1. 项目背景与核心价值
小电流接地系统是6-35kV配电网最常见的接地方式，主要包括中性点不接地和经消弧线圈接地两种形式。这类系统在发生单相接地故障时，故障电流较小（通常不超过10A），系统可继续运行1-2小时，但需要及时定位并排除故障点。实际电力系统中，直接进行故障试验成本高、风险大，因此Simulink仿真成为研究接地故障特性的重要手段。
这个仿真模型的价值在于：
- 可直观展示两种接地方式下的零序电流、电压波形差异
- 能模拟不同过渡电阻（0-10kΩ）下的故障特征
- 为继电保护算法开发提供测试平台
- 帮助理解消弧线圈的补偿原理
提示：仿真时建议采用Matlab R2020b及以上版本，其对电力系统模块库（SimPowerSystems）的支持更完善。
2. 模型搭建关键步骤
2.1 系统主电路构建
在Simulink中新建空白模型，从Simscape/Electrical/Specialized Power Systems库中拖拽以下元件：
- 三相电压源（参数：10.5kV，50Hz）
- 三相双绕组变压器（Δ/Y接线，10.5kV/0.4kV）
- 分布式参数线路模块（设置正序阻抗Z1=0.17+j0.38Ω/km）
- 负载模块（阻感性负载，功率因数0.85）
接地方式通过切换开关选择：
- 不接地系统：直接悬空中性点
- 消弧线圈接地：中性点经L=0.33H电感接地（对应10kV系统补偿电流10A）
2.2 故障模拟模块
使用三相故障模块（Three-Phase Fault）实现：
- 故障类型选择A-G（单相接地）
- 故障起始时间设为0.1s
- 过渡电阻可调（默认设置为100Ω）
- 通过接地开关控制是否投入消弧线圈
2.3 测量与显示配置
关键测量点包括：
- 三相电压、电流（使用V-I Measurement模块）
- 零序电压（3U0）：通过三相电压求和获得
- 零序电流（3I0）：中性点电流或三相电流矢量和
- 消弧线圈电流（仅补偿模式）
示波器建议分组显示：
- Scope1：三相电压/电流波形
- Scope2：零序分量对比
- Scope3：消弧线圈电流与故障点电流相位关系
3. 核心参数计算与设置
3.1 消弧线圈补偿度计算
补偿电流公式：
[ I_L = \frac{U_{ph}}{X_L} = \frac{10.5kV/\sqrt{3}}{2πfL} ]
以10A补偿电流为例：
[ L = \frac{10.5kV/\sqrt{3}}{2π×50×10} ≈ 0.33H ]
3.2 线路对地电容估算
电缆线路每相对地电容：
[ C_0 ≈ \frac{0.055ε_r}{\ln(D/r)} [μF/km] ]
其中D为线间距，r为导体半径，ε_r为绝缘介电常数
典型10kV电缆参数：
- 导体半径8mm
- 绝缘厚度4.5mm
- ε_r=2.3
计算得C0≈0.2μF/km
3.3 故障电阻设置原则
仿真时应覆盖典型场景：
- 金属性接地：0.1-1Ω
- 电弧接地：100-1000Ω
- 高阻接地：1-10kΩ
4. 仿真结果对比分析
4.1 中性点不接地系统特性
- 故障相电压降为0，非故障相电压升高至线电压
- 零序电压3U0=相电压（约6kV）
- 零序电流3I0=3ωC0Uph（容性电流）
- 故障点电流超前零序电压90°
典型波形特征：
code复制故障相电压 [V]
|‾‾‾‾‾‾‾‾‾‾
| |
|__________|
非故障相电压 [V]
__________
| |
| |‾‾‾‾‾‾‾‾
|__________|
4.2 消弧线圈接地系统特性
- 零序电压大小与补偿度相关（欠补偿时略低于相电压）
- 故障点电流=│Ic-IL│（残余电流通常<5A）
- 电流相位关系：
- 全补偿时：故障点电流最小
- 欠补偿时：残余电流呈容性
- 过补偿时：残余电流呈感性
注意：实际系统中通常采用过补偿5-10%，避免全补偿时的谐振风险。
5. 工程应用中的关键问题
5.1 故障选线方法验证
可通过仿真测试不同选线算法的有效性：
- 零序电流幅值比较法
- 五次谐波法
- 暂态能量法
- 群体比幅比相法
5.2 电弧建模技巧
使用非线性电阻模拟电弧：
matlab复制R_arc = R0 * exp(-t/τ) // τ为时间常数
或采用Cassie/Mayer电弧模型
5.3 仿真加速技巧
- 使用变步长求解器ode23tb
- 对线性网络启用"加速"选项
- 合理设置最大步长（建议≤1e-4s）
6. 模型扩展方向
-
多分支线路仿真：
- 添加3-5条出线
- 设置不同长度的电缆/架空线混合线路
- 验证故障定位算法
-
间歇性电弧接地：
- 使用可控开关模拟电弧重燃
- 观察过电压现象（可达3.5p.u.）
-
铁磁谐振模拟：
- 在PT一次侧添加非线性电感
- 激发分频谐振（1/3次谐波）
-
保护装置测试：
- 接入自定义保护逻辑（S-function）
- 测试动作准确性与速动性
7. 常见问题排查
7.1 仿真不收敛
- 检查是否存在拓扑冲突（如电压源短路）
- 尝试更改求解器为ode15s
- 减小初始步长（1e-6s起步）
7.2 波形异常
- 确认所有元件参数单位一致（kV vs V）
- 检查接地连接是否完整
- 验证变压器接线组别设置
7.3 零序测量误差
- 确保三相测量同步采样
- 使用Simulink的"Three-Phase Sequence Analyzer"模块
- 检查滤波器设置（通常需要低通滤波）
8. 实操经验分享
-
参数扫描技巧：
matlab复制for R_fault = [0.1 100 1000 10000] set_param('model/ThreePhaseFault', 'Ron', num2str(R_fault)); simout = sim('model'); % 分析数据... end -
波形标注最佳实践：
- 使用Simulink的"Signal Label"功能
- 添加图例说明曲线对应关系
- 导出数据到MATLAB workspace进行FFT分析
-
模型版本控制：
- 使用Git管理不同配置版本
- 通过Model Properties记录修改历史
- 对关键参数使用变量而非硬编码
-
性能优化记录：
- 某次仿真将步长从1e-5调整为5e-5后，运行时间从43s降至17s
- 启用"加速器"模式后，100次蒙特卡洛仿真耗时从2.1h降至38min
这个模型经过多次迭代，最初版本只能模拟金属性接地，现在已扩展支持：
- 可变过渡电阻（0-50kΩ）
- 动态补偿调节（自动调谐消弧线圈）
- 故障录波数据导出（兼容COMTRADE格式）
在实际教学中发现，学生最容易忽略的是：
- 忘记设置变压器的磁化曲线（导致铁磁谐振模拟不准确）
- 混淆线路的集中参数与分布参数模型适用场景
- 未考虑测量设备的带宽限制（特别是高频暂态分析时）