CN111815077A - 一种水库洪水调度优化方法及*** - Google Patents

一种水库洪水调度优化方法及*** Download PDF

Info

Publication number
CN111815077A
CN111815077A CN202010892101.XA CN202010892101A CN111815077A CN 111815077 A CN111815077 A CN 111815077A CN 202010892101 A CN202010892101 A CN 202010892101A CN 111815077 A CN111815077 A CN 111815077A
Authority
CN
China
Prior art keywords
reservoir
water level
flood
flow
scheduling
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.)
Granted
Application number
CN202010892101.XA
Other languages
English (en)
Other versions
CN111815077B (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202010892101.XA priority Critical patent/CN111815077B/zh
Publication of CN111815077A publication Critical patent/CN111815077A/zh
Application granted granted Critical
Publication of CN111815077B publication Critical patent/CN111815077B/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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of 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)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供一种水库洪水调度优化方法及***，其中方法包括获取水库相关数据，还包括以下步骤：设定初始时段t=1，获取时段初水位Z1；计算时段平均入库流量QIAvg；计算时段末库容V2；判断V2是否大于Vlim；判断Z1在调洪规程规定的水位级别，如果不是最高级，则按照调洪规程调度计算洪水调度；反之，则按照最大出库流量调度计算洪水调度；判断t<=N是否成立，输出洪水调度量。本发明提出的一种水库洪水调度优化方法及***，按照库水位及出库流量规则设置水库调度规程，该规则明确易实现，易于编程实现，通用性强。

Description

