Since the appearance of wave-equation migration, many have tried to improve the resolution and effectiveness of this technology. Least-squares wave-equation migration is one of those attempts that tries to fill the gap between the migration assumptions and reality in an iterative manner. However, these iterations do not come cheap. A proven solution to limit the number of least-squares iterations is to correct the gradient direction within each iteration via the action of a preconditioner that approximates the Hessian inverse. However, the Hessian computation, or even the Hessian approximation computation, in large-scale seismic imaging problems involves an expensive computational bottleneck, making it unfeasible. Therefore, we propose an efficient computation of the Hessian approximation operator, in the context of one-way wave-equation migration (WEM) in the space-frequency domain. We build the Hessian approximation operator depth by depth, considerably reducing the operator size each time it is calculated. We prove the validity of our proposed method with two numerical examples. We then extend our proposal to the framework of full-wavefield migration, which is based on WEM principles but includes interbed multiples. Finally, this efficient preconditioned least-squares full-wavefield migration is successfully applied to a dataset with strong interbed multiple scattering.

@en

Since seismic imaging creates an image of the subsurface structure based on information received from the measured wavefield, it is essential to fully utilize the reflected waves. Full Wavefield Modeling (FWMod) was developed with recursive and iterative up-and-down wavefield propagation, using one-way wave propagation, to model both primary and multiple reflections. Using FWMod as the modeling engine, Full Wavefield Migration (FWM) has been introduced to directly image data including internal multiples, where internal multiple crosstalk is suppressed automatically via an inversion-based data-fitting process. This avoids the need for applying internal multiple removal, which is often challenging. Conventional one-way wave propagators calculated in the wavenumber domain, like the phase shift (PS) operator, have limitations when applied to strongly inhomogeneous media. Even when computing a new operator at each lateral grid point, they still suffer difficulties because the medium is assumed to be locally homogeneous. In the past, matrix eigendecomposition has been proposed as a way to create accurate, local velocity-based one-way propagation operators. In this article, an accurate propagator based on eigendecomposition is incorporated into FWMod and FWM. In the numerical examples, four models with strong lateral velocity variations were used to test the propagator. With a comparison of the conventional FWM based on the PS operator with input data including FWMod and a finite-difference (FD) approach, the numerical examples demonstrated that the proposed method has the potential to significantly enhance image reflectivity, suppress internal multiples, and maintain convergence speed during the least-squares inversion.

@en

Enhanced pre-stack depth migration, characterized by improved resolution and amplitudes, ensures a more accurate representation of the subsurface, proving essential for reducing the likelihood of geological misinterpretations and facilitating informed decision-making in seismic exploration. However, obtaining high-resolution images with preserved amplitudes through standard depth migration could face several hurdles known as migration artifacts. Iterative least-squares migration (LSM) was developed to address these migration artifacts. However, the convergence rate of LSM using a gradient descent approach tends to be slow. Several researchers have attempted to achieve computational efficiency in linearized LSM through gradient preconditioning. In the context of iterative least-squares wave-equation migration, this extended abstract compares three minimization approaches that differ in error functions and gradient preconditioning—including the depth-dependent Hessian approximation inverse—through two numerical examples, one with an inverse-crime scenario and the other with a non-inverse-crime scenario.@en

Seismic imaging is crucial for subsurface exploration and monitoring, with a focus on deep and complex structures. Seismic wave migration solves the wave equation, and an accurate propagator is essential. Full Wavefield Modeling (FWMod) was developed based on recursive and iterative up/down wavefield propagation, modeling both primaries and multiples. Embedded within Full Wavefield Migration (FWM) it can be used to image data including multiples, resulting in better illumination in case primary illumination is not sufficient. FWM can be efficient and effective, but conventional one-way wave operators, such as Phase Shift Plus Interpolation Migration, have limitations in strongly inhomogeneous media. Local velocity-based one-way operator based on eigen decomposition was proposed and integrated within FWMod and FWM in this study, improving image amplitudes and fidelity and improving converage speed in the least-squares inversion process.@en

