scholarly journals Comparison of Two Ensemble Kalman-Based Methods for Estimating Aquifer Parameters from Virtual 2-D Hydraulic and Tracer Tomographic Tests

Geosciences ◽  
2020 ◽  
Vol 10 (7) ◽  
pp. 276 ◽  
Author(s):  
Emilio Sánchez-León ◽  
Daniel Erdal ◽  
Carsten Leven ◽  
Olaf A. Cirpka

We compare two ensemble Kalman-based methods to estimate the hydraulic conductivity field of an aquifer from data of hydraulic and tracer tomographic experiments: (i) the Ensemble Kalman Filter (EnKF) and (ii) the Kalman Ensemble Generator (KEG). We generated synthetic drawdown and tracer data by simulating two pumping tests, each followed by a tracer test. Parameter updating with the EnKF is performed using the full transient signal. For hydraulic data, we use the standard update scheme of the EnKF with damping, whereas for concentration data, we apply a restart scheme, in which solute transport is resimulated from time zero to the next measurement time after each parameter update. In the KEG, we iteratively assimilate all observations simultaneously, here inverting steady-state heads and mean tracer arrival times. The inversion with the dampened EnKF worked well for the transient pumping-tests, but less for the tracer tests. The KEG produced similar estimates of hydraulic conductivity but at significantly lower costs. We conclude that parameter estimation in well-defined hydraulic tests can be done very efficiently by iterative ensemble Kalman methods, and ambiguity between state and parameter updates can be completely avoided by assimilating temporal moments of concentration data rather than the time series themselves.

Geosciences ◽  
2020 ◽  
Vol 10 (11) ◽  
pp. 462
Author(s):  
Emilio Sánchez-León ◽  
Carsten Leven ◽  
Daniel Erdal ◽  
Olaf A. Cirpka

Pumping and tracer tests are site-investigation techniques frequently used to determine hydraulic conductivity. Tomographic test layouts, in which multiple tests with different combinations of injection and observation wells are performed, gain a better insight into spatial variability. While hydraulic tomography has repeatedly been applied in the field, tracer tomography lags behind. In a previous publication, we presented a synthetic study to investigate whether the ensemble Kalman Filter (EnKF) or the Kalman Ensemble Generator (KEG) performs better in inverting hydraulic- and tracer-tomographic data. In this work, we develop an experimental method for solute-tracer tomography using fluorescein as a conservative tracer. We performed hydraulic- and tracer-tomographic tests at the Lauswiesen site in Germany. We analyzed transient drawdown and concentration data with the EnKF and steady-state hydraulic heads and mean tracer arrival times with the KEG, obtaining more stable results with the KEG at lower computational costs. The spatial distribution of the estimated hydraulic conductivity field agreed with earlier descriptions of the aquifer at the site. This study narrows the gap between numerical studies and field applications for aquifer characterization at high resolution, showing the potential of combining ensemble-Kalman filter based methods with data collected from hydraulic and solute-tracer tomographic experiments.


2012 ◽  
Vol 9 (11) ◽  
pp. 13083-13115
Author(s):  
E. Crestani ◽  
M. Camporese ◽  
D. Baú ◽  
P. Salandin

Abstract. The significance of estimating the spatial variability of the hydraulic conductivity K in natural aquifers is relevant to the possibility of defining the space and time evolution of a non-reactive plume, since the transport of a solute is mainly controlled by the heterogeneity of K. At the local scale, the spatial distribution of K can be inferred by combining the Lagrangian formulation of the transport with a Kalman filter-based technique and assimilating a sequence of time-lapse concentration C measurements, which, for example, can be evaluated on-site through the application of a geophysical method. The objective of this work is to compare the ensemble Kalman filter (EnKF) and the ensemble smoother (ES) capabilities to retrieve the hydraulic conductivity spatial distribution in a groundwater flow and transport modeling framework. The application refers to a two-dimensional synthetic aquifer in which a tracer test is simulated. Moreover, since Kalman filter-based methods are optimal only if each of the involved variables fit to a Gaussian probability density function (pdf) and since this condition may not be met by some of the flow and transport state variables, issues related to the non-Gaussianity of the variables are analyzed and different transformation of the pdfs are considered in order to evaluate their influence on the performance of the methods. The results show that the EnKF reproduces with good accuracy the hydraulic conductivity field, outperforming the ES regardless of the pdf of the concentrations.


