Characterization of layered anisotropic media from prestack PS-wave-reflection data

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. D171-D182 ◽  
Author(s):  
Jason E. Gumble ◽  
James E. Gaiser
Keyword(s):  
Gulf Of Mexico ◽  
Seismic Data ◽  
Equivalent Model ◽  
Time Shift ◽  
Layer By Layer ◽  
Converted Wave ◽  
Data Set ◽  
Fracture Characterization ◽  
Anisotropic Models ◽  
Hti Media

Anisotropy and fracture characterization in individual layers is realized through iterative layer stripping corrections of four, converted-wave (PS-wave) synthetic reflection seismic data sets, generated from azimuthally anisotropic (HTI and TTI) models, and a four component (4-C) data set from the Teal South, Gulf of Mexico. The corrections were applied on a layer-by-layer basis to evaluate the efficacy of constant polarization rotation and time-shift operators. Equivalent isotropic models were compared to anisotropic models after layer-stripping corrections using rms amplitude and shear-wave-splitting time-difference maps to quantify and identify inherent errors in estimating seismic polarization parameters. For HTI media radial and transverse components of PS data that have had layer-stripping corrections applied, exhibit incorrect symmetry and orientations. This may adversely affect inversion and/or amplitude-variation with angle offset (AVO) and amplitude versus azimuth (AVA)analysis. Layer-stripping corrections applied to fast and slow ([Formula: see text] and [Formula: see text], respectively) components exhibit the correct symmetry and orientation. Time differences between PS1 and PS2 are computed using crosscorrelation. Previous studies have addressed some of the problems associated with layer-stripping corrections for the case of vertical fractures (HTI media) and poststack layer-stripping analyses. This study includes an equivalent model with dipping fractures (TTI media) and extends the scope to encompass the effects of anisotropy on prestack data. The results from an application of the same technique are also applied to a limited set of 4-C data from the Teal South project in the Gulf of Mexico. Results are consistent with those of previous studies involving solely poststack 4-C rotation analysis in terms of average, or zero offset, time differences and symmetry orientation. Offset and azimuth amplitude/traveltime variations, however, indicate that there is more information contained in prestack seismic data than 4-C rotation can comprehend.

scholarly journals Converted-wave azimuth moveout

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. S99-S110
Author(s):  
Daniel A. Rosales ◽  
Biondo Biondi
Keyword(s):  
Seismic Data ◽  
Extreme Case ◽  
Ocean Bottom ◽  
Sampled Data ◽  
Computationally Efficient ◽  
Converted Wave ◽  
Data Set ◽  
Prestack Migration ◽  
Conversion Point ◽  
The Common

A new partial-prestack migration operator to manipulate multicomponent data, called converted-wave azimuth moveout (PS-AMO), transforms converted-wave prestack data with an arbitrary offset and azimuth to equivalent data with a new offset and azimuth position. This operator is a sequential application of converted-wave dip moveout and its inverse. As expected, PS-AMO reduces to the known expression of AMO for the extreme case when the P velocity is the same as the S velocity. Moreover, PS-AMO preserves the resolution of dipping events and internally applies a correction for the lateral shift between the common-midpoint and the common-reflection/conversion point. An implementation of PS-AMO in the log-stretch frequency-wavenumber domain is computationally efficient. The main applications for the PS-AMO operator are geometry regularization, data-reduction through partial stacking, and interpolation of unevenly sampled data. We test our PS-AMO operator by solving 3D acquisition geometry-regularization problems for multicomponent, ocean-bottom seismic data. The geometry-regularization problem is defined as a regularized least-squares-objective function. To preserve the resolution of dipping events, the regularization term uses the PS-AMO operator. Application of this methodology on a portion of the Alba 3D, multicomponent, ocean-bottom seismic data set shows that we can satisfactorily obtain an interpolated data set that honors the physics of converted waves.

Time-lapse acoustic impedance variations during CO2 injection in Weyburn oilfield, Canada

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. M1-M13 ◽  
Author(s):  
Yichuan Wang ◽  
Igor B. Morozov
Keyword(s):  
Oil Recovery ◽  
Seismic Data ◽  
Acoustic Impedance ◽  
Time Lapse ◽  
Time Shift ◽  
Co2 Injection ◽  
Inversion Method ◽  
Data Sets ◽  
Storage Site ◽  
Data Set

