Navigation – Plan du site

AccueilNuméros33 (suppl.)Methods and innovations2D Electrical Resistivity Tomogra...

Methods and innovations

2D Electrical Resistivity Tomography survey optimisation of solute transport in porous media

Gregory Lekmine, Marc Pessel et Harold Auradou
p. 309-312

Entrées d’index

Haut de page

Texte intégral

1The purpose of this study is to quantify experimentally the evolution of dissolved species in porous media from 2D resistivity models. Transport experiments are done at the lab level by performing flow in a model porous medium obtained by filling a transparent container with mono disperse glass beads. A tracer made by mixing a dissolved blue dye and a known NaCl concentration is injected with a constant flow rate through the porous medium already saturated by a transparent fluid. ERT measurements are acquired during the fluids flow. The measurement conditions and the inversion parameters are estimated so that the relation between spatial and temporal resolutions is optimised. A video follow-up is also carried out during the flow. A comparison of the temporal evolution of the NaCl concentration distribution estimated from ERT models with video analysis shows remarkable agreement.

2Electrical resistivity tomography is used to detect the contaminant’s plumes in saturated porous media using the resistivity contrasts between the fluids saturating initially the pore spaces and the moving fluids (Binley, 1996; Oldenborger et al., 2007). Archie’s law (1942) relates the measured bulk apparent conductivities to the fluids resistivity in the poral network via a priori information about the porosity and the particles cementation state. As previously shown by (Slater et al., 2000), it is possible to use this method to estimate the tracer spatial reparation in terms of concentrations. This method is non-invasive and non-destructive, but it is based on indirect measurements of the electrical properties. The electrical conductivity distribution is obtained through data inversion. However the non-unicity of solution can lead to misinterpretations of the geological reality. We test the model’s reliability in laboratory conditions by experiments in porous media at mesoscopic level where the hydrodynamic parameters can be easily estimated. We optimise the protocols of measurements and inversions according to the medium’s geometry for tracer follow-up. The built models are then expressed in terms of tracer concentration distribution on several time steps and compared to the estimated concentrations by a parallel video follow-up.

Experiment protocol

3The flow experiments are carried out in a transparent Plexiglas container (with a volume of V=HxLxE=233.75 cm3) filled with glass beads of 166 microns diameter. Porosity is F = 0.365, close to the theoretical value estimated by (Gondret, 1994): the porous medium may be considered homogeneous and isotropic. Hydrodynamic parameters are supposed constant, and the edge effects on the porosity are neglected since the pores diameter is small in comparison to the wall dimensions (Gondret, 1994). A permeable nylon fabric of 100 μm mesh covers the lower opening of the container. This part is in contact with a tank which contains the tracer, whereas three injectors in the sealed upper cover ensure a constant flow rate (Fig. 1a). Two lines of twenty one electrodes of 1.5 mm diameter, spaced every 1 cm are installed along the edges of the cell (Fig. 1b). The first and the last electrode of each line are located respectively at 3.5 cm of the lower and upper opening of the container. This kind of electrodes configuration reproduces the measurements conditions used in a borehole.

Figure 1a: Schematic of the experimental Plexiglas cell filled with spherical glass beads of 125 to 166 micron diameter. H=27.5 cm, L=8.5 cm, E=1 cm.

Figure 1a: Schematic of the experimental Plexiglas cell filled with spherical glass beads of 125 to 166 micron diameter. H=27.5 cm, L=8.5 cm, E=1 cm.

Figure 1b: Side view of the experimental device. s=1 cm, a=3.5 cm, b=21 cm.

Figure 1b: Side view of the experimental device. s=1 cm, a=3.5 cm, b=21 cm.

4Four experimental conditions must be respected in order to optimise the ERT follow-up of tracer propagation (Slater et al., 2000):

  1. all the pore space is saturated;

  2. the matrix conductivity is low in comparison to tracer conductivity;

  3. temporal variations of electric conductivity are only due to conductivity changes of interstitial fluids;

  4. electrical conductivity is linearly related to tracer concentration.

