龙门山断裂带中段深地震宽角反射/折射剖面研究

李佳欣, 洪启宇, 郑需要

李佳欣, 洪启宇, 郑需要. 2017: 龙门山断裂带中段深地震宽角反射/折射剖面研究. 地震学报, 39(2): 188-206. DOI: 10.11939/jass.2017.02.003
引用本文: 李佳欣, 洪启宇, 郑需要. 2017: 龙门山断裂带中段深地震宽角反射/折射剖面研究. 地震学报, 39(2): 188-206. DOI: 10.11939/jass.2017.02.003
Li Jiaxin, Hong Qiyu, Zheng Xuyao. 2017: Study on the deep seismic wide-angle reflection/refraction profile in the middle of Longmenshan fault zone. Acta Seismologica Sinica, 39(2): 188-206. DOI: 10.11939/jass.2017.02.003
Citation: Li Jiaxin, Hong Qiyu, Zheng Xuyao. 2017: Study on the deep seismic wide-angle reflection/refraction profile in the middle of Longmenshan fault zone. Acta Seismologica Sinica, 39(2): 188-206. DOI: 10.11939/jass.2017.02.003

龙门山断裂带中段深地震宽角反射/折射剖面研究

基金项目: 

国家自然科学基金 40974033

国家自然科学基金 (41090292,40974033) 资助

国家自然科学基金 41090292

详细信息
    通讯作者:

    郑需要, e-mail: zhengxy@cea-igp.ac.cn

  • 中图分类号: P315.3+1

