长江中下游江湖水沙调控数值模拟 Numerical simulation of flow-sediment transport on the river-lake in the middle and lower reaches of the Yangtze River.pdfVIP

第2l卷第6期 水科学进展 V01．21．No．6 2010年11月 ADVANCESINWATERSCIENCE Nov．．20lO 长江中下游江湖水沙调控数值模拟 施 勇，栾震宇，陈炼钢，金 秋 (南京水利科学研究院，江苏南京210029) 摘要：以长江中下游防洪系统为对象，在大型复杂防洪系统洪水运动数值模拟基础上，成功地将面向长江中下游防 洪规划论证需求的洪水演进数学模型转化为面向长江防洪系统江湖水沙调控需求的长江中下游江湖水沙调控数学 模型。为适应江湖水沙调控和评估的要求，提出了基于水动力学的闸坝调度计算模式。此外，还对河网分汉泥沙分 配模式进行了深入研究。通过长江中下游防洪规划及其洞庭湖区“控支强干”方案论证模拟计算，较好地解决了 防洪措施蓄泄后效评估和工程优化调控模式等关键难题，为防洪规划方案的制定提供了定量依据，主要成果已应 用于长江中下游防洪规划方案。 关键词：水沙调控；数值模拟；防洪规划；长江中下游 中图分类号：P333．9 文献标志码：A 文章编号：1001-679l(2010)06-0823-09 利用水流数值模拟试验评价工程对水流的短期影响，该方法广泛应用于江河湖泊的防洪规划论证中¨…。 针对工程对河道水沙的长期影响及其河势的变化，仅依靠水流模型进行试验是不够的，不能反映河床演变带 来的影响一·，需要采用水沙数学模型开展工程运行后的水沙演变模拟试验。闸坝调控改变了原河道水沙演 进规律，改变了河网区中各水道的水量沙量的分配，利用水库或闸坝控制水沙资源，使下游河道朝着有利的 河床演变方向发展”·。水沙调控数值模拟是解决闸控条件下河网分汉河道的水沙优化配置的重要手段之一， 但其中两方面的技术难点需要解决：一是水库、闸坝调控方式的数值模拟；二是闸坝水沙计算模式与坝下河 道的水沙优化配置优化模式的耦合求解。近年来，笔者针对这一技术问题，开展了相关研究。在长江中下游 洪水演进数学模型的基础上p。，围绕水沙输运、河床演变和闸坝水沙调控等相关的关键问题开展深入研究， 同时针对三峡建成后江湖关系变化对区域防洪形势的影响预测、防洪对策措施及调度预案的后效评估，开展 面向应用的研究，取得了预期成果∽。。本文以长江中下游防洪系统为对象，概述了近年来在防洪系统水沙 调控数值模拟关键技术方面的主要成果。 1 江湖水沙数值模拟算法 河网水沙计算采用在同一时步内先计算水流后计算泥沙的非耦合联解，其中水流控制方程是圣维南方程 沙及河床变形计算采用隐式迎风格式旧1，即将泥沙连续方程和河床变形方程在相邻断面间离散成微段方程， 结合汉点处沙量平衡方程，进行逐微段计算，得到河段各断面上的含沙量解。然后利用河床变形方程，根据 河段内泥沙冲淤量与该河段河底高程变化量的关系，计算河床冲淤变化量。 收稿日期：2009-10．16 2007Shz02(1-9)) 作者简介：施 勇(1964一)，男，江苏南通人，教授级高级工程师，主要从事水力学及河流动力学研究。 E—mail：mrshiyong@263．net 万方数据 824 水科学进展 第2I卷 1．1 一维河网泥沙输运算法 河网泥沙算法包括河道泥沙算法、河网汊点分沙模式、河网区内边界点含沙量计算以及河床变形计算。 其中，一维河道泥沙连续方程的离散采用隐式逆风差分格式，即当河道流速大于0时(设沿河道下泄方向流 速为正，反之为负)采用隐式后差，当河道流速小于0时采用隐式前差。该格式形式简单，计算稳定性好， 符合特征理论，能较好地模拟涨落潮不同流态下含沙量的时空变化和水沙输运的往复特征一。。 对于汉点泥沙分配算法，在原计算模式一3的基础上进行了改进。首先采用泥沙垂线分布公式1计算主 支汊河段首断面的含沙量分布，再根据汊点处沙量平衡方程(进出每一汉点的沙量与该汉点水体中的沙量增 减率相平衡)对主支汊首断面的计算含沙量进行校正。具体计算步骤如下： 计算断面的平均含沙量。 t