一种水库洪水调度优化方法及***
技术领域
本发明涉及水库流量计算的技术领域，特别是一种水库洪水调度优化方法及***。
背景技术
水库可以对洪水进行调蓄，削减洪峰流量，避免或减轻洪水灾害。在正常情况下，水库洪水调度是指按照水库调洪规程，对入库洪水进行调度，使得在满足水库安全运行的前提下，出库流量满足下游防护对象安全流量，保证下游人民群众生命财产安全。在非正常情况下，应首选保证水库的安全。
洪水调度的基本要求是：在进行洪水调度时，当库水位超过控制水位（一般为汛限水位）时，开始进行调度，并使得洪水过程结束后，库水位仍然保持在控制水位之下。
在2004年的《水电与新能源》杂志上刊登了彭天波、梁年生、袁林山等人的文章《黄龙滩水电厂洪水调度方案生成模型与计算方法》，针对黄龙滩水电厂洪水调度的具体情况，在指定的目标要求和内外环境下，运用定量、定性或两者相结合的方法，提出了常规调度、优化调度和自拟调度三种洪水调度模型，产生了相应的洪水调度方案，这些方案是电厂运行人员分析、比较、评判和选择洪水调度方式的基础。然而，在实际调度中，水库调度规程是水库调度的指导文件。但是水库调度规程与库水位、入库流量、闸门启闭顺序、梯级水库的联合调度、上下游安全防护对象的要求等有关，各个水库的调洪规程差异较大，难以统一规范设置，因此在开发计算机程序时，往往需要根据水库调度规程对程序进行特殊处理，难以通用推广应用。
发明内容
为了解决上述的技术问题，本发明提出的一种水库洪水调度优化方法及***，按照库水位及出库流量规则设置水库调度规程，该规则明确易实现，易于编程实现，通用性强。
本发明的第一目的是提供一种水库洪水调度优化方法，包括获取水库相关数据并设定初始时段t=1，还包括以下步骤：
步骤1：获取时段初水位Z1=Z(t)，其中，Z(t)为调度期内水库水位过程；
步骤2：计算时段平均入库流量QIAvg；
步骤3：计算时段末库容V2=V1+QIAvg*DeltT，其中，V1为时段初库容，DeltT为时段长；
步骤4：判断V2是否大于Vlim，如果V2大于Vlim，则执行步骤5；如果V2小于等于Vlim，则QO(t)=0，其中，Vlim为控制库容，QO(t)为出库流量；
步骤5：判断Z1在调洪规程规定的水位级别，如果不是最高级，则按照调洪规程调度计算洪水调度；反之，则按照最大出库流量调度计算洪水调度；
步骤6：判断t<=N是否成立，如果t<=N，则输出洪水调度量，所述洪水调度量包括库水位和出库流量，如果t>N，则t=t+1，执行步骤1，其中，N为时段的数量。
优选的是，所述水库相关数据包括水位库容关系曲线Z~V、水位出库流量关系曲线Z~Q、调洪规程、长度为N的入库流量数组QI()、控制水位Zlim、控制库容Vlim、第1时段的时段初水位ZB、时段长DeltT和水库的最大可能出库流量QOmax中至少一种，其中，入库流量为时段初瞬时流量。
在上述任一方案中优选的是，当t=1时，Z1=ZB，查询所述水位库容关系曲线Z~V获取所述时段初库容V1，其中，ZB为调度开始时刻的库水位。
在上述任一方案中优选的是，所述时段平均入库流量QIAvg的计算公式为QIAvg=(QI(t)+QI(t+1))/2，其中，QI(t)为t时段的时段初入库流量，QI(t+1)为t+1时段的时段初入库流量。
在上述任一方案中优选的是，所述步骤4还包括如果V2小于等于Vlim，根据V2查所述水位库容关系曲线，获取时段末水位Z(t+1)=Z2，t=t+1，重新执行步骤1。
在上述任一方案中优选的是，所述调洪规程调度的计算方法包括以下子步骤：
步骤501：根据所述时段初水位Z1查询对应的调洪规程水位级别，获取出库流量QO(t)；
步骤502：计算时段末库容V2=V1+(QIAvg-QO(t))*DeltT；
步骤503：判断V2是否小于控制库容Vlim，如果是，则执行步骤504，如果不是，则执行步骤505；
步骤504：按照时段初库容V1与控制库容Vlim之差，重新计算出库流量；
步骤505：根据V2查所述水位库容关系曲线Z~V，获取时段末水位Z2；令Z(t+1)=Z2。
在上述任一方案中优选的是，所述出库流量的计算公式为：QO(t)=(V1-Vlim)/DeltT，令时段末库容V2=Vlim。
在上述任一方案中优选的是，所述最大出库流量调度的计算方法包括以下子步骤：
步骤511：根据所述时段初水位Z1查所述水位出库流量关系曲线Z~Q，获取时段初出库流量QO1；
步骤512：利用水量平衡方程计算得到时段末库容和时段末出流的组合值.；
步骤513：计算库容变化范围，最高值Vmax=V1+QIAvg*DeltT，最低值Vmin=V1-QOmax*DeltT，其中，QOmax为 水库最大可能出库流量；
步骤514：在库容变化范围[Vmin,Vmax]内，采用0.618法计算时段末库容V2和时段末出库流量QO2，使得水量平衡方程成立；
步骤515：计算时段平均出库流量QO(t)=(QO1+QO2)/2；
步骤516：根据V2查所述水位库容关系曲线Z~V，求得时段末水位Z2，令Z(t+1)=Z2。
在上述任一方案中优选的是，所述水量平衡方程为V2-V1=(QIAvg-(QO1+QO2)/2)*DeltT，方程变形后得到V2+(QO2/2)*DeltT=V1+(QIAvg-QO1/2)*DeltT。
本发明的第二目的是提供一种水库洪水调度优化***，包括用于获取水库相关数据的数据获取模块，还包括以下模块：
初始化模块：用于设定初始时段t=1，获取时段初水位Z1=Z(t)，其中，Z(t)为调度期内库水位过程；
计算模块：用于计算时段平均入库流量QIAvg和时段末库容V2=V1+QIAvg*DeltT，其中，V1为时段初库容，DeltT为时段长；
判断模块：用于判断V2是否大于Vlim或Z1是否在调洪规程规定的水位之内；
洪水调度计算模块：用于按照调洪规程调度或最大出库流量调度计算洪水调度；
所述***按照如第一目的所述的方法对进行水库洪水调度优化。
优选的是，所述水库相关数据包括水位库容关系曲线Z~V、水位出库流量关系曲线Z~Q、调洪规程、长度为N的入库流量数组QI()、控制水位Zlim、控制水位对应库容Vlim、第1时段的时段初水位ZB、时段长DeltT和水库的最大可能出库流量QOmax中至少一种，其中，入库流量为时段初瞬时流量。
在上述任一方案中优选的是，当t=1时，Z1=ZB，查询所述水位库容关系曲线Z~V获取所述时段初库容V1，其中，ZB为调度开始时刻的库水位。
在上述任一方案中优选的是，所述时段平均入库流量QIAvg的计算公式为QIAvg=(QI(t)+QI(t+1))/2，其中，QI(t)为t时段的时段初入库流量，QI(t+1)为t+1时段的时段初入库流量。
在上述任一方案中优选的是，所述判断模块还用于如果V2小于等于Vlim，根据V2查所述水位库容关系曲线，获取时段末水位Z(t+1)=Z2，t=t+1。
在上述任一方案中优选的是，所述最大出库流量调度的计算方法包括以下子步骤：
步骤501：根据所述时段初水位Z1查询对应的调洪规程水位级别，获取出库流量QO(t)；
步骤502：计算时段末库容V2=V1+(QIAvg-QO(t))*DeltT；
步骤503：判断V2是否小于控制库容Vlim，如果是，则执行步骤504，如果不是，则执行步骤505；
步骤504：按照时段初库容V1与控制库容Vlim之差，重新计算出库流量；
步骤505：根据V2查所述水位库容关系曲线Z~V，获取时段末水位Z2；令Z(t+1)=Z2。
在上述任一方案中优选的是，所述出库流量计算公式为：QO(t)=(V1-Vlim(t))/DeltT，令时段末库容V2=Vlim。
在上述任一方案中优选的是，所述最大出库流量调度的计算方法包括以下子步骤：
步骤511：根据所述时段初水位Z1查所述水位出库流量关系曲线Z~Q，获取时段初出库流量QO1；
步骤512：利用水量平衡方程计算得到时段末库容和时段末出流的组合值；
步骤513：计算库容变化范围，最高值Vmax=V1+QIAvg*DeltT，最低值Vmin=V1-QOmax*DeltT，其中，QOmax为水库最大可能出库流量；
步骤514：在库容变化范围[Vmin,Vmax]内，采用0.618法计算时段末库容V2和时段末出库流量QO2，使得水量平衡方程成立；
步骤515：计算时段平均出库流量QO(t)=(QO1+QO2)/2；
步骤516：根据V2查所述水位库容关系曲线Z~V，求得时段末水位Z2，令Z(t+1)=Z2。
在上述任一方案中优选的是，所述水量平衡方程为V2-V1=(QIAvg-(QO1+QO2)/2)*DeltT，方程变形后得到V2+(QO2/2)*DeltT=V1+(QIAvg-QO1/2)*DeltT。
本发明提出了一种水库洪水调度优化方法及***，提供了一种通用的水库调度规程设置及计算方法，易于编程实现。
附图说明
图1为按照本发明的水库洪水调度优化方法的一优选实施例的流程图。
图2为按照本发明的水库洪水调度优化***的一优选实施例的模块图。
图3为按照本发明的水库洪水调度优化方法的水库调度规程的一实施例的计算流程图。
图4为按照本发明的水库洪水调度优化方法的最大出库流量调度的一实施例的计算流程图。
图5为按照本发明的水库洪水调度优化方法的全过程调度的一实施例的流量结果对比图。
图6为按照本发明的水库洪水调度优化方法的全过程调度的一实施例的库水位过程图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的阐述。
实施例一
如图1、2所示，一种水库洪水调度优化方法，执行步骤1010，数据获取模块2000获取水库相关数据，初始化模块2100设定初始时段t=1，当t=1时，Z1=ZB，查询所述水位库容关系曲线获取所述时段初库容V1。水库相关数据包括水位库容关系曲线Z~V、水位出库流量关系曲线Z~Q、调洪规程、长度为N的入库流量数组QI()、控制水位Zlim、控制库容Vlim、控制水位对应库容为时段末库容、第1时段的时段初水位ZB、时段长DeltT和水库的最大可能出库流量QOmax中至少一种，其中，入库流量为时段初瞬时流量， ZB为调度开始时刻的库水位，可以通过实际测量或人工设置得到。在步骤1010中包括制定调洪规程表示规则，所述调洪规程表示规则为从控制水位以上，对水库水位进行分级，各级水位对应出库流量，各级水位及出库流量依次增大。第1级表示库水位在ZL1之下时，出库流量为Q1，第2级表示库水位在ZL2之下时，出库流量为Q2，以此类推，在第m级水位ZLm之下时，出库流量为Qm，当水位超过最高级水位ZLm时，按照最大出库流量泄洪。当按照给定出库流量泄洪时，若时段末水位低于控制水位，则对给定出库流量进行调整，减小出库流量，使得时段末水位等于控制水位。
执行步骤1020，初始化模块2100获取时段初水位Z1=Z(t)，其中，Z(t)为调度期内库水位过程。
执行步骤1030，计算模块2200计算时段平均入库流量QIAvg，时段平均入库流量QIAvg的计算公式为QIAvg=(QI(t)+QI(t+1))/2，其中，QI(t)为t时段的时段初入库流量，QI(t+1)为t+1时段的时段初入库流量。
执行步骤1040，计算模块2200计算时段末库容V2=V1+QIAvg*DeltT，其中，V1为时段初库容，DeltT为时段长。
执行步骤1050，判断模块2300判断V2是否大于Vlim。如果V2小于等于Vlim，则顺序执行步骤1060和步骤1020，QO(t)=0，根据V2查所述水位库容关系曲线，获取时段末水位Z(t+1)=Z2，t=t+1，重新获取时段初水位Z1=Z(t)，其中，Vlim(t)为t时段的控制水位对应库容数组，QO(t)为出库流量。
如果V2大于Vlim，则执行步骤1070，判断模块2300判断Z1在调洪规程规定的水位级别。
如果Z1不是调洪规程规定的水位最高级别，则执行步骤1080，洪水调度计算模块2400按照调洪规程调度计算洪水调度。如图3所示，执行步骤1081，根据所述时段初水位Z1查询对应的调洪规程水位级别，获取出库流量QO(t)。执行步骤1082，计算时段末库容V2=V1+(QIAvg-QO(t))*DeltT。执行步骤1083，判断V2是否小于控制水位ZC对应的库容Vlim，如果V2大于等于控制水位ZC对应的库容Vlim，则直接执行步骤1085，根据V2查所述水位库容关系曲线，获取时段末水位Z2；令Z(t+1)=Z2。如果V2小于控制水位ZC对应的库容Vlim，则顺序执行步骤1084，按照时段初库容V1与控制库容Vlim之差，重新计算出库流量，QO(t)=(V1-Vlim)/DeltT，令时段末库容V2=Vlim。执行步骤1085，根据V2查所述水位库容关系曲线，获取时段末水位Z2；令Z(t+1)=Z2。
如果Z1是调洪规程规定的水位最高级别，如图4所示，则执行步骤1090，洪水调度计算模块2400按照最大出库流量调度计算洪水调度。执行步骤1091，根据所述时段初水位Z1查所述水位出库流量关系曲线Z~Q，获取时段初出库流量QO1。执行步骤1092，利用水量平衡方程计算得到时段末库容和时段末出流的组合值，水量平衡方程为V2-V1=(QIAvg-(QO1+QO2)/2)*DeltT，方程变形后得到V2+(QO2/2)*DeltT=V1+(QIAvg-QO1/2)*DeltT。执行步骤1093，计算库容变化范围，最高值Vmax=V1+QIAvg*DeltT，最低值Vmin=V1-QOmax*DeltT，其中，QOmax为可能最大出库流量。执行步骤1094，在库容变化范围[Vmin,Vmax]内，采用0.618法计算时段末库容V2和时段末出库流量QO2，使得水量平衡方程成立。执行步骤1095，计算时段平均出库流量QO(t)=(QO1+QO2)/2。执行步骤1096，根据V2查所述水位库容关系曲线Z~V，求得时段末水位Z2，令Z(t+1)=Z2。
执行完步骤1080或步骤1090后，执行步骤1100，判断t<=N是否成立,其中，N为时段的数量。如果t>N，则顺序执行步骤1110和步骤1020，t=t+1，重新获取时段初水位Z1=Z(t)。
如果t<=N，则执行步骤1120，洪水调度计算模块2400输出洪水调度量。
实施例二
本发明提出一种通用的调洪规程表示规则，并提出该调洪规程下的洪水调度方法。
1、调洪规程表示规则
从控制水位Zlim以上，对水库水位进行分级，各级水位对应出库流量，各级水位及出库流量依次增大，如表1所示：
表1 水库调洪规程
Figure 672345DEST_PATH_IMAGE001
表1中，第1级表示库水位在ZL1之下时，出库流量为Q1，第2级表示库水位在ZL2之下时，出库流量为Q2，以此类推，在第m级水位ZLm之下时，出库流量为Qm，当水位超过最高级水位ZLm时，按照最大出库流量泄洪。
当按照给定出库流量泄洪时，若时段末水位低于控制水位，则对给定出库流量进行调整，减小出库流量，使得时段末水位等于控制水位。
2、水库调度计算方法
已知条件：水位库容关系曲线（Z~V曲线）、水位出库流量关系曲线（Z~Q曲线）、调洪规程（各级水位及出库流量）、长度为N的入库流量数组QI()、控制水位、控制库容Vlim，其中入库流量为时段初瞬时流量、第1时段的时段初水位ZB、DeltT-时段长（分钟），水库的最大可能出库流量QOmax。
求解内容：Z()-库水位数组（时段初水位），QO()-出库流量数组（时段平均流量），
2.1基本原理
对于某一个调度时段，如果时段入库流量全部入库后，时段末水位（库容）仍低于控制水位（库容），则出库流量为0；如果时段末水位（库容）超过控制水位（库容），则需要开闸放水，如果时段初水位在调洪规程规定的分级水位以下，则按照调洪规程的出库流量调度，如果时段初水位高于调洪规程规定的最高分级水位，则按照最大出库流量调度。
2.2计算步骤
从时段t=1开始计算到t=N，计算步骤如下：
1、获取时段初水位Z1=Z(t)，当t=1时，Z1=ZB，查询水位库容曲线Z~V获取时段初库容V1；
2、计算时段平均入库流量QIAvg=(QI(t)+QI(t+1))/2，当t=N时，令QIAvg=QI(t)；
3、计算时段末库容V2=V1+QIAvg*DeltT；
4、判断V2是否大于Vlim，如果是，进行第5步；如果不是，则QO(t)=0，，根据V2查水位库容曲线Z~V，获取时段末水位Z(t+1)=Z2，t=t+1，返回第1步。
5、判断Z1是否在调洪规程规定的水位之内，如果是，则按照规程调度；如果不是，则按照最大出库流量调度。
6、判断是否是最后一个时段，如果是，则结束，如果不是，则t=t+1，返回1。
2.3按调洪规程计算步骤
按照调洪规程调度计算t时段的洪水调度步骤如下：
1、根据时段初水位Z1查询对应的调洪规程水位级别，获取出库流量QO(t)；
2、计算时段末库容V2=V1+(QIAvg-QO(t))*DeltT；
3、判断V2是否小于控制库容Vlim，如果是，则进入第4步，如果不是，则进入第5步；
4、按照时段初库容V1与控制库容Vlim之差，重新计算出库流量，QO(t)=(V1-Vlim)/DeltT，令时段末库容V2=Vlim;
5、根据V2查水位库容曲线Z ~V，获取时段末水位Z2；令Z(t+1)=Z2。
2.4按照最大出库流量调度计算步骤
按照最大出库流量调度计算t时段的洪水调度步骤如下：
1、根据时段初水位Z1，查水位出库流量曲线Z~Q，获取时段初出库流量QO1；
2、水量平衡方程V2-V1=(QIAvg-(QO1+QO2)/2)*DeltT，其中QO2为时段末出库流量，方程变形为V2+(QO2/2)*DeltT=V1+(QIAvg-QO1/2)*DeltT，则方程的等号右半部分可以由已知条件计算；
3、计算库容变化范围，最高值Vmax=V1+QIAvg*DeltT，最低值Vmin=V1-QOmax*DeltT；
4、在库容变化范围[Vmin,Vmax]内，采用0.618法计算时段末库容V2和时段末出库流量QO2，使得步骤2中水量平衡方程成立；
5、计算时段平均出库流量QO(t)=(QO1+QO2)/2；
6、根据V2查水位库容关系曲线Z~V，求得时段末水位Z2，令Z(t+1)=Z2。
实施例三
某水库的水位库容关系曲线Z~V，水位出库流量关系曲线Z~Q如表2所示，调洪规程如表3所示。
表2 某水库水位-库容-出库流量关系曲线表
Figure 154273DEST_PATH_IMAGE002
表3 水库调洪规程
Figure 475533DEST_PATH_IMAGE003
该水库的调洪控制水位为27.6m，即当库水位超过27.6m时，进行洪水调度。该水库调洪规程表示，当27.6<Z<=28时，出库流量为100m3/s；当28<Z<=28.3时，出库流量为200m3/s；当Z>28.3时，出库流量为最大出库流量。
某场洪水的计算时段数为72小时，计算时段长为1小时，时段初水位为27.6m，最大可能下泄流量为1000m3/s，计算结果如表4所示，入库流量和出库流量过程如图5所示，库水位过程如图6所示。
如表4所示，第1-33时段的时段初水位低于调洪规程的第1级水位28m，因此按照第1级出库流量进行调度。第1-22时段，各时段出库流量等于入库流量，库水位保持在控制水位27.6m，这是因为按照调洪规程的第1级水位调度时，如果出库流量为100m3/s，则时段末水位会低于27.6m，因此需要重新计算出库流量，使得库水位保持在27.6m；第23-25时段，按照第1级出库流量100m3/s进行调度；第26-29时段，按照第1级出库流量调度，需要调整出库流量，使得时段末水位等于27.6m；第30-33时段，按照第1级出库流量100m3/s调度。
第34时段的时段初水位28.22m高于第1级水位28m，低于第2级水位28.3m，因此按照第2级出库流量200m3/s调度；
第35-43时段，各时段初水位高于第2级水位28.3m，按照最大出库流量调度；
第44-46时段，各时段初水位28.22m高于第1级水位28m，低于第2级水位28.3m，按照第2级出库流量200m3/s调度；
第47-52时段，各时段初水位低于第1级水位28m，按照第1级出库流量100m3/s调度；
第53-72时段，时段初水位低于第1级水位28m，按照第1级出库流量调度时，时段末水位低于27.6m，因此对出库流量进行调整，使得出库流量等于入库流量，时段末水位等于27.6m。m3/s。
表4 洪水调度计算成果表
Figure 511753DEST_PATH_IMAGE004
为了更好地理解本发明，以上结合本发明的具体实施例做了详细描述，但并非是对本发明的限制。凡是依据本发明的技术实质对以上实施例所做的任何简单修改，均仍属于本发明技术方案的范围。本说明书中每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似的部分相互参见即可。对于***实施例而言，由于其与方法实施例基本对应，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

