CN103345663B - 考虑爬坡速率约束的电力系统机组组合优化方法 - Google Patents

考虑爬坡速率约束的电力系统机组组合优化方法 Download PDF

Info

Publication number
CN103345663B
CN103345663B CN201310301711.8A CN201310301711A CN103345663B CN 103345663 B CN103345663 B CN 103345663B CN 201310301711 A CN201310301711 A CN 201310301711A CN 103345663 B CN103345663 B CN 103345663B
Authority
CN
China
Prior art keywords
unit
individual
population
individuality
sit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310301711.8A
Other languages
English (en)
Other versions
CN103345663A (zh
Inventor
张景瑞
林爽
邱卫霞
郭逸然
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xiamen University
Original Assignee
Xiamen University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xiamen University filed Critical Xiamen University
Priority to CN201310301711.8A priority Critical patent/CN103345663B/zh
Publication of CN103345663A publication Critical patent/CN103345663A/zh
Application granted granted Critical
Publication of CN103345663B publication Critical patent/CN103345663B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

考虑爬坡速率约束的电力系统机组组合优化方法,涉及电力系统调度。以煤耗成本和启动成本之和最小为目标,在离散粒子群框架中引入差分加速技术以提高求解速度,对不可行个体进行修复以提高可行性,并采用随机等效λ迭代法进行负荷分配并处理爬坡速率约束。具体包括以下步骤:收集相关数据并设置算法参数;建立机组组合优化数学模型;随机生成初始种群;对不满足约束的个体进行修复;采用随机等效λ迭代法进行经济负荷分配并处理爬坡速率约束;计算适应度值,更新种群最好和个体最好;执行差分加速搜索;更新种群;达到迭代结束条件则输出最优机组组合方案。可以显著提高求解速度,有效降低运行成本。

Description

考虑爬坡速率约束的电力系统机组组合优化方法
技术领域
本发明涉及电力系统调度,尤其是涉及一种考虑爬坡速率约束的电力系统机组组合优化方法。
背景技术
随着电力系统规模的增加和日益加剧的能源危机,人们对电力系统经济性的运行要求越来越高。电力系统机组组合优化对电力系统经济运行、安全调度具有重要影响,它不仅可以节省大量的经济成本,而且能通过保持一定的旋转备用来提高电力系统的可靠性。机组组合优化问题涉及两个子问题,其一是机组组合,用于确定有哪些机组出力;其二为经济负荷分配,用于确定需要这些机组出多少力。机组组合优化问题的决策变量不仅涉及表示机组运行状态(离线、在线分别用0、1表示)的离散变量,而且涉及表示机组出力的连续变量,需要考虑包括电量平衡、旋转备用等在内的大量线性、非线性的等式或不等式约束。
目前,机组组合优化的方法有优先顺序法、动态规划法、拉格朗日松弛法等传统数学规划方法和遗传算法、神经网络等智能进化算法。优先顺序法计算简单,速度快,但很难得到最优解;动态规划法能得到最优解,但计算量大,存在维数灾问题;拉格朗日松弛法算子选择不易,且其优劣直接影响到算法的收敛性;遗传算法容易早熟收敛,且后期易出现震荡现象;神经网络的训练需要大量的样本数据。
Kennedy和Eberhart于1995年提出的粒子群算法是模拟鸟群觅食的基于群体智能的仿生算法,便于在连续变量空间处理全局约束优化问题,随后他们于1997年将所提出的粒子群算法扩展为离散粒子群算法,用于处理含离散变量的优化问题,具有原理简单、鲁棒性好和易于实现的优点。Zhao等人(B.Zhao,C.X.Guo,B.R.Bai and Y.J.Cao.An improvedparticle swarm optimization algorithm for unit commitment.InternationalJournal of Electrical Power&Energy Systems,28(7):482-490,2006)在2006年构建了一个增强粒子群算法求解电力系统机组组合优化问题,但该法的求解速度慢,容易陷入局部最优,而且该法所求得的机组组合优化的结果并不理想。中国专利201110243050.9(浙江大学)公开一种基于离散粒子群和拉格朗日乘数算法的机组组合优化方法,但该法不一定能获得满足各种运行约束的最优解(其实施案例所给出的验证性结果未满足最小在线时间约束)且未考虑爬坡速率约束,求解速度慢,不适于实际应用。差分加速技术是Zhang等人(Jingrui Zhang,Jian Wang,and Chaoyuan Yue,Small population-based particleswarm optimization for short-term hydrothermal scheduling,IEEE Transactionson Power Systems,27(1):142-152,2012)于2012年提出的短期水火联调问题的粒子群算法求解框架中用于加速搜索的一种方法,但该法仅针对连续变量,对含离散变量的机组组合优化无能为力。
发明内容
本发明的目的在于为进一步提高粒子群算法求解机组组合优化问题的速度,提供一种考虑爬坡速率约束的电力系统机组组合优化方法。
本发明针对含爬坡速率约束的电力系统机组组合优化问题,将差分加速技术扩展于处理离散变量的加速搜索,在经济负荷分配上引入适于智能进化算法的随机等效Lambda(在本说明书剩余部分采用λ表示)迭代法。
本发明包括以下步骤:
(1)收集电力系统相关数据并设置算法参数;
(2)根据步骤(1)所收集的电力系统相关数据,建立考虑各种运行约束的电力系统机组组合优化数学模型;
(3)根据步骤(2)所建立的电力系统机组组合优化问题数学模型和步骤(1)所设置的算法参数,初始化种群个体的飞行速度,并确定种群个体的位置也即是机组的运行状态;
(4)判断种群个体是否满足最小在线/离线时间约束,若满足,则转步骤(5);若不满足,则对种群个体进行最小在线/离线时间约束的修复;
(5)判断满足最小在线/离线时间约束的种群个体是否满足旋转备用约束,若满足,则转步骤(6);若不满足,则对种群个体进行旋转备用约束的修复;
(6)对种群个体执行过度盈余机组删除操作;
(7)根据步骤(1)所收集的不同时段的预测负荷情况对种群个体进行经济负荷分配;
(8)根据步骤(2)所建立的优化问题目标函数和步骤(7)的经济负荷分配情况,计算种群个体的目标函数值,设置该值为相应个体的适应度值,并更新个体最好和种群最好;
(9)采用差分加速技术对种群个体最好进行加速搜索;
(10)若达到迭代结束条件,则转步骤(12),否则转下一步;
(11)根据个体最好和种群最好更新种群个体的飞行速度和位置,转步骤(4);
(12)输出电力系统机组组合优化方案,求解结束。
在步骤(1)中,所述电力系统相关数据包括参与调度机组的特性数据、调度周期及其时段数、调度周期内各时段的预测负荷和旋转备用要求;所述特性数据包括煤耗函数、启动成本、最小在线/离线时间、爬坡速率约束、机组出力上下限及机组初始运行状态;所述算法参数包括粒子群算法相关参数、差分加速方法相关参数和随机等效λ迭代法的相关参数。
在步骤(2)中,所述电力系统机组组合优化数学模型以机组煤耗成本和启动成本之和最小为目标,优化的目标函数优化的约束条件为:
1)电量平衡
t=1,2,...,T
2)旋转备用约束
t=1,2,...,T
3)最小在线时间
uit=1∧uit+1=0
4)最小离线时间
uit=0∧uit+1=1
5)机组出力约束
i=1,2,...,Ns,t=1,2,...,T
6)火电机组初始运行状态
i=1,2,...,Ns
7)机组爬坡速率约束
Psit-Psit-1≤URi,Psit-1-Psit≤UDii=1,2,...,Ns,t=1,2,...,T
在上述公式中,F是优化的目标,t=1,2,...,T代表调度时段,i=1,2,...,Ns代表参与调度的机组,Psit代表机组i在时段t的出力,机组i的煤耗成本可以表示为相应出力的二次多项式函数:fi(Psit)=ai+biPsit+ciPsit 2,aibi和ci为机组i的煤耗成本函数的系数,uit为机组i在t时段工作状态,1表示在线,0表示离线,机组i在t时段的启动费用SUit可以表示为:
Chsi为机组i的热启动成本,Ccsi为机组i的冷启动成本,Hcsi为机组i的冷启动时间,为机组i的最小离线时间,机组i在t时段末的连续离线时间更新,为机组i的最小在线时间,机组i在t时段末的连续在线时间PL(t)和PD(t)分别为t时段的输电损耗和系统负荷,R(t)为t时段的旋转备用要求,分别为机组i出力的上下限,分别为机组i的初始连续在线和离线时间,URi和DRi分别为机组i出力增加和降低的速度上限也即是反映爬坡速率约束的两个变量。在不加说明情况下,本发明说明书剩余部分相同数学符号均代表上述意义。
在步骤(4)中,所述最小在线/离线时间约束修复的步骤如下:
1):获取待修复个体,初始化令t=1,i=1;
2):若uit=1,则转步骤3),否则转步骤5);
3):若uit-1=0,则转步骤4),否则转步骤7);
4):若则令uit=0转步骤7),否则直接转步骤7);
5):若uit-1=1,则转步骤6),否则转步骤7);
6):若则令uit=1转步骤7),否则直接转步骤7);
7):更新
8):i=i+1,若i≤Ns,则转步骤2),否则转步骤9);
9):t=t+1,若t≤T,则i=1转步骤2),否则转步骤10);
10):最小在线/离线时间约束修复完毕。
在步骤(5)中,所述旋转备用约束的修复的步骤如下:
1):获取待修复个体,将机组按照装机容量排成降序,令t=1;
2):按下式计算当前个体在t时段的旋转备用情况令g=1,转步骤3),否则转步骤7);
3):若ugt=1,则g=g+1直至ugt=0,计算
4):令ugt=1,若 否则令
5):令ugl=1, l=l+1;
6):若l≤t,则转步骤5),否则令则转步骤7),否则g=g+1转步骤3);
7):t=t+1,若t≤T,则转步骤2),否则转步骤8);
8):旋转备用约束修复完毕;
在上述步骤中,g代表排序后的机组序号,l代表时段号,代表t时段的实际旋转备用。
在步骤(6)中,所述过度盈余机组删除操作的具体步骤如下:
1):获取待修复个体,按容量大小将机组排成升序,并令g=1,t=1;
2):若uit=1,则转步骤3),否则转步骤7);
3):若则转步骤4),否则转步骤8);
4):若则转步骤6),否则转步骤5);
5):若则转步骤6),否则转步骤7);
6):ugt=0,更新
7):g=g+1,若g≤Ns,则转步骤2),否则转步骤8);
8):t=t+1,若t≤T,则令g=1转步骤2),否则转步骤9);
9):过度盈余机组删除完毕;
在上述步骤中,g代表排序后的机组序号。
在步骤(7)中,所述对种群个体进行经济负荷分配可采用随机等效λ迭代法,并在分配中处理爬坡速率约束,可按时段t=1,2,...,T进行,其具体步骤如下:
1):获取待分配个体和时段t,令iteration=1,i=1,λ=rand(λminmax);
2):求的解为
3):令Psit=min{max{Pi',Psit-1-UDi},Psit-1+URi};
4):若i<Ns,则i=i+1,转步骤2),否则计算
5):若ε≥τ,则转步骤6),否则转步骤10);
6):若iteration=1,则λ1=λ,ε1=ε转步骤7),否则λ2=λ1,ε2=ε1,λ1=λ,ε1=ε,转步骤9);
7):若ε>0,则λ=(1+η)*λ,否则λ=(1-η)*λ;
8):iteration=iteration+1,i=1转步骤2);
9):若|ε12|<ξ,则转步骤7),否则令λ=λ2+(λ12)*ε2/(ε21)转步骤8);
10):待分配个体t时段的经济负荷分配结束;
在上述步骤中,iteration代表迭代次数,和Pi'是为获得Psit的中间值,λ1和λ2是λ在前两次迭代中的值,ε1和ε2是前两次迭代中电量平衡的违背量,λmin和λmax分别是λ的最小、最大值,rand(λminmax)产生在[λminmax]上服从均匀分布的随机数,τ为电量平衡允许误差,η为λ的更新步长,ξ为近两次迭代误差的阈值。
在步骤(9)中,所述采用差分加速技术对种群个体最好进行加速搜索,可直接以种群个体最好为差分加速的初始种群,以pbest表示个体最好,gbest表示种群最好,表示在第m次迭代中个体r1和个体r2的个体最好在第j维的距离,则差分加速搜索的详细步骤如下:
1):以种群当前个体最好为差分加速的初始种群,令k=1,i=1;
2):按式产生个体i的尝试个体,式中,tmp为尝试个体,mut反映其变异部分,且 为位反转函数,
3):计算尝试个体tmp的适应度值;
4):若tmp的适应度值优于个体i的个体最好适应度值,则将尝试个体更新为个体i的个体最好,并转下一步,否则转步骤6);
5):将个体i的个体最好适应度值与种群最好适应度值进行比较,若优于种群最好适应度值,则更新种群最好;
6):i=i+1,判断是否遍历种群,若已遍历则转步骤7),否则转步骤2);
7):k=k+1,判断是否达到加速次数,若达到则转步骤8),否则,令i=1,转步骤2);
8):差分加速搜索结束;
在上述步骤中,i代表个体序号,j代表维数,k代表差分加速迭代次数,m代表粒子群算法迭代次数,Facc和CRacc分别为差分加速搜索方法的收缩因子和交叉因子。
本发明要解决的技术问题是,克服现有技术的不足,在求解考虑爬坡速率约束的电力系统机组组合优化问题的粒子群算法框架中引入差分加速技术加速搜索,提高机组组合优化问题求解的速度;在机组组合中引入对不可行种群个体进行修复的一系列方法,提高机组组合方案的可行性;在经济负荷分配中采用适于智能进化算法的随机等效λ迭代法并处理爬坡速率约束。总之,本发明为考虑爬坡速率约束的电力系统机组组合优化问题提供了一种基于差分加速和不可行个体修复的粒子群优化方法,该方法根据所预测的电力系统负荷数据和参与调度机组的特性参数,以煤耗成本和启动成本之和最小为目标,给出在一个调度周期内满足各种约束的运行方案。
本发明旨在提供一种考虑爬坡速率约束的电力系统机组组合优化方法。本发明以煤耗成本和启动成本之和最小为目标,在离散粒子群框架中引入差分加速技术以提高求解速度,对不可行个体进行修复以提高可行性,并采用随机等效λ迭代法进行负荷分配并处理爬坡速率约束。具体包括以下步骤:收集相关数据并设置算法参数;建立机组组合优化数学模型;随机生成初始种群;对不满足约束的个体进行修复;采用随机等效λ迭代法进行经济负荷分配并处理爬坡速率约束;计算适应度值,更新种群最好和个体最好;执行差分加速搜索;更新种群;达到迭代结束条件则输出最优机组组合方案。采用本发明可以显著提高求解速度,有效降低运行成本。
附图说明
图1为本发明实施例流程图。
图2为本发明方法与其他方法的收敛曲线比较图。在图2中,a为本发明方法,b为离散粒子群,c为增强粒子群。
具体实施方式
下面结合附图,以一个含10机组的典型电力系统为例,对本发明所述的考虑爬坡速率约束的电力系统机组组合优化方法的实施作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。如图1所示是本发明的流程图,所提供方法包括如下步骤:
(1)收集电力系统相关数据并设置算法参数。
电力系统相关数据包括:参与调度机组的特性数据如表1所示(该表中机组初始状态为负值表示初始连续离线时间,为正表示初始连续在线时间);调度周期为1天并分为24时段;调度周期内各时段的预测负荷如表2所示,旋转备用要求为相应时段负荷的10%;设置微粒群算法、差分加速搜索方法及随机等效λ迭代法的相关参数数据。
表1机组特性参数
表2时段负荷(MW)
(2)根据步骤(1)所收集的电力系统相关数据,建立考虑各种运行约束的电力系统机组组合优化数学模型。
(3)根据步骤(2)所建立的电力系统机组组合优化数学模型和步骤(1)所设置的算法参数,初始化种群个体的飞行速度,并确定种群个体的位置也即是机组的运行状态。
(4)判断种群个体是否满足最小在线/离线时间约束,若满足则转步骤(6);
(5)对不满足最小在线/离线时间约束的种群个体进行最小在线/离线时间约束的修复;
(6)判断满足最小在线/离线时间约束的种群个体是否满足旋转备用约束,若满足,则转步骤(8);
(7)对不满足旋转备用约束的种群个体进行旋转备用约束的修复;
(8)对种群个体执行过度盈余机组删除操作;
(9)根据步骤(1)所收集的不同时段的预测负荷情况对种群个体进行经济负荷分配;
(10)根据步骤(2)所建立的优化问题目标函数和步骤(9)的经济负荷分配情况,计算种群个体的目标函数值,设置该值为相应个体的适应度值,并更新个体最好和种群最好;
(11)采用差分加速方法对种群个体最好进行加速搜索;
(12)若达到迭代结束条件,则转步骤(14),否则转下一步;
(13)根据个体最好和种群最好更新种群个体的飞行速度和位置,转步骤(4);
(14)输出电力系统机组组合优化方案,求解结束。
经上述方法优化后的10机组24时段的负荷及相应的成本如表3所示。
表3优化后的各机组不同时段出力安排及成本
在表3中,机组时段出力为0表示该机组在该时段处于离线状态即停机状态,从表3可见,所得分配方案满足所有约束条件,总运行成本为$563937。本发明所述方法与其他方法(包括离散粒子群算法和增强粒子群算法)的收敛性比较曲线见图2。经对不同方法的优化结果进行对比分析和比较,可知本发明所述方法能够快速收敛。可见本发明所述方法在电网调度计划编制中能有效地分配资源,提高电网运行的经济性,速度快,效率高,具有广泛推广的应用前景。

