2015年四川乐山金口河M5.0地震孕震构造条件分析

2015年1月14日乐山金口河M5.0地震发生在历史地震强度较低的川南山区与四川盆地交界一带。基于四川区域地震台网的震相报告与波形资料,采用双差定位法对地震序列进行重新定位,同时,采用CAP波形反演方法及HASH方法反演了主震及序列中8次M_L≥2.0地震的震源机制解。另外,利用Coulomb3计算了主震发生后库仑应力改变量,得到的结果如下：①重新定位结果显示,金口河M5.0地震位于（103.18°E,29.32°N）,震源深度16.6km,略深于波形反演结果（12km）。序列分布在NNW向天全-荥经断裂和NE向西河-美姑断裂的交汇部位,余震序列在空间上呈NE向展布。②M5.0主震的机制解为节面Ⅰ：走向350°/倾角46°/滑动角107°,节面Ⅱ：走向146°/倾角47°/滑动角73°,表现为走向NW（NNW）、中等倾角的逆冲型运动方式。序列中其余8次M_L≥2.0余震大多以走向NE的逆冲型地震为主,个别为走滑或正断层类型。主震和大部分余震的节面方向不一致,主震节面方向与余震长轴方向也不一致。③主震后库仑应力改变量显示,余震主要发生在主震引起的库仑破裂应力增加的区域。综合分析推测,NNW向天全-荥经断裂为本次地震主震的发震构造,倾向NE的机制解节面Ⅰ指出了该断裂的几何产状;M5.0主震发生后,立即触发了其旁侧的NE向西河-美姑断裂,并激发了多次余震。     

网格内不满足均匀性假设对应力场反演结果的影响——以喜马拉雅东构造结及其周边地区应力场研究为例

本研究搜集了GCMT中自1976年1月1日至2017年4月30日喜马拉雅东构造结及其周边区域的102个震源机制解数据,应用MSATSI软件反演了该区域的构造应力场.反演结果给出了该区域应力场的总体特征:最大主压应力轴自西向东呈现NNE-NE偏转,最小主压应力轴自西向东呈现NWW-NW-NNW扇形偏转;整个研究区域σ_3轴倾角较小接近水平,σ_1轴倾角情况复杂,大体上自西向东由水平转向直立;R值自西北至东南出现规律性递增现象.根据反演结果,部分网格点中应力轴倾角的置信区间过大,且不同应力轴的置信区间存在重叠现象,同时该部分网格内震源机制解类型多样.针对上述现象,本文利用两种反演方法对单个网格基于不同类型震源机制解数据进行单独反演、利用人工合成震源机制数据进行验证以及各震源机制解相对于参考应力张量的残差分析等方法对其进行了深入的验证和分析,分析结果表明导致这一现象的原因可能是单个网格内应力场不均匀所致,此时应力场反演结果为多个应力场的综合.本研究结果可为后续应力场研究中网格的合理划分以及应力场结果评价提供参考.     

安徽地区震源参数的相关性研究

基于安徽数字地震台网2010年1月至2017年12月记录到的ML2.5~5.0级地震,根据中小地震震源参数测定原理,利用多台多地震联合反演Atkinson方法和Moya方法分别计算了安徽地区地壳平均非弹性衰减因子Q值和安徽24个省属专业地震的台站场地。在此基础上,运用遗传算法求得安徽地区96个地震事件的震源谱参数,进而根据Brune中小地震圆盘模型计算其拐角频率、地震矩、应力降、矩震级、震源尺度等震源参数,并分析其特征及相互之间的关系。研究结果表明:安徽地区中小地震的ML震级与其他震源参数之间存在一定相关关系,而地震矩M0与应力降Δσ、震源尺度r和拐角频率fc之间并未表现出明显的相关关系。     

基于sPL震相测定山西地区的震源深度

本文利用sPL震相计算了山西及周边地区50 km之内的39次ML≥3.0地震的震源深度。结果显示:震源深度分布在6~30 km,优势震源深度为11~25 km,呈现出山西北部地震的震源深度比中、南部浅的趋势。将本文得到的39个地震的震源深度与中国地震台网中心统一编目的深度结果相比较,发现除极个别的地震事件外,两种方法得出的深度结果差值不大。     

2017年8月9日精河MS6.6地震余震序列精定位及发震构造分析

