CN113947315A - 一种水库分级分期旱限水位设计与控制方法 - Google Patents
一种水库分级分期旱限水位设计与控制方法 Download PDFInfo
- Publication number
- CN113947315A CN113947315A CN202111230658.8A CN202111230658A CN113947315A CN 113947315 A CN113947315 A CN 113947315A CN 202111230658 A CN202111230658 A CN 202111230658A CN 113947315 A CN113947315 A CN 113947315A
- Authority
- CN
- China
- Prior art keywords
- water
- reservoir
- drought
- water level
- water supply
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 532
- 238000000034 method Methods 0.000 title claims abstract description 90
- 238000013461 design Methods 0.000 title claims abstract description 40
- 230000008569 process Effects 0.000 claims abstract description 43
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 230000001932 seasonal effect Effects 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000003860 storage Methods 0.000 claims description 13
- 239000008235 industrial water Substances 0.000 claims description 7
- 239000003621 irrigation water Substances 0.000 claims description 7
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical group [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 6
- 230000018109 developmental process Effects 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 239000011780 sodium chloride Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 2
- 238000001704 evaporation Methods 0.000 claims description 2
- 238000012946 outsourcing Methods 0.000 claims description 2
- 238000004321 preservation Methods 0.000 claims 1
- 238000011160 research Methods 0.000 description 7
- 238000007689 inspection Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000009776 industrial production Methods 0.000 description 4
- 238000003973 irrigation Methods 0.000 description 4
- 230000002262 irrigation Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000012271 agricultural production Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 239000008400 supply water Substances 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/06313—Resource planning in a project environment
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Sewage (AREA)
- Alarm Systems (AREA)
Abstract
本发明提供一种水库分级分期旱限水位设计与控制方法，属于抗旱调度技术领域。首先，分析供水保障目标及设计供水需求。其次，通过审查长系列水文资料的可靠性、一致性、代表性，设计不同预警级别的入库径流过程；综合考虑入库径流、水库运行水位、供水需求等因素的季节性特点，对水文年进行干旱预警分期。再次，根据设计需水过程、供水保障优先级，确定分级分期旱限水位供水规则，并初步确定分级分期的旱限水位。最后，基于长系列径流过程模拟进行分级分期旱限水位合理性分析与调整。本发明提出的基于长系列模拟的合理性分析方法及水库旱限水位反馈修正方法，能够综合考虑工程应用中对系统供水保证率与破坏深度的要求，计算快速、便捷、适于工程推广应用。
Description
技术领域
本发明属于抗旱调度技术领域，涉及到一种水库分级分期旱限水位设计与控制方式。
背景技术
受气候变化与强人类活动影响，近年来我国极端干旱事件趋多趋强。同时，随着经济社会的快速发展，用水需求迅速增加，水资源供需矛盾愈发突出，严重制约了我国经济社会发展(屈艳萍,吕娟,苏志诚,等.抗旱减灾研究综述及展望[J].水利学报,2018,49(01):115-125.)。水库是重要的蓄水工程，制定水库抗旱特征指标，指导水资源合理分配，有助于降低缺水事件的发生频率、减轻缺水的严重程度，对实现抗旱减灾具有重要意义。2011年，国家防汛抗旱总指挥部办公室在《旱限水位(流量)确定办法》(国家防洪抗旱总指挥部办公室.旱限水位(流量)的确定办法[Z].2011.)中首次提出以旱限水位作为水库抗旱关键性控制水位，并给出了相应的计算方法。该办法首先通过逐月滑动计算得到各月水库应供水量，再将各月水库应供水量与死库容之和最大值对应的水库水位确定为旱限水位。
该方法全年设置单一旱限水位。而水库来水、蓄水、用水均存在季节特性，年内不同时期枯水具有不同特点。全年设置相同的旱限水位对于年内某些时期可能偏高而造成不必要的供水破坏，对于年内某些时期可能偏低使得预蓄水量不足。为此，刘攀等(刘攀,李立平,吴荣飞,等.论水库旱限水位分期控制的必要性与计算方法探讨[J].水资源研究,2012(1):51-55.)首次提出对旱限水位进行分期控制的理念，之后宋树东(宋树东,朱文才.水库旱限水位分期确定的研究[J].长春理工大学学报(自然科学版),2014(3):160-163.)提出了分期旱限水位滑动计算方法。该方法简单易用，但在计算旱限水位后，未能进行合理性分析，忽略了设置旱限水位对供水系统的影响，不能确保初步设置旱限水位的调控效果。近年来，一些学者从不同角度定义了旱限水位，基于模拟-优化思想提出了相应的计算方法，从而考虑了设置旱限水位对供水系统的影响，保证了设置旱限水位后系统供水性能符合预期，进行有效的旱灾防御。张礼兵等(张礼兵,伍露露,金菊良,等.大型灌区骨干水库分期旱限水位研究[J].水利学报,2018,49(6):757-766.)将旱限水位定义为水库水位较低时，应削减供水的临界水位，提出在满足各行业供水保证率前提下，使系统供水经济效益最大化的旱限水位优化计算方法。彭少明等(彭少明,王煜,张永永,等.多年调节水库旱限水位优化控制研究[J].水利学报,2016,47(4):552-559.)针对多年调节水库跨年度补水问题，将旱限水位定义为水库年末应蓄水量对应的水位，提出年际缺水均衡控制的旱限水位优化计算方法，有效减小了连枯水年供水破坏深度。Chang等人(Chang J,Guo A,Wang Y,etal.Reservoir Operations to Mitigate Drought Effects With a Hedging PolicyTriggered by the Drought Prevention Limiting Water Level[J].Water ResourcesResearch,2019,55(2):904-922.)将旱限水位定义为不同程度干旱发生时水库的最低运行水位，综合历史水文气象情况与使系统总缺水程度最小的水位过程，确定旱限水位。
然而，以上研究在确定旱限水位时仅考虑了供水保证率或系统整体缺水程度中某一方面供水要求，没有全面考虑工程应用中对系统供水保证率、破坏深度的要求。同时，上述旱限水位的确定均依赖优化方法，在推广应用时有较高的技术门槛与难度。为此，本发明全面考虑系统供水保证率、破坏深度要求，提出一种快速、便捷，适于工程推广的旱限水位简便计算方法，为水利部门指导水库抗旱运行，制定抗旱预案提供技术支撑。
发明内容
针对现有技术的不足，本发明提供了一种水库分级分期旱限水位设计与控制方法。
本发明采用的技术方案如下：
一种水库分级分期旱限水位设计与控制方法，流程图见附图1，其主要包括以下步骤：
步骤一：分析供水保障目标及设计供水需求。
(1)分析规划水平年水库生态环境保护和经济社会保障目标，包括供水对象、供水时段、供水范围及其供水保障优先级。
(2)计算河道内、河道外各用户规划水平年需水量和需水过程，并明确各用水户设计供水保证率、供水破坏深度要求。河道外需水包括工业、农业、生活经济社会发展需水与生态环境需水，总量为各用户需水量之和；河道内需水包括生态、航运、生产等，总量按一水多用的原则，取各用户需水量的外包值。河道外经济社会发展需水计算参考GBT 51051-2014《水资源规划规范》，生态环境需水计算参考SL 429-2008《水资源供需预测分析技术规范》。河道内生态环境需水计算参考SL/Z 712《河湖生态环境需水计算规范》，交通航运需水计算参考GB 50139《内河通航标准》，其他生产需水量根据生产过程对流量、流速、水位等要求综合确定。
步骤二：设计不同预警级别的入库径流过程。
审查长系列水文资料的可靠性、一致性、代表性，设计不同预警级别的设计枯水年(系列)径流过程，包括水文频率为75％的一般枯水年(系列)径流过程和水文频率为95％的特枯年(系列)径流过程。
上述设计枯水年径流过程推求方法分为以下两种情况：
(1)对于年调节、季调节水库，采用水文频率分析方法设计枯水年径流量及其年内分配，径流频率曲线线型、统计参数及代表年的选取参考SL/T 278—2020《水利水电工程水文计算规范》。
(2)对于多年调节水库，采用扣除允许破坏年数法确定设计枯水系列。具体为：
首先，按式(1)计算设计保证率条件下正常工作允许破坏的年数Tf。然后，在实测资料中选出最枯连续枯水年组，从该最枯连续枯水年组最后一年起逆时序扣除允许破坏年数，余下的即为初步确定的设计枯水系列。本发明采用连续枯水年组的平均年径流量、年径流量变差系数、连续枯水年数三个指标评估枯水的严重程度，三个指标按重要性排序依次是：平均年径流量、连续枯水年数、年径流量变差系数。
Tf＝n-Pd(n+1) (1)
式中，n表示水文系列总年数，Pd为设计保证率。
步骤三：确定干旱预警分期。
综合考虑入库径流、水库运行水位(库容)、供水需求等因素的季节性特点，采用Fisher最优分割法、成因分析法、模糊分析法、分形分析法等方法对水文年进行干旱预警分期，如图2所示。
步骤四：确定分级分期的旱限水位供水规则。
根据步骤一确定的设计需水过程、供水保障优先级，确定分级分期旱限水位供水规则，如附图3所示。
水库旱警水位指水库水位持续偏低，影响城乡生活、工农业生产、生态环境等用水安全，应采取一般抗旱应急措施的水位。该水位是能够在一般枯水年份(一般取年径流频率75％)保障城乡生活需水量、工业生产需水量、重要农业生产需水量、河道生态环境需水量的水库水位。
(1)当水库水位高于旱警水位时，水库采取“充分供水”的供水措施，不限制供水，“充分供水”供水量的计算公式如下所示：
Wa,t＝Dagr,t+Dind,t+Dlife,t+max{Deco,t,Dship,t} (2)
式中，Wa,t为水库t时段“充分供水”供水量；Dagr,t为水库t时段正常灌溉需水量；Dind,t为水库t时段正常工业需水量；Dlife,t为水库t时段正常生活需水量；Deco,t为水库t时段生态需水量；Dship,t为水库t时段航运需水量。航运需水和生态需水都属于河道内需水，可共用，两者取大。
水库旱保水位指水库水位显著偏低，严重影响城乡生活、工农业生产、生态环境等用水安全，应采取高级别抗旱应急措施的水位。该水位指能够在特枯年份(一般取年径流频率95％)保障城乡生活需水量、重要工业需水量、重要农业需水量、河道基本生态环境需水量的水库水位。
(2)当水库水位在旱警水位以下，旱保水位以上时，水库采取“一级限供”的供水措施，轻度限制供水，“一级限供”供水量的计算公式如下所示：
Wb,t＝Dagr,t·vagr,t+Dind,t+Dlife,t+max{Deco,t,Dship,t} (3)
式中，vagr,t为水库t时段农业限制供水系数。
(3)当水库水位在旱保水位以下时，水库采取“二级限供”的供水措施，进一步限制供水，“二级限供”供水量的计算公式如下所示：
Wc,t＝Dagr,t·vagr,t+Dind,t·vind,t+Dlife,t+max{Deco,t,Dship,t} (4)
式中，vind,t为水库t时段工业限制供水系数。
上式中的工业、农业限制供水系数均需结合各用水户的允许破坏深度与当地抗旱应急预案，综合分析确定，也可直接采用用水户允许破坏深度，对用水需求进行折减。
步骤五：初步确定分级分期旱限水位。
基于步骤一确定的设计需水过程、步骤二确定的设计入库径流过程，以及步骤四确定的旱限水位供水规则，采用逆序递推方法初步确定分级分期旱限水位。计算旱警水位时，以一般枯水年(系列)径流过程及“充分供水”供水过程作为输入，计算旱保水位时，以特枯年(系列)径流过程及“一级限供”供水过程作为输入。
以供水期末水位恰好达到水库死水位或允许消落的最低水位(见公式6)，以旬或月为计算时段，依据水库兴利调节原理，逆序递推得到各月(旬)初水位，递推公式如公式(5)所示。
逆序递推过程中，水库水位要满足基本约束：库水位在死水位和兴利蓄水上限水位之间，其中汛期兴利蓄水上限水位为水库汛限水位，非汛期兴利蓄水上限水位为正常蓄水位，见式(7)；各月(旬)旱警水位不应低于同期旱保水位。满足以上约束的各月(旬)初水位即为各月(旬)旱限水位和旱警水位。
ZT+1＝Zdead (6)
式中，Zhj,t、Zhj,t+1分别为水库t、t+1月(旬)的旱警水位；f()为水库库容-水位曲线；f'()为水库水位-库容曲线；Wa,t为t月(旬)“充分供水”供水量；Wb,t为t月(旬)“一级限供”供水量；Wloss,t为水库t月(旬)蒸发、渗漏损失水量；Wp1,t、Wp2,t分别为一般枯水年和特枯水年t月(旬)的水库来水量；Zhb,t、Zhb,t+1别为水库t、t+1月(旬)的旱保水位；ZT+1为水文年末水位；Zdead为水库死水位或考虑水质要求允许取水的最低水位；Zlimit为水库汛限水位；Znormal为水库正常蓄水位。
在上一步推求的各月(旬)旱限水位和旱警水位的基础上，基于步骤三确定的干旱预警分期，取各分期内最低月(旬)旱限水位、旱警水位作为该分期旱限水位、旱警水位。
步骤六：分级分期旱限水位合理性分析与调整。
具体包括如下内容：
(1)基于步骤四确定的旱限水位供水规则以及步骤五确定的分级分期旱限水位，对历史长系列径流过程进行模拟，统计各用户供水保证率Rel、最大供水破坏深度Vul，计算公式如下：
Vul＝max(Deft/Demt) (9)
式中，N为总的统计时段数；Zt为t时段供水状态，若当前时段满足供水需求，则Zt＝1，若当前时段未满足供水需求，即发生供水破坏，Zt＝0；Deft为时段t的缺水量；Demt为时段t的需水量。
(2)判别供水保证率、破坏深度是否符合规划设计要求：若各用水户供水保证率、供水破坏深度指标均满足规划设计要求，则说明初步确定的旱限水位科学合理。若有任意指标不满足规划设计要求，则需要对初步确定的旱限水位进行反馈修正与调整，直至指标达标。
(3)分级分期旱限水位反馈修正方法：
若供水保证率不达标，首先采用更换典型枯水年(系列)方式进行反馈修正，接着按照“据旱限水位年内变化规律，调整对系统供水保证率、破坏深度影响大的旱限水位”思路进行反馈修正。本发明提出提升高用水时期旱警水位、降低枯水期旱保水位的反馈修正方法。
若供水保证率达标而供水破坏深度不达标，则采用提升高用水时期旱警水位、降低枯水期旱保水位的反馈修正方法。
若无法通过修正使得各用水户供水保证率、供水破坏深度指标达标，则说明水库供水系统无法承担当前设计需水量，应对设计需水量进行调整，直至各用水户供水保证率、供水破坏深度达到设计要求。
与现有技术相比，本发明具有以下优点和效果：
本发明综合考虑各行业供水的保障优先级、保证率、破坏深度以及入库径流特性，结合干旱径流设计技术、水库模拟调度技术，制定了科学合理的水库分级分期旱限水位确定方法，提出了基于长系列径流过程模拟的合理性分析方法及水库分级分期旱限水位反馈修正方法。该方法综合考虑了工程应用中对系统供水保证率与破坏深度的要求，计算快速、便捷、适于工程推广应用，为水利部门制定水库旱限水位，指导水库抗旱调度提供了一套科学便捷的方法。
附图说明
图1是水库分级分期旱限水位计算流程图。
图2是水库干旱预警分期确定示意图。
图3是水库分级分期旱限水位供水规则示意图。
图4是清河水库旱警(保)水位反馈调整前后对比图。
具体实施方式
本发明以清河水库为例，结合技术方案和附图详细叙述具体实施方式。具体包括以下步骤：
第一步，供水保障目标及供水需求分析。
清河水库近景规划水平年供水任务包括工业生产供水、城市生活供水、农业灌溉及生态环境供水，年需水量分别为30百万m3、20百万m3、300百万m3、32百万m3。其中，农业灌溉设计保证率为75％，工业供水设计保证率为95％。农业灌溉供水破坏深度一般情况下应不超过0.3，极端枯水时应不超过0.5，工业、城乡生活供水破坏深度应不超过0.1。各用户规划水平年需水过程如表1所示。
第二步，设计不同预警级别的入库径流。
采用Mann-Kendall趋势检验对清河水库1950-2017共68年入库年径流量序列作趋势检验，采用Pettitt变点检验法作突变点检验。检验结果表明，清河水库入库径流符合一致性要求。采用多年调节水库设计枯水系列方法确定清河水库不同预警级别的枯水系列过程。扣除17年最严重的连续枯水年组，对应水文频率为75％的一般枯水系列为2006～2007年(水文年)。扣除3年最严重的连续枯水年组，对应水文频率为95％的特枯水系列为1999～2001年(水文年)。
第三步，确定干旱预警分期。
综合考虑清河流域多年月平均降水量、入库流量、库水位、供水量，本发明采用Fisher最优分割法对清河水库进行干旱预警分期，将一个水文年分为5个时期，6月为灌溉期，7～8月为主汛期，9～10月为汛后蓄水期，11月～翌年4月为枯水期，5月为灌溉用水高峰期。
第四步，确定分级分期旱限水位供水规则。
清河水库农业灌溉取水口在水库下游，农业灌溉用水与河道生态环境用水可共用，故当清河水库农业灌溉供水流量高于生态流量时，清河水库不需要额外给生态环境供水，否则按生态流量要求进行生态补水。当水库采取“一级限供”的供水措施时，限制农业供水，保障生活、工业、生态供水。当水库采取“二级限供”的供水措施时，限制农业、工业供水，保障生活、生态供水。各用水户限制供水量由允许破坏深度确定，农业供水折减系数取0.7，工业供水折减系数取0.9。清河水库不同预警级别设计供水过程如表1所示。
表1清河水库不同预警级别设计供水过程表单位：百万m3
第五步，初步确定分级分期旱限水位。
各月旱警水位：一般枯水系列2006～2007年(水文年)，采用“充分供水”的供水过程，由死水位逆序递推得到各月初水位，枯水系列中各月最低水位为该月旱警水位。各月旱保水位：特枯水系列1999～2001(水文年)，采用为“一级限供”供水过程，由死水位逆序递推得到各月初水位，特枯枯水系列中各月最低水位为该月旱保水位。初步确定的清河水库各月旱警(保)水位如表2所示。
表2初步确定的清河水库各月旱警(保)水位表单位：m
各分期旱警(保)水位为该分期最低月旱警(保)水位，初步确定的清河水库分期旱警(保)水位如表3所示。
表3初步确定的清河水库各分期旱警(保)水位表单位：m
第六步，分级分期旱限水位合理性分析与调整。
在初步确定的旱警(保)水位控制运行下，对清河水库进行长系列径流过程进行供水模拟，各项供水指标如表4所示。虽然各用水户供水保证率满足设计要求，但农业供水最大破坏深度与工业供水最大破坏深度均不达标。按照本发明所提方法，当供水保证率达标而破坏深度不达标时，采用提高需水高峰期旱警水位方式进行反馈修正。同时，为了便于工程实际操作，将其余不参与修正的旱限水位进行四舍五入取整。经过试算，按该方法适当提高旱警水位后，各用水户供水保证率与最大破坏深度均满足要求。调整前后清河水库旱警(保)水位如图4所示。
表4长系列模拟后各项供水指标
以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。
Claims (5)
1.一种水库分级分期旱限水位设计与控制方法，其特征在于，包括以下步骤：
步骤一：分析供水保障目标及设计供水需求；
步骤二：设计不同预警级别的设计枯水年径流过程，包括一般枯水年径流过程和特枯年径流过程，设计分为以下两种情况：
(1)对于年调节、季调节水库，采用水文频率分析方法设计枯水年径流量及其年内分配，径流频率曲线线型、统计参数及代表年的选取参考SL/T 278—2020《水利水电工程水文计算规范》；
(2)对于多年调节水库，采用扣除允许破坏年数法确定设计枯水系列；具体为：
首先，按式(1)计算设计保证率条件下正常工作允许破坏的年数Tf；然后，在实测资料中选出最枯连续枯水年组，从该最枯连续枯水年组最后一年起逆时序扣除允许破坏年数，余下的即为初步确定的设计枯水系列；采用连续枯水年组的平均年径流量、年径流量变差系数、连续枯水年数三个指标评估枯水的严重程度，三个指标按重要性排序依次是：平均年径流量、连续枯水年数、年径流量变差系数；
Tf＝n-Pd(n+1) (1)
式中，n表示水文系列总年数，Pd为设计保证率；
步骤三：确定干旱预警分期；
综合考虑入库径流、水库运行水位、供水需求或其他因素的季节性特点，对水文年进行干旱预警分期；
步骤四：确定分级分期的旱限水位供水规则；
根据步骤一确定的设计需水过程、供水保障优先级，确定分级分期旱限水位供水规则；
(1)当水库水位高于旱警水位时，水库采取“充分供水”的供水措施，不限制供水，“充分供水”供水量的计算公式如下所示：
Wa,t＝Dagr,t+Dind,t+Dlife,t+max{Deco,t,Dship,t} (2)
式中，Wa,t为水库t时段“充分供水”供水量；Dagr,t为水库t时段正常灌溉需水量；Dind,t为水库t时段正常工业需水量；Dlife,t为水库t时段正常生活需水量；Deco,t为水库t时段生态需水量；Dship,t为水库t时段航运需水量；航运需水和生态需水都属于河道内需水，可共用，两者取大；
(2)当水库水位在旱警水位以下，旱保水位以上时，水库采取“一级限供”的供水措施，轻度限制供水，“一级限供”供水量的计算公式如下所示：
Wb,t＝Dagr,t·vagr,t+Dind,t+Dlife,t+max{Deco,t,Dship,t} (3)
式中，vagr,t为水库t时段农业限制供水系数；
(3)当水库水位在旱保水位以下时，水库采取“二级限供”的供水措施，进一步限制供水，“二级限供”供水量的计算公式如下所示：
Wc,t＝Dagr,t·vagr,t+Dind,t·vind,t+Dlife,t+max{Deco,t,Dship,t} (4)
式中，vind,t为水库t时段工业限制供水系数；
步骤五：初步确定分级分期旱限水位；
基于步骤一确定的设计需水过程、步骤二确定设计入库径流过程，以及步骤四确定的旱限水位供水规则，采用逆序递推方法初步确定分级分期的旱限水位；计算旱警水位时，以一般枯水年径流过程及“充分供水”供水过程作为输入；计算旱保水位时，以特枯年径流过程及“一级限供”供水过程作为输入；
以供水期末水位恰好达到水库死水位或允许消落的最低水位，以旬或月为计算时段，依据水库兴利调节原理，逆序递推得到各月/旬初水位，递推公式如公式(5)所示；
逆序递推过程中，水库水位要满足基本约束：库水位在死水位和兴利蓄水上限水位之间，其中汛期兴利蓄水上限水位为水库汛限水位，非汛期兴利蓄水上限水位为正常蓄水位，见式(7)；各月/旬旱警水位不应低于同期旱保水位；满足以上约束的各月/旬初水位即为各月/旬旱限水位和旱警水位；
ZT+1＝Zdead (6)
式中，Zhj,t、Zhj,t+1分别为水库t、t+1月/旬的旱警水位；f()为水库库容-水位曲线；f'()为水库水位-库容曲线；Wa,t为t月/旬“充分供水”供水量；Wb,t为t月/旬“一级限供”供水量；Wloss,t为水库t月/旬蒸发、渗漏损失水量；Wp1,t、Wp2,t分别为一般枯水年和特枯水年t月/旬的水库来水量；Zhb,t、Zhb,t+1别为水库t、t+1月/旬的旱保水位；ZT+1为水文年末水位；Zdead为水库死水位或考虑水质要求允许取水的最低水位；Zlimit为水库汛限水位；Znormal为水库正常蓄水位；
在上一步推求的各月/旬旱限水位和旱警水位的基础上，基于步骤三确定的干旱预警分期，取各分期内最低月/旬旱限水位、旱警水位作为该分期旱限水位、旱警水位；
步骤六：分级分期旱限水位合理性分析与调整；
(1)基于步骤四确定的旱限水位供水规则以及步骤五确定的分级分期旱限水位，对历史长系列径流过程进行模拟，统计各用户供水保证率Rel、最大供水破坏深度Vul，计算公式如下：
Vul＝max(Deft/Demt) (9)
式中，N为总的统计时段数；Zt为t时段供水状态，若当前时段满足供水需求，则Zt＝1，若当前时段未满足供水需求，即发生供水破坏，Zt＝0；Deft为时段t的缺水量；Demt为时段t的需水量；
(2)判别供水保证率、破坏深度是否符合规划设计要求：若各用水户供水保证率、供水破坏深度指标均满足规划设计要求，则说明初步确定的旱限水位科学合理；若有任意指标不满足规划设计要求，则需要对初步确定的旱限水位进行反馈修正与调整，直至指标达标；
(3)分期旱限水位反馈修正的方法：
若供水保证率不达标，首先采用更换典型枯水年方式进行反馈修正，再按照“据旱限水位年内变化规律，调整对系统供水保证率、破坏深度影响大的旱限水位”思路进行反馈修正；
若供水保证率达标而供水破坏深度不达标，则采用提升高用水时期旱警水位、降低枯水期旱保水位的反馈修正方法；
若无法通过修正使得各用水户供水保证率、供水破坏深度指标达标，则说明水库供水系统无法承担当前设计需水量，应对设计需水量进行调整，直至各用水户供水保证率、供水破坏深度达到设计要求。
2.根据权利要求1所述的一种水库分级分期旱限水位设计与控制方法，其特征在于，所述步骤一具体为：
(1)分析规划水平年水库生态环境保护和经济社会保障目标，包括供水对象、供水时段、供水范围及其供水保障优先级；
(2)计算河道内、河道外各用户规划水平年需水量和需水过程，并明确各用水户设计供水保证率、供水破坏深度要求；河道外需水包括工业、农业、生活经济社会发展需水与生态环境需水，总量为各用户需水量之和；河道内需水包括生态、航运、生产等，总量按一水多用的原则，取各用户需水量的外包值；河道外经济社会发展需水计算参考GBT 51051-2014《水资源规划规范》，生态环境需水计算参考SL 429-2008《水资源供需预测分析技术规范》；河道内生态环境需水计算参考SL/Z 712《河湖生态环境需水计算规范》，交通航运需水计算参考GB 50139《内河通航标准》，其他生产需水量根据生产过程对流量、流速、水位等要求综合确定。
3.根据权利要求1所述的一种水库分级分期旱限水位设计与控制方法，其特征在于，所述步骤三中，采用Fisher最优分割法、成因分析法、模糊分析法、分形分析法对水文年进行干旱预警分期。
4.根据权利要求1所述的一种水库分级分期旱限水位设计与控制方法，其特征在于，所述步骤五中，采用逆序递推法初步确定旱限水位。
5.根据权利要求1所述的一种水库分级分期旱限水位设计与控制方法，其特征在于，所述步骤六中，基于长系列径流过程模拟，统计分析供水保证率与破坏深度，对初步确定的旱限水位进行和理性分析与调整。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111230658.8A CN113947315B (zh) | 2021-10-22 | 2021-10-22 | 一种水库分级分期旱限水位设计与控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111230658.8A CN113947315B (zh) | 2021-10-22 | 2021-10-22 | 一种水库分级分期旱限水位设计与控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113947315A true CN113947315A (zh) | 2022-01-18 |
| CN113947315B CN113947315B (zh) | 2024-07-12 |
Family
ID=79332018
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111230658.8A Active CN113947315B (zh) | 2021-10-22 | 2021-10-22 | 一种水库分级分期旱限水位设计与控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113947315B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114783164A (zh) * | 2022-06-20 | 2022-07-22 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种全类型供水水源地旱情预警与测算更新方法 |
| CN114997683A (zh) * | 2022-06-15 | 2022-09-02 | 中国水利水电科学研究院 | 一种基于旱限水位的水库精细化调度方法 |
| CN115018335A (zh) * | 2022-06-15 | 2022-09-06 | 中国水利水电科学研究院 | 一种控制性供水湖泊旱警水位确定方法 |
| CN115099477A (zh) * | 2022-06-15 | 2022-09-23 | 中国水利水电科学研究院 | 一种水库旱限水位优化及抗旱调度方法 |
| CN115169825A (zh) * | 2022-06-15 | 2022-10-11 | 中国水利水电科学研究院 | 一种确定不同调节能力水库的旱警水位的方法 |
| CN116703134A (zh) * | 2023-08-10 | 2023-09-05 | 长江勘测规划设计研究有限责任公司 | 一种大型跨流域调水水库的多目标调度方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140350737A1 (en) * | 2013-05-24 | 2014-11-27 | Judd Goodman | Optimized hydromodification management with active stormwater controls |
| CN112288183A (zh) * | 2020-11-12 | 2021-01-29 | 华中科技大学 | 一种湖泊分级分期旱限水位确定方法及系统 |
| CN113177189A (zh) * | 2021-05-20 | 2021-07-27 | 中国水利水电科学研究院 | 一种梯级水库分级分期旱限库容的计算方法 |
-
2021
- 2021-10-22 CN CN202111230658.8A patent/CN113947315B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140350737A1 (en) * | 2013-05-24 | 2014-11-27 | Judd Goodman | Optimized hydromodification management with active stormwater controls |
| CN112288183A (zh) * | 2020-11-12 | 2021-01-29 | 华中科技大学 | 一种湖泊分级分期旱限水位确定方法及系统 |
| CN113177189A (zh) * | 2021-05-20 | 2021-07-27 | 中国水利水电科学研究院 | 一种梯级水库分级分期旱限库容的计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张礼兵;伍露露;金菊良;吴成国;周玉良;朱文礼;: "大型灌区骨干水库分期旱限水位研究", 水利学报, no. 06, 26 June 2018 (2018-06-26) * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114997683A (zh) * | 2022-06-15 | 2022-09-02 | 中国水利水电科学研究院 | 一种基于旱限水位的水库精细化调度方法 |
| CN115018335A (zh) * | 2022-06-15 | 2022-09-06 | 中国水利水电科学研究院 | 一种控制性供水湖泊旱警水位确定方法 |
| CN115099477A (zh) * | 2022-06-15 | 2022-09-23 | 中国水利水电科学研究院 | 一种水库旱限水位优化及抗旱调度方法 |
| CN115169825A (zh) * | 2022-06-15 | 2022-10-11 | 中国水利水电科学研究院 | 一种确定不同调节能力水库的旱警水位的方法 |
| CN114997683B (zh) * | 2022-06-15 | 2022-11-29 | 中国水利水电科学研究院 | 一种基于旱限水位的水库精细化调度方法 |
| CN115018335B (zh) * | 2022-06-15 | 2023-10-27 | 中国水利水电科学研究院 | 一种控制性供水湖泊旱警水位确定方法 |
| CN114783164A (zh) * | 2022-06-20 | 2022-07-22 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种全类型供水水源地旱情预警与测算更新方法 |
| CN116703134A (zh) * | 2023-08-10 | 2023-09-05 | 长江勘测规划设计研究有限责任公司 | 一种大型跨流域调水水库的多目标调度方法及系统 |
| CN116703134B (zh) * | 2023-08-10 | 2023-11-10 | 长江勘测规划设计研究有限责任公司 | 大型跨流域调水水库的多目标调度方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113947315B (zh) | 2024-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113947315A (zh) | 一种水库分级分期旱限水位设计与控制方法 | |
| WO2021073192A1 (zh) | 一种考虑预报误差降低水库洪水起调水位的预报调度方法 | |
| Xun et al. | Evaluation of ecological sustainability based on a revised three-dimensional ecological footprint model in Shandong Province, China | |
| CN113177189B (zh) | 一种梯级水库分级分期旱限库容的计算方法 | |
| CN105913146B (zh) | 南方湿润地区的水资源优化配置系统 | |
| CN104047258B (zh) | 一种面向环境保护的联合调度水库群的生态库容确定方法 | |
| Zhang et al. | Risk analysis of dynamic control of reservoir limited water level by considering flood forecast error | |
| Diao et al. | Risk analysis of flood control operation mode with forecast information based on a combination of risk sources | |
| CN108133331B (zh) | 一种黄河下游宽滩区滞洪沉沙功能和减灾效应的评价方法 | |
| CN108681844A (zh) | 一种调水工程洪水资源利用风险评价方法 | |
| CN112819284B (zh) | 以流域为基础的生态环境演变评价方法及装置 | |
| CN109558626A (zh) | 基于时变设计洪水的梯级水库汛期运行水位动态控制方法 | |
| Yang et al. | Method for quantitatively assessing the impact of an inter-basin water transfer project on ecological environment-power generation in a water supply region | |
| Pan et al. | Analysis of the variation characteristics of runoff and sediment in the Yellow River within 70 years | |
| CN113177727B (zh) | 一种河道分级分期旱限水位或流量的计算方法 | |
| CN111126693B (zh) | 基于上游水库运行对下游电站发电能力影响的调度方法 | |
| CN116362552A (zh) | 一种评估小型水库安全风险等级的方法 | |
| CN116703134A (zh) | 一种大型跨流域调水水库的多目标调度方法及系统 | |
| CN111415103A (zh) | 一种基于psr的流域水安全度量方法 | |
| CN116402282A (zh) | 一种基于水量平衡的水库水位调控方法及系统 | |
| CN110348692A (zh) | 一种大型混联水库群多目标蓄能调度图计算方法 | |
| Zhao et al. | Research on energy-environment-economy-ecology coupling development in the yellow river basin | |
| CN115358587A (zh) | 区域多部门协同的基础设施规划方法及系统 | |
| CN109118073B (zh) | 一种水库洪水资源利用风险决策方法 | |
| CN110162848B (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 |