李 棟，郭 曦，王 迪,，馬鵬博，孫偉棟，王明吉

（1.東北石油大學(xué) 土木建筑工程學(xué)院，黑龍江 大慶 163318；2.東北石油大學(xué) 物理與電子工程學(xué)院，黑龍江 大慶 163318；3.中國(guó)石化管道儲(chǔ)運(yùn)有限公司，江蘇 徐州 221008）

引言

中央排水管是大型浮頂罐必需的排水裝置，其有效防止浮頂雨雪堆積，但由于儲(chǔ)液腐蝕、儲(chǔ)罐沉降、受力不均等原因?qū)е缕淦茡p發(fā)生油品微泄漏，其泄漏通常難以察覺，若不及時(shí)發(fā)現(xiàn)則將導(dǎo)致環(huán)境污染，甚至引發(fā)燃爆事故[1-2]。因此，開展中央排水管泄漏油品擴(kuò)散研究對(duì)監(jiān)控系統(tǒng)設(shè)計(jì)和安全預(yù)警措施制定具有重要工程意義。

目前，關(guān)于浮頂罐中央排水管油氣泄漏的相關(guān)研究較少，而對(duì)輸油輸氣管道及儲(chǔ)罐泄漏的研究較多[3-5]。中央排水管泄漏油氣擴(kuò)散不同于現(xiàn)有管道及儲(chǔ)罐泄漏研究，其需要考慮半封閉空間內(nèi)多組分氣體擴(kuò)散過程。謝飛等人采用FLACS軟件對(duì)浮頂罐密封圈油氣泄漏狀態(tài)進(jìn)行分析，建議在容易泄漏處安裝可燃?xì)怏w監(jiān)測(cè)裝置[6]。而中央排水管管徑較小，無(wú)法在泄漏口安裝氣體監(jiān)測(cè)裝置?？烧{(diào)諧二極管激光吸收光譜技術(shù)（tunable diode laser absorption spectroscopy，TDLAS）不同于傳統(tǒng)的氣體檢測(cè)方法，具有檢測(cè)靈敏度高、精度高、響應(yīng)快等優(yōu)點(diǎn)，成為氣體泄漏安全監(jiān)測(cè)領(lǐng)域研究的熱點(diǎn)[7-10]。然而該技術(shù)是否適用于對(duì)浮頂罐中央排水管泄漏油氣濃度的檢測(cè)，目前尚未發(fā)現(xiàn)相關(guān)研究。

基于以上研究成果，本文建立了浮頂罐中央排水管泄漏油氣擴(kuò)散模型，以中央排水管出口監(jiān)測(cè)點(diǎn)的濃度為對(duì)象，模擬分析了泄漏位置對(duì)中央排水管泄漏油氣擴(kuò)散的影響，并通過實(shí)驗(yàn)驗(yàn)證了中央排水管出口油氣泄漏激光檢測(cè)的可行性。

1 模型建立

1.1 中央排水管泄漏問題

圖1所示浮頂罐中央排水管出現(xiàn)破損點(diǎn)后，受壓力及重力影響罐內(nèi)儲(chǔ)液（通常為原油）從破損點(diǎn)流出。油品泄漏后中央排水管內(nèi)會(huì)聚集大量的輕烴組分（簡(jiǎn)稱油氣）。當(dāng)排水閥處于關(guān)閉狀態(tài)時(shí)，浮盤排水入口在中央排水管泄漏狀態(tài)下變成了油氣擴(kuò)散出口，故出口處油氣濃度檢測(cè)結(jié)果可作為判斷中央排水管是否發(fā)生泄漏的重要依據(jù)。

圖1 浮頂罐中央排水管結(jié)構(gòu)示意圖Fig.1 Structure diagram of central drain pipe of floatingroof tank

1.2 數(shù)理模型

如圖2所示為中央排水管泄漏油氣擴(kuò)散物理模型。

圖2 物理模型Fig.2 Physical model

原油揮發(fā)的油氣成分復(fù)雜，根據(jù)文獻(xiàn)[11]，揮發(fā)過程的油氣組分中丙烷、丁烷及戊烷所占體積分?jǐn)?shù)較高，因此本文選用了丙烷、丁烷、戊烷的混合氣作為泄漏油氣。

中央排水管內(nèi)油氣泄漏擴(kuò)散滿足以下方程[12]：

連續(xù)性方程

式中：ρ為油氣密度；u為油氣在中央排水管內(nèi)的流速。

動(dòng)量守恒方程

