Image-domain wavefield tomography for VTI media

Abstract

Processing algorithms for transversely isotropic (TI) media are widely used in depth imaging and typically bring substantial improvements in reflector focusing and positioning. Here, we develop acoustic image-domain tomography (IDT) for reconstructing VTI (TI with a vertical symmetry axis) models from P-wave reflection data. The modeling operator yields an integral wave-equation solution, which is based on a separable dispersion relation and contains only P-waves. The zero-dip NMO velocity ([Formula: see text]) and anellipticity parameter [Formula: see text] are updated by focusing energy in space-lag images obtained by least-squares reverse-time migration (LSRTM). Application of LSRTM helps mitigate aperture- and illumination-induced artifacts in space-lag gathers and improve the robustness of [Formula: see text]-estimation. The impact of the trade-off between [Formula: see text] and [Formula: see text] is reduced by a three-stage inversion algorithm that gradually relaxes the constraints on the spatial variation of [Formula: see text]. Assuming that the depth profile of the Thomsen parameter [Formula: see text] is known at two or more borehole locations, we employ image-guided interpolation to constrain the depth scale of the parameter fields and of the migrated image. Image-guided smoothing is also applied to the IDT gradients to facilitate convergence towards geologically plausible models. The algorithm is tested on synthetic reflection and borehole data from the structurally complicated elastic VTI Marmousi-II model. Although the initial velocity field is purely isotropic and substantially distorted, all three relevant parameters ([Formula: see text], [Formula: see text], and [Formula: see text]) are estimated with sufficient accuracy. The algorithm is also applied to a line from a 3D ocean-bottom-node data set acquired in the Gulf of Mexico.