Reflection waveform inversion (RWI) is a method that relies on primary pure reflection data to recover the subsurface background velocity based on the associated evolving seismic images. Background velocity updates estimated by conventional RWI are nonoptimal, which is partly attributed to low-resolution tomographic wavepaths and migration isochrones. Preconditioning RWI sensitivity kernels using Hessian information solves this problem but is not practical for a large number of model parameters. One-way reflection waveform inversion (ORWI) is a reflection waveform tomography technique in which the forward modeling scheme operates in one direction (downward and then upward) via virtual parallel data levels in the medium. The ORWI framework allows us to break down the Hessian matrix into smaller operators, which makes the preconditioning operation more efficient and less computationally expensive. This extended abstract turns conventional ORWI into a high-resolution but computationally feasible ORWI (Gauss-Newton ORWI) to improve the nonoptimal background velocity updates.@en

Conventional full waveform inversion has proven ineffective in recovering background velocity models for the targets out of the reach of refracted and diving waves. Therefore, different reflection waveform inversion (RWI) tools have been presented so far, one of which is joint migration inversion (JMI). JMI, unlike most similar tools, is parameterized by two classes of parameters, reflectivity and velocity, and equipped with decoupled imaging and tomographic sensitivity kernels. The fundamental characteristic shared by JMI or any RWI technique is the arrival-time consistency between the anchor parts of the modeled and observed data. This characteristic helps and keeps any reflection waveform tomography tool robust against cycle-skipping in near-offset data. In this paper, through examining the tomographic gradient of JMI, we show that conventional JMI suffers from cycle skipping in near-offset traces, which degenerates the fidelity of the tomographic gradient of JMI. Next, we demonstrate how the degeneracy can be lifted from two different perspectives, either by muting the misguiding tomographic wavepaths or migration isochrones built by near-offset residuals. Coupled with this, we show how excluding the cycle-skipped far-offset data significantly improves the validity of the tomographic gradient of JMI.

@en

Since the appearance of wave-equation migration, many have tried to push this technology to the limits of its capacity. Least-squares wave-equation migration is one of those attempts that tries to fill the gap between the migration assumptions and reality in an iterative manner. However, these iterations do not come cheap. A proven solution to limit the number of iterations is to correct the update direction in every iteration via a proper but cheap preconditioner that approximates the second-order partial derivatives, known as Hessian, of the data-error functional employed in the least-squares formulation. In this study, we will address this issue, namely the cheap computation of the Hessian inverse in the context of our one-way wave-equation-based migration method, cited as full wavefield migration. To do so, we will build and apply the Hessian inverse matrix depth by depth, reducing the Hessian matrix size considerably each time it is calculated. We prove the validity of our proposed method by two numerical examples.@en

Considering the steadily declining prices in the oil and gas industry, nowadays, the requirement for geophysical information becomes more important in order to get the most out of available reservoirs. Electromagnetic measurements could complement seismic data where it lacks information. The EM response to fluid fill complements the resolving power of seismic data. In the current study, we use a probabilistic method to estimate reservoir parameters individually and jointly through two simplistic, synthetic, 2D reservoir models which can be considered as the geometrical limits of water-oil-contacts in oil and gas fields. We demonstrate a constructive contribution of the measurements with different physical natures in the estimation of reservoir parameters.@en

In this paper we compared time-domain FWI problem vs. frequency-domain in terms of the behavior of cost function and model misfit error on a synthetic model with transmission geometry. Synthetic data was generated by solution of 2D acoustic wave equation and optimization was performed using L-BFGS algorithm for both approaches in time and frequency-domain. Although FWI problem has a better expression in frequency-domain but inversion of selected frequency components is drawback of this approach despite its lower computational costs in term of memory. This paper shows the difference of full band data inversion in the time-domain with sequential inversion of some selected frequency components in the frequency-domain numerically. According to the results timedomain FWI leads to better results due to its full band characteristics when computational resources permits us to use it.@en