邢翠巧 尹洪軍 李興科 劉岢鑫 付京

摘要：中國已開發(fā)的大多數(shù)油田，經(jīng)過多年的加密調(diào)整和注水開發(fā)，目前已處于高含水或高采出程度階段，開采和穩(wěn)產(chǎn)難度越來越大。因此，對注水開發(fā)油田進(jìn)行水驅(qū)面積波及系數(shù)的研究意義重大。根據(jù)注水油田的基本特點(diǎn)，基于流管法建模基本原理，建立了五點(diǎn)法流管模型，利用Beckley-Leverett方程，結(jié)合單元分析法和流線積分法，對注水油田五點(diǎn)法井網(wǎng)面積波及效率的計(jì)算公式進(jìn)行了詳細(xì)推導(dǎo)，繪制并分析了面積波及系數(shù)和采出程度隨時(shí)間的變化曲線。結(jié)果表明：當(dāng)無流管突破時(shí)，注水油田的面積波及系數(shù)隨時(shí)間的增加而線性增大；當(dāng)主流管突破見水后，波及系數(shù)增大程度越來越慢，最后趨于平緩。當(dāng)油田被注入水完全波及后，繼續(xù)利用注水提高采收效果的作用很小，水驅(qū)效率明顯降低。研究結(jié)果對于注水油田面積波及系數(shù)的求解及水驅(qū)動(dòng)態(tài)開發(fā)規(guī)律的研究具有借鑒意義。

關(guān)鍵詞：滲流力學(xué)；流管法；面積波及效率；水驅(qū)動(dòng)態(tài)開發(fā)；流線積分法

中圖分類號：TE344文獻(xiàn)標(biāo)志碼：Adoi： 10.7535/hbgykj.2018yx04001

水驅(qū)面積波及系數(shù)是描述油田開發(fā)效果的重要參數(shù)之一，提高水驅(qū)波及面積是提高原油采收率的重要手段，代表了水驅(qū)過程中注水驅(qū)替的好壞程度，對于研究、計(jì)算不同注采井網(wǎng)類型的水驅(qū)驅(qū)油效率和最終采收率等參數(shù)具有重要意義。所以，雖然面積波及系數(shù)的計(jì)算難度大，但是國內(nèi)外已有很多學(xué)者對其進(jìn)行了研究[1-3]。

用來研究水驅(qū)規(guī)律的最早的實(shí)驗(yàn)技術(shù)方法是電解模型，由MUSKAT等[4]于1934年最先提出，用電流來演示穩(wěn)定滲流。隨著對水驅(qū)滲流規(guī)律的不斷研究，吸墨紙式模型和明膠模型相繼出現(xiàn)，但是這些模型仍是電解模型的延續(xù)。DOYLE等[5]利用流線位置相似的特性，結(jié)合電位模型研究獲得的技術(shù)，生成了一個(gè)簡化網(wǎng)絡(luò)，從流線定義出發(fā)，擴(kuò)展了流動(dòng)通道的形狀。他們利用這種方法進(jìn)行了經(jīng)濟(jì)可行性研究，其中對于儲層的最佳模式和密度可以有多個(gè)選擇，并通過與實(shí)驗(yàn)室模型相比較，建立了流度比的計(jì)算程序。尹洪軍等[6]對基于流管法的注水油田開發(fā)進(jìn)行了數(shù)值模擬研究，通過建立“一注一采”的概念模型，求解并繪制該概念模型的流管法滲流圖版，給出了水驅(qū)油田動(dòng)態(tài)開發(fā)分析方法，對于準(zhǔn)確描述和評價(jià)小層的動(dòng)用程度和水驅(qū)波及效率具有重要的應(yīng)用意義。目前，對于注水油藏水驅(qū)動(dòng)態(tài)及流體滲流規(guī)律的研究不斷深入，研究方法主要包括油藏工程方法、實(shí)驗(yàn)方法和數(shù)值模擬方法[7-8]。牛彥良等[9]通過對井網(wǎng)系數(shù)、驅(qū)油效率和波及系數(shù)等因素的分析，結(jié)合非達(dá)西滲流理論，研究了注水油藏水驅(qū)采收率計(jì)算方法。范江等[10]運(yùn)用概率論和量綱分析方法建立了波及系數(shù)計(jì)算模型。該方法中有大量的統(tǒng)計(jì)學(xué)公式，難以應(yīng)用到礦場分析中，不利于推廣。近年來國內(nèi)外學(xué)者針對注水油藏提出了一些計(jì)算面積波及系數(shù)的公式，給出了面積波及系數(shù)的計(jì)算方法[11-15]，但計(jì)算公式應(yīng)用的大多是標(biāo)定值，是通過與地質(zhì)或井網(wǎng)參數(shù)回歸得到的，并沒有考慮到其他參數(shù)對面積波及系數(shù)的影響，因此，在一定程度上存在不合理性。

