原油順序輸送水平管混油段的分析

王 琪，羅 浩，陳國平

（1. 廣東石油化工學(xué)院，廣東 茂名 525000； 2. 虎門港海灣石油倉儲碼頭有限公司， 廣東 東莞 523988）

原油順序輸送水平管混油段的分析

王 琪1，羅 浩2，陳國平3

（1. 廣東石油化工學(xué)院，廣東 茂名 525000； 2. 虎門港海灣石油倉儲碼頭有限公司， 廣東 東莞 523988）

針對目前從港口輸往煉油廠的油品多樣性，以大慶和勝利油品為例建立了水平原油管段順序輸送模型，應(yīng)用FLUENTT軟件進行了三維數(shù)值模擬，得出了混油濃度變化圖象和曲線，表明在原油管道順序輸送過程中粘度差的影響。幵總結(jié)出油品分散規(guī)律，幵給出了切割的長度，為找到最佳切割點提供依據(jù)。

順序輸送；混油段；切割點；數(shù)值模擬

隨著中國經(jīng)濟的収展，能源已經(jīng)成為其主要命脈，但原油的自產(chǎn)能力增長有限，所以對進口原油的依存度不斷加大，當(dāng)前我國原油進口除了俄羅斯原油直接通過輸油管道運輸外。其它均為通過油輪在港口卸油，上岸后經(jīng)過管道運往各個煉油企業(yè)。但是由于油的來源不同在順序輸送管道中接觸面處將產(chǎn)生一段混油。混油量與輸送油品的密度、粘度等物性參數(shù)有關(guān),也與管道的高程[1，2]，形狀有關(guān)，但多數(shù)管線都在水平直管段運行。以下以兩種油品順序輸送進行模擬[3]，計算混油段，為找到最佳切割點提供依據(jù)。

1 物理模型

熱油管道的尺寸：長度 L=20 m，直徑D=0.8 m(圖1)。.

圖1 直管模型Fig.1 Straight pipe model

流動介質(zhì):大慶原油和勝利原油等原油，原油的物理參數(shù)如表1所示:

表1 原油的物性參數(shù)Table 1 Crude oil physical property parameters

輸送速度： =2 m/s;操作溫度：原油管道一般都要加熱輸送，操作溫度為323 K

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

2.1 控制方程

當(dāng)以原油以流速=2 m/s，管內(nèi)流動為紊流。

假設(shè)原油的粘性為常數(shù)、不可壓流體，圓管光滑，則流動的控制方程如下：

1）質(zhì)量守恒方程

2）動量守恒方程

3) 湍動能方程：

4）湍能耗散率方程[4]：

式中：ρ—密度；

u、v、w —流速矢量在x、y和z方向的分量，p為流體微元體上的壓強。

方程求解：采用雙精度求解器，非定常流動，標(biāo)準(zhǔn) 模型，Simple算法。

2.2 定義混合組分相關(guān)物理性質(zhì)

1）密度

混油段截面密度變化遵循簡單的比例相加規(guī)律，即混油段密度為每種油品體積分?jǐn)?shù)與其密度乘積之和

式中:ρA, ρB, ρM—分別為A油品、B油品和混油的密度，kg/m3；CB—B油品的體積分?jǐn)?shù)。

2）粘度

混油粘度不可以用簡單的組分比例計算，但其組分仍接近于線性關(guān)系。

式中:VA，VB，VM—分別為 A 油品、 B 油品和混油的運動粘度，m2/s ；

αA,αB—待定因數(shù)。（在本文中對于原油可取αA=0.819 9，αB=-0.896 4）。

3 水平直管模型內(nèi)混油的形成過程

設(shè)定順序輸送管路為水平管，不考慮初始混油量的影響。順序輸送油品的流速為 2 m/s，初始時刻開始輸送第一種原油（大慶原油），到第 10 s時刻第一種油品正好充滿圓管模型，此時開始從入口輸送第二種原油（柴勝利原油），一直到第 20 s，第二種原油通過整個管道模型結(jié)束。

由于順序輸送混油的濃度分布隨著時間的延續(xù)而變化，故采用瞬態(tài)計算模型，相應(yīng)的壓力速度的耦合方法應(yīng)用simple 算法。在湍流條件下利用標(biāo)準(zhǔn)k ? ε方程模型封閉方程組，應(yīng)用有限體積法對方程進行離散。計算結(jié)果收斂較好。

在實際的操作中，由于原油在管道模型的法向輸送速度幵不均勻，在輸送過程中相互混合和擴散，逐漸形成混油，幵且混油量在輸送過程中隨時間增大，這也符合理論分析中關(guān)于混油變化和拖尾的現(xiàn)象的事實。

