35kV电力系统三段式电流保护Matlab仿真指南-代码聚汇网

35kV电力系统三段式电流保护Matlab仿真指南

光启元

1. 35kV电力系统三段式电流保护仿真实战

作为一名在电力系统保护领域摸爬滚打十年的工程师，我深刻理解三段式电流保护作为电力系统第一道防线的重要性。今天就用Matlab/Simulink带大家完整走一遍35kV线路的仿真建模过程，这个电压等级在工业园区和分布式电网中非常典型。

先说说为什么选择这个案例：35kV系统既有足够的复杂度（需要考虑分布式电容电流、线路阻抗等参数），又不像110kV以上系统那样需要引入距离保护等更复杂的保护方案。通过这个仿真，你不仅能掌握Simulink电力系统建模的核心技巧，更能透彻理解I段速断、II段限时速断和III段过电流保护的配合逻辑。

2. 仿真模型构建与参数设定

2.1 系统单线图与等效模型

我们先构建一个典型的35kV放射状供电网络：

电源等效阻抗：Xs=0.2Ω（系统短路容量约200MVA）
线路参数：LGJ-240导线，长度15km
- 正序阻抗Z1=0.17+j0.38Ω/km
- 零序阻抗Z0=0.23+j1.72Ω/km
负载：8MW工业负荷（功率因数0.85）

在Simulink中打开Simscape Electrical库，按这个顺序搭建：

三相电压源（设置35kV/√3相电压）
串联RL支路模拟系统阻抗
分布参数线路模型（Bergeron模型）
三相并联负载

关键技巧：线路模型一定要用Frequency-Dependent模式，否则高频暂态过程会失真。双击线路模块，在"Model"选项中选择"Frequency-dependent (phase)"。

2.2 保护CT/PT配置

电流互感器选择需要注意：

变比：400/5（根据负载电流8MW/35kV≈132A，取1.5倍裕度）
准确限值系数：10（确保短路电流下不饱和）
在Simulink中用Current Measurement模块实现，记得勾选"谐波测量"选项。

电压互感器直接用Voltage Measurement模块，变比35000/100。保护装置采样率建议设为每周波64点（对应3200Hz采样率）。

3. 三段式保护算法实现

3.1 I段速断保护建模

I段保护的核心是躲过本线路末端最大短路电流：

matlab复制function [trip] = Stage1_Protection(Ia, Ib, Ic)
    I_set = 1800; % 动作值设为1800A（约4.5倍额定电流）
    I_fault = max([abs(Ia), abs(Ib), abs(Ic)]);
    trip = (I_fault > I_set); 
end

这个函数可以直接封装成Matlab Function模块。动作时间设为0秒，但实际会有约20ms的继电器固有延时。

3.2 II段限时速断保护

II段要保证与下级保护配合，典型配置：

动作值：I段定值的80%（约1440A）
延时：0.3-0.5秒（这里取0.4秒）

Simulink实现要点：

用Timer模块实现延时逻辑
加入方向元件判断（用电压电流相位比较）
配置闭锁逻辑（通过二进制信号交互）

3.3 III段过电流保护

作为后备保护，III段需要：

动作值：1.2倍最大负荷电流（约160A）
反时限特性：使用IEC标准反时限曲线

matlab复制function [trip] = Stage3_Protection(I, TMS)
    % TMS（时间倍数）设为0.1
    k = 0.14 / ( (I/160)^0.02 - 1 );
    trip_delay = k * TMS; 
    % 用persistent变量实现计时器
    persistent timer;
    if isempty(timer), timer = 0; end
    
    if I > 160
        timer = timer + 1/3200; % 假设采样率3200Hz
        trip = (timer >= trip_delay);
    else
        timer = 0;
        trip = false;
    end
end

4. 典型故障仿真与分析

4.1 三相短路测试

在线路5km处设置0.1秒时三相金属性短路：

故障电流波形显示峰值约2.8kA
I段保护在22ms后动作
断路器分闸时间建模为60ms

关键观察点：

用Powergui做FFT分析，可见二次谐波含量<5%（验证CT未饱和）
电压跌落至15%额定值（符合35kV系统特性）

4.2 高阻接地故障模拟

设置经100Ω过渡电阻的单相接地：

零序电流约85A（需在保护配置中设3I0>40A启动）
III段保护在1.2秒后动作
用Sequence Analyzer模块查看序分量

重要经验：高阻故障时，传统过流保护可能失效，这时需要配置零序方向保护作为补充。在Simulink中可以通过比较零序电压和电流相位来实现。

5. 保护配合验证与优化

5.1 时间-电流特性曲线绘制

通过批量仿真获取动作时间数据：

matlab复制current_levels = 100:100:3000;
trip_times = zeros(size(current_levels));

for i = 1:length(current_levels)
    set_param('ProtectionModel/Fault', 'R', num2str(35000/sqrt(3)/current_levels(i)));
    simOut = sim('ProtectionModel');
    trip_times(i) = simOut.tripTime(end);
end

semilogx(current_levels, trip_times); 
grid on;
xlabel('Fault Current (A)');
ylabel('Trip Time (s)');

这条曲线应该显示清晰的三个阶段特性，且与理论计算误差<5%。

5.2 保护灵敏度校验

验证保护启动值：

最小运行方式下（系统阻抗最大），线路末端短路电流应仍大于I段定值
校验公式：Ksen = I_fault.min / I_set > 1.2
在Simulink中通过逐步减小电源阻抗来寻找临界点

6. 工程实践中的坑与经验

CT饱和问题：当短路电流超过4000A时（对应二次侧50A），仿真显示波形畸变严重。解决方案：
- 改用TPY级CT
- 在算法中加入二次谐波制动
分布式电容影响：长线路充电电流可能导致保护误动。实测数据：
- 15km线路电容电流约3.5A
- 需要在零序保护中设置>5A的启动值
Simulink仿真提速技巧：
- 使用变步长solver：ode23tb
- 关闭所有scope的"Limit data points"选项
- 对非关键部件用简化模型
实际工程与仿真的差异补偿：
- 断路器动作时间增加20ms余量
- 考虑CT传变误差（在算法中加入±5%的裕度）
- 气候条件影响（通过调整线路参数±15%做边界测试）

7. 模型扩展与进阶应用

这个基础模型还可以进一步开发：

加入自动重合闸逻辑（Simulink中用Stateflow实现）
构建双端电源系统研究方向保护
与RTDS等实时仿真器做联合仿真
生成C代码部署到实际保护装置（用Simulink Coder）

对于想深入研究的同行，我特别推荐用这个模型做以下实验：

研究过渡电阻对保护性能的影响规律
分析分布式电源接入对保护配合的影响
开发基于波形畸变识别的故障类型判别算法