陳晶華， 付 璇

（1.中海石油（中國(guó)）有限公司北京研究中心，北京100027；

2.中國(guó)石油大學(xué)（北京）城市油氣輸配北京市重點(diǎn)實(shí)驗(yàn)室，北京102249）

近年來，海上油田不斷發(fā)展，對(duì)海底輸油管道停輸過程傳熱問題的研究也迫在眉睫。研究海底輸油管道停輸過程的傳熱問題，不僅可降低輸油成本、減少能耗、保護(hù)環(huán)境，而且可以為實(shí)際生產(chǎn)管理提供科學(xué)的依據(jù)。對(duì)于指導(dǎo)海上油田的輸油生產(chǎn)、管道安全運(yùn)行和節(jié)能降耗也具有重要意義。

管道運(yùn)行過程中，不可避免地會(huì)遇到計(jì)劃性或事故性停輸問題。而管內(nèi)原油溫度不斷下降，使油品黏度增大、形成結(jié)構(gòu)并膠凝［1］，管道再啟動(dòng)所需壓力顯著增加。如果停輸時(shí)間過長(zhǎng)，管道再啟動(dòng)所需壓力超過泵的揚(yáng)程或管道強(qiáng)度，則無法進(jìn)行正常的再啟動(dòng)過程。所以為了防止易凝高黏原油管道出現(xiàn)凝管現(xiàn)象，需要確定管道內(nèi)的溫度分布情況。

針對(duì)管道的停輸溫降，國(guó)內(nèi)外許多學(xué)者已經(jīng)做了大量研究工作，克利文推導(dǎo)出埋地管道周圍土壤溫度場(chǎng)的計(jì)算公式［1］，J.A.Wheeler［2］應(yīng)用有限元、有限差分理論研究了埋地管道溫度分布。我國(guó)羅塘湖［3］、李長(zhǎng)?。?］、邢曉凱等［5］也從數(shù)值模擬的角度對(duì)原油管道的溫降進(jìn)行了模擬計(jì)算。目前，對(duì)原油管道停輸降溫過程的研究主要采用數(shù)值計(jì)算的方法，通過控制網(wǎng)格的劃分，計(jì)算得到較高精度的解［6］。本文利用Fluent對(duì)原油管道停輸溫降過程進(jìn)行了模擬，分析研究原油管道在停輸后徑向溫度分布隨時(shí)間的變化規(guī)律。

1 管道概況

由海上平臺(tái)到陸上輸油站的原油管線全長(zhǎng)58.78km，包括平臺(tái)到海岸的海底管段和陸上埋地管段兩部分。其中海底管段長(zhǎng)44.2km，保溫層厚50mm，總傳熱系數(shù)為0.85，管道管徑為Φ406.4 mm×14.3mm；陸上管段長(zhǎng)14.58km，保溫層厚60mm，總傳熱系數(shù)為0.9，管道管徑為Φ406.4 mm×11.1mm。管道輸量在一定范圍內(nèi)變化，海底管段冬季平均環(huán)境溫度為3.7℃，夏季平均環(huán)境溫度21.7℃，陸管段環(huán)境溫度按極端氣溫考慮，最高為34.6℃，最低為－17.2℃。管道所在位置的冬夏溫差比較大，為保證管道安全運(yùn)行，需要確定管內(nèi)原油停輸一定時(shí)間后的管道內(nèi)的溫度分布，這對(duì)大口徑海底管道再啟動(dòng)壓力計(jì)算有很大影響。

2 原油及管材基本物性

原油20℃下的密度895.95kg／m3，凝點(diǎn)30℃，比熱容2 000J／（kg·K），初始油溫80℃，臨界溫度為35℃，原油含蠟質(zhì)量分?jǐn)?shù)為23.2%。原油的黏溫?cái)?shù)據(jù)如表1所示，管材性質(zhì)如表2所示。

表1 原油黏溫?cái)?shù)據(jù)Table1 The viscosity of the crude oil

表2 管材性質(zhì)Table2 The nature of the pipeline

3 管道停輸后管內(nèi)的傳熱過程及特點(diǎn)

熱油管道停輸后，由于管壁溫度始終低于管內(nèi)液態(tài)原油的溫度，因此形成由管中心指向管壁的徑向溫度梯度［7］。由于管內(nèi)的自然對(duì)流，管道截面上的原油溫度以一個(gè)較快的速度比較均勻地下降。當(dāng)油溫低于析蠟點(diǎn)時(shí)，原油中的蠟將逐漸析出，并和油相互作用形成凝油層。由于管壁處溫度最低，因而該處原油最先發(fā)生膠凝。隨著停輸時(shí)間的不斷延長(zhǎng)，管內(nèi)油溫不斷下降，凝油層也會(huì)不斷增厚，液相區(qū)面積不斷減?。?］。自然對(duì)流逐漸減弱，凝油區(qū)的熱傳導(dǎo)將占主要地位［9］。由于石蠟結(jié)晶釋放潛熱以及凝油層的保溫作用，會(huì)使得管道的停輸溫降變慢。根據(jù)停輸管道的徑向溫度分布和所輸原油滯流點(diǎn)可以確定自然對(duì)流區(qū)和導(dǎo)熱區(qū)的邊界，這是停輸管道溫降計(jì)算的關(guān)鍵。

