陳珊珊 時(shí)培忠 馮義 蓋靖安 樓一珊 翟曉鵬 史寶成

(1.長(zhǎng)江大學(xué)石油工程學(xué)院 2.新疆油田分公司工程技術(shù)研究院3. 中石油吐哈油田分公司工程技術(shù)研究院 4. 中石油西部鉆探工程有限公司國(guó)際工程公司)

0 引 言

防砂失效會(huì)導(dǎo)致油井出砂、產(chǎn)液低，大幅降低油田的生產(chǎn)效益。篩網(wǎng)故障通常發(fā)生在完井與油井生產(chǎn)的過(guò)程中，但是故障原因卻難以準(zhǔn)確確定，一般是堵塞與沖蝕共同作用下導(dǎo)致的失效[1]。因此，解決因沖蝕磨損引發(fā)的防砂失效，準(zhǔn)確預(yù)測(cè)防砂壽命成為保證油田穩(wěn)定生產(chǎn)的關(guān)鍵舉措。

當(dāng)前對(duì)于防砂篩管的沖蝕研究多集中在篩管的性能評(píng)價(jià)上。S.MONDAL等[2-3]采用試驗(yàn)研究與數(shù)值模擬相結(jié)合的方式來(lái)評(píng)估篩網(wǎng)性能，優(yōu)選篩網(wǎng)。G.GILLESPIE等[4]以Troika油田為例，開(kāi)發(fā)出防砂篩管沖蝕模型，用于驗(yàn)證篩管的實(shí)際應(yīng)用情況，但若要使用該模型預(yù)測(cè)篩管壽命，則仍需進(jìn)行更多的測(cè)試評(píng)估。J.A.CAMERON等[5]以基本沖蝕理論為基礎(chǔ)，考慮砂粒大小和沖擊角度等因素，通過(guò)試驗(yàn)結(jié)果獲得新的沖蝕數(shù)學(xué)模型，并以2%質(zhì)量損失極限作為防砂失效標(biāo)準(zhǔn)，推導(dǎo)得到篩管壽命預(yù)測(cè)公式。但現(xiàn)有研究并沒(méi)有考慮井中使用的特定篩管系統(tǒng)，無(wú)法給出現(xiàn)有流動(dòng)條件下預(yù)期油井壽命的定量指標(biāo)，或者油井壽命會(huì)隨著流量的增加而縮短等問(wèn)題。

針對(duì)防砂篩管沖蝕壽命難以預(yù)測(cè)的問(wèn)題，本文以金屬網(wǎng)布復(fù)合篩管[6]為例，建立了全結(jié)構(gòu)篩管模型，基于計(jì)算流體力學(xué)(CFD)進(jìn)行數(shù)值模擬，根據(jù)沖蝕率的物理意義預(yù)測(cè)沖蝕壽命，以期為其現(xiàn)場(chǎng)應(yīng)用提供理論指導(dǎo)。

1 仿真模擬計(jì)算模型

油井生產(chǎn)過(guò)程中夾帶的砂粒會(huì)對(duì)防砂篩管造成沖蝕磨損。防砂篩管沖蝕過(guò)程屬于固液兩相流作用下的沖蝕問(wèn)題。由于涉及到的固體顆粒體積分?jǐn)?shù)遠(yuǎn)小于10%，計(jì)算區(qū)域?yàn)橄∠嗫臻g，所以對(duì)固體顆粒采用離散相模型進(jìn)行計(jì)算。相關(guān)計(jì)算模型如下。

1.1 流體控制方程

流體控制方程包括連續(xù)性方程和動(dòng)量方程。連續(xù)性方程為：

(1)

動(dòng)量方程為：

(2)

式中：ρ為流體密度,kg/m3；u為流體速度,m/s；x為坐標(biāo)軸，i，j為空間方位；p為壓力,Pa；υ為連續(xù)相的運(yùn)動(dòng)黏度,m2/s；fi為質(zhì)量力,N；f為流場(chǎng)中受到的其他力,N。

