首页 / 专利库 / 燃料种类 / 能源 / 可再生能源 / 包括电动汽车充电负荷的微电网经济调度优化方法

包括电动汽车充电负荷的微电网经济调度优化方法

阅读:665发布:2020-05-11

专利汇可以提供包括电动汽车充电负荷的微电网经济调度优化方法专利检索,专利查询,专利分析的服务。并且为了解决 现有技术 中确定性优化方法所建立的优化模型没有将电动 汽车 入网考虑其中而造成的综合 能源 系统灵活性不够的缺点,本 发明 公开了一种包括电动汽车充电负荷的微 电网 经济调度优化方法,同时考虑电动汽车充电负荷、 风 光出 力 以及CCHP系统冷热电负荷的不确定性和随机性,利用随机规划中的期望值模型来描述充电负荷的不确定性和随机性,同时采用机会约束形式来表征充电设施的投建容量能够满足电动汽车的充电需求,采用多场景随机规划的方法来描述各随机变量的不确定性和随机性。本发明期望目标是提高综合能源系统灵活性、增强 可再生能源 消纳能力,同时降低系统在运行周期内的综合运行成本,为微电网经济调度运行提供必要的技术 支撑 。,下面是包括电动汽车充电负荷的微电网经济调度优化方法专利的具体信息内容。

