利用宽角反射/折射地震剖面揭示芦山MS7.0地震震区深部孕震环境

王帅军， 王夫运， 张建狮， 刘宝峰， 张成科， 赵金仁， 段玉玲，宋向辉， 邓晓果， 马策军， 孙一男， 臧怡然， 李怡青

中国地震局地球物理勘探中心， 郑州 450002



利用宽角反射/折射地震剖面揭示芦山MS7.0地震震区深部孕震环境

王帅军， 王夫运*， 张建狮， 刘宝峰， 张成科， 赵金仁， 段玉玲，宋向辉， 邓晓果， 马策军， 孙一男， 臧怡然， 李怡青

中国地震局地球物理勘探中心， 郑州 450002

2013年4月芦山地震发生后，中国地震局迅速成立了芦山地震科学考察指挥部，要求查明芦山地震的深部构造环境和孕震背景.为此，中国地震局地球物理勘探中心于2013年9月至11月在芦山震源区布设了一条长约410 km的人工地震高分辨宽角反射/折射探测剖面，获得了信噪比较高的人工地震探测数据，采用地震射线走时正演拟合构建了该区的地壳及上地幔二维P波速度结构模型，结果显示：扬子块体和松潘—甘孜块体显示出迥异的速度结构特征，地壳厚度由南向北逐渐加厚.沉积盖层在四川盆地厚达7.8 km，而进入松潘—甘孜块体沉积层最薄处只有几百米厚，几乎出露地表；在中上地壳，扬子块体平均速度比松潘—甘孜块体高0.2 km·s-1，在盆地与高原耦合部位(构造转换带)以北深度大约20 km左右有一厚度为8.0 km的软弱层(低速层)，该层内的速度为5.80 km·s-1，明显低于周围介质的平均速度6.00～6.10 km·s-1；构造转换带内，震相显示紊乱、不清晰、不能连续对比，由地表至上地幔顶部壳内界面不连续、速度结构异常紊乱且呈现低速异常特征；在中下地壳，沿剖面速度呈现正梯度垂向增大变化；壳内界面在扬子块体内部起伏变化不大，但在构造转换带以北呈现急速加深的趋势，特别是Moho界面起伏变化较为明显，界面深度在距离50 km范围内由扬子块体的36.2 km迅速变化至松潘—甘孜块体下方的45.8 km，形成一陡变带.芦山MS7.0级地震震源位置位于二维速度结构异常紊乱和界面起伏变化的地带，研究表明，壳内界面及速度结构差异、起伏变化的特征与该区域的地震活动性关系密切.

芦山地震； 人工地震剖面； 速度结构； 壳内低速层； 构造转换带

According to the 2D crust-mantle P wave velocity structure model and previous research results of the study area, this research has obtained the following accomplishments: (1) The 2D crustal velocity structure has obvious zoning characteristics in the transverse direction, and the velocity structure shows low velocity anomaly in the tectonic transformation zone where the velocity structure varies significantly on both sides. The average velocity of the Yangtze block is obviously 0.10 km·s-1higher than that of Songpan-Garzê block, and the crustal thickness increases gradually from south to north, which presents rapid changes with the discontinuous interface near the tectonic transition zone. (2) The basement along the profile in Songpan-Garzê block is thinner, almost exposed to the surface near Xiaojin, and the thickest part is approximately 2.5 km near Yongfu County. The basement of Sichuan Basin is relatively thicker, about 7.8 km; clear changes of the upper crustal velocity structure are located between C1 and C2 under the Songpan-Garzê block. There is a negative gradient velocity structure, whose velocity varies from 5.80 to 5.90 km·s-1, 0.10～0.20 km·s-1lower than the average velocity around it. The depth range is from 13 to 21 km and the thickness is about 8 km. (3) The Pn wave of the upper mantle in the profile is difficult to identify, only seen when the SP7 shot point reflects velocity structure of the upper mantle beneath the Yangtze block. There is no Pn wave phase in other shot points. The velocity structure of the upper mantle is clearly different beneath the two blocks. It is 8.00～8.05 km·s-1beneath the Yangtze block while it is 7.95～7.90 km·s-1beneath the Songpan-Garzê block. (4) The 2D crust-mantle velocity model results show that along the profile the 2D contoured velocity is disordered with obvious interface undulation changes in the vicinity of the profile distance 150 km, 260 km, 280 km and 300 km. This phenomenon exactly matches with positions of Longquanshan fault, Jiangyou-Dujiangyan fault, Maowen-Tianquan fault and Xiaojin arc fracture in the geological structure map.

In the light of the projection results of aftershocks in this section, the main shock and aftershocks are mainly distributed in the high-low speed transition and negative velocity gradient areas. From the perspective of abnormal thickness in the negative gradient velocity structure, this layer appears to act as a detachment surface in the process of eastward movement of material or the stress on the Qinghai Tibet plateau, which leads to local stress accumulation around the block edge and tectonic activity. In terms of the geological tectonic background of LushanMS7.0 earthquake, the regional tectonic structure of modern earthquakes is controlled by nearly NS local stress field and EW regional stress field, and 2D model shows that the velocity structure has a great reversal in the vicinity of Lushan earthquake epicenter (i.e. about 13～14 km), which presents deep medium conditions leading to strong earthquakes. Therefore, LushanMS7.0 earthquake may be a sudden release of the stress on the southwest section in the Longmenshan fault zone.Keywords Lushan earthquake；Artificial earthquake profile；Velocity structure；Crustal low-speed layer；Tectonic transformation zone

1 引言

2013年4月20日，四川省雅安市芦山县境内发生MS7.0级强烈地震，主震及余震精定位的震源深度为13 km左右(赵博等，2013；吕坚等，2013)，震中位于龙门山断裂带南段，2008年5月12日的汶川MS8.0级地震也发生在该条断裂带上，两次地震的震中相距仅为80多公里(陈运泰等，2013；徐锡伟等，2013；高原等2013)．地震发生后，国内外学者对该地震进行了空间定位、震源机制和震源破裂过程的研究(王卫民等，2013；刘成利等，2013；张勇等，2013).结果表明，芦山地震是发生在龙门山断裂带南端的一次逆冲型地震，这次地震在震源性质上与2008年汶川地震同为逆冲破裂(Zhang et al.，2009a).由于先后两次地震均发生在龙门山推覆构造带上，在时空展布、发震构造和震源机制等方面存在一定的关联性，因此，关于芦山地震是否为汶川地震强余震引起了很大的争议(徐锡伟等，2013；单斌等，2013；杜方等，2013；颜照坤等，2014).

