CN117335449B - 一种高比例新能源电力系统最低调频能力需求计算方法 - Google Patents

Abstract

本发明公开了一种高比例新能源电力系统最低调频能力需求计算方法，属于电力系统频率控制技术领域。高比例新能源电力系统最低调频能力需求计算方法包括以下步骤：确定电力系统内的同步机实时开机数量及同步机参数；考虑同步机调速器死区及限幅，构建预想扰动下电力系统频率支撑动态特性模型；基于有限差分法将电力系统频率动态模型时域离散化；根据设定的频率稳定需求指标，建立目标函数及约束条件；使用商用求解器求解模型并得到高比例新能源电力系统频率稳定的最小功率需求轨迹。本发明采用上述高比例新能源电力系统最低调频能力需求计算方法，能够实现预想扰动后的电力系统最小调频功率需求计算，提高系统运行稳定性。

Description

一种高比例新能源电力系统最低调频能力需求计算方法

技术领域

本发明涉及电力系统频率控制技术领域，尤其是涉及一种高比例新能源电力系统最低调频能力需求计算方法。

背景技术

随着“碳达峰、碳中和”目标的进一步推进，高比例可再生能源、高比例电力电子设备的“双高”电网格局将进一步突显，电力系统频率稳定问题也将愈发突出。传统同步发电机与电网直接耦合，当电网发生扰动时，同步发电机自发地释放转子惯性动能以补偿系统功率不平衡，并能够通过调速器快速响应系统频率变化，阻止系统频率持续下降。与传统同步发电机相比，新能源电源和其他并网电力电子设备对电网动态事件表现出无惯性或低惯性，系统惯量水平下降，抵抗扰动的能力降低；并且为了获得最大的发电功率，新能源电源一般运行在最大功率跟踪模式，传统控制模式下不响应电网的频率变化，不具备有功支撑能力。

新能源发电在电网中占比越来越高，其参与系统调频将会成为未来电网发展的一种必然要求。随着电力电子技术的发展，虚拟惯量控制、虚拟同步机技术得到了广泛发展与应用。电力系统受到低频扰动后，通过换流器接口并入电网的新能源可以通过虚拟惯量控制、虚拟同步机控制等方式，以快速、灵活的策略短时间内提供大量能量，平缓系统频率的下跌趋势。

因此，当预想极端故障发生后，如何对新能源的转子动能、有功备用与附加储能协调控制，同时考虑储能参与系统调频的作用，实现对系统综合惯量的提升，成了新能源参与系统调频亟待解决的问题。

发明内容

本发明的目的是提供一种高比例新能源电力系统最低调频能力需求计算方法，实现预想扰动后的电力系统最小调频功率需求计算，帮助系统调度人员实时掌握系统频率安全水平并指导调整系统运行方式，提高系统运行稳定性。

为实现上述目的，本发明提供了一种高比例新能源电力系统最低调频能力需求计算方法，包括以下步骤：

S1、确定电力系统内的同步机实时开机数量及同步机参数；

S2、考虑同步机调速器死区及限幅，构建预想扰动下电力系统频率支撑动态特性模型；

S3、基于有限差分法将电力系统频率支撑动态特性模型时域离散化；

S4、根据设定的频率稳定需求指标，建立目标函数及约束条件；

S5、使用商用求解器求解模型并得到高比例新能源电力系统频率稳定的最小功率需求轨迹。

优选的，所述S2具体包括以下步骤：

S21、将电力系统中已投入运行的同步机组聚合为一台等值同步机；

S22、计及调速器死区与限幅非线性环节，建立预想扰动下电力系统频率支撑动态特性模型。

优选的，所述S21具体包括以下步骤：

S211、对电力系统中已投入运行的同步机组的参数进行等值聚合，将同步机等值聚合成单台机组，将其有功功率施加于一个等值转子上：

系统等效惯性时间常数为：

式中：S_B1，……，S_Bn分别表示n台同步机组的额定容量；H₁，……,H_n分别表示n台同步机组的惯性时间常数；S_N表示系统总额定容量；

S212、以调差系数表示同步机组调速器作用，系统等效调差系数按下式计算：

式中：R_i表示第i台机组的调差系数；

式中：S_N表示系统总额定容量；

S213、等值机组再热器时间常数计算方式为：

式中：T_RHi表示第i台机组的再热器时间常数；

