一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优化方法
【技术领域】
[0001] 本发明涉及电力系统的运行、分析与调度领域，特别涉及一种考虑动态频率约束 的含高渗透率光伏电源的孤立电网机组组合优化方法。
【背景技术】
[0002] 由直流换流站联结，形成了直流孤立电网送电的网架结构成为一种趋势。互联系 统中各区的惯性响应、一次调频互为支持，系统频率调整能力较强。而孤立电网内机组数 少，惯性系数小，一次调频响应能力有限。由于光伏电站无旋转部件以及光伏并网逆变器的 隔离，光伏电站无惯性以及一、二次调频能力，高渗透率光伏电源接入孤立电网后，无论其 替代部分常规电源，还是作为新增电源，均会进一步削弱系统动态调频能力。含高渗透率光 伏电源的孤立电网遭受紧急事故，频率的过低或者过高容易导致低频减载或者高频切机动 作，严重时将导致全网频率崩溃。且机组组合优化结果直接决定系统的惯性系数以及一次 调频响应能力。相同的故障扰动下，不同的机组组合优化结果对应的系统动态频率最低点 不同。
[0003] 而目前广泛采用的机组组合优化模型采用的约束条件为系统的有功平衡约束、机 组最大与最小出力约束、机组最小运行与停运持续时间约束，部分文献考虑了稳态情况下 的网络安全约束，即潮流约束，而机组组合对动态频率特性的影响，以及考虑动态频率约束 的机组组合已有的文献尚未涉及。
[0004] 此外，高渗透率光伏电源接入孤立电网面临动态调频能力不足的风险。已有不少 文献针对光伏电源参与系统调频展开讨论，得到光伏电源在频率动态过程中能快速响应改 变出力，提高系统安全性的结论。而以上研究均是针对光伏电源调频的控制策略或是通过 仿真对光伏电源调频进行探讨，并未针对光伏电源通过机组组合所需预留的调频容量进行 探讨。然而，光伏电源通过机组组合所预留的调频容量直接影响系统的频率稳定性以及系 统的发电成本。
【发明内容】
[0005] 本发明针对现有机组组合技术中缺少考虑动态频率约束的问题，从而提供一种考 虑动态频率约束的机组组合方法。且现有的光伏电源减载参与调频的技术未有确定光伏调 频容量的方法，本发明提供一个在机组组合中预留合理的光伏调频容量的方法。从而提供 一种兼顾安全性和经济性的考虑动态频率约束的孤立电网机组组合方法。
[0006] 本发明提供了一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组 合优化方法，其特征在于，基于一个光伏电源并网下考虑动态频率约束机组组合的优化规 划模型，该模型基于以下目标函数以及约束条件：
[0007] 所述目标函数定义为系统发电成本最低，基于：
[0009] 式中N为常规机组的台数，T为研究周期内小时数，Plh为常规机组i在时段h的 输出有功功率，f (Plh)为常规机组i的运行成本，Ulh为常规机组i在时段h的启停机状态， Uih= 1表示运行，U ih= 0表示停机，用式表示为
为成本函 数的系数；Slh为机组i在时段h的启动成本；
[0010] 所述约束条件为：
[0011] 常规机组发电、光伏发电与负荷功率平衡约束：
式二；系统 旋转备用容量
式四；机组有功出力上 下限约束P1minS P lh< P Γ5式五；机组有功功率爬坡速率约束D P ^P1 L i式六；最 小启停时间约束DT, ST;f Λ L、UT； 式八；动态频率最低点约束f_ 49. OHz式九； 光伏调频备用容量约束
分别为第j个光伏电源、 第h时段的预测值和调度值，Rh为第h时段负荷预测值大小；P Γ\ ΡΓη分别为机组i的 出力上、下限；为第j个光伏电源预测误差百分比为第i台机组的爬坡上、下限；
分别为第i台机组在h时段连续开、停时间；UTp DT1*别为第i台机组最小开、 停时间；f_ h为第h时段，系统发电功率损失后动态频率的最小值；
为第j个光伏电 源、在h时段的最小调频容量；
[0012] 优化方法具体包括以下步骤：
[0013] 步骤1 :获取电力系统的每个发电机组的机组特性数据ΡΓΧ、ΡΓη、Lp Dp UTp DT1; 负荷预测数据Plh;光伏电源出力预测数据
[0014] 步骤2 :将优化规划模型分解成3个相互制约的子问题，分别为无动态频率约束以 及光伏调频容量约束的机组组合问题SPl、动态频率校验问题SP2、光伏调频备用容量校验 问题SP3 ;
[0015] 步骤3 :计算无动态频率约束以及光伏调频容量约束的机组组合问题SP1，问题 SPl的目标函数为式一，约束条件为式二至式八，求解子问题SPl获得当前机组启停方案；
[0016] 步骤4 :子问题SP2为动态频率最低点即约束式九的校验；将步骤3所得机组启停 方案、以及出力代入检测问题SP2 ;如果SP2问题无越限情况，则进入下一步；如果检测到某 个时段有动态频率最低点越限，且该时段系统发电机单位调节功率有增加的可能，则产生 优化割，并在SPl中增加一个新的约束条件，然后计算SP1，重复步骤4和步骤5使SPl和 SP2反复相互迭代，直至SP2问题满足；
[0017] 步骤5 :将步骤3、步骤4迭代收敛最终所得机组启停方案、机组出力作为光伏调频 容量检测问题SP3运算的初始值，计算SP3子问题的目标函数，若SP3子问题目标函数为0， 则代表光伏调频容量校验子问题SP3无越限，所得机组启停方案、机组出力为系统最优方 案；若越限，则进入下一步；
[0018] 步骤6 :将步骤5中得到的非0惩罚变量值进行标记，并在SPl中增加一个新的约 束条件，然后计算SP1，重复步骤4、5、6 ;使SP1、SP2、SP3反复相互迭代，直至得到系统最优 方案。
[0019] 在上述的一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优 化方法，所述步骤4中计算动态频率最低点的方法如下：
[0024] 其中M为光伏电站个数，&为发电机i的单位调节功率，T $发电机i的时间常 数，Kv]为光伏电站j的单位调节功率，T v]为光伏j的时间常数，H屬h时段系统的惯性时 间常数；
[0025] 将步骤3中计算得到的各机组启停方案作为步骤4中已知的初始条件，代入式 十二中，并基于式十一至式十四，求得各时段在故障下，最大频降f_ h;若某些时段f _ 满足条件，且该时段仍然有机组处于关机状态，则标记这些时段为h%该时段产生相应的优 化割元>〇 ;其中
为系统的单位调节功率，鳥代表上一轮迭代 计算问题SPl所产生的系统单位调节功率大小。
[0026] 在上述的一种考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合优 化方法，步骤5中对应的子问题SP3为：
[0030] 其中为第j个光伏电站在时段h的预测值、通过步骤4、5迭代收敛后得 到的第j个光伏电站在时段h的调度值、Δ^Γ1为通过故障后的最大频降A f_ M及最大 频降出现的时间t_ ^十算得到的光伏电源最小调频容量、s ]h为当光伏电站j在h时段通 过机组组合所预留的调频容量小于其最小调频容量时引入的松弛变量；其中第j个光伏电 站在时段h的最小调频容量计算式为：
[0032] 在上述的考虑动态频率约束的含高渗透率光伏电源的孤立电网机组组合方法，所 述步骤6中，返回的Benders割为：
其中jh为拉格朗日乘子。
[0033] 本发明与现有技术相比，具有以下优点和效果：本发明在机组组合模型中考虑动 态频率约束，并且让光伏电源减载参与调频来增强系统的调频能力，并确定机组组合在优 化过程留有充足但不过量的光伏电源调频容量，克服了高渗透率光伏接入系统使得系统调 频能力下降的问题，以及传统机组组合模型优化结果因不满足动态频率约束而无法实施的 问题。且本方明所提的求解方法能针对所提模型进行有效的求解，本方明具有良好的推广 应用价值和前景。
【附图说明】
[0034] 图1为系统故障后动态频率响应图。
[0035] 图2为光伏频率响应系统的等效模型。
[0036] 图3为系统的动态频率响应模型。
[0037] 图4本发明方法的计算流程图。
【具体实施方式】
[0038] 下面通过实施例，并结合数据分析，对本发明的技术方案作进一步具体的说明。
[0039] 实施例：
[0040] 本专利所提方法在多个算例模型下进行了验证，限于篇幅，本实施例针对以改进 的十机系统算例为例，对本文所提方法的可行性及有效性进行分析及验证。具体情况如 下：
[0041] 在系统中接入2个光伏电站，1号光伏电源额定容量为500丽，2号光伏电源额定容 量为300MW，两