scholarly journals Fluid pressure diffusion in fractured media: the role played by the geometry of real fractures

Author(s):  
Simón Lissa ◽  
Nicolás D. Barbosa ◽  
Beatriz Quintal
2021 ◽  
Author(s):  
Samuel Chapman ◽  
Jan V. M. Borgomano ◽  
Beatriz Quintal ◽  
Sally M. Benson ◽  
Jerome Fortin

<p>Monitoring of the subsurface with seismic methods can be improved by better understanding the attenuation of seismic waves due to fluid pressure diffusion (FPD). In porous rocks saturated with multiple fluid phases the attenuation of seismic waves by FPD is sensitive to the mesoscopic scale distribution of the respective fluids. The relationship between fluid distribution and seismic wave attenuation could be used, for example, to assess the effectiveness of residual trapping of carbon dioxide (CO2) in the subsurface. Determining such relationships requires validating models of FPD with accurate laboratory measurements of seismic wave attenuation and modulus dispersion over a broad frequency range, and, in addition, characterising the fluid distribution during experiments. To address this challenge, experiments were performed on a Berea sandstone sample in which the exsolution of CO2 from water in the pore space of the sample was induced by a reduction in pore pressure. The fluid distribution was determined with X-ray computed tomography (CT) in a first set of experiments. The CO2 exosolved predominantly near the outlet, resulting in a heterogeneous fluid distribution along the sample length. In a second set of experiments, at similar pressure and temperature conditions, the forced oscillation method was used to measure the attenuation and modulus dispersion in the partially saturated sample over a broad frequency range (0.1 - 1000 Hz). Significant P-wave attenuation and dispersion was observed, while S-wave attenuation and dispersion were negligible. These observations suggest that the dominant mechanism of attenuation and dispersion was FPD. The attenuation and dispersion by FPD was subsequently modelled by solving Biot’s quasi-static equations of poroelasticity with the finite element method. The fluid saturation distribution determined from the X-ray CT was used in combination with a Reuss average to define a single phase effective fluid bulk modulus. The numerical solutions agree well with the attenuation and modulus dispersion measured in the laboratory, supporting the interpretation that attenuation and dispersion was due to FPD occurring in the heterogenous distribution of the coexisting fluids. The numerical simulations have the advantage that the models can easily be improved by including sub-core scale porosity and permeability distributions, which can also be determined using X-ray CT. In the future this could allow for conducting experiments on heterogenous samples.</p>


2012 ◽  
Vol 117 (B5) ◽  
pp. n/a-n/a ◽  
Author(s):  
Luca Malagnini ◽  
Francesco Pio Lucente ◽  
Pasquale De Gori ◽  
Aybige Akinci ◽  
Irene Munafo'

2018 ◽  
Vol 115 (8) ◽  
pp. E1720-E1729 ◽  
Author(s):  
Maria Kozłowska ◽  
Michael R. Brudzinski ◽  
Paul Friberg ◽  
Robert J. Skoumal ◽  
Nicholas D. Baxter ◽  
...  

Understanding the causes of human-induced earthquakes is paramount to reducing societal risk. We investigated five cases of seismicity associated with hydraulic fracturing (HF) in Ohio since 2013 that, because of their isolation from other injection activities, provide an ideal setting for studying the relations between high-pressure injection and earthquakes. Our analysis revealed two distinct groups: (i) deeper earthquakes in the Precambrian basement, with larger magnitudes (M > 2), b-values < 1, and many post–shut-in earthquakes, versus (ii) shallower earthquakes in Paleozoic rocks ∼400 m below HF, with smaller magnitudes (M < 1), b-values > 1.5, and few post–shut-in earthquakes. Based on geologic history, laboratory experiments, and fault modeling, we interpret the deep seismicity as slip on more mature faults in older crystalline rocks and the shallow seismicity as slip on immature faults in younger sedimentary rocks. This suggests that HF inducing deeper seismicity may pose higher seismic hazards. Wells inducing deeper seismicity produced more water than wells with shallow seismicity, indicating more extensive hydrologic connections outside the target formation, consistent with pore pressure diffusion influencing seismicity. However, for both groups, the 2 to 3 h between onset of HF and seismicity is too short for typical fluid pressure diffusion rates across distances of ∼1 km and argues for poroelastic stress transfer also having a primary influence on seismicity.


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.


