程猛猛 吳 明 趙 玲 劉 柯 邵大偉 楊明悅

(1.遼寧石油化工大學石油天然氣工程學院)(2.沈陽市特種設備檢測研究院)

安全與環(huán)保

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

程猛猛1吳 明1趙 玲1劉 柯2邵大偉2楊明悅1

(1.遼寧石油化工大學石油天然氣工程學院)(2.沈陽市特種設備檢測研究院)

針對城市埋地天然氣管道穿孔泄漏擴散問題，結合有限容積法，利用Gambit 2.4建立了天然氣管道不同泄漏位置的CFD仿真模型，利用Fluent 6.3分別對天然氣管道上部、下部及背風側3種泄漏工況下，氣體在土壤中和空氣中的擴散規(guī)律進行了數(shù)值模擬。研究結果表明，下部泄漏在土壤和空氣中的危險范圍最大，關閉泄漏管段兩端閥門以后，氣體擴散危害范圍逐漸變小。研究結果為城市埋地天然氣管道泄漏事故現(xiàn)場人員疏散及安全搶修提供了理論依據(jù)。

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

隨著城市的發(fā)展和人民生活水平的提高，城市燃氣事業(yè)得到迅速的發(fā)展，燃氣作為重要的能源形式之一，應用越來越廣泛。城市燃氣管道應采用埋地方式鋪設，但埋地管道可能在人為破壞或自然災害時發(fā)生泄漏事故[1]。燃氣具有易燃、易爆的性質，燃氣泄漏爆炸事故一般會造成大量人員傷亡和嚴重的財產(chǎn)損失，具有較大的危害性。因此，進行天然氣管道泄漏擴散數(shù)值模擬具有重大意義。

本文針對城市埋地天然氣管道，考慮風速隨高度變化和泄漏后閥門自動關閉或者人工關閉工況，利用C語言編寫與Fluent軟件兼容的UDF計算程序，修正風速隨高度的變化和閥門關閉后氣體泄漏率隨時間的變化[2-3]。由于腐蝕、施工等不同的損傷過程，造成泄漏口在管道上的位置會發(fā)生變化，故選取管道泄漏口分別位于管道上部、下部和背風側3種工況，建立泄漏模型來模擬計算，得出了天然氣管道不同損傷下的天然氣泄漏擴散規(guī)律，為安全施工、預防事故及制定天然氣管路泄漏應急預案提供了理論依據(jù)。

1 模型建立

1.1數(shù)學模型的建立

天然氣在泄漏過程中遵循質量守恒定律、動量守恒定律和能量守恒定律[4],其控制方程如下：

(1) 連續(xù)性方程：

(1)

式中，ρ為流體密度，kg/m3；t為時間，s；ui為(x、y、z)3個方向上的速度，m/s。

(2) 動量方程：

(2)

式中，f為單位質量力矢量，m/s2；t為時間，s；u為速度，m/s；μ為動力黏度，Pa·s；p為流體微元上的壓力，Pa。

(3) 能量方程

三組學生網(wǎng)上收集“汶川地震17歲少年馬健用雙手救同學”的故事，談談勇氣、堅毅、信心、意志和干勁是人的精神狀態(tài)，如何理解“精神不是萬能的，但沒有精神是萬萬不能的”？

(3)

式中，E為流體微元的總能，J；hj為組分j的焓，J/kg；keff為有效熱傳導系數(shù)；Jj為組分j的擴散通量，mol/(m2·s)；Sh為體積熱源項，kg·J/(m3·s)。

1.2泄漏擴散控制方程

天然氣的擴散可認為是多組分氣體相互作用而形成的湍流。綜合考慮計算條件和工程要求，采用耦合隱式求解器，其基本方程為定常Reynolds時均方程[5]。由于高壓管道泄漏初速度很大，故選用湍流模型中的可實現(xiàn)k-ε模型。

湍流脈動動能方程(k方程)：

ρε-Ym

(4)

湍流動能耗散方程(ε方程)：

(5)

(6)

