青樹(shù)勇代 卿郭鴻雁杜明俊崔 奇揚(yáng)子寧

(1.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司北京分公司,北京100000;2.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司西南分公司,四川成都610000;3.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司華北分公司,河北任丘062552)

冷熱原油長(zhǎng)距離順序輸送動(dòng)態(tài)傳熱特性數(shù)值研究

青樹(shù)勇1,代 卿2,郭鴻雁1,杜明俊3,崔 奇3,揚(yáng)子寧3

(1.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司北京分公司,北京100000;2.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司西南分公司,四川成都610000;3.中國(guó)石油集團(tuán)工程設(shè)計(jì)有限責(zé)任公司華北分公司,河北任丘062552)

冷熱原油順序輸送對(duì)土壤溫度場(chǎng)的要求及其嚴(yán)格,這也是確保冷油過(guò)后,熱油能否安全進(jìn)站的關(guān)鍵?；趥鳠釋W(xué)和流體動(dòng)力學(xué),建立了埋地管道流動(dòng)與傳熱控制方程,數(shù)值模擬了冷熱原油順序輸送過(guò)程中管道沿線(xiàn)不同位置軸向油溫及土壤溫度場(chǎng)的動(dòng)態(tài)變化過(guò)程。研究表明,隨著出站油溫的冷熱交替周期運(yùn)行,管道沿線(xiàn)不同位置的管內(nèi)油溫及周?chē)欢ǚ秶鷥?nèi)的土壤溫度場(chǎng)呈現(xiàn)周期性變化,并存在一定的時(shí)間或空間滯后性,且對(duì)于低輸量運(yùn)行的管道來(lái)說(shuō),超過(guò)一定輸送距離后,輸送溫度趨于一致;在熱油-冷油交替輸送過(guò)程中,熱油受前端冷油的影響,熱油頭進(jìn)站溫度最低,這是管道安全運(yùn)行方案應(yīng)考慮的主要因素。

長(zhǎng)輸管道; 冷熱原油; 順序輸送; 動(dòng)態(tài)傳熱

隨著國(guó)際石油貿(mào)易的快速發(fā)展,國(guó)內(nèi)進(jìn)口原油總量不斷提高,煉化企業(yè)的油源也形成了多元化供給模式。與國(guó)產(chǎn)原油不同,進(jìn)口原油多為低黏度、低凝點(diǎn)、流動(dòng)性好的油品,與國(guó)產(chǎn)原油摻混加熱輸送顯然是不經(jīng)濟(jì)的[1],且不同性質(zhì)的原油直接摻混容易造成管內(nèi)油品分層,析出瀝青質(zhì),造成局部沉積,直接影響下游儲(chǔ)煉設(shè)施的安全[2]。冷熱原油順序輸送將成為今后原油管道的主要運(yùn)行方式。因此,對(duì)于工程設(shè)計(jì)單位來(lái)說(shuō),開(kāi)展冷熱原油順序輸送工藝的技術(shù)研究至關(guān)重要。

與成品油管道順序輸送不同,冷熱原油順序輸送對(duì)土壤溫度場(chǎng)的要求及其嚴(yán)格,這也是確保冷油過(guò)后,熱油能否安全進(jìn)站的關(guān)鍵。目前,國(guó)內(nèi)外已有部分學(xué)者對(duì)冷熱原油順序輸送技術(shù)進(jìn)行了較為深入的研究,并取得一定成果[3-5]。崔秀國(guó)等[6]建立了冷熱原油交替輸送過(guò)程非穩(wěn)態(tài)水力、熱力數(shù)學(xué)模型。李傳憲等[7]建立了冷熱原油順序輸送水力準(zhǔn)穩(wěn)態(tài)、熱力非穩(wěn)態(tài)數(shù)學(xué)模型,數(shù)值模擬了冷熱原油交替輸送過(guò)程中土壤溫度場(chǎng)的變化規(guī)律。王凱等[8]在一定假設(shè)的前提下建立了順序輸送非穩(wěn)態(tài)傳熱數(shù)學(xué)模型,并采用熱力特征線(xiàn)法和有限容積法相結(jié)合的數(shù)值算法對(duì)模型進(jìn)行了求解。筆者從工程設(shè)計(jì)角度出發(fā),基于有限差分和有限單元法,研究了冷熱原油順序輸送過(guò)程的動(dòng)態(tài)傳熱特性,所得結(jié)果可為實(shí)際工程設(shè)計(jì)提供一定的理論指導(dǎo)。

