[发明专利]多孔介质全阶渗透率张量的预测方法

权利要求书

1.一种多孔介质全阶渗透率张量的预测方法，其特征在于，所述预测方法包括：

利用稳态斯托克斯方程建立多孔介质内流体流动的控制模型，并将所述多孔介质划分成多个交错的网格；

使流体在所述多孔介质内沿第一方向流动，利用有限差分法对所述控制模型进行差分求解，以确定所述网格中不同点处的局部速度分量u_1i,j、v_1i,j；

根据所述局部速度分量u_1i,j、v_1i,j，利用分段二次抛物线数值积分法确定多孔介质内流体的第一达西速度u₁^D、v₁^D；

使流体在所述多孔介质内沿与所述第一方向垂直的第二方向流动，利用有限差分法对所述控制模型进行差分求解，以确定所述网格中不同点处的局部速度分量u_2i,j、v_2i,j；

根据所述局部速度分量u_2i,j、v_2i,j，利用分段二次抛物线数值积分法确定多孔介质内流体的第二达西速度u₂^D、v₂^D；

根据所述第一达西速度u₁^D、v₁^D及第二达西速度u₂^D、v₂^D，确定所述多孔介质内流体的全阶渗透率张量K：

$K=\frac{\mathrm{\μ}}{\mathrm{\ρ}g}\left[\begin{array}{cc}{u_1}^D& {u_2}^D\\ {v_1}^D& {v_2}^D\end{array}\right];$

其中，μ为多孔介质内流体的动力黏度；ρ为多孔介质内流体的密度；g为重力加速度。

2.根据权利要求1所述的预测方法，其特征在于，如果所述网格中待求解点位于多孔介质的流体区域且其近邻点也位于流体区域，利用有限差分法对所述控制模型进行差分求解，以确定所述网格中不同点处的局部速度分量u_1i,j、v_1i,j，包括：

利用位于多孔介质流体区域中的所述近邻点处的速度分量对所述控制模型进行差分求解，以确定所述网格中不同点处的局部速度分量u_1i,j、v_1i,j。

3.根据权利要求1所述的预测方法，其特征在于，如果所述网格中待求解点位于多孔介质的流体区域且其近邻点位于固体区域，利用有限差分法对所述控制模型进行差分求解，以确定所述网格中不同点处的局部速度分量u_1i,j、v_1i,j，包括：

将位于多孔介质固体区域的所述近邻点的速度分量替换为距离其最近的流固界面处的速度分量；

沿与所述流固界面处的速度分量相垂直的方向往流体区域移动1/4个网格；

对所述控制模型进行差分求解，以确定所述待求解点处的局部速度分量u_1i,j、v_1i,j。