Estimating the probability of CO2 leakage using rare event simulation

A.h. Elsheikh, S. Oladyshkin, W. Nowak, M. Christie

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Estimating the probability of rare events is an extremely challenging task. For example, estimating the probability of leakage of CO2 from a saline aquifer using a direct Monte-Carlo approach would in general require a number of simulations proportional to the inverse of the rare event probability. Since it is likely that any action will require a probability of failure of less than $10^{-6}$, requiring $10^7$ to $10^8$ simulations, it is understandable that few such studies have been published. In this paper, we propose a means of simulating such rare events using a multilevel splitting algorithm called subset simulation (SS)[Au and Beck, 2001]. SS is an iterative algorithm that introduces intermediate events which are easier to sample from, and then iteratively samples within each constrained region in the probability space until the rare event threshold is reached. We show results for a standard benchmark for CO2 leakage through a leaky well [Class et al., 2009]. In this test case, CO2 is injected into a deep aquifer, then spreads within the aquifer and, upon reaching an abandoned well, it rises to a shallower aquifer. We show that subset simulation is an effective algorithm for estimating the probability of rare events with significant computational advantages over a direct Monte-Carlo approach. The SS algorithm relies on two parameters that have to be adjusted at the start of the simulation; we show the effect of these parameters on the quality of estimated rare event probabilities and propose guidelines on how to select these parameters.
Original languageEnglish
Title of host publication14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV)
PublisherEAGE Publishing BV
ISBN (Print)9781634391689
DOIs
Publication statusPublished - 8 Sep 2014
Event14th European Conference on the Mathematics of Oil Recovery 2014 - Catania, Italy
Duration: 8 Sep 201411 Sep 2014

Conference

Conference14th European Conference on the Mathematics of Oil Recovery 2014
Abbreviated titleECMOR 2014
CountryItaly
CityCatania
Period8/09/1411/09/14

Fingerprint

leakage
Set theory
Aquifers
simulation
aquifer
Abandoned wells
parameter

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geotechnical Engineering and Engineering Geology
  • Energy Engineering and Power Technology

Cite this

Elsheikh, A. H., Oladyshkin, S., Nowak, W., & Christie, M. (2014). Estimating the probability of CO2 leakage using rare event simulation. In 14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV) EAGE Publishing BV. https://doi.org/10.3997/2214-4609.20141876
Elsheikh, A.h. ; Oladyshkin, S. ; Nowak, W. ; Christie, M. / Estimating the probability of CO2 leakage using rare event simulation. 14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV). EAGE Publishing BV, 2014.
@inproceedings{8273e11d8be945acbe183527080e129a,
title = "Estimating the probability of CO2 leakage using rare event simulation",
abstract = "Estimating the probability of rare events is an extremely challenging task. For example, estimating the probability of leakage of CO2 from a saline aquifer using a direct Monte-Carlo approach would in general require a number of simulations proportional to the inverse of the rare event probability. Since it is likely that any action will require a probability of failure of less than $10^{-6}$, requiring $10^7$ to $10^8$ simulations, it is understandable that few such studies have been published. In this paper, we propose a means of simulating such rare events using a multilevel splitting algorithm called subset simulation (SS)[Au and Beck, 2001]. SS is an iterative algorithm that introduces intermediate events which are easier to sample from, and then iteratively samples within each constrained region in the probability space until the rare event threshold is reached. We show results for a standard benchmark for CO2 leakage through a leaky well [Class et al., 2009]. In this test case, CO2 is injected into a deep aquifer, then spreads within the aquifer and, upon reaching an abandoned well, it rises to a shallower aquifer. We show that subset simulation is an effective algorithm for estimating the probability of rare events with significant computational advantages over a direct Monte-Carlo approach. The SS algorithm relies on two parameters that have to be adjusted at the start of the simulation; we show the effect of these parameters on the quality of estimated rare event probabilities and propose guidelines on how to select these parameters.",
author = "A.h. Elsheikh and S. Oladyshkin and W. Nowak and M. Christie",
year = "2014",
month = "9",
day = "8",
doi = "10.3997/2214-4609.20141876",
language = "English",
isbn = "9781634391689",
booktitle = "14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV)",
publisher = "EAGE Publishing BV",
address = "Netherlands",

}

