Detection of Natural Fractures from Observed Surface Seismic Data Based on a Linear-Slip Model

2018 ◽  
Vol 175 (8) ◽  
pp. 2769-2784 ◽  
Author(s):  
Huaizhen Chen ◽  
Guangzhi Zhang
Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. R85-R98 ◽  
Author(s):  
Xinpeng Pan ◽  
Guangzhi Zhang

Detection of fracture and fluid properties from subsurface azimuthal seismic data improves our abilities to characterize the saturated porous reservoirs with aligned fractures. Motivated by the fracture detection and fluid identification in a fractured porous medium, we have developed a feasible approach to perform a rock physics model-based amplitude variation with offset and azimuth (AVOAz) inversion for the fracture and fluid parameters in a horizontal transversely isotropic (HTI) medium using the PP-wave angle gathers along different azimuths. Based on the linear-slip model, we first use anisotropic Gassmann’s equation to derive the expressions of saturated stiffness components and their perturbations of first-order approximation in terms of elastic properties of an isotropic porous background and fracture compliances induced by a single set of rotationally invariant fractures. We then derive a linearized PP-wave reflection coefficient in terms of fluid modulus, dry-rock matrix term, porosity, density, and fracture compliances or quasi-compliances for an interface separating two weakly HTI media based on the Born scattering theory. Finally, we solve the AVOAz inverse problems iteratively constrained by the Cauchy-sparse regularization and the low-frequency regularization in a Bayesian framework. The results demonstrate that the fluid modulus and fracture quasi-compliances are reasonably estimated in the case of synthetic and real seismic data containing moderate noise in a gas-filled fractured porous reservoir.


Geophysics ◽  
2021 ◽  
pp. 1-97
Author(s):  
kai lin ◽  
Bo Zhang ◽  
Jianjun Zhang ◽  
Huijing Fang ◽  
Kefeng Xi ◽  
...  

The azimuth of fractures and in-situ horizontal stress are important factors in planning horizontal wells and hydraulic fracturing for unconventional resources plays. The azimuth of natural fractures can be directly obtained by analyzing image logs. The azimuth of the maximum horizontal stress σH can be predicted by analyzing the induced fractures on image logs. The clustering of micro-seismic events can also be used to predict the azimuth of in-situ maximum horizontal stress. However, the azimuth of natural fractures and the in-situ maximum horizontal stress obtained from both image logs and micro-seismic events are limited to the wellbore locations. Wide azimuth seismic data provides an alternative way to predict the azimuth of natural fractures and maximum in-situ horizontal stress if the seismic attributes are properly calibrated with interpretations from well logs and microseismic data. To predict the azimuth of natural fractures and in-situ maximum horizontal stress, we focus our analysis on correlating the seismic attributes computed from pre-stack and post-stack seismic data with the interpreted azimuth obtained from image logs and microseismic data. The application indicates that the strike of the most positive principal curvature k1 can be used as an indicator for the azimuth of natural fractures within our study area. The azimuthal anisotropy of the dominant frequency component if offset vector title (OVT) seismic data can be used to predict the azimuth of maximum in-situ horizontal stress within our study area that is located the southern region of the Sichuan Basin, China. The predicted azimuths provide important information for the following well planning and hydraulic fracturing.


2019 ◽  
Vol 38 (4) ◽  
pp. 268-273
Author(s):  
Maheswara Phani ◽  
Sushobhan Dutta ◽  
Kondal Reddy ◽  
Sreedurga Somasundaram

Raageshwari Deep Gas (RDG) Field is situated in the southern part of the Barmer Basin in Rajasthan, India, at a depth of 3000 m. With both clastic and volcanic lithologies, the main reservoirs are tight, and hydraulic fracturing is required to enhance productivity, especially to improve permeability through interaction of induced fractures with natural fractures. Therefore, optimal development of the RDG Field reservoirs requires characterization of faults and natural fractures. To address this challenge, a wide-azimuth 3D seismic data set over the RDG Field was processed to sharply define faults and capture anisotropy related to open natural fractures. Anisotropy was indicated by the characteristic sinusoidal nature of gather reflection events processed using conventional tilted transverse imaging (TTI); accordingly, we used orthorhombic imaging to correct for these, to quantify fracture-related anisotropy, and to yield a more correct subsurface image. During prestack depth migration (PSDM) processing of the RDG data, TTI and orthorhombic velocity modeling was done with azimuthal sectoring of reflection arrivals. The resultant PSDM data using this velocity model show substantial improvement in image quality and vertical resolution at the reservoir level compared to vintage seismic data. The improved data quality enabled analysis of specialized seismic attributes like curvature and thinned fault likelihood for more reliable characterization of faults and fractures. These attributes delineate the location and distribution of probable fracture networks within the volcanic reservoirs. Interpreted subtle faults, associated with fracture zones, were validated with microseismic, production, and image log data.