式中，Gk為平均速度梯度引起的湍流動能產(chǎn)生項；Gb為浮力引起的湍流動能產(chǎn)生項；Ym為可壓縮湍流中脈動擴張貢獻；μt為湍動黏度，Pa·s；σk、σε分別為k方程和ε方程的湍流Prandtl數(shù)；C1ε=1.44、C2ε=1.92、C3ε=0.09、Cμ=0.084 5為經(jīng)驗常數(shù)。

1.3物理模型

燃氣管道泄漏時，由于外界壓力較小，泄漏氣體迅速膨脹，氣體擴散的過程極為復雜。為了研究管道泄漏后氣體擴散規(guī)律，確定危害范圍，以城市天然氣管道泄漏為例，模擬泄漏過程及天然氣流場變化。天然氣組分中以甲烷(CH4)為主(體積分數(shù)為99.88%)，還包含有少量其他氣體，故以CH4代替天然氣進行模擬計算。假設發(fā)生泄漏后從管內到接觸大氣的過程中，摩阻減速的作用不大。天然氣的傷害下限為1%(φ)，爆炸極限范圍為5%(φ)～15%(φ)[6]，故取天然氣的體積分數(shù)為1%、5%和15%的區(qū)域為評定泄漏氣體危害分界濃度。

某天然氣管道運行壓力為0.2 MPa，根據(jù)GB 50028-2006《城市燃氣設計規(guī)范》，屬于中壓輸氣管道，管道直徑325 mm、壁厚8 mm，管道埋深1.7 m，氣體溫度為20 ℃。由于所選天然氣管道在城鎮(zhèn)人口密集地區(qū)，考慮管道所處的實際位置，取建筑物高度為20 m，管道中心距建筑物距離為3 m[7]，用較大空間 (80 m×80 m)進行仿真模擬計算，具體模型及邊界條件如圖1所示。劃分網(wǎng)格時，由于管道泄漏口附近氣體各項參數(shù)變化梯度較大，為滿足計算精度，對此區(qū)域進行網(wǎng)格加密。本文以管道穿孔泄漏為例，根據(jù)API 581-2000《基于風險的檢驗基礎資源文件(Risk-Based Inspection Base Resource Document)》選取管道泄漏口徑為25.4 mm的泄漏模型來模擬計算。

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

設定模擬的環(huán)境溫度和氣體溫度均為300 K，距離地面高度10 m處風速為4 m/s，利用UDF程序對風速隨高度變化進行修正。由于動力因素和熱力因素，在近地層中，風速會隨高度發(fā)生變化，并且服從普朗特亂流中經(jīng)驗理論公式[8]。求高度為H處的風速的公式為：

v=v0(H/H0)n

(7)

式中：v為高度在H處的風速，m/s；v0為高度在H0處的風速，m/s；H0一般為10 m；n為風速隨高度的變化系數(shù)，最常用的為0.142。

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

(8)

(亞音速泄漏)(9)

Qm(t)=Qm0·[F(t)](k+1)/(k-1)

(10)

(11)

(12)

m0=L·Ac·ρ

(13)

(14)

式中，m0為管道內氣體質量，kg；Qm為質量流量，kg/s；k為比熱容；t為關閉閥門后泄漏時間，s；Ac為管道橫截面積，m2；Aor為泄漏口面積，m2；M為氣體摩爾質量，kg/kmol；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。

泄漏口設為質量流量進口；左側設置為風入口，設置為速度進口；空間上邊界設為壓力出口；土壤兩側設為壓力出口，坡面和土壤層表面設為多孔介質，土壤物性為：密度1 500 kg/m3、導熱系數(shù)1.5 W/ (m·K)、比熱容2 200 J/(kg·K )、孔隙度0.43[11]、土壤平均直徑0.2 mm，根據(jù)文獻[12]計算出滲透率及慣性阻力系數(shù)。管道、建筑物及地面的邊界條件設為壁面。穩(wěn)態(tài)泄漏120 s內質量流量進口設置為1 077.95 kg/(m2·s)，關閉閥門后，根據(jù)質量流量隨時間變化擬合出公式(15)：