S214、等值机组高压缸功率比例计算方式为：

式中：F_HPi为第i台机组的高压缸功率比例，表示高压缸稳态输出功率与汽轮机总输出功率的比值。

优选的，所述S22具体包括以下步骤：

S221、由摇摆方程刻画电网扰动后的频率动态响应过程：

式中：f₀为发电机额定频率，Δf为发电机频率偏差，H为发电机等效惯性时间常数，D为发电机阻尼系数，ΔP_m为调速器功率变化量，ΔP_e为电磁功率变化量，为时域微分算子；

S222、基于SFR模型刻画火电机组的频率响应动态，包括再热式火电机组和非再热式火电机组，将上级电网内各火电机组等值聚合为一台机组，得到理想情况下等值火电机组的调速器-原动机模型：

式中：ΔP_k为等值机组的有功调整量，K为等值机组的调速器增益；T_RH为等值机组的再热器时间常数；F_HP为等值机组的高压缸功率比例；

S223、实际电力系统中，还需要考虑调频机组的调速器死区与限幅非线性环节的约束，等值火电机组的有功功率实际调整量如下式：

式中：Δf_db为频率死区，为等值火电机组的调频有功限幅值，ΔP_k为等值火电机组经调速器后的有功调整量。

优选的，所述S3具体包括以下步骤：

S31、采用前向差分法建立系统频率响应的离散化模型，扰动后的系统功率缺额如下式：

P_dis,n＝P_dis,0-P_i,n-P_k,n (9)

式中，P_dis,n为系统第n个步长的有功功率缺额，P_dis,0为扰动初始有功功率缺额，P_i,n为系统第n个步长时的新能源增发有功功率，P_k,n为系统第n个步长时的等值同步机组增发有功功率；

S32、扰动后的系统频率偏差如下式：

式中，Δf_n为系统第n个步长的频率偏差，Δn为差分步长，H为系统等效惯性时间常数，S_b为系统基准容量。

S33、扰动后的等值同步机组增发有功功率如下式：

式中，P_k,n为系统第n个步长时等值同步机组增发有功功率，b_k,n为系统第n个步长时等值同步机组经调速器死区后的增发有功功率，a_k,n为系统第n个步长时等值同步机组经调速器限幅后的增发有功功率，Δf_db为频率死区。

优选的，所述S4具体包括以下步骤：

S41、确定优化目标并根据扰动发生后的调频时段设立不同权重，建立目标函数；

S42、根据系统频率稳定要求建立约束条件。

优选的，所述S41具体包括以下步骤：

S411、将扰动发生后的调频过程划分为不同阶段，根据时段设立不同的新能源调频权重；

S412、根据扰动发生后的不同划分时段及其权重建立目标函数，具体表示为：

式中，N_a、N_b、N_c表示扰动发生后不同时段的差分起点；Δ_n为差分步长，表示第n差分步长时系统中新能源调频资源增发的有功功率；a、b、c表示扰动发生后不同时段内新能源调频资源分配到的权重。

优选的，所述S42具体包括以下步骤：

S421、基于系统频率稳定最低需求，建立扰动后最大频率变化率约束：

式中，ROCOF^extre为设定的满足系统频率稳定最低需求的最大频率变化率，Δf_n为第n个差分步长时的系统频率偏差。

S422、基于系统频率稳定最低需求，建立扰动后最大频率偏差约束：

|Δf_n|≤Δf^extre (14)

式中，Δf^extre为设定的满足系统频率稳定最低需求的最大频率偏差，Δf_n为第n差分步长时的系统频率偏差。

S423、考虑各机组暂态出力可能性、防止新能源增发有功功率出现频繁的突降及突增，建立频率变化率约束：

式中，不等式左边和右边分别表示第n+1个差分步长和第n个差分步长的系统频率变化率。

优选的，所述S5具体包括以下步骤：

调用商用求解器对目标函数进行优化，计算满足约束条件的各变量最优解并输出各个差分步长的计算结果，最终得到满足系统频率稳定要求的最小有功功率边界曲线。

本发明所述的一种高比例新能源电力系统最低调频能力需求计算方法的优点和积极效果是：

本发明能够实现预想扰动后的电力系统最小调频功率需求计算，帮助系统调度人员实时掌握系统频率安全水平并指导调整系统运行方式，提高系统运行稳定性。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明一种高比例新能源电力系统最低调频能力需求计算方法实施例的流程图；

