日本1923年关东地震（M7.9）近场加速度合成

本文用统计一经验格林函数法估计了日本1923年关东地震(M7.9)的近场加速度时程,合成的加速度峰值、持时与他人的研究结果符合较好,说明统计-经验格林函数法是在没有小震记录场点估计近场地震动的一个有效方法。合成结果也说明1971年Kanamori推定的关东地震断层模型2适合用来估计该次地震的近场地震动。根据九个场地的合成结果推测,关东地震近场基岩加速度最大值约为300Gal,地表加速度最大值约为600Gal。     

用经验格林函数方法研究澜沧-耿马MS=7.6地震的破裂过程x,auto,auto,415px);} h

运用经验格林函数反褶积方法,用中国数字地震台网(CDSN)6个台记录的MS=6.3余震的地震图作为经验格林函数,研究了1988年11月6日中国云南澜沧-耿马MS=7.6地震的破裂过程．反褶积结果表明,MS=7.6地震是一次比较简单的事件．分析了破裂的方向性,得出MS=7.6地震的破裂开始于微观震中,然后沿西北、东南两个方向近似对称地扩展,破裂总长度为70 km,持续时间19 s,破裂速度约为2 km/s．      

用经验格林函数方法从长周期数字波形资料中提取共和地震的震源时间函数

1990年4月26日，在青海省共和县发生了一次震级MS=6.9地震．这次地震之后，于1990年5月7日、1994年1月3日和1994年2月16日在原震区又发生了MS分别为5.5，6.0和5.7的余震．我们分别以这几次余震的中国数字地震台网的长周期记录作为经验格林函数对该地震的长周期记录进行反褶积，提取了该地震的远场震源时间函数．无论以哪次余震的记录作为经验格林函数，反褶积的结果都相当一致．结果表明，发生于1990年4月26日的MS=6.9地震至少由两次规模相当的事件构成，两次事件的发震时刻相差约30 s．第一次事件的持续时间较短，约12 s，震源时间函数的上升时间约5 s；第二次事件的持续时间较长，约17 s，上升时间约8 s．分析分别从P波和SH波提取的震源时间函数，我们注意到，从SH波中提取的震源时间函数比从P波中提取的震源时间函数复杂，表明这次地震的震源过程除了上述两次规模相当的事件外，还存在规模较小的事件．这一结果与我们由矩张量反演得到的结果以及用宽频带波形资料借助于经验格林函数方法所得的结果一致．我们以较小余震对较大余震做反褶积得到较大余震的震源时间函数．结果表明，发生在1994年1月3日的MS=6.0和1994年2月16日的MS=5.7地震的震源过程相当简单，其震源时间函数为一简单脉冲．对于1990年4月26日MS=6.9地震，我们从不同台站的资料中得到的震源时间函数具有方位依赖性．然而，这种方位依赖性在3个余震的震源时间函数中几乎没有显示．从不同台站或同一台站的不同震相求得的震源时间函数的幅度存在差异，但震源时间函数在其持续时间上的积分却相当稳定，亦即从任何一个台站测得的标量地震矩相当一致． 以3个余震的记录作为经验格林函数得到的1990年4月26日地震的相对标量地震矩分别为22(以1994年1月3日MS=6.0地震的记录为经验格林函数)、26(以1994年2月16日MS=5.7地震的记录为经验格林函数)和66(以1990年5月7日MS=5.5地震的记录为经验格林函数)．据此推算，MS=6.0地震与MS=5.7地震的相对标量地震矩为1.18；MS=6.0地震与MS=5.5地震的相对标量地震矩为3.00；MS=5.7地震与MS=5.5地震的相对标量地震矩为2.54．这些结果与我们从资料中直接得到的相对标量地震矩(分别为1.15，3.43和3.05)一致，表明了作为提取震源时间函数的经验格林函数方法副产品的相对标量地震矩是地震大小的一种很稳定的度量.      

基于改进的经验格林函数方法的地震动模拟

改进的经验格林函数方法弃掉以往完全用经验方法确定用于合成大震的小震数目的方法,考虑了大地震断层各个子断层位错的不均匀性,使得经验格林函数方法的物理意义更加清楚。用改进的经验格林函数方法,以汶川MS6.4余震记录作为经验格林函数模拟了九寨白河台、松潘安宏台、皋兰台和成都中和4个台站的加速度时程,与以往的经验格林函数方法比较,模拟结果的反应谱与记录符合得更好,证明了改进的经验格林函数方法可以反映震源位错的复杂分布,可以在工程地震中用来估计场地地震动过程。     

