最大可能地震相关参数的统计估计

提出了一种最大可能震级Ｍｍａｘ及一些相关参数估计的新型统计方法,推导了Ｍｍａｘ无偏估计公式及该参数的任意函数。这种估计已有明确的公式,且方差很小。亦推导了这种预示性估计的极限分布,并将所得公式通过几个例子进行了说明．     

用不完整的数据文件估计地震危险性参数 第一部分：具有不同临界震级的极值

地震危险性参数（最大区域震级mmax、地震活动速率λ和古登堡－里克特方程式中的参数b）的最大似然估算方法已被扩展用于处理包括历史大震和近代完整观测资料在内的混合资料。该方法可以收取利用具有不同临界震级值的地震目录中各种质量不一的完整资料。本程序作为一个案列证被用来估计卡拉布里亚和西西里东部地区的地震活动性参数。     

地震震级—频度关系b值参数的稳健估计

文中依据古登堡－里克特震级－频度公式，结合实例说明了参数ｂ值的稳健估计方法，并将其与最小二乘法及最大似然法比较，结果表明，一般情况下，稳健估算结果与实际数据分布较吻合，并且可以在不剔除低震级段和大震级的偏离分布值的情况下，取得与选取合理震级上下限最小二乘法一致的结果，此外，稳健估计结果的余差分析同样也得到地震震级－频度关系不是简单的对数线性关系的结果，而原始Ｇ－Ｒ公式需加二以上的高次项才能较好地押     

完整资料的地震危险性参数估计及敏感性分析

应用考虑震级不确定笥和资料完整性估计地震危险性参数的极大似然法，分别采用１９４０年和１０７０年至１９９３年完整的地震资料，计算了汾渭地震带及陕南地区的地震危险性参数。结果表明：利用近几十年来的完整地震资料结合研究区历史上发生的最大地震，可以得到较为可靠的地震危险性参数；参数敏感性分析表明，完整震级不确定性对β的影响较小；随震级不确定性的增大地震年发生率减小；在资料完整的前提下，震级下限在一定范围内     

哈萨克斯坦境内地震活动的主要特点

哈萨克斯坦地震活动带从西往东延伸长约２０００ｋｍ，宽１５０ｋｍ。在这个区域内曾发生过世界闻名的灾难性地震。外伊力阿拉套，昆格阿拉套和准噶尔阿拉套是这里地震最活跃的地区。在这里进行综合地震研究可以更具体、更准确地说明：哈萨克斯坦境内地震活动的时空分布特征及其详细部位；主要地震活动带长期的地震动态参数平均值；根据不同参数进行境内的地震区划；估计最大能级Ｋｍａｘ（Ｍｍａｘ）值并以此为依据论证和评估哈萨克斯坦境内的地震危险性。     

地震危险性分析中的参数敏感性研究

强调地震危险性分析是包含多参数的复杂系统，其参数的敏感性应作为一个整体来研究．考察了参数连续变化对结果的影响以及超越概率和潜在震源区震级上限对参数敏感性的影响，得到了一些有益的认识．全面计算了包括地震空间分布函数在内的参数之间的交互作用，提出了相对交互作用的概念．得出在50年超越概率10%时相对交互作用不超过10%；在年超越概率10－4时，高震级潜在震源区的震级上限与年平均发生率、空间分布函数的相对交互作用分别可达到17%和22%，为实际应用时判断是否要考虑参数的交互作用提供了方便．指出了求解结果的概率分布的重要性和必要性，分析了各参数取不同水平的权重对结果分布的影响，揭示多参数敏感性的全貌，为进行合理的不确定性校正作了必要的准备．      

丹江口水库及邻区地震震源参数研究

在讨论震源参数和Ｍ震级关系以及论述利用地方震尾波估算地震矩方法的基础上，尝试性地系统估算了丹江口水库及邻区内４３次ＭＬ≥２．５级地震的震源参数，并讨论了丹江口水库地震的特征。丹江口水库及邻区地震震源参数的初步分析表明：丹江口及邻区属低应力降地区。水库蓄水的附加引张应力效应可能使得丹江口水库中局部以引张应变为主的地区应变迅速增加，地震应力降高于该区一般地震原有水平、从而诱发地震。此外，丹江口水库地震的破裂长度可能和震源机制的破裂类型有关。     

