天然氣混輸管網(wǎng)熱值波動過程模擬

潘國君，郭開華，王冠培，梁金鳳

（中山大學工學院，廣東 廣州510006）

引 言

天然氣作為清潔能源在我國得到迅猛發(fā)展，2012年我國天然氣消費量達1438億立方米，同比增長10.2%［1］。天然氣發(fā)電作為一種新型、清潔、高效的能源在國內(nèi)外受到極大重視，特別是在燃氣發(fā)電當中。在國外天然氣發(fā)電裝機容量已經(jīng)接近或者超過常規(guī)的水、火電裝機容量；在國內(nèi)天然氣發(fā)電雖剛剛起步，但天然氣發(fā)電的優(yōu)點也得到了國家和民眾的認識與重視，特別是廣東地區(qū)，隨著多方位多氣源局面的形成［2］，LNG項目逐步發(fā)展，燃氣發(fā)電用氣量也在逐步增加［3］。

國內(nèi)天然氣氣源種類繁多，多氣源接入同一管網(wǎng)系統(tǒng)的格局逐步形成，對下游用戶的安全、高效運行提出了新的挑戰(zhàn)。相較于民用用戶，工商業(yè)用戶對氣質(zhì)質(zhì)量具有更高的要求。以燃氣電廠為例，一般燃氣電廠常用的燃氣機組有E級和F級，E級機組運行中熱值允許變化速率要求在0.1%·s－1之內(nèi)且熱值波動范圍≤±5%，F(xiàn)級機組運行中熱值允許變化速率要求在4%·min－1之內(nèi)且熱值波動范圍≤±2%。目前研究表明，國內(nèi)現(xiàn)階段的氣源之間均要求滿足熱值范圍的互換要求［4］。但對于兩種滿足互換要求的氣源進行相互替換時造成的熱值波動速率能否達到燃氣機組的運行要求，則需要進一步研究。

本研究是在中山大學－BP液化天然氣教育培訓與研究中心燃氣混輸管網(wǎng)動態(tài)仿真平臺［5］的基礎(chǔ)上建立了天然氣混輸管網(wǎng)氣質(zhì)質(zhì)換情況下的仿真模型，可以計算氣質(zhì)質(zhì)換后下游各位置的熱值變化及熱值變化率情況，為下游工商業(yè)用戶提供參考。

1 混輸管網(wǎng)熱值波動仿真模型

混輸管網(wǎng)熱值波動一般是由不同熱值的氣體質(zhì)換造成的。氣質(zhì)質(zhì)換指的是兩種不同氣質(zhì)的可互換氣源替換時的氣質(zhì)瞬變過程。對于質(zhì)換型混輸管網(wǎng)，熱值的波動與質(zhì)換氣體的熱值變化、管道輸氣流量、管道的幾何尺寸、流速、氣質(zhì)傳遞特性等相關(guān)。對于一個混輸管網(wǎng)，任何一點的熱值波動取決于該點到氣質(zhì)質(zhì)換點路由的各段管道的流動和氣質(zhì)傳遞特性。以圖1所示的支狀管網(wǎng)為例，對于任意節(jié)點k，該點到質(zhì)換點的路由管線由k段管道組成，各管道之間通過節(jié)點進行連接，管網(wǎng)中任意節(jié)點的熱值變化由氣質(zhì)質(zhì)換點的熱值波動及i＝1～k的k段管道的流動性質(zhì)和氣質(zhì)傳遞特性決定。

圖1 支狀混輸管網(wǎng)仿真模型平面拓撲Fig.1 Topology of mixing pipeline network

文獻 ［6］通過CFD數(shù)值模擬和理論分析相結(jié)合的方法給出了氣源質(zhì)換條件下單根管道的氣質(zhì)擴散及熱值波動規(guī)律，其中包括管道內(nèi)不同位置處熱值變化率及熱值變化的數(shù)學表達式。對于定常流動的單根氣質(zhì)質(zhì)換管道而言，管內(nèi)天然氣氣質(zhì)濃度及其燃氣熱值的變化規(guī)律遵循沿軸向一維瞬變熱質(zhì)擴散的一般規(guī)律和基本方程［7］。氣質(zhì)質(zhì)換時，流動介質(zhì)中組分沿管的軸向的瞬變質(zhì)量擴散問題可以等效為一個勻速運動的觀測者在混氣交界面處所觀測到的平均流動為零的軸向瞬變質(zhì)量對流擴散問題［8－9］。管道內(nèi)的氣質(zhì)濃度及燃氣熱值隨時間變化規(guī)律滿足正態(tài)概率積分函數(shù)，即誤差函數(shù)［6］。

