CN113177727B - 一种河道分级分期旱限水位或流量的计算方法 - Google Patents

一种河道分级分期旱限水位或流量的计算方法 Download PDF

Info

Publication number
CN113177727B
CN113177727B CN202110551124.9A CN202110551124A CN113177727B CN 113177727 B CN113177727 B CN 113177727B CN 202110551124 A CN202110551124 A CN 202110551124A CN 113177727 B CN113177727 B CN 113177727B
Authority
CN
China
Prior art keywords
water
drought
flow
river
demand
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
CN202110551124.9A
Other languages
English (en)
Other versions
CN113177727A (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 CN202110551124.9A priority Critical patent/CN113177727B/zh
Publication of CN113177727A publication Critical patent/CN113177727A/zh
Application granted granted Critical
Publication of CN113177727B publication Critical patent/CN113177727B/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/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
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression 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
    • 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
    • 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

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Data Mining & Analysis (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Databases & Information Systems (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Primary Health Care (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种河道分级分期旱限水位或流量的计算方法，包括根据江河的多种功能和干旱期江河断面的保障目标，具体分析河道的供水保障目标；根据河道的供水保障目标，计算河道内需水量和河道外需水量；综合河道内、河道外保障目标，对旱限水位或流量进行分级分期控制，即分级设置旱限水位或流量，具体分为旱警水位或流量和旱保水位或流量；江河断面划分有多个干旱预警分期，获取每个预警分期的分期旱限水位或流量；通过重现期分析法分析旱限水位或流量的合理性。优点是：有效的解决了抗旱减灾指挥决策时机把握不准或应急响应过度的现象，建立河道水位等干旱特征指标体系，为抗旱指挥决策提供科学依据和技术支撑。

Description

一种河道分级分期旱限水位或流量的计算方法
技术领域
本发明涉及水资源调度技术领域，尤其涉及一种河道分级分期旱限水位或流量的计算方法。
背景技术
干旱发生最直接的原因是区域内工农业生产生活及生态环境用水需求得不到满足，缺水程度越重，干旱灾害危害就越大，抗旱减灾的关键就是通过抗旱应急调水、科学分配水量、强化节约用水等措施，最大程度地缓解旱区缺水矛盾，减轻旱灾损失，因此，及时掌握及合理判断区域内水资源布局和供需缺口、影响程度就成为抗旱指挥决策的基本前提。一般而言，区域内工农业生产、生活及生态环境用水需求相对明确，江河来水、水库蓄水等水源可供水量往往成为抗旱应急水量调度最为关注的指标之一。江河湖库水位下降或者水量减少到某一临界值，区域内将会发生缺水，持续下降到某种程度，将会对城乡居民生活、工农业生产、航运及生态环境等各种用水需求造成影响，这些指标都是做好抗旱应急工作必须研究的重要课题。
河流是所在流域内自然地理背景下的产物，是以不同形态和经过不同转化途径的降水为补给来源的。河流对于人类社会的发展具有重要意义，干旱则会对沿岸工农业生产、城市供水、居民生活、航运交通和自然生态的安全产生巨大影响。开展江河断面旱限水位(流量)确定方法及试点应用研究，能够为各级决策部门和社会大众及时提供干旱预警，对全国抗旱应急管理工作将是一个有力的推动，具有重要的现实意义。
目前，在防汛领域，我国江河断面已经建立了较为完善的警戒水位、保证水位等防汛指标体系，在防洪减灾中发挥了重要作用，但在抗旱领域，我国至今尚未建立可以用于抗旱应急管理工作的水位(流量)等干旱特征指标，水文报旱、旱情预警、抗旱会商、应急响应、水量调度等抗旱应急管理工作往往缺乏科学依据，存在抗旱减灾指挥决策时机把握不准或应急响应过度的现象，一定程度上影响了抗旱工作科学有序地开展。因此，迫切需要建立河道水位(流量)等干旱特征指标体系，为抗旱指挥决策提供科学依据和技术支撑。
发明内容
本发明的目的在于提供一种河道分级分期旱限水位或流量的计算方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明采用的技术方案如下：
一种河道分级分期旱限水位或流量的计算方法，包括如下步骤，
S1、根据江河的多种功能和干旱期江河断面的保障目标，具体分析河道的供水保障目标；
S2、根据河道的供水保障目标，计算河道内需水量和河道外需水量；
S3、综合河道内、河道外保障目标，由于年内不同时期江河断面水位或流量存在差异，应针对性设置不同的旱限水位或流量，即分期设置旱限水位或流量；在年内同一时期，江河断面水位或流量不同时，保障未来供水能力不同，应设置不同预警级别的旱限水位或流量，对旱限水位或流量进行分级分期控制，即分级设置旱限水位或流量，具体分为旱警水位或流量和旱保水位或流量；
S4、江河断面划分有多个干旱预警分期，在同一预警分期内存在多个不同的旱限水位或流量，取同一预警分期中所有旱限水位或流量中的最大值作为该预警分期的分期旱限水位或流量；
S5、通过重现期分析法分析旱限水位或流量的合理性，当重现期结果与历史上旱情出现的重现期吻合，则表示确定的旱限水位或流量合理。
优选的，河道的供水保障目标包括社会经济供水保障目标、生态环境保障目标和航运保障目标；步骤S1具体包括如下内容，
S11、通过调查收集江河取水口及供水信息，分析确定旱限期需要重点保障的社会经济供水保障目标；所述江河取水口及供水信息包括取水口位置、取水口高程以及各取水口的受水区范围、用水行业、用水时段和用水量；
S12、通过调查收集江河生态环境信息，分析确定旱限期需要重点保障的相关生态环境保障目标；所述江河生态环境信息包括水生栖息及洄游物种的生态需水量、冲沙和污染物输移需水量以及环境容量；
S13、通过调查收集江河航运信息，分析确定旱限期需要重点保障的航运保障目标；所述江河航运信息包括通航水域、航道等级、通航时间、通航水位和航运调度。
优选的，所述河道外需水量为社会经济需水量；所述河道内需水量包括生态需水量和航运需水量；所述生态需水量包括基本生态环境需水量和适宜生态环境需水量；步骤S2具体包括如下内容，
S21、利用用水量调查统计法或定额计算法计算社会经济需水量；
S22、以当地供水规则为准，对比采用7Q10法、最小连续30d平均流量法、流量百分数法和Tennant法的计算结果，或采用近10年最枯月平均流量法确定基本生态环境需水量；以当地供水规则为准利用Tennant法确定适宜生态环境需水量；
S23、利用保障率频率法计算航运需水量。
优选的，利用用水量调查统计法计算社会经济需水量的具体过程为，
调查统计断面下游10年分行业用水统计数据；
依据用水类型及方式，重点调查统计城乡供水、企业生产、农业灌溉及环境生态用水数据；
选取来水频率为75％的多个年份作为一般枯水年组，对各个行业，计算一般枯水年组中逐月平均用水量，作为该行业在一般干旱年份社会经济需水量；同时，选取来水频率为95％的多个年份作特枯水年组，对各个行业，计算特枯水年组中逐月平均用水量，作为该行业在特大干旱年份社会经济需水量。
优选的，利用定额计算法计算社会经济需水量的具体过程为，
根据规划水平年居民生活用水定额和人口发展预测成果计算生活需水量；
根据规划水平年工业用水定额及工业生产总值计算工业需水量；
根据规划水平年灌溉面积、种植结构、灌溉制度和灌溉利用系数计算灌溉需水量；
综合生活需水量、工业需水量和灌溉需水量，获取社会经济需水量。
优选的，在计算基本生态环境需水量时，对于长系列水文资料河流控制断面，采用流量历时曲线法、7Q10法或Qp法计算基本生态环境需水量，通过对比计算结果合理确定基本生态环境需水量；对于缺乏长系列水文资料的河流控制断面，采用近10年最枯月平均流量法计算基本生态环境需水量。
优选的，利用Tennant法计算适宜生态环境需水量具体为，
Figure BDA0003075371040000041
其中，WR为多年平均条件下维持河道一定功能的需水量；Mi为第i月天数；Qi为i月多年平均流量；Pi为第i月适宜生态环境需水百分比，少水期通常选取多年平均流量的10％～20％作为河道适宜生态环境需水量，多水期选取多年平均流量的30％～40％作为河道适宜生态环境需水量。
优选的，步骤S3中利用外包线法以及逐级控制、逐级叠加旱限水位或流量法，计算各分期各级别的江河旱限水位或流量；具体过程为，
S31、利用外包线法取生态需水量和航运需水量各分期的最大值；
Wrs,i＝max{Wr，i，Ws，i}
其中，Wr，i为第i月的河道生态需水量；Ws,i为第i月的河道航运需水量；Wrs,i为生态需水量和航运需水量第i月的最大值；
S32、利用获取的生态需水量和航运需水量各分期的最大值，使用逐级控制、逐级叠加旱限水位或流量法，获取各分期各级别的江河旱限水位或流量；
Qi＝μ{max{Wr,i,Ws,i}+(Wul,i+WF，i+WLI，i)}
其中，Wul,,i、WF,i和WLI,i分别为第i月一般枯水年或特枯年的月生活需水量、工业需水量和灌溉需水量；Qi为第i月旱警或旱保流量；μ()为河道需水量-流量转换函数；
Zi＝max{H，z(Qi)}
其中，Zi为月旱警或旱保水位；z()为河道流量-水位转换函数；H为水口高程。
优选的，步骤S4中旱警/旱保水位或流量分期的计算方式分别如下，
QT＝max{Q1，Q2…Qi，}
ZT＝max{Z1，Z2…Zi，}
其中，QT为分期T的旱警或旱保流量；ZT为分期T的旱警或旱保水位；Qi为分期T内第i个月的河道旱警或旱保流量；Zi为分期T内第i个月的河道旱警或旱保水位。
优选的，步骤S5具体包括如下内容，
S51、明确旱限流量分期(T1,T2…)，获得各分期的旱限流量；
S52、从长系列数据第一年开始遍历所有年份分期内的月流量，统计分期内流量低于旱限流量的年数n；
S53、获得分期实际流量低于河道旱限流量的重现期R；
Figure BDA0003075371040000051
其中，N为长系列流量资料的年数；旱限水位重现期与旱限流量重现期计算方法相同；
S54、判断获取的重现期R与历史上旱情出现的重现期是否吻合，若是，则表示确定的旱限水位或流量合理；否则，则表示确定的旱限水位或流量不合理。
本发明的有益效果是：综合考虑江河断面的社会经济供水保障目标、生态环境保障目标和航运保障目标等主要用水需求，结合河道一定的设计来水过程，通过逐级控制、逐级叠加旱限水位(流量)法及外包线法，对旱限水位(流量)进行分期，建立可以用于抗旱应急管理工作的水位等干旱特征指标，有效的解决了抗旱减灾指挥决策时机把握不准或应急响应过度的现象，建立河道水位等干旱特征指标体系，为抗旱指挥决策提供科学依据和技术支撑。
附图说明
图1是本发明实施例中计算方法的流程示意图；
图2是本发明实施例中水文站A干旱预警分期示意图；
图3是本发明实施例中外包线、逐级叠加和逐级控制的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例一
如图所示，本实施例中，提供了一种河道分级分期旱限水位或流量的计算方法，包括如下步骤，
S1、根据江河的多种功能和干旱期江河断面的保障目标，具体分析河道的供水保障目标；
S2、根据河道的供水保障目标，计算河道内需水量和河道外需水量；
S3、综合河道内、河道外保障目标，由于年内不同时期江河断面水位(流量)存在差异，应针对性设置不同的旱限水位(流量)，即分期设置旱限水位(流量)；在年内同一时期，江河断面水位(流量)不同时，保障未来供水能力不同，应设置不同预警级别的旱限水位(流量)，对旱限水位(流量)进行分级分期控制，即分级设置旱限水位(流量)，具体分为旱警水位(流量)和旱保水位(流量)；
S4、江河断面划分有多个干旱预警分期，在同一预警分期内存在多个不同的旱限水位(流量)，取同一预警分期中所有旱限水位(流量)中的最大值作为该预警分期的分期旱限水位(流量)；
S5、通过重现期分析法分析旱限水位或流量的合理性，当重现期结果与历史上旱情出现的重现期吻合，则表示确定的旱限水位或流量合理。
本实施例中，所述计算方法具体包括五部分内容，分别为分析供水保障目标、计算河道内及河道外用水需水量、计算各分期各级别的旱限水位(流量)、计算各个预警分期的旱限水位(流量)、判断计算获取的旱限水位(流量)的合理性。下面分别针对这五部分内容进行详细的解释说明。
一、分析供水保障目标
该部分对应上述步骤S1的内容，河道的供水保障目标包括社会经济供水保障目标、生态环境保障目标和航运保障目标；步骤S1具体包括如下内容，
S11、社会经济供水保障目标；通过调查收集江河取水口及供水信息，分析确定旱限期需要重点保障的社会经济供水保障目标；所述江河取水口及供水信息包括取水口位置、取水口高程以及各取水口的受水区范围、用水行业、用水时段和用水量等资料；
S12、生态环境保障目标；通过调查收集江河生态环境信息，分析确定旱限期需要重点保障的相关生态环境保障目标；所述江河生态环境信息包括水生栖息及洄游物种的生态需水量、冲沙和污染物输移需水量以及环境容量等资料；
S13、航运保障目标；通过调查收集江河航运信息，分析确定旱限期需要重点保障的航运保障目标；所述江河航运信息包括通航水域、航道等级、通航时间、通航水位和航运调度等资料。
二、计算河道内及河道外用水需求量
该部分对应步骤S2的内容，所述河道外需水量为社会经济需水量；所述河道内需水量包括生态需水量和航运需水量；所述生态需水量包括基本生态环境需水量和适宜生态环境需水量；步骤S2具体包括如下内容，
S21、利用用水量调查统计法或定额计算法计算社会经济需水量；
S22、以当地供水规则为准，对比采用7Q10法、最小连续30d平均流量法、流量百分数法和Tennant法的计算结果，或采用近10年最枯月平均流量法确定基本生态环境需水量；以当地供水规则为准利用Tennant法确定适宜生态环境需水量；
S23、利用保障率频率法计算航运需水量。
本实施例中，计算社会经济需水量包括两种方式，分别为：
1、利用用水量调查统计法计算社会经济需水量的具体过程为，
调查统计断面下游近10年分行业用水统计数据；
依据用水类型及方式，重点调查统计城乡供水、企业生产、农业灌溉及环境生态用水数据；
选取来水频率为75％的多个年份作为一般枯水年组，对各个行业，计算一般枯水年组中逐月平均用水量，作为该行业在一般干旱年份社会经济需水量；同时，选取来水频率为95％的多个年份作特枯水年组，对各个行业，计算特枯水年组中逐月平均用水量，作为该行业在特大干旱年份社会经济需水量。
2、利用定额计算法计算社会经济需水量的具体过程为，(采用定额计算法时，社会经济需水量计算参考GBT 51051-2014《水资源规划规范》)
根据规划水平年居民生活用水定额和人口发展预测成果计算生活需水量；
根据规划水平年工业用水定额及工业生产总值计算工业需水量；
根据规划水平年灌溉面积、种植结构、灌溉制度和灌溉利用系数计算灌溉需水量；
综合生活需水量、工业需水量和灌溉需水量，获取社会经济需水量。
生态需水量尽量以当地供水规则为准，参考《河湖生态环境需水计算规范(SL/Z712-2014)》，所述生态需水量包括基本生态环境需水量和适宜生态环境需水量。
本实施例中，在计算基本生态环境需水量时，可在7Q10法、最小连续30d平均流量法、流量百分数法、Tennant法等基础上，比较多种方法计算结果，合理确定基本生态环境需水量。应按下列规定执行：
1、对于长系列(n>30年)水文资料河流控制断面，采用流量历时曲线法、7Q10法或Qp法计算基本生态环境需水量；
2、对于缺乏长系列水文资料的河流控制断面，采用近10年最枯月平均流量法计算基本生态环境需水量。
3、应比较分析多种方法计算结果，合理确定基本生态环境需水量最小值。
本实施例中，适宜生态环境需水量可按照保护目标对应的生态环境功能维持在正常水平的需水量要求采取Tennant法进行计算；
利用Tennant法计算适宜生态环境需水量具体为，
Figure BDA0003075371040000081
其中，WR为多年平均条件下维持河道一定功能的需水量；Mi为第i月天数；Qi为i月多年平均流量；Pi为第i月适宜生态环境需水百分比，少水期通常选取多年平均流量的10％～20％作为河道适宜生态环境需水量，多水期选取多年平均流量的30％～40％作为河道适宜生态环境需水量，如表1。
表1在河道内不同生境情况下推荐的河道流量百分比
Figure BDA0003075371040000082
本实施例中，航运需水量的计算参考《内河通航标准》(GB 40139-2014)，在不受潮汐影响和湖汐影响不明显的河段，采用保证率频率法计算确定通航水位。
三、计算各分期各级别的江河旱限水位(流量)
该部分对应步骤S3的内容，步骤S3中利用外包线法以及逐级控制、逐级叠加旱限水位(流量)法，计算各分期各级别的江河旱限水位(流量)；具体过程为，
S31、利用外包线法取生态需水量和航运需水量各分期(各月)的最大值；
Wrs，i＝max{Wr，i，Ws，i}
其中，Wr，i为第i月的河道生态需水量；Ws，i为第i月的河道航运需水量；Wrs，i为生态需水量和航运需水量第i月的最大值；
S32、利用获取的生态需水量和航运需水量各分期的最大值，使用逐级控制、逐级叠加旱限水位(流量)法，获取各分期各级别的江河旱限水位(流量)；
Qi＝μ{max{Wr,i,Ws,i}+(Wul，i+WF,i+WLI,i)}
其中，Wul,,i、WF,i和WLI,i分别为第i月一般枯水年或特枯年的月生活需水量、工业需水量和灌溉需水量；Qi为第i月旱警(旱保)流量；μ()为河道需水量-流量转换函数；
Zi＝max{H，z(Qi)}
其中，Zi为月旱警(旱保)水位；z()为河道流量-水位转换函数；H为水口高程。
四、计算各个预警分期的旱限水位(流量)
该部分对应步骤S4的内容，步骤S4中旱警(旱保)的水位(流量)分期的计算方式分别如下，
QT＝max{Q1,Q2…Qi,}
ZT＝max{Z1，Z2…Zi，}
其中，QT为分期T的旱警(旱保)流量；ZT为分期T的旱警(旱保)水位；Qi为分期T内第i个月的河道旱警(旱保)流量；Zi为分期T内第i个月的河道旱警(旱保)水位。
江河断面一般会划分为几个干旱预警分期，在同一预警分期内可能会有几个不同的旱限水位(流量)，此时取这些旱限水位中的最大值作为该分期的旱限水位(流量)。
五、判断计算获取的旱限水位(流量)的合理性
该部分对应步骤S5的内容，步骤S5具体包括如下内容，
S51、明确旱限水位(流量)分期(T1,T2…)，获得各分期的旱限水位(流量)；
S52、从长系列数据第一年开始遍历所有年份分期内的月流量，统计分期内流量低于旱限水位(流量)的年数n；
S53、获得分期实际流量低于河道旱限水位(流量)的重现期R；
Figure BDA0003075371040000101
其中，N为长系列流量资料的年数；旱限水位重现期与旱限流量重现期计算方法相同；
S54、判断获取的重现期R与历史上旱情出现的重现期是否吻合，若是，则表示确定的旱限水位(流量)合理；否则，则表示确定的旱限水位(流量)不合理。
通过重现期分析法以及历史干旱一致性分析法判断旱限水位(流量)的合理性，其中重现期分析法建议南方地区旱警水位(流量)10年一遇，旱保水位(流量)20年一遇；北方地区旱警水位(流量)5年一遇，旱保水位(流量)10年一遇。
实施例二
本实施例中，具体举例说明本发明的计算方法的实施过程。如图1所示，某水文站A下游沿河取水用户为一大型灌区，取水口位于河床底部，对断面水位没有要求。
结合河道特点，参考步骤S1的内容，本实施例中计算主要考虑河道的社会经济需水量、生态需水量和航运需水量。
结合步骤S2的内容，本实施例中，社会经济需水量的计算分别采用用水量调查统计法和定额法计算法得到灌区一般枯水年和特枯水年的逐月需水过程。从偏安全的角度考虑，对两种方法得到的用水量逐月取最大值得到综合需水量，作为旱限水位(流量)计算的基础。这两种方法计算获取的各月份的灌区抽水量如下表2。
表2灌区抽水量(m3/s)
Figure BDA0003075371040000111
生态需水量中以Tennant法计算河道适宜生态流量，20.1m3/s；根据水文站所在河流水量调度细则，该水文站需要保障低限流量12m3/s，预警流量25m3/s。为了与现行调度管理工作相结合，采用12m3/s作为最小生态流量(基本生态环境需水量)，25m3/s作为适宜生态流量。获取的各鱼粉的河道生态需水量具体如下表3所示。
表3河道生态需水量(m3/s)
Figure BDA0003075371040000112
Figure BDA0003075371040000121
航运需水量的计算需要参考《内河通航标准》(GB 40139-2014)，不受潮汐影响和湖汐影响不明显的河段，采用保证率频率法计算确定通航水位。获取的各月份的河道航运需水量如下表4所示。
表4河道航运需水量(m3/s)
月份 75％排序年灌区需水量 95％排序年灌区需水量
1月 6 7
2月 3 4
3月 8 9
4月 10 10
5月 9 10
6月 10 11
7月 19 21
8月 20 22
9月 0 0
10月 0 0
11月 3 4
12月 6 7
根据步骤S3，采用逐月累计的方法计算逐月旱警水位和逐月旱保水位，计算结果如下表5所示。
表5水文站A逐月旱警、旱保流量(m3/s)
Figure BDA0003075371040000122
Figure BDA0003075371040000131
根据步骤S4，综合考虑水文站A多年月平均降水量、径流量、水位，以及灌区用水量规律，对该水文站进行干旱预警分期，将一个水文年分为4个时期，3-5月为农业灌溉期，6-8月为主汛期，9-10月为汛末，11～次年2月为非汛期。采用分期内各月旱警(保)流量最大值作为分期旱警流量和分期旱保流量，结果如下表6所示。
表6水文站A分期旱警、旱保流量(m3/s)
Figure BDA0003075371040000132
根据步骤S5，采用重现期法进行旱限水位(流量)合理性分析。根据1961～2018年58年逐月流量资料统计，旱警流量中枯水期共有14年流量低于旱警流量，灌溉期有9年流量低于旱警流量，汛期有9年流量低于旱警流量，重现期在4-6年以内；旱保流量中枯水期共有7年流量低于旱保流量，灌溉期有6年流量低于旱保流量，汛期有6年流量低于旱保流量，重现期在8-10年以内。合理性验证的结果如下表7所示。
表7水文站A旱限流量合理性检验
Figure BDA0003075371040000133
Figure BDA0003075371040000141
通过采用本发明公开的上述技术方案，得到了如下有益的效果：
本发明提供了一种河道分级分期旱限水位或流量的计算方法，综合考虑江河断面的社会经济供水保障目标、生态环境保障目标和航运保障目标等主要用水需求，结合河道一定的设计来水过程，通过逐级控制、逐级叠加旱限水位(流量)及外包线法，对旱限水位(流量)进行分期，建立可以用于抗旱应急管理工作的水位等干旱特征指标，有效的解决了抗旱减灾指挥决策时机把握不准或应急响应过度的现象，建立河道水位等干旱特征指标体系，为抗旱指挥决策提供科学依据和技术支撑。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。

Claims (7)

1.一种河道分级分期旱限水位或流量的计算方法，其特征在于：包括如下步骤，
S1、根据江河的多种功能和干旱期江河断面的保障目标，具体分析河道的供水保障目标；
S2、根据河道的供水保障目标，计算河道内需水量和河道外需水量；
S3、综合河道内、河道外保障目标，由于年内不同时期江河断面水位或流量存在差异，应针对性设置不同的旱限水位或流量，即分期设置旱限水位或流量；在年内同一时期，江河断面水位或流量不同时，保障未来供水能力不同，应设置不同预警级别的旱限水位或流量，对旱限水位或流量进行分级分期控制，即分级设置旱限水位或流量，具体分为旱警水位或流量和旱保水位或流量；
S4、江河断面划分有多个干旱预警分期，在同一预警分期内存在多个不同的旱限水位或流量，取同一预警分期中所有旱限水位或流量中的最大值作为该预警分期的分期旱限水位或流量；
S5、通过重现期分析法分析旱限水位或流量的合理性，当重现期结果与历史上旱情出现的重现期吻合，则表示确定的旱限水位或流量合理；
所述河道外需水量为社会经济需水量；所述河道内需水量包括生态需水量和航运需水量；所述生态需水量包括基本生态环境需水量和适宜生态环境需水量；
步骤S2具体包括如下内容，
S21、利用用水量调查统计法或定额计算法计算社会经济需水量；
S22、以当地供水规则为准，对比采用7Q10法、最小连续30d平均流量法、流量百分数法和Tennant法的计算结果，或采用近10年最枯月平均流量法确定基本生态环境需水量；以当地供水规则为准利用Tennant法确定适宜生态环境需水量；
S23、利用保障率频率法计算航运需水量；
步骤S3中利用外包线法以及逐级控制、逐级叠加旱限水位或流量法，计算各分期各级别的江河旱限水位或流量；具体过程为，
S31、利用外包线法取生态需水量和航运需水量各分期的最大值；
Wrs,i＝max{Wr,i,Ws,i}
其中，Wr,i为第i月的河道生态需水量；Ws,i为第i月的河道航运需水量；Wrs,i为生态需水量和航运需水量第i月的最大值；
S32、利用获取的生态需水量和航运需水量各分期的最大值，使用逐级控制、逐级叠加旱限水位或流量法，获取各分期各级别的江河旱限水位或流量；
Qi＝μ{max{Wr,i,Ws,i}+(Wul,i+WF,i+WLI,i)}
其中，Wul,i、WF,i和WLI,i分别为第i月一般枯水年或特枯年的月生活需水量、工业需水量和灌溉需水量；Qi为第i月旱警或旱保流量；μ()为河道需水量-流量转换函数；
Zi＝max{H,z(Qi)}
其中，Zi为月旱警或旱保水位；z()为河道流量-水位转换函数；H为水口高程；
步骤S5具体包括如下内容，
S51、明确旱限流量分期(T1,T2…)，获得各分期的旱限流量；
S52、从长系列数据第一年开始遍历所有年份分期内的月流量，统计分期内流量低于旱限流量的年数n；
S53、获得分期实际流量低于河道旱限流量的重现期R；
Figure FDA0003354420210000021
其中，N为长系列流量资料的年数；旱限水位重现期与旱限流量重现期计算方法相同；
S54、判断获取的重现期R与历史上旱情出现的重现期是否吻合，若是，则表示确定的旱限水位或流量合理；否则，则表示确定的旱限水位或流量不合理。
2.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：河道的供水保障目标包括社会经济供水保障目标、生态环境保障目标和航运保障目标；步骤S1具体包括如下内容，
S11、通过调查收集江河取水口及供水信息，分析确定旱限期需要重点保障的社会经济供水保障目标；所述江河取水口及供水信息包括取水口位置、取水口高程以及各取水口的受水区范围、用水行业、用水时段和用水量；
S12、通过调查收集江河生态环境信息，分析确定旱限期需要重点保障的相关生态环境保障目标；所述江河生态环境信息包括水生栖息及洄游物种的生态需水量、冲沙和污染物输移需水量以及环境容量；
S13、通过调查收集江河航运信息，分析确定旱限期需要重点保障的航运保障目标；所述江河航运信息包括通航水域、航道等级、通航时间、通航水位和航运调度。
3.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：利用用水量调查统计法计算社会经济需水量的具体过程为，
调查统计断面下游10年分行业用水统计数据；
依据用水类型及方式，重点调查统计城乡供水、企业生产、农业灌溉及环境生态用水数据；
选取来水频率为75％的多个年份作为一般枯水年组，对各个行业，计算一般枯水年组中逐月平均用水量，作为该行业在一般干旱年份社会经济需水量；同时，选取来水频率为95％的多个年份作特枯水年组，对各个行业，计算特枯水年组中逐月平均用水量，作为该行业在特大干旱年份社会经济需水量。
4.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：利用定额计算法计算社会经济需水量的具体过程为，
根据规划水平年居民生活用水定额和人口发展预测成果计算生活需水量；
根据规划水平年工业用水定额及工业生产总值计算工业需水量；
根据规划水平年灌溉面积、种植结构、灌溉制度和灌溉利用系数计算灌溉需水量；
综合生活需水量、工业需水量和灌溉需水量，获取社会经济需水量。
5.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：在计算基本生态环境需水量时，对于长系列水文资料河流控制断面，采用流量历时曲线法、7Q10法或Qp法计算基本生态环境需水量，通过对比计算结果合理确定基本生态环境需水量；对于缺乏长系列水文资料的河流控制断面，采用近10年最枯月平均流量法计算基本生态环境需水量。
6.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：利用Tennant法计算适宜生态环境需水量具体为，
Figure FDA0003354420210000041
其中，WR为多年平均条件下维持河道一定功能的需水量；Mi为第i月天数；Qi为i月多年平均流量；Pi为第i月适宜生态环境需水百分比，少水期通常选取多年平均流量的10％～20％作为河道适宜生态环境需水量，多水期选取多年平均流量的30％～40％作为河道适宜生态环境需水量。
7.根据权利要求1所述的河道分级分期旱限水位或流量的计算方法，其特征在于：步骤S4中旱警/旱保水位或流量分期的计算方式分别如下，
QT＝max{Q1,Q2...Qi,}
ZT＝max{Z1,Z2...Zi,}
其中，QT为分期T的旱警或旱保流量；ZT为分期T的旱警或旱保水位；Qi为分期T内第i个月的河道旱警或旱保流量；Zi为分期T内第i个月的河道旱警或旱保水位。
CN202110551124.9A 2021-05-20 2021-05-20 一种河道分级分期旱限水位或流量的计算方法 Active CN113177727B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110551124.9A CN113177727B (zh) 2021-05-20 2021-05-20 一种河道分级分期旱限水位或流量的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110551124.9A CN113177727B (zh) 2021-05-20 2021-05-20 一种河道分级分期旱限水位或流量的计算方法

Publications (2)

Publication Number Publication Date
CN113177727A CN113177727A (zh) 2021-07-27
CN113177727B true CN113177727B (zh) 2022-02-01

Family

ID=76929811

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110551124.9A Active CN113177727B (zh) 2021-05-20 2021-05-20 一种河道分级分期旱限水位或流量的计算方法

Country Status (1)

Country Link
CN (1) CN113177727B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115018335B (zh) * 2022-06-15 2023-10-27 中国水利水电科学研究院 一种控制性供水湖泊旱警水位确定方法
CN116222705B (zh) * 2023-03-16 2023-11-24 中国水利水电科学研究院 一种地下水旱限水位确定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104047259A (zh) * 2014-06-23 2014-09-17 清华大学 一种面向生态环境保护的待建水库的生态库容确定方法
CN112288183A (zh) * 2020-11-12 2021-01-29 华中科技大学 一种湖泊分级分期旱限水位确定方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140343855A1 (en) * 2013-05-15 2014-11-20 The Regents Of The University Of California Drought Monitoring and Prediction Tools

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104047259A (zh) * 2014-06-23 2014-09-17 清华大学 一种面向生态环境保护的待建水库的生态库容确定方法
CN112288183A (zh) * 2020-11-12 2021-01-29 华中科技大学 一种湖泊分级分期旱限水位确定方法及系统

Also Published As

Publication number Publication date
CN113177727A (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
WO2021109848A1 (zh) 一种考虑提升量的生态流量确定方法
Fattahi et al. A compromise programming model to integrated urban water management
CN113177727B (zh) 一种河道分级分期旱限水位或流量的计算方法
CN113177189B (zh) 一种梯级水库分级分期旱限库容的计算方法
Li et al. Multistage integrated water security assessment in a typical region of Northwestern China
CN111401732B (zh) 一种水资源供需形势分层次识别方法
Sivakumar Strategy to be adopted in Preparation of National Water Resource Master Plan
O’riordan et al. Choice in water use
Tursunova et al. Water resources of Kazakhstan in conditions of uncertainty
Ho et al. Apply MIKE 11 model to study impacts of climate change on water resources and develop adaptation plan in the Mekong Delta, Vietnam: a case of Can Tho city
Fellizar Urban water management in metropolitan Manila, Philippines
Ivanova et al. Optimization model of a tandem water reservoir system management
CN114897396A (zh) 评价河湖水系连通水安全保障能力的方法
Zhou Vulnerability and adaptation to climate change in North China: the water sector in Tianjin
Tarjuelo et al. Water resources deficit and water engineering
LU502097B1 (en) Calculation Method for Grading and Staged Drought-limited Water Level or Flow Rate in River Course
Statzu et al. Water demand for residential uses in a mediterranean region: Econometric analysis and policy implications
Kampragou et al. Towards policy recommendations for future drought risk reduction
Dinar et al. The Evolving Nature of California’s Water Economy
DeMaagd Water Demand in the Residential and Tourism Sectors: Evidence and Implications for Efficient Management
Baran Integrated Model-Based Impact Assessment of Climate Change and Land Use Change on the Occoquan Watershed
Rhode Using the standardised precipitation index to determine the severity, duration, and frequency of drought events in the western cape water supply system area
Ait-Aoudia et al. Algiers water issue: Key indicators for a diagnosis of water resources carrying capacity
CN116090846A (zh) 基于污染物允许排放量的水环境管控分区划分方法
Shahraki Sustainable development in hydro-drought regions by improving hydro-indicators

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