地震体波走时层析成像方法研究综述

赵烽帆; 张明辉; 徐涛

doi:10.6038/pg20140312

地震体波走时层析成像方法研究综述

- 1.
  中国地震台网中心, 北京 100045;
- 2.
  中国科学院地质与地球物理研究所, 岩石圈演化国家重点实验室 100029;
- 3.
  中国科学院大学, 北京 100039

详细信息

作者简介:
赵烽帆，女，1983年生，工程师，主要从事震源物理及地震活动性研究.（E-mail：weina22@163.com）

通讯作者: 徐涛，男，1978年生，副研究员，主要从事地震波传播和壳幔结构成像研究.（E-mail：xutao@mail.iggcas.ac.cn）

中图分类号: P631

收稿日期: 2014-02-20

修回日期: 2014-03-20

刊出日期: 2014-06-20

A review of body wave traveltime tomography methods

- 1.
  China Earthquake Networks Center, Beijing 100045, China;
- 2.
  State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 10029, China;
- 3.
  University of Chinese Academy of Sciences, Beijing 100039, China

More Information

Corresponding author: XU Tao

Received Date: 20 February 2014

Revised Date: 20 March 2014

Publish Date: 20 June 2014

摘要

摘要: 地球介质参数是地震深部探测的重要目标，蕴含了丰富的地球结构和演化信息.全波形反演技术由于计算量大的原因，目前阶段在实际壳幔结构成像应用中还存在较大的难度.基于P波和S波到时的体波成像仍然是壳幔精细速度结构成像的重要途径.本文主要介绍了三种具有代表性的地震体波走时层析成像方法，即基于高频近似射线理论的成像方法、基于程函方程数值解的初至波成像方法及有限频成像方法.基于射线理论的成像方法计算量小，计算效率高，长期以来是体波成像的主要方法之一；基于程函方程数值解的成像方法，在初至波走时成像方面有优势，但是要求网格节点较密，计算量大；有限频成像方法基于波动方程理论，比射线理论包含更丰富的介质信息，但是由于计算量大，目前在实现过程中引入简化的假设，成像效果还有较大的提升空间.
- 体波 /
- 走时层析成像 /
- 射线理论 /
- 程函方程 /
- 有限频
Abstract: The medium parameter of the Earth is an important target of the deep seismic exploration, which contains rich information of the structure and evolution of the Earth. Full waveform inversion technique is restricted in the application to the crust and mantle imaging due to the time-consuming computation. The body wave tomography based on the arrival time of P-wave and S-wave plays a key role in imaging fine velocity structure of the crust and mantle. This paper mainly introduces three typical seismic body wave traveltime tomography methods, which are traveltime tomography method based on ray theory, first-arrival traveltime tomography method based on eikonal equation numerical solver and finite frequency tomography method. The first one requires a small calculation amount and is efficient; the second one has the advantage in the first arrival traveltime imaging, while is more time-consuming; the last one which is based on the wave equation theory, has higher resolution than the method based on ray theory, while is more time-consuming and less efficient.
- body wave /
- traveltime tomography /
- ray theory /
- eikonal equation /
- finite frequency

HTML全文

参考文献(86)

