一种基于水库调度的河流生态流量管理方法.pdf

发明内容

针对现有基于水库调度的河流生态流量管理方法的不足，本发明要解决的技术问题是建立一种新的生态流量管理方法，该方法能适合不同的供水保证率要求，同时能提供尽可能优化的生态流量过程，有效的保护河流生态系统。本发明解决这些技术问题所采用的技术方案如下：首先，分别建立生态流量管理基本方案、折衷方案和理想方案。这三种方案分别适合高、中和低供水保证率要求，并分别能维持基本、良好和理想的河流生态系统健康状态。然后，将水库的库容空间分为三个区，从上到下这三个区依次对应生态流量管理理想、折衷和基本方案，从而通过水库的水位将这三个方案耦合为一个统一的方案。当规划的供水保证率高时，水库管理者扩大水库库容最下面的一个区；当水库规划的供水保证率低时，管理者扩大水库库容最上面的一个区，从而该耦合方案可以适合不同的供水保证率，也能提供尽可能优化的生态流量过程。其次，采用水库入流和水库调度曲线相结合的方法指导水库向人类供水，将人类供水方案和以上建立的耦合的生态流量管理方案相结合，构成完整的水库调度方案，指导水库向人类和河流供水。最后，运用遗传算法，以满足规划的供水保证率为基本约束，以减小河流水文情势的扰动为优化目标，对水库库容分区、调度曲线等参数同时进行优化，从而得出最优的基于水库调度的河流生态流量管理方案。其具体的步骤如下：

1.水库入流流量区间分割

采用参数I₁，I₂和I₃将水库入流(R)分为四个区间：R≥I₁，I₁＞R≥I₂，I₂＞R＞I₃，R≤I₃。I₁设为河流平滩流量和水库规划供水量之和；I₂设为水库断面多年日径流量的第75百位数；I₃设为河流最小生态需水量。

2.河流生态流量管理基本方案建立

河流生态流量管理基本方案只提供河流最小生态需水和低于该需水的流量，以保证高的人类用水保证率，具体的方案如下：

a)当水库的入流流量R大于I₃时，河流的生态流量设为河流最小生态需水量；

b)当水库的入流流量R小于或等于I₃时，河流的生态流量设为水库入流量。

3.河流生态流量管理折衷方案建立

河流生态流量管理折衷方案提供河流最小生态需水、低于该需水的流量以及限定数量的高流量水文事件(包括洪水和高脉冲事件)。水库管理者对每个月规定可以下泄的高流量水文事件最大数目[M_a(a＝1，2，...，12)]。M_a的大小取决于规划的供水保证率和河流的生态需求。具体的方案如下：

a)当水库的入流流量R大于或等于I₂时，1)如果已经下泄的高流量水文事件数目小于规划的最大数目M_a，河流的生态流量设为河流入流和规划的供水量之差，同时该生态流量要不大于河流平滩流量；

2)如果已经下泄的高流量水文事件数目等于或大于规划的最大数目M_a，河流的生态流量设为河流最小生态需水量；

b)当水库的入流流量R位于I₂和I₃之间时，河流的生态流量设为河流最小生态需水量；

c)当水库的入流流量R小于或等于I₃时，河流的生态流量设为水库入流量。

4.河流生态流量管理理想方案建立

河流生态流量管理基本方案提供河流最小生态需水、低于该需水的流量以及所有的高流量水文事件。该方案能有效的维持河流的主要生态功能，维持河流生态系统健康，但只能提供较低的人类供水保证率，具体的方案如下：

a)当水库的入流流量R大于或等于I₁时，河流的生态流量设为河流平滩流量；

b)当水库的入流流量R位于I₁和I₂之间时，河流的生态流量设为河流入流和规划的供水量之差；

c)当水库的入流流量R位于I₂和I₃之间时，河流的生态流量设为河流最小生态需水量；

d)当水库的入流流量R小于等于I₃时，河流的生态流量设为水库入流量。

5.河流生态流量管理耦合方案建立

采用两个水库水位H₁和H₂将水库库容分为三个区，这三个区从上到下依次称为理想生态库容区、折衷生态库容区和基本生态库容区。当水库水位位于不同生态库容区时，采取不同的生态流量管理方案，具体如下：

a)当水库水位位于理想生态库容区时，采取河流生态流量管理理想方案；

b)当水库水位位于折衷生态库容区时，采取河流生态流量管理折衷方案；

c)当水库水位位于基本生态库容区时，采取河流生态流量管理基本方案。

6.人类用水供给方案建立

图1表示的是M-5型水库调度曲线，它包括三条曲线：上调度曲线(upper limit curve)，下调度曲线(lower limit curve)和关键调度曲线(critical limit curve)。为了保持水库的防洪功能，上调度曲线可事先通过模拟的方式得到，在本发明中保持不变，另外两条曲线根据实际需要优化确定。这两条曲线都可以通过6个参数来表示，其中两个参数用来表示曲线的上、下水位高度(X₁，X₂：下调度曲线；X₃，X₄：关键调度曲线)，另外四个参数用来表示曲线转折的时间(T₁，T₂，T₃，T₄：下调度曲线；T₅，T₆，T₇，T₈：关键调度曲线)。本发明采用水库入流和水库调度曲线相结合的方法，指导水库向人类供水，具体的方案如下：

a)当水库入流流量R大于或等于I₂时，如果水库水位位于理想生态库容区，或者水库水位位于折衷生态库容区且已经下泄的高流量水文事件数目小于规划的最大数目M_a，此时向人类供给规划的供水量；如果水库水位位于基本生态库容区，或者水库水位位于折衷生态库容区且已经下泄的高流量水文事件数目大于或等于规划的最大数目M_a，此时水库按照以下规则向人类供给供水量：