西安及周围地区小地震的震源参数

利用中日合作微震观测台网１９８８年３月至１９９１年３月间的地震记录计算了西安及周围地区小地震的波谱和震源参数,并探讨了震源参数与震级ＭL之间的关系.     

多随机变量回归及其参数的误差估计

本文从理论上比较严格地推导了多随机变量回归参数及其误差估计的一般表达式,从而从理论上完善了多随机变量回归方法。本文推导得到的多随机变量回归参数及其误差估计表达式具有计算步骤简捷、变量之间的函数关系可以任意的特点。最后,文中还利用导得的公式对断层长度与震级关系以及震级与震中烈度关系中的参数及其误差进行了具体的计算。     

用于地震作用预测的弱震数字记录分析

本文中使用了麦迪奥台９０个小震的近震数字记录，这些小震分布在阿拉木图市附近，处在东经７６°４０’～７７°４０’和北纬４２°５０’～４３°２０’座标所局限的区域范围内。分析了ＣＫＭ信道水平分量Ｓ波振相的记录。对于所有被研究的数字地震记录，得到了如下动力学参数：基底最大位移的振幅Ａｍａｘ，与Ａｍａｘ相应的周期Ｔ和０．５Ａｍａｘ水平上的相对持续时间ｄ。这些参数是对于各个分量按每个滤波信道来测定的。编制了原始地震记录图以及最大振动分量的参数Ｔ和ｄ的整体分布。确定了７０％置信区间内的中值和误差。它们的分布接近于正态分布。对Ｔ和ｄ的研究表明；这些参数彼此互相关：周期愈大，相对持续时间的值也愈大。得到了Ｔ和ｄ与地震能级及震级的初步关系。     

余震对地震危险性评估的影响：以北安纳托利亚断裂带为例

为了调查研究余震对地震危险性评估的影响,根据主震目录和原始地震目录用最大似然方法估计了北安纳托利亚断裂带（ＮＡＦＺ）的地震危险性参数（λｍ、ｂ和Ｍｍａｘ）。主震目录是使用窗口方法朋原始地震目录中去除余震编制的,原始地震目录是由１９００￣、９９２年仪器检测地震和１０００￣１９００年发生的历史地质组成的。对于主震目录的事件而言,泊松过程是适用的,而对于原始地震目录则不适用。本文论证了使用主震目录和原始     

中国大陆活动地块边界带最大震级分布特征研究

中国大陆活动地块边界带是强震活动的集中带,也是地震预测研究与防灾减灾的主要目标区,而最大地震震级是描述区域地震活动性的一个基本参数,被广泛应用于地震危险性分析研究中,因此有必要开展各活动地块边界带最大震级分布特征研究.本文利用广义极值(GEV)分布对各边界带作极值统计分析,对于分布参数进行估计,比较极值分布99％分位数和b值截距法最大震级、历史最大震级,并计算各震级水平下地震平均复发周期与发震次数,探讨活动地块各边界带未来强震活动的危险性.     

地震震级－频度关系b值参数的稳健估计

文中依据古登堡－里克特（Ｇ－Ｒ）震级－频度公式，结合实例说明了参数ｂ值的稳健估计方法，并将其与最小二乘法及最大似然法比较。结果表明，一般情况下，稳健估算结果与实际数据分布较吻合，并且可以在不剔除低震级段和大震级段的偏离分布值的情况下，取得与选取合理震级上下限后最小二乘法一致的结果。此外，稳健估计结果的余差分析同样也得到地震震级－频度关系不是简单的对数线性关系的结果，而原始Ｇ－Ｒ公式需加二次以上的高次项才能较好地描述该特性。     

单值地震危险性图的参数不确定性和可变性的处理