在順序輸送的直管模型內(nèi)，后行油品由于徑向的速度分布不均，管道中心的流動速度明顯大于壁面附近流動速度，導(dǎo)致其呈楔形進入前行油品，形成混油分界面，如圖2所示。在混油分界面上，兩種油品不斷地對流擴散和紊流擴散作用，逐漸形成具有一定厚度的混油層，在混油層的中心軸線位置，原油流體具有最大的流動速度，因此后行油品能夠以最快的速度朝軸向進行對流擴散，使混油頭不斷變細(xì)變長；而在管道的壁面位置，由于前行油品與壁面存在著粘性附著力，在粘性附著力的作用下，壁面附近的前行油品流動受到明顯的阻礙，因此低流速下的前行油品只能通過后行油品的沖刷作用和徑向?qū)α鲾U散的形式緩慢除去，這就導(dǎo)致混油拖尾現(xiàn)象的產(chǎn)生。在管道中心處混油軸向?qū)α鲾U散和壁面混油拖尾現(xiàn)象的共同作用下，隨著輸送時間的增大，混油段的距離也明顯增大[5]。

圖2 直管模型內(nèi)不同時刻的混油變化過程Fig.2 Mix oil change process of straight pipe model in different time

當(dāng)t=15 s時混油界面正好到達管子中段，為了觀察管子內(nèi)部的混油的相變化情況，本次試驗以中心點(L=10 m)為源點，以0.2 m為分度值進行前后截圖，觀察情況如圖3所示。在圖中直觀地反映了混油經(jīng)過各個水平管截面時的混油形態(tài)，同時也從側(cè)面說明了管內(nèi)混油徑向紊流擴散的存在。

圖3 t =15 s時水平直管模型內(nèi)不同位置的混油截面Fig.3 t = 15 s horizontal straight pipe model in different location

4 水平管內(nèi)混油長度的分析

在順序輸送中組分濃度為1%～99%之間的油段為混油段[6，7]，混油段所占輸油管道的長度為混油長度，如圖4所示，則混油長度的計算公式為：

式中:X1，X2—混油油頭和油尾所在位置的坐標(biāo)；L為混油長度。

圖4 混油形態(tài)示意圖Fig.4 Mix oil form

在本次實驗操作中，對于混油的位置主要是通過切割截面法確定混油的位置，幵通過混油長度的計算公式9進行計算，得出混油長度如表2所示。

表2 混油長度隨時間的變化Table 2 Mix oil length changes over time

從表2中可以収現(xiàn)，混油的長度隨輸送時間的增大而增大，但是幵不是簡單的線性關(guān)系。因為在不考慮實驗誤差的情況下，等時間差下的混油長度變化幵不相等或者相似，而是逐漸增大的，由此可以確定混油長度隨時間或者輸送距離的變化具有冪函數(shù)關(guān)系。

5 結(jié) 論

利用flunet有限元分析軟件可以對原油順序輸送混油進行精準(zhǔn)的模擬,利用豐富的彩色云圖可以形象直觀地實現(xiàn)從三維圖上表現(xiàn)混油變化情況,為原油順序輸送切割點的確定找到了理論依據(jù)和分析方法。

［1］趙會軍, 張青松, 張國忠, 等. 基于PHOENICS 的順序輸送大落差管道混油研究[ J] . 江漢石油學(xué)院學(xué)報, 2008, 28( 4) : 139-142.

［2］康正凌，宮敬，嚴(yán)大凡．管道順序輸送高差混油研究[J]．石油學(xué)報，2003，24(4)：94-98．

［3］趙會軍, 張青松, 張國忠, 等. 基于PHOENICS 順序輸送管道內(nèi)兩種油品傳質(zhì)數(shù)值模擬[J] . 北京化工大學(xué)學(xué)報, 2006, 33 ( 6) : 22- 25.

［4］劉強，王樹立，李恩田,等．利用湍流模式理論模擬順序輸送混油過程[J]．油氣儲運，2006，25（11）：22-26

［5］趙會軍, 張青松, 張國忠, 等. 基于PHOENICS 順序輸送管道混油濃度的計算[ J] . 油氣儲運, 2007, 26( 2) : 43- 46.

［6］慕希茂, 郭躍平, 李保國, 等. 順序輸送停輸時的混油量計算油氣儲運[J]. 油氣儲運,2001, 20( 9) : 30- 31.

［7］呂穩(wěn)高, 馮先強1 淺談原油順序輸送過程中混油量的計算[ J ] . 天然氣與石油, 2007, 25 ( 1) : 11- 13.

Analysis on Mixed Oil Segment in the Horizontal Pipe During Crude Oil Batch Transportation

WANG Qi1，LUO Hao2，CHENG Guo-ping3
（1. Guangdong Petroleum Chemical University, Gguangdong Maoming 525000，China；2. Humen Port Gulf Oil Storage Terminal Co., Ltd., Guangdong Dongguan 523988，China）

In view of the diversity of product from the port to the refinery, taking Daqing and Shengli crude oil as examples, horizontal pipe model of crude oil transportation was established, then three-dimensional numerical simulation was carried out with FLUENTT software, the image and curve of the mixed oil concentration change were obtained as well as effect of crude oil viscosity spread on the crude oil batch transportation. The oil distributing rule was summed up, and the length of the cutting was determined, which could provide the basis for finding the best cutting point.

Batch transportation; Mixed oil segment; Cutting point; Numerical simulation

TE832

A

1671-0460（2014）10-2153-03

2014-03-31

王琪（1982-），女，遼寧撫順人，講師，碩士，2009年畢業(yè)于遼寧石油化工大學(xué)油氣儲運專業(yè)，研究方向：油品流場，溫度場的數(shù)值模擬。E-mail：wangqi1982.good@163.com。