甘泉泉 劉少胡 管鋒 李明紅

1.長江大學(xué)機(jī)械工程學(xué)院；2.中國石油西南管道分公司

水平井技術(shù)以其低成本、高產(chǎn)能等特點(diǎn)而廣泛應(yīng)用于油氣藏的開發(fā)，但由于水平井井身結(jié)構(gòu)特殊性，大量滲透進(jìn)井筒的砂粒會(huì)沉積在水平井段，形成大面積的砂床，造成采油困難甚至油井停產(chǎn)等嚴(yán)重井下事故［1-2］。

水平井段沖砂清洗較困難，目前主要采用機(jī)械撈砂和水力沖砂清除砂床，水力沖砂因其效率高、效果好而得到廣泛應(yīng)用，目前水力沖砂洗井技術(shù)主要包括常規(guī)密閉連續(xù)沖砂洗井技術(shù)、連續(xù)管沖砂洗井技術(shù)、真空清砂技術(shù)和泡沫流體沖砂洗井技術(shù)［3-5］。常規(guī)密閉連續(xù)沖砂洗井技術(shù)解決了常規(guī)油管沖砂時(shí)沖砂液回流、卡鉆風(fēng)險(xiǎn)、沖洗液用量大、沖洗時(shí)間過長和無法實(shí)現(xiàn)密閉沖砂作業(yè)等問題；連續(xù)管沖砂洗井工藝無需連接管柱，可以帶壓連續(xù)作業(yè)，對(duì)地層傷害小且洗井效率高；同心連續(xù)管真空清砂技術(shù)對(duì)洗井液的要求較低，可根據(jù)井況優(yōu)選，大幅度降低了作業(yè)成本，避免了因洗井液滲入地層造成的地層傷害；泡沫流體沖砂洗井技術(shù)其攜砂能力遠(yuǎn)高于空氣、清水甚至普通鉆井液，同時(shí)對(duì)附著在油管內(nèi)外壁、套管內(nèi)壁及井下工具上的石蠟及粘結(jié)物有較好的剝離清除作用［6-7］。

上述水力沖砂多采用的是電機(jī)帶動(dòng)沖砂清洗工具旋轉(zhuǎn)實(shí)現(xiàn)脈沖射流，本文提出了一種工具內(nèi)脈沖射流沖砂洗井方法，該方法不需要電機(jī)帶動(dòng)工具旋轉(zhuǎn)產(chǎn)生射流，而是通過工具殼體小孔間歇式開關(guān)實(shí)現(xiàn)脈沖射流，可對(duì)砂粒床進(jìn)行快速高頻破壞清除。采用工具內(nèi)脈沖對(duì)砂床清理的流體動(dòng)力特性發(fā)生了變化，故需論證該脈沖射流方法的可行性，探尋工作參數(shù)對(duì)沖洗砂床的影響規(guī)律。

1 工具內(nèi)脈沖射流沖砂數(shù)值模擬

脈沖射流是非連續(xù)射流，這種間斷脈沖射流在水錘作用下產(chǎn)生較高瞬態(tài)能量，來達(dá)到?jīng)_散堆積砂床的目的［8］。圖1 所示為工具內(nèi)脈沖射流示意圖，當(dāng)液體流入進(jìn)口，工具殼體小孔打開時(shí)分流液體，產(chǎn)生了間歇射流(圖1a)，形成局部高速。當(dāng)工具殼體小孔關(guān)閉時(shí)，工具間歇射流關(guān)閉(圖1b)，此時(shí)射流速度為0，上述兩種工況不間斷轉(zhuǎn)化實(shí)現(xiàn)了工具內(nèi)脈沖射流。

1.1 數(shù)值模型建立及邊界條件

圖1 內(nèi)脈沖沖砂射流原理及工具照片F(xiàn)ig.1 Principle of internal pulsed jet sand washing and the picture of the tool

根據(jù)現(xiàn)場水平井井眼參數(shù)，建立數(shù)值計(jì)算模型。連續(xù)管外徑為50.8 mm，內(nèi)脈沖射流沖砂工具外徑為88.9 mm，井眼環(huán)空內(nèi)徑為127 mm，工具射流小孔尺寸根據(jù)現(xiàn)有井下工具設(shè)計(jì)分布。對(duì)模型進(jìn)行四面體網(wǎng)格劃分，并在工具壁面處和射流孔部位網(wǎng)格插入膨脹層，細(xì)化網(wǎng)格，提高求解精度，具體模型與網(wǎng)格分布如圖2 所示。

