一种面向环境保护的单一供水水库的生态库容确定方法
【专利摘要】本发明涉及一种面向环境保护的单一供水水库的生态库容确定方法，属于环境保护与资源综合利用【技术领域】。首先根据河流管理部门规定的河流生态流量要求，确定水库生态用水调度的优化目标，结合优化变量与约束条件，建立水库生态用水调度优化模型，构造多目标的权重集合，将多目标化成单目标进行求解；然后根据经济社会和生态缺水率的分析，确定优化供水方案；最后根据优化的经济社会和生态供水量，按比例进行生态库容的计算。本发明弥补了传统供水水库库容设置未能考虑生态供水的弊端，符合未来水库进行生态用水调度的需求；根据优化模型结果计算，能够得到更为合理的生态库容，具有理论意义明确、可操作性强、容易得到实际应用的优点。
【专利说明】
【技术领域】
[0001] 本发明涉及一种面向环境保护的单一供水水库的生态库容确定方法，属于环境保 护与资源综合利用【技术领域】。 一种面向环境保护的单一供水水库的生态库容确定方法
【背景技术】
[0002] 传统水库的运行与调度，都是以防洪、供水、发电等社会功能和经济功能极大化为 目标。这些水库的调度运用，改变了河流水文的自然过程及规律，影响到鱼类等水生生物的 栖息生存，对河流生态系统造成威胁。随着生态环境问题的日益突出，河流生态环境需水 也就成为水库调度中越来越重要的任务，一些学者纷纷提出了生态调度的概念（董哲仁， 2007)，并将生态流量与水库调度结合，发展了考虑生态需水量的水库调度方案和规则（杨 志峰，2010 ;戴会超，2012 ;刘攀，2013)。
[0003] 这些方法都是从水库的传统库容和调度方法角度出发，按下游河道内生态需水要 求约束下泄生态水量，在具体实施中，容易产生经济用水和生态用水的矛盾，难以协调。近 年有学者提出了水库生态库容概念（廖四辉，2011)，并讨论了其对保障生态流量的必要性 (雍婷，2013 ;吕孙云，2013)，其确定方法是对设置的生态限制供水线进行优化，仍是在传 统库容和调度方法上的约束调度，容易使枯水季节只有生态供水，而人类生活生产需水缺 口较大，矛盾依然突出。此外，传统的水库设计运行中库容设置仅有死库容、兴利库容和防 洪库容确定方法（顾圣平，2009)，根据现实需要，有必要在传统库容的基础上专设生态库 容，发展生态库容确定方法。目前，在供水水库中，还未有生态库容的确定计算方法，水库运 行中还难以明晰和协调人类社会经济、防洪和生态目标的水量调度。
【发明内容】
[0004] 本发明的目的是提出一种面向环境保护的单一供水水库的生态库容确定方法，结 合当前生态流量需求和优化模型工具，建立已有单一供水水库的生态库容确定方法，使供 水水库在发挥河川径流调蓄和经济社会供水作用的同时，动用专门库容进行生态供水，保 障河道内环境流量，实施环境保护和生态友好的生态调度。
[0005] 本发明提出的面向环境保护的单一供水水库的生态库容确定方法，其流程框图如 图1所示，包括以下步骤：
[0006] (1)根据河流管理部门的要求，分别确定供水水库下游河流的生态控制断面各月 份的生态需水量其中k为月份，k = 1，2,…，12,将各月份的生态需水量》;作为生 态流量目标；
[0007] (2)建立一个水库生态用水调度优化模型，优化模型的目标函数为：
[0008] Μ?ηΖ = + + ^W：ihlck) + ±w4W^k /=1 /.=1
[0009] 上式中为供水水库的第i个供水区的生活缺水量，》为第i个供水区 的工业及城市缺水量，PLrf为第i个供水区的农业缺水量，为供水水库下游河流的 第j个生态控制断面的生态缺水量，Wp w2、w3、w4为权重系数，通过随机或人工生成，满足 = 1，η为由供水水库供水的供水区总数，m为水库下游河流的生态控制断面总 数；
[0010] (3)设定关于供水水库、供水水库所在河流、供水区、水库下游河流的生态控制断 面、汇流节点和渠道的优化模型的约束条件，如下：
[0011] 供水水库：
[0012] V (t+1) = V (t) +ffresin (t) -Wresout (t) -Wres_sup (t) -Lres (t)
[0013] Vds^ V(t) ^ Vfflax
[0014] Wres-out (t) < Qres max
[0015] Wres_sup(t) ^
[0016] 供水水库所在河流：
[0017] Wriv in (t) = ffriv out (t) +Lriv (t)
[0018] ffriv in(t) ^ Qriv max
[0019] Wriv out(t)彡 Qriv max
[0020] 供水区：
[0021 ] wsuP(t) = ffcon(t)+wret(t)
[0022] Wsup(t) = Wdem(t)-fflack(t)
[0023] 生态控制断面：
[0024] Weco (t) = ffedem (t) -ffelack (t) +ffeinc (t)
[0025] 汇流节点： N
[0026] WcdYj￥:m,」n(t) /-1
[0027] 渠道：
[0028] Wsup(t) = (l-a)Xffres_sup(t)
[0029] 其中，t为优化模型的时段指标，t = 1，2，···，T，T为模型计算的总时间，V(t)、 V(t+1)分别为第t时段和第t+1时段的供水水库的蓄水量，WMS in(t)、WMS 、WMS sup(t)、 L_(t)分别为第t时段的供水水库的入流量、下泄流量、向供水区的供水量和损失量，Vds为 供水水库的下限库容，为供水水库的上限库容，Q MS _为供水水库的下泄能力，Qd_为 供水水库的引水能力，WHv in(t)、WHv ^(t)、Liv(t)分别为第t时段的供水水库的上游河 流入流量、下游河流出流量和所在河流的损失量，Q&_ max为下游河流的过流能力，Wsup(t)、 Wcon(t)、Wret (t)、Wdem(t)、Wlack(t)分别为第t时段供水区得到的供水量、耗水量、回归下游河 流的水量、需水量和缺水量，W_(t) 、W eladi(t)、Weinc;(t)分别为第t时段供水水库下 游河流的生态控制断面的生态供水量、需水量、缺水量和加大量，W_ _α)为第t时段汇流 节点的出流量，表示第t时段汇流节点的第1个分支入流量，N为流入同一汇流节 点的分支总数，α为渠道的损失系数，取值范围为〇-1的实数，上述约束条件中，供水水库 的下限库容、供水水库的上限库容、供水水库的下泄能力、供水水库的引水能力、下游河流 的过流能力和渠道的损失系数由河流管理部门提供；
[0030] (4)采用线性规划方法求解，以10天为时间步长，根据30年水文历史资料中的供 水水库的上游河流入流量W& in(t)，求解上述步骤（2)和步骤（3)构成的优化模型，得到优 化的供水水库向供水区的供水量W Ms sup(t)和供水水库的生态供水量W_(t);
[0031] (5)随机或人工生成多组权重系数，满足A+W2+W3+W4 = 1，得到权重集 合，，重复步骤（4)，计算得到多组优化的供水水库向供水区的供水量W_sup(t)和供水水 库的生态供水量1。。(〇,形成一个由多组优化的供水水库向供水区的供水和 供水水库的生态供水量W_ (t)组成的非劣解集f (WMS_SUP (t)，W_ (t));
[0032] (6)根据上述多组优化的供水水库向供水区的供水量WMs sup(t)和供水水库的生 态供水量W_ (t)，通过下式计算人类社会经济缺水率Rs ladt和生态缺水率R_ ladt :
【权利要求】
1. 一种面向环境保护的单一供水水库的生态库容确定方法，其特征在于该计算方法包 括以下各步骤： (1) 根据河流管理部门的要求，分别确定供水水库下游河流的生态控制断面各月份的 生态需水量其中k为月份，k = 1，2，…，12,将各月份的生态需水量d,作为生态流 量目标； (2) 建立一个水库生态用水调度优化模型，优化模型的目标函数为： /=1 /=1 上式中为供水水库的第i个供水区的生活缺水量，为第i个供水区的工业 及城市缺水量，为第i个供水区的农业缺水量，KLr为供水水库下游河流的第j个生 态控制断面的生态缺水量，为权重系数，通过随机或人工生成，满足Wi+WfwfWj =1，η为由供水水库供水的供水区总数，m为水库下游河流的生态控制断面总数； (3) 设定关于供水水库、供水水库所在河流、供水区、水库下游河流的生态控制断面、汇 流节点和渠道的优化模型的约束条件，如下： 供水水库： V (t+1) = V (t) +ffres in (t) -Wres out (t) -ffres_sup (t) -Lres (t) Vds ^ V(t) ^ Vfflax Wres-out (t) Qres-max Wres-sup (t) Qdmax 供水水库所在河流： Wrivjn(t) = Wriv out (t)+Lriv (t) Wriv-in ⑴彡 Q riv-max Wfiv-out (t) Qriv-max 供水区： WsuP(t) = Wcon(t)+ffret(t) Wsup (t) = Wdem (t) -fflack (t) 生态控制断面： weco (t) = ffedem (t) -welack (t) +ffeinc (t) 汇流节点： υ)=?υ /=1 渠道： wsup ⑴=(l-a)XWres-sup(t) 其中，t为优化模型的时段指标，t = 1，2，一,111为模型计算的总时间，V(t)、V(t+l) 分别为第t时段和第t+1时段的供水水库的蓄水量，WMS in(t)、WMS wt(t)、WMS sup(t)、Lres(t) 分别为第t时段的供水水库的入流量、下泄流量、向供水区的供水量和损失量，V ds为供水 水库的下限库容，为供水水库的上限库容，Qras _为供水水库的下泄能力，为供水 水库的引水能力，WHv in(t)、WHv Mt(t)、L,iv(t)分别为第t时段的供水水库的上游河流入流 量、下游河流出流量和所在河流的损失量，Q&_max为下游河流的过流能力，Wsup(t)、W_(t)、 WMt(t)、Wd"(t)、Wladi(t)分别为第t时段供水区得到的供水量、耗水量、回归下游河流的水 量、需水量和缺水量，W_(t) 、W eladi(t)、Weinc;(t)分别为第t时段供水水库下游河流 的生态控制断面的生态供水量、需水量、缺水量和加大量，w__α)为第t时段汇流节点的 出流量，ML」,#)表示第t时段汇流节点的第1个分支入流量，N为流入同一汇流节点的分 支总数，α为渠道的损失系数，取值范围为〇-1的实数，上述约束条件中，供水水库的下限 库容、供水水库的上限库容、供水水库的下泄能力、供水水库的引水能力、下游河流的过流 能力和渠道的损失系数由河流管理部门提供； (4) 采用线性规划方法求解，以10天为时间步长，根据30年水文历史资料中的供水水 库的上游河流入流量W& in(t)，求解上述步骤（2)和步骤（3)构成的优化模型，得到优化的 供水水库向供水区的供水量W Ms sup(t)和供水水库的生态供水量W_(t); (5) 随机或人工生成多组权重系数％、w2、w3、w4，满足Wi+w 2+w3+w4 = 1，得到权重集合 ?Ω，重复步骤（4)，计算得到多组优化的供水水库向供水区的供水量W_ sup(t)和供水水库 的生态供水量W_ (t)，形成一个由多组优化的供水水库向供水区的供水量sup (t)和供 水水库的生态供水量W_ (t)组成的非劣解集f (WMS_SUP (t)，W_ (t)); (6) 根据上述多组优化的供水水库向供水区的供水量WMs sup(t)和供水水库的生态供 水量W_ (t)，通过下式计算人类社会经济缺水率Rs ladt和生态缺水率R_ ladt :
当人类社会经济缺水率Rs_ladt与生态缺水率R__ladt的值相等时（或某一个特定的协调 值），从非劣解集W" (WMS sup (t)，W_ (t))中检索出与人类社会经济缺水率Rs ladt和生态缺水 率&。。」^相等条件下（或某一个特定的协调值）相应的Ls_sup(t)和^。。(〇 ; (7) 根据上述_和W_，利用下式计算生态库容入。。：
其中，Vactive为从河流管理部门获取的当前供水水库的兴利库容。
【文档编号】E02B1/02GK104047257SQ201410283712
【公开日】2014年9月17日 申请日期:2014年6月23日 优先权日:2014年6月23日
【发明者】王忠静, 朱金峰, 尚文绣, 赵建世, 郑航, 黄草 申请人:清华大学