考虑凹凸体理论的经验格林函数方法

以15次地震资料为基础,回归得到断层上所有凹凸体地震矩与凹凸体面积的关系;研究经验格林函数方法中凹凸体子源的应力降,提出基于断层凹凸体理论的经验格林函数方法;以唐山地震余震（Ms=6.9）作为经验格林函数,合成1976年唐山地震（Ms=7.8）的地震动,比较合成结果与观测记录表明,该方法符合断层破裂过程,计算结果可靠．      

汶川地震破裂过程联合反演及高频辐射研究

汶川地震是发生在叠瓦状曲面断层上的一次复杂破裂过程,北川—映秀断层和灌县—江油断层以及与其垂直相交的小鱼洞断层均有地表破裂发生,已有动力学和运动学研究表明断层破裂扩展至3个断层相交处发生了重要的转换。因而,北川断层、彭灌断层和小鱼洞断层的破裂顺序是反演汶川地震破裂过程的关键。本文首先建立合理的三维复杂断层模型,采用3种可能的破裂方式,基于并行非负最小二乘法和多时间视窗技术,联合远场、近场、GPS和地表破裂资料,反演汶川地震精细破裂过程,给出了合理的破裂方式。进一步采用近场加速度记录和差分进化方法,反演断层面上高频辐射分布规律,对比分析高低频地震波辐射的差异。主要成果如下:1.综合考虑三维发震构造模型、余震分布和地表破裂调查,建立更符合实际的曲面断层模型;采用远场36个台站垂直向P波位移记录,利用并行非负最小二乘法结合多时间视窗技术,反演了汶川地震破裂过程。研究表明:(1)目前震源破裂过程反演广泛采用的单侧破裂,在北川断层虹口—映秀近地表区域不产生与地表破裂相符的位错。如果彭灌断层与小鱼洞断层相交处发生双侧破裂,则会在断层南段近地表处产生高达4 m的位错,且在北川断层虹口—映秀近地表区域不产生位错,这与地表破裂矛盾。而北川断层浅部区域从与小鱼洞断层相交处发生双侧破裂,则会在虹口—映秀近地表产生与地表破裂相符的位错,明显优于单侧破裂。综合3种破裂方式的结果,本文认为北川断层在与小鱼洞断层相交处发生双侧破裂更符合汶川地震的实际情况。(2)北川断层南段和彭灌断层空间上接近,两者在远场台站的格林函数波形相似,基于远场记录的反演不能区分两者的位错,且北川断层南段位错分布的可靠性大于彭灌断层的位错分布。2.选用近场方位角覆盖较均匀的43个台站三分向速度记录,采用并行非负最小二乘算法结合多时间视窗技术,得到了汶川地震破裂过程。研究表明:(1)与远场记录反演结果相似,单侧破裂和彭灌断层在与小鱼洞断层相交处发生双侧破裂,北川断层虹口—映秀浅部都不产生位错。要使虹口—映秀区域发生与地表破裂吻合的位错,只有北川断层在与小鱼洞断层相交处发生双侧破裂才能满足,所以北川断层需发生双侧破裂。北川断层西南侧高倾角区域有明显的破裂停顿和二次破裂,部分区域破裂持时可达15 s左右。(2)51MZQ、51SFB、51MXN和51MXT等台站水平向较大的PGV与断层面相邻区域的滑动速率具有很好的一致性。说明断层面上发生较大滑动速率的区域,其临近台站往往伴随有较大的PGV产生。3.为了克服单一数据分辨率不足的缺陷,联合远场和近场资料以及联合远场、近场、GPS和同震位移观测资料反演了汶川地震破裂过程。研究表明:(1)联合远场和近场资料,提高了北川断层南段高倾角部分、PGF南半段区域以及北川断层北川附近滑动分布的识别能力。(2)反演中加入GPS资料能很好地控制断层浅部和北川断层北段的滑动分布。(3)联合反演结果显示,汶川地震破裂持续时间达100 s,释放地震矩为1.058×1021 N·m,断层面上存在5个凹凸体,表明此次地震至少由5个子事件组成。滑动主要分布在北川断层上,说明北川断层是主要的破裂面。在北川断层南段上,龙门山镇下侧以及虹口—映秀近地表区域的位错以逆冲错动为主,最大滑动量达12 m,位于虹口下侧;在岳家山到清平近地表附近错动以逆冲为主兼有走滑错动,最大滑动量约为10 m。北川断层北段上,北川附近滑动以逆冲为主,最大滑动量10 m;南坝到青川区域以走滑错动为主,最大滑动量10 m。在彭灌断层上,白鹿下方区域的位错也以逆冲为主,断层深部位错达8 m。4.利用芦山地震记录建立的加速度包络衰减关系和汶川地震近场30个台站的加速度包络,基于线源模型采用差分进化方法反演了汶川地震断层面上高频(1 Hz)辐射区域分布。结果表明:(1)断层面上高频辐射分布很不均匀,辐射较强的区域主要位于:产生较大地表破裂的映秀、北川和南坝区域;映秀和北川等凹凸体的周边区域,包括震中东北侧60—90 km区域、北川和南坝东北侧30 km处;断层破裂停止的东北端约30 km长的区域。其中,破裂贯穿到地表的映秀、北川和南坝是高频和低频辐射都很强的区域。(2)对于无观测记录场点,选择其临近且场地条件类似的台站加速度提取平稳随机过程,结合高频辐射分布和衰减关系得到的包络,合成了加速度时程,可为汶川地震结构震害分析提供地震动输入。     

