中长期地震危险性概率预测中的统计检验方法Ⅰ:Molchan图表法

对中长期地震危险性概率预测结果的统计检验是地震预测预报的重要研究内容,采用通用的统计检验方法是促进地震预测理论、模型"无障碍"交流和发展的必要条件。本文通过对青藏高原东北缘地区PI算法和RI算法回溯性预测结果统计检验的实例,介绍了如何运用CSEP计划中已常规采用的Molchan图表法进行地震概率预测统计检验和分析,并分别考虑了网格权重和地震活动权重两种算法计算"时空占有率"的情况。结果表明,Molchan图表法可通过显著性水平α和概率增益Gain有效地评价概率预测模型的预测能力。此外,对于同样的预测结果,在考察不同的预测策略时得到的统计检验结果也可能不同。其中地震活动权重相对网格权重采用更为严格的统计检验。

前兆异常区地震危险性概率的评估方法

讨论了根据前兆观测估计地震危险性概率的问题，给出了单学科观测异常区和多学科观测异常区地震的危险性概率的评估方法。利用３ 个先验概率值，即地震背景概率—Ｐ（Ｂ）、孕震条件下某学科观测发生异常的概率—Ｐ（Ａｉ／Ｂ）、非孕震条件下某学科观测发生异常的概率—Ｐ（Ａｉ／Ｂ），就可计算：①异常区有震概率（也称有震报准概率）Ｐ（Ｂ／Ａ） 或Ｐ（Ｂ／ ∩ｎｉ＝ １Ａｉ） 和异常区无震概率（也称虚报概率）Ｐ（Ｂ／Ａ） 或Ｐ（Ｂ／ ∩ｎｉ＝ １Ａｉ）；②非异常区有震概率（也称漏报概率）Ｐ（Ｂ／Ａ）或Ｐ（Ｂ／ ∩ｎｉ＝ １Ａｉ）和无异常无震概率（也称无震报准概率）Ｐ（Ｂ／Ａ） 或Ｐ（Ｂ／ ∩ｎｉ＝ １Ａｉ）。通过几组模拟数据的计算，可初步得出结论：①对于单学科观测，降低非孕震异常概率，即减小干扰，可提高报准概率；②多学科观测同时出现异常，将提高报准概率，但漏报概率也将显著增加。

地震危险性概率分析方法:存在的问题和纠正(英文)

地震危险性概率分析(PSHA)是目前最广泛应用于地震灾害与风险性评估的方法。然而它在计算中却存在着一个错误:把强地面运动衰减关系(一个函数)的条件超越概率等同于强地面运动误差(一个变量)的超越概率。这个错误导致了运用强地面运动误差(空间分布特征)去外推强地面运动的发生(时间分布特征)或称之为遍历性假设,同时也造成了对PSHA理解和应用上的困难。本文推导出新的灾害计算方法(称之为KY-PSHA)来纠正这种错误。

基于ETAS模型的余震概率性地震危险性分析

余震概率性地震危险性分析对主震发生后的地震风险评估至关重要。经典的余震概率性地震危险性分析框架基于Reasenberg-Jones(R-J)模型建立,但R-J模型对复杂的地震序列可能存在预测失效的问题,因此,将对真实余震序列描述较好的传染型余震序列模型(epidemic-type aftershock sequence,ETAS)结合到余震概率性地震危险性分析框架中,以2008年汶川7.9级余震序列和2014年鲁甸6.2级余震序列为例,评估余震危险性的重要性以及探讨ETAS模型相对于R-J模型在余震概率性地震危险性分析中的差异。结果表明:余震危险性水平与主震震级大小有关,主震震级越大,余震危险性水平越高,持续时间越久,对于震级较大的强主震事件,主震发生后早期的余震危险性不容忽略。对于汶川7.9级和鲁甸6.2级地震事件,在余震发生率预测效能上ETAS模型优于R-J模型,而在余震危险性评估中,余震危险性的大小取决于余震发生率,基于R-J模型计算的余震发生率大于ETAS模型计算的余震发生率,因此,在余震危险性评估中使用R-J模型相比于ETAS模型会造成对余震危险性水平的高估,ETAS模型可用于改进余震概率性地震危险性评估框架。

