CN115935843A - 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法 - Google Patents

一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法 Download PDF

Info

Publication number
CN115935843A
CN115935843A CN202211578528.8A CN202211578528A CN115935843A CN 115935843 A CN115935843 A CN 115935843A CN 202211578528 A CN202211578528 A CN 202211578528A CN 115935843 A CN115935843 A CN 115935843A
Authority
CN
China
Prior art keywords
fracture
seepage
model
water
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211578528.8A
Other languages
English (en)
Inventor
朱争
党海龙
王强
康胜松
李留杰
梁卫卫
白璞
凃兴平
崔鹏兴
王继超
石立华
张天龙
李东
王敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Yanchang Petroleum Group Co Ltd
Original Assignee
Shaanxi Yanchang Petroleum Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaanxi Yanchang Petroleum Group Co Ltd filed Critical Shaanxi Yanchang Petroleum Group Co Ltd
Priority to CN202211578528.8A priority Critical patent/CN115935843A/zh
Priority to CN202310366098.1A priority patent/CN116579201A/zh
Publication of CN115935843A publication Critical patent/CN115935843A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，步骤如下：步骤1：建立模型；在油水两相渗流的基本方程基础上，考虑渗吸驱替对油水渗流的双重作用及注水诱导动态裂缝，分别将驱替和渗吸过程的毛管力模型、相对渗透率模型及基质裂缝窜流量带入到渗流控制方程，建立渗吸驱替双重作用下注水诱导动态裂缝油水两相渗流数学模型；步骤2：模型求解，采用全隐式有限差分方法获得考虑渗吸驱替双重作用的动态裂缝油水两相渗流场数据。本发明更能精确表征储层实际渗流环境，开发指标预测精度更高。

Description

