湄洲湾潮流特性的数值研究.pdfVIP

第 28卷 　第 1期 台 湾 海 峡 Vo l. 28, 　No. 1 　2009年 2 月 JOURNAL OF OCEANO GRA PHY IN TA WI AN STRA IT Feb. , 2009 湄洲湾潮流特性的数值研究 1 2 2 3 刘 　浩 ,潘伟然 ,骆智斌 ,石晓枫 ( 1. 中国科学院烟台海岸带可持续发展研究所 ,山东 烟台 264003; 2. 厦门大学海洋学系 ,福建 厦门 361005; 3. 厦门大学环境科学研究中心 ,福建 厦门 36 1005) 摘要 :本文运用一个含动边界的二维河口海岸动力模型模拟了湄洲湾潮汐潮流的基本特征 ,模拟结 果与实测数据吻合较好. 在此基础上估算了湄洲湾大 、小潮过程的纳潮量 ,并根据湾内保守示踪物 ( ) 的质量 浓度 随涨 、落潮流周期性的变化 ,进一步估算了湄洲湾的水交换周期 ,其半交换和 80%的 交换周期分别为 5d和 15d. 同时可以看出 ,主航道深水区的水交换特性明显强于湾顶浅水区. 关键词 :湄洲湾 ;模型 ;潮汐潮流 ;水交换 中图分类号 : P 7 文献标识码 : A 文章编号 : 1000 8 160 (2 009 ) 0 10 118 05 湄洲湾位于台湾海峡西岸 ,湾内已建有多个船坞码头 ,构成了海上交通运输的重要枢纽. 此外 ,三大石化 基地 ———泉港 、泉惠和东吴石化基地的规划和建设必将使湄洲湾在福建未来的区域经济发展中扮演重要角 色. 另一方面 ,这些港口和石化基地的建设发展又可能给当地的水环境资源造成潜在的危害. 因此了解湄洲 湾的水动力特性对于灾害事故的预防治理具有重要意义. 本文正是基于这一 目的 ,决定采用数值模型的方法 系统地了解湄洲湾的潮汐 、潮流特性 ,并对整个湾的纳潮量和水交换周期给出一个定量的估算. 图 1给出了本文的研究区域. 该区域不仅包括了湄洲湾 ,还包括了相邻的两个小湾 ———平海湾和大港. 图 1 中的 C 1、F和 C2表示 3个潮流测量站位 ,用以模拟结果的验证. 图 1中虚线为假定的湄洲湾的分界线 , 用以估算湾内的纳潮量和水交换特性. 1　数值模型 1. 1　模型介绍 本研究采用的水动力模型是由 B lum berg和 M ellor开发的 三维斜压原始方程海流模型 [ 1 ] 发展而来的 ,适于河口及沿岸海 域的水动力环境的模拟 ,简称为 ECOM. 该模型基本特征如下 : 水平和垂直方向分别采用正交曲线坐标系和随地形变化的 sig ( ) m a坐标系 图 2 ; 自由表面可以模拟水位变化 ;垂直和水平方 向的混合扩散分别采用 2. 5 阶的 M ellorYam ada 湍流闭合模式 和 Sm agorin sk i模式 ; 内外模态分别处理速度较慢的内重力波和 速度较快的外重力波 , 以提高整个模式计算效率等. 1. 2 　动边界的处理方法 ECOM 模型的一个显著特征就是在计算过程中所有网格点 图 1　研究区域 ( ) ( ) 都被看作是湿点 水点 ,即便是干点 陆地点 也被赋予一个很 F ig. 1　M ap s of region s in th is study ( ) 小的水深 小于实际的最小水深 . 这样 ,当每一时间步长的计算 ( ) 完成后再通过一个转换器 FSM 其值在干点处定义为 0,湿点处定义为 1 将干点处的计算结