G=1 047.4×e-0.001 3t

(15)

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

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

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

當管道受到不同的損傷，泄漏口在管道上出現(xiàn)的位置是不同的，而對于埋地燃氣管道泄漏口位置不同，氣體擴散的規(guī)律也不盡相同。根據(jù)泄漏口在管道上的不同位置修改幾何模型，上部泄漏和背風側泄漏地面距離幾何模型下邊界為5 m，下部泄漏工況下地面距離幾何模型下邊界為10 m。以泄漏口直徑為25.4 mm為例說明，氣體擴散云圖如圖2～圖4所示，表1更詳細地說明了CH4在土壤中和空氣中擴散的范圍。

表1 不同泄漏位置天然氣危害范圍Table1 Gashazardareaofdifferentleakageposition泄漏位置上部泄漏下部泄漏背風側泄漏泄漏時間/s120434120434120434傷害下限距離/m342838383329.6傷害下限距地面高度/m>70>63>65>57.8>70>62爆炸下限距離/m181325.520.518.815.9爆炸下限距地面高度/m6347584357.847土壤中爆炸下限距離/m2016.53326.520.817.7土壤中傷害下限距離/m>3829.8>38>3835.531

由圖2～圖4可以看出，管道泄漏位置不同，泄漏天然氣在土壤和空氣中擴散的范圍也不同，其中管道下部泄漏在土壤中擴散的范圍最大，并且在泄漏率接近零時，甲烷的危害范圍最大。從表1可以看出，當泄漏120 s 時，上部泄漏情況下，傷害下限和爆炸下限在下風向貼近地面處距離泄漏口分別為34 m和18 m，擴散高度分別為70 m以上和63 m，在土壤中沿下風向擴散至距泄漏口38 m和20 m；下部泄漏情況，傷害下限在下風向貼近地面處和土壤中距離泄漏口均超過38 m，擴散高度為65 m，爆炸下限在下風向貼近地面處和土壤中至泄露口距離分別為25.5 m和33 m，擴散高度為58 m；背風側泄漏情況，傷害下限和爆炸下限在下風向貼近地面處距離泄漏口分別為33 m和18.8 m，擴散高度分別為70 m以上和57.8 m，在土壤中沿下風向擴散至距泄漏口35.5 m和20.8 m。由此可知，當泄漏口位于埋地管道下部時，甲烷首先向下部土壤中擴散，在土壤中擴散的范圍最大，在地面附近擴散距離最遠，并且泄漏的CH4不易擴散。

結合圖2～圖4及表1可以看出，關閉閥門以后，泄漏氣體質量流量逐漸變小，由于管道壓力較低，當管道內壓力與環(huán)境壓力相等時，泄漏的質量流量接近于零，受風和濃度差擴散作用CH4危險區(qū)域減小，例如管道上部泄漏工況下，當泄漏120 s 時，上部泄漏情況，傷害下限和爆炸下限在下風向貼近地面處距離泄漏口分別為34 m和18 m，擴散高度分別為70 m以上和63 m，在土壤中沿下風向擴散至距泄漏口38 m和20 m；在約434 s泄漏口處質量流量接近于零，故氣體危險濃度范圍也在變小，傷害下限和爆炸下限在下風向貼近地面處距離泄漏口分別為28 m和13 m，擴散高度分別為63 m以上和47 m，在土壤中沿下風向擴散至距泄漏口29.8 m和16.5 m。由此可知，埋地燃氣管道一旦發(fā)生泄漏，盡快關閉泄漏管段兩端的閥門，能夠減小CH4的危險濃度范圍。

3 結 論

城市埋地天然氣管道發(fā)生泄漏后，在土壤中由于濃度差作用，CH4在管道周圍呈發(fā)散狀擴散。CH4進入空氣以后，在浮力和風作用下，分別向高空和下風向迅速擴散，由于建筑物的阻礙作用，在建筑物的迎風側形成高濃度區(qū)；在土壤中氣體繞過建筑物擴散到建筑物背風側，在浮力與濃度差的驅動下，天然氣從土壤中擴散到空氣中，并在建筑物背風側形成高濃度危險區(qū)。對天然氣管道上部、下部及背風側3種泄漏工況下，下部泄漏在土壤和空氣中的危險范圍最大，關閉泄漏管段兩端閥門以后，氣體擴散危害范圍逐漸變小。