式中：P為油氣微團(tuán)上的壓力；μ為油氣動(dòng)力黏度；τ為分子粘性作用而產(chǎn)生的作用在微元體表面上的粘性應(yīng)力；ρa(bǔ)為空氣密度。

能量守恒方程

式中：T為油氣溫度；σT、σc分別為油氣表面拉應(yīng)力和表面壓應(yīng)力；ω為各組分的質(zhì)量分?jǐn)?shù)；cpv、cpa、cp分別為油氣定壓比熱容、空氣定壓比熱容、混合氣體定壓比熱容。

標(biāo)準(zhǔn)k-ε方程

式中：k、ε分別為湍動(dòng)能和耗散率；μt為湍流黏度；Gk、Gb分別為由平均速度梯度和浮力影響引起的湍動(dòng)能產(chǎn)生；YM為可壓湍流項(xiàng)中脈動(dòng)擴(kuò)張的貢獻(xiàn)；Sk、Sε為源項(xiàng)；C1ε=1.44、C2ε=1.92、C3ε=0.09。

組分運(yùn)輸方程

1.3 邊界條件及模型求解

泄漏口采用速度入口邊界，泄漏速度為0.5 m/s，泄漏孔徑為5 mm；出口邊界采用壓力出口邊界條件，出口壓力采用標(biāo)準(zhǔn)大氣壓；其余均采用壁面邊界且無(wú)滑移?；趬毫η蠼馄鬟M(jìn)行求解，壓力與速度耦合選用SIMPLE算法，除能量方程殘差收斂標(biāo)準(zhǔn)設(shè)定為10?6，其余均設(shè)為10?3，迭代時(shí)間步長(zhǎng)設(shè)為0.1 s。

1.4 TDLAS技術(shù)濃度測(cè)量原理

根據(jù)比爾-朗伯定律，推導(dǎo)得出直接吸收光譜法待測(cè)油氣濃度表達(dá)式為[13]：

式中：C為待測(cè)油氣體積分?jǐn)?shù)，%；A為油氣光譜吸光度；P為環(huán)境壓強(qiáng)，atm；S（T）為待測(cè)油氣在該波段的吸收譜線線強(qiáng)，僅與環(huán)境溫度有關(guān)，cm?2·atm?1；L為吸光光程，cm。

2 模型驗(yàn)證

為驗(yàn)證本文模型，與文獻(xiàn)[14]天然氣管道泄漏擴(kuò)散模型進(jìn)行了對(duì)比。文獻(xiàn)中二維模型尺寸為200 m×2 m，泄漏口為速度入口，泄漏速度為100 m/s，泄漏孔徑為10 mm。本文建立了與文獻(xiàn)中尺寸相同的模型，其余設(shè)置參考1.3節(jié)模擬條件。

如圖3所示，通過對(duì)比本文和文獻(xiàn)中模擬10 s的結(jié)果，發(fā)現(xiàn)氣體有著相似的形狀和擴(kuò)散規(guī)律，說明本文模型及求解方法的正確性。

圖3 模型驗(yàn)證Fig.3 Model validation

3 結(jié)果與分析

3.1 模擬結(jié)果分析

當(dāng)泄漏口在Leak1、Leak2、Leak3位置時(shí)（見圖2），報(bào)警時(shí)間分別為561.5 s、1 023.5 s、1 869.5 s（報(bào)警時(shí)間取泄漏初始時(shí)刻至監(jiān)測(cè)點(diǎn)a體積分?jǐn)?shù)到達(dá)0.3%的時(shí)間[15]）。

通過模擬得到監(jiān)測(cè)點(diǎn)a、b、c、d、e、f濃度變化曲線（監(jiān)測(cè)點(diǎn)位置見圖2）。觀察圖4的油氣濃度變化曲線，發(fā)現(xiàn)雖然3幅圖的曲線變化規(guī)律明顯不同，但結(jié)合各監(jiān)測(cè)點(diǎn)的位置，依然可以找到相似的擴(kuò)散規(guī)律：當(dāng)監(jiān)測(cè)點(diǎn)位于泄漏口同一側(cè)時(shí)，各監(jiān)測(cè)點(diǎn)的曲線斜率接近，而不同側(cè)的監(jiān)測(cè)點(diǎn)曲線斜率相差較大。因此可以通過分析管內(nèi)各點(diǎn)濃度變化曲線的斜率，大致分析出泄漏點(diǎn)的位置。

圖4 不同泄漏位置各監(jiān)測(cè)點(diǎn)濃度變化曲線Fig.4 Concentration change curves of each monitoring point at different leak locations

3.2 實(shí)驗(yàn)研究