流線模擬最大的優(yōu)勢就是計(jì)算速度快，且精度高，具有良好的可視性。流線模擬的計(jì)算過程都是沿流線進(jìn)行的，在每一個(gè)網(wǎng)格尺寸內(nèi)可以直觀地描述流體從何處來，將到何處去，這就相當(dāng)于加密網(wǎng)格的情形，但是并沒有增加大量的計(jì)算量。因此，本文考慮將流線方法結(jié)合到五點(diǎn)法井網(wǎng)的注水驅(qū)替問題中，給出水驅(qū)過程中面積波及效率等動(dòng)態(tài)參數(shù)的變化特征。

1流管法

流管法就是假設(shè)在注水驅(qū)替過程中，油藏的物理驅(qū)替機(jī)理是可以用注采井之間的一根根不相交的流管來表示。流線表示的就是流體過程中的流動(dòng)軌跡，兩根流線組成1根流管，在研究區(qū)塊內(nèi)的注采井之間有規(guī)律的分布著多根流管。本文用形狀因子來考慮每一根流管的形狀，其大小與地層性質(zhì)等有關(guān)。

面積注水是指將注水井和采油井按一定的幾何形狀和密度均勻地布置在整個(gè)開發(fā)區(qū)上進(jìn)行注水和采油。五點(diǎn)井網(wǎng)是當(dāng)注采比為1時(shí)的井位布置形式，如圖1所示。

模型假設(shè)條件：

1）油層是均質(zhì)、等厚的單一油層；

2）油井和水井之間由流管組成，流管中的流體為水驅(qū)油過程；

3）水驅(qū)油過程仍然遵循已經(jīng)確定的多孔介質(zhì)中均質(zhì)流體穩(wěn)定滲流時(shí)的流線；

4）不考慮地層巖石及流體的壓縮性；

5）不考慮毛管力和重力的作用；

6）油井以定產(chǎn)量生產(chǎn)。

根據(jù)以上假設(shè)條件，在每一口油井與水井之間建立多根流管，通過圖1及五點(diǎn)法井網(wǎng)流體滲流特征及對稱性的分析，取五點(diǎn)法井網(wǎng)中的1/8單元作為計(jì)算單元，見圖2。

在圖2中，油、水井半徑為rw，油、水井距為l，油藏厚度為h。任意一根三角流管的夾角為Δα，流管角平分線由對稱的L1和L2組成，它與主流線的夾角為α。沿三角流管中線上的位置表示的是注水驅(qū)替過程中流管內(nèi)前緣移動(dòng)的位置，用坐標(biāo)ξ表示。

2數(shù)學(xué)模型

根據(jù)上述所建立的流管模型，研究對應(yīng)的數(shù)學(xué)模型，可以將水驅(qū)油藏五點(diǎn)法井網(wǎng)的二維兩相驅(qū)替問題簡化為若干個(gè)一維兩相滲流問題，以降低計(jì)算求解的難度，且有利于實(shí)現(xiàn)后期的編程計(jì)算。

