1. 项目背景与核心价值
小电流接地系统是配电网中最常见的接地方式之一，主要包括中性点不接地和经消弧线圈接地两种形式。这类系统在发生单相接地故障时，故障电流较小，系统可以继续运行1-2小时，大大提高了供电可靠性。但这也带来了故障检测和定位的难题。
通过Simulink搭建这两种接地方式的仿真模型，我们可以：
- 直观观察不同接地方式下的故障特征
- 验证继电保护算法的有效性
- 为实际系统保护配置提供理论依据
- 缩短现场试验周期和成本
2. 模型搭建关键步骤
2.1 系统参数设计
典型10kV配电网参数设置示例：
matlab复制% 系统基准值
Sb = 10e6; % 基准容量10MVA
Vb = 10e3; % 线电压10kV
Zb = Vb^2/Sb; % 基准阻抗
% 线路参数（π型等效）
R = 0.17; % Ω/km
X = 0.38; % Ω/km
C = 0.0076; % μF/km
线路长度 = 10; % km
% 变压器参数
变比 = 110/10.5; % 110kV/10.5kV
短路阻抗 = 0.1; % 10%
2.2 中性点不接地模型搭建
关键模块配置：
- 电源模块：使用Three-Phase Source，设置内阻和相角
- 线路模块：采用Distributed Parameter Line或π型等效电路
- 负荷模块：用Three-Phase Parallel RLC Load实现
- 故障模块：配置Three-Phase Fault实现单相接地
重要提示：不接地系统仿真时，需将电源中性点直接接地电阻设置为1e6Ω以上来模拟浮空状态。
2.3 消弧线圈接地模型实现
消弧线圈关键参数计算：
matlab复制% 系统对地电容电流计算
Ic = 3 * ω * C * 线路长度 * Vb/sqrt(3) * 1e-6;
% 消弧线圈补偿度取5%-10%
L = 1/(ω^2 * C_total * (1.05~1.10));
在Simulink中使用Three-Phase Series RLC Branch模块实现消弧线圈，接线方式为：
code复制电源中性点 → 消弧线圈 → 接地
3. 故障特征对比分析
3.1 中性点不接地系统特征
当A相发生金属性接地时：
- 故障相电压降为0
- 非故障相电压升高至线电压
- 零序电压U0 ≈ 相电压
- 故障点电流仅为容性电流（一般<30A）
仿真波形特征：
python复制# 伪代码展示波形特征
if 接地故障发生:
Ua = 0
Ub = Uc = 线电压
Ia = 容性电流(相位超前U0约90°)
3.2 消弧线圈接地系统特征
采用过补偿方式时：
- 故障点残流为感性电流
- 中性点位移电压降低
- 故障相电压不完全为零
- 接地点电流进一步减小（通常<10A）
关键参数关系：
code复制残流 = 电容电流 - 消弧线圈电流
补偿度 = (IL - Ic)/Ic × 100%
4. 保护算法实现与验证
4.1 零序电流保护实现
典型逻辑框图：
code复制零序电流采集 → 滤波处理 → 幅值比较 → 延时元件 → 出口动作
Simulink实现要点：
- 使用Three-Phase Sequence Analyzer提取零序分量
- 配置Moving Average滤波器消除高频噪声
- Relay模块设置动作阈值（通常0.2-0.5A）
4.2 五次谐波选线法
MATLAB函数示例：
matlab复制function [fault_phase] = fifth_harmonic_selection(I0)
[pxx,f] = periodogram(I0,[],[],采样频率);
五次谐波能量 = pxx(find(f==5*基频));
[~, fault_phase] = max(五次谐波能量);
end
5. 仿真技巧与问题排查
5.1 常见报错解决
-
代数环问题：
- 现象：仿真无法启动，提示"Algebraic loop"
- 解决：在适当位置添加Unit Delay模块
- 原理：打破信号直接反馈形成的代数环
-
波形失真：
- 检查步长设置：最大步长设为1/1000基波周期
- 验证求解器：推荐使用ode23tb或ode15s
5.2 提高仿真效率的方法
-
对不关注的部分使用简化模型：
- 用Current Source代替详细发电机模型
- 用集中参数线路代替分布参数模型
-
并行计算设置：
matlab复制parpool('local',4); % 开启4个worker
spmd
sim('model_name');
end
6. 模型扩展应用
6.1 弧光接地故障模拟
实现方法：
- 使用Controlled Voltage Source模拟动态电弧
- 电弧模型采用Cassie或Mayr方程：
code复制d(g)/dt = (1/τ)*(i·v/P - g)
其中g为电弧电导，τ为时间常数
6.2 与继电保护装置联合测试
硬件在环(HIL)实现步骤：
- 将保护装置通过IO板卡与仿真机连接
- 使用Simulink Real-Time生成实时代码
- 配置信号调理接口：
- 电压输出：0-10V对应0-1.2Un
- 电流输出：0-5A对应0-20In
7. 工程应用建议
-
参数整定原则：
- 零序电流保护：Iset > 1.5倍最大不平衡电流
- 动作延时：0.5-1.5s（躲过瞬时干扰）
-
现场试验对照：
- 仿真结果应与注入试验数据偏差<15%
- 重点关注相位关系一致性
-
新型保护算法验证：
- 暂态能量法
- 小波变换法
- 人工智能算法