随着时间的推移，新资料不断涌现，尤其是活动断裂研究、地球物理探测、GPS观测等方面新成果不断出现，这些成果对于进一步深入认识该区的构造变形、地震活动背景等都有着重要影响.但是，现有的研究由于受到资料本身的研究精度和技术手段的限制，难以精确地给出该区的壳幔速度结构模型(王椿镛等，2003a；徐涛等，2014).因此，通过跨青藏高原东南缘和四川盆地布设了一条金川—芦山—乐山人工地震剖面.利用多种计算方法及综合研究工作，重新构建该区的高精度的壳幔速度结构模型，不仅有助于进一步正确理解青藏高原的最新地壳变形过程及机制，科学评价该区未来的大震活动趋势，还可为区域国土合理规划与开发利用，以及重大工程建设的地壳稳定性评价提供更为可靠的基础数据(王运生等，2013).

芦山地震发生后，亟需全面了解青藏高原东缘地区大震活动的新构造与活动构造背景及成因问题(陈立春等，2013).在此背景下中国地震局迅速成立了科学考察指挥部，由中国地震局地球物理研究所牵头，设立了多个科考专题任务组，要求从多个方面查明地震发生的背景，其中人工地震探测任务组要完成一条长约410 km、通过芦山震区的宽角反射/折射剖面，获取穿过震区的高分辨二维地壳速度结构和发震区精细结构，探明芦山地震的发震构造、孕震环境及深浅部构造关系，构建芦山地震发震构造模型，探讨芦山地震深部动力学背景，分析和对比龙门山断裂带不同区段的深部结构和构造差异特征，研究芦山地震与汶川地震的发震构造关系，对芦山震区的地震活动趋势和地震危险性判定提供深浅部构造依据，并对地球物理其他方法(楼海等，2008)的研究带来更为精细的地壳速度结构模型作为参考、相互验证与补充(丁志峰等，2008；Bai et al.，2011； Robert et al.，2010)，这一研究不仅有助于深化认识地震破裂机制，对地震间相互触发和空间迁移或未来地震发展趋势等也有十分重要的科学意义.

2 区域地质构造背景

研究区位于青藏高原的东南缘，是青藏高原现今地壳形变和地震活动最强烈的地区之一，在大地构造上是扬子地块、松潘—甘孜地块等多个微陆块体汇聚的地带，新生代期间经历了多期岩浆活动和强烈的构造变形作用，形成了众多规模与性质不同的断裂构造(Bai et al.，2010；Wang et al.，2009).进入新构造期间，特别是第四纪或晚第四纪期间，在印度板块与欧亚板块持续碰撞的动力学背景下，该区的地壳运动十分活跃，发育有众多规模、类型、活动性等各不相同的活动断裂(Lei and Zhao，2009；徐锡伟等，2013；刘保金等，2009；Zhang et al.，2010)，龙门山断裂带位于中国大陆中部南北地震带的中段,其东南侧为扬子准地台,西侧是松潘一甘孜地块.印支运动尤其是喜山期运动以来,伴随青藏高原的强烈抬升,发生了不同程度的大陆壳的相互碰撞、挤压,在龙门山前形成了一条以逆冲为主,走滑为辅的大型断裂带(石玉涛等，2009；孙长青等，2011；雷建设等，2009；Chen et al.，2013)，龙门山断裂带总体走向N45°E，倾向NW.自北西向南东，龙门山逆冲推覆构造带由3条西倾的主干断裂带组成，即龙门山后山断裂、中央断裂及前山断裂.每条断裂又分别由几个不同的段落组成.地质考察结果表明(石玉涛等，2009；孙长青等，2011)，这几条主断裂带自晚第四纪以来均显示由北西向南东的逆冲运动，并伴有显著的右旋走滑分量.在芦山地震发生的龙门山断裂带南段，叠瓦状次级逆断层发育.自西向东，排列着盐井—五龙、双石—大川、新开店、大邑、浦江—新津等几条NE向的断裂.其中，盐井—五龙断裂是主中央断裂的南段，双石—大川断裂属于前山断裂的南段，为龙门山构造带东南侧的边界断裂带.大邑断裂则是四川盆地内部的地表倾角较陡的隐伏断层.该区的新构造变形过程及其动力学机制、现今的活动构造变形样式、与大震活动关系等都是国内外地球物理学科领域极为关心的问题(Lei et al.，2008；Zhang et al.，2010).因此，在该区建立高精度的壳幔速度结构模型对于理解青藏高原东南缘地区现今地壳变形方式及其动力学机制、地震发生的背景都具有重要的意义(Xu et al.，2010；Wang et al.，2007).前人围绕上述问题，在该地区进行了大量的地震地质、地球物理等方面的工作(王椿镛等，2003a； Royden et al.，1997；Unsworth et al.，2005；Wei et al.，2001；Zheng et al.，2010；Zhang et al.，2009b；Yao et al.，2009)，这些工作对本研究提供了重要的参考.

