(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111389397.4 (22)申请日 2021.11.23 (71)申请人 吉林化工学院 地址 132022 吉林省吉林市龙潭区承 德街 45号 (72)发明人 董如意　马龙　邢雪　石聪　 刘亚男　 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 10/04(2012.01) G06Q 50/06(2012.01) (54)发明名称 一种基于混合核搜索和黏菌优化的电力经 济排放调度方法 (57)摘要 本发明涉及一种混合的核搜索和黏菌算法 优化的电力经济排放调度方法， 包括以下步骤： S1.采用权重加和法将燃 料成本和污染排放转化 为单目标函数， 建立电力经济排放调度模型； S2. 随机初始化发电机组有功出力调度矩阵； S3.根 据牛顿迭代法， 由潮流等式求出最后一组发电机 有功出力； S4.计算最后一组发电机有功出力超 出约束的罚函数值； S5.计算调度方案中发电机 组的燃料成本和污染排放之和， 加上罚函数值共 同归一化； S6.执行混合核搜索和黏菌优化算法 更新发电机组有功出力； S7.若达到最大迭代次 数， 则输出最优发电机组出力； 否则转S3 。 权利要求书2页 说明书6页 附图1页 CN 114066282 A 2022.02.18 CN 114066282 A 1.一种基于混合核搜索和黏菌优化的电力经济排放调度方法， 其特征在于， 包括以下 步骤： S1.采用权重加和法将燃料成本和污染排放两个目标函数转化为单目标优化函数， 并 结合电力系统潮流等式约束和发电机组有功出力约束， 共同建立电力系统经济排放调 度模 型； S2.设有N组发电机组有功出力需要调度， 随机初始化 M个调度方案， 根据公式 （1） 计算 N‑1组有功出力 Pi，i=1,2,…N‑1， 共同构成发电机组有功出力方案矩阵 A； ， （1） 其中Pmin和Pmax分别为第i组发电机组有功出力的上 下限， rand为[0,1]的随机数； S3.根据牛 顿迭代法， 由潮流 等式约束求出 A中第N组发电机有功出力 PN； S4.计算第 N组发电机有功出力 PN超出发电机组约束范围的罚函数值； S5.计算调度方案 A中N组发电机组的燃料成本和污染排放之和， 并与罚函数值共同构 成适应度函数值 yi， 然后归一 化至[0,1]； S6.执行混合核搜索和黏菌优化方法更新发电机组有功出力 Pi； S7.若达到最大迭代次数， 则输出最优发电机组出力 Pgbest； 否则转S3 。 2.如权利1要求所述的方法， 其特 征在于， 所述 步骤S6具体包括： 步骤S6.1： 将适应度 yi从小到大排序， 形 成新的适应度序列 y， 其最小适应度值为 ybest， 最大适应度值 为yworst； 步骤S6.2： 计算 位置更新速度 ， 如公式 （2） 所示；  (2) 其中，t为当前迭代次数， M为总迭代次数； 步骤S6.3： 计算 黏菌优化 算法的权 重 ， 如公式 （3） 所示；    (3) 其中， rand表示[0,1]的随机数,  median(y)表示 适应度y的中位数； 步骤S6.4： 计算 N‑1组有功出力 ， 如公式(4)所示；  (4) 其中， 为最优发电机组有功出力， 和 为随机挑选的两个 体； 步骤S6.5： 依次将发电机组有功出力 P’i替换A中每个调度方案相应的发电机组有功出 力Pi， 形成新的 M* （N‑1） 个调度方案矩阵 A’， 如公式(5)所示；权　利　要　求　书 1/2 页 2 CN 114066282 A 2(5) 步骤S6.6： 根据牛 顿迭代法， 由潮流 等式约束求出 A’第N组发电机有功出力 P’N； 步骤S6.7： 计算第 N组发电机有功出力 P’N超出发电机组约束范围的罚函数值； 步骤S6.8： 计算调度方案 A’中N组发电机组的燃料成本和污染排放之和， 并与罚函数值 共同构成适应度函数值 y’i； 步骤S6.9： 计算核向量 a， 具体公式如 （6） 所示； （6） 步骤S6.10： 计算近似最优调度方案 Pbest， 具体公式如 （7） 所示； （7） 步骤S6.1 1： 计算新的发电机组有功出力Pnew， 具体公式如 （8） 所示； （8） 其中， rand为[0,1]的随机数。权　利　要　求　书 2/2 页 3 CN 114066282 A 3