陳 平， 梅華鋒， 趙海林

（北京科技大學機械工程學院，北京 100083）

目前世界上95%以上的天然氣都采用管道輸送方式[1]，但是輸氣管道在運行過程中可能會受到各種因素的干擾，致使管道發(fā)生泄漏事故[2-3]。根據(jù)泄漏孔徑的大小，通常將管道泄漏模型分為小孔模型、管道模型和大孔模型 3種類型[4-5]。國內(nèi)外學者針對上述3種不同失效模式下的氣體泄漏率的計算模型展開了一系列研究，并取得了可喜的成果[3-4,6-12]。但是對于大孔泄漏模型的研究仍然不太成熟，比如不同壓力下的狀態(tài)分析與計算模型等，本文結(jié)合上述文獻，對上面的3種模型的適用范圍進行了研究，對大孔泄漏模型泄漏速率的變化進行了理論分析，并且對大孔泄漏過程中存在的3種狀態(tài)在高壓和低中壓情況下進行詳細的分析討論，得出更加符合實際情況的泄漏計算模型，并且通過實例驗證其正確性，可為長輸管道的定量風險評價提供依據(jù)。

1 基本方程

圖1為管道氣體泄漏示意圖[4]。如圖1所示，管道泄漏點處距離輸氣起始點的間距為Le，1為起始點處中心位置，2為與泄漏孔在同一截面上的管內(nèi)某點，3為泄漏口處位置，Pi、ρi、Ti、ωi（i=1、2、3）分別表示1、2、3點處的氣體壓力、密度、溫度、氣體流速，Pa、ρa、Ta分別表示外界環(huán)境的壓力、密度和溫度。

圖1 管道氣體泄漏示意圖

泄漏孔口處泄漏率的大小取決于泄漏孔處的流動狀態(tài)，其判別標準由臨界壓力比CPR值[13]來確定，即：

式中，Pa：大氣環(huán)境的壓力，Pa（絕）； P2：氣體在點2處的壓力，Pa（絕）；γ：泊松比。

（1）當 Pa/P2＜CPR時，孔口氣體為臨界泄漏，氣體泄漏率m為：

式中，m：氣體泄漏率(kg/s)；Cd：流量系數(shù)，在雷諾數(shù)Re＞30000的非臨界流條件下取0.61，其他情況取 1.0；Ah：泄漏孔口面積(m2)；Wg：氣體的摩爾質(zhì)量(kg·mol-1)；Z：壓縮因子；R：理想氣體常數(shù)，為 8.314(J·mol-1·K-1)；T2：氣體在點2處的溫度(K)。

2 大孔泄漏模型

大孔泄漏模型是一種孔徑介于小孔和管道之間的模型，Montiel等[7]認為在泄漏過程中存在 3種可能的情況：管內(nèi)亞臨界流、泄漏孔口臨界流狀態(tài)；管內(nèi)、泄漏孔口亞臨界流狀態(tài)；管內(nèi)、泄漏孔口臨界流狀態(tài)。董玉華等[3-4]認為第三種情況不存在。而王兆芹等[14-15]則認為在高壓和低壓情況下管內(nèi)、泄漏孔口的狀態(tài)會有所不同。結(jié)合上述文獻，下面分別對 3種情況進行詳細的論述。

（1）管內(nèi)亞臨界流、泄漏孔口臨界流狀態(tài)：當氣體的運動滿足式(4)的條件時，氣體為管內(nèi)亞臨界流、泄漏孔口臨界流狀態(tài)：

氣體泄漏率的計算公式由式(2)來表示。

管內(nèi)點 1和點 2各參數(shù)間的關系式(5)表示[6,16]：

式中，Mai：點 i處的馬赫數(shù)(i=1,2)；ρi：點 i處的氣體密度kg/m3(i=1,2)；Ti：點i處的氣體溫度(K)(i=1,2)。

（2）管內(nèi)、泄漏孔口亞臨界流狀態(tài)：當氣體的運動滿足式(6)的條件時，氣體在管內(nèi)、泄漏孔口均為亞臨界流狀態(tài)。

氣體泄漏率的計算公式由式(3)來表示。

