Full-waveform inversion FWI is a challenging data-fitting procedure based on full-wavefield modeling to extract quantitative information from seismograms. High-resolution imaging at half the propagated wavelength is expected. Recent advances in high-performance computing and multifold/multicomponent wide-aperture and wide-azimuth acquisitions make 3D acoustic FWI feasible today. Key ingredients of FWI are an efficient forward-modeling engine and a local differential approach, in which the gradient and the Hessian operators are efficiently estimated. Local optimization does not, however, prevent convergence of the misfit function toward local minima because of the limited accuracy of the starting model, the lack of low frequencies, the presence of noise, and the approximate modeling of the wave-physics complexity. Different hierarchical multiscale strategiesaredesignedtomitigatethenonlinearityandill-posedness of FWI by incorporating progressively shorter wavelengths in the parameter space. Synthetic and real-data case studies address reconstructing various parameters, from VP and VS velocities to density, anisotropy, and attenuation. This review attempts to illuminate the state of the art of FWI. Crucial jumps, however, remain necessary to make it as popular as migration techniques. The challenges can be categorized as 1 building accurate starting models with automatic procedures and/or recording low frequencies, 2 defining new minimization criteria to mitigate the sensitivity of FWI to amplitude errors and increasing the robustness of FWI when multiple parameter classes are estimated, and 3 improving computational efficiency by data-compression techniquestomake3DelasticFWIfeasible.

/pdf/an-overview-of-full-waveform-inversion-in-exploration-59sei7fzvh.pdf

An overview of full-waveform inversion in exploration geophysics

Migration of stacked or zero-offset sections is based on deriving the wave amplitude in space from wave field observations at the surface. Conventionally this calculation has been carried out through a depth extrapolation. We examine the alternative of carrying out the migration through a reverse time extrapolation. This approach may offer improvements over existing migration methods, especially in cases of steeply dipping structures with strong velocity contrasts. This migration method is tested using appropriate synthetic data sets.

/pdf/reverse-time-migration-46dbsrafae.pdf

Reverse time migration

Ground-penetrating radar and its use in sedimentology: principles, problems and progress

A space-variant interpolation is required to compensate for the migration of signal energy through range resolution cells when processing synthetic aperture radar (SAR) data, using either the classical range/Doppler (R/D) algorithm or related frequency domain techniques. In general, interpolation requires significant computation time, and leads to loss of image quality, especially in the complex image. The new chirp scaling algorithm avoids interpolation, yet performs range cell migration correction accurately. The algorithm requires only complex multiplies and Fourier transforms to implement, is inherently phase preserving, and is suitable for wide-swath, large-beamwidth, and large-squint applications. This paper describes the chirp scaling algorithm, summarizes simulation results, presents imagery processed with the algorithm, and reviews quantitative measures of its performance. Based on quantitative comparison, the chirp scaling algorithm provides image quality equal to or better than the precision range/Doppler processor. Over the range of parameters tested, image quality results approach the theoretical limit, as defined by the system bandwidth. >

Precision SAR processing using chirp scaling

Accurate methods for the solution of the migration of zero-offset seismic records have been developed. The numerical operations are defined in the frequency domain. The source and recorder positions are lowered by means of a phase shift, or a rotation of the phase angle of the Fourier coefficients. For applications with laterally invariant velocities, the equations governing the migration process are solved very accurately by the phase-shift method. The partial differential equations considered include the 15 degree equation, as well as higher order approximations to the exact migration process. The most accurate migration is accomplished by using the asymptotic equation, whose dispersion relation is the same as that of the full wave equation for downward propagating waves. These equations, however, do not account for the reflection and transmission effects, multiples, or evanescent waves. For comparable accuracy, the present approach to migration is expected to be computationally more efficient than finite-difference methods in general.

/pdf/wave-equation-migration-with-the-phase-shift-method-34wz4vawgx.pdf

Wave equation migration with the phase-shift method

Wave equation migration is known to be simpler in principle when the horizontal coordinate or coordinates are replaced by their Fourier conjugates. Two practical migration schemes utilizing this concept are developed in this paper. One scheme extends the Claerbout finite difference method, greatly reducing dispersion problems usually associated with this method at higher dips and frequencies. The second scheme effects a Fourier transform in both space and time; by using the full scalar wave equation in the conjugate space, the method eliminates (up to the aliasing frequency) dispersion altogether. The second method in particular appears adaptable to three‐dimensional migration and migration before stack.

/pdf/migration-by-fourier-transform-2srmg05k7s.pdf

R. H. Stolt

Papers

Migration by fourier transform