摘 """""要：基于“碳達峰、碳中和”國家戰(zhàn)略目標，分析現(xiàn)有儲罐油氣蒸發(fā)研究方法，歸納現(xiàn)場實測、風洞實驗和數(shù)值模擬等方法的應用范圍，闡述不同模擬方法的研究對象和優(yōu)缺點，提出油品蒸發(fā)損耗的主要影響因素及改進措施，指出了環(huán)境溫度與太陽輻射對其影響的重要性，為優(yōu)化儲油設備及提高降耗效率提供了改進思路與借鑒。

關(guān) "鍵 "詞：儲罐；蒸發(fā)損耗；現(xiàn)場實測；風洞實驗；數(shù)值模擬

中圖分類號：TE85"""""文獻標識志碼：A """"文章編號：1004-0935（20202024）0×12-1887-05

隨著我國油罐群數(shù)量不斷增多，一系列的環(huán)境安全問題也越發(fā)突出?！半p碳”目標提出以來，我國加快了揮發(fā)性有機物的排放綜合整治，力爭到2025年，實現(xiàn)揮發(fā)性有機物排放總量下降10%以上的目標^[1-2]。儲罐內(nèi)油氣因蒸發(fā)量大、濃度高等特點，受到諸多學者關(guān)注。在早期研究中，認為油氣蒸發(fā)擴散過程是氣液相變、油氣傳質(zhì)和油氣排放的綜合效應，歸納了儲罐油品蒸發(fā)擴散的整體變化規(guī)律，但與實際結(jié)果存在一定偏差^[3]。目前，最有效的研究方法主要是現(xiàn)場實測、風洞實驗和數(shù)值模擬^[4]，通過開展不同外界條件下的現(xiàn)場實測、風洞實驗，可以直接掌握油品的蒸發(fā)擴散機理；利用數(shù)值模擬軟件，可實現(xiàn)擴散物質(zhì)的大范圍追蹤，節(jié)省成本和時間，因此數(shù)值模擬方法在儲罐油品蒸發(fā)擴散機理研究中被廣泛采用。

1 "蒸發(fā)損耗研究方法

現(xiàn)有的儲罐油氣蒸發(fā)研究方法包括：現(xiàn)場實測、風洞實驗與數(shù)值模擬等^[5]。

1.1 "現(xiàn)場實測

現(xiàn)場實測可直接觀察到油品蒸發(fā)損耗的真實情況。劉澤陽^[6]基于現(xiàn)場實測修正了API公式，使聯(lián)合站沉降罐的油品蒸發(fā)擴散規(guī)律的結(jié)論相對誤差從80%降低到10%以下。Tamaddoni等^[7]對碼頭油輪的油氣蒸發(fā)情況進行實驗研究，發(fā)現(xiàn)當油輪裝載30%的原油時，排出氣體中的碳氫化合物含量比裝載前增加15%；當裝載60%原油時，含量減少到9.73%。Hata等^[8]分析溫度和氣體空間等參數(shù)對蒸發(fā)速度的影響，得到了油氣蒸發(fā)速度與晝夜溫差和油罐內(nèi)氣體空間成正比。Lu 等^[9]搭建了油氣的工作損耗和靜止損耗的評估實驗平臺，有效分析了拱頂罐排放的廢氣。

現(xiàn)場實測所得結(jié)論既可用于油氣蒸發(fā)機理的研究，也可為其他研究方法的修正和理論模型的建立提供依據(jù)。

1.2""風洞實驗

風洞實驗因具有優(yōu)越的可控性與穩(wěn)定性，在模型評估和數(shù)據(jù)驗證方面具有獨特優(yōu)勢。ZHANG等^[10]借助風洞平臺搭建內(nèi)浮頂罐縮比模型研究了正己烷蒸氣的擴散規(guī)律，提出儲罐內(nèi)的氣體空間越大，儲罐內(nèi)外的氣流交換越弱，隨著浮盤高度增加，艙內(nèi)形成較大的渦流，加劇氣體擾動。許雪等^[11]建立罐區(qū)模型，通過風洞實驗和數(shù)值模擬，發(fā)現(xiàn)當儲罐高度相等時，儲罐間油氣主要積聚在前排儲罐背風側(cè)下方；當前排儲罐高于后排儲罐時，氣流會挾卷罐間油氣至后排儲罐罐頂，使其油氣質(zhì)量分數(shù)較高于同等高度區(qū)域。為對儲罐油品蒸發(fā)擴散進行全面深入的研究，常州大學研究人員建立了直流式風洞平臺，討論了內(nèi)浮頂罐外風場及風壓分布規(guī)律、風速對內(nèi)浮頂罐油氣流場分布及油氣擴散濃度的影響，得到不同風速下，罐內(nèi)油氣分布整體呈現(xiàn)對稱狀態(tài)，且浮盤位置越低、風速越大，蒸發(fā)速率越快^[12]。

風洞實驗研究主要集中于罐內(nèi)油品的風荷載量、風向和儲罐結(jié)構(gòu)的變化等方面，而對于儲罐內(nèi)部氣體空間的蒸發(fā)實質(zhì)分析以及光照輻射對內(nèi)部氣體空間的影響等研究仍需進一步深入。

1.3""數(shù)值模擬

數(shù)值模擬技術(shù)成本低，周期短，效率高，可進行全尺寸模擬。目前，依據(jù)求解尺度的不同，數(shù)值模擬方法主要分為三種：直接數(shù)值模擬（DNS）、雷諾時均方法（RANS）和大渦模擬（LES）。

DNS方法用于求解低雷諾數(shù)的流體模擬問題。早期DNS方法建立的氣液兩相混輸工藝分析模型，對于氣體擴散的模擬應用較少^[13]。隨著計算機性能的不斷提升，DNS方法現(xiàn)已能夠通過開展儲罐的油氣蒸發(fā)過程研究，對罐體內(nèi)部的傳熱、傳質(zhì)和化學反應運動規(guī)律進行總結(jié)^[14]。通過DNS方法來全面認識油品蒸發(fā)和油氣擴散的機理，有助于優(yōu)化罐體設備結(jié)構(gòu)，提高過程效率。

RANS方法計算效率高，解的精度可以滿足工程實際需要。Tan等^[15]將氣體通過簡單障礙物時擴散的CFD預測與風洞測量結(jié)果進行比較，探究了障礙物引起的流場特性與氣體濃度之間的相關(guān)性。Liu等^[16]建立罐內(nèi)油品的傳熱和流動數(shù)學模型，繪制儲罐油溫場的中心溫度隨時間的動態(tài)曲線，分析了儲罐的軸向和徑向溫度。

LES方法能夠解決了DNS求解時間過長問題，又彌補了 RANS 無法計算非穩(wěn)態(tài)、非平衡的大尺度效應不足^[17]。Sun等^[18]利用LES方法分析了外浮頂儲罐周圍的湍流流動和風荷載，得到LES方法在計算儲罐周圍非定常風荷載和渦流結(jié)構(gòu)方面優(yōu)于RANS方法的結(jié)論。針對強風作用下儲罐的繞流特性和密封圈的油氣耗散特點，李文欣^[19]基于儲罐繞流流場旋渦結(jié)構(gòu)與風洞實驗結(jié)果，分析風速和液位對儲罐內(nèi)外壁風荷載特性的影響規(guī)律，提出了基于LES的計算公式。

基于已有研究成果，對三種模擬方法的優(yōu)勢、不足及適用對象等進行總結(jié)，見表1。

2 "儲罐蒸發(fā)損耗影響因素

2.1 "儲罐內(nèi)部影響因素

2.1.1 "油品影響

儲罐內(nèi)油品蒸發(fā)損耗會受到油品自身因素的影響。Sharma等^[20]對臥式儲罐中的15種汽油進行實驗分析，發(fā)現(xiàn)正丁烷和異戊烷的濃度越高，呼吸損耗越大。此外，儲罐收發(fā)油的裝載因素不同，油氣損耗率也不同。Hassanvand^[21]等利用VOF模型進行數(shù)值模擬，分析了裝油速度、油罐初始油氣濃度對油罐損失率的影響。Wang等^[22]通過模擬表明裝油速度增加，裝油蒸發(fā)損耗率變化不大，約為0.022%；初始油氣質(zhì)量分數(shù)增加，蒸發(fā)損耗率上升明顯，見圖1。

針對收發(fā)油過程中造成的油氣損失，可采用氮封系統(tǒng)維持儲罐內(nèi)壓力平衡，有效減少儲存介質(zhì)蒸發(fā)損耗、防止空氣污染，保證儲存油品質(zhì)量和儲罐安全。

2.1.2 "儲罐影響

我國工業(yè)儲罐結(jié)構(gòu)類型眾多，儲罐的不同構(gòu)造對油氣蒸發(fā)影響不一。Huang、Jing、黃維秋、李飛等基于儲罐自身結(jié)構(gòu)對油品蒸發(fā)的影響進行了大量的工作，取得了一定的科研成果^[23-26]，見表2。

由表2可知，通風口位置、輸油口高度和浮盤位置等不同儲罐構(gòu)造是影響油品蒸發(fā)損耗的主要儲罐因素?；诖?，可以對油罐結(jié)構(gòu)進行改造升級，以減少油品蒸發(fā)損耗。王建偉^[27]通過改變擋板的位置，發(fā)現(xiàn)當呼吸閥擋板的吊桿長度為呼吸閥接合口直徑的2倍時，對降低儲罐內(nèi)油品蒸發(fā)損耗的效果最好。

2.2 "儲罐外部影響因素

大氣溫度是影響儲罐蒸發(fā)損耗的因素之一。Wang等^[28]和Li等^[29]基于數(shù)值模擬結(jié)果與現(xiàn)場監(jiān)測數(shù)據(jù)，提出了有效的物理-數(shù)學模型。張佩宇等^[30]建立了拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型，利用自制的儲罐實測罐內(nèi)氣相溫度變化，對拱頂罐一天內(nèi)的蒸發(fā)損耗過程進行了數(shù)值模擬，測得小呼吸損耗實驗值與模擬值最大誤差不超過4.7%。Ginestet等^[31]提出一種評估蒸發(fā)流量的熱分析方法，討論了溫度、壓力和蒸發(fā)量隨時間的變化規(guī)律，并對蒸發(fā)損耗進行了評估。

輻射是影響儲罐蒸發(fā)損耗的另一重要因素。建立拱頂罐非穩(wěn)態(tài)傳熱傳質(zhì)模型，見圖2，太陽輻射熱量通過大氣經(jīng)罐壁傳遞到油品，由于油品與大氣之間存在較大的溫差，使得靠近儲罐邊界部分的油品通過與內(nèi)壁之間的自然對流、罐壁導熱、強制對流及輻射的方式向外界散失熱量。

Huang等^[32]基于CFD數(shù)值模擬方法，分析不同季節(jié)下太陽輻射對拱頂罐呼吸損耗的影響，得到了不同VOCs的日損失率的估算公式。張佩宇^[30]對拱頂罐進行數(shù)值模擬，發(fā)現(xiàn)總傳熱系數(shù)中與液相接觸的罐壁的總傳熱系數(shù)最大，與氣相接觸的罐壁的次之，罐頂?shù)淖钚　un等^[18]對浮頂儲罐溫度場進行數(shù)值模擬，指出罐頂溫降速率隨著太陽輻射強度的降低而增大；罐壁溫降速率隨大氣溫度的降低而增大；罐底近似于絕熱，溫降速率受外界環(huán)境影響較小。

除了以上兩種因素外，風速對儲罐蒸發(fā)的影響也不能忽視。周寧等^[33]考慮到LNG儲罐蒸發(fā)泄漏過程中的特點，結(jié)合相平衡原理和歐拉多相流模型進行數(shù)值模擬，開展了風速、泄漏速率和地面粗糙度對LNG氣云擴散影響的研究。Fang等^[34]通過分析儲罐的風速場和濃度場的疊加效應，發(fā)現(xiàn)在不同的環(huán)境風速下，油罐之間的干擾是不同的，外浮頂儲罐數(shù)量越多，其在順風方向的疊加效應就越明顯。Kountouriotis等^[35]通過分析汽油擴散規(guī)律，發(fā)現(xiàn)風速可以加速蒸汽的擴散，在高溫和低風速的情況下，罐體泄漏處上方存在可燃云，通過控制揮發(fā)性有機物的排放可將可燃云最小化，從而降低爆燃的風險。

大氣溫度，輻射和環(huán)境風速是影響儲罐蒸發(fā)損耗的重要外部因素。針對儲罐外部因素造成的蒸發(fā)損耗，Chen等^[36]研發(fā)了一種應用于儲罐外壁的聚合物被動輻射冷卻涂層（PRCC），該涂層對陽光的反射率達96.2%，可以減少儲罐內(nèi)50%以上的油品損耗。劉風章^[37]進一步研究發(fā)現(xiàn)，涂覆隔熱反射涂層的罐體全天外表面平均溫度可降低21.1%，罐頂外表面平均溫度可降低 37.1%，罐內(nèi)呼吸氣體平均溫度可降低 27.4%。因此，開發(fā)儲罐隔熱涂層，可以降低罐體所受輻射，達到理想的降耗效果。

3 "結(jié)論與建議

通過蒸發(fā)損耗的產(chǎn)生規(guī)律分析與綜合探討，得到以下結(jié)論與建議：

1）環(huán)境溫度和太陽輻射是影響儲罐蒸發(fā)損耗的主控因素。使用儲罐隔熱涂層控制溫度以降低油品的蒸發(fā)損耗是最為簡單有效的方法。當使用隔熱涂層時，罐體外表面平均溫度可降低20%以上，大幅度減少油品的蒸發(fā)損耗。

2）結(jié)合現(xiàn)有儲罐的研究方法，得到各自的應用范圍?，F(xiàn)場實測的數(shù)據(jù)既可用于油氣蒸發(fā)機理的研究，也可為其他實驗方法的修正和理論模型的建立提供依據(jù)；數(shù)值模擬能夠有效分析罐內(nèi)流場特征，節(jié)約分析時間和成本；風洞實驗的研究主要集中于罐內(nèi)油品的風荷載量、風向和儲罐結(jié)構(gòu)的變化等方面，而對于儲罐內(nèi)部氣體空間的蒸發(fā)實質(zhì)分析等研究仍需進一步深入。

3）從綜合效能出發(fā)，應優(yōu)化呼吸閥結(jié)構(gòu)，減少油面上的氣體空間；當儲罐進行收發(fā)油作業(yè)時，可以采用氮封系統(tǒng)維持儲罐內(nèi)壓力平衡，防止油氣外泄，降低呼吸損耗。

參考文獻：

[1] 王陶， 張志智， 孫瀟磊. 低碳背景下碳市場的發(fā)展及其對石油石化行業(yè)的啟示[J]. 遼寧化工， 2020， 49（6）： 647-650.

[2] SIMAYI M， SHI Y Q， XI Z Y， et al. Emission trends of industrial VOCs in China since the clean air action and future reduction perspectives[J]. Science of the Total Environment， 2022， 826： 153994.

[3] CHEN W H， CHEN Z B， YUAN C S， et al. Investigating the differences between receptor and dispersion modeling for concentration prediction and health risk assessment of volatile organic compounds from petrochemical industrial complexes[J]. Journal of Environmental Management， 2016， 166： 440-449.

[4] 石俊杰， 趙東旭， 王德陽， 等. 儲罐油品蒸發(fā)擴散規(guī)律研究[J]. 遼寧化工， 2022， 51（12）： 1725-1729.

[5] CHEN L， WANG X T， ZHAO Y Q， et al. Research on application of oil and gas recovery technology in storage and transportation system， 2022

[6] 劉澤陽. 聯(lián)合站典型源VOCs排放規(guī)律研究及應用[D]. 東營： 中國石油大學（華東）， 2019.

[7] TAMADDONI M， SOTUDEH-GHAREBAGH R， NARIO S， et al. Experimental study of the VOC emitted from crude oil tankers[J]. Process Safety and Environmental Protection， 2014， 92（6）： 929-937.

[8] HATA H， YAMADA H， KOKURYO K， et al. Estimation model for evaporative emissions from gasoline vehicles based on thermodynamics[J]. Science of the Total Environment， 2018， 618： 1685-1691.

[9] LU C， HUANG H， CHANG S T， et al. Emission characteristics of VOCs from three fixed-roof p-xylene liquid storage tanks[J]. Environmental Monitoring and Assessment， 2013， 185（8）： 6819-6830.

[10] ZHANG G， HUANG F Y， HUANG W Q， et al. Analysis of influence of floating-deck height on oil-vapor migration and emission of internal floating-roof tank based on numerical simulation and wind-tunnel experiment[J]. Processes， 2020， 8（9）： 1026.

[11] 許雪， 陳風， 黃維秋， 等. 基于風洞平臺實驗的大型罐區(qū)溢油事故后的油氣擴散模擬[J]. 環(huán)境工程學報， 2021， 15（12）： 3946-3956.

[12] HUANG W Q， HUANG F Y， FANG J， et al. A calculation method for the numerical simulation of oil products evaporation and vapor diffusion in an internal floating-roof tank under the unsteady operating state[J]. Journal of Petroleum Science and Engineering， 2020， 188： 106867.

[13] BRANFORD S， COCEAL O， THOMAS T G， et al. Dispersion of a point-source release of a passive scalar through an urban-like array for different wind directions[J]. Boundary-Layer Meteorology， 2011， 139（3）： 367-394.

[14] Alizadeh MohammadReza，F(xiàn)atemi Mobeen. Pore-doublet computa- tional fluid dynamic simulation of the effects of dynamic contact angle and interfacial tension alterations on the displacement mechanisms of oil by low salinity water[J].International Journal of Multiphase Flow， 2003， 29（6）： 1039.

[15] TAN W， WANG K， LI C J， et al. Experimental and numerical study on the dispersion of heavy gases in urban environments[J]. Process Safety and Environmental Protection， 2018， 116： 640-653.

[16] LIU Y， FAN J W， CHENG Q L. Mathematical modeling of large floating roof reservoir temperature arena[J]. Polish Journal of Chemical Technology， 2018， 20（1）： 67-74.

[17] SALOMONE T， PIOMELLI U， DE S G. WALL-MOD-SALOMONE T， et"al. Wall-modeled and hybrid large-eddy simulations of the flow over roughness strips[J]. Fluids， 2022， 8（1）： 10.

[18] SUN W， QINGLIN C， PEIDI W， et al. Effect of periodic boundary conditions on heat transfer process of crude oil in floating roof storage tank， 2016

[19] 李文欣. 強風作用下大型外浮頂罐密封圈油氣耗散研究[D]. 北京： 中國石油大學（北京）， 2019.

[20] SHARMA Y K， MAJHI A， KUKRETI V S， et al. Stock loss studies on breathing loss of gasoline[J]. Fuel， 2010， 89（7）： 1695-1699.

[21] HASSANVAND A， HASHEMABADI S H， BAYAT M. Evaluation of gasoline evaporation during the tank splash loading by CFD techniques[J]. International Communications in Heat and Mass Transfer， 2010， 37（7）： 907-913.

[22] WANG Z L，HUANG W Q，JI H. Numerical simulation of vapor diffusion and emission in loading gasoline into dome roof tank[J]. Acta. Petrol. Sin.，2017， 33：26-33.

[23] HUANG W Q， FANG J， LI F， et al. Numerical simulation and applications of equivalent film thickness in oil evaporation loss evaluation of internal floating-roof tank[J]. Process Safety and Environmental Protection， 2019， 129： 74-88.

[24] JING HB，HUANG W Q. Study on similarity criteria number of wind tunnel experiment for oil vapor diffusion based on numerical simulation technology[J]. Changzhou Univ.，2019， 31：25-34.

[25] 黃維秋， 盧建國， 趙書華. 油罐注油蒸發(fā)損耗的實驗測定[J]. 實驗技術(shù)與管理， 2009， 26（1）： 19-21.

[26] 李飛， 黃維秋， 紀虹， 等. 浮盤和氣窗位置對內(nèi)浮頂罐正己烷蒸發(fā)損耗的影響[J]. 環(huán)境工程學報， 2018， 12（2）： 410-416.

[27] 王建偉.立式儲油罐油品蒸發(fā)損耗及外干涉流場研究[D]. 西安：西安石油大學，2020.

[28] WANG M， YU B， ZHANG X Y， et"al. Experim-ental and numerical study on the heat transfer character-istics of waxy crude oil in a 100，000 m3 double-plate floating roof oil tank[J]. Applied Thermal Engineering， 2018，136.

[29] LI W， SHAO Q Q， LIANG J. Numerical study on oil temperature field during long storage in large floating roof tank[J]. International Journal of Heat and Mass Transfer，2019，130.

[30] 張佩宇， 黃維秋， 景海波， 等. 拱頂罐內(nèi)油品蒸發(fā)損耗的數(shù)值模擬和實驗研究[J]. 安全與環(huán)境工程， 2020， 27（4）： 143-151.

[31] GINESTET S， LE B C. Evaporation flow assessment from pet-roleum product storage tanks exposed to fire conditions[J]. Oil amp; Gas Science and Technology，2018，73.

[32] HUANG W Q， WANG S， JING H B， et al. A calculation method for simulation and evaluation of oil vapor diffusion and breathing loss in a dome roof tank subjected to the solar radiation[J]. Journal of Petroleum Science and Engineering， 2020， 195： 107568.

[33] 周寧， 陳力， 呂孝飛， 等. 大型LNG儲罐連續(xù)泄漏擴散過程影響因素的分析[J]. 化工進展， 2019， 38（10）： 4423-4436.

[34] FANG J， HUANG W Q， HUANG F Y， et"al. Inv-estigation of the Superposition Effect of Oil Vapor Leakage and Diffusion from External Floating-Roof Tanks Using CFD Numerical Simulations and Wind-Tunnel Experiments[J]. Processes，2020，8（3）.

[35] KOUNTOURIOTIS A， ALEIFERIS P G， CHARALAMBIDES A G. Numerical investigation of VOC levels in the area of petrol stations[J]. Science of the Total Environment， 2014， 470： 1205-1224.

[36] CHEN S R， LIN K X， PAN A Q， et al. Study of a passive radiative cooling coating on chemical storage tanks for evaporative loss control[J]. Renewable Energy， 2023， 211： 326-335.

[37] 劉風章. 恩威爾特CC-100復合型絕熱涂層在吐哈油田原油儲罐上的應用研究[J]. 石油石化節(jié)能， 2011， 1（1）： 4-6， 15.

Research Progress oin Evaporation Loss of Oil in Storage Tanks

ZHANG Guang-jian¹，"ZHAO Dong-xu²，"QIN Li-guang²，"YU Bing²，"LI Yao²，"WANG Wei-qiang¹

（1."College of Petroleum Engineering， Liaoning Petrochemical University， Fushun Liaoning 113001， China;

2."Sinopec Marketing Liaoning Branch， Shenyang Liaoning 110031， China）

Abstract："Based on the national strategic goal of \"carbon peaking and carbon neutrality\"， this paper analyzes the existing research methods for oil and gas evaporation in storage tanks"were"analyzed，"summarizes"the application scopes"of on-site measurement， wind tunnel experiments，"and numerical simulation methods"were"summarized， elaborates on the research objects and advantages and disadvantages of different simulation methods"were"discussed， proposes the main influencing factors and improvement measures for oil evaporation loss"were"proposed， and"points"out"the importance of environmental temperature and solar radiation on its impact"was"pointed"out， provideing"improvement ideas and references"for optimizing oil storage equipment and improving consumption reduction efficiency.

Key words："Storage tank;"Evaporation loss;"Field measurement;"Wind tunnel experiments;"Numerical simulation