【讨论】小型水库洪水调节计算

现在一般的调洪计算，大多采用VB或者VBA配合EXCEl来做，因为EXCEL应用比较广泛，大家都比较熟悉。如果应用SQL或者Oracle的话，在大计算量情况下比较有优势。一般的计算似乎有点得不偿失，因此本文仅对一般小型水库洪水调节计算应用VBA和EXCEL进行说明。
首先说明一下，大型水库尤其是洪量大，泄流能力大水库，这样的水库调节起来比较容易。如果有朋友计算过不同类型的水库调洪，会深有感触。大型水库一般只要过了最高水位时段，以后时段水位大多平稳回落到起调水位。小型水库不然，在洪水过程线最开始的若干时段，属于控泄时段，随着来水量加大，水库按泄流能力泄流。达到最高水位后，水位开始回落，此时入库逐渐减小，出库也减小，但出库减小幅度小于入库减小幅度。在整个洪水过程中，一个水位变化的完整过程是：1、起调水位-2、最高水位->3、起调水位。那么小水库的洪水时段划分的应适当小一些，而且3、起调水位出现之前，出库是大于入库的，因此往往存在调节到最后，水位低于起调水位，这种情况是不允许，因为那意味这你可能放走了一部分兴利库容。即使你简单控制低于3、起调水位也不是很容易，这就是与大型水库调节计算之间的差别。大型水库小频率洪水流量相当大，每个洪水时段（h）相对而言就较小，因此调节起来就符合水位变化过程。小型水库，尤其是小型水库的大频率洪水，比如33%，20%，洪水流量一般都不大，但是由于库区堤防等需要这样频率的洪水水位，仍要进行计算。当然，常规的计算能够得到最高水位，但是最高水位出现后的其余洪水过程不对特征水位有影响，但是如果看到一个这样过程的话（即最终水位低于起调水位），这样的计算过程是不完善的，也是不合理的。
因此，我们需要对计算过程进行修正。一般计算过程见下：


H（1）为起调水位，HΔ为某一微量
V（1）为起调水位对应库容，VΔ为某一微量

洪水时段 入 库 出 库 水库水位 库 容
1 Qr（1） Qc（1） H（1） V（1）
2 Qr（2） Qc（2） H（2） V（2）
................................
最高水位出现在两虚线之间
................................
n Qr（n） Qc（n） H（n） V（n）
n+1 Qr（n+1） Qc（n+1） H（n+1） V（n+1）

假设第n个时段 H（n）>H（1）,第n+1个时段H（n+1） 那么n+1个时段的出库流量Qc（n+1）=（V（n）-V（1）+Qr（n+1））/时段长，
这样，出库流量才是正确的，第n+1时段末水位及库容均与起调时相应，其以后时段按入库泄流。整个过程完全正常！

附件为某一小水库的洪水调节计算过程，在小频率的时候特别明显，即时段末无法正常回到起调水位。因为以前做大水库的时候，根本遇不到这样的问题，因此，把这个提出来希望和大家一起探讨。
附件的excel程序中带有宏，没有病毒请放心，另外，那个excel都是比较简单的VBA，我不是跟大家探讨编程序的问题，相信论坛很多都是编程高手，因此那个excel也不是通用性很好，仅作交流用，希望不要把它用在实际工程计算上，不然出问题本人概不负责！
呵呵！


[ 本帖最后由 nidr_zys 于 2011-2-13 16:17 编辑 ]