輸氣管道工程放空系統(tǒng)的模擬計(jì)算

趙文利

(中國石油工程建設(shè)有限公司華北分公司，河北 任丘 062550)

隨著國民生活質(zhì)量的提高，人們對綠色節(jié)能的天然氣資源需求量日益增加，保障輸氣管道的安全運(yùn)行已成為國民經(jīng)濟(jì)生產(chǎn)運(yùn)行的基石[1-4]。然而，天然氣屬于易燃易爆氣體，且輸送氣體的長輸管道具有高壓力、長距離、大管徑的特點(diǎn)，因此聚焦天然氣管道的運(yùn)行安全是十分必要的[5-9]。放空系統(tǒng)是輸氣管道運(yùn)行安全的重要保障，通過設(shè)置安全閥對壓力容器超壓情況進(jìn)行泄壓保護(hù)和對故障管道做出迅速反應(yīng)的緊急截?cái)嗟却胧楣艿谰S搶修爭取時間，更重要的是在緊急工況下達(dá)到自動放空、快速切斷的作用[10]。該措施能夠在特殊情況下迅速作出反應(yīng)，為預(yù)防輸氣管道出現(xiàn)放空系統(tǒng)安全問題，減輕放空對環(huán)境不利影響做基礎(chǔ)和鋪墊[11]。

1 安全閥的選型與計(jì)算

根據(jù)國家推薦標(biāo)準(zhǔn)GB/T 12241《安全閥 一般要求》中定義安全閥是一種自動控制閥，能夠保護(hù)系統(tǒng)不超過壓力設(shè)定值，一旦壓力超過規(guī)定值，安全閥會自動打開，將系統(tǒng)中規(guī)定范圍內(nèi)的流體排出，當(dāng)壓力恢復(fù)到正常值后，該閥門自動關(guān)閉防止系統(tǒng)超壓造成事故傷害[12]。因此，安全閥在輸氣站場中至關(guān)重要，對靈敏度也有著較高的要求[13-14]。在計(jì)算安全閥口徑前，應(yīng)首先確定該事故屬于哪種工況，如閥門誤關(guān)斷、控制閥故障及外部火災(zāi)等，從而確定在不同工況下的最大泄放量，進(jìn)而計(jì)算安全閥的最小泄放面積，通過核算安全閥的入口壓力降及安全閥背壓可得到安全閥的喉徑和管徑，最終確定安全閥的型號[15-16]。圖1為安全閥的選型與計(jì)算步驟示意圖。

根據(jù)HG/T 20570.2-95《工藝系統(tǒng)工程設(shè)計(jì)技術(shù)規(guī)定》中的相關(guān)規(guī)定，在不同工況下最大泄放量有所不同。本文中只列出與天然氣輸送管道相關(guān)的事故工況，其他類型在此不在贅述[17]。

1.1 不同工況下的最大泄放量

1.1.1 閥門誤關(guān)閉

①入口截?cái)嚅y打開，出口截?cái)嚅y關(guān)閉時，其泄放量為被封閉內(nèi)管道的最大流量。

②當(dāng)安全閥口徑較小(≤DN25)，管道兩側(cè)的截?cái)嚅y關(guān)閉時，安全閥的泄放量為管道內(nèi)液體的膨脹量。 對于大口徑的長距離管道，液體膨脹量按式 (1)計(jì)算：

Q=B·H/(ρ1·Cp)

(1)

式中：Q為體積泄放量，m3/h；B為體積膨脹系數(shù)，1/℃；H為正常工作時的最大傳熱量，kJ/h；ρ1為液相密度，kg/m3；Cp為定壓比熱，kJ/(kg·℃)。

圖1 安全閥的選型與計(jì)算步驟示意圖

1.1.2 控制閥故障

安裝在設(shè)備入口的控制調(diào)節(jié)閥出現(xiàn)故障，閥門位置為全關(guān)狀態(tài)，泄放量為流過該閥門的最大流量；相反，當(dāng)發(fā)生故障時閥門位置處于全開狀態(tài)，則安全閥的泄放量應(yīng)按式(2)計(jì)算，若為輸氣管道，則應(yīng)滿足低壓側(cè)的設(shè)計(jì)壓力小于三分之二的高壓側(cè)的設(shè)計(jì)壓力。

Qm=3171.3(Cv1-Cv2)Ph(Gg/T)1/2

(2)

式中：Qm為質(zhì)量泄放流量，kg/h；Cv1為控制閥的Cv值；Cv2為控制閥最小流量下的Cv值；Ph為高壓側(cè)工作壓力，MPa；Gg為氣相密度，kg/m3；T為泄放溫度，K。

1.2 判斷流動條件

