Seong H. Lee

Multi-scale finite-volume method for elliptic problems in subsurface flow simulation

In this paper we present a multi-scale finite-volume (MSFV) method to solve elliptic problems with many spatial scales arising from flow in porous media. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of the differential operator. This leads to a multi-point discretization scheme for the finite-volume solution algorithm. Transmissibilities for the MSFV have to be constructed only once as a preprocessing step and can be computed locally. Therefore this step is perfectly suited for massively parallel computers. Furthermore, a conservative fine-scale velocity field can be constructed from the coarse-scale pressure solution. Two sets of locally computed basis functions are employed. The first set of basis functions captures the small-scale heterogeneity of the underlying permeability field, and it is computed in order to construct the effective coarse-scale transmissibilities. A second set of basis functions is required to construct a conservative fine-scale velocity field. The accuracy and efficiency of our method is demonstrated by various numerical experiments.

Hierarchical modeling of flow in naturally fractured formations with multiple length scales

This paper describes hierarchical approach to modeling flow in a naturally fractured formation. Our model is based on calculating the effective permeability of a fractured formation, as a function of grid block size, and using the results in a conventional finite difference flow simulator. On the basis of their length (lf) relative to the finite difference grid size (lg), fractures are classified as belonging to one of three groups: (1) short fractures (lf≪ lg), (2) medium-length fractures (lf ∼ lg), and (3) long fractures (lf≫ lg). The effects of the fractures belonging to each class are computed in a hierarchical manner. The permeability contribution from short fractures is derived in an analytical expression and used as an enhanced matrix permeability for the next-scale (medium-length) calculation. The effective matrix permeability associated with medium-length fractures is numerically solved using a boundary element method. The long fractures are modeled explicitly as major fluid conduits. As numerical examples, tracer transport in fractured formations was illustrated. The numerical results clearly indicated that effective tensor permeability well represented directional, enhanced permeability in fractured formations. The fluid-conduit formulation captured the efficient fluid transport by long fractures.

Adaptive Multiscale Finite-Volume Method for Multiphase Flow and Transport in Porous Media

We present a multiscale finite-volume (MSFV) method for multiphase flow and transport in heterogeneous porous media. The approach extends our recently developed MSFV method for single-phase flow. We use a sequential scheme that deals with flow (i.e., pressure and total velocity) and transport (i.e., saturation) separately and differently. For the flow problem, we employ two different sets of basis functions for the reconstruction of a conservative fine-scale total velocity field. Our basis functions are designed to have local support, and that allows for adaptive computation of the flow field. We use a criterion based on the time change of the total mobility field to decide when and where to recompute our basis functions. We show that at a given time step, only a small fraction of the basis functions needs to be recomputed. Numerical experiments of difficult two-dimensional and three-dimensional test cases demonstrate the accuracy, computational efficiency, and overall scalability of the method.

Adaptive fully implicit multi-scale finite-volume method for multi-phase flow and transport in heterogeneous porous media

We describe a sequential fully implicit (SFI) multi-scale finite volume (MSFV) algorithm for nonlinear multi-phase flow and transport in heterogeneous porous media. The method extends the recently developed multiscale approach, which is based on an IMPES (IMplicit Pressure, Explicit Saturation) scheme [P. Jenny, S.H. Lee, H.A. Tchelepi, Adaptive multi-scale finite volume method for multi-phase flow and transport, Multiscale, Model. Simul. 3 (2005) 50-64]. That previous method was tested extensively and with a series of difficult test cases, where it was clearly demonstrated that the multiscale results are in excellent agreement with reference fine-scale solutions and that the computational efficiency of the MSFV algorithm is much higher than that of standard reservoir simulators. However, the level of detail and range of property variability included in reservoir characterization models continues to grow. For such models, the explicit treatment of the transport problem (i.e. saturation equations) in the IMPES-based multiscale method imposes severe restrictions on the time step size, and that can become the major computational bottleneck. Here we show how this problem is resolved with our sequential fully implicit (SFI) MSFV algorithm. Simulations of large (million cells) and highly heterogeneous problems show that the results obtained with the implicit multi-scale method are in excellent agreement with reference fine-scale solutions. Moreover, we demonstrate the robustness of the coupling scheme for nonlinear flow and transport, and we show that the MSFV algorithm offers great gains in computational efficiency compared to standard reservoir simulation methods.

Modeling of subgrid effects in coarse-scale simulations of transport in heterogeneous porous media

A methodology for incorporating subgrid effects in coarse-scale numerical models of flow in heterogeneous porous media is presented. The method proceeds by upscaling the deterministic fine-grid permeability description and then solving the pressure equation over the coarse grid to obtain coarse-scale velocities. The coarse-grid saturation equation is formed through a volume average of the fine-scale equations and includes terms involving both the average and fluctuating components of the velocity field. The terms involving the fluctuating components are subgrid effects that appear as length- and time-dependent dispersivities. A simplified model for the coarse-scale dispersivity, in terms of these subgrid velocity fluctuations, is proposed, and a numerical scheme based on this model is implemented. Results using the new method are presented for a variety of two- dimensional heterogeneous systems characterized by moderate permeability correlation length in the dominant (horizontal) flow direction and small correlation length in the vertical direction. The new method is shown to provide much more accurate results than comparable coarse-grid models that do not contain the subgrid treatment. Extensions of the overall methodology to handle more general systems are discussed.

