珠江河口二维整体数值模拟及其西北江主要分流比研究.pdfVIP

珠江河口二维整体数值模拟及其西北江 主要分流比研究1 许炜铭 包芸 （中山大学力学系 广州510275，） 摘 要：珠江河口是世界是最复杂的河口之一，它由河口湾和错综复杂的网河组成。本文生成了整个珠江 河口包括网河区和河口湾区的整体计算网格。在此计算网格的基础上应用自主开发的计算程序 ccost，对珠 江河口水动力特性进行了二维整体模拟计算，得到了整个珠江河口合理的水动力场。通过二维模拟计算流场 容易显示的优点，可以清晰的看到在整个珠江河口中的流动情况，尤其是网河中各个主要河道内在涨落潮过 程中非同步性。这种在河口区网河河道中涨落的非同步性与潮波在复杂网河中的传播特性间的关联有待进一 步研究。同时对两个主要西、北江汊道的枯季分流比进行了讨论。 关键词：珠江河口 整体数值模拟 潮波传播 1 引言 珠江河口作为典型的复合水域，三江汇流、八口入海决定了河口各区域在地貌形态和动 力特征上均存在着显著的时空差异。在河网区纵横交错的河道主要以一维流动为特征，在河 口区则以一、二维过渡流动为特征，口 门外则是典型的二维和三维流动特征。 河口作为一个整体系统, 其动力及物质 输移等各种特性在空间和时间上是相互 影响的, 兼具陆域和海洋交汇地带的特 殊性质，系统整体模拟是珠江河口数值 模拟的必然趋势[1,2] 。但珠江河口由于其 边界的复杂性，给河口的整体数值模拟 工作带来了一定挑战。国内对珠江口的 数值模型研究近年来从一维、二维、三 维的分维计算模式逐渐向不同维之间连 接计算的模式发展[3,4]，其中河网与口门 外海域的一二维连接计算模型以及一三 维连接计算模型[5] ，使得不同区域水动 力特征的整合分析这一目标逐步实现， 为模拟复杂网状三角洲及河口区动力及 物质输移提供了条件和方法。 图1 研究区域和计算网格 1 国家基金项目资助 针对珠江河口的实际地形边界，本文设计并生成了包含网河区和河口湾的整体计算网格。 在此基础上应用中山大学近岸海洋研究中心自主研究开发的二维河口海洋水沙特性数值计算 模式ccost-2d，对该区域内的流动进行了二维水动力模拟计算，并通过计算得到的整个河口的 枯季水动力场来分析不同区域的潮波传播特征和涨落潮非同步性。 2 研究区域概况和计算网格 本文研究区域（图 1）的上边界至西江高要站、北江石角站、流溪河老鸦岗站、东江博罗 站、潭江长沙站；外海边界延伸至南海大陆架水深约 60m 处。 网河区域与河口湾整体联解的重点和难点在于网河区河道纵横交错，如何分配网格数并 与河口湾相衔接从而生成适用的计算网格。这是一项技术性非常强的工作，也是网河区域与 河口湾整体计算的关键点。本文经过多次的计算和调整，生成了珠江河口覆盖网河区和河口 湾的贴体网格。图 1 所示曲线网格共 352 行、273 列，参与计算网格数 23579 个，外海网格 最宽处 6.5km，河网最窄处网格宽 70m 。该网格为基础网格，未来根据研究的需要，可针对各 重点分析区域进行加密。 3 二维水动力模型 目前常用的二维海洋水动力模式有ADI方法、汉堡模式、POM模式、ECOM模式以及 ECOMSED模式等。这些模式在珠江河口的水动力计算中各有优缺点。在综合考虑原有几种模 型的基础上，中山大学近岸海洋研究中心自主开发了二维河口海洋水沙特性数值计算模式 ccost －2d[6] ，其主要特点包括：采用平面正交曲线网格；采用平面有限体积法求解水动力方 程组；通过旋转矩阵克服科氏力产生的计算弱不稳定问题；水位计算类似于ECOM采用隐式技 术；采用冻结网格法模拟滩地出没水面的动边界技术；对水点记数，只计算求解水点的物理 量；计算高程，即负水深。 从 2005 年发展至今，在珠江口的水沙特性数值模拟中，ccost-2d 模型已经用于伶仃洋、 磨刀门和崖门水域的流场计算，并作为基础水动力、泥沙计算程序，引进到珠江三角洲 6000 年以来地貌环境演变的长周期计算模型和珠江河口 10-100 年环境演变的长周期计算模型中， 均取得了令人满意的效果。 4 模型验证