Kirchhoff elastic wave migration for the case of noncoincident source and receiver

Geophysics ◽  
1984 ◽  
Vol 49 (8) ◽  
pp. 1223-1238 ◽  
Author(s):  
John T. Kuo ◽  
Ting‐fan Dai
Keyword(s):  
Elastic Wave ◽  
Seismic Data ◽  
Signal To Noise Ratio ◽  
Horizontal Layer ◽  
P Wave ◽  
Seismic Section ◽  
S Waves ◽  
Surface Data ◽  
Marine Seismic ◽  
Wave Migration

In taking into account both compressional (P) and shear (S) waves, more geologic information can likely be extracted from the seismic data. The presence of shear and converted shear waves in both land and marine seismic data recordings calls for the development of elastic wave‐migration methods. The migration method presently developed consists of simultaneous migration of P- and S-waves for offset seismic data based on the Kirchhoff‐Helmholtz type integrals for elastic waves. A new principle of simultaneously migrating both P- and S-waves is introduced. The present method, named the Kirchhoff elastic wave migration, has been tested using the 2-D synthetic surface data calculated from several elastic models of a dipping layer (including a horizontal layer), a composite dipping and horizontal layer, and two layers over a half‐space. The results of these tests not only assure the feasibility of this migration scheme, but also demonstrate that enhanced images in the migrated sections are well formed. Moreover, the signal‐to‐noise ratio increases in the migrated seismic section by this elastic wave migration, as compared with that using the Kirchhoff acoustic (P-) wave migration alone. This migration scheme has about the same order of sensitivity of migration velocity variations, if [Formula: see text] and [Formula: see text] vary concordantly, to the recovery of the reflector as that of the Kirchhoff acoustic (P-) wave migration. In addition, the sensitivity of image quality to the perturbation of [Formula: see text] has also been tested by varying either [Formula: see text] or [Formula: see text]. For varying [Formula: see text] (with [Formula: see text] fixed), the migrated images are virtually unaffected on the [Formula: see text] depth section while they are affected on the [Formula: see text] depth section. For varying [Formula: see text] (with [Formula: see text] fixed), the migrated images are affected on both the [Formula: see text] and [Formula: see text] depth sections.

Real data results of Kirchhoff elastic wave migration

Geophysics ◽  
1986 ◽  
Vol 51 (4) ◽  
pp. 1006-1011 ◽  
Author(s):  
Ting‐Fan Dai ◽  
John T. Kuo
Keyword(s):  
Seismic Data ◽  
Real Data ◽  
Horizontal Layer ◽  
S Wave ◽  
P Waves ◽  
S Waves ◽  
Surface Data ◽  
And Migration ◽  
Wave Migration ◽  
Kirchhoff Integral

Although Kirchhoff integral migration has attracted considerable attention for seismic data processing since the early 1970s, it, like all other seismic migration methods, is only applicable to compressional (P) waves. Because of a recent surge of interest in shear (S) waves, Kuo and Dai (1984) developed the Kirchhoff elastic (P and S) wave migration (KEWM) formulation and migration principle for the case of source and receiver noncoincidence. They obtained encouraging results using two‐dimensional (2-D) synthetic surface data from various geometric elastic models, including a dipping layer, a composite dipping and horizontal layer, and two layers over a half‐space.

scholarly journals PENERAPAN METODE F-K DEMULTIPLE DALAM KASUS ATENUASI WATER-BOTTOM MULTIPLE

2016 ◽  
Vol 11 (1) ◽  
pp. 33
Author(s):  
Subarsyah Subarsyah ◽  
Sahudin Sahudin
Keyword(s):  
Seismic Data ◽  
Final Section ◽  
Signal To Noise Ratio ◽  
Multiple Reflections ◽  
Seismic Section ◽  
Multiple Attenuation ◽  
Seismic Acquisition ◽  
Marine Seismic ◽  
Second Category ◽  
Water Bottom

