Joint iterative inversion of tectonic stress field and fault identification in Yingjiang area
-
摘要: 搜集了盈江及其邻区5组地震序列震源机制解并反演得到该区域的构造应力场。研究结果显示：盈江地区整体主压应力以NNE向为主，主张应力以ESE向为主；但其局部应力场不完全一致，沿苏典断裂分布的主压应力轴走向随着断裂走向由北向南延伸角度逐渐向北偏移，而盈江地区的西南部，其主压应力走向更偏向于东，这可能与大盈江断裂的横向拉伸有关。此外，本研究通过应力场反演识别出了5组地震序列震源机制解的主发震断层节面的走向、倾角、滑动角及发震断层的摩擦系数，为今后该区域的地震研究及地壳动力学变迁提供了参考。Abstract: There has always been a problem in inversion of tectonic stress field by focal mechanism solution, it is impossible to determine which nodel plane is the correct seismogenic fault plane. The joint iterative stress inversion method can effectively avoid the error caused by incorrect fault selection by identifying the fault instability. In this paper, the focal mechanism solutions of five groups of earthquake sequences in Yingjiang and its adjacent areas are collected, and the tectonic stress field in this area is inversed. The results show that, in Yingjiang area, the main compressive stress is NNE, and the main tensile stress is ESE, which is basically consistent with the previous research results. However, the local stress field is not completely consistent. The strike of the principal stress axis along the Sudian fault extends from north to south, and the angle gradually shifts to the north. In the southwest of Yingjiang, the strike of the principal stress tends to the East, which may be related to the transverse extension of the Dayingjiang fault. In addition, the strike, dip angle, slip angle and friction coefficient of the main seismogenic fault nodal planes of the focal mechanism solutions of five groups of earthquake sequences are identified by stress field inversion. This provides a valuable reference for the future seismic and crustal dynamic changes research in this region.
-
引言
由于水库蓄水或排水而引发的库区附近地震活动性（地震的频次和地震的震级）的明显增强，称为水库诱发地震，或简称水库地震（陈颙，2009）。自二十世纪六十年代全球接连发生M6以上水库诱发强震后，水库诱发地震在全世界范围内引起了广泛关注。广西龙滩水库横跨广西天峨县和贵州罗甸县，地处云贵高原南缘向桂西山区与丘陵过渡的斜坡地带，坝高220 m，正常蓄水位为400 m，总库容为273亿m3，是典型的高山峡谷高坝大库容水电站（周斌等，2014；阎春恒等，2015）。该水库自2006年9月蓄水以来，库区地震活动明显增多，截止到2019年底，已经记录到近5 000余次地震（图1），其中ML＞4.0地震3次，对库区周边人们的生产生活造成较大影响。库区断裂发育（图1），九条主要断层中除罗甸—望谟断裂（F1）、达恒—达良断裂（F5）、望谟—逻西断裂（F7）外，其余均在早—中更新世有过不同程度的活动（周斌等，2014；阎春恒等，2015）。在库区的出露地层中（图1），泥盆系—二叠系地层绝大部分为碳酸盐岩建造，岩溶、裂隙发育程度较高，渗透性较强（阎春恒等，2015；叶庆东等，2018）。作为国内最重要的水利枢纽工程之一，龙滩水库在带来巨大经济价值的同时，也持续吸引着学者们对库区地震地质构造条件（潘建雄，1989）、危险性评估（郭培兰等，2006）、活动性与孕震环境 （陈翰林等，2009a，b；詹艳等，2012；阎春恒等，2015，2020；叶庆东等，2018）、诱发的物理机制（刘耀炜等，2011；周斌等，2014）等方面的关注并开展相关研究，取得了许多重要的认识。
![]()
图 1 龙滩水库地质构造及2006年9月至2019年10月的地震事件分布F1：罗甸—望谟断裂；F2：高圩—八茂断裂；F3：凤亭—下老断裂；F4：马耳—拉浪断裂；F5：达恒—达良断裂；F6：党明—桂花断裂；F7：望谟—逻西断裂；F8：长里—八南断裂；F9：龙凤—八腊断裂Figure 1. Structural outlines and the distribution of seismic events from September of 2006 to October of 2019 in the Longtan reservoirF1：Luodian-Wangmo fault；F2：Gaoyu-Bamao fault；F3：Fengting-Xialao fault；F4：Maer-Lalang fault；F5：Daheng-Daliang fault；F6：Dangming-Guihua fault；F7：Wangmo-Luoxi fault；F8：Changli-Ba’nan fault；F9：Longfeng-Bala fault水库地震主要受库水载荷与流体作用两方面的影响（蒋海昆等，2014）。水库蓄水或排水受降水量及生产需要的影响一般呈年度和季节性变化，载荷与流体作用响应也呈现相同的变化规律，进而导致地震活动性也呈现相同的变化趋势（Nur，Booker，1972）。基于此，近年来一些学者（Saar，Manga，2003；Christiansen et al，2005；Telesca，2010；Telesca et al，2015，2020）直接对地震活动与水位进行相关分析，获取水库地震活动相对于水位的延迟时间，并通过延迟时间估算库区的孔隙压力扩散系数。该系数是一个非常重要的物理参数，其数值大小直接关系到地震活动的机理，并影响库区地震活动的时空尺度，一直是地震和水文地质工作者关注的重点。然而，这方面的研究在国内并不多见，鉴于此，本文拟对龙滩水库2013年3月1日至2019年10月31日的地震活动与水位变化进行相关分析，计算地震活动对库区蓄水响应的延迟时间，进而根据延迟时间估算龙滩水库库区的孔隙压力扩散系数。
1. 数据
龙滩水库从2006年9月开始蓄水，水位在3个月左右上升100余m，几乎与此同时就触发了大量的地震。这些地震属于对库区蓄水作用的即时响应，受到库区自身的应力状态、蓄水量短时间内大量增加所致的加载和孔隙压力扩散等因素的影响。一般认为，库水加载对诱发地震的作用非常有限甚至是抑制地震发生的（Roeloffs，1988；蒋海昆等，2014）；蓄水量短时间内大量增加引起孔隙压力增加的影响时段通过分析误差余补函数的变化来获得，通常为几个月到几年（Talwani，Acree，1984；Roeloffs，1988；Turcotte，Schubert，2002；蒋海昆等，2014）。本文仅考虑蓄水完成后水位波动与地震活动的关系。为了削弱即时响应的影响，我们收集了龙滩水库完成蓄水6年以后，2013年3月1日至2019年10月31日共六年半的地震目录。该目录共包含516个地震事件，其N-t和M-t分布见图2，这些地震均已由广西地震台网进行了精确定位。采用最大曲率法（Wiemer，Wyss，2000）计算得到这个时段内的最小完整性震级为0.5，为此我们仅保留震级不小于0.5的事件信息。图2b还给出了水位的日均值变化，可以看出水位基本呈周年变化，并且对地震活动具有一定的调控作用，在高水位时期地震活动强度较大，反之在低水位时期地震活动相对较弱。
![]()
图 2 龙滩水库2013年3月1日至2019年10月31日地震活动和水位变化（a） 地震活动N-t图；（b） M-t图和水位变化Figure 2. The seismicity and the variation of water level from March 1，2013 to October 31，2019 in the Longtan reservoir（a） The N-t map of seismicity；（b） The M-t map of seismicity and the change of water level2. 地震活动与水位变化的互相关分析
Nur和Booker （1972）在研究流体孔隙压力触发余震活动时提出，余震活动水平与孔隙压力变化成正比。对水库诱发地震而言，孔隙压力的变化表现为水库水位的变化，地震活动水平表现为一定时段内的地震数目或这些地震释放的地震矩（Saar，Manga，2003；Christiansen et al，2005），本文选择前者作为地震活动水平的表征。我们先将筛选出的地震目录整理成以天为单位的地震次数时间序列N，然后对每天的地震次数N与水位W进行标准化处理，即减去均值后再除以标准差，分别得到新的时间序列n和w （Telesca et al，2012）。类似于Saar和Manga （2003）、Telesca （2010）、Telesca等（2015）的做法，对序列n与w作互相关后，再采用替代数据来检验。假设地震发生与水位不存在关系，则n与w的互相关系数没有明显的极大值，或者极大值仅仅是因为原始数据包含噪声而产生的假象。为了检验这种假设，通过替代数据生成方法产生一系列与w具有相同功率谱但相位谱随机的时间序列wR来“重新洗牌”，尽可能破坏掉原序列的原始动力结构（Little et al，2006）。如果地震活动与水位不相关，则n与w和n与wR的互相关系数统计上处于同一水平，在n与w最大互相关系数对应的延迟时间附近，其数值相对于n与wR的互相关系数没有明显优势。反之，如果地震活动与水位相关，则在n与w最大互相关系数对应的延迟时间附近，其值显著地大于n与wR的互相关系数。
本文中采用傅里叶变换法（Saar，Manga，2003；刘鑫，鲍长春，2014）来生成替代数据，总共产生了1 000组替代时间序列，每一组w的替代数据wR都与n作互相关，得到1 000组互相关系数。由于互相关系数的每一个时间点对应一个延迟时间，1 000组数据意味着每一个延迟时间处都有1 000个互相关系数值与之对应，在这些值中取出不小于序列中99%值的最小值（1 000个数据的序列，从小到大排列取第990个点）作为该点的置信水平。对每一个延迟时间对应的1 000个互相关系数作同样处理，得到99%置信水平的曲线。图3a中灰色曲线给出了n与w的互相关系数及基于替代数据得到的99%的置信水平，可以看出最大互相关系数出现在延迟时间37天处，最大值为0.186，尽管相关性不高，但却高于99%的置信水平。
![]()
图 3 标准化地震数目与水位的互相关（a） 以天为单位的标准化后的地震数目与水位的互相关，其中灰色代表滤波之前的结果，黑色为滤波后的结果；（b） 以月为单位的标准化后的地震数目与水位的互相关Figure 3. Cross-correlation between the standardized number of earthquakes and water level（a） Cross-correlation between standardized number of earthquakes and water level in day，where the gray line represents the results before filtering，and the black line represents the results after filtering；（b） Cross-correlation between standardized number of earthquakes and water level in month与水位相比，地震活动显然随机性更大，频谱成分的差异可能是导致n与w互相关系数小的直接原因。为此，我们对时间序列n和w进行了功率谱分析，其结果见图4。可以看出：w含有比较明显的年度、半年度和季节性变化；而n除了年度变化比较明显之外，半年度及季节性变化不明显，高频成分较多。考虑到库区地震活动可能主要受水位长周期变化的影响，我们对序列w和n进行低通滤波，对于w保留周期大于一个季度的成分，而对于n仅保留周期大于150天的成分。w和n滤波后的时间序列ws和ns见图5a，在箭头1—3处，可以看出ws与ns相关性明显，地震活动的变化ns相对于水位ws略为滞后。此外，我们以月为单位统计地震数目（nm）和水位数据（wm），其中水位数据取每个月的平均值，其结果见图5b。相对于以日为单位的统计，以月为单位的统计也有低通滤波的作用。由于时间间隔的拉大和序列长度的缩短，在图5b中地震数目相对于水位变化的滞后更加明显，其中在箭头1，2，4和5处的地震活动相对于水位滞后1个月，在箭头6处滞后3个月。
![]()
图 4 标准化后的水位w （a）和地震数目n （b）的功率谱Figure 4. The power spectra of standardized water level w （a） and number n of earthquakes （b）![]()
图 5 （a） 滤波后的标准化地震数目ns和水位ws；（b） 以月为单位的标准化地震数目nm和水位wmFigure 5. （a） Standardized number of earthquakes ns and water level ws after filtering；（b） Standardized number of earthquakes nm and water level wm in month同处理n和w一样，我们对ns与ws和nm与wm也进行互相关计算，并用替代数据进行检验。图3a中黑色曲线为ns与ws的互相关系数及其99%置信水平曲线，对比滤波前的灰色曲线，可见互相关系数形态基本相同，但是滤波后的最大相关系数达到了0.67，是滤波前的3.6倍，表明水位的变化对地震活动产生的调控确实主要来自于其长周期波动。滤波前后两者的延迟时间均为37天，说明高频变化仅改变了相关系数的数值大小，掩盖不了水位长周期变化是引起库区地震活动的主要因素这一实质。图3b为nm与wm的互相关系数及其99%置信水平曲线，最大相关系数为0.58，对应的最大延迟时间为1个月，与ns和ws互相关系数的结果一致，进一步确认了延迟时间为37天是可靠的。
3. 扩散系数估算
孔隙压力的增加降低了岩石的有效正应力和岩石摩擦系数，从而使介质的抗剪强度降低。当岩石处于破裂的临界状态时，孔隙压力微小的增加都可能触发地震，因此孔隙压力扩散被认为是触发水库地震的一个主要机制（Nur，Booker，1972；Talwani，Acree，1984；Roeloffs，1996）。孔隙压力扩散由扩散方程描述，在扩散方程中，扩散系数是衡量扩散快慢的物理参数，与介质的渗透系数、流体的动力学黏度及孔隙弹性模量等参数的大小相关。地表附近的孔隙压力扩散系数可以通过相关试验获取，深部的孔隙压力扩散系数获取非常困难。水库库区的水位变化与地震活动的关系提供了一个探测库区孔隙压力扩散系数的“窗口”，是孔隙压力扩散系数原位测试的试验场（Talwani，Acree，1984；龚钢延，谢原定，1991）。
3.1 基于自由扩散模式的扩散系数估算
在基于扩散方程估算孔隙压力扩散系数的已有研究中，大多忽略孔隙压力扩散与应力的耦合（Talwani，Acree，1984；Saar，Manga，2003；Telesca，2010；Telesca et al，2015），仅考虑孔隙压力的自由扩散。在此情形下仅考虑垂向扩散，一维孔隙压力扩散满足的方程为：
$$ D\frac{{{\partial ^2}P}}{{\partial {h^2}}} ＝ \frac{{\partial P}}{{\partial t}} \text{，} $$ (1) 式中，P为孔隙压力，D为孔隙压力扩散系数，h为深度，t为时间。在水库坝基处孔隙压力主要受库水周年性变化调控，则边界条件为P（t，h＝0）＝P0exp（iωt），其中P0为坝基处（h=0）由水位变化引起的压力周年性变化的振幅，ω＝2π/T，T为周期，T＝1 a，i为虚数单位。式（1）的解为（Turcotte，Schubert，2002）：
$$ \begin{split} \\[-6pt] \frac{P}{{{P_0}}} ＝ \exp \left( { - h\sqrt {\frac{\omega }{{2D}}} } \right)\exp \left[ {{\rm{i}}\omega \left( {t － \frac{h}{{\sqrt {2\omega D} }}} \right)} \right] {\text{.}} \end{split}$$ (2) 式（2）表明扩散方程在周期性的边界条件下，其解可表示为一个振幅衰减的平面波，衰减系数为$ \sqrt {{\omega / {2D}}}$，波速为$ \sqrt {2\omega D}$。指数项$\exp （ - h\sqrt {{\omega \mathord{\left/ {\vphantom {\omega {2D}}} \right. } {2D}}} ） $描述了孔隙压力振幅随深度h的衰减，而谐波项$\exp ［ {\rm{i}}\omega （ t －{h / {\sqrt {2\omega D} }} ） ］ $描述了孔隙压力作为时间t与深度h的函数周期性变化的规律，两者均与扩散系数D相关。在某一深度岩石已处于破裂的临界应力状态，孔隙压力扰动前锋传播至此触发地震，则应满足$t －h / {\sqrt {2\omega D} ＝ 0}$，整理得
$$ D ＝ \frac{{{h^2}}}{{2\omega {t^2}}} \text{，} $$ (3) 式中，t＝37 d为延迟时间，h取趋肤深度（Saar，Manga，2003）或者地震分布的优势深度（Telesca，2010；Telesca et al，2015），本文取后者。图6给出了库区地震的深度分布图，可见深度主要集中在5—7 km，均值μ为5.94 km，标准偏差σ为1.44 km。将ω，t和h的值代入式（3），可得D＝（8.66±4.11） m2/s，其中，不确定度由地震深度分布的一倍标准偏差确定。
在自由扩散模式下，基于扩散方程在周期性边界条件下的解和库区地震平均深度得到的龙滩水库孔隙压力扩散系数D＝（8.66±4.11） m2/s，处于0.1—10 m2/s的合理范围内（Talwani，Acree，1984）。然而，该结果与刘耀炜等（2011）基于水压致裂中使用的单点触发模型采用抛物线包络拟合方法所得到的龙滩库区蓄水初期的孔隙压力扩散系数（0.04—0.65 m2/s）差异较大。我们认为差异可能源于以下三方面。其一，根据地质资料（广西壮族自治区地质局，1972；潘建雄，1989），龙滩库区以泥盆系—二叠系碳酸盐岩建造为主，岩溶、裂隙发育程度较高，渗透性较好，受地震定位精度的影响，加之抛物线包络拟合具有一定的主观性，刘耀炜等（2011）估算的值可能偏低。其二，刘耀炜等（2011）是在水位快速上升阶段，假定孔隙压力边界条件为阶跃函数得到的结果；而本文结果是在蓄水完成六年以后，基于水位总体稳定仅有周年扰动变化得到的结果，不同时段库区的介质状态不一样，因而扩散系数也不一样。其三，自由扩散忽略孔隙压力与应力的耦合，模式过于简单。
3.2 基于耦合模式的扩散系数估算
虽然自由扩散模式下得到的孔隙压力扩散系数在合理范围内，但物理上孔隙压力与应力是耦合的，自由扩散模式仍过于简化。根据Biot （1941）的固结理论，Roeloffs （1988）得到了周期性边界条件下一维孔隙压力与应力耦合模式扩散方程：
$$ \frac{1}{\alpha }\frac{{{\partial ^2}}}{{\partial {h^2}}}\frac{P}{{{P_0}}} ＋ \frac{{{\rm{i}}\omega }}{D} ＝ \frac{{{\rm{i}}\omega P}}{{\alpha D{P_0}}} \text{，} $$ (4) 式中，$\alpha ＝ {{B （ 1 ＋{\nu _{\rm{u}}} ） } / { ［ {3 （ 1 － {\nu _{\rm{u}}} ） } ］ }}$，B为Skempton系数，νu为无排水条件下孔隙介质的泊松比。式（4）的解为
$$ \frac{P}{{{P_0}}} ＝ （ {1 － \alpha } ） \exp \left[ { - （ {1 ＋ {\rm{i}}} ） \sqrt {{{{\frac{\omega }{{2D}}}} }}h} \right] ＋ \alpha {\text{.}} $$ (5) 式（5）表明深度h处的孔隙压力除了与扩散率D和角频率ω有关，还与系数α有关，α为解耦模式与无排水响应模式的加权平均。当α＝0时，式（5）与式（2）等价，对应于孔隙压力与应力解耦的模式；当α＝1时，对应于无排水响应模式。
式（5）中，对于给定$ \sqrt{\omega /2D}h$，相位仅与α相关，可以通过相位的变化分析时间的超前与滞后，相位正代表时间的超前，相位负对应着时间的滞后，最小值对应着最大延迟时间。图7给出了α＝0，0.05，0.30，0.70，1.00时，孔隙压力及相位与$ \sqrt{\omega /2D}h$ 的关系，可见在相同的$ \sqrt{\omega /2D}h$ 情况下，α越大，孔隙压力越大，延迟时间越小。已知延迟时间为37 d，通过搜索相位的最小值求得对应的α＝0.219，对应的$ \sqrt{\omega /2D}h$＝1.43，ω=2π/T，周期T取1 a，h取地震的平均深度，可得扩散系数D＝（1.72±0.82） m2/s。可见，是否考虑孔隙压力扩散与应力耦合，估算出来的扩散系数差异较大。这种差异可以由延迟时间直观地判断，当延迟时间较大（接近180 d）时，是否考虑孔隙压力扩散与应力的耦合对估算D影响不大；当延迟时间较小（远小于180 d）时，是否考虑孔隙压力与应力的耦合则可能对估算D有一定影响。由$ \sqrt{\omega /2D}h$＝1.43得到对应的$ P/{P}_{0}＝0.307 $，根据图1，龙滩水库水位年变化在330—370 m之间，水位变化Δz＝40 m，则${P}_{0}＝\rho g\Delta {\textit{z}}/2$，其中ρ为水的密度，取1 000 kg/m3，g为重力加速度，取9.8 N/kg，则P＝60 kPa，孔隙压力的年变化为120 kPa。假定岩石强度的降低完全由孔隙压力的变化引起，120 kPa的孔隙压力增加会导致近100 kPa的库仑应力变化，在临界状态下足以触发地震（Hardebeck et al，1998；Stein，1999）。
![]()
图 7 α取不同值时孔隙压力幅值（a）及相位（b）随$ \sqrt{\omega /2D}h $ 变化Figure 7. Variation of the magnitude （a） and phase （b） of pore pressure as a function of $ \sqrt{\omega /2D}h $ for different α4. 讨论与结论
通过对龙滩水库的水位与地震数目的互相关分析，得到了地震活动相对于水位的延迟时间为37天。与一些蓄水多年的水库（Talwani，Acree，1984；Ferreira et al，1995；Telesca，2010）相比，37天的延迟时间是一个较小的值，由于孔隙压力扩散范围正比于延迟时间（Talwani，Acree，1984；龚钢延，谢原定，1991），这暗示库水的影响主要局限于库区较小的范围内，这与库区地震精定位的结果（阎春恒等，2015，2020）一致。由于水库地震震源深度较浅，通常M3以上地震均会有震感，引发社会舆论，对人们的生活生产造成影响，地震活动强度相对于水位峰值37天的延迟时间为相关部门制定震情跟踪、舆情应对等预案提供了时间节点信息。
在不考虑孔隙压力扩散与应力耦合的情况下，基于周期性边界条件下孔隙压力扩散方程的解和库区地震的深度分布，利用37天的延迟时间得到扩散系数D＝（8.66±4.11） m2/s；在考虑孔隙压力扩散与应力耦合的情况下，得到扩散系数D＝（1.72±0.82） m2/s，两者均处于合理的范围内。数值上，基于孔隙压力与应力耦合模式得到的扩散系数与刘耀炜等（2011）的结果比较接近，也与近年来部分学者（Saar，Manga，2003；Telesca，2010；Telesca et al，2015）利用地震活动与水位变化的互相关分析得到的其它库区的扩散系数接近，但由于这些研究大多采用了自由扩散的假定，因此无法说明D＝（1.72±0.82） m2/s比D＝（8.66±4.11） m2/s更为合理。从物理上以及近年来对断裂带的孔隙压力扩散系数原位测量的结果（Doan et al，2006；Xue et al，2016）来看支持孔隙压力扩散系数取较小的值，这可能说明对龙滩水库而言，孔隙压力扩散与应力耦合是现阶段诱发地震的主要因素。
-
![]()
图 1 盈江区域地震序列震源机制解空间分布图
F1：怒江断裂；F2：龙陵—瑞丽断裂；F3：瓦德龙断裂；F4：大盈江断裂；F5：昔马—盘龙山断裂；F6：苏典断裂；F7：槟榔断裂；F8：猴桥—中和断裂；F9：固东—腾冲断裂图中A−E区分别代表2008年3月21日、8月21日及9月3日盈江地震序列、2011年3月10日盈江地震序列和2011年8月9日腾冲地震序列的震源机制解分区
Figure 1. Spatial distribution of focal mechanism solutions of earthquake sequences in Yingjiang region
F1：Nujiang fault；F2：Longling-Ruili fault；F3：Wadelon fault；F4：Dayingjiang fault；F5：Xima-Panlongshan fault；F6：Sudian fault；F7：Binglang fault；F8：Houqiao-Zhonghe fault；F9：Gudong-Tengchong fault.A−E in the figure represent the focal mechanism solution partition of Yingjiang earthquake sequences on 21 March 2008，21 August 2008，3 September 2008，10 March 2011 and Tengchong earthquake sequence on 9August 2011，respectively
![]()
图 2 断层不稳定性的莫尔圆定义（Vavryčuk，2014）
红色圆点标出了以不稳定性为特征I＝1的主断层上的牵引力，蓝色圆点标出了以任意不稳定性为特征I的主断层上的牵引力
Figure 2. Definition of the fault instability in Mohr’s diagram （Vavryčuk，2014）
The red solid circle marks the traction on the main fault characterized by instability I＝1，and the blue solid circle marks the traction on the main fault characterized by arbitrary instability
![]()
图 3 A—E区（a—e）5组地震序列震源机制解联合迭代反演的应力场结果
左侧为地震序列震源机制解反演的断层摩尔圆分布，右侧为主轴方向及P，T轴
Figure 3. Stress field results of joint iterative inversion for focal mechanism solutions of five seismic sequences in areas A−E （a−e）
The left side is the fault Mohr’s circle distribution inversed by focal mechanism solutions of seismi sequences， and the right side is the main axis direction and P，T axis
![]()
图 4 5组地震序列构造应力场反演的应力形因子R的分布（左）及主应力轴方向（右）
Figure 4. Stress shape factor R distribution （left） and principal stress axis direction （right） of tectonic stress field inversion results for five earthquake sequences
（a） 2008-13-21；（b） 2008-08-21；（c） 2008-09-03；（d） 2011-03-10；（e） 2011-08-09
![]()
图 5 5组地震序列震源机制解应力反演识别出的主断层节面图
Ⅰ ：地震序列震源机制解的双节面下半球投影图；Ⅱ ：最优断层节面下半球投影；Ⅲ ：两组震源丛最优发震断层的主断层节面；Ⅳ ：主断层节面的P，T轴
Figure 5. Nodal map of main fault identified by stress inversion of focal mechanism solutions of five groups of earthquake sequences
Ⅰ ：Double-nodal lower hemispherical projection of focal mechanism solutions of the earthquake sequence；Ⅱ ：Lower hemispherical projection of optimal fault nodal plane；Ⅲ ：Optimal occurrence of the main fault nodal plane for the two sets of focal clusters；Ⅳ ：The P and T axis of the main fault nodal plane （a） 2008-03-21；（b） 2008-08-21；（c） 2008-09-03；（d） 2011-03-10；（e） 2011-08-09
![]()
图 6 主压应力轴应力场分布模型图
Figure 6. Stress field distribution model diagram of principal compressive stress axis
表 1 构造应力场应力主轴反演结果
Table 1 Inversion results of stress principal axis of tectonic stress field
地震序列主震
发震日期分区 震源机制解
个数σ1 σ2 σ3 应力形因子R 方位角/° 倾角/° 方位角/° 倾角/° 方位角/° 倾角/° 2008-03-21 A区 17 44 3 258 86 135 2 0.78 2008-08-21 B区 36 215 33 40 57 306 2 0.40 2008-09-03 C区 48 34 3 212 87 304 0 0.75 2011-03-10 D区 61 193 15 333 70 99 12 0.57 2011-08-09 E区 23 195 1 294 83 105 7 0.65 
下载: 导出CSV
表 2 迭代应力反演识别出的最优主断层节面
Table 2 Optimal main fault nodal plane identified by iterative stress inversion
地震序列主震
发震日期反演识别的两组主断层节面 地震序列主震
发震日期反演识别的两组主断层节面 走向/° 倾角/° 滑动角/° 走向/° 倾角/° 滑动角/° 2008-03-21 250 89 4 2011-03-10 158 88 161 199 86 −178 223 72 −7 2008-08-21 241 77 −31 2011-08-09 164 85 176 12 73 −151 225 84 3 2008-09-03 187 89 −177 60 88 −3
下载: 导出CSV
-
邓菲,刘杰. 2014. 2008年盈江地震序列的震源参数和震源机制相关系数研究[J]. 地震,34(2):22–34. doi: 10.3969/j.issn.1000-3274.2014.02.003 Deng F,Liu J. 2014. Source parameters and correlation coefficients of focal mechanisms for the 2008 Yingjiang earthquake sequence[J]. Earthquake,34(2):22–34 (in Chinese).
房立华,吴建平,张天中,黄静,王长在,杨婷. 2011. 2011年云南盈江MS5.8地震及其余震序列重定位[J]. 地震学报,33(2):262–267. doi: 10.3969/j.issn.0253-3782.2011.02.013 Fang L H,Wu J P,Zhang T Z,Huang J,Wang C Z,Yang T. 2011. Relocation of mainshock and aftershocks of the 2011 Yingjiang MS5.8 earthquake in Yunnan[J]. Acta Seismologica Sinica,33(2):262–267 (in Chinese).
郭祥云,陈学忠,王生文,王恒信. 2014. 川滇地区中小地震震源机制解及构造应力场的研究[J]. 地震工程学报,36(3):599–607. doi: 10.3969/j.issn.1000-0844.2014.03.0599 Guo X Y,Chen X Z,Wang S W,Wang H X. 2014. Focal mechanism of small and moderate earthquakes and tectonic stress field in Sichuan-Yunnan areas[J]. China Earthquake Engineering Journal,36(3):599–607 (in Chinese).
胡晓辉. 2020. 云南及邻区构造应力场研究[D]. 三河: 防灾科技学院: 45–46. Hu X H. 2020. The Study of Tectonic Stress Field in Yunnan and Its Adjacent Areas[D]. Sanhe: Institute of Disaster Prevention: 45–46 （in Chinese）.
黄小龙,吴中海,赵小艳,吴坤罡,黄小巾,杜锦锦. 2015. 2014年5月云南盈江MS5.6、MS6.1地震发震构造分析[J]. 地球学报,36(6):760–769. Huang X L,Wu Z H,Zhao X Y,Wu K G,Huang X J,Du J J. 2015. Seismogenic structure of 2014 MS5.6 and MS6.1 earthquakes in Yingjiang,Yunnan Province[J]. Acta Geoscientica Sinica,36(6):760–769 (in Chinese).
季建清,钟大赉,张连生. 2000. 滇西南新生代走滑断裂运动学、年代学、及对青藏高原东南部块体运动的意义[J]. 地质科学,35(3):336–349. doi: 10.3321/j.issn:0563-5020.2000.03.009 Ji J Q,Zhong D L,Zhang L S. 2000. Kinematics and dating of cenozoic strike-slip faults in the Tengchong area,West Yunnan:Implications for the block movement in the southeastern Tibet Plateau[J]. Scientia Geologica Sinica,35(3):336–349 (in Chinese).
罗艳, 包丰, 倪四道, 陈颙, 陈棋福, 付虹. 2011. 2008年云南盈江MS5.9地震序列初步研究[C]//中国地球物理学会第二十七届年会论文集. 长沙: 中国科学技术大学出版社: 422. Luo Y, Bao F, Ni S D, Chen Y, Chen Q F, Fu H. 2011. Study on the 2008 MS5.9 earthquake sequence in Yingjiang, Yunnan Province[C]//Proceedings of the 27th Annual Meeting of Chinese Geophysical Society. Changsha: China University of Science and Technology Press: 422 （in Chinese）.
孙尧,吴中海,安美建,龙长兴. 2014. 川滇地区主要活动断裂的活动特征及其近十年的地震活动性[J]. 地震工程学报,36(2):320–330. doi: 10.3969/j.issn.1000-0844.2014.02.0320 Sun Y,Wu Z H,An M J,Long C X. 2014. Activity characteristics of primary active faults in Yunnan-Sichuan area and their seismic activity in the past[J]. China Earthquake Engineering Journal,36(2):320–330 (in Chinese).
万永革. 2012. 根据震源机制求解构造应力场的网格搜索法[J]. 国际地震动态,(6):19. Wan Y G. 2012. The grid search algorithm of tectonic stress tensor based on focal mechanism data[J]. Recent Developments in World Seismology,(6):19 (in Chinese).
谢富仁,苏刚,崔效锋,舒赛兵,赵建涛. 2001. 滇西南地区现代构造应力场分析[J]. 地震学报,23(1):17–23. doi: 10.3321/j.issn:0253-3782.2001.01.003 Xie F R,Su G,Cui X F,Shu S B,Zhao J T. 2001. Modern tectonic stress field in southwestern Yunnan,China[J]. Acta Seismologica Sinica,14(1):17–23 (in Chinese).
徐彦. 2013. 云南地区ML3.0级以上中小地震震源机制解汇编2008—2012[M]. 云南: 云南科技出版社: 168–186. Xu Y. 2013. Decollection of Focal Mechanism of Moderate and Small Earthquakes With ML3.0 or Higher in Yunnan Area （2008−2012）[M]. Yunnan: Yunnan Science and Technology Press: 168–186 （in Chinese）.
张丽娜,韩立波,罗艳. 2016. 2014年5月云南盈江两次中强地震震源参数研究[J]. 地震,36(1):59–68. doi: 10.3969/j.issn.1000-3274.2016.01.007 Zhang L N,Han L B,Luo Y. 2016. Source parameters of May 24,2014 MS5.6 and May 30,2014 MS6.1 earthquakes in Yingjiang,Yunnan Province[J]. Earthquake,36(1):59–68 (in Chinese).
张艳凤,常祖峰,彭永中,张亚宏. 2018. 2014年云南盈江6.1级地震宏观异常及其与构造关联性讨论[J]. 地震研究,41(1):148–156. doi: 10.3969/j.issn.1000-0666.2018.01.019 Zhang Y F,Chang Z F,Peng Y Z,Zhang Y H. 2018. Discuss on macro-anomalies of the Yunnan Yingjiang M6.1 earthquake and its relation of tectonics[J]. Journal of Seismological Research,41(1):148–156 (in Chinese).
赵小艳,韩立波,龙锋. 2012. 2011年盈江MS≥4.0地震序列震源机制解与发震构造研究[J]. 地震研究,35(4):477–484. doi: 10.3969/j.issn.1000-0666.2012.04.006 Zhao X Y,Han L B,Long F. 2012. Focal mechanism solutions and seismogenic structures of Yingjiang MS≥4.0 earthquake sequences in 2011[J]. Journal of Seismological Research,35(4):477–484 (in Chinese).
Gephart J W,Forsyth D W. 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data:Application to the San Fernando earthquake sequence[J]. J Geophys Res,89(B11):9305–9320. doi: 10.1029/JB089iB11p09305
Hardebeck J L,Michael A J. 2006. Damped regional-scale stress inversions:Methodology and examples for southern California and the Coalinga aftershock sequence[J]. J Geophys Res,111(B11):B11310.
Harris R A. 1998. Introduction to special section:Stress triggers,stress shadows,and implications for seismic hazard[J]. J Geophys Res,103(B10):24347–24358. doi: 10.1029/98JB01576
Jacques A. 2002. Inversion of earthquake focal mechanisms to obtain the seismotectonic stress IV:A new method free of choice among nodal planes[J]. Geophys J Int,150(3):588–609. doi: 10.1046/j.1365-246X.2002.01713.x
Lund B,Slunga R. 1999. Stress tensor inversion using detailed microearthquake information and stability constraints:Application to Őlfus in southwest Iceland[J]. J Geophys Res,104(B7):14947–14964. doi: 10.1029/1999JB900111
Martínez-Garzón P,Kwiatek G,Ickrath M,Bohnhoff M. 2014. MSATSI:A MATLAB package for stress inversion combining solid classic methodology,a new simplified user-handling,and a visualization tool[J]. Seismol Res Lett,85(4):896–904. doi: 10.1785/0220130189
Michael A J. 1984. Determination of stress from slip data:Faults and folds[J]. J Geophys Res,89(B13):11517–11526. doi: 10.1029/JB089iB13p11517
Michael A J. 1987. Use of focal mechanisms to determine stress:A control study[J]. J Geophys Res,92(B1):357–368. doi: 10.1029/JB092iB01p00357
Vavryčuk V,Bouchaala F,Fischer T. 2013. High-resolution fault image from accurate locations and focal mechanisms of the 2008 swarm earthquakes in West Bohemia,Czech Republic[J]. Tectonophysics,590:189–195. doi: 10.1016/j.tecto.2013.01.025
Vavryčuk V. 2014. Iterative joint inversion for stress and fault orientations from focal mechanisms[J]. Geophys J Int,199(1):69–77. doi: 10.1093/gji/ggu224
-
期刊类型引用(2)
1. 向亚红，张峰，邓念武，谢雨航. 基于分形插值与支持向量机混合模型的大坝变形分析. 中国农村水利水电. 2023(01): 185-188+193 . 
百度学术
2. 李友平，王汉，程广兵，王晓波. 基于映射算法的高精度交流钳形表设计. 安徽电气工程职业技术学院学报. 2023(02): 53-58 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 735
- HTML全文浏览量: 368
- PDF下载量: 64
- 被引次数: 4
