雷啟盟，賈海洋，付曉明

?

熱油管道停輸再啟動數(shù)值模擬

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

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

通過建立管道停輸再啟動非穩(wěn)態(tài)數(shù)學(xué)模型，采用數(shù)值模擬的方法進行求解，考察不同輸量、輸油溫度下停輸再啟動過程管道起點壓力隨時間的變化規(guī)律。研究表明，一定條件下，管道起點壓力峰值隨管道輸量的增加而增大，隨輸油溫度的增加而減?。磺译S管道輸量的增加，再啟動過程壓力峰值恢復(fù)到穩(wěn)定狀態(tài)所需時間減小。

熱油管道；數(shù)值模擬；停輸再啟動

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

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

1 數(shù)學(xué)模型與計算方法

1.1 正常運行的數(shù)學(xué)模型

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

圖1 埋地?zé)嵊凸艿朗疽鈭D

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

連續(xù)性方程：

動量方程：

能量方程：

——貫流斷面面積，m2；

——原油比內(nèi)能，J/kg；

——高程，m；

——原油比焓，J/kg；

——時間，s；

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

——管道軸向位置，m；

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

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

——達西摩阻系數(shù)；

——管道內(nèi)徑，m；

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

根據(jù)連續(xù)性方程、動量方程和能量方程，推導(dǎo)出管道正常運行數(shù)值模型：

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

——油溫，℃；

——原油膨脹系數(shù)，℃-1。

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

——導(dǎo)熱系數(shù)，W/(m·℃)；

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

（3）土壤導(dǎo)熱方程：

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

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

T——土壤溫度，℃；

——土壤導(dǎo)熱系數(shù)，W/(m·℃)；

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

——深度，m。

（4）管內(nèi)原油、結(jié)蠟層、管壁、防腐層以及土壤的傳熱過程相互關(guān)聯(lián)，滿足：

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

0——管道內(nèi)壁溫度，℃。

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

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

T——恒溫層溫度，℃；

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

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

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

1.2 停輸過程的數(shù)學(xué)模型

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

在方程（9）中，對于滯流狀態(tài)，導(dǎo)熱系數(shù)取真實的導(dǎo)熱系數(shù)；對于自然對流狀態(tài)，取當(dāng)量導(dǎo)熱系數(shù)。結(jié)蠟層、管壁和防腐層的導(dǎo)熱方程和土壤導(dǎo)熱方程與式（5）和式（6）相同。

1.3 再啟動過程的數(shù)學(xué)模型

當(dāng)原油沒表現(xiàn)出觸變性時，再啟動過程的數(shù)值模型與正常輸送相同。當(dāng)原油表現(xiàn)出觸變性時，再啟動過程的數(shù)值模型為：

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

1.4 計算區(qū)域離散化及數(shù)值計算方法

為保證計算精度和準(zhǔn)確性，采取DELAUNAY三角化方法、極坐標(biāo)網(wǎng)格和有限差分法在管道熱力影響區(qū)內(nèi)對土壤、結(jié)蠟層、鋼管壁、防腐層和管道分別進行離散，結(jié)果見圖2和圖3。

圖2 結(jié)蠟層、管壁和防腐層及其周圍土壤的網(wǎng)格劃分

圖3 管道離散示意圖

2 分析與討論

2.1 模擬工況

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

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

2.2 停輸再啟動壓力變化規(guī)律

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

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

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

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

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

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

3 結(jié)束語

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

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

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

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

［3］李長俊，駱建武, 等. 埋地?zé)嵊凸艿绬⑤敓崃δM數(shù)值模擬[J].油氣儲運，2002, 21 (12)：16-19.

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

［5］崔秀國, 張勁軍. 埋地?zé)岷炘凸艿赖姆欠€(wěn)態(tài)傳熱問題[J]. 中國工程科學(xué), 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-），男，研究生，黑龍江省佳木斯市人，研究方向：油氣儲運工程。

研究生科技創(chuàng)新項目，項目號：xm173411。

TQ 794.2

A

1004-0935（2017）12-1182-04