一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法
技术领域
本发明涉及一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法。
背景技术
目前的研究普遍聚焦在静态裂缝与基质的渗吸作用，带压动态裂缝与基质渗吸双重作用(毛管力与相对渗透率渗吸驱替)下的渗流机理的研究较少。同时对动态裂缝的表征目前只是通过改变近井地带导流能力来等效裂缝，而没有考虑在注水过程中动态裂缝的几何尺度开启及延伸特征。在渗流尺度上，基于有效应力与地层压力，简单地考虑裂缝渗透率随压力的变化。目前鲜有考虑动态裂缝与渗吸耦合作用的渗流模型，导致注水诱导裂缝扩展过程中渗吸置换的渗流规律尚不明确，无法进一步指导现场开发实践。针对目前相关数值模型存在的问题，有必要研究一种同时考虑动态扩展裂缝及其与基质渗吸作用、嵌入双重毛管力曲线建立渗流数学模型来精准表征诱导裂缝发生时流体的渗流规律及其对开发效果的影响规律。
发明内容
本发明旨在针对现有技术存在问题，提出一种同时考虑渗吸驱替双重作用的注水诱导动态裂缝油水两相的渗流数学模型。
本发明的技术方案在于：
一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，步骤如下：
步骤1：建立模型；
在油水两相渗流的基本方程基础上，考虑渗吸驱替对油水渗流的双重作用及注水诱导动态裂缝，分别将驱替和渗吸过程的毛管力模型、相对渗透率模型及基质裂缝窜流量带入到渗流控制方程，建立渗吸驱替双重作用下注水诱导动态裂缝油水两相渗流数学模型；
步骤2：模型求解，采用全隐式有限差分方法获得考虑渗吸驱替双重作用的动态裂缝油水两相渗流场数据。
其中，建立模型的具体过程如下：
(1)建立油水两相渗流的基本微分方程；
在考虑产量项的情况下，油水两相渗流的基本微分方程组为
Figure BDA0003989942780000021
Figure BDA0003989942780000022
(2)建立带压渗吸双重毛管力模型
考虑驱替过程的毛管力计算模型为
Figure BDA0003989942780000023
式中：Pcq为考虑驱替过程的毛管压力，MPa；
Pd为驱替入口压力，MPa；
Sw为含水饱和度，无量纲；Swr为水相对渗透率趋近于0时的含水饱和度，
无量纲；λ为孔隙大小分布指数，无量纲；
考虑渗吸过程的毛管力计算模型为
Figure BDA0003989942780000024
式中：Pcs为考虑渗吸过程的毛管压力，MPa；
Pneg为渗吸曲线上，水饱和度达到最大值时的最大毛管压力，MPa；
m、n、α为模型拟合参数；
Siw为束缚水饱和度，无量纲；Sor为残余油饱和度，无量纲；
考虑驱替过程的相渗模型为
Figure BDA0003989942780000025
Figure BDA0003989942780000026
式中：Krwq为考虑驱替过程的水相相对渗透率，无量纲；
Kroq为考虑驱替过程的油相相对渗透率，无量纲；
考虑渗吸过程的相渗模型为:
Figure BDA0003989942780000031
Figure BDA0003989942780000032
式中：Kros为考虑渗吸过程的水相相对渗透率，无量纲；
Krws为考虑渗吸过程的油相相对渗透率，无量纲；
在渗流基本模型基础上，考虑渗吸驱替对油水渗流的双重作用，分别将驱替和渗吸过程的毛管力模型和相对渗透率模型，代入到渗流控制方程即可；
(3)建立油水井人工压裂裂缝数值刻画模型EDFM
嵌入式离散裂缝数值模型的控制方程为：
Figure BDA0003989942780000033
Figure BDA0003989942780000034
式中：
Figure BDA0003989942780000035
为基质孔隙度，无量纲；
Figure BDA0003989942780000036
为裂缝孔隙度，无量纲；
Sl,m为基质中l相饱和度，无量纲；Sl,f为裂缝中l相饱和度，无量纲；
km为基质渗透率，mD；Kf为裂缝渗透率，mD；
krl,m为基质中l相的相对渗透率，无量纲；krl,f为裂缝中l相的相对渗透率，无量纲；
l为油相或者水相，o、w；μl为l相的粘度，mpa.s；ρlg为l相的重度，MPa/m；Z为深度，m；
pl,m为基质中l相的目前地层压力，MPa；pl,f为裂缝中l相的目前地层压力，MPa；
qm为油藏基质与嵌入裂缝间的窜流量的窜流量，m3/d；
qf为嵌入裂缝与油藏基质间的窜流量的窜流量，m3/d；
qm＝-qf；
其中，qf的计算过程包括四种情况；
四种情况分别为三个类型非相邻连接NNC及裂缝与井眼之间的传导率系数；三个类型非相邻连接NNC定义如下所示：
①NNC类型Ⅰ：裂缝段与其穿过的基质网格间的连接；
②NNC类型Ⅱ：某一裂缝内相邻裂缝段间的连接；
③NNC类型Ⅲ：相交的不同裂缝段之间的连接；
对某一类NNC所连接的两个网格间的裂缝窜流量表示为：
qf＝λlTNNCΔP(15)
式中：
λl为l相流度，mD·mPa-1·s-1；
TNNC为NNC对传导率系数，mD·m；
ΔP为相邻网格压力差值，MPa；
Figure BDA0003989942780000041
式中：kNNC为NNC对渗透率，mD；
ANNC为NNC对接触面积，m2；
dNNC为与NNC连接对特征距离，m；
①NNC类型Ⅰ
NNC类型Ⅰ的基质与嵌入网格裂缝传导率系数为：
Figure BDA0003989942780000042
式中：Tmf为基质与嵌入网格裂缝传导率系数，mD·m；
Amf为基质与嵌入网格裂缝的接触面积，m2；
kmf为基质渗透率和裂缝渗透率的调和平均，mD；
dmf为基质到裂缝平面的平均法向距离，m；
Figure BDA0003989942780000051
基质到裂缝平面的平均法向距离dmf为：
Figure BDA0003989942780000052
式中：dv为网格块内体积微元，m3；
xn为该体积微元到裂缝法向距离，m；
V为网格块体积，m3；
②NNC类型Ⅱ
NNC类型Ⅱ相邻裂缝段间传导率系数为：
Figure BDA0003989942780000053
Figure BDA0003989942780000054
Figure BDA0003989942780000055
式中：
Figure BDA0003989942780000056
为相邻裂缝段间传导率系数，mD·m；
kf1为相邻裂缝段1的渗透率，mD；wf1为相邻裂缝段1的开度，m；
kf2为相邻裂缝段2的渗透率，mD；wf2为相邻裂缝段2的开度，m；
Lint为裂缝段交线的长度，m；
df1为相邻裂缝段1中心到交线的平均距离，m；
df2为相邻裂缝段2中心到交线的平均距离，m；
相交裂缝计算式为：
Figure BDA0003989942780000057
Figure BDA0003989942780000058
式中：dSi为面积微元；
Si为裂缝段i面积，m2；xn为面积微元到交线距离，m；
③NNC类型Ⅲ
NNC类型Ⅲ的相交裂缝段间传导率系数为：
Figure BDA0003989942780000061
Figure BDA0003989942780000062
Figure BDA0003989942780000063
式中：
Figure BDA0003989942780000064
为相交裂缝段间传导率系数，mD·m；
Ac为两个相交裂缝段的公共面，m2；
kf1′为相交裂缝段1的渗透率，mD；kf2′为相交裂缝段1的渗透率，mD；
dseg1为相交裂缝段1形心到公共面的距离，m；
dseg2为相交裂缝段2形心到公共面的距离，m；
④裂缝与井眼之间的传导率系数
基于Peaceman数学模型推导出裂缝与井眼之间的传导率系数及裂缝网格的等效井系数分别为：
Figure BDA0003989942780000065
Figure BDA0003989942780000066
式中：WIf为裂缝与井眼之间的传导率系数，mD·m；
re为裂缝网格的等效井系数，m；rw为井眼的半径，m；
Δθ为径向井包含在裂缝内的圆心角角度，即Δθ＝2π；
Lf、hf分别为裂缝的长度及高度，m；
考虑注水诱发动态裂缝的问题，引入动态的裂缝参数
Figure BDA0003989942780000067
及kf′：
其中，动态裂缝渗透率kf′为
Figure BDA0003989942780000068
式中：kfi为初始裂缝渗透率，mD；
γf裂缝渗透率应力敏感模量，MPa-1；
pfi初始裂缝内压力，MPa；pf裂缝内压力，MPa；
kmi为基质渗透率，mD；
其中，动态裂缝孔隙度
Figure BDA0003989942780000071
为
Figure BDA0003989942780000072
式中：
Figure BDA0003989942780000073
为裂缝孔隙度，无量纲；
Figure BDA0003989942780000074
为初始裂缝孔隙度，无量纲；
Cf为裂缝压缩系数，MPa-1；
pf为裂缝内压力，MPa；
另外，模型的边界条件和初始条件为
1)边界条件
①外边界条件
定压力梯度外边界条件：
Figure BDA0003989942780000075
式中：
Figure BDA0003989942780000076
为边界方向上的法向压力梯度，MPa/m；
const为常数；
②内边界条件
定井底压力：
P|rw＝Pwf(x,y,z)＝const      (35)
式中：P|rw为定井底压力，MPa；
Pwf为井底流压，MPa；const为常数；
2)初始条件
假设数值模拟初始时刻油藏的压力为某一已知函数P0(x,y,z)，油藏数值模拟压力的初始条件为
P(x,y,z,0)＝P0(x,y,z)(36)在开发过程中，油藏中油水两相发生流动时，需要确定开发初始时刻的某一相的饱和度场，含水饱和度的初始条件为：
Sw＝Sw0(x,y,z)(37)
式中：Sw为油藏各介质***中任意一点含水饱和度的参数，无量纲；
P为油藏各介质***中任意一点压力的参数，MPa。
其中，所述模型求解的具体过程为：
在对渗流控制方程进行差分后，对离散模型进行数值求解计算，计算流程如下：(1)网格划分和初始化：对基质进行网格划分并对基质和裂缝物性进行初始化，并计算传导率；
(2)裂缝网格划分、嵌入：根据裂缝网格坐标，将裂缝嵌入至基质网格中；
(3)计算裂缝宏观渗流参数：根据裂缝不同状态分别计算裂缝的渗透率、孔隙度、相对渗透率等参数；
(4)计算窜流量：分别计算裂缝-基质，裂缝-裂缝的传导率和窜流量，构建线性方程组系数矩阵；
(5)求解渗流微分方程：将计算得到的窜流量带入渗流微分方程中，全隐式进行求解计算；
(6)输出结果：计算至最终时刻，输出压力和饱和度场。
本发明的技术效果在于：
本发明在通过考虑动态诱导裂缝扩展及其与基质的渗吸作用，同时利用嵌入吸吮、驱替的两条毛管力曲线，建立了低渗透油藏诱导动态裂缝渗流数学模型。模型在几何尺度上将裂缝的动态开启与闭合过程等效成裂缝网格随时间的激活和湮灭过程，进而在渗流尺度上通过引入裂缝动态孔隙度和动态渗透率，并对建立渗流模型进行离散以及数值求解，揭示了产生诱导裂缝时渗吸作用规律及其对开发效果的影响特征，从而在油藏尺度上实现了注水诱导缝与渗吸耦合作用机制的数值表征，更能精确表征储层实际渗流环境，开发指标预测精度更高。
附图说明
图1为基质-裂缝的窜流量计算示意图。
图2为嵌入式离散裂缝动态表征过程。
图3为动态裂缝计算流程图。
图4为注水诱导形成动态裂缝过程图。
图5为岩心有效压力与渗透率关系示意图。
图6为岩心有效压力与渗透率关系示意图。
图7为不同裂缝尺寸岩心采出程度随注入量的变化关系曲线。
图8为考虑诱导裂缝及双重毛管力作用的渗流规律模拟计算生产制度。
图9为基质压力场云图。
图10为基质含水饱和度场云图。
图11为是否考虑渗吸作用的井底压力随时间变化图。
图12为是否考虑渗吸作用的产水速率和累计产水量随时间变化图。
图13为是否考虑带压渗吸双重作用的产油速率和累计产油量随时间变化图。
具体实施方式
实施例1—模型验证(注水诱导裂缝起裂实验结果对比)
(1)数值模型构建
为了验证上述构建的注水诱导裂缝模型准确性，基于低渗油藏注水诱导裂缝起裂实验，构建实验尺度下的岩心数值模型，分别对岩心的注水诱导裂缝实验过程进行模拟。
构建的数值模型，模型网格为22×22×33。模型孔隙度为7.9％，渗透率为0.845mD，原油粘度为1mPa·s，压力根据实验所得突破压力设置为4MPa；水的总矿化度为30000mg/L，模型温度设置为25℃，其他参数设置参考黑油模型默认值，因此在其他属性设置时考虑该因素，对模型参数设置做出相应的小幅度调整。在进行岩样设置时，将岩样设置为方岩心，岩样尺寸为2.5cm×2.5cm×5cm，为后续对注入及产出乘以相应的系数以匹配实验室的数据条件。岩样前端设置一口注水井，进行裂缝模拟；
(2)注水诱导裂缝起裂实验数值模拟验证
数值模型形成动态裂缝的过程如图4所示，裂缝随着注水的进行而不断扩展，最终形成贯穿缝；
根据数值模型的模拟结果，可得到注入压力与渗透率的关系如图5所示；
分别与岩心的注水诱导裂缝实验进行拟合对比，可以发现上述数值模型在注水时压力曲线同实验压力曲线一样表现出频繁的波动，出现多个突破特征，前者突破压力在2.8MPa～5Mpa之间，岩心第一次突破；
对于模型，注水憋压后注入压力急剧上升，在有效压力为3.0MPa～3.3MPa时接近峰值，压力达到裂缝开启压力值但还未开启，随着注水的进行，当模拟到2.3min时，注入压力增大到4.5MPa，达到模型突破压力，产生微裂缝。在达到突破压力后，压力有个释放降低的过程，后续随着水注入量与注入时长增加，压力波动上升，裂缝继续扩展，渗透率随之上升，最终在10min裂缝延伸至岩心表面，此时注入压力保持平稳状态。即上述数注水诱导裂缝起裂实验数值模拟结果与实验所得结论能够精准拟合，再一次证明了压力与诱导裂缝开启的关系；
(3)商业软件与数值模性精确度验证
为了验证本文建立的注水诱导裂缝模型较现有商业数值模拟软件在动态裂缝模拟方面的优越性，本节利用商业软件(Tnavigator)建立裂缝岩心数值模型。进行注水诱导裂缝起裂实验模拟，可得到如图6所示的压力图版。通过拟合对比，对于Tnavigator所建模型，在进行注水诱导裂缝模拟时，只能产生静态裂缝，因此无法精确描述裂缝开启，产生微裂缝的时机状态等。静态裂缝同上述模型所建动态裂缝相比，准确度不够高，只能简单模拟注水生产裂缝的过程，无法刻画裂缝的动态扩展，因此上文所建数值模型可作为验证实验准确性的模型。
实施例2—模型验证(带压动态渗吸实验结果对比)
(1)数值模型构建
为验证上述构建考虑渗吸机理的注水诱导裂缝模型准确性，基于低渗储层带压动态渗吸与裂缝耦合作用影响因素实验，在上述所建的数值模型的基础上，构建岩心数值模拟，将实验结果与数值模拟结果进行对比从而验证数值模型在带压动态渗吸的准确性。数值模型保留原有的网格属性，设置模型1孔隙度为11.9％，渗透率为0.421×10-3μm2；设置模型2地层孔隙度为10.9％，渗透率为0.381×10-3μm2；设置模型3孔隙度为11.3％，渗透率为0.357×10-3μm2；设置岩样4孔隙度为11.6％，渗透率为0.244×10-3μm2。统一设置地层温度为50℃，岩样前端注水介质选择盐水，地层水总矿化度为78000mg/L。设置模型1为基质，对模型2、3、4分别进行裂缝模拟，设置裂缝大小为1/3裂缝、2/3裂缝和贯通缝，其他参数不变，4种数值模型。模型构建完成后进行生产模拟，前端为注水井，后端为生产井。当生产井含水率达到98％停止注水。随后焖井12h，继续开始注水，记录不同时刻的采出程度与注水量的关系。
(2)裂缝岩心动态渗析实验数值模拟验证
不同裂缝尺寸岩心采出程度随注入量的变化关系的数值模型结果与实验结果如图7所示。与裂缝岩心动态渗析实验结果进行拟合对比后可以看出，动态渗析的数值模型的结果与实验结果拟合度非常高。从数模结果来看，基质模型1在焖井前采出程度为48.68％，焖井后最终采出程度为50.65％，焖井渗吸采出程度提高了1.97％。设置1/3裂缝的模型2在焖井前采出程度为56.4％，焖井后最终采出程度为59.93％，焖井渗吸采出程度提高了3.53％。设置2/3裂缝的模型3在焖井前采出程度为61.58％，焖井后最终采出程度为67.15％，焖井渗吸采出程度提高了5.57％。设置贯通缝的模型4在焖井前采出程度为40％，焖井后最终采出程度为45.89％，焖井渗吸采出程度提高了5.89％。
该模型验证了实验所得到的结论，即裂缝可以提高驱替采出程度，且采出程度随着裂缝的增加而增加，但裂缝贯穿后，形成了高渗通道，驱替的采出程度会降低，焖井渗吸对采出程度的贡献随着裂缝的增长而增加，这是由于模型越长，其与注入水接触面积远大，增加了渗吸作用。
具体实例
在考虑渗吸、驱替两条毛管力曲线基础上，探究是否考虑渗吸作用对渗流场和生产动态的影响，建立了考虑渗吸作用和不考虑渗吸作用的数值模型算例，裂缝设置为动态裂缝；生产制度以及油藏数值模拟相关参数与前面一致，工作制度为注水吞吐，即定压生产-注水-闷井这三个工作制度不断循环，生产制度如下图8所示。
通过对比基质压力场(图9)可以发现，考虑和不考虑渗吸作用对压力的影响较大，考虑渗吸作用以后，水更多的进入基质中，使基质维持较高的压力，因此考虑渗吸作用注水后，对储层有补充能量的作用。
通过对比基质饱和度场的云图可以发现(图10)，考虑和不考虑渗吸作用对含水饱和度的影响。这是因为考虑渗吸作用以后，由于毛管力的原因，注入的水可以进入基质的更深处，因此含水饱和度较高。
图11为是否考虑带压渗吸双重作用的井底压力随时间变化图，通过分析井底压力曲线可以发现，注水时井底压力逐渐上升，闷井后由于注入的水逐渐进入基质，因此井底压力下降。对于是否考虑渗吸作用二者的井底压力差距不明显，是因为注入压力远大于渗吸的毛管力，因此渗吸对井底压力的作用不明显；
图12为是否考虑渗吸作用的产水速率和累计产水量随时间变化，第一次定压生产时，未考虑渗吸作用的产水速率高于考虑渗吸作用，这是由于在第一次开采时，考虑渗吸作用的毛管力为开采时的阻力，因此产水量小。对于整个生产过程的累计产水量而言，未考虑渗吸的累计产水量是高于考虑渗吸作用的累计产水量。
图13为是否考虑渗吸作用的产油速率和累计产油量随时间变化，通过对比是否考虑渗吸作用的产油速度和累计产油量可以发现：随着生产的持续，无论是否考虑渗吸作用，二者的产油速度下降，但通过注水蓄能作用，使产油速度下降速度降低。

