基于页岩油水两相渗流特性的油井产能模拟研究

孙 鑫，刘礼军，侯树刚，戴彩丽，杜焕福，王春伟

（1.中石化经纬有限公司地质测控技术研究院,山东青岛 266003；2.中国石油大学（华东）石油工程学院,山东青岛 266580；3.成都理工大学能源学院,四川成都 610059）

我国油气对外依存度持续上升，常规油气田进入开发后期，页岩油气等非常规油气资源的重要性日益显著[1-2]。2022 年，中国石油的页岩油产量突破300×104t，页岩油作为油气资源的后起之秀，其规模化开发正加速推进[3]。作为实现页岩油藏规模效益开发的关键技术，水平井和体积压裂技术受页岩复杂孔隙结构和固液相互作用影响[4-5]，目前在页岩油渗流机理、压裂液返排规律及产能预测等方面仍面临诸多问题与挑战。

页岩油藏压裂过程中，高压泵注的压裂液使主裂缝和次生裂缝延伸，形成复杂人工裂缝网络；同时，压裂液会通过渗吸作用进入并滞留在基质中，返排过程中少量排出，影响页岩油后续产能[6-8]。实践表明，不同页岩油藏储层的压裂液返排率与产能差异较大，其返排特征主要受页岩储层中油水两相渗流特性影响，由页岩孔隙结构特征控制。页岩油储层孔隙结构通常具有纳米级孔隙发育、孔径分布范围广的特点，前人基于页岩孔隙内流体流动规律、孔径分布等提出了多种两相相对渗透率计算方法。Wang Jinxun 等人[9]将多孔介质概念化为不同尺寸管道串并联的毛细管模型，考虑孔隙尺寸分布及孔隙形状，推导了储层相渗的经典计算方法。Li Ran 等人[10]综合考虑单孔内两相流动特征和孔隙分形结构，建立了页岩两相相渗计算方法，并分析了孔隙尺寸和结构对相渗特征的影响。Su Yuliang 等人[11]考虑页岩有机和无机孔隙内油水赋存特征，建立了页岩油水两相相渗计算方法。数值模拟是常用的油藏产能评价手段，但体积压裂后的页岩油藏多尺度孔隙和裂缝储渗空间发育[12-13]，传统双重介质模型、局部网格加密等模拟方法具有计算量大、难以处理复杂结构裂缝等局限，同时，页岩油流动受多种机理影响[14-17]，数值模拟难度大。

嵌入式离散裂缝模型是将复杂裂缝几何形态直接嵌入正交背景网格中，简化了裂缝的几何剖分过程，极大地降低了计算量和复杂度[18]。国内外学者综合页岩油渗流机理和嵌入式离散裂缝模型，开展了页岩油井产能模拟分析[19-20]，但目前尚未实现微观页岩油水相渗计算与宏观页岩油井产能的耦合。因此，笔者结合页岩油藏相渗计算方法、嵌入式离散裂缝模型和油水两相渗流数学模型，提出了考虑页岩孔隙结构作用下油水两相渗流特性的页岩油井产能数值模拟方法，分析了体积压裂后页岩油藏压裂液空间分布特征和油井产能，实现了页岩微观油水两相渗流特性与宏观油井产能的一体化评价。

1 数学模型及求解

1.1 页岩油藏相渗计算方法

基于毛细管相渗计算模型，结合实际页岩孔隙形状和孔径分布，建立页岩油藏油水两相相渗计算方程。考虑页岩储层中复杂的孔隙形状，采用三角形毛细管模型表征页岩油藏储层。根据三角形毛细管中油水分布状态，单个毛细管中的油水两相流动规律可表示为：

式中：Qo和Qw分别为毛细管中油相和水相的流量，m3/s；μo和μw分别为油相和水相黏度，Pa·s；L为毛细管长度，m；Δp为施加在毛细管上的压差，Pa；ζ为水与孔隙壁面相互作用的无因次阻力系数（用于表征孔隙表面性质对流体流动的影响）；reff为毛细管有效油相半径，m；G为三角形毛细管的形状因子；rd为油水稳定状态下的界面曲率半径，m。

