何顯榮 張引弟 楊建平 伍麗娟

1．長江大學(xué)石油工程學(xué)院，湖北武漢430100;

2．中國石油遼河油田公司SAGD開發(fā)項(xiàng)目管理部，遼寧盤錦124010

不同傾角下成品油順序輸送混油數(shù)值模擬

何顯榮1張引弟1楊建平2伍麗娟1

1．長江大學(xué)石油工程學(xué)院，湖北武漢430100;

2．中國石油遼河油田公司SAGD開發(fā)項(xiàng)目管理部，遼寧盤錦124010

為了能更好地認(rèn)識管道成品油順序輸送停輸時混油段變化特性，應(yīng)用Fluent軟件組分輸送模型，對不同傾角下的成品油順序輸送停輸進(jìn)行了數(shù)值模擬，并對模擬結(jié)果和導(dǎo)出數(shù)據(jù)進(jìn)行分析。研究結(jié)果表明:停輸時，重力和油品間的密度差是造成混油段特性變化的主要因素，管道傾角越大，柴油在重力勢能高處時產(chǎn)生的自然對流使混油段長度明顯增大且所需時間更短，汽油在重力勢能高處時，管道傾角的增大對混油段的增加有抑制作用且混油段濃度分層與管道傾角一致。研究結(jié)果對減少和預(yù)防停輸時混油段的產(chǎn)生和增加有一定借鑒作用。

順序輸送;停輸;不同傾角;混油特性;數(shù)值模擬

0 前言

成品油順序輸送過程中由于種種原因會造成管道停輸［1－2］，而輸送管道常經(jīng)過地形復(fù)雜有較大落差的地段［3］，一般為減少混油段產(chǎn)生，管道會在紊流狀態(tài)下運(yùn)行，其雷諾數(shù)一般大于104［4－5］。而當(dāng)管道停輸時，管道內(nèi)油品容易出現(xiàn)層流［6］，管道中心區(qū)域與管壁附近流速差增加，使混油段增加［7］，如蘇聯(lián)的古比雪夫—勃克斯克成品油管段，停輸幾小時的混油量比連續(xù)輸送216 h形成的混油量還要多1倍［8］。而在不同傾角和輸送順序下，油品間由于物性參數(shù)不同會使混油特性發(fā)生怎樣變化。在前人研究中，對管道順序輸送模擬大多采用VOF模型［9－11］，此模型多適用于兩種不同相態(tài)流體存在自由分界面的流動情況模擬，對兩種相同相態(tài)油品間混油特性的模擬效果并不理想［12］。針對以上問題，本文采用Fluent軟件組分輸送(Species Transport)模型，建立了成品油順序輸送控制方程，就不同傾角和輸送順序下管道停輸對混油段特性進(jìn)行模擬分析，對管道停輸時混油段的控制提供理論指導(dǎo)。

1 模型建立

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

管道中混油濃度分布隨時間而變化，故對管道紊流計(jì)算采用雷諾時均法［13］，通過Boussinesq假設(shè)，引入脈動造成的Reynolds應(yīng)力和質(zhì)量輸送流率。選用物質(zhì)擴(kuò)散模型，在考慮重力條件下，建立管道輸送與停輸?shù)臅r均控制方程組求得。

1.1.1 連續(xù)性方程

1.1.3 質(zhì)量輸送方程

式中:xi、xj為空間坐標(biāo);ui、uj為空間時均速度變量，m/s; u'i、u'j為脈動速度分量，m/s;ρ為密度，kg/m3;gi為質(zhì)量力矢量，m/s2;P為壓力，Pa;μ為層流動力黏度，Pa· s;c'為前行油品脈動濃度，%;c為前行油品的體積時均濃度，%;Pri(c)為層流Schmidt數(shù)。

1.1.4 物質(zhì)擴(kuò)散模型

Fluent中對物質(zhì)擴(kuò)散守恒方程求解i組分的質(zhì)量擴(kuò)散速率的方程［14］為:

式中:μt為紊流黏度，Pa·s;Di，m為物質(zhì)i的擴(kuò)散系數(shù); Sct為紊流Schmidt數(shù)，默認(rèn)值取0.7;γi為物質(zhì)i的質(zhì)量分?jǐn)?shù)，%。

