日本熊本M_W 7.0地震的长周期地震动

利用日本K-NET和KiK-net强震动台网获取的距离发震断层100 km以内136个强震动台站的三分量加速度记录,研究熊本M_W7.0地震地震动的长周期特性.基于残差分析研究不同周期地震动的空间分布差异,将观测分析结果与美国NGA经验模型、汶川和芦山地震观测结果进行对比,揭示此次熊本地震近场强震动的长周期特点及其形成机理.研究结果表明:(1)虽然总体上此次地震的近场地震动水平与美国NGA-West2经验模型的预测结果接近,但周期2 s以上地震动的分布在断层不同方位有系统性差异,在断层的北东方位,周期2.0～10.0 s的反应谱高于NGA-West2经验模型的预测结果,在西南方位,谱值低于经验预测模型.(2)我们认为此次地震2.0～10.0 s的长周期地震动的空间分布差异主要受破裂方向性的影响,在破裂传播的正前方,周期T=2.0 s,3.0 s,5.0 s,7.5 s和10.0 s的加速度谱被放大到整体观测平均水平的1 4～2.0倍.从周期T=2.0 s到10.0 s,破裂向前方向的放大作用和破裂反方向的减弱作用均有所增强,此次地震观测到的速度大脉冲记录均位于断层的东北方位,这与方向性脉冲的产生机理相吻合,速度大脉冲对加速度反应谱有显著的长周期放大作用,放大倍数值可以超过4.0,放大作用的影响主要位于脉冲的特征周期T_p附近.(3)近断层记录在建筑结构敏感的周期(0.5～2.0 s)的反应谱达到芦山地震的3～6倍,虽然与芦山地震震级接近,此次地震近断层地震动破坏力大大超过了芦山M_W6.8地震,甚至超过了汶川W_W7.9地震,这种长周期特点应该引起工程抗震设计和相关研究人员的重视.     

2016年日本熊本MW7.0地震强震动记录特征分析

对2016年日本熊本MW7. 0地震K-net和Kik-net强震动台网获取的强震动记录进行校正处理,分析了峰值加速度与峰值速度幅值特征,将观测结果与日本常用的3组衰减关系进行了对比,最后分析了强震动记录的持时特征和频谱特征。结果表明:①本次地震EW,NS及UD向最大峰值加速度分别为1 156. 95 gal,-786. 60 gal,874. 40 gal,与日本地区峰值加速度和峰值速度衰减公式预测值比较吻合;②5%-75%和5%-95%的2种能量持时整体上符合随距离增大而增大的规律;③KMMH16台站3个方向的加速度反应谱在0. 5～1. 2 s都超过1995年阪神地震中JR Takaroti台站的反应谱,本次地震对木结构房屋和低矮及中层钢筋混凝土房屋破坏严重。     

2021年日本福岛海域MW7.1地震S-net海底地震动特征

2021年2月13日日本福岛县东部海域(Off-Fukushima)发生的M_W7.1地震是目前世界最大的海底观测系统——日本海沟地震和海啸海底观测网(S-net)——迄今为止记录到最大的俯冲带地震.为认识海底地震动并探究海陆差异,本文介绍了S-net海域台网及其数据处理的特殊性和必要性,分析了Off-Fukushima地震的地震动幅值、频谱和持时特征,比较了海底与陆地震动衰减特性的异同,对比了近海岸埋置台站、近海沟未埋台站及陆地台站地震动的差异,并与俯冲带地震动模型(GMM)进行了比较.结果表明：Off-Fukushima地震海底水平向PGV和周期大于0.5 s的加速度反应谱值显著大于陆地记录;相比陆地,海底动力放大系数谱向长周期方向偏移,表明海底记录长周期成分更丰富;受到海水的抑制作用,海底记录V/H反应谱比值小于陆地记录;由于海底记录与陆地GMM在路径和场地方面存在差异,水平向PGA及反应谱的残差分布出现明显偏移;海底近海岸埋置台站与近海沟未埋台站的地震动特征存在系统性差异,水平向差异很可能与海底地形有关,而竖向可能受到布设方式对台站与海底耦合程度的影响.本文...     

九寨沟MS 7.0地震强震动特征及构造影响

从余震序列分布、震源机制、断裂构造、地震烈度和发震构造,简要论述九寨沟M_S7.0地震及构造背景,在此基础上,整理此次地震198条强震动记录,分析强震动峰值加速度分布特征（衰减规律）及构造影响因素,发现峰值加速度在发震断层两侧的衰减具有明显的不对称性,其中NE向衰减缓慢,SW向衰减迅速,该现象不仅与发震断层的倾向、倾角有关,还与不同构造块体的介质差异密切相关。     

