致密油藏分段多簇压裂水平井非稳态产能模型

2019-06-25姬靖皓席家辉曾凤凰杨啟桂

岩性油气藏 2019年4期

姬靖皓，席家辉，曾凤凰，杨啟桂

（1.中国石油长庆油田分公司第四采油厂，陕西榆林719000；2.中国石油长庆油田分公司第十二采油厂，甘肃庆阳745000）

0 引言

随着能源需求的不断增大和压裂改造技术的长足进步，致密油作为一种储量丰富的非常规油气［中国致密油资源总量为（106.7～111.5）亿 t］，是目前我国较为现实的石油接替资源［1-3］。致密油藏由于低孔低渗，其开发实践表明，常规压裂改造很难满足工业生产需要［4-6］，目前通常采用水平井钻井结合分段多簇射孔压裂的方式形成多条高导流能力裂缝，并利用缝间干扰产生复杂裂缝网络，其大小用储层改造体积（SRV）表示，改善近井地带储层的渗透性，增大泄流面积，从而提高单井产能［7-9］。

国内外众多学者为了研究考虑SRV的压裂水平井生产动态，主要应用数值模拟方法和解析/半解析方法。数值模拟方法［10-12］可灵活处理复杂裂缝形态，但很难考虑非常规储层的渗流机理，而解析/半解析方法因其严格的数学推导和较高的计算效率成为目前的研究热点，其中Brown等［13］提出的三线性流模型和Stalgorova等［14］提出的五线性流模型因其简单有效的特点而被众多学者采用。姚军等［15］针对低渗油藏非达西渗流特征，在三线性流模型中引入了启动压力梯度，分析了低渗透油藏压裂水平井压力动态特征。苏玉亮等［16］利用三线性流模型研究了致密油藏水平井多级压裂布缝策略，分析了SRV与累积产油量的关系。赵二猛等［17］基于双重介质模型，建立了窜流系数随着距主裂缝距离逐渐变化的致密油体积压裂水平井线性流模型。苏玉亮等［18］、姬靖皓等［19］、JI等［20］结合致密储层低速非达西流动和裂缝应力敏感特征建立了致密油藏体积压裂水平井五区复合流动模型，分析了其相应的压力动态和产能变化特征。为了考虑段间和簇内渗流规律的差异性，郭肖等［21］建立了页岩气藏分段多簇压裂水平井复合线性流模型，然而，目前还没有适合致密油藏分段多簇压裂水平井的产能模型。因此，基于前人研究成果，结合致密储层低速非达西流动、裂缝应力敏感特征，并考虑段间渗流场和簇间渗流场差异性，建立致密油藏分段多簇压裂水平井非稳态产能模型，以期为致密油藏分段多簇压裂水平井产能评价和裂缝参数优化提供理论依据。

1 物理模型及假设

致密油藏经水平井分段多簇压裂改造后，每个压裂段内会形成多簇高导流能力裂缝，簇间裂缝由于缝间干扰形成复杂裂缝网络，经过简化的物理模型如图1所示。图1中：由于分段多簇压裂形成的裂缝区别于多段单簇压裂缝，段间渗流场改造程度弱，簇间渗流场改造程度强，储层渗透性得到改善，2个渗流场具有明显的差异性，不能简单地应用三线性流模型或五线性流模型来进行表征。如图2所示，段间渗流场为每段改造体间的渗流区域，可划分为7个区域：区域1、区域1 c、区域2为未改造区，可采用单重介质模型并引入启动压力梯度来表征其低速非达西流动特征；区域3为部分受效区（改造程度弱），区域4、区域4 c为有效改造区（改造程度强），这3个区域可利用双重介质模型来进行表征，并考虑裂缝渗透率应力敏感性；区域F为压裂裂缝区，考虑裂缝渗透率应力敏感性。簇间渗流场为每段改造体内部的渗流区域，可划分为3个区域：区域1为未改造区；区域2为有效改造区；区域F为人工裂缝区。簇间渗流场各区域表征方法同段间渗流场。图1和图2中：LR为油藏长度，m；WR（2 xe）为油藏宽度，m；ns为压裂段数；nc为每段压裂簇数；Δ yc为簇间距，m；Δ yi为段间距，m；Δ ye为两端裂缝到储层边界的距离之和，m；l为段裂缝间改造体积的延伸长度，m；xF为裂缝半长，m；wF为裂缝宽度，m。

