Interpretation of fractures and stress anisotropy in Marcellus Shale using multicomponent seismic data

2014 ◽  
Vol 2 (2) ◽  
pp. SE105-SE115 ◽  
Author(s):  
Mehdi E. Far ◽  
Bob Hardage
Keyword(s):  
Seismic Data ◽  
Marcellus Shale ◽  
P Wave ◽  
Accurate Information ◽  
Amplitude Analysis ◽  
Data Set ◽  
Regional Stress ◽  
Stress Anisotropy ◽  
Multicomponent Seismic ◽  
Principal Directions

Using a data set from the Marcellus Shale, we evaluated the advantages of multicomponent seismic data for fracture and anisotropy studies over conventional P-wave data. Using traveltime and amplitude analysis on pre- and poststack seismic data, we concluded that PS-waves can provide more accurate information about the location, orientation, and intensity of natural fractures and stress anisotropy than P-waves. Our analysis indicated that regional stress was the main cause of velocity anisotropy. Amplitude variation with offset and azimuth appeared to be more useful for fracture studies, whereas traveltime variations (especially PS-waves) provided a better indication of regional stress orientations. Principal directions for amplitudes and traveltimes of PP- and PS-waves were different. Misalignment of PP- and PS-waves principal directions suggested that the simplest, most realistic anisotropy model for the fractured Marcellus is monoclinic symmetry.

scholarly journals Semiglobal viscoacoustic full-waveform inversion

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. R271-R293 ◽  
Author(s):  
Nuno V. da Silva ◽  
Gang Yao ◽  
Michael Warner
Keyword(s):  
Physical Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Waveform Inversion ◽  
P Wave ◽  
Full Waveform Inversion ◽  
Data Set ◽  
Intrinsic Attenuation ◽  
Full Waveform ◽  
P Wave Velocity

Full-waveform inversion deals with estimating physical properties of the earth’s subsurface by matching simulated to recorded seismic data. Intrinsic attenuation in the medium leads to the dispersion of propagating waves and the absorption of energy — media with this type of rheology are not perfectly elastic. Accounting for that effect is necessary to simulate wave propagation in realistic geologic media, leading to the need to estimate intrinsic attenuation from the seismic data. That increases the complexity of the constitutive laws leading to additional issues related to the ill-posed nature of the inverse problem. In particular, the joint estimation of several physical properties increases the null space of the parameter space, leading to a larger domain of ambiguity and increasing the number of different models that can equally well explain the data. We have evaluated a method for the joint inversion of velocity and intrinsic attenuation using semiglobal inversion; this combines quantum particle-swarm optimization for the estimation of the intrinsic attenuation with nested gradient-descent iterations for the estimation of the P-wave velocity. This approach takes advantage of the fact that some physical properties, and in particular the intrinsic attenuation, can be represented using a reduced basis, substantially decreasing the dimension of the search space. We determine the feasibility of the method and its robustness to ambiguity with 2D synthetic examples. The 3D inversion of a field data set for a geologic medium with transversely isotropic anisotropy in velocity indicates the feasibility of the method for inverting large-scale real seismic data and improving the data fitting. The principal benefits of the semiglobal multiparameter inversion are the recovery of the intrinsic attenuation from the data and the recovery of the true undispersed infinite-frequency P-wave velocity, while mitigating ambiguity between the estimated parameters.

A robust workflow for performing joint impedance inversion and its applications

2020 ◽  
Vol 8 (1) ◽  
pp. T141-T149
Author(s):  
Ritesh Kumar Sharma ◽  
Satinder Chopra ◽  
Larry R. Lines
Keyword(s):  
Time Domain ◽  
Seismic Data ◽  
Sedimentary Basin ◽  
Vertical Component ◽  
Data Sets ◽  
Data Set ◽  
Matching Process ◽  
Multicomponent Seismic ◽  
Impedance Inversion ◽  
Straightforward Approach

Multicomponent seismic data offer several advantages for characterizing reservoirs with the use of the vertical component (PP) and mode-converted (PS) data. Joint impedance inversion inverts both of these data sets simultaneously; hence, it is considered superior to simultaneous impedance inversion. However, the success of joint impedance inversion depends on how accurately the PS data are mapped on the PP time domain. Normally, this is attempted by performing well-to-seismic ties for PP and PS data sets and matching different horizons picked on PP and PS data. Although it seems to be a straightforward approach, there are a few issues associated with it. One of them is the lower resolution of the PS data compared with the PP data that presents difficulties in the correlation of the equivalent reflection events on both the data sets. Even after a few consistent horizons get tracked, the horizon matching process introduces some artifacts on the PS data when mapped into PP time. We have evaluated such challenges using a data set from the Western Canadian Sedimentary Basin and then develop a novel workflow for addressing them. The importance of our workflow was determined by comparing data examples generated with and without its adoption.

