CN110348600B - 一种水库群防洪效益最大的库容分配方法 - Google Patents

一种水库群防洪效益最大的库容分配方法 Download PDF

Info

Publication number
CN110348600B
CN110348600B CN201910478191.5A CN201910478191A CN110348600B CN 110348600 B CN110348600 B CN 110348600B CN 201910478191 A CN201910478191 A CN 201910478191A CN 110348600 B CN110348600 B CN 110348600B
Authority
CN
China
Prior art keywords
reservoir
flood control
flood
storage capacity
maximum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910478191.5A
Other languages
English (en)
Other versions
CN110348600A (zh
Inventor
周建中
周华艳
卢程伟
沈柯言
陈潇
杨鑫
朱锦干
曾德晶
金倩芳
曾昱
顿晓晗
李力
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201910478191.5A priority Critical patent/CN110348600B/zh
Publication of CN110348600A publication Critical patent/CN110348600A/zh
Application granted granted Critical
Publication of CN110348600B publication Critical patent/CN110348600B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种防洪库容效益最大化的库容分配方法，属于水文学中的水库防洪调度库容分配利用领域，包括以下步骤：(1)建立各水库对下游防洪控制点的防洪优化调度模型；(2)根据防洪调度模型计算的水库消耗防洪库容和下游削减洪峰的对应关系；(3)进行函数拟合得水库防洪库容效益函数；(4)建立基于水库群防洪效益最大的联合防洪调度模型并求解库容分配方案。本发明考虑了洪水过程、河道演进和水库自身特性的影响，通过模拟防洪调度建立了上游水库群对下游防洪控制点的防洪库容效益关系，基于水库群总的防洪效益最大，求解在保证下游防洪控制站点安全时上游水库的合理分配，分配方法直观，操作简单，适用性强。

Description

一种水库群防洪效益最大的库容分配方法
技术领域
本发明属于水文学中的水库防洪调度库容分配利用领域，更具体地，涉及一种水库群防洪效益最大的库容分配方法。
背景技术
水库防洪调度是水库对入库洪水进行拦蓄、削减洪峰、补偿调度等对下游防洪安全进行的出入库调度操作。对于一场洪水过程，当洪水量级不大，水库可以存储部分或全部的洪水保证下游的防洪安全，而量级较大时即使利用了水库的所有防洪库容，下游亦会出现超安全流量的洪水即产生洪灾。
对于一个防洪安全控制点，常有多个串并联耦合的水库等水库工程共同参与防洪调度来保证流域的防洪安全，且不同的水库由于水库自身的库容大小、距离防洪控制点的远近、防洪任务的不同会产生不同的防洪效益和调度规则。因此，在流域水库群联合防洪调度中，如何解析水库各自的防洪任务并合理的进行水库各自防洪库容分配，使得水库防洪库容充分利用在保证下游防洪控制点安全的条件下达到最大的防洪效益是重要的研究内容。
现有的水库群联合防洪调度库容分配利用主要通过水库群剩余防洪库容最大模型或引入权重系数实现不同防洪控制点的区分来进行库容分配，这些方法在实际使用存在以下缺陷：
(1)建立联合防洪优化调度模型时，约束耦合条件复杂，求解难度大，难以运用于实际问题中。
(2)现有研究一般将库容分配转化为求解的约束或通过优化调度的结果来分析出库容分配方案，不够直观。
(3)现有的库容分配方法在使用中只考虑了洪水形状的影响，暂时忽略了河道传播对洪水的坦化作用，但河道演进对洪水形状过程影响较大，进而影响库容分配结果。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种防洪库容效益最大化的库容分配方法，其目的在于综合考虑水库自身特性、河道演进特性和洪水形状过程的影响因素确定水库群防洪库容效益最大的各水库库容，优化水库群防洪库容分配方案，提高水库群流域防洪的整体效益，有效降低水库群防洪风险。
为实现上述目的，本发明提供了一种水库群防洪效益最大的库容分配方法，其特征在于，包括如下步骤：
(S1)以水库消耗防洪库容最小为目标，建立单个水库对下游防洪控制点的防洪优化调度模型；
(S2)根据所述防洪优化调度模型多次模拟调度计算得到单水库消耗防洪库容和下游防洪控制点削减洪峰流量的多组数据，利用该多组数据拟合得到各水库防洪库容效益的函数关系；所述水库防洪库容效益FSE是指单位水库的防洪库容产生的对应防洪控制站点的洪峰削减量值，表达式为：
FSE＝ΔQ/Vf (1)
式中，ΔQ为下游防洪控制点削减洪峰流量，Vf为水库消耗的防洪库容；
(S3)以水库群防洪库容效益最大为目标建立基于防洪库容效益最大的联合防洪调度模型；
(S4)利用步骤(S2)得到的各水库的防洪库容效益的函数关系，对建立的联合防洪调度模型进行求解；根据求解的计算结果得到防洪库容效益最大时各水库库容的优化分配方案。
优选地，所述单个水库对下游防洪控制点的防洪优化调度模型如下：
Figure BDA0002082962740000031
式中，t为时段，T为调度周期，Tm为取最大消耗防洪库容的时段点，It为水库t时段入库流量，Qt为水库t时段出库流量，Δt为调度时段间隔长度，Vf为水库消耗的防洪库容；
约束条件如下：
防洪库容约束：Vf≤V总 (3)
河道演进约束：Qt'＝C0Qt+C1Qt-1+C2Q't-1 (4)
削减洪峰流量：ΔQ＝Qpeak-qt-Qt' (5)
水库最大最小下泄：Qmin≤Qt≤Qmax (6)式中，V总为水库最大可用防洪库容，Q't为t时段水库出库流量演进到防洪控制点的流量，Qt为t时段出库流量，Qt-1为t-1时段出库流量，Q't-1为t-1时段出库流量演进到防洪控制点的流量，C0,C1,C2为河道演进系数，ΔQ为水库削减洪峰流量，Qpeak为防洪控制点最大洪峰流量，qt为区间流量，Qmax为水库最大下泄流量，Qmin为水库最小下泄流量。
优选地，根据所述防洪优化调度模型多次模拟调度计算得到多组单水库消耗防洪库容V和下游防洪控制点削减洪峰流量ΔQ的多组数据，采用多项式V(ΔQ)＝a1ΔQn+a2ΔQn -1+…+anΔQ+b对该多组数据进行拟合，得到防洪库容效益的函数关系；其中a1,a2,…,an,b为多项式中的待求解参数。
优选地，所述基于防洪库容效益最大的联合防洪调度模型如下：
Figure BDA0002082962740000032
式中，Z为目标指代量，n为水库群中水库的总个数，Vf,i为水库i在保障下游防洪控制点安全时消耗的防洪库容，Vi为水库i消耗的防洪库容，t为调度时段，T为调度周期，Tm为调度单水库消耗最大库容时段点，Ii,t为水库i在t时段的入库流量，Qi,t为水库i在t时段的出库流量，Δt为时段间隔；
约束条件如下：
防洪库容约束：Vf≤V总 (3)
河道演进约束：Qt'＝C0Qt+C1Qt-1+C2Q't-1 (4)
削减洪峰流量：ΔQ＝Qpeak-qt-Qt' (5)
水库最大最小下泄：Qmin≤Qt≤Qmax (6)
式中，V总为水库最大可用防洪库容，Q't为t时段水库出库流量演进到防洪控制点的流量，Qt为t时段出库流量，Qt-1为t-1时段出库流量，Q't-1为t-1时段出库流量演进到防洪控制点的流量，C0,C1,C2为河道演进系数，ΔQ为水库削减洪峰流量，Qpeak为防洪控制点最大洪峰流量，qt为区间流量，Qmax为水库最大下泄流量，Qmin为水库最小下泄流量。
优选地，根据所述基于防洪库容效益最大的联合防洪调度模型，对下游最大洪峰削减ΔQ，使上游水库群消耗的防洪库容最小且保证下游防洪控制点的防洪安全，抽象出数学模型如下:
Figure BDA0002082962740000041
式中，V为水库群消耗的库容，V1为一号水库使下游防洪控制点最大洪峰削减ΔQ1所消耗的库容，V2为二号水库使下游防洪控制点最大洪峰削减ΔQ2所消耗的库容，...，Vn为n号水库使下游防洪控制点最大洪峰削减ΔQn所消耗的防洪库容，n个水库共同使下游防洪控制点最大洪峰削减ΔQ，使下游防洪控制点最大洪峰流量不超过安全流量；
利用Lagrange乘子法，构造辅助函数，
F＝V1+V2+…+Vn+λ(ΔQ1+ΔQ2+…+ΔQn) (9)
式中，F为辅助函数，λ为Lagrange乘子，其取极值条件为
Figure BDA0002082962740000042
联立公式(8)、(9)和(10)，利用步骤(S2)得到的各水库的防洪库容效益的函数关系，求解得到各水库防洪库容效益关系曲线微增率相等时的各水库消耗的防洪库容V1,V2,…,Vn，即为防洪库容效益最大时各水库库容的优化分配方案。
优选地，由公式(1)可知，水库群下游防洪控制点削减指定的洪峰流量ΔQ时，水库群消耗防洪库容最小时水库群防洪库容效益最大；
将水库的防洪库容效益函数转换表达为消耗库容对下游防洪控制点削减的洪峰量值的偏导数相等时求得各水库消耗的防洪库容V1,V2,…,Vn，此时各水库的库容微增率相等，水库群消耗防洪库容最小，即为水库群防洪库容效益最大的各水库库容的优化分配方案。
优选地，建立所述单个水库对下游防洪控制点的防洪优化调度模型之前，先分析并选择典型入库洪水。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
(1)本发明提出的基于防洪库容效益最大化的库容分配方法，充分考虑了水库自身特性、河道传播特性和洪水过程特性等各种因素的影响，能够有效避免影响因素不全面对库容分配的影响，得到的库容分配方案更加合理准确。
(2)本发明提出的基于防洪库容效益最大化的库容分配方法，在建立了防洪库容效益曲线的基础上进行库容分配求解，可运用于大规模水库群联合防洪调度求解、水库预报防洪调度及水库防洪库容规划设计等各防洪调度领域中，易于实际运用，操作简单方便。
附图说明
图1是本发明较佳实施例中一种防洪库容效益最大化的库容分配方法流程图；
图2是本发明较佳实施例中水库库容对应下游控制点削减峰量值关系曲线；
图3(a)是本发明较佳实施例中1998年上游溪向梯级水库和下游防洪控制点削峰量值关系曲线；
图3(b)是本发明较佳实施例中1998年上游瀑布沟水库和下游防洪控制点削峰量值关系曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供一种水库群防洪效益最大化的库容分配方法，综合考虑水库自身特性、河道演进特性和洪水形状过程等影响因素，提供了一种简单直观操作方便的水库群在保障下游防洪控制点安全时的库容分配途径。包括以下步骤：
(1)以防洪控制站点为控制的设计洪水分析，确定具有代表性的典型洪水过程为研究洪水。入库洪水过程的不同会使水库对防洪控制站点产生不同的防洪效益，因此需针对不同的入库洪水过程进行具体分析，需首先收集防洪控制站点、研究水库入库洪水的长系列洪水资料，对该资料洪水进行辨析和聚类，确定对下游危害较大、不易控制的大洪水典型过程为研究洪水对象。所述的典型入库洪水，一般按照洪峰和洪量划分，选取峰高量大的历史洪水位典型洪水。
(2)研究分析水库的防洪库容效益关系。综合考虑水库自身特性、洪水过程和河道演进影响，建立各研究水库对应防洪控制站点的模拟调度模型，计算消耗不同防洪库容下对应的防洪控制站点的洪峰削减量，得到水库的防洪库容消耗和下游防洪控制站点削峰量之间的关系点，对这些数据关系进行函数拟合，得到水库的防洪库容效益函数。
水库防洪库容效益FSE是指单位水库的防洪库容产生的对应防洪控制站点的洪峰削减量值，表达式为：
FSE＝ΔQ/Vf (1)
式中，ΔQ为下游防洪控制点削减洪峰流量，Vf为水库消耗的防洪库容；
1)模型建立：
目标函数为：
Figure BDA0002082962740000071
式中，It、Qt分别为水库t时段初入库流量和出库流量，Δt为调度时段间隔长度，Vf为水库消耗的防洪库容。
2)约束条件为：
①防洪库容约束：Vf≤V总 (3)
②河道演进约束：Qt'＝C0Qt+C1Qt-1+C2Q't-1 (4)
③削减洪峰流量：ΔQ＝Qpeak-qt-Qt' (5)
④水库最大最小下泄：Qmin≤Qt≤Qmax (6)
式中，V总为水库最大可用防洪库容，Qt-1为t-1时段出库流量，Q't-1为t-1时段出库流量演进到防洪控制点的流量，Qt为t时段出库流量，Qt'为水库出库流量演进到防洪控制点的流量，qt为区间流量，Qpeak为防洪控制点最大洪峰流量，C0,C1,C2为河道演进系数，Qmax,Qmin为水库最大最小下泄流量。
(3)防洪库容效益FSE函数计算
根据上述建立水库对下游防洪控制点的调度模型，可得到水库消耗的不同的防洪库容对下游的防洪控制点的削峰量的对应关系。对这系列数据关系，可采用多项式V(ΔQ)＝a1ΔQn+a2ΔQn-1+…+anΔQ+b进行拟合，其中a1,a2,…,an,b为多项式中的待求解参数，因大部分的函数关系都可以通过多项式进行展开，且可对拟合的函数进行检验，设定达到评定标准的为可适用的函数表达。
(4)水库群防洪库容分配计算
对选取的每一种洪水类型，通过模拟计算得到了各个水库的防洪库容效益曲线，为使防洪库容效益最大化，进行防洪库容的合理分配，以保证下游防洪控制点安全为前提，构建水库群防洪库容效益最大为目标的模型。既有模型为：
目标函数：
Figure BDA0002082962740000081
式中，Z为目标指代量，n为水库群数量，Vf,i为水库i在保障下游防洪安全时消耗的防洪库容，Vi为水库消耗的防洪库容，t为调度时段，T为调度周期，Tm为调度单水库消耗最大库容时段点，Ii,t为i水库t时段的入库流量，Qi,t为i水库t时段的出库流量，Δt为时段间隔。
约束条件：同上述表达。
根据上述模型可知，要使上游水库消耗的防洪库容最小，且保证下游防洪控制点的防洪安全，即在下游最大洪峰流量削减到安全流量下的情况下，上游水库库容消耗最小，而对下游削减固定的洪水峰量，上游参与调度的水库不同、各水库分别储存的洪水不同都会产生不同的水库库容消耗。由此可抽象出简单的数学模型：
Figure BDA0002082962740000082
式中，V为水库群消耗的库容，V1为一号水库使下游防洪控制点最大洪峰削减ΔQ1所消耗的库容V2为二号水库使下游防洪控制点最大洪峰削减ΔQ2所消耗的库容，Vn为n号水库使下游防洪控制点最大洪峰削减ΔQn所消耗的防洪库容，n个水库共同达到是下游防洪控制点最大洪峰流量不超过安全流量，即需使最大洪峰削减ΔQ。则可利用Lagrange乘子法，构造一个辅助函数，
F＝V1+V2+…+Vn+λ(ΔQ1+ΔQ2+…+ΔQn) (9)
式中，F为辅助函数，λ为Lagrange乘子，其取极值条件为
Figure BDA0002082962740000091
即当各水库防洪库容效益关系曲线微增率相等时为防洪库容效益最大时各水库库容分配取值点。
本发明考虑了洪水过程、河道演进和水库自身特性的影响，创新性地定义了单位库容对下游防洪控制站点削减洪峰流量作用的防洪库容效益概念，通过模拟防洪调度建立了上游水库群对下游防洪控制点的防洪库容效益关系，基于水库群总的防洪效益最大，求解在保证下游防洪控制站点安全时上游水库的合理分配，此方法较传统库容分配方法的考虑的影响因素更加全面，分配方法更加直观，操作更加简单，适用性强。
以下结合附图和实例对本发明提供的防洪库容效益最大化的库容分配方法进一步详细说明。
如附图1所示，本发明包括以下流程：
(1)选取典型入库洪水；
(2)分别建立各水库对下游防洪控制站点的防洪优化调度模型；
(3)根据单库对下游防洪控制点的调度模型，计算防洪库容效益函数关系；
(4)在求解各水库对下游防洪控制点防洪库容效益关系的基础上，建立以水库防洪库容效益最大的联合防洪调度模型，并求解库容分配结果。
以长江流域川江河段水库群(金沙江下游溪洛渡、向家坝水库和岷江瀑布沟水库)为研究对象，综合考虑水库自身特性、河道演进特性和洪水过程影响罂粟，求解在保证下游宜宾城市(防洪控制点为李庄)达到50年一遇防洪标准(洪峰流量不超过51000m3/s)时上游水库群库容分配方案。限于研究资料和洪水类型的分析缺乏，本文选取1998年以李庄为控制的50年一遇的设计洪水为例进行分析，分别建立溪向梯级对李庄的防洪调度模型和瀑布沟水库对李庄的防洪调度模型，并求解对应防洪库容效益关系。
建立溪向梯级和瀑布沟水库以防洪库容效益最大为目标的联合防洪调度模型，要使上游水库群消耗的防洪库容最小，且保证下游防洪控制点的防洪安全，即在下游最大洪峰流量削减到安全流量下的情况下，上游水库库容消耗最小，而对下游削减固定的洪水峰量，上游参与调度的水库不同、各水库分别储存的洪水不同都会产生不同的水库库容消耗。由此可抽象出简单的数学模型：
Figure BDA0002082962740000101
式中，V为水库群消耗的库容，这里简单的考虑两个水库便于表达，V1为一号水库使下游防洪控制点最大洪峰削减ΔQ1所消耗的库容，V2为二号水库使下游防洪控制点最大洪峰削减ΔQ2所消耗的库容，两者共同使下游防洪控制点最大洪峰流量不超过安全流量，即需使最大洪峰削减ΔQ。如附图2所示，当一号水库削减ΔQ1时，二号水库削减ΔQ2时分别位于图线上的a、b点，若一号水库增加削减量，则二号水库需减少对应量以保证下游防洪控制点最大洪峰流量刚好不超过安全流量。根据曲线特征可知，随着削减峰量的增加，曲线的斜率增加，若改变两水库间的削峰量分配可使得总的消耗库容减少则是需要的。即当dV1/dΔQ1＝dV2/dΔQ2有dV＝dV1+dV2＝0，若继续增加一号水库削减的峰量值，并减少二号水库的峰量值，由于dV1/dΔQ1＞dV2/dΔQ2会使得dV＝dV1+dV2＞0，消耗的总库容会增加。综上所述：在确定下游防洪控制点削减的峰量值时，水库群消耗防洪库容最小的条件是两个水库的库容微增率相等，因此将水库的防洪库容效益函数转换表达为消耗库容对下游削峰量值的偏导数相等即为合理分配的库容点。
拟合得1998年防洪库容效益关系曲线如附图3(a)和图3(b)所示。其中溪向梯级防洪库容效益曲线相关系数为0.9915，瀑布沟水库防洪库容效益关系曲线相关系数为0.9978，相关系数趋近于1，相关性强，拟合函数关系可靠。由此建立确保李庄站点防洪安全的防洪库容效益最大化联合防洪调度模型，利用前述求解方法得到库容分配结果。考虑相同约束条件，建立联合优化调度模型，将两种方法所得结果对比如下表所示。
表1 1998年防洪效益库容分配与联合防洪优化调度结果对比表
Figure BDA0002082962740000111
由表1可知在确保防洪控制站点李庄安全的条件下得到的溪向梯级和瀑布沟水库的库容分配和联合防洪优化调度所得的库容分配结果基本一致，总库容相差0.01亿立方米，相对误差5％，表明所提基于防洪库容效益最大化的库容分配方法合理有效。且本发明提供的方法也可运用于大规模水库群联合防洪调度时的库容分配求解，简化了复杂水库群联合防洪调度库容分配求解，由于该方法考虑洪水形状过程、水库自身特性及河道传播特性等影响因素，因此所得库容分配方案相对可靠，准确度高。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种水库群防洪效益最大的库容分配方法，其特征在于，包括如下步骤：
(S1)以水库消耗防洪库容最小为目标，建立单个水库对下游防洪控制点的防洪优化调度模型；
(S2)根据所述防洪优化调度模型多次模拟调度计算得到单水库消耗防洪库容和下游防洪控制点削减洪峰流量的多组数据，利用该多组数据拟合得到各水库防洪库容效益的函数关系；所述水库防洪库容效益FSE是指单位水库的防洪库容产生的对应防洪控制站点的洪峰削减量值，表达式为：
FSE＝ΔQ/Vf (1)
式中，ΔQ为下游防洪控制点削减洪峰流量，Vf为水库消耗的防洪库容；
(S3)以水库群防洪库容效益最大为目标建立基于防洪库容效益最大的联合防洪调度模型；
(S4)利用步骤(S2)得到的各水库的防洪库容效益的函数关系，对建立的联合防洪调度模型进行求解；根据求解的计算结果得到防洪库容效益最大时各水库库容的优化分配方案。
2.根据权利要求1所述的一种水库群防洪效益最大的库容分配方法，其特征在于，所述单个水库对下游防洪控制点的防洪优化调度模型如下：
Figure FDA0002082962730000011
式中，t为时段，T为调度周期，Tm为取最大消耗防洪库容的时段点，It为水库t时段入库流量、Qt为水库t时段出库流量，Δt为调度时段间隔长度，Vf为水库消耗的防洪库容；
约束条件如下：
防洪库容约束：Vf≤V总 (3)
河道演进约束：Qt'＝C0Qt+C1Qt-1+C2Q't-1 (4)
削减洪峰流量：ΔQ＝Qpeak-qt-Qt' (5)
水库最大最小下泄：Qmin≤Qt≤Qmax (6)
式中，V总为水库最大可用防洪库容，Q't为t时段水库出库流量演进到防洪控制点的流量，Qt为t时段出库流量，Qt-1为t-1时段出库流量，Q't-1为t-1时段出库流量演进到防洪控制点的流量，C0,C1,C2为河道演进系数，ΔQ为水库削减洪峰流量，Qpeak为防洪控制点最大洪峰流量，qt为区间流量，Qmax为水库最大下泄流量，Qmin为水库最小下泄流量。
3.根据权利要求1或2所述的一种水库群防洪效益最大的库容分配方法，其特征在于，根据所述防洪优化调度模型多次模拟调度计算得到多组单水库消耗防洪库容V和下游防洪控制点削减洪峰流量ΔQ的多组数据，采用多项式V(ΔQ)＝a1ΔQn+a2ΔQn-1+…+anΔQ+b对该多组数据进行拟合，得到防洪库容效益的函数关系。
4.根据权利要求1或2所述的一种水库群防洪效益最大的库容分配方法，其特征在于，所述基于防洪库容效益最大的联合防洪调度模型如下：
Figure FDA0002082962730000021
式中，Z为目标指代量，n为水库群中水库的总个数，Vf,i为水库i在保障下游防洪控制点安全时消耗的防洪库容，Vi为水库i消耗的防洪库容，t为调度时段，T为调度周期，Tm为调度单水库消耗最大库容时段点，Ii,t为水库i在t时段的入库流量，Qi,t为水库i在t时段的出库流量，Δt为时段间隔；
约束条件如下：
防洪库容约束：Vf≤V总 (3)
河道演进约束：Qt'＝C0Qt+C1Qt-1+C2Q't-1 (4)
削减洪峰流量：ΔQ＝Qpeak-qt-Qt' (5)
水库最大最小下泄：Qmin≤Qt≤Qmax (6)式中，V总为水库最大可用防洪库容，Q't为t时段水库出库流量演进到防洪控制点的流量，Qt为t时段出库流量，Qt-1为t-1时段出库流量，Q't-1为t-1时段出库流量演进到防洪控制点的流量，C0,C1,C2为河道演进系数，ΔQ为水库削减洪峰流量，Qpeak为防洪控制点最大洪峰流量，qt为区间流量，Qmax为水库最大下泄流量，Qmin为水库最小下泄流量。
5.根据权利要求4所述的一种水库群防洪效益最大的库容分配方法，其特征在于，根据所述基于防洪库容效益最大的联合防洪调度模型，对下游最大洪峰削减ΔQ，使上游水库群消耗的防洪库容最小且保证下游防洪控制点的防洪安全，抽象出数学模型如下:
Figure FDA0002082962730000031
式中，V为水库群消耗的库容，V1为一号水库使下游防洪控制点最大洪峰削减ΔQ1所消耗的库容，V2为二号水库使下游防洪控制点最大洪峰削减ΔQ2所消耗的库容，...，Vn为n号水库使下游防洪控制点最大洪峰削减ΔQn所消耗的防洪库容，n个水库共同使下游防洪控制点最大洪峰削减ΔQ，使下游防洪控制点最大洪峰流量不超过安全流量；
利用Lagrange乘子法，构造辅助函数，
F＝V1+V2+…+Vn+λ(ΔQ1+ΔQ2+…+ΔQn) (9)
式中，F为辅助函数，λ为Lagrange乘子，其取极值条件为
Figure FDA0002082962730000032
联立公式(8)、(9)和(10)，利用步骤(S2)得到的各水库的防洪库容效益的函数关系，求解得到各水库防洪库容效益关系曲线微增率相等时的各水库消耗的防洪库容V1,V2,…,Vn，即为防洪库容效益最大时各水库库容的优化分配方案。
6.根据权利要求5所述的一种水库群防洪效益最大的库容分配方法，其特征在于，由公式(1)可知，水库群下游防洪控制点削减指定的洪峰流量ΔQ时，水库群消耗防洪库容最小时水库群防洪库容效益最大；
将水库的防洪库容效益函数转换表达为消耗库容对下游防洪控制点削减的洪峰量值的偏导数相等时求得各水库消耗的防洪库容V1,V2,…,Vn，此时各水库的库容微增率相等，水库群消耗防洪库容最小，即为水库群防洪库容效益最大的各水库库容的优化分配方案。
7.根据权利要求1-6任一所述的一种水库群防洪效益最大的库容分配方法，其特征在于，建立所述单个水库对下游防洪控制点的防洪优化调度模型之前，先分析并选择典型入库洪水。
CN201910478191.5A 2019-06-03 2019-06-03 一种水库群防洪效益最大的库容分配方法 Active CN110348600B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910478191.5A CN110348600B (zh) 2019-06-03 2019-06-03 一种水库群防洪效益最大的库容分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910478191.5A CN110348600B (zh) 2019-06-03 2019-06-03 一种水库群防洪效益最大的库容分配方法