標(biāo)準(zhǔn)k-ε模型[7]是一種半經(jīng)驗(yàn)?zāi)Ｐ?，基于湍?dòng)能k和耗散率ε二階相關(guān)量的輸送方程，具有精度高、適用范圍廣的優(yōu)點(diǎn)，故采用標(biāo)準(zhǔn)k-ε模型。其控制方程為：

(3)

(4)

式中：k為湍流動(dòng)能，J；μt為湍流黏度，Pa·s；σk、σε分別為k方程和ε方程的湍流Prandtl數(shù)，σk取1.0，σε取1.3；Gk為層流速度梯度產(chǎn)生的湍流動(dòng)能，J；Gb為浮力產(chǎn)生的湍流能項(xiàng)，J；ε為湍流耗散率，W/m3；C1ε、C2ε、C3ε分別為經(jīng)驗(yàn)常數(shù)。

1.2 砂粒運(yùn)動(dòng)方程

砂粒在流體中主要受到流體作用力與顆粒碰撞作用力。顆粒碰撞作用力主要為顆粒間的碰撞作用力和顆粒與壁面的碰撞作用力。因?yàn)殡x散相體積分?jǐn)?shù)遠(yuǎn)小于10%，計(jì)算區(qū)域?yàn)橄∠嗫臻g，所以可忽略顆粒間的相互作用。固體顆粒在兩相流中受到虛擬質(zhì)量力、壓力梯度力、布朗力和重力等作用力，具體表達(dá)式如下[8]：

(5)

(6)

(7)

(8)

式中：up為砂粒速度，m/s；ρp為砂粒密度，kg/m3；d為砂粒直徑，m；Fy為y方向的其他作用力，包括虛擬質(zhì)量力、壓力梯度力、布朗力和Saffman升力，N；CD為曳力系數(shù)；Rep為相對(duì)Reynolds數(shù)；α1、α2和α3為常數(shù)；FD為流體黏性阻力，N；gy為重力加速度，m/s2。

離散相模型可以忽略顆粒間的相互作用。顆粒與壁面材料發(fā)生碰撞時(shí)，將碰撞前、后速度分量的比值定義為恢復(fù)系數(shù)。采用G.GRANT和W.TABAKOFF提出的恢復(fù)系數(shù)[9]進(jìn)行數(shù)值計(jì)算，其表達(dá)式為：

(9)

式中：εN為法相恢復(fù)系數(shù)；εT為切向恢復(fù)系數(shù)；α為顆粒軌跡與壁面的沖擊角度，(°)。

1.3 沖蝕模型

Fluent所使用的DPM沖蝕模型[10]為：

(10)

式中：Re為沖蝕率，kg/(m2·s)；C(dp)為顆粒直徑函數(shù)，C(dp)=1.8×10-9；f(α)為沖蝕角度函數(shù)，f=1；v為固體顆粒的相對(duì)速度；b(v)為顆粒相對(duì)速度函數(shù)；A為沖蝕面積，m2；mb為質(zhì)量流率，kg/s。

2 篩管沖蝕模擬

2.1 篩管全結(jié)構(gòu)建模

通過(guò)對(duì)實(shí)際篩管進(jìn)行切片分析，研究金屬網(wǎng)布復(fù)合篩管的結(jié)構(gòu)布局。金屬網(wǎng)布復(fù)合篩管由保護(hù)殼、多層篩網(wǎng)和基管組成。根據(jù)現(xiàn)場(chǎng)實(shí)際沖蝕破壞結(jié)果，主要沖蝕破壞位置為第一層篩網(wǎng)，考慮到模型計(jì)算工作量，建立了保護(hù)殼、雙層篩網(wǎng)和基管的三維模型。金屬網(wǎng)布的編織方法為麻花狀，通過(guò)測(cè)量計(jì)算，最終確定金屬網(wǎng)布的單元尺寸。

