**Peer Reviewed**

**Source**

**SPE Disciplines**

**Journal**

**Publisher**

**Theme**

**Concept Tag**

- basis function (1)
- boundary condition (1)
- coarse block (1)
- coarse grid (1)
- Computation (1)
- construction (1)
- convergence (1)
- correction function (1)
- equation (2)
- finite-volume restriction operator (1)
- flow in porous media (1)
- Fluid Dynamics (1)
- formulation (1)
- full construction (1)
- grid (1)
- iteration (2)
- iteration matrix (1)
- june 2012 (1)
- march 2012 (1)
- mass conservation (1)
- matrix (1)
- Modeling & Simulation (2)
- msfvm (1)
- Multiscale Method (1)
- operator (2)
- permeability (1)
- permeability field (1)
- preconditioner (1)
- prolongation operator (2)
- reservoir simulation (2)
- restriction operator (1)
- saturation (1)
- saturation equation (1)
- scaling method (2)
- TAM (1)
- Tchelepi (1)
- Timestep (1)
- transport equation (1)
- Upstream Oil & Gas (2)
- velocity field (1)

**Country**

**Industry**

**Oilfield Places**

**Technology**

**File Type**

Goto

Layer | Fill | Outline |
---|

Theme | Visible | Selectable | Appearance | Zoom Range (now: 0) |
---|

Fill | Stroke |
---|---|

Zhou, Hui (ConocoPhillips Subsurface Technology) | Tchelepi, Hamdi A. (Stanford University)

An efficient two-stage algebraic multiscale solver (TAMS) that converges to the fine-scale solution is described. The first (global) stage is a multiscale solution obtained algebraically for the given fine-scale problem. In the second stage, a local preconditioner, such as the Block ILU (BILU) or the Additive Schwarz (AS) method, is used. Spectral analysis shows that the multiscale solution step captures the low-frequency parts of the error spectrum quite well, while the local preconditioner represents the high-frequency components accurately. Combining the two stages in an iterative scheme results in efficient treatment of all the error components associated with the fine-scale problem. TAMS is shown to converge to the reference fine-scale solution. Moreover, the eigenvalues of the TAMS iteration matrix show significant clustering, which is favorable for Krylov-based methods. Accurate solution of the nonlinear saturation equations (i.e., transport problem) requires having locally conservative velocity fields. TAMS guarantees local mass conservation by concluding the iterations with a multiscale finite-volume step. We demonstrate the performance of TAMS using several test cases with strong permeability heterogeneity and large-grid aspect ratios. Different choices in the TAMS algorithm are investigated, including the Galerkin and finite-volume restriction operators, as well as the BILU and AS preconditioners for the second stage. TAMS for the elliptic flow problem is comparable to state-of-the-art algebraic multigrid methods, which are in wide use. Moreover, the computational time of TAMS grows nearly linearly with problem size.

boundary condition, coarse grid, convergence, equation, finite-volume restriction operator, flow in porous media, Fluid Dynamics, grid, iteration, iteration matrix, june 2012, mass conservation, matrix, Modeling & Simulation, msfvm, operator, permeability, permeability field, preconditioner, prolongation operator, reservoir simulation, restriction operator, scaling method, TAM, Upstream Oil & Gas

Country:

- North America > United States (1.00)
- Europe (0.93)

Oilfield Places:

- Europe > Norway > North Sea > Tarbert Formation (0.99)
- Europe > Germany > North Sea > Tarbert Formation (0.99)

SPE Disciplines:

Zhou, Hui (ConocoPhilips) | Lee, Seong H. (Chevron Energy Technology Company) | Tchelepi, Hamdi A. (Stanford University)

Recent advances in multiscale methods have shown great promise in modeling multiphase flow in highly detailed heterogeneous domains. Existing multiscale methods, however, solve for the flow field (pressure and total velocity) only. Once the fine-scale flow field is reconstructed, the saturation equations are solved on the fine scale. With the efficiency in dealing with the flow equations greatly improved by multiscale formulations, solving the saturation equations on the fine scale becomes the relatively more expensive part. In this paper, we describe an adaptive multiscale finite-volume (MSFV) formulation for nonlinear transport (saturation) equations. A general algebraic multiscale formulation consistent with the operator-based framework proposed by Zhou and Tchelepi (*SPE Journal*, June 2008, pages 267-273) is presented. Thus, the flow and transport equations are solved in a unified multiscale framework. Two types of multiscale operators--namely, restriction and prolongation--are used to construct the multiscale saturation solution. The restriction operator is defined as the sum of the fine-scale transport equations in a coarse gridblock. Three adaptive prolongation operators are defined according to the local saturation history at a particular coarse block. The three operators have different computational complexities, and they are used adaptively in the course of a simulation run. When properly used, they yield excellent computational efficiency while preserving accuracy. This adaptive multiscale formulation has been tested using several challenging problems with strong heterogeneity, large buoyancy effects, and changes in the well operating conditions (e.g., switching injectors and producers during simulation). The results demonstrate that adaptive multiscale transport calculations are in excellent agreement with fine-scale reference solutions, but at a much lower computational cost.

basis function, coarse block, Computation, construction, correction function, equation, formulation, full construction, iteration, march 2012, Modeling & Simulation, Multiscale Method, operator, prolongation operator, reservoir simulation, saturation, saturation equation, scaling method, Tchelepi, Timestep, transport equation, Upstream Oil & Gas, velocity field

SPE Disciplines:

Thank you!