云南景谷MS6.6地震地震动特征分析

2014年8月3日云南省鲁甸发生MS6.5地震后,10月7日云南省内再次发生景谷MS6.6地震。对我国数字强震台网捕获的39组主震记录进行处理后,分析幅值参数、Arias强度以及持时参数的衰减特征,发现两次地震触发台站远场记录的地震动参数均显著低于预测方程,鲁甸地震的近场地震动参数略高于景谷地震;而鲁甸地震地震动参数衰减快于景谷地震,这可能与云南东部地区Q值低于西部地区的区域差异性特征相关。分析结果可以为研究该区域6.5级左右地震动衰减特征提供参考。     

基于强震动记录测定近场地震参数方法研究

本文基于大理州的23个强震动台2007年以来记录的11组3.5级以上近场地震共51条强震记录,根据数字地震仪台网系统测定地震参数的基本原理,讨论了通过强震动台网测定近场地震参数的方法,并对测定近场地震参数的计算机软件MDEPS系统的编制进行了介绍。     

四川九寨沟MS7.0地震的地震动分布和衰减特征初步分析

2017年8月8日在四川阿坝州九寨沟县发生了M_S 7.0地震,这是继2008年四川汶川M_S 8.0地震和四川芦山M_S 7.0地震后在此地区发生的又一强震。本研究对此次地震中的66组主震记录进行了常规校正处理,基于此分析了此次地震中典型台站——九寨百合台站的地震动特征,并通过地震动的峰值加速度、反应谱参数分析了近场效应特征。研究结果表明:九寨百合台站的三分量峰值加速度都超过了100 gal,而地震动持时相对较短,并且三分量反应谱在短周期大于Ⅷ度多遇地震设计谱,小于Ⅷ度罕遇地震设计谱;由于九寨沟地震属于典型的走滑型地震,因此地震动参数的方向性特征不明显;九寨沟地震的地震动的峰值加速度以及反应谱谱值小于同等震级的芦山地震,并且小于几个典型衰减关系的预测值。对九寨沟地震的地震动参数特征研究表明,此次地震是一次震级较大、但灾害损毁程度较低的地震。     

山地与平原地形对地震动影响的对比分析

山地和平原地形在破坏性地震中对震害和地震动有很大影响,不同地形对地震动有不同的地形效应。为研究山地和平原地形对地震动的影响规律,基于集集地震后续M L6.6余震的强震观测记录,对位于不同地形且与震中位于同一直线上的四个强震台记录处理计算,从波形特征、频谱分析、加速度反应谱计算、与规范谱比较等方面进行对比分析,完成长周期地震动反应谱分布的计算。研究表明,山地地形对地震动的高频成分有地形放大效应,沉积平原对地震动的低频长周期部分有放大作用,并且长周期反应谱的峰值区域始终位于沉积平原中部。平原上强震记录产生的长周期放大系数谱已超出规范谱,沉积平原上长周期建筑的抗震设防应引起重视。     

利用近场强震记录反演2016年日本熊本县地震震源破裂过程

2016年4月15日16时25分（UTC）,日本熊本县发生M_W7.1强烈地震,给当地人员、建筑及经济造成严重灾难和巨大损失.日本地震观测网F-net给出的震源机制解显示此次地震的震源位置为130.7630°E,32.7545°N,深度12.45 km,节面Ⅰ：走向N131°E、倾角53°、滑动角-7°;节面Ⅱ：走向N226°E、倾角84°、滑动角-142°.与此同时,余震的震中分布及其震源机制结果显示主震的震源机制在破裂过程中有可能发生了变化,单一的震源机制不足以充分解释观测数据.本文依据GNSS和InSAR地表形变反演结果为约束,并结合活动构造资料为参考,构建了震源机制变化的有限断层模型,采用水平层状介质模型,利用日本强震观测台网K-NET和KiK-net的近场加速度观测记录,通过多时间窗线性波形反演方法反演了此次地震的震源破裂过程.研究结果显示,这是一次沿Futagawa-Hinagu断层带发生的右旋走滑破裂事件,发震断层分为南北两段,其中北段走向N235°E、倾角60°,南段走向N205°E、倾角72°,断层深度范围和余震深度分布基本一致,断层面上滑动主要集中于断层北段,最大滑动量约7.9 m,整个断层的破裂过程持续约18 s,释放地震矩5.47×10¹⁹ N·m（M_W7.1）.

     

基于汶川地震远场强震动记录的厚覆盖土层对长周期地震动影响分析

