References of "Nguyen, Frédéric"
     in
Bookmark and Share    
Full Text
Peer Reviewed
See detailQuantitative temperature monitoring of a heat tracing experiment using cross-borehole ERT
Hermans, Thomas ULg; Wildemeersch, Samuel ULg; Jamin, Pierre ULg et al

in Geothermics (in press)

The growing demand for renewable energy leads to an increase in the development of geothermal energy projects and heat has become a common tracer in hydrology and hydrogeology. Designing geothermal ... [more ▼]

The growing demand for renewable energy leads to an increase in the development of geothermal energy projects and heat has become a common tracer in hydrology and hydrogeology. Designing geothermal systems requires a multidisciplinary approach including geological and hydrogeological aspects. In this context, electrical resistivity tomography (ERT) can bring relevant, qualitative and quantitative information on the temperature distribution in operating shallow geothermal systems or during heat tracing experiments. We followed a heat tracing experiment in an alluvial aquifer using cross-borehole time-lapse ERT. Heated water was injected in a well while water of the aquifer was extracted at another well. An ERT section was set up across the main flow direction. The results of ERT were transformed into temperature using calibrated petrophysical relationships. These ERT-derived temperatures were then compared to direct temperature measurements in control piezometers collected with distributed temperature sensing (DTS) and groundwater temperature loggers. Spatially, it enabled to map the horizontal and vertical extent of the heated water plume, as well as the zones where maximum temperatures occurred. Quantitatively, the temperatures and breakthrough curves estimated from ERT were in good agreement with the ones observed directly during the rise and maximum of the curve. An overestimation, likely related to 3D effects, was observed for the tail of the heat breakthrough curve. The error made on temperature can be estimated to be between 10 to 20 %, which is a fair value for indirect measurements. From our data, we estimated a quantification threshold for temperature variation of 1.2°C. These results suggest that ERT should be considered when designing heat tracing experiments or geothermal systems. It could help also to assess the geometrical complexity of the concerned reservoirs. It also appears that ERT could be a useful tool to monitor and control geothermal systems once they are in operation. [less ▲]

Detailed reference viewed: 13 (7 ULg)
Full Text
Peer Reviewed
See detailCase studies of incorporation of prior information in electrical resistivity tomography: comparison of different approaches
Hermans, Thomas ULg; Caterina, David ULg; Nguyen, Frédéric ULg

in Near Surface Geophysics (in press)

Many geophysical inverse problems are ill-posed and their solution non-unique. It is thus important to reduce the amount of mathematical solutions to more geologically plausible models by regularizing the ... [more ▼]

Many geophysical inverse problems are ill-posed and their solution non-unique. It is thus important to reduce the amount of mathematical solutions to more geologically plausible models by regularizing the inverse problem and incorporating all available prior information in the inversion process. We compare three different ways to incorporate prior information for electrical resistivity tomography (ERT): using a simple reference model or adding structural constraints to Occam's inversion and using geostatistical constraints. We made the comparison on four real cases representing different field applications in terms of scales of investigation and level of heterogeneities. In those cases, when electromagnetic logging data are available in boreholes to control the solution, it appears that incorporating prior information clearly improves the correspondence with logging data compared to the standard smoothness constrain. However, the way to incorporate it may have a major impact on the solution. A reference model can often be used to constrain the inversion; however, it can lead to misinterpretation if its weight is too strong or the resistivity values inappropriate. When the computation of the vertical and/or horizontal correlation length is possible, the geostatistical inversion gives reliable results everywhere in the section. However, adding geostatistical constraints can be difficult when there is not enough data to compute correlation lengths. When a known limit between two layers exists, the use of structural constrain seems to be more indicated particularly when the limit is located in zones of low sensitivity for ERT. This work should help interpreters to include their prior information directly into the inversion process through an appropriate way. [less ▲]

Detailed reference viewed: 33 (16 ULg)
Peer Reviewed
See detailInverting Hydraulic Heads In An Alluvial Aquifer Constrained With Electrical Resistivity Tomography Data Through Multiple-Point Statistics And Probability Perturbation Method: A Case Study
Hermans, Thomas ULg; Scheidt, Celine; Caers, Jef et al

Conference (2014, July)

Solving spatial inverse problems in the Earth Sciences remains a considerable challenge given the large number of parameters to invert for, the non-linearity of forward models and as a result the ill ... [more ▼]