1.一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:同时考虑电动汽车充电负荷、光出以及CCHP系统冷热电负荷的不确定性和随机性,利用随机规划中的期望值模型来描述充电负荷的不确定性和随机性,同时采用机会约束形式来表征充电设施的投建容量能够满足电动汽车的充电需求,采用多场景随机规划的方法来描述各随机变量的不确定性和随机性。
2.根据权利要求1所述的一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:具体步骤包括:
S1.获取数据集,构建电动汽车充电负荷数据矩阵;
S2.根据S1步骤获得的电动汽车充电负荷数据矩阵构建计及电动汽车充电负荷经济调度模型的目标函数;
S3.构建S2步骤所述计及电动汽车充电负荷经济调度模型的约束条件;
S4.对S3步骤中的约束条件中的机会约束条件进行转化和处理;
S5.对S2步骤中的计及电动汽车充电负荷经济调度模型的进行求解;
S5步骤中具体的过程是:
S501.读取S1步骤所述的系统数据,设置S3步骤中机会约束条件下的置信平,输入改进万有引力搜索算法所需参数;
S502.获取S3步骤所述的目标函数和约束条件;
S503.利用蒙特卡洛方法和后向缩减法进行场景生成和削减,并初始化各随机变量在各场景下相应的发生概率;
S503.采用混沌方法初始化粒子,检验粒子是否在取值范围内:随机产生一个n维、每个分量数值在0~1之间的向量Z1=[Z11,Z12,…,Z1n],利用混沌算法得出Z2、Z3,…,ZM;
设n维优化变量ai≤xi≤bi,将Zi的各个分量载波到优化变量的取值范围得:
其中,xi为第i维优化变量,ai、bi为常数,对应第i维优化变量的上下
限;
S505.计算适应度及粒子受到的引力和加速度:根据大数定律,检验优化变量是否满足约束条件,对于满足约束条件的粒子,以目标函数值作为其适应度值;对于不满足约束条件的粒子,采用惩罚函数方法计算其适应度值;其中, 表示t时刻i粒子的更新加速度,可表示为:
式中:Fik(t)表示在t时刻第k维粒子的适应度值,mi(t)表示t时刻粒子的质量;其中,表示在t时刻第k维粒子a受到粒子b的引力,可表示为:
式中:k1、k2是与万有引力计算公式相关的两个常数;||Na(t),Nb(t)||2表示粒子Na(t)和Nb(t)之间的欧氏距离; 表示在t时刻第k维粒子a、b的位置
S506.更新个体的位置和速度:从M个初始群体中择优选取m个粒子作为初始解,采用随机模拟方法从中选出群体最优值,利用PSO算法中的速度更新公式改写GSA中的速度更新公式:
式中:ω为PSO中的惯性权重;k3、k4、k5为3个[0,1]中的随机数;c1、c2为PSO中的学习因子; 分别为粒子的个体最优和群体最优;利用载波得到的优化变量取值改写GSA中的位置更新公式:
其中,aj对应第j维优化变量上限、Zij对应Zi
各个分量, 为低k维粒子i在t+1时刻的速度;
然后检验是否满足约束条件,得出适应度值;
S507.对优化变量进行混沌扰动,并更新个体最优值和群体最优值;
S508.检验是否满足终止条件:是,输出结果;否,更新各随机变量在各场景下相应的生成概率,然后转到步骤S505。
3.根据权利要求2所述的一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:S1步骤的具体过程如下:
从微电网能量管理系统MGEMS获取某地某段时间数据集,并利用蒙特卡洛方法生成S个充电负荷场景,将获取到的一天中电动汽车充电桩i在场景s下时段t的充电功率,构建电动汽车充电负荷数据矩阵;
其中PiEV(s,t)表征电动汽车充电桩i在场景s下时段t的充电功率, 为j
型充电桩的规划数目上限,j=1,2分别代表快、慢充电桩;S、T分别为生成的场景总数、每天划分的时段总数;单个电动汽车充电桩i每日充电负荷数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T;
从微电网能量管理系统MGEMS获取某地某段时间历史和实时数据集,并利用蒙特卡洛方法生成S个风机出力场景,将获取到的每日风速数据,通过处理、计算和分析,构建风机在场景s下时段t的出力数据矩阵:
其中 表征第d天在场景s下时段t的风机出力,d为该段时间的天数,d=1,
2,...,;S、T分别为生成的场景总数、每天划分的时段总数;单日风机出力数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
从微电网能量管理系统MGEMS获取某地某段时间历史和实时数据集,并利用蒙特卡洛方法生成S个光伏出力场景,将获取到的每日光照强度数据,构建光伏发电站在场景s下时段t的出力数据矩阵:
其中 表征第d天在场景s下时段t的光伏板出力,d为该时段的天数,d=1,
2,...,;S、T分别为生成的场景总数、每天划分的时段总数;单日光伏板出力数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T;
从微电网能量管理系统MGEMS获取某地某段时间的历史数据和实时数据集,并利用蒙特卡洛方法生成S个电负荷功率场景,将获取到的电负荷功率数据,构建系统在场景s下时段t的电负荷功率数据矩阵:
其中 表征第d天系统在场景s下时段t的电负荷功率,d为该时段的天数,d=1,
2,...,;S、T分别为生成的场景总数、每天划分的时段总数。单日系统电负荷功率数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T;
从微电网能量管理系统MGEMS获取某地某段时间的历史数据和实时数据集,并利用蒙特卡洛方法生成S个冷/热负荷功率场景,将获取到的冷/热负荷功率数据,构建系统在场景s下时段t的冷/热负荷功率数据矩阵:
其中 表征第d天在场景s下时段t的冷/热负荷功率,d为该时段的天数,d=1,
2,...,;j=c,h分别对应冷、热负荷;S、T分别为生成的场景总数、每天划分的时段总数;单日冷/热负荷功率数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
4.根据权利要求2所述的一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:S2步骤的具体过程如下:
综合能源系统中综合运行成本目标函数的数学表达式为:
minCS=min(C1+C2+C3+C4+C5+C6)
经济调度模型中各成本分量具体表述如下:
1)C1表示系统中电动汽车充电基础设施总的规划成本:
式中:C11为充电基础设施的建设投资成本;C12为充电基础设施的运行成本;C13为充电基础设施的维护成本; 为j型充电桩的规划数目上限,j=1,2分别代表快、慢充电桩;Yij是充电桩i是否投建的0-1变量,j=1,2分别代表快、慢充电桩; 是充电桩i的最大充电容量,j=1,2分别代表快、慢充电桩;kij是充电桩i的容量投资单价,j=1,2分别代表快、慢充电桩;R0为贴现率; 是充电桩i的折旧年限,j=1,2分别代表快、慢充电桩;T为时段总数;S为场景总数;ρs是充电负荷场景s发生的概率, 是场景s下时段t充电桩i的充电功率,j=1,2分别代表快、慢充电桩;k(t)为时段t的购电价格; 为每年充电桩i的单位容量维护成本,j=1,2分别代表快、慢充电桩;
2)C2表示系统中常规发电机组即运行成本、启动成本、停机成本的和:
式中:Nc为常规发电机组数量; 分别为常规发电机组的运行成本、
启动成本和停机成本;Δt为时间间隔;
3)C3表示系统中快速调节机组即运行成本、启动成本、停机成本的和:
式中:Nk为快速调节机组数量; 分别为快速调节机组的运行成本、
启动成本和停机成本;
4)C4表示系统中燃气轮机即运行成本、启动成本、停机成本的和:
式中:Nr为燃气轮机数量; 分别为燃气轮机的运行成本、启动成本
和停机成本;
5)C5表示系统中燃气锅炉的运行成本:
式中: 为燃气锅炉的运行成本。
6)C6表示系统与主网之间的购售电成本:
式中: 分别为场景s下时段t联络线购电和售电的功率;kB、kS分别为购电和售电单价。
5.根据权利要求2所述的一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:S3步骤中所述的约束条件包括:
1)充电桩的充电功率约束:
是场景s下时段t充电桩i的充电功率,j=1,2分别代表
快、慢充电桩;Yij是充电桩i是否投建的0-1变量,j=1,2分别代表快、慢充电桩; 是充电桩i的最大充电容量,j=1,2分别代表快、慢充电桩;
2)充电站总容量的机会约束:
式中:kCAP为充电站容量的冗余度系数;Q(s,t)为场景s下时段t电动汽车集群的总充电负荷;Prob{·}为{·}中事件成立的概率,α1为给定约束条件成立的置信水平;
3)充电站充电功率平衡约束:
Q(s,t)为场景s下时段t电动汽车集群的总充电负荷;
4)充电站电量约束:
式中:Tij为充电桩i一次工作的平均充电时间,j=1,2分别代表快、慢充电桩; 为充电桩i全天的最大可充电次数,j=1,2分别代表快、慢充电桩;RENER为充电站电量的冗余度系数;
5)充电设施建设的总数目限制:
式中:nmax为充电站内可以容纳的慢充桩与快充桩的最大数目;
6)系统功率平衡约束:
式中: 分别为常规发电机组、快速调节机组、燃气轮机场景s下时段t
的出力; 分别为风机和光伏场景s下时段t的出力; 分别为场景s下
时段t系统的电、冷、热负荷功率; 为电制冷机场景s下时段t消耗的电功率; 为吸收式制冷机场景s下时段t的制冷功率; 为电制冷机场景s下时段t的制冷功率; 为余热锅炉场景s下时段t提供的热功率; 为燃气锅炉场景s下时段t提供的热功率;
7)系统机组出力上下限约束:
式中: 分别为常规发电机组、快速调节机组、燃气轮机、燃气锅炉
的开停机状态,1为开机,0为停机; 分别为常规发电机组出力的最小值和最大值; 分别为快速调节机组出力最小值和最大值; 分别为燃气轮机出
力的最小值和最大值; 分别为吸收式制冷机和电制冷机的额定功率;
分别为燃气锅炉的最小出力和最大出力; 分别为燃气锅炉用于制热
和供吸收式制冷机制冷的功率;
8)发电机组爬坡率约束:
式中: 分别为常规发电机组、快速调节机组出力增加和减少的爬
坡率; 分别为常规发电机组、快速调节机组、燃气轮机、燃气锅炉的开停机状态;
9)最小开停机时间约束:
式中: 分别为常规发电机组、燃气锅炉的最小开机时间
和最小停机时间; 分别为常规发电机组、快速调节机组、燃气轮机、燃气锅炉的开停机状态;
10)联络线功率约束:
式中:Pex,max为联络线最大传输功率; 为购售电状态,1为购电,0为售电;
分别为场景s下时段t联络线购电和售电的功率。

说明书全文

包括电动汽车充电负荷的微电网经济调度优化方法

技术领域

[0001] 本发明属于电系统及其自动化技术领域,涉及一种计及电动汽车充电负荷的CCHP微电网经济调度优化方法。

背景技术

[0002] 随着科学技术的发展,使得全球工业迅速崛起,经济平不断提高,但随之所带来的能源危机和环境污染问题也日益严重。传统的能源结构及利用方式难以为继,因此,推动能源结构转型、提高能源利用效率、发展清洁能源和可再生能源势在必行。
[0003] 如今,微网研究与应用越来越引起世界各国的重视。微网是通过分布式微型电源,向用户端提供各类负荷的小型供能系统,它有孤岛运行和并网运行两种运行方式。而CCHP型微网可以将电、冷、热、气、可再生能源等多种能源有机耦合,对区域内多种分布式资源进行群管群控,建立多能综合利用平台,整合利用区域内的可再生能源、可控供能机组、冷热电三联供等分布式资源,结合能量的存储及转换技术,实现电、气、冷、热可再生能源等多种能源互补供应、协调运行,提高系统运行效率,降低运行成本。
[0004] 微网经济分配就是根据电价和负荷需求等因素灵活地调节各个发电单元的输出功率,在保证供需平衡的同时,使经济效益达到最高。然而以、光为代表的可再生能源、其出力受自然资源影响较大,具有间歇性、波动性和随机性,易对电力系统的运行控制造成不良影响,因此需要提高能源供给系统灵活性,以增强可再生能源消纳能力。
[0005] 而且近年来,我国持续加大对于电动汽车发展的支持力度,以促进燃油替代,减少尾气排放,保障能源安全、防治大气污染,不断推动我国从汽车大国向汽车强国转变。为了进一步促进电动汽车的普及应用,完善充电基础设施的规划发展,有必要在CCHP系统中考虑电动汽车入网对系统经济分配的影响。
[0006] 针对新能源入网问题,以往通常采用确定性的优化方法进行建模和求解,且建立的CCHP系统大都仅考虑了风光等可再生能源,而很少将电动汽车充电基础设施的规划考虑其中。且由于可再生能源出力及负荷功率具有不确定性,传统的确定性优化方法存在一定的局限性。所以,传统的确定性优化方法所建立的模型大都没有将电动汽车入网考虑其中,更无法确定电动汽车充电负荷对CCHP系统经济分配的影响。

