[发明专利]一种非均质压裂后储层流动数值模拟方法及系统

权利要求书

1.一种非均质压裂后储层流动数值模拟方法，其特征在于，该方法包括：

获取非均质压裂后的离散裂缝网络和改造区储层范围，所述改造区储层范围包括近井区储层和远井区储层；

将所述离散裂缝网络划分为若干裂缝微元，以及分别将所述近井区储层的边界和远井区储层的边界划分为若干边界单元；

分别将所述近井区储层的流动微分方程和远井区储层的流动微分方程变换为边界积分方程，根据所述裂缝微元和所述边界单元，对变换后的边界积分方程进行求解，得到近井区储层的边界元解和远井区储层的边界元解；

将离散裂缝流动方程、近井区储层的边界元解和远井区储层的边界元解进行耦合，得到离散裂缝网络-两区流动耦合矩阵；

根据所述离散裂缝网络-两区流动耦合矩阵，对所述近井区储层和远井区储层进行流动数值模拟。

2.根据权利要求1所述的非均质压裂后储层流动数值模拟方法，其特征在于，在划分的步骤中包括：

将所述离散裂缝网络划分为若干裂缝微元，通过“岔型”变化法消去所述裂缝微元的交汇点；

利用线性边界元法分别将所述近井区储层的边界和远井区储层的边界划分为若干边界单元。

3.根据权利要求1所述的非均质压裂后储层流动数值模拟方法，其特征在于，所述分别将所述近井区储层的流动微分方程和远井区储层的流动微分方程变换为边界积分方程，包括：

分别确定近井区储层的流动微分方程和远井区储层的流动微分方程的基本解；

根据所述基本解，分别将近井区储层的流动微分方程和远井区的流动微分方程变换为边界积分方程。

4.根据权利要求3所述的非均质压裂后储层流动数值模拟方法，其特征在于，所述确定近井区储层的流动微分方程和远井区储层的流动微分方程的基本解包括：

对所述近井区储层的流动微分方程和远井区储层的流动微分方程进行Laplace变换，得到Hemholz方程；对所述Hemholz方程进行求解，得到所述基本解。

5.根据权利要求1所述的非均质压裂后储层流动数值模拟方法，其特征在于，在求解近井区储层的边界元解和远井区储层的边界元解的步骤中包括：

使用线性边界元法对所述边界积分方程进行离散化；

根据所述裂缝微元和所述边界单元，对离散后的边界积分方程进行求解，得到所述近井区储层的边界元解和远井区储层的边界元解。

6.根据权利要求1所述的非均质压裂后储层流动数值模拟方法，其特征在于，在将离散裂缝流动方程、近井区储层的边界元解和远井区储层的边界元解进行耦合的步骤中包括：

根据离散裂缝流动方程，建立离散裂缝网络流动的Laplace空间矩阵方程，求解所述Laplace空间矩阵方程，得到所述Laplace空间矩阵方程的求解矩阵；

基于离散裂缝面上的流量和压力的连续性条件，耦合所述Laplace空间矩阵方程的求解矩阵、所述近井区储层的边界元解、远井区储层的边界元解和生产边界条件，得到所述离散裂缝网络-两区流动耦合矩阵。

7.根据权利要求1-6任一项所述的非均质压裂后储层流动数值模拟方法，其特征在于，所述方法还包括：

根据所述近井区储层和远井区储层流动数值模拟结果，确定不同时刻储层流动的压力分布和储层动用情况；

对所述不同时刻储层流动的压力分布和储层动用情况进行可视化处理。