CN110852525A - 基于预报误差越级设防的水库防洪限制水位动态调节方法 - Google Patents

基于预报误差越级设防的水库防洪限制水位动态调节方法 Download PDF

Info

Publication number
CN110852525A
CN110852525A CN201911136441.3A CN201911136441A CN110852525A CN 110852525 A CN110852525 A CN 110852525A CN 201911136441 A CN201911136441 A CN 201911136441A CN 110852525 A CN110852525 A CN 110852525A
Authority
CN
China
Prior art keywords
rain
forecast
rainfall
reservoir
grade
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
CN201911136441.3A
Other languages
English (en)
Other versions
CN110852525B (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.)
Changan University
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN201911136441.3A priority Critical patent/CN110852525B/zh
Publication of CN110852525A publication Critical patent/CN110852525A/zh
Application granted granted Critical
Publication of CN110852525B publication Critical patent/CN110852525B/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)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Control Of Non-Electrical Variables (AREA)
  • Alarm Systems (AREA)

Abstract

本发明公开了一种基于预报误差越级设防的水库防洪限制水位动态调节方法，包括步骤：一、降水预报等级划分；二、确定越级控制雨量阈值；三、获取未来24小时降水预报等级信息并判断未来24小时降水预报等级是否为小雨及小雨以下降水预报等级；四、小雨及小雨以下降水预报等级的水库防洪限制水位动态调节；五、中雨及中雨以上降水预报等级的水库防洪限制水位动态调节。本发明按照在流域上发生平均面雨量相当于预报越级降水量来对水库进行设防，再由水库防洪库容和安全控泄流量来反推对应的安全限制水位，基于降水预报推求防洪限制水位，考虑预报精度的同时还将可能的降水量取值范围进行越级拓宽，增加水库对洪水资源的利用，降低了水库的调洪风险。

Description

