Gladys Gonzalez

Multipathing in three-parameter common-image gathers from reverse-time migration

ABSTRACT Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region.

Comparison of Irregular Cartesian Finite Difference Methods For Acoustic RTM

We compare three Irregular Cartesian Finite Difference (ICFD) methods towards deployment into Forward Modeling and RTM. The ICFD methods considered are based on: Fornberg (1988), Gamet et al. (1999) and Pitarka et al. (1999). These methods consist of providing a two-way acoustic wave propagation adapted to the P-wave velocity field. While respecting stability and low numerical dispersion criteria, they avoid the spatial oversampling introduced by a uniform spacing. Thus they reduce in a significant manner the required computer memory and execution time. We show in examples that ICFD methods allow to obtain results as accurate as high-order Regular Cartesian Finite Difference methods. Also, we evaluate the accuracy and the performance of the non-uniform stencil computation for each of the three ICFD methods. Finally, we validate that ICFD methods provide significant speed-up for FD propagators-based applications.