|
| A mimetic finite difference based novel streamline simulation method for two-phase porous flow |
| Received:June 11, 2023 Revised:August 27, 2023 |
| View Full Text View/Add Comment Download reader |
| DOI：10.7511/jslx20230611001 |
| KeyWord:mimetic finite difference method streamline-based simulation discontinuous Galerkin method oil-water two-phase porous flow |
| Author | Institution |
| 饶翔 |
油气钻采工程湖北省重点实验室长江大学, 武汉 , 中国;阿卜杜拉国王科技大学 物理科学与工程部Ali I. Al-Naimi石油工程研究中心, 吉达 23955-6900, 沙特阿拉伯;沙特阿美国家石油公司 EXPEC高级研究中心, 达曼 3, 沙特阿拉伯 |
| 何旭鹏 |
阿卜杜拉国王科技大学 物理科学与工程部Ali I. Al-Naimi石油工程研究中心, 吉达 23955-6900, 沙特阿拉伯;沙特阿美国家石油公司 EXPEC高级研究中心, 达曼 3, 沙特阿拉伯 |
| 赵辉 |
油气钻采工程湖北省重点实验室长江大学, 武汉 , 中国 |
| 刘怡娜 |
油气钻采工程湖北省重点实验室长江大学, 武汉 , 中国 |
|
| Hits: 16 |
| Download times: 12 |
| Abstract: |
| In this paper,a new streamline simulation method based on mimetic finite difference method (denoted as MFD-SL) is proposed for the first time.This method uses mimetic finite difference method applicable for full tensor permeability and unstructured meshes to obtain high-accuracy pressure and flow velocity distributions and a simple streamline tracking method to quickly obtain the streamline distribution in a two-dimensional triangular or a three-dimensional tetrahedral mesh and calculate the flight time distribution on the streamlines.The discontinuous Galerkin method based on the simplified Runge-Kutta scheme is used to obtain the water saturation distribution of low numerical dissipation on each one-dimensional streamline.Theoretically,compared with traditional streamline simulation methods,this method can be applied to both complex-shaped reservoirs and anisotropic reservoirs.Compared with finite volume and mimetic finite difference methods,which only calculates the grid-average saturation,this method can significantly reduce the numerical dissipation of saturation calculation.Compared with high-order methods such as the discontinuous Galerkin method,which calculates the vertex saturation of the grid,this method can avoid the time step restriction of the strong stability condition,improve the computational efficiency,and further reduce the numerical dissipation.Numerical examples in this paper also show that the proposed method has the advantages of both accuracy and efficiency in the case of anisotropic permeability,unstructured meshing and complex geometry,and can provide intuitive flow diagnosis results. |
|
|
|