Field-scale application of a semi-analytical model for estimation of CO2 and brine leakage along old wells

Celia, Michael A.; Nordbotten, Jan Martin; Court, Benjamin; Dobossy, M.; Bachu, Stefan

doi:10.1016/j.ijggc.2010.10.005

Cited by 138 publications

(123 citation statements)

References 34 publications

Supporting

Mentioning

118

Contrasting

Order By: Relevance

“…Therefore, the first part of this section provides a detailed description of work presented in [40] and [7].…”

Section: Methodsmentioning

confidence: 99%

“…Reference [7] expresses this radial overpressure response, p, at the bottom of a confined aquifer for a single well injecting CO 2 as: (1) where p 0 and p are the initial and resulting fluid pressures at the bottom of the aquifer,  is fluid density, g is gravitational acceleration, H is aquifer thickness, and subscripts b and c denote phase types brine and CO 2 , respectively. In addition, p' is a dimensionless function defined as:…”

Section: The Estimating Leakage Semi-analytically (Elsa) Algorithmmentioning

confidence: 99%

“…diffuse leakage [10]), this work focuses upon the multi-phase subsurface flow model proposed by [40] and further developed by [7] because it is the only semi-analytical model able to simulate multi-phase flow in domains having multiple injection wells and multiple aquifer and aquitard (i.e. caprock) layers.…”

Section: Introductionmentioning

confidence: 99%

