CN109149566A - 一种大功率缺失下频率最低点预测的仿真模型的建模方法 - Google Patents

Abstract

一种大功率缺失下频率最低点预测的仿真模型的建模方法属于电力***分析领域，包括以下步骤：S1：对各机组调速器阶跃响应曲线进行离散微分运算得到其冲激响应曲线；S2：利用多项式对冲激响应曲线进行拟合；S3：利用最小二乘法确定多项式拟合系数；S4：对步骤S2所得多项式进行Laplace变换，求解调速器的近似传递函数；S5：计算***等效调速器的传递函数；S6：基于步骤S5所得的***等效调速器传递函数，建立大功率缺失下频率最低点预测的仿真模型。本发明基于各调速器阶跃响应曲线，利用所提技术方案，可建立大功率缺失下频率最低点预测的仿真模型，从而准确快速计算大功率缺失下的最大频率偏移和频率最低点时间。

Description

一种大功率缺失下频率最低点预测的仿真模型的建模方法

技术领域

本发明属于电力***分析领域，涉及一种大功率缺失下频率最低点预测的仿真模型的建模方法。

背景技术

随着高比例可再生能源的接入，风力、光伏等新能源机组将逐步成为电力***的主力电源。相比传统发电机组，新能源机组的惯性响应和频率响应能力先天不足，其大规模并网将使得***频率稳定问题日益凸显。另一方面，随着特高压交直流输电工程的发展，由于直流闭锁故障造成的大功率缺失，易引发直流馈入受端电网的低频事故。上述两方面因素的作用，将导致频率控制“攻防”两端出现此消彼长的态势，从而使电力***的频率稳定面临严峻的挑战。

针对大功率缺失引发的电网频率骤降，快速精准的频率响应资源调用和低频减载控制是维持电力***频率稳定的重要手段。最大频率偏移和频率最低点时间作为频率紧急控制的决策参量，对其进行准确地预测和估计尤为关键。

发明内容

根据现有技术存在的问题，本发明公开了一种大功率缺失下频率最低点预测的仿真模型的建模方法。

本发明的技术方案为：

一种大功率缺失下频率最低点预测的仿真模型的建模方法，包括以下步骤：

S1：在大功率缺失下的频率响应过程中，各机组依据本地频率偏差进行功率调整。为模拟各机组的频率响应动态过程，需建立各机组调速器的传递函数，通常可由其冲激响应进行Laplace变换获得。为此，收集各机组的调门扰动试验曲线，亦即调速器的阶跃响应曲线。对各机组调速器阶跃响应曲线进行离散微分运算得到其冲激响应曲线其中，t为响应时间，i为机组编号。

S2：利用n次多项式对冲激响应曲线在t∈[0,t_fit]内进行拟合，从而获取各机组调速器冲激响应的时域解析表达式。

所述的n次多项式如公式(1)所示：

式中：为多项式拟合系数向量；t为响应时间；i为机组编号；t_fit为时域拟合时长；n为拟合阶数。

S3：为求解各机组调速器冲激响应曲线的多项式拟合参数，定义拟合误差并计算冲激响应曲线中各采样点的均方误差。通过最小二乘法最小化累积拟合误差，求解多项式拟合系数

所述的拟合误差如公式(2)所示：

式中：m为时域拟合区间内的样本数量；k为样本数据编号；i为机组编号；t_k为采样时间；t_fit为时域拟合时长。

S4：基于各机组调速器冲激响应的时域解析表达式，利用Laplace变换求解其复频域内的传递函数。对步骤S2所得n次多项式进行Laplace变换，求解调速器i的近似传递函数为：

式中：t为响应时间；s为复频域算子；n为拟合阶数；l为拟合阶数编号；i为机组编号；为机组i的n次多项式拟合系数；为Laplace变换。

S5：***频率响应是所辖区域内各机组频率响应的综合。在发生大功率缺失故障后，所有机组共同参与频率响应过程，其复频域内等效的调速器模型是前述各机组调速器模型的叠加。因此，包含N_G台机组的***等效调速器传递函数为：

其中：

Γ＝diag{0！,1！,…,n！} (6)

式中：K_i为标幺值转换系数，即从机组容量标幺值转化为***容量标幺值；N_G为***内频率响应机组的数量；i为机组编号；s为复频域算子；n为拟合阶数；l为拟合阶数编号；为机组调速器i冲激响应的多项式拟合系数向量；为***等效调速器冲激响应的多项式拟合系数向量。

S6：忽略***网络结构的影响，保留各机组调速器动态，基于***等效惯性和步骤S5所得的***等效调速器传递函数建立大功率缺失下频率最低点预测的仿真模型。

与现有技术相比，本发明的有益效果为：

(1)数据需求少。大功率缺失下频率最低点预测的仿真模型无须已知各调速器结构及参数，仅对调门扰动试验曲线进行多项式拟合便可建立***的频率响应模型。

