1. 项目缘起当“快”与“慢”在电力仿真中相遇在电力系统仿真这个行当里干了十几年我越来越觉得我们很多时候不是在和物理定律较劲而是在和时间尺度与计算资源做一场永无止境的博弈。传统上我们有两类主要的仿真工具一类是电磁暂态EMT仿真它像一台高速摄像机能捕捉到电力系统中微秒级的快速动态比如雷电冲击、开关操作引起的过电压、电力电子器件的快速开关过程。另一类是机电暂态TS仿真它更像一个延时摄影关注的是秒到分钟级的慢速动态比如发电机的功角稳定、频率波动、负荷的缓慢变化。这两者各有各的“主场”EMT擅长细节但算得慢、范围小TS算得快、范围大但对快速过程“视而不见”。然而现代电力系统尤其是随着新能源高比例接入、直流输电和柔性交流输电系统FACTS的广泛应用让“快”与“慢”的界限变得模糊。一个风电场的并网既涉及电力电子变流器毫秒级的开关动态EMT领域又关系到整个区域电网的功角稳定TS领域。这时候如果只用EMT仿真整个大电网计算量会大到令人绝望如果只用TS仿真又会漏掉那些关键的快速动态导致分析结果失真。于是“混合仿真”应运而生它试图让EMT和TS这对“快慢兄弟”联手干活EMT负责仿真需要精细建模的局部快速动态区域TS负责仿真其余的大电网慢速动态部分。我这次要聊的“EMT-TS混合仿真精度评估与三序接口优化研究”正是这个领域里一个既基础又核心的难题。说白了就是当我把EMT和TS两个仿真器“粘”在一起时这个“粘合剂”——也就是接口——到底靠不靠谱它传递信息时会不会“失真”整个混合仿真结果的精度到底如何又该怎么去优化这个接口让它传得更准、更稳这不仅仅是学术问题更是工程实践中决定混合仿真能否真正替代部分全EMT仿真、成为大型电网分析可靠工具的关键。2. 混合仿真的核心接口处的“信息翻译官”要让EMT和TS两个世界顺利对话接口是唯一的桥梁。但问题在于它们俩说的“语言”完全不同。EMT仿真是在三相abc坐标系下进行的它能看到每一相电压、电流的瞬时波形是真实物理世界的直接映射。而TS仿真则通常在正序d-q旋转坐标系或相量域下进行它处理的是经过派克变换或傅里叶变换后的“平均化”或“基波分量”信息本质上是一种简化模型。因此接口的核心任务就是充当“信息翻译官”。它需要从TS侧获取正序电压相量幅值和相位然后“合成”出EMT侧所需的三相瞬时电压边界条件反过来又需要从EMT侧的三相瞬时电流中“提取”出正序电流相量反馈给TS侧。这个“合成”与“提取”的过程就是混合仿真精度问题的根源所在。目前最主流的接口方法是“理想变压器模型ITM”及其各种改进型。它的思想很直观在接口处虚拟一个三相理想变压器TS侧提供正序电压作为这个虚拟变压器的原边电压从而在EMT侧感应出三相电压EMT侧的三相电流则“流回”这个虚拟变压器在TS侧等值为正序电流注入。听着很完美对吧但魔鬼藏在细节里。2.1 接口处的“失真”从何而来在实际操作中我发现接口的失真主要来自以下几个层面它们环环相扣数据交换的“时间差”EMT和TS的仿真步长天差地别。EMT步长通常在几十微秒而TS步长在几毫秒到几十毫秒。接口数据交换不可能每个EMT步长都进行一次那样通信开销太大。通常是在每个TS步长处进行交换。这就意味着在一个TS步长内EMT侧是在一个“恒定”的边界电压下运行的这个电压是上个TS时刻TS侧提供的而TS侧接收的电流反馈也是EMT侧在一个完整TS步长内动态过程的“总体效果”。这个“保持”和“平均”的过程必然会引入误差尤其是在EMT侧动态剧烈变化时。序分量转换的“理想化假设”三序接口正序、负序、零序的理论基础是系统对称或基本对称。它默认EMT侧网络在接口处是三相平衡的或者不平衡分量可以完美分解为正、负、零序。但在真实的电力电子化电网中接口处很可能因为非线性元件、不对称故障或控制策略而产生丰富的谐波和不平衡分量。此时强行用基波正序分量去“代表”EMT侧的复杂波形就像用一个人的平均体温去描述他发烧时每分钟的体温变化一样会丢失大量关键信息。负序和零序分量的引入就是为了部分弥补这种不平衡但其建模和传递的准确性又是一大挑战。网络等值的“简化伤”为了解耦TS侧看到的EMT部分通常被等值为一个诺顿或戴维南等效电路一个电流源并联阻抗或电压源串联阻抗。这个等效阻抗的参数通常是次暂态或暂态电抗如何选取它是线性的还是非线性的在实际仿真中EMT侧系统的等效阻抗会随着运行点变化特别是含有电力电子设备时。用一个固定阻抗去代表在系统大扰动时必然会产生误差。数值振荡的“幽灵”这是最棘手的问题之一。由于接口处电压和电流的相互依赖关系以及数据交换的延时很容易引发数值不稳定性。表现为仿真波形中出现高频振荡这种振荡没有物理意义纯粹是数值计算耦合产生的“幽灵”。它可能淹没真实的物理动态导致仿真失败。3. 精度评估不能只看“长得像”要问“为什么不像”当我们拿到一份混合仿真的结果比如一条母线电压的波形我们如何判断它的精度很多初入行的朋友会习惯性地把混合仿真的结果和全EMT仿真的结果叠在一起看如果两条曲线重合度很高就认为精度好。这没错但不够。作为深度实践者我认为精度评估必须是一个多层次、多维度的“体检”过程。3.1 建立可靠的“金标准”首先必须有一个无可争议的参照系即“金标准”。对于混合仿真最理想的金标准是对同一个完整系统进行全EMT仿真。但全EMT仿真计算代价极高有时甚至不可行。因此在实践中我们可以采用一种“阶梯式”验证策略层级一简单测试系统。构建一个包含典型元件的简单系统如同步机、变压器、线路、负荷分别进行全EMT仿真和混合仿真。这个系统要能涵盖我们关心的主要动态过程。这里的全EMT仿真是可行的结果作为金标准。层级二与专业TS仿真结果对比。在TS侧将混合仿真中TS部分的动态响应如发电机功角、频率与用专业TS仿真软件如PSASP、BPA对同一TS网络进行独立仿真的结果进行对比。这可以验证接口信息传递没有破坏TS侧原有的动态特性。层级三局部精细化验证。对于大规模系统可以对EMT侧关注的局部区域将其从大系统中“切割”出来赋予一个合理的边界条件可从混合仿真或TS仿真中获取进行独立的、更精细的全EMT仿真然后将这个局部结果与混合仿真中EMT侧的结果进行对比。3.2 量化误差不止于RMSE对比波形时除了肉眼观察必须引入量化指标。均方根误差RMSE和平均绝对百分比误差MAPE是常用的但它们主要反映整体偏差。对于电力系统仿真我们更应关注关键特征点的精度暂态峰值误差过电压、冲击电流的最大值偏差是多少这直接关系到设备绝缘设计和保护定值。上升/下降时间误差波形的陡度是否一致这影响了对于波前时间敏感的器件如避雷器的分析。稳定值误差暂态过程结束后稳态值的偏差。这关系到系统最终的运行状态。主要振荡频率与阻尼对于功率振荡或谐振其频率和衰减系数是否吻合这关乎系统稳定性判断。我通常会制作一个误差分析表格将上述指标在全EMT结果和混合仿真结果之间一一对比。评估指标全EMT仿真值混合仿真值绝对误差相对误差可接受阈值是否达标故障点电压暂态峰值 (kV)1.251.18-0.07-5.6% ±10%是故障切除后振荡频率 (Hz)1.321.28-0.04-3.0% ±5%是发电机功角最大摆幅 (度)45.648.12.55.5% ±8%是联络线功率稳态值 (MW)520508-12-2.3% ±3%是3.3 压力测试在极端场景下暴露问题精度评估绝不能只在“风平浪静”的稳态下进行。必须设计一系列压力测试场景把接口推到极限大扰动测试模拟三相短路、大型发电机跳闸等严重故障。观察接口在功率和电压剧烈突变时能否保持稳定和准确。不对称故障测试模拟单相接地、两相短路等。这是检验三序接口处理不平衡能力的关键。看负序和零序分量的传递是否准确会不会引起错误放大。非线性元件测试在EMT侧引入饱和变压器、电弧模型等强非线性元件。观察接口算法能否处理由此产生的谐波和非周期分量。弱电网连接测试让EMT侧系统通过一个高阻抗长线路或小容量变压器连接到TS主网。这种场景下接口点的电压支撑弱对接口算法的鲁棒性要求极高。提示在进行压力测试时一定要同步监测接口交换点的数据。我习惯将每个交换时刻TS侧发出的电压相量幅值、相位和EMT侧反馈的电流相量都记录下来绘制成随时间变化的曲线。观察这些曲线是否平滑有无跳变或振荡这往往是发现接口算法潜在问题的直接窗口。4. 三序接口的优化实战从理论到参数调校当我们通过精度评估发现了接口存在的问题接下来就是优化。优化不是推倒重来而是一个从模型、算法到参数的精细化调校过程。4.1 接口模型的改进超越理想变压器基础的ITM模型假设接口变压器是理想的、无损耗的且TS侧提供的正序电压能无畸变地产生完美三相正弦波。这对于强电网、对称情况尚可但对于现代电网远远不够。实践中我主要从以下几个方向进行模型改进引入耦合阻抗在虚拟变压器原副边之间加入一个串联阻抗通常用TS侧系统的短路阻抗或更精细的等值阻抗。这个阻抗起到了“缓冲”作用可以一定程度上吸收由于两侧步长不一致和网络不等值带来的功率不平衡减少数值振荡。其参数需要根据TS侧系统的戴维南等效阻抗进行仔细整定。改进序分量计算传统的基于单周期基波提取的正负零序计算方法在波形畸变时误差大、延时高。可以采用更先进的信号处理技术如瞬时对称分量法通过一个变换矩阵直接从abc瞬时值计算瞬时正、负、零序分量响应更快但受高频噪声影响大。锁相环PLL增强型提取利用高性能PLL如基于二阶广义积分器的SOGI-PLL在非理想电网条件下快速、准确地提取基波正序分量。对于负序和零序也可以设计相应的提取器。移动窗口DFT/FFT采用一个滑动数据窗口进行傅里叶分析可以在每个仿真步长都更新序分量平衡了精度和速度。考虑电磁耦合对于同塔多回线路等存在电磁耦合的情况简单的三序解耦可能不充分。此时需要考虑更详细的相域模型与序域模型之间的转换关系或者在接口处保留部分相域模型特性。4.2 数据交换策略的优化寻找节奏感数据交换的节奏何时交换、交换什么直接影响精度和稳定性。插值法TS侧在每个大步长给出一个电压相量。如果EMT侧直接在整个TS步长内使用这个恒定值误差大。改进方法是让TS侧提供未来一个步长内的电压预测值例如线性预测或者让EMT侧在内部采用插值算法如线性插值、三次样条插值根据TS侧当前和上一个时刻的电压推测出本TS步长内每个EMT小步长的边界电压。这能显著平滑边界条件提高精度。迭代法或预报-校正法在每个TS步长内进行多次数据交换迭代直到两侧接口处的功率或电压误差小于某个阈值。这类似于一个微型的协同仿真迭代过程精度最高但通信开销和计算量也最大通常用于对精度要求极高的关键研究而非大规模实时仿真。多速率接口这不是简单地固定交换频率而是根据系统动态的剧烈程度自适应调整。当系统处于稳态或慢动态时降低交换频率当检测到故障或快速扰动时立即提高交换频率。这需要设计灵敏的“事件检测”逻辑。4.3 阻尼与滤波给接口装上“减震器”数值振荡是混合仿真的顽疾。直接硬抗往往适得其反。我的经验是适当引入“阻尼”或“滤波”像给接口装上减震器。接口阻尼电阻在EMT侧的接口节点与地之间或在虚拟变压器支路上并联一个高值电阻。这个电阻不显著影响工频特性但可以为高频的数值振荡能量提供一个泄放通道有效抑制振荡。电阻值的选择是关键太小会影响系统动态太大则阻尼效果不足通常需要经过多次测试确定。一阶滞后滤波对从TS侧接收的电压指令或向TS侧发送的电流反馈施加一个一阶低通滤波。公式类似于y(k) α * x(k) (1-α) * y(k-1)。其中α是滤波系数介于0和1之间。这相当于给交换信号增加了一个惯性环节平滑了突变抑制了高频分量但会引入相位滞后。需要权衡稳定性和动态响应速度。前馈补偿为了抵消滤波带来的相位滞后可以采用前馈补偿。例如在已知系统主要动态频率的情况下设计一个超前补偿器对信号进行相位超前校正使整体相频特性在关注频段内保持平坦。4.4 参数整定一场精细的“微调手术”所有上述改进模型和策略都离不开参数整定。这是一个经验与科学结合的过程。我通常遵循以下步骤确定调参目标是优先保证精度还是优先保证稳定性对于不同的研究目的过电压研究、稳定性研究侧重点不同。单参数扫描在简单测试系统上保持其他参数不变逐个调整关键参数如耦合阻抗值、阻尼电阻值、滤波系数α观察其对精度指标如峰值误差和稳定性仿真能否顺利完成波形有无振荡的影响。记录下参数变化与效果之间的趋势。多参数协同优化关键参数之间往往相互影响。例如增大阻尼电阻可能允许使用更大的滤波系数。这时可以采用实验设计DOE的方法或者基于优化算法如粒子群、遗传算法以综合误差最小或稳定裕度最大为目标进行多参数寻优。虽然计算量大但对于一个要反复使用的标准接口模型这个投入是值得的。在压力测试下验证将初步整定好的参数放到第3章设计的各种压力测试场景中去跑。观察在极端条件下参数是否依然鲁棒。往往需要根据压力测试结果进行微调。5. 一个典型优化案例风电并网系统的混合仿真接口调校让我用一个简化但典型的案例把上面的理论串起来。假设我们研究一个通过长线路接入主网的风电场EMT侧对主网TS侧暂态稳定性的影响。5.1 初始设置与问题浮现最初我们使用标准的三序ITM接口无阻尼无滤波TS步长20msEMT步长50μs。在模拟主网发生三相短路故障时混合仿真出现了问题现象故障期间及切除后风电场出口电压和电流波形出现明显的高频毛刺数值振荡且风电场输出的有功功率振荡衰减很慢与全EMT仿真结果相比阻尼偏小。精度评估峰值电压误差在8%以内但功率振荡频率误差达12%阻尼时间常数误差超过30%。这会导致对系统稳定性的误判。5.2 分步优化过程第一步抑制数值振荡保稳定。首先在EMT侧并网点的每相对地之间并联一个阻尼电阻。从1kΩ开始尝试发现振荡有所减弱但未消除。逐步减小到500Ω振荡基本消失但通过录波发现在故障瞬间引入了一个很小的额外有功冲击约0.5%。权衡后选择600Ω作为折中。第二步改进电压边界平滑度提精度。观察到TS侧提供的电压指令在故障瞬间是阶跃变化的。我们在EMT侧接口程序中加入了线性插值模块。利用当前TS时刻和上一TS时刻的电压相量为中间每一个EMT步长插值计算边界电压。这一步后电压电流波形的平滑度显著提升峰值误差减小到5%。第三步优化序分量提取针对不平衡。为了更精确地模拟单相接地故障我们将接口从仅正序扩展为完整三序。但发现传统的基于整周期DFT的负序提取在故障初期响应慢。我们换用了基于双二阶广义积分器DSOGI的瞬时对称分量法虽然计算量稍增但负序分量的提取速度和精度在动态过程中明显改善单相故障下的仿真精度大幅提高。第四步精细调参协同优化。现在我们有三个可调参数阻尼电阻R_d、插值算法的启用/禁用、序分量提取器的带宽。我们设计了一个简单的网格搜索在单相故障和三相故障两种场景下组合不同的参数计算综合误差指标RMSE和稳定时间误差的加权和。最终找到一组较优参数R_d600Ω启用插值DSOGI带宽设为30Hz在动态响应和抗噪间平衡。5.3 优化效果对比优化后重新进行压力测试数值振荡完全消除。精度在各种故障类型下关键电气量的峰值误差稳定在5%以内振荡频率和阻尼误差控制在8%以内达到工程实用要求。稳定性在更弱的电网条件下短路比更小仿真也能稳定进行。这个案例说明接口优化是一个系统工程需要先解决稳定性振荡再解决精度误差并且要根据具体的仿真对象如风电场的电力电子特性选择合适的模型和算法最后通过系统性的参数整定达到最佳效果。没有一劳永逸的“万能接口”只有针对特定场景不断打磨的“专用接口”。每一次成功的混合仿真应用背后都离不开对接口精度与优化技术的深刻理解和细致实践。