Stress estimation in reservoirs using an integrated inverse method

Antoine Mazuyer and Paul Cupillard and Richard Giot and Marianne Conin and Y. Leroy and Pierre Thore. ( 2018 )
in: Computers & Geosciences, 114 (30-40)

Abstract

Estimating the stress in reservoirs and their surroundings prior to the production is a key issue for reservoir management planning. In this study, we propose an integrated inverse method to estimate such initial stress state. The 3D stress state is constructed with the displacement-based finite element method assuming linear isotropic elasticity and small perturbations in the current geometry of the geological structures. The Neumann boundary conditions are defined as piecewise linear functions of depth. The discontinuous functions are determined with the CMA-ES (Covariance Matrix Adaptation Evolution Strategy) optimization algorithm to fit wellbore stress data deduced from leak-off tests and breakouts. The disregard of the geological history and the simplified rheological assumptions mean that only the stress field, statically admissible and matching the wellbore data should be exploited. The spatial domain of validity of this statement is assessed by comparing the stress estimations for a synthetic folded structure of finite amplitude with a history constructed assuming a viscous response.

Download / Links

BibTeX Reference

@ARTICLE{,
    author = { Mazuyer, Antoine and Cupillard, Paul and Giot, Richard and Conin, Marianne and Leroy, Y. and Thore, Pierre },
     title = { Stress estimation in reservoirs using an integrated inverse method },
   journal = { Computers & Geosciences },
    volume = { 114 },
      year = { 2018 },
     pages = { 30-40 },
       url = { https://www.sciencedirect.com/science/article/pii/S0098300417305010 },
       doi = { 10.1016/j.cageo.2018.01.004 },
  abstract = { Estimating the stress in reservoirs and their surroundings prior to the production is a key issue for reservoir management planning. In this study, we propose an integrated inverse method to estimate such initial stress state. The 3D stress state is constructed with the displacement-based finite element method assuming linear isotropic elasticity and small perturbations in the current geometry of the geological structures. The Neumann boundary conditions are defined as piecewise linear functions of depth. The discontinuous functions are determined with the CMA-ES (Covariance Matrix Adaptation Evolution Strategy) optimization algorithm to fit wellbore stress data deduced from leak-off tests and breakouts. The disregard of the geological history and the simplified rheological assumptions mean that only the stress field, statically admissible and matching the wellbore data should be exploited. The spatial domain of validity of this statement is assessed by comparing the stress estimations for a synthetic folded structure of finite amplitude with a history constructed assuming a viscous response. }
}