Daniela Amazonas

Wide-angle FD and FFD migration using complex Padé approximations

Seismic migration by downward continuation using the one-way wave-equation approximations has two shortcomings: imaging steep-dip reflectors and handling evanescent waves. Complex Pade approximations allow a better treatment of evanescent modes, stabilizing finite-difference migration without requiring special treatment for the migration domain boundaries. Imaging of steep-dip reflectors can be improved using several terms in the Pade expansion. We discuss the implementation and evaluation of wide-angle complex Pade approximations for finite-difference and Fourier finite-difference migration methods. The dispersion relation and the impulsive response of the migration operator provide criteria to select the number of terms and coefficients in the Pade expansion. This ensures stability for a prescribed maximum propagation direction. The implementations are validated on the Marmousi model data set and SEG/EAGE salt model data.

Anisotropic complex Padé hybrid finite-difference depth migration

Standard real-valued finite-difference (FD) and Fourier finite-difference (FFD) migrations cannot handle evanescent waves correctly, which can lead to numerical instabilities in the presence of strong velocity variations. A possible solution to these problems is the complex Pade approximation, which avoids problems with evanescent waves by rotating the branch cut of the complex square root. We have applied this approximation to the acoustic wave equation for vertical transversely isotropic media to derive more stable FD and hybrid FD/FFD migrations for such media. Our analysis of the dispersion relation of the new method indicates that it should provide more stable migration results with fewer artifacts and higher accuracy at steep dips. Our studies lead to the conclusion that the rotation angle of the branch cut that should yield the most stable image is 60° for FD migration, as confirmed by numerical impulse responses and work with synthetic data.

Including lateral velocity variations into true-amplitude common-shot wave-equation migration

In heterogeneous media, standard one-way wave equations describe only the kinematic part of one-way wave propagation correctly. For a correct description of amplitudes, the one-way wave equations must be modified. In media with vertical velocity variations only, the resulting true-amplitude one-way wave equations can be solved analytically. In media with lateral velocity variations, these equations are much harder to solve and require sophisticated numerical techniques. We present an approach to circumvent these problems by implementing approximate solutions based on the one-dimensional analytic amplitude modifications. We use these approximations to show how to modify conventional migration methods such as split-step and Fourier finite-difference migrations in such a way that they more accurately handle migration amplitudes. Simple synthetic data examples in media with a constant vertical gradient demonstrate that the correction achieves the recovery of true migration amplitudes. Applications to the SEG/EAGE salt model and the Marmousi data show that the technique improves amplitude recovery in the migrated images in more realistic situations.

Noise Reduction on Numerical Solutions of Partial Differential Equations using Fuzzy Transform

Fuzzy Transform (F-transform) has been introduced as an approximation method which encompasses both classical transforms as well as approximation methods studied in fuzzy modelling and fuzzy control. It has been proved that, under some conditions, Ftransform can remove a periodical noise and it can significantly reduce random noise. In this work we apply the F-transform methodology on the study of numerical solutions of partial differential equations with noisy initial conditions.

Regularization in Slope Tomography

Seismic imaging in depth is limited by the accuracy of velocity model estimation. Slope tomography uses the slowness components and traveltimes of picked reflection or diffraction events for velocity model building. The unavoidable data incompleteness requires additional information to assure stability to inversion. One natural constraint for ray based tomography is a smooth velocity model. We propose a new, reflection-angle-based kind of smoothness constraint as regularization in slope tomography and compare its effects to three other, more conventional constraints. The effect of these constraints are evaluated through angle domain common image gathers, computed with wave-equation migration using the estimated velocity model. We find the smoothness constraints to have a distinct effect on the velocity model but a weaker effect on the migrated data. In numerical tests on synthetic data, the new constraint leads to geologically more consistent models.