统计地震学在地震危险性概率预测方法研究中的应用与讨论

地震是给人民生命财产造成重大威胁的自然灾害之一,地震预测是推动现代地震学学科建设与发展的原动力.在地震的确定性预测迄今仍面临诸多困难的情况下,寻找地震发生的时空统计特征,构建地震危险性统计预测模型成为了评估地震风险并应用于城市与生命线工程等抗震规划设计实际需求的重要途径. 21世纪以来,以2008年汶川8级地震为代表的系列强震在一些低风险预测区发生,引起了世界各国地震学家对传统的地震危险性预测方法的反思.文章对统计地震学的发展进行了简要综述,着重分析了统计地震学在地震危险性预测中的研究成果及存在的问题,对其发展方向进行了讨论.分析指出:在大多数地区以现代地震目录为基础资料的地震危险性预测是有效的,尤其是在大型城市群地区及厚沉积、弱震区的地震危险性区划图的编制中,引入基于传染型余震序列模型的地震危险性概率预测应该是最为简便和有效的方法.

概率地震危险性分析、分解与设定地震及其在西安地区的应用

本文以Arc GIS为平台,以西安地区为研究对象,以西安地震局提供的信息为依托,给出西安地区的潜在震源区划分以及地震活动性参数,根据随机有限断层模拟得到西安地区地震动强度预测模型,基于离散算法平台,以潜在震源区为基本单位,将潜在震源区离散成栅格并看成点源,分别计算每个栅格对场地的地震动超越概率,即概率地震危险性分析。在此基础上进行概率地震危险性分解去确定某一最可能造成所需危险水平下的地震动强度参数的贡献率,最后根据加权方法确定设定地震。研究结果表明:在0. 2 s、1. 0 s和2. 4 s这3个周期的谱加速度值,相同危险性水平下,对场地贡献较大的潜在震源基本相同,但是各潜在震源的相对贡献率依然有差别;对于不同的危险性水平,对场地贡献较大的潜在震源区的数目区别很大,但基本都是靠近场地附近的几个潜在震源,且比较集中。

基于Newmark模型的概率地震滑坡危险性分析方法研究——以甘肃天水地区为例

本文以天水地区为研究区,结合地震潜在震源区模型和Newmark位移预测方程,采用概率地震危险性分析方法,计算了该地区50年超越概率10%水平下的Newmark位移。同时,根据天水地区50年超越概率10%下的阿里亚斯烈度,并结合Newmark位移与阿里亚斯烈度的关系式,计算了天水地区在遭受50年超越概率10%下的阿里亚斯烈度影响时,潜在滑坡体产生的Newmark位移分布。通过比较上述两种方法得到的天水地区不同Newmark位移的分布特征,本文认为二者虽然存在较大差异,但其空间分布特征均能反映天水地区每个场点处的相对滑坡危险性。对滑坡危险性水平进行分区的结果显示,天水地区60%以上的区域具有高地震滑坡危险性,50%以上的区域具有甚高地震滑坡危险性。本文的研究结果可以作为天水地区地震危险性及风险评估的参考资料,也可以作为天水地区城市规划、土地使用规划、地震应急准备以及其它公共政策制定的参考资料。

余震时空丛集对概率地震危险性分析的影响