基于预报误差越级设防的水库防洪限制水位动态调节方法
技术领域
本发明属于防洪限制水位调节技术领域，具体涉及一种基于预报误差越级设防的水库防洪限制水位动态调节方法。
背景技术
受气候变化及人类活动影响，水资源时空分配更加趋于不均，加之我国社会经济的高速发展，人们对水资源需求的缺口愈加增大。通过大江大河上的水库对地表径流进行调节，是满足社会用水需求的重要措施。水库洪水调度是水库对河川径流过程调节的一种重要类型。水库洪水调度长期以来实行基于不考虑预报信息的静态调度模式，即，基于洪水频率分析及水文统计，设定固定的汛限水位和确定的调度规则来调节洪水。这种静态调度模式造成洪水资源难以充分利用，具体表现为：在汛期，为确保防洪安全不敢蓄水，要严格控制水库水位不超过汛限水位，留出较大的防洪库容准备调蓄可能会发生的设计洪水和校核洪水，较低的库水位影响供水、发电效益；而汛期过后河道来水量很少，水库在汛限水位以上的库容又比较大，常常造成库水位难以达到正常蓄水位，严重影响到水库的兴利和洪水资源利用。
对于水库防洪调度而言，防洪限制水位是一个控制性的关键水位，其关系到水库兴利和防洪安全两个方面的重要参数。防洪限制水位低则发电水头低，水库蓄水不足、降低兴利效益，但防洪相对偏于安全；限制水位高则恰恰相反，增加了蓄水量和水头，减小了调洪库容，防洪风险上升。长期以来，水库防洪限制水位的定制仍是基于泄流能力和由历史洪水频率分析得到入库设计洪水过程经调洪演算确定得到，频率分析时采用的洪水要素均为年最大值系列。由历史洪水数据系列确定得到的汛限水位显然是要求水库在汛期时刻以防备发生设计或校核标准(发生概率很小)的洪水为目标，其防范级别显然较高，影响蓄水发电和洪水资源利用。在气候变化和人类活动等多重作用影响下，河川水文系列变异明显，以历史数据作为设计依据指导水库未来调度弊端更加显现。
随着遥测、遥感与卫星技术的飞速发展，包括基于数值分析、卫星云图、测雨雷达的天气形势与降水预报信息，已在气象预报上发挥了重要作用，降水预报的精度提高，预见期延长。我国天气预报每日发布全国和各地区24小时、48小时及72小时降水预报等级分为无雨、小雨、中雨、大雨、暴雨和大暴雨，由于预报的不确定性，实际天气的降水预报等级在给定的降水预报等级条件下，有可能会发生相邻等级的降水事件，但越级的降水事件概率极低，据有关数据分析表明，我国未来24小时预报的准确率已经超过80％，定性分级预报准确率已经达到90％以上，这些信息对于水库防洪调度、充分利用洪水资源意义重大。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于预报误差越级设防的水库防洪限制水位动态调节方法，按照在流域上发生平均面雨量相当于预报越级降水量来对水库进行设防，再由水库防洪库容和安全控泄流量来反推对应的安全限制水位，基于未来短期降水预报推求防洪限制水位考虑预报精度的同时，还将可能的降水量取值范围进行越级拓宽，一方面增加了水库对洪水资源的利用，另一方面在更大程度上降低了水库的调洪风险，便于推广使用。
为解决上述技术问题，本发明采用的技术方案是：基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于，该方法包括以下步骤：
步骤一、降水预报等级划分：对天气预报播报的降雨信息，按照流域或区域面雨量大小进行降水预报等级划分，所述降水预报等级由低至高依次包括无雨、小雨、中雨、大雨、暴雨和大暴雨，所述无雨按预报降水概率由低至高依次包括晴天、多云和阴天，所述无雨的实际降水区间为[0，0.1mm)，所述小雨的实际降水区间为[0.1mm，6mm)，所述中雨的实际降水区间为[6mm，15mm)，所述大雨的实际降水区间为[15mm，30mm)，所述暴雨的实际降水区间为[30mm，60mm)，所述大暴雨的实际降水区间为[60mm，∞)；
步骤二、确定越级控制雨量阈值：流域内未来24小时降水预报为晴天间多云时，其对应的越级为小雨及小雨以上降水预报等级，其对应的越级控制雨量阈值取3mm；
流域内未来24小时降水预报为多云间阴天时，其对应的越级为中雨及中雨以上降水预报等级，其对应的越级控制雨量阈值取10mm；
流域内未来24小时降水预报为小雨时，其对应的越级为大雨及大雨以上降水预报等级，其对应的越级控制雨量阈值取15mm；
流域内未来24小时降水预报为中雨时，其对应的越级为无雨或暴雨及暴雨以上降水预报等级，其对应的越级控制雨量上限阈值取30mm，其对应的越级控制雨量下限阈值取3mm；
流域内未来24小时降水预报为大雨及大雨以上降水预报等级时，其对应的越级为中雨及中雨以下降水预报等级，其对应的越级控制雨量下限阈值取10mm；
步骤三、获取未来24小时降水预报等级信息并判断未来24小时降水预报等级是否为小雨及小雨以下降水预报等级：根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为小雨及小雨以下降水预报等级时，执行步骤四；否则，执行步骤五；
步骤四、小雨及小雨以下降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤401、根据公式
Figure BDA0002279712650000031
计算小雨及小雨以下降水预报等级下未来24小时水量其中，W'为水库24小时汇流区河道槽蓄水量，
Figure BDA0002279712650000033
为未来24小时内对应降水预报等级的越级控制雨量阈值，单位为mm，α为流域内径流系数，F为水库24小时汇流区面积，单位为km2；
步骤402、根据公式
Figure BDA0002279712650000041
计算小雨及小雨以下降水预报等级下第一限制水位值
Figure BDA0002279712650000042
其中，f[·]为水库水位库容关系函数，Zh为防洪高水位，V(Zh)为防洪高水位Zh下对应的防洪库容，Qu为水库生态基流常流量；
步骤403、根据公式
Figure BDA0002279712650000043
计算小雨及小雨以下降水预报等级下未来24小时越级控制雨量对应的平均来流
Figure BDA0002279712650000044
步骤404、根据公式
Figure BDA0002279712650000045
计算小雨及小雨以下降水预报等级下未来24小时内平均可附加下泄流量其中，Qs为水库安全允许泄流量；
步骤405、根据公式
Figure BDA0002279712650000047
计算小雨及小雨以下降水预报等级下未来24小时内可预泄水库水量Wout；
步骤406、根据公式计算小雨及小雨以下降水预报等级下第二限制水位值其中，Zl为水库的原防洪限制水位，V(Zl)为水库的原防洪限制水位Zl对应的库容；
步骤407、比较小雨及小雨以下降水预报等级下第一限制水位值
Figure BDA00022797126500000410
和小雨及小雨以下降水预报等级下第二限制水位值
Figure BDA00022797126500000411
的数值大小，当
Figure BDA00022797126500000412
时，小雨及小雨以下降水预报等级的水库防洪限制水位
Figure BDA00022797126500000413
当
Figure BDA00022797126500000414
时，小雨及小雨以下降水预报等级的水库防洪限制水位Zx取
Figure BDA00022797126500000415
和
Figure BDA00022797126500000416
中的较小值；
步骤五、中雨及中雨以上降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤501、根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，根据公式
Figure BDA00022797126500000417
计算当前预报量级下未来24小时越级水量
Figure BDA00022797126500000418
其中，
Figure BDA00022797126500000419
为未来24小时内对应降水预报等级的越级控制雨量上限阈值，单位为mm；根据公式
Figure BDA00022797126500000420
计算当前预报量级下第一限制水位值
Figure BDA0002279712650000051
当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，且当前预报量级下第一限制水位值
Figure BDA0002279712650000052
高于防洪限制水位Zl时，取中雨及中雨以上降水预报等级下第一限制水位值
Figure BDA0002279712650000053
步骤502、根据公式
Figure BDA0002279712650000054
计算中雨及中雨以上降水预报等级下未来24小时水量
Figure BDA0002279712650000055
其中，
Figure BDA0002279712650000056
为未来24小时内对应降水预报等级的越级控制雨量下限阈值，单位为mm；
步骤503、根据公式
Figure BDA0002279712650000057
计算中雨及中雨以上降水预报等级下未来24小时越级控制雨量对应的平均来流
Figure BDA0002279712650000058
步骤504、根据公式计算中雨及中雨以上降水预报等级下未来24小时内超出兴利用水的平均流量Qa；
步骤505、根据公式Wa＝Qa×86400，计算中雨及中雨以上降水预报等级下未来24小时内超出流量部分对应的有效水量Wa；
步骤506、根据公式
Figure BDA00022797126500000510
计算中雨及中雨以上降水预报等级下第二限制水位值
Figure BDA00022797126500000511
步骤507、天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，比较
Figure BDA00022797126500000512
和
Figure BDA00022797126500000513
的数值大小，当
Figure BDA00022797126500000514
与
Figure BDA00022797126500000515
相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure BDA00022797126500000516
当与
Figure BDA00022797126500000518
不相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure BDA00022797126500000519
取
Figure BDA00022797126500000520
与
Figure BDA00022797126500000521
中的较大值；再比较中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure BDA00022797126500000522
和原防洪限制水位Zl的数值大小，当
Figure BDA00022797126500000523
时，中雨及中雨以上降水预报等级的水库防洪限制水位
Figure BDA00022797126500000524
当
Figure BDA00022797126500000525
时，中雨及中雨以上降水预报等级的水库防洪限制水位Zzd＝Zl。
上述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述水库为暴雨洪水入流类型且承担防洪任务的水库，所述暴雨洪水入流类型水库包括供水水库、灌溉水库、发电水库及多目标综合利用水库；
当水库为供水水库或灌溉水库时，水库生态基流常流量Qu为满足河道下游生态及航运要求的下泄流量；
当水库为发电水库或多目标综合利用水库时，水库生态基流常流量Qu取发电下泄等常流量和下游生态流量两者中的大值。
上述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述水库为单一水库。
上述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述流域为扇形流域或汇流面积不大于50000km2的羽形流域。
上述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述流域内径流系数α取值范围为0.6～0.8。
本发明与现有技术相比具有以下优点：
1、本发明对天气预报播报的降雨信息，按照流域或区域面雨量大小进行降水预报等级划分，考虑天气预报的不确定性，通过基于未来短期降水预报推求防洪限制水位考虑预报精度的同时，还将可能的降水量取值范围进行越级拓宽，进而确定各个降水预报等级的越级控制雨量阈值，为水库防洪限制水位动态调节提供参考的基础，便于推广使用。
2、本发明针对不同降水预报等级制定不同的水库防洪限制水位动态调节方式，根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为小雨及小雨以下降水预报等级时，以洪水资源利用为主要目标，按预设的越级控制雨量阈值防控洪水风险；当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，以控制和降低洪水风险为主要目标，按预设的越级控制雨量阈值上限降低调洪风险，按预设的越级控制雨量阈值下限减少兴利损失，在更大程度上降低了水库的调洪风险，使用效果好。
3、本发明方法步骤简单，针对不同降水预报等级，均计算两个限制水位值，其中，不同降水预报等级下的第一限制水位值的制定主要是站在洪量控制的角度上，计算过程均以防洪高水位或设计洪水位为基础进行；不同降水预报等级下的第二限制水位值的拟定则主要是基于洪水资源利用或减少兴利库容损失的要求，计算过程以常规防洪高水位为基础进行，便于推广使用。
综上所述，本发明按照在流域上发生平均面雨量相当于预报越级降水量来对水库进行设防，再由水库防洪库容和安全控泄流量来反推对应的安全限制水位，基于未来短期降水预报推求防洪限制水位考虑预报精度的同时，还将可能的降水量取值范围进行越级拓宽，一方面增加了水库对洪水资源的利用，另一方面在更大程度上降低了水库的调洪风险，便于推广使用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明方法的方法流程框图。
具体实施方式
如图1所示，本发明的基于预报误差越级设防的水库防洪限制水位动态调节方法，包括以下步骤：
步骤一、降水预报等级划分：对天气预报播报的降雨信息，按照流域或区域面雨量大小进行降水预报等级划分，所述降水预报等级由低至高依次包括无雨、小雨、中雨、大雨、暴雨和大暴雨，所述无雨按预报降水概率由低至高依次包括晴天、多云和阴天，所述无雨的实际降水区间为[0，0.1mm)，所述小雨的实际降水区间为[0.1mm，6mm)，所述中雨的实际降水区间为[6mm，15mm)，所述大雨的实际降水区间为[15mm，30mm)，所述暴雨的实际降水区间为[30mm，60mm)，所述大暴雨的实际降水区间为[60mm，∞)；
本实施例中，所述水库为单一水库，仅基于单一水库考虑其汛期限制水位的动态规划，不针对梯级和水库群；所述流域为扇形流域或汇流面积不大于50000km2的羽形流域；且参照水文情报预报规范，采用的天气预报产品的合格率大于85％，预报精度不低于甲级水平。
步骤二、确定越级控制雨量阈值：流域内未来24小时降水预报为晴天间多云时，其对应的越级为小雨及小雨以上降水预报等级，其对应的越级控制雨量阈值取3mm，取小雨中位数；
流域内未来24小时降水预报为多云间阴天时，其对应的越级为中雨及中雨以上降水预报等级，其对应的越级控制雨量阈值取10mm，取中雨中位数；
流域内未来24小时降水预报为小雨时，其对应的越级为大雨及大雨以上降水预报等级，其对应的越级控制雨量阈值取15mm，取大雨下限；
流域内未来24小时降水预报为中雨时，其对应的越级为无雨或暴雨及暴雨以上降水预报等级，其对应的越级控制雨量上限阈值取30mm，取暴雨下限，其对应的越级控制雨量下限阈值取3mm，取小雨中位数；
流域内未来24小时降水预报为大雨及大雨以上降水预报等级时，其对应的越级为中雨及中雨以下降水预报等级，其对应的越级控制雨量下限阈值取10mm，取中雨中位数；
需要说明的是，对天气预报播报的降雨信息，按照流域或区域面雨量大小进行降水预报等级划分，考虑天气预报的不确定性，通过基于未来短期降水预报推求防洪限制水位考虑预报精度的同时，还将可能的降水量取值范围进行越级拓宽，进而确定各个降水预报等级的越级控制雨量阈值，为水库防洪限制水位动态调节提供参考的基础。
步骤三、获取未来24小时降水预报等级信息并判断未来24小时降水预报等级是否为小雨及小雨以下降水预报等级：根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为小雨及小雨以下降水预报等级时，执行步骤四；否则，执行步骤五；
实际使用时，考虑未来72小时降水预报等级信息，根据未来24小时内的降水预报等级来进行推导计算，即24小时降水预报是否为小雨及小雨以下等级，确定并修正水库防洪限制水位，将未来48小时和72小时的降水预报等级信息作为参考，辅助作业人员做好水位调节准备工作。
步骤四、小雨及小雨以下降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤401、根据公式计算小雨及小雨以下降水预报等级下未来24小时水量
Figure BDA0002279712650000092
其中，W'为水库24小时汇流区河道槽蓄水量，
Figure BDA0002279712650000093
为未来24小时内对应降水预报等级的越级控制雨量阈值，单位为mm，α为流域内径流系数，F为水库24小时汇流区面积，单位为km2；
本实施例中，所述流域内径流系数α取值范围为0.6～0.8。
实际实施时，流域内径流系数α取值参考前期降雨情况而定。
步骤402、根据公式
Figure BDA0002279712650000094
计算小雨及小雨以下降水预报等级下第一限制水位值
Figure BDA0002279712650000095
其中，f[·]为水库水位库容关系函数，Zh为防洪高水位，V(Zh)为防洪高水位Zh下对应的防洪库容，Qu为水库生态基流常流量；
本实施例中，所述水库为暴雨洪水入流类型且承担防洪任务的水库，所述暴雨洪水入流类型水库包括供水水库、灌溉水库、发电水库及多目标综合利用水库；
当水库为供水水库或灌溉水库时，水库生态基流常流量Qu为满足河道下游生态及航运要求的下泄流量；
当水库为发电水库或多目标综合利用水库时，水库生态基流常流量Qu取发电下泄等常流量和下游生态流量两者中的大值。
需要说明的是，水库防洪限制水位的制定考虑未来24小时降水预报，将预报等级越级拓宽，按照越级后相应的雨量折算入库水量，再从防洪库容中扣除对应水量，作为防洪安全储备，并以扣除后的库容基于水位库容关系求得限制水位。
步骤403、根据公式
Figure BDA0002279712650000096
计算小雨及小雨以下降水预报等级下未来24小时越级控制雨量对应的平均来流
Figure BDA0002279712650000097
步骤404、根据公式
Figure BDA0002279712650000101
计算小雨及小雨以下降水预报等级下未来24小时内平均可附加下泄流量
Figure BDA0002279712650000102
其中，Qs为水库安全允许泄流量；
步骤405、根据公式
Figure BDA0002279712650000103
计算小雨及小雨以下降水预报等级下未来24小时内可预泄水库水量Wout；
步骤406、根据公式
Figure BDA0002279712650000104
计算小雨及小雨以下降水预报等级下第二限制水位值
Figure BDA0002279712650000105
其中，Zl为水库的原防洪限制水位，V(Zl)为水库的原防洪限制水位Zl对应的库容；
步骤407、比较小雨及小雨以下降水预报等级下第一限制水位值
Figure BDA0002279712650000106
和小雨及小雨以下降水预报等级下第二限制水位值
Figure BDA0002279712650000107
的数值大小，当时，小雨及小雨以下降水预报等级的水库防洪限制水位当
Figure BDA00022797126500001010
时，小雨及小雨以下降水预报等级的水库防洪限制水位Zx取和
Figure BDA00022797126500001012
中的较小值；
步骤五、中雨及中雨以上降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤501、根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，根据公式计算当前预报量级下未来24小时越级水量
Figure BDA00022797126500001014
其中，为未来24小时内对应降水预报等级的越级控制雨量上限阈值，单位为mm；根据公式
Figure BDA00022797126500001016
计算当前预报量级下第一限制水位值
当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，且当前预报量级下第一限制水位值
Figure BDA00022797126500001018
高于防洪限制水位Zl时，取中雨及中雨以上降水预报等级下第一限制水位值
Figure BDA00022797126500001019
步骤502、根据公式
Figure BDA00022797126500001020
计算中雨及中雨以上降水预报等级下未来24小时水量
Figure BDA00022797126500001021
其中，
Figure BDA00022797126500001022
为未来24小时内对应降水预报等级的越级控制雨量下限阈值，单位为mm；
步骤503、根据公式
Figure BDA0002279712650000111
计算中雨及中雨以上降水预报等级下未来24小时越级控制雨量对应的平均来流
步骤504、根据公式
Figure BDA0002279712650000113
计算中雨及中雨以上降水预报等级下未来24小时内超出兴利用水的平均流量Qa；
步骤505、根据公式Wa＝Qa×86400，计算中雨及中雨以上降水预报等级下未来24小时内超出流量部分对应的有效水量Wa；
步骤506、根据公式
Figure BDA0002279712650000114
计算中雨及中雨以上降水预报等级下第二限制水位值
Figure BDA0002279712650000115
步骤507、天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，比较
Figure BDA0002279712650000116
和
Figure BDA0002279712650000117
的数值大小，当
Figure BDA0002279712650000118
与
Figure BDA0002279712650000119
相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure BDA00022797126500001110
当与
Figure BDA00022797126500001112
不相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位取
Figure BDA00022797126500001114
与
Figure BDA00022797126500001115
中的较大值；再比较中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure BDA00022797126500001116
和原防洪限制水位Zl的数值大小，当时，中雨及中雨以上降水预报等级的水库防洪限制水位
Figure BDA00022797126500001118
当
Figure BDA00022797126500001119
时，中雨及中雨以上降水预报等级的水库防洪限制水位Zzd＝Zl。
需要说明的是，针对不同降水预报等级确定不同的水库防洪限制水位动态调节方式，根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为小雨及小雨以下降水预报等级时，以洪水资源利用为主要目标，按预设的越级控制雨量阈值防控洪水风险；当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，以控制和降低洪水风险为主要目标，按预设的越级控制雨量阈值上限降低调洪风险，按预设的越级控制雨量阈值下限减少兴利损失，在更大程度上降低了水库的调洪风险；针对不同降水预报等级，均计算两个限制水位值，其中，不同降水预报等级下的第一限制水位值的制定主要是站在洪量控制的角度上，计算过程均以防洪高水位或设计洪水位为基础进行，考虑漏报上限，未来可能入库水量有多少，就将库容从防洪高水位或设计洪水位向下浮动多少；不同降水预报等级下的第二限制水位值的拟定则主要是基于洪水资源利用或减少兴利库容损失的要求，计算过程以常规防洪高水位为基础进行，在小雨及小雨以下降水预报等级时，不确定性主要考虑漏报，在不影响下游防洪目标的前提下水库未来24小时可安全下泄多少水量，就将库蓄水量向上抬升多少；在中雨及中雨以上降水预报等级时，不确定性主要考虑空报，水库未来24小时最少可能会回充多少水量不会显著影响到水库蓄水兴利，就将库蓄水量向下浮动多少。
实际实施时，不同降水预报等级下的水库防洪限制水位均选小值，从而实现对小雨及小雨以下降水预报等级的情况以兴利为主要目的，控制防洪风险；对中雨及中雨以上降水预报等级的情况以防洪为主要目的，控制兴利损失。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (5)