1 數(shù)學(xué)模型

基于一維流動(dòng)理論,假設(shè)管內(nèi)油溫在同一截面上均勻分布,僅為時(shí)間和管道軸向距離的函數(shù)。由于冷熱原油的黏度差異較大,故二者交匯處假設(shè)為活塞式驅(qū)油(即不考慮混油段),通過(guò)求解連續(xù)性方程、動(dòng)量方程、能量方程及流體面積方程可得出管道瞬態(tài)水力、熱力參數(shù)。

建立的控制方程:

式中,A為流體面積,m2;z為管道軸向距離,m;v為軸向Z方向的流體速度,m/s;ρo為流體密度, kg/m3;t為管道運(yùn)行時(shí)間,s;gc為量綱常數(shù);p為管道內(nèi)流體的壓力,Pa;f為摩阻系數(shù);h為熱膨脹系數(shù),1/℃;d為管道內(nèi)徑,m;U為流體的內(nèi)能,J;h為流體與管內(nèi)壁的對(duì)流換熱系數(shù),W/(m2·℃);T為流體溫度,℃;Tk為管內(nèi)壁溫度,℃;A0為管道橫截面積,m2;Am為面積系數(shù);pm為壓力系數(shù);p0為密閉輸送基準(zhǔn)壓力,Pa;T0為密閉輸送基準(zhǔn)溫度,℃。

忽略土壤的軸向傳熱,只考慮徑向傳熱。埋地?zé)嵊凸艿赖膹较騻鳠嵊扇糠纸M成,即管內(nèi)油流以對(duì)流方式將熱量傳給管道內(nèi)壁,經(jīng)管壁、保溫層將熱量傳遞給周?chē)寥?再經(jīng)地表與大氣進(jìn)行對(duì)流換熱[6-8]。這里忽略水分遷移和冰水相變(計(jì)算時(shí)采用實(shí)測(cè)的含水土壤導(dǎo)熱系數(shù))對(duì)土壤溫度場(chǎng)的影響,建立的管道徑向非穩(wěn)態(tài)傳熱數(shù)學(xué)模型如下:

管內(nèi)壁對(duì)流換熱方程:

式中,T為流體溫度,℃;λ1為油品的導(dǎo)熱系數(shù),W/ (m·℃);h為流體與管內(nèi)壁的對(duì)流換熱系數(shù),W/ (m2·℃);Tk為管內(nèi)壁溫度,℃

管道壁、保溫層傳熱方程:

式中,i=1、2分別表示鋼管及保溫層;ρi為第i層的密度,kg/m3;ci為第i層的比熱容J/(kg·℃);Ti為各層的溫度,℃;λi為第i層的導(dǎo)熱系數(shù),W/ (m·℃);r為徑向位置,m;θ為環(huán)向弧度。

土壤導(dǎo)熱方程:

式中,ρc為土壤的密度,kg/m3;cc為土壤的比熱容, J/(kg·℃);Tc為土壤的溫度,℃;λc為土壤的導(dǎo)熱系數(shù),W/(m·℃)。

2 數(shù)值計(jì)算及結(jié)果分析

以東北某輸油管道為例,管線(xiàn)全長(zhǎng)150 km,管徑D406 mm×7 mm,設(shè)計(jì)壓力6.3 MPa,管道中心埋深1.6 m,保溫層厚30 mm,輸送介質(zhì)為俄羅斯原油和國(guó)產(chǎn)原油,輸送方式為冷熱順序輸送,冬季最冷月管道沿線(xiàn)埋深處平均地溫4℃,地表環(huán)境平均溫度-15℃,平均風(fēng)速2 m/s。油品、土壤及保溫層的物性參數(shù)見(jiàn)表1,油品50 s-1剪切速率下的黏溫曲線(xiàn)見(jiàn)圖1。