式中：P为三角形毛细管截面周长，m；rin为毛细管内切圆半径，m。

结合单个毛细管中油水流动规律和页岩孔径分布，可得页岩储层油水相对渗透率计算公式：

式中：n为不同尺寸孔隙总数，其中，1～m为中心含水的孔隙尺寸数量，m+1～n为边缘含水的孔隙尺寸数量；Ak为第k个尺寸的毛细管截面积，m2；fk为第k个尺寸的毛细管所占比例。

因此，已知页岩储层的孔径分布后，便可根据式（5）—式（7）计算出页岩油水相对渗透率

1.2 页岩油藏油水两相渗流数学模型

考虑页岩油藏中的油水两相渗流过程，其基质和裂缝中油水两相流体质量守恒关系可统一表达为连续性方程。

式中：β为o 或w，代表油相或水相；ϕ为孔隙度；ρβ为β相流体的密度，kg/m3；Sβ为β相流体饱和度；qβ为β相流体的源汇项，kg/(m3·s)；vβ为β相流体的渗流速度，m/s。

考虑页岩中流体流动的最小启动压力梯度效应，用非线性渗流模型描述基质中油水两相流动[21]。裂缝中通常不存在启动压力梯度效应，因此采用常规达西定律描述裂缝内的油水两相流动过程：

式中：K为绝对渗透率，m2；Krβ为β相流体的相对渗透率；b为拟启动压力梯度的倒数，(Pa/m)-1；a为非线性渗流凹形曲线段的影响因子；ψβ为β相流体的流动势，Pa；pβ为β相流体压力，Pa；h为深度，m。

1.3 页岩油藏渗流模型求解方法

为了高效求解页岩油藏油水两相流体流动，基于嵌入式离散裂缝模型对体积压裂后页岩油藏中复杂裂缝进行几何离散和网格剖分（见图1）。对于给定的体积压裂页岩油藏模型，采用结构化网格对基质区域进行剖分，将水力压裂缝和天然裂缝网络嵌入至剖分后的结构化网格中，利用结构化网格边界切割裂缝网络，形成离散裂缝网格单元，综合形成页岩油藏数值模拟的网格单元系统。

图1 嵌入式离散裂缝模型示意Fig.1 Embedded discrete fracture model

基于网格单元系统，采用有限体积法对油水两相渗流模型进行数值离散，并推导得到离散方程的残差形式：

式中：dij为单元i和j间的距离，m。

采用牛顿-拉夫森方法求解离散的残差方程。

式中：k为迭代层次；c为主变量向量元素；x为主变量向量，选取油相压力和含水饱和度为主变量。

在每个时间步中，采用上述求解格式进行迭代计算，并更新主变量至残差向量的范数小于设定的允许误差，进入下一个时间步进行计算。

2 页岩油藏压裂液分布及产能模拟分析

2.1 页岩油藏油水相对渗透率计算

选取单峰型孔径分布和双峰型孔径分布2 种典型孔径分布页岩（见图2），在孔隙形状参数相同的基础上，采用页岩油藏相渗计算方法计算油水相对渗透率，结果见图3。由图3 可知，相比于单峰型孔径分布，双峰型孔径分布的孔隙尺寸更大，油相的流动能力更强。因此，双峰型孔径分布的页岩储层油相相对渗透率更大，水相相对渗透率更小。

图2 两种典型页岩孔径分布概率曲线Fig.2 Two typical probability curves for shale pore size distribution

图3 两种孔径分布计算的油水两相相渗曲线Fig.3 Oil-water two-phase relative permeability curves calculated with two pore size distributions

2.2 页岩油藏压裂液分布特征分析

