深地震测深剖面的层析成像研究及其应用
|
摘要
本文首先回顾了深地震测深(DSS)方法的发展。在分析深地震测深剖面层析成像基本理论的基础上提出了利用走时、振幅和重力数据，通过有限差分反演和射线反演方法对地壳上地幔结构进行层析成像的研究思路。
本文第三章研究了相对振幅比的获取方法，探讨了利用振幅和重力资料模拟来约束地壳速度结构的有效性，从而形成本文层析成像工作的总体结构。
本文的第四章通过一系列模型试验详细探讨了利用射线反演方法反演速度和深度参数的实际步骤，认为走时反演对速度层内顶部和底部速度变化的分辨能力很弱，仅仅依靠走时数据难以反演出层内顶部和底部的速度变化，故一般情况下应采用速度层底部和顶部速度相同的纵向均匀模型。在参数反演过程中应采用先反演速度参数、再反演深度参数的方案，这样的步骤使得反演结果较稳定，得到的解更为可信。
利用检测板试验的方法以及与其它二维和三维成像结果的对比验证了本文成像结果的可靠性。检测板试验的结果表明，HQ-13线速度地壳速度结构横向变化的最佳绝对分辨应在30km左右，该测线上界面形态横向变化的绝对分辨可以达到10km以内；取速度节点间距和界面深度节点间距均为20km的检测板试验表明，云南地区四条测线地壳内界面形态的横向分辨良好。
本文第五章利用有限差分反演和射线反演的方法获得了符离集—奉贤地震测深剖面(HQ-13线)精细的地壳上地幔结构剖面，揭示出下扬子地区一些重要的结构和动力学特征：
地壳速度结构在纵向上大致可分上地壳、中地壳和下地壳三部分，横向上可划分为6个块体，各块体的P波速度沿测线方向呈现高速—低速—高速的组合。上部地壳高速异常和低速异常分别与地表的隆起和坳陷高度一致，高速异常和低速异常的接触带往往与地表的断裂位置有良好的对应关系，沉积盖层底部的基底界面横向上起伏变化较大。中地壳的中间层速度分的基本特征类似上地壳底部，横向均匀性显著；中地壳下部速横向不均匀性明显，在灵壁隆起和湖苏隆起的下部有显著的高速异常。下地壳的两构造层速度差异不明显。莫霍界面深度为30～36km。上地幔顶部速度值为8.06～8.29km／s。
成像结果显示符离集大断裂、嵩沟坳陷南界断裂、洪泽-流均沟断裂等均延伸至上地幔，反之郯庐断裂带缺乏深大断裂构造的迹象。就地壳深部结构而言，苏北坳陷区板块的划分办法、各板块的接触边界以及相应的动力学背景，有进一步求证的必要。
根据所得结构剖面，分析了该测线经过区域的地震构造环境。镇江块体容易发生中强地震的解释是：由于地壳结构偏刚性及壳内界面整体上隆，块体中活动性显著的大断裂很可能将地壳深部或上地幔的能量源源不断地传输至中上地壳，在有利部位积聚并诱发地震。苏北坳陷、无锡块体和上海块体的下部存在低速、塑性的中间层。这些块体内部结构之间互动的位能容易被塑性中间层吸收，且由于区域构造活动不显著，从而不容易发生大的地震。
本文第六章对云南地区四条DSS剖面进行层析成像研究，得到以下一些重要认
识：
云南地区地壳大致可分为包括沉积层和基底层的上地壳、包含两个弱反射层的
中地壳以及下地壳构成。各个反射层的界面形态有相似性，各层的速度分布密切联
系，存在显著的深浅构造的一致性，意味着浅部物质活动有深部背景，同时说明云
南地区是典型的构造活动区。云南地区莫霍界面总体形态南端浅、北端深，地壳厚
度约33刀～54刀km。上地慢速度明显偏低，可能与贯穿整个新生代的热过程相联系。
根据界面形态和速度异常分布特征，判别出云南地区断裂构造的深部特征。作
为一级构造单元分界线的红河断裂、怒江断裂是超壳断裂，而昌宁－双江断裂显示为
低角度铲式断层，中下地壳无与之对应的速度异常，可能意味着该断裂切割不深。
澜沧江断裂南段为基底层断裂。中甸断裂、剑q断裂和丽江－宁落断裂组成的断裂
系延伸至下地壳。
通过对龙陵地震区、思茅一普洱地震区、建水附近地震带以及剑川］附近等地震
构造环境的分析，提出一个云南地区强烈地震孕育发生的模式：在历史上火山活动
剧烈的地区，现今深部活动仍未停止，存在较强的应力作用。由于地表火山活动的
停止，来自中下地壳和上地慢的大部分能量得不到释放，从而在有利的构造部位积
聚，进而在流体的作用下，或在其它外动力诱发下，导致孕震系统失稳而发生地震。
本文第七章总结了本文的主要工作和意义，并指出一些需要深入研究的问题
This paper reviews the development of the Deep Seismic Sounding (DSS) method firstly. Based on the analyzing of the DSS's basic theory of tomography,an inversion scheme that can be described as the combination of the finite-difference inversion and ray inversion using traveltime,amplitude and gravity data is proposed.
The third chapter of this paper explored the method of acquiring amplitude-ratio,and also the significance of the amplitude and gravity data for the probing of the crust and upper mantle's structure is investigated. Thus the framework of inversion adopted by this paper is formed.
Several numerical experiments are carried out to detect the actual steps and limitation for the velocity and depth parameter's ray inversion in the fourth chapter,and the we find that:a. the resolution capability of the layer's velocity change along both the top and bottom boundary is very poor. Thus generally the vertical homogeneous layer model whose top velocity equal to that of its bottom boundary should be considered. During the inversion process of velocity and depth parameter,the velocity parameter should be inversed firstly and then do the depth so as to assure the stability of the final results and the reliability of the solution.
Take advantage of the checker-board-test method to test the horizontal resolution of velocity and depth parameter on HQ-13 profile and the four profiles of Yunnan region. It's found that the best absolute lateral resolution of the velocity on the HQ-13 profile is about 30 kilometer,while that of the depth parameter on this profile can be smaller than 10 kilometer. Both of the velocity and depth parameter points are of 20 kilometers interval on the four profiles in Yunnan region,their lateral resolution are also detected by the checker-board-test method. It shows that the resolution of the crust interface's shape on this four profiles is very good,nevertheless that of the velocity structure is relatively poor.
The crust and upper mantle's fine configuration on the HQ-13 profile is obtained by our tomography,which reveals some important structure and dynamic features of this region:
The velocity structure of the crust consists of three layers,the upper crust,the middle crust and the lower crust. However it also can be divided into six lateral blocks,each of them 's velocity is always higher or lower than that of its neighboring blocks. The high or low velocity anomalies's locations are consistent with the uplift and depression respectively,and their contacting boundary always correspond the faults exposed on the ground's surface. The basement interface that is the bottom interface of the sedimentary cover undulates strongly. The characters of velocity's distribution within the mid-layer of the middle crust resemble that of the upper crust's bottom layer,and the velocity strongly disturbed. The prominent higher velocity anomalies exist below the Linbi uplift and the
Hu-Su uplift within this layer. The velocity difference between the two layers of the lower crust is not obvious. The depth of Moho interface is about 30-36 km,and the velocity of the upper mantle is nearly 8.06-8.29 km/s.
The final structure also shows that the Fuliji fault,the Songgou basin south fault,the Hongze-Liujungou fault among with other deep faults stretch into the upper mantle. While there is no distict evidence to prove the Tanlu fault's noted deep and dynamic behavior within the lower crust or even on the upper mantle. Thus there still exists a necessary to research and check the understanding about the boundary's location and the dynamic background.
Base on the velocity profile,the seismotectonic environment along this profile is studied. The reason for the moderate and large earthquake of Zhengjiang block could be attribute as such features:
Because of the structure's more rigidity and the rising of the several interfaces in the crust,the energy from the deep crust of upper mantle is transported to the middle or upper crust persistently. Thus it always accumulates at some locations and this may
本文第三章研究了相对振幅比的获取方法，探讨了利用振幅和重力资料模拟来约束地壳速度结构的有效性，从而形成本文层析成像工作的总体结构。
本文的第四章通过一系列模型试验详细探讨了利用射线反演方法反演速度和深度参数的实际步骤，认为走时反演对速度层内顶部和底部速度变化的分辨能力很弱，仅仅依靠走时数据难以反演出层内顶部和底部的速度变化，故一般情况下应采用速度层底部和顶部速度相同的纵向均匀模型。在参数反演过程中应采用先反演速度参数、再反演深度参数的方案，这样的步骤使得反演结果较稳定，得到的解更为可信。
利用检测板试验的方法以及与其它二维和三维成像结果的对比验证了本文成像结果的可靠性。检测板试验的结果表明，HQ-13线速度地壳速度结构横向变化的最佳绝对分辨应在30km左右，该测线上界面形态横向变化的绝对分辨可以达到10km以内；取速度节点间距和界面深度节点间距均为20km的检测板试验表明，云南地区四条测线地壳内界面形态的横向分辨良好。
本文第五章利用有限差分反演和射线反演的方法获得了符离集—奉贤地震测深剖面(HQ-13线)精细的地壳上地幔结构剖面，揭示出下扬子地区一些重要的结构和动力学特征：
地壳速度结构在纵向上大致可分上地壳、中地壳和下地壳三部分，横向上可划分为6个块体，各块体的P波速度沿测线方向呈现高速—低速—高速的组合。上部地壳高速异常和低速异常分别与地表的隆起和坳陷高度一致，高速异常和低速异常的接触带往往与地表的断裂位置有良好的对应关系，沉积盖层底部的基底界面横向上起伏变化较大。中地壳的中间层速度分的基本特征类似上地壳底部，横向均匀性显著；中地壳下部速横向不均匀性明显，在灵壁隆起和湖苏隆起的下部有显著的高速异常。下地壳的两构造层速度差异不明显。莫霍界面深度为30～36km。上地幔顶部速度值为8.06～8.29km／s。
成像结果显示符离集大断裂、嵩沟坳陷南界断裂、洪泽-流均沟断裂等均延伸至上地幔，反之郯庐断裂带缺乏深大断裂构造的迹象。就地壳深部结构而言，苏北坳陷区板块的划分办法、各板块的接触边界以及相应的动力学背景，有进一步求证的必要。
根据所得结构剖面，分析了该测线经过区域的地震构造环境。镇江块体容易发生中强地震的解释是：由于地壳结构偏刚性及壳内界面整体上隆，块体中活动性显著的大断裂很可能将地壳深部或上地幔的能量源源不断地传输至中上地壳，在有利部位积聚并诱发地震。苏北坳陷、无锡块体和上海块体的下部存在低速、塑性的中间层。这些块体内部结构之间互动的位能容易被塑性中间层吸收，且由于区域构造活动不显著，从而不容易发生大的地震。
本文第六章对云南地区四条DSS剖面进行层析成像研究，得到以下一些重要认
识：
云南地区地壳大致可分为包括沉积层和基底层的上地壳、包含两个弱反射层的
中地壳以及下地壳构成。各个反射层的界面形态有相似性，各层的速度分布密切联
系，存在显著的深浅构造的一致性，意味着浅部物质活动有深部背景，同时说明云
南地区是典型的构造活动区。云南地区莫霍界面总体形态南端浅、北端深，地壳厚
度约33刀～54刀km。上地慢速度明显偏低，可能与贯穿整个新生代的热过程相联系。
根据界面形态和速度异常分布特征，判别出云南地区断裂构造的深部特征。作
为一级构造单元分界线的红河断裂、怒江断裂是超壳断裂，而昌宁－双江断裂显示为
低角度铲式断层，中下地壳无与之对应的速度异常，可能意味着该断裂切割不深。
澜沧江断裂南段为基底层断裂。中甸断裂、剑q断裂和丽江－宁落断裂组成的断裂
系延伸至下地壳。
通过对龙陵地震区、思茅一普洱地震区、建水附近地震带以及剑川］附近等地震
构造环境的分析，提出一个云南地区强烈地震孕育发生的模式：在历史上火山活动
剧烈的地区，现今深部活动仍未停止，存在较强的应力作用。由于地表火山活动的
停止，来自中下地壳和上地慢的大部分能量得不到释放，从而在有利的构造部位积
聚，进而在流体的作用下，或在其它外动力诱发下，导致孕震系统失稳而发生地震。
本文第七章总结了本文的主要工作和意义，并指出一些需要深入研究的问题
This paper reviews the development of the Deep Seismic Sounding (DSS) method firstly. Based on the analyzing of the DSS's basic theory of tomography,an inversion scheme that can be described as the combination of the finite-difference inversion and ray inversion using traveltime,amplitude and gravity data is proposed.
The third chapter of this paper explored the method of acquiring amplitude-ratio,and also the significance of the amplitude and gravity data for the probing of the crust and upper mantle's structure is investigated. Thus the framework of inversion adopted by this paper is formed.
Several numerical experiments are carried out to detect the actual steps and limitation for the velocity and depth parameter's ray inversion in the fourth chapter,and the we find that:a. the resolution capability of the layer's velocity change along both the top and bottom boundary is very poor. Thus generally the vertical homogeneous layer model whose top velocity equal to that of its bottom boundary should be considered. During the inversion process of velocity and depth parameter,the velocity parameter should be inversed firstly and then do the depth so as to assure the stability of the final results and the reliability of the solution.
Take advantage of the checker-board-test method to test the horizontal resolution of velocity and depth parameter on HQ-13 profile and the four profiles of Yunnan region. It's found that the best absolute lateral resolution of the velocity on the HQ-13 profile is about 30 kilometer,while that of the depth parameter on this profile can be smaller than 10 kilometer. Both of the velocity and depth parameter points are of 20 kilometers interval on the four profiles in Yunnan region,their lateral resolution are also detected by the checker-board-test method. It shows that the resolution of the crust interface's shape on this four profiles is very good,nevertheless that of the velocity structure is relatively poor.
The crust and upper mantle's fine configuration on the HQ-13 profile is obtained by our tomography,which reveals some important structure and dynamic features of this region:
The velocity structure of the crust consists of three layers,the upper crust,the middle crust and the lower crust. However it also can be divided into six lateral blocks,each of them 's velocity is always higher or lower than that of its neighboring blocks. The high or low velocity anomalies's locations are consistent with the uplift and depression respectively,and their contacting boundary always correspond the faults exposed on the ground's surface. The basement interface that is the bottom interface of the sedimentary cover undulates strongly. The characters of velocity's distribution within the mid-layer of the middle crust resemble that of the upper crust's bottom layer,and the velocity strongly disturbed. The prominent higher velocity anomalies exist below the Linbi uplift and the
Hu-Su uplift within this layer. The velocity difference between the two layers of the lower crust is not obvious. The depth of Moho interface is about 30-36 km,and the velocity of the upper mantle is nearly 8.06-8.29 km/s.
The final structure also shows that the Fuliji fault,the Songgou basin south fault,the Hongze-Liujungou fault among with other deep faults stretch into the upper mantle. While there is no distict evidence to prove the Tanlu fault's noted deep and dynamic behavior within the lower crust or even on the upper mantle. Thus there still exists a necessary to research and check the understanding about the boundary's location and the dynamic background.
Base on the velocity profile,the seismotectonic environment along this profile is studied. The reason for the moderate and large earthquake of Zhengjiang block could be attribute as such features:
Because of the structure's more rigidity and the rising of the several interfaces in the crust,the energy from the deep crust of upper mantle is transported to the middle or upper crust persistently. Thus it always accumulates at some locations and this may
引文
陈炳蔚，王铠元，刘万熹，等．1987．怒江—澜沧江—金沙江地区大地构造[M]．北京：地质出版社，204
陈步云，高文海．1988．贾家湾—沙园剖面地壳上地幔的深部结构[A]．见：国家地震局科技监测司编．中国大陆深部构造的研究与进展[C]．北京：地质出版社，154～168
陈国英，宋仲和，安昌强，等．1995．中国北部及其邻区地壳上地幔三维速度结构[J]．地球物理学报，38：321～328
陈沪生．1988．下扬子地区 HQ-13 线的综合地球物理调查及其地质意义[J]．石油与天然气地质，9(3)：211～222
陈沪生，周雪清，李道琪，等．1992．中国东部灵壁—奉贤(HQ—13)地学断面图1：1，000，000说明书[M]．北京：地质出版社．
陈培善，刘福田，李强，等．1990．云南地区速度结构的横向不均匀性[J]．中国科学，B辑(4)：431～439
陈天培，秦德厚，韩文华，等．1969．中华人民共和国区域地质调查报告．1：20万，建水幅[M]．昆明：云南省地质局，119
陈希孺．1992．概率论与数理统计[M]．合肥：中国科学技术大学出版社，138．
陈学波，吴跃强、杜平山等．1988．龙门山构造带两侧地壳速度结构特征 [A]．见：国家地震局科技监测司编．中国大陆深部构造的研究与进展[C]．北京：地质出版社，97～113
陈学波，李金森，赵静娴，等．1994．三江构造带基底构造特征与强震活动关系[A]．见：陈运泰，阚荣举，滕吉文主编．中国固体地球物理学进展[C]．北京：海洋出版社，41～54
陈运泰，滕吉文，张中杰，等．2001．地球物理学的回顾与展望[J]．地球科学进展，16(5)：634～642
崔作舟，卢德源，陈纪平，等．1987．攀西地区的深部地壳结构与构造[J]．地球物理学报，30：566～580
丁志峰，何正勤，孙为国，等．1999a．青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J]．地球物理学报，42(2)：197～205
丁志峰．1999b．近震层析成像的理论及应用[D]．40
冯锐．1985．中国地壳厚度及上地幔密度分布(三维重力反演结果)[J]．地震学报，7(2)：143～157
冯锐，郭强绪，江月昌，等．1986．地震测深的观测与回放系统[J]．地震学报，8：57～71
冯锐，郑书真．1987．对华北岩石圈构造结构的综合研究[J]．科学通报，22：1723～1727
冯锐，陶裕录．1990．地震—重力联合反演[A]．见：中国地球物理学会年刊(1990)[C]．北京：地震出版社，44
冯锐，陶裕录．1993．地震—重力联合反演中的非块状一致性模型[J]．地球物理学报，36：463～474
傅淑芳，朱仁益．1998．地球物理反演问题[M]．北京：地震出版社，10．
国家地震局《深部物探成果》编写组．1986．中国地壳上地幔地球物理探测成果[M]．北京：地震出版社．
胡鸿翔，陆涵行，王椿镛，等．1996a．云南地区地壳结构的爆破地震研究[J]．地球物理学报，29(2)：133～144
胡鸿翔，陈学波，张碧秀，等．1986b．中国中原地区随县—安阳剖面深地震测深资料的解释[J]．地震学报，8(1)：38～49
胡鸿翔，高世玉．1993．云南地区地壳浅部基底速度细结构的研究[J]．中国地震，19(4)：356～363
胡鸿翔，李学清．1994．云南地区孟连—思茅—马龙剖面的基底速度系结构[A]．见：陈运泰，阚荣举，滕吉文主编．中国固体地球物理学进展[C]．北京：海洋出版社，100～106
黄怀勇，张碧秀，朱建亚．1989．安阳—泗水DSS剖面二维速度结构[A]．见：汤永安主编，华北平原南部地区深部地球物理综合探测方法与成果[C]．北京：海洋出版社，64～78
江苏省地震局深部构造组．1986．HQ—13测线人工地震测深资料图集．
江苏省重工业局区测队革命委员会．1970．中华人民共和国地质部区域地质测量报告书(1：200，000)，扬州幅，汤山幅，上党幅[M]．
江苏省地质局．1976．中华人民共和国区域地质调查报告(1：200，000)，无锡幅，苏州幅[M]．
江苏省地质局．1979．中华人民共和国区域地质调查报告(1：200，000)，盱眙幅[M]．
金安蜀，刘福田，孙永智．1980．北京地区地壳和上地幔三维P波速度结构[J]．地球物理学报，23：172～182
阚荣举，林中洋．1986．云南地壳上地幔构造的初步研究[J]．中国地震，2(4)：50～61
阚荣举，韩源．1992．云南遮放至马龙地学断面说明书[M]．北京：地震出版社。
阚荣举，赵晋明．1994．腾冲火山地区的壳幔构造[A]．见：陈运泰，阚荣举，滕吉文，等主编．中国固体地球物理学进展[C]．北京：海洋出版社，23～30
赖晓玲，张先康，郑需要．1997．利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态[J]．地震学报，19(5)：506～516
李成猛．1996．地震资料的加权约束反演方法[J]．石油物探，35(4)：29～36
李清河，郭建康，李流玉，等．1989．清丰—汶上DSS剖面的Q值结构[A]．见：深部地球物理综合探测方法与成果[C]．北京：海洋出版社，78～100
李庆忠．1993．走向精确勘探的道路[M]．北京：石油工业出版社，38．
李松林，樊计昌，吴宁远．1995．利用地震层析成像技术确定横向非均匀介质中速度分布和界面位置[J]．地震研究，18：430～444
李松林，吴宁远，宋占隆，等．1997．速度分布和界面位置的联合反演[J]．地震学报，19(4)：383～392
李雄．1991．地震—重力联合反演—不同方式的实现与对比[A]．见：中国地球物理学会年刊(1991)[C]．北京：地震出版社，7．
林中洋，胡鸿翔，张文彬，等．1993．云南地区地壳上地幔速度结构特征的研究[J]．地震学报，15(4)：427～440
刘宝诚，佟伟，等．1986．腾冲火山地热区的微震观测[J]．地球物理学报，29：547～555
刘昌铨，祝治平，李捍东，等．1983．连云港—临沂—泗水测深剖面及临沂8.5级地震深部构造背景[J]．地震，3：172～183
刘福田．1984．震源位置和速度结构的联合反演(Ⅰ)—理论和方法[J]．地球物理学报，27(2)：167～175
刘福田．1985．我国地球物理反演问题研究的回顾[J]．31(增刊)：144～149
刘福田，曲克信，吴华，等．1986．华北地区的地震层面成象[J]．地球物理学报，29(5)：442～449
刘福田，李强，吴华，等．1989a．用于速度图像重建的层析成像法[J]．地球物理学报，32：46～61
刘福田，曲克信，吴华，等。1989b．中国大陆及其邻近地区的地震层析成象[J]．地球物理学报，32(3)：281～293
刘福田．1995．地球物理学中的反演问题[M]．讲义，17
刘建华，刘福田，吴华，等．1989．中国南北带地壳和上地幔的三维速度图象[J]．地球物理学报，32(2)：143～153
刘清林．1993．地震初至波射线路径的追踪[J]．石油物探，32(2)：14～20
刘清林．1995．程函方程差分法层析成像[J]．石油物探，34(4)：14～26
刘清林．1996．初始模型和约束条件对层析反演结果的影响[J]．石油物探，35(2)：35～43
刘瑞丰，陈培善，李强．1993．云南及其邻近地区三维速度图象[J]．地震学报，15(1)：61～67
刘天佑．1992．重磁异常反演理论与方法[M]．中国地质大学出版社，21
楼海，王椿镛．1999．三维连续密度分布的重力计算及应用[J]．地震学报，21(3)：298～304
楼海，王椿镛，皇甫岗，等．2002．云南腾冲火山上部地壳三维地震速度层析成像[J]．地震学报，24(3)：243～251
马杏垣等．中国岩石圈动力学图集．北京：中国地图出版社，51
卢军，潘卫平，曾晓献．1991．人工爆炸地震野外数据处理装置[J]．地球物理学报，34：717～723
穆治国，佟伟．1987．腾冲火山活动的时代和岩浆来源问题[J]．地球物理学报，30：261～270
秦大正，刘昌森，丁颂华，等．1983．长江断裂带东延问题及裂谷特性的讨论[J]．地震地质，5(2)：70～77
全幼黎，滕吉文．1986．混合法计算GB体波理论地震图[J]．地球物理学报，31：165～175
任纪舜，姜春发，张正坤，等．1980．中国大地构造及其演化[M]．北京：科学出版社．
孙洁等．1989．滇西地区地壳上地幔电性结构与地壳构造运动的关系[J]．地震地质，11：35～45
孙武城，刘昌铨，祝治平，等．1984．揭示郯庐大断裂地震深部的某些特征—连云港—临沂—泗水人工地震测深资料的初步解释[A]．见：构造地质论丛[C]．北京：地质出版社，238～251
孙武城，祝治平，张利，等．1988．对华北地壳上地幔的探测与研究[A]．见：国家地震局《深部物探成果》编写组．中国地壳上地幔地球物理探测成果[C]．北京：地震出版社：19～37
孙武城，徐杰，杨主恩，等．1992．上海奉贤至内蒙古阿拉善左旗地学断面说明书[M]．北京：地震出版社．
汤永安．1989．华北平原南部地区深部地球物理综合探测方法与成果[M]．北京：海洋出版社．
滕吉文，熊绍柏，尹周勋，等．1983．喜马拉雅北部地区的地壳结构模型和速度分布特征[J]．地球物理学报，26：525～540
滕吉文，孙克忠，熊绍柏，等．1985．中国东部马鞍山—常熟—启东地带地壳与上地幔结构和速度分布的爆炸地震研究[J]．地球物理学报，28(2)：155～169
滕吉文，刘福田，全幼黎，等，1994a．中国西北造山带与沉积盆地地区的地壳和地幔的层析成像[A]．见：陈运泰等主编．中国固体地球物理学进展[C]．北京：海洋出版社，66～80
滕吉文．1994b．全球地球科学大断面与地球动力学(I)[J]．地球物理学进展，9 (1)：40～47
田宗勇，吴建成，刘家琦，等．1997．地震反射层析成像方法的研究及其在工程中的应用[J]．物探与化探，21(5)：394～397
王椿镛．1982．层状不均匀介质中合成地震图的反射率法[J]．地球物理学报，25：424～433
王椿镛．1985．二维射线合成地震图的计算及平凉—渭南地震测深剖面的解释[J]．地震研究，8(6)：607～616
王椿镛．1986．层状不均匀衰减介质中的反射率法[J]．地震学报，8(增刊)：85～94
王椿镛，陆原，陈光英，等．1989．地震测深数据处理软件包[J]．地震地磁观测和研究，10(5)：50～58
王椿镛，陆原，谢剑波，等．1991．地震测深数据文件系统[J]．地球物理学报，34：455～464
王椿镛，楼海，宋亦青．1994．人工地震测深数据库系统[J]．地震学报，16：89～95
王椿镛．1994c．中国深地震反射剖面探测与研究[A]．见：陈运泰等主编．中国固体地球物理学进展[C]．北京：海洋出版社，120～127
王椿镛，林中洋．1995a．地壳地震学研究进展[A]．见：陈运泰主编．地球与空间科学观测技术进展[C]．北京：地震出版社，272～282
王椿镛，林中洋，陈学波．1995b．青海门源—福建宁德地学断面综合地球物理研[J]．地球物理学报，38(5)：590～598
王椿镛，张先康，丁志峰，等．1997a．大别造山带上部地壳结构的有限差分层析成像[J]．地球物理学报，40(4)：496～501
王椿镛．1997b．中国岩石层结构研究的回顾与展望[J]．地球物理学报，40(增刊)：82～109
王椿镛．1998．地球内部构造和动力学的地震学研究[J]．地学前缘(中国地质大学，北京)，5(1～2)：91～98
王椿镛，楼海，吴建平等．2001a．腾冲火山地热区地壳结构的地震学研究[J]．地震学报，24(3)：
王椿镛，楼海，魏修成，等．2001b．天山北缘的地壳结构和1906年玛纳斯地震的地震构造[J]．地震学报，23(5)：460～470
王椿镛，Mooney W D, 王溪莉，等．2002a．川滇地区地壳上地幔三维速度结构研究[J]．地震学报，24(1)：1～16
王椿镛，楼海，吴建平，等．2002b．腾冲火山地热区地壳结构的地震学研究[J]．地震学报，24(3)：231～242
王华忠，谢海兵，马在田．1997．有限差分法地震波走时计算[J]．同济大学学报，25(3)：318～321
王家映．1998．地球物理反演理论[M]．武汉：中国地质大学出版社，
谢海兵，马在田，王燕琨．1997．一种新程函方程差分法层析成像[J]．石油地球物理勘探，32(5)：623～627
熊绍柏，滕吉文，尹周勋，等．1986．攀西构造带南部地壳上地幔结构和爆破地震研究[J]．地球物理学报，29(3)：235～244
杨文采．1989．地球物理反演和地震层析成像[M]，北京：地质出版社
杨文采，李幼铭，等，1992．应用地震层析成像[M]，北京：地质出版社
杨晓春，李小凡，张美根．2001．地震波反演方法研究的某些进展及其数学基础[J]．地球物理学进展，16(4)：96～109
云南省地质局．1976．中华人民共和国区域地质调查报告(1：200，000)，普洱幅[M]．
曾融生，阚荣举，何传大．1960．柴达木盆地低频地震探测结晶基底的工作方法[J]．科学通报，10：313～316
曾融生，阚荣举．1961a．柴达木盆地西部地壳深界面反射波[J]．地球物理学报，10：120～125
曾融生，阚荣举，何传大．1961b．柴达木盆地地震探测的基岩首波和大角度反射[J]．地球物理学报，10：54～66
曾融生，朱介寿，周兵，等．1992．西藏高原及其东部邻区的三维地震波速度结构与大陆碰撞模型．地震学报，14(增刊)：523～533
张少泉，俞玉琪，高世玉，等．1983．横向非均匀介质的地震射线追踪及其应用[J]．地震地磁观测与研究，4(6)：1～8
张少泉，武利均，郭建明，等．1985．中国西部地区门源—平凉—渭南地震测深剖面资料的分析解释[J]．地球物理学报，28：349～361
张四维，张锁喜，唐荣余，等．1988．下扬子地区符离集—奉贤地震测深资料解释[J]．地球物理学报，31(6)：637～647
张文生，何樵登．1996．有限差分法精确计算地震走时[J]．物探化探计算技术，18 (4)：331～335
张先康，杨卓欣，杨玉春，等．1995．地壳三维结构的层析成像方法—爆炸和地震资料的联合反演[J]．地震学报，17(4)：422～431
张先康，张成科，赵金仁，等．2002．长白山天池火山区岩浆系统深部结构的深地震测深研究[J]．地震学报，24(2)：135～143
张元生，李清河，徐果明．1998．联合利用走时与波形反演技术研究地壳三维速度结构(I)—理论与方法[J]．西北地震学报，20(2)：8～51
张中杰，滕吉文，张霖斌，等．1996．地球深部研究的复杂性及其对策[J]．地球物理学进展，11(3)：11～20
张中杰，滕吉文，李英康．2002．藏南地壳速度结构与地壳物质东西向“逃逸”—以佩枯错—普莫雍错宽角反射剖面为例[J]．中国科学，32(10)：793～798
赵爱华，张中杰，王光杰，等．2000．非均匀介质中地震波走时与射线路径快速计算技术[J]．地震学报，22(2)：151～157
周兵，朱介寿，S. A. Greenhalgh, 等．1994．一种新的地震射线层析成像计算方法[J]．石油物探，33(1)：45～54
朱介寿，严忠琼，曹俊兴，等．1998．勘探地球物理层析成像软件系统及其应用[J]．物探与化探，22(2)：90～98
朱金明，王丽燕．1992．地震波走时的有限差分法计算[J]．地球物理学报．35(1)：86～92
朱良保．1989．计算横向非均匀介质体波理论地震图的 Maslov 渐近法[J]．地球物理学报，32：163～172
朱仁益．1990．用WKBJ法计算衰减介质中的理论地震图．地球物理学报．33(增刊)：212～224．
朱天飞，刘福田．1982．用地震资料测定岩石层的三维速度结构[J]．地球物理学报．25：500～508．
Ammon C J, Vidale J E. 1993.Tomography without rays[J]. Bull Seism Soc Am, 83:509～528
Anderson D I, Dziewonski A M. 1984.Seismic tomography[J], Scientific American., 251(4):58～66
Benz H M, Smith R B, Mooney W D.1990.Crustal structure of the northwestern basin and range province from the 1986 program for the array seismic studies of the continental lithosphere seismic experiment [J]. J Geophys Res,95:21,823～21,842
Birch F. 1961. The velocity of compression waves in rocks to 10 Kilobars [J], part 2. J Geophys Res, 66:2199～2224
Bishop T N, Bube K P, Cutler R T, et al. 1985.Tomographic determination of velocity and depth in laterally varying media [J]. Geophysics, 50:903～923
Boschetti F, Dentith D M, List D R. 1996. A fractal-based algorithm for detecting first arrivals on seismic traces [J]. Geophysics, 61:1905～1102
Braile L W. 1977.Interpretation of crustal velocity gradients and Q structure using amplitude-corrected seismic refraction profiles[A]. In: Heacock J G. The Earth's Crust, Geophys Monogr Ser[C]. AGU Washington D C, 20:427～439
Braun J A, Smith R B.1989. Two-dimensional inversion for crustal structure using refraction/wide-angle reflection data from the 1986 PASSCAL Basin-Range experiment [J]. EOS Trans Am geophys Un., 70:1206
Bregman N D, Chapman C H,Bailey R C. 1989.Travel time and Amplitude Analysis in Seismic Tomography[J]. J Geophys Res,94(B6):7577～7587
Catherine R D, Woodgold 1990.Estimation of Q In Eastern Canada Using Coda Waves[J].Bull Seisin Soc Am,80(2):411～429
Cerveny V, Molotkov I A, Psencik I.1977. Ray Method in Seismology [M], Univerzita Karlova, Prague.
Cerveny V, Hron F. 1980. The ray series method and dynamic ray tracing system for three-dimensional inhomogeneous media[J].Bull Seism Soc Am,70(1):47～77
Cerveny V., Propov M. M., Psencik I. .1982. Computation of wave fields in inhomogeneous media-Gaussian beam approach [J]. Geophys. J. R. astr. Soc., 70:109～128
Cerveny, V. and Psencik, I. 1984. SEIS83-numerical modeling of seismic wavefield in 2-D laterally varying layered structures by the ray method, In: E.R. Engdahl (ed,), Documentation of Earthquake Algorithm, World Data Center (A) for Solid Earth Geophysics, Boulder, Colo,, Rep. SE-35, pp.36～40.
Chang X, Liu Y K, Ashida Y. 1999. Hausdorff fractal algorithm for picking first break in seismic traces [J]. BUTSURI-TANSA 52:316～322
Chapman C H, Drummond R. 1982. Body-wave seismograms in heterogeneous media using Maslov asymptotic theory[J], Bull Seism Soc Am, 72: 277～417.
Chapman C H. 1985. Ray theory and its extensions: WKBJ and Maslov seismograms[J], J. Geophys., 58:27～43
Cheng Ning-ya,House L. 1996. (Short Note) Minimum traveltime calculation in 3-D graph theory[J].Geophys,61(6): 1895～1898
Coppens F. 1985. First arrival picking on common-offset trace collections for automatic estimation of static corrections [J]. Geophys. Prosp., 33:1212～1213
Cotte N,Pedersen H,Campillo M,et al. 1999.Determination of the crustal structure in southern Tibet by dispersion and amplitude analysis of Rayleigh waves[J].Geophys J Int,138:809～819
Daily W D. 1984. Underground oil-shale retort monitoring using geotomography [J]. Geophysics, 49:1701～1711
Douma H,Kabir Roy-Chowdhury. 2001. Amplitude effects due to multiscale impedance contrasts and multiple scattering: implications for Ivrea-type continental lower crust [J].Geophys J Int, 147:435～448
Du Zhi-jun,Foulger G R. 2001. Variation in the crustal structure across central Iceland[J]. Geophys J Int, 145:246～264
Dziewonski A M. 1984. Mapping the lower mantle: determination of lateral heterogeneity in P velocity up to degree and order 6. J Geophys Res., 89: 5929～5952.
Eppstein M J,Dougherty D E.1998.Optimal 3-D traveltime tomography [J] . Geophysics, 63 (3): 1053～1061
Firbas P.1987. Tomography from seismic profiles, in: Seismic Tomography. Nolet G., Reidel, Dordrecht. 189～202.
Funck T, Louden K E,Reid I D. 2000. Wide-angle seismic imaging of a Mesoproterozoic anorthosite complex:The Nain Plutonic Suite in Labrador, Canada [J]. J Geophys Res, 105 (B11): 25,693～25,707
Giese P, Prodehl C. 1976. Problems of evaluation of seismic-refracrion data for crustal and upper-mantle studies[A]. In: Giese p,et al. Explosion Seismology in Central Europe[C]. Springer, 146～147.
Hole J A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography [J]. J Geophys Res, 97: 6553～6562.
Huang Huai-yong, Spencer C,Green A. 1986.A Method for the Inversion of Refraction and Reflection Travel Times for Laterally Varying Velocity Structures[J]. Bull Seism Soc Am.,76(3):837～846
Humphreys E,Clayton R W. 1988. Adaptation of back projection tomography to seismic travel time problems[J]. J Geophys Res, 93:1 073～1 085
Inoue H, Fukao Y, Tanabe K, et al. 1990. Whole mantle P-wave travel time tomography [J]. Phys Earth Planet Inter, 59:294～328
Jin An-shu,Aki K. 1988. Spatial and Temporal Correlation Between Coda Q And Seismicity in China[J]. Bull Seism Soc Am,78(2):741～769
Kan Rong-ju J, Hu Hong-xiang, Zeng Rong-sheng, Mooney W D, et al. 1986. Crustal structure of Yunnan Province, People's Republic of China, from seismic refraction profiles [J]. Science, 234: 433～437.
Konstantin N. Poplavskii,Yuri Y. Podladchikov, Randell A. Stephenson. 2001. Twodimensional inverse modeling of sedimentary basin subsidence [J]. J Geophys Res, 106 (B4): 6657～6671
Kruger J M, Keller G R, 1986. Interpretation of crustal structure from regional gravity anomalies, Ouachita Mountains area and adjacent Gulf coastal plain[J], Am Assoc Pet Geol Bull, 70:667～689
Lees J M, Vandecar J C. 1991.Seismic tomography constrained by Bouger gravity anomalies: applications in western Washington [J],pure appli Geophys, 135:31～51.
Ludwing W J, Nafe J E, Drake C L. 1970. seismic refraction,in the sea[M], part 1. A. E. Maxwell, wisely interscience,New York. 4:53～84 Luetgert J.1992.Macray-interactive two-dimensional seismic raytracing for the Macintosh [M]. USGS Open File Rept., 92～356
Lutter W J,Nowack R L, Braile L W. 1990. Inversion for crustal structure using reflections from the PASSCAL Ouachita experiment [J]. J Geophys Res, 95:4633～4646
McCaughey M, Singh S C. 1997. Simultaneous velocity and interface tomography of normal-incidence and wide-aperture traveltime data [J]. Geophys. J. Int., 131: 87～99.
McCaughey M,Barton P J, Singh S C. 2000. Joint traveltime inversion of wide-angle seismic data and a deep reflection profile from the central North Sea [J]. Geophys J Int, 141:100～114
McMechan G. A., Mooney W. D. 1980.Asymptotic ray theory and synthetic seismograms for laterally varying structures: theory and application to the Imperial Valley, California [J]. Bull. Seism. Soc. Am., 70:2021～2035
Menke W. 1984.Geophysical Data Analysis: Discrete Inverse Theory[M], Academic Press, Orlando.
Michelena R J,Harris J M. 1991.Tomographic traveltime inversion using natural pixels [J]. Geophysics, 56 (5):635～644
Mochizuki E. 1997.Nonuniqueness of Two-Dimensional Anisotropic Tomography[J]. Bull Seism Soc Am, 87(1):261～264
Mooney W D, Brocher T M. 1987. Coincident seismic reflection refraction studies of the continental lithosphere: a global review[J]. Rev Geophys.,25:723～732.
Mooney W D.1989.Seismic methods for determining earthquake source parameters and lithospheric structure[A]. In:Pakiser L C, Mooney W D. Geophysical Framework of the Continental United States[C]. Boulder, Colorado: Geophysical Society of America Memoir, 172:11～34
Morgan R P L, Barton P J, Warner M,et al .2000.Lithospheric structure north of Scotland-I. P-wave modelling, deep reflection profiles and gravity[J].Geophys J Int,142:716～736
Moser T J. 1991. Shortest path calculation of seismic rays[J]. Geophys,56(1):59～67
Nakajima J, Matsuzawa T, Hasegawa A, et al. 2001. Three-dimensional structure of Vp, Vs, and Vp/Vs beneath northeastern Japan: Implications for arc magmatism and fluids [J]. J. Geophys. Res, 106 (B10):21,843～21,857
Neves F A, Singh S C, Priestley K.2001.Velocity structure of the upper mantle discontinuities beneath North America from waveform inversion of broadband seismic data using a genetic algorithm [J]. J Geophys Res, 106 (B10): 21,883～21,895
Oprsal I, Zahradnik J. 1999. Elastic finite-difference method for irregular grids [J]. Geophysics, 64(1):240～250
Parker R L. 1972.The Rapid Calculation of Potential Anomalies[J]. Geophys J R astr Soc, 31:447～455
Pratt R G, Worthington M H. 1988. The application of diffraction tomograply to crossbole seismic data [J]. Geophysics, 53:1284～1294
Pratt R G.. 1999.Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model [J]. Geophys, 64(3): 888～901
Qin Fu-hao, Luo Yi, Olsen K B, et al. 1992. Finite-difference solution of the eikonal equation along expanding wavefronts[J]. Geophysics,57(3):478～487
Ramananantoandro R, Bemitsas N. 1987. A computer algorithm for automatic picking of refraction first-arrival time [J].Geoexploration, 24:147～151
Ruhl T. 1996. Finite-difference migration derived from the Kirchhoff-Helmholtz integral [J]. (Short Note).Geophys,61(5): 1394～1399
Saito M. 1974.Hilbert Transforms for Sampled Data [J].J Phys Earth, 22:313～324
Sambridge M, Braun J, McQueen H. 1995. Geophysical parametrization and interpolation of irregular data using natural neighbours[J].Geophys J Int,122:837～857
Sambridge M, Gudmundsson O. 1998. Tomographic systems of equations with irregular cells [J]. J Geophys Res, 103(B1): 773～781
Sethian J A, Popovici A M. 1999. 3-D traveltime computation using the fast marching method [J]. Geophysics, 64(2):516～523
Sheriff R E, Geldart L P. 1983. Exploration Seismology Volume 2: Data-processing and interpretation [M]. Cambridge University Press
Somerstein. 1984. Radio-frequencey geotomography for remotely probing the interior of operating mini and commercial-sized oil-shale retorts [J]. Geophysics, 49:1288～1300
Spence G. D., Whittall K. P., Clowes R. M. 1984. Practical synthetic seismograms for laterally varying media calculated by asymptotic ray theory [J]. Bull. Seism. Soc. Am., 74:1209～1223
Spence G. D., Clowes R. M., Ellis R. M. 1985. Seismic structure across the active subduction zone of western Canada[J].J. geophys. Res., 90:6754～6772
Starostenko V I, Kostynkevich A S, Kozlenko V G. 1988.seismogravimetric method: principles, algorithms, results[J].Geophys J,93:295～309.
Talwani M, Worzel J L, Landisman M.1959.Rapid Gravity Computations for two-dimensional bodies with application to the mendocino submarine fracture zone [J]. J Geophys Res,64(1):49～59
Talwani M, Worzel J L,Landisman M.1959a. Rapid computaions for two-dimentional bodies with application to the Mendocino submarine fracture zone [J], J Geophys Res., 64 (1): 49～59
Talwani M, Sutton G H, Worzel J L. 1959b. A crustal section across the Puerto Rico trench, J Geophys Res., 64(10): 1547～1555
Tanimoto T, Anderson D L. 1984. Mapping convection in the mantle [J]. Geophys Res Lett., 11:287～290
Tarantola A. 1987. Inverse Problem Theory [M]. Elsevier. Amsterdam.
Vesnaver A V. 1996.Irregular grids in seismic tomography and minimum-time ray tracing [J]. Geophys J Int, 126:147～165
Vidale J E. 1988. Finite-difference traveltime calculation[J]. Bull Seism Soc. Am, 78:2062～2076
Vidale J E. 1990. Finite-difference calculation of traveltimes in three dimentions [J]. Geophys., 55:521～526
Vladimir Y G, McMechan George A. 1997. Analysis of reflection traveltimes in 3-D transversely isotropic heterogeneous media [J]. Geophys., 62(6): 1884～1895
Weng S., Chen H., Zhou X., et al. 1990. Deep seismic probing of continental crust in the Lower Yangtze region, eastern China [J]. Tectonophysics, 173:299～305
Wang Yang-hua,Gregory A, Houseman. 1994. Inversion of reflection seismic amplitude data for interface geometry[J]. Geophys J Int,117: 92～110
Wesson R L. 1971. Travel time inversion for laterally inhomogeneous crustal velocity models [J]. Bull Seism Soc Am., 61:729～746
White D J. 1989.Two-dimensional seismic refraction tomography[J]. Geophys J Int,97:223～245
White D. J., Clowes R. M. 1990.Shallow crustal structure beneath the Juan de Fuca Ridge from 2-D seismic refraction tomography [J]. Geophys. J. Int., 100:349～367
Whittall K. P., Clowes R. M. 1979. A simple, efficient method for the calculation of traveltimes and ray paths in laterally inhomogeneous media[J], J Can Soc expl Geophys, 15:21～29
Woodgold C R D. 1990.Estimation of Q in eastern Canada using coda waves[J]. Bull Seismol Soc Am, 80:411～429
Woodhouse J H, Dziewonski A M. 1984. Mapping the upper mantle: three-dimensional modelling of earth structure by inversion of seismic waveforms[J], J Geopbys Res., 89:5953～5986
Yao Z S,Osypov K S,Roberts R G. 1998.Traveltime tomography using regularized recursive least squares[J]. Geophys J Int., 134:545～553
Zelt B, Ellis R M, Clowes R M, et al.1996.Inversion of three-dimensional wide-angle scismic data from the southwestern Canadian Cordillera [J]. J Geophys Res, 101: 8503～8529
Zelt C A, Ellis R M. 1988. Practical and efficient ray tracing in two-dimensional media for rapid traveltime and amplitude forward modeling[J], Can J Explor Geophys, 24: 16-31,
Zelt C A, Smith R B.1992. Scismic traveltime inversion for 2-D crustal velocity structure [J]. Geophys J Int, 108:16～34
Zelt C A,Forsyth D A.1994.Modeling wide-angle seismic data for crustal structure: Southeastern Grenville Province [J]. J Geophys Res,99(B6): 11,687～11,704
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]. J Geophys Res, 103 (B4): 7187～7210
Zelt C A. 1999.Modelling stragies and model assessment for wide-angle seismic traveltime data[J]..Geophys J Int,139:183～204
Zhang J ten, Brink U S,Toksoz M N. 1998. Nonlinear refraction and reflection traveltime tomography [J]. J geophys Res, 103:29,743～29,757
Zhang Jie,Toksoz M N. 1998. Nonlinear refraction traveltime tomography [J]. Geophys.,63 (5): 1726～1737
陈步云，高文海．1988．贾家湾—沙园剖面地壳上地幔的深部结构[A]．见：国家地震局科技监测司编．中国大陆深部构造的研究与进展[C]．北京：地质出版社，154～168
陈国英，宋仲和，安昌强，等．1995．中国北部及其邻区地壳上地幔三维速度结构[J]．地球物理学报，38：321～328
陈沪生．1988．下扬子地区 HQ-13 线的综合地球物理调查及其地质意义[J]．石油与天然气地质，9(3)：211～222
陈沪生，周雪清，李道琪，等．1992．中国东部灵壁—奉贤(HQ—13)地学断面图1：1，000，000说明书[M]．北京：地质出版社．
陈培善，刘福田，李强，等．1990．云南地区速度结构的横向不均匀性[J]．中国科学，B辑(4)：431～439
陈天培，秦德厚，韩文华，等．1969．中华人民共和国区域地质调查报告．1：20万，建水幅[M]．昆明：云南省地质局，119
陈希孺．1992．概率论与数理统计[M]．合肥：中国科学技术大学出版社，138．
陈学波，吴跃强、杜平山等．1988．龙门山构造带两侧地壳速度结构特征 [A]．见：国家地震局科技监测司编．中国大陆深部构造的研究与进展[C]．北京：地质出版社，97～113
陈学波，李金森，赵静娴，等．1994．三江构造带基底构造特征与强震活动关系[A]．见：陈运泰，阚荣举，滕吉文主编．中国固体地球物理学进展[C]．北京：海洋出版社，41～54
陈运泰，滕吉文，张中杰，等．2001．地球物理学的回顾与展望[J]．地球科学进展，16(5)：634～642
崔作舟，卢德源，陈纪平，等．1987．攀西地区的深部地壳结构与构造[J]．地球物理学报，30：566～580
丁志峰，何正勤，孙为国，等．1999a．青藏高原东部及其边缘地区的地壳上地幔三维速度结构[J]．地球物理学报，42(2)：197～205
丁志峰．1999b．近震层析成像的理论及应用[D]．40
冯锐．1985．中国地壳厚度及上地幔密度分布(三维重力反演结果)[J]．地震学报，7(2)：143～157
冯锐，郭强绪，江月昌，等．1986．地震测深的观测与回放系统[J]．地震学报，8：57～71
冯锐，郑书真．1987．对华北岩石圈构造结构的综合研究[J]．科学通报，22：1723～1727
冯锐，陶裕录．1990．地震—重力联合反演[A]．见：中国地球物理学会年刊(1990)[C]．北京：地震出版社，44
冯锐，陶裕录．1993．地震—重力联合反演中的非块状一致性模型[J]．地球物理学报，36：463～474
傅淑芳，朱仁益．1998．地球物理反演问题[M]．北京：地震出版社，10．
国家地震局《深部物探成果》编写组．1986．中国地壳上地幔地球物理探测成果[M]．北京：地震出版社．
胡鸿翔，陆涵行，王椿镛，等．1996a．云南地区地壳结构的爆破地震研究[J]．地球物理学报，29(2)：133～144
胡鸿翔，陈学波，张碧秀，等．1986b．中国中原地区随县—安阳剖面深地震测深资料的解释[J]．地震学报，8(1)：38～49
胡鸿翔，高世玉．1993．云南地区地壳浅部基底速度细结构的研究[J]．中国地震，19(4)：356～363
胡鸿翔，李学清．1994．云南地区孟连—思茅—马龙剖面的基底速度系结构[A]．见：陈运泰，阚荣举，滕吉文主编．中国固体地球物理学进展[C]．北京：海洋出版社，100～106
黄怀勇，张碧秀，朱建亚．1989．安阳—泗水DSS剖面二维速度结构[A]．见：汤永安主编，华北平原南部地区深部地球物理综合探测方法与成果[C]．北京：海洋出版社，64～78
江苏省地震局深部构造组．1986．HQ—13测线人工地震测深资料图集．
江苏省重工业局区测队革命委员会．1970．中华人民共和国地质部区域地质测量报告书(1：200，000)，扬州幅，汤山幅，上党幅[M]．
江苏省地质局．1976．中华人民共和国区域地质调查报告(1：200，000)，无锡幅，苏州幅[M]．
江苏省地质局．1979．中华人民共和国区域地质调查报告(1：200，000)，盱眙幅[M]．
金安蜀，刘福田，孙永智．1980．北京地区地壳和上地幔三维P波速度结构[J]．地球物理学报，23：172～182
阚荣举，林中洋．1986．云南地壳上地幔构造的初步研究[J]．中国地震，2(4)：50～61
阚荣举，韩源．1992．云南遮放至马龙地学断面说明书[M]．北京：地震出版社。
阚荣举，赵晋明．1994．腾冲火山地区的壳幔构造[A]．见：陈运泰，阚荣举，滕吉文，等主编．中国固体地球物理学进展[C]．北京：海洋出版社，23～30
赖晓玲，张先康，郑需要．1997．利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态[J]．地震学报，19(5)：506～516
李成猛．1996．地震资料的加权约束反演方法[J]．石油物探，35(4)：29～36
李清河，郭建康，李流玉，等．1989．清丰—汶上DSS剖面的Q值结构[A]．见：深部地球物理综合探测方法与成果[C]．北京：海洋出版社，78～100
李庆忠．1993．走向精确勘探的道路[M]．北京：石油工业出版社，38．
李松林，樊计昌，吴宁远．1995．利用地震层析成像技术确定横向非均匀介质中速度分布和界面位置[J]．地震研究，18：430～444
李松林，吴宁远，宋占隆，等．1997．速度分布和界面位置的联合反演[J]．地震学报，19(4)：383～392
李雄．1991．地震—重力联合反演—不同方式的实现与对比[A]．见：中国地球物理学会年刊(1991)[C]．北京：地震出版社，7．
林中洋，胡鸿翔，张文彬，等．1993．云南地区地壳上地幔速度结构特征的研究[J]．地震学报，15(4)：427～440
刘宝诚，佟伟，等．1986．腾冲火山地热区的微震观测[J]．地球物理学报，29：547～555
刘昌铨，祝治平，李捍东，等．1983．连云港—临沂—泗水测深剖面及临沂8.5级地震深部构造背景[J]．地震，3：172～183
刘福田．1984．震源位置和速度结构的联合反演(Ⅰ)—理论和方法[J]．地球物理学报，27(2)：167～175
刘福田．1985．我国地球物理反演问题研究的回顾[J]．31(增刊)：144～149
刘福田，曲克信，吴华，等．1986．华北地区的地震层面成象[J]．地球物理学报，29(5)：442～449
刘福田，李强，吴华，等．1989a．用于速度图像重建的层析成像法[J]．地球物理学报，32：46～61
刘福田，曲克信，吴华，等。1989b．中国大陆及其邻近地区的地震层析成象[J]．地球物理学报，32(3)：281～293
刘福田．1995．地球物理学中的反演问题[M]．讲义，17
刘建华，刘福田，吴华，等．1989．中国南北带地壳和上地幔的三维速度图象[J]．地球物理学报，32(2)：143～153
刘清林．1993．地震初至波射线路径的追踪[J]．石油物探，32(2)：14～20
刘清林．1995．程函方程差分法层析成像[J]．石油物探，34(4)：14～26
刘清林．1996．初始模型和约束条件对层析反演结果的影响[J]．石油物探，35(2)：35～43
刘瑞丰，陈培善，李强．1993．云南及其邻近地区三维速度图象[J]．地震学报，15(1)：61～67
刘天佑．1992．重磁异常反演理论与方法[M]．中国地质大学出版社，21
楼海，王椿镛．1999．三维连续密度分布的重力计算及应用[J]．地震学报，21(3)：298～304
楼海，王椿镛，皇甫岗，等．2002．云南腾冲火山上部地壳三维地震速度层析成像[J]．地震学报，24(3)：243～251
马杏垣等．中国岩石圈动力学图集．北京：中国地图出版社，51
卢军，潘卫平，曾晓献．1991．人工爆炸地震野外数据处理装置[J]．地球物理学报，34：717～723
穆治国，佟伟．1987．腾冲火山活动的时代和岩浆来源问题[J]．地球物理学报，30：261～270
秦大正，刘昌森，丁颂华，等．1983．长江断裂带东延问题及裂谷特性的讨论[J]．地震地质，5(2)：70～77
全幼黎，滕吉文．1986．混合法计算GB体波理论地震图[J]．地球物理学报，31：165～175
任纪舜，姜春发，张正坤，等．1980．中国大地构造及其演化[M]．北京：科学出版社．
孙洁等．1989．滇西地区地壳上地幔电性结构与地壳构造运动的关系[J]．地震地质，11：35～45
孙武城，刘昌铨，祝治平，等．1984．揭示郯庐大断裂地震深部的某些特征—连云港—临沂—泗水人工地震测深资料的初步解释[A]．见：构造地质论丛[C]．北京：地质出版社，238～251
孙武城，祝治平，张利，等．1988．对华北地壳上地幔的探测与研究[A]．见：国家地震局《深部物探成果》编写组．中国地壳上地幔地球物理探测成果[C]．北京：地震出版社：19～37
孙武城，徐杰，杨主恩，等．1992．上海奉贤至内蒙古阿拉善左旗地学断面说明书[M]．北京：地震出版社．
汤永安．1989．华北平原南部地区深部地球物理综合探测方法与成果[M]．北京：海洋出版社．
滕吉文，熊绍柏，尹周勋，等．1983．喜马拉雅北部地区的地壳结构模型和速度分布特征[J]．地球物理学报，26：525～540
滕吉文，孙克忠，熊绍柏，等．1985．中国东部马鞍山—常熟—启东地带地壳与上地幔结构和速度分布的爆炸地震研究[J]．地球物理学报，28(2)：155～169
滕吉文，刘福田，全幼黎，等，1994a．中国西北造山带与沉积盆地地区的地壳和地幔的层析成像[A]．见：陈运泰等主编．中国固体地球物理学进展[C]．北京：海洋出版社，66～80
滕吉文．1994b．全球地球科学大断面与地球动力学(I)[J]．地球物理学进展，9 (1)：40～47
田宗勇，吴建成，刘家琦，等．1997．地震反射层析成像方法的研究及其在工程中的应用[J]．物探与化探，21(5)：394～397
王椿镛．1982．层状不均匀介质中合成地震图的反射率法[J]．地球物理学报，25：424～433
王椿镛．1985．二维射线合成地震图的计算及平凉—渭南地震测深剖面的解释[J]．地震研究，8(6)：607～616
王椿镛．1986．层状不均匀衰减介质中的反射率法[J]．地震学报，8(增刊)：85～94
王椿镛，陆原，陈光英，等．1989．地震测深数据处理软件包[J]．地震地磁观测和研究，10(5)：50～58
王椿镛，陆原，谢剑波，等．1991．地震测深数据文件系统[J]．地球物理学报，34：455～464
王椿镛，楼海，宋亦青．1994．人工地震测深数据库系统[J]．地震学报，16：89～95
王椿镛．1994c．中国深地震反射剖面探测与研究[A]．见：陈运泰等主编．中国固体地球物理学进展[C]．北京：海洋出版社，120～127
王椿镛，林中洋．1995a．地壳地震学研究进展[A]．见：陈运泰主编．地球与空间科学观测技术进展[C]．北京：地震出版社，272～282
王椿镛，林中洋，陈学波．1995b．青海门源—福建宁德地学断面综合地球物理研[J]．地球物理学报，38(5)：590～598
王椿镛，张先康，丁志峰，等．1997a．大别造山带上部地壳结构的有限差分层析成像[J]．地球物理学报，40(4)：496～501
王椿镛．1997b．中国岩石层结构研究的回顾与展望[J]．地球物理学报，40(增刊)：82～109
王椿镛．1998．地球内部构造和动力学的地震学研究[J]．地学前缘(中国地质大学，北京)，5(1～2)：91～98
王椿镛，楼海，吴建平等．2001a．腾冲火山地热区地壳结构的地震学研究[J]．地震学报，24(3)：
王椿镛，楼海，魏修成，等．2001b．天山北缘的地壳结构和1906年玛纳斯地震的地震构造[J]．地震学报，23(5)：460～470
王椿镛，Mooney W D, 王溪莉，等．2002a．川滇地区地壳上地幔三维速度结构研究[J]．地震学报，24(1)：1～16
王椿镛，楼海，吴建平，等．2002b．腾冲火山地热区地壳结构的地震学研究[J]．地震学报，24(3)：231～242
王华忠，谢海兵，马在田．1997．有限差分法地震波走时计算[J]．同济大学学报，25(3)：318～321
王家映．1998．地球物理反演理论[M]．武汉：中国地质大学出版社，
谢海兵，马在田，王燕琨．1997．一种新程函方程差分法层析成像[J]．石油地球物理勘探，32(5)：623～627
熊绍柏，滕吉文，尹周勋，等．1986．攀西构造带南部地壳上地幔结构和爆破地震研究[J]．地球物理学报，29(3)：235～244
杨文采．1989．地球物理反演和地震层析成像[M]，北京：地质出版社
杨文采，李幼铭，等，1992．应用地震层析成像[M]，北京：地质出版社
杨晓春，李小凡，张美根．2001．地震波反演方法研究的某些进展及其数学基础[J]．地球物理学进展，16(4)：96～109
云南省地质局．1976．中华人民共和国区域地质调查报告(1：200，000)，普洱幅[M]．
曾融生，阚荣举，何传大．1960．柴达木盆地低频地震探测结晶基底的工作方法[J]．科学通报，10：313～316
曾融生，阚荣举．1961a．柴达木盆地西部地壳深界面反射波[J]．地球物理学报，10：120～125
曾融生，阚荣举，何传大．1961b．柴达木盆地地震探测的基岩首波和大角度反射[J]．地球物理学报，10：54～66
曾融生，朱介寿，周兵，等．1992．西藏高原及其东部邻区的三维地震波速度结构与大陆碰撞模型．地震学报，14(增刊)：523～533
张少泉，俞玉琪，高世玉，等．1983．横向非均匀介质的地震射线追踪及其应用[J]．地震地磁观测与研究，4(6)：1～8
张少泉，武利均，郭建明，等．1985．中国西部地区门源—平凉—渭南地震测深剖面资料的分析解释[J]．地球物理学报，28：349～361
张四维，张锁喜，唐荣余，等．1988．下扬子地区符离集—奉贤地震测深资料解释[J]．地球物理学报，31(6)：637～647
张文生，何樵登．1996．有限差分法精确计算地震走时[J]．物探化探计算技术，18 (4)：331～335
张先康，杨卓欣，杨玉春，等．1995．地壳三维结构的层析成像方法—爆炸和地震资料的联合反演[J]．地震学报，17(4)：422～431
张先康，张成科，赵金仁，等．2002．长白山天池火山区岩浆系统深部结构的深地震测深研究[J]．地震学报，24(2)：135～143
张元生，李清河，徐果明．1998．联合利用走时与波形反演技术研究地壳三维速度结构(I)—理论与方法[J]．西北地震学报，20(2)：8～51
张中杰，滕吉文，张霖斌，等．1996．地球深部研究的复杂性及其对策[J]．地球物理学进展，11(3)：11～20
张中杰，滕吉文，李英康．2002．藏南地壳速度结构与地壳物质东西向“逃逸”—以佩枯错—普莫雍错宽角反射剖面为例[J]．中国科学，32(10)：793～798
赵爱华，张中杰，王光杰，等．2000．非均匀介质中地震波走时与射线路径快速计算技术[J]．地震学报，22(2)：151～157
周兵，朱介寿，S. A. Greenhalgh, 等．1994．一种新的地震射线层析成像计算方法[J]．石油物探，33(1)：45～54
朱介寿，严忠琼，曹俊兴，等．1998．勘探地球物理层析成像软件系统及其应用[J]．物探与化探，22(2)：90～98
朱金明，王丽燕．1992．地震波走时的有限差分法计算[J]．地球物理学报．35(1)：86～92
朱良保．1989．计算横向非均匀介质体波理论地震图的 Maslov 渐近法[J]．地球物理学报，32：163～172
朱仁益．1990．用WKBJ法计算衰减介质中的理论地震图．地球物理学报．33(增刊)：212～224．
朱天飞，刘福田．1982．用地震资料测定岩石层的三维速度结构[J]．地球物理学报．25：500～508．
Ammon C J, Vidale J E. 1993.Tomography without rays[J]. Bull Seism Soc Am, 83:509～528
Anderson D I, Dziewonski A M. 1984.Seismic tomography[J], Scientific American., 251(4):58～66
Benz H M, Smith R B, Mooney W D.1990.Crustal structure of the northwestern basin and range province from the 1986 program for the array seismic studies of the continental lithosphere seismic experiment [J]. J Geophys Res,95:21,823～21,842
Birch F. 1961. The velocity of compression waves in rocks to 10 Kilobars [J], part 2. J Geophys Res, 66:2199～2224
Bishop T N, Bube K P, Cutler R T, et al. 1985.Tomographic determination of velocity and depth in laterally varying media [J]. Geophysics, 50:903～923
Boschetti F, Dentith D M, List D R. 1996. A fractal-based algorithm for detecting first arrivals on seismic traces [J]. Geophysics, 61:1905～1102
Braile L W. 1977.Interpretation of crustal velocity gradients and Q structure using amplitude-corrected seismic refraction profiles[A]. In: Heacock J G. The Earth's Crust, Geophys Monogr Ser[C]. AGU Washington D C, 20:427～439
Braun J A, Smith R B.1989. Two-dimensional inversion for crustal structure using refraction/wide-angle reflection data from the 1986 PASSCAL Basin-Range experiment [J]. EOS Trans Am geophys Un., 70:1206
Bregman N D, Chapman C H,Bailey R C. 1989.Travel time and Amplitude Analysis in Seismic Tomography[J]. J Geophys Res,94(B6):7577～7587
Catherine R D, Woodgold 1990.Estimation of Q In Eastern Canada Using Coda Waves[J].Bull Seisin Soc Am,80(2):411～429
Cerveny V, Molotkov I A, Psencik I.1977. Ray Method in Seismology [M], Univerzita Karlova, Prague.
Cerveny V, Hron F. 1980. The ray series method and dynamic ray tracing system for three-dimensional inhomogeneous media[J].Bull Seism Soc Am,70(1):47～77
Cerveny V., Propov M. M., Psencik I. .1982. Computation of wave fields in inhomogeneous media-Gaussian beam approach [J]. Geophys. J. R. astr. Soc., 70:109～128
Cerveny, V. and Psencik, I. 1984. SEIS83-numerical modeling of seismic wavefield in 2-D laterally varying layered structures by the ray method, In: E.R. Engdahl (ed,), Documentation of Earthquake Algorithm, World Data Center (A) for Solid Earth Geophysics, Boulder, Colo,, Rep. SE-35, pp.36～40.
Chang X, Liu Y K, Ashida Y. 1999. Hausdorff fractal algorithm for picking first break in seismic traces [J]. BUTSURI-TANSA 52:316～322
Chapman C H, Drummond R. 1982. Body-wave seismograms in heterogeneous media using Maslov asymptotic theory[J], Bull Seism Soc Am, 72: 277～417.
Chapman C H. 1985. Ray theory and its extensions: WKBJ and Maslov seismograms[J], J. Geophys., 58:27～43
Cheng Ning-ya,House L. 1996. (Short Note) Minimum traveltime calculation in 3-D graph theory[J].Geophys,61(6): 1895～1898
Coppens F. 1985. First arrival picking on common-offset trace collections for automatic estimation of static corrections [J]. Geophys. Prosp., 33:1212～1213
Cotte N,Pedersen H,Campillo M,et al. 1999.Determination of the crustal structure in southern Tibet by dispersion and amplitude analysis of Rayleigh waves[J].Geophys J Int,138:809～819
Daily W D. 1984. Underground oil-shale retort monitoring using geotomography [J]. Geophysics, 49:1701～1711
Douma H,Kabir Roy-Chowdhury. 2001. Amplitude effects due to multiscale impedance contrasts and multiple scattering: implications for Ivrea-type continental lower crust [J].Geophys J Int, 147:435～448
Du Zhi-jun,Foulger G R. 2001. Variation in the crustal structure across central Iceland[J]. Geophys J Int, 145:246～264
Dziewonski A M. 1984. Mapping the lower mantle: determination of lateral heterogeneity in P velocity up to degree and order 6. J Geophys Res., 89: 5929～5952.
Eppstein M J,Dougherty D E.1998.Optimal 3-D traveltime tomography [J] . Geophysics, 63 (3): 1053～1061
Firbas P.1987. Tomography from seismic profiles, in: Seismic Tomography. Nolet G., Reidel, Dordrecht. 189～202.
Funck T, Louden K E,Reid I D. 2000. Wide-angle seismic imaging of a Mesoproterozoic anorthosite complex:The Nain Plutonic Suite in Labrador, Canada [J]. J Geophys Res, 105 (B11): 25,693～25,707
Giese P, Prodehl C. 1976. Problems of evaluation of seismic-refracrion data for crustal and upper-mantle studies[A]. In: Giese p,et al. Explosion Seismology in Central Europe[C]. Springer, 146～147.
Hole J A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography [J]. J Geophys Res, 97: 6553～6562.
Huang Huai-yong, Spencer C,Green A. 1986.A Method for the Inversion of Refraction and Reflection Travel Times for Laterally Varying Velocity Structures[J]. Bull Seism Soc Am.,76(3):837～846
Humphreys E,Clayton R W. 1988. Adaptation of back projection tomography to seismic travel time problems[J]. J Geophys Res, 93:1 073～1 085
Inoue H, Fukao Y, Tanabe K, et al. 1990. Whole mantle P-wave travel time tomography [J]. Phys Earth Planet Inter, 59:294～328
Jin An-shu,Aki K. 1988. Spatial and Temporal Correlation Between Coda Q And Seismicity in China[J]. Bull Seism Soc Am,78(2):741～769
Kan Rong-ju J, Hu Hong-xiang, Zeng Rong-sheng, Mooney W D, et al. 1986. Crustal structure of Yunnan Province, People's Republic of China, from seismic refraction profiles [J]. Science, 234: 433～437.
Konstantin N. Poplavskii,Yuri Y. Podladchikov, Randell A. Stephenson. 2001. Twodimensional inverse modeling of sedimentary basin subsidence [J]. J Geophys Res, 106 (B4): 6657～6671
Kruger J M, Keller G R, 1986. Interpretation of crustal structure from regional gravity anomalies, Ouachita Mountains area and adjacent Gulf coastal plain[J], Am Assoc Pet Geol Bull, 70:667～689
Lees J M, Vandecar J C. 1991.Seismic tomography constrained by Bouger gravity anomalies: applications in western Washington [J],pure appli Geophys, 135:31～51.
Ludwing W J, Nafe J E, Drake C L. 1970. seismic refraction,in the sea[M], part 1. A. E. Maxwell, wisely interscience,New York. 4:53～84 Luetgert J.1992.Macray-interactive two-dimensional seismic raytracing for the Macintosh [M]. USGS Open File Rept., 92～356
Lutter W J,Nowack R L, Braile L W. 1990. Inversion for crustal structure using reflections from the PASSCAL Ouachita experiment [J]. J Geophys Res, 95:4633～4646
McCaughey M, Singh S C. 1997. Simultaneous velocity and interface tomography of normal-incidence and wide-aperture traveltime data [J]. Geophys. J. Int., 131: 87～99.
McCaughey M,Barton P J, Singh S C. 2000. Joint traveltime inversion of wide-angle seismic data and a deep reflection profile from the central North Sea [J]. Geophys J Int, 141:100～114
McMechan G. A., Mooney W. D. 1980.Asymptotic ray theory and synthetic seismograms for laterally varying structures: theory and application to the Imperial Valley, California [J]. Bull. Seism. Soc. Am., 70:2021～2035
Menke W. 1984.Geophysical Data Analysis: Discrete Inverse Theory[M], Academic Press, Orlando.
Michelena R J,Harris J M. 1991.Tomographic traveltime inversion using natural pixels [J]. Geophysics, 56 (5):635～644
Mochizuki E. 1997.Nonuniqueness of Two-Dimensional Anisotropic Tomography[J]. Bull Seism Soc Am, 87(1):261～264
Mooney W D, Brocher T M. 1987. Coincident seismic reflection refraction studies of the continental lithosphere: a global review[J]. Rev Geophys.,25:723～732.
Mooney W D.1989.Seismic methods for determining earthquake source parameters and lithospheric structure[A]. In:Pakiser L C, Mooney W D. Geophysical Framework of the Continental United States[C]. Boulder, Colorado: Geophysical Society of America Memoir, 172:11～34
Morgan R P L, Barton P J, Warner M,et al .2000.Lithospheric structure north of Scotland-I. P-wave modelling, deep reflection profiles and gravity[J].Geophys J Int,142:716～736
Moser T J. 1991. Shortest path calculation of seismic rays[J]. Geophys,56(1):59～67
Nakajima J, Matsuzawa T, Hasegawa A, et al. 2001. Three-dimensional structure of Vp, Vs, and Vp/Vs beneath northeastern Japan: Implications for arc magmatism and fluids [J]. J. Geophys. Res, 106 (B10):21,843～21,857
Neves F A, Singh S C, Priestley K.2001.Velocity structure of the upper mantle discontinuities beneath North America from waveform inversion of broadband seismic data using a genetic algorithm [J]. J Geophys Res, 106 (B10): 21,883～21,895
Oprsal I, Zahradnik J. 1999. Elastic finite-difference method for irregular grids [J]. Geophysics, 64(1):240～250
Parker R L. 1972.The Rapid Calculation of Potential Anomalies[J]. Geophys J R astr Soc, 31:447～455
Pratt R G, Worthington M H. 1988. The application of diffraction tomograply to crossbole seismic data [J]. Geophysics, 53:1284～1294
Pratt R G.. 1999.Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model [J]. Geophys, 64(3): 888～901
Qin Fu-hao, Luo Yi, Olsen K B, et al. 1992. Finite-difference solution of the eikonal equation along expanding wavefronts[J]. Geophysics,57(3):478～487
Ramananantoandro R, Bemitsas N. 1987. A computer algorithm for automatic picking of refraction first-arrival time [J].Geoexploration, 24:147～151
Ruhl T. 1996. Finite-difference migration derived from the Kirchhoff-Helmholtz integral [J]. (Short Note).Geophys,61(5): 1394～1399
Saito M. 1974.Hilbert Transforms for Sampled Data [J].J Phys Earth, 22:313～324
Sambridge M, Braun J, McQueen H. 1995. Geophysical parametrization and interpolation of irregular data using natural neighbours[J].Geophys J Int,122:837～857
Sambridge M, Gudmundsson O. 1998. Tomographic systems of equations with irregular cells [J]. J Geophys Res, 103(B1): 773～781
Sethian J A, Popovici A M. 1999. 3-D traveltime computation using the fast marching method [J]. Geophysics, 64(2):516～523
Sheriff R E, Geldart L P. 1983. Exploration Seismology Volume 2: Data-processing and interpretation [M]. Cambridge University Press
Somerstein. 1984. Radio-frequencey geotomography for remotely probing the interior of operating mini and commercial-sized oil-shale retorts [J]. Geophysics, 49:1288～1300
Spence G. D., Whittall K. P., Clowes R. M. 1984. Practical synthetic seismograms for laterally varying media calculated by asymptotic ray theory [J]. Bull. Seism. Soc. Am., 74:1209～1223
Spence G. D., Clowes R. M., Ellis R. M. 1985. Seismic structure across the active subduction zone of western Canada[J].J. geophys. Res., 90:6754～6772
Starostenko V I, Kostynkevich A S, Kozlenko V G. 1988.seismogravimetric method: principles, algorithms, results[J].Geophys J,93:295～309.
Talwani M, Worzel J L, Landisman M.1959.Rapid Gravity Computations for two-dimensional bodies with application to the mendocino submarine fracture zone [J]. J Geophys Res,64(1):49～59
Talwani M, Worzel J L,Landisman M.1959a. Rapid computaions for two-dimentional bodies with application to the Mendocino submarine fracture zone [J], J Geophys Res., 64 (1): 49～59
Talwani M, Sutton G H, Worzel J L. 1959b. A crustal section across the Puerto Rico trench, J Geophys Res., 64(10): 1547～1555
Tanimoto T, Anderson D L. 1984. Mapping convection in the mantle [J]. Geophys Res Lett., 11:287～290
Tarantola A. 1987. Inverse Problem Theory [M]. Elsevier. Amsterdam.
Vesnaver A V. 1996.Irregular grids in seismic tomography and minimum-time ray tracing [J]. Geophys J Int, 126:147～165
Vidale J E. 1988. Finite-difference traveltime calculation[J]. Bull Seism Soc. Am, 78:2062～2076
Vidale J E. 1990. Finite-difference calculation of traveltimes in three dimentions [J]. Geophys., 55:521～526
Vladimir Y G, McMechan George A. 1997. Analysis of reflection traveltimes in 3-D transversely isotropic heterogeneous media [J]. Geophys., 62(6): 1884～1895
Weng S., Chen H., Zhou X., et al. 1990. Deep seismic probing of continental crust in the Lower Yangtze region, eastern China [J]. Tectonophysics, 173:299～305
Wang Yang-hua,Gregory A, Houseman. 1994. Inversion of reflection seismic amplitude data for interface geometry[J]. Geophys J Int,117: 92～110
Wesson R L. 1971. Travel time inversion for laterally inhomogeneous crustal velocity models [J]. Bull Seism Soc Am., 61:729～746
White D J. 1989.Two-dimensional seismic refraction tomography[J]. Geophys J Int,97:223～245
White D. J., Clowes R. M. 1990.Shallow crustal structure beneath the Juan de Fuca Ridge from 2-D seismic refraction tomography [J]. Geophys. J. Int., 100:349～367
Whittall K. P., Clowes R. M. 1979. A simple, efficient method for the calculation of traveltimes and ray paths in laterally inhomogeneous media[J], J Can Soc expl Geophys, 15:21～29
Woodgold C R D. 1990.Estimation of Q in eastern Canada using coda waves[J]. Bull Seismol Soc Am, 80:411～429
Woodhouse J H, Dziewonski A M. 1984. Mapping the upper mantle: three-dimensional modelling of earth structure by inversion of seismic waveforms[J], J Geopbys Res., 89:5953～5986
Yao Z S,Osypov K S,Roberts R G. 1998.Traveltime tomography using regularized recursive least squares[J]. Geophys J Int., 134:545～553
Zelt B, Ellis R M, Clowes R M, et al.1996.Inversion of three-dimensional wide-angle scismic data from the southwestern Canadian Cordillera [J]. J Geophys Res, 101: 8503～8529
Zelt C A, Ellis R M. 1988. Practical and efficient ray tracing in two-dimensional media for rapid traveltime and amplitude forward modeling[J], Can J Explor Geophys, 24: 16-31,
Zelt C A, Smith R B.1992. Scismic traveltime inversion for 2-D crustal velocity structure [J]. Geophys J Int, 108:16～34
Zelt C A,Forsyth D A.1994.Modeling wide-angle seismic data for crustal structure: Southeastern Grenville Province [J]. J Geophys Res,99(B6): 11,687～11,704
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]. J Geophys Res, 103 (B4): 7187～7210
Zelt C A. 1999.Modelling stragies and model assessment for wide-angle seismic traveltime data[J]..Geophys J Int,139:183～204
Zhang J ten, Brink U S,Toksoz M N. 1998. Nonlinear refraction and reflection traveltime tomography [J]. J geophys Res, 103:29,743～29,757
Zhang Jie,Toksoz M N. 1998. Nonlinear refraction traveltime tomography [J]. Geophys.,63 (5): 1726～1737