(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111591693.2 (22)申请日 2021.12.23 (71)申请人 山东理工大 学 地址 255086 山东省淄博市高新 技术产业 开发区高创园A座313室 (72)发明人 陈佳佳　戎泽坤　陈文刚　徐丙垠　 肖传亮　 (74)专利代理 机构 淄博市众朗知识产权代理事 务所(特殊普通 合伙) 37316 代理人 程强强　张宁 (51)Int.Cl. G06Q 30/02(2012.01) G06Q 10/06(2012.01) G06Q 50/06(2012.01) (54)发明名称 一种多能微电网非线性综合需求响应的管 理方法 (57)摘要 一种多能微电网非线性综合需求响应的管 理方法， 属于智 能微电网技术领域。 初始化每小 时负荷需求， 电价、 天然气价格和预测的风机出 力； 确定风机出力， 电力和天然气价格的波动范 围； 初始化进化捕食算法中的种群规模以及最大 迭代次数； 根据最短距离算法确定最佳峰谷平时 段， 获得电力和燃气负荷的最优分时价格， 并计 算电力弹性系数矩阵； 依次执行电与气的幂、 指 数、 对数以及线性需求响应模型， 以获得每小时 电负荷与气负荷的需求响应； 根据可信度理论确 定风机出力的置信度模型； 评估多能微电网中的 约束， 计算目标函数并通过进化捕食算法更新决 策变量。 本发 明确定多能微电网的最优分时价格 及不同时刻最优需求响应量， 降低了多能微电网 的运行费用。 权利要求书4页 说明书11页 附图3页 CN 114266600 A 2022.04.01 CN 114266600 A 1.一种多能微电网非线性综合需求响应的管理方法， 其特 征在于： 包括如下步骤： 初始化每小时负荷需求， 电价、 天然气价格和预测的风机出力； 确定风机出力， 电力和天然气价格的波动范围； 初始化进化捕食算法中的种群规模以及最大迭代次数； 根据最短距离算法确定最佳峰谷平时段， 获得电力和燃气负荷的最优分时价格， 并计 算电力弹性系数矩阵； 依次执行电与气的幂、 指数、 对数以及线性需求响应模型， 以获得每小时电负荷与气负 荷的需求响应； 根据可信度理论确定风机出力的置信度模型； 评估多能微电网中的约束， 计算目标函数并通过进化捕食算法更新决策变量。 2.根据权利要求1所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 所 述的风机出力PW为： 其中， V、 Vci、 Vra和Vco 分别是实际风速、 切入风速、 额定风速和 切出风速， Pra为风电场的额定有功 功率。 3.根据权利要求1所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 所 述的种群规模为200， 所述的最大迭代次数为10000， 进化捕食算法的策略捕食者比例为 30％。 4.根据权利要求1所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 所 述的峰谷平时段为： t∈G1， if|dt‑m1|＝min{|dt‑mj|j＝1， 2， 3} t∈G2， if|dt‑m2|＝min{|dt‑mj|j＝1， 2， 3}； t∈G3， if|dt‑m3|＝min{|dt‑mj|j＝1， 2， 3} 其中， {d(t)， t＝1， ...， T}为每个时刻配电网系统的初始负荷需求， 令d ′(t)＝max{d (t)‑|| μ||1， 0}； 记m1、 m2、 m3分别为d(t)最小值、 中间值和最大值， G1、 G2、 G3分别为谷时段、 平 时段和峰时段； 对{d ′D(t)}进行排序， 得到d1＜d2＜...＜dT∈{d′D(t)}。 5.根据权利要求4所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 电 力弹性系数矩阵E(t， t ′)为： 其中， t和t ’表示时间， ρ0和ρ 分别为需求响应前后的电价。权　利　要　求　书 1/4 页 2 CN 114266600 A 26.根据权利要求5所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 电 的幂、 指数、 对数以及线性需求响应模型分别为： 其中， dpow(t)、 dexp(t)、 dlog(t)、 dlin(t)分别为电的幂、 指数、 对数以及线性需求响应模 型， d(t)和d0(t)表示t时刻需求响应之后与之前的电负荷需求量， ρ0(t)和ρ(t)分别表示初 始电价和实时电价。 7.根据权利要求5所述的多能微电网非线性综合需求响应的管理方法， 其特征在于： 气 的幂、 指数、 对数以及线性需求响应模型分别为： 其中， gpow(t)、 gexp(t)、 glog(t)、 glin(t)分别为气的幂、 指数、 对数以及线性需求响应模 型， g(t)和g0(t)表示t时刻需求响应之后与之前的气负荷需求量， ρg0(t)和ρg(t)分别表示权　利　要　求　书 2/4 页 3 CN 114266600 A 3