Generalized Polynomial Chaos Expansion for Fast and Accurate Uncertainty Quantification in Geomechanical Modelling

Geomechanical modelling of the processes associated to the exploitation of subsurface resources, such as land subsidence or triggered/induced seismicity, is a common practice of major interest. The prediction reliability depends on different sources of uncertainty, such as the parameterization of the constitutive model characterizing the deep rock behaviour. In this study, we focus on a Sobol’-based sensitivity analysis and uncertainty reduction via assimilation of land deformations. A synthetic test case application on a deep hydrocarbon reservoir is considered, where land settlements are predicted with the aid of a 3-D Finite Element (FE) model. Data assimilation is performed via the Ensemble Smoother (ES) technique and its variation in the form of Multiple Data Assimilation (ES-MDA). However, the ES convergence is guaranteed with a large number of Monte Carlo (MC) simulations, that may be computationally infeasible in large scale and complex systems. For this reason, a surrogate model based on the generalized Polynomial Chaos Expansion (gPCE) is proposed as an approximation of the forward problem. This approach allows to efficiently compute the Sobol’ indices for the sensitivity analysis and greatly reduce the computational cost of the original ES and MDA formulations, also enhancing the accuracy of the overall prediction process.

[1]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

[2]  Alberto Guadagnini,et al.  Moment-based metrics for global sensitivity analysis of hydrological systems , 2017 .

[3]  Holger Class,et al.  Bayesian updating via bootstrap filtering combined with data-driven polynomial chaos expansions: methodology and application to history matching for carbon dioxide storage in geological formations , 2013, Computational Geosciences.

[4]  Y. Marzouk,et al.  A stochastic collocation approach to Bayesian inference in inverse problems , 2009 .

[5]  Michael S. Eldred,et al.  Sparse Pseudospectral Approximation Method , 2011, 1109.2936.

[6]  Nicola Castelletto,et al.  II cycle compressibility from satellite measurements , 2013 .

[7]  Carlo Janna,et al.  A two‐invariant pseudoelastic model for reservoir compaction , 2017 .

[8]  Takeshi Kodaka,et al.  An elasto‐viscoplastic model for diatomaceous mudstone and numerical simulation of compaction bands , 2011 .

[9]  M. Ferronato,et al.  Modeling fault activation due to fluid production: Bayesian update by seismic data , 2019, Computational Geosciences.

[10]  A. Reynolds,et al.  History matching time-lapse seismic data using the ensemble Kalman filter with multiple data assimilations , 2012, Computational Geosciences.

[11]  D. Xiu.  Efficient collocational approach for parametric uncertainty analysis , 2007 .

[12]  Bruno Sudret,et al.  Global sensitivity analysis using polynomial chaos expansions , 2008, Reliab. Eng. Syst. Saf..

[13]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[14]  Habib N. Najm,et al.  Uncertainty Quantification and Polynomial Chaos Techniques in Computational Fluid Dynamics , 2009 .

[15]  Dongbin Xiu,et al.  A generalized polynomial chaos based ensemble Kalman filter with high accuracy , 2009, J. Comput. Phys..

[16]  G. Evensen,et al.  Data assimilation and inverse methods in terms of a probabilistic formulation , 1996 .

[17]  Domenico Baù,et al.  Basin-scale compressibility of the northern Adriatic by the radioactive marker technique , 2002 .

[18]  Roger Ghanem,et al.  Characterization of reservoir simulation models using a polynomial chaos‐based ensemble Kalman filter , 2009 .

[19]  Gioacchino Viggiani,et al.  A general formulation of hypoplasticity , 2004 .

[20]  Ali Rezaei,et al.  A global sensitivity analysis and reduced-order models for hydraulically fractured horizontal wells , 2018, Computational Geosciences.

[21]  Youssef M. Marzouk,et al.  Adaptive Smolyak Pseudospectral Approximations , 2012, SIAM J. Sci. Comput..

[22]  Geir Evensen,et al.  Analysis of iterative ensemble smoothers for solving inverse problems , 2018, Computational Geosciences.

[23]  N. Wiener.  The Homogeneous Chaos , 1938 .

[24]  Bruno Sudret,et al.  Using sparse polynomial chaos expansions for the global sensitivity analysis of groundwater lifetime expectancy in a multi-layered hydrogeological model , 2015, Reliab. Eng. Syst. Saf..

[25]  Carlo Janna,et al.  Robust numerical implementation of a 3D rate‐dependent model for reservoir geomechanical simulations , 2019, International Journal for Numerical and Analytical Methods in Geomechanics.

[26]  R. Hanssen,et al.  Application of an ensemble smoother with multiple data assimilation to the Bergermeer gas field, using PS-InSAR , 2016 .

[27]  Domenico Baù,et al.  On the importance of the heterogeneity assumption in the characterization of reservoir geomechanical properties , 2016 .

[28]  Alexandre A. Emerick.  Analysis of the performance of ensemble-based assimilation of production and seismic data , 2016 .

[29]  A. Brovelli,et al.  A strain-rate-dependent modified Cam-Clay model for the simulation of soil/rock compaction , 2017 .

[30]  Francesca Bottazzi,et al.  A Functional Data Analysis Approach to Surrogate Modeling in Reservoir and Geomechanics Uncertainty Quantification , 2017, Mathematical Geosciences.

[31]  Giuseppe Gambolati,et al.  Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations , 2001 .

[32]  Albert C. Reynolds,et al.  Ensemble smoother with multiple data assimilation , 2013, Comput. Geosci..

[33]  Domenico Baù,et al.  Data assimilation of surface displacements to improve geomechanical parameters of gas storage reservoirs , 2016 .

[34]  Emanuele Borgonovo,et al.  A new uncertainty importance measure , 2007, Reliab. Eng. Syst. Saf..

[35]  Alberto Guadagnini,et al.  Global sensitivity analysis through polynomial chaos expansion of a basin-scale geochemical compaction model , 2013, Computational Geosciences.

[36]  Olivier P. Le Maître,et al.  Polynomial chaos expansion for sensitivity analysis , 2009, Reliab. Eng. Syst. Saf..

[37]  Vinay Prasad,et al.  Uncertainty quantification of the factor of safety in a steam-assisted gravity drainage process through polynomial chaos expansion , 2020, Comput. Chem. Eng..

[38]  Mark D. Zoback,et al.  Viscous creep in room-dried unconsolidated Gulf of Mexico shale (II): Development of a viscoplasticity model , 2010 .

[39]  Dongxiao Zhang,et al.  Stochastic analysis of unsaturated flow with probabilistic collocation method , 2009 .

[40]  Alberto Guadagnini,et al.  Use of global sensitivity analysis and polynomial chaos expansion for interpretation of nonreactive transport experiments in laboratory‐scale porous media , 2011 .

[41]  Alberto Guadagnini,et al.  Uncertainty Quantification and Global Sensitivity Analysis of Subsurface Flow Parameters to Gravimetric Variations During Pumping Tests in Unconfined Aquifers , 2018 .

[42]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[43]  Rainer Helmig,et al.  An integrative approach to robust design and probabilistic risk assessment for CO2 storage in geological formations , 2011 .

[44]  Olivier Le Maitre,et al.  Numerical approximation of poroelasticity with random coefficients using Polynomial Chaos and Hybrid High-Order methods , 2020 .

[45]  Roger Ghanem,et al.  Stochastic Finite Elements with Multiple Random Non-Gaussian Properties , 1999 .

[46]  Mary F. Wheeler,et al.  Efficient Bayesian inference of subsurface flow models using nested sampling and sparse polynomial chaos surrogates , 2014 .