发明内容

[0007] 为了解决现有技术中确定性优化方法所建立的优化模型没有将电动汽车入网考虑其中而造成的综合能源系统灵活性不够的缺点,本发明公开了一种包括电动汽车充电负荷的微电网经济调度优化方法。
[0008] 本发明为了解决上述技术体所采用的技术方案是:一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:同时考虑电动汽车充电负荷、风光出力以及CCHP系统冷热电负荷的不确定性和随机性,利用随机规划中的期望值模型来描述充电负荷的不确定性和随机性,同时采用机会约束形式来表征充电设施的投建容量能够满足电动汽车的充电需求,采用多场景随机规划的方法来描述各随机变量的不确定性和随机性。
[0009] 本发明的有益效果是:本发明同时考虑电动汽车充电负荷、风光出力以及CCHP系统冷热电负荷的不确定性和随机性,通过微网能量管理系统MGEMS获取某地微电网运行的数据,在考虑电动汽车充电负荷的不确定性和随机性时,采用多场景随机规划的方法,利用随机规划中的期望值模型来描述充电负荷的不确定性和随机性,同时采用机会约束形式来表征充电设施的投建容量能够在一定概率下满足电动汽车的充电需求,在考虑风光出力以及CCHP系统冷热电负荷的不确定性和随机性时,同样采用多场景随机规划的方法来描述各随机变量的不确定性和随机性。其期望目标是提高综合能源系统灵活性、增强可再生能源消纳能力,同时降低系统在日、月、年、多年等运行周期内的综合运行成本,为微电网经济调度运行提供必要的技术支撑附图说明
[0010] 图1是本发明所针对的风-光-气-蓄联合调度系统结构图。
[0011] 图2是本发明的流程框图
[0012] 图3是本发明的求解过程流程图

具体实施方式