Solving spatial inverse problems in the Earth Sciences remains a considerable challenge given the large number of parameters to invert for, the non-linearity of forward models and as a result the ill-posedness of the problem. Geostatistics is therefore needed to specify prior models, more particularly, information to control the spatial features of the inverse solutions. We used multiple-point statistics (MPS) to build models of pre-defined hydrofacies: clay, sand and gravel facies constrained to geological data (hard data) and geophysical data (soft data). The electrical resistivity tomography method was chosen to bring relevant spatially distributed information on the presence of the facies, given its sensitivity to variations in lithology and porosity. The comparison of the geophysical signature of the deposits with direct observations in boreholes enables to derive the conditional probability of observing a facies given its electrical resistivity. This is used to produce probability maps for each facies and constrain stochastic simulations of the alluvial aquifer. Then, the probability perturbation method (PPM) is used to integrate hydraulic heads data, using MPS to generate models. This process enables us to obtain calibrated models of the aquifer. The PPM algorithm will automatically seek solutions fitting both hydrogeological data and training-image based geostatistical constraints. Only geometrical features of the model are affected by the perturbation, i.e. we do not attempt to directly find the optimal value of hydrogeological parameters (chosen a priori), but the optimal spatial distribution of facies whose prior distribution is quantified in a training image. The methodology is first tested with a synthetic benchmark. The tests performed show that the choice of the training image is a major source of uncertainty. Therefore, one first needs to select those training images consistent with the geophysical data (and hence reject the inconsistent ones). Then, we proceed with them to hydrogeological inversions. Geophysical data (soft constraints) acts as an accelerator of convergence by reducing prior uncertainty. The hydraulic conductivity of each facies is a sensitive parameter, but it can be easily optimized prior to the PPM process. The stochastic method is then successfully applied within the context of an alluvial aquifer submitted to a pumping experiment. We show how the integration of various sources of data (borehole logs, geophysics, hydraulic heads) aids in calibrating hydrogeological models, locating high hydraulic conductivity zones and reducing uncertainty. The developed methodology proposes a common framework (multiple-point statistics) to integrate various information sources with variable resolutions relevant for hydrogeology: geological, geophysical and hydrogeological data. The method can be extended to integrate tracer tests to enable the calibration of transport parameters as well. The originality of the method is to use geophysical data both to refine the choice of the training image and to constrain the inversion of hydrogeological models. [less ▲]

Detailed reference viewed: 7 (1 ULg)
See detailUtilisation de la géostatistique multi-points pour l'intégration de données de tomographie de résistivité électrique aux modèles hydrogéologiques
Hermans, Thomas ULg; Scheidt, Céline; Caers, Jef et al

Scientific conference (2014, January 15)

Detailed reference viewed: 30 (13 ULg)
Full Text
Peer Reviewed
See detailDétermination des propriétés de couches dans le béton à l'aide d'un géoradar commercial à hautes fréquences: approche pic-à-pic et analyse fréquentielle du coefficient de réflexion
Van der Wielen, Audrey ULg; Nguyen, Frédéric ULg; Courard, Luc ULg

in Annales du Bâtiment et des Travaux Publics (2014)

The Ground Penetrating Radar (GPR) is an efficient tool for the non-destructive inspection of concrete structures. It is widely used for the detection of rebars or humid zones or for evaluating the ... [more ▼]

The Ground Penetrating Radar (GPR) is an efficient tool for the non-destructive inspection of concrete structures. It is widely used for the detection of rebars or humid zones or for evaluating the thickness of elements. But when an element contains a thin layer, the radar waves are submitted to multiple reflections on the interfaces and the layer appears in the radargram as a single reflection, whose detailed analysis can allow determining the thickness and the permittivity of the thin layer. Two approaches were considered in this paper. In the first one, the analysis is based on the peak-to-peak reflection amplitude. The second approach uses a frequency analysis of the reflection coefficient, whose amplitude and phase can then be calculated for several frequencies. With this method, the thickness and permittivity of the layer can in theory be simultaneously determined. Both methods were numerically validated through finite difference simulations and experimentally tested on concrete samples containing an air layer of variable thickness. We showed that the frequency analysis allowed to reach a higher precision in the parameters estimation for a limited additional computing cost. The method efficiency depends on the conditions and is optimal for layers with a high permittivity presenting a large contrast with the matrix. [less ▲]