Publications (2)

Publication Number Publication Date
CN110348600A CN110348600A (zh) 2019-10-18
CN110348600B true CN110348600B (zh) 2021-08-31

Family

ID=68181420

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910478191.5A Active CN110348600B (zh) 2019-06-03 2019-06-03 一种水库群防洪效益最大的库容分配方法

Country Status (1)

Country Link
CN (1) CN110348600B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112418491B (zh) * 2020-11-06 2021-08-03 黄河勘测规划设计研究院有限公司 一种水库剩余拦沙库容动态配置方法
CN112785151B (zh) * 2021-01-22 2023-11-07 黄河勘测规划设计研究院有限公司 串联水库联合防凌补偿调度方法及其专用系统
CN113626754B (zh) * 2021-07-12 2023-01-06 中国长江三峡集团有限公司 基于防洪对象安全的水库群决策实时动态调整方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102817335A (zh) * 2012-07-12 2012-12-12 贵州乌江水电开发有限责任公司 一种梯级水库群联合防洪优化调度的方法及系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102817335A (zh) * 2012-07-12 2012-12-12 贵州乌江水电开发有限责任公司 一种梯级水库群联合防洪优化调度的方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
" 基于B/S架构的实时校正洪水预报模型研究";樊青松等;《浙江水利水电学院学报》;20190415;第32-34页 *
"流域实时防洪调度关键技术研究与应用";钟平安;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20060815(第8期);第C037-11页 *
"综合利用水库调度风险分析理论与方法研究";张验科;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20121115(第11期);第C037-8页 *

