CN113256028B - 一种面向大范围长历时干旱的应急供水调配方法 - Google Patents
一种面向大范围长历时干旱的应急供水调配方法 Download PDFInfo
- Publication number
- CN113256028B CN113256028B CN202110695430.XA CN202110695430A CN113256028B CN 113256028 B CN113256028 B CN 113256028B CN 202110695430 A CN202110695430 A CN 202110695430A CN 113256028 B CN113256028 B CN 113256028B
- Authority
- CN
- China
- Prior art keywords
- water
- drought
- water supply
- emergency
- allocation
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 561
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000003860 storage Methods 0.000 claims abstract description 32
- 238000011156 evaluation Methods 0.000 claims abstract description 9
- 238000002156 mixing Methods 0.000 claims description 54
- 238000005457 optimization Methods 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012546 transfer Methods 0.000 claims description 15
- 238000012384 transportation and delivery Methods 0.000 claims description 13
- 239000002352 surface water Substances 0.000 claims description 11
- 238000010586 diagram Methods 0.000 claims description 8
- 238000004519 manufacturing process Methods 0.000 claims description 7
- 238000005065 mining Methods 0.000 claims description 6
- 230000001105 regulatory effect Effects 0.000 claims description 6
- 238000011084 recovery Methods 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 3
- 238000005086 pumping Methods 0.000 claims description 2
- 238000001556 precipitation Methods 0.000 description 6
- 238000005192 partition Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000001704 evaporation Methods 0.000 description 3
- 230000008020 evaporation Effects 0.000 description 3
- 230000008595 infiltration Effects 0.000 description 3
- 238000001764 infiltration Methods 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000013468 resource allocation Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000009189 diving Effects 0.000 description 2
- 239000008235 industrial water Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000013316 zoning Methods 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/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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种面向大范围长历时干旱的应急供水调配方法,包括对干旱区内各调配单元常规水源和应急水源分类识别;针对不同干旱程度选取应急供水调配模式;计算分析各调配单元不同干旱程度下供需水量;制定各调配单元不同调配模式下各水源和用水户干旱应急供水调配规则;以整体缺水量最小、整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为目标构建双层大范围长历时干旱应急供水调配模型,所述模型中包括数据输入模块、干旱程度评价模块、供需水模块和优化调配模块。本发明能提高干旱条件下应急供水调配效率、公平性及可持续性。
Description
技术领域
本发明涉及水资源调配的技术领域,特别是一种面向大范围长历时干旱的应急供水调配方法。
背景技术
干旱可能发展为重大的环境灾难,会导致重大的经济损失、社会问题和生态影响。随着全球气候变化加剧,极端气候事件发生的频率不断增加,旱灾发生的频率及规模也不断加大,不仅影响粮食安全,还会引发供水危机和生态危机。干旱是最主要的自然灾害之一。水资源调配是缓解旱情的有效手段之一,干旱发生时,通过水利工程群的调配,不仅能在一定程度上增加整个干旱期的可供水量,还能灵活地控制供水过程,适当压缩干旱初期供水量,保证干旱发展后期基本的用水需要,从而降低重度干旱在整个干旱期的比重,减少旱灾损失。目前,已有众多学者进行相关研究,例如利用多目标线性规划优化干旱期间城市供水;构建双层动态优化模型解决多水库之间干旱应急调度问题等。以上研究大部分集中于封闭区域内部各用水户之间的水量分配,另有一部分着眼于小范围内单个水库干旱应急调水决策问题,均未在大范围长历时干旱情况下的配水和调水问题进行统一考虑。
2014年2月的《华北水利水电大学学报》的第35卷第1期公开了邵东国、李旭东、唐明和黄文红的《干旱条件下城市水资源应急调配模型》的文章,该文章以南昌市为例,构建了基于城市供水水源管理与行政分区的二层大系统分解协调应急调配模型,分别以2010年、2020年、2030年为水平年,在不同干旱条件下进行水资源供需分析,结果表明,该市南昌县与进贤县7-9月最易发生干旱缺水;不同干旱条件下,全市缺水率将达23.70%-45.85%;如采用水资源统一调配供水方案,则可使全市缺水率下降到0.00%-41.63%,基于水资源统一管理理念及实际情况,限制性协调供水方案可用来降低缺水率,减轻干旱带来的损失,论证了所提模型方法的有效性。该方法的缺点是:①该方法只考虑不同水平年、不同来水频率下的干旱情景,没有考虑到干旱在年内的动态演变过程;②没有对区域内的常规水源及应急水源进行划分,不便于指导决策者对是否启用应急水源等相关措施进行决策;③在调配情景中只考虑了供水端的不同方案,没有考虑不同的干旱阶段在供需两侧应采取怎样的抗旱措施。
发明内容
为了解决上述的技术问题,本发明提出的一种面向大范围长历时干旱的应急供水调配方法,针对大范围长历时干旱特点,行政分区及水资源分区对干旱区进行划分,对各调配单元常规水源和应急水源进行分类识别,并对不同干旱程度下各水源可供水量及各行业需水量进行评估计算,在基于时段旱情判断的基础上,制定调配模式与调配规则,以整体缺水量最小、整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为目标,构建应急供水调配模型,求解得到调配方案。
本发明的目的是提供一种面向大范围长历时干旱的应急供水调配方法,包括以下步骤:
步骤1:根据行政分区和水资源分区对大范围干旱区进行划分,得到调配单元;
步骤2:对所述调配单元中的常规水源和应急水源进行分类统计;
步骤3:对干旱区利用ArcGIS软件及DEM数据对不同干旱程度下各水库之间水力联系进行识别,确定上下游关系,并对水库进行编码;
步骤4:建立各水源、调配单元之间的拓扑关系,绘制干旱区水资源系统网络图;
步骤5:针对不同干旱程度选定相应的应急供水调配模式;
步骤6:对各调配单元不同干旱程度下可供水量及需水量进行分析计算;
步骤7:制定不同调配模式下各水源供水及各行业用水优先序;
步骤8:构建双层大范围长历时干旱应急供水调配模型,优化求解得到调配方案。
优选的是,所述步骤2包括对两种水源中地表水、地下水、非常规水、外调水数量和基本参数进行分类统计,用于对不同干旱程度不同调配模式下可供水量进行计算。
在上述任一方案中优选的是,所述分类统计的结果还用于在不同干旱程度下,提取相应水源的可供水量进行调配,包括:
1)发生轻度干旱时使用常规水源可供水量进行调配;
2)发生中度干旱时在所述轻度干旱的调配基础上增添当地地表水、地下水应急水源供水量;
3)发生重度干旱时启用应急水源工程进行供水;
4)连续多年干旱时启用外调应急水源供水。
在上述任一方案中优选的是,所述步骤3为利用ArcGIS软件对DEM数据进行分析,生成干旱区河网,结合各水库经纬度坐标,确定各水库上下游关系,并对水库进行编码,识别干旱区各水库之间河道在轻旱、中旱、重旱时期内水流连通性,判断不同干旱程度下水库之间的水力联系。
在上述任一方案中优选的是,所述步骤3包括以下子步骤:
步骤31:根据DEM数据利用ArcGIS软件生成干旱区河网,确定水库之间的上下游关系;
步骤32:对水库进行编码;
步骤33:依据历史资料,判断不同干旱程度下各水库之间是否连通。
在上述任一方案中优选的是,所述步骤4包括以下子步骤:
步骤41:建立干旱区内不同干旱程度下各调配单元、供水节点、用水节点、输水系统和退水系统中的点、线、面要素之间的拓扑关系,绘制干旱区各要素关系图;
步骤42:对调配单元内引提水及应急水源工程统一建立一个虚拟供水节点;
步骤43:将所述虚拟供水节点融入到所述干旱区各要素关系图中,形成所述水资源系统网络图。
在上述任一方案中优选的是,所述供水节点包括水库、引水工程、提水工程和应急水源工程。
在上述任一方案中优选的是,所述用水节点包括生活、生产和生态用水。
在上述任一方案中优选的是,所述应急供水调配模式包括:
1)调配模式一,适用于轻度干旱,压缩需水量;
2)调配模式二,适用于中度干旱,在所述调配模式一的基础上进行调配单元内供水挖潜;
3)调配模式三,适用于重度干旱,在所述调配模式二的基础上启用应急水源工程;
4)调配模式四,适用于连续多年干旱,在所述调配模式三的基础上启用应急外调水。
在上述任一方案中优选的是,所述分析计算包括以下内容:
1)在需水量计算方面,确定不同干旱程度下生活、生产、生态抗旱用水定额,根据各时段t社会经济指标,采用定额法对需水量进行分析计算;
2)在可供水量计算方面,根据不同干旱程度及调配模式,分别进行计算。
在上述任一方案中优选的是,所述根据不同干旱程度及调配模式的计算方法如下:
1)轻度干旱可供水量为对应时段轻度干旱条件下的常规水源可供水量;
2)中度干旱可供水量包括对应时段中度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量,其中:
水库应急水量为水库死水位对应水库蓄水量;
地下水应急开采量应建立地下水流数值模型模拟不同开采程度下对地下水位变化及恢复情况,选取不产生区域水位持续下降情况下的可开采量;
3)重度干旱可供水量包括对应时段重度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量及应急水源地储存量;
4)连续多年干旱可供水量包括对应时段连续多年干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量、应急水源地储存量、应急外调水量,应急外调水量即调配单元外水库满足该水库所在单元用水需求后的可供水量,计算公式如下:
在上述任一方案中优选的是,所述用水优先序包括:
1)水库向各调配单元供水优先序依次为:水库所在调配单元用水、受灾程度较高调配单元用水、受灾程度较低调配单元用水;
2)水源供水优先序依次为:常规地表水、常规再生水、常规地下水、常规外调水、应急地表水、应急地下水、应急水源供水和应急外调水;
3)行业用水优先序依次为:生活用水、生产用水、生态用水。
在上述任一方案中优选的是,双层大范围长历时干旱应急供水调配模型包括数据输入模块、干旱等级评价模块、供需水模块和优化调配模块。
在上述任一方案中优选的是,所述数据输入模块用于输入t时段供需水、降水预测数据及各水库信息。
在上述任一方案中优选的是,所述干旱程度评价模块用于根据降水数据计算SPI指数,判断各调配单元t时段干旱程度,选取相应调配模式与调配规则。
在上述任一方案中优选的是,所述供需水模块用于根据干旱程度评价结果调取相应干旱程度及调度模式下所需的供需水数据,并依据方案需求对供需水量进行处理。
在上述任一方案中优选的是,所述优化调配模块用于以整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为下层优化目标,以整体缺水量最小作为上层优化目标,实现双层大范围长历时干旱应急供水调配。
在上述任一方案中优选的是,下层优化得到各水库向调配单元供水的最优方案,下层优化目标函数:
其中,表示t时段k调配单元i水库需要满足的最小需水量,表示i水库向k调配单元的供水量,表示i水库向k调配单元供水的单位效益,表示t时段i水库向干旱水库输水的单位成本,表示t时段从i水库调水的水量,表示t时段i水库是否调水,表示t时段i水库实际蓄水量和目标蓄水量之间的差值,T表示干旱期内总时段数。
在上述任一方案中优选的是,上层优化将下层优化得到的最优解作为参数,优化得到各调配单元内向用水节点供水的最优方案,上层优化目标函数:
在上述任一方案中优选的是,所述双层大范围长历时干旱应急供水调配模型的模型约束条件包括水库水量平衡约束、水库可供水量约束、水库输水能力约束、需水约束、可供水量约束和非负约束中至少一种。
在上述任一方案中优选的是,所述水库水量平衡约束的公式为:
在上述任一方案中优选的是,所述水库可供水量约束的公式为:
在上述任一方案中优选的是,所述水库输水能力约束的公式为:
在上述任一方案中优选的是,所述需水约束的公式为:
在上述任一方案中优选的是,所述可供水量约束的公式为:
在上述任一方案中优选的是,所述非负约束的公式为
本发明提出了一种面向大范围长历时干旱的应急供水调配方法,针对大范围长历时干旱情景,分区对常规水源和应急水源进行分类识别,根据长历时干旱从轻度干旱到连续多年干旱发展过程,对不同程度干旱情况采用不同的调配模式、调配规则,并进行相应供需水计算。采用整体缺水量最小、整体受灾损失与及经济成本最小、水库实际蓄水量与目标蓄水量最小作为目标函数构建双层优化调配模型,降低各区域之间及当前时段利益与未来利益之间的供水利益冲突,提高了干旱水资源调配的效率、公平性及可持续性,使调配方案更具合理性。
DEM数据是指数字高程模型数据。
ArcGIS软件是一款美国ESRI公司研发的地理信息系统软件。
附图说明
图1为按照本发明的面向大范围长历时干旱的应急供水调配方法的一优选实施例的流程图。
图2为按照本发明的面向大范围长历时干旱的应急供水调配方法的水源分类统计的一实施例的示意图。
图3为按照本发明的面向大范围长历时干旱的应急供水调配方法的不同干旱程度下应急供水调配模式的一实施例的示意图。
图4为按照本发明的面向大范围长历时干旱的应急供水调配方法的大范围长历时干旱应急供水调配模型构建方法的一实施例的流程图。
图5为按照本发明的面向大范围长历时干旱的应急供水调配方法的一优选实施例的系统网络示意图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的阐述。
实施例一
本发明提出一种面向大范围长历时干旱的应急供水调配模型,来解决大范围长历时干旱下的水资源配置问题,降低干旱带来的社会经济损失。
如图1所示,执行步骤110,针对大范围干旱区,利用ArcGIS软件对地市级行政分区和水资源分区进行剖分叠加,形成“水资源分区+行政分区”的调配单元,便于分析计算和政府部门决策管理。
执行步骤120,区分各调配单元内常规供水水源及应急供水水源,同时对两种水源中地表水、地下水、非常规水、外调水数量和基本参数(例如水库的水位库容曲线、水利工程供水能力)等进行分类统计,便于对不同干旱程度不同调配模式下可供水量进行计算,水源统计分类如图2所示。
进一步的便于在不同干旱程度下,提取相应水源的可供水量进行调配,包括:
1)在轻度干旱时使用常规水源可供水量进行调配;
2)中度干旱时在轻度干旱的基础上增添当地地表水、地下水应急水源供水量;
3)重度干旱时启用应急水源工程进行供水;
4)连续多年干旱时启用外调应急水源供水。
执行步骤130,利用ArcGIS软件对DEM数据进行分析,生成干旱区河网,结合各水库经纬度坐标,确定各水库上下游关系,并对水库进行编码,识别干旱区各水库之间河道在轻旱、中旱、重旱时期内水流连通性,判断不同干旱程度下水库之间的水力联系,包括以下内容:
1、根据DEM数据利用ArcGIS软件生成干旱区河网,确定水库之间的上下游关系;
2、对水库进行编码;
3、依据历史资料,判断不同干旱程度下各水库之间是否连通。
执行步骤140,建立各水源、调配单元之间的拓扑关系,绘制干旱区水资源系统网络图,包括以下步骤:
步骤141,建立干旱区内不同干旱程度下各调配单元、供水节点、用水节点、输水系统、退水系统等点、线、面要素之间的拓扑关系,绘制干旱区各要素关系图,其中供水节点包括水库、引水工程、提水工程、应急水源工程等;
步骤142,对调配单元内引提水及应急水源工程统一建立一个虚拟供水节点;
步骤143,用水节点包括生活、生产、生态用水等,隐含在调配单元中,形成水资源系统网络图。
执行步骤150,如图3所示,针对不同干旱程度选定相应的应急供水调配模式,应急供水调配模式包括:
1)调配模式一,适用于轻度干旱,压缩需水量;
2)调配模式二,适用于中度干旱,在所述调配模式一的基础上进行调配单元内供水挖潜;
3)调配模式三,适用于重度干旱,在所述调配模式二的基础上启用应急水源工程;
4)调配模式四,适用于连续多年干旱,在所述调配模式三的基础上启用应急外调水。
执行步骤160,对各调配单元不同干旱程度下可供水量及需水量进行分析计算,分析计算包括以下内容:
1、在需水量计算方面,确定不同干旱程度下生活、生产、生态抗旱用水定额,根据各时段t社会经济指标,采用定额法对需水量进行分析计算;
2、在可供水量计算方面,根据不同干旱程度及调配模式,分别进行计算,方法如下:
①轻度干旱可供水量为对应时段轻度干旱条件下的常规水源可供水量;
②中度干旱可供水量包括对应时段中度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量,其中:
水库应急水量为水库死水位对应水库蓄水量;
地下水应急开采量应建立地下水流数值模型模拟不同开采程度下对地下水位变化及恢复情况,选取不产生区域水位持续下降情况下的可开采量;
③重度干旱可供水量包括对应时段重度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量及应急水源地储存量;
④连续多年干旱可供水量包括对应时段连续多年干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量、应急水源地储存量、应急外调水量,应急外调水量即调配单元外水库满足该水库所在单元用水需求后的可供水量,计算公式如下:
执行步骤170,根据各调配单元不同干旱程度下选取的调配模式,针对干旱区社会经济发展及生态保护需求,制定各水库向各调配单元供水优先序,制定各调配单元内水源供水优先序及行业用水优先序,包括以下内容:
1、水库向各调配单元供水优先序依次为:水库所在调配单元用水、受灾程度较高调配单元用水、受灾程度较低调配单元用水;
2、水源供水优先序依次为:常规地表水、常规再生水、常规地下水、常规外调水、应急地表水、应急地下水、应急水源供水、应急外调水;
3、行业用水优先序依次为:生活用水、生产用水、生态用水。
执行步骤180,构建双层大范围长历时干旱应急供水调配模型,模型优化求解得到大范围长历时干旱最优调配方案,双层大范围长历时干旱应急供水调配模型包括数据输入模块、干旱等级评价模块、供需水模块和优化调配模块,模型构建流程见图4,各模块包括以下内容:
1、数据输入模块用于输入t时段供需水、降水预测数据及各水库信息。
2、干旱程度评价模块用于根据降水数据计算SPI指数,判断各调配单元t时段干旱程度,选取相应调配模式与调配规则。
3、供需水模块用于根据干旱程度评价结果调取相应干旱程度及调度模式下所需的供需水数据,并依据方案需求对供需水量进行处理。
4、优化调配模块用于以整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为下层优化目标,以整体缺水量最小作为上层优化目标,实现双层大范围长历时干旱应急供水调配。
下层优化得到各水库向调配单元供水的最优方案;上层优化将下层优化得到的最优解作为参数,优化得到各调配单元内向用水节点供水的最优方案。模型目标函数及约束条件如下:
上层优化目标函数:
下层优化目标函数:
其中,表示t时段k调配单元i水库需要满足的最小需水量,表示i水库向k调配单元的供水量,表示i水库向k调配单元供水的单位效益,表示t时段i水库向干旱水库输水的单位成本,表示t时段从i水库调水的水量,表示t时段i水库是否调水,表示t时段i水库实际蓄水量和目标蓄水量之间的差值,T表示干旱期内总时段数。
模型约束条件包括水库水量平衡约束、水库可供水量约束、水库输水能力约束、需水约束、可供水量约束和非负约束中至少一种,具体约束表达式如下:
①、水库水量平衡约束:
②、水库可供水量约束;
③、水库输水能力约束:
④、需水约束:
⑤、可供水量约束:
⑥、非负约束:
实施例二
某一干旱区内划分了A、B、C三个调配单元,区域内有编号为1、2、3的三座水库,调配过程如下:
①“水库-水库”及“水库-调配单元”之间拓扑关系如表1所示。
表1“水库-水库”及“水库-调配单元”之间拓扑关系
②绘制如图5所示的系统网络图。
③经判断,干旱区内某时段t时,调配单元A为轻度干旱,选取压缩需水量的调配模式;调配单元B为中度干旱,选取“压缩需水+调配单元内供水挖潜”的调配模式;调配单元C为重度干旱,选取“压缩需水+调配单元内供水挖潜+应急水源工程”的调配模式。
④各水库常规可供水量、应急可供水量、水库目标蓄水量及最小应满足需水量如表2所示。
表2水库调配数据参数
⑤经下层优化调配后,各水库向各调配单元供水量如表3所示。
表3时段t水库优化调配方案
⑥根据上层优化调配结果,将各调配单元常规供水量计入常规地表水源可供水量中,将应急供水量计入地表水挖潜供水量中,对各调配单元可供水量及需水量进行计算,如表4所示。
表4调配单元供需水计算
⑦经过下层优化调配后,各调配单元不同行业供水量如表5所示。
表5各调配单元行业用水优化调配方案
实施例三
双层大范围长历时干旱应急供水调配模型的构建求解方法包括以下步骤:
步骤1、输入t时段供需水、降水预测数据及各水库信息;
步骤2、根据降水数据计算SPI指数,判断干旱区t时段干旱程度,选取相应调配模式与调配规则;
步骤3、根据步骤2中干旱程度评价结果调取相应干旱程度及调配模式下所需的供需水数据,并依据方案需求对供需水量进行处理;
步骤4、以整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为下层优化目标,以整体缺水量最小作为上层优化目标,构建双层大范围长历时干旱应急供水调配模型。
下层优化得到各水库水源向调配单元供水的最优方案;上层优化将下层优化得到的最优解作为参数,优化得到各调配单元内各水源向各用水户供水的最优方案。模型目标函数及约束条件如下:
上层优化目标函数:
下层优化目标函数:
其中,表示t时段k调配单元i水库需要满足的最小需水量,表示i水库向k调配单元的供水量,表示i水库向k调配单元供水的单位效益,表示t时段i水库向干旱水库输水的单位成本,表示t时段从i水库调水的水量,表示t时段i水库是否调水,用0或1表示;表示t时段i水库实际蓄水量和目标蓄水量之间的差值,T表示干旱期内总时段数
模型约束条件包括水库水量平衡约束、水库可供水量约束、水库输水能力约束、需水约束、可供水量约束等,具体约束表达式如下:
①、水库水量平衡约束:
②、水库可供水量约束;
③、水库输水能力约束:
④、需水约束:
⑤、可供水量约束:
⑥、非负约束:
步骤5、模型求解得到大范围长历时干旱最优调配方案。
为了更好地理解本发明,以上结合本发明的具体实施例做了详细描述,但并非是对本发明的限制。凡是依据本发明的技术实质对以上实施例所做的任何简单修改,均仍属于本发明技术方案的范围。本说明书中每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似的部分相互参见即可。对于系统实施例而言,由于其与方法实施例基本对应,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
Claims (8)
1.一种面向大范围长历时干旱的应急供水调配方法,包括以下步骤:
步骤1:根据行政分区和水资源分区对大范围干旱区进行划分,得到调配单元;
步骤2:对所述调配单元中的常规水源和应急水源进行分类统计;
步骤3:对干旱区利用ArcGIS软件及DEM数据对不同干旱程度下各水库之间水力联系进行识别,确定上下游关系,并对水库进行编码;
步骤4:建立各水源、调配单元之间的拓扑关系,绘制干旱区水资源系统网络图;
步骤5:针对不同干旱程度选定相应的应急供水调配模式;
步骤6:对各调配单元不同干旱程度下可供水量及需水量进行分析计算;
步骤7:制定不同调配模式下各水源供水及各行业用水优先序;
步骤8:构建双层大范围长历时干旱应急供水调配模型,优化求解得到调配方案;
其中,上层优化将下层优化得到的最优解作为参数,优化得到各调配单元内向用水节点供水的最优方案,上层优化目标函数:
下层优化得到各水库向调配单元供水的最优方案,下层优化目标函数:
2.如权利要求1所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述分类统计的结果还用于在不同干旱程度下,提取相应水源的可供水量进行调配,包括:
1)发生轻度干旱时使用常规水源可供水量进行调配;
2)发生中度干旱时在所述轻度干旱的调配基础上增添当地地表水、地下水应急水源供水量;
3)发生重度干旱时启用应急水源工程进行供水;
4)连续多年干旱时启用外调应急水源供水。
3.如权利要求2所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述步骤3包括以下子步骤:
步骤31:根据DEM数据利用ArcGIS软件生成干旱区河网,确定水库之间的上下游关系;
步骤32:对水库进行编码;
步骤33:依据历史资料,判断不同干旱程度下各水库之间是否连通。
4.如权利要求3所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述步骤4包括以下子步骤:
步骤41:建立干旱区内不同干旱程度下各调配单元、供水节点、用水节点、输水系统和退水系统中的点、线、面要素之间的拓扑关系,绘制干旱区各要素关系图;
步骤42:对调配单元内引提水及应急水源工程统一建立一个虚拟供水节点;
步骤43:将所述虚拟供水节点融入到所述干旱区各要素关系图中,形成所述水资源系统网络图。
5.如权利要求4所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述分析计算包括以下内容:
1)在需水量计算方面,确定不同干旱程度下生活、生产、生态抗旱用水定额,根据各时段t社会经济指标,采用定额法对需水量进行分析计算;
2)在可供水量计算方面,根据不同干旱程度及调配模式,分别进行计算。
6.如权利要求5所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述根据不同干旱程度及调配模式的计算方法如下:
1)轻度干旱可供水量为对应时段轻度干旱条件下的常规水源可供水量;
2)中度干旱可供水量包括对应时段中度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量,其中:
水库应急水量为水库死水位对应水库蓄水量;
地下水应急开采量应建立地下水流数值模型模拟不同开采程度下对地下水位变化及恢复情况,选取不产生区域水位持续下降情况下的可开采量;
3)重度干旱可供水量包括对应时段重度干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量及应急水源地储存量;
4)连续多年干旱可供水量包括对应时段连续多年干旱条件下的常规水源可供水量、各调配单元内水库应急水量、各调配单元内在不产生区域水位持续下降为约束条件下的地下水应急开采量、应急水源地储存量、应急外调水量,应急外调水量即调配单元外水库满足该水库所在单元用水需求后的可供水量,计算公式如下:
7.如权利要求6所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述双层大范围长历时干旱应急供水调配模型包括数据输入模块、干旱等级评价模块、供需水模块和优化调配模块。
8.如权利要求7所述的面向大范围长历时干旱的应急供水调配方法,其特征在于,所述优化调配模块用于以整体受灾损失和经济成本最小、水库实际蓄水量与目标蓄水量偏差最小作为下层优化目标,以整体缺水量最小作为上层优化目标,实现双层大范围长历时干旱应急供水调配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110695430.XA CN113256028B (zh) | 2021-06-23 | 2021-06-23 | 一种面向大范围长历时干旱的应急供水调配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110695430.XA CN113256028B (zh) | 2021-06-23 | 2021-06-23 | 一种面向大范围长历时干旱的应急供水调配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113256028A CN113256028A (zh) | 2021-08-13 |
CN113256028B true CN113256028B (zh) | 2021-10-01 |
Family
ID=77189285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110695430.XA Active CN113256028B (zh) | 2021-06-23 | 2021-06-23 | 一种面向大范围长历时干旱的应急供水调配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113256028B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115099477B (zh) * | 2022-06-15 | 2023-02-03 | 中国水利水电科学研究院 | 一种水库旱限水位优化及抗旱调度方法 |
CN115797101A (zh) * | 2022-11-09 | 2023-03-14 | 江苏省水利科学研究院 | 湿润地区山丘区现状工程条件下历史可供水量折算方法 |
CN117669978A (zh) * | 2023-12-14 | 2024-03-08 | 中国水利水电科学研究院 | 一种面向不同干旱程度的生活需水分级阈值的调控方法 |
CN118037120A (zh) * | 2024-02-27 | 2024-05-14 | 青海九零六工程勘察设计院有限责任公司 | 一种河谷地下水源地应急供水能力评价方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106545048A (zh) * | 2016-12-08 | 2017-03-29 | 中国水利水电科学研究院 | 一种基于山区的层次化农村供水水源优化配置方法 |
CN107657329A (zh) * | 2017-07-04 | 2018-02-02 | 江苏鸿利智能科技有限公司 | 基于极端气候条件下防汛防旱的智能调度决策方法 |
CN108797511A (zh) * | 2018-06-19 | 2018-11-13 | 中国水利水电科学研究院 | 一种面向河道生态需水的水资源配置的方法 |
CN111105170A (zh) * | 2019-12-31 | 2020-05-05 | 张旭 | 水资源仿真配置计算方法及水资源的配置方法 |
CN112443010A (zh) * | 2019-08-29 | 2021-03-05 | 刘璐 | 一种城市多功能应急储水和供水系统 |
CN112443023A (zh) * | 2019-08-29 | 2021-03-05 | 刘璐 | 一种城市多功能防涝储水设施 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7736557B2 (en) * | 2003-12-02 | 2010-06-15 | Evans Brothers Investments | Aged roofing tile system |
US20080162047A1 (en) * | 2007-01-03 | 2008-07-03 | Nepal Chandra Mondal | Cross correlation technique to delineate groundwater recharge potential zone in hard rock terrain |
CN106950613A (zh) * | 2017-01-10 | 2017-07-14 | 于洋 | 一种干旱区灌溉水资源平衡分析系统 |
CN108681974B (zh) * | 2018-05-29 | 2021-11-02 | 中国水利水电科学研究院 | 一种面向内陆干旱区的水资源多维均衡配置方法 |
CN108764573B (zh) * | 2018-05-29 | 2021-11-02 | 中国水利水电科学研究院 | 一种面向内陆干旱区的水资源多维均衡配置系统 |
-
2021
- 2021-06-23 CN CN202110695430.XA patent/CN113256028B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106545048A (zh) * | 2016-12-08 | 2017-03-29 | 中国水利水电科学研究院 | 一种基于山区的层次化农村供水水源优化配置方法 |
CN107657329A (zh) * | 2017-07-04 | 2018-02-02 | 江苏鸿利智能科技有限公司 | 基于极端气候条件下防汛防旱的智能调度决策方法 |
CN108797511A (zh) * | 2018-06-19 | 2018-11-13 | 中国水利水电科学研究院 | 一种面向河道生态需水的水资源配置的方法 |
CN112443010A (zh) * | 2019-08-29 | 2021-03-05 | 刘璐 | 一种城市多功能应急储水和供水系统 |
CN112443023A (zh) * | 2019-08-29 | 2021-03-05 | 刘璐 | 一种城市多功能防涝储水设施 |
CN111105170A (zh) * | 2019-12-31 | 2020-05-05 | 张旭 | 水资源仿真配置计算方法及水资源的配置方法 |
Non-Patent Citations (6)
Title |
---|
Study on water resources allocation and planning in circular economy pilot zones;Fu-xin Chai;《IEEE》;20110616;全文 * |
基于遗传算法的多目标水资源优化配置模型在岩溶地区的应用;马兴华等;《中国防汛抗旱》;20180831;第28卷(第8期);全文 * |
多水源供给模式下灌区旱灾诱发危机诊断与预警——以泾惠渠灌区为例;王斌;《中国优秀硕士学位论文全文数据库 (基础科学辑)》;20180430;D046-2 * |
干旱条件下城市水资源应急调配模型;邵东国等;《华北水利水电大学学报(自然科学版)》;20140228;第35卷(第1期);全文 * |
引汉济渭跨流域调水工程运行调度模式研究及实现;张晓;《中国博士学位论文全文数据库 (工程科技Ⅱ辑)》;20190831;C037-14 * |
江苏省南水北调受水区水资源配置;宋丹丹等;《 南水北调与水利科技》;20150331;第13卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113256028A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113256028B (zh) | 一种面向大范围长历时干旱的应急供水调配方法 | |
Kucukmehmetoglu et al. | International water resources allocation and conflicts: the case of the Euphrates and Tigris | |
Hedden et al. | Parched prospects-the emerging water crisis in South Africa | |
Liu et al. | A model for the optimal allocation of water resources in a saltwater intrusion area: a case study in Pearl River Delta in China | |
CN105297827A (zh) | 一种考虑多用户需水和多水源供水的水资源配置方法 | |
CN106485366A (zh) | 一种复杂梯级水库群蓄水期优化调度方法 | |
CN108985585B (zh) | 一种考虑预报不确定性影响的水库洪水资源利用风险对冲决策方法 | |
CN104047258B (zh) | 一种面向环境保护的联合调度水库群的生态库容确定方法 | |
Dimova et al. | Complementary use of the WEAP model to underpin the development of SEEAW physical water use and supply tables | |
CN114663254B (zh) | 一种水资源-粮食-能源-生态协同调控方法 | |
CN104636834B (zh) | 一种改进的联合概率规划模型系统优化方法 | |
CN110533326A (zh) | 一种与水相关的生态环境承载力预估调控系统 | |
Huang et al. | An optimization model for water resources allocation in Dongjiang River Basin of Guangdong-Hong Kong-Macao Greater Bay Area under multiple complexities | |
CN110633849B (zh) | 一种基于启发式大系统分解协调思想的水资源配置方法 | |
Rak et al. | Assessment of water supply diversification using the Pielou index | |
CN111915065A (zh) | 一种河流枯水期多目标动态水资源优化配置系统及方法 | |
Heydari Kushalshah et al. | Hybrid modelling for urban water supply system management based on a bi-objective mathematical model and system dynamics: A case study in Guilan province. | |
CN113762774A (zh) | 一种矿区所在地矿井水配置方法 | |
CN117391334A (zh) | 一种适用于大型旱区的水资源统一优化配置方法 | |
Shijie et al. | Quantifying the degree of water resource utilization polarization: a case study of the Yellow River Basin | |
Arena et al. | A simulation/optimization model for selecting infrastructure alternatives in complex water resource systems | |
Huang et al. | A stochastic optimization approach for integrated urban water resource planning | |
Efstratiadis et al. | Minimizing water cost in water resource management of Athens | |
CN113887073A (zh) | 流域水土的资源优化配置方法、系统、电子设备及存储介质 | |
Davidsen et al. | The cost of ending groundwater overdraft on the North China Plain. |
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 |