估计震源破裂模型的经验格林函数法

用经验格林函数法模拟断层破裂过程，合成近场地面运动，并通过合成结果和真实记录的时程曲线、反应谱、傅氏谱的对比，估计断层破裂模型。     

2003年10月25日甘肃省民乐-山丹MS 6.1地震强震近场记录和分析

2003年10月25日甘肃民乐-山丹间发生Ms6．1地震．震前在兰州地震研究所明确的预报意见指导下就派出4台数字流动强震仪到甘肃省河西地区进行流动观测，其中民乐流动台获得近场主震记录．震中距为12．6km，EW向加速度峰值为33lcm／s^2，NS向加速度峰值为332cm／s^2，垂直向加速度峰值为209cm／s^2．其它固定自由场数字强震台山丹、肃南、武威也获得记录．本文分析了近场主震频谱以及远场加速度衰减规律．     

青海门源5.1级地震强震动记录及初步分析

对西北强震动台网中心收集的2013年9月20甘肃肃南与青海门源交界5.1级地震的44组三分向加速度记录进行了处理和初步分析,给出了126条零基线校正后的三分向峰值加速度值。距震中37 km的63 MEY台加速度峰值约为50 gal,可得出该次地震震中烈度约为Ⅴ度,与震后的实际调查结果基本一致;选取了近场典型台站的加速度记录进行分析,绘制出了校正后的加速度时程曲线以及加速度反应谱,反应谱幅度及卓越周期的差别体现了震中距对地震传播和地震在不同场地环境中的影响。     

2021年5月21日云南漾濞MS6.4地震随机有限断层三维地震动模拟

2021年5月21日云南省大理州漾濞县发生MS6.4 地震,震中附近遭受了强烈地震破坏.为预测此次地震的地震动影响场,利用震源运动学破裂随机模型,基于随机有限断层三维地震动模拟方法,给出了此次地震中 28个触发强震动台站的三分量加速度时程模拟记录,并结合强震动观测记录,估计了此次地震的地震应力降及震源破裂过程,进一步模拟给出了此次地震中 2823 个虚拟观测点的三分量加速度时程.结果表明,模拟记录的峰值地面加速度(PGA)、峰值地面速度(PGV)与观测值接近,并体现了地震动峰值的衰减规律、近场饱和效应和破裂方向性效应;模拟与观测记录的 5%阻尼比拟加速度反应谱(PSA)的幅值接近、谱形相似,在0.05～10 s周期段,模拟记录可以很好地预测地震动.基于三分量模拟记录给出了漾濞MS 6.4 地震的仪器测定地震烈度图,与云南省地震局发布的烈度图接近,极震区烈度最高可达Ⅷ度,震源破裂方向性导致震中 SE方向的烈度普遍高于 NW方向,受局部场地条件影响沿洱海西侧出现高烈度异常区.     

台湾费丛峡谷对强地面运动的影响研究──利用单台资料及2．5维SH混合方法

本文通过格林函数反褶积方法，由台湾峡谷附近的记录资料预测峡谷区的强地面运动的时程曲线．峡谷区的理论格林函数应用２．５维ＳＨ混合方法求解．通过格林函数反褶积方法得到的峡谷区费丛１及费丛２台的位移、速度、加速度和实际资料对比，取得了满意结果．计算结果表明，峡谷底部的峰值加速度相对峡谷边缘为最小；在靠近震中的一侧，峡谷的加速度的最大振幅相对比另一侧大．还给出了其余３个台的预测结果，研究了它们的加速度傅里叶谱和反应谱．     

