干配电网继电保护的同仁应该都有过这种体验母线发单相接地报警选线装置给出一个“疑似故障线路”结果现场试拉了一轮真正故障的却是另一条线。尤其是在中性点经消弧线圈接地的小电流系统里故障残流小、相位关系飘忽不定选线这个老问题硬是坑了一批又一批人。我的做法是与其靠经验去猜不如直接用MATLAB/Simulink把不同中性点接地方式下的单相接地故障仿真模型搭出来用波形图和零序分量数据说话。这套模型不仅能复现“故障线路零序电流最大”“方向相反”这些理论结论还能帮我们验证选线判据在暂态和稳态情形下的实际表现算得上选线方案开发前最划算的一笔投入。本文就围绕这套模型展开适合正在做配电网故障分析、选线装置算法验证或电力系统课程设计的读者参考。1. 单相接地选线这么难仿真模型能帮上什么忙1.1 选线难题的根源故障量实在太小小电流接地系统中性点不接地、经消弧线圈接地发生单相接地后故障点流过的仅是系统对地电容电流最小的几安培大的也就几十安培比负荷电流小一两个数量级。而且三相线电压仍保持对称负荷照常运行继电保护不像相间短路那样立即跳闸。规程允许带故障运行一段时间可值班员心里都清楚非故障相电压已经升高到线电压绝缘薄弱点随时可能发展成两相短路电弧还会引起过电压。选线的目的就是在众多出线里找到发生接地的那一条然后安排转移负荷或者停电处理。难点在哪儿呢稳态情况下零序电流的幅值和相位在消弧线圈补偿后变得很不明显。特别是过补偿运行时故障线路的零序电流方向跟不接地系统相反如果选线装置还沿用稳态幅值或方向判据就容易选错。暂态过程虽然信息丰富但衰减快、频率高传统电磁式PT/CT难以准确传变也给工程实现增加了不少难度。说到底选线本质上是一个“信号弱、干扰大、判据适应性差”的问题需要把故障机理吃透才能解决好。1.2 仿真模型能提供什么价值搭建仿真模型不是简单去验证课本结论而是给自己提供一个可以反复折腾的“试验台”。比如说我想研究消弧线圈补偿度从欠补偿到过补偿零序电流相位到底怎么翻转想看看故障发生时刻在电压过零和电压峰值时暂态分量有多大区别想验证五次谐波判据在长线路占比很高的电网上还灵不灵——这些在真实电网里动一次刀很贵但在Simulink里改一个参数、点一下运行波形和录波数据立刻就能出来。另一个关键价值在于算法验证。选线算法的开发往往需要大量样本数据来整定阈值和评估鲁棒性单靠现场录波肯定不够。仿真模型配合Matlab脚本可以批量修改故障线路、故障位置、故障电阻、故障时刻、系统运行方式等参数一次跑出几百个案例导出的零序电流零序电压数据直接喂给选线判据做统计分析。这种“仿真算法验证”的工作流是目前我做过的最省力的选线研究方式。2. 中性点三种接地方式在Simulink里的建模区别2.1 中性点不接地最容易被忽略的对地电容建模前先把系统参数定下来。以最常见的10kV配电网为例三相电源10kV/50Hz变压器容量选6.3MVA或10MVA高压侧Y接线中性点引出。不接地就是变压器中性点悬空不做任何连接Simulink里对应的是三线制或者带中性点但不连接。不接地系统最容易忽略的是线路对地电容。单相接地电流完全由线路和负载侧对地电容形成所以线路模型里的零序电容参数必须给准。我一般用集中参数π型线路每条出线取典型单位参数正序电阻0.27Ω/km正序电感1.33mH/km正序电容0.339μF/km零序电阻2.7Ω/km零序电感4.13mH/km零序电容0.28μF/km架空线和电缆差别很大具体以实际线路手册为准。线路长度按5km、8km、12km、20km设置几回线这样各线路的对地电容电流大小就有明显区分度后面波形对比才好看。2.2 经消弧线圈接地补偿度是建模的灵魂消弧线圈接在变压器中性点和地之间模型里就是一个电感串联一个阻尼电阻通常还会带一个小电阻模拟零序阻尼。电感值要根据系统对地电容电流来算需要补偿的容性电流I_C约等于所有线路对地电容电流之和消弧线圈在工频下的感抗X_LU_ph/I_C。如果按过补偿10%设计实际补偿电流I_L1.1×I_C感抗X_LU_ph/(1.1×I_C)电感LX_L/(2πf)。我举个例子假设系统对地电容电流总和为30A10kV系统相电压约5774V过补偿10%后补偿电流33AX_L≈175ΩL≈0.557H。阻尼电阻通常取感抗的10%~15%约17.5Ω到26Ω。仿真时把这个RLC串联支路放到中性点与地之间即可Simulink里用Series RLC Branch或者Three-Phase Series RLC Branch注意电感后必须可靠接地。如果要研究五次谐波选线最好给变压器或电感做饱和特性不然谐波源缺失后面验证判据会落空。2.3 经电阻接地选线判据反而最简单中性点经电阻接地分为小电阻接地几十欧以下和中电阻接地100~1000Ω。小电阻接地时单相接地故障电流可达数百安培故障特征非常明显零序电流保护可以直接动作跳闸选线压力不大。这个时候中性点与地之间只接一个三相电阻支路电阻值按工程实际取6Ω、10Ω或20Ω都行。中电阻接地则介于两者之间主要用于限制过电压同时提供一定的零序电流。我把三种接地方式搭成同一个系统的变子部分也就是用三个不同的中性点支路模块通过手动开关或者变子系统Variant Subsystem切换。这样做的好处是线路和故障参数完全一致只改中性点支路就能在同一个Simulink环境下对比三种接地方式的零序特征不需要复制三个大模型后面做参数扫描也清爽得多。3. 从零开始搭一套选线仿真模型关键模块与参数3.1 系统级拓扑电源、变压器、线路、负载我习惯用Simscape Electrical里的Specialized Power Systems库老版本叫SimPowerSystems。搭建顺序三相电源→变压器→母线→多条π型线路→负载故障模块放在需要考察的出线上测量模块放在母线出线断路器处。电源参数比较容易踩坑。Three-Phase Source里填的是线电压有效值比如10kV系统填10e3频率50Hz接线方式选Y接中性点不接地。变压器用Three-Phase TransformerTwo Windings高压侧Yg低压侧Δ高压侧中性点引出后接入中性点支路。如果变压器模型选Yg就会自带一个中性点端子方便接消弧线圈或电阻。线路我用Three-Phase Pi Section Line每段按5km来串联两段就是10km这样既能近似长线的分布效应又不至于让仿真慢到没法忍。线路末端接三相并联RL负载参数按每条线路所带负荷设置比如5MWj0.5Mvar。3.2 故障与测量模块的放置故障用Three-Phase Fault模块选择A相接地在Phase-A Fault那栏打勾。故障时刻我习惯用外部控制信号接一个Step模块0.2s之前为00.2s时跳变到1触发故障。故障电阻设成0.1Ω金属性接地也可以做成变量在Matlab脚本里循环赋值专门跑高阻接地场景。测量方面母线处接Three-Phase V-I Measurement把三相电压电流引出来。为了得到零序电压和零序电流最干净的做法是在Matlab脚本里合成i0(iaibic)/3u0(uaubuc)/3。注意三相信号必须时间对齐离散仿真里如果不同数据点没有对齐会有误差所以我把To Workspace的采样时间设成跟模型步长一致并保存为Structure With Time格式后面处理比较方便。Simulink里也有零序分量提取模块Sequence Analyzer但我觉得脚本合成更灵活、不占模型资源。3.3 求解器与仿真参数的设置这是很多模型跑不出理想波形的关键。我推荐用离散仿真Powergui的仿真类型选Discrete采样时间设成10μs或20μs这样才能捕捉到几百赫兹到几千赫兹的暂态高频分量如果只关心稳态基波波形50μs也够。求解器选ode23tb或ode15s容差保持默认。仿真时长设置0.4s左右前0.2s建立稳态0.2s发生故障后面0.2s用于观察暂态和稳态。如果要做暂态能量判据建议把波形记录窗口覆盖故障前至少一个周波方便对齐零序电压过零点。核心参数参考如下表元件参数典型值三相电源线电压有效值10kV / 50Hz变压器高压侧中性点Yg引出架空线正序R1/L1/C10.27Ω/km, 1.33mH/km, 0.339μF/km架空线零序R0/L0/C02.7Ω/km, 4.13mH/km, 0.28μF/km消弧线圈过补偿度10%中性点电阻小电阻接地10Ω仿真采样时间离散Ts10~20μs4. 故障波形里藏着哪些选线特征重点看零序4.1 不接地系统的波形幅值最大方向相反以4条出线为例线路1长20km电容电流最大线路2长12km线路3长8km线路4长5km。假设线路4在0.2s发生A相金属性接地。母线零序电压从0跃变到接近相电压峰值约5774V的幅值。观察各线路首端零序电流非故障线路1/2/3的零序电流是各自线路的对地电容电流幅值约几安到十几安相位超前零序电压90°容性方向故障线路4的零序电流等于其余三条线路零序电流之和幅值明显最大相位滞后零序电压90°从线路流向母线。示波器上看起来故障线路4的零序电流和线路1的零序电流波形几乎反相这就是稳态幅值加方向选线的基础。这里要提醒一点如果故障线路恰好是短线路而某条正常线路很长故障线路电容电流相对非故障线路优势就不大。比如线路4只有5km线路1有20km故障线路零序电流幅值理论上是线路123之和仍然大于任意单条非故障线路所以幅值比较在理论上成立但现场CT误差大实际还是倾向用方向判据或者暂态判据。4.2 经消弧线圈接地后的相位翻转现象同样的故障场景把中性点换成消弧线圈过补偿10%后故障线路零序电流变成补偿后的残流。因为过补偿时电感电流大于系统总电容电流故障线路的零序电流方向变成与不接地系统相反超前零序电压90°感性方向。从波形上看故障线路4的零序电流相位与正常线路1-3的相位相同都是超前如果还套用“方向相反”的判据必然误判。幅值上故障线残流可能很小有时候甚至小于某条长正常线路的电容电流稳态幅值判据也不可靠。如果改成欠补偿故障线路零序电流方向又会变回容性方向但系统运行方式一变补偿度就会漂移稳态特征非常不稳。这时候必须看暂态故障开始的最初几个毫秒零序电压突变量和零序电流突变量里包含丰富的暂态分量其中故障线路的暂态零序电流首半波极性与零序电压相反非故障线路相同。在Simulink波形里能看到故障线路的高频毛刺明显更大尤其在相电压峰值附近接地时。4.3 经小电阻接地时波形清晰得令人感动中性点经10Ω电阻接地单相接地故障电流主要由中性点电阻提供的零序电流构成数值可达几百安。各线路零序电流中故障线路的零序电流幅值显著大于非故障线路方向相反而且基波分量非常稳定。波形上几乎不用滤波直接用基波有效值幅值比较就能选出故障线。这也是为什么很多城市配电网改成小电阻接地——代价是接地故障必须立即跳闸供电可靠性确实降低了但保护的选择性和灵敏性都大幅改善。三种接地方式波形特征对比如下表接地方式故障线路零序电流幅值方向特征稳态判据可靠性不接地等于其余线路电容电流之和与非故障线路相反较好消弧线圈过补偿残流可能很小与非故障线路相同差需暂态/谐波小电阻接地远大于电容电流与非故障线路相反非常好5. 选线判据在Simulink/Matlab里的实际验证手法5.1 零序电流幅值/方向判据的脚本化实现模型跑完后我把To Workspace导出的三相电流数据放到Matlab脚本里合成零序电流然后计算基波幅值和相位。幅值判据就是比较各线路零序电流的基波有效值方向判据需要结合零序电压相位看零序电流是超前还是滞后。代码框架大致这样% 假设 ts 为采样时间ia1,ib1,ic1 为线路1的三相电流信号 i0_1 (ia1.data ib1.data ic1.data) / 3; Fs 1 / ts; N length(i0_1); Y fft(i0_1); % 基波对应的频率分量索引忽略直流 f (0:N-1) * Fs / N; idx50 find(abs(f - 50) 1); amp1 abs(Y(idx50)) / N * 2; % 基波幅值 phase1 angle(Y(idx50)); % 基波相位不过实际我更推荐用Powergui的FFT Analysis工具它可以直接看选定信号的基波幅值相位不用自己在脚本里算。批量验证判据时再写循环脚本批量提取数据统一计算。核心是数据要干净负序、零序分量没有混叠所以采样时间务必要小。5.2 五次谐波判据的仿真条件五次谐波判据针对消弧线圈系统原理是消弧线圈只在基波下补偿对五次谐波感抗变成基波的5倍所以五次谐波电流基本不补偿仍按容性分布故障线路的五次谐波电流方向与非故障线路相反。要在仿真中造出五次谐波需要变压器饱和或者在线路中注入谐波源。我用的是让变压器模型勾选Enable saturation设置饱和特性参数让励磁电流在过电压时产生畸变自然带出5次谐波。我实测发现一个问题只靠变压器饱和产生的5次谐波非常弱只有基波百分之几加上测量噪声后不太容易分析。更可控的做法是用一个并联的谐波电流源Simulink里用Three-Phase Harmonic Suppressor或者直接加电流源在母线上注入1-2A的5次谐波零序电流模拟实际系统里的背景谐波。这样各线路5次谐波零序电流分布就很清晰故障线路方向相反的特征也明显。5.3 暂态首半波与暂态能量判据首半波原理接地瞬间的暂态电容电流远大于稳态电容电流且频率高故障线路与零序电压的极性关系在首半波时区别于正常线路。Simulink里要看到准确的暂态过程采样时间必须小于20μs。从波形上手动找首半波极性不难但批量验证就要写检测脚本先找零序电压故障后第一个过零点再看零序电流在这个区间里是正还是负。暂态能量法更抗干扰定义暂态零序能量为W ∫v0(t)·i0(t)dt在故障后一个或半个周波内积分。正常线路W为正或很小故障线路W显著偏大且符号相反。我通常用梯形积分函数trapz实现然后对大量故障案例统计W的分布整定选线阈值。仿真模型的优势在这里体现得很充分故障位置、故障电阻、故障起始角全都能扫样本量比现场录波不知道多到哪里去了。6. 我在这套模型上踩过的坑与调参经验6.1 线路测量极性搞反选线方向全反头一回搭模型时我把线路零序电流互感器极性接反导致正常线路和故障线路的方向特征完全颠倒选线判据跟着错位。后来逐个检查测量模块的正方向定义Three-Phase V-I Measurement在母线出线侧测线路电流正方向默认是母线流向线路故障点设置在线路上短路电流从两侧流向故障点。实际工程中零序CT的极性也有统一要求仿真阶段如果不把方向理清楚后面所有分析都白搭。检查方法很简单不接地系统单相接地时非故障线路零序电流是容性的超前零序电压90°如果画出来是滞后大概率极性反了。另外可以对比两条不同长度非故障线路的零序电流幅值应该大致跟线路长度电容成正比这一点也能帮助确认接线正确。6.2 消弧线圈电感整定不能只算线路电容我一开始计算消弧线圈补偿电流时只把各条线路线路电容电流加起来忽略了配电变压器低压侧对地电容以及母线侧杂散电容结果补偿度偏小故障残流巨大波形完全失真。后来用Powergui的Impedance vs Frequency功能扫了母线零序阻抗找到工频附近的谐振点从谐振点反推出系统实际总电容电流再回头调整电感参数模型才靠谱。具体做法是在母线处加一个很小的扰动电流源频率从1Hz到1000Hz扫频观察母线零序电压的幅值响应。谐振峰对应的频率和系统电容、电感参数直接关联。这个功能不复杂但非常实用建议搭完模型先跑一次扫频再设置消弧线圈参数。6.3 故障时刻和故障电阻对暂态波形的影响极大故障时刻选在相电压峰值与过零暂态分量差出一个数量级。峰值附近接地时du/dt最大暂态电容电流可达数十安培过零附近接地时暂态很小甚至看不清。故障电阻提升到1000Ω时暂态分量被明显阻尼首半波极性几乎淹没在噪声里。很多论文里给出的漂亮波形都是金属性接地加电压峰值时刻的极端工况实做时不能只看这种理想结果。所以我做验证时会用parfor批量扫描故障时刻从0.18s到0.22s步长0.002s故障电阻从0.1Ω到1000Ω按对数均匀分布每条线路都轮流做故障线。跑完一批数据统计各判据的准确率和失效工况再针对失效区域决定要不要加辅助判据。这套流程我已经重复了很多次是选线研究中性价比最高的一步。6.4 数据保存与内存优化的两个小技巧To Workspace模块如果同时保存十几路三相波形采样时间10μs仿真0.4s就是十几条×3×40000点内存压力不小。我的习惯是只保存需要分析的三相电流和母线三相电压其他中间量用Terminator或者干脆不接入记录模块。其次To Workspace输出格式建议选Structure With Time这样在脚本里直接取.data和.time字段不用再处理bus数据结构。另外离散仿真数据默认是均匀采样但要注意故障时刻0.2s和采样步长10μs不一定正好对齐导致暂态起始点有小偏差。影响不大但如果做首半波算法验证建议把故障触发信号也记录到WorkSpace确保取数据窗口时能精确对齐故障瞬间。最后分享一个让我效率翻倍的办法把三种中性点接地支路封装成Variant Subsystem用SS中性点选择变量控制到底接入消弧线圈、电阻还是悬空。我把这个变量设成Simulink参数配合脚本循环一套模型就能把三种接地方式、几十条线路参数、上百组故障工况全部跑完。选线仿真这件事模型本身不是核心核心是你能不能快速批量制造出足够多样的故障样本数据来检验你的判据。这套模型的价值就在于把原本需要现场蹲守才能积累的故障录波变成了按一下运行就能得到的东西。