Nonlinear refraction traveltime tomography
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
Cai, W., and Qin, F., 1994, Three‐dimensional refraction imaging: 63rd Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstract, 629–632.
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.
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.
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.
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.
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
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.