干了这么多年配网故障分析我最大的感受是单相接地故障里真正让人头疼的从来不是金属性接地而是电弧性接地。你用固定电阻去模拟接地点出来的电压电流波形永远差口气一旦把电弧的动态特性加进去故障相电压恢复、非故障相过电压、零序分量这些现场录波里常见的现象才一个一个对得上。要把电弧行为与系统响应之间的关系彻底摸透绕不开ATP-EMTP这套老而弥坚的电磁暂态仿真工具。这篇博文就是从一次完整的“配电网单相接地电弧模型建模与仿真分析”项目里沉淀下来的从系统搭建、电弧模型实现到结果解读和避坑经验全部按实操顺序讲清楚。适合刚接触ATP的研究生、做配网故障分析的工程师以及准备做小电流接地选线仿真的同学参考。1. 为什么非要做这个仿真1.1 单相接地故障占了七八成还最难缠的故障配电网跟输电线路不一样中压系统绝大多数是中性点非有效接地包括不接地、经消弧线圈接地、经大电阻接地这几种方式。在这种系统里单相接地故障占到全部故障的70%到80%但故障电流并不大稳态值往往只有几安到几十安。这个特点决定了它“难发现、难处理、后果又很严重”。一个中压馈线发生单相接地规程允许继续带电运行一到两个小时目的是给运维人员留出寻线时间。但这只是针对稳定接地的说法一旦接地点是间歇性电弧情况就完全不同了。电弧反复熄弧、重燃会在健全相上堆积电荷产生几倍于相电压的暂态过电压既威胁电缆绝缘也容易把单相接地诱骗成两相短路。很多配网事故的录波最后追根溯源都是某个不起眼的电弧接地先烧起来的。所以对配电网做故障分析不能只算基波稳态必须看暂态过程。仿真研究要回答的问题也不是“接地后电压变成多少伏”这种简单问题而是故障相电压怎么恢复、非故障相过电压到底多高、零序电压电流里哪些特征量能拿来选线、消弧线圈补偿度对电弧自熄有什么影响。这些问题在现场很难大规模做试验只能靠电磁暂态仿真来支撑。1.2 电弧不能等效成固定电阻很多人第一次做单相接地仿真图省事直接在故障点放一个固定电阻比如10欧、100欧觉得接地阻抗差不多就得了。这个做法在算工频稳态时可以大致对付但一涉及暂态和电弧行为结果完全不能用。电弧的本质是气体放电通道它的电阻不是一个常量而是一个高度非线性的时变参数。电流大时弧道电离充分电阻可能只有几欧到几十欧电流过零附近电离减弱电阻迅速恢复到几十千欧甚至兆欧级。这个“从导通到截流”的过程直接决定了熄弧瞬间系统储能怎么释放、电压怎么恢复、会不会再次击穿重燃。你把这段动态过程压成一颗固定电阻等于把问题的核心部分直接删掉了。所以正确的路只有一条把电弧用电网模型能够接受的动态数学方程表达出来再耦合进ATP-EMTP的求解回路里。这也是“电弧行为与系统响应关系”研究的核心难点。电弧行为变了系统电压电流波形就变反过来系统参数变了电弧的燃弧时间、重燃次数也跟着变。两者是互相耦合的不是简单的先有鸡还是先有蛋而是必须放在同一个暂态求解框架里一起算。1.3 仿真要回答的工程问题清单动手之前我建议先把研究目标列清楚。目标是“过电压评估”还是“选线特征量提取”仿真方案是完全不同的。按我这次项目经验常见目标可以分成四类过电压水平评估统计非故障相暂态过电压峰值、弧光接地过电压倍数、不同故障相角下的最恶劣工况。消弧线圈参数验证调谐度、脱谐度对残流和熄弧效果的影响全补偿和过补偿的区别。故障选线特征分析零序电压、零序电流的幅值、相位、谐波含量暂态首半波方向。保护配合校验反时限保护、接地距离保护的动作边界间歇性接地时保护会不会误动。不同的目标对电弧模型的精度要求也不同。做选线特征重点看零序分量的相对关系和暂态主频电弧模型可以简化一些做过电压定量评估就必须把熄弧和重燃的细节模拟准。先明确要交什么结论再决定模型复杂度这个顺序不能反。2. 系统建模ATP里的10kV配电网怎么搭2.1 为什么用ATP-EMTP可能有人会问现在PSCAD、Matlab Simulink那么方便为什么还要用ATP这种界面老派的工具我的回答是ATP是老牌电磁暂态程序EMTP的免费替代实现1980年代开始流行到现在电力系统科研和工程里还在大量使用它不是被淘汰的老古董而是经过几十年工程检验的稳定工具。单说单相接地电弧仿真这件事ATP有几个优势是实打实的。第一它在开关动作、非线性元件、输电线路模型这些电磁暂态的看家领域处理得很扎实数值稳定性好第二模型自由度极高电弧这种非标准元件可以很方便地用TACS或MODELS写进去不强迫你用现成封装第三全网有大量文献和测试数据可以直接参考你做出来的结果容易跟同行对标。PSCAD的优势是集成度高、上手快但自定义微分方程模型的灵活性反倒不如ATP顺手。Matlab做后处理强但电磁暂态求解的底层细节需要自己操心仿真步长一大就发散。我整理过一个简单的工具选型对比可以参考对比项ATP-EMTPPSCAD/EMTDCMatlab/Simulink电磁暂态求解专业性强算法积累深厚强界面友好一般需自己搭架构自定义电弧模型非常灵活MODELS/TACS灵活可写自定义组件灵活但连接电网模型麻烦学习门槛偏高需要适应输入习惯中等中等文献参考丰富度高电弧仿真大量采用较高一般后处理较弱配合PlotsXY/外部脚本较强最强做配网单相接地电弧研究我的习惯始终是ATP-EMTP作为仿真核心波形数据导出后用Matlab或者Python做特征提取。前后端分离各用各的强项效率最高。2.2 系统拓扑与参数仿真系统不用搭得很复杂但也不能太简单。我这次用的是典型10kV配电网缩比模型一条10kV母线带三条馈线其中两条纯架空线一条架空电缆混合线路故障点设置在其中一条馈线的中段。电源按50Hz、线电压10.5kV考虑用内阻抗模拟系统等值短路容量一般按变电站母线短路容量500MVA来处理。线路参数采用典型10kV架空线参数我这里给出常用参考值参数架空线路典型值电缆线路典型值正序电阻 r10.27 Ω/km0.26 Ω/km正序电感 l1约1.1 mH/km约0.8 mH/km正序电容 c1约9.5 nF/km约250 nF/km零序电阻 r00.35 Ω/km0.75 Ω/km零序电感 l0约3.2 mH/km约3.4 mH/km零序电容 c0约6.0 nF/km约200 nF/km注意这些参数只是典型参考值具体工程仿真必须按实际导线型号和敷设方式查表或实测。不过对于机理研究和相对比较用典型值完全够用。负荷用恒阻抗模型每条馈线挂一个R-L并联支路有功和无功按实际线路容量比例设置。负荷建得过于精细对单相接地暂态影响不大没必要一开始就上潮流迭代。中性点接地方式直接在电源变压器中性点处理悬空、经消弧线圈接地、经小电阻接地三个方案都建出来方便后面对比。2.3 ATPdraw操作要点我一般习惯在ATPdraw里搭图搭完再生成ATP文件运行。几个关键操作点说一下第一电源和接地。电源用Type 14电压源模型设置幅值要算清楚线电压有效值10.5kV对应相电压峰值10.5kV×√2/√3≈8.57kV。频率50Hz内阻抗按短路容量反算。第二故障点。用Fault switch最方便可以直接指定故障相、接地方式和投入时间。这里并不直接接固定电阻而是预留一个节点对把故障相经过一个受控电阻元件接地。这个受控电阻就是后面MODELS电弧模型的输出接口。第三线路。架空线用LCC模型里的JMarti或Bergeron模型对地电容和相间耦合都能自动算出来。千万不要用纯串联R-L代替线路做暂态仿真那样对地电容通路就没有了零序特征完全出不来。第四仿真设置。步长和时长这个必须单独说。电弧时间常数在几十到几百微秒量级想捕捉到电弧电流过零附近的熄弧和重燃细节步长至少取1到5微秒我习惯取1微秒最大1微秒再大波形就失真了。仿真时长取0.2到0.4秒就够了前面留出系统稳定时间故障设置在0.04秒左右投入能看4到10个周期的燃弧过程。这里有一个很关键的习惯每次改步长先看同一工况下固定电阻接地的波形有没有变化。如果固定电阻工况的波形都随步长漂移那是系统建模有问题不是电弧模型的问题。3. 电弧模型建模从微分方程到ATP可执行模型3.1 电弧的两种经典数学描述电弧模型的建立核心思想不是去描述电弧内部每一个粒子的运动而是把电弧当成一个黑盒电阻用能量守恒关系描述它的动态特性。这个思路最早来自电弧的热平衡过程电弧等离子体从电场吸收能量同时向周围散发热量电导率的变化取决于能量平衡和热惯性。工程上最常用的两个模型是Cassie模型和Mayr模型。Cassie模型假设电弧通道半径随电流变化而电弧电压基本恒定适合描述燃弧阶段电流大、弧道电压变化不大的情况。Mayr模型则假设电弧通道半径不变用散热功率P来描述能量损失适合描述电流过零附近弧道温度剧变、电阻急剧变化的阶段。两者的微分方程形式非常相似区别只在于对电弧电压和散热功率的处理方式不同。Mayr模型方程的基本形式是(1/g) · dg/dt (1/τ) · ( (u·i)/P - 1 )其中g是电弧电导u是电弧电压i是电弧电流τ是电弧时间常数P是散热功率。方程的含义非常直观当输入功率u·i大于散热功率P时电导增大电弧加强当输入功率不足时电导减小电弧趋向熄灭。时间常数τ则决定这个变化的快慢τ越小电导变化越快电弧越容易不稳定。实际配电网单相接地研究中很少只用单一模型跑全程。更常见的做法是Mayr-Cassie组合模型燃弧期间电流大让Cassie模型起作用电流过零前后切换成Mayr模型描述熄弧到重燃的完整过程。切换逻辑一般按电弧电流的瞬时幅值自动判断电流大于某个阈值用Cassie低于阈值用Mayr并设置滞回避免在阈值附近来回跳。3.2 在ATP里写电弧模型ATP里没有现成的电弧元件需要借助MODELS语言自己写。MODELS是ATP内嵌的通用描述语言可以在每个仿真步长内求解微分方程并把计算结果输出给电网模型。我这次用的方案是电弧模型用MODELS求解输出动态电阻值接到电网故障点一个受控电阻元件上。这样只要改MODELS里的方程和参数就能自由切换Mayr、Cassie、组合模型不用动主电网图。MODELS伪代码大致是这个样子MODEL ARC_MAYR INPUT v_arc -- 电弧两端电压 OUTPUT r_arc -- 动态电弧电阻 VAR g, tau, P, i_arc INIT g : 1e-4 -- 初始电导给一个很小的值避免分母为零 tau : 100e-6 -- 时间常数初始取值100微秒 P : 5e3 -- 散热功率初始取值5千瓦 ENDINIT EXEC i_arc : v_arc * g IF ABS(i_arc) 1e-3 THEN dgdt : (1/tau) * ( (v_arc * i_arc / P) - 1 ) * g ELSE dgdt : 0 ENDIF g : MAX(g dgdt * tstep, 1e-6) r_arc : 1.0 / g ENDEXEC ENDMODEL这段代码里最值得注意的就是电导下限的处理。电弧电导出现负值或者零值是仿真发散的头号原因我习惯在每步更新后限制电导最小值比如1e-6西门子这样电弧熄灭时电阻不会变成无穷大导致数值故障而是一个很大的有限值。MODELS里tstep是当前仿真步长不需要自己显式取直接在积分更新里引用就可以。如果步长设置过大MODELS积分会不稳定电弧波形会出现异常锯齿这个问题后面第5节还会专门说。MODELS写好后在ATPdraw里把MODELS元件和受控电阻元件连接起来受控电阻的控制信号来自MODELS输出电阻本体接在故障相和地之间。这样电网求解器每走一步都会用上一时刻的MODELS输出电阻去计算电网状态再返回给MODELS更新电弧电阻实现两个系统耦合推进。3.3 参数怎么定电弧模型参数标定是电弧仿真里最主观、最容易被喷的地方。不同文献中τ和P的取值差异很大这不是模型不对而是电弧本身受电流等级、电极材料、弧道长度、介质条件影响太大了。配电网单相接地电流小几安到几十安与高压断路器开断大电流的电弧参数不在一个数量级不能直接抄大电流试验的τ值。我整理过一个配网接地电弧参数的初值参考范围参数取值范围说明时间常数 τ50~500 μs电流等级越小τ越大电弧越不容易熄灭散热功率 P0.5~5 kW与弧道长度、风冷条件相关初始电导 g01e-6~1e-4 S越小越接近冷态起始重燃阈值电压500 V~2 kV取决于间隙距离一般通过外施电压判断是否重燃这里必须说清楚上述参数只能作为仿真研究的起点不是“标准答案”。严格的做法是用实际现场录波数据或人工接地试验数据进行参数拟合以故障相电压恢复曲线和电弧电流过零时间作为目标反向调整τ和P。如果没有实测数据就至少在论文或报告里把参数依据写清楚再对τ和P做敏感性分析说明结论在参数偏差下依然稳定。我个人的经验是先把τ定在100微秒量级、P按故障电流峰值×电弧电压估算然后看电弧是否在一个周波内熄弧、是否出现合理的重燃现象。熄弧太快就加大P重燃过于频繁就减小P让波形先定性对上再定量微调。4. 仿真分析电弧行为与系统响应到底怎么互相影响4.1 金属性接地与电弧接地波形对比模型跑通之后第一件事不是急着调参数而是做一个对照同一系统、同一故障时刻分别用固定电阻接地和动态电弧接地做仿真把波形摆在一起看差异。固定电阻接地时故障相电压直接掉落到接近零非故障相电压稳态升到线电压整个过程干净利落电压波形是标准正弦波没有畸变。而动态电弧接地时故障相在电弧燃烧期间电压不是零而是呈现一个近似方波的小幅电压数值通常在几百伏到一千多伏之间这是弧道压降。更明显的是故障相电流正半波和负半波不对称电流过零附近不是快速线性过零而是出现明显的“零休”平台电弧在此短暂熄灭之后电压恢复并可能重燃。这个零休现象千万别当成数值误差它就是电弧非线性的直接表现。电弧在电流过零时输入功率为零但散热还在继续弧道温度快速下降电导率骤减电流无法立即跟随外电路电压变化。这段平台持续几十到几百微秒之后如果恢复电压高于介质击穿强度电弧又重新点燃。从系统侧看非故障相电压不再是规律的正弦波而是叠加了明显的谐波畸变尤其是3次和5次谐波。这些谐波就是后面选线算法里“5次谐波法”“谐波分量比相法”的基础。仿真如果不建电弧模型这些特征一个都不会出现。4.2 零序分量变化与接地方式的关系配电网单相接地最核心的系统响应是零序电压和零序电流。我分别对不接地系统和经消弧线圈接地系统做了故障仿真结果差异非常大。中性点不接地系统单相接地后零序电压3U0稳态幅值接近相电压的三倍并且因电弧非线性叠加了明显的谐波分量。零序电流主要是健全线路对地电容电流的叠加接地线路的零序电流方向与健全线路相反、幅值是全网电容电流之和这就是经典稳态选线的幅值方向比较依据。经消弧线圈接地系统情况完全不同。消弧线圈给接地故障点提供了一个感性电流补偿支路故障点电流从容性为主变成感性为主残流变小电弧更容易自然熄灭。仿真里能看到的是故障相电流波形峰值明显降低从几十安降到十安以内零序电压虽然仍然较高但谐波含量明显下降电弧重燃次数也减少。但消弧线圈接入也带来一个隐患如果补偿度接近全补偿即脱谐度接近零系统可能会出现串联谐振零序电压反而异常升高。我在仿真中把脱谐度从-5%扫到-15%清楚看到零序电压和故障点电流的谐振趋势。工程上坚持过补偿运行、脱谐度不要太小就是避免这个谐振风险。仿真结果和运行经验对上了模型的可信度就上来了。4.3 间歇性电弧重燃与弧光过电压电弧接地最危险的一面是间歇性重燃引起的弧光过电压。这个过电压的机理并不复杂电弧第一次熄灭瞬间故障相和非故障相上还有残余电荷系统对地电容上的电荷无法立刻泄放电弧重燃时相当于把故障相重新接地非故障相对地电容上的电荷被强制重新分配产生高频电荷转移电压出现强烈振荡。每一次重燃都是一次新的冲击如果重燃时刻刚好落在系统恢复电压峰值附近过电压会叠加得非常高。仿真中我故意把电弧参数设置为“易于重燃”的模式让故障相电流过零后不久就重新击穿。结果显示非故障相电压峰值可以达到正常运行相电压峰值的2.8到3.5倍个别恶劣相角条件甚至更高。这个倍数和教科书上弧光接地过电压的理论范围是一致的也是仿真结果可信度的有力佐证。过电压的严重程度跟电弧在电流过零后什么时刻重燃、外电路阻抗特性、故障相初始相位都有关系。我在分析中还做了一组故障相角扫描从0到360度每10度一个工况看不同起始相角下的过电压分布。结果发现故障发生时刻的相角对暂态过电压峰值影响很大最恶劣情况往往出现在相电压接近峰值附近起弧时这提醒做仿真研究的人只跑一个固定起弧相位就下结论是不严谨的必须做相角扫描找到包络。4.4 离脱谐度与故障距离的敏感性分析电弧行为和系统参数耦合得最明显的地方就是消弧线圈脱谐度。我把脱谐度依次设为0、-5%、-10%、-15%观察故障电流和零序电压变化发现脱谐度从过补偿方向增大时残流进一步减小电弧重燃次数减少非故障相过电压峰值也有下降。这个结论对运行部门选择补偿度很有指导意义也符合工程低残流运行的理念。故障距离的影响也值得单独说。同一馈线上故障点离母线越近健全线路对地电容放电回路越集中暂态电流幅值越大、频率越高故障点远离母线时暂态过程相对平缓但零序电流中各回线路的方向特征依然可辨识。从保护角度看近端故障特征突出但更容易饱和远端故障特征弱但是方向关系稳定选线判据的灵敏度和可靠性在这两种情形下需要折中评估。负荷水平对单相接地暂态的影响相对较小但也不是完全没有。重负荷时系统阻抗略微降低故障电流和暂态振荡幅值会有几个百分点的变化。如果研究目标是高精度定量评估建议做重载、轻载、空载三个工况的对比如果只是定性研究选线特征恒阻抗负荷模型加典型负荷率就够了。5. 常见问题与排查技巧实录5.1 仿真发散和不收敛这是电弧仿真里遇到最多的问题现象很吓人波形突然变成巨大毛刺要么飞到几千千伏要么直接NaN。我排查的顺序是固定的。先看步长。MODELS微分方程稳定性对步长极其敏感电弧时间常数100微秒步长如果取50微秒等于一个时间常数内只走两步积分基本废了。把步长降到1微秒大部分发散问题立刻消失。再看电弧电导限制确认代码里有没有给电导设置下限没有的话在过零瞬间电导可能冲成负值或者零导致整条支路阻抗非法。最后看受控电阻连接MODELS输出的是电阻而不是电导时注意1/g在g0时会除以零必须先限幅再求倒数。如果这些都没问题还发散把故障点上并联一个1纳法的小数值稳定电容试一下有时候只是LCC模型和时变电阻直接相连时开关切换带来的数值振荡不是电弧方程本身的问题。5.2 波形毛刺太多不光滑波形能跑完但是布满高频锯齿这种问题最容易让人误判成出现了真实的高频振荡。判断方法很简单把放电弧模型的仿真跑一遍用固定电阻接地看同样的位置有没有毛刺。如果固定电阻工况也毛刺那基本是系统模型问题比如线路模型选择不当、步长过大、或者Fault switch动作时刻与计算节点恰好冲突。具体解决办法有几个线路改用JMarti模型替代简单的Bergeron或π型模型JMarti在高频段拟合更平滑降低输出时采集数据的间隔用平均值滤波避免在元件参数突变的同时又让MODELS电阻快速变化把故障投入时刻和电弧起弧时刻分开比如0.04秒接入故障MODELS初始电导从很小值开始逐渐增大而不是直接突变。5.3 电弧不熄弧或者瞬间就熄灭这个问题几乎每次调参数都会遇到。电弧一直燃着不熄在波形上表现为故障相电流始终连续、过零后没有零休平台稍微一调参数又变成起弧一个周波后直接熄灭后面全是开路状态根本没有间歇性重燃过程。本质是τ和P不匹配。电弧不熄P设置过大散热功率远大于电弧输入功率弧道始终冷却不下来或者τ设置过大电导变化太慢跟不上电流过零。先减小P一个数量级试试。电弧瞬间熄灭P过小电流过零后弧道温度衰减太快恢复电压根本撑不起重燃。把P调回去并把τ适当增大给弧道增加一点热惯性。这里有个从录波里总结出来的经验真实配电网电弧很少是固定周期重燃重燃时刻与外电路恢复电压和弧道热状态联合决定所以仿真结果里重燃间隔不整齐反而正常。如果重燃严格按固定周期出现多半是模型里某个逻辑写成了定时刻开关不是真实的电弧行为。5.4 仿真结果和现场录波对不上仿真波形和现场故障录波对不上是电弧研究中不可避免的一步关键在于知道哪些偏差可以接受、哪些不能接受。我的判断标准是电弧电流峰值和过零时间的量级必须对得上故障相电压恢复过程的波形形态必须对得上但如果高频振荡的幅值和频率有几个百分点的差异可以接受因为现场系统参数本身就不完全清楚线路对地电容尤其难测准。调参顺序建议先对工频稳态调系统等值阻抗和线路电容把零序电压稳态幅值对齐再对暂态过零过程调电弧τ和P把电流过零附近的零休持续时间和重燃间隔对齐最后才对高频振荡调线路模型和杂散参数。反过来调会非常痛苦因为绕来绕去不知道是哪个参数在起作用。参数敏感性做扎实之后结论即使不能精确复现某一条现场波形也能保证它的趋势和相对关系是可信的。这一点在写论文和做工程结论时很有用评审专家问起来你可以明确说出哪些结论受参数影响大、哪些不受影响。结束前的一点个人体会这类电弧仿真做多了我越来越觉得最关键的不是模型公式多高级而是耐心做对比和校验。先和固定电阻工况对比证明电弧模型在起作用再和理论过电压范围对比证明参数合理最后和现场录波对比证明结论可信。每一步都留好波形和数据后面分析才有底气。下一步如果有精力可以把这套电弧模型应用到小电流接地选线算法验证里把零序电压、零序电流的特征量提取出来做量化评估甚至可以叠加分布式光伏和电缆混接场景看看电弧接地特性在新型配网里又会出什么新花样。仿真的边界拓展空间还很大但基础一定要先扎实。