基于区域宽频带和近震波形数据的多点源反演——以2016 M_(W)7.1日本熊本地震为例  

Multiple point source inversion with regional broadband and local waveform data-a case study of the 2016 M_(W)7.1 Kumamoto,Japan earthquake

作　　者：施其斌 韦生吉 SHI QiBin;WEI ShengJi(Asian School of the Environment,Nanyang Technological University,Singapore 639798;Earth Observatory of Singapore,Nanyang Technological University,Singapore 639798)

机构地区：[1]南洋理工大学亚洲环境学院,新加坡639798 [2]南洋理工大学新加坡地球观测中心,新加坡639798

出　　处：《地球物理学报》2021年第12期4378-4393,共16页Chinese Journal of Geophysics

基　　金：新加坡教育部项目(MOE2019-T2-1-182)资助.

摘　　要：中强度(M 6—7)地震可以对人类社会基础设施和生命财产造成威胁,并且其发震频率要远高于M 8+地震.从震源物理的角度来说,此类地震的发震断层的几何形态和破裂过程往往较为复杂.因此,准确地确定此类地震的震源参数对地震防治和了解震源物理都具有重大意义.作为承接地震单点源解和有限破裂模型的桥梁,地震多点源解可以为有限破裂模型反演提供可靠的断层几何,并有效地帮助约束破裂的时空演化过程.本文基于马尔可夫链蒙特卡洛(MCMC)非线性反演算法,将一种多点源波形反演方法的适用范围从近震距离推广到了区域距离,并验证了利用余/前震(或历史地震)对地震波传播路径进行校正的有效性.我们将改进后的方法应用于具有代表性的2016 M_(W)7.1熊本地震的研究中.结果发现,利用一个M_(W)6.0和一个M_(W)5.4的前震,可以对半径100 km范围内强震记录和1000 km内的宽频带记录进行有效的路径校正,筛选出现有速度模型所适用的地震台和震相,并且可以确定波形拟合中高质量Pnl波和面波所能利用的频率上限,即高达0.3 Hz(Pnl)和0.2 Hz(面波).研究结果还显示区域宽频带数据反演的结果和强震+高频GPS数据结果有高度的一致性,验证了将该方法推广到区域数据反演的可行性.对熊本M_(W)7.1主震的反演结果显示,在反演所用的频段,我们需要至少四个子事件来表示该地震的运动学过程.第一个M_(W)6.7子事件发生在地下12 km深的Hinagu断层上,并具有高倾角的纯走滑位错解;第二(M_(W)6.7)和第三个(M_(W)6.7)子事件的深度均为7 km,其走向和Futagawa断层一致并有一定的正断滑移成分;最后一个子事件(M_(W)6.6)最浅(1~2 km),也具有最强(滑动角=-137°)的正断滑移成分.这些子事件揭示了发震层的厚度随着主震向东北方向的破裂传播距离增加而变薄,并且其深度随着破裂位置靠近Aso火山而变得越浅.当�Moderate-size(M 6-7)earthquakes possess great threats to human society,as they occur much more frequent than M 8+earthquakes and many take place inland.The rupture processes and fault geometry of this kind of earthquake are usually very complex.Precise determination of their kinematic rupture parameters is a fundamental task in seismology,which is critical to understanding the earthquake physics as well as to better prepare for seismic hazard.As an earthquake source model with a precision between point source and finite fault model,Multiple Point Source(MPS)solution of large earthquakes provide reliable fault geometry and first order rupture evolution of earthquakes,which is extremely helpful for more detailed finite fault inversion and modeling.Here we generalize a Markov-Chain-Monte-Carlo(MCMC)based MPS inversion algorithm to the regional broadband waveform inversion and apply to the 2016 M_(W)7.1 Kumamoto mainshock.We first calibrate the path and select the frequency ranges for the strong motion(<100 km)and regional(<1000 km)broadband stations through the point source modeling of the M_(W)5.4 and M_(W)6.0 foreshocks,which results in up to 0.3 Hz for Pnl waves and 0.2 Hz for surface waves.The inversion results from regional,strong motion and high-rate GPS datasets show highly consistent solutions,validating the generalization of our MPS inversion method.Our preferred model show that the earthquake is composed of at least four sub-events.The first(M_(W)6.7)sub-event is a high-dip-angle purely strike-slip event,taking place near the hypocenter at the depth of 12 km.The strike of the first sub-event is well consistent with the strike of the Hinagu fault.The second(M_(W)6.7)and the third(M_(W)6.7)sub-events all took place at the depth of 7 km,with strike more consistent with the Futagawa fault and with slight normal slip component.The last sub-event(M_(W)6.6)is the shallowest(1~2 km)and its focal mechanism is composed of almost equal percentage of strike-slip and normal-slip components.Our five point source inversi

关 键 词：地震破裂过程 非线性反演 多点源 非对称震源时间函数 熊本地震 

分 类 号：P315[天文地球—地震学]