陸上輸油管道溢油擴(kuò)散的三維態(tài)勢(shì)推演模型

摘 要 為了評(píng)估陸上輸油管道溢油的危害范圍，以陸上輸油管道溢油擴(kuò)散事故為研究對(duì)象，設(shè)計(jì)三維態(tài)勢(shì)推演模型。將計(jì)算得到的陸上輸油管道溢油速度、溢油量等參數(shù)作為模型的輸入項(xiàng)，設(shè)置質(zhì)量守恒和動(dòng)量守恒作為三維態(tài)勢(shì)推演模型的約束條件。模擬陸上輸油管道的溢油過(guò)程，在相關(guān)條件的約束下，實(shí)現(xiàn)橫向擴(kuò)散和縱向滲透的態(tài)勢(shì)推演，綜合考慮蒸發(fā)和吸附現(xiàn)象，調(diào)整態(tài)勢(shì)推演結(jié)果，輸出橫向擴(kuò)散范圍、縱向滲透深度等參數(shù)的推演結(jié)果。通過(guò)模型分析實(shí)驗(yàn)得出結(jié)論：與傳統(tǒng)推演模型相比，優(yōu)化設(shè)計(jì)模型的橫向擴(kuò)散面積推演誤差和縱向擴(kuò)散深度推演誤差分別降低約3.57 m2和0.2 m，即優(yōu)化設(shè)計(jì)模型在態(tài)勢(shì)推演精度性能方面具有明顯優(yōu)勢(shì)。

關(guān)鍵詞 陸上輸油管道 溢油擴(kuò)散 三維態(tài)勢(shì) 推演模型

中圖分類號(hào) X55 " 文獻(xiàn)標(biāo)志碼 A " 文章編號(hào) 0254?6094（2024）04?0613?07

輸油管道由油管及其附件構(gòu)成，并根據(jù)工藝流程的要求，配置與之相適應(yīng)的油泵單元，將其設(shè)計(jì)并安裝為一個(gè)完整的管路系統(tǒng)，以實(shí)現(xiàn)油料的接卸和輸轉(zhuǎn)。輸油管道的管子通常是鋼管，通過(guò)焊接、法蘭等連接方式將長(zhǎng)距離管道連接起來(lái)，并通過(guò)閥門(mén)進(jìn)行開(kāi)閉控制和流量調(diào)節(jié)。由于輸油管道的長(zhǎng)期使用和接口位置的腐蝕，導(dǎo)致輸油管道發(fā)生溢油現(xiàn)象［1］。溢油是指在石油勘探、開(kāi)發(fā)、煉制及儲(chǔ)運(yùn)過(guò)程中，因意外事故或操作失誤，導(dǎo)致原油或油品從作業(yè)現(xiàn)場(chǎng)或儲(chǔ)器里外泄。根據(jù)輸油管道的溢油位置，可以分為陸上溢油和水下溢油兩種類型，陸上溢油指的是輸油管道中的傳輸石油通過(guò)缺口泄漏至陸地表面。溢油造成了巨大的石油資源的流失，并且由于殘余原油的黏性、稠度顯著增加，使得溢油的收集與處理變得十分困難。石油泄漏不僅會(huì)對(duì)大氣、生態(tài)及岸灘等造成嚴(yán)重影響，而且還會(huì)對(duì)人類身體造成嚴(yán)重威脅，限制了經(jīng)濟(jì)與社會(huì)的可持續(xù)發(fā)展。

輸油管道溢油會(huì)在陸上向水平、豎直等多個(gè)維度擴(kuò)散，為了實(shí)現(xiàn)對(duì)輸油管道溢油的處理，需要對(duì)溢油的擴(kuò)散形式進(jìn)行推演，最大程度上降低陸上輸油管道溢油產(chǎn)生的負(fù)面影響。為此，筆者設(shè)計(jì)了陸上輸油管道溢油擴(kuò)散的三維態(tài)勢(shì)推演模型。現(xiàn)階段發(fā)展較為成熟的三維態(tài)勢(shì)推演模型包括：基于風(fēng)洞平臺(tái)實(shí)驗(yàn)的三維態(tài)勢(shì)推演模型、基于貝葉斯網(wǎng)絡(luò)的三維態(tài)勢(shì)推演模型和基于動(dòng)態(tài)摩阻模型的三維態(tài)勢(shì)推演模型，上述傳統(tǒng)模型在實(shí)際工作過(guò)程中存在推演誤差大的問(wèn)題，降低了推演模型的應(yīng)用性能，為此以陸上輸油管道的溢油擴(kuò)散行為作為研究對(duì)象，優(yōu)化設(shè)計(jì)三維態(tài)勢(shì)推演模型，以期能夠提升三維態(tài)勢(shì)推演模型的推演與應(yīng)用性能。