研究区域地质构造十分复杂，多组深大断裂纵横交错.近年来，在该区域发生的多个中等强度的地震研究表明，地震的发生与该区的深大断裂带、复杂的壳幔过渡带、中上地壳高低速相间介质环境等密切相关，鉴于该区特殊的地震构造环境，一些研究者在该地区已经开展了大量的地球物理研究工作，为该区的深部物理场特性、地震孕育环境等提供了重要的深部资料，王椿镛等(2010)根据龙门山及其周边多个台站的宽频带远震记录，使用H-k叠加方法计算了该区的地壳厚度和波速比.结果表明该区域的地壳厚度总体变化是：从东向西增加，其中横跨龙门山断裂带的地壳厚度变化较大，松潘-甘孜地体北部和西秦岭造山带具有低泊松比，松潘—甘孜块体南部的上地壳物质向东运动，受刚性强度较大的扬子地台的阻挡，导致沿龙门山断裂带产生应变积累.当断层被地壳流体弱化，积累的应变能量快速释放，产生汶川MS8.0地震；赵博等(2013)利用双差定位法对芦山地震主震及余震序列进行重新定位.结果表明，地震深度主要集中在10～20 km．余震序列沿近似NE方向分布，长约35 km.沿近NW剖面，余震大多分布在大川—双石断裂的SE方向．随着流动台站震相数据和余震数据的增加，重新定位得到的余震分布特征与断层构造产状相似；郑勇等(2013)借助全国地震台网连续波形数据,使用背景噪声层析成像方法和远震接收函数分析方法,获得了震区及其周边地区精细的S波速度结构和地壳厚度、泊松比分布情况,研究结果发现，龙门山断裂带两侧剪切波速度和地壳厚度有非常显著的差异;芦山和汶川地震均位于地壳厚度和波速结构变化剧烈之处,断层的破裂面和余震的分布均处于地震波横向速度梯度和地壳厚度的横向梯度跳变的地区;地震深度处于从均匀波速结构向非均匀波速结构过渡区域.为了更好地理解该地区的深部地球物理场特征，2013年，我们横跨芦山MS7.0震源区实施了一条金川—芦山—乐山人工源宽角反射/折射地震探测剖面，为该区域的地震精确定位、地震孕育机制等提供深部速度模型.

3 剖面位置及观测系统

研究剖面东南起宜宾市王场镇附近(其坐标为104°39′15″E，29°05′40″N，剖面桩号50 km)，由东南向西北依次经自贡市的古文镇、留佳镇、乐山市的土主镇、汉阳镇、仁美镇、张场镇、雅安市的车岭镇、新店镇、龙门镇、宝兴、永富、小金、名山县北、芦山县的龙门、宝兴、陇东、永富、小金、金川等地，剖面的终点位于金川县太阳河乡且斯都村附近(其坐标为101°46′30″E，31°41′25″N，剖面桩号460 km)，全长410 km(图1).物探中心于2013年9月至11月完成了该条人工地震测深剖面的野外打井、炮点激发、观测点安置及数据采集等工作.

该地震探测剖面穿过了不同的地形条件和地理环境，跨越了多条断裂构造带和不同的地质构造单元，沿剖面地形条件和地理环境比较复杂，基本上全部为山区地形，剖面的东南端大约100 km区段位于四川盆地.在280 km桩号以北测线基本上位于川西高原，海拔较高，一般都在2500～3500 m，局部地段达到4000 m以上；在280 km桩号以南地区基本上位于扬子块体，海拔一般在1000～2000 m，观测点基本上都是沿山路或者乡村小路布设，观测点背景噪声水平较低.

在长约410 km的剖面上总共布设了268台三分向PDS-2型地震仪器，观测点距采取了疏密相间的布设方案，观测点距在0.8～2.0 km之间变化，在断裂带及特殊构造部位附近点距适当加密为0.8 km，使得沿测线变化剧烈的地段，保证浅部结构能够很好地成像.在川西高原及四川盆地内部观测点距为2.0 km.根据地震波场信息对地壳与上地幔深、浅结构的有效反映和主要断裂构造带的控制，确定炮点位置.沿剖面共布置了8个炮点，各个炮点的具体参数见表1，单炮药量的大小根据研究的目标、炮点所在位置周围环境情况、激发条件的不同和观测区间的大小而定，最大药量为3.5 t，最小药量为1.0 t，沿剖面8个炮点的激发和268台PDS-2型地震仪器的共同接收构成了较为完善的观测系统.

图1 剖面位置及地质构造图Fig.1 Location of the seismic sounding profile and geological structures

表1 炮点参数表Table 1 Table of shots parameters

图2 SP1炮点折合时间剖面图Fig.2 Converted time sections of shot (SP1)

图3 SP7炮点折合时间剖面图Fig.3 Converted time sections of shot (SP7)

4 震相分析

对获得的野外探测数据用折合速度6.00 km·s-1进行了折合(图2—3)，(图2—3中的纵坐标为折合到时，记录截面的横坐标为炮检距，右侧负炮检距对应测线的南东方向，横坐标正炮检距对应测线的北西方向)并对记录截面进行了2～12 Hz滤波，滤波后有效信号的信噪比大大提高，各组震相清晰、连续、能够可靠追踪与对比.根据获得的8炮记录截面的各组地震波震相特征进行了综合分析、研究，主要对比出了7组震相，它们分别为结晶基底以上强梯度带内回折波Pg、壳内反射波P1、P2、P3、P4、莫霍面强反射波Pm以及上地幔顶部回折波Pn.

通过对上述7组震相的走时拾取，利用WH、PLUCH等计算方法(Michel and Hirn，1980；国家地震局科技监测司，1988)求取了各波组的视速度、平均速度以及沿测线两个不同构造地质单元的大致地壳厚度及平均深度，这些参数为二维壳幔速度结构模型的构建提供了重要的参考和依据.

5 二维壳幔速度结构模型构建

在获得一维壳幔速度结构模型的基础上，依据有限差分反演方法(段永红等，2002；徐涛等，2014)和WH反演得出的各炮点附近地表速度值，确定剖面结晶基底以上的浅部速度分布，壳内各层界面沿剖面起伏变化和层内速度变化可参照一维地壳速度结构模型参数，再综合考虑用其他方法所计算的结果以及研究区内已有的钻孔资料、浅部地震探测成果、地质资料、其他地球物理探测资料(徐锡伟等，2013；李德威等，2013；谢祖军等2013，郑勇等，2013；房立华等，2013)和实际地震观测资料的特征，综合设计出剖面的初始二维地壳结构模型(Wang et al.，2014；王帅军等，2007).利用改进后的深地震测深资料处理Seis程序包，对该剖面获得的8炮观测资料进行了二维非均匀介质中动力学射线追踪(Çerveny et al.,1997；Çerveny and Hron,1980；Çerveny and Pençík，1984；Çerveny,1984)、 走时拟合的计算(徐涛等，2004；Xu et al.，2006，2010，2014；王夫运和张先康，2004)，经过反复正演计算，最终使理论走时、各波组的振幅特征与实测资料的记录特征达到了最佳的拟合，如图4.与此同时，获得了沿剖面的地壳二维P波速度结构和构造模型(图5).结果表明，在松潘—甘孜块体与四川盆地速度结构差异特征显著，我们把金川—芦山—乐山地震测深剖面所揭示的地壳结构分为上地壳、中地壳和下地壳三层.

