Search

Predicting the permeability of porous media in saturated and partially saturated conditions is of crucial importance in many geo-engineering areas, from water resources to vadose zone hydrology or contaminant transport predictions. Many models have been proposed in the literature to estimate the permeability from properties of the porous media such as porosity, grain size or pore size. In this study, we develop a model of the permeability for porous media saturated by one or two fluid phases with all physically-based parameters using a fractal upscaling technique. The model is related to microstructural properties of porous media such as fractal dimension for pore space, fractal dimension for tortuosity, porosity, maximum radius, ratio of minimum pore radius and maximum pore radius, water saturation and irreducible water saturation. The model is favorably compared to existing and widely used models from the literature. Then, comparison with published experimental data for both unconsolidated and consolidated samples, we show that the proposed model estimate the permeability from the medium properties very well., Comment: Hydrogeology Journal, Springer Verlag, 2021

SUMMARY Geophysical imaging using the inversion procedure is a powerful tool for the exploration of the Earth's subsurface. However, the interpretation of inverted images can sometimes be difficult, due to the inherent limitations of existing inversion algorithms, which produce smoothed sections. In order to improve and automate the processing and interpretation of inverted geophysical models, we propose an approach inspired from data mining. We selected an algorithm known as DBSCAN (Density-Based Spatial Clustering of Applications with Noise) to perform clustering of inverted geophysical sections. The methodology relies on the automatic sorting and clustering of data. DBSCAN detects clusters in the inverted electrical resistivity values, with no prior knowledge of the number of clusters. This algorithm has the advantage of being defined by only two parameters: the neighbourhood of a point in the data space, and the minimum number of data points in this neighbourhood. We propose an objective procedure for the determination of these two parameters. The proof of concept described here is applied to simulated ERT (electrical resistivity tomography) sections, for the following three cases: two layers with a step, two layers with a rebound, and two layers with an anomaly embedded in the upper layer. To validate this approach, sensitivity studies were carried out on both of the above parameters, as well as to assess the influence of noise on the algorithm's performance. Finally, this methodology was tested on real field data. DBSCAN detects clusters in the inverted electrical resistivity models, and the former are then associated with various types of earth materials, thus allowing the structure of the prospected area to be determined. The proposed data-mining algorithm is shown to be effective, and to improve the interpretation of the inverted ERT sections. This new approach has considerable potential, as it can be applied to any geophysical data represented in the form of sections or maps.

Precipitation and dissolution are prime processes in carbonate rocks and being able to monitor them is of major importance for aquifer and reservoir exploitation or environmental studies. Electrical conductivity is a physical property sensitive both to transport phenomena of porous media and to dissolution and precipitation processes. However, its quantitative use depends on the effectiveness of the petrophysical relationship to relate the electrical conductivity to hydrological properties of interest. In this work, we develop a new physically-based model to estimate the electrical conductivity by upscaling a microstructural description of water-saturated fractal porous media. This model is successfully compared to published data from both unconsolidated and consolidated samples, or during precipitation and dissolution numerical experiments. For the latter, we show that the permeability can be linked to the predicted electrical conductivity.

Self-potential signals can be generated by different sources and can be decomposed in various contributions. Streming potential is the contribution due to the water flux in the subsurface and is of particular interest in hydrogeophysics and reservoir characterization. Being able to estimate water fluxes in porous and fractured media using streaming potential data relies on our understanding of the electrokinetic coupling at the mineral-solution interface and our capacity to understand, model, and upscale this phenomenon. Two main approaches have been proposed to predict streaming potential generation in geological media. One of these approaches is based on determining the excess charge which is effectively dragged in the medium by water flow. In this chapter, we describe how to model the streaming potential by considering this effective excess charge density, how it can be defined, calculated and upscaled. We provide a short overview of the theoretical basis of this approach and we describe different applications to both water saturated and partially saturated soils and fractured media., Comment: 36 pages, 14 figures, book chapter