對于圖1所示的支狀質(zhì)換型燃氣混輸管網(wǎng)而言，任意節(jié)點k處的熱值變化及其熱值變化速率可分別由式 （1）和式 （2）決定

式中，Hk為節(jié)點k處的熱值；Hs為質(zhì)換氣的熱值，Hr為被質(zhì)換氣的熱值，H0為Hs和Hr的平均值；uk為第k段管的氣體流速；xk為節(jié)點k到氣質(zhì)質(zhì)換點的距離，由k點到氣質(zhì)質(zhì)換點的路由管線決定；t為從氣質(zhì)質(zhì)換開始到當前的時間；t＊為等效預混時間，其大小與氣質(zhì)質(zhì)換點的管路結(jié)構(gòu)和切換時間等因素相關(guān)；Dm，eff，k為節(jié)點k處的平均有效質(zhì)擴散系數(shù)，與節(jié)點k到氣質(zhì)質(zhì)換點的路由管線及各管道流動特性相關(guān)。

式中，Dm，eff，i為節(jié)點前路由第i根管道的有效質(zhì)擴散系數(shù)，與 Reynolds數(shù)相關(guān)［10－12］，文獻 ［6］給出了具體表達式：

式中，ν為氣體運動黏度。

任意時刻氣質(zhì)質(zhì)換交界面的位置x0（t）與到氣質(zhì)質(zhì)換點的管道路由中不同管路的流速和幾何尺寸相關(guān)。氣質(zhì)交界面到達任意節(jié)點j的時間可由式（6）確定

對于節(jié)點k，當j≤k且tj－1＜t≤tj時，氣質(zhì)交界面到氣質(zhì)質(zhì)換點的距離x0（t）可由式 （7）確定

當t＞tk時則

2 實例仿真分析及結(jié)果

某干線管網(wǎng)由20個分輸站組成，全長441km，運行壓力8MPa以上，燃氣經(jīng)過各個分輸站輸送給電廠及商業(yè)用戶，氣源入口為管網(wǎng)唯一供氣氣源，其管網(wǎng)拓撲如圖2所示。氣源供應(yīng)兩種氣質(zhì)的天然氣，分別為熱值Hs＝38.8MJ·m－3（標況狀況，下同）的澳大利亞氣和Hr＝37.6MJ·m－3的卡塔爾氣，在t＝0時刻發(fā)生由卡塔爾氣到澳大利亞氣的質(zhì)換過程。管網(wǎng)用氣負荷一般是在平均負荷附近波動，因此管網(wǎng)給定各節(jié)點日平均小時流量作為邊界條件進行計算，見表1，并取預混時間t＊＝20s。

表1 各節(jié)點外輸流量Table 1 Outflow of nodes

當某節(jié)點熱值與被質(zhì)換氣熱值不同時，認為該節(jié)點熱值開始發(fā)生變化，以每一節(jié)點開始變化的時間為起始點，令這一時刻為t′＝0，得到圖3的各節(jié)點熱值變化曲線?？梢钥闯?，經(jīng)過一段時間之后，各節(jié)點均完成由卡塔爾氣質(zhì)換到澳大利亞氣的過程，在氣質(zhì)質(zhì)換開始和結(jié)束階段變化較緩，中間較陡。各節(jié)點完成質(zhì)換的時間不一樣，熱值變化速率也不一樣。完成質(zhì)換的時間越長，熱值變化越慢，反之越快。

圖2 某市管網(wǎng)拓撲圖Fig.2 Topology of city pipeline network

圖3 各節(jié)點熱值變化曲線Fig.3 Heating value of each node

每個節(jié)點由于與氣源距離和各管道流量的不同，氣質(zhì)變化速率最大的時間也不一樣。表2給出了每個節(jié)點從氣源開始質(zhì)換到熱值變化速率最大時的時間。

表2 各節(jié)點最大熱值變化率時間Table 2 Minutes to reach max heating value rate

圖4給出了各節(jié)點的熱值變化率曲線?？梢钥闯觯S著節(jié)點與氣源點距離的增加，氣質(zhì)變換的時間隨之增加，節(jié)點熱值變化率先逐漸增加后逐漸減小。最大熱值變化率同節(jié)點與氣質(zhì)質(zhì)換點的距離有關(guān)，一般而言距離越遠最大熱值變化率越小。

