水力壓裂管柱彎管處的固體顆粒沖蝕研究

嚴(yán)淳鳀

水力壓裂管柱彎管處的固體顆粒沖蝕研究

嚴(yán)淳鳀

（長(zhǎng)江大學(xué)機(jī)械工程學(xué)院，湖北 荊州 434023）

為了研究固/液兩相流體對(duì)壓裂管柱彎管處的沖蝕問(wèn)題，運(yùn)用DPM沖蝕模型對(duì)彎管處在不同固相顆粒體積分?jǐn)?shù)條件下的沖蝕進(jìn)行計(jì)算與分析，追蹤固體顆粒粒子路徑，研究結(jié)果表明：彎管內(nèi)壁管徑最大區(qū)域最易產(chǎn)生沖蝕；在不同固相顆粒體積分?jǐn)?shù)條件下，彎管最大沖蝕率隨著固體顆粒含量的增加而增加，且固相顆粒體積分?jǐn)?shù)與彎管的沖蝕率呈線性變化趨勢(shì)。

壓裂彎管；固體顆粒；沖蝕；液/固兩相流

在非常規(guī)油氣井開(kāi)采中，水力壓裂技術(shù)是提高采收率和可采儲(chǔ)量的重要手段之一；而且水力壓裂在低滲透油氣田中是目前最有效的開(kāi)發(fā)技術(shù)。目前，水力壓裂工藝向大的施工排量和大的加砂量趨勢(shì)發(fā)展，這就對(duì)管柱和工具的耐磨性和安全性提出了新的要求[1]。

水力壓裂就是利用壓裂車將壓裂液以高壓力和高流速?gòu)木沧⑷刖诘倪^(guò)程，在壓裂過(guò)程中，油氣層產(chǎn)生裂縫后還會(huì)注入含有支撐劑的混砂液；在水力壓裂施工工序中，注入支撐劑、泵入頂替液和反洗都將會(huì)產(chǎn)生固體顆粒，凡是與運(yùn)動(dòng)流體產(chǎn)生直接接觸的所有類型的設(shè)備，都將會(huì)產(chǎn)生沖蝕磨損[2]；因此，本文對(duì)井下管柱的彎管處進(jìn)行沖蝕研究。

管道產(chǎn)生沖蝕是一個(gè)極為復(fù)雜的過(guò)程，當(dāng)流體運(yùn)動(dòng)為湍流時(shí)，其運(yùn)動(dòng)狀態(tài)非常無(wú)序，流體可從各個(gè)方向?qū)饘俦诿孢M(jìn)行沖擊，金屬壁面受沖蝕較為嚴(yán)重，因此，在管道的變徑處和彎曲段壁面受到的沖蝕更加明顯。而且，管道沖蝕與流體內(nèi)固體顆粒含量、管線材料、管道幾何結(jié)構(gòu)、管道彎曲處彎曲角度、流體沖擊角度等許多因素都有一定的關(guān)系[3]。

文中利用SolidWorks對(duì)彎管模型進(jìn)行構(gòu)建，再將模型導(dǎo)入Ansys的計(jì)算流體模型Fluent中，開(kāi)展固/液兩相流體流動(dòng)特性仿真模擬計(jì)算，研究混砂液在彎管處的內(nèi)部流場(chǎng)分布規(guī)律，對(duì)不同固相顆粒體積分?jǐn)?shù)條件下彎管的沖蝕進(jìn)行計(jì)算與分析，得出了相關(guān)結(jié)論。

1 模型建立

1.1 沖蝕磨損模型

影響液/固兩相流沖蝕的因素有很多，Mixture- discrete phase model（DPM）分析法模型不僅考慮了沖擊速度和角度的影響，還將管線的材料管線的壁厚、顆粒粒徑、顆粒形狀和顆粒材料性質(zhì)一起考慮進(jìn)去[4]，文中采用此模型完成沖蝕計(jì)算[5]：

式中：R—彎管沖蝕速率，kg·（m2·s）-1；

（d）—離散相顆粒直徑函數(shù)；

N—顆粒數(shù)目；

m—固體顆粒質(zhì)量流量，kg·s-1；

—固體顆粒與彎管壁面的碰撞角，（°）；

A—彎管壁面計(jì)算單元面積，m2；

（）—侵入角函數(shù)。

1.2 幾何模型及網(wǎng)格

模型如圖1所示，由上游直管入口段、下游直管出口段和彎管段組成，管道采用3 1/16″彎頭規(guī)格，彎管角度為90°，曲率半徑=3D，管徑=100 mm，壁厚=15 mm，==300 mm，含顆粒的流體從入流端以30 m·s-1的速度垂直進(jìn)入；在Ansys中建立流道模型后，并進(jìn)行網(wǎng)格劃分，采用六面體網(wǎng)格單元，得到網(wǎng)格數(shù)為1 441 872，如圖2所示。

圖1 彎管幾何模型

圖2 彎管模型網(wǎng)格

2 計(jì)算結(jié)果及比較

2.1 計(jì)算前處理

入口邊界條件為速度入口，設(shè)置固相顆粒和沖蝕磨損模型，固相顆粒的入射速度和流體相同；彎管內(nèi)連續(xù)相介質(zhì)為壓裂液，分散相為砂粒；出口邊界采用outflow。