The self-potential (SP) method is a passive geophysical technique, which may offer insights about water and ionic fluxes in the vadose zone. The main obstacles presently prohibiting its routine use in quantitative vadose zone hydrology are the superposition of signals arising from various source mechanisms, difficult-to-predict electrode polarization effects that depend on electrode design and age, as well as water saturation, pore water chemistry, clay content, and temperature in the immediate vicinity of the electrodes. We present a unique long-term SP monitoring experiment focusing on the first four years of data acquired at different depths in the vadose zone within the HOBE hydrological observatory in Denmark. Using state-of-the-art SP theory combined with flow and transport simulations, we attempt to replicate the observed data and suggest reasons for observed discrepancies. The predictions are overall satisfactory during the first six months of monitoring after which both the patterns and magnitudes of the observed data change drastically. Our main observations are (1) that predicted SP magnitudes are strongly sensitive to how the effective excess charge scales with water saturation implying that continued research is needed to build more accurate models of electrokinetic phenomena in unsaturated conditions, (2) that significant changes in electrode polarization occur in the shallowest electrodes at time scales of a year, suggesting that electrode effects cannot be ignored and that explicit electrode modeling should be considered in future monitoring studies, and (3) that multi-rate mass transfer and reactive transport modeling are needed to better predict salinity and pore water conductivity. We hope to stimulate other researchers to test new SP modeling approaches and interpretation strategies against these data by making the SP and complimentary data time-series available., Comment: 17 pages, 12 figures, 2 tables

Electrical conductivity is one of the most commonly used geophysical method for reservoir and environmental studies. Its main interest lies in its sensitivity to key properties of storage and transport in porous media. Its quantitative use therefore depends on the efficiency of the petrophysical relationship to link them. In this work, we develop a new physically based model for estimating electrical conductivity of saturated porous media. The model is derived assuming that the porous media is represented by a bundle of tortuous capillary tubes with a fractal pore-size distribution. The model is expressed in terms of the porosity, electrical conductivity of the pore liquid and the microstructural parameters of porous media. It takes into account the interface properties between minerals and pore water by introducing a surface conductivity. Expressions for the formation factor and hydraulic tortuosity are also obtained from the model derivation. The model is then successfully compared with published data and performs better than previous models. The proposed approach also permits to relate the electrical conductivity to other transport properties such as the hydraulic conductivity., Comment: 24 pages, 10 figures, 3 tables

Movement of soil moisture associated with tree root-water uptake is ecologically important, but technically challenging to measure. Here, the self-potential (SP) method, a passive electrical geophysical method, is used to characterize water flow in situ. Unlike tensiometers, which use a measurement of state (i.e. matric pressure) at two locations to infer fluid flow, the SP method directly measures signals generated by water movement. We collected SP measurements in a two-dimensional array at the base of a Douglas-fir tree (Pseudotsuga menziesii) in the H.J. Andrews Experimental Forest in western Oregon over five months to provide insight on the propagation of transpiration signals into the subsurface under variable soil moisture. During dry conditions, SP data appear to show downward unsaturated flow, while nearby tensiometer data appear to suggest upward flow during this period. After the trees enter dormancy in the fall, precipitation-induced vertical flow dominates in the SP and tensiometer data. Diel variations in SP data correspond to periods of tree transpiration. Changes in volumetric water content occurring from soil moisture movement during transpiration are not large enough to appear in volumetric water content data. Fluid flow and electrokinetic coupling (i.e. electrical potential distribution) were simulated using COMSOL Multiphysics to explore the system controls on field data. The coupled model, which included a root-water uptake term, reproduced components of both the long-term and diel variations in SP measurements, thus indicating that SP has potential to provide spatially and temporally dense measurements of transpiration-induced changes in water flow. This manuscript presents the first SP measurements focusing on the movement of soil moisture in response to tree transpiration., Comment: 18 pages, 7 figures

Among the different contributions generating self-potential, the streaming potential is of particular interest in hydrogeology for its sensitivity to water flow. Estimating water flux in porous media using streaming potential data relies on our capacity to understand, model, and upscale the electrokinetic coupling at the mineral-solution interface. Different approaches have been proposed to predict streaming potential generation in porous media. One of these approaches is the flux averaging which is based on determining the excess charge which is effectively dragged in the medium by water flow. In this study, we develop a physically-based analytical model to predict the effective excess charge in saturated porous media using a flux-averaging approach in a bundle of capillary tubes with a fractal pore-size distribution. The proposed model allows the determination of the effective excess charge as a function of pore water ionic concentration and hydrogeological parameters like porosity, permeability and tortuosity. The new model has been successfully tested against different set of experimental data from the literature. One of the main findings of this study is the mechanistic explanation to the empirical dependence between the effective excess charge and the permeability that has been found by several researchers. The proposed model also highlights the link to other lithological properties and it is able to reproduce the evolution of effective excess charge with electrolyte concentrations., Comment: 22 pages, 7 figures

