埋地天然氣管道泄漏擴散數(shù)值模擬

王 健

（新疆油田公司油氣儲運公司，新疆 昌吉 831100）

埋地天然氣管道泄漏擴散數(shù)值模擬

王 健

（新疆油田公司油氣儲運公司，新疆 昌吉 831100）

針對埋地天然氣管道穿孔泄漏擴散問題，結合有限容積法，利用Gambit2.4建立了天然氣管道不同泄漏口徑和不同環(huán)境溫度的CFD仿真模型，利用Fluent6.3分別對不同泄漏口直徑（6.35、25.4、101.6 mm）和不同環(huán)境溫度（0、10、30 ℃）泄漏工況下，氣體在土壤中和空氣中擴散規(guī)律進行了數(shù)值模擬。研究結果表明，隨著泄漏口直徑增加天然氣危害范圍逐漸增大，關閉泄漏管段兩端閥門以后，氣體擴散危害范圍逐漸變?。浑S著環(huán)境溫度的升高氣體危害范圍在豎直方向上有明顯增加，在地面處危險區(qū)域也增加。研究結果為城市埋地天然氣管道泄漏事故現(xiàn)場人員疏散及安全搶修提供了理論依據(jù)。

天然氣管道；泄漏；擴散；數(shù)值模擬

隨著我國國民經(jīng)濟的不斷發(fā)展和人民生活水平的提高，天然氣的應用得到迅速的發(fā)展，燃氣作為重要的能源形式之一應用越來越廣泛。城市燃氣管道一般采用埋地方式鋪設，但由于受腐蝕和人為破壞等因素的作用，天然氣管道泄漏事故時有發(fā)生[1,2]。天然氣爆炸極限范圍較大，具有易燃、易爆的性質，天然氣泄漏爆炸事故一般會造成大量人員傷亡和嚴重的財產損失，具有較大的危害性。因此，對天然氣管道泄漏后，確定其危害范圍具有重大意義。

本文針對城市埋地天然氣管道，利用C語言編寫與FLUENT軟件兼容的UDF計算程序，修正風速隨高度的變化，并且控制閥門由于人工或自動關閉后氣體泄漏率[3]隨時間的變化。由于腐蝕、施工等不同的損傷過程，造成泄漏口的尺寸不同，根據(jù)API581-2000—基于風險的檢驗基礎資源文件（Risk-Based Inspection Base Resource Document）選取管道泄漏口直徑分別為6.35 mm、25.4 mm、 101.6mm，并且根據(jù)不同的環(huán)境溫度，建立泄漏模型來模擬計算，得出了天然氣管道泄漏后天然氣擴散規(guī)律，為制定天然氣管道泄漏事故應急預案及保障措施提供了理論依據(jù)。

1 模型建立

1.1 數(shù)學及物理模型

以天然氣管道泄漏為例，模擬泄漏過程及天然氣流場變化。天然氣為多種氣體的混合物，主要成分為甲烷(CH4)（體積分數(shù)為99.88%），故可以用CH4代替天然氣進行模擬計算。天然氣的傷害下限為1%，爆炸極限范圍為5%～15%[4]，故以天然氣的體積分數(shù)為1%、5%和15%為分界線為評定泄漏氣體危害濃度區(qū)域。

天然氣泄漏過程遵循連續(xù)方程、動量守恒方程以及能量守恒方程[5]。綜合考慮計算條件和工程要求，采用耦合隱式求解器，其基本方程為定常Reynolds時均方程[6]。同時天然氣泄漏為湍流流動，故選用湍流模型中的可實現(xiàn)k-ε模型。

某天然氣管道運行壓力為0.2 MPa，根據(jù)GB 50028-2006《城市燃氣設計規(guī)范》，屬于中壓輸氣管道，管道直徑325 mm，壁厚8 mm，管道埋深1.7 m，氣體溫度為20 ℃，由于所選天然氣管道在城鎮(zhèn)人口密集地區(qū)，考慮管道所處的實際位置，取建筑物高度為20 m，管道中心距建筑物距離為3 m[7]，土壤物性為：密度1 500 kg/m3，導熱系數(shù)1.5 W/ (m?K)，比熱2 200 J/(kg?K )，孔隙度0.43[8,9]，土壤平均直徑0.2 mm，根據(jù)文獻[10]計算出滲透率及慣性阻力系數(shù)。用較大空間 (80 m×80 m)進行仿真模擬計算，具體模型及邊界條件如圖1所示。劃分網(wǎng)格時，由于管道泄漏口附近氣體各項參數(shù)變化梯度較大，為滿足計算精度，對此區(qū)域進行網(wǎng)格加密。

