馬貴陽(yáng)，王　新*，李　丹，劉兵兵，黃騰龍

埋深對(duì)平坦地區(qū)天然氣管道泄漏的影響研究

馬貴陽(yáng)1，王新1*，李丹2，劉兵兵3，黃騰龍3

1.遼寧石油化工大學(xué)石油天然氣工程學(xué)院，遼寧 撫順 113001；2.中國(guó)石油天然氣管道工程有限公司，河北 廊坊 065000；3.中國(guó)石油天然氣管道局，河北 廊坊 065000

針對(duì)三維埋地輸氣管道泄漏擴(kuò)散問題，對(duì)不同埋深的平坦地區(qū)天然氣管道泄漏情況進(jìn)行數(shù)值模擬。根據(jù)單一變量原則在相同氣候條件下，對(duì)于不同工況只針對(duì)埋深作為單一變量，對(duì)忽略埋深的準(zhǔn)確性進(jìn)行論證，并分別研究了埋深為1.4 m和2.0 m工況埋深對(duì)地下、地表、和空氣中泄漏的影響。研究結(jié)果表明：埋深對(duì)泄漏的影響非常大，忽略埋深的工況與埋深為1.4 m和2.0 m的工況相比所得出的各項(xiàng)結(jié)果都有很大的誤差，忽略埋深是不準(zhǔn)確的。埋深與擴(kuò)散范圍、泄漏速度、質(zhì)量分?jǐn)?shù)、高濃度范圍成反比，埋深越小擴(kuò)散范圍、泄漏速度、質(zhì)量分?jǐn)?shù)、高濃度范圍越大。

天然氣泄漏；三維埋地管道；危險(xiǎn)區(qū)域；埋深；數(shù)值模擬

網(wǎng)絡(luò)出版地址：http：//www.cnki.net/kcms/detail/51.1718.TE.20150330.1719.014.html

馬貴陽(yáng)，王 新，李 丹，等.埋深對(duì)平坦地區(qū)天然氣管道泄漏的影響研究［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2015，37（2）：165-172.

Ma Guiyang，Wang Xin，Li Dan，et al.Influence of Burial Depth for Natural Gas Pipeline Leakage in Flat Country［J］.Journal of Southwest Petroleum University：Science&Technology Edition，2015，37（2）：165-172.

引言

管道運(yùn)輸是天然氣最主要的運(yùn)輸方式，而且多以埋地管道為主。由于管道所受到的腐蝕作用，輸氣管網(wǎng)本身的一些問題（主要的是質(zhì)量問題），人類的破壞，以及來自地震的毀滅性的破壞，天然氣輸氣管網(wǎng)的泄漏問題有時(shí)候是不可避免的。一旦發(fā)生泄漏事故，考慮到天然氣的燃燒擴(kuò)散特點(diǎn)（易爆炸性和易擴(kuò)散性），將給人類和人類賴以生存的環(huán)境帶來嚴(yán)重的影響，因此對(duì)于天然氣輸氣管網(wǎng)泄漏后的泄漏范圍（危險(xiǎn)區(qū)域范圍和安全區(qū)域范圍）的研究計(jì)算顯得尤為重要［15］。國(guó)內(nèi)外學(xué)者已對(duì)天然氣泄漏后擴(kuò)散的一般規(guī)律和模擬結(jié)果進(jìn)行分析。得出一些編程相對(duì)簡(jiǎn)單、計(jì)算快的半經(jīng)驗(yàn)半理論計(jì)算公式，但其結(jié)果的準(zhǔn)確性很大程度上依賴于有害氣體的擴(kuò)散條件，如泄漏類型、自然風(fēng)向和強(qiáng)度、大氣壓力梯度分布、周圍地形建筑等與模型假設(shè)的接近程度［68］。國(guó)內(nèi)學(xué)者朱紅鈞等對(duì)平坦地區(qū)天然氣泄漏進(jìn)行二維數(shù)值模擬［9］，得出了泄漏后一定范圍內(nèi)的安全區(qū)域。馬存棟研究了高壓天然氣輸氣管道發(fā)生泄漏事故時(shí)氣體移動(dòng)擴(kuò)散的規(guī)律［10］。文獻(xiàn)［11］采用特征線法進(jìn)行天然氣管道泄漏模擬，采用儲(chǔ)罐和管道的零維模型來建立，由于孔洞導(dǎo)致的壓力下降和氣體泄漏過程的模型，這些儲(chǔ)罐和管道的模型適合小孔洞或完全斷裂工況下的仿真。前人對(duì)天然氣泄漏問題的研究適用條件有局限性，且多為二維情況，雖有少數(shù)學(xué)者對(duì)三維問題進(jìn)行研究，但所建立的模型只有地上的部分，無埋地管道模型，研究的結(jié)果都是在忽略了埋深的條件得出的，其結(jié)果不夠精確，與實(shí)際情況相差較大。