“…This algorithm creates a set of linear equations describing the pressure distribution throughout the domain by superimposing pressure changes caused by each source or sink in each aquifer. The general algorithm presented in [38] in conjunction with the development of a multiphase pressure response function [33,34,35,38] has led to a semi-analytical CO 2 leakage algorithm, presented in [40] and expounded upon in [7] and [36], which estimates both brine and CO 2 flux across confining layers resulting from the injection of CO 2 . While there are multiple pathways for the leakage of sequestered CO 2 from subsurface storage reservoirs (e.g.…”

Section: Introductionmentioning

confidence: 99%

“…Therefore, it is important to be continually developing the efficiency of simulation tools. Herein, a fixed point type iterative global pressure solution modification to the semi-analytical CO 2 leakage algorithm [7,36,40] is presented and explored. This work first includes a detailed description of the original semi-analytical leakage algorithm then presents the methodology for applying the proposed modification.…”

Section: Introductionmentioning

confidence: 99%

See 4 more Smart Citations

An iterative global pressure solution for the semi-analytical simulation of geological carbon sequestration

Baù

Cody²,

González‐Nicolás

2015

Comput Geosci

View full text Add to dashboard Cite

Successful large-scale implementation of geological CO 2 sequestration (GCS) will require the preliminary assessment of multiple potential injection sites. Risk assessment and optimization tools used in this effort typically require large numbers of simulations. This makes it important to choose the appropriate level of complexity when selecting the type of simulation model. A promising multiphase semi-analytical method proposed by [40] to estimate key system attributes (i.e. pressure distribution, CO 2 plume extent, and fluid migration) has been found to reduce computational run times by three orders of magnitude when compared to other standard numerical techniques. The premise of the work presented herein is that the existing semi-analytically leakage algorithm proposed by [40] may be further improved in computational efficiency by applying a fixed point type iterative global pressure solution to eliminate the need to solve large sets of linear equations at each time step. Results show that significant gains in computational efficiency are obtained with this new methodology. In addition, this modification provides the same enhancement to similar semi-analytical algorithms that simulate singe-phase injection into multi-layer domains.

show abstract

“…Therefore, the first part of this section provides a detailed description of work presented in [40] and [7].…”

Section: Methodsmentioning

confidence: 99%

Section: The Estimating Leakage Semi-analytically (Elsa) Algorithmmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

An iterative global pressure solution for the semi-analytical simulation of geological carbon sequestration

Baù

Cody²,

González‐Nicolás

2015

Comput Geosci

View full text Add to dashboard Cite

show abstract

Flow in horizontally anisotropic multilayered aquifer systems with leaky wells and aquitards

et al. 2014

View full text Add to dashboard Cite

[1] Flow problems in an anisotropic domain can be transformed into ones in an equivalent isotropic domain by coordinate transformations. Once analytical solutions are obtained for the equivalent isotropic domain, they can be back transformed to the original anisotropic domain. The existing solutions presented by Cihan et al. (2011) for isotropic multilayered aquifer systems with alternating aquitards and multiple injection/pumping wells and leaky wells were modified to account for horizontal anisotropy in aquifers. The modified solutions for pressure buildup distribution and leakage rates through leaky wells can be used when the anisotropy direction and ratio (K x =K y ) are assumed to be identical for all aquifers alternating with aquitards. However, for multilayered aquifers alternating with aquicludes, both the principal direction of the anisotropic horizontal conductivity and the anisotropy ratio can be different in each aquifer. With coordinate transformation, a circular well with finite radius becomes an ellipse, and thus in the transformed domain the head contours in the immediate vicinity of the well have elliptical shapes. Through a radial flow approximation around the finite radius wells, the elliptical well boundaries in the transformed domain are approximated by an effective well radius expression. The analytical solutions with the effective radius approximations were compared with exact solutions as well as a numerical solution for elliptic flow. The effective well radius approximation is sufficiently accurate to predict the head buildup at the well bore of the injection/pumping wells for moderately anisotropic systems (K x =K y 25). The effective radius approximation gives satisfactory results for predicting head buildup at observation points and leakage through leaky wells away from the injection/pumping wells even for highly anisotropic aquifer systems ðK x =K y 1000Þ.Citation: Cihan, A., Q. Zhou, J. T. Birkholzer, and S. R. Kraemer (2014), Flow in horizontally anisotropic multilayered aquifer systems with leaky wells and aquitards, Water Resour. Res., 50,[741][742][743][744][745][746][747]

show abstract

Analytical solutions for two‐phase subsurface flow to a leaky fault considering vertical flow effects and fault properties

Kang

Nordbotten

Doster

et al. 2014

Water Resources Research

Self Cite

View full text Add to dashboard Cite

Conductive faults in geologic formations can serve as leakage pathways for fluids such as CO 2 and methane, which would otherwise remain trapped beneath low-permeability layers. To estimate leakage rates and the associated pressure effects on adjacent aquifers, modeling of the flow in the vicinity of leaky faults must be performed. The flow from an aquifer to a leaky fault is controlled by aquifer properties as well as fault properties, including the fault permeability, fault width, and anisotropy in the fault permeability. Depending on aquifer properties and leakage rates, vertical flow equilibrium may not be reached in the aquifer in the vicinity of the leaky fault. Therefore, analytical solutions for leakage through faults that allow for vertical flow need to be developed. To do this, the vertically integrated form of the two-phase flow formulation based on vertical equilibrium is expanded by assuming a linearly structured vertical flow field, and steady state analytical solutions for single-phase and two-phase leakage are derived. In the case of twophase leakage, anisotropic aquifer permeabilities are shown to produce stronger vertical flow effects than in the case of single-phase leakage. Using numerical simulations, a model representing fault properties that is compatible with the analytical solutions is developed. The combination of the fault model and the analytical solutions captures the effects of leakage through faults with different properties and vertical flow effects. The incorporation of this fault model/analytical solution in larger basin-wide multiscale models can enable subgrid-scale effects due to leakage through faults to be captured with improved efficiency.

show abstract

Field-scale application of a semi-analytical model for estimation of CO2 and brine leakage along old wells

Cited by 138 publications

References 34 publications

An iterative global pressure solution for the semi-analytical simulation of geological carbon sequestration

An iterative global pressure solution for the semi-analytical simulation of geological carbon sequestration

Flow in horizontally anisotropic multilayered aquifer systems with leaky wells and aquitards

Analytical solutions for two‐phase subsurface flow to a leaky fault considering vertical flow effects and fault properties

Contact Info

Product

Resources

About