台湾海峡台风浪的数值模拟.docVIP

台湾海峡台风浪的数值模拟 袁凯瑞1,商少平1,2,3*,谢燕双2,3,张莉1,张余得1,张峰1 1.厦门大学海洋与地球学院, 2.厦门大学海洋观测技术研发中心,3.水声通信与海洋信息技术教育部重点实验室厦门大学,福建 厦门 361102 摘要:本文采用第3代海浪模式SWAN,在考虑了波浪折射、底摩擦、破碎、白浪、风能输入、非线性波波相互作用等物理过程的基础上,分别应用Jelesnianski模型风场、藤田气压公式计算的梯度风场以及考虑台湾海峡及台湾岛地形影响的台风风场陈德文台风风场模型,数值模拟了0813号“森拉克”、0815号“蔷薇”和0908号“莫拉克”3场台风过程在台湾海峡产生的台风浪,并与浮标观测结果对比表明,陈德文台风风场模型数值模拟的台风浪结果与实测符合较好,其中模拟结果有效波高平均绝对误差为0.65 m,有效波高峰值平均绝对误差为0.16 m,说明本文对台风浪的模拟尤其是灾害性台风浪的模拟是成功的。本研究可为台湾海峡台风浪的数值预报提供参考。 关键词:台湾海峡;台风浪;SWAN;数值模拟 中图分类号:P 731.2文献标志码: A我国海岸线长,台风活动频繁,由台风引起的台风浪是我国主要海洋灾害之一,台风浪造成的经济损失平均每年超过7亿元。因此对台风浪的研究和预报显得尤为重要。Ou等[1]首次用SWAN模型对台湾东岸与西岸浅水区域的台风浪进行过模拟,但结果并不理想。陈希等[2,3]运用SWAN模型初步分析了台湾岛邻近海域的台风浪分布特征。张军等[5]用SWAN模型并考虑波流相互作用模拟了兴化湾的台风浪,结果与实测符合良好。姬厚德等[6-8]用SWAN模型通过大小网格相互嵌套的方法对台湾海峡及厦门湾台风浪进行模拟,验证了该模型在近岸浅海良好的适应性。姬厚德等[9]通过SWAN模型与缓坡方程联合应用,较好的模拟了“莲花”台风在泉州湾的台风浪。上述研究中,模拟台湾海峡近岸台风浪时通过大小网格嵌套、考虑波流相互作用以及与缓坡方程联合运用等方法有效的提高了SWAN在近岸或者湾内台风浪模拟的准确性,而对于SWAN在台湾海峡中部海域模拟台风浪研究还较少。本文分别采用Jelesnianski台风风场模型、藤田气压公式梯度风风场模型以及带台湾海峡附加风场的台风风场模型 [10]模拟台湾海峡中部海域3场台风的台风浪,并且与实测数据进行对比分析。 1 数值模拟 1.1 海浪模式 本文采用荷兰Delft大学开发的第3代海浪模式SWANSimulation Waves Nearshore的最新版本4085版,该模式全面考虑了波浪浅化、折射、底摩擦、破碎、白浪、风能输入、非线性波波相互作用,并且采用了全隐式有限差分格式,使得计算空间网格和时间步长不会受到限制。SWAN模式采用二维动谱密度Eσ,θ表示随机波,主要是为计算浅水与近岸海浪而设计的,但Rogers等[11]运用了新的数值计算方法,减少了在计算尺度比较大的海域时会产生的数值扩散,从而使得SWAN适用于浅水到大洋各种尺度。 1.2 台风风场本文采用3种不同的台风风场模型进行数值模拟。 1.2.1 台风风场模型IJelesnianski1966台风风场模型由圆形台风风场和移行台风风场叠加组成,即: 。其中圆形台风风场为: 1 其中为计算点风速,圆形风场的风速,移行风场风速,R为最大风速半径,为台风最大风速,为入流角这里取30°,r为计算点到台风中心的距离,,分别为x,y方向的单位向量;移行风场为:2其中为计算点的移动风速,为台风中心的移动速度。 1.2.2 台风风场模型Ⅱ采用藤田气压公式计算的梯度风场叠加Jelesnianski移行风场。其中藤田气压公式为: 3 4 其中为计算点的气压,为台风外围环境气压,台风中心气压,为科氏参量,为空气密度,其他参数与上相同。 1.2.3台风风场模型Ⅲ采用陈德文[6]考虑了台湾海峡及台湾岛地形影响下的台风风场模型,该模型是在圆形对称风场模型Holland台风风场模型及Jelesnianski移行风场的基础上,考虑台湾岛和台湾海峡地形影响设计的台风风场模型。 圆对称风场由Holland提出的飓风径向风速剖面方程简化得到的,它忽略了柯氏力的影响,其径向风速剖面方程如下式: 5其中为空气密度,B为热带气旋参数,为台风中心气压降,其他参数同上; ,6 7 其中,为风速,为圆形风场模型风速,为地形作用产生的附加风速,为近台风中心最大风速,为登陆台风风速衰减因子,为台风中心位置权重因子,为海峡内各点位置权重因子,为风速增强极大因子,具体表达见相应文献。 1.3 数值模拟计算域如图1110°~130°E,18°~30°N,网格分辨率为6′×6′,时间步长为900 s,谱分布在360°范围内,方向步长3°,离散频率为0.125 Hz。选用的物理过程包括了波浪浅化、折射