基于主观判断和有限资料解释的概率地震危险性研究的输入(地震震源区、地震发生率、衰减函数等)都是不确定的。面对这些不确定性,我们要考虑:(a)如何“合理地”确定单值概率地震危险性图中所标示的地震运动水平;(b)计算得到的地面运动水平的不确定性在这种图上可以明确地表达到何种程度。如果单个震源区的地震发生率、古登堡-里克特公式中的b值,以及最大震级的最佳猜测估计值被看作是不相关的并且每一参数中的不确定性可以看作是关于估计值对称的,那么,使用这些最优估计计算的概率地面运动水平表达了最可能的值和近似的均值。另一方面,根据不同震源区方案计算的地面运动大小,可能集中于每一方案的不同的中心值周围,而且,根据最可能方案计算出的地面运动水平可显著不同于所有方案计算出的地面运动水平的平均。我们得出结论(因为至少在美国东部和中部对地震过程是缺乏了解的),要得到单一的危险性估计结果,假定未来地震将趋向于发生在过去已经发生过地震的地区,并且具有与过去相当的发生率是合理的(在美国西部当单条断层上存在古地震资料时,这些信息是可以结合使用的)。在有多种震源区方案的情况下,使用不同的地面运动衰减关系得到的危险性估计,彼此会显著不同。人们或许考虑采用数个衰减关系的平均,但是,这样一来技术上的问题(例如,如何平均那些采用不同距离量度的关系式,如何估计σ,等等)会使操作非常困难,而且当包含一个与其它衰减关系离散较大的衰减关系时,它对均值会产生很大的影响。为了简化和可操作起见,使用单一的、不被认为是“极值”的衰减关系似乎是合理的,并且这一做法已获得危险性分析人员的接受。     

青藏高原东北缘地震活动性广义帕累托模型的全域敏感性分析

基于广义帕累托分布构建地震活动性模型,因其输入参数取值难以避免不确定性,导致依据该模型所得的地震危险性估计结果具有不确定性。鉴于此,本文选取青藏高原东北缘为研究区,提出了基于全域敏感性分析的地震危险性估计的不确定性分析流程和方法。首先,利用地震活动性广义帕累托模型,进行研究区地震危险性估计;然后,选取地震记录的起始时间和震级阈值作为地震活动性模型的输入参数,采用具有全域敏感性分析功能的E-FAST方法,对上述两个参数的不确定性以及两参数之间的相互作用对地震危险性估计不确定性的影响进行定量分析。结果表明:地震危险性估计结果（不同重现期的震级重现水平、震级上限及相应的置信区间）对两个输入参数中的震级阈值更为敏感;不同重现期的地震危险性估计结果对震级阈值的敏感程度不同;对不同的重现期而言,在影响地震危险性估计结果的不确定性上,两个输入参数之间存在非线性效应,且非线性效应程度不同。本文提出的不确定性分析流程和方法,可以推广应用于基于其它类型地震活动性模型的地震危险性估计不确定性分析。      

基于芦山主震及余震的地震预警参数与快速震级估算研究

我国四川地区继08年汶川地震后,在2013年再次发生了7级强烈大地震,造成了重大人员伤亡并诱发了多处灾难性滑坡灾害.地震预警系统是防震减灾中最为行之有效的方法之一,其中对地震震级的快速估算是非常重要且困难的.目前世界上提出较早、应用较为广泛的震级估计的方法主要有:基于最大卓越周期(τ_p~(max))、特征周期(τ_c)和最大位移幅值(Pd)方法.本文利用2013年芦山主震及余震P波初期部分的信息,研究了地震震级快速估算中三个预警参数与震级的相关性,提出了四川地区的震级估计模型,然后对其进行验证分析,并与前人得出的震级估计模型进行对比和评价.结果表明:三种方法均能在短时间内(3 s)有效地进行震级估算,总体而言,P_d方法估算效果最优,τ_c方法次优,τ_p~(max)方法较弱.在震级较大的主震震级估计中,三种方法均出现明显的震级低估(震级饱和)现象.对于τ_p~(max)方法,基于芦山地震的估计模型与基于汶川地震的估计模型较为接近,但与南加州地区模型区别较大,可能与四川地区龙门山断裂带的复杂地质条件有关;而τ_c方法的估计模型大体均较为接近,更具有普适性和稳定性.在地震预警系统的实际应用中,由于目前我国仍处于刚起步阶段,尚未建立密集的地震监测台网系统,因此在短时间内难以得到较为准确的震源距,所以与震源距相独立的τ_p~(max)和τ_c两种算法则显得较为实用,其中τ_c方法略优于τ_p~(max)方法,同时能较好地满足地震预警系统的精度要求,因此推荐使用τ_c方法应用于四川地区地震预警系统中的快速震级估算.需要注意的是,地震预警参数具有一定的区域适用性,即在不同震源机制、地形和地质条件影响下,不同地方的震级估算回归模型存在一定的差别.     