圖1 物理模型及局部放大圖Fig.1 Physical model and partial enlargement

1.2 參數(shù)確定及邊界條件設置

假設管道兩個閥門之間的距離為1 000 m[7]。泄漏開始為穩(wěn)態(tài)泄漏，泄漏壓力為0.2 MPa，此壓力大于臨界壓力Pc，泄漏氣體流動為臨界流動[11]，泄漏氣體質量流率Qm0按式（7）計算得Qm0=0.181 964 kg/s，G=1 077.95 kg/（m2?s）。泄漏120 s時閥門自動關閉[11]，關閉閥門后管道內壓力減小，當壓力大于臨界壓力Pc時，泄漏氣體質量流率Qm0按式（7）、（9）、（10）、（11）計算；當壓力小于臨界壓力Pc時，泄漏氣體流動為亞臨界流，泄漏氣體質量流率Qm0按式（8）、（9）、（10）、（11）計算。

式中：

Qm— 質量流量，kg/s；

k — 比熱容；

t — 關閉閥門后泄漏時間，s；

Aor— 泄漏口面積，m2；

R — 氣體常數(shù)，R=8.314 J/(mol?K)；

P2— 管道內壓力，MPa；

Pa— 大氣壓，0.101 325 MPa；

G — 質量流率，kg/（m2?s）；

T — 溫度，K；

L — 管線長度，m；

ρ— 密度，kg/m3。

穩(wěn)態(tài)泄漏120 s內質量流量進口設置為1 077.95 kg/（m2?s），關閉閥門后，根據(jù)質量流量流率隨時間變化擬合出公式（14）：

將公式（14）編寫成UDF程序并按照風速邊界設置方法設置泄漏口邊界條件。

2 數(shù)值模擬結果和分析

埋地天然氣管道發(fā)生泄漏擴散后，由于受土壤阻力作用，泄漏氣體的動能衰減程度較大，速度大大降低，泄漏氣體不足以沖破管道上方的埋土，氣體會在土壤中擴散，到達地表后向大氣中擴散[12]。

根據(jù)不同泄漏口直徑修改幾何模型，以泄漏口在天然氣管道上表面為例說明，氣體擴散云圖如圖2-圖4所示。

2.1 不同泄漏孔大小下中壓天然氣管道泄漏

中壓天然氣管道泄漏不同泄漏孔徑質量流量如表1。

表1 中壓天然氣管道泄漏不同泄漏孔徑質量流量Table 1 Medium pressure natural gas pipeline leakage mass flow of different aperture

圖2-圖4分別為泄漏口為6.35、25.4和101.6 mm下天然氣泄漏擴散云圖(a—泄漏120 s，b—泄漏停止時)，可以看出當泄漏120 s時，泄漏孔徑為 6.35 mm，由于泄漏孔徑很小，天然氣泄漏速度和質量流量也很小，天然氣和大氣對流擴散比較慢，天然氣危險范圍較??；當泄漏孔徑為 25.4 mm 時，甲烷危險范圍增大，傷害下限范圍和爆炸下限范圍在下風向貼近地面處分別為36 m和22 m，擴散高度分別為大于75 m和60 m，在土壤中沿下風向擴散距離分別為38 m和20 m；當泄漏孔徑為101.6 mm時，由于泄漏率的增加，氣體危險濃度范圍進一步擴大，傷害下限范圍已經(jīng)超過了計算區(qū)域即下風向擴散距離超過了38 m，而高度遠遠超過了計算區(qū)域的上邊界，爆炸下限在下風向貼近地面處擴散距離為22 m，高度超過了計算區(qū)域上邊界，在土壤中沿下風向擴散距離為25 m；可見泄漏孔徑變大，泄漏質量流率變大，天然氣危害范圍也較大。

