致密油藏动态裂缝扩展机理及应用1)

邸士莹 程时清 白文鹏 魏 操 汪 洋 秦佳正

* (中国石油大学(北京)石油工程教育部重点实验室,北京 102249)

† (西南石油大学油气藏地质及开发工程国家重点实验室,成都 610500)

引 言

致密油藏水平井体积压裂投产后,随着地层能量不断释放,采油能力衰减,产量递减快.在注水吞吐、水驱、CO2驱等能量补充方式中,长庆、吐哈等油田采用注水吞吐,初期增油效果好.注水吞吐成为接替衰竭开采、补充地层能量的有效方式[1-2].Hagoort 等[3]、Fan 等[4]和Wang 等[5]提出注水易诱导部分天然裂缝扩展.田虓丰等[6]和蒲春生等[7]根据储层应力阴影效应,从数值模拟角度分析了水力裂缝同步扩展到多裂缝的动态延伸特征.Rountree 等[8]、王勃等[9]和唐巨鹏等[10]认为周期注水应力改造可以减小水平地应力差对裂缝扩展的约束能力,促进微裂缝大量发育、沟通天然裂缝形成缝网.王友净等[11]、汪洋等[12]、严谨等[13]、Wang 等[14]注意到天然裂缝扩展形成的不规则复杂缝网极大的增加了焖井期间裂缝与基质的接触面积,注入水通过缝网与基质发生静态渗吸.杨正明等[15-16]、Gao 等[17]和汪勇等[18]发现裂缝性致密储层与普通油藏不同的渗流机理使得逆向渗吸、正向渗吸共同作用置换原油.蔡建超和郁伯铭[19]、许建红和马丽丽[20]明确了注水吞吐初期主要是逆向渗吸作用置换采油.程志林等[21]和Ghasemi 等[22]研究了致密砂岩油水系统自发渗吸特征及规律,通过T2谱反映不同边界条件对渗吸采收率的影响.随着吞吐轮次的增加,地层能量不断释放,导致每一轮次的采油能力都在衰减,又因为致密油藏注水波及范围小,逆向渗吸减弱,注水吞吐效果变差,产量递减快.王家禄等[23]认为在一定的驱替速度范围内,油水黏度比越小,动态渗吸效果越好.徐中一等[24-25]从油藏尺度及实例分析致密储层中注水吞吐效果,以期利用逆向渗吸作用提高开发效果,但效果不明显.王向阳等[26]、Liu 等[27]、马剑等[28]和李帅等[29]建立了不同尺度岩心渗吸物理模拟实验方法,研究了致密储集层渗吸过程的影响因素,并构建了水驱油时渗吸作用的定量评价方法,提出了渗吸及驱替采油的新思想.计秉玉等[30]为更有效地增大驱替作用,针对高含水油井提出改变液流方向的开发方式,提高整体水驱油效率.吴忠宝等[31]和李伟才[32]采用缩小井距及油水井互换的思想提高水驱动用储量,应用于低渗油藏大港油田某区块,初步实施效果显著.

目前裂缝性致密油藏动态裂缝延伸规律、逆向渗吸及驱替机理不清楚,注水吞吐转变开发方式尚未成熟.笔者根据艾尔文理论、弹性力学及渗流力学分析裂缝尖端附近的应力场分布,结合注水诱导裂缝扩展机理及致密油藏发生逆向渗吸原理,提出将注水吞吐转为不稳定脉冲注水开发方式.利用逆向渗吸及线性驱替开发方式,模拟实例井生产10年,预测累计采油、压力及剩余油分布情况,探讨了此方法的可行性.

1 注水诱导裂缝扩展机理

致密油藏大规模压裂水平井注水吞吐后,注水期间井底压力高于裂缝扩展压力,天然微细缝被激发、扩展及延伸,或者天然微细裂缝的充填物被冲刷.相互沟通的天然微细缝形成了新的渗流空间,储层渗透率增加,导流能力增加.

1.1 天然裂缝扩展机理