采用10%濃度的CH4代替成分復(fù)雜的油氣進(jìn)行TDLAS系統(tǒng)監(jiān)測(cè)中央排水管出口濃度的可行性研究，實(shí)驗(yàn)檢測(cè)原理如圖5所示。

圖5 激光檢測(cè)實(shí)驗(yàn)系統(tǒng)示意圖Fig.5 Schematic diagram of laser detection experimental system

圖6所示為本實(shí)驗(yàn)系統(tǒng)的實(shí)物圖。本實(shí)驗(yàn)在壓力為1 atm、溫度為20 ℃的室內(nèi)進(jìn)行。具體操作如下：開啟系統(tǒng)所有設(shè)備，校準(zhǔn)激光器和光電探測(cè)器位置，使數(shù)字示波器輸出信號(hào)為鋸齒波的形狀。關(guān)閉控制閥1并打開控制閥2，調(diào)節(jié)甲烷輸配控制器，設(shè)定流量為100 ml/min后打開氣瓶。CH4在一段時(shí)間后擴(kuò)散至模擬管道外，當(dāng)看到數(shù)字示波器吸收峰穩(wěn)定后開始記錄20組波形數(shù)據(jù)。繼續(xù)將流量分別調(diào)為150 ml/min、200 ml/min，各測(cè)20組數(shù)據(jù)。

圖6 實(shí)驗(yàn)系統(tǒng)圖Fig.6 Picture of experimental system

本文采用直接吸收光譜法進(jìn)行濃度反演，將實(shí)驗(yàn)數(shù)據(jù)導(dǎo)入至計(jì)算機(jī)中進(jìn)行篩選得到如圖7所示的CH4分子原始吸收譜線波形。通過對(duì)該原始信號(hào)進(jìn)行S-G卷積平滑降噪處理以及Voigt線型函數(shù)擬合得到CH4光譜吸光度曲線積分面積A，特定環(huán)境溫度下吸收譜線線強(qiáng)S（T）可根據(jù)HITRAN數(shù)據(jù)庫(kù)查出，在環(huán)境壓強(qiáng)P和測(cè)量光程L已知的情況下，根據(jù)（8）式即可計(jì)算得到待測(cè)CH4體積分?jǐn)?shù)。

圖7 原始吸收譜線波形Fig.7 Original absorption spectral-line waveform

圖8所示Voigt線型函數(shù)擬合殘差，可以看出Voigt線型函數(shù)的擬合殘差在±5%以內(nèi)，曲線擬合優(yōu)度為0.997 9。

圖8 Voigt線型函數(shù)擬合殘差Fig.8 Voigt linear function fitting residual

計(jì)算結(jié)果如圖9所示，圖中CH4的濃度出現(xiàn)了較大的波動(dòng)，這是室內(nèi)空氣的擾動(dòng)和氣體擴(kuò)散的不確定性造成的結(jié)果。排水管出口CH4的濃度隨流量的增大而增大，原因是當(dāng)氣體流量較小時(shí)，氣體動(dòng)能較小，運(yùn)移至排水管出口時(shí)立即向周圍空氣中擴(kuò)散，導(dǎo)致濃度較低。當(dāng)流量為200 ml/min時(shí)，個(gè)別實(shí)驗(yàn)組CH4濃度高于10%，可能因?yàn)闈舛仍诜囱莸倪^程中，平滑降噪、Voigt擬合等處理均會(huì)產(chǎn)生誤差，對(duì)反演結(jié)果造成影響。由此實(shí)驗(yàn)說明，TDLAS技術(shù)可以用于油氣擴(kuò)散檢測(cè)實(shí)驗(yàn)，并具有適用于中央排水管泄漏油氣檢測(cè)的前景，但尚需進(jìn)一步完善其具體測(cè)量方法和求解算法。

圖9 不同流量下CH4濃度反演結(jié)果Fig.9 Inversion results of CH4 concentration at different flow rates

4 結(jié)論

本文通過數(shù)值模擬分析了不同泄漏位置對(duì)油氣擴(kuò)散特性的影響，找到了中央排水管內(nèi)各點(diǎn)濃度變化曲線的斜率與泄漏口間的規(guī)律性。開展了TDLAS技術(shù)檢測(cè)原理的激光檢測(cè)實(shí)驗(yàn)，對(duì)排水管出口的CH4濃度進(jìn)行連續(xù)監(jiān)測(cè)，驗(yàn)證了浮頂罐中央排水管泄漏激光檢測(cè)的可行性，為中央排水管監(jiān)控系統(tǒng)設(shè)計(jì)提供了一定的參考價(jià)值。