上地壳：C2界面以上的部分为上部地壳，上地壳结构又分为三个亚层，G界面以上为第一亚层，G界面在四川盆地沉积盖层较厚，厚度在3.0～7.8 km之间变化.在松潘—甘孜块体盖层较薄，厚度在1.0～2.5 km之间变化.速度随深度的增加迅速加大，由近地表的4.70 km·s-1变化至结晶基底底部的5.75 km·s-1，可以看出，G界面以上为一强速度梯度层.第二亚层为G界面与C1界面之间的层位，C1界面在四川盆地内部缺失，上部地壳结构显示双层结构特征(Clark and Royden，2000).构造转换带以北逐渐出现C1界面.在松潘—甘孜块体上部地壳结构显示三层结构.C1界面深度大约为10 km左右，C1界面之上为一弱的梯度层， 之下为一低速层，其速度有0.10～0.20 km·s-1的负跳跃，速度在5.80～5.90 km·s-1之间变化，在两块体耦合位置地区下方速度等值线紊乱，反映了该区为强、弱速度梯度层的复杂结构，横向变化强烈.第三亚层为C1界面与C2界面之间的层位，C2界面沿剖面由南向北逐渐加深，在两块体交界位置附近界面不连续，两侧速度结构差异显著，且向北有明显的陡变加深，C2界面深度在16.0～19.8 km之间变化.

图4 Sp1和Sp7炮的射线追踪(b)与走时拟合(a)图Fig.4 Synthetic seismograms, ray tracing (b) and travel time fitting curves (a) for shots SP1 and SP7

图5 金川—芦山—乐山地震探测剖面二维速度结构图Fig.5 2-D crustal section along Jinchuan-Lushan-Leshan seismic profileF1: Longquan mountain fault; F2：Jiangyou—Dujiangyan fault; F3：Maowen—Tianquan fault; F4：Xiaojin arcuate fault.

中地壳：C3界面与C2界面之间为中地壳，C3界面的形态与C2界面的起伏变化较为一致，层间速度在6.15～6.40 km·s-1之间变化，在四川盆地该层速度高于松潘—甘孜块体0.05～0.10 km·s-1，由上至下呈现随深度加深而增加的正梯度变化，速度变化相对较为缓和，只有在两块体交汇处显示出速度结构紊乱等特征，C3界面在四川盆地埋深大约在26～28 km之间变化，在松潘—甘孜块体界面埋深大约在29～34 km之间变化，在两块体耦合位置附近界面显示不连续，由南向北呈现急剧加深趋势.

下地壳：C3界面之下与M界面之间的层位为下地壳，在松潘—甘孜块体内部C3至M界面之间有一C4界面，该界面深度大约在42～48 km之间变化，该层位之间速度变化相对较小，由C3界面下方的6.40 km·s-1变化至C4界面上方的6.50 km·s-1，在构造转换带下方速度结构有局部低速块体存在，低于周围介质速度0.05～0.10 km·s-1，该界面在构造转换带部位同样显示界面不连续和陡变特征，在松潘—甘孜块体下方下地壳显示双层结构特征，C4界面进入四川盆地逐渐消失，下地壳在四川盆地下方显示单层结构特征.从获得的8炮记录截面上几乎都能看到较强能量和大振幅的PmP反射震相.从获得的二维剖面可以看出Moho界面变化较为剧烈，该界面埋置深度在四川盆地约为42 km，在青藏高原下方最深处达到了62 km，特别是在扬子与松潘—甘孜块体耦合部位Moho界面显示迅速加深8 km，陡变幅度较壳内其他界面大，在构造转换带内从上至下反射波显示能量弱、波形紊乱、界面不连续，也暗示该区域为应力汇集区和两大块体碰撞的脆弱带，芦山MS7.0地震震源位置刚好位于构造转换带、高、低速结构变化剧烈区域.研究表明(王琼和高原，2014)强震往往发生在高低速变化剧烈、结构陡变地带内.Moho界面之上的速度在松潘—甘孜块体内速度值为6.60～6.70 km·s-1，在扬子块体内速度为6.70～6.80 km·s-1.沿剖面Moho界面显示出南浅北深的特征，Moho界面起伏变化特征与基底呈镜像关系.

总体来看，在四川盆地与松潘—甘孜块体的过渡带附近，显示出了速度结构紊乱、界面起伏变化明显的壳幔速度结构特征，这可能与该区段特殊的接触耦合关系及复杂的地质运动构造背景关系极为密切(滕吉文等，2014；高锐等，2004；吴建平等2001).

6 结论与讨论

(1) 从计算获得的二维壳幔速度结构模型可以看出，地壳速度结构在横向上存在明显的分区特征，构造转换带内速度结构显示低速异常，构造转换带两侧速度结构差异显著，扬子块体平均速度明显高于松潘—甘孜块体0.10 km·s-1左右，地壳厚度由南向北逐渐加厚，其中在构造转换带附近为界面急剧变化带，且界面不连续.

(2) 沿剖面基底在松潘—甘孜块体厚度较薄，在小金附近基底几乎出露地表，永富县附近最厚处约为2.5 km左右，四川盆地基底相对较厚，最厚约为7.8 km；上地壳速度结构变化较明显的区域位于松潘—甘孜块体下方C2与C1之间.有一负梯度速度结构体存在，速度变化在5.80～5.90 km·s-1之间，低于周围平均速度的0.10～0.20 km·s-1；深度范围为13～21 km，厚度大约8 km左右，赵博等(2013)利用芦山余震双差定位结果显示，芦山余震震源深度大多位于10～20 km之间，对应本剖面震源深度刚好位于剖面的高速和低速交接处.本文利用芦山地震122个3级以上的芦山地震余震在本剖面上的投影结果来看，主震及余震主要集中分布在本剖面的构造转换带内(图5).研究表明,中等强度的地震往往发生在低速与高速相间的位置(王琼和高原，2014).从负梯度速度结构异常厚度来看，该层在青藏高原物质或者应力东向运移过程中似乎充当了一个滑脱面，并导致周边块体边缘及活动构造局部应力不断积累.

