摘 要

摘　要：對天然氣輸氣管道開裂泄漏，考慮管道有連續(xù)供氣以及破裂口泄漏符合按時間指數(shù)下降規(guī)律的邊界條件，按管道不定常流動推導連續(xù)泄漏模型，得到關于管道出現(xiàn)較大破裂時，管道流量

摘　要：對天然氣輸氣管道開裂泄漏，考慮管道有連續(xù)供氣以及破裂口泄漏符合按時間指數(shù)下降規(guī)律的邊界條件，按管道不定常流動推導連續(xù)泄漏模型，得到關于管道出現(xiàn)較大破裂時，管道流量和壓力沿長度的分布，以及破裂口泄漏流量計算結果。

關鍵詞：輸氣管道　開裂泄漏　不定常流動

Cracking Leakage Model of Natural Gas Transmission Pipeline

Abstract：For cracking leakage of natural gas transmission pipeline，considering the boundary conditions that pipeline supplies gas continuously and broken leakage is in accordance with exponential decline rule according to time，the continuous leakage model is derived from the pipeline unsteady flow．The flow and pressure distribution along the length and the broken leakage flow calculation results are obtained when a large pipe burst occurs．

Key words：gas transmission pipeline；crackingleakage；unsteady flow

 

1　概述

對于一定長度的天然氣管道，若管道出現(xiàn)較大破裂口或完全斷裂時，發(fā)生較大量的天然氣泄漏，即為開裂泄漏，本文簡稱泄漏。管道開裂泄漏屬于不定常連續(xù)泄漏。

關于天然氣管道泄漏，有很多討論事故原因與對策、檢測方法、事故評價等的文獻，關于輸氣管道開裂泄漏計算研究的報道較少。已有對于天然氣管道泄漏的討論，多限于小孔泄漏。而小孔泄漏，可以足夠近似地看為定常流動，完全可以直接套用現(xiàn)有的孔口臨界出流公式計算。對開裂泄漏，顯然小能簡單地運用孔口臨界出流公式，因為隨泄漏進行，破裂口前的壓力等參數(shù)隨時間變化，是不定常流動(傳統(tǒng)稱為不穩(wěn)定流動)。

對于高壓天然氣輸送管道，當出現(xiàn)開裂泄漏造成較大量的天然氣泄漏時，破裂口天然氣流動將會較長時間都處于臨界流動狀態(tài)，即壓力比b≤bc。

 

式中b——壓力比

pa——大氣壓力(絕對)，Pa

pL——破裂口前天然氣壓力(絕對)，Pa

 

式中bc——臨界壓力比

k——天然氣等熵指數(shù)

對于泄漏點上游持續(xù)供氣的管道，當管道發(fā)生泄漏時，泄漏流量會導致管道內壓力顯著降低，同時管道內沿程阻力和局部阻力也會對破裂口前壓力降低產(chǎn)生影響。總的效果是管道泄漏孔口前的壓力會隨著泄漏的持續(xù)而降低。

本文按上游持續(xù)供氣的高壓天然氣管道開裂泄漏造成天然氣不定常流動推導連續(xù)泄漏模型，得到管道流量和壓力沿長度的分布，以及破裂口泄漏流量計算結果。

2　天然氣管道開裂泄漏方程推導

依據(jù)管道天然氣流動方程式，并考慮管道阻力，通過對下游用氣負荷周期性變動的天然氣管道不定常流動的分析可知，對長度超過一定距離的管道，用氣負荷流量的波動傳導到管道進口會有很大衰減^[1]，管道供給的流量基本保持定值。因此，對高壓天然氣管道開裂泄漏計算模型，以相應管道的進口流量為常量、破裂口泄漏流動為臨界流動為邊界條件，構成偏微分方程定解問題，分3個步驟進行推導。

2.1　破裂口天然氣泄漏臨界流量

在破裂口，天然氣大量泄出，致使管道壓力發(fā)生顯著下降。破裂口面積為：

Ac＝zA    　　      　　　(3)

式中Ac——破裂口面積，m²

z——破裂口面積系數(shù)

