An exploration of the IGA method for efficient reservoir simulation

Novel numerical methods present exciting opportunities to improve the efficiency of reservoir simulators. Because potentially significant gains to computational speed and accuracy may be obtained, it is worthwhile explore alternative computational algorithms for both general and case-by-case application to the discretization of the equations of porous media flow, fluid-structure interaction, and/or production. In the present work, the fairly new concept of isogeometric analysis (IGA) is evaluated for its suitability to reservoir simulation via direct comparison with the industry standard finite difference (FD) method and 1st order standard finite element method (SFEM). To this end, two main studies are carried out to observe IGA’s performance with regards to geometrical modeling and ability to capture steep saturation fronts. The first study explores IGA’s ability to model complex reservoir geometries, observing L2 error convergence rates under a variety of refinement schemes. The numerical experimental setup includes an 'S' shaped line sink of varying curvature from which water is produced in a 2D homogenous domain. The accompanying study simplifies the domain to 1D, but adds in multiphase physics that traditionally introduce difficulties associated with modeling of a moving saturation front. Results overall demonstrate promise for the IGA method to be a particularly effective tool in handling geometrically difficult features while also managing typically challenging numerical phenomena.Petroleum and Geosystems Engineerin

Large scale cavity dissolution: From the physical problem to its numerical solution

Dissolution of underground cavities by ground water (or solutions) may cause environmental problems and geological hazards. Efficient modeling and numerical solving of such phenomena are critical for risk analysis. To solve the cavity dissolution problems, we propose to use a porous medium based local non-equilibrium diffuse interface method (DIM) which does not need to track the dissolution fronts explicitly as the sharp front methods (such as ALE). To reduce the grid blocks when using the DIM method, an adaptive mesh refinement (AMR) method is used to have higher resolutions following the moving fronts. An efficient fully implicit scheme is used by taking care of the velocities across the gridblock interfaces on the AMR grid. Numerical examples of salt dissolution under different flow conditions were performed to validate the modeling and numerical solving. Core-scale and reservoir-scale cases were carried out to study the mass transport and the evolution of the profiles of the dissolution fronts. Gravity-driven physical instabilities are found to be more strong in the infinite channel with upper and lower planes than in the 3D tube configuration under the same condition. The implementations with the AMR method also showed a very good computational efficiency, while obtaining good agreement with the finest-grid solutions