中国海啸研究-概率海啸危险性分析(PTHA)

中国海啸

概率海啸危险性分析

早在20世纪七十年代，国外就已经开始了概率海啸危险性分析的研究。在1978年Houston和Garcia(1978)进行了开创性的海啸危险性研究，他借鉴Cornell(1968)提出的PSHA方法，将概率方法与简单的地震发生模型结合，计算了夏威夷群岛附近海啸波高的发生频率。但当时的研究存在很多的局限性，只考虑的远场地震源模型，没有考虑近场地震源模型，采用的水深地形数据精度很低。随后Lin和Tung(1982)借鉴Cornell提出的概率地震危险性分析方法，提出概率海啸危险性分析的概念。Rikitake和Aida(1988)提出了使用海啸历史记录和典型地震断层模型估计海啸波高在海岸线附近的超越概率。Downes和Stirling(2001)1提出在概率海啸危险性分析中使用经验衰减关系近似海啸波在传播过程中的衰减，但由于这种经验衰减关系只适用于特定的震源和场点，许多地区没有足够的数据让研究人员获取经验衰减关系，这种方法并未被广泛采纳。

随着流体力学的发展和计算机计算能力的提升，海啸数值模拟技术得到了长足的发展，越来越多的研究人员采用数值模拟的方式模拟海啸生成和传播的全过程，计算海啸危险性，解决了很多地区海啸记录匮乏的问题，使PTHA研究进入了新的时代。2006年Geist和Parsons(2006)将前人的研究进行了总结，提出将海啸历史记录与数值模拟结果相结合，分析海啸风险，并采用逻辑树方法解决PTHA中不确定性问题，使PTHA得到进一步发展。随后，在全世界范围内，针对不同地区的PTHA研究不断涌现。

概率地震危险性分析(PSHA)与概率海啸危险性分析(PTHA)对比

传统的PTHA，把震级范围分档，然后通过震中随机采样，计算出每一档地震情况下的波高超越概率。温瑞智根据我国的地震潜源特点和潜源参数，提出了我国地震海啸危险性分析的主要步骤。杨智博针对我国面临的海啸潜源，利用Monte Carlo对震级和震中随机取样，改进了传统的PTHA。总结改进后的PTHA方法流程如下：

第一步：确定地震海啸潜源，选取研究区域周边的地震潜源。结合其震级上限、地质构造背景及历史地震活动性等信息，选择可能诱发海啸并对研究区域产生影响的地震潜源，作为地震海啸潜源。

第二步：对所选取的地震海啸潜源进行活动性参数分析。针对局地海啸潜源和区域海啸潜源采用不同的方法计算能够诱发海啸的地震年发生率。并得到每个地震海啸潜源震级的累积分布函数。

第三步：使用Monte Carlo技术采样制定地震海啸目录。假设第i个地震海啸潜源随机发生N_i次地震，针对每个地震海啸潜源使用Monte Carlo技术进行N_i次采样，使样本的震级分布满足第二步求得的震级累积分布函数，震中位置满足在地震海啸潜源上均匀分布(将地震海啸潜源考虑为线源)。断层破裂的长度、宽度、滑移量等根据与震级相关的经验公式确定。将Monte Carlo技术采样出的包含震级、震中位置、破裂长度、宽度、滑移量等参数的样本作为地震海啸目录

第四部：数值模拟。对第三步制定的地震海啸目录中N_i次地震，进行数值模拟，模拟每一次地震诱发海啸从生成到传播到我们关注场点的全过程，记录每一次模拟，关注场点海啸波高的最大值。

第五步：计算每个场点最大海啸波高超越概率。统计i个地震海啸潜源发生N_i次地震海啸情况下，我们所关注的每个场点的最大海啸波分布情况，使用相应统计分布函数进行拟合，结合第二步求出的每个地震海啸潜源的地震年发生率，最终计算出每个场点最大海啸波高超越概率。

传统PTHA流程

改进PTHA流程