1 三維態(tài)勢(shì)推演模型設(shè)計(jì)

優(yōu)化設(shè)計(jì)陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演模型的基本運(yùn)行思路為：根據(jù)陸上輸油管道的破口尺寸，確定輸油管道的溢油量、溢油速度等參數(shù)，將其作為三維態(tài)勢(shì)推演模型的輸入?yún)?shù)，考慮陸上輸油管道中輸送石油的物理性質(zhì)，設(shè)置三維態(tài)勢(shì)推演模型的約束條件，具體包括動(dòng)量守恒、質(zhì)量守恒等，模擬陸上輸油管道的溢油過(guò)程，最終從橫向擴(kuò)展、縱向向上蒸發(fā)、縱向向下滲透等多個(gè)維度，推演溢油的擴(kuò)散態(tài)勢(shì)，并得出陸上輸油管道溢油區(qū)域內(nèi)任意位置在任意時(shí)刻的溢油量及存在形式，以此作為三維態(tài)勢(shì)推演模型的量化輸出。

1.1 計(jì)算陸上輸油管道溢油參數(shù)

陸上輸油管道溢油參數(shù)是三維態(tài)勢(shì)推演模型的輸入項(xiàng)，具體的溢油參數(shù)包括溢油速度、溢油量、溢油溫度及空間位置等，陸上輸油管道與外界環(huán)境之間的壓強(qiáng)差直接決定輸油管道的溢油速度，管道內(nèi)外的壓強(qiáng)差Δp可以表示為：

Δp=p-p （1）

p=ρgD （2）

式中 D——陸上輸油管道的直徑；

"g——重力加速度；

" "p——陸上輸油管道外的壓強(qiáng)值；

p——陸上輸油管道內(nèi)壓強(qiáng)值；

"ρ——陸上輸油管道內(nèi)輸送的石油量。

將式（2）的計(jì)算結(jié)果代入式（1）中，即可得出壓強(qiáng)差的計(jì)算結(jié)果［2］。根據(jù)輸油管道壓強(qiáng)差的計(jì)算結(jié)果，確定陸上輸油管道的溢油速度v為：

v=Δp××W（3）

式中 S——陸上輸油管道上的破口面積；

S——陸上輸油管道的總面積；

" W——單位體積內(nèi)輸送的石油總量。

當(dāng)溢油量不變時(shí)，由于在溢油點(diǎn)附近有少量的油滴，在溢油量不變的情況下，溢油率增加，油滴擴(kuò)散范圍擴(kuò)大。在溢油模型中，表征溢油特征的溢油量和溢油率是影響溢油擴(kuò)散的重要因素［3］。根據(jù)溢油速度的計(jì)算結(jié)果，確定陸上輸油管道的溢油量，溢油量W的計(jì)算公式為：

W=v×S×W " "（4）

將相關(guān)數(shù)據(jù)代入式（4）中，即可得出溢油量參數(shù)的計(jì)算結(jié)果。另外溢油溫度參數(shù)的計(jì)算結(jié)果為：

T=κT " " "（5）

式中 T——輸油管道內(nèi)石油的傳輸溫度［4］；

κ——溫度衰減系數(shù)，該參數(shù)的具體取

值與環(huán)境溫度有關(guān)。

同理，可以在硬件設(shè)備的支持下，得出陸上輸油管道其他溢油參數(shù)的計(jì)算結(jié)果。除了陸上輸油管道的溢油參數(shù)外，溢油位置區(qū)域的環(huán)境參數(shù)也會(huì)對(duì)溢油擴(kuò)散效果產(chǎn)生影響，因此還需要收集區(qū)域環(huán)境內(nèi)地面粗糙度、地面坡度等參數(shù)，并將上述所有參數(shù)以數(shù)據(jù)的形式輸入到三維態(tài)勢(shì)推演模型中。

1.2 設(shè)置三維態(tài)勢(shì)推演模型約束條件