(3) 上地幔Pn波在整条剖面上难以识别对比(Liang et al.，2004)，只在SP7炮点反映扬子块体下方的上地幔速度结构时显示出来，其他炮点均无Pn波震相表现，上地幔速度结构沿剖面两个块体下方变化较为明显，在扬子块体下方速度为8.00～8.05 km·s-1，而在松潘—甘孜块体下方的速度为7.90～7.95 km·s-1.扬子块体上地幔顶部速度高于松潘—甘孜块体，从金川—芦山—乐山地震探测剖面的结果显示四川盆地的刚性特征强于松潘—甘孜块体的尺度至少延伸到上地幔上部.暗示松潘—甘孜块体下方上地幔物质活动性较扬子块体下方明显活跃的特征.

(4) 研究区二维壳幔速度模型结果显示，沿剖面分别在150、260、280 km及300 km桩号附近二维速度等值线较为紊乱、壳内界面起伏变化明显，这一现象刚好与地质结构图上的龙泉山断裂、江油—都江堰断裂、茂汶—天全断裂及小金弧形断裂位置相互吻合.断裂的具体深度范围、倾角及倾向等有待进一步的工作研究.该区是我国强震多发地区之一，从芦山7.0级地震发生的地质构造背景看，该地区的现代地震构造受控于近南北向的区域构造应力场和东西向的局部应力场，二维模型结果显示在芦山地震震中位置附近大约在13.0～14.0 km左右速度结构发生较大的逆转，具备了强烈地震发生的深部介质条件，因此，芦山MS7.0地震的发生可能是龙门山断裂带西南段上应力的一次突然释放.

致谢 高信噪比的地震数据是后续解释工作的基础和保证.物探中心野外工作人员在数据采集过程中克服了种种困难，三位匿名审稿专家对本文提出了具有建设性的意见和建议，在此一并表示感谢.另外，谨以此文向中国地震局地球物理勘探中心成立60周年献礼！

Bai D H, Unsworth M J, Meju M A, et al. 2010. Crustal deformation of the eastern Tibetan Plateau revealed by magnetotelluric imaging.NatureGeoscience, 3(5): 358-362.Bai Z M, Tian X B, Tian Y. 2011. Upper mantle P-wave tomography across the Longmenshan fault belt from passive-source seismic observations along Aba-Longquanshan profile.J.AsianEarthSci., 40(4): 873-882.

Çerveny V, Molotkov I A, Pençík L. 1977. Ray Method in Seismology.Prague:University of Karlova.

Çerveny V, Hron F. 1980. The ray series method and dynamic ray tracing system for three-dimensional inhomogeneous media.BulletinoftheSeismologicalSocietyofAmerica, 70(1): 47-77.

Çerveny V. 1984. SEIS83-Numerical modeling of seismic wave fields in 2-D laterally varying layered structures by the ray method. Documentation of Earthquake Algorithms.

Chen Y, Zhang Z J, Sun C Q, et al. 2013. Crustal anisotropy from Moho converted Ps wave splitting analysis and geodynamic implications beneath the eastern margin of Tibet and surrounding regions.GondwanaResearch, 24(3-4): 946-957.

Chen Y T, Yang Z X, Zhang Y, et al. 2013. From 2008 Wenchuan earthquake to 2013 Lushan earthquake.ScientiaSinicaTerrae(in Chinese), 43(6): 1064-1072.

China Earthquake Administration Technology Monitoring Division Compiles. 1988. Deep Tectonics Research and Progress of China Mainland. Beijing: Geology Pess, 18-46.

Clark M K, Royden L H. 2000. Topographic ooze: Building the eastern margin of Tibet by lower crustal flow.Geology, 28(8): 703-706.

Ding Z F, Wu Y, Wang H, et al. 2008. Variations of shear wave splitting in the 2008 Wenchuan earthquake region.ScienceinChinaSeriesD:EarthSciences, 51(12): 1712-1716.

Du F, Long F, Ruan X, et al. 2013. TheMS7.0 Lushan earthquake and the relationship with theM8.0 Wenchuan earthquake in Sichuan, China.ChineseJ.Geophys.(in Chinese), 56(5): 1772-1783,doi:10.6038/cjg2013535. Duan Y H, Zhang X K, Fang S M. 2002. Three-dimensional finite-difference tomography of velocity structure of the upper crustal in North China.ChineseJ.Geophys. (in Chinese), 45(3): 362-369.

Fang L H, Wu J P, Wang W L, et al. 2013. Relocation of the mainshock and aftershock sequences ofMS7.0 Sichuan Lushan earthquake.Chin.Sci.Bull., 58(28-29): 3451-3459.

Gao R, Dong S W, He R Z, et al. 2004. Subduction process of the Yangtze continental block from Moho reflection image, south China.EarthScienceFrontiers(in Chinese), 11(3): 43-49.

Gao Y, Wang Q, Zhao B, et al. 2014. A rupture blank zone in middle south part of Longmenshan Faults: Effect after LushanMS7.0 earthquake of 20 April 2013 in Sichuan, China.ScienceChina:EarthSciences, 57(9): 2036-2044, doi: 10.1007/s11430-013-4646-x.

Lei J S, Zhao D P. 2009. Structural heterogeneity of the Longmenshan fault zone and the mechanism of the 2008 Wenchuan earthquake (MS8.0).Geochem.Geophys.Geosyst., 10(10), doi: 10.1029/2009GC002590.

Lei J S, Zhao D P, Su J R, et al. 2009. Fine seismic structure under the Longmenshan fault zone and the mechanism of the large Wenchuan earthquake.ChineseJ.Geophys. (in Chinese), 52(2): 339-345.

Li D W, Chen G F, Chen J L, et al. 2013. Earthquake prediction: From Lushan earthquake to continental earthquakes.EarthScienceFrontiers(in Chinese), 20(3): 1-10.

Liang CT, Song X D, Huang J L. 2004. Tomographic inversion of Pn travel times in China.J.Geophys.Res., 109: B11304.

Liu B J, Zhang X K, Feng S Y, et al. 2009. High-resolution seismic reflection profile across Pengzhou buried fault in the frontal areas of Longmen Shan.ChineseJ.Geophys. (in Chinese), 52(2): 538-546. Liu C L, Zheng Y, Ge C, et al. 2013. Rupture process of theMS7.0 Lushan earthquake,2013.ScienceChina:EarthSciences, 56(7):1187-1192, doi: 10.1007/s11430-013-4639-9.