此外，根據(jù)管道內(nèi)徑向溫度分布還可以判斷出凝油層界面位置和析蠟量等重要參數(shù)，對(duì)停輸管道水力分析具有重要意義［10］。凝油層界面的位置不能通過實(shí)驗(yàn)來測(cè)得，故采用數(shù)值模擬的方法進(jìn)行分析，能夠?yàn)閷?shí)際生產(chǎn)和理論研究提供依據(jù)。

4 停輸管道徑向溫度分布計(jì)算

在研究過程中，由于水下原油管道的傳熱方式、熱效應(yīng)范圍與埋地管道不同［11］，海水的熱容量較大，且與管道最外層之間的放熱系數(shù)很大，所以海底原油管道在溫降模擬過程中可假設(shè)環(huán)境溫度保持不變。由于原油密度變化和重力作用，管內(nèi)原油與管壁、保溫層與管外海水進(jìn)行熱交換的方式屬于自然對(duì)流換熱。由于停輸后原油密度的變化不大，可選擇采用Boussinesq模型解決自然對(duì)流問題。原油與管壁之間、管壁與保溫層之間均采用耦合邊界。

在管道上取一部分管段作為研究對(duì)象，建立三維立體模型并對(duì)整個(gè)圓柱區(qū)域進(jìn)行網(wǎng)格劃分，用Fluent流體分析軟件模擬管道停輸過程，并計(jì)算出所取管段的流場(chǎng)分布及溫度分布云圖。在計(jì)算時(shí)，假設(shè)管道內(nèi)部完全充滿原油，即任何截面的介質(zhì)都是均勻分布的。環(huán)境溫度冬季取3.7℃，初始油溫為80、40℃時(shí)，停輸時(shí)間分別取12、20h。夏季環(huán)境溫度取21.7℃，初始油溫分別為80、40℃時(shí)，停輸時(shí)間分別取12、20h。停輸開始時(shí)，可認(rèn)為管內(nèi)油溫分布均勻，再以此作為初始條件，開始非穩(wěn)態(tài)計(jì)算。計(jì)算結(jié)果如圖1所示。

由圖1可以看出，隨著停輸時(shí)間的延長(zhǎng)，初溫80℃時(shí)冬季管內(nèi)溫降比夏季快很多，且在相同的停輸時(shí)間下，冬季管中心處與管壁處的溫差比夏季大。根據(jù)不同時(shí)間段的斷面溫度分布，結(jié)合管輸原油的凝點(diǎn)、溶解度曲線（DSC）可以進(jìn)一步確定停輸管道的凝油層界面位置及原油析蠟量等參數(shù)。同時(shí)，通過初溫80℃與40℃的溫降對(duì)比可以看出，提高初始油溫可延長(zhǎng)油品的停輸時(shí)間，但需考慮管材物性和熱能動(dòng)力消耗等因素。

圖1 不同季節(jié)和初溫停輸12h和20h徑向溫度分布Fig.1 The radial temperature distribution of different seasons and the initial temperature with the pipeline 12 h and 20 h shutdown

為了更直觀的看出初始油溫對(duì)停輸后管內(nèi)溫降的影響，用Fluent軟件計(jì)算得到的不同初始油溫距管中心100mm處的油溫隨停輸時(shí)間變化的曲線，如圖2所示。

圖2 不同初始油溫距管中心100mm處溫度隨停輸時(shí)間變化曲線Fig.2 The pipeline temperature curve of different initial temperature at the position where is 100 mm from the center over the shutdown time

由圖2可知，在其他條件不變的情況下，較高的初始油溫可延長(zhǎng)油品的停輸時(shí)間。但完全依靠過度提高初始油溫的方法并不可取，不僅增加燃料費(fèi)用也對(duì)管材性質(zhì)要求比較高。在0～20h溫降速度很快，主要是因?yàn)樵撾A段管內(nèi)原油的自然對(duì)流較強(qiáng)烈。停輸20h后的一段時(shí)間內(nèi)溫降緩慢，降溫在5℃以內(nèi)，這是因?yàn)楣軆?nèi)原油接近臨界溫度，原油黏度增大及蠟晶析出，使得自然對(duì)流強(qiáng)度減弱。

海底原油管道停輸時(shí)，海水溫度對(duì)管道影響很大，所以計(jì)劃停輸盡量選在夏季海水溫度較高時(shí)。為了清楚的了解夏季和冬季海水溫度對(duì)管內(nèi)溫降的影響，所以給出所選管道冬夏兩季距管中心100 mm處的溫度隨停輸時(shí)間變化的曲線，如圖3所示。

圖3 不同季節(jié)距管中心100mm處溫度隨停輸時(shí)間變化曲線Fig.3 The pipeline temperature curve of different seasons at the position where is 100 mm from the center over the shutdown time

