泄漏原油在河流表面漂移的數值模擬*

曹 權，王 吉，張 行，余東亮，蔣 毅

(1.中國石油大學(北京) 機械與儲運工程學院，北京 102249； 2.國家管網集團西南管道有限責任公司技術中心，四川 成都 610000)

0 引言

管道運輸具有運量大、密封性好、能耗小、長期運行穩定等優點。截至2017年底，中國油氣長輸管道總里程位列全球第3名，其中原油管道總里程達到3.2萬km[1]。我國地形復雜，西南地區地貌以高原與河流為主，長距離鋪設輸油管道常需要穿越河流。穿河管道長期受到河床泥沙沖蝕，河水沖刷與腐蝕，極易發生泄漏，泄漏后原油進入河流，在河水的攜帶下會快速順流而下，短時間內造成大范圍的環境污染[2]。2003年8月5日，停泊在吳涇電廠的中海“長陽”輪被1艘小船碰撞，油箱受損導致85 t燃油泄漏，在江面造成近15萬m2的污染[3]。如何快速阻止泄漏原油順河而下的漂移，設置圍油欄回收泄漏原油，對于相關部門做出應急決策，減輕事故風險具有重要意義。而指導應急決策的核心環節就是對原油順河漂移過程的準確把握。

關于水面溢油，國內外學者開展相關研究，大多數研究針對海洋表面展開。針對海洋溢油，已開發出一系列溢油模擬系統，最新1代溢油軟件如GNOME[4]，OILMAP[5]，MIKE21/3SA[6]等，基于油膜擴展理論[7]結合油粒子模型[8]，采用隨機游走法和拉格朗日粒子追蹤法預測油粒子時空分布，同時加入海面溢油的蒸發、溶解、乳化、沉降等過程[9]。杜文強[10]采用泊松圓盤采樣算法，在油膜達到最大面積時將油膜離散為油粒子，解決油膜分散問題。現有研究多集中于海洋，港口等區域，對淺窄河道溢油研究極少。部分學者利用MIKE軟件中的水動力模型及溢油模塊對河流溢油進行了二維數值模擬，得到河道溢油的漂移軌跡及濃度分布等信息[11-12]。

隨著計算流體力學技術發展，眾多學者利用計算流體力學方法建立一系列水下原油泄漏運移模型。Yapa等[13]、Johansen[14]基于拉格朗日積分法建立水下原油噴射三維模型。Tkalich[15]采用歐拉法對原油海底擴散過程進行數值模擬，并預測溢出點油膜厚度。曹學文等[16]建立水下二維水箱模型，采用多相流VOF(Volume of Fluid)模型對水下原油管道泄漏進行二維非穩態數值模擬，研究不同漏油速率、流速等因素對水下原油擴散的影響。Cloete等[17]基于離散相DPM模型(Discrete Phase Model)，結合VOF模型，對水下天然氣管道泄漏擴散行為進行模擬。

本文選用國內某一原油管道穿越的真實河流模擬原油在河面的漂移過程，采用Mixture(混合)多相流模型，計算獲得原油泄漏后的漂移范圍，為管道泄漏后制定及時準確的攔截方案提供理論依據。

1 數值計算模型

研究對象選取中國西南地區某條河(記為S江)，有1段原油管道穿河而過，穿越段地處平原，河流落差小，河水較淺，可以假設穿越段河水深度不變。原油在河流表面漂移主要受到河水流動的影響，忽略蒸發、溶解、乳化和沉降等因素的作用。水下原油管道發生泄漏后，上浮至河流表面的過程占據較少的時間，本文忽略原油上浮過程，從原油到達水面形成油膜之后開始考慮。由于原油被河水裹挾流動，原油在河面的漂移軌跡可近似為泄漏點位置的河水表面流體微團的流動軌跡，因此本文將在泄漏源位置處添加被標記的水分子質量源項模擬原油泄漏到達河流表面的過程，通過記錄這些被標記的水分子在河流表面的移動得到被主流河水裹挾的原油在河面擴散漂移的過程。

1.1 多相流模型

多相流模型Mixture允許相間滲透，與標記水分子進入主流過程原理一致。河水與泄漏源位置產生的被標記的水均視為不可壓縮流體，混合模型的連續性方程如式(1)～(3)所示：

(1)

(2)

(3)

式中：t為時間，s；ρm為混合密度，kg/m3；為哈密頓算子；為質量平均流速，m/s；k為相數；αk為第k相的體積分數；ρk為第k相密度，為第k相滑移速度，m/s。

混合相的動量方程如式(4)～(5)所示：

(4)

(5)

第k相的體積分數計算公式如式(6)所示：

(6)