●当水库水位高于上调度曲线时，通过溢洪道泄流使水库的水位降到上调度曲线水位，并向用水户提供规划的水量；

●当水库水位位于上调度曲线和下调度曲线之间时，水库正常供水；

●当水库水位位于下调度曲线和关键调度曲线之间时，水库供水减小α％；

●当水库水位位于关键调度曲线和死库容线之间时，水库供水减小β％；

●当水库水位低于死库容线时，水库不供水；

其中α和β为参数，并且0＜α＜β＜100。

b)当水库入流流量R小于I₂时，水库按照以下规则向人类供水：

●当水库水位高于上调度曲线时，通过溢洪道泄流使水库的水位降到上调度曲线水位，并向用水户提供规划的水量；

●当水库水位位于上调度曲线和下调度曲线之间时，水库正常供水；

●当水库水位位于下调度曲线和关键调度曲线之间时，水库供水减小α％；

●当水库水位位于关键调度曲线和死库容线之间时，水库供水减小β％；

●当水库水位低于死库容线时，水库不供水。

7.水库调度目标量化及优化参数范围确定

本发明采用变化范围法(Range of Variability Approach，RVA)衡量河流水文情势扰动程度。该方法首先选取了32个具有重要生态意义的水文指标(表1)，并运用公式(1)确定各个指标的水文扰动程度D_m，然后运用公式(2)求出这32个指标的水文情势扰动程度的算术平均值D。

Dm=|Na,m-Nb,mNb,m|×100%---(1)]]>

D=132Σm=132Dm---(2)]]>

式中D_m是第m个水文指标的扰动程度；N_a，m是人为扰动后第m个水文指标的数值落在其目标变化范围的次数；N_b，m是人为扰动前该水文扰动指标的数值落在其目标变化范围的次数；D为水文情势的平均扰动程度。

水库的调度需满足人类所需的供水保证率。水库供水的供水保证率S可通过公式(3)计算。

S=1N×Σb=1NΣa=1TRRabPRab×100%---(3)]]>

其中，N表示水库供水的年数，RR_ab表示第b年中第a天的实际供水量，PR_ab表示第b年中第a天的规划供水量，T表示第b年的天数。

水库调度需要同时考虑人类的需求和河流生态保护的需求。本发明将水库调度的目标量化为：在满足规划的供水保证率的前提下，使河流水文情势扰动值D最小，可用式(4)和(5)表示。根据河流生态需水理论，河流的水文情势扰动越小，生态流量维持的效果越好，河流生态系统越健康。

L＝Min(D)                                                        (4)

s.t.：R≥R₀                                                      (5)

本发明的优化变量包括每个月规定可以下泄的高流量水文事件最大数目[M_a(a＝1，2，...，12)]，水库生态库容分割参数H_k(k＝1，2)以及水库调度曲线的水位高度参数X_i(i＝1，2，3，4)、时间参数T_j(j＝1，2，...，8)。确定M_a的约束的方法为：通过多年自然日径流水文序列，确定各个月历史最大的高流量水文事件发生次数[MT_a(a＝1，2，...，12)]，M_a约束设为0≤M_a≤MT_a。H_k(k＝1，2)、X_i(i＝1，2，3，4)和T_j(j＝1，2，...，8)的约束条件为式(6)-(12)。

MAXlevel≥H₁≥H₂≥MINlevel                                                (6)

MAXlevel＞X₁＞X₂                                                          (7)

MAXlevel＞X₁＞X₃                                                          (8)

X₂＞X₄＞MINlevel                                                          (9)

X₃＞X₄＞MINlevel                                                          (10)

1≤T₁＜T₂＜T₃＜T₄≤36                                                   (11)

1≤T₅＜T₆＜T₇＜T₈≤36                                                     (12)

式中MAXlevel和MINlevel分别为水库最高和最低允许水位。

8水库调度参数优化

采用遗传算法，以河流水文扰动程度D最小化为目标，对每个月规定可以下泄的高流量水文事件最大数目[M_a(a＝1，2，...，12)]，水库生态库容分割参数H_k(k＝1，2)以及水库调度图的水位高度参数X_i(i＝1，2，3，4)、时间参数T_j(j＝1，2，...，8)同时进行优化，按照建立的生态流量管理耦合方案和人类用水供给方案，对不同参数下的水库运行情况进行模拟，得出不同参数下的实际供水保证率R和水文情势扰动D值，进而得出实际供水保证率不小于R₀且能使D值最小的调度参数值。遗传算法中个体的形式为I＝[M₁，...，M₁₂，H₁，H₂，X₁，...X₄，T₁，...，T₈]。

本发明具有两个优点：1)适合较宽的水库供水保证率需求。当规划的供水保证率高时，水库管理者可以扩大基本生态库容区，减小生态流量的供给；当规划的供水保证率低时，水库管理者可以扩大理想生态库容区，增加生态流量的供给。随着三个生态库容区大小的变化，水库可以实现不同的供水保证率，从而提高了本发明的广泛适用性。2)更能实现河流生态系统的保护。现有生态流量管理方案相当于本发明中河流生态流量管理基本方案的第一条(只提供最小生态需水)，本发明根据水库入流量和水库水位的不同，对生态流量管理方案做了大量的扩展，更有利于满足河流生态流量的需求，保持河流生态系统健康；另外，受规划的供水保证率的约束，水库难以长期采用生态调度理想方案或折衷方案。采用本发明建立的生态流量管理耦合方案，在水库水量充足时，水库启动生态调度折衷或者理想方案，从而比只采用基本方案或折衷方案更能保护河流生态系统健康。