交城断裂带北段最大潜在地震发震概率评估1

基于太原市目标区交城断裂带的定量研究,特别是对活动断裂上的古地震资料进行的系统、详细的分析与总结,建立了反映该断裂地震地质特点和运动学属性的复发模式和概率模型.引入震级-地表破裂长度、震级-震源破裂长度、震级-断层破裂面积以及震级-地震矩的经验关系进行震级估计,最后进行综合评估以确定交城断裂带北段潜在地震的最大震级.复发模式的建立兼顾了泊松和准周期两种模式,利用专家意见法组合相应的Poisson模型和BPT模型,计算活动断裂最大潜在地震的复发概率.结果表明,交城断裂带北段潜在地震最大震级为Ms7.2级,而未来50a、100a、200a发生Ms7.2级地震的概率分别为2.1％、4.0％和7.9％.     

丽江7.0级地震的余震震源参数研究

本从ＰＤＲ－２数字化近源（６．２ｋｍ≤△≤４２ｋｍ）台网记录的丽江余震中初选了７４个地震，初定了震中位置，震中方位角。用相应的台网数据处理技术，研究了丽江地震序列的体波谱。对震级２．５≤Ｍ≤５．７、地震矩２１．４０≤ｌｏｎｇＭ０≤２３．２８的地震，得出了丽江地区地震系列的震源参数；以及Ｐ波、Ｓ波的ｌｏｇＭ０与Ｍｄ的关系式分别为：ｌｏｇＭ０＾ｐ＝０．６２Ｍｄ＋１９．９３；ｌｏｇＭ０＾ｓ＝０．５９Ｍ     

应用贝叶斯方法估算场地参数

地震小区划工作的目的是对一定区域未来地震可能遭遇的地震破坏分布进行估计，而这项工作中对场地参数的估计是一项非常重要的基础性工作，也往往是在场地区划中投入最大的一项工作。如何能在投入较少的情况下较正确合理估计场地物理力学参数，是非常重要的。本文介绍应用后验贝叶斯方法进行场地物理参数的估计，这不仅可以充分利用工程地质资料，节约大量的地质勘探的投入，同时也使得场地物理参数的估计更加合理。     

巴颜喀拉地块及其周边地震危险性分析

巴颜喀拉地块及其周边因地震动活动频繁成为中国地震预测研究的重要目标区,我们利用广义极值(GEV)分布对该地区的最大震级进行了极值统计分析.本文首先对于分布参数进行估计,然后计算出重现水平、地震平均复发周期与发震概率和震级危险率等参数,并比较了GEV分布99％分位数和b值截距法最大震级、历史最大震级.最后基于Monte Carlo模拟方法验证了运用GEV分布对研究区域做极值统计分析的稳定性.经计算,该地区6.0级地震的平均复发周期为1.8年,8.0级以上巨震的复发周期仅为76.8年.百年的重现水平高达8.08级.另外,该地块未来发生5.0～7.0级地震的概率较大,说明该地块未来仍然会比较活跃.Monte Carlo模拟结果表明,利用GEV分布探讨该地块的地震危险性在一定程度上是稳定的.