Lou H, Wang C Y, Lü Z Y, et al. 2009. Deep tectonic setting of the 2008 WenchuanMS8.0 earthquake in southwestern China—Joint analysis of teleseismic P-wave receiver functions and Bouguer gravity anomalies.ScienceinChina(SeriesD):EarthScience,52(2):166-179.

Lü J, Wang X S, Su J R, et al. 2013. Hypocentral location and source mechanism of theMS7.0 Lushan earthquake sequence.ChineseJ.Geophys. (in Chinese), 56(5): 1753-1763,doi:10.6038/cjg20130533.

Michel B, Hirn A. 1980. Velocity-depth estimation from wide angle seismic reflection arrivals.AnnalsdeGeophysique, 36(1): 107-118.

Robert A, Pubellier M, de Sigoyer J, et al. 2010. Structural and thermal characters of the Longmen Shan (Sichuan, China).Tectonophysics, 491(1-4): 165-173.

Royden L H, Burchfiel B C, King R W, et al. 1997. Surface deformation and lower crustal flow in Eastern Tibet.Science, 276(5313): 788-790.

Shan B, Xiong X, Zheng Y, et al. 2013. Stress changes on major faults caused by 2013 Lushan earthquake and its relationship with 2008 Wenchuan earthquake.ScienceChina:EarthSciences, 56(7):1169-1176, doi: 10.1007/s11430-013-4642-1. Shi Y T, Gao Y, Zhao C P, et al. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence.ChineseJ.Geophys. (in Chinese), 52(2): 398-407. Sun C Q, Chen Y, Gao E G. 2011. The crustal anisotropy and its geodynamical significance of the strong basin-range interaction zone beneath the east margin of Qinghai-Tibet plateau.ChineseJ.Geophys. (in Chinese), 54(5): 1205-1214,doi:10.3969/j.issn.0001-5733.2011.05.009.

Teng J W, Li S L, Zhang Y Q, et al. 2014. Fine velocity structures and deep processes in crust and mantle of the Qinling orogenic belt and the adjacent North China craton and Yangtze craton.ChineseJ.Geophys. (in Chinese), 57(10): 3154-3175,doi:10.6038/cjg2014006.

Unsworth M J, Jones A G, Wei W, et al. 2005. Crustal rheology of the Himalaya and Southern Tibet inferred from magnetotelluric data.Nature, 438(7064): 78-81.

Wang C Y, Wu J P, Lou H, et al. 2003a. P-wave crustal structure in western Sichuan and eastern Tibetan region.ScienceChina(SeriesD), 46(Supp.): 254-265.

Wang C Y, Han W B, Wu J P, et al. 2003b. Crustal structure beneath the Songpan-Garze orogenic belt.ActaSeismologicaSinica(in Chinese),25(3) :229-241.

Wang C Y, Han W B, Wu J P, et al. 2007. Crustal structure beneath the eastern margin of the Tibetan Plateau and its tectonic implications.J.Geophys.Res., 112(B7): B07307, doi: 10.1029 /2005 JB 003873.

Wang C Y, Lou H, Yao Z X, et al. 2010. Crustal thicknesses and poisson′s ratios in Longmenshan mountains and adjacent regions.QuaternarySciences(in Chinese), 30(4): 652-661.

Wang F Y, Zhang X K. 2004. Genetic algorithms in seismic waveform inversion and its application in seismic wide-angle reflection and refraction data.∥Chen Yuntai, Teng Jiwen,Kan Rongjv,et al.,eds.Advances in Seismology and Physics of the Earth′s Interior in China: in Honor of Academician Zeng Rongsheng′s 80th Birthday (in Chinese) Beijing: Seismological Press, 256-269.

Wang Q, Gao Y. 2014. Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the Tibetan Plateau.ScienceChina:EarthSciences, 57(16): 2532-2542, doi: 10.1007/s11430-014-4908-2.

Wang S J, Zhang X K, Zhang C K, et al. 2007. 2-D crustal structures along Wuqing-Beijing-Chicheng deep seismic sounding profile.ChineseJ.Geophys. (in Chinese), 50(6): 1769-1777.

Wang S J, Wang F Y, Zhang J S, et al. 2014. The P-wave velocity structure of the lithosphere of the North China Craton—Results from the Wendeng-Alxa Left Banner deep seismic sounding profile.ScienceChina:EarthSciences, 57(9):2053-2063.

Wang W M, Hao J L, Yao Z X. 2013. Preliminary result for rupture process of Apr. 20, 2013,MW6.7 earthquake, Lushan, Yaan, China.ChineseJ.Geophys. (in Chinese), 56(4): 1412-1417,doi:10.6038/cjg20130436. Wang Y S, Quan Q, Luo Y H, et al. 2013. Geological environmental influence of LushanMS7.0 earthquake in Sichuan.JournalofEarthSciencesandEnvironment(in Chinese), 35(2): 92-98.

Wang Z, Fukao Y, Pei S P. 2009. Structural control of rupturing of theMw7.9 2008 Wenchuan Earthquake, China.EarthPlanet.Sci.Lett., 279(1-2): 131-138. Wei W B, Unsworth M, Jones A, et al. 2001. Detection of widespread fluids in the Tibetan crust by magnetotelluric studies.Science, 292(5517): 716-718. Wu J P, Ming Y H, Wang C Y. 2001. The S wave velocity structure beneath digital seismic stations of Yunnan province inferred from teleseismic receiver function modeling.ChineseJ.Geophys. (in Chinese), 44(2): 228-237.

Xie Z J, Jin B K, Zheng Y, et al. 2013. Source parameters inversion of the 2013 Lushan earthquake by combining teleseismic waveforms and local seismograms.ScienceChina:EarthSciences, 56(7):1177-1186, doi: 10.1007/s11430-013-4640-3.

Xu T, Xu G M, Gao E G, et al. 2004. Block modeling and shooting ray tracing in complex 3-D media.ChineseJ.Geophys.(in Chinese), 47(6): 1118-1126.Xu T, Xu G M Gao E G, et al. 2006. Block modeling and segmentally iterative ray tracing in complex 3D media.Geophysics, 71(3): T41-T51. Xu T, Zhang Z, Gao E, et al. 2010a. Segmentally iterative ray tracing in complex 2D and 3D heterogeneous block models.Bull.Seismol.Soc.Am., 100(2): 841-850.

