Traveltime tomography in anisotropic media—I. Theory

SUMMARY In principle, crosshole, traveltime tomography is ideal for directly detecting and measuring seismic anisotropy. The traveltimes of multiple rays with wide angular coverage will be sensitive both to inhomogeneities and to anisotropy. In practice, the traveltimes will depend only on a limited number of the anisotropic velocity parameters, and the data may not be adequate even to determine these parameters uniquely. In addition, trade-offs may exist between anisotropy and inhomogeneities. In this paper, we use the linear perturbation theory for traveltimes in general, weakly anisotropic media to discuss the dependence of traveltimes in 2-D crosshole tomographic experiments on the anisotropic parameters. In a companion paper, we apply the results to synthetic and real data examples. We show that when measurements are restricted to a 2-D plane, the qP and qS traveltimes depend on subsets of the complete set of 21 anisotropic velocity parameters. Formulae are developed for the differential coefficients of the traveltimes with respect to these parameters in piecewise homogeneous and in linearly interpolated models. It is shown how in a generally oriented model element, the local parameters are related to the same parameters in the global model. The parameters that can be determined from 2-D tomographic data do not in general determine the full nature of the anisotropy. Rather, these parameters serve only to describe the intersection of the slowness sheet with the 2-D plane. Since many models may fit this description, additional information on symmetry properties and orientations is required. For example, if a priori information suggests that the anisotropy is transversely isotropic (TI), then we can determine some of the TI parameters and some information on the orientation of the axis of symmetry. Formulae are given relating the general parameters to those of a TI system with general orientation of the symmetry axis. The general formulae for qS traveltimes are intrinsically more complicated than those for qP. In the qS case, the traveltime perturbation depends on the polarization, which in turn depends on the perturbation. This makes the general problem non-linear even for small perturbations. However, the mean qS traveltime and the traveltime dependence on various subsets of parameters are linear. Although linear perturbation theory is invalid for qS rays, degenerate perturbation theory is valid for the calculation of the traveltimes and could be used in a non-linear inversion scheme.

[1]  M. Fehler,et al.  CROSS-HOLE SEISMIC SURVEYS: APPLICATIONS FOR STUDYING SUBSURFACE FRACTURE SYSTEMS AT A HOT DRY ROCK GEOTHEMAL SITE , 1983 .

[2]  G. F. West,et al.  Crosshole seismology and seismic imaging in crystalline rocks , 1983 .

[3]  R. G. Pratt,et al.  Traveltime tomography in anisotropic media—II. Application , 1992 .

[4]  L. Thomsen Weak elastic anisotropy , 1986 .

[5]  R. Lytle,et al.  Computerized geophysical tomography , 1979, Proceedings of the IEEE.

[6]  A. Love,et al.  The Mathematical Theory of Elasticity. , 1928 .

[7]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[8]  A Seismic-Ray Tracing Method for Two-Dimensional Inhomogeneous Media , 1976 .

[9]  V. Červený Seismic Rays and Ray Intensities in Inhomogeneous Anisotropic Media , 1972 .

[10]  Robert R. Stewart,et al.  An algebraic reconstruction technique for weakly anisotropic velocity , 1988 .

[11]  J. Virieux,et al.  Ray tracing for earthquake location in laterally heterogeneous media , 1988 .

[12]  Sven Ivansson,et al.  Tomographic velocity estimation in the presence of low velocity zones , 1984 .

[13]  M. Worthington,et al.  Some Sources of Distortion in Tomographic Velocity IMAGES1 , 1988 .

[14]  Donald F. Winterstein,et al.  Velocity anisotropy in shale determined from crosshole seismic and vertical seismic profile data , 1990 .

[15]  C. H. Chapman,et al.  Crosshole seismic tomography , 1989 .

[16]  John E. Peterson,et al.  Applications of algebraic reconstruction techniques to crosshole seismic data , 1985 .

[17]  G. Backus A geometrical picture of anisotropic elastic tensors , 1970 .

[18]  Donald F. Winterstein,et al.  Velocity anisotropy terminology for geophysicists , 1990 .

[19]  G. West,et al.  Cross-hole seismic scanning and tomography , 1987 .

[20]  G. West,et al.  Crosshole audiofrequency seismology in granitic rocks using piezoelectric transducers as sources and detectors , 1984 .

[21]  Chris Chapman,et al.  Ray theory and its extensions; WKBJ and Maslov seismograms , 1985 .

[22]  P. Bois,et al.  Well-to-well seismic measurements , 1972 .

[23]  I. M. Mason Algebraic reconstruction of a two‐dimensional velocity inhomogeneity in the High Hazles seam of Thoresby colliery , 1981 .

[24]  Jiří Jech,et al.  First-order perturbation method for anisotropic media , 1989 .