本文以东北、华北及川滇地区为例,系统研究了余震时空丛集对概率地震危险性分析的影响.采用基于传染型余震序列模型(ETAS)的蒙特卡罗模拟方法,模拟了包含余震和不包含余震的两套地震序列,然后以模拟地震目录为基础输入,采用基于空间光滑地震活动性模型的地震危险性分析方法计算了两套地震危险性结果——PGA(Peak Ground Acceleration,峰值加速度),通过分析比较这两套PGA的绝对差值和相对差值来研究余震时空丛集对概率地震危险性分析的影响.研究结果表明余震对50年超越概率10%地震危险性计算结果的影响均值为6%左右,最大可达10%,并且随着超越概率水平的提高,余震影响也越大.弱地震活动区余震对概率地震危险性分析的影响要高于强地震活动区.研究结果还进一步揭示两套PGA结果绝对差值的最大值约为15 cm·s^(-2),且出现在高PGA区,这意味着余震对概率地震危险性计算结果不会产生显著影响.因此在地震区划或一般性地震危险性分析中可考虑不用删除余震.

基于MAPGIS的概率性地震危险性分析

对工程场地的地震危险性分析是地震安全性评价的主要方法,并且为工程师提供抗震参数.在分析国内概率性地震危险性分析(PSHA)方法基础上,提出了基于GIS的概率性地震危险性分析的可行性方案.通过MAPGIS二次开发编写了地震危险性分析程序,并以三峡坝区某工程场地为例对程序进行了测试.基于MAPGIS的概率性地震危险性分析程序提供友好的人机交互界面,提高了地震危险性分析的可操作性,更重要的是帮助用户从空间数据中挖掘更多的信息.

基于阿里亚斯烈度估值的概率性地震危险性分析——以四川丹棱县及其周缘为例

以编制《中国地震动参数区划图》时所构建的地震潜在震源区模型和地震活动性模型为基本输入,在四川省丹棱县及其周缘地区开展基于阿里亚斯烈度的概率地震危险性分析,计算50年超越概率10%的阿里亚斯烈度(Arias Intensity,Ia)。结果表明:丹棱县及其周缘绝大部分地区的Ia值都在0.11m/s以上,部分地区在0.32m/s甚至0.54m/s以上,具有较高的地震诱发滑坡风险,应当加强人员密集区的地震诱发滑坡危险性评估。根据不同Ia预测方程计算得到的Ia分布有较大差异,因此在计算Ia时应考虑多个Ia预测方程,对最终结果进行加权平均,以减小Ia结果的不确定性。同时还发现Ia值与1.0s的谱加速度具有很好的相关性,这也印证了Ia和1.0s谱加速度与砂土液化的相关性。

基于特征地震模型含时间的概率地震危险性分析方法及其应用研究

本文介绍了特征地震的对数正态分布模型、正态分布模型和布朗过程时间模型,提出了使用地震破裂面源模型的特征地震含时间的概率地震危险性分析理论和方法.通过具体算例对不同的特征地震模型进行了比较,并对特征地震危险性分析方法进行了系统探索.研究结果表明,特征地震含时间模型在复发周期早期的地震危险性低于不含时间模型,而在后期其地震危险性则高于不含时间模型.特征地震复发周期的对数正态分布模型与布朗过程时间模型计算得出的地震危险性差别不大.在未到期望复发时间时,正态分布模型与前两种模型计算的地震危险性差别不大;而接近期望复发时间及之后时段,正态分布模型计算的地震危险性则迅速增大.

潜在地震破裂面源模型及在概率地震危险性分析中的应用方法

目前国内外常用的概率地震危险性分析方法中,一般把未来可能发生的地震震源(潜在震源)简化为点源或者是线源.但对于较大震级的地震,仍然采用点源或线源模型来描述潜在震源显然是不合理的.为此,在概率地震危险性分析方法中本文提出了潜在地震破裂面源模型,并以1999年9月21日台湾集集7.6级地震的车笼埔断层为例,探讨了其在地震危险性分析中的应用方法.结果表明:①采用潜在地震破裂面源模型是合理的,因为它可以模拟地震破裂面与地震动影响场的三维展布特征,尤其适用于较大震级地震的近场区域;②潜在地震破裂面源的大小、产状,对近震源场点的地震危险性分析和地震区划结果有明显的控制影响.

基于GIS的概率性地震危险性分析(英文)

