Direct traveltime inversion of VSP data for elliptical anisotropy in layered media
This paper presents an exact method for inversion of VSP data for anisotropic parameters that uses traveltimes directly, rather than phase slownesses computed from traveltime differences. There is a tradeoff in the present traveltime-based method compared with the slowness-based method. In our method we have to assume some function representing the vertical velocity variation, something that is not necessary in the phase-slowness method. There is, however, a higher numerical accuracy in the direct traveltime method. This is due partly to our using observed traveltimes directly in the computations rather than taking differences between values of comparable size, which greatly magnifies the relative errors. Error analysis shows that there also are other intrinsic reasons why our technique has less error. A numerical example representing a typical VSP yielded errors from the phase-slowness method (requiring two sources and two receivers) that were about 18 times larger than for a single determination by the direct traveltime method. Furthermore, in cases where the VSP has been acquired with a single source offset, e.g., for some offshore wells, the traveltime method will yield results, whereas the slowness method will not. INTRODUCTION In the last decade or two there has been a greatly increased recognition that significant seismic anisotropy is rather common in sedimentary rocks. It has also been shown that neglecting even a seemingly modest degree of anisotropy can lead to significant degradation in processed seismic images (e.g., Winterstein, 1986; Wright, 1987; Larner, 1993; Tsvankin, 1996; Chen and Castagna, 2000) or to significant errors in vertical or lateral positioning of subsurface features (e.g. Banik, 1984, 1987; Isaac and Lawton, 1999; Vestrum et al., 1999). Gaiser (1990) determined how to estimate the anisotropic parameters of a transversely isotropic medium from VSP data using vertical and horizontal phaseslowness measurements. Miller and Spencer (1994) and Miller et al. (1994) presented methods for inverting phase slowness estimates from walkaway VSPs for anisotropy parameters. Using a similar method, Leaney et al. (1999) have inverted multiazimuthal multioffset VSP P-wave data over a marine carbonate reservoir to determine azimuthal anisotropy in the carbonate and transverse isotropy in the overlying shale. Leslie and Lawton (1999) present a method for inverting refraction seismic data for the anisotropic parameters of steeply dipping strata. Tsvankin and Thomsen (1994) and Alkhalifah and Tsvankin (1995) deal with the problem of inverting for longer spreads, beyond the near-vertical elliptic-anisotropic approximation, in order to recover velocity fields for processing. * Authors listed in alphabetical order.
Analysis of multiazimuthal VSP data for anisotropy and AVO
Multioffset vertical seismic profile (VSP) experiments, commonly referred to as walkaways, enable anisotropy to be measured reliably in the field. The results can be fed into modeling programs to study the impact of anisotropy on velocity analysis, migration, and amplitude versus offset (AVO). Properly designed multioffset VSPs can also provide the target AVO response measured under optimum conditions, since the wavelet is recorded just above the reflectors of interest with minimal reflection point dispersal. In this paper, the multioffset VSP technique is extended to include multioffset azimuths, and a multiazimuthal multiple VSP data set acquired over a carbonate reservoir is analyzed for P-wave anisotropy and AVO. Direct arrival times down to the overlying shale and reflection times and amplitudes from the carbonate are analyzed. Data analysis involves a three‐term fit to account for nonhyperbolic moveout, dip, and azimuthal anisotropy. Results indicate that the overlying shale is transversely isotropi...
Determination of anisotropic velocity models from walkaway VSP data acquired in the presence of dip
Wide‐aperture walkaway vertical seismic profile (VSP) data acquired through transversely isotropic horizontal layers can be used to determine the P phase‐slowness surface, local to a receiver array in a borehole. In the presence of dip, errors in the slowness surface may occur if the medium is assumed to be layered horizontally. If the acquisition plane is oriented parallel to the dip direction, the derived slowness is too large for sources offset from the well in the down‐dip direction and too small for sources offset from the well in the up‐dip direction. For acquisition parallel to the strike of the layers, the recovery of the P phase‐slowness in the vicinity of the receiver array is excellent. It is therefore preferable to orient the walkaway VSP in the strike direction to estimate the anisotropic parameters of the medium in the vicinity of a receiver array. However, this may not be possible. If the dip direction of all layers has the same azimuth, the variation of walkaway traveltimes with azimuth ha...