做继电保护的人绕不开小电流接地系统。本科《电力系统继电保护》课上会讲现场做选线装置调试的工程师更要天天面对它。原因很实在中性点不接地和经消弧线圈接地这两种方式占了10kV配电网的绝对主流而单相接地故障又是配电网里发生概率最高的故障类型。可是真要把这两种方式下的故障特征讲清楚、调明白光靠公式推导不够光靠现场试验更不行——现场想人为制造一次单相接地风险和成本都太高。所以用Simulink做故障仿真是最稳妥、也最能出效果的路子。这篇文章完整记录一下我做“小电流接地系统单相接地故障仿真”的整个过程包含中性点不接地和经消弧线圈接地两套模型的搭建、参数计算、波形分析和避坑经验适合电力专业学生、继保方向的研究生以及正在做配网选线项目的工程师参考。1. 项目背景小电流接地系统为什么值得搞明白1.1 小电流接地系统是怎么定义的电力系统的中性点接地方式本质上就两类大电流接地和小电流接地。直接接地、经低电阻接地故障电流很大属于大电流接地中性点不接地、经消弧线圈接地、经高电阻接地故障电流被限制到很小属于小电流接地。国外把这两种方式叫有效接地和非有效接地叫法不同内核一样。我国3~66kV配电网主要采用小电流接地方式其中10kV系统最常见的就是中性点不接地和经消弧线圈接地。选择这两种方式的核心考虑是供电可靠性单相接地时故障电流小电弧一般可以自熄系统还能带故障继续运行1~2小时给调度留出倒负荷、查故障的时间。代价就是选线和定位变难了这就是继电保护里最经典的矛盾——可靠性上去了选择性下来了。仿真模型里我按10kV系统来做因为这是现场最常见的电压等级而且作为学习案例10kV的参数比较好理解。真要换成6kV或35kV原理不变只要换电压等级和线路参数就行。1.2 单相接地故障占比高、难处理、必须深入理解配电网故障里单相接地占了80%以上剩下才是两相短路、三相短路这些。两相短路、三相短路特征明显电流激增、电压骤降保护容易动作单相接地恰恰相反故障电流只有电容电流量级几安到几十安电压变化又不至于让负荷停摆所以保护装置很难直观捕捉。更麻烦的是选线问题。中性点不接地系统里非故障线路也有零序电流故障线路的零序电流反而更大方向也相反这还算有规律可循。一旦中性点经消弧线圈接地故障点残流被补偿得很小故障线路和非故障线路的零序电流差别变得非常微弱传统稳态选线方法直接失效。现场很多选线装置误报、漏报根子就在这里。做这个仿真目的就是把这两种接地方式下的故障特征量化、可视化从波形和数值上理解“为什么选线难”“保护要怎么整定”再去研究解决方案就有抓手了。1.3 这次仿真要实现什么目标这个项目的核心目标有三个。第一搭建一个可复现的10kV小电流接地系统Simulink仿真模型包含电源、变压器、多条出线、负载和故障模块能模拟A相单相接地金属性故障。第二实现中性点接地方式的灵活切换——不接地和经消弧线圈接地两套方案在同一个模型框架内完成便于对比。第三提取母线零序电压、各出线零序电流、消弧线圈电流等关键电气量分析故障前后波形特征和相量关系为理解零序保护原理和选线方法提供直观依据。这三个目标做完读者不仅能得到能跑的模型还能把“为什么这么搭、为什么这样算参数”搞清楚。下面按我的实际搭建顺序来写。2. 两种中性点接地方式的故障特征与方案选型2.1 中性点不接地系统的故障特征中性点不接地系统在正常运行时三相对地电压对称中性点电位为零。A相发生金属性接地后A相电压被强制拉低到接近零中性点发生位移B相和C相对地电压升高到线电压水平。如果用有效值说话正常时相电压是5.77kV故障后非故障相电压变成10kV也就是升高了根号3倍。这个电压变化是所有小电流接地系统共有的特征也是绝缘监测装置能发出接地告警的物理基础。零序电压3U0在故障后会由接近零变成约等于相电压这个量可以从母线电压互感器的开口三角绕组直接取到做保护的启动判据非常方便。电流侧的规律更关键。非故障线路的零序电流本质上就是本线路的对地电容电流方向从母线流向线路数值等于3ωC₀Uφ。故障线路的零序电流等于整个系统中除本线路外所有线路对地电容电流之和方向从线路流向母线。也就是说故障线路的零序电流最大且方向与正常线路相反。这就是稳态零序电流幅值法和方向法的原理基础。2.2 经消弧线圈接地系统是如何工作的中性点不接地系统最大的隐患是间歇性电弧接地时容易产生过电压甚至发展成相间短路。消弧线圈的作用就是利用中性点接入的电感电流去抵消单相接地时流过故障点的容性电流让电弧电流减小到能自熄的程度。补偿方式有三种欠补偿电感电流小于电容电流、全补偿两者相等、过补偿电感电流大于电容电流。全补偿会引发串联谐振中性点位移电压很高现场绝对禁止。欠补偿在系统运行方式变化时可能失去补偿作用也不可靠。所以工程上统一采用过补偿脱谐度一般在5%~10%也就是电感电流比电容电流大5%~10%。消弧线圈投入后故障点残流等于电容电流和电感电流的差值方向由两者大小决定。过补偿状态下故障点残流呈感性故障线路零序电流方向会反转这正是消弧线圈系统选线困难的直接原因。理解这一点仿真时看波形方向就自然明白了。2.3 仿真建模思路与系统参数选型模型整体思路用一个110kV等效电源加一台Yd11变压器把电压降到10kV10kV侧引出4条出线故障点设在其中一条出线上。变压器低压侧Y中性点引出后串联一个可投切开关再接消弧线圈——开关闭合就是经消弧线圈接地开关断开就是中性点不接地。这样一个模型两套系统切换方便对比直观。为什么不直接用一个10kV三相电源因为Three-Phase Source模块本身不提供中性点引出接不了消弧线圈。实际变电站里10kV系统的中性点也是通过主变低压侧Y绕组引出的。仿真要贴近实际就得把变压器加进来。线路选型上用Three-Phase PI Section Line来模拟架空线和电缆的等值参数。为了让单相接地电容电流落在典型范围本模型按30A左右设计在10kV母线上并联了一组集中等值对地电容用来等效较长电缆网络的对地电容效应。这个做法在工程仿真里很常见避免为了加大电容电流而无限制加长线路导致模型跑起来又慢又笨。3. Simulink建模实操从库选型到参数设置3.1 模型顶层架构和模块清单我的模型搭建顺序是先放电源和变压器再引出母线和出线接着加负载和故障模块最后加测量环节。全部使用Simscape Electrical库里的专用电力系统模块不要用基础Simulink库搭电路那样既不准确调试也麻烦。核心模块清单如下模块库内名称作用110kV电源Three-Phase Source提供系统电源主变Three-Phase Transformer110kV/10kV降压Yd11接线输电线路Three-Phase PI Section Line模拟10kV出线等值参数集中对地电容Three-Phase Parallel RLC Branch等效系统对地电容故障模块Three-Phase Fault模拟A相单相接地消弧线圈Series RLC Branch补偿电容电流测量模块Three-Phase V-I Measurement测量电压电流零序提取Sequence Analyzer提取零序电压、零序电流显示Scope观察波形模块数量不多但接线需要仔细。高压侧电源接变压器高压侧三角形绕组变压器低压侧星形绕组引出10kV母线中性点经断路器和消弧线圈接地。四条出线从母线上引出每条出线末端接三相负载其中一条出线上串入故障模块。3.2 电源与变压器配置110kV电源的参数我设置如下线电压有效值110kV频率50Hz三相短路容量取1000MVAX/R比取10。短路容量大说明系统比较强故障前后母线电压基本由主变和线路决定仿真结果更干净。变压器选择Three-Phase Transformer接线方式Yd11额定容量50MVA。这里把高压侧接成Y型、低压侧接成Δ型还是反过来有讲究对Yd11变压器来说低压侧是三角形高压侧是星形而我们需要从低压10kV侧引出中性点所以应该让10kV侧做成星形接线。标准Yd11变压器一般是高压侧星形、低压侧三角形所以这里反过来用选Yd11时要注意把两侧电压等级设对10kV侧选星形Y110kV侧选三角形D这在实际工程中对应的是“高压侧三角形、低压侧星形”的变压器只是为了引出中性点的合理设定。正常运行时10kV侧相电压为5774V线电压10kV。变压器短路阻抗设为10.5%负载损耗和空载损耗按默认值即可它们对故障特征影响不大。3.3 输电线路与对地电容参数设置线路参数直接影响零序电流大小是全模型最需要较真的地方。Three-Phase PI Section Line模块里需要填正序电阻、正序电感、正序电容以及零序电阻、零序电感、零序电容单位分别是Ω/km、H/km、F/km还要设置线路长度。我的4条出线设定如下三条10kV架空线每条长度15km正序电阻0.27Ω/km正序电感1.3mH/km正序电容0.0096μF/km零序电阻0.35Ω/km零序电感3.2mH/km零序电容0.006μF/km。第四条是电缆线路长度5km正序电阻0.12Ω/km正序电感0.31mH/km正序电容0.4μF/km零序电阻0.35Ω/km零序电感1.2mH/km零序电容0.3μF/km。单靠这些线路的分布电容系统对地电容电流并不大主要来自电缆线路。为了把系统电容电流撑到典型工业现场的量级模型里在母线上并联了一组集中对地电容等效零序电容约5μF。这样系统单相接地电容电流大致为I_c 3 × ω × C₀Σ × Uφ 3 × 314 × 5.5×10⁻⁶ × 5774 ≈ 29.9A这个30A左右的电容电流在10kV电缆网络的现场是常态。消弧线圈参数也基于这个值来设计逻辑才成立。3.4 故障模块与消弧线圈参数计算故障模块用Three-Phase Fault这是模拟单相接地最直接的组件。需要勾选Phase A和Ground表示A相经故障点接地其他相不要勾选。Transitional resistance设为0.001Ω模拟金属性接地。Transition Time设成[0.2 0.4]意思是0.2秒故障发生0.4秒故障切除仿真时长0.6秒这样能清楚看到故障前、故障中、故障后的三段波形。消弧线圈参数需要实际算一遍不能瞎填。按过补偿10%设计补偿后的电感电流IL取33A系统相电压Uφ等于5774V那么消弧线圈的感抗为X_L Uφ / IL 5774 / 33 ≈ 175Ω电感量就是L X_L / ω 175 / 314 ≈ 0.557H消弧线圈并非理想电感内部有绕组电阻和铁芯损耗等值电阻一般取感抗的3%~10%。取R6Ω这样品质因数比较高仿真中不会因为电阻过大影响补偿效果也不会因为电阻过小产生数值振荡。在Series RLC Branch里填Resistance6ΩInductance0.557H串联接在中性点和地之间。这里有个关键点消弧线圈支路必须串接在变压器中性点和地之间而不是随便并在母线上。中性点不接地时把连接消弧线圈的断路器断开即可经消弧线圈接地时断路器合上。用Breaker模块加Step信号控制信号为1合闸0分闸。3.5 测量环节与波形采集测量环节是理解仿真结果的关键。在母线处和每条出线首端分别放置Three-Phase V-I Measurement模块电压测量用“Phase-to-ground”方式电流测量选三相电流输出类型用Complex。零序电压3U0的获取有两个办法。一是用Sequence Analyzer模块直接把三相电压分解成正序、负序、零序分量输出零序电压二是用加法器把三相电压瞬时值加在一起再乘以1/3就是零序电压瞬时值。我推荐Sequence Analyzer它输出的是相量幅值和相角观察相位关系非常方便。零序电流3I0同理在每条出线首端的电流测量后接Sequence Analyzer得到各出线的零序电流相量。这里要特别注意测量点是出线首端、靠近母线侧这样才能正确反映流向量故障线路和非故障线路的差别。从波形上可以直接对比零序电流的幅值和相位配合故障线路识别方法一起验证。示波器Scope建议分组一组看母线三相电压一组看3U0一组看4条出线的3I0一组看消弧线圈电流。这样每个物理量单独显示波形清楚不互相干扰。如果想导出数据做进一步分析可以在Scope里点右上角齿轮把“Log data to workspace”勾上仿真后通过base workspace直接拿到变量。4. 仿真结果分析与继电保护应用4.1 中性点不接地系统的故障波形解读跑完不接地模型最直观的现象是母线三相电压的变化。故障前三相电压对称幅值相等0.2秒A相接地后A相电压跌落趋近于零B相和C相电压升高。用有效值看B、C相电压从原来的5.77kV升到10kV左右正好是线电压水平。这个波形和理论完全对应说明模型搭得对。零序电压3U0的变化更关键。正常时近似为零故障后零序电压有效值跳到接近相电压5.77kV。这个量是绝缘监测或零序保护启动的判据仿真里看到3U0稳定出现就知道接地信号是可靠的。再看各出线零序电流的对比。非故障线路的3I0约等于本线路的电容电流只有几安培故障线路的3I0则要大得多接近整个系统的电容电流30A左右。更重要的是方向非故障线路零序电流超前零序电压90度是容性故障线路零序电流滞后零序电压90度看起来两者方向相反。用相量图观察时3U0与3I0的相位关系一览无余这就是零序方向保护的核心依据。4.2 经消弧线圈接地系统的补偿效果把中性点断路器合上消弧线圈接入后再次仿真波形会发生明显变化。最显著的是故障点残流大幅下降。这是因为消弧线圈产生的电感电流抵消了大部分电容电流系统电容电流30A左右消弧线圈补偿电流33A左右残流只有大约3A。从零序电流波形看故障线路的3I0幅值明显减小不再像不接地时那样“一枝独秀”。非故障线路的零序电流依然是各自电容电流故障线路零序电流在过补偿状态下呈感性方向发生了反转。这就是现场选线装置在消弧线圈系统里经常失灵的原因——特征量被削弱方向还变了传统的幅值比较和方向比较都不可靠。从保护的角度讲消弧线圈带来的“副作用”是选线难但换个角度看它也有收益故障点电流小了电弧自熄概率提高系统可以继续带故障运行供电可靠性上升。到底选不选消弧线圈本质是可靠性和选择性的权衡仿真结果把这种权衡摆在了眼前。4.3 从仿真波形看选线保护的整定思路看完波形再谈整定就有依据了。中性点不接地系统里零序电流幅值选线法可以奏效设定一个启动值大于该值的线路判定为故障线路。启动值要躲过非故障线路的最大电容电流又要低于故障线路的最小零序电流整定时需要根据系统实际出线长度和电容分布计算不能拍脑袋。经消弧线圈接地系统里稳态零序电流幅值法基本失效这时候工程上通常用无功功率方向法、五次谐波法、暂态量法、注入法这些更复杂的手段。比如五次谐波法消弧线圈对五次谐波的感抗会变大补偿效果变弱故障线路五次谐波电流特征反而更明显。仿真模型可以直接把这些方法做成算法模块接入测量结果进行验证。我在模型里就预留了接线端口把三条非故障线路和一条故障线路的3I0输出到同一个Scope里后面接一个简单的比较模块就能验证幅值法的局限。这也是初学者最容易上手的研究切入点。5. 实操中的常见问题与避坑指南5.1 仿真不收敛步长与求解器设置Simulink电力电子仿真最容易翻车的地方就是求解器设置。用默认变步长跑我这个模型大概率会出现波形毛刺甚至“Simulink cannot solve the algebraic loop”的报错。原因在于电力系统模型里有大量的快速暂态过程默认步长跟不上。我的做法是solver选择ode23tb最大步长设为1e-5秒相对容差改到1e-4。这样跑0.6秒仿真时间要计算6万步看起来慢但模型规模不大实际也就几十秒的事。如果发现波形还是振铃就把最大步长再压到5e-6。另外模型里的Three-Phase Fault模块自带过渡电阻参数如果设成0可能造成数值刚性建议至少给0.001Ω。消弧线圈支路的阻尼电阻也别设成0否则电感电流会出现不衰减的振荡波形没法看。5.2 故障模块和断路器为什么“没反应”这是初学者问得最多的问题故障模块放上去了时间参数也设了但波形一点变化都没有。排查顺序就三件事。第一确认Transition Time设置的时间在仿真时长范围内。比如仿真时长0.6秒故障时间设0.2到0.4秒中间时间段要覆盖住别设成0.8秒开始那跑到波形结束都看不到故障。第二确认故障相勾选正确。Three-Phase Fault默认可能只选了Phase A没选Ground单相接地要同时勾选对应相和Ground少一个等于没接地。第三确认模块的输出端已经接入被测线路而且故障模块不是悬空的——这点容易忽略模块画上去了却忘记串进主电路电气上相当于没接。断路器模块类似。带Breaker的模型控制信号必须是标量0/1而且要单独用Step模块给定不能直接用三相电压信号去控。另外Breaker内部有导通电阻默认值很大可能导致流过电流很小需要把导通电阻调到0.01Ω以下。5.3 消弧线圈参数错误导致的谐振问题消弧线圈参数如果按全补偿来设计仿真里可能能看到灾难性的一幕中性点位移电压异常升高三相电压严重不平衡零序电压比正常相电压还大。这就是串联谐振的表现。实际工程中绝对禁止全补偿运行仿真里也要避开。我之前用L0.524H、对电容电流30A的补偿刚好接近全补偿跑出来波形一度让我怀疑模型错了后来重新核算脱谐度才发现问题。脱谐度的计算很简单v (I_C - I_L) / I_C按过补偿设计I_L取33A、I_C取30A脱谐度就是(30-33)/30 -10%负号代表过补偿。符合这个范围波形才能物理正确。5.4 模型扩展还能往哪些方向玩两套基础模型跑通后后续扩展空间很大。我列几个自己试过或正在做的方向给读者参考。第一个方向是间歇性电弧接地仿真。用两个三相断路器串联让其中一个按特定频率反复通断模拟电弧重燃和熄灭过程。配电网单相接地里最危险的就是这种间歇性电弧它能产生很高的过电压用静态故障模块模拟不了必须上这种动态控制方案。第二个方向是选线算法验证。在模型输出端接数字信号处理模块尝试五次谐波法、首半波法、小波变换法对比每种方法在两套系统中的可靠性。这个方向适合做毕业设计或者科技论文内容很丰富。第三个方向是把仿真模型和代码生成、硬件在环打通。Simulink模型本身可以做C代码生成配合嵌入式控制器验证选线算法的实时性也可以用外部模式在Simulink和真实采集卡之间在线调试。这个模型作为算法测试床非常合适。第四个方向是对比经高电阻接地。调节消弧线圈支路为纯电阻观察零序电流特征的变化分析高电阻接地系统的选线差异性这在部分新能源并网场景里值得研究。我个人在实际操作中最大的体会是这个仿真项目的价值不在“能跑”而在于把抽象的原理变成了看得见、摸得着的信号关系。第一次看到非故障线路和故障线路零序电流方向相反时我对选线原理的理解比看十遍教科书都深刻。后面调试选线算法、分析录波数据、跟现场工程师讨论故障案例脑子里都是有画面的。最后再分享一个小技巧把消弧线圈支路和中性点断路器封装成子系统用手动开关切换接地方式整个模型就能一键对比两种系统。平时做仿真时我还会把Scope的波形导出成图片存档积累各种故障条件下的波形库。仿真做多了对故障特征的敏感度自然会上去这比死记硬背任何结论都管用。