介绍了地理信息系统的主要特征、概率性地震危险性分析的一般步骤及其在MapGIS软件上的实现方法。并以三峡坝区某工程场地为例对程序进行了测试。GIS的用户界面使用户易于掌握概率性地震危险性分析的流程和数据。更重要的是它提高了用户对复杂地震危险性分析过程的认知能力,可以帮助用户在地震危险性分析中从空间数据中挖掘更多的信息。

拟负刚度控制隔震结构的概率地震危险性评估

考虑地震的不确定性,运用概率地震危险性评估的方法对拟负刚度隔震结构的性能进行了评估.以隔震Benchmark模型为研究对象,将拟负刚度(PNS)控制系统与普通双线性隔震(BIS)系统的性能从3个方面进行了对比:增量动力分析、概率地震易损性和概率地震危险性,多方位探讨PNS控制对隔震结构的减震效果并定量分析其控制性能的优越程度.分析结果表明,PNS控制系统楼层最大加速度立即使用极限状态的年超越概率较BIS系统小30%,隔震支座最大位移倒塌极限状态的年超越概率较BIS系统小20%.由此可见,PNS控制更易保证隔震结构在小震下的功能性和大震下的安全性.

融合三维断层源和二维潜在震源区的随机抽样概率地震危险性分析算法研发

文中采用蒙特卡罗随机抽样方法,研发了一套融合传统二维潜在震源区(下文简称“潜源”)和三维断层源的概率地震危险性算法。该算法不仅适用于传统的区域面源,同时还能考虑地震的破裂尺度并兼容三维断层源的概率地震危险性计算。文中研发的算法可高效实现断层源地震事件集的三维模拟,并将地震破裂尺度引入到概率地震危险性计算中,显著提高了近断层地区地震危险性计算的合理性。为了提高程序的执行效率,算法采用预先在平面潜源中充填网格点的方式随机模拟地震事件在潜源内的均匀分布。对于椭圆衰减的地震危险性计算,算法采用了预先构建不同震级、距离及不同场点与潜源长轴方向夹角下的短轴距的三维矩阵,通过查表和插值方式直接获得相应场点的短轴距,避免了循环迭代逼近短轴距计算效率低下的问题。分别利用五代图的概率地震危险性程序和文中研发的算法,计算了湖南长-株-潭(长沙-株洲-湘潭)城市群所处的中强地震活动环境的区域地震危险性以及近断层源的常德、株洲2个场点在不同概率水平下(重现期分别为50a、 475a和2 475a)的地震危险性。比较研究表明,五代图的程序低估了三维断层源附近的地震危险性,且随着概率水平的降低,低估的程度越来越高。最后,利用太平洋地震工程中心(Pacific Earthquake Engineering Research Center, PEER)验证概率地震危险性程序的算例(数据集1案例10)验证了文中算法的可靠性。

地震动参数之阿里亚斯烈度和累积绝对速度的特征及其在概率地震危险性分析中的应用