第i根流管的角平分線與主流線的夾角為α（i），則

α（i）=Δα·（i-12）=π4N（i-12）。 （1）

流管角平分線的長度為

L（i）=lcos α（i） 。 （2）

以下將L（i）簡記為L。

角平分線上某一位置ξ處的橫截面積為

A（ξ）=2tanΔα2h·ξ，rw≤ξ≤L2；

2tanΔα2h·（L-ξ）， L2<ξ≤L-rw。（3）

在確定了流管的基本參數(shù)之后，給出了水驅(qū)油藏五點(diǎn)法井網(wǎng)面積波及系數(shù)的求解方法：模型中有N根流管，并且每根流管都是由2個(gè)對稱的等腰三角形拼接而成，流管流動(dòng)區(qū)域所占的面積表示的是去除注水井和采油井所占的面積后其余的面積。每一根流管之間均互不交叉，第i根流管可流動(dòng)區(qū)域的面積記為Stube（i），該流管中劃分成的每個(gè)小塊的面積記為Scell（i），所有流管可流動(dòng)區(qū)的總面積記為Stotal。將每根流管等體積的分為2M份，根據(jù)幾何關(guān)系，可得：

Stube（i）=2tanΔα2（14L2-r2w），（4）

Scell（i）=Stube（i）2M， （5）

Stotal=∑Ni=1Stube（i）。 （6）

記第i根流管內(nèi)已經(jīng)波及到的區(qū)域的面積為SAC（i），前緣Lf（i）處于第k（i）個(gè)小塊內(nèi)（以下將Lf（i）簡記為Lf），計(jì)算SAC（i）需要分不同的情況進(jìn)行討論，具體如下。

1）當(dāng)前緣未過流管中線半長時(shí)

SAC（i）=Scell（i）（k-1）+Strapezoid， （7）

其中

Strapezoid=2tan（Δα2）（yk-1+Lf）·Lf-yk-12。

2）當(dāng)前緣過流管中線半長但未突破時(shí)

SAC（i）=Scell（i）（k-1）+Strapezoid， （8）

其中

Strapezoid=2tan（Δα2）[（L-yk-1）+

（L-Lf）]·Lf-yk-12。

3）當(dāng)流管突破時(shí)

SAC（i）=Stube（i）。 （9）

綜上所述，整理得到

SAC（i）=

Scell（i）（k-1）+tan（Δα2）（L2f-y2k-1），rw≤Lf≤L2；

Scell（i）（k-1）+tan（Δα2）（2L-yk-1-Lf）·

（Lf-yk-1），L2Stube（i），Lf≥L-rw。（10）

因此可得：

φ=∑Ni=1SAC（i）Stotal， （11）

式中φ為面積波及效率。

利用式（11）計(jì)算不同時(shí)刻對應(yīng)的波及面積，進(jìn)而可以計(jì)算得到面積波及系數(shù)的值，但是計(jì)算過程中涉及到滲流阻力、不同時(shí)刻的前緣位置及前緣含水飽和度值等參數(shù)的計(jì)算。

計(jì)算滲流阻力：

R（i）=1K（i）[∫LfrwdξA（ξ）（Krwμw+Kroμo）+

∫L-rwLfdξA（ξ）Kro（Swc）μo] ，（12）

式中：Swc為束縛水飽和度；μw為水相黏度，mPa·s；μo為油相黏度，mPa；Krw為水相相對滲透率；Kro為油相相對滲透率；rw為井筒半徑，m。

計(jì)算各根流管分配到的注入量：

q（i）=Q1R（i）∑Ni=11R（i），（13）

式中Q為整個(gè)單元的總注入量，m3。

計(jì)算水驅(qū)前緣位置：