5Prior to any experiments, the tank in contact with the lower opening is filled by a transparent fluid, which is pumped upwards by the injectors until total saturation of the medium. This fluid is produced by dissolution of a small amount of NaCl in degassed and demineralised water.The pumping rate (1.5 ml/min) is sufficiently low to obtain a perfect saturation of the medium. The tracer is obtained by dissolution of 0.1 g/L blue dye mixed with 1.0 g/L NaCl in demineralised water.

6The two fluids are miscible, have the same viscosity and their different NaCl concentration creates a small density difference. During the experiments a mixing zone appeared between the two fluids: its characteristics depend upon the densities contrasts, pore size and flow rate.

7The geometry evolution of this zone is observed, at the same time, by two independent techniques: the video follow-up which links transmitted light intensity through the cell measured on each pixel of the camera to the dye concentration, and the electrical measurements which reveal the NaCl dispersion during transport. The experiments are reproduced at several constant flow rates with different NaCl concentrations between the two fluids and different density differences.

Electrical measurement protocol optimisation

8A sensitivity analysis (Furman et al., 2007) is necessary to estimate the best series of quadripoles to describe the dispersion front on each time-step, while minimizing the acquisition time. First we did a set of experiments where the pumping was stopped when the tracer filled half of the porous medium. In these conditions gravity stabilizes the front and its position is fixed. We have tested several standard protocol measurements (Wenner-Schlumberger, gradient, dipole-dipole). The first results show that the measurements are strongly perturbed by the resistive edges of the plexiglas tank. In order to remove this effect, data have to be normalised (Daily et al., 1995; Nimmer et al., 2008) and the inversion parameters are optimised to produce reliable models of conductivity distribution. Archie’s law and linear relationship between conductivity and concentration (Fig. 2) allows to relate the bulk conductivity of each model to NaCl concentration. Vertical variations of dye concentrations estimated by video (red curve) and NaCl estimated by ERT (black curve) are compared. For all experiments, the two concentration curves are shifted: the mean position of the NaCl mixing front appears to be in advance in comparison to the dye. This shift is probably due to dye trapping in dead end pore zones and adsorption on particle surfaces.

Figure 2: Vertical variations concentrations of NaCl (C/C0) and dye obtained by ERT (black curve) and video analysis (red curve).

Figure 2: Vertical variations concentrations of NaCl (C/C0) and dye obtained by ERT (black curve) and video analysis (red curve).

Dynamic monitoring of the mixing front

9Tests carried out for a static front position allow the efficiency of each measurement protocol to be evaluated in description on the fluids interface. We found that a transversal dipole-dipole configuration with 190 quadripoles is best adapted to the geometry of our medium. The temporal resolution is also optimised: each time step is performed in 5 minutes while 30 minutes were needed for the measurements in static conditions. Yet, the continuous tracer propagation causes a temporal shift between the first and the last measurements of each time step. This drawback is corrected assuming a linear conductivity variation with time. Data analyses and inversion parameters are similar as in the case of static fluids (Fig. 3), vertical variations of concentrations estimated by both methods present the same tendencies. Thickness and mean front position on first time steps are similar (Fig. 3), but the tracer propagation creates a shift of the mean front position which increases upward during the experiment.

Figure 3: Dynamic tracer follow-up in terms of electrical resistivity distribution and comparison with relative concentrations estimated on each time step.

Figure 3: Dynamic tracer follow-up in terms of electrical resistivity distribution and comparison with relative concentrations estimated on each time step.

Conclusions

Optimisation of the electrical measurements protocol provides an acceptable compromise between spatial and temporal resolutions.

Dynamics monitoring gives a reliable estimation of position and thickness of the mixing front.

The shift observed between video and ERT is caused by adsorption and dye trapping, which may be quantified by our method.

10In the future, it will be necessary to check on the repetitiveness of the results with several density differences and flow rates. The system will also be complexed, while adding permeability heterogeneities and a 2D vertical fracture.

Haut de page

Bibliographie