地震动参数是表征地震动特性(振幅、频率和持时)的重要度量参数。目前使用最多的地震动参数主要有峰值加速度和谱加速度等,然而,任何单一地震动参数都不能全面有效地反映所有的地震动破坏特征(如峰值加速度仅反映振幅)。为了更加全面地描述地震动特征及其潜在的破坏作用,满足工程应用的需求,科学家们一直在寻找能够表征地震动多种属性和特定地震破坏的地震动参数。经研究发现,阿里亚斯烈度Ia(Arias Intensity)、累积绝对速度CAV(Cumulative Absolute Velocity)及Newmark位移在地震诱发滑坡预测、沙土液化评估、地震灾害损失估计以及地震危险性分析等方面优于峰值加速度,因此,开展Ia、CAV和Newmark位移的特征研究及与其相关的地震危险性分析方法研究对于拓展我们对地震动特性的认识和提高地震灾害预测水平具有重要的科学意义和应用价值。本文我们根据汶川地震的177个台站的528条可用记录的强震数据,在对地震记录进行基线校正后计算了地震动参数Ia、CAV及Newmark位移,并研究了这些地震动参数的相关特征,最后我们还开展了基于Ia、CAV及Newmark位移的概率地震危险性分析方法研究。我们研究了汶川地震la的空间分布特征、衰减特征、近断层特征以及la与其他地震动参数的相关性。在空间分布特征研究中,我们比较了竖向和水平向la值的差异,统计了二者的比值分布及随距离的变化特征;对于垂直断层走向和平行断层走向的la值也做了同样的工作。在衰减特征方面,我们主要研究了3种场地类型下的la衰减特征及其与现有la预测模型的比较,指出了他们之间的差异并分析了导致差异的原因。在近断层特征方面,主要研究了上盘效应和方向性效应对1a的影响。最后我们研究了la与PGA和Housner烈度的相关性,建立了la与PGA和Housner烈度之间的经验关系式,并分析了场地条件对地震动参数之间相关性的影响。我们对CAV也同样研究了其空间分布特征、衰减特征、近断层特征以及CAV与其他地震动参数的相关性。特别地,我们还研究了标准CAV的衰减特征和其与PGA的相关性,分析了场地条件对它们的影响。另外,我们简要分析了汶川地震Newmark位移的衰减特征及其与PGA和la的相关性。本文中我们还开展了基于Ia、CAV的概率地震危险性分析方法研究,开展了基于Newmark位移的概率地震滑坡危险性分析方法研究,梳理了相关的理论、方法,建立了相关的技术框架,编制了计算程序并在典型地区进行了试算。本文的研究对我们认识汶川地震乃至我国西部地震地质构造环境下地震的Ia、CAV及Newmark位移的相关特征具有重要的科学意义。建立的相关概率地震危险性分析方法的技术框架对于提高我国的概率地震危险性分析水平和地震灾害预测水平具有重要的应用价值。

基于自适应空间光滑模型和三维断层模型的概率地震危险性分析方法研究

概率地震危险性分析(PSHA)直接为工程建设提供抗震设计参数,其重要性不言而喻。进行相关研究的核心科学问题在于,如何优化参数、改进模型来更好地描述地震的中长期活动规律和破坏特征,从而提高结果的可靠性。目前在该领域,我国现行的方法体系与美国国家地质调查局(USGS)提出的基于空间光滑模型和三维断层模型的方法存在较大差异,两种方法的特征和适用性也存在一定的差异。因此,本人在各位老师的指导下开始学习国外先进的理论和技术,将其用于国内研究区,并与国内常用的方法进行对比分析,总结各自的特点。在此基础上尝试做一些改进和创新,目前已经取得部分进展和成果。

对概率地震危险性分析方法中地震动离散截断位置影响的讨论

概率地震危险性分析方法（PSi）衰减关系不确定（即地震动离散）校正，采用的是截断的正态分布模型，在平均值正负ε倍标准差范围内考虑其对场址地震动影响。ε取值对极低年超越概率（APE），如104或更低的地震动参数具有重要影响。本文将从数学分析、场址分析、地震动离散分析等方面分析ε值对PSHA结果的影响，并给出ε值取值的建议。

时间相依的概率地震危险性分析方法研究及应用

对使用不同模型的郯庐断裂带莒县-郯城段和山西交城断裂的地震危险性进行了计算，并比较了时间独立模型和不同时间相依模型的计算结果。传统的时间独立（time-independent）的概率地震危险性分析（probabilistic seismic hazard analysis ， PSHA）模型假设地震的发生服从泊松分布。

基于地震带细致划分为地震构造区的概率地震危险性分析

在我国当前地震安全性评价中,普遍使用考虑地震活动时空不均匀性的概率地震危险性分析方法(CP-SHA),它规定以地震带为统计区域计算地震活动性参数b值和v4值。正在编制的中国地震动区划图(五代图)提出了针对潜在震源区进行三级划分原则:划分地震带、地震带上划分地震构造区、地震构造区内再划分潜在震源区。本文提出以地震构造区为统计区域回归统计方法计算b值和v4值,然后进行概率危险性分析计算,这样得到的结果可能更为合理。