從圖2-圖4可以看出，隨著泄漏孔徑的變大，距離泄漏口相同距離處，天然氣濃度會加大，而泄漏擴散危險區(qū)域也會變大。

當泄漏口直徑為6.35 mm時，泄漏口直徑很小可以看做穩(wěn)態(tài)泄漏，即氣體以恒定泄漏率泄漏，由于風的擴散作用，天然氣危害范圍變化不大。而對于泄漏口直徑分別為25.4 mm和101.6 mm時，當關閉泄漏口管道兩端閥門后，管道內壓力逐漸降低，從而氣體泄漏率逐漸減小，氣體擴散危險范圍逐漸減小，當管道內壓與大氣壓力相等時，泄漏率接近于零，如圖3(b)、4(b)所示。

當泄漏口徑為25.4 mm時，泄漏口處氣體質量流量逐漸減小，在434 s左右泄漏口處質量流量接近于零，故氣體危險濃度范圍也在變小，傷害下限范圍和爆炸下限范圍在下風向貼近地面處分別為28和13 m，擴散高度分別為大于68和52 m，在土壤中沿下風向擴散距離分別為29.8和16.5 m；而泄漏口為101.6 mm時，由于泄漏口較大，氣體質量流量很大，在19 s左右管內氣體質量流量已接近零，由于從關閉閥門到泄漏停止之間時間較短，氣體危險濃度范圍變化不大。

圖2 泄漏口直徑d=6.35 mm CH4濃度分布圖Fig.2 CH4Concentration distribution when leakage aperture is 6.35 mm

2.2 不同環(huán)境溫度下中壓天然氣管道泄漏

天然氣溫度恒定為20 ℃，泄漏口直徑取25.4 mm，泄漏口在管道上表面處，環(huán)境溫度分別選取0、10及30 ℃，在平坦地區(qū)對天然氣泄漏擴散過程進行數(shù)值模擬，模擬結果如圖5-圖7可以得到表2。

圖3 泄漏口直徑d=25.4 mm CH4濃度分布圖Fig.3 CH4Concentration distribution when leakage aperture is 25.4 mm

圖4 泄漏口直徑d=101.6 mm CH4濃度分布圖Fig.4 CH4Concentration distribution when leakage aperture is 101.6 mm

表2 中壓天然氣管道泄漏不同環(huán)境溫度下氣體危害范圍Table 2 Hazard range of medium pressure gas pipeline leakage under different ambient temperature

從圖5-7和表2中可以看出，可以看出，當環(huán)境溫度從0 ℃增加到15 ℃時，浮升力影響程度增大，氣體上升高度增加，同時在空氣中擴散距離也增大；而當環(huán)境溫度為30 ℃時，氣體上升進一步增大。因此由圖5-圖7同一時刻不同溫度天然氣變化趨勢得：隨著環(huán)境溫度的升高，天然氣受到浮升力也隨著增大，上升高度亦增加。

Numerical Simulation of Natural Gas Pipeline Leakage and Diffusion

WANG Jian
（PetroChina Xinjiang Oilfield Company Oil & Gas Storage & Transportation Branch，Xinjiang Changji 831100，China）

To study the diffusion of gas after piercing leakage of buried gas pipeline, combined with the finite volume method, CFD simulation model of the gas pipeline leakage under different leak caliber and different environmental temperature was established based on Gambit 2.4. By using Fluent6.3,gas diffusion rules in soil and air were respectively simulated under different leak caliber (6.35, 25.4, 101.6 mm) and different environmental temperature (0, 10, 30 ℃). The research results show that, with increasing of the leakage hole diameter, gas hazard scope gradually increases. After shutting off the valves at both ends of the leak pipeline, the gas diffusion hazard area gradually decreases. With the increase of environmental temperature, the gas hazard scope in the vertical direction significantly increases, the gas hazard scope on the ground also increases at the same time.

Gas pipeline; Different leak location; Diffusion; Numerical simulation

TQ 018

： A

： 1671-0460（2015）10-2494-03

2015-04-20

王?。?987-），男，遼寧錦州人，助理工程師，碩士研究生，2012年畢業(yè)于遼寧石油化工大學，主要從事石油與天然氣運行工作。E-mail: 759439811@qq.com。