Keberadaan water-bottom multiple merupakan hal yang tidak bisa dihindari dalam akuisisi data seismik laut, tentu saja hal ini akan menurunkan tingkat perbandingan sinyal dan noise. Beberapa metode atenuasi telah dikembangkan dalam menekan noise ini. Metode atenuasi multiple diklasifikasikan dalam tiga kelompok meliputi metode dekonvolusi yang mengidentifikasi multiple berdasarkan periodisitasnya, metode filtering yang memisahkan refleksi primer dan multiple dalam domain tertentu (F-K,Tau-P dan Radon domain) serta metode prediksi medan gelombang. Penerapan metode F-K demultiple yang masuk kategori kedua akan diterapkan terhadap data seismik PPPGL tahun 2010 di perairan Teluk Tomini. Atenuasi terhadap water-bottom multiple berhasil dilakukan akan tetapi pada beberapa bagian multiple masih terlihat dengan amplitude relatif lebih kecil. F-K demultiple tidak efektif dalam mereduksi multiple pada offset yang pendek dan multiple pada zona ini yang memberikan kontribusi terhadap keberadaan multiple pada penampang akhir. Kata kunci : F-K demultiple, multiple, atenuasi The presence of water-bottom multiple is unavoidable in marine seismic acquisition, of course, this will reduce signal to noise ratio. Several attenuation methods have been developed to suppress this noise. Multiple attenuation methods are classified into three groups first deconvolution method based on periodicity, second filtering method that separates the primary and multiple reflections in certain domains (FK, Tau-P and the Radon domain) ang the third method based on wavefield prediction. Application of F-K demultiple incoming second category will be applied to the seismic data in 2010 PPPGL at Tomini Gulf waters. Attenuation of the water-bottom multiple successful in reduce multiple but in some parts of seismic section multiple still visible with relatively smaller amplitude. FK demultiple not effective in reducing multiple at near offset and multiple in this zone contribute to the existence of multiple in final section. Key words : F-K demultiple, multiple, attenuation

Least-squares data-to-data migration: An approach for migrating free-surface-related multiples

Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S83-S94 ◽  
Author(s):  
Yikang Zheng ◽  
Yibo Wang ◽  
Xu Chang
Keyword(s):  
Free Surface ◽  
Least Squares ◽  
Seismic Data ◽  
Signal To Noise Ratio ◽  
Data Migration ◽  
Subsurface Imaging ◽  
Full Data ◽  
Recording Time ◽  
Numerical Tests ◽  
Marine Seismic

Free-surface-related multiples can provide extra illumination of the subsurface and thus can be usefully included in migration procedures. However, most multiple migration approaches require separation of primaries and free-surface-related multiples or at least prediction of multiples in advance, which is time consuming and prone to errors. The data-to-data migration (DDM) method migrates free-surface-related multiples by forward and backward propagating the recorded full data (containing primaries and free-surface-related multiples). For DDM, there is no need to predict or separate multiples, but the migration results suffer from the crosstalk generated by crosscorrelations of undesired seismic events, e.g., primaries and second-order free-surface-related multiples. We have developed least-squares DDM (LSDDM) for marine data to eliminate the crosstalk generated by DDM. In each iteration, the forward-propagated primaries and free-surface-related multiples are crosscorrelated with the backward-propagated primary and free-surface-related multiple residuals to form the reflectivity gradient. We use a three-layer model and the Marmousi model for numerical tests. The results validate that LSDDM can provide a migrated image with higher signal-to-noise ratio and more balanced amplitudes than DDM. The LSDDM approach might be valuable for general subsurface imaging for marine seismic data when the migration velocity is accurate, and the acquired data have sufficient recording time.

scholarly journals Deep learning for low-frequency extrapolation from multioffset seismic data

Geophysics ◽  
2019 ◽  
Vol 84 (6) ◽  
pp. R989-R1001 ◽  
Author(s):  
Oleg Ovcharenko ◽  
Vladimir Kazei ◽  
Mahesh Kalita ◽  
Daniel Peter ◽  
Tariq Alkhalifah
Keyword(s):  
Deep Learning ◽  
Seismic Data ◽  
Signal To Noise Ratio ◽  
Waveform Inversion ◽  
Low Frequency ◽  
Low Frequencies ◽  
Seismic Surveys ◽  
Full Waveform ◽  
Frequency Components ◽  
Marine Seismic

