Low-dimensional modelling of numerical groundwater flow

Numerical models are often used for simulating groundwater flow. Written in state-space form, the dimension of these models is of the order of the number of grid cells used and can be very high (more than a million). As a result, these models are computationally very demanding, especially if many different scenarios have to be simulated. In this paper we introduce a model reduction approach to develop an approximate model with a significantly reduced dimension. The reduction method is based upon several simulations of the large-scale numerical model. By computing the covariance matrix of the model results, we obtain insight into the variability of the model behaviour. Moreover, by selecting the leading eigenvectors of this covariance matrix, we obtain the spatial patterns that represent the directions in state space where the model variability is dominant. These patterns are also called empirical orthogonal functions. We can project the original numerical model onto those dominant spatial patterns. The result is a low-dimensional model that is still able to reproduce the dominant model behaviour. Copyright © 2004 John Wiley & Sons, Ltd.

[1]  Andrew J. Newman,et al.  Model Reduction via the Karhunen-Loeve Expansion Part I: An Exposition , 1996 .

[2]  Ioannis G. Kevrekidis,et al.  Alternative approaches to the Karhunen-Loève decomposition for model reduction and data analysis , 1996 .

[3]  M. A. Hooimeijer.  Reduction of Complex Computational Models , 2001 .

[4]  Andrew J. Newman,et al.  Model Reduction via the Karhunen-Loeve Expansion: An Exposition , 1996 .

[5]  A. Segers,et al.  Variance reduced ensemble Kalman filtering , 2001 .

[6]  Gene H. Golub,et al.  Matrix computations , 1983 .

[7]  Andrés Sahuquillo,et al.  An eigenvalue numerical technique for solving unsteady linear groundwater models continuously in time , 1983 .

[8]  Roger Ghanem,et al.  Scales of fluctuation and the propagation of uncertainty in random porous media , 1998 .

[9]  Allan D. Woodbury,et al.  The Arnoldi reduction technique for efficient direct solution of radionuclide decay chain transport in dual-porosity media , 2000 .

[10]  Arthur Veldman,et al.  Proper orthogonal decomposition and low-dimensional models for driven cavity flows , 1998 .

[11]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[12]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[13]  Raymond A. Adomaitis.  RTCVD Model Reduction: A Collocation on Empirical Eigenfunctions Approach , 1995 .

[14]  H. Park,et al.  On the solution of inverse heat transfer problem using the Karhunen–Loève Galerkin method , 1999 .

[15]  H. M. Park,et al.  Low dimensional modeling of flow reactors , 1996 .

[16]  Allan D. Woodbury,et al.  Application of the Lanczos algorithm to the simulation of groundwater flow in dual-porosity media , 2000 .

[17]  Allan D. Woodbury,et al.  Application of the Lanczos Algorithm to the solution of the groundwater flow equation , 1989 .

[18]  Ian T. Jolliffe,et al.  Applied Factor Analysis in the Natural Sciences. , 1994 .