安徽二期成品油管道順序輸送的運行分析

2018-05-22 08:00:49李鑫偉

石油庫與加油站 2018年1期

李鑫偉

〔中國石化銷售有限公司華中分公司湖北武漢 430023〕

1 安徽二期管道概況

中國石化安徽二期成品油管道是將安慶石化的成品油輸送到合肥油庫、淮南油庫、阜陽油庫和蚌埠油庫，以滿足管道沿線合肥、淮南、阜陽和蚌埠等地區的成品油市場需求。管道全長519.4 km，管道設計壓力9.5 MPa，管道入口設計輸量為410×104t/a。全線采用密閉工藝順序輸送92號汽油、95號汽油、0號車用柴油和0號普通柴油兩大類4個品種。各段間距、管徑、管容及設計壓力見表1。

表1 安徽二期管道沿線管段一覽表

管道沿線設置18座線路截斷閥室(其中4座為遠控閥室)、9座陰極保護站和10座光通信站。全線自動控制采用SCADA系統，沿線各站場與調度中心采用光纜傳輸方式進行數據交換，調度中心通過各站的站控系統對全線進行數據采集及控制，并負責編排批次、批量計劃、輸送方案、跟蹤批量、泄漏檢測和優化運行等。輸送流程為：

本文以安徽二期成品油管道為例，依據成品油管道水力特性和順序輸送理論，對該管道各管段的臨界雷諾數、臨界流量、水力摩阻系數等數據進行計算和分析，并總結出適合各管段的水力計算公式，為配泵和高效節能運行提供理論依據。

2 臨界雷諾數Re、混油臨界雷諾數Rec和臨界流量Qc的計算

2.1 臨界雷諾數Re的計算

管線和所輸成品油的流動狀態用臨界雷諾數表示式(1)：

(1)

式中：Q——油品在管路中的體積流量，m3/s；

d——管道直徑，m；

ν——油品的運動粘度，m2/s。

管壁相對當量粗糙度ε計算，如式(2)：

(2)

式中：e——管壁的絕對當量粗糙度，m。

光滑區與混合摩擦區的臨界雷諾數計算，見式(3)：

(3)

混合摩擦區與粗糙區臨界的雷諾數計算，見式(4)：

(4)

2.2 混油臨界雷諾數Rec的計算

采用奧斯汀和柏爾弗萊由雷諾數Re與無因

次量C2/Ld的關系總結出的經驗公式(5)：

(5)

2.3 臨界流量Qc的計算

臨界流量的計算公式(6)如下：

Qc=3 600πdvRec/4

(6)

式中：QC——臨界流量，m3/h；

d——管道內徑，m；

ν——油品的運動粘度，m2/s；

Rec——混油臨界雷諾數。

式(5)、(6)的計算結果見表2。當雷諾數ReRe>Re1時流態處于混合摩擦區，當Re>Re2時流態處于粗糙區，此時混油長度則隨雷諾數的減小增長較為緩慢。故參考表2計算結果，順序輸送過程中當管道存在混油頭時，應保證各管段的混油界面雷諾數Re>Rec，最小流量Qmin>QC，以防因管輸流量小而產生大量混油。

表2 混油臨界雷諾數Rec和臨界流量Qc的計算結果

3 沿程水力計算

管道輸油過程中壓力能的消耗主要包括兩部分，一是用于克服地形高差所需的位能；二是克服油品沿管路流動過程中的摩擦及撞擊產生的能量損失轉換成的液柱高度，通常稱為摩阻損失。長輸管道的摩阻損失主要由通過直管段所產生的沿程摩阻和通過各閥件、管件所產生的局部摩阻組成。而長輸管道站間管路的摩阻損失主要是沿程摩阻，局部摩阻只占1 %～2 %。

因此，為了計算各控制點在不同工況下的壓力控制值，首先需用達西公式計算推導出列賓宗公式，再由列賓宗公式計算出沿程摩阻大小，最后由管道全線壓力供需平衡關系式計算出壓力控制值。

水力摩阻系數計算，如式(7)：

(7)

式中：A、m根據雷諾試驗可知水力摩阻系數λ與不同流態下流體的雷諾數有關，因此得到不同流態下與雷諾數相關的常數系數A、m。

達西公式如式(8)：

(8)

式中：V—油品在管路中的流速，m/s；

hl—管路的沿程摩阻損失，m；

L—管道長度，m；

d—管道直徑，m；

g—重力加速度，取9.81 m/s2。

(9)

式中：Q——油品在管路中的體積流量，m3/s；

d——管道直徑，m；

ν——油品的運動粘度，m2/s；

L——管道長度，m；

A、m——不同流態下與雷諾數相關的常數系數；

β——列賓宗公式(9)含A、m的簡化替代量。

各流態區的A、m、β值及沿程摩阻計算式見表3。

得出各流態下流體對應的沿程摩阻計算公式后，將其代入管道全線壓力供需平衡關系式中即可求得各工況下各管段沿程摩阻損失的計算值。通過實際值與計算值比較可以得到誤差系數，具體記錄數據見表4。

假設一條管道有N座泵站，且泵站特性相同，全線管徑相同，無分支，首站進站壓頭和各站內摩阻均為常量，可寫出全線的壓力供需平衡關系式(10)如下

Hs1+NHY-N×hm=HI+(Z2-Z1)+Hs2

(10)

式中：Q——全線工作流量，m3/s；

N——全線泵站數；

Hl——管道沿程摩阻值，m液柱；

Hs1——管道首站進站壓頭，m液柱；

Hs2——管道終點剩余壓頭，m液柱；

HY——泵站揚程，m；

Z1、Z2——管道起、終點高程，m；

hm——每個泵站的站內損失，m液柱。

表3 不同流態時的A、m、β值及列賓宗公式

表4 安徽二期各管段不同工況下運行的記錄數據

通過采集數據將其代入公式(10)即可得到實際各管段在不同工況下沿程摩阻損失值。再通過式(7)、(8)、(9)及表3對應公式計算出相應的沿程摩阻損失值。由計算值和實際值比較得到的各管段在不同工況下的誤差系數見表5。由此值可以更加準確地計算各管段的摩阻損失，為配泵和高效節能運行提供了理論依據。各段管道在不同工況下的誤差系數見表5。

表5 各段管道在不同工況下的誤差系數

由表5可知，安徽二期管道各管段在單輸汽油時計算摩阻損失值與實際摩阻損失值相差過大，因此在計算時不可以用列賓宗公式。由于安合二期管道各管段在單輸汽油時流態均由水力光滑區進入混合摩擦區，可見列賓宗公式在混合摩擦區計算時誤差較大。

Moody圖是由自然粗糙管實驗得出，Cole-brook-White(以下簡稱C-W)公式(11)為：

(11)

該公式綜合了光滑管路普朗特公式和粗糙管尼古拉茲公式得出的，又經實驗修正。因此C-W公式被公認為計算水力摩阻系數λ的精確性標準。遺憾的是，公式(11)為一隱函數，必須用迭代方法才能求解，不便于工程應用。下面介紹幾個混合摩擦區的摩阻系數公式均為顯函數形式，以此來確定安徽二期管道各管段混合摩擦區計算公式。

Jain公式(12)：

(12)

Swamee&Jain公式(13):

(13)

Haaland公式(14)：

(14)

Manadilli公式(15)：

(15)

將表4安徽二期各管段單輸汽油數據代入式(12)、(13)、(14)和(15)中得到表6、表7、表8、表9。