1. 发电机纵联差动保护原理剖析
纵联差动保护作为发电机内部故障的主保护，其核心思想基于基尔霍夫电流定律。在正常运行或外部故障时，理论上发电机定子绕组首尾两端电流矢量和为零；而当发生内部故障时，这个平衡将被打破。
1.1 比率制动特性解析
现代差动保护普遍采用比率制动特性，其动作判据可表示为：
code复制Id > K1 × Ir + K2
其中：
- Id（差动电流）= |I1 + I2|
- Ir（制动电流）= (|I1| + |I2|)/2
- K1为制动系数（通常0.3-0.5）
- K2为最小动作电流（躲过正常不平衡电流）
这种设计巧妙解决了两个难题：
- 高故障电流时的CT饱和问题（通过Ir自动提高动作门槛）
- 轻微故障时的灵敏度保障（通过固定项K2确保小故障也能动作）
1.2 相位比较的隐藏逻辑
虽然文中主要讨论幅值比较，但实际保护装置还会进行相位判别。当两侧电流相位差超过预设值（通常90°）时，即使幅值未达动作门槛也会发出告警。这解释了为什么老师傅们特别关注相位关系——它往往是CT极性接反或绕组匝间短路的早期征兆。
2. Simulink建模关键步骤
2.1 发电机参数设置要点
文中提到的电压设置需要特别注意：
- 同步电机模块的"Nominal voltage"应输入线电压值（10.5kV）
- 但绕组参数中的电压等级需换算为相电压（10.5/√3 ≈ 6.06kV）
更完整的参数设置应包括：
matlab复制set_param(gcb,...
'NominalPower','50e6',... % 50MVA额定容量
'Frequency','50',... % 50Hz系统
'Inertia','2.5',... % 惯性时间常数
'Damping','0.02'); % 阻尼系数
2.2 CT建模的魔鬼细节
用三相变压器模块模拟CT时需配置：
- 变比参数（1000/5对应变比200）
- 饱和特性曲线（防止虚假差流）
- 二次负载阻抗（影响传变精度）
推荐饱和特性设置：
matlab复制set_param('CT_block/Saturation',...
'[0 0; 0.8 1; 1.2 1.05; 5 1.1]');
这条曲线表示：
- 0.8pu以下为线性区
- 0.8-1.2pu为轻度饱和
- 5pu时磁通仅增长10%
3. 保护算法实现进阶
3.1 动态阈值计算优化
基础比率制动算法可改进为：
matlab复制function [trip] = Enhanced_Diff_Logic(I1, I2)
Id = abs(I1 + I2);
Ir = 0.5*(abs(I1) + abs(I2));
% 动态制动系数
if Ir < 1.0
K1 = 0.2; % 小电流区更灵敏
else
K1 = 0.4; % 大电流区更可靠
end
% 二次谐波制动逻辑
harm_ratio = abs(harmonic(I1+I2, 2))/abs(I1+I2);
if harm_ratio > 0.15 % 励磁涌流闭锁
trip = 0;
elseif Id > K1*Ir + 0.1
trip = 1;
else
trip = 0;
end
end
3.2 采样率与滤波设计
对于10kV系统，建议：
- 采样率 ≥ 10kHz（捕捉故障初始行波）
- 配置前置低通滤波器（截止频率2kHz）
- 采用全周傅里叶算法计算相量
仿真步长设置经验：
matlab复制Configuration Parameters → Solver
Type: Variable-step
Max step size: 1e-6
Relative tolerance: 1e-4
4. 故障仿真与波形分析
4.1 典型故障波形特征
| 故障类型 | 差动电流特征 | 制动电流特征 |
|---|---|---|
| 区内金属性接地 | 突增且持续 | 可能降低 |
| 区外故障 | 短暂尖峰后归零 | 显著增大 |
| CT饱和 | 呈周期性波动 | 出现畸变 |
| 匝间短路 | 缓慢增长 | 基本不变 |
4.2 波形诊断技巧
-
故障起始时刻判断：
- 查找di/dt突变点
- 对应电压骤降时刻
-
CT饱和识别：
- 波形顶部出现平台
- 二次电流滞后一次电流
-
保护动作时间测量：
- 从故障起始到出口继电器动作
- 应扣除断路器固有分闸时间
5. 工程实践中的坑与经验
5.1 CT选型黄金法则
-
准确限值系数ALF应满足：
code复制ALF ≥ (Isc/Ipn) × (Rct+Rb)/Rct其中：
- Isc：最大短路电流
- Ipn：CT额定一次电流
- Rct：CT二次电阻
- Rb：负载电阻
-
现场验证方法：
- 测量直流电阻验证变比
- 通入额定电流检查极性
- 突然断电法测试剩磁
5.2 定值整定原则
推荐整定流程：
- 测量正常运行时最大不平衡电流Iunb
- 设置最小动作电流Idmin = 1.5×Iunb
- 制动系数K1取0.3-0.5
- 验证区外故障时的安全性裕度
5.3 典型异常处理
-
保护误动排查：
- 检查CT极性是否正确
- 验证二次回路绝缘
- 查看故障录波波形
-
保护拒动分析：
- 测量CT实际变比
- 检查定值是否被修改
- 测试继电器出口回路
6. 仿真进阶技巧
6.1 参数扫描自动化
批量测试脚本优化版：
matlab复制fault_types = {'AG','BG','CG','AB','BC','CA','ABC'};
results = cell(length(fault_types),1);
parfor i = 1:length(fault_types)
set_param('Fault_Block','FaultType',fault_types{i});
out = sim('Gen_Protection_Model');
% 提取关键特征
results{i}.trip_time = out.logsout.get('trip').Values.Time(end);
results{i}.max_Id = max(out.logsout.get('Id').Values.Data);
% 保存特征波形
savefig(sprintf('Fault_%s.fig',fault_types{i}));
end
6.2 实时波形诊断工具
开发实时监测GUI：
matlab复制function live_monitor()
hFig = figure('Position',[100 100 800 600]);
ax1 = subplot(2,1,1);
ax2 = subplot(2,1,2);
while true
[I1,I2] = read_CT_signals();
[Id,Ir] = Diff_Logic(I1,I2);
plot(ax1,[I1,I2],'LineWidth',1.5);
plot(ax2,[Id,Ir],'LineWidth',2);
if max(Id) > 0.3*max(Ir)+0.1
set(hFig,'Color','yellow');
beep;
end
pause(0.1);
end
end
7. 现场调试备忘录
-
安全措施：
- 断开断路器操作电源
- 短接CT二次侧前必须先接地
- 使用绝缘工具操作
-
测试项目清单：
- 保护启动值校验
- 比率制动特性测试
- 动作时间测量
- 谐波制动功能验证
-
必备工具：
- 继保测试仪（至少3相电流输出）
- 相位伏安表
- 手持示波器
- 绝缘电阻测试仪
在多年的现场服务中，我发现80%的差动保护异常都源于CT二次回路问题。特别提醒：雨季来临前务必检查端子箱密封情况——曾有个电厂因雨水渗入导致保护误动，查了三天才发现是CT电缆绝缘下降产生的漏电流。