燃氣長輸管線動態模擬分析

李曉鳳

（天津海運職業學院，天津 300350）

?

燃氣長輸管線動態模擬分析

李曉鳳

（天津海運職業學院，天津 300350）

摘 要：長輸管線針對燃氣長輸管道流動系統，應用連續性方程、運動方程和氣體狀態方程建立數學模型，聯合初始條件和邊界條件，采用中心隱式差分法進行數值計算，可得到燃氣管道中各點每個時刻的流量和壓力，隨之根據管道流動狀態可以計算管道儲氣量，同時應用模型對文獻中某實例進行計算。結果與文獻數據對比，誤差在5%以內。

關鍵詞：燃氣 長輸管線 動態模擬 中心隱式差法 儲氣能力

引言

能源與環境是困擾人類的兩大問題。煤和石油的燃燒造成環境污染，產生大量的CO2，致使氣候發生了顯著變化。近年來，隨著“霧霾”天氣的頻頻出現，原本最常見的藍天白云成為奢望。天然氣是一種高效、清潔、優質的能源，對改善環境污染起到了積極作用。隨著“西氣東輸”“俄氣中輸”等工程的建成，逐漸形成地區性或全國性的天然氣管網供氣系統。燃氣輸送管道與儲存設備是燃氣進行輸配的載體，為保證合理地向用戶供氣，降低運行費用，需掌握燃氣管道的流動規律。燃氣輸送管道流動動力系統中的很多狀態參數是未知的。掌握這些參數，了解它們之間的規律，對于實際生產、科研都很有裨益［1］。

1 模型的建立

燃氣輸氣管道末端的壓力就是城市儲配站的壓力，流量就是配氣站向城市的供氣量。這個壓力和流量隨著城市的耗氣量在一定范圍內時刻變化著。由于氣體的可壓縮性，壓力、流速和密度這些參數的變化隨管道長度衰減，因而輸氣管工況始終處于動態流動。

輸氣管的直徑在很長距離上是不變的，管路曲率半徑比直徑大很多，垂直于流線的氣體特性變化率同沿著流線的變化率相比可以忽略不計，所以流動可以假定為一維的［2］。對于一維流動，流動運動參數壓力p、密度ρ、速度v，只是時間t和沿管道軸線長度x的函數。

1.1 基本方程

描述氣體管流狀態的參數有三個：壓力p、密度ρ和速度v。為求解這些參數有三個基本方程：連續性方程、運動方程和氣體狀態方程［2-4］。

運動方程中對流向影響小，可以忽略。當標高的差值不太大時，運動方程中的重力項一般也可忽略不計［5］，并與式（1）和（3）聯立簡化為：

1.2 中心隱式差法計算

由式（4）與式（5）組成的方程組為非線性方程組。解析法雖然計算簡單［6］，但是存在較大誤差，主要用于計算機不很發達的年代。本文采用中心隱式差分進行計算。隱式差分法同時考慮所有節點和管道內部網格點，建立統一非線性方程組，運算增加，所需計算機存儲量增加，但是可以保證數值穩定性。

1.2.1 差分方程

圖1 中心隱式差分網格圖

式中：A，b為變量。

將式（6）、（7）代入式（4）和式（5），可得每個管段單元的差分格式。

1.2.2 邊界條件和初始條件

（1）邊界條件

①管道起點（一般為起點）的壓力為定值或為時間函數（已知）；

②已知管道末端流量為時間的函數，可根據城市用氣量的實時參數即可。

（2）初始條件

理論上，初始條件對動態問題的數值計算結果影響不大。但是，對一實際的數值計算格式，若所取初邊值條件不恰當，將影響格式的收斂速度，有時甚至會導致計算發散。在沒有實測數據或預測數據作為初始條件下，按穩態計算結果作為初值。

（3）方程組求解

由前面分析可知，當規定邊界值后，由各網格方程式以及邊界條件和初始條件可組成一封閉的方程組。對整個管道系統的節點和剖面上的變量統一編號，待求變量以向量X表示，其增量用?X表示，各方程用F表示。迭代方程可寫為：

式中：?X為修正向量；F 為函數向量；J（X）為雅可比矩陣。

本文采用MATLAB進行編程，可計算出燃氣管道中各點每個時刻的流量、壓力。

2 結果驗證

本文以文獻［7］中所列數據（某輸氣管道儲氣不穩定流實測數據）為依據進行驗證。應用本文計算方法進行計算，計算結果跟文獻中的末端壓力相近，誤差在±5%以內，可以認為本文的計算是正確的。

計算結果跟原數據有差距，產生這些誤差的主要原因可以歸結為如下幾點：（1）本文在計算時采用了簡化模型，對其影響不大的慣性項、對流項及高差項忽略不計；（2）每個小時內始端的壓力及末端的流量均按一定的數值計算，但是實際情況是這兩個參數時刻在變化；（3）管道中，實際上每個時刻每點內的摩阻系數均不相同，本文計算時采用了一個固定的摩阻系數。（4）管道中燃氣的溫度是變化的，但為了計算方便，本文采用了定溫處理。

3 結論

燃氣輸送管道是燃氣輸配的載體。為了合理供氣、降低運行費用，需掌握管道的流動規律。本文針對管道流動系統，應用連續性方程、運動方程和氣體狀態方程建立模型，進行合理簡化，采用穩定性、準確性較好的中心隱式差分法進行數值計算，可得到管道中任一點的壓力、流量等參數，并可應用模擬結果計算得到管道儲氣量。最后，文章應用該模型以文獻［7］所列數據為依據進行計算，計算結果與文獻數據對比，誤差在5%以內，吻合性較好。

參考文獻

［1］徐彥峰.燃氣管網仿真、優化的研究與開發［M］.哈爾濱：哈爾濱建筑大學，1999.

［2］姚光鎮.輸氣管道設計與管理［M］.青島：石油大學出版社，1991.

［3］江茂澤，徐羽鏜，王壽喜，曾自強.輸配氣管網的模擬與分析［M］.北京：石油工業出版社，1995.

［4］李長俊.天然氣管道輸送［M］.青島：石油工業出版社，2003.

［5］周游，田貫三，張增剛.數值解法模擬長輸管道末端儲氣規律［J］.油氣儲運，2004，（7）：29-34.

［6］孫建國，王壽喜.氣體管網的動態仿真［J］.油氣儲運，2001，（8）：18-21.

［7］李猷嘉.長輸管道末段儲氣量的計算與分析［J］.煤氣與熱力，2002，（1）：8-11.

Dynamic Simulation Analysis on the Long Transmission Gas Line

LI Xiaofeng
（TianjinMarine Career Academy， Tianjin 300350）

Abstract:For the long trans mission gas line system， this paper applied equation of continuity， motion and state of gas to build the mathematical model. Then the flow condition of pipeline can be obtained from using implicit central difference method and the gas storage capacity of line can be received. In the end， this model was used to calculate s ome example and the consequence was compared with the literature data with the error less than 5%.

Key words:long transmission gas line， dynamic simulation，implicit central difference method， gas storage capacity