台湾费丛峡谷对强地面运动的影响研究──利用单台资料及2．5维SH混合方法

本文通过格林函数反褶积方法，由台湾峡谷附近的记录资料预测峡谷区的强地面运动的时程曲线．峡谷区的理论格林函数应用２．５维ＳＨ混合方法求解．通过格林函数反褶积方法得到的峡谷区费丛１及费丛２台的位移、速度、加速度和实际资料对比，取得了满意结果．计算结果表明，峡谷底部的峰值加速度相对峡谷边缘为最小；在靠近震中的一侧，峡谷的加速度的最大振幅相对比另一侧大．还给出了其余３个台的预测结果，研究了它们的加速度傅里叶谱和反应谱．     

震源深度误差对矩张量反演的影响

在用矩张量反演方法确定地震的震源机制和震源的破裂过程时，由于测定震源深度的误差，很难保证计算格林函数时所采用的震源深度(理论震源深度)和实际震源深度完全相等．理论震源深度和实际震源深度的差异(震源深度误差)对矩张量的反演或多或少会造成影响．借助于合成地震图较系统地讨论了震源深度误差对3种基本类型(正断层、逆断层和走滑断层)的震源反演造成的影响．对于正断层和逆断层，震源深度误差主要影响各向同性成分和补偿线性向量偶极成分．在理论震源深度大于实际震源深度时:对于正断层，震源深度误差导致出现虚假的正的各向同性分量和负的补偿线性向量偶极分量；随着差异的增大其绝对值逐渐增大；对于逆断层，震源深度误差导致出现虚假的负的各向同性分量和正的补偿线性向量偶极分量，其绝对值随着震源深度误差的增加而增加．在理论震源深度小于实际震源深度时，结果恰好相反:对于正断层，震源深度误差产生了虚假的负的各向同性分量和正的补偿线性向量偶极分量；对于逆断层，震源深度误差产生了虚假的正的各向同性分量和负的补偿线性向量偶极分量．同样，随着深度误差的增大，它们的绝对值也逐渐增大．对于走滑断层，震源深度误差对矩张量反演的影响则不同于正断层和逆断层两种情况，受影响最大的是震源时间函数的形状．无论理论震源深度小于还是大于实际震源深度，震源时间函数的形状都发生不同程度的拖尾现象．数字试验表明，当震源深度误差小于20 km时对于地震的总体机制的反演没有明显影响．另外，当理论震源深度大于实际震源深度时，震源深度误差对矩张量反演的影响相对较小．      

基于强震和测震数据对鲁甸6.5级地震发震断层方向的研究

2014年8月3日16时30分,云南省昭通市鲁甸县发生M6.5地震,宏观震中位于鲁甸县龙头山镇(27.1°N,103.3°E),深度约12km。发震断层为NW向包谷垴-小河断裂,是NE向昭通-鲁甸断裂系的NW向次级走滑断层。文中利用鲁甸地震50个强震台站观测记录与135个测震台站观测记录,分别进行峰值加速度计算和S波辐射图案匹配,对此次地震发震断层方向进行讨论。结果表明,"NW向"的断层方向更符合强震、测震仪器观测值,同时与地震灾害现场调查结果一致。     

日本东京城区四个场地的土层放大效应分析(英文)

本文用统计一经验格林函数法合成１９２３年关东地震近场加速度时程，对比分析东京城区四个相近场地的土层放大效应，说明在结构抗震分析中考虑土层放大效应的重要性。     

日本东京城区四个场地的土层次放大效应分析

本用统计－经济格林函数法合成１９２３年关东地震近场加速度时程，对比分析东京城区四个相近场地的土层放大效应，说明在结构抗震分析中考虑土层放大效主尖的重要性。     

云南漾濞6.4级地震强地面运动的模拟和空间分布特征分析

北京时间 2021 年05月21日 21 时48分云南大理州漾濞县(北纬 25.67°,东经 99.87°)发生了 6.4 级地震.本文利用震源破裂过程控制的 NNSIM随机有限断层方法,模拟对比了近场 6 个强震台的加速度、速度、反应谱记录,据此确定了开展近场地震动模拟所需要的参数的大小,进而建立了近断层 200 km范围内的地面运动.分析了漾濞地震近场地面运动的空间分布特征,结果显示此次地震地面运动的峰值速度PGV、峰值加速度PGA以及反应谱分布均表现为圆形分布,无明显的上下盘特征和走向特征,震中极震区峰值加速度超过了 400 cm·s－2 ,对应国家标准烈度Ⅸ度.     