由圖3可以看出，在停輸32h時(shí)，距管中心100 mm處夏季溫降為22℃，冬季由于海水溫度低，距管中心100mm處溫降比夏季大7℃左右。

5 結(jié)論

通過Fluent模擬，結(jié)合管道運(yùn)行工況，預(yù)測(cè)了輸油管線不同停輸時(shí)間管道溫降的情況，得到主要結(jié)論：

（1）隨著管道停輸時(shí)間的延長(zhǎng)，中心液相區(qū)的流動(dòng)半徑在不斷減小。

（2）提高初始油溫可延長(zhǎng)管道的安全停輸時(shí)間，但需考慮管材物性和熱能動(dòng)力消耗等因素。

（3）管道停輸0～20h溫降速度很快，主要是因?yàn)樵撾A段管內(nèi)原油的自然對(duì)流較強(qiáng)烈。停輸20 h后的一段時(shí)間內(nèi)溫降緩慢，降溫在5℃以內(nèi)，這是因?yàn)楣軆?nèi)原油接近臨界溫度，原油黏度增大及蠟晶析出，使得自然對(duì)流強(qiáng)度減弱。

（4）通過與現(xiàn)場(chǎng)采集的數(shù)據(jù)進(jìn)行對(duì)比，可知利用Fluent軟件對(duì)海底原油管道的溫度場(chǎng)進(jìn)行數(shù)值計(jì)算有較高的準(zhǔn)確性。

［1］楊筱蘅.輸油管道設(shè)計(jì)與管理［M］.東營(yíng)：中國(guó)石油大學(xué)出版社，2006：239－241.

［2］Wheeler J A.Simulation of heat transfer from a warm pipeline Buried in permafrost［C］.AICHE 74th National Meetings，New Orleans，La：［s.n.］，1973：102－105.

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

［4］李長(zhǎng)俊，駱建武，陳玉寶.埋地?zé)嵊凸艿绬⑤敓崃?shù)值模擬［J］.油氣儲(chǔ)運(yùn)，2002，21（12）：16－19.Li Changjun，Luo Jianwu，Chen Yubao.Numerical simulation of heated crude oil pipeline at starting［J］.Oil ＆ Gas Storage and Transportation，2002，21（12）：16－19.

［5］邢曉凱，張國(guó)忠.埋地?zé)嵊凸艿勒＿\(yùn)行溫度場(chǎng)的確定［J］.油氣儲(chǔ)運(yùn)，1999，18（12）：28－30.Xing Xiaokai，Zhang Guozhong.Determination of temperature field of crude oil pipeline［J］.Oil ＆ Gas Storage and Transportation，1999，18（12）：28－30.

［6］張煜，朱紅鈞，陳小榆.熱油管道停輸溫降過程的模擬研究［J］.石油工業(yè)計(jì)算機(jī)應(yīng)用，2010：17（1）：38－40.Zhang Yu，Zhu Hongjun，Chen Xiaoyu.Process simulation of hot oil pipeline temperature drop after shutdown［J］.Computer Applications of Petroleum，2010：17（1）：38－40.

［7］楊晶，趙興民，蘆靜，等.熱油管道停輸溫降過程的數(shù)值模擬［J］.科學(xué)技術(shù)與工程，2011，11（16）：3797－3800.Yang Jing，Zhao Xingmin，Lu Jing，et al.Numerical simulation of the course of temperature drop in the hot oil pipeline［J］.Science Technology and Engineering，2011，11（16）：3797－3800.

［8］王敏，于遠(yuǎn)洋.埋地?zé)嵊凸艿劳］敎亟档腃FD 模擬［J］.科學(xué)技術(shù)與工程，2011，11（22）：5281－5285.Wang Min，Yu Yuanyang.Buried oil pipeline in the CFD simulation of shutdown temperature drop ［J］.Science Technology and Engineering，2011，11（22）：5281－5285.

［9］龍安厚，張帆.基于Fluent的海底輸油管道停輸溫降數(shù)值模擬［J］.科學(xué)技術(shù)與工程，2011，11（34）：8474－8477.Long Anhou，Zhang Fan.The temperature drop numerical simulation of the submarine oil pipeline based on Fluent［J］.Science Technology and Engineering，2011，11（34）：8474－8477.

［10］趙虎，王為民，王雷，等.超稠油埋地?zé)彷敼艿辣厥У臄?shù)值模擬［J］.當(dāng)代化工，2012，41（1）：85－87.Zhao Hu，Wang Weimin，Wang Lei，et al.Numberical simulation for the buried super heavy oil heating transportation pipeline of invalidation of heat insulator［J］.Contemporary Chemical Industry，2012，41（1）：85－87.

［11］韓佳燁，丁魯振，丁寧.埋地?zé)嵊凸艿罍亟导巴寥罍囟葓?chǎng)［J］.遼寧化工，2013，42（12）：1437－1440.Han Jiaye，Ding Luzhen，Ding Ning.Temperature drop and soil temperature field of buried oil pipelines［J］.Liaoning Chemical Industry，2013，42（12）：1437－1440.