A——管道截面積，m²

考慮到天然氣管道的原有高壓狀態(tài)，破裂口處的天然氣泄漏流動會較長一段時間保持臨界流動狀態(tài)(b≤bc)，隨后轉到亞l臨界流動狀態(tài)(b＞bc)。

其臨界流量按式(4)計算：

 

式中qL——破裂口處天然氣泄漏質量流量，kg／s

m——考慮破裂口環(huán)境對泄漏流動阻礙的流量系數(shù)

c0——破裂口前的天然氣滯止聲速，m／s

pL——破裂口前的天然氣滯止密度，kg／m³

亞臨界流量按式(5)計算：

 

pL、c0隨發(fā)生泄漏的管道內壓力、溫度的降低而減小，因而臨界流量會隨時間有所減小?？紤]到管道內流速有限，故在以下方程的推導中視破裂口前密度為滯止密度pL。

2.2　管道天然氣泄漏流量衰減指數(shù)

不考慮流動阻力，忽略破裂口處管道內流向下游的流量，對管道起點到破裂口的管段，建立天然氣質量平衡方程，用以估計破裂口天然氣泄漏質量流量衰減指數(shù)bb。

-dm＝(qL－qin)dt 　　　　　　　   (6)

式中m——管段內天然氣質量，kg

qin——管段進口天然氣流入質量流量，kg/s

t——時間，s

以泄漏開始為時間起點，t=0，考慮足夠長的管段(例如管段長度≥200km)，管段進口天然氣流入質量流量qin?；颈３殖Ａ?，管道內天然氣密度采用平均值，則有：

 

式中V——管段幾何容積，m³

pin——管段進口天然氣密度，kg／m³

D——管道內徑，m

L——管段起點到破裂口的管段長度，m

采用綜合參數(shù)Zc，在開始b≤bc時，有：

 

式中Zc——綜合參數(shù)，m³／s

在隨后b＞bc時，有：

 

對亞臨界狀態(tài)b＞bc，式(10)中b不是常量，因而Zc不是常量。對天然氣，bc＝0.544，臨界壓力(絕對)為0.18MPa。近似認為輸氣管道天然氣泄漏始終保持臨界流動狀態(tài)。

將式(9)代入式(4)、將式(10)代入式(5)，都可得到式(11)：

qL＝ZcpL  　　　　　　　　　  (11)

方程(6)的解為：

 

式中pL0——泄漏開始時刻破裂口前的天然氣滯止密度，kg／m³

e——自然對數(shù)的底

b0——管段進口天然氣流入質量流量    qin為零時的破裂口泄漏質量流量衰減指數(shù)

W——參數(shù)，kg／(m³·s)

由式(11)，對泄漏開始時刻，有：

qL0＝ZcpL0 　　　　　　　　　　　(15)

式中qL0——泄漏開始時刻破裂口泄漏質量流量，kg／s

將式(13)～(15)代入式(12)，有：

qL＝qL0e^-b0t+qin(1-e^-b0t) 　　　　   (16)

設想qL可用更簡單的表達式給出：

qL＝qL0e^-b0t 　　　　　　　　　　  (17)

式中bb——破裂口泄漏質量流量衰減指數(shù)

在實際計算時，由式(16)計算得到ti時刻的qLi。qLi與ti用最小二乘法按式(17)求出bb：

 

式中qLi——第i時刻計算的泄漏質量流量，kg／s

ti——第i時刻的時間，s

n——計算時刻數(shù)

設bb初值，由式(16)計算得到qLi，用n組qLi與ti由式(18)求出bb，經(jīng)迭代計算逼近。

2．3　管段壓力及泄漏量

研究管道氣體連續(xù)泄漏的目的在于求得管段壓力及泄漏量用于實際問題計算或用于理論分析。

考慮過余壓力：

p(x,t)＝pr(x,t)-p(x,0)   　　　 (19)

p(x,0)＝¦(x)　　　　　　　　    (20)

式中p(x，t)——管道過余壓力(絕對)，Pa

x——管道軸向坐標，管段起點為原點，m