(2)模型阶数低。大功率缺失下频率最低点预测的仿真模型阶数不随机组数量增多而进一步升高，有效提高了模型的求解计算速度。

(3)计算精度高。多项式拟合阶数和拟合时长对模型求解结果影响不大，采用低阶模型便可获取较高的计算精度。

附图说明

图1为大功率缺失下频率最低点预测仿真建模的流程图；

图2为实施例频率最低点预测的仿真模型。

具体实施方式

为使本发明的技术方案和优点更加清楚，下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

如图1所示的一种大功率缺失下频率最低点预测的仿真建模流程，包括以下步骤：

S1：在大功率缺失下的频率响应过程中，各机组依据本地频率偏差进行功率调整。为模拟各机组的频率响应动态过程，需建立各机组调速器的传递函数，通常可由其冲激响应进行Laplace变换获得。为此，收集各机组的调门扰动试验曲线，亦即调速器的阶跃响应曲线。在本实施例中，共有10台发电机组，其中G1～G9为火电机组，采用IEEEG1型调速器，G10为水电机组，采用IEEEG3型调速器。对各机组调速器阶跃响应曲线进行离散微分运算得到其冲激响应曲线其中，t为响应时间，i为机组编号。

S2：利用2次多项式对10台发电机组调速器的冲激响应曲线在0～5s内进行拟合，从而获取各机组调速器冲激响应的时域解析表达式。

所述的2次多项式如公式(1)所示：

式中：为多项式拟合系数向量；t为响应时间；i为机组编号；时域拟合时长t_fit＝5。

S3：为求解各机组调速器冲激响应曲线的多项式拟合参数，定义拟合误差并计算冲激响应曲线中各采样点的均方误差。通过最小二乘法最小化累积拟合误差，求解多项式拟合系数

所述的拟合误差如公式(2)所示：

式中：m为时域拟合区间内的样本数量；k为样本数据编号；i为机组编号；t_k为采样时间；时域拟合时长t_fit＝5。

利用最小二乘法计算所得G1～G10调速器5s的二阶多项式拟合结果如表1所示。

表1各机组二阶多项式拟合结果

S4：基于各机组调速器冲激响应的时域解析表达式，利用Laplace变换求解其复频域内的传递函数。对步骤S2所得2次多项式进行Laplace变换，求解调速器i的近似传递函数为：

式中：s为复频域算子；为机组i的多项式拟合系数；i为机组编号。

S5：***频率响应是所辖区域内各机组频率响应的综合。在G6发生机组跳闸故障后，其余机组共同参与频率响应过程，其复频域内等效的调速器模型是前述各机组调速器模型的叠加。因此，包含10台发电机组的***等效调速器传递函数为：

其中，标幺值转换系数K＝[0.5263 0.0526 0.0526 0.0526 0.0526 0.05260.0526 0.0526 0.0526 0.0526]。

S6：忽略***网络结构的影响，保留各机组调速器动态，基于***等效惯性和步骤S5所得的***等效调速器传递函数建立大功率缺失下频率最低点预测的仿真模型如图2所示。

Claims (1)

1.一种大功率缺失下频率最低点预测的仿真模型的建模方法，其特征在于以下步骤：

S1：建立各机组调速器的传递函数，模拟各机组的频率响应动态过程；收集各机组的调门扰动试验曲线，即调速器的阶跃响应曲线，对各机组调速器阶跃响应曲线进行离散微分运算得到其冲激响应曲线其中，t为响应时间，i为机组编号；

S2：利用n次多项式对冲激响应曲线在t∈[0,t_fit]内进行拟合，获取各机组调速器冲激响应的时域解析表达式；

所述的n次多项式如公式(1)所示：

式中：为多项式拟合系数向量，t_fit为时域拟合时长，n为拟合阶数；

S3：定义拟合误差求解各机组调速器冲激响应曲线的多项式拟合参数，并计算冲激响应曲线中各采样点的均方误差；通过最小二乘法最小化累积拟合误差，求解多项式拟合系数

所述的拟合误差如公式(2)所示：

式中：m为时域拟合区间内的样本数量，k为样本数据编号，t_k为采样时间；

S4：对步骤S2所得n次多项式进行Laplace变换，求解调速器i的近似传递函数为：

式中：s为复频域算子，l为拟合阶数编号，p为多项式拟合系数；

S5：大功率缺失故障发生后，所有机组共同参与频率响应过程，其复频域内等效的调速器模型是前述各机组调速器模型的叠加，包含N_G台机组的***等效调速器传递函数为：

其中：

Γ＝diag{0！,1！,…,n！} (6)

式中：K_i为标幺值转换系数，即从机组容量标幺值转化为***容量标幺值；N_G为***内频率响应机组的数量，为机组调速器i冲激响应的多项式拟合系数向量；为***等效调速器冲激响应的多项式拟合系数向量；

S6：忽略***网络结构的影响，保留各机组调速器动态，基于***等效惯性和步骤S5所得的***等效调速器传递函数建立大功率缺失下频率最低点预测的仿真模型。