scholarly journals On the link between stress field and small-scale hydraulic fracture growth in anisotropic rock derived from microseismicity

Author(s):  
Valentin S. Gischig ◽  
Joseph Doetsch ◽  
Hansruedi Maurer ◽  
Hannes Krietsch ◽  
Florian Amann ◽  
...  

Abstract. To characterize the stress field at the Grimsel Test Site (GTS) underground rock laboratory a series of hydrofracturing test and overcoring test were performed. Hydrofracturing was accompanied by seismic monitoring using a network of highly sensitive piezo sensors and accelerometers that were able to record small seismic events associated with decimeter-sized fractures. Due to potential discrepancies between the hydro-fracture orientation and stress field estimates from overcoring, it was essential to obtain high-precision hypocenter locations that reliably illuminate fracture growth. Absolute locations were improved using a transverse isotropic P-wave velocity model and by applying joint hypocenter determination that allowed computation of station corrections. We further exploited the high degree of waveform similarity of events by applying cluster analysis and relative relocation. Resulting clouds of absolute and relative located seismicity showed a consistent east-west strike and 70° dip for all hydro-fractures. The fracture growth direction from microseismicity is consistent with the principal stress orientations from the overcoring stress tests provided an anisotropic elastic model for the rock mass is used in the data inversions. σ1 is significantly larger than the other two principal stresses, and has a reasonably well-defined orientation that is subparallel to the fracture plane. σ2 and σ3 are almost equal in magnitude, and thus lie on a circle defined by the standard errors of the solutions. The poles of the microseismicity planes also lie on this circle towards the north. The trace of the hydraulic fracture imaged at the borehole wall show that they initiated within the foliation plane, which differs in orientation from the microseismicity planes. Thus, fracture initiation was most likely influenced by a foliation-related strength anisotropy. Analysis of P-wave polarizations suggested double-couple focal mechanisms with both thrust and normal faulting mechanisms present, whereas strike-slip and thrust mechanisms would be expected from the overcoring-derived stress solution. The reasons for these discrepancies are not well understood, but may involve stress field rotation around the propagating hydrofracture. Our study demonstrates that microseismicity monitoring along with high-resolution event locations provides valuable information for interpreting stress characterization measurements.

Solid Earth ◽  
2018 ◽  
Vol 9 (1) ◽  
pp. 39-61 ◽  
Author(s):  
Valentin Samuel Gischig ◽  
Joseph Doetsch ◽  
Hansruedi Maurer ◽  
Hannes Krietsch ◽  
Florian Amann ◽  
...  

Abstract. To characterize the stress field at the Grimsel Test Site (GTS) underground rock laboratory, a series of hydrofracturing and overcoring tests were performed. Hydrofracturing was accompanied by seismic monitoring using a network of highly sensitive piezosensors and accelerometers that were able to record small seismic events associated with metre-sized fractures. Due to potential discrepancies between the hydrofracture orientation and stress field estimates from overcoring, it was essential to obtain high-precision hypocentre locations that reliably illuminate fracture growth. Absolute locations were improved using a transverse isotropic P-wave velocity model and by applying joint hypocentre determination that allowed for the computation of station corrections. We further exploited the high degree of waveform similarity of events by applying cluster analysis and relative relocation. Resulting clouds of absolute and relative located seismicity showed a consistent east–west strike and 70° dip for all hydrofractures. The fracture growth direction from microseismicity is consistent with the principal stress orientations from the overcoring stress tests, provided that an anisotropic elastic model for the rock mass is used in the data inversions. The σ1 stress is significantly larger than the other two principal stresses and has a reasonably well-defined orientation that is subparallel to the fracture plane; σ2 and σ3 are almost equal in magnitude and thus lie on a circle defined by the standard errors of the solutions. The poles of the microseismicity planes also lie on this circle towards the north. Analysis of P-wave polarizations suggested double-couple focal mechanisms with both thrust and normal faulting mechanisms present, whereas strike-slip and thrust mechanisms would be expected from the overcoring-derived stress solution. The reasons for these discrepancies can be explained by pressure leak-off, but possibly may also involve stress field rotation around the propagating hydrofracture. Our study demonstrates that microseismicity monitoring along with high-resolution event locations provides valuable information for interpreting stress characterization measurements.


1981 ◽  
Vol 21 (01) ◽  
pp. 21-29 ◽  
Author(s):  
Gordon D. Anderson