图1 分段多簇压裂裂缝与改造体积示意图Fig.1 Schematic diagram of segmented multi-cluster fractures and SRV

图2 段间渗流场和簇间渗流场示意图Fig.2 Seepage field schematic diagram of segments and clusters

物理模型的其他假设条件为：①油藏外边界封闭，压裂水平井位于油藏中心以定井底流压生产；②各裂缝性质相同且每条裂缝都完全穿透储层；③单相微可压缩液体等温渗流，忽略重力、毛管压力和井筒中的流动阻力；④流体流动为线性流动，由未改造区流向改造区进而从压裂裂缝流入水平井筒。

2 产能模型的建立

首先定义相关无因次变量如表1所列，Pi为原始地层压力，MPa；Pwf为井底流压，MPa；ηref为参考导压系数，μm2/（mPa·s·MPa-1）；qF为单条裂缝产量，m3/d；t为时间，d；μ 为流体黏度，mPa·s；B 为原油体积系数；Kref为参考渗透率，D；h为储层厚度，m；α为形状因子；Km为基质渗透率，D；Kfi为原始状态下的裂缝渗透率，D；Lref为参考长度，m；Ф为孔隙度，%；Ct为综合压缩系数，MPa-1；KFi为原始状态下的人工裂缝渗透率，D；CL为流体压缩系数，MPa-1；λm为启动压力梯度，MPa/m；γ为渗透率模量，MPa-1。其中的下标，j代表各个渗流区域；i代表原始状态；m，f代表双重介质的基质和裂缝系统；F代表人工压裂裂缝；ref代表参考值；D代表无因次。

表1 无因次变量的定义Table 1 Dimensionless variables

2.1 簇间渗流场模型

对于簇间渗流场，选取1/4裂缝控制区域进行研究时可利用三线性流模型进行表征，此时区域F的无因次化渗流数学模型可表示为

为了消除式（1）的非线性，可采用 Pedrosa［22］变换式如下：

将式（2）带入式（1），并对式（1）进行 Pedrosa［22］变换。通常γD非常小，因此零阶摄动即可满足计算精度的要求。对变换后的式子取零阶摄动，而后进行 Laplace变换（t→ u）得到

据姚军等［15］的求解过程，可求得本文模型单条簇裂缝在拉氏空间的产量解为

2.2 段间渗流场模型

对于段间渗流场，选取1/2裂缝控制区域进行研究，区域F在拉氏空间中的无因次化渗流数学模型可表示为

段间渗流场的左半边类似于一个三线性流模型，因此据姚军等［15］的求解过程，区域4 c的解可表示为

由物理模型假设可知，段间渗流场的区域4 c与簇间渗流场的区域2都属于有效改造区，因此，式（6）中的 β4c和K4fi等于式（4）中的 β2和K2fi。段间渗流场的右半边类似于一个五线性流模型，根据姬靖皓等［19］、JI等［20］的研究，区域 4 的解可表示为

2.3 多裂缝叠加处理

在得到上述段间渗流场和簇间渗流场的单条裂缝解后，利用Stehfest［23］数值反演方法得到其时空间的解，而后采用Meyer等［24］的方法进行产量叠加处理，分别定义不同位置裂缝所对应的形状因子：端部裂缝形状参数为δe=WR/Δ ye，段间裂缝形状参数为 δi=WR/Δ yi，簇间裂缝形状参数为 δc=WR/Δ yc。则分段多簇压裂水平井的产量解可表示为

3 产能动态特征分析

3.1 矿场实例验证