Claims (3)

1.一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：步骤如下：
步骤1：建立模型；
在油水两相渗流的基本方程基础上，考虑渗吸驱替对油水渗流的双重作用及注水诱导动态裂缝，分别将驱替和渗吸过程的毛管力模型、相对渗透率模型及基质裂缝窜流量带入到渗流控制方程，建立渗吸驱替双重作用下注水诱导动态裂缝油水两相渗流数学模型；
步骤2：模型求解，采用全隐式有限差分方法获得考虑渗吸驱替双重作用的动态裂缝油水两相渗流场数据。
2.根据权利要求1所述考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：所述建立模型的具体过程如下：
(1)建立油水两相渗流的基本微分方程；
在考虑产量项的情况下，油水两相渗流的基本微分方程组为
Figure FDA0003989942770000011
Figure FDA0003989942770000012
(2)建立带压渗吸双重毛管力模型
考虑驱替过程的毛管力计算模型为
Figure FDA0003989942770000013
式中：Pcq为考虑驱替过程的毛管压力，MPa；
Pd为驱替入口压力，MPa；
Sw为含水饱和度，无量纲；Swr为水相对渗透率趋近于0时的含水饱和度，无量纲；λ为孔隙大小分布指数，无量纲；
考虑渗吸过程的毛管力计算模型为
Figure FDA0003989942770000014
式中：Pcs为考虑渗吸过程的毛管压力，MPa；
Pneg为渗吸曲线上，水饱和度达到最大值时的最大毛管压力，MPa；
m、n、α为模型拟合参数；
Siw为束缚水饱和度，无量纲；Sor为残余油饱和度，无量纲；
考虑驱替过程的相渗模型为
Figure FDA0003989942770000021
Figure FDA0003989942770000022
式中：Krwq为考虑驱替过程的水相相对渗透率，无量纲；
roq为考虑驱替过程的油相相对渗透率，无量纲；
考虑渗吸过程的相渗模型为：
Figure FDA0003989942770000023
Figure FDA0003989942770000024
式中：ros为考虑渗吸过程的水相相对渗透率，无量纲；
Krws为考虑渗吸过程的油相相对渗透率，无量纲；
在渗流基本模型基础上，考虑渗吸驱替对油水渗流的双重作用，分别将驱替和渗吸过程的毛管力模型和相对渗透率模型，代入到渗流控制方程即可；
(3)建立油水井人工压裂裂缝数值刻画模型EDFM
嵌入式离散裂缝数值模型的控制方程为：
Figure FDA0003989942770000025
Figure FDA0003989942770000026
式中：
Figure FDA0003989942770000027
为基质孔隙度，无量纲；
Figure FDA0003989942770000028
为裂缝孔隙度，无量纲；
Sl,m为基质中l相饱和度，无量纲；l,f为裂缝中l相饱和度，无量纲；
m为基质渗透率，mD；f为裂缝渗透率，mD；
krl,为基质中l相的相对渗透率，无量纲；krl,为裂缝中l相的相对渗透率，无量纲；
l为油相或者水相，o、w；μl为l相的粘度，mpa.s；ρlg为l相的重度，MPa/m；
为深度，m；
pl,m为基质中l相的目前地层压力，MPa；pl,f为裂缝中l相的目前地层压力，MPa；
qm为油藏基质与嵌入裂缝间的窜流量的窜流量，m3/d；
qf为嵌入裂缝与油藏基质间的窜流量的窜流量，m3/d；
qm＝-qf；
其中，qf的计算过程包括四种情况；
四种情况分别为三个类型非相邻连接NNC及裂缝与井眼之间的传导率系数；
三个类型非相邻连接NNC定义如下所示：
①NNC类型Ⅰ：裂缝段与其穿过的基质网格间的连接；
②NNC类型Ⅱ：某一裂缝内相邻裂缝段间的连接；
③NNC类型Ⅲ：相交的不同裂缝段之间的连接；
对某一类NNC所连接的两个网格间的裂缝窜流量表示为：
qf＝λlTNNCΔP(15)式中：
λl为l相流度，mD·mPa-1·s-1；
TNNC为NNC对传导率系数，mD·m；
ΔP为相邻网格压力差值，MPa；
Figure FDA0003989942770000031
式中：kNNC为NNC对渗透率，mD；
ANNC为NNC对接触面积，m2；
dNNC为与NNC连接对特征距离，m；
①NNC类型Ⅰ
NNC类型Ⅰ的基质与嵌入网格裂缝传导率系数为：
Figure FDA0003989942770000032
式中：Tmf为基质与嵌入网格裂缝传导率系数，mD·m；
Amf为基质与嵌入网格裂缝的接触面积，m2；
kmf为基质渗透率和裂缝渗透率的调和平均，mD；
dmf为基质到裂缝平面的平均法向距离，m；
Figure FDA0003989942770000033
基质到裂缝平面的平均法向距离dmf为：
Figure FDA0003989942770000041
式中：dv为网格块内体积微元，m3；
xn为该体积微元到裂缝法向距离，m；
V为网格块体积，m3；
②NNC类型Ⅱ
NNC类型Ⅱ相邻裂缝段间传导率系数为：
Figure FDA0003989942770000042
Figure FDA0003989942770000043
Figure FDA0003989942770000044
式中：
Figure FDA0003989942770000045
为相邻裂缝段间传导率系数，mD·m；
kf1为相邻裂缝段1的渗透率，mD；wf1为相邻裂缝段1的开度，m；kf2为相邻裂缝段2的渗透率，mD；wf2为相邻裂缝段2的开度，m；
Lint为裂缝段交线的长度，m；
df1为相邻裂缝段1中心到交线的平均距离，m；
df2为相邻裂缝段2中心到交线的平均距离，m；
相交裂缝计算式为：
Figure FDA0003989942770000046
Figure FDA0003989942770000047
式中：dSi为面积微元；
Si为裂缝段i面积，m2；xn为面积微元到交线距离，m；
③NNC类型Ⅲ
NNC类型Ⅲ的相交裂缝段间传导率系数为：
Figure FDA0003989942770000048
Figure FDA0003989942770000049
Figure FDA0003989942770000051
式中：
Figure FDA0003989942770000052
为相交裂缝段间传导率系数，mD·m；
Ac为两个相交裂缝段的公共面，m2；
kf1′为相交裂缝段1的渗透率，mD；kf2′为相交裂缝段1的渗透率，mD；
dseg1为相交裂缝段1形心到公共面的距离，m；
dseg2为相交裂缝段2形心到公共面的距离，m；
④裂缝与井眼之间的传导率系数
基于Peaceman数学模型推导出裂缝与井眼之间的传导率系数及裂缝网格的等效井系数分别为：
Figure FDA0003989942770000053
Figure FDA0003989942770000054
式中：WIf为裂缝与井眼之间的传导率系数，mD·m；
re为裂缝网格的等效井系数，m；w为井眼的半径，m；
Δθ为径向井包含在裂缝内的圆心角角度，即Δθ＝2π；
Lf、hf分别为裂缝的长度及高度，m；
考虑注水诱发动态裂缝的问题，引入动态的裂缝参数
Figure FDA0003989942770000055
及kf′：
其中，动态裂缝渗透率kf′为
Figure FDA0003989942770000056
式中：fi为初始裂缝渗透率，mD；
γf裂缝渗透率应力敏感模量，MPa-1；
pfi初始裂缝内压力，MPa；pf裂缝内压力，MPa；
kmi为基质渗透率，mD；
其中，动态裂缝孔隙度
Figure FDA0003989942770000057
为
Figure FDA0003989942770000058
式中：
Figure FDA0003989942770000059
为裂缝孔隙度，无量纲；
Figure FDA00039899427700000510
为初始裂缝孔隙度，无量纲；
Cf为裂缝压缩系数，MPa-1；
pf为裂缝内压力，MPa；
另外，模型的边界条件和初始条件为1)边界条件
①外边界条件
定压力梯度外边界条件：
Figure FDA0003989942770000061
式中：
Figure FDA0003989942770000062
为边界方向上的法向压力梯度，MPa/m；
const为常数；
②内边界条件
定井底压力：
P|rw＝Pwf(x,y,z)＝const(35)式中：P|rw为定井底压力，MPa；
Pwf为井底流压，MPa；const为常数；
2)初始条件
假设数值模拟初始时刻油藏的压力为某一已知函数P0(x,y,z)，油藏数值模拟压力的初始条件为P(x,y,z,0)＝P0(x,y,z)(36)在开发过程中，油藏中油水两相发生流动时，需要确定开发初始时刻的某一相的饱和度场，含水饱和度的初始条件为：
Sw＝Sw0(x,y,z)(37)式中：Sw为油藏各介质***中任意一点含水饱和度的参数，无量纲；
P为油藏各介质***中任意一点压力的参数，MPa。
3.根据权利要求2所述考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：所述模型求解的具体过程为：
在对渗流控制方程进行差分后，对离散模型进行数值求解计算，计算流程如下：
(1)网格划分和初始化：对基质进行网格划分并对基质和裂缝物性进行初始化，并计算传导率；
(2)裂缝网格划分、嵌入：根据裂缝网格坐标，将裂缝嵌入至基质网格中；
(3)计算裂缝宏观渗流参数：根据裂缝不同状态分别计算裂缝的渗透率、孔隙度、相对渗透率等参数；
(4)计算窜流量：分别计算裂缝-基质，裂缝-裂缝的传导率和窜流量，构建线性方程组系数矩阵；
(5)求解渗流微分方程：将计算得到的窜流量带入渗流微分方程中，全隐式进行求解计算；
(6)输出结果：计算至最终时刻，输出压力和饱和度场。
CN202211578528.8A 2022-12-09 2022-12-09 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法 Pending CN115935843A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211578528.8A CN115935843A (zh) 2022-12-09 2022-12-09 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法
CN202310366098.1A CN116579201A (zh) 2022-12-09 2023-04-07 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211578528.8A CN115935843A (zh) 2022-12-09 2022-12-09 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法