1.基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于，该方法包括以下步骤：
步骤一、降水预报等级划分：对天气预报播报的降雨信息，按照流域或区域面雨量大小进行降水预报等级划分，所述降水预报等级由低至高依次包括无雨、小雨、中雨、大雨、暴雨和大暴雨，所述无雨按预报降水概率由低至高依次包括晴天、多云和阴天，所述无雨的实际降水区间为[0，0.1mm)，所述小雨的实际降水区间为[0.1mm，6mm)，所述中雨的实际降水区间为[6mm，15mm)，所述大雨的实际降水区间为[15mm，30mm)，所述暴雨的实际降水区间为[30mm，60mm)，所述大暴雨的实际降水区间为[60mm，∞)；
步骤二、确定越级控制雨量阈值：流域内未来24小时降水预报为晴天间多云时，其对应的越级为小雨及小雨以上降水预报等级，其对应的越级控制雨量阈值取3mm；
流域内未来24小时降水预报为多云间阴天时，其对应的越级为中雨及中雨以上降水预报等级，其对应的越级控制雨量阈值取10mm；
流域内未来24小时降水预报为小雨时，其对应的越级为大雨及大雨以上降水预报等级，其对应的越级控制雨量阈值取15mm；
流域内未来24小时降水预报为中雨时，其对应的越级为无雨或暴雨及暴雨以上降水预报等级，其对应的越级控制雨量上限阈值取30mm，其对应的越级控制雨量下限阈值取3mm；
流域内未来24小时降水预报为大雨及大雨以上降水预报等级时，其对应的越级为中雨及中雨以下降水预报等级，其对应的越级控制雨量下限阈值取10mm；
步骤三、获取未来24小时降水预报等级信息并判断未来24小时降水预报等级是否为小雨及小雨以下降水预报等级：根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为小雨及小雨以下降水预报等级时，执行步骤四；否则，执行步骤五；
步骤四、小雨及小雨以下降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤401、根据公式
Figure FDA0002279712640000021
计算小雨及小雨以下降水预报等级下未来24小时水量
Figure FDA0002279712640000022
其中，W'为水库24小时汇流区河道槽蓄水量，
Figure FDA0002279712640000023
为未来24小时内对应降水预报等级的越级控制雨量阈值，单位为mm，α为流域内径流系数，F为水库24小时汇流区面积，单位为km2；
步骤402、根据公式
Figure FDA0002279712640000024
计算小雨及小雨以下降水预报等级下第一限制水位值
Figure FDA0002279712640000025
其中，f[·]为水库水位库容关系函数，Zh为防洪高水位，V(Zh)为防洪高水位Zh下对应的防洪库容，Qu为水库生态基流常流量；
步骤403、根据公式
Figure FDA0002279712640000026
计算小雨及小雨以下降水预报等级下未来24小时越级控制雨量对应的平均来流
Figure FDA0002279712640000027
步骤404、根据公式
Figure FDA0002279712640000028
计算小雨及小雨以下降水预报等级下未来24小时内平均可附加下泄流量
Figure FDA0002279712640000029
其中，Qs为水库安全允许泄流量；
步骤405、根据公式计算小雨及小雨以下降水预报等级下未来24小时内可预泄水库水量Wout；
步骤406、根据公式计算小雨及小雨以下降水预报等级下第二限制水位值其中，Zl为水库的原防洪限制水位，V(Zl)为水库的原防洪限制水位Zl对应的库容；
步骤407、比较小雨及小雨以下降水预报等级下第一限制水位值
Figure FDA00022797126400000213
和小雨及小雨以下降水预报等级下第二限制水位值的数值大小，当
Figure FDA00022797126400000215
时，小雨及小雨以下降水预报等级的水库防洪限制水位当
Figure FDA00022797126400000217
时，小雨及小雨以下降水预报等级的水库防洪限制水位Zx取
Figure FDA00022797126400000218
和
Figure FDA00022797126400000219
中的较小值；
步骤五、中雨及中雨以上降水预报等级的水库防洪限制水位动态调节，过程如下：
步骤501、根据天气预报获取未来24小时降水预报等级信息，当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，根据公式
Figure FDA0002279712640000031
计算当前预报量级下未来24小时越级水量其中，
Figure FDA0002279712640000033
为未来24小时内对应降水预报等级的越级控制雨量上限阈值，单位为mm；根据公式计算当前预报量级下第一限制水位值
Figure FDA0002279712640000035
当天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，且当前预报量级下第一限制水位值
Figure FDA0002279712640000036
高于防洪限制水位Zl时，取中雨及中雨以上降水预报等级下第一限制水位值
Figure FDA0002279712640000037
步骤502、根据公式
Figure FDA0002279712640000038
计算中雨及中雨以上降水预报等级下未来24小时水量
Figure FDA0002279712640000039
其中，
Figure FDA00022797126400000310
为未来24小时内对应降水预报等级的越级控制雨量下限阈值，单位为mm；
步骤503、根据公式
Figure FDA00022797126400000311
计算中雨及中雨以上降水预报等级下未来24小时越级控制雨量对应的平均来流
Figure FDA00022797126400000312
步骤504、根据公式
Figure FDA00022797126400000313
计算中雨及中雨以上降水预报等级下未来24小时内超出兴利用水的平均流量Qa；
步骤505、根据公式Wa＝Qa×86400，计算中雨及中雨以上降水预报等级下未来24小时内超出流量部分对应的有效水量Wa；
步骤506、根据公式
Figure FDA00022797126400000314
计算中雨及中雨以上降水预报等级下第二限制水位值
Figure FDA00022797126400000315
步骤507、天气预报播报未来24小时降水预报等级为中雨及中雨以上降水预报等级时，比较和
Figure FDA00022797126400000317
的数值大小，当
Figure FDA00022797126400000318
与
Figure FDA00022797126400000319
相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure FDA00022797126400000320
当
Figure FDA00022797126400000321
与
Figure FDA00022797126400000322
不相等时，中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure FDA00022797126400000323
取
Figure FDA00022797126400000324
与
Figure FDA00022797126400000325
中的较大值；再比较中雨及中雨以上降水预报等级的估计水库防洪限制水位
Figure FDA00022797126400000326
和原防洪限制水位Zl的数值大小，当时，中雨及中雨以上降水预报等级的水库防洪限制水位
Figure FDA0002279712640000041
当
Figure FDA0002279712640000042
时，中雨及中雨以上降水预报等级的水库防洪限制水位Zzd＝Zl。
2.按照权利要求1所述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述水库为暴雨洪水入流类型且承担防洪任务的水库，所述暴雨洪水入流类型水库包括供水水库、灌溉水库、发电水库及多目标综合利用水库；
当水库为供水水库或灌溉水库时，水库生态基流常流量Qu为满足河道下游生态及航运要求的下泄流量；
当水库为发电水库或多目标综合利用水库时，水库生态基流常流量Qu取发电下泄等常流量和下游生态流量两者中的大值。
3.按照权利要求2所述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述水库为单一水库。
4.按照权利要求1所述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述流域为扇形流域或汇流面积不大于50000km2的羽形流域。
5.按照权利要求1所述的基于预报误差越级设防的水库防洪限制水位动态调节方法，其特征在于：所述流域内径流系数α取值范围为0.6～0.8。
CN201911136441.3A 2019-11-19 2019-11-19 基于预报误差越级设防的水库防洪限制水位动态调节方法 Active CN110852525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911136441.3A CN110852525B (zh) 2019-11-19 2019-11-19 基于预报误差越级设防的水库防洪限制水位动态调节方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911136441.3A CN110852525B (zh) 2019-11-19 2019-11-19 基于预报误差越级设防的水库防洪限制水位动态调节方法

