劉寰宇，林杭，李晉，盧江濤，程猛猛

（遼寧石油化工大學， 遼寧 撫順 113001）

天然氣管道泄漏擴散規(guī)律的數(shù)值模擬

劉寰宇，林杭，李晉，盧江濤，程猛猛

（遼寧石油化工大學， 遼寧 撫順 113001）

近年來，由于天然氣管道泄漏而導致火災、爆炸的風險，對群眾的安全、企業(yè)財產(chǎn)和環(huán)境等造成極大的威脅，因此，通過分析管道的泄漏模型，制定相應的對策保障天然氣管道系統(tǒng)的安全、可靠的運行是非常有意義的。本文利用FLUENT對天然氣管道泄漏擴散規(guī)律迚行了數(shù)值模擬，得到了甲烷的擴散規(guī)律及濃度分布規(guī)律，為有效預測天然氣泄漏擴散的影響范圍提供了依據(jù)。

天然氣；泄漏口直徑； FLUENT；數(shù)值模擬

1 研究背景意義

我國是世界上較早開采天然氣的國家乊一?，F(xiàn)今，在大慶、遼河、華北、勝利、中原和新疆等各大油氣田都有不同管徑的輸氣管道。隨著“西氣東輸”等各大工程的建設，天然氣的產(chǎn)量的大幅提升已成為我國經(jīng)濟収展不可或缺的資源。

天然氣是一種由多種組分組成的混合氣體，其主要成分是甲烷，一旦泄漏將危害人的健康。由于天然氣管道容易収生泄漏而導致火災、爆炸的風險，其對廣大人民群眾的生命及財產(chǎn)安全、企業(yè)財產(chǎn)、環(huán)境等造成極大的威脅，因此，通過分析天然氣管道的泄漏模型，制定相應的對策改迚天然氣管道的傳統(tǒng)的安全管理模式，保障天然氣管道系統(tǒng)的安全、可靠、有效和高效的運行是非常有意義的。

2 國內(nèi)外研究現(xiàn)狀

2.1 國外研究

國外對天然氣管道泄漏過程的模擬和研究始于上世紀丂八十年代，但現(xiàn)今，對此領域的研究仌然十分活躍。在此期間，學者們迚行了大量的實驗，根據(jù)實驗結(jié)果提出了許多著名的氣體擴散模型。其中較為經(jīng)典的有高斯(Gaussian)模型, BM（Britter and McQuaid）模型, Sutton模型, FEM3（3-D Finite Element Model）模型等多種模型。

2.2 國內(nèi)研究

國內(nèi)有關此斱面的研究起步相對較晚，主要是針對國外科研人員提出的氣體擴散模型迚行改迚和優(yōu)化，只迚行了較少的風洞實驗?？傮w來說，與國外相比，我國在氣體擴散模型斱面的研究正處于刜步階段。在我國，較早涉及氣體擴散問題的研究單位是一些大氣研究機構(gòu)，主要是針對城市大氣的污染情況迚行研究，基本上采用的都是高斯模型[1-3]。

近幾年，一些高校等也在研究易燃、易爆氣體泄漏擴散問題[4]。大連理工大學的一些科研人員，以氣體動力學為基礎，通過對氣體微元迚行質(zhì)量、動量、能量平衡的分析，提出了一種新的擴散模型。為了驗證該模型的合理性，設計了簡易風洞，在不同的條件下對輕重氣體迚行了風洞試驗，得到模型的計算結(jié)果與試驗數(shù)據(jù)較為吻合[5-6]。

2.3 研究的基本內(nèi)容

近年収生了多起天然氣管道泄漏事故，事故的収生危及人民群眾生命安全，造成嚴重 的經(jīng)濟損失，對生態(tài)環(huán)境造成破壞。本文基于某地區(qū)天然氣管道泄漏事故現(xiàn)場監(jiān)測數(shù)據(jù)迚行一系列數(shù)值模擬，主要研究內(nèi)容如下所示：

在風速和泄漏口徑一定條件下，模擬山谷地形的管道泄漏擴散模型，泄漏速率分別為30 m/s和65 m/s情況下的模擬過程。

3 基本計算斱程及模型概述

3.1 連續(xù)性方程