陸上輸油管道溢油擴(kuò)散的三維推演需遵循質(zhì)量守恒和動(dòng)量守恒規(guī)則，溢油控制單元體在陸上擴(kuò)散過(guò)程中，周圍環(huán)境對(duì)溢油的溶解、滲透及分解等作用引起的液體質(zhì)量變化滿足質(zhì)量守恒定律，質(zhì)量守恒定律約束條件的設(shè)置如下：

=-∑- " （6）

其中，和分別為滲透到陸上地表下的溢油質(zhì)量損失和水平擴(kuò)散引起的溢油質(zhì)量損失［5］。溢油在陸地環(huán)境輸移和擴(kuò)散過(guò)程中，周圍環(huán)境元素不斷摻混進(jìn)入溢油內(nèi)，并滿足如下動(dòng)量守恒方程：

其中，κ為拖拽系數(shù)，將模型輸入的陸上輸油管道溢油初始參數(shù)代入式（7）中，得出動(dòng)量守恒方程的構(gòu)建結(jié)果［6］。在溢油擴(kuò)散三維態(tài)勢(shì)推演過(guò)程中，要求推演結(jié)果均滿足上述守恒條件。

1.3 模擬陸上輸油管道溢油過(guò)程

在陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演模型前，首先需要確定溢油在陸上環(huán)境中的基本擴(kuò)散規(guī)律，圖1為陸上輸油管道的溢油過(guò)程。

輸油管道的滲漏過(guò)程可劃分為3個(gè)階段。在第1個(gè)階段，石油在初動(dòng)量的驅(qū)動(dòng)下，以固定的速度被噴入陸面，并形成了一道油柱，這時(shí)石油的動(dòng)量很大，初動(dòng)量占主導(dǎo)地位。第2個(gè)階段是在溢油運(yùn)動(dòng)到一定距離后，在這一階段，受外部環(huán)境的影響，溢油的初始噴流動(dòng)能不斷減小，油柱逐漸破裂為無(wú)數(shù)粒徑大小不一的微小油珠，這些微小油珠與海水混合并一起向外擴(kuò)散［7］。油污在被卷吸的同時(shí)，會(huì)與周圍的其他物質(zhì)發(fā)生混雜，改變油污的組成和性質(zhì)。當(dāng)溢油持續(xù)擴(kuò)散至某一區(qū)域后，油滴將徹底喪失原有的動(dòng)力，從而使溢油的運(yùn)動(dòng)狀態(tài)進(jìn)入第3個(gè)階段。在此期間，油滴的初始動(dòng)量已全部消失，且與周圍的物質(zhì)發(fā)生了同化，因此，大多數(shù)油滴在橫向和縱向上都是依賴于周圍的動(dòng)量來(lái)進(jìn)行擴(kuò)散的［8］。在陸上輸油管道溢油擴(kuò)散過(guò)程中，存在外部阻力，導(dǎo)致擴(kuò)散現(xiàn)象的停止，若存在水流等其他影響因素，則擴(kuò)散現(xiàn)象、走向?qū)?huì)發(fā)生改變。當(dāng)陸上輸油管道溢油在水平和垂直方向均不再發(fā)生位移時(shí)，溢油擴(kuò)散現(xiàn)象結(jié)束。

1.4 實(shí)現(xiàn)陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演

綜合考慮陸上輸油管道溢油過(guò)程的3個(gè)擴(kuò)散階段，在三維環(huán)境下進(jìn)行態(tài)勢(shì)推演，圖2表示的是輸油管道溢油擴(kuò)散的三維環(huán)境。

陸上輸油管道溢油在水平維度中的擴(kuò)散包括慣性擴(kuò)散、黏性擴(kuò)散、表面張力擴(kuò)散和擴(kuò)散停止4個(gè)階段，在模型約束條件的控制作用下，綜合考慮影響擴(kuò)散現(xiàn)象的多個(gè)因素，得出三維態(tài)勢(shì)的推演結(jié)果。

1.4.1 溢油橫向擴(kuò)散態(tài)勢(shì)推演

