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.
-
Keywords:
- Longmenshan fault zone /
- deep seismic sounding /
- SIRT /
- velocity structure /
- travel time inversion
-
引言
2022年1月8日1时45分，青海省海北州门源县发生MS6.9地震，震中位置为（37.8°N，101.3°E），这是中国大陆2022年第一个接近M7.0的强震。该次地震发生在青藏高原东北缘祁连地块中的托莱山—冷龙岭断裂附近（图1）。据中国地震台网中心（2022）的速报结果，该地震震中距离门源县城大约54 km，震源深度为10 km，震源深度重定位结果为12.9 km （Fan et al，2022）。截至2022年1月19日8时整，共记录到M≥3.0余震20次，其中MS4.0—4.9地震5次，MS≥5.0地震2次，最大为1月12日18时20分MS5.2地震。
![]()
图 1 青藏高原东北缘区域构造背景与地震台站分布Figure 1. Tectonic settings of the northeastern Tibetan Plateau and distribution of seismic stations门源地区位于祁连地块中东部，由于受到印度板块和欧亚板块持续碰撞、挤压的远程作用，区域内活动断裂发育、构造变形强烈。祁连地块内部发育了一系列以挤压逆冲为主兼具左旋走滑特征的活动断裂带，其中主要受到阿尔金和祁连—海原两条大型边界走滑断裂控制（袁道阳等，2004）。2022年门源MS6.9地震发生在冷龙岭断裂西端、托莱山断裂东端，基本上将冷龙岭断裂与托莱山断裂连接起来（图2），这两条断裂与海原断裂、老虎山断裂、毛毛山断裂、金强河断裂等构成了长约1 000 km的祁连—海原断裂带。祁连—海原断裂带不仅控制着青藏高原东北缘地区的几何和构造格局（Zheng et al，2013；Shi et al，2020），同时也调节着青藏地块相对于阿拉善地块的向东运动（Gaudemer et al，1995；Lasserre et al，2002）。
![]()
图 2 门源MS6.9地震及余震分布2022年地震序列定位结果引自Fan等（2022），2016年门源地震位置引自梁姗姗等（2017），1986年门源地震位置引自兰州地震研究所青海省地震局联合考察队（1987）。F1：龙首山南缘断裂东段；F2：河西堡—四道山断裂；F3：榆木山东缘断裂；F4：民乐—永昌断裂；F5：肃南—祁连断裂；F6：民乐—大马营断裂；F7：皇城—双塔断裂；F8：托莱山断裂；F9：冷龙岭断裂；F10：金强河断裂；F11：毛毛山断裂；F12：天桥沟—黄羊川断裂；F13：玛雅雪山北缘断裂；F14：大通山北缘断裂；F15：木里—江仓断裂；F16：日月山断裂；F17：门源断裂；F18：达坂山断裂；F19：黑林河断裂Figure 2. The distribution of Menyuan MS6.9 earthquake and its aftershocksThe location of 2022 seismic sequence is from Fan et al （2022），the location of 2016 Menyuan earthquake is from Liang et al （2017），and the location of 1986 Menyuan earthquake is from Lanzhou Institute of Seismology and Seismological Bureau of Qinghai Province （1987）. F1：Eastern segment of southern Longshoushan fault；F2：Hexipu-Sidaoshan fault；F3：Yumushan eastern marginal fault；F4：Minyue-Yongchang fault； F5：Su’nan-Qilian fault；F6：Minyue-Damaying fault；F7：Huangcheng-Shuangta fault；F8：Tuolaishan fault；F9：Lenglongling fault；F10：Jinqianghe fault；F11：Maomaoshan fault；F12：Tianqiaogou-Huangyangchuan fault；F13：Mayaxueshan northern marginal fault；F14：Datongshan northern marginal fault；F15：Muli-Jiangcang fault；F16：Riyueshan fault； F 17：Menyuan fault；F18：Dabanshan fault；F19：Heilinhe fault1900年以来，祁连—海原断裂带附近多次发生M6.0以上地震，位置相近的有1986年MS6.4地震和2016年MS6.4地震（图2），其中2016年门源MS6.4地震发生在冷龙岭北侧断裂（胡朝忠等，2016），与此次门源地震序列中的1月12日MS5.2余震位置相差约9 km。两次门源地震虽然相距不远，但其发震机制表现不同，即2016年门源地震为冷龙岭北侧断裂上的逆冲型地震，2022年门源地震则为道沟断裂上的走滑型地震（中国地震局地质研究所，2022），表明该地区的地质构造特征比较复杂。相关研究表明，在剪切变形较大的区域更易发生地震（Jin et al，2019）。位于青藏东北缘的祁连地块是新生代再次活跃的早古生代岛弧、微陆块拼合体（Yin，Harrison，2000），当前仍以12 mm/a的速率沿NNE−SSW方向水平缩短（Zhang et al，2004）。这种NE方向的推挤对阿尔金断裂带东部和河西走廊产生一定的加载作用，可能是导致2016年门源MS6.4地震发生的原因（胡朝忠等，2016）。2022年门源MS6.9地震发生后，中国地震局组织了考察工作，初步结果显示：此次门源MS6.9地震极震区的烈度可能达Ⅸ以上（中国地震局地球物理研究所，2022）；地震序列在托莱山—冷龙岭断裂附近产生了四条地表破裂带，其中断裂带南、北两侧破裂程度差异较大，北侧明显高于南侧（中国地震局地质研究所，2022）。
结构成像研究显示，强震更容易发生在高、低速变化的过渡带，例如青藏东南缘（王琼，高原，2014），青藏东北缘地区的岷漳MS6.7地震和九寨沟MS7.0地震的震源也都位于高、低速的过渡带（夏思茹等，2021）。为解剖强震发生的深部构造背景，本文收集了速度结构等资料，拟通过分析地壳结构与地震分布之间的关系，探讨门源MS6.9地震的深部孕震构造背景。
1. 接收函数揭示的门源地区地壳厚度和vP/vS不均匀分布
地壳厚度和泊松比是描述地壳结构和物质成分的两个重要参数，可以为地壳的地质演化过程提供重要的约束，同时有益于对地震孕震环境的研究。Wang等（2016）使用来自中国地震台网2007年8月至2013年10月间的三分量宽频带地震数据，对震级M＞5.0，震中距为30°—90°且震相清晰、初动尖锐的地震记录进行接收函数计算，得到了青藏高原东北缘地区的地壳厚度和波速比。Wang等（2017）使用来自中国科学台阵Ⅱ期和鄂尔多斯地块内流动台阵的资料（共724个台站）得到了青藏高原东北缘的地壳厚度和波速比。为了分析此次门源MS6.9地震的发震位置与深部结构之间的关系，本文根据两种数据结果重绘了门源地区的地壳厚度和vP/vS的分布图像（图3）。
![]()
从图3可以看到，研究区域内地壳厚度约为30—70 km，变化非常剧烈，青藏高原东北部的地壳比周缘区域厚。地壳厚度与地形起伏变化比较一致，从东向西呈梯度加深的趋势。门源台（MEY）的地壳厚度为59.3 km，与王椿镛等（1995）的基于门源—平凉—渭南人工地震测深剖面得到的门源地壳厚度（58 km）较为一致。吴立辛等（2011）运用小波多尺度分析也观察到青藏高原东北缘地壳厚度呈现由东北向西南增厚的趋势，且门源地区位于梯度变化较大的区域。通过本文的数据再分析，门源地区处于地壳厚度由西向东剧烈减薄的地带，而门源MS6.9地震的震中位于地壳厚度南厚北薄的局部急剧变化过渡带上，同时又处于波速比（vP/vS）东高西低快速变化的过渡带，揭示了门源MS6.9地震与深部结构的不均匀分布有关，特别是与物性参数的快速变化有关。
由于使用不同的资料和方法，Wang等（2016）和Wang等（2017）的结果在门源地区存在一些差异，也揭示了该区域地壳物性在小范围内的不均匀分布，表明关键地域密集台阵观测的必要性。这些结果虽然呈现不均匀分布的具体形态差异，但它们都揭示了门源MS6.9地震震源区深部结构存在剧烈的（或快速的）物性变化。
2. 体波成像揭示的P波、S波速度和泊松比分布
基于2009年至2017年青藏高原东北缘区域固定台网记录到的3万9 971次地震的初至P波和S波走时数据，利用双差成像方法对2万9 491次地震进行了重新定位，获取了青藏高原东北缘地壳的波速和泊松比结构（肖卓，高原，2017）。为了探究门源MS6.9地震的深部发震结构，对门源地震震源区的深部结构进行了重绘放大，得到了门源地震震源区的P波速度、S波速度以及泊松比在不同深度的横向和垂向分布图像，结果如图4和图5所示。
![]()
图 4 门源地区P波速度（a）、S波速度（b）和泊松比（c）红色星形表示2022年MS6.9门源主震，红色圆点表示其余震，地壳深度标于子图的左下角，下同Figure 4. P-wave velocity （a），S-wave velocity （b） and Poisson’s ratio （c） in the Menyuan regionThe red star represents the 2022 Menyuan MS6.9 earthquake，the red solid circle represents its aftershock，and the crustal depth is labled at the lower-left corner of the subfigure，the same below![]()
图 5 穿过门源地震震源区的P波速度（a）、S波速度（b）和泊松比（c）的垂向剖面图剖面方向根据图4白色实线绘制，剖面数据来自肖卓和高原（2017）Figure 5. Vertical profiles of P-wave velocity （a），S-wave velocity （b） and Poisson’s ratio （c） through the source region of Menyuan earthquakeThe section direction is drawn according to the white line in Fig. 4，and the profile data is from Xiao and Gao （2017）P波速度分布显示，门源地震震源区的结构存在明显的垂向变化和横向变化，特别是在垂向上，上地壳顶部（0—10 km）的P波速度较高，但随着深度的增加（10—20 km），P波速度逐渐转变为低速异常（图4，5）。根据门源地震序列重新定位结果（Fan et al，2022），MS6.9主震的震源深度为12.9 km。在震源区，P波速度在5—15 km深度由高速变为低速（图5）。S波速度和泊松比分布显示，门源地震震源区的结构存在明显的横向变化。在10—15 km深度，门源地震震源区东西两侧分别呈现为高S波速度、低泊松比和低S波速度、高泊松比异常（图4）。门源地震震源区泊松比（图4）在东西方向的变化特征与接收函数得到的vP/vS变化（图3）具有很好的一致性，揭示了地壳介质物性在东西方向上的快速变化很可能主要分布在10—20 km深度范围。
夏思茹等（2021）使用青藏高原东北缘71个固定台和418个流动台站得到了更高分辨率的三维P波速度结构，本文使用其数据绘制了两个跨过震中位置的垂向剖面图（图6）。从图6可以看到，门源地震发生在上地壳高低速变化过渡区，震源下方存在明显低速体。从这一点上看，图6与图5的P波速度垂向变化结果基本一致。
![]()
图 6 穿过门源地震震源区P波速度垂向剖面图剖面方向根据图1中白色实线绘制，剖面数据来自夏思茹等（2021）Figure 6. Vertical P-wave velocity profile through the source region of Menyuan earthquakeThe section direction is drawn according to the white line in Fig. 1，and the prolife data is from Xia et al （2021）人工地震测深剖面（王椿镛等，1995）和远震面波成像（Li et al，2017）均揭示，祁连造山带下方中地壳存在S波低速层；地震与速度结构之间关系的研究结果表明，地震的发生往往与壳内低速度层存在明显的关联性（Zhao et al，2002；李永华等，2014）；重力资料多阶小波分析也显示，重力异常分布与大地震发生的位置有关联性（Wu，Gao，2019）；王新胜等（2013）利用分离的布格重力异常揭示大地震多发生在高、低密度异常边界区域或者下地壳低密度层之上，结果显示这次门源MS6.9震源区下方下地壳表现为强烈的低密度异常；赵凌强等（2019）通过大地电磁探测剖面发现冷龙岭断裂下方存在高导层，认为该断裂下方可能形成了明显的力学强度软弱区，这种力学强度软弱区的存在易导致地震蠕动、滑移与发生。由此可见，地震带下方壳内薄弱层的存在很可能使其上覆脆性上地壳物质易于形成应力集中而发生强震。
3. 背景噪声成像揭示的相速度和方位各向异性分布
利用青藏高原东北缘地区固定地震台网三分量连续波形数据，采用Yao等（2006）发展起来的背景噪声数据处理技术，使用瑞雷波能量较强的垂向记录，反演得到了青藏高原东北缘地区相速度和方位各向异性（王琼，高原，2018）。为了分析门源地震与深部结构的关系，本文重绘了门源地震震源区的区域相速度和方位各向异性结构（图7）。
![]()
图 7 门源地震震源区相速度和方位各向异性分布Figure 7. Phase velocity and azimuthal anisotropy distribution in the source region of the Menyuan earthquake（a） T＝8 s；（b） T＝12 s根据门源MS6.9地震的震源深度结果，把门源地震标示在周期8 s和12 s的相速度图（图7）上。可以看到，MS6.9门源地震发生在高低速的过渡区，方位各向异性也呈现出方向和幅值的变化。郭瑛霞（2017）利用布设在祁连山断裂带的40个短周期流动台站和38个宽频带固定台站的记录，采用背景噪声成像得到了祁连地块S波速度结构，结果同样显示门源震源区处于高速与低速的过渡区域。王琼和高原（2014）分析青藏高原东南缘地区强震活动与速度结构关系时发现，M5.0—6.9地震主要发生在高低速过渡区和低速异常区，M≥7.0地震则主要发生在高低速过渡区但更深入到高速异常区里。房立华等（2009）认为，速度变化强烈的地区存在应力集中，介质比较脆弱，更易于释放应力而触发较大的地震。
地震各向异性可以揭示壳幔内部构造变形等信息（高原，滕吉文，2005）。连续记录的背景噪声数据不依赖于地震信号，更有利于获得地壳的结构特征。从图7可以看到，以冷龙岭断裂为界，各向异性快波方向呈现不同特征：断裂以北，快波方向主要呈北东方向，与区域平均最大主压应力方向一致（许忠淮，2001）；断裂以南，快波方向主要呈北西方向，表明各向异性形态可能受到区域断裂带的影响，揭示门源地区地壳（至少上地壳）介质各向异性可能受到区域构造应力和局部断裂（构造）的双重约束。
祝意青等（2016）关于2011—2015年流动重力观测资料的分析结果显示，重力变化在冷龙岭断裂两侧出现正负边界带，重力变化梯度带与冷龙岭断裂的走向一致，冷龙岭北侧的重力变化以正值为主，南侧则以负值为主，表明冷龙岭断裂有地震孕育的构造背景。赵凌强等（2019）结合现今水准、重力、GPS速度场及大地电磁测深结果，认为冷龙岭断裂所处的祁连—海原断裂带是青藏高原东北缘的边界断裂，目前正承受着巨大的NE向挤压力。Gao等（1996）和高锐等（1998）根据深地震反射剖面结果认为青藏高原在向北扩展过程中处于双向挤压应力作用之下，同时受到印度板块向北俯冲和阿拉善地块沿宽滩山断裂插入到祁连山之下的动力作用。在双向挤压力作用下，上地壳岩片发生逆冲叠覆，岩层发生脆性变形，使上地壳缩短增厚，因此对于祁连造山带内低速高导层的存在，可能与地壳增厚和地表抬升有关。Wang （2001）提出，由于壳幔流变性质的差异，地壳的增厚升温可以降低岩石层的整体强度。新生代祁连造山带地壳的增厚和构造演化可能导致放射性生热元素含量增加，产生相对较高的热流值。羌塘地块和松潘—甘孜地块的热流值较低（44—45 mW/m2），柴达木盆地的热流值为54 mW/m2，而祁连造山带的热流值高达66 mW/m2 （Wang，2001），说明祁连山的局部低速异常可能由上地壳增厚、生热元素含量增加引起，而门源地震同时也处于地壳厚度快速变化的区域（图3），容易造成应力的集中，具有发生大地震的深部构造背景。
4. 讨论与结论
本文根据地壳厚度、vP/vS比值、P波和S波速度结构、面波相速度及方位各向异性结果，结合密集地震台阵探测、人工地震测深、深地震反射、大地电磁探测、区域重力场等地球物理资料，初步探讨了门源MS6.9地震的深部构造背景。
此次门源MS6.9地震发生在地壳厚度和vP/vS值都出现快速空间变化的区域。门源MS6.9地震震源位于P波速度从浅到深由高速变低速的垂向过渡区，也是S波速度和泊松比分布呈现明显横向变化的过渡区域，该过渡区域大致处于10—20 km深度范围内。背景噪声成像结果显示，冷龙岭断裂两侧相速度和方位各向异性变化明显。这些现象表明，地震活动与地壳结构有较强的对应关系。通常情况下，地震波速度的变化可能代表岩石强度的变化，而大的破裂或强烈的地壳变形往往都集中在地壳强度（或地震波速度）反差较大的地方（Tian et al，2021）。
1月12日的MS5.2余震震中与2016年MS6.4地震震中很近，表明这次门源MS6.9地震及其余震导致冷龙岭断裂的破裂比较充分，两次门源地震之间及邻近地区短时间内难以积累更大的能量，即意味着短时间内发生更大地震的可能性不大。
本文结果显示，不同资料、不同方法及不同数据处理得到的结果可能存在差异，需要甄别大的形态一致性与小的局部不同。重要的是强震震源区处于地壳介质物性快速变化的过渡区，这是不同的研究结果给出了共同的认识。对于小区域深部结构的准确成像，需要密集台阵观测资料的支持。此外，P波与S波速度分布不同，可能表明岩性、组分或流体充填状态等的不同。
受到印度板块和欧亚板块持续碰撞、挤压的远程作用，青藏高原东北缘祁连地块的地壳显著增厚，但陇西盆地没有明显增厚（Tian et al，2021），而祁连地块东部正处于地壳厚度横向变化较为剧烈的区域，容易形成应力集中和积累，具有发生大地震的深部构造背景。
中国地震局地球物理研究所王兴臣副研究员提供了地壳厚度和波速比数据，石磊副研究员和硕士研究生夏思茹提供了P波速度剖面图件，审稿专家提出的宝贵意见使本文的质量有了明显改善，作者在此一并表示衷心的感谢。
-
![]()
图 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
-
期刊类型引用(3)
1. 赵庆旭，王延伟，莫红艳，曹振中. 基于多输入高斯过程回归的震级快速估算方法. 地震学报. 2024(05): 806-824 . 
本站查看
2. 李子昊，赵国为，郭攀，赵亚楠，王晓睿. 人工智能在地震学中的应用. 防灾减灾学报. 2024(04): 84-90 . 百度学术
3. 徐繁树，王保云，韩俊. 一种沟谷型潜在泥石流危险性评价方法：基于多源数据融合的卷积神经网络. 地球信息科学学报. 2023(03): 588-605 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 717
- HTML全文浏览量: 440
- PDF下载量: 30
- 被引次数: 5
