大型LNG儲罐儲液翻滾特性的多階段研究

王　萍，彭文山，曹學文，王　慶，靳學堂.勝利油田森諾勝利工程有限公司，山東東營 57000.中國石油大學（華東）儲運與建筑工程學院，山東青島 66580

大型LNG儲罐儲液翻滾特性的多階段研究

王萍1，彭文山2，曹學文2，王慶2，靳學堂2
1.勝利油田森諾勝利工程有限公司，山東東營 257000
2.中國石油大學（華東）儲運與建筑工程學院，山東青島 266580

一般認為相鄰兩層LNG間密度差大于0.5 kg/m3、溫度差大于0.2℃就會引起LNG翻滾，導致大量LNG快速蒸發，罐內壓力急劇升高影響存儲安全，不僅浪費資源且造成極大安全隱患。在罐內LNG液體已經形成分層且可引發翻滾的基礎上，利用Fluent軟件建立罐內LNG翻滾過程中氣相空間數值計算模型。得出主要結論如下：（1）儲罐內LNG分層后翻滾過程可以分為四個過程：界面擾動階段、擾動發展階段、劇烈翻滾增壓階段和平穩恢復階段；（2）儲罐內LNG發生翻滾時，罐內LNG迅速蒸發，罐內壓力急劇上升至儲罐最大工作壓力時，安全泄放閥打開泄放氣體以防止罐內超壓；（3）研究表明16萬m3儲罐充滿率為70%、上層重質LNG厚度為1 m、LNG分層間密度差為1 kg/m3時，翻滾發生時的平均蒸發速率是靜態蒸發速率的32.3倍，安全泄放閥在翻滾發生約16 h后打開泄壓。

LNG儲罐；分層與翻滾；翻滾特性；速度；壓力

隨著國家能源結構的調整以及環保要求的提高，液化天然氣作為一種無色、無味、無毒、無腐蝕性的清潔能源在國內應用越來越廣泛［1］。LNG儲罐是LNG接收站的核心設備也是LNG存儲的主要設備［2-3］。LNG在罐內存儲過程中受到外部熱源的侵擾而發生靜態蒸發，LNG受熱蒸發罐內壓力升高，影響儲罐的安全存儲；上層LNG液體中沸點較低的組分（如氮、甲烷等）優先吸熱蒸發［4-5］，形成上下兩層或多層的密度、溫度等不同的LNG分層；新舊LNG的混裝也會引起罐內LNG分層。當相鄰兩層LNG之間的密度差、溫度差達到一定值則會導致其翻滾［6-7］，并大量快速蒸發［8］，造成極大的安全威脅和經濟損失。

本文以16萬m3LNG儲罐為例，在罐內LNG已經分層且會發生翻滾的基礎上，建立LNG翻滾過程中氣相部分數值模型，利用Fluent軟件進行計算，研究翻滾過程中氣相空間運動規律和壓力變化規律。這是首次建立LNG儲罐內儲液翻滾時氣相空間的研究模型，對罐內的翻滾過程進行分階段分析，對罐內翻滾過程中的壓力變化進行研究，為生產中的翻滾預防提供理論支持。

1　模型建立

1.1物理模型

以16萬m3LNG儲罐為原型，建立密閉大型LNG儲罐分層與翻滾的物理模型，對翻滾過程中氣相空間運動規律和壓力變化進行研究。LNG儲罐內的LNG液體發生分層是儲罐翻滾的必要條件。圖1所示為儲罐內的LNG液體已發生分層與翻滾。

圖1　L NG儲罐分層與翻滾示意

當兩層LNG之間的密度差、溫度差達到一定值時，下層LNG與上層LNG之間發生劇烈的傳熱傳質現象。下層LNG上升到頂面，由于壓力的降低使得其成為過熱飽和液體，蓄積的能量急劇釋放，LNG蒸發產生大量的BOG氣體，儲罐內的壓力急劇上升，罐內發生翻滾現象。

