李慶杰,陳保東,洪麗娜,馬金鵬
(1. 遼寧石油化工大學石油天然氣工程學院, 遼寧 撫順 113001; 2. 遼河油田公司錦州采油廠,遼寧 錦州 121209;3. 中國石油管道大慶輸油氣分公司, 黑龍江 大慶 163458)
LNG管道輸送工藝計算
李慶杰1,陳保東1,洪麗娜2,馬金鵬3
(1. 遼寧石油化工大學石油天然氣工程學院, 遼寧 撫順 113001; 2. 遼河油田公司錦州采油廠,遼寧 錦州 121209;3. 中國石油管道大慶輸油氣分公司, 黑龍江 大慶 163458)
以某LNG公司建設的LNG接收站為依托背景,對該地區建設一條6 km長管道的運行參數進行了計算。利用有限容積法對管道截面的溫度場分布進行了數值計算,得出了保冷層厚度為95 mm時的外表面溫度。通過保冷層厚度和外管道直徑的選取,計算出輸送到終點的壓降和溫升,驗證了 LNG在“過冷”狀態下輸送的可行性。
液化天然氣; 管中管輸送; 保冷層; 過冷
天然氣液化后是原來體積的1/600左右,采用液態輸送較氣態相比,輸送相同體積的天然氣,LNG的管徑要小,泵的能耗也相對較低,輸送成本較為理想;相對于居民用戶來講,如果把LNG直接輸送到化工企業,這將為企業帶來巨大的冷量,節約企業制冷費用。所以采用管道輸送LNG具有巨大的優勢和發展潛力。
據相關文獻介紹,國際上實行LNG管道輸送的實例較少而在國內目前尚處于起步階段。加拿大曾在此方面做過研究,結果表明當LNG輸量越大,距離越長,地溫越低時,液化輸送的經濟效果較好。目前制約 LNG管道發展除了輸送技術及設備等因素外,在理論上仍缺少相應的研究。
管道敷設方式一般有3種形式:架空敷設、埋地敷設和管溝敷設。由于LNG溫度極低,管道輸送方式中架空較為理想,但是考慮到建設在人口較為稠密的城區時,此方案較為不可行。其次由于埋地敷設方式可以降低地面的安全隱患,減少對周邊人和物的冷損傷,可靠性較好,但是LNG的保冷層收到土壤的擠壓后保冷效果迅速下降,不利于管道的無氣化輸送,所以此方案也有一定的缺點。為了能夠解決這一問題,本文選用了管中管敷設的方式,這樣既能夠保證保冷效果又可以進行直埋敷設降低地面安全隱患。使用管道輸送時,為了保持天然氣輸送到接受終端仍以液態方式存在,需要盡量減少管道與外界的換熱量;確保LNG處于一定的壓力范圍內,使其溫度低于相應壓力下的飽和溫度,讓其處于“過冷”狀態。本文利用“過冷”態原理進行了無氣化管道輸送計算。
本文以某LNG公司為依托背景,在其附近建設一條約6 km的埋地LNG管道。已知參數如下:以進口澳大利亞的LNG成分為例(其成分如表1),氣候條件可參照當地氣象參數。在計算過程中,泵的揚程根據國外進口 LNG泵的性能提供的 1.194 MPa壓力,輸送距離為6 km,高程差48.5 m。其保冷層厚度選95 mm,管道材質選取0Cr18Ni9,其導熱系數=11.8 W/(m·K),管徑為159 mm,壁厚為4.5 mm,空氣導熱系數為 W/(m·K)。保溫材料選取LT/Armaflex(福樂斯)系列產品,根據計算得到該材質平均導熱系數為W/(m·K),泵出口溫度為-158℃。

表1 LNG組分Table 1 LNG composition
為了能夠更合理計算管道輸送過程中沿程的溫降,本文將外管道至于年平均溫度為9.3 ℃的土壤中。為了盡量減少管道對土壤的溫度場影響,并考慮到建設成本,盡量在投資較小的情況下選擇保護管道外表面溫度與地溫近似,通過計算,本文選取的參數相對較為合理。其參數選擇如下:輸送介質管道的外徑為159 mm,保冷層厚度為95 mm,保護管道直徑為450 mm,保冷層外徑與最外層管道之間充滿了空氣。利用有限容積法對管道截面處的溫度場分布進行了數值計算,得出了保冷層和保護管道的外表面溫度。在計算過程中做出如下模型假設:
(1) 輸送過程中管道介質溫度與管道內壁溫度相同;
(2) 輸送介質管道放置在保護管中,保冷層與保護管之間的空氣為對流換熱;
(3) 管道在運行過程中處于穩態導熱。