t——時間，泄漏開始時為起點，s

pr(x，t)——管道實際壓力(絕對)，Pa

pr(x，0)——管道初始壓力(絕對)，Pa

¦(x)——管道初始壓力(絕對)分布，Pa

 

管道始點壓力pr(0，t)記為pr0，按實際管道運行情況，一般pr0＝C，C為常量。

采用長輸管道流動拋物型方程：

 

式中x——拋物型方程系數(shù)，m²／s

c——聲速，m／s

k——線性化系數(shù)，s^-1

l——摩阻系數(shù)

`w——管道內天然氣線性化平均流速，m／s

初始條件：p(x，0)=pr(x，0)- ¦(x)＝0

邊界條件：

 

由于泄漏質量流量衰減指數(shù)bb隱含了待求函數(shù)p(x，t)，因此破裂口邊界條件未知。為解決此問題，采取了式(18)估計bb，使得破裂口邊界條件成為已知函數(shù)。

采用拉普拉斯變換法求解，得到：

 

式中`t——卷積積分的積分自變量，s

對泄漏口處，x=L，代入式(25)，則有：

 

對式(28)作時間積分還可得到累計泄漏質量：

 

式中mb——累計泄漏質量，kg。

2．4　算例

天然氣壓氣站后一段長輸管道，L＝200km處發(fā)生管道破裂泄漏。管道內徑D=0.7m，壓氣站出口壓力為6.4MPa，管道運行平均壓力為4.5MPa，破裂口處初始壓力為5.5MPa，管道年供氣量為25×10⁸m³／a。天然氣參數(shù)：溫度T=293 K，密度(標準狀態(tài))P=0.7471kg／m³，等熵指數(shù)k=1.309，氣體常數(shù)Rg=517.1J／(kg·K)，壓縮因子Z=0.95。摩阻系數(shù)l=0.02，破裂口面積系數(shù)z=0.10，破裂口流量系數(shù)m=0.43，破裂口泄漏初始流量70.09kg／s。

對式(26)～(28)編程計算得到結果，天然氣長輸管道破裂口泄漏流量及壓力見圖1。計算結果bb=0.6843。

 

p(L，t)計算迭代過程表明，bb值經(jīng)過0.1872、0.4366、0.5615、0.6474、0.6972、0.6843等迭代值變化，收斂到bb=0.6843。p(L，t)曲線的形狀改變十分微小。

3　討論

①對計算模型中的一些參數(shù)，如流量系數(shù)m等可通過實驗進行研究。

②由于管道中天然氣的高壓，在泄漏的前一段時間，破裂口前壓力和泄漏質量流量改變都很小。經(jīng)歷一段時間后，兩者都急劇下降。

③bb經(jīng)迭代計算逼近，由初值0.1872收斂到0.6843，隨時間變化的破裂口壓力p(L，t)曲線改變很小，表明p(L，t)對bb十分不敏感。其實質是，破裂口泄漏流量由臨界流量所控制，是實際的邊界條件。給定的指數(shù)函數(shù)(指數(shù)為bb值)不過是形式上的。但屆。值最后收斂到0.6843，即臨界流量也可由指數(shù)為0.6843的指數(shù)函數(shù)表達。這表明設定破裂口泄漏流量按時間指數(shù)規(guī)律下降是合理的。

④設管道進口流量為0，即是管道放散工況。

4　結論

考慮在發(fā)生開裂泄漏時輸氣管道有繼續(xù)供氣以及破裂口泄漏流量按時間指數(shù)規(guī)律下降的邊界條件，考慮管道流動阻力等因素，采取合理的簡化條件推導的天然氣管道開裂泄漏壓力和泄漏流量方程，可用于對天然氣高壓輸氣管道開裂泄漏事故工況進行預測或實時估算。也可用于管道放散計算。

 

參考文獻：

[1]博布羅夫斯基C A．天然氣管路輸送[M]．陳祖澤，譯．北京：石油工業(yè)出版社，1985：208-209．

 

本文作者：嚴銘卿

作者單位：中國市政工程華北設計研究總院