一种基于概率方法的南海海啸危险性预测方法
【技术领域】：
[0001] 本发明属于灾难预测，特别涉及一种灾难预警仿真应用中基于概率方法的南海海 啸危险性预测方法。
【背景技术】
[0002] 海啸危险性，是指靠近海岸的特定地段遭受特定规模的海啸袭击的可能性，这种 可能性可以表示为特定地点在特定时间段内波高超过特定数值的海啸的发生概率。目前， 地震危险性分析已广泛应用于地震区划分和具体地震的估计。可以说，我国地震危险性分 析的研究已趋于成熟，产生的研究价值也很高，不仅能为城市规划、重大工程建设提供技术 支持，还能为政府部门在某些重大决策的讨论中提供重要的技术参考，为我国地震时减少 人员伤亡和财产损失及地震防灾减灾体系的构建做出了很大贡献。
[0003] 相对于地震，海啸的危害也不可忽视，2004年苏门答腊海啸就是最好的说明。基于 概率方法的海啸危险性是通过充分借鉴地震危险性的方法和思路，并从历史数据、数值模 拟的角度提出的海啸分析方法。
[0004] 地震危险性方法在上世纪六十年代被提出来，而我国在这方面的研究始于上世纪 八十年代末。基于地震危险性分析方法，学者们把该理论思想应用于海啸方面的研究，这样 就形成了概率性海啸危险性分析方法。该分析方法的优点是考虑了地震海啸发生时间的不 确定性和它的可能性。地震危险性分析方法目前在世界范围已经广泛使用，主要思路是基 于Cornell提出的方法，模型是地震危险性分析模型。地震活动发生在一定的空间范围内， 即这里指的就是潜在震源区：该区内的大小地震比例关系符合古登堡-里克特的震级-频 度（G-R)关系，在时间上满足泊松分布。潜在震源区域内各点地震发生的概率相等，符合均 匀分布。我国目前地震危险性分析理论是胡聿贤基于以上理论进行了一定的修正。尽管此 方法思路明确，但地震和海啸灾害的本质有所不同，以及历史海啸资料也不够完善，所以海 啸危险性分析方法的步骤也不完全相同。
【发明内容】
[0005] 本发明的目的在于现有南海海啸历史资料不完善等原因，所以将提出一种基于概 率方法的南海海啸危险性预测方法，即给出特定地点在波高超过特定数值的海啸发生概 率。
[0006] 本发明具体内容：
[0007] -种基于概率方法的海啸危险性预测方法，该方法包括以下具体步骤：
[0008] 1)确定南海潜在的海啸源区。根据南海板块地质情况、历年海啸记录、南海水深情 况等划分潜在海啸源区。结合地形条件、地质构造条件等，马尼拉海沟被认为是南海风险最 商的区域；
[0009] 2)基于概率方法，通过统计南海历史海啸记录，给出相关地震参数。基于概率方 法，通过历年来南海海啸地震的记录数据，统计出潜在海啸源区的古登堡-里克特震级-频 度关系，然后给出相关参数0、b值等。因此，可得到马尼拉海沟海啸源区的不同震级的累 积频率分布为：
【主权项】
1. 一种基于概率方法的南海海啸危险性预测方法，该方法有以下具体步骤： 1) 确定南海潜在的海啸源区。根据南海板块地质情况、历年海啸记录、南海水深情况等 划分潜在海啸源区。结合地形条件、地质构造条件等，马尼拉海沟被认为是南海风险最高的 区域； 2) 基于概率方法，通过统计南海历史海啸记录，给出相关地震参数。基于概率方法，通 过历年来南海海啸地震的记录数据，统计出潜在海啸源区的古登堡-里克特震级-频度关 系。因此，可得到马尼拉海沟海啸源区的不同震级的累积频率分布； 3) 给出南海区域产生地震的样本。根据马尼拉海沟海啸源区的不同震级的累积频率分 布的关系，在马尼拉沟海啸源区内空间上均匀的产生地震。低震级可以由单一断层作为震 源，而高震级地震均由组合断层作为震源。高震级地震均由组合断层作为震源，组合的规律 是组合断层长度在日本规范和Papazachos B. C.回归公式计算的两者之间； 4) 通过数学模型计算海啸传播到近岸的波高。在潜在海啸源区内，对每个地震事件运 用海啸数值模式计算沿岸波高值； 5) 计算马尼拉海沟海啸源区的累积频率。以特定波高值为界，针对潜在海啸源区内地 震事件进行数值模拟的结果，统计超过特定波高个数，给出该场点波高的累积频率； 6) 计算马尼拉海沟海啸源区的超越概率和重现期。通过指定场点的波高对应的累积频 率可以得到相对应的超越概率和重现期，最后给出海啸波高对应的超越概率曲线，即给出 特定场点在特定时间段内波高超过特定数值的海啸的发生概率。
2. 根据权利要求1所述的一种基于概率方法的海啸危险性预测方法，该方法有以下具 体步骤： 1) 确定南海潜在的海啸源区。根据南海板块地质情况、历年海啸记录、南海水深情况等 划分潜在海啸源区。结合地形条件、地质构造条件等，马尼拉海沟被认为是南海风险最高的 区域； 2) 基于概率方法，通过统计南海历史海啸记录，给出相关地震参数。基于概率方法，通 过历年来南海海啸地震的记录数据，统计出潜在海啸源区的古登堡-里克特震级-频度关 系，然后给出相关参数P、b值等。因此，可得到马尼拉海沟海啸源区的不同震级的累积频 率分并T -
(1) 式中，M。和Mt分别为震级上限和震级下限； 3) 给出南海区域产生地震的样本。根据马尼拉海沟海啸源区的不同震级的累积频率分 布的关系，在马尼拉沟海啸源区内空间上均匀的产生地震，震级范围在M t与M。之间。低震 级可以由单一断层作为震源，而高震级地震均由组合断层作为震源。高震级地震均由组合 断层作为震源，组合的规律是组合断层长度在日本规范和Papazachos B.C.回归公式计算 的两者之间； 4) 通过数学模型计算海啸传播到近岸的波高。在马尼拉海沟海啸源区内，对于每一个 地震事件给出由经验公式计算的断层面参数，例如断层面长度、宽度、深度、滑移量等参数。 基于以上参数，通过C0MC0T海啸模型计算每个地震事件产生的沿岸波高值，或者给出指定 场点K个沿岸波高值。COMCOT海啸模型可以计算从海啸产生到海啸传播至近岸这一过程。 为了考虑到计算的精度以及效率，该海啸模型采用了网格嵌套方法。 海啸波在大洋传播时，如果计算区域较大，可以用球面坐标下的线性浅水方程来模拟， 公式如下：
式中，n代表水面高度，（)是纬度和经度，R是地球半径，g是重力加速度，h是 水深。（P，Q)表示在x和y方向上的体积通量，他们是速度和水深的乘积，也就是P = hu， ch Q = hv。表示的是瞬时海底运动，可以用在计算滑坡引起的海啸。f代表科氏力系数or
式中，Q是地球转速。 当计算区域较小时，应该采用笛卡尔坐标下的浅水方程。公式如下：
当海啸传到大陆架，接近海岸时，球面坐标下的非线性浅水方程如下：
笛卡尔坐标下的非线性浅水方程：
式中，H是总水深，H= n+h。FjPFy分别表示x和y方向上的底摩擦，可以通过曼宁 公式计算：
式中，n是曼宁系数。 COMCOT采用显式蛙跳有限差分数值方法。水面高度和体积通量（P，?在时间和空间上 交错排列。球面坐标下，显式蛙跳有限差分格式的浅水方程可以表示为：
笛卡尔坐标下，可以表示为：
5) 计算马尼拉海沟海啸源区的累积频率。通过计算海啸源区内地震事件，以某一场点 波高值hu为界，统计超过波高值h u的个数为k，该场点波高h u的累积频率为：
6) 计算马尼拉海沟海啸源区的超越概率和重现期。指定场点T年波高hu的超越概率 可表不为：
那么，波高hu的重现期可由下式计算：
最后给出海啸波尚对应的超越概率曲线。
【专利摘要】一种基于概率方法的南海海啸危险性预测方法，对于没有足够历史海啸观测记录的地区，只能采用基于数值模拟的分析方法来评价该地区的海啸危险性。概率性海啸危险性分析是通过充分借鉴地震危险性分析的方法和思路，并从数值模拟的角度提出概率性海啸危险性分析的方法。数值计算采用COMCOT海啸模型，使用多层大中小网格嵌套的方法计算海啸从远海传至近海。以风险区域最高的马尼拉海沟海啸源区为例计算地震引发海啸对南海海域的影响，给出该海域某场点波面升高对应的年累积频率、超越概率和重现期。所以本发明可以给出特定地点在特定时间段内波高超过特定数值的海啸的发生概率，即基于概率方法的海啸危险性预测方法。
【IPC分类】G01V1-30, G06F19-00
【公开号】CN104615847
【申请号】CN201410814810
【发明人】黄明汉, 张文忠, 吕迎雪, 黄宣军, 毕琳, 孔丛颖, 郭立栋
【申请人】中交天津港湾工程研究院有限公司, 中交第一航务工程局有限公司, 天津港湾工程质量检测中心有限公司
【公开日】2015年5月13日
【申请日】2014年12月24日