Lf=r2w+f′w（Swf）Witan（Δα2）h，rw≤Lf≤L2；

L-L2-C，L2

1）根據(jù)流管法的基本原理，建立水驅(qū)油藏的五點(diǎn)井網(wǎng)三角流管模型，可以將平面的二維兩相驅(qū)替問題轉(zhuǎn)換為一維兩相驅(qū)替問題，這樣就從計(jì)算方法上大大降低了求解難度。

2）在沒有流管突破時(shí)，注入水在油藏中均勻波及，波及系數(shù)隨時(shí)間的增加而線性增大；當(dāng)主流管突破后，油井見水，波及系數(shù)增大的程度越來越慢，直至趨于平緩，表明注入水完全波及整個(gè)井網(wǎng)單元，注水對提高原油采收率的作用已經(jīng)十分微弱，水驅(qū)效率明顯降低。

3）流度比越大，注水開發(fā)油藏的驅(qū)油效果越差，通過提高注入液黏度來降低流度比，可以更有效地提高水驅(qū)開發(fā)的效果，進(jìn)而增加原油采收率。

4）當(dāng)油田進(jìn)行注水驅(qū)替時(shí)，在注入壓力不超過巖石破裂壓力的情況下，可以通過適當(dāng)增大注入量的方式來提高原油采油率，加快油田采收速度。

利用流管法求解油藏動(dòng)態(tài)參數(shù)的方法快速有效，未來會得到廣泛的應(yīng)用。本文考慮的是規(guī)則五點(diǎn)法井網(wǎng)，對于油田現(xiàn)場的復(fù)雜不規(guī)則井網(wǎng)，還需進(jìn)一步研究流線的分布，以便運(yùn)用流管法求解，獲得有關(guān)的油田動(dòng)態(tài)參數(shù)，指導(dǎo)油田后期開發(fā)，有效提高原油采收率。

參考文獻(xiàn)/References：

[1]沈非， 程林松， 黃世軍， 等. 基于流管法的普通稠油水驅(qū)波及系數(shù)計(jì)算方法[J]. 石油鉆采工藝， 2016， 38（5）： 645-649.

SHEN Fei， CHENG Linsong， HUANG Shijun， et al. Calculation of sweep efficiency for water flooding development of conventional heavy oil using the stream-tube method[J]. Oil Drilling and Production Technology， 2016， 38（5）： 645-649.

[2]計(jì)秉玉， 李莉， 王春艷. 低滲透油藏非達(dá)西滲流面積井網(wǎng)產(chǎn)油量計(jì)算方法[J]. 石油學(xué)報(bào)， 2008， 29（2）： 256-261.

JI Bingyu， LI Li， WANG Chunyan. Oil production calculation for areal well pattern of low permeability reservoir with non-Darcy seepage flow[J]. Acta Petrolei Sinica， 2008， 29（2）： 256-261.

[3]HIGGINS R V， LEIGHTON A J. A computer method to calculate two-phase flow in any irregularly bounded porous medium [J]. Journal of Petroleum Technology， 1962，14（6）： 679-683.

[4]MUSKAT M， WYCKOFF R D. A theoretical analysis of water flooding networks[J]. Transactions of the AIME， 1934，107（1）：62-76.

[5]DOYLE R E， WURL T M. Stream channel concept applied to water flood performance calculations for multi well， multi zone， three-component cases [J]. Soc Pet Eng AIME Pap， 1971，spe2653（3）： 373-380.

[6]尹洪軍， 付瑩， 王美楠. 基于流管法低滲透油藏開發(fā)數(shù)值模擬研究[J]. 石油化工高等學(xué)校學(xué)報(bào)， 2015， 28（3）： 61-64.

YIN Hongjun， FU Ying， WANG Meinan. Numerical simulation research for low-permeability reservoirs based on stream tube model [J]. Journal of Petrochemical Universities， 2015， 28（3）： 61-64.

