1. 项目背景与核心价值
小电流接地系统是6-35kV配电网最常见的运行方式，其中性点不接地或经消弧线圈接地的设计能有效限制单相接地故障电流。但这类系统故障特征隐蔽，传统保护装置易出现误判。通过Simulink搭建精确的仿真模型，我们可以获得：
- 不同接地方式下的零序电流/电压波形特征
- 故障暂态过程的电磁特性可视化呈现
- 保护算法验证的安全实验环境
某供电局2022年故障统计显示，小电流系统单相接地故障占比达83%，其中因保护误动/拒动导致的停电事故占37%。精确的仿真模型对提高保护可靠性具有直接工程价值。
2. 模型架构设计要点
2.1 系统主电路建模
采用三相电压源（10.5kV/50Hz）经π型线路等效供电，关键参数包括：
matlab复制R1 = 0.17 Ω/km; % 正序电阻
L1 = 1.2 mH/km; % 正序电感
C1 = 9.5 nF/km; % 正序电容
线路长度建议设置为5-10km以模拟典型配网场景，过长会导致电容电流过大超出小电流系统范围。
注意：电压源内阻需设为0.1%以下，否则会影响中性点位移电压精度
2.2 接地方式实现
中性点不接地模型：
直接保持变压器二次侧中性点悬空，此时系统对地电容将决定故障电流通路。
消弧线圈接地模型：
需使用Three-Phase Series RLC Branch模块，典型参数配置：
matlab复制L = 0.5 H; % 电感量（按过补偿5%设计）
R = 10 Ω; % 线圈等效电阻
建议通过Powergui的Load Flow计算功能校核补偿度是否合理。
2.3 故障注入设计
使用Three-Phase Fault模块实现：
- A相金属性接地（过渡电阻0Ω）
- 高阻接地（500-1000Ω）
- 间歇性电弧接地（通过Timer控制通断）
实操技巧：故障触发时间建议设为0.1-0.2s，避开系统启动暂态过程
3. 关键测量与信号处理
3.1 零序信号提取
采用三相电压/电流求和法：
matlab复制V0 = (Va + Vb + Vc)/3;
I0 = Ia + Ib + Ic;
需在测量环节添加Butterworth低通滤波器（截止频率150Hz）抑制高频噪声。
3.2 特征量计算
- 零序电压突变量：ΔU0(k) = |U0(k)| - |U0(k-1)|
- 五次谐波含量：通过FFT模块提取100Hz分量
- 暂态能量积分：∫(i0²(t) + u0²(t))dt
3.3 示波器布局建议
创建多视图显示界面：
- 顶层：三相电压/电流波形
- 中层：零序量时域波形
- 底层：谐波分析频谱图
4. 仿真参数配置要点
4.1 求解器选择
推荐采用ode23tb（刚性系统专用）：
matlab复制Max step size: 50e-6 s
Relative tolerance: 1e-4
Absolute tolerance: 1e-6
4.2 采样率设置
满足Nyquist定理前提下：
- 基波分析：10kHz足够
- 暂态过程研究：建议≥100kHz
- 电弧建模：需提高到1MHz
4.3 常见报错处理
代数环问题：
在电压测量端添加1e-6s延时模块
发散振荡：
检查线路参数单位是否统一（Ω/km还是Ω）
5. 典型仿真结果分析
5.1 金属性接地特征对比
| 特征量 | 不接地系统 | 消弧线圈接地 |
|---|---|---|
| 零序电压(V) | 6.06kV | 5.8kV |
| 零序电流(A) | 8.7 | 1.2 |
| 五次谐波占比% | 3.2 | 15.8 |
5.2 高阻接地难点
当过渡电阻>700Ω时：
- 不接地系统U0可能低于15%Un（保护启动阈值）
- 消弧线圈系统会出现欠补偿现象
解决方案：增加暂态量启动元件
6. 模型验证方法
6.1 理论校验
- 不接地系统：验证Ic = 3ωC0Uph
- 消弧线圈：检查IL/IC比值是否在0.95-1.05间
6.2 实物对比
将仿真波形与PT/CT录波数据对比，重点关注：
- 故障初相角对暂态过程的影响
- 电流过零点的相位偏移
7. 高级应用扩展
7.1 保护算法测试
在模型中接入自定义保护逻辑（通过MATLAB Function模块），测试：
- 基波方向元件灵敏度
- 谐波制动特性
- 小波变换暂态检测
7.2 弧光接地模拟
通过非线性电阻模块实现：
matlab复制R_arc = 1000/(1 + exp(-(t-0.1)*100)) + 50*rand();
7.3 参数优化设计
使用Simulink Design Optimization工具箱自动调整消弧线圈抽头位置，目标函数设为：
matlab复制min( |IL - IC| + 0.1*THD )
8. 工程经验实录
-
某110kV变电站10kV出线仿真发现：当电缆占比>70%时，消弧线圈补偿电流需增加12-15%
-
实测表明：间歇性接地故障的二次电弧建立时间通常在0.5-2ms范围内，仿真时需精确设置Timer参数
-
对于风电汇集线路，建议在模型中加入6nF/km的并联电抗器模拟电缆分布电容
-
发现Simulink 2021a版本存在零序测量模块相位偏移bug，可通过添加π/6相位补偿解决