For seismic monitoring injected fluids during enhanced oil recovery or geologic [Formula: see text] sequestration, it is useful to measure time-lapse (TL) variations of acoustic impedance (AI). AI gives direct connections to the mechanical and fluid-related properties of the reservoir or [Formula: see text] storage site; however, evaluation of its subtle TL variations is complicated by the low-frequency and scaling uncertainties of this attribute. We have developed three enhancements of TL AI analysis to resolve these issues. First, following waveform calibration (cross-equalization) of the monitor seismic data sets to the baseline one, the reflectivity difference was evaluated from the attributes measured during the calibration. Second, a robust approach to AI inversion was applied to the baseline data set, based on calibration of the records by using the well-log data and spatially variant stacking and interval velocities derived during seismic data processing. This inversion method is straightforward and does not require subjective selections of parameterization and regularization schemes. Unlike joint or statistical inverse approaches, this method does not require prior models and produces accurate fitting of the observed reflectivity. Third, the TL AI difference is obtained directly from the baseline AI and reflectivity difference but without the uncertainty-prone subtraction of AI volumes from different seismic vintages. The above approaches are applied to TL data sets from the Weyburn [Formula: see text] sequestration project in southern Saskatchewan, Canada. High-quality baseline and TL AI-difference volumes are obtained. TL variations within the reservoir zone are observed in the calibration time-shift, reflectivity-difference, and AI-difference images, which are interpreted as being related to the [Formula: see text] injection.

Fully automated near-surface analysis by surface-consistent refraction method

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U39-U49 ◽  
Author(s):  
Daniele Colombo ◽  
Federico Miorelli ◽  
Ernesto Sandoval ◽  
Kevin Erickson
Keyword(s):  
Surface Analysis ◽  
Seismic Data ◽  
Linear Equations ◽  
Time Shift ◽  
Data Sets ◽  
Time Lags ◽  
Data Set ◽  
Near Surface ◽  
Seismic Acquisition ◽  
Key Aspects

Industry practices for near-surface analysis indicate difficulties in coping with the increased number of channels in seismic acquisition systems, and new approaches are needed to fully exploit the resolution embedded in modern seismic data sets. To achieve this goal, we have developed a novel surface-consistent refraction analysis method for low-relief geology to automatically derive near-surface corrections for seismic data processing. The method uses concepts from surface-consistent analysis applied to refracted arrivals. The key aspects of the method consist of the use of common midpoint (CMP)-offset-azimuth binning, evaluation of mean traveltime and standard deviation for each bin, rejection of anomalous first-break (FB) picks, derivation of CMP-based traveltime-offset functions, conversion to velocity-depth functions, evaluation of long-wavelength statics, and calculation of surface-consistent residual statics through waveform crosscorrelation. Residual time lags are evaluated in multiple CMP-offset-azimuth bins by crosscorrelating a pilot trace with all the other traces in the gather in which the correlation window is centered at the refracted arrival. The residuals are then used to build a system of linear equations that is simultaneously inverted for surface-consistent shot and receiver time shift corrections plus a possible subsurface residual term. All the steps are completely automated and require a fraction of the time needed for conventional near-surface analysis. The developed methodology was successfully performed on a complex 3D land data set from Central Saudi Arabia where it was benchmarked against a conventional tomographic work flow. The results indicate that the new surface-consistent refraction statics method enhances seismic imaging especially in portions of the survey dominated by noise.

Building tilted transversely isotropic depth models using localized anisotropic tomography with well information

Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. D27-D36 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Synthetic Data ◽  
Rock Physics ◽  
Transversely Isotropic ◽  
Data Set ◽  
Anisotropic Models ◽  
Well Data

Tilted transverse isotropy (TTI) is increasingly recognized as a more geologically plausible description of anisotropy in sedimentary formations than vertical transverse isotropy (VTI). Although model-building approaches for VTI media are well understood, similar approaches for TTI media are in their infancy, even when the symmetry-axis direction is assumed known. We describe a tomographic approach that builds localized anisotropic models by jointly inverting surface-seismic and well data. We present a synthetic data example of anisotropic tomography applied to a layered TTI model with a symmetry-axis tilt of 45 degrees. We demonstrate three scenarios for constraining the solution. In the first scenario, velocity along the symmetry axis is known and tomography inverts for Thomsen’s [Formula: see text] and [Formula: see text] parame-ters. In the second scenario, tomography inverts for [Formula: see text], [Formula: see text], and velocity, using surface-seismic data and vertical check-shot traveltimes. In contrast to the VTI case, both these inversions are nonunique. To combat nonuniqueness, in the third scenario, we supplement check-shot and seismic data with the [Formula: see text] profile from an offset well. This allows recovery of the correct profiles for velocity along the symmetry axis and [Formula: see text]. We conclude that TTI is more ambiguous than VTI for model building. Additional well data or rock-physics assumptions may be required to constrain the tomography and arrive at geologically plausible TTI models. Furthermore, we demonstrate that VTI models with atypical Thomsen parameters can also fit the same joint seismic and check-shot data set. In this case, although imaging with VTI models can focus the TTI data and match vertical event depths, it leads to substantial lateral mispositioning of the reflections.

Focusing analysis of seismic data with peg‐leg multiples