Claims (10)

1.一种水库洪水调度优化方法，包括获取水库相关数据并设定初始时段t=1，其特征在于，还包括以下步骤：
步骤1：获取时段初水位Z1=Z(t)，其中，Z(t)为调度期内库水位过程；
步骤2：计算时段平均入库流量QIAvg；
步骤3：计算时段末库容V2=V1+QIAvg*DeltT，其中，V1为时段初库容，DeltT为时段长；
步骤4：判断V2是否大于Vlim，如果V2大于Vlim，则执行步骤5；如果V2小于等于Vlim，则QO(t)=0，其中，Vlim为控制水位对应库容，QO(t)为出库流量；
步骤5：判断Z1在调洪规程规定的水位级别，如果不是最高级，则按照调洪规程调度计算洪水调度；反之，则按照最大出库流量调度计算洪水调度；
步骤6：判断t<=N是否成立，如果t<=N，则输出洪水调度量，所述洪水调度量包括库水位、出库流量，如果t>N，则t=t+1，执行步骤1，其中，N为时段的数量。
2.如权利要求1所述的水库洪水调度优化方法，其特征在于，所述水库相关数据包括水位库容关系曲线Z~V、水位出库流量关系曲线Z~Q、调洪规程、长度为N的入库流量数组QI()、控制水位Zlim、控制库容Vlim、第1时段的时段初水位ZB、时段长DeltT和水库的最大可能出库流量QOmax中至少一种，其中，入库流量为时段初瞬时流量。
3.如权利要求2所述的水库洪水调度优化方法，其特征在于，当t=1时，Z1=ZB，查询所述水位库容关系曲线Z~V获取所述时段初库容V1，其中，ZB为调度开始时刻的库水位。
4.如权利要求2所述的水库洪水调度优化方法，其特征在于，所述时段平均入库流量QIAvg的计算公式为QIAvg=(QI(t)+QI(t+1))/2，其中，QI(t)为t时段的时段初入库流量，QI(t+1)为t+1时段的时段初入库流量。
5.如权利要求4所述的水库洪水调度优化方法，其特征在于，所述步骤4还包括如果V2小于等于Vlim，根据V2查所述水位库容关系曲线，获取时段末水位Z(t+1)=Z2，t=t+1，重新执行步骤1。
6.如权利要求5所述的水库洪水调度优化方法，其特征在于，所述调洪规程调度的计算方法包括以下子步骤：
步骤501：根据所述时段初水位Z1查询对应的调洪规程水位级别，获取出库流量QO(t)；
步骤502：计算时段末库容V2=V1+(QIAvg-QO(t))*DeltT；
步骤503：判断V2是否小于控制库容Vlim，如果是，则执行步骤504，如果不是，则执行步骤505；
步骤504：按照时段初库容V1与Vlim之差，重新计算出库流量；
步骤505：根据V2查所述水位库容关系曲线，获取时段末水位Z2；令Z(t+1)=Z2。
7.如权利要求6所述的水库洪水调度优化方法，其特征在于，所述出库流量的计算公式为：QO(t)=(V1-Vlim)/DeltT，令时段末库容V2=Vlim。
8.如权利要求5所述的水库洪水调度优化方法，其特征在于，所述最大出库流量调度的计算方法包括以下子步骤：
步骤511：根据所述时段初水位Z1查所述水位出库流量关系曲线Z~Q，获取时段初出库流量QO1；
步骤512：利用水量平衡方程计算得到时段末库容和时段末出流的组合值；
步骤513：计算库容变化范围，最高值Vmax=V1+QIAvg*DeltT，最低值Vmin=V1-QOmax*DeltT，其中，QOmax为水库最大可能出库流量；
步骤514：在库容变化范围[Vmin,Vmax]内，采用0.618法计算时段末库容V2和时段末出库流量QO2，使得水量平衡方程成立；
步骤515：计算时段平均出库流量QO(t)=(QO1+QO2)/2；
步骤516：根据V2查所述水位库容关系曲线Z~V，求得时段末水位Z2，令Z(t+1)=Z2。
9.如权利要求8所述的水库洪水调度优化方法，其特征在于，所述水量平衡方程为V2-V1=(QIAvg-(QO1+QO2)/2)*DeltT，方程变形后得到V2+(QO2/2)*DeltT=V1+(QIAvg-QO1/2)*DeltT。
10.一种水库洪水调度优化***，包括用于获取水库相关数据的数据获取模块，其特征在于，还包括以下模块：
初始化模块：用于设定初始时段t=1，获取时段初水位Z1=Z(t)，其中，Z(t)为调度期内库水位过程；
计算模块：用于计算时段平均入库流量QIAvg和时段末库容V2=V1+QIAvg*DeltT，其中，V1为时段初库容，DeltT为时段长；
判断模块：用于判断V2是否大于Vlim或Z1是否在调洪规程规定的水位之内；
洪水调度计算模块：用于按照调洪规程调度或最大出库流量调度计算洪水调度；
所述***按照如权利要求1所述的方法对进行水库洪水调度优化。
CN202010892101.XA 2020-08-31 2020-08-31 一种水库洪水调度优化方法及*** Active CN111815077B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010892101.XA CN111815077B (zh) 2020-08-31 2020-08-31 一种水库洪水调度优化方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010892101.XA CN111815077B (zh) 2020-08-31 2020-08-31 一种水库洪水调度优化方法及***