图2为本发明实施例采用的改进IEEE 39节点系统拓扑示意图；

图3为本发明实施例扰动发生后系统频率偏差计算结果示意图；

图4为本发明实施例扰动发生后节点有功功率缺额计算结果示意图；

图5为本发明实施例扰动发生后等值同步机组增发的有功功率计算结果示意图；

图6为本发明实施例扰动发生后节点新能源机组增发的有功功率计算结果示意图。

具体实施方式

以下通过附图和实施例对本发明的技术方案作进一步说明。

如图1所示，一种高比例新能源电力系统最低调频能力需求计算方法，包括以下步骤：

S1、确定电力系统内的同步机实时开机数量及同步机参数。

步骤S1中需要统计的同步机参数主要包括：额定容量S_n、高压缸功率比例F_HP、再热器时间常数T_RH、惯性时间常数H、一次调频下垂系数K。

S2、考虑同步机调速器死区及限幅，构建预想扰动下电力系统频率支撑动态特性模型。

步骤S2具体包括以下步骤：

S21、将电力系统中已投入运行的同步机组聚合为一台等值同步机。

步骤S21具体包括以下步骤：

S211、系统中同步机有很多台，需要对其参数进行等值聚合，实际系统包含多台发电机组，以再热式机组为例，为简化模型需要将同步机等值聚合成单台机组，将其有功功率施加于一个等值转子上：

系统等效惯性时间常数为：

式中：S_B1，……，S_Bn分别表示n台同步机组的额定容量；H₁，……,H_n分别表示n台同步机组的惯性时间常数；S_N表示系统总额定容量。

S212、以调差系数表示同步机组调速器作用，系统等效调差系数按下式计算：

式中：R_i表示第i台机组的调差系数；

式中：S_N表示系统总额定容量。

S213、等值机组再热器时间常数计算方式为：

式中：T_RHi表示第i台机组的再热器时间常数。

S214、等值机组高压缸功率比例计算方式为：

式中：F_HPi为第i台机组的高压缸功率比例，表示高压缸稳态输出功率与汽轮机总输出功率的比值。

S22、计及调速器死区与限幅非线性环节，建立预想扰动下电力系统频率支撑动态特性模型。

步骤S22具体包括以下步骤：

S221、由摇摆方程刻画电网扰动后的频率动态响应过程：

式中：f₀为发电机额定频率，Δf为发电机频率偏差，H为发电机等效惯性时间常数，D为发电机阻尼系数，ΔP_m为调速器功率变化量，ΔP_e为电磁功率变化量，为时域微分算子。

S222、基于SFR模型刻画火电机组的频率响应动态，包括再热式火电机组和非再热式火电机组。我国火电机组以再热式火电机组为主，将上级电网内各火电机组等值聚合为一台机组，得到理想情况下等值火电机组的调速器-原动机模型：

式中：ΔP_k为等值机组的有功调整量，K为等值机组的调速器增益；T_RH为等值机组的再热器时间常数；F_HP为等值机组的高压缸功率比例。

S223、实际电力系统中，还需要考虑调频机组的调速器死区与限幅等非线性环节的约束，等值火电机组的有功功率实际调整量如下式：

式中：Δf_db为频率死区，为等值火电机组的调频有功限幅值，ΔP_k为等值火电机组经调速器后的有功调整量。

S3、基于有限差分法将电力系统频率支撑动态特性模型时域离散化。

S3具体包括以下步骤：

S31、采用前向差分法建立系统频率响应的离散化模型，以低频扰动为例，得到扰动后的系统功率缺额如下式：

P_dis,n＝P_dis,0-P_i,n-P_k,n (9)

式中，P_dis,n为系统第n个步长的有功功率缺额，P_dis,0为扰动初始有功功率缺额，P_i,n为系统第n个步长时的新能源增发有功功率，P_k,n为系统第n个步长时的等值同步机组增发有功功率。

S32、采用前向差分法建立系统频率响应的离散化模型，以低频扰动为例，得到扰动后的系统频率偏差如下式：

式中，Δf_n为系统第n个步长的频率偏差，Δn为差分步长，H为系统等效惯性时间常数，S_b为系统基准容量。

S33、采用前向差分法建立系统频率响应的离散化模型，以低频扰动为例，得到扰动后的等值同步机组增发有功功率如下式：