固相顆粒體積分?jǐn)?shù)分別設(shè)置為1%、3%、5%、7%、9%、11%；入口速度為30 m·s-1；砂粒顆粒視密度為=1 500 kg·m-3，顆粒粒徑為= 0.35 mm；沖擊角函數(shù)定義為線性數(shù)值，如表2所示。

表1 沖擊角函數(shù)參數(shù)設(shè)置

2.2 仿真結(jié)果分析

利用上述仿真方法在固相顆粒體積分?jǐn)?shù)分別為1%、3%、5%、7%、9%、11%時(shí)進(jìn)行模擬計(jì)算，得到固體顆粒體積分?jǐn)?shù)對(duì)彎管沖蝕的影響。經(jīng)沖蝕仿真計(jì)算后，觀察固體顆粒通過(guò)彎管時(shí)的運(yùn)動(dòng)軌跡，如圖3所示。固體顆粒在0 s時(shí)刻，從下方入口處進(jìn)入彎管中，在彎管彎曲處于管壁發(fā)生碰撞，隨后在0.052 6 s時(shí)刻從出口處流出。

在Fluent軟件中計(jì)算不同固體顆粒體積分?jǐn)?shù)下的彎管沖蝕率，其云圖如圖4所示。不同固體顆粒體積分?jǐn)?shù)與最大沖蝕率的曲線圖如圖5所示。

圖3 固體顆粒運(yùn)動(dòng)軌跡圖

圖4 不同固體顆粒含量下彎管沖蝕云圖

圖5 最大沖蝕率-固相顆粒體積分?jǐn)?shù)曲線圖

由圖5相關(guān)結(jié)果可知，隨著固相顆粒含量的增加，彎管內(nèi)壁面最大沖蝕速率呈近似線性增長(zhǎng)。壓裂液質(zhì)量濃度增加，使得在單位時(shí)間、單位體積內(nèi)壓裂液中固體顆粒增加，則導(dǎo)致在單位時(shí)間、單位面積上與彎管內(nèi)壁的碰撞次數(shù)增多，從而使得最大沖蝕速率增大。

3 結(jié) 論

本文主要分析水力壓裂管柱彎管處在不同的固體顆粒含量下對(duì)管柱的沖蝕影響。通過(guò)Fluent的仿真模擬，可以得到以下結(jié)論：

1）根據(jù)彎管沖蝕云圖可知，易產(chǎn)生沖蝕區(qū)域在彎管內(nèi)壁管徑最大區(qū)域，在上游直管段和下游直管段產(chǎn)生沖蝕較小，應(yīng)重點(diǎn)關(guān)注易產(chǎn)生沖蝕區(qū)域，以免造成重大安全事故。

2）彎管沖蝕率隨著固體顆粒體積分?jǐn)?shù)的增加而增加，固相顆粒體積分?jǐn)?shù)在1%、3%、5%、7%、9%、11%時(shí)，最大沖蝕率分別為0.000 064 7、0.000 065 3、0.000 352、0.001 07、0.001 67、0.002 01 kg·m-2·s-1，且固相顆粒體積分?jǐn)?shù)與彎管的沖蝕呈線性變化趨勢(shì)。

［1］路鑫, 常宗瑜, 趙林, 等. 大口徑海洋輸氣管道彎管處的固體顆粒沖蝕效果分析[J]. 機(jī)械設(shè)計(jì), 2019, 36 (S1): 34-39.

［2］周蘭, 張紅, 陳文康, 等. 頁(yè)巖氣壓裂管匯彎頭的沖蝕磨損影響分析[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù), 2020, 16 (10): 53-58.

［3］董春朋. 分層壓裂工藝下管柱及工具易沖蝕部件研究[D]. 大慶：東北石油大學(xué), 2015.

［4］王健剛, 孫巧雷, 嚴(yán)淳鳀, 等. 測(cè)試地面流程彎管沖蝕磨損的影響研究[J]. 石油機(jī)械, 2021, 49 (1): 88-94.

［5］ZHANG Y, REUTERFORS E P, MCLAURY B S, et al. Comparison of computed and measured particle velocities and erosion in water and air flows[J]., 2006, 263 (1)：330-338.

Research on Solid Particle Scouring in Bends of Hydraulic Fracturing Pipe Columns

（School of Mechanical Engineering, Yangtze University, Jingzhou Hubei 434023, China）

In order to investigate the erosion of fracture tubes at bends of fracture tube columns by solid/liquid two-phase fluids, the erosion model of DPM was used to calculate and analyze the erosion of bends at different solid-phase particle volume fractions and to trace the path of solid particles. The results showed that, the area with the largest pipe diameter on the inner wall of the elbow was most likely to be eroded, the maximum erosion rate of the elbow increased with the increase of the solid particle content, and the volume fraction of solid particles and the erosion rate of the elbow showed a linear change trend.

Fracture bends; Solid particles; Erosion; Liquid/solid two-phase flow

2020-11-30

嚴(yán)淳鯷（1994-），女，湖南常德人，碩士研究生，研究方向：機(jī)械系統(tǒng)仿真與診斷技術(shù)。

TQ050.9

A

1004-0935（2021）04-0526-03