(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111483569.4 (22)申请日 2021.12.07 (71)申请人 沈阳理工大 学 地址 110158 辽宁省沈阳市 浑南新区南屏 中路6号 申请人 沈阳工程学院 (72)发明人 张铁岩　王启民　赵琰　姜河　 李兆滢　宋世巍　王东来　李昱材　 王健　林盛　杨君宝　魏莫杋　 (74)专利代理 机构 沈阳之华 益专利事务所有限 公司 21218 代理人 韩凌宇 (51)Int.Cl. G06F 30/20(2020.01) G06Q 10/04(2012.01)G06Q 10/06(2012.01) G06Q 50/06(2012.01) G06F 111/06(2020.01) (54)发明名称 一种含复合储能装置的冷热电联供系统优 化方法 (57)摘要 本发明公开了一种含复合储能装置的冷热 电联供系统优化方法， 包括以下步骤： 建立含有 复合储能装置的冷热电联供系统模 型； 确定含复 合储能装置的冷热电联供系统的多目标优化模 型和约束条件； 使用蜻蜓优化算法进行含复合储 能装置的冷热电联供系统优化调度求解问题。 本 发明采用热电复合的冷热电联供系统并采用了 蜻蜓优化算法有更少的迭代次数、 更准确的全局 搜索能力以及更好的经济效益， 实现了含复合储 能装置的冷热电联供系统的最优求 解。 权利要求书2页 说明书5页 附图4页 CN 114117823 A 2022.03.01 CN 114117823 A 1.一种含复合储能装置的冷热电联 供系统优化方法， 其特 征在于： 包括以下步骤： 步骤1： 建立含有复合储能装置的冷热电联 供系统模型； 步骤2： 确定含复合储能装置的冷热电联 供系统的多目标优化模型和约束条件； 步骤3： 使用蜻蜓优化 算法进行含复合储能装置的冷热电联 供系统优化调度求 解问题。 2.根据权利要求1所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于： 步骤1所述的含有复合储能装置的冷热电联供系统包括微型燃气轮机、 光伏系统、 燃气 锅炉以及含电热的复合储能系统。 3.根据权利要求2所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于： 所述微型燃气轮机建模为： PMT＝ ηMTPMTN 其中， PMT为燃气轮机发出的电功率， ηMT为燃气轮机的发电效率， PMTN为燃气轮机额定电 功率； 其中， QMT为燃气轮机的余热， ηi为燃气轮机的散热衰减系数， PMT为燃气轮机发出的电功 率， ξ 为常数， ηMT为燃气轮机的发电效率。 4.根据权利要求2所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于： 所述光伏系统建模为： 其中， PPV为光伏系统输出的功率， K为光伏系统的功率温度系数， T为光伏系统电池外表 的温度， 为标准测试环境下的温度， 为标准测试环境下光伏系统发出的功率， G为 太阳对光伏系统的光照强度， 为标准测试环境下太阳对光伏系统的光照强度。 5.根据权利要求2所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于： 所述燃气锅炉建模为： 其中， QGB为燃气锅炉的燃气量， PGB为燃气锅炉发出的热功率， ηGB为燃气锅炉的热效率， θGB为常数热值。 6.根据权利要求2所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于： 所述含电热的复合储能系统包括储电系统和储热系统。 7.根据权利要求1所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于步骤2所述的含复合储能装置的冷热电联供系统的多目标优化模型的总运行成本用下式 表示： 权　利　要　求　书 1/2 页 2 CN 114117823 A 2其中， FC为总运行成本， τ为系数， Fy(t)为第t时间段内人工成本， Fq(t)为第t时间段内 燃烧废料处理成本， Fp(t)为复合储能装置中储电系统和储热系统复合成本， t在1到T的时 间范围内取值， Fa为第a种燃料的燃烧成本， a在1到N之间取值； 含复合储能装置的冷热电联 供系统总体的燃料消耗 量如下式： 其中， 为总体燃料消耗量， PMT为燃气轮机发出的电功率， ωMT为燃气轮机的燃料消耗 量， QGB为燃气锅炉的燃气量， ωGB为燃气锅炉发燃料消耗量， Pe为电力系统总体的燃料消耗 量。 8.根据权利要求1所述的一种含复合储能装置的冷热电联供系统优化方法， 其特征在 于步骤3所述的蜻蜓优化 算法包括以下步骤： (1)通过对分离度Ai、 对齐度Bi、 聚集度Ci、 吸引度Ei以及排斥度Pi的组合， 来描述多目标 函数优化问题， 引入步长来代 表方向： ΔXt+1＝ μΔXt+( α Ai+β Bi+χCi+γEi+λPi) 其中， ΔXt+1为更新后的步长， ΔXt为更新前的步长， μ、 α 、 β 、 χ、 γ、 λ为分离度因子； (2)若在设定的Ω范围内， 存在两个变量的欧氏距离小于Ω， 则更新步长直到满足设定 的终止条件时输出程序： Xt+1＝Xt+Levy(l)Xt+1 其中， Levy(l)为飞行 方程。权　利　要　求　书 2/2 页 3 CN 114117823 A 3