1. 项目背景与核心价值
小电流接地系统是6-35kV配电网最常用的运行方式，其故障仿真对继电保护设计至关重要。中性点不接地和经消弧线圈接地这两种典型方式，在单相接地故障时呈现完全不同的暂态特性。通过Simulink搭建精确的仿真模型，可以直观观察到：
- 中性点不接地系统产生的高幅值暂态过电压
- 消弧线圈对接地电流的补偿效果
- 不同接地电阻下的故障特征差异
这种仿真能力对于保护装置定值校验、选线算法验证以及教学演示都具有不可替代的价值。笔者在电力系统仿真领域有8年实战经验，本文将分享从模型搭建到结果分析的全套方法论。
2. 模型架构设计要点
2.1 系统主电路建模
采用三相电压源（10.5kV/50Hz）通过π型等效线路连接负载，关键参数包括：
matlab复制R_line = 0.17Ω/km # 线路电阻
L_line = 1.24mH/km # 线路电感
C_line = 9.8nF/km # 对地电容
提示：线路长度建议设置为5-10km，过短会导致故障特征不明显，过长会增大计算量。
2.2 接地方式实现
- 中性点不接地：直接保持变压器次级绕组中性点悬空
- 消弧线圈接地：需在MATLAB Function模块中实现动态调谐算法：
matlab复制function L = auto_tuning(U0, I0)
% U0: 中性点位移电压
% I0: 接地电容电流
w = 2*pi*50;
L = 1/(3*w^2*C) * (1 + 0.1*sign(U0)); % 带5%过补偿
end
2.3 故障模拟模块
使用Three-Phase Fault模块实现：
- 故障类型：A相金属性接地/经电阻接地
- 触发方式：定时触发或手动控制
- 过渡电阻范围：0-10kΩ可调
3. 关键仿真技术实现
3.1 零序电流检测电路
在每条出线装设零序CT，通过Simulink的Current Measurement模块获取信号后，用以下处理链：
code复制三相电流 → 克拉科夫变换 → 0.5Hz高通滤波 → 有效值计算
3.2 暂态特征提取
配置Powergui模块进行FFT分析，重点关注：
- 故障初瞬的5次谐波含量
- 零序电压/电流相位差
- 暂态过程持续时间（通常10-20ms）
3.3 仿真步长选择
推荐采用变步长ode23t算法，最大步长设置为：
$$
\Delta t_{max} = \frac{1}{10f_{highest}} = 0.1ms
$$
对于需要观察电弧重燃的场景，建议局部加密至0.01ms。
4. 典型仿真结果对比
| 特征参数 | 中性点不接地 | 消弧线圈接地 |
|---|---|---|
| 零序电压倍数 | 1.73Un | <0.3Un |
| 接地电流 | 3-30A | <5A |
| 暂态振荡频率 | 300-1500Hz | 50-100Hz |
| 故障相电压恢复 | 无 | 快速恢复 |
5. 工程应用中的注意事项
-
模型验证：先用理论公式计算电容电流，与仿真结果偏差应<5%
$$
I_C = 3U_phωC_{total}
$$ -
参数敏感性：线路对地电容误差对结果影响最大，建议实测获取
-
实时性优化：
- 关闭所有Scope的"Limit data points"选项
- 将非关键模块设为刚性系统（Stiff）
-
教学演示技巧：
- 在故障点添加Spark Gap可视化效果
- 用XY Graph展示电压电流轨迹图
6. 常见问题解决方案
Q1：仿真出现代数环错误
- 检查所有测量模块与电源的电气连接
- 在适当位置插入Unit Delay模块
Q2：消弧线圈补偿效果不理想
- 调整调谐算法的过补偿度（建议8-12%）
- 检查线圈电阻参数是否过大
Q3：暂态过程采样不足
- 修改Powergui的Solver为Discrete
- 设置固定步长0.01ms
笔者在多个变电站改造项目中验证过该模型的准确性。特别提醒：实际系统中消弧线圈的调谐精度会显著影响仿真结果，建议配合PSCAD进行交叉验证。对于需要获取完整模型文件的研究者，可通过文末联系方式获取经过工业验证的仿真模板。