2016 ◽  
Vol 140 (4) ◽  
pp. 2554-2570 ◽  
Author(s):  
Nicolás D. Barbosa ◽  
J. Germán Rubino ◽  
Eva Caspari ◽  
Marco Milani ◽  
Klaus Holliger

2020 ◽  
Vol 8 (2) ◽  
pp. SG21-SG31
Author(s):  
Yadong Zhang ◽  
Reza Rezaee ◽  
Tobias M. Müller ◽  
Guangjie Zheng ◽  
Jimmy X. Li ◽  
...  

We have predicted the flow permeability and its spatial distribution for the Longmaxi shale gas reservoir using microseismicity induced during hydraulic-fracturing stimulation. In the time-of-occurrence versus distance-from-injector plot, we find that microseismic points exhibit a parabolic envelope, which we interpret as a triggering front. This reveals that fluid pressure diffusion is at least one of the underlying mechanisms of microseismicity generation. We derive the large-scale equivalent diffusivity from the triggering front plot and thereafter obtain a 3D diffusivity map of the heterogeneous reservoir by solving an eikonal-like equation suggested previously. During this process, we apply kriging interpolation to increase the density of sparsely distributed microseismic points. The resulting diffusivity ranges between 1.0 and [Formula: see text] with the peak probability attained at [Formula: see text], which is consistent with the estimate we obtain from the triggering front analysis. We transform the diffusivity map into a permeability map using three different theories of fluid pressure diffusion in porous media: the seismicity-based reservoir characterization method (SBRC) based on Biot’s theory of poroelasticity, the quasirigid medium approximation (QRMA), and the deformable medium approximation (DMA) based on the de la Cruz-Spanos theory. The permeability according to QRMA is slightly higher than that from SBRC, yet we observe no significant difference. However, these estimates are by one order of magnitude higher compared with the permeability estimate from DMA. Furthermore, the permeability from all three theories is much higher than that from previously reported core sample measurements. We interpret this as the difference between large-scale equivalent and matrix permeability and therefore lend weight to the hypothesis that there exist highly conducting fluid pathways, such as natural fractures.


2020 ◽  
Vol 91 (6) ◽  
pp. 3171-3181 ◽  
Author(s):  
Maomao Wang ◽  
Hongfeng Yang ◽  
Lihua Fang ◽  
Libo Han ◽  
Dong Jia ◽  
...  

Abstract Human activity-induced earthquakes are emerging as a global issue, and revealing its underlying mechanisms is essential for earthquake hazard mitigation and energy development. We investigated the relationship between the seismotectonic model and seismic sequences from moderate Mw 4.3 and Mw 5.2 earthquakes that occurred in February and September 2019, respectively, in the Weiyuan anticline of Sichuan basin, China. We found that the Mw 5.2 earthquake ruptured a back thrust of structural wedges and released most aftershocks near the wedge tip. However, the two foreshocks of the Mw 4.3 earthquake sequence occurred in hydrofractured Silurian shale at depth of 2.5–3 km, and the mainshock ruptured the overlying oblique tear fault at a depth of ∼1  km. Hydraulic fracturing in the sedimentary cover of this block may induce earthquakes through fluid pressure diffusion in the Silurian shale and through poroelastic effects on back thrusts within structural wedges, respectively. We assessed the hazard potential of four seismic sources in the Weiyuan block and suggest it is critical to conduct a coupled flow-geomechanics assessment and management on induced seismicity and related cascading effects in the densely inhabited and seismically active Sichuan basin.


Sign in / Sign up

Export Citation Format

Share Document