企业官网建设流程全解析

干了这么多年配网故障分析，我最大的感受是：单相接地故障里真正让人头疼的从来不是金属性接地，而是电弧性接地。你用固定电阻去模拟接地点，出来的电压电流波形永远差口气；一旦把电弧的动态特性加进去，故障相电压恢复、非故障相过电压、零序分量这些现场录波里常见的现象，才一个一个对得上。要把电弧行为与系统响应之间的关系彻底摸透，绕不开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-EMTP作为仿真核心，波形数据导出后用Matlab或者Python做特征提取。前后端分离，各用各的强项，效率最高。

2.2 系统拓扑与参数

仿真系统不用搭得很复杂，但也不能太简单。我这次用的是典型10kV配电网缩比模型：一条10kV母线带三条馈线，其中两条纯架空线，一条架空电缆混合线路，故障点设置在其中一条馈线的中段。

电源按50Hz、线电压10.5kV考虑，用内阻抗模拟系统等值短路容量，一般按变电站母线短路容量500MVA来处理。线路参数采用典型10kV架空线参数，我这里给出常用参考值：

注意这些参数只是典型参考值，具体工程仿真必须按实际导线型号和敷设方式查表或实测。不过对于机理研究和相对比较，用典型值完全够用。负荷用恒阻抗模型，每条馈线挂一个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的取值差异很大，这不是模型不对，而是电弧本身受电流等级、电极材料、弧道长度、介质条件影响太大了。配电网单相接地电流小（几安到几十安），与高压断路器开断大电流的电弧参数不在一个数量级，不能直接抄大电流试验的τ值。

我整理过一个配网接地电弧参数的初值参考范围：

这里必须说清楚，上述参数只能作为仿真研究的起点，不是“标准答案”。严格的做法是用实际现场录波数据或人工接地试验数据进行参数拟合，以故障相电压恢复曲线和电弧电流过零时间作为目标，反向调整τ和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在g=0时会除以零，必须先限幅再求倒数。

如果这些都没问题还发散，把故障点上并联一个1纳法的小数值稳定电容试一下，有时候只是LCC模型和时变电阻直接相连时开关切换带来的数值振荡，不是电弧方程本身的问题。

5.2 波形毛刺太多，不光滑

波形能跑完，但是布满高频锯齿，这种问题最容易让人误判成出现了真实的高频振荡。判断方法很简单：把放电弧模型的仿真跑一遍，用固定电阻接地看同样的位置有没有毛刺。如果固定电阻工况也毛刺，那基本是系统模型问题，比如线路模型选择不当、步长过大、或者Fault switch动作时刻与计算节点恰好冲突。

具体解决办法有几个：线路改用JMarti模型替代简单的Bergeron或π型模型，JMarti在高频段拟合更平滑；降低输出时采集数据的间隔，用平均值滤波；避免在元件参数突变的同时又让MODELS电阻快速变化，把故障投入时刻和电弧起弧时刻分开，比如0.04秒接入故障，MODELS初始电导从很小值开始逐渐增大而不是直接突变。

5.3 电弧不熄弧或者瞬间就熄灭

这个问题几乎每次调参数都会遇到。电弧一直燃着不熄，在波形上表现为故障相电流始终连续、过零后没有零休平台，稍微一调参数又变成起弧一个周波后直接熄灭，后面全是开路状态，根本没有间歇性重燃过程。本质是τ和P不匹配。

电弧不熄：P设置过大，散热功率远大于电弧输入功率，弧道始终冷却不下来；或者τ设置过大，电导变化太慢，跟不上电流过零。先减小P一个数量级试试。电弧瞬间熄灭：P过小，电流过零后弧道温度衰减太快，恢复电压根本撑不起重燃。把P调回去，并把τ适当增大，给弧道增加一点热惯性。

这里有个从录波里总结出来的经验：真实配电网电弧很少是固定周期重燃，重燃时刻与外电路恢复电压和弧道热状态联合决定，所以仿真结果里重燃间隔不整齐反而正常。如果重燃严格按固定周期出现，多半是模型里某个逻辑写成了定时刻开关，不是真实的电弧行为。

5.4 仿真结果和现场录波对不上

仿真波形和现场故障录波对不上，是电弧研究中不可避免的一步，关键在于知道哪些偏差可以接受、哪些不能接受。我的判断标准是：电弧电流峰值和过零时间的量级必须对得上；故障相电压恢复过程的波形形态必须对得上；但如果高频振荡的幅值和频率有几个百分点的差异，可以接受，因为现场系统参数本身就不完全清楚，线路对地电容尤其难测准。

调参顺序建议：先对工频稳态，调系统等值阻抗和线路电容，把零序电压稳态幅值对齐；再对暂态过零过程，调电弧τ和P，把电流过零附近的零休持续时间和重燃间隔对齐；最后才对高频振荡，调线路模型和杂散参数。反过来调会非常痛苦，因为绕来绕去不知道是哪个参数在起作用。

参数敏感性做扎实之后，结论即使不能精确复现某一条现场波形，也能保证它的趋势和相对关系是可信的。这一点在写论文和做工程结论时很有用，评审专家问起来，你可以明确说出哪些结论受参数影响大、哪些不受影响。

结束前的一点个人体会

这类电弧仿真做多了，我越来越觉得最关键的不是模型公式多高级，而是耐心做对比和校验。先和固定电阻工况对比证明电弧模型在起作用，再和理论过电压范围对比证明参数合理，最后和现场录波对比证明结论可信。每一步都留好波形和数据，后面分析才有底气。下一步如果有精力，可以把这套电弧模型应用到小电流接地选线算法验证里，把零序电压、零序电流的特征量提取出来做量化评估，甚至可以叠加分布式光伏和电缆混接场景，看看电弧接地特性在新型配网里又会出什么新花样。仿真的边界拓展空间还很大，但基础一定要先扎实。