Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
This is the second paper of a series where we introduce a control volume based finite element method （CVFEM） to simulate multiphase flow in porous media. This is a fully conservative method able to deal with unstruc...This is the second paper of a series where we introduce a control volume based finite element method （CVFEM） to simulate multiphase flow in porous media. This is a fully conservative method able to deal with unstructured grids which can be used for representing any complexity of reservoir geometry and its geological objects in an accurate and efficient manner. In order to deal with the inherent heterogeneity of the reservoirs, all operations related to discretization are performed at the element level in a manner similar to classical finite element method （FEM）. Moreover, the proposed method can effectively reduce the so-called grid orientation effects. In the first paper of this series, we presented this method and its application for incompressible and immiscible two-phase flow simulation in homogeneous and heterogeneous porous media. In this paper, we evaluate the capability of the method in the solution of highly nonlinear and coupled partial differential equations by simulating hydrocarbon reservoirs using the black-oil model. Furthermore, the effect of grid orientation is investigated by simulating a benchmark waterflooding problem. The numerical results show that the formulation presented here is efficient and accurate for solving the bubble point and three-phase coning problems.展开更多
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金Iranian Offshore OilCompany (IOOC) for financial support of this work
文摘This is the second paper of a series where we introduce a control volume based finite element method （CVFEM） to simulate multiphase flow in porous media. This is a fully conservative method able to deal with unstructured grids which can be used for representing any complexity of reservoir geometry and its geological objects in an accurate and efficient manner. In order to deal with the inherent heterogeneity of the reservoirs, all operations related to discretization are performed at the element level in a manner similar to classical finite element method （FEM）. Moreover, the proposed method can effectively reduce the so-called grid orientation effects. In the first paper of this series, we presented this method and its application for incompressible and immiscible two-phase flow simulation in homogeneous and heterogeneous porous media. In this paper, we evaluate the capability of the method in the solution of highly nonlinear and coupled partial differential equations by simulating hydrocarbon reservoirs using the black-oil model. Furthermore, the effect of grid orientation is investigated by simulating a benchmark waterflooding problem. The numerical results show that the formulation presented here is efficient and accurate for solving the bubble point and three-phase coning problems.
