杨子由,陈民锋,屈 丹,杨金欣
(1.油气资源与探测国家重点实验室中国石油大学(北京),北京 102249;2.中国石油大学(北京)石油工程学院,北京 102249)
在特低渗透油藏中,流体在多孔介质的渗流过程中往往存在启动压力梯度,且由于高渗透条带的发育,通常使得油藏渗透率在不同方向差异化分布,导致特低渗透油藏出现渗透率各向异性特征。高渗透条带的存在使得水驱开发过程中出现单方向水窜,因此,常规井网形式不再适合此类油藏的开发[1-5]。
在对各向异性油藏的渗流理论研究中,往往使用流管模型对个别渗流方向上的渗流规律进行分析[6-10],研究渗透率在不同方向上的差异对整个注采井网系统渗流规律的影响较少;在储量动用规律研究中,多基于油藏注水波及情况及油井产能等指标进行井网调整[11-17],较少在油田注采井网模型的基础上,对特低渗透各向异性油藏井网加密方式进行适应性评价,并以此来论证井网加密方式对此类油藏的合理性。为此,基于实际油田注采井网形式,考虑油藏各向异性的渗流特点,得到不同注采井网条件下渗流场分布规律,论证排状加密井网对各向异性油藏的适应性,建立储量动用评价方法,以研究井网加密后的储量动用规律。
油藏的非均质性一般包括层间、层内、平面和微观非均质性,其中平面上渗透率的各向异性最为普遍,直接影响水驱波及系数和驱油效率。对于重点解决油藏各向异性问题的注采井网设计来说,为保证各方向上的均衡驱替,需要重点考虑渗透率在不同位置及方向上的差异。
油藏在投入开发时,由于各类资料相对较少,油藏的各向异性特征和主渗流方向不甚明确,基础井网系统往往是正方形反九点井网。注水开发过程中,油藏水驱动态越来越表现出单方向水窜的特点,反映了油藏渗流非均匀性对水驱开发的不利影响。
由于正方形反九点井网越来越表现出不适应地下水驱流动的非均匀性,后续开发过程中力图通过井网加密改变注水井网形式,使调整后的井网形式与各向异性油藏的渗流特点相适应,以最大可能地提高油藏水驱波及系数,进而提高注采井间储量动用程度。
若均质油藏中有1 口井,其产液量或注入量为q,则油藏平面上任意点处的势为:
若油藏中有多口井同时注入或生产,根据叠加原理,假定各投产井的坐标为(xi,yi),则任意点处(x,y)的势为:
若为各向异性油藏,任意点处的势为:
其中:
假定Kx=nKy(n=1 时为均质油藏),将(4)式代入(3)式得:
根据势可确定驱动压力梯度。压力梯度为矢量,其沿x与y轴正方向的分量为:
在CZ特低渗透各向异性油藏注水开发过程中,针对沿主渗流方向水驱突进的问题,一般将初始的正方形反九点基础井网整体加密为交错排状井网,使注水井排方向与油藏主渗流方向一致,以改善水驱效果。初始注采井网为井距与排距相等的正方形反九点井网(图1a);在初始反九点井网的两井排中间,生产井交错位置处加密生产井,初始井网中的生产井转为注水井,形成排状加密井网(图1b)。
图1 加密前后注采井网形式Fig.1 Well patterns before and after infilling
基于各向异性油藏渗流方程,可得到2 种井网形式下注采井间的势和驱动压力梯度分布。
许某某是程瀚的老朋友,2006年,许某某在合肥注册成立安徽某置业投资公司。2008年年底,许某某在程瀚帮助下,其女朋友杨某某顺利进入了交警支队工作。2012年年底,其公司出现纠纷,便找到了程瀚,程瀚问包河公安分局打招呼。2009年6月份的一天,许某某请求程瀚向合肥市公安局车管所打招呼,为其办理一套好的车牌照。朋友的请托之事,程瀚动动嘴就能轻松地予以解决。来而不往非礼也,程瀚把他的商人朋友们当成“提款机”。
初始反九点井网 在初始反九点井网中取1个正方形注采单元(图1a中深黄区域)作为研究单元,以注水井作为坐标原点,考虑配产配注原则,设注水井的注入量为qw;各角生产井的产量为qoA,其坐标为(xoAi,yoAi);各边生产井的产量为qoE,其坐标为(xoEi,yoEi)。
根据势叠加原理,可以得到研究单元中任意点M(x,y)处的势为:
根据(6)式,可以得到任意点M处的驱动压力梯度为:
排状加密井网 在排状加密井网中(图1b)取1个注采单元作为研究区域。以中心生产井为坐标原点,考虑配产配注原则,设中心生产井的产量为qo;各注水井的注入量为qw,其坐标为(xwti,ywti);各角生产井的产量为qoAt,其坐标为(xoAti,yoAti);各边生产井的产量为qoEt,坐标为(xoEti,yoEti)。同理,可以得到排状加密井网任意点M处的势和驱动压力梯度分别为:
驱动压力梯度可反映流体在油藏中流动的难易程度,通过统计注采井间驱动压力梯度的分布,定量表征注采井间的储量动用情况。为此,提出无因次动用范围与无因次动用提高程度2个概念。
无因次动用范围 在注采井间,设研究范围的总面积为RT。在此范围内,若驱动压力梯度大于某值的面积为R,则该条件下无因次动用范围定义为:
无因次动用提高程度 井网加密后,若某位置处的驱动压力梯度更大,则视该位置处的储量动用程度提高。在研究范围内,设储量动用程度提高范围的总面积为NT,在此范围内,若井网加密后驱动压力梯度增量大于某值的面积为N,定义无因次动用提高程度为:
在特低渗透油藏中,若某区域的驱动压力梯度小于油藏启动压力梯度,流体不流动,此区域称为难动用区,反之则为储量有效动用区。由于多孔介质在微观尺度上存在差异,不同特低渗透油藏的启动压力梯度不同,判别储量有效动用的标准也不同。依据CZ特低渗透各向异性油藏的物性,将驱动压力梯度小于0.015 MPa/m 的区域定义为难动用区,驱动压力梯度为0.015~0.03 MPa/m 的区域定义为低-中动用区,驱动压力梯度大于0.03 MPa/m 的区域定义为中-高动用区。
特低渗透油藏井网加密的目的是缩小注采井间难动用区范围,扩大低-中动用区范围,进而提高全油田范围内的储量动用程度。
根据CZ 特低渗透各向异性油藏开发实践及调整经验,设置各向异性油藏在不同方向上渗透率的关系为Kx=5Ky,此条件下可以明显体现不同方向上渗透率的差异。研究范围为各自井网形式下2个相邻的注采单元(图1中深黄和浅黄区域)。各井定产液量或注入量,初始注采井距为200~600 m。
排状井网的线性排驱和注采沿线 在排状井网中,注水井排和生产井排交替排列,在注水井排附近形成高压(势)区,生产井排附近形成低压(势)区。若油藏为各向异性油藏,且井排布置方向与高渗透方向一致时,注水井附近的高压区交汇,形成高压水线。该高压水线的延伸方向与注水井排一致,并与生产井排附近形成低压带交替分布,在注采井间形成线状排驱。
在注水井排上,以不同位置为起点,分别向相邻生产井排作垂线,将此类垂线定义为注采沿线。当各向异性油藏使用排状井网进行开发时,注水井排不同位置为起点的诸多注采沿线上,其势分布规律相近,说明不同注采沿线上流体的流动状态相似,体现各向异性油藏在排状井网形式下的开发特点。由此可知,对于各向异性油藏,线状排驱可以使注入流体沿低渗透方向整体排驱,有效提高低渗透方向的储量动用程度,进而提高储量动用程度。
井网加密前后注采井间势分布规律 根据(7)式和(13)式,可分别得到均质油藏和各向异性油藏在井网加密前后的势分布(初始井距为400 m)。由图2可以看出:①在反九点井网中,均质油藏等势线呈圆形分布,而各向异性油藏中势场呈椭圆形分布,二者形态差异反映了流体在各自油藏条件下的渗流特征。②对于各向异性油藏,井网加密后,注水井排附近高势带呈现线状分布,形成高压水线,与生产井排处的低势带交错分布,形成线状排驱。③相比于均质油藏,各向异性油藏在井网加密后,注采井间势等值线近似线状分布,体现注水井排整体排驱的特性;井排间等值线分布更密集,说明沿低渗透方向势梯度更大,注入流体更多地沿低渗透方向驱替,发挥线状排驱优势。
注采沿线势分布规律 注采主沿线定义为以注水井为起点的注采沿线,注采非主沿线定义为起点非注水井的注采沿线。为使论证更具代表性和充分性,取注采主沿线、距注水井1/4 井距和1/2 井距处为起点的非主沿线作为研究对象。井网加密前后不同油藏注采沿线上势分布规律(图3)表明:①在注采主沿线上,近井区域(距井约40 m 内的区域)势呈漏斗状下降趋势,后渐变为近线性下降的趋势;在非主沿线上,势随距注水井排距离的增加呈相对平缓下降的趋势。②在均质油藏中,3 条注采沿线中任意2 条在确定距离条件下,势始终存在一定的差值,该差值在超过近井区域范围内基本趋于一致,反映了注采井间径向驱动的特点。③在各向异性油藏中,注采主沿线、非主沿线上,势的变化趋势在超过近井区域基本趋于一致,且势的大小接近,说明在此条件下形成了注水井排的高压水线,势分布等值线横向轮廓几乎与井排相互平行,产生整体向生产井排线性排驱的态势。
综合分析图2 和图3 后发现:在各向异性油藏中,将初始反九点井网转化为排状加密井网,可以在注采井间形成线性排驱;不同注采沿线上势相近的变化规律,表明注入水可以对油藏中流体沿低渗透方向(注采沿线方向)整体性排驱,使该方向的水驱效果更为均衡。因此,排状加密井网更适合各向异性油藏的注采开发。
图2 井网加密前后不同油藏中的势分布Fig.2 Potential distribution in different reservoirs before and after well pattern infilling
注采井间驱动压力梯度分布 基于(12)式和(14)式,可以得到在确定条件下井网加密前后2 个相邻注采单元的驱动压力梯度分布。由初始井距为400 m时的驱动压力梯度分布(图4)可以看出:①井网加密后,在注采井间的较大范围内,驱动压力梯度有一定程度的增大,有效缩小了井网加密前难动用区范围,储量动用程度有所提高。②动用提高范围主要集中在加密井附近,由此可知井网加密对该区域储量动用程度提升效果明显,而两相邻注水井之间区域动用效果较差,后期可采取其他措施以改善该区域的储量动用效果。
图3 井网加密前后注采沿线势分布规律Fig.3 Potential distribution along injector and producer line before and after well pattern infilling
图4 注采井间驱动压力梯度分布Fig.4 Driving pressure gradient distribution between injectors and producers
不同条件下储量动用变化规律 根据(21)式和(22)式,可得到在不同初始井距条件下,井网加密前后的无因次动用范围和无因次动用提高程度。由无因次动用范围曲线(图5a,5b)可以看出:①当井网形式一定时,相同驱动压力梯度下,随着井距的减小,无因次动用范围均增加,储量动用程度有所提高;在任意曲线上,随驱动压力梯度的升高,均出现曲线斜率突增的拐点(图中虚线与各曲线交点)。在注采井间,小于拐点处驱动压力梯度的区域,可以看作是此条件下注采井间主要动用区域。②相比于初始井网,井网加密后,相同驱动压力梯度下的无因次动用范围更大,说明井网加密后,注采井间的储量动用效果更好。③井网加密前,对于井距大于400 m 的曲线,其拐点处对应的驱动压力梯度较小(<0.015 MPa/m),此时注采井间的主要动用区域为难动用区,说明在此条件下,注采井间主要动用范围的驱动压力梯度较低,初始井距越大,其井网加密后储量动用提高潜力越大。
根据图5c 可以看出:①当驱动压力梯度一定时,井距越小,动用提高区域的范围越大,储量动用程度越高。②在无因次动用提高程度曲线中,亦有拐点的出现。在井距一定的条件下,小于此拐点以及大于反九点井网无因次动用曲线(图5a)中拐点对应的驱动压力梯度范围,即为该条件下井网加密主要提升的动用范围。③当井距大于400 m 时,井网加密主要缩小难动用区范围(<0.015 MPa/m),提升低-中动用区(0.015~0.03 MPa/m)动用范围,而井距在小于400 m 的条件下对中-高动用区提升效果较好。
图5 无因次动用范围和无因次动用提高程度变化曲线Fig.5 Dimensionless producing range and dimensionless incremental percentage of producing range
基于典型油藏注采井网的加密方式,建立反映油藏各向异性的注采井间渗流场计算方法,可以得到不同条件下注采井间渗流场分布。
基于特低渗透各向异性油藏有效动用条件,建立了定量分析此类油藏井网加密前后储量动用变化规律的评价方法。
加密后的排状井网更适合各向异性油藏的水驱开发,注采井间的储量有效动用程度显著提高;当初始反九点井网井距大于400 m 时,井网加密对储量动用提升效果更为明显。
符号解释
q——井产液量或注入量,m3/d;ϕ——势,无量纲;r——任意点与井之间的距离,m;C——常数;xi,yi——第i口井横、纵坐标,m;i——井编号;m——井总数,口;ϕi——第i口井在任意点产生的势,无量纲;qi——第i口井产液量或注入量,m3/d;x,y——任意点的横、纵坐标,m;Ke——油藏平均渗透率,mD;Kx,Ky——油藏横、纵向渗透率,mD;n——渗透率级差,无量纲;p——地层中任意处的压力,MPa;μo——地层原油黏度,mPa·s;qw——注水井注入量,m3/d;qoA——初始井网中各角生产井的产量,m3/d;xoAi,yoAi——初始井网中各角生产井的横、纵坐标,m;qoE——初始井网中各边生产井的产量,m3/d;xoEi,yoEi——初始井网中各边生产井的横、纵坐标,m;ϕM——任意点M处的势,无量纲;xai,yai——初始井网中任意点距各角生产井在x和y方向的距离,m;xbi,ybi——初始井网中任意点距各边生产井在x和y方向的距离,m;qo——排状加密井网中中心生产井产量,m3/d;xwti,ywti——排状加密井网中各注水井的横、纵坐标,m;qoAt——排状加密井网中各角生产井的产量,m3/d;xoAti,yoAti——排状加密井网中各角生产井的横、纵坐标,m;qoEt——排状加密井网中各边生产井的产量,m3/d;xoEti,yoEti——排状加密井网中各边生产井的横、纵坐标,m;xci,yci——排状加密井网中任意点距各注水井在x和y方向的距离,m;xdi,ydi——排状加密井网中任意点距各角生产井在x和y方向的距离,m;xei,yei——排状加密井网中任意点距各边生产井在x和y方向的距离,m;RT——研究区域的总面积,m2;R——研究区域内驱动压力梯度高于某值时的面积,m2;RD——无因次动用范围,无量纲;ND——无因次动用提高程度,无量纲;N——研究范围内驱动压力梯度增量高于某值时的面积,m2;NT——研究范围内储量动用程度提高区域的总面积,m2。