表1 油品、土壤及保溫層物性參數(shù)表Table 1 Oil product,soil and insulation physical property parameter table

圖1 冷熱原油的黏溫特性曲線(xiàn)Fig.1 Viscosity temperature characteristic curve of cold&hot crude oil

目前,管道處于低輸量運(yùn)行狀態(tài),輸送流量300 m3/h,熱油出站溫度55℃,冷油出站溫度20℃,熱油頭進(jìn)站油溫不低于18℃,順序輸送周期為8 d(即冷熱原油各輸送4 d)。

2.1 冷熱原油順序輸送管道軸向溫降變化

圖2(a)給出了熱油過(guò)后,輸送冷油不同時(shí)刻管內(nèi)軸向油溫變化規(guī)律。由圖2(a)可知,隨著冷油的不斷注入,冷油頭受前端熱油和周?chē)寥罍囟葓?chǎng)的作用很快被加熱。根據(jù)管輸量核算,輸送冷油12 h后,冷油頭應(yīng)到達(dá)管道29.8 km處,從軸向油溫來(lái)看,此時(shí)冷油溫度與輸送熱油時(shí)的溫度基本相同,冷油被加熱的現(xiàn)象十分明顯,這主要是由于冷油出站時(shí),前端的熱油溫度(55℃)較高,冷熱原油溫差大,傳熱快,加之周?chē)寥佬罘e的熱量多,進(jìn)一步加快了油品的傳熱,使冷油前緣很快升高了一定的溫度。冷油輸送初期由于冷熱原油接觸時(shí)間短,雖然傳熱速率快,但傳遞熱量有限。因此,輸送初期管道軸向油溫過(guò)度段變化較為劇烈。隨著后續(xù)冷油的不斷注入,冷油前緣繼續(xù)吸收熱油和周?chē)寥赖臒崃?并向冷油尾端傳遞,被加熱的冷油段不斷延長(zhǎng),使得管道軸向油溫過(guò)度段變化越來(lái)越平緩。管道輸送60.3 h后,冷油頭到達(dá)末站,進(jìn)站油溫39.9℃;此時(shí)整個(gè)管道內(nèi)超過(guò)20℃的原油長(zhǎng)度約62.5 km,進(jìn)一步表明了熱油及管道周?chē)寥罍囟葓?chǎng)對(duì)冷油動(dòng)態(tài)傳熱特性的影響。

圖2 輸送冷/熱油時(shí)的軸向油溫變化Fig.2 The axial temperature variation of hot or cold crude oil transportation

冷油輸送96 h后,開(kāi)始輸送熱油,輸送熱油時(shí)不同時(shí)刻管內(nèi)軸向油溫變化規(guī)律見(jiàn)圖2(b)。由圖2 (b)可以看出,管內(nèi)軸向輸送熱油與輸送冷油的油溫變化正好相反。輸送熱油時(shí),熱油頭不斷被前端的冷油尾冷卻,加之熱油需要不斷的向土壤散熱,故熱油前緣溫降較快,輸送12 h后,熱油頭到達(dá)管道29.8 km處,油溫降至20.4℃。隨著熱油注入時(shí)間的延長(zhǎng),熱油前緣不斷向土壤散熱并被冷油尾冷卻降溫。60.3 h后熱油頭到達(dá)末站,進(jìn)站油溫18.2℃,即冷油尾溫度18.2℃。若該管道只輸送熱油,則進(jìn)站油溫39.8℃,與順序輸送熱油進(jìn)站油溫相差21.6℃,這部分熱損失主要用于加熱冷油尾及向土壤散熱。結(jié)合圖2可知,在冷油-熱油順序輸送過(guò)程中,冷油運(yùn)行96 h后,溫降只有1.8℃,進(jìn)一步說(shuō)明了冷油對(duì)熱油輸送的影響。