Low-frequency seismic data are crucial for convergence of full-waveform inversion (FWI) to reliable subsurface properties. However, it is challenging to acquire field data with an appropriate signal-to-noise ratio in the low-frequency part of the spectrum. We have extrapolated low-frequency data from the respective higher frequency components of the seismic wavefield by using deep learning. Through wavenumber analysis, we find that extrapolation per shot gather has broader applicability than per-trace extrapolation. We numerically simulate marine seismic surveys for random subsurface models and train a deep convolutional neural network to derive a mapping between high and low frequencies. The trained network is then tested on sections from the BP and SEAM Phase I benchmark models. Our results indicate that we are able to recover 0.25 Hz data from the 2 to 4.5 Hz frequencies. We also determine that the extrapolated data are accurate enough for FWI application.

Iterative P-wave velocity inversion using Markov chain Monte Carlo method

2020 ◽  
Author(s):  
Hyunggu Jun ◽  
Hyeong-Tae Jou ◽  
Han-Joon Kim ◽  
Sang Hoon Lee
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Velocity Structure ◽  
Low Frequency ◽  
P Wave ◽  
Velocity Analysis ◽  
Seismic Section ◽  
Traveltime Tomography ◽  
P Wave Velocity ◽  
Frequency Components

<p>Imaging the subsurface structure through seismic data needs various information and one of the most important information is the subsurface P-wave velocity. The P-wave velocity structure mainly influences on the location of the reflectors during the subsurface imaging, thus many algorithms has been developed to invert the accurate P-wave velocity such as conventional velocity analysis, traveltime tomography, migration velocity analysis (MVA) and full waveform inversion (FWI). Among those methods, conventional velocity analysis and MVA can be widely applied to the seismic data but generate the velocity with low resolution. On the other hands, the traveltime tomography and FWI can invert relatively accurate velocity structure, but they essentially need long offset seismic data containing sufficiently low frequency components. Recently, the stochastic method such as Markov chain Monte Carlo (McMC) inversion was applied to invert the accurate P-wave velocity with the seismic data without long offset or low frequency components. This method uses global optimization instead of local optimization and poststack seismic data instead of prestack seismic data. Therefore, it can avoid the problem of the local minima and limitation of the offset. However, the accuracy of the poststack seismic section directly affects the McMC inversion result. In this study, we tried to overcome the dependency of the McMC inversion on the poststack seismic section and iterative workflow was applied to the McMC inversion to invert the accurate P-wave velocity from the simple background velocity and inaccurate poststack seismic section. The numerical test showed that the suggested method could successfully invert the subsurface P-wave velocity.</p>

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.

Converted-wave model building and imaging based on common-focus-point methodology

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo
Keyword(s):  
Elastic Wave ◽  
Wave Velocity ◽  
Seismic Data ◽  
Model Building ◽  
P Wave ◽  
Rock Properties ◽  
S Wave ◽  
Converted Wave ◽  
Seismic Images ◽  
Focus Point

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

A plane‐wave decomposition for elastic wave fields applied to the separation of P‐waves and S‐waves in vector seismic data

Geophysics ◽  
1986 ◽  
Vol 51 (2) ◽  
pp. 419-423 ◽  
Author(s):  
A. J. Devaney ◽  
M. L. Oristaglio
Keyword(s):  
Plane Wave ◽  
Elastic Wave ◽  
Seismic Data ◽  
Seismic Profile ◽  
P Waves ◽  
Wave Fields ◽  
S Waves ◽  
Plane Wave Decomposition ◽  
Wave Decomposition ◽  
Seismic Measurements

We describe a method to decompose a two‐dimensional (2-D) elastic wave field recorded along a line into its longitudinal and transverse parts, that is, into compressional (P) waves and shear (S) waves. Separation of the data into P-waves and S-waves is useful when analyzing vector seismic measurements along surface lines or in boreholes. The method described is based on a plane‐wave expansion for elastic wave fields and is illustrated with a synthetic example of an offset vertical seismic profile (VSP) in a layered elastic medium.