Geophysics ◽  
1998 ◽  
Vol 63 (1) ◽  
pp. 139-148 ◽  
Author(s):  
Jeroen Groenenboom ◽  
Jacob T. Fokkema

During scaled hydraulic fracturing experiments in our laboratory, the fracture growth process is monitored in a time‐lapse experiment with ultrasonic waves. We observe dispersion of compressional waves that have propagated across the hydraulic fracture. This dispersion appears to be related to the width of the hydraulic fracture. This means that we can apply the dispersion measurements to monitor the width of the hydraulic fracture in an indirect manner. For a direct determination of the width, the resolution of the signal is required to distinguish the reflections that are related with two distinct fluid/solid interfaces delimiting the hydraulic fracture from its solid embedding. To make this distinction, the solid/fluid interfaces must be separated at least one eighth of a wavelength and represent sufficient impedance contrast. The applicability of the indirect dispersion measurement method however, extends to a fracture width that is in the order of 1% of the incident wavelength. The time‐lapse ultrasonic measurements allow us to relate the small difference in arrival time and amplitude between two measurements solely to the small changes in the width of the fracture. Additional experimental data show that shear waves are completely shadowed by hydraulic fractures, indicating that there is no acoustic contact mechanism at the fracture interface. Therefore we think it is appropriate to use a thin fluid‐filled layer model for these hydraulic fractures instead of the standard empirically oriented linear slip model. Nevertheless, the thin layer model is consistent with the linear‐slip model, if interpreted correctly. A comparison of width measurements inside the wellbore and width estimates by means of dispersion measurements close to the wellbore shows that the method can be successfully applied, at least under laboratory conditions, and that small changes in the width of the fracture are directly expressed in the dispersion of the transmitted signal. This opens the way for the important new application of width monitoring of hydraulic fractures.


Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. T57-T62 ◽  
Author(s):  
Chunling Wu ◽  
Jerry M. Harris ◽  
Kurt T. Nihei ◽  
Seiji Nakagawa

Within the context of seismic wave propagation, fractures can be described as thin layers or linear-slip interfaces. In this paper, numerical simulations of elastic wave propagation in a medium with a single fracture represented by these two models are performed by 2D finite-difference codes: a variable-grid isotropic code for the thin-layer model and a regular-grid anisotropic code for the linear-slip model. Numerical results show excellent agreement between the two models for wavefields away from the fracture; the only discrepancy between the two is the presence of a slow wave traveling primarily within the fracture fluid of the thin-layer model. The comparison of the computational cost shows that modeling of the linear-slip model is more efficient than that of the thin-layer model. This study demonstrates that the linear-slip model is an efficient and accurate modeling approach for the remote seismic characterization of fractures.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. WB165-WB176 ◽  
Author(s):  
Tatiana Chichinina ◽  
Irina Obolentseva ◽  
Leonid Gik ◽  
Boris Bobrov ◽  
Gerardo Ronquillo-Jarillo

