[发明专利]一种高效快速模拟水力压裂的方法

权利要求书

1.一种高效快速模拟水力压裂的方法，该方法包括如下模拟步骤：

步骤1)、根据有限元变分原理，得出水力压裂方程组的弱形式，包括固体应力平衡方程的弱形式和流体压力方程的弱形式；建模过程具体如下：

首先，对于固体应力平衡方程部分：

σ_e＝Dε (2)

其中，上述式(1)-式(3)为固体方程及相关关系式，σ为地层岩石的总应力张量，单位为MPa；ε为地层岩石的应变张量，其为百分数或小数；σ_e为岩石骨架的有效应力张量，单位为MPa；u为地层岩石的位移向量，单位为m；D为岩石骨架的线弹性本构张量，单位为MPa；T表示矩阵的转置符号；

固体方程边界条件定义如下：

其中σ_H、σ_h分别为远场最大、最小水平主应力，单位为MPa；p_f(s,t)为作用在裂缝面上的流体压力大小，单位为MPa；n_t、n_f为远场外边界、裂缝内边界的法向量；s为裂缝位置，单位为m；t为时间，单位为s；

根据有限元变分原理，以上固体方程对应的弱形式为：

其次，对于流体方程部分：

水力裂缝内流压方程如下：

其中w为裂缝宽度，单位为m；Q₀为流入液体的排量，单位为m³/s；δ(s)为Dirac函数；μ为液体粘度，单位为Pa·s；

流体方程边界条件定义如下：

其中，2L为裂缝长度，单位为m；-L和L分别对应裂缝的两个尖端部；

根据有限元变分原理，以上流体方程对应的弱形式为：

综上，式(5)和式(8)分别是得到的水力压裂问题的弱形式；

步骤2)、根据步骤1)得到的固体应力平衡方程的弱形式与流体压力方程的弱形式，根据伽辽金(Galerkin)有限元离散化方法，并对上一步中得到的流体压力方程的弱形式中关于时间的一阶导数采用有限差分近似，得出有限元离散格式；具体如下：

其中，式(9)中的矩阵定义为：

其中式(10)中的矩阵定义为：M^f＝(M^s)^T、N^u、N^p为有限元形函数，为每个有限元网格节点值组成的向量；

步骤3)、对流体方程的刚度矩阵K_f进行奇异值分解(SVD)，分解时间域和空间域；由于流体方程中的立方项具有非线性，对矩阵进行奇异值分解(SVD)，分解时间域和空间域，得到如下：

其中

步骤4)、将PGD(Proper Generalized Decomposition，最佳广义分解法)解代入到有限元离散格式中，分别得到时间域和空间域分解的离散方程；具体如下：

设所求压力场的PGD解具有如下形式：

由式(9)可知，如果求解压力场p，很容易得到位移场u的PGD解，因此只需对流体压力方程部分进行PGD求解即可；

将上式(12)所求的PGD解代入到步骤2)的式(10)中分别得到：

首先固定T，得到关于X的线性代数方程：

其中式中各系数定义如下：

再固定式(13)刚求出的X，得到关于T的一阶常微分方程：

其中式中各系数定义如下：

步骤5)、用交替方向迭代法进行求解即首先固定T，得到关于X的线性代数方程，再利用刚求出的X，得到关于T的一阶常微分方程，当满足一定收敛准则时，即可得到X和T的大小，迭代停止；具体如下：

假定初始项X₀、T₀，应用交替方向迭代法对方程(13)和方程(15)进行交替迭代，当满足下述收敛准则时，即可得到X和T的大小，迭代停止准则如下：

其中ε为允许误差；

步骤6)、重复上述步骤4)和步骤5)，得到上式(17)中每一项X_m和T_m，当满足以下迭代停止准则时，退出循环，否则m＝m+1，直到满足条件(18)时结束；

||X_mT_m-X_m+1T_m+1||＜ε (18)

步骤7)、将步骤6)得到的每一项X和T相乘并求和，得到当前压力场的PGD解；

步骤8)、对流体方程进行不动点迭代，迭代格式如下：

当满足以下迭代准则时，停止迭代：

||p^δ+1-p^δ||＜ε (20)

步骤9)、输出最终压力场p，再根据式(9)的固体离散方程得到的位移场u；

步骤10)、根据得到的位移场u，计算裂缝的几何参数。

2.如权利要求1所述的高效快速模拟水力压裂的方法，其特征在于：所述裂缝几何参数包括缝宽和缝长；其中，缝宽为裂缝面两侧位移的差值，缝长为两个裂缝尖端点间的距离。