劉玉晗，劉寶玉，蘇 鵬

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

基于ANSYS的LNG管道保冷結(jié)構(gòu)分析

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

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

介紹了LNG管道保冷結(jié)構(gòu)以及保冷層厚度的計算方法，應(yīng)用有限元分析軟件對LNG低溫長輸管道進行了穩(wěn)態(tài)熱力分析。在計算管道截面溫度場時，采用保冷層材料導(dǎo)熱系數(shù)隨溫度變化和平均導(dǎo)熱系數(shù)兩種方法進行模擬。研究結(jié)果表明溫度分布基本相同，而熱流密度值不同。

LNG低溫管道；溫度分布；數(shù)值模擬；保冷

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

1 LNG管道保冷結(jié)構(gòu)及保冷層厚度計算方法

1.1 保冷結(jié)構(gòu)

LNG低溫管道的保冷結(jié)構(gòu)如圖1所示，由內(nèi)而外主要有：防銹層，過渡層，保冷層，防潮層以及外保護層[2]。

圖1 LNG管道保冷結(jié)構(gòu)徑向剖面圖Fig.1 LNG pipeline insulation structure profiles

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

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

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

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

1.2 保冷層厚度計算方法

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

（1）外表面溫度法

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

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

Di—絕熱層內(nèi)直徑，m；

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

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

t, ta, ts—分別為介質(zhì)（管道的外表面）、環(huán)境、外表面溫度，℃。

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

式中：td—當(dāng)?shù)貧庀髼l件下最熱月的露點溫度，℃。

（2）允許冷損失法

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

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

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

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

（3）保冷層經(jīng)濟厚度法

保冷層經(jīng)濟厚度法計算公式為：

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

T—運行時間，h；

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

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

S—絕熱工程投資年分?jǐn)偮省?/p>

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

2 管道溫度分布模擬

本算例是基于文獻[5]中實例進行的模擬計算分析。1條約6 km的LNG長輸管道，輸送壓力1.194 MPa，管道材質(zhì)選取304 L不銹鋼（0Cr18Ni9），導(dǎo)熱系數(shù)為11.8 W/(m·K)，管道外徑159 mm，管道壁厚4.5 mm，管道架空，由于當(dāng)?shù)刈顭嵩缕骄鄬穸?3%，容易在管道外表面結(jié)露，為了防止結(jié)露滲水發(fā)生所以保冷層厚度的計算方法采用外表面溫度法，保冷材料按照LT/Armaflex的平均導(dǎo)熱系數(shù)0.015W/(m·K)進行計算，得到保冷層厚度為145 mm，管道輸送LNG介質(zhì)溫度為-158 ℃，當(dāng)?shù)叵募韭饵c溫度26 ℃，環(huán)境溫度按照28.4 ℃計算，隨溫度變化的導(dǎo)熱系數(shù)值如圖2所示，管外壁至大氣的放熱系數(shù)為23.753 W/(m2·℃)，此處放熱系數(shù)包含管道外壁與大氣的對流換熱和輻射過程，在計算過程中將輻射視為對流換熱一部分。

圖2 LT/Armaflex導(dǎo)熱系數(shù)與溫度的關(guān)系曲線Fig.2 Relationship curve of LT/Armaflex conductivity coefficient with temperature

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

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

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

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

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

3 結(jié) 論

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

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

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

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

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

［5］ 李慶杰. LNG管道輸送技術(shù)[D]. 撫順：遼寧石油化工大學(xué)石油天然氣工程學(xué)院，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-），女，遼寧凌源人，遼寧石油化工大學(xué)碩士研究生，LN G管道及設(shè)備保冷專業(yè)。E-m ail：okliuyuhan@163.com。