注水诱导裂缝扩展之前天然裂缝已经存在,第i次裂缝扩展压力计算不同于第一次.艾尔文理论能够用裂缝尖端附近的应力场强度分析第i次裂缝扩展问题,该理论是艾尔文等在格里菲斯脆性断裂理论基础上经过不断丰富和发展而形成的,主要适用于小范围裂缝尖端应力与应变场分析[33-34].在注水后原本闭合的天然裂缝发生扩展,扩展的天然微细裂缝可以视为小范围裂缝,分为3 种类型:I 型裂缝张开、II 型裂缝划开、III 型裂缝撕开(图1).注水诱发裂缝扩展的形式可视为I 型、II 型及交叉裂缝扩展.由于II 型模拟难度系数比较大,本文只对I 型及交叉扩展型开展研究.

图1 裂缝3 种类型Fig.1 Three types of fractures

1.1.1 I 型裂缝扩展规律

假设闭合的天然裂缝只发生I 型破坏,采用Westergaard 应力函数求解,根据拉氏变换和广义胡克定律得到解析函数,进而求裂缝尖端附近的应变分量及应力场强度因子KI.

致密油藏在开发扰动之前就已经承受应力的作用,这个应力一般称为原地应力(或地应力).原地应力状态由上覆岩层压力和两个水平方向的主应力组成,x方向为最大水平主应力方向,y方向为最小水平主应力方向,z方向为铅垂方向(图2(a)).

对于二维平面裂缝,假设“无限大”板,板上有一个长为2a的矩形中心贯穿裂纹(图2(b)),这个板在无限远处受双向等值拉伸应力的作用,忽略流动过程中流体剪切强度的影响、注液系统内各种摩擦力影响、流体与岩石表面间的摩擦阻力，假设在任一时刻,缝内各点的压力均相等.针对平面裂缝体,为研究方便取极坐标系,坐标原点取在裂缝尖端(图2(c)).其中,σ为板两端承受的均匀拉应力,MPa;a为贯穿的矩形切口的半长轴,m;r为裂缝尖端控制区极半径,m;θ为以裂缝尖端为圆心区域上的方位角,(°).

图2 应力分布示意图Fig.2 Schematic diagram of stress distribution

基于裂缝尖端应力场强度的观点,运用线弹性理论和复变函数理论可以求得裂缝尖端附近任意一点A(r,θ)处的应力分布如式(1).裂缝尖端附近的应力场分布[35-36]形式为

其中,fij(θ)为角分布函数;K为裂缝尖端应力场强度因子,其国际单位为MPa∙m1/2.应力场强度因子的临界值即断裂韧度,它反映了裂缝扩展即抵抗脆断的能力[37].针对不同的地层,K值可以通过实验测得.

根据艾尔文理论及裂缝尖端的扩展准则[38],K表达式为

其中,Pnet(t)裂缝张开的净内压,MPa;L为裂缝半长,m.将式(2)代入式(1)得到裂缝尖端附近的应力场与净内压的关系表达式为

由于储层应力场计算过程较繁琐,本节仅利用式(3)定性地分析裂缝尖端的应力与裂缝面受内压的关系.注水初期,裂缝净内压随着注水量的增加而升高,由式(2)可知,应力场强度因子也随之升高.当应力场强度因子K达到断裂韧度KIC,在裂缝尖端会发生扩展.断裂韧度可以实验获得或者查询已有成果[39].

该简化模型有一定的局限性,适用于单一裂缝的判据,且对Ⅰ型裂缝问题比较适用.KI是表征裂缝尖端区域应力场强弱程度的参量,而且是唯一的参量,局限于裂缝尖端附近区域的应力场分析.

1.1.2 多裂缝交叉裂缝扩展渗流模型

致密油藏注水吞吐后,地层中因存在天然裂缝及非均质等因素,还会发生裂缝交叉扩展(图3(a)和图3(b)).交叉扩展的天然裂缝沟通压裂裂缝形成了动态缝网(图3(c)).但由于交叉裂缝扩展从断裂力学分析难度系数较大,因此,从渗流力学的角开展研究.