Archie, G. E., 1942. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics. J. Petroleum Technol. Tech. Pub., 5 p.

Binley, A., Henry-Poulter, S. and Shaw, B., 1996. Examination of solute transport in an undisturbed soil column using electrical resistance tomography. Water Resources Research, 32: 763-769.

Flowers, T. C. and Hunt, J. R., 2007. Viscous and gravitational contributions to mixing during vertical brine transport in water-saturated porous media. Water Resources Research, 43, W01407.

Furman, A., Ferré, T. P. A. and Health, G. L., 2007. Spatial focusing of electrical resistivity surveys considering geologic and hydrologic layering. Geophysics, 72: 65-73.

Kemna, A., Vanderborght, J., Kulessa, B. and Vereecken, H., 2002. Imaging and characterisation of subsurface solute transport using electrical resistivity (ERT) tomography and equivalent transport models. Journal of Hydrology, 267: 125-165.

Oldenborger, G. A., Knoll, M. D., Routh, P. S. and LaBrecque, D. J., 2007. Time-lapse ERT monitoring of an injection/withdrawal experiment in a shallow unconfined aquifer. Geophysics, 72: 177-187.

Slater, L., Binley, A. M., Daily, W. and Johnson, R., 2000. Cross-hole electrical imaging of a controlled saline tracer injection. Journal of Applied Geophysics, 44: 85-102.

Haut de page

Table des illustrations

Titre Figure 1a: Schematic of the experimental Plexiglas cell filled with spherical glass beads of 125 to 166 micron diameter. H=27.5 cm, L=8.5 cm, E=1 cm.
URL http://journals.openedition.org/archeosciences/docannexe/image/1756/img-1.png
Fichier image/png, 63k
Titre Figure 1b: Side view of the experimental device. s=1 cm, a=3.5 cm, b=21 cm.
URL http://journals.openedition.org/archeosciences/docannexe/image/1756/img-2.png
Fichier image/png, 92k
Titre Figure 2: Vertical variations concentrations of NaCl (C/C0) and dye obtained by ERT (black curve) and video analysis (red curve).
URL http://journals.openedition.org/archeosciences/docannexe/image/1756/img-3.png
Fichier image/png, 71k
Titre Figure 3: Dynamic tracer follow-up in terms of electrical resistivity distribution and comparison with relative concentrations estimated on each time step.
URL http://journals.openedition.org/archeosciences/docannexe/image/1756/img-4.jpg
Fichier image/jpeg, 406k
Haut de page

Pour citer cet article

Référence papier

Gregory Lekmine, Marc Pessel et Harold Auradou, « 2D Electrical Resistivity Tomography survey optimisation of solute transport in porous media »ArcheoSciences, 33 (suppl.) | 2009, 309-312.

Référence électronique

Gregory Lekmine, Marc Pessel et Harold Auradou, « 2D Electrical Resistivity Tomography survey optimisation of solute transport in porous media »ArcheoSciences [En ligne], 33 (suppl.) | 2009, mis en ligne le 30 octobre 2011, consulté le 29 mars 2024. URL : http://journals.openedition.org/archeosciences/1756 ; DOI : https://doi.org/10.4000/archeosciences.1756

Haut de page

Auteurs

Gregory Lekmine

Laboratoire IDES- UMR 8148, Université Paris Sud bâtiment 504 – 91405 Orsay Cedex

Marc Pessel

Laboratoire IDES- UMR 8148, Université Paris Sud bâtiment 504 – 91405 Orsay Cedex

Harold Auradou

Laboratoire FAST-UMR 7608, Universités Paris Sud & Paris VI bâtiment 502 – 91405 Orsay Cedex

Haut de page

Droits d’auteur

CC-BY-NC-ND-4.0

Le texte seul est utilisable sous licence CC BY-NC-ND 4.0. Les autres éléments (illustrations, fichiers annexes importés) sont « Tous droits réservés », sauf mention contraire.

Haut de page
Rechercher dans OpenEdition Search

Vous allez être redirigé vers OpenEdition Search