Fast stokes flow simulations for geophysical-geodynamic inverse problems and sensitivity analyses based on reduced order modeling
Visualitza/Obre
Cita com:
hdl:2117/335920
Tipus de documentArticle
Data publicació2020-03-14
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
ProjecteMATHROCKS - Multiscale Inversion of Porous Rock Physics using High-Performance Simulators: Bridging the Gap between Mathematics and Geophysics (EC-H2020-777778)
ASIMILACION DE DATOS PARA UNA SIMULACION INGENIERIL CREIBLE (AEI-DPI2017-85139-C2-2-R)
ASIMILACION DE DATOS PARA UNA SIMULACION INGENIERIL CREIBLE (AEI-DPI2017-85139-C2-2-R)
Abstract
Markov chain Monte Carlo (MCMC) methods have become standard in Bayesian inference and multi-observable inversions in almost every discipline of the Earth sciences. In the case of geodynamic and/or coupled geophysical-geodynamic inverse problems, however, the computational cost associated with the solution of large-scale 3-D Stokes forward problems has rendered probabilistic formulations impractical. Here we present a novel and extremely efficient method to produce ultrafast solutions of the 3-D Stokes problem for MCMC simulations. Our approach combines the individual benefits of Reduced Basis techniques, goal-oriented error formulations, and MCMC algorithms to produce an accurate and computationally efficient surrogate for the forward problem. Importantly, the surrogate adapts itself during the MCMC simulation according to the history of the chain and the goals of the inversion. This maximizes the efficiency of the forward problem and removes the need for preinversion off-line computations to build a surrogate. We demonstrate the benefits and limitations of the method with several numerical examples and show that in all cases the computational cost is of the order of <1% compared to a traditional MCMC approach. The method is general enough to be applied to a range of problems, including uncertainty quantification/propagation, adjoint-based geodynamic inversions, sensitivity analyses in mantle convection problems, and in the creating surrogate models for complex forward problems (e.g., heat transfer, seismic tomography, and magnetotellurics).
CitacióOrtega, O. [et al.]. Fast stokes flow simulations for geophysical-geodynamic inverse problems and sensitivity analyses based on reduced order modeling. "Journal of geophysical research B: solid earth", 14 Març 2020, vol. 125, núm. 3, p. e2019JB018314:1-e2019JB018314:25.
ISSN2169-9356
Versió de l'editorhttps://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2019JB018314
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
2020-JGR-OZAD.pdf | 6,644Mb | Visualitza/Obre |