[1]	Aki K, Lee W H K. 1976. Determination of three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model[J]. Journal of Geophysical Research, 81(23): 4381-4399.
[2]	Bai C, Li X, Huang G, et al. 2013. Simultaneous Inversion for Velocity and Reflector Geometry Using Multi-phase Fresnel Volume Rays[J]. Pure and Applied Geophysics, 1-17.
[3]	Cerveny V. 2001.Seismic Ray Theory[M].Cambridge University Press.
[4]	Cerveny V, Soares J E P. 1992. Fresnel volume ray tracing[J]. Geophysics, 57(7): 902-915.
[5]	Chen L H, Song Z H. 1990. Crust-upper mantle P wave velocity structure beneath Northern China[J]. Geophysics. (in Chinese), 33(5): 540-546.Dahlen F A, Hung S H, Nolet G. 2000. Fréchet kernels for finite-frequency traveltimes-I.Theory[J]. Geophysical Journal International, 141(1): 157-174.
[6]	Dahlen F A, Tromp J. 1998. Theoretical global seismology[M]. Princeton university press.
[7]	Ding Z F, He Z Q, Sun W G, et al. 1999. 3-D crust and upper mantle velocity structure in Eastern Tibetan plateau and its surrounding areas[J]. Geophysics. (in Chinese), 42(2): 197-205.
[8]	Duan Y H, Zhang X K, Fang S M. 2002. Three-dimensional finite-difference tomography of velocity structure of the upper crustal in North China[J]. Chinese J. Geophys. (in Chinese), 45(3): 362-369.
[9]	Engdahl E R, Lee W H K. 1976. Relocation of local earthquakes by seismic ray tracing[J]. Journal of Geophysical Research, 81(23): 4400-4406.
[10]	Favier N, Chevrot S. 2003. Sensitivity kernels for shear wave splitting in transverse isotropic media[J]. Geophysical Journal International, 153(1): 213-228.
[11]	Guo B, Liu Q Y, Chen J H, et al. 2009. Teleseismic P-wave tomography of the crust and upper mantle in Longmenshan area, west Sichuan[J]. Chinese J. Geophys. (in Chinese), 52(2): 346-355.
[12]	Hole, J.A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography[J]. Journal of Geophysical Research, 97(B5): 6553-6562.
[13]	Hung S H, Dahlen F A, Nolet G. 2000.Fréchet kernels for finite-frequency traveltimes-II.Examples[J]. Geophysical Journal International, 141(1): 175-203.
[14]	Jiang G M, Zhao D P, Zhang G B. 2008. Detailed structure of the subducting Pacific slab beneath the Japan islands and Japan sea[J]. Earth Science Frontiers (in Chinese), 15(1): 222-231.
[15]	Jiang G M, Zhao D P, Zhang G B. 2009. Crustal correction in teleseismic tomography and its application[J]. Chinese J. Geophys. (in Chinese), 52(6): 1508-1514.
[16]	Jiang Y, Chen X F. 2011. Review on the comparative study between finite-frequency tomography and ray-theoretical tomography[J]. Progress in Geophys. (in Chinese), 26(5): 1566-1575.
[17]	Jin A S, Liu F T, Sun Y Z. 1980. Three-dimensional P velocity structure of the crust and upper mantle under Beijing region[J]. Chinese J. Geophys. (in Chinese), 23(2), 172-182.
[18]	Julian B R, Gubbins D. 1977.Three-dimensional seismic ray tracing[J]. J. geophys, 43(1): 95-114.
[19]	Kayal J R, Zhao D. 1998.Three-dimensional seismic structure beneath Shillong Plateau and Assam Valley, northeast India[J]. Bulletin of the Seismological Society of America, 88(3): 667-676.
[20]	Komatitsch D, Liu Q, Tromp J, et al. 2004. Simulations of ground motion in the Los Angeles basin based upon the spectral-element method[J]. Bulletin of the Seismological Society of America, 94(1): 187-206.
[21]	Lan H, Zhang Z. 2013. Topography-dependent eikonal equation and its solver for calculating first-arrival traveltimes with an irregular surface[J]. Geophysical Journal International, 193(2): 1010-1026.
[22]	Lan H, Zhang Z. 2013. A High-Order Fast-Sweeping Scheme for Calculating First-Arrival Travel Times with an Irregular Surface[J]. Bulletin of the Seismological Society of America, 103(3): 2070-2082.
[23]	Langan, R.T., Lerche, I., Cutler, R.T. 1985. Tracing of rays through heterogeneous media: Anaccurate and efficient procedure[J]. Geophysics, 50(9): 1456-1465.
[24]	Lei J S, Zhou H L, Zhao D P. 2002. 3-D velocity structure of P-wave in the crust and upper-mantle beneath Pamir and adjacent region[J]. Chinese J. Geophys. (in Chinese), 45(6): 802-811.
[25]	Lei J S, Zhao D P, Su J R, et al. 2009.Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuanearthquake[J]. Chinese J. Geophys. (in Chinese), 52(2): 339-345.
[26]	Li, C., and van der Hilst, R.D. 2010.Structure of the upper mantle and transition zone beneath South East Asia from travel time tomography[J].Journal of Geophysical Research, 115, B07308, doi:10.1029/2009JB006882.
[27]	Li F, Xu T, Wu Z B, et al. 2013. Segmentally iterative ray tracing in 3-D heterogeneous geological models[J]. Chinese J. Geophys. (in Chinese), 56(10): 3514-3522.
[28]	Li F, Xu T, Zhang M, et al. 2013. Seismic traveltime inversion of 3D velocity model with triangulated interfaces[J]. Earthquake Science, 1-10.
[29]	Li Z W, Xu Y, Hao T Y, et al. 2006. Seismic tomography and velocity structure in the crust and upper mantle around Bohai Sea area[J]. Chinese J. Geophys. (in Chinese), 49(3): 797-804.
[30]	Liang X, Sandvol E, Chen Y J, et al. 2012.A complex Tibetan upper mantle: A fragmented Indian slab and no south-verging subduction of Eurasian lithosphere[J]. Earth and Planetary Science Letters, 333: 101-111.
[31]	Liu Q, Gu Y J. 2012. Seismic imaging: from classical to adjointtomography[J]. Tectonophysics, 566: 31-66.
[32]	Lv Q T, Jiang M, Ma K Y. 1996. 3-D traveltime inversion and the deep structure in the southern Tibetan plateau[J]. ActaSeismologicasinica (in Chinese), 18(4): 451-459.
[33]	Maarten V, van Der Hilst R D. 2005. On sensitivity kernels for 'wave-equation'transmissiontomography[J]. Geophysical Journal International, 160(2): 621-633.
[34]	Marquering, H., Dahlen, F.A., Nolet, G. 1999. Three-dimensionalsensitivity kernels for finite-frequency traveltimes: the banana-doughnut paradox[J]. Geophys. J. Int., 137(3): 805-815.
[35]	Montelli R, Nolet G, Dahlen F A, et al. 2004. Finite-frequency tomography reveals a variety of plumes in the mantle[J]. Science, 303(5656): 338-343.
[36]	Nolet G. 2008. A breviary of seismic tomography: Imaging the Interior[M].Cambridge University Press.
[37]	Pereyra V. 1992.Two-point ray tracing in general 3D media[J]. Geophysical Prospecting, 40(3): 267-287.
[38]	Pratt R G. 1999. Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model[J]. Geophysics, 64(3): 888-901.
[39]	Qin F, Luo Y, Olsen K B, et al. 1992. Finite-difference solution of the eikonal equation along expanding wavefronts[J]. Geophysics, 57(3): 478-487.
[40]	Rawlinson N, Pozgay S, Fishwick S. 2010. Seismic tomography: a window into deep Earth[J]. Physics of the Earth and Planetary Interiors, 178(3): 101-135.
[41]	Rawlinson N, Reading A M, Kennett B L N. 2006. Lithospheric structure of Tasmania from a novel form of teleseismictomography[J]. Journal of Geophysical Research: Solid Earth(1978-2012), 111(B2).
[42]	Rawlinson N, Sambridge M. 2004. Multiple reflection and transmission phases in complex layered media using a multistage fast marching method[J]. Geophysics, 69(5): 1338-1350.
[43]	Rawlinson N., ReadingA M, Kennett B L N. 2006. Lithospheric structure of Tasmania from a novel form of teleseismictomography[J]. Journal of Geophysical Research: Solid Earth(1978-2012),111(B2).
[44]	Rawlinson N, Pozgay S, Fishwick S. 2010. Seismic tomography: A window into deep Earth[J]. Physics of the Earth and Planetary Interiors, 178(3): 101-135.
[45]	Sethian J A. 1996. A fast marching level set method for monotonically advancing fronts[J]. Proceedings of the National Academy of Sciences, 93(4): 1591-1595.
[46]	Sethian J A, Popovici A M. 1999.3-D traveltime computation using the fast marching method[J]. Geophysics, 64(2): 516-523.
[47]	Shearer P M. 2009. Introduction to Seismology[M]. Cambridge University Press.
[48]	Teng J W, Liu F T, Quan Y L, et al. 1994. The crust and upper mantle seismic tomography in the orogenic belts and sedimentary basins of the Northwest China[J].Progress in solid Geophysics of China (in Chinese), 66-80.
[49]	Teng J W, Zhang Z J, Bai W M. 2003. Lithosphere Physics[J]. Beijing: Science Press.
[50]	Tian Y, Liu C, Feng X. 2011. P-wave velocity structure of crust and upper mantle in Northeast China and its control on the formation of mineral and energy[J]. Chinese J. Geophys. (in Chinese), 54(2): 407-414.
[51]	Tian Y, Zhao D P, Liu C, et al. 2009.A review of body-wave tomography and its applications to studying the crust and mantle structure in China[J]. Earth Science Frontiers (in Chinese), 16(2): 347-360.
[52]	Tian Y, Zhao D P, Sun R M, et al. 2007. The 1992 Landers earthquake: effect of crustal heterogeneity on earthquake generation[J]. Chinese J. Geophys. (in Chinese), 50(5): 1488-1496.
[53]	ThurberC H, Ellsworth W L. 1980. Rapid solution of ray tracing problems in heterogeneous media[J].Bulletin of the Seismological Society of America, 70(4): 1137-1148.
[54]	Thurber C H. 1983. Earthquake locations and three-dimensional crustal structure in the Coyote Lake area, central California[J]. Journal of Geophysical Research: Solid Earth (1978–2012), 88(B10): 8226-8236.
[55]	Tromp J, Komattisch D, Liu Q. 2008. Spectral-element and adjoint methods in seismology[J]. Communications in Computational Physics, 3(1): 1-32.
[56]	Um J, Thurber C. 1987.A fast algorithm for two-point seismic ray tracing[J]. Bulletin of the Seismological Society of America, 77(3): 972-986.
[57]	Van Der Hilst R D, Maarten V. 2005. Banana-doughnut kernels and mantle tomography[J]. Geophysical Journal International, 163(3): 956-961.
[58]	Van der Hilst R D, Widiyantoro S, Engdahl E R. 1997. Evidence for deep mantle circulation from global tomography[J]. Nature, 386(6625): 578-584.
[59]	Vidale J. 1988. Finite-difference calculation of travel times[J]. Bulletin of the Seismological Society of America, 78(6): 2062-2076.
[60]	Vidale J E. 1990. Finite-difference calculation of traveltimes in three dimensions[J]. Geophysics, 55(5): 521-526.
[61]	Virieux J, Farra V. 1991. Ray tracing in 3-D complex isotropic media: An analysis of the problem. Geophysics, 56(12): 2057-2069.
[62]	Wang H, Huang D C. 2000. Seismic tomography and its application to rock mass structure study[J]. Journal of engineering geology (in Chinese), 8(1): 109-117.
[63]	Woodward M J. 1992.Wave-equation tomography[J]. Geophysics, 57(1): 15-26.
[64]	Wu Z B, Xu T, Li F, et al. 2014. Acoustic full-waveform inversion in the finite frequency-domain[J]. See: Chen Y T, eds, the physics of the earth's interior and dynamics in China-to celebrate the 60th anniversary of engaging in the geophysical research of the academician TengJiwen (in Chinese), Beijing: Science press, 306-314.
[65]	Xu T, Li F, Wu Z B, et al. 2014. A successive three-point perturbation method for fast ray tracing in complex 2D and 3D geological models[J]. Tectonophysics, doI:10.1016/j.tecto.2014.02.012.
[66]	Xu T, Xu G M, Gao E G, et al. 2004. Block modeling and shooting ray tracing in complex 3-D media[J]. Chinese J. Geophys. (in Chinese), 47(6): 1118-1126.
[67]	Xu T, Xu G, Gao E, Li Y, Jiang X, Luo K. 2006. Block modeling and segmentally iterative ray tracing in complex 3D media[J]. Geophysics, 71(3): T41-T51.
[68]	Xu T, Zhang Z, Gao E, Xu G, Sun L. 2010. Segmentally iterative ray tracing in complex 2D and 3D heterogeneous block models[J]. Bull. Seismol. Soc. Am., 100(2): 841-850.
[69]	Xu Y, Liu F T, Liu J H, et al. 2001. Deep features of continental collision belts in northwestern China and their dynamic implications[J]. Chinese J. Geophys. (in Chinese), 44(1): 41-47.
[70]	Xu Y, Liu J H, Hao T Y, et al. 2006. P wave velocity structure and tectonics analysis of lithospheric mantle in eastern China seas and adjacent regions[J]. Chinese J. Geophys. (in Chinese), 49(4): 1053-1061.
[71]	Yang G H. 2009. Travel time tomography method with Fresnel volume and its applications[D]. Xiamen University.
[72]	Yang T, Shen Y. 2006.Frequency-dependent crustal correction for finite-frequency seismic tomography[J]. Bulletin of the Seismological Society of America, 96(6): 2441-2448.
[73]	Zelt C A, Smith R B. 1992. Seismic traveltime inversion for 2-D crustal velocity structure[J]. Geophysical journal international, 108(1): 16-34.
[74]	Zhang F X. 2011. Finite-frequency traveltime tomography by body wave and its application in the North China[D]. Institute of Geophysics, China Earthquake Administration.
[75]	Zhang F X, Li Y H, Wu Q J, et al. 2011. The P wave velocity structure of upper mantle beneath the North China and surrounding regions from FMTT[J]. Chinese J. Geophys. (in Chinese), 54(5): 1233-1242.
[76]	Zhang F X, Wu Q J, Li Y H. 2013. The traveltime tomography study by teleseismic P wave data in the Northeast China area[J]. Geophys. (in Chinese), 56(8): 2690-2700.
[77]	Zhao D, Christensen D, Pulpan H. 1995. Tomographic imaging of the Alaska subduction zone[J]. Journal of Geophysical Research: Solid Earth (1978–2012), 100(B4): 6487-6504.
[78]	Zhao D, Hasegawa A, Horiuchi S. 1992. Tomographic imaging of P and S wave velocity structure beneath northeastern Japan[J]. Journal of Geophysical Research: Solid Earth (1978–2012), 97(B13): 19909-19928.
[79]	Zhao D, Hasegawa A, Kanamori H. 1994. Deep structure of Japan subduction zone as derived from local, regional, and teleseismicevents[J]. Journal of Geophysical Research: Solid Earth (1978–2012), 99(B11): 22313-22329.
[80]	Zhao D, Kanamori H, Humphreys E. 1996. Simultaneous inversion of local and teleseismic data for the crust and mantle structure of southern California[J]. Physics of the earth and planetary interiors, 93(3): 191-214.
[81]	Zhao D, Wang Z, Umino N, et al. 2007. Tomographic imaging outside a seismic network: Application to the northeast Japan arc[J]. Bulletin of the Seismological Society of America, 97(4): 1121-1132.
[82]	Zhao L, Allen R M, Zheng T, et al. 2009. Reactivation of an Archean craton: Constraints from P-and S-wave tomography in North China[J]. Geophysical Research Letters, 36: L17306.
[83]	Zhao L, Allen R M, Zheng T, et al. 2012. High-resolution body wave tomography models of the upper mantle beneath eastern China and the adjacent areas[J]. Geochemistry, Geophysics, Geosystems, 13(6).
[84]	Zhao L, Jordan T H, Chapman C H. 2000. Three-dimensional Fréchet differential kernels for seismicdelaytimes[J]. Geophysical Journal International, 141(3): 558-576.
[85]	Zheng H W, Li T D, Gao R, et al. 2007. Teleseismic P-wave tomography evidence for the Indian lithospheric mantle subducting northward beneath the Qiangtangterrane[J]. Chinese J. Geophys. (in Chinese), 50(5): 1418-1426.
[86]	Zheng H W, Li T D, Gao R. 2011. The progress of seismic tomography research in the Tibetan plateau[J]. Geophysical and geochemical exploration, 35(2): 160-164.

施引文献

资源附件(0)

访问统计

计量

文章访问数:
PDF下载数:
施引文献: 0

地震体波走时层析成像方法研究综述

作者简介: 赵烽帆，女，1983年生，工程师，主要从事震源物理及地震活动性研究.（E-mail：weina22@163.com）

通讯作者: 徐涛，男，1978年生，副研究员，主要从事地震波传播和壳幔结构成像研究.（E-mail：xutao@mail.iggcas.ac.cn）

A review of body wave traveltime tomography methods

Corresponding author: XU Tao

计量

出版历程

目录

作者简介:
赵烽帆，女，1983年生，工程师，主要从事震源物理及地震活动性研究.（E-mail：weina22@163.com）