图3 多裂缝交叉扩展形成动态缝网Fig.3 Multi-fractures cross and expand to form a dynamic fracture network

裂缝发生交叉扩展,假设油、水两相流动且等温,考虑地层岩石及流体的可压缩性,建立注水吞吐油水两相连续性方程、渗流模型以及传导率计算方程.

(1) 注水吞吐油水两相连续性方程

对于注水吞吐油水两相渗流,注入水和原油产出导致储层含水及含油饱和度是动态变化的,则油相及水相的连续方程表示为

其中,po及pw分别为油相和水相压力,MPa;ρo及ρw分别为油层条件下原油密度和地层水的密度,103t/m3;kro和krw分别为油相和水相对渗透率,10−3μm2;k是基质绝对渗透率,10−3μm2;μo及μw分别为油和水的黏度,cp;qo及qw分别为产液量和注入量,m3;so及sw分别为含油饱和度和含水饱和度,%.

随着注水吞吐开发,油层流体的采出,地层压力不断下降,导致岩石发生变形、流体被压缩.基质及流体的状态方程为

其中,po,0为油相初始压力,MPa;pw,0为水相初始压力,MPa;ρo,0为油层条件下原油初始密度,103kg/m3;ρw,0为油层条件下地层水初始密度,t/m3;φ 为油层条件下孔隙度,无量纲;φo为原始油层孔隙度,无量纲;co为原油压缩系数,MPa−1;cr为岩石压缩系数,MPa−1;cw为地层水压缩系数,MPa−1;p0为原始地层压力,MPa;p为地层压力,MPa.

(2)渗流模型

致密油藏基质与裂缝的接触面主要发生逆向渗吸,毛管力占主导作用.则毛管力辅助方程为

由于致密油藏基质渗透率极低,原油从基质和天然裂缝以及压裂裂缝流动到水平井筒中,其产量方程为

其中,pc为毛管力,MPa−1;pwf为井底流压,MPa;WIo及WIw为井指数公式.采用Peaceman 模型[40-41]概念并予以修正,其井指数计算公式如下

其中,re为供给半径,m;rw为井半径,m;h为基质中裂缝高度,m;Bo及Bw分别为原油及水的体积系数,无因次;s为表皮系数,无因次.

由于致密油藏储层及流体的复杂性,考虑启动压力,渗流速度方程为

其中v表示渗流速度,m/s;∇p为地层压力变化量,无因次;k为储层渗透率,10−3μm2,μ为流体黏度,cp;G为启动压力梯度,MPa/m.

基于储层存在的应力敏感效应,根据基质和裂缝的敏感程度设置模型

其中,k是有效渗透率,10−3μm2;ko为初始压力下基质的渗透率,10−3μm2;αk为应力敏感系数,无因次.

考虑裂缝导流系数的动态变化,得到计算公式

其中,kf(t)是动态裂缝渗透率,10−3μm2;wf(t)表示动态裂缝宽度,m.

(3)传导率计算方程

针对交叉裂缝的情况,将交叉裂缝单元之间的传导率近似处理,将裂缝网格可看作普通网格,此时裂缝网格与基质网格传导率表示为

式中Ti为裂缝网格与基质网格传导率,kfi是裂缝和基质渗透率的调和平均值,10−3μm2;Afi是裂缝与基质接触面积,m2;dfi是两条裂缝中心距离,m.

当两条裂缝相交时,裂缝网格与裂缝网格传导率[42-43]为

其中,kf1和kf2是裂缝段的渗透率,10−3μm2;df1和df2是裂缝段到相交线的平均距离,Af1及Af2是裂缝段接触面积,m2.

建立压裂水平井M56-152H 常规数值模型及交叉裂缝数值模型,常规模型模拟日产液与实际数据有差异(图4(a)),交叉裂缝扩展形成的动态缝网模型模拟日产液,历史拟合较好(图4(b)).

