一种基于响应面的海底滑坡海啸的快速预报方法技术

本发明专利技术涉及滑坡海啸技术领域，具体涉及一种基于响应面的海底滑坡海啸的快速预报方法；本发明专利技术使用最小二乘法对滑坡深度、滑坡体体积量和海啸最大波高三个参数进行了曲面拟合，构建了以滑坡深度和滑坡体体积量为自变量，以海啸最大波高为因变量的公式，完成了响应面的构建，使对滑坡海啸的预报能够较为准确且高效的完成，采用构建响应面的方法，只需对特定的数十个滑坡海啸情景进行模拟，便可求得响应面范围内所有情景的海啸最大波高，使得到的滑坡海啸的数据量由点变为面，大大降低了模拟滑坡海啸的计算成本和时间成本。啸的计算成本和时间成本。啸的计算成本和时间成本。

【技术实现步骤摘要】
一种基于响应面的海底滑坡海啸的快速预报方法


[0001]本专利技术涉及滑坡海啸
，具体涉及一种基于响应面的海底滑坡海啸的快速预报方法。

技术介绍

[0002]除海底地震这一引发海啸的最主要因素外，海底滑坡是引发海啸又一重要因素，从近期的记录中发现，在所有的海啸事件中，将近30%的海啸和海底滑坡等海底沉积物运动有关，其中大型的滑坡海啸会对沿岸人类活动和海洋生态造成巨大破坏，也会对人类的海底工程设备造成不可挽回的损失。因为海底滑坡不易观测且海底地质和地形复杂，故人们对海底滑坡的认识较少，随着近几年观测技术的发展，滑坡海啸才逐渐引起人们的重视。
[0003]我国陆架宽阔且存在大量富含天然气水合物的区域，极有可能发生海底滑坡，且滑坡位置离我国沿岸较近，迫切需要开展滑坡海啸的评估，构建一种滑坡海啸的快速预报方法，为滑坡海啸的危险性分析和灾害预警工作提供关键的参考和指导。
[0004]随着计算机算力与海啸数值模式的发展，预测和评估滑坡海啸变得愈发的可能。
[0005]现有技术对滑坡海啸的危险性评估多为使用数值模式对滑坡海啸进行特定的情景式模拟并进行分析，选定的滑坡参数来源于前人对滑坡海啸危险区域的研究，所以进行数值模拟的滑坡海啸情景的数量少且参数固定，还原现实中的滑坡海啸的概率非常小，无法起到对滑坡海啸的预警作用。
[0006]综上所述，研发一种基于响应面的海底滑坡海啸的快速预报方法，是滑坡海啸
中急需解决的关键问题。

技术实现思路