在建立金屬網(wǎng)布復(fù)合篩管三維模型的過(guò)程中，發(fā)現(xiàn)由于篩網(wǎng)的編織結(jié)構(gòu)，使得Fluent進(jìn)行布爾計(jì)算時(shí)多個(gè)界面線性相切，這對(duì)計(jì)算機(jī)性能提出了更高要求，無(wú)法開(kāi)展后續(xù)數(shù)值模擬計(jì)算。因此，需要對(duì)金屬網(wǎng)布三維模型進(jìn)行簡(jiǎn)化并驗(yàn)證其可行性。根據(jù)金屬網(wǎng)布的防砂機(jī)制，保留篩網(wǎng)的多孔流道結(jié)構(gòu)，將其簡(jiǎn)化為圓柱孔布置，并與相同尺寸的原始篩網(wǎng)建模，在附加保護(hù)殼建模的情況下進(jìn)行數(shù)值模擬分析。經(jīng)模擬結(jié)果驗(yàn)證，簡(jiǎn)化后模型與篩網(wǎng)編織模型的沖蝕區(qū)域位置吻合，故簡(jiǎn)化后篩網(wǎng)可以應(yīng)用于后續(xù)數(shù)值模擬。金屬篩網(wǎng)模型簡(jiǎn)化前、后沖蝕云圖分布如圖1所示。

圖1 金屬篩網(wǎng)模型簡(jiǎn)化前、后沖蝕云圖分布

確定金屬網(wǎng)布簡(jiǎn)化模型后，建立復(fù)合篩管的全結(jié)構(gòu)三維模型。其中保護(hù)殼流道采用錯(cuò)位布置，金屬網(wǎng)布網(wǎng)孔直徑為120 μm。金屬網(wǎng)布復(fù)合篩管的三維模型如圖2所示。

圖2 金屬網(wǎng)布復(fù)合篩管三維模型

2.2 網(wǎng)格劃分及邊界條件

2.2.1 網(wǎng)格劃分

為保證計(jì)算結(jié)果精度，對(duì)主體計(jì)算區(qū)域這類(lèi)復(fù)雜的工程結(jié)構(gòu)采取非結(jié)構(gòu)化的四面體網(wǎng)格，并進(jìn)行加密，剩余流體域采用六面體網(wǎng)格劃分方法，并進(jìn)行無(wú)關(guān)性驗(yàn)證，最終確定當(dāng)網(wǎng)格數(shù)為12 877 212個(gè)時(shí)，符合計(jì)算精度要求。最終網(wǎng)格劃分結(jié)果如圖3所示。

圖3 網(wǎng)格劃分結(jié)果

2.2.2 邊界條件

基于CFD模擬金屬網(wǎng)布復(fù)合篩管在固液兩相流作用下的沖蝕過(guò)程，仿真模擬計(jì)算主體區(qū)域是保護(hù)殼和兩層篩網(wǎng)。固液混合物以一定速度沖擊篩管，通過(guò)保護(hù)殼進(jìn)入篩網(wǎng)，由于篩網(wǎng)的擋砂機(jī)制，大量砂粒與篩網(wǎng)的直接碰撞與觸底反彈碰撞使防砂篩網(wǎng)產(chǎn)生塑性變形；部分砂粒進(jìn)入篩網(wǎng)，對(duì)篩網(wǎng)造成切削磨損。固液混合物沖蝕入口為保護(hù)殼外界面，出口為壓力outflow，離散相與沖蝕壁面為彈性反彈。

2.3 壽命預(yù)測(cè)

G.GILLESPIE等[11]通過(guò)對(duì)不同篩管類(lèi)型進(jìn)行完井模擬，提出了基于質(zhì)量損失的篩管防砂失效判斷依據(jù)，認(rèn)為繞絲篩管質(zhì)量損失2%，縫寬變大0.05 mm；金屬絲篩管質(zhì)量損失8%，出砂量會(huì)嚴(yán)重增多，篩管防砂失效。其中8%的質(zhì)量損失極限為保守值，本研究通過(guò)大量模擬與試驗(yàn)對(duì)比，最終取5%為質(zhì)量損失極限。

CFD獲得的沖蝕率定義為單位時(shí)間和單位面積下金屬材料損失的質(zhì)量，故以5%作為質(zhì)量損失極限，以第一層金屬網(wǎng)布的平均沖蝕率為計(jì)算依據(jù)，推算出金屬網(wǎng)布復(fù)合篩管的沖蝕壽命[12]。沖蝕壽命預(yù)測(cè)公式為：

(11)

