Thermal conductivity of unsaturated clay-rocks

The parameters used to describe the electrical conductivity of a porous material can be used to describe also its thermal conductivity. A new relationship is developed to connect the thermal conductivity of an unsaturated porous material to the thermal conductivity of the different phases of the composite, and two electrical parameters called the first and second Archie's exponents. A good agreement is obtained between the new model and thermal conductivity measurements performed using packs of glass beads and core samples of the Callovo-Oxfordian clay-rocks at different saturations of the water phase. We showed that the three model parameters optimised to fit the new model against experimental data (namely the thermal conductivity of the solid phase and the two Archie's exponents) are consistent with independent estimates. We also observed that the anisotropy of the effective thermal conductivity of the Callovo-Oxfordian clay-rock was mainly due to the anisotropy of the thermal conductivity of the solid phase

Impact of water saturation on seismoelectric transfer functions: a laboratory study of coseismic phenomenon

Seismic waves propagating in a porous medium, under favourable conditions, generate measurable electromagnetic fields due to electrokinetic effects. It has been proposed, following experimental and numerical studies, that these so-called ‘seismoelectromagnetic' couplings depend on pore fluid properties. The theoretical frame describing these phenomena are based on the original Biot's theory, assuming that pores are fluid-filled. We study here the impact of a partially saturated medium on amplitudes of those seismoelectric couplings by comparing experimental data to an effective fluid model. We have built a 1-m-length-scale experiment designed for imbibition and drainage of an homogeneous silica sand; the experimental set-up includes a seismic source, accelerometers, electric dipoles and capacitance probes in order to monitor seismic and seismoelectric fields during water saturation. Apparent velocities and frequency spectra (in the kiloHertz range) are derived from seismic and electrical measurements during experiments in varying saturation conditions. Amplitudes of seismic and seismoelectric waves and their ratios (i.e. transfer functions) are discussed using a spectral analysis performed by continuous wavelet transform. The experiments reveal that amplitude ratios of seismic to coseismic electric signals remain rather constant as a function of the water saturation in the Sw=[0.2-0.9] range, consistently with theoretically predicted transfer function

Streaming potential modeling in fractured rock: Insights into the identification of hydraulically active fractures

Numerous field experiments suggest that the self-potential (SP) geophysical method may allow for the detection of hydraulically active fractures and provide information about fracture properties. However, a lack of suitable numerical tools for modeling streaming potentials in fractured media prevents quantitative interpretation and limits our understanding of how the SP method can be used in this regard. To address this issue, we present a highly efficient two-dimensional discrete-dual-porosity approach for solving the fluid flow and associated self-potential problems in fractured rock. Our approach is specifically designed for complex fracture networks that cannot be investigated using standard numerical methods. We then simulate SP signals associated with pumping conditions for a number of examples to show that (i) accounting for matrix fluid flow is essential for accurate SP modeling and (ii) the sensitivity of SP to hydraulically active fractures is intimately linked with fracture-matrix fluid interactions. This implies that fractures associated with strong SP amplitudes are likely to be hydraulically conductive, attracting fluid flow from the surrounding matrix.Comment: 8 pages, 2 figure

Un modelo fractal para estimar la permeabilidad a partir de la porosidad

En este trabajo se presenta un modelo fractal para estimar la permeabilidad de una roca en función de su porosidad. El modelo asume que el medio poroso se puede representar mediante un conjunto de tubos capilares constrictivos con una distribución fractal de tamaño de poro. La relación obtenida es una ley de potencias que depende de tres parámetros independientes con significado físico o geométrico. Cabe resaltar que el exponente de la relación propuesta depende de la dimensión fractal y resulta siempre mayor a 3. En el caso límite de un exponente igual a 3, la relación propuesta es similar al modelo semi-empírico de Kozeny-Carman (KC). El modelo fractal se validó con datos experimentales de la literatura obteniéndose un buen ajuste para distintas texturas de suelo. La comparación con datos experimentales muestra que la relación propuesta predice mejor los valores de permeabilidad que la ecuación de KC para todo el rango de magnitudes.Eje: Ciencias Hidrológicas y Criósfera.Facultad de Ciencias Astronómicas y Geofísica

An analytical study of seismoelectric signals produced by 1-D mesoscopic heterogeneities

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 figure

Toward characterization of organic matter rich in aromatic compounds by spectral-induced polarization: preliminary investigation and perspectives

Spectral induced polarization (SIP) has been suggested as a non-invasive and cost-effective tool to detect and monitor aromatic rich organic matter such as biochar. In our study, we show that SIP can track biochar concentration up to 10 per cent (wt.) in a soil with a clay content of 20 per cent. Assessment of changes in the concentration of biochar was conducted according to double Pelton parameters and the maximal phase determined at 11.7 Hz, a frequency at which a polarization peak is observed in the presence of biochar. All SIP-derived parameters were correlated with the biochar content, with the exception of the relaxation time of the polarization peak occurring at 11.7 Hz, which was related to soil water saturation in previous investigations. Among studied parameters, the phase value that we measured at 11.7 Hz may therefore consist in a simple and reliable methodology to evaluate the biochar content on SIP in our experiment. Several steps are still necessary before a widespread field application notably by considering how modifications in the chemistry of biochar with time can interact with biochar concentration and water saturation to modify polarization processes shaping SIP curves. Beyond the scope of tracking changes in the content of highly aromatic OM—such as biochar—in soils, this study suggests that the degree of aromaticity of OM can play a key role in the SIP response paving the way for wider use of SIP in soil science

Self-Potential as a Predictor of Seawater Intrusion in Coastal Groundwater Boreholes

This work was supported by the Natural Environment Research Council in the UK, as part of the Science and Solutions for a Changing Planet Doctor Training Partnership, run by the Grantham Institute for Climate Change at Imperial College London. We thank Southern Water for access to the boreholes at Saltdean and Balsdean. We thank Southern Water and Atkins Global for funding the installation of the equipment. We also thank Dr Amadi Ijioma for providing a prototype of the electrodynamic modelling code in MATLAB, which has since been adapted for use in a coastal chalk aquifer. Three anonymous reviewers are thanked for their comments, which greatly helped to improve the manuscript. The data used in this paper are in the tables, figures and cited information. The authors have no conflicts of interest to declare.Peer reviewedPublisher PDFPublisher PD

Diffusion of ions in unsaturated porous materials

International audienceIn a salinity gradient, the diffusion of ions through the connected porosity of a porous and charged material is influenced by the charged nature of the interface between the pore water and the solid. This influence is exerted through the generation of a macroscopic electrical field termed the diffusion or membrane potential. This electrical field depends on the excess of counterions located in the pore space counterbalancing the charge density of the surface of the solid. In unsaturated porous materials, we have to consider (1) the effect of the charged nature of the air/water interface, (2) the increase of the counterion density as the counterions are packed in a smaller volume when the saturation of the nonwetting phase (air) increases, and (3) the influence of the water saturation upon the tortuosity of the water phase. The volume average of the Nernst–Planck equation is used to determine the constitutive equations for the coupled diffusion flux and current density of a multicomponent electrolyte in unsaturated conditions. We assume that water is the wetting phase for the solid phase. We neglect the electro-osmotic flow in the coupled constitutive equations and the deformation of the medium (the medium is assumed to be both isotropic and rigid). This model explains well the observed tendency of strong decreases of the apparent diffusion coefficient of ions with the decrease of the saturation of the water phase under steady-state conditions. This decrease is mainly due to the influence of the saturation upon the tortuosity of the water phase