(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211442287.4 (22)申请日 2022.11.18 (71)申请人 海南浙江大学研究院 地址 572025 海南省三 亚市崖州区崖州湾 科技城用友产业园1 1号楼 (72)发明人 白晔斐　刘金伟　周一帆　房欣　 刘英男　魏笑然　支泓欢　 (74)专利代理 机构 北京专赢专利代理有限公司 11797 专利代理师 李斌 (51)Int.Cl. G06F 30/20(2020.01) G01C 13/00(2006.01) G06F 111/10(2020.01) (54)发明名称 一种基于响应面的海底滑坡海啸的快速预 报方法 (57)摘要 本发明涉及滑坡海啸技术领域， 具体涉及一 种基于响应面的海底滑坡海啸的快速预报方法； 本发明使用最小二乘法对滑坡深度、 滑坡体体积 量和海啸最大波高三个参数进行了曲面拟合， 构 建了以滑坡深度和滑坡体体积量为自变量， 以海 啸最大波高为因变量的公式， 完成了响应面的构 建， 使对滑坡海啸的预报能够较为准确且高效的 完成， 采用 构建响应面的方法， 只需对特定的数 十个滑坡海啸情景进行模拟， 便可求得响应面范 围内所有情景的海啸最大波高， 使得到的滑坡海 啸的数据量由点变为面， 大大降低了模拟滑坡海 啸的计算成本和时间成本 。 权利要求书2页 说明书5页 附图1页 CN 115495941 A 2022.12.20 CN 115495941 A 1.一种基于响应面的海底 滑坡海啸的快速预报方法， 其特 征在于， 包括以下步骤： （1） 判断会发生海底 滑坡的位置； （2） 获得包 含滑坡位置的研究区域的高精度地形文件； （3） 根据海底滑坡所在的深度范围及滑坡体体积量的区间， 将深度和滑坡体体积量进 行离散化， 离散的区间个数分别为 和 ； （4） 建立 一个包含所有滑坡参数的目录； （5） 对 个滑坡情景进行滑坡海啸的模拟， 得到 个研究范围内的海啸最大 波高数据； （6） 对研究范 围内的所有像素点的最大波高进行统计， 每个像素点受到深度 和滑坡 体体积量 两个参数的影响， 会产 生 个海啸的最大波高 ， 对这三个变量的数 据进行拟合； （7） 构建出每 个像素点的滑坡 海啸最大波高响应面； （8） 当有海底滑坡发生时， 判断该滑坡发生的深度 和滑坡体体积量 是否在构建的 响应面内； （9） 获取整个区域的海啸最大波高， 对海啸波高超过阈值的高风险区域进行快速预警， 并对滑坡 海啸的最大波高数据进行 可视化处 理。 2.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （1） 中， 判断会发生海底滑坡的位置的具体方法为： 根据对富含天然气水合物区 域位置计算得到的滑坡安全系 数分布云图， 对安全系 数小于设定阈值的区域绘制轮廓线， 轮廓线内的区域即为会发生海底 滑坡的位置 。 3.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （2） 中， 在GEBCO网站中获得高精度地形文件。 4.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （4） 中， 在目录中深度和滑坡体体积量的大小取值于各区间的中间数， 海底滑坡 的情景数量为 。 5.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （5） 中， 采用Navier ‑Stokes海啸三维数值模型对 个滑坡情景进行滑坡海 啸的模拟。 6.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （6） 中， 拟合公式为： ； 式中， 和 均为待定系数。 7.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （7） 中， 构建出每 个像素点的滑坡 海啸最大波高响应面的过程 为： 计算获得每个像素点的拟合公式的系数； 将系数代入拟合公式， 获得每个像素点的海权　利　要　求　书 1/2 页 2 CN 115495941 A 2啸最大波高公式， 自变量为滑 坡深度 和滑坡体体积量 ； 构建出每个像素点的滑坡海啸 最大波高响应面。 8.根据权利要求7所述的一种基于响应面的海底滑坡海啸的快速预报方法， 其特征在 于， 在步骤 （8） 中， 当海底滑坡的深度 和滑坡体体积量 在响应面的影响范围内时， 将海底滑坡的发 生的深度 和滑坡体体积量 输入到步骤 （7） 构建的响应面中， 求得整 个区域的海啸最 大 波高数据； 当海底滑坡的深度 和滑坡体体积量 不在响应面范围内， 使用快速建模的方法对 滑坡海啸进行预报， 并扩大响应面滑 坡的深度 和滑坡体体积量 的范围， 并继续进行滑 坡海啸的情景模拟， 将响应面的影响范围进一 步扩大。权　利　要　求　书 2/2 页 3 CN 115495941 A 3