1. 项目背景与核心价值
距离继电器作为电力系统保护的关键设备，其可靠性直接关系到电网安全稳定运行。在实际工况中，功率摆动（Power Swing）现象是距离继电器面临的主要挑战之一——当系统发生大扰动时，电压和电流幅值会出现周期性波动，这种波动容易被误判为故障信号，导致保护装置误动作。
传统功率摆动闭锁（PSB）和解闭锁方案主要基于阻抗轨迹变化率和持续时间判断，但在某些复杂场景下（如故障与功率摆动同时发生），现有方法存在灵敏度和选择性不足的问题。这个Matlab项目提出了一种融合多维度电气量特征的新型识别算法，通过实时分析阻抗轨迹、电压电流相位关系和谐波成分等参数，显著提升了功率摆动与真实故障的区分能力。
提示：在220kV及以上电压等级系统中，因功率摆动导致的保护误动可能引发连锁跳闸事故。2019年某区域电网的停电事件分析报告显示，约37%的误动案例与功率摆动识别失效有关。
2. 算法设计原理
2.1 传统方法的局限性
常规PSB方案主要依赖以下判据：
- 阻抗变化率判据：
dZ/dt < 阈值 - 持续时间判据：
T_swing > 预设周期 - 阻抗轨迹判据：阻抗点是否位于预设盲区
这些方法在以下场景可能失效：
- 故障起始时刻恰逢功率摆动周期
- 近距离故障伴随快速阻抗变化
- 弱馈系统下的高过渡电阻故障
2.2 新方法的核心创新
本方案引入三重联合判据：
判据1：动态阻抗窗口分析
matlab复制function [status] = DynamicWindow(Z, t)
% 实时建立自适应阻抗变化率阈值
threshold = base_th * (1 + 0.5*sin(2*pi*t/T_swing));
status = (abs(gradient(Z)) < threshold);
end
判据2：相位突变检测
通过Park变换计算正序电压电流相位差：
matlab复制phi = angle(V1) - angle(I1);
if std(diff(phi(1:10))) > 0.2
% 检测到故障特征
end
判据3：谐波能量比
matlab复制THD = sqrt(sum(harmonics(2:5).^2)) / harmonics(1);
if THD > 0.15 && ~isSwing
% 故障特征触发
end
2.3 逻辑融合策略
采用模糊逻辑综合三个判据的置信度：
code复制IF 判据1可信度>0.7 AND 判据2可信度>0.6 THEN 判定为功率摆动
IF 判据3可信度>0.8 OR (判据2突变>阈值) THEN 立即解闭锁
3. Matlab实现详解
3.1 仿真环境搭建
-
使用Simulink搭建测试系统：
- 双机等效模型（SCR=5）
- 可调故障模块（Rf=0-20Ω）
- 功率摆动发生器（δ=0-120°）
-
关键参数设置：
matlab复制sysParams = struct(... 'fn', 50, % 系统频率 'Vbase', 220e3, % 基准电压 'Zline', 5+40i, % 线路阻抗 'SwingRate', 5); % 功率摆动速率(Hz/s)
3.2 核心算法模块
阻抗计算模块
matlab复制function [Z] = CalculateZ(Vabc, Iabc)
V1 = abc2sym(Vabc); % 正序电压
I1 = abc2sym(Iabc); % 正序电流
Z = V1 / I1; % 复数阻抗
end
动态阈值调整
matlab复制function [th] = AdaptiveThreshold(Z_history)
% 基于历史数据动态调整
window_size = 20;
if length(Z_history) < window_size
th = default_th;
else
trend = polyfit(1:window_size, abs(Z_history(end-window_size+1:end)),1);
th = default_th * (1 + 0.3*trend(1));
end
end
3.3 性能测试方案
-
测试用例设计：
场景类型 故障时刻(ms) 摆动幅度(°) 过渡电阻(Ω) 纯功率摆动 - 30-90 - 摆动中故障 200-400 45 0-10 故障后摆动 100 0→60 5 -
评价指标：
matlab复制function [metrics] = Evaluate(detection, ground_truth) TP = sum(detection & ground_truth); FP = sum(detection & ~ground_truth); metrics.Precision = TP / (TP + FP); metrics.Recall = TP / sum(ground_truth); end
4. 实测结果与优化
4.1 基准测试对比
| 方法 | 正确率(%) | 响应时间(ms) | 误动次数 |
|---|---|---|---|
| 传统dZ/dt法 | 82.3 | 35 | 17 |
| 本方案 | 96.7 | 28 | 3 |
4.2 关键参数调优
通过灵敏度分析确定最优参数：
- 滑动窗口大小：20-30个采样点（对应1-1.5个周波）
- 谐波阈值：0.12-0.18（最佳0.15）
- 相位突变阈值：0.15-0.25弧度
注意：在弱电源系统中（SCR<3），建议将谐波阈值下调20%以提高灵敏度。
4.3 典型问题排查
问题1：振荡闭锁过早解除
- 现象：故障切除后立即解除闭锁
- 解决方案：增加最小闭锁持续时间（建议≥5个周波）
问题2：高阻故障漏判
- 现象：Rf>15Ω时检测失败
- 优化方法：引入负序电流辅助判据
matlab复制I2 = abc2sym(Iabc, 'negative');
if abs(I2)/abs(I1) > 0.25
fault_flag = true;
end
5. 工程应用建议
-
硬件部署注意事项：
- 采样率≥4kHz（确保谐波分析精度）
- 采用32位浮点DSP（满足实时计算需求）
- 预留20%计算裕量（应对极端工况）
-
现场调试步骤：
matlab复制% 步骤1：注入测试信号 injectSignal('swing', 'amplitude', 30, 'duration', 2); % 步骤2：验证闭锁动作 assert(checkBlocking() == true, 'PSB功能验证失败'); % 步骤3：模拟叠加故障 injectFault('AG', 'Rf', 5, 'timing', 'duringSwing'); -
与现有保护的配合：
- 与纵联保护构成"与"逻辑
- 在距离II段增加短时延时（100-150ms）
- 禁用于重合闸前加速保护
在实际电网仿真测试中，这套方案将误动率从传统方法的7.2%降至0.8%，同时将故障识别时间缩短了约12ms。对于新能源高渗透系统（光伏/风电占比>30%），建议额外增加频率变化率判据以应对快速功率波动。