Recent advances in seismic lithologic analysis

Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 42-46 ◽  
Author(s):  
John P. Castagna
Keyword(s):  
Seismic Data ◽  
Recent Progress ◽  
Seismic Analysis ◽  
Rock Physics ◽  
Pore Fluid ◽  
Analysis Tool ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Fluid Content ◽  
Key Aspects

An objective of seismic analysis is to quantitatively extract lithology, porosity, and pore fluid content directly from seismic data. Rock physics provides the fundamental basis for seismic lithology determination. Beyond conventional poststack inversion, the most important seismic lithologic analysis tool is amplitude‐variation‐with‐offset (AVO) analysis. In this paper, I review recent progress in these two key aspects of seismic lithologic analysis.

Seismic reservoir characterization of a class-1 amplitude variation with offset turbiditic system located offshore Cote d’Ivoire, West Africa

2018 ◽  
Vol 6 (2) ◽  
pp. SD115-SD128
Author(s):  
Pedro Alvarez ◽  
William Marin ◽  
Juan Berrizbeitia ◽  
Paola Newton ◽  
Michael Barrett ◽  
...  
Keyword(s):  
West Africa ◽  
Seismic Data ◽  
Reservoir Characterization ◽  
Rock Physics ◽  
Rock Properties ◽  
Well Log ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Fluid Content ◽  
Class 1

We have evaluated a case study, in which a class-1 amplitude variation with offset (AVO) turbiditic system located offshore Cote d’Ivoire, West Africa, is characterized in terms of rock properties (lithology, porosity, and fluid content) and stratigraphic elements using well-log and prestack seismic data. The methodology applied involves (1) the conditioning and modeling of well-log data to several plausible geologic scenarios at the prospect location, (2) the conditioning and inversion of prestack seismic data for P- and S-wave impedance estimation, and (3) the quantitative estimation of rock property volumes and their geologic interpretation. The approaches used for the quantitative interpretation of these rock properties were the multiattribute rotation scheme for lithology and porosity characterization and a Bayesian lithofluid facies classification (statistical rock physics) for a probabilistic evaluation of fluid content. The result indicates how the application and integration of these different AVO- and rock-physics-based reservoir characterization workflows help us to understand key geologic stratigraphic elements of the architecture of the turbidite system and its static petrophysical characteristics (e.g., lithology, porosity, and net sand thickness). Furthermore, we found out how to quantify and interpret the risk related to the probability of finding hydrocarbon in a class-1 AVO setting using seismically derived elastic attributes, which are characterized by having a small level of sensitivity to changes in fluid saturation.

An integrated broadband preprocessing method for towed-streamer seismic data

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. V201-V221 ◽  
Author(s):  
Mehdi Aharchaou ◽  
Erik Neumann
Keyword(s):  
Seismic Data ◽  
High Frequency ◽  
Pressure Data ◽  
Amplitude Variation ◽  
Field Source ◽  
Amplitude Variation With Offset ◽  
Ringing Artifacts ◽  
Joint Application ◽  
Amplitude Compensation ◽  
Unified Method

Broadband preprocessing has become widely used for marine towed-streamer seismic data. In the standard workflow, far-field source designature, receiver and source-side deghosting, and redatuming to mean sea level are applied in sequence, with amplitude compensation for background [Formula: see text] delayed until the imaging or postmigration stages. Thus, each step is likely to generate its own artifacts, quality checking can be time-consuming, and broadband data are only obtained late in this chained workflow. We have developed a unified method for broadband preprocessing — called integrated broadband preprocessing (IBP) — which enables the joint application of all the above listed steps early in the processing sequence. The amplitude, phase, and amplitude-variation-with-offset fidelity of IBP are demonstrated on pressure data from the shallow, deep, and slanted streamers. The integration allows greater sparsity to emerge in the representation of seismic data, conferring clear benefits over the sequential application. Moreover, time sparsity, full dimensionality, and early amplitude [Formula: see text] compensation all have an impact on broadband data quality, in terms of reduced ringing artifacts, improved wavelet integrity at large crossline angles, and fewer residual high-frequency multiples.

Value of information analysis and Bayesian inversion for closed skew-normal distributions: Applications to seismic amplitude variation with offset data

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. R151-R163 ◽  
Author(s):  
Javad Rezaie ◽  
Jo Eidsvik ◽  
Tapan Mukerji
Keyword(s):  
Seismic Data ◽  
Value Of Information ◽  
Bayesian Inversion ◽  
Information Analysis ◽  
Amplitude Variation ◽  
Mean Values ◽  
Amplitude Variation With Offset ◽  
Seismic Amplitude ◽  
Value Of Information Analysis ◽  
Skew Normal

