西北太平洋海浪的特性.docx

《西北太平洋海浪的特性.docx》由会员分享，可在线阅读，更多相关《西北太平洋海浪的特性.docx（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、西北太平洋海浪的特性应用中国海洋大学海洋遥感研究所基于Hasselmanns提出的MPI方法用C语言开发的SAR反演海浪方向谱的软件，对2003年1月2012年1月间西北太平洋海域Envisat波形式数据进行海浪方向谱反演，共得到西北太平洋海域观测的海浪方向谱146796个。统计由反演的海浪方向谱得到的海浪有效波高、平均波向、平均波周期等数据，分析了西北太平洋海浪场分布的特点，得到了一些有参考意义的结果。关键词：SAR海浪方向谱；海浪特性；西北太平洋早期的海浪观测主要依靠于现场定点浮标获得的时间序列海浪有效波高数据，在时间和空间上都有一定的局限性。空间卫星技术和传感器技术的发展为海浪观测提供了

2、新的技术手段。星载合成孔径雷达SyntheticApertureRadar，SAR是一种主动式微波成像雷达，它对海浪的独特的观测能力始于1978年Seasat-A/SAR数据1。20世纪90年代后，以欧空局发射的ERS-1、2/SAR和Envisat/ASAR为代表，世界各国发射的星载SAR为海浪研究和预报提供了强有力的支持。1991年，德国著名科学家Hasselmanns首先提出了卫星SAR海浪图像谱与海浪方向谱之间的非线性理论关系，并基于数据同化的概念提出了从SAR图像反演海浪方向谱的方法(MaxPlanckInstitute，MPI)1。随后，众多学者亦展开了SAR反演海浪方向谱的研究2

3、-8。太平洋处于亚洲大陆、大洋洲、南极洲、南美洲、北美洲之间，是四大洋中面积最大、深度最深的。我国接近西北太平洋，20世纪60年代以来，中国海及西北太平洋等海域，都遭受太多次重大和特大灾祸性海浪的冲击，致使沿海海上和沿岸近海渔业的生产、港口码头、交通运输、盐业生产等造成多起宏大经济损失和人员伤亡事件。因而，深化研究其海浪场的时间变化规律及空间上的分布特性，对我国有极其重大的意义与应用前景。诸多学者利用各种数据源现场、高度计、形式研究过西北太平洋的海浪场9-15，结果表明：(1)西北太平洋季平均和月平均有效波高，冬季最大，高达2.0m以上，有效波高2.54m的5级大浪频率小于20%；春季、秋季，

4、高达1.5m左右，大浪频率小于13%；夏季最小，小于1.0m，大浪频率小于10%；(2)整个太平洋的有效波高呈现明显的上升趋势，大约每年增加0.53cm。目前为止，未见利用SAR海浪方向谱数据进行海浪场时空分布特征研究的文献。本文应用中国海洋大学海洋遥感研究所基于Hasselmanns提出的MPI方法用C语言开发的SAR反演海浪方向谱的软件16-17反演得到的海浪方向谱数据，对西北太平洋海域海浪场的分布特征进行了研究。2方法与数据2.1MPI反演方法星载SAR并非对海浪直接成像，SAR仅与海面短重力波或毛细波互相作用进而成像。由于海浪对海面短重力波的调制作用，因此在SAR图像上可观测到作为调制

