东京湾温盐长期变化数值计算.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 东京湾温盐长期变化的数值计算· 杨洁余锡平 (水沙科学与水利水电工程国家重点实验室，清华大学水利水电工程系，北京，100084) 摘要：东京湾是一个典型的半封闭型海湾，湾内尤其在湾口处地形起伏较大。研究东京 湾内温度、盐度场空间分布的长期变化及其与水流变化的关系，是进一步研究湾内生态环境 规律的基础．本研究计算采用c0HERENs水动力模型，垂向盯坐标采用混合分层，外海开边界 计入M2，Kl，01，s2等八个主要分潮，河流开边界考虑主要河流的日径流量输入，并采用实 测气象数据估算大气一海表面间的热、盐通量和动量通量等。通过对2005年2—11月份的计算 结果进行分析，并和实测数据相比较，较好地再现了湾内的流场和温盐场的长期变化规律。 关键词：东京湾；CoHBRENS；混合分层；数值计算 1研究背景 按照狭义的定义，东京湾一般是指观音崎一富津岬连线以北的水域，亦称内湾，它通过浦 贺水道与太平洋相连，水域面积约1380km2。东京湾内湾平均水深约15m，而外湾湾口处水深 可达700m左右。湾口处两岬间距约7I锄，为典型的封闭型海湾，内湾与外海间的海水交换情 况较差，海水的贮留时间约为30~45d【11。东京湾背靠人类活动非常频繁的东京都，在内湾有 多摩川、江户川和荒川等主要河流汇入，潮余流、密度流、风生流等是内湾与外海间物质交 换的主要动力因素。因此，潮汐、河川径流等对湾内物质循环都有不可忽视的影响。数值计 算是研究东京湾水域水环境过程的重要手段。Shimizi等12J根据实测的雷达数据，对东京湾的大、 小潮期间潮余流等物理特性进行了数值计算及分析，GuopJ等人也先后对东京湾余流特征进行 了详细的计算和探讨。这些研究也表明，温盐的时空分布与河流淡水流入，大气一海水间的热 通龟交换，降水、蒸发等形成的营养盐通量有着密切的关系。 2计算方法 计算采用水动力生态耦合模型CoHERENS【4】。该模型控制方程基于垂向静力平衡和 Boussinesq近似，可同时求解表征温度、盐度以及污染物等变量的对流扩散过程。模型在垂直 。国家自然科学基金资助 ．918． 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 方向上采用。分层，按内外模态分离的方式求解。 东京湾内湾地形变化幅度较大。湾顶到湾口底部落差达百米。若在整个计算域内采用。 均匀分层，与密度分层的实际构造相差较远，计算结果的精度将受到很大的影响。为了在水 深较大的海域也能较好地获取混合层信息，需要在海面和海底附近加大分层的密度。为此， 本y删l引入垂向。混合分层【51。若按通常的均匀分层方法在垂向分Ⅳz层，则各层界面处的 口坐标为： 七=1～Ⅳ么+1 (1)(1) 吼’：等，|i}：l～ⅣZ+1 吼=面， 其中，七=O时为海底，I|}=ⅣZ+1时为水面。在表层、海层实现局部加密的口非均匀分层可采 用以下形式： 肛型‰黑筹竺一·一¨ 圆 其中，以和刃均为正值，分别为控制表、底层加密的参数，当珥=吃=O时，屏’=％’，非均 匀分层退化为均匀分层。混合分层实质上是a坐标的均匀分层《和非均匀分层群的加权平均： 吼=哆《+【l一哆)群 (3) 其中，权函数 一一m；。l(厦+。一群)一告(《+。一《) q硼ml可再虿商石二丁 权函数的值取决于阎值水深口与总水深D之比值，当总水深小于阀值水深时，混合分层等价 于沿整个水深均匀分层；反之，则由Z，和Z参数控制垂向的分层。 计算中水平网格间距为1km，垂向分16层。对横须贺测站2005年的验潮数据进行调和 摩川、江户川、荒川和鹤见川等四条一级河流的日流量输入过程。温盐的外海开边界结合实 测数据给定或零梯度给定，河流边