摘　要：輸氣管道氣相泄漏速率計(jì)算是泄漏風(fēng)險(xiǎn)評(píng)估的前提和基礎(chǔ)。為此，通過(guò)搭建氣相管道小孔泄漏實(shí)驗(yàn)系統(tǒng)，開(kāi)展圓孔、周向矩形方孔、軸向矩形方孔的泄漏速率實(shí)驗(yàn)，獲取了不同裂口幾何形態(tài)的氣相泄漏速率。在此基礎(chǔ)上，建立了基于FLUENT的管道孔口泄漏CFD仿真模型，并用其研究了裂口幾何形態(tài)對(duì)中低壓管道小孔泄漏速率的影響機(jī)理、泄漏孔口附近的氣體動(dòng)力學(xué)特征量，包括速度分布、馬赫數(shù)分布等。仿真結(jié)果表明：最大速度發(fā)生在泄漏孔口截面中心處，矩形方孔的最大速度明顯高于圓孔，而裂口方向?qū)ζ溆绊懖伙@著；從臨界壓力比來(lái)看，孔口面積一定，矩形方孔更容易在孔口處達(dá)到臨界流，圓孔、周向矩形方孔、軸向矩形方孔的臨界壓力比模擬值均低于理論計(jì)算值。該實(shí)驗(yàn)結(jié)果有助于氣相管道小孔泄漏研究的深化，也為氣相管道泄漏事故的應(yīng)急處置提供了參考依據(jù)。?

關(guān)鍵詞：輸氣管道??小孔泄漏??實(shí)驗(yàn)??數(shù)值模擬??幾何形態(tài)??泄漏速率??速度分布??馬赫數(shù)分布

Abstract：Gas leakage rate calculation is the premise of and foundation for leakage risk assessment in natural gas pipelines．An experimental system was established for the small leakage of natural gas pipelines to evaluate the respective leakage rates respectively of?round holes，circumferential rectangular square holes and axial rectangular square holes，and to obtain the gas leaking rate of holes with different geometrical characteristics．On this basis，F(xiàn)LUENT-based pipe hole leak CFD simulation was performed to investigate?the influence mechanism of hole geometrical characteristics on medium-and low-pressure pipe hole-leaking rate and gas dynamic characteristic quantities，including speed distribution and Mach number distribution．Simulation results showed that the maximum rate was at the cross-section center of a leaking hole，and that the maximum speed of a rectangular hole was significantly higher thanthat of a square hole，but the influence of hole direction was insignificant；and in terms of critical pressure ratio，if the hole area was determined．it is easier for a rectangular square hole to achieve the critical flow；but as for the critical pressure of a round hole，a circumferential rectangular square hole，and an axial rectangular square hole，their analogous values were all lower than theoretical values．The experimental results are of great help to the deepening of gas pipe hole leaking research and provide scientific guidance for?the emergency response to gas pipe leakage．

Keywords：gas pipeline，hole-leaking，experiment，numerical simulation，geometrical characteristics，leaking rate，speed distribution，Mach number distribution

隨著天然氣管網(wǎng)規(guī)模的不斷擴(kuò)大，由于各種原因?qū)е碌男孤┦鹿蕰r(shí)有發(fā)生。泄漏速率計(jì)算是天然氣管網(wǎng)安全領(lǐng)域的重要問(wèn)題，也是泄漏后擴(kuò)散、風(fēng)險(xiǎn)評(píng)估的前提和基礎(chǔ)。

較早的泄漏速率模型研究主要基于理論推導(dǎo)和數(shù)值計(jì)算，假設(shè)氣體在管內(nèi)為絕熱流動(dòng)，在泄漏點(diǎn)為等熵流動(dòng)，遵循理想氣體狀態(tài)方程和泊松方程，根據(jù)能量守恒和動(dòng)量守恒方程定量描述管內(nèi)氣體的流動(dòng)過(guò)程。Montiel^[1]提出了適用于不同泄漏孔徑的管孔綜合模型，討論了管內(nèi)為亞臨界流、孔口為臨界流或亞臨界流以及管內(nèi)和孔口均為臨界流的泄漏速率計(jì)算方法。國(guó)內(nèi)學(xué)者^[2-4]。通過(guò)對(duì)Montiel模型的研究發(fā)現(xiàn)，管內(nèi)和孔口均出現(xiàn)臨界流的情形在工程實(shí)際中很少見(jiàn)。因此在氣體泄漏速率計(jì)算過(guò)程中僅保留前兩種情形。這類(lèi)模型本身存在一定的缺陷，如沒(méi)有考慮裂口幾何形態(tài)對(duì)泄漏速率的影響，對(duì)于理論推導(dǎo)的泄漏速率缺少相應(yīng)的驗(yàn)證等。