Determination of the principal directions of azimuthal anisotropy from P-wave seismic data

Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 692-706 ◽  
Author(s):  
Subhashis Mallick ◽  
Kenneth L. Craft ◽  
Laurent J. Meister ◽  
Ronald E. Chambers
Keyword(s):  
Anisotropic Medium ◽  
Seismic Data ◽  
Principal Direction ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Amplitude Variation ◽  
S Waves ◽  
Wave Data ◽  
Principal Directions ◽  
Qualitative Measure

In an azimuthally anisotropic medium, the principal directions of azimuthal anisotropy are the directions along which the quasi-P- and the quasi-S-waves propagate as pure P and S modes. When azimuthal anisotropy is induced by oriented vertical fractures imposed on an azimuthally isotropic background, two of these principal directions correspond to the directions parallel and perpendicular to the fractures. S-waves propagating through an azimuthally anisotropic medium are sensitive to the direction of their propagation with respect to the principal directions. As a result, primary or mode‐converted multicomponent S-wave data are used to obtain the principal directions. Apart from high acquisition cost, processing and interpretation of multicomponent data require a technology that the seismic industry has not fully developed. Anisotropy detection from conventional P-wave data, on the other hand, has been limited to a few qualitative studies of the amplitude variation with offset (AVO) for different azimuthal directions. To quantify the azimuthal AVO, we studied the amplitude variation with azimuth for P-wave data at fixed offsets. Our results show that such amplitude variation with azimuth is periodic in 2θ, θ being the orientation of the shooting direction with respect to one of the principal directions. For fracture‐induced anisotropy, this principal direction corresponds to the direction parallel or perpendicular to the fractures. We use this periodic azimuthal dependence of P-wave reflection amplitudes to identify two distinct cases of anisotropy detection. The first case is an exactly determined one, where we have observations from three azimuthal lines for every common‐midpoint (CMP) location. We derive equations to compute the orientation of the principal directions for such a case. The second case is an overdetermined one where we have observations from more than three azimuthal lines. Orientation of the principal direction from such an overdetermined case can be obtained from a least‐squares fit to the reflection amplitudes over all the azimuthal directions or by solving many exactly determined problems. In addition to the orientation angle, a qualitative measure of the degree of azimuthal anisotropy can also be obtained from either of the above two cases. When azimuthal anisotropy is induced by oriented vertical fractures, this qualitative measure of anisotropy is proportional to fracture density. Using synthetic seismograms, we demonstrate the robustness of our method in evaluating the principal directions from conventional P-wave seismic data. We also apply our technique to real P-wave data, collected over a wide source‐to‐receiver azimuth distribution. Computations using our method gave an orientation of the principal direction consistent with the general fracture orientation in the area as inferred from other geological and geophysical evidence.

Elastic impedance parameterization and inversion with Young’s modulus and Poisson’s ratio

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. N35-N42 ◽  
Author(s):  
Zhaoyun Zong ◽  
Xingyao Yin ◽  
Guochen Wu
Keyword(s):  
Parameter Estimation ◽  
Young’S Modulus ◽  
Seismic Data ◽  
Poisson’S Ratio ◽  
Young's Modulus ◽  
Incident Angle ◽  
P Wave ◽  
Poisson's Ratio ◽  
Data Set ◽  
Elastic Impedance

Young’s modulus and Poisson’s ratio are related to quantitative reservoir properties such as porosity, rock strength, mineral and total organic carbon content, and they can be used to infer preferential drilling locations or sweet spots. Conventionally, they are computed and estimated with a rock physics law in terms of P-wave, S-wave impedances/velocities, and density which may be directly inverted with prestack seismic data. However, the density term imbedded in Young’s modulus is difficult to estimate because it is less sensitive to seismic-amplitude variations, and the indirect way can create more uncertainty for the estimation of Young’s modulus and Poisson’s ratio. This study combines the elastic impedance equation in terms of Young’s modulus and Poisson’s ratio and elastic impedance variation with incident angle inversion to produce a stable and direct way to estimate the Young’s modulus and Poisson’s ratio, with no need for density information from prestack seismic data. We initially derive a novel elastic impedance equation in terms of Young’s modulus and Poisson’s ratio. And then, to enhance the estimation stability, we develop the elastic impedance varying with incident angle inversion with damping singular value decomposition (EVA-DSVD) method to estimate the Young’s modulus and Poisson’s ratio. This method is implemented in a two-step inversion: Elastic impedance inversion and parameter estimation. The introduction of a model constraint and DSVD algorithm in parameter estimation renders the EVA-DSVD inversion more stable. Tests on synthetic data show that the Young’s modulus and Poisson’s ratio are still estimated reasonable with moderate noise. A test on a real data set shows that the estimated results are in good agreement with the results of well interpretation.

