Seismic imaging using lateral adaptive windows

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S73-S79
Author(s):  
Ørjan Pedersen ◽  
Sverre Brandsberg-Dahl ◽  
Bjørn Ursin

One-way wavefield extrapolation methods are used routinely in 3D depth migration algorithms for seismic data. Due to their efficient computer implementations, such one-way methods have become increasingly popular and a wide variety of methods have been introduced. In salt provinces, the migration algorithms must be able to handle large velocity contrasts because the velocities in salt are generally much higher than in the surrounding sediments. This can be a challenge for one-way wavefield extrapolation methods. We present a depth migration method using one-way propagators within lateral windows for handling the large velocity contrasts associated with salt-sediment interfaces. Using adaptive windowing, we can handle large perturbations locally in a similar manner as the beamlet propagator, thus limiting the impact of the errors on the global wavefield. We demonstrate the performance of our method by applying it to synthetic data from the 2D SEG/EAGE [Formula: see text] salt model and an offshore real data example.

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.


2003 ◽  
Vol 46 (6) ◽  
pp. 1176-1185 ◽  
Author(s):  
Shengchang CHEN ◽  
Jingzhong CAO ◽  
Zaitian MA

2003 ◽  
Vol 2003 (2) ◽  
pp. 1-4
Author(s):  
James Sun ◽  
Carl Notfors ◽  
Zhang Yu ◽  
Gray Sam ◽  
Young Jerry

Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. S17-S26 ◽  
Author(s):  
Daniel A. Rosales ◽  
Sergey Fomel ◽  
Biondo L. Biondi ◽  
Paul C. Sava

Wavefield-extrapolation methods can produce angle-domain common-image gathers (ADCIGs). To obtain ADCIGs for converted-wave seismic data, information about the image dip and the P-to-S velocity ratio must be included in the computation of angle gathers. These ADCIGs are a function of the half-aperture angle, i.e., the average between the incidence angle and the reflection angle. We have developed a method that exploits the robustness of computing 2D isotropic single-mode ADCIGs and incorporates both the converted-wave velocity ratio [Formula: see text] and the local image dip field. It also maps the final converted-wave ADCIGs into two ADCIGs, one a function of the P-incidence angle and the other a function of the S-reflection angle. Results with both synthetic and real data show the practical application for converted-wave ADCIGs. The proposed approach is valid in any situation as long as the migration algorithm is based on wavefield downward continuation and the final prestack image is a function of the horizontal subsurface offset.


Geophysics ◽  
1999 ◽  
Vol 64 (5) ◽  
pp. 1535-1545 ◽  
Author(s):  
Lian‐Jie Huang ◽  
Michael C. Fehler ◽  
Peter M. Roberts ◽  
Charles C. Burch

We develop a novel depth‐migration method termed the extended local Rytov Fourier (ELRF) migration method. It is based on the scalar wave equation and a local application of the Rytov approximation within each extrapolation interval. Wavefields are Fourier transformed back and forth between the frequency‐space and frequency‐wavenumber domains during wavefield extrapolation. The lateral slowness variations are taken into account in the frequency‐space domain. The method is efficient due to the use of a fast Fourier transform algorithm. Under the small angle approximation, the ELRF method leads to the split‐step Fourier (SSF) method that is unconditionally stable. The ELRF method and the extended local Born Fourier (ELBF) method that we previously developed can handle wider propagation angles than the SSF method and account for the phase and amplitude changes due to the lateral variations of slowness, whereas the SSF method only accounts for the phase changes. The stability of the ELRF method is controlled more easily than that of the ELBF method.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB203-WB210 ◽  
Author(s):  
Gilles Hennenfent ◽  
Lloyd Fenelon ◽  
Felix J. Herrmann

We extend our earlier work on the nonequispaced fast discrete curvelet transform (NFDCT) and introduce a second generation of the transform. This new generation differs from the previous one by the approach taken to compute accurate curvelet coefficients from irregularly sampled data. The first generation relies on accurate Fourier coefficients obtained by an [Formula: see text]-regularized inversion of the nonequispaced fast Fourier transform (FFT) whereas the second is based on a direct [Formula: see text]-regularized inversion of the operator that links curvelet coefficients to irregular data. Also, by construction the second generation NFDCT is lossless unlike the first generation NFDCT. This property is particularly attractive for processing irregularly sampled seismic data in the curvelet domain and bringing them back to their irregular record-ing locations with high fidelity. Secondly, we combine the second generation NFDCT with the standard fast discrete curvelet transform (FDCT) to form a new curvelet-based method, coined nonequispaced curvelet reconstruction with sparsity-promoting inversion (NCRSI) for the regularization and interpolation of irregularly sampled data. We demonstrate that for a pure regularization problem the reconstruction is very accurate. The signal-to-reconstruction error ratio in our example is above [Formula: see text]. We also conduct combined interpolation and regularization experiments. The reconstructions for synthetic data are accurate, particularly when the recording locations are optimally jittered. The reconstruction in our real data example shows amplitudes along the main wavefronts smoothly varying with limited acquisition imprint.


Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 251-260 ◽  
Author(s):  
Gary F. Margrave