圖4 各節(jié)點熱值變化率Fig.4 Heating value rate of each node

表3 各節(jié)點最大相對熱值變化率Table 3 Max relative heating value rate of nodes

對于一般燃氣電廠常用的E級和F級機組，允許的熱值 （或沃泊數(shù)）變化范圍分別為±5%和±2%，同時允許的變化速率分別為0.1%·s－1和4%·min－1，由此可得相對熱值允許變化速率分別為0.02s－1和0.033s－1。從圖4可以得出每個節(jié)點的最大相對熱值變化率，見表3。

由表3可以看出，除了福華德電廠外，其他節(jié)點均在燃氣機組熱值變化速率允許范圍之內(nèi)，氣質(zhì)質(zhì)換不會對其產(chǎn)生影響。福華德電廠由于距離氣質(zhì)質(zhì)換點較近，在該處的相對熱值變化速率已經(jīng)大于E級機組所允許的熱值變化速率，因此在氣質(zhì)質(zhì)換時要采取一定的措施。

3 結(jié) 論

混輸管網(wǎng)熱值波動一般是由質(zhì)換兩種不同氣質(zhì)的天然氣引起的。本研究建立了氣質(zhì)質(zhì)換情況下混輸管網(wǎng)熱值波動的仿真模型，模擬了燃氣氣質(zhì)過程中混輸管網(wǎng)不同位置處的熱值及熱值變化率的情況，通過實例工況的模擬驗證此仿真模型能很好地滿足設(shè)計的要求，可以分析質(zhì)換型天然氣混輸管網(wǎng)在不同流量、不同位置處的熱值變化情況，為下游工商業(yè)用戶提供參考，具有實際應(yīng)用價值。

［1］ BP.BP Statistical Review of World Energy［R］.2012.

［2］ Xiang Xiaomin（香小敏），Chen Maohao（陳茂濠），Huang He（黃何），Guo Xianming（郭賢明）.Guangdong province natural gas power generation prospects analysis［J］.Energy Technology（能源技術(shù)），2007，28 （4）：214－217.

［3］ Chen Liang （陳亮），Li Jialong（李嘉龍），Wu Xiaomei（武小梅），Peng Xiangang （彭顯剛）.Problems and solutions of natural gas power generation in Guangdong power grid［J］.ElectricPower（中國電力），2011，（9）：76－79.

［4］ Guo Kaihua（郭開華），Wang Guanpei（王冠培），Huangpu Lixiang（皇甫立霞），Li Ning （李寧）.Standards of natural gas specifications and interchangeability in China ［J］.Natural GasIndustry（天然氣工業(yè)），2011，31 （3）：97－101.

［5］ Wang Guanpei（王冠培），Su Qingbo（蘇清博），Guo Kaihua（郭開華），Li Ning（李寧）.Gas supply program for pipeline network with multiple gas sources［J］.Oil＆GasStorageand Transportation（油氣儲運），2013，32 （10）：63－67.

［6］ Wang Guanpei （王 冠 培 ）. Quality and mixing－flow properties for multi－source pipelines ［D］. Guangzhou：SYSU，2014.

［7］ Eckert E R G，Drake R M Jr.Analysis of Heat and Mass Transfer ［M ］. New York： McGraw－Hill Book Company，1972.

［8］ Freitas Rachid F B，Carneiro de Araujo J H，Baptista R M.Predicting mixing volumes in serial transport in pipelines［J］.JournalofFluidsEngineering，2002，124：528－534.

［9］ Geoffrey Taylor.Dispersion of soluble matter in solvent flowing slowly through a tube［J］.ProcedingoftheRoyal SocietyofLondon（SeriesA：MathematicalandPhysical Sciences），1953，219 （1137）：186－203.

［10］ Geoffrey Taylor.Dispersion of matter in turbulent flow through a pipe ［J］.ProceedingoftheRoyalSocietyof London（SeriesA：MathematicalandPhysicalSciences），1954，223：446－468.

［11］ James B Rawlings，John G Ekerdt.Chemical Reactor Analysis and Design Fundamentals［M］.Madison：Nob Hill Publishing，2006.

［12］ Sittel C N，Threadgill W D，Schnelle K B.Longitudinal dispersion for turbulent flow in pipes［J］.Industrialand EngineeringChemistryFundamentals，1968，7：39－43.