式中：t為沖蝕破壞時(shí)間，s；Re′為第一層篩網(wǎng)的平均沖蝕率，kg/(m2·s)；ΔGmax為第一層篩網(wǎng)的5%質(zhì)量損失量，kg。

2.4 模擬可行性驗(yàn)證

對(duì)金屬網(wǎng)布篩管試件進(jìn)行不同速度的沖蝕試驗(yàn)，得到不同速度下防砂篩管的沖蝕壽命，進(jìn)而驗(yàn)證全結(jié)構(gòu)金屬網(wǎng)布復(fù)合篩管仿真模擬結(jié)果的可行性。

2.4.1 試驗(yàn)裝置

采用水射流噴射含砂液體沖擊試件表面的方法進(jìn)行沖蝕試驗(yàn)，測(cè)試篩管試件的沖蝕破壞時(shí)間及形貌，分析不同流速、不同篩管結(jié)構(gòu)對(duì)篩管沖蝕壽命的影響規(guī)律。篩管試件沖蝕試驗(yàn)裝置流程如圖4所示。

圖4 沖蝕試驗(yàn)裝置流程圖

2.4.2 試驗(yàn)方法

試驗(yàn)試件采用金屬網(wǎng)布復(fù)合篩管切片，試件尺寸為200 mm×115 mm，擋砂精度為120 μm。砂粒質(zhì)量分?jǐn)?shù)為40%～60%，砂子粒徑為100～120 μm，噴嘴距離為40 mm，沖擊角度為90°，設(shè)定沖擊速度v分別為20、30、40和50 m/s。

試驗(yàn)過(guò)程中，根據(jù)沖擊速度不同，分為多個(gè)輪次，每個(gè)輪次加砂約75 kg，每輪次加砂沖蝕持續(xù)破壞時(shí)長(zhǎng)約35 min。當(dāng)試件過(guò)濾層完全破壞、底部出現(xiàn)明顯漏砂時(shí)，沖蝕破壞試驗(yàn)結(jié)束。累計(jì)篩管的有效沖蝕時(shí)長(zhǎng)，即為該速度下的金屬網(wǎng)布復(fù)合篩管沖蝕壽命。

2.4.3 試驗(yàn)結(jié)果與模擬結(jié)果

沖蝕試驗(yàn)結(jié)束后，拆解并觀察試件的破壞形態(tài)。試件由保護(hù)殼、多層金屬網(wǎng)布和基管3部分組成。通過(guò)觀察發(fā)現(xiàn)，保護(hù)殼出現(xiàn)明顯沖蝕孔，第一層金屬網(wǎng)布出現(xiàn)圓形沖蝕孔，且破壞面積普遍遠(yuǎn)大于外保護(hù)罩的破壞面積。第二層金屬網(wǎng)布破壞形態(tài)與第一層類(lèi)似，但沖蝕的面積減小，內(nèi)部導(dǎo)流層出現(xiàn)一定程度的沖蝕破壞。這是由于試驗(yàn)采取的沖蝕角度為90°，砂礫撞擊試件表面幾乎完全是垂直的應(yīng)力，造成嚴(yán)重的塑形變形[13]。

統(tǒng)計(jì)不同流速下的金屬網(wǎng)布復(fù)合篩管沖蝕壽命，與仿真模擬壽命預(yù)測(cè)結(jié)果相差不大，誤差在10%以內(nèi)，說(shuō)明仿真模擬預(yù)測(cè)金屬網(wǎng)布復(fù)合篩管沖蝕壽命具有可行性。隨著流速增加，沖蝕壽命也在不斷縮短。仿真模擬與物理試驗(yàn)壽命預(yù)測(cè)對(duì)比如圖5所示。

圖5 不同沖擊速度下模擬與試驗(yàn)預(yù)測(cè)壽命對(duì)比

2.5 模擬結(jié)果分析

