熱油管道停輸再啟動數值模擬

雷啟盟，賈海洋，付曉明

?

熱油管道停輸再啟動數值模擬

雷啟盟1，賈海洋2，付曉明3

（1. 東北石油大學 石油工程學院， 黑龍江 大慶 163318； 2. 大慶油田工程建設有限公司油建公司物資管理中心，黑龍江 大慶 163453； 3. 大慶油田有限責任公司(儲運銷售分公司)，黑龍江 大慶 163517）

通過建立管道停輸再啟動非穩態數學模型，采用數值模擬的方法進行求解，考察不同輸量、輸油溫度下停輸再啟動過程管道起點壓力隨時間的變化規律。研究表明，一定條件下，管道起點壓力峰值隨管道輸量的增加而增大，隨輸油溫度的增加而減小；且隨管道輸量的增加，再啟動過程壓力峰值恢復到穩定狀態所需時間減小。

熱油管道；數值模擬；停輸再啟動

我國盛產易凝高粘原油[1]，在常溫下流動性較差，主要采用加熱工藝降低原油粘度的方式輸送。隨含蠟原油逐漸進入事故高發期，不可避免的導致管道停輸，如在停輸時間較長的情況下，極易發生凝管事故，嚴重影響管道的安全運行。目前，針對熱油管道停輸再啟動過程已有一定的研究報道。羅塘湖[2]、李長俊[3]等從解析計算和數值模擬的角度模擬計算了熱油管道的溫降，得到了一定條件下管內原油的溫度分布。安家榮、史秀敏等[4]通過建立管道停輸再啟動過程數學模型，開發熱油管道停輸與再啟動過程的數學模擬軟件，解決了熱油管道停輸與再啟動過程預測的技術難題。

目前，針對熱油管道停輸再啟動過程問題的研究主要采用數值計算的方法[5]，采用有限元法和有限體積法，通過控制網格的劃分，計算得到較高精度的解。因此，本文針對熱油管道停輸再啟動過程進行數值模擬，分析不同輸量、不同輸油溫度管道起點壓力的變化規律，為熱油管道的生產管理提供一定的科學依據。

1 數學模型與計算方法

1.1 正常運行的數學模型

由于埋地管道真實的熱力變化情況比較復雜，在建立管道正常運行數值計算模型時，假設管道周圍土壤為各向同性的均勻介質，且管道在土壤內的熱力影響區域取為方向5 m，土壤恒溫層為方向12.5 m。

圖1 埋地熱油管道示意圖

（1）基于上述簡化與假設, 建立熱油管道正常運行的數學模型。

連續性方程：

動量方程：

能量方程：

——貫流斷面面積，m2；

——原油比內能，J/kg；

——高程，m；

——原油比焓，J/kg；

——時間，s；

——原油平均速率，m/s；

——管道軸向位置，m；

——管道軸向與水平方向的夾角，°；

——油流截面平均壓力，Pa；

——達西摩阻系數；

——管道內徑，m；

——原油在單位時間內單位管壁面積上的散熱量，W/m2。

根據連續性方程、動量方程和能量方程，推導出管道正常運行數值模型：

式中：c——原油定壓比熱容，J/(kg·℃)；

——油溫，℃；

——原油膨脹系數，℃-1。

（2）結蠟層、管壁和防腐層的傳熱方程：

——導熱系數，W/(m·℃)；

——取1、2、3，分別代表結蠟層、管壁和防腐層。

（3）土壤導熱方程：

式中：——土壤密度，kg/m3；

C——土壤比熱容，J/(kg·℃)；

T——土壤溫度，℃；

——土壤導熱系數，W/(m·℃)；

——垂直于軸向的水平位置，m；

——深度，m。

（4）管內原油、結蠟層、管壁、防腐層以及土壤的傳熱過程相互關聯，滿足：

式中：0——油流對管內壁的放熱系數，W/(m·℃)；

0——管道內壁溫度，℃。

（5）由于計算區域的對稱性，僅取管道的右半部分進行研究，邊界條件為：

式中：T——大氣溫度，℃；

T——恒溫層溫度，℃；

——地表向大氣的放熱系數，W/(m·℃)；

——熱力影響區域水平方向邊界位置，m；

——熱力影響區域深度方向邊界位置，m。

1.2 停輸過程的數學模型

在建立停輸溫降過程數值模型時，將液態原油的自然對流問題當成導熱問題來處理，引入當量導熱系數；同時，假設管道和蠟晶析出形成凝油層的增長同心，將停輸溫降問題簡化為一個純導熱問題。停輸溫降過程管內原油的導熱方程為：

在方程（9）中，對于滯流狀態，導熱系數取真實的導熱系數；對于自然對流狀態，取當量導熱系數。結蠟層、管壁和防腐層的導熱方程和土壤導熱方程與式（5）和式（6）相同。

1.3 再啟動過程的數學模型

當原油沒表現出觸變性時，再啟動過程的數值模型與正常輸送相同。當原油表現出觸變性時，再啟動過程的數值模型為：

式中：用Houska觸變模型來計算[6]。

1.4 計算區域離散化及數值計算方法

為保證計算精度和準確性，采取DELAUNAY三角化方法、極坐標網格和有限差分法在管道熱力影響區內對土壤、結蠟層、鋼管壁、防腐層和管道分別進行離散，結果見圖2和圖3。

圖2 結蠟層、管壁和防腐層及其周圍土壤的網格劃分

圖3 管道離散示意圖

2 分析與討論

2.1 模擬工況