圖2 井眼環(huán)空模型網(wǎng)格圖Fig.2 Grid map of borehole annulus model

攜砂液在井筒內(nèi)流動(dòng)可視為不可壓縮兩相流紊態(tài)流動(dòng)［9］，在水平段砂粒初始運(yùn)移速度在一定條件下會(huì)達(dá)到一種動(dòng)平衡狀態(tài)。根據(jù)水力學(xué)中的泥砂瞬時(shí)起動(dòng)流速獲得砂粒沉降速度，計(jì)算公式為［10］

式中，vs為砂粒瞬時(shí)起動(dòng)流速，m/s；d′為油層套管當(dāng)量直徑，m； ρs為 砂粒密度，kg/m3； ρl為 沖洗液密度，kg/m3；g為 重力加速度，m/s2；ds為 砂粒直徑，m；r0為砂床面泥砂干容重；r′為泥砂顆粒的穩(wěn)定干容重；

0hw為薄膜水厚度，mm； δw為 薄膜水厚度參數(shù)，取2.31×10-5cm。

水平段沖砂洗井所需的最小工作排量為［10］

式中，D為井眼環(huán)空外徑，m；d為井眼環(huán)空內(nèi)徑，m。

邊界條件：計(jì)算得到進(jìn)口沖砂液排量Qmin=500 L/min，出口壓力根據(jù)工作井深3 000 m 進(jìn)行預(yù)算，得到其工作壓力為30.2 MPa。井壁和連續(xù)管、沖砂工具壁面假設(shè)為無滑移固定邊界，流相設(shè)置為攜砂液和占比80%的砂床，攜砂液密度為1 200 kg/m3，黏度為26 mPa · s，砂粒初始直徑設(shè)置為0.8 mm，密度為2 500 kg/m3。

1.2 控制方程

采用Realizablek-ε模型對(duì)脈沖射流流場進(jìn)行模擬計(jì)算，質(zhì)量守恒方程為［11］

式中， div(a)= ?ax/?x+?ay/?y+?az/?z； ρ為流體密度，kg/m3；t為時(shí)間，s。

動(dòng)量守恒方程為

式中，u、v、w分別為速度矢量u在x、y、z方 向的分量，m/s；p為流體微元體上的壓力，Pa； τi j為剪切應(yīng)力τ的分量，Pa；Fx、Fy、Fz分 別為微元體在x、y、z方 向上的體力，kg/(m · s)。

1.3 砂床顆粒運(yùn)移受力分析

圖3 為井筒內(nèi)砂粒的受力分析簡圖，砂粒在井筒內(nèi)運(yùn)移受到重力、流體舉升力、拖曳力和射流沖擊力［12］，脈沖射流沖擊力周期性變化。

圖3 水平井筒內(nèi)砂粒受力分析Fig.3 Force analysis on the sand particles in the borehole of horizontal well

砂粒受到的重力

砂粒受到舉升力

砂粒受到拖曳力

射流沖擊力

式中，F(xiàn)g為砂粒重力，N； ρf為沖砂液密度，kg/m3；CL為舉升系數(shù)，無因次；AL為舉升力作用面積，m2；uf為 摩阻流速，m/s；CD為 拖曳系數(shù)，無因次；AD為流體流動(dòng)方向上砂粒投影面積，m2；v為工具射流出口流速，m/s。

2 內(nèi)脈沖射流沖砂與旋轉(zhuǎn)射流沖砂效果對(duì)比分析

根據(jù)內(nèi)脈沖沖砂工具的工作情況，取脈沖幅度比為 42.86% 、占空比比例為1 ∶1［9］，采用水平井段沉積的砂床完全運(yùn)移所用的時(shí)間、使用沖砂工具后沉積在工具下的砂床清潔程度和在沖砂工具尾部平面的砂粒完全流過所需的時(shí)間(即清除砂床所需的時(shí)間)為標(biāo)準(zhǔn)進(jìn)行評(píng)價(jià)。

為對(duì)比內(nèi)脈沖射流與旋轉(zhuǎn)射流沖砂的運(yùn)移規(guī)律，以長度1 000 mm、高度18.5 mm、濃度 80%砂床為研究對(duì)象，采用Eulerian 兩相流模型求解沖砂液-砂床兩相流［13］，以此來研究砂床的運(yùn)移規(guī)律。

