多孔介质全阶渗透率张量的预测方法

摘要

本发明提供了一种多孔介质全阶渗透率张量的预测方法，该预测方法包括：建立多孔介质内流体流动的控制模型，并将多孔介质划分成多个交错的网格；使流体在多孔介质内沿第一方向流动，利用有限差分法对控制模型进行求解，以确定网格中不同点处的局部速度分量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_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{\mu }{\rho g}\left[\begin{array}{cc}{u_1}^D& {u_2}^D\\ {v_1}^D& {v_2}^D\end{array}\right];$其中，μ为多孔介质内流体的动力黏度；ρ为多孔介质内流体的密度；g为重力加速度。