架空天然氣管道泄漏數值模擬

（中石化石油工程設計有限公司, 山東 東營 257026）

架空天然氣管道泄漏數值模擬

高 永 強

（中石化石油工程設計有限公司, 山東 東營 257026）

以計算流體力學軟件為基礎，利用組分輸運模型，建立了天然氣泄漏擴散控制方程，對高含硫架空天然氣管道泄漏數值模擬，研究穩態泄漏和非穩態泄漏兩種情況。分析了風速、重力、泄漏量、工況、輸送壓力等因素對天然氣泄漏后擴散過程的影響，得到了硫化氫在不同工況下的擴散規律及安全區域云圖。結合模擬結果，分析了高含硫天然氣的泄漏擴散規律，得到了不同風速條件對架空天然氣管道泄漏的影響，且其模擬結果可以為石油天然氣行業制定相關應急預案及制定安全規章提供指導意義。

天然氣；泄漏；穩態；非穩態；數值模擬

天然氣是重要的清潔能源，是世界一次能源的三大支柱之一，2007年天然氣占世界一次能源消費總量的23.76%。近年來，在世界石油價格不斷攀升的推動下，天然氣勘探開發的速度明顯加快，其產量年均增長率為3.1%～3.4%，是石油產量增長率的3倍[1]。目前，世界管道總長度已超過了鐵路總里程，達到230多萬公里，其中輸氣管道占近60%[2]。對于長距離輸送油氣，管輸是最經濟最高效的方式。但考慮到油氣的易爆、易燃和有毒等特點，發生事故易引起爆炸，污染和中毒等嚴重后果，尤其發生在人口稠密區，同時還會造成惡劣的社會和政治影響[3]。所以，含硫天然氣管道泄漏擴散模擬具有現實意義。

1 天然氣泄漏擴散理論及硫化氫氣體危害

1.1 天然氣泄漏量

孔口泄漏和裂縫泄漏是天然氣管道泄漏主要的兩種方式，由于氣體具有可壓縮性，致使管道泄漏口附近一定范圍之內的較長一段時期仍處于高危險區。泄漏口大小決定了天然氣管道泄漏事故危害等級，在分析過去的事故實例中發現，小孔徑泄漏事故比大孔徑泄漏事故發生的概率大得多。對于裂縫等不規則開孔的微小泄漏源可將其當量為孔口出流情況進行分析[4]。

假設天然氣以小孔出流方式泄漏，且泄漏過程為絕熱擴散過程，則伯努利方程和絕熱方程可得到下列計算式：

實際上由于孔口摩阻損失的存在，泄漏速度往往小于理論值，需要對流速進行修正，故得到的泄漏速度如下：

考慮到泄漏孔形態的多樣性，這里引入流出系數，對計算公式進行修訂：

式中：Φ—為流速系數；

k—為絕熱指數；

Q—為體積流量。

1.2 硫化氫的危害

硫化氫是一種無色、劇毒、強酸性氣體，一旦高含硫化氫氣井發生井噴失控等造成含硫天然氣泄漏，可能導致災難性的后果。曾經，某油田的一口油井，在試油作業起電纜時發生井噴時空，高濃度硫化氫氣體大量噴出，致使7人死亡，數百人中毒，22.6萬人大疏散[5]。

2 架空管道泄漏擴散數值模擬研究

管線腐蝕穿孔是天然氣管道泄漏最常見的工況之一。文章以含硫天然氣管道泄漏為例，對短時穩定泄漏和非穩定泄漏過程進行數值模擬計算，分析了地表風速、管內壓力對擴散過程影響，得出模擬區域內天然氣云團的擴散特性。由于甲烷和硫化氫的云圖分布規律具有相似性，故只給出了硫化氫質量濃度安全區域云圖。

2.1 架空天然氣管道穩態泄漏分析

2.1.1 穩態泄漏擴散控制方程

（1）物質傳輸方程：

湍流脈動動能方程（K方程）

湍流動能耗散方程

湍流粘度方程：

其中：GK—平均速度梯度造成湍動能k產生項；

GB—浮力造成湍動能k產生項；

YM—可壓縮湍流脈動擴張項；

C1ε，C3ε，Cμ—經驗常數，分別為1.44，1，0.09。

（2）物質擴散方程：

其中：Si—第i種物質速率；

Ji—第i種物質質量擴散速率。

湍流中的質量擴散：

其中：Sct—湍流施密特數。

2.1.2 邊界條件

天然氣入口邊界：velocity-inlet(速度入口)；

左側為風速入口邊界:velocity-inlet(速度入口)；上側和右側邊界：pressure-outlet(壓力出口)。

2.1.3 穩定泄漏擴散數值模擬

以某架空管道為例，該管道管外徑610 mm，計算段管道長間距5 km，起始壓力3.5 MPa，終端壓力3.35 MPa，輸量20×108m3/a，硫化氫含量20mg/L，外界環境溫度15 ℃，天然氣溫度26 ℃，泄漏口直徑10 cm，泄漏處位于管段中心，泄漏點處于人口稠密區。模擬區域定為100 m×100 m。經折算，天然氣泄漏速度為213.42 m/s。

