王新濤， 盛磊祥， 鞠少棟， 王世強(qiáng)， 邢洪憲， 劉傳剛

（1.中海油能源發(fā)展股份有限公司 工程技術(shù)分公司，天津 300452；2.中海石油（中國(guó)）有限公司 北京研究中心，北京 100020）

0 引言

海上油田防砂井大部分采用了篩管簡(jiǎn)易防砂和篩管礫石充填的防砂方式，油氣井開(kāi)發(fā)至后期，儲(chǔ)層近井地帶地層顆粒、膠質(zhì)等向井筒周?chē)\(yùn)移聚集，堵塞篩管，造成其滲透率降低，最終導(dǎo)致油井產(chǎn)量下降，甚至停產(chǎn)[1]。目前解除地層堵塞主要采用的是物理解堵法和化學(xué)解堵法，物理解堵法污染較小，解堵效果較差；化學(xué)解堵方法投資較大、工藝復(fù)雜，易對(duì)地層造成二次污染。

設(shè)計(jì)井下解堵輔助配套工具，解決物理解堵方法中解堵效果較差的問(wèn)題[2]，主要是通過(guò)提供周向高頻振動(dòng)，清理篩管外表面沉積砂及篩管網(wǎng)布堵塞的細(xì)粉砂，從而提高油氣井解堵效率。通過(guò)對(duì)井下解堵輔助配套工具進(jìn)行理論計(jì)算、流場(chǎng)模擬分析及強(qiáng)度校核，驗(yàn)證其提供周向高頻振動(dòng)、輔助解堵的可行性。

1 井下解堵輔助配套工具流場(chǎng)分析

1.1 工具流體計(jì)算理論

井下解堵輔助配套工具采用流體作為驅(qū)動(dòng)介質(zhì)。流體進(jìn)入工具內(nèi)部，分成3部分：經(jīng)過(guò)工具中部流道下端的噴嘴繼續(xù)向下流動(dòng)；從工具中部流道分流一部分，進(jìn)入沖擊腔，推動(dòng)沖擊錘周向轉(zhuǎn)動(dòng)；進(jìn)入換向通道。與其他兩個(gè)通道相比，換向通道體積很小，計(jì)算工具的水力學(xué)參數(shù)時(shí)可以忽略。

工具內(nèi)部流體的分流情況如圖1所示，中部流道的流體和沖擊腔內(nèi)的流體為并聯(lián)的關(guān)系。根據(jù)流體并聯(lián)關(guān)系的特點(diǎn)，沖擊腔內(nèi)的流體和中部流道所造成的壓降相等。

由于整個(gè)工具的長(zhǎng)度較短，同時(shí)工具表面較為光滑，因此工具的摩擦阻力可以忽略，整個(gè)輔助配套工具壓降的主要來(lái)源是工具內(nèi)部流體流向改變和截面積突變的局部壓降[3]。

1）中部流道。進(jìn)入中部流道的流體通過(guò)中部噴嘴處，截面積突減，造成的壓降為：

式中：Δp11為噴嘴處所造成的壓降；ξ11為工具局部壓降系數(shù)；ρ為流道中流體密度；v11為中部流道噴嘴出口處液流流速；Q1為從中部流道經(jīng)過(guò)的流體的流量；A11為噴嘴的過(guò)流面積。

2）沖擊腔。工具內(nèi)部有2個(gè)完全相同且對(duì)稱(chēng)的沖擊腔，這2個(gè)沖擊腔是并聯(lián)關(guān)系，分析時(shí)只需要計(jì)算1個(gè)沖擊腔。流體從軸向運(yùn)動(dòng)的方向進(jìn)入沖擊腔內(nèi)，流體的運(yùn)動(dòng)方向改變90°，同時(shí)流道面積也有突變，產(chǎn)生的壓降為：

式中：Δp31為該處所造成的壓降；ξ31為局部壓降系數(shù)；ξ32為截面積突變時(shí)的局部壓降系數(shù)；v31為從沖擊腔通道內(nèi)排出流體的流速；Q2為進(jìn)入單個(gè)沖擊腔內(nèi)流體的流量；A31為沖擊腔排出通道的過(guò)流面積。

排出后的流體從下端的通道流出并和中部通道流出的流體匯合，這部分的流體截面積發(fā)生突變，產(chǎn)生的壓降為：2

中部噴嘴的過(guò)流面積越小，工具的壓降越大，從噴嘴處流到下端的流體的流量越小，進(jìn)入沖擊腔內(nèi)分流的液流排量越大，沖擊錘振動(dòng)的頻率越高，所做的沖擊功就越大。