式(1)～(4)應(yīng)用于Fluent軟件中，對不同傾角下管道混油停輸工況進(jìn)行模擬計(jì)算。

1.2 物理模型

“哦，翠絲，僵尸人吧？”艾瑞克說著，突然咧嘴假笑了一下，這笑扯動唇環(huán)，環(huán)孔一下子被拉寬了，那樣子讓我畏縮了一下，“我倒要看看，你能撐多久?！?/p>

本文主要討論在不同傾角下停輸管道中混油段變化特性的影響，以管徑為0.5 m，管長為10 m的管道作為研究對象，在傾角下管道內(nèi)油品的重力分解圖見圖1。

圖1 傾角下管道內(nèi)油品重力分解圖

1.3 油品參數(shù)

以常見的0#柴油和90#汽油作為輸送油品［15］，其物性參數(shù)見表1。

表1 油品物性參數(shù)

2 數(shù)值模擬

利用Fluent軟件模擬組分為成品汽油和柴油的順序輸送，應(yīng)用PISO算法實(shí)現(xiàn)對壓力速度耦合數(shù)值計(jì)算求解，選擇雙方程湍流模型及標(biāo)準(zhǔn)固壁無滑移邊界條件，默認(rèn)操作壓力，忽略溫度影響。管道進(jìn)口條件設(shè)為速度進(jìn)口，速度均設(shè)為1 m/s;管道出口條件設(shè)為自由出口。不同傾角下重力因素的影響由x、y軸的重力分量替代。在前行油品充滿管道的情況下，輸送后行油品。油品初始速度均為1 m/s，考慮重力因素(不同傾角下)。當(dāng)混油段行進(jìn)至管內(nèi)3 m時，油品停止輸送(設(shè)此時t=0 s)。

2.1 有無重力下管道停輸模擬

首先研究重力對混油特性在水平管內(nèi)的變化情況，先輸汽油后輸柴油。管道充滿汽油后，柴油在管道內(nèi)行進(jìn)3 s時停輸(t=0 s)，混油變化情況和混油濃度分布見圖2。

圖2 有無重力下水平管道各時間點(diǎn)的停輸模擬

由圖2知，不考慮重力的情況下未停輸t1=0 s時，在混油段任一截面的中心部位，濃度分布都較均勻，管壁處的層流邊界層存在較大濃度梯度，與理論描述一致［16－17］，停輸后，由于濃度差導(dǎo)致分子擴(kuò)散作用明顯，混油沿徑向和軸向具有濃度差的地方進(jìn)行擴(kuò)散，管壁處的層流邊界層濃度梯度逐漸平衡，混油段徑向濃度差逐漸減小，整個混油段變得均勻，混油段長度逐漸增長一段時間后停止;考慮重力的情況下未停輸t2=0 s時，由于重力和不同油品間密度差，使得密度大的后行柴油下沉，密度小的前行汽油上浮，停輸t2=5 s時，這種對流效果更加明顯，很快就在管內(nèi)形成較長的混油段。

為了更清晰地體現(xiàn)混油濃度沿管線的分布情況，有無重力因素管道停輸時，混合油品沿管道中心軸向分布濃度從模擬結(jié)果中導(dǎo)出。

圖3為混合油品沿管道中心軸向分布濃度，圖3中無重力停輸t1=5 s時，混油段由于慣性作用，混油段濃度還往前運(yùn)行了一小段，而在停輸t1=10 s后，管道內(nèi)剩余紊流波動消失，油品逐漸穩(wěn)定，依靠濃度差的分子擴(kuò)散作用，兩種油品之間開始向管內(nèi)四周進(jìn)行分子擴(kuò)散，t1=20 s之后，混油沿管道中心軸向濃度曲線基本不再發(fā)生變化，擴(kuò)散速率隨時間延長逐漸平緩，這是因?yàn)榉肿訑U(kuò)散所引起的濃度變化已近極限;有重力停輸t2=0 s時，混油沿管道中心軸向濃度相比同時刻無重力情況時，混油段前部濃度大，后部濃度小，這是因?yàn)橹亓κ骨靶衅蜕细『笮胁裼拖鲁?，混油段延長，混油濃度特性傾斜，在停輸后，這種依靠重力作用產(chǎn)生的對流效果愈發(fā)明顯，在t2=3、5 s時，管道內(nèi)混油沿管道中心軸向濃度迅速擴(kuò)散且增幅明顯，說明停輸時重力對混油量的影響很大，有必要進(jìn)行詳細(xì)分析。