設(shè)置入口速度為40 m/s、含砂質(zhì)量分?jǐn)?shù)為50%進(jìn)行數(shù)值模擬，得到金屬網(wǎng)布復(fù)合篩管的沖蝕速率分布，如圖6所示。從圖6可以看出，沖蝕速率最大值分布在第一層篩網(wǎng)，且集中在與保護(hù)殼流入通道下的交錯(cuò)區(qū)域。這是由于保護(hù)殼的特殊結(jié)構(gòu)，使得流體流經(jīng)保護(hù)殼進(jìn)入金屬網(wǎng)布的過(guò)程中，流體方向發(fā)生轉(zhuǎn)變，固體顆粒在交錯(cuò)處形成局部熱點(diǎn)，造成該區(qū)域更加容易受到?jīng)_蝕破壞。第一、二層篩網(wǎng)的沖蝕率分布相同，第二層的沖蝕率小于第一層，這說(shuō)明流體從第一層篩網(wǎng)流經(jīng)第二層篩網(wǎng)的過(guò)程中，運(yùn)動(dòng)方向沒(méi)有發(fā)生太多改變，但由于第一層篩網(wǎng)阻隔了大量砂粒，使得第二層篩網(wǎng)的沖蝕程度遠(yuǎn)小于第一層。隨即分別統(tǒng)計(jì)了20～60 m/s沖擊速度下的平均沖蝕率，并分析各個(gè)元件之間的沖蝕規(guī)律，結(jié)果如表1和表2所示。

圖6 沖蝕速率分布云圖

表1 不同沖擊速度下篩管各元件的平均沖蝕速率

表2 不同沖擊速度下篩管各元件之間平均沖蝕速率比值

由表1和表2可以看出：第一、二層篩網(wǎng)與保護(hù)殼的平均沖蝕速率比值分別為2.476和1.053；第一層與第二層篩網(wǎng)的平均沖蝕率比值為2.351，即第一層篩網(wǎng)在沖蝕過(guò)程中最先發(fā)生沖蝕破壞。

3 實(shí)例應(yīng)用

A區(qū)塊油井井深在1 400～2 800 m之間，完井方式為裸眼完井，采用優(yōu)質(zhì)金屬網(wǎng)布防砂篩管，內(nèi)置金屬網(wǎng)布的擋砂精度均為120 μm，井身結(jié)構(gòu)為水平井。根據(jù)各井的產(chǎn)油和產(chǎn)液情況，結(jié)合含水體積分?jǐn)?shù)，對(duì)各井口的產(chǎn)出液密度進(jìn)行計(jì)算，進(jìn)而計(jì)算各井口的質(zhì)量流量、體積流量以及流速。根據(jù)已有的出砂井現(xiàn)場(chǎng)數(shù)據(jù)，使用Rosin-Rammler模型[15]對(duì)現(xiàn)場(chǎng)出砂井砂粒的粒徑分布進(jìn)行分析，并提取出砂井現(xiàn)場(chǎng)油水混合物密度和黏度等信息，綜合分析模擬參數(shù)的設(shè)定方式，得出粒徑中值和最頻粒徑為240 μm左右。含砂體積分?jǐn)?shù)設(shè)定為0.5%。仿真模擬與現(xiàn)場(chǎng)監(jiān)測(cè)壽命對(duì)比見(jiàn)表3。

由表3可以看出，本文建立的沖蝕壽命預(yù)測(cè)模型在實(shí)際應(yīng)用中誤差小于10%，符合工程壽命預(yù)測(cè)要求。

表3 仿真模擬與現(xiàn)場(chǎng)監(jiān)測(cè)壽命對(duì)比

4 結(jié)論及認(rèn)識(shí)

(1)通過(guò)金屬網(wǎng)布復(fù)合篩管仿真模擬獲得篩管沖蝕率，以5%質(zhì)量損失極限作為防砂失效標(biāo)準(zhǔn)，推算出篩管防砂壽命，與試驗(yàn)結(jié)果符合率達(dá)到90%。這為篩管沖蝕壽命預(yù)測(cè)提供了一種新的思路和方法。

(2)不考慮砂粒堵塞，流速與沖蝕壽命呈現(xiàn)負(fù)相關(guān)；第一層篩網(wǎng)為最先被沖蝕破壞的部件。

(3)本次防砂篩管模擬是不考慮砂粒堵塞工況的沖蝕預(yù)測(cè)模型，考慮砂粒堵塞工況的沖蝕模擬仍需進(jìn)一步研究。