CN108491983A - 一种多堆燃料电池发电系统大规模机组组合优化方法 - Google Patents
一种多堆燃料电池发电系统大规模机组组合优化方法 Download PDFInfo
- Publication number
- CN108491983A CN108491983A CN201810318214.1A CN201810318214A CN108491983A CN 108491983 A CN108491983 A CN 108491983A CN 201810318214 A CN201810318214 A CN 201810318214A CN 108491983 A CN108491983 A CN 108491983A
- Authority
- CN
- China
- Prior art keywords
- fuel cell
- pile
- scene
- cell pile
- extreme value
- 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
- 239000000446 fuel Substances 0.000 title claims abstract description 164
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005457 optimization Methods 0.000 title claims abstract description 51
- 230000005611 electricity Effects 0.000 title claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 31
- 230000003044 adaptive effect Effects 0.000 claims abstract description 16
- 238000013459 approach Methods 0.000 claims abstract description 12
- 238000010248 power generation Methods 0.000 claims abstract description 12
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000012804 iterative process Methods 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000002485 combustion reaction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012913 prioritisation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000012528 membrane Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- H02J3/387—
-
- 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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Power Engineering (AREA)
- Fuel Cell (AREA)
Abstract
本发明公开一种多堆燃料电池发电系统大规模机组组合优化方法,包括步骤:通过K均值聚类法和蒙特卡罗模拟法,根据场景将燃料电池电堆聚类并建立概率分布模型;运用自适应惯性权值算法,使各燃料电池电堆根据概率分布模型在解空间中搜索追随当前的最优解,实现各燃料电池电堆出力功率寻优分配。本发明保证燃料电池发电阵列系统大规模机组的合理出力,通过自适应惯性权值算法控制每个电堆的启停和电堆之间的轮换,从而提高燃料电池的使用寿命,提高系统稳定性,降低系统成本,实现燃料经济性的目的。
Description
技术领域
本发明属于燃料电池技术领域,特别是涉及一种多堆燃料电池发电系统大规模机组组合优化方法。
背景技术
能源问题如今是不可忽视的严重问题,传统的化石能源日益减少甚至枯竭,已经无法满足人类的日常需求。于是越来越多新的优质能源涌现,同时需要掌握并运用先进的能源技术,新型、清洁和可再生能源的开发利用势在必行。燃料电池是一种清洁、高效的发电装置,其作为分布式电源的重要组成部分,以其诸多优点受到广泛关注。尤其是质子交换膜燃料电池(proton exchange membrane fuel cell,PEMFC)具有运行温度低、噪音低、效率高、响应速度快等特点,在燃料电池中应用最为广泛,是一种极具吸引力的发电装置。燃料电池发电系统在分布式发电、车辆、轨道交通和移动设备等领域均有广泛应用。
但是将燃料电池大规模应用于轨道交通领域仍受到诸多限制,尤其是燃料电池单体功率等级不高、可靠性低,难以满足轨道交通车辆实时功率需求,并且存在其耐久性不足、成本相对较高等问题。当燃料电池系统以阵列模式进行运行时,其系统效率能够得到较大提升。因此,多电堆轮换投入工作,综合考虑燃料电池各个电堆寿命,有必要针对燃料电池发电阵列系统大规模电堆组合的工作情况开展深入研究。燃料电池发电机组组合问题难点之一在于存在大量的整数变量,并且这些整数变量在时段上具有紧密的耦合性,因此对整数变量的处理是亟需重点解决的问题。
目前国内外对于燃料电池发电阵列的研究相对较少,针对多电堆大规模机组出力的控制没有相对较为详细的解释和说明,无法提高燃料电池寿命较短、降低成本。
发明内容
为了解决上述问题,本发明提出了一种多堆燃料电池发电系统大规模机组组合优化方法,保证燃料电池发电阵列系统大规模机组的合理出力,通过自适应惯性权值算法控制每个电堆的启停和电堆之间的轮换,从而提高燃料电池的使用寿命,提高系统稳定性,降低系统成本,实现燃料经济性的目的。
为达到上述目的,本发明采用的技术方案是:一种多堆燃料电池发电系统大规模机组组合优化方法,包括步骤:
步骤S100,通过K均值聚类法和蒙特卡罗模拟法,根据场景将燃料电池电堆聚类并建立概率分布模型;
步骤S200,运用自适应惯性权值算法,使各燃料电池电堆根据概率分布模型在解空间中搜索追随当前的最优解,实现各燃料电池电堆出力功率寻优分配。
进一步的是,步骤S100中,包括步骤:
将N个燃料电池电堆分成K个聚类,通过将聚类簇中具有较大的相似度的点,来完成各个场景的分类;所述场景为燃料电池发电系统中各个燃料电池电堆的使用情况;
生成情景的概率分布模型;根据燃料电池电堆情景建立非线性化模型,再通过求解非线性化模型得到情景树所有节点的情景及其概率,构成概率分布模型。
进一步的是,所述K均值聚类法包括步骤:首先,选取K个类中心;然后,将样本按规则归类并计算类的中心,不断变换类的中心重新归类;最后,重新得到K个聚类。
进一步的是,所述K均值聚类法中,每个场景到质心距离的概率距离计算公式为:
对于燃料电池发电系统的聚类目标函数为:
其中K为聚类数目即为燃料电池电堆群的数目,ck为k簇的中心,xj为k簇中的场景,优化目标为K个质心到簇内其它消减场景距离的总和最小;Xn为簇内 N个场景,ck为选出的K个质心场景。
进一步的是,所述蒙特卡罗模拟法包括步骤:首先,建立燃料电池发电阵列的概率分布模型;然后,确定模拟次数并随机抽样,经过数学计算得到随机解,随机解即为燃料电池电堆群的各个情景;最后,处理随机解数据,生成情景的概率分布模型。
进一步的是,所述生成情景的的概率分布模型为:
其中R表示所用统计特征的集合,包括燃料电池的使用情况以及整体系统的经济性;ωi表示统计特征i的权重;s表示情景树中各个节点的情景所构成的集合;Fi是事先已知的统计特征i的值;B是只由0和1构成的矩阵,它的每一行都表示一个节点的分布情况;p是按节点顺序将各情景的发生概率组合而成的列向量;e表示一个全为1的N维列向量。
进一步的是,在所述步骤S200中,综合考虑燃料电池的寿命和整体系统的成本因素,将燃料电池机组按空间划分为几个区域;在各个区域内的各次迭代中,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度;
所述极值包括个体极值和全局极值;所述个体极值为燃料电池电堆在迭代过程中找到的自身最优解,自身最优解所对应的燃料电池电堆是极值电堆pi,d;所述全局极值是各区域存在的局部最优解,局部最优解所对应的电堆是局部极值电堆pb,d;
基于所述极值的局部优化:开始时段随机分布燃料电池电堆的初始状态,将燃料电池电堆按使用程度划分为几个区域;在各次迭代中,每个燃料电池电堆通过跟踪2个极值来更新自身在解空间中的使用情况;
基于所述极值的全局优化:为了避免搜索过程中算法过早地陷入局部最优解,经过一定的迭代计算后,寻优过程改用全局优化算法以加速收敛;多个燃料电池电堆通过跟踪个体极值电堆pid和一个全局极值电堆pg,d更新自身的使用情况和内部状态,在解空间中搜寻全局最优解。
进一步的是,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数为:
其中,需要满足以下等式:
F(P)=aP+bP+c;
约束条件:It,j·Pmin,j<Pt,j<It,j·Pmax,j,
式中,NG和T分别为燃料电池电堆数目以及总时段数;F(Pt,j)和St,j分别为 t时段第j个燃料电池电堆的发电费用函数和启动耗量函数;It,j为0、1整数变量,0表示t时段第j个燃料电池电堆停机,1表示t时段第j个燃料电池电堆开机;Dt和Pt,L分别为t时段燃料电池系统负荷和整体系统损耗;Rt为t 时段燃料电池系统备用需求;Pmin,j和Pmax,j为第j个燃料电池电堆的最小和最大输出功率;zzt-1,j为燃料电池电堆j在时段t-1连续停运的时间;
每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数中对应每一个时段下的燃料电池机组启停计划即变量It,j的选择,都存在机组有功出力即连续变量Pt,j的优化选择;
再将燃料电池机组启停状态对应的0、1整数变量松弛为区间[0,1]内的连续变量,通过构造罚函数将目标函数对应的问题转换成单层的连续变量优化问题;再通过在算法迭代过程中对松弛后的0、1变量与机组有功出力变量并行地进行优化,避免了由于决策变量过多造成的维数灾难问题。引入变动阈值,解决了在寻优过程中各个电堆的值出现振荡时可能会丢失机组有效启停状态的问题。
通过引入罚函数方法和惯性权值自适应过程将约束下的机组组合问题转化为无约束问题求解。
进一步的是,将约束下的机组组合问题转化为无约束问题求解的目标函数为:
式中
其中,Δω为机组出力上升或下降速度限制;M,λ,η,v为罚因子。
进一步的是,在燃料电池阵列系统优化算法中引入自适应惯性权值算法,有效地处理了机组启、停时间的约束并提高了优化算法的全局收敛能力;所述自适应惯性权值算法,包括步骤:
首先,初始化各个燃料电池电堆的使用情况;
然后,进入循环判断程序,若达到迭代次数,计算燃料电池电堆的适值并求出个体极值与全局极值,并根据燃料电池电堆的个体极值与全局极值求出迭代后的燃料电池电堆使用情况;若没达到一定的迭代次数,计算燃料电池电堆的适值求出个体极值与局部极值,并根据个体极值与局部极值求出迭代后的燃料电池电堆使用情况;
随后,继续循环迭代进程,若达到一定的迭代次数,则进行整数运算;若没达到迭代次数,则进行惯性权值自适应调整过程,并更新所有燃料电池电堆的使用情况;
最后,判断迭代次数和寻优次数的完成情况,若完成寻优次数,则结束循环并输出数据;若没有完成寻优次数,返回循环进程。
采用本技术方案的有益效果:
本发明在一定运行周期内通过合理地开、停机组并满足各种运行条件的情况下达到燃料电池寿命损耗最低,并有效提高整体系统的经济性;该优化算法能够在达到最佳工作性能的同时延长多堆燃料电池的工作寿命,系统具有高度可扩展性和可靠性,能够在负载功率需求环境发生变化后进行自我调节,无需人为干预。
本发明以自适应惯性权值作为系统选择/动作等判定行为的依据,并能够根据运行工况的变化进行实时调整;同时在用于产生和修正处理惯性权值的函数中应有基于当前算法运行效能表现对权值进行修正的变量,当前算法运行效能表现的判定应由另一套预设的反映当前系统最优运行工况标准的函数实现,通过对所述惯性权值函数的不断运行和修正,可使系统始终运行于最优工况,并在所处外界环境发生变化之后通过自修正对最优工况进行持续跟踪,从而提高整体系统的自适应性和稳定性。
本发明采用的寻优方法能够避免搜索过程中算法过早地陷入局部最优解;经过一定的迭代计算后,寻优过程改用全局优化算法以加速收敛而且该优化算法采用实数编码并且没有许多参数需要调整,因此算法的实现更容易。机组组合问题就是在一定运行周期内通过合理地开、停机组并满足各种运行条件的情况下达到燃料电池的寿命损耗最小且燃料电池发电阵列系统总的运行成本最低。
本发明在算法迭代过程中并行优化松弛后的机组启停状态变量与机组有功出力变量,将两层优化问题转化为单层优化问题,降低了算法的时间复杂度;该算法有效地提高了大规模机组优化算法搜索全局最优解的能力和处理机组组合问题的收敛效率,能有效的节省计算时间、扩大计算规模;采用实数编码并且没有许多参数需要调整,因此该算法的实现更容易;该优化控制方法有较为广泛的应用范围,具有极高的通用性,能够应用于大规模机组组合问题的优化求解中;该优化方法针对多堆燃料电池系统投入工作有较大的改善,有效延长各个电堆的使用寿命,并提高整体系统的经济性。
附图说明
图1为本发明的一种多堆燃料电池发电系统大规模机组组合优化方法流程示意图;
图2为本发明实施例中K均值聚类法的流程示意图;
图3为本发明实施例中蒙特卡罗模拟法的流程示意图;
图4为本发明实施例中自适应惯性权值算法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步阐述。
在本实施例中,参见图1所示,本发明提出了一种多堆燃料电池发电系统大规模机组组合优化方法,包括步骤:
步骤S100,通过K均值聚类法和蒙特卡罗模拟法,根据场景将燃料电池电堆聚类并建立概率分布模型;
步骤S200,运用自适应惯性权值算法,使各燃料电池电堆根据概率分布模型在解空间中搜索追随当前的最优解,实现各燃料电池电堆出力功率寻优分配。
基于上述实施例的优化方案,步骤S100中,包括步骤:
将N个燃料电池电堆分成K个聚类,通过将聚类簇中具有较大的相似度的点,来完成各个场景的分类;所述场景为燃料电池发电系统中各个燃料电池电堆的使用情况;
生成情景的概率分布模型;根据燃料电池电堆情景建立非线性化模型,再通过求解非线性化模型得到情景树所有节点的情景及其概率,构成概率分布模型。
其中,如图2所示,所述K均值聚类法包括步骤:首先,选取K个类中心;然后,将样本按规则归类并计算类的中心,不断变换类的中心重新归类;最后,重新得到K个聚类。
所述K均值聚类法中,每个场景到质心距离的概率距离计算公式为:
对于燃料电池发电系统的聚类目标函数为:
其中K为聚类数目即为燃料电池电堆群的数目,ck为k簇的中心,xj为k簇中的场景,优化目标为K个质心到簇内其它消减场景距离的总和最小;Xn为簇内 N个场景,ck为选出的K个质心场景。
其中,如图3所示,所述蒙特卡罗模拟法包括步骤:首先,建立燃料电池发电阵列的概率分布模型;然后,确定模拟次数并随机抽样,经过数学计算得到随机解,随机解即为燃料电池电堆群的各个情景;最后,处理随机解数据,生成情景的概率分布模型。
所述生成情景的的概率分布模型为:
其中R表示所用统计特征的集合,包括燃料电池的使用情况以及整体系统的经济性;ωi表示统计特征i的权重;s表示情景树中各个节点的情景所构成的集合;Fi是事先已知的统计特征i的值;B是只由0和1构成的矩阵,它的每一行都表示一个节点的分布情况;p是按节点顺序将各情景的发生概率组合而成的列向量;e表示一个全为1的N维列向量。
基于上述实施例的优化方案,在所述步骤S200中,综合考虑燃料电池的寿命和整体系统的成本因素,将燃料电池机组按空间划分为几个区域;在各个区域内的各次迭代中,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度;
所述极值包括个体极值和全局极值;所述个体极值为燃料电池电堆在迭代过程中找到的自身最优解,自身最优解所对应的燃料电池电堆是极值电堆pi,d;所述全局极值是各区域存在的局部最优解,局部最优解所对应的电堆是局部极值电堆pb,d;
基于所述极值的局部优化:开始时段随机分布燃料电池电堆的初始状态,将燃料电池电堆按使用程度划分为几个区域;在各次迭代中,每个燃料电池电堆通过跟踪2个极值来更新自身在解空间中的使用情况;
基于所述极值的全局优化:为了避免搜索过程中算法过早地陷入局部最优解,经过一定的迭代计算后,寻优过程改用全局优化算法以加速收敛;多个燃料电池电堆通过跟踪个体极值电堆pid和一个全局极值电堆pg,d更新自身的使用情况和内部状态,在解空间中搜寻全局最优解。
其中,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数为:
其中,需要满足以下等式:
F(P)=aP+bP+c;
约束条件:It,j·Pmin,j<Pt,j<It,j·Pmax,j,
式中,NG和T分别为燃料电池电堆数目以及总时段数;F(Pt,j)和St,j分别为 t时段第j个燃料电池电堆的发电费用函数和启动耗量函数;It,j为0、1整数变量,0表示t时段第j个燃料电池电堆停机,1表示t时段第j个燃料电池电堆开机;Dt和Pt,L分别为t时段燃料电池系统负荷和整体系统损耗;Rt为t 时段燃料电池系统备用需求;Pmin,j和Pmax,j为第j个燃料电池电堆的最小和最大输出功率;zzt-1,j为燃料电池电堆j在时段t-1连续停运的时间;
每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数中对应每一个时段下的燃料电池机组启停计划即变量It,j的选择,都存在机组有功出力即连续变量Pt,j的优化选择;
再将燃料电池机组启停状态对应的0、1整数变量松弛为区间[0,1]内的连续变量,通过构造罚函数将目标函数对应的问题转换成单层的连续变量优化问题;再通过在算法迭代过程中对松弛后的0、1变量与机组有功出力变量并行地进行优化,避免了由于决策变量过多造成的维数灾难问题。引入变动阈值,解决了在寻优过程中各个电堆的值出现振荡时可能会丢失机组有效启停状态的问题。
通过引入罚函数方法和惯性权值自适应过程将约束下的机组组合问题转化为无约束问题求解。
其中,将约束下的机组组合问题转化为无约束问题求解的目标函数为:
式中
其中,Δω为机组出力上升或下降速度限制;M,λ,η,v为罚因子。
作为上述实施例的优化方案,在燃料电池阵列系统优化算法中引入自适应惯性权值算法,有效地处理了机组启、停时间的约束并提高了优化算法的全局收敛能力;如图4所示,所述自适应惯性权值算法,包括步骤:
首先,初始化各个燃料电池电堆的使用情况;
然后,进入循环判断程序,若达到迭代次数,计算燃料电池电堆的适值并求出个体极值与全局极值,并根据燃料电池电堆的个体极值与全局极值求出迭代后的燃料电池电堆使用情况;若没达到一定的迭代次数,计算燃料电池电堆的适值求出个体极值与局部极值,并根据个体极值与局部极值求出迭代后的燃料电池电堆使用情况;
随后,继续循环迭代进程,若达到一定的迭代次数,则进行整数运算;若没达到迭代次数,则进行惯性权值自适应调整过程,并更新所有燃料电池电堆的使用情况;
最后,判断迭代次数和寻优次数的完成情况,若完成寻优次数,则结束循环并输出数据;若没有完成寻优次数,返回循环进程。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (10)
1.一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,包括步骤:
步骤S100,通过K均值聚类法和蒙特卡罗模拟法,根据场景将燃料电池电堆聚类并建立概率分布模型;
步骤S200,运用自适应惯性权值算法,使各燃料电池电堆根据概率分布模型在解空间中搜索追随当前的最优解,实现各燃料电池电堆出力功率寻优分配。
2.根据权利要求1所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,步骤S100中,包括步骤:
将N个燃料电池电堆分成K个聚类,通过将聚类簇中具有较大的相似度的点,来完成各个场景的分类;所述场景为燃料电池发电系统中各个燃料电池电堆的使用情况;
生成情景的概率分布模型;根据燃料电池电堆情景建立非线性化模型,再通过求解非线性化模型得到情景树所有节点的情景及其概率,构成概率分布模型。
3.根据权利要求2所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,所述K均值聚类法包括步骤:首先,选取K个类中心;然后,将样本按规则归类并计算类的中心,不断变换类的中心重新归类;最后,重新得到K个聚类。
4.根据权利要求3所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,所述K均值聚类法中,每个场景到质心距离的概率距离计算公式为:
对于燃料电池发电系统的聚类目标函数为:
其中K为聚类数目即为燃料电池电堆群的数目,ck为k簇的中心,xj为k簇中的场景,优化目标为K个质心到簇内其它消减场景距离的总和最小;Xn为簇内N个场景,ck为选出的K个质心场景。
5.根据权利要求4所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,所述蒙特卡罗模拟法包括步骤:首先,建立燃料电池发电阵列的概率分布模型;然后,确定模拟次数并随机抽样,经过数学计算得到随机解,随机解即为燃料电池电堆群的各个情景;最后,处理随机解数据,生成情景的概率分布模型。
6.根据权利要求5所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,所述生成情景的的概率分布模型为:
其中R表示所用统计特征的集合,包括燃料电池的使用情况以及整体系统的经济性;ωi表示统计特征i的权重;s表示情景树中各个节点的情景所构成的集合;Fi是事先已知的统计特征i的值;B是只由0和1构成的矩阵,它的每一行都表示一个节点的分布情况;p是按节点顺序将各情景的发生概率组合而成的列向量;e表示一个全为1的N维列向量。
7.根据权利要求1所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,在所述步骤S200中,综合考虑燃料电池的寿命和整体系统的成本因素,将燃料电池机组按空间划分为几个区域;在各个区域内的各次迭代中,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度;
所述极值包括个体极值和全局极值;所述个体极值为燃料电池电堆在迭代过程中找到的自身最优解,自身最优解所对应的燃料电池电堆是极值电堆pi,d;所述全局极值是各区域存在的局部最优解,局部最优解所对应的电堆是局部极值电堆pb,d;
基于所述极值的局部优化:开始时段随机分布燃料电池电堆的初始状态,将燃料电池电堆按使用程度划分为几个区域;在各次迭代中,每个燃料电池电堆通过跟踪2个极值来更新自身在解空间中的使用情况;
基于所述极值的全局优化:为了避免搜索过程中算法过早地陷入局部最优解,经过一定的迭代计算后,寻优过程改用全局优化算法以加速收敛;多个燃料电池电堆通过跟踪个体极值电堆pid和一个全局极值电堆pg,d更新自身的使用情况和内部状态,在解空间中搜寻全局最优解。
8.根据权利要求7所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数为:
其中,需要满足以下等式:
F(P)=aP+bP+c;
约束条件:It,j·Pmin,j<Pt,j<It,j·Pmax,j,
式中,NG和T分别为燃料电池电堆数目以及总时段数;F(Pt,j)和St,j分别为t时段第j个燃料电池电堆的发电费用函数和启动耗量函数;It,j为0、1整数变量,0表示t时段第j个燃料电池电堆停机,1表示t时段第j个燃料电池电堆开机;Dt和Pt,L分别为t时段燃料电池系统负荷和整体系统损耗;Rt为t时段燃料电池系统备用需求;Pmin,j和Pmax,j为第j个燃料电池电堆的最小和最大输出功率;zzt-1,j为燃料电池电堆j在时段t-1连续停运的时间;
每个燃料电池电堆通过跟踪极值来更新自身在解空间中的位置与速度的目标函数中对应每一个时段下的燃料电池机组启停计划即变量It,j的选择,都存在机组有功出力即连续变量Pt,j的优化选择;
再将燃料电池机组启停状态对应的0、1整数变量松弛为区间[0,1]内的连续变量,通过构造罚函数将目标函数对应的问题转换成单层的连续变量优化问题;
通过引入罚函数方法和惯性权值自适应过程将约束下的机组组合问题转化为无约束问题求解。
9.根据权利要求8所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,将约束下的机组组合问题转化为无约束问题求解的目标函数为:
式中
其中,Δω为机组出力上升或下降速度限制;M,λ,η,v为罚因子。
10.根据权利要求1-9中任一所述的一种多堆燃料电池发电系统大规模机组组合优化方法,其特征在于,所述自适应惯性权值算法,包括步骤:
首先,初始化各个燃料电池电堆的使用情况;
然后,进入循环判断程序,若达到迭代次数,计算燃料电池电堆的适值并求出个体极值与全局极值,并根据燃料电池电堆的个体极值与全局极值求出迭代后的燃料电池电堆使用情况;若没达到一定的迭代次数,计算燃料电池电堆的适值求出个体极值与局部极值,并根据个体极值与局部极值求出迭代后的燃料电池电堆使用情况;
随后,继续循环迭代进程,若达到一定的迭代次数,则进行整数运算;若没达到迭代次数,则进行惯性权值自适应调整过程,并更新所有燃料电池电堆的使用情况;
最后,判断迭代次数和寻优次数的完成情况,若完成寻优次数,则结束循环并输出数据;若没有完成寻优次数,返回循环进程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810318214.1A CN108491983B (zh) | 2018-04-11 | 2018-04-11 | 一种多堆燃料电池发电系统大规模机组组合优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810318214.1A CN108491983B (zh) | 2018-04-11 | 2018-04-11 | 一种多堆燃料电池发电系统大规模机组组合优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108491983A true CN108491983A (zh) | 2018-09-04 |
CN108491983B CN108491983B (zh) | 2021-09-03 |
Family
ID=63315618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810318214.1A Active CN108491983B (zh) | 2018-04-11 | 2018-04-11 | 一种多堆燃料电池发电系统大规模机组组合优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108491983B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109685259A (zh) * | 2018-12-12 | 2019-04-26 | 聊城大学 | 基于多模态特性感知的混合流水调度进化优化方法 |
CN110190628A (zh) * | 2019-06-12 | 2019-08-30 | 西南交通大学 | 一种基于综合能源的复线全并联at牵引供电系统及方法 |
CN111755719A (zh) * | 2020-05-15 | 2020-10-09 | 广东鸿力氢动科技有限公司 | 一种燃料电池系统集群的功率分配方法及装置 |
CN111953035A (zh) * | 2020-07-17 | 2020-11-17 | 华中科技大学鄂州工业技术研究院 | 一种多模块固体氧化物燃料电池能量管理系统 |
CN116504069A (zh) * | 2023-06-26 | 2023-07-28 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110244353A1 (en) * | 2010-04-02 | 2011-10-06 | Chung-Hsin Electric And Machinery Manufacturing Corp. | Method of obtaining optimal design for a header of fuel cell stack and fuel cell stack with an optimally designed header |
CN105719015A (zh) * | 2016-01-19 | 2016-06-29 | 国网河北省电力公司电力科学研究院 | 一种基于pepso的分布式电源优化选址定容方法 |
CN106230028A (zh) * | 2016-09-08 | 2016-12-14 | 安徽电气工程职业技术学院 | 一种风电—抽水蓄能联合系统的多目标优化方法 |
CN106651059A (zh) * | 2017-01-13 | 2017-05-10 | 国网山西省电力公司 | 一种电动汽车充电站优化配置的方法 |
US20170301934A1 (en) * | 2012-10-16 | 2017-10-19 | Bloom Energy Corporation | Energy load management system |
CN107425521A (zh) * | 2017-06-28 | 2017-12-01 | 国网江苏省电力公司电力科学研究院 | 一种考虑负荷波动性的动态聚合建模及参与调频方法 |
CN107634518A (zh) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | 一种“源‑网‑荷”相协调的主动配电网经济调度方法 |
-
2018
- 2018-04-11 CN CN201810318214.1A patent/CN108491983B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110244353A1 (en) * | 2010-04-02 | 2011-10-06 | Chung-Hsin Electric And Machinery Manufacturing Corp. | Method of obtaining optimal design for a header of fuel cell stack and fuel cell stack with an optimally designed header |
US20170301934A1 (en) * | 2012-10-16 | 2017-10-19 | Bloom Energy Corporation | Energy load management system |
CN105719015A (zh) * | 2016-01-19 | 2016-06-29 | 国网河北省电力公司电力科学研究院 | 一种基于pepso的分布式电源优化选址定容方法 |
CN106230028A (zh) * | 2016-09-08 | 2016-12-14 | 安徽电气工程职业技术学院 | 一种风电—抽水蓄能联合系统的多目标优化方法 |
CN106651059A (zh) * | 2017-01-13 | 2017-05-10 | 国网山西省电力公司 | 一种电动汽车充电站优化配置的方法 |
CN107425521A (zh) * | 2017-06-28 | 2017-12-01 | 国网江苏省电力公司电力科学研究院 | 一种考虑负荷波动性的动态聚合建模及参与调频方法 |
CN107634518A (zh) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | 一种“源‑网‑荷”相协调的主动配电网经济调度方法 |
Non-Patent Citations (1)
Title |
---|
朱亚男 等: ""2kW 燃料电池分布式发电系统设计"", 《2017中国自动化大会(CAC2017)暨国际智能制造创新大会(CIMIC2017)论文集》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109685259A (zh) * | 2018-12-12 | 2019-04-26 | 聊城大学 | 基于多模态特性感知的混合流水调度进化优化方法 |
CN110190628A (zh) * | 2019-06-12 | 2019-08-30 | 西南交通大学 | 一种基于综合能源的复线全并联at牵引供电系统及方法 |
CN110190628B (zh) * | 2019-06-12 | 2022-04-15 | 西南交通大学 | 一种基于综合能源的复线全并联at牵引供电系统及方法 |
CN111755719A (zh) * | 2020-05-15 | 2020-10-09 | 广东鸿力氢动科技有限公司 | 一种燃料电池系统集群的功率分配方法及装置 |
CN111755719B (zh) * | 2020-05-15 | 2021-10-12 | 广东鸿力氢动科技有限公司 | 一种燃料电池系统集群的功率分配方法及装置 |
CN111953035A (zh) * | 2020-07-17 | 2020-11-17 | 华中科技大学鄂州工业技术研究院 | 一种多模块固体氧化物燃料电池能量管理系统 |
CN111953035B (zh) * | 2020-07-17 | 2021-11-30 | 华中科技大学鄂州工业技术研究院 | 一种多模块固体氧化物燃料电池能量管理系统 |
CN116504069A (zh) * | 2023-06-26 | 2023-07-28 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
CN116504069B (zh) * | 2023-06-26 | 2023-09-05 | 中国市政工程西南设计研究总院有限公司 | 城市道路网络容量优化方法、装置、设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108491983B (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108491983A (zh) | 一种多堆燃料电池发电系统大规模机组组合优化方法 | |
Wang et al. | Optimal capacity allocation of standalone wind/solar/battery hybrid power system based on improved particle swarm optimisation algorithm | |
Welch et al. | Energy dispatch fuzzy controller for a grid-independent photovoltaic system | |
Yang et al. | A novel parallel-series hybrid meta-heuristic method for solving a hybrid unit commitment problem | |
CN108471143A (zh) | 基于正负反馈粒子群算法的微电网多能源调度优化方法 | |
CN104036329A (zh) | 一种基于多智能体协同寻优的含光伏微源主动配网拓扑重构方法 | |
CN109034587B (zh) | 一种协调多种可控单元的主动配电系统优化调度方法 | |
CN106684913B (zh) | 一种基于多代理的储能电站跟踪发电计划控制系统和方法 | |
Hu et al. | Short-term photovoltaic power prediction based on similar days and improved SOA-DBN model | |
CN107273968A (zh) | 一种基于动态模糊混沌粒子群的多目标调度方法及装置 | |
CN112163711A (zh) | 一种基于mopso/cc的综合能源系统优化调度方法 | |
Yang et al. | A Novel Binary Jaya Optimization for Economic/Emission Unit Commitment. | |
Deng et al. | Hybrid short-term load forecasting using the Hadoop MapReduce framework | |
CN109950933A (zh) | 一种基于改进粒子群算法的风光储联合调峰优化方法 | |
CN111461478B (zh) | 一种大规模水光能源互补调度方法与系统 | |
Li et al. | Distributed power dispatching solution for a future economic and environment-friendly energy internet | |
Sedighizadeh et al. | Parameter identification of proton exchange membrane fuel cell using a Hybrid Big Bang-Big Crunch optimization | |
CN108022055A (zh) | 一种基于粒子群模型的微电网经济调度方法 | |
Wu et al. | Brain Storm Optimization Algorithm based on adaptive inertial selection strategy for the RCPSP | |
Li et al. | A review of multi-objective particle swarm optimization algorithms in power system economic dispatch | |
Zhao et al. | A novel binary social learning particle swarm optimizer for power system unit commitment | |
Hui-chcng et al. | Solving on agile satellites mission planning based on tabu search-parallel genetic algorithms | |
Dong et al. | Dynamic reconfiguration strategy based on partition of time intervals with improved fuzzy C-means clustering | |
Cao et al. | Predictive model based on graph convolutional neural network | |
Ghiassi et al. | On the use of AI as a requirement for improved insolation forecasting accuracy to achieve optimized PV utilization |
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 |