圖3 混合油品沿管道中心軸向分布濃度

2.2 不同傾角下未停輸?shù)臄?shù)值模擬

不同傾角下，油品順序輸送未停輸時，混油段特性也有所不同。模擬了傾角為30°、45°、60°、90°下管道未停輸?shù)那闆r，混合油品濃度沿混油管道中心徑向分布，模擬結(jié)果見圖4～5。

前行汽油后行柴油未停輸時，后行柴油向下嵌入前行汽油中，隨著傾角增大，重力在管道徑向的作用力逐漸減弱，使得柴油沿管道徑向下沉嵌入汽油的趨勢減小，混油段濃度在管道中心徑向上越均勻，混油段濃度分層與管道傾角角度大體一致，混油段長度隨傾角增大略有減小;前行柴油后行汽油未停輸時，后行汽油上浮嵌入到前行柴油中，混油頭偏離管道中心，在不同傾角下，混油尾長度和濃度變化并不明顯，說明重力對紊流輸送時的影響很小，這與理論一致［18－20］，傾角為90°時，混油濃度沿管道截面分布均勻，與前行汽油后行柴油的情況相比，由于柴油的黏度比汽油大得多，在管壁邊界層處濃度變化梯度明顯。

兩種輸送順序下對比可知，在有傾角的紊流輸送時，前行柴油后行汽油的混油比前行汽油后行柴油形成的混油量要多，這是由于前行柴油后行汽油時，后行汽油嵌入前行柴油中，管道紊流輸送加劇了混油段的增長。

圖4 不同傾角下管道停輸前混油特性模擬結(jié)果

圖5 混油段中心濃度沿管道徑向分布圖

2.3 不同傾角下停輸后的數(shù)值模擬

1)當(dāng)前行汽油后行柴油的管道停輸時，不同傾角下管道內(nèi)混油段混油特性模擬見圖6，混合油品濃度沿管道中心徑向分布見圖7。

處于高重力勢能的是質(zhì)量較輕的汽油，較重的柴油處于低重力勢能處，不同傾角下混油段沒有明顯的由于重力引發(fā)的對流情況，但管道內(nèi)剩余的紊流脈動和分子間的擴(kuò)散作用使混油段擴(kuò)散有明顯波動，隨著傾角增大，重力分量偏向軸向方向，混油濃度在管道中心徑向上趨于平穩(wěn)，混油段長度也隨之減小?；煊投坞S著停輸時間增長，混油段濃度逐漸分層且分層角度與管道傾斜角度基本一致，此時混油段長度維持穩(wěn)定，管道傾角增加對混油段增長起抑制作用。

圖6 前行汽油后行柴油時不同傾角下管道停輸后混油特性模擬

圖7 前行汽油后行柴油時不同傾角下停輸管道徑向混油濃度

2)當(dāng)前行柴油后行汽油的管道停輸時，不同傾角下管道內(nèi)混油特性模擬見圖8?；煊蜐舛妊毓艿乐行妮S向分布見圖9。

前行柴油處于重力勢能高處，后行汽油處于重力勢能低處，未停輸時，混油段受重力的影響使柴油滲流到汽油液層之下，混油段長度比前行輸汽油后行輸柴油時長(見圖4)，不同傾角下管道中混油段長度大致相同。停輸1 s時，不同傾角下混油段長度差距不明顯，停輸3 s時，傾角越大，下沉到后行汽油中的混油頭越大，混油段擴(kuò)散速度也越快，擴(kuò)散范圍越廣，當(dāng)傾角90°時，混油擴(kuò)散方式略有不同，混油沿管壁向下延伸，在管道中心軸處，混油濃度相比其他傾角下反而較小，起始混油段長度比其他傾角時小，但混油段擴(kuò)散速度明顯比其他傾角時快。因此在前行柴油后行汽油時，傾角增加對停輸管道混油段增長起促進(jìn)作用。