Detailed reference viewed: 8 (1 ULg)
Full Text
Peer Reviewed
See detailA modified DOI-based method to statistically estimate the depth of investigation of dc resistivity surveys
Deceuster, J.; Etienne, A.; Robert, Tanguy et al

in Journal of Applied Geophysics (2014), 103

Several techniques are available to estimate the depth of investigation or to identify possible artifacts in dc resistivity surveys. Commonly, the depth of investigation (DOI) is mainly estimated by using ... [more ▼]

Several techniques are available to estimate the depth of investigation or to identify possible artifacts in dc resistivity surveys. Commonly, the depth of investigation (DOI) is mainly estimated by using an arbitrarily chosen cut-off value on a selected indicator (resolution, sensitivity or DOI index). Ranges of cut-off values are recommended in the literature for the different indicators. However, small changes in threshold values may induce strong variations in the estimated depths of investigation. To overcome this problem, we developed a new statistical method to estimate the DOI of dc resistivity surveys based on a modified DOI index approach. This method is composed of 5 successive steps. First, two inversions are performed by using different resistivity reference models for the inversion (0.1 and 10 times the arithmetic mean of the logarithm of the observed apparent resistivity values). Inversion models are extended to the edges of the survey line and to a depth range of three times the pseudodepth of investigation of the largest array spacing used. In step 2, we compute the histogram of a newly defined scaled DOI index. Step 3 consists of the fitting of the mixture of two Gaussian distributions (G1 and G2) to the cumulative distribution function of the scaled DOI index values. Based on this fitting, step 4 focuses on the computation of an interpretation index (II) defined for every cell j of the model as the relative probability density that the cell j belongs to G1, which describes the Gaussian distribution of the cells with a scaled DOI index close to 0.0. In step 5, a new inversion is performed by using a third resistivity reference model (the arithmetic mean of the logarithm of the observed apparent resistivity values). The final electrical resistivity image is produced by using II as alpha blending values allowing the visual discrimination between well-constrained areas and poorly-constrained cells. The efficiency of the proposed methodology is assessed on synthetic and field data. By using synthetic benchmark analysis, we demonstrate that the selected well-constrained cells are well-reconstructed in size and shape as well as in resistivity contrasts. Compared to the existing image appraisal tools, the proposed statistical method allows the identification of the statistically well-constrained cells of the model without using any arbitrary cut-off value. Using this statistical method in combination with the resolution, when interpreting dc resistivity surveys, provides the geophysicist valuable information to avoid over- or misinterpretation of ERT images. [less ▲]

Detailed reference viewed: 8 (2 ULg)
Peer Reviewed
See detailCoupling of hydrogeological models with hydrogeophysical data to characterize seawater intrusion and shallow geothermal systems
Beaujean, Jean ULg; Kemna, Andreas; Engesgaard, Peter et al

Conference (2013, December 12)

While coastal aquifers are being stressed due to climate changes and excessive groundwater withdrawals require characterizing efficiently seawater intrusion (SWI) dynamics, production of geothermal energy ... [more ▼]

While coastal aquifers are being stressed due to climate changes and excessive groundwater withdrawals require characterizing efficiently seawater intrusion (SWI) dynamics, production of geothermal energy is increasingly being used to hinder global warming. To study these issues, we need both robust measuring technologies and reliable predictions based on numerical models. SWI models are currently calibrated using borehole observations. Similarly, geothermal models depend mainly on the temperature field at few locations. Electrical resistivity tomography (ERT) can be used to improve these models given its high sensitivity to TDS and temperature and its relatively high lateral resolution. Inherent geophysical limitations, such as the resolution loss, can affect the overall quality of the ERT images and also prevent the correct recovery of the desired hydrochemical property. We present an uncoupled and coupled hydrogeophysical inversion to calibrate SWI and thermohydrogeologic models using ERT. In the SWI models, we demonstrate with two synthetic benchmarks (homogeneous and heterogeneous coastal aquifers) the ability of cumulative sensitivity-filtered ERT images using surface-only data to recover the hydraulic conductivity. Filtering of ERT-derived data at depth, where resolution is poorer, and the model errors make the dispersivity more difficult to estimate. In the coupled approach, we showed that parameter estimation is significantly improved because regularization bias is replaced by forward modeling only. Our efforts are currently focusing on applying the uncoupled/coupled approaches on a real life case study using field data from the site of Almeria, SE Spain. In the thermohydrogeologic models, the most sensitive hydrologic parameters responsible for heat transport are estimated from surface ERT-derived temperatures and ERT resistance data. A real life geothermal experiment that took place on the Campus De Sterre of Ghent University, Belgium and a synthetic case are tested. They consist in a thermal injection and storage of water in a shallow sandy aquifer. The use of a physically-based constraint accounting for the difference in conductivity between the formation and the tap injected water and based on the hydrogeological model calibrated first on temperatures is necessary to improve the parameter estimation. Results suggest that time-lapse ERT data may be limited but useful information for estimating groundwater flow and transport parameters for both the convection and conduction phases. [less ▲]