Publications (2)

Publication Number Publication Date
CN111815077A true CN111815077A (zh) 2020-10-23
CN111815077B CN111815077B (zh) 2020-12-11

Family

ID=72859742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010892101.XA Active CN111815077B (zh) 2020-08-31 2020-08-31 一种水库洪水调度优化方法及***

Country Status (1)

Country Link
CN (1) CN111815077B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113177189A (zh) * 2021-05-20 2021-07-27 中国水利水电科学研究院 一种梯级水库分级分期旱限库容的计算方法
CN114118554A (zh) * 2021-11-18 2022-03-01 武汉大学 基于台风路径的水库汛期运行水位动态控制方法及***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102298731A (zh) * 2010-06-25 2011-12-28 华东电网有限公司 考虑顶潮供水综合要求的梯级水库短期发电优化调度方法
CN102776870A (zh) * 2012-07-12 2012-11-14 贵州乌江水电开发有限责任公司 基于水位流量的梯级水电站的无预报洪水调度方法
KR101670918B1 (ko) * 2015-10-29 2016-10-31 대한민국 지상과 지하의 실시간 침수 해석방법
CN107180313A (zh) * 2017-06-05 2017-09-19 中国水利水电科学研究院 基于风险对冲规则的大型水库汛末蓄水方案的编制方法
CN110232479A (zh) * 2019-06-13 2019-09-13 福州市规划设计研究院 一种城市水库防洪补偿优化调度方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102298731A (zh) * 2010-06-25 2011-12-28 华东电网有限公司 考虑顶潮供水综合要求的梯级水库短期发电优化调度方法
CN102776870A (zh) * 2012-07-12 2012-11-14 贵州乌江水电开发有限责任公司 基于水位流量的梯级水电站的无预报洪水调度方法
KR101670918B1 (ko) * 2015-10-29 2016-10-31 대한민국 지상과 지하의 실시간 침수 해석방법
CN107180313A (zh) * 2017-06-05 2017-09-19 中国水利水电科学研究院 基于风险对冲规则的大型水库汛末蓄水方案的编制方法
CN110232479A (zh) * 2019-06-13 2019-09-13 福州市规划设计研究院 一种城市水库防洪补偿优化调度方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113177189A (zh) * 2021-05-20 2021-07-27 中国水利水电科学研究院 一种梯级水库分级分期旱限库容的计算方法
CN113177189B (zh) * 2021-05-20 2022-02-01 中国水利水电科学研究院 一种梯级水库分级分期旱限库容的计算方法
CN114118554A (zh) * 2021-11-18 2022-03-01 武汉大学 基于台风路径的水库汛期运行水位动态控制方法及***
CN114118554B (zh) * 2021-11-18 2024-05-14 武汉大学 基于台风路径的水库汛期运行水位动态控制方法及***

