[发明专利]一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法

权利要求书

1.一种考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：步骤如下：

步骤1：建立模型；

在油水两相渗流的基本方程基础上，考虑渗吸驱替对油水渗流的双重作用及注水诱导动态裂缝，分别将驱替和渗吸过程的毛管力模型、相对渗透率模型及基质裂缝窜流量带入到渗流控制方程，建立渗吸驱替双重作用下注水诱导动态裂缝油水两相渗流数学模型；

步骤2：模型求解，采用全隐式有限差分方法获得考虑渗吸驱替双重作用的动态裂缝油水两相渗流场数据。

2.根据权利要求1所述考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：所述建立模型的具体过程如下：

(1)建立油水两相渗流的基本微分方程；

在考虑产量项的情况下，油水两相渗流的基本微分方程组为

(2)建立带压渗吸双重毛管力模型

考虑驱替过程的毛管力计算模型为

式中：P_cq为考虑驱替过程的毛管压力，MPa；

P_d为驱替入口压力，MPa；

S_w为含水饱和度，无量纲；S_wr为水相对渗透率趋近于0时的含水饱和度，无量纲；λ为孔隙大小分布指数，无量纲；

考虑渗吸过程的毛管力计算模型为

式中：P_cs为考虑渗吸过程的毛管压力，MPa；

P_neg为渗吸曲线上，水饱和度达到最大值时的最大毛管压力，MPa；

m、n、α为模型拟合参数；

S_iw为束缚水饱和度，无量纲；S_or为残余油饱和度，无量纲；

考虑驱替过程的相渗模型为

式中：K_rwq为考虑驱替过程的水相相对渗透率，无量纲；

_roq为考虑驱替过程的油相相对渗透率，无量纲；

考虑渗吸过程的相渗模型为：

式中：_ros为考虑渗吸过程的水相相对渗透率，无量纲；

K_rws为考虑渗吸过程的油相相对渗透率，无量纲；

在渗流基本模型基础上，考虑渗吸驱替对油水渗流的双重作用，分别将驱替和渗吸过程的毛管力模型和相对渗透率模型，代入到渗流控制方程即可；

(3)建立油水井人工压裂裂缝数值刻画模型EDFM

嵌入式离散裂缝数值模型的控制方程为：

式中：为基质孔隙度，无量纲；为裂缝孔隙度，无量纲；

S_l,m为基质中l相饱和度，无量纲；_l,f为裂缝中l相饱和度，无量纲；

_m为基质渗透率，mD；_f为裂缝渗透率，mD；

k_rl,为基质中l相的相对渗透率，无量纲；k_rl,为裂缝中l相的相对渗透率，无量纲；

l为油相或者水相，o、w；μ_l为l相的粘度，mpa.s；ρ_lg为l相的重度，MPa/m；

为深度，m；

p_l,m为基质中l相的目前地层压力，MPa；p_l,f为裂缝中l相的目前地层压力，MPa；

q_m为油藏基质与嵌入裂缝间的窜流量的窜流量，m³/d；

q_f为嵌入裂缝与油藏基质间的窜流量的窜流量，m³/d；

q_m＝-q_f；

其中，q_f的计算过程包括四种情况；

四种情况分别为三个类型非相邻连接NNC及裂缝与井眼之间的传导率系数；

三个类型非相邻连接NNC定义如下所示：

①NNC类型Ⅰ：裂缝段与其穿过的基质网格间的连接；

②NNC类型Ⅱ：某一裂缝内相邻裂缝段间的连接；

③NNC类型Ⅲ：相交的不同裂缝段之间的连接；

对某一类NNC所连接的两个网格间的裂缝窜流量表示为：