以某致密油藏1口分段多簇压裂水平井基本参数为例，原始地层压力为15 MPa，井底流压为5 MPa，水平井长度为1 400 m，8段2簇压裂（簇间距15 m），每段改造体两端裂缝各向外延展的改造体宽度为20 m，可根据图3进行SRV大小的计算。裂缝半长为100 m，裂缝导流能力为5 D·cm，油藏厚度为20 m，基质渗透率为0.1 mD，孔隙度为14%，原油黏度为 4 mPa·s，体积系数为 1.2，改造区基质和裂缝压缩系数分别为0.000 75 MPa-1和0.001 MPa-1，启动压力梯度为0.005 MPa/m，裂缝渗透率模量为0.02 MPa-1。

由图4可以看出，油井生产初期存在压裂液反排现象，产量波动大，中后期模型计算结果与实际生产数据拟合程度较高，进一步证明了本文所推导的非稳态产能模型具有一定的准确性。

图3 改造体积计算方法示意图Fig.3 Schematic diagram of calculating SRV

3.2 低速非达西流动的影响

低速非达西流动的程度可利用启动压力梯度值的大小来进行表征［25-26］，启动压力梯度越大，储层低速非达西流动越明显，其对油井产能动态特征的影响越大。如图5所示，启动压力梯度主要影响油井中后期的产能，这是因为分段多簇压裂会在近井地带形成复杂裂缝网络，有效提高了储层的渗透性，从而减小了非达西流动对油井前期产能的影响。随启动压力压力梯度的增大，油井中后期的产量越小，累积产油量也越小。

图4 计算数据与实际数据对比图Fig.4 Comparison of calculated data and actual production data

图5 启动压力梯度对油井产能的影响Fig.5 Effect of threshold pressure gradient on oil well productivity

3.3 裂缝渗透率应力敏感的影响

图6表明了裂缝渗透率应力敏感对致密油藏分段多簇压裂水平井产能动态的影响。由于主裂缝和SRV区域的裂缝是流体主要的渗流通道，裂缝应力敏感性的存在会导致其渗透性随开发过程中地层压力的减小而降低，因此，裂缝的应力敏感性影响油井整个生产阶段。应力敏感系数越大，油井的产油量和累积产油量越低。

3.4 裂缝条数的影响

现场分段多簇压裂每段簇数一般为2～4簇。分别计算22种布缝形式下（4/6/8/10/12/14/16段2簇缝共7种；4/5/6/7/8/9/10/11段3簇缝共8种；3/4/5/6/7/8/9段4簇缝共7种）油井生产5 a的累积产油量。从图7可以看出，相同裂缝分布形式下，裂缝条数的增加可有效提高油井的产能，但累积产油量的增幅会随着裂缝条数的增加增幅逐渐变缓。不同布缝方式对比可以看出，当裂缝总数较少时，裂缝的段簇比对累积产油量影响较大，而当裂缝条数较多时，段簇比对累积产油量影响不明显。相同裂缝条数下，段簇比越大，累积产油量越大，这是因为压裂后油井生产效率很大程度上取决于单条裂缝的储层改造范围，簇数越多，簇间干扰越严重。

图6 裂缝渗透应力敏感性对油井产能的影响Fig.6 Effect of fracture stress sensitivity on oil well productivity

图7 裂缝条数对油井产能的影响Fig.7 Effect of fracture number on oil well productivity

3.5 SRV的影响

图8反映了22种裂缝分布方案对应的5 a累积产油量和SRV的关系，随SRV增大，油井累积产油量增幅逐渐变缓，因此现场压裂并不能一味地追求较高的SRV。这是因为，SRV较小时，裂缝条数较少，不同裂缝段之间的区域向改造体的线性流供液时间越长，而当SRV增大时，对应的裂缝条数增多，缝间干扰加剧，累积产油量增幅变缓。因此，矿场实际生产时为了提高生产效率，除了增大裂缝条数提高SRV外，不能使裂缝分布过密，导致缝间干扰而降低开发效果。