大华北及其邻区地壳上地幔三维速度结构的地震层析成像研究
|
摘要
本文对地球物理反演的发展历史和研究现状作了简单回顾,交待了大华北及其邻区的构造环境,阐述了本文所采用的反演方法的基本原理,分析了资料的可靠性和准确性,选择了适当的反演模型,得出了较好的反演结果,并对本地区的反演分辨率进行了讨论。
本文的研究区域是大华北及其邻区(107°E-124°E,28°N-42°N),即通常所指的大华北地区及其以南部分地区。这个区域是我国地震活动强烈的地区之一,曾多次发生强烈破坏性地震,区内既有大陆内部的碰撞造山带(秦岭-大别造山带),又有贯通本区南北的大型断裂带(郯庐断裂带)和大型重力梯级带(太行山),地质构造十分复杂。
本文利用区域地震到时资料、部分人工地震资料和远震的P波资料,使用赵大鹏博士提供的程序,用体波地震成像技术重建了本地区地壳上地幔的三维速度结构。对所使用的资料的可靠性作了充分的分析,对资料在本地区的层析成像分辨能力作了比较,使用了适合本次工作的分辨率0.5°×0.5°,反演出了与区域大地构造相对应的较好的结果。
本文利用大华北及其邻区475个地震台站的区域地震资料(共15714个区域地震,P波震相114971个,S波震相102142个)、三条人工地震剖面初至P波走时数据(江苏8401,山东8509,大别山9511,共51炮的数据)和4975个远震(76620个远震P波)资料。以水平方向0.5°×0.5°,垂直方向以1km、10km、25km、42km、60km为间隔划分网格,并建立初始三维速度模型。用近似弯曲射线追踪方法计算走时和射线路径,用LSQR算法进行反演。在反演之前,用检测方法分析了所用地震资料的层析成像分辨率。层析成像结果较好地反映了本地区的地壳上地幔速度结构特征。
结果表明:
(1) 使用区域地震资料、远震资料和人工地震资料联合进行层析成像能够取得好的结果。
(2) 大华北及其邻区大的构造如燕山构造带、太行山前构造带、汾渭构造带、郯庐断裂带及秦岭-大别构造带等都与该区地壳上地幔的低速或高、低速度相邻的部位相对应,华北断块与扬子断块速度异常沿秦岭-大别碰撞带存在明显的分区特征,扬子断块波速异常较少,没有明显的横向不均匀性,而华北断块多
The history of geophysical inversion and seismic tomography is reviewed in this paper. The tectonic of study region and the method of inversion are explained simply also. In order to get good result, the reliability and preciseness of the seismic data are analyzed carefully and the good inversion model is selected. The resolving power in this area is discussed in the end.The research region of this dissertation is 107°E-124°E and 28°N-42°N. In China geological and tectonic map, this region is referred to as the North China and some south area of it. This region is located in a seismically very active region in China. A lot of large earthquakes have occurred in it. The present study area includes many active tectonics such as Yanshan uplift, Qinling-Dabie uplift, the fault in Taihangshan piedmont, Tanlu fault ect.In total, 114971 P arrival times, 102142 S arrival times from 15714 local earthquakes, 76620 P arrive times from 4975 teleseismic events and some data from seismic sounding profile have been collected. A detailed three-dimensional(3-D) velocity model of the crust and uppermost mantle under the region is determined with a spatial resolution of 0.5° by 0.5° in the horizontal direction and 9-20 kilometer in depth. To analyze the arrival time data, we have used the tomographic method of Zhao et al.(1992). This method is adaptable to a velocity structure which includes several complex-shaped velocity discontinuities and allows 3-D velocity variations everywhere in the model. The discontinuities represent known geological boundaries, like the Moho discontinuity. A 3-D grid net is set up in model to express the 3-D structure. The tomograhic inversion was calculated using this program. According to the data distribution and several check-board tests, reliable inversion results are obtained for a model with grid size of 0.5° by 0.5° in the horizontal directions and 5 layers of 1,10,25,42 and 60 km deep from the surface to 60 km in the depth.The tomographic inversion results show that: Yanshan uplift, Qinling-Dabie uplift, the fault in Taihangshan piedmont, Tanlu fault etc are corresponding with low velocity or between low and high velocity in the crust and uppermost mantle in this region. Depth around 300 kilometer strong P-wave velocity variations of up to 4% are found in some area such as Yanshan uplift and Qinling-Dabie uplift, indicating the
本文的研究区域是大华北及其邻区(107°E-124°E,28°N-42°N),即通常所指的大华北地区及其以南部分地区。这个区域是我国地震活动强烈的地区之一,曾多次发生强烈破坏性地震,区内既有大陆内部的碰撞造山带(秦岭-大别造山带),又有贯通本区南北的大型断裂带(郯庐断裂带)和大型重力梯级带(太行山),地质构造十分复杂。
本文利用区域地震到时资料、部分人工地震资料和远震的P波资料,使用赵大鹏博士提供的程序,用体波地震成像技术重建了本地区地壳上地幔的三维速度结构。对所使用的资料的可靠性作了充分的分析,对资料在本地区的层析成像分辨能力作了比较,使用了适合本次工作的分辨率0.5°×0.5°,反演出了与区域大地构造相对应的较好的结果。
本文利用大华北及其邻区475个地震台站的区域地震资料(共15714个区域地震,P波震相114971个,S波震相102142个)、三条人工地震剖面初至P波走时数据(江苏8401,山东8509,大别山9511,共51炮的数据)和4975个远震(76620个远震P波)资料。以水平方向0.5°×0.5°,垂直方向以1km、10km、25km、42km、60km为间隔划分网格,并建立初始三维速度模型。用近似弯曲射线追踪方法计算走时和射线路径,用LSQR算法进行反演。在反演之前,用检测方法分析了所用地震资料的层析成像分辨率。层析成像结果较好地反映了本地区的地壳上地幔速度结构特征。
结果表明:
(1) 使用区域地震资料、远震资料和人工地震资料联合进行层析成像能够取得好的结果。
(2) 大华北及其邻区大的构造如燕山构造带、太行山前构造带、汾渭构造带、郯庐断裂带及秦岭-大别构造带等都与该区地壳上地幔的低速或高、低速度相邻的部位相对应,华北断块与扬子断块速度异常沿秦岭-大别碰撞带存在明显的分区特征,扬子断块波速异常较少,没有明显的横向不均匀性,而华北断块多
The history of geophysical inversion and seismic tomography is reviewed in this paper. The tectonic of study region and the method of inversion are explained simply also. In order to get good result, the reliability and preciseness of the seismic data are analyzed carefully and the good inversion model is selected. The resolving power in this area is discussed in the end.The research region of this dissertation is 107°E-124°E and 28°N-42°N. In China geological and tectonic map, this region is referred to as the North China and some south area of it. This region is located in a seismically very active region in China. A lot of large earthquakes have occurred in it. The present study area includes many active tectonics such as Yanshan uplift, Qinling-Dabie uplift, the fault in Taihangshan piedmont, Tanlu fault ect.In total, 114971 P arrival times, 102142 S arrival times from 15714 local earthquakes, 76620 P arrive times from 4975 teleseismic events and some data from seismic sounding profile have been collected. A detailed three-dimensional(3-D) velocity model of the crust and uppermost mantle under the region is determined with a spatial resolution of 0.5° by 0.5° in the horizontal direction and 9-20 kilometer in depth. To analyze the arrival time data, we have used the tomographic method of Zhao et al.(1992). This method is adaptable to a velocity structure which includes several complex-shaped velocity discontinuities and allows 3-D velocity variations everywhere in the model. The discontinuities represent known geological boundaries, like the Moho discontinuity. A 3-D grid net is set up in model to express the 3-D structure. The tomograhic inversion was calculated using this program. According to the data distribution and several check-board tests, reliable inversion results are obtained for a model with grid size of 0.5° by 0.5° in the horizontal directions and 5 layers of 1,10,25,42 and 60 km deep from the surface to 60 km in the depth.The tomographic inversion results show that: Yanshan uplift, Qinling-Dabie uplift, the fault in Taihangshan piedmont, Tanlu fault etc are corresponding with low velocity or between low and high velocity in the crust and uppermost mantle in this region. Depth around 300 kilometer strong P-wave velocity variations of up to 4% are found in some area such as Yanshan uplift and Qinling-Dabie uplift, indicating the
引文
Aki K, Christoffersson A, Husebye E S, et al. 1974. Three-dimensional seismic velocity anomalies in the crust and upper-mantle under the U.S.GS. California seismic array (abstract), EOS, Trans. Am. Geophys. Union, 56:1145.
Aki K, Christoffersson A, Husebye E S. 1976. Three-dimensional seismic structure of the lithosphere under Montana LASA[J]. Bull. Seismol. Soc. Am.,66:501~524.
Aki K, Lee W H K. 1976. Determination of three dimensional velocity anomalies under a seismic array using first P arrival times form local earthquakes 1 .A homogeneous initial model[J]. Journal of Geophysical Research, 81: 4381-4399.
Aki K, Christoffersson A, Husebye E S. 1977. Determination of the three-dimensional seismic structure of the lithosphere[J]. Journal of Geophysical Research , 82: 277-296.
Ankeny L A, Braile L W. 1986. Upper crustal structure beneath the Jemez Mountains Volcanic Filed , New Mexica , Determined by Three-Dimensional Simultaneous inversion of seismic refraction and earthquake data[J]. Journal of Geophysical Research , 91(B6): 6188-6198.
Anderson D I, Dziewonski A M. 1984. Seismic tomography[J]. Scientific American, 251(4):58-66.
Backus G, Gilbert F. 1967. Numerical applications of a formalism for geophysical inverse problem[J]. Geophys. J. R. astr. Soc, 13:247-276,
Backus G, Gilbert F. 1968. The resolving power of gross earth data[J]. Geophys. J.RAstr.Soc, 16:169-205,
Backus G, Gilbert F. 1970. Uniqueness in the inversion of inaccurate gross earth data, Phil[J]. Trans. Roy. Soc, London, Ser. A266: 123-196,
Bishop T N, Bube K P, Citler R T, et al. 1985. Tomogrophic determination of velocity and depth in laterallg varying media[J]. Geophysics, 50:903~923,
Blumling P, Prodehl C. 1983. Crustal structure beneath the eastern part of the Coast Ranges(Diablo Range) of central California from explosion seismic and near earthquake data[J]. Physics of the Earth and Planetary Interios, 31: 313-326.
Benz H M, Smith R B. 1984. Simultaneous inversion for lateral velocity variations and hypocenters in the Yellowstone region using earthquake and refraction data[J]. Journal of Geophysical Research, 89(B2): 1208-1220.
Bullock A D, Minshull T A. 2004. Poisson's Ratio Structure Through a Zone of Exhumed Mantle at the Goban Spur Rifted Margin, Southwest of the UK[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Crosson R S. 1976. Crustal structure modeling of earthquake data 1. Simultaneous least squares estimation of hypocenter and velocity parameter[J]. Journal of Geophysical Research, 81:3036-3064.
Chun-Yong Wang, Chan W W, Mooney W D. 2000. A crustal model of the ultrahigh-pressure Dabie Shan orogenic belt, China, derived from deep seismic refraction profiling[J]. Journal of Geophysical Research, 105(B5): 10857-10869.
Chun-Yong WANG , Gui-Mei WANG , Zhong-Yang LIN , et al. 1993 . Study of fine crustal structure in Xintai earthquake region by deep seismic reflection method[J]. Chinese Journal of Geophysics , 36(5): 445-451 (in Chinese).
Christoph H A, Knackstedt M A, Pinczewskiz W V., et al. 2002. Computation of linear elastic properties from microtomographic images: Methodology and agreement between theory and experiment[J]. Geophysics, 67 (5): 1396-1405.
Dziewonski A M, Gilbert F. 1976. The effect of small aspherical perturbations on travel times and re-examination of the corrections for ellipticity[J]. Geophys. J. R. Astr. Soc, 44: 7-17.
Dziewonski A M. 1984. Mapping the lower mantle: determination of lateral heterogeneity in P velocity up to degree and order 6[J]. Journal of Geophysical Research, 89:5929-5952.
Dziewonski A M, Anderson D L. 1984 Seismic tomography of the Earth's interior[J]. Am.Sci. ,721:483~494.
DvorkinJ, PrasadM, Sakai A, et al. 1999. Elasticity of marine sediments: Rock hysics modeling[J]. Geophys. Res. Lett., 26(12): 1781-1784.
DeShon H R, Schwartz S Y, Dorman L M, et al. 2004. Seismogenic zone structure along the Middle America Trench, Nicoya Peninsula, Costa Rica, from 3D local earthquake tomography[A]. In : 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Eberhart-Philips D. 1986. Three-dimensional velocity structure in north California Coast Ranges from inversion of local earthquakes arrival times[J]. Bull Seism Soc Amer, 76: 1 025-1 052.
Farra V, Madariaga R. 1988. Non-linear reflection tomography[J]. Geophys. J. Int., 95:135-147.
Geiger L. 1912. Probalibility method for the determination of earthquake epicenters from arrival time only[J]. Bull.St.Louis. Univ, 8:60-71.
Huang J, Zhao D. 2004. Crustal heterogeneity and seismotectonics of the region around Beijing[J]. Tectonophysics ,385:159-180.
Hole J A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography[J]. Journal of Geophysical Research, 97:6553-6562.
Humphreys E, Clayton R W. 1988. Adaptation of back projection tomography to seismic travel time problems[J]. Journal of Geophysical Research, 93:1073-1085.
Husen S Smith R B, Waite G P. 2004. Evidence for gas and magmatic sources beneath the Yellowstone volcanic field from seismic tomographic imaging[J]. J. Volcanol. Geotherm. Res. 131: 397-410.
Ivanson S. 1985. A study of methods for tomographic velocity estimation in the presence of low-velocity zones[J]. Geophysics, 50:969-988,.
Ivanson S. 1986. Some remarks concerning seismic reflection tomography and velocity ananlysis[J]. Geophys. J. R.Astr.Soc,87.539-557,.
Julian B R, Gubbins D. 1977. Three-dimensional sesmic ray tracing[J]. J. Geophys, 43:95-113.
Kamei R, Minami Y, Xue Z. 2004. Evaluation of Seismic Properties and Their Dependence on CO_2 Phases in a Porous Sandstone; Insights From Laboratory Injection Experiments[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Krajewski J, Roecker S, Thurber C. 2004. Combined Teleseismic and Local Earthquake Tomography of the SAFOD Drill Site[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Koch M. 1985. A numerical study on th e determination of the 3-D structure of th lithosphere by linear and non-liner inversion of teleseismic traveltimes[J]. Geophys. J.R.Astr.Soc,80:73~93.
Koch M. 1993a. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-I.evelopment,validation and optimal regularization of the method[J]. Geophys. J.Int., 112:385-412.
Koch M. 1993b. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-II.Application to the northern Rhine Grabern/Rhanish Massif Region Germany[J]. Geophys. J.Int., 112:413-428,.
Koch M. 1993c. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-III.Application to the southern Rhine Grabern seismic regin, Germany[J]. Geophys. J.Int., 112:429-447,
Lees J M, Crosson R S. 1989. Tomographic inversion for three-dimensional velocity structure at Mount St. Helens using earthquake data[J]. Journal of Geophysical Research, 94: 5716-5728.
Lees J M, Wu H. 2000. Possion's ratio and porosity at Coso geothermal area, Californica[J]. J. Volcanol. Geotherm. Res., 95: 157-163.
Mishra O P, Zhao D. 2003. Crack density, saturation rate and porosity at the 2001 Bhuj, India, earthquake hypocenter; a fluid-driven earthquake[J]. Earth Plant. Sci. Letters, 212: 393-405.
O'Connell R J, Budiansky B. 1974. Seismic velocities in dry and saturated cracked solids[ J]. Journal of Geophysical Research,79: 5412-5426.
Pavlis G L, Booker J R. 1980. The Mixed Discrete-Continuous Inverse Problem: Application to the Simultaneous Determination of Earthquake Hypocenters and Velocity Structure[J]. Journal of Geophysical Research . 85(B9) 4801-4810.
Pickett G R. 1963. Acoustic character logs and their applications in formation evaluation[J]. J. Petrol. Tech, 15: 650~667.
Shedloc. 1987. Elastic wave velocity structure of the crust and upper mantle beneath the north China basin[J]. Journal of Geophysical Research, 92(B9):9327~9350.
Thurber C H. 1983. Earthquake Locations and Three-Dimensional Crustal structure in the Coyote Lake Area, Central California[J]. Journal of Geophysical Research, 88(B10): 8226-8236.
Um J, Thurber C. 1987. A fast algorithm for two-point seismic ray tracing[J]. Bull. Seism. Soc.Am., 77(3):972~986.
Vanorio T, Virieux J, and Lotorre D. 2004. Beyond Delay Time Tomography the Estimation of Rock physics properties from P- and S- Velocity models, In:2004 Fall Meeting, American Geophysical Union, Abs.[C], 2004.
Wagner L S, Zandt G, Beck S, et al. 2004. Black upper mantle beneath the Sierras Pampeanas, Argentina? Results from seismic tomography[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Wang Qing-Chen, Sun Shu, Li Ji-Liang et al. 1989. Tectonic evolution of QinLing mountains[J]. Scientia Geologica Sinica(in Chinese),2:129~141
Wiggins, R. A, 1972. The generalized linear inverse problem: Implication of surface waves and free oscillations for earth structure[J]. Rev. Geophys. Space Phys., 10:251-285.
Zhao D, Horiuchi S, Takgi A. 1990. 3-D seismic velocity structure of the crust and the uppermost mantle in the northeastern Japan arc, Tectonopgysics, 181:135-149.
Zhao D. 1991. A tomographic study of seismic velocity structure in the Japan Islands,Ph.D.Theis, Tohoku University.
Zhao D, Hasegawa A, Horiuchi S. 1992. Tomographic imaging of P and S wave velocity structure eneath northeastern Japan[J]. Journal of Geophysical Research, 97(B 13): 19909-19928.
Zhao D, Hasegawa A. 1993. P wave tomographic imaging of the crust and upper mantle beneath the Japan lslands[J] Journal of Geophysical Research, 98(B3):4333~4353.
Zhao D.1994. Deep structure of Japan subduction zone as deriveved from local, regional, and teleseismic events[J]. Journal of Geophysical Research, 99(B11):22313-22329.
Zhao D. 1998. The 1995 Kobe earthquake: Seismic image of the source zone and its implications for rupture nucleation[J]. Journal of Geophysical Research, 103(B5): 9967～9986.
Zhao D. 2001. Seismic structure and origin of hotspots and mantle plumes[J]. Earth. Planet. Sci. Lett., 192: 251～265.
Zhao D. 2004. Seismic imaging of North China: insight into intraplate volcanism and seismoteetonics[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Zhao D. 2004. Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics[J]. Phys. Earth Planet. Int., 146: 3～34.
Zcng Rongsheng, Wang Chunyong, Zhang Dongning. 1995a. On the dynamics of extensional basin, PAGEOPH, 145(3/4): 579-603.
陈辉,刘福田,游庆瑜,等.1999.地震体波成像中的几个问题[J].地球物理学进展,14(4):9～17.
丁国瑜.1991.中国岩石圈动力学概论[M].北京.地震出版社.
丁志峰.1999.近震层析成像的理论及应用[学位论文].北京:中国地震局地球物理研究所.
丁志峰,曾融生.1994.用近震资料反演京津唐地区的地壳三维速度结构[J].华北地震科学,12(2):14～20.
丁志峰.2001.青藏高原地震波三维速度结构的研究[J].中国地震,17(2).
地质部航空物探大队904队.1957.华北平原南部及周围山区航空磁测结果报告.
方盛明,张先康,嘉世旭,等.2002.华北地区布格重力异常的多尺度分解特征与地震活动性[J].大地测量与地球动力学,22 (1):34～39.
付英祺,杨季楷.1987.地史学简明教程[M].北京:地质出版社.
金安蜀,刘福田,孙永智.1980.北京地区地壳上地幔的三维P波速度结构[J].地球物理学报,23(2):172～182.
嘉世旭,张先康,方盛明.2001.华北裂陷盆地不同块体地壳结构及演化研究[J].地学前缘,8 (2):259～265.
嘉世旭,张先康.2005.华北不同构造块体地壳结构及其对比研究[J].地球物理学报[J],48 (3):611～620.
国家地震局地球物理研究所.1978.近震分析[M].北京:地震出版社.
国家地震局《深部物探探测成果》编写组.1986.中国地壳上地幔地球物理探测成果[M].北京:地震出版社.
国家地震局地质研究所.1987.郯庐断裂[M].北京:地震出版社.
国家地震局鄂尔多斯周缘活动断裂系课题组.1988.鄂尔多斯周缘活动断裂系[M].北京:地震出版社.
郭增建,秦保燕.1979.震源物理[M].北京:地震出版社.
何正勤,丁志峰等.2001.中国大陆及其邻域地壳上地幔速度结构的面波层析成像研究[J].地球物理学报,(06):32～39.
李松林,张先康,张成科等.2002.玛沁-兰州-靖边地震测深剖面地壳速度结构的初步研究[J].地球物理学报,45(2):210～216.
李强,王椿镛,刘瑞丰,王溪莉,陈光英,陈光,李桂银.1999.应用层析成像技术研究华北地壳速度结构[J].地震地磁观测与研究.
刘瑞丰,陈培善,李强.1993.云南及其邻近地区三维速度图象[J].地震学报,15(1):61～67.
刘福田.1984.震源位置和速度结构的联合反演(Ⅰ)-理论和方法.地球物理学报[J],27(2):167～175.
刘建华,刘福田,孙若昧,等.1995.秦岭-大别造山带及其南北缘地震层析 成像,地球物理学报[J],38 (1):46～53.
刘福田,曲克信,吴华等.1986.华北地区的地震层析成像[J].地球物理学报,29(5):442～449.
刘福田,曲克信,吴华等.1989.中国大陆及其邻近地区的地震层析成象[J].地球物理学报[J],32(3):281～291.
刘福田.1991.三维速度结构的研究现状和展望.地球物理学报[J],34 (6):788～196.
刘建忠,李三忠,周立宏,高振平,郭晓玉.2004.华北板块东部中生代构造变形与盆地格局.海洋地质与第四纪地质[J],24(4):46～54.
刘伊克,常旭.2000.地震成析成像中解的定量评价及其应用.地球物理学报[J],43(2):252～256.
马杏垣,吴正文,谭应佳等.1979.华北地台基底构造[J].地质学报,53 (4):293～304.
马杏垣.1989.中国岩石圈动力学地图集.北京:中国地图出版社.
裴顺平等.2002.新疆及邻区Pn速度层析成像.地球物理学报[J],45(2):218-224.
孙若昧,刘福田.1995.京津地区地壳结构与强震的发生—Ⅰ:P波速度结构[J].地球物理学报,38(5):599～607.
孙若昧,赵燕来,吴丹.1996,京津唐地区地壳结构与强震的发生—Ⅱ:S波速度结构[J].地球物理学报,39(3):347～355.
孙若昧,赵燕来,梅世蓉.1993.勃海及其邻区的地震成析成像[J].地球物理学报,36(1):44～54.
宋仲和,陈国英,安昌强等.1993发中国大陆及其相邻海域瑞利波群速度分布特征[J].地震学报,15(1):32～38.
滕吉文,曾融生等.2002.东亚大陆及周边海域Moho界面深度分布和基本构 造格局[J].中国科学D辑,32(2):89～100.
滕吉文,张中杰,白武明等.2004.岩石圈物理学[M].北京:科学出版社.
王椿镛,王溪莉,颜其中.1993.昆明地震台网多事件定位问题的初步研究[J].地球物理学报,15 (2):136～145.
王椿镛,丁志峰等.1997.大别山造山带地壳S波速度结构[J].地球物理学报,40(3):337～345.
王椿镛,丁志峰等.1997.大别造山带地壳S波分裂和介质各向异性[J].科学通报.
王椿镛,张先康,丁志峰,阮红,邓宏钊.1997.大别造山带上部地壳结构的有限差分层析成像[J].地球物理学报,40(4):495～502.
王椿镛,W. D. Mooner,王溪莉,吴建平,楼海,王飞.2002.川滇地区地壳地幔三维速度结构研究[J].地震学报,24(1):1～16.
王椿镛,楼海,王飞.1999.大别山超高压变质带地壳结构及其构造意义[J].地震学报,21(5):534～544.
王夫运,张先康等.2002.用地震走时反演长白山天池火山地区的二维地壳结构[J].地震学报,24(2):144～152.
王夫运,张先康等.2004.张渤构造带中东段二维P波速度结构安新—香河—宽城剖面[J].地震学报,26增刊:31～41.
徐果明.2003.反演理论及其应用[M].北京:地震出版社.
徐锡伟,吴卫民,张先康等.2002.首都圈地区地壳最新构造变动与地震[M].北京.科学出版社.
姚振兴,李白基,梁尚鸿,朱培定,张立敏,卢善声.1981.青藏高原地区瑞利波群速度和地壳构造[J].地球物理学报,24(3):287～295.
杨景春.1993.中国地貌特征与演化[M].北京:海洋出版社.
杨文采.1997.地球物理反演的理论与方法[M].北京:地质出版社.
Aki K, Christoffersson A, Husebye E S. 1976. Three-dimensional seismic structure of the lithosphere under Montana LASA[J]. Bull. Seismol. Soc. Am.,66:501~524.
Aki K, Lee W H K. 1976. Determination of three dimensional velocity anomalies under a seismic array using first P arrival times form local earthquakes 1 .A homogeneous initial model[J]. Journal of Geophysical Research, 81: 4381-4399.
Aki K, Christoffersson A, Husebye E S. 1977. Determination of the three-dimensional seismic structure of the lithosphere[J]. Journal of Geophysical Research , 82: 277-296.
Ankeny L A, Braile L W. 1986. Upper crustal structure beneath the Jemez Mountains Volcanic Filed , New Mexica , Determined by Three-Dimensional Simultaneous inversion of seismic refraction and earthquake data[J]. Journal of Geophysical Research , 91(B6): 6188-6198.
Anderson D I, Dziewonski A M. 1984. Seismic tomography[J]. Scientific American, 251(4):58-66.
Backus G, Gilbert F. 1967. Numerical applications of a formalism for geophysical inverse problem[J]. Geophys. J. R. astr. Soc, 13:247-276,
Backus G, Gilbert F. 1968. The resolving power of gross earth data[J]. Geophys. J.RAstr.Soc, 16:169-205,
Backus G, Gilbert F. 1970. Uniqueness in the inversion of inaccurate gross earth data, Phil[J]. Trans. Roy. Soc, London, Ser. A266: 123-196,
Bishop T N, Bube K P, Citler R T, et al. 1985. Tomogrophic determination of velocity and depth in laterallg varying media[J]. Geophysics, 50:903~923,
Blumling P, Prodehl C. 1983. Crustal structure beneath the eastern part of the Coast Ranges(Diablo Range) of central California from explosion seismic and near earthquake data[J]. Physics of the Earth and Planetary Interios, 31: 313-326.
Benz H M, Smith R B. 1984. Simultaneous inversion for lateral velocity variations and hypocenters in the Yellowstone region using earthquake and refraction data[J]. Journal of Geophysical Research, 89(B2): 1208-1220.
Bullock A D, Minshull T A. 2004. Poisson's Ratio Structure Through a Zone of Exhumed Mantle at the Goban Spur Rifted Margin, Southwest of the UK[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Crosson R S. 1976. Crustal structure modeling of earthquake data 1. Simultaneous least squares estimation of hypocenter and velocity parameter[J]. Journal of Geophysical Research, 81:3036-3064.
Chun-Yong Wang, Chan W W, Mooney W D. 2000. A crustal model of the ultrahigh-pressure Dabie Shan orogenic belt, China, derived from deep seismic refraction profiling[J]. Journal of Geophysical Research, 105(B5): 10857-10869.
Chun-Yong WANG , Gui-Mei WANG , Zhong-Yang LIN , et al. 1993 . Study of fine crustal structure in Xintai earthquake region by deep seismic reflection method[J]. Chinese Journal of Geophysics , 36(5): 445-451 (in Chinese).
Christoph H A, Knackstedt M A, Pinczewskiz W V., et al. 2002. Computation of linear elastic properties from microtomographic images: Methodology and agreement between theory and experiment[J]. Geophysics, 67 (5): 1396-1405.
Dziewonski A M, Gilbert F. 1976. The effect of small aspherical perturbations on travel times and re-examination of the corrections for ellipticity[J]. Geophys. J. R. Astr. Soc, 44: 7-17.
Dziewonski A M. 1984. Mapping the lower mantle: determination of lateral heterogeneity in P velocity up to degree and order 6[J]. Journal of Geophysical Research, 89:5929-5952.
Dziewonski A M, Anderson D L. 1984 Seismic tomography of the Earth's interior[J]. Am.Sci. ,721:483~494.
DvorkinJ, PrasadM, Sakai A, et al. 1999. Elasticity of marine sediments: Rock hysics modeling[J]. Geophys. Res. Lett., 26(12): 1781-1784.
DeShon H R, Schwartz S Y, Dorman L M, et al. 2004. Seismogenic zone structure along the Middle America Trench, Nicoya Peninsula, Costa Rica, from 3D local earthquake tomography[A]. In : 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Eberhart-Philips D. 1986. Three-dimensional velocity structure in north California Coast Ranges from inversion of local earthquakes arrival times[J]. Bull Seism Soc Amer, 76: 1 025-1 052.
Farra V, Madariaga R. 1988. Non-linear reflection tomography[J]. Geophys. J. Int., 95:135-147.
Geiger L. 1912. Probalibility method for the determination of earthquake epicenters from arrival time only[J]. Bull.St.Louis. Univ, 8:60-71.
Huang J, Zhao D. 2004. Crustal heterogeneity and seismotectonics of the region around Beijing[J]. Tectonophysics ,385:159-180.
Hole J A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography[J]. Journal of Geophysical Research, 97:6553-6562.
Humphreys E, Clayton R W. 1988. Adaptation of back projection tomography to seismic travel time problems[J]. Journal of Geophysical Research, 93:1073-1085.
Husen S Smith R B, Waite G P. 2004. Evidence for gas and magmatic sources beneath the Yellowstone volcanic field from seismic tomographic imaging[J]. J. Volcanol. Geotherm. Res. 131: 397-410.
Ivanson S. 1985. A study of methods for tomographic velocity estimation in the presence of low-velocity zones[J]. Geophysics, 50:969-988,.
Ivanson S. 1986. Some remarks concerning seismic reflection tomography and velocity ananlysis[J]. Geophys. J. R.Astr.Soc,87.539-557,.
Julian B R, Gubbins D. 1977. Three-dimensional sesmic ray tracing[J]. J. Geophys, 43:95-113.
Kamei R, Minami Y, Xue Z. 2004. Evaluation of Seismic Properties and Their Dependence on CO_2 Phases in a Porous Sandstone; Insights From Laboratory Injection Experiments[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Krajewski J, Roecker S, Thurber C. 2004. Combined Teleseismic and Local Earthquake Tomography of the SAFOD Drill Site[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Koch M. 1985. A numerical study on th e determination of the 3-D structure of th lithosphere by linear and non-liner inversion of teleseismic traveltimes[J]. Geophys. J.R.Astr.Soc,80:73~93.
Koch M. 1993a. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-I.evelopment,validation and optimal regularization of the method[J]. Geophys. J.Int., 112:385-412.
Koch M. 1993b. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-II.Application to the northern Rhine Grabern/Rhanish Massif Region Germany[J]. Geophys. J.Int., 112:413-428,.
Koch M. 1993c. Simultaneous inversion for 3-D crustal structure and hypocenters including direct,refracted and reflected phase-III.Application to the southern Rhine Grabern seismic regin, Germany[J]. Geophys. J.Int., 112:429-447,
Lees J M, Crosson R S. 1989. Tomographic inversion for three-dimensional velocity structure at Mount St. Helens using earthquake data[J]. Journal of Geophysical Research, 94: 5716-5728.
Lees J M, Wu H. 2000. Possion's ratio and porosity at Coso geothermal area, Californica[J]. J. Volcanol. Geotherm. Res., 95: 157-163.
Mishra O P, Zhao D. 2003. Crack density, saturation rate and porosity at the 2001 Bhuj, India, earthquake hypocenter; a fluid-driven earthquake[J]. Earth Plant. Sci. Letters, 212: 393-405.
O'Connell R J, Budiansky B. 1974. Seismic velocities in dry and saturated cracked solids[ J]. Journal of Geophysical Research,79: 5412-5426.
Pavlis G L, Booker J R. 1980. The Mixed Discrete-Continuous Inverse Problem: Application to the Simultaneous Determination of Earthquake Hypocenters and Velocity Structure[J]. Journal of Geophysical Research . 85(B9) 4801-4810.
Pickett G R. 1963. Acoustic character logs and their applications in formation evaluation[J]. J. Petrol. Tech, 15: 650~667.
Shedloc. 1987. Elastic wave velocity structure of the crust and upper mantle beneath the north China basin[J]. Journal of Geophysical Research, 92(B9):9327~9350.
Thurber C H. 1983. Earthquake Locations and Three-Dimensional Crustal structure in the Coyote Lake Area, Central California[J]. Journal of Geophysical Research, 88(B10): 8226-8236.
Um J, Thurber C. 1987. A fast algorithm for two-point seismic ray tracing[J]. Bull. Seism. Soc.Am., 77(3):972~986.
Vanorio T, Virieux J, and Lotorre D. 2004. Beyond Delay Time Tomography the Estimation of Rock physics properties from P- and S- Velocity models, In:2004 Fall Meeting, American Geophysical Union, Abs.[C], 2004.
Wagner L S, Zandt G, Beck S, et al. 2004. Black upper mantle beneath the Sierras Pampeanas, Argentina? Results from seismic tomography[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Wang Qing-Chen, Sun Shu, Li Ji-Liang et al. 1989. Tectonic evolution of QinLing mountains[J]. Scientia Geologica Sinica(in Chinese),2:129~141
Wiggins, R. A, 1972. The generalized linear inverse problem: Implication of surface waves and free oscillations for earth structure[J]. Rev. Geophys. Space Phys., 10:251-285.
Zhao D, Horiuchi S, Takgi A. 1990. 3-D seismic velocity structure of the crust and the uppermost mantle in the northeastern Japan arc, Tectonopgysics, 181:135-149.
Zhao D. 1991. A tomographic study of seismic velocity structure in the Japan Islands,Ph.D.Theis, Tohoku University.
Zhao D, Hasegawa A, Horiuchi S. 1992. Tomographic imaging of P and S wave velocity structure eneath northeastern Japan[J]. Journal of Geophysical Research, 97(B 13): 19909-19928.
Zhao D, Hasegawa A. 1993. P wave tomographic imaging of the crust and upper mantle beneath the Japan lslands[J] Journal of Geophysical Research, 98(B3):4333~4353.
Zhao D.1994. Deep structure of Japan subduction zone as deriveved from local, regional, and teleseismic events[J]. Journal of Geophysical Research, 99(B11):22313-22329.
Zhao D. 1998. The 1995 Kobe earthquake: Seismic image of the source zone and its implications for rupture nucleation[J]. Journal of Geophysical Research, 103(B5): 9967～9986.
Zhao D. 2001. Seismic structure and origin of hotspots and mantle plumes[J]. Earth. Planet. Sci. Lett., 192: 251～265.
Zhao D. 2004. Seismic imaging of North China: insight into intraplate volcanism and seismoteetonics[A]. In: 2004 Fall Meeting, American Geophysical Union, Abs.[C].
Zhao D. 2004. Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics[J]. Phys. Earth Planet. Int., 146: 3～34.
Zcng Rongsheng, Wang Chunyong, Zhang Dongning. 1995a. On the dynamics of extensional basin, PAGEOPH, 145(3/4): 579-603.
陈辉,刘福田,游庆瑜,等.1999.地震体波成像中的几个问题[J].地球物理学进展,14(4):9～17.
丁国瑜.1991.中国岩石圈动力学概论[M].北京.地震出版社.
丁志峰.1999.近震层析成像的理论及应用[学位论文].北京:中国地震局地球物理研究所.
丁志峰,曾融生.1994.用近震资料反演京津唐地区的地壳三维速度结构[J].华北地震科学,12(2):14～20.
丁志峰.2001.青藏高原地震波三维速度结构的研究[J].中国地震,17(2).
地质部航空物探大队904队.1957.华北平原南部及周围山区航空磁测结果报告.
方盛明,张先康,嘉世旭,等.2002.华北地区布格重力异常的多尺度分解特征与地震活动性[J].大地测量与地球动力学,22 (1):34～39.
付英祺,杨季楷.1987.地史学简明教程[M].北京:地质出版社.
金安蜀,刘福田,孙永智.1980.北京地区地壳上地幔的三维P波速度结构[J].地球物理学报,23(2):172～182.
嘉世旭,张先康,方盛明.2001.华北裂陷盆地不同块体地壳结构及演化研究[J].地学前缘,8 (2):259～265.
嘉世旭,张先康.2005.华北不同构造块体地壳结构及其对比研究[J].地球物理学报[J],48 (3):611～620.
国家地震局地球物理研究所.1978.近震分析[M].北京:地震出版社.
国家地震局《深部物探探测成果》编写组.1986.中国地壳上地幔地球物理探测成果[M].北京:地震出版社.
国家地震局地质研究所.1987.郯庐断裂[M].北京:地震出版社.
国家地震局鄂尔多斯周缘活动断裂系课题组.1988.鄂尔多斯周缘活动断裂系[M].北京:地震出版社.
郭增建,秦保燕.1979.震源物理[M].北京:地震出版社.
何正勤,丁志峰等.2001.中国大陆及其邻域地壳上地幔速度结构的面波层析成像研究[J].地球物理学报,(06):32～39.
李松林,张先康,张成科等.2002.玛沁-兰州-靖边地震测深剖面地壳速度结构的初步研究[J].地球物理学报,45(2):210～216.
李强,王椿镛,刘瑞丰,王溪莉,陈光英,陈光,李桂银.1999.应用层析成像技术研究华北地壳速度结构[J].地震地磁观测与研究.
刘瑞丰,陈培善,李强.1993.云南及其邻近地区三维速度图象[J].地震学报,15(1):61～67.
刘福田.1984.震源位置和速度结构的联合反演(Ⅰ)-理论和方法.地球物理学报[J],27(2):167～175.
刘建华,刘福田,孙若昧,等.1995.秦岭-大别造山带及其南北缘地震层析 成像,地球物理学报[J],38 (1):46～53.
刘福田,曲克信,吴华等.1986.华北地区的地震层析成像[J].地球物理学报,29(5):442～449.
刘福田,曲克信,吴华等.1989.中国大陆及其邻近地区的地震层析成象[J].地球物理学报[J],32(3):281～291.
刘福田.1991.三维速度结构的研究现状和展望.地球物理学报[J],34 (6):788～196.
刘建忠,李三忠,周立宏,高振平,郭晓玉.2004.华北板块东部中生代构造变形与盆地格局.海洋地质与第四纪地质[J],24(4):46～54.
刘伊克,常旭.2000.地震成析成像中解的定量评价及其应用.地球物理学报[J],43(2):252～256.
马杏垣,吴正文,谭应佳等.1979.华北地台基底构造[J].地质学报,53 (4):293～304.
马杏垣.1989.中国岩石圈动力学地图集.北京:中国地图出版社.
裴顺平等.2002.新疆及邻区Pn速度层析成像.地球物理学报[J],45(2):218-224.
孙若昧,刘福田.1995.京津地区地壳结构与强震的发生—Ⅰ:P波速度结构[J].地球物理学报,38(5):599～607.
孙若昧,赵燕来,吴丹.1996,京津唐地区地壳结构与强震的发生—Ⅱ:S波速度结构[J].地球物理学报,39(3):347～355.
孙若昧,赵燕来,梅世蓉.1993.勃海及其邻区的地震成析成像[J].地球物理学报,36(1):44～54.
宋仲和,陈国英,安昌强等.1993发中国大陆及其相邻海域瑞利波群速度分布特征[J].地震学报,15(1):32～38.
滕吉文,曾融生等.2002.东亚大陆及周边海域Moho界面深度分布和基本构 造格局[J].中国科学D辑,32(2):89～100.
滕吉文,张中杰,白武明等.2004.岩石圈物理学[M].北京:科学出版社.
王椿镛,王溪莉,颜其中.1993.昆明地震台网多事件定位问题的初步研究[J].地球物理学报,15 (2):136～145.
王椿镛,丁志峰等.1997.大别山造山带地壳S波速度结构[J].地球物理学报,40(3):337～345.
王椿镛,丁志峰等.1997.大别造山带地壳S波分裂和介质各向异性[J].科学通报.
王椿镛,张先康,丁志峰,阮红,邓宏钊.1997.大别造山带上部地壳结构的有限差分层析成像[J].地球物理学报,40(4):495～502.
王椿镛,W. D. Mooner,王溪莉,吴建平,楼海,王飞.2002.川滇地区地壳地幔三维速度结构研究[J].地震学报,24(1):1～16.
王椿镛,楼海,王飞.1999.大别山超高压变质带地壳结构及其构造意义[J].地震学报,21(5):534～544.
王夫运,张先康等.2002.用地震走时反演长白山天池火山地区的二维地壳结构[J].地震学报,24(2):144～152.
王夫运,张先康等.2004.张渤构造带中东段二维P波速度结构安新—香河—宽城剖面[J].地震学报,26增刊:31～41.
徐果明.2003.反演理论及其应用[M].北京:地震出版社.
徐锡伟,吴卫民,张先康等.2002.首都圈地区地壳最新构造变动与地震[M].北京.科学出版社.
姚振兴,李白基,梁尚鸿,朱培定,张立敏,卢善声.1981.青藏高原地区瑞利波群速度和地壳构造[J].地球物理学报,24(3):287～295.
杨景春.1993.中国地貌特征与演化[M].北京:海洋出版社.
杨文采.1997.地球物理反演的理论与方法[M].北京:地质出版社.