A comparison of deterministic and stochastic approaches for regional scale inverse modeling on the Mar del Plata aquifer

被引:36
|
作者
Pool, M. [1 ,2 ]
Carrera, J. [1 ,2 ]
Alcolea, A. [3 ]
Bocanegra, E. M. [4 ]
机构
[1] CSIC, Inst Environm Assessment & Water Res IDAEA, ES-08034 Barcelona, Spain
[2] UPC, CSIC, Hydrogeol Grp, Associated Unit, Barcelona, Spain
[3] TK Consult AG, Zurich, Switzerland
[4] Univ Nacl Mar del Plata, Inst Geol Costas & Cuaternario, Mar Del Plata, Buenos Aires, Argentina
关键词
Real regional aquifer; Deterministic-stochastic approaches; Pilot points; Plausible calibration; 3-DIMENSIONAL NUMERICAL INVERSION; GROUND-PENETRATING-RADAR; CROSS-HOLE TESTS; HYDRAULIC TOMOGRAPHY; STEADY-STATE; GEOSTATISTICAL INVERSION; CONDITIONAL SIMULATION; SUBSURFACE HETEROGENEITY; TRANSMISSIVITY FIELDS; SPATIAL VARIABILITY;
D O I
10.1016/j.jhydrol.2015.09.064
中图分类号
TU [建筑科学];
学科分类号
0813 ;
摘要
Inversion of the spatial variability of transmissivity (T) in groundwater models can be handled using either stochastic or deterministic (i.e., geology-based zonation) approaches. While stochastic methods predominate in scientific literature, they have never been formally compared to deterministic approaches, preferred by practitioners, for regional aquifer models. We use both approaches to model groundwater flow and solute transport in the Mar del Plata aquifer, where seawater intrusion is a major threat to freshwater resources. The relative performance of the two approaches is evaluated in terms of (i) model fits to head and concentration data (available for nearly a century), (ii) geological plausibility of the estimated T fields, and (iii) their ability to predict transport. We also address the impact of conditioning the estimated fields on T data coming from either pumping tests interpreted with the Theis method or specific capacity values from step-drawdown tests. We find that stochastic models, based upon conditional estimation and simulation techniques, identify some of the geological features (river deposit channels and low transmissivity regions associated to quartzite outcrops) and yield better fits to calibration data than the much simpler geology-based deterministic model, which cannot properly address model structure uncertainty. However, the latter demonstrates much greater robustness for predicting sea water intrusion and for incorporating concentrations as calibration data. We attribute the poor performance, and underestimated uncertainty, of the stochastic simulations to estimation bias introduced by model errors. Qualitative geological information is extremely rich in identifying large-scale variability patterns, which are identified by stochastic models only in data rich areas, and should be explicitly included in the calibration process. (C) 2015 Elsevier B.V. All rights reserved.
引用
收藏
页码:214 / 229
页数:16
相关论文
共 50 条