Résultat de la recherche
Can one identify karst conduit networks geometry and properties from hydraulic and tracer test data?
Karst aquifers are characterized by extreme heterogeneity due to the presence of karst conduits embedded in a fractured matrix having a much lower hydraulic conductivity. The resulting contrast in the physical properties of the system implies that the system reacts very rapidly to some changes in the boundary conditions and that numerical models are extremely sensitive to small modifications in properties or positions of the conduits. Furthermore, one major issue in all those models is that the location and size of the conduits is generally unknown. For all those reasons, estimating karst network geometry and their properties by solving an inverse problem is a particularly difficult problem.
In this paper, two numerical experiments are described. In the first one, 18,000 flow and transport simulations have been computed and used in a systematic manner to assess statistically if one can retrieve the parameters of a model (geometry and radius of the conduits, hydraulic conductivity of the conduits) from head and tracer data. When two tracer test data sets are available, the solution of the inverse problems indicate with high certainty that there are indeed two conduits and not more. The radius of the conduits are usually well identified but not the properties of the matrix. If more conduits are present in the system, but only two tracer test data sets are available, the inverse problem is still able to identify the true solution as the most probable but it also indicates that the data are insufficient to conclude with high certainty.
In the second experiment, a more complex model (including non linear flow equations in conduits) is considered. In this example, gradient-based optimization techniques are proved to be efficient for estimating the radius of the conduits and the hydraulic conductivity of the matrix in a promising and efficient manner.
These results suggest that, despite the numerical difficulties, inverse methods should be used to constrain numerical models of karstic systems using flow and transport data. They also suggest that a pragmatic approach for these complex systems could be to generate a large set of karst conduit network realizations using a pseudo-genetic approach such as SKS, and for each karst realization, flow and transport parameters could be optimized using a gradient-based search such as the one implemented in PEST.
Stochastic forecasts of seawater intrusion in the Korba aquifer (Tunisia) towards a sustainable groundwater management
Truncated plurigaussian simulations to characterize aquifer heterogeneity
Evaluation of reactive transport parameters to assess specific vulnerability in karst systems
Modelling of the predictive hydrogeological impacts of the radwaste geological repository construction on limestone aquifers of the Meuse/Haute-Marne site (France)
Analysis and evaluation of hydraulic impacts of the shafts and the tunnel access to radwaste repository on the Meuse/Haute-Marne site by the use of finite element variably saturated flow model
The influence of paleoclimatic events on the functioning of an alpine thermal system (France): the contribution of hydrodynamic–thermal modeling
Numerical models of the Aix-les-Bains thermal aquifer (France) were used to investigate the influence of Quaternary paleoclimatic events on the current thermal state of the groundwater. Initial numerical tests were successful in that present-day fluid flows (heads and flow rates) and the resulting velocities were compatible with residence time data. Water flowing through an aquifer cools the rock mass; therefore, the rate of water flow governs the outlets temperature. For the Aix-les-Bains aquifer, applying present-day flow rates to the entire history of the aquifer leads to much more substantial cooling of the rock mass than is indicated by the outlets temperature (i.e. present-day flow rates are 10 times too high). This suggests that the aquifer may have gone through alternating functioning phases, during which the rock mass cooled, and blocked phases, during which the aquifer reheated. Other results indicate that the main parameters affecting thermal behavior during a functioning phase are the total inflow volume, rather than individual inflow rates, and the initial heat field. As phenomena linked to glaciation can lead to the blocking of infiltration zones and aquifer outlets, the findings suggest that the hypothesis of intermittent aquifer functioning related to glaciations is compatible with the current thermal field.
A method for the stochastic modeling of karstic systems accounting for geophysical data: an example of application in the region of Tulum, Yucatan Peninsula (Mexico)
Analytical modelling of stable isotope fractionation of volatile organic compounds in the unsaturated zone
Analytical models were developed that simulate stable isotope ratios of volatile organic compounds (VOCs) near a point source contamination in the unsaturated zone. The models describe diffusive transport of VOCs, biodegradation and source ageing. The mass transport is governed by Fick's law for diffusion. The equation for reactive transport of VOCs in the soil gas phase was solved for different source geometries and for different boundary conditions. Model results were compared to experimental data from a one-dimensional laboratory column and a radial-symmetric field experiment. The comparison yielded a satisfying agreement. The model results clearly illustrate the significant isotope fractionation by gas phase diffusion under transient state conditions. This leads to an initial depletion of heavy isotopes with increasing distance from the source. The isotope evolution of the source is governed by the combined effects of isotope fractionation due to vaporisation, diffusion and biodegradation. The net effect can lead to an enrichment or depletion of the heavy isotope in the remaining organic phase, depending on the compound and element considered. Finally, the isotope evolution of molecules migrating away from the source and undergoing degradation is governed by a combined degradation and diffusion isotope effect. This suggests that, in the unsaturated zone, the interpretation of biodegradation of VOC based on isotopic data must always be based on a model combining gas phase diffusion and degradation.