基于ANSYS軟件的熱油管路保溫的數值模擬

曹先慧，馬貴陽，于麗麗, 高艷波，孟向楠

（遼寧石油化工大學 石油天然氣工程學院， 遼寧 撫順 113001）

基于ANSYS軟件的熱油管路保溫的數值模擬

曹先慧，馬貴陽，于麗麗, 高艷波，孟向楠

（遼寧石油化工大學 石油天然氣工程學院， 遼寧 撫順 113001）

應用ANSYS軟件對埋地熱油管道沿徑向溫度進行了數值模擬，得到了不同保溫層厚度時管道徑向溫度及熱流量的變化。分析了含蠟層的徑向溫度變化，隨著保溫層厚度的增加，保溫層內部溫降變化減小，管道向外傳遞的熱流密度逐漸減小，保溫效果更好。熱流量減小的速度隨著保溫層厚度的增加變得緩慢，模擬結果與編制計算機程序計算結果相吻合。

熱油管道；徑向溫降；熱流密度；數值模擬

對埋地輸油管道進行保溫，能有效地阻止管道內油品溫度向周圍土壤的擴散，減少管內油品的熱量損失，達到確保管道安全運行、節約能源、降低輸送成本的目的。對管道施加保溫時，除了要選擇良好的保溫材料，采用合理的保溫結構之外，保溫的效果還取決于保溫層的厚度。本文借助ANSYS軟件對不同保溫層厚度的埋地熱油管道進行了數值模擬，對埋地熱油管道保溫層的厚度對管道徑向溫度即熱流量的傳遞影響規律進行了研究。

1 數值模擬分析的基本原理[1]

本論文采用ANSYS作為數值模擬分析的工具，從熱分析模塊展開研究。熱分析模塊是基于能量守恒原理的熱平衡方程，用有限元法計算各節點的溫度，并導出其他熱物理參數。

熱分析主要包括穩態傳熱分析和瞬態傳熱分析。在穩態傳熱分析中任一節點的溫度不隨時間變化，它主要用于分析穩態的熱載荷對系統或部件的影響；瞬態傳熱過程指一個系統的加熱或冷卻過程，主要用于計算一個系統隨時間變化的溫度場以及其他熱參數。本論文采用瞬態熱分析進行分析計算。

瞬態傳熱過程指一個系統的加熱或冷卻過程。在這個過程中系統的溫度、熱流率、熱邊界條件及系統內能都隨時間明顯變化，根據能量守恒，瞬態熱分析的有限平衡方程為（以矩陣形式表示）:

式中： [K]—傳導矩陣，包含導熱系數、對流系數及輻射率和形狀系數；[C]—比熱矩陣，考慮系統內能的增加；{T}—節點溫度向量；—溫度對時間的導數；[Q]—節點熱流率向量，包括熱生成。

熱分析的控制方程：

其中，Vx，Vy，Vz為媒介傳導速率。

2 計算公式[2]

熱流量計算公式：

保溫層外壁溫度計算公式：

放熱系數a2：

其中：Q — 熱流量，W/m；

t0，t，tw— 分別為土壤溫度，管內油品溫度，保溫層外壁溫度，℃；

λ0，λla，λg，λb— 分別為土壤導熱系數，蠟層的導熱系數，鋼管的導熱系數，保溫層的導熱系數，W/(m·℃)；

a1— 油流至管內壁放熱系數，W/（m2·℃）；

a2— 管外壁至土壤的放熱系數，W/（m2·℃）；

rla,rn,rw,rb— 分別為蠟層內徑，鋼管內徑，鋼管外徑，保溫層外徑m；

h，Db— 分別為管道埋深和保溫層外直徑，m。

3 模擬計算與分析