Detailed reference viewed: 31 (5 ULg)
Full Text
See detailMonitoring temperature changes during heat tracing experiments using electrical resistivity tomography
Hermans, Thomas ULg; Wildemeersch, Samuel ULg; Nguyen, Frédéric ULg

Conference (2013, December 06)

Thermal tracing experiments are becoming common in hydrogeology to estimate parameters governing heat transport processes and to study geothermal reservoirs. Electrical resistivity tomography (ERT) has ... [more ▼]

Thermal tracing experiments are becoming common in hydrogeology to estimate parameters governing heat transport processes and to study geothermal reservoirs. Electrical resistivity tomography (ERT) has proven its ability to monitor salt tracer tests, but few studies have investigated its performances, both qualitatively and quantitatively, in thermal tracing experiments. In this study, we monitored a heat injection and pumping experiment in an alluvial aquifer using both surface and crosshole ERT. The data sets of the surface profile, located along the main direction of flow, are distorted during injection by an electrical short-circuit through the external pumping-heating-injection experimental set-up. Current is flowing outside the subsurface leading to bad data for electrode dipoles located near the pumping and injection wells. The crosshole ERT panel is perpendicular to the main direction of flow. Difference inversion time-lapse images clearly show a preferential flow path in the bottom of the aquifer related to the presence of a coarse and clean gravel layer. Direct temperature measurements are available in control piezometers during the experiment to validate the ERT-derived temperatures and confirm the spatial pattern of temperature observed with ERT. Breakthrough curves are correctly retrieved in time and difference of 10 to 20% are observed for temperature estimation. The latter requires site-specific petrophysical laws and chemical stability assumptions that must be carefully verified. Our study proves that ERT, especially crosshole ERT, is a reliable tool to follow thermal tracing experiments but also to characterize heat transfer in the subsurface and to monitor geothermal resource exploitations. We also show that surface ERT may be impacted by the survey layout in unsuspected ways. [less ▲]

Detailed reference viewed: 15 (1 ULg)
Full Text
See detailReliability of resistivity-derived temperature: insights from laboratory measurements
Robert, Tanguy ULg; Hermans, Thomas ULg; Dumont, Gaël ULg et al

Conference (2013, December 06)

This contribution consists in studying the reliability of resistivity-derived temperature, for example from time-lapse electrical resistivity tomography (ERT) surveys. The idea of using temperature as a ... [more ▼]

This contribution consists in studying the reliability of resistivity-derived temperature, for example from time-lapse electrical resistivity tomography (ERT) surveys. The idea of using temperature as a quantitative tracer is growing in the hydrogeophysical community, especially to simulate geo/hydrothermal systems. However, plenty of physico-chemical processes are influenced by temperature and most of them impact directly resistivity measurements. Therefore, one needs to take them into account to retrieve quantitative temperature estimates from resistivity measurements but, up to now, it is seldom the case. The experiment we conducted consisted in simulating an ERT monitoring of heat storage in a sandy aquifer. We show that using experimental relationships between fluid electrical conductivity and temperature alone does not allow reliable temperature estimates, simply because rock-water interactions are neglected. Worst, from a certain temperature (45°C here), the bulk resistivity starts to increase with temperature although this is not expected from the experimental law. Chemical analyses made on water samples collected during the experiment highlight the importance of accounting chemical reactions (e.g. calcite precipitation with increasing temperature) occurring when temperature changes as well as their kinetics. Finally, other parameters as surface conductivity cannot always be neglected when estimating temperature from resistivity measurements. This means that retrieving reliable temperatures from bulk resistivity measurements (e.g. time-lapse ERT) requires the knowledge of water mineralization as well as the rock / soil mineralogy in order to fully integrate physico-chemical reactions between groundwater and the host rock, for example with a joint inversion scheme. [less ▲]