SPE Journal ◽  
2018 ◽  
Vol 23 (02) ◽  
pp. 449-466 ◽  
Author(s):  
Siavash Hakim Elahi ◽  
Behnam Jafarpour

Summary Hydraulic fracturing is performed to enable production from low-permeability and organic-rich shale-oil/gas reservoirs by stimulating the rock to increase its permeability. Characterization and imaging of hydraulically induced fractures is critical for accurate prediction of production and of the stimulated reservoir volume (SRV). Recorded tracer concentrations during flowback and historical production data can reveal important information about fracture and matrix properties, including fracture geometry, hydraulic conductivity, and natural-fracture density. However, the complexity and uncertainty in fracture and reservoir descriptions, coupled with data limitations, complicate the estimation of these properties. In this paper, tracer-test and production data are used for dynamic characterization of important parameters of hydraulically fractured reservoirs, including matrix permeability and porosity, planar-fracture half-length and hydraulic conductivity, discrete-fracture-network (DFN) density and conductivity, and fracture-closing (conductivity-decline) rate during production. The ensemble Kalman filter (EnKF) is used to update uncertain model parameters by sequentially assimilating first the tracer-test data and then the production data. The results indicate that the tracer-test and production data have complementary information for estimating fracture half-length and conductivity, with the former being more sensitive to hydraulic conductivity and the latter being more affected by fracture half-length. For characterization of DFN, a stochastic representation is adopted and the parameters of the stochastic model along with matrix and hydraulic-fracture properties are updated. Numerical examples are presented to investigate the sensitivity of the observed production and tracer-test data to fracture and matrix properties and to evaluate the EnKF performance in estimating these parameters.


2013 ◽  
Vol 17 (4) ◽  
pp. 1517-1531 ◽  
Author(s):  
E. Crestani ◽  
M. Camporese ◽  
D. Baú ◽  
P. Salandin

Abstract. Estimating the spatial variability of hydraulic conductivity K in natural aquifers is important for predicting the transport of dissolved compounds. Especially in the nonreactive case, the plume evolution is mainly controlled by the heterogeneity of K. At the local scale, the spatial distribution of K can be inferred by combining the Lagrangian formulation of the transport with a Kalman-filter-based technique and assimilating a sequence of time-lapse concentration C measurements, which, for example, can be evaluated on site through the application of a geophysical method. The objective of this work is to compare the ensemble Kalman filter (EnKF) and the ensemble smoother (ES) capabilities to retrieve the hydraulic conductivity spatial distribution in a groundwater flow and transport modeling framework. The application refers to a two-dimensional synthetic aquifer in which a tracer test is simulated. Moreover, since Kalman-filter-based methods are optimal only if each of the involved variables fit to a Gaussian probability density function (pdf) and since this condition may not be met by some of the flow and transport state variables, issues related to the non-Gaussianity of the variables are analyzed and different transformation of the pdfs are considered in order to evaluate their influence on the performance of the methods. The results show that the EnKF reproduces with good accuracy the hydraulic conductivity field, outperforming the ES regardless of the pdf of the concentrations.


2012 ◽  
Vol 428-429 ◽  
pp. 152-169 ◽  
Author(s):  
Liangping Li ◽  
Haiyan Zhou ◽  
J. Jaime Gómez-Hernández ◽  
Harrie-Jan Hendricks Franssen

1999 ◽  
Vol 30 (4-5) ◽  
pp. 333-360 ◽  
Author(s):  
Larry McKay ◽  
Johnny Fredericia ◽  
Melissa Lenczewski ◽  
Jørn Morthorst ◽  
Knud Erik S. Klint

