彎管順序輸送影響因素研究

王 升，王 岳，劉麗艷

（遼寧石油化工大學(xué)石油天然氣工程學(xué)院， 遼寧 撫順 113001）

由于在成品長輸管道的線路上會出現(xiàn)各種工況及各種復(fù)雜管道，使得無法建立和研究它的確切模型。而采用數(shù)值模擬計算以其投資少、見效快的優(yōu)勢在研究及模擬方面占據(jù)了越來越重要的地位。其中，順序輸送過程中的油品經(jīng)過彎管是一種比較普遍的現(xiàn)象[1]。因為彎管在不同情況下其影響因素有很多，本文以90°彎管為例，利用CFD軟件進行了模擬計算，分析了流速、粘度和彎徑比在彎管中混油段的影響。

1 物理模型

1.1 物理模型及邊界條件

采用 CFD軟件建模并對彎管條件下混油進行數(shù)值模擬計算。對彎管上游和下游長度均為5 m，管徑600 mm，彎徑比為3的90°水平彎管進行建模（如圖1所示）。管內(nèi)區(qū)域采用結(jié)構(gòu)化六面體網(wǎng)格進行單元劃分，在貼近管壁出的邊界層區(qū)將網(wǎng)格細化（如圖2所示）。且本文列出的所有圖形均是以模型左端為入口，模型上端為出口。

本文是以汽油先行柴油后行的輸送過程為研究條件，采用多相流（VOF）進行模型計算[2]，入口邊界設(shè)置為速度入口（velocity）,出口邊界設(shè)置為自由出口（outflow）。為了使分析更加準確，先將入口處的柴油體積分數(shù)設(shè)置為 0，等前行的汽油運行一段時間直至前行油品在管道內(nèi)充分運動后，再將后行的柴油的體積分數(shù)改為 1，從而實現(xiàn)了這種交替輸送的過程。

圖1 彎管順序輸送的三維幾何模型Fig.1 The three-dimensional geometric model of batch pipelining transportation in an angle branch conditions

圖2 管道入口處邊界層放大示意圖Fig.2 The enlarged view of entrance to the boundary layer pipes

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

由于管道內(nèi)截面流速分布不均勻，管道內(nèi)紊流流動時混油量要比層流流動時小很多。故為了減少混油損失，應(yīng)使管道在紊流狀態(tài)下運行。該紊流計算采用雷諾時均法[3,4]，建立成品油順序輸送管道混油計算的時均方程組。本文不考慮重力因素，并且因為輸送過程為等溫輸送，所以不引入重力項和能量方程。非穩(wěn)態(tài)的油品順序輸送的數(shù)學(xué)模型如下：

連續(xù)性方程：

動量方程（N-S方程）：

組分方程：

紊流模型采用k-e 模型。

k方程：

e方程：

式中：ρ—油品密度；

ui、uj—油品速度；

xi、xj—距離；

p—壓力；

Prt（C）—為紊流Schmidt數(shù)；

Prl（C）—層流Schmidt數(shù)；

ut、ul—紊流和層流運動粘度；

t —時間；

gi—重力加速度；

c —前行油品體積分數(shù)；

K —紊流脈動動能；

e —紊流脈動動能耗散率。

1.3 物性參數(shù)

常溫下0#柴油和93#汽油的密度分別取800、720 kg/m3，運動粘度分別取 6×10-6、4×10-6m2/2。混油段截面密度變化遵循簡單的比例相加規(guī)律，但是混油粘度不可以用簡單的組分比例計算，但其組分仍接近于線性關(guān)系。

式中：ρA、ρB、ρm—分別為A油、B油和混油的密度；

νA、νB、νm—分別為A油、B油和混油的運動粘度；

CA、CB—分別為A油、B油的體積分數(shù)；

αA、αB—待定系數(shù)。

2 模擬計算結(jié)果及分析

2.1 流速變化

在前行油品為93#汽油，后行油品為0#柴油，流速分別為1、3 m/s的情況下使用CFD軟件建模并對混油的特性數(shù)值作模擬計算。經(jīng)計算，在突擴管中雷諾數(shù)最小值為100 000，為紊流流動。當混油頭處于同一位置時混油段的俯視圖如圖3所示。

圖3 速度不同時混油段濃度分布示意圖Fig.3 The concentration distribution diagram of contaminated segments at different speeds

可以看到，當油品按照不同的速度輸送時，混油段在管內(nèi)的形態(tài)不一致。輸送速度小的混油油頭更加尖銳。這是由于雷諾數(shù)隨速度的增大而增大，雷諾數(shù)越大混油截面上的分布就越均勻，強烈的紊流擴散使得液體質(zhì)點發(fā)生劇烈的橫向交換，雷諾數(shù)越大，紊流擴散就越強烈，而此時分子擴散的作用則可以忽略不計[5,6]。同時，由于受到彎管曲率的影響，在離心力的作用下流體向著曲率半徑較大的彎曲部分外側(cè)方向移動，從而導(dǎo)致大量前行油品堆積在外側(cè)彎管的壁面附近。位于管道外壁面的前行油品相較于靠近彎管內(nèi)側(cè)的油品流動要相對滯后。而在流速不同的情況下，速度較大的輸送方式相對于速度較小的在管道彎曲處產(chǎn)生的向心力的差值更大。所以輸送速度較大的輸送方式相比速度小的混油油頭更加靠近彎管內(nèi)側(cè)部分。