圖8 前行柴油后行汽油時不同傾角下管道停輸后混油特性模擬

圖9 前行柴油后行汽油不同傾角下管道中心軸向混油濃度

3 結(jié)論

1)導(dǎo)致順序輸送管道停輸混油段增長的主要因素在于重力作用下油品間的密度差引起的自然對流，不考慮重力時，停輸管道混油量的增加主要是由于油品間濃度差導(dǎo)致的分子擴(kuò)散造成的，這種混油量的增速緩慢，且在一段時間后，分子擴(kuò)散達(dá)到極限，混油段長度將不再變化。

2)管道未停輸，前行汽油后行柴油時，柴油嵌入前行汽油中，混油段濃度分層基本與管道傾角一致，混油段長度隨傾角增大而略有減小;前行柴油后行汽油時，混油段長度和濃度變化不明顯，混油段濃度分布隨傾角增大而愈發(fā)均勻。傾角管道紊流輸送時，前行汽油后行柴油所形成的混油量比前行柴油后行汽油所形成混油量要多，這是由于前行輸柴油后行輸汽油時，后行汽油嵌入前行柴油中，管道紊流輸送又加劇了混油段的增長。

3)管道停輸后，前行汽油后行柴油時，傾角越大，混油段長度越小，混油段濃度分布越均勻，紊流脈動越不明顯，混油段濃度分層角度與管道傾角基本一致，傾角增加對停輸管道混油段增長起抑制作用;在前行柴油后行汽油時，柴油向下嵌入汽油中，混油形成自然對流，管道傾角越大，下嵌的混油頭濃度越大，混油段增速越快。

4)如需對輸送管道停輸時，使密度較小的油品處于重力勢能高處，密度較大的油品處于重力勢能低處，可以有效地減小停輸造成的混油量增大。

［1］楊筱蘅．輸油管道設(shè)計(jì)與管理［M］．東營:中國石油大學(xué)出版社，2011．Yang Xiaoheng．Design and Management of Oil Pipeline［M］．Dongying:China University of Petroleum Press，2011．

［2］張其敏，陳寧．成品油順序輸送管道優(yōu)化設(shè)計(jì)［J］．油氣儲運(yùn)，2004，23(3):5－7．Zhang Qimin，Chen Ning．Optimization Design on Products Batch Transportation Pipeline［J］．Oil＆Gas Storage and Transportation，2004，23(3):5－7．

［3］康正凌，宮敬，嚴(yán)大凡．管道順序輸送高差混油研究［J］．石油學(xué)報(bào)，2003，24(4):94－98．Kang ZhengLin，Gong Jing，Yan Dafan．Research on Contamination Caused by Elevation Difference in Batching Transportation［J］．ACTA Petrolei Sinica，2003，24(4):94－98．

［4］尚鳳山．原油管道的順序輸送［J］．油氣儲運(yùn)，2004，23(1):9－10．Shang Fengshan．Batch Transportation of Crude Oil［J］．Oil＆Gas Storage and Transportation，2004，23(1):9－10．

［5］郭小倩．慶哈輸油管道安全停輸再啟動技術(shù)研究［D］．大慶:東北石油大學(xué)，2014．Guo Xiaoqian．A Technical Research on the Safety Shutdown and Restart of Qingha Oil Pipeline［D］．Daqing:Northeast Petroleum University，2014．

［6］蔣新國，劉愛虢，丁啟敏．熱油管道停輸再啟動過程［J］．天然氣與石油，2005，23(2):25－27．Jiang Xinguo，Liu Aiguo，Ding Qimin．Restarting Process of Hot Oil Pipeline after Shutdown［J］．Natural Gas and Oil，2005，23(2):25－27．