[1] 馮慶善,陳健峰,艾慕陽,等.管道完整性管理在應對地震災害中的應用[J].石油學報, 2010,31(1): 139-143.

[2] 霍春勇,董玉華,余大濤,等.長輸管線氣體泄漏率的計算方法研究[J].石油學報,2004,25(1):101-105.

[3] 張寶柱,王秀紅,高軍,等.障礙物不同坡度對H2S擴散影響研究[J].石油與天然氣化工,2012,41(2):229-231.

[4] 朱紅鈞,林元華,謝龍漢.FLUENT流體分析及仿真實用教程[M].北京:人民郵電出版社,2010.

[5] 劉勇峰,吳明,趙玲,等.城市天然氣管道泄漏數(shù)值模擬[J].工業(yè)安全與環(huán)保,2012,38(8):51-53.

[6] Rae M, Lee J H. Acoustic emission technique for pipeline leak detection [J]. Key Engineering Materials, 2000, 186(4): 555-892.

[7] 中華人民共和國建設部.GB 50028-2006城市燃氣設計規(guī)范[S].北京:中國建筑工業(yè)出版社,2006-11-01.

[8] 潘陽，許國良，中山顕，等.計算傳熱學理論及其在擴孔介質中的應用[M].北京：科學出版社，2011.

[9] 張舜德,高文元,王現(xiàn)青,等.風場風速特性的研究[J].機械設計與制造,2010(8):150-151.

[10] 張紅兵,李長俊.長距離輸氣管線泄漏監(jiān)測技術研究[J].石油與天然氣化工,2005,34(2):146-147.

[11] 李朝陽，馬貴陽，劉亮.埋地輸油管道泄漏油品擴散模擬[J].油氣儲運，2011，30(9)：674-676.

[13] 李柯.天然氣長輸管道泄漏工況數(shù)值模擬[J].油氣儲運,http://www.cnki.net/kcms/detail/13.1093.TE.20130802.0958.007.html.

[14] 楊昭,賴建波,韓金麗. 天然氣管道孔口泄漏危險域的研究[J]. 天然氣工業(yè), 2006, 26(11):156-159.

Numericalsimulationofurbanburiedgaspipelineleakageanddiffusion

ChengMengmeng1，WuMing1，ZhaoLing1，LiuKe2，ShaoDawei2，YangMingyue1

(1.CollegeofPetroleumEngineering，LiaoningShihuaUniversity，F(xiàn)ushun113001，Liaoning,China；2.ShenyangInstituteofSpecialEquipmentInspection&Research，Shenyang110000,Liaoning，China)

To study the diffusing of gas piercing leakage in urban buried gas pipeline, combining the finite volume method, CFD simulation models are established by Gambit 2.4 at different leakage locations. The diffusion rule of gas in soil and air are simulated by Fluent 6.3 under the three leakage conditions of the upper, lower and leeward side of the gas pipeline. The research result shows that the lower leakage in soil and air is the largest diffusion hazard area. If the valves at both ends of the pipeline are shut off, the gas diffusion hazard area will gradually get smaller. The studies can also provide the theoretical foundation for the personnel evacuation and security maintenance of pipelines in the scene of natural gas release.

gas pipeline, leakage, diffusion, numerical simulation

TE88

：ADOI: 10.3969/j.issn.1007-3426.2014.01.020

2013-10-21；

2013-11-21；編輯：鐘國利

國家質檢總局科研項目“質檢公益性行業(yè)科研專項項目”(201010025)。

程猛猛(1987-)，男，河南商丘人，遼寧石油化工大學在讀碩士研究生，主要從事輸氣管道技術方面的科研工作。地址：(113001)遼寧省撫順市望花區(qū)丹東路西段一號。E-mail: cm07131042@126.com