Claims (1)

1.考虑爬坡速率约束的电力系统机组组合优化方法,其特征在于包括以下步骤:
(1)收集电力系统相关数据并设置算法参数;
所述电力系统相关数据包括参与调度机组的特性数据、调度周期及其时段数、调度周期内各时段的预测负荷和旋转备用要求;所述特性数据包括煤耗函数、启动成本、最小在线/离线时间、爬坡速率约束、机组出力上下限及机组初始运行状态;所述算法参数包括粒子群算法相关参数、差分加速方法相关参数和随机等效λ迭代法的相关参数;
(2)根据步骤(1)所收集的电力系统相关数据,建立考虑各种运行约束的电力系统机组组合优化数学模型;
所述电力系统机组组合优化数学模型以机组煤耗成本和启动成本之和最小为目标,优化的目标函数优化的约束条件为:
1)电量平衡
Σ i = 1 N s u i t P s i t - P L ( t ) = P D ( t ) t = 1 , 2 , ... , T
2)旋转备用约束
Σ i = 1 N s u i t P s i max ≥ P D ( t ) + R ( t ) t = 1 , 2 , ... , T
3)最小在线时间
4)最小离线时间
5)机组出力约束
u i t P s i min ≤ P s i t ≤ u i t P s i max i = 1 , 2 , ... , N s , t = 1 , 2 , ... , T
6)火电机组初始运行状态
T i 0 o n = T i i n i t o n , T i 0 o f f = T i i n i t o f f i = 1 , 2 , ... , N s
7)机组爬坡速率约束
Psit-Psit-1≤URi,Psit-1-Psit≤UDi i=1,2,...,Ns,t=1,2,...,T
在上述公式中,F是优化的目标,t=1,2,...,T代表调度时段,i=1,2,...,Ns代表参与调度的机组,Psit代表机组i在时段t的出力,机组i的煤耗成本表示为相应出力的二次多项式函数:fi(Psit)=ai+biPsit+ciPsit 2,ai,bi和ci为机组i的煤耗成本函数的系数,uit为机组i在t时段工作状态,1表示在线,0表示离线,机组i在t时段的启动费用SUit表示为:
Chsi为机组i的热启动成本,Ccsi为机组i的冷启动成本,Hcsi为机组i的冷启动时间,Ti MD为机组i的最小离线时间,机组i在t时段末的连续离线时间更新,Ti MU为机组i的最小在线时间,机组i在t时段末的连续在线时间PL(t)和PD(t)分别为t时段的输电损耗和系统负荷,R(t)为t时段的旋转备用要求,分别为机组i出力的上下限,Ti initon和Ti initoff分别为机组i的初始连续在线和离线时间,URi和UDi分别为机组i出力增加和降低的速度上限也即是反映爬坡速率约束的两个变量;
(3)根据步骤(2)所建立的电力系统机组组合优化问题数学模型和步骤(1)所设置的算法参数,初始化种群个体的飞行速度,并确定种群个体的位置也即是机组的运行状态;
(4)判断种群个体是否满足最小在线/离线时间约束,若满足,则转步骤(5);若不满足,则对种群个体进行最小在线/离线时间约束的修复;
所述最小在线/离线时间约束修复的步骤如下:
1):获取待修复个体,初始化令t=1,i=1;
2):若uit=1,则转步骤3),否则转步骤5);
3):若uit-1=0,则转步骤4),否则转步骤7);
4):若则令uit=0转步骤7),否则直接转步骤7);
5):若uit-1=1,则转步骤6),否则转步骤7);
6):若则令uit=1转步骤7),否则直接转步骤7);
7):更新
8):i=i+1,若i≤Ns,则转步骤2),否则转步骤9);
9):t=t+1,若t≤T,则i=1转步骤2),否则转步骤10);
10):最小在线/离线时间约束修复完毕;
(5)判断满足最小在线/离线时间约束的种群个体是否满足旋转备用约束,若满足,则转步骤(6);若不满足,则对种群个体进行旋转备用约束的修复;
所述旋转备用约束的修复的步骤如下:
1):获取待修复个体,将机组按照装机容量排成降序,令t=1;
2):按下式计算当前个体在t时段的旋转备用情况令g=1,转步骤3),否则转步骤7);
r ‾ ( t ) = Σ i = 1 N s u i t P s i max P D ( t ) - R ( t )
3):若ugt=1,则g=g+1直至ugt=0,计算
4):令ugt=1,若否则令
5):令ugl=1,l=l+1;
6):若l≤t,则转步骤5),否则令则转步骤7),否则g=g+1转步骤3);
7):t=t+1,若t≤T,则转步骤2),否则转步骤8);
8):旋转备用约束修复完毕;
在上述步骤中,g代表排序后的机组序号,l代表时段号,代表t时段的实际旋转备用;
(6)对种群个体执行过度盈余机组删除操作;
所述过度盈余机组删除操作的具体步骤如下:
1):获取待修复个体,按容量大小将机组排成升序,并令g=1,t=1;
2):若uit=1,则转步骤3),否则转步骤7);
3):若则转步骤4),否则转步骤8);
4):若则转步骤6),否则转步骤5);
5):若则转步骤6),否则转步骤7);
6):ugt=0,更新
7):g=g+1,若g≤Ns,则转步骤2),否则转步骤8);
8):t=t+1,若t≤T,则令g=1转步骤2),否则转步骤9);
9):过度盈余机组删除完毕;
在上述步骤中,g代表排序后的机组序号;
(7)根据步骤(1)所收集的不同时段的预测负荷情况对种群个体进行经济负荷分配;
所述对种群个体进行经济负荷分配是采用随机等效λ迭代法,并在分配中处理爬坡速率约束,按时段t=1,2,...,T进行,其具体步骤如下:
1):获取待分配个体和时段t,令iteration=1,i=1,λ=rand(λminmax);
2):求的解为Pi *,令
3):令Psit=min{max{Pi',Psit-1-UDi},Psit-1+URi};
4):若i<Ns,则i=i+1,转步骤2),否则计算
5):若ε≥τ,则转步骤6),否则转步骤10);
6):若iteration=1,则λ1=λ,ε1=ε转步骤7),否则λ2=λ1,ε2=ε1,λ1=λ,ε1=ε,转步骤9);
7):若ε>0,则λ=(1+η)*λ,否则λ=(1-η)*λ;
8):iteration=iteration+1,i=1转步骤2);
9):若|ε1-ε2|<ξ,则转步骤7),否则令λ=λ2+(λ12)*ε2/(ε21)转步骤8);
10):待分配个体t时段的经济负荷分配结束;
在上述步骤中,iteration代表迭代次数,Pi *和Pi'是为获得Psit的中间值,λ1和λ2是λ在前两次迭代中的值,ε1和ε2是前两次迭代中电量平衡的违背量,λmin和λmax分别是λ的最小、最大值,rand(λminmax)产生在[λminmax]上服从均匀分布的随机数,τ为电量平衡允许误差,η为λ的更新步长,ξ为近两次迭代误差的阈值;
(8)根据步骤(2)所建立的优化问题目标函数和步骤(7)的经济负荷分配情况,计算种群个体的目标函数值,设置该值为相应个体的适应度值,并更新个体最好和种群最好;
(9)采用差分加速技术对种群个体最好进行加速搜索;所述差分加速技术处理离散变量的加速搜索;
所述采用差分加速技术对种群个体最好进行加速搜索,是直接以种群个体最好为差分加速的初始种群,以pbest表示个体最好,gbest表示种群最好,表示在第m次迭代中个体r1和个体r2的个体最好在第j维的距离,则差分加速搜索的详细步骤如下:
1):以种群当前个体最好为差分加速的初始种群,令k=1,i=1;
2):按式产生个体i的尝试个体,式中,tmp为尝试个体,mut反映其变异部分,且 为位反转函数,
f l i p ( gbest y k ) = 1 i f gbest j k = 0 0 i f gbest j k = 1 ;
3):计算尝试个体tmp的适应度值;
4):若tmp的适应度值优于个体i的个体最好适应度值,则将尝试个体更新为个体i的个体最好,并转下一步,否则转步骤6);
5):将个体i的个体最好适应度值与种群最好适应度值进行比较,若优于种群最好适应度值,则更新种群最好;
6):i=i+1,判断是否遍历种群,若已遍历则转步骤7),否则转步骤2);
7):k=k+1,判断是否达到加速次数,若达到则转步骤8),否则,令i=1,转步骤2);
8):差分加速搜索结束;
在上述步骤中,i代表个体序号,j代表维数,k代表差分加速迭代次数,m代表粒子群算法迭代次数,Facc和CRacc分别为差分加速搜索方法的收缩因子和交叉因子;
(10)若达到迭代结束条件,则转步骤(12),否则转下一步;
(11)根据个体最好和种群最好更新种群个体的飞行速度和位置,转步骤(4);
(12)输出电力系统机组组合优化方案,求解结束。
CN201310301711.8A 2013-07-18 2013-07-18 考虑爬坡速率约束的电力系统机组组合优化方法 Active CN103345663B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310301711.8A CN103345663B (zh) 2013-07-18 2013-07-18 考虑爬坡速率约束的电力系统机组组合优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310301711.8A CN103345663B (zh) 2013-07-18 2013-07-18 考虑爬坡速率约束的电力系统机组组合优化方法

