Three-dimensional seismic velocity structure beneath the M6.4 Yangbi，Yunnan earthquake region
-
摘要: 利用2015年1月至2021年5月28日期间我国云南省漾濞县及周边地区固定台站和漾濞地震后布设的流动台站所记录到的近震资料，使用双差层析成像方法获得了该地震震区的高分辨率地壳三维速度结构和震源位置。重定位结果显示，漾濞M6.4地震序列主要沿NW−SE向展布，与维西—乔后—巍山断裂走向一致，地震主要集中在4—10 km的深度范围，呈约80°高倾角分布。结合定位结果与三维速度结构显示：漾濞M6.4地震序列的空间分布与速度结构变化具有相关性，主震位于P波、S波高低速异常交界处，这种介质物性变化的交界地带可能有利于中强地震的孕育和发生，余震主要分布在低P波速度、高S波速度和低波速比的脆性区域；沿漾濞地震序列的分布走向，主震两侧呈现完全不同的速度结构，其西北部具有明显的高P波速度、低S波速度特征，该地区高密度、强韧性的地层可能是阻挡漾濞地震的NW向破裂而呈单向破裂特征的原因。Abstract: This paper collected the seismic travel-time data both from temporary stations employed after Yangbi M6.4 earthquake and the seismic networks in Yangbi region and its vicinity from January of 2015 to 28 May, 2021 and performed the high-resolution inversion for three-dimensional velocity structure and accurate hypocentral locations by using double-difference seismic tomography method. The relocation results show that the sequence spread in the NW−SE direction along the Weixi-Qiaohou-Weishan fault. The focal depth in the area is generally in the range of 2−5 km with high dip angle of 80°. The three-dimensional velocity structure shows that the spatial distribution characteristics of the Yangbi M6.4 earthquake sequence are closely related to the velocity structure. The epicenter of the Yangbi M6.4 main shock was located near the P- and S-wave high-to-low-velocity anomaly transitional zones, which are favorable for occurrence of moderate-strong earthquakes, and the aftershocks are mainly distributed in the brittle region with low vP, high vS and low vP/vS. In addition, along the strike of Yangbi earthquake sequence, there are totally different velocity structures on both sides of main earthquake, and significant higher vP and lower vS anomalies are observed in the northwest of the Yangbi M6.4 main shock compared with the southeast part, which may obstruct the northwestward slipping of the seismogenic fault of Yangbi M6.4 earthquake, leading to the striking unilateral source rupture.
-
引言
地震计是地震动波形数据采集系统中最重要的部分，不同地震计之间的对比研究是很多研究人员感兴趣的课题。目前国内用于地震观测的地震计有多种（刘瑞丰等，2008），这些地震计所记录数据之间的一致性值得深入研究。
一般来说，地震计的记录特性可以用传递函数来表述，因此，仪器参数的对比基本上可以反映出不同地震计之间的性能差别。但不同地震计的实际观测差别却不容易通过仪器参数的对比直观地看出或者定量地得到。地震计记录的信号包括地震、爆破信号、地脉动信号、环境噪声（气流、海浪、人文等噪声）、环境参数（气压、温度等）变化引起的噪声以及地震计内部噪声，有时环境参数改变引起的噪声也被划归为仪器噪声。地震计之间的实际观测差异一方面来自于仪器噪声的差异，另一方面，因为对相同外部信号记录的不一致性，在实际观测中这些差异耦合在一起，使不同地震计在观测频带范围内呈现出各自的记录特性，从而可能影响到波形分析结果。例如，对于长周期地震计来说，低频范围内的数据常被用于基于接收函数、噪声成像等方法的地球自由振荡研究，而这一频率范围内的记录由于常受到气流、温度和气压的变化及地倾斜等因素的影响，导致不同地震计的记录存在一定的差别。在数据处理时，如果对不同地震计的观测差别有直观的、量化的认识，则可以更合理地选择适宜的数据处理方法。
本文拟对比两种宽频带地震仪器所记录波形的噪声概率密度函数（probability density functions，简写为PDFs），以期定量地观察不同仪器观测性能之间的差异，以便深入理解仪器参数与观测性能之间的关系。
1. PDF方法
1.1 PDF简介
PDF统计法是由Mcnamara和Buland （2004）于2004年首次提出，是基于传统的噪声功率谱密度（power spectral density，简写为PSD）分析方法，利用长时间内完整连续的波形记录分析地震观测台站噪声水平的方法。该方法区别于传统PSD方法之处在于，其在计算中不需要排除包括地震在内的突发事件，也不必刻意地选取外界噪声较小的时段来计算PSD，而是对所有记录数据进行相同处理，因此可以在保持数据连续性的同时，将各种对背景噪声的影响体现在PDF的概率值中，使得该方法包含了地震计对所有地表振动的记录，全面地反映出台站的噪声水平。这一方法近年来已被越来越多的国际同行应用，在美国该方法已被应用于美国地震学合作研究会（Incorporated Research Institutions for Seismology，简写为IRIS）和美国国家地震监测台网系统（Advanced National Seismic System，简写为ANSS）的数据管理和流动台阵观测数据质量控制（Mcnamara，Boaz，2005；Diaz et al，2010）。国内，PDF方法越来越多地被用来评估测震台站数据的观测质量，例如：吴建平等（2012）利用PDF方法分析了华北地区的背景噪声，得到了华北不同地区的背景噪声特征，并提出基于PDF分析的环境噪声特征可以对流动台阵的数据观测质量进行定量评估；葛洪魁等（2013）通过对不同台基的台站背景噪声进行PDF分析，提出了流动台阵观测的有效降噪方法。
PDF的计算方法可概述为：首先将一个较长时段的数据记录分成连续的若干段，以方便对每一段数据进行PSD计算；然后将所有的PSD计算结果进行概率统计，即统计出在某一频率，各功率谱密度值在所有可能出现的功率谱密度值中所占的比例，其计算方法为
$P\left( {f,d} \right) {\text{＝}} \frac{{{N_{df}}}}{{{N_f}}}{\text{，}}$
(1) 式中：P（f，d）表示频率为f、功率谱密度为d时的概率密度；Ndf表式频率为f时，功率谱密度在[d，d＋δd] dB范围内出现的次数；Nf表示频率为f时，各功率谱密度值出现的总次数。
1.2 数据处理过程
具体计算步骤为：① 用SAC软件（黄金莉，顾小虹，2001）去除地震波形记录的均值和线性趋势，再去除仪器响应，得到速度记录（单位：m/s）；② 将每小时的波形数据分为14段，每段数据的长度为1 000 s （即数据段重合率为80%），用以减小PSD的方差；③ 对每一段数据进行汉宁（Hanning）窗加窗处理，以压制PSD计算时的旁瓣影响，并利用功率谱修正周期图法进行PSD计算；④ 对每小时14段PSD结果进行平均，得到该时段的PSD均值；⑤ 对PSD均值结果进行1/10倍频程滤波，得到对数坐标下均匀平滑分布的速度PSD；⑥ 将各频率对应的PSD值转换为以dB为单位，统计各频率PSD的概率分布，得到两套仪器所记录波形的PDF结果；⑦ 绘制PDF结果图，并与全球地震背景噪声模型（Peterson，1993）的新低噪声模型（new low noise model，简写为NLNM）和新高噪声模型（new high noise model，简写为NHNM）进行对比，以便直观地了解两套仪器在各频段所记录噪声的分布情况。这两个模型是Peterson在定量地分析了全球75个固定地震台站的背景噪声功率谱密度后得到的。
在实际资料处理过程中，本文选取了两套仪器自2011年11月1日至2012年10月31日一整年的记录波形，按照上述数据处理过程对数据进行处理，并绘制PDF结果。
2. 利用PDF对比地震计特性
地震计在工作中会受到各种因素的干扰和影响，两台地震计在同一观测环境下记录到的波形也会存在一定的差别，这种差别可以体现出地震计不同的观测性能，还可以用来评估某一地震计是否能够维持正常工作。一般来说，不同类型、不同厂家的地震计，或者同一类型、不同工作年限的地震计，其记录波形均可能出现一定的差别。仪器参数的对比往往不能够细致地体现出这些差别，而PDF对比则可以给出定量的对比结果，这种对比结果可以帮助我们了解观测环境对不同地震计的影响程度，进而对观测数据进行评估，并可根据不同的观测条件和需求选择相应的地震计。
2.1 资料
因为宽频地震观测实验和壳幔速度研究项目“深部探测技术与实验研究”（SinoProbe）的开展，中国地质科学院地质研究所和中国地震局兰州地震研究所在甘肃省合作市测震台网的固定台站上架设了宽频带流动观测仪器，架设周期为2011年10月至2012年12月。在此期间，固定台和流动台的仪器并行观测，二者均位于基岩山洞内，且放置在同一水泥摆墩上，均采取了防风保温等措施，观测环境非常一致，良好的观测条件保证了其观测数据的质量。表1列出了观测仪器的主要参数，可以看出，这两套仪器均适合于长周期观测，且其记录特性优良；两套设备从仪器参数到观测参数虽然存在一定的差别，但这些差别恰好可以帮助我们衡量PDF方法能否体现出两套仪器观测性能的细微差别，从而体现出仪器参数不能提供的信息，进而对仪器观测性能有更为具体的认识。
表 1 观测仪器的主要参数Table 1. Main parameters of the observation instruments台站类型 数采型号 数采字长/bit 数采动态范围/dB 地震计型号 地震计频率范围 地震计动态范围/dB 采样率/sps 固定台 EDAS-24 24 ＞135 BBVS-60 40 Hz—60 s ＞130 100 流动台 REFTEK-130 24 ＞135 CMG-3T 50 Hz—120 s ＞140 50 2.2 仪器波形记录
为确保两套仪器正常工作，我们选取一次远震和一次近震的波形记录分别进行对比，对所有波形记录均去中值和去仪器响应。考虑到仪器可能存在方位角不完全一致的问题（陈继峰等，2016），为避免这种不一致带来的水平分量的波形差异，我们先将两台地震计各自的东西分量和南北分量合并为不考虑矢量方向的水平向记录。从两次地震的原始记录波形对比中可以看出，两套仪器所记录的波形走时曲线和幅值的一致性均非常高，差别细微，各震相的相位无明显差别（图1）。
![]()
图 1 流动台和固定台不同分量的地震波形记录对比（a）和（b）分别为远震的垂直分量和水平合并分量； （c）和（d）分别为近震的垂直分量和水平合并分量Figure 1. Comparisons of seismic waveforms for different component between the mobile station and the fixed station（a） and （b） are vertical component and synthesized horizontal component of distant earthquake, respectively；（c） and （d） are vertical component and synthesized horizontal component of near earthquake, respectively进一步对上述两套仪器的近震和远震波形记录进行互相关计算，在计算前将所有波形均进行0.02—20 Hz的带通滤波处理。近震波形的长度为20 s，选取相关窗长为1 s；远震波形的长度为2 400 s，选取相关窗长为10 s。从互相关结果（图2）中可以看出：无论垂直分量还是水平分量，近震波形在S波到达之前均有一定的差异，而在S波到达之后差异较小，这种差异可能是由于地震波传播过程中微小的路径差异所造成的；对于远震记录，其在面波部分存在一定差异，水平向的差异较垂直向小。通过近震和远震波形资料的对比可以看出，两套仪器运行正常，波形记录效果均很好。
![]()
图 2 流动台和固定台不同分量地震波形记录的相关分析结果（a）和（b）分别为远震垂直分量和水平合并分量； （c）和（d）分别为近震垂直分量和水平合并分量Figure 2. Cross-correlations of different component seismic waveforms of the mobile station and the fixed station（a） and （b） are vertical and combined horizontal component of distant earthquake respectively；（c） and （d） are vertical and combined horizontal component in near earthquake respectively2.3 PDF对比分析
由于两套观测仪器的观测环境非常一致，因此用PDF方法可以去除单次波形记录对比中存在的偶然性，从而可以看出两套仪器在对比频率范围内的记录差异。图3为利用两套仪器在相同观测时段（2011年11月1日至2012年10月31日）的波形数据绘制的PDF对比图，周期范围介于20 Hz—50 s之间，以确保两套仪器均处于速度平坦响应的范围内。
![]()
图 3 两套仪器不同分量的PDF对比（a）和（b）分别为流动台和固定台的垂直分量PDF；（c）和（d）分别为流动台和固定台的南北分量PDFFigure 3. The comparisons of PDF of different components（a） and （b） show vertical component PDFs of mobile station and fixed station； （c） and （d） show north-south component PDFs of mobile station and fixed station图3中对各台站的PDF结果均计算了均值和中值，可以看出，中值更接近实际的功率谱概率密度最大值。对比两个台站的垂直分量和南北分量的PDF中值，同时对两套仪器的仪器噪声（Ringler，Hutt，2010）进行对比（图4），可以看出：两套仪器的PDF中值相似度非常高，仅在周期小于0.1 s （10 Hz）以及大于10 s范围内略微有差别；两套设备的仪器噪声均非常低，BBVS-60的仪器噪声相对更低，两套仪器的环境噪声在全频段均远远高于仪器噪声。
![]()
图 4 流动台和固定台垂直分量（a）和南北分量（b）的PDF中值及仪器自噪声对比Figure 4. The comparison of the PDF median and instrument self noise between the mobile station and the fixed station in vertical component （a） and north-south component （b）2.4 PDF差异及讨论
地震仪器记录的信息包括地面运动信息和仪器噪声，在PDF对比分析时，也可以采用“环境背景噪声（包括环境噪声、干扰和地震事件）＋仪器噪声”的模型来进行解释。在地震计速度平坦输出的频率范围内，仪器所记录地震动的振幅远大于仪器噪声；在速度平坦频率范围之外，这种差别将明显减小（Mark et al，1990）。在本文计算PDF所涉及的频率范围内，仪器噪声也远低于背景噪声，因此噪声功率谱得到的结果主要反映了仪器对环境背景噪声的记录情况。
大体上，可以将仪器记录的背景噪声的频率范围分为高频（1 Hz以上）、低频（1—10 s）和长周期（10—50 s） 3个部分。前人的研究表明：高频部分的干扰主要以人类活动等为主，这种干扰随距离的衰减最快，另外气流以及地方震等自然因素也会形成高频干扰（Mcnamara，Buland，2004）；低频部分的干扰主要来自微震噪声，微震噪声与海底压力扰动和气流及气流引起的植物根茎扰动等有关，其主要峰值的形成一般认为与海洋波与海底或海岸线的非线性相互作用引起的海底压力扰动有关（Longuet-Higgins，1950；Frontera et al，2010 ），此外中远震也会影响到这一频率范围的背景噪声；长周期部分的噪声则主要与自然因素相关（Bonnefoy-Claudet et al，2006 ），如风、急流（瀑布和河流）、温度变化、气压变化和地倾斜等，其中地倾斜使重力耦合至水平分量中，因此水平分量的长周期噪声会高于垂直分量（Wielandt，2002）；此外各类远震波也在低频和长周期范围影响背景噪声。
从两种仪器的背景噪声PDF中值和仪器噪声对比来看：当周期在0.1 s以上时，两者的仪器噪声差异超过了背景噪声差异；当周期处于10 s以上的低频范围时，两种仪器的背景噪声PDF中值出现一定的差别，即固定台站的中值相对略低，这一频率范围内PDF概率密度分布则出现较为明显的差别，流动台站的概率密度分布相对较为聚拢；在周期小于0.1 s的频率范围内，二者的仪器噪声水平非常相近，但固定台站的PDF中值相对较低，二者的PDF概率密度分布并未出现明显差别，因此可以推断，这种差别源于仪器记录背景噪声幅值的差异。
为了更好地观察不同频率处两套仪器记录背景噪声的差别，我们以垂直分量为例，观察PDF中0.06 （16.7 Hz），0.1 （10 Hz），1，10，30和40 s等频率处的概率分布细节（图5）。可以看出：① 在0.06 s处，两套仪器的PSD分布的整体形态一致，但二者的PSD值存在一定偏差，说明两套仪器在这一频段记录的环境噪声是一致的。由于两套设备的仪器噪声大小在这一频率处相当，因此，PSD分布的偏移则反映出两套观测仪器在该频率处的灵敏度的实际细微差异，这一差异从仪器参数中并不能直接看出；② 在0.1，1，10 s等3处，两套仪器的PSD分布几乎完全一致，说明两套仪器的记录特性在这一频率范围内的一致性非常高。反观二者的仪器噪声则存在一定差异，考虑到仪器噪声的幅值远小于背景噪声的幅值，说明仪器噪声的差异并不影响到背景噪声的记录；③ 在30 s和40 s处，两套仪器的仪器噪声存在差异，固定台的相对较低。背景噪声也存在一定差异，不过，相较于背景噪声PDF中值（接近PSD峰值）的差别，其PSD分布差别较大，主要体现在PSD分布形态，固定台的分布更宽，流动台则较为聚拢。结合仪器噪声幅值，可以看出在这一频率，虽然固定台的仪器有更低的仪器噪声，但其对于外界因素的影响更为“敏感”，从而表现出更大的PSD变化范围，而流动台站的PSD值分布则较为集中。
![]()
图 5 不同频率处流动台和固定台PSD概率分布对比（a）—（f）分别为0.06，0.1，1，10，30和40 s处固定台和流动台PSD概率分布的对比Figure 5. PSD probability distribution comparison between the mobile station and the fixed station at certain frequencies（a）−（f） are the comparisons at frequencies of 0.06，0.1，1，10，30 and 40 s，respectively在地震观测系统中，数据采集器对于记录背景噪声的影响很小，而本文中两台数据采集器的工作参数也非常接近：均为24位高精度AD转换，动态范围大；工作温度均介于−20—60 ℃之间；GPS精度均小于±100 μs；数据采集器增益为1时，自身工作引起的噪声水平均约为1个计数（1个计数对应地动速度小于0.001 6 μm/s），因此，数据采集器对于PDF结果的影响极小，两套仪器记录的差异主要来自于地震计。地震计是受外界干扰影响的主要部分，诸如温度、气流、地倾斜等干扰因素对于地震观测系统的影响主要作用于地震计。如前文所述，地震计的仪器噪声大致可以分为地震计内部噪声和仪器受环境参数（气压、温度等）变化引起的噪声。对于速度地震计来说，地震计内部噪声主要包括换能器放大器的电子噪声和机械摆体的布朗热噪声（崔庆谷，2003）。本文选取的两种地震计均采用电容位移换能器，这是由于相较于速度换能器，位移换能器在噪声抑制方面有优良的性能（崔庆谷，2003），因此，两种地震计在换能器电子抑制噪声方面均有很好的效果。机械摆体的布朗热噪声大小由8πKbTe/QT0ms决定，其中Kb为玻尔兹曼常数，Te为开尔文温度，Q为摆体的品质因子，T0为摆体的自振周期，ms为摆体的等效质量。Q，T0和ms共同决定了机械摆体的布朗热噪声的大小。相对于地震计的内部噪声，环境参数变化引起的噪声往往要大得多，摆体受气压变化引起的噪声振幅有时比地面震动的振幅还要高几个数量级（Wielandt，Streckeisen，1982）；另外，虽然诸如地震的计恒温观测系统在不断改进（孙宏志，2016），但温度和地倾斜引起的噪声同样也是不容忽视的（崔庆谷，2003），这些环境参数变化引起的噪声也与摆体的Q，T0和ms等地震计参数有关。水平向的摆体特别容易受到地面形变和地倾斜的影响（段天山，袁顺，2011），垂直方向的摆体由于处在非真空环境中，更容易受到气压的影响。在30—40 s频率范围内，气压、温度变化等正是影响观测的主要因素，因此，两种仪器的PDF在低频范围的差别与各自不同的Q，T0和ms等地震计参数有很大关系。结合图5e，f可以看出，固定台BBVS-60型地震计在30—50 s范围内受自然因素的干扰而引起的背景噪声变化范围较流动台CMG-3T型地震计更大，CMG-3T型地震计在此频率范围内的记录变化范围较为集中。
从两台地震计零极点描述的传递函数来看，BBVS-60型地震计给出的零点为0，0和−200.5，极点p1，p2，···，p6分别为−0.074 5＋0.074 5i，−0.074 5−0.074 5i，−100.65＋264.63i，−100.65−264.63i，−134.3，−493.7。CMG-3T型地震计给出的零点为0和0，极点为p1＝−0.037 01＋0.037 008i，p2＝−0.037 01−0.037 008i。根据各极点的具体值可知，p1和p2为闭环主导极点，系统特性主要由零点及p1和p2组成的二阶系统描述，因此，二者的传递函数均可表示为
$H\left( s \right) {\text{＝}}\frac{{k{s^2}}}{{\left( {s - p_1} \right)\left( {s - p_2} \right)}} {\text{＝}} \frac{{k{s^2}}}{{{s^2} {\text{＋}} {{{K}}_{11}}s {\text{＋}} {{{K}}_{12}}}},$
(2) 式中，s为复变量，k为灵敏度和归一化因子，K11和K12为频率常数，K11与阻尼及闭环周期相关，K12与闭环周期相关。闭环主导极点主要反映了地震计的低频特性（蔡亚先等，2004），一组共轭复数闭环主导极点的负实部反应了地震计低频振荡的衰减特性。式（2）中两台地震计极点值的实部和虚部均相差近一倍，这一差别可能导致两台地震计对于噪声响应的衰减不同，也可能导致PDF低频区域的差别，同时不同的极点也决定了二者有不同的闭环周期。
为了定量地得到两种地震计在本文分析频率范围内的记录差异，本文提出一种新的PDF对比方法，即将两个台站的PDF结果进行差值运算，得到的结果可以显示出二者在各频率值的PSD分布及PDF中值差异（图6）。从图6中的PDF中值差曲线可以明显地看出两套仪器背景噪声的最大概率差值，根据色条可以更为细致地得到不同频率范围内两套仪器背景噪声的具体差值。总体来看，两套仪器在周期20 s以下的记录性能，除了低于0.1 s的周期范围有一定的中值差之外，总体非常一致。周期20 s以上差异则较为明显，从图6中可以直观地、量化地得到两种仪器背景噪声的具体差值。
![]()
图 6 流动台和固定台的PDF分布定量差值Figure 6. The quantitative difference of PDFs between the mobile station and the fixed station3. 讨论与结论
PDF是了解地震观测台站环境噪声水平的方法，也是评估台站记录波形质量的重要方法，它不仅能够提供环境噪声评估，而且从本文的对比分析来看，PDF方法可以有效地判断出不同观测仪器在同样观测条件下的细微差别。这种差别来自于仪器的结构和具体工作参数，诸如仪器频率工作范围、灵敏度、零极点等参数可以直观地从仪器参数中得到，但是仪器在一定观测环境下的具体观测性能却不容易从仪器参数中量化得出，因此利用PDF分析方法，我们可以对不同仪器观测性能的差别给出量化的、直观的结果，并可结合仪器的具体参数，更为细致地了解不同地震计的差异，进而对观测资料作出更为合理的评估。尽管如此，但是由于条件所限，本研究中采用了不同的数据采集器，虽然对于观测及分析结果影响很小，但在今后的研究中，可以避免这一类的影响因素以得到更为合理的结果。此外，人为地改变温度、湿度、气流等环境观测参数，能得到更为全面的不同地震计的性能差异。
几位审稿专家对本文提出了诸多宝贵意见，作者在此表示感谢！
-
![]()
图 9 不同深度z的vP （左）和vS （右）分布
黑色圆圈代表地震， 图（b）中AA′和BB′分别为平行、垂直于漾濞地震序列走向的剖面
Figure 9. vP （left） and vS （right） distribution at the depth z of 5，10 and 15 km
The circles represent earthquakes，and in Fig. （b） AA′ and BB′ are the two profiles parallel to and perpendicular with the strike of the Yangbi sequence，respectively
![]()
图 1 漾濞地区构造背景、台站和2010年以来M4.0以上地震分布图
Figure 1. Tectonic settings，seismic stations and M≥4.0 earthquakes since 2010 in Yangbi region
![]()
图 2 射线路径、台站及地震分布
Figure 2. Geographical distribution of seismic stations （triangles），earthquake hypocenters （circles） and ray paths （blue lines） considered in this study
![]()
图 4 本文选取的最佳阻尼因子（a）和平滑因子（b）的折中曲线图
Figure 4. Trade-off curves for determining optimal damping （a） and smoothing （b） parameters
![]()
图 5 深度为5 km （a），10 km （b）和15 km （c）时P波（左）和S波（右）检测板的测试结果（0.1°×0.1°）
Figure 5. P （left） and S （right） wave checkboard resolution tests at the depth of 5 km （a），10 km （b） and 15 km （c） with the grid 0.1°×0.1°
![]()
图 6 漾濞M6.4地震序列重定位前（a）、后（b）的空间分布变化
Figure 6. Spatial distribution variation of the Yangbi M6.4 earthquake sequence before （a） and after （b） relocation
![]()
图 7 漾濞M6.4地震序列重定位前（a）、后（b）的震源深度变化
Figure 7. Focal depth variation of the Yangbi M6.4 earthquake sequence before （a） and after （b） relocation
![]()
图 8 重定位前（红）、后（蓝）走时残差直方图
Figure 8. Histograms of travel time residual before （red） and after （blue） relocation
![]()
图 11 vP，vS和vP/vS沿纵向剖面的变化及漾濞地震序列震中分布图
Figure 11. vP，vS and vP/vS images along different vertical profiles and epicenter distribution of Yangbi earthquake sequence
表 1 本文所用的P波初始速度模型
Table 1 Initial P wave velocity model used in this study
深度/km vP/（km·s−1） 深度/km vP/（km·s−1） 0 4.4 20 6.0 5 5.6 25 6.1 10 5.8 30 6.3 15 5.9 50 8.0 
下载: 导出CSV
-
常祖峰,张艳凤,周青云,虎雄林,臧阳. 2014. 2013年洱源MS5.5地震烈度分布及震区活动构造背景研究[J]. 中国地震,30(4):560–570. doi: 10.3969/j.issn.1001-4683.2014.04.009 Chang Z F,Zhang Y F,Zhou Q Y,Hu X L,Zang Y. 2014. Intensity distribution characteristics and active tectonic background in area of the 2013 Eryuan MS5.5 earthquake[J]. Earthquake Research in China,30(4):560–570 (in Chinese).
常祖峰,常昊,臧阳,代博洋. 2016. 维西—乔后断裂新活动特征及其与红河断裂的关系[J]. 地质力学学报,22(3):517–530. doi: 10.3969/j.issn.1006-6616.2016.03.009 Chang Z F,Chang H,Zang Y,Dai B Y. 2016. Recent active features of Weixi-Qiaohou fault and its relationship with the Honghe fault[J]. Journal of Geomechanics,22(3):517–530 (in Chinese).
胡亚平,王志,刘冠男,柳存喜,伏毅. 2017. 南北地震带地壳结构多参数成像及强震触发机制研究[J]. 地球物理学报,60(6):2113–2129. doi: 10.6038/cjg20170608 Hu Y P,Wang Z,Liu G N,Liu C X,Fu Y. 2017. Crustal structure imaging of multi-geophysical parameters and generating mechanisms of large earthquakes in North-South Seismic Zone[J]. Chinese Journal of Geophysics,60(6):2113–2129 (in Chinese).
李大虎,詹艳,丁志峰,高家乙,吴萍萍,孟令媛,孙翔宇,张旭. 2021. 四川长宁MS6.0地震震区上地壳速度结构特征与孕震环境[J]. 地球物理学报,64(1):18–35. doi: 10.6038/cjg2021O0241 Li D H,Zhan Y,Ding Z F,Gao J Y,Wu P P,Meng L Y,Sun X Y,Zhang X. 2021. Upper crustal velocity and seismogenic environment of the Changning MS6.0 earthquake region in Sichuan,China[J]. Chinese Journal of Geophysics,64(1):18–35 (in Chinese).
李姣,姜金钟,杨晶琼. 2020. 2017年漾濞MS4.8和MS5.1地震序列的微震检测及重定位[J]. 地震学报,42(5):527–542. doi: 10.11939/jass.20190161 Li J,Jiang J Z,Yang J Q. 2020. Microseismic detection and relocation of the 2017 MS4.8 and MS5.1 Yangbi earthquake sequence,Yunnan[J]. Acta Seismologica Sinica,42(5):527–542 (in Chinese).
李敏娟,沈旭章,张元生,刘旭宙,梅秀苹. 2018. 基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析[J]. 地球物理学报,61(5):2075–2087. doi: 10.6038/cjg2018L0720 Li M J,Shen X Z,Zhang Y S,Liu X Z,Mei X P. 2018. Fine crustal structures of northeast margin of the Tibetan Plateau and structural features of Jiuzhaigou earthquake focal area constrained by the data from a high-density seismic array[J]. Chinese Journal of Geophysics,61(5):2075–2087 (in Chinese).
李永华,吴庆举,安张辉,田小波,曾融生,张瑞青,李红光. 2006. 青藏高原东北缘地壳S波速度结构与泊松比及其意义[J]. 地球物理学报,49(5):1359–1368. doi: 10.3321/j.issn:0001-5733.2006.05.015 Li Y H,Wu Q J,An Z H,Tian X B,Zeng R S,Zhang R Q,Li H G. 2006. The Poisson ratio and crustal structure across the NE Tibetan Plateau determined from receiver functions[J]. Chinese Journal of Geophysics,49(5):1359–1368 (in Chinese).
刘伟,吴庆举,张风雪. 2019. 利用双差层析成像方法反演青藏高原东南缘地壳速度结构[J]. 地震学报,41(2):155–168. doi: 10.11939/jass.20180083 Liu W,Wu Q J,Zhang F X. 2019. Crustal structure of southeastern Tibetan Plateau inferred from double-difference tomography[J]. Acta Seismologica Sinica,41(2):155–168 (in Chinese).
马宏生,张国民,闻学泽,周龙泉,邵志刚. 2008. 川滇地区三维P波速度结构反演与构造分析[J]. 地球科学:中国地质大学学报,33(5):591–602. Ma H S,Zhang G M,Wen X Z,Zhou L Q,Shao Z G. 2008. 3-D P wave velocity structure tomographic inversion and its tectonic interpretation in Southwest China[J]. Earth Science:Journal of China University of Geosciences,33(5):591–602 (in Chinese).
缪思钰,张海江,陈余宽,谭玉阳,苗园园,黄振华,王飞,谢庆明. 2019. 基于微地震定位和速度成像的页岩气水力压裂地面微地震监测[J]. 石油物探,58(2):262–271. doi: 10.3969/j.issn.1000-1441.2019.02.012 Miao S Y,Zhang H J,Chen Y K,Tan Y Y,Miao Y Y,Huang Z H,Wang F,Xie Q M. 2019. Surface microseismic monitoring of shale gas hydraulic fracturing based on microseismic location and tomography[J]. Geophysical Prospecting for Petroleum,58(2):262–271 (in Chinese).
孙权,裴顺平,苏金蓉,刘雁冰,薛晓添,李佳蔚,李磊,左洪. 2021. 2019年6月17日四川长宁MS6.0地震震源区三维速度结构[J]. 地球物理学报,64(1):36–53. doi: 10.6038/cjg2021O0246 Sun Q,Pei S P,Su J R,Liu Y B,Xue X T,Li J W,Li L,Zuo H. 2021. Three-dimensional seismic velocity structure across the 17 June 2019 Changning MS6.0 earthquake,Sichuan,China[J]. Chinese Journal of Geophysics,64(1):36–53 (in Chinese).
王椿镛,Mooney W D,王溪莉,吴建平,楼海,王飞. 2002. 川滇地区地壳上地幔三维速度结构研究[J]. 地震学报,24(1):1–16. Wang C Y,Mooney W D,Wang X L,Wu J P,Lou H,Wang F. 2002. Study on 3-D velocity structure of crust and upper mantle in Sichuan-Yunnan region,China[J]. Acta Seismologica Sinica,24(1):1–16 (in Chinese).
王伟平. 2016. 鲜水河、安宁河和龙门山断裂带交汇区双差层析成像研究[D]. 北京: 中国地震局地球物理研究所: 39−70. Wang W P. 2016. The Application of Double-Difference Tomography Method in the Intersection Zone of Xianshuihe Fault, Anninghe Fault and Longmenshan Fault[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 39−70 （in Chinese）.
韦伟,孙若昧,石耀霖. 2010. 青藏高原东南缘地震层析成像及汶川地震成因探讨[J]. 中国科学:地球科学,40(7):831–839. Wei W,Shi R M,Shi Y L. 2010. P-wave tomographic images beneath southeastern Tibet:Investigating the mechanism of the 2008 Wenchuan earthquake[J]. Science China Earth Science,53(9):1252–1259.
吴建平,明跃红,王椿镛. 2006. 川滇地区速度结构的区域地震波形反演研究[J]. 地球物理学报,49(5):1369–1376. doi: 10.3321/j.issn:0001-5733.2006.05.016 Wu J P,Ming Y H,Wang C Y. 2006. Regional waveform inversion for crustal and upper mantle velocity structure below Chuan-Dian region[J]. Chinese Journal of Geophysics,49(5):1369–1376 (in Chinese).
肖卓,高原. 2017. 利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构[J]. 地球物理学报,60(6):2213–2225. doi: 10.6038/cjg20170615 Xiao Z,Gao Y. 2017. Crustal velocity structure beneath the northeastern Tibetan Plateau and adjacent regions derived from double difference tomography[J]. Chinese Journal of Geophysics,60(6):2213–2225 (in Chinese).
徐锡伟, 韩竹军, 杨晓平, 张世民, 于贵华, 周本刚, 李峰, 马保起, 陈桂华, 冉永康. 2016. 中国及邻区地震构造图[M]. 北京: 地震出版社: 1. Xu X W, Han Z J, Yang X P, Zhang S M, Yu G H, Zhou B G, Li F, Ma B Q, Chen G H, Ran Y K. 2016. Seismotectonic Map in China and Its Adjacent Regions[M]. Beijing: Seismological Press: 1 （in Chinese）.
胥颐,杨晓涛,刘建华. 2013. 云南地区地壳速度结构的层析成像研究[J]. 地球物理学报,56(6):1904–1914. doi: 10.6038/cjg20130613 Xu Y,Yang X T,Liu J H. 2013. Tomographic study of crustal velocity structures in the Yunnan region southwest China[J]. Chinese Journal of Geophysics,56(6):1904–1914 (in Chinese).
云南省地震局. 2021. 云南漾濞6.4级地震烈度图发布[EB/OL]. [2021-05-25]. http://yndzj.gov.cn/yndzj/_300559/_300651/629959/index.html. Yunnan Earthquake Agency. 2021. Seismic intensity map of the M6.4 Yangbi earthquake[EB/OL]. [2021-05-25]. http://yndzj.gov.cn/yndzj/_300559/_300651/629959/index.html （in Chinese）.
Eberhart-Phillips D. 1986. Three-dimensional velocity structure in northern California Coast Ranges from inversion of local earthquake arrival times[J]. Bull Seismol Soc Am,76(4):1025–1052.
Lei J S,Zhao D P,Su Y J. 2009. Insight into the origin of the Tengchong intraplate volcano and seismotectonics in Southwest China from local and teleseismic data[J]. J Geophys Res:Solid Earth,114(B5):B05302. doi: 10.1029/2008JB005881
Lei J S,Zhao D P. 2016. Teleseismic P-wave tomography and mantle dynamics beneath eastern Tibet[J]. Geochem Geophys Geosyst,17(5):1861–1884. doi: 10.1002/2016gc006262
Lei J S,Zhao D P,Xu X W,Xu Y G,Du M F. 2019. Is there a big mantle wedge under eastern Tibet?[J]. Phys Earth Planet Inter,292:100–113.
Liu Y, Yao H J, Zhang H J, Fang H J. 2021. The community velocity model V.1.0 of Southwest China, constructed from joint body- and surface-wave travel-time tomography[J]. Seismol Res Lett: 1–16. doi: 10.1785/0220200318.
Nakajima J,Matsuzawa T,Hasegawa A,Zhao D P. 2001. Three-dimensional structure of vP,vS and vP/vS beneath northeastern Japan:Implications for arc magmatism and fluids[J]. J Geophys Res:Solid Earth,106(B10):21843–21857. doi: 10.1029/2000JB000008
Shen W S,Ritzwoller M H,Kang D,Kim Y H,Lin F C,Ning J Y,Wang W T,Zhang Y,Zhou L Q. 2016. A seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion[J]. Geophys J Int,206(2):954–979. doi: 10.1093/gji/ggw175
Waldhauser F,Ellsworth W. 2000. A double-difference earthquake location algorithm:Method and application to the northern Hayward fault,California[J]. Bull Seismol Soc Am,90(6):1353–1368. doi: 10.1785/0120000006
Xiao X,Cheng S H,Wu J P,Wang W L,Sun L,Wang X X,Wen L X. 2021. Shallow seismic structure beneath the continental China revealed by P-wave polarization,Rayleigh wave ellipticity and receiver function[J]. Geophys J Int,225(2):998–1019. doi: 10.1093/gji/ggab022
Zhang H J,Thurber C H. 2003. Double-difference tomography:The method and its application to the Hayward fault,California[J]. Bull Seismol Soc Am,93(5):1875–1889. doi: 10.1785/0120020190
Zhang H J,Thurber C H. 2006. Development and applications of double-difference seismic tomography[J]. Pure Appl Geophys,163(2/3):373–403. doi: 10.1007/s00024-005-0021-y
Zheng C,Zhang R Q,Wu Q J. 2019. Variations in crustal and uppermost mantle structures across eastern Tibet and adjacent regions:Implications of crustal flow and asthenospheric upwelling combined for expansions of the Tibetan Plateau[J]. Tectonics,38(8):3167–3181. doi: 10.1029/2018TC005276
-
期刊类型引用(8)
1. 谢张迪，于湘伟，章文波. 2022年青海门源M_W6.6地震震源动力学破裂过程. 地震学报. 2025(01): 21-36 . 
本站查看
2. Qiong Wang，ShuYu Li，XinYi Li，Yue Wu，PanPan Zhao，Yuan Gao. Deep tectonics and seismogenic mechanisms of the seismic source zone of the Jishishan M_s 6.2 earthquake on December 18, 2023, at the northeast margin of the Tibetan Plateau. Earth and Planetary Physics. 2024(03): 514-521 . 必应学术
3. 高原，李心怡，李抒予，夏新宇，杨逸文，王琼. 2023年12月18日积石山6.2级地震的深浅变形构造分析. 地震. 2024(01): 160-166 . 百度学术
4. 刘会龙，陈楚桐，范俊杰，马德锡，彭伍胥，陈亮. 青藏高原东北缘祁连至阿拉善地块深部电性结构特征及其地质意义. 地球物理学进展. 2024(03): 916-929 . 百度学术
5. 于书媛，黄显良，郑海刚，李玲利，骆佳骥，丁娟，范晓冉. 2022年门源M_W6.7地震的同震破裂模型及应力研究. 地震地质. 2023(01): 286-303 . 百度学术
6. 万悦，苏鹤军，李晨桦，周慧玲. 祁连山断裂带中东段断层土壤气氡浓度强度时空分布特征及其地震危险性分析. 地震工程学报. 2023(02): 491-500 . 百度学术
7. 李莹，田建慧，李心怡，李抒予，王琼，高原. 2022年9月5日四川泸定M_S6.8地震深部构造特征. 地球物理学报. 2023(04): 1385-1396 . 百度学术
8. 赵云峰，祝意青，隗寿春，刘芳，梁伟锋，孙和平. 2022年1月8日青海门源M_S6.9地震前重力场动态变化. 地球物理学报. 2023(06): 2337-2351 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 1074
- HTML全文浏览量: 354
- PDF下载量: 267
- 被引次数: 10
