3D Kirchhoff prestack time migration in average illumination-azimuth and incident-angle domain for isotropic and vertical transversely isotropic media

Geophysics ◽  
2011 ◽  
Vol 76 (1) ◽  
pp. S15-S27 ◽  
Author(s):  
Jiubing Cheng ◽  
Jianhua Geng ◽  
Huazhong Wang ◽  
Zaitian Ma
Keyword(s):  
Incident Angle ◽  
Fractured Reservoirs ◽  
Plane Waves ◽  
Transversely Isotropic ◽  
Image Point ◽  
Amplitude Variation ◽  
Time Migration ◽  
Prestack Time Migration ◽  
Ray Bending ◽  
Average Illumination

Conventional offset domain prestack migration tends to bring ambiguity and migration artifacts because it smears energy from different angles at the image point. To avoid this, prestack depth migration implementations in angle domain have been investigated in the past decades. As an efficient imaging tool, angle domain Kirchhoff prestack time migration is still useful and was proposed recently. However, existing algorithms cannot handle ray bending and anisotropy correctly. Practically, azimuth analysis for fractured reservoirs should be carried out after migration for most geological settings. Unfortunately, the existing migration algorithm implicitly involves some kind of binning to source-receiver azimuth, which may not be the real wave-pathazimuth, especially for side-scattering or out-of-plane waves. In this paper, we present an algorithm for 3D Kirchhoff prestack time migration in average illumination azimuth and incident angle domain, which matches true wave path naturally and more accurately. To handle ray bending and vertical transversely isotropy, we propose several approaches to estimate two-way traveltime and the corresponding angular attributes through extended offset-to-angle mapping. Based upon these approaches, our 3D prestack time migration can provide high-quality common-image gathers for amplitude variation with incident angle and/or amplitude variation with offset and azimuth analyses, even in media with slight to moderate lateral heterogeneity. The 2D and 3D synthetic examples prove the validity of our methods.

Improving lateral and vertical resolution of seismic images by correcting for wavelet stretch in common-angle migration

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. C95-C104 ◽  
Author(s):  
Gabriel Perez ◽  
Kurt J. Marfurt
Keyword(s):  
Signal To Noise Ratio ◽  
Incident Angle ◽  
Large Angle ◽  
Acoustic Medium ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Velocity Anisotropy ◽  
Time Migration ◽  
Reflection Angle ◽  
Seismic Reflections

Long-offset or high-incident-angle seismic reflections provide us with improved velocity resolution, better leverage against multiples, less contamination by ground roll, and information that is often critical when estimating lithology and fluid product. Unfortunately, high-incident-angle seismic reflections suffer not only from nonhyperbolic moveout but also from wavelet stretch during imaging, resulting in lower-resolution images that mix the response from adjacent lithologies. For an arbitrary acoustic medium, wavelet stretch from prestack migration depends only on the cosine of the reflection angle, such that the amount of wavelet stretch will be the same for all samples of a common-reflection-angle migrated trace. Thus, we are able to implement a wavelet stretch correction by applying a simple stationary spectral shaping operation to common-angle migrated traces. We obtain such traces directly by a prestack Kirchhoff migration algorithm. Correcting for stretch effectively increases the fold of imaged data, far beyond that achieved in conventional migration, resulting in improved signal-to-noise ratio of the final stacked section. Increasing the fidelity of large incident angles results in images with improved vertical and lateral resolution and with increased angular illumination, valuable for amplitude variation with angle (AVA) and amplitude variation with offset (AVO) analysis. Finally, such large-angle images are more sensitive to and therefore provide increased leverage over errors in velocity and velocity anisotropy. These ideas were applied to prestack time migration on seismic data from the Fort Worth basin, in Texas.

Enhancements to prestack frequency‐wavenumber (f-k) migration

Geophysics ◽  
1991 ◽  
Vol 56 (1) ◽  
pp. 27-40 ◽  
Author(s):  
Z. Li ◽  
W. Lynn ◽  
R. Chambers ◽  
Ken Larner ◽  
Ray Abma
Keyword(s):  
Computational Effort ◽  
Velocity Variation ◽  
Lateral Velocity ◽  
Time Migration ◽  
Migration Velocity Analysis ◽  
Prestack Time Migration ◽  
And Migration ◽  
Definition Of ◽  
Ray Bending ◽  
Excellent Source

