黄河泥沙数学模型的研究进展与展望.pdfVIP

354 新世纪水利工程科技前沿 黄河泥沙数学模型的研究进展与展望 张红武钟德钰 (清华大学 黄河研究中心北京 100084) 摘要黄河泥沙数学模型是同泥沙学科发展关系最为密切且难度最大的课题之一。本文论述 了半个世纪以来黄河数学模型的研究进展和展望，认为河流泥沙运动理论实质性的深化及对黄 河河性的深刻把握是黄河泥沙数学模型得到发展的动力，也只有如此，才能使数学模型的研究水 平不断提高。 关键词黄河泥沙数学模型进展展望 河流泥沙数学模型是建立在流体力学、计算数学、泥沙运动力学等学科基础上，用于模拟河流水沙运 动、河床演变过程的重要工具之一，在河流工程中广泛应用。近年来，一些优秀的数值方法、格式在应用过 程中被保留和继承下来，同时又有一些新方法、新格式不断被引入，加上计算机技术的飞速发展以及泥沙 学科取得的进展，促进了水沙数学模型的发展，目前主要工作逐渐转向模型功能的扩展、前后处理技术、计 算可视化等的开发研制。还需指出的是，数学模型的发展水平不仅取决于求解相关方程的计算方法和计 算机的计算能力，更莺要的是取决于泥沙学科的研究进展。 黄河泥沙计算方法研究起步较早，早在20世纪50年代已有许多学者通过计算来解决黄河的实际问 题，由于黄河泥沙问题极其复杂，泥沙学科的发展速度不快，相对于长江等河流的数学模型发展较为缓慢。 黄河数学模型真正得到快速发展是近10多年来的事情，逐渐在治黄研究中发挥了莺大的作用。本文在简 要回顾黄河数学模型初期与发展阶段的研究状况的基础上，着重介绍了黄河数学模型中的困难和最新进 展，并对其未来发展做了展望。 1 黄河数学模型初期与发展阶段 计算方法是数学模型的基础和核心，因而初期的研究主要以黄河泥沙计算方法为主要目标。1954 年，黄河规划委员会在进行流域规划时，开展了三门峡水库的淤积计算和三门峡水库建成后下游河道冲刷 的计算。计算中虽区别于以前假定泥沙全部淤积在死库容内的方法，但存在许多粗糙的假定和其他问题。 其后前苏联专家罗辛斯基根据中苏技术合作协定来华指导工作，又引进了新的方法。其计算的总冲淤量 与实际情况比较接近，但滩地和主槽的冲淤量与实际相差甚远。 1962年，三门峡水库因库区严重淤积而改变运用方式后，需要定性判断河床变形发展趋势，而当时的 计算方法不能反映挟沙能力在河床冲淤过程中不断调整的规律。麦乔威、赵业安、潘贤娣等通过分析各水 文测站实测资料发现，黄河下游河床的冲淤变化取决于上游来沙量与水流挟沙能力问的对应关系，并建立 了著名的多来多排计算公式。黄委水科所与水科院河渠所用此方法合作研究了三门峡水库低水位运用后 不同防洪排沙设施下黄河下游河床冲淤变化。 、 1964～1965年，黄委水科所与水科院河渠所一起进行了多沙河流拦洪水库下游河床变形计算。 1969—1972年，张启舜针对三门峡水库的改建和确定改建后的运用方式，进行水库泥沙冲淤过程计算，提 出“三门峡工程改建规模分析”、“三门峡水库(潼关以下)冲淤计算方法”。1972年，清华大学水利系用三 门峡水库历年冲淤资料，参考有关水库的情况，对三门峡库区冲淤的基本规律进行了研究，探讨了库区淤 作者简介：张红武(195卜)．男．教授．博士生导师，黄河水利委员会副总工程师，中国水利学会泥沙专业委员会副主任。 第4部分工程水文、泥沙与水力学 355 积及冲刷的计算方法。1972年黄委规划设计大队根据三门峡水库滞洪及蓄水运用时的排沙资料，建立了 水库排沙比与较多复合因素相关的定量关系式。 1975～19r77年，在进行治黄规划过程中，为研究各种规划方案对黄河下游的减淤作用，黄委及其有关 科研、设计部门的科技人员，使用两种河道冲淤计算的经验模型开展了大量的计算工作，麦乔威、李保如两 位学者进行了总结。其中，在规划方案的初步比较阶段，由于方案较多，工作量很大，使用的是李保如提出 的框算模型，只考虑影响输沙能力的最主要因素(来水、来沙)，计算时段分汛期及非汛期两个时段，计算中 龙门以下干流河段划分为四个河段。对于问题比较复杂和初选的方案，则使用较为详细的经验模型。该 模型与实测差别不大，且两种模型计算的结果也相差不多。 随着计算机技术的快速发展和广泛应用，