CN106969755B - 一种基于水压力反演海浪波高与周期的方法 - Google Patents
一种基于水压力反演海浪波高与周期的方法 Download PDFInfo
- Publication number
- CN106969755B CN106969755B CN201710152430.9A CN201710152430A CN106969755B CN 106969755 B CN106969755 B CN 106969755B CN 201710152430 A CN201710152430 A CN 201710152430A CN 106969755 B CN106969755 B CN 106969755B
- Authority
- CN
- China
- Prior art keywords
- pressure
- wave
- water
- water pressure
- period
- 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
- G01C13/002—Measuring the movement of open water
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)
- Measuring Fluid Pressure (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于水压力反演海浪波高与周期的方法,包括:用水压式波高仪测量得到水下压力强度;得到水压力强度时间序列数据;通过对水压力强度时间序列数据进行傅里叶变换得到压力频谱;考虑水压力的线性趋势变化和压力波在深度上、水层过滤作用的衰减,计算出表面波浪的压力校正系数,对水压力进行校正;依据谱矩法或者对压力频谱进行傅里叶反变换利用上跨零点分波法,反演出高时空分辨率的海浪波高与周期。
Description
技术领域
本发明涉及海洋中海浪要素测量领域,具体涉及一种基于水压力反演海浪波高与周期的方法。
背景技术
海浪的波高与周期是海洋观测研究的重点内容之一,现场观测海浪的波高与周期过程最基本的就是要明确波浪的时空变化规律,为海洋波动规律提供基础认识。早期的海浪观测是对海浪施行目测,易产生较大误差、数据资料不连续、费时费力。之后采用数字记录,即以数字记录波面的高低,精度上欠缺。利用水压式波高仪的水压力强度信号反演海浪的波高与周期的测量技术具有较高精度,实施方便,并且能够得到高时空分辨率的波浪要素资料,为海洋波动观测奠定基础。
发明内容
有鉴于此,本发明提供了一种基于水压力反演海浪波高与周期的方法及其操作方法,目的在于提高海浪要素测量精度与避免资料处理复杂的问题,实现高时空分辨率海浪波高与周期的测量,详见下文描述。
一种基于水压力反演海浪波高与周期的方法,其特征在于,包括以下步骤:
1)将压力式波高仪放置在水下z处进行定点长期观测,通过测量水体压力变化得到压力强度信息。
2)对压力强度数据进行分段,去除其线性变化的趋势,得到水压力强度时间序列数据。
3)通过对水压力强度时间序列数据进行傅里叶变换得到压力频谱;
4)考虑水压力的线性趋势变化和压力波在深度上、水层过滤作用的衰减,计算出表面波浪的压力校正系数,对水压力进行校正;
5)依据谱矩法或者对压力频谱进行傅里叶反变换利用上跨零点分波法,反演出高时空分辨率的海浪波高与周期。
本发明实施例提供的技术方案的有益效果是:
本发明的基于水压力反演海浪波高与周期的方法,在接收到水压力强度和小振幅与波动振幅理论技术的基础上,综合考虑压力波的作用随着深度在衰减,水层的过滤作用是非线性且随频率变异,对压力频谱进行校正和水压力线性趋势变化的去除等过程,得到较为精确的海浪波高与周期。同时,此方法对于水体干扰小,具有较强的抗附着能力。并且可以定点连续长时间高频率测量水压力强度,最终实现高时空分辨率的海浪波高与周期的测量。
附图说明
图1为一种基于水压力反演海浪波高与周期的方法示意图;
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
本发明实施例介绍了一种基于水压力反演海浪波高与周期的方法,目的在于实现高时空分辨率海浪波高与周期的测量。具体实施步骤如下:
1、将压力式波高仪放置在水下z处(海底最佳)进行定点长期观测,通过测量水体压力变化得到压力强度信息。
2、对压力强度数据进行分段,去除其线性变化的趋势;比如2分钟一段,利用每段测量数据减去该段的平均值之后,即可去除其线性变化趋势。
3、压力频谱的计算。通过接收到的水压力强度时间序列数据进行傅里叶变换得到压力频谱。
4、对压力频谱进行校正,需要计算出表面波浪的压力校正系数。
依小振幅理论,水下z处压力随时间t的变化为:
其中,z为深度,波沿x轴传播,k为波数,σ为圆频率,a为自由表面上的振幅,h为水深。
表面波浪的波动振幅η随时间t的变化为:
η(x,t)=asin(kx-σt)
综合上述关系式,可得:
η(x,t)=p(x,z,t)/Kp(z)
其中,即为表面波浪的压力校正系数。下面的关键问题是如何求得kh,kh的计算方法如下:
其中,根据频散关系σ2=gktanh(kh),并且d1=0.6666,d2=0.3555,d3=0.1608465608,d4=0.0632098765,d5=0.0217540484,d6=0.0065407983。为防止校正系数对波谱在高频段的过度校正,需要对校正频率范围设置一个限制点,通过对压力数据和波浪的观测对比,得出限制频率为F=0.39-0.0094DHz。
5、利用压力校正系数校正压力频谱,能够换算成波谱表示,可根据谱矩法计算海浪的波高与周期。
6、利用校正后的压力频谱,通过傅里叶反变换将压力频谱转换成压力信号时间序列表达,根据上跨零点分波法计算出海浪的波高与周期。
7、整个反演方法可通过Matlab编程来实现,具体编程与操作步骤可按如下方法实现:
1)将观测的长时间序列数据进行分段,比如2分钟一段。
2)利用分段数据减去该段数据的平均值,去掉其线性变化的趋势,再按顺序合成新的时间序列数据。
3)对新的时间序列数据进行傅里叶变换,得到压力频谱信号。
4)对压力频谱进行校正,这里需要计算出表面波浪的压力校正系数。
5)校正后的压力频谱信号,可以根据谱矩法计算海浪的波高与周期。
6)或者从第4步开始,对校正后的压力频谱信号进行傅里叶反变换,将其转换成压力信号的时间序列。
7)根据上跨零点分波法,计算出海浪的波高与周期。
Claims (1)
1.一种基于水压力反演海浪波高与周期的方法,其特征在于,包括以下步骤:
1)将压力式波高仪放置在水下z处进行定点长期观测,通过测量水体压力变化得到压力强度信息;
2)对压力强度数据进行分段,去除其线性变化的趋势,得到水压力强度时间序列数据;
3)通过对水压力强度时间序列数据进行傅里叶变换得到压力频谱;
4)考虑水压力的线性趋势变化和压力波在深度上、水层过滤作用的衰减,计算出表面波浪的压力校正系数,对水压力进行校正,
所述的表面波浪的压力校正系数其中,z为深度,k为波数,a为自由表面上的振幅,h为水深,式中,其中,根据频散关系:σ2=gktanh(kh),而且d1=0.6666,d2=0.3555,d3=0.1608465608,d4=0.0632098765,d5=0.0217540484,d6=0.0065407983;
5)依据谱矩法或者对压力频谱进行傅里叶反变换利用上跨零点分波法,反演出高时空分辨率的海浪波高与周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710152430.9A CN106969755B (zh) | 2017-03-15 | 2017-03-15 | 一种基于水压力反演海浪波高与周期的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710152430.9A CN106969755B (zh) | 2017-03-15 | 2017-03-15 | 一种基于水压力反演海浪波高与周期的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106969755A CN106969755A (zh) | 2017-07-21 |
CN106969755B true CN106969755B (zh) | 2019-08-13 |
Family
ID=59328580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710152430.9A Expired - Fee Related CN106969755B (zh) | 2017-03-15 | 2017-03-15 | 一种基于水压力反演海浪波高与周期的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106969755B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860372B (zh) * | 2017-10-30 | 2020-02-07 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 海底摩擦导致波浪衰减的实地测量方法 |
CN110006409B (zh) * | 2019-01-17 | 2021-06-08 | 广东省水利水电科学研究院 | 一种用于智慧海堤感知系统的风暴潮波高采集设备 |
CN113155408B (zh) * | 2020-10-24 | 2023-04-18 | 大连理工大学 | 一种基于点压力反馈浪高显示的造波机实验装置及其方法 |
CN112683245B (zh) * | 2020-12-31 | 2021-11-30 | 广州海洋地质调查局 | 一种海洋孤立内波预警强度校正方法 |
CN113281757B (zh) * | 2021-04-22 | 2023-07-11 | 浙江大学 | 一种倒置式多波束回声仪及海面波浪测量方法 |
CN113408401B (zh) * | 2021-06-16 | 2022-02-22 | 中国科学院南海海洋研究所 | 一种基于机器学习的船行波快速自动识别方法及装置 |
CN115435757B (zh) * | 2022-11-09 | 2023-03-14 | 海南浙江大学研究院 | 一种利用水下动力压力反演表面波浪波高的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10141953A (ja) * | 1996-11-08 | 1998-05-29 | Kaijo Corp | 圧力式波高計 |
CN101706275A (zh) * | 2009-11-09 | 2010-05-12 | 哈尔滨工程大学 | 动态海浪测量装置及测量方法 |
CN101806594A (zh) * | 2010-04-23 | 2010-08-18 | 中国海洋石油总公司 | 一种基阵式波浪监测装置及其波浪测量方法 |
-
2017
- 2017-03-15 CN CN201710152430.9A patent/CN106969755B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10141953A (ja) * | 1996-11-08 | 1998-05-29 | Kaijo Corp | 圧力式波高計 |
CN101706275A (zh) * | 2009-11-09 | 2010-05-12 | 哈尔滨工程大学 | 动态海浪测量装置及测量方法 |
CN101806594A (zh) * | 2010-04-23 | 2010-08-18 | 中国海洋石油总公司 | 一种基阵式波浪监测装置及其波浪测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106969755A (zh) | 2017-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106969755B (zh) | 一种基于水压力反演海浪波高与周期的方法 | |
AU2008263602B2 (en) | Method of representing seismic signals | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡系统频率估计方法 | |
CN110135291B (zh) | 一种低信噪比信号的参数估计方法 | |
CN104833937A (zh) | 一种基于mir-rsd高精度余弦窗插值fft算法的谐波测量通道校准方法 | |
CN103454494A (zh) | 一种高精度的谐波分析方法 | |
CN106990402A (zh) | 一种基于波浪理论的导航x波段雷达波群检测方法 | |
CN103148790A (zh) | 一种锂电池薄膜微位移测厚方法 | |
CN111458678B (zh) | 一种基于时频干涉谱和辐射噪声声强测量的被动测距方法 | |
CN104143031B (zh) | 一种基于小波多尺度分解的植被指数时序数据重构方法 | |
CN113505479A (zh) | 密度反演方法、装置及电子设备 | |
Song et al. | Utilization of multimode surface wave dispersion for characterizing roadbed structure | |
CN109188542B (zh) | 一种波区相关性检测的远参考大地电磁阻抗计算方法 | |
CN113866836A (zh) | 一种基于归一化磁异常导数标准差的多目标边界识别方法 | |
CN114089416B (zh) | 一种利用薛定谔方程进行地震波衰减梯度估计的方法 | |
CN114280366A (zh) | 基于改进型频率插值算法的正弦信号频率估计方法 | |
CN109446552B (zh) | 多轴相关随机激励下结构疲劳寿命时域计算方法 | |
CN105372492B (zh) | 基于三条dft复数谱线的信号频率测量方法 | |
CN104914439A (zh) | 一种超声波测距的双相位测量方法 | |
CN107870034B (zh) | 一种基于相位差的水声声速测量方法 | |
CN108152363B (zh) | 一种抑端固有时间尺度分解的管道缺陷识别方法 | |
CN112505640B (zh) | 基于参数自适应的扩展b分布脉冲信号时频分析方法 | |
CN110297199B (zh) | 一种基于全相位fft的铯光泵磁力仪频率测量方法及系统 | |
RU2578126C1 (ru) | Способ формирования радиолокационных изображений | |
CN102928874B (zh) | 相对震级类比反演方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for 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 |
Granted publication date: 20190813 Termination date: 20210315 |
|
CF01 | Termination of patent right due to non-payment of annual fee |