1. 项目背景与核心目标
在电力系统运行中，暂态稳定性分析是确保电网安全的关键环节。当系统发生短路故障时，发电机的机械功率与电磁功率平衡被打破，转子转速会出现波动。这个项目聚焦于两相接地短路故障场景下，发电机转速变化的动态过程仿真。
我选择Simulink作为仿真平台，主要考虑到它在电力系统建模方面的独特优势：直观的图形化界面、丰富的电力元件库、以及强大的数值计算能力。通过搭建单机无穷大系统模型，可以清晰地观察到故障发生、切除及后续恢复过程中，发电机转速的振荡特性和稳定状态。
2. 系统建模与参数设置
2.1 单机无穷大系统架构
典型的单机无穷大系统由以下核心组件构成：
- 同步发电机（采用二阶实用模型）
- 变压器（额定容量匹配发电机）
- 输电线路（π型等效电路）
- 无穷大母线（电压幅值相位恒定）
- 故障模块（设置在两相线路间）
在Simulink中搭建模型时，我特别注意了几个关键连接点：
- 发电机出口电压需要经过变压器升压
- 线路阻抗参数需换算为标幺值
- 测量模块需正确接入转速信号输出端
2.2 发电机参数配置示例
% 同步发电机参数（标幺值） H = 3.5; % 惯性常数 D = 2; % 阻尼系数 Xd = 1.8; % d轴同步电抗 Xq = 1.7; % q轴同步电抗 X'd = 0.3; % d轴暂态电抗 T'do = 8.0; % d轴开路暂态时间常数重要提示：参数设置必须采用标幺值系统，否则仿真结果将失去可比性。建议先在工作区定义基值，再计算各元件参数。
3. 故障场景实现与仿真
3.1 两相接地短路设置
在Simulink中实现AB相接地短路需要：
- 使用Three-Phase Fault模块
- 设置Fault A和Fault B为接地状态
- 配置故障时序：
- 起始时间：1.0秒（系统已进入稳态）
- 持续时间：0.1秒（典型故障切除时间）
我特别添加了故障电阻参数（通常取0.001pu），以模拟实际电弧电阻。过小的电阻值可能导致数值计算困难，出现仿真发散问题。
3.2 仿真参数配置要点
| 参数项 | 推荐值 | 说明 |
|---|---|---|
| Solver | ode23tb | 适合电力系统刚性方程 |
| Max step size | 0.01s | 捕捉转速快速变化 |
| Relative error | 1e-4 | 平衡精度与计算效率 |
| Stop time | 10s | 观察完整动态过程 |
4. 结果分析与典型波形
4.1 转速变化动态过程
故障期间转速变化呈现典型的三阶段特征：
- 加速阶段（故障发生瞬间）
- 电磁功率突降导致机械功率过剩
- 转速以近似线性增长
- 振荡阶段（故障切除后）
- 出现衰减振荡（频率约0.5-2Hz）
- 振幅受阻尼系数显著影响
- 稳定阶段（约5秒后）
- 恢复同步运行或失步
4.2 关键指标测量
通过Simulink示波器可获取：
- 最大转速偏差（Δω_max）
- 振荡周期（T_osc）
- 稳定时间（t_settling）
我习惯使用To Workspace模块导出数据，再用MATLAB脚本进行定量分析：
[peak,loc] = findpeaks(omega); damping_ratio = log(peak(1)/peak(2))/sqrt(4*pi^2 + log(peak(1)/peak(2))^2);5. 工程经验与调试技巧
5.1 常见问题排查
- 仿真发散问题：
- 检查发电机初始条件是否合理
- 尝试减小步长或更换求解器
- 异常振荡：
- 确认变压器接线组别正确
- 检查线路参数单位是否统一
5.2 参数敏感性分析
通过参数扫描发现：
- 惯性常数H每增加1秒，最大转速偏差减小约12%
- 阻尼系数D超过4后对改善稳定性效果有限
- 故障切除时间每增加0.01秒，恢复时间延长约15%
建议搭建参数化模型，使用Simulink的Batch Simulation功能进行多场景分析。
6. 模型扩展与应用
基础模型可通过以下方式增强：
- 添加AVR（自动电压调节器）
- 引入PSS（电力系统稳定器）
- 考虑多机系统交互影响
我在实际项目中发现，加入调速器模型后，转速的二次波动会明显减弱。这可以通过修改原动机功率输出来实现：
function Pm = governor(omega) omega_ref = 1.0; R = 0.05; % 调差系数 Pm = max(0, 1 - (omega-omega_ref)/R); end对于需要更高精度的场景，建议采用六阶发电机模型，但会显著增加计算复杂度。在大多数暂态稳定性分析中，二阶模型已经能够满足工程精度要求。