管內(nèi)點1和點2各參數(shù)間的關系由式(5)表示。

（3）管內(nèi)、泄漏孔口臨界流狀態(tài)：當氣體的運動滿足式(7)的條件時，氣體在管內(nèi)和泄漏孔口均為臨界流狀態(tài)。

氣體泄漏率的計算公式由式(8)來表示：

管內(nèi)點1和點2各參數(shù)間的關系由式(5)表示。

3 計算案例與分析

某天然氣長輸管道長度為872 km，管道內(nèi)徑為660 mm，管道內(nèi)壁粗糙度為46 μm，管線的起始點壓力為5 MPa，起始點氣體溫度為293 K。天然氣的摩爾質(zhì)量為 17.1 kg/kmol，壓縮因子為0.9，黏度為 1.01×10-5Pa·s，理想氣體常數(shù)為8.314 J·mol-1·K-1，取泊松比為 1.3，管道在距離起始點等效長度Le=126 km處發(fā)生泄漏。

假設氣體在泄漏過程中，管道起始點的壓力，即點 1處的壓力 P1為定值。根據(jù)文獻中的公式計算小孔泄漏、管道泄漏情況下的泄漏速率以及采用上述討論的大孔泄漏模型得出的計算結(jié)果見圖2～6。

圖2 不同模型的氣體泄漏速率與孔徑比的關系曲線

圖3 小孔模型和大孔模型泄漏速率比值與孔徑比的關系曲線

圖4 大孔模型和管道模型泄漏速率比值與孔徑比的關系曲線

從圖2～4可知，當管道的泄漏孔徑較小，即在其孔徑比（d / D）小于0.15時，小孔模型和大孔模型計算所得的結(jié)果相近，相對誤差小于5%，可以用小孔模型代替大孔模型來計算。隨著泄漏孔徑的不斷增大，在孔徑比（d / D）大于0.9時，大孔泄漏模型所得計算結(jié)果與管道模型相近，相對誤差小于5%，可以用管道模型代替大孔模型來計算。

圖5 不同起始壓力下的氣體泄漏速率與孔徑比的關系曲線（Le=126 km）

從圖5可以看出，隨著起始壓力的增大，同一孔徑比所對應的氣體泄漏速率也逐漸增大，在相同的起始壓力下，氣體的泄漏速率隨著孔徑比的增大而逐漸增大，最終趨于平緩。在P1≤0.184 MPa時，泄漏孔口處的泄漏狀態(tài)一直為亞臨界流狀態(tài)；在0.184 MPa＜P1＜7 MPa時，當孔徑比增大到某一點時，泄漏孔口處的泄漏狀態(tài)將由臨界流狀態(tài)轉(zhuǎn)變?yōu)閬喤R界流狀態(tài)，如P1=5 MPa時，在孔徑比為0.87（泄漏孔徑大約為0.57 m）左右時氣體的泄漏狀態(tài)發(fā)生了變化；在P1≥7 MPa時，當孔徑比從0逐漸增大到1的過程中，氣體的泄漏一直處于臨界流狀態(tài)。從而驗證了上面對大孔模型的論述的正確性，即在高壓時只有管內(nèi)為亞臨界流狀態(tài)、泄漏孔口為臨界流狀態(tài)和管內(nèi)、泄漏孔口均為臨界流狀態(tài)兩種情況，低中壓時只有管內(nèi)、泄漏孔口均為亞臨界流狀態(tài)和管內(nèi)為亞臨界流狀態(tài)、泄漏孔口為臨界流狀態(tài)兩種情況。

圖6 不同泄漏點距離下的氣體泄漏率與孔徑比的關系曲線

從圖6可知，在壓力一定的情況下，泄漏點處離起始點距離越遠，同一孔徑比下的泄漏速率就越小，這主要是管道內(nèi)壁摩擦致使氣體壓力下降所造成的。從圖中還可以看出當P1=5 MPa時，雖然泄漏點距離不同，但是氣體在剛開始時均處于臨界流狀態(tài)，隨著孔徑比的增大，泄漏狀態(tài)都轉(zhuǎn)變?yōu)閬喤R界流狀態(tài)；當P1=7 MPa時，無論泄漏距離為多少，氣體都一直處于臨界流狀態(tài)。這也從另一方面驗證了前面大孔模型討論的正確性。