The signal band of reflection seismic data is that portion of the temporal Fourier spectrum which is dominated by reflected source energy. The signal bandwidth directly determines the spatial and temporal resolving power and is a useful measure of the value of such data. The realized signal band, which is the signal band of seismic data as optimized in processing, may be estimated by the interpretation of appropriately constructed f-x spectra. A temporal window, whose length has a specified random fluctuation from trace to trace, is applied to an ensemble of seismic traces, and the temporal Fourier transform is computed. The resultant f-x spectra are then separated into amplitude and phase sections, viewed as conventional seismic displays, and interpreted. The signal is manifested through the lateral continuity of spectral events; noise causes lateral incoherence. The fundamental assumption is that signal is correlated from trace to trace while noise is not. A variety of synthetic data examples illustrate that reasonable results are obtained even when the signal decays with time (i.e., is nonstationary) or geologic structure is extreme. Analysis of real data from a 3-C survey shows an easily discernible signal band for both P-P and P-S reflections, with the former being roughly twice the latter. The potential signal band, which may be regarded as the maximum possible signal band, is independent of processing techniques. An estimator for this limiting case is the corner frequency (the frequency at which a decaying signal drops below background noise levels) as measured on ensemble‐averaged amplitude spectra from raw seismic data. A comparison of potential signal band with realized signal band for the 3-C data shows good agreement for P-P data, which suggests the processing is nearly optimal. For P-S data, the realized signal band is about half of the estimated potential. This may indicate a relative immaturity of P-S processing algorithms or it may be due to P-P energy on the raw radial component records.


Geosciences ◽  
2018 ◽  
Vol 8 (12) ◽  
pp. 497
Author(s):  
Fedor Krasnov ◽  
Alexander Butorin

Sparse spikes deconvolution is one of the oldest inverse problems, which is a stylized version of recovery in seismic imaging. The goal of sparse spike deconvolution is to recover an approximation of a given noisy measurement T = W ∗ r + W 0 . Since the convolution destroys many low and high frequencies, this requires some prior information to regularize the inverse problem. In this paper, the authors continue to study the problem of searching for positions and amplitudes of the reflection coefficients of the medium (SP&ARCM). In previous research, the authors proposed a practical algorithm for solving the inverse problem of obtaining geological information from the seismic trace, which was named A 0 . In the current paper, the authors improved the method of the A 0 algorithm and applied it to the real (non-synthetic) data. Firstly, the authors considered the matrix approach and Differential Evolution approach to the SP&ARCM problem and showed that their efficiency is limited in the case. Secondly, the authors showed that the course to improve the A 0 lays in the direction of optimization with sequential regularization. The authors presented calculations for the accuracy of the A 0 for that case and experimental results of the convergence. The authors also considered different initialization parameters of the optimization process from the point of the acceleration of the convergence. Finally, the authors carried out successful approbation of the algorithm A 0 on synthetic and real data. Further practical development of the algorithm A 0 will be aimed at increasing the robustness of its operation, as well as in application in more complex models of real seismic data. The practical value of the research is to increase the resolving power of the wave field by reducing the contribution of interference, which gives new information for seismic-geological modeling.


2017 ◽  
Vol 5 (1) ◽  
pp. T1-T9 ◽  
Author(s):  
Rui Zhang ◽  
Kui Zhang ◽  
Jude E. Alekhue

More and more seismic surveys produce 3D seismic images in the depth domain by using prestack depth migration methods, which can present a direct subsurface structure in the depth domain rather than in the time domain. This leads to the increasing need for applications of seismic inversion on the depth-imaged seismic data for reservoir characterization. To address this issue, we have developed a depth-domain seismic inversion method by using the compressed sensing technique with output of reflectivity and band-limited impedance without conversion to the time domain. The formulations of the seismic inversion in the depth domain are similar to time-domain methods, but they implement all the elements in depth domain, for example, a depth-domain seismic well tie. The developed method was first tested on synthetic data, showing great improvement of the resolution on inverted reflectivity. We later applied the method on a depth-migrated field data with well-log data validated, showing a great fit between them and also improved resolution on the inversion results, which demonstrates the feasibility and reliability of the proposed method on depth-domain seismic data.


Sign in / Sign up

Export Citation Format

Share Document