图4 M56-152H 产液量历史拟合结果Fig.4 M56-152H Production history matching results

致密油藏在不同区域均能使用相同形式的运动方程构建统一的控制方程.因此,可以对模型进行离散化求解,降低耦合问题的复杂性,使模型求解变得简单统一.

下面研究主要以I 型扩展为主,推导裂缝扩展长度、刻画动态缝网渗流规律.

1.2 裂缝扩展长度

注水诱发裂缝扩展后,裂缝扩展长度及裂缝半长是主要的关键参数.根据艾尔文裂缝扩展原理,当K>KIC,裂缝尖端压力达到地层破裂压力,天然裂缝发生扩展.由图2(b)模型假设裂缝扩展缝宽为矩形,第i次扩展长度Li,长度增量为ΔLi(图5).

图5 裂缝扩展长度示意图Fig.5 Schematic diagram of crack propagation length

假设10 级压裂水平井存在60 条天然裂缝,天然裂缝扩展瞬间,根据格里菲斯准则采用断裂力学的能量平衡原理[44-46],从旧裂缝生成且还未充填流体的瞬间开始,到新裂缝产生且还未充填流体的瞬间结束为单元,对整个系统进行线弹性求解为

其中,ΔEi表示地面流入的能量,ΔEl表示滤失的能量,ΔEf表示缝内弹性能量,ΔErs表示基质内弹性能,ΔEic表示新的裂缝的能量

可根据格里菲斯准则,克服基质内储存的弹性能

产生新裂缝表面所吸收的能量Erci+1为

裂缝扩展的整个系统中能量平衡方程为

整理式(26)得到裂缝扩展长度的计算公式

其中,Q为注入量,m3;P(t)注入压力,MPa;Li为第i次裂缝扩展的裂缝半长,m;PP地层孔隙压力,MPa;ρ注入流体的密度,kg/m3;PH井筒液柱产生的压力,MPa;P为裂缝扩展压力,MPa;Ti为扩展周期,d;Ai+1第i+1 次扩展中半裂缝的侧面积,m2;σh地层最小水平主应力,MPa;E地层岩石弹性模量,MPa;ω为裂缝宽度,m;H油层厚度,m;γ 表面能密度,MPa/m2;CL扩展后的滤失系数,m/d1/2;Pio第i次裂缝产生时,地面压力,MPa;Ti为扩展周期,d;Ti0扩展周期初始时间,d;C(t)滤失系数,无因次.

根据某致密油藏M56 块注水吞吐期间的5 几口井的注入量数据,计算得到滤失系数C(t)为

随着注水时间延长,扩展的天然裂缝延伸,开发过程中需将缝长控制在合理的范围,防止井间裂缝水窜.由式(2)可知储层中裂缝尖端方位角θ是随机分布的,天然裂缝会扩展、延伸并相互沟通,呈现不规则复杂缝网(图6(b)),此时储层中流体的渗流机理比较复杂.因此,揭示致密油藏动态裂缝的渗流规律是转变开发方式的关键.

图6 注水诱发天然裂缝扩展Fig.6 Water-induced the expansion of natural fractures

1.3 裂缝扩展过程中的逆向渗吸判别

自发渗吸是低渗透裂缝水驱油藏的重要开采机理[47].衰竭开采及注水吞吐都会有大量注入水通过不规则复杂缝网进入地层,与基质大面积接触,发生正向及逆向渗吸作用.为了确定两种渗吸主控作用,根据渗吸物理机理的判别标准[48],计算NB−1(重力与毛管力的比值Bond 数的倒数)的数值.利用NB−1数值的范围判定基质渗吸过程中毛细管力和重力的贡献.当NB−1>5 时,毛管力支配渗吸过程,与基质的接触面主要发生逆向渗吸;当NB−1<0.2 时,重力支配渗吸过程,主要发生正向渗吸;当0.2