PROCESSING AND IMAGING OF MARINE SEISMIC DATA FROM THE JEQUITINHONHA BASIN (BAHIA, BRAZIL)

2015 ◽  
Vol 33 (3) ◽  
pp. 403
Author(s):  
Lourenildo W.B. Leite ◽  
J. Mann ◽  
Wildney W.S. Vieira
Keyword(s):  
Seismic Data ◽  
Oil And Gas ◽  
Signal To Noise Ratio ◽  
Sedimentary Basins ◽  
Data Driven ◽  
Geological Interpretation ◽  
Static Correction ◽  
Time Migration ◽  
Study Results ◽  
Marine Seismic

ABSTRACT. The present study results from a consistent processing and imaging of marine seismic data from a set collected over sedimentary basins of the East Brazilian Atlantic. Our general aim is first to subsidize geological interpretations with plausible subsurface images for oil and gas exploration. In second place, to verify published schematic geological interpretation for these basins by underlying the sediment/basement contact, from where subvertical faults are projected upwards through the basin followed by folded structures. The data-driven results can be used to trace the reflector boundaries in the time sections, submitted to time-to-depth axis transformation, and to be used as a first model for further basin pressure prediction, where natural pumps necessarily develop for the mechanism of oil and gas accumulation. The applied fundamental techniques are mainly based on the data-driven common reflection surface stack, where it is shown the improvement of the signal-to-noise ratio, the lateral continuity of the reflection events, the resolution, and that time migration collapses the diffraction events. The CRS migration strongly collapsed the diffraction events, allowing some subsurface structures be more evident. The free surface and some shallow internal multiples can be clearly traced for further processing aiming at their attenuation. The interpretation lines are meant to show the geometry of selected reflectors, and to help comparing the results with other similar sections. One can trace some subvertical fault systems starting from the lower part of the section (interpreted as the basement), and their extension upwards through the sedimentary sequence.Keywords: CRS stack, CRS migration, residual static correction, NIP wave tomography. RESUMO. O presente artigo resulta de um processamento e imageamento consistentes de dados sísmicos marinhos de levantamento realizado em bacias sedimentares do Atlântico do Nordeste brasileiro. Nossos objetivos gerais são em primeiro lugar subsidiar as interpretações geológicas com imagens plausíveis do subsolo, e voltadas à exploração de óleo e gás. Em segundo lugar, verificar as interpretações geológicas esquemáticas publicadas para estas bacias, para conferir o delineamento do contato sedimento/embasamento, de onde falhas subvertical são projetadas através da bacia, seguidas de estruturas dobradas. Os resultados baseados em dados reais podem ser usados para delinear interfaces refletoras contidas nas seções tempo, submetidos à transformação da coordenada tempo para profundidade, e que podem ser usados posteriormente como um primeiro modelo para a predição de pressão em bacias sedimentares, onde se desenvolve um bombeamento natural necessário para a acumulação de óleo e gás. As técnicas fundamentais aplicadas baseiam-se principalmente no denominado empilhamento de superfície de reflexão comum, baseado em dados observados, onde se mostra a evolução da relação sinal-ruído, da continuidade lateral dos eventos de reflexão, da resolução, e o colapso dos eventos de difração nas seções de migração do tempo. A migração CRS colapsa fortemente os eventos de difração permitindo que algumas estruturas do subsolo sejam mais evidentes. Múltiplas da superfície livre, e algumas internas rasas, podem ser claramente traçadas para processamento adicional que visam a atenuação. As linhas de interpretação trac¸adas visam mostrar a geometria dos refletores selecionados, e ajudar na comparação com outros resultados de seções semelhantes. Pode-se traçar um sistema de falhas subvertical a partir da base inferior (interpretada como o embasamento) da seção escolhida como referência, e os seus prolongamentos através da sequência sedimentar.Palavras-chave: empilhamento CRS, migração CRS, correção estática residual, tomografia NIP.

Sign in / Sign up

Export Citation Format

Share Document