封　猛，丁心魯，王金禮，徐偉紅，張曉文，盧德唐

(1.中國石油西部鉆探工程有限公司試油公司，新疆 克拉瑪依 834000；2.中國科學技術(shù)大學近代力學系，安徽 合肥 230026)

在油氣田勘探與開發(fā)中，針對壓力系數(shù)低的非自噴稠油井、含硫化氫井、氣油比大的井，多采用螺桿泵求產(chǎn)工藝[1-2]。為了獲取螺桿泵求產(chǎn)井的地層參數(shù)，現(xiàn)場形成了螺桿泵+偏心井口環(huán)空測試、封隔器+MFE(APR)多流測試器+螺桿泵多聯(lián)作等技術(shù)。但由于受到井下、井口生產(chǎn)工藝及其他多種因素制約[3]，難以推廣應(yīng)用。采用螺桿泵求產(chǎn)，獲取地層參數(shù)是現(xiàn)場生產(chǎn)急需解決的問題。

螺桿泵試產(chǎn)分為開泵求產(chǎn)階段和停泵液面恢復(fù)階段，井底流動特征符合間歇試產(chǎn)分析條件[4-8]，采用數(shù)值方法，求解地層壓力分布及井底壓力，可實現(xiàn)地層參數(shù)解釋。

1　間歇試產(chǎn)試井分析原理

間歇方式采油井有以下幾方面特征：①無論是抽油還是停抽階段，地層都不斷地有液體產(chǎn)出，地層的壓力分布都不斷地變化，地層壓降漏斗是由多次抽油形成壓力波的迭加；②和段塞流動期間相同，這種生產(chǎn)方式無法確定井底產(chǎn)量，即產(chǎn)量和井底壓力都是未知的。井筒中的流動和地層滲流是偶合的，求出各時刻的壓力分布、井底壓力，就可以實現(xiàn)對地層進行解釋求參。

1.1　流體流量的計算

螺桿泵求產(chǎn)的回收量(V1)由兩部分組成，一部分為恢復(fù)期間已存在的液體，另一部分是地層流入井筒中的流體，設(shè)生產(chǎn)Δt時間后的井底壓力為p1，生產(chǎn)Δt時間后井底壓力為p2，則回收量為V2:

(1)

式中：A為井筒截面積，m2；ρ為液體密度，g/cm3；其中地層流入井筒的部分為：

(2)

利用達西公式，可以得到地層流入井筒里總量：

(3)

1.2　井底壓力的計算

以均質(zhì)無限大地層為例，停泵期間的方程及其定解條件可以表示成：

(4)

p(Δt=0,r)=f(r)

(5)

p(r→∞,t)=pi

(6)

(7)

則井底壓力為：

(8)

式中f(r)為地層初始時刻的地層壓力分布(即壓降漏斗)，可由(3)得到。

1.3　壓力歷史擬合

由于存在多個求產(chǎn)及停泵過程，全壓力計算較復(fù)雜，采用數(shù)值方法進行計算。

考慮到油井為圓形，包含井的網(wǎng)格，常稱之為井網(wǎng)格。當網(wǎng)格是井網(wǎng)格，其它網(wǎng)格與井網(wǎng)格的傳導系數(shù)要特殊處理，如圖1、圖2所示。

圖1井網(wǎng)格與相鄰網(wǎng)格

圖2等效外半徑，rw為井半徑

井網(wǎng)格與其相鄰網(wǎng)格的流動仍是徑向流，相鄰網(wǎng)格j從rj處按徑向流方式流向井網(wǎng)格i。顯然要計算出井網(wǎng)格的平均壓力所對應(yīng)的半徑rc，即等效半徑rc。對井網(wǎng)格，傳導系數(shù)應(yīng)是從rj處流到rc處。下面給出rc的計算方法。

整個井網(wǎng)格i的面積A可用等效圓表示，為區(qū)別等效半徑，稱等效圓的半徑為等效外接圓半徑，如圖2所示，即ro為：

(9)

其中A為井網(wǎng)格i的面積，β為射開的角度，當β=2π為全射開，則井網(wǎng)格的等效半徑為