陸上輸油管道溢油地表橫向擴(kuò)展主要用于描述溢油在不同地表狀態(tài)下的流動(dòng)擴(kuò)散情況，確定溢油的擴(kuò)散范圍和到達(dá)指定區(qū)域的時(shí)間。從水文地質(zhì)的角度，分析了不同的地形條件、不同的地表粗糙度和不同的表土微孔隙結(jié)構(gòu)對(duì)石油的擴(kuò)散有不同的影響。不同地形對(duì)溢油流態(tài)的影響也有很大差別［9］。與地面的水流過(guò)程相似，地面上的油膜也會(huì)受到地面地形的影響。油膜在坡面更容易產(chǎn)生溢油，但過(guò)高的地表粗糙度和過(guò)多的微細(xì)孔隙將使溢油擴(kuò)散區(qū)域變小。擴(kuò)展過(guò)程受到重力、表面張力和黏滯力的共同作用，在不同的受力情況下，陸上輸油管道溢油的擴(kuò)散面積S可以表示為：

S=

，Y

，Y

，Y

（8）

式中 k——經(jīng)驗(yàn)系數(shù)；

" k——運(yùn)動(dòng)黏滯系數(shù)；

" t——溢油在陸上的擴(kuò)散時(shí)間；

Y、Y、Y——溢油在陸上地表的3個(gè)擴(kuò)散階

段［10］；

γ——溢油擴(kuò)散區(qū)域內(nèi)的地面粗糙度；

"δ——溢油的凈表面張力；

ρ——環(huán)境密度；

ρ——溢油密度。

隨著時(shí)間的推移，利用式（8）可以得出不同時(shí)段溢油在橫向維度上擴(kuò)散直徑的推演結(jié)果。

1.4.2 溢油縱向滲透態(tài)勢(shì)推演

陸上輸油管道溢油在豎直方向上的態(tài)勢(shì)推演目標(biāo)就是確定溢油向下的滲透深度，溢油在土體中擴(kuò)散的過(guò)程中，因其周圍土體的孔隙中常含有水和空氣，油品在土壤顆粒中的流動(dòng)過(guò)程，通常伴隨水和空氣的同時(shí)流動(dòng)。在陸上輸油管道溢油擴(kuò)散過(guò)程中，石油和水是互不相溶的，因此溢油在縱向維度中的擴(kuò)散過(guò)程就是一個(gè)多孔介質(zhì)多相流過(guò)程［11］。陸上輸油管道溢油在任意時(shí)刻的縱向滲透深度h的推演結(jié)果可以表示為：

h=W×λ×ρ×t " （9）

式中 λ——地面坡度；

ρ——土壤密度。

將相關(guān)參數(shù)代入式（9）中，得出溢油縱向滲透深度的推演結(jié)果，將其作為縱向維度態(tài)勢(shì)的推演輸出結(jié)果。

1.4.3 溢油蒸發(fā)態(tài)勢(shì)推演

在陸上輸油管道溢油擴(kuò)散過(guò)程中，受到溢油溫度和環(huán)境溫度的影響，導(dǎo)致溢油可能出現(xiàn)蒸發(fā)情況，溢油的蒸發(fā)直接導(dǎo)致溢油量的變化，進(jìn)而影響陸上輸油管道溢油在三維環(huán)境中的擴(kuò)散效果。輸油管道溢油的蒸發(fā)率μ為：

μ=κ×p× （10）

式中 M——溢油的分子質(zhì)量；

p——溢油的蒸氣壓；

T——溢油溫度；

α——?dú)怏w常數(shù)；

κ——蒸發(fā)系數(shù)。

從式（10）中可以看出，溢油溫度越高，溢油的蒸發(fā)率越高，即轉(zhuǎn)換為氣體的溢油數(shù)量越多，此時(shí)溢油擴(kuò)散至空氣中，同時(shí)地表橫向與縱向擴(kuò)散維度中的初始溢油數(shù)量發(fā)生變化，剩余的溢油量[W][′]可以表示為：

[W][′]=（1-μ）W "（11）

利用式（11）對(duì)上述橫向和縱向擴(kuò)展態(tài)勢(shì)推演中的溢油量參數(shù)進(jìn)行調(diào)整，得到考慮蒸發(fā)條件下的態(tài)勢(shì)推演調(diào)整結(jié)果。

1.4.4 溢油吸附態(tài)勢(shì)推演

土壤對(duì)石油類污染物有一定的吸附作用，土壤對(duì)溢油的吸附能力Z可以表示為：

Z=L×K×H× χ （12）

式中 H——能夠?qū)崿F(xiàn)溢油吸附土壤區(qū)域的高度；

