Publicado 6 números por año
ISSN Imprimir: 2152-5080
ISSN En Línea: 2152-5099
Indexed in
USING PARALLEL MARKOV CHAIN MONTE CARLO TO QUANTIFY UNCERTAINTIES IN GEOTHERMAL RESERVOIR CALIBRATION
SINOPSIS
We introduce a parallel rejection scheme to give a simple but reliable way to parallelize the Metropolis-Hastings algorithm. This method can be particularly useful when the target density is computationally expensive to evaluate and the acceptance rate of the Metropolis-Hastings is low. We apply the resulting method to quantify uncertainties of inverse problems, in which we aim to calibrate a challenging nonlinear geothermal reservoir model using real measurements from well tests. We demonstrate the parallelized method on various well-test scenarios. In some scenarios, the sample-based statistics obtained by our scheme shows clear advantages in providing robust model calibration and prediction compared with those obtained by nonlinear optimization methods.
-
Higdon, D., Lee, H., and Zhuoxin, B., A Bayesian Approach to Characterizing Uncertainty in Inverse Problems Using Coarse and Fine-Scale Information, IEEE Trans. Signal Proc, 50(2):389-399, 2002.
-
Andersen, K.E., Brooks, S.P., and Hansen, M.B., Bayesian Inversion of Geoelectrical Resistivity Data, J. R. Stat. Soc., Ser. B, 65:619-644,2004.
-
Cui, T., Fox, C., and O'Sullivan, M.J., Bayesian Calibration of a Large-Scale Geothermal Reservoir Model by a New Adaptive Delayed Acceptance Metropolis Hastings Algorithm, Water Res. Res., 47:W10521, 2011.
-
Nicholls, G.K. and Fox, C., Prior Modelling and Posterior Sampling in Impedance Imaging, Proc. SPIE, 3459:116-127,1998.
-
Vauhkonen, P., Vauhkonen, M., Savolainen, T., and Kaipio, J., Three-Dimensional Electrical Impedance Tomography based on the Complete Electrode Model, IEEE Trans. Biomed. Eng., 46:1150-1160, 1999.
-
Watzenig, D. and Fox, C., A Review of Statistical Modelling and Inference for Electrical Capacitance Tomography, Meas. Sci. Technol., 20(5):052002, 2009.
-
McKeague, I.W., Nicholls, K., G., Speer, K., and Herbei, R., Statistical Inversion of South Atlantic Circulation in an Abyssal Neutral Density Layer, J. Marine Res, 63(4):683-704, 2005.
-
Haario, H., Laine, M., Lehtinen, M., Saksman, E., and Tamminen, J., Markov Chain Monte Carlo Methods for High Dimensional Inversion in Remote Sensing, J. R. Stat. Soc., Ser. B, 66:591-608, 2004.
-
Martin, J., Wilcox, L.C., Burstedde, C., and Ghattas, O., A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion, SIAMJ. Sci. Comput., 34(3):A1460-A1487, 2012.
-
Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., and Teller, E., Equation of State Calculations by Fast Computing Machines, J. Chem. Phys, 21:1087-1092, 1953.
-
Hastings, W., Monte Carlo Sampling Using Markov Chains and Their Applications, Biometrika, 57:97-109, 1970.
-
Geyer, C.J., Markov Chain Monte Carlo Maximum Likelihood, in Computing Science And Statistics: Proc. of 23rd Symposium on the Interface, E. Keramigas, Ed., 5:156-163, Fairfax Station, VA: Interface Foundation of North America, 1991.
-
Chen, Y., Keyes, D., Law, K.J.H., and Ltaief, H., Accelerated Dimension-Independent Adaptive Metropolis, SIAM J. Sci. Comput, 38(5):S539-S565, 2016.
-
Craiu, R.V., Rosenthal, J., and Yang, C., Learn from Thy Neighbor: Parallel-Chain and Regional Adaptive MCMC, J. Am. Stat. Assoc, 104(488):1454-1466, 2009.
-
Solonen, A., Ollinaho, P., Laine, M., Haario, H., Tamminen, J., and Jarvinen, H., Efficient MCMC for Climate Model Parameter Estimation: Parallel Adaptive Chains and Early Rejection, Bayesian Anal., 7(3):715-736,2012.
-
Calderhead, B., A General Construction for Parallelizing Metropolis-Hastings Algorithms, Proc. Natl. Acad. Sci, 111(49):17408-17413, 2014.
-
Finsterle, S., ITOUGH2 User's Guide, Berkeley, CA: Lawrence Berkeley National Laboratory, 1993.
-
Doherty, J., Groundwater Model Calibration Using Pilot Points and Regularization, Ground Water, 41:170-177,2003.
-
Marquardt, D.W., An Algorithm for Least Square Estimation of Nonlinear Parameters, J. Soc. Ind. Appl. Math., 11(2):431- 441, 1963.
-
Roberts, G.O. and Stramer, S., Langevin Diffusions and Metropolis-Hastings Algorithms, Methodol. Comput. Appl. Probl., 4(4):337-357, 2002.
-
Braak, C.J.F.T., A Markov Chain Monte Carlo Version of the Genetic Algorithm Differential Evolution: Easy Bayesian Computing for Real Parameter Spaces, Stat. Comput, 16:239-249, 2006.
-
Finsterle, S., Pruess, K., Bullivant, D., and O'Sullivan, M.J., Application of Inverse Modeling to Geothermal Reservoir Simulation, in Proc. of Twenty-Second Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, CA, pp. 309-316, 1997.
-
Liu, J.S., Monte Carlo Strategies in Scientific Computing, Berlin: Springer-Verlag, 2001.
-
Roberts, G.O. and Rosenthal, J.S., Optimal Scaling for Various Metropolis-Hastings Algorithms, Stat. Sci., 16:351-367,2001.
-
Roberts, G.O., Gelman, A., and Gilks, W.R., Weak Convergence and Optimal Scaling of Random Walk Metropolis Algorithms, Ann. Appl. Probab, 7:110-120, 1997.
-
Sohn, A., Parallel Speculative Computation of Simulated Annealing, IEEE Trans. Parallel Distributed Syst., 6(10):997-1005, 1995.
-
Mascagni, M. and Srinivasan, A., Algorithm 806: SPRNG: A Scalable Library for Pseudorandom Number Generation, ACM Trans. Math. Software, 26:436-461, 2000.
-
Efendiev, Y., Thomas, H., and Luo, W., Preconditioning Markov Chain Monte Carlo Simulations Using Coarse-Scale Models, SIAMJ Sci. Comput., 28:776-803,2006.
-
van Genuchten, M.T., A Closed-Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils, Soil Sci. Soc. Am. J., 44:892-898,1980.
-
International Formulation Committee, A Formulation of the Thermodynamic Properties of Ordinary Water Substance, Dusseldorf, Germany: IFC Secretariat, 1967.
-
Pruess, K., TOUGH2-A General-Purpose Numerical Simulator for Multiphase Fluid and Heat Flow, Berkeley, CA: Lawrence Berkeley National Laboratory, 1991.
-
Jaynes, E.T., Prior Probability, IEEE Trans. Syst. Sci. Cybern., 4(3):227-241, 1968.
-
Haario, H., Saksman, E., and Tamminen, J., An Adaptive Metropolis Algorithm, Bernoulli, 7:223-242, 2001.
-
Narayan, A., Jakeman, J., and Zhou, T., A Christoffel Function Weighted Least Squares Algorithm for Collocation Approximations, Math. Comput., 86:1913-1947, 2017.
-
Yan, L. and Zhou, T., An Adaptive Multi-Fidelity PC-Based Ensemble Kalman Inversion for Inverse Problems, Int. J. Uncertainty Quantif, 9:205-220, 2019.
-
Yan, L. and Zhou, T., Adaptive Multi-Fidelity Polynomial Chaos Approach to Bayesian Inference in Inverse Problems, J. Comput. Phys., 381:110-128,2019.
-
Cui, T., Marzouk, Y.M., and Willcox, K.E., Data-Driven Model Reduction for the Bayesian Solution of Inverse Problems, Int. J. Numer. Methods Eng., 102:966-990, 2015.
-
Cui, T., Marzouk, Y.M., and Willcox, K.E., Scalable Posterior Approximations for Large-Scale Bayesian Inverse Problems via Likelihood-Informed Parameter and State Reduction, J. Comput. Phys., 315:363-387, 2016.
-
Maclaren Oliver J., Nicholson Ruanui, Bjarkason Elvar K., O'Sullivan John P., O'Sullivan Michael J., Incorporating Posterior‐Informed Approximation Errors Into a Hierarchical Framework to Facilitate Out‐of‐the‐Box MCMC Sampling for Geothermal Inverse Problems and Uncertainty Quantification, Water Resources Research, 56, 1, 2020. Crossref
-
Scott S. W., O’Sullivan J. P., Maclaren O. J., Nicholson R., Covell C., Newson J., Guðjónsdóttir M. S., Bayesian Calibration of a Natural State Geothermal Reservoir Model, Krafla, North Iceland, Water Resources Research, 58, 2, 2022. Crossref