圖3給出了管道沿線(xiàn)不同位置冷熱原油順序輸送過(guò)程中管道截面油溫變化規(guī)律。由圖3可知,隨著出站冷熱原油的周期交替輸送,管道沿線(xiàn)不同位置的油溫呈現(xiàn)周期性變化,但是振幅不同,并表現(xiàn)出一定的時(shí)間滯后性。這主要是受空間距離的影響,管內(nèi)冷油頭或熱油頭向前運(yùn)行需要一定的時(shí)間。由圖3不難看出,距出站距離越遠(yuǎn),管內(nèi)油溫變化幅度相對(duì)越小,這主要是由于油流從出站口開(kāi)始向前運(yùn)行過(guò)程中不斷與前端油品及周?chē)寥肋M(jìn)行熱交換,熱交換的程度越來(lái)越大,強(qiáng)度相對(duì)減弱。仔細(xì)觀(guān)察圖3還會(huì)發(fā)現(xiàn),順序輸送過(guò)程中,沿線(xiàn)不同位置高溫油品流經(jīng)管道截面的時(shí)間要長(zhǎng)于低溫油品,即冷油段被加熱的長(zhǎng)度大于熱油段被冷卻的長(zhǎng)度。這主要是由于冷油-熱油順序輸送過(guò)程中,前端熱油及土壤傳遞給冷油的熱量多,使冷油前緣溫度升高較快,形成較長(zhǎng)的高溫區(qū)。

圖3 順序輸送管道不同位置的油溫變化Fig.3 The oil temperature change in different positions of pipeline in sequential transportation

2.2 管道周?chē)寥罍囟葓?chǎng)的變化

圖4、5分別給出了管道順序輸送4個(gè)周期后的首末站土壤溫度場(chǎng)云圖。結(jié)合計(jì)算的數(shù)據(jù)分析可知,管道周?chē)欢ǚ秶鷥?nèi)的等溫線(xiàn)隨著出站油溫的周期變化,呈現(xiàn)周期性移動(dòng),一定范圍以外的等溫線(xiàn)受輸送周期及土壤傳熱滯后性的影響,等溫線(xiàn)波動(dòng)不明顯。

圖4 輸送熱油后的首末站土壤溫度場(chǎng)Fig.4 The first and last station soil temperature field after the hot oil transportation

圖5 輸送冷油后的首末站土壤溫度場(chǎng)Fig.5 The first and last station soil temperature field after the cold oil transportation

輸送冷油前,上一批次熱油形成了較高的土壤溫度場(chǎng),冷油出站后油流溫度小于土壤溫度,土壤向油流放熱,隨著輸送時(shí)間的延長(zhǎng),土壤溫度場(chǎng)逐漸降低,高溫等溫線(xiàn)不斷向管中心方向移動(dòng)。

當(dāng)輸送熱油時(shí),上一批次的冷油形成低溫溫度場(chǎng),熱油出站后油流溫度大于土壤溫度,油流向土壤放熱,隨著輸送時(shí)間的延長(zhǎng),土壤溫度場(chǎng)重新建立,高溫等溫線(xiàn)不斷向遠(yuǎn)離管中心方向移動(dòng)。綜上所述:隨著出站油溫的循環(huán)周期變化,沿線(xiàn)不同位置的管內(nèi)油溫及周?chē)欢ǚ秶鷥?nèi)的土壤溫度場(chǎng)也呈現(xiàn)周期性循環(huán)變化。

圖6給出了管道首末站正上方0.5 m處,土壤溫度隨冷熱原油順序輸送過(guò)程的變化規(guī)律。由圖6可知,隨著冷熱原油周期順序輸送的進(jìn)行,首末站管道正上方0.5 m處的土壤溫度呈周期性變化。但受保溫層及土壤傳熱滯后性的影響,冷熱原油輸送過(guò)程中該處的土壤溫度波動(dòng)不大,且隨著循環(huán)周期的增加,該處的溫度波動(dòng)幅度進(jìn)一步減小,并逐漸趨于平穩(wěn)。循環(huán)4個(gè)周期后,首站管道正上方0.5 m處的最大波動(dòng)溫差只有1.4℃。末站管道正上方0.5 m處的最大波動(dòng)溫差只有1.1℃,這也進(jìn)一步說(shuō)明了保溫層對(duì)熱油管道傳熱特性的影響是非常明顯的。

圖6 管道正上方0.5 m處的土壤溫度變化Fig.6 The soil temperature changes above the pipeline 0.5 m

3 結(jié)論及建議

(1)當(dāng)管道低輸量運(yùn)行時(shí),超過(guò)一定輸送距離后,無(wú)論是冷油-熱油還是熱油-冷油順序輸送,原油溫度都趨于一致,這主要是由于隨著輸送距離和時(shí)間的延長(zhǎng),熱交換程度充分的緣故。(2)在分析了各類(lèi)制約因素后,發(fā)現(xiàn)熱油-冷油交替輸送過(guò)程中,受前端冷油的影響,熱油頭進(jìn)站溫度最低,這是管道安全經(jīng)濟(jì)運(yùn)行方案應(yīng)考慮的主要因素。合理確定冷油的出站溫度對(duì)整個(gè)管道安全平穩(wěn)運(yùn)行及其重要。(3)從計(jì)算的土壤溫度場(chǎng)來(lái)看,隨著出站油溫的交替運(yùn)行,沿線(xiàn)管道不同位置周?chē)寥赖牡葴鼐€(xiàn)呈周期性移動(dòng),且振幅不同,并存在空間滯后性,這主要是受不同距離處管內(nèi)油溫與土壤換熱的影響。(4)目前,關(guān)于冷熱原油順序停輸再啟動(dòng)的研究還不多,且考慮的多為停輸過(guò)程的溫降變化,對(duì)其啟動(dòng)壓力及啟動(dòng)過(guò)程對(duì)混油的影響研究還很少。(5)對(duì)于已建的管道若想進(jìn)行冷熱原油順序輸送必須測(cè)定管道周?chē)耐寥罍囟葓?chǎng),并進(jìn)行冷熱原油順序輸送數(shù)值計(jì)算,從而判斷其輸送可行性、冷熱油最佳出站溫度及順序輸送周期。6)冷熱原油順序輸送因其復(fù)雜性,致使計(jì)算過(guò)程仍需要進(jìn)行多方面的假設(shè)。因此,開(kāi)展先導(dǎo)性試驗(yàn),進(jìn)一步優(yōu)化計(jì)算模型是必要的。

[1] 吳玉國(guó).冷熱原油順序輸送技術(shù)研究[D].東營(yíng):中國(guó)石油大學(xué)(華東),2010.

[2] 吳明,周詩(shī)崠王雷,等.冷熱原油交替順序輸送優(yōu)化模型研究[J].石油化工高等學(xué)校學(xué)報(bào),2012,25(5):64-67. Wu Ming,Zhou Shidong,Wang Lei,et al.Research on optimization model of hot&cold crude oil sequential transportation[J].Journal of Petrochemical Universities,2012,25(5):64-67.

[3] 楊云鵬,劉寶玉,張玉廷,等.冷熱原油順序輸送溫度場(chǎng)波動(dòng)規(guī)律[J].遼寧石油化工大學(xué)學(xué)報(bào),2013,33(1):53-56. Yang Yunpeng,Liu Baoyu,Zhang Yuting,et al.Temperature field undulate rule of hot&cold crude oil sequential transportation[J].Journal of Liaoning Shihua University,2013,33(1):53-56.

[4] 邢曉凱,孫瑞艷,王世剛.西部原油管道冷熱油交替輸送過(guò)程溫度場(chǎng)[J].北京工業(yè)大學(xué)學(xué)報(bào),2013,39(4):614-618. Xing Xiaokai,Sun Ruiyan,Wang Shigang.Temperature field of hot&cold crude oil alternate transportation in western crude oil pipeline[J].Journal of Beijing University of Technology,2013,39(4):614-618.

