Nonlinear refraction traveltime tomography

Geophysics - Tập 63 Số 5 - Trang 1726-1737 - 1998
Jie Zhang1,2, M. Nafi Toksöz3,4
1Earth Resources Lab., Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts 02142.
2GeoTomo LLC, Louisville, CO, United States
3Earth Resources Lab., Department of Earth, Atmospheric, and Planetary Sciences, E34-462 MIT, 42 Carleton Street, Cambridge, Massachusetts 02142.
4Massachussets Institute of Technology, United States

Tóm tắt

A few important issues for performing nonlinear refraction traveltime tomography have been identified. They include the accuracy of the traveltime and raypath calculations for refraction, the physical information in the refraction traveltime curves, and the characteristics of the refraction traveltime errors. Consequently, we develop a shortest path ray‐tracing method with an optimized node distribution that can calculate refraction traveltimes and raypaths accurately in any velocity model. We find that structure ambiguity caused by short and long rays in the seismic refraction method may influence the inversion solution significantly. Therefore, we pose a nonlinear inverse problem that explicitly minimizes the misfits of the average slownesses (ratios of traveltimes to the corresponding ray lengths) and the apparent slownesses (derivatives of traveltimes with respect to distance). As a result, we enhance the resolution as well as the convergence speed. To regularize our inverse problem, we use the Tikhonov method to avoid solving an ill‐posed inverse problem. Errors in refraction traveltimes are characterized in terms of a common‐shot error, a constant deviation for recordings from the same shot, and a relative traveltime‐gradient error with zero mean with respect to the true gradient of the traveltime curve. Therefore, we measure the uncertainty of our tomography solution using a nonlinear Monte Carlo approach and compute the posterior model covariance associated with two different types of random data vectors and one random model vector. The nonlinear uncertainty analysis indicates that the resolution of a tomography solution may not correspond to the ray coverage. We apply this tomography technique to image the shallow velocity structure at a coastal site near Boston, Massachusetts. The results are consistent with a subsequent drilling survey.

Từ khóa


Tài liệu tham khảo

Ammon C. J., 1993, Bull. Seis. Soc. Am., 83, 509, 10.1785/BSSA0830020509

10.1098/rsta.1970.0005

10.1098/rsta.1970.0005

Cai, W., and Qin, F., 1994, Three‐dimensional refraction imaging: 63rd Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstract, 629–632.

10.1029/92JB02441

10.1190/1.1443489

10.1111/j.1365-246X.1995.tb05723.x

10.1111/j.1365-246X.1994.tb03293.x

Kutrubes, D. L., Zhang, J., and Hager, J., 1996, Comparison of conventional processing techniques and nonlinear refraction traveltime tomography for surveys at eastern Massachusetts coastal site: Proc., SAGEEP, Environ. Eng. Geophys. Soc., 135–138.

10.1190/1.1443407

Matarese, J. R., 1993, Nonlinear traveltime tomography: Ph.D. thesis, Mass. Inst. Tech.

McMechan G. A., 1987, Bull. Seis. Soc. Am., 77, 1945

Morgan, F. D. O., 1981, Electronics of sufide minerals: implications for induced polarization: Ph.D. thesis, Mass. Inst. Tech.

Moser, T. J., 1989, Efficient seismic ray tracing using graph theory: 59th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstract, 1106–1031.

1991, Shortest path calculation of seismic rays: Geophysics,56, 59–67.

10.4294/jpe1952.34.195

Nolet, G., 1987, Waveform tomography,inG. Nolet, Ed., Seismic tomography: D. Reidel Publ. Co., 301–323.

Palmer, D., 1980, The generalized reciprocal method of seismic refraction interpretation: Soc. Expl. Geophys.

Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., 1992, Numerical recipes in Fortran: the art of scientific computing, 2nd edition: Cambridge Univ. Press.

10.1190/1.1443263

Rodi, W. L., 1989, Regularization and Backus‐Gilbert estimation in nonlinear inverse problems: application to magnetotellurics and surface waves: Ph.D. thesis, Pennsylvania State Univ.

Saito, H., 1989, Traveltimes and raypaths of first arrival seismic waves: Computation method based on Huygens’ Principle: 59th Ann. Internat. Mtg., Soc Expl. Geophys., Expanded Abstract, 244–247.

1990, 3-D ray‐tracing method based on Huygens’ principle: 60th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstract, 1024–1027.

10.1190/1.1442293

10.1088/0266-5611/6/1/011

10.1190/1.1443923

Tarantola, A., 1987, Inverse problem theory: Elsevier.

Tikhonov, A. N., and Arsenin, V. Y., 1977. Solutions of ill‐posed problems: W. H. Winston and Sons.

Vidale J. E., 1988, Bull. Seis. Soc. Am., 78, 2062

10.1190/1.1443499

10.1111/j.1365-246X.1989.tb00498.x

10.1190/1.1443868

Zhang, J., Rodi, W., Mackie, R. L., and Shi, W., 1996, Regularization in 3-D dc resistivity tomography: Proc. SAGEEP, Environ. Eng. Geophy. Soc., 130–135.

10.1002/ima.1850010103

Thông tin
Thông tin xuất bản

Geophysics

Tập 63 Số 5

1726-1737

Thông tin tác giả