以某埋地輸油管道為例，管道規格為D426× 7，管中心埋深H=1.8 m，平均地溫為-10 ℃，土壤導熱系λ0=1.5 W/(m·℃),油品溫度70 ℃，結蠟厚度為8 mm，蠟層導熱系數2.5 W/(m·℃)，鋼管導熱系數λg=46 W/(m·℃)，保溫材料選取聚氨酯硬質泡沫塑料，材質的參數為：密度ρ=50 kg/m3，導熱系數λb=0.035 W/(m·℃)，a1=556.432 W/（m2·℃），設置保溫層厚度分別為δ=10 ，20 ，30 ，40 ，50 ，60 ，80 ，100 mm時其對應的管外壁放熱系數a2值分別為2.454, 2,418, 2.384，2.352，2.320， 2.290，2.233， 2.18 W/（m2·℃）,油品流速0.5 m/s。

3.1 模型的建立

目前國內采用的管道保溫結構基本形式為[3]：鋼管—防腐層—保溫層—防水保護層，防腐層和防水保護層相對于保溫層的厚度對保溫效果的影響可以忽略不計，得到簡化后的埋地熱油管道的物理模型如圖1。

圖1埋地熱油管道的簡化模型Fig.1 simplified model of buried hot oil pipeline

3.2 數值模擬及結果分析

數值模擬及結果分析如圖2-7和表1。

圖2 保溫層沿徑向溫度分布曲線Fig.2 The radial temperature of the thermal insulation layer

圖3 保溫層管道三維溫度分布云圖Fig.3 The three dimensional distributed cloud chart of the pipeline

圖4 100 mm保溫層徑向溫度分布云圖Fig.4 The cloud chart of the radial temperature of the thermal insulation layer of 100 mm

圖5 管道不同位置處溫度隨保溫層厚度的變化Fig.5 The change of the temperature of different position of the pipeline with the insulation layer thickness

表1 不同保溫層厚度的熱流量Table 1 The heat flow of thermal insulation layer with different thickness

圖6不同保溫層厚度下的熱流密度Fig.6 The heat flow of thermal insulation layer with different thickness

埋地熱油管道的傳熱過程分三部分組成[4]，即管內熱油以對流方式將熱量傳給結蠟層內壁，經過結蠟層、管壁、保溫層將熱量傳給周圍土壤，再經地面與大氣換熱。通常在傳熱計算當中我們近似將油溫作為蠟層內壁溫度，由圖5可以看出，油溫為70 ℃時，保溫層厚度分別為10，20，30，40，50，60，80，100 mm時對應的結蠟層內壁溫度分別為69.770，69.834，69.868，69.890，69.906，69.917，69.932，69.942 ℃,其值非常接近油溫，但是有較小溫差存在，并且結蠟層內壁溫度隨著保溫層厚度的增加而升高。

圖7保溫層外壁溫度隨保溫層厚度的變化Fig.7 The change of the temperature of the outer wall of the thermal insulation layer with its thickness

對于計算結果分析可知，保溫層厚度不同，其散熱量也不同，通過用VB編制計算機程序計算得到的散熱量的理論解為Q1，由ANSYS軟件模擬得到的散熱量為Q2，數據比較見表1。

從表1中可以看出，Q1與Q2的誤差非常小，在允許范圍內，由此ANSYS軟件得到了驗證。并且保溫層向外傳遞熱流密度值隨著保溫層厚度的增大而明顯減少[5]。這是由于保溫層越厚，保溫層的熱阻增大，導熱系數越小，管道散失的熱量越少，對于一定管道長度的管線，增加保溫層厚度可以提高出口油品的溫度，節約能源。

將模擬計算數據繪成圖6可知：隨管道保溫厚度δ的不同，其管道的散熱量變化很大。當保溫層厚度δ較小時，散熱量Q值隨δ變化激烈；當保溫層厚度δ達到一定值后，散熱量Q值隨δ變化較緩慢。

取不同的保溫層半徑，進行運算后，得到不同保溫層厚度下的沿徑向溫降如圖8所示。

