Study on the Static Stress Triggering Effect of the 2023 Turkey Doublet Earthquakes
-
摘要: 为了研究2023年土耳其双震间是否存在应力触发作用以及双震对周边断裂和余震的影响，本文基于USGS给出的双震破裂模型、全球震源机制解（GCMT）目录和土耳其海峡大学坎迪利天文台与地震研究所区域地震海啸监测中心实时地震资料，利用Coulomb 3.3软件从静态应力触发角度对土耳其双震序列进行了研究。研究结果表明：第1次主震在第2次主震震源处产生的库仑应力量值为0.033 MPa，超过应力触发阈值0.01 MPa，反映出第2次主震的发生受到第1次主震的触发作用。两次主震在第1次主震所在的东安纳托利亚断裂破裂段的东北部和西南部有应力加载作用，且加载的库仑应力量值较大；在7.8级地震破裂段上的作用为应力卸载，即发震段应力得到释放。两次主震在第2次主震所在的卡达克断裂的破裂段1和3交汇部位产生了应力卸载作用。余震库仑应力计算结果表明2次主震对余震存在明显的触发作用。上述研究结果可以为后续地震危险性分析等相关研究提供参考。Abstract: In order to investigate the stress triggering between the doublet Earthquakes in Turkey in 2023 and the impact of the doublet earthquakes on surrounding faults and aftershocks, this study utilized the doublet earthquakes rupture models provided by the US Geological Survey (USGS), the Global Centroid Moment Tensor (GCMT) catalog, and real-time seismic data from the Bogazici University Kandilli Observatory and Earthquake Research Institute Regional Earthquake-Tsunami Monitoring Center to calculate the static Coulomb stress by using Coulomb3.3. The results of this study show that the Coulomb stress value generated by the first mainshock at the hypocenter of the second mainshock is 0.033 MPa, exceeding the stress triggering threshold of 0.01 MPa. This indicates that the occurrence of the second mainshock is triggered by the first mainshock. The two mainshocks exert stress loading on the northeast and southwest of the East Anatolian fault, where the first mainshock ruptured. In addition, the loading Coulomb stress value is relatively large. On the rupture segment of the magnitude 7.8 earthquake, there is a stress unloading, indicating the stress release in the fault rupture segment. The two mainshocks exert stress unloading at the intersection of rupture’s segment 1 and segment 3 of the Cardak fault, where the second mainshock occurred. The calculated Coulomb stress of aftershocks indicates that the two mainshocks have a significant triggering effect on the aftershocks. The results of this study can provide references for subsequent research on earthquake hazard analysis and related studies.
-
Key words:
- Turkey earthquake /
- Doublet earthquakes /
- Static stress triggering /
- Focal mechanism /
- Seismogenic fault
1)1 2 震情简介参考https://www.scieau.com/articles/2023028636 2)2 3https://earthquake.usgs.gov/earthquakes/browse/significant.php?year=2023 3)3 4http://www.koeri.boun.edu.tr/sismo/2/latest-earthquakes/automatic-solutions/ -
引言
地形变测量中的定点潮汐形变观测主要以地倾斜、地应变等测量指标为主，其目的在于记录地壳周期性潮汐形变，并能捕捉地震中短临前兆异常信息（孙伶俐等，2013）。但实际形变观测结果受多方面因素影响，主要包括仪器系统状态、固体潮汐、观测环境、人为因素、自然环境、地震波、地下介质变化等（蔡佩蕊等，2020）。为分析和判识观测中的各类干扰信号、潮汐波及地震前兆异常信息等问题，近年来，诸多数字信号处理技术被引入形变资料分析领域，如小波分析方法（侯跃伟等，2015）、模态分解类方法（孔祥瑞等，2018）、超限率分析方法（李娜等，2020）、奇异谱分析方法（赵佳佳等，2017）、时频分析方法（王宁等，2014）等。其中，时频分析方法以其能够表征信号频谱随时间变化规律的优势，在形变观测干扰分析（张小艳等，2019）、前兆异常判定（武善艺等，2018）、日常跟踪分析（孟庆筱等，2018）等方面得到了应用与发展。目前，制约时频分析方法进一步在台站形变数据处理中的应用主要问题为：①传统的连续小波变换（CWT）、短时傅里叶变换（STFT）、S变换（ST）、希尔伯特-黄变换（HHT）等方法时频分辨率较低，对于波形类别多、分布频段广、振幅差异大的各类形变观测信号的泛化适用能力较差；②台站缺乏相关集成处理软件程序（平建军等，2013），数据导入、函数传递、方法计算、成果绘图等功能和步骤缺乏统一设计规划，不便于操作与推广。
为解决上述问题，本研究基于MATLAB科学计算平台研发了定点形变观测数据时频分析软件，该软件通过集成数据导入、插值补全、小波分析、同步挤压时频变换、自动绘图等功能，处理过程采用命令窗参数绘图交互方式，能够灵活计算特定频段、特定时段形变观测数据的高精度时频谱，以期为台站工作人员深入挖掘形变观测资料中的各类潜在变化特征提供简便实用的处理工具。
1. 方法原理
本研究采用信号分解方法中的小波分析（DWT），时频处理方法中基于CWT和广义S变换（GST）的同步挤压小波变换（SSCWT）和同步挤压广义S变换（SSGST），以及时频辅助分析方法中的叠加幅值特征函数（SCCF）、叠加系数包络函数（SCEF）和信息熵（Renyi），全面刻画不同频段分布的复杂形变观测数据时频响应特征。
DWT是利用选定小波基和分解层数以二进塔式分解方式将原始信号分解到一系列相互独立的频带上（Mallat，1989），凸显信号中不同频段成分信息，并可通过筛选重构得到特定优势频段子信号分量的方法。
CWT是将原始信号与经过伸缩平移后的母小波作褶积运算的方法，具有多尺度表达信号频率分辨率特征，其时频谱结果频率轴尺度为指数型分布，对于中低频段信号具有较好的时频响应能力。
GST是对传统ST方法高斯窗函数引入2个可调参数（魏学强等，2016），使其窗函数随频率呈现出更多非线性变化特征而提高时频分辨率的方法，其时频谱结果频率轴尺度为线性分布，对于高频段信号具有较好的时频响应能力。
SSCWT、SSGST是分别在CWT和GST“前处理”结果的基础上，进一步利用基于微分算子的挤压“后处理”操作，沿频率轴重排时频谱瞬时频带，将特定频率区间内的时频系数叠加到理论频率点上（严海滔等，2019；潘晓等，2020），可得到时频聚集性更精确的同步挤压时频结果的方法。
SCCF、SCEF是分别描述刻画同步挤压时频谱在频率尺度上的能量分布特征和在时间尺度上的强度分布特征（Mousavi等，2016）的指标。
Renyi熵是定量检测时频谱所含信息量和时频谱能量集中度的指标（黄昱丞等，2018）。
整个方法计算流程为：（1）利用DWT方法对原始形变信号作多尺度分解；（2）根据设定的研究目标，对分解出的子信号及时间范围作选择性重构，以提取待分析的目标信号分量；（3）根据目标分量的频段分布特征，基于分频处理策略，对于短周期高频目标信号和长周期低频目标信号分别采用高频响应较好的GST方法和中低频响应较好的CWT方法进行时频计算，在此基础上利用同步挤压操作得到更高分辨率的SSGST和SSCWT时频结果；（4）进一步对计算出的同步挤压时频谱结果，采用SCCF方法、SCEF方法和Renyi熵分别从频率尺度、时间尺度、能量复杂度方面进行辅助分析。
上述计算方法基本公式和控制参数如表1所示。
表 1 软件包所用方法及基本定义Table 1. Summary of methods and basic definitions of the software package类型 方法名称 基本定义 备注 信号分解方法 小波分析
（DWT）${A}_{j}s\left(t\right)={A}_{j+1}s\left(t\right)+{D}_{j+1}s\left(t\right),j=\mathrm{0,1},\cdots ,N$ $ s\left(t\right) $为原始信号，$ N $为分解层数，$ {A}_{j}s\left(t\right) $为第$ j $阶低频趋势分量，$ {D}_{j}s\left(t\right) $为第$ j $阶高频细节分量 时频前处理
方法连续小波变换
（CWT）$ C W T\left(a,\tau \right)=\left\langle{s,{\psi }_{a,\tau }}\right\rangle={\displaystyle\int }_{-\infty }^{+\infty }s\left(t\right){a}^{-\frac{1}{2}}{\psi }^{*}\left(\dfrac{t-\tau }{a}\right){\rm{d}}t $ $ {a}^{-\frac{1}{2}}\psi \left(\dfrac{t-\tau }{a}\right) $为母小波族，$ a $为小波尺度因子，$ \tau $为时间平移因子 广义S变换
（GST）$ G S T\left(f,t\right)={\displaystyle\int }_{-\infty }^{+\infty }s\left(\tau \right)\dfrac{{\left|f\right|}^{\lambda }}{\sqrt{2{\text{π}} }\rho }{\rm{e}}^{-\dfrac{{\left(\tau -t\right)}^{2}{f}^{2\lambda }}{2{\rho }^{2}}}{\rm{e}}^{-j2{\text{π}} \tau }{\rm{d}}\tau $ $ \dfrac{\left|f\right|}{\sqrt{2{\text{π}} }}{\rm{e}}^{-\dfrac{{t}^{2}{f}^{2}}{2}} $为Gauss窗，$ \lambda $为时宽调节参数，$ \rho $为衰减趋势调节参数 时频后处理
方法同步挤压小波变换（SSCWT） $ S S C W T\left({\omega }_{l},\tau \right)={\Delta \omega }^{-1}\displaystyle\sum _{{a}_{k}:\left|C W T\left(a,\tau \right)-{\omega }_{l}\right|\leqslant \Delta \omega /2}C W T\left(a,\tau \right){{a}_{k}}^{-3/2}\left({\Delta a}_{k}\right) $ $ C W T\left(a,\tau \right) $为CWT时频谱，$ {a}_{k} $为第$ k $个离散化尺度，$ {\omega }_{l} $为点$ l $处离散化频率值 同步挤压广义S变换（SSGST） $ S S G S T\left({\tilde f}_{l},t\right)={(\Delta {\tilde f}_{l})}^{-1}\displaystyle\sum _{{f}_{k}:\left|{\tilde f}_{l}\left({f}_{k},t\right)-{\tilde f}_{l}\right|\leqslant \Delta {\tilde f}_{l}/2}\left|GST\left({f}_{k},t\right)\right|{f}_{k}\Delta {f}_{k} $ $ GST\left({f}_{k},t\right) $为GST时频谱，$ {f}_{k} $为离散频点，$ {\tilde f}_{l} $为挤压中心频率，$ \Delta {\tilde f}_{l} $为挤压带宽 辅助分析
方法叠加幅值特征函数（SCCF） $ S C C F\left(t\right)=\displaystyle\sum _{t=1}^{N}\left|TFR(f,t)\right|f=1,\cdots ,M $ $ TFR(f,t)\in {\boldsymbol{C}}^{M\cdot N} $，为时频谱矩阵 叠加系数包络函数（SCEF） $ S C E F\left(f\right)=\displaystyle\sum _{f=1}^{M}\left|E(f,t)\right|t=1,\cdots ,N $ $ E(f,t) $为频率点$ f $处时频谱系数的包络函数 信息熵（Renyi） $ {H}_{a}\left(T F R\right)=\dfrac{1}{1-a}{\mathrm{log}}_{2}\displaystyle\iint T F R(f,t){\rm{d}}t{\rm{d}}f $ $ a $为Renyi熵的阶次，一般取3，$ {H}_{a} $为熵值 2. 软件设计与功能实现
2.1 总体框架
软件包总体由8个函数（脚本）程序组成（表2），按照功能可划分为：1个封装有全部功能模块函数，执行全计算与绘图流程的主脚本程序；4个分别实现数据导入、预处理、小波分解与重构、同步挤压时频分析的二级功能模块函数；3个被相应功能模块函数调用，封装有核心计算方法程序的三级子调用函数。通过在MATLAB环境下运行主脚本程序Main_func.m，依次自动执行各计算模块（图1），为增强软件包处理不同信号的效果，对于信号分解、提取重构、时频变换等关键步骤，会在命令行窗口设置关键参数输入选项和相关信息提示，用户可根据设定研究目标不同，灵活输入和选择合适参数，各模块计算完成后会立即产出中间绘图图件进行反馈，便于调参及后续操作，产出最终成果图。软件整体架构和程序代码简明高效。
图 1 软件总体处理框架流程Figure 1. General processing framework and flow of the software package表 2 软件包所含函数程序与功能Table 2. Summary of programs and functions of the software package函数（脚本）名称 所属功能模块 主要功能描述 Main_func.m 主脚本程序 封装全部功能模块函数，执行全计算流程 deformation_txt_read.m 数据导入模块 导入TXT文件中原始数据及测项信息 deformation_data_interp.m 数据预处理模块 利用DCT-PLS方法平滑插值缺失数据 deformation_mallat_wavelet.m 小波分解与重构模块 利用DWT分解并提取待分析目标信号 deformation_timefrequency_analysis.m 同步挤压时频分析模块 分频计算SSCWT或SSGST时频谱，并计算SCCF函数、
SCEF函数、Renyi熵等辅助分析项dctpls_smoothn.m 模块子调用函数 封装DCT-PLS插值方法实现算法 sscwt_tfrs.m 模块子调用函数 封装SSCWT时频方法实现算法 ssgst_tfrs.m 模块子调用函数 封装SSGST时频方法实现算法 2.2 数据导入模块
为方便统一数据接口和在台站间进行应用推广，软件包读取数据格式统一为中国地震前兆台网数据处理系统（QZProcess）软件下载的单测项TXT格式原始数据文件。用户可通过QZProcess软件选择和导出指定仪器、测项、时间范围和采样率下的观测数据，并进行时频分析。主脚本程序执行数据导入模块函数deformation_txt_read.m后，会自动打开并显示导入数据对话框，利用交互操作将数据导入MATLAB中，且相应文件信息会自动显示在命令行窗口，包括导入文件位置、台站信息、测项信息、采样率、数据类型、采样点长度、采样时间范围等。
2.3 数据预处理模块
由于时频计算的数据序列必须完整，因此执行数据预处理模块deformation_data_interp.m会对导入的原始数据进行完整性检查。若存在数据缺失，则利用离散余弦表示的惩罚最小二乘回归（DCT-PLS）方法对缺失数据进行无参数自动平滑插值补全（张桂欣等，2016），处理完成后会在命令行窗口自动显示是否缺数及缺失点位置序号，同时产出导入原始数据及插值结果绘图图件。
2.4 小波分解与重构模块
为提取观测数据中需分析的目标信号分量，对预处理后完整数据执行小波分解与重构模块deformation_mallat_wavelet.m，将信号中不同频率成分进行分离与重组。首先，用户需根据导入数据绘图和显示信息特征，指定小波分解参数，即在命令行窗口选择和输入小波分解层数和小波基函数，执行小波分解计算，产出原始数据小波分解结果绘图图件，在定点形变数据分析中，常将分解层数设为4～10层（张中旭等，2020），小波基函数常选用db族小波、sym族小波等（刘建明等，2016）。然后，根据小波分解结果中目标信号所在频段和时间范围分布特征，指定小波重构参数，即以首尾序号方式输入小波重构（叠加）区间和提取目标时段范围，提取目标信号，并产出信号提取结果绘图图件。
2.5 同步挤压时频计算模块
最后执行同步挤压时频计算模块deformation_timefrequency_analysis.m，对小波提取的目标信号进行处理。基于分频处理策略，软件包内置2种同步挤压时频分析方法，用户可根据目标信号频段特征自主选择采用高频分辨率较好的SSGST方法或低频分辨率较好的SSCWT方法进行时频谱计算。选择SSGST方法时，需选择和输入决定SSGST方法时频成像分辨率的高斯窗参数值，该参数由时宽调节参数和衰减趋势调节参数共同决定，一般取值范围为0～1（杨千里等，2020），可根据实际处理信号特征进行灵活调整，达到最优时频成像效果。选择SSCWT方法时，需选择和输入主要影响SSCWT方法结果的母小波类型参数和小波尺度控制参数，常用母小波包括morlet小波、cmor小波等（张维辰等，2019），小波尺度控制参数决定了结果在尺度方向上采样点数和分辨率。为统一不同采样率数据和不同计算方法的单位，时频结果中频率轴尺度采用了归一化频率。进一步利用3种辅助分析方法对计算出的同步挤压时频谱进行统计分析，最终产出目标信号时频分析成果绘图图件。
成果图件上部为显示导入原始数据及提取目标信号的时间范围（红色虚线内），下部为显示提取目标信号及其同步挤压时频谱、SCCF特征谱和SCEF包络谱，可更全面地描述观测数据中研究目标对象的时频响应特征及其在频率、时间尺度上的变化趋势，如图2所示（以西昌小庙台2021年3月DSQ型水管仪NS向观测数据为例）。
3. 应用实例
为验证本文所提定点形变观测数据时频分析软件包的实用性，利用该软件对2020年1月1日至2021年6月30日西昌小庙台DSQ型水管仪、SS-Y型伸缩仪4个测项整时值采样观测数据进行处理，提取原始数据中固体潮频段分量信号进行时频分析，以评估观测质量。小波分解层数设为4层，小波基函数采用db5小波。小波重构区间设为细节1～细节4，即提取周期为2～32 h的信号，并将提取目标时段范围设为整个采样时间。时频计算选择SSCWT方法，并将母小波类型设为morlet小波。通过绘制的最终成果图进行对比分析。
3.1 DSQ型水管仪固体潮分析
从时频谱角度考察观测数据精度，主要依据为观测条件好、仪器运行稳定的测项可获得信噪比强的固体潮信号，其反映在时频谱上为能量主要集中在优势波群频段（吕品姬等，2011）。一般通过时频谱视觉上的定性对比确定时频能量集中性的优劣，本研究为定量评估时频谱所含信息量和复杂度，引入了Renyi熵作为参考指标，熵值越小，表明信号在时频域的能量集中度越高，反之越低。
由处理结果可知，研究时段内西昌小庙台DSQ型水管仪NS分量观测固体潮数据（图3（a））与EW分量观测固体潮数据（图3（b））在时域变化形态上均较规律，潮汐波形态清晰，时频谱中半日波、全日波等信号频率能量分布集中，除部分由地震波、降雨造成的突跳信号外，基本不存在非固体潮频段的杂波信号能量，将观测固体潮时频结果与理论固体潮时频结果（图3（c）、图3（d））进行对比，发现西昌小庙台DSQ型水管仪观测数据在时域信号形态、时频谱能量分布、SCCF特征谱频率分布、SCEF包络谱强度变化等方面一致性较好，观测固体潮时频谱Renyi熵值较小，且接近理论固体潮熵值（表3），总体上，西昌小庙台DSQ型水管仪观测精度较好，观测质量较高。
图 3 小庙台2020年1月至2021年6月DSQ型水管仪观测固体潮与理论固体潮时频结果对比Figure 3. Comparison of time-frequency results between observed and theoretical earth tide of DSQ data of Xiaomiao seismic station from Jan, 2020 to June,2021表 3 观测固体潮与理论固体潮时频谱Renyi熵值Table 3. The Renyi entropy of observed earth tide and theoretical earth tide项目 DSQ（NS向） DSQ（EW向） SS-Y（NS向） SS-Y（EW向） 观测固体潮时频谱Renyi熵值 6.26 5.56 8.44 8.17 理论固体潮时频谱Renyi熵值 5.50 5.50 5.13 5.13 Renyi熵误差 0.76 0.06 3.31 3.04 3.2 SS-Y型伸缩仪固体潮分析
西昌小庙台SS-Y型伸缩仪NS分量固体潮（图4（a））与EW分量固体潮（图4（b））在研究时段内信号存在较大噪声干扰，时频谱中主要波群信号能量较分散，且被强杂波能量覆盖。由SCCF特征谱可知，相较于理论固体潮（图4（c）、图4（d））的频率分布，观测数据频率能量明显分散在固体潮中心频率以外的其他频带上，且时频谱Renyi熵值偏大，说明西昌小庙台SS-Y型伸缩仪观测精度较差，数据质量有待提高。这是因为近年来西昌小庙台应变仪器相比倾斜仪器更易受到台站周边施工、交通、机井抽水等观测环境干扰，且传感器老化严重，导致日常观测数据常出现高频噪声和波形畸变。
图 4 小庙台2020年1月至2021年6月SS-Y型伸缩仪观测固体潮与理论固体潮时频结果对比Figure 4. Comparison of time-frequency results between observed and theoretical earth tide of SS-Y data of Xiaomiao seismic station from Jan,2020 to June,20214. 结论
本研究所提基于MATLAB的定点形变观测数据时频分析软件包采用了小波分析、同步挤压时频分析和时频辅助分析等多种方法，程序基于模块化集成封装原则，以绘图-参数输入-绘图交互运行方式，实现了提取原始形变观测数据中特定频段、特定时段的目标信号分量，并对其进行高精度时频计算与自动绘图。将该软件包应用于形变观测资料的固体潮数据分析，通过对比理论固体潮，可从时频谱角度评估观测数据质量，具有一定推广应用价值。
-
图 1 研究区地质构造与地震分布图
Figure 1. Geological structure and earthquake distribution of the study area
图 2 M7.8地震在M7.5地震破裂面上产生的库仑应力分布
Figure 2. The Coulomb stress change distribution on the rupture plane of the M7.5 earthquake caused by the M7.8 mainshock
图 3 第1次主震在周边断裂上产生的库仑应力变化图
Figure 3. The coulomb stress change on the surrounding active faults caused by the first mainshock
图 4 双震在周边断裂上产生的库仑应力变化图
Figure 4. The coulomb stress change on the surrounding active faults caused by the doublet earthquakes
图 5 双震产生的库仑应力及余震分布图
Figure 5. The distribution of aftershocks and coulomb stress caused by the doublet earthquakes
表 1 震源机制解参数表
Table 1. The parameter table of focal mechanisms
序号 发震时间/
年-月-日纬度/（°） 经度/（°） 震级MW 深度/km 节面1 节面2 走向/（°） 倾角/（°） 滑动角/（°） 走向/（°） 倾角/（°） 滑动角/（°） 1 1979-12-28 37.47 35.85 5.4 41.0 231 90 0 141 90 180 2 1986-05-05 37.97 37.77 6.0 10.0 260 54 9 164 82 144 3 1986-06-06 37.97 37.88 5.8 10.0 250 90 0 160 90 180 4 1989-06-24 36.70 35.88 5.1 41.0 27 62 −88 203 28 −93 5 1991-04-10 37.21 36.01 5.3 33.0 29 72 −70 160 27 −136 6 1997-01-22 36.25 35.95 5.7 10.0 243 39 −15 345 81 −128 7 1998-05-09 38.28 38.99 5.1 10.0 251 83 −7 341 83 −173 8 1998-06-27 36.88 35.31 6.3 5.8 53 81 15 321 75 171 9 1998-07-04 36.877 35.32 5.4 33.0 72 55 8 338 84 145 10 2001-06-25 37.24 36.21 5.4 5.0 1 75 −92 189 15 −83 11 2003-07-13 38.29 38.96 5.5 10.0 72 89 1 342 89 179 12 2005-11-26 38.26 38.81 5.1 8.5 237 51 −20 339 75 −139 13 2006-03-29 35.25 35.43 5.0 27.3 219 43 −10 317 83 −132 14 2008-09-03 37.51 38.50 5.0 5.7 219 79 −10 311 80 −169 15 2008-11-12 38.84 35.52 5.1 10.0 227 70 −13 321 78 −160 16 2010-11-14 36.58 36.01 4.9 2.5 24 53 −94 211 37 −84 17 2012-07-22 37.55 36.38 4.8 7.6 38 53 −78 198 38 −106 18 2012-09-19 37.31 37.10 5.0 7.0 210 48 −11 307 82 −138 19 2014-02-14 36.74 36.08 4.9 10.0 35 70 −59 155 36 −144 20 2014-06-09 36.74 36.05 4.8 17.6 34 65 −63 164 36 −135 21 2015-11-29 38.82 37.74 5.1 10.0 338 72 161 74 72 19 22 2017-03-02 37.62 38.43 5.6 10.0 225 78 −21 319 69 −167 23 2018-04-24 37.60 38.51 5.2 10.0 212 63 −3 304 87 −153 24 2018-10-02 37.67 37.40 4.7 5.0 242 90 −42 332 48 −180 25 2019-03-25 38.69 38.07 4.9 10.0 343 64 −157 242 70 −28 26 2020-02-25 38.34 38.80 5.0 10.0 232 40 −29 345 72 −127 27 2020-04-03 35.94 35.49 4.8 12.7 345 50 −117 204 47 −61 28 2020-04-15 35.86 35.53 4.9 10.0 219 47 −10 316 83 −137 29 2020-06-05 38.24 38.76 5.1 10.0 234 49 −24 340 72 −137 30 2020-08-04 38.19 38.70 5.6 10.0 235 75 −15 329 76 −165 31 2020-09-08 38.06 38.78 4.8 10.0 238 63 −7 332 83 −153 32 2021-11-12 38.20 38.78 5.0 7.0 237 82 −25 330 65 −171 33 2022-04-09 38.11 38.67 5.3 10.0 248 83 −14 340 77 −173 34 2022-10-11 37.26 36.23 5.0 10.0 17 52 −80 181 39 −103 35 2023-02-06 37.17 37.03 7.8 17.9 54 70 11 320 80 160 36 2023-02-06 37.13 36.94 6.8 14.5 211 66 −13 306 78 −156 37 2023-02-06 38.02 37.20 7.7 10.0 261 42 −8 358 84 −132 38 2023-02-07 37.76 37.74 5.5 10.0 204 54 10 108 80 144 39 2023-02-08 37.95 37.65 5.5 7.5 206 74 177 297 87 16 40 2023-02-20 36.11 36.02 6.3 16.0 227 45 −16 329 79 −134 
下载: 导出CSV
-
程佳, 徐锡伟, 2018. 巴颜喀拉块体周缘强震间应力作用与丛集活动特征初步分析. 地震地质, 40（1）: 133—154 doi: 10.3969/j.issn.0253-4967.2018.01.011Cheng J. , Xu X. W. , 2018. Features of earthquake clustering from calculation of Coulomb stress around the Bayan Har Block, Tibetan Plateau. Seismology and Geology, 40(1): 133—154. (in Chinese) doi: 10.3969/j.issn.0253-4967.2018.01.011 郝平, 傅征祥, 田勤俭等, 2004. 昆仑山口西8.1级地震强余震库仑破裂应力触发研究. 地震学报, 26（1）: 30—37Hao P. , Fu Z. X. , Tian Q. J. , et al. , 2004. Large aftershocks triggering by Coulomb failure stress following the 2001 MS=8.1 great Kunlun earthquake. Acta Seismologica Sinica, 26(1): 30—37. (in Chinese) 何金, 许鑫, 吴彪等, 2022. 汶川地区区域构造应力场特征及汶川MS8.0地震对周围主要断层面的影响. 地震科学进展, 52（3）: 97—115He J. , Xu X. , Wu B. , et al. , 2022. Characteristics of regional tectonic stress field in Wenchuan area and the effect of Wenchuan MS8.0 earthquake on surrounding major fault surfaces. Progress in Earthquake Sciences, 52(3): 97—115. (in Chinese) 贾若, 蒋海昆, 2014. 基于同震库仑应力变化的汶川地震余震频次研究. 中国地震, 30（1）: 74—90Jia R. , Jiang H. K. , 2014. Aftershock sequence frequency research on the Wenchuan aftershock sequence based on the coseismic Coulomb stress changes. Earthquake Research in China, 30(1): 74—90. (in Chinese) 靳志同, 崔华伟, 刘佳璐等, 2023.2023年土耳其两次强震对周围地区的静态应力影响. 防灾科技学院学报, 25（2）: 1—12Jin Z. T. , Cui H. W. , Liu J. L. , et al. , 2023. Impact of two strong earthquakes in Turkey in 2023 on the static stress in the surrounding areas of the epicenters. Journal of Institute of Disaster Prevention, 25(2): 1—12. (in Chinese) 李健, 詹文欢, 朱俊江等, 2016. 南海东部俯冲洋脊区段地震静态库伦应力分析. 海洋地质与第四纪地质, 36（5）: 63—73Li J. , Zhan W. H. , Zhu J. J. , et al. , 2016. Static Coulomb stress analysis for earthquakes in the spreading ridge of the South China Sea. Marine Geology & Quaternary Geology, 36(5): 63—73. (in Chinese) 李健, 詹文欢, 朱俊江等, 2017.1990年菲律宾Mw7.7级强震对马尼拉俯冲带静态应力触发影响. 海洋地质与第四纪地质, 37（6）: 93—99Li J. , Zhan W. H. , Zhu J. J. , et al. , 2017. A preliminary study on static stress triggering effects on Manila subduction zone by the Philippine Mw7.7 earthquake 1990. Marine Geology & Quaternary Geology, 37(6): 93—99. (in Chinese) 李玉江, 石富强, 张辉等, 2020. 川滇地区主要断裂带上的库仑应力变化及其对地震危险性的指示. 地震地质, 42（2）: 526—546Li Y. J. , Shi F. Q. , Zhang H. , et al. , 2020. Coulomb stress change on active faults in Sichuan-Yunnan region and its implications for seismic hazard. Seismology and Geology, 42(2): 526—546. (in Chinese) 刘盼, 李平恩, 廖力, 2017. 从库仑破裂应力和余震分布角度探讨汶川地震和芦山地震的关系. 震灾防御技术, 12（1）: 40—55Liu P. , Li P. E. , Liao L. , 2017. Discussion of relationship between the Wenchuan earthquake and Lushan earthquake from the viewpoint of Coulomb failure stress change and spatial distribution of aftershocks. Technology for Earthquake Disaster Prevention, 12(1): 40—55. (in Chinese) 刘强, 倪四道, 秦嘉政等, 2007.2007年宁洱6.4级地震强余震库仑破裂应力触发研究. 地震研究, 30（4）: 331—336Liu Q. , Ni S. D. , Qin J. Z. , et al. , 2007. Triggered strong aftershock by Coulomb failure stress change caused by the 2007 Ning'er, Yunnan, MS6.4 earthquake. Journal of Seismological Research, 30(4): 331—336. (in Chinese) Seismology小组, 2023. 2023年2月6日土耳其双强震的震源机制中心解、地震触发关系及产生的地表变形. （2023-02-08）[2023-07-06]. https://mp.weixin.qq.com/s/Au_fjta-kqAuHmwkp4LUw.Seismology Group, 2023. The center focal mechanism solution of the 2023 Türkiye Doublet Earthquakes on February 6, 2023, and their static stress triggering relationship and surface deformation. （2023-02-08）[2023-07-06].https://mp.weixin.qq.com/s/Au_fjta-kqAuHmwkp4LUw. (in Chinese) 盛书中, 万永革, 程佳等, 2012.2011年日本9.0级大地震的应力触发作用初步研究. 地震地质, 34（2）: 325—337Sheng S. Z. , Wan Y. G. , Cheng J. , et al. , 2012. Primary research on the Coulomb stress triggering of the 2011 Mw 9.0 Tohoku earthquake. Seismology and Geology, 34(2): 325—337. (in Chinese) 盛书中, 万永革, 蒋长胜等, 2015.2015年尼泊尔MS8.1强震对中国大陆静态应力触发影响的初探. 地球物理学报, 58（5）: 1834—1842Sheng S. Z. , Wan Y. G. , Jiang C. S. , et al. , 2015. Preliminary study on the static stress triggering effects on China mainland with the 2015 Nepal MS8.1 earthquake. Chinese Journal of Geophysics, 58(5): 1834—1842. (in Chinese) 盛书中, 万永革, 徐锡伟等, 2019. 用大量地震震源机制解检验2008年汶川地震对后续地震的触发作用[J]. 地球物理学报, 62（12）: 4588—4603 doi: 10.6038/cjg2019M0198Sheng S. Z. , Wan Y. G. , Xu X. W. , et al. , 2019. Using a large number of focal mechanism solutions to examine the Coulomb stress triggering effect of the 2008 Wenchuan earthquake on its subsequent earthquakes. Chinese Journal of Geophysics, 62(12): 4588—4603. (in Chinese) doi: 10.6038/cjg2019M0198 石富强, 张辉, 邵志刚等, 2020. 华北地区库仑应力演化与强震活动关系. 地球物理学报, 63（9）: 3338—3354Shi F. Q. , Zhang H. , Shao Z. G. , et al. , 2020. Coulomb stress evolution and stress interaction among strong earthquakes in North China. Chinese Journal of Geophysics, 63(9): 3338—3354. (in Chinese) 宋金, 周龙泉, 2014.2014年于田MS7.3地震产生的静态库仑应力变化及对周边断层的影响. 中国地震, 30（2）: 168—177Song J. , Zhou L. Q. , 2014. The static stress triggering effects related with the Yutian MS7.3 earthquake. Earthquake Research in China, 30(2): 168—177. (in Chinese) 万永革, 吴忠良, 周公威等, 2000. 几次复杂地震中不同破裂事件之间的“应力触发”问题. 地震学报, 22（6）: 568—576Wan Y. G. , Wu Z. L. , Zhou G. W. , et al. , 2000. "Stress triggering" between different rupture events in several earthquakes. Acta Seismologica Sinica, 22(6): 568—576. (in Chinese) 汪建军, 许才军, 2017.2017年Mw6.5九寨沟地震激发的同震库仑应力变化及其对周边断层的影响. 地球物理学报, 60（11）: 4398—4420Wang J. J. , Xu C. J. , 2017. Coseismic Coulomb stress changes associated with the 2017 Mw6.5 Jiuzhaigou earthquake (China) and its impacts on surrounding major faults. Chinese Journal of Geophysics, 60(11): 4398—4420. (in Chinese) 解朝娣, 王英楠, 闫如玉等, 2021.2014年四川康定MS6.3和MS5.8地震的应力触发研究. 地震研究, 44（1）: 9—14Xie C. D. , Wang Y. N. , Yan R. Y. , et al. , 2021. Study on stress triggering of 2014 Kangding MS6.3 and MS5.8 earthquakes. Journal of Seismological Research, 44(1): 9—14. (in Chinese) 熊维, 谭凯, 刘刚等, 2015.2015年尼泊尔Mw7.9地震对青藏高原活动断裂同震、震后应力影响. 地球物理学报, 58（11）: 4305—4316Xiong W. , Tan K. , Liu G. , et al. , 2015. Coseismic and postseismic Coulomb stress changes on surrounding major faults caused by the 2015 Nepal Mw7.9 earthquake. Chinese Journal of Geophysics, 58(11): 4305—4316. (in Chinese) 尹迪, 董培育, 石耀霖, 2022. 紫坪铺库区小震产生的库仑应力变化及其与汶川地震的关系. 地球物理学报, 65（1）: 256—267Yin D. , Dong P. Y. , Shi Y. L. , 2022. Coulomb stress changes induced by small earthquakes in the Zipingpu reservoir area and its significance to the Wenchuan earthquake. Chinese Journal of Geophysics, 65(1): 256—267. (in Chinese) 尹凤玲, 蒋长胜, 韩立波等, 2018. 红河断裂带库仑应力演化及未来地震危险性估计. 地球物理学报, 61（1）: 183—198Yin F. L. , Jiang C. S. , Han L. B. , et al. , 2018. Seismic hazard assessment for the Red River fault: insight from Coulomb stress evolution. Chinese Journal of Geophysics, 61(1): 183—198. (in Chinese) 张迎峰, 张国宏, 单新建等, 2017.2015年尼泊尔Gorkha Mw 7.9地震与Kodari Mw 7.3地震InSAR数据反演及其应力触发分析. 地震地质, 39（1）: 104—116Zhang Y. F. , Zhang G. H. , Shan X. J. , et al. , 2017. The coseismic source slip and Coulomb stress triggering of 2015 Nepal Gorkha Mw 7.9 and Kodari Mw 7.3 earthquake based on InSAR measurements. Seismology and Geology, 39(1): 104—116. (in Chinese) 周龙泉, 马宏生, 夏红等, 2008.2007年苏门答腊8.5级、8.3级地震强余震库仑破裂应力触发研究. 地震, 28（1）: 40—46Zhou L. Q. , Ma H. S. , Xia H. , et al. , 2008. Large aftershocks triggering by coulomb failure stress following the 2007 MS8.5 and 8.3 Sumatra great earthquakes. Earthquake, 28(1): 40—46. (in Chinese) 周云, 潘正洋, 王卫民等, 2021.1998年以来伽师地震（Mw≥6.0）应力相互作用与强震活动的关系. 地震地质, 43（2）: 280—296Zhou Y. , Pan Z. Y. , Wang W. M. , et al. , 2021. Relationship between stress interaction and strong earthquake activity of Jiashi earthquakes (Mw≥6.0) since 1998. Seismology and Geology, 43(2): 280—296. (in Chinese) Cotton F. , Coutant O. , 1997. Dynamic stress variations due to shear faults in a plane-layered medium. Geophysical Journal International, 128(3): 676—688. doi: 10.1111/j.1365-246X.1997.tb05328.x Hardebeck J. L. , Nazareth J. J. , Hauksson E. , 1998. The static stress change triggering model: constraints from two southern California aftershock sequences. Journal of Geophysical Research: Solid Earth, 103(B10): 24427—24437. doi: 10.1029/98JB00573 Harris R. A. , 1998. Introduction to special section: stress triggers, stress shadows, and implications for seismic hazard. Journal of Geophysical Research: Solid Earth, 103(B10): 24347—24358. doi: 10.1029/98JB01576 King G. C. P. , Stein R. S. , Lin J. , 1994. Static stress changes and the triggering of earthquakes. Bulletin of the Seismological Society of America, 84(3): 935—953. Lin J. , Stein R. S. , 2004. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults. Journal of Geophysical Research: Solid Earth, 109(B2): B02303. Nalbant S. S. , Hubert A. , King G. C. P. , 1998. Stress coupling between earthquakes in northwest Turkey and the north Aegean Sea. Journal of Geophysical Research: Solid Earth, 103(B10): 24469—24486. doi: 10.1029/98JB01491 Papadimitriou E. E. , Karakostas V. G. , Papazachos B. C. , 2001. Rupture zones in the area of the 17.08. 99 Izmit (NW Turkey) large earthquake (Mw7.4) and stress changes caused by its generation. Journal of Seismology, 5(2): 269—276. doi: 10.1023/A:1011463420557 Pauchet H. , Rigo A. , Rivera L. , et al. , 1999. A detailed analysis of the February 1996 aftershock sequence in the eastern Pyrenees, France. Geophysical Journal International, 137(1): 107—127. doi: 10.1046/j.1365-246x.1999.00776.x Pinar A. , Honkura Y. , Kuge K. , 2001. Seismic activity triggered by the 1999 Izmit earthquake and its implications for the assessment of future seismic risk. Geophysical Journal International. , 146(1): F1—F7. doi: 10.1046/j.0956-540x.2001.01476.x Stein R. S. , Barka A. A. , Dieterich J. H. , 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering. Geophysical Journal International, 128(3): 594—604. doi: 10.1111/j.1365-246X.1997.tb05321.x Toda S. , Stein R. S. , Reasenberg P. A. , et al. , 1998. Stress transferred by the 1995 Mw=6.9 Kobe, Japan, shock: effect on aftershocks and future earthquake probabilities. Journal of Geophysical Research: Solid Earth, 103(B10): 24543—24565. doi: 10.1029/98JB00765 Toda S. , Stein R. S. , Richards-Dinger K. , et al. , 2005. Forecasting the evolution of seismicity in southern California: animations built on earthquake stress transfer. Journal of Geophysical Research: Solid Earth, 110(B5): B05 S16. Toda S., Lin J., Stein R. S., 2011. Using the 2011 Mw 9.0 off the pacific coast of Tohoku earthquake to test the Coulomb stress triggering hypothesis and to calculate faults brought closer to failure. Earth, Planets and Space, 63(7): 725—730. Toda S., Stein R. S., Özbakir A. D., et al., 2023. Stress change calculations provide clues to aftershocks in 2023 Türkiye earthquakes. (2023-02-08)[2023-07-06]. http://doi.org/10.32858/temblor.295. Wan Y. G., Wu Z. L., Zhou G. W., 2004. Focal mechanism dependence of static stress triggering of earthquakes. Tectonophysics, 390(1—4): 235—243. 期刊类型引用(1)
1. 侯跃伟，刘孝峰，杜存鹏. 定点形变数据可视化分析系统的研制及应用. 地震地磁观测与研究. 2023(06): 154-162 . 
百度学术其他类型引用(0)
-
下载:
点击查看大图
计量
- 文章访问数: 143
- HTML全文浏览量: 36
- PDF下载量: 12
- 被引次数: 1