The self-potential (SP) method is a passive geophysical method that relies on the measurement of naturally occurring electrical field. One of the contributions to the SP signal is the streaming potential, which is of particular interest in hydrogeophysics as it is directly related to both the water flow and porous medium properties. The streaming current is generated by the relative displacement of an excess of electrical charges located in the electrical double layer surrounding the minerals of the porous media. In this study, we develop a physically based analytical model to estimate the effective excess charge density dragged by the water flow under partially saturated conditions. The proposed model is based on the assumption that the porous media can be represented by a bundle of tortuous capillary tubes with a fractal pore size distribution. The excess charge that is effectively dragged by the water flow is estimated using a flux averaging approach. Under these hypotheses, this new model describes the effective excess charge density as a function of saturation and relative permeability while also depending on the chemical and interface properties, and on petrophysical parameters of the media. The expression of the model has an analytical single closed-form which is consistent with a previous model developed from a different approach. The performance of the proposed model is then tested against previous models and different sets of laboratory and field data from the literature. The predictions of the proposed model fits fairly well the experimental data and shows improvements to estimate the magnitude of the effective excess charge density over the previous models. This new model proposes a simple and efficient way to model the streaming current generation for partially saturated porous media., Comment: 30 pages, 10 figures, 7 tables

In this paper we present a constitutive model to describe unsaturated flow that considers the hysteresis phenomena. This constitutive model provides simple mathematical expressions for both saturation and hydraulic conductivity curves, and a relationship between permeability and porosity. The model is based on the assumption that the porous media can be represented by a bundle of capillary tubes with throats or "ink-bottles" and a fractal pore size distribution. Under these hypotheses, hysteretic curves are obtained for saturation and relative hydraulic conductivity in terms of pressure head. However, a non-hysteretic relationship is obtained when relative hydraulic conductivity is expressed as a function of saturation. The proposed relationship between permeability and porosity is similar to the well-known Kozeny-Carman equation but depends on the fractal dimension. The performance of the constitutive model is tested against different sets of experimental data and previous models. In all of the cases the proposed expressions fit fairly well the experimental data and predicts values of permeability and hydraulic conductivity better than others models., Comment: 15 pages, 5 figures

We implement a procedure to simulate the drainage and imbibition in random, two-dimensional, square networks. We compute the resistivity index, the relative permeability, and the characteristic lengths of a correlated network at various saturation states, under the assumption that the surface conductivity can be neglected. These parameters exhibit a hysteretic behaviour. Then, we calculate the Spectral Induced Polarization (SIP) response of the medium, under the assumption that the electrical impedance of each tube follows a local Warburg conductivity model, with identical DC conductivity and chargeability for all the tubes. We evidence that the shape of the SIP spectra depends on the saturation state. The analysis of the evolution of the macroscopic Cole-Cole parameters of the spectra in function of the saturation also behaves hysteretically, except for the Cole-Cole exponent. We also observe a power-law relationship between the macroscopic DC conductivity and time constant and the relative permeability. We also show that the frequency peak of the phase spectra is directly related to the characteristic length and to the relative permeability, underlining the potential interest of SIP measurements for the estimation of the permeability of unsaturated media., Comment: 47 pages, 19 figures, 2 tables

Zeta potential is a physicochemical parameter of particular importance in describing the surface electrical properties of charged porous media. However, the zeta potential of calcite is still poorly known because of the difficulty to interpret streaming potential experiments. The Helmholtz-Smoluchowski (HS) equation is widely used to estimate the apparent zeta potential from these experiments. However, this equation neglects the influence of surface conductivity on streaming potential. We present streaming potential and electrical conductivity measurements on a calcite powder in contact with an aqueous NaCl electrolyte. Our streaming potential model corrects the apparent zeta potential of calcite by accounting for the influence of surface conductivity and flow regime. We show that the HS equation seriously underestimates the zeta potential of calcite, particularly when the electrolyte is diluted (ionic strength < 0.01 M) because of calcite surface conductivity. The basic Stern model successfully predicted the corrected zeta potential by assuming that the zeta potential is located at the outer Helmholtz plane, i.e. without considering a stagnant diffuse layer at the calcite-water interface. The surface conductivity of calcite crystals was inferred from electrical conductivity measurements and computed using our basic Stern model. Surface conductivity was also successfully predicted by our surface complexation model., Comment: 14 pages, 9 figures

