Impact of geostatistical reconstruction approaches on model calibration for flow in highly heterogeneous aquifers
- PDF / 4,135,665 Bytes
- 16 Pages / 595.276 x 790.866 pts Page_size
- 70 Downloads / 211 Views
ORIGINAL PAPER
Impact of geostatistical reconstruction approaches on model calibration for flow in highly heterogeneous aquifers Martina Siena1
•
Monica Riva1
The Author(s) 2020
Abstract Our study is aimed at assessing the extent at which relying on differing geostatistical approaches may affect characterization of the connectivity of geomaterials (or facies) and, in turn, model calibration outputs in highly heterogeneous aquifers. We set our study within a probabilistic framework, by relying on a numerical Monte Carlo (MC) approach. The reconstruction of the spatial distribution of geomaterials and flow simulations are patterned after a field scenario corresponding to the aquifer system serving the city of Bologna (Northern Italy). Two collections of MC realizations of facies distributions, conditional on available lithological data, are generated through two alternative geostatistically-based techniques, i.e., Sequential Indicator and Transition-Probability simulation. Hydraulic conductivity values of the least- and most-conductive facies are estimated within each MC simulation in the context of a Maximum Likelihood (ML) approach by considering available piezometric data. We provide evidence that the choice of the facies reconstruction technique (1) impacts the degree of connectivity of facies whose proportions are close to the percolation threshold while (2) is not sensibly affecting the connectivity associated with facies whose proportions are much larger than the percolation threshold. By relying on the unique (lithological and hydrological) data-set at our disposal, we also explore the performance of MLbased model identification criteria to (1) discriminate amongst competitive facies reconstruction geostatistical models and (2) quantify the (posterior probabilistic) weight associated with each model. We then show that ML-based model averaging provides estimates of hydraulic heads which are slightly more in agreement with available data when compared to the bestperforming realization in the T-PROGS set than considering its counterpart associated with the SISIM-based collection. Keywords Geostatistical reconstruction Connectivity metrics Groundwater flow model calibration Bayesian model averaging Transition probability
1 Introduction Probabilistic reconstruction of complex aquifer systems is nowadays considered a common practice to account for uncertainty resulting from our lack of knowledge of highlyheterogeneous subsurface structures and spatial distribution of attributes of aquifer systems. In a stochastic framework, aquifer heterogeneity can be conceptualized by considering system attributes, such as hydraulic conductivity, as random functions. Widely employed geostatistical methods (e.g., Sequential Gaussian Simulation; Deutsch and Journel & Martina Siena [email protected] 1
Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano, Piazza L. Da Vinci 32, 20133 Milan, Italy
1998) describe hydraulic properties as multivariate Gaussian random fields and quantify the degree
Data Loading...