［7］王昆，陳保東，郭淑娟，等．管道順序輸送產(chǎn)生混油的問題研究［J］．油氣儲運(yùn)，2007，26(8):8－12．Wang Kun，Chen Baodong，Guo Shujuan，et al．A Study on the Product Oil Contamination Caused by Batching Transportation［J］．Oil＆Gas Storage and Transportation，2007，26 (8):8－12．

［8］康正凌，宮敬，嚴(yán)大凡．成品油管道輸送高差混油模型研究［J］．石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2003，27(6):65－67．Kang Zhenglin，Gong Jing，Yan Dafan．Development of Model for Contaminated Oil Product Cased by Elevation Difference in Multiproduct Pipelines［J］．Journal of the University of Petroleum:Edition of Natural Science，2003，27 (6):65－67．

［9］呂坦．成品油管道順序輸送混油分析與模擬研究［D］．西安:西安石油大學(xué)，2014．Lü Tan．The Analysis and Simulation Study About Batching Transportation of Mix Oil by Products Pipeline［D］．Xi'an: Xi'an Shiyou University，2014．

［10］鹿欽禮，杜明俊，馬貴陽，等．變徑管道對冷熱原油順序輸送混油的影響［J］．油氣儲運(yùn)，2011，30(3):183－186．Lu Qinli，Du Mingjun，Ma Guiyang，et al．Influence of Reduced Pipeline on the Batch Transportation of Cold-Hot Contaminated Crude Oil［J］．Oil＆Gas Storage and Transportation，2011，30(3):183－186．

［11］杜明俊，馬貴陽，陳笑寒．凍土區(qū)埋地?zé)嵊凸艿劳］敎亟禂?shù)值模擬［J］．天然氣與石油，2010，28(4):54－57．Du Minjun，Ma Guiyang，Chen Xiaohan．Numerical Simulation for Temperature Drop of Buried Hot Oil Pipeline during Shutdown in Permafrost Region［J］．Natural Gas and Oil，2010，28(4):54－57．

［12］劉強(qiáng)，王樹立，李恩田，等．重力場對管內(nèi)液液間傳質(zhì)影響的數(shù)值模擬［J］．化學(xué)工程，2007，35(5):9－12．Liu Qiang，Wang Shuli，Li Entian，et al．Numerical Simulation of the Effects of Gravitational Field on Mass Transfer BetweenTwoLiquidsinPipeline［J］．Chemical Engineering，2007，35(5):9－12．

［13］趙會軍，張青松，張國忠，等．停輸對順序輸送管道混油影響的PHOENICS模擬［J］．西南石油大學(xué)學(xué)報(bào)，2007，29(4):149－151．Zhao Huijun，Zhang Qingsong，Zhang Guozhong，et al．Phoenics Simulation for the Effect of Shut-Down to Contamination in Batch Transportation［J］．Journal of Southwest Petroleum University，2007，29(4):149－151．

［14］陶文銓．?dāng)?shù)值傳熱學(xué)(第二版)［M］．西安:西安交通大學(xué)出版社，2001．Tao Wenquan．Numerical Heat Transfer(EditionⅡ)［M］．Xi'an:Xi'an Jiaotong University Press，2001．

［15］趙海燕．順序輸送混油的CFD模擬［D］．大慶:大慶石油學(xué)院，2010．Zhao Haiyan．CFD Simulation of Contamination on Batch Transportation Pipeline［D］．Daqing:Northeast Petroleum University，2010．

［16］李欣澤．成品油管道的混油長度計(jì)算方法［J］．油氣儲運(yùn)，2015，34(5):497－499．Li Xinzhe．Calculation Method for Length of Oil Mixtures in ProductPipelines［J］．Oil＆GasStorageand Transportation，2015，34(5):497－499．

［17］劉強(qiáng)，王樹立，趙會軍，等．停輸對混油量影響數(shù)值模擬研究［J］．油氣儲運(yùn)，2007，26(3):13－15．Liu Qiang，Wang Shuli，Zhao Huijun，et al．Study on Numerical Simulation of Contaminated Volume for Shutdown BatchingPipeline［J］．Oil＆GasStorageand Transportation，2007，26(3):13－15．