Publications (2)

Publication Number Publication Date
CN110852525A true CN110852525A (zh) 2020-02-28
CN110852525B CN110852525B (zh) 2023-04-18

Family

ID=69602654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911136441.3A Active CN110852525B (zh) 2019-11-19 2019-11-19 基于预报误差越级设防的水库防洪限制水位动态调节方法

Country Status (1)

Country Link
CN (1) CN110852525B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111476676A (zh) * 2020-04-02 2020-07-31 中国水利水电科学研究院 一种城镇小流域分类分级智能调度方法
CN111612292A (zh) * 2020-04-03 2020-09-01 贵州乌江水电开发有限责任公司 基于关键水位控制的梯级水电站调度控制系统及方法
CN111650572A (zh) * 2020-07-23 2020-09-11 浪潮云信息技术股份公司 一种减小短时降水估测偏差的方法及系统

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1610156A1 (fr) * 2004-03-31 2005-12-28 Centre National Du Machinisme Agricole, Du Genie Rural, Des Eaux Et Des Forets (Cemagref) Dispositif d'alerte pour au moins un phénomène météorologique tel que les pluies, les écoulements de crues de cours d'eau, pour une zone géographique donnée
CN102599118A (zh) * 2012-04-11 2012-07-25 宁德市富发水产有限公司 采捕海水池塘中人工培养桡足类的方法
CN104179148A (zh) * 2014-04-29 2014-12-03 河海大学 一种基于不可控洪水特性的洪水调度安全度评价方法
CN105825309A (zh) * 2016-03-24 2016-08-03 大连理工大学 一种协调水库汛期防洪风险与兴利效益的调度方法
CN106682789A (zh) * 2017-01-09 2017-05-17 泰华智慧产业集团股份有限公司 规划防汛物资储备数量的方法及系统
CN106874657A (zh) * 2017-01-18 2017-06-20 贵州黔源电力股份有限公司 一种汛期常遇洪水下的水库水位控制方法
CN107059761A (zh) * 2017-06-19 2017-08-18 武汉大学 水库群防洪库容时空分配设计方法
CN108681848A (zh) * 2018-08-25 2018-10-19 黄河水利委员会黄河水利科学研究院 一种“三位一体”小型水库或淤地坝汛期预警方法
CN108985585A (zh) * 2018-06-28 2018-12-11 河海大学 一种考虑预报不确定性影响的水库洪水资源利用风险对冲决策方法
CN109829611A (zh) * 2018-12-24 2019-05-31 长江勘测规划设计研究有限责任公司 基于防洪库容动态分配的梯级优化调度方法
CN109829584A (zh) * 2019-02-02 2019-05-31 中国水利水电科学研究院 一种山洪风险动态评估方法
CN110288239A (zh) * 2019-06-27 2019-09-27 华中科技大学 基于自适应分级预泄法的汛限水位动态控制方法及系统
CN110363356A (zh) * 2019-07-17 2019-10-22 长安大学 一种面向生态的水土资源时空优化配置方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1610156A1 (fr) * 2004-03-31 2005-12-28 Centre National Du Machinisme Agricole, Du Genie Rural, Des Eaux Et Des Forets (Cemagref) Dispositif d'alerte pour au moins un phénomène météorologique tel que les pluies, les écoulements de crues de cours d'eau, pour une zone géographique donnée
CN102599118A (zh) * 2012-04-11 2012-07-25 宁德市富发水产有限公司 采捕海水池塘中人工培养桡足类的方法
CN104179148A (zh) * 2014-04-29 2014-12-03 河海大学 一种基于不可控洪水特性的洪水调度安全度评价方法
CN105825309A (zh) * 2016-03-24 2016-08-03 大连理工大学 一种协调水库汛期防洪风险与兴利效益的调度方法
CN106682789A (zh) * 2017-01-09 2017-05-17 泰华智慧产业集团股份有限公司 规划防汛物资储备数量的方法及系统
CN106874657A (zh) * 2017-01-18 2017-06-20 贵州黔源电力股份有限公司 一种汛期常遇洪水下的水库水位控制方法
CN107059761A (zh) * 2017-06-19 2017-08-18 武汉大学 水库群防洪库容时空分配设计方法
CN108985585A (zh) * 2018-06-28 2018-12-11 河海大学 一种考虑预报不确定性影响的水库洪水资源利用风险对冲决策方法
CN108681848A (zh) * 2018-08-25 2018-10-19 黄河水利委员会黄河水利科学研究院 一种“三位一体”小型水库或淤地坝汛期预警方法
CN109829611A (zh) * 2018-12-24 2019-05-31 长江勘测规划设计研究有限责任公司 基于防洪库容动态分配的梯级优化调度方法
CN109829584A (zh) * 2019-02-02 2019-05-31 中国水利水电科学研究院 一种山洪风险动态评估方法
CN110288239A (zh) * 2019-06-27 2019-09-27 华中科技大学 基于自适应分级预泄法的汛限水位动态控制方法及系统
CN110363356A (zh) * 2019-07-17 2019-10-22 长安大学 一种面向生态的水土资源时空优化配置方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SU-XIANG ZHANG; MING MEI; XIN-FA XU; LI-MO TANG; NIAN-HUA WU; SHA-SHA FU; XI LI: "Numerical modeling of water level and current in very shallow large Lake Poyang" *
刘招，黄强，于兴杰，王丽霞: "基于6h预报径流深的安康水库防洪预报调度方案研究" *
王国利;梁国华;何斌;王本德;: "基于分级预泄条件的水库汛限水位动态控制域与应用" *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111476676A (zh) * 2020-04-02 2020-07-31 中国水利水电科学研究院 一种城镇小流域分类分级智能调度方法
CN111476676B (zh) * 2020-04-02 2023-11-21 中国水利水电科学研究院 一种城镇小流域分类分级智能调度方法
CN111612292A (zh) * 2020-04-03 2020-09-01 贵州乌江水电开发有限责任公司 基于关键水位控制的梯级水电站调度控制系统及方法
CN111650572A (zh) * 2020-07-23 2020-09-11 浪潮云信息技术股份公司 一种减小短时降水估测偏差的方法及系统
CN111650572B (zh) * 2020-07-23 2023-02-24 浪潮云信息技术股份公司 一种减小短时降水估测偏差的方法及系统

