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.

S-S imaging with vertical-force sources

2014 ◽  
Vol 2 (2) ◽  
pp. SE29-SE38 ◽  
Author(s):  
Bob A. Hardage ◽  
Donald Wagner
Keyword(s):  
Seismic Data ◽  
P Wave ◽  
Energy Sources ◽  
Minimal Amount ◽  
Equal Weight ◽  
Vertical Force ◽  
3D Seismic ◽  
S Wave ◽  
Receiver Arrays ◽  
3D Seismic Data

We show examples of S-S images created from multicomponent seismic data generated by vertical-force sources that can be quite useful to seismic interpreters. Two source types are used: vertical vibrators and shot-hole explosives. We first discuss S-S images made from data generated by a vertical vibrator and recorded with vertical receiver arrays of 3C geophones. We next show images extracted from surface-based 3C geophones deployed around this VSP well as a 3D seismic grid. The energy sources used to generate these surface 3D seismic data were shot-hole explosives. In all data examples, we observe that each type of vertical-force source (vertical vibrator and shot-hole explosive) produces abundant direct-S energy on radial and transverse geophones. We find only minimal amounts of P-wave energy on transverse-receiver data. In contrast, radial-receiver data have significant P-wave events intermingled with radial-S events. The minimal amount of P-wave noise on transverse-receiver data makes it easier to study S-S wave physics and to create S-S images with transverse-S data. The data examples focus on transverse-S data created by vertical-force sources because interpreters will find it more convenient to process and use this S-mode. Subsequent publications will assign equal weight to radial-S and transverse-S data.

Robust vector median filtering with a structure-adaptive implementation

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. V407-V414
Author(s):  
Yanghua Wang ◽  
Xiwu Liu ◽  
Fengxia Gao ◽  
Ying Rao
Keyword(s):  
Seismic Data ◽  
Impulse Noise ◽  
Reservoir Characterization ◽  
Median Filter ◽  
Original Data ◽  
Thin Layers ◽  
3D Seismic ◽  
Data Set ◽  
Vector Median ◽  
3D Seismic Data

The 3D seismic data in the prestack domain are contaminated by impulse noise. We have adopted a robust vector median filter (VMF) for attenuating the impulse noise from 3D seismic data cubes. The proposed filter has two attractive features. First, it is robust; the vector median that is the output of the filter not only has a minimum distance to all input data vectors, but it also has a high similarity to the original data vector. Second, it is structure adaptive; the filter is implemented following the local structure of coherent seismic events. The application of the robust and structure-adaptive VMF is demonstrated using an example data set acquired from an area with strong sedimentary rhythmites composed of steep-dipping thin layers. This robust filter significantly improves the signal-to-noise ratio of seismic data while preserving any discontinuity of reflections and maintaining the fidelity of amplitudes, which will facilitate the reservoir characterization that follows.

Application of sparse-layer inversion and harmonic bandwidth extension for a channel system in Southern Alberta, Canada

2020 ◽  
Vol 8 (2) ◽  
pp. T217-T229
Author(s):  
Yang Mu ◽  
John Castagna ◽  
Gabriel Gil
Keyword(s):  
Seismic Data ◽  
Reflection Coefficients ◽  
3D Seismic ◽  
Bandwidth Extension ◽  
Data Set ◽  
Frequency Spectra ◽  
Channel System ◽  
Southern Alberta ◽  
Reflectivity Inversion ◽  
3D Seismic Data

Sparse-layer reflectivity inversion decomposes a seismic trace into a limited number of simple layer responses and their corresponding reflection coefficients for top and base reflections. In contrast to sparse-spike inversion, the applied sparsity constraint is less biased against layer thickness and can thus better resolve thin subtuning layers. Application to a 3D seismic data set in Southern Alberta produces inverted impedances that have better temporal resolution and lateral stability and a less blocky appearance than sparse-spike inversion. Bandwidth extension harmonically extrapolated the frequency spectra of the inverted layers and nearly doubled the usable bandwidth. Although the prospective glauconitic sand tunes at approximately 37 m, bandwidth extension reduced the tuning thickness to 22 m. Bandwidth-extended data indicate a higher correlation with synthetic traces than the original seismic data and reveal features below the original tuning thickness. After bandwidth extension, the channel top and base are more evident on inline and crossline profiles. Lateral facies changes interpreted from the inverted acoustic impedance of the bandwidth-extended data are consistent with observations in wells.

Amplitude and frequency anomalies in regional 3D seismic data surrounding the Mallik 5L-38 research site, Mackenzie Delta, Northwest Territories, Canada

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. B183-B191 ◽  
Author(s):  
M. Riedel ◽  
G. Bellefleur ◽  
S. R. Dallimore ◽  
A. Taylor ◽  
J. F. Wright
Keyword(s):  
Northwest Territories ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Ratio Method ◽  
Unfrozen Water ◽  
3D Seismic ◽  
Mackenzie Delta ◽  
Data Set ◽  
Seismic Amplitude ◽  
3D Seismic Data