[0007]针对现有技术所存在的上述缺点，本专利技术在于提供一种基于响应面的海底滑坡海啸的快速预报方法，本专利技术所提供的方法，使用滑坡深度和滑坡体体积量作为自变量，海啸最大波高作为因变量进行曲面拟合，构建响应面，通过构建响应面，降低滑坡海啸模拟次数，且还能准确且高效地对滑坡海啸进行预报。
[0008]为实现上述目的，本专利技术提供了如下技术方案：本专利技术提供了一种基于响应面的海底滑坡海啸的快速预报方法，包括以下步骤：（1）根据研究资料，判断会发生海底滑坡的位置。
[0009]（2）获得包含滑坡位置的研究区域的高精度地形文件。在具体操作中，可以采用在GEBCO网站中获得高精度地形文件。
[0010]（3）根据海底滑坡所在的深度范围及滑坡体体积量的区间，将深度和滑坡体体积量进行离散化，离散的区间个数分别为和。
[0011]（4）建立一个包含所有滑坡参数的目录，优选的，在目录中深度和滑坡体体积量的大小取值于各区间的中间数，海底滑坡的情景数量为。
[0012]（5）采用Navier
‑
Stokes海啸三维数值模型对个滑坡情景进行滑坡海啸的模拟，对个滑坡情景进行滑坡海啸的模拟，得到个研究范围内的海啸最大波高数据。
[0013]（6）对研究范围内的所有像素点的最大波高进行统计，每个像素点受到深度和滑坡体体积量两个参数的影响，会产生个海啸的最大波高，对这三个变量的数据进行拟合；所述拟合公式为：；式中，和均为待定系数。
[0014]（7）构建出每个像素点的滑坡海啸最大波高响应面，具体过程为：计算获得每个像素点的拟合公式的系数；将系数代入拟合公式，获得每个像素点的海啸最大波高公式，自变量为滑坡深度和滑坡体体积量；构建出每个像素点的滑坡海啸最大波高响应面。
[0015]（8）当有海底滑坡发生时，判断该滑坡发生的深度和滑坡体体积量是否在构建的响应面内；若海底滑坡的深度和滑坡体体积量在响应面的影响范围内时，则将海底滑坡的发生的深度和滑坡体体积量输入到步骤（7）构建的响应面中，求得整个区域的海啸最大波高数据。
[0016]若海底滑坡的深度和滑坡体体积量不在响应面范围内，则使用快速建模的方法对滑坡海啸进行预报，并扩大响应面滑坡的深度和滑坡体体积量的范围，并继续进行滑坡海啸的情景模拟，将响应面的影响范围进一步扩大。
[0017]（9）获取整个区域的海啸最大波高，对海啸波高超过阈值的高风险区域进行快速预警，并对滑坡海啸的最大波高数据进行可视化处理。
[0018]有益效果采用本专利技术提供的技术方案，与已知的公有技术相比，具有如下有益效果：（1）本专利技术使用最小二乘法对滑坡深度、滑坡体体积量和海啸最大波高三个参数进行了曲面拟合，构建了以滑坡深度和滑坡体体积量为自变量，以海啸最大波高为因变量的公式，完成了响应面的构建，使对滑坡海啸的预报能够较为准确且高效的完成。
[0019]（2）本专利技术采用构建响应面的方法，只需对特定的数十个滑坡海啸情景进行模拟，便可求得响应面范围内所有情景的海啸最大波高，使得到的滑坡海啸的数据量由点变为面，大大降低了模拟滑坡海啸的计算成本和时间成本。
附图说明
[0020]图1为本专利技术一种基于响应面的海底滑坡海啸的快速预报方法的流程图。
具体实施方式
[0021]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将对本专利技术实施例中
的技术方案进行清楚、完整地描述。显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0022]下面结合实施例对本专利技术作进一步的描述。
[0023]如图1所示，本专利技术提供了一种基于响应面的海底滑坡海啸的快速预报方法，包括以下步骤：（1）根据研究资料，判断会发生海底滑坡的位置。
[0024]在本步骤中，需要说明的是，在判断会发生海底滑坡的位置时，通常海底滑坡易发生在富含天然气水合物的区域。根据专家对富含天然气水合物区域位置计算得到的滑坡安全系数分布云图，对安全系数小于1.1的区域绘制轮廓线，轮廓线内的区域即为发生海底滑坡的位置。
[0025]（2）获得包含滑坡位置的研究区域的高精度地形文件。
[0026]其中，在GEBCO网站中获得高精度地形文件。
[0027]在本步骤中，需要说明的是，GEBCO（General Bathymetric Chart of the Oceans）网站是由国际海道测量组织和政府间海洋学委员会协调有关国家联合编制的覆盖全球18幅1:1000万（两极为1:600万）的地势图。
[0028]（3）根据专家的研究得到海底滑坡所在的深度范围及滑坡体体积量的区间，并将深度和滑坡体体积量进行离散化，离散的区间个数分别为和，建议滑坡发生深度的离散区间为，滑坡体体积量离散区间为，则离散的区间个数分别为：；。
[0029]在本步骤中，需要说明的是，离散间隔视计算机算力而定，间隔越小，计算量则越庞大，但对滑坡海啸的预报则越准确。
[0030]（4）建立一个包含所有滑坡参数的目录。
[0031]其中，在目录中深度和滑坡体体积量的大小取值于各区间的中间数，海底滑坡的情景数量为。
[0032]（5）对个滑坡情景进行滑坡海啸的模拟，得到个研究范围内的海啸最大波高数据。
[0033]其中，采用Navier
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于响应面的海底滑坡海啸的快速预报方法，其特征在于，包括以下步骤：（1）判断会发生海底滑坡的位置；（2）获得包含滑坡位置的研究区域的高精度地形文件；（3）根据海底滑坡所在的深度范围及滑坡体体积量的区间，将深度和滑坡体体积量进行离散化，离散的区间个数分别为和；（4）建立一个包含所有滑坡参数的目录；（5）对个滑坡情景进行滑坡海啸的模拟，得到个研究范围内的海啸最大波高数据；（6）对研究范围内的所有像素点的最大波高进行统计，每个像素点受到深度和滑坡体体积量两个参数的影响，会产生个海啸的最大波高，对这三个变量的数据进行拟合；（7）构建出每个像素点的滑坡海啸最大波高响应面；（8）当有海底滑坡发生时，判断该滑坡发生的深度和滑坡体体积量是否在构建的响应面内；（9）获取整个区域的海啸最大波高，对海啸波高超过阈值的高风险区域进行快速预警，并对滑坡海啸的最大波高数据进行可视化处理。2.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法，其特征在于，在步骤（1）中，判断会发生海底滑坡的位置的具体方法为：根据对富含天然气水合物区域位置计算得到的滑坡安全系数分布云图，对安全系数小于设定阈值的区域绘制轮廓线，轮廓线内的区域即为会发生海底滑坡的位置。3.根据权利要求1所述的一种基于响应面的海底滑坡海啸的快速预报方法，其特征在于，在步骤（2）中，在GEBCO网站中获得高精度地形文件。4.根据权利要求1所述的一种基于响应面的海底滑坡...

【专利技术属性】
技术研发人员：白晔斐，刘金伟，周一帆，房欣，刘英男，魏笑然，支泓欢，
申请(专利权)人：海南浙江大学研究院，
类型：发明
国别省市：