Detailed reference viewed: 12 (4 ULg)
Full Text
See detailMinimum gradient support and geostatistics regularization approaches for inverting time-lapse data
Nguyen, Frédéric ULg; Hermans, Thomas ULg; Robert, Tanguy ULg

Conference (2013, December 05)

Inversion of time-lapse resistivity data allows obtaining ‘snapshots’ of changes occurring in monitored systems for applications such as aquifer storage, site remediation or tracer tests. Based on these ... [more ▼]

Inversion of time-lapse resistivity data allows obtaining ‘snapshots’ of changes occurring in monitored systems for applications such as aquifer storage, site remediation or tracer tests. Based on these snapshots, one can infer qualitative information on the location and morphology of changes occurring in the subsurface but also quantitative estimates on the degree of changes in certain property such as temperature or total dissolved solid content. Analysis of these changes can provide direct insight into flow and transport processes and controlling parameters. However, the reliability of the analysis is dependent on survey geometry, measurement schemes, data error, or regularization. Except regularization, survey design parameters may be optimized prior to the monitoring survey. Regularization, on the other hand, may be chosen depending on available information collected during the monitoring. Common approaches consider smoothing model changes both in space and/or time. We here propose to use two alternative regularization approaches which may be better suited to invert time-lapse data. The first approach is the minimum gradient support (MGS) regularization, which focus the changes in tomograms snapshots. MGS will limit the occurrences of changes in electrical resistivity but will also restrict the variations of these changes inside the different zones. The second approach is based on geostatistics and requires first to derive variogram parameters for the model changes. In this contribution, we demonstrate the benefits and limitations of these regularization approaches to time-lapse data on numerical benchmarks and three case studies. [less ▲]

Detailed reference viewed: 21 (6 ULg)
See detail3D ERT monitoring of the reactivation of waste biodegradation with fresh leachate injection
Dumont, Gaël ULg; Robert, Tanguy ULg; Pilawski, Tamara et al

Conference (2013, December 04)

The aim of this study is to monitor (bio) physical processes occurring in a landfill. The experiment consists in injecting leachate towards a drain in unsaturated and not yet digested waste to reactivate ... [more ▼]

The aim of this study is to monitor (bio) physical processes occurring in a landfill. The experiment consists in injecting leachate towards a drain in unsaturated and not yet digested waste to reactivate (or activate) waste biodegradation. The target is the first 15 meters of the studied landfill subsurface. The visualization of the wet front arrival (short term effect) is crucial because we want to ensure that waste is entirely humidified to allow the reactivation of waste digestion. The second process is a long term effect consisting in the increase of the internal temperature of the landfill which is synonymous of the reactivation of biodegradation processes. We use 3D time-lapse ERT on a monthly basis to capture the decrease of electrical resistivity related to the increasing temperature. We also collect ground truth data, including distributed temperatures in a borehole to validate results. For short term effects, we monitored the wet front arrival with three 2D ERT profiles composing the 3D image, during an entire day. Preliminary results, corroborated by ground truth data, show that leachate flow in anisotropic (more rapid horizontally than vertically). So far, waste was completely humidified and slight changes of temperature occurred. [less ▲]

Detailed reference viewed: 8 (0 ULg)
Peer Reviewed
See detailGeophysical characterisation of a former waste disposal site in the context of landfill mining
Dumont, Gaël ULg; Robert, Tanguy ULg; Pilawski, Tamara et al

in EarthDoc - Near Surface Geoscience 2013 – 19th European Meeting of Environmental and Engineering Geophysics (2013, September 11)

Detailed reference viewed: 7 (0 ULg)
Peer Reviewed
See detail3D ERT Monitoring of the Reactivation of Waste Biodegradation with Fresh Leachate Injection
Robert, Tanguy ULg; Dumont, Gaël ULg; Pilawski, Tamara ULg et al

in EarthDoc - Near Surface Geoscience 2013 – 19th European Meeting of Environmental and Engineering Geophysics (2013, September 11)

Detailed reference viewed: 21 (5 ULg)
Full Text
Peer Reviewed
See detailReliability of ERT-derived Temperature - Insights from Laboratory Measurements
Robert, Tanguy ULg; Hermans, Thomas ULg; Dumont, Gaël ULg et al