(10)

(11)

在網(wǎng)格i中，對方程(11)作體積分：

(12)

應(yīng)用高斯定理，則方程(12)左邊的體積分可以變?yōu)槊娣e分，

(13)

相鄰網(wǎng)格中心點的連線垂直兩個網(wǎng)格的相鄰面，因而每個面的法線與流動方向平行，這使面積分大大簡化。由于認為壓力、粘度等參數(shù)在相鄰面上是常數(shù)，因而有

(14)

(15)

其中dij為網(wǎng)格中心距離，如圖3所示。

圖3網(wǎng)格參數(shù)示意

對累積項的積分，有

(16)

聯(lián)立方程(15)與方程(16)，單相流體的滲流方程(11)變?yōu)椋?/p>

(17)

其中pj，pi為網(wǎng)格i與網(wǎng)格j的壓力,Vi為i網(wǎng)格的體積，Dj=γZj為地下原油重度與從某一基準面算起的高度Zj的乘積。Tij為傳導系數(shù)，是網(wǎng)格流動系數(shù)γij與其幾何因子Gij的乘積。

建立網(wǎng)格后，對數(shù)值離散后形成的線性方程組進行求解，可以得到地層壓力分布及井底壓力，由井底壓力即可以實現(xiàn)壓力歷史擬合。

2　應(yīng)用實例

C21井區(qū)位車排子油田北段，位于克拉瑪依市以南約40 km處，該區(qū)塊石炭系地面原油密度為0.902～0.937 g/cm3；50℃時粘度平均為363.14 mPa·s。開發(fā)區(qū)塊采用螺桿泵機抽方式求產(chǎn)，經(jīng)三個月的試采，平均日產(chǎn)油7.5 t，取得較好效果。

C222井位于C21井區(qū)內(nèi)，目的層位石炭系，巖性為綠灰色熒光凝灰質(zhì)細砂巖、深灰色熒光凝灰?guī)r。射孔井段1 311.00～1 341.00 m，射后采用胍膠液458.0 m3，石英砂60.0 m3進行措施改造，自噴退液時出口發(fā)現(xiàn)硫化氫，采取除硫措施后下螺桿泵(底帶電子壓力計)機抽求產(chǎn)，平均日抽出油2.12 m3，日抽出壓裂液2.81 m3，(螺桿泵轉(zhuǎn)速30 r/min)，壓力計監(jiān)測全過程井下壓力、溫度(如圖4)。

圖4實測壓力、溫度曲線

螺桿泵求產(chǎn)期間，具有多個間歇停泵階段，560小時以后，間歇求產(chǎn)井底壓力曲線相對平滑，選擇該階段進行地層參數(shù)解釋。利用全壓力歷史擬合軟件，在壓力的雙對數(shù)導數(shù)圖可以看出，1/2斜率線特征明顯，選擇“無限導流垂直裂縫+均質(zhì)油藏+無限大邊界”模型進行擬合，見圖5、圖6。

圖5雙對數(shù)導數(shù)曲線

圖6壓力歷史擬合

求得地層參數(shù)如下：地層壓力為14.612 3 MPa，地層滲透率為0.131 1×10-3μm2，表皮系數(shù)為-6.57，裂縫半長為44.52 m。。分析結(jié)果反映儲層物性較差，措施改造增大了泄油面積，供液能力有所改善，表皮系數(shù)反映井周不存在堵塞。

3　結(jié)論

(1)采用全壓力歷史擬合方法，可以對螺桿泵間歇求產(chǎn)井的壓力資料進行擬合求參，解決了氣油比大的井、非自噴稠油井的地層參數(shù)解釋，為下步開發(fā)方案提供依據(jù)。

(2)全壓力歷史擬合解釋方法，可以不需要關(guān)井壓力數(shù)據(jù)進行地層解釋，采用螺桿泵底帶電子壓力計即可實現(xiàn)數(shù)據(jù)采集，比多流測試器+螺桿泵多聯(lián)作螺桿泵測試及螺桿泵+偏井口環(huán)空測試技術(shù)優(yōu)勢明顯，優(yōu)化了螺桿泵測試工藝。