Prestack frequency‐wavenumber (f-k) migration is a particularly efficient method of doing both full prestack time migration and migration velocity analysis. Conventional implementations of the method, however, can encounter several drawbacks: (1) poor resolution and spatial aliasing noise caused by insufficient sampling in the offset dimension, (2) poor definition of steep events caused by insufficient sampling in the velocity dimension, and (3) inadequate handling of ray bending for steep events. All three of these problems can be mitigated with modifications to the prestack f-k algorithm. The application of linear moveout (LMO) in the offset dimension prior to migration reduces event moveout and hence increases the bandwidth of non‐spatially aliased signals. To reduce problems of interpolation for steep events, the number of constant‐velocity migrations can be economically increased by performing residual poststack migrations. Finally, migration with a dip‐dependent imaging velocity addresses the issue of ray bending and thereby improves the positioning of steep events. None of these enhancements substantially increases the computational effort of f-k migration. Prestack f-k migration possesses a limitation for which no solution is readily available. Where lateral velocity variation is modest, steep events (such as fault‐plane reflections in sediments) may not be imaged as well as by other migration approaches. This shortcoming results from the restriction that, in the prestack f-k approach, a single velocity field must serve to perform two different functions: imaging and stacking. Nevertheless, in areas of strong velocity variation and gentle to moderate dip, the detailed velocity control afforded by the prestack f-k method is an excellent source of geologic information.

Common-image gathers in the incident phase-angle domain from reverse time migration in 2D elastic VTI media

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. S197-S206 ◽  
Author(s):  
Qunshan Zhang ◽  
George A. McMechan
Keyword(s):  
Phase Angle ◽  
Incident Angle ◽  
Image Point ◽  
Common Source ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wave Phase ◽  
Time Migration ◽  
Imaging Time ◽  
Vti Media

Reverse time migration (RTM) was implemented with a modified crosscorrelation imaging condition for data from 2D elastic vertically transversely isotropy (VTI) media. The computation cost was reduced because scalar qP- and qS-wavefield separations are performed in VTI media, for the source and receiver wavefields only at the RTM imaging time, to calculate the migrated qP and qS images. Angle-domain common-image gathers (CIGs) were extracted from qPqP and qPqS common-source RTM images. The local incident angle was produced as the difference between the qP-wave phase angle, obtained directly from the source wavefield polarization, and the normal to the reflector, calculated as the instantaneous wavenumber direction via a directional Hilbert transform of the stacked image. Angle-domain CIGs were extracted by reordering the prestack-migrated images by local incident phase angle, source by source. Vector decomposition of the source qP-wavefield was required to calculate the qP-wave phase polarization direction for each image point at its imaging time. RTM and angle-domain CIG extraction were successfully implemented and illustrated with a synthetic 2D elastic VTI example.

scholarly journals Reection and refraction of plane waves at the interface between magnetoelectroelastic and liquid media

2013 ◽  
Vol 40 (3) ◽  
pp. 427-439 ◽  
Author(s):  
Rong Zhang
Keyword(s):  
Liquid Medium ◽  
Incident Angle ◽  
Plane Waves ◽  
Transversely Isotropic ◽  
Liquid Media ◽  
Reflection And Transmission ◽  
Transmission Coefficients ◽  
Reflection And Refraction ◽  
Acoustic Device ◽  
Wave Incidence

This paper analyzes the reflection and refraction of plane wave incidences at the interface between magnetoelectroelastic (MEE) and liquid media. The MEE medium is assumed to be transversely isotropic and the liquid medium to be nonviscous. Three cases, i.e., the coupled quasipressure wave incidence from the MEE medium, the coupled quasi-shear vertical wave incidence from the MEE medium, and the pressure wave incidence from the liquid medium, are discussed. The expressions of reflection and transmission coefficients varying with the incident angle are obtained. This investigation would be useful to the MEE acoustic device field.

Azimuth-preserved local angle-domain prestack time migration in isotropic, vertical transversely isotropic and azimuthally anisotropic media

Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. S51-S64 ◽  
Author(s):  
Jiubing Cheng ◽  
Tengfei Wang ◽  
Chenlong Wang ◽  
Jianhua Geng
Keyword(s):  
Seismic Imaging ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Azimuthal Anisotropy ◽  
Velocity Estimation ◽  
Energy Propagation ◽  
Directional Information ◽  
Time Migration ◽  
Prestack Time Migration ◽  
Local Angle