Information analysis can be used in the context of reservoir decisions under uncertainty to evaluate whether additional data (e.g., seismic data) are likely to be useful in impacting the decision. Such evaluation of geophysical information sources depends on input modeling assumptions. We studied results for Bayesian inversion and value of information analysis when the input distributions are skewed and non-Gaussian. Reservoir parameters and seismic amplitudes are often skewed and using models that capture the skewness of distributions, the input assumptions are less restrictive and the results are more reliable. We examined the general methodology for value of information analysis using closed skew normal (SN) distributions. As an example, we found a numerical case with porosity and saturation as reservoir variables and computed the value of information for seismic amplitude variation with offset intercept and gradient, all modeled with closed SN distributions. Sensitivity of the value of information analysis to skewness, mean values, accuracy, and correlation parameters is performed. Simulation results showed that fewer degrees of freedom in the reservoir model results in higher value of information, and seismic data are less valuable when seismic measurements are spatially correlated. In our test, the value of information was approximately eight times larger for a spatial-dependent reservoir variable compared with the independent case.

Traveltime and relative geometric spreading approximation in elastic orthorhombic medium

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. C153-C162 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas ◽  
Hitoshi Mikada
Keyword(s):  
Seismic Data ◽  
Real Data ◽  
P Wave ◽  
Form Expression ◽  
Amplitude Variation ◽  
Rational Form ◽  
Amplitude Variation With Offset ◽  
Practical Applications ◽  
Numerical Tests ◽  
Geometric Spreading

Wavefield properties such as traveltime and relative geometric spreading (traveltime derivatives) are highly essential in seismic data processing and can be used in stacking, time-domain migration, and amplitude variation with offset analysis. Due to the complexity of an elastic orthorhombic (ORT) medium, analysis of these properties becomes reasonably difficult, where accurate explicit-form approximations are highly recommended. We have defined the shifted hyperbola form, Taylor series (TS), and the rational form (RF) approximations for P-wave traveltime and relative geometric spreading in an elastic ORT model. Because the parametric form expression for the P-wave vertical slowness in the derivation is too complicated, TS (expansion in offset) is applied to facilitate the derivation of approximate coefficients. The same approximation forms computed in the acoustic ORT model also are derived for comparison. In the numerical tests, three ORT models with parameters obtained from real data are used to test the accuracy of each approximation. The numerical examples yield results in which, apart from the error along the y-axis in ORT model 2 for the relative geometric spreading, the RF approximations all are very accurate for all of the tested models in practical applications.

Processing, inversion, and interpretation of a 2D seismic data set from the North Viking Graben, North Sea

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 837-848 ◽  
Author(s):  
Gislain B. Madiba ◽  
George A. McMechan
Keyword(s):  
Seismic Data ◽  
Dynamic Range ◽  
Velocity Ratio ◽  
Total Porosity ◽  
Shear Velocity ◽  
Seismic Section ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Fluid Content ◽  
The North

Simultaneous elastic impedance inversion is performed on the 2D North Viking Graben seismic data set used at the 1994 SEG workshop on amplitude variation with offset and inversion. P‐velocity (Vp), S‐velocity (Vs), density logs, and seismic data are input to the inversion. The inverted P‐impedance and S‐impedance sections are used to generate an approximate compressional‐to‐shear velocity ratio (Vp/Vs) section which, in turn, is used along with water‐filled porosity (Swv) derived from the logs from two wells, to generate fluid estimate sections. This is possible as the reservoir sands have fairly constant total porosity of approximately 28 ± 4%, so the hydrocarbon filled porosity is the total porosity minus the water‐filled porosity. To enhance the separation of lithologies and of fluid content, we map Vp/Vs into Swv using an empirical crossplot‐derived relation. This mapping expands the dynamic range of the low end of the Vp/Vs values. The different lithologies and fluids are generally well separated in the Vp/Vs–Swv domain. Potential hydrocarbon reservoirs (as calibrated by the well data) are identified throughout the seismic section and are consistent with the fluid content estimations obtained from alternative computations. The Vp/Vs–Swv plane still does not produce unique interpretation in many situations. However, the critical distinction, which is between hydrocarbon‐bearing sands and all other geologic/reservoir configurations, is defined. Swv ≤ 0.17 and Vp/Vs ≤ 1.8 are the criteria that delineate potential reservoirs in this area, with decreasing Swv indicating a higher gas/oil ratio, and decreasing Vp/Vs indicating a higher sand/shale ratio. As these criteria are locally calibrated, they appear to be valid locally; they should not be applied to other data sets, which may exhibit significantly different relationships. However, the overall procedure should be generally applicable.

Fracture detection via correlating P-wave amplitude variation with offset and azimuth analysis and well data in eastern central Saudi Arabia

