Azimuth Of Maximum Horizontal Stress Research Articles

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi2 (62 km2) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Experimental study of directional propagation of hydraulic fracture guided by multi-radial slim holes

The conventional hydraulic fracturing fails in the target oil/gas zone (remaining oil, closed reservoir, etc.) which is not located in the azimuth of maximum horizontal stress of available wellbores. The technology of directional propagation of hydraulic fracture guided by vertical multi-radial slim holes is innovatively developed. In order to verify this technology, lots of true triaxial hydraulic fracturing simulation experiments were carried out with artificial cores, aiming at investigating the influence of in-situ stress, fracturing fluid displacement, and azimuth, diameter, number, and spacing of radial slim holes on fracture propagation. The results show that directional propagation of hydraulic fractures guided by vertical multi-radial slim holes is feasible. Regardless of varied radial slim hole and in-situ stress parameters, hydraulic fracture always initiates in the heel of radial slim hole. For hydraulic fracturing assisted by single radial slim hole, smaller horizontal stress difference (3 MPa, σH/σh = 1.25) and smaller radial slim hole azimuth (15°) may guide propagation of hydraulic fractures along radial slim hole, and larger horizontal stress difference (6 MPa, σH/σh = 1.67) can sharply reduce guidance of radial slim hole. Affected by mutual interference from guidance of radial slim holes and controlling of horizontal in-situ stress, fractures tend to distort along fracture length and fracture height. For hydraulic fracturing assisted by vertical multi-radial slim holes, horizontal stress difference is one of key factors influencing the directional propagation effect of hydraulic fracture. Larger horizontal in-situ stress difference (>6 MPa, σH/σh > 1.67) hardly results in directional propagation of hydraulic fracture along radial slim hole row, and creates bigger diversion angle of fracture, and larger azimuth (>45°) of radial slim hole row reduces the guidance strength, resulting small fracture diversion radius. In conditions of larger angle between target zone and maximum horizontal stress, and requirement of large azimuth of radial slim hole row, the effective hydraulic fracture propagation along radial slim holes and ideal fracture height will be achieved through human intervention, e.g. increment of the number and diameter of radial slim holes, reduction of radial slim holes spacing, increment of fracturing fluid displacement, etc. The simulation results provide methods for designing directional propagation of hydraulic fractures guided by vertical multi-radial slim holes.

Seismic azimuthal anisotropy analysis after hydraulic fracturing

Many tight sandstone, limestone, and shale reservoirs require hydraulic fracturing to provide pathways that allow hydrocarbons to reach the well bore. Most of these tight reservoirs are now produced using multiple stages of fracturing through horizontal wells drilled perpendicular to the present-day azimuth of maximum horizontal stress. In a homogeneous media, the induced fractures are thought to propagate perpendicularly to the well, parallel to the azimuth of maximum horizontal stress, thereby efficiently fracturing the rock and draining the reservoir. We evaluated what may be the first anisotropic analysis of a Barnett shale-gas reservoir after extensive hydraulic fracturing and focus on mapping the orientation and intensity of induced fractures and any preexisting factures, with the objective being the identification of reservoir compartmentalization and bypassed pay. The Barnett Shale we studied has near-zero permeability and few if any open natural fractures. We therefore hypothesized that anisotropy is therefore due to the regional northeast–southwest maximum horizontal stress and subsequent hydraulic fracturing. We found the anisotropy to be highly compartmentalized, with the compartment edges being defined by ridges and domes delineated by the most positive principal curvature [Formula: see text]. Microseismic work by others in the same survey indicates that these ridges contain healed natural fractures that form fracture barriers. Mapping such heterogeneous anisotropy field could be critical in planning the location and direction of any future horizontal wells to restimulate the reservoir as production drops.

Determination of magnitude and orientation of the in-situ stress from borehole breakout and effect of pore pressure on borehole stability — Case study in Cheshmeh Khush oil field of Iran

The detailed profile of the orientations of borehole breakouts as a function of depth were investigated in wells CK-8 and CK-9 in the Cheshmeh Khush oilfield in south-west of Iran. Borehole breakouts without tensile induced fractures were detected from UBI logs. By analyzing statistically of borehole breakout orientation in two wells, it was found that while a mean azimuth of maximum horizontal stress in CK-8 is N–S (s.d = 0.4°); the azimuth of breakout in CK-9 varies as a function of depth with a mean azimuth 345° + 26°; which since no 3D seismic survey has been run upon the field, assumed to be caused by active tectonic state. By utilizing presence and absence of borehole breakout and tensile induced fracture, as well as density, rock strength, pore and mud pressure in deviated borehole, in situ stresses magnitude in CK-8 were constrained. The experiment showed that the state of stress is normal faulting regime and maximum horizontal stress magnitude is closer to SV comparing with S hmin. Pore pressure was predicted by deploying Eaton's method on observed sonic data deviation from normal compaction trend. Considering the behavior of breakout and pore pressure in CK-8 in accordance with rock strength changes, demonstrated that the borehole gets more stable as pore pressure increases.

Some influences of stratigraphy and structure on reservoir stress orientation

The azimuth of maximum horizontal stress in a reservoir can vary significantly with depth and with position on a subsurface structure. We present and discuss evidence from field data for such variation and demonstrate both analytically and with finite‐element modeling how such changes might take place. Under boundary conditions of uniform far‐field displacement, changes in stratigraphic layering can reorient the principal stress direction if the formation is intrinsically anisotropic. If the formation stiffness is lower perpendicular to bedding than parallel to bedding (as is often the case in layered geologic media), an increase in dip will reduce the component of compressive stress in the dip azimuth direction. Folds can reorient principal stresses because flexural strain varies with depth and position. Compressive stress perpendicular to a fold axis increases with depth at the crest of an anticline and decreases with depth at the limb. When the regional stress anisotropy is weak, this change in stress magnitude can reorient the local principal stress directions. Numerical simulations of such effects gave results consistent with changes in stress orientation at the Cymric and Lost Hills oil fields in California as observed via shear‐wave polarization analyses and tiltmeter surveys of hydraulic fracturing. Knowledge of such variation of stress direction with depth and structural position is critical for drilling, completions, hydraulic fracture, and well pattern designs.