做配电网继电保护的人十有八九都得跟“小电流接地系统”打交道。这类系统最典型的问题就是单相接地故障——占了配电网故障的八成左右可故障电流只有几安到几十安常规过流保护根本没法可靠动作。我刚接触这个方向的时候最头疼的就是现场没法随便拿真系统做接地试验所以Simulink仿真几乎成了研究这类问题的标配。这篇文章就围绕一套实际跑过的仿真模型来拆解中性点不接地和经消弧线圈接地两种接地方式下的单相接地故障仿真从系统原理、建模参数到波形分析、常见坑位尽量一次讲明白。适合正在做配电网继电保护课题的学生、刚接触小电流接地选线的工程师以及所有想用Simulink把10kV配电网故障过程“看个清楚”的人。仿真这件事看着是搭模型实际上是先把“选线为什么难”这个问题想清楚。故障电流小、方向规律不明显、消弧线圈投入后特征更复杂这些都是小电流接地选线的痛点。用Simulink把两种接地方式的故障响应仿真出来就能把痛点一个个变成看得见摸得着的波形后面再谈选线算法、消弧线圈优化才有底气。1. 小电流接地系统故障仿真的整体思路1.1 两种接地方式的本质区别先明确概念。所谓小电流接地系统指的是中性点不接地、经消弧线圈接地、经高阻接地这几类系统。它们本质上都是中性点不直接接地所以单相接地时故障电流很小。本文讨论的两种方式——中性点不接地和经消弧线圈接地——是最常见、也最有对比价值的研究对象。中性点不接地很好理解变压器中性点没有任何人为接地整个系统对地只有分布电容这条回路。单相接地时故障点流过的电流来自系统内所有线路的对地电容电流大小取决于系统规模、线路长度和线路类型。架空线路电容小电流可能只有几安电缆线路电容大电流可以到几十安甚至上百安。经消弧线圈接地则是在变压器中性点接入一个电感线圈。单相接地时消弧线圈产生感性电流与系统的容性电流在故障点相互抵消从而把故障电流降到更小电弧更容易自行熄灭避免故障扩大成相间短路。这个思路本质上和“给故障点的电容电流找个感性对头”是一样的道理。但这两种方式在仿真里的差别可不止“多加一个线圈”这么简单。中性点的引出方式、零序回路的结构、故障后的电压分布、零序电流的方向规律全都会变。建模型之前如果脑子里没有一张完整的系统拓扑图很容易出现“波形跑出来但不知道对不对”的情况。1.2 单相接地故障后的电气量变化规律无论哪种接地方式单相接地后的基本规律是一致的故障相电压降为0非故障相电压升高为线电压也就是原来相电压的√3倍中性点电压升高为相电压。这个规律是判断仿真结果是否正确的第一道标准。区别在于零序电流。中性点不接地系统里单相接地时所有非故障线路的对地电容电流都会通过母线流向故障线路再从故障点流回系统。所以非故障线路的零序电流等于本线路自身的对地电容电流数值小故障线路的零序电流等于所有非故障线路对地电容电流之和数值大且方向与非故障线路相反。这就是“零序电流幅值最大、方向相反”选线判据的物理基础。经消弧线圈接地之后消弧线圈的感性电流会叠加到零序回路里。如果感性电流与容性电流完全抵消就是全补偿故障点电流理论上为零但实际为了避免串联谐振一般都设计成过补偿也就是感性电流略微大于容性电流。这时故障线路的零序电流变成感性分量占主导方向反而可能与非故障线路相同幅值也变得很小。传统幅值比较法在这个场景下基本失效这也是为什么消弧线圈系统选线特别麻烦。1.3 为什么要用Simulink做对比仿真选线算法要验证不能只靠理论推导。拿真实系统做单相接地试验要么涉及停电要么需要专门的安全措施成本很高。Simulink的价值就在于用几分钟就能搭出一个10kV配电系统自由切换接地方式、故障位置、故障电阻甚至过渡电阻从0.01Ω到1000Ω都能模拟波形、向量、谐波分析一应俱全。另一个好处是可控性强。仿真里可以把消弧线圈设置成全补偿、欠补偿、过补偿观察不同补偿度下的零序电流变化这在真实系统里要反复调整线圈分接头才能做到。先仿真摸清规律再去现场看数据效率会高很多。所以很多论文和课题的第一步都是先建一个可靠的接地故障仿真模型这篇文章讲的就是这一步。2. Simulink建模前的系统参数设计2.1 仿真模型整体架构与模块选型模型整体架构必须在一开始想清楚不能一边搭一边改。按照10kV配电网的典型结构我习惯把模型分成四块电源与变压器、两条及以上出线、故障设置模块、测量与分析模块。如果是经消弧线圈接地系统还要在变压器中性点加一条消弧线圈支路。Simulink里用的库主要是Simscape Electrical的Specialized Power Systems旧版本叫Power System Blockset。核心模块如下Three-Phase Source系统电源用来模拟110kV侧等效电源Three-Phase Transformer (Two Windings)降压变压器110kV/10kVThree-Phase PI Section Line输电线路模型参数按单位长度输入Three-Phase Fault单相接地故障模块Series RLC Branch消弧线圈支路用R、L串联模拟Three-Phase V-I Measurement测量断路器处三相电压电流Voltage Measurement、Current Measurement单相或零序测量RMS、Fourier、Scope后处理与波形观察模块选型看似简单其实有个关键点线路一定要用三相PI Section Line不要用三个单相模块拼。因为零序参数的设定在单相模块里很容易搞错而三相PI Section Line可以直接输入正序和零序R、L、C参数符合实际线路参数的习惯。2.2 10kV配电网典型参数设定参数设计是仿真能不能反映真实系统的核心。我以一个典型10kV辐射状配电网为例参数如下系统电压10kV频率50Hz变压器110kV/10kV5000kVA短路电压百分比Uk%10.5连接组别Dyn11出线110km架空线正序参数R10.17Ω/kmL11.2mH/kmC19.7nF/km零序参数R00.23Ω/kmL04.1mH/kmC06nF/km出线25km电缆线路正序参数R10.13Ω/kmL10.25mH/kmC1300nF/km零序参数R00.28Ω/kmL01.1mH/kmC0180nF/km故障设置在出线1上距母线5km处A相接地这里有个容易被忽略的点Simulink PI Section Line模块里输入的C参数是每相的对地电容单位F/km。但很多资料里给的线路电容是“线间电容”或者“正序电容”要换算清楚。三相PI Section Line模块直接用每相对地等效电容即可正序电容和零序电容分别在对应参数栏里输入。架空线和电缆的零序电容差异很大这直接决定了系统电容电流的大小。电缆的C0比架空线大一个数量级以上所以含电缆的系统里单相接地电容电流主要由电缆贡献。做配电网仿真时不能为了省事把所有线路都设成架空线那样算出来电流会很假。2.3 消弧线圈电感值的计算方法消弧线圈参数是经消弧线圈接地模型的关键不能随便填一个电感值就开始仿真。计算思路是先算系统总的对地电容电流再根据补偿度确定消弧线圈的感性电流最后反推电感值。单相接地时系统总电容电流Ic的估算公式Ic √3 × ω × C_eq × U_N × L_total其中ω是角频率314rad/sC_eq是单位长度等效对地电容U_N是线电压L_total是线路总长度。这个公式推导自故障后非故障相电压升高为线电压、对地电容电流矢量合成的结果建议直接记。按上面的参数估算出线1架空线10kmC06nF/km出线2电缆5kmC0180nF/km。总电容电流约为Ic 1.732 × 314 × (10×6×10⁻⁹ 5×180×10⁻⁹) × 10000 ≈ 5.1A取过补偿度8%消弧线圈感性电流IL Ic × 1.08 ≈ 5.5A。10kV系统相电压Uφ 10000/√3 ≈ 5774V消弧线圈电感值L Uφ / (ω × IL) 5774 / (314 × 5.5) ≈ 3.34H。实际模型中Series RLC Branch按此计算值填入L同时串联一个小电阻模拟线圈内阻一般取电感感抗的2%~5%也就是几欧到十几欧。这个内阻很重要它决定了残流的有功分量是很多选线算法比如有功分量法的信号来源不能设成0。消弧线圈的全补偿、过补偿、欠补偿对比仿真也很有价值。全补偿时L Uφ/(ω×Ic)此时系统处于谐振状态零序电压会异常放大实际工程中要避免欠补偿时L比计算值大过补偿时L比计算值小。通过修改L值就能在一套模型里观察不同补偿度下选线特征量的变化。3. 模型搭建与关键环节实现3.1 电源、变压器与中性点引出电源用Three-Phase SourceInternal connection选Y中性点接地或不接地视上层系统而定。我这里模拟110kV系统设为中性点直接接地比较合理接入110kV母线。变压器是模型里最容易出问题的地方。用Three-Phase Transformer (Two Windings)高压侧110kV接线方式选Delta低压侧10kV选Yg也就是星型中性点引出。为什么低压侧必须选Yg因为无论是消弧线圈还是中性点不接地系统的中性点测量都要从中性点端子引线。如果选成Delta或者不带g的Y中性点根本没引出后续接消弧线圈、测中性点电压都是空谈。这是新手特别容易卡住的地方仿真前一定检查变压器对话框里的接线方式设置。变压器参数里还要正确填入容量、频率、绕组电压、短路电阻和短路电感。这些参数由前面的铭牌参数换算得到Simulink支持直接输入标幺值电压和标幺值电阻填起来比较方便。3.2 线路模型与单相接地故障设置两条出线分别用两个Three-Phase PI Section Line搭建设置对应的正序、零序参数和线路长度。PI Section Line内部默认用集中参数π模型多条π段串联可以提高分布特性精度但对单相接地分析来说单段或两三段已经够用。故障模块用Three-Phase Fault在Fault Selection里勾选Phase A Fault和Ground FaultTransition Time设为0.05s也就是仿真开始后0.05s发生A相接地。故障持续0.2s到0.25s结束方便观察故障稳态和切除后恢复两个阶段。故障电阻的设置在Fault Resistance栏里。金属性接地取0.001Ω过渡接地根据场景可以取10Ω、100Ω甚至1000Ω。高阻接地时零序电压显著降低、电流很小更贴近实际现场难选线的情况建议把不同故障电阻都跑一遍对比差异。故障位置放在线路1距离母线5km处。想要观察不同位置对选线的影响只需要修改线路模型段数分配。比如把10km线路拆成两段故障接在两段之间调整前段长度就能模拟不同故障距离这种方式很灵活。3.3 消弧线圈支路的连接与参数录入如果是中性点不接地系统变压器低压侧中性点直接悬空什么都不接。如果是经消弧线圈接地系统则需要在中性点和地之间接一个Series RLC Branch。Series RLC Branch的Branch Type选RLC串联R填5ΩL填3.34HC填inf无穷大表示不接电容。这里注意C值默认是1e-6F必须改成inf否则它默认带一个电容仿真结果会错得离谱。还有一个隐藏细节消弧线圈支路最好并联一个高阻电阻作为阻尼电阻。Simulink里电感支路直接接入地数值积分时容易产生高频数值振荡并联一个高阻电阻比如消弧线圈感抗的10倍以上可以有效抑制数值振荡同时不影响稳态特性。我一般并联一个500Ω左右的高阻电阻波形会干净很多。测量模块方面在变压器低压侧出口用Three-Phase V-I Measurement测母线电压和电源侧电流在线路1、线路2末端各加一个Three-Phase V-I Measurement测线路电流便于分别观察故障线路和非故障线路的零序电流。3.4 零序电压电流的提取与显示零序量的提取是波形分析的基础。Simulink里最直接的方法是用相加模块对三相电压或三相电流求和。零序电压U0 (Ua Ub Uc) / 3 零序电流I0 (Ia Ib Ic) / 3但实际工程里经常看3U0和3I0也就是三相之和本身。这个不影响分析只要保持一致即可。具体实现上把Three-Phase V-I Measurement的Uabc输出接一个Demux取Ua、Ub、Uc三路用两个Sum模块相加再乘1/3增益就得到U0。电流一样处理。如果嫌麻烦也可以用Power MeasurementDemux的“Sequence Analyzer”模块直接输出正序、负序、零序分量但那个模块在部分Simulink版本里位置比较深入门先用Sum加增益的方式最稳妥。波形观察用Scope就能满足需求。如果想进一步做频谱分析把零序信号接FFT或Fourier模块提取基波幅值和相位。选线算法里经常用零序电流的幅值、相位、有功分量这些量都可以通过Fourier模块得到。我建议Scope里同时显示以下信号三相电压、中性点电压、故障线路零序电流、非故障线路零序电流、消弧线圈电流。一次性全部铺开故障前后的变化一目了然。4. 仿真结果与波形分析4.1 中性点不接地系统的单相接地故障响应先跑中性点不接地模型。求解器用ode23tb最大步长设为1e-4s仿真时长0.3s故障从0.05s开始到0.25s结束。故障后波形应当符合以下几个特征第一A相电压迅速跌到接近0B相和C相电压升高到线电压水平也就是从约5.77kV升到约10kV。故障瞬间非故障相电压会有轻微的暂态过冲随后进入稳态这是对地电容充放电造成的属于正常现象。第二中性点电压从0升到相电压5.77kV和理论一致。第三零序电流的对比是关键。故障线路出线1的零序电流明显大于非故障线路出线2大约等于非故障线路电容电流之和方向相反。这个方向怎么判断看波形的相位关系故障线路零序电流与非故障线路零序电流相位差接近180°。这个结果对应了“比幅比相”选线法的物理基础。我每次仿真完先看这一点如果方向不对或者幅值关系不对先检查测量极性再看线路参数而不是急着调消弧线圈。4.2 经消弧线圈接地系统的单相接地故障响应接着改模型把变压器中性点经消弧线圈接地支路接入保留原来的故障设置。同样的求解器和步长设置再跑一遍。最直观的变化是故障点电流明显减小。中性点不接地时故障线路零序电流约为系统总电容电流5.1A左右接入消弧线圈后变成残流只有零点几安到一两安幅值差距非常大。更值得注意的是方向。在过补偿状态下故障线路零序电流的方向和非故障线路零序电流变成同向而不是反相。原因前面说过过补偿时感性电流大于容性电流故障线路零序电流由感性分量主导整体相位翻转。这个现象直接影响选线判断。很多基于“零序电流方向和幅值比较”的传统选线方法在消弧线圈系统里会失效或者出现误判。仿真波形把这个过程展示得非常清楚比干讲理论要直观得多。4.3 两种接地方式的关键特征对比我把两种方式的仿真结果整理成一个对比表方便各位在做课题报告或论文时直接参考对比项中性点不接地经消弧线圈接地过补偿故障相电压降低到接近0降低到接近0非故障相电压升到线电压升到线电压中性点电压升到相电压升到相电压故障线路零序电流大约等于系统电容电流之和小等于残流非故障线路零序电流小约等于本线路电容电流小约等于本线路电容电流故障线路与非故障线路零序电流相位反相同相过补偿时选线判据效果比幅比相法有效传统比幅比相法失效需要功率方向或有功分量法这个表格基本就是模型的“结论页”。仿真的价值不在于跑出一堆波形截图而在于把这些特征梳理清楚。实际项目里我通常先把这张表做出来再决定选线算法往哪个方向设计。5. 仿真中的常见问题与排查技巧5.1 仿真不收敛或波形高频振荡这是Simulink电力电子和电力系统模型里最常遇到的一类问题症状是波形充满毛刺或者干脆报错“Simulink cannot solve the algebraic loop”。排查方向有三个。第一换求解器。小电流接地系统里既有电感又有电容属于典型的刚性系统用ode45很容易跑不动或精度差。我一般用ode23tb或ode15s这两种变步长求解器对刚性系统比较友好。第二限制最大步长。变步长求解器默认最大步长可能太大导致故障瞬间的暂态过程被跳过或产生数值振荡。把Max Step Size设为1e-4s甚至1e-5s波形会明显改善。第三消弧线圈支路并联阻尼电阻。如果用了消弧线圈后出现持续振荡大概率是纯电感接入地导致的数值问题。并联一个高阻电阻稳压不影响稳态结果但能让数值计算稳定很多。5.2 零序电流方向判断反复出错新手最常踩的坑是零序电流相位看着不对。原因往往很简单测量模块里电流方向设置反了。Three-Phase V-I Measurement模块的电流正方向默认从端流入、-端流出。如果电流方向接反相位整体翻转180°故障线路看起来就变成了“非故障线路”。判断方向之前先确认所有测量点的电流正方向都是一致的最好都朝向负荷侧或者都朝向母线侧再比较相位关系。还有一点经验不要直接在Scope波形图上用肉眼判断方向而是用FFT或Fourier模块提取零序基波分量的相位角用数值对比。相位差接近0°就是同相接近180°就是反相这样判断准确、效率高。5.3 变压器中性点引不出来调模型到经消弧线圈接地时发现中性点没有端子可以接线基本可以断定变压器接线方式设错了。反复强调一下低压绕组必须选Yg也就是带接地的星形接法而不是Y或Delta。如果变压器选的是Three-Winding类型也要确认绕组配置。另外要注意变压器中性点引出后不要直接接地否则就变成中性点直接接地系统跟小电流接地完全不搭边。中性点要么悬空要么经消弧线圈接地别搞混。5.4 过渡电阻太大导致选线特征消失用小电阻做故障仿真一切正常把故障电阻改成100Ω甚至1000Ω后零序电压和零序电流的幅值都会大幅下降这时候波形看起来几乎和正常状态没什么区别。这不是模型坏了而是高阻接地的真实情况。工程现场最难处理的就是这类高阻接地故障因为故障特征弱保护装置启动都困难更别说选线。如果要做高阻接地研究建议在测量环节加装“启动判据”模块比如零序电压越限触发录波然后专门观察越限后的暂态零序电流。这部分内容属于选线算法层面但仿真模型必须先把高阻场景跑出来算法才有验证基础。5.5 故障切除后的恢复过程不用忽略很多仿真只盯着故障期间的波形忽略了0.25s故障切除后的恢复过程。实际上故障切除瞬间系统对地电容重新充放电非故障相电压会出现一个短暂的过冲零序电流也会有一个衰减振荡。这个过程对研究故障熄弧、系统绝缘配合也有参考价值建议Scope时间轴拉长到0.3s以上把恢复过程一起录下来。我在实际项目中往往会再往后设置几个不同故障距离和故障电阻的仿真case把波形数据导出到MATLAB工作区批量计算零序电流的幅值和相位做成表格或散点图。这样做的好处是后面写选线算法时数据都是现成的不用反复重新跑模型。最后再分享一个操作习惯建模时给每个测量信号起好名字用Goto标签把关键信号送到同一个Scope里显示而不是在模型里到处拉线。模型复杂之后拉线越多越容易接错信号经过Goto/From传递要清晰得多排查起来也快。这套模型改造成本很低但能让你在后续换线路参数、换补偿度、加算法时省下大量时间。