1.2假設條件

根據所研究的內容，進行了以下的假設：

（1）以16萬m3大型LNG儲罐為例，建立三維模型，直徑取80 m，忽略內罐壁厚。

（2）儲罐的漏熱已經造成了LNG液體分層，分層形式為上重下輕結構。

（3）LNG儲罐與外界環境之間的熱交換全部用于液相蒸發和溫度提升。

（4）儲罐內液相、氣相部分的溫度是均勻的。

（5）整個儲罐的內部空間分為三部分：氣相部分、液相主體部分和氣液交界面部分，其中氣液交界面部分只考慮其存在而忽略其厚度。

1.3數學模型

建立儲罐內部翻滾的數學計算模型。

（1） 氣相空間質量守恒方程：

式中：τ為時間，s；ρg為氣相密度，kg/m3；vg為氣相體積，m3；Gw為單位時間儲罐內壁凝結的氣體質量，kg/s；Gf為單位時間氣液界面上凝結的氣體質量，kg/s。

（2） 氣相空間能量守恒方程：

式中：hg為飽和氣體焓，J/kg；αwg為氣相空間與儲罐內壁間對流換熱系數，J/（m2·K）；Fwg為氣相空間與儲罐內壁接觸面積，m2；Twg為氣相空間壁面溫度，K；Ts為儲罐內壓力對應的飽和溫度，K；hs為飽和溫度Ts對應的飽和氣體焓，J/kg；αf為氣相空間與氣液界面對流換熱系數，J/（m2·K）；Ff為氣液界面面積，m2；Tg為氣相空間溫度，K。

（3）氣相空間內壁面單位時間氣體冷凝質量：

本文的研究假設氣相空間的氣體溫度均勻一致，氣液界面只考慮其存在而忽略其厚度，則可以認為在氣相空間的儲罐壁面上氣體冷凝量是均勻的，則氣相空間儲罐壁面單位時間氣體冷凝質量為：

式中：λs為氣體冷凝膜導熱系數，W/（m·K）；δ為氣體冷凝膜的厚度，m；γs為飽和溫度Ts對應的氣化潛熱，J/kg。

（4）氣液界面單位時間氣體冷凝質量：

式中：Tl為液相溫度，K；λl為液相導熱系數，W/ （m·K）；Cpl為液相比熱容，J/（kg·K）；ρs為飽和溫度Ts對應的氣相密度，kg/m3；ρl為液相密度，kg/m3。

（5） 液相質量守恒方程：

式中：vl為液相體積，m3。

（6） 液相能量守恒方程：

式中：hl為液體焓，J/kg；αw l為液相與壁面間對流換熱系數，J/（m2·K）；Fw l為液相與儲罐內壁接觸面積，m2；Tw l為液相壁面溫度，K。

（7） 體積守恒方程：

式中：V為儲罐總容積，m3。

（8） 氣體狀態方程：

式中：P為儲罐內壓力，Pa。

1.4模擬計算

本文采用計算流體動力學方法（Computational Fluid Dynamics，CFD）建立大型LNG儲罐內分層與翻滾模型。針對大型LNG儲罐內LNG液體在分層形成之后的翻滾過程進行數值模擬，模擬時選用壓力基（Pressure-Based）求解器以及瞬態（Transient）模型，模擬選用mixture模型、k-ε方程，儲罐頂部設為pressure-out類型，其他邊界條件設為wall類型。在模型的設置中，相間作用選擇 evaporation-condensation模型。

依據以上的物理模型、數學模型以及數值計算方法的分析，利用Fluent軟件進行模型的求解。儲罐初始充滿率為70%，初始壓力為101.325 kPa，罐頂出口背壓為130.325 kPa；罐內液相主體為輕質LNG，其上層覆蓋1 m厚重質LNG，氣相空間的初始時刻充滿甲烷。研究采用的介質物性見表1。

