The effect of Biot coefficient and elastic moduli stress–pore pressure dependency on poroelastic response to fluid injection: laboratory experiments and geomechanical modeling

2020 ◽  
Vol 10 (5) ◽  
pp. 980-998 ◽  
Author(s):  
Samin Raziperchikolaee ◽  
Vivek Singh ◽  
Mark Kelley
Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-21
Author(s):  
Xi Chen ◽  
Wei Wang ◽  
Yajun Cao ◽  
Qizhi Zhu ◽  
Weiya Xu ◽  
...  

The study on hydromechanical coupling properties of rocks is of great importance for rock engineering. It is closely related to the stability analysis of structures in rocks under seepage condition. In this study, a series of conventional triaxial tests under drained condition and hydrostatic compression tests under drained or undrained condition on sandstones were conducted. Moreover, complex cyclic loading and unloading tests were also carried out. Based on the experimental results, the following conclusions were obtained. For conventional triaxial tests, the elastic modulus, peak strength, crack initiation stress, and expansion stress increase with increased confining pressure. Pore pressure weakened the effect of the confining pressure under drained condition, which led to a decline in rock mechanical properties. It appeared that cohesion was more sensitive to pore pressure than to the internal friction angle. For complex loading and unloading cyclic tests, in deviatoric stress loading and unloading cycles, elastic modulus increased obviously in first loading stage and increased slowly in next stages. In confining pressure loading and unloading cycles, the Biot coefficient decreased first and then increased, which indicates that damage has a great impact on the Biot coefficient.


2021 ◽  
Vol 40 (6) ◽  
pp. 413-417
Author(s):  
Chunfang Meng ◽  
Michael Fehler

As fluids are injected into a reservoir, the pore fluid pressure changes in space and time. These changes induce a mechanical response to the reservoir fractures, which in turn induces changes in stress and deformation to the surrounding rock. The changes in stress and associated deformation comprise the geomechanical response of the reservoir to the injection. This response can result in slip along faults and potentially the loss of fluid containment within a reservoir as a result of cap-rock failure. It is important to recognize that the slip along faults does not occur only due to the changes in pore pressure at the fault location; it can also be a response to poroelastic changes in stress located away from the region where pore pressure itself changes. Our goal here is to briefly describe some of the concepts of geomechanics and the coupled flow-geomechanical response of the reservoir to fluid injection. We will illustrate some of the concepts with modeling examples that help build our intuition for understanding and predicting possible responses of reservoirs to injection. It is essential to understand and apply these concepts to properly use geomechanical modeling to design geophysical acquisition geometries and to properly interpret the geophysical data acquired during fluid injection.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. KS105-KS118 ◽  
Author(s):  
Himanshu Barthwal ◽  
Mirko van der Baan

Hydraulic fracturing in low-permeability hydrocarbon reservoirs creates/reactivates a fracture network leading to microseismic events. We have developed a simplified model of the evolution of the microseismic cloud based on the opening of a planar fracture cavity and its effect on elastic stresses and pore pressure diffusion during fluid injection in hydraulic fracturing treatments. Using a material balance equation, we compute the crack tip propagation over time assuming that the hydraulic fracture is shaped as a single penny-shaped cavity. Results indicate that in low-permeability formations, the crack tip propagates much faster than the pore pressure diffusion front thereby triggering the microseismic events farthest from the injection domain at any given time during fluid injection. We use the crack tip propagation to explain the triggering front observed in distance versus time plots of published microseismic data examples from hydraulic fracturing treatments of low-permeability hydrocarbon reservoirs. We conclude that attributing the location of the microseismic triggering front purely to pore pressure diffusion from the injection point may lead to incorrect estimates of the hydraulic diffusivity by multiple orders of magnitude for low-permeability formations. Moreover, the opening of the fracture cavity creates stress shadow zones perpendicular to the principal fracture walls in which microseismic triggering due to the elastic stress perturbations is suppressed. Microseismic triggering in this stress shadow region may be attributed mainly to pore pressure diffusion. We use the width, instead of the longest size, of the microseismic cloud to obtain an enhanced diffusivity measure, which may be useful for subsequent production simulations.


2010 ◽  
Vol 13 (02) ◽  
pp. 265-274 ◽  
Author(s):  
Ashraf Al-Tahini ◽  
Younane Abousleiman