式中,B为结构常数,无因次;Δρ为油水密度差,kg/m3;g为重力加速度,取9.8 m/s2;H为多孔介质高度,m;K为渗透率,10−3μm2;Φ为孔隙度,无因次;δ为表面张力,mN/m.式(29)没有考虑基质的润湿性,而润湿性对基质渗吸有重要的影响.Iffly 等[49]对上述渗吸机理判别方程进行修正,刘卫东等[50]和姚同玉等[51]也做过类似研究,并做了大量实验.改进后的NB−1的表达式为

式中,α为润湿接触角,(°).Standnes 等[52]进一步实验研究了润湿性对渗吸机理的影响.润湿性越弱,润湿接触角越大,界面张力越大.当接触角大于90°时,界面张力成为渗流阻力.其中B的数值是由孔隙结构常数的表达式[53]计算得出,该公式为

式中,τ是流体流线的迂曲度,无因次;Df为孔隙分形维数,且0

式中,β为拟合常数.裂缝性致密油藏具有低孔、低渗的特征,致密储层渗吸过程毛管力比常规油藏大.某致密油藏M56 块,孔隙度φ为0.165,β取值0.51.在天然裂缝扩展后形成的动态缝网中,裂缝视为海岸线科赫曲线类型,根据文献[54],Df取值为1.2618,由式(31)计算得到B1= 0.45.动态缝网中的基质视为Sierpinski 垫片和Sierpinski 地毯型,Df分别取值1.585 及1.892,计算得到B2=0.394 及B3= 0.22.本文的结构系数近似取三类型的平均值,即B≈ 0.4.

某致密油藏M56 块,Δρ=1.9 × 102kg/m3,重力加速度g取9.8 m/s2;H=20 m,K=1×10−5μm2,Φ= 16.5%,δ= 30 mN/m,α= 33.6°.由式(30)计算得到致密油藏M56 块NB−1= 35.7.

根据Schechter 提出的判别标准,当NB−1>5 时,毛管力支配渗吸过程,与基质的接触面主要发生逆向渗吸.在毛管力作用下,致密油藏M56 块发生大规模的逆向渗吸,将原油置换到复杂缝网中,这为注水吞吐转变为脉冲注水开发方式提供了大量可驱替的原油(图7).

图7 逆向渗吸作用Fig.7 Reverse imbibition

注水初期,毛管力是逆向渗吸的动力.随着注水吞吐轮次的增加,地层能量不断释放,注水波及面积范围降低,导致岩石润湿性减弱.润湿接触角α变大,毛管力变成了逆向渗吸阻力.逆向渗吸作用减弱,由逆向渗吸作用置换原油数量将减少,注水吞吐效果变差.

2 致密油藏线性驱替方式

裂缝性致密油藏注水吞吐主要发生径向驱替及逆向渗吸作用,多轮次注水吞吐后,渗吸置换作用减弱,产量降低明显.合理利用天然裂缝扩展形成的高导流通道,尽可能地发挥逆向渗吸加驱替作用,是目前裂缝性致密油藏亟待解决的难题.

下面研究将注水吞吐转为脉冲注水的可行性.

2.1 脉冲注水原理

脉冲注水充分发挥注入水在动态裂缝中的线性驱替作用,增加流体的流动,促使注水波及面积迅速增加,逆向渗吸范围增加.脉冲注水与注水吞吐主要不同点在于,注水吞吐只发生在本井附近,多轮次吞吐后,致密油藏注入水推进速度较为缓慢,两井之间部分区域的原油未能被波及.脉冲注水是通过注水井周期性地改变注水量,造成地层压力周期性升高和降低,从而在裂缝和基质之间产生压差.压力升高到高于裂缝开启压力时,天然裂缝扩展、延伸并沟通压裂裂缝而形成的复杂缝网.注入水驱替原油,向前推进,在生产井被采出.

从宏观上分析,注水吞吐在注水阶段驱替作用方向是径向的(图8(a)),脉冲注水主要是线性驱替作用占主导地位,主要沿着井间裂缝线性向前驱替(图8(b)).