Conventional prestack migration does not preserve local directional information of the seismic waves at the image points. New attempts such as sectored migration of azimuth-limited or common-offset-vector data only concern source-receiver azimuth and offset on the surface, which can be poor representation of subsurface wavepath direction. Moreover, they could result in inaccurate imaging because they do not account for the energy propagation between azimuths or offset-vectors. In the past decade, local angle-domain seismic imaging has been highly advocated to avoid migration artifacts and to improve velocity estimation in complex media. Considering prestack time migration (PSTM) is still widely used in seismic imaging and seismic data preconditioning for amplitude variations with offset or incident-angle (AVO/AVA) analysis, fracture detection, and reservoir characterization, we present an azimuth-preserved local angle-domain Kirchhoff PSTM approach for such purposes. We apply a seismic imaging condition in 3D local angle domain and use extended superposition of impulse responses retaining subsurface angular attributes, which are evaluated through the incident and scattering phase slowness vectors using classical-diffraction moveout equations in isotropic, vertical transversely isotropic (VTI) and azimuthally anisotropic media. Two-dimensional synthetic examples demonstrate what the migrated results look like in local angle domain. A wide-azimuth synthetic example with horizontal transversely isotropy (HTI) proves the necessity of azimuthal migration for reliable imaging and azimuthal analysis when azimuthal anisotropy exists in the overburden. Real data examples show the advantages of imaging in subsurface angle domain for properly focusing and revealing azimuth- and angle-dependent variations of residual moveout and migrated amplitudes.

Propagation and Reflection of Plane Waves in a Rotating Magneto-Elastic Fibre-Reinforced Semi Space with Surface Stress

2020 ◽  
Vol 22 (4) ◽  
pp. 939-958
Author(s):  
Indrajit Roy ◽  
D. P. Acharya ◽  
Sourav Acharya
Keyword(s):  
Surface Stress ◽  
Incident Angle ◽  
Plane Waves ◽  
Three Dimensional ◽  
Elastic Fibre ◽  
Amplitude Ratios ◽  
Governing Equations ◽  
Rotation Surface ◽  
Wave Velocities ◽  
Magnetic Field Rotation

AbstractThe present paper investigates the propagation of quasi longitudinal (qLD) and quasi transverse (qTD) waves in a magneto elastic fibre-reinforced rotating semi-infinite medium. Reflections of waves from the flat boundary with surface stress have been studied in details. The governing equations have been used to obtain the polynomial characteristic equation from which qLD and qTD wave velocities are found. It is observed that both the wave velocities depend upon the incident angle. After imposing the appropriate boundary conditions including surface stress the resultant amplitude ratios for the total displacements have been obtained. Numerically simulated results have been depicted graphically by displaying two and three dimensional graphs to highlight the influence of magnetic field, rotation, surface stress and fibre-reinforcing nature of the material medium on the propagation and reflection of plane waves.

Reverse time migration in tilted transversely isotropic media

2020 ◽  
Vol 38 (2) ◽  
Author(s):  
Razec Cezar Sampaio Pinto da Silva Torres ◽  
Leandro Di Bartolo
Keyword(s):  
Seismic Anisotropy ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Computer Clusters ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Tti Media

ABSTRACT. Reverse time migration (RTM) is one of the most powerful methods used to generate images of the subsurface. The RTM was proposed in the early 1980s, but only recently it has been routinely used in exploratory projects involving complex geology – Brazilian pre-salt, for example. Because the method uses the two-way wave equation, RTM is able to correctly image any kind of geological environment (simple or complex), including those with anisotropy. On the other hand, RTM is computationally expensive and requires the use of computer clusters. This paper proposes to investigate the influence of anisotropy on seismic imaging through the application of RTM for tilted transversely isotropic (TTI) media in pre-stack synthetic data. This work presents in detail how to implement RTM for TTI media, addressing the main issues and specific details, e.g., the computational resources required. A couple of simple models results are presented, including the application to a BP TTI 2007 benchmark model.Keywords: finite differences, wave numerical modeling, seismic anisotropy. Migração reversa no tempo em meios transversalmente isotrópicos inclinadosRESUMO. A migração reversa no tempo (RTM) é um dos mais poderosos métodos utilizados para gerar imagens da subsuperfície. A RTM foi proposta no início da década de 80, mas apenas recentemente tem sido rotineiramente utilizada em projetos exploratórios envolvendo geologia complexa, em especial no pré-sal brasileiro. Por ser um método que utiliza a equação completa da onda, qualquer configuração do meio geológico pode ser corretamente tratada, em especial na presença de anisotropia. Por outro lado, a RTM é dispendiosa computacionalmente e requer o uso de clusters de computadores por parte da indústria. Este artigo apresenta em detalhes uma implementação da RTM para meios transversalmente isotrópicos inclinados (TTI), abordando as principais dificuldades na sua implementação, além dos recursos computacionais exigidos. O algoritmo desenvolvido é aplicado a casos simples e a um benchmark padrão, conhecido como BP TTI 2007.Palavras-chave: diferenças finitas, modelagem numérica de ondas, anisotropia sísmica.

