机载波谱仪斑点噪声去除方法对比分析

针对斑点噪声降低机载波谱仪海浪谱反演精度问题,基于STORM机载波谱仪VALPARESO试验数据,分别使用数学法、阈值法、交叉谱法和小波分析法四种斑点噪声去除方法对其进行去噪处理,得到了去噪前后调制密度谱、海浪谱和噪声谱,并将反演所得海浪谱有效波高与同步观测的浮标结果进行了对比,最后对四种方法的去噪效果以及噪声谱特性进行了分析与评价。结果表明:四种去噪方法均可以有效去除斑点噪声,按照反演精度排序依次为交叉谱法、阈值法、小波分析法和数学法,其中交叉谱方法去噪效果最理想,RMSE仅为0.008,有效波高绝对误差为0.25m,相对误差为7%。     

波谱仪调制功率谱估计算法性能的仿真研究

基于星载波谱仪海浪方向谱探测原理，仿真了不同海况、风速下的海浪波谱仪接收信号，并采用周期图法、Welch法、AR模型法以及最小方差法共4种不同的调制谱估计方法反演出海浪谱，比较各种调制谱估计方法的海浪方向谱反演性能。仿真结果表明:对于一定方位向下的一维海浪谱反演，不同调制谱估计方法反演海浪谱性能优劣没有绝对的顺序。对于二维海浪谱反演，在成长中海浪条件下周期图法反演性能最差，其他3种方法的反演性能没有绝对优劣顺序；对于成熟风浪，最小方差法在积分能量误差、有效波高误差两个指标上的反演性能最好，在主波波向、主波波长误差上，周期图法反演性能最差，其他3种方法没有绝对优劣顺序。在涌浪条件下，AR模型法反演性能优于其他3种方法。在不同海况下，随着波谱仪中心入射角的下降，反演性能会下降。基于这些仿真结果，本文推荐最小方差法作为充分成长海浪的海浪方向谱反演的调制谱估计方法，AR模型法作为涌浪海浪方向谱反演的调制谱估计方法。     

基于同步风场的机载SAR海浪参数反演方法

针对机载合成孔径雷达(SAR)对海探测特点,采用多入射角法从SAR数据本身得到与海浪参数反演区域时空匹配的同步海面风速和风向,并结合线性变换关系,计算得到海浪初猜谱对应的仿真SAR图像谱,将仿真SAR图像谱和观测SAR图像谱输入代价函数中进行迭代运算,通过非线性方程的解算得到最适海浪谱;采用交叉谱法去除海浪传播180°方向模糊,最终得到海浪参数。论文提出的基于同步风场的机载SAR海浪参数反演方法,充分利用了机载SAR海洋环境探测的优势,解决了传统SAR海浪参数反演中初猜谱构造依赖外部风场的问题,机载同步飞行试验的海浪参数反演结果与浮标观测值的有效波高、波向的均方根误差分别为0.23 m和13.23°,验证了该方法的有效性,可为机载SAR海浪参数反演业务化提供支持。     

GPS浮标数据反演海浪谱的理论仿真与试验验证

GPS浮标作为一种新型的海洋测量设备,近年来在海面高度现场测量和星载高度计定标方面取得了重要应用。通过仿真试验对反演海浪谱的方法和流程进行研究,旨在探索从GPS浮标测量的海面高度序列中提取海浪谱的方法。首先,使用Longuest模型生成了海浪波面位移时间序列,并通过Pierson-Moscowitz风浪谱对波面位移的统计特性进行约束,其随机性由相位引入。结合典型潮汐和GPS浮标仪器噪声的仿真时间序列,合成了仿真时间长度1h的1Hz(每秒1次)随机海面高度序列。然后,利用自相关函数法,进行高通滤波和数据压缩,得到了仿真的海浪谱。该仿真结果和理论海浪谱非常接近,可满足海浪谱反演的需求。最后,通过山东石岛外海的GPS浮标现场试验,验证了本文提出的反演方法的适用性。本文的研究解决了GPS浮标反演海浪谱的关键问题,丰富了海浪谱反演的手段,拓展了GPS浮标的应用领域,有望为未来我国的星载波谱仪定标服务。     

基于GPS高频数据的海浪测量方法

研究了基于GPS高频数据进行海浪测量的方法。利用TRACK解算得海上载体高精度的垂向位移,经浪潮分离提取海浪信号,分别采用相关函数法和周期图法估计海浪信号的功率谱,并计算海浪要素。利用实测数据进行试验,结果表明,周期图法推算的平均波高和平均周期精度较高,与测波仪结果差异分别小于2cm和0.25s,基于GPS高频数据的海浪测量方法可有效反演海浪要素。     

基于SAR海浪方向谱的中国海海浪特性研究

基于Hasselmanns提出的SAR反演海浪方向谱的MPI算法,对2003-2012年间的中国海Envisat ASAR波模式数据进行了海浪方向谱反演。统计由反演的海浪方向谱得到的海浪有效波高数据,依据海浪浪级的划分,分析了中国海海浪浪级的分布特点,获得一些有参考价值的结果:(1)中国海3-4级海况的年出现概率最高,达85%;(2)中国海累月不同浪级的海浪出现概率分布符合高斯分布函数:f(x)=a*exp(-((x-1-b)/c)^2);(3)分析上述高斯分布函数的拟合系数a、b、c,发现其分布也有一定规律性。同时应用反演得到的海浪有效波高、平均波周期、平均波向等数据,分析了中国海的海浪时空分布特性,得到一些可供参考的结果。     

星载高级合成孔径雷达波模式算法及其反演数据精度验证

搭载在欧洲环境卫星(ENVISAT)上的高级合成孔径雷达(Advanced Synthetic Aperture Radar,ASAR)二级波模式数据提供了诸多海浪信息包括有效波高、波向、波长和二维海浪谱等,在海浪预报模式中具有重要作用。本文拟利用浮标观测数据对ASAR波模式算法及其反演数据精度进行对比验证。由于SAR卫星在海面的特殊成像机制,不同海况下会有不同的测量结果,通过与美国国家浮标中心(NDBC)的浮标数据对比,显示ASAR有效波高在高海况下低估和在低海况下高估的现象,在中等海况下的测量结果较优。通过研究ASAR数据集中对应的海浪谱,按照能量与方向分布可分为四种类型:单一方向海浪谱(Ⅰ类谱),180°方向模糊海浪谱(Ⅱ类谱),海浪两个方向且能量分布杂乱(Ⅲ类谱),多个传播方向且谱型杂乱海浪谱(Ⅳ类谱)。探究在不同类型下的海浪参数的精度,结果表明在单一波向正常海浪谱情况下,有效波高、波向与浮标数据一致性较好,存在180°方向模糊的对称海浪谱仅有效波高精度较高,谱型杂乱的海浪谱海浪有效波高和波向反演结果均较差。     

X-波段雷达近海海浪频谱反演的神经网络模型

X-波段雷达作为国内海浪观测的一种新工具,在海浪频谱获取和有效波高反演方面仍存在较多问题.本文利用非线性回归方法,将现场实测浮标数据频谱和雷达一维图像谱分别与标准频谱模型进行拟合,发现浮标频谱和一维图像谱具有标准频谱的特征,能够较准确地获取相应的谱参数.提出了建立由雷达一维图像谱参数反演海浪频谱参数的神经网络模型,同时在模型中加入影像序列信噪比,进而反演有效波高,并将反演结果与现场实测数据和传统算法(建立影像序列信噪比与有效波高之间的线性回归方程)进行了对比,结果表明,获取谱参数的误差和反演有效波高的平均误差在20％以内,而传统算法计算有效波高平均误差在20％以上.     

基于SAR图像速度聚束调制的海浪反演研究

本文首先对合成孔径雷达（SAR）海浪成像中的3种调制（倾斜调制、流体力学调制与速度聚束调制）的影响进行了对比分析,结果显示:速度聚束调制对SAR图像的影响最为显著。另外,由于SAR图像中固有相干斑噪声的存在,较低波数范围的噪声难以滤除或抑制,利用经典MPI方法反演海浪谱会造成低波数范围谱值偏大。基于此,本文借鉴经典MPI海浪谱反演算法,建立了基于速度聚束调制的海浪方位向斜率谱和有效波高的反演算法。通过将经典MPI方法、同极化调制法及本文算法等3种海浪反演方法所得有效波高与浮标数据进行比较,结果显示:本文方法反演得到的海浪有效波高与浮标数据获得的有效波高之间的均方误差为0.79 m,为3种方法中最小。     

基于X波段海洋雷达的风浪联合反演方法研究

提出了一种新的利用X波段海洋雷达联合反演海面风速与海浪谱的方法,该方法不需要额外的信息输入来反演海浪谱。通过利用风速与雷达后向散射强度的经验关系获得海表风速,然后将反演的风速输入风浪谱,通过求解该模拟风浪谱与雷达观测图像谱的约束函数的最小值来确定海浪谱。利用实验数据对反演方法进行了验证,风速、有效波高、主波周期以及主波波向反演的均方根误差分别为1.9 m/s,0.4 m,1.2 s和9.6°,证明了该方法的可行性。     

Comparison of inversion algorithms for HF radar wave measurements

All ocean wave components contribute to the second-order scattering of a high-frequency (HF) radio wave by the sea surface. It is therefore theoretically possible to estimate the ocean wave spectrum from the radar backscatter. To extract the wave information, it is necessary to solve the nonlinear integral equation that describes the relationship between the backscatter spectrum and the ocean wave directional spectrum. Different inversion techniques have been developed for this problem by different researchers, but there is at present no accepted “best” method. This paper gives an assessment of the current status of two methods for deriving sea-state information from HF radar observations of the sea surface. The methods are applied to simulated data and to an experimental data set with sea-truth being provided by a directional wave buoy     

基于海面多普勒谱模型的海浪反演精度分析

Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.     

Measurements of ocean wave spectrum from airborne radar at small incidence angles

