CN105136126B - 利用深海海底压力数据进行海啸波侦测的方法 - Google Patents
利用深海海底压力数据进行海啸波侦测的方法 Download PDFInfo
- Publication number
- CN105136126B CN105136126B CN201510535091.3A CN201510535091A CN105136126B CN 105136126 B CN105136126 B CN 105136126B CN 201510535091 A CN201510535091 A CN 201510535091A CN 105136126 B CN105136126 B CN 105136126B
- Authority
- CN
- China
- Prior art keywords
- level value
- time
- tidal level
- tsunami
- moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用深海海底压力数据进行海啸波侦测的方法,采用以下步骤:S1:采用压力传感器获取一段时间内的深海海底压力数据,从而获得相应时间段内的潮位值;S2:采用插值多项式拟合天文潮波形,以时间为插值点横坐标,外推下一时刻t′的潮位值H′(t′);S3:实测t′时刻深海海底压力数据,从而获得t′时刻的实测潮位值,从t′时刻的实测潮位值中分离由步骤S2计算得到的t′时刻潮位值H′(t′),得到差值,依据连续的多个差值数据判断是否有海啸波传播。本发明利用实测的、离岸较远的深海海底压力数据对海啸波进行分析和判别,能够及早准确地发现海啸波,从而提高海啸预报准确率,获得更多的应急响应时间。
Description
技术领域
本发明属于海洋监测技术和海洋防灾减灾技术领域,涉及海啸的监测预警预报,特别是利用在深海海底实时监测的压力数据进行海啸波检测识别的方法。
背景技术
我国以往的海啸预警系统是以陆上地震台网监测的海底地震数据为主要依据,应用海啸预警模式进行计算来产生预警信息的,并由沿岸验潮站数据提供验证。然而,海啸预警模式的开发是以历史海啸观测数据为基础,我国相对缺乏;模式参数需要实测数据加以修正;模式本质上是对海啸过程的模拟计算,也离不开实测验证;沿岸验潮站受位置所限,无法测量开阔海洋上的海啸波数据,仍只能应用于海啸发生后近岸数据的验证。以往的海啸预警存在的误报和漏报较多,代价非常高昂。
对于我国来说,地震海啸源主要位于马尼拉群岛、琉球群岛、日本列岛等位置,距离大陆很远。海啸波从其发源地传播到海岸边往往需要一段时间,如果在海啸波传播过程中,在离岸较远的深水海洋中,就能识别出海啸波并发出海啸预警,可为灾害响应提供足够时间,及时疏散民众、降低生命财产损失。
发明内容
本发明为解决公知技术中存在的技术问题而提供一种利用深海海底压力数据进行海啸波侦测的方法,采用该方法可以依据实时监测到的远岸深海海底压力数据及早准确地发现海啸波,从而获得更多的应急响应时间。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种利用深海海底压力数据进行海啸波侦测的方法,采用以下步骤:
S1:采用压力传感器获取一段时间内的深海海底压力数据,从而获得相应时间段内的潮位值;
S2:采用插值多项式拟合天文潮波形,以时间为插值点横坐标,外推下一时刻t′的潮位值H′(t′):
其中:ai为插值系数,当插值点的间隔时间Δt和要推算的下一个潮位值时间t′确定时,插值系数ai可由拉格朗日插值基函数计算得出;
为插值点的函数值,通过读取步骤S1中对应插值时间点的一段时间内的潮位值并采用算术平均法计算获得;
t是当前时间,也是步骤S1的结束时间;
p是算术平均法取值的时长;
Δt是两个相邻插值点的间隔时间;
S3:实测t′时刻深海海底压力数据,从而获得t′时刻的实测潮位值,从t′时刻的实测潮位值中分离由步骤S2计算得到的t′时刻潮位值H′(t′),得到差值数据;
S4:依据连续的多个差值数据判断是否有海啸波传播。
本发明具有的优点和积极效果是:依据实测的海底压力数据,采用插值多项式拟合天文潮波形,以时间为插值点横坐标,外推下一时刻的潮位值,并将其视作天文潮波潮位,从实测的下一时刻的海底压力数据中分离出去,得到差值数据,依据连续的多个差值数据判别是否有海啸波在传播。本发明利用实测的、离岸较远的深海海底压力数据对海啸波进行分析和判别,能够及早准确地发现海啸波,从而提高海啸预报准确率,获得更多的应急响应时间。
附图说明
图1为本发明采用三次多项式拟合潮波时插值点的示意图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
一种利用深海海底压力数据进行海啸波侦测的方法,采用以下步骤:
S1:采用压力传感器获取一段时间内的深海海底压力数据,从而获得相应时间段内的潮位值。
实测的海底压力信号反映的是海平面高度变化,而这一变化的最主要表现形式是波。根据波在海洋中垂直方向的传播规律,在深海海底压力实测信号中,压力值波动的主要分量为海啸波和潮波,除此外,还存在一些由高频波和海洋动物活动等引起的环境杂波。
S2:采用插值多项式拟合天文潮波形,以时间为插值点横坐标,外推下一时刻t′的潮位值H′(t′):
其中:ai为插值系数,当插值点的间隔时间Δt和要推算的下一个潮位值时间t′确定时,插值系数ai可由拉格朗日插值基函数计算得出;
为插值点的函数值,通过读取步骤S1中对应插值时间点的一段时间内的潮位值并采用算术平均法计算获得;采用算术平均法对一段时间内连续采集的多个压力数据取平均值,得到H*,同时可以滤除海底压力波动中高于海啸波频率的高频杂波。
t是当前时间,也是步骤S1的结束时间;
p是算术平均法取值的时长;
Δt是两个相邻插值点的间隔时间;
下一个潮位值对应的时间t′由压力数据的采样频率确定。
天文潮的波形为正弦波,可用多项式来拟合,多项式的阶数n是综合考虑曲线拟合精度和硬件计算性能来确定。以三阶插值多项式拟合天文潮为例,请参见图1,在步骤S1中,设定监测的结束时间为当前时刻t=0,将当前时刻t前、Δt时间前、2Δt时间前、3Δt时间前的前p分钟内的压力数据进行平均,得到四个潮位值则插值点已知函数值对应的插值时间点为(-p/2)、(-p/2-Δt)、(-p/2-2Δt)、(-p/2-3Δt),下一时刻潮位值H′(t′)外推计算公式即可写为:
S3:实测t′时刻深海海底压力数据,从而获得t′时刻的实测潮位值,从t′时刻的实测潮位值中分离由步骤S2计算得到的t′时刻潮位值H′(t′),得到差值数据;
将步骤S2推算出的下一时刻潮位值H′(t′),视作天文潮波潮位;从实测压力数据中分离天文潮波潮位数据,得到的剩余信号为海平面波动的幅度,其大小和变化特征可用于判断是否有海啸波。
S4:依据连续的多个差值数据判断是否有海啸波传播。
实测t′时刻的深海海底压力数据和由步骤S2计算得到的t′时刻潮位值H′(t′)按照压力传感器的采样间隔不断更新,依据连续的多个差值数据就能够准确判断是否有海啸波。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种利用深海海底压力数据进行海啸波侦测的方法,其特征在于,采用以下步骤:
S1:采用压力传感器获取一段时间内的深海海底压力数据,从而获得相应时间段内的潮位值;
S2:采用插值多项式拟合天文潮波形,以时间为插值点横坐标,外推下一时刻t′的潮位值H′(t′):
<mrow>
<msup>
<mi>H</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>t</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>H</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>p</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mi>i</mi>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
其中:ai为插值系数,当插值点的间隔时间Δt和要推算的下一个潮位值时间t′确定时,插值系数ai可由拉格朗日插值基函数计算得出;
为插值点的函数值,通过读取步骤S1中对应插值时间点的一段时间内的潮位值并采用算术平均法计算获得;
t是当前时间,也是步骤S1的结束时间;
p是算术平均法取值的时长;
Δt是两个相邻插值点的间隔时间;
S3:实测t′时刻深海海底压力数据,从而获得t′时刻的实测潮位值,从t′时刻的实测潮位值中分离由步骤S2计算得到的t′时刻潮位值H′(t′),得到差值数据;
S4:依据连续的多个差值数据判断是否有海啸波传播。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510535091.3A CN105136126B (zh) | 2015-08-27 | 2015-08-27 | 利用深海海底压力数据进行海啸波侦测的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510535091.3A CN105136126B (zh) | 2015-08-27 | 2015-08-27 | 利用深海海底压力数据进行海啸波侦测的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105136126A CN105136126A (zh) | 2015-12-09 |
CN105136126B true CN105136126B (zh) | 2017-10-10 |
Family
ID=54721566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510535091.3A Expired - Fee Related CN105136126B (zh) | 2015-08-27 | 2015-08-27 | 利用深海海底压力数据进行海啸波侦测的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105136126B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6858415B2 (ja) * | 2019-01-11 | 2021-04-14 | 学校法人福岡工業大学 | 海面計測システム、海面計測方法および海面計測プログラム |
CN111323809B (zh) * | 2020-03-17 | 2021-09-28 | 河海大学 | 一种监测海底地震引起海啸的装置及方法 |
CN112233387A (zh) * | 2020-10-12 | 2021-01-15 | 中国海洋大学 | 一种沿海风暴潮监测装置及在线监测预警系统 |
CN114572347A (zh) * | 2022-03-23 | 2022-06-03 | 国家海洋技术中心 | 一种海啸预警监测系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288147A (zh) * | 2000-11-01 | 2001-03-21 | 华东师范大学 | 微波潮位传感器及其用途 |
CN101441078A (zh) * | 2008-12-25 | 2009-05-27 | 杭州电子科技大学 | 江河涌潮分段实时预警方法 |
CN102221389A (zh) * | 2011-04-11 | 2011-10-19 | 国家海洋信息中心 | 结合统计模型与动力模型的乘潮水位预报方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3471701B2 (ja) * | 2000-03-15 | 2003-12-02 | エヌイーシーテレネットワークス株式会社 | 海底津波計システム、海底津波計装置及びその方法 |
JP2012058062A (ja) * | 2010-09-08 | 2012-03-22 | Nippon Telegr & Teleph Corp <Ntt> | 津波規模予測装置、方法、およびプログラム |
-
2015
- 2015-08-27 CN CN201510535091.3A patent/CN105136126B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1288147A (zh) * | 2000-11-01 | 2001-03-21 | 华东师范大学 | 微波潮位传感器及其用途 |
CN101441078A (zh) * | 2008-12-25 | 2009-05-27 | 杭州电子科技大学 | 江河涌潮分段实时预警方法 |
CN102221389A (zh) * | 2011-04-11 | 2011-10-19 | 国家海洋信息中心 | 结合统计模型与动力模型的乘潮水位预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105136126A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105136126B (zh) | 利用深海海底压力数据进行海啸波侦测的方法 | |
CN104180873A (zh) | 一种单波束测深仪水深粗差检测修正方法及系统 | |
CN102901560A (zh) | 一种海洋导管架平台结构安全综合监测系统 | |
Cheng et al. | Field measurements of inhomogeneous wave conditions in Bjørnafjorden | |
CN110765686B (zh) | 利用有限波段海底地形进行船载声呐测深测线设计的方法 | |
CN104075734A (zh) | 水下组合导航故障智能诊断方法 | |
CN105651265A (zh) | 基于波浪压力的跨海桥梁海域的波浪要素及潮位测量方法 | |
Wang et al. | An automated procedure to calculate the morphological parameters of superimposed rhythmic bedforms | |
Christou et al. | Examining a comprehensive dataset containing thousands of freak wave events: Part 1—description of the data and quality control procedure | |
CN102087107B (zh) | 绳系多传感器协同优化近海测波浮标及其滤波融合方法 | |
CN112632868B (zh) | 高频地波雷达观测径向流缺失值的填补修正方法及系统 | |
CN117930206A (zh) | 用于水面无人艇的海底地形测量系统及测量数据校正方法 | |
Bressan et al. | Detecting the 11 March 2011 Tohoku tsunami arrival on sea-level records in the Pacific Ocean: application and performance of the Tsunami Early Detection Algorithm (TEDA) | |
CN116738375A (zh) | 基于单条带测深数据的诱导升沉误差探测消除方法及系统 | |
Wang et al. | Detection and Elimination of Bathymetric Outliers in Multibeam Echosounder System Based on Robust Multi-quadric Method and Median Parameter Model. | |
Zhao et al. | Detection method for submarine oil pipeline leakage under complex sea conditions by unmanned underwater vehicle | |
CN106125100A (zh) | 一种空间信号连续性评估方法 | |
CN116499532B (zh) | 基于水文模型构建的复杂海洋环境深水群桩施工监测系统 | |
Pouliquen | Recommendations for in-situ data Real Time Quality Control | |
Bell | Determination of bathymetry using marine radar images of waves | |
CN107064875B (zh) | 基于一步状态估计和协方差估计的距离野值剔除方法 | |
Cao et al. | Travel time processing for LBL positioning system | |
Goring | Detecting datum changes using tide gauge records | |
CN103278851A (zh) | 一种海底管道掩埋状况检测方法 | |
Taylor et al. | Maximising data return: Towards a quality control strategy for managing and processing TRDI ADCP data sets from moored instrumentation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171010 Termination date: 20200827 |