實(shí)驗(yàn)和計(jì)算流體力學(xué)(CFD)可以更詳盡地研究管道泄漏的流場(chǎng)特征。趙金輝等吲提出了大孔亞臨界流泄漏速率的計(jì)算公式，并搭建氣體管道泄漏模擬實(shí)驗(yàn)臺(tái)驗(yàn)證公式的正確性和適用性。Kostowski等^[6]利用實(shí)驗(yàn)和CFD方法比較絕熱模型和等溫模型的適用性。目前，基于實(shí)驗(yàn)和CFD方法研究輸氣管道泄漏主要針對(duì)圓孔，考慮到不同的泄漏部位和破壞原因可能造成不同的裂口形狀和方向^[7]。因此無(wú)法全面描述裂口幾何形態(tài)對(duì)泄漏速率的影響。

筆者的主要研究對(duì)象為中低壓輸氣管道因腐蝕等原因引發(fā)的小孔泄漏，這主要是因?yàn)橄鄬?duì)而言，小孔泄漏在實(shí)際情況中更難發(fā)現(xiàn)與定位^[8]。利用相似原理搭建氣相管道小孔泄漏實(shí)驗(yàn)系統(tǒng)，設(shè)計(jì)并加工不同的泄漏模塊以模擬不同的泄漏場(chǎng)景；以FLUENT為平臺(tái)建立仿真模型，并對(duì)模型的可行性展開(kāi)全面的理論分析和實(shí)驗(yàn)驗(yàn)證，在此基礎(chǔ)上研究裂口幾何形態(tài)對(duì)泄漏速率的影響機(jī)理并得到泄漏孔口附近的氣體動(dòng)力學(xué)特征量，包括速度分布、馬赫數(shù)分布等，使研究趨于精細(xì)化，為管道泄漏事故的應(yīng)急救援和事故調(diào)查提供理論依據(jù)。

利用相似原理搭建氣相管道小孔泄漏實(shí)驗(yàn)系統(tǒng)，實(shí)驗(yàn)過(guò)程中以空氣為介質(zhì)，主要是出于安全和成本方面的考慮，盡管天然氣與空氣性質(zhì)不同，但并不會(huì)影響研究裂口幾何形態(tài)對(duì)泄漏速率的影響。

圖1為氣相管道小孔泄漏實(shí)驗(yàn)系統(tǒng)流程圖，系統(tǒng)主要由動(dòng)力單元、泄漏模塊與測(cè)量?jī)x表單元、數(shù)據(jù)采集單元組成。其中：

?

動(dòng)力單元：包括空氣壓縮機(jī)、緩沖罐、減壓閥，為系統(tǒng)提供壓力恒定的氣源。

泄漏模塊與測(cè)量?jī)x表單元：包括一套可拆卸更換的泄漏模塊以及安裝在泄漏模塊兩端的壓力表、渦街流量計(jì)，模擬泄漏過(guò)程發(fā)生及監(jiān)測(cè)過(guò)程參數(shù)的變化。

數(shù)據(jù)采集單元：包括PLC控制器、計(jì)算機(jī)，實(shí)現(xiàn)過(guò)程參數(shù)的實(shí)時(shí)采集與處理。

實(shí)驗(yàn)過(guò)程中，空氣通過(guò)壓縮機(jī)增壓至緩沖罐，經(jīng)減壓閥減壓后，在泄漏模塊處發(fā)生泄漏，未泄漏的空氣對(duì)大氣進(jìn)行放空；通過(guò)安裝在泄漏模塊兩端的壓力表和渦街流量計(jì)監(jiān)測(cè)壓力、流量參數(shù)的實(shí)時(shí)變化；通過(guò)泄漏模塊兩端渦街流量計(jì)的差值換算得到泄漏速率^[9]。

實(shí)驗(yàn)通過(guò)調(diào)節(jié)減壓閥、更換泄漏模塊等手段構(gòu)建不同的泄漏場(chǎng)景，研究管內(nèi)壓力、泄漏孔口形狀、裂口方向?qū)π孤┧俾实挠绊憴C(jī)理。綜合考慮實(shí)驗(yàn)的可操作性及相關(guān)標(biāo)準(zhǔn)規(guī)范，確定具有代表性的孔口尺度與管道壓力范圍。為盡可能降低讀數(shù)的觀測(cè)誤差及舍入誤差，采用多次實(shí)驗(yàn)取平均值的方法。

