Hydro-mechanical coupling in porous rocks: hidden dependences to the microstructure?

2020 ◽  
Vol 224 (2) ◽  
pp. 973-984
Author(s):  
Lucas Pimienta ◽  
Beatriz Quintal ◽  
Eva Caspari

SUMMARY While hydro-mechanical coupling in rocks is generally well understood, exotic rock poroelastic responses—such as unexpected dependence to fluid diffusion time and low skeleton moduli—have been reported. Hydro-mechanical coupling, or poroelasticity, explains how fluid-saturated rocks respond to either confining or fluid pressure variations. This coupling is usually inferred from the apparent mechanical and hydraulic properties: mechanical properties determine the strain level experienced by the rock when submitted to pressures at a timescale when fluid pressure is equilibrated, which is in turn ruled by hydraulic properties. However, the coupling between properties might not always be straightforward, particularly for rocks in which two distinct families of pore types coexist: spherical pores and cracks. Comparing it with reported laboratory data sets on pressure-dependent hydraulic and elastic properties in sandstones of different porosity confirms that the well-known simple concepts of networks in parallel apply and yield opposite dependencies for the two properties on the two pore families. Because hydro-mechanical coupling implies that the two properties—that depend differently on the pore families—will be interdependent, we further apply the same concept of parallel network. It yields that, although under apparent drained conditions, typical poroelasticity experiments could underestimate the rock compressibility ${C_{\mathrm{ bp}}}$, measured as a response to fluid pressure variation, and underestimate the related skeleton (or unjacketed) bulk modulus ${K_\mathrm{ s}} = \ 1/{C_\mathrm{ s}}$.

Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1173-1181 ◽  
Author(s):  
S. Richard Taylor ◽  
Rosemary J. Knight

Our new method incorporates fluid pressure communication into inclusion‐based models of elastic wave velocities in porous rocks by defining effective elastic moduli for fluid‐filled inclusions. We illustrate this approach with two models: (1) flow between nearest‐neighbor pairs of inclusions and (2) flow through a network of inclusions that communicates fluid pressure throughout a rock sample. In both models, we assume that pore pressure gradients induce laminar flow through narrow ducts, and we give expressions for the effective bulk moduli of inclusions. We compute P‐wave velocities and attenuation in a model sandstone and illustrate that the dependence on frequency and water‐saturation agrees qualitatively with laboratory data. We consider levels of water saturation from 0 to 100% and all wavelengths much larger than the scale of material heterogeneity, obtaining near‐exact agreement with Gassmann theory at low frequencies and exact agreement with inclusion‐based models at high frequencies.


Minerals ◽  
2020 ◽  
Vol 11 (1) ◽  
pp. 33
Author(s):  
Valérie Laperche ◽  
Bruno Lemière

Portable X-ray fluorescence spectroscopy is now widely used in almost any field of geoscience. Handheld XRF analysers are easy to use, and results are available in almost real time anywhere. However, the results do not always match laboratory analyses, and this may deter users. Rather than analytical issues, the bias often results from sample preparation differences. Instrument setup and analysis conditions need to be fully understood to avoid reporting erroneous results. The technique’s limitations must be kept in mind. We describe a number of issues and potential pitfalls observed from our experience and described in the literature. This includes the analytical mode and parameters; protective films; sample geometry and density, especially for light elements; analytical interferences between elements; physical effects of the matrix and sample condition, and more. Nevertheless, portable X-ray fluorescence spectroscopy (pXRF) results gathered with sufficient care by experienced users are both precise and reliable, if not fully accurate, and they can constitute robust data sets. Rather than being a substitute for laboratory analyses, pXRF measurements are a valuable complement to those. pXRF improves the quality and relevance of laboratory data sets.


Minerals ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 162
Author(s):  
A.A. Jameei ◽  
S. Pietruszczak

This paper provides a mathematical description of hydromechanical coupling associated with propagation of localized damage. The framework incorporates an embedded discontinuity approach and addresses the assessment of both hydraulic and mechanical properties in the region intercepted by a fracture. Within this approach, an internal length scale parameter is explicitly employed in the definition of equivalent permeability as well as the tangential stiffness operators. The effect of the progressive evolution of damage on the hydro-mechanical coupling is examined and an evolution law is derived governing the variation of equivalent permeability with the continuing deformation. The framework is verified by a numerical study involving 3D simulation of an axial splitting test carried out on a saturated sample under displacement and fluid pressure-controlled conditions. The finite element analysis incorporates the Polynomial-Pressure-Projection (PPP) stabilization technique and a fully implicit time integration scheme.


2018 ◽  
Vol 11 (11) ◽  
pp. 6203-6230 ◽  
Author(s):  
Simon Ruske ◽  
David O. Topping ◽  
Virginia E. Foot ◽  
Andrew P. Morse ◽  
Martin W. Gallagher