in EarthDoc - Near Surface Geosciences 2013 - 19th European Meeting of Environmental and Engineering Geophysics (2013, September)

We performed laboratory measurements on fully saturated sand samples in the context of deriving reliable temperature from time-lapse electrical resistivity tomography (ERT). The experiment consisted in ... [more ▼]

We performed laboratory measurements on fully saturated sand samples in the context of deriving reliable temperature from time-lapse electrical resistivity tomography (ERT). The experiment consisted in monitoring an increase of temperature in sand samples with electrical resistivity measurements. We neglected the effect of surface conductivity since experiments showed two orders of magnitude between surface and fluid conductivities. We show that using simple linear relationship between fluid electrical conductivity and temperature alone does not allow reliable temperature estimates. Indeed, chemical analyses highlight the importance of accounting chemical reactions occurring when temperature changes, including dissolution/precipitation processes. We performed two experiments based on typical in-situ conditions. We first simulated the injection of a less conductive tap water and second, the injection of heated formation water. In the second case, minerals solubility decreases and precipitation occurs, leading to an increase of bulk resistivity. This mechanism competes with dissolution of minerals when tap water is injected, since tap water is not in equilibrium with the medium. In any case, further research is needed to fully understand the mechanisms and to develop a fully integrated law to derive better temperature estimates. [less ▲]

Detailed reference viewed: 13 (3 ULg)
Full Text
Peer Reviewed
See detailAssessing the Probability of Training Image-Based Geological Scenarios Using Geophysical Data
Hermans, Thomas ULg; Caers, Jef; Nguyen, Frédéric ULg

in Pardo-Iguzquiza, Eulogio; Guardiola-Albert, Carolina; Heredia, Javier (Eds.) et al Mathematics of Planet Earth - Proceedings of the 15th Annual Conference of the International Association for Mathematical Geosciences (2013, September)

In multiple-point statistics (MPS), the construction of training im-ages (TIs) is one of the most critical steps. Reliable geological studies may not always be available to depict with certainty what ... [more ▼]

In multiple-point statistics (MPS), the construction of training im-ages (TIs) is one of the most critical steps. Reliable geological studies may not always be available to depict with certainty what geological patterns or heterogeneity are present. In this context, geophysical techniques may provide additional information to reduce the possible large uncertainty in the understanding of prior geological scenarios. To overcome this problem, we developed a methodology to verify the consistency of geophysical data with independently-built TIs representing different plausible geological scenarios. If a TI is deemed consistent with the field geophysical survey, then in a sec-ond step we calculate a likelihood probability for each consistent TI. Our methodology starts by creating subsurface models with each TI. From these models we create synthetic geophysical data and from this synthetic data, synthetic inverted models. These models are now compared with a single inverted model obtained from the field sur-vey, allowing for our definition of what is “consistent”. To that ex-tent, we calculate the Euclidean distance between any two inverted models as well as field data and visualize the results in a 2D or 3D space using multidimensional scaling (MDS). With this technique, it is possible to verify if field cases fall in the distribution represented by synthetic cases, and thus are consistent with them. In a second step, we present a cluster analysis on the MDS-map to highlight which parameters are the most sensitive for the construction of TI. Based on this analysis, a probability of each geological scenario is computed through kernel smoothing of the densities in reduced pro-jected metric space. This approach was tested using electrical resistivity tomography as geophysical data to analyze TI scenarios for the Meuse alluvial aqui-fer (Belgium), where the lack of reliable sedimentological data lead to the definition of a multitude of geological scenarios, hence TIs. [less ▲]

Detailed reference viewed: 36 (14 ULg)
Full Text
Peer Reviewed
See detailA heat injection and pumping experiment in a gravel aquifer monitored with crosshole electrical resistivity tomography
Hermans, Thomas ULg; Wildemeersch, Samuel ULg; Jamin, Pierre ULg et al

in EarthDoc - Near Surface Geosciences 2013 - 19th European Meeting of Environmental and Engineering Geophysics (2013, September)

Thermal tracing experiments are becoming common in hydrogeology to estimate parameters governing heat transport processes and to study geothermal reservoirs. Electrical resistivity tomography (ERT) has ... [more ▼]