在工具尾部建立切片面Plane 1，脈沖前入口處建立了切面Plane 2，兩脈沖入口之間建立了切面Plane 3，脈沖后入口處建立了切面Plane 4，脈沖后入口和射流入口之間建立了切面Plane 5，射流入口處建立切面Plane 6，具體分布如圖4 所示(其中紅色部分表示砂床的分布，藍(lán)色部分為沖砂液的分布)。

圖4 環(huán)空模型切片面分布圖Fig.4 Distribution of slice planes in the annulus model

圖5 所示為射流攜砂和內(nèi)脈沖射流攜砂對(duì)比云圖，在相同工況下，射流2.5 s 后旋轉(zhuǎn)射流沖砂中只有井底靠近噴孔處砂床被運(yùn)移(圖5a 中Plane 6)，其余區(qū)域含砂體積分?jǐn)?shù)約60%；而采用內(nèi)脈沖射流沖砂方法后環(huán)空砂床體積分?jǐn)?shù)全部小于6%，基本被沖洗干凈。由此得出內(nèi)脈沖射流沖砂對(duì)砂床的破壞能力要明顯高于常規(guī)射流沖砂。

圖5 射流沖砂和內(nèi)脈沖沖砂各切片面含砂體積分?jǐn)?shù)Fig.5 Sand volume fraction of each slice plane of jet sand washing and internal pulsed sand washing

由圖6 可知，旋轉(zhuǎn)射流沖砂時(shí)，砂粒被完全運(yùn)移通過Plane 1 所需要的時(shí)間為5.6 s，而內(nèi)脈沖射流沖砂時(shí)，砂粒被完全運(yùn)移通過Plane 1 所需要的時(shí)間為3.5 s，運(yùn)移時(shí)間減少了約37.5%。內(nèi)脈沖清砂的砂粒平均運(yùn)移速度比射流沖砂要快。

圖6 Plane 1 處的含砂體積分?jǐn)?shù)Fig.6 Sand volume fraction at Plane 1

如圖7 所示為內(nèi)脈沖射流時(shí)，水平井眼中砂床運(yùn)移情況。初始時(shí)刻砂床堆積在沖砂工具下方，在t=0.50 s時(shí) 脈沖射流使得砂床在井眼中完全紊亂；t=1.00 s時(shí) 脈沖射流的速度處于波谷，射流速度較小，砂床出現(xiàn)重新堆積；在t=1.50 s時(shí) ，波峰速度來臨，砂床被再次破壞。

圖7 不同時(shí)刻砂床破壞情況Fig.7 Damage situation of sand bed at different moments

3 內(nèi)脈沖射流沖砂參數(shù)的確定

3.1 射流頻率對(duì)沖砂效果的影響

采用控制變量法進(jìn)行對(duì)比分析，由于脈沖頻率過高會(huì)導(dǎo)致小孔射流不具備連續(xù)性，對(duì)砂床的破壞能力不能正常體現(xiàn)，因此設(shè)計(jì)在脈沖射流頻率范圍0.67～2.5 Hz(脈沖沖砂周期為 0.4～1.6 s)內(nèi)優(yōu)選工具脈沖射流頻率。砂粒直徑1 mm、排量500 L/min、砂床高度18.5 mm、砂床長度1 000 mm、砂床體積分?jǐn)?shù)80%，沖砂工具在井眼內(nèi)居中。在0.67～2.5 Hz區(qū)間內(nèi)，由0.4 s 開始，每隔0.2 s 設(shè)置一個(gè)周期點(diǎn)，觀察水平井段沉積在脈沖工具下的砂床完全運(yùn)移離開工具下部所用的時(shí)間。由圖8 可知，在頻率為0.71 Hz 時(shí)，沖砂所用的時(shí)間最少，且砂床的運(yùn)移距離最長，即在該工況參數(shù)下，沖砂效果最好，更有利于砂床的破壞和運(yùn)移。綜上所述，0.71 Hz 為該脈沖射流沖砂的最佳工作頻率。

圖8 不同脈沖頻率時(shí)砂床的運(yùn)移情況Fig.8 Migration situation of sand bed at different pulse frequencies

3.2 沖砂液排量對(duì)沖砂效果的影響