沖擊腔入口流道尺寸越小，進(jìn)入沖擊腔和從沖擊腔排出的阻力越大，損耗的能量越大。流道面積越大越好，但是過(guò)大會(huì)降低沖擊錘強(qiáng)度，因此在設(shè)計(jì)時(shí)需要綜合考慮壓降和工具強(qiáng)度。

1.2 工具內(nèi)部流場(chǎng)模擬

為分析工具內(nèi)部流體在流道中運(yùn)動(dòng)及壓力狀態(tài)，使用流體力學(xué)軟件FLUENT進(jìn)行數(shù)值模擬計(jì)算。建立工具三維的內(nèi)部流道結(jié)構(gòu)，如圖2和圖3所示。

圖2 井下解堵輔助配套工具內(nèi)部流道三維結(jié)構(gòu)

圖3 井下解堵輔助配套工具不同位置處橫截剖面

建立CFD流體動(dòng)力學(xué)模型，通過(guò)Gambit軟件對(duì)輔助工具流體域建立網(wǎng)格模型[4]。

計(jì)算時(shí)工具內(nèi)部通道液流取水的密度值，液流的動(dòng)力黏度取0.001 kg/（m·s）。示例中，井下解堵輔助配套工具取中心孔的直徑為30 mm，從中心孔進(jìn)入的流速取30 L/s，可以看到此時(shí)為強(qiáng)湍流的流動(dòng)狀態(tài)，因?yàn)檩o助配套工具內(nèi)部的雷諾數(shù)高達(dá)2.48×105。進(jìn)口條件：速度進(jìn)口取0.6114 m/s，湍流強(qiáng)度取3.39%，水力直徑取0.154 m。出口條件：壓力出口。壁面條件：取無(wú)滑移壁面的邊界條件。

沖擊錘在流體壓力的推動(dòng)下運(yùn)動(dòng)，從開(kāi)始運(yùn)動(dòng)到滿(mǎn)程時(shí)總共所需時(shí)間為0.011 72 s，工具的沖擊頻率為42.66 Hz，模擬計(jì)算整套井下解堵輔助配套工具的進(jìn)出口壓力差為1.2～2.7 MPa。

井下解堵輔助配套工具的沖擊錘在運(yùn)動(dòng)過(guò)程中的壓差衰減曲線如圖4所示，井下解堵輔助配套工具的沖擊錘在流體壓力的推動(dòng)下運(yùn)動(dòng)時(shí)，剛開(kāi)始的流阻是最大的，沖擊錘啟動(dòng)之后，內(nèi)部流體進(jìn)行分流，工具的流阻逐漸減小；在結(jié)構(gòu)一定的情況下，輔助配套工具總體的壓差損失隨排量的減小而降低。

圖4 井下解堵輔助配套工具壓差隨時(shí)間變化曲線

不同流量下井下解堵輔助配套工具在中心孔處和分流孔處的液流的速度-時(shí)間變化規(guī)律曲線如圖5所示，對(duì)比不同的排量，在一定的時(shí)間范圍內(nèi)的同一時(shí)間點(diǎn)時(shí)，排量增加，配套工具在中心孔處和分流孔處的液流的速度是變大的。

圖5 井下解堵輔助配套工具不同排量下不同位置速度-時(shí)間變化曲線

分析不同排量下井下解堵輔助配套工具的運(yùn)動(dòng)特性的變化規(guī)律，需要重點(diǎn)研究沖擊錘轉(zhuǎn)矩和沖擊錘加速度的特性變化規(guī)律。如圖6所示，在剛開(kāi)始啟動(dòng)的時(shí)刻，輔助配套工具的沖擊錘所受轉(zhuǎn)矩值及加速度值均為最大，沖擊錘啟動(dòng)之后，內(nèi)部流體進(jìn)行分流，沖擊錘所受的轉(zhuǎn)矩值及加速度值均逐漸減小。對(duì)比不同排量，在一定的時(shí)間范圍內(nèi)的同一時(shí)間點(diǎn)時(shí)，排量增加，配套工具的沖擊錘轉(zhuǎn)矩和加速度均增加。

圖6 井下解堵輔助配套工具不同流量下沖擊錘轉(zhuǎn)矩及沖擊錘加速度變化曲線

2 井下解堵輔助配套工具強(qiáng)度分析

2.1 沖擊錘結(jié)構(gòu)碰撞力分析

