珠江三角洲河网与河口夏季水沙通量的模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 珠江三角洲河网与河口夏季水沙通量的模拟 胡嘉镗 李适宇 中山大学环境科学与工程学院 广东广州 摘要 本文利用一 三维水动力与泥沙耦合模型 计算珠江三角洲河网与河口夏季水沙通量 构建其收支模式 并分 析水沙迁移路径与泥沙沉积特征 研究表明 西江为最主要的水沙输入源 磨刀门为最主要的水沙输出口门 蕉门 次之 夏季泥沙以淤积为主 上游汇入的泥沙 有 沉积于河网区 其余 经八大口门输入珠江口后 有 发生沉积 另外 输入外海 河网区的水沙输送由径流控制 而河口区则由径流 潮汐 季风等因素控 制 河网区各区域的沉积特点因动力条件的差异而呈现不同的规律 大量泥沙在西江干流 虎门水道淤积 珠江 口中以内伶仃洋与磨刀门海域的沉积量最大 泥沙在西滩周边和磨刀门海域快速沉积 其中磨刀门海域淤积最为 强烈 关键词 珠江三角洲 水 悬浮泥沙 通量 数值模拟 中图分类号 文献标识码 在河口系统中 泥沙与水是重金属 营养盐和有机物等多种物质的重要载体 其输移是初级生产 生 物地球化学循环 污染物迁移转化过程的关键环节 珠江三角洲是一个极为复杂的大尺度河口系 统 涵盖一个巨大 独特的感潮河网区和河口湾 三江汇流 八口出海 携带大量的水 泥沙和污染物流 经河网区 经多通道从不同口门汇入珠江口 输往南海 其间经历一系列复杂的生物地球化学反应 据统计 珠江三角洲多年平均径流量和悬沙量分别为 和 其中丰水期 月 水 沙量各占全年的 呈现明显的季节特征 珠江三角洲水系错综复杂 动力因子众多 泥沙行为复杂 为全面掌握水沙的输送通量及收支过程 需借助于数学模型 这不但可以弥补统计上的 不足 而且能够提高通量估算的精度 近年来 珠江三角洲的水沙数值模拟已取得较多成果 但多数研究着重于模型构建以及水动力与泥沙输移特征的分析 并未涉及通量 再者 研究多偏重于河 口区 而面向河网区 尤其是将河网区与河口区视作整体的模拟研究仅为少数 实际上 河网区 河口区 与近岸海域之间水流相互贯通 多种动力交互作用剧烈 难以将其截然分开 此外 河网区的物质输送形 态非常独特 其对物质的滞留 积蓄作用对于通量的确定以及对河口海岸环境的影响至关重要 因此 本文将河网区与河口区作为一个整体河口系统 建立一 三维水动力与泥沙耦合模型 计算河网区与河 口区夏季的水沙通量 构建水沙的收支模式 并分析水沙的迁移路径以及泥沙沉积特征 一 三维泥沙耦合数学模型 一维悬移质泥沙模型 汇入珠江三角洲的推移质不及悬移质的 作用较小 故模型中未纳 入推移质输移过程 一维悬沙输运方程 表述如下 悬沙输运方程 收稿日期 基金项目 广东省近岸海洋综合调查与评价 专项 子课题 珠江口主要环境问题分析与对策 作者简介 胡嘉镗 男 广东深圳人 博士生 主要从事水环境数值模拟与水环境管理研究 汊点连接方程 式中 为过水面积 为水面宽度 为流量 为纵向分散系数 为泥沙浓度 分别为恢 复饱和系数 沉降速度 水流挟沙力 为水位 为汊点水面面积 为汊点编号 为与汊点 相联接 的河道编号 为汊点连接断面个数 模型采用隐式差分迎风格式离散 求解方法见文献 水流挟沙力 泥沙沉降速度公式 如下 式中 为重力加速度 为水力半径 ? 为平均流速 为经验系数 为水体运动黏性系 数 为悬沙中值粒径 分别为泥沙 水的容重 当水体泥沙浓度大于水流挟沙力时 泥沙发生淤积 当水体泥沙浓度小于水流挟沙力 且流速大于 泥沙起动速度 与扬动速度 时 发生冲刷 起动速度公式 扬动速度公式 如下 式中 为水