5、信号的海浪信息，表现为明暗相间的波纹图像。从SAR图像能够反演获得海浪方向谱。1991年，Hasselmann等1根据SAR工作原理和流体力学理论提出包含倾斜调制、流体力学调制和速度聚束的SAR海浪成像机制，将二维海浪谱同SAR图像谱联络起来，推导出从海浪方向谱到SAR图像谱的非线性变换关系。对于垂直极化和右视的SAR，各调制函数表示如下：在此只给出HasselmannsSAR对海浪成像的完全非线性关系的式式6及其一阶近似准线性式式9(具体推导经过见Hasselmanns等的文章1：式中：n表示非线性阶数，m为速度聚束参数的阶数，kx表示方位向波数。式中：PRql为纯RAR谱，Pvbql为纯准

6、线性速度聚束谱，Pintql为准线性干预项。Krogstad根据Hasselmanns推导出的海浪方向谱到SAR图像谱的非线性变换关系，通过迭代求逆方法得到最优海浪方向谱的MPI反演算法2。MPI反演流程图见图1。图1能够看出MPI算法的框架核心是一个迭代求逆经过。该经过需要SAR图像谱和第一猜想谱作为输入项。通过迭代来不断修改海浪方向谱，使得第一猜想谱通过非线性变换得到的模拟SAR图像谱尽量接近观测的SAR图像谱。定义了价值函数J式12，当J取最小值时所得的海浪方向谱F(k)为最适海浪方向谱。迭代初始所用的海浪方向谱称为第一猜想谱，它由海浪数值预报形式WAM计算得来。WAM的运行需要海面风场

7、和海底地形作为输入。式中：P(k)为观测SAR图像谱，由观测的卫星SAR图像作FFT变换得到。P(k)为最优SAR图像谱，由SAR图像谱和海浪方向谱的非线性变换关系的前向模拟得到。F(k)是由海浪数值预报形式得到的第一猜想谱。F(k)代表最适海浪方向谱。由J的定义能够看到，价值函数由两部分组成，一部分是代表前向模拟SAR图像谱与卫星观测SAR图像谱之间的差，另一部分是海浪方向谱与猜想海浪方向谱之间的差。价值函数最小化意味着既要知足反演得到的海浪方向谱和猜想海浪方向谱最接近，又要知足前向模拟SAR图像谱与卫星观测SAR图像谱最接近。为代表从最优SAR谱到一阶估计谱的置信区间的权重系数，通常取0.

8、1P2max；B为一个小正数，以避免F(k)=0时价值函数的计算值出现奇异值无穷大，通常取0.01Fmax。因而问题变成求一般非线性变分问题的解：此方程的解采用准线性变换关系式(9)进行迭代求出。谱的初值为第一猜想谱。假设Fn(k)为经过n次迭代后的近似解，Pn(k)可用式(6)进行计算。若改良解：由Fn求得Fn+1，然后进行下一次迭代，直至知足一定的条件为止。假设：在大多数的情况下，选取：中国海洋大学海洋遥感研究所基于上述MPI方法应用C语言开发了SAR反演海浪方向谱的软件。软件构成框架图如图2所示。应用该反演软件分别在全球以及中国海和西北太平洋海域内，对反演的海浪有效波高数据与现场浮标有效

9、波高数据进行了同步印证。结果显示：在全球海域，反演误差为5.2%。在中国海和西北太平洋反演误差为4.9%。比照国际上的反演精度，ORSI基于MPI方法的海浪方向谱反演软件在中国海和西北太平洋的表现性能与其相当9-10。2.2数据介绍2.2.1SAR波形式数据Hasselmanns提出的MPI迭代循环需要观测的SAR图像谱作为输入项。本文用于反演海浪方向谱的数据是Envisat_ASAR波形式单视复数据Level1B产品ASA_WVI_1P，来源于欧洲太空局EuropeanSpaceAgency，ESA。其数据是一个个的小图像，大小为10km5km，小图像在轨道方向的间距为100km，极化方式为

10、VV或者HH，天天在全球海域可得到大约2500个小图像。2.2.2数值预报形式数据MPI反演算法需要海浪数值预报形式提供初猜谱信息。本文第一猜想谱使用的是WAMcycle4.5的数值形式结果。WAM的运行需要风场和地形数据驱动。本文采用ECMWFERA40再分析形式风场数据，其空间分辨率为11，时间间隔为6h。，采用全球地形数据：，其空间分辨率为5。3反演结果图3是SAR反演软件的输出结果示例图。输出的产品是原始SAR图像和它对应的反演结果。图3a为SAR图像；图3b为SAR图像谱，由SAR图像通过FFT变换得到；图3c为第一猜想谱，是WAM形式预报结果；图3d为最优谱，即反演得到的海浪方向谱

11、,并且能够输出海浪的有效波高Sig.Height、平均波长Avg.Length、平均周期Avg.Period、平均波向Avg.Dir等海浪参数。本文利用该SAR反演海浪方向谱软件对2003年1月2012年1月间西北太平洋海域的Envisat波形式数据进行海浪方向谱反演，共获得海浪方向谱观测数据146796个。统计由反演的海浪方向谱得到的海浪有效波高、平均波向、平均波周期等数据，研究了西北太平洋海浪场的各季节性分布特征。4分析与讨论利用反演得到的146796组海浪有效波高及平均波向数据，将有效波高按00.5m，0.51m，11.5m，1.52m，22.5m，2.53m，33.5m，3.54m，4

12、m，波向按N,NNE,NE,ENE,E,ESE,SE，SSE，S，SSW，SW，WSW，W，WNW，NW，NNW进行频率统计，以2030N为过渡区域，图4中B区域研究中发现该海区南北海域海浪场分布特征差异较大，故由北至南分3个区域进行分析，详见图4，得到了西北太平洋海域春夏秋冬四季的海浪统计玫瑰图见图58。综合图58以及统计的有效波高、平均波周期数据分析，得西北太平洋的海浪场各季节分布特征如下注：本文浪向均代表海浪传播的方向。春季，西北太平洋A、B、C海区海浪统计玫瑰图如图5所示。西北太平洋30N以北海域，浪向分布较为杂乱，但以E-SSE浪向为主，频率达20%以上。2030N海域内，浪向由以北

13、海域盛行的东偏南浪向向20N以南海域盛行的西偏南浪向过度，总体以偏南浪向为主。春季西北太平洋的平均有效波高为1.8m左右，平均周期在7s左右。夏季，西北太平洋A、B、C海区海浪统计玫瑰图如图6所示。西北太平洋30N以北海域西南季风盛行，浪向总体以东偏北向E-N为主，整体发生频率可达40%以上。2030N海域则以西偏北浪向W-N为主，整体频率达55%。20N以南海域则以西偏南浪向W-S为主，发生频率48%左右。夏季西北太平洋的平均有效波高在1.2m左右，平均周期在6s左右。秋季，西北太平洋A、B、C海区海浪统计玫瑰图如图7所示。西北太平洋30N以北海域浪向分布较为平均，但以偏东浪向NE-SE为多

14、，频率在33%左右，其中以E为主向浪。2030N海域则是西偏南浪向W-S为主，整体发生频率可达65%。20N以南海域盛行偏西浪WNW-SW，发生频率72%。秋季西北太平洋平均有效波高为1.6m左右，平均周期在67s左右。冬季，西北太平洋A、B、C海区海浪统计玫瑰图如图8所示。这是西北太平洋海域海浪场最强的季节。30N以北海域盛行西风带作用突出，浪向为东偏南向E-S，频率在46%以上。2030N的过渡海域，则整体以偏南浪向SE-SW为主，整体发生频率可达55%。20N以南海域则以西南浪向WSW-SW为主，整体频率60%以上。冬季西北太平洋平均有效波高在2.4m左右，平均波周期约78s。鉴于近10

15、a的数据积累，本文统计了西北太平洋海浪有效波高在20032011年间平均年有效波高的变化趋势见图9。发现西北太平洋的年平均有效波高在这几年间整体呈现出微弱的上升趋势。通过拟合得到其上升趋势约为1.8cm/a，这与之前的学者研究的北太平洋的海浪有效波高上升趋势0.53cm/a9-10的结论不冲突。其中2005年的异常情况可能是由于整年间有8个强台风和热带风爆发生，其中有6次强台风的中心风速大于等于45m/s所致；而2010年是我国气候异常、极端、罕见创世纪之最的一年，海浪有效波高也发生异常。5结论与瞻望本文利用ORSI基于MPI方法应用C语言开发的SAR反演海浪方向谱的软件，反演得到了西北太平洋20032011年间146796组海浪方向谱数据。利用这些数据分析了西北太平洋的海浪场分布特点，发现以2030N为过渡海域，其南北海域海浪场分布特征差异较大。同时发现西北太平洋海浪有效波高在20032011年间整体呈现出微弱的上升趋势，通过拟合得到其上升趋势约为1.8cm/a。但本文应用的是基于MPI方法的海浪方向谱反演方法，其反演精度相比于PARSA反演算法略低，况且应用观测的SAR图像反演海浪方向谱以研究海浪场的分布特征，必须收集大量的SAR数据。以后的工作中，我们会致力于提高反演算法的精度，继续收集长时间序列的SAR数据，愈加深化的分析西北太平洋的海浪场特性。