Amplitude and frequency anomalies associated with lakes and drainage systems were observed in a 3D seismic data set acquired in the Mallik area, Mackenzie Delta, Northwest Territories, Canada. The site is characterized by large gas hydrate deposits inferred from well-log analyses and coring. Regional interpretation of the gas hydrate occurrences is mainly based on seismic amplitude anomalies, such as brightening or blanking of seismic energy. Thus, the scope of this research is to understand the nature of the amplitude behavior in the seismic data. We have therefore analyzed the 3D seismic data to define areas with amplitude reduction due to contamination from lakes and channels and to distinguish them from areas where amplitude blanking may be a geologic signal. We have used the spectral ratio method to define attenuation (Q) over different areas in the 3D volume and subsequently applied Q-compensation to attenuate lateral variations ofdispersive absorption. Underneath larger lakes, seismic amplitude is reduced and the frequency content is reduced to [Formula: see text], which is half the original bandwidth. Traces with source-receiver pairs located inside of lakes show an attenuation factor Q of [Formula: see text], approximately half of that obtained for source-receiver pairs situated on deep, continuous permafrost outside of lakes. Deeper reflections occasionally identified underneath lakes show low-velocity-related pull-down. The vertical extent of the washout zones is enhanced by acquisition with limited offsets and from processing parameters such as harsh mute functions to reduce noise from surface waves. The strong attenuation and seismic pull-down may indicate the presence of unfrozen water in deeper lakes and unfrozen pore water within the sediments underlying the lakes. Thus, the blanking underneath lakes is not necessarily related to gas migration or other in situ changes in physical properties potentially associated with the presence of gas hydrate.

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>

P‐wave anisotropy from azimuthal velocity and AVO using wide‐azimuth 3D seismic data

2004 ◽  
Author(s):  
Ahmed M. Al‐Marzoug ◽  
Fernando A. Neves ◽  
Jung J. Kim ◽  
Edgardo L. Nebrija
Keyword(s):  
Seismic Data ◽  
Azimuthal Velocity ◽  
P Wave ◽  
3D Seismic ◽  
3D Seismic Data

Bayesian linearized AVO inversion

Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 185-198 ◽  
Author(s):  
Arild Buland ◽  
Henning Omre
Keyword(s):  
Wave Velocity ◽  
Posterior Distribution ◽  
Acoustic Impedance ◽  
Velocity Ratio ◽  
P Wave ◽  
Inversion Method ◽  
Inversion Algorithm ◽  
S Wave ◽  
Data Set ◽  
Avo Inversion

A new linearized AVO inversion technique is developed in a Bayesian framework. The objective is to obtain posterior distributions for P‐wave velocity, S‐wave velocity, and density. Distributions for other elastic parameters can also be assessed—for example, acoustic impedance, shear impedance, and P‐wave to S‐wave velocity ratio. The inversion algorithm is based on the convolutional model and a linearized weak contrast approximation of the Zoeppritz equation. The solution is represented by a Gaussian posterior distribution with explicit expressions for the posterior expectation and covariance; hence, exact prediction intervals for the inverted parameters can be computed under the specified model. The explicit analytical form of the posterior distribution provides a computationally fast inversion method. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved when the noise approached zero. With realistic noise levels, acoustic impedance was the best determined parameter, while the inversion provided practically no information about the density. The inversion algorithm has also been tested on a real 3‐D data set from the Sleipner field. The results show good agreement with well logs, but the uncertainty is high.

scholarly journals Use of Attributes on 3D Seismic Data for Studying Mishrif Formation in East Abu-Amoud Field, South-Eastern Iraq

2021 ◽  
pp. 4802-4809
Author(s):  
Mohammed H. Al-Aaraji ◽  
Hussein H. Karim
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Seismic Attributes ◽  
3D Seismic ◽  
Seismic Section ◽  
Reflected Waves ◽  
South Eastern ◽  
Seismic Sections ◽  
Sonic Logs ◽  
3D Seismic Data

      The seismic method depends on the nature of the reflected waves from the interfaces between layers, which in turn depends on the density and velocity of the layer, and this is called acoustic impedance. The seismic sections of the East Abu-Amoud field that is located in Missan Province, south-eastern Iraq, were studied and interpreted for updating the structural picture of the major Mishrif Formation for the reservoir in the field. The Mishrif Formation is rich in petroleum in this area, with an area covering about 820 km2. The horizon was calibrated and defined on the seismic section with well logs data (well tops, check shot, sonic logs, and density logs) in the interpretation process to identify the upper and lower boundaries of the Formation.  Seismic attributes were used to study the formation, including instantaneous phase attributes and relative acoustic impedance on time slice of 3D seismic data . Also, relative acoustic impedance was utilized to study the top of the Mishrif Formation. Based on these seismic attributes, karst features of the formation were identified. In addition, the nature of the lithology in the study area and the change in porosity were determined through the relative acoustic impedance The overlap of the top of the Mishrif Formation with the bottom of the Khasib Formation was determined because the Mishrif Formation is considered as an unconformity surface.

Sign in / Sign up

Export Citation Format

Share Document