摘 要

摘　要：為了防止LNG全容罐溫度升高、LNG沸騰氣化導致罐內壓力升高，常采用安全閥釋放罐內氣體的方法來控制壓力。在LNG蒸氣釋放過程中，若被點燃會發生火災，LNG全容罐在火災產生的

摘　要：為了防止LNG全容罐溫度升高、LNG沸騰氣化導致罐內壓力升高，常采用安全閥釋放罐內氣體的方法來控制壓力。在LNG蒸氣釋放過程中，若被點燃會發生火災，LNG全容罐在火災產生的溫度應力和LNG氣壓等荷載作用下，有可能開裂、坍塌導致LNG泄漏，進而造成重大的安全事故。目前國內尚無相關規范規定，而國外的標準、規范也僅給出了一般性的條款，沒有具體的計算公式。為此，闡述了安全閥火災熱輻射研發的必要性；探討和建立了LNG全容罐的傳導、強制對流、火災輻射和內部空間自輻射，以及相互耦合的復雜三維瞬態熱分析有限元模型，獲得了瞬態熱分析各個燃燒時間點的計算結果，并給出了合理的解釋；最后將熱學模型轉換為力學模型，施加火災產生的溫度，組合氣壓和自重，經計算獲得了LNG全容罐在火災燃燒時限內的最大應變。結果表明，與混凝土的極限應變對比，儲罐不會發生壓碎坍塌。該研究成果為后續國內規范的編制和LNG全容罐結構優化設計奠定了基礎。

關鍵詞：LNG全容罐  傳熱邊界條件  安全閥火災  瞬態熱分析  熱力耦合  最大壓應變

Transient thermal-mechanical coupling analysis of full-containment LNG tanks with relief valves against fire and radiation

Abstract：In order to prevent the rising temperature or the increasing pressure caused by LNG boiling gasification in a full-containment LNG tank，pressure control is usually done through a relief valve releasing the LNG vapor．In this process，if such easily ignited fire occurs，the resulted heat stress，LNG vapor pressure and other load conditions will lead to crackin9，collapsin9，or spilling of the LNG tank，thus resulting in some major accidents．Unfortunately，no related codes or regulations can be found yet in China and in some foreign codes only general clauses but no specific calculating formula can be found at present．In view of this，we first explained the necessity of this study，and then established a complicated 3D transient thermal analysis FEA model，reflecting the conduction，forced convection，fire radiation，and self radiation in the inner Space．With this model，the calculation results were thus obrained at each burning time point，and thus reasonable explanation were given．Finally，we transformed this thermal model into a mechanical one，and with the fire temperature exerted，LNG vapor pressure，and self weight as the applied load，we calculated the maximum strain of a full containment LNG tank，which is lower than its limiting strain．This study provides solid foundation for the subsequent work like the programming of relevant domestic codes，or optimal structural design of full-containment LNG tanks．

Keywords：full containment l．NG tank，heat transfer boundary conditions，fire relief valve，transient thermal analysis，thermal-mechanical coupling analysis，maximum compressive strain

隨著石油資源的日趨緊張，LNG的應用前景越來越廣闊^[1-4]。LNG儲罐是整個LNG接收站的核心，其結構形式主要有單壁型、雙壁雙頂型、雙壁單頂敞口型和全容式。其中全容罐具有容量大、結構性能好的優點而得到廣泛應用^[5-6]。本文瞬態熱力學分析的有限元模型為LNG全容罐(以下簡稱儲罐)。

在外界溫度變化下，儲罐內的LNG達到沸點蒸騰，從而產生氣壓。為保護混凝土外罐免受不斷增加的氣壓作用，需通過穹頂上部的安全閥將超壓的LNG蒸氣排出。在此過程中，LNG蒸氣或被點燃發生火災。在這種偶然工況下，國內現有的混凝土規范尚無規定，而國外的ENl473等規范僅給出一般性的條款，沒有具體的計算公式。所以發生安全閥火災時，儲罐各部位的溫度怎么分布?火災導致的混凝土熱應變是否達到了極限應變等，這些需借助大型有限元軟件進行詳細的熱力學分析才能得到，指導后續設計工作。

筆者基于熱力學基本理論和大型有限元軟件Ansys，考慮內部材料的保溫熱傳導、外部有風情況下的強制對流以及安全閥火災的熱輻射作用，進行基于時間進程的瞬態熱分析，獲得儲罐溫度場云圖分布，然后通過熱—固耦合，將熱學模型轉換為固體力學模型，獲取儲罐各個部分的力學結果，從而判斷儲罐在火災作用下的安全性，對國內相關規范的編制具有一定的意義。

1　安全閥火災輻射熱模型

1．1　熱力學理論^[7-11]

火災作用下儲罐的熱力學分析，需要考慮傳導、對流和輻射3種傳熱方式。熱傳導理論較為簡單，此不贅述。

1．1．1強制對流

對流有自由對流和強制對流。自由對流和溫度緊密相連，而強制對流還需考慮風速的影響。牛頓依據固體與周圍流體之間的熱量交換，提出了牛頓冷卻公式：

Q＝h(tw-tf)            (1)

式中q為熱流密度，W／m²；h為對流換熱的表面傳熱系數，W／(m²·K)；tw為固體表面溫度，℃；tf為周圍流體溫度，℃。

對流換熱系數的計算需要考慮固體表面與流體之間邊界層的導熱性質。

邊界層可以理解為通過熱傳導傳遞熱量：

 

結合以上兩式，可以得到：

 

結合本文參考文獻[10]可得：

 

所以

 

式中ReL為雷諾系數，Pr是普朗特常數；N為努賽爾系數，與雷諾和普朗特常數相關；k為空氣熱導率；L為幾何體的有效長度。

通過式(4)和(5)將對流換熱系數(h)的問題，轉為解決努賽爾系數的問題。

1．1．2輻射^[12]

熱輻射指物體發射電磁能，并被其他物體吸收轉變為熱的熱量交換過程。熱輻射無須任何介質，在真空中熱輻射效率最高。

在工程中通常考慮兩個或兩個以上物體之間的輻射，系統中每個物體同時輻射并吸收熱量。它們之間的凈熱量傳遞可以用斯蒂芬～波爾茲曼方程來計算：

Q=esAiFij(Ti⁴-Tj⁴)                (6)

式中q為熱流率；e為實際物體的輻射率，或稱為黑度，它的數值處于0～l之間；s為斯蒂芬—波爾茲曼常數，約為5.67×10^-8w／(m²·K⁴)；Ai為輻射面i的面積；Fij為由輻射面i到輻射面j的形狀系數(又稱角系數)，是關于表面面積、面的取向及面間距離的函數，它的計算是熱輻射分析中的重點和難點；Ti為輻射面i的絕對溫度；Tj為輻射面j的絕對溫度。

可以看出，熱輻射引起的熱流與物體表面絕對溫度的4次方成正比，所以包含輻射的熱分析是高度非線性的。

在儲罐熱分析中要考慮兩部分輻射，即內部空間的相互輻射和火災燃燒輻射。其中內部環境的相互輻射是罐內吊頂和穹頂之間的輻射，有3000多個輻射面，角系數的計算工作量非常大。

1．2　有限元模型^[13]

儲罐瞬態熱力分析模型包含穹頂和混凝土外罐墻。局部還有結構的鋼內襯板、玻璃纖維彈性氈、珍珠巖保溫材料和鋼筋混凝土承壓環。其穹頂厚度約為400mm。

其中承壓環及襯板用shelll31單元建立、其他部分可用實體單元solid70。有限元模型求解的是連續體問題，但在此模型中，安全閥的位置離穹頂約16m，是非連續孤立的節點。為了能讓安全閥的熱輻射能傳遞到穹頂。必須借助一類特殊的單元即表面效應單元，實現孤立火災點和儲罐各個部分的熱輻射。

有限元模型見圖1、2，其中圖2是儲罐內部的情況，紅色是儲罐內部的保溫材料。儲罐節點組成見圖3。

 

 

 

穹頂表面覆蓋了表面效應單元，通過此特殊單元，實現非連續介質有限元模型熱輻射的傳遞，見圖4。

 

1．3　材料熱工性質

儲罐中部分材料的物理性質見表1。

 

1．4　瞬態熱分析的邊界條件

安全閥的熱輻射條件是根據相關國外規范^[14]，熱輻射能量為32kW／m²，安全閥高度離穹頂約l6m，持火燃燒6h。火災熱輻射是瞬態熱分析，即熱輻射隨著時間發生的。為了建立有限元模型的初始條件，在火災將要發生的ls時，執行穩態熱分析，得到儲罐的溫度分布情況。然后從1s～6h考慮安全閥火災的瞬態熱輻射，給安全閥孤立節點施加熱輻射能量，使儲罐獲得相應的背景輻射溫度。模型還要考慮通風情況下的強制對流，施加根據努賽爾系數計算得出的對流換熱系數，并考慮外部環境溫度的約束條件。

1．5　瞬態熱分析的結果

1s時的溫度分布可以看到穹頂已經獲得了初始的熱邊界條件，穹頂外表面是夏季最熱天的環境溫度41℃(圖5)。

 

隨著火災燃燒時間的持續增加，在火災發生lh后，穹頂外表面的溫度已經升高到了373℃(圖6)。

 

火災發生6h時后，穹頂的外表面溫度已經升高到了500℃左右，同時高溫區域也增大了很多(圖7)。

 

穹頂混凝土的外表面溫度無論是徑向還是環向，最高溫度均在450℃左右。內外表面溫度相差l00℃左右，將產生較大的溫度應力(圖8、9)。

 

 

2　安全閥火災輻射力模型

2．1　有限元模型

利用Ansys的“ETCHG，TTS”命令，將安全閥火災輻射的熱模型轉換為力學模型。熱單元solid70和shelll31分別轉換為solidl85和shelll81，而整個模型的幾何尺寸并沒有發生改變。

2．2　材料物理性質

表2是部分材料隨著溫度變化的物理性質。

 

2．3　荷載與邊界條件