图8 径向和线性驱替作用示意图Fig.8 Diagram of radial and linear displacement

从油藏尺度分析,注水吞吐方式在基质中的注水波及主要是拟径向的(图9(a)),脉冲注水波及主要是沿着注水诱导裂缝扩展方向发生线性驱替作用(图9(b)).

图9 裂缝扩展前后径向和线性驱替作用示意图Fig.9 Diagram of radial and linear displacement

动态缝网中注入水受毛细管力作用和水动力效应,加之天然裂缝扩展形成的高导流通道,注水波及面积迅速增加,逆向渗吸范围扩大(图10(a)和图10(b)).逆向渗吸作用加强,缝网与基质之间的流体交换速度和质量发生质变,置换出更多的原油.

图10 注水吞吐和脉冲注水逆向渗吸作用范围对比Fig.10 Comparison of water-injection huff and puff and pulse water injection imbibition range

图10 注水吞吐和脉冲注水逆向渗吸作用范围对比(续)Fig.10 Comparison of water-injection huff and puff and pulse water injection imbibition range (continued)

在下一个脉冲注水期间,裂缝中被置换的原油及无渗吸作用波及区域的原油,在线性驱替作用下流动到井筒.水平井脉冲注水利用了逆向渗吸及线性驱替的采油原理(图11),与注水吞吐开发方式相比极大地提高了驱油效率.下面以某实际区块为例,讨论该开发方式合理的油井工作制度.

图11 脉冲注水逆向渗吸及线性驱替作用Fig.11 Reverse imbibition and linear displacement

2.2 脉冲注水可行性探讨

某致密油藏属于凝灰岩裂缝性致密油藏,是国内外规模较大的致密油藏水平井注水吞吐试验区.该区块天然微细裂缝比较发育,地层压力系数为1.013,井距为100 m～200 m,断裂韧度KIC为0.908 MPa·m1/2.2015年开始注水吞吐,截止2020年9月注水吞吐120 井次,其中56 口井吞吐3～6 轮次.该油藏M56 块为例,模拟实际井注水吞吐转化为脉冲注水的累计采油及平均地层压力.

M56-151H 井(图12(a)) 于2017年5月开始注水吞吐,建立该井数值模型,有关地层参数、初始地层条件等参数等见表1.

图12 M56-151H 井组井位图及裂缝发育情况Fig.12 Well location and fracture development

表1 数值模拟参数表Table 1 Numerical simulation parameter table

模拟结果显示,注水初期裂缝净内压随着注水量的增加而升高,应力场强度因子也随之升高.当应力场强度因子K达到0.908 MPa·m1/2,天然微细缝发生开启及扩展现象.模拟结果还表明,采用定产液方法拟合的实际产液,历史拟合效果较好(图13).

图13 M56-151H 产液量历史拟合结果Fig.13 M56-151H Production history matching results

2.3 注水吞吐转脉冲注水

注水吞吐转为脉冲注水前,在设计井距及段间距离时,要着重考虑裂缝扩展长度,以期确保不会发生水窜.图14(a)中A点至B点为第i次裂缝扩展长度,根据式(20) 及式(27) 计算得到Li扩= 0.67 m,Lif= 7.66.注水120 天,天然裂缝扩展长度为13.33 m(图14(b)).每次天然裂缝扩展,Li扩及Lif均小于区块的井距和段距,水平井段和邻井之间没有发生水窜现象(图15).

图14 模拟裂缝扩展长度结果Fig.14 Simulated fracture propagation length

图15 裂缝扩展长度与井距段距对比Fig.15 Result of simulated fracture propagation length

随着脉冲注水周期的增加,被置换的原油大幅度增多.分析原因,裂缝扩展形成的动态缝网加强了脉冲注水的驱替作用,提高了注入水与基质接触的面积.逆向渗吸范围扩大,逆向渗吸作用加强,置换原油的速度及数量发生了质变.