Thermal tracing experiments are becoming common in hydrogeology to estimate parameters governing heat transport processes and to study geothermal reservoirs. Electrical resistivity tomography (ERT) has proven its ability to monitor salt tracer tests, but few studies have investigated its performances in thermal tracing experiments. In this study, we monitor the injection and pumping of heated water using crosshole ERT in a panel crossing the main flow direction. Difference inversion time-lapse images clearly show the heterogeneous pattern of resistivity changes, and thus temperature changes, highlighting the existence of preferential flow paths in the aquifer. Comparison of temperature estimates from ERT and direct measurements in boreholes show the ability of ERT to quantify the temperatures in the aquifer and to draw the breakthrough curves of the thermal tracer with a relative accuracy. Such resistivity data may provide important information to improve hydrogeological models. Our study proves that ERT, especially crosshole ERT, is a reliable tool to follow thermal tracing experiments. It also confirms that ERT should be included to in situ techniques to characterize heat transfer in the subsurface and to monitor geothermal resources exploitation. [less ▲]

Detailed reference viewed: 36 (17 ULg)
Full Text
See detailProbability perturbation method applied to the inversion of groundwater flow models using HydroGeoSphere
Hermans, Thomas ULg; Scheidt, Céline; Caers, Jef et al

Conference (2013, April 04)

Solving spatial inverse problems in Earth Sciences remains a big challenge given the high number of parameters to invert for and the complexity of non-linear forward models. Techniques were developed to ... [more ▼]

Solving spatial inverse problems in Earth Sciences remains a big challenge given the high number of parameters to invert for and the complexity of non-linear forward models. Techniques were developed to reduce the number of parameters to invert for or to produce geologically consistent simulations from an initial guess. These techniques ask for a prior model to constrain the spatial distribution of the solution. Geostatistical models contain, by nature, information to control the spatial features of the inverse solutions, but the integration of dynamic data into such models remains difficult. We adapted, the “probability perturbation algorithm” (PPM) using Matlab® to invert hydrogeological data using multiple-point geostatistics to build models of pre-defined hydrofacies. The algorithm uses HydroGeoSphere (HGS) to compute the forward response of the model and SGems to produce geostatistical realizations. The algorithm only needs the proper definition of all the parameters to be used by HydroGeoSphere (grid matching with SGems, position of the wells, pumping rate, facies properties, boundary conditions, etc.). The PPM algorithm will automatically seek solutions fitting both hydrogeological data and geostatistical constraints. Through the inversion process, the initial geostatistical realization is perturbed. Only geometrical features of the model are affected, i.e. we do not attempt to directly find the optimal value of hydrogeological parameters, but the optimal spatial distribution of facies whose prior distribution is quantified in a training image. The algorithm can be divided in three steps. In the first step, we use SGems to generate an initial facies model with the multiple-point geostatistical algorithm SNESIM (single normal equation simulation). The facies model is composed of several categories representing hydrological facies (e.g. gravel, sand and clay). It can be conditioned using hard data (borehole data) and/or soft data (e.g. geophysical data). We then run a first flow simulation with HydroGeoSphere. This requires defining hydrogeological parameters (porosity, hydraulic conductivity, etc.) for each category of the facies model to create a hydrogeological model. The response of the latter model is compared to the expected one through an objective function. In the second step, a perturbation to the facies model is computed using a single parameter called rD. This perturbation is used to generate a new facies model with SGems and calculate a new objective function value via HGS, as done in the first step. An inner loop optimizes the value of rD. In the third step, we verify if the objective function of the best fitting model is smaller than a predefined value. If it is the case, we stop the algorithm, otherwise we go back to step 2 until convergence. We illustrate the methodology with a synthetic example in an alluvial aquifer. The model is based on a training image depicting gravel channels and clay lenses in a coarse sand aquifer. We simulate a pumping test and inverse water level data recorded at 9 wells using our implementation of the PPM algorithm. Using this method, it is possible to generate multiple solutions and to derive a posterior probability of the facies distribution. [less ▲]

Detailed reference viewed: 104 (8 ULg)
See detailHydrogeological processes in fractured and porous media: insights from geophysical case studies
Robert, Tanguy ULg; Hermans, Thomas ULg; Nguyen, Frédéric ULg

Conference (2013, January 18)

This presentation focuses on geophysical case studies with the aim to highlight the possibilities to study and monitor hydrogeological processes in the subsurface, including transport processes in ... [more ▼]