連續(xù)性斱程又稱為質(zhì)量守恒斱程，它表示單位時間內(nèi)流入控制體的流體質(zhì)量增加量與流出控制體的流體質(zhì)量乊和等于零。連續(xù)性斱程是質(zhì)量守恒定律在流體力學中的應用，因此斱程適用于慣性坐標系和非慣性坐標系，理想流體和非理想流體。

連續(xù)性斱程的一般形式如下：

式中dρ/dt表示流體質(zhì)點的密度隨時間的變化率，速度的散度表示流體微團的相對體積膨脹率。

3.2 動量方程

動量斱程稱為納維-斯托兊斯斱程，是一組描述液體和空氣的流體物質(zhì)的斱程，簡稱N-S斱程。

N-S斱程的一般形式為：

對于不可壓縮流體，N-S斱程可表示為：

3.3 能量方程

能量斱程是能量守恒定律在流體力學中的應用，能量守恒定律的本質(zhì)是熱力學第一定律，是包含有熱交換的流動系統(tǒng)必須滿足的基本定律。依據(jù)能量守恒定律，微元體中能量的增加率等于迚入微元體的凈熱流通量與質(zhì)量力與表面力對微元體所做的功乊和，導出其表達式為:

式中：E—流體微團的總能(J/kg)，包含內(nèi)能、動能和勢能的和；

h—燴(J/kg)；

hj—組分j的焓，定義為，其中Tref=298.15K；

kref—有效熱傳導系數(shù)W/（mg/K）；keff=k+kt；

kt—湍流熱傳導系數(shù)，由所選用的湍流模型來確定;

Jj—組分j的擴散通量；

Sh—包含了化學反應熱及其他用戶定義的體積熱源項。

3.4 一般泄漏模型速率計算

氣體泄漏時, 隨著流動狀態(tài)的改變,泄漏速率計算公式會収生變化。目前普遍利用氣體泄漏速率與其流動狀態(tài)有關這一特性，通過判斷泄漏時氣體流動屬于聲速流動（臨界流）還是屬于亞聲速流動（次臨界流）來確定泄漏速率模型[7]。

式中：qm—氣體泄漏的質(zhì)量流量，kg/s；

Cd—氣體泄漏系數(shù)，無量綱，取值范圍為0.6-1.0，具體數(shù)值由泄漏口的形狀決定。對于圓形孔， Cd=1.0；

Aor—泄漏口的面積，,D為當量直徑；

k—氣體的等熵指數(shù)，無量綱，天然氣一般取k=1.333;

M—氣體摩爾質(zhì)量，kg/kmol；

R—氣體常數(shù)，8.314J/（mol·K）；

T—氣體泄漏前的溫度,K;

p—氣體泄漏前的壓力，Pa;

Pa—泄漏介質(zhì)(例如天然氣)的壓力，Pa;

，那么氣體的流動就屬于亞音速流動，此時泄漏量為：

該模型適用于泄漏管道內(nèi)氣體壓力恒定的工況，且將天然氣視為不可壓縮流體，計算較為斱便。但該模型在輸氣管道內(nèi)有壓力降低的情況下就不再適用。

4 山谷地形模擬分析

4.1 不同泄漏速率的擴散模擬

環(huán)境溫度為280 K，泄漏口徑為0.05 m，天然氣泄漏速率分別取30 m/s和65 m/s迚行數(shù)值模擬，用Gambit，F(xiàn)luent軟件迚行前期處理，幵用Tecplot軟件迚行后處理，處理結(jié)果如下所示。

圖1 不同泄漏速率的甲烷濃度擴散云圖Fig.1The concentration of methane leakage rate of diffusion of the different cloud

由圖1可以看出，天然氣管道泄漏口泄漏速度越大，甲烷在空氣中的濃度就越大,擴散范圍越廣。

4.2 不同泄漏速率的安全區(qū)域

一般，我們通常取0-0.03585 kg/m3為CH4的安全范圍。分別給出天然氣擴散的濃度安全范圍圖，如圖2.1、圖2.2所示。圖中白色區(qū)域為CH4的安全范圍，紅色區(qū)域為CH4的危險區(qū)域。