將單相的標準k-ε模型擴展至多相流模型，將單相流的壓力速度耦合SIMPLEC算法擴展至多相流中，利用Fluent求解器對方程進行數值求解。

1.2 邊界條件

河流入口流速計算公式如式(7)所示：

(7)

式中：v為河流入口流速，m/s；Q為河水流量，m3/s；A為河流入口截面積，m2。

管道泄漏質量流速計算公式如式(8)～(10)所示：

(8)

(9)

(10)

式中：G為泄漏質量流速，kg/(m2·s)；As為單個泄漏質量源面積，m2；vl為管道泄漏質量流量，kg/s；Al為管道裂口面積，m2；Cv為流體泄漏系數，為經驗常數，Cv=0.60；g為重力加速度，取9.81 m/s2；Δh為泄漏點處流體壓力水頭，m；p為輸送管道內介質壓力，Pa；p0為環境壓力，Pa；ρ為泄漏流體密度，kg/m3；v0為管內介質流速，m/s。

河岸為無滑移固體邊界，河道入口、出口分別為速度入口邊界、自由出流邊界。泄漏口定義為泄漏質量源項。

1.3 物理模型

穿越段河流二維物理模型如圖1所示。在得到原始河流地圖后，對原始圖像進行截取和裁剪，只保留河道信息，得到僅包括河道區域的目標圖像，河道信息地圖如圖1(a)所示。為獲得河流邊緣坐標信息，對目標圖像進行二值化處理，二值目標圖像如圖1(b)所示。二值目標圖像邊緣坐標利用Canny算子提取，為重建連續通道邊界，將這些離散點連接起來形成NURBS曲線。由于原始圖片分辨率較高，直接重建的河道細節過多，存在尖銳邊緣。為平滑河道形狀，在重建過程中，沿邊界均勻間隔選取1/10的離散數據點，以較少的數據表示主要幾何特征，平滑插值獲得可用于后續數值計算的河道形狀。河道二維模型如圖1(c)所示。

圖1 穿越段河流二維物理模型Fig.1 Two-dimensional physical model of crossing section river

1.4 網格生成

自西向東在管道穿越位置上取10個泄漏點，即在水面等距建立10個泄漏質量源項，分別編號為1～10，代表原油溢出區域。對物理模型進行結構化網格劃分，對河岸附近及泄漏口附近進行網格加密處理，最終生成網格總量277 952。建立的泄漏口模型及附近的網格如圖2所示。

圖2 泄漏口模型Fig.2 Model of leakage holes

1.5 模擬工況

河流、原油及管道參數見表1。

穿越河流管道發生泄漏位置不同，導致河面溢油區域位置不同，考慮河流表面河水速度分布以及河岸對原油漂移的影響，將泄漏位置劃分為河岸邊緣、近岸處及河流中部3個區域。S江全年平均流量為250.1 m3/s，計算段河流入口平均流速為2 m/s，考慮S江在枯水期及豐水區河水流速不同，將河流入口流速分別設置為1.58，2.58，3.58 m/s。水下原油管道工作環境惡劣，本文泄漏量按照管道發生完全斷裂計算得到河流表面泄漏質量流速范圍。模擬工況設置見表2。

2 計算結果分析

2.1 泄漏位置對原油漂移過程的影響

《地表水環境質量標準》(GB 3838—2002)[18]規定：地表水中石油類污染物體積分數大于5.638×10-8為污染區。河道邊緣位置發生泄漏后，原油在河流表面的漂移過程如圖3所示。原油在河流邊緣處發生泄漏時，河流變窄與河流曲折對漂移過程影響不大，原油會緊貼河岸向下游漂移。河流擴張段對原油漂移影響較大，在西側(泄漏口1)發生泄漏后，原油漂移至河道擴張段，會在擴張段西側形成漩渦。與西側發生泄漏相比，在東側(泄漏口10)發生泄漏后，原油漂移進入河流擴張段未出現旋渦，但會在擴張段東側末端近岸處出現堆積。經過擴張段后，原油漂移橫向擴散范圍變寬，基本充斥整個河道繼續向下游漂移。

表1 河流、原油及管道參數Table 1 Parameters of river, crude oil and pipeline

河流近岸處發生泄漏時原油在河流表面的漂移過程如圖4所示。原油向下游漂移過程中會逐漸向河岸靠近，最終緊貼河岸。在西側近岸處(泄漏口3)發生泄漏時，原油漂移進入河流擴張段，會出現旋渦，在河面大量堆積，幾乎充滿整個擴張段，東側近岸處(泄漏口9)發生泄漏時，原油漂移進入河流擴張段的角度偏向河流中部，會在擴張段出現微小旋渦。經過擴張段后，原油會充斥整個河流表面繼續向下游漂移。

