地震体波走时层析成像方法研究综述
详细信息    查看全文 | 推荐本文 |
  • 英文篇名:A review of body wave traveltime tomography methods
  • 作者:赵烽帆 ; 张明辉 ; 徐涛
  • 英文作者:ZHAO Feng-fan;ZHANG Ming-hui;XU Tao;China Earthquake Networks Center;State Key Laboratory of Lithospheric Evolution,Institute of Geology and Geophysics,Chinese Academy of Sciences;University of Chinese Academy of Sciences;
  • 关键词:体波 ; 走时层析成像 ; 射线理论 ; 程函方程 ; 有限频
  • 英文关键词:body wave;;traveltime tomography;;ray theory;;eikonal equation;;finite frequency
  • 中文刊名:DQWJ
  • 英文刊名:Progress in Geophysics
  • 机构:中国地震台网中心;中国科学院地质与地球物理研究所岩石圈演化国家重点实验室;中国科学院大学;
  • 出版日期:2014-06-15
  • 出版单位:地球物理学进展
  • 年:2014
  • 期:v.29;No.125
  • 基金:中国地震局地震科技星火计划青年项目(XH14055Y);; 中国地震局监测预报司震情跟踪合同制定向工作任务(2014020412);; 国家自然科学基金(41174075)联合资助
  • 语种:中文;
  • 页:DQWJ201403012
  • 页数:12
  • CN:03
  • ISSN:11-2982/P
  • 分类号:126-137
摘要
地球介质参数是地震深部探测的重要目标,蕴含了丰富的地球结构和演化信息.全波形反演技术由于计算量大的原因,目前阶段在实际壳幔结构成像应用中还存在较大的难度.基于P波和S波到时的体波成像仍然是壳幔精细速度结构成像的重要途径.本文主要介绍了三种具有代表性的地震体波走时层析成像方法,即基于高频近似射线理论的成像方法、基于程函方程数值解的初至波成像方法及有限频成像方法.基于射线理论的成像方法计算量小,计算效率高,长期以来是体波成像的主要方法之一;基于程函方程数值解的成像方法,在初至波走时成像方面有优势,但是要求网格节点较密,计算量大;有限频成像方法基于波动方程理论,比射线理论包含更丰富的介质信息,但是由于计算量大,目前在实现过程中引入简化的假设,成像效果还有较大的提升空间.
        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.
