珠江河口整体数值模拟及潮波传播特征研究.pdfVIP

第九届全国水动力学学术会议暨第二二十二届全国水动力学研讨会文集 珠江河口整体数值模拟及潮波传播特征研究1 许炜铭陈祖辉包芸 (中山大学力学系广州510275，E豫il：stsby锄ail。sysu。edu．cn) 摘要：珠江河口是世界是最复杂的河口之一，它由河口湾和错综复杂的网河组成．本研究 生成了整个珠江河口包括网河区和河口湾区的整体计算网格。在此计算网格的基础上应用自 主开发的计算程序ccost，对珠江河口水动力特性进行了二维整体模拟计算，得到了整个珠江 河口合理的水动力场。通过二维模拟计算流场容易显示的优点，可以清晰地看到在整个珠江 河口中的流动情况。对网河中局部区域主要河道在涨落潮过程中的非同步性进行分析，同时 对西，北江主要汉道的枯季分流比进行了讨论。 关键词：珠江河口；整体数值模拟；潮波传播 1引言 珠江河口作为典型的复合水域，三江汇流、八口入海决定了河口各区域在地貌形态和动 力特征上均存在着显著的时空差异。在河网区纵横交错的河道主要以一维流动为特征，在河 口区则以一、-二维过渡流动为特征，口门外则是典型的二维和三维流动特征。河口作为一个 整体系统，其动力及物质输移等各种特性在空间和时间上是相互影响的，兼具陆域和海洋交 汇地带的特殊性质，系统整体模拟是珠江河口数值模拟的必然趋势fl埘。但珠江河口由于其边 界的复杂性，给河口的整体数值模拟工作带来了一定挑战。国内对珠江口的数值模型研究近 年来从一维、_二维、二维的分维计算模式逐渐向不同维之间连接计算的模式发展f3圳，其中河 网与口I’J外海域的一_二维连接计算模型以及一二维连接计算模型【5J，使得不同I又：域水动力特征 的整合分析这一目标逐步实现，为模拟复杂网状三角洲及河口区动力及物质输移提供了条件 和方法。 针对珠江河口的实际地形边界，本研究设计并生成了包含网河区和河口湾的整体计算网 格。在此基础上应用中山大学近岸海洋研究中心自主研究开发的二维河口海洋水沙特性数值 计算模式ccost-2d，对该区域内的流动进行了_二维水动力模拟计算，并通过计算得到的整个河 口的枯季水动力场来分析不同区域的潮波传播特征和涨落潮非同步性。 国家皋会项目资助 ．852． 第九届全国水动力学学术会议暨第二十一二膈全国水动力学研讨会文集 2研究区域概况和计算网格 研究区域(图1)的上边界至硝江 高要站、北江右角站、流溪河老鸦岗站、 东江博罗站、潭扛长沙站：外海边界延 伸至南海人陆架水深约60m处。 网河区域与河口湾整体联解的重点和难 点在丁网闸医河道纵横交错，如何分再己 网格数升与河口湾相衔接从而生成适圳 的计算网格。这是一项技术性非常强的 一『一作，也是网河区域与河口湾整体计算 的关键点。经过多次的计算和调整，生 成了珠江河口覆盖网河医和河口湾的贴 体网格。图l所示曲线网格共352行、 273列，参与计算网格数23579个．外 海网格最宽处65km，河网晟窄处网格 宽70m。该网格为基础网格．未来根据 研究的需要，可针对各重点分析区域进 行加密。 圈i研究区域和计算罔格 3二维水动力模型 目前常用的二维海洋水动力模式有ADI方法、汉堡模式、PoM模式、Ec0M模式以及 EcOMsED模式等。这些模式在珠江河口的水动力计算中并有优缺点。在综合考虑原有几种 模型的基础上，中山大学近岸海洋研究中心自主开发了二维河口海洋水沙特性数值计算模式 程组i通过旋转矩阵克服科氏力产生的计算弱不稳定问题；水位计算类似丁EcoM采用隐式技 术：采用冻结网格法模拟濉地出没水面的动边界技术；对水点记数，只计算求解水点的物理 量；计算高科，即负水深。 从2005年发展至今．在珠江口的水沙特性数值模拟中，ccos-2d模型已经用丁伶仃洋、 磨刀门和崖fj水域的流场计算．井作为基础水动力、泥沙计算程序，引进到珠江三角洲6000 年以米地貌环境演变的长周期计算模型和珠江河口10—100年环境演变的氏周期计算模型中， 均取得了令人满意的效果。 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 4模型验证 模型模拟了枯季期间的水动力特征，模拟时段为2001年2月7一16日共196