Solving iTOUGH2 simulation and optimization problems using the PEST protocol

被引:59
作者
Finsterle, Stefan [1 ]
Zhang, Yingqi [1 ]
机构
[1] Lawrence Berkeley Natl Lab, Div Earth Sci, Berkeley, CA 94720 USA
关键词
Optimization; Sensitivity analysis; Inverse modeling; Uncertainty quantification; iTOUGH2; PEST; UNSATURATED FLOW; AQUIFER PARAMETERS; MODEL; SEEPAGE; INVERSE; FIELD; ALGORITHM; TRANSIENT; SITE;
D O I
10.1016/j.envsoft.2011.02.008
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The PEST protocol has been implemented into the iTOUGH2 code, allowing the user to link any simulation program (with ASCII-based inputs and outputs) to iTOUGH2's sensitivity analysis, inverse modeling, and uncertainty quantification capabilities. These application models can be pre- or post-processors of the TOUGH2 non-isothermal multiphase flow and transport simulator, or programs that are unrelated to the TOUGH suite of codes. PEST-style template and instruction files are used, respectively, to pass input parameters updated by the iTOUGH2 optimization routines to the model, and to retrieve the model-calculated values that correspond to observable variables. We summarize the iTOUGH2 capabilities and demonstrate the flexibility added by the PEST protocol for the solution of a variety of simulation optimization problems. In particular, the combination of loosely coupled and tightly integrated simulation and optimization routines provides both the flexibility and control needed to solve challenging inversion problems for the analysis of multiphase subsurface flow and transport systems. (C) 2011 Elsevier Ltd. All rights reserved.
引用
收藏
页码:959 / 968
页数:10
相关论文
共 51 条
[31]  
Kowalsky MB, 2005, WATER RESOUR RES, V41, DOI [10.1029/2005wr004237, 10.1029/2005WR004237]
[32]   Estimating flow parameter distributions using ground-penetrating radar and hydrological measurements during transient flow in the vadose zone [J].
Kowalsky, MB ;
Finsterle, S ;
Rubin, Y .
ADVANCES IN WATER RESOURCES, 2004, 27 (06) :583-599
[33]  
Levenberg K, 1944, Q Appl Math, V2, P164, DOI [10.1090/QAM/10666, 10.1090/qam/10666, DOI 10.1090/QAM/10666, DOI 10.1090/QAM/1944-02-02]
[34]   Inversion of tracer test data using tomographic constraints [J].
Linde, N ;
Finsterle, S ;
Hubbard, S .
WATER RESOURCES RESEARCH, 2006, 42 (04)
[35]   A note on unsaturated flow in two-dimensional fracture networks [J].
Liu, HH ;
Bodvarsson, GS ;
Finsterle, S .
WATER RESOURCES RESEARCH, 2002, 38 (09)
[36]   AN ALGORITHM FOR LEAST-SQUARES ESTIMATION OF NONLINEAR PARAMETERS [J].
MARQUARDT, DW .
JOURNAL OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS, 1963, 11 (02) :431-441
[37]   EQUATION OF STATE CALCULATIONS BY FAST COMPUTING MACHINES [J].
METROPOLIS, N ;
ROSENBLUTH, AW ;
ROSENBLUTH, MN ;
TELLER, AH ;
TELLER, E .
JOURNAL OF CHEMICAL PHYSICS, 1953, 21 (06) :1087-1092
[38]   The coprecipitation of Sr into calcite precipitates induced by bacterial ureolysis in artificial groundwater: Temperature and kinetic dependence [J].
Mitchell, AC ;
Ferris, FG .
GEOCHIMICA ET COSMOCHIMICA ACTA, 2005, 69 (17) :4199-4210
[39]   Evaluation of alternative designs for an injectable subsurface barrier at the Brookhaven National Laboratory site, Long Island, New York [J].
Moridis, GJ ;
Finsterle, S ;
Heiser, J .
WATER RESOURCES RESEARCH, 1999, 35 (10) :2937-2953
[40]   A SIMPLEX-METHOD FOR FUNCTION MINIMIZATION [J].
NELDER, JA ;
MEAD, R .
COMPUTER JOURNAL, 1965, 7 (04) :308-313