對(duì)埋地天然氣管道泄漏進(jìn)行三維數(shù)值模擬，考慮到天然氣泄漏后會(huì)經(jīng)過土壤而泄漏到大氣當(dāng)中，而土壤為多孔介質(zhì)，在多孔介質(zhì)的內(nèi)部，流體的流動(dòng)以及相變過程很復(fù)雜，而有限容積法是處理這類問題較為常用的方法［1217］，采用有限容積法為原理的計(jì)算流體力學(xué)軟件研究埋深地泄漏擴(kuò)散的影響，并只針對(duì)埋深作為單一變量，其他氣候條件地理?xiàng)l件均相同。首先對(duì)以往忽略埋深的工況進(jìn)行研究，對(duì)其準(zhǔn)確與否進(jìn)行討論，并對(duì)不同埋深管道泄漏進(jìn)行模擬，分別研究了埋深對(duì)地下、地表和空氣中的泄漏情況，且定量證明了埋深對(duì)泄漏的重要性。為事故后確定泄漏范圍和危險(xiǎn)區(qū)域提供了理論依據(jù)。

1　計(jì)算模型

流體流動(dòng)遵循動(dòng)量守恒定律，忽略溫度對(duì)流動(dòng)的影響。模型主要包含質(zhì)量守恒方程、動(dòng)量守恒方程和組分守恒方程，其中，動(dòng)量守恒方程為

式中：i=x，y，z；

Sx=0；

Sy=0

Sz=gβT（T0?T）

ρ—流體密度，kg/m3；

ux、uy、uz—x、y，z向的速度分量，m/s；

t—時(shí)間，s；

u—速度矢量，m/s；

p—流體壓強(qiáng)，Pa；

μeff—有效黏度，Pa·s；

g—重力加速度，g=9.8 m/s2；

βT—體積膨脹系數(shù)，無因次；

T0—參考溫度，K；

T—溫度，K。

天然氣輸氣管道發(fā)生泄漏事故規(guī)模的大小可以按照泄漏口的大小進(jìn)行劃分，而且孔徑較小的泄漏事故發(fā)生的概率較大一些，孔徑較大的泄漏事故發(fā)生的概率則相對(duì)較小。假設(shè)所研究的泄漏過程為泄漏口孔徑為較小的情況，而且從孔徑中排出氣體的過程看做是絕熱過程，認(rèn)為天然氣為理想氣體。由伯努利方程以及絕熱方程可以得到，泄漏口的氣體流速

式中：

v0—泄漏速度，m/s；

φ—流速系數(shù)，無因次，代表實(shí)際流速與理論流速之比，一般為0.97～0.98；

k—絕熱指數(shù)，無因次，對(duì)于天然氣一般取1.3；

R—?dú)怏w常數(shù)，R=8.314 J/（mol·K）；

T1—?dú)怏w溫度，K；

p0—環(huán)境壓力，Pa；

p1—?dú)怏w泄漏前壓力（管道內(nèi)介質(zhì)壓力），Pa。

由于孔口形狀不同，引入流量系數(shù)，則體積流量為

式中：

Q—體積流量，m3/s；

CD—流量系數(shù)，無因次；

A—泄漏孔面積，m2；

M—分子量，天然氣取0.017 kg/mol；

σ—速度系數(shù)，無因次。

式（3）的流量系數(shù)CD，亦稱泄漏系數(shù)。一般取值在0.6～1.0，對(duì)圓形孔，CD=1.0；三角孔，CD=0.95；長(zhǎng)形孔，CD=0.90；孔口為內(nèi)層腐蝕形成的減縮孔，0.9＜CD＜1.0；孔口為外力機(jī)械損傷形成漸擴(kuò)孔，0.6＜CD＜0.9。本文選圓形孔，CD=1.0。

泄漏口的氣體流速與流體的流動(dòng)狀態(tài)有很大關(guān)系，所以在研究泄漏過程前要先確定氣體的流動(dòng)狀態(tài)（音速或亞音速）。在長(zhǎng)輸氣體管道中由于介質(zhì)壓力較高，氣體處于流速臨界狀態(tài)，流動(dòng)屬于音速流動(dòng)，即，此時(shí)速度系數(shù)為

氣體從泄漏口出來，可以按射流處理，沿射流中心方向上的速度可表示為

式中：vm—圓斷面射流的軸心速度，m/s；dl—泄漏孔半徑，m；a—紊流系數(shù)，無因次；s—到射流出口的距離，m。

由于土壤對(duì)泄漏氣體的阻礙作用，從泄漏口到地表這段距離，氣體的泄漏速度可以看做是分層遞減的，為簡(jiǎn)化分析計(jì)算，將遞減過程近似為均勻遞減，則當(dāng)氣體泄漏到地表時(shí)會(huì)形成一個(gè)等半徑的圓形的泄漏區(qū)域，此時(shí)的氣體泄漏速度才是更為準(zhǔn)確的，為簡(jiǎn)化計(jì)算，實(shí)際泄漏速度應(yīng)該為射流乘以一個(gè)小于1的速度調(diào)整常數(shù)α

式中：α—沿程速度修正系數(shù)，無因次，與通過土壤性質(zhì)有關(guān)。

根據(jù)天然氣的泄漏特點(diǎn)，天然氣通過管道泄漏口射出管外，對(duì)于埋地天然氣輸氣管道，泄漏剛開始時(shí)，氣體射入土壤中，氣體在流動(dòng)過程中會(huì)受到來自土壤的孔隙阻力以及毛細(xì)管壓力，其泄漏速度將會(huì)在一定的程度上減小。由式（5）可知，地表速度的大小在不計(jì)其他耗省的工況中受射流距離的影響很大；由式（6），對(duì)于埋地天然氣泄漏射流距離即管道埋深對(duì)地表泄漏速度的影響更大，而地表泄漏速度的大小將決定氣體擴(kuò)散到大氣中的初速度，進(jìn)而影響擴(kuò)散范圍。

2　計(jì)算實(shí)例

2.1實(shí)例參數(shù)

選取平坦地區(qū)泄漏工況，有埋深的情況模擬區(qū)域尺寸為：40.0 m×40.0 m×40.0 m，其中y方向范圍為（-10，30），即空氣域高度為30 m；不計(jì)埋深的時(shí)模擬區(qū)域尺寸為40.0 m×30.0 m×40.0 m，其中y方向范圍為與前面一致。定義管道與地表的最近距離為埋深。管道外徑1.0 m，泄漏口直徑0.3 m，管內(nèi)和環(huán)境溫度300 K，土壤含水率9%。漏口氣體流速為5 m/s，泄漏口方向向上，天然氣中甲烷和硫化氫的濃度分別為93.5%和6.5%。網(wǎng)格模型如圖1所示，網(wǎng)格中紅色為空氣域，藍(lán)色為土壤域。

圖1　網(wǎng)格模型Fig.1　Mesh model

其中，圖1a為忽略埋深時(shí)仰視網(wǎng)格圖，黑色密集中心位置為泄漏口，紅色密集中心區(qū)域?yàn)樾孤┛谡戏剑瑘D1b、圖1c分別為有埋深時(shí)管道附近網(wǎng)格圖、管道及泄漏口的計(jì)算網(wǎng)格圖。

2.2忽略埋深的不準(zhǔn)確性

同一工況的天然氣泄漏事故，當(dāng)環(huán)境因素和泄漏速度相同時(shí)，管道的埋深對(duì)天然氣的泄漏擴(kuò)散范圍影響很大。先對(duì)忽略埋深的不準(zhǔn)確性進(jìn)行驗(yàn)證。

圖2中給出了泄漏10 s和泄漏60 s時(shí)不同埋深工況地表以上1.0 m高處平面中同樣位置各點(diǎn)的甲烷和硫化氫質(zhì)量分?jǐn)?shù)。對(duì)于泄漏10 s的情況，由圖2a、圖2b可以得出，在平面內(nèi)同樣位置忽略埋深的情況甲烷和硫化氫的質(zhì)量分?jǐn)?shù)都比埋深為1.4 m和2.0 m的情況甲烷的質(zhì)量分?jǐn)?shù)大的多，而埋深為1.4 m和2.0 m的情況甲烷和硫化氫的質(zhì)量分?jǐn)?shù)相差不大。對(duì)于泄漏60 s的情況，由圖2c、圖2d可以得出，泄漏60 s時(shí)忽略埋深的情況甲烷和硫化氫的質(zhì)量分?jǐn)?shù)仍然比埋深為1.4 m和2.0 m的情況大，但大小差距相對(duì)泄漏10 s時(shí)有所下降。

所以，無論是同一時(shí)刻還是不同時(shí)刻甲烷和硫化氫的質(zhì)量分?jǐn)?shù)忽略埋深的情況總是比埋深為1.4 m和2.0 m的情況大，可以得出模擬過程中忽略埋深即忽略埋深所得出的各項(xiàng)結(jié)果都有很大的誤差，忽略埋深是不準(zhǔn)確的，所以下文對(duì)埋深為1.4 m和2.0 m的工況進(jìn)行研究，在考慮埋深的基礎(chǔ)上研究埋深對(duì)泄漏的影響。

圖2　不同埋深工況下，地表以上1.0 m高處平面上同一位置天然氣擴(kuò)散速度Fig.2　The diffusion velocity of natural gas at the same location of the same plane under different burial depth conditions

2.3埋深對(duì)天然氣泄漏擴(kuò)散的影響

為研究埋深對(duì)天然氣泄漏擴(kuò)散的影響，分別研究了地下、地表和空氣中天然氣的泄漏情況。圖3、圖4分別給出了該平面內(nèi)埋深為1.4 m和2.0 m時(shí)不同時(shí)刻天然氣質(zhì)量分?jǐn)?shù)流場(chǎng)圖，圖5為泄漏過程中該平面速度散點(diǎn)圖。

由圖3、圖4可以得出，天然氣的擴(kuò)散范圍在一定時(shí)間內(nèi)與時(shí)間成正比，分別對(duì)比圖3a和圖4a，圖3b和圖4b，圖3c和圖4c可以看出，相同時(shí)刻埋深1.4 m比埋深2.0 m硫化氫的擴(kuò)散范圍大，又由圖5a中可以看出泄漏過程中埋深1.4 m比埋深2.0 m泄漏速度大，尤其是位置0處（泄漏口上方對(duì)應(yīng)的位置），速度越大相同時(shí)間內(nèi)泄漏的天然氣就越多，其質(zhì)量分?jǐn)?shù)就會(huì)越大，驗(yàn)證了圖3、圖4。圖5b與5c中速度相差無幾，說明此時(shí)該平面內(nèi)速度已經(jīng)穩(wěn)定。綜合圖3～圖5可以得出，此平面內(nèi)無論是擴(kuò)散范圍、泄漏速度、還是質(zhì)量分?jǐn)?shù)埋深1.4 m都比埋深2.0 m大，即埋深為擴(kuò)散范圍、泄漏速度、質(zhì)量分?jǐn)?shù)的減函數(shù)，埋深越小擴(kuò)散范圍、泄漏速度、質(zhì)量分?jǐn)?shù)越大。

圖3　埋深1.4 m地表下1.0 m深處平面內(nèi)天然氣泄漏后的流場(chǎng)圖Fig.3　The velocity of natural gas below the surface of the deep 1.0 m under 1.4 m burial depth conditions

圖4　埋深2.0 m地表下1.0 m深處平面內(nèi)天然氣泄漏后的流場(chǎng)圖Fig.4　The velocity of natural gas below the surface of the deep 1.0 m under 2.0 m burial depth conditions

圖5　不同埋深地表以下1.0 m深處平面內(nèi)天然氣泄漏后的速度散點(diǎn)圖Fig.5　The leakage natural gas velocity below the surface of the deep 1.0 m under different burial depth conditions

對(duì)泄漏過程中地表處的天然氣平均速度和質(zhì)量分?jǐn)?shù)進(jìn)行跟蹤監(jiān)測(cè)得到圖6。其中圖6a、圖6b為地表處天然氣平均速度隨時(shí)間的變化，圖6c為地表處天然氣質(zhì)量分?jǐn)?shù)隨時(shí)間的變化。泄漏30 s之前同一部分面域內(nèi)同樣位置泄漏速度散點(diǎn)圖，由圖6a可以看出泄漏10 s時(shí)兩種埋深工況下泄漏速度都是中間大兩邊小，同樣的位置埋深1.4 m都比埋深2.0 m速度大，圖6b、圖6c中為了看的更清晰減少了位置節(jié)點(diǎn)，可以看出這兩個(gè)時(shí)刻泄漏速度中間高兩邊低的特點(diǎn)不是很明顯但同樣的位置埋深1.4 m還是都比埋深2.0 m速度大。分析其原因可能為：高速氣流從泄漏口噴出，大量高濃度天然氣迅速充滿泄漏口附近地帶，由于埋深的存在，天然氣速度會(huì)受土壤阻礙作用速度會(huì)減小，由此可知經(jīng)過越大范圍的土壤層其速度會(huì)減小的越大?？梢缘贸雎裆钆c泄漏速度成反比，埋深越小泄漏地表速度越大。

圖6　不同埋深地表以下1.0 m深處平面內(nèi)天然氣泄漏后的速度、質(zhì)量分?jǐn)?shù)變化圖Fig.6　The changes of leakage natural gas velocity and mass fraction below the surface of the deep 1.0 m under different burial depth conditions

圖7～圖8為不同埋深整個(gè)模擬范圍內(nèi)空氣域中天然氣的泄漏擴(kuò)散范圍3D顯示。

可以看出10 s時(shí)天然氣還未擴(kuò)散到模擬范圍邊界，隨著泄露的進(jìn)行擴(kuò)散范圍越來越大，60 s時(shí)已經(jīng)擴(kuò)散到邊界，其原因同圖6中結(jié)論：埋深越深，地表泄漏速度越小。綜合圖7、圖8可以得出埋深與擴(kuò)散范圍成、高濃度范圍反比，埋深越小擴(kuò)散范圍越大、高濃度范圍越大。

圖7　埋深1.4 m，空氣中天然氣泄漏情況Fig.7　The leakage situation of natural gas in the air at 1.4 m burial depth

圖8　埋深2.0 m，空氣中天然氣泄漏情況Fig.8　The leakage situation of natural gas in the air at 2.0 m burial depth

為了更加清楚直觀地分析管道埋深對(duì)天然氣泄漏的影響，對(duì)相關(guān)數(shù)據(jù)進(jìn)行分析處理，做出不同埋深危險(xiǎn)邊界的覆蓋范圍3D圖，如圖9所示。由3D圖可以清楚地看出埋深為1.4 m的工況危險(xiǎn)范圍大于埋深2.0 m的工況。其原因?yàn)椋汗艿缆裼诘叵?，泄漏初期大量氣體首先擴(kuò)散到土壤中，受土壤孔隙阻力、毛細(xì)管壓力的影響，埋于地下的深度越大，氣體湍動(dòng)能的損失越大。

圖9　不同埋深，空氣中天然氣泄漏準(zhǔn)確范圍Fig.9　The accurate leakage range of natural gas in the air at different burial depth conditions

3　結(jié)　語(yǔ)

（1）管道泄漏后相同條件下硫化氫的危險(xiǎn)區(qū)域總是比甲烷的危險(xiǎn)區(qū)域大，建議事故處理時(shí)首先考慮硫化氫的危險(xiǎn)區(qū)域。

（2）忽略埋深將導(dǎo)致硫化氫軸向危險(xiǎn)區(qū)域擴(kuò)大化，地表危險(xiǎn)區(qū)域比實(shí)際情況小，得出的危險(xiǎn)區(qū)域結(jié)果不精確。對(duì)于軸向，埋深對(duì)靜風(fēng)天氣下硫化氫的危險(xiǎn)區(qū)域影響較大，忽略埋深的比有埋深的危險(xiǎn)區(qū)域大；埋深對(duì)有風(fēng)天氣硫化氫的危險(xiǎn)區(qū)域影響不大，但對(duì)擴(kuò)散范圍影響較大，擴(kuò)散范圍忽略埋深的比有埋深的大。對(duì)于地表，埋深對(duì)靜風(fēng)天氣下硫化氫的危險(xiǎn)區(qū)域影響較大，忽略埋深的比有埋深的小但擴(kuò)散范圍忽略埋深的比有埋深的大；埋深對(duì)有風(fēng)天氣硫化氫危險(xiǎn)區(qū)域影響不大但擴(kuò)散范圍有埋深的比忽略埋深的大。