Abstract Small-scale laboratory experiments were performed to study the growth of hydraulically driven fractures in the vicinity of an unbonded interface in rocks. The purpose was to evaluate under which conditions the hydraulic fractures would cross the interface. The materials used in these studies were Nugget sandstone from Utah (3 to 6% porosity) and Indiana limestone (12 to 15% porosity). The fracturing fluid was oil (viscosity appx. 300 cp) injected into the rock through high-pressure steel tubing. Prismatic blocks of the rock materials to be studied were held adjacent to one another in a hydraulic press so that a normal stress was set up across their mutual interface. Lubricants and surface roughening were used to vary the frictional properties of the interfaces. It was found that as the interface surface friction coefficient was decreased, the normal stress had to be increased for a hydraulic fracture to cross the interface. The frictional shear stress that the interface can support without slippage appears to be critical in determining fracture growth across the interface. Additional experiments were performed to evaluate the coefficient of friction for the different interface surface preparations used. These experiments demonstrated that a variation in the frictional properties along an interfacial surface in the vicinity of hydraulic fracture growth can alter the path of the fracture. The experiments also demonstrated that cracks, which intersect the interface from the side opposite the approaching hydraulic fracture, can impede fracture growth across the interface. Introduction Hydraulic fracturing and a variant - massive hydraulic fracturing (MHF) - are primary candidates for stimulating production from the tight-gas reservoirs in the U.S. Hydraulic fracturing has been used widely as a well completion technique for about 30 years. MHF is a more recent application that differs from hydraulic fracturing in that larger quantities of fluid and proppant are pumped to create more extensive fractures in the reservoirs. Application of MHF to increase production from the tight reservoirs has provided mixed and, in many cases, disappointing results, especially in lenticular reservoirs. For MHF to be successful in enhancing the production of gas from tight reservoirs, it is important that the fractures be emplaced in productive reservoir rock providing large drainage surfaces in the low-permeability material and conductive channels back to the wellbore. We then are faced with the problem of containing fractures in a given formation.Under the U.S. DOE'S Unconventional Gas Recovery program, Lawrence Livermore Natl. Laboratory is conducting a research program to study the hydraulic fracture process. The general goal of this research is to determine if and to what extent the reservoir parameters control the geometry of the created fractures. These reservoir parameters include (1) the mechanical properties of the rock (i.e., elastic moduli, mechanical strength, etc.), (2) the physical state of the rock (i.e., presence of pre-existing cracks or faults, porosity, pore fluid, etc.), (3) presence of layering or interfaces between different rock strata, and (4) stress field on the rock. In addition to reservoir parameters, the growth of a hydraulically driven crack will be influenced by (1) the manner in which the driving fluid is injected into the rock, (2) the characteristics of the fracturing fluid (i.e., viscosity, presence of proppant, etc.), and (3) any chemical reaction between the fluid and rock. Previous work has shown that crack orientation is controlled primarily by the in-situ or applied stress field, with crack growth oriented perpendicular to the least principal stress. SPEJ P. 21^


Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 603-611 ◽  
Author(s):  
Jeroen Groenenboom ◽  
Dirkjan B. van Dam

We carry out small‐scale hydraulic fracture experiments to investigate the physics of hydraulic fracturing. The laboratory experiments are combined with time‐lapse ultrasonic measurements with active sources using both compressional and shear‐wave transducers. For the time‐lapse measurements we focus on ultrasonic measurement changes during fracture growth. As a consequence we can detect the hydraulic fracture and characterize its shape and geometry during growth. Hence, this paper deals with fracture characterization using time‐lapse acoustic data. During fracture growth the acoustic waves generate diffractions at the tip of the fracture. The direct compressional and shear diffractions are used to locate the position of the tip of the fracture. More detailed analysis of these diffractions can be used to obtain information on the geometry and configuration of the fracture tip, including the creation of a zone that is not penetrated by fluid. Furthermore, it appears that the acoustic diffraction is generated mainly at the fluid front and only weakly at the dry tip. In addition, the wavefield that has been transmitted through the hydraulic fracture is measured. Shear‐wave transmissions are shadowed because the shear modulus vanishes inside the fluid‐filled fracture. From this observation we conclude that the fracture is mechanically open. In other words, no friction occurs related to the movement of fracture faces that are in mechanical contact. Compressional transmissions show a distinctive dispersion relative to the measurement in the unfractured medium. This dispersion can be used to determine the width (or aperture) of the fracture by fitting the measured dispersion with the theoretical prediction as a function of the unknown fracture width. We show that the width profile of the fracture can be reconstructed by using a set of transmission records with different source and receiver locations. By performing a validation experiment, we show that the width determination method is reliable, although the estimated fracture width is only a few percent of the incident wavelength. The strength of the method relies on time‐lapse measurements combined with fitting the changes in the measured waveforms during the experiment. The combination of diffractions and transmissions helps us visualize the dynamic process of hydraulic fracture growth. Hence, acoustic measurements with active sources prove their usefulness for fracture characterization.


Solid Earth ◽  
2014 ◽  
Vol 5 (1) ◽  
pp. 275-297 ◽  
Author(s):  
W. Friederich ◽  
A. Brüstle ◽  
L. Küperkoch ◽  
T. Meier ◽  
S. Lamara ◽  
...  

Abstract. The lateral variation of the stress field in the southern Aegean plate and the subducting Hellenic slab is determined from recordings of seismicity obtained with the CYCNET and EGELADOS networks in the years from 2002 to 2007. First motions from 7000 well-located microearthquakes were analysed to produce 540 well-constrained focal mechanisms. They were complemented by another 140 derived by waveform matching of records from larger events. Most of these earthquakes fall into 16 distinct spatial clusters distributed over the southern Aegean region. For each cluster, a stress inversion could be carried out yielding consistent estimates of the stress field and its spatial variation. At crustal levels, the stress field is generally dominated by a steeply dipping compressional principal stress direction except in places where coupling of the subducting slab and overlying plate come into play. Tensional principal stresses are generally subhorizontal. Just behind the forearc, the crust is under arc-parallel tension whereas in the volcanic areas around Kos, Columbo and Astypalea tensional and intermediate stresses are nearly degenerate. Further west and north, in the Santorini–Amorgos graben and in the area of the islands of Mykonos, Andros and Tinos, tensional stresses are significant and point around the NW–SE direction. Very similar stress fields are observed in western Turkey with the tensional axis rotated to NNE–SSW. Intermediate-depth earthquakes below 100 km in the Nisyros region indicate that the Hellenic slab experiences slab-parallel tension at these depths. The direction of tension is close to east–west and thus deviates from the local NW-oriented slab dip presumably owing to the segmentation of the slab. Beneath the Cretan sea, at shallower levels, the slab is under NW–SE compression. Tensional principal stresses in the crust exhibit very good alignment with extensional strain rate principal axes derived from GPS velocities except in volcanic areas, where both appear to be unrelated, and in the forearc where compressional principal stresses are very well aligned with compressional principal strain rates. This finding indicates that, except for volcanic areas, microseismic activity in the southern Aegean is not controlled by small-scale local stresses but rather reflects the regional stress field. The lateral and depth variations of the stress field reflect the various agents that influence tectonics in the Aegean: subduction of the Hellenic slab, incipient collision with continental African lithosphere, roll back of the slab in the southeast, segmentation of the slab, arc volcanism and extension of the Aegean crust.


2019 ◽  
Vol 6 ◽  
pp. 50-53
Author(s):  
A.R. Davletova ◽  
◽  
A.I. Fedorov ◽  
G.A. Shchutsky ◽  
◽  
...  

1976 ◽  
Vol 13 (3) ◽  
pp. 212-217 ◽  
Author(s):  
Han-Shou Liu ◽  
Edward S. Chang ◽  
George H. Wyatt

Author(s):  
Suresh Dande ◽  
◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
◽  
...  

Laboratory physical models play an important role in understanding rock properties and wave propagation, both theoretically and at the field scale. In some cases, 3D-printing technology can be adopted to construct complex rock models faster, more inexpensively, and with more specific features than previous model-building techniques. In this study, we use 3D-printed rock models to assist in understanding the effects of various fluids (air, water, engine oil, crude oil, and glycerol) on the models’ elastic properties. We first used a 3D-printed, 1-in. cube-shaped layered model. This model was created with a 6% primary porosity and a bulk density of 0.98 g/cc with VTI anisotropy. We next employed a similar cube but with horizontal inclusions embedded in the layered background, which contributed to its total 24% porosity (including primary porosity). For air to liquid saturation, P-velocities increased for all liquids in both models, with the highest increase being with glycerol (57%) and an approximately 45% increase for other fluids in the inclusion model. For the inclusion model (dry and saturated), we observed a greater difference between two orthogonally polarized S-wave velocities (Vs1 and Vs2) than between two P-wave velocities (VP0 and VP90). We attribute this to the S2-wave (polarized normal to both the layering and the plane of horizontal inclusions), which appears more sensitive to horizontal inclusions than the P-wave. For the inclusion model, Thomsen’s P-wave anisotropic parameter (ɛ) decreased from 26% for the air case to 4% for the water-saturated cube and to 1% for glycerol saturation. The small difference between the bulk modulus of the frame and the pore fluid significantly reduces the velocity anisotropy of the medium, making it almost isotropic. We compared our experimental results with theory and found that predictions using Schoenberg’s linear slip theory combined with Gassmann’s anisotropic equation were closer to actual measurements than Hudson’s isotropic calculations. This work provides insights into the usefulness of 3D-printed models to understand elastic rock properties and wave propagation under various fluid saturations.


Author(s):  
Behrouz Arash ◽  
Quan Wang

Free vibration of single- and double-layered graphene sheets is investigated by employing nonlocal continuum theory and molecular dynamics simulations. Results show that the classical elastic model overestimated the resonant frequencies of the sheets by a percentage as high as 62%. The dependence of small-scale effects, sizes of sheets, boundary conditions, and number of layers on vibrational characteristic of single- and double-layered graphene sheets is studied. The resonant frequencies predicted by the nonlocal elastic plate theory are verified by the molecular dynamics simulations, and the nonlocal parameter is calibrated through the verification process. The simulation results reveal that the calibrated nonlocal parameter depends on boundary conditions and vibrational modes. The nonlocal plate model is found to be indispensable in vibration analysis of grapheme sheets with a length less than 8 nm on their sides.


Sign in / Sign up

Export Citation Format

Share Document