Xu T, Li F, Wu Z B, et al. 2014. A successive three-point perturbation method for fast ray tracing in complex 2D and 3D geological model.Tectonophysics, 627: 72-81.

Xu T, Zhang M H, Tian X B, et al. 2014. Upper crustal velocity of Lijiang-Qingzhen profile and its relationship with the seismogenic environment of theMS6.5 Ludian earthquake.ChineseJ.Geophys. (in Chinese), 57(9): 3069-3079,doi:10.6038/cjg20140932.Xu X W, Chen G H, Yu G H, et al. 2013. Seismogenic structure of Lushan earthquake and its relationship with Wenchuan earthquake.EarthScienceFrontiers(in Chinese), 20(3): 11-20.

Xu Y, Li Z W, Huang R, et al. 2010b. Seismic structure of the Longmen Shan region from S-wave tomography and its relationship with the WenchuanMS8.0 earthquake on 12 May 2008, southwestern China.Geophys.Res.Lett., 37: L02304.

Yan Z K, Li Y, Zhao G H, et al. 2014. The relationship between Lushan earthquake and Wenchuan earthquake by segmentation of geology and geomorphology of Longmen Shan.ChineseJournalofNature(in Chinese), 36(1): 51-57.

Yao H J, Campman X, de Hoop M V, et al. 2009. Estimation of surface wave Green′s functions from correlation of direct waves, coda waves, and ambient noise in SE Tibet.Phys.EarthPlanet.Int., 177(1-2): 1-11. Zhang P Z, Wen X Z, Shen Z K, et al. 2010. Oblique, high-angle, listric-reverse faulting and associated development of strain: The Wenchuan earthquake of May 12, 2008, Sichuan, China.Annu.Rev.EarthPlanet.Sci., 38(1): 353-382.

Zhang Y, Xu L S, Chen Y T. 2013. Rupture process of the Lushan 4. 20 earthquake and preliminary analysis on the disaster-causing mechanism.ChineseJ.Geophys.(in Chinese), 56(4): 1408-1411,doi:10.6038/cjg20130435.

Zhang Z J, Teng J W, Badal J, et al. 2009a. Construction of regional and local seismic anisotropic structures from wide-angle seismic data: crustal deformation in the southeast of China.J.Seismol., 13(2):241-252, doi: 10.1007/s10950-008-9124-0.

Zhang Z J, Wang Y H, Chen Y, et al. 2009b. Crustal structure across Longmenshan fault belt from passive source seismic profiling.Geophys.Res.Lett., 36(17).

Zhao B, Gao Y, Huang Z B, et al. 2013. Double difference relocation, focal mechanism and stress inversion of LushanMS7.0 earthquake sequence.ChineseJ.Geophys. (in Chinese), 56(10): 3385-3395,doi:10.6038/cjg20131014.

Zheng Y, Yang Y J, Ritzwoller M H, et al. 2010. Crustal structure of the northeastern Tibetan Plateau, the Ordos block and the Sichuan basin from ambient noise tomography.Earthq.Sci., 23(5): 465-476.

Zheng Y, Ge C, Xie Z J, et al. 2013. Crustal and upper mantle structure and the deep seismogenic environment in the source regions of the Lushan earthquake and the Wenchuan earthquake.ScienceChina:EarthSciences, 56(7):1158-1168.

附中文参考文献

陈立春, 冉勇康, 王虎等. 2013. 芦山地震与龙门山断裂带南段活动性. 科学通报, 58(20): 1925-1932.

陈运泰, 杨智娴, 张勇等. 2013. 从汶川地震到芦山地震. 中国科学D辑: 地球科学, 43(6): 1064-1072.

丁志峰, 武岩, 王辉等. 2008. 2008年汶川地震震源区横波分裂的变化特征. 中国科学D辑:地球科学, 38(12): 1600-1604.

杜方, 龙锋, 阮祥等. 2013. 四川芦山7.0级地震及其与汶川8.0级地震的关系. 地球物理学报, 56(5): 1772-1783,doi:10.6038/cjg20130535.

段永红, 张先康,方盛明.2002.华北地区上部地壳结构的三维有限差分层析成像. 地球物理学报, 45(3): 362-369.

房立华, 吴建平, 王未来等. 2013. 四川芦山MS7.0级地震及其余震序列重定位. 科学通报, 58(20): 1901-1909.

高锐, 董树文, 贺日政等. 2004. 莫霍面地震反射图像揭露出扬子陆块深俯冲过程. 地学前缘, 11(3): 43-49.

高原, 王琼, 赵博等. 2013. 龙门山断裂带中南段的一个破裂空段—芦山地震的震后效应. 中国科学D辑: 地球科学, 43(6): 1038-1046.

国家地震局科技监测司.1988.中国大陆深部构造的研究与进展. 北京: 地质出版, 18-46.

雷建设, 赵大鹏, 苏金蓉等. 2009. 龙门山断裂带地壳精细结构与汶川地震发震机理. 地球物理学报, 52(2): 339-345.

李德威, 陈桂凡, 陈继乐等. 2013. 地震预测从芦山地震到大陆地震. 地学前缘, 20(3): 1-10.

刘保金, 张先康, 酆少英等. 2009. 龙门山山前彭州隐伏断裂高分辨率地震反射剖面. 地球物理学报, 52(2): 538-546.

刘成利, 郑勇, 葛粲等. 2013. 2013年芦山7.0级地震的动态破裂过程. 中国科学 地球科学, 43(6): 1020-1026.

楼海, 王椿镛, 吕智勇等. 2008. 2008年汶川MS8.0级地震的深部构造环境—远震P波接收函数和布格重力异常的联合解释. 中国科学D辑: 地球科学, 38(10): 1207-1220.

吕坚, 王晓山, 苏金蓉等. 2013. 芦山7.0级地震序列的震源位置与震源机制解特征. 地球物理学报, 56(5): 1753-1763,doi:10.6038/cjg20130533.

单斌, 熊熊, 郑勇等. 2013. 2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系. 中国科学D辑: 地球科学, 43(6): 1002-1009.

