1. 单机无穷大系统暂态稳定性仿真研究概述
电力系统暂态稳定性分析是保障电网安全运行的重要技术手段。单机无穷大系统作为研究电力系统动态行为的经典模型，通过将复杂电网简化为单台发电机与无穷大母线相连的系统，能够清晰地展现发电机在故障扰动下的动态响应特性。本次研究聚焦于两相接地短路故障场景，采用理论计算与Simulink仿真相结合的方法，深入分析故障切除时间对发电机转速变化的影响规律。
在实际工程应用中，这种仿真分析具有重要价值。以某区域电网为例，2020年曾发生过因线路故障导致的多台机组失稳事件，事后分析发现故障切除时间超出临界值是主要原因。通过此类仿真研究，可以预先确定各类故障下的极限切除时间，为继电保护整定提供理论依据。
2. 理论基础与关键参数计算
2.1 系统建模原理
单机无穷大系统的核心假设在于"无穷大母线"具有无限大的惯性和恒定的电压频率。这种简化处理使得我们可以将关注点集中在待研究发电机的动态行为上。系统等效电路如图1所示，其中：
- Xd' 表示发电机暂态电抗
- XT 为变压器阻抗
- XL 是线路电抗
- V∞∠0° 代表无穷大母线电压
重要提示：在实际建模时，所有参数都需要换算到统一基准值下，通常选择系统额定容量SB=100MVA，电压基准取短路点平均额定电压。
2.2 暂态稳定判据
等面积定则是判断暂态稳定性的核心理论工具。其物理意义是：故障期间转子获得的加速面积（A1）必须小于故障切除后的减速面积（A2）。通过求解以下方程可以得到极限切除角δcr：
∫δ0^δcr (Pm - Pe1)dδ = ∫δcr^δmax (Pe2 - Pm)dδ
其中：
- Pm 为机械功率（恒定值）
- Pe1 是故障期间电磁功率
- Pe2 为故障后电磁功率
- δ0 是初始功角
2.3 两相接地故障特性分析
两相接地短路（假设为BC相接地）时，系统将出现不对称运行状态。采用对称分量法分析，需要建立正序、负序和零序网络。故障点的边界条件为：
- VB = VC = 0
- IA = 0
这种故障类型会产生较大的负序电流，导致发电机转子产生附加的发热和振动，对机组安全构成威胁。仿真时需要特别注意故障期间的电抗变化：
Xeq = X1 + (X2 || X0)
其中X1、X2、X0分别为正序、负序和零序电抗。
3. Simulink建模与参数设置
3.1 整体模型架构
基于Simulink搭建的仿真模型包含以下核心模块：
- 同步发电机模块（Synchronous Machine pu Standard）
- 励磁系统（Excitation System）
- 变压器模块（Three-Phase Transformer）
- 输电线路（Three-Phase PI Section Line）
- 故障模块（Three-Phase Fault）
- 测量与显示模块（Scope, Display）
模型连接顺序为：发电机→变压器→线路→无穷大母线，在变压器高压侧设置故障点。
3.2 关键参数配置示例
以某300MW机组为例，典型参数设置如下：
matlab复制% 发电机参数（标幺值）
H = 3.5; % 惯性常数(s)
Xd = 1.8; % 直轴同步电抗
Xq = 1.7; % 交轴同步电抗
Xd' = 0.3; % 直轴暂态电抗
Tdo' = 5.0; % 直轴暂态开路时间常数
% 变压器参数
XT = 0.12; % 阻抗电压百分比(换算为标幺值)
% 线路参数
XL = 0.2; % 100km线路电抗(标幺值)
3.3 仿真步长设置技巧
暂态稳定仿真对步长选择非常敏感，建议采用变步长求解器ode23tb（适用于刚性系统），按以下原则设置：
- 最大步长不超过0.01s
- 相对容差1e-4
- 绝对容差1e-6
对于包含电力电子元件的系统，可能需要将最大步长减小到0.001s以获得准确结果。
4. 仿真结果分析与讨论
4.1 正常工况验证
在投入故障模块前，首先验证系统稳态运行参数：
- 发电机端电压：1.02pu
- 输出功率：0.95pu（有功）
- 功角：约35°
- 转速：严格保持1pu（3000rpm）
这个步骤常被初学者忽视，但却是确保模型正确的关键检查点。
4.2 故障场景对比分析
设置两相接地短路发生在t=1s，比较两种切除时间下的转速响应：
| 工况类型 | 切除时间(s) | 最大转速偏差(%) | 稳定时间(s) | 最终状态 |
|---|---|---|---|---|
| 及时切除 | 0.15 | +1.2 | 3.5 | 稳定 |
| 延迟切除 | 0.25 | +8.7 | - | 失步 |
典型转速响应曲线特征：
- 及时切除：呈现衰减振荡，3-5个周期后恢复同步
- 延迟切除：单调递增，最终触发过速保护
4.3 灵敏度分析
通过参数扫描研究关键因素的影响程度：
-
惯性常数H：
- H增大→极限切除时间延长
- 但H>5s后改善效果递减
-
故障位置：
- 近发电机端故障最严重
- 线路中点故障比末端轻微约15%
-
初始负载率：
- 重载工况（P>0.9pu）下稳定裕度降低30-40%
5. 工程应用与问题排查
5.1 保护整定参考
根据仿真结果，建议继电保护设置：
- 两相接地故障：切除时间≤0.85tcr
- 考虑1.2倍安全裕度
- 配合断路器固有动作时间（通常50-80ms）
5.2 常见问题解决方案
问题1：仿真出现代数环错误
- 检查所有测量模块是否都加了滤波器
- 在适当位置插入Transport Delay模块
问题2：转速曲线异常平滑
- 确认是否启用了详细的发电机模型
- 检查阻尼系数D是否设置合理
问题3：故障后电压无法恢复
- 验证励磁系统参数
- 检查AVR的PID调节器设置
5.3 模型验证方法
为确保仿真结果可信，建议采用三重验证：
- 理论计算验证：对比等面积法则计算结果
- 商业软件验证：如PSASP、PSS/E等
- 现场数据对比：利用PMU录波数据校验
6. 扩展研究方向
本次基础研究可以进一步深化：
-
考虑励磁系统动态响应
- 添加IEEE AC1A型励磁模型
- 研究电压调节对稳定的影响
-
引入调速器模型
- 模拟一次调频过程
- 分析机械功率变化效应
-
复杂故障场景
- 发展性故障仿真
- 重合闸过程模拟
-
多机系统扩展
- 构建3机9节点系统
- 研究机组间的动态交互
在实际工程应用中，我们发现当系统初始功角超过70°时，即使很小的扰动也可能导致失稳。这提示我们在电网运行中需要严格控制传输功率，保留足够的稳定裕度。对于重要输电通道，建议配置双重化快速保护，确保故障切除时间控制在100ms以内。