HR: 17:30h
AN: H14C-06    [Abstracts]
TI: An Inverse Model of Three-Dimensional Flow and Transport in Heterogeneous Porous Media
AU: * Robinson, B A
EM: robinson@lanl.gov
AF: Los Alamos National Laboratory, PO Box 1663 LANL, Los Alamos, NM 87545, United States
AU: Vrugt, J A
EM: vrugt@lanl.gov
AF: Los Alamos National Laboratory, PO Box 1663 LANL, Los Alamos, NM 87545, United States
AU: Yoon, H
EM: hyoon3@uiuc.edu
AF: University of Illinois at Urbana Champaign, 205 N. Mathew Ave MC-250, Urbana, IL 61801, United States
AU: Zhang, C
EM: czhang@uiuc.edu
AF: University of Illinois at Urbana Champaign, 205 N. Mathew Ave MC-250, Urbana, IL 61801, United States
AU: Werth, C J
EM: werth@uiuc.edu
AF: University of Illinois at Urbana Champaign, 205 N. Mathew Ave MC-250, Urbana, IL 61801, United States
AU: Kitanidis, P K
EM: peterk@stanford.edu
AF: Stanford University, 380 Panama Mall Terman Building, Room M-19, Stanford, CA 94305, United States
AU: Lichtner, P C
EM: lichtner@lanl.gov
AF: Los Alamos National Laboratory, PO Box 1663 LANL, Los Alamos, NM 87545, United States
AU: Lu, C
EM: clu@lanl.gov
AF: Los Alamos National Laboratory, PO Box 1663 LANL, Los Alamos, NM 87545, United States
AB: A three-dimensional flow and transport model was developed to simulate the results of a laboratory-scale experiment in which snapshots of concentration were obtained using magnetic resonance imaging (MRI) during the displacement of tracer through a 14 by 8 by 8 cm flow cell. The medium was deliberately constructed to be heterogeneous with a known spatial correlation structure using sand of five different grain-size distributions. The extremely well characterized flow cell and large, high-precision data set of concentrations during displacement make this a unique experiment for examining the validity of flow and transport models, and for exploring new methods for interpreting large data sets using advanced optimization algorithms. A transport model was constructed by solving the steady state flow equations using the Finite Element Heat and Mass (FEHM) code, using FEHM's particle tracking transport model for simulating tracer migration. The particle tracking model was selected so that precise estimates of the transport parameters could be obtained that are not corrupted by numerical dispersion; a large number of particles (typically one million) were required to provide accuracy. The inverse model included nine uncertain parameters, the five permeability values of the individual sand units, and four dispersion/diffusion parameters. The inverse problem was solved with AMALGAM and DREAM, two recently developed self-adaptive multimethod optimization algorithms. The computations were enabled by performing both the transport model and the optimization loop on a high-performance computing cluster. Computational results indicate that parameter estimates and increased understanding of the behavior of the system can be obtained, and significant improvements in the fit to the data over hand calibration can be achieved, using this inverse modeling approach. The study also illustrates that numerical methods that make effective use of high- performance computing resources and advanced optimization algorithms are crucial in enabling the interpretation of very large data sets using large-scale, distributed parameter models.
DE: 0550 Model verification and validation
DE: 1805 Computational hydrology
DE: 1832 Groundwater transport
DE: 1846 Model calibration (3333)
SC: Hydrology [H]
MN: 2007 Fall Meeting