5 小浪底水库异重流研究与利用展望 - 黄河水利科学研究院.doc

黄河小浪底水库异重流研究与塑造 张俊华 李书霞 陈书奎 马怀宝 王艳平 （黄河水利科学研究院） 摘 要小浪底水库自1954年工程规划至今，历经半个多世纪。伴随着工程规划设计的不同阶段，均不可避免地涉及水库异重流问题。水库施工阶段开展的小浪底水库实体模型及数学模型研究，针对小浪底水库异重流进行更为全面系统的研究。小浪底水库投入运用以来，历次调水调沙不仅将长期的研究成果付诸实施，而且使得对异重流运行规律的研究在认识－实践－再认识的过程中不断深化。进一步加强基本规律研究、物理模型试验及原型资料观测是深化异重流研究的重点。 关键词 小浪底水库 异重流塑造 调水调沙 1 工程规划及设计阶段 黄委设计院在进行小浪底水库淤积过程计算时所涉及到的水库异重流问题，主要是依据三门峡等已建水库的实测资料所建立的经验关系及范家骅提出的异重流排沙计算方法[1,2]。 范家骅异重流排沙计算的基本思路是：根据水库蓄水位及地形特点将库区划分为若干库段，分析确定不同库段及不同高程相应的阻力系数、平均河宽、底坡等；按均匀流流速公式计算第二断面流速，利用建立的流速与悬沙的关系曲线，推算出第二断面的平均含沙量，依此类推，可得到异重流运行至坝前的平均含沙量。 黄委设计院通过对已建三门峡等水库壅水排沙资料的分析认为，水库在蓄水状态，异重流与壅水眀流排沙比可用统一的表达式描述，即： （1） 式中，为浑水容重t/m3；为出库流量m3/s；为水库蓄水容积m3；为悬移质泥沙群体沉速；、根据水沙及边界条件而取不同的值。因此，作为简化方法，利用式（1）进行水库异重流排沙计算。 2 “八五”攻关期间 作为“八五”国家重点科技攻关项目“黄河治理与水资源开发利用”的专题之——黄河泥沙冲淤数学模型的应用[3]，在模型研发的过程中，不可避免地涉及水库异重流计算问题。涉及异重流计算模型包括黄科院曲少军模型、中国水科院孙卫东模型、清华大学王兴奎模型。 黄科院曲少军模型判断异重流潜入条件时，在库底坡降较陡时，采用关系式（式中，、分别为异重流潜入断面平均流速和水深；、分别为异重流与清水的容重）；若底坡缓时，采用韩启为等人的研究成果（为异重流均匀流水深）。 采用异重流均匀流运动水深公式试算得到异重流水深、面积及流速。异重流挟沙力、含沙量及冲淤量等计算与眀流同，其水力因素采用异重流相应值。 中国水科院孙卫东模型对异重流的潜入条件采用表示，其中及分别采用范家骅及异重流均匀运动的水深公式计算。异重流挟沙能力、含沙量、悬移质级配等计算与明流完全相同，其水力因素采用异重流相应值；对于支流异重流倒灌问题，采用李义天在范家骅工作基础上得出的异重流流速及含沙量沿程变化表达式计算。 清华大学王兴奎模型对异重流的运动规律进行了系统的试验研究，得出了异重流时均流速分分布、含沙浓度分布、水流挟沙力等公式，采用黄河三门峡水库的实测资料对所建的公式进行回归计算，得出： 垂线时均流速分布 （2） 其中， （3） 含沙浓度分布 （4） 式中，为垂线最大时均流速，为其相应高度；为距潜入点距离；为潜入断面平均浓度。 在水文水动力学泥沙数学模型中考虑了干支流倒灌淤积。在计算中，按近似水平计算倒灌淤积体修正地形。 3 工程施工期 黄科院在小浪底水库施工期分别采用实体模型与数学模型同时开展了小浪底水库异重流研究，两者达到了互为补充与验证的作用。中国水科院、清华大学、西安理工大学、黄委设计院等单位利用数学模型对小浪底水库异重流进行了研究。 3.1 水库实体模型研究 小浪底水库拦沙初期，库区大多为异重流排沙。为保证水库模型排沙与原型相似，黄科院专门进行了异重流相似条件研究。文献[4]及[5]通过推导非恒定异重流运动方程式及非恒定二维非均匀条件下的扩散方程： （5） （6） 运用相似转化原理得出的异重流发生相似条件，及异重流流连续相似条件分别为： (7) （8） （9） 式中：为浑水容重分布修正系数，若=1则式（7）可转化为异重流发生相似条件的特殊形式；为泥沙容重比尺；角标、分别表示原型及模型。 小浪底库区模型经过两个5年系列试验，重点研究了水库运用初期库区水沙输移及排沙特性