K——能夠?qū)崿F(xiàn)溢油吸附土壤區(qū)域的寬度；

L——能夠?qū)崿F(xiàn)溢油吸附土壤區(qū)域的長(zhǎng)度；

χ——吸附土壤的有效孔隙率［12］。

根據(jù)土壤的吸附能力，可以確定輸油管道溢油在土壤區(qū)域內(nèi)被吸附的概率，從而對(duì)向下滲透的溢油量進(jìn)行調(diào)整，得出考慮土壤吸附作用下的溢油擴(kuò)散態(tài)勢(shì)推演結(jié)果。

綜合考慮蒸發(fā)和吸附溢油，得出任意時(shí)刻溢油的擴(kuò)散推演結(jié)果，并將推演結(jié)果以量化形式輸出，完成模型的態(tài)勢(shì)推演功能。

2 模型分析

為了測(cè)試優(yōu)化設(shè)計(jì)陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演模型的態(tài)勢(shì)推演性能，設(shè)計(jì)模型分析實(shí)驗(yàn)，此次實(shí)驗(yàn)的基本執(zhí)行思路為：在某實(shí)驗(yàn)環(huán)境中在人為控制作用下產(chǎn)生輸油管道溢油現(xiàn)象，并通過(guò)對(duì)輸電管線工作參數(shù)與破口大小的控制，得出輸油管道溢油參數(shù)的設(shè)置結(jié)果。分別設(shè)置多個(gè)不同的態(tài)勢(shì)推演目標(biāo)時(shí)刻，利用優(yōu)化設(shè)計(jì)模型得出對(duì)應(yīng)時(shí)刻的態(tài)勢(shì)推演結(jié)果，并在對(duì)應(yīng)時(shí)刻利用可視化觀測(cè)設(shè)備和測(cè)量設(shè)備，得出水平擴(kuò)散范圍和縱向擴(kuò)散深度的實(shí)際值，通過(guò)模型輸出推演結(jié)果數(shù)據(jù)與實(shí)際測(cè)定數(shù)據(jù)之間的差距，得出反映優(yōu)化設(shè)計(jì)模型推演精度的測(cè)試結(jié)果。

2.1 選擇陸上輸油管道研究對(duì)象

此次模型分析實(shí)驗(yàn)選擇某輸油管道的陸上部分作為研究對(duì)象，該輸油管道直徑為1.2 m，能夠同時(shí)實(shí)現(xiàn)5 t石油資源的運(yùn)輸，最快輸送速度為20 t/h，由于石油具有較高的易燃性，為保證石油資源的輸送安全，選擇輸油管道一般將輸送速度控制在10 t/h。整條輸油管道有2/3埋設(shè)在地面以下，僅1/3在陸地表面，選擇陸上輸油管道研究對(duì)象的長(zhǎng)度約為6 km。

2.2 設(shè)定陸上輸油管道溢油現(xiàn)象

在選擇的陸上輸油管道區(qū)域選擇任意位置，采用人為破壞的方式生成一個(gè)10 mm×10 mm的破口，實(shí)現(xiàn)陸上輸油管道的溢油現(xiàn)象。在溢油現(xiàn)象生成狀態(tài)下，輸油管道中石油資源的輸送速度為10 t/h，輸送總量為20 t，并平均分布整個(gè)管道區(qū)域內(nèi)，確定溢油總量為0.32 t，初始溢油速度與輸送速度相同。輸油管道溢油區(qū)域土壤的孔隙率約為0.2，地面粗糙度和坡度分別為0.1和0.6°，以此作為輸入模型中的初始數(shù)據(jù)。

2.3 描述模型測(cè)試實(shí)驗(yàn)過(guò)程

利用程序轉(zhuǎn)換工具，實(shí)現(xiàn)對(duì)優(yōu)化設(shè)計(jì)陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演模型的開(kāi)發(fā)，并保證推演模型能夠在主測(cè)計(jì)算機(jī)中正常運(yùn)行。輸入準(zhǔn)備的初始溢油數(shù)據(jù)，以陸上輸油管道溢油現(xiàn)象的發(fā)生時(shí)刻作為初始時(shí)刻，分別設(shè)置5、10、15、20、25、30、45、60 min為模型的態(tài)勢(shì)推演目標(biāo)時(shí)刻，通過(guò)模型的運(yùn)行與推演，輸出水平（橫向）和豎直（縱向）兩個(gè)方向上的態(tài)勢(shì)推演結(jié)果。圖3表示的是5 min時(shí)刻模型輸出的態(tài)勢(shì)推演結(jié)果。

