(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111389352.7 (22)申请日 2021.11.23 (71)申请人 吉林化工学院 地址 132022 吉林省吉林市龙潭区承 德街 45号 (72)发明人 董如意　王通　孙立勋　邢雪　 高兴泉　朱建军　 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 50/06(2012.01) G06Q 10/04(2012.01) G06F 17/18(2006.01) G06F 17/16(2006.01) (54)发明名称 一种基于改进蜜獾优化的电力经济排放调 度方法 (57)摘要 本发明涉及一种改进蜜獾优化的电力经济 排放调度方法， 包括以下步骤： S1.采用线 性加权 法将燃料成本和污染排放的多目标优化问题转 化为单目标函数， 建立电力经济排放调度模型； S2.随机初始化发电机组有功出力调度矩阵； S3. 随机初始化另一调度方案； S4.根据牛顿迭代法， 由潮流等式求出两个方案中最后一组发电机有 功出力； S5.计算两个方案中最后一组发电机有 功出力超出电力约束的罚函数值； S6.计算调度 方案中发电机组的燃料成本和污染排放之和， 加 上罚函数值的适应度值； S7.利用改进蜜獾优化 更新发电机组有功出力； S8.若达到最大迭代次 数， 则输出最优发电机组出力； 否则转S3 。 权利要求书2页 说明书5页 附图1页 CN 114066280 A 2022.02.18 CN 114066280 A 1.一种基于改进 蜜獾优化的电力经济排 放调度方法， 其特 征在于， 包括以下步骤： S1.采用线性加权法将燃料成本和污染排放两个目标函数转化为单目标优化函数， 并 结合电力系统潮流等式约束和发电机组有功出力约束， 共同建立电力系统经济排放调 度模 型； S2.设有N组发电机组有功出力需要调度， 随机初始化M个调度方案， 根据公式 （1） 计算 N‑1组有功出力Pi， i=1,2,…N‑1， 共同构成发电机组有功出力方案矩阵A； ， （1） 其中Pmax和Pmin分别为第i组发电机组有功出力的上 下限， r1为[0,1]的随机数； S3.随机初始化另一个调度方案， N ‑1组有功出力为Pi’， i=1,2,…N‑1， 形成新的发电组 有功出力调度方案矩阵A ’， 如公式 （2） 所示； , （2） 其中Pmax和Pmin分别为另一组发电机组有功出力的上 下限， r2为[0,1]的随机数； S4.根据牛 顿迭代法， 由潮流 等式约束求出A与A ’第N组发电机有功出力PN和PN’； S5.计算第N组发电机有功出力PN和PN’超出发电机组约束范围的罚函数值； S6.计算调度方案A与A ’中N组发电机组的燃料成本和污染排放之和， 并与罚函数值共 同构成适应度函数值yi与yi’； S7.比较yi和yi’， 取最小值作为 ybest， 对应的发电机组有功出力为Pbest， 执行改进蜜獾优 化更新发电机组有功出力Pi； S8.若达到最大迭代次数tmax， 则输出最优发电机组出力Pgbest； 否则转S3 。 2.如权利1要求所述的方法， 其特 征在于， 所述 步骤S7具体包括： 步骤S7.1:计算自适应步长 ， 如公式 （3） 所示； （3） 其中， 为调节因子， 通常为0.5； t为当前迭代次数， tmax为最大迭代次数； 步骤S7.2： 计算强度I， 如公式 （4） 所示； （4） 其中， r3为[0,1]的随机数； 步骤S7.3： 更新当前迭代的发电机组有功出力Pnew， 如公式 （5） 所示；权　利　要　求　书 1/2 页 2 CN 114066280 A 2（5） 其中， r4、 r5、 r6、 r7、 r8为[0,1]的随机数； β 为成功系数， 通常取6； 步骤S7.4： 计算Pnew对应的适应度值ynew，更新当前发电机组最优有功出力Pbest和下一次 循环的Pi， 如公式 （6） 所示； （6） 。权　利　要　求　书 2/2 页 3 CN 114066280 A 3