A skeptic's view of VVAz and AVAz

2020 ◽  
Vol 39 (2) ◽  
pp. 128-134 ◽  
Author(s):  
Norbert Van De Coevering ◽  
Klaas Koster ◽  
Rob Holt
Keyword(s):  
Phase Ii ◽  
Seismic Data ◽  
Wave Mode ◽  
P Wave ◽  
Noise Levels ◽  
Reservoir Properties ◽  
Data Set ◽  
Lower Amplitude ◽  
Software Note ◽  
Dense Sampling

We have applied a modern amplitude- and azimuth-preserving seismic data processing workflow to the SEG Advanced Modeling Program (SEAM) Phase II Barrett classic data set — an orthorhombic synthetic seismic model that has extremely dense sampling of all azimuths and offsets. We analyze the resulting prestack depth-migrated offset vector tiles with a variety of methods and software. Note that we only analyze the P-P wave mode, which is the focus of our study. We demonstrate that observed azimuthal changes cannot be correlated with the model's reservoir properties. We have made the migrated data available through SEAM. Compared to modeled data, real onshore seismic data have significantly lower amplitude fidelity, higher noise levels, and more uncertainty in the migration velocity field used for processing. Since we are unable to relate the anisotropy measured from the fully sampled clean SEAM Phase II Barrett synthetic seismic data to the model's known anisotropy, we conclude that it is highly unlikely that azimuthal variations observed on real onshore seismic data will be predictive of reservoir fracture properties.

Simultaneous inversion of PP and PS seismic data

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. R1-R10 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Martin Landrø
Keyword(s):  
Seismic Data ◽  
Reservoir Characterization ◽  
Model Building ◽  
Joint Inversion ◽  
P Wave ◽  
Inversion Method ◽  
System Of Equations ◽  
S Wave ◽  
Data Set ◽  
Simultaneous Inversion

Elastic parameters derived from seismic data are valuable input for reservoir characterization because they can be related to lithology and fluid content of the reservoir through empirical relationships. The relationship between physical properties of rocks and fluids and P-wave seismic data is nonunique. This leads to large uncertainties in reservoir models derived from P-wave seismic data. Because S- waves do not propagate through fluids, the combined use of P-and S-wave seismic data might increase our ability to derive fluid and lithology effects from seismic data, reducing the uncertainty in reservoir characterization and thereby improving 3D reservoir model-building. We present a joint inversion method for PP and PS seismic data by solving approximated linear expressions of PP and PS reflection coefficients simultaneously using a least-squares estimation algorithm. The resulting system of equations is solved by singular-value decomposition (SVD). By combining the two independent measurements (PP and PS seismic data), we stabilize the system of equations for PP and PS seismic data separately, leading to more robust parameter estimation. The method does not require any knowledge of PP and PS wavelets. We tested the stability of this joint inversion method on a 1D synthetic data set. We also applied the methodology to North Sea multicomponent field data to identify sand layers in a shallow formation. The identified sand layers from our inverted sections are consistent with observations from nearby well logs.

The structural architecture of the Whataroa Valley at the Alpine Fault (New Zealand) from first-arrival tomography and reflection imaging using an extended 3D VSP survey

2020 ◽  
Author(s):  
Vera Lay ◽  
Stefan Buske ◽  
Sascha Barbara Bodenburg ◽  
Franz Kleine ◽  
John Townend ◽  
...  
Keyword(s):  
New Zealand ◽  
Seismic Data ◽  
Average Velocity ◽  
Velocity Model ◽  
P Wave ◽  
3D Seismic ◽  
Data Set ◽  
Alpine Fault ◽  
P Wave Velocity ◽  
First Arrival

