Abstract

Velocity macromodel building is an essential step of the seismic imaging workflow. Indeed, obtaining acceptable results through migration or full waveform inversion is highly dependent on the kinematic accuracy of the background/initial velocity model. Two decades ago, stereotomography was proposed as an alternative to reflection traveltime tomography, the first relying on semi-automatic picking of locally coherent events associated with small reflection or diffraction segments tied to scatterers in depth by a pair of rays, while the latter on interpretive picking of laterally continuous reflections. The flexibility of stereotomography paved the way for many developments that have shown the efficiency of the method whilst emphasizing on the complementary information carried out by traveltimes and slopes of locally coherent events. A recent formulation recast stereotomography under a matrix-free formulation based on eikonal solvers and the adjoint-state method. In the latter, like in the previous works, the scatterer positions and the velocity field are updated jointly to tackle the ill-famed velocity–position coupling in reflection tomography. Following on from this adjoint-state formulation, we propose a new parsimonious formulation of slope tomography that offers the chance to restrain the problem to minimizing the residuals of a single data class being a slope, in search of a sole parameter class being the subsurface velocity field. This parsimonious formulation results from a variable projection, which is implemented by enforcing a consistency between the scatterer coordinates and the velocity macromodel through migration of kinematic attributes. We explain why the resulting reduced-parametrization inversion is more suitable for tomographic problems than the most common joint inversion strategy. We benchmark our method against the complex Marmousi model along with a validation through time domain full waveform inversion and then present the results of a field data case study.

Highlights

  • The key purpose of seismic imaging methods is the retrieval of the subsurface properties like for instance wave speeds, density, attenuation or anisotropy

  • The non-uniqueness of the solution in first-arrival traveltime tomography motivated the use of additional kinematic attributes as in reflection tomography (Bishop et al 1985; Farra & Madariaga 1988) or even higher order attributes relying on directional reception as in polarization tomography (Hu et al 1994; Farra & Le Begat 1995) where the wholeness of the slowness vector is exploited in a transmission regime

  • We review the parametrization of the classical stereotomography and the slope tomography based on eikonal solvers and the adjoint-state method (AST)

Read more

Summary

Introduction

The key purpose of seismic imaging methods is the retrieval of the subsurface properties like for instance wave speeds, density, attenuation or anisotropy. Building a kinematically accurate smooth velocity model of the subsurface from acquired seismic data is essential for obtaining reliable depth-migrated images (Etgen et al 2009) or adequate starting models for full-waveform inversion (FWI; Tarantola 1984; Virieux & Operto 2009). Several authors addressed this ill-posed inverse problem; most adopted the asymptotic high-frequency approximation (Cerveny 2001) whilst utilizing different types of data and methods. In reflection settings, Controlled Directional Reception (CDR; Rieber 1936; Riabinkin 1957; Sword 1987), relying on locally coherent events defined by traveltime and its first-order derivative, proposed an interesting approach in the sense that

Objectives
Methods
Findings
Discussion
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call