06.08.2013 Views

Sequential Methods for Coupled Geomechanics and Multiphase Flow

Sequential Methods for Coupled Geomechanics and Multiphase Flow

Sequential Methods for Coupled Geomechanics and Multiphase Flow

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

1.5. FINITE-VOLUME AND FINITE-ELEMENT METHODS 9<br />

algebraic constraint (Brenan et al., 1996). They claimed that the differential algebraic equa-<br />

tions have first-order accuracy when the first step in the algebraic equation is redundant.<br />

However, this is not applicable to the coupled flow <strong>and</strong> quasi-static mechanics because the<br />

mechanical problem is not redundant <strong>for</strong> the first time step. For example, consolidation<br />

problems are often driven by instant loading in the mechanical problem during the first time<br />

step. Rather, the work of Vijalapura et al. (2005) supports the fact that typical sequential<br />

methods can provide zeroth order accuracy in time, when the first step in the algebraic<br />

equation is not redundant.<br />

1.5 Finite-Volume <strong>and</strong> Finite-Element <strong>Methods</strong><br />

Numerical experiments are per<strong>for</strong>med using the finite-volume <strong>and</strong> finite-element methods<br />

(the FVM <strong>and</strong> FEM) <strong>for</strong> flow <strong>and</strong> mechanics, respectively. The benefits from this discretiza-<br />

tion strategy are as follows.<br />

1. The choice <strong>for</strong> the discretization <strong>and</strong> the primary unknowns are practical <strong>and</strong> natural<br />

in order to make use of a reservoir simulator <strong>and</strong> a geomechanics simulator in a<br />

sequential fashion.<br />

2. Reservoir simulation with the finite-volume method yields local mass conservation,<br />

while that with the finite-element method is not locally conservative.<br />

3. This mixed-space discretization strategy avoids spatial oscillations at early time, which<br />

are typically encountered in consolidation problems.<br />

There are two reasons <strong>for</strong> the early time spatial oscillations in the nodal based finite-<br />

element methods. One is the violation of the LBB condition (Fortin <strong>and</strong> Brezzi, 1991) <strong>and</strong><br />

the other is the discontinuity of pressure at the drainage boundary. At early time, the<br />

coupled problem converges to the undrained mechanical response. The LBB condition is<br />

required to solve the incompressible mechanical problem (Fortin <strong>and</strong> Brezzi, 1991), <strong>and</strong> it<br />

must be satisfied when the fluid <strong>and</strong> solid grains are incompressible (Murad <strong>and</strong> Loula,<br />

1992, 1994; Truty <strong>and</strong> Zimmermann, 2006; White <strong>and</strong> Borja, 2008).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!