式中，P_k,n为系统第n个步长时等值同步机组增发有功功率，b_k,n为系统第n个步长时等值同步机组经调速器死区后的增发有功功率，a_k,n为系统第n个步长时等值同步机组经调速器限幅后的增发有功功率，Δf_db为频率死区。

S4、根据设定的频率稳定需求指标，建立目标函数及约束条件。

S4具体包括以下步骤：

S41、确定优化目标并根据扰动发生后的调频时段设立不同权重，建立目标函数。

S41具体包括以下步骤：

S411、将扰动发生后的调频过程划分为不同阶段，根据时段设立不同的新能源调频权重。例如，扰动发生后0-0.5s内，储能系统可以迅速释放大量能量以改善频率的大幅度跌落，调频价值较高，为鼓励快速调频资源参与辅助服务，应设置最高权重。扰动发生后0.5-5s内风光储新能源系统可以提供快速灵活的虚拟惯量支撑，可设置次高权重。扰动发生后5s之后系统中同步机组一次调频功率得到逐渐提高并作为调频主导，因此新能源可设置较低权重。具体设置情况根据实际调频需求决定。

S412、根据扰动发生后的不同划分时段及其权重建立目标函数，具体表示为：

式中，N_a、N_b、N_c表示扰动发生后不同时段的差分起点，Δ_n为差分步长，表示第n差分步长时系统中新能源调频资源增发的有功功率，a、b、c表示扰动发生后不同时段内新能源调频资源分配到的权重。

该目标函数的物理意义为：使扰动发生后的加权新能源调频资源需求最小化。

S42、根据系统频率稳定要求建立约束条件。

S42具体包括以下步骤：

S421、基于系统频率稳定最低需求，建立扰动后最大频率变化率约束：

式中，ROCOF^extre为设定的满足系统频率稳定最低需求的最大频率变化率，Δf_n为第n个差分步长时的系统频率偏差。

目前，我国电网虽然尚未对频率变化率做出明确要求，但在惯量相对减小的新能源高占比电力系统，扰动下较大的频率变化率可能导致机组内部结构损坏以及分布式电源脱网，威胁机组运行安全，未来考虑会将其作为频率约束指标。此外，国外许多电网规定了系统的频率变化率保护限值。

S422、基于系统频率稳定最低需求，建立扰动后最大频率偏差约束：

|Δf_n|≤Δf^extre (14)

式中，Δf^extre为设定的满足系统频率稳定最低需求的最大频率偏差，Δf_n为第n差分步长时的系统频率偏差。

根据GB 38755-2019《电力系统安全稳定导则》规定，频率稳定是指系统在受到严重扰动后，发电与负荷需求出现大的不平衡情况下，系统频率能够保持或恢复到允许的范围内、不发生频率奔溃的能力。作为频率响应极值，最大频率偏差是调速器参数优化，低频减载、低频解列、高频切机(对应于频率动态上升场景)、超速保护控制、直流调制等安全稳定措施动作值整定的重要依据。

S423、考虑各机组暂态出力可能性、防止新能源增发有功功率出现频繁的突降及突增，建立频率变化率约束：

式中，不等式左边和右边分别表示第n+1个差分步长和第n个差分步长的系统频率变化率。该约束的物理意义是使整个调频过程中的系统频率变化率保持单调不增状态，保证计算结果的调频有效性。

S5、使用商用求解器求解模型并得到高比例新能源电力系统频率稳定的最小功率需求轨迹。

调用商用求解器对目标函数进行优化，计算满足约束条件的各变量最优解并输出各个差分步长的计算结果，最终得到满足系统频率稳定要求的最小有功功率边界曲线。

实施例

如图2所示，基于该改进IEEE 39节点系统进行系统频率稳定的有功功率需求计算，操作步骤如下：

1、在改进IEEE 39节点系统母线25位置处设置负荷突增50MW扰动，系统频率越过同步机调速器一次调频动作死区。统计系统中各台同步发电机的额定容量S_n、高压缸功率比例F_HP、再热器时间常数T_RH、惯性时间常数H、一次调频下垂系数K。

2、将改进IEEE 39节点系统中的同步机组聚合为一台等值同步机，计及调速器死区与限幅等非线性环节，建立预想扰动下电力系统频率支撑动态特性模型。在本实施例中，同步机组调速器调频死区均设置为0.033Hz，有功限幅值设置为容量的10％。

3、采用前向差分法建立系统频率响应的离散化模型，得到扰动后的系统功率缺额、系统频率偏差、等值同步机组增发有功功率离散表达式。在本实施例中，差分步长为0.01s。

在调用商用求解器进行优化计算时，需要对各变量进行初始化。在本实施例中，初始系统频率偏差为0，等值同步机组初始增发有功功率为0。

4、将扰动发生后的调频过程划分为不同阶段，以根据时段设立不同的新能源调频权重，并根据扰动发生后的不同划分时段及其权重建立目标函数。在本实施例中，扰动发生后0-0.5s为高权重，扰动发生后0.5-5s为次高权重，扰动发生5s后为低权重。

5、根据系统频率稳定要求建立约束条件，包括最大系统频率变化率约束、最大系统频率偏差约束、频率变化率单调不增约束等。

根据电力系统安全稳定导则GB38755-2019，频率稳定的判据是系统频率能迅速恢复到额定频率附近继续运行，不发生频率持续震荡或频率崩溃，也不使系统频率长期悬浮于某一过高或过低的数值。根据电能质量电力系统频率偏差GB15945-2008，正常运行条件下的频率偏差为±0.2～0.5Hz。因此，在本实施例中，对于最大系统频率偏差约束设置为0.5Hz。

我国电网虽然尚未对频率变化率做出明确要求，但在惯量相对减小的新能源高占比电力系统，扰动下较大的频率变化率可能导致机组内部结构损坏以及分布式电源脱网，威胁机组运行安全，未来考虑会将其作为频率约束指标。此外，许多电网规定了系统的频率变化率保护限值，如表1所示。参考国外电网频率变化率保护限值，在本实施例中，对于最大频率变化率约束设置为1.0Hz/s。

表1不同电网的频率变化率保护限值

电网名称	频率变化率保护限值/(Hz/s)
		Texas	1.0
Great Britain	0.5
		Ireland	1.0
Nordic system	0.5
		South Australia	1.5～3.0

6、调用商用求解器对目标函数进行优化，计算满足约束条件的各变量最优解并输出各个差分步长的计算结果。最终得到扰动发生后的系统频率偏差计算结果曲线，如图3所示；同样地，可以得到扰动发生后节点有功功率缺额计算结果曲线、扰动发生后等值同步机组增发的有功功率计算结果曲线、满足系统频率稳定要求的最小有功功率边界曲线，分别如图4、图5、图6所示。

因此，本发明采用上述高比例新能源电力系统最低调频能力需求计算方法，实现预想扰动后的电力系统最小调频功率需求计算，帮助系统调度人员实时掌握系统频率安全水平并指导调整系统运行方式，提高系统运行稳定性。

最后应说明的是：以上实施例仅用以说明本发明的技术方案而非对其进行限制，尽管参照较佳实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对本发明的技术方案进行修改或者等同替换，而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。

Claims (7)

1.一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，包括以下步骤：

S1、确定电力系统内的同步机实时开机数量及同步机参数；

S2、考虑同步机调速器死区及限幅，构建预想扰动下电力系统频率支撑动态特性模型；

S3、基于有限差分法将电力系统频率支撑动态特性模型时域离散化；

S4、根据设定的频率稳定需求指标，建立目标函数及约束条件；

S5、使用商用求解器求解模型并得到高比例新能源电力系统频率稳定的最小功率需求轨迹；

所述S2具体包括以下步骤：

S21、将电力系统中已投入运行的同步机组聚合为一台等值同步机；

S22、计及调速器死区与限幅非线性环节，建立预想扰动下电力系统频率支撑动态特性模型；

所述S22具体包括以下步骤：

S221、由摇摆方程刻画电网扰动后的频率动态响应过程：

式中：f₀为发电机额定频率，Δf为发电机频率偏差，H为发电机等效惯性时间常数，D为发电机阻尼系数，ΔP_m为调速器功率变化量，ΔP_e为电磁功率变化量，为时域微分算子；

S222、基于SFR模型刻画火电机组的频率响应动态，包括再热式火电机组和非再热式火电机组，将上级电网内各火电机组等值聚合为一台机组，得到理想情况下等值火电机组的调速器-原动机模型：