Publications (2)

Publication Number Publication Date
CN103345663A CN103345663A (zh) 2013-10-09
CN103345663B true CN103345663B (zh) 2017-07-04

Family

ID=49280458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310301711.8A Active CN103345663B (zh) 2013-07-18 2013-07-18 考虑爬坡速率约束的电力系统机组组合优化方法

Country Status (1)

Country Link
CN (1) CN103345663B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104036320B (zh) * 2014-05-30 2017-02-08 国家电网公司 一种基于改进粒子群算法的微网系统动态经济调度方法
CN106655153B (zh) * 2016-10-11 2018-11-13 广西大学 二类状态变量的电力系统机组组合模型及其建立方法
CN109390981B (zh) * 2018-11-30 2022-05-27 国家电网公司西南分部 新能源参与电量平衡的送端电网机组组合运行控制方法
CN109636046A (zh) * 2018-12-17 2019-04-16 广东电网有限责任公司 一种基于天牛须算法的智能经济调度方法及设备
CN112491049B (zh) * 2020-11-26 2022-10-11 贵州电网有限责任公司 一种考虑在线煤耗曲线的多能源接入电网优化调度方法
CN113032991B (zh) * 2021-03-22 2022-03-18 华中科技大学 一种热电联产机组模型的建模方法及装置
CN114135983B (zh) * 2021-11-19 2023-05-02 国网山东省电力公司电力科学研究院 一种水冷机组与风冷热泵机组的联合优化调度方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101876988A (zh) * 2009-11-25 2010-11-03 天津工业大学 电力系统经济负荷分配的混沌搜索方法
CN102684190B (zh) * 2012-05-25 2014-07-09 浙江大学 一种高效求解含有交流潮流的机组组合方法
CN102789542B (zh) * 2012-06-27 2015-03-04 西安电子科技大学 基于多目标分布式估计的水火电系统负荷分配方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
粒子群优化算法及其应用研究;李丹;《万方数据》;20110630;1-128 *

