基于ANSYS的LNG管道保冷結構分析

劉玉晗，劉寶玉，蘇 鵬

（1. 遼寧石油化工大學 石油天然氣工程學院， 遼寧 撫順 113001； 2. 中海油銷售惠州有限責任公司， 廣州 惠州 516002）

基于ANSYS的LNG管道保冷結構分析

劉玉晗1，劉寶玉1，蘇 鵬2

（1. 遼寧石油化工大學 石油天然氣工程學院， 遼寧 撫順 113001； 2. 中海油銷售惠州有限責任公司， 廣州 惠州 516002）

介紹了LNG管道保冷結構以及保冷層厚度的計算方法，應用有限元分析軟件對LNG低溫長輸管道進行了穩態熱力分析。在計算管道截面溫度場時，采用保冷層材料導熱系數隨溫度變化和平均導熱系數兩種方法進行模擬。研究結果表明溫度分布基本相同，而熱流密度值不同。

LNG低溫管道；溫度分布；數值模擬；保冷

液化天然氣（Liquefied Natural Gas，簡稱LNG），是天然氣經過脫硫、脫碳、脫水、脫重烴和雜質等預處理后，再經液化工藝冷卻降溫至-162 ℃而成[1]。天然氣液化后的體積是原體積的1/625左右，與氣態相比采用液態輸送相同體積的所需輸送動力費用優勢明顯，基本建設成本低。同時，采用LNG輸送安全、環保，可以大大降低氣體輸送帶來爆炸、環境污染的可能性。由于LNG在輸送過程中溫度很低，需要采取相應的保冷措施，合理的保冷絕熱結構能夠減少冷損失和介質溫度升高，阻止保冷層表面結霜凝露，從而節約能源。通常設計計算保冷層厚度時，采用平均導熱系數進行計算，實際當中大部分材料的導熱系數隨著溫度的變化而變化。本文基于ANSYS軟件對LNG低溫管道結構進行了溫度場分析，計算對比分析了保冷材料導熱系數隨溫度變化和平均導熱系數兩種情況下的溫度分布和管道徑向熱流率情況，為LNG管道保冷結構的優化提供了相應的依據。

1 LNG管道保冷結構及保冷層厚度計算方法

1.1 保冷結構

LNG低溫管道的保冷結構如圖1所示，由內而外主要有：防銹層，過渡層，保冷層，防潮層以及外保護層[2]。

圖1 LNG管道保冷結構徑向剖面圖Fig.1 LNG pipeline insulation structure profiles

防銹層主要是防止不銹鋼管道腐蝕，但是本文以文獻[3]中管道實例為研究對象不需涂刷防銹漆；

過渡層材料主要一具有較高強度和彈性的無機纖維制品；

保冷層常用結構為捆扎式，由于管道和保冷材料的線脹系數不一樣，為避免材料收縮后應力過大導致保冷結構損壞，要求每層保冷材料每隔4 m左右設置伸縮縫，各層伸縮縫錯開，本文保冷材料選用LT/Armaflex，該制品在低溫下具有較小的導熱系數；

防潮層主要作用是為了防止外界濕氣滲透到保冷絕熱層中，產生“冷橋”使保冷材料導熱系數變大，以致破壞保冷結構，但文中所選用的LT/Armaflex材料本身透濕系數為≤1.31× 10-11g/(m· s·Pa)，達到防潮標準，所以可以不設防潮層；外保護層一般由0.25～0.50 mm厚的鋁材、鍍鋅鐵皮或不銹鋼構成主要作用是為了保護保冷絕熱材料免受外界的影響和損傷。

1.2 保冷層厚度計算方法

低溫管道保冷層的計算方法主要有外表面溫度法，允許冷損失法，保冷層經濟厚度法等[4]。本文僅對以下常用的3種計算方法進行分析。

（1）外表面溫度法

若保冷層為多層時的總厚度的計算公式為

其中：Do—絕熱層外直徑，m；

Di—絕熱層內直徑，m；

α—絕熱層外表面對大氣的放熱系數，W/(m2·℃)；

λi—絕熱層第i層材料的導熱系數，W/(m·℃)；

t, ta, ts—分別為介質（管道的外表面）、環境、外表面溫度，℃。

若管道的保冷材料僅用一種則視為單層保冷，同時防止管道外表面結露，公式（1）可簡化為：

式中：td—當地氣象條件下最熱月的露點溫度，℃。

（2）允許冷損失法

多層保冷絕熱情況下保冷層厚度計算公式：

其中：q’—允許最大冷損失值，W/m2；

K —對于冷損失控制較嚴格的管道要求對q’進行修正，一般取K=0.8～0.9；

單層保冷情況下的保冷厚度：

（3）保冷層經濟厚度法

保冷層經濟厚度法計算公式為：

其中：fn—熱能價格，元/106 kJ；

T—運行時間，h；

?t—?t=ts-ta，℃；

Pi—絕熱材料單位價格，元/m3；

S—絕熱工程投資年分攤率。

上述3種絕熱保冷層厚度的計算方法各有利弊，應該按照具體實際工程來選擇最佳的計算方法，外表面溫度法（防止表面結露）法和允許冷損失法應互為校核計算，選取符合要求的保冷層厚度。

2 管道溫度分布模擬

本算例是基于文獻[5]中實例進行的模擬計算分析。1條約6 km的LNG長輸管道，輸送壓力1.194 MPa，管道材質選取304 L不銹鋼（0Cr18Ni9），導熱系數為11.8 W/(m·K)，管道外徑159 mm，管道壁厚4.5 mm，管道架空，由于當地最熱月平均相對濕度83%，容易在管道外表面結露，為了防止結露滲水發生所以保冷層厚度的計算方法采用外表面溫度法，保冷材料按照LT/Armaflex的平均導熱系數0.015W/(m·K)進行計算，得到保冷層厚度為145 mm，管道輸送LNG介質溫度為-158 ℃，當地夏季露點溫度26 ℃，環境溫度按照28.4 ℃計算，隨溫度變化的導熱系數值如圖2所示，管外壁至大氣的放熱系數為23.753 W/(m2·℃)，此處放熱系數包含管道外壁與大氣的對流換熱和輻射過程，在計算過程中將輻射視為對流換熱一部分。

圖2 LT/Armaflex導熱系數與溫度的關系曲線Fig.2 Relationship curve of LT/Armaflex conductivity coefficient with temperature

用ANSYS計算分析LNG低溫管道的溫度分布結果如圖3所示。圖3(a)是取保冷材料導熱系數為平均數計算得到，計算過程是線性的。圖3(b)是按照保冷材料的導熱系數隨溫度變化進行的計算分析，計算的過程是非線性的。但是從計算結果可以看出溫度分布基本一致，取保冷材料導熱系數的平均值進行設計即可滿足外表面溫度法的要求。圖4所示的溫度梯度在整個徑向方向分布趨勢一樣，具體數值上差別也較小。

圖3 LNG管道徑向截面溫度分布Fig.3 Temperature distribution of LNG pipeline profiles

圖4 LNG管道徑向方向溫度梯度Fig.4 Temperature grads of LNG pipeline in the radial direction

LNG管道截面溫度能夠滿足設計要求，但是從LNG管道徑向熱流率曲線來看，發現二者隨著計算時采用的導熱系數的不同得到的結果顯示出差異。圖5(b)中輸液管管壁處的熱流率明顯大于圖5 (a)中的輸液管道管壁處的熱流率；同樣在保冷層外壁處(b)圖中的熱流率值是(a)圖中相同位置的熱流率值的2倍多，如果考慮到限制冷損失的情況應用平均導熱系數計算在某些情況下可能會不滿足條件。圖5(b)中保冷材料的導熱系數是隨著溫度的變化而變化的，所以計算出來的結果較為接近實際情況。

圖5 LNG管道徑向熱流率Fig.5 Heat flux of LNG pipeline in the radial direction

3 結 論

LNG管道保冷層設計過程中，設計條件是滿足外表面溫度時，為方便計算可以采用平均導熱系數進行計算；如果需要滿足限制最大冷損失要考慮用保冷材料導熱系數隨溫度變化進行計算。

［1］ 揚帆，陳保東，姜文全,等.液化天然氣技術在調峰領域的應用[J].油氣儲運，2006，25(10):26.

［2］ 李鈺.低溫管道的保冷設計[J].石油化工設計，2009,26(1)：52-55.

［3］ 李慶杰，陳保東，洪麗娜,等. LNG管道輸送工藝計算[J]. 當代化工，2011，40(2):163-164;167.

［4］ 張德姜，王懷義，劉紹葉. 石油化工裝置工藝管道安裝設計手冊[M]. 北京：中國石化出版社，1994.

［5］ 李慶杰. LNG管道輸送技術[D]. 撫順：遼寧石油化工大學石油天然氣工程學院，2010.

Analysis of LNG Pipeline Insulation Structure Based on ANSYS Numerical Simulation

LIU Yu-han1，LIU Bao-yu1，SU Peng2
（1. College of Petroleum Engineering, Liaoning Shihua University, Liaoning Fushun 113001, China；2. CNOOC Marketing Huizhou Co.,Ltd., Guangzhou Huizhou 516002, China）

LNG pipeline insulation structure and calculation methods of insulation thickness were introduced. Static temperature distribution of LNG cryogenic long-distance pipeline was analyzed by using finite software. During calculation of pipeline cross-section temperature distribution, two methods including changing thermal conductivity with temperature and average thermal conductivity of the same insulation material were used to simulate the temperature distribution,separately. The results show that temperature distributions simulated by above two methods are similar, but heat flow densities are different.

LNG cryogenic pipeline; Temperature distribution; Numerical simulation; Insulation

TE 832

A

1671-0460（2012）02-0185-04

2011-00-00

劉玉晗（1986-），女，遼寧凌源人，遼寧石油化工大學碩士研究生，LN G管道及設備保冷專業。E-m ail：okliuyuhan@163.com。