圖1 LNG管道截面溫度場分布等值線圖Fig.1 The contour map of LNG pipeline cross-section temperature distribution
通過計算得出管道保冷層厚度為95 mm時其表面溫度為-20.55 ℃,保護管道外表面溫度為 6.69℃。
3.1 管道的水力計算
水力計算由于粘度較小,運行流速較大,其流態一般為混合摩擦區,根據列賓宗公式可得其水力坡降為

沿程摩阻為

式中:i— 水力坡降;
Q— 體積流量,m3/s;
d— 管道內徑,m;
v— LNG的運動粘度;
L—管道輸送的長度,m。
3.2 管道的熱力計算
管道的熱力計算可有下式得出:

式中:G— 油品的質量流量,kg/s;
c— 平均溫度下輸送介質的比熱容,J/(kg·℃);
D— 管道外徑,m;
TR— 管道起點輸送介質溫度,℃;
TL— 距起點L處輸送介質溫度,℃;
T0— 周圍介質溫度,℃;
i— 水力坡降;
g— 重力加速度,m/s2。
3.3 總傳熱系數K

式中:Dw— 管道的結構外徑,m;
Di、Di+1— 鋼管、防腐絕緣層及絕熱層的內徑和外徑,m;
λi— 與上述各層相對應的導熱系數W/(m·℃);
K— 管道總傳熱系數,W/(m2·℃);
α1— 介質至管內壁的放熱系數,W/(m2·℃);
α2— 管外壁至大氣的放熱系數,W/(m2·℃);
D— 計算直徑,m。對于保冷管道,可取保冷層內外直徑平均值。 最終通過計算得出輸送到6 km的溫度為-153.39 ℃,壓降為 0.338 MPa。由于LNG管道輸送到終端的壓力降到0.856 MPa,相對應的飽和壓力為-128 ℃,而此時的溫度上升為-153.39℃,故由此得出在1.194 MPa下輸送,LNG仍處于液態,從而實現無氣化輸送。
(1) 通過管道敷設方式的分析得出了在人口較為稠密地區敷設 LNG管線較為理想的方式為管中管埋地敷設;
(2) 通過有限容積法計算得出95 mm的保冷層厚度結果較為理想,在該厚度條件下,能夠達到無氣化輸送;
(3) 通過實例分析,在1.194 MPa下輸送,由于輸送介質粘度較小,保溫層熱阻較大,致使整個管道溫升較小,其介質仍處于“過冷”狀態。
[1] 顧安忠,魯雪生,汪榮順,等.液化天然氣技術[M].北京:機械工業出版社,2003:18-20.
[2] K.S.佩德森.石油和天然氣的性質[M].北京:中國石化出版社,1989:204-207.
[3] B.E.波林,J.M.普勞斯尼茨,J.P.奧康奈爾,等.氣液物性估算手冊[M].北京:化學工業出版社,2006:157.
[4] 童景山.流體熱物性學基本理論與計算[M].北京:中國石化出版社,2008:155.
[5] 楊筱蘅.輸油管道設計與管理[M].東營:中國石油大學出版,2006:33-77.
Calculation of Conveying Process of LNG Pipeline
LI Qing-jie1,CHEN Bao-dong1,HONG Li-na2,MA Jin-peng3
(1. College of Petroleum Engineering,Liaoning Shihua University, Liaoning Fushun 113001, China; 2 .Liaohe Oil Field Company Jinzhou Oil Product Plant, Liaoning Jinzhou 121209, China; 3. PetroChina Pipeline Daqing Oil & Gas,Heilongjiang Daqing 163458,China)
Taking LNG terminal construction of a LNG company as background, operation parameters of a 6 km pipeline were calculated. Numerical calculation of the temperature field distribution of pipeline cross-section was carried out by the finite volume method ,and the outward surface temperature was gained when the insulation thickness was 95 mm. After selection of the insulation thickness and outer pipe diameter, the pressure drop and temperature rise at the terminal was calculated, the feasibility of “supercooling” state transportation was verified in the LNG transportation.
LNG; Tube in tube Conveying; Insulation thickness; Supercooling
TE 832
A
1671-0460(2011) 02-0163-02
遼寧省高等學校優秀人才支持計劃(2009R37)。
2010-07-16
李慶杰(1984-),男,碩士研究生,江蘇贛榆人,遼寧石油化工大學油氣儲運專業在讀,研究方向:從事 LNG管道輸送研究。E-mail:liqingjie01@yahoo.com.cn。