［18］宮敬，王琴，王衛(wèi)東，等．成品油管道的混油計(jì)算方法研究［J］．管道技術(shù)與設(shè)備，2013，(1):1－3．Gong Jing，Wang Qin，Wang Weidong，et al．Study of the Calculation Method for Mixing Volume in Products Pipeline［J］．Pipeline Technology and Equipment，2013，(1):1－3．

［19］劉勇峰，吳明，張宏濤，等．裸露管線停輸溫降規(guī)律數(shù)值模擬［J］．天然氣與石油，2011，29(5):18－20．Liu Yongfeng，Wu Ming，Zhang Hongtao，et al．Numerical SimulationofTemperatureDropofExposedPipeline Shutdown［J］．Natural Gas and Oil，2011，29(5):18－20．

［20］賈永英．埋地含蠟原油管道運(yùn)行及停輸再啟動過程研究［D］．大慶:大慶石油學(xué)院，2009．Jia Yongying．The Research on Operation and Shutdown and Restart-Up of Buried Waxy Crude Oil Pipeline［D］．Daqing:Northeast Petroleum University，2009．

中國石油工程建設(shè)公司海外再攬大單

阿爾及利亞首都阿爾及爾當(dāng)?shù)貢r間2016年11月6日，中國石油工程建設(shè)公司與阿爾及利亞國家石油公司正式簽約阿爾及爾煉廠改擴(kuò)建項(xiàng)目，項(xiàng)目合同簽約價5.6億美元。這是繼10月在伊拉克攬得殼牌2.8億美元工程大單之后，中國石油工程建設(shè)公司在國際市場開發(fā)上取得的又一亮眼成績。

阿爾及利亞國家石油公司是非洲第一大石油公司，其子公司阿爾及爾煉油廠始建于1964年，主要處理來自阿國南部沙漠地區(qū)的原油。2010年，隨著油價高企，阿爾及利亞為增加對歐洲及國內(nèi)市場的成品油供應(yīng)，擬對煉廠的常壓蒸餾裝置、氣體廠、催化重整裝置等設(shè)施進(jìn)行翻新改建，同時新建一些裝置，擴(kuò)建后的年處理能力將達(dá)360×104t。

石油石化產(chǎn)品收入占阿國國家財(cái)政預(yù)算的95%、國際貿(mào)易收入的65%，是阿國重要經(jīng)濟(jì)命脈。在當(dāng)前世界經(jīng)濟(jì)復(fù)蘇乏力、國際油價低谷徘徊的大環(huán)境下，阿爾及利亞也遭遇嚴(yán)重的經(jīng)濟(jì)發(fā)展問題。因此，阿國希望加快煉廠改擴(kuò)建工程進(jìn)度，盡快投產(chǎn)，增加石油石化產(chǎn)品的對外出口，減緩國際石油價格下跌對阿國經(jīng)濟(jì)的沖擊。

中國石油工程建設(shè)公司自2004年進(jìn)入阿爾及利亞市場，先后成功實(shí)施了凝析油煉廠、圖瓦油田項(xiàng)目和水利部的水泵站項(xiàng)目，目前正在執(zhí)行阿爾及利亞國家石油公司油泵站項(xiàng)目。通過4個項(xiàng)目的實(shí)施，中國石油工程建設(shè)公司在阿國積累了豐富的EPC項(xiàng)目經(jīng)驗(yàn)，憑借在工期、成本和質(zhì)量安全等方面的優(yōu)異業(yè)績和顯著優(yōu)勢，在阿國樹立了良好的品牌形象，成為受邀參加項(xiàng)目投標(biāo)并最終戰(zhàn)勝韓國的三星、韓華等競爭對手，成功中標(biāo)。

(曾妍摘自中國石油新聞網(wǎng))

10.3969/j．issn．1006－5539.2016.06.004

2016－05－27

中國石油科技創(chuàng)新基金(2015 D－5006－0603);長江青年科技創(chuàng)新團(tuán)隊(duì)基金(2015 cqt 01)

何顯榮(1992－)，男，湖北荊州人，碩士研究生，主要從事多相輸送及流體傳質(zhì)傳熱。