丁國玉

（北京市城市規(guī)劃設(shè)計研究院，北京　100045）

環(huán)狀燃氣管網(wǎng)非等溫穩(wěn)態(tài)流動計算分析

丁國玉

（北京市城市規(guī)劃設(shè)計研究院，北京100045）

針對當前數(shù)學模型不適用于燃氣環(huán)狀管網(wǎng)的問題，基于燃氣管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型建立了匯流三通非管元件數(shù)學模型，對管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型的方程封閉性進行了分析，給出了算例并進行計算。其計算結(jié)果通過與Pipeline軟件計算結(jié)果對比，最大相對誤差僅為0.29%，驗證了匯流三通數(shù)學模型的適用及可操作性。

環(huán)狀燃氣管網(wǎng)非等溫穩(wěn)態(tài)模型能量方程匯流三通

0　引言

燃氣管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型由于引入了能量方程，計算結(jié)果更加接近實際運行工況。同時，燃氣管網(wǎng)非等穩(wěn)數(shù)學模型又可作為燃氣管網(wǎng)動態(tài)模型離線仿真的初始條件，因此其解法的正確與否直接影響到離線仿真的計算結(jié)果。當前，國外對燃氣管網(wǎng)非等溫穩(wěn)態(tài)模型的研究已經(jīng)成熟，并開發(fā)了相應的軟件［1-3］；國內(nèi)燃氣管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型的研究主要集中在枝狀管網(wǎng)［4-10］，對環(huán)狀管網(wǎng)的研究涉及較少。為此，筆者擬就國內(nèi)環(huán)狀燃氣管網(wǎng)的非等溫穩(wěn)態(tài)數(shù)學模型計算作一分析。

1　非等溫穩(wěn)態(tài)數(shù)學模型的建立

式中，ρ為燃氣密度，kg/m3；p為燃氣絕對壓力，Pa；A為管段的流通橫截面積，m2；M為燃氣質(zhì)量流量，kg/s；x為管段軸向長度，m；λ為摩阻系數(shù)；D為管道內(nèi)徑，m；g為重力加速度，m/s2；θ為燃氣管道對水平面的傾斜角，用弧度表示；h為燃氣比焓，J/kg；K為傳熱系數(shù)，W/（m2·K）；T為燃氣溫度，K；T0為土壤溫度，K。

非管元件模型主要為非管元件進出口間的壓力p、流量M和溫度T的關(guān)系。在城市燃氣管網(wǎng)中常見的非管元件主要為三通［11］。文獻給出的三通數(shù)學模型為：

根據(jù)流向的不同，三通分為分流三通和匯流三通，見圖1和圖2。燃氣流經(jīng)匯流三通時，由于氣源溫度或管道對外傳熱的不同，匯流三通進口處燃氣溫度也不相同，即T1i≠T2k，此時上述模型中的溫度方程不適用于匯流三通。因此，筆者對匯流三通的數(shù)學模型進行推導。

燃氣流經(jīng)匯流三通時為多股燃氣混合過程，忽略局部阻力及位能、動能的變化，則進出口間的流量、壓力和溫度所滿足的方程式如下：

圖1　分流三通示意圖

圖2　匯流三通示意圖

式中，M為質(zhì)量流量，kg/s；p為絕對壓力，Pa；T為溫度，K；h為比焓，J/kg；i，j，k為管段剖分面；1，2，3為管段編號。

2　數(shù)學模型求解

2.1方程封閉性分析

下面結(jié)合一典型環(huán)狀燃氣管網(wǎng)，對非等溫穩(wěn)態(tài)數(shù)學模型封閉性進行分析，見圖3。在圖中，管網(wǎng)由6個節(jié)點、6個管段和6個非管元件組成。節(jié)點Ⅰ、Ⅱ為氣源點，節(jié)點Ⅲ～Ⅵ為分氣點，非管元件1～3、6為分流三通，4、5為匯流三通。

圖3　燃氣環(huán)狀管網(wǎng)示意圖

表1　三通數(shù)學模型所列方程個數(shù)表

2.2計算結(jié)果校驗

在燃氣環(huán)狀管網(wǎng)中，一般給定氣源點的壓力、溫度和各分氣點的流量。此時，非等溫穩(wěn)態(tài)數(shù)學模型的求解為多點邊界問題，可采用Newton-Rapshan迭代法進行求解。