Also Published As

Publication number Publication date
CN110348600A (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN110348600B (zh) 一种水库群防洪效益最大的库容分配方法
WO2021196552A1 (zh) 一种基于互馈关系解析的梯级水库风险评估方法及系统
CN111125969A (zh) 一种跨水库流域河道径流演算方法和系统
CN111160430A (zh) 一种基于人工智能算法的水资源优化配置方法
CN111126847B (zh) 耦合河道水动力过程的梯级水库短期优化调度方法和系统
CN108985585B (zh) 一种考虑预报不确定性影响的水库洪水资源利用风险对冲决策方法
CN109063901B (zh) 一种省级电网水电系统中长期发电能力分析方法
CN105069527A (zh) 一种基于数据挖掘技术的台区合理线损预测方法
CN109948869B (zh) 基于有序聚类离散化的脱硫系统so2出口浓度预测方法
CN112633578A (zh) 一种引调水工程影响下梯级水库群优化调度方法
CN111080157B (zh) 一种梯级水电站排磷量的调度方法和系统
CN104636834B (zh) 一种改进的联合概率规划模型系统优化方法
CN107565585A (zh) 储能装置调峰回报时间预测方法及其模型创建方法
CN113449993B (zh) 一种城市水源供水调度方法
CN109403425B (zh) 供水管网的优化调度系统
CN112036632A (zh) 一种梯级水库生态发电多目标中长期随机调度模型的构建方法以及优化调度方法
CN111428970A (zh) 一种大规模水电站群跨省外送能力分析模型及求解方法
CN110633849B (zh) 一种基于启发式大系统分解协调思想的水资源配置方法
CN109494810B (zh) 间歇性能源电站调度评价方法及该评价的限电分配方法
CN105225000A (zh) 一种基于模糊序优化的风功率概率模型非参数核密度估计方法
CN104636831B (zh) 一种面向多电网的水电站短期调峰特征值搜索方法
CN104504455A (zh) 一种梯级蓄能控制下水电站群长期优化调度方法
Yan et al. A multi-objective optimal control model of cascade pumping stations considering both cost and safety
CN113256028B (zh) 一种面向大范围长历时干旱的应急供水调配方法
CN107392460A (zh) 一种水库群多目标调度风险分析最佳均衡解的获取方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant