舰船科学技术  2019, Vol. 41 Issue (8): 131-134   PDF    
基于船载北斗导航仪测量波浪参数方法研究
周红进1, 李方能2, 黄谦1, 李伟1, 兰国辉1     
1. 海军大连舰艇学院 航海系,辽宁 大连 116001;
2. 海军工程大学 导航工程系,湖北 武汉 430032
摘要: 获取舰船波浪参数对于舰船的航行安全和武器使用十分重要。针对舰员缺乏有效手段实时获取舰船当前航行局部海域的波浪信息的问题,提出利用船载北斗接收机输出的天线高度信息反演波浪场信息的新方法。该方法认为北斗测量天线相对海平面高度的信息变化,反映了波浪的运动状态,通过对天线高度信息进行谱分析,根据波浪频谱特征,进行带通滤波,从而在频域提取出波浪信息,经过傅里叶反变换后,得到时域波浪信息,然后计算波高和波周期。海上实船试验结果表明,与浮标测量的准确波浪参数相比,通过该方法测量得到的SWH波高误差小于6%,波周期误差小于10%。
关键词: 波浪     北斗卫星导航系统     谱分析     带通滤波     傅里叶变换    
Research on wave parameters measurements via vessel mounted compass navigator
ZHOU Hong-jin1, LI Fan-neng2, HUANG Qian1, LI Wei1, LAN Guo-hui1     
1. Department of Marine Navigation, Dalian Naval Academy, Dalian 116001, China;
2. Department of Navigation, Naval University of Engineering, Wuhan 430032, China
Abstract: It's very important to get wave parameters information for vessels' safety. But lack of vessels mounted equipments to measure wave parameters in real time and real spot make it difficult for crews. One new method is presented to solve the problem by dealing with vessel mounted compass navigator’s output data. This method deems that the change of altitude of antenna reflects the wave movements above mean sea level. First, to make a spectrum analysis of antenna altitude data, second, to design a band-pass filter according to the wave frequency-spectrum character to extract the signal in frequency-domain, third, to extract the data in time-domain via inverse Fourier transformation, at last, the wave parameters such as SWH and period can be resolved. The test results on a ship shows the methods is feasible, its SWH determination accuracy is better than 6%, and the period determination accuracy is better than 10%.
Key words: wave     compass satellite navigation system     spectrum analysis     band-pass filter     fourier transformation    
0 引 言

波浪参数(波高、波周期等)对于舰船航行安全、拟定航行计划至关重要。现代舰船获取波浪参数信息的途径主要通过接收气象传真图或卫星云图获取实况分析和短期预报,但这种方法获取的波浪参数尺度一般都很大,是关于整个海区的波浪参数信息,舰船自身所处海况则无从获知。导致这一现状的原因是缺乏船载波浪参数测量设备。

船载测量波高对于远洋运输和军事行动意义重大,但一直都没有很好解决,主要有两方面的原因:一是传统的波高测量设备大多是接触式波高计或浪潮仪,如浮标式波高计、压力传感式波高计/浪潮仪、超声波波高计等等,它们或浮于海面,或沉于海底,不能跟随船体移动,因此难以用于船载测量;是船体在海洋中受到海浪的影响而产生颠簸,测量海浪波高时没有静止参照物,因此不能正确地测量海浪波高[15]

北斗卫星导航系统是我国独立自主发展和运行的卫星导航系统。现代舰船,尤其是军舰普遍安装了北斗卫星导航仪。北斗导航仪能够实时测量天线的高度(相对海平面的高度),在短时间的航行中(如10-20 min),舰船所处局部位置的海平面高度是不变的。但舰船会随着波浪升沉,天线的实际高度也随之变化。通过对北斗导航仪输出的天线高度数据进行频谱分析,结合波浪频谱特征,进行带通滤波,在频域内提取波浪信号,然后进行傅里叶逆变换,反演得到时域波浪数据,从而计算波高、波周期等参数。

1 测量方法 1.1 波高测量统计

波高与周期的测量包括最大波高、1/10最大波高、1/3最大波高、平均周期等特征值[6]

不规则波的波形中具有一定统计特征的波高。由于外观波形的不规则性和随机性,目前多数采用“上跨零点法”对外观波形(表视波)及其要素给予定义。对一个不规则波的波高仪记录,首先确定其无趋势项的水位零点,波动水位过程线自波谷上升与零线交点称为上跨零点,自波峰下降与零线交点称为下跨零点,2个相邻上跨零点之间的波形称为一个表视波,2个相邻上跨零点的跨离称为周期,波峰顶和波谷底的垂直距离称为波高。

最大波高是波浪连续记录中波高的最大值。将波浪之波高大小按照顺序排列,所有波高个数中最大之1/10个波高的平均值,称为1/10最大波高,记为 ${H_{1/10}}$ ;所有波高个数中最大之1/3个波高的平均值,称为1/3最大波高,或称为有意义波高、有效波高,记为 $SWH$ 。平均波高是波浪连续记录中所有波高的平均值,其对应周期的平均值为平均周期。

1.2 波浪信息提取

本文利用北斗导航仪输出的天线高度信息研究海浪波形参数。安装在舰船上的北斗导航仪输出信息包括:北斗天线随舰船运动的速度信息,实时反映舰船对地航速;北斗天线的位置信息,实时反映舰船的位置;北斗天线相对海平面的高度信息,既反映了舰船所处位置的海平面高度,也实时反映了舰船随波浪升沉的信息[710]

速度信息虽然也包含了波浪波速信息,但主要包含了舰船在螺旋桨推动下的对地航速,难以从中提取波浪信息。位置信息通常是经纬度数据,主要包含了舰船的平面坐标,由于舰船是主动力推进运动的,不是浮标,也难以从中提取波浪信息。

在10~20 min,舰船位置变化小,海平面高度基本不变,此时天线高度信息则直接反映了舰船随波浪升沉的信息。因此适合提取波浪信息。

对北斗导航仪输出的天线高度序列进行傅离叶变换获得高度谱。自相关函数定义为:

${R_{ii}}(\tau ) = \mathop {\lim }\limits_{T \to \infty } \frac{1}{T}\int^{T/2}_{ - T/2} {{S_i}(t){S_i}(t + \tau )} {\rm{d}}t\text{,}$ (1)

其中, ${S_i}(t)$ 为北斗导航仪输出的天线高度序列。对式进行傅里叶变换得到:

$S(f) = \int_{ - \infty }^\infty {{R_{ii}}(t){e^{( - i2{\text{π}} ft)}}} {\rm{d}}t\text{。}$ (2)

$S(f)$ 为高度谱,表示天线高度在频域中某个频段的能量分布。由于波浪的运动位于一个特定的频带内,通过带通滤波可以获得波浪的频谱。滤波方程为:

${S_W}(f) = S(f) \cdot H(f)\text{,}$ (3)

其中: ${S_W}(f)$ 为波浪运动某个频段内的谱, $H(f)$ 为带通滤波函数。

$H(f) = \left\{ {\begin{array}{*{20}{l}} {0\;\;f < A,f > B}\text{,} \\ {1\;\;A \leqslant f \leqslant B}\text{。} \end{array}} \right.$ (4)

$A = 0.03,B = 0.4$ ,比较吻合波浪运动谱性质。

${S_W}(f)$ 进行逆傅里叶变换:

${S_W}(t) = \frac{1}{{2{\text{π}} }}\int_{ - \infty }^\infty {{S_W}(\omega ){e^{(j\omega t)}}} {\rm{d}}\omega \text{。}$ (5)

其中, ${S_W}(t)$ 为时域波浪信号表示。

${S_W}(t)$ 提取前1/10波高最大值,记为 ${S_W}{(n)_{1/10}}$ ,从 ${S_W}(t)$ 提取前1/3波高最大值,记为 ${S_W}{(n)_{1/3}}$ ,根据2.1节的定义有:

${H_{1/10}} = \frac{1}{N}\sum\limits_{n = 1}^N {{S_W}{{(n)}_{1/10}}} \text{,}$ (6)
$SWH = \frac{1}{N}\sum\limits_{n = 1}^N {{S_W}{{(n)}_{1/3}}} \text{,}$ (7)
$T = \frac{{{f_s} \cdot {N_W}}}{N}\text{。}$ (8)

其中: $T$ 为波周期; ${f_s}$ 为采样频率; $N$ 为采样点数。

2 实验及结果

某船在东海海域执行任务期间,采集记录了北斗数据,图1为北斗天线安装位置。通过国家海洋局水文气象中心,获得了当时船所在位置的准确波浪信息,将该波浪信息作为基准衡量本文提出的测量波高方法的可行性和准确度。

图 1 北斗天线安装位置 Fig. 1 Beidou antenna mounting position

利用本文测量波浪的方法,对采集的4个位置和时间段的数据分别进行了处理,得到波高和波周期,图2图9分别为同时间4个位置的波浪高度谱和波形参数曲线。

图 2 波浪谱(20140912-0600) Fig. 2 Wave spectrum (20140912-0600)

图 3 波形(20140912-0600) Fig. 3 Wave height (20140912-0600)

图 4 波浪谱(20140912-1200) Fig. 4 Wave spectrum (20140912-1200)

图 5 波形(20140912-1200) Fig. 5 Wave height (20140912-1200)

图 6 波浪谱(20140915-0000) Fig. 6 Wave spectrum (20140915-0000)

图 7 波形(20140915-0000) Fig. 7 Wave height (20140915-0000)

图 8 波浪谱(20140918-0600) Fig. 8 Wave spectrum (20140918-0600)

图 9 波形(20140918-0600) Fig. 9 Wave height (20140918-0600)

与国家海洋局水文气象中心实测数据比较结果如表1所示(水温气象中心波高为1/3波高,即Significant Wave Height-SWH)。

表 1 试验结果 Tab.1 Test results

从数据处理结果来看,利用本文提出的方法,波高测量的精度优于6.0%,总体而言,波浪较大的时候,测量精度要把波浪小的时候精度要高。从数据处理过程来看,每个位置点的数据持续时间大约为15 min,短的时候只有5 min,比如2014年9月12日1 200的数据,由于北斗发生了信号丢失,有效数据只有5 min。对于一个波周期而言,5 min的数据也足够,在进行数据处理时,尤其是最后一步针对提取处理的波浪信号进行逆傅里叶变换后,计算波高和波周期时,由于滤波稳定时间的存在,通常应在1~2波周期后,开始计算波高,这样比较准确。

3 结 语

利用船载北斗导航仪测量的天线高度信息,通过谱分析,结合波浪的频谱分布,进行带通滤波,可以从天线高度信息提取波浪信号,将频域信号通过逆傅里叶变换,就可以获取波浪在时域内的波形变化规律,从而计算出波高和波周期。实测试验结果表明,该方法可以准确的测量波浪参数,波高测量误差可以达到6.0%以内,波周期测量误差可以达到10%以内。这种方法为远航舰船或执行军事任务的军舰提供了一种实时获取舰船所在位置的实时波浪信息的有效手段。

在使用本文提出的方法提取波浪信息时,有3点需要注意:

1)需要近似获取波浪的频谱特征。可以从气象传真图或卫星云图的大尺度海域波浪分析图中获取。

2)本文的方法目前还只是经过了4批海上实测数据的检验,其普遍适用性还有待进一步的验证。

3)波浪信息参数不仅包括波高、波周期,还包括波向分布。获取波向谱对于舰船操纵具有重要意义。这也是下一步的研究方向。

参考文献
[1]
YE Xiao-mm, LIN Mings-en, XU Ying. Validation of Chinese HY-2 satellite radar ahimeter significant wave height[J]. Acta Oceanologica Sinica, 2015, 34(5): 60-67. DOI:10.1007/s13131-015-0667-y
[2]
PARK K A, WOO H J, I. EE E Y, et al. Validation of significant wave height from satellite ahimeter in the seas around Korea and Error Charactensncs[J]. Korean Journal of Remote Sensing, 2013, 29(6): 631-644. DOI:10.7780/kjrs.2013.29.6.6
[3]
SHANAS P R, SANIL KUMAR V, HITHIN N K. Companson of gridded multi-mission and along-track mono-mission satellite ahimetry wave heights with in sun nearshore buoy data[J]. Ocean Engineering, 2014, 83: 24-35. DOI:10.1016/j.oceaneng.2014.03.014
[4]
韩伟孝, 杨俊钢, 王际朝. 基于浮标数据的卫星雷达高度计海浪波高数据评价与校正[J]. 海洋学报, 2016, Vol.38(11): 73-89. DOI:10.3969/j.issn.0253-4193.2016.11.007
[5]
JANSSEN P A E M, ABDALLA S, HERSBACH H, et al. Error estimation of buoy, satellite, and model wave height data[J]. Journal of Atmospheric and Technology, 2007, 24(9): 1665-1677. DOI:10.1175/JTECH2069.1
[6]
李静. 船模通航波高采集测量系统及其应用[J]. 水科学与工程技术, 2007(1): 52-55. DOI:10.3969/j.issn.1672-9900.2007.01.020
[7]
尹聪, 曹云昌, 朱彬, 等. GNSS-R信号测量海面有效波高的应用[J]. 热带海洋学报, 2012, Vol.31(6): 36-40. DOI:10.3969/j.issn.1009-5470.2012.06.006
[8]
李颖, 朱雪瑗, 崔璨, 等. 船载GNSS-R有效波高测量的初步研究[J]. 海洋环境科学, 2016, 35(2): 180-183.
[9]
沈继红, 宋静怡, 李众, 戴运桃, 王淑娟. 导航雷达反演海面波高的方法研究[J]. 中国航海, 2013, 36(1): 9-13. DOI:10.3760/cma.j.issn.1009-6906.2013.01.004
[10]
朱雪松, 严硕华, 张训械, 等. 基于导航卫星载波信号的海浪高度反演系统[J]. 电波科学学报, 2012, 27(6): 1070-1075.