Inversion of walkaway VSP data in the presence of lateral velocity heterogeneity

Abstract

Multi-azimuth walkaway vertical seismic profiling is an established technique for the estimation of in situ slowness surfaces and inferring anisotropy parameters. Normally, this technique requires the assumption of lateral homogeneity, which makes the horizontal slowness components at depths of downhole receivers equal to those measured at the surface. Any violations of this assumption, such as lateral heterogeneity or nonzero dip of intermediate interfaces, lead to distortions in reconstructed slowness surfaces and, consequently, to errors in estimated anisotropic parameters. In this work, we relax the assumption of lateral homogeneity and discuss how to correct vertical seismic profile data for weak lateral heterogeneity. We describe a procedure of downward continuation of recorded traveltimes that accounts for the presence of both vertical inhomogeneity and weak lateral heterogeneity, which produces correct slowness surfaces at depths of downhole receivers, noticing that sufficiently dense receiver coverage along a borehole is required to separate influences of vertical and lateral heterogeneity on measured traveltimes and obtain accurate estimates of the slowness surfaces. Once the slowness surfaces are found and a desired type of anisotropic model to be inverted is selected, the corresponding anisotropic parameters, providing the best fit to the estimated slownesses, can be obtained. We invert the slowness surfaces of P-waves for parameters of the simplest anisotropic model describing dipping fractures (transversely isotropic medium with a tilted symmetry axis). Five parameters of this model, namely, the P-wave velocity V0 in the direction of the symmetry axis, Thomsen's anisotropic coefficients ϵ and δ, the tilt ν, and the azimuth β of the symmetry axis, can be estimated in a stable manner when maximum source offset is greater than half of receiver depth.