選用三維數值模擬，在計算區域中采用結構化網格進行網格劃分，網格劃分數目為2 193 457。圖2 為LNG儲罐的網格劃分圖。

表1　L NG和甲烷物性

圖2　L NG儲罐網格劃分

2　結果分析

LNG儲罐內液體發生分層后，其密度差、溫度差達到一定值時罐內會發生LNG層間的翻滾。通過對整個過程的模擬研究分析，罐內的翻滾過程可以分為四個階段：界面擾動階段、擾動發展階段、劇烈翻滾增壓階段和平穩恢復階段。利用三維模型研究LNG儲罐分層與翻滾的氣相空間運動規律及壓力變化規律，以下取儲罐中心剖面進行結果分析。

2.1界面擾動階段

圖3為儲罐內LNG液相在分層后發生翻滾時界面擾動階段速度云圖。

從圖3（a）可以看出，儲罐內的翻滾現象從罐壁開始，這是因為壁面無滑移，且罐壁與外界環境之間存在熱交換，分層界面首先在壁面處發生彎曲開始翻滾。從圖3（b）可以看出LNG分層間的翻滾現象由罐壁向罐中心發展，最大速度出現在氣液交界面處，這是因為LNG層間翻滾而蒸發大量氣體，氣液交界面處是蒸發氣體的發源地，因而此處的氣體速度最大。從圖3（b）、（c）、（d）可以看出，在罐中心處的氣相速度較小，形成一個旋渦，且以中心軸線為分界；兩側速度方向相反，背向而行。這是由于罐內LNG層間翻滾的發生是由罐壁向罐中心發展的，在氣液交界面處貼近液相的一層氣體速度方向是由罐壁到罐中心，兩側氣體在罐中心相遇且相互沖撞排斥，由于對稱性使得氣體在罐中心處向上運動。儲罐頂部壓力較大，罐中心兩側的氣體向上運動一定距離后方向改變，返回罐壁方向。這樣在液相上方，罐中心兩側各自形成一個氣體旋渦。

圖3　L NG儲罐內翻滾時界面擾動階段速度云圖（單位：m/s）

2.2擾動發展階段

圖4為儲罐內翻滾的擾動發展階段的速度云圖。

從圖4可以看出，儲罐內氣體的最大速度呈現增長趨勢，從圖4（a）中的1.44 m/s增加到圖4（d）中的1.78 m/s。由于儲罐內翻滾的發生，大量的LNG蒸發成為氣體，使得儲罐內氣體速度呈現增長趨勢。從圖中可以明顯看出，儲罐內的氣體由氣液界面向罐頂方向擴展，氣體速度及波動范圍越來越大。這也說明，儲罐內的翻滾處于擾動發展階段時，不僅表現在氣體速度的增長，也表現在速度波動范圍的增大趨勢。從圖4中還可看出，氣體旋渦在罐中心處兩側反向氣體的作用下，氣體旋渦向罐壁方向運動。

2.3劇烈增壓階段

圖5為儲罐內翻滾的劇烈增壓階段的速度云圖。

從圖5可以看出，氣體旋渦中心處具有最大的氣相速度。在氣體旋渦由罐中心向罐壁運動以及由罐壁向罐中心運動的過程中，氣體具有的最大速度，比氣體旋渦與儲罐罐壁碰撞以及兩氣旋在罐中心處碰撞的速度要大。

在氣體旋渦由罐中心向罐壁運動以及由罐壁向罐中心運動的過程中，氣體具有的最大速度約為1.1 m/s；在氣體旋渦與罐壁碰撞以及兩氣旋在罐中心處碰撞過程中，氣體具有的最大速度約為0.7 m/s。這是由于氣旋與罐壁之間的作用以及氣旋碰撞的相互作用，使得氣旋速度降低。

圖6為儲罐內翻滾劇烈增壓階段的另外一次碰撞循環過程。