Summary In this study, we determine experimentally the effect of inherent and stress-induced anisotropy on stiffness components, elastic moduli, and Biot's pore-pressure coefficients (PPCs) for Lyons outcrop Colorado sandstone, which exhibits a clear transverse isotropic rock structure. Both dynamic and quasistatic methods were used under a nonhydrostatic state of stress to perform the measurements on dry core samples. Our assumption of apparent transverse anisotropy was confirmed initially with acoustic velocity measurements and at a later stage in the loading with experimental transverse anisotropic failure analysis. The objective of this study is to identify and isolate the effect of stress-induced anisotropy from the inherent transverse anisotropy on the measured stiffness components, elastic moduli, and Biot's PPCs. The effect of stress-induced anisotropy appears to have significant control on measured stiffness components, elastic moduli, and Biot's PPCs in comparison to the inherent-transverse-anisotropy effect. Our work shows that the stiffness components, Mij and thus the computed elastic moduli, are highly influenced by the stress-induced anisotropy, especially the off-diagonal stiffness components, M12 and M13, where the increase in their magnitudes from the dynamic measurements before failure is determined to be 100 and 81%, respectively. The difference in the magnitude between the axial and lateral Biot's PPCs in line with bedding planes and perpendicular to them is measured to be 24 and 16% from the quasistatic and dynamic methods, respectively; whereas, the effect of stress-induced anisotropy reduced the dynamic average magnitude of the Biot's PPCs along the bedding planes and transverse to these planes by 63% across a stress range of 145 MPa.


1982 ◽  
Vol 22 (03) ◽  
pp. 333-340 ◽  
Author(s):  
Norman R. Warpinski ◽  
James A. Clark ◽  
Richard A. Schmidt ◽  
Clarence W. Huddle

Abstract Laboratory experiments have been conducted to determine the effect of in-situ stress variations on hydraulic fracture containment. Fractures were initiated in layered rock samples with prescribed stress variations, and fracture growth characteristics were determined as a function of stress levels. Stress contrasts of 300 to 400 psi (2 to 3 MPa) were found sufficient to restrict fracture growth in laboratory samples of Nevada tuff and Tennessee and Nugget sandstones. The required stress level was found not to depend on mechanical rock properties. However, permeability and the resultant pore pressure effects were important. Tests conducted at biomaterial interfaces between Nugget and Tennessee sandstones show that the resultant stresses set up near the interface because of the applied overburden stress affect the fracture behavior in the same way as the applied confining stresses. These results provide a guideline for determining the in-situ stress contrast necessary to contain a fracture in a field treatment. Introduction An under-standing of the factors that influence and control hydraulic fracture containment is important for the successful use of hydraulic fracturing technology in the enhanced production of natural gas from tight reservoirs. Optimally, this understanding would provide improved fracture design criteria to maximize fracture surface area in contact with the reservoir with respect to volume injected and other treatment parameters. In formations with a positive containment condition (i.e., where fracturing out of zone is not anticipated), long penetrating fractures could be used effectively to develop the resource. For the opposite case, the options would beto use a small treatment so that large volumes are not wasted in out-of-zone fracturing and to accept a lower productivity improvement, orto reject the zone as uneconomical. These decisions cannot be made satisfactorily unless criteria for vertical fracture propagation are developed and techniques for readily measuring the important parameters are available. Currently, both theoretical and experimental efforts are being pursued to determine the important parameters and their relative effects on fracture growth. Two modes of fracture containment are possible. One is the situation where fracture growth is terminated at a discrete interface. Examples of this include laboratory experiments showing fracture termination at weak or unbonded interfaces and theoretical models that predict that fracture growth will terminate at a material property interface. The other mode may occur when the fracture propagates into the bounding layer, but extensive growth does not take place and the fracture thus is restricted. An example is the propagation of the fracture into a region with an adverse stress gradient so that continued propagation results in higher stresses on the fracture, which limits growth, as suggested by Simonson et al. and as seen in mineback experiments. Another example is the possible restriction caused by propagation into a higher modulus region where the decreased width results in increased pressure drop in the fracture, which might inhibit extensive growth into that region relative to the lower modulus region. Other parameters, such as natural fractures, treatment parameters, pore pressure, etc., may affect either of these modes. Laboratory and mineback experiments have shown that weak interfaces and in-situ stress differences are the most likely factors to contain the fracture, and weak interfaces are probably effective only at shallow depths. Thus, our experiments are being performed to determine the effect of in-situ stresses on fracture containment, both in a uniform rock sample and at material properly interfaces. SPEJ P. 333^


Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. WA157-WA168 ◽  
Author(s):  
Osni Bastos de Paula ◽  
Marina Pervukhina ◽  
Dina Makarynska ◽  
Boris Gurevich

