配电网里的单相接地故障这几年在仿真交流群里被反复问起也是不少电气学生和运行人员最想搞明白的问题同样是接地为什么有的系统马上跳闸有的系统却能带故障撑一两个小时同样是故障电流为什么有些只有十几安培看起来一点都不像“短路”关键就在中性点的接地方式上。最近我做了一个Simulink环境下的仿真项目在同一个10kV系统模型里分别搭建了中性点不接地系统和经消弧线圈接地系统设置单相接地故障把两种接地方式下电压、电流、零序量的变化完整跑出来做了对比。这篇文章就把这个项目从建模思路、参数计算、实操步骤到踩坑心得完整拆一遍给需要做同类仿真对比的朋友做个参考。这个仿真适合谁看电力系统方向的研究生、配电网运行管理人员、搞继电保护整定的工程师还有做毕业设计的本科生。出发点很单纯用肉眼看得见的波形把教科书里那几句话变成脑子里的物理画面。中性点不接地系统的电容电流到底怎么流消弧线圈补偿的又是什么故障瞬间零序电压为什么突然飘高——跑完这几个仿真这些问题基本就通了。1. 项目思路拆解把两种接地系统放进同一个试验台1.1 两种系统在物理本质上差在哪中性点不接地系统的三相绕组中性点是悬空的或者说电位是漂浮的。正常运行时三相对称中性点电位为零各相对地电压就是相电压整个系统对地只有分布电容这一条通路。单相接地以后故障相通过接地点、大地、另外两相对地电容构成一个闭合回路这个回路里的电流就是电容反复充放电形成的容性电流。它的幅值一般不大几十安培以内但分布在整个网络上所有非故障线路的对地电容都会向故障点馈入电流。消弧线圈系统则是在中性点与大地之间串入了一个电感也就是消弧线圈。故障以后中性点电压发生位移消弧线圈两端承受的是位移电压于是电感中流过感性电流。这个电流方向和电容电流正好相反在故障点处相互抵消把接地电流压到很小。换句话说消弧线圈的本质不是堵住故障电流而是让感性电流和容性电流在故障点“碰面抵消”把故障点的能量降下来让电弧更容易熄灭。用大白话打个比方电容电流像一根水管渗漏出来的细流消弧线圈等于在旁边修了一个蓄水池把渗出来的水往回引漏水点的水量自然就小了。1.2 对比研究能回答哪些实际问题把这两种系统放在同一个模型里做故障对比能直接回答几个平时只靠概念很难想透的问题。第一为什么消弧线圈系统的零序电压保护整定思路和不接地系统不完全一样。两种系统在稳态电压特征上确实很像都是故障相电压下降、非故障相电压升高、中性点位移但电流的量级和相位不同导致可供检测的故障特征完全不同。第二为什么故障选线在不接地系统里相对好办到消弧线圈系统里就变难。不接地系统的故障线路和非故障线路零序电流方向相反、幅值差异明显很容易分辨。消弧线圈系统里残流太小方向又被补偿度影响传统的稳态选线判据经常失效。第三为什么消弧线圈系统对绝缘配合的要求更严格。消弧线圈把稳态电流压小了但故障消失瞬间和暂态过程中可能激发谐振过电压非故障相的电压幅值可能远超线电压这个在仿真里设置不同补偿度就能亲眼看到。选择Simulink做这件事是因为它元件库完整、波形可视化好、参数可以任意修改重跑。故障暂态过程靠手算几乎没法算仿真却能一帧一帧看清故障瞬间发生了什么——这对理解“单相接地故障时网络状态如何变化”特别有效。2. 建模前的关键准备参数计算与元件规划2.1 系统电压等级与模型边界设定仿真模型以10kV配电网为对象这是国内最常见的配电电压等级中性点非有效接地方式主要应用在这一层级。建模的第一步是确定模型边界也就是哪些部分要建、哪些部分可以省略。这里要研究的是中性点接地方式对单相故障特征的影响核心在于中压网络的对地电容回路低压侧负荷不是重点。所以模型可以做一个合理简化用一个Three-Phase Source模块设定10kV/50Hz等效系统电源线路用π型集中参数等效末端挂一个轻载电阻性负荷保证系统有功潮流正常最重要的对地电容参数集中在母线处设置。越简单的模型干扰因素越少越能看清故障的本质特征。如果要做多馈线对比可以在母线上挂两条或者三条出线一条设置故障另外一条正常运行这样才能观察到故障线路与非故障线路零序电流方向相反的现象。如果只建单条线路基本电压电流特征都能看到但零序方向选线的逻辑验证是不完整的。2.2 对地电容整个仿真的灵魂参数为什么说对地电容是整个仿真的灵魂因为中性点不接地系统的单相接地电流本质就是由系统各相对地电容决定的。公式很直接Ic 3ωC0Uph其中ω是角频率50Hz系统取314rad/sUph是相电压C0是系统单相等效对地电容。这个公式算出来的就是系统发生金属性单相接地时流过故障点的稳态电容电流。实际线路的对地电容分散在整个线路上仿真中通常做集中等效处理。架空线路单位长度对地电容大约是0.005到0.01μF/km交联聚乙烯电缆则要高得多每公里可以达到0.2到0.5μF。如果是一个包含架空线和电缆的混合配电网总等效C0取2到10μF都是正常的。我在这个项目里取C05μF这是很典型的配电网量级。有些初学者习惯把对地电容取得很小比如0.1μF结果故障电流只有零点几安培到一安培波形几乎看不出什么特征取得太大比如50μF电容电流上百安培又和中压配电网小电流接地系统的实际情况严重不符。C0取5μF左右能兼顾波形可见度和物理合理性。2.3 消弧线圈参数整定与计算脚本消弧线圈的电感量不是拍脑袋定的必须按补偿度计算。工程上消弧线圈一般运行在过补偿状态也就是让电感电流比系统电容电流略大通常取105%到110%的补偿度。选择过补偿的原因很实际系统运行方式变化时对地电容会改变如果运行在欠补偿状态某条线路断开导致电容减小系统可能偶然进入接近完全补偿的状态容易激发谐振过电压。过补偿则可以避免这个问题。计算流程是这样的先算系统电容电流再根据过补偿度算出需要的电感电流然后计算感抗和电感值最后配上阻尼电阻。可以直接把这个逻辑写成MATLAB脚本U_line 10e3; % 线电压单位V U_ph U_line / sqrt(3); % 相电压 f 50; w 2 * pi * f; C0 5e-6; % 单相对地等效电容单位F Ic 3 * w * C0 * U_ph; % 系统电容电流单位A compensation 1.10; % 过补偿度10% IL compensation * Ic; % 消弧线圈目标电流 XL U_ph / IL; % 消弧线圈感抗 L XL / w; % 电感值单位H Rd 0.1 * XL; % 阻尼电阻取感抗的10% fprintf(系统电容电流: %.2f A\n, Ic); fprintf(消弧线圈电流: %.2f A\n, IL); fprintf(消弧线圈电感: %.3f H\n, L); fprintf(阻尼电阻: %.1f Ω\n, Rd);以C05μF为例算出来系统电容电流大约是27安培过补偿10%后消弧线圈电流约30安培电感值约0.61H阻尼电阻约20Ω。这几个参数直接填进Simulink的RLC支路里就可以用。为什么不加阻尼电阻不行因为消弧线圈的电感在工频下感抗和系统对地电容的容抗非常接近本质上构成了一个并联谐振回路。如果没有阻尼任何一点扰动都可能激发持续振荡实际工程里这是谐振过电压的隐患仿真里则表现为波形剧烈震荡、数值发散。阻尼电阻相当于给谐振回路加了损耗把振荡能量消耗掉这也是实际消弧线圈装置必然串联阻尼电阻的原因。参数规划完成后做一个速查表方便对照参数中性点不接地系统经消弧线圈接地系统系统线电压10kV10kV对地等效电容C05μF5μF消弧线圈电感L不接入0.61H阻尼电阻Rd不接入20Ω补偿度无110%过补偿初始运行时间0~0.05s0~0.05s故障时间窗0.05~0.25s0.05~0.25s3. Simulink模型搭建实操全流程3.1 模型骨架元件库选择和连接顺序Simulink里做电力电子和电力系统仿真现在统一在Simscape Electrical早期叫SimPowerSystems元件库下。搭建这个模型需要用到的模块有Three-Phase Source三相电压源、Three-Phase PI Section Line三相π型线路、Three-Phase V-I Measurement三相电压电流测量、Three-Phase Fault三相故障、RLC Branch电阻电感电容支路、Powergui电力系统图形用户界面模块以及必要的Scope和To Workspace存储模块。模型连接顺序从电源开始三相电压源输出端接线路π型模型线路末端接负荷和测量模块测量模块后面接故障模块故障模块的输出端接地。这里有一个前提性操作必须从Simscape Electrical库拖一个Powergui模块放到模型里否则整个电气仿真无法运行。Powergui可选Continuous连续仿真或Discrete离散仿真本项目选Continuous模式、变步长求解器即可。线路π型模块的参数比较关键。Three-Phase PI Section Line需要填正序阻抗R和L、零序阻抗以及每相对地电容。如果觉得这个模块参数填写比较繁琐可以用三个RLC Branch并联在母线处代替集中对地电容效果一样且参数更直观三个支路分别接A、B、C三相到地每个支路的电容值填C0。3.2 故障模块与测量通道配置Three-Phase Fault模块是故障仿真的核心开关。双击模块进入参数设置在Fault Type中选择A相To Ground也就是A相接地故障。Switching Times填[0.05 0.25]意思是0.05秒时刻故障发生0.25秒时刻故障切除。Transition Time是故障过渡时间金属性接地时把这个值设成1e-6秒或者更小模拟瞬态开关动作如果设置得太大比如1ms波形上会看到非物理的渐变过程故障相的电压是缓慢下降而不是骤跌和实际金属性接地差距很大初学者经常在这个参数上翻车。测量配置方面Three-Phase V-I Measurement模块可以同时输出三相电压和三相电流接到Scope就能直接观察。为了后续做FFT分析或者画对比图建议同时接To Workspace模块把数据存到MATLAB工作区。不同信号的采样率可以在To Workspace里设置一般不需要太高20μs到50μs的采样周期足够覆盖暂态过程的主要频率。零序电压的提取有两种常见方法。一种是纯模块计算把三相电压测量信号送入Demux拆成三路单相信号用加法器求UAUBUC这就是3U0如果要求U0就再除以3。另一种是使用Sequence Analyzer序分量分析模块直接输出正序、负序、零序分量。前一种方法更透明能看到计算过程适合学习阶段用后一种方法更省事适合模型已经成熟以后用。3.3 从“不接地”切换到“消弧线圈接地”的关键操作中性点不接地系统建好之后切换成消弧线圈接地系统只需要在模型上做一处改动在电源中性点与地之间串入一个RLC Branch把电感值和电阻值填成0.61H和20Ω。问题是中性点在哪引出如果用的是Three-Phase Source模块参数对话框里有一个Neutral connection选项选择引出中性点后模块会多出一个中性点端子从这个端子接线即可。如果用的是Yg联结变压器中性点端子本来就在变压器图标上直接用Three-Phase Breaker连接后接地。这里特别容易出错的点是如果把变压器中性点直接接地而忘了串消弧线圈那就变成了中性点直接接地系统故障电流会是几千安培的短路电流整个仿真性质就完全变了。我实际操作中还尝试过一种更灵活的做法在消弧线圈支路上串联一个Three-Phase Breaker或者理想开关通过控制开关的合闸时间可以在同一个模型里先跑中性点不接地工况再在某个时刻合上消弧线圈支路观察系统切换的过程。不过如果目标是做两种情况的结果对比建议还是复制出两个独立模型分别仿真数据更干净分析也更省心。3.4 仿真求解器设置与运行时长规划仿真是否稳定和求解器设置关系很大。Simulink默认的ode45对简单的电气系统通常没问题但单相故障涉及开关动作信号存在刚性跳变我试过几次ode45在这种场景下容易发散或者仿真速度极慢。建议改用ode23tb刚性求解器并且把最大步长限制在1e-4秒以内。这个最大步长数值很关键如果设得太大比如0.01秒高频暂态分量完全采不到波形会严重失真设成1e-4秒左右既能捕捉到故障前后几十个周波的完整过程计算量也不会太大。仿真总时长设为0.3秒比较合适。前0.05秒系统从零状态起跑到稳态0.05秒时故障发生到0.25秒故障切除最后留0.05秒观察系统恢复到什么状态。这样一个时间轴安排稳态、暂态、故障持续时间、恢复过程全都能看清楚。4. 仿真结果对比与特征分析4.1 稳态电压特征两种系统殊途同归的部分把两种模型的仿真结果放在一起看首先注意到的是稳态电压特征几乎一致。金属性单相接地发生后A相电压跌落至接近零B相和C相电压升高到线电压量级也就是从5.77kV升到10kV左右。中性点电压从接近零的位置运动到负的相电压值也就是大约-5774V。这个现象可以从零序网络的角度解释。单相接地故障在零序等效电路中相当于在故障点接入了一个电压源这个电压源把整个系统的对地电压抬升了于是非故障相承受的电压从相电压变成线电压。这就是为什么中压配电网的设备绝缘水平要按线电压来考核——不是没有道理的单相接地时非故障相的过电压是实实在在的。消弧线圈对工频稳态电压的影响很小因为它改变的是电流回路而不是电压分布所以在电压波形上不接地系统和消弧线圈系统的稳态表现非常接近。如果只看电压波形很难判断系统到底是哪种接地方式。4.2 稳态电流走向容性主导与感性补偿的区别电流波形是两种系统差别最明显的地方。中性点不接地系统的故障相电流幅值大约在27安培左右和前面计算的理论值吻合相位超前电压90度标准容性电流特征。这个电流是由全网非故障相对地电容涌向故障点的所以它的值直接反映系统对地电容的大小。消弧线圈系统的故障点电流变化非常直观波形幅值从27安培骤降到2到3安培相位也改变了。这是因为感性电流滞后电压90度和超前90度的容性电流在故障点正好反相两者向量相加后基波分量几乎被抵消剩下的残流主要是有功分量和谐波分量。有功分量来自阻尼电阻谐波分量来自系统非线性因素。有一点需要特别注意消弧线圈补偿后故障电流并不是零。因为补偿度设定在110%留了一点电流裕度加上阻尼电阻带来的有功分量残流一般在2到5安培范围。这个数值虽然小但仍然存在选线保护需要在这个量级上想办法辨识故障线路难度确实不小。4.3 暂态过程的秘密故障瞬间的波形信息量稳态波形只能说明“故障发生以后系统变成什么样”而暂态过程告诉我们的却是“故障发生的那一瞬间系统经历了什么”。故障初瞬的暂态电流波形两种系统差异非常大。不接地系统在故障发生的第一个周波内高频分量很丰富非故障相的对地电容通过故障相迅速放电形成幅值高、衰减快的振荡电流振荡频率在几百赫兹到上千赫兹持续时间只有几个毫秒然后快速衰减到稳态电容电流值。这个高频分量是故障相电压骤降引起的本质上是电容的充放电过程。消弧线圈系统除了高频振荡之外还叠加了一个明显的衰减直流分量。原因在于电感电流不能突变故障瞬间消弧线圈中还没有电流它必须按照时间常数τL/Rd慢慢上升直到达到稳态补偿值。所以故障点电流从零开始爬升经过大约0.05到0.1秒才稳定到残流值。这个“从零爬升到残流”的过程就是电感储能逐渐建立的过程在波形上非常明显。这个暂态差异对保护原理有直接影响。基于稳态零序电流幅值和方向的选线方法在消弧线圈系统里经常失效因为残流太小了。实际工程中要利用暂态首半波的方向、零序电压变化率等特征来做选线这正是消弧线圈系统暂态量蕴含的工程价值。如果把消弧线圈的补偿度临时调到100%并且把阻尼电阻设得很小还能在仿真里观察到有趣的谐振现象故障消失后电压电流不会立刻衰减而是持续振荡很长时间。这就是谐振过电压的危险状态。仿真里观察到这个现象没有安全风险但能非常直观地理解为什么工程上必须采用过补偿并且加装阻尼电阻。5. 常见问题与调试经验速查5.1 仿真发散和波形毛刺的排查仿真发散是电力电子和电力系统仿真最常见的灾难。症状表现为波形突然出现巨大的尖峰数值跑到1e10量级仿真直接停止。排查方向有三个。第一是求解器不匹配。故障仿真包含开关动作属于刚性系统ode45效率低还容易崩。换成ode23tb或者ode15s以后稳定性提升非常明显。这是见效最快的一步。第二是模型参数刚度问题。如果线路电感和对地电容构成的谐振频率非常高而求解器的最大步长又不够小数值计算就会失稳。把最大步长限制到1e-4秒以内大部分毛刺问题都能解决。如果还是不行检查线路参数有没有量级错误比如把0.01Ω/km的电阻填成了0.01Ω/m。第三是故障模块的过渡时间。Transition Time设置过大故障过程本身就被拉长波形看起来像“缓坡”而不是“断崖”容易误判为系统特性。金属性接地必须把过渡时间设到微秒量级。5.2 新手最容易踩的坑第一个坑是电源模块的中性点设置。Three-Phase Source默认情况下中性点是悬空的很多新手以为已经实现了“中性点不接地”实际上是误打误撞。但当你想改造成消弧线圈系统时如果不知道怎么引出中性点就会卡住。一定要在电源参数对话框里找到Neutral connection选项正确引出后再接线。第二个坑是消弧线圈参数算反了方向。过补偿的意思是电感电流略大于电容电流IL Ic。有些人在脚本里把补偿度写成0.95变成欠补偿虽然也能运行但系统在一部分运行方式下有进入谐振的风险。仿真的结果和理论分析的边界条件就对不上了。第三个坑是故障时间太短。有人把故障持续时间设为0.02秒一个周波波形上看不到稳态特征全都是暂态过程拿这个波形去和理论公式对比自然对不上。故障时间窗口至少留0.1到0.2秒保证暂态衰减完成后还有足够的稳态数据。5.3 波形解读中的易混淆概念最后说几个波形解读时容易混淆的概念这些是我自己走过弯路的地方。零序电压不等于中性点电压。零序电压3U0是三相电压的向量和中性点电压是中性点对大地的位移电压。两者数值上有关联但物理含义不同报文和图表里不要把标签写混。残流不等于最小的故障电流。消弧线圈把基波容性电流补偿掉以后残流主要由有功分量和谐波分量构成它的幅值不是零相位也受阻尼电阻影响。看残流波形时不要期望看到一个幅值接近零的干净正弦波。单相接地故障和三相短路是两回事。三相短路不涉及零序网络中性点接地方式对它几乎没有影响单相接地才走零序通路才是我们研究的对象。在Three-Phase Fault模块里选择故障类型时一定要确认选的是A相To Ground而不是Three-Phase。为了排查方便把常见问题整理成速查表现象可能原因解决方法仿真中途发散求解器不合适或步长过大ode23tb 最大步长1e-4s故障相电压缓慢下降Transition Time设置过大设为1e-6s故障电流高达千A中性点直接接地而非经线圈检查消弧线圈支路是否串入故障后持续振荡不退补偿度接近100%且阻尼不足改为110%过补偿并增加阻尼电阻找不到中性点端子电源模块未导出中性点勾选Neutral connection选项暂态波形采不到高频求解器步长太大减小最大步长或降低采样周期把中性点不接地和消弧线圈接地两种系统从稳态到暂态完整跑一遍以后我最大的体会是教科书上那些一句话带过的概念在波形里全都有具体的模样。“电容电流超前电压90度”这句话看波形比背公式直观一百倍“消弧线圈的感性电流补偿容性电流”这件事看过残流从27安培降到3安培的瞬间你就再也不会忘掉。这个模型后续的扩展空间很大。可以在母线两侧挂多条馈线模拟真实配电网的选线环境可以设置不同过渡电阻对比金属性接地和弧光接地的特征差异还可以把零序电流互感器和选线算法用Simulink模块搭进去做一个闭环的故障选线测试平台。如果大家在搭建过程中遇到波形异常、参数调不通这类问题把模型里的线路参数和消弧线圈设置整理出来对照排查大部分问题都可以在几个常见检查点里找到答案。