最底部的節點固結約束；珍珠巖和纖維玻璃由于其彈性模量小，不考慮它對數值模型剛度的貢獻，約束其全部自由度，在Ansys中讀入熱分析模型的溫度場分布，施加溫度荷載。然后施加儲罐的自重、LNG的蒸氣壓進行求解。

2．4　安全閥火災力學分析結果

施加各種荷載后進行求解，獲得混凝土的壓縮應蠻(圖10)。

 

可以看到，模型中最大的壓應變是-0.002，沒有達到混凝土的極限壓應變0.0033，所以穹頂不會發生混凝土壓碎破壞，滿足儲罐安全性要求。

3　結論

1)通過傳熱學理論，闡述了3大傳熱方式的基本理論和計算方法以及怎么合適地施加到有限元模型上。

2)建立了LNG全容罐的瞬態熱分析有限元模型，獲得了熱分析中各個部分以及各個時間段上的溫度結果，此溫度分析結果是后續力學分析的基礎。

3)通過熱一固耦合，將熱學模型轉換為力學模型，并施加熱荷載、自重等荷載在模型上，獲得了儲罐混凝土穹頂各個部分的熱應變值。通過此結果，可以證明在安全閥火災熱輻射下，LNG全容罐是安全的。

4)研究結果為后期LNG全容罐結構的優化和國內規范的編制工作打下了較好的基礎。

 

參考文獻

[1]羅偉中，鄭洪搜，孟勐．美國出口LNG及對世界LNG市場的影響分析[J]．天然氣工業，2012，32(6)：93-99．

LUO Weizhon,ZHENG Hongtao，MENG Meg．LNG exports from the United States and their impact on the global LNG market[J]．Natural Gas Industry，2012，32(6)：93-99．

[2]李健胡，蕭彤．日本LNG接收站的建設[J]．天然氣工業，2010，30(1)：109-115．

LI Jianhu．XIAO Tong．Experiences learned from Japan’s LNG receiving terminals[J]．Natural Gas Industry，2010，30(1)：109-115．

[3]黃群，夏芳．LNG儲罐國產化的可行性[J]．天然氣工業，2010，30(7)：80-86．

HUANG Qun，XIA Fang．A feasibility of China made LNG storage tanks[J]．Natural Gas Industry，2010，30(7)：80-86．

[4]劉小麗．中國天然氣市場發展現狀與特點[J]．天然氣工業，2010，30(7)：1-6．

LIU xiaoli．Present situation and characteristics of the development of natural gas market in China[J]．Natural Gas Industry，2010，30(7)：1-6．

[5]付子航，單彤文．大型l．NG儲罐完整性管理初探[J]．天然氣工業，2012，32(3)：86-93．

FU Zihang，SHAN Tongwen．Integrity management system for large LNG tanks[J]．Natural Gas Industry，2012，32(3)：86-93．

[6]王冰，陳學東，王國平．大型低溫LNG儲罐設計與建造技術的新進展[J]．天然氣工業，2010，30(5)：108-112．

WANG Bing，CHEN Xuedong，WANG Guoping．Design of large low-temperature LNG storage tanks and new progress in its construction technology[J]．Natural Gas Industry，2010，30(5)：108-112．

[7]章偉星，李科浚，周吳，等．薄膜式LNG運輸船溫度場研究[J]．天然氣工業，2005，25(10)：110-112．

ZHANG Weixing，LI Kejun，ZHOU Hao，et al．Study on temperature field for membrane type LNG carriers[J]．Natural Gas Industry，2005，25(10)：110-112．

[8]彭明，丁乙．全容式LNG儲罐絕熱性能及保冷系統研究[J]．天然氣工業，2012，32(3)：94-97．

PENG Ming，DING Yi．Thermal insulation performance and cold insulation system of full containment LNG storage tanks[J]．Natural Gas Industry，2012，32(3)：94-97．

[9]陶文栓．數值傳熱學[M]．西安：西安交通大學出版社，2001．

TAO Wenshuan．Numerical heat transfer[M]．Xi¢an：Xi¢an Jiaotong University Press，2001．

[10]THEODORE，ADRIENNE，FRANK P，et al．Funda mentals of heat and mass transfer[M]．New Jersey：John Wiley&Sons Inc，1996．

[11]陶文栓，楊世銘．傳熱學[M]．北京：高等教育出版社，2006．

TAO Wenshuan，YANG Shiming．Heat transfer[M]．Beijing：Higher Education Press，2006．

[12]余其錚．輻射換熱原理[M]．哈爾濱：哈爾濱工業大學出版社，2000．

Yu Qizheng．Radiation heat transfer principle[M]．Harbin：Harbin Institute of Technology Press，2000．

[13]王勖成．有限單元法[M]．北京：清華大學出版社，2003．

WANG Xucheng．Finite element method[M]．Beijing：Tshinghua University Press，2003．

[14]European Committee for Standardization．ENl473 Installation and equipment for liquefied natural gas-Design of onshore installations[S]．London：Britain Standardization Institution，2007．

 

本文作者：余曉峰  王松生  蘇軍偉

作者單位：中石化洛陽工程有限公司