Sparse data integration for the interpolation of concentration measurements using kriging in natural coordinates

Summary Groundwater contaminant plumes often display a curvilinear anisotropy which conventional kriging and geostatistical simulation approaches generally fail to reproduce. In this paper, we use a physically relevant coordinate transformation in order to improve the kriging of contaminant plumes in 2D. The proposed coordinate transformation maps the Cartesian grid into the natural coordinates of flow: the hydraulic head and the stream function. This simplifies the specification of a nonlinear anisotropy by modifying the spatial relationships between grid points. The computation of the natural coordinates, however, requires the availability of a flow model, which often needs to be defined using limited data. For this reason, a data integration procedure is included in the methodology to support the use of the natural coordinate transformation in real cases. The performance of the approach is tested on two synthetic test cases where it produces concentration maps that are more accurate than those obtained with Cartesian coordinates kriging. These test cases also highlight some limitations of the approach. Whereas the transformation enables kriging to account for advection, results show the need to consider dispersion as well. Further work is also required to generalize the approach to 3D cases and to groundwater flow including wells and sources.

[1]  Giansilvio Ponzini,et al.  Identification of aquifer transmissivities: The comparison model method , 1982 .

[2]  Philippe Lagacherie,et al.  Modelling Spatial Variability Along Drainage Networks with Geostatistics , 2006 .

[3]  J. Hoef,et al.  Spatial statistical models that use flow and stream distance , 2006, Environmental and Ecological Statistics.

[4]  The Principal Direction Technique: A New Approach to Groundwater Contaminant Transport Modeling , 1982 .

[5]  P Goovaerts,et al.  Geostatistical assessment and validation of uncertainty for three-dimensional dioxin data from sediments in an estuarine river. , 2001, Environmental science & technology.

[6]  Phaedon C. Kyriakidis,et al.  Forward and Inverse Transformations between Cartesian and Channel-fitted Coordinate Systems for Meandering Rivers , 2007 .

[7]  Olaf A Cirpka,et al.  Interpolation of Steady‐State Concentration Data by Inverse Modeling , 2009, Ground water.

[8]  A. Michalak,et al.  A Geostatistical Data Assimilation Approach for Estimating Groundwater Plume Distributions from Multiple Monitoring Events , 2013 .

[9]  Graham E. Fogg,et al.  Automatic Generation of Flow Nets with Conventional Ground-Water Modeling Algorithmsa , 1985 .

[10]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[11]  Clayton V. Deutsch,et al.  Programs for kriging and sequential Gaussian simulation with locally varying anisotropy using non-Euclidean distances , 2011, Comput. Geosci..

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

[13]  A. Michalak,et al.  A geostatistical framework for incorporating transport information in estimating the distribution of a groundwater contaminant plume , 2005 .

[14]  Philippe Pasquier,et al.  Hydraulic head field estimation using kriging with an external drift: A way to consider conceptual model information , 2008 .

[15]  J. Snepvangers,et al.  Mapping Curvilinear Structures with Local Anisotropy Kriging , 2005 .

[16]  M. Dagbert,et al.  Computing Variograms in Folded Strata-Controlled Deposits , 1984 .

[17]  Clayton V. Deutsch,et al.  Kriging in the Presence of Locally Varying Anisotropy Using Non-Euclidean Distances , 2009 .

[18]  Anders Løland,et al.  Spatial covariance modelling in a complex coastal domain by multidimensional scaling , 2003 .

[19]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[20]  C. Fouquet,et al.  Modèles géostatistiques de concentrations ou de débits le long des cours d'eau , 2006 .

[21]  A. Dassargues,et al.  Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium) , 2009 .

[22]  Philippe Pasquier,et al.  Steady- and transient-state inversion in hydrogeology by successive flux estimation , 2006 .

[23]  P. Guttorp,et al.  Nonparametric Estimation of Nonstationary Spatial Covariance Structure , 1992 .

[24]  Numerical generation of flow nets - The FLOWNS model. , 1990 .

[25]  Frank C. Curriero,et al.  On the Use of Non-Euclidean Distance Measures in Geostatistics , 2007 .

[26]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[27]  G. Matanga,et al.  Stream functions in three‐dimensional groundwater flow , 1993 .

[28]  B. Minsker,et al.  Spatial Interpolation Methods for Nonstationary Plume Data , 2004, Ground water.

[29]  L. Feyen,et al.  Quantifying geological uncertainty for flow and transport modeling in multi-modal heterogeneous formations , 2006 .

[30]  Stephen L. Rathbun,et al.  Spatial modelling in irregularly shaped regions: Kriging estuaries , 1998 .

[31]  H. T. Buxton,et al.  Patterns and Rates of Ground‐Water Flow on Long Island, New York , 1992 .

[32]  Clayton V. Deutsch,et al.  Geostatistical Reservoir Modeling , 2002 .

[33]  Stan Matwin,et al.  Machine Learning for the Detection of Oil Spills in Satellite Radar Images , 1998, Machine Learning.

[34]  Wenlong Xu,et al.  Conditional curvilinear stochastic simulation using pixel-based algorithms , 1996 .

[35]  U. Meissner,et al.  Finite elements in water resources. , 1982 .

[36]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[37]  Rainer Helmig,et al.  Streamline-oriented grid generation for transport modelling in two-dimensional domains including wells , 1999 .

[38]  Clayton V. Deutsch,et al.  Hierarchical object-based stochastic modeling of fluvial reservoirs , 1996 .

[39]  Dionissios T. Hristopulos,et al.  On the physical geometry concept at the basis of space/time geostatistical hydrology , 2000 .

[40]  Emil O. Frind,et al.  The Dual Formulation of Flow for Contaminant Transport Modeling: 1. Review of Theory and Accuracy Aspects , 1985 .

[41]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.