This study attempts to validate a mathematical formalism of introducing attenuation into Schoenberg’s linear slip model. This formalism is based on replacing the real-valued weaknesses by complex-valued ones. During an ultrasonic experiment, performed at a central frequency of [Formula: see text] on a plate-stack model with [Formula: see text]-thick Plexiglas™ plates, the velocity and attenuation (inverse of the quality factor [Formula: see text]) of P-, SH-, and SV-waves are measured in directions from 25° to 90° with the symmetry axis for dry and oil-saturated models and loading uniaxial pressures of 2 and [Formula: see text]. The velocity and attenuation data are fitted by the derived theoretical functions. The values of the real and imaginary parts of the complex-valued weaknesses are estimated. Thereal parts of the weaknesses, which have a clear physical meaning (they affect the weakening of the material), are three times larger for the dry model than for the oil-saturated one. The imaginary parts of the weaknesses are responsible for attenuation; their values are an order of magnitude smaller than the real parts. The derived expressions for angle-dependent velocities and attenuations can be used to distinguish between dry and oil-saturated fractures. In particular, the P-wave attenuation function in the symmetry-axis direction (normal to fracture planes) is different in dry and saturated media. The experiment shows that the plate-stack model is inhomogeneous because of the nonuniform pressure distribution, which degrades the experimental results and creates difficulties in the inversion for the complex-valued weaknesses — particularly in joint inversion of P- and S-wave data.


2016 ◽  
Vol 43 (15) ◽  
pp. 7973-7981 ◽  
Author(s):  
Konstantinos Chousianitis ◽  
A. Ozgun Konca ◽  
G.-Akis Tselentis ◽  
Gerassimos A. Papadopoulos ◽  
Michail Gianniou

2012 ◽  
Vol 52 (1) ◽  
pp. 455 ◽  
Author(s):  
Adam Bailey ◽  
Rosalind King ◽  
Guillaume Backé

Understanding natural fracture networks has increasingly been recognised as an important factor for the prospectivity of a geothermal play, as they commonly exert a prime control over permeability at depth. The onshore Northern Perth Basin provides a good example of how fracture stimulation, and subsequent enhancement of the structural permeability, during hydrocarbon production can enhance flow rate from original tight gas reservoirs. Low primary porosity and permeability values have been initially recorded in the Northern Perth Basin due to silica-rich groundwater infiltration and consequent quartz cementation. Geothermal energy prospectivity in the region will therefore depend heavily on similar engineering techniques or on the presence of secondary permeability due to interconnected natural fractures. The existence and extent of these natural fractures are verified in this study through an integrated analysis of geophysical logs (including wellbore image logs), wells tests, and 3D seismic data. Wellbore image logs from 11 petroleum wells in the Northern Perth Basin are used to identify borehole failure (such as borehole breakout and drilling-induced tensile fractures) to give a present-day maximum horizontal stress orientation of N076°E (with an s.d. of 13°). Density logs and leak off tests from 13 petroleum wells are used to constrain the present-day stress magnitudes, giving a transitional strike-slip fault to reverse-fault stress regime in the Northern Perth Basin. 870 fractures are identified in image logs from 13 petroleum wells in the Northern Perth Basin, striking roughly north to south and northwest to northeast. Fractures aligned in the present-day stress field are optimally oriented for reactivation, and are hence likely to be open to fluid flow. Electrically resistive and conductive natural fractures are identified on the wellbore image logs. Resistive fractures are considered to be cemented with electrically resistive cement (such as quartz or calcite) and thus closed to fluid-flow. Conductive fractures are considered to be uncemented and open to fluid-flow, and are thus important to geothermal exploration. Fracture susceptibility diagrams constructed for the identified fractures illustrate that the conductive fractures are optimally oriented for reactivation in the present-day strike-slip fault to reverse-fault stress regime, and so are likely to be open to fluid flow. This is reinforced by the correlation of drilling fluid loss and conductive natural fractures in three wells in the Northern Perth Basin. To gain an understanding of the extent and interconnectedness of these fractures, it is necessary to look at more regional data, such as 3D seismic surveys. It is, however, well-documented that fault and fracture networks like those generally observed in image logs lie well below seismic amplitude resolution, making them difficult to observe directly on amplitude data. Seismic attributes can be calculated to provide some information on sub-seismic scale structural and stratigraphic features. Using a 3D seismic cube acquired over the Dongara North gas field, attribute maps of complex multi-trace dip-steered coherency and most positive curvature were used to document the presence of natural fractures and to best constrain the likely extent of the fracture network. The resulting fracture network model displays relatively good connectivity, which is likely to extend across much of the basin. These optimally oriented fractures are therefore likely to form a secondary permeability network throughout the cemented sediments of the Northern Perth Basin, offering potential deep fluid flow conduits, which may be exploited for the production of geothermal energy.


Sign in / Sign up

Export Citation Format

Share Document