實(shí)驗(yàn)所用管道為不銹鋼材質(zhì)，內(nèi)徑50mm，壁厚3.5mm，壓力的選取主要依據(jù)《城鎮(zhèn)燃?xì)庠O(shè)計(jì)規(guī)范》(GB 50028—2006)，結(jié)合減壓閥的工作特性，確定本實(shí)驗(yàn)的管道壓力范圍，如表1所示。

?

設(shè)計(jì)并加工一套可拆卸更換的泄漏模塊，即3根并聯(lián)連接的管道，每根管道長(zhǎng)1.2m，在管道頂部開(kāi)圓孔和矩形方孔，圓孔和矩形方孔具有等水力半徑的幾何特征，矩形方孔的長(zhǎng)邊沿周向或軸向分布，通過(guò)閥門(mén)切換管道模擬不同裂口形狀及方向的泄漏過(guò)程。圓孔孔徑為5mm，主要參照《基于風(fēng)險(xiǎn)的檢驗(yàn)》(API?581—2008)關(guān)于小孔泄漏孔徑的推薦范圍。圖2為泄漏模塊流程圖。

?

圖3為圓孔、周向矩形方孔、軸向矩形方孔的泄漏速率實(shí)驗(yàn)值對(duì)比。由圖3可以看出，對(duì)于圓孔、周向矩形方孔、軸向矩形方孔，管內(nèi)壓力增大時(shí)泄漏速率隨之增加，呈近似線性關(guān)系，因?yàn)闅怏w通過(guò)孔流出的過(guò)程可假定為等熵膨脹，孔口面積一定，隨著管內(nèi)壓力增大，更多的壓力能轉(zhuǎn)化為氣體的動(dòng)能，泄漏速率隨之增加。相同的管內(nèi)壓力下，矩形方孔的泄漏速率顯著高于圓孔，從實(shí)驗(yàn)數(shù)據(jù)來(lái)看，裂口方向?qū)π孤┧俾实挠绊懖伙@著，壓力為0.245MPa時(shí)，周向矩形方孔與軸向矩形方孔泄漏速率的差值僅為0.0323×10^-3kg／s，但總體上，相同的管內(nèi)壓力下，周向矩形方孔泄漏速率略高于軸向矩形方孔。

?

流體的流動(dòng)應(yīng)遵循質(zhì)量守恒定律、動(dòng)量守恒定律和能量守恒定律，而控制方程是上述守恒定律的數(shù)學(xué)描述。建立氣體管道泄漏數(shù)值模擬的控制方程組，包括可壓縮理想流體的質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程，其通用形式是^[10]：

?

從左至右，方程中各項(xiàng)依次為瞬態(tài)項(xiàng)、對(duì)流項(xiàng)、擴(kuò)散項(xiàng)、源項(xiàng)。j是通用變量；G是廣義擴(kuò)散系數(shù)；S是廣義源項(xiàng)。對(duì)于特定的方程，j、G、S具有特定的形式。此外，控制方程組還應(yīng)包括附加的湍流輸運(yùn)方程。u是速度失量。

湍流模型采用標(biāo)準(zhǔn)k-e雙方程模型，具有較高的適用性和精度。標(biāo)準(zhǔn)k-e雙方程模型假定湍動(dòng)黏度是各向同性的，基本輸運(yùn)方程為：

?

式中Gk是由平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng)；Gb是由于浮力引起湍動(dòng)能k的產(chǎn)生項(xiàng)；YM代表可壓湍流中脈動(dòng)擴(kuò)張的貢獻(xiàn)；C1e、C2e和C3e為經(jīng)驗(yàn)常數(shù)；ak和ae分別為湍動(dòng)能是和湍流耗散率e對(duì)應(yīng)的Prandtl數(shù)；Sk和Se是用戶定義的源項(xiàng)。

求解控制方程組采用基于壓力的穩(wěn)態(tài)求解器，選擇該方法是因?yàn)楸疚闹饕P(guān)心近似穩(wěn)定狀態(tài)下的泄漏速率和氣體動(dòng)力學(xué)特征量分布，對(duì)于中低壓輸氣管道，選擇基于壓力的求解器可提高計(jì)算效率同時(shí)可以保證較高的精度。