當(dāng)液流通過(guò)中部通道分流進(jìn)入沖擊腔內(nèi)，液流推動(dòng)沖擊錘開(kāi)始轉(zhuǎn)動(dòng)。沖擊錘的運(yùn)動(dòng)關(guān)系為：1）初始狀態(tài)時(shí)，沖擊錘為靜止?fàn)顟B(tài)。當(dāng)入口通道打開(kāi)，液流進(jìn)入工具的沖擊腔，沖擊錘開(kāi)始轉(zhuǎn)動(dòng)。這時(shí)沖擊錘的轉(zhuǎn)速較低，沖擊錘低壓腔內(nèi)的流體排出的速度較低，沖擊錘的阻尼力較小，沖擊錘的加速度很大。2）沖擊錘的轉(zhuǎn)速越高，沖擊錘低壓腔內(nèi)的流體排出的速度越高，沖擊錘的阻尼力越大。當(dāng)阻尼力與驅(qū)動(dòng)力大小相等時(shí)，此時(shí)沖擊錘受力為零，沖擊錘做勻速轉(zhuǎn)動(dòng)，直至撞擊沖擊腔壁面。3）以恒定速度撞擊到錘外殼上時(shí)，該沖擊功會(huì)通過(guò)錘外殼下端的方形塊傳遞到下部。同時(shí)沖擊錘的速度變?yōu)榱?，這時(shí)啟動(dòng)器繼續(xù)轉(zhuǎn)動(dòng)，改變高低壓腔通道方向，驅(qū)動(dòng)沖擊錘向著相反的方向轉(zhuǎn)動(dòng)。

沖擊錘的臨界轉(zhuǎn)速為：

式中：f為沖擊頻率；A24為擺錘葉片的側(cè)面積；m為擺錘質(zhì)量。

仿真計(jì)算的邊界條件設(shè)定為通過(guò)Fluent動(dòng)網(wǎng)格模擬計(jì)算出的液流對(duì)沖擊錘的沖擊速度，以此計(jì)算沖擊錘的撞擊力。在極限狀態(tài)下，最大碰撞力發(fā)生在沖擊錘與壁面撞擊后反彈的時(shí)候。沖擊錘在液流作用下與沖擊腔的碰撞分為彈性碰撞和阻尼碰撞，即撞擊后沖擊錘反彈和速度為零的2種情況[5]。

經(jīng)過(guò)模擬計(jì)算，輔助配套工具的沖擊錘彈性碰撞和阻尼碰撞轉(zhuǎn)矩曲線如圖7所示，工具沖擊錘彈性、阻尼碰撞的轉(zhuǎn)矩及應(yīng)力數(shù)據(jù)如表1所示。

圖7 輔助配套工具沖擊錘彈性碰撞和阻尼碰撞轉(zhuǎn)矩曲線

表1 輔助配套工具沖擊錘彈性碰撞和阻尼碰撞數(shù)據(jù)

2.2 工具結(jié)構(gòu)強(qiáng)度校核

通過(guò)模擬計(jì)算來(lái)分析輔助配套工具本體的疲勞強(qiáng)度。輔助配套工具本體材質(zhì)為4145H，彈性模量為210 GPa，泊松比為0.3，密度為7850 kg/m3，通過(guò)計(jì)算本體最大等效應(yīng)力值為149 MPa，最大靜變形為0.201 mm，如圖8所示。工具整體滿(mǎn)足強(qiáng)度要求。

圖8 輔助配套工具本體等效應(yīng)力和整體變形云圖

沖擊錘等效應(yīng)力分布如圖9所示，在沖擊錘的邊界處及換向槽的棱角處出現(xiàn)了一定的應(yīng)力集中，最大應(yīng)力值為128 MPa，而沖擊錘所采用的材料為20CrMoTi，最大應(yīng)力遠(yuǎn)小于屈服強(qiáng)度，因此沖擊錘的強(qiáng)度滿(mǎn)足需求。

圖9 沖擊錘等效應(yīng)力分布

3 結(jié)語(yǔ)

1）在方案設(shè)計(jì)的基礎(chǔ)上進(jìn)行工具內(nèi)部流場(chǎng)模擬，工具沖擊錘在流體壓力的推動(dòng)下運(yùn)動(dòng)時(shí)，剛開(kāi)始的流阻最大，沖擊錘啟動(dòng)之后，內(nèi)部流體進(jìn)行分流，工具的流阻逐漸減小；結(jié)構(gòu)一定的情況下工具總體的壓差損失隨排量的減小而降低。2）研究流量對(duì)工具運(yùn)動(dòng)特性影響，在剛開(kāi)始啟動(dòng)的時(shí)刻，輔助配套工具的沖擊錘所受轉(zhuǎn)矩及加速度均為最大，沖擊錘啟動(dòng)之后，內(nèi)部流體分流，轉(zhuǎn)矩及加速度均逐漸減小；在一定的時(shí)間范圍內(nèi)，排量增加，配套工具的沖擊錘轉(zhuǎn)矩和加速度均增加。3）對(duì)工具進(jìn)行碰撞力和強(qiáng)度分析，工具整體滿(mǎn)足強(qiáng)度要求，沖擊錘強(qiáng)度亦滿(mǎn)足需求。