結(jié)合管道過程，就ESD閥關(guān)斷順序和延遲對管道流量和泄漏速率的影響進(jìn)行了分析。結(jié)果表明，隨著泄漏孔徑的增加，管道內(nèi)將出現(xiàn)逆流現(xiàn)象，管道ESD閥關(guān)斷延遲順序?qū)π】仔孤┖痛罂?a target="_blank" href="http://www.13house.cn/leakhunting/">泄漏速率的控制存在明顯不同的影響，對釋放持續(xù)時間的影響較小。

1、天然氣管道泄漏速率

　　天然氣擴(kuò)散是一種非常復(fù)雜的大氣污染物擴(kuò)散過程，在對其模擬評價過程中存在大量不確定性因素，這些不確定性因素直接影響了模擬結(jié)果的正確性，泄漏速率或泄漏量即為其中之一。天然氣泄漏速率是模擬評價天然氣管線泄漏事故過程中的重要因素之一，對其研究將有利于提高模擬評價結(jié)果的正確性。應(yīng)用熱力學(xué)和氣體動力學(xué)理論，結(jié)合理想氣體狀態(tài)方程、絕熱方程和能量守恒方程，研究分析了天然氣管道的泄漏過程，分別給出了在臨界泄漏階段以及亞臨界泄漏階段的泄漏速率計算公式，為天然氣管輸安全工程的相關(guān)計算提供了可靠的理論依據(jù)。

2、基于HYSYS的動態(tài)天然氣集輸管道泄漏模型

　　天然氣管道泄漏模型與適用性分析：

　　1)穩(wěn)態(tài)泄漏模型。這類模型最常見的是一般泄漏速率模型，利用氣體在泄漏點流動狀態(tài)有關(guān)的特性，認(rèn)為泄漏發(fā)生時泄漏處管內(nèi)壓力不變，通常根據(jù)各種守恒方程可得出解析解的泄漏速率方程。該方法適用于泄漏時管道內(nèi)氣體壓力恒定工況。也有很多研究將小孔泄漏近似為穩(wěn)態(tài)泄漏進(jìn)行計算。

　　2)動態(tài)非穩(wěn)定泄漏模型。目前動態(tài)泄漏模型均主要對大孔泄漏和管道泄漏情景進(jìn)行計算和研究，因難以確定實際泄漏量，基本都假設(shè)管道發(fā)生泄漏后ESD閥同時關(guān)閉，認(rèn)為泄漏量即為兩ESD閥間氣體量，并根據(jù)能量守恒方程等得到泄漏速率方程。實際上，天然氣管道泄漏后，由于受泄漏檢測技術(shù)、人員干預(yù)、執(zhí)行機(jī)構(gòu)響應(yīng)等因素的制約，泄漏處上下游ESD閥關(guān)斷動作并不能保證同步，即存在不同程度的關(guān)斷延遲。在此借鑒加拿大阿爾伯特能源與公共事業(yè)局(AEUB)推薦的天然氣管道有效泄漏長度概念，建立管道泄漏的等效開口系統(tǒng)模型，結(jié)合過程模擬技術(shù)，研究ESD閥關(guān)斷延遲對泄漏過程的影響。

3、ESD閥關(guān)斷動作對泄漏率的影響

3.1、計算算例條件

　　由于高壓天然氣管道泄漏速度較高，為清晰對比泄漏后的動態(tài)平衡過程和ESD閥關(guān)斷延遲影響，基于前述理論和模型，建立20000m長、內(nèi)徑為254mm水平天然氣管道HYSYS泄漏等效模型。計算管道泄漏常見的小孔、大孔范圍內(nèi)的泄漏口面積，小孔取內(nèi)管道內(nèi)橫截面積的1%(50615mm²)，大孔取10%(5065mm²)。對于全尺寸斷裂，泄漏視為管道泄漏模式，泄漏速率為上下游管道單獨泄漏速率之和，2個ESD閥延遲相互影響很小，篇幅所限，這里不再分析。為求不同泄漏孔徑的氣體泄漏速度系數(shù)，將等效模型內(nèi)介質(zhì)指定為純甲烷并與美國環(huán)保署和國家海洋大氣管理局推薦的ALOHA計算數(shù)據(jù)進(jìn)行對比驗證，雖然ALOHA只能對純物質(zhì)提供ESD閥關(guān)閉后泄漏量的分鐘內(nèi)均值，但是可以近似得到不同泄漏孔徑的泄漏系數(shù)，提高模型計算泄漏速率準(zhǔn)確度。在此小孔和大孔速度系數(shù)分別取0.65和0.85。設(shè)定的條件為：

　　1)因管道處于集輸管網(wǎng)的背景壓力環(huán)境，指定管段的入口壓力(本文取95.2kPa)，下游出口壓力(穩(wěn)態(tài)模式求解后指定);

　　2)入口氣體溫度為45℃;

　　3)指定實際的氣體成分，相關(guān)參數(shù)如流經(jīng)上下游ESD閥的流量，泄漏點處壓力、溫度等由HYSYS計算。

　　基于前述構(gòu)建的泄漏等效模型，在驗證其有效性后，本文利用HYSYS的過程計算優(yōu)勢，計算管道摩擦壓力損失變化和沿程熱量損失，將管道氣流改為實際天然氣成分(本文為相對分子質(zhì)量為19.41的高含硫天然氣)，分析管道泄漏段上下游ESD閥關(guān)閉時機(jī)對泄漏的影響。為分析和對比ESD閥關(guān)斷順序?qū)π】准按罂仔孤┧俾实挠绊懀�定義以下動作時間點：

