Analytical and numerical investigation of multicomponent multiphase WAG displacements

In this paper the ability of analytical solutions for four-component three-phase flow to predict displacement efficiency in water alternating gas (WAG) injection processes is studied. First analytical solutions for Riemann problems with injection compositions that are the average water and gas mixture for various WAG injection schemes are presented. These solutions are compared to numerical calculations with variable slug sizes and used to explore the effect of slug size, injecting water vs gas first, and the average injection composition on displacement efficiency in compositional WAG schemes. The example model is partially miscible WAG injection of water and CO₂ into an oil reservoir containing C₁₀ and CH₄ with and without a mobile aqueous phase present initially. The trailing end of the water and gas profiles are sensitive to whether water or gas is injected first, but the magnitude of the oil bank and the breakthrough time of the injected fluids are accurately predicted by the analytical solutions, even for displacements where large water and gas slugs are injected. Fluctuations in the saturation and composition profiles resulting from the alternating injection sequence in the WAG simulations appear as super-imposed on top of the sequence of rarefaction and shock waves predicted by analytical solutions. As the number of slugs increases, the effect of alternating boundary conditions diminishes and the displacements predicted by numerical calculations converge to the analytical solutions.     

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mathematical Geosciences - Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the...     

Multiple Well Systems with Non‐Darcy Flow

Optimization of groundwater and other subsurface resources requires analysis of multiple‐well systems. The usual modeling approach is to apply a linear flow equation (e.g., Darcy's law in confined aquifers). In such conditions, the composite response of a system of wells can be determined by summating responses of the individual wells (the principle of superposition). However, if the flow velocity increases, the nonlinear losses become important in the near‐well region and the principle of superposition is no longer valid. This article presents an alternative method for applying analytical solutions of non‐Darcy flow for a single‐ to multiple‐well systems. The method focuses on the response of the central injection well located in an array of equally spaced wells, as it is the well that exhibits the highest pressure change within the system. This critical well can be represented as a single well situated in the center of a closed square domain, the width of which is equal to the well spacing. It is hypothesized that a single well situated in a circular region of the equivalent plan area adequately represents such a system. A test case is presented and compared with a finite‐difference solution for the original problem, assuming that the flow is governed by the nonlinear Forchheimer equation.     

Insight from analytical solutions for improved simulation of miscible WAG flooding in one dimension

In this work, the analytical and numerical solutions for modeling miscible gas and water injection into an oil reservoir are presented. Conservation laws with three levels of complexity are considered. Only the most complex model has the correct phase behavior for the example system, which is a multicontact miscible condensing gas drive with simultaneous water and gas injection. Example displacements in which one or both of the simpler models result in accurate simulations in a fraction of the computation time are presented, along with an example in which neither simplified thermodynamic model achieves a truly satisfactory result. A methodology is presented that can be used to establish the accuracy of simplified models in 1-D simulation based on convergence to analytical solutions for the full three-phase system.