按照上述計(jì)算得出所需的泄放量，根據(jù)國家標(biāo)準(zhǔn)GB/T 36588《過壓保護(hù)安全裝置 通用數(shù)據(jù)》中第5條安全閥性能的相關(guān)規(guī)定從而確定安全閥的泄放面積[18]。當(dāng)氣體介質(zhì)流經(jīng)安全閥孔道的流量隨著下游壓力的減小將不再持續(xù)增加，則稱該流體為臨界流動，其對應(yīng)的流量為臨界流量，相反則稱為亞臨界流動，其不同的流動條件對應(yīng)不同的流道面積、排放系數(shù)計(jì)算。因此首先判斷該流體是臨界流動還是亞臨界流動。 滿足式(3)為臨界流動，滿足式(4)為非臨界流動。

(3)

(4)

式中：Pb為介質(zhì)流經(jīng)安全閥出口處形成的壓力(簡稱背壓力)，MPa；Pc為排放壓力，MPa；k為絕熱系數(shù)，1。

其中根據(jù)現(xiàn)行國家標(biāo)準(zhǔn)GB50251-2015《輸氣管道工程設(shè)計(jì)規(guī)范》的要求，通過第3.4條輸氣管道的安全泄放中安全閥的定壓要求計(jì)算得到安全閥的定壓值P0[19]，在確定流動條件下，根據(jù)計(jì)算得到的定壓值從而計(jì)算出安全閥的流道面積。

當(dāng)在臨界流動條件下的氣體介質(zhì)的排量按式(5)和式(6)計(jì)算：

(5)

(6)

當(dāng)在亞臨界流動條件下的氣體介質(zhì)的排量按式(7)和式(8)計(jì)算：

(7)

(8)

式中：A為安全閥的閥座與閥瓣間最小面積即流道面積，mm2；C為絕熱指數(shù)k的函數(shù)，1；Kb為亞臨界流動下的理論排量修正系數(shù)，1；Kd為排量系數(shù)，1；Kdr為額定排量系數(shù)，1；M為摩爾質(zhì)量，g/mol；ν0為實(shí)際排放壓力和溫度下的比容，m3/kg；T0為排放溫度，K；Z為壓縮因子，1。

1.3 確定最小泄放面積

根據(jù)化工行業(yè)推薦標(biāo)準(zhǔn)HG/T 20570《工藝系統(tǒng)工程設(shè)計(jì)技術(shù)規(guī)定》的第8.0.2條規(guī)定：

1)對于氣體、蒸汽在臨界條件下的最小泄放面積按式(9)計(jì)算：

(9)

2)對于氣體或蒸汽在亞臨界條件下的最小泄放面積按式(10)計(jì)算：

(10)

式中：A為最小泄放面積，mm2；W為質(zhì)量泄放流量，kg/h；X為氣體特性系數(shù)；P為泄放壓力，MPa；Z為氣體壓縮因子；T為泄放溫度，K；M為分子量；kf為亞臨界流動系數(shù)；C0為流量系數(shù)。

1.4 確定安全閥口徑

根據(jù)上述公式得到的最小泄放面積，計(jì)算安全閥喉徑(d)或閥座口徑(D)。

1)對于全啟式安全閥，其計(jì)算公式為(11)：

(11)

2)對于平面密封型微啟式安全閥，其計(jì)算公式為(12)：

a=πDh

12)

3)對于錐面密封型微啟式安全閥，其計(jì)算公式為(13)：

a=πDhsinφ

(13)

式中：d為安全閥喉徑，mm；h為開啟高度，mm；D為安全閥的閥座口徑，mm；φ為密封面的半錐角，度。

通過不同工況計(jì)算得到的泄放量、判斷流體是臨界流動或非臨界流動從而確定最小的泄放面積最終選擇其合適的口徑，當(dāng)滿足入口壓力降及安全閥背壓時，可確定最終的選型。

2 計(jì)劃性放空

線路段的計(jì)劃放空一般是因線路可能發(fā)生的維檢修、改線、換管、動火連頭等需進(jìn)行放空操作[20-21]。在放空前，需關(guān)斷上下游的截?cái)嚅y，緩慢打開放空閥，通過調(diào)節(jié)閥門開度的大小從而改變放空流量，在規(guī)定時間內(nèi)，完成一定氣量的泄放。本文中利用Pipeline Studio(TG NET)軟件進(jìn)行瞬態(tài)模擬計(jì)算，觀察線路的壓力、放空時間以及放空流量的變化。同時根據(jù)石油化工標(biāo)準(zhǔn)規(guī)范SH3009-2013《石油化工可燃性氣體排放系統(tǒng)設(shè)計(jì)規(guī)范》：“排放系統(tǒng)管網(wǎng)的馬赫數(shù)不應(yīng)大于0.7；可能出現(xiàn)凝結(jié)液的可燃性氣體排放管道末端的馬赫數(shù)不宜大于0.5”，因此根據(jù)放空需滿足環(huán)境要求的條件，下面項(xiàng)目中按照滿足馬赫數(shù)不大于0.5考慮[22]。

