《北部湾潮汐潮流的三维数值模拟》.pdf

第 23 卷 　第 2 期 海 　　洋 　　学 　　报 Vol . 23 , No . 2 200 1 - 03 AC TA OCEAN OL O GICA SIN ICA March , 200 1 北部湾潮汐潮流的三维数值模拟 1 1 孙洪亮 ,黄卫民 ( 1 国家海洋局 第一海洋研究所 ,山东 青岛 26606 1) 摘要 : 基于二阶湍流闭合模型计算涡动粘性系数的 POM 三维水动力模式 ,采用细网 格 ,考虑 6 个岛屿、海底摩擦系数进行划片取值 ,模拟北部湾潮汐潮流. 所得潮汐调 和常数与 81 个实测站比较 ,绝对平均误差 : K1 分潮振幅为 4 6 cm ,迟角为 9°;O1 分 潮振幅为 56 cm ,迟角为 7°;M2 分潮振幅为 62 cm ,迟角为 15°. 由模拟结果分析出 该海区潮汐 、潮流 、余水位和潮余流 , 以及水平速度垂直分布等特征. 关键词 : 北部湾 ;潮汐潮流 ;POM ( ) 中图分类号 : P7227 ; P73123 　　文献标识码 : A 　　文章编号 : 0253 - 4 193 2001 02 - 0001 - 08 1 　引言 1) 20 世纪 60 年代初 , 中国和越南 曾经合作对北部湾潮汐和潮流进行观测和分析 . 郑文 振 、方国洪等在 1964 和 1965 年最早对北部湾的潮汐潮流进行系统研究[ 1 ] . 对北部湾潮汐潮 流的数值研究已有 30 多年历史 ,但多数研究将其作为整个南海的一个组成部分. 专门对北部 湾潮汐潮流进行研究始于周朦和方国洪. 曹德明和方国洪用有限差分法模拟北部湾 m1 和 M2 分潮的潮汐潮流[2 ] . 他们采用的均为二维模式. 夏华永等依经 Sigma 坐标变换的三维非线性 潮波 N avier - Stokes 方程 ,用分裂算子法对北部湾 K1 和 M2 分潮进行三维潮流数值模拟[3 ] . 笔者认为基于二阶湍流闭合模型计算涡动粘性系数的 POM 三维水动力模式能够获得更好的 结果. 与以往北部湾潮汐潮流数值模拟相比,本文有以下不同与改进 : ( 1) 把东开边界选在 111°12 ′E ; (2) 选用 3 ′×3 ′细网格. 在模拟范围内,除海南岛外 ,另考虑了涠洲岛等大小 5 个岛 屿 ; (3) 对海底摩擦系数进行划片取值 ; (4) 开边界潮位输入为 O1 + K1 + M2 ; ( 5) 收集了 81 个 潮汐观测站资料. 2 　计算海区及数值模式 计算海区为整个北部湾 ,范围为 16°00 ′～2 1°54 ′N , 105°36 ′～111°12 ′E. 东边界在 111°12 ′ 收稿 日期 : 1999 - 12 - 26 ; 修订日期 : 2000 - 05 - 08 . ( ) 基金项 目: 广西壮族 自治区资助项 目 9734 . ( ) 作者简介 : 孙洪亮 1943 — , 男 , 山东省潍坊市人 , 副研究员 , 现主要从事物理海洋学研究. 1) 中越北部湾海洋