1997年中国西藏玛尼Mem>Ssub>7.9地震的时空破裂过程l

用中国数字地震台网（CDSN）的长周期波形资料反演了1997年11月8日中国西藏玛尼地区MS7.9地震的地震矩张量,用频率域里反褶积方法从P波和S波震相中分别提取了震源时间函数,并经反演依赖于方位的震源时间函数获取了断层面上破裂随时空变化的图象．矩张量反演结果表明:玛尼地震发震应力场的P轴和T轴均接近于水平,P轴在NNE方向（方位角29,倾角7）,T轴在SEE方向(方位角122,倾角23),断层错动以走滑为主;标量地震矩为3.41020 Nm,矩震级MW=7.6．由矩张量反演得到的震源时间函数显示,这次地震是由一次较小事件和较大事件组成的,较小事件大约持续5 s,较大事件持续约10 s．由余震分布可推断出玛尼地震的发震断层是走向为250、以走滑为主的左旋-逆断层,断层面的倾角比较陡,约88．根据反演结果计算了理论格林函数,然后用反褶积方法提取了震源时间函数．从不同台站的P波和S波中分别提取的震源时间函数一致表明这次地震破裂的时间历史比较简单,可用一宽度约10 s的正弦形的函数近似表示．进一步反演从不同台站上得到的、依赖于方位的P波和S波震源时间函数,获得了断层面上滑动的时空分布图象．从破裂的记忆式快照看,破裂开始于断层的西端,然后向东向下发展,总体上具有单侧破裂的特征．破裂面由3个破裂子区构成．一个在断层西端,深度约10 km(西区);另一个距断层西端约55 km,深度约35 km(东区);第3个距断层西端约30 km,深度约40 km(中区)．３个破裂子区构成约长70 km,宽60 km的破裂面．从破裂的遗忘式快照看,这次地震的破裂过程是相当复杂的,在不同时刻断层面上发生错动的地点并不相同,显示出这次地震的破裂过程具有愈合脉冲的特征,而且在断层面上的某些部位发生了多次错动;另一特征是最先和最后破裂的部位都不是主要的破裂区．根据标量地震矩计算了断层面上静态位错的分布,位错最高的3处（西区、东区和中区）的位错值分别为956 cm,743 cm和1 060 cm．由断层面上位错的分布推知,破裂主要集中在震中以东长约70 km的断层上;从余震的分布看,震中以西余震稀疏而震中以东余震密集．这些都表明这次玛尼MS7.9地震是北东东-南西西向至近东-西向断层向东扩展的结果．      

1976年唐山地震震源动力过程的数值模拟

用新LDDA(Lagrangian DiscontinuousDeformation Analysis)方法模拟了唐山地震断层的破裂、错动和应力释放的整个动力过程．模拟结果表明,唐山地震的震源滑动过程在发震断层上各处不一样．近场位移受断层的曲率影响,断层凹侧的位移大于断层凸侧的位移．滑动过冲现象在震中处最大,并向断层两端衰减．我们发现,唐山地震断层的破裂速度和应力降与断层上的初始剪应力大小有关．唐山发震断层的最大动态、准静态位错量和剪应力降均发生在中间部位,分别是7.1 m、6.2 m和8.1 MPa、5.4 MPa,发震断层的平均准静态位错量和剪应力降分别为4.5 m和3.3 MPa,断层破裂的传播速度从震中向东南和西北方向分别为3.08 km/s和1.18 km/s．      

华北地区地震危险性分析和地面运动预测的一致性方法

利用华北地区地震活动性资料,建立了地震危险性计算的一致性模型.在此模型的基础上,得出了北京、天津、唐山和济南等7个城市未来2500年内地震的时空强度分布,并计算了2500年回复周期的地震动峰值加速度(PGA).结果表明,唐山和太原的PGA最大(>0.2g),石家庄和北京次之(≈0.17g).对华北地区2500年地震记录的正演计算结果表明,太原和唐山地区的潜在地震危险最有可能来源于震级在6.0～7.0、震中距离在12～15km的地震活动;而北京、天津和石家庄地区则可能来源于震级在5.5～6.0、震中距离在10km左右的地震活动.采用IBC(International Building Code)方法计算后的结果显示,太原、唐山等地区的PGA与2001年我国地震动峰值加速度值基本一致,与此地区的较高地震活动性特征相符.利用随机震源模型,还给出了影响此7个城市的最大地震记录的加速度、速度及位移时程曲线,这对本区工程建筑的抗震性设计以及对救援设施的选址等有重要作用.