Mayr与Cassie电弧模型：配电网故障仿真双引擎解析

电弧仿真这件事干电力系统的人迟早要碰。不管是算断路器开断能力、研究故障弧光接地还是做继电保护装置的动作逻辑验证电弧模型都是绕不过去的坎。而只要一提到电弧模型Mayr和Cassie这两个名字几乎永远是并列出现的。它俩一个管大电流燃弧一个管小电流过零搭配起来几乎覆盖了工程上90%以上的电弧仿真场景。这篇文章我不打算念教科书。我从模型推导的物理直觉讲起结合我自己搭过的真实仿真实例把两个模型的脾气、适用范围、参数坑、和工程选型逻辑一次说清楚。内容主要面向做配网仿真、保护装置研发、或者刚接触电磁暂态仿真的工程师当然学校里做科研的朋友想快速建立概念也能直接用得上。1. 两个模型到底在描述什么从能量平衡到黑盒参数很多人把Mayr和Cassie当成两个拿来就用的现成封装模块参数照着论文抄一遍就开始跑仿真。这样做的结果就是波形看起来挺像那么回事但换一个工况马上就不准了。问题的根源在于我们根本不理解这两个模型背后的物理假设是什么。1.1 Mayr模型热惯性决定的小电流电弧行为Mayr模型的出发点是能量平衡。电弧通道本质上是一个气体放电通道电流通过时焦耳发热同时通道向周围散热。电弧是否维持取决于热量的净积累。Mayr假设电弧通道是一个半径固定、能量均匀分布的圆柱体散热功率近似恒定于是可以推导出电弧电导的动态方程[ \frac{1}{g} \frac{dg}{dt} \frac{1}{\tau} \left( \frac{u i}{P_0} - 1 \right) ]这里 ( g ) 是电弧电导( u ) 是弧压( i ) 是弧流( P_0 ) 是电弧散热功率常数( \tau ) 是时间常数。注意这个方程的核心逻辑当注入功率 ( u i ) 大于散热功率 ( P_0 ) 时电导上升电弧变强反过来电导衰减电弧趋向熄灭。整个模型管的是“电弧通道的热动态过程”它天生适合描述电流过零附近、电弧接近熄灭和重燃临界状态的动态行为。之所以说它适合小电流是因为Mayr模型假定电弧电压由热游离主导散热量恒定。当电流很小、弧柱截面收缩明显时这一假设和真实电弧的物理过程非常贴近。1.2 Cassie模型对流换主导的电压梯度恒定假设Cassie模型的思路完全不同。它描述的是大电流燃弧阶段此时电弧通道截面大、温度高弧柱中单位体积储存的能量不是主导因素真正重要的是气流或电磁力作用下的对流换能力。Cassie假设电弧电压梯度是常数即弧压不随电流变化只与弧长有关。方程形式如下[ \frac{1}{g} \frac{dg}{dt} \frac{1}{\tau} \left( \frac{u^2}{E_0^2} - 1 \right) ]其中 ( E_0 ) 是静态电弧电压梯度。这个模型把电弧视作一个可变的电阻通道通道截面积随电流大小伸缩而单位长度弧柱的电压基本不变。实际物理图像就是大电流时弧柱粗小电流时弧柱细可弧柱单位长度上维持导电所需的电场强度基本恒定。基于这个假设Cassie模型能比较准地模拟短路电流大燃弧阶段的弧压波形弧压表现为接近方波状、幅值基本和电流幅度解耦的特征这与高压断路器开断试验中录到的弧压波形相当吻合。1.3 一个表看懂两者定位差异对比维度Mayr模型Cassie模型物理假设通道半径固定散热近似恒定弧柱电压梯度恒定主导物理过程热游离与热惯性对流换能力和通道截面积变化适用电流范围小电流近零区几十安培以下大电流数十至数十千安典型应用弧光接地、消弧线圈、小电流选线断路器开断、短路燃弧、保护配合弧压波形特征电流过零时弧压尖峰明显燃弧区弧压接近方波参数敏感度对散热常数和时间常数敏感对电压梯度和时间常数敏感这个表基本就是选型时的第一层过滤网。拿到一个仿真任务先判断关注的是燃弧大电流阶段还是熄弧临界阶段模型就定了一半。2. 仿真实例设计为什么组合使用是工程常态单纯把一个模型拿出来单独用往往达不到工程精度。因为真实的电弧过程是连续演化的短路初期电流大Cassie模型好使电流过零后进入弧后临界区热惯性和散热变成核心Mayr模型更贴近实际。所以工程上最常见的做法是用分段切换或加权组合的方式把两个模型拼成一个“复合电弧模型”。这也解释了为什么很多商业软件比如EMTP-ATP、PSCAD/EMTDC里的标准电弧模块内部结构都是MayrCassie的组合形式。我在一个配电网单相接地故障仿真实例中实测过两套方案一套只用Mayr另一套用MayrCassie组合。故障点设在10kV中性点不接地系统的A相过渡电阻为0.01欧姆模拟的是金属性弧光接地。整个仿真时长0.2秒故障在0.05秒时投入。单用Mayr模型时燃弧阶段的弧压峰值被明显高估因为Mayr模型默认散热功率恒定大电流下发热功率远超散热电导增长速度被数学限制卡住导致弧压畸变。而改用组合模型后大电流燃弧段由Cassie主导弧压被限制在一个合理的梯度范围过零附近再切换给Mayr模型处理弧压尖峰和熄弧重燃的过渡过程都更接近实测录波。2.1 仿真模型搭建参数选择是第一个坑参数选择直接影响结果的真实性。Mayr模型中两个关键参数时间常数 ( \tau_M ) 和散热功率 ( P_0 )。工程上( \tau_M ) 一般在0.5到2微秒之间取( P_0 ) 需要结合系统短路电流水平和预估弧压来反推。经验公式大致是[ P_0 \approx 0.4 \times U_{arc} \times I_{arc(max)} ]即散热功率约为弧压和最大燃弧电流乘积的百分之四十。当然这只是初值后面要走参数整定。Cassie模型的时间常数 ( \tau_C ) 一般取得比Mayr大典型在5到20微秒之间反映对流换的过程比导热慢。弧柱电压梯度 ( E_0 ) 则和介质、气压、弧长有关。空气中的一个标准大气压下经验值在每米10到20伏每厘米也就是1000到2000伏每米。如果仿真的是SF6断路器这个梯度会高很多可以到每米上万伏的量级。参数设置的顺序也讲究。先根据故障电流幅值定 ( E_0 )再根据燃弧时间内的弧压强度定 ( P_0 )最后根据观测到的弧压过零尖峰宽度反推 ( \tau_M ) 和 ( \tau_C )。整个过程是“先宏观后微观”一上来就抠时间常数往往越调越乱。2.2 分段切换逻辑什么时候该听谁的在组合模型里切换逻辑是灵魂。最简单的做法是根据电流瞬时值判断当 ( |i| I_{threshold} ) 时启用Cassie低于阈值时启用Mayr。阈值一般是额定电流的10%到20%。但我实测下来这个硬切换在阈值边界处容易造成电导跳变引发数值振荡。更平滑的做法是使用加权函数过渡[ \frac{1}{g_{eq}} \frac{\alpha}{g_{Cassie}} \frac{1-\alpha}{g_{Mayr}} ]其中 ( \alpha ) 是与电流幅值相关的平滑权重参数比如[ \alpha \frac{1}{1 e^{-k(|i|-I_{th})}} ]这样在阈值附近有一个连续的过渡带而不是突然从一个模型切到另一个模型。( k ) 控制过渡带宽度我一般取 ( k I_{th}/10 )效果比较理想电导变化连续波形没有毛刺。2.3 我实测的波形特征记录用组合模型跑出来的故障相电压波形有几个典型特征。燃弧期间弧压稳定在一个较低水平相对相电压过零时弧压出现明显的尖峰突起尖峰幅值可以达到燃弧电压的数倍这对应着电弧熄灭后介质强度快速恢复的过程。故障相电压不再是简单接地零电压而是在过零附近呈现出明显的“双峰”结构这一细节对保护装置的零序电压门槛整定影响很大很多保护误动案例仿真阶段没把这个细节做出来导致装置整定值偏差。3. 实操过程一个完整的配电网弧光接地仿真计算下面我拆解一个完整的仿真工程案例用的平台是MATLAB/Simulink配合自编程的电弧动态模块。这个方法的好处是完全透明不会黑箱每条公式和参数都摊开在面前。工业上当然可以用PSCAD现成模块但对理解模型本质来讲手写一遍的价值不可替代。3.1 系统参数与仿真场景设定以一个10kV中性点不接地配电系统为背景。电源经降压变压器提供10kV母线电压馈线长度设为15公里线路参数采用架空线典型值正序阻抗 ( R_1 0.27\Omega/km )( L_1 1.28mH/km )零序参数酌情放大。在馈线末端设置一个单相接地故障故障点前方加装电弧模型。仿真总时间0.3秒故障触发时刻0.08秒持续到0.15秒后由保护动作切除。步长设定为1微秒。这里一定要说清楚电弧模型方程的时间常数在微秒级仿真步长若取太大比如50微秒以上电弧动态过程会被严重平滑掉过零尖峰会失真保护逻辑验证就失去了意义。3.2 核心代码状态方程与仿真流程电弧模型在Simulink里通过S-Function实现。下面给出一段核心残差计算代码展示Mayr和Cassie模型的嵌入方式MATLAB语法function [sys,x0,str,ts] arc_model(t,x,u,flag,param) % u(1): 电弧电流 i % u(2): 弧长或控制信号 % x(1): ln(g)取对数是为了保证电导始终为正 switch flag case 0 [sys,x0,str,ts] mdlInitializeSizes(); case 3 sys mdlOutputs(t,x,u,param); case 1 sys mdlDerivatives(t,x,u,param); case {2,4,9} sys []; otherwise error([Unhandled flag ,num2str(flag)]); end function sys mdlDerivatives(t,x,u,param) g exp(x(1)); i u(1); % 根据电流幅值计算加权系数alpha abs_i abs(i); alpha 1 / (1 exp(-param.k * (abs_i - param.Ith))); % Mayr模型电流 dMayr (1/param.tauM) * ((param.Uarc * i / param.P0) - 1); % Cassie模型电导导数 dCassie (1/param.tauC) * ((param.Uarc^2 / param.E0^2) - 1); % 等效电导导数 加权和 dlnG alpha * dCassie (1 - alpha) * dMayr; sys dlnG; function sys mdlOutputs(t,x,u,param) g exp(x(1)); u_arc u(1) / g; % 弧压 sys [g; u_arc];这里有个关键设计状态变量选的是 ( \ln(g) ) 而不是 ( g ) 本身。原因是电导值跨越几个数量级直接用 ( g ) 做状态变量在过零附近容易出现负值导致数值崩溃。对数化处理后无论电导多小状态变量都有明确的数学意义积分过程中不会因为舍入误差跌破物理下限。这个小技巧在多次工程实践中都很管用。仿真流程方面每个步长内电源、线路、电弧模块联立求解。电弧模块接收上一时刻的电流输出电弧电导和弧压再回传给电网方程。接口上用希尔伯特变换做延时校正避免单步长数据错位这个细节大家正式做课题时需要注意。3.3 参数整定过程记录首轮参数取值Mayr时间常数0.8微秒散热功率 ( P_0 15000W )Cassie时间常数8微秒( E_0 1500V/m )阈值电流50A。跑完一轮后发现燃弧区弧压偏低。对照录波数据弧压峰值应该在700到900伏区间实际只有500伏。检查发现是 ( E_0 ) 取值偏小把弧长设为0.8米时 ( E_0 ) 应该在2200伏每米附近。调整后再跑弧压抬升到850伏波形形状和实测对上了。第二处调整是Mayr模型中的 ( P_0 )。初始值取值偏大导致过零后电弧重燃不积极叠加上 ( \tau_C ) 影响后重燃延迟了大约80微秒。把 ( P_0 ) 从15000降到8500瓦重燃时刻回摆到范围内故障相的电压过零尖峰宽度也从原来的160微秒压缩到90微秒接近实测值。3.4 波形数据解读提取保护相关的关键特征量仿真跑完后不能只盯着波形看个大概。要输出电弧电压峰值、燃弧持续时间、重燃时刻、弧流过零前后的 ( di/dt ) 等特征参数。这些量直接关系到保护装置的整定。比如我在这轮仿真里提取到燃弧期间弧压峰值862伏是相电压峰值8165伏的大约10.5%表现为非接地故障特征过零后首个电压尖峰幅值2400伏持续时间约90微秒故障电流从50安培向上重燃的速率约3.2安培每微秒。这几个参数拿去跟保护装置的零序电压启动门槛做匹配时就能发现常规固定值整定法在这类高阻接地场景会漏判动态门槛或时频特征提取才更有优势。这就是电弧模型仿真在保护研发中的实际价值。4. 调试与问题排查电弧仿真里我踩过的五个坑长时间调电弧模型总能攒下一些刻骨铭心的教训。这里集中整理五个高频问题每一条都是实际跑仿真时真实遇到过的。4.1 数值刚性步长小于时间常数的必要性电弧模型的时间常数在微秒级而电网机械暂态往往在毫秒到百毫秒级两者天然形成刚性系统。最开始我试过用变步长求解器自动调节结果发现电弧过零附近步长缩到纳秒级而大电流燃弧段步长放宽到几十微秒来回切换导致收敛性极差。后来学乖了固定步长直接设为最大时间常数的三分之一以下。比如 ( \tau_M 0.8\mu s )步长就定在0.2微秒稳得很。代价是仿真时间拉长但至少结果可信。4.2 时间常数的量级错位判断收敛性的第一道关判断电弧仿真有没有跑疯第一眼看时间常数。如果Mayr模型的 ( \tau ) 被整到了几十微秒那电弧过零行为早就变味了因为真实燃气电弧的热时间常数就是亚微秒到几个微秒的量级。Cassie的 ( \tau ) 相对可以大一些但也不能超过50微秒量级太多。超了就说明算法在掩盖数值问题而不是在模拟物理过程。这个检查原则在论文和工程报告里都适用审稿人和领导看到这个数据都会认可的。4.3 参数突变后的高频振荡处理组合模型在切换点附近因为电导加权过渡有时候会在切换点附近激发高频数值振荡频率可达几十千赫。处理手段是设置一个限幅滤波器对输出的弧压做一阶低通截止频率设为燃弧电压基频的10到20倍。但注意这个滤波器会轻微影响过零尖峰幅值低通要轻用。我在调试末期把滤波截止频率从80kHz调到160kHz才在波光滑度和过零特征保真之间找到平衡。4.4 模型参数与线路分布参数的匹配关系这个坑比较隐蔽。电弧模型的弧压工况依赖于实际流过电弧的电流而分布参数线路模型在高频下会改变故障点的电流幅值。如果你用的是集中参数线路模型故障电流偏大弧压会被压到偏低用分布参数模型则更贴近实际。所以参数整定前先确认线路模型的选择。不要套用论文里的电弧参数时忽略对方线路建模方式。4.5 仿真时长陷阱把暂态跑成稳态的假象有的工程师为了省时间把故障持续时间缩短到几个周期只看前几个周期的波形。这样做在大燃弧阶段问题不大但涉及过零重燃的统计特性时就完全劝退了。电弧过零行为天然存在分散性没有一个电弧是每次过零行为完全一样的。至少跑20个工频周期以上统计重燃时刻和电压尖峰幅值的分布范围得到的中位数和极值才有工程参考意义。5. 从电弧模型到数字化系统与“电力系统CDT规约解析”的关联思考不少做配网自动化或变电站监控的工程师看到电弧模型会觉得这是“电磁暂态仿真圈”的事跟自己的监控后台关系不大。但电弧故障特征的获取最终要落到保护动作和监控告警上。比如弧光接地产生的高频暂态信号经过电压互感器传到保护装置或故障录波器由监控系统从前端设备汇集数据并上送调度。这里就绕不开“电力系统CDT规约解析”这类传统而关键的通信协议工作。CDT循环式远动规约在配网终端中依然大量存在尤其是早期投运的厂站。从电弧模型仿真中得到故障波形通过电力系统CDT规约解析把故障时刻的暂态特征量电压有效值、零序电压、故障电流量化打包上送是连接“仿真分析”和“现场监控”的一条现实路径。我在做仿真时会有意识地把电弧过零区间的特征量按规约格式整理为后续接入现场监控系统做准备。这里建议做相关工作的朋友别只盯着波形图乐呵将仿真输出与传统远动规约的字段设计匹配起来能让研究成果更快落地。另一个关联方向是电力系统自动化模板的标准化。电弧仿真的临时脚本越来越多代码风格千奇百怪。如果团队里有标准化模板把电弧模型参数的命名规范、输出变量格式、仿真报告框架做成模板对项目交付和审计都大有帮助。我自己的做法是做一个统一的电弧仿真结果清单模板包含故障条件、模型参数、关键波形截图、特征量汇总四大部分任何人接手都能快速复现。最后再说点实在的我自己的体会是搞电弧仿真的核心其实是给物理过程选择合适的“数学代理人”。Mayr和Cassie都不是全能的甚至在今天的标准里看有点粗糙但它们是经典模型体系里最经得起工程检验的那一档。用得好不好取决于你是否了解它们各自守护的物理边界。如果你正准备做一个电弧仿真项目别急着写代码先回答三个问题仿真关注的是燃弧大电流段还是过零熄弧段故障点的电流幅值区间覆盖什么范围实测录波数据有没有能不能用来校准参数这三个问题想清楚了选择Mayr还是Cassie、或者怎么组合自然会有答案。如果这篇文章里的某个坑正好帮你少走了一段弯路或者你有更离奇的调试经历想交换欢迎在评论区聊。毕竟电弧这玩意儿理论上有边界工程上总是出其不意。