简介这份PPT面向电气工程专业学生及电力系统运行人员系统讲解电力系统控制的核心理论与工程实现帮助读者理解如何通过有功与无功的协调控制维持系统频率和电压稳定。资源为单个PPT文件压缩包约21.64MB内容以章节讲义形式组织涵盖概述、基本发电控制环、负荷频率控制、自动发电控制、无功与电压控制、现代控制应用等模块并配有发电机、负荷、原动机与调速系统的传递函数框图。读者可从中掌握LFC与AVR的独立分析思路、AGC对联络线功率的调节机制以及基于摆动方程的小扰动建模方法同时了解MATLAB控制工具箱与SCADA在现代能量控制中心的应用。目前已有166人学习适合作为电力系统分析课程的复习与拓展材料。1. 从一张 CHAPTER 12 的 PPT 说起电力系统控制到底在控什么很多人第一次翻到《电力系统分析》的 CHAPTER 12 电力系统控制会以为这章是讲继电保护或者开关操作其实它盯的是另一件事当负荷突然涨了 200MW或者一台机组跳了系统频率和联络线功率怎么在几十秒内被拉回来。这个动作的核心就是 AGC自动发电控制在北美教材里常写成 LFC负荷频率控制。它要解决的是有功功率的实时平衡而不是故障切除。这一章通常会把控制目标拆成三层一次调频靠调速器本地响应秒级二次调频就是 AGC分钟级负责把频率偏差和联络线交换功率偏差同时压到零三次调频是经济调度十几分钟往上。PPT 里那些框图和传递函数落到工程上就是一台调度主站给各电厂下发调节指令。适合正在学电力系统分析、准备做课程设计或者要从 MATLAB 仿真入手理解 AGC 的人。2. AGC/LFC 的控制模型与 MATLAB 建模准备2.1 单区域与多区域 LFC 的框图差异单区域模型里调速器、汽轮机、发电机和负荷构成一个闭环输入是负荷扰动 ΔP_L输出是频率偏差 Δf。多区域则要在每个区域里再加一条联络线功率偏差 ΔP_tie控制目标是 ACEArea Control ErrorACE ΔP_tie B × ΔfB 是频率偏差系数单位 MW/0.1Hz工程上一般取区域自然频率响应系数的倒数再乘一个系数。PPT 上常画成两个区域用一条联络线连起来每个区域有自己的 ACE 计算点和 PI 控制器。理解这张图的关键是AGC 不是把频率调回 50Hz 就完事它还要保证联络线交换功率回到计划值否则两个区域会互相“甩锅”。2.2 用 MATLAB 搭一个两区域 LFC 最小模型常见做法是用 Simulink 搭框图但如果你想快速验证参数用脚本写状态空间更直接。下面是一个两区域 LFC 的简化模型参数取自典型 1000MW 区域% 两区域 LFC 简化模型参数 M1 10; D1 1.0; % 区域1 惯性时间常数与阻尼 M2 10; D2 1.0; R1 0.05; R2 0.05; % 调差系数 Tg1 0.2; Tg2 0.2; % 调速器时间常数 Tt1 0.5; Tt2 0.5; % 汽轮机时间常数 B1 21; B2 21; % 频率偏差系数 MW/0.1Hz T12 0.2; % 联络线同步系数 % 状态变量: [df1, dPm1, dPg1, df2, dPm2, dPg2, dPtie] A zeros(7,7); A(1,1) -D1/M1; A(1,2) 1/M1; A(1,7) -1/M1; A(2,3) 1/Tt1; A(2,2) -1/Tt1; A(3,1) -1/(R1*Tg1); A(3,3) -1/Tg1; A(4,4) -D2/M2; A(4,5) 1/M2; A(4,7) 1/M2; A(5,6) 1/Tt2; A(5,5) -1/Tt2; A(6,4) -1/(R2*Tg2); A(6,6) -1/Tg2; A(7,1) 2*pi*T12; A(7,4) -2*pi*T12; B [0 0; 0 0; 1/Tg1 0; 0 0; 0 0; 0 1/Tg2; 0 0]; C [B1 0 0 0 0 0 1; 0 0 0 B2 0 0 -1]; % ACE1 和 ACE2 D zeros(2,2); sys ss(A,B,C,D); step(sys); % 观察 ACE 阶跃响应这段代码的逻辑是前六个状态描述两个区域的频率、机械功率和调速器阀位第七个状态是联络线功率偏差。A 矩阵里 A(1,7) 和 A(4,7) 的符号相反体现了一个区域多送电、另一个区域少送电的耦合关系。C 矩阵输出的是两个区域的 ACEB1 和 B2 的单位是 MW/0.1Hz所以 Δf 要用标幺值或 Hz 统一量纲。跑完 step 你会看到 ACE 不会自己回零因为还没有加 PI 控制器。2.3 参数表建模时最容易填错的几个量参数典型范围单位常见错误R 调差系数0.04~0.06pu误填成百分数 5B 频率偏差系数15~25MW/0.1Hz与 R 的倒数混淆Tt 汽轮机时间常数0.3~1.5s再热机组要加再热时间常数T12 联络线同步系数0.1~0.5pu忘记乘 2πACE 采样周期2~4s设成 0.1s 导致指令抖动注意B 的取值直接影响 ACE 中频率项的权重。B 偏大频率偏差被放大AGC 会频繁动作B 偏小联络线功率偏差得不到纠正。3. 从模型到指令AGC 控制器实现与 MATLAB 仿真3.1 PI 控制器的离散化写法AGC 主站里跑的是数字 PI不是连续传递函数。把上面的模型接上 PI 控制器用离散方式实现Ts 2; % AGC 采样周期 2s Kp 0.3; Ki 0.05; % 比例、积分增益 N 600; % 仿真 1200s ace1 zeros(1,N); u1 zeros(1,N); ace2 zeros(1,N); u2 zeros(1,N); int1 0; int2 0; for k 2:N % 简化用上一时刻 ACE 近似当前值 int1 int1 ace1(k-1)*Ts; int2 int2 ace2(k-1)*Ts; u1(k) -Kp*ace1(k-1) - Ki*int1; u2(k) -Kp*ace2(k-1) - Ki*int2; % 限幅机组调节容量通常 ±100MW u1(k) max(min(u1(k),100),-100); u2(k) max(min(u2(k),100),-100); % 这里应调用区域状态更新省略具体差分方程 end逻辑说明积分项用累加代替连续积分Ts 是 AGC 下发周期工程上 2~4 秒。限幅不能省否则负荷大扰动时指令会超出机组可调范围。Kp 和 Ki 的整定顺序是先加 Kp 让 ACE 快速下降再加 Ki 消除稳态误差。如果 ACE 出现等幅振荡先把 Ki 减半。3.2 用 MATLAB 画阶跃响应看控制效果把 PI 输出反馈到状态空间模型可以用lsim看 ACE 和频率偏差t 0:Ts:N*Ts; % 假设已构建闭环 Acl, Bcl, Ccl, Dcl sys_cl ss(Acl, Bcl, Ccl, Dcl, Ts); [ace, t] step(sys_cl, N); plot(t, ace(:,1), b, t, ace(:,2), r); xlabel(时间 (s)); ylabel(ACE (MW)); legend(区域1 ACE,区域2 ACE); grid on;参数说明step的第二个参数是采样点数不是时间。如果曲线在 30 秒内回到零附近说明 Kp/Ki 合适如果 100 秒还在振荡检查 T12 是否填得过大。MATLAB 的pidtune也可以直接对线性模型整定 PI但要注意它默认是连续域离散化后要重新验证。3.3 多区域 ACE 计算中的量纲坑ACE 的公式里 ΔP_tie 单位是 MWΔf 单位是 HzB 的单位是 MW/0.1Hz。很多人直接把 Δf 按标幺值代入结果 ACE 小了 50 倍。正确做法是统一成 MWΔf 用 Hz 数值乘以 10 再乘 B或者把 B 换算成 MW/Hz。仿真时可以在 C 矩阵里把 B1 写成 210 而不是 21对应 MW/Hz。这个坑在课程设计里非常常见PPT 上不会专门讲但一跑仿真就露馅。4. AGC 参数整定与常见故障排查4.1 三个必调参数Kp、Ki、BKp 决定 ACE 偏差出现时的即时调节力度典型值 0.2~0.5Ki 决定消除稳态误差的速度典型值 0.02~0.1B 决定频率项在 ACE 中的权重典型值 15~25 MW/0.1Hz。整定顺序是先设 Ki0逐步加大 Kp 直到 ACE 响应出现轻微超调然后固定 Kp加大 Ki 直到 ACE 在 60 秒内回零且不振荡。如果两个区域互相影响导致联络线功率来回摆把 B 适当调大让本区域更关注频率而不是联络线。4.2 仿真不收敛时先查这四处第一检查 A 矩阵中联络线项的符号两个区域符号必须相反。第二检查调速器时间常数 Tg 是否小于采样周期 Ts如果 Tg0.2s 而 Ts2s离散化会失真需要把 Tg 调大或减小 Ts。第三检查限幅是否设得太小导致积分饱和。第四检查 ACE 输出是否在 C 矩阵里正确组合了 Δf 和 ΔP_tie。用eig(A)看闭环极点如果有正实部说明反馈极性反了。4.3 和一次调频的配合关系AGC 指令下发后调速器先动发电机机械功率变化频率才开始回升。如果仿真里把调速器旁路掉直接让 AGC 输出控制机械功率响应会快得不真实。正确做法是保留调速器一阶惯性环节让 AGC 输出作为调速器的负荷参考值。这样频率曲线会先跌后升符合实际录波。PPT 上的框图通常把调速器和汽轮机分开画仿真时不要合并成一个增益。5. 把 CHAPTER 12 的模型用到课程设计与工程估算课程设计里常见的要求是给一个两区域系统负荷突增 100MW要求频率偏差不超过 0.2Hz联络线功率偏差在 120 秒内回零。用上面的模型先把 B 设成 21Kp0.3Ki0.05跑一遍看最大频率偏差。如果超过 0.2Hz把 Kp 加到 0.4 再跑。记录每次的调节时间和超调量画成表格对比。工程估算时可以用一个简化公式ACE 的稳态误差与 Ki 成反比调节时间约为 4/(Ki×B)。比如 Ki0.05、B21调节时间约 3.8 秒但实际因为调速器和汽轮机惯性会拉长到 60~90 秒。这个差距就是 PPT 上那些时间常数存在的意义。提示MATLAB 的step和lsim对离散系统默认零阶保持如果你用c2d转换连续模型采样周期要和 AGC 实际下发周期一致否则频率响应会对不上。最后给一个快速验证技巧把 ACE 曲线和频率偏差曲线画在同一张图上如果 ACE 回零时频率偏差还没回零说明 B 偏小如果频率偏差先回零而 ACE 还在说明联络线功率项权重不够。调 B 比调 Kp/Ki 更快让两条曲线同步收敛。本文还有配套的精品资源点击获取 电力系统AGC自动发电控制：MATLAB两区域LFC建模与仿真
简介这份PPT面向电气工程专业学生及电力系统运行人员系统讲解电力系统控制的核心理论与工程实现帮助读者理解如何通过有功与无功的协调控制维持系统频率和电压稳定。资源为单个PPT文件压缩包约21.64MB内容以章节讲义形式组织涵盖概述、基本发电控制环、负荷频率控制、自动发电控制、无功与电压控制、现代控制应用等模块并配有发电机、负荷、原动机与调速系统的传递函数框图。读者可从中掌握LFC与AVR的独立分析思路、AGC对联络线功率的调节机制以及基于摆动方程的小扰动建模方法同时了解MATLAB控制工具箱与SCADA在现代能量控制中心的应用。目前已有166人学习适合作为电力系统分析课程的复习与拓展材料。1. 从一张 CHAPTER 12 的 PPT 说起电力系统控制到底在控什么很多人第一次翻到《电力系统分析》的 CHAPTER 12 电力系统控制会以为这章是讲继电保护或者开关操作其实它盯的是另一件事当负荷突然涨了 200MW或者一台机组跳了系统频率和联络线功率怎么在几十秒内被拉回来。这个动作的核心就是 AGC自动发电控制在北美教材里常写成 LFC负荷频率控制。它要解决的是有功功率的实时平衡而不是故障切除。这一章通常会把控制目标拆成三层一次调频靠调速器本地响应秒级二次调频就是 AGC分钟级负责把频率偏差和联络线交换功率偏差同时压到零三次调频是经济调度十几分钟往上。PPT 里那些框图和传递函数落到工程上就是一台调度主站给各电厂下发调节指令。适合正在学电力系统分析、准备做课程设计或者要从 MATLAB 仿真入手理解 AGC 的人。2. AGC/LFC 的控制模型与 MATLAB 建模准备2.1 单区域与多区域 LFC 的框图差异单区域模型里调速器、汽轮机、发电机和负荷构成一个闭环输入是负荷扰动 ΔP_L输出是频率偏差 Δf。多区域则要在每个区域里再加一条联络线功率偏差 ΔP_tie控制目标是 ACEArea Control ErrorACE ΔP_tie B × ΔfB 是频率偏差系数单位 MW/0.1Hz工程上一般取区域自然频率响应系数的倒数再乘一个系数。PPT 上常画成两个区域用一条联络线连起来每个区域有自己的 ACE 计算点和 PI 控制器。理解这张图的关键是AGC 不是把频率调回 50Hz 就完事它还要保证联络线交换功率回到计划值否则两个区域会互相“甩锅”。2.2 用 MATLAB 搭一个两区域 LFC 最小模型常见做法是用 Simulink 搭框图但如果你想快速验证参数用脚本写状态空间更直接。下面是一个两区域 LFC 的简化模型参数取自典型 1000MW 区域% 两区域 LFC 简化模型参数 M1 10; D1 1.0; % 区域1 惯性时间常数与阻尼 M2 10; D2 1.0; R1 0.05; R2 0.05; % 调差系数 Tg1 0.2; Tg2 0.2; % 调速器时间常数 Tt1 0.5; Tt2 0.5; % 汽轮机时间常数 B1 21; B2 21; % 频率偏差系数 MW/0.1Hz T12 0.2; % 联络线同步系数 % 状态变量: [df1, dPm1, dPg1, df2, dPm2, dPg2, dPtie] A zeros(7,7); A(1,1) -D1/M1; A(1,2) 1/M1; A(1,7) -1/M1; A(2,3) 1/Tt1; A(2,2) -1/Tt1; A(3,1) -1/(R1*Tg1); A(3,3) -1/Tg1; A(4,4) -D2/M2; A(4,5) 1/M2; A(4,7) 1/M2; A(5,6) 1/Tt2; A(5,5) -1/Tt2; A(6,4) -1/(R2*Tg2); A(6,6) -1/Tg2; A(7,1) 2*pi*T12; A(7,4) -2*pi*T12; B [0 0; 0 0; 1/Tg1 0; 0 0; 0 0; 0 1/Tg2; 0 0]; C [B1 0 0 0 0 0 1; 0 0 0 B2 0 0 -1]; % ACE1 和 ACE2 D zeros(2,2); sys ss(A,B,C,D); step(sys); % 观察 ACE 阶跃响应这段代码的逻辑是前六个状态描述两个区域的频率、机械功率和调速器阀位第七个状态是联络线功率偏差。A 矩阵里 A(1,7) 和 A(4,7) 的符号相反体现了一个区域多送电、另一个区域少送电的耦合关系。C 矩阵输出的是两个区域的 ACEB1 和 B2 的单位是 MW/0.1Hz所以 Δf 要用标幺值或 Hz 统一量纲。跑完 step 你会看到 ACE 不会自己回零因为还没有加 PI 控制器。2.3 参数表建模时最容易填错的几个量参数典型范围单位常见错误R 调差系数0.04~0.06pu误填成百分数 5B 频率偏差系数15~25MW/0.1Hz与 R 的倒数混淆Tt 汽轮机时间常数0.3~1.5s再热机组要加再热时间常数T12 联络线同步系数0.1~0.5pu忘记乘 2πACE 采样周期2~4s设成 0.1s 导致指令抖动注意B 的取值直接影响 ACE 中频率项的权重。B 偏大频率偏差被放大AGC 会频繁动作B 偏小联络线功率偏差得不到纠正。3. 从模型到指令AGC 控制器实现与 MATLAB 仿真3.1 PI 控制器的离散化写法AGC 主站里跑的是数字 PI不是连续传递函数。把上面的模型接上 PI 控制器用离散方式实现Ts 2; % AGC 采样周期 2s Kp 0.3; Ki 0.05; % 比例、积分增益 N 600; % 仿真 1200s ace1 zeros(1,N); u1 zeros(1,N); ace2 zeros(1,N); u2 zeros(1,N); int1 0; int2 0; for k 2:N % 简化用上一时刻 ACE 近似当前值 int1 int1 ace1(k-1)*Ts; int2 int2 ace2(k-1)*Ts; u1(k) -Kp*ace1(k-1) - Ki*int1; u2(k) -Kp*ace2(k-1) - Ki*int2; % 限幅机组调节容量通常 ±100MW u1(k) max(min(u1(k),100),-100); u2(k) max(min(u2(k),100),-100); % 这里应调用区域状态更新省略具体差分方程 end逻辑说明积分项用累加代替连续积分Ts 是 AGC 下发周期工程上 2~4 秒。限幅不能省否则负荷大扰动时指令会超出机组可调范围。Kp 和 Ki 的整定顺序是先加 Kp 让 ACE 快速下降再加 Ki 消除稳态误差。如果 ACE 出现等幅振荡先把 Ki 减半。3.2 用 MATLAB 画阶跃响应看控制效果把 PI 输出反馈到状态空间模型可以用lsim看 ACE 和频率偏差t 0:Ts:N*Ts; % 假设已构建闭环 Acl, Bcl, Ccl, Dcl sys_cl ss(Acl, Bcl, Ccl, Dcl, Ts); [ace, t] step(sys_cl, N); plot(t, ace(:,1), b, t, ace(:,2), r); xlabel(时间 (s)); ylabel(ACE (MW)); legend(区域1 ACE,区域2 ACE); grid on;参数说明step的第二个参数是采样点数不是时间。如果曲线在 30 秒内回到零附近说明 Kp/Ki 合适如果 100 秒还在振荡检查 T12 是否填得过大。MATLAB 的pidtune也可以直接对线性模型整定 PI但要注意它默认是连续域离散化后要重新验证。3.3 多区域 ACE 计算中的量纲坑ACE 的公式里 ΔP_tie 单位是 MWΔf 单位是 HzB 的单位是 MW/0.1Hz。很多人直接把 Δf 按标幺值代入结果 ACE 小了 50 倍。正确做法是统一成 MWΔf 用 Hz 数值乘以 10 再乘 B或者把 B 换算成 MW/Hz。仿真时可以在 C 矩阵里把 B1 写成 210 而不是 21对应 MW/Hz。这个坑在课程设计里非常常见PPT 上不会专门讲但一跑仿真就露馅。4. AGC 参数整定与常见故障排查4.1 三个必调参数Kp、Ki、BKp 决定 ACE 偏差出现时的即时调节力度典型值 0.2~0.5Ki 决定消除稳态误差的速度典型值 0.02~0.1B 决定频率项在 ACE 中的权重典型值 15~25 MW/0.1Hz。整定顺序是先设 Ki0逐步加大 Kp 直到 ACE 响应出现轻微超调然后固定 Kp加大 Ki 直到 ACE 在 60 秒内回零且不振荡。如果两个区域互相影响导致联络线功率来回摆把 B 适当调大让本区域更关注频率而不是联络线。4.2 仿真不收敛时先查这四处第一检查 A 矩阵中联络线项的符号两个区域符号必须相反。第二检查调速器时间常数 Tg 是否小于采样周期 Ts如果 Tg0.2s 而 Ts2s离散化会失真需要把 Tg 调大或减小 Ts。第三检查限幅是否设得太小导致积分饱和。第四检查 ACE 输出是否在 C 矩阵里正确组合了 Δf 和 ΔP_tie。用eig(A)看闭环极点如果有正实部说明反馈极性反了。4.3 和一次调频的配合关系AGC 指令下发后调速器先动发电机机械功率变化频率才开始回升。如果仿真里把调速器旁路掉直接让 AGC 输出控制机械功率响应会快得不真实。正确做法是保留调速器一阶惯性环节让 AGC 输出作为调速器的负荷参考值。这样频率曲线会先跌后升符合实际录波。PPT 上的框图通常把调速器和汽轮机分开画仿真时不要合并成一个增益。5. 把 CHAPTER 12 的模型用到课程设计与工程估算课程设计里常见的要求是给一个两区域系统负荷突增 100MW要求频率偏差不超过 0.2Hz联络线功率偏差在 120 秒内回零。用上面的模型先把 B 设成 21Kp0.3Ki0.05跑一遍看最大频率偏差。如果超过 0.2Hz把 Kp 加到 0.4 再跑。记录每次的调节时间和超调量画成表格对比。工程估算时可以用一个简化公式ACE 的稳态误差与 Ki 成反比调节时间约为 4/(Ki×B)。比如 Ki0.05、B21调节时间约 3.8 秒但实际因为调速器和汽轮机惯性会拉长到 60~90 秒。这个差距就是 PPT 上那些时间常数存在的意义。提示MATLAB 的step和lsim对离散系统默认零阶保持如果你用c2d转换连续模型采样周期要和 AGC 实际下发周期一致否则频率响应会对不上。最后给一个快速验证技巧把 ACE 曲线和频率偏差曲线画在同一张图上如果 ACE 回零时频率偏差还没回零说明 B 偏小如果频率偏差先回零而 ACE 还在说明联络线功率项权重不够。调 B 比调 Kp/Ki 更快让两条曲线同步收敛。本文还有配套的精品资源点击获取