Well Modeling in the Multiscale Finite Volume Method for Subsurface Flow Simulation

A multiscale method for effective handling of wells (source/sink terms) in the simulation of multiphase flow and transport processes in heterogeneous porous media is developed. The approach extends the multiscale finite volume (MSFV) framework. Our multiscale well model allows for accurate reconstruction of the fine‐scale pressure and velocity fields in the vicinity of wells. Accurate and computationally efficient modeling of complex wells is a prerequisite for field applications, and the ability to model wells within the MSFV framework makes it possible to solve large‐scale heterogeneous problems of practical interest. Our approach consists of removal of the well singularity from the multiscale solution via a local change of variables and the computation of a smoothly varying background field instead. The well effects are computed using a separate basis function, which is superposed on the background solution to yield accurate representation of the flow field. The multiscale well treatment accounts for b...

Adaptive multiscale finite-volume method for nonlinear multiphase transport in heterogeneous formations

In the previous multiscale finite-volume (MSFV) method, an efficient and accurate multiscale approach was proposed to solve the elliptic flow equation. The reconstructed fine-scale velocity field was then used to solve the nonlinear hyperbolic transport equation for the fine-scale saturations using an overlapping Schwarz scheme. A coarse-scale system for the transport equations was not derived because of the hyperbolic character of the governing equations and intricate nonlinear interactions between the saturation field and the underlying heterogeneous permeability distribution. In this paper, we describe a sequential implicit multiscale finite-volume framework for coupled flow and transport with general prolongation and restriction operations for both pressure and saturation, in which three adaptive prolongation operators for the saturation are used. In regions with rapid pressure and saturation changes, the original approach, with full reconstruction of the velocity field and overlapping Schwarz, is used to compute the saturations. In regions where the temporal changes in velocity or saturation can be represented by asymptotic linear approximations, two additional approximate prolongation operators are proposed. The efficiency and accuracy are evaluated for two-phase incompressible flow in two- and three-dimensional domains. The new adaptive algorithm is tested using various models with homogeneous and heterogeneous permeabilities. It is demonstrated that the multiscale results with the adaptive transport calculation are in excellent agreement with the fine-scale solutions. Furthermore, the adaptive multiscale scheme of flow and transport is much more computationally efficient compared with the previous MSFV method and conventional fine-scale reservoir simulation methods.

An iterative multiscale finite volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

An Operator Formulation of the Multiscale Finite-Volume Method with Correction Function

The multiscale finite-volume (MSFV) method has been derived to efficiently solve large problems with spatially varying coefficients. The fine-scale problem is subdivided into local problems that can be solved separately and are coupled by a global problem. This algorithm, in consequence, shares some characteristics with two-level domain decomposition (DD) methods. However, the MSFV algorithm is different in that it incorporates a flux reconstruction step, which delivers a fine-scale mass conservative flux field without the need for iterating. This is achieved by the use of two overlapping coarse grids. The recently introduced correction function allows for a consistent handling of source terms, which makes the MSFV method a flexible algorithm that is applicable to a wide spectrum of problems. It is demonstrated that the MSFV operator, used to compute an approximate pressure solution, can be equivalently constructed by writing the Schur complement with a tangential approximation of a single-cell overlapping grid and incorporation of appropriate coarse-scale mass-balance equations.

A finite-volume method with hexahedral multiblock grids for modeling flow in porous media

This paper presents a finite-volume method for hexahedral multiblock grids to calculate multiphase flow in geologically complex reservoirs. Accommodating complex geologic and geometric features in a reservoir model (e.g., faults) entails non-orthogonal and/or unstructured grids in place of conventional (globally structured) Cartesian grids. To obtain flexibility in gridding as well as efficient flow computation, we use hexahedral multiblock grids. These grids are locally structured, but globally unstructured. One major advantage of these grids over fully unstructured tetrahedral grids is that most numerical methods developed for structured grids can be directly used for dealing with the local problems. We present several challenging examples, generated via a commercially available tool, that demonstrate the capabilities of hexahedral multiblock gridding. Grid quality is discussed in terms of uniformity and orthogonality. The presence of non-orthogonal grid and full permeability tensors requires the use of multi-point discretization methods. A flux-continuous finite-difference (FCFD) scheme, previously developed for stratigraphic hexahedral grid with full-tensor permeability, is employed for numerical flow computation. We extend the FCFD scheme to handle exceptional configurations (i.e. three- or five-cell connections as opposed to the regular four), which result from employing multiblock gridding of certain complex objects. In order to perform flow simulation efficiently, we employ a two-level preconditioner for solving the linear equations that results from the wide stencil of the FCFD scheme. The individual block, composed of cells that form a structured grid, serves as the local level; the higher level operates on the global block configuration (i.e. unstructured component). The implementation uses an efficient data structure where each block is wrapped with a layer of neighboring cells. We also examine splitting techniques [14] for the linear systems associated with the wide stencils of our FCFD operator. We present three numerical examples that demonstrate the method: (1) a pinchout, (2) a faulted reservoir model with internal surfaces and (3) a real reservoir model with multiple faults and internal surfaces.