Efficient parameter estimation for a methane hydrate model with active subspaces
arXiv:1801.09499 · doi:10.1007/s10596-018-9769-x
Abstract
Methane gas hydrates have increasingly become a topic of interest because of their potential as a future energy resource. There are significant economical and environmental risks associated with extraction from hydrate reservoirs, so a variety of multiphysics models have been developed to analyze prospective risks and benefits. These models generally have a large number of empirical parameters which are not known a priori. Traditional optimization-based parameter estimation frameworks may be ill-posed or computationally prohibitive. Bayesian inference methods have increasingly been found effective for estimating parameters in complex geophysical systems. These methods often are not viable in cases of computationally expensive models and high-dimensional parameter spaces. Recently, methods have been developed to effectively reduce the dimension of Bayesian inverse problems by identifying low-dimensional structures that are most informed by data. Active subspaces is one of the most generally applicable methods of performing this dimension reduction. In this paper, Bayesian inference of the parameters of a state-of-the-art mathematical model for methane hydrates based on experimental data from a triaxial compression test with gas hydrate-bearing sand is performed in an efficient way by utilizing active subspaces. Active subspaces are used to identify low-dimensional structure in the parameter space which is exploited by generating a cheap regression-based surrogate model and implementing a modified Markov chain Monte Carlo algorithm. Posterior densities having means that match the experimental data are approximated in a computationally efficient way.
arXiv comment21 pages, 15 figures
References in corpus (11)
- Scikit-learn: Machine Learning in Python
- MCMC using Hamiltonian dynamics
- Active subspace methods in theory and practice: applications to kriging surfaces
- Dimension-independent likelihood-informed MCMC
- Spectral gaps for a Metropolis-Hastings algorithm in infinite dimensions
- Non-Isothermal, Multi-phase, Multi-component Flows through Deformable Methane Hydrate Reservoirs
- Accelerating MCMC with active subspaces
- A Consistent Bayesian Formulation for Stochastic Inverse Problems Based on Push-forward Measures
- Solving Large-Scale PDE-constrained Bayesian Inverse Problems with Riemann Manifold Hamiltonian Monte Carlo
- Computing active subspaces with Monte Carlo
- Testing a coupled hydro-thermo-chemo-geomechanical model for gas hydrate bearing sediments using triaxial compression lab experiments