用户名: 密码: 验证码:
近地表地震层析成像方法综述
详细信息    查看全文 | 推荐本文 |
  • 英文篇名:Review of seismic tomography methods in near-surface structures reconstruction
  • 作者:张明辉 ; 刘有山 ; 侯爵 ; 徐涛 ; 白志明
  • 英文作者:ZHANG Ming-hui;LIU You-shan;HOU Jue;XU Tao;BAI Zhi-ming;Institute of Information Engineering, Binzhou University;State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Science;University of Chinese Academy of Sciences;Institute of Geophysics,China Earthquake Administration;CAS Center for Excellence in Tibetan Plateau Earth Science;
  • 关键词:近地表 ; 地震层析成像 ; 高频近似 ; 有限频 ; 波动方程
  • 英文关键词:Near-surface;;Seismic tomography;;High frequency approximation;;Finite frequency;;Wave equation
  • 中文刊名:DQWJ
  • 英文刊名:Progress in Geophysics
  • 机构:滨州学院信息工程学院;中国科学院地质与地球物理研究所岩石圈演化国家重点实验室;中国科学院大学;中国地震局地球物理研究所;中国科学院青藏高原地球科学卓越创新中心;
  • 出版日期:2019-02-15
  • 出版单位:地球物理学进展
  • 年:2019
  • 期:v.34;No.153
  • 基金:国家重点研发计划“深地资源勘查开采”重点专项(2016YFC0600201,2016YFC0600101,2016YFC0600302);; 国家自然科学基金(41574082,41874065,41604076,41804060,41774097)联合资助
  • 语种:中文;
  • 页:DQWJ201901008
  • 页数:16
  • CN:01
  • ISSN:11-2982/P
  • 分类号:54-69
摘要
近地表地层与人类生产生活密切相关,利用地震层析成像方法准确重建浅部地壳速度结构有助于开展高精度地震勘探、探查浅部矿产资源、规避潜在自然灾害,并利于城市地下空间建设.中国大陆地表条件复杂,尤其中西部盆岭结合带地形起伏剧烈,对浅部地壳精确速度建模构成严重挑战.本文系统论述了地震层析成像领域基于高频近似理论的走时成像方法和有限频层析成像方法,阐明两类方法的基本原理、存在问题和发展方向等.依据正演走时有无显式射线追踪,基于高频近似理论的走时成像方法分为传统走时层析成像方法和无射线路径的走时层析成像方法.基于射线追踪的传统走时层析成像方法,在浅层速度强烈变化时,因存在阴影区或多路径现象引起成像失真,严重影响成像效率;而无射线路径的层析成像方法通过程函方程走时场的正传和逆传直接计算敏感核,并利用伴随状态法获得目标函数的梯度,具有快速、稳健的优点.以上两种基于地震射线高频近似理论的走时成像方法由于未考虑地震波频率的带限性,存在波散射、波前愈合及反演约束差等问题.有限频层析成像方法克服了射线理论"无限高频"假设所带来的弊端,已成为重要的研究方向之一.该类方法主要分为射线有限频层析成像方法和基于波动方程的有限频层析成像方法.射线有限频层析成像方法能够提高成像的分辨率,但在方法本质上仍依赖于射线理论,较难处理较复杂的波现象问题;基于波动方程的有限频层析成像方法能准确处理复杂地质问题,提高成像可靠性并能以图像形式直观展示地球内部地震波的速度结构分布,但是该方法在实际应用中强烈依赖于数据中的低频信息及较精确的初始速度模型,其推广应用仍需进一步探索.
        The near-surface crustal structure is closely related to human production and life. Accurate reconstruction of the velocity structure of shallow crust by seismic tomography can help us to carry out high-precision seismic exploration,explore shallow resources and minerals,avoid potential natural disasters,and construct urban underground space. Complex surface conditions are widely found in China mainland, especially in the basin-range junction zone of central and western China, which poses a serious challenge to accurately reconstruct the velocity model of shallow crust. This paper discusses systematically the basic principle,existing problems and development direction of two types of current seismic tomography methods,namely,traveltime tomography method based on high frequency approximation theory and on the finite frequency tomography theory. The former is based on the high-frequency approximation theory of seismic rays. According to whether there is explicit ray tracing in the forward calculation,it can be divided into the conventional traveltime tomography method and the traveltime tomography method without ray paths. The conventional traveltime tomography method utilizes the traditional ray tracing methods.When the near-surface velocity changes strongly, the imaging distortion is usually caused by the existence of shadow zone or multipaths,which seriously affects the imaging efficiency. While,the traveltime tomography method without ray paths calculates the sensitive kernel directly by forward and backward propagation of the traveltime field of the eikonal equation,and obtains the gradient of the objective function by the adjoint state method,which has the advantages of fast and robustness. Since the band-limited property of the seismic wave frequency is not considered,there are some problems of the traveltime tomography based on high frequency approximation theory,such as wave scattering,wave front healing and poor inversion constraint. The latter, that is, the finite frequency tomography overcomes the drawbacks of the assumption of"infinite high frequency" in ray theory and has become more important. The methods are mainly divided into ray finite frequency tomography method and the finite frequency tomography methodbased on wave equation. Ray finite frequency tomography can improve the resolution of imaging,but in essence,it still depends on ray theory and is difficult to deal with complex wave phenomena.The finite frequency tomography method based on wave equation can accurately deal with complex geological problems,improve imaging reliability and visually display the velocity structure of seismic waves in the earth's interior in the form of images. However,this method strongly depends on sufficient low-frequency information in data or a more accurate initial velocity model in practical applications,and its wide application requires further exploration and development.