Geophysics ◽  
1997 ◽  
Vol 62 (1) ◽  
pp. 177-182
Author(s):  
Einar Maeland
Keyword(s):  
Seismic Data ◽  
Synthetic Data ◽  
Velocity Estimation ◽  
Time Shift ◽  
Data Set ◽  
Seismic Migration ◽  
Imaging Condition ◽  
Multilayered Medium ◽  
Velocity Information ◽  
Zero Offset

Seismic migration with an erroneous velocity field produces an “image” that must be interpreted to obtain reliable velocity information. Inclusion of multiples in common‐shot data makes velocity estimation more difficult. Zero‐offset migration, based on the exploding reflector model, can be used to identify peg‐leg multiples by application of an extra time shift in the imaging condition. The time and the corresponding position when reflected energy focuses must be detected by inspection of the migrated data set. Formulas are derived and the method is tested on synthetic data from a multilayered medium.

The effect of geologic structure on fracture characterization using azimuthal VSPs

Geophysics ◽  
2020 ◽  
pp. 1-62
Author(s):  
Ali Sayed ◽  
Robert R. Stewart ◽  
Dhananjay Kumar
Keyword(s):  
Transverse Energy ◽  
Structural Model ◽  
Transverse Isotropy ◽  
Structural Response ◽  
Fracture Density ◽  
Geologic Structure ◽  
Converted Wave ◽  
Fracture Characterization ◽  
Flat Interface ◽  
Hti Media

Azimuthal VSP (AzVSP) surveys have been commonly used for fracture characterization by analyzing the P-to-S converted wave response across fractured zones. Fractured media produce characteristic two-cycle patterns on AzVSP gathers that are disrupted in the presence of complex structures. Aiming to characterize the complexity of AzVSP response in the presence of geologic structure, we derive anisotropy parameters for effective horizontal transverse isotropy (HTI) media, generate AzVSP signatures for the flat-interface case as the baseline, and compare the baseline signatures to AzVSP signatures for the dipping-interface cases. Fracture fill and fracture intensity are varied to capture the effect of intrinsic fracture parameters on AzVSP signatures. We find that dry fractures produce a stronger response compared to fluid-filled fractures for the same fracture density. The structural response on AzVSP signatures is isolated by generating synthetic seismograms across structurally equivalent isotropic models. Transverse energy, resulting from structures, could be misinterpreted as evidence for fracturing. AzVSP signatures for a fractured dipping-interface two-layer model show significant distortion of the fracture response for the 10° dip case as compared to the flat-interface signatures. As a possible solution to the structural problem, we use the arrival azimuths of the transmitted P-to-S event in synthetic signatures generated with the isotropic structural model, to orient the radial and transverse components parallel and perpendicular to the isotropic P-to-S event, respectively. Such structurally consistent orientation negates the effect of structure on azimuthal VSP gathers and uncovers the underlying fracture response. This methodology can be extended to complex overburdens that are structurally well constrained.

Localized anisotropic tomography with well information in VTI media

Geophysics ◽  
2010 ◽  
Vol 75 (5) ◽  
pp. D37-D45 ◽  
Author(s):  
Andrey Bakulin ◽  
Marta Woodward ◽  
Dave Nichols ◽  
Konstantin Osypov ◽  
Olga Zdraveva
Keyword(s):  
Vertical Velocity ◽  
Seismic Data ◽  
Model Building ◽  
Synthetic Data ◽  
Resolution Limit ◽  
Data Set ◽  
Anisotropic Models ◽  
Velocity Models ◽  
First Case ◽  
Well Data

We develop a concept of localized seismic grid tomography constrained by well information and apply it to building vertically transversely isotropic (VTI) velocity models in depth. The goal is to use a highly automated migration velocity analysis to build anisotropic models that combine optimal image focusing with accurate depth positioning in one step. We localize tomography to a limited volume around the well and jointly invert the surface seismic and well data. Well information is propagated into the local volume by using the method of preconditioning, whereby model updates are shaped to follow geologic layers with spatial smoothing constraints. We analyze our concept with a synthetic data example of anisotropic tomography applied to a 1D VTI model. We demonstrate four cases of introducing additionalinformation. In the first case, vertical velocity is assumed to be known, and the tomography inverts only for Thomsen’s [Formula: see text] and [Formula: see text] profiles using surface seismic data alone. In the second case, tomography simultaneously inverts for all three VTI parameters, including vertical velocity, using a joint data set that consists of surface seismic data and vertical check-shot traveltimes. In the third and fourth cases, sparse depth markers and walkaway vertical seismic profiling (VSP) are used, respectively, to supplement the seismic data. For all four examples, tomography reliably recovers the anisotropic velocity field up to a vertical resolution comparable to that of the well data. Even though walkaway VSP has the additional dimension of angle or offset, it offers no further increase in this resolution limit. Anisotropic tomography with well constraints has multiple advantages over other approaches and deserves a place in the portfolio of model-building tools.

Sign in / Sign up

Export Citation Format

Share Document