2017 ◽  
Vol 5 (4) ◽  
pp. T531-T544
Author(s):  
Ali H. Al-Gawas ◽  
Abdullatif A. Al-Shuhail
Keyword(s):  
Saudi Arabia ◽  
Seismic Data ◽  
Production Data ◽  
Fracture Zones ◽  
Fracture Detection ◽  
Amplitude Variation ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Well Data ◽  
Patch Geometry

The late Carboniferous clastic Unayzah-C in eastern central Saudi Arabia is a low-porosity, possibly fractured reservoir. Mapping the Unayzah-C is a challenge due to the low signal-to-noise ratio (S/N) and limited bandwidth in the conventional 3D seismic data. A related challenge is delineating and characterizing fracture zones within the Unayzah-C. Full-azimuth 3D broadband seismic data were acquired using point receivers, low-frequency sweeps down to 2 Hz, and 6 km patch geometry. The data indicate significant enhancement in continuity and resolution of the reflection data, leading to improved mapping of the Unayzah-C. Because the data set has a rectangular patch geometry with full inline offsets to 6000 m, using amplitude variation with offset and azimuth (AVOA) may be effective to delineate and characterize fracture zones within Unayzah-A and Unayzah-C. The study was undertaken to determine the improvement of wide-azimuth seismic data in fracture detection in clastic reservoirs. The results were validated with available well data including borehole images, well tests, and production data in the Unayzah-A. There are no production data or borehole images within the Unayzah-C. For validation, we had to refer to a comparison of alternative seismic fracture detection methods, mainly curvature and coherence. Anisotropy was found to be weak, which may be due to noise, clastic lithology, and heterogeneity of the reservoirs, in both reservoirs except for along the western steep flank of the study area. These may correspond to some north–south-trending faults suggested by circulation loss and borehole image data in a few wells. The orientation of the long axis of the anisotropy ellipses is northwest–southeast, and it is not in agreement with the north–south structural trend. No correlation was found among the curvature, coherence, and AVOA in Unayzah-A or Unayzah-C. Some possible explanations for the low correlation between the AVOA ellipticity and the natural fractures are a noisy data set, overburden anisotropy, heterogeneity, granulation seams, and deformation.

scholarly journals AVO analysis for high amplitude anomalies using 2D pre-stack seismic data

2022 ◽  
Author(s):  
Lamees N. Abdulkareem ◽  
Keyword(s):  
Seismic Data ◽  
High Amplitude ◽  
Rock Properties ◽  
Fluid Substitution ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
The Real ◽  
North West ◽  
Avo Analysis ◽  
Controlling Parameter

Amplitude variation with offset (AVO) analysis is an 1 efficient tool for hydrocarbon detection and identification of elastic rock properties and fluid types. It has been applied in the present study using reprocessed pre-stack 2D seismic data (1992, Caulerpa) from north-west of the Bonaparte Basin, Australia. The AVO response along the 2D pre-stack seismic data in the Laminaria High NW shelf of Australia was also investigated. Three hypotheses were suggested to investigate the AVO behaviour of the amplitude anomalies in which three different factors; fluid substitution, porosity and thickness (Wedge model) were tested. The AVO models with the synthetic gathers were analysed using log information to find which of these is the controlling parameter on the AVO analysis. AVO cross plots from the real pre-stack seismic data reveal AVO class IV (showing a negative intercept decreasing with offset). This result matches our modelled result of fluid substitution for the seismic synthetics. It is concluded that fluid substitution is the controlling parameter on the AVO analysis and therefore, the high amplitude anomaly on the seabed and the target horizon 9 is the result of changing the fluid content and the lithology along the target horizons. While changing the porosity has little effect on the amplitude variation with offset within the AVO cross plot. Finally, results from the wedge models show that a small change of thickness causes a change in the amplitude; however, this change in thickness gives a different AVO characteristic and a mismatch with the AVO result of the real 2D pre-stack seismic data. Therefore, a constant thin layer with changing fluids is more likely to be the cause of the high amplitude anomalies.

The Shuey-Radon transform

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. V197-V206 ◽  
Author(s):  
Ali Gholami ◽  
Milad Farshad
Keyword(s):  
Radon Transform ◽  
Seismic Data ◽  
Matching Pursuit ◽  
Real Data ◽  
Transform Domain ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Linear System Of Equations ◽  
Matching Pursuit Algorithm ◽  
Zero Offset