以M56-152H 井组为例,分别设计3 种方案,注入时间设置为1 d,2 d,3 d.设计注入量设置为300 m3/d,但在模拟含水饱和度场图时发现,水平井间发生明显水窜现象(图16).

图16 注水量300 m3/d 时发生水窜Fig.16 Water breakthrough occurs when the water injection volume is 300 m3/d

经过反复调整参数,当注入量设置为100 m3/d时,井底压力达到裂缝开启压力41.7 MPa,低于地层破裂压力60 MPa.停注时间设置为1 d,2 d,3 d,两口采油井一直以50 m3/d 的采油量生产10年(表2).

表2 脉冲注水3 种方案Table 2 Three schemes of pulse water injection

M56-151H 井3 种设计方案生产10年的累计采油在1.02 × 105～1.05 × 105m3,方案3 累计采油最高,预测10年累计采油量为1.05 × 105m3(图17(a)).该方案的工作制度是注3 天停3 天,注入量100 m3/d,产油量50 m3/d.预测10年后地层压力降至20 MPa,但仍然保持较高的水平(图17(b)).

图17 模拟脉冲注水3 种方案生产10年预测产量Fig.17 Predicted output of the three schemes of simulated pulse water injection in 10 years

模拟M56-152H 井注水吞吐4 轮次后转变为脉冲注水,并按照方案3 开采制度,模拟马56-152H井4 轮次后仍以注水吞吐方式生产至2029年.对比两者的累计采油量,转为脉冲注水预测10年累计采油量高于注水吞吐(图18).

图18 注水吞吐与脉冲注水生产10年预测产量Fig.18 Predicted output of water-injection huff and puff and pulsed water injection production for 10 years

M56-152H 在脉冲注水前,储层中部分天然裂缝处于闭合状态(图19(a),图19(b)).脉冲注水后天然裂缝形态发生了扩展延伸、并沟通了压裂裂缝形成了动态缝网(图19(c),图19(d)).

图19 脉冲注水前后裂缝分布形态Fig.19 Fracture distribution

图19 脉冲注水前后裂缝分布形态 (续)Fig.19 Fracture distribution (continued)

模拟M56-151H 井底剩余油分布,脉冲注水方式的剩余油饱和度较低.说明注水吞吐4 轮次后转为脉冲注水方式生产8年,剩余油充分动用(图20).

图20 注水吞吐与脉冲注水剩余油分布对比Fig.20 Comparison of remaining oil distribution

3 结论

(1)裂缝性致密油藏注水诱导裂缝扩展可视为张开裂缝及交叉扩展型.注水期间,裂缝净内压随着注水量的增加而升高时,应力场强度因子也随之升高.当应力场强度因子K达到断裂韧度KIC,在裂缝尖端会发生扩展、延伸并相互沟通,呈现不规则复杂缝网.提出的裂缝扩展长度计算方法可用于确定裂缝延伸位置,预防裂缝水窜.改进的逆向渗吸的判别公式确定M56 区块动态裂缝发生了大规模的逆向渗吸.

(2)周期性地改变注水量,从而产生裂缝和基质压差,造成地层中压力场的不稳定分布.在压力上升过程中,注入水通过天然裂缝扩展形成的高导流通道进入到基质,极大地增加了注水波及面积.促使逆向渗吸范围扩大,逆向渗吸作用加强,缝网与基质之间的流体交换速度和数量均发生质变.更多的原油被置换到缝网中,在下一个脉冲注水周期的线性驱替作用下流向采油井.

(3)讨论了裂缝性致密油藏水平井注水吞吐转变为脉冲注水方式,即由逆向渗吸与径向驱替作用转变为逆向渗吸及线性驱替作用.脉冲注水累计产油高、逆向渗吸作用强,充分发挥逆向渗吸及线性驱替作用,能够实现有效驱油的目的.

(4)实际注水时,应当控制注水量,使得注入期间井底压力接近天然裂缝开启压力,但要低于地层破裂压力,其目的是增加地层有效渗流能力,同时防止注入水沿裂缝水窜.