[0013] 下面参照附图并结合实例对本发明的具体实施方式做进一步详细描述。
[0014] 如图1,一种包括电动汽车充电负荷的微电网经济调度优化方法,其特征在于:同时考虑电动汽车充电负荷、风光出力以及CCHP系统冷热电负荷的不确定性和随机性,利用随机规划中的期望值模型来描述充电负荷的不确定性和随机性,同时采用机会约束形式来表征充电设施的投建容量能够满足电动汽车的充电需求,采用多场景随机规划的方法来描述各随机变量的不确定性和随机性。
[0015] 图2中的步骤1描述构建数据矩阵的过程和方法。从微电网能量管理系统MGEMS获取某地一年历史和实时数据集,并利用蒙特卡洛方法生成S个充电负荷场景,将获取到的一天中电动汽车充电桩i在场景s下时段t的充电功率,通过处理、计算和分析,构建电动汽车充电负荷数据矩阵:
[0016]
[0017] 其中PiEV(s,t)表征电动汽车充电桩i在场景s下时段t的充电功率, (为j型充电桩的规划数目上限,j=1,2分别代表快、慢充电桩);S、T分别为生成的场景总数、每天划分的时段总数。单个电动汽车充电桩i每日充电负荷数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
[0018] 从微电网能量管理系统MGEMS获取某地一年历史和实时数据集,并利用蒙特卡洛方法生成S个风机出力场景,将获取到的每日风速数据,通过处理、计算和分析,构建风机在场景s下时段t的出力数据矩阵:
[0019]
[0020] 其中 表征第d天在场景s下时段t的风机出力,d=1,2,...,365(对应一年当中相应的日期);S、T分别为生成的场景总数、每天划分的时段总数。单日风机出力数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
[0021] 从微电网能量管理系统MGEMS获取某地一年历史和实时数据集,并利用蒙特卡洛方法生成S个光伏出力场景,将获取到的每日光照强度数据,通过处理、计算和分析,构建光伏发电站在场景s下时段t的出力数据矩阵:
[0022]
[0023] 其中 表征第d天在场景s下时段t的光伏板出力,d=1,2,...,365(对应一年当中相应的日期);S、T分别为生成的场景总数、每天划分的时段总数。单日光伏板出力数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
[0024] 从微电网能量管理系统MGEMS获取某地一年历史和实时数据集,并利用蒙特卡洛方法生成S个电负荷功率场景,将获取到的电负荷功率数据,通过处理、计算和分析,构建系统在场景s下时段t的电负荷功率数据矩阵:
[0025]
[0026] 其中 表征第d天系统在场景s下时段t的电负荷功率,d=1,2,...,365(对应一年当中相应的日期);S、T分别为生成的场景总数、每天划分的时段总数。单日系统电负荷功率数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
[0027] 从微电网能量管理系统MGEMS获取某地一年历史和实时数据集,并利用蒙特卡洛方法生成S个冷/热负荷功率场景,将获取到的冷/热负荷功率数据,通过处理、计算和分析,构建系统在场景s下时段t的冷/热负荷功率数据矩阵:
[0028]
[0029] 其中 表征第d天在场景s下时段t的冷/热负荷功率,d=1,2,...,365(对应一年当中相应的日期);j=c,h分别对应冷、热负荷;S、T分别为生成的场景总数、每天划分的时段总数。单日冷/热负荷功率数据样本的数量N=S*T,数据集的数量n=T,s=1,2,...,S,t=1,2,...,T。
[0030] 图2中的步骤2描述计及电动汽车充电负荷经济调度模型的目标函数构建的过程和方法。综合能源系统中综合运行成本要满足各场景下系统中电动汽车充电基础设施总的规划成本、常规发电机组(运行成本+启动成本+停机成本)、快速调节机组(运行成本+启动成本+停机成本)、燃气轮机(运行成本+启动成本+停机成本)、燃气锅炉的运行成本、与主网之间的购售电成本六者期望之和最小化,其目标函数的数学表达式为:
[0031] minCS=min(C1+C2+C3+C4+C5+C6)
[0032] 经济分配模型中各成本分量具体表述如下:
[0033] 1)C1表示系统中电动汽车充电基础设施总的规划成本,其具体数学表达式为:
[0034] C1=C11+C12+C13
[0035]
[0036] 式中:C11为充电基础设施的建设投资成本;C12为充电基础设施的运行成本;C13为充电基础设施的维护成本; 为j型充电桩的规划数目上限,j=1,2分别代表快、慢充电桩;Yij是充电桩i是否投建的0-1变量,j=1,2分别代表快、慢充电桩; 是充电桩i的最大充电容量,j=1,2分别代表快、慢充电桩; 是充电桩i的容量投资单价,j=1,2分别代表快、慢充电桩;R0为贴现率; 是充电桩i的折旧年限,j=1,2分别代表快、慢充电桩;T为时段总数;S为场景总数;ρs是充电负荷场景s发生的概率, 是场景s下时段t充电桩i的充电功率,j=1,2分别代表快、慢充电桩;k(t)为时段t的购电价格; 为每年充电桩i的单位容量维护成本,j=1,2分别代表快、慢充电桩。
[0037] 2)C2表示系统中常规发电机组(运行成本+启动成本+停机成本),其具体数学表达式为:
[0038]
[0039] 式中:Nc为常规发电机组数量; 分别为常规发电机组的运行成本、启动成本和停机成本;Δt为时间间隔。
[0040] 3)C3表示系统中快速调节机组(运行成本+启动成本+停机成本),其具体数学表达式为:
[0041]
[0042] 式中:Nk为快速调节机组数量; 分别为快速调节机组的运行成本、启动成本和停机成本。
[0043] 4)C4表示系统中燃气轮机(运行成本+启动成本+停机成本),其具体数学表达式为:
[0044]
[0045] 式中:Nr为燃气轮机数量; 分别为燃气轮机的运行成本、启动成本和停机成本。
[0046] 5)C5表示系统中燃气锅炉的运行成本,其具体数学表达式为:
[0047]
[0048] 式中: 为燃气锅炉的运行成本。
[0049] 6)C6表示系统与主网之间的购售电成本,其具体数学表达式为:
[0050]
[0051] 式中: 分别为场景s下时段t联络线购电和售电的功率;kB、kS分别为购电和售电单价。
[0052] 图2中的步骤3描述计及电动汽车充电负荷经济调度模型的约束条件构建的过程和方法。计及电动汽车充电负荷经济调度模型的约束条件包括:充电桩的充电功率约束、充电站总容量的机会约束、充电站充电功率平衡约束、充电站电量约束、充电设施建设的总数目限制、系统功率平衡约束、系统机组出力上下限约束、发电机组爬坡率约束、最小开停机时间约束、联络线功率约束等。
[0053] 各约束条件具体表述如下:
[0054] 1)充电桩的充电功率约束:
[0055]
[0056] 2)充电站总容量的机会约束:
[0057]
[0058] 式中:kCAP为充电站容量的冗余度系数;Q(s,t)为场景s下时段t电动汽车集群的总充电负荷;Prob{·}为{·}中事件成立的概率,α1为给定约束条件成立的置信水平。
[0059] 3)充电站充电功率平衡约束:
[0060]
[0061] 4)充电站电量约束:
[0062]
[0063] 式中:Tij为充电桩i一次工作的平均充电时间,j=1,2分别代表快、慢充电桩;为充电桩i全天的最大可充电次数,j=1,2分别代表快、慢充电桩;RENER为充电站电量的冗余度系数。
[0064] 5)充电设施建设的总数目限制:
[0065]
[0066] 式中:nmax为充电站内可以容纳的慢充桩与快充桩的最大数目。
[0067] 6)系统功率平衡约束:
[0068]
[0069]
[0070]
[0071] 式中: 分别为常规发电机组、快速调节机组、燃气轮机场景s下时段t的出力; 分别为风机和光伏场景s下时段t的出力; 分别为场
景s下时段t系统的电、冷、热负荷功率; 为电制冷机场景s下时段t消耗的电功率; 为吸收式制冷机场景s下时段t的制冷功率; 为电制冷机场景s下时段t的制冷功率; 为余热锅炉场景s下时段t提供的热功率; 为燃气锅炉场景s下时段t提供的热功率。
[0072] 7)系统机组出力上下限约束:
[0073]
[0074]
[0075]
[0076]
[0077]
[0078] 式中: 分别为常规发电机组、快速调节机组、燃气轮机、燃气锅炉的开停机状态,1为开机,0为停机; 分别为常规发电机组出力的最小值和最大值; 分别为快速调节机组出力最小值和最大值; 分别为燃气
轮机出力的最小值和最大值; 分别为吸收式制冷机和电制冷机的额定功率;
分别为燃气锅炉的最小出力和最大出力; 分别为燃气锅炉用于制热
和供吸收式制冷机制冷的功率。
[0079] 8)发电机组爬坡率约束:
[0080]
[0081]
[0082] 式中: 分别为常规发电机组、快速调节机组出力增加和减少的爬坡率。
[0083] 9)最小开停机时间约束:
[0084]
[0085]
[0086] 式中: 分别为常规发电机组、燃气锅炉的最小开机时间和最小停机时间。
[0087] 10)联络线功率约束:
[0088]
[0089] 式中:Pex,max为联络线最大传输功率; 为购售电状态,1为购电,0为售电。
[0090] 图2中的步骤4描述机会约束条件处理的过程和方法。在一定置信水平下,电动汽车充电站的总容量满足充电负荷的机会约束条件可以等价为如下形式:
[0091]
[0092] 式中:Kα1是一个较大的常数;Φ-1(·)为随机变量KCAP·Q(s,t)的累积概率分布函数。
[0093] 图2中的步骤5描述计及电动汽车充电负荷经济调度模型的求解过程和方法。采用智能算法中的万有引力搜索算法,并利用混沌粒子群算法对引力搜索算法的速度和位置更新公式进行改进,结合各负荷场景s所对应的概率进行优化调整,可以求解出电动汽车总规划成本的最小方案、综合能源系统中各机组或设备的开停机计划、联络线购售电计划使得经济分配模型所得到的解最优。
[0094] 具体步骤如下:
[0095] 1)读取系统数据,设置机会约束条件下的置信水平,输入改进万有引力搜索算法所需参数;
[0096] 2)获取目标函数和约束条件;
[0097] 3)利用蒙特卡洛方法和后向缩减法进行场景生成和削减,并初始化各随机变量在各场景下相应的发生概率;
[0098] 4)采用混沌方法初始化粒子,检验粒子是否在取值范围内:随机产生一个n维、每个分量数值在0~1之间的向量Z1=[Z11,Z12,…,Z1n],利用混沌算法得出Z2、Z3,…,ZM。设n维优化变量ai≤xi≤bi,将Zi的各个分量载波到优化变量的取值范围得:
[0099] 5)计算适应度及粒子受到的引力和加速度:根据大数定律,检验优化变量是否满足约束条件,对于满足约束条件的粒子,以目标函数值作为其适应度值;对于不满足约束条件的粒子,采用惩罚函数方法计算其适应度值;其中, 表示t时刻i粒子的更新加速度,可表示为:
[0100]
[0101] 式中:Fik(t)表示在t时刻第k维粒子的适应度值,mi(t)表示t时刻粒子的质量;其中, 表示在t时刻第k维粒子a受到粒子b的引力,可表示为:
[0102]
[0103]
[0104] 式中:k1、k2是与万有引力计算公式相关的两个常数,且k2非常小;||Na(t),Nb(t)||2表示粒子Na(t)和Nb(t)之间的欧氏距离; 表示在t时刻第k维粒子a、b的位置;
[0105] 6)更新个体的位置和速度:从M个初始群体中择优选取m个粒子作为初始解,它们的位置和目标函数值即为个体最优值,采用随机模拟方法从中选出群体最优值,利用PSO算法中的速度更新公式改写GSA中的速度更新公式:
[0106]
[0107] 式中:ω为PSO中的惯性权重;k3、k4、k5为3个[0,1]中的随机数;c1、c2为PSO中的学习因子; 分别为粒子的个体最优和群体最优;利用载波得到的优化变量取值改写GSA中的位置更新公式:
[0108]
[0109] 然后检验是否满足约束条件,计算出适应度值;
[0110] 7)对优化变量进行混沌扰动,并更新个体最优值和群体最优值;
[0111] 8)检验是否满足终止条件:是,输出结果;否,更新各随机变量在各场景下相应的生成概率,然后转到步骤5)。
[0112] 以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易变化或替换,都属于本发明的保护范围之内。因此本发明的保护范围所述以权利要求的保护范围为准。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