利用2008年5月12日汶川8．0级地震的远场强震动记录,研究太原盆地厚覆盖土层对长周期地震动的影响,并将其结果与地震安全性评价工作中的土层反应分析结果进行比较。结果表明：厚覆盖土层对长周期地震动有很强的放大作用,覆盖层厚度越厚,峰值出现的周期越大,而相对于基岩最大放大可达9倍左右。另外,对于厚覆盖层的场地,观测结果与场地地震安全性评价计算结果在一定的长周期范围内差别相当大,相对于基岩,观测结果的放大系数至少高于场地地震安全性评价计算结果的2倍以上,最大达到近9倍。     

日本熊本MW7.0地震同震位移场和震源滑动模型反演

利用2016年4月16日日本熊本M_W7.0地震震中周围94个近场强震动台的观测资料和新近改进的强震经验基线校正方法SMBLOC,尝试解算并绘制了一个内陆M7左右走滑型地震的同震位移场全貌,并反演了其震源滑动模型.与日本国土地理院(GSI)公布的该地震57个GPS同震位移结果的比较显示,两种完全不同资料、不同解算方法给出的水平同震位移场的最大幅值均为100 cm左右,均呈右旋走滑为主兼具部分正断分量的震源机制.强震最大水平和垂直永久位移分别为104.5 cm和58.0 cm, 分别出现在震中东北侧的KMMH162台和KMM005台.两种资料单独以及联合反演的震源滑动模型均表明,此次地震为北东侧破裂为主并呈双事件特征,且主要滑动均不在初始破裂点附近, 而是集中于第二次事件周围,即距离初始破裂点东北侧约20 km处的走向长约40 km、倾向宽约20 km的范围内.基于强震和GPS模型所得的最大滑动量分别为5.10 m和5.87 m,量级一致,反演矩震级均为M_W7.1左右;主破裂区近地表滑动量比野外调查结果略微偏大,可能与数值效应有关.此外,还利用不同方法得到的解算结果比较了熊本地震特有的12组台间距在3 km以内的GPS-强震台站对各自的三分量同震位移,其结果表明对于M7左右的地震而言,SMBLOC方法解算同震位移时方向和幅值的可靠性下限约为2 cm.      

2010年新西兰MW7.0主震与2011年MW6.1余震近场强地面运动特征比较分析

2010年9月3日16时35分46秒新西兰南岛Greendale附近发生了M_W7.0地震, 震源深度约10.0 km. 2011年2月21日新西兰南岛又发生了M_W6.1地震, 为2010年M_W7.0主震后最大的一次余震, 震源深度约5.0 km, 发震断层为Christchurch南约9 km一条近东西走向逆冲的隐伏断层, 该地震造成Christchurch城内多处建筑物严重损毁. 本文分析了2010年新西兰地震事件M_W7.0主震与M_W6.1余震强地面运动的特征. 新西兰M_W6.1余震近场强地面运动整体高于M_W7.0主震. 将主震和余震的强震观测记录分别与新一代衰减关系(NGA)进行对比, 发现余震强震观测数据整体高于其震级对应的NGA. 分别选取距离主震和余震震中最近且强震观测记录最高的两个台站(GDLC台站和HVSC台站)作为参照台站, 建立动态复合震源模型(DCSM)及有限断层随机振动模型(SFFM)进行强地面运动的模拟计算, 分析两种模型的模拟结果并对比二者的优势及局限, 以便在未来工作中更好地通过模型计算强地面运动特征, 实现区域化特征快速、 实时分析及局部重点、 细致分析相结合的目标.      

综合InSAR和应变张量估计2016年MW7.0熊本地震同震三维形变场

利用合成孔径雷达(Synthetic Aperture Radar,SAR)影像提取与地质活动相关的三维地表形变场,对深入理解地质灾害的形成机制及其潜在灾害风险评估非常重要.目前,利用SAR影像的同震三维形变场提取主要利用单个像素点的多次观测构建观测方程,然后基于加权最小二乘(Weighted Least Squares,WLS)方法分解从而获得同震三维形变场,因此该方法缺乏对相邻像素点空间相关性的约束.考虑相邻同震位移点的应力连续性,研究学者提出了顾及大地测量应变张量和卫星形变观测的SAR同震三维形变场方法(Extended Simultaneous and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measurements,ESISTEM).本文以2016年MW7.0熊本地震为例,收集了覆盖此次地震的ALOS-2卫星升降轨影像,利用传统差分InSAR(DInSAR)方法和子孔径雷达干涉测量(Multiple Aperture InSAR,MAI)方法分别对升降轨SAR影像对进行处理,得到视线向(LOS)形变和方位向形变,最后利用ESISTEM方法获取此次地震的三维同震形变场.此外,利用GPS和野外考察观测对本文的三维形变场结果进行结果精度分析.研究结果表明,与传统WLS方法相比,ESISTEM方法不仅能有效抑制奇异像素点对形变结果的干扰,同时对近断层的失相干信号能进行较好的恢复,更有助于解释地表破裂区的地震形变特征和掌握地震发生机制.本文确定的三维同震形变场结果显示主形变区发生在Futagawa断层中部和Hinagu断层最北端,最大水平位移为2m,抬升为0.55m.断层破裂以NE-SW走向的右旋走滑为主兼有部分正断成分.应变张量分析表明发震断层处受到了明显的收缩力和剪切力的作用.     