石玉涛, 高原, 赵翠萍等. 2009. 汶川地震余震序列的地震各向异性. 地球物理学报, 52(2): 398-407.

孙长青, 陈赘, 高尔根. 2011. 青藏高原东缘强烈盆山相互作用区的地壳各向异性特征及其动力学意义探讨. 地球物理学报, 54(5): 1205-1214，doi:10.3969/j.issn.0001-5733.2011.05.009.

滕吉文, 李松岭, 张永谦等. 2014. 秦岭造山带与邻域华北克拉通和扬子克拉通的壳、幔精细速度结构与深层过程. 地球物理学报, 57(10): 3154-3175, doi: 10.6038/cjg20141006.

王椿镛, 吴建平, 楼海等. 2003a. 川西藏东地区的地壳P波速度结构. 中国科学D辑: 地球科学, 33(B04): 181-189.

王椿镛, 韩渭宾, 吴建平等. 2003b. 松潘—甘孜造山带地壳速度结构. 地震学报, 25(3): 229-241.

王椿镛, 楼海, 姚志祥等. 2010. 龙门山及其邻区的地壳厚度和泊松比. 第四纪研究, 30(4): 652-661.

王夫运, 张先康. 2004. 波形反演的遗传算法—地震宽角反射/折射数据.∥陈运泰，腾吉文，阚荣举等.中国大陆地震学与地球内部物理学进展—庆贺曾融生院士八十寿诞. 北京: 地震出版社,256-269.

王琼, 高原. 2014. 青藏东南缘背景噪声的瑞利波相速度层析成像及强震活动. 中国科学D辑: 地球科学, 44(11):2440-2450.

王帅军, 张先康, 张成科等. 2007. 武清—北京—赤城二维地壳结构和构造. 地球物理学报, 50(6): 1769-1777.

王卫民, 郝金来, 姚振兴. 2013. 2013年4月20日四川芦山地震震源破裂过程反演初步结果. 地球物理学报, 56(4):1412-1417,doi:10.6038/cjg20130436.

王运生, 全清, 罗永红等. 2013. 四川芦山MS7.0级地震的地质环境影响分析. 地球科学与环境学报, 35(2): 92-98.

吴建平, 明跃红,王椿镛. 2001. 云南数字地震台站下方的S波速度结构研究. 地球物理学报, 44(2): 228-237.

谢祖军, 金笔凯, 郑勇等. 2013. 近远震波形反演2013年芦山地震震源参数. 中国科学D辑: 地球科学, 43(6): 1010-1019.

徐涛, 徐果明, 高尔根等. 2004. 三维复杂介质的块状建模和试射射线追踪. 地球物理学报, 47(6): 1118-1126.

徐涛, 张明辉, 田小波等. 2014. 丽江—清镇剖面上地壳速度结构及其与鲁甸MS6.5级地震孕震环境的关系. 地球物理学报, 57

(9): 3069-3079,doi:10.6038/cjg20140932.

徐锡伟, 陈桂华, 于贵华等. 2013. 芦山地震发震构造及其与汶川地震关系讨论. 地学前缘, 20(3): 11-20.

颜照坤, 李勇, 赵国华等. 2014. 从龙门山地质地貌分段性探讨芦山地震与汶川地震的关系. 自然杂志, 36(1): 51-57.

张勇, 许力生, 陈运泰. 2013. 芦山4.20地震破裂过程及其致灾特征初步分析. 地球物理学报, 56(4): 1408-1411,doi:10.6038/cjg20130435.

赵博, 高原, 黄志斌等. 2013. 四川芦山MS7.0地震余震序列双差定位、震源机制及应力场反演. 地球物理学报, 56(10): 3385-3395,doi:10.6038/cjg20131014.

郑勇, 葛粲, 谢祖军等. 2013. 芦山与汶川地震震区地壳上地幔结构及深部孕震环境. 中国科学D辑: 地球科学, 43(6): 1027-1037.

(本文编辑 胡素芳)

The deep seismogenic environment of LushanMS7.0 earthquake zone revealed by a wide-angle reflection/refraction seismic profile

WANG Shuai-Jun， WANG Fu-Yun*， ZHANG Jian-Shi， LIU Bao-Feng， ZHANG Cheng-Ke, ZHAO Jin-Ren, DUAN Yu-Ling, SONG Xiang-Hui， DENG Xiao-Guo，MA Ce-Jun， SUN Yi-Nan， ZANG Yi-Ran， LI Yi-Qing

GeophysicalExplorationCenterofChinaEarthquakeAdministration,Zhengzhou450002,China

To study the deep tectonic environment and the seismogenic background of the Lushan earthquake in November, 2013, we completed an about 410 km-long high resolution wide angle reflection / refraction profile, which went across the Yangtze block and the Songpan-Garzê block and stretched over a number of major faults. By the well-desighed and thorough spread geometry, we obtained the artificial seismic recording data with higher SNR along the profile. Through 7 groups of clear crustal P wave phase from 8 blasting sections in this paper, we applied the seismic ray travel time forward fitting to construct 2D P wave velocity model of the crust and the upper mantle. Combined with the reliable tracking interval and amplitude variations in different interface phases of the upper crust, we analyzed and discussed the characteristics of crust-mantle velocity distribution, crustal interface distribution and the deep structure of fractures along the profile.

王帅军，王夫运，张建狮等. 2015. 利用宽角反射/折射地震剖面揭示芦山MS7.0地震震区深部孕震环境.地球物理学报，58(9):3193-3204,

10.6038/cjg20150915.

Wang S J, Wang F Y, Zhang J S, et al. 2015. The deep seismogenic environment of LushanMS7.0 earthquake zone revealed by a wide-angle reflection/refraction seismic profile.ChineseJ.Geophys. (in Chinese),58(9):3193-3204,doi:10.6038/cjg20150915.

10.6038/cjg20150915

P315

2015-05-20，2015-07-08收修定稿

中国地震局四川省芦山“4.20”7.0级强烈地震科学考察项目和国家自然科学基金项目(41474076)共同资助.

王帅军,男,1978年生,副研究员,主要从事地震测深资料处理与解释等方面的研究. E-mail：wsjdzj@126.com

*通讯作者 王夫运，男，1962年生，研究员，主要从事地震波走时和波形反演理论研究. E-mail：fuyunwang@x263.net