CN110135627A - 水资源优化方法及装置 - Google Patents
水资源优化方法及装置 Download PDFInfo
- Publication number
- CN110135627A CN110135627A CN201910328856.4A CN201910328856A CN110135627A CN 110135627 A CN110135627 A CN 110135627A CN 201910328856 A CN201910328856 A CN 201910328856A CN 110135627 A CN110135627 A CN 110135627A
- Authority
- CN
- China
- Prior art keywords
- water
- subregion
- resource optimization
- water resource
- target
- 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 306
- 238000005457 optimization Methods 0.000 title claims abstract description 165
- 238000000034 method Methods 0.000 title claims abstract description 95
- 239000002352 surface water Substances 0.000 claims abstract description 132
- 230000008901 benefit Effects 0.000 claims abstract description 87
- 238000013439 planning Methods 0.000 claims abstract description 79
- 238000003973 irrigation Methods 0.000 claims abstract description 48
- 230000002262 irrigation Effects 0.000 claims abstract description 48
- 230000006870 function Effects 0.000 claims description 65
- 238000007726 management method Methods 0.000 claims description 14
- 230000007704 transition Effects 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 7
- 230000001052 transient effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 230000036642 wellbeing Effects 0.000 claims description 3
- 206010049040 Weight fluctuation Diseases 0.000 claims description 2
- 230000008859 change Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000003621 irrigation water Substances 0.000 description 4
- 230000001154 acute effect Effects 0.000 description 2
- 238000010612 desalination reaction Methods 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 238000012913 prioritisation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- JXASPPWQHFOWPL-UHFFFAOYSA-N Tamarixin Natural products C1=C(O)C(OC)=CC=C1C1=C(OC2C(C(O)C(O)C(CO)O2)O)C(=O)C2=C(O)C=C(O)C=C2O1 JXASPPWQHFOWPL-UHFFFAOYSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009313 farming Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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/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/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例提供一种水资源优化方法及装置。其中,方法包括:建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为第二水资源优化模型;根据第二水资源优化模型和预设的违背风险,获取水资源优化方案以使得目标区域规划时间内的综合效益最大。本发明实施例提供的水资源优化方法及装置,能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
Description
技术领域
本发明涉及农业技术领域,更具体地,涉及一种水资源优化方法及装置。
背景技术
水资源管理问题往往需要满足多方面、不同利益部门的需求,涉及多个目标之间的权衡问题。并且,水资源管理中还普遍存在着许多的不确定性因素,如可用水资源量等。
线性加权方法是最常用适用于多目标优化问题的方法。使用线性加权方法的过程中,若采用客观赋权法确定权重,仅利用数据本身的信息,则忽略了决策者偏好和社会影响等外界因素的影响;若采用主观赋权法确定权重,结果会受到决策者主观偏好的影响,存在着不确定性;若采用主客观结合的方法确定权重,权重的赋值同样会受到决策者主观因素的影响。由于多目标权重存在着不确定性信息,而现有的赋权方法忽略了这样的不确定信息。
由于水资源优化包括多目标,多目标权重具有不确定性以及可用水资源量等参数具有不确定性,在多目标权重和参数具有不确定性的情况下获得水资源优化方案是亟待解决的问题。
发明内容
本发明实施例提供一种水资源优化方法及装置,用以解决或者至少部分地解决现有技术中无法在多目标权重和参数具有不确定性的情况下获得水资源优化方案的缺陷。
第一方面,本发明实施例提供一种水资源优化方法,包括:
获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,获取每一所述子区域的种植面积上限和种植面积下限,获取每一所述子区域内各种作物的产量需求上限和产量需求下限;
根据所述每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,每一所述子区域的种植面积上限和种植面积下限,以及每一所述子区域内各种作物的产量需求上限和产量需求下限,建立以最大化所述目标区域所述规划时间内的经济目标和生态目标的第一水资源优化模型;其中,所述第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;
对于每一采用地表水灌溉的子区域,根据所述采用地表水灌溉的子区域的地表水量历史数据,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;
根据线性加权方法和鲁棒优化方法将所述第一水资源优化模型转化为以最大化所述目标区域所述规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以所述标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中;所述第二水资源优化模型以转化后的所述第一水资源优化模型中的约束条件为约束条件;
根据所述第二水资源优化模型和预设的违背风险,获取每一所述子区域内在所述规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得所述目标区域所述规划时间内的综合效益最大。
第二方面,本发明实施例提供一种水资源优化装置,包括:
数据获取模块,用于获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,获取每一所述子区域的种植面积上限和种植面积下限,获取每一所述子区域内各种作物的产量需求上限和产量需求下限;
模型建立模块,用于根据所述每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,每一所述子区域的种植面积上限和种植面积下限,以及每一所述子区域内各种作物的产量需求上限和产量需求下限,建立以最大化所述目标区域所述规划时间内的经济目标和生态目标的第一水资源优化模型;其中,所述第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;
条件转化模块,用于对于每一采用地表水灌溉的子区域,根据所述采用地表水灌溉的子区域的地表水量历史数据,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;
模型转化模块,用于根据线性加权方法和鲁棒优化方法将所述第一水资源优化模型转化为以最大化所述目标区域所述规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以所述标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中;所述第二水资源优化模型以转化后的所述第一水资源优化模型中的约束条件为约束条件;
方案获取模块,用于根据所述第二水资源优化模型和预设的违背风险,获取每一所述子区域内在所述规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得所述目标区域所述规划时间内的综合效益最大。
第三方面,本发明实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,执行所述程序时实现如第一方面的各种可能的实现方式中任一种可能的实现方式所提供的水资源优化方法的步骤。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如第一方面的各种可能的实现方式中任一种可能的实现方式所提供的水资源优化方法的步骤。
本发明实施例提供的水资源优化方法及装置,通过将优化目标的标准化权重扩展为随机分布于名义值为中心,以波动半径为半径的对称区间中来消除目标权重的不确定性,并根据地表水量历史数据消除子区域地表水可用量的不确定性,从而能消除多目标权重和参数的不确定性,能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为根据本发明实施例提供的水资源优化方法的流程示意图;
图2为根据本发明实施例提供的水资源优化装置的结构示意图;
图3为根据本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了克服现有技术的上述问题,本发明实施例提供一种水资源优化方法及装置,其发明构思是,通过将优化目标的标准化权重扩展为随机分布于名义值为中心,以波动半径为半径的对称区间中来消除目标权重的不确定性,并根据地表水量历史数据消除子区域地表水可用量的不确定性,从而能消除多目标权重和参数的不确定性,获得水资源优化方案。
图1为根据本发明实施例提供的水资源优化方法的流程示意图。如图1所示,方法包括:步骤S101、获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,获取每一子区域的种植面积上限和种植面积下限,获取每一子区域内各种作物的产量需求上限和产量需求下限。
具体地,目标区域包括多个子区域,子区域可以称为灌区。目标区域内种植有多种作物,每一子区域内可以种植上述多种作物中的一种或多种。目标区域由多种灌溉水源对作物进行灌溉,上述多种灌溉水源至少包括地表水和地下水,还可以包括再生水、淡化后的水资源等其它非常规水源,本发明实施例对地表水和地下水之外的灌溉水源不作具体限制。对于每一子区域,可以采用上述多种灌溉水源中的至少一种进行灌溉。
本发明实施例提供的水资源优化方法,是为了获取规划时间(例如未来十年)内每一单位时间(例如每一年)对目标区域的各子区域进行配水的方案。
可以用i表示作物种类的编号;j表示灌溉水源的编号(j=1表示地下水,j=2表示地表水);k表示子区域的编号;t表示规划时间内单位时间的编号;I表示作物种类的数量;J表示灌溉水源的数量;K表示子区域的数量;T表示规划时间包括的单位时间的数量;Aijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的灌溉面积。
可以根据目标区域的各子区域的历史数据,获取每一子区域内每一种作物的单位经济净效益NBijkt、单位产量YIit和单位生态净效益CFijk,获取目标区域在规划时间内每一单位时间的总可用水量TAt和每一子区域的地下水可用量GWk,获取每一子区域的种植面积下限和种植面积上限获取每一子区域内每一种作物的产量需求下限和产量需求上限
其中,NBijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的单位经济净效益;YIit表示第t个单位时间第i种作物的单位产量;CFijk表示第k个子区域内用第j种灌溉水源灌溉第i种作物的单位生态净效益;TAt表示目标区域在第t个单位时间的总可用水量;GWk表示第k个子区域的地下水可用量;和分别表示第k个子区域内第i种作物的产量需求下限和上限。
单位经济净效益NBijkt,可以根据第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物历年的单位经济净效益进行预测获取。
单位产量YIit,可以根据目标区域历年的历史单位产量进行预测获取。
单位生态净效益CFijk,可以根据第k个子区域内用第j种灌溉水源灌溉第i种作物的历史单位生态净效益进行预测获取。
总可用水量TAt,可以目标区域历年的历史总可用水量进行预测获取。
地下水可用量GWk,可以根据第k个子区域历年的地下水可用量进行预测获取或直接获取。由于地下水可用量具有一定的稳定性,可以认为第k个子区域每一单位时间的地下水可用量相同。例如,可以将第k个子区域历年的地下水可用量的最大值作为GWk。
种植面积下限和种植面积上限可以根据第k个子区域历年各种作物种植面积的总和进行预测获取或直接获取。由于如果耕地等农用地不发生剧烈变化,总种植面积具有一定的稳定性,可以认为第k个子区域每一单位时间的种植面积下限相同,每一单位时间的种植面积上限相同。例如,可以将第k个子区域内历年各种作物种植面积的总和中的最小值作为历年各种作物种植面积的总和中的组最大值作为
产量需求下限和产量需求上限可以根据第k个子区域内第i种作物历年的产量需求进行预测获取或直接获取。由于如果人口、市场等因素不发生剧烈变化,作物的产量需求具有一定的稳定性,可以认为第k个子区域内第i种作物每一单位时间的产量需求下限相同,第k个子区域内第i种作物每一单位时间的产量需求上限相同。例如,可以将第k个子区域内第i种作物历年的产量需求中的最小值作为历年的产量需求中的最大值作为
可以采用机器学习的方法基于历史数据分别获得多个预测模型,并根据多个预测模型分别预测出NBijkt、YIit、CFijk、TAt、GWk、和
步骤S102、根据每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,每一子区域的种植面积上限和种植面积下限,以及每一子区域内各种作物的产量需求上限和产量需求下限,建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;其中,第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件。
具体地,根据每一组i,j,k,t取值下的NBijkt、YIit、CFijk、TAt、GWk、和可以建立第一水资源优化模型。
第一水资源优化模型包括两个目标函数,分别为经济目标函数和生态目标函数。
经济目标函数为
生态目标函数为
第一水资源优化模型中的约束条件至少包括总水量约束、地表水约束、地下水约束、种植面积约束和产量需求约束。
总水量约束为
若第k个子区域采用了地表水灌溉,则地表水约束为
其中,SW(u)表示地表水可用量SW在概率空间U上的随机参数,u∈U。
若第k个子区域采用了地下水灌溉,则地下水约束为
若第k个子区域不利用地表水,则该子区域的Aijkt=0,j=2,
第k个子区域的种植面积约束为
第k个子区域的产量需求约束为
需要说明的是,对于再生水、淡化后的水资源等其它非常规水源,其约束条件与地下水约束类似,此处不再赘述。
步骤S103、对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件。
具体地,由于地表水可用量具有不确定性,所以在地表水约束中用概率分布的形式表示地表水可用量。
为了消除地表水可用量的不确定性,对地表水可用量进行机会约束规划。
进行机会约束规划的具体步骤包括:对于每一采用地表水灌溉的子区域(假设为第k个子区域),可以根据该子区域地表水量历史数据,获取某一置信水平p下该子区域的地表水可用量,将该子区域的地表水约束中的SW(u)替换为某一置信水平p下该子区域的地表水可用量SW(p),从而将该子区域的地表水量约束条件转化为预设的置信水平p下的确定性约束条件。转化后的地表水量约束条件为
该子区域地表水量历史数据,为该子区域历年地表水量的观测值。
置信水平p的取值范围为p∈[0,1]。
可以理解的是,除了地表水约束之外,第一水资源优化模型中的其他约束条件均为确定性约束条件,不需要进行转化。
步骤S104、根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;第二水资源优化模型以转化后的第一水资源优化模型中的约束条件为约束条件。
具体地,由于第一水资源优化模型包括两个目标函数(经济目标函数和生态目标函数),最大化经济目标和最大化生态目标之间存在一定的矛盾,可以根据线性加权方法对经济目标函数和生态目标函数进行加权求和,获得综合效益目标函数,通过最大化综合效益,达到获最大化经济目标和最大化生态目标的平衡。
加权求和之后,对标准化权重进行不确定性处理消除不确定性。
标准化权重,指相对权重,即以其中一个目标函数的权重为1,将另一个目标函数的权重与该目标函数的权重的比值,作为上述另一个目标函数的标准化权重。
对综合效益目标函数中经济目标函数或生态目标函数的标准化权重进行不确定性处理,将该目标函数的标准化权重W处理为表示W随机分布在对称区间中。其中,表示W的名义值(即根据预设的经济目标函数和生态目标函数的权重值获取的二者的比值);表示预设的W的波动半径。
根据鲁棒优化方法将对标准化权重进行不确定性处理后的综合效益目标函数进行线性转化,获得第二水资源优化模型。
第一水资源模型中的约束条件进行转化后,直接作为第二水资源模型中的约束条件。
需要说明的是,步骤S103和步骤S104之间的时序关系不限于先执行步骤S103后执行步骤S104,也可以同时执行步骤S103和步骤S104,或者先执行步骤S104后执行步骤S103。
步骤S105、根据第二水资源优化模型和预设的违背风险,获取每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得目标区域规划时间内的综合效益最大。
通过步骤S103和步骤S104获得的第二水资源优化模型,是约束条件均为确定性约束条件、目标函数为线性单目标函数的模型;由于鲁棒优化方法引入了违背风险(Probability bound of unattainability,PBU),可以根据预设的违背风险,很方便地获得第二水资源优化模型的一组最优解作为水资源优化方案。
一组最优解为目标区域的每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积。在该组最优解下,目标区域规划时间内的综合效益最大。
综合效益将随着目标权重(包括标准化权重的名义值和波动半径)和地表水可用量的变化而变化,能够清晰地反映出综合效益对随机参数和目标权重两种不确定性信息变化的响应情况。具体而言,综合效益会随着标准化权重波动范围的减小、目标可行性风险的增加或者地表水可用量约束违背风险的增加而增加。第二水资源优化模型还可以提供各种作物最适宜的灌溉水源及最佳的配水量,从而为区域管理者提供更为明确的管理方案。
本发明实施例通过将优化目标的标准化权重扩展为随机分布于名义值为中心,以波动半径为半径的对称区间中来消除目标权重的不确定性,并根据地表水量历史数据消除子区域地表水可用量的不确定性,从而能消除多目标权重和参数的不确定性,能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。进一步地,获取的水资源优化方案具有稳健性,获取水资源优化方案更加简便、耗时更短。
基于上述各实施例的内容,根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型的具体步骤包括:根据线性加权方法,将第一水资源优化模型中的两个目标函数以预设的经济目标和生态目标的原始权重进行加权,转化为以最大化目标区域规划时间内的经济目标和生态目标的加权和为目标的过渡目标函数。
具体地,根据预设的经济目标的原始权重和预设的生态目标的原始权重,对第一水资源优化模型中的经济目标函数和生态目标函数进行加权,获得过渡目标函数。
过渡目标函数为
其中,Wl表示过渡目标函数Fl(x)的原始权重。
l=1,Fl(x)为l=2,Fl(x)为或者l=1,Fl(x)为l=2,Fl(x)为
将经济目标和生态目标的原始权重进行不确定性处理,将过渡目标函数转化为单目标函数。
具体地,将过渡目标函数中经济目标和生态目标的原始权重进行标准化,并对标准化后的原始权重进行不确定性处理,消除原始权重的不确定性,将过渡目标函数转化为单目标函数。
根据鲁棒优化方法将单目标函数线性化,获得第二水资源优化模型中的目标函数,并将转化后的第一水资源优化模型中的约束条件作为第二水资源优化模型中的约束条件。
根据鲁棒优化方法将对标准化权重进行不确定性处理后获得的单目标函数线性化,获得第二水资源优化模型中的目标函数。
可以将第一水资源模型中的约束条件进行转化后,直接作为第二水资源模型中的约束条件。
获得第二水资源优化模型中的目标函数和约束条件之后,即可获得为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型。
本发明实施例通过线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型,能消除多目标权重的不确定性,从而能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
基于上述各实施例的内容,将经济目标和生态目标的原始权重进行不确定性处理的具体步骤包括:将经济目标的原始权重和生态目标的原始权重,分别除以经济目标的原始权重,获得生态目标的标准化权重的名义值,将生态目标的标准化权重随机分布于以生态目标的标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;
或者,将经济目标的原始权重和生态目标的原始权重,分别除以生态目标的原始权重,获得经济目标的标准化权重的名义值,将经济目标的标准化权重随机分布于以经济目标的标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中。
具体地,获得过渡目标函数,可以将过渡目标函数中经济目标的原始权重和生态目标的原始权重,均除以经济目标的原始权重,获得经济目标的标准化权重的名义值为1、生态目标的标准化权重的名义值
对生态目标的原始权重的名义值进行不确定性处理,使生态目标的标准化权重随机分布在以名义值为中心、以波动半径为半径的对称区间中。
同理,也可以将过渡目标函数中经济目标的原始权重和生态目标的原始权重,均除以生态目标的原始权重,获得生态目标的标准化权重的名义值为1、经济目标的标准化权重的名义值
对经济目标的原始权重的名义值进行不确定性处理,使经济目标的标准化权重随机分布在以名义值为中心、以波动半径为半径的对称区间中。
本发明实施例通过将经济目标或生态目标的标准化权重扩展为随机分布于名义值为中心,以波动半径为半径的对称区间中来消除权重的不确定性,从而能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
基于上述各实施例的内容,第二水资源优化模型中的目标函数为
max z
s.t.
q≥0
rijkt≥0,
其中,z表示目标区域规划时间内的综合效益;i表示作物种类的编号;j表示灌溉水源的编号;k表示子区域的编号;t表示规划时间内单位时间的编号;I表示作物种类的数量;J表示灌溉水源的数量;K表示子区域的数量;T表示规划时间包括的单位时间的数量;Aijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的灌溉面积;表示Aijkt的经济收益;表示Aijkt的生态收益;Γ表示保护度水平;q表示第一辅助变量;rijkt表示Aijkt对应的第二辅助变量;表示标准化权重的名义值;表示标准化权重的波动半径。
具体地,进行原始权重标准化时,除以经济目标的原始权重,根据线性加权方法和鲁棒优化方法进行转化获得的第二水资源优化模型中的目标函数为
maxz
s.t.
q≥0
rijkt≥0,
其中,z表示目标区域规划时间内的综合效益;i表示作物种类的编号;j表示灌溉水源的编号;k表示子区域的编号;t表示规划时间内单位时间的编号;I表示作物种类的数量;J表示灌溉水源的数量;K表示子区域的数量;T表示规划时间包括的单位时间的数量;Aijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的灌溉面积;表示Aijkt的经济收益;表示Aijkt的生态收益;Γ表示保护度水平;q表示第一辅助变量;rijkt表示Aijkt对应的第二辅助变量;表示标准化权重的名义值;表示标准化权重的波动半径。
第一辅助变量和各第二辅助变量,可以在对第二水资源优化模型进行最优化求解的过程中确定。
对于每一Aijkt,和分别为目标区域经济目标和生态目标的子目标,即分别为经济子目标和生态子目标。
可以理解的是,经济子目标和生态子目标的个数均为物种类的数量、灌溉水源的数量、子区域的数量和规划时间包括的单位时间的数量四者的乘积。
本发明实施例中第二水资源优化模型中的目标函数为线性单目标函数,可以很方便地获得该目标函数在第二水资源优化模型中的约束条件下的最优解,能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
基于上述各实施例的内容,保护度水平根据预设的违背风险获取。
具体地,保护度水平用于表示可以保护的子目标标准化权重的最大数量。保护度水平与违背风险相关。违背风险,用于表示第二水资源优化模型中的目标函数的可达性。
本发明实施例引入保护度水平,可以获取不同违背风险下的水资源优化方案,使得水资源优化结果具有更好的适应性。
基于上述各实施例的内容,根据预设的违背风险获取保护度水平的公式为
其中,PBU表示预设的违背风险;v=(Γ+n)/2;
n=I×J×K×T。
具体地,可以根据预设的违背风险确定该违背风险下的保护度水平Γ。具体公式为
其中,PBU表示预设的违背风险;v=(Γ+n)/2;
n=I×J×K×T。
B表示贝塔分布;表示v向下取整。
本发明实施例根据预设的违背风险获取保护度水平,从而可以获取不同违背风险下的水资源优化方案,使得水资源优化结果具有更好的适应性。
基于上述各实施例的内容,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件的具体步骤包括:根据采用地表水灌溉的子区域的地表水量历史数据,获取采用地表水灌溉的子区域的地表水量的分布模型。
具体地,对于任一采用地表水灌溉的子区域,可以根据该子区域的地表水量历史数据,获取采用地表水灌溉的子区域的地表水量的分布模型。
例如,可以根据该子区域的地表水量历史数据绘制累计分布曲线,获得该子区域的地表水量随时间的分布模型。
根据分布模型,获取预设的置信水平下采用地表水灌溉的子区域单位时间的地表水可用量。
具体地,根据预设的置信水平,可以获得该分布模型中预设的置信水平下采用地表水灌溉的子区域单位时间的地表水可用量。
根据预设的置信水平下采用地表水灌溉的子区域单位时间的地表水可用量,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件。
具体地,将该子区域的地表水约束中的SW(u)替换为某一置信水平p下该子区域的地表水可用量SW(p),从而将该子区域的地表水量约束条件转化为预设的置信水平p下的确定性约束条件。
本发明实施例根据地表水量历史数据获取某一置信水平下的地表水可用量,替换概率空间中的随机参数,能消除子区域地表水可用量的不确定性,从而能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。
图2为根据本发明实施例提供的水资源优化装置的结构示意图。基于上述各实施例的内容,如图2所示,该装置包括数据获取模块201、模型建立模块202、条件转化模块203、模型转化模块204和方案获取模块205,其中:
数据获取模块201,用于获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,获取每一子区域的种植面积上限和种植面积下限,获取每一子区域内各种作物的产量需求上限和产量需求下限;
模型建立模块202,用于根据每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,每一子区域的种植面积上限和种植面积下限,以及每一子区域内各种作物的产量需求上限和产量需求下限,建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;其中,第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;
条件转化模块203,用于对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;
模型转化模块204,用于根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;第二水资源优化模型以转化后的第一水资源优化模型中的约束条件为约束条件;
方案获取模块205,用于根据第二水资源优化模型和预设的违背风险,获取每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得目标区域规划时间内的综合效益最大。
具体地,数据获取模块201,可以根据目标区域的各子区域的历史数据,获取每一子区域内每一种作物的单位经济净效益NBijkt、单位产量YIit和单位生态净效益CFijk,获取目标区域在规划时间内每一单位时间的总可用水量TAt和每一子区域的地下水可用量GWk,获取每一子区域的种植面积下限和种植面积上限获取每一子区域内每一种作物的产量需求下限和产量需求上限
其中,NBijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的单位经济净效益;YIit表示第t个单位时间第i种作物的单位产量;CFijk表示第k个子区域内用第j种灌溉水源灌溉第i种作物的单位生态净效益;TAt表示目标区域在第t个单位时间的总可用水量;GWk表示第k个子区域的地下水可用量;和分别表示第k个子区域内第i种作物的产量需求下限和上限。
模型建立模块202,可以根据每一组i,j,k,t取值下的NBijkt、YIit、CFijk、TAt、GWk、和可以建立第一水资源优化模型。
第一水资源优化模型包括两个目标函数,分别为经济目标函数和生态目标函数。第一水资源优化模型中的约束条件至少包括总水量约束、地表水约束、地下水约束、种植面积约束和产量需求约束。
条件转化模块203,对于每一采用地表水灌溉的子区域(假设为第k个子区域),可以根据该子区域地表水量历史数据,获取某一置信水平p下该子区域的地表水可用量,将该子区域的地表水约束中的SW(u)替换为某一置信水平p下该子区域的地表水可用量SW(p),从而将该子区域的地表水量约束条件转化为预设的置信水平p下的确定性约束条件。
模型转化模块204,根据线性加权方法对经济目标函数和生态目标函数进行加权求和,获得综合效益目标函数,通过最大化综合效益,达到获最大化经济目标和最大化生态目标的平衡;对标准化权重进行不确定性处理消除不确定性;根据鲁棒优化方法将对标准化权重进行不确定性处理后的综合效益目标函数进行线性转化,获得第二水资源优化模型;第一水资源模型中的约束条件进行转化后,直接作为第二水资源模型中的约束条件。
方案获取模块205,可以根据预设的违背风险,很方便地获得第二水资源优化模型的一组最优解作为水资源优化方案。
一组最优解为目标区域的每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积。在该组最优解下,目标区域规划时间内的综合效益最大。
本发明实施例提供的水资源优化装置,用于执行本发明上述各实施例提供的水资源优化方法,该水资源优化装置包括的各模块实现相应功能的具体方法和流程详见上述水资源优化方法的实施例,此处不再赘述。
该水资源优化装置用于前述各实施例的水资源优化方法。因此,在前述各实施例中的水资源优化方法中的描述和定义,可以用于本发明实施例中各执行模块的理解。
本发明实施例通过将优化目标的标准化权重扩展为随机分布于名义值为中心,以波动半径为半径的对称区间中来消除目标权重的不确定性,并根据地表水量历史数据消除子区域地表水可用量的不确定性,从而能消除多目标权重和参数的不确定性,能在目标权重和地表水可用量具有不确定性的情况下获得水资源优化方案。进一步地,获取的水资源优化方案具有稳健性,获取水资源优化方案更加简便、耗时更短。
图3为根据本发明实施例提供的电子设备的结构框图。基于上述实施例的内容,如图3所示,该电子设备可以包括:处理器(processor)301、存储器(memory)302和总线303;其中,处理器301和存储器302通过总线303完成相互间的通信;处理器301用于调用存储在存储器302中并可在处理器301上运行的计算机程序指令,以执行上述各方法实施例所提供的水资源优化方法,例如包括:获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,获取每一子区域的种植面积上限和种植面积下限,获取每一子区域内各种作物的产量需求上限和产量需求下限;根据每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,每一子区域的种植面积上限和种植面积下限,以及每一子区域内各种作物的产量需求上限和产量需求下限,建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;其中,第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;第二水资源优化模型以第一水资源优化模型转化后的约束条件为约束条件;根据第二水资源优化模型和预设的违背风险,获取每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得目标区域规划时间内的综合效益最大。
本发明另一实施例公开一种计算机程序产品,计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序,计算机程序包括程序指令,当程序指令被计算机执行时,计算机能够执行上述各方法实施例所提供的水资源优化方法,例如包括:获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,获取每一子区域的种植面积上限和种植面积下限,获取每一子区域内各种作物的产量需求上限和产量需求下限;根据每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,每一子区域的种植面积上限和种植面积下限,以及每一子区域内各种作物的产量需求上限和产量需求下限,建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;其中,第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;第二水资源优化模型以第一水资源优化模型转化后的约束条件为约束条件;根据第二水资源优化模型和预设的违背风险,获取每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得目标区域规划时间内的综合效益最大。
此外,上述的存储器302中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本发明另一实施例提供一种非暂态计算机可读存储介质,非暂态计算机可读存储介质存储计算机指令,计算机指令使计算机执行上述各方法实施例所提供的水资源优化方法,例如包括:获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,获取每一子区域的种植面积上限和种植面积下限,获取每一子区域内各种作物的产量需求上限和产量需求下限;根据每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,目标区域在规划时间内每一单位时间的总可用水量和每一子区域的地下水可用量,每一子区域的种植面积上限和种植面积下限,以及每一子区域内各种作物的产量需求上限和产量需求下限,建立以最大化目标区域规划时间内的经济目标和生态目标的第一水资源优化模型;其中,第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;对于每一采用地表水灌溉的子区域,根据采用地表水灌溉的子区域的地表水量历史数据,将采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;根据线性加权方法和鲁棒优化方法将第一水资源优化模型转化为以最大化目标区域规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以标准化权重的名义值为中心,以预设的标准化权重的波动半径为半径的对称区间中;第二水资源优化模型以第一水资源优化模型转化后的约束条件为约束条件;根据第二水资源优化模型和预设的违背风险,获取每一子区域内在规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得目标区域规划时间内的综合效益最大。
以上所描述的装置实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行上述各个实施例或者实施例的某些部分的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种水资源优化方法,其特征在于,包括:
获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,获取每一所述子区域的种植面积上限和种植面积下限,获取每一所述子区域内各种作物的产量需求上限和产量需求下限;
根据所述每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,每一所述子区域的种植面积上限和种植面积下限,以及每一所述子区域内各种作物的产量需求上限和产量需求下限,建立以最大化所述目标区域所述规划时间内的经济目标和生态目标的第一水资源优化模型;其中,所述第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;
对于每一采用地表水灌溉的子区域,根据所述采用地表水灌溉的子区域的地表水量历史数据,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;
根据线性加权方法和鲁棒优化方法将所述第一水资源优化模型转化为以最大化所述目标区域所述规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以所述标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中;所述第二水资源优化模型以转化后的所述第一水资源优化模型中的约束条件为约束条件;
根据所述第二水资源优化模型和预设的违背风险,获取每一所述子区域内在所述规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得所述目标区域所述规划时间内的综合效益最大。
2.根据权利要求1所述的水资源优化方法,其特征在于,所述根据线性加权方法和鲁棒优化方法将所述第一水资源优化模型转化为以最大化所述目标区域所述规划时间内的综合效益为目标的第二水资源优化模型的具体步骤包括:
根据线性加权方法,将所述第一水资源优化模型中的两个目标函数以预设的经济目标和生态目标的原始权重进行加权,转化为以最大化所述目标区域所述规划时间内的经济目标和生态目标的加权和为目标的过渡目标函数;
将所述经济目标和生态目标的原始权重进行不确定性处理,将所述过渡目标函数转化为单目标函数;
根据鲁棒优化方法将所述单目标函数线性化,获得所述第二水资源优化模型中的目标函数,并将转化后的所述第一水资源优化模型中的约束条件作为所述第二水资源优化模型中的约束条件。
3.根据权利要求2所述的水资源优化方法,其特征在于,所述将所述经济目标和生态目标的原始权重进行不确定性处理的具体步骤包括:
将所述经济目标的原始权重和所述生态目标的原始权重,分别除以所述经济目标的原始权重,获得所述生态目标的标准化权重的名义值,将所述生态目标的标准化权重随机分布于以所述生态目标的标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中;
或者,将所述经济目标的原始权重和所述生态目标的原始权重,分别除以所述生态目标的原始权重,获得所述经济目标的标准化权重的名义值,将所述经济目标的标准化权重随机分布于以所述经济目标的标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中。
4.根据权利要求1所述的水资源优化方法,其特征在于,所述第二水资源优化模型中的目标函数为
max z
s.t.
q≥0
其中,z表示目标区域规划时间内的综合效益;i表示作物种类的编号;j表示灌溉水源的编号;k表示子区域的编号;t表示规划时间内单位时间的编号;I表示作物种类的数量;J表示灌溉水源的数量;K表示子区域的数量;T表示规划时间包括的单位时间的数量;Aijkt表示第k个子区域内在第t个单位时间用第j种灌溉水源灌溉第i种作物的灌溉面积;表示Aijkt的经济收益;表示Aijkt的生态收益;Γ表示保护度水平;q表示第一辅助变量;rijkt表示Aijkt对应的第二辅助变量;W表示标准化权重的名义值;表示标准化权重的波动半径。
5.根据权利要求4所述的水资源优化方法,其特征在于,所述保护度水平根据所述预设的违背风险获取。
6.根据权利要求5所述的水资源优化方法,其特征在于,根据所述预设的违背风险获取所述保护度水平的公式为
其中,PBU表示预设的违背风险;v=(Γ+n)/2;
n=I×J×K×T。
7.根据权利要求1至6任一所述的水资源优化方法,其特征在于,所述根据所述采用地表水灌溉的子区域的地表水量历史数据,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件的具体步骤包括:
根据所述采用地表水灌溉的子区域的地表水量历史数据,获取所述采用地表水灌溉的子区域的地表水量的分布模型;
根据分布模型,获取所述预设的置信水平下所述采用地表水灌溉的子区域单位时间的地表水可用量;
根据所述预设的置信水平下所述采用地表水灌溉的子区域单位时间的地表水可用量,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件。
8.一种水资源优化装置,其特征在于,包括:
数据获取模块,用于获取目标区域的每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,获取所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,获取每一所述子区域的种植面积上限和种植面积下限,获取每一所述子区域内各种作物的产量需求上限和产量需求下限;
模型建立模块,用于根据所述每一子区域内各种作物的单位经济净效益、单位产量和单位生态净效益,所述目标区域在规划时间内每一单位时间的总可用水量和每一所述子区域的地下水可用量,每一所述子区域的种植面积上限和种植面积下限,以及每一所述子区域内各种作物的产量需求上限和产量需求下限,建立以最大化所述目标区域所述规划时间内的经济目标和生态目标的第一水资源优化模型;其中,所述第一水资源优化模型以总水量、各种灌溉水源的水量、种植面积和产量需求为约束条件;
条件转化模块,用于对于每一采用地表水灌溉的子区域,根据所述采用地表水灌溉的子区域的地表水量历史数据,将所述采用地表水灌溉的子区域的地表水量约束条件转化为预设的置信水平下的确定性约束条件;
模型转化模块,用于根据线性加权方法和鲁棒优化方法将所述第一水资源优化模型转化为以最大化所述目标区域所述规划时间内的综合效益为目标的第二水资源优化模型;其中,线性加权方法中的标准化权重,随机分布于以所述标准化权重的名义值为中心,以预设的所述标准化权重的波动半径为半径的对称区间中;所述第二水资源优化模型以转化后的所述第一水资源优化模型中的约束条件为约束条件;
方案获取模块,用于根据所述第二水资源优化模型和预设的违背风险,获取每一所述子区域内在所述规划时间的每一单位时间内用每一种灌溉水源灌溉每种作物的灌溉面积,以使得所述目标区域所述规划时间内的综合效益最大。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7任一项所述的水资源优化方法的步骤。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至7任一项所述的水资源优化方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910328856.4A CN110135627B (zh) | 2019-04-23 | 2019-04-23 | 水资源优化方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910328856.4A CN110135627B (zh) | 2019-04-23 | 2019-04-23 | 水资源优化方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110135627A true CN110135627A (zh) | 2019-08-16 |
CN110135627B CN110135627B (zh) | 2021-07-16 |
Family
ID=67570823
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910328856.4A Active CN110135627B (zh) | 2019-04-23 | 2019-04-23 | 水资源优化方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110135627B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112215432A (zh) * | 2020-10-21 | 2021-01-12 | 河南大学 | 一种水资源可达性预测方法和预测*** |
CN112465206A (zh) * | 2020-11-20 | 2021-03-09 | 广东工业大学 | 一种水资源优化配置及优化结果评价反馈方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020111758A1 (en) * | 2000-10-18 | 2002-08-15 | Qing-Guo Wang | Robust process identification and auto-tuning control |
US20080022114A1 (en) * | 1996-07-02 | 2008-01-24 | Wistaria Trading, Inc. | Optimization methods for the insertion, protection, and detection of digital watermarks in digital data |
CN104657591A (zh) * | 2015-01-19 | 2015-05-27 | 华北电力大学 | 水资源配置方法和装置 |
CN105787501A (zh) * | 2015-12-17 | 2016-07-20 | 武汉大学 | 输电线路走廊区域自动选择特征的植被分类方法 |
CN105956718A (zh) * | 2016-05-31 | 2016-09-21 | 中国农业大学 | 一种灌区水资源优化配置方法和*** |
-
2019
- 2019-04-23 CN CN201910328856.4A patent/CN110135627B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080022114A1 (en) * | 1996-07-02 | 2008-01-24 | Wistaria Trading, Inc. | Optimization methods for the insertion, protection, and detection of digital watermarks in digital data |
US20020111758A1 (en) * | 2000-10-18 | 2002-08-15 | Qing-Guo Wang | Robust process identification and auto-tuning control |
CN104657591A (zh) * | 2015-01-19 | 2015-05-27 | 华北电力大学 | 水资源配置方法和装置 |
CN105787501A (zh) * | 2015-12-17 | 2016-07-20 | 武汉大学 | 输电线路走廊区域自动选择特征的植被分类方法 |
CN105956718A (zh) * | 2016-05-31 | 2016-09-21 | 中国农业大学 | 一种灌区水资源优化配置方法和*** |
Non-Patent Citations (3)
Title |
---|
QIAN TAN等: "Robust fractional programming approach for improving agricultural wateruse efficiency under uncertainty", 《JOURNAL OF HYDROLOGY》 * |
R. TIMOTHY MARLER等: "The weighted sum method for multi-objective optimization: new insights", 《STRUCTURAL MULTIDISCIPLINARY OPTIMIZATION》 * |
张珊等: "基于模糊优选和可信性的农业水资源多目标优化配置模型", 《南水北调与水利科技》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112215432A (zh) * | 2020-10-21 | 2021-01-12 | 河南大学 | 一种水资源可达性预测方法和预测*** |
CN112465206A (zh) * | 2020-11-20 | 2021-03-09 | 广东工业大学 | 一种水资源优化配置及优化结果评价反馈方法 |
CN112465206B (zh) * | 2020-11-20 | 2023-01-20 | 广东工业大学 | 一种水资源优化配置及优化结果评价反馈方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110135627B (zh) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Montazar et al. | Conjunctive water use planning in an irrigation command area | |
Bekchanov et al. | Pros and cons of adopting water-wise approaches in the lower reaches of the Amu Darya: a socio-economic view | |
Loch et al. | The role of water markets in climate change adaptation | |
Dubovyk et al. | Land suitability assessment for afforestation with Elaeagnus angustifolia L. In degraded agricultural areas of the lower amudarya river basin | |
CN106779255B (zh) | 流域农业面源污染控制的最佳管理措施优化方法及装置 | |
Reddy et al. | Evolving strategies for crop planning and operation of irrigation reservoir system using multi-objective differential evolution | |
Schindler | A multi-agent system for simulating land-use and land-cover change in the Atankwidi catchment of Upper East Ghana | |
CN108022023A (zh) | 一种农作物的种植产量预测方法、装置及存储介质 | |
McConnell et al. | Precision conservation: A geospatial decision support tool for optimizing conservation and profitability in agricultural landscapes | |
CN104732426A (zh) | 一种农产品产销决策方法、装置及*** | |
CN113785676A (zh) | 一种盐碱地土壤改良体系及其构建方法 | |
Safavi et al. | Conjunctive use of surface water and groundwater resources under deficit irrigation | |
CN110135627A (zh) | 水资源优化方法及装置 | |
Niles et al. | Farmers share their perspectives on California water management and the Sustainable Groundwater Management Act | |
Bjornlund et al. | Institutional innovation and smart water management technologies in small-scale irrigation schemes in southern Africa | |
CN116307191A (zh) | 一种基于人工智能算法的水资源配置方法、装置和设备 | |
Sonneveld et al. | The impact of climate change on crop production in west Africa: An assessment for the oueme river basin in Benin | |
CN109767046A (zh) | 一种土地空间优化配置方法及*** | |
Yu et al. | Model based decision support system for land use changes and socio-economic assessments | |
Ciancio et al. | Systemic forest management and operational perspectives for implementing forest conservation in Italy under a pan-European framework | |
Knapp et al. | ECONOMIC POLICIES FOR REGULATING AGRICULTURAL DRAINAGE WATER 1 | |
Malano et al. | Ranking and classification of irrigation system performance using fuzzy set theory: case studies in Australia and China | |
Flichman et al. | Economic policy and water pollution | |
NA et al. | Optimal regeneration regime under continuous crown cover requirements in cork oak woodlands | |
Lacerda et al. | Multi-objective forest harvesting under sustainable and economic principles |
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 |