q_f＝λ_lT_NNCAP     (15

式中：

λ_l为l相流度，mD·mPa^-1·s^-1；

T_NNC为NNC对传导率系数，mD·m；

ΔP为相邻网格压力差值，MPa；

式中：k_NNC为NNC对渗透率，mD；

A_ANC为NNC对接触面积，m²；

d_NNC为与NNC连接对特征距离，m；

①NNC类型Ⅰ

NNC类型Ⅰ的基质与嵌入网格裂缝传导率系数为：

式中：T_mf为基质与嵌入网格裂缝传导率系数，mD·m；

A_mf为基质与嵌入网格裂缝的接触面积，m²；

k_mf为基质渗透率和裂缝渗透率的调和平均，mD；

d_mf为基质到裂缝平面的平均法向距离，m；

基质到裂缝平面的平均法向距离d_mf为：

式中：dv为网格块内体积微元，m³；

x_n为该体积微元到裂缝法向距离，m；

V为网格块体积，m³；

②NNC类型Ⅱ

NNC类型Ⅱ相邻裂缝段间传导率系数为：

式中：为相邻裂缝段间传导率系数，mD·m；

k_f1为相邻裂缝段1的渗透率，mD；w_f1为相邻裂缝段1的开度，m；k_f2为相邻裂缝段2的渗透率，mD；w_f2为相邻裂缝段2的开度，m；

L_int为裂缝段交线的长度，m；

d_f1为相邻裂缝段1中心到交线的平均距离，m；

d_f2为相邻裂缝段2中心到交线的平均距离，m；

相交裂缝计算式为：

式中：dS_i为面积微元；

S_i为裂缝段i面积，m²；x_n为面积微元到交线距离，m；

③NNC类型Ⅲ

NNC类型Ⅲ的相交裂缝段间传导率系数为：

式中：为相交裂缝段间传导率系数，mD·m；

A_c为两个相交裂缝段的公共面，m²；

k_f1′为相交裂缝段1的渗透率，mD；k_f2′为相交裂缝段1的渗透率，mD；

d_seg1为相交裂缝段1形心到公共面的距离，m；

d_seg2为相交裂缝段2形心到公共面的距离，m；

④裂缝与井眼之间的传导率系数

基于Peaceman数学模型推导出裂缝与井眼之间的传导率系数及裂缝网格的等效井系数分别为：

式中：WI_f为裂缝与井眼之间的传导率系数，mD·m；

r_e为裂缝网格的等效井系数，m；_w为井眼的半径，m；

Δθ为径向井包含在裂缝内的圆心角角度，即Δθ＝2π；

L_f、h_f分别为裂缝的长度及高度，m；

考虑注水诱发动态裂缝的问题，引入动态的裂缝参数及k_f′：

其中，动态裂缝渗透率k_f′为

式中：_fi为初始裂缝渗透率，mD；

γ_f裂缝渗透率应力敏感模量，MPa^-1；

p_fi初始裂缝内压力，MPa；p_f裂缝内压力，MPa；

k_mi为基质渗透率，mD；

其中，动态裂缝孔隙度为

式中：为裂缝孔隙度，无量纲；为初始裂缝孔隙度，无量纲；

C_f为裂缝压缩系数，MPa^-1；

p_f为裂缝内压力，MPa；

另外，模型的边界条件和初始条件为

1)边界条件

①外边界条件

定压力梯度外边界条件：

式中：为边界方向上的法向压力梯度，MPa/m；

const为常数；

②内边界条件

定井底压力：

P|_rw＝P_wf(x，y，z)＝const    (35

式中：P|_rw为定井底压力，MPa；

P_wf为井底流压，MPa；const为常数；

2)初始条件

假设数值模拟初始时刻油藏的压力为某一已知函数P₀(x，y，z)，油藏数值模拟压力的初始条件为

P(x，y，z，0)＝P₀(x，y，z)    (36)

在开发过程中，油藏中油水两相发生流动时，需要确定开发初始时刻的某一相的饱和度场，含水饱和度的初始条件为：

S_w＝S_w0(x，y，z)      (37)

式中：S_w为油藏各介质系统中任意一点含水饱和度的参数，无量纲；

P为油藏各介质系统中任意一点压力的参数，MPa。

3.根据权利要求2所述考虑渗吸机理的注水诱导动态裂缝渗流数值模拟方法，其特征在于：所述模型求解的具体过程为：

在对渗流控制方程进行差分后，对离散模型进行数值求解计算，计算流程如下：

(1)网格划分和初始化：对基质进行网格划分并对基质和裂缝物性进行初始化，并计算传导率；

(2)裂缝网格划分、嵌入：根据裂缝网格坐标，将裂缝嵌入至基质网格中；

(3)计算裂缝宏观渗流参数：根据裂缝不同状态分别计算裂缝的渗透率、孔隙度、相对渗透率等参数；

(4)计算窜流量：分别计算裂缝-基质，裂缝-裂缝的传导率和窜流量，构建线性方程组系数矩阵；

(5)求解渗流微分方程：将计算得到的窜流量带入渗流微分方程中，全隐式进行求解计算；

(6)输出结果：计算至最终时刻，输出压力和饱和度场。