High-resolution imaging: An approach by incorporating stationary-phase implementation into deabsorption prestack time migration

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. S317-S331 ◽  
Author(s):  
Jianfeng Zhang ◽  
Zhengwei Li ◽  
Linong Liu ◽  
Jin Wang ◽  
Jincheng Xu
Keyword(s):  
Stationary Phase ◽  
Fresnel Zone ◽  
Computational Cost ◽  
Computational Effort ◽  
Propagation Path ◽  
Time Migration ◽  
Resolution Imaging ◽  
Prestack Time Migration ◽  
Dip Angle ◽  
Absorption And Dispersion

We have improved the so-called deabsorption prestack time migration (PSTM) by introducing a dip-angle domain stationary-phase implementation. Deabsorption PSTM compensates absorption and dispersion via an actual wave propagation path using effective [Formula: see text] parameters that are obtained during migration. However, noises induced by the compensation degrade the resolution gained and deabsorption PSTM requires more computational effort than conventional PSTM. Our stationary-phase implementation improves deabsorption PSTM through the determination of an optimal migration aperture based on an estimate of the Fresnel zone. This significantly attenuates the noises and reduces the computational cost of 3D deabsorption PSTM. We have estimated the 2D Fresnel zone in terms of two dip angles through building a pair of 1D migrated dip-angle gathers using PSTM. Our stationary-phase QPSTM (deabsorption PSTM) was implemented as a two-stage process. First, we used conventional PSTM to obtain the Fresnel zones. Then, we performed deabsorption PSTM with the Fresnel-zone-based optimized migration aperture. We applied stationary-phase QPSTM to a 3D field data. Comparison with synthetic seismogram generated from well log data validates the resolution enhancements.

3D constraints and finite-difference modeling of massive sulfide deposits: The Kristineberg seismic lines revisited, northern Sweden

Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. WC69-WC79 ◽  
Author(s):  
Mahdieh Dehghannejad ◽  
Alireza Malehmir ◽  
Christopher Juhlin ◽  
Pietari Skyttä
Keyword(s):  
Finite Difference ◽  
Seismic Data ◽  
Massive Sulfide ◽  
Northern Sweden ◽  
Seismic Reflection Data ◽  
Sulfide Deposits ◽  
Time Migration ◽  
Massive Sulfide Deposits ◽  
Prestack Time Migration ◽  
Seismic Lines

The Kristineberg mining area in the western part of the Skellefte ore district is the largest base metal producer in northern Sweden and currently the subject of extensive geophysical and geologic studies aimed at constructing 3D geologic models. Seismic reflection data form the backbone of the geologic modeling in the study area. A geologic cross section close to the Kristineberg mine was used to generate synthetic seismic data using acoustic and elastic finite-difference algorithms to provide further insight about the nature of reflections and processing challenges when attempting to image the steeply dipping structures within the study area. Synthetic data suggest processing artifacts manifested themselves in the final 2D images as steeply dipping events that could be confused with reflections. Fewer artifacts are observed when the data are processed using prestack time migration. Prestack time migration also was performed on high-resolution seismic data recently collected near the Kristineberg mine and helped to image a high-amplitude, gently dipping reflection occurring stratigraphically above the extension of the deepest Kristineberg deposit. Swath 3D processing was applied to two crossing seismic lines, west of the Kristineberg mine, to provide information on the 3D geometry of an apparently flat-lying reflection observed in both of the profiles. The processing indicated that the reflection dips about 30° to the southwest and is generated at the contact between metasedimentary and metavolcanic rocks, the upper part of the latter unit being the most typical stratigraphic level for the massive sulfide deposits in the Skellefte district.

Sign in / Sign up

Export Citation Format

Share Document