設(shè)定脈沖頻率0.71 Hz、砂粒直徑1 mm、砂床高度18.5 mm。為保證沖砂液在豎直井段將砂粒正常攜帶出井，排量需不小于500 L/min［8］。實(shí)驗(yàn)研究沖砂液排量分別為550、600、650、700 和750 L/min時(shí)射流沖砂效果和砂床的運(yùn)移情況。

從圖9 可以看出，井筒內(nèi)的含砂體積分?jǐn)?shù)隨著沖砂液排量V的增大而減小，當(dāng)V=750 L/min 時(shí)，截面內(nèi)固定砂床體積減少到0，總的含砂體積分?jǐn)?shù)不大于15%，同時(shí)固定砂床的體積明顯減少。原因在于當(dāng)沖砂液排量增大，而沖砂工具所處的環(huán)空面積沒有變化，使得環(huán)空內(nèi)和射流孔處的沖砂液流速增大，同時(shí)射流沖擊力正比于射流速度的平方，最終增大套管內(nèi)砂床的運(yùn)移距離，提高沖砂效率。

圖9 不同排量下截面含砂體積分?jǐn)?shù)Fig.9 Sand volume fraction of the section plane at different displacements

3.3 砂床高度對(duì)沖砂效果的影響

脈沖射流參數(shù)為：頻率0.71 Hz、排量500 L/min、砂粒直徑1 mm，初始時(shí)刻砂床高度H分別為6.5、9.5、12.5、15.5、18.5 mm 時(shí)套管中砂床及工具的截面分布如圖10 所示。

圖10 初始時(shí)刻套管中砂床及工具的截面分布Fig.10 Distribution of section plane of tool and sand bed in the casing at the initial moment

從圖11 可以看出，砂床高度為6.5 mm 和18.5 mm時(shí)，砂床的運(yùn)移距離較大，而其他砂床高度時(shí)運(yùn)移距離較短。這是由于砂床高度小于9.5 mm 時(shí)，井段的含砂量較低，尚未形成固定砂床，在脈沖射流作用下運(yùn)移較快。砂床高度大于9.5 mm 時(shí)，隨著砂床高度的增加，井筒環(huán)空過流面積減少，沖砂液流速增大，作用在砂床上的沖擊力也隨之增加，進(jìn)而增大了砂床的運(yùn)移距離，提高了沖砂效率。因此，在實(shí)際操作過程中為保證沖砂效率，建議在砂床高度小于9.5 mm時(shí)沖砂作業(yè)。

圖11 不同砂床高度運(yùn)移情況Fig.11 Migration situation of different sand bed heights

3.4 砂粒直徑對(duì)沖砂效果的影響

設(shè)沖砂液排量為500、550、600 L/min，砂床高度H為18.5 mm，脈沖頻率為0.71 Hz，研究砂粒直徑對(duì)砂床運(yùn)移影響規(guī)律［14］。從圖12 可以看出，在3 種排量下，在砂粒直徑小于0.8 mm 時(shí)砂床運(yùn)移距離較長，直徑大于0.8 mm 時(shí)，砂床運(yùn)移距離變化較短，其中在排量為500 L/min 和550 L/min 時(shí)砂床的運(yùn)移距離基本沒變化。粒徑小于0.8 mm 時(shí)，砂粒的舉升力足以將其“舉起”并沿軸向運(yùn)移。因此得出砂粒直徑小于0.8 mm 時(shí)沖砂洗井效率較優(yōu)。

圖12 不同砂粒直徑運(yùn)移情況Fig.12 Migration situation of different sand particle diameters

4 結(jié)論

(1)研制出了內(nèi)脈沖射流工具，對(duì)比分析了內(nèi)脈沖射流沖砂和旋轉(zhuǎn)射流沖砂的砂床運(yùn)移規(guī)律，研究認(rèn)為采用內(nèi)脈沖射流沖砂比旋轉(zhuǎn)射流沖砂效率要高，同樣工況下，內(nèi)脈沖射流沖砂對(duì)砂床的破壞力更大，所需時(shí)間可縮短37.5%。

(2)通過內(nèi)脈沖射流沖砂對(duì)沖砂效果的影響分析，得到以下認(rèn)識(shí)：優(yōu)選脈沖頻率為0.71 Hz；選用較大沖砂液排量，有利于高效沖砂；砂床高度小于9.5 mm和砂粒直徑小于0.8 mm 時(shí)沖砂洗井效率較優(yōu)。