为确定2017年8月9日精河M_S6.6地震的发震构造,本文使用双差定位方法对发震时刻至2017年10月震源区所发生的余震进行了精定位,同时利用CAP波形反演方法,得到了主震的震源机制解,同时使用GPAT方法反演得到了部分余震的震源机制解,并基于两者对本次地震的发震构造予以分析。结果显示:精定位后主震位于（44.27°N,82.85°E）,震源深度为17 km;主震最佳双力偶解对应的节面Ⅰ的走向为260°、倾角为51°、滑动角为84°,节面Ⅱ的走向为89.5°、倾角为39.4°、滑动角为97.4°;余震序列位于主震东侧,并向东展布约30 km,在3—18 km深度范围内均有分布,其优势方向为近EW向,次优势方向为SW向。结果表明,本次地震是一次逆冲型地震,通过反演得到的大量小震震源机制解的结果与主震震源机制解结果相一致。结合余震震中分布、主震及余震的震源机制解以及震源区的地质构造,本文推断近EW走向具有逆冲性质的库松木楔克山前断裂为精河主震的发震构造。      

四川长宁MS6.0地震区域构造应力场特征分析

2019年6月17日在四川宜宾市长宁县（28.34°N,104.90°E）发生M_S6.0地震,余震发育。本文利用区域测震台网的地震观测数据基于CAP方法计算了28°~29°N,104°~105°E范围内的14个M_S>3.0以上地震的震源机制解,结合全球矩心矩张量目录和部分前人研究结果中该区域的共27个震源机制解数据,应用MSATSI软件反演了研究区域的应力场。将研究区域按0.1°×0.1°划分成25个应力网格,最终得到9个网格的应力分布结果,大多数应力场方向稳定,根据主震所在应力网格点得到主震的断层类型为主逆冲型。本文研究成果为四川长宁地区的孕震机理、活动构造以及地震趋势判定提供了可靠的参考依据。     

2016年青海门源Ms6.4级地震多个震源机制中心解的确定

2016年1月21日,青海省门源县发生Ms6.4级地震。不同机构、学者依据不同方法和资料计算得到该地震的震源机制解,由于采用的方法和选择的资料不同,所求解的震源机制解存在着一定的差别。为了在多个震源机制解中确定一个震源机制解进行后续分析,采用最小空间旋转角表达了2个震源机制的差别,将待定震源机制与所有震源机制的最小空间旋转角的平方和作为目标函数,并且采用Levenberg-Marquardt方法求解,使其与所有收集的震源机制解的差别的平方最小,这些震源机制结果就是多个震源机制解的中心解,同时给出标准差,选取标准差最小的中心解作为最终结果。根据这种方法,已知多个该地震的震源机制中确定中心解,从而对其他诸如应力触发、地球动力学分析等提供后续分析。     

晋冀蒙交界地区最小一维速度模型联合反演

文章使用内蒙古地震台网记录到的晋冀蒙交界地区517个定位地震,在参考原始编目报告的基础上,用和达法拟合走时曲线删除错误数据,使用Kissling方法联合反演得到晋冀蒙交界地区最小一维P波速度模型和震源位置。结果显示,研究区P波速度随着深度增加而递增,无明显间断面,联合反演得到的震中位置较原始报告分布更集中,簇状分布更明显。     

2017年9月4日河北临城M_L4.4地震震源深度测定

震源深度是地震定位中的重要参数之一。对于一个地震,sPn与Pn震相到时差是一个常量,与震中距无关,只与震源深度和地壳模型有关。基于该理论,以河北数字台网数据为基础,识别sPn震相和Pn震相,计算2017年9月4日河北临城ML4.4地震震源深度,结果为6.86km。双差定位法是目前确定震源深度比较精确的方法,双差定位得到的震源深度为7.01km。双差定位法所用数据量比较大,计算过程繁杂。利用sPn与Pn到时差计算震源深度的结果与双差定位法结果接近,计算快速、简单,所需数据量小。因此,将该方法应用到实际工作中,可以快速、准确地测定近震震源深度。     

基于三维应变格林函数反演中小震震源机制

随着复杂速度结构反演的发展和高性能计算能力的提升,基于高精度3D介质模型计算格林函数反演震源机制更具可行性。中小地震因具有更好的覆盖和近似点源效应,在区域结构成像中有着广泛的应用。基于波形类的反演方法如波动方程层析成像\,全波形反演都需要震源机制解,而传统的震源机制反演方法不能很好地应用于中小地震。采用有限差分法构建应变格林张量(Strain Green Tensor,SGT)数据库,将合成波形和实际波形按震相截窗并滤波到不同的频带范围,先通过最小化互相关走时差来进行震源重定位,再通过最小化波形残差反演震源机制。通过合成数据测试验证方法的正确性,随后将该方法应用于青藏高原东部边缘龙门山地区,反演一系列M_W3.4～5.7的中小地震震源机制。由于应变格林张量数据库可预先构建,该方法可以应用于(近)实时震源机制解反演。