假定天然氣組分為甲烷，氣源點Ⅰ的壓力為4.0 MPa、溫度為293.15 K，氣源點Ⅱ的壓力為4.0 MPa、溫度為288.15 K；各分氣點Ⅲ～Ⅵ對應流量分別為105 kg/s、210 kg/s、140 kg/s、70 kg/s。各管段長度為20 km，管段內(nèi)徑為685.6 mm；管材為鋼管，當量粗糙度為0.04 mm；管段埋深處的溫度為278.15 K，傳熱系數(shù)取2.88 W/（m2·K）。空間計算網(wǎng)絡(luò)步長取1 000 m。經(jīng)求解得出管段各剖分面的壓力及溫度。由于條件限制無法對上述模型求解結(jié)果進行實例驗證，為使計算具有可比性，筆者同時采用Pipeline軟件在相同的條件對上述管網(wǎng)進行了核算。該軟件是一套成熟的商業(yè)軟件，具有一定的參考價值。計算結(jié)果對比見表2。

經(jīng)對比可知，數(shù)學模型計算結(jié)果與Pipeline軟件的計算結(jié)果基本一致。其中，氣源點流量的最大偏差為0.002 kg/s，最大相對誤差為0.001%；分氣點壓力的最大偏差為0.001 MPa，最大相對誤差為0.04%，溫度的最大偏差為-0.808 K，最大相對誤差為0.29%。因此從總體上看，建立的數(shù)學模型計算方法是可以接受的，能滿足一般的工程需求。

3　結(jié)束語

基于高壓燃氣管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型提出了匯流三通的數(shù)學模型，解決了當前數(shù)學模型不適用于燃氣環(huán)狀管網(wǎng)的問題，對環(huán)狀燃氣管網(wǎng)非等溫穩(wěn)態(tài)數(shù)學模型的方程封閉性進行了分析，并給出了算例。通過與Pipeline軟件計算結(jié)果的對比，發(fā)現(xiàn)最大相對誤差僅為0.29%，驗證了匯流三通數(shù)學模型的適用及可操作性。

表2　氣源點流量及分氣點壓力、溫度計算結(jié)果對比表

［1］W.O.Tao，H.C.Ti.Transient analysis of gas pipeline network［J］.Chemica1 engineeringJournal，1998，69（3）：47-52.

［2］Guang-yan Zhun，Michael A.Hensonn，Lawrence Megan.Dynamic modeling and linear model predictive control of gas pipeline network［J］.Journal of process control，2001，11（2）：129-148.

［3］常大海，王善坷，肖尉.國外管道仿真技術(shù)發(fā)展狀況［J］.油氣儲運，1997，16（10）：9-13.

［4］嚴明卿.燃氣輸配工程分析［M］.北京：石油工業(yè)出版社，2007.

［5］李長俊.天然氣管道輸送［M］.北京：石油工業(yè)出版社，2006.

［6］高蘭周.燃氣管道非穩(wěn)態(tài)工況的分析方法與軟件研究［D］.濟南：山東建筑大學，2008.

［7］梁炎明，鄭崗，劉涵，等.天然氣長輸管網(wǎng)穩(wěn)態(tài)計算方法研究［J］.天然氣技術(shù)，2008，2（3）：43-45.

［8］丁國玉.變邊界條件下高壓燃氣管網(wǎng)動態(tài)仿真［D］.濟南：山東建筑大學，2009.

［9］沈孝風.輸氣干線管網(wǎng)瞬態(tài)模擬仿真與優(yōu)化技術(shù)研究［D］.北京：中國石油大學，2010.

［10］劉慶堂.高壓燃氣管網(wǎng)動態(tài)模擬與壅塞流動特性分析［D］.濟南：山東大學，2010.

［11］李興泉，鞠秀峰，田貫三.燃氣管道局部阻力的當量長度法計算［J］.煤氣與熱力，2009，29（6）：9-13.

（編輯：蔣龍）

B

2095-1132（2016）04-0050-03

10.3969/j.issn.2095-1132.2016.04.014

修訂回稿日期：2016-07-12

丁國玉（1983-），碩士，工程師，從事燃氣輸配及利用研究工作。E-mail：dinggy_gh@126.com。