芦山7.0级地震近断层地震动的方向性

基于芦山7.0级地震中断层距小于100 km自由场台站的强震动记录观测数据,研究此次地震近断层地震动的方向性特性,并探讨方向性特性与震源破裂机制、断层距离和空间方位的关系.研究结果表明：（1）与距断层较远记录不同,近断层地震动在不同的观测方向上表现出显著的强度差异,存在明显的极大和极小作用方向.在不同的方向上,最大加速度反应可达最小值的4倍以上;（2）这种方向性差异在T=1.0 s以上的长周期段更为明显,在T=0.1 s以下的短周期段,地震动随方向变化的差异较小.地震动强度随方向变化的差异随周期增大而增大,不同方向上的加速度反应谱值的最大值与最小值之比从周期T=0.01 s时的约1.7增大到周期T=10 s时的约2.4;（3）在距离断层约35 km以内,地震动具有明显方向性,地震动卓越方向具有垂直断层走向的特征,随断层距的增大,这种方向性不明显.从不同方向上地震动强度的差异来看,随断层距增大,地震动强度在不同方向上的差异在减小,表现为各个周期的最大值/中值和最大值/最小值比值均随断层距离增大缓慢减小;（4）近断层地震动的方向性特性主要受断层上、下盘的相对运动所控制,其在长周期的卓越方向与水平同震位移方向一致,且该卓越方向上的地震动强度绝对大小与地震破裂造成的静态位移明显相关,表现为地震动强度随水平同震位移的增大而增大.

     

九寨沟MS7.0地震强地震动模拟及漳扎镇地震动强度预测

破坏性地震强度预测可用于工程领域抗震设防以及地震危险性分析评估,是防震减灾中一项很重要的基础工作.为了再现九寨沟地震的地震动强度,评估缺失强震记录的九寨章扎台站的地震动强度,本文用经验格林函数法对九寨沟地震进行了数值模拟.选取了震源周边地震动峰值加速度超过10 Gal的10个强震台站进行模拟.因未得到九寨沟地震的余震,初次尝试将汶川地震和定西地震的余震作为格林函数模拟九寨沟地震.模拟结果整体上可以反映各台站地震动的强度特征,尤其是地震动高频成份拟合较好.模拟值的地震动峰值加速度、时程数据、反应谱等与观测值拟合较好.预测结果显示漳扎镇的地震动峰值加速度值约为180~200 Gal.预测结果也表明在缺少大震的余震记录时,经验格林函数法使用其他大震的余震同样可以再现目标地震的强度特征.本研究也为经验格林函数方法在缺乏小震记录地区的使用积累了经验.最后总结了格林函数的选取标准,为经验格林函数方法来预测未来强震动特征积累了经验.

     

智利M_W8.3地震与M_W8.8地震的震源过程及其引起的库仑应力特征

2015年9月17日6时54分32秒(北京时间)智利中部伊拉佩尔附近(震中31.57°S,71.67°W)发生了一次M_w8.3大地震,在此次地震震中以南约500 km处的马乌莱地区曾于2010年2月27日14时34分11秒发生过一次M_w8.8强震(震中36.12°S,72.90°W),两次地震余震分布区之间有约75 km的地震空区.本文利用远场体波与面波波形,基于有限断层模型,反演了这两次地震的震源破裂过程.结果显示这两次地震均为逆冲型大地震,2015年伊拉佩尔M_w8.3地震的平均滑动角度为107°,平均滑动量为2.43 m,平均破裂速度为1.82 km·s~(-1),标量地震矩为3.28×10~(21)Nm,95%的标量地震矩在104 s内得到了释放.最大滑动量约8 m,位于沿走向75 km,深度8 km处.2010年马乌莱M_w8.8地震的平均滑动角度为109°,平均滑动量为4.95 m,平均破裂速度1.90 km·s~(-1),标量地震矩为1.86×10~(22)Nm,95%的标量地震矩在121 s内得到了释放.最大滑动量约12.5 m,位于沿走向100 km,深度21 km处.2015年伊拉佩尔M_w8.3地震浅部更大的滑动量应该是其引起了较大海啸的一个原因.基于破裂滑动分布,我们计算了这两次地震引起的周边俯冲带上静态库仑应力变化,结果显示两次地震均显著增加了周边俯冲带上的库仑应力,2010年马乌莱地震使得2015.年伊拉佩尔地震震源区附近的库仑应力增加了(0.01~0.15)×10~5Pa,从应力积累的角度看,2010年马乌莱地震有利于2015年伊拉佩尔地震的发生,对后者的发生起到了促进作用.     

