[发明专利]基于迎风GFDM的多孔介质油水两相流计算方法

权利要求书

1.一种基于迎风GFDM的多孔介质油水两相流计算方法，其特征在于，包括以下步骤：

S1、建立多孔介质油水两相流模型；

所述多孔介质油水两相流模型包括油相物质守恒方程、水相物质守恒方程和辅助方程，其中多孔介质中油水两相流互不相溶；

其中，所述油相物质守恒方程为：

上式中：k是渗透率，mD；k_ro＝k_ro(S_w)是与含水饱和度S_w有关的油相相对渗透率；μ_o是油相粘度，mPa·s；p_o是油相压力，MPa；q_o是油相源汇项，1/天；t是时间，天；S_o是油相饱和度；φ＝φ(p)是与压力有关的储层孔隙度；

所述水相物质守恒方程为：

上式中：k_rw＝k_rw(S_w)是与含水饱和度有关的水相相对渗透率；μ_w是水相粘度，mPa·s；p_w是水相压力，MPa；q_w是水相源汇项，1/天；

所述辅助方程为：

S_o＝1-S_w；

p_w＝p_o-P_c(S_w)；

k_rw＝k_rw(S_w)；

k_ro＝k_ro(S_w)

                                (3)

上式中：p_o为油相压力；S_w为含水饱和度；C_t＝C_t(S_w)是综合压缩系数，1/MPa；p_o⁰是储层原始油相压力；φ₀是p_o⁰对应的孔隙度；

S2、基于迎风GFDM对所述多孔介质油水两相流模型进行离散，得到所述多孔介质油水两相流模型的全隐式离散格式；

S3、利用基于牛顿迭代和自动微分的非线性求解器对所述多孔介质油水两相流模型的全隐式离散格式进行求解，实现压力和含水饱和度的高精度无网格计算。