The potential of time-lapse GPR full-waveform inversion as high resolution imaging technique for salt and ethanol transport

Crosshole GPR full-waveform inversion (FWI) has shown a high potential to characterize the near surface at a decimeter scale which is crucial for flow and transport. GPR FWI provide high-resolution tomograms of dielectric permittivity and electrical conductivity, which can be linked lithological properties. This study tests the potential of time-lapse GPR FWI to monitor tracers of different geophysical properties (salt, heat, ethanol). Synthetic and preliminary field results show that both properties can resolve major transport processes

Advancing measurements and representations of subsurface heterogeneity and dynamic processes: towards 4D hydrogeology

Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability generally leads to effective behaviors and emerging phenomena that cannot be predicted from conventional approaches based on homogeneous assumptions and models. However, it is not always clear when, why, how, and at what scale the 4D (3D + time) nature of the subsurface needs to be considered in hydrogeological monitoring, modeling, and applications. In this paper, we discuss the interest and potential for the monitoring and characterization of spatial and temporal variability, including 4D imaging, in a series of hydrogeological processes: (1) groundwater fluxes, (2) solute transport and reaction, (3) vadose zone dynamics, and (4) surface–subsurface water interactions. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. We then highlight recent innovations that have led to significant breakthroughs in high-resolution space–time imaging and modeling the characterization, monitoring, and modeling of these spatial and temporal fluctuations. We finally propose a classification of processes and applications at different scales according to their need and potential for high-resolution space–time imaging. We thus advocate a more systematic characterization of the dynamic and 3D nature of the subsurface for a series of critical processes and emerging applications. This calls for the validation of 4D imaging techniques at highly instrumented observatories and the harmonization of open databases to share hydrogeological data sets in their 4D components

Individual and joint inversion of head and flux data by geostatistical hydraulic tomography

Hydraulic tomography is a state-of-the-art method for inferring hydraulic conductivity fields using head data. We employed geostatistical inversion using synthetically generated head and flux data individually and jointly in a steady-state experiment. We designed 96 inversion scenarios to better understand the relative merits of each data type. For the typical case of a small number of observation points, we find that flux data provide a better resolved hydraulic conductivity field compared to head data when considering data with similar signal-to-noise ratios. This finding is further confirmed by a resolution analysis. When considering a high number of observation points, the estimated fields are of similar quality regardless of the data type. In terms of borehole boundary conditions, the best setting for flux and head data are constant head and constant rate, respectively, while joint inversion results are insensitive to the borehole boundary type. When considering the same number of observations, the joint inversion of head and flux data does not offer advantages over individual inversions. When considering the same number of observation points and, hence, twice as many observations, the joint inversion performs better than individual inversions. The findings of this paper are useful for future planning and design of hydraulic tomography tests comprising flux and head data

An assessment of the relative information content of ground water flux and pressure data in the context of geostatistical inversion

International audienceSubsurface characterization often relies on inversion of either pressure or tracer data. Unless data from many pumping and observation wells are available, the inversion process only resolves smooth low-resolution images of subsurface properties, which leads to less accurate subsurface flow and reactive transport predictions. Furthermore, tracer tomography can be very challenging and convergence to a global minimum is difficult. Active-distributed temperature sensing technology opens up the prospect of replacing tracer test data with estimates of subsurface groundwater flux [1]. Here, the value of using estimated subsurface groundwater fluxes as a data source to reconstruct subsurface hydraulic properties is explored using a sequence of synthetic multivariate Gaussian aquifers with different measurement configurations. These results are compared to inversion of pressure data and joint inversion of the two data types with the inversions being based on the Principal Component Geostatistical Approach [2]. Inversion of pressure data resulted in a smoothed reconstruction of aquifer heterogeneity capturing approximately high and low conductivity regions while ground water flux data inversion leads to higher-resolution estimates. This is reflected, for one of the considered examples, by a correlation coefficient that increases from 0.57 for the pressure data to 0.65 for the ground water flux data. The complimentary nature of the data sets is represented by a correlation coefficient that increases to 0.74 for the joint inversion of the two data types.To conclude, inversion of ground water flux whether individually or jointly with pressure data, can provide enhanced information about the heterogeneity of subsurface media compared with using pressure data alone

Use of Groundwater Flux Data Measured by Active Fiber Optic DTS for Reconstruction of Aquifers Heterogeneity