Elsheikh, AH, Oladyshkin, S, Nowak, W & Christie, M 2014, Estimating the probability of CO2 leakage using rare event simulation. in 14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV). EAGE Publishing BV, 14th European Conference on the Mathematics of Oil Recovery 2014, Catania, Italy, 8/09/14. https://doi.org/10.3997/2214-4609.20141876

Estimating the probability of CO2 leakage using rare event simulation. / Elsheikh, A.h.; Oladyshkin, S.; Nowak, W.; Christie, M.

14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV). EAGE Publishing BV, 2014.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

TY - GEN

T1 - Estimating the probability of CO2 leakage using rare event simulation

AU - Elsheikh, A.h.

AU - Oladyshkin, S.

AU - Nowak, W.

AU - Christie, M.

PY - 2014/9/8

Y1 - 2014/9/8

N2 - Estimating the probability of rare events is an extremely challenging task. For example, estimating the probability of leakage of CO2 from a saline aquifer using a direct Monte-Carlo approach would in general require a number of simulations proportional to the inverse of the rare event probability. Since it is likely that any action will require a probability of failure of less than $10^{-6}$, requiring $10^7$ to $10^8$ simulations, it is understandable that few such studies have been published. In this paper, we propose a means of simulating such rare events using a multilevel splitting algorithm called subset simulation (SS)[Au and Beck, 2001]. SS is an iterative algorithm that introduces intermediate events which are easier to sample from, and then iteratively samples within each constrained region in the probability space until the rare event threshold is reached. We show results for a standard benchmark for CO2 leakage through a leaky well [Class et al., 2009]. In this test case, CO2 is injected into a deep aquifer, then spreads within the aquifer and, upon reaching an abandoned well, it rises to a shallower aquifer. We show that subset simulation is an effective algorithm for estimating the probability of rare events with significant computational advantages over a direct Monte-Carlo approach. The SS algorithm relies on two parameters that have to be adjusted at the start of the simulation; we show the effect of these parameters on the quality of estimated rare event probabilities and propose guidelines on how to select these parameters.

AB - Estimating the probability of rare events is an extremely challenging task. For example, estimating the probability of leakage of CO2 from a saline aquifer using a direct Monte-Carlo approach would in general require a number of simulations proportional to the inverse of the rare event probability. Since it is likely that any action will require a probability of failure of less than $10^{-6}$, requiring $10^7$ to $10^8$ simulations, it is understandable that few such studies have been published. In this paper, we propose a means of simulating such rare events using a multilevel splitting algorithm called subset simulation (SS)[Au and Beck, 2001]. SS is an iterative algorithm that introduces intermediate events which are easier to sample from, and then iteratively samples within each constrained region in the probability space until the rare event threshold is reached. We show results for a standard benchmark for CO2 leakage through a leaky well [Class et al., 2009]. In this test case, CO2 is injected into a deep aquifer, then spreads within the aquifer and, upon reaching an abandoned well, it rises to a shallower aquifer. We show that subset simulation is an effective algorithm for estimating the probability of rare events with significant computational advantages over a direct Monte-Carlo approach. The SS algorithm relies on two parameters that have to be adjusted at the start of the simulation; we show the effect of these parameters on the quality of estimated rare event probabilities and propose guidelines on how to select these parameters.

UR - http://www.scopus.com/inward/record.url?scp=84908213257&partnerID=8YFLogxK

U2 - 10.3997/2214-4609.20141876

DO - 10.3997/2214-4609.20141876

M3 - Conference contribution

SN - 9781634391689

BT - 14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV)

PB - EAGE Publishing BV

ER -

Elsheikh AH, Oladyshkin S, Nowak W, Christie M. Estimating the probability of CO2 leakage using rare event simulation. In 14th European Conference on the Mathematics of Oil Recovery 2014 (ECMOR XIV). EAGE Publishing BV. 2014 https://doi.org/10.3997/2214-4609.20141876