2.1.4 結果分析

由圖1分析可知：在泄漏口流速不變的情況下，隨著地表風速的不斷增大，周圍硫化氫分布規律不同。一旦發生泄漏后，高含硫天然氣會以高壓射流的方式射向高空，近地面除泄漏口附近氣體濃度較高外，一定范圍內非危險區。隨著風速的增大，泄漏口處硫化氫不斷向順風側飄移，氣體射出高度降低，受空氣阻力、重力和氣體濃度影響，導致泄漏氣體在下風側一定位置處形成高濃云團，近地面處危險區域增大。因此，在建立應急預案時，應予以考慮，提供相應的解決方案。

圖1 天然氣以213.42 m/s速度泄漏，H2S質量濃度安全區域云圖Fig.1 H2S concentration cloud security area during gas leak at 213.42 m/s speed

2.2 架空天然氣管道非穩定泄漏分析

2.2.1 非穩態泄漏擴散控制方程

（1）物質傳輸方程：

湍流脈動動能方程（K方程）

湍流動能耗散方程

湍流粘度：

其中：GK——平均速度梯度造成湍動能k產生項；

GB——浮力造成湍動能k產生項；

YM——可壓縮湍流脈動擴張項；

C1ε，C3ε，Cμ—經驗常數，分別為1.44，1，0.09。

（2）物質擴散方程：

其中： Si—第i種物質擴散速率；

Ji—第i種物質質量擴散速率。

湍流中的質量擴散：

式中：Sct為湍流施密特數。

2.2.2 邊界條件

天然氣入口邊界：velocity-inlet(速度入口)；左側為風速入口邊界：velocity-inlet(速度入口)；上側和右側邊界：pressure-outlet(壓力出口)。

2.2.3 數值模擬結果

本小節所研究的天然氣管道泄漏過程為管道穿孔后兩端氣源沒有截斷，管道繼續輸送天然氣。由于輸送壓力大，天然氣泄漏湍動能大，且實際中輸送會有壓力脈動，認為泄漏時為非穩態泄漏。研究含硫天然氣隨時間推移的泄漏規律。并分析了風速對天然氣管道非穩定泄漏擴散過程影響。泄漏速速仍為213.42 m/s。

2.2.4 結果分析

圖2 風速0 m/s，泄漏1 s、5 s、10 s，H2S質量濃度分布安全云圖Fig.2 H2S concentration distribution security cloud when wind speed 0 m/s, leak 1 s, 5 s, 10 s

圖3 風速2 m/s，泄漏1 s、5 s、10 s，H2S質量濃度分布安全云圖Fig.3 H2S concentration distribution security cloud when wind speed 2 m/s, leak 1 s, 5 s, 10 s

圖4 風速5 m/s，泄漏1 s、5 s、10 s，H2S質量濃度分布安全云圖Fig.4 H2S concentration distribution security cloud when wind speed 5 m/s, leak 1 s, 5 s, 10 s

由云圖2-4分析可知：天然氣管道泄漏后，如果沒有及時關閉閥門，天然氣繼續流動，則天然氣的非穩態泄漏過程受風速影響較大。當地表無風時，泄漏的天然氣會由泄漏口向上擴散，泄漏1 s時，天然氣氣團達到20 m左右高度；泄漏5 s時，天然氣氣團達到65 m左右高度；泄漏10 s時，氣團已經超過了100 m高度范圍。隨著時間的推移，擴散范圍不斷擴大，而地面附近除了泄漏口附近一定區域范圍，沒有危險區域。當風速達到2 m/s時，泄漏1 s時，天然氣氣團高度達到將近20 m，向下風側有一個小角度的傾斜；泄漏時間達到5 s時，氣團高度達到30 m左右，向下風側進一步擴散、傾斜；當時間達到10 s時，氣團高度已經沒有明顯的增大了，泄漏天然氣進一步向下風側擴散、傾斜，范圍已經超過了50 m。當風速達到5 m/s時，泄漏時間為1 s時，天然氣氣團擴散高度將近20 m，傾斜程度相對于風速2 m/s時幅度變大；當泄漏時間為5 s時，泄漏氣團達到25 m左右高度，天然氣在下風側進一步擴散、傾斜，已經超出了50 m的區域；當泄漏時間為10 s時，氣團高度已經沒有進一步上升了，而是繼續向下風側飄移、擴散。有風速時，危險區域主要集中在近地面下風側區域。針對模擬結果，可以看到，若要采取救援，在無風時，應帶好防毒面具，實施近地救援；當有風速時，應在上風側進行救援，以降低危險隱患。

3 結論與建議

3.1 總結