Abstract. Primary biological aerosol including bacteria, fungal spores and pollen have important implications for public health and the environment. Such particles may have different concentrations of chemical fluorophores and will respond differently in the presence of ultraviolet light, potentially allowing for different types of biological aerosol to be discriminated. Development of ultraviolet light induced fluorescence (UV-LIF) instruments such as the Wideband Integrated Bioaerosol Sensor (WIBS) has allowed for size, morphology and fluorescence measurements to be collected in real-time. However, it is unclear without studying instrument responses in the laboratory, the extent to which different types of particles can be discriminated. Collection of laboratory data is vital to validate any approach used to analyse data and ensure that the data available is utilized as effectively as possible. In this paper a variety of methodologies are tested on a range of particles collected in the laboratory. Hierarchical agglomerative clustering (HAC) has been previously applied to UV-LIF data in a number of studies and is tested alongside other algorithms that could be used to solve the classification problem: Density Based Spectral Clustering and Noise (DBSCAN), k-means and gradient boosting. Whilst HAC was able to effectively discriminate between reference narrow-size distribution PSL particles, yielding a classification error of only 1.8 %, similar results were not obtained when testing on laboratory generated aerosol where the classification error was found to be between 11.5 % and 24.2 %. Furthermore, there is a large uncertainty in this approach in terms of the data preparation and the cluster index used, and we were unable to attain consistent results across the different sets of laboratory generated aerosol tested. The lowest classification errors were obtained using gradient boosting, where the misclassification rate was between 4.38 % and 5.42 %. The largest contribution to the error, in the case of the higher misclassification rate, was the pollen samples where 28.5 % of the samples were incorrectly classified as fungal spores. The technique was robust to changes in data preparation provided a fluorescent threshold was applied to the data. In the event that laboratory training data are unavailable, DBSCAN was found to be a potential alternative to HAC. In the case of one of the data sets where 22.9 % of the data were left unclassified we were able to produce three distinct clusters obtaining a classification error of only 1.42 % on the classified data. These results could not be replicated for the other data set where 26.8 % of the data were not classified and a classification error of 13.8 % was obtained. This method, like HAC, also appeared to be heavily dependent on data preparation, requiring a different selection of parameters depending on the preparation used. Further analysis will also be required to confirm our selection of the parameters when using this method on ambient data. There is a clear need for the collection of additional laboratory generated aerosol to improve interpretation of current databases and to aid in the analysis of data collected from an ambient environment. New instruments with a greater resolution are likely to improve on current discrimination between pollen, bacteria and fungal spores and even between different species, however the need for extensive laboratory data sets will grow as a result.


Geophysics ◽  
1986 ◽  
Vol 51 (3) ◽  
pp. 780-787 ◽  
Author(s):  
Kai Hsu ◽  
Arthur B. Baggeroer

Modern digital sonic tools can record full waveforms using an array of receivers. The recorded waveforms are extremely complicated because wave components overlap in time. Conventional beamforming approaches, such as semblance processing, while robust, sometimes do not resolve the interfering wave components propagating at similar speeds, such as multiple compressional arrivals due to the formation alteration surrounding the borehole. Here the maximum‐likelihood method (MLM), a high‐resolution array processing algorithm, is modified and applied to process borehole array sonic data. Extensive modifications of the original MLM algorithm were necessary because of the transient character of the sonic data and its effect upon the spectral covariance matrix. We applied MLM to several array sonic data sets, including laboratory data, synthetic waveforms, and field data taken by a Schlumberger array sonic tool. MLM’s slowness resolution is demonstrated in resolving a secondary compressional arrival from the primary compressional arrival in an altered formation, and the formation compressional arrival in the presence of a stronger casing arrival in an unbonded cased hole. In comparison with the semblance processing results, the MLM results clearly show a better slowness resolution. However, in the case of a weak formation arrival, the semblance processing tends to enhance and resolve the weak arrival by the semblance normalization procedure, while the MLM, designed to estimate the signal strength, does not. The heavy computational requirement (mainly, many matrix inversions) in the MLM makes it much slower than semblance processing, which may prohibit implementation of the MLM algorithm in a real‐time environment.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA167-WA183 ◽  
Author(s):  
Patrick N. J. Rasolofosaon