为分析压裂过程中压裂液的分布特征，基于胜利油田某页岩油井地质及压裂设计资料，结合嵌入式离散裂缝模型，建立体积压裂页岩油藏模型（见图4）。该页岩油藏基质孔隙度7.0%，渗透率0.5 μD；水力裂缝开度为4 mm，渗透率为5 D；天然裂缝开度为0.3 mm，渗透率0.1 D；初始油藏压力40 MPa，初始含水饱和度为0.05，油相和水相的黏度分别为0.40 和0.25 mPa·s，压裂液注入量为1.0×104m3，压裂后闷井时间为30 d，油水相对渗透率曲线采用图3中单峰型孔径分布的计算结果。注入压裂液后的页岩油藏基质、裂缝中的压力和含水饱和度分布模拟结果如图5 所示，压裂结束闷井30 d 后的基质、裂缝中的压力和含水饱和度分布则如图6 所示。

图4 体积压裂页岩油藏模型Fig.4 Shale oil reservoir model by volume fracturing

图5 压裂结束时页岩油藏压力和含水饱和度分布模拟结果Fig.5 Simulation results of pressure and water saturation distribution in shale oil reservoir after fracturing

图6 闷井30 d 后页岩油藏压力和含水饱和度分布模拟结果Fig.6 Simulation results of pressure and water saturation distribution in shale oil reservoir after 30 days of shut-in

从图5 和图6 可以看出，压裂液主要进入压裂缝及其周边天然裂缝和基质中，引起近水力裂缝周边区域压力升高，该区域裂缝内含水饱和度显著上升，近水力裂缝基质内含水饱和度有所提升。进入闷井阶段后，裂缝和基质中的压力逐渐向周围区域耗散，近水力裂缝高压区域内的压力逐渐降低。同时，裂缝内的压裂液在毛细管力作用下渗吸进入基质，裂缝内含水饱和度降低，对基质中的原油产生一定的渗吸置换作用。

2.3 页岩油藏压后产能分析

在压裂液注入和闷井的模拟结果基础上，开展页岩油藏压后产能数值模拟，对页岩油藏衰竭开发动用范围和产油量进行评价。衰竭开发1 000 d 后的储层基质和裂缝中的压力和含油饱和度分布如图7 所示。经体积压裂后的页岩油藏裂缝网络发育，衰竭开发过程中油藏动用程度高，天然裂缝发育范围内基本可以动用开发。此外，衰竭开发后裂缝内含水量低，但基质内的含水饱和度分布与生产前差异不大，其原因在于压裂液在毛管力作用下滞留在基质中，生产压差难以克服毛管阻力，这也解释了实际页岩储层压裂后压裂液返排率低的现象。页岩油藏生产1 000 d 的日产油量和累计产油量曲线如图8 所示。体积压裂页岩油藏衰竭开发产量递减速度快，经1 000 d 开发后油井累计产油量可达61 145 m3。此外，开发过程中累计产水量为3 335 m3，忽略地层水产出，计算得出压裂液反排率为33%，与现场实际基本符合。由此可见，受页岩油水两相渗流特性及毛管力作用影响，页岩储层压裂后压裂液返排率较低，但体积压裂后的页岩油藏动用程度较好。

图7 生产1 000 d 后页岩油藏压力和含水饱和度分布模拟结果Fig.7 Simulation results of pressure and water saturation distribution in shale oil reservoir after 1 000 days of production

图8 生产 1 000 d 后页岩油藏日产油量和累积产油量曲线Fig.8 Daily oil production and cumulative oil production curves of shale oil reservoir after 1 000 days of production

3 结论

1）考虑页岩孔隙结构作用下油水两相渗流特性，压裂页岩油藏产能数值模拟方法可实现页岩油藏油水相对渗透率、压裂液分布和返排以及油井产能的全流程评价。

2）基于页岩储层孔径分布以及孔隙结构参数，采用毛细管模型可得到页岩油藏油水相对渗透率，不同页岩孔径分布下油水相对渗透率存在较大差异。

3）压裂液主要分布于压裂缝及其周边的天然裂缝和基质中，闷井阶段进入周边基质，基质毛管阻力作用导致压裂液返排率较低，但体积压裂后的页岩油藏动用范围和程度较好。