引文
Aki K,Lee W H K.1976.Determination of three-dimensionalvelocity anomalies under a seismic array using first P arrivaltimes from local earthquakes:1.A homogeneous initial model[J].Journal of Geophysical Research,81(23):4381-4399.
    Bai C,Li X,Huang G,et al.2013.Simultaneous Inversion forVelocity and Reflector Geometry Using Multi-phase FresnelVolume Rays[J].Pure and Applied Geophysics,1-17.
    Cerveny V.2001.Seismic Ray Theory[M].Cambridge UniversityPress.
    Cerveny V,Soares J E P.1992.Fresnel volume ray tracing[J].Geophysics,57(7):902-915.
    Chen L H,Song Z H.1990.Crust-upper mantle P wave velocitystructure beneath Northern China[J].Geophysics.(inChinese),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.
    Dahlen F A,Tromp J.1998.Theoretical global seismology[M].Princeton university press.
    Ding Z F,He Z Q,Sun W G,et al.1999.3-D crust and uppermantle velocity structure in Eastern Tibetan plateau and itssurrounding areas[J].Geophysics.(in Chinese),42(2):197-205.
    Duan Y H,Zhang X K,Fang S M.2002.Three-dimensional finitedifference tomography of velocity structure of the upper crustalin North China[J].Chinese J.Geophys.(in Chinese),45(3):362-369.
    Engdahl E R,Lee W H K.1976.Relocation of local earthquakes byseismic ray tracing[J].Journal of Geophysical Research,81(23):4400-4406.
    Favier N,Chevrot S.2003.Sensitivity kernels for shear wavesplitting in transverse isotropic media[J].Geophysical JournalInternational,153(1):213-228.
    Guo B,Liu Q Y,Chen J H,et al.2009.Teleseismic P-wavetomography of the crust and upper mantle in Longmenshanarea,west Sichuan[J].Chinese J.Geophys.(in Chinese),52(2):346-355.
    Hole,J.A.1992.Nonlinear high-resolution three-dimensionalseismic travel time tomography[J].Journal of GeophysicalResearch,97(B5):6553-6562.
    Hung S H,Dahlen F A,Nolet G.2000.Fréchet kernels for finitefrequency traveltimes-II.Examples[J].Geophysical JournalInternational,141(1):175-203.
    Jiang G M,Zhao D P,Zhang G B.2008.Detailed structure of thesubducting Pacific slab beneath the Japan islands and Japan sea[J].Earth Science Frontiers(in Chinese),15(1):222-231.
    Jiang G M,Zhao D P,Zhang G B.2009.Crustal correction inteleseismic tomography and its application[J].Chinese J.Geophys.(in Chinese),52(6):1508-1514.
    Jiang Y,Chen X F.2011.Review on the comparative study betweenfinite-frequency tomography and ray-theoretical tomography[J].Progress in Geophys.(in Chinese),26(5):1566-1575.
    Jin A S,Liu F T,Sun Y Z.1980.Three-dimensional P velocitystructure of the crust and upper mantle under Beijing region[J].Chinese J.Geophys.(in Chinese),23(2),172-182.
    Julian B R,Gubbins D.1977.Three-dimensional seismic ray tracing[J].J.geophys,43(1):95-114.
    Kayal J R,Zhao D.1998.Three-dimensional seismic structurebeneath Shillong Plateau and Assam Valley,northeast India[J].Bulletin of the Seismological Society of America,88(3):667-676.
    Komatitsch D,Liu Q,Tromp J,et al.2004.Simulations of groundmotion in the Los Angeles basin based upon the spectral-elementmethod[J].Bulletin of the Seismological Society of America,94(1):187-206.
    Lan H,Zhang Z.2013.Topography-dependent eikonal equation andits solver for calculating first-arrival traveltimes with anirregular surface[J].Geophysical Journal International,193(2):1010-1026.
    Lan H,Zhang Z.2013.A High-Order Fast-Sweeping Scheme forCalculating First-Arrival Travel Times with an Irregular Surface[J].Bulletin of the Seismological Society of America,103(3):2070-2082.
    Langan,R.T.,Lerche,I.,Cutler,R.T.1985.Tracing of raysthrough heterogeneous media:Anaccurate and efficientprocedure[J].Geophysics,50(9):1456-1465.
    Lei J S,Zhou H L,Zhao D P.2002.3-D velocity structure of Pwave in the crust and upper-mantle beneath Pamir and adjacentregion[J].Chinese J.Geophys.(in Chinese),45(6):802-811.
    Lei J S,Zhao D P,Su J R,et al.2009.Fine seismic structure underthe Longmenshan fault zone and the mechanism of the largeWenchuanearthquake[J].Chinese J.Geophys.(in Chinese),52(2):339-345.
    Li,C.,and van der Hilst,R.D.2010.Structure of the uppermantle and transition zone beneath South East Asia from traveltime tomography[J].Journal of Geophysical Research,115,B07308,doi:10.1029/2009JB006882.
    Li F,Xu T,Wu Z B,et al.2013.Segmentally iterative ray tracingin 3-D heterogeneous geological models[J].Chinese J.Geophys.(in Chinese),56(10):3514-3522.
    Li F,Xu T,Zhang M,et al.2013.Seismic traveltime inversion of3Dvelocity model with triangulated interfaces[J].EarthquakeScience,1-10.
    Li Z W,Xu Y,Hao T Y,et al.2006.Seismic tomography andvelocity structure in the crust and upper mantle around BohaiSea area[J].Chinese J.Geophys.(in Chinese),49(3):797-804.
    Liang X,Sandvol E,Chen Y J,et al.2012.A complex Tibetanupper mantle:A fragmented Indian slab and no south-vergingsubduction of Eurasian lithosphere[J].Earth and PlanetaryScience Letters,333:101-111.
    Liu Q,Gu Y J.2012.Seismic imaging:from classical toadjointtomography[J].Tectonophysics,566:31-66.
    Lv Q T,Jiang M,Ma K Y.1996.3-D traveltime inversion and thedeep structure in the southern Tibetan plateau[J].ActaSeismologicasinica(in Chinese),18(4):451-459.
    Maarten V,van Der Hilst R D.2005.On sensitivity kernels for‘wave-equation’transmissiontomography[J].GeophysicalJournal International,160(2):621-633.
    Marquering,H.,Dahlen,F.A.,Nolet,G.1999.Threedimensionalsensitivity kernels for finite-frequency traveltimes:the banana-doughnut paradox[J].Geophys.J.Int.,137(3):805-815.
    Montelli R,Nolet G,Dahlen F A,et al.2004.Finite-frequencytomography reveals a variety of plumes in the mantle[J].Science,303(5656):338-343.
    Nolet G.2008.A breviary of seismic tomography:Imaging theInterior[M].Cambridge University Press.
    Pereyra V.1992.Two-point ray tracing in general 3D media[J].Geophysical Prospecting,40(3):267-287.
    Pratt R G.1999.Seismic waveform inversion in the frequencydomain,Part 1:Theory and verification in a physical scalemodel[J].Geophysics,64(3):888-901.
    Qin F,Luo Y,Olsen K B,et al.1992.Finite-difference solution ofthe eikonal equation along expanding wavefronts[J].Geophysics,57(3):478-487.
    Rawlinson N,Pozgay S,Fishwick S.2010.Seismic tomography:awindow into deep Earth[J].Physics of the Earth and PlanetaryInteriors,178(3):101-135.
    Rawlinson N,Reading A M,Kennett B L N.2006.Lithosphericstructure of Tasmania from a novel form ofteleseismictomography[J].Journal of Geophysical Research:Solid Earth(1978-2012),111(B2).
    Rawlinson N,Sambridge M.2004.Multiple reflection andtransmission phases in complex layered media using a multistagefast marching method[J].Geophysics,69(5):1338-1350.
    Rawlinson N.,ReadingA M,Kennett B L N.2006.Lithosphericstructure of Tasmania from a novel form ofteleseismictomography[J].Journal of Geophysical Research:Solid Earth(1978-2012),111(B2).
    Rawlinson N,Pozgay S,Fishwick S.2010.Seismic tomography:Awindow into deep Earth[J].Physics of the Earth and PlanetaryInteriors,178(3):101-135.
    Sethian J A.1996.A fast marching level set method formonotonically advancing fronts[J].Proceedings of the NationalAcademy of Sciences,93(4):1591-1595.
    Sethian J A,Popovici A M.1999.3-D traveltime computation usingthe fast marching method[J].Geophysics,64(2):516-523.
    Shearer P M.2009.Introduction to Seismology[M].CambridgeUniversity Press.
    Teng J W,Liu F T,Quan Y L,et al.1994.The crust and uppermantle seismic tomography in the orogenic belts andsedimentary basins of the Northwest China[J].Progress in solidGeophysics of China(in Chinese),66-80.
    Teng J W,Zhang Z J,Bai W M.2003.Lithosphere Physics[J].Beijing:Science Press.
    Tian Y,Liu C,Feng X.2011.P-wave velocity structure of crustand upper mantle in Northeast China and its control on theformation of mineral and energy[J].Chinese J.Geophys.(inChinese),54(2):407-414.
    Tian Y,Zhao D P,Liu C,et al.2009.A review of body-wavetomography and its applications to studying the crust and mantlestructure in China[J].Earth Science Frontiers(in Chinese),16(2):347-360.
    Tian Y,Zhao D P,Sun R M,et al.2007.The 1992 Landersearthquake:effect of crustal heterogeneity on earthquakegeneration[J].Chinese J.Geophys.(in Chinese),50(5):1488-1496.
    ThurberC H,Ellsworth W L.1980.Rapid solution of ray tracingproblems in heterogeneous media[J].Bulletin of theSeismological Society of America,70(4):1137-1148.
    Thurber C H.1983.Earthquake locations and three-dimensionalcrustal structure in the Coyote Lake area,central California[J].Journal of Geophysical Research:Solid Earth(1978–2012),88(B10):8226-8236.
    Tromp J,Komattisch D,Liu Q.2008.Spectral-element and adjointmethods in seismology[J].Communications in ComputationalPhysics,3(1):1-32.
    Um J,Thurber C.1987.A fast algorithm for two-point seismic raytracing[J].Bulletin of the Seismological Society of America,77(3):972-986.
    Van Der Hilst R D,Maarten V.2005.Banana-doughnut kernels andmantle tomography[J].Geophysical Journal International,163(3):956-961.
    Van der Hilst R D,Widiyantoro S,Engdahl E R.1997.Evidencefor deep mantle circulation from global tomography[J].Nature,386(6625):578-584.
    Vidale J.1988.Finite-difference calculation of travel times[J].Bulletin of the Seismological Society of America,78(6):2062-2076.
    Vidale J E.1990.Finite-difference calculation of traveltimes in threedimensions[J].Geophysics,55(5):521-526.
    Virieux J,Farra V.1991.Ray tracing in 3-D complex isotropicmedia:An analysis of the problem.Geophysics,56(12):2057-2069.
    Wang H,Huang D C.2000.Seismic tomography and its applicationto rock mass structure study[J].Journal of engineering geology(in Chinese),8(1):109-117.
    Woodward M J.1992.Wave-equation tomography[J].Geophysics,57(1):15-26.
    Wu Z B,Xu T,Li F,et al.2014.Acoustic full-waveform inversionin the finite frequency-domain[J].See:Chen Y T,eds,thephysics of the earth’s interior and dynamics in China-tocelebrate the 60th anniversary of engaging in the geophysicalresearch of the academician TengJiwen(in Chinese),Beijing:Science press,306-314.
    Xu T,Li F,Wu Z B,et al.2014.A successive three-pointperturbation method for fast ray tracing in complex 2Dand 3Dgeological models[J].Tectonophysics,doI:10.1016/j.tecto.2014.02.012.
    Xu T,Xu G M,Gao E G,et al.2004.Block modeling and shootingray tracing in complex 3-D media[J].Chinese J.Geophys.(inChinese),47(6):1118-1126.
    Xu T,Xu G,Gao E,Li Y,Jiang X,Luo K.2006.Block modelingand segmentally iterative ray tracing in complex 3Dmedia[J].Geophysics,71(3):T41-T51.
    Xu T,Zhang Z,Gao E,Xu G,Sun L.2010.Segmentally iterativeray tracing in complex 2Dand 3Dheterogeneous block models[J].Bull.Seismol.Soc.Am.,100(2):841-850.
    Xu Y,Liu F T,Liu J H,et al.2001.Deep features of continentalcollision belts in northwestern China and their dynamicimplications[J].Chinese J.Geophys.(in Chinese),44(1):41-47.
    Xu Y,Liu J H,Hao T Y,et al.2006.P wave velocity structureand tectonics analysis of lithospheric mantle in eastern Chinaseas and adjacent regions[J].Chinese J.Geophys.(inChinese),49(4):1053-1061.
    Yang G H.2009.Travel time tomography method with Fresnelvolume and its applications[D].Xiamen University.
    Yang T,Shen Y.2006.Frequency-dependent crustal correction forfinite-frequency seismic tomography[J].Bulletin of theSeismological Society of America,96(6):2441-2448.
    Zelt C A,Smith R B.1992.Seismic traveltime inversion for 2-Dcrustal velocity structure[J].Geophysical journal international,108(1):16-34.
    Zhang F X.2011.Finite-frequency traveltime tomography by bodywave and its application in the North China[D].Institute ofGeophysics,China Earthquake Administration.
    Zhang F X,Li Y H,Wu Q J,et al.2011.The P wave velocitystructure of upper mantle beneath the North China andsurrounding regions from FMTT[J].Chinese J.Geophys.(inChinese),54(5):1233-1242.
    Zhang F X,Wu Q J,Li Y H.2013.The traveltime tomographystudy by teleseismic P wave data in the Northeast China area[J].Geophys.(in Chinese),56(8):2690-2700.
    Zhao D,Christensen D,Pulpan H.1995.Tomographic imaging ofthe Alaska subduction zone[J].Journal of GeophysicalResearch:Solid Earth(1978–2012),100(B4):6487-6504.
    Zhao D,Hasegawa A,Horiuchi S.1992.Tomographic imaging of Pand S wave velocity structure beneath northeastern Japan[J].Journal of Geophysical Research:Solid Earth(1978–2012),97(B13):19909-19928.
    Zhao D,Hasegawa A,Kanamori H.1994.Deep structure of Japansubduction zone as derived from local,regional,andteleseismicevents[J].Journal of Geophysical Research:SolidEarth(1978–2012),99(B11):22313-22329.
    Zhao D,Kanamori H,Humphreys E.1996.Simultaneous inversionof local and teleseismic data for the crust and mantle structure ofsouthern California[J].Physics of the earth and planetaryinteriors,93(3):191-214.
    Zhao D,Wang Z,Umino N,et al.2007.Tomographic imagingoutside a seismic network:Application to the northeast Japanarc[J].Bulletin of the Seismological Society of America,97(4):1121-1132.
    Zhao L,Allen R M,Zheng T,et al.2009.Reactivation of anArchean craton:Constraints from P-and S-wave tomography inNorth China[J].Geophysical Research Letters,36:L17306.
    Zhao L,Allen R M,Zheng T,et al.2012.High-resolution bodywave tomography models of the upper mantle beneath easternChina and the adjacent areas[J].Geochemistry,Geophysics,Geosystems,13(6).
    Zhao L,Jordan T H,Chapman C H.2000.Three-dimensionalFréchet differential kernels for seismicdelaytimes[J].Geophysical Journal International,141(3):558-576.
    Zheng H W,Li T D,Gao R,et al.2007.Teleseismic P-wavetomography evidence for the Indian lithospheric mantlesubducting northward beneath the Qiangtangterrane[J].Chinese J.Geophys.(in Chinese),50(5):1418-1426.
    Zheng H W,Li T D,Gao R.2011.The progress of seismictomography research in the Tibetan plateau[J].Geophysicaland geochemical exploration,35(2):160-164.
    陈立华,宋仲和.1990.华北地区地壳上地幔P波速度结构[J].地球物理学报,33(5):540-546.
    丁志峰,何正勤,孙为国,等.1999.青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J].地球物理学报,42(2):197-205.
    段永红,张先康,方盛明.2002.华北地区上部地壳结构的三维有限差分层析成像[J].地球物理学报,45(3):362-369.
    郭飚,刘启元,陈九辉,等.2009.川西龙门山及邻区地壳上地幔远震P波层析成像.地球物理学报,52(2):346-355.
    江国明,赵大鹏,张贵宾.2008.日本列岛下太平洋俯冲板块的精细结构[J].地学前缘,15(1):222-231.
    江国明,赵大鹏,张贵宾.2009.远震层析成像中的地壳校正研究及应用[J].地球物理学报,52(6):1508-1514.
    江燕,陈晓非.2011.有限频与射线层析成像比较研究综述[J].地球物理学进展,26(5):1566-1575.
    金安蜀,刘福田,孙永智.1980.北京地区地壳和上地幔的三维P波速度结构[J].地球物理学报,23(2):172-182.
    雷建设,周蕙兰,赵大鹏.2002.帕米尔及邻区地壳上地幔P波三维速度结构的研究[J].地球物理学报,45(6):802-811.
    雷建设,赵大鹏,苏金蓉,等.2009.龙门山断裂带地壳精细结构与汶川地震发震机理[J].地球物理学报,52(2):339-345.
    李飞,徐涛,武振波,等.2013.三维非均匀地质模型中的逐段迭代射线追踪[J].地球物理学报,56(10):3514-3522.
    李志伟,胥颐,郝天珧,等.2006.环渤海地区的地震层析成像与地壳上地幔结构[J].地球物理学报,49(3):797-804.
    吕庆田,姜枚,马开义.1996.三维走时反演与青藏高原南部深部构造[J].地震学报,18(4):451-459.
    滕吉文,刘福田,全幼黎.1994.中国西北地区造山带及沉积盆地的地壳和上地幔地震层析成像[J].见:陈运泰等主编,中国固体地球物理进展,66-80.
    滕吉文,张中杰,白武明.2003.岩石圈物理学[M].北京:科学出版社.
    田有,刘财,冯晅.2011.中国东北地区地壳、上地幔速度结构及其对矿产能源形成的控制作用[J].地球物理学报,54(2):407-414.
    田有,赵大鹏,刘财,等.2009.体波走时层析成像方法及其在中国壳幔结构研究中的应用评述[J].地学前缘,16(2):347-360.
    田有,赵大鹏,孙若昧,等.2007.1992年美国加州兰德斯地震-地壳构造不均匀性对地震发生的影响[J].地球物理学报,50(5):1488-1496.
    王辉,黄鼎成.2000.地震层析成像方法及其在岩体结构研究中的应用[J].工程地质学报,8(1):109-117.
    武振波,徐涛,李飞,等.2014.声波方程的有限频率域全波形反演[J].见:陈运泰等主编,中国大陆地球内部物理学与动力学研究-庆贺滕吉文院士从事地球物理学研究60周年,北京:科学出版社[M],306-314.
    徐涛,徐果明,高尔根,等.2004.三维复杂介质的块状建模和试射射线追踪[J].地球物理学报,47(6):1118-1126.
    胥颐,刘福田,刘建华,等.2001.中国西北大陆碰撞带的深部构造特征和地球动力学意义[J].地球物理学报,44(1):40-47.
    胥颐,刘建华,郝天珧,等.2006.中国东部海域及邻区岩石层地幔的P波速度结构与构造分析[J].地球物理学报,49(4):1053-1061.
    杨国辉.2009.菲涅耳体旅行时层析成像方法及应用研究[D].厦门大学.
    张风雪.2011.有限频体波走时层析成像及其在华北地区的应用[D].中国地震局地球物理研究所.
    张风雪,李永华,吴庆举,等.2011.FMTT方法研究华北及邻区上地幔P波速度结构[J].地球物理学报,54(5):1233-1242.
    张风雪,吴庆举,李永华.2013.中国东北地区远震P波走时层析成像研究[J].地球物理学报,56(8):2690-2700.
    郑洪伟,李廷栋,高锐,等.2007.印度板块岩石圈地幔向北俯冲到羌塘地体之下的远震P波层析成像证据[J].地球物理学报,50(5):1418-1426.
    郑洪伟,李廷栋,高锐.2011.青藏高原地震层析成像的研究进展[J].物探与化探,35(2):160-164.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700