[5] 韓東旭,王鵬,邵倩倩.冷熱油交替輸送熱力快速計(jì)算研究[J].工程熱物理學(xué)報(bào),2015,36(1):130-135. Han Dongxu,Wang Peng,Shao Qianqian.Fast thermal calculation research on hot&cold crude oil alternate transportation[J].Journal of Engineering Thermophysics,2015,36(1):130-135.

[6] 崔秀國(guó),張勁軍.冷熱油交替順序輸送過(guò)程熱力問(wèn)題的研究[J].油氣儲(chǔ)運(yùn),2004,23(11):15-19. Cui Xiuguo,Zhang Jinjun.Research on thermal problems of hot&cold crude oil sequential transportation[J].Oil&Gas Storage and Transportation,2004,23(11):15-19.

[7] 李傳憲,施靜.冷熱原油順序輸送過(guò)程的熱力分析[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,37(2):112-118. Li Chuanxian,Shi Jing.Thermal analysis on batch pipelining of cold and hot crude oil with different temperatures[J]. Journal of China University of Petroleum(Edition of Natural Sciences),2013,37(2):112-118.

[8] 王凱,張勁軍,宇波.原油管道差溫順序輸送水力—熱力耦合計(jì)算模型[J].油氣儲(chǔ)運(yùn),2013,32(2):143-151. Wang Kai,Zhang Jinjun,Yu Bo.Hydro-thermal coupling computational model of difference temperature of crude oil sequential transportation[J].Oil&Gas Storage and Transportation,2013,32(2):143-151.

(編輯 王戩麗)

Numerical Study on Dynamic Heat Transfer Characteristics of Cold and Hot Crude Oil in Long Distance Transportation

Qing Shuyong1,Dai Qing2,Guo Hongyan1,Du Mingjun3,Cui Qi3,Yang Zining3
(1.China Petroleum Engineering Co.Ltd.Beijing Company,Beijing100000,China; 2.China Petroleum Engineering Co.Ltd.Southwest Company,Chengdu Sichuan610000,China; 3.China Petroleum Engineering Co.Ltd.North Company,Renqiu Hebei062552,China)

Cold and hot crude oil orderly transportation is very strict to the temperature field of the soil,which is also the key for the hot oil to enter the station safely after the cold oil is transported.Based on heat transfer and fluid dynamics,the governing equations of flow and heat transfer of buried pipelines have been established,and the numerical simulation of the dynamic change process of cold and hot crude oil orderly transportation process along the pipeline in different axial position of the oil temperature and soil temperature field have been done.The results show that with the alternating cycle operation of temperature of the outbound crude oil,there is a periodic variation in the oil temperature of different locations along the pipeline and a certain range around the soil temperature field but a certain time or space lag.For the low throughput operation of the pipeline,more than a certain distance,the transmission temperature tends to be consistent.There is a main factor which should be considered in the pipeline safe operation plan that in the process of hot-cold oil batch transportation,hot oil inlet temperature drops to a minimum which is affected by the previous cold oil.

Long pipeline;Cold and hot crude oil;Order delivery;The dynamic heat transfer

TE832

:A

10.3969/j.issn.1006-396X.2017.01.018

1006-396X(2017)01-0092-05投稿網(wǎng)址:http://journal.lnpu.edu.cn

2016-11-08

:2016-12-18

中國(guó)石油工程設(shè)計(jì)公司管道技術(shù)改進(jìn)項(xiàng)目(2016JG-13)。

青樹(shù)勇(1983-),男,工程師,從事油氣管道及油氣田地面工程設(shè)計(jì)與研究;E-mail:cpebj-qingsy@163.com。

杜明俊(1983-),男,碩士,工程師,從事石油天然氣工程設(shè)計(jì)及系統(tǒng)仿真研究;E-mail:dmj260750009@163.com。