[发明专利]利用深海海底压力数据进行海啸波侦测的方法

说明书

技术领域

本发明属于海洋监测技术和海洋防灾减灾技术领域，涉及海啸的监测预警预报，特别是利用在深海海底实时监测的压力数据进行海啸波检测识别的方法。

背景技术

我国以往的海啸预警系统是以陆上地震台网监测的海底地震数据为主要依据，应用海啸预警模式进行计算来产生预警信息的，并由沿岸验潮站数据提供验证。然而，海啸预警模式的开发是以历史海啸观测数据为基础，我国相对缺乏；模式参数需要实测数据加以修正；模式本质上是对海啸过程的模拟计算，也离不开实测验证；沿岸验潮站受位置所限，无法测量开阔海洋上的海啸波数据，仍只能应用于海啸发生后近岸数据的验证。以往的海啸预警存在的误报和漏报较多，代价非常高昂。

对于我国来说，地震海啸源主要位于马尼拉群岛、琉球群岛、日本列岛等位置，距离大陆很远。海啸波从其发源地传播到海岸边往往需要一段时间，如果在海啸波传播过程中，在离岸较远的深水海洋中，就能识别出海啸波并发出海啸预警，可为灾害响应提供足够时间，及时疏散民众、降低生命财产损失。

发明内容

本发明为解决公知技术中存在的技术问题而提供一种利用深海海底压力数据进行海啸波侦测的方法，采用该方法可以依据实时监测到的远岸深海海底压力数据及早准确地发现海啸波，从而获得更多的应急响应时间。

本发明为解决公知技术中存在的技术问题所采取的技术方案是：一种利用深海海底压力数据进行海啸波侦测的方法，采用以下步骤：

S1：采用压力传感器获取一段时间内的深海海底压力数据，从而获得相应时间段内的潮位值；

S2：采用插值多项式拟合天文潮波形，以时间为插值点横坐标，外推下一时刻t′的潮位值H′(t′)：

其中：a_i为插值系数，当插值点的间隔时间Δt和要推算的下一个潮位值时间t′确定时，插值系数a_i可由拉格朗日插值基函数计算得出；

为插值点的函数值，通过读取步骤S1中对应插值时间点的一段时间内的潮位值并采用算术平均法计算获得；

t是当前时间，也是步骤S1的结束时间；

p是算术平均法取值的时长；

Δt是两个相邻插值点的间隔时间；

S3：实测t′时刻深海海底压力数据，从而获得t′时刻的实测潮位值，从t′时刻的实测潮位值中分离由步骤S2计算得到的t′时刻潮位值H′(t′)，得到差值数据；

S4：依据连续的多个差值数据判断是否有海啸波传播。

本发明具有的优点和积极效果是：依据实测的海底压力数据，采用插值多项式拟合天文潮波形，以时间为插值点横坐标，外推下一时刻的潮位值，并将其视作天文潮波潮位，从实测的下一时刻的海底压力数据中分离出去，得到差值数据，依据连续的多个差值数据判别是否有海啸波在传播。本发明利用实测的、离岸较远的深海海底压力数据对海啸波进行分析和判别，能够及早准确地发现海啸波，从而提高海啸预报准确率，获得更多的应急响应时间。

附图说明

图1为本发明采用三次多项式拟合潮波时插值点的示意图。

具体实施方式

为能进一步了解本发明的发明内容、特点及功效，兹例举以下实施例，并配合附图详细说明如下：

一种利用深海海底压力数据进行海啸波侦测的方法，采用以下步骤：

S1：采用压力传感器获取一段时间内的深海海底压力数据，从而获得相应时间段内的潮位值。

实测的海底压力信号反映的是海平面高度变化，而这一变化的最主要表现形式是波。根据波在海洋中垂直方向的传播规律，在深海海底压力实测信号中，压力值波动的主要分量为海啸波和潮波，除此外，还存在一些由高频波和海洋动物活动等引起的环境杂波。

S2：采用插值多项式拟合天文潮波形，以时间为插值点横坐标，外推下一时刻t′的潮位值H′(t′)：

其中：a_i为插值系数，当插值点的间隔时间Δt和要推算的下一个潮位值时间t′确定时，插值系数a_i可由拉格朗日插值基函数计算得出；

为插值点的函数值，通过读取步骤S1中对应插值时间点的一段时间内的潮位值并采用算术平均法计算获得；采用算术平均法对一段时间内连续采集的多个压力数据取平均值，得到H^*，同时可以滤除海底压力波动中高于海啸波频率的高频杂波。

t是当前时间，也是步骤S1的结束时间；

p是算术平均法取值的时长；

Δt是两个相邻插值点的间隔时间；