基于Matlab的电压依赖电晕效应输电线建模与电磁暂态仿真
1. 项目概述1.1 核心需求解析做电力系统仿真的人尤其是研究高压输电线路特性的工程师大概率都绕不开一个东西——电晕效应。我这次在Matlab环境下把一个带电压依赖特性的电晕效应输电线模型完整实现了出来模型是基于电磁暂态仿真思路搭建的重点刻画了电压变化如何直接影响电晕损耗、导线对地电容修正以及波过程传播特性。这个模型能解决什么问题一句话说清楚传统输电线路仿真里要么干脆忽略电晕要么把电晕损耗当成固定常数叠加进去这样在正常电压下误差还能忍受但一旦遇到过电压、操作暂态、谐振抬升这类电压剧烈波动的场景固定参数的模型就会严重失真。电压依赖型电晕模型就是让线路参数跟随瞬时电压动态调整把电晕放电那个非线性过程用工程上够用的精度表达出来。适合谁来参考两类人最对口一是做电力系统过电压分析与绝缘配合的工程师需要在仿真平台里评估操作过电压、雷电过电压下线路的衰减特性二是在校研究生或做输电线路在线监测算法验证的开发者需要一套可复现、可修改的电压依赖线路模型作为算法测试床。我自己在模型调通之后拿它跑了几组典型的电压波动工况对比恒定参数模型波前陡度、峰值衰减这些指标的差异肉眼可见可以说这个模型的价值就在于把“电晕会让电压波形变形”这个经验性认知变成可量化的仿真结果。底层思路其实不复杂电晕放电本质上是一个场强阈值触发的非线性过程导线表面电场强度低于起始场强时线路表现接近理想一旦超过起始场强电晕电流开始出现等效来看就是导线周围空气被局部击穿形成导电通道宏观表现就是线路出现额外损耗、等效半径增大电气上相当于导线变粗、对地电容变化。这三个宏观表现全部跟电压瞬时值强相关因此模型的核心就是构建电压到这些参数的映射关系。1.2 模型适用场景电压依赖电晕模型的应用面比想象中宽我自己梳理下来至少有这些典型场景值得用操作过电压仿真合空线、切空线、重合闸过程中暂态电压幅值很容易冲到额定值的2倍以上这个区间电晕损耗对波头削平作用非常明显不考虑电晕的话过电压峰值计算结果会偏保守绝缘子串选型就可能过度设计。雷电过电压分析雷电波侵入变电站时波前极陡电晕效应会导致波在传播过程中明显变形衰减电压依赖模型能让雷电过电压的沿线分布计算贴合实测规律。长距离输电线路的工频过电压研究空载长线容升效应让末端电压升高电晕损耗相当于一个随电压变化的负荷影响容升幅值。电磁暂态程序算法验证自己写FDTD、时域有限元或者改进贝杰龙法的时候需要有一个带非线性环节的线路模型来验证算法处理非线性边界的能力。还需要提醒一点这个模型面向的是宏观电磁暂态仿真维度不是用来做导线表面电场微观分析或者放电等离子体机理研究的。如果要做电晕放电的物理机制研究那需要PIC仿真或者流体模型那完全是另一个世界观的东西。我们这个模型定位是在保证暂态精度可接受的前提下把电晕效应的宏观电压依赖特性嵌入到输电线模型里让系统级仿真更接近真实。2. 电晕效应的物理基础与建模思路2.1 电晕放电的物理图景电晕放电说穿了就是导线周围空气在高场强下发生的局部击穿。导线表面是不光滑的加上污秽、水滴这些因素局部电场会远高于平均场强。当导线表面电场强度超过空气的击穿阈值工程上常取20~30kV/cm跟气压、温度、湿度都有关时导线附近极薄一层空气里的自由电子被电场加速撞击中性分子产生新的电子和正离子——这就是电子雪崩大量雪崩叠加就形成了自持放电。注意一个关键物理事实电晕放电发生在导线附近的强场区不会像火花放电那样贯穿整个空气间隙。所以电晕本质上是一种“局部自持放电”它不导致闪络但会持续消耗能量、产生臭氧和氮氧化物、发出可听噪声和无线电干扰。对输电线路仿真而言真正的归属是损耗模型和波形畸变模型。电晕损耗跟很多因素纠缠在一起导线半径、导线表面状况、大气条件、电压幅值和频率其中电压的影响最直接也最根本——因为电晕起始与否完全由电压决定的表面场强说了算。正半周和负半周的电晕特性还有差异负极性电晕起始电压通常更低这就是后面模型里需要区分电压极性至少对直流线路必须区分对交流系统可以按有效值走简化路径的原因。工程上估算导线表面场强常用皮克公式。皮克公式的典型形式为E₀ 30.3 × m × δ × (1 0.298 / √(r_eff × δ))单位kV/cm其中m是导线表面粗糙系数光滑导线取1但实际工程导线取0.8左右因为绞线表面不光滑δ是相对空气密度r_eff是导线等效半径单位cm。从公式就能看出场强阈值跟半径成反比趋势——细导线更容易发生电晕这也符合物理直觉。那么有了起始场强阈值电压依赖模型的思路就清晰了正常运行电压下导线表面场强低于阈值线路无电晕或者电晕极弱线路参数接近恒定当电压升高导致场强超过阈值进入电晕区这时候就应该计算出一个随电压变化的电晕电导或者视在损耗功率并联到线路模型里。这里“电压依赖”四个字的含义就落地了。2.2 电晕电压依赖特性的数学刻画要让模型可计算必须把电晕电压依赖特性用具体数学形式表达出来。工程上积累了大量电晕损耗的经验公式常见的有这么几类幂函数型P_c K × (U / U₀)^n其中n通常落在3~5之间K是比例系数U是实际电压有效值U₀是电晕起始电压。这个形式简洁但精度一般。分段线性型电压低于起始电压时损耗为零或极小常数超过之后损耗按线性或二次增长。这种方法在暂态仿真里容易被接受因为导数连续性问题不严重。基于V-I特性曲线拟合实测或仿真得到电晕电流与电压的关系曲线比如经典的I C × (U - U₀)^1.5对应空间电荷限制电流再用损耗P U × I折算成电导。这个更接近物理机制。我在模型里采用的是第三种思路的简化版把电晕电流和电压的V-I特性用带阈值和指数项的函数描述同时考虑电压极性和相位。具体来说对交流输电线路电晕损耗正比于电压瞬时值与起始电压之差的某个幂次且损耗在电压峰值附近的贡献最大零交叉附近基本可以忽略。需要特别说明纯交流输电线路里由于电晕损耗要在一个周期内积分所以对瞬时电压的依赖体现在每一个时间步的损耗计算上而如果做的是基波相量域仿真那需要把周期平均损耗和电压有效值之间建立关系。两种域下的模型表达形式不同但底层物理是同一套。我在Simulink里建的是时域模型所以直接处理瞬时值。再补一个容易被忽略的工程细节交流线路电晕损耗在雨天会大幅上升因为雨滴和水膜改变了导线表面状态m系数会往下降等效起始场强随之下降。我的模型里把天气因素折算进粗糙系数m和起始电压U₀的修正系数中留了一个可调参数的入口这样模型既能做晴天基线仿真也能做雨天保守仿真灵活性高很多。2.3 为什么不能忽略电压依赖特性有人可能会问贝杰龙输电线路模型本身已经能处理线路的分布参数特性再加个电晕依赖是不是多此一举这个问题问到了点子上。贝杰龙模型确实把线路的波过程、反射折射都考虑进去了但它假设线路参数——电阻、电感、电容——都是常数不随电压变化。这个假设在电晕起始电压以下的工况里误差很小。但一旦电压冲高电晕出现情况就变了电晕等效于在线路对地之间并联了一个非线性电导而且这个电导随电压增长非常快相当于一条额外的泄放通道把波前能量消耗掉一部分。电晕还让导线等效半径变大导线对地电容跟着变化波的传播速度被改变波形前沿被拉缓。电晕损耗对操作过电压的峰值有抑制作用若忽略它算出来的过电压峰值偏高约10%~20%具体跟线路长度、电压等级、气象条件有关这对绝缘配合影响很大。所以电压依赖模型的价值在于它不是对线路模型的简单修补而是把电晕这个强非线性物理过程以工程上足够精确的方式耦合进暂态仿真主循环里让仿真的电压和时间尺度都能逼近真实物理过程。这正是暂态仿真工具里被反复强调的“模型保真度”问题——线性模型的误差不是叠加在细节上而是叠加在量级和趋势上的。3. Matlab实现方案整体设计3.1 技术选型对比Simulink vs 纯M脚本拿到这个项目第一个要决策的问题就是实现载体。我对比了三条技术路线纯M脚本自定义数值积分优点是完全可控每一步的物理过程都是透明的适合算法验证与教学缺点是写起来繁琐非线性迭代和事件检测比如电压穿越起始阈值都要自己处理后期改模型结构成本高。Simulink Simscape Electrical模块库优点是图形化建模效率高工具箱自带的输电线路模型和元件库丰富适合系统级集成缺点是输电线路库里的模型大多是固定参数版本要改成电压依赖需要自己封装模块调试时看不到内部细节出了问题不好定位。Simulink S函数/MATLAB Function混合方案兼顾两者优点——线路主模型用Simulink框架搭建电压依赖电晕模块用MATLAB Function或Level-2 S函数实现既能图形化集成又能保留自定义非线性逻辑的灵活性。这是我最终选择的路线。配套电气参数计算我放在初始化脚本里用脚本算好线路单位长度参数、贝杰龙模型需要的特征阻抗和波速然后灌入Simulink模型中的参数变量。这一步的好处是参数修改只用改脚本顶部的物理参数完全不用动模型结构工程上复现和调参都方便。3.2 模型架构分区解读整个模型由三大块构成线路主体用贝杰龙模型表达输电线分布参数特性包括单位长度电阻R、电感L、电容C、电导G特征阻抗Zc sqrt((R jωL)/(G jωC))波速v 1/sqrt(LC)不计损耗的近似或按精确公式算衰减系数和传播常数。电晕模块这是核心新增模块。对每一相或每一极线路实时计算导线表面场强判断当前瞬时电压是否超过电晕起始阈值若超过则计算电晕电导和等效附加电容并以受控电流源或可变电导的形式并入线路两端。源与负载部分电源模块用可控电压源模拟不同工况的激励负载端根据研究目标接不同阻抗或开路来模拟空载线路。3.3 提升仿真效率的三个设计节拍第一个设计节拍是解耦。电压依赖电晕模块和贝杰龙波动模型的耦合在数学上是强耦合的——电晕电流会影响线路端电压端电压又反过来决定电晕强度。强行解这个强耦合在每个时间步里做牛顿迭代仿真速度会非常难看。我的处理是用上一个时间步的电压值来估算当前时刻的电晕电导再代入线路方程。因为电晕效应的时间常数相对电磁波传播时间尺度来说不算极快电晕建立和消退有微秒到毫秒量级的过程这种显式解耦在步长短于10微秒时精度完全够用。实测下来在固定步长5微秒、仿真时长20毫秒的工况里解耦误差不超过2%。第二个设计节拍是查表加速。电晕V-I曲线虽然非线性强但它是单变量函数瞬时电压到等效电导完全可以用预计算的查找表加线性插值来替代每个步长的公式计算。我在初始化脚本里生成一条从0到2.5倍额定电压的映射表Simulink里用Lookup Table模块直接查表单步开销几乎可以忽略。第三个设计节拍是参数条件编译。模型里做了一个开关enabled电晕模型的开关变量置0时模型退化为纯贝杰龙模型置1时激活电压依赖模块。这样在调试阶段可以先关掉电晕跑通线路主模型再开电晕做对比验证不搞两套模型省时省力。4. 核心模块实现细节4.1 电压依赖电晕模块的数学方程模块输入是线路对地瞬时电压u(t)输出是电晕电流i_c(t)或等效并联电导g_c(t)。判断依据是起始电压U₀而U₀本身由导线参数和气象条件决定U₀ E₀ × r_eff × ln(D / r_eff)其中D是导线对地等效距离或者几何平均距离针对分裂导线要换算r_eff是导线等效半径E₀由皮克公式得到。注意这个公式来自圆柱导体同轴电容的场强分布表面场强E_surf U / (r_eff × ln(D/r_eff))让E_surf等于E₀反解出来的电压就是起始电压。对于分裂导线等效半径要做修正r_eq r_s × (n × r_s / R_b)^(1/n)其中n是分裂数r_s是子导线半径R_b是分裂间距子导线所在圆的半径。这是工程上常用的几何等效法虽然忽略了子导线之间场强的详细分布但对宏观电晕损耗计算足够准确。电晕电流-电压关系我采用如下形式i_c 0当u ≤ U₀时i_c k × sgn(u) × (|u| - U₀)^m当u U₀时。m取值在1.5到2之间空间电荷限制区更接近1.5强电晕区间更接近2k是比例系数需根据线路实测或者已有文献数据校准。对交流线路分段校准有一个技巧把仿真得到的V-I曲线和参考曲线对比调整k和m最小化误差即可。不需要做物理放电的详细模拟这就是工程模型的取舍。等效电导g_c i_c / u注意这个电导是瞬时值且只在超过阈值的半周期内起作用正负半周对称。如果要做周期平均损耗则对P_c u × i_c在工频周期上积分再除以周期时间。4.2 贝杰龙线路模型与电晕模块的耦合方式贝杰龙模型的核心思想是把分布参数线路等效成两端带电流源的电阻网络线路两端各有一个等效电阻Zc以及一个由对端历史值决定的电流源。用公式表示就是i_k(t) u_k(t) / Zc I_hist_k(t)其中I_hist_k(t)是来自对端的历史波贡献它基于上一时刻的电压电流递归计算。那么电晕模块怎么耦合进去在物理上电晕是线路对地之间的并联导纳所以在贝杰龙等效电路里我把它放到端节点的对地支路上。具体操作是在贝杰龙线路模块的每一端节点并联一个受控电流源i_c它的值由该节点当前电压通过电晕模块实时计算得出。这样一来节点电压方程变为u_k(t) / Zc I_hist_k(t) i_c 外部注入电流这样电晕损耗就直接成为节点电流平衡的一部分不需要改动贝杰龙模型的内部递归结构耦合干净利落。需要提醒一个数值细节i_c的计算要用到当前时刻的u_k但u_k本身就是待求量所以必须做解耦用上一个步长的电压值去查表估算这个步长的电晕电流否则形成代数环。Simulink里如果不处理代数环会导致仿真变慢甚至报错。解决方案就是MATLAB Function模块里把输入配置为前一步电压u_pre输出为当前步电晕电流i_c用unit delay模块做一个步长的延迟即可打破代数环。4.3 Simulink模块搭建步骤整个模型在Simulink里搭建的步骤我整理如下纯图形操作但每一步都有明确的物理含义新建模型仿真步长设为固定步长建议5微秒或者更小。这里必须用固定步长因为变步长求解器在非线性开关类模块下容易漏掉事件点。放置两个贝杰龙输电线路模块对应线路两端中间用分布式参数线段连接或者直接用Simscape的分布式线路模块按线路长度、单位长度参数配置。在每个端节点处添加MATLAB Function模块编写电晕电流计算函数输入为端电压延迟一步后的值输出为电晕电流。将电晕电流通过受控电流源并入节点电流方向设为从节点流向地对应损耗方向。电源侧接交流电压源要模拟操作过电压就加一个阶跃或冲击信号发生器负载侧根据工况接开路或匹配电阻。在模型中添加示波器和To Workspace模块记录端电压和电晕电流波形用于后处理。整个搭建过程在熟练操作下大概一小时可以完成剩下的大部分时间都花在参数校准和验证调试上。4.4 核心M脚本框架在Simulink外部需要一个加载参数的初始化脚本。我给它命名叫init_corona_model.m核心内容如下% 输电线路基本参数 V_rated 500e3; % 额定电压单位V线电压有效值 L_line 300e3; % 线路长度单位m r_s 1.5e-2; % 子导线半径单位m n_split 4; % 分裂导线数 R_bundle 0.45; % 分裂间距单位m D_equiv 12; % 导线对地等效距离单位m m_surf 0.8; % 导线表面粗糙系数 delta_air 1.0; % 相对空气密度 k_corona 2.5e-8; % 电晕电流比例系数需校准 m_exp 1.7; % 电晕电流指数 % 计算等效半径 r_eq r_s * (n_split * r_s / R_bundle)^(1/n_split); % 皮克公式计算起始场强 E0 30.3 * m_surf * delta_air * (1 0.298 / sqrt(r_eq * 100 * delta_air)); % 注意r_eq要换算成cm % 计算电晕起始电压相电压幅值 U0_phase E0 * r_eq * log(D_equiv / r_eq); U0_amp U0_phase * sqrt(2) / sqrt(3); % 折算到相电压幅值粗略估算工程取整 % 建立电晕V-I查找表 u_table linspace(0, 2.5 * V_rated / sqrt(3) * sqrt(2), 2000); i_table zeros(size(u_table)); for idx 1:length(u_table) if u_table(idx) U0_amp i_table(idx) k_corona * (u_table(idx) - U0_amp)^m_exp; end end这段脚本把所有物理参数集中在顶部修改起来非常直观。我建议工程人员拿到代码后第一件事就是调U0_amp和k_corona这两个参数对结果影响最大而且需要结合你实际线路的参数来设置。5. 模型参数校准与结果验证5.1 参数敏感性分析电晕模型里参数不少但不是每个都对结果有同等影响力。我实际做敏感性分析发现按影响排序大概是U0_amp起始电压、k_corona比例系数、m_exp指数、r_eq等效半径、m_surf表面系数。尤其是U0_amp它直接决定电晕在什么电压水平下“起燃”起燃点偏差10%损耗就可能差30%以上。所以校准参数时的第一优先级就是U0_amp。怎么校准U0_amp最稳妥的方式是用电晕起始场强的实测报告数据反推参考同一电压等级、同类型导线的文献值。几个典型数值可以供参考330kV线路四分裂导线电晕起始电压大致在相电压峰值附近500kV线路八分裂导线晴天起始电压约在250kV相电压有效值附近。不同气象条件差异明显雨天起始电压比晴天低15%~25%这是因为水膜降低了表面粗糙系数。5.2 仿真结果对比我拿一个典型算例做对比演示500kV线路、300公里、末端开路合闸空载线路仿真时间20毫秒。分别跑纯贝杰龙模型和电压依赖电晕模型对比末端电压波形。纯贝杰龙模型下合闸暂态末端电压峰值可以达到额定相电压峰值的1.8倍以上波前陡峭反射波多次叠加后峰值高且衰减缓慢。加入电晕模型后末端电压峰值下降约15%波前明显变缓并且后续的振铃幅度衰减更快。原因很直观电晕在电压冲高区间提供了一个额外损耗通道而且电压越高损耗越大相当于一个非线性限幅器把高幅值脉冲削平了。从数值上电晕损耗贡献大约占总注入能量的8%到12%这个量级与工程经验吻合。若天气设置为雨天把m_surf从0.8降到0.6起始电压下降电晕损耗占比会进一步提升到15%以上电压峰值继续下压。这个趋势与实际观测一致。5.3 模型验证要点验证模型对错不能只看波形好看要看物理量是否自洽。我建议做三个检查第一起始电压阈值测试慢慢升高电源电压观察电晕电流出现时刻对应的电压值必须和初始化的U0_amp一致。如果不一致查单位换算——这是最容易出错的地方kV、V、相电压、线电压、峰值、有效值稍不注意就差了根号3或根号2。第二功率守恒检查统计电源注入功率、线路末端吸收功率、电晕模块消耗功率三者的代数和应为零忽略数值误差。如果电晕模块消耗功率为负或异常偏大多半是电晕电流方向设反了或者电压解耦出了问题。第三收敛性检查逐步加密步长从10微秒往1微秒加密观察末端电压峰值是否趋于稳定。如果峰值波动超过设定容许范围比如3%说明步长还不够小需要对模型步长做细化。电晕模块的非线性指数m_exp带来高增益区段太粗的步长会在这个区段引入数值振荡。6. 常见问题与调试实录6.1 代数环与仿真发散第一次搭完模型跑仿真最可能撞见的是Simulink报代数环错误或者仿真直接发散。我在调试过程中遇到过一次原因就是电晕电流计算直接引用了当前步的节点电压形成代数环。Simulink对代数环会尝试迭代求解但非线性强、高增益的环节很容易导致迭代不收敛。解决思路有两个一是在反馈路径上插入Unit Delay模块打破代数环二是把MATLAB Function模块的输入改成上一个采样时刻的电压。方法二更物理也更稳定。我自己最后采用的是方法二并且在延迟环节前做了一次低通滤波进一步抑制高频数值噪声。6.2 参数单位错误导致电晕永不触发或永远触发这个坑太经典了。第一次跑的时候我把U0_amp按相电压有效值计算了但Simulink里加的电压源是线电压有效值峰值一对比模型以为始终超过起始电压电晕电流一直很大波形被削得面目全非。最后逐行核对才发现是单位约定不一致。总结一个经验所有参数在初始化脚本里统一采用V伏特和m米作为基本单位并且用注释明确标注是线电压还是相电压、是峰值还是有效值。一定要多写一行注释不要默认自己记得住。6.3 电晕模块在电压过零附近产生数值抖动电晕电流公式里有|u| - U0_amp这个差值在电压过零附近差值从负到正跳变如果直接用分段函数硬切换导数不连续容易在零点附近产生振荡这个振荡还会沿着贝杰龙模型的历史波传递到对端。我的处理方式是在阈值附近做平滑过渡引入一个过渡宽度参数delta_u比如U0_amp的2%在这个区间内用余弦函数插值让电流从零平滑上升到公式值。这样处理后数值抖动完全消失代价是模型在阈值附近的精度有小幅损失但对暂态研究的精度要求来说完全无关紧要。这也是工程取舍的常规操作——硬切换换成软切换稳定性换局部精度。6.4 查表精度与内存的取舍查找表步长越细插值误差越小但表越长占内存越多。2000个点的表对应2.5倍额定电压范围插值误差已经远小于模型本身的物理不确定性。如果你在嵌入式环境或要批量跑蒙特卡洛仿真可以适当把表压缩到500个点线性插值带来的误差依然可控。但要注意表的区间必须覆盖仿真全过程的电压范围超出表范围会触发Simulink默认的饱和外推这个行为可能掩盖真实的过电压幅值有隐患。6.5 参数速查表我把调试过程中常用到的参数取值范围整理成一张表方便快速试参参数符号典型范围影响起始电压U0_amp0.6~0.95倍相电压峰值决定电晕触发点影响最大电晕比例系数k_corona1e-9 ~ 1e-7决定电晕电流幅值影响损耗量电晕指数m_exp1.5~2.0决定电压升高时损耗增速表面粗糙系数m_surf0.5~0.9修正起始场强雨天取低值等效半径r_eq0.1~0.5m分裂导线影响起始场强和电容参数7. 模型扩展与应用建议7.1 从单相模型到三相模型的扩展很多实际研究需要三相模型因为操作过电压往往是三相联动的。我目前的实现是在单相模型基础上做验证的扩展到三相有两个思路第一个思路是三套独立的单相电晕模型并行各相之间通过线路模块的相间耦合互电感、互电容自然关联。电晕模块本身是各相对地的因此互不影响实现起来只需要把单相模块复制三份然后接入三相贝杰龙线路模块。第二个思路是考虑相间电晕的交互效应——一相电晕产生的空间电荷可能畸变邻近相的表面场强影响邻近相的电晕起始电压。这个效应在紧凑型输电线路里更明显但常规分裂导线间距下可以忽略。所以工程上我建议先用第一个思路除非你很明确要研究相间交互效应否则不要轻易加这个复杂度。7.2 用于过电压峰值预测的实证技巧我做了一组批量仿真把电源电压从0.9倍额定值扫描到2.0倍额定值记录末端电压峰值。结果可以拟合出一条非线性曲线在1.2倍额定电压以下电晕对峰值影响有限峰值基本随电压线性上升超过1.2倍之后曲线斜率明显下降说明电晕损耗的“限幅”作用开始显现。这个规律可以用来快速评估线路的过电压水平不需要每一条线路都跑完整仿真先用归一化曲线做粗略预估再针对关键工况做细致仿真。7.3 从仿真到工程应用的距离作为从业者我最后提醒一句任何仿真模型的结论要用于工程决策比如绝缘配合、避雷器选型都必须经过实验室或现场实测数据的校准和验证。电晕损耗的实测受天气、导线表面老化状态的影响非常大模型参数必须在实际运行数据面前低头而不是反过来让现实去迁就模型参数。这个模型的定位是提供一个稳定可靠的平台把电压依赖电晕效应从“拍脑袋经验”变成“可迭代修正的数字化表达”。有了这个基础后续无论是加装线路避雷器的策略仿真还是高海拔地区线路修正方案验证都能在这个模型上做增量式演进不至于每次从零搭一套。我个人的实际操作体会是电晕模型的最大价值不在于精确复现某一组实验数据而在于把电压依赖这个非线性维度纳入系统级仿真视野。当你习惯了用这种视角看输电线路再去审视那些不做电晕修正的仿真结果你会自动多问一句“电压冲高的时候这个结果还靠得住吗”。就冲这一句花在模型搭建上的时间就值回来了。