We propose a method for filling gaps and removing interferences in time series for applications involving continuous monitoring of environmental variables. The approach is non-parametric and based on an iterative pattern-matching between the affected and the valid parts of the time series. It considers several variables jointly in the pattern matching process and allows preserving linear or non-linear dependences between variables. The uncertainty in the reconstructed time series is quantified through multiple realizations. The method is tested on self-potential data that are affected by strong interferences as well as data gaps, and the results show that our approach allows reproducing the spectral features of the original signal. Even in the presence of intense signal perturbations, it significantly improves the signal and corrects bias introduced by asymmetrical interferences. Potential applications are wide-ranging, including geophysics, meteorology and hydrology., Comment: 6 pages, 4 figures

The presence of mesoscopic heterogeneities in fluid-saturated porous rocks can produce measurable seismoelectric signals due to wave-induced fluid flow between regions of differing compressibility. The dependence of these signals on the petrophysical and structural characteristics of the probed rock mass remains largely unexplored. In this work, we derive an analytical solution to describe the seismoelectric response of a rock sample, containing a horizontal layer at its center, that is subjected to an oscillatory compressibility test. We then adapt this general solution to compute the seismoelectric signature of a particular case related to a sample that is permeated by a horizontal fracture located at its center. Analyses of the general and particular solutions are performed to study the impact of different petrophysical and structural parameters on the seismoelectric response. We find that the amplitude of the seismoelectric signal is directly proportional to the applied stress, to the Skempton coefficient contrast between the host rock and the layer, and to a weighted average of the effective excess charge of the two materials. Our results also demonstrate that the frequency at which the maximum electrical potential amplitude prevails does not depend on the applied stress or the Skempton coefficient contrast. In presence of strong permeability variations, this frequency is rather controlled by the permeability and thickness of the less permeable material. The results of this study thus indicate that seismoelectric measurements can potentially be used to estimate key mechanical and hydraulic rock properties of mesoscopic heterogeneities, such as compressibility, permeability, and fracture compliance., Comment: 14 pages, 8 figures

The self-potential (SP) method is sensitive to water fluxes in saturated and partially saturated porous media, such as those associated with rainwater infiltration and groundwater recharge. We present a field-based study at the Voulund agricultural test site, Denmark, that is, to the best of our knowledge, the first to focus on the vertical self-potential distribution prior to and during a saline tracer test. A coupled hydrogeophysical modeling framework is used to simulate the SP response to precipitation and saline tracer infiltration. A layered hydrological model is first obtained by inverting dielectric and matric potential data. The resulting model that compares favorably with electrical resistance tomography models is subsequently used to predict the SP response. The electrokinetic contribution (caused by water fluxes in a charged porous soil) is modeled by an effective excess charge approach that considers both water saturation and pore water salinity. Our results suggest that the effective excess charge evolution prior to the tracer injection is better described by a recent flux-averaged model based on soil water retention functions than by a previously proposed volume-averaging model. This is the first time that raw vertically distributed SP measurements have been explained by a physically based model. The electrokinetic contribution cannot alone reproduce the SP data during the tracer test and an electro-diffusive contribution (caused by concentration gradients) is needed. The predicted amplitude of this contribution is too small to perfectly explain the data, but the shape is in accordance with the field data. This discrepancy is attributed to imperfect descriptions of electro-diffusive phenomena in partially saturated soils, unaccounted soil heterogeneity, and discrepancies between the measured and predicted electrical conductivities in the tracer infiltration area., Comment: 45 pages, 9 figures, 2 tables

Water movement in unsaturated soils gives rise to measurable electrical potential differences that are related to the flow direction and volumetric fluxes, as well as to the soil properties themselves. Laboratory and field data suggest that these so-called streaming potentials may be several orders of magnitudes larger than theoretical predictions that only consider the influence of the relative permeability and electrical conductivity on the self potential (SP) data. Recent work has partly improved predictions by considering how the volumetric excess charge in the pore space scales with the inverse of water saturation. We present a new theoretical approach that uses the flux-averaged excess charge, not the volumetric excess charge, to predict streaming potentials. We present relationships for how this effective excess charge varies with water saturation for typical soil properties using either the water retention or the relative permeability function. We find large differences between soil types and the predictions based on the relative permeability function display the best agreement with field data. The new relationships better explain laboratory data than previous work and allow us to predict the recorded magnitudes of the streaming potentials following a rainfall event in sandy loam, whereas previous models predict three orders of magnitude too small values. We suggest that the strong signals in unsaturated media can be used to gain information about fluxes (including very small ones related to film flow), but also to constrain the relative permeability function, the water retention curve, and the relative electrical conductivity function.