An efficient surrogate-based simulation-optimization method for calibrating a regional MODFLOW model

Abstract Simulation-optimization method entails a large number of model simulations, which is computationally intensive or even prohibitive if the model simulation is extremely time-consuming. Statistical models have been examined as a surrogate of the high-fidelity physical model during simulation-optimization process to tackle this problem. Among them, Multivariate Adaptive Regression Splines (MARS), a non-parametric adaptive regression method, is superior in overcoming problems of high-dimensions and discontinuities of the data. Furthermore, the stability and accuracy of MARS model can be improved by bootstrap aggregating methods, namely, bagging. In this paper, Bagging MARS (BMARS) method is integrated to a surrogate-based simulation-optimization framework to calibrate a three-dimensional MODFLOW model, which is developed to simulate the groundwater flow in an arid hardrock-alluvium region in northwestern Oman. The physical MODFLOW model is surrogated by the statistical model developed using BMARS algorithm. The surrogate model, which is fitted and validated using training dataset generated by the physical model, can approximate solutions rapidly. An efficient Sobol’ method is employed to calculate global sensitivities of head outputs to input parameters, which are used to analyze their importance for the model outputs spatiotemporally. Only sensitive parameters are included in the calibration process to further improve the computational efficiency. Normalized root mean square error (NRMSE) between measured and simulated heads at observation wells is used as the objective function to be minimized during optimization. The reasonable history match between the simulated and observed heads demonstrated feasibility of this high-efficient calibration framework.

[1]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[2]  Ilya M. Sobol,et al.  Theorems and examples on high dimensional model representation , 2003, Reliab. Eng. Syst. Saf..

[3]  R. Dennis Cook,et al.  Cross-Validation of Regression Models , 1984 .

[4]  Hari S. Viswanathan,et al.  An Integrated Framework for Optimizing CO2 Sequestration and Enhanced Oil Recovery , 2014 .

[5]  J. Friedman.  Multivariate adaptive regression splines , 1990 .

[6]  Farrokh Mistree,et al.  Kriging Models for Global Approximation in Simulation-Based Multidisciplinary Design Optimization , 2001 .

[7]  Rajesh J. Pawar,et al.  A response surface model to predict CO 2 and brine leakage along cemented wellbores , 2015 .

[8]  M. Powell.  The BOBYQA algorithm for bound constrained optimization without derivatives , 2009 .

[9]  Bryan A. Tolson,et al.  Review of surrogate modeling in water resources , 2012 .

[10]  Javier Samper,et al.  Inverse problem of multicomponent reactive chemical transport in porous media: Formulation and applications , 2004 .

[11]  Robert J. Mellors,et al.  An efficient optimization of well placement and control for a geothermal prospect under geological uncertainty , 2015 .

[12]  Rajesh J. Pawar,et al.  Reduced order models for assessing CO 2 impacts in shallow unconfined aquifers , 2016 .

[13]  Dongxiao Zhang,et al.  A sparse grid based Bayesian method for contaminant source identification , 2012 .

[14]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[15]  Abelardo Ramirez,et al.  An efficient Bayesian inversion of a geothermal prospect using a multivariate adaptive regression spline method , 2014 .

[16]  J. Samper,et al.  Inverse modeling of water flow and multicomponent reactive transport in coastal aquifer systems , 2006 .

[17]  Wenke Wang,et al.  An ecology-oriented exploitation mode of groundwater resources in the northern Tianshan Mountains, China , 2016 .

[18]  A. Izady,et al.  An efficient methodology to design optimal groundwater level monitoring network in Al-Buraimi region, Oman , 2017, Arabian Journal of Geosciences.

[19]  Yunwei Sun,et al.  Surrogate-based optimization of hydraulic fracturing in pre-existing fracture networks , 2013, Comput. Geosci..

[20]  Christine A. Shoemaker,et al.  A Stochastic Radial Basis Function Method for the Global Optimization of Expensive Functions , 2007, INFORMS J. Comput..

[21]  Leah L. Rogers,et al.  Solving Problems in Environmental Engineering and Geosciences with Artificial Neural Networks , 1996 .

[22]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[23]  B. Wagner.  Recent advances in simulation-optimization groundwater management modeling (95RG00394) , 1995 .

[24]  Virginia M. Johnson,et al.  Applying soft computing methods to improve the computational tractability of a subsurface simulation–optimization problem , 2001 .

[25]  Raghavan Srinivasan,et al.  Approximating SWAT Model Using Artificial Neural Network and Support Vector Machine 1 , 2009 .

[26]  Agus Sudjianto,et al.  Computer Aided Reliability and Robustness Assessment , 1998 .

[27]  G. Gary Wang,et al.  Review of Metamodeling Techniques in Support of Engineering Design Optimization , 2007 .

[28]  T. Simpson,et al.  Comparative studies of metamodelling techniques under multiple modelling criteria , 2001 .

[29]  A. W. Harbaugh.  MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[30]  M. Wang,et al.  Optimal Remediation Policy Selection under General Conditions , 1997 .

[31]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[32]  Richard G. Niswonger,et al.  MODFLOW-NWT, A Newton Formulation for MODFLOW-2005 , 2014 .

[33]  Donald D. Lucas,et al.  The parametric sensitivity of CAM5's MJO , 2014 .

[34]  P. Bühlmann.  Bagging, subagging and bragging for improving some prediction algorithms , 2003 .

[35]  Hsien-Chie Cheng,et al.  Assessing a Response Surface-Based Optimization Approach for Soil Vapor Extraction System Design , 2009 .

[36]  Leo Breiman,et al.  Bagging Predictors , 1996, Machine Learning.

[37]  I. Sobola,et al.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .