CN106529166B - 一种基于maepso算法的区域水资源优化配置方法 - Google Patents
一种基于maepso算法的区域水资源优化配置方法 Download PDFInfo
- Publication number
- CN106529166B CN106529166B CN201610973576.5A CN201610973576A CN106529166B CN 106529166 B CN106529166 B CN 106529166B CN 201610973576 A CN201610973576 A CN 201610973576A CN 106529166 B CN106529166 B CN 106529166B
- Authority
- CN
- China
- Prior art keywords
- particle
- value
- water
- optimal
- algorithm
- 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.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000002245 particle Substances 0.000 claims description 108
- 238000005457 optimization Methods 0.000 claims description 39
- 230000008901 benefit Effects 0.000 claims description 25
- 230000006870 function Effects 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000005086 pumping Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000002068 genetic effect Effects 0.000 description 3
- 238000013468 resource allocation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000001363 water suppression through gradient tailored excitation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于MAEPSO算法的区域水资源优化配置方法，步骤如下：首先获取水资源系统的基本信息数据；其次建立水资源多目标配置模型；然后执行基于MAEPSO算法，求解区域水资源系统多目标配置的Pareto最优解集；最后按照一定的规则从最优解集中选取最终解。本发明实现全局寻优，提高计算效率，满足选择水资源系统多目标最优配置方案要求。
Description
技术领域
本发明属于水利水电领域的水资源配置技术领域，特别是一种基于MAEPSO算法的区域水资源优化配置方法。
背景技术
区域水资源优化配置是一个高度复杂的多层次、多阶段、多主体、多目标、多关联、非线性的风险决策问题。近年来，越来越多的学者对水资源优化配置的模型构建和求解问题进行了深入研究，其中，智能优化算法已在模型求解中得到广泛应用，如遗传算法、模糊退化、人工神经网络、混沌优化算法、蚁群算法等进化算法逐渐广泛地应用于水资源优化配置中。但是由于水资源优化配置的复杂性和多目标性，常规优化算法在收敛性、计算速度、初值敏感性等方面效果不是很理想。为实现总体综合效益的最大化，需要对各目标进行权衡择优，很多研究者开始致力于多目标进化算法(MOEA)的研究。近年来有关MOEA的研究大多是在群体进化过程中引入非支配解集(Pareto)思想，为水资源的多目标优化配置提供了坚实的基础。
粒子群优化算法(PSO)是1995年Kennedy和Eberhart提出的基于群体智能的全局随机搜索算法，源于鸟群觅食过程中的迁徙和群聚行为的研究。该方法基于群体迭代，同遗传算法类似，但是没有遗传算法的交叉、变异操作，在单目标优化、约束优化、动态优化、多目标优化以及动态多目标优化等问题上得到广泛的应用。现实生活中需要对多个目标同时进行优化，单目标粒子群算法无法直接用于多目标问题的求解，需要修正单目标粒子群算法以便能够生成多个解，维持找到的多个解并且保持其多样性。多目标粒子群(MOPSO)算法因其概念简单、容易实现以及需要调整的参数较少等优点，很多学者加入到该算法的研究中来，使其成为解决多目标优化问题的有力工具。
近年来，国内外不少学者将MOPSO算法应用到区域水资源优化配置领域并得到了广泛的应用，现阶段提出的基于MOPSO算法的区域水资源优化配置存在易陷于局部最优解、非劣解分布不均等缺陷。同时，区域水资源优化配置领域在采用多目标决策方法处理非劣解集时，需要基于主客观因素相结合确定评价指标权重。
发明内容
发明目的：针对基本MOPSO算法容易陷入局部收敛，外部档案集非劣解分布不均，个体难收敛到非劣解前沿等缺陷，提供一种基于MAEPSO算法(Multi-objective AdaptiveExpanded Particle Swarm Optimization，多目标自适应扩展粒子群算法)的区域水资源优化配置方法，并在外部档案非劣解集的基础上采用主客观权重相结合的多目标决策理论选择最优水资源配置方案，实现区域水资源多目标优化配置。
技术方案：一种基于MAEPSO算法的区域水资源优化配置方法，包括以下步骤：
步骤一：获取区域水资源系统工程的基本信息数据，包括：泵、闸站以及水库水利工程过流能力值q，初末库容限制V，正常蓄水位Z正、防洪限制水位Z防、死水位Z死，容积-水位关系曲线S～Z，下游水位-下泄流量关系曲线Z～Q，发电机组出力约束值N，来水量W；
步骤二：建立以社会效益、经济效益以及生态环境等的综合效益最大的目标函数，考虑水量平衡、机组出力、过流能力等约束条件的多目标水资源优化配置数学模型：
opt F(x)＝{f1(x),f2(x),…,fn(x)} (1)
s.t x∈G(x) (2)
式中，opt表示优化方向，包括最大方向和最小方向，n表示水资源系统优化调度的目标数；F(x)目标函数集；fn(x)表示为社会效益、经济效益以及生态环境等的综合效益最大的目标函数；G(x)表示约束条件集，包括水量平衡约束、水源可供水量约束、水源输水能力约束、用户需水能力约束、变量非负约束。
步骤三：执行改进的MAEPSO算法；
步骤四：基于外部归档集EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案。
进一步的，所述步骤三：执行改进的MAEPSO算法具体为：
步骤3-1：确定以水资源工程时段末蓄水量S、下泄流量u或工程过水流量q为决策变量，划分水资源调度周期时段T；
步骤3-2：设定参数；确定由决策变量组成的初始种群规模G，搜索空间为D维，全局迭代次数MAXGEN，学习因子c1和c2，惯性因子w，粒子速度[vmin,vmax]，其中，vmin为离子速度最小值，vmax为离子速度最大值，定义第j个粒子位置表示为向量Xj＝(xj1,xj2,…,xjD)，速度表示为向量Vj＝(vj1,vj2,…,vjD)，历史最优位置表示为Qj＝(qj1,qj2,…,qjD)，外部归档集规模NEA；
步骤3-3，生成初始种群，包括初始位置向量与初始速度向量；初始全局迭代次数，GEN＝0，外部归档集EA赋空集；
步骤3-4，计算每个个体的各个目标函数值，设置初始种群的个体极值及全局极值gbest(Pg＝min(P1 1,P2 1,…,PG 1))，其中，Pj 1表示第1次迭代第j个粒子极值；将其非劣解加入到外部档案EA中；
步骤3-5：分别利用和更新粒子当前的速度和位置，并且保证
其中，c1和c2为学习因子，c1是用来调节粒子飞向其最优位置的步长，c2是用来调节粒子飞向整个群体最优位置的步长；r1和r2为[0，1]范围之间的随机数；w为权重；表示第t次迭代第j个粒子第n决策变量的飞行速度；表示第t次迭代第j个粒子第n决策变量位置；表示第j个粒子第n决策变量的历史最优位置；qgn表示第n决策变量的历史全局最优位置；
惯性因子采用依据粒子群群体适应度值变化而进行调整的自适应策略(self-adaptive strategy)如下：
式中，fi表示粒子i的适应度值，favg1表示适应度值大于种群粒子平均适应度值的全部粒子的平均适应度值，favg2表示适应度值小于种群粒子平均适应度值的全部粒子的平均适应度值；采用线性调整学习因子策略，即c1先大后小，c2先小后大的方法：
式中，c1s,c1e为学习因子c1的初末值，分别取2.5、0.5；c2s,c2e为学习因子c2的初末值，分别取0.5、2.5；t为目前的迭代次数，tmax是最大迭代次数；
步骤3-6：计算当前个体目标函数值，根据更新当前个体至最优pbest，同时根据更新群体至最优gbest，其中t表示第t代循环，Xj t、Xj t+1分别表示第t、t+1次迭代第j个粒子，f(Xj t)、f(Xj t+1)表示第t、t+1次迭代第j个粒子适应性函数值；
步骤3-7：根据新的非劣解采用基于距离的方法维护外部档案NEA：每次迭代都计算每个非劣解距离左右两个粒子的距离，通过比较该非劣解的两个距离是否小于某个限值Dis，如果均小于限值Dis，则从外部档案中删除此非劣解，并重新计算左右两边粒子的距离；
步骤3-8：基于Sigma选取策略来选择全局向导，为各粒子选择新的pbest，以双目标问题为例，f1和f2分别代表两个目标；
步骤3-9：判断是否达到全局迭代次数MAXGEN，若没有达到，
GEN＝GEN+1，转到步骤25，继续下一轮的全局搜索，否则，算法结束，输出EA。
进一步的，所述步骤3-7根据新的非劣解采用基于距离的方法维护外部档案NEA：每次迭代都计算每个非劣解距离左右两个粒子的距离，通过比较该非劣解的两个距离是否小于某个限值Dis，如果均小于限值Dis，则从外部档案中删除此非劣解，并重新计算左右两边粒子的距离具体为：
步骤3-7-1，将找到的所有非劣解集按照某个目标适应度值从小到大进行排序；
步骤3-7-2，计算非劣解与非劣解在目标空间的欧氏距离：
其中，n为目标的个数；
步骤3-7-3，判断上一步中计算的每个距离是否小于限值Dis，如果disi-1＜Dis且disi＜Dis，则删除非劣解i并重新计算非劣解i-1与非劣解i+1的欧氏距离作为disi-1；临界距离Dis设计为根据非劣解的情况自适应调整：
步骤3-7-4，判断非劣解的个数是否超过外部档案的大小，若超过转S2，否则转S5；
步骤3-7-5，返回外部归档集EA。
进一步的，所述步骤3-8，基于Sigma选取策略来选择全局向导，为各粒子选择新的pbest，以双目标问题为例，f1和f2分别代表两个目标，具体为：
步骤3-8-1，计算外部归档集EA中所有成员的Sigma值σa，如下式计算：
其中:σa就是Sigma值，f1,f2为对应的目标函数值，a＝1,2,…,NEA；
步骤3-8-2，从种群中选择一个粒子，计算该粒子的Sigma值σi；
步骤3-8-3，计算该粒子的Sigma值σi与外部档案所有成员的Sigma值之间的欧氏距离，选择与该粒子Sigma值欧氏距离最小的外部档案成员作为该粒子的pbest，更新粒子当前的速度和位置。
进一步的，所述步骤四基于外部归档集EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案，包括以下步骤：
步骤4-1，利用模糊层次分析法确定属性主观权重W'；
步骤4-2，利用CRITIC法确定属性客观权重W”；
步骤4-3，确定属性组合权重W：当得到各评价指标的主观权重wi'和客观权重wj”后，按式计算各指标组合权重wj:
步骤4-4，基于TOPSIS法方案确定最终优化配置方案。
进一步的，所述步骤4-1利用模糊层次分析法确定属性主观权重W'具体如下：
步骤4-1-1，根据所选取的各个目标函数，建立系统递阶层次结构，并进行各层次之间指标间的两两比较，采用0.1～0.9标度法进行定量描述，建立模糊一致判断矩阵R：
其中rij表示元素ai比aj重要的隶属度，rij越大，ai比aj越重要；n为指标数目；
步骤4-1-2，主观权重wj'计算
其中，参数a取值大小的选择反映了决策者的个人偏好，a越小则表明决策者偏好权重之间差异程度越大；k＝1,2,…,n。
进一步的，所述步骤4-2利用CRITIC法确定属性客观权重W”具体如下：
步骤4-2-1，计算指标j对应的信息量Cj，指标i与其他指标之间的冲突性量化指标而用式表示：
式中，δj为指标j的标准差，pjk为指标j,k的相关系数；
步骤4-2-2，所以指标j对应的归一化的权重wj”为：
进一步的，所述步骤4-4基于TOPSIS法方案确定最终优化配置方案具体如下：
步骤4-4-1，设多目标水资源优化配置问题的决策矩阵A为
式中，fij表示第i个可行方案第j个目标值，i＝1,2,…NEA,j＝1,2,…,n；
由该矩阵A构成规范决策矩阵Z'，其元素为
步骤4-4-2，构造加权的规范决策矩阵Zij，其中的元素zij为zij＝wjrij；
步骤4-4-3，在加权规范化决策矩阵Zij中选取各个属性都达到最优的方案为正理想方案Z+，各个属性都达到最劣的方案为负理想方案Z-；
其中，n为效益型指标数目，n”为成本型指标数目，n＝n'+n”；
步骤4-4-4，采用欧几里得范数来计算各方案与正、负理想方案距离，不同方案到正、负理想方案的距离分别为：
式中，分别为第i个方案与正、负理想方案的欧几里得范数距离；zij为第i个方案第j个指标的规范化值；分别为第j个指标的正、负理想方案属性值；
步骤4-4-5，计算各个方案到正、负理想方案的相对接近度
步骤4-4-6，根据相对接近度指数的大小对样本方案进行优选排序；当方案越靠近正理想方案时，越接近1，取最大值对应的配置方案作为最优方案。
有益效果：本发明提出的一种基于MAEPSO算法的区域水资源优化配置方法，首先获取水资源系统的基本信息数据；其次建立水资源多目标配置模型；然后执行基于MAEPSO算法，求解区域水资源系统多目标配置的Pareto最优解集；最后按照一定的规则从最优解集中选取最终解。本发明实现全局寻优，提高计算效率，满足选择水资源系统多目标最优配置方案要求。相对于现有技术，本发明具有以下优势：
(1)满足区域水资源系统多目标优化配置的要求；
(2)惯性因子采用依据粒子群群体适应度值变化而进行调整的自适应策略，更好的实现全局搜索与局部搜索之间的平衡，防止陷入局部收敛；
(3)采用线性调整学习因子策略，实现搜索初期粒子飞行主要参考粒子本身的历史信息，后期注重社会信息的功能，改善算法收敛性；
(4)采用基于距离的方法对多目标优化算法的外部档案进行维护，确保非劣解个体分布均匀，具有良好的多样性，加快全局收敛。
(5)基于Sigma选取策略来选择全局向导，使得种群里的各个粒子被各自离得近的非支配解导向，能够促使种群个体尽快收敛到非劣解前沿；
(6)所采用的多目标决策方法，主客观权重相结合，增强了水资源系统多目标配置系统的自主决策能力，又避免决策过程中过度带入主观偏好。
附图说明
图1是本发明的方法流程图；
图2是本发明中MAEPSO算法的流程图；
图3是本发明中多目标决策方法流程图；
图4是本发明中平水年南水北调东线工程湖泊群EA解集空间分布图。
具体实施方式
下面结合具体实施例，进一步阐明本发明。
本发明针对基本MOPSO算法容易陷入局部收敛，外部档案集非劣解分布不均，个体难收敛到非劣解前沿等缺陷，提供一种基于MAEPSO算法的区域水资源优化配置方法。该方法采用惯性因子以及学习因子自适应调整策略，利用基于距离的方法对外部档案进行维护，基于Sigma选取策略选择全局向导，并在非劣解集的基础上采用主客观权重相结合的多目标决策理论选择最优水资源配置方案，实现区域水资源多目标优化配置。
如图1-3所示，一种基于MAEPSO算法的区域水资源优化配置方法，包括以下步骤：
步骤一：获取区域水资源系统工程的基本信息数据，包括：泵、闸站以及水库等水利工程过流能力值q，初末库容限制V，正常蓄水位Z正、防洪限制水位Z防、死水位Z死，容积-水位关系曲线S～Z，下游水位-下泄流量关系曲线Z～Q，发电机组出力约束值N，来水量W；
步骤二：建立以社会效益、经济效益以及生态环境等的综合效益最大等目标函数，考虑水量平衡、机组出力、过流能力等约束条件的多目标水资源优化配置数学模型；
minF(x)＝{f1(x),f2(x),…,fn(x)} (1)
式中，n—表示水资源系统优化调度的目标数，n＝1,2,…,N；F(x)—目标函数集；fn(x)—表示为社会效益、经济效益以及生态环境等的综合效益最大的目标函数；Δt—表示计算时段区间；Wt,ut—表示水库Δt时段内水库、泵站、水闸等单元入流量、下泄流量，m3/s；St,St+1—表示水库、泵站、水闸等单元t时段、t+1时段末蓄水量，m3；It—表示Δt时段内的损失量，m3；St,min,St,max—表示t时段末允许的最低库容、最高库容，m3；qt—表示t时段内的发电流量，m3/s；qt,min,qt,max—表示t时段末允许的最小下泄流量、最大下泄流量，m3/s；Nt—表示水库Δt时段内水电站的出力，kW；Nt,min,Nt,max表示水电站t时段末允许的最小出力值、最大出力值，kW；S,S'—表示水库初末库容限制值，m3。
步骤三：执行MAEPSO算法；
步骤四：基于外部归档集EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案。
考虑多粒子影响和自适应控制参数(惯性因子以及学习因子)的多目标自适应扩展粒子群优化算法，步骤三：执行MAEPSO算法，包括以下步骤：
步骤21，确定以水资源工程时段末蓄水量(末水位)S、下泄流量u或工程过水流量q为决策变量，划分水资源调度周期时段T；
步骤22，设定参数。确定由决策变量组成的初始种群规模G，搜索空间为D维，全局迭代次数MAXGEN，学习因子c1和c2，惯性因子w，粒子速度[vmin,vmax]，其中，vmin为离子速度最小值，vmax为离子速度最大值，定义第j个粒子位置表示为向量Xj＝(xj1,xj2,…,xjD)，速度表示为向量Vj＝(vj1,vj2,…,vjD)，历史最优位置表示为Qj＝(qj1,qj2,…,qjD)，外部归档集规模NEA；
步骤23，生成初始种群，包括初始位置向量与初始速度向量；初始全局迭代次数，GEN＝0，外部归档集(EA)赋空集；
步骤24，计算每个个体的各个目标函数值，设置初始种群的个体极值及全局极值gbest(Pg＝min(P1 1,P2 1,…,PG 1))，其中，Pj 1表示第1次迭代第j个粒子极值；将其非劣解加入到外部档案(EA)中。
步骤25，分别利用和更新粒子当前的速度和位置，并且保证
其中，c1和c2为学习因子，c1是用来调节粒子飞向其最优位置的步长，c2是用来调节粒子飞向整个群体最优位置的步长；r1和r2为[0，1]范围之间的随机数；w为权重；表示第t次迭代第j个粒子第n决策变量的飞行速度；表示第t次迭代第j个粒子第n决策变量位置；表示第j个粒子第n决策变量的历史最优位置；qgn表示第n决策变量的历史全局最优位置；
惯性因子采用依据粒子群群体适应度值变化而进行调整的自适应策略，如下：
式中，fi表示粒子i的适应度值，favg1表示适应度值大于种群粒子平均适应度值的全部粒子的平均适应度值，favg2表示适应度值小于种群粒子平均适应度值的全部粒子的平均适应度值。
本发明采用线性调整学习因子策略，即c1先大后小，c2先小后大的思想：
式中，c1s,c1e为学习因子c1的初末值，分别取2.5、0.5；c2s,c2e为学习因子c2的初末值，分别取0.5、2.5；t为目前的迭代次数，tmax是最大迭代次数。
步骤26，计算当前个体目标函数值，根据更新当前个体至最优pbest，同时根据更新群体至最优gbest，其中t表示第t代循环，Xj t、Xj t+1分别表示第t、t+1次迭代第j个粒子，f(Xj t)、f(Xj t+1)表示第t、t+1次迭代第j个粒子适应性函数值；
步骤27，根据新的非劣解采用基于距离的方法维护外部档案EA：每次迭代都计算每个非劣解距离左右两个粒子的距离，通过比较该非劣解的两个距离是否小于某个限值Dis，如果均小于限值Dis，则从外部档案中删除此非劣解，并重新计算左右两边粒子的距离，步骤27具体包括：
步骤S1，将找到的所有非劣解集按照某个目标适应度值从小到大进行排序。
步骤S2，计算非劣解与非劣解在目标空间的欧氏距离：
其中，n为目标的个数。
步骤S3，判断上一步中计算的每个距离是否小于限值Dis，如果disi-1＜Dis且disi＜Dis，则删除非劣解i并重新计算非劣解i-1与非劣解i+1的欧氏距离作为disi-1。临界距离Dis设计为根据非劣解的情况自适应调整：
步骤S4，判断非劣解的个数是否超过外部档案的大小，若超过转S2，否则转S5。
步骤S5，返回外部档案。
步骤28，基于Sigma选取策略来选择全局向导，为各粒子选择新的pbest，以双目标问题为例，f1和f2分别代表两个目标，步骤28具体包括：
步骤S1，计算档案EA中所有成员的Sigma值σa，如下式计算：
其中:σa就是Sigma值，f1,f2为对应的目标函数值，a＝1,2,…,NEA。
步骤S2，从种群中选择一个粒子，计算该粒子的Sigma值σi；
步骤S3，计算该粒子的Sigma值σi与外部档案所有成员的Sigma值之间的欧氏距离，选择与该粒子Sigma值欧氏距离最小的外部档案成员作为该粒子的pbest，更新粒子当前的速度和位置。
步骤29，判断是否达到全局迭代次数MAXGEN，若没有达到，GEN＝GEN+1，转到步骤25，继续下一轮的全局搜索，否则，算法结束，输出EA。
基于EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案，包括以下步骤：
步骤31，利用模糊层次分析法确定属性主观权重W'，具体步骤如下：
步骤S1，根据所选取的各个目标函数，建立系统递阶层次结构，并进行各层次之间指标间的两两比较，采用0.1～0.9标度法进行定量描述，建立模糊一致判断矩阵R：
其中rij表示元素ai比aj重要的隶属度，rij越大，ai比aj越重要；n为指标数目；
步骤S2，主观权重wj'计算
其中，参数a取值大小的选择反映了决策者的个人偏好，a越小则表明决策者偏好权重之间差异程度越大；k＝1,2,…,n；
步骤32，利用CRITIC法确定属性客观权重W”，具体如下：
步骤S1，计算指标j对应的信息量Ci，指标i与其他指标之间的冲突性量化指标而可用式表示：
式中，δj为指标j的标准差，pjk为指标j,k的相关系数。
步骤S2，所以指标j对应的归一化的权重wj”为：
步骤33，确定属性组合权重W：当得到各评价指标的主观权重wi'和客观权重wj”后，可按式计算各指标组合权重wj:
步骤34，基于TOPSIS法方案确定最终优化配置方案，具体步骤如下：
步骤S1，设多目标水资源优化配置问题的决策矩阵A为
式中，fij表示第i个可行方案第j个目标值，i＝1,2,…NEA,j＝1,2,…,n。
由该矩阵A构成规范决策矩阵Z'，其元素为
Step2构造加权的规范决策矩阵Zij，其中的元素zij为zij＝wjrij；
步骤S3，在加权规范化决策矩阵Zij中选取各个属性都达到最优的方案为正理想方案Z+，各个属性都达到最劣的方案为负理想方案Z-。
其中，n为效益型指标数目，n”为成本型指标数目，n＝n'+n”。
步骤S4，采用欧几里得范数来计算各方案与正、负理想方案距离，不同方案到正、负理想方案的距离分别为：
式中，分别为第i个方案与正、负理想方案的欧几里得范数距离；zij为第i个方案第j个指标的规范化值；分别为第j个指标的正、负理想方案属性值。
步骤S5，计算各个方案到正、负理想方案的相对接近度
步骤S6，根据相对接近度指数的大小对样本方案进行优选排序。当方案越靠近正理想方案时，越接近1，取最大值对应的配置方案作为最优方案。
现以南水北调东线工程江苏段水资源多目标优化配置为例，说明发明方法的有效性与合理性。南水北调东线工程江苏段工程沿线分布有众多天然湖泊，自江都站起依次连接洪泽湖、骆马湖、南四湖，总调蓄库容45.25亿m3，邻湖泊间的水位差都在l0m左右，若以湖泊为节点，则从长江至东平湖下共可分为三个大段，每段设3级提水泵站，共计9个提水梯级，可划分为16个计算受水区。
本发明将以平水年泵站抽水量设为决策变量，采用改进的粒子群算法进行优化调度，实现缺水量最小、抽江水量最小最大两目标。经反复测试计算，确定改进粒子群算法的具体参数设置为：粒子群种群规模设置为100，迭代次数为1000，外部档案Np设置为100，wstart＝0.9，wend＝0.2，c1s＝2.5，c1e＝0.5，c2s＝0.5，c2e＝2.5，调度方案集空间分布如图4所示，调度方案集在空间分布呈现为一条非凸曲线，调度方案分布广泛且均匀，抽江水量最小和缺水量最小两目标之间相互制约、相互冲突，存在明显的反比关系，因此改进的粒子群求解的水资源配置方案集是合理有效的。基于Pareto最优解，采用多目标决策方法确定水库最优配置方案，优化调度方案多属性模糊一致矩阵及权重见表1。在此基础上，选择使受水区缺水量目标值最小的前10个方案作为待选方案，经过方案优选，所得调度方案所对应的受水区缺水量为0.177亿m3，抽江水量为8.059亿m3，需水满足度达到99.9％，湖泊汛末库容变化量为3.303亿m3。输水线路运河线、运西线总抽水量分别为160.232亿m3、47.351亿m3，具体年调水情况见表2。系统全年总供水量为125.685亿m3，其中高峰期(6月～9月)供水62.084亿m3，非高峰期(10月～次年5月)供水63.601亿m3，各受水区全年实际供水情况见表3。
表1优化调度方案多属性模糊一致矩阵及权重
表2各线路具体年调水结果 单位：亿m3
表3各受水区全年实际供水情况 单位：亿m3
Claims (5)
1.一种基于MAEPSO算法的区域水资源优化配置方法，其特征在于，所述MAEPSO算法为多目标自适应扩展粒子群算法，所述方法包括以下步骤：
步骤一：获取区域水资源系统工程的基本信息数据，包括：泵、闸站以及水库水利工程过流能力值q，初末库容限制V，正常蓄水位Z正、防洪限制水位Z防、死水位Z死，容积-水位关系曲线S～Z，下游水位-下泄流量关系曲线Z～Q，发电机组出力约束值N，来水量W；
步骤二：建立以社会效益、经济效益以及生态环境的综合效益最大的目标函数，考虑水量平衡、机组出力、过流能力约束条件的多目标水资源优化配置数学模型：
opt F(x)＝{f1(x),f2(x),…,fn(x)} (1)
其中x∈G(x) (2)
式中，opt表示优化方向，包括最大方向和最小方向，n表示水资源系统优化调度的目标数；F(x)目标函数集；fn(x)表示为社会效益、经济效益以及生态环境的综合效益最大的目标函数；G(x)表示约束条件集，包括水量平衡约束、水源可供水量约束、水源输水能力约束、用户需水能力约束、变量非负约束；
步骤三：执行改进的MAEPSO算法，具体包括以下步骤：
步骤3-1：确定以水资源工程时段末蓄水量S、下泄流量u或工程过水流量q为决策变量，划分水资源调度周期时段T；
步骤3-2：设定参数；确定由决策变量组成的初始种群规模G，搜索空间为D维，全局迭代次数MAXGEN，学习因子c1和c2，惯性因子w，粒子速度[vmin,vmax]，其中，vmin为离子速度最小值，vmax为离子速度最大值，定义第j个粒子位置表示为向量Xj＝(xj1,xj2,…,xjD)，速度表示为向量Vj＝(vj1,vj2,…,vjD)，历史最优位置表示为Qj＝(qj1,qj2,…,qjD)，外部归档集规模NEA；
步骤3-3，生成初始种群，包括初始位置向量与初始速度向量；初始全局迭代次数，GEN＝0，外部归档集EA赋空集；
步骤3-4，计算每个个体的各个目标函数值，设置初始种群的个体极值及全局极值gbest(Pg＝min(P1 1,P2 1,…,PG 1))，其中，Pj 1表示第1次迭代第j个粒子极值；将其非劣解加入到外部档案EA中；
步骤3-5：分别利用和更新粒子当前的速度和位置，并且保证
其中，c1和c2为学习因子，c1是用来调节粒子飞向其最优位置的步长，c2是用来调节粒子飞向整个群体最优位置的步长；r1和r2为[0，1]范围之间的随机数；w为权重；表示第t次迭代第j个粒子第n决策变量的飞行速度；表示第t次迭代第j个粒子第n决策变量位置；表示第j个粒子第n决策变量的历史最优位置；xn,min为允许的最小值；xn,max为允许的最大值；qgn表示第n维决策变量的历史全局最优位置；
惯性因子采用依据粒子群群体适应度值变化而进行调整的自适应策略如下：
式中，wstart表示惯性因子的最大值，wend表示惯性因子最小值，fi表示粒子i的适应度值，favg1表示适应度值大于种群粒子平均适应度值的全部粒子的平均适应度值，favg2表示适应度值小于种群粒子平均适应度值的全部粒子的平均适应度值；采用线性调整学习因子策略，即c1先大后小，c2先小后大的方法：
式中，c1s,c1e为学习因子c1的初末值，分别取2.5、0.5；c2s,c2e为学习因子c2的初末值，分别取0.5、2.5；t为目前的迭代次数，tmax是最大迭代次数；
步骤3-6：计算当前个体目标函数值，根据更新当前个体至最优pbest，同时根据更新群体至最优gbest，其中t表示第t代循环，Xj t、Xj t+1分别表示第t、t+1次迭代第j个粒子，f(Xj t)、f(Xj t+1)表示第t、t+1次迭代第j个粒子适应性函数值；
步骤3-7：根据新的非劣解采用基于距离的方法维护外部档案EA：每次迭代都计算每个非劣解距离左右两个粒子的距离，通过比较该非劣解的两个距离是否小于某个限值Dis，如果均小于限值Dis，则从外部档案中删除此非劣解，并重新计算左右两边粒子的距离；
步骤3-8：基于Sigma选取策略来选择全局向导，为各粒子选择新的pbest，以双目标问题为例，f1和f2分别代表两个具有竞争关系的目标，即f1和f2不会同时达到最优值；
步骤3-9：判断是否达到全局迭代次数MAXGEN，若没有达到，GEN＝GEN+1，转到步骤3-5，继续下一轮的全局搜索，否则，算法结束，输出EA；
步骤四：基于外部归档集EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案。
2.如权利要求1所述的基于MAEPSO算法的区域水资源优化配置方法，其特征在于，所述步骤3-7根据新的非劣解采用基于距离的方法维护外部档案EA：每次迭代都计算每个非劣解距离左右两个粒子的距离，通过比较该非劣解的两个距离是否小于某个限值Dis，如果均小于限值Dis，则从外部档案中删除此非劣解，并重新计算左右两边粒子的距离具体为：
步骤3-7-1，将找到的所有非劣解集按照某个目标适应度值从小到大进行排序；
步骤3-7-2，计算非劣解与非劣解在目标空间的欧氏距离：
其中，n为目标的个数；
步骤3-7-3，判断上一步中计算的每个距离是否小于限值Dis，如果disi-1＜Dis且disi＜Dis，则删除非劣解i并重新计算非劣解i-1与非劣解i+1的欧氏距离作为disi-1；临界距离Dis设计为根据非劣解的情况自适应调整：
其中，NEA为外部档案EA中个体的数量；
步骤3-7-4，判断非劣解的个数是否超过外部档案的大小，若超过转3-7-2，否则转3-7-5；
步骤3-7-5，返回外部归档集EA。
3.如权利要求1所述的基于MAEPSO算法的区域水资源优化配置方法，其特征在于，所述步骤3-8，基于Sigma选取策略来选择全局向导，为各粒子选择新的pbest，以双目标问题为例，f1和f2分别代表两个目标，具体为：
步骤3-8-1，计算外部归档集EA中所有成员的Sigma值σa，如下式计算：
其中:σa就是Sigma值，f1,f2为对应的目标函数值，a＝1,2,…,NEA；
步骤3-8-2，从种群中选择一个粒子，计算该粒子的Sigma值σi；
步骤3-8-3，计算该粒子的Sigma值σi与外部档案所有成员的Sigma值之间的欧氏距离，选择与该粒子Sigma值欧氏距离最小的外部档案成员作为该粒子的pbest，更新粒子当前的速度和位置。
4.如权利要求1所述的基于MAEPSO算法的区域水资源优化配置方法，其特征在于，所述步骤四基于外部归档集EA，采用基于组合权重的多目标决策方法确定水资源系统最优配置方案，包括以下步骤：
步骤4-1，利用模糊层次分析法确定属性主观权重，包括以下步骤：
步骤4-1-1，根据所选取的各个目标函数，建立系统递阶层次结构，并进行各层次之间指标间的两两比较，采用0.1～0.9标度法进行定量描述，建立模糊一致判断矩阵R：
其中rij表示元素ai比aj重要的隶属度，rij越大，ai比aj越重要；n为指标数目；
步骤4-1-2，主观权重wj'计算
其中，参数a取值大小的选择反映了决策者的个人偏好，a越小则表明决策者偏好权重之间差异程度越大；rjk表示元素aj比ak重要的隶属度，j＝1,2,…,n，k＝1,2,…,n；
步骤4-2，利用CRITIC法确定属性客观权重，包括以下步骤：
步骤4-2-1，计算指标j对应的信息量Cj，指标i与其他指标之间的冲突性量化指标而用下式表示：
式中，δj为指标j的标准差，pjk为指标j,k的相关系数，j＝1,2,…,n，k＝1,2,…,n；
步骤4-2-2，所以指标j对应的归一化的权重wj”为：
步骤4-3，确定属性组合权重W：当得到各评价指标的主观权重wi'和客观权重wj”后，按下式计算各指标组合权重wj:
步骤4-4，基于TOPSIS法方案确定最终优化配置方案。
5.如权利要求4所述的基于MAEPSO算法的区域水资源优化配置方法，其特征在于，所述步骤4-4基于TOPSIS法方案确定最终优化配置方案具体如下：
步骤4-4-1，设多目标水资源优化配置问题的决策矩阵A为
式中，fij表示第i个可行方案第j个目标值，i＝1,2,…NEA,j＝1,2,…,n，NEA为外部档案EA中个体的数量；
由该矩阵A构成规范决策矩阵Z'，其元素为
步骤4-4-2，构造加权的规范决策矩阵Zij，其中的元素zij为zij＝wjrij；
步骤4-4-3，在加权规范化决策矩阵Zij中选取各个属性都达到最优的方案为正理想方案Z+，各个属性都达到最劣的方案为负理想方案Z-；
其中，n为效益型指标数目，n”为成本型指标数目，n＝n'+n”；
步骤4-4-4，采用欧几里得范数来计算各方案与正、负理想方案距离，不同方案到正、负理想方案的距离分别为：
式中，分别为第i个方案与正、负理想方案的欧几里得范数距离；zij为第i个方案第j个指标的规范化值；分别为第j个指标的正、负理想方案属性值；
步骤4-4-5，计算各个方案到正、负理想方案的相对接近度
步骤4-4-6，根据相对接近度指数的大小对样本方案进行优选排序；当方案越靠近正理想方案时，越接近1，取最大值对应的配置方案作为最优方案。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610973576.5A CN106529166B (zh) | 2016-11-04 | 2016-11-04 | 一种基于maepso算法的区域水资源优化配置方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610973576.5A CN106529166B (zh) | 2016-11-04 | 2016-11-04 | 一种基于maepso算法的区域水资源优化配置方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106529166A CN106529166A (zh) | 2017-03-22 |
| CN106529166B true CN106529166B (zh) | 2019-03-12 |
Family
ID=58349549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610973576.5A Expired - Fee Related CN106529166B (zh) | 2016-11-04 | 2016-11-04 | 一种基于maepso算法的区域水资源优化配置方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106529166B (zh) |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107229971B (zh) * | 2017-06-06 | 2019-10-08 | 西安电子科技大学 | 基于gapso算法的最优自适应策略决策方法 |
| CN107423564A (zh) * | 2017-07-26 | 2017-12-01 | 中国环境科学研究院 | 流域生态修复策略的决策分析方法 |
| CN107480813A (zh) * | 2017-07-27 | 2017-12-15 | 河海大学 | 基于多目标混沌遗传算法的流域水资源优化配置方法 |
| CN107748930B (zh) * | 2017-10-19 | 2020-06-02 | 扬州大学 | 一种充分灌溉条件下直接补渠的单库-多站系统水资源优化配置方法 |
| CN108049453A (zh) * | 2017-10-31 | 2018-05-18 | 江苏省水利勘测设计研究院有限公司 | 一种南水北调调水系统水资源配置装置 |
| CN108005016B (zh) * | 2017-10-31 | 2023-11-17 | 江苏省水利勘测设计研究院有限公司 | 一种南水北调调水系统水资源配置装置及配置方法 |
| CN108363303B (zh) * | 2018-03-09 | 2020-09-04 | 重庆科技学院 | 基于ar偏好信息的差分进化铝电解多目标优化方法 |
| CN110888732B (zh) * | 2018-09-10 | 2023-04-25 | 中国移动通信集团黑龙江有限公司 | 一种资源配置方法、设备、装置和计算机可读存储介质 |
| CN109446630A (zh) * | 2018-10-23 | 2019-03-08 | 宁夏大学 | 基于混合粒子群算法的多目标水资源优化配置方法 |
| CN109460624A (zh) * | 2018-11-30 | 2019-03-12 | 中国矿业大学 | 一种新的圆筒开关磁阻直线电机多目标优化设计方法 |
| CN110033164B (zh) * | 2019-03-04 | 2021-07-02 | 华中科技大学 | 一种水库群联合防洪调度的风险评估与决策方法 |
| CN110084428B (zh) * | 2019-04-26 | 2021-07-02 | 中国水利水电科学研究院 | 基于决策者偏好方案计算的水资源配置方法及系统 |
| CN110428099B (zh) * | 2019-07-30 | 2022-09-27 | 中国水利水电科学研究院 | 一种基于粒子群算法水库多目标供水能力优化方法 |
| CN111598447B (zh) * | 2020-05-15 | 2021-03-23 | 浙江大学 | 一种基于hmaqga的水库群联合优化调度方法 |
| CN112149246A (zh) * | 2020-09-11 | 2020-12-29 | 华中科技大学 | 一种永磁电机多目标优化设计方法和系统 |
| CN113240163B (zh) * | 2021-04-29 | 2024-02-27 | 中铁第一勘察设计院集团有限公司 | 复杂山区铁路线路多目标优化方法 |
| CN113312735B (zh) * | 2021-05-19 | 2022-06-03 | 太原理工大学 | 一种城市供水管网dma分区方法 |
| CN113539383B (zh) * | 2021-06-02 | 2022-10-04 | 中南大学 | 基于动态多目标优化和模糊多属性决策的除铜过程控制方法 |
| CN114204551A (zh) * | 2021-11-25 | 2022-03-18 | 哈尔滨工业大学 | 一种含光伏热泵的多能互补型微电网多元化优化运行方法、设备 |
| CN114358379B (zh) * | 2021-12-02 | 2024-01-23 | 武汉大学 | 梯级水库旱限水位调度线优化方法、装置及系统 |
| CN114492566A (zh) * | 2021-12-20 | 2022-05-13 | 西南科技大学 | 一种可调权重的高维数据降维方法及系统 |
| CN114997619B (zh) * | 2022-05-24 | 2023-06-20 | 中国标准化研究院 | 一种区域水资源动态调控平衡方法及系统 |
| CN115271437B (zh) * | 2022-07-27 | 2023-04-07 | 中国水利水电科学研究院 | 一种基于多决策主体的水资源配置方法及系统 |
| CN116777037B (zh) * | 2023-04-18 | 2024-01-30 | 中国农业大学 | 一种生态用水配水策略优化方法 |
| CN117649315A (zh) * | 2023-12-22 | 2024-03-05 | 中国水利水电科学研究院 | 基于水资源配置的空间均衡程度评价方法及装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103971174A (zh) * | 2014-05-06 | 2014-08-06 | 大连理工大学 | 基于改进量子粒子群算法的水电站群优化调度方法 |
| CN104156945A (zh) * | 2014-07-16 | 2014-11-19 | 西安电子科技大学 | 基于多目标粒子群算法的灰度图像分割方法 |
| CN105956718A (zh) * | 2016-05-31 | 2016-09-21 | 中国农业大学 | 一种灌区水资源优化配置方法和系统 |
-
2016
- 2016-11-04 CN CN201610973576.5A patent/CN106529166B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103971174A (zh) * | 2014-05-06 | 2014-08-06 | 大连理工大学 | 基于改进量子粒子群算法的水电站群优化调度方法 |
| CN104156945A (zh) * | 2014-07-16 | 2014-11-19 | 西安电子科技大学 | 基于多目标粒子群算法的灰度图像分割方法 |
| CN105956718A (zh) * | 2016-05-31 | 2016-09-21 | 中国农业大学 | 一种灌区水资源优化配置方法和系统 |
Non-Patent Citations (4)
| Title |
|---|
| MOPSO: a proposal for multiple objective particle swarm optimization;Coello C A C et al.;《IEEE Computer Society》;20021231;全文 |
| Strategies for finding good local guides in multi － objective particle swarm optimization;MOSTAGHIM S et al.;《Piscataway: IEEE Press》;20031231;全文 |
| 一种多尺度协同变异的粒子群优化算法;陶新民等;《软件学报》;20121231;全文 |
| 滨海地区水资源多目标优化调度模型研究;闫堃 等;《南 水 北 调 与 水 利 科 技》;20160229;第14卷(第1期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106529166A (zh) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106529166B (zh) | 一种基于maepso算法的区域水资源优化配置方法 | |
| WO2022193642A1 (zh) | 一种基于图卷积神经网络和nsga-ⅱ算法的水库调度多目标优化方法 | |
| CN105243458B (zh) | 一种基于多目标混合蛙跳差分算法的水库调度方法 | |
| CN109886473B (zh) | 一种考虑下游生态的流域风光水系统多目标优化调度方法 | |
| CN102629106B (zh) | 供水控制方法及系统 | |
| CN107527119A (zh) | 基于改进的多目标量子遗传算法的水资源优化调度方法 | |
| CN107370188A (zh) | 一种计及风电出力的电力系统多目标调度方法 | |
| CN109670650B (zh) | 基于多目标优化算法的梯级水库群调度模型的求解方法 | |
| Mohammadi et al. | Investigation of a new hybrid optimization algorithm performance in the optimal operation of multi-reservoir benchmark systems | |
| WO2022012093A1 (zh) | 一种基于"分区-分级"理论的复杂水资源系统优化配置方法 | |
| CN112184479B (zh) | 一种水库群调度规则型式和参数对气候变化的适应性研究方法 | |
| Yong | An improved artificial fish swarm algorithm for optimal operation of cascade reservoirs | |
| CN111598447B (zh) | 一种基于hmaqga的水库群联合优化调度方法 | |
| CN104933245A (zh) | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 | |
| CN106779177A (zh) | 基于粒子群优化的多分辨率小波神经网络用电量预测方法 | |
| CN110837915B (zh) | 一种基于混合集成深度学习的电力系统低压负荷点预测及概率预测方法 | |
| CN111724064B (zh) | 一种基于改进免疫算法的含储能配电网规划方法 | |
| CN103440535A (zh) | 基于免疫优化和模糊决策的多目标厂级负荷优化方法 | |
| CN111008790A (zh) | 一种水电站群发电调度规则提取方法 | |
| CN113191086A (zh) | 一种基于遗传算法的电采暖热负荷需求优化方法及系统 | |
| CN108376294A (zh) | 一种供能反馈与气象因素的热负荷预测方法 | |
| Akpinar et al. | Forecasting natural gas consumption with hybrid neural networks—Artificial bee colony | |
| CN106845696B (zh) | 一种智能优化水资源配置方法 | |
| CN115940267A (zh) | 一种区域配电网分布式光伏集群划分方法 | |
| CN112711896A (zh) | 一种考虑多源预报误差不确定性的复杂水库群优化调度方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190312 |