針對架空含硫天然氣管道發生穿孔泄漏后的擴散過程進行系統分析[6]，建立了泄漏擴散數學模型，并借助計算流體力學軟件進行求解模擬計算，分析了不同工況下天然氣泄漏擴散規律，主要結論如下:

（1）系統分析了風速對天然氣泄漏擴散的影響規律，結合多相流體力學、傳熱學、工程熱力學等相關理論知識，建立了架空天然氣管道泄漏后擴散數學模型，應用有限體積法對模型偏微分方程組離散、轉化為總體有限元方程，同時加上邊界條件進行求解方程。

（2）利用計算流體力學軟件對架空天然氣管道泄漏過程模擬計算，給出了硫化氫質量濃度安全區域云圖，計算結果符合擴散、著火爆炸理論，進一步驗證了所建立的泄漏擴散模型是正確的。

（3）架空天然氣管道泄漏后形成高壓射流，進入大氣中開始迅速膨脹，濃度和速度逐漸減小，為高速噴射擴散過程，且天然氣擴散規律受到氣候條件因素的影響較大。

（4）影響天然氣擴散規律的因素主要包括輸送壓力、泄漏速度、風速、周圍環境大氣溫度、管道沿線地形、泄漏處位置等，不同因素對含硫天然氣擴散影響程度不同。

（5）在重力作用下，各個組分會出現分層現象，原因是多組分天然氣發生泄漏而擴散時，重組分會發生下沉，含硫天然氣泄漏后，硫化氫氣體會在低洼地區發生沉積，造成人員中毒傷亡，使得危險性增大。

（6）為了有效防止天然氣泄漏后引起中毒、爆炸事故的發生，可以通過控制或人為改變某些因素，例如從擴散和爆炸的影響因素著手，從而降低其危險性。

3.2 建議

天然氣管道泄漏擴散和爆炸影響因素眾多，機理復雜，實驗研究還尚存一定的困難，為了進一步的研究其擴散規律和爆炸機理，提出以下幾點建議：

（1）本文所建立的天然氣管道泄漏擴散模型與實際工況下的泄漏擴散過程有一定的差距，因為它是在一定的假設條件下建立的，為了進一步提高模型的普遍適用性，則需要結合大量的現場數據和實驗做進一步的分析和研究，。

（2）由于天然氣是一種多組分混合物，其泄漏后與一般的工業氣體有較大的區別。天然氣管道發生泄漏后，各種因素疊加，錯綜復雜，受到氣候和地形等條件的綜合影響。雖然可以借助計算流體軟件進行計算，但是受到物理模型抽象和假定一些因素不變所約束，還不能完全分析天然氣管道的實際狀況及多種因素綜合作用下的擴散規律。

（3）可以根據本文的基礎條件，建立三維模型，進一步對天然氣管道泄漏過程進行模擬，在實踐中，還需要多多觀察、記錄，確定影響泄漏擴散的主要因素，忽略次要因素，進一步完善天然氣的泄漏擴散模擬，要想完美的模擬出實際的擴散情況，還有很多的工作需要做。

［1］賀永德．天然氣應用技術手冊[M]．北京：化學工業出版社，2009.

［2］張其敏，孟江．油氣管道輸送技術[M]．北京：中國石化出版社，2008.

［3］楊莜蘅.油氣管道安全工程[M]．北京：中國石化出版社，2005.

［4］胡棟，龍屬川，吳傳謙．可爆性氣體爆炸極限和爆燃轉變成爆轟的研究[J]．爆炸與沖擊，1959，9(3)：266-275.

［5］易俊，王以朗，朱俊，謝代安．天然氣采輸作業硫化氫防護[M]．重慶：西南師范大學出版社，2010.

［6］朱紅鈞，林元華，謝龍漢．FLUENT流體分析及仿真實用教程[M]．北京：人民郵電出版社，2010.

Numerical Simulation on Leakage of Overhead Natural Gas Pipeline

GAO Yong-qiang
（Sinopec Petroleum Engineering Design Co., Ltd., Shandong Dongying 257026，China）

Based on the CFD software, making use of components transport models, the control equation of leakage and diffusion of natural gas was established, numerical simulation for leakage of overhead high sulfur natural gas pipeline was carried out, two conditions including steady-state leakage and unsteady-state leakage were studied. Effect of wind speed, gravity, leakage amount, working conditions, transport pressure and other factors on the diffusion process of natural gas after the leakage was analyzed, diffusion law and safe regional cloud of hydrogen sulfide in different working conditions were obtained. Combined with the simulation results, the diffusion law of high sulfur natural gas was analyzed, and effect of different wind speed on the leakage of natural gas was obtained, the simulation results can provide guidance on formulating related emergency plan and safe practice for oil and natural gas industry.

Natural gas; Leakage; Steady state; Unsteady state；Numerical simulation

TQ 018

： A

： 1671-0460（2014）04-0655-04

2013-10-14

高永強（1981-），男，河北石家莊人，研究方向：油氣儲運設計相關工作。E-mail：gaoyongqiang77@163.com。