Modeling dispersion and attenuation of elastic waves in fluid-saturated rocks due to squirt flow requires the knowledge of a number of geometrical parameters of the pore space, in particular, the characteristic aspect ratio of the pores. These parameters are usually inferred by fitting measurements on saturated rocks to model predictions. To eliminate such fitting and thus make the model more predictive, we propose to recover the geometrical parameters of the pore space from the pressure dependency of elastic moduli on dry samples. Our analysis showed that the pressure dependency of elastic properties of rocks (and their deviation from Gassmann’s prediction) at ultrasonic frequencies is controlled by the squirt flow between equant, stiff, and so-called intermediate pores (with aspect ratios between [Formula: see text]). Such intermediate porosity is expected to close at confining pressures of between 200 and 2000 MPa, and thus cannot be directly obtained from ultrasonic experiments performed at pressures below 50 MPa. However, the presence of this intermediate porosity is inferred from the significant linear trend in the pressure dependency of elastic properties of the dry rock and the difference between the bulk modulus of the dry rock computed for spherical pores and the measured modulus at 50 MPa. Moreover, we can infer the magnitude of the intermediate porosity and its characteristic aspect ratio. Substituting these parameters into the squirt model, we have computed elastic moduli and velocities of the water-saturated rock and compared these predictions against laboratory measurements of these velocities. The agreement is good for a number of clean sandstones, but not unexpectedly worse for a broad range of shaley sandstones. Our predictions showed that dispersion and attenuation caused by the squirt flow between compliant and stiff pores may occur in the seismic frequency band. Confirmation of this prediction requires laboratory measurements of elastic properties at these frequencies.


Energies ◽  
2021 ◽  
Vol 14 (17) ◽  
pp. 5514
Author(s):  
Dariusz Knez ◽  
Mohammad Ahmad Mahmoudi Zamani

In this research, two empirical correlations have been introduced to calculate the dynamic Biot coefficients of low-porosity and high-porosity sandstone samples from two open pit mines located in South-West Poland. The experiments were conducted using an acoustic velocity measurement apparatus. Under the undrained condition, firstly, the confining pressure was increased in increments of 200 psi, and the values of pore pressure and dynamic elastic modulus were recorded. This experiment was continued until the Skempton coefficient remained in the range of 0.98–1. Secondly, an experiment on the same sample was conducted under drained conditions, and the corresponding dynamic elastic moduli were calculated. Then, using the calculated dynamic elastic moduli, the dynamic Biot coefficient was determined for each sample under different confining pressure. Finally, two empirical correlations were formulated for each sandstone category. The results demonstrate that, as the confining pressure increases, the Biot coefficient decreases from 0.79 to 0.50 and from 0.84 to 0.45 for low-porosity and high-porosity samples, respectively. Furthermore, as the porosity increases, the sandstone behavior increasingly approaches that of soil. The empirical correlations can be used for sandstone formations with the same porosity in projects where there is not a measurement method for the Biot coefficient.


2020 ◽  
Vol 39 (12) ◽  
pp. 893-900
Author(s):  
Inga Berre ◽  
Ivar Stefansson ◽  
Eirik Keilegavlen

Hydraulic stimulation of geothermal reservoirs in low-permeability basement and crystalline igneous rock can enhance permeability by reactivation and shear dilation of existing fractures. The process is characterized by interaction between fluid flow, deformation, and the fractured structure of the formation. The flow is highly affected by the fracture network, which in turn is deformed because of hydromechanical stress changes caused by the fluid injection. This process-structure interaction is decisive for the outcome of hydraulic stimulation, and, in analysis of governing mechanisms, physics-based modeling has potential to complement field and experimental data. Here, we show how recently developed simulation technology is a valuable tool to understand governing mechanisms of hydromechanical coupled processes and the reactivation and deformation of faults. The methodology fully couples flow in faults and matrix with poroelastic matrix deformation and a contact mechanics model for the faults, including dilation because of slip. Key elements are high aspect ratios of faults and strong nonlinearities in highly coupled governing equations. Example simulations using our open-source software illustrate direct and indirect hydraulic fault reactivation and corresponding permeability enhancement. We investigate the effect of the fault and matrix permeability and the Biot coefficient. A higher matrix permeability leads to more leakage from a permeable fault and thus suppresses reactivation and slip of the fault compared to the case with a lower matrix permeability. If a fault is a barrier to flow, increase of pressure because of the fluid injection results in stabilization of the fault; the situation is opposite if the fault is highly permeable compared to the matrix. For the given setup, lowering the Biot coefficient results in more slip than the base case. While conceptually simple, the examples illustrate the strong hydromechanical couplings and the prospects of physics-based numerical models in investigating the dynamics.


Sign in / Sign up

Export Citation Format

Share Document