The traditional hyperbolic Radon transform (RT) decomposes seismic data into a sum of constant amplitude basis functions. This limits the performance of the transform when dealing with real data in which the reflection amplitudes include the amplitude variation with offset (AVO) variations. We adopted the Shuey-Radon transform as a combination of the RT and Shuey’s approximation of reflectivity to accurately model reflections including AVO effects. The new transform splits the seismic gather into three Radon panels: The first models the reflections at zero offset, and the other two panels add capability to model the AVO gradient and curvature. There are two main advantages of the Shuey-Radon transform over similar algorithms, which are based on a polynomial expansion of the AVO response. (1) It is able to model reflections more accurately. This leads to more focused coefficients in the transform domain and hence provides more accurate processing results. (2) Unlike polynomial-based approaches, the coefficients of the Shuey-Radon transform are directly connected to the classic AVO parameters (intercept, gradient, and curvature). Therefore, the resulting coefficients can further be used for interpretation purposes. The solution of the new transform is defined via an underdetermined linear system of equations. It is formulated as a sparsity-promoting optimization, and it is solved efficiently using an orthogonal matching pursuit algorithm. Applications to different numerical experiments indicate that the Shuey-Radon transform outperforms the polynomial and conventional RTs.

A successful qualitative and quantitative integrated interpretation of Lower Miocene wells and seismic data in Salina del Istmo Basin, Mexico

2021 ◽  
Vol 40 (12) ◽  
pp. 897-904
Author(s):  
Manuel González-Quijano ◽  
Gregor Baechle ◽  
Miguel Yanez ◽  
Freddy Obregon ◽  
Carmen Vito ◽  
...  
Keyword(s):  
Seismic Data ◽  
Rock Physics ◽  
Reservoir Rock ◽  
Amplitude Variation With Offset ◽  
Volcanic Material ◽  
Qualitative And Quantitative ◽  
Rock Types ◽  
Deep Waters ◽  
Well Data ◽  
Gulf Of Mexico Basin

The study area is located in middepth to deep waters of the Salina del Istmo Basin where Repsol operates Block 29. The objective of this work is to integrate qualitative and quantitative interpretations of rock and seismic data to predict lithology and fluid of the Early Miocene prospects. The seismic expression of those prospects differs from age-equivalent well-studied analog fields in the U.S. Gulf of Mexico Basin due to the mineralogically complex composition of abundant extrusive volcanic material. Offset well data (i.e., core, logs, and cuttings) were used to discriminate lithology types and to quantify mineralogy. This analysis served as input for developing a new rock-physics framework and performing amplitude variation with offset (AVO) modeling. The results indicate that the combination of intercept and gradient makes it possible to discriminate hydrocarbon-filled (AVO class II and III) versus nonhydrocarbon-filled rocks (AVO class 0 and IV). Different lithologies within hydrocarbon-bearing reservoirs cannot be discriminated as the gradient remains negative for all rock types. However, AVO analysis allows discrimination of three different reservoir rock types in water-bearing cases (AVO class 0, I, and IV). These conclusions were obtained during studies conducted in 2018–2019 and were used in prospect evaluation to select drilling locations leading to two wildcat discoveries, the Polok and Chinwol prospects, drilled in Block 29 in 2020.

Anisotropic 3D elastic full-wavefield inversion to directly estimate elastic properties and its role in interpretation

2021 ◽  
Vol 40 (4) ◽  
pp. 277-286
Author(s):  
Haiyang Wang ◽  
Olivier Burtz ◽  
Partha Routh ◽  
Don Wang ◽  
Jake Violet ◽  
...  
Keyword(s):  
Elastic Properties ◽  
Seismic Data ◽  
Reservoir Characterization ◽  
P Wave ◽  
Quality Estimation ◽  
Amplitude Variation ◽  
Linear Inversion ◽  
Amplitude Variation With Offset ◽  
P Wave Velocity ◽  
Wavefield Inversion

Elastic properties from seismic data are important to determine subsurface hydrocarbon presence and have become increasingly important for detailed reservoir characterization that aids to derisk specific hydrocarbon prospects. Traditional techniques to extract elastic properties from seismic data typically use linear inversion of imaged products (migrated angle stacks). In this research, we attempt to get closer to Tarantola's visionary goal for full-wavefield inversion (FWI) by directly obtaining 3D elastic properties from seismic shot-gather data with limited well information. First, we present a realistic 2D synthetic example to show the need for elastic physics in a strongly elastic medium. Then, a 3D field example from deepwater West Africa is used to validate our workflow, which can be practically used in today's computing architecture. To enable reservoir characterization, we produce elastic products in a cascaded manner and run 3D elastic FWI up to 50 Hz. We demonstrate that reliable and high-resolution P-wave velocity can be retrieved in a strongly elastic setting (i.e., with a class 2 or 2P amplitude variation with offset response) in addition to higher-quality estimation of P-impedance and VP/VS ratio. These parameters can be directly used in interpretation, lithology, and fluid prediction.

Sign in / Sign up

Export Citation Format

Share Document