從圖6（a）可以看出，在氣液界面上方的氣相空間里，儲罐中心軸線兩側的氣體旋渦均由罐中心向罐壁方向運動，與罐壁碰撞后，部分氣體沿罐壁向上運動至儲罐頂部，剩余氣體上下翻滾形成旋渦，旋渦的旋轉方向不改變。從圖6（b）可知，罐壁的碰撞作用使氣體旋渦整體運動方向改變，由罐壁返回罐中心。圖6（c）顯示的是儲罐中心處的氣體由罐頂向下運動，并在中心兩側再次形成氣體旋渦，但旋渦的旋轉方向與碰撞之前相反。圖6（d）顯示的是氣體旋渦在儲罐的中心碰撞中，由于相互作用力的作用，使得氣體旋渦整體改變運動方向，再次向儲罐的罐壁移動。LNG儲罐內液相之間的密度差、溫度差，以及儲罐與外界環境之間的熱交換為翻滾中氣體的蒸發提供了氣旋持續運動的動力。這正是儲罐內的劇烈翻滾增壓階段。

圖4　L NG儲罐內翻滾時擾動發展階段速度云圖（單位：m/s）

圖5　L NG儲罐翻滾時劇烈增壓階段速度云圖（單位：m/s）

圖6　L NG儲罐翻滾時劇烈增壓階段速度云圖（單位：m/s）

2.4平穩恢復階段

當LNG儲罐內液相之間的密度差、溫度差消除后，儲罐與外界環境之間熱交換時儲罐處于靜態蒸發狀態。儲罐內的LNG翻滾逐漸減弱直至平穩，此時為翻滾的平穩恢復階段，儲罐內的LNG翻滾現象結束。

2.5罐內壓力研究

通過對LNG儲罐儲液分層與翻滾的模擬，對儲罐氣液界面上方氣相的壓力進行密切關注，可以得到儲罐內氣相的壓力隨時間變化的規律。該模擬的監測點為儲罐中心頂端的壓力出口處。圖7為氣相出口處的壓力隨時間的變化曲線。

從圖7可以看出，儲罐氣相出口處的壓力隨時間呈現增長趨勢，且增長速度越來越快，即儲罐內發生翻滾，大量LNG急劇蒸發，儲罐內壓力迅速增長。從圖中標注的虛線處可以推算，在分層發生翻滾到罐內壓力達到儲罐的最大工作壓力（表壓29 kPa），時間經歷約16 h。即從儲罐內發生分層與翻滾到罐內壓力達到儲罐的最大工作壓力共經歷約16 h，此時儲罐的安全泄放閥打開，釋放儲罐內過多的氣體，維持儲罐的安全。

圖7　L NG儲罐翻滾時罐氣相空間壓力變化曲線

通過計算，該工況下的儲罐罐內壓力增長率約為1.81 kPa/h，這與文獻［9］中的數據比較，約為該工況下靜態日蒸發時的罐內壓力增長率0.056 kPa/h的32.3倍。分層平衡被破壞后的LNG出現翻滾沸騰狀態，其蒸發量可達到靜態蒸發量的10～50倍甚至更高。

3　結論

翻滾過程是一種急劇的相變過程，對整個氣相空間的氣相運動和壓力變化影響巨大，嚴重威脅儲罐的安全存儲。對翻滾過程中氣相空間運動規律和壓力變化的研究，得出如下主要結論：

（1）儲罐內分層發生后，翻滾過程可以分為四個階段：界面擾動階段、擾動發展階段、劇烈翻滾增壓階段和平穩恢復階段。

（2）當儲罐內發生翻滾時，罐內LNG迅速蒸發，罐內壓力急劇上升，罐內壓力達到儲罐的最大工作壓力時，儲罐的安全泄放閥打開，防止罐內超壓。