International audienceSubsurface characterization often relies on inversion of either pressure or tracer data. Unless data from many pumping and observation wells are available, the inversion process only resolves smooth low-resolution images of subsurface properties, which leads to less accurate subsurface flow and reactive transport predictions. Furthermore, tracer tomography can be very challenging and convergence to a global minimum is difficult. Active-distributed temperature sensing technology opens up the prospect of replacing tracer test data with estimates of subsurface groundwater flux. Here, the value of using estimated subsurface groundwater fluxes as a data source to reconstruct subsurface hydraulic properties is explored using a sequence of synthetic multivariate Gaussian aquifers with different measurement configurations. These results are compared to inversion of pressure data and joint inversion of the two data types with the inversions being based on the Principal Component Geostatistical Approach. Inversion of pressure data resulted in a smoothed reconstruction of aquifer heterogeneity capturing approximately high and low conductivity regions while ground water flux data inversion leads to higher-resolution estimates. To conclude, inversion of ground water flux whether individually or jointly with pressure data, can provide enhanced information about the heterogeneity of subsurface media compared with using pressure data alone

Modelling borehole flows from Distributed Temperature Sensing data to monitor groundwater dynamics in fractured media

International audienceFractured aquifers are known to be very heterogeneous with complex flow path geometries. Their characterization and monitoring remain challenging despite the importance to better understand their behavior at all spatial and temporal scales. Heat and correspondingly temperature data have gained much interest in recent years and are often used as a tracer for characterizing groundwater flows. In the current work, a fast computer code is developed using Ramey and Hassan and Kabir analytical solutions which converts the temperature profile to the flow rate profile along the borehole. The method developed is validated through numerical simulations. A global sensitivity study recognizes the media thermal properties as the most influential parameters. For testing the method in the field, fiber-optic distributed temperature sensing (FO-DTS) data were used to monitor the dynamic behavior of fractured aquifers at the borehole scale at the Ploemeur-Guidel field site in Brittany, France. DTS data are used to infer the flow rates in the different sections of a fractured wellbore (flow profile) and calculate the contribution of each fracture to the total flow. DTS data were acquired for about three days in three different hydraulic conditions corresponding to two different ambient flow conditions and one pumping condition. Flow profiling using distributed temperature data matches satisfactorily with results from heat-pulse flow metering performed in parallel for cross-checking. Moreover, flow profiling reveals the daily variations of ambient flow in this fractured borehole. Furthermore, it shows that during ambient flowing conditions, shallow and deep fractures contribute roughly equally to the total flow while during the pumping condition, the deepest fractures contribute more to the total flow, suggesting a possible reorganization of flow and hydraulic heads depending on the hydraulic conditions. Thus, although the proposed method (DTS data and proposed framework) may be costlier and is based on indirect characterization through temperature measurements, it provides real-time monitoring of complex fracture interactions and recharge processes in fractured media. Thus, this method allows for a full analysis of the temporal behavior of the system with a simple and fast analytical model. Furthermore, thanks to its narrow width, DTS can be used and installed in boreholes for long-term monitoring while heat-pulse flow metering may lead to head losses in the borehole and may not be always possible depending on some borehole conditions. One of the limitations the approach proposed is the proper knowledge of the thermal properties of media required to infer the flow rate from the temperature. Nevertheless, surface rate measurement can be useful to constrain these properties and reduce the flow profiling uncertainty. Thus, the method proposed appears to be an interesting and complementary method for characterizing borehole flows and groundwater dynamics in fractured media such as for instance, monitoring the recharge dynamic

Testing in sandbox experiments the potentialities of active-Distributed Temperature Sensing to quantify distributed groundwater fluxes in porous media