A field experiment shows that rapid downward migration of solutes and microorganisms can occur in a fractured till. A solute tracer, chloride, and a bacteriophage tracer, PRD-1, were added to groundwater and allowed to infiltrate downwards over a 4 × 4 m area. Chloride was detected in horizontal filters at 2.0 m depth within 3-40 days of the start of the tracer test, and PRD-1 was detected in the same filters within 0.27 - 27 days. At 2.8 m depth chloride appeared in all the filters, but PRD-1 appeared in only about one-third of the filters. At 4.0 m depth chloride appeared in about one-third of the filters and trace amounts of PRD-1 were detected in only 2 of the 36 filters. Transport rates and peak tracer concentrations decreased with depth, but at each depth there was a high degree of variability. The transport data is generally consistent with expectations based on hydraulic conductivity measurements and on the observed density of fractures and biopores, both of which decrease with depth. Transport of chloride was apparently retarded by diffusion into the fine-grained matrix between fractures, but the rapid transport of PRD-1, with little dispersion, indicates that it was transported mainly through the fractures.


1985 ◽  
Vol 17 (2-3) ◽  
pp. 197-208 ◽  
Author(s):  
H. Bode ◽  
C. F. Seyfried

The interrelationship between mixing characteristics and tracer response curves in activated sludge tanks is explained. In some cases the return sludge cycle has a strong influence on the tracer response curves. Results from tracer tests in the field are hard to interpret because the tracer in the return sludge interferes with the initial tracer. Therefore a special evaluation procedure has to be applied. The paper closes with results from a field tracer test study.


2021 ◽  
Author(s):  
Arezou Dodangeh ◽  
Mohammad Mahdi Rajabi ◽  
Marwan Fahs

<p>In coastal aquifers, we face the problem of salt water intrusion, which creates a complex flow field. Many of these coastal aquifers are also exposed to contaminants from various sources. In addition, in many cases there is no information about the characteristics of the aquifer. Simultaneous identification of the contaminant source and coastal aquifer characteristics can be a challenging issue. Much work has been done to identify the contaminant source, but in the complex velocity field of coastal aquifer, no one has resolved this issue yet. We want to address that in a three-dimensional artificial coastal aquifer.</p><p>To achieve this goal, we have developed a method in which the contaminant source can be identified and the characteristics of the aquifer can be estimated by using information obtained from observation wells. First, by assuming the input parameters required to simulate the contaminant transfer to the aquifer, this three-dimensional coastal aquifer that is affected by various phenomena such as seawater intrusion, tides, shore slope, rain, discharge and injection wells, is simulated and the time series of the output parameters including head, salinity and contaminant concentration are estimated. In the next step, with the aim of performing inverse modeling, random values ​​are added to the time series of outputs obtained at specific points (points belonging to observation wells) in order to rebuilt the initial conditions of the problem to achieve the desired unknowns (contaminant source and aquifer characteristics). The unknowns estimated in this study are the contaminant source location (x, y, z), the initial contaminant concentration, the horizontal and vertical hydraulic conductivity of the aquifer. SEAWAT model in GMS software environment has been used to solve the equations of flow and contaminant transfer and simulate a three-dimensional coastal aquifer. Next, for reverse modeling, one of the Bayesian Filters subset (ensemble Kalman filter) has been used in the Python programming language environment. Also, to reduce the code run time, the neural network model is designed and trained for the SEAWAT model.</p><p>This method is able to meet the main purpose of the study, namely estimating the value ​​of unknown input parameters, including the contaminant source location, the initial contaminant concentration, the horizontal and vertical hydraulic conductivity of the aquifer. In addition, that makes it possible to achieve a three-dimensional numerical model of the coastal aquifer that can be used as a benchmark to examine more accurately the impact of different phenomena simultaneously. In conclusion, we have developed an algorithm which can be used in the world's coastal aquifers to identify the contaminant source and estimate its characteristics.</p><p> </p><p>Key words: coastal aquifer, seawater intrusion, contaminants, groundwater, flow field, parameter estimation, ensemble kalman filter</p>


Sign in / Sign up

Export Citation Format

Share Document