This presentation focuses on geophysical case studies with the aim to highlight the possibilities to study and monitor hydrogeological processes in the subsurface, including transport processes in fractured or in porous media. The presentation emphasizes two geoelectrical methods, namely electrical resistivity tomography (ERT) which images the electrical resistivity distribution of the subsurface and self-potential (SP) whose measured signal is directly sensitive to groundwater fluxes. The first case study concerns the geophysical identification and characterization of large hydraulically-active fractured areas in calcareous synclines and in particular the assessment of the joint use of ERT and SP to set up new piezometers in fractured limestone. This assessment shows that piezometers drilled inside less resistive areas and/or in negative SP anomalies presented high hydraulic capacities. Inversely, piezometers drilled inside more resistive zones and/or outside an SP anomaly presented low hydraulic capacities. The SP anomaly related to preferential flow in fractures was thus demonstrated for the first time. All these fractures information, obtained with geophysics, improved the conceptualization and calibration of the groundwater flow model of the calcareous valley. A seasonal monitoring of SP signals proved to be a successful methodology to better understand the hydrodynamics of calcareous aquifers and in particular to follow the seasonal drawdown of the water table in the calcareous valley. Different methodologies to delineate the main groundwater flow direction were also tested. The latter can be achieved for example by drawing an SP map showing the main hydraulic gradients or by monitoring a salt tracer test with ERT to highlight preferential flow in fractures. The second case study concerns the ERT monitoring of a shallow geothermal test conducted in a porous medium (sand). The main objective of this study was to derive temperature from a series of electrical resistivity images since the electrical resistivity is directly sensitive to temperature changes. This field work demonstrates that surface electric resistivity tomography can monitor heat injection and storage experiments in shallow aquifers providing a number of practical applications, such as the monitoring or the design of shallow geothermal systems or the use of heated water to replace salt water in tracer tests. Through these two different case studies, this presentation also emphasizes in a practical way on the importance of data inversion and image appraisal since these issues are crucial to quantitatively study hydrogeological processes. [less ▲]

Detailed reference viewed: 59 (11 ULg)
Full Text
Peer Reviewed
See detailA comparison study of image appraisal tools for electrical resistivity tomography
Caterina, David ULg; Beaujean, Jean ULg; Robert, Tanguy ULg et al

in Near Surface Geophysics (2013)

To date, few studies offer a quantitative comparison of the performance of image appraisal tools. Moreover, there is no commonly accepted methodology to handle them even though it is a crucial aspect for ... [more ▼]

To date, few studies offer a quantitative comparison of the performance of image appraisal tools. Moreover, there is no commonly accepted methodology to handle them even though it is a crucial aspect for reliable interpretation of geophysical images. In this study, we compare quantitatively different image appraisal indicators to detect artefacts, estimate depth of investigation, address parameters resolution and appraise ERT-derived geometry. Among existing image appraisal tools, we focus on the model resolution matrix (R), the cumulative sensitivity matrix (S) and the depth of investigation index (DOI) that are regularly used in the literature. They are first compared with numerical models representing different geological situations in terms of heterogeneity and scale and then used on field data sets. The numerical benchmark shows that indicators based on R and S are the most appropriate to appraise ERT images in terms of the exactitude of inverted parameters, DOI providing mainly qualitative information. In parallel, we test two different edge detection algorithms – Watershed’s and Canny’s algorithms – on the numerical models to identify the geom-etry of electrical structures in ERT images. From the results obtained, Canny’s algorithm seems to be the most reliable to help practitioners in the interpretation of buried structures. On this basis, we propose a methodology to appraise field ERT images. First, numerical bench¬mark models representing simplified cases of field ERT images are built using available a priori information. Then, ERT images are produced for these benchmark models (all simulated acquisition and inversion parameters being the same). The comparison between the numerical benchmark mod¬els and their corresponding ERT images gives the errors on inverted parameters. These discrepan¬cies are then evaluated against the appraisal indicators (R and S) allowing the definition of threshold values. The final step consists in applying the threshold values on the field ERT images and to validate the results with a posteriori knowledge. The developed approach is tested successfully on two field data sets providing important information on the reliability of the location of a contamina¬tion source and on the geometry of a fractured zone. However, quantitative use of these indicators remains a difficult task depending mainly on the confidence level desired by the user. Further research is thus needed to develop new appraisal indicators more suited for a quantitative use and to improve the quality of inversion itself. [less ▲]

Detailed reference viewed: 50 (27 ULg)