某兩個站場間間距為 30.0 km，進(jìn)站壓力為 5.0 MPa，輸氣總量為40×104m3/d，設(shè)計(jì)壓力 6.3 MPa，溫度為 0 ℃，總傳熱系數(shù)為 3.25 W/(m2·℃)，放空閥后的管道長度為 0.1 km，放空口徑為DN50，在滿足最大放空流速(馬赫數(shù))要求的前提下，放空總時間不超過 12 h，對其進(jìn)行線路段放空的動態(tài)模擬分析計(jì)算，在TG NET模型(如圖2所示)中輸入表1的線路段參數(shù)，經(jīng)過計(jì)算，得到表2中線路段管道壓力與放空時間的變化關(guān)系。

圖2 線路段放空示意圖

通過管道壓力與放空時間的變化關(guān)系可看出，計(jì)劃在第 1.5 h 開始放空，計(jì)算結(jié)果表明，線路段放空時間小于 5 h，滿足放空時間 12 h 要求，同時通過該方法可核算出放空立管的口徑是正確的。

表1 線路段參數(shù)表

表2 線路段管道壓力與放空時間計(jì)算結(jié)果表

3 站內(nèi)緊急放空

站內(nèi)緊急放空主要是輸氣站內(nèi)發(fā)生緊急事故(如管道泄漏、發(fā)生火災(zāi)等特殊工況)時利用緊急截?cái)嚅y(ESDV)自動執(zhí)行處置操作的過程[23-25]。根據(jù)GB50251《輸氣管道工程設(shè)計(jì)規(guī)范》中的規(guī)定：“當(dāng)站場出現(xiàn)事故需緊急放空時，應(yīng)能滿足在 15 min 內(nèi)將站內(nèi)設(shè)備及管道內(nèi)壓力從最初的壓力降到設(shè)計(jì)壓力的一半”，同時SH3009-2013《石油化工可燃性氣體排放系統(tǒng)設(shè)計(jì)規(guī)范》對環(huán)境的相關(guān)規(guī)定，由HYSYS軟件進(jìn)行瞬態(tài)模擬，實(shí)現(xiàn)對放空系統(tǒng)中放空量的計(jì)算及限流孔板尺寸的確定。

聯(lián)合被放空管道和容器(可簡化為與實(shí)際容積、內(nèi)部表面積、金屬壁重量相等的一個容器)、節(jié)流閥進(jìn)行模擬，根據(jù)減壓速率要求試算節(jié)流閥(在此認(rèn)為限流孔板)的流通面積，若初選的限流孔板流通面積較小導(dǎo)致 15 min 無法將站內(nèi)管道或容器的壓力降至設(shè)計(jì)壓力的50%，需擴(kuò)大限流孔板尺寸直至滿足減壓速率的要求。

某輸氣站場設(shè)計(jì)壓力 6.3 MPa，最高操作壓力為 6.0 MPa，操作溫度為 20 ℃，輸氣總量為40×104m3/d，設(shè)置有1具放空立管。站內(nèi)設(shè)置有ESD自動放空系統(tǒng)，ESD閥門至放空設(shè)施管道長度 100 m，站內(nèi)管存的有效容積為 60 m3。請核算本工程ESD放空管網(wǎng)限流孔板口徑及放空量。

將站場ESD放空簡化為一個等容積的壓力容器(Vessel)，添加緊急泄放閥(VapourValve)，設(shè)定物流組分、壓力、溫度，緊急泄放閥的流道面積(即限流孔板)等參數(shù)，建立如下模型，如圖3所示。

圖3 站場緊急放空模型圖

圖4 壓力隨放空時間變化結(jié)果圖

通過圖4的壓力隨放空時間變化的結(jié)果圖和表3中的泄放數(shù)據(jù)，發(fā)現(xiàn)當(dāng)時間為 900 s (15 min)時，壓力降至 0.52 MPa，滿足在 15 min 內(nèi)將站內(nèi)設(shè)備及管道內(nèi)壓力從最初的壓力降到設(shè)計(jì)壓力的50%的要求，同時經(jīng)過馬赫數(shù)的計(jì)算[26]，馬赫數(shù)小于0.5，符合對噪音控制的要求，因此該限流孔板的尺寸是合適的。

表3 泄放數(shù)據(jù)表

表3(續(xù))

4 結(jié)論

本文系統(tǒng)地闡述了安全閥的選型與計(jì)算過程，同時利用TG NET軟件和HYSYS軟件模擬了線路段計(jì)劃性放空及站場緊急放空的動態(tài)過程，詳細(xì)地介紹了兩種軟件的操作步驟。研究表明，在長距離輸氣管道放空設(shè)計(jì)中，根據(jù)管道設(shè)計(jì)壓力、溫度、管徑等具體數(shù)據(jù)搭建基礎(chǔ)模型，通過計(jì)算的放空氣量、放空時間及泄放后的壓力等關(guān)鍵數(shù)據(jù)，為能解決長輸管道在放空過程遇到的實(shí)際問題做基礎(chǔ)，同時為指導(dǎo)運(yùn)行管理和維搶修工作提供決策依據(jù)。