引文
Aki K,Christofferson A,Husebye E S. 1977. Determination of the three-dimensional seismic structure of the lithosphere [J]. J. Geophys. Res.,82(2): 277-296.
    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.
    Bai C Y,Greenhalgh S. 2006. 3-D local earthquake hypocenter determination with an ‘irregular’ shortest-path method [J]. Bull. Seism. Soc. Am.,96(6): 2257-2268,doi: 10.1785/0120040178.
    Bai C Y,Greenhalgh S,Zhou B. 2007. 3D ray tracing using a modified shortest-path method [J]. Geophysics,72(4): T27-T36,doi: 10.1190/1.2732549.
    Bai C Y,Huang G J,Li Z S. 2011. Simultaneous inversion combining multiple-phase travel times within 3D complex layered media [J]. Chinese Journal of Geophysics (in Chinese),54(1): 182-192,doi: 10.3969/j.issn.0001-5733.2011.01.019.白超英,黄国娇,李忠生. 2011. 三维复杂层状介质中多震相走时联合反演成像[J]. 地球物理学报,54(1): 182-192,doi: 10.3969/j.issn.0001-5733.2011.01.019.
    Bai C Y,Huang G J,Zhao R. 2010. 2D/3D irregular shortest-path ray tracing for multiple arrivals and its applications [J]. Geophys. J. Int.,183(3): 1596-1612,doi: 10.1111/j.1365-246x.2010.04817.x.
    Bona A,Slawinski M A,Smith P. 2009. Ray tracing by simulated annealing: Bending method [J]. Geophysics,74(2): T25-T32,doi: 10.1190/1.3068006.
    Brenders A J,Pratt R G. 2007. Full waveform tomography for lithospheric imaging: Results from a blind test in a realistic crustal model [J]. Geophys. J. Int.,168(1): 133-151,doi: 10.1111/j.1365-246X.2006.03156.x.
    Brossier R,Operto S,Virieux J. 2009. Seismic imaging of complex onshore structures by 2D elastic frequency-domain full-waveform inversion [J]. Geophysics,74(6): WCC105-WCC118,doi: 10.1190/1.3215771.
    Brzostowski M A,Mcmechan G A. 1992. 3-D tomographic imaging of near-surface seismic velocity and attenuation [J]. Geophysics,57(3): 396- 403.
    Bunks C,Saleck F M ,Zaleski S,et al. 1995. Multiscale seismic waveform inversion [J]. Geophysics,60(5): 1457-1473.
    Capdeville Y. 2005. An efficient Born normal mode method to compute sensitivity kernels and synthetic seismograms in the Earth [J]. Geophys. J. Int.,163(2): 639- 646,doi: 10.1111/j.1365-246X.2005.02765.x.
    Cerveny V. 2001. Seismic ray theory [M]. Cambridge University Press.
    Cerveny V,Soares J E P. 1992. Fresnel volume ray tracing [J]. Geophysics,57(7): 902-915.
    Chavent G. 1974. Identification of functional parameters in partial differential equations [J]. Joint Automatic Control Conference,12: 155-156.
    Christensen N I,Mooney W D. 1995. Seismic velocity structure and composition of the continental crust: A global view [J]. J. Geophys. Res.,100(B6): 9761-9788.
    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,doi: 10.1046/j.1365-246X.2000.00070.x.
    Dahlen F A,Tromp J. 1998. Theoretical global seismology [M]. Princeton University Press.
    Ding J C,Sun W B,Huang X G,et al. 2017. The strategies of FWI realization for marine seismic data [J]. Geophysical Prospecting for Petroleum (in Chinese),56(1): 75-80,doi: 10.3969/j.issn.1000-1441.2017.01.009.丁继才,孙文博,黄小刚,等. 2017. 海上地震数据全波形反演实际应用[J]. 石油物探,56(1): 78-80,doi: 10.3969/j.issn.1000-1441.2017.01.009.
    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]. Chinese Journal of Geophysics (in Chinese),42(2): 197-205,doi: 10.3321/j.issn:0001-5733.1999.02.007.丁志峰,何正勤,孙为国,等. 1999. 青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J]. 地球物理学报,42(2): 197-205,doi: 10.3321/j.issn:0001-5733.1999.02.007.
    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 Journal of Geophysics (in Chinese),45(3): 362-369,doi: 10.3321/j.issn:0001-5733.2002.03.008.段永红,张先康,方盛明. 2002. 华北地区上部地壳结构的三维有限差分层析成像[J]. 地球物理学报,45(3): 362-369,doi: 10.3321/j.issn:0001-5733.2002.03.008.
    Fichtner A,Bunge H P,Igel H. 2006a. The adjoint method in seismology-I. Theory [J]. Physics of the Earth and Planetary Interiors,157(1-2): 86-104,doi: 10.1016/j.pepi.2006.03.016.
    Fichtner A,Bunge H P,Igel H. 2006b. The adjoint method in seismology-II. Applications: traveltimes and sensitivity functionals [J]. Physics of the Earth and Planetary Interiors,157(1-2): 105-123,doi: 10.1016/j.pepi.2006.03.018.
    Fichtner A,Kennett B L N,Igel H,et al. 2009. Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint Methods [J]. Geophys. J. Int.,179(3): 1703-1725,doi: 10.1111/j.1365-246X.2009.04368.x.
    Fichtner A,Tramper J. 2011. Resolution analysis in full waveform inversion [J]. Geophys. J. Int.,187(3): 1604-1624,doi: 10.1111/j.1365-246X.2011.05218.x.
    Fischer R,Lees J M. 1993. Shortest path ray tracing with sparse graphs [J]. Geophysics,58(7): 987-996.
    Gao E G,Xu G M. 1996. A new kind of step by step iterative ray-tracing method [J]. Chinese Journal of Geophysics (in Chinese),39(S1): 302-308.高尔根,徐果明. 1996. 二维速度随机分布逐步迭代射线追踪方法[J]. 地球物理学报,39(S1): 302-308.
    Gao E G,Xu G M,Jiang X Y,et al. 2002. Iterative ray-tracing method segment by segment under 3-D construction [J]. Oil Geophysical Prospecting (in Chinese),37(1): 11-16,doi: 10.3321/j.issn:1000-7210.2002.01.003.高尔根,徐果明,蒋先艺,等. 2002. 三维结构下逐段迭代射线追踪方法[J]. 石油地球物理勘探,37(1): 11-16,doi: 10.3321/j.issn:1000-7210.2002.01.003.
    Gao E G,Xu G M,Zhao Y. 1998. Segmentally-iterative ray tracing method for any interface [J]. Oil Geophysical Prospecting (in Chinese),33(1): 54- 60.高尔根,徐果明,赵焱. 1998. 一种任意界面的逐段迭代射线追踪方法[J]. 石油地球物理勘探,33(1): 54- 60.
    Gauthier O,Virieux J,Tarantola A. 1986. Two-dimensional nonlinear inversion of seismic waveforms: numerical results [J]. Geophysics,51(7): 1387-1403.
    Gautier S,Nolet G,Virieux J. 2008. Finite-frequency tomography in a crustal environment: Application to the western part of the Gulf of Corinth [J]. Geophysical Prospecting,56(4): 493-503,doi: 10.1111/j.1365-2478.2007.00683.x.
    Hole J A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography [J]. Journal of Geophysical Research,97(B5): 6553- 6562.
    Hole J A,Clowes R M. Ellis R M. 1993. Interpretation of three-dimensional seismic refraction data from western Hecate Strait,British Columbia: structure of the crust [J]. Canadian Journal of earth sciences,30(7): 1440-1452.
    Hole J A,Zelt B C. 1995. 3-D finite-difference reflection travel times [J]. Geophys. J. Int.,121(2): 427- 434.
    Hou H S,Gao R,Lu Z W,et al. 2010. Reflection seismic first-arrival wave tomography of Longqiao iron deposit and concealed deposit forecast in Luzong iron-polymetallic ore concentrated area [J]. Acta Petrologica Sinica (in Chinese),26(9): 2623-2629.侯贺晟,高锐,卢占武,等. 2010. 庐枞铁多金属矿集区龙桥铁矿反射地震初至波层析成像与隐伏矿床预测[J]. 岩石学报,26(9): 2623-2629,doi: 1000- 0569/2010/026(09)-2323-29.
    Huang G J,Bai C Y. 2010. Simultaneous inversion with multiple traveltimes within 2-D complex layered media [J]. Chinese Journal of Geophysics (in Chinese),53(12): 2972-2981,doi: 10.3969/j.issn.0001-5733.2010.12.021.黄国娇,白超英. 2010. 二维复杂层状介质中地震多波走时联合反演成像[J]. 地球物理学报,53(12): 2972-2981,doi: 10.3969/j.issn.0001-5733.2010.12.021.
    Huang G J,Bai C Y,Qian W. 2016. Multi-phase traveltime tomography in general anisotropic TI media: an example of crosshole tomography [J]. Oil Geophysical Prospecting (in Chinese),51(1): 115-126,doi: 10.13810/j.cnki.issn.1000-7210.2016.01.015.黄国娇,白超英,钱卫. 2016. 一般TI介质中多震相走时层析成像-以井间成像为例[J]. 石油地球物理勘探,51(1): 115-126,doi: 10.13810/j.cnki.issn.1000-7210.2016.01.015.
    Huang G J,Sun J B,Bai C Y,et al. 2018. Seismic multi-wave traveltime tomography in 3D TI media [J]. Oil Geophysical Prospecting (in Chinese),53(1): 63-72,doi: 10.13810/j.cnki.issn.1000-7210.2018.01.008.黄国娇,孙江兵,白超英,等. 2018. 三维TI介质中多波走时层析成像[J]. 石油地球物理勘探,53(1): 63-72,doi: 10.13810/j.cnki.issn.1000-7210.2018.01.008.
    Huang J W,Bellefleur G. 2012. Joint transmission and reflection traveltime tomography using the fast sweeping method and the adjoint-state technique [J]. Geophysical Journal International,188(2): 570-582,doi: 10.1111/j.1365-246x.2011.05273.x.
    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.
    Julian B R,Gubbins D. 1977. Three dimensional seismic ray tracing [J]. Journal of Geophysics,43: 95-113.
    Komatitsch D,Tromp J. 2002a. Spectral-element simulations of global seismic wave propagation-I. Validation [J]. Geophysical Journal International,149(2): 390- 412,doi: 10.1046/j.1365-246X.2002.01653.x.
    Komatitsch D,Tromp J. 2002b. Spectral-element simulations of global seismic wave propagation-II. Three-dimensional models,oceans,rotation and self-gravitation [J]. Geophysical Journal International,150(1): 303-318,doi: 10.1046/j.1365-246X.2002.01716.x.
    Lailly P. 1984. The seismic inverse problem as a sequence of before-stack migrations [C]. Conference on inverse scattering: Theory and application: SIAM,206-220.
    Lan H Q,Chen J Y,Zhang Z J. 2014. A fast sweeping scheme for calculating p wave first-arrival travel times in transversely isotropic media with an irregular surface [J]. Pure and Applied Geophysics,171(9): 2199-2208,doi: 10.1007/s00024- 014- 0836-5.
    Lan H Q,Zhang Z J. 2013a. Topography-dependent eikonal equation and its solver for calculating first-arrival traveltimes with an irregular surface [J]. Geophysical Journal International,193(2): 1010-1026,doi: 10.1093/gji/ggt036.
    Lan H Q,Zhang Z J. 2013b. 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,doi: 10.1785/0120120199.
    Lan H Q,Zhang Z,Xu T,et al. 2012. A comparative study on the fast marching and fast sweeping methods in the calculation of first-arrival traveltime field [J]. Progress in Geophysics (in Chinese),27(5): 1863-1870,doi: 10.6038/j.issn.1004-2903.2012.05.005.兰海强,张智,徐涛,等. 2012. 地震波走时场模拟的快速推进法和快速扫描法比较研究[J]. 地球物理学进展,27(5): 1863-1870,doi: 10.6038/j.issn.1004-2903.2012.05.005.
    Langan R T,Lerche I,Cutler R T. 1985. Tracing of rays through heterogeneous media: An accurate and efficient procedure [J]. Geophysics,50(9): 1456-1465.
    Leung S Y,Qian J L. 2006. An adjoint state method for three-dimensional transmission traveltime tomography using first-arrivals [J]. Communications in Mathematical Sciences,4(1): 249-266,doi: 10.4310/CMS.2006.v4.n1.a10.
    Li F,Xu T,Wu Z B,et al. 2013. Segmentally iterative ray tracing in 3-D heterogeneous geological models [J]. Chinese Journal of Geophysics (in Chinese),56(10): 3514-3522,doi: 10.6038/cjg20131026.李飞,徐涛,武振波,等. 2013. 三维非均匀地质模型中的逐段迭代射线追踪[J]. 地球物理学报,56(10): 3514-3522,doi: 10.6038/cjg20131026.
    Li F,Xu T,Zhang M H,et al. 2014. Seismic traveltime inversion of 3D velocity model with triangulated interfaces [J]. Earthquake Science,27(2): 127-136,doi: 10.1007/s11589- 013- 0025- 0.
    Li J G,Zhao B,Sun S W,et al. 2012. Application of seismic first arrival tomography inversion on the exploration of polymetallic mine [J]. Geoscience,26(6): 1218-1224.李建国,赵斌,孙少伟,等. 2012. 地震初至波速度层析反演在多金属矿探测中的应用[J]. 现代地质,26(6): 1218-1224,doi: 10.3969/j.issn.1000-8527.2012.06.013.
    Lions J L. 1971. Optimal control of systems governed by partial differential equations [M]. Springer Verlag,New York.
    Liu H,Meng F L,Li Y M. 1995. The interface grid method for seeking global minimum travel-time and the correspondent raypath [J]. Chinese Journal of Geophysics (in Chinese),38(6): 823-832.刘洪,孟凡林,李幼铭. 1995. 计算最小走时和射线路径的界面网全局方法[J]. 地球物理学报,38(6): 823-832.
    Liu Q,Gu Y J. 2012. Seismic imaging: From classical to adjoint tomography [J]. Tectonophysics,566-567: 31- 66,doi: 10.1016/j.tecto.2012.07.006.
    Liu Y S,Teng J W,Xu T,et al. 2017a. Robust time-domain full waveform inversion with normalized zero-lag cross-correlation objective function [J]. Geophys. J. Int.,209(1): 106-122,doi: 10.1093/gji/ggw485.
    Liu Y S,Teng J W,Xu T,et al. 2017b. Effects of conjugate gradient methods and step-length formulas on the multiscale full waveform inversion in time domain: numerical experiments [J]. Pure and Applied Geophysics,174(5): 1983-2006,doi: 10.1007/s00024- 017-1512-3.张云姝,顾汉明,师学明. 2005. 基于Moser曲射线追踪的SIRT声波层析成像[J]. 工程地球物理学报,2(3): 167-176,doi: 10.3969/j.issn.1672-7940.2005.03.001.
    Liu Z K,Wu Y G,Liu Y J,et al. 1995. Near-surface static correction method using first-arrival traveltimes velocity tomography [J]. Chinese J. Geophys. (in Chinese),38(Supp.1): 153-159.刘振宽,吴永刚,刘英杰,等. 1995. 利用初至波走时速度成像进行近地表静校正的方法[J]. 地球物理学报,38(增刊I): 153-159.
    Luo Y,Schuster G T. 1991. Wave-equation traveltime inversion[J]. Geophysics,56(5): 645- 653.
    Menke W. 1984. Geophysical data analysis: Discrete inverse theory [M]. Academic Press,Florida.
    Mora P. 1987. Nonlinear two-dimensional elastic inversion of multi-offset seismic data [J]. Geophysics,52(9): 1211-1228.
    Moser T J. 1991. Shortest path calculation of seismic rays [J]. Geophysics,56(1): 59- 67.
    Mulder W A. 2006. The perturbed traveltime equation and the adjoint-states gradient of the traveltime error [J]. Geophysical Journal International,167(2): 679- 683,doi: 10.1111/j.1365-246X.2006.02932.x.
    Nakanishi L,Yamaguchi K. 1986. A numerical experiment on nonlinear image reconstruction from first-arrival time for two-dimensional island arc structure [J]. Journal of Physics of the Earth,34(2): 195-201.
    Noble M,Thierry P,Taillandier C,et al. 2010. High-performance 3D first-arrival traveltime tomography [J]. Leading Edge,29(1): 86-93,doi: 10.1190/1.3284057.
    Nolet G. 2008. A breviary of seismic tomography: Imaging the Interior [M]. Cambridge University Press.
    Operto S,Ravaut C,Improta L,et al. 2004. Quantitative imaging of complex structures from dense wide aperture seismic data by multiscale traveltime and waveform inversion: A case study [J]. Geophysical Prospecting,52(6): 625- 651,doi: 10.1111/j.1365-2478.2004.00452.x.
    Pageot D,Operto S,Vallée M,et al. 2013. A parametric analysis of two-dimensional elastic full waveform inversion of teleseismic data for lithospheric imaging [J]. Geophys. J. Int.,193(3): 1479-1505,doi: 10.1093/gji/ggs132.
    Pereyra V,Lee W H K,Keller H B. 1980. Solving two-point seismic-ray tracing problems in a heterogeneous medium. Part I. A general adaptive finite difference method [J]. Bull. Seismol. Soc. Am.,70(1): 79-99.
    Phillips W S,Fehler M C. 1991. Traveltime tomography: A comparison of popular methods [J]. Geophysics,56(10): 1639-1649.
    Plessix R-é. 2006. A review of the adjoint-state method for computing the gradient of a functional with geophysical applications [J]. Geophysical Journal International,167(2): 495-503,doi: 10.1111/j.1365-246X.2006.02978.x.
    Plessix R-é,Baeten G,Maag J W D,et al. 2012. Full waveform inversion and distance separated simultaneous sweeping: a study with a land seismic data set [J]. Geophysical Prospecting,60(4): 733-747,doi: 10.1111/j.1365-2478.2011.01036.x.
    Plessix R é,Mulder W A,ten Kroode A P E. 2000. Automatic cross-well tomography by semblance and differential semblance optimization: Theory and gradient computation [J]. Geophysical Prospecting,48(5): 913-935,doi: 10.1046/j.1365-2478.2000.00217.x.
    Podvin P,Lecomte I. 1991. Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools [J]. Geophysical Journal International,105(1): 271-284.
    Pollitz F F. 2007. Finite-frequency traveltime tomography of San Francisco Bay regional crustal velocity structure [J]. Geophys. J. Int.,171(2): 630- 656,doi: 10.1111/j.1365-246X.2007.03532.x.
    Pratt R G. 1999. Seismic waveform inversion in the frequency domain,Part I: Theory and verification in a physical scale model [J]. Geophysics,64(3): 888-901.
    Pratt R G,Shin C,Hick G J. 1998. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion [J]. Geophys. J. Int.,133(2): 341-362.
    Pratt R G,Song Z M,Williamson P,et al. 1996. Two-dimensional velocity models from wide angle seismic data by wavefield inversion [J]. Geophys. J. Int.,124(2): 323-340.
    Qian J L,Zhang Y T,Zhao H K. 2007a. A fast sweeping method for static convex Hamilton-Jacobi equations [J]. Journal of Scientific Computing,31(1-2): 237-271,doi: 10.1007/s10915- 006-9124- 6.
    Qian J L,Zhang Y T,Zhao H K. 2007b. Fast sweeping methods for eikonal equations on triangular meshes [J]. SIAM Journal on Numerical Analysis,45(1): 83-107,doi: 10.1137/050627083.
    Rawlinson N,Hauser J,Sambridge M. 2008. Seismic ray tracing and wavefront tracking in laterally heterogeneous media [J]. Advances in Geophysics,49: 203-273,doi: 10.1016/S0065-2687(07)49003-3.
    Rawlinson N,Houseman G A,Collins C D N,et al. 2001. New evidence of Tasmania’s tectonic history from a novel seismic experiment [J]. Geophys. Res. Lett.,28(17): 3337-3340,doi: 10.1029/2001GL013342.
    Rawlinson N,Pozgay S,Fishwick S. 2010. Seismic tomography: a window into deep Earth [J]. Physics of the Earth and Planetary Interiors,178(3- 4): 101-135,doi: 10.1016/j.pepi.2009.10.002.
    Sambridge M S. 1990. Non-linear arrival time inversion: constraining velocity anomalies by seeking smooth models in 3-D [J]. Geophys. J. Int.,102(3): 653- 677.
    Scales J A,Docherty P,Gersztenkorn A. 1990. Regularisation of nonlinear inverse problems: Imaging the near-surface weathering layer [J]. Inverse Problems,6(1): 115-131.
    Sei A,Symes W W. 1994. Gradient calculation of the travel time cost function without ray-tracing [C]. 64th Annual International Meeting,SEG,Expanded Abstracts,1351-1354.
    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.
    Sethian J A. 1999. Level set methods and fast marching methods [M]. Cambridge University Press,Cambridge.
    Sethian J A,Popovici A M. 1999. 3-D traveltime computation using the fast marching method [J]. Geophysics,64(2): 516-523.
    Shen P,Symes W W,Stolk C C. 2003. Differential semblance velocity analysis by wave-equation migration [C]. 73rd Annual International Meeting,SEG,Expanded Abstracts,2132-2135.
    Sirgue L,Barkved O I,Dellinger J,et al. 2010. Full waveform inversion: the next leap forward in imaging at Valhall [J]. First Break,28(4): 65-70,doi: 10.3997/1365-2397.2010012.
    Sun M Y,Zhang J,Zhang W. 2017. Alternating first-arrival traveltime tomography and waveform inversion for near-surface imaging [J]. Geophysics,82(4): R245-R257,doi: 10.1190/geo2016- 0576.1.
    Symes W W. 2008. Migration velocity analysis and waveform inversion [J]. Geophysical Prospecting,56(6): 765-790,doi: 10.1111/j.1365-2478.2008.00698.x.
    Taillandier C,Noble M,Chauris H,et al. 2009. First-arrival traveltime tomography based on the adjoint-state method [J]. Geophysics,74(6): WCB57-WCB66,doi: 10.1190/1.3250266.
    Tang X P,Bai C Y. 2009. Multiple ray tracing within 3-D layered media with the shortest path method [J]. Chinese J Geophys (in Chinese),52(10): 2635-2643,doi: 10.3969/j.issn.0001-5733.2009.10.024.唐小平,白超英. 2009. 最短路径算法下三维层状介质中多次波追踪[J]. 地球物理学报,52(10): 2635-2643,doi: 10.3969/j.issn.0001-5733.2009.10.024.
    Tape C,Liu Q Y,Maggi A,et al. 2010. Seismic tomography of the southern California crust based on spectral-element and adjoint methods [J]. Geophys. J. Int.,180(1): 433- 462,doi: 10.1111/j.1365-246X.2009.04429.x.
    Tape C,Liu Q Y,Tromp J. 2007. Finite-frequency tomography using adjoint methods-Methodology and examples using membrane surface waves [J]. Geophys. J. Int.,168(3): 1105-1129,doi: 10.1111/j.1365-246X.2006.03191.x.
    Tarantola A. 1984. Inversion of seismic reflection data in the acoustic approximation [J]. Geophysics,49(8): 1259-1266.
    Teng J W,Zhang Z J,Bai W M,et al. 2004. Physics of lithosphere [M]. Beijing: Science Press (in Chinese).滕吉文,张中杰,白武明,等. 2004. 岩石圈物理学[M]. 北京: 科学出版社.
    Thurber C 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.
    Tong P,Zhao D,Yang D,et al. 2014a. Wave-equation-based travel-time seismic tomography-Part 1: Method [J]. Solid Earth,5: 1151-1168,doi: 10.5194/se-5-1151-2014.
    Tong P,Zhao D,Yang D,et al. 2014b. Wave-equation-based travel-time seismic tomography-Part 2: Application to the 1992 Landers earthquake (Mw7.3) area [J]. Solid Earth,5(2): 1169-1188,doi: 10.5194/se-5-1169-2014.
    Tong P,Zhao D P,Yang D H. 2011. Tomography of the 1995 Kobe earthquake area: comparison of finite-frequency and ray approaches [J]. Geophys. J. Int.,187(1): 278-302,doi: 10.1111/j.1365-246x.2011.05139.x.
    Tromp J,Komatitsch D,Hj?rleifsdóttir V,et al. 2010. Near real-time simulations of global CMT earthquakes [J]. Geophysical Journal International,183(1): 381-389,doi: 10.1111/j.1365-246X.2010.04734.x.
    Tromp J,Komatitsch D,Liu Q Y. 2008. Spectral-element and adjoint methods in seismology [J]. Communications in Computational Physics,3(1): 1-32.
    Tromp J,Tape C,Liu Q Y. 2005. Seismic tomography,adjoint methods,time reversal and banana-doughnut kernels [J]. Geophys. J. Int.,160(1): 195-216,doi: 10.1111/j.1365-246X.2004.02453.x.
    Um J,Thurber C. 1987. A fast algorithm for two-point ray tracing [J]. Bulletin of the Seismological Society of America,77(3): 972-986.
    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 three dimensions [J]. Geophysics,55(5): 521-526.
    Vinje V,Iversen E,?steb?l K,et al. 1996a. Estimation of multivalued arrivals in 3D models using wavefront construction-Part I [J]. Geophysical Prospecting,44(5): 819-842.
    Vinje V,Iversen E,?steb?l K,et al. 1996b. Estimation of multivalued arrivals in 3D models using wavefront construction-Part II [J]. Geophysical Prospecting,44(5): 843-858.
    Vinje V,Iversen E,Gjoystdal H. 1993. Traveltime and amplitude estimation using wavefront construction [J]. Geophysics,58(8): 1157-1166.
    Virieux J,Farra V. 1991. Ray tracing in 3-D complex isotropic media: An analysis of the problem [J]. Geophysics,56(12): 2057-2069.
    Waheed U,Flagg G,Yarman C. 2016. First-arrival traveltime tomography for anisotropic media using the adjoint-state method [J]. Geophysics,81(4): R147-R155,doi: 10.1190/segam2016-13822039.1.
    Wang C Y,Zhang X K,Ding Z F,et al. 1997. Finite-difference tomography of upper crustal structure in Dabieshan orogenic belt [J]. Chinese Journal of Geophysics (in Chinese),40(4): 459-501.王椿镛,张先康,丁志峰,等. 1997. 大别造山带上部地壳结构的有限差分层析成像[J]. 地球物理学报,40(4): 459-501.
    Wang Y H,Rao Y. 2009. Reflection seismic waveform tomography [J]. Journal of Geophysical Research: Atmospheres,114(B3): 2008JB005916,doi: 10.1029/2008JB005916.
    Warner M,Ratcliffe A,Nangoo T,et al. 2013. Anisotropic 3D full-waveform inversion [J]. Geophysics,78(2): R59-R80,doi: 10.1190/GEO2012- 0338.1.
    White D J. 1989. Two-dimensional seismic refraction tomography [J]. Geophys. J. Int.,97(2): 223-245.
    Williamson P R. 1991. A guide to the limits of resolution imposed by scattering in ray tomography [J]. Geophysics,56(2): 202-207.
    Woodward M J. 1992. Wave-equation tomography [J]. Geophysics,57(1): 15-26.
    Xie C,Liu Y Z,Dong L G,et al. 2014. First arrival traveltime tomography based on the adjoint state method [J]. Oil Geophysical Prospecting (in Chinese),49(5): 877-883.谢春,刘玉柱,董良国,等. 2014. 伴随状态法初至波走时层析[J]. 石油地球物理勘探,49(5): 877-883.
    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,627: 72-81,doi: 10.1016/j.tecto.2014.02.012.
    Xu T,Ning J R,Liu C C,et al. 2007. Influence of the self-organization of the Earth interior upon the traveltime and amplitude of seismic wave [J]. Chinese Journal of Geophysics (in Chinese),50(4): 1174-1181,doi: 10.3321/j.issn:0001-5733.2007.04.026.徐涛,宁俊瑞,刘春成,等. 2007. 地球介质自组织性对地震波走时和振幅的影响[J]. 地球物理学报,50(4): 1174-1181,doi: 10.3321/j.issn:0001-5733.2007.04.026.
    Xu T,Xu G M,Gao E G,et al. 2004. Block modeling and shooting ray tracing in complex 3-D media [J]. Chinese Journal of Geophysics (in Chinese),47(6): 1118-1126,doi: 10.3321/j.issn:0001-5733.2004.06.027.徐涛,徐果明,高尔根,等. 2004. 三维复杂介质的块状建模和试射射线追踪[J]. 地球物理学报,47(6): 1118-1126,doi: 10.3321/j.issn:0001-5733.2004.06.027.
    Xu T,Xu G M,Gao E G,et al. 2006. Block modeling and segmentally iterative ray tracing in complex 3D media [J]. Geophysics,71(3): T41-T51,doi: 10.1190/1.2192948.
    Xu T,Zhang M H,Tian X B,et al. 2014. Upper crustal velocity of Lijiang-Qingzhen profile and its relationship with the seismogenic environment of the Ms6.5 Ludian earthquake [J]. Chinese J. Geophys. (in Chinese),57(9): 3069-3079,doi: 10.6038/cjg20140932.徐涛,张明辉,田小波,等. 2014. 丽江-清镇剖面上地壳速度结构及其与鲁甸Ms6.5级地震孕震环境的关系[J]. 地球物理学报,57(9): 3069-3079,doi: 10.6038/cjg20140932.
    Xu T,Zhang Z J,Gao E G,et al. 2010. Segmentally iterative ray tracing in complex 2D and 3D heterogeneous block models [J]. Bulletin of the Seismological Society of America,100(2): 841-850,doi: 10.1785/0120090155.
    Xu T,Zhang Z J,Liu B F,et al. 2015. Crustal velocity structure in the Emeishan large igneous province and evidence of the Permian mantle plume activity [J]. Science China: Earth Science,58(7): 1133-1147,doi: 10.1007/s11430- 015-5094- 6.徐涛,张忠杰,刘宝峰,等. 2015. 峨眉山大火成岩省地壳速度结构与古地幔柱活动遗迹:来自丽江-清镇宽角地震资料的约束[J]. 中国科学,45(5): 561-576,doi: 10.1007/s11430- 015-5094- 6.
    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,doi: 10.1785/0120060038.
    Zelt B C,Ellis R M,Zelt C A,et al. 2001. Three-dimensional crustal velocity structure beneath the Strait of Georgia,British Columbia [J]. Geophys. J. Int.,144(3): 695-712,doi: 10.1046/j.0956-540X.2000.01364.x.
    Zelt C A,Azara A,Levander A. 2006. 3D seismic refraction traveltime tomography at a groundwater contamination site [J]. Geophysics,71(5): H67-H78.
    Zelt C A,Barton P J. 1998. Three-dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin [J]. Journal of Geophysical Research,103(B4): 7187-7210.
    Zelt C A,Hojka A M,Flueh E R,et al. 1999. 3D simultaneous seismic refraction and reflection tomography of wide-angle data from the central Chilean margin [J]. Geophys. Res. Lett.,26(16): 2577-2580.
    Zelt C A,Smith R B. 1992. Seismic traveltime inversion for 2-D crustal velocity structure [J]. Geophys. J. Int.,108(1): 16-34.
    Zhang J,Toks?z M N. 1998. Nonlinear refraction traveltime tomography [J]. Geophysics,63(5): 1726-1737.
    Zhang Y S,Gu H M,Shi X M. 2005. SIRT tomography inversion of two-dimensional sonic data based on Moser’s curved ray tracing method [J]. Chinese Journal of Engineering Geophysics (in Chinese),2(3): 167-176,doi: 10.3969/j.issn.1672-7940.2005.03.001.
    Zhao D P,Hasegawa A,Horiuchi S. 1992. Tomographic imaging of P and S wave velocity structure beneath northeastern Japan [J]. Journal of Geophysical Research,97(B13): 19909-19928.
    Zhao D P,Hasegawa A,Kanamori H. 1994. Deep structure of Japan subduction zone as derived from local,regional,and teleseismic events [J]. Journal of Geophysical Research,99(B11): 22313-22329.
    Zhao F F,Ma T,Xu T. 2014a. A review of the travel-time calculation methods of seismic first break [J]. Progress in Geophysics (in Chinese),29(3): 1102-1113,doi: 10.6038/pg20140313.赵烽帆,马婷,徐涛. 2014a. 地震波初至走时的计算方法综述[J]. 地球物理学进展,29(3): 1102-1113,doi: 10.6038/pg20140313.
    Zhao F F,Zhang M H,Xu T. 2014b. A review of body wave traveltime tomography methods [J]. Progress in Geophysics (in Chinese),29(3): 1090-1101,doi: 10.6038/pg20140312.赵烽帆,张明辉,徐涛. 2014b. 地震体波走时层析成像方法研究综述[J]. 地球物理学进展,29(3): 1090-1101,doi: 10.6038/pg20140312.
    Zhao H. 2005. A fast sweeping method for eikonal equations [J]. Mathematics of computation,74(250): 603- 628,doi: 10.1090/S0025-5718- 04- 01678-3.
    Zhao L,Jordan T H,Chapman C H. 2000. Three-dimensional Fréchet differential kernels for seismicdelay times [J]. Geophysical Journal International,141(3): 558-576,doi: 10.1046/j.1365-246x.2000.00085.x.
    Zhao R,Bai C Y. 2010. Fast multiple ray tracing within complex layered media: the shortest path method based on irregular grid cells [J]. Acta Seismologica Sinica (in Chinese),32(4): 433- 444,doi: 10.3969/j.issn.0253-3782.2010.04.006.赵瑞,白超英. 2010. 复杂层状模型中多次波快速追踪—一种基于非规则网格的最短路径算法[J]. 地震学报,32(4): 433- 444,doi: 10.3969/j.issn.0253-3782.2010.04.006.
    Zheng T Y,He Y M,Yang J H,et al. 2015. Seismological constraints on the crustal structures generated by continental rejuvenation in northeastern China [J]. Scientific Reports,5: 14995,doi: 10.1038/srep14995.
    Zhou B,Greenhalgh S. 2008. Nonlinear traveltime inversion for 3D seismic tomography in strongly anisotropic media [J]. Geophysical Journal International,172(1): 383-394,doi: 10.1111/j.1365-246X.2007.03649.x.
    Zhou B,Greenhalgh S A. 2005. ‘Shortest path’ ray tracing for most general 2D/3D anisotropic media [J]. Journal of Geophysics & Engineering,2(1): 54- 63,doi: 10.1088/1742-2132/2/1/008.

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

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

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