This paper proposes the retrieval method of ocean wave spectrum for airborne radar observations at small incidence angles, which is slightly modified from the method developed by Hauser. Firstly, it makes use of integration method to estimate total mean square slope instead of fitting method, which aims to reduce the affects of fluctuations superposed on normalized radar cross-section by integration. Secondly, for eliminating the noise spectrum contained in signal spectrum, the method considers the signal spectrum in certain look direction without any long wave components as the assumed noise spectrum, which would be subtracted from signal spectrum in any look direction for linear wave spectrum retrieval. Estimated ν from the integration method are lower than the one from fitting method and have a standard deviation of 0.004 between them approximately. The assumed noise spectrum energy almost has no big variations along with the wave number and is slightly lower to the high wave number part of signal spectrum in any look direction, which follows that the assumption makes sense. The retrieved directional spectra are compared with the buoy records in terms of peak wavelength, peak direction and the significant wave height. Comparisons show that the retrieved peak wavelength and significant wave height are slightly higher than the buoy records but don’t differs significantly (error less than 10%). For peak direction, the swell waves in first case basically propagate in the wind direction 6 hours ago and the wind-generated waves in second case also propagate in the wind direction, but the 180? ambiguity remains. Results show that the modified method can carry out the retrieval of directional wave spectrum.     

Directional wave spectra measured with a cloverleaf buoy during ARSLOE

The cloverleaf buoy is designed to determine the directional wave spectra with high directional resolution by measuring the vertical acceleration, surface slope, and curvature of the ocean wave surface. This paper describes the properties of the directional wave spectra measured with the cloverleaf buoy during the Atlantic Remote Sensing Land Ocean Experiment (ARSLOE). It is shown that the directional wave spectra measured under relatively constant wind agree fairly well with the similarity spectrum reported previously, but some differences are found in the spectral parameters. The differences in the scale parameters are attributed to unstable atmospheric conditions, though reasonable explanations for those in the shape parameters are difficult now.     

Tsunami Wave Interaction with Data Buoys

To plan for proper mitigation measures, one should have an advanced knowledge of the phenomenon of tsunami propagation from the deep ocean to coastal waters. There are a few methods to predict tsunamis in the ocean waters; one method is the effective use of data buoy measurements. Although data buoys have been used along the Indian waters there has been a tremendous growth in the number of buoy deployment recently. Under the National Data Buoy Programme (NDBP) of India, the 2.2 m diameter discus data buoys were deployed along the east and west coasts of India for measuring meteorological and ocean parameters. It would be advantageous if these buoys could be efficiently used to measure rare events such as tsunamis. Understanding the dynamic behavior of the buoy is of prime importance if a tsunami warning system is to be successful. This may be accomplished through experimental or numerical studies. A comprehensive experimental study has been conducted to understand the dynamic behavior of a wave rider buoy exposed to a variety of waves. It is common that tsunami waves are represented in terms of shallow water waves, namely solitary and cnoidal waves. Hence, in the present study, the discus type data buoy is scale modeled and tested under the action of solitary and cnoidal waves in the laboratory. The time histories of wave elevations, as well as heave and pitch motions of the buoy model, were analyzed through a spectral approach as well as through wavelet transformations. The wavelet approach gives more detailed insight into the spectral characteristics of the buoy motion in the time scale. The harmonic analyses were performed for the cnoidal wave elevations and subsequent motion characteristics that give an insight into the energy variations. The details of the model, instrumentation, testing conditions and the results are presented in this paper.     

数值模式与统计模型相耦合的近岸海浪预报方法

针对数值模式和统计模型预报近岸海浪存在的局限性,构建了数值模式和统计模型相耦合的近岸海浪预报框架,在模式计算格点和近岸预报目标点之间定义一个海浪能量密度谱传递系数,通过经验正交函数分解和卡尔曼滤波方法建立传递系数的统计预报模型并与数值模式进行耦合。经过对近岸波浪观测站1a的预报试验表明:该方法能够提高近岸海浪有效波高预报精度,有效波高的均方根误差降低了约0.16m,平均相对误差降低约9%。进一步试验和分析发现,该方法的预报有效时间小于24h,将海浪能量密度谱经过分解后得到的基本模态反映了近岸波侯的主要特征,海浪能量密度谱传递系数的变化体现了波侯的季节变化特点。     

直驱式波浪能转换系统的研究和性能优化

A direct-drive wave energy conversion system based on a three-phase permanent magnet tubular linear generator （PMTLG） and a heaving buoy is proposed to convert wave energy into electrical energy. Sufficient experimental methods are adopted to compare the computer simulations, the validity of which is verified by the experiment results from a wave tank laboratory. In the experiment, the motion curves of heaving buoy are with small fluctuations, mainly caused by the PMTLG＇s detent force. For the reduction of these small fluctuations and a maximum operational efficiency of the direct-drive wave energy conversion system, the PMTLG＇s detent force minimization technique and the heaving buoy optimization will be discussed. It is discovered that the operational efficiency of the direct-drive wave energy conversion system increases dramatically after optimization. The experiment and optimization results will provide useful reference for the future research on ocean wave energy conversion system.