本发明涉及微电网领域，尤其是涉及考虑冷热电联供的微电网优化调度方法。
背景技术：
冷热电联供(combinedcoolingheatingandpower，cchp)系统存在冷、热、电三种能量，能够提高能源的综合利用率并减少污染物排放，具有良好的社会和经济效益。微电网解决了分布式电源大规模接入电网的问题，能够有效、灵活地利用各种分布式能源。cchp型微电网存在冷、热、电三种能量之间的平衡关系，而微电网具有并网和孤岛运行两种方式。基于能源利用率和经济效益的考虑，在满足系统冷热负荷的前提下，如何制定最佳的优化调度策略这一问题受到越来越多的关注。
现有的微电网系统中风电、光伏发出功率存在间歇性、随机性等问题。蒸汽轮机在发出电能的同时也会产生大量热能及污染气体。微网的能源综合利用率较低，风电、光伏的消纳能力较差。现有微网的需求响应往往通过负荷削减或转移参与需求响应，会影响用户的用电体验，影响其舒适度。
技术实现要素：
本发明的目的就是为了克服上述现有技术存在微电网运行中的波动性和不稳定性的缺陷而提供一种考虑冷热电联供的微电网优化调度方法。
本发明的目的可以通过以下技术方案来实现：
一种考虑冷热电联供的微电网优化调度方法，用于优化冷热电联供的微电网的运行，所述微电网包括风机、光伏发电机和微型燃气轮机，所述方法采用改进蝴蝶算法求解预建立的微电网优化模型，获取最优解，所述微电网中设有蓄电池，用以在风机和光伏发电不足时，为微电网供电，所述改进蝴蝶算法通过反向学习策略对精英蝴蝶种群进行优化，通过混沌函数对随机数进行优化。
进一步地，所述微电网优化模型的目标函数的表达式为：
f(x)＝jo(x) jf(x) jb(x)
式中，jo(x)为设备的运行维护成本，jf(x)为燃气轮机的燃料成本，jb(x)为蓄电池的折旧成本。
进一步地，所述微电网优化模型包括电储能模型，该电储能模型的表达式为：
式中，ees(t)为t时刻蓄电池的储能容量，τ为蓄电池的自放电率，pes_ch(t)为蓄电池在t时刻的充电功率，pes_dis(t)为蓄电池在t时刻的放电功率，ηsch为蓄电池在t时刻的充电效率，ηsdis为蓄电池在t时刻的放电效率。
进一步地，所述微电网优化模型包括蓄电池约束，该蓄电池约束的表达式为：
soci 1-soci＞δ
式中，pes_ch_max为蓄电池最大允许充电功率，pes_dis_max为蓄电池最大允许放电功率，pes为蓄电池的功率。soci 1为i 1时刻蓄电池的soc值，soci为i 1时刻蓄电池的soc值，δ为两相邻时刻蓄电池的标准soc值。
进一步地，所述改进蝴蝶算法包括以下步骤：
s1：初始化输入参数，载入微电网优化模型，所述输入参数包括风光发电功率、负荷量、微型燃气轮机、蓄电池的出电量以及微电网运行成本函数；
s2：根据反向学习策略，计算出反向解蝴蝶；
s3：计算精英蝴蝶种群中每只蝴蝶产生的香味浓度，根据原始蝴蝶和反向解蝴蝶，并舍弃香味浓度低于预设的第一阈值的蝴蝶，构建精英蝴蝶种群；
s4：计算步骤s3更新后的蝴蝶种群的切换概率，并获取一随机数，若切换函数大于所述随机数，则依次进行步骤s5和s6，否则进行步骤s7，所述随机数的取值范围在0至1之间，所述随机数根据混沌函数获取；
s5：进行全局位置搜索，获取全局最优解；
s6：对全局最优解进行柯西变异，然后进行步骤s8；
s7：进行局部位置搜索，获取局部最优解，然后进行步骤s8；
s8：更新局部最优解和全局最优解；
s9：对步骤s3更新后的精英蝴蝶种群中所有的蝴蝶个体，进行正余弦操作，优化蝴蝶个体进行位置更新；
s10：判断是否达到最大迭代次数，若达到，则结束迭代，输出最优值和最优解，否则执行步骤s3至s10。
进一步地，所述反向解蝴蝶的计算表达式为：
x′ij＝m(aij bij)-xij
式中，xij为普通蝴蝶xi在j维上的值，x′ij为普通蝴蝶的反向解，精英反向系数m为(0，1)的随机数，aij和bij为x′i在j维的最大最小值，j为算法空间的纬度。
进一步地，所述步骤s3中，所述构建精英蝴蝶种群还包括对不在精英蝴蝶种群区间内的蝴蝶赋值，更新精英蝴蝶种群。
首先引入反向学习策略优化蝴蝶算法，将当前迭代的最优个体视为精英个体，通过利用精英个体反向学习生成精英反向解，从当前解和精英反向解中挑选优异个体作为下一代种群，扩大搜索区域的范围，在保证种群多样性的同时，可以增强算法的全局搜索能力和寻优精度。
进一步地，所述对不在初始精英蝴蝶种群区间内的蝴蝶赋值的表达式为：
x′ij＝aij；x′ij＞aij
x′ij＝bij；xij′＜aij
式中，xij为普通蝴蝶xi在j维上的值，
进一步地，步骤s4中，所述随机数的计算表达式为：
rk 1＝sin(πrk),r0＝0.7
式中，rk为第k次迭代时的随机数，r0为随机数的初始值。
进一步地，步骤s6中，所述对全局最优解进行柯西变异的计算表达式为：
xnewbest＝xbest xbest×cauchy(0,1)
式中，xnewbest为柯西变异后的全局最优解，xbest为全局最优解。
与现有技术相比，本发明具有以下优点：
(1)本发明针对楼宇型冷热电联供微电网，利用发电侧的能源互补，以系统运行成本最低为优化目标，优先对微网系统中的风电、光伏进行消纳，提高系统的综合能源利用率，并利用改进的蝴蝶算法进行求解，通过反向学习策略和混沌函数，加快了算法的收敛速度，提高了算法求解的可靠性。
(2)本发明微电网优化模型充分考虑了微电网中风机、光伏的出力，引入cchp系统，提高了微电网对风机、光伏的消纳能力，充分考虑cchp系统中，冷热负荷的潜在储能，利用多能互补提高了微电网系统的综合能源利用率，降低了系统中各设备的储能容量，能有效减少微网的运行成本，提高微网内能源的利用率，减少环境污染。
(3)本发明将混沌算法和正余弦算法引入传统蝴蝶算法，利用柯西变异提高了种群的多样性，通过反向学习策略对精英蝴蝶种群进行优化，提高了蝴蝶算法的全局搜索能力，加快了算法的收敛速度，提高收敛的准确性。
附图说明
图1为含冷热电联供的微电网系统的结构示意图；
图2为本发明改进的蝴蝶算法的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
针对目前微电网运行中的波动性和不稳定性等问题，本实施例拟从发电侧出发，主要考虑冷热电联供系统中存在的冷热负荷的惯性因素，提高微电网对风电和光伏的消纳能力和微电网运行的经济性。cchp的能源利用率较高，能够提高微电网对风能和太阳能的消纳能力。cchp除自身配置的储能设备以外，利用其冷热负荷中的惯性因素构建多元储能系统，进一步优化系统设备的容量配置，提高整体的经济效益。
本实施例提供一种考虑冷热电联供的微电网优化调度方法，应用在含冷热电联供的微电网中。针对楼宇型冷热电联供微电网，利用发电侧的能源互补，以系统运行成本最低为优化目标，优先对微网系统中的风电、光伏进行消纳。重点考虑风机、光伏和燃气轮机的运维成本，建立单目标的微电网优化模型。将cchp系统引入微电网，充分考虑cchp系统中冷热负荷的潜在储能，提高系统的综合能源利用率。采用改进的蝴蝶算法进行求解，求得微网中各设备的最佳运行方式。
下面对本实施例提供的方法进行具体介绍。
一、微电网优化模型
微电网优化模型包括各设备的模型、目标函数以及约束条件，下面进行具体描述。
1、各设备的模型
1.1、微型燃气轮机模型
微型燃气轮机时cchp系统的核心，在发电的同时会产生高温废烟，高温费烟进入余热锅炉转化为热能使用。本发明假设废烟在转化过程中温度不发生变化，忽略外界环境对微型蒸汽轮机的影响。其数学模型为：
其中，
式中cmt(t)为t时段燃气轮机的燃料成本；pmt(t)、ηmt(t)为t时段微型气轮机发出有功功率和发电效率系数；δτ为微型燃气轮机运行时间；rgas为天然气单位价格；lh为天然气低位热值。
1.2、余热锅炉模型
式中hhe(t)为t时段余热锅炉的制热功率输出；η1为燃气轮机散热损失系数；ηhe为热回收效率；cophe为余热锅炉的能效比。x(t)、1-x(t)分别为t时段微型燃气轮机烟气通入余热锅炉的流量比和烟气通入吸收式制冷机的流量比。
1.3、吸收式制冷机模型
式中qac(t)、hac(t)为t时段吸收式制冷机的制冷功率输出和吸收的热功率；copac为吸收式制冷机的能效比；qac_mt(t)为t时段微型燃气轮机余热通过吸收式制冷机的制冷功率输出。
1.4、电制冷机模型
qec(t)＝pec(t)×copec
式中qec(t)、pec(t)分别为t时段电制冷机的制冷功率输出和消耗的电功率；copec为电制冷机的能效比。
1.5、电锅炉模型
heh(t)＝peh(t)×copeh
式中heh(t)、peh(t)为t时段电锅炉的制热功率输出和消耗的电功率；copeh为电锅炉的能效比。
1.6、电储能模型
电储能系统选取蓄电池作为储能设备。其储能容量与蓄电池充放电功率应满足：
式中ees(t)为t时刻蓄电池的储能容量；τ为蓄电池的自放电率；pes_ch(t)、pes_dis(t)和ηsch、ηsdis分别为蓄电池在t时刻的充/放电功率和充/放电效率。
2、确定系统优化目标
2.1、优化目标
选取24h内，蓄电池和微型燃气轮机在每小时内的出力，为cchp型微电网优化调度问题的优化变量。在满足各单元物理约束和系统运行约束的条件下，运行期间优先使用风机、光伏、微型燃气轮机的发电满足电负荷，使系统的运行成本最低。
2.2、目标函数
本实施例中系统的目标函数为系统的运行成本f(x)，其计算表达式为：
f(x)＝jo(x) jf(x) jb(x)
式中:jo(x)为设备的运行维护成本；jf(x)为燃气轮机的燃料成本；jb(x)为蓄电池的折旧成本；目标函数的表达式可表示为：f(x)＝min(f(x))
3、约束条件
3.1、蓄电池
蓄电池应满足充放电功率约束，即：
-pes_ch_max≤pes≤pes_dis_max
式中，pes_ch_max为蓄电池最大允许充电功率，pes_dis_max为蓄电池最大允许放电功率，pes为蓄电池的功率。
通过等价转换可用蓄电池荷电状态表示充放电功率，即
soci 1-soci＞δ
式中，soci 1和soci分别为两相邻时刻的蓄电池soc值。在不同运行状态下(充电或放电状态)，δ值不同。
3.2、微型燃气轮机
微型燃气轮机的发电功率应满足功率的上下限，即：
pgen-min≤pgen(t)≤pgen-max
式中，pgen_min为发电机的最小启动功率；pgen_max为最大发电功率，一般取额定功率值。
此外，燃气轮机还应满足爬坡率的约束：
式中，pup、pdown分别为发电机爬坡率的上、下限。
3.3、功率和能力平衡
系统运行时，应满足电功率平衡，即：
pload(t)＝pes(t) ppv(t) pwt(t)
式中，pload(t)、ppv(t)、pwt(t)分别为负荷预测功率、光伏和风力预测发电功率；pes(t)为蓄电池充放电功率。
同时还需要考虑热能平衡和冷能平衡，即：
qdac(t)×cophe≥qhe(t)
qdac(t)×copco≥qco(t)
式中，qdac(t)表示t时段微燃机在发电过程中产生的余热；cophe和copco分别表示吸收式机组的制热系数和制冷系数；qhe(t)和qco(t)分别表示t时段热负荷和冷负荷水平，两种负荷通常不同时存在。
二、通过改进的烟花算法对系统运行进行优化
本发明中，以优先使用风电、光伏为调度原则，利用改进的蝴蝶算法对该问题进行求解。取微型燃气轮机、蓄电池、的出电量和风机、光伏的预测功率为优化变量，输入改进蝴蝶算法进行优化。
图2所示为改进蝴蝶算法的流程图。具体步骤如下：
s1：初始化输入参数，包括预测风光发电功率、负荷量、微型燃气轮机、蓄电池的出电量以及微电网运行成本函数、约束条件，建立优化调度模型；
s2：应用反向学习策略，通过公式(1)计算出反向解蝴蝶，比较原始蝴蝶和反向解蝴蝶并舍弃香味浓度较低的蝴蝶；利用公式(2)将不再精英区间内的蝴蝶赋值，得到更好的精英蝴蝶种群；
s3：利用混沌函数代替随机函数，减少选择全局搜索和局部搜索的随机性，在每次迭代优化时，采用公式(3)计算rand(随机数)的值。若切换概率p>rand，则按照公式(4)、(5)进行全局位置更新；
s4：若切换概率p<rand，则按照公式(6)进行局部位置更新；
s5：得到局部和全局最优解后，通过公式(7)对所有的蝴蝶个体进行正余弦操作，进一步优化蝴蝶个体进行位置更新；
s6：判断是否达到最大迭代次数，若满足最大迭代次数，则结束程序，输出最优值和最优解，否则转步骤(3)。
此算法中，以微电网运行成本最低为目标函数，优先使用风机、光伏的出力，对微电网内各个设备单元进行优化调度，在满足约束条件的作用下，实现微电网经济最优。对于传统的微电网而言，引入冷热电联供系统后，降低了微电网运行成本，提高了系统的经济性。与传统蝴蝶算法相比，改进蝴蝶算法加快了收敛速度，提高了全局搜索能力。
此算法中各公式的计算表达式为：
公式(1)：
x＇ij＝m(aij bij)-xij
公式(2)：
x′ij＝aij；x′ij＞aij
x′ij＝bij；xij＇＜aij
式中，xij为普通蝴蝶xi在j维上的值，x＇ij为普通蝴蝶的反向解，精英反向系数m为(0，1)的随机数，aij和bij为x＇i在j维的最大最小值，j为算法空间的纬度。
公式(3)：
rk 1＝sin(πrk),r0＝0.7
公式(4)：
式中，
公式(5)：
xnewbest＝xbest xbest×cauchy(0,1)
公式(6)：
式中，
公式(7)：
上式中含有四个参数，r1、r2、r3、r4，r1决定在下一次迭代第i个个体的位置更新方向，r2是[0,2π]之间的随机数，它决定下一次迭代中个体的移动距离；r3是随机权重的取值范围是[0,2]，r3>1时，
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。
1.一种考虑冷热电联供的微电网优化调度方法，用于优化冷热电联供的微电网的运行，所述微电网包括风机、光伏发电机和微型燃气轮机，其特征在于，所述方法采用改进蝴蝶算法求解预建立的微电网优化模型，获取最优解，所述微电网中设有蓄电池，用以在风机和光伏发电不足时，为微电网供电，所述改进蝴蝶算法通过反向学习策略对精英蝴蝶种群进行优化，通过混沌函数对随机数进行优化。
2.根据权利要求1所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述微电网优化模型的目标函数的表达式为：
f(x)＝jo(x) jf(x) jb(x)
式中，jo(x)为设备的运行维护成本，jf(x)为燃气轮机的燃料成本，jb(x)为蓄电池的折旧成本。
3.根据权利要求1所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述微电网优化模型包括电储能模型，该电储能模型的表达式为：
式中，ees(t)为t时刻蓄电池的储能容量，τ为蓄电池的自放电率，pes_ch(t)为蓄电池在t时刻的充电功率，pes_dis(t)为蓄电池在t时刻的放电功率，ηsch为蓄电池在t时刻的充电效率，ηsdis为蓄电池在t时刻的放电效率。
4.根据权利要求1所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述微电网优化模型包括蓄电池约束，该蓄电池约束的表达式为：
soci 1-soci＞δ
式中，pes_ch_max为蓄电池最大允许充电功率，pes_dis_max为蓄电池最大允许放电功率，pes为蓄电池的功率。soci 1为i 1时刻蓄电池的soc值，soci为i 1时刻蓄电池的soc值，δ为两相邻时刻蓄电池的标准soc值。
5.根据权利要求1所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述改进蝴蝶算法包括以下步骤：
s1：初始化输入参数，载入微电网优化模型，所述输入参数包括风光发电功率、负荷量、微型燃气轮机、蓄电池的出电量以及微电网运行成本函数；
s2：根据反向学习策略，计算出反向解蝴蝶；
s3：计算精英蝴蝶种群中每只蝴蝶产生的香味浓度，根据原始蝴蝶和反向解蝴蝶，并舍弃香味浓度低于预设的第一阈值的蝴蝶，构建精英蝴蝶种群；
s4：计算步骤s3更新后的蝴蝶种群的切换概率，并获取一随机数，若切换函数大于所述随机数，则依次进行步骤s5和s6，否则进行步骤s7，所述随机数的取值范围在0至1之间，所述随机数根据混沌函数获取；
s5：进行全局位置搜索，获取全局最优解；
s6：对全局最优解进行柯西变异，然后进行步骤s8；
s7：进行局部位置搜索，获取局部最优解，然后进行步骤s8；
s8：更新局部最优解和全局最优解；
s9：对步骤s3更新后的精英蝴蝶种群中所有的蝴蝶个体，进行正余弦操作，优化蝴蝶个体进行位置更新；
s10：判断是否达到最大迭代次数，若达到，则结束迭代，输出最优值和最优解，否则执行步骤s3至s10。
6.根据权利要求5所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述反向解蝴蝶的计算表达式为：
x′ij＝m(aij bij)-xij
式中，xij为普通蝴蝶xi在j维上的值，x′ij为普通蝴蝶的反向解，精英反向系数m为(0，1)的随机数，aij和bij为x′i在j维的最大最小值，j为算法空间的纬度。
7.根据权利要求5所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述步骤s3中，所述构建精英蝴蝶种群还包括对不在精英蝴蝶种群区间内的蝴蝶赋值，更新精英蝴蝶种群。
8.根据权利要求7所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，所述对不在初始精英蝴蝶种群区间内的蝴蝶赋值的表达式为：
x′ij＝aij；x′ij＞aij
x′ij＝bij；xij′＜aij
式中，xij为普通蝴蝶xi在j维上的值，x′ij为普通蝴蝶的反向解，精英反向系数m为(0，1)的随机数，aij和bij为x′i在j维的最大最小值，j为算法空间的纬度。
9.根据权利要求5所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，步骤s4中，所述随机数的计算表达式为：
rk 1＝sin(πrk),r0＝0.7
式中，rk为第k次迭代时的随机数，r0为随机数的初始值。
10.根据权利要求5所述的一种考虑冷热电联供的微电网优化调度方法，其特征在于，步骤s6中，所述对全局最优解进行柯西变异的计算表达式为：
xnewbest＝xbest xbest×cauchy(0,1)
式中，xnewbest为柯西变异后的全局最优解，xbest为全局最优解。
技术总结