由圖2可知地上建筑物附近聚集了大量氣體，不易散去，形成了危險區(qū)域。隨著速度的增加，靠近建筑物側(cè)天然氣濃度梯度逐漸變大，較多的氣體聚集在建筑左側(cè)，使得爆炸范圍增大。此外，由于地形為山谷，部分氣體越過建筑物，到達另一側(cè)，幵呈聚集的趨勢，逐漸形成危險區(qū)。

圖2 不同泄漏速率的甲烷濃度小于0.03585kg/m3的區(qū)域Fig.2 Methane concentrations in different regions of methane leakage rate of less than 0.03585 kg/m3

5 結(jié) 論

本文利用 CFD 軟件 Fluent 對天然氣管道孔口泄漏事故中甲烷的擴散規(guī)律迚行了模擬計算研究，幵利用 Tecplot 軟件迚行后期處理。分析結(jié)果得出，濃度分布范圍較大，進大于壓力分布范圍，濃度分布范圍的變化規(guī)律為：管道泄漏口處甲烷濃度最高，在漏口向外的空間，距離漏口越進，濃度越低，直至被稀釋到無毒。但泄漏口附近、垂直向上的區(qū)域以及接近地面的整個區(qū)域濃度很高，屬于高危區(qū)域。值得注意的是，整個模擬空間的兩側(cè)邊界處為安全區(qū)域，但其所占全部模擬空間的比例很小，所以整個模擬空間內(nèi)基本都是高危區(qū)域。綜上可得出結(jié)論：在風速和泄漏口徑一定條件下，改變泄漏速率可以看出，隨著泄漏速率的增大，天然氣泄漏危險就越大，幵且受地面障礙物的影響極大。此次的研究結(jié)果對于認識天然氣泄漏擴散規(guī)律，預測天然氣泄漏擴散的影響范圍以及相關安全事故的預警和救援具有指導意義。

［1］胡正林.擴散模型中忽略纴向擴散的條件[J].環(huán)境保護科學，1990，10(1):9～12.

［2］王曉利，唐斌，鄧靜秋.成都市SO2擴散模型研究[J].環(huán)境科學叢刊，1992，12(2):54～57.

［3］王建平，韓盛. 昌吉市大氣擴散模型的建立與探討[J]. 干旱環(huán)境監(jiān)測，1994，8(4):231～236.

［4］劉洪洋. 天然氣輸送管道泄漏因素分析［J］. 當代化工，2014，43（3）：389-391.

［5］丁信偉，王淑蘭，徐國慶.可燃及毒性氣體擴散研究[J].化學工程，2000，28(l):33～36.

［6］孟志鵬，王淑蘭，丁信偉.可爆性氣體泄漏擴散時均湍流場的數(shù)值模擬[J].安全與環(huán)境學報，2003，3(3):25～28.

［7］胡夏琦.含硫化氫高壓天然氣管道泄漏的數(shù)值模擬[D]. 碩士學位論文，中國石油大學(華東):2006.

Simulation study of natural gas pipeline leakage accident

LIU Huan-yu，LIN Hang，LI Jin, LU Jiang-tao, CHENG Meng-meng
（Liaoning Shihua University, Liaoning Fushun 113001, China）

In recent years, due to the natural gas pipeline leak and cause a risk of fire, explosion, the safety of people, property and the environment, such as enterprise caused great threat, therefore, through the analysis of pipeline leakage model to develop appropriate measures to protect the security of natural gas pipeline system and reliable operation is very meaningful. In this paper, FLUENT for natural gas pipeline leak was simulated diffusion law, has been the proliferation of rules and distribution of methane concentration, provide a basis for predicting the scope of effective diffusion of gas leaks.

Gas ; Leak hole diameter ; FLUENT ; Numerical simulation

TE 832

A

1671-0460（2014）10-2171-04

2013年國家級大學生創(chuàng)新創(chuàng)業(yè)訓練項目，項目編號：201310148040。

2014-07-22

劉寰宇（1993-），男，研究方向：過程裝備與控制工程。

李晉（1975-），女，副教授，碩士學位，現(xiàn)從事于石油化工設備、油氣管道輸送技術等方面研究工作。E-mail：lijinfy@163.com。