做配电网继电保护研究这段时间我花了不少精力在小电流接地系统的仿真上。这类课题在电力系统里非常典型6kV到35kV的配电网中性点大多采用不接地或经消弧线圈接地的方式一旦发生单相接地故障电流小、特征弱保护装置很难准确判别所以很多论文和工程验证都绕不开Simulink仿真。我先后搭过中性点不接地和经消弧线圈接地两套故障模型从电压源参数设定、线路模型选择到故障模块配置踩了不少坑也积累了一些可以直接复用的实操经验。这篇就围绕这个仿真模型把原理、建模、参数算例和典型问题一次讲透。先说结论如果你正在做小电流接地系统故障仿真或者论文里需要对比不同接地方式下的故障特征用Simulink的Simscape Electrical库老版本叫SimPowerSystems完全够用。关键是把电容电流路径理清楚把消弧线圈的补偿度算对模型不会太复杂但每一步的物理概念必须正确。1. 为什么小电流接地系统的仿真要单独建模型1.1 小电流接地系统的定义与应用背景小电流接地系统学术上叫中性点非有效接地系统包括中性点不接地、经消弧线圈接地和经高电阻接地三种典型方式。它主要应用在66kV及以下的配电网。这个“小电流”是有对比对象的在中性点直接接地的大电流接地系统里单相接地故障时的短路电流可以达到几千安到几十千安但在小电流接地系统里单相接地时流过故障点的电流只是系统对地电容电流一般只有几安到几十安故障特征非常微弱。我在搭建仿真模型之前最先确认的就是系统所在的电压等级和线路规模因为这直接决定了对地电容电流的大小。实际工程中一条10kV出线的对地电容电流大概在几安到十几安量级整个变电站所有出线加起来可能几十安。这个数值直接影响消弧线圈的容量选择也决定了故障后非故障相电压的升高程度所以不是随便拍脑袋定的。很多初接触这个方向的人会问直接接地系统故障电流大保护不是更好动作吗为什么配电网非要搞小电流接地最核心的原因是供电可靠性。小电流接地系统发生单相接地故障时三相之间的线电压仍然保持对称所以允许继续带故障运行1到2小时运行人员有充足时间去查找故障点而不是立即跳闸停电。这种“带病运行”能力对末端用户来说很宝贵代价就是故障检测难度大。1.2 单相接地故障的核心特征与仿真目标单相接地最常见的是A相接地是小电流接地系统中最主要的故障类型占比超过80%。故障发生后系统的电气量变化有几个非常鲜明的特征故障相电压变成零非故障相电压升高到线电压即相电压的根号3倍母线零序电压3U0也从正常时的很小值跃升到接近相电压的值。线电压保持不变这是系统能继续运行的根本原因。我在仿真模型里要实现的核心目标就是准确复现这几个特征。具体来说要能清楚地看到故障相电压、非故障相电压、零序电压和故障点电流的波形变化同时要能在两种接地方式下对比同一个故障场景分析消弧线圈对故障电流和电压恢复过程的抑制作用。这两个目标实现之后后续做故障选线、保护判据研究才有基础。这个仿真适合谁来参考一种是电力系统相关专业的研究生论文方向涉及配电网单相接地故障检测或消弧线圈控制策略另一种是做配网自动化或继电保护整定的工程师需要用仿真数据辅助分析现场波形。如果你只是初学者第一次接触这类模型这篇文章里的参数计算方法和建模步骤也可以直接照搬不用自己去啃一堆参考书。2. 两种接地方式的故障机理拆解2.1 中性点不接地时电容电流与电压抬升过程中性点不接地系统在正常运行时三相对地电容基本对称所以中性点电位接近零系统各相对地电压就是相电压。但一旦A相发生金属性接地情况就完全不同了。A相的对地电容被短路A相电压直接变成零而B相和C相的对地电容两端电压从原来的相电压升高到线电压对地电容电流也相应增大。这里最需要理解的是故障点电流的路径。A相接地后电流从B相、C相的电源端流出经过线路对地电容回到大地再经过故障点流入A相最后回到电源。所以故障点电流等于全系统B相和C相对地电容电流的向量和而不是某一条线路单独贡献的。实际就是所有非故障线路的对地电容电流叠加到故障线路上这就是为什么故障线路的零序电流反而比非故障线路大而且是方向相反的关系。我建模型时特别验证了这一点单一出线模型里故障点电流等于该线路非故障相电容电流的根号3倍。如果接多回出线故障电流就是各回线电容电流的向量和所以系统规模越大、出线越多、电缆线路越长单相接地时的电容电流就越大。这也是设计中需要考虑的当系统电容电流超过一定阈值中性点不接地方式就不能满足电弧自灭的要求了需要升级为经消弧线圈接地。电压方面中性点不接地系统单相接地后会面临一个很实际的运行问题非故障相的对地绝缘长期承受线电压如果这个状态持续时间过长可能引发第二点接地形成相间短路。所以规程允许带故障运行1到2小时但运行人员必须尽快处理。仿真波形里非故障相电压稳定在1.732倍额定值这是判断单相接地故障的一个关键特征。2.2 消弧线圈的补偿逻辑与脱谐度整定消弧线圈本质上是一个带铁芯的电感线圈接在变压器中性点与大地之间。配电网中性点是通过接地变压器引出的这也是一个容易忽略的工程细节——普通的配电变压器一次侧没有中性点引出需要单独加接地变压器来构造中性点仿真模型里可以直接用接地的电压源模拟这个效果。消弧线圈的工作原理其实用一句话就能说明白单相接地发生后故障点流过系统对地电容电流消弧线圈提供了电感电流它恰好与电容电流方向相反两者在故障点相互抵消把故障点电流降到很小电弧就容易自己熄灭。这里的关键是补偿度的整定。补偿度K (IL - Ic) / Ic当IL大于Ic时叫过补偿小于Ic叫欠补偿相等叫全补偿。实际工程中几乎都采用过补偿也就是消弧线圈的电感电流稍大于电容电流脱谐度1-K通常控制在5%到10%。为什么不用全补偿因为全补偿时系统处于串联谐振状态正常运行时中性点会因三相不对称而产生很高的位移电压极容易引发过电压。欠补偿在系统部分线路退出运行后可能变为全补偿同样危险。所以只有过补偿才能确保任何运行方式下都不进入谐振状态。这一点在做仿真参数选择的时候如果直接按文献中的经验公式计算很容易忽略背后的逻辑。我仿真时采用的脱谐度是8%也就是IL取1.08倍Ic。计算电感值时先算系统单相接地时的电容电流再根据消弧线圈端电压即相电压算所需电感量。10kV系统相电压约5.77kV如果电容电流是14A那么消弧线圈的感抗需要满足IL Uφ / (ωL)算出来L大约0.237H。消弧线圈的损耗通过串联一个小电阻模拟Q值取10左右仿真波形里的衰减特性会更接近实际。2.3 两种接地方式在保护检测上的本质差异把两种方式放在一起对比保护检测的难点就清晰了。中性点不接地系统单相接地时零序电流的分布规律是故障线路零序电流大且方向与非故障线路相反这是方向型选线保护的基础。但因为故障电流只有几安培CT电流互感器误差、信号干扰都会对判别带来很大影响加上电缆线路和架空线路的电容电流差别很大传统基于零序电流幅值比较的选线方法并不可靠。经消弧线圈接地后情况更复杂。故障点流过的是补偿后的残流通常只有几安培甚至更小零序电流的幅值特征基本消失故障线路与非故障线路的零序电流方向也不再具有清晰的区分度。这也是消弧线圈系统里“故障选线难”的根源。不少研究人员会通过改变消弧线圈的脱谐度或注入特殊频率信号来实现选线仿真里验证这些策略时就需要把消弧线圈模型做得足够精细包括铁芯饱和特性、补偿度变化对零序电流的影响等。从保护角度来看两种系统最基本的共同判据是零序电压越限即3U0超过整定值通常取额定相电压的15%左右就启动报警或选线程序。仿真模型里的零序电压可以直接从三相电压测量模块中提取零序分量也可以自己搭建3U0计算回路两种方式我都试过后者在调试波形时更灵活。3. Simulink建模仿真从参数计算到模型搭建3.1 系统参数与算例设计仿真模型的电气参数是我在实际搭建之前就计算好的建议你也先做这一步而不是在模型里边改边试。算例的系统参数如下系统电压等级10kV中性点不接地/经消弧线圈接地电源侧三相理想电压源线电压10kV频率50Hz变压器10kV/10kV隔离变压器用于构造中性点引出容量5MVA线路两条10kV出线长度分别为5km和10km用π型集中参数模型故障设置单相A相接地故障时间0.1s后触发持续0.1s消弧线圈按过补偿8%整定串联电阻模拟损耗线路参数也是关键。我采用的是典型10kV架空线路参数正序电阻0.27Ω/km正序电感1.02mH/km正序电容9.6nF/km零序电阻2.7Ω/km零序电感4.2mH/km零序电容6nF/km。如果你手里有具体的电缆参数可以直接替换但要注意电缆的对地电容通常比架空线大一个量级电容电流会大很多故障特征也会更强。两条出线的总对地电容电流计算10km线路的零序电容电流按每公里0.07A估算5km线路就减半。这样总电容电流Ic约10.5A。如果使用电缆线路单公里电容电流可以达到1.4A以上差异非常明显。我做的是架空线路算例消弧线圈按Ic的1.08倍补偿线圈电流IL约11.3A。这些参数都会直接写入模型里。有一点要提醒Simulink模型的初值是影响仿真稳定性的重要因素最好让系统先运行0.1s无故障等电气量稳定后再投入故障模块这样故障波形的起振时刻才是准确的也不会出现开头一大段暂态振荡干扰分析。3.2 模型搭建步骤详解打开Simulink后我在Simscape Electrical库中按下面这个顺序搭建系统整体搭下来也就20分钟。先放一个Three-Phase Source作为电源线电压设为10kV内阻取默认值。然后是Three-Phase Transformer Two Windings用于把10kV系统降到0.4kV供负荷使用同时在其Yg/Δ接线的Y侧引出中性点——这一步直接对应实际变电站的接地变功能。线路模型方面我选的是Three-Phase PI Section Line。这个模块自带正序和零序参数输入框直接把前面计算的R、L、C填进去就行。每公里一条线还是每5公里一条线完全取决于你要观察的暂态细节。我仿真时把10km线路分成2段串联5km线路用1段这个做法在保证精度的同时不会让模型变得太复杂运行速度也快。故障模块用的是Three-Phase Fault它有A、B、C三个相和一个接地端子可以分别配置各相的故障类型。我把故障相选为A相故障类型选为“A相接地短路”过渡电阻设置为0.01Ω模拟金属性接地。要特别留意的是Transition Times的设置它有两个时间点第一个是故障投入时间第二个是故障切除时间。我设置为[0.1 0.2]即0.1s发生故障0.2s故障切除持续0.1s。如果设置反了模型会报错或者表现完全不符合预期。测量部分我用了两组Three-Phase V-I Measurement一组接在线路末端一组接在电源出口用来录三相电压和电流。故障发生后查看非故障相电压的波形分析零序电压变化都靠这两组测量数据。零序电压的计算可以直接在测量模块后接一个Zero-Order Hold再进FFT Analysis或者更简单地用数学运算模块把三相电压相加后除以3得到的就是3U0的标幺值这个信号直接反映了接地故障的严重程度。整个模型建立后需要在solver设置中先把仿真算法和步长选对。我建议选择ode23tb默认适用于电力电子和电力系统模型或ode15s最大步长设置为1e-4秒即0.1ms这样能捕捉到故障开关动作时的快速暂态过程。采样时间设得太大会出现波形锯齿太小的步长会拖慢仿真速度0.1ms对这类系统已经是比较合适的平衡点。3.3 消弧线圈建模与中性点引出细节消弧线圈在Simulink里实现起来不难一个电感串联一个电阻就够了但有几个细节容易忽略。我先添加了一个Series RLC Branch模块设成“串联R、L”模式R取0.1Ω模拟线圈铜损L按计算结果填入0.237H。然后把这个支路接在变压器中性点与大地的Ground之间。关键点在于中性点引出方式。我在模型里用了一个Yg/Δ接线的变压器Y侧中性点直接引出来接消弧线圈。这是仿真模型与真实工程高度一致的地方。如果你手头的模型里变压器只是简单的Δ/Y接线没有中性点引出就需要额外加一个Zigzag Transformer或Grounding Transformer模块来构造中性点。这也是很多初学者容易卡壳的地方——明明加了消弧线圈却发现系统还是找不到中性点故障电流路径建立不起来。消弧线圈投退的仿真可以这样设计在消弧线圈支路里串联一个理想开关用Breaker模块实现然后在不同的仿真中控制这个开关的开合状态实现“不接地”和“经消弧线圈接地”两种工况的切换。我在实际建模时是直接复制了一套系统模型一套不带消弧线圈一套带分开运行可以避免频繁改参数。还要提一下消弧线圈的过补偿校验。我根据上文的电容电流和脱谐度计算出IL11.3A然后在仿真里测量了故障点电流的实际波形残流应在1A以内证明补偿是有效的。如果残流明显偏大说明脱谐度没设对或者电感量算错了。排查方法很简单在故障点前面串联一个电流表看故障期间电流有效值再和理论计算值对比。4. 仿真结果解读两种工况的波形对比与保护意义4.1 中性点不接地系统故障波形特征先看不接地系统的仿真结果。故障发生前三相电压是对称的中性点电压基本为零。大约在0.1s时刻A相电压突降到零B相和C相电压从原来的5.77kV相电压峰值约8.16kV同步升高到10kV线电压水平峰值大约14.1kV。这个升高过程有一个短暂的过渡振荡因为线路电容和电感在故障瞬间发生了能量再分配之后很快进入稳态。零序电压的波形最值得看。正常运行时由于三相负载和线路参数的不完全对称零序电压有一个很小的底值通常不超过额定相电压的3%到5%。故障之后零序电压迅速上升到接近额定相电压波形基本是50Hz正弦波。这个特征就是“零序电压越限启动”的判据来源——我设置的阈值是额定相电压的15%即约866V模型运行到0.1s之后该判据立即满足。故障电流方面A相接地点的电流等于B相和C相对地电容电流的向量和仿真测得约为11.2A和理论估算的10.5A基本吻合。这个电流幅值不大但相位超前于零序电压约90度体现出容性电流的特征。故障线路首端的零序电流比非故障线路大相位差接近180度这条规律在仿真波形中非常明显也是后续做故障选线的基础。如果你用FFT分析工具看故障电流的频谱会发现除了基波50Hz之外还存在一定的高次谐波分量特别是当过渡电阻不是零时故障电流的波形会轻微畸变。实际工程中的小电流接地故障常伴随间歇性电弧仿真里用固定过渡电阻只能模拟稳定接地如果要做弧光接地研究就需要用Arc模块或者自定义状态方程这个复杂度更高但是是后续的一个自然扩展方向。4.2 经消弧线圈系统的补偿效果验证切换到带消弧线圈的模型后同样在0.1s设置A相接地故障点电流的波形发生了明显变化。由于消弧线圈提供了感性电流叠加在容性电流上两者的向量差就是残流。仿真结果显示残流有效值约1.0A到1.5A比不接地系统下降了约90%。这验证了消弧线圈的补偿效果——电弧能够自然熄灭。电压波形方面消弧线圈系统在故障期间非故障相电压同样升高到线电压水平这与不接地系统一致。但差别体现在故障切除后的过渡过程不接地系统故障切除后电压很快恢复稳定消弧线圈系统则会多出一个低频衰减振荡因为消弧线圈储存的磁能需要时间来释放。这个暂态过程如果谐振频率接近工频可能会出现短时的过电压这是设计消弧线圈时必须考虑的。零序电压在消弧线圈系统里出现了明显的相位变化。不接地系统的零序电流超前零序电压近90度经消弧线圈后零序电流滞后于零序电压一个角度这个角度的正切值和脱谐度相关。换句话说通过测量零序电流的相位可以判断系统当前的补偿状态很多消弧线圈自动跟踪补偿装置就是利用这个原理调节绕组分接头的。补偿度改变对波形的影响也值得一看。我把脱谐度分别调成3%、8%和15%做过对比脱谐度越小残流越小但故障后的过渡振荡越持久过电压幅值也越高。8%在工程上是比较折中的选择既能有效灭弧又不会让系统在故障切除后陷入长时间的不稳定状态。4.3 保护启动判据与选线逻辑的仿真验证思路模型的最终价值要落在保护策略验证上。我仿真时重点验证了三套逻辑零序电压越限启动、零序电流方向比较、零序电流幅值比较。零序电压越限是最简单可靠的一层任何单相接地故障必伴随这个特征但它的缺点是只能告诉你“发生了故障”不能告诉你“故障在哪条线路上”。零序电流方向比较则稍微复杂。在中性点不接地系统中故障线路的零序电流是从线路流向母线的方向非故障线路则是从母线流向线路两者呈现明显的180度相位差。在仿真中我从两条出线的测量模块中分别提取零序电流经过零序滤波器后比较它们的相位。相位差接近180度的就是故障线路这个方法在模型里能准确选线。但换到经消弧线圈接地系统后这个方法就不灵了。由于消弧线圈补偿故障线路和非故障线路的零序电流相位差不再保持180度幅值差异也不明显这就是实际消弧线圈系统选线困难的根源。在仿真里我把这两种情况做了对比用图形清晰地展现给读者不接地系统的方向判据清晰消弧线圈系统的判据变得模糊。这直接说明为什么消弧线圈系统需要更高级的选线算法比如注入信号法、零序电流有功分量法或暂态量法。基于这个仿真模型你还可以验证一种改进思路在消弧线圈支路中临时控制补偿度增大或减小脱谐度观察故障线路零序电流方向或幅值是否出现特殊变化。这种方法本质上是通过主动改变系统参数来“制造”故障特征比被动测量更可靠。我在另一个课题里用这个方式做了仿真模型的灵活性足够支撑这类控制策略的验证。5. 常见仿真问题排查与避坑经验5.1 仿真不收敛与步长设置问题Simulink电力系统仿真最常见的错误就是“Solver failed to converge”或者系统运行到一半直接发散。这类问题我排查过很多次一般有三个原因。第一是仿真步长太大。电力电子开关器件和故障断路器在动作瞬间会产生陡峭的波形如果最大步长大于us级别求解器可能无法精确捕获这些快速变化导致迭代发散。我在模型中把最大步长设为1e-4秒后基本就没再出现发散问题如果模型里还有更快的开关器件建议继续收紧到1e-5秒或1e-6秒。第二是电感电流和电容电压的初始条件不一致。Simulink电力系统模块通常默认从零初始状态开始仿真但一个变压器和线路组成的网络如果在零时刻突然接入10kV电源会产生很大的冲击电流。我在变压器和电源之间加了一个初始相位控制并用solver configuration里的“Start simulation from steady state”选项进行初始化让系统从稳态解开始运行这一步效果非常明显。第三是纯电感或纯电容回路导致的代数环问题。消弧线圈直接并联在变压器中性点和大地之间时如果不带任何串联电阻理想电感在直流稳态下会形成短路路径求解器很容易报错。给消弧线圈串一个0.01到0.1Ω的小电阻问题自然解决。5.2 故障模块和测量模块的配置陷阱Three-Phase Fault模块的配置有几个容易踩的坑。首先故障类型选择时一定要确认是否勾选了“Ground”故障。如果不勾选故障只在相与相之间短接不会产生零序电流波形完全不对。其次故障阻抗不要设成绝对零用0.001Ω到0.01Ω之间的金属性接地电阻更符合实际也避免数值奇异。在测量侧很多人直接用Three-Phase V-I Measurement测出来的电压除以根号3当相电压用这是错误的。该模块输出的是线电压还是相电压取决于你勾选的“Voltage measurement”是“phase-to-ground”还是“phase-to-phase”。我习惯选“phase-to-ground”这样三相电压直接对应实际相电压非故障相电压抬升效果一目了然。如果选“phase-to-phase”看到的波形会一直是线电压相电压升高的特征就被掩盖了。零序电压的提取也不要指望直接在测量模块里找到。需要把三相电压通过一个Fourier分解模块或者直接搭建数学运算单元三个电压求和除以3来得到。我建议用后一种方式因为可以实时观察3U0的时域波形比FFT更直观。5.3 消弧线圈参数调整后的重新校验消弧线圈的电感值不是一次算好就可以一直用的。如果你改变了系统电压、线路长度或负荷大小系统的对地电容电流会变化原来整定的消弧线圈补偿度就失效了。我在换了一条更长的电缆线路后就发现故障点电流从1A又涨到了5A明显是欠补偿。排查思路是这样的先在新工况下测量单相接地故障时的故障点电流不经过消弧线圈得到新的Ic值。然后按过补偿8%计算新的IL再反解出电感量。最后把消弧线圈的电感参数替换掉重新跑一遍验证残流是否降到合理范围。这个流程要养成习惯否则仿真结果很容易与实际脱节。另外不同文献里提供的线路零序电容参数差别可能很大架空线和电缆差距是一个数量级。如果你是在做对比研究建议统一采用同一套线路参数并在文中说明参数来源这对论文评审或者工程交流都是负责任的做法。做电力系统仿真的这几年我最大的感受是模型永远是为物理概念服务的。小电流接地系统的仿真看似不复杂但如果对电容电流路径、消弧线圈补偿度、零序电压特征这些基本原理理解不透模型搭得再漂亮结果也是空中楼阁。我这个模型在参数上做了一些简化比如忽略了负载动态特性、单相接地用固定过渡电阻模拟但对大多数故障特征研究和保护策略验证来说精度已经完全足够。后续如果你需要做谐振过电压、弧光接地或选线算法对比只需要在这个模型基础上做增量扩展不需要推倒重来。