2.2 粘度變化

在前行油品為 93#汽油，后行油品的粘度分別為 和 的 0#柴油，流速都為 1 m/s 的情況下使用CFD軟件建模并對混油的特性做模擬計算。當運行13 s時混油段的俯視圖如圖4所示。

可以看到，在相同的輸送時間下，后行油品粘度較大的比黏度較小的混油油頭的位置相對滯后，這是由于粘度大的油品分子與管壁的附著力較大使得油品輸送相對緩慢所致[7]，此處不詳細論述。

2.3 彎徑比變化

將彎徑比為1、3和5的這3種管道以前行汽油后行柴油的輸送順序，在1 m/s的流速下輸送。當油品輸送至下游管段的同一位置時，三種情況下的俯視圖如圖5所示。

圖5 彎徑比不同時混油濃度分布示意圖Fig.5 The concentration distribution diagram of different bend diameter ratio

可以看到，在彎徑比不同的情況下混油油頭的尖銳程度是隨著彎徑比的增加而更加尖銳。這是因為彎徑比增加，則管道彎曲處所產(chǎn)生的向心力和油品所經(jīng)歷的路程都會隨之增大，

致使紊流擴散程度減弱，液體質(zhì)點所發(fā)生的橫向交換減弱，由頭則更加尖銳。同時，觀察這幾種情況，彎徑比為3的情況下混油油頭相較于其他兩種情況更加更加靠近彎管內(nèi)側(cè)部分。這是因為，當有品流經(jīng)彎管時，大量前行油品堆積在外側(cè)彎管的壁面附近，彎徑比較大的彎管中油品前行油品所經(jīng)歷的路程較長，后行油品還沒有足以填補彎曲部分，致使混油油頭更向彎曲部分偏移。但是，在彎徑比為5的情況下，后行油品基本上完全替代了前行油品在彎管中的位置。后行油品此時已經(jīng)大量堆積在彎曲部分，使原先偏向于彎管內(nèi)側(cè)的油頭重新調(diào)整回來（如圖6所示）。

圖6 彎徑比為5時后行油品運行20s混油濃度分布示意圖Fig.6 The concentration distribution diagram of bend diameter ratio of 5 after trailing oil product running 20 seconds

在彎徑比不同的情況下，由于混油的油頭和油尾的位置相對于下游管道的軸線位置一直在發(fā)生變化，所以混油在同樣時間范圍內(nèi)的長度變化同樣值得考慮，為了使下面的研究結(jié)果更加明顯，將彎徑比為1、3和5的管道的下油管道分別加長15、10和5 m。當彎徑比為5的情況下，混油段進入下游管段的時間最晚，故以其進入下游管段的時間為標準，其時間為第19 s，在19～32 s內(nèi)記錄結(jié)果如圖7所示：

圖7 19s～32 s三種不同彎徑比管道內(nèi)的混油段長度Fig.7 The contaminated length of three different bend diameter ratio in 19 second to 32 second

可以看出，三種情況下混油段長度的變化趨勢基本一致，且隨著彎徑比的增加，混油段長度越長。這是因為流線是一條光滑曲線，它既不能轉(zhuǎn)折也不能相交，所以在油品流經(jīng)彎管彎曲處時，壁面附近會產(chǎn)生兩個漩渦區(qū)，彎曲部分內(nèi)側(cè)的漩渦區(qū)的形成是由于流體流動時的慣性使其與內(nèi)壁脫離造成的，而彎曲部分外側(cè)的漩渦的形成主要是因為流體與壁面撞擊形成的。隨著彎徑比的增大，在彎管彎曲部分的油品的流向變化越平緩，則使得漩渦區(qū)增大。從而維持這兩個漩渦區(qū)所需的能量增大而帶動流體主流束的速度減小，繼而混油段長度增加。

3 結(jié) 論

建立了紊流順序輸送模型，采用了CFD軟件對汽油前行柴油后行情況下針對彎管順序輸送條件的流速、粘度和彎徑比的改變進行模擬計算。結(jié)果表明，流速越小混油油頭越尖銳，并且流速越大混油油頭越靠近管道彎曲處；而后行油品粘度較大的比粘度較小的混油油頭的位置相對滯后；彎徑比較大的管道混油流經(jīng)下游管段時混油段長度越長。與理論結(jié)果基本相符。為以后研究順序輸送混油段計算和檢測方面提供了幫助和參考。

［1］ 楊筱蘅.輸油管道設(shè)計與管理[M].東營：中國石油大學(xué)出版社，2006.

［2］ 蔣仕章，蒲家寧．成品油順序輸送時混油粘度的計算與誤差分析[J]．油氣儲運，2003，22(2)：16-19．

［3］陶文銓.數(shù)值傳熱學(xué)[M].2版.，西安：西安交通大學(xué)出版社，2001:332-409.

［4］ 吳子牛.計算流體力學(xué)[M]．北京：科學(xué)出版社，2001：70-76．

［5］ 趙海燕.不同工況下順序輸送混油的CFD模擬[J].化學(xué)工業(yè)與工程技術(shù)，2012,33（2）：54-57.

［6］ 常艷兵，王為民.蘭咸成品油管道順序輸送混油量分析與計算.當代化工，2011,40（11）：1195-1197.

［7］ 王岳、馬躍，張國軍,等.順序輸送粘度差對混油量影響分析[J]．遼寧石油化工大學(xué)學(xué)報，2012,32(4):37-40.．