管道全長30.44 km，管道外徑630 mm，壁厚為8 mm，根據管道實際運行數據的統計結果，選取不同輸油溫度和輸量，計算管道停輸再啟動管道起點壓力隨時間的變化。在給定環境溫度下，分別選取最高、最低和平均的輸量和輸油溫度，對不同工況下的停輸再啟動過程進行數值模擬，具體模擬工況見表1所示。

表1 管道停輸再啟動模擬工況

2.2 停輸再啟動壓力變化規律

不同輸量下，熱油管道停輸再啟動過程管道起點壓力隨時間的變化如圖4所示。

圖4 不同輸量的停輸再啟動起點壓力隨時間變化

由圖4可以看出，其他條件相同時，穩定運行時隨管道輸量的增加，管道起點壓力逐漸增加，是因為以恒流量的方式輸送原油時，原油通過外輸泵源源不斷地打入管內，輸量越大，原油沿管道流動過程的摩擦阻力越大，需要的起點壓力就越高。再啟動過程中，隨管道輸量的增加，管道起點峰值逐漸增加，輸量分別為40、80、120 kg/s時的起點壓力峰值為964.72、1 322.03、1 714.74 kPa，這是因為相同環境溫度和輸油溫度下，輸量越大，增加了原油與管道內壁之間的作用力，要想使管道在短時間內快速啟動，故起點的出站壓力迅速升高。啟動一定時間后，膠凝原油從屈服到穩定運行，隨著管道輸量增加，起點壓力由峰值到穩定值所需時間逐漸減小，這是由于輸量越大，管道內的剪切應力越大，由于原油的剪切稀釋性，原油粘度降低，導致壓力下降速率快，故達到穩定運行所需時間少。

不同輸油溫度下，熱油管道停輸再啟動過程管道起點壓力隨時間的變化如圖5所示。

圖5 不同輸油溫度的停輸再啟動起點壓力隨時間變化

由圖5可以看出，其他條件相同時，管道穩定運行時，隨輸油溫度的增加，管道起點壓力逐漸減小，是因為輸油溫度越高，原油粘度越低，管道的壓降減小。再啟動過程中，隨管道輸油溫度的增加，管道起點峰值逐漸減小，輸油溫度分別為40、50、60 ℃時的起點壓力峰值分別為1714.74、751.7、493.04 kPa。這是因為輸送油溫較高時，原油的膠凝結構弱，只需提供的壓頭可以克服管道的高程差，管內原油就可以恢復流動；隨著油溫的降低，原油將呈現出越來越強的屈服性，在施壓啟動過程中原油將經歷屈服、裂降過程，在這種情況下，當施壓而產生的剪切應力大于原油的屈服應力時，原油才會流動。此外，當油溫從40 ℃升高至50 ℃時，原油粘度大幅度降低，再啟動壓力迅速下降，而當油溫從50 ℃升高至60 ℃時，管道的再啟動壓力則變化不大，導致升高油溫對降低再啟動壓力的作用很有限。

3 結束語

（1）其他條件一定時，管道輸量越大，管道穩定運行時起點壓力越大，再啟動時壓力峰值越大，由峰值恢復到穩定值所需時間越小，起點壓力達到峰值后的壓降速率越大。

（2）其他條件一定時，管道輸油溫度越高，管道穩定運行時起點壓力越小，再啟動時壓力峰值越小，且油溫從40 ℃升至50 ℃時，再啟動壓力下降顯著，油溫從50 ℃升高至60 ℃時，再啟動壓力變化不大。

［1］楊筱衡，張國忠． 輸油管道設計與管理［M］．東營: 石油大學出版社，1996.

［2］羅塘湖. 含蠟原油流變特性及其管道輸送［M］. 北京：石油工業出版社，1991: 93-95.

［3］李長俊，駱建武, 等. 埋地熱油管道啟輸熱力模擬數值模擬[J].油氣儲運，2002, 21 (12)：16-19.

［4］安家榮, 史秀敏, 張國忠. 熱油管道停輸與再啟動過程模擬計算軟件[J]. 油氣儲運, 1998, 17 (3): 12-14.

［5］崔秀國, 張勁軍. 埋地熱含蠟原油管道的非穩態傳熱問題[J]. 中國工程科學, 2003 (07): 77-83.

［6］Houska M.Engineering aspects of the rheology of thixotropic liquids[D]. Prague: Czech Technical University of Prague, 1981.

Numerical Simulation of Shutdown and Restarting of Hot Oil Pipeline

1,2,3

（1. College of Petroleum Engineering, Northeast Petroleum University, Heilongjiang Daqing 163318,China ; 2. Oilfield Construction Company of Daqing Oilfield Construction Engineering Co., Ltd., Heilongjiang Daqing 163453,China; 3.Daqing Oilfield Limited Company(Storage,Transportation and Marketing Company), Heilongjiang Daqing 163517,China）

The unsteady mathematical model of pipeline shutdown and restartingwas established and solved by numerical simulation method. The variation law of the starting pressure of pipeline with time was studied in the process of shutdown and restarting with different throughput and different oil temperature. The results showed that, under certain conditions, the peak value of the starting pressure of pipeline increased with the increase of pipeline throughput and decreased with the increase of the oil temperature; and with the increase of the pipeline throughput, the time required for the pressure peak to resume to the steady state during the restart process decreased.

waxy crude oil pipeline; numerical simulation; restarting process

2017-11-20

雷啟盟（1993-），男，研究生，黑龍江省佳木斯市人，研究方向：油氣儲運工程。

研究生科技創新項目，項目號：xm173411。

TQ 794.2

A

1004-0935（2017）12-1182-04