International audienceActive-Distributed Temperature Sensing is a new method that has been recently developed for quantifying groundwater fluxes in the sub-surface along fibre-optic cables with a great spatial resolution. It consists in measuring and modelling the increase of temperature due to a heat source, dissipated through heat conduction and heat advection, depending on groundwater fluxes. Here, we propose to estimate the applicability and limitations of the method using sandbox experiments where flow rate and temperature are well controlled. For doing so, active-DTS experiments have been achieved under different flow rates and experimental conditions. In addition, we compare three different and complementary methods to estimate in practice the spatial resolution of DTS measurements. Active-DTS experiments have been conducted by deploying a fiber optic cable in a large PVC tank (1.6m long; 1.2 m width and 0.3 m height) and filled with 0.4-1.3 mm diameter sand. The height of water in water reservoirs on either side of the sandbox can be adjusted to control the head gradient and the flow rate through the sand. Heating was done by injecting during at least 8 hours for each experiment, a well-controlled electrical current along the steel armouring of the fiber optic cable. The three methods for estimating spatial resolution were applied and compared using FO-DTS measurements obtained on the same fiber-optic cable but with two different DTS units having different spatial resolution. Results show that a large range of groundwater fluxes may be estimated with a very good accuracy. Finally, we compare the advantages and complementarities of the different methods proposed for estimating the spatial resolution of measurements. In particular, the spatial resolution estimated using a temperature step change is both dependent on the effective spatial resolution of the DTS unit but also on heat conduction induced because of the high thermal conductivity of the cable. By showing the applicability of the method for a large range of flow rates and with an excellent spatial resolution, these experiments demonstrate the potentialities of the method for quantifying fluid fluxes in porous media for a large range of applications

A Comparison of Different Methods to Estimate the Effective Spatial Resolution of FO-DTS Measurements Achieved during Sandbox Experiments

International audienceFor many environmental applications, the interpretation of fiber-optic Raman distributed temperature sensing (FO-DTS) measurements is strongly dependent on the spatial resolution of measurements, especially when the objective is to detect temperature variations over small scales. Here, we propose to compare three different and complementary methods to estimate, in practice, the “effective” spatial resolution of DTS measurements: The classical “90% step change” method, the correlation length estimated from experimental semivariograms, and the derivative method. The three methods were applied using FO-DTS measurements achieved during sandbox experiments using two DTS units having different spatial resolutions. Results show that the value of the spatial resolution estimated using a step change depends on both the effective spatial resolution of the DTS unit and on heat conduction induced by the high thermal conductivity of the cable. The correlation length method provides an estimate much closer to the value provided by the manufacturers, representative of the effective spatial resolutions along cable sections where temperature gradients are small or negligible. Thirdly, the application of the derivative method allows for verifying the representativeness of DTS measurements all along the cable, by localizing sections where measurements are representative of the effective temperature. We finally show that DTS measurements could be validated in sandbox experiments, when using devices with finer spatial resolution

Numerical and experimental validation of the applicability of active‐DTS experiments to estimate thermal conductivity and groundwater flux in porous media

International audienceGroundwater flow depends on the heterogeneity of hydraulic properties whose field characterization is challenging. Recently developed active‐Distributed Temperature Sensing (DTS) experiments offer the possibility to directly measure groundwater fluxes resulting from heterogeneous flow fields. Here, based on fundamental principles and numerical simulations, two interpretation methods of active‐DTS experiments are proposed to estimate both the porous media thermal conductivities and the groundwater fluxes in sediments. These methods rely on the interpretation of the temperature increase measured along a single heated fiber optic (FO) cable and consider heat transfer processes occurring both through the FO cable itself and through the porous media. The first method relies on the Moving Instantaneous Line Source (MILS) model that reproduces the temperature increase and provides estimates of thermal conductivity and groundwater flux as well as an evaluation of the temperature rise due to the FO cable. The second method, based on the graphical identification of three characteristic times, provides complementary estimates of flux, fully independent of the effect of the FO cable. Sandbox experiments provide an experimental validation of the interpretation methods, demonstrate the excellent accuracy of groundwater flux estimates (< 5%) and highlight the complementarity of both methods. Active‐DTS experiments allow investigating groundwater fluxes over a large range spanning 1x10‐6 to 5x10‐2 m/s, depending on the duration of the experiment. Considering the applicability of active‐DTS experiments in different contexts, we propose a general experimental framework for the application of both interpretation methods in the field, making active‐DTS field experiments especially promising for many subsurface applications

The potential of time-lapse GPR full-waveform inversion as high resolution imaging technique for salt and ethanol transport

Cross-hole GPR full-waveform inversion (FWI) has shown a high potential to characterize the near surface at a decimeter scale which is crucial for flow and transport. GPR FWI provide high-resolution tomograms of dielectric permittivity (ε) and electrical conductivity (EC), which can be linked lithological properties. This study tests the potential of time-lapse GPR FWI to monitor tracers of different geophysical properties (salt, heat, ethanol). Synthetic and field results show that both properties can resolve the plume at decimeter- channelized transport- scale