Also Published As

Publication number Publication date
CN111815077B (zh) 2020-12-11

Similar Documents

Publication Publication Date Title
CN111815077B (zh) 一种水库洪水调度优化方法及***
CN107818385B (zh) 一种梯级水电站群实时运行趋势预测方法
CN105760213B (zh) 云环境下虚拟机资源利用率的预警***及方法
CN112184070B (zh) 一种协同生态流量需求的梯级水电站多目标优化调度方法及***
CN110599363A (zh) 一种计及梯级水电站优化调度的电力***可靠性评估方法
CN105243438A (zh) 一种考虑径流不确定性的多年调节水库优化调度方法
CN110739726B (zh) 一种考虑海上风电接入的多类型电源容量长期规划方法
CN113530793B (zh) 一种空压气站用智能调节***
CN105335561A (zh) 一种基于指标排序的梯级水电站群超短期调度方法
WO2016000133A1 (en) Method and system of efficiency evaluation of rcaes system
CN110912200A (zh) 一种梯级水电站优化调度***、方法及安全电网***
CN115619180A (zh) 一种确定多能互补***发电方案的方法、装置及电子设备
CN113592281A (zh) 电站的性能分析方法、装置、设备及存储介质
CN112950096A (zh) 厂网河一体化分类智能调度方法
CN104063810A (zh) 一种基于大数据的用电因素影响模型
CN110729721A (zh) 电力***全局备用容量计算方法
CN109409569A (zh) 一种考虑直流输电约束的水库群中长期调度的离散梯度逐步优化算法
CN111009927A (zh) 一种风电配合压缩空气储能容量优化方法及***
CN113506185A (zh) 梯级水电站的发电优化调度方法、装置和计算机设备
CN117910711B (zh) 均衡风险与效益的汛期水位智能动态控制模型的构建方法
CN118352989A (zh) 一种多场景下的电力电量综合平衡裕度计算方法及***
CN117584790B (zh) 一种无增容充电桩控制***
CN118278770A (zh) 一种基于机组数据挖掘的机组agc性能排序与关键机组甄别方法
CN108241912B (zh) 抽油机井盈亏平衡周期确定方法及装置
CN115099006A (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