圖8 不同保溫層厚度的管道徑向溫降Fig.8 The drop of the radial temperature of thermal insulation layer with different thickness

取圖8中的0 mm至15 mm（結蠟層內壁為徑向距離原點）段作圖得到蠟層和鋼管內的溫度變化圖9。

從結蠟層內壁至鋼管外壁的范圍內，徑向溫度接近油溫70 ℃，這是由于保溫層厚度為8 mm比較小，并且鋼管的導熱系數較大，對管道的徑向溫度的影響小于1 ℃。在0 mm至8 mm即結蠟層內，保溫層越厚，溫度降低越慢。在8 mm至15 mm即鋼管壁內，曲線接近水平，溫度變化非常小，這是由鋼管較大的導熱系數所致。當徑向距離大于15 mm時，即從鋼管外壁至保溫層外壁，管道沿徑向的溫度是呈弧線均勻降低的。保溫層厚度為10 mm時的斜率最大，保溫層厚度為100 mm時的斜率最小，即保溫層厚度越大，管道徑向溫度降低的速度越慢。熱油管道運行過程中，在徑向溫降的推動下，油流所攜帶的熱量不斷往管外散失。保溫層厚度為10 mm的管道保溫層外壁溫度為36.206 ℃，保溫層厚度為100 mm的管道保溫層外壁溫度為-0.612℃。

圖9 蠟層和鋼管內的溫度變化圖Fig.9 The temperature variation of the wax layer and steel pipe

4 結束語

應用ANSYS軟件對不同保溫層下的埋地熱油管道沿徑向溫降的溫度分布進行數值模擬，結果與公式編制的VB程序求得的理論解相吻合，驗證了ANSYS軟件在模擬管道傳熱方面的可行性和準確性。

應用ANSYS模擬了不同保溫層厚度的埋地熱油管道沿向的熱流量與溫降變化，得到了保溫層厚度對管道傳熱的影響，為埋地熱油管道結構參數的合理設計提供了理論依據和分析方法。

［1］王澤鵬，張秀輝,等，ANSYS 12.0 熱力學有限元分析從入門到精通[M].2010:112-129.

［2］吳國忠，曲洪權，龐麗萍,等，埋地輸油管道非穩態熱力計算求解方法[J].油氣田地面工程，2001（6）：6-7．

［3］楊嘉瑜,譯.國外輸油管道保溫技術概況及發展[J].國外油田工程，1997（5）:42-45．

［4］賈雪松.基于CFD的熱油管道數值模擬[D].大慶石油學院,2010.

［5］王常斌，趙月,等：熱油管道溫度場的數值模擬[J].管道技術與設備，2011（3）:9-12.

Numerical Simulation on Thermal Insulation of Hot Oil Pipelines Based on ANSYS

CAO Xian-hui，MA Gui-yang，YU Li-li，GAO Yan-bo，MENG Xiang-nan
（Liaoning Shihua University, Liaoning Fushun 113001，China）

Through the numerical simulation of radial temperature of buried hot oil pipelines with ANSYS software, effect of the thermal insulation layer thickness on the radial temperature and heat flow were analyzed and calculated. Change of the radial temperature of different wax layers was analyzed. The results show that along with the increase of the insulation layer thickness, the reducing of temperature in the thermal insulation layer is slow and the heat flux that transfers between pipes and outside reduces gradually，so the heat preservation effect is better；With the increase of the insulation layer thickness, the decline rate of the heat flow is slow. The simulation result of numerical simulation is consistent with theory calculation result, which shows the feasibility and accuracy of ANSYS in heat transfer calculation.

Hot oil pipeline; Radial temperature drop; Heat flow; Numerical simulation

TE 832

A

1671-0460（2012）02-0181-04

2011-12-19

曹先慧（1986-），女，山東濟寧人，現為遼寧石油化工大學油氣儲運工程專業在讀碩士。研究方向：計算流體力學數值仿真。E-m ail：caoxianhui@yeah.net。