（2）對(duì)埋深1.8 m，泄漏口速度243.42 m/s，風(fēng)速為5 m/s工況進(jìn)行數(shù)值模擬。研究結(jié)果表明，含硫化氫6.5%的高含硫天然氣泄漏690 s時(shí)達(dá)到穩(wěn)定，徑向危險(xiǎn)區(qū)域約為（-170，80），軸向危險(xiǎn)區(qū)域約為（-435，122），高度達(dá)到40 m左右?？蔀槭鹿屎髶岆U(xiǎn)救災(zāi)提供較為準(zhǔn)確的危險(xiǎn)區(qū)域預(yù)測(cè)。

［1］秦政先.天然氣管道泄漏擴(kuò)散及爆炸數(shù)值模擬研究［D］.南充：西南石油學(xué)院，2006：1-5.

［2］畢靜，張默，何利民.天然氣管道泄漏擴(kuò)散規(guī)律CFD研究現(xiàn)狀分析［J］.石油地質(zhì)與工程，2012（6）：1291-131.

［3］趙洪華，艾長(zhǎng)勝.管道傳聲特性及典型破壞聲的識(shí)別［J］.中國(guó)石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2006，30（3）：101-105.

Zhao Honghua，Ai Changsheng.Pipeline acoustic propagation features and recognition of typical damage sounds［J］.Journal of China University of Petroleum：Edition of Natural Science，2006，30（3）：101-105.

［4］李振林，姚孝庭，張永學(xué)，等.基于FLUENT的高含硫天然氣管道泄漏擴(kuò)散模擬［J］.油氣儲(chǔ)運(yùn)，2008，27（5）：38-41.

Li Zhenlin，Yao Tingxiao，Zhang Yongxue，et al.Simulation study on leakage and dispersion of high sulfur nature gas pipeline with FLUENT CFD-based software［J］.Oil& Gas Storage and Transportation，2008，27（5）：38-41.

［5］李朝陽(yáng)，馬貴陽(yáng).埋地與架空輸氣管道泄漏數(shù)值模擬對(duì)比分析［J］.天然氣工業(yè)，2011，31（7）：90-93.

Li Zhaoyang，Ma Guiyang.A comparative study of numerical simulation on underground and aerial pipeline leakage at instantaneous and continuous states［J］.Natural Gas Industry，2011，31（7）：90-93.

［6］錢新明，馮長(zhǎng)根.易燃易爆危險(xiǎn)物質(zhì)泄漏擴(kuò)散仿真及其應(yīng)用的研究［J］.中國(guó)安全科學(xué)學(xué)報(bào)，1997，7（3）：30-33.

Qian Xinming，F(xiàn)eng Changgen.Study on simulationand application of the release and dispersion of flammableand explosivematerials［J］.Journal of China Safety Science，1997，7（3）：30-33.

［7］Sagrado A P G，van Beeck J，Rambaud P，et al.Numerical and experimental modelling of pollutant dispersion in a street canyon［J］.Journal of Wind Engineering and Industrial Aerodynamics，2002，90（4）：321-339.

［8］Berkowicz R.Modelling street canyon pollution：model requirements and expectations［J］.International Journal of Environment and Pollution，1997，8（3）：609-619.

［9］朱紅鈞，林元華，馬成學(xué)，等.平坦地區(qū)含硫化氫集輸管道的泄漏擴(kuò)散模擬［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，31（6）：156-157.

ZhuHongjun，LinYuanhua，MaChengxue，etal.Numerical smulation of sulfured hydrogen-contaning gas pipeline leak and diffusion on flat areas［J］.Joumal of Southwest Petroleum University：Edition of Natural Science，2009，31（6）：156-157.