按照上述方式可以得出實(shí)驗(yàn)中設(shè)置所有推演目標(biāo)時(shí)刻模型輸出的三維態(tài)勢(shì)推演結(jié)果。分別在目標(biāo)時(shí)刻利用硬件設(shè)備測(cè)量輸油管道溢油的實(shí)際擴(kuò)散半徑與深度，得出溢油擴(kuò)散的真實(shí)數(shù)據(jù)。其中5 min時(shí)刻溢油水平擴(kuò)散情況的實(shí)際測(cè)量結(jié)果如圖4所示。

按照上述方式得出所有時(shí)刻陸上輸油管道溢油在水平和豎直方向上擴(kuò)散的實(shí)際數(shù)據(jù)。為體現(xiàn)出優(yōu)化設(shè)計(jì)模型在推演精度方面的優(yōu)勢(shì)，設(shè)置傳統(tǒng)的基于風(fēng)洞平臺(tái)實(shí)驗(yàn)的三維態(tài)勢(shì)推演模型和基于貝葉斯網(wǎng)絡(luò)的三維態(tài)勢(shì)推演模型作為實(shí)驗(yàn)的對(duì)比模型，在相同的實(shí)驗(yàn)環(huán)境下實(shí)現(xiàn)對(duì)比模型的開(kāi)發(fā)與運(yùn)行，并得出相應(yīng)的模型輸出結(jié)果數(shù)據(jù)。

2.4 設(shè)置模型測(cè)試分析指標(biāo)

為了實(shí)現(xiàn)對(duì)模型推演性能的量化測(cè)試，設(shè)置橫向擴(kuò)散面積推演誤差和縱向擴(kuò)散深度推演誤差作為模型性能的量化測(cè)試指標(biāo)，其中橫向擴(kuò)散面積推演誤差ε的數(shù)值結(jié)果為：

ε=S

-S

（13）

式中 S——橫向擴(kuò)散面積的實(shí)測(cè)結(jié)果；

" "S——橫向擴(kuò)散面積的推演結(jié)果。

縱向擴(kuò)散深度推演誤差ε的測(cè)試結(jié)果如下：

ε=h

-h

（14）

式中 h——縱向擴(kuò)散深度的模型的

實(shí)測(cè)結(jié)果；

"h——縱向擴(kuò)散深度的模型輸

出結(jié)果。

最終計(jì)算得出推演誤差越小，證明對(duì)應(yīng)模型的推演精度性能越優(yōu)。

2.5 模型分析實(shí)驗(yàn)結(jié)果

通過(guò)相關(guān)數(shù)據(jù)的統(tǒng)計(jì)，得出反映陸上輸油管道溢油擴(kuò)散三維態(tài)勢(shì)推演模型精度性能的測(cè)試結(jié)果，見(jiàn)表1。

將表1中的數(shù)據(jù)代入式（13）、（14）中，計(jì)算得出兩種傳統(tǒng)模型橫向擴(kuò)散面積推演誤差的平均值分別為5.75、2.88 m2，縱向擴(kuò)散深度推演誤差的平均值分別為0.34、0.16 m，另外優(yōu)化設(shè)計(jì)模型橫向擴(kuò)散面積和縱向擴(kuò)散深度的推演誤差平均值分別為0.75 m2和0.05 m。

3 結(jié)束語(yǔ)

由于受到腐蝕、自然災(zāi)害或第三方損害等因素的影響，導(dǎo)致輸油管道出現(xiàn)安全事故的現(xiàn)象并不少見(jiàn)。管道輸送原油的易燃、易爆和劇毒特性決定了管道輸送原油的燃爆危險(xiǎn)性高，波及范圍廣，對(duì)環(huán)境造成嚴(yán)重的污染。在此次研究中，設(shè)計(jì)并開(kāi)發(fā)了陸上輸油管道溢油擴(kuò)散的三維態(tài)勢(shì)推演模型，能夠?qū)斢凸艿酪缬偷臄U(kuò)散趨勢(shì)進(jìn)行推演與預(yù)測(cè)，為管道泄漏擴(kuò)散應(yīng)急方案的制定與執(zhí)行提供有效參考，降低陸上輸油管道溢油對(duì)環(huán)境的影響。