<p>The Alpine Fault along the West Coast of the South Island (New Zealand) is a major plate boundary that is expected to rupture in the next 50 years, likely as a magnitude 8 earthquake. The Deep Fault Drilling Project (DFDP) aims to deliver insight into the geological structure of this fault zone and its evolution by drilling and sampling the Alpine Fault at depth.  </p><p>Here we present results from a 3D seismic survey around the DFDP-2 drill site in the Whataroa Valley where the drillhole penetrated almost down to the fault surface. Within the glacial valley, we collected 3D seismic data to constrain valley structures that were obscured in previous 2D seismic data. The new data consist of a 3D extended vertical seismic profiling (VSP) survey using three-component receivers and a fibre optic cable in the DFDP-2B borehole as well as a variety of receivers at the surface.</p><p>The data set enables us to derive a reliable 3D P-wave velocity model by first-arrival travel time tomography. We identify a 100-460 m thick sediment layer (average velocity 2200±400 m/s) above the basement (average velocity 4200±500 m/s). Particularly on the western valley side, a region of high velocities steeply rises to the surface and mimics the topography. We interpret this to be the infilled flank of the glacial valley that has been eroded into the basement. In general, the 3D structures implied by the velocity model on the upthrown (Pacific Plate) side of the Alpine Fault correlate well with the surface topography and borehole findings.</p><p>A reliable velocity model is not only valuable by itself but it is also required as input for prestack depth migration (PSDM). We performed PSDM with a part of the 3D data set to derive a structural image of the subsurface within the Whataroa Valley. The top of the basement identified in the P-wave velocity model coincides well with reflectors in the migrated images so that we can analyse the geometry of the basement in detail.</p>

AVO inversion and poroelasticity with P- and S-wave moduli

Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. N17-N24 ◽  
Author(s):  
Zhaoyun Zong ◽  
Xingyao Yin ◽  
Guochen Wu
Keyword(s):  
Seismic Data ◽  
Real Data ◽  
P Wave ◽  
Model Parameters ◽  
Initial Model ◽  
S Wave ◽  
Data Set ◽  
Statistical Probability ◽  
Avo Inversion ◽  
Poroelasticity Theory

The fluid term in the Biot-Gassmann equation plays an important role in reservoir fluid discrimination. The density term imbedded in the fluid term, however, is difficult to estimate because it is less sensitive to seismic amplitude variations. We combined poroelasticity theory, amplitude variation with offset (AVO) inversion, and identification of P- and S-wave moduli to present a stable and physically meaningful method to estimate the fluid term, with no need for density information from prestack seismic data. We used poroelasticity theory to express the fluid term as a function of P- and S-wave moduli. The use of P- and S-wave moduli made the derivation physically meaningful and natural. Then we derived an AVO approximation in terms of these moduli, which can then be directly inverted from seismic data. Furthermore, this practical and robust AVO-inversion technique was developed in a Bayesian framework. The objective was to obtain the maximum a posteriori solution for the P-wave modulus, S-wave modulus, and density. Gaussian and Cauchy distributions were used for the likelihood and a priori probability distributions, respectively. The introduction of a low-frequency constraint and statistical probability information to the objective function rendered the inversion more stable and less sensitive to the initial model. Tests on synthetic data showed that all the parameters can be estimated well when no noise is present and the estimated P- and S-wave moduli were still reasonable with moderate noise and rather smooth initial model parameters. A test on a real data set showed that the estimated fluid term was in good agreement with the results of drilling.

Inversion and interpretation of a 3D seismic data set from the Ouachita Mountains, Oklahoma

Geophysics ◽  
2009 ◽  
Vol 74 (2) ◽  
pp. B37-B45 ◽  
Author(s):  
Abuduwali Aibaidula ◽  
George McMechan
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
3D Visualization ◽  
Wave Impedance ◽  
P Wave ◽  
3D Seismic ◽  
S Wave ◽  
Ouachita Mountains ◽  
Data Set ◽  
3D Seismic Data

Acoustic impedance inversion (AI) and simultaneous angle-dependent inversion (SADI) of a 3D seismic data set characterize reservoirs of Mississippian Morrowan age in the triangle zone of the frontal Ouachita Mountains, Oklahoma. Acoustic impedance of the near-angle seismic data images the 3D spatial distributions of Wapanucka limestone and Cromwell sandstone. Lamé [Formula: see text] ([Formula: see text] and [Formula: see text]) and [Formula: see text] sections are derived from the P-wave and S-wave impedance ([Formula: see text] and [Formula: see text]) sections produced by the SADI. Lithology is identified from the gamma logs and [Formula: see text]. The [Formula: see text], [Formula: see text], and [Formula: see text] are interpreted in terms of a hydrocarbon distribution pattern. The [Formula: see text] is used to identify high [Formula: see text] regions that are consistent with high sand/shale ratio. The estimated impedances and derived Lamé parameter sections are consistent with the interpretation that parts of the Wapanucka limestone and Cromwell sandstone contain potential gas reservoirs in fault-bounded compartments. The Cromwell sandstone contains the main inferred reservoirs; the two largest of these are each [Formula: see text] in pore volume. The inversion results also explain the observed low production in previous wells because those did not sample the best compartments. We propose a single new well location that would penetrate both reservoirs; 3D visualization facilitates this recommendation.

Sign in / Sign up

Export Citation Format

Share Document