(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110163295.4 (22)申请日 2021.02.0 5 (65)同一申请的已公布的文献号 申请公布号 CN 112949909 A (43)申请公布日 2021.06.11 (73)专利权人 清华大学深圳国际研究生院 地址 518055 广东省深圳市南 山区西丽 街 道深圳大 学城清华校区A栋二楼 (72)发明人 石立宝　翟放　 (74)专利代理 机构 深圳新创友知识产权代理有 限公司 4 4223 代理人 江耀纯 (51)Int.Cl. G06Q 10/04(2012.01) G06Q 50/06(2012.01)G06F 30/27(2020.01) G06N 3/00(2006.01) G06K 9/62(2022.01) H02J 3/46(2006.01) G06F 113/04(2020.01) (56)对比文件 CN 107330546 A,2017.1 1.07 CN 108448650 A,2018.08.24 CN 109785183 A,2019.0 5.21 CN 109118024 A,2019.01.01 CN 106940869 A,2017.07.1 1 陈云翔等.考虑可靠性的时域证据组合方 法. 《控制与决策》 .2018, 审查员 高静 (54)发明名称 一种考虑多种不确定性并存的机组组合优 化方法 (57)摘要 本发明公开了一种考虑多种不确定性并存 的机组组合优化方法， 包括以下步骤： S1、 输入电 网、 火电机组和确定性负荷基本 数据； S2、 输入新 能源机组的不确定性变量； S3、 将不确定性变量 转化为DS结构后再生成噪声符号ε； S4、 基于优 先排序法生成决策变量的初始群体 和 S5、 基于扩展仿射算术理论和直流潮流B系数法 计算电网的网损； S6、 建立最小化电网调度总成 本的目标函数及约束条件， 并计算求解得到适应 度函数值和总 约束违反TCV； S7、 对适应度函数值 和总约束违 反TCV进行适应度比较得到前三名较 优解； S8、 输出前三名较优 解。 本发明有利于电力 系统调度者进行 更安全和经济的优化调度。 权利要求书3页 说明书8页 附图4页 CN 112949909 B 2022.03.08 CN 112949909 B 1.一种考虑多种不确定性并存的机组 组合优化方法， 其特 征在于， 包括以下步骤： S1、 获取电网、 火电机组和 确定性负荷基本数据； S2、 获取新能源机组的不确定性变量； S3、 将所述不确定性变量表示为列变量κ， 并基于证据理论将列变量转化为Dempster – Shafer结构， 然后将转 化为Dempster –Shafer结构的列变量κ 进行线性变化 生成噪声符号 ε； S4、 设定优先排序法的迭代次数， 并基于优先排序法生成决策变量的初始群体 和 其中 为t时刻第i个火电机组的启停变量， 启停变 量启动为1、 关停为0， 为t时 刻第i个火电机组的有功出力； S5、 利用步骤S3得到的列变量κ、 噪声符号ε， 以及步骤S4得到的初始群体 和 并基于扩展仿射 算术理论和直 流潮流B系数法计算电网的网损； S6、 建立最小化电网调度总成本的目标函数及约束条件， 并对所述目标函数及约束条 件多次迭代计算 求解得到适应度函数值和总约束违反TCV； 建立的最小化电网调度总成本的目标函数如公式(6)、 (7)、 (8)所示： 公式(6)、 (7)、 (8)中： t为时间变量， T为调度时间总数， NG为火电机组总数， 为火 电机组启停成本函数， 为火电机组燃料成本函数， 为t时刻第i个火电机组的启 动成本， 为t时刻第i个火电机组的关停成本， ξp为功率偏 移的惩罚成本， 为t时刻 电网的功率偏差， ai、 bi、 ci为第i个火电机组的成本系数， DS表示基于二次多项式表示的输 出量通过Dempster –Shafer结构的卷积运 算得到的Dempster –Shafer结构； 建立的约束条件包括火力发电机组出力上下限约束、 最小开停机时间约束、 爬坡约束、 功率平衡约束、 备用容 量约束、 网络安全约束， 其中： 火力发电机组出力上 下限约束如公式(9)所示： 最小开停机时间约束如公式(10)、 公式(1 1)所示： 爬坡约束如公式(12)、 (13)所示： 功率平衡约束如公式(14)所示：权　利　要　求　书 1/3 页 2 CN 112949909 B 2备用容量约束如公式(15)所示： 网络安全约束如公式(16)所示： 公式(9)‑公式(16)中： NW为新能源机组总数， NB为电网支路总数， LG(i)和UG(i)分别为第i 个火电机组的出力下限和上限， 和 分别为t时刻第i个火电机组已经启动或者 关停的时间， MUTG(i)和MDTG(i)分别为第i个火电机组最小开机和关停时间， URG(i)和DRG(i)分 别为第i个火电机组最大和最小爬坡速率， 为t时刻b支路上负荷的二次多项式， Δ为允 许的功率偏差量， rm为备用系统容量系数， 为t时刻第 i个火电机组最大调 整功率， Cl 为l支路的线路最大容量， 为b支路上、 连接l和k支路的线路转移因子， 为连接在b支 路上的火电机组集合， 为连接在b支路上的风电场集合， DS表示基于二次多项式表示的 输出量通过Dempster –Shafer结构的卷积运算得到的Dempster –Shafer结构， 为 二次多项式表示， 为t时刻第i个风电场的有功出力， 为 的二次多项式表 示， 为t时刻系统网损的二次多 项式表示； 所述总约束违反TCV的计算方式如下： 记FM和 为变量M的概率上界和下界， m为任意实数， Pr为事件的概率， 则满足不等式条 件为： 设X为决策变量的列向量， gc(X)为第c个DS结构下的不等式约束， 则第c个约束满足不等 式的条件为： a、 最小概率没有低于预先设定的概率阈值σc； b、 否则， 计算约束违反CVc， 公式 如下所示： CVc＝max{σc‑inf Pr[gc(x)≤0]， 0}， 总约束违反TCV的公式如下 所示： S7、 对步骤S6迭代计算得到的适应度函数值和总约束违反TCV进行适应度比较， 并得到 前三名较优解；权　利　要　求　书 2/3 页 3 CN 112949909 B 3