水库调洪计算理论方法及研究应用进展.docVIP

水库调洪计算理论方法及研究应用进展 　　 【摘要】在径流调节计算中，起着重要作用的是水库调洪计算。关于水库调洪计算的理论方法，按对水库入流过程的描述方式和调节计算所采用的数学方法，大体可以分为时历法、数理统计法和随机模拟法。文章由此对水库调洪计算的理论方法及研究、应用进展进行了综述。 　　 【关键词】水库调洪计算；理论方法；研究应用进展 　　 　　降水到达地面产生径流后，径流在流域的流动过程中必然会受到坡地、河网、湖泊及人工水利工程等的调节作用，使其在时间、空间上的分配发生了变化，因此，在研究径流分配以及相关问题时，必须考虑到径流调节所起到的作用。其中起到径流调节作用较为突出的就是水库调洪。 　　一、水库调洪计算的理论方法 　　在引言中已经提到，水库调洪计算的理论方法可以分为三类，以下做出介绍： 　　1. 时历法[1]。即“先调节计算后频率统计”的方法。根据过去按时历顺序的流量资料进行调节计算，再将调节后的调节流量、水库水位、库容的多年变化情况，绘制成相应的频率曲线，将供设计选择。 　　2. 数理统计法[1]。即“先频率统计后调节计算”的方法。其过程是先将流量系列进行频率统计，用统计特征值来描述系列的多年变化，再按概率论原理和水量平衡进行调节计算，获得调节流量、多余、不足水量及水库蓄水等多年变化情况，以供决策。 　　3. 随机模拟法[2]。先根据实际径流，应用水文时间序列分析的理论和方法（如水文随机模拟，水文系列资料的生成），建立径流序列的随机模型并据此模拟出足够长的径流系列，然后按照时历法进行调节计算。 　　二、水库调洪计算研究应用进展情况 　　 （一）水库调洪计算的传统方法 　　从上面介绍的三种水库调洪计算的理论方法中，可知时历法和数理统计法是比较传统的方法，在这两种方法的理论基础上，结合现在不断发展的计算机技术、通信技术等发展起随机模拟法，这是当今研究水库调洪计算理论方法的新趋势。所以，首先对水库调洪计算的传统方法作一下介绍。 　　水库调洪计算的理论方法是基于水力学的圣维南方程组进行了简化，忽略了洪水入库至泄洪建筑物间的行进时间、沿程流速变化及动库容等的影响[3]。经过简化后得到专用于水库调洪计算的实用公式，即水量平衡方程为： 　　（1） 　　式中 Q1 ，Q2 ─分别为计算时段初、末的入库流量 （m3/s）； 　　q1 ，q2 ─分别为计算时段初、末的下泄流量 （m3/s）； 　　V1 ，V2 ─分别为计算时段初、末水库的蓄水量（m3）； 　　ΔV ─ V1 和V2 之差； 　　Δt ─计算时段。 　　再结合蓄泄流量关系式：（2） 　　如果要考虑动库容的影响，需在公式（2）的右边加上一个参数Qe（入库流量）。不管水库调洪计算采用何种方法，都是用以上两个方程式或由二者变形而来。 　　对每一计算时段来说，公式中的V2 ， q2 为未知数根据水力学知识[3]，可推出泄流量q 与库容V 的关系曲线公式为（2）。将该公式与水量平衡方程组成方程组，即可进行求解。根据上述方程组进行调洪计算的具体方法常用的是列表试算法[1]、图解分析法[1]（即半图解法）和图解法[1]。 　　1.列表试算法。列表试算法是根据已知条件对每一时段的q2和V2 进行试算，一次次地假定数值、再查表或查图。计算结果不合理时要重新计算，直到符合要求为止。通常要试算多次才能得到接近值，精度还是能够保证，其最大的缺点是运算量较大，而且洪水经水库调节后的最大下泄流量qm。 　　2.图解分析法。图解分析法是由公式 （3）变形得 ： 　　 （3）其中 =（Q1+Q2）/2； 　　作两个辅助函数： （4）； （5） 　　并据此作图，然后进行量线分析，得到洪水下泄过程线。 　　图解分析法相对于试算法来说，计算过程简便，但必须借助于辅助曲线，对每一时段的q2 和V2 通过量取、查找得到计算结果。此法对制图的精度尤其重要，受图形的限制，使欲量的线段在开始时会特别短，这样就会造成一开始就会量取不方便，且有较大的误差，因此，所得结果的精度受人为因素的影响较大。 　　3. 图解法。图解法是由公式（3）变形得 　　（6） 　　 　　作辅助线： （7） 　　再结合入库洪水过程线，可以得到水库下泄过程线。 　　图解法可以克服列表试算法和图解分析法的一些缺点，但它对工作人员的制图能力和熟练程度要求很高。 　　4.综合法。因为以上三种传统的水库调节方法都存在着不同程度的不足，所以针对它们各自的优缺点，进行取长补短。如图解分析法受图形的限制，使得欲量的线段在开始时会特别短，这样就会造成一开始就有了误差，为了解决这个问题，可以在前一两个时段Δt内使用列表分析法，也能较快地算出前一两个时段Δt内的q2、V2，然后再使用图解分析法再