Also Published As

Publication number Publication date
CN110852525B (zh) 2023-04-18

Similar Documents

Publication Publication Date Title
CN110852525B (zh) 基于预报误差越级设防的水库防洪限制水位动态调节方法
CN108108838B (zh) 一种高水量利用率的季调节水库优化调度方法
CN110288239B (zh) 基于自适应分级预泄法的汛限水位动态控制方法及系统
CN106202978B (zh) 多泥沙河流综合利用水库分期汛限水位优化方法及其系统
CN107992989B (zh) 一种水库群洪水资源化预蓄预泄风险决策方法
CN108596417A (zh) 一种基于库容补偿法的梯级水库洪水资源化利用方法
CN106779232B (zh) 一种城市内涝的建模预测方法
CN110991688B (zh) 基于气象数值预报的水库调度预警方法
CN109272186B (zh) 一种水资源调度方法
Chang et al. Floodwater utilization for cascade reservoirs based on dynamic control of seasonal flood control limit levels
CN107679021A (zh) 一种河流入湖库流量计算方法
CN107180318B (zh) 基于变权重剩余防洪库容最大的水库群防洪库容分配方法
CN116149187B (zh) 一种用于水库库尾鱼类栖息地生境保护的生态调度方法
US20220279741A1 (en) Decision-making Method for Variable Rate Irrigation Management
CN108537449A (zh) 计及河道度汛需求的水库协调调度策略获取方法
CN110414051A (zh) 一种抑制河流水华的生态需水核算方法
CN109002946B (zh) 一种河湖补水的“两库—两站”系统水资源优化调度方法
CN117436619B (zh) 一种基于等效防洪作用的梯级水库防洪库容联合预留方法
CN108090671B (zh) 一种多阶段域间-域内反馈迭代的跨流域负荷分配方法
CN110378805B (zh) 梯级水电站发电调度分级调峰方法
CN116402282A (zh) 一种基于水量平衡的水库水位调控方法及系统
CN111915173A (zh) 一种梯级小水电站优化调度关键因素的筛选方法
CN115994602A (zh) 水库汛期水位动态控制方法及水库汛期水位动态控制设备
CN114707879A (zh) 一种基于预泄规则下的汛期分期动态汛限水位确定方法
CN111126693B (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