z-logo
Premium
Inverse modeling for locating dense nonaqueous pools in groundwater under steady flow conditions
Author(s) -
Sciortino Antonella,
Harmon Thomas C.,
Yeh William WG.
Publication year - 2000
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2000wr900047
Subject(s) - inverse problem , inverse , minification , mathematics , least squares function approximation , flow (mathematics) , sensitivity (control systems) , constant (computer programming) , mathematical optimization , statistics , computer science , mathematical analysis , geometry , engineering , estimator , electronic engineering , programming language
In this work we develop an inverse modeling procedure to identify the location and the dimensions of a single‐component dense nonaqueous phase liquid (DNAPL) pool in a saturated porous medium under steady flow conditions. The inverse problem is formulated as a least squares minimization problem and solved by a search procedure based on the Levenberg‐Marquardt method. Model output is calculated by an existing three‐dimensional analytical model describing the transport of solute from a dissolving distributed noise upon the forward model‐generated concentration field. We further test the algorithm's ability to predict the location and size of a DNAPL pool placed in a controlled three‐dimensional bench‐scale experiment. In this case we apply the Levenberg‐Marquardt algorithm to the minimization of three types of residuals: ordinary residuals, weighted residuals with weights equal to the square of the inverse of the observations, and weighted residuals with weights obtained by adding a constant term to the observed concentrations. The results are sensitive to the location of the observation wells and to the type of residuals minimized. In general, better results in terms of pool location and dimensions were obtained by the minimization of weighted residuals with weights obtained by adding a constant term to the observed concentrations. The results also indicate that the inverse problem is nonunique and nonconvex even in the absence of observation errors. Finally, the sensitivity of the inverse modeling scheme to transport parameter uncertainty was addressed. The inverse solution was found to be extremely sensitive to errors in the dispersion coefficients and relatively insensitive to errors in the mass transfer coefficient.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here