24. Caustics in qSV Rayfields of Transversely Isotropic and Vertically Inhomogeneous Media

2001 ◽  
pp. 409-418 ◽  
Author(s):  
Andrzej Hanyga ◽  
M. A. Slawinski
Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1139-1150 ◽  
Author(s):  
Tariq Alkhalifah

I develop an efficient modeling technique for transversely isotropic, inhomogeneous media using a mix of analytical equations and numerical calculations. The analytic equation for the raypath in a factorized transversely isotropic (FTI) media with linear velocity variation, derived by Shearer and Chapman, is used to trace rays between two points. In addition, I derive an analytical equation for geometrical spreading in FTI media that aids in preserving program efficiency; however, the traveltimes are calculated numerically. I then generalize the method to treat general transversely isotropic (TI) media that are not factorized anisotropic inhomogeneous by perturbing the FTI traveltimes, following the perturbation ideas of Červený and Filho. A Kirchhoff‐summation‐based program relying on Trorey’s diffraction method is used to generate synthetic seismograms for such a medium. For the type of velocity models treated, the program is much more efficient than finite‐difference or general ray‐trace modeling techniques.


Geophysics ◽  
1993 ◽  
Vol 58 (10) ◽  
pp. 1454-1467 ◽  
Author(s):  
Ken L. Larner ◽  
Jack K. Cohen

Given the sensitivity of imaging accuracy to the velocity used in migration, migration founded (as in practice) on the erroneous assumption that a medium is isotropic can be expected to be inaccurate for steep reflectors. Here, we estimate errors in interpreted reflection time and lateral position as a function of reflector dip for transversely isotropic models in which the axis of symmetry is vertical and the medium velocity varies linearly with depth. We limit consideration to media in which ratios of the various elastic moduli are independent of depth. Tests with reflector dips up to 120 degrees on a variety of anisotropic media show errors that are tens of wavelengths for dips beyond 90 degrees when the medium (unrealistically) is homogeneous. For a given anisotropy, the errors are smaller for inhomogeneous media; the larger the velocity gradient, the smaller the errors. For gradients that are representative of the subsurface, lateral‐position errors tend to be minor for dips less than about 60 degrees, growing to two to five wavelengths as dip passes beyond 90 degrees. These errors depend on reflector depth and average velocity to the reflector only through their ratio, i.e., migrated reflection time. Migration error, which is found to be unrelated to the ratio of horizontal to vertical velocity, is such that reflections with later migrated reflection times tend to be more severely overmigrated than are those with earlier times. Over a large range of dips, migration errors that arise when anisotropy is ignored but inhomogeneity is honored tend to be considerably smaller than those encountered when inhomogeneity is ignored in migrating data from isotropic, inhomogeneous media.


Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1079-1092 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin

We present a new equation for normal‐moveout (NMO) velocity that describes azimuthally dependent reflection traveltimes of pure modes from both horizontal and dipping reflectors in arbitrary anisotropic inhomogeneous media. With the exception of anomalous areas such as those where common‐midpoint (CMP) reflection time decreases with offset, the azimuthal variation of NMO velocity represents an ellipse in the horizontal plane, with the orientation of the axes determined by the properties of the medium and the direction of the reflector normal. In general, a minimum of three azimuthal measurements is necessary to reconstruct the best‐fit ellipse and obtain NMO velocity in all azimuthal directions. This result provides a simple way to correct for the azimuthal variation in stacking velocity often observed in 3-D surveys. Even more importantly, analytic expressions for the parameters of the NMO ellipse can be used in the inversion of moveout data for the anisotropic coefficients of the medium. For homogeneous transversely isotropic media with a vertical axis of symmetry (VTI media), our equation for azimuthally dependent NMO velocity from dipping reflectors becomes a relatively simple function of phase velocity and its derivatives. We show that the zero‐dip NMO velocity Vnmo(0) and the anisotropic coefficient η are sufficient to describe the P-wave NMO velocity for any orientation of the CMP line with respect to the dip plane of the reflector. Using our formalism, Vnmo(0) and η (the only parameters needed for time processing) can be found from the dip‐dependent NMO velocity at any azimuth or, alternatively, from the azimuthally dependent NMO for a single dipping reflector. We also apply this theory to more complicated azimuthally anisotropic models with the orthorhombic symmetry used to describe fractured reservoirs. For reflections from horizontal interfaces in orthorhombic media, the axes of the normal moveout ellipse are aligned with the vertical symmetry planes. Therefore, azimuthal P-wave moveout measurements can be inverted for the orientation of the symmetry planes (typically determined by the fracture direction) and the NMO velocities within them. If the vertical velocity is known, symmetry‐plane NMO velocities make it possible to estimate two anisotropic parameters equivalent to Thomsen’s coefficient δ for transversely isotropic media.


Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 662-675 ◽  
Author(s):  
Tariq Alkhalifah

The first and most important step in processing data in transversely isotropic (TI) media for which velocities vary with depth is parameter estimation. The multilayer normal‐moveout (NMO) equation for a dipping reflector provides the basis for extending the TI velocity analysis of Alkhalifah and Tsvankin to vertically inhomogeneous media. This NMO equation is based on a root‐mean‐square (rms) average of interval NMO velocities that correspond to a single ray parameter, that of the dipping event. Therefore, interval NMO velocities [including the normal‐moveout velocity for horizontal events, [Formula: see text]] can be extracted from the stacking velocities using a Dix‐type differentiation procedure. On the other hand, η, which is a key combination of Thomsen's parameters that time‐related processing relies on, is extracted from the interval NMO velocities using a homogeneous inversion within each layer. Time migration, like dip moveout, depends on the same two parameters in vertically inhomogeneous media, namely [Formula: see text] and η, both of which can vary with depth. Therefore, [Formula: see text] and ε estimated using the dip dependency of P‐wave moveout velocity can be used for TI time migration. An application of anisotropic processing to seismic data from offshore Africa demonstrates the importance of considering anisotropy, especially as it pertains to focusing and imaging of dipping events.


Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1474-1484 ◽  
Author(s):  
Tariq Alkhalifah

Gaussian beam migration (GBM), as it is implemented today, efficiently handles isotropic inhomogeneous media. The approach is based on the solution of the wave equation in ray‐centered coordinates. Here, I extend the method to work for 2-D migration in generally anisotropic inhomogeneous media. Extension of the Gaussian‐beam method from isotropic to anisotropic media involves modification of the kinematics and dynamics in the required ray tracing. While the accuracy of the paraxial expansion for anisotropic media is comparable to that for isotropic media, ray tracing in anisotropic media is much slower than that in isotropic media. However, because ray tracing is just a small portion of the computation in GBM, the increased computational effort in general anisotropic GBM is typically only about 40%. Application of this method to synthetic examples shows successful migration in inhomogeneous, transversely isotropic media for reflector dips up to and beyond 90°. Further applications to synthetic data of layered anisotropic media show the importance of applying the proper smoothing to the velocity field used in the migration. Also, tests with synthetic data show that the quality of anisotropic migration of steep events in a medium with velocity increasing with depth is much more sensitive to the Thomsen anisotropy parameter ε than to the parameter δ. Thus, a good estimate of ε is needed to apply anisotropic migration with confidence.


Sign in / Sign up

Export Citation Format

Share Document