參 考 文 獻(xiàn)

［1］ 許雪，陳風(fēng)，黃維秋，等.基于風(fēng)洞平臺(tái)實(shí)驗(yàn)的大型罐區(qū)溢油事故后的油氣擴(kuò)散模擬［J］.環(huán)境工程學(xué)報(bào)，2021，15（12）：3946-3956.

［2］ 屈靜，張建彬，李旭芳，等.基于貝葉斯網(wǎng)絡(luò)的輸油管道泄漏事故情景推演［J］.中國(guó)安全科學(xué)學(xué)報(bào)，2021，31（1）：192-198.

［3］ 顧一凡，王彬諭，倪云林，等.綠華島周邊海域溢油擴(kuò)散風(fēng)險(xiǎn)數(shù)值模擬研究［J］.海洋湖沼通報(bào)，2022，44（2）：17-26.

［4］ 孫琰，劉大召，許源興，等.湛江灣海域溢油漂移擴(kuò)散數(shù)值模擬研究［J］.海洋湖沼通報(bào)，2022，44（2）：73-83.

［5］ 顧靖華，朱建榮，金智.長(zhǎng)江口寶鋼碼頭溢油事故油膜漂移擴(kuò)散數(shù)值模擬［J］.熱帶海洋學(xué)報(bào)，2022，41（6）：159-170.

［6］ 李旺，馮亮，黃立宇，等.基于動(dòng)態(tài)摩阻模型的大落差輸油管道不穩(wěn)定流動(dòng)模擬［J］.科學(xué)技術(shù)與工程，2022，22（9）：3512-3517.

［7］ 李巖松，丁鼎倩，韓東，等.起伏輸油管道臨界完全攜積水油速數(shù)值模擬［J］.上海交通大學(xué)學(xué)報(bào)，2021，55（7）：878-890.

［8］ 劉雪光，張偉浩，吳東容，等.中緬原油管道瀾滄江跨越段原油泄漏擴(kuò)散數(shù)值模擬［J］.油氣儲(chǔ)運(yùn)，2021，40（1）：96-106.

［9］ 王鵬，申霞，王船海，等.復(fù)雜河網(wǎng)水系油粒子模型開(kāi)發(fā)及溢油污染模擬［J］.水利學(xué)報(bào)，2021，52（6）：750-757.

［10］ 陳勤思，胡松，張春玲，等.海洋溢油溯源數(shù)值模擬方法的研究進(jìn)展［J］.海洋通報(bào)，2022，41（5）：481-490.

［11］ 陳秋良，蔡學(xué)石，姜倩，等.曹妃甸甸頭海域溢油行為數(shù)值模擬研究［J］.海洋環(huán)境科學(xué)，2022，41（1）：85-93.

［12］ 熊海濱，孫昭華，陳立，等.碼頭密度對(duì)長(zhǎng)江中游典型河段內(nèi)溢油運(yùn)移影響數(shù)值模擬［J］.安全與環(huán)境學(xué)報(bào)，2021，21（5）：2252-2261.

（收稿日期：2023-09-05，修回日期：2024-07-12）

A 3D Situation Prediction Model for Oil Spills and Diffusion from

Onshore Oil Pipelines

SHI Ying?he

（China Petroleum Longway Engineering Project Management Co.，Ltd.）

Abstract " For purpose of evaluating the hazard range of oil spills from the onshore oil pipeline， the oil spill diffusion accident of onshore oil pipeline was taken as the object of study， and a 3D situation simulation model was designed， in which， having parameters of oil spills’velocity and amount there calculated taken as the input items of the model， and the mass conservation and momentum set as the constraints of this model. In simulating the oil?spilling process， the situation of lateral diffusion and vertical penetration was realized under the constraints of relevant conditions， including having the evaporation and adsorption phenomena considered comprehensively， and the situation results adjusted to output the simulation results of both lateral diffusion range and vertical penetration depth. The model analysis experiment indicates that， compared with the traditional simulation model， this optimal model can reduce the simulation error of lateral diffusion area and the longitudinal diffusion depth by about 3.57 m2 and 0.2 m respectively， which means that the optimized design model has obvious advantages in accurately simulating the oil spills and diffusion.

Key words " "onshore oil pipeline， oil spill and diffusion， 3D situation， simulation model