表2 模擬工況Table 2 Simulation conditions

圖3 邊緣處泄漏漂移過程Fig.3 Drift process of leaking at edges

圖4 近岸處泄漏漂移過程Fig.4 Drift process of leaking at nearshore

河流中部(泄漏口5，泄漏口7)發生泄漏時原油在河流表面的漂移過程如圖5所示。分析可知：原油在河流中部發生泄漏后，河流變窄與河流曲折對漂移過程影響不大，原油基本沿著河流中部向下游漂移，漂移至河流擴張段，原油會充斥整個擴張段。此外泄漏位置越靠近河流中部，原油漂移進入擴張段形成的污染面積越大。

圖5 中部泄漏漂移過程Fig.5 Drift process of leaking at middle

不同泄漏位置發生泄漏后原油在計算域中的污染面積隨漂移時間變化如圖6所示。泄漏初期，原油在河流表面污染面積與漂移時間變化曲線呈近似線性關系，在約0.85 h時，原油漂移至河流擴張段，河水流速與方向變化較大，原油在河流表面污染面積的增長速率變大。在約2.3 h時，原油漂移到計算域中的河流出口，此后計算域內污染面積基本穩定下來，隨時間變化不大。此外，泄漏位置越靠近河流中部，河流污染面積增長速率越大，最終造成的污染面積越大。

圖6 不同泄漏位置原油漂移污染面積變化Fig.6 Variation of pollution area by crude oil drift under different leakage locations

2.2 河流入口流速對原油漂移過程的影響

圖7 不同河流入口流速下原油漂移過程(2.78 h)Fig.7 Drift process of crude oil under different river inlet flow velocities (2.78 h)

發生泄漏2.78 h后，不同河流入口流速下原油在河流表面漂移過程如圖7所示。原油在河流表面漂移造成的污染面積與漂移時間變化如圖8所示。泄漏位置越靠近中部，造成污染面積越大。對不同泄漏口，河流入口流速一致時，漂移至河流出口所需時間基本相同，隨著河流入口流速增加，原油在河流表面的漂移速度會明顯變大，在河流入口流速為1.58 m/s時，原油在約2.3 h時漂移至河流出口，在河流入口流速為2.58 m/s時，原油在約1.7 h時漂移至河流出口，在河流入口流速為3.58 m/s時，原油在約1.3 h時漂移至河流出口。對同一泄漏口，原油漂移至河流擴張段時，河水流速的變大導致油膜變薄，污染面積隨漂移時間增大的速率變大，但原油漂移至河流出口后，最終在計算域內河流表面形成的污染面積會趨于一致，說明河流入口流速對原油漂移在計算域內最終形成的污染面積影響不大。

圖8 不同河流入口流速下原油漂移污染面積變化Fig.8 Variation of pollution area by crude oil drift under different river inlet velocities

2.3 泄漏質量流速對原油漂移過程的影響

發生泄漏2.78 h后，不同泄漏質量流速下原油在河流表面漂移過程如圖9所示。原油在河流表面漂移造成的污染面積與漂移時間變化如圖10所示。原油均在1.3 h左右漂移至計算區域出口，泄漏位置越靠近河流中部，造成污染面積越大。對同一泄漏口，隨著泄漏質量流速增加，原油污染面積略有增加，油膜變厚。

圖9 不同泄漏質量流速下原油漂移過程(2.78 h)Fig.9 Drift process of crude oil under different leakage mass flow rates (2.78 h)

圖10 不同泄漏質量流速下原油漂移污染面積變化Fig.10 Variation of pollution area by crude oil drift under different leakage mass flow rates

3 結論

1)采用CFD軟件建立原油河面漂移預測模型，對不同泄漏位置、不同河水入口流速、不同泄漏質量流速情況下發生泄漏后原油在河流表面的擴散漂移過程進行模擬分析，從而得到泄漏原油的漂移范圍及污染面積等重要信息。

2)穿越河流管道泄漏位置對原油在河流表面漂移過程影響較大。若河流形狀未出現擴張，在河流邊緣及近岸處發生泄漏時，漂移過程基本貼近河岸；在河流中部發生泄漏時，基本不會漂移至河道邊緣；突遇河流形狀出現擴張，原油會在擴張段出現漩渦，泄漏位置越靠近河流中部，在擴張段形成的污染面積越大。

3)原油在河流表面漂移速度主要受到河水入口流速影響，與泄漏位置與泄漏質量流量關系不大，河流入口流速越大，漂移速度越快。泄漏質量流速大小對原油在河流表面污染范圍影響不大。