（3）翻滾發生時，LNG的蒸發速率是正常情況下的10～50倍甚至上百倍，對于本文的研究，16萬m3儲罐充滿率為70%、在上層重質LNG厚度為1 m、LNG分層間密度差為1 kg/m3時，其翻滾的平均蒸發速率是正常情況下的32.3倍，儲罐的安全泄放閥在翻滾發生約16 h時打開。

［1］劉小麗.中國天然氣市場發展現狀與特點［J］.天然氣工業，2010，30（7）：1-6.

［2］彭文山，程旭東，林楠，等.LNG儲罐環向預應力筋設計方法［J］.低溫建筑技術，2013，35（5）：38-40.

［3］程旭東，彭文山，孫連方，等.熱應力作用下LNG儲罐外罐裂縫及失效時間分析［J］.天然氣工業，2015，35（3）：102-109.

［4］汪順華，魯雪生，趙紅霞.低溫液體容器無損存儲傳熱模型［J］.低溫工程，2001（6）：37-41.

［5］徐烈，孫恒，李兆慈，等.溫度與充滿率對低溫容器無損貯存性能的影響［J］.化工學報，2001，52（10）：891-895.

［6］王武昌，李玉星，孫法峰，等.大型LNG儲罐內壓力及蒸發率的影響因素分析［J］.天然氣工業，2010，30（7）：87-92.

［7］喬國發.影響LNG儲存容器蒸發率因素的研究［D］.青島：中國石油大學（華東），2007.

［8］候金偉.LNG儲罐分層與翻滾的數值模擬［D］.哈爾濱：哈爾濱工業大學，2009.

［9］王小尚，王萍，曹學文，等.大型LNG儲罐安全存儲影響因素分析［J］.石油工程建設，2015，41（5）：5-11.

國家自然科學基金（51274232）和中央高校基本科研業務費專項資金（13CX06074A）資助。

Multi-stage Studyon Rolling Characteristics ofL iquid in L arge L NGStorage Tank

WANG Ping1，PENG Wenshan2，CAO Xuewen2，WANG Qing2，JIN Xuetang2
1.ShengliOilfield Sennuo ShengliEngineering Co.，Ltd.，Dongying 257100，China
2.College of Pipeline and CivilEngineering，China University of Petroleum，Qingdao 266580，China

It's generally believed that the rolling of liquid in LNG storage tank happens as the density difference of two adjacent layers of LNG gets greater than 0.5 kg/m3and the temperature difference gets greater than 0.2℃.Large amounts of LNG evaporate quickly，which raises the pressure in the tank and shortens the secure storage time，even may lead to stratification and rolling causing massive evaporation of LNG.This phenomenon is not only a waste of resources but also a potential safety hazard.On the basis of condition that LNG liquid in the tank has already become layered and has the possibility for rolling，this paper uses Fluent software to establish numerical model for gas phase space in the rolling process of LNG in the tank.The main results are as follows：（1） Following the stratification，the process of rolling can be divided into four periods：occurrence of interfacial disturbance，development of disturbance，violent rolling with increasing pressure and smooth recovery.（2） The rolling inside the tank makes LNG evaporate at a high speed which sharply increases the pressure in the tank.When the pressure reaches the maximum working pressure of tank，the relief valve would open to release the excessive gas in case of overpressure.（3） For the conditions that the volume of tank is 160 000 m3，the filling rate is 70%，thickness of top heavy LNG is 1 m and the density difference in stratified layers is 1 kg/m3，the average evaporation rate is 32.3 times as high compared to the normal case，and the relief valve would open about 16 h later after the beginning of rolling.

LNG storage tank；stratification and rolling；rolling characteristics；velocity；pressure

10.3969/j.issn.1001-2206.2016.01.004

王萍（1989-），女，山東德州人，2015年畢業于中國石油大學（華東）油氣儲運專業，碩士，主要從事大型LNG儲罐干燥置換、分層與翻滾以及海管完整性等方面的研究。

Email：wangpflyhigh@126.com

2015-11-18