2018年2月6日花莲M_W6.4地震近场地震动方向性效应

利用中国台湾省内222个强震动台站以及Palert地震预警系统520个台站所观测的三分量加速度记录,研究此次花莲M_W6.4地震近场强地震动空间分布和衰减特征,将观测结果与美国NGA-West2地震动经验预测模型进行对比,揭示此次台湾花莲地震近场地震动的长周期特点,基于回归残差分析研究地震动峰值加速度(PGA)、峰值速度(PGV)和不同周期地震动的空间分布差异,定量考察近场地震动的方向性效应.研究结果表明:(1)整体上此次地震的近场PGV观测值和周期1.0s以上的长周期加速度谱值与美国NGA-West2地震动预测模型结果接近,PGA观测值和周期小于1.0s的加速度反应谱略低于预测模型结果.从空间分布来看,周期1.0s以上的长周期地震动在断层的不同方位有系统性差异,在破裂传播前方(震中西南方位),周期大于1.0s时的反应谱明显高于美国NGA-West2地震动经验预测模型,在破裂传播后方(震中东北方位),周期大于1.0s时的反应谱低于经验预测模型,表明此次地震近场地震动具有显著的方向性效应.(2)破裂传播的方向性效应主要影响周期超过1.0s的长周期,而对PGA以及周期小于1.0s的短周期地震动影响较弱.在破裂传播前方,周期1.0～10.0s的加速度反应谱值被增强到整体观测平均水平的1.16～1.52倍;在破裂传播后方,周期1.0～10.0s的加速度反应谱值被减弱到整体观测平均水平的0.36～0.70倍.(3)此次地震破裂方向性效应的影响表现出明显的窄带效应,破裂方向性的影响(包括破裂传播前方的增强作用和破裂传播后方的减弱作用)在周期T=3.0s时达到最大,在该周期破裂传播前方的增强系数为1.52,破裂传播后方的减弱系数为0.36.从周期T=3.0s到10.0s,破裂方向性效应的影响随周期增大总体上呈减弱趋势,这与2016年日本熊本M_W7.0地震破裂方向性效应的影响特点显著不同.     

2016年熊本MJ7.3地震的工程地震动参数模拟及分布特征分析

本文选取2016年4月16日日本熊本县M_J7.3 （M_W7.0）地震近场区域内K-net地震台网的47个强震台站所记录的加速度数据,运用经验格林函数法模拟分析此次地震主要的工程地震参数并给出了各工程参数的空间分布。通过对比分析得到结论如下:① 地震动时程的基本频谱的模拟结果较好,尤其是1—15 Hz的高频段内;② 在震源距小于50 km范围内,峰值加速度和阿里亚斯强度的观测值与模拟值拟合较好,二者在近场区域均以椭圆形向周围扩散衰减,且阿里亚斯强度的模拟值整体大于观测值;③ 卓越周期拟合整体较好,但在局部区域存在较大差异,模拟结果难以表征场地环境的复杂性。      

2008年青海大柴旦Mw6.3地震地面运动模拟

2008年11月10日在青海柴达木盆地北缘发生了大柴旦M_W6.3地震,为了研究该地震的区域地震波传播与地面运动特征,本文利用地质资料和地壳速度结构研究成果,构建了柴达木盆地及周边区域三维传播介质模型,采用有限差分方法模拟了大柴旦地震波场传播过程以及地面运动分布特征.结果表明,柴达木盆地对波场传播有明显影响,表现为地震波传入盆地后在边界产生次生面波,盆地沉积物对地震波具有围陷作用,地震地面运动在盆地内振幅增大、持时延长.模拟结果给出的地震地面运动峰值速度分布以及理论地震图均和观测结果符合较好,反映数值模拟较好地给出了观测地面运动的主要特征以及传播介质模型的合理性.