4 結(jié) 論

（1）當管道的泄漏孔徑較小，即在其相對孔徑比（d / D）小于0.15時，小孔模型和大孔模型計算所得的結(jié)果相近，誤差小于 5%，可用小孔模型來代替大孔模型計算氣體的泄漏速率；當泄漏孔徑接近管徑，即在其相對孔徑比（d / D）大于0.90時，大孔模型計算所得結(jié)果與管道模型趨于相等，誤差小于 5%，可用管道模型來代替大孔模型；在處于中間范圍時，則只能采用大孔模型來計算氣體的泄漏速率。

（2）對于大孔模型，在高壓情況下，只存在管內(nèi)亞臨界流、孔口臨界流和管內(nèi)、孔口均為臨界流兩種情況；在低中壓情況下，只存在管內(nèi)、孔口亞臨界流和管內(nèi)亞臨界流、孔口臨界流兩種情況。

（3）在壓力一定的情況下，泄漏點處離起始點距離越遠，同一孔徑比下的泄漏速率就越小，而且不同孔徑下的泄漏速率變化不同。另外在不同的壓力情況下，氣體泄漏所處的狀態(tài)也會不一樣。

[1] 霍桂霞. 管線球閥的設計特點研究[J]. 寧夏機械,2010,1(2): 19-21.

[2] Fearnehough G A. The control of risk in gas transmission pipeline. Institute of Chemical Engineers,Symposium [S].1985,93: 25-44.

[3] 霍春勇,董玉華,余大濤,高惠臨,周敬恩,馮耀榮.長輸管線氣體泄漏率的計算方法研究[J]. 石油學報,2004,25(1): 101-105.

[4] 董玉華,周敬恩,高惠臨,霍春勇， 馮耀榮. 長輸管道穩(wěn)態(tài)氣體泄漏率的計算[J]. 油氣儲運,2002,21(8): 11-15.

[5] European Gas Pipeline Incident Data Group. Gas Pipeline Incident [R]. European Gas Pipeline Incident Data Group,2005.

[6] Levenspiel O. Engineering flow and heat exchange [M].New York: Plenum Press,1986: 1-10.

[7] Montiel H,Vilchez J A,Casal J,Arnaldos J.Mathematical modeling of accidental gas releases [J].Journal of Hazardous Materials,1998,59(2/3): 211-233.

[8] 肖建蘭,呂保和,王明賢,周加倍. 氣體管道泄漏模型的研究進展[J]. 煤氣與熱力,2006,26(2): 7-9.

[9] 楊 昭,張甫仁,賴建波. 非等溫長輸管線穩(wěn)態(tài)泄漏計算模型[J]. 天津大學學報,2005,38(12):1115-1121.

[10] 賴建波,楊 昭. 高壓燃氣管道破裂的定量風險分析[J]. 天津大學學報,2007,40(5): 589-593.

[11] 吳 起,姜 平,冷世榮,郭培杰. 長輸天然氣管道事故泄漏速率算法改進研究[J]. 中國安全生產(chǎn)科學技術,2012,8(9): 38-42.

[12] 劉中良,羅志云,王皆騰,張 建， 孫曉春. 天然氣管道泄漏速率的確定[J]. 化工學報,2008,59(8):2121-2126.

[13] 蔡鳳英,談宗山,孟 赫,蔡仁良. 化工安全工程[M]. 北京: 科學出版社,2001: 201-204.

[14] 王兆芹. 高壓輸氣管道泄漏模型研究及后果影響區(qū)域分析[D]. 北京: 中國地質(zhì)大學(北京),2009.

[15] 王兆芹,馮文興,李在蓉,李保吉,程五一. 高壓輸氣管道泄漏模型[J]. 油氣儲運,2009,28(12): 28-30.

[16] 王大慶,霍春勇,高惠臨. 長輸管線氣體泄漏率簡化計算方法[J]. 天然氣工業(yè),2008,28(1): 116-118.