　　1)50min為管道由正常平穩(wěn)運行狀態(tài)轉(zhuǎn)為泄漏事故狀態(tài)的開始時刻;

　　2)150min時上游(或下游)ESD閥立即執(zhí)行關(guān)斷;

　　3)250min時剩余的下游(或上游)ESD閥執(zhí)行關(guān)斷。

　　需說明的是，所設(shè)各動作時間點只為清晰對比ESD閥關(guān)斷影響而設(shè)定，并不代表管道ESD閥關(guān)斷動作的實際執(zhí)行時間。

3.2、不同關(guān)斷條件下的管道流量變化

　　泄漏的發(fā)生使管道出入口流量不再相等，為計算有效泄漏質(zhì)量，在分析泄漏事故狀態(tài)時，了解天然氣管道流量變化有助于分析泄漏速率的變化影響關(guān)系。對于小孔泄漏，泄漏發(fā)生后(50～150min)泄漏段管道入口流量增加，出口流量降低，泄漏量由泄漏段的上游管道所提供，即在泄漏段管道出口無回流。

　　當(dāng)上游ESD閥先關(guān)斷后，下游管道出口將立即出現(xiàn)逆流，補充泄漏量。當(dāng)下游ESD閥先行關(guān)閉時，上游管道入口處流量降低，但持續(xù)補充流量明顯大于上游先關(guān)斷段時的逆向補充流量，將造成較多的氣體釋放。當(dāng)泄漏為大孔泄漏時，泄漏量則由ESD閥上游管道和ESD閥出口的下游管道共同提供，上游入口流量顯著增加，下游管道出口出現(xiàn)大的逆向回流，造成更多的氣體釋放。當(dāng)上游ESD閥先行關(guān)閉后，泄漏管段下游出口的逆向回流量大幅增加，補充泄漏量。當(dāng)下游ESD閥先行關(guān)閉后，管內(nèi)上游流量大幅增加，但相比低于上游ESD閥先關(guān)斷后的出口逆向補充流量。不同孔徑和ESD閥不同關(guān)斷順序下管內(nèi)流量存在差異，流量變化存在動態(tài)平衡過程，這些因素都決定實際的氣體釋放質(zhì)量和速率。顯然，如不能保證ESD閥同時關(guān)斷，應(yīng)優(yōu)先關(guān)斷管道入口ESD閥，此時無論大小孔補充流量都相對較低。

3.3、不同關(guān)斷條件下的泄漏速率變化

　　泄漏段管道出入口流量在不同關(guān)斷順序和時差條件下的變化決定著泄漏率的變化，泄漏發(fā)生后(50min時)由于管內(nèi)壓力存在動態(tài)平衡過程，小孔由于泄漏速率低，平衡時間很短，初始時刻流量稍高，十幾秒內(nèi)達(dá)到壓力平衡后的穩(wěn)態(tài)泄漏流量(由600.215kg·min^-1降至585.59kg·min^-1)。而大孔泄漏后初始時刻流量較高，需經(jīng)過10余分鐘系統(tǒng)方達(dá)到壓力平衡后的穩(wěn)態(tài)泄漏，降幅較高(由5565.21kg·min^-1降至4377.33kg·min^-1)。實際過程中當(dāng)上游ESD閥先關(guān)斷時，大孔和小孔泄漏速率都會降低，但是大孔泄漏速率降低程度明顯高于小孔泄漏，泄漏速率在上游ESD閥關(guān)斷達(dá)到穩(wěn)定后泄漏速率降低34.88%小孔泄漏速率僅僅降低3.45%。當(dāng)下游ESD閥先行關(guān)閉時，大孔泄漏速率也會大幅下降27.53%而小孔泄漏速率并未降低反而有所增加。

4、結(jié)論

　　利用本文方法可對不同距離和梯度變化的管道，定義不同的泄漏位置和尺寸，結(jié)合管道過程控制參數(shù)建立泄漏速率計算模型，為天然氣管道泄漏控制提供有益參考。通過本文研究可得到以下結(jié)論。

　　1)隨著泄漏孔徑的增加，管道內(nèi)將出現(xiàn)回流現(xiàn)象，管道ESD閥關(guān)斷延遲順序?qū)π】仔孤┖痛罂仔孤┧俾实目刂菩Ч�存在不同影響。先關(guān)斷上游ESD閥時，無論小孔還是大孔泄漏，都使下游出口ESD閥處出現(xiàn)逆向流，補充泄漏。而下游先關(guān)斷時，小孔泄漏使上游入口處管道流量有所降低，而大孔泄漏卻使入口處流量大幅增加。

　　2)ESD閥的關(guān)斷延遲和關(guān)斷順序，主要表現(xiàn)在對泄漏速率的影響，即影響單位時間內(nèi)釋放到大氣的天然氣量，而相同間隔不同順序的關(guān)斷對釋放持續(xù)時間的影響較小。

　　3)在現(xiàn)有天然氣管道關(guān)斷多采用壓力、流量監(jiān)控等自動關(guān)斷的現(xiàn)狀下，為最大限度地控制事故天然氣泄漏速率和釋放時間，管道的緊急關(guān)斷應(yīng)優(yōu)先采用聯(lián)動的同時關(guān)斷方式，其次為上游先關(guān)斷方式。