Also Published As

Publication number Publication date
CN103345663A (zh) 2013-10-09

Similar Documents

Publication Publication Date Title
CN103345663B (zh) 考虑爬坡速率约束的电力系统机组组合优化方法
CN103971181B (zh) 一种用于虚拟电厂的日前经济调度方法
CN112865190A (zh) 计及光伏和充电需求的光储充电站优化调度方法和系统
CN104795828A (zh) 基于遗传算法的风储容量配置方法
CN113780688B (zh) 一种电热联合系统的优化运行方法、系统、设备及介质
CN103440535A (zh) 基于免疫优化和模糊决策的多目标厂级负荷优化方法
CN111553750A (zh) 一种计及电价不确定性和损耗成本的储能竞价策略方法
CN114696351A (zh) 一种电池储能系统动态优化方法、装置、电子设备和存储介质
CN116247648A (zh) 一种考虑源荷不确定性下微电网能量调度的深度强化学习方法
Juul et al. Strategies for charging electric vehicles in the electricity market
CN106355980A (zh) 一种基于限定记忆极限学习机的电网调节能力预测方法
Liu et al. Multi-agent quantum-inspired deep reinforcement learning for real-time distributed generation control of 100% renewable energy systems
Cui et al. Demand response method considering multiple types of flexible loads in industrial parks
Li et al. Dynamic economic dispatch with spinning reserve constraints considering wind power integration
Abasi et al. A novel metaheuristic approach to solve unit commitment problem in the presence of wind farms
Velamuri et al. Reserve constrained economic dispatch incorporating solar farm using particle swarm optimization
CN114884101B (zh) 一种基于自适应模型控制预测的抽水蓄能调度方法
Zhang et al. Distributed Training and Distributed Execution Based Stackelberg Multi-Agent Reinforcement Learning for EV Charging Scheduling
CN117060386A (zh) 一种基于值分布深度q网络的微电网储能调度优化方法
CN112132379A (zh) 考虑经济性的新能源跨区域消纳评估方法和存储介质
CN115912367A (zh) 一种基于深度强化学习的电力系统运行方式智能生成方法
CN113344316A (zh) 一种电网与气网协调运行调度策略确定方法和系统
CN112182835A (zh) 一种考虑风电不确定性和储能调节的电力系统可靠性评估方法及系统
CN113705067B (zh) 一种微网优化运行策略生成方法、系统、设备及存储介质
Zhuang et al. Online Scheduling of PV and Energy Storage System Based on Deep Reinforcement Learning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant