勝利海上油田的累計開發時間長達26年，部分平臺采油井已進入開采后期，出砂嚴重，含水率高達90%以上，平臺間大量采用油氣水混輸，這導致油井至計量間的單井管線容易受到沖蝕，形成溝槽，彎頭位置尤為嚴重， 且隨著采出液含水率不斷上升以及化學驅影響，油井采出液攜帶砂的能力增強，管道的沖蝕風險加劇，管線壽命縮短，原油泄漏事故頻發。

固體顆粒沖蝕是結構壁面受到離散微粒沖擊產生損傷的一種現象。影響沖蝕的因素多達三十余種，其中比較重要的影響因素有顆粒速度，沖擊角，顆粒尺寸、形狀、材質，以及被沖擊表面的材料屬性等。含固體顆粒的多相流管道的沖蝕非常復雜，管道內的流型多變。這導致離散顆粒的運動更為復雜，而顆粒的運動及分布又直接決定了管道的沖蝕特性。

近年來，國內外學者對于管道的沖蝕進行了一系列研究，但關于油氣水砂三相流環境中的沖蝕研究尚不多見。為此，西安石油大學石油工程學院和中石化勝利油田分公司海洋采油廠的研究人員利用計算流體力學(CFD)數值分析方法，連接海底與海上平臺相連的輸送管道內的砂粒流動特性及其對管壁的腐蝕特性，考慮更為真實的油氣水三相黏性湍流態流體與砂粒之間的相互作用。

1 研究方法

針對海上油田的典型管道布置方式，選取沖蝕風險較大的海上平臺至海底的管段進行分析研究，管道三維模型如圖1所示。管道流動計算模型的網格剖分使用六面體單元進行劃分，最終得到的模型網格數量為670萬。計算網格模型如圖2所示。

圖1 管道結構示意

圖2 網格劃分示意

使用CFD程序對輸送管道內的油氣水砂多相流的物理現象進行求解，分別使用歐拉-歐拉多相流體系和歐拉-拉格朗日多相流體系求解連續相及離散相。在計算中，CFD多相流控制方程主要包含質量守恒和動量守恒方程，使用Edwards等給出的一種可應用于CFD的磨損預測模型預測顆粒腐蝕，該模型考慮了顆粒的碰撞速度、侵入角、形狀、直徑和質量流量等因素，其表達式為：

式中：V為磨損量；p為靶材流動應力；K為常數，經試驗驗證；a為沖擊角；m為顆粒質量；v為顆粒速度。

本工作并不需要考慮流體的化學屬性，因此近似將任意一相流體介質考慮為物性接近的純凈物。選取油田的典型輸量和含砂率進行分析。

2 試驗結果

流場計算結果

管道各處橫截面上的流速分布情況如圖3所示。可以發現，在重力和浮升力的作用下，油氣水混合物的高速區向下方偏移。進入立管彎頭時，在離心力的作用下在速度分布上體現為偏心程度更強。總體而言，速度分布在水平管段主要受重力作用，在彎頭附近主要受離心力作用。

圖3 截面流速分布云圖

油氣兩相的相含率分布情況如圖4所示?？梢园l現，各相在進入平臺水平管段的初始分布接近于均勻混合，隨著在水平管段的流動，在重力和浮升力的作用下氣相和油相逐漸向管頂聚集。進入立管彎頭后，在離心力的作用下氣相迅速向彎頭外側移動。在立管段，油氣兩相又逐漸恢復均勻混合分布，當由立管段進入海底水平段時氣相又向彎外側偏移。

（a）油相

（b）氣相

圖4 截面相分布云圖

沖蝕計算結果

管道兩個彎頭附近的含砂率分布如圖5所示?？梢园l現，砂粒在重力和自身慣性的作用下在平臺進入立管的彎頭附近時主要向彎心聚集。而在立管進入海底管的彎頭附近時，重力不再影響砂粒的偏轉，在重力的作用下砂粒向彎頭外側聚集?？梢姡瑳_蝕將主要發生在彎頭外側的壁面。

圖5 截面顆粒相分布云圖

圖6和圖7所示為不同條件下管道彎頭附近壁面當量沖蝕量，表示因沖蝕而產生的壁面單位面積材料損失量?？梢园l現，由于在彎頭附近離心力的作用，砂粒在彎頭外側的運動速度更高。而隨著輸液量的增加，砂粒對下彎管處管壁的沖蝕隨之增大，沖蝕范圍也有所增加，尤其是在輸液量達到1萬立方米/天時，上部彎管也開始出現砂粒沖蝕較明顯的區域。

（a） 0.2萬方/天  （b）0.5萬方/天

（c）1萬方/天 （d）2萬方/天

圖6 不同輸量下管道彎頭處腐蝕量分布圖

（a）0.005%    （b）0.01%

（c）0.02%     （d）0.04%

圖7 不同含砂率下管道彎頭處腐蝕量分布圖

圖8 不同條件下的管道最大腐蝕量

結合圖6~8可知，隨著輸液量的增加，沖蝕范圍有所增加，沖蝕率也隨之增長。相同輸液量的情況下，含砂率的增加雖然不會對沖蝕范圍造成影響，但卻會導致沖蝕率的增長。

3 結 論

(1) 使用CFD數值計算方法，選取適當的物理及數學模型能夠良好捕捉到管道中的流動形態以及相分布規律。選取領域內經過驗證的顆粒沖蝕模型能夠良好捕捉到管道中的腐蝕分布規律；

(2) 在慣性和連續流體介質共同作用下，砂粒堆積主要發生在管道彎頭的外側區域。而在本工作研究的海上油田低含砂量范圍下砂粒在管道內的總體沖蝕量仍然較低；

(3) 沖蝕率隨著輸液量的增加而增大，沖蝕范圍也同時增加。相同輸液量下，沖蝕率隨著含砂率的增加而增大，但對沖蝕范圍無明顯影響。