The great diversity of the microstructures of rocks impedes the use of a universal rock physics model with idealized geometry to correctly describe the mechanical behavior of rocks. In this quest for universality, by ignoring the detailed description of the causes of the observed phenomenon and only focusing on the empirical relation between the cause (applied stress) and the effect (resulting strain), phenomenological models such as the linear elastic Hooke’s law roughly describe the mechanical behavior of rocks of contrasted microstructures. However, in detail, numerous laboratory experiments covering broad frequency and strain ranges (both typically more than eight orders of magnitude) on various types of rocks have also shown deviations from Hooke’s law due to anisotropy, frequency dependence, nonlinearity, possibly with the presence of hysteresis, and poroelasticity. A phenomenological model has been recently proposed that synthesizes all these behaviors in a single model, but unfortunately does not integrate the porous nature of rocks. The new model is based on a reformulation in modified spectral decomposition of the previous work using the 7D poroelastic tensor linking the dynamic parameters (i.e., the six stress components and fluid pressure) and the kinematic parameters (i.e., the six strain components and the local increase of fluid content ζ). In addition to the elastic hysteresis of the stress-strain curves, the model also predicts the existence of a second hysteresis, or hydraulic hysteresis, of the curve fluid pressure p versus fluid content ζ, qualitatively similar to the first one. Indeed, the elastic hysteresis is due to the opening and the closure of some compliant pores at different stress levels. These pores represent possible access radii for the saturating fluid; the hysteresis in the geometry of the porous network also induces the hydraulic hysteresis in the p-ζ curves.


2020 ◽  
Vol 222 (1) ◽  
pp. 715-733
Author(s):  
Gabriel A Castromán ◽  
Nicolás D Barbosa ◽  
J Germán Rubino ◽  
Fabio I Zyserman ◽  
Klaus Holliger

SUMMARY The presence of sets of open fractures is common in most reservoirs, and they exert important controls on the reservoir permeability as fractures act as preferential pathways for fluid flow. Therefore, the correct characterization of fracture sets in fluid-saturated rocks is of great practical importance. In this context, the inversion of fracture characteristics from seismic data is promising since their signatures are sensitive to a wide range of pertinent fracture parameters, such as density, orientation and fluid infill. The most commonly used inversion schemes are based on the classical linear slip theory (LST), in which the effects of the fractures are represented by a real-valued diagonal excess compliance matrix. To account for the effects of wave-induced fluid pressure diffusion (FPD) between fractures and their embedding background, several authors have shown that this matrix should be complex-valued and frequency-dependent. However, these approaches neglect the effects of FPD on the coupling between orthogonal deformations of the rock. With this motivation, we considered a fracture model based on a sequence of alternating poroelastic layers of finite thickness representing the background and the fractures, and derived analytical expressions for the corresponding excess compliance matrix. We evaluated this matrix for a wide range of background parameters to quantify the magnitude of its coefficients not accounted for by the classical LST and to determine how they are affected by FPD. We estimated the relative errors in the computation of anisotropic seismic velocity and attenuation associated with the LST approach. Our analysis showed that, in some cases, considering the simplified excess compliance matrix may lead to an incorrect representation of the anisotropic response of the probed fractured rock.


2020 ◽  
Author(s):  
Yury Alkhimenkov ◽  
Lyudmila Khakimova ◽  
Ludovic Raess ◽  
Beatriz Quintal ◽  
Yury Podladchikov

<p> <span>Elastodynamic hydro-mechanical coupling based on Biot’s theory describes an upscaling of the fluid-solid deformation at a porous scale. Examples of applications of this theory are near surface geophysics, CO2 monitoring, induced seismicity, etc. The dynamic response of a coupled hydro-mechanical system can produce fast and slow compressional waves and shear waves. In many earth materials, a propagating slow wave degenerates into a slow diffusion mode on orders of magnitude larger time scales compared to wave propagation. In the present work, we propose a new approach to accelerate the numerical simulation of slow diffusion processes. We solve the coupled Biot elastodynamic hydro-mechanical equations for particle velocity and stress in the time domain using the finite volume method on a rectangular grid in three dimensions. The MPI-based multi-GPU code is implemented using CUDA-C programming language. We prescribe a fluid injection at the center of the model that generates a fast propagating wave and a significantly slower fluid-diffusion event. The fast wave is attenuated due to absorbing boundary conditions after what the slow fluid-diffusion process remains active. A Courant stability condition for the fast wave controls the time-step in the entire simulation, resulting in a suboptimal short time step for the diffusion process. Once fast waves are no longer present in the model domain, the hydro-mechanical coupling vanishes in the inertial terms allowing for an order of magnitude larger time steps. We accelerate the numerical simulation of slow diffusion processes using a pseudo-transient method that permits to capture the transition in time step restrictions. This latest development enables us to simulate quasi-static and dynamic responses of two-phase media. We present benchmarks confirming the numerical efficiency and accuracy of the novel approach. The further development of the code will capture inelastic physics starting from the dynamic events (earthquake modeling) to quasi-static faulting. </span></p>


Sign in / Sign up

Export Citation Format

Share Document