系统频率响应模型(SFR)推导与安全指标

黄昏ivi

已于 2025-04-17 15:13:28 修改

阅读量549

点赞数 9

分类专栏：电气基本概念文章标签：算法机器学习人工智能

于 2025-04-17 15:12:59 首次发布

本文链接：https://blog.csdn.net/m0_69522810/article/details/147306430

版权

51 篇文章

订阅专栏

电力系统频率响应分析的核心是建立扰动下的动态方程，结合同步机组和逆变器的动态特性：

$\sum\left(H_{gi}+H_{vj}\right)\frac{d\Delta f(t)}{dt}=-D\Delta f(t)+\Delta P_L-\sum\Delta P_{gi}(t)-\sum\Delta P_{vj}(t)$

参数说明：
$H_{gi}$ 、 $H_{vj}$ 为机组惯性时间常数； $D$ 为负载阻尼系数； $\Delta P_{L}$ 为有功扰动； $\Delta P_{gi}(t)$ 、 $\Delta P_{vj}(t)$ 为机组调频功率调整量。

推导关键步骤：

$\Delta f(s)=\frac{\Delta P_L-\sum\Delta P_{gi}(s)-\sum\Delta P_{vj}(s)}{s\sum(H_{gi}+H_{vj})+D}$

2、调频环节建模：
- 同步机组调频功率： $\Delta P_{gi}(s)=\frac{R_{gi}}{1+T_{ri}s}\Delta f(s)$
- 虚拟同步机调频功率： $\Delta P_{vj}(s)=\frac{R_{vj}}{1+T_{vj}s}\Delta f(s)$
  其中 $R_{gi},R_{vj}$ 为调速器增益， $T_{ri},T_{vj}$ 为时间常数。
3、参数聚合：将多机组等效为系统级参数：

$H_t=\sum\left(\frac{S_i}{S_{sys}}u_{i,t}H_i\right),\quad K_{sys}=\sum\left(\frac{S_i}{S_{sys}}u_{i,t}R_i\right)$

其中 $S_{i}$ 为机组容量， $u_{i,t}$ 为运行状态变量。

通过传递函数极点分析，得到频率偏差的时域表达式：

$\Delta f(t)=\frac{\Delta P}{D+K_{sys}}\left(1-e^{-\frac{D+K_{sys}}{2H_t}t}\cos\omega_dt\right)$

关键参数：
- 阻尼比： $\zeta=\frac{D+K_{sys}}{2\sqrt{H_t(D+K_{sys})}}$
- 阻尼振荡频率： $\omega_n=\sqrt{\frac{D+K_{sys}}{H_t}}$
- 自然振荡频率： $\omega_d=\omega_n\sqrt{1-\zeta^2}$

频率最低点（Nadir）计算：
通过求导确定极值时间 $t_{nadir}$ ，代入时域方程得：

$\Delta f_{nadir}=\frac{\Delta P}{D+K_{sys}}\left(1-\frac{\zeta}{\sqrt{\zeta^2+1}}\right)$

此过程需结合系统传递函数的动态特性分析。

$\Delta f_{ss}=\frac{\Delta P}{D+K_{sys}}$

$RoCoF_{max}=\left|\frac{d\Delta f}{dt}\right|_{max}=\frac{\Delta P}{2H_t}$

$FSM=\frac{\Delta f_{nadir}-\Delta f_{threshold}}{\Delta f_{threshold}}$

指标名称	表达式	物理意义
系统等效惯性时间常数	$H_t=\sum\left(\frac{S_i}{S_{sys}}u_{i,t}H_i\right)$	反映系统整体惯性
频率最低点	$\begin{aligned} & \Delta f_{nadir}= \\ & \frac{\Delta P}{D+K_{sys}}\left(1-\frac{\zeta}{\sqrt{\zeta^2+1}}\right) \end{aligned}$	最大瞬时频率跌落
频率变化率（RoCoF）	$RoCoF_{max}=\frac{\Delta P}{2H_t}$	频率跌落速度，与惯性成反比
频率安全边际指数	$FSMI=\frac{\Delta f_{nadir}-\Delta f_{threshold}}{\Delta f_{threshold}}$	量化频率安全裕度