Crustal thickness and vP/vS ratio of Arxan volcanic group in Xing’an-Mongolia orogenic belt
-
摘要: 利用架设在我国东北地区阿尔山火山区的宽频带流动地震台站记录的远震波形数据，采用P波接收函数H-κ叠加扫描方法，得到了阿尔山火山区的地壳厚度和平均波速比。结果显示：阿尔山地区的地壳厚度范围为33.9—37.9 km，整体呈西北厚东南薄的特点，火山带附近地壳较薄；地壳平均波速比范围为1.73—1.83，主要有柴河镇—明水河镇、伊敏德仁北部和天池镇三处高波速比区。结合前人研究结果推断， 阿尔山火山区薄的地壳和高波速比值可能是由地幔物质上涌、玄武岩浆底侵下地壳所致。Abstract: The Arxan volcanic group is one of the most important volcanic groups in the Quaternary in Northeast China. It is of great significance to obtain information of crustal structure and vP/vS ratios in this area for studying the physical properties of the crust and tectonic evolution beneath the volcanoes. Original waveform data came from temporary broadband seismic stations in the Arxan volcanic area. The crustal thickness and average vP/vS in the study area are obtained by the P-wave receiver function and H-κ stacking method. The results show that the crust thickness is thicker in the northwest and thinner in the southeast, ranging from 33.9 km to 37.9 km, with thinner crust near the volcanic belt. There is a positive correlation between crustal thickness and elevation. The vP/vS ratio ranges from 1.73 to 1.83. There are three high vP/vS ratio areas in the study area: Chaihe town-Mingshuihe town, the north of Yimin Doren and Tianchi town. The thin crust and high wave velocity ratio areas beneath Arxan volcanic region may be caused by the upwelling of mantle material and underplating of basaltic magma.
-
Keywords:
- Arxan volcanic region /
- receiver function /
- H-κ stacking /
- crustal thickness /
- vP/vS ratio
-
引言
紫坪铺水库位于四川省都江堰市岷江上游的麻溪乡，于2005年9月开始蓄水，其正常蓄水位为877 m，当水库水位达到850 m的汛限水位时，水库水深最大约为100 m。大地构造位置分布显示，紫坪铺水库库区及邻近区域（103.3°E—103.7°E，30.85°N—31.15°N）位于青藏高原东缘的松潘—甘孜地块与扬子陆块过渡带上的龙门山断裂带内中南段，处于中央北川—映秀断裂与前山安县—灌县断裂之间（周宏伟等，2007）。
紫坪铺水库距离2008年汶川MS8.0地震震中仅十几千米，地震发生后，国内外许多学者对紫坪铺水库与汶川地震之间可能的关系展开了研究（Kerr，Stone，2009，2010；Deng et al，2010；Ruan et al，2017）。先前的研究大多使用四川区域地震台网（程万正等，2010）或当地紫坪铺水库地震台网（雷兴林等，2008；马文涛等，2011）数据资料，但是一些震级较小的地震由于信噪比较低往往会遗失在原始地震目录中，使地震目录的完整性降低。地震目录是区域地震活动性研究重要的基础资料，编目的完整性直接影响后续研究结果的可信性（Woessner，Wiemer，2005）。
地震的识别一般利用两类方法（Grigoli et al，2018；Zhang et al，2019）：基于震相拾取（Allen，1982；Withers et al，1998）方法和基于波形（Kao，Shan，2007；Shelly et al，2007；Zhang，Wen，2015）识别方法。在检测小震事件时，基于波形的方法通常比基于震相拾取的方法表现更好（Pesicek et al，2014；Grigoli et al，2018），其中基于波形互相关的匹配滤波技术（Match Filter）通过模板波形与连续波形中潜在小震信号之间的叠加互相关来检测小震事件，已得到了广泛的应用（Shelly et al，2007；Peng，Zhao，2009）。然而，该方法只能检测那些紧邻参考模板的地震，不能对检测到的地震事件进行定位。Zhang和Wen （2015）提出了匹配定位方法（Match and Locate），该方法采用了与匹配滤波技术中波形互相关检测类似的概念，但考虑了小震事件与模板可能存在的位置差异，同时提供了高精度的事件位置信息。但由于该方法需要对模板周围的三维空间进行搜索，因此运行起来非常耗时。为了提高匹配定位的计算效率和探测能力，Liu等（2020）发展了基于图像处理器加速的匹配定位（graphics processing unit-based match and locate，缩写为GPU-M&L）方法，综合考虑了几何扩散和地震波形的信噪比对模板事件各个分量赋予不同的权重系数，在通过图像处理器加速的同时也提高了地震事件的探测能力。
本文基于紫坪铺水库地震台网所提供的2005—2008年连续波形数据，首先以四川省地震局提供的地震目录重定位后的地震事件作为模板，使用GPU-M&L （Liu et al，2020）对紫坪铺水库及邻区进行遗漏地震事件的检测。然后，利用基于深度学习算法的去噪技术Deep Denoiser （Zhu et al，2019）对检测出的事件进行去噪处理，并设计了实现噪声和事件分类的卷积神经网络模型对检测出来的事件进行筛选。最后，再次使用相对定位方法hypoDD （Waldhauser，Ellsworth，2000）对筛选后的检测事件精定位，依据得到的地震目录和重定位结果进行地震时空分布特征、地震活动性与深部地质构造等分析，以期为后续的研究工作提供一定的数据支撑。
1. 数据和方法
本文使用四川省地震局提供的地震目录和紫坪铺水库地震台网2005—2008年的连续波形资料，研究区范围为（30.8°N—31.2°N，103.3°E—103.7°E），涉及7个地震台站，台站分布如图1所示。地震目录中记录研究区域内的地震事件1 241个，连续波形数据记录到了1 017次地震事件，其中最大震级为ML5.2，M＞4地震一共有17次。
![]()
图 1 研究区断裂、台站及模板事件重定位前后的震中分布图F1：汶川—茂汶断裂；F2北川—映秀断裂；F3安县—灌县断裂Figure 1. Distribution of the faults，stations，and epicenters of template events before and after relocation in the study areaF1：Wenchuan-Maowen fault；F2：Beichuan-Yingxiu fault；F3：Anxian-Guanxian fault数据处理主要分为以下几个步骤：
1） 利用深度学习算法PhaseNet对地震目录中事件进行P波、S波震相拾取，依据拾取结果挑选模板事件，并利用相对定位法hypoDD对其进行精定位；
2） 利用GPU-M&L方法进行微震检测；
3） 利用基于深度学习的去噪技术对检测到的事件进行去噪处理，并依据去噪后波形各分量的时频图，训练卷积神经网络模型对检测事件进行筛选；
4） 再次利用PhaseNet和hypoDD对筛选后的检测事件进行震相拾取和重定位。
1.1 模板事件选取与精定位
Zhu和Beroza （2019）发展了基于深度学习的震相拾取算法PhaseNet，该算法使用三分量地震波作为输入，输出P波、S波和噪声的概率分布。与现有方法相比，该算法在震相拾取的准确度方面有显著的改进，特别是对于S波。本文利用该算法对四川省地震局所提供地震目录中的事件进行P波、S波初动拾取，根据经验设置P波和S波的拾取可信度阈值为0.5，并舍去拾取台站数小于3的PhaseNet拾取结果。当在同一事件中拾取到多个相同震相时，保留可信度较高的震相到时，并利用TauP计算理论到时标记未被PhaseNet拾取到的震相，最终通过人工读取波形进行校正。基于上述步骤拾取的P波、S波初动，本文以S波到时前1 s为起始的6 s时间窗作为信号窗口，P波到时前6 s为起始的6 s时间窗作为噪声窗口，计算了地震事件每个分量的信噪比，并依据信噪比大于5的分量至少为9个的原则，共挑选了718个高质量的地震事件。
为提高模板事件震源位置的分辨率，本文基于PhaseNet得到的P波、S波到时记录，利用相对定位方法hypoDD对挑选的718个高质量事件进行了精定位。在该过程中，利用波形互相关技术对所有的地震对走时差进行校正，用于校正的波形窗口长度为P波到时前0.2 s至后1 s，滑动窗口为0.6 s；S波到时前0.5 s至后1.5 s，滑动窗口为1 s。地震对之间的最大震源位置差为20 km，最少震相连接个数为4，并选用前人层析成像结果作为初始速度模型。最终得到了578个模板事件的高分辨率震源位置（图1）。
1.2 基于图像处理器加速的匹配定位技术
基于上述重定位后的模板事件，本文利用基于图像处理器加速的匹配定位法（GPU-M&L）进行微震识别，以提高匹配定位的计算效率和探测能力，其程序与步骤主要遵循先前的研究（Liu et al，2020），仅在此简要描述。首先，对连续波形和模板事件都进行2—8 Hz带通滤波，来放大有用的地震信号，减轻背景噪声和大型事件尾波的影响，同时为减小运行内存提升运行速度，将数据从100 Hz减采样到20 Hz。然后，扫描时间窗口为模板事件S波到时前1 s至后5 s，并设置当在6 s的时间窗口内检测到多个事件时，只保留互相关系数（CC）值最高的事件。在扫描过程中，匹配定位法会将每一个网格点内各个台站分量的滑动互相关系数进行叠加，当平均滑动互相关系数大于所设阈值时，这一时刻被设置为一个可能事件的发震时刻，这一空间位置被设置为一个可能事件的震源位置，并根据所有分量中检测事件与模板事件之间峰值幅度比的中值计算其震级大小。
在以往的相关研究中，常常使用9—15倍的中值绝对偏差（MAD）作为检测阈值（Meng et al，2013；Walter et al，2015；Meng，Peng，2016；Ruan et al，2017）。为确定适用于本研究的阈值条件，表1统计了分别以9，12和15倍的MAD作为阈值条件下的检测事件数量和低CC （＜0.2）检测事件所占的比例。综合考虑检测数量与错误检测率，本文选取了12倍MAD值作为检测阈值。同时，由于低CC （＜0.2）值表明其可能是一个错误的检测（Meng et al，2013；Huang et al，2017），本文以CC＝0.2作为另一个条件对结果进行约束，最终得到2万零746个检测事件。图2显示了当地时间2008年6月6日02:21:33检测到的事件，这一事件未在当地地震目录中列出。该检测模板为2008年6月6日ML1.4地震事件，检测的平均CC值为0.792 1，推断震级为ML2.10。
表 1 不同阈值下的检测事件Table 1. The detected events under different thresholds阈值 检测事件个数 互相关系数＜0.2的事件占比 9MAD 109 822 21.15% 12MAD 21 213 2.20% 15MAD 8 643 2.13% 注：MAD为中值绝对偏差。 ![]()
图 2 基于图像处理器加速的匹配定位技术检测结果示例 （起始时间为2008−06−06 00:00:00）（a） 平均互相关系数CC值，红色虚线表示检测阈值为0.3，红色圆点表示检测到的地震事件；（b） 对应图（a）中所示窗口互相关系数的分布统计图；（c） 连续波形（灰色）与模板事件（红色）在时间窗口中的比较，台网台站和各分量CC值分别标记于波形左右两侧Figure 2. An example of detection based on the GPU-M&L technique（a） Mean CC value trace，where the red dashed line marks the detection threshold of 0.3，and the positive detection is shown by the red dot；（b） Histogram of the CC value distribution corresponding to the window shown in Fig. （a）；（c） Comparisons of continuous waveforms （gray） with templates （red） in a zoom-in time window. The stations of seismic networks and CC values are given to the left and right，respectively1.3 检测事件筛选
为进一步提高检测事件的真实性，本文对上述获得的检测事件进行了如下的筛选。首先，通过利用基于深度学习的去噪方法DeepDenoiser （Zhu et al，2019）对检测事件进行了去噪处理，并依据短时傅里叶变换对去噪后各分量时间域波形绘制了发震时刻后30 s的时频谱图，共得到36万7079幅时频谱图。然后，结合去噪后的时间域波形图，将时频谱图人工分类标记为“事件”和“噪声”两类标签，得到了分类后的事件和噪声的时频图各5 000幅，图3为去噪后事件和去噪后噪声的时频谱图示例。
![]()
图 3 去噪后地震事件（a，b）和去噪后噪声（c，d）的时频谱图Figure 3. Time-frequency spectra of events （a，b） and noises （c，d） after denoising将上述1万幅标记好的时频谱图以3 ∶ 1 ∶ 1的比例分为训练集、验证集和测试集，训练卷积神经网络模型，实现噪声和事件分类。卷积神经网络结构如图4所示，具体包括3个卷积层，卷积核个数分别为32，64，128，卷积核个数的增加保证了对图像特征的充分提取；卷积核大小均为3×3，移动步长均为1；每一个卷积层后增加Relu层进行非线性激活，并在卷积层与非线性部分之间使用批量归一化层来加速网络训练并降低对网络初始化的敏感度（Ioffe，Szegedy，2015）。3层池化层均采用最大池化的方法，池化窗口大小均为2×2，滑动窗口步长均为2；两个全连接层的隐藏单元数分别为512和2，分别使用Relu和Softmax激活函数进行激活，并在第一个全连接层前加入Dropout层防止过拟合，提升模型的泛化能力（Srivastava et al，2014）。最终得到的模型在测试集上的准确率为96%。将该模型应用于上述得到的检测事件的各分量时频谱图，并将被分类为“事件”的时频谱图总数大于9的检测事件筛选为一个真实的地震事件，最终得到了1万6 836个筛选后的地震事件。
1.4 检测事件精定位
小震活动包含了丰富的构造信息，其位置的精确程度直接关系到对研究区构造活动的认识。对于水库地震，震源深度也是判断水库触发地震的重要基础（姜金钟等，2016）。因此，本文再次按照与1.1节中相似的步骤和参数，利用深度学习算法PhaseNet对筛选后的1万6836个检测事件进行震相拾取，并利用波形互相关技术得到了互相关系数大于0.5的96万1899条P波和33万1331条S波互相关走时差数据。双差定位法可以同时使用绝对走时和互相关走时差数据进行定位，但当一个地震对在同一台站的绝对走时差和互相关走时差相差较大时，容易造成定位结果不稳定。因此，为了避免这种情况，本文以精度更高的互相关走时差数据为主进行双差定位，最终得到了1万5745个事件的高分辨率震源位置。
2. 结果分析
2.1 地震事件检测结果
本文以台网目录中经过重定位的578个事件作为模板事件，利用基于图像处理器加速的匹配定位技术进行事件检测，筛选后得到的事件高达1万6836个，约为台网地震目录的地震记录数量的13倍。检测目录的M-t图 （图5）显示检测出来的地震事件在2005年至2008年4月期间分布较为零散，集中分布在2008年5月之后的时段内。图6给出了台网地震目录与检测目录的震级-频率关系，可以看出，GPU-M&L识别出了大量台网地震目录中遗漏的地震事件，且为震级集中于ML2.0以下的微震事件。同时，本文使用地震活动性分析软件ZMAP （Wiemer，2001），基于最大曲率法来估算地震目录的完备震级Mc，相较于台网震级目录，基于GPU-M&L获得的检测目录使得完备震级从1.4降低到−0.1。
![]()
图 6 台网地震目录和检测目录的震级-频率关系图Figure 6. The magnitude-frequency relationship for the local catalog and detected catalog根据古登堡-里希特（Gutenberg-Richter，缩写为G-R）的震级-频次关系式（Gutenberg，Richter，1944），b值直接反映了大小地震之间的比例关系，也是用于分析水库触发地震的重要参数。实验室及大量的实际观测（Scholz，1968；Wiemer，Benoit，1996；Murru et al，2007）都充分表明了b值与应力呈负相关关系，b值越低，表明应力越高。与同震级天然构造地震相比，由于水库蓄水造成地下介质孔隙压力增大或者水的润滑作用使断层摩擦力减小等影响，水库触发地震的应力降值明显偏低，即b值较高（do Nascimento et al，2004；Mekkawi et al，2004）。小震级事件的遗漏可能会导致计算的b值偏低，大地震后的低b值也可能指示了地震台网无法检测到许多小震级的地震事件（Chang et al，2007）。因此，本文利用MPU-M&L补充了更多小震级事件的检测目录，统计自2005年9月蓄水以来研究区大于完备震级（Mc−0.1）的震级-频度关系（图7），并通过最大似然法（Aki，1965）计算得到b值及拟合误差为0.632±0.006，这一结果与Zhang和Zhou （2016）利用1985—2012年中国地震台网中心地震目录及刘雁冰和裴顺平（2017）利用1992—2014年龙门山地区全部M＞2.0地震目录计算所得的研究区所属区域b值相比较小，指示较高的区域应力，即研究区蓄水后的地震活动可能主要是构造应力积累导致的天然地震活动。
![]()
图 7 水库蓄水后的检测目录震级-频率关系图Figure 7. The magnitude-frequency relationship for the detected catalog after reservoir impoundment2.2 检测事件精定位结果
本文再次利用PhaseNet对筛选后的检测目录中地震事件进行P波、S波初动拾取，联合波形互相关技术使用相对定位法hypoDD对其进行精定位，定位结果水平和垂直方向的平均误差分别为42 m，51 m和38 m，发震时刻误差为7 ms，最终得到了1万5745个事件的高分辨率震源位置（图8），其水平方向误差分布如图9所示。重定位后的地震事件集中分布在北川—映秀断裂和汶川—茂汶断裂之间，延展分布趋势与区域内三条大断裂的北东走向基本一致。将重定位后的事件分别沿平行断裂走向方向（剖面AB）和垂直断裂走向方向（剖面CD）进行投影。剖面AB （图8a）显示，重定位后的事件震源集中分布在0—20 km深度范围内，且沿剖面方向距剖面A端约20—30 km位置处存在一个地震空区，与三维空间分布中（图10）绿色圆圈所指示的位置一致，对应于水库库体覆盖区以下的区域。剖面CD （图8b）显示，重定位后震源深度在0—5 km和10—15 km存在优势分布，两层大致平行，呈水平状。在深度7 km左右存在一地震空区，与龙门山地区三叠系雷口坡组和嘉陵江组主滑脱面位置一致（周斌等，2010）。重定位后检测目录的P波和S波走时分布如图11所示，其数据离散程度和整体趋势较为合理，但考虑到研究区范围过小且深度剖面上存在的空区较小，在走时分布图中未能呈现明显的与震源深度分布相对应的“两支现象”。
![]()
图 8 检测事件重定位后的震中分布（a）及沿剖面AB和CD方向的事件投影分布（b，c）F1：汶川—茂汶断裂； F2：北川—映秀断裂；F3：安县—灌县断裂Figure 8. The distribution of relocated earthquakes （a） and projection of the detected events along AB and CD （b，c）F1：Wenchuan-Maowen fault；F2：Beichuan-Yingxiu fault；F3：Anxian-Guanxian fault![]()
图 9 检测事件的水平方向重定位误差分布Figure 9. The distribution of horizontal relocation errors of the detected events![]()
图 10 检测事件重定位后的三维分布图Figure 10. 3-D distribution of the detected events after relocation![]()
图 11 检测目录重定位后的P波、S波走时分布图Figure 11. Travel time distribution of P wave and S wave of the detected events after relocation震源深度也是判断水库触发地震的重要参数，水库触发地震的震源深度一般不超10 km （马文涛等，2013）。杨智娴等（2004）对中国中西部地区（21°N—36°N，98°E—112°E）1992—1999年间的6 459次地震进行了重定位，精定位后的震源深度在0—10 km，0—15 km和0—20 km范围内的地震事件比例分别为52%，71%和82%。而Yang等（2005）利用双差定位法对中国中西部地区（21°N—36°N，98°E—112°E）1992—1999年间的1万零57次地震进行重定位的结果显示，平均震源深度为11 km，震源深度在0—15 km和0—20 km范围内的地震事件占地震总数的比例分别为77%和91%。统计本文双差定位结果显示，水库蓄水后库区平均震源深度为5.14 km，震源深度小于10 km的地震数占地震总数的80.9%。将其与天然构造地震深度进行对比（杨智娴等，2004；Yang et al，2005），指示库区内存在着水库触发地震活动。
3. 讨论与结论
汶川地震之后，为探究紫坪铺水库与汶川地震之间可能的关系，国内外许多研究者开展了有关紫坪铺水库蓄水后地震活动性的研究。Deng等（2010）认为根据库仑应力改变，紫坪铺水库的蓄水仅能造成水库附近地区浅层地震（＜5 km）的频率增加。Ruan等（2017）通过对汶川地震前5个月微震活动的研究，未发现当地地震活动性与紫坪铺水库水位之间可能的相关性。张致伟等（2009）和程万正等（2010）则认为，紫坪铺水库蓄水以来增加的地震活动从表面上看与水库的加卸载过程有关，但从长期序列来看，其频次和强度仍在区域地震活动的正常起伏范围之内。
本文基于紫坪铺水库地震台网提供的2005—2008年连续波形数据，利用GPU-M&L技术进行了地震检测。相较于前人的研究（Ruan et al，2017；Huang et al，2020；Liu et al，2020）， 本文为了确保检测事件的真实性，增加了对检测事件进行筛选的步骤。通过使用深度学习算法DeepDenoiser对检测事件进行去噪处理，并基于去噪后波形绘制的时频谱图，设计了一个卷积神经网络模型对其进行筛选，最终得到的地震目录事件增加到1万6836个，极大地丰富了原始地震目录，提高了地震活动性分析的可靠性。
将最终得到的检测目录与台网地震目录给出的重定位前后的震源分布进行对比，可以发现，台网地震目录重定位前（图12中灰色圆点）在平面和深度剖面上的分布较为分散，无明显的优势震源分布；重定位后的台网地震目录（图12中红色圆点）在平面分布上的特征与检测目录的结论类似，在沿断层走向和垂直断层走向两个方向的深度剖面中存在一个深度处于5—10 km范围的地震空区。检测目录中大量地震事件重定位后的结果将该地震空区的范围缩小，在沿断层走向深度剖面上勾勒出了更明显的水库库体覆盖区下方的地震空区，在垂直断层走向的深度剖面上将上述地震空区的范围缩小到7 km深度左右。
![]()
图 12 台网目录重定位前后震源分布（a）及沿剖面AB和CD方向的事件投影分布（b，c）F1：汶川—茂汶断裂； F2：北川—映秀断裂；F3：安县—灌县断裂Figure 12. The distribution of earthquakes before and after relocation （a） and projections of the events along AB and CD profiles，respectively （b，c）F1：Wenchuan-Maowen fault；F2：Beichuan-Yingxiu fault；F3：Anxian-Guanxian fault龙门山断裂带发育多层次的滑脱层，历史地震的优势震源深度可能反映龙门山断裂带的推覆逆冲面的铲式变化及其下伏的滑脱层位置（李勇等，2013）。滑脱层作为一个塑性软弱区域，易发生形变，且较难积累应变能，应力容易在其上部的脆性地层中积累（李大虎等，2021）。检测目录重定位后的结果显示，地震活动集中分布在0—20 km深度范围内脆性程度高、渗透稳定性低的岩层中，而7 km左右深度处为岩性较软弱、渗透稳定性高的三叠系雷口坡组和嘉陵江组白云岩和膏盐，少有地震发生，表现为一个明显的地震空区。依据本文的震源分布的优势深度，推测在研究区壳内20 km左右的深度存在一滑脱层顶面，与朱介寿（2008）对跨龙门山爆破地震剖面研究得出的松潘−甘孜地块及龙门山推覆体的中地壳内（20—30 km深）普遍存在一个厚约10 km的滑脱层的结果相一致。这一结果也证实了本文定位结果的可靠性。
前人研究总结出水库触发地震的震源深度一般不超过10 km，震中分布一般距库区不超过10 km，b值相较于区域构造地震的b值较高（胡毓良，陈献程，1979；丁原章等，1983；华卫等，2012；马文涛等，2013）。本文的重定位结果显示紫坪铺水库库体覆盖区下方的地震活动非常少，并在水库覆盖区下方呈现出明显的地震空区，与周斌等（2010）的研究结果一致。利用地震检测目录统计自2005年9月蓄水后研究区的震级−频度关系计算得到的b值结果显示，与前人计算的研究区所属区域b值相比本研究结果较小，指示了在研究时段内区域整体的平均应力水平较高。然而重定位结果显示，80%左右的地震震源深度不超过10 km，结合水库蓄水后研究区地震活动与水库加卸载过程有一定程度吻合的结果（程万正等，2010；张致伟等，2009；周斌等，2010），推测研究区蓄水后的地震主要是天然构造地震活动，同时混杂发生水库触发地震。
总之，本文利用紫坪铺水库地震台网7个台站2005—2008年的波形记录，采用基于图像处理器加速的匹配定位技术（GPU-M&L）对其遗漏的地震事件进行检测和筛选，最终得到的检测目录事件高达1万6 836个，约为台网目录记录的地震数量的13倍，并使完备震级从台网目录的1.4降低到−0.1。然后，通过再次使用结合波形互相关的双差定位法进行重定位，结果显示地震事件沿北东方向展布于北川—映秀断裂和汶川—茂汶断裂之间，水库覆盖区下方存在一个明显的地震空区，且地震活动集中分布在0—20 km深度范围内脆性程度高、渗透稳定性低的岩层中。结合研究区域地质构造背景，优势震源深度较好地指示了研究区地壳内滑脱层位置。最后，根据检测目录计算的b值和震源深度分析结果推测研究区蓄水后的地震主要是构造应力积累导致的天然地震活动，并伴随水库触发地震混杂发生。
四川省地震局为本文提供了数据支持，香港中文大学刘敏博士为本文提供了GPU-M&L代码，加州理工学院朱尉强博士为本文提供了DeepDenoiser以及PhaseNet代码，审稿专家为本文提供了宝贵的修改意见及建议，本文大部分图件采用 GMT 绘制，作者在此一并表示感谢。
-
![]()
图 5 AR19 （a）和AR35 （b）台站的接收函数波形（上）和H-κ叠加结果（下）
Figure 5. Receiver functions （upper panels） and H-κ stacking results （lower panels） of the stations AR19 （a） and AR35 （b）
![]()
图 2 本文所用远震事件分布
红色三角形为研究区域中心，蓝色圆点为地震事件
Figure 2. Distribution of teleseismic events used in this study
Red triangle is the center of study area，blue dots are seismic events
![]()
图 4 地壳厚度H （a）和波速比κ （b）的误差分布
Figure 4. The error distribution of crustal thickness H （a） and vP/vS ratio κ （b）
![]()
图 6 研究区内地壳厚度H （a）与波速比κ （b）分布（白色圈为高波速比区域）
Figure 6. Distribution of crustal thickness H （a） and vP/vS ratio κ （b） beneath the stations in the study area （The white circles denote the areas with high vP/vS ratio）
![]()
图 7 地壳厚度H与高程E （a）、波速比κ （b）的关系图
Figure 7. Correlations between crustal thickness H and elevations E （a） and vP/vS ratios κ （b）
表 1 各台站基本信息及相应的H-κ叠加结果
Table 1 Basic information of each station and corresponding H-κ stacking results
台站 高程/m 接收函数
条数tPs /s 地壳厚度/km vP/vS 台站 高程/m 接收函数
条数tPs /s 地壳厚度/km vP/vS AR03 833 72 4.7 37.9±1.9 1.73±0.04 AR18 971 92 4.7 36.9±2.4 1.74±0.06 AR04 809 29 4.8 36.5±3.1 1.77±0.08 AR19 1129 47 4.8 34.0±1.6 1.83±0.04 AR05 481 10 4.7 34.4±2.1 1.79±0.05 AR22 663 43 4.8 34.5±2.3 1.81±0.06 AR06 450 47 4.4 33.9±2.4 1.76±0.06 AR24 878 35 4.8 35.5±2.5 1.77±0.06 AR07 529 27 4.4 34.4±2.6 1.75±0.06 AR27 986 87 4.7 36.5±2.9 1.76±0.07 AR08 588 42 4.5 35.5±2.9 1.75±0.07 AR28 1038 79 4.8 35.5±2.8 1.80±0.07 AR09 883 37 4.6 35.5±1.8 1.76±0.05 AR29 1108 56 5.0 36.0±3.4 1.80±0.09 AR10 935 27 4.7 36.5±2.6 1.76±0.06 AR30 1179 71 4.7 35.9±3.7 1.77±0.08 AR11 1124 28 4.9 36.9±1.8 1.77±0.04 AR31 1193 23 4.7 37.5±2.8 1.74±0.06 AR12 688 83 4.9 34.5±4.0 1.83±0.10 AR32 881 105 5.0 36.0±2.2 1.81±0.06 AR14 710 90 4.7 34.4±1.9 1.78±0.05 AR33 1042 94 4.7 36.6±2.8 1.77±0.06 AR16 1254 61 4.7 36.0±2.5 1.76±0.06 AR35 1071 18 4.7 37.0±1.4 1.74±0.04 AR17 816 43 4.6 35.9±1.7 1.77±0.04 
下载: 导出CSV
-
樊祺诚,赵勇伟,李大明,武颖,隋建立,郑德文. 2011. 大兴安岭哈拉哈河—绰尔河第四纪火山分期:K-Ar年代学与火山地质特征[J]. 岩石学报,27(10):2827–2832. Fan Q C,Zhao Y W,Li D M,Wu Y,Sui J L,Zheng D W. 2011. Studies on Quaternary volcanism stages of Halaha river and Chaoer river area in the Great Xing’an Range:Evidence from K-Ar dating and volcanic geology features[J]. Acta Petrologica Sinica,27(10):2827–2832 (in Chinese).
高延光,李永华. 2014. 中国东北—华北地区地壳厚度与泊松比及其地质意义[J]. 地球物理学报,57(3):847–857. Gao Y G,Li Y H. 2014. Crustal thickness and vP/vS in the Northeast China−North China region and its geological implication[J]. Chinese Journal of Geophysics,57(3):847–857 (in Chinese).
高占永. 2015. 中国东北地区地壳上地幔结构的接收函数研究[D]. 北京: 中国地震局地球物理研究所: 35–42. Gao Z Y. 2015. The Study of the Crustal and Upper Mantle Structure in Northeast China From Teleseismic Receiver Function[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 35–42 （in Chinese）.
韩江涛,王天琪,刘文玉,刘国兴,韩松,刘立家. 2018. 阿尔山火山群深部“拱桥式”岩浆系统及其稳定性分析[J]. 地震地质,40(3):590–610. Han J T,Wang T Q,Liu W Y,Liu G X,Han S,Liu L J. 2018. Deep “arch-bridge” magmatic system of the Aershan volcanic group and its stability analysis[J]. Seismology and Geology,40(3):590–610 (in Chinese).
嵇少丞,王茜,杨文采. 2009. 华北克拉通泊松比与地壳厚度的关系及其大地构造意义[J]. 地质学报,83(3):324–330. Ji S C,Wang Q,Yang W C. 2009. Correlation between crustal thickness and Poisson’s ratio in the North China Craton and its implication for lithospheric thinning[J]. Acta Geologica Sinica,83(3):324–330 (in Chinese).
李卫东, 颜雷, 刘俊斌, 于宝洵, 赵玉琢, 王迺宸, 马志通, 苏本聪, 李飞达, 李小鹏, 赵超, 张廷绣. 2014. 伊敏林场幅L51E001003: 1/5万区域矿产地质调查报告[DS]. 全国地质资料馆, doi: 10.35080/n01.c.135496. Li W D, Yan L, Liu J B, Yu B X, Zhao Y Z, Wang N C, Ma Z T, Su B C, Li F D, Li X P, Zhao C, Zhang T X. 2014. Yimin forest field area L51E001003: 1/50 000 regional mineral geological survey report[DS]. National Geological Data Museum, doi: 10.35080/n01.c.135496 （in Chinese）.
李英康,高锐,姚聿涛,米胜信,李文辉,熊小松,高建伟. 2014. 大兴安岭造山带及两侧盆地的地壳速度结构[J]. 地球物理学进展,29(1):73–83. Li Y K,Gao R,Yao Y T,Mi S X,Li W H,Xiong X S,Gao J W. 2014. The crust velocity structure of Da Hinggan Ling orogenic belt and the basins on both sides[J]. Progress in Geophysics,29(1):73–83 (in Chinese).
梁宏达,高锐,侯贺晟,金胜,韩江涛,韩松,刘国兴. 2016. 大兴安岭与两侧盆地结合地带深部电性结构与岩石圈尺度构造关系[J]. 地球物理学报,59(5):1696–1704. Liang H D,Gao R,Hou H S,Jin S,Han J T,Han S,Liu G X. 2016. Deep electrical structure beneath the Da Hinggan Ling and the junction zone with adjacent basins and their tectonic relationship at a lithospheric scale[J]. Chinese Journal of Geophysics,59(5):1696–1704 (in Chinese).
刘嘉麒. 1999. 中国火山[M]. 北京: 科学出版社: 13–28. Liu J Q. 1999. Volcano in China[M]. Beijing: Science Press: 13–28 （in Chinese）.
潘桂棠,肖庆辉,陆松年,邓晋福,冯益民,张克信,张智勇,王方国,邢光福,郝国杰,冯艳芳. 2009. 中国大地构造单元划分[J]. 中国地质,36(1):1–28. Pan G T,Xiao Q H,Lu S N,Deng J F,Feng Y M,Zhang K X,Zhang Z Y,Wang F G,Xing G F,Hao G J,Feng Y F. 2009. Subdivision of tectonic units in China[J]. Geology in China,36(1):1–28 (in Chinese).
潘佳铁,吴庆举,李永华,余大新. 2014. 中国东北地区噪声层析成像[J]. 地球物理学报,57(3):812–821. Pan J T,Wu Q J,Li Y H,Yu D X. 2014. Ambient noise tomography in Northeast China[J]. Chinese Journal of Geophysics,57(3):812–821 (in Chinese).
汤吉,王继军,陈小斌,赵国泽,詹燕. 2005. 阿尔山火山区地壳上地幔电性结构初探[J]. 地球物理学报,48(1):196–202. Tang J,Wang J J,Chen X B,Zhao G Z,Zhan Y. 2005. Preliminary investigation for electric conductivity structure of the crust and upper mantle beneath Aershan volcano area[J]. Chinese Journal of Geophysics,48(1):196–202 (in Chinese).
谢振新,吴庆举,周仕勇,朱敏. 2018. 兴蒙造山带诺敏河火山群地壳厚度与波速比研究[J]. 地球物理学报,61(12):4805–4816. Xie Z X,Wu Q J,Zhou S Y,Zhu M. 2018. Study of crustal thickness and vP/vS ratio beneath the Nuomin River volcanoes[J]. Chinese Journal of Geophysics,61(12):4805–4816 (in Chinese).
许文良,孙晨阳,唐杰,栾金鹏,王枫. 2019. 兴蒙造山带的基底属性与构造演化过程[J]. 地球科学:中国地质大学学报,44(5):1620–1646. Xu W L,Sun C Y,Tang J,Luan J P,Wang F. 2019. Basement nature and tectonic evolution of the Xing’an-Mongolian orogenic belt[J]. Earth Science:Journal of China University of Geosciences,44(5):1620–1646 (in Chinese).
张风雪,吴庆举,李永华. 2013. 中国东北地区远震P波走时层析成像研究[J]. 地球物理学报,56(8):2690–2700. Zhang F X,Wu Q J,Li Y H. 2013. The traveltime tomography study by teleseismic P wave data in the Northeast China area[J]. Chinese Journal of Geophysics,56(8):2690–2700 (in Chinese).
张风雪,吴庆举,李永华. 2014. 中国东北地区远震S波走时层析成像研究[J]. 地球物理学报,57(1):88–101. Zhang F X,Wu Q J,Li Y H. 2014. A traveltime tomography study by teleseismic S wave data in the Northeast China area[J]. Chinese Journal of Geophysics,57(1):88–101 (in Chinese).
张广成. 2012. 用接收函数方法研究中国东北地区地壳上地幔结构[D]. 北京: 中国地震局地球物理研究所: 45–46. Zhang G C. 2012. The Structure of the Crust and Upper Mantle in Northeast China From Teleseismic Receiver Function[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 45–46 （in Chinese）.
张广成,吴庆举,潘佳铁,张风雪,余大新. 2013. 利用H-κ叠加方法和CCP叠加方法研究中国东北地区地壳结构与泊松比[J]. 地球物理学报,56(12):4084–4094. doi: 10.6038/cjg20131213 Zhang G C,Wu Q J,Pan J T,Zhang F X,Yu D X. 2013. Study of crustal structure and Poisson ratio of NE China by H-κ stack and CCP stack methods[J]. Chinese Journal of Geophysics,56(12):4084–4094 (in Chinese).
赵海玲,邓晋福,陈发景,胡泉,赵世柯. 1996. 东北地区新生代火山作用、深部作用与大陆裂谷型盆地[J]. 地球科学,21(6):51–55. Zhao H L,Deng J F,Chen F J,Hu Q,Zhao S K. 1996. Cenozoic volcanism,deep interior processes and continental rift basin formation in the northeastern China[J]. Earth Science:Journal of China University of Geosciences,21(6):51–55 (in Chinese).
赵勇伟,樊祺诚. 2012. 大兴安岭哈拉哈河—绰尔河第四纪火山岩地幔源区与岩浆成因[J]. 岩石学报,28(4):1119–1129. Zhao Y W,Fan Q C. 2012. Mantle sources and magma genesis of Quaternary volcanic rocks in the Halaha river and Chaoer river area,Great Xing’an Range[J]. Acta Petrologica Sinica,28(4):1119–1129 (in Chinese).
朱洪翔,田有,刘财,冯晅,杨宝俊,刘才华,刘廷,马锦程. 2017. 中国东北地区高分辨率地壳结构:远震接收函数[J]. 地球物理学报,60(5):1676–1689. Zhu H X,Tian Y,Liu C,Feng X,Yang B J,Liu C H,Liu T,Ma J C. 2017. High-resolution crustal structure of Northeast China revealed by teleseismic receiver functions[J]. Chinese Journal of Geophysics,60(5):1676–1689 (in Chinese).
朱勤文,路凤香,谢意红,郑建平. 1997. 大陆边缘扩张型活动带火山岩组合:松辽盆地周边中生代火山岩研究[J]. 岩石学报,13(4):551–562. doi: 10.3321/j.issn:1000-0569.1997.04.009 Zhu Q W,Lu F X,Xie Y H,Zheng J P. 1997. Volcanic rocks assemblages in active belt of spreading type in continental margin:Study on Mesozoic volcanic rocks around Songliao basin[J]. Acta Petrologica Sinica,13(4):551–562 (in Chinese).
Guo Z,Chen Y J,Ning J Y,Feng Y G,Grand S P,Niu F L,Kawakatsu H,Tanaka S,Obayashi M,Ni J. 2015. High resolution 3-D crustal structure beneath NE China from joint inversion of ambient noise and receiver functions using NECESSArray data[J]. Earth Planet Sci Lett,416:1–11. doi: 10.1016/j.jpgl.2015.01.044
Guo Z,Chen Y J,Ning J Y,Yang Y J,Afonso J C,Tang Y C. 2016. Seismic evidence of on-going sublithosphere upper mantle convection for intra-plate volcanism in Northeast China[J]. Earth Planet Sci Lett,433:31–43. doi: 10.1016/j.jpgl.2015.09.035
Herrmann R B. 2013. Computer programs in seismology:An evolving tool for instruction and research[J]. Seismol Res Lett,84(6):1081–1088. doi: 10.1785/0220110096
Kikuchi M,Kanamori H. 1982. Inversion of complex body waves[J]. Bull Seismol Soc Am,72(2):491–506.
Ligorría J P,Ammon C J. 1999. Iterative deconvolution and receiver-function estimation[J]. Bull Seismol Soc Am,89(5):1395–1400. doi: 10.1785/BSSA0890051395
Liu J Q,Han J T,Fyfe W S. 2001. Cenozoic episodic volcanism and continental rifting in northeast China and possible link to Japan Sea development as revealed from K-Ar geochronology[J]. Tectonophysics,339(3/4):385–401.
Liu Y N,Niu F L,Chen M,Yang W C. 2017. 3-D crustal and uppermost mantle structure beneath NE China revealed by ambient noise adjoint tomography[J]. Earth Planet Sci Lett,461:20–29. doi: 10.1016/j.jpgl.2016.12.029
Tao K,Niu F L,Ning J Y,Chen Y J,Grand S,Kawakatsu H,Tanaka S,Obayashi M,Ni J. 2014. Crustal structure beneath NE China imaged by NECESSArray receiver function data[J]. Earth Planet Sci Lett,398:48–57. doi: 10.1016/j.jpgl.2014.04.043
Tarkov A P, Vavakin V V. 1982. Poisson’s ratio behaviour in crystalline rocks: Application to the study of the Earth’s interior[J]. Physics Earth Planet Inter, 29（1）: 24-29.
Wessel P, Smith W H F, Scharroo R, Luis J, Wobbe F. 2013. Generic mapping tools: Improved version released[J]. Eos Trans Am Geophys Union, 94（45）: 409–410.
Zandt G,Ammon C J. 1995. Continental crust composition constrained by measurements of crustal Poisson’s ratio[J]. Nature,374(6518):152–154. doi: 10.1038/374152a0
Zhu L P,Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions[J]. J Geophys Res:Solid Earth,105(B2):2969–2980. doi: 10.1029/1999JB900322
-
期刊类型引用(1)
1. 黄春梅，吴朋，苏金蓉，王宇航，魏娅玲，李大虎，颜利君. 2017年8月九寨沟M_S7.0地震序列S波分裂特征. 华北地震科学. 2020(01): 29-37 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 326
- HTML全文浏览量: 171
- PDF下载量: 113
- 被引次数: 2