［10］馬存棟.高壓天然氣管道破裂氣體擴(kuò)散規(guī)律模擬結(jié)果分析［J］.油氣田地面工程，2005，24（3）：8-10.

［11］Liu Y，Zheng J.Numerical simulation of the diffusion of natural gas due to pipeline failure［C］//ASME 2007 pressure vessels and piping conference.American Society of Mechanical Engineers，2007：319-325.

［12］盧濤，姜培學(xué).多孔介質(zhì)融化相變自然對(duì)流數(shù)值模擬［J］.工程熱物理學(xué)報(bào)，2005，26（1）：167-170.

［13］盧濤，佟德斌.飽和含水土壤埋地原油管道冬季停輸溫降［J］.北京化工大學(xué)學(xué)報(bào)，2006，33（4）：37-40.

［14］Beckermann C，Viskanta R.Natural convection solid/liquid phase change in porous media［J］.International Journal of Heat and Mass Transfer，1988，31（1）：35-46.

［15］Wang B X，Du J H.Forced convective heat transfer in a vertical annulus filled with porous media［J］.International Journal of Heat and Mass Transfer，1993，36（17）：4207-4213.

［16］Hamdami N，Monteau J Y，Le Bail A.Simulation of coupled heat and mass transfer during freezing of a porous humid matrix［J］.International journal of refrigeration，2004，27（6）：595-603.

［17］李又綠，姚安林，李永杰.天然氣管道泄漏擴(kuò)散模型研究［J］.天然氣工業(yè)，2004，24（8）：102-104.

Li Youl¨u，Yao Anlin，Li Yongjie.Study on diffusion of gas pipeline leaking［J］.Natural Gas Industry，2004，24（8）：102-104.

馬貴陽(yáng)，1965年生，男，漢族，內(nèi)蒙古寧城人，教授，博士，主要從事多孔介質(zhì)傳熱理論及計(jì)算流體力學(xué)的研究。E-mail：guiyangma1@163.com王新，1988年生，男，滿族，遼寧開原人，碩士，主要從事油氣管道輸送技術(shù)。E-mail：wangxin_0414@163.com

李 丹，1988年生，女，漢族，河北衡水人，助理工程師，主要從事油氣管道設(shè)計(jì)工作。E-mail：danli0117@163.com

劉兵兵，1986年生，男，漢族，遼寧盤錦人，碩士，主要從事石油天然氣管道安全施工技術(shù)工作。E-mail：113992871@qq.com

黃騰龍，1988年生，男，漢族，遼寧臺(tái)安人，碩士，主要從事給排水和消防的研究。E-mail：276671985@qq.com

編輯：王旭東

編輯部網(wǎng)址：http：//zk.swpuxb.com

Influence of Burial Depth for Natural Gas Pipeline Leakage in Flat Country

Ma Guiyang1，Wang Xin1*，Li Dan2，Liu Bingbing3，Huang Tenglong3
1.College of Petroleum Engineering，Liaoning Shihua University，F(xiàn)ushun，Liaoning 113001，China 2.China Petroleum Pipeline Engineering Corporation，Langfang，Hebei 065000，China 3.China Petroleum Pipeline Bureau，Langfang，Hebei 065000，China

Aiming at the leak and diffusion problem of the buried gas pipeline，the three-dimensional numerical simulation is carried out for the natural gas pipeline leakage in flat area with different buried depth.According to the principle of single variable in the same climate conditions，we take the buried depth as a single variable for different conditions，and study the effect of buried depth of 1.4 m and 2.0 m on the leakage underground，the surface，and in the air.The results show that the buried depth has a great influence on leakage.The results from 0 conditions are great error compared with those of buried depth of 1.4 m and 2.0 m conditions.The buried depth is inversely proportional to the diffusion range，leakage rate，mass fraction and high concentration range.The shallower the buried depth，the bigger the diffusion range，leakage rate，mass fraction and high concentration range.

natural gas leak；three dimensional buried pipeline；deathtrap；burial depth；numerical simulation

10.11885/j.issn.1674-5086.2013.03.14.01

1674-5086（2015）02-0165-08

TE99

A

2013-03-14網(wǎng)絡(luò)出版時(shí)間：2015-03-30

王 新，E-mail：wangxin_0414@163.com

遼寧省自然科學(xué)基金（201102118）。