Publications (1)

Publication Number Publication Date
CN115935843A true CN115935843A (zh) 2023-04-07

Family

ID=86648660

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202211578528.8A Pending CN115935843A (zh) 2022-12-09 2022-12-09 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法
CN202310366098.1A Pending CN116579201A (zh) 2022-12-09 2023-04-07 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202310366098.1A Pending CN116579201A (zh) 2022-12-09 2023-04-07 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法

Country Status (1)

Country Link
CN (2) CN115935843A (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117077577B (zh) * 2023-10-17 2024-02-02 中国石油大学(华东) 一种适用于低渗透裂缝性油藏的快速模拟及优化方法

Also Published As

Publication number Publication date
CN116579201A (zh) 2023-08-11

Similar Documents

Publication Publication Date Title
Jing et al. Mechanism simulation of oil displacement by imbibition in fractured reservoirs
CN106194145B (zh) 一种多级暂堵深度网络酸压方法
CN107313759B (zh) 低渗稠油油藏直井热采压力分布预测方法及***
Chong et al. Numerical investigation of hydraulic fracturing in transversely isotropic shale reservoirs based on the discrete element method
CN114048695B (zh) 一种基于返排数据的页岩气有效缝网体积反演方法
Furui et al. A Comprehensive Model of High-Rate Matrix-Acid Stimulation for Long Horizontal Wells in Carbonate Reservoirs: Part II—Wellbore/Reservoir Coupled-Flow Modeling and Field Application
CN112541287B (zh) 疏松砂岩压裂充填防砂增产调剖一体化设计方法
CN116306385B (zh) 一种油藏压裂渗吸增能数值模拟方法、***、设备及介质
US20230046288A1 (en) New foamed diverter/sand control model for fluid diversion in integrated wellbore-reservoir system
CN111577269B (zh) 一种基于离散元流固耦合的多簇压裂裂缝形态预测方法
CN111219176A (zh) 一种注水井扩容储层改造方法
CN112392472A (zh) 确定页岩与邻近油层一体化开发方式的方法及装置
CN115935843A (zh) 一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法
Zheng et al. Coupling a geomechanical reservoir and fracturing simulator with a wellbore model for horizontal injection wells
WO2012109191A1 (en) A quantitative method of determining safe steam injection pressure for enhanced oil recovery operations
CN114592840B (zh) 暂堵压裂方法及其应用
CN116629165B (zh) 储层压裂改造区和未改造区参数反演方法、***及设备
CN115081352B (zh) 深水高温高压疏松砂岩端部脱砂工艺设计优化方法及装置
CN113738351B (zh) 一种断缝体油藏物理模型的制作方法和实验方法
CN110566174A (zh) 一种水锁伤害物理模拟方法及装置
CN104392131A (zh) 一种水驱沙过程中破碎岩石渗流场计算方法
Lu et al. A polymer injectivity model establishment and application for early polymer injection
CN115841083A (zh) 确定注水井压驱配注量的方法
Liu et al. Effect of capillary force on performance of shale rock fracturing
CN106529199B (zh) 一种砾岩油藏化学驱井距的确定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20230407