Study on the deep seismic wide-angle reflection/refraction profile in the middle of Longmenshan fault zone

  • 摘要: 射线追踪方法是研究地震波在横向非均匀地壳介质中传播的重要方法.本文推导了理论走时对网格化节点速度的偏导数公式,提出了针对深地震宽角反射/折射剖面数据反演的联合迭代法,并使用该方法对横跨龙门山断裂带中段的一条深地震宽角反射/折射剖面进行了反演和解释.首先,对每一炮的观测数据进行一维反演,在此基础上插值出一个粗略的二维速度模型;然后,使用射线追踪方法计算理论走时,再根据理论走时与观测走时的拟合程度对二维模型进行调整,以获得更加接近实际的二维速度模型;最后,利用联合迭代法对观测走时进行反演,经反复迭代使所有接收点理论走时与实测走时的残差平方和最小,最终获取该剖面的二维地壳速度结构.反演结果表明:测线东段的沉积盖层明显厚于中段褶皱带和西部高原,中部褶皱带部分地区出现基岩裸露;构造转换带两侧的地层分界面近于水平层状分布,其西侧的中、下地壳内各存在一个层间速度间断面;构造转换带内存在薄厚不等的低速层,自西向东有增厚趋势.此外,龙门山断裂带的3条主断裂向下深切结晶基底,这是由于西部松潘—甘孜地块自西向东运动,受到刚性扬子地块的阻挡,沿铲式断裂向上爬升所致;而在断层上盘距地表约15 km深处出现的最大剪应力极值区,正是发生汶川MS8.0地震的震源位置.
    Abstract: Ray tracing is an important method to study seismic wave propagating in the lateral inhomgeneous media of the crust. Partial derivatives of theoretical trival time with respect to the velocity at the meshing point of the model are derived. Simultaneous iterative reconstruction tomography (SIRT for short) method is used to invert the data of a deep seismic wide-angle reflection/refraction profle across the middle of Longmenshan fault zone, which can cover the shortage of time consuming and non-quantitative error analysis for trial-and-error method. Observed data for each shot are firstly inverted to get a one-dimensional velocity model, and then a simple two-dimensional velocity model is obtained by interpolating method. The two-dimensional model is adjusted by fitting observed and theoretical travel times, which makes the new model more close to the real one. SIRT method is further used to travel time inversion, and thus reaches final two-dimensional crustal velocity structure. The inverted result shows that the deposit in the east part of the profile is thicker than that in the middle folded zone and the western plateau, and there is basement rock appearing in the part area of the middle folded zone. The interfaces of the strata in the two sides of the transitional zone are almost horizontal, and in the western side of the transitional zone, there is a discontinuity interface in both the middle and lower crust. In addition, there are also lower velocity layers with different thickness in the transitional zone, and their thickness shows an increase tendency from east to west. It should be mentioned, that the three main faults in Longmenshan fault zone cut deeply the bottom of basement, which was resulted from the upward climbing movement of Songpan-Garze block to the eastern direction along a shovel-like fault and prevented from Yangtze block. There is a maximum shear stress region in the footwall of the fault inside the upper crust at the depth of 15 km, which is just the source location of the MS8.0 Wenchuan earthquake.
  • 龙门山断裂带位于四川盆地西部,青藏高原东缘,南起四川天全,总体以北东40°—50°方向延伸至广元,全长绵延约500 km,宽约30—40 km (李勇等,2009a)(图 1).整个龙门山断裂带主要由3条北东走向的断裂组成:汶川—茂县断裂 (后山断裂)、北川—映秀断裂 (中央断裂) 和江油—灌县断裂 (前山断裂); 3条断裂呈大致平行延伸,经过长期且多阶段的构造演化作用,该区地质构造也愈加复杂 (徐杰等,2010李志伟等,2011).作为稳定的扬子地块与活动的松潘—甘孜地块之间的过渡带,龙门山断裂带不仅在地形上显示出剧烈的变化,其两侧的地壳厚度 (含沉积盖层,下同) 也存在着相似的变化模式,即地壳厚度从松潘—甘孜地块之下的57—58 km减少至四川盆地之下的43—44 km,地形和地壳厚度的巨大差异表明该区发生了强烈的构造变形 (李志伟等,2011).

    图  1  龙门山区域主要活动构造及深地震测深剖面位置图
    F1:江油—灌县断裂;F2:北川—映秀断裂;F3:汶川—茂县断裂
    Figure  1.  Main active tectonic settings and location of the deep seismic sounding profile in Longmenshan region
    F1: Jiangyou-Guanxian fault; F2: Beichuan-Yingxiu fault; F3: Wenchuan-Maoxian fault

    据历史地震资料记载,中等强度及以上的地震主要集中在龙门山断裂带的南段及中南段,且龙门山断裂带的3条断裂上均已发生过M6以上的中等强度地震. 李勇等 (2009b)根据历史地震记录分析推测,龙门山断裂带近几十年来已进入地震高发期,3条主断裂均具备发生更大强度地震的潜在危险.

    龙门山断裂带长期以来一直是地球科学家们关注的焦点. 2008年5月12日,位于龙门山中央的映秀—北川断裂上的汶川县映秀镇发生MS8.0地震 (杨智娴等,2012).此次地震引起了世界范围内地球科学界的极大关注.在此之前,关于龙门山地区的地下地质构造,地学界已开展了大量的研究工作. 蒋福珍和方剑 (2001)使用重力场分离和密度反演方法对康滇地区的地壳构造进行了研究; 李勇等 (2005)给出了龙门山断裂带的均衡重力异常及其对青藏高原东缘山脉地壳隆升的约束; 王谦身等 (2009)研究了四川中西部地区的地壳结构与重力均衡; 李立和金国元 (1987)以及孙洁等 (2003)使用大地电磁测深方法得到了攀西裂谷带及龙门山断裂带地壳上地幔的电性结构;张先等 (1996)获得了四川盆地及其西部边缘震区的居里等温面;张培震等 (2008)徐锡伟等 (2008)对龙门山断裂带的构造演化过程以及地质成因进行了深入的探讨;刘建华等 (1989)Huang等 (2002)Wang等 (2003)雷建设等 (2009)利用大量天然地震走时资料确定了该区的三维速度结构;刘启元等 (2009)李志伟等 (2011)利用地震层析成像方法和远震接收函数反演了地下速度结构、地壳厚度分布、平均泊松比和波速比,并且通过接收函数的偏振分析研究了龙门山断裂带两侧地壳的各向异性差异 (楼海等,2008Zhang et al,2009Robert et al,2010郑勇等,2013);刘启元等 (2009)邓文泽等 (2014)通过对汶川地震余震序列的研究确定了龙门山断裂带中上地壳速度结构和断裂构造,从而探讨了该区域深部构造和汶川地震的形成原因 (李志伟等,2011).

    自20世纪80年代中期起,龙门山地区实施了一系列壳幔构造的探测计划,布设了多条人工地震剖面,仅横穿龙门山断裂带中段的测线就有两条.但是,由于当时观测系统不完善,地震仪器落后,一直未能得到翔实可靠的地下深部构造数据资料,难以对龙门山地区的地壳结构进行细致的研究.此外,李志伟等 (2011)关于该区天然地震的观测研究显示,低速层分布和形态分布等方面存在着不一致.而汶川地震的发生,推翻了之前对于龙门山断裂带的一些认知,提出了一系列新的科学问题.为了更加深入细致地研究该区的壳幔精细速度结构,中国地震局地球物理勘探中心布设了一条横穿龙门山断裂带的人工深地震测深剖面.

    本文拟对横穿龙门山断裂带的这条人工深地震测深剖面数据进行反演计算和解释.首先用每一炮的多个震相反演一维速度模型,利用多炮一维模型插值出一个较为粗略的二维初始速度模型;然后使用射线追踪法计算反射、折射波的理论到时,并且根据理论走时与观测走时的拟合程度 (试错法),调整不同层位的速度和深度,得到一个较为接近真实地壳上地幔顶部的二维速度模型.在试错法建立的速度模型基础上,本文提出一种对深地震测深剖面数据进行联合迭代反演 (simultaneous iterative reconstruction tomography,简写为SIRT) 的方法, 并使用该方法对横跨龙门山断裂带中段的这条深地震宽角反射/折射剖面进行反演,以期获取该剖面的二维地壳速度结构.

    2008年汶川MS8.0地震发生后,中国地震局地球物理勘探中心布设了一条横穿龙门山断裂带的人工深地震测深剖面.该剖面东起四川遂宁 (105°30′02″E,30°36′51″N),西至阿坝 (101°48′49″E,32°55′08″N),全长约500 km (图 1).由于整条测线横穿四川盆地、龙门山断裂带和川西北高原,因此对该测线的分析研究能够加深对该区域速度结构特征和深部构造的认识,推动对龙门山断裂带发展演化过程的研究.

    该测线横穿龙门山断裂带中段,途经绵竹、茂县等极震区,向西北经过黑水,最终到达阿坝州.整条测线布设12个炮点,炸药量为600—2800 kg (嘉世旭等,2014),采用井下爆破的方式,井深范围为15—61 m.龙门山断裂带及其附近的炮点较密集,最小炮检距约为10 km,而在远离龙门山断裂带的四川盆地和高原地区炮检距最大约为110 km.

    图 2给出了测线所采用的宽角反射/折射观测系统,通过该系统可以实现对地下介质进行多次覆盖和交叉采样,保证壳内波组可以连续对比追踪 (张恩会等,2013),从而压制部分干扰波,提高地震记录的信噪比.本文选取12炮中震相清晰、容易识别且较为典型的6炮数据进行处理分析,具体列于表 1.

    图  2  龙门山断裂带中段及两侧宽角反射/折射观测系统及检波器高程
    Figure  2.  Layout of wide-angle reflection/ refraction observatory system of the middle Longmenshan fault zone and elevation of receivers along the profile
    表  1  炮点参数一览表
    Table  1.  Shot parameters list
    炮点
    编号
    桩号
    /km
    炮点坐标 药量
    /kg
    平均井深
    /m
    高程
    /m
    地理位置
    东经/° 北纬/°
    SP1 113.478 105°30′02″ 30°36′51″ 1800 61 280 四川省遂宁市新桥镇黄连沱村
    SP5 255.298 104°18′26″ 31°23′14″ 800 20 486 四川省绵竹县绵远乡双狮村
    SP8 290.191 104°03′13″ 31°37′07″ 1000 10 1310 四川省绵竹县清平乡长河坝村
    SP9 320.638 103°48′10″ 31°47′26″ 600 26.5 1 712 四川省茂县沟口乡五家村
    SP12 467.784 102°24′14″ 32°24′44″ 2688 20.8 3596 四川省红原县龙日坝乡
    SP13 545.360 101°48′49″ 32°55′08″ 2814 42.5 3468 四川省阿坝县麦昆乡
    下载: 导出CSV 
    | 显示表格

    整条测线的原始数据分为12个组段,首先,从其中挑选出同一炮点的数据,重新分组形成共炮点数据文件;然后,对每一炮的地震数据进行分析,识别出清晰可靠的震相;最后,进行走时拾取,得到观测走时数据.

    识别震相时,不仅要对多炮的对应波组进行对比,而且要分析震相的运动学和动力学特征.震相识别是所有工作的基础,它不仅是构建初始速度模型的依据,同时还对反演的可靠性有决定作用.因此,对于不太清晰的或者无法确定的震相应采取宁缺毋滥的原则予以剔除.最终分析识别得到的主要震相为:结晶基底折射波Pg,莫霍面反射波PM,折射波Pn和地壳内反射波Pi(图 3).测线以龙门山断裂带为界,东侧为扬子地块 (四川盆地),西侧为松潘—甘孜地块 (川西北高原),两大地块之间存在着不同的地下地质构造和速度结构,因而使得两侧的地震波震相、走时和振幅也存在较大差异.

      3  深地震测深地震记录截面及走时拾取图 (折合速度为6.0 km/s)
    黑色短线表示拾取的走时. (a) SP1; (b) SP5; (c) SP8
      3.  Seismic record sections and picked travel times for the shots SP1 (a), SP5 (b) and SP8 (c)
    Reduced travel time is 6.0 km/s, and black short lines represent picked travel times
      3  深地震测深地震记录截面及走时拾取图 (折合速度为6.0 km/s)
    黑色短线表示拾取的走时. (d) SP9; (e) SP12; (f) SP13
      3.  Seismic record sections and picked travel times for the shots SP9 (d), SP12 (e) and SP13 (f)
    Reduced travel time is 6.0 km/s, and black short lines represent picked travel times

    Pg震相是地壳内结晶基底的折射波.测线东侧 (100—270 km桩号段) 追踪区间为10—100 km,以清晰可靠的初至波出现 (图 3).在炮点附近视速度迅速增加,远离炮点时视速度逐渐趋于稳定,大约为5.8—6.2 km/s,折合到时约为0—1.8 s,且其能量随着炮检距的增大而逐渐减弱.这些特征体现了扬子地块稳定的结构和较厚的沉积层.测线中部 (270—460 km桩号段) 位于龙门山断裂带,复杂的地下地质构造使得Pg波组的追踪范围大大缩小,仅在炮点附近可以追踪到,折合到时降至约0.1—0.4 s.测线西侧 (460—600 km桩号段) 的Pg波折合到时约为0.3—1.0 s,显示了较浅的沉积层厚度和高速的地下速度结构. Pg震相波组的上述特征是由地表沉积盖层厚度和地壳结晶层顶部介质速度结构所决定的,而局部折合到时的超前、滞后一般是由地表局部的隆起和凹陷构造所引起.

    作为壳幔分界面莫霍面的反射波,由于上下岩性差异较大,速度差异也较大,PM波在整个可追踪段震相清晰且振幅能量强,自炮检距80 km处出现,最远可达300 km,可连续对比追踪.与其它壳内反射波震相的比较表明,PM波的时距曲线更为陡峭,视速度更大.

    Pn波是来自上地幔的折射波,能够反映上地幔顶部的速度结构.由于该波组的能量较弱,连续对比追踪比较困难,除了SP8可以在炮点两侧追踪到外,其余炮仅能在单侧追踪. Pn波出现的最小距离为140 km,最远可达300 km,追踪长度最小仅有十几千米,显示了复杂的上地幔地质构造.以龙门山断裂带为界,测线西段川西北高原的折合到时在约190 km处为0;而在东段的四川盆地,其折合到时在约260 km处为0,显示了高速稳定的上地幔速度结构,视速度明显大于川西北高原.

    Pi波是地壳内各层的反射波组,其中P1和P3分别为上地壳底部C1界面和中地壳底部C3界面的反射波震相.在所选6炮地震记录截面 (图 3) 上可以看出,地壳内反射波Pi震相较弱;川西北高原可识别的波组震相比四川盆地明显要多,说明川西北高原的地下地质结构更为复杂.

    位于四川盆地的炮点SP1西侧和SP5东侧的震相显示,有3组反射波震相较明显,为P1,P3和P4波组.这3组反射波出现在炮检距为50—180 km范围内,随着炮检距的增大,折合到时逐渐减小,其时距曲线相对龙门山断裂带西侧较为陡峭,与四川盆地较为简单的层位分布和较高的地下介质速度相对应.位于龙门山断裂带的SP8和SP9炮点,其附近的反射波震相较凌乱,震相追踪距离较短,折合到时不是单纯地随炮检距增大,这反映了龙门山断裂带复杂的地下速度结构和层内界面展布.位于川西北高原的反射波组震相增多,主要的震相有P1,P2,P3,P4和P5,追踪范围为50—270 km.相较于四川盆地,其时距曲线较平缓,说明壳内介质速度比四川盆地要低.

    构建地下速度模型需要自上而下确定地壳各层的埋深和速度值.地表的速度值利用炮点附近的Pg波走时得到,界面深度及其上覆层介质的平均速度由x2-t2方法反演得到.在此基础上,构造每一层的速度梯度模型,然后通过走时拟合获取各炮点下方的地壳上地幔一维速度模型.

    龙门山断裂带构造变形剧烈.横跨该断裂带的地壳和上地幔速度结构存在强烈的横向不均匀性,因此使用简单的一维地壳模型不可能很好地描述真实的地下速度结构.从图 4可以看出,位于该区的SP5的右支 (图 4b),SP8的两支 (图 4c) 和SP9的左支 (图 4d) 均较难拟合,比较突出的是SP8左支和SP9左支的P1震相,其走时曲线出现凹凸弯曲,表明其下方存在比较大的界面起伏.此外,SP5和SP8炮点的PM震相理论走时与实测走时的残差较大,说明该区下方莫霍面起伏较大. SP1位于龙门山东侧的四川盆地中部,SP12和SP13位于川西北高原,一维拟合即可取得较好的结果,说明这两个地区的地下结构呈水平层状分布.在所有的实测走时中,会出现局部走时点的突跳,这与地下存在的不规则高速体、低速体有关,一维拟合无法处理,需要建立二维地壳速度结构模型进行研究.

      4  单炮基于一维速度模型的走时拟合图. (a) SP1; (b) SP5;(c) SP8
    小矩形框内曲线表示一维速度模型;实线代表理论走时,红点代表拾取走时
      4.  Computed travel-time data of single shot based on 1-D velocity model
    (a) SP1; (b) SP5; (c) SP8. The curves in the small frame represent the 1-D velocity model. Solid lines represent theory travel time, red dots represent observed travel time
      4  单炮基于一维速度模型的走时拟合图. (d) SP9; (e) SP12;(f) SP13
    小矩形框内曲线表示一维速度模型;实线代表理论走时,红点代表拾取走时
      4.  Computed travel-time data of single shot based on 1-D velocity model
    (d) SP9; (e) SP12; (f) SP13. The curves in the small frame represent the 1-D velocity model. Solid lines represent theory travel time, red dots represent observed travel time

    首先利用每一炮的一维速度模型,使用插值法构造出一个较为粗略的二维初始速度模型;然后应用“剥皮法”,由浅入深地对沉积层、上地壳、中地壳、下地壳和上地幔顶部,采用射线追踪法计算反射、折射波的理论到时 (Červený V, Pšenčík,1984);最后根据理论走时与观测走时的拟合程度,调整不同层位的速度和深度,即可得到一个较为接近真实地壳上地幔的二维速度模型.

    由于试错法存在耗时和不能给出定量误差的缺点,仅仅借助试错法,很难得到高精度的地壳速度结构.为此,本文将使用SIRT方法对观测走时进行反演.首先对二维初始模型进行网格化并建立射线路径上任意节点处的速度与周边网格节点速度之间的关系,然后推导出走时相对于网格节点速度的偏导数,最后建立走时反演方程组.

    首先使用矩形网格对二维速度模型进行网格化.每一层介质为一个独立的网格单元,单元的左、右边界与研究区域的左、右边界重合,网格的上、下边界将所在层位的上、下两个界面包含在内.每个层位的矩形网格尺度可以不一致,网格的大小取决于二维初始模型速度的纵向和横向变化,速度梯度较大时采用较密的网格,反之,则采用较稀疏的网格.

    对之前得到的二维初始速度模型进行横向和纵向线性插值,可以得到用于反演的二维网格化速度模型.每个网格有4个节点,节点速度用vij表示,其中i为沿x方向 (水平方向) 的节点序号,j为沿深度方向的节点序号.取其中一个矩形网格进行分析,网格的4个节点的速度分别为v11, v12, v21, v22(图 5).不同射线可能穿过网格的不同部位,射线在该网格内的长度也不同,如图 5中射线1和射线2.由于正演射线追踪得到的射线可能在某网格内的射线节点过少,常常会出现射线2的情况 (V1V2,…, Vn代表沿射线的节点),这可能导致反演时所求的方程组欠定.在这种情况下,需要对网格内的射线点进行线性插值,以保证最终反演的方程组为适定或超定.

    图  5  速度网格单元
    空心圈表示射线点,红色圆点表示插值点
    Figure  5.  Velocity grid cell
    Open circles represent ray points, and red dots represent interpolation points

    对网格内的射线节点采用反距离平方法进行线性插值.根据射线追踪结果可知,射线2在图 5所示的矩形网格内没有节点,它在网格外的两个节点分别为V1V2.为了得到射线2在矩形网格内任意一段射线两个节点 (图中任意两个相邻红色圆点, 例如AB) 的速度参数,首先要在V1V2的连线上确定射线段两节端点AB的坐标,然后分别求出每一个节点到矩形4个顶点的距离.设节点A到矩形4个顶点的距离为d1, d2, d3d4,节点B到矩形4个顶点的距离为e1, e2, e3e4,4个节点的速度依次为V11, V12, V21, V22(图 5),那么,如果用4个网格节点速度表示射线节点速度,AB两节点的速度为

    (1)

    (2)

    射线路径上任意相邻两点AB的走时差tB-tA对矩形节点处速度的偏导数为

    (3)

    式中

    其中s1为射线节点AB之间的路径长度.同样可得∂tB/∂vA∂tB/∂vB. ∂vA/∂vij∂vB/∂vij可由式 (1) 和 (2) 对vij求偏导获取.

    对于每一条射线,沿射线路径周围节点速度的改变会导致射线终点走时的改变.用Δt1, Δt2, …, Δtn分别表示第1条,第2条,…, 第n条射线的走时增量,可得

    (4)

    上式中每一条射线的走时是该条射线上所有射线段走时的总和.将方程组 (4) 等号左右两边互换,可用矩阵形式表示为

    (5)

    此即为用于反演的方程组,其右边为观测走时与理论走时的差,左边系数矩阵的每个元素是理论走时对模型节点速度的偏导数.在完成射线追踪后, 系数矩阵可由式 (1),(2) 和 (3) 求出;等号左边的矢量 (ΔV11, ΔV12, …, Δvmn)T是模型节点速度的改变量,可以通过反演计算求得.方程组 (5) 的系数矩阵是一个稀疏矩阵,只有那些有射线穿过的网格的相关节点的偏导数不为零,其余节点的偏导数均为零.

    大型稀疏矩阵方程组 (5) 可通过许多方法求解,例如,直接解法,雅克比共轭梯度法,不完全柯列斯基 (Cholesky) 共轭梯度法, 预条件共轭梯度法等.这些方法有的要求高内存,有的要求计算速度快,有的可能存在收敛和稳定性问题.本文将采用层析成像理论中的代数迭代法和联合迭代法 (SIRT).利用代数迭代法计算一条射线之后就可以将走时残差所对应的速度变化分配到射线路径上的相关节点;利用SIRT法则是计算完所有射线后将走时残差所对应的速度变化量分配到相关的节点上.由于SIRT法要求模型更新的次数较少,所以比较稳健.本文主要使用SIRT进行迭代反演.

    将走时拟合后得到的二维速度模型作为初始反演模型,使用上述反演方法对2668个走时数据 (见表 2) 进行反演,可以进一步缩小理论走时与观测走时之间的误差. Pg波组的到时记录用来约束地壳顶部速度结构,P1-P5反射波组用来约束对应的壳内C1-C5界面,PM和Pn分别用来约束莫霍面和上地幔顶部速度结构.迭代反演进行13次后,所有检波器实测走时与理论走时残差的平方和呈减小趋势,如图 6所示.可见,在反演的初期,残差平方和急剧减小,随后出现震荡,并最终趋于收敛,其中在第11次反演过程中残差平方和最小,为25.7257 s2.

    表  2  用于反演的全部震相及其走时数据统计
    Table  2.  Statistics of the phases and travel times used in inversion
    炮点 Pg P1 P2 P3 P4 P5 PM Pn
    SP1 31 27 - 26 - - 48 12
    SP5 136 83 63 35 56 - 72 6
    SP8 105 106 94 58 69 40 63 13
    SP9 155 101 89 73 84 37 - 13
    SP12 81 73 88 94 91 77 85 11
    SP13 65 36 26 55 33 55 61 10
    合计 573 426 360 347 359 209 329 65
    注:“-”表示无此震相.
    下载: 导出CSV 
    | 显示表格
    图  6  迭代反演的走时残差平方和与迭代次数变化关系曲线
    Figure  6.  Relationship between the sum of travel time residual squares and iteration number of inversion

    图 7给出了反演前与反演迭代第11次的对比结果,可以看出,速度模型有了较大的改善. SP1(图 7a) 和SP13(图 7f) 炮点模型反演前实测走时与理论走时的残差较小,拟合效果比较好;反演后速度模型变化不大. SP5(图 7b) 反演前左支走时残差较小,右支P2-PM震相走时残差较大; 反演后右支走时残差得到了压制. SP8(图 7c) 反演前右支走时残差较小,左支残差主要出现在P1和P2震相,理论走时偏小;反演后走时残差变小. SP9(图 7d) 反演前左支P1和P2震相及右支的P3-PM震相残差较大,左支理论走时偏小,右支理论走时偏大;反演后拟合效果变好. SP12(图 7e) 反演前左支P2震相理论走时偏小,P3-PM震相理论走时偏大, 反演后走时残差变小.

      7  基于二维速度模型的走时拟合图. (a) SP1; (b) SP5
    左图为反演迭代前结果,右图为反演迭代第11次时结果.黑色符号表示理论走时,红色符号表示实测走时,箭头表示炮点位置
      7.  Travel-times fitting based on 2-D velocity model. (a) SP1; (b) SP5
    Left panels represent the results befor inversion, right panels represent the results after 11 iteration inversions. Black symbols represent theory travel times, red symbols represent observed travel times, arrows represent shot location
      7  基于二维速度模型的走时拟合图. (c) SP8; (d) SP9; (e) SP12; (f) SP13
    左图为反演迭代前结果,右图为反演迭代第11次时结果.黑色符号表示理论走时,红色符号表示实测走时,箭头表示炮点位置
      7.  Travel-times fitting based on 2-D velocity model. (c) SP8; (d) SP9; (e) SP12; (f) SP13
    Left panels represent the results befor inversion, right panels represent the results after 11 iteration inversions. Black symbols represent theory travel times, red symbols represent observed travel times, arrows represent shot location

    通过反演得到了二维P波速度结构剖面图, 如图 8所示.可以看出,整条测线下方的地壳结构以龙门山为界,两侧的厚度和速度结构均存在较大的差异.测线东段 (100—270 km桩号段) 位于四川盆地,地壳厚度约为43—44 km;测线西段 (460—600 km桩号段) 位于川西北高原,地壳厚度约为57—58 km;测线中段 (270—460 km桩号段) 位于松潘—甘孜构造转换带,强烈的构造活动使得该区界面起伏明显并伴有错断,地壳厚度介于东、西两段之间.

    图  8  宽角反射/折射剖面二维P波速度结构图
    红色星形表示汶川地震震源在剖面上的投影,白色水平箭头表示块体运动方向F1:江油—灌县断裂;F2:北川—映秀断裂;F3:汶川—茂县断裂
    Figure  8.  2-D P wave velocity model along the wide-angle reflection/refraction profile The red star represents the projection of hypocenter of the Wenchuan earthquake
    The red star represents the projection of hypocenter of the Wenchuan earthquake on the profile, and the white horizontal arrows represent the direction of block motion F1: Jiangyou-Guanxian fault; F2: Beichuan-Yingxiu fault; F3: Wenchuan-Maoxian fault

    沿着测线自西向东,东、西两段沉积盖层的特征存在巨大差异.川西北高原沉积盖层的厚度为2.5—3.0 km,介质速度为3.60—5.25 km/s;沿测线向东经过龙日坝后,沉积盖层逐渐变薄;而从桩号280 km左右开始,盖层厚度急剧增大;当测线到达德阳市时,盖层厚度达到5.5—6.5 km,并延伸至遂宁,介质速度为3.55—5.34 km/s.龙门山断裂带两侧沉积盖层厚度的差异与扬子板块和松潘—甘孜板块的地质年龄相吻合.

    上地壳位于沉积盖层GC1界面之间,由基底回折波Pg和P1反射波组约束控制.川西北高原上地壳厚度为15—19 km,速度为5.94—6.16 km/s,在10 km深处存在约3 km厚的低速层;四川盆地上地壳厚度为12—15 km,速度为5.87—6.18 km/s.在松潘—甘孜构造转换带的东部,上地壳界面自西向东逐渐抬升,层厚变小,层内出现不规则低速块体,速度在横向上存在差异;其西部层厚逐渐加大,界面趋于水平,层上部存在高速块体.

    中地壳位于C1界面与C3界面之间,由反射波组P2和P3控制约束.西部高原地区中地壳厚度为12—14 km,速度为6.24—6.31 km/s;其间存在一个明显的速度间断面C2,速度由6.16 km/s突变为6.31 km/s,该间断面将中地壳分为两层,且均为高低速互层;龙门山褶皱带下方存在一长150 km厚40 km的低速体.四川盆地中地壳厚度为6—7 km,相较于川西北高原大幅度变薄,速度为6.29—6.40 km/s,其埋深起伏较大,在龙门山前陆盆地其埋深为18—25 km,在成都平原埋深为21—28 km.

    下地壳位于C3界面与莫霍面之间,由反射波组P4,P5和莫霍面反射波PM共同控制约束.川西北高原下地壳厚22—25 km,速度值偏低,约为6.46—6.91 km/s;四川盆地下地壳厚16—18 km,速度值较高,为6.66—7.20 km/s.下地壳中存在两个速度间断面:一个为C4间断面,该间断面横贯东西,并且其西部的速度跃变值为0.15 km/s,其东部的速度跃变值为0.25 km/s;另一个速度间断面C5仅存在于龙门山构造带西部之下,其速度从6.38 km/s跃变为6.61 km/s.

    莫霍面作为壳幔边界,其深度沿测线变化很大,由川西北高原的53—55 km,自西向东经龙门山褶皱带,逐渐抬升至四川盆地的43—44 km.上地幔顶部的速度Pn在横向上存在比较剧烈的变化,川西北高原速度约为7.95 km/s,而四川盆地高达8.10 km/s.

    将汶川地震的震源垂直投影至测量剖面上,如图 8红色星形所示,其在地面的投影坐标为 (103.7830°E,31.6879°N),震源深度为 (15.5±0.3) km,桩号为332.9 km,可见其位于上地壳底部的地层抬升区.

    本文研究的测线横跨西部松潘—甘孜地块、中部龙门山断裂带和东部扬子地块,在地壳厚度、壳内速度结构方面存在着巨大的差异,低速层、高速体的分布特征也有其独特的地质成因.测线剖面显示:沉积盖层在东部四川盆地的分布明显厚于中部褶皱带和西部高原,而中部褶皱带部分地区出现基岩裸露,表明扬子板块的地质年龄要老于松潘—甘孜板块,中部褶皱带处于活跃的地质构造期;地壳内构造转换带两侧的地层分界面近于水平层状分布,转换带西侧的中、下地壳内各存在一个层间速度间断面,而东侧则不存在类似的间断面.对应上、中、下地壳的分界面,西侧比东侧的速度要小,而且自上而下逐渐增大,下地壳底部莫霍面的速度西侧比东侧要低0.3 km/s.在厚度方面,西侧比东侧的上、中、下地壳厚度分别厚3—7 km,6—8 km,6—9 km,可见地壳的增厚主要由被改造的中下地壳介质形变增厚所致 (李志伟等,2011嘉世旭等,2014),这有别于蔡学林等 (1999)一文中地壳厚度基本不变的结果.在构造转换带内,存在着薄厚不等的低速层,自西向东有增厚的趋势.川西北高原的松潘—甘孜地块在向东挤压刚性的扬子地块时受到阻挡,使得转换带内的低速层堆积、增厚,中下地壳物质向上抬升,形成了上陡下缓、向北西倾斜的构造,从而使测线东、西两侧的莫霍面产生了较大的高差,地壳厚度由四川盆地的43—44 km增厚至川西北高原的57—58 km.

    龙门山断裂带的3条主断裂即前山断裂F1、中央断裂F2和后山断裂F3,分别向下深切结晶地壳,而断裂F1F3正是自上而下并向西侧倾斜的巨型铲式低速通道的包络 (图 8).这是由于西部松潘—甘孜地块自西向东运动,受到刚性扬子地块的阻挡,沿铲式断裂向上爬升,便在断层上盘距地表 15 km左右深处出现最大剪切应力的极值区,最终应变能以地震的形式得以释放,与杨智娴等 (2012)对汶川MW7.9地震震源位置的测定结果相吻合. 陶玮等 (2011)利用有限元方法对该区的铲式断裂进行模拟,其结果也显示最大剪应力变化的极值区位于深约14—18 km地壳内.

    龙门山断裂带由南西向北东延伸约500 km,不同地段的深部构造存在着差异.以往的深地震测深测线主要布测于龙门山南段,且观测系统布设不够完善.横穿龙门山断裂带的这条测线完善了龙门山深地震测深资料,改善了观测系统的布设,从而得到了精细的地下速度结构,同时为其它地质和地球物理研究提供了参考和依据.

    本文提出的反演方法是在试错法建立的速度模型的基础上,对试错法的一种补充.试错法在试错的后期残差收敛较慢,新的反演方法大大减少了工作量和时间;试错法无法对走时残差进行量化,本文使用的方法则弥补了试错法的缺点,使误差处于可控的范围内.然而,该方法尚有很多需要改进的地方,例如网格内射线点速度的插值方法,可以考虑使用更多的网格节点进行插值;对射线上的点进行样条插值,以期更精确地表示射线追踪得到的射线.

    中国地震局地球物理勘探中心为本文提供了人工地震数据,参与该项目野外施工的工作人员克服困难,取得了高质量的野外地震资料,中国地震局地球物理勘探中心嘉世旭研究员和刘志研究员为初期的资料处理给予了有益的指导帮助,审稿专家提出了诸多改进意见,作者在此一并表示感谢.

  • 图  1   龙门山区域主要活动构造及深地震测深剖面位置图

    F1:江油—灌县断裂;F2:北川—映秀断裂;F3:汶川—茂县断裂

    Figure  1.   Main active tectonic settings and location of the deep seismic sounding profile in Longmenshan region

    F1: Jiangyou-Guanxian fault; F2: Beichuan-Yingxiu fault; F3: Wenchuan-Maoxian fault

    图  2   龙门山断裂带中段及两侧宽角反射/折射观测系统及检波器高程

    Figure  2.   Layout of wide-angle reflection/ refraction observatory system of the middle Longmenshan fault zone and elevation of receivers along the profile

    3   深地震测深地震记录截面及走时拾取图 (折合速度为6.0 km/s)

    黑色短线表示拾取的走时. (a) SP1; (b) SP5; (c) SP8

    3.   Seismic record sections and picked travel times for the shots SP1 (a), SP5 (b) and SP8 (c)

    Reduced travel time is 6.0 km/s, and black short lines represent picked travel times

    3   深地震测深地震记录截面及走时拾取图 (折合速度为6.0 km/s)

    黑色短线表示拾取的走时. (d) SP9; (e) SP12; (f) SP13

    3.   Seismic record sections and picked travel times for the shots SP9 (d), SP12 (e) and SP13 (f)

    Reduced travel time is 6.0 km/s, and black short lines represent picked travel times

    4   单炮基于一维速度模型的走时拟合图. (a) SP1; (b) SP5;(c) SP8

    小矩形框内曲线表示一维速度模型;实线代表理论走时,红点代表拾取走时

    4.   Computed travel-time data of single shot based on 1-D velocity model

    (a) SP1; (b) SP5; (c) SP8. The curves in the small frame represent the 1-D velocity model. Solid lines represent theory travel time, red dots represent observed travel time

    4   单炮基于一维速度模型的走时拟合图. (d) SP9; (e) SP12;(f) SP13

    小矩形框内曲线表示一维速度模型;实线代表理论走时,红点代表拾取走时

    4.   Computed travel-time data of single shot based on 1-D velocity model

    (d) SP9; (e) SP12; (f) SP13. The curves in the small frame represent the 1-D velocity model. Solid lines represent theory travel time, red dots represent observed travel time

    图  5   速度网格单元

    空心圈表示射线点,红色圆点表示插值点

    Figure  5.   Velocity grid cell

    Open circles represent ray points, and red dots represent interpolation points

    图  6   迭代反演的走时残差平方和与迭代次数变化关系曲线

    Figure  6.   Relationship between the sum of travel time residual squares and iteration number of inversion

    7   基于二维速度模型的走时拟合图. (a) SP1; (b) SP5

    左图为反演迭代前结果,右图为反演迭代第11次时结果.黑色符号表示理论走时,红色符号表示实测走时,箭头表示炮点位置

    7.   Travel-times fitting based on 2-D velocity model. (a) SP1; (b) SP5

    Left panels represent the results befor inversion, right panels represent the results after 11 iteration inversions. Black symbols represent theory travel times, red symbols represent observed travel times, arrows represent shot location

    7   基于二维速度模型的走时拟合图. (c) SP8; (d) SP9; (e) SP12; (f) SP13

    左图为反演迭代前结果,右图为反演迭代第11次时结果.黑色符号表示理论走时,红色符号表示实测走时,箭头表示炮点位置

    7.   Travel-times fitting based on 2-D velocity model. (c) SP8; (d) SP9; (e) SP12; (f) SP13

    Left panels represent the results befor inversion, right panels represent the results after 11 iteration inversions. Black symbols represent theory travel times, red symbols represent observed travel times, arrows represent shot location

    图  8   宽角反射/折射剖面二维P波速度结构图

    红色星形表示汶川地震震源在剖面上的投影,白色水平箭头表示块体运动方向F1:江油—灌县断裂;F2:北川—映秀断裂;F3:汶川—茂县断裂

    Figure  8.   2-D P wave velocity model along the wide-angle reflection/refraction profile The red star represents the projection of hypocenter of the Wenchuan earthquake

    The red star represents the projection of hypocenter of the Wenchuan earthquake on the profile, and the white horizontal arrows represent the direction of block motion F1: Jiangyou-Guanxian fault; F2: Beichuan-Yingxiu fault; F3: Wenchuan-Maoxian fault

    表  1   炮点参数一览表

    Table  1   Shot parameters list

    炮点
    编号
    桩号
    /km
    炮点坐标 药量
    /kg
    平均井深
    /m
    高程
    /m
    地理位置
    东经/° 北纬/°
    SP1 113.478 105°30′02″ 30°36′51″ 1800 61 280 四川省遂宁市新桥镇黄连沱村
    SP5 255.298 104°18′26″ 31°23′14″ 800 20 486 四川省绵竹县绵远乡双狮村
    SP8 290.191 104°03′13″ 31°37′07″ 1000 10 1310 四川省绵竹县清平乡长河坝村
    SP9 320.638 103°48′10″ 31°47′26″ 600 26.5 1 712 四川省茂县沟口乡五家村
    SP12 467.784 102°24′14″ 32°24′44″ 2688 20.8 3596 四川省红原县龙日坝乡
    SP13 545.360 101°48′49″ 32°55′08″ 2814 42.5 3468 四川省阿坝县麦昆乡
    下载: 导出CSV

    表  2   用于反演的全部震相及其走时数据统计

    Table  2   Statistics of the phases and travel times used in inversion

    炮点 Pg P1 P2 P3 P4 P5 PM Pn
    SP1 31 27 - 26 - - 48 12
    SP5 136 83 63 35 56 - 72 6
    SP8 105 106 94 58 69 40 63 13
    SP9 155 101 89 73 84 37 - 13
    SP12 81 73 88 94 91 77 85 11
    SP13 65 36 26 55 33 55 61 10
    合计 573 426 360 347 359 209 329 65
    注:“-”表示无此震相.
    下载: 导出CSV
  • 蔡学林, 曹家敏, 刘援朝, 魏显贵. 1999.青藏高原多向碰撞—楔入隆升地球动力学模式[J].地学前缘, 6(3): 181-189. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY199903024.htm

    Cai X L, Cao J M, Liu Y C. Wei X G. 1999. Geodynamic models of multidirectional collision-wedging uplift of the Qinghai-Tibet Plateau[J]. Earth Science Frontiers, 6(3) 181-189 (in Chinese).

    邓文泽, 陈九辉, 郭飚, 刘启元, 李顺成, 李昱, 尹昕忠, 齐少华. 2014.龙门山断裂带精细速度结构的双差层析成像研究[J].地球物理学报, 57(4): 1101-1110. doi: 10.6038/cjg20140408

    Deng W Z, Chen J H, Guo B, Liu Q Y, Li S C, Li Y, Yin X Z, Qi S H. 2014. Fine velocity structure of the Longmenshan fault zone by double-difference tomography[J]. Chinese Journal of Geophysics, 57(4): 1101-1110 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract10289.shtml

    嘉世旭, 刘保金, 徐朝繁, 刘志, 酆少英, 张建狮, 林吉焱, 田晓峰, 刘巧霞, 郭文斌. 2014.龙门山中段及两侧地壳结构与汶川地震构造[J].中国科学:地球科学, 44(3): 497-509. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201403010.htm

    Jia S X, Liu B J, Xu Z F, Liu Z, Feng S Y, Zhang J S, Lin J Y, Tian X F, Liu Q X, Guo W B. 2014. The crustal structures of the central Longmenshan along and its margins as related to the seismotectonics of the 2008 Wenchuan earthquake[J]. Science China Earth Sciences, 57(4): 777-790. doi: 10.1007/s11430-013-4744-9

    蒋福珍, 方剑. 2001.康滇地区重力场分离、密度反演与地壳构造[J].地震学报, 23(4): 391-397. http://www.dzxb.org/Magazine/Show?id=26963

    Jiang F Z, Fang J. 2001. Gravity field separation, density inversion and crustal tectonics in Kang-Dian region[J]. Acta Seismologica Sinica, 23(4): 391-397 (in Chinese). doi: 10.1007%2Fs11589-001-0119-y.pdf

    雷建设, 赵大鹏, 苏金蓉, 张光伟, 李凤. 2009.龙门山断裂带地壳精细结构与汶川地震发震机理[J].地球物理学报, 52(2): 339-345. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200902006.htm

    Lei J S, Zhao D P, Su J R, Zhang G W, Li F. 2009. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake[J]. Chinese Journal of Geophysics, 52(2): 339-345 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract913.shtml

    李立, 金国元. 1987.攀西裂谷带及龙门山断裂带地壳上地幔的大地电磁测深研究[J].物探与化探, 11(3): 161-169. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH198703000.htm

    Li L, Jin G Y. 1987. Telluric electromagnetic sounding study of crust and upper mantle in the Panxi "rift zone" and the Longmenshan faulted zone[J]. Geophysical and Geochemical Exploration, 11(3): 161-169 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-WTYH198703000.htm

    李勇, 徐公达, 周荣军, Densmore A L, Ellis M A. 2005.龙门山均衡重力异常及其对青藏高原东缘山脉地壳隆升的约束[J].地质通报, 24(12): 1162-1168. doi: 10.3969/j.issn.1671-2552.2005.12.012

    Li Y, Xu G D, Zhou R J, Densmore A L, Ellis M A. 2005. Isostatic gravity anomalies in the Longmen mountains and their constraints on the crustal uplift below the mountains on the eastern margin of the Qinghai-Tibet Plateau[J]. Geological Bulletin of China, 24(12): 1162-1168 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZQYD200512013.htm

    李勇, 黄润秋, Densmore A L, 周荣军, 曹叔尤. 2009a.汶川8.0级地震的基本特征及其研究进展[J].四川大学学报:工程科学版, 41(3): 7-25. http://www.cnki.com.cn/Article/CJFDTOTAL-SCLH200903006.htm

    Li Y, Huang R Q, Densmore A L, Zhou R J, Cao S Y. 2009a. Basic features and research progresses of Wenchuan MS8.0 earthquake[J]. Journal of Sichuan University: Engineering Science Edition, 41(3): 7-25 (in Chinese). https://www.researchgate.net/publication/286044531_Basic_features_and_research_progresses_of_Wenchuan_Ms_80_earthquake

    李勇, 黄润秋, 周荣军, Densmore A L, Ellis M A, 闫亮, 董顺利, Richardson N, 张毅, 何玉林, 陈浩, 乔宝成, 马博琳. 2009b.龙门山地震带的地质背景与汶川地震的地表破裂[J].工程地质学报, 17(1): 3-18. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ200901003.htm

    Li Y, Huang R Q, Zhou R J, Densmore A L, Ellis M A, Yan L, Dong S L, Richardson N, Zhang Y, He Y L, Chen H, Qiao B C, Ma B L. 2009b. Geological background of Longmen shan seismic belt and surface ruptures in Wenchuan earthquake[J]. Journal of Engineering Geology, 17(1): 3-18 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-GCDZ200901003.htm

    李志伟, 胥颐, 黄润秋, 郝天珧, 徐亚, 刘劲松, 刘建华. 2011.龙门山地区的P波速度结构与汶川地震的深部构造特征[J].中国科学:地球科学, 41(3): 283-290. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201103001.htm

    Li Z W, Xu Y, Huang R Q, Hao T Y, Xu Y, Liu J S, Liu J H. 2011. Crustal P-wave velocity structure of the Longmenshan region and its tectonic implications for the 2008 Wenchuan earthquake[J]. Science China Earth Sciences, 54(9): 1386-1393. doi: 10.1007/s11430-011-4177-2

    刘建华, 刘福田, 吴华, 李强, 胡戈. 1989.中国南北带地壳和上地幔的三维速度图象[J].地球物理学报, 32(2): 143-152. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX198902003.htm

    Liu J H, Liu F T, Wu H, Li Q, Hu G. 1989. Three dimensional velocity images of the crust and upper mantle beneath North-South Zone in China[J]. Chinese Journal of Geophysics, 32(2): 143-152 (in Chinese). https://www.researchgate.net/profile/Jiankun_He/publication/228373448_Three-Dimensional_Velocity_Images_of_the_Crust_and_Upper_Mantle_beneath_the_North-South_Zone_in_China/links/0deec53363d7c9dc3f000000.pdf

    刘启元, 李昱, 陈九辉, 郭飚, 李顺成, 王峻, 张绪奇, 齐少华. 2009.汶川MS8.0地震:地壳上地幔S波速度结构的初步研究[J].地球物理学报, 52(2): 309-319. http://www.cqvip.com/QK/94718X/2009002/29656012.html

    Liu Q Y, Li Y, Chen J H, Guo B, Li S C, Wang J, Zhang X Q, Qi S H. 2009. Wenchuan MS8.0 earthquake: Preliminary study of the S-wave velocity structure of the crust and upper mantle[J]. Chinese Journal of Geophysics, 52(2): 309-319 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWX200902003.htm

    楼海, 王椿镛, 吕智勇, 姚志祥, 戴仕贵, 尤惠川. 2008. 2008年汶川MS8.0级地震的深部构造环境:远震P波接收函数和布格重力异常的联合解释[J].中国科学: D辑, 38(10): 1207-1220. http://www.cnki.com.cn/Article/CJFDTotal-JDXK200810005.htm

    Lou H, Wang C Y, Lü Z Y, Yao Z X, Dai S G, You H C. 2009. Deep tectonic setting of the 2008 Wenchuan MS8.0 earthquake in southwestern China: Joint analysis of teleseismic P-wave receiver functions and Bouguer gravity anomalies[J]. Science in China: Series D, 52(2): 166-179. doi: 10.1007/s11430-009-0009-z

    孙洁, 晋光文, 白登海, 王立凤. 2003.青藏高原东缘地壳、上地幔电性结构探测及其构造意义[J].中国科学: D辑, 33(增刊): 173-180. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGEM200209001122.htm

    Sun J, Jin G W, Bai D H, Wang L F. 2003. Sounding of electrical structure of the crust and upper mantle along the eastern border of Qinghai-Tibet Plateau and its tectonic significance[J]. Science in China: Series D, 46(S2): 243-253. doi: 10.1360/03dz0019

    陶玮, 胡才博, 万永革, 沈正康, 王康. 2011.铲形逆冲断层地震破裂动力学模型及其在汶川地震研究中的启示[J].地球物理学报, 54(5): 1260-1269. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201105017.htm

    Tao W, Hu C B, Wan Y G, Shen Z K, Wang K. 2011. Dynamic modeling of thrust earthquake on listric fault and its inference to study of Wenchuan earthquake[J]. Chinese Journal of Geophysics, 54(5): 1260-1269 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract7959.shtml

    王谦身, 滕吉文, 张永谦, 张雪梅, 杨辉. 2009.四川中西部地区地壳结构与重力均衡[J].地球物理学报, 52(2): 579-583. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200902031.htm

    Wang Q S, Teng J W, Zhang Y Q, Zhang X M, Yang H. 2009. The crustal structure and gravity isostasy in the middle western Sichuan area[J]. Chinese Journal of Geophysics, 52(2): 579-583 (in Chinese). https://www.researchgate.net/publication/286887745_The_crustal_structure_and_gravity_isostasy_in_the_middle_western_Sichuan_area

    徐杰, 高祥林, 周本刚, 计凤桔, 张进, 白玉柱, 陈国光. 2010. 2008年汶川8.0级地震的发震构造:沿龙门山断裂带新生的地壳深部断裂[J].地学前缘, 17(5): 117-127. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201005012.htm

    Xu J, Gao X L, Zhou B G, Ji F J, Zhang J, Bai Y Z, Chen G G. 2010. The seismogenic structure of the 2008 Wenchuan MS8.0 earthquake: A newly generated deep fault in crust along the Longmen mountains fault zone in the shallow subsurface[J]. Earth Science Frontiers, 17(5): 117-127 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DXQY201005012.htm

    徐锡伟, 闻学泽, 叶建青, 马保起, 陈杰, 周荣军, 何宏林, 田勤俭, 何玉林, 王志才, 孙昭民, 冯希杰, 于贵华, 陈立春, 陈桂华, 于慎鄂, 冉勇康, 李细光, 李陈侠, 安艳芬. 2008.汶川MS8.0地震地表破裂带及其发震构造[J].地震地质, 30(3): 597-629. http://www.wenkuxiazai.com/doc/7b7719d14028915f804dc264.html

    Xu X W, Wen X Z, Ye J Q, Ma B Q, Chen J, Zhou R J, He H L, Tian Q J, He Y L, Wang Z C, Sun Z M, Feng X J, Yu G H, Chen L C, Chen G H, Yu S E, Ran Y K, Li X G, Li C X, An Y F. 2008. The MS8.0 Wenchuan earthquake surface ruptures and its seismogenic structure[J]. Seismology and Geology, 30(3): 597-629 (in Chinese). http://www.dz-dz.com.cn/EN/abstract/abstract8419.shtml

    杨智娴, 陈运泰, 苏金蓉, 陈天长, 吴朋. 2012. 2008年5月12日汶川MW7.9地震的震源位置与发震时刻[J].地震学报, 34(2): 127-136. http://www.dzxb.org/Magazine/Show?id=28738

    Yang Z X, Chen Y T, Su J R, Chen T C, Wu P. 2012. The hypocenter and origin time of the MW7.9 Wenchuan earthquake of May 12, 2008[J]. Acta Seismologica Sinica, 34(2): 127-136 (in Chinese). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZXB201202001.htm

    张恩会, 楼海, 嘉世旭, 李永华. 2013.云南西部地壳深部结构特征[J].地球物理学报, 56(6): 1915-1927. doi: 10.6038/cjg20130614

    Zhang E H, Lou H, Jia S X, Li Y H. 2013. The deep crust structure characteristics beneath western Yunnan[J]. Chinese Journal of Geophysics, 56(6): 1915-1927 (in Chinese). http://manu39.magtech.com.cn/Geophy/EN/abstract/abstract9576.shtml

    张培震, 徐锡伟, 闻学泽, 冉勇康. 2008. 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因[J].地球物理学报, 51(4): 1066-1073. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200804017.htm

    Zhang P Z, Xu X W, Wen X Z, Ran Y K. 2008. Slip rates and recurrence intervals of the Longmen Shan active fault zone, and tectonic implications for the mechanism of the May 12 Wenchuan earthquake, 2008, Sichuan, China[J]. Chinese Journal of Geophysics, 51(4): 1066-1073 (in Chinese). https://www.researchgate.net/publication/279545184_Slip_rates_and_recurrence_intervals_of_the_Longmen_Shan_active_fault_zone_and_tectonic_implications_for_the_mechanism_of_the_May_12_Wenchuan_earthquake_2008_Sichuan_China

    张先, 虎喜凤, 沈京秀, 田长征, 赵丽, 刘敏, 陈秀文. 1996.四川盆地及其西部边缘震区居里等温面的研究[J].地震学报, 18(1): 83-88. http://www.dzxb.org/Magazine/Show?id=27644

    Zhang X, Hu X F, Shen J X, Tian C Z, Zhao L, Liu M, Chen X W. 1996. Study of Curie isothermal surface in Sichuan basin and the seismic area on its western margin[J]. Acta Seismologica Sinica, 18(1): 83-88 (in Chinese). doi: 10.1007/BF02650629

    郑勇, 葛粲, 谢祖军, Yang Y J, 熊熊, 许厚泽. 2013.芦山与汶川地震震区地壳上地幔结构及深部孕震环境[J].中国科学:地球科学, 43(6): 1027-1037. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201306011.htm

    Zheng Y, Ge C, Xie Z J, Yang Y J, Xiong X, Hsu H T. 2013. Crustal and upper mantle structure and the deep seismogenic environment in the source regions of the Lushan earthquake and the Wenchuan earthquake[J]. Science China Earth Sciences, 56(7): 1158-1168. doi: 10.1007/s11430-013-4641-2

    Červený V, Pšenčík I. 1984. SEIS83: Numerical modelling of seismic wave fields in 2-D laterally varying layered structure by the ray method[G]// Documentation of Earthquake Algorithms. Boulder Colo: World Data Center for Solid Earth Geophysics, SE-35: 36-40.

    Huang J L, Zhao D P, Zheng S H. 2002. Lithospheric structure and its relationship to seismic and volcanic activity in southwest China[J]. J Geophys Res, 107(B10): ESE 13-1-ESE 13-14. https://www.researchgate.net/publication/259450893_Lithospheric_structure_and_its_relationship_to_seismic_and_volcanic_activity_in_southwest_China

    Robert A, Zhu J, Vergne J, Cattin R, Chan L S, Wittlinger G, Herquel G, de Sigoyer J, Pubellier M, Zhu L D. 2010. Crustal structures in the area of the 2008 Sichuan earthquake from seismologic and gravimetric data[J]. Tectonophysics, 491(1/2/3/4): 205-210. http://www.academia.edu/8282440/Crustal_structures_in_the_area_of_the_2008_Sichuan_earthquake_from_seismologic_and_gravimetric_data

    Wang C Y, Chan W W, Mooney W D. 2003. Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications[J]. J Geophys Res, 108(B9): 2442. https://www.researchgate.net/profile/Walter_Mooney/publication/236271821_Three-dimensional_velocity_structure_of_Crust_and_Upper_Mantle_in_southwestern_China_and_its_tectonic_implications/links/00b4951780a3f50bc1000000.pdf

    Zhang Z J, Wang Y H, Chen Y, Houseman G A, Tian X B, Wang E, Teng J W. 2009. Crustal structure across Longmenshan fault belt from passive source seismic profiling[J]. Geophys Res Lett, 36(17): L17310. doi: 10.1029/2009GL039580

  • 期刊类型引用(2)

    1. Yingkang LI,Jianwei GAO,Jian HAN,Yu'e YANG. Geophysical evidence for thrusting of crustal materials from orogenic belts over both sides of the Yangtze Block and its geological significance. Science China(Earth Sciences). 2019(05): 812-831 . 必应学术
    2. 李英康,高建伟,韩健,杨予鄂. 扬子块体两侧造山带地壳推覆的地球物理证据及其地质意义. 中国科学:地球科学. 2019(04): 687-705 . 百度学术

    其他类型引用(1)

图(11)  /  表(2)
计量
  • 文章访问数:  717
  • HTML全文浏览量:  440
  • PDF下载量:  30
  • 被引次数: 3
出版历程
  • 收稿日期:  2016-05-30
  • 修回日期:  2016-06-19
  • 发布日期:  2017-02-28

目录

/

返回文章
返回