强潮海湾中的湍流与混合：观测与数值模拟.pdfVIP

1 强潮海湾中的湍流与混合：观测与数值模拟 2 3 连强①③ ，刘志宇②③* 4 ① 厦门大学环境与生态学院，福建 厦门 361102 ； 5 ② 厦门大学海洋与地球学院，福建 厦门 361102 ； 6 ③ 近海海洋环境科学国家重点实验室（厦门大学），福建 厦门 361102 ； 7 * 联系人, E-mail: zyliu@xmu.edu.cn 8 9 国家自然科学基金（编号41476006 ）资助 10 11 摘要：对垂直剪切应力（雷诺应力）、湍动能剪切生成率与湍动能耗散率以及垂直涡粘性系数等 12 湍流参数的现场观测与数值模拟，是理解强潮海湾中垂直混合过程及其动力与生物地球化学效应 13 的基础. 利用坐底 ADCP 流速观测资料、湍流剖面观测资料以及二阶湍流封闭模型，本文对一典 14 型强潮海湾（厦门湾）中的湍流混合过程进行了观测与数值模拟研究. 观测结果表明：在 M2 分 15 潮占优的厦门湾中，湍动能剪切生成率与湍动能耗散率等湍流参数呈现明显的 M4 分潮周期，并 16 具有显著的涨-落潮不对称性，涨潮时强，落潮时弱. 在涨潮期间，整个水体在垂向上近似混合 -6 2 -3 -2 2 -1 17 均匀，湍动能剪切生成率和垂直涡粘性系数的垂向平均值分别达到 510 m s 和 210 m s ; 18 在落潮期间，强湍流则主要局限于厚度为 5～8 m 的底边界层内，且湍流强度随离底距离的增大 19 而逐渐减小. 基于重力位能异常 （potential energy anomaly ）控制方程的分析表明，潮流拉伸 （tidal 20 straining ）是造成湍流强度涨-落潮不对称性的主要原因，潮混合的贡献较小. 基于二阶湍流封闭 21 模型的模拟结果与实测结果的对比分析表明，在弱层化强潮海湾中，二阶湍流封闭模型能够较好 22 地模拟湍动能耗散率与垂直涡粘性系等湍流参数的垂向结构与潮内变化；模拟结果与实测值之间 23 的差异主要源于一些理想化的模型假设和数学近似. 24 关键词：强潮海湾 湍流混合 二阶湍流封闭模型 湍动能耗散率 垂直涡粘性系数 25 26 湍流混合是潮驱海湾中水体垂直交换的主要物理过程. 混合过程引起的物质、动量与能量输 27 运不仅影响着湾内水体的垂向结构，同时也是湾内众多物理与生物地球化学过程的重要控制因子. 28 对湍流通量的准确量化是湍流混合现场观测与数值模拟研究的核心. 近年来，随着高频宽带 29 ADCP 的发展与商业化湍流剖面仪的出现与广泛应用，陆架海及近岸河口海湾中湍流混合的观测 30 研究取得了重要进展，已逐步实现对湍动能剪切生成率、湍动能耗散率等湍流参量的准同步观测 31 [1,2,3]. 这使得从根本上理解强潮海湾中湍流混合特征与机理，并从湍流能量平衡与混合的角度直 32 接检验湍流封闭模型的适用性成为可能. 33 高频宽带ADCP 利用声学多普勒效应测量 beam 坐标下的流速分量. 根据这些原始流速数据， 1 [4] [5] 34 由“方差法”可以计算雷诺应力和垂直涡粘性系数等湍流参数 . Lohrmann 等 首先运用该方法对 35 一海湾中的湍流混合特征进行了观测研究，随后 Lu 和 Lueck[6,7], Stacey 等[8,9] 以及 Rippeth 等[3] 36 进一步发展了该方法，并将其运用到陆架浅海与河口海湾中. 37 湍流剖面仪（如 FLY[10]和 MSS[11]等）利用流速剪切探头测量微尺度流速剪切. 根据这些流 38 速数据，基于湍流的“各向同性假设”可以计算湍动能耗散率等湍流参数. Lueck 等[12]系统地回 39 顾了这类湍流观测技术的发展历史，并详细介绍了相关仪器的工作原理与使用方法.