式中：ΔP_k为等值机组的有功调整量，K为等值机组的调速器增益，T_RH为等值机组的再热器时间常数，F_HP为等值机组的高压缸功率比例；

S223、实际电力系统中，还需要考虑调频机组的调速器死区与限幅非线性环节的约束，等值火电机组的有功功率实际调整量如下式：

式中：Δf_db为频率死区，为等值火电机组的调频有功限幅值，ΔP_k为等值火电机组经调速器后的有功调整量。

2.根据权利要求1所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S21具体包括以下步骤：

S211、对电力系统中已投入运行的同步机组的参数进行等值聚合，将同步机等值聚合成单台机组，将其有功功率施加于一个等值转子上：

系统等效惯性时间常数为：

式中：S_B1，……，S_Bn分别表示n台同步机组的额定容量；H₁，……,H_n分别表示n台同步机组的惯性时间常数；S_N表示系统总额定容量；

S212、以调差系数表示同步机组调速器作用，系统等效调差系数按下式计算：

式中：R_i表示第i台机组的调差系数；

式中：S_N表示系统总额定容量；

S213、等值机组再热器时间常数计算方式为：

式中：T_RHi表示第i台机组的再热器时间常数；

S214、等值机组高压缸功率比例计算方式为：

式中：F_HPi为第i台机组的高压缸功率比例，表示高压缸稳态输出功率与汽轮机总输出功率的比值。

3.根据权利要求2所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S3具体包括以下步骤：

S31、采用前向差分法建立系统频率响应的离散化模型，扰动后的系统功率缺额如下式：

P_dis,n＝P_dis,0-P_i,n-P_k,n (9)

式中，P_dis,n为系统第n个步长的有功功率缺额，P_dis,0为扰动初始有功功率缺额，P_i,n为系统第n个步长时的新能源增发有功功率，P_k,n为系统第n个步长时的等值同步机组增发有功功率；

S32、扰动后的系统频率偏差如下式：

式中，Δf_n为系统第n个步长的频率偏差，Δn为差分步长，H为系统等效惯性时间常数，S_b为系统基准容量；

S33、扰动后的等值同步机组增发有功功率如下式：

式中，P_k,n为系统第n个步长时等值同步机组增发有功功率，b_k,n为系统第n个步长时等值同步机组经调速器死区后的增发有功功率，a_k,n为系统第n个步长时等值同步机组经调速器限幅后的增发有功功率，Δf_db为频率死区。

4.根据权利要求3所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S4具体包括以下步骤：

S41、确定优化目标并根据扰动发生后的调频时段设立不同权重，建立目标函数；

S42、根据系统频率稳定要求建立约束条件。

5.根据权利要求4所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S41具体包括以下步骤：

S411、将扰动发生后的调频过程划分为不同阶段，根据时段设立不同的新能源调频权重；

S412、根据扰动发生后的不同划分时段及其权重建立目标函数，具体表示为：

式中，N_a、N_b、N_c表示扰动发生后不同时段的差分起点，Δ_n为差分步长，ΔP_i,n表示第n差分步长时系统中新能源调频资源增发的有功功率，a、b、c表示扰动发生后不同时段内新能源调频资源分配到的权重。

6.根据权利要求5所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S42具体包括以下步骤：

S421、基于系统频率稳定最低需求，建立扰动后最大频率变化率约束：

式中，ROCOF^extre为设定的满足系统频率稳定最低需求的最大频率变化率，Δf_n为第n个差分步长时的系统频率偏差；

S422、基于系统频率稳定最低需求，建立扰动后最大频率偏差约束：

|Δf_n|≤Δf^extre (14)

式中，Δf^extre为设定的满足系统频率稳定最低需求的最大频率偏差，Δf_n为第n差分步长时的系统频率偏差；

S423、考虑各机组暂态出力可能性、防止新能源增发有功功率出现频繁的突降及突增，建立频率变化率约束：

式中，不等式左边和右边分别表示第n+1个差分步长和第n个差分步长的系统频率变化率。

7.根据权利要求6所述的一种高比例新能源电力系统最低调频能力需求计算方法，其特征在于，所述S5具体包括以下步骤：

调用商用求解器对目标函数进行优化，计算满足约束条件的各变量最优解并输出各个差分步长的计算结果，最终得到满足系统频率稳定要求的最小有功功率边界曲线。