[7]王強(qiáng)， 計(jì)秉玉， 聶俊. 聚合物驅(qū)油過程中不同粘度比情況下波及系數(shù)計(jì)算方法 [J]. 石油與天然氣地質(zhì)， 2014， 35（4）： 551-555.

WANG Qiang， JI Bingyu， NIE Jun. Calculation methods of sweep efficiency under different viscosity ratios for polymer flooding process[J]. Oil and Gas Geology， 2014， 35（4）： 551-555.

[8]NELSON R A. Geologic analysis of naturally fractured reservoirs [M]. Tesas： Goulf Publishing Company， 1985.

[9]牛彥良， 李莉， 韓德金， 等. 低滲透油藏水驅(qū)采收率計(jì)算新方法[J]. 石油學(xué)報(bào)， 2006， 27（2）： 77-79.

NIU Yanliang， LI Li， HAN Dejin， et al. A new method for calculation recovery factor of water drive in low permeability reservoir[J]. Acta Petrolei Sinica， 2006， 27（2）： 77-79.

[10]范江， 張子香，扎拉杜新А Б. 非均質(zhì)油層波及系數(shù)計(jì)算模型[J]. 石油學(xué)報(bào)， 1993， 14（1）： 92-97.

FAN Jiang， ZHANG Zixiang， ZALADUXIN А Б. Calculation model of sweep efficiency for heterogeneous reservoir [J]. Acta Petrolei Sinica， 1993， 14（1）： 92-97.

[11]郭粉轉(zhuǎn)， 唐海， 呂棟梁， 等. 滲流啟動(dòng)壓力梯度對低滲透油田五點(diǎn)井網(wǎng)面積波及效率的影響[J]. 大慶石油學(xué)院學(xué)報(bào)， 2010， 34（3）： 65-68.

GUO Fenzhuan， TANG Hai， LYU Dongliang， et al. Effects of seepage threshold pressure gradient on areal sweep efficiency for five-spot pattern of low permeability reservoirs [J]. Journal of Daqing Petroleum Institute， 2010， 34（3）： 65-68.

[12]朱圣舉， 朱潔， 安小平， 等. 低滲透油藏菱形反九點(diǎn)井網(wǎng)面積波及系數(shù)研究[J]. 重慶科技學(xué)院學(xué)報(bào)（自然科學(xué)版）， 2013， 15（2）： 80-84.

ZHU Shengju， ZHU Jie， AN Xiaoping， et al. Research on areal sweep efficiency for rhombus invert 9-spot areal well pattern of low permeability reservoir [J]. Journal of Chongqing University of Science and Technology （Natural Science Edition）， 2013， 15（2）： 80-84.

[13]喻秋蘭， 唐海， 呂棟梁， 等. 平面非均質(zhì)性對面積波及系數(shù)的影響[J]. 重慶科技學(xué)院學(xué)報(bào)（自然科學(xué)版）， 2011， 13（5）： 59-61.

YU Qiulan， TANG Hai， LYU Dongliang， et al. On the Effect of areal heterogeneous on areal sweep efficiency[J]. Journal of Chongqing University of Science and Technology （Natural Science Edition）， 2011， 13（5）： 59-61.

[14]蒲軍， 劉傳喜， 尚根華. 基于流線積分法的注水井網(wǎng)非穩(wěn)態(tài)產(chǎn)量模型[J]. 西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2013， 15（2）： 80-84.

PU Jun， LIU Chuanxi， SHANG Genhua. Non-steady water flooding production model of areal well pattern based on flow line integral method [J]. Journal of Southwest Petroleum University （Science and Technology Edition）， 2013， 15（2）： 80-84.

[15]DING D Y， FARAH N， BOURBIAUX B， et al. Simulation of matrix-fracture interaction in low permeability fractured unconventional reservoirs[C]//SPE Reservoir Simulation Conference. [S.l.]：[s.n.]， 2017： 1-29.