变压器纵联差动保护的仿真,是电力系统方向里一个特别容易“看起来简单、跑起来打脸”的活儿。我见过太多人拿着教材上的三段式比率制动特性曲线,在Simulink里搭完模型,波形一跑,区外故障时保护动作了,区内故障时反而纹丝不动。然后就开始怀疑人生:公式没错啊,整定值也是按规程算的,怎么就误动拒动了?问题往往不在保护原理本身,而在于仿真模型里那些教材不会告诉你、但实际搭建时处处埋雷的细节——CT饱和特性的建模方式、两侧电流相位的折算、制动量的取法、故障初相角的影响、离散化步长的选择,每一个都能让你的仿真结果偏离预期。这篇内容就是把这些年我在变压器差动保护仿真里踩过的坑、排查过的异常、以及最后总结出来的一套可复现的建模思路,完整地摊开讲一遍。不管你是刚接触继电保护仿真的学生,还是需要快速验证保护逻辑的工程人员,下面这些内容应该都能帮你少走不少弯路。
1. 先搞清楚:你的模型到底在“仿”什么
很多人搭差动保护仿真,一上来就打开Simulink拖模块,变压器、CT、继电器、断路器一顿连,跑出来不对再回头改。这种“先搭后想”的路子,在差动保护这个场景里几乎必然翻车。因为差动保护的仿真本质上是在复现一个“比较”过程——比较两侧电流的幅值和相位关系,而这个比较的准确性,取决于你对每一个环节的物理建模是否到位。所以在动手之前,得先把“仿什么”这件事想清楚。
1.1 差动保护的判据本质:一个被很多人误解的“和”
纵联差动保护的核心判据,说穿了就是比较变压器各侧电流的相量和。理想情况下,正常运行和区外故障时,流入变压器的电流等于流出的电流,差动电流为零;区内故障时,故障点提供了额外的电流通路,差动电流等于故障电流。这个逻辑本身没有任何问题,但问题出在“相量和”这三个字上。
变压器各侧的电流,由于绕组接线方式不同(比如Y/Δ接法),两侧电流之间存在30度的相位差,幅值也因为变比关系不相等。你在仿真里如果直接把两侧CT二次电流拿来做减法,得到的差动电流在正常运行时就有一个很大的稳态值,保护自然误动。所以模型里必须有一个相位和幅值的折算环节,把两侧电流归算到同一侧、同一相位基准下。这个折算在真实保护装置里是由CT接线组别和装置内部算法共同完成的,在仿真里则需要你显式地建模出来。
我见过一个典型的错误做法:有人在CT二次侧直接加了一个移相模块,把相位硬掰30度,幅值用理想变比折算。乍一看没问题,但一跑区外故障,差动电流还是不为零。原因在于,他用的移相模块是理想移相,而实际变压器在暂态过程中,两侧电流的相位关系并不是恒定30度,铁芯饱和、励磁涌流都会让这个关系发生偏移。所以相位折算环节的建模,必须和变压器的电磁暂态模型耦合起来考虑,不能孤立地做一个理想移相。
1.2 仿真精度和模型复杂度的取舍
另一个在动手前要想清楚的问题是:你这个仿真到底要验证什么?如果只是验证保护判据的逻辑正确性,那变压器模型可以用简化的理想变压器加集中参数,CT用线性变比模型,重点放在继电器逻辑上。但如果你要复现误动、拒动这类异常工况,那就必须把CT饱和、铁芯非线性、励磁涌流这些因素建模进去,否则你永远复现不出真实装置在现场遇到的那些边界情况。
这里有个经验:差动保护仿真里,最容易被简化掉、但又最影响结果的两个环节,一个是CT的饱和特性,一个是变压器的励磁支路。前者决定了区外故障穿越性电流下的差动不平衡电流大小,后者决定了空载合闸时的涌流波形。这两个环节如果都用理想模型代替,你的仿真结果会“过于干净”,干净到和现场完全对不上。所以我的建议是,哪怕其他环节简化,这两个地方也要用带非线性特性的模型。
2. 电流折算环节:差动保护仿真里最容易埋雷的地方
如果说差动保护仿真有一个“万恶之源”,那一定是电流折算环节。我排查过的误动案例里,超过一半最后都定位到这个环节的建模错误上。这个环节涉及相位折算、幅值折算、以及两者在暂态过程中的耦合,任何一个处理不当,都会在差动电流里引入虚假的不平衡量。
2.1 相位折算:不是简单移相30度就完事
变压器Y/Δ接线带来的30度相位差,是差动保护必须处理的问题。在仿真里处理这个相位差,常见的有三种做法,每种都有各自的适用场景和坑。
第一种是CT二次侧接线组别补偿,也就是模拟真实装置里通过CT的Δ接法来抵消相位差。这种做法在仿真里实现起来比较麻烦,因为你要在CT二次侧搭一个Δ接法,涉及到电流的相减和再分配。好处是它最接近真实装置的物理过程,暂态响应也最真实。缺点是建模复杂,而且Δ接法CT在内部故障时会有灵敏度问题。
第二种是装置内部软件移相,也就是在保护算法里用数学方法把两侧电流折算到同一相位。这是目前数字保护装置的主流做法,仿真里也最常用。具体实现上,可以用全周傅里叶算法提取基波相量后做旋转,也可以在时域里用延时采样来移相。这里有个坑:时域移相如果用的是固定延时,在频率偏移或者暂态过程中会引入误差,所以最好用相量域的旋转,配合一个能跟踪频率的算法。
第三种是直接在三相电流上做Clarke变换或者对称分量分解,把三相电流分解成正序、负序、零序分量,然后在序分量域里做折算。这种做法在处理不对称故障时比较优雅,但对建模者的理论功底要求高,而且序分量域的判据和相分量域的判据之间的对应关系容易搞混。
我个人的建议是,如果你的仿真重点是保护判据逻辑,用第二种(相量域软件移相)最省事也最可控;如果重点是暂态过程复现,用第一种(CT接线补偿)更真实。但无论用哪种,都要注意一个细节:相位折算的基准侧要固定,不能两侧都折算,否则会引入额外的相位误差。
2.2 幅值折算:变比、CT变比、平衡系数三者别搞混
幅值折算这块,混乱程度比相位折算有过之而无不及。变压器有变比,CT有变比,保护装置里还有平衡系数,这三个东西在仿真里如果混在一起处理,很容易算错。
正确的逻辑链条是这样的:变压器两侧的一次电流,先经过各自CT变比变换成二次电流,然后保护装置用平衡系数把两侧二次电流归算到同一个基准。平衡系数的计算,要同时考虑变压器变比和两侧CT变比的差异。具体来说,如果以高压侧为基准,低压侧的平衡系数应该等于(变压器变比 × 低压侧CT变比)/(高压侧CT变比),再取倒数或者按装置定义调整。
仿真里常见的错误是:直接把变压器变比当成平衡系数用,忽略了CT变比的影响。或者反过来,把CT变比折算做了一次,又在平衡系数里重复折算了一次,导致差动电流被放大或缩小了一个系数。这种错误在正常运行工况下可能看不出来(因为差动电流本来就接近零),但一跑区外故障,差动电流的数值就会明显偏离预期。
提示:在仿真里做幅值折算时,建议把每一步的折算系数单独列出来,用一个表格记录每个环节的输入输出关系,跑仿真前先手算一遍额定工况下的差动电流,确认接近零再往下走。这一步花十分钟,能省后面几个小时的排查。
2.3 暂态过程中的折算误差:一个容易被忽略的细节
相位和幅值折算在稳态下做对了,不代表暂态下也对。变压器在暂态过程中,两侧电流的相位关系会随着铁芯饱和程度变化而偏移,这时候如果折算环节用的是固定参数,就会在差动电流里引入暂态不平衡量。
这个问题的根源在于,真实的差动保护装置在处理相位折算时,用的算法本身就有一定的自适应能力(比如跟踪频率变化、根据电流大小调整算法参数),而仿真里如果用一个固定的移相模块,就丢失了这种自适应能力。解决办法有两个:一是把折算环节做成参数可调的,根据仿真工况手动调整;二是在折算环节后面加一个滤波或者补偿环节,把暂态不平衡量压下去。前者适合做机理分析,后者适合做保护逻辑验证。
3. CT饱和建模:区外故障误动的头号嫌疑人
区外故障时差动保护误动,十有八九和CT饱和有关。这个结论在仿真里同样成立。但CT饱和的建模,恰恰是很多人偷懒的地方——直接用理想变比模型代替,或者用一个简单的限幅环节模拟饱和,结果就是仿真里永远复现不出区外故障误动的现象,然后你就以为自己的保护逻辑没问题,到了现场才发现不是那么回事。
3.1 CT饱和的物理本质和建模要点
CT饱和的本质是铁芯磁通密度超过饱和点后,励磁阻抗急剧下降,导致一次电流中有很大一部分变成了励磁电流,二次电流不再按变比线性跟随一次电流。在仿真里要复现这个过程,核心是建立一个能反映铁芯非线性磁化特性的模型。
常见的建模方式有三种。第一种是用Jiles-Atherton磁滞模型,这个模型能比较准确地反映铁芯的磁滞回线,但参数辨识比较麻烦,需要实际的铁芯材料数据。第二种是用分段线性化的励磁电感模型,把励磁电感随磁通变化的曲线用几段折线近似,实现起来简单,精度也够用,是仿真里最常用的做法。第三种是用查表法,直接把实测的励磁特性曲线做成查找表,精度最高但需要实测数据。
对于大多数差动保护仿真场景,我推荐用分段线性化模型。具体做法是:把CT的等效电路搭出来,一次侧是电流源(代表一次电流),二次侧是负载,中间是励磁支路。励磁支路的电感用一个受磁通控制的可变电感,磁通通过对励磁电压积分得到,当磁通超过饱和阈值时,电感值切换到饱和区的小电感值。这个模型在Simulink里用Simscape或者基本的积分器和查表模块都能搭出来。
这里有个关键参数:饱和阈值和饱和后的电感值。这两个参数直接决定了CT饱和的深度和二次电流的畸变程度。如果饱和阈值设得太高,CT永远不饱和,区外故障误动就复现不出来;设得太低,CT过早饱和,差动电流畸变过大,又会掩盖其他问题。我的经验是,按照实际CT的准确限值系数(ALF)来设定饱和阈值,饱和后的电感值取正常励磁电感的几十分之一到几百分之一,具体看CT的剩磁情况。
3.2 剩磁:一个让仿真结果“随机”的变量
CT铁芯的剩磁,是差动保护仿真里一个特别讨厌的变量。它的大小和方向取决于CT上一次断电时的状态,在仿真里如果不显式设定,默认是零。但实际运行中,CT铁芯往往带有不同程度的剩磁,剩磁会显著改变CT的饱和时机和饱和深度。
在仿真里处理剩磁,有两种方式。一种是把它作为一个初始条件,在仿真开始前给励磁支路的积分器设一个初始磁通值。另一种是在仿真过程中通过一个预充电环节来模拟剩磁的建立过程。前者简单直接,适合做参数扫描分析;后者更真实,但实现起来复杂。
我做过一组对比仿真:同样的区外故障工况,剩磁为零时差动电流的不平衡量很小,保护不误动;剩磁为正饱和方向时,CT在故障暂态过程中提前饱和,差动电流出现明显的畸变分量,保护误动。这个对比说明,如果你的仿真里不设定剩磁,你可能会得出“保护不会误动”的错误结论。所以做区外故障误动分析时,剩磁这个变量一定要纳入考虑,至少要做几个不同剩磁水平的对比。
3.3 二次负载对CT饱和的影响:别把负载设成纯电阻
CT二次负载的建模,很多人图省事直接用一个电阻代替。这在稳态分析里问题不大,但在暂态分析里会引入误差。实际的CT二次负载包括继电器线圈、电缆阻抗等,是感性的,而且在不同频率下阻抗特性不同。
仿真里如果只用纯电阻负载,CT饱和后的二次电流波形会和实际有差异,尤其是在高次谐波分量上。这个差异会进一步影响差动保护算法里傅里叶滤波的结果,导致差动电流的计算出现偏差。所以建议二次负载用电阻串联电感的模型,电感值根据实际电缆和继电器参数估算。如果实在没有参数,用一个小的电感值(比如几毫亨)也比纯电阻强。
4. 励磁涌流:空载合闸时躲不过去的坎
变压器空载合闸产生的励磁涌流,是差动保护必须躲过的工况。仿真里要复现涌流,并且验证保护能不能可靠躲过,是差动保护仿真的另一个重头戏。但涌流的仿真,坑一点也不比CT饱和少。
4.1 涌流波形对变压器模型参数的敏感性
励磁涌流的波形特征——幅值大小、衰减速度、二次谐波含量——对变压器模型的参数非常敏感。尤其是铁芯的饱和特性和剩磁,这两个参数稍微变一点,涌流波形就明显不同。
仿真里搭变压器模型,如果用的是线性变压器模型,那根本产生不了涌流,因为线性模型没有饱和特性。必须用带饱和特性的变压器模型,把励磁支路的非线性伏安特性建进去。Simulink的Simscape Electrical里有现成的饱和变压器模型,但默认参数不一定合适,需要根据实际变压器的铭牌数据和空载试验数据来调整。
这里有个实操经验:涌流的二次谐波含量,是差动保护二次谐波制动判据的关键输入。如果变压器模型的饱和特性设得太“硬”(饱和后电感下降太快),涌流的二次谐波含量会偏高;设得太“软”,二次谐波含量偏低。这两种情况都会导致你对二次谐波制动系数的整定产生误判。所以建议在仿真里做一组饱和特性参数的扫描,看看二次谐波含量在什么范围内变化,再结合保护定值来判断。
4.2 合闸初相角:决定涌流大小的关键变量
涌流的大小和合闸初相角密切相关。当合闸时刻电压过零时,磁通积分从零开始,最容易产生最大的涌流;当合闸时刻电压峰值时,涌流最小。仿真里如果不做初相角的扫描,只在一个固定相角下跑,很可能错过最严重的涌流工况。
我的做法是,在合闸控制环节里加一个可调的延时,让合闸时刻相对于电压波形的相位可以扫描。一般扫0到180度,步长15度,看涌流幅值和二次谐波含量的变化范围。这个扫描不需要跑很多次,十几次就能看出趋势。然后取最严重工况下的涌流波形,去验证保护的二次谐波制动能不能可靠闭锁。
4.3 涌流衰减和CT饱和的耦合:一个复合工况
涌流工况下,CT也可能饱和,尤其是涌流幅值大、衰减慢的时候。涌流导致的CT饱和,和故障电流导致的CT饱和,特性不一样:涌流是单极性的,CT饱和后二次电流的畸变更偏向于单方向,这会影响差动保护算法里的一些对称性假设。
仿真里要复现这个复合工况,需要把涌流模型和CT饱和模型耦合起来跑。这时候计算量会大一些,步长要设得小一点,否则积分误差会累积。我一般把仿真步长设在微秒级,用变步长求解器,在饱和发生的时刻自动缩小步长。
5. 保护算法实现:离散化带来的那些“意外”
前面的建模环节都做对了,最后一步是把保护判据用代码或者模块实现出来。这一步看似简单,但离散化实现里藏着不少能让结果跑偏的细节。
5.1 傅里叶滤波器的窗长和采样率匹配
差动保护算法里,提取基波相量最常用的是全周傅里叶算法。这个算法需要一个工频周期的数据窗,窗长和采样率必须匹配。仿真里如果采样率设得和窗长不匹配,比如一个周期采了20个点但窗长设成了24个点,滤波结果就会出错。
正确的做法是:先确定仿真采样率(比如每个周期64点),然后傅里叶窗长严格等于一个周期的采样点数(64点)。如果仿真步长和采样率不一致,还要加一个采样保持环节,把仿真步长上的连续值采样成离散序列再送进傅里叶算法。这个细节如果处理不好,差动电流的计算会有一个随机的偏差,时大时小,很难排查。
5.2 制动量的取法:标积制动和比率制动的区别
差动保护的制动量,常见的有比率制动和标积制动两种取法。比率制动取两侧电流幅值的最大值或者平均值,标积制动取两侧电流相量的点积。这两种取法在仿真里的实现不一样,对故障类型的响应也不一样。
比率制动实现简单,在区外故障时制动量随穿越性电流增大而增大,能有效防止误动。但在区内故障时,如果两侧电流都有贡献(比如变压器带负荷时发生区内故障),制动量也会比较大,可能影响灵敏度。标积制动在区内故障时制动量小,灵敏度高,但在区外故障严重饱和时,制动量的计算容易受畸变分量影响。
仿真里选哪种,取决于你要验证的保护装置用的是哪种判据。如果装置说明书里写的是比率制动,你仿真里用标积制动去验证,那结果对不上是必然的。所以这一步一定要和实际装置的判据对齐,不能想当然。
5.3 二次谐波制动:闭锁逻辑的时序问题
二次谐波制动是躲涌流的主要手段,仿真里实现这个逻辑,除了要算对二次谐波含量,还要注意闭锁逻辑的时序。涌流工况下,二次谐波含量是随时间变化的,合闸初期含量高,随着涌流衰减逐渐降低。如果闭锁逻辑的延时设得太短,可能在涌流后期二次谐波含量降低时误开放;设得太长,又可能在涌流已经衰减到不足以引起误动时还闭锁着,影响保护的快速性。
仿真里验证这个逻辑,要把涌流全过程跑出来,看二次谐波含量随时间的变化曲线,再结合闭锁逻辑的延时参数,判断闭锁和开放的时机是否合理。这个分析用波形图看最直观,建议把二次谐波含量、差动电流、制动量、闭锁信号画在一张图上,时序关系一目了然。
6. 排查链路:当仿真结果不对时,按这个顺序查
前面讲了建模和算法实现的要点,但实际做仿真,结果不对是常态。关键是要有一套系统的排查链路,而不是东改一下西改一下。下面这个排查顺序,是我这些年总结下来效率最高的。
6.1 第一步:确认稳态额定工况下差动电流接近零
这是最基础的一步,也是最能暴露折算环节错误的一步。在变压器带额定负荷、无故障的工况下跑仿真,看差动电流的有效值。如果差动电流明显不为零(比如超过额定电流的百分之几),那一定是相位折算或者幅值折算出了问题。
排查方法:把两侧折算后的电流波形画出来,看幅值是否相等、相位是否一致。如果幅值不等,查幅值折算系数;如果相位不一致,查相位折算环节。这一步不要跳过,因为稳态下的折算错误会在所有后续工况里持续存在,掩盖其他问题。
6.2 第二步:区外故障下看差动不平衡电流是否合理
稳态确认无误后,跑一个区外故障工况(比如变压器出口处短路),看差动电流的不平衡量。理想情况下,区外故障时差动电流应该很小,只有CT传变误差和饱和引起的不平衡量。如果差动电流很大,先查CT饱和模型是否设置合理,再查两侧电流的暂态折算是否匹配。
这一步的关键是区分“合理的不平衡量”和“错误的不平衡量”。合理的不平衡量通常在额定电流的百分之几到十几之间,随故障电流增大而增大,波形上有明显的饱和畸变特征。错误的不平衡量往往数值异常大,或者波形特征不对(比如没有饱和畸变但数值很大,那可能是折算系数错了)。
6.3 第三步:区内故障下看保护动作是否正确
区外工况确认后,跑区内故障,看保护是否可靠动作。如果拒动,先查制动量是否过大(可能是制动量取法不对,或者两侧电流都贡献导致制动量偏大),再查故障电流是否足够大(可能是故障点设置或者系统阻抗参数问题)。
如果误动或者动作时间不对,查傅里叶滤波器的窗长和采样率匹配,查二次谐波闭锁逻辑是否误闭锁。区内故障时二次谐波含量通常很低,不应该触发闭锁,如果触发了,说明谐波计算有问题。
6.4 第四步:涌流工况下看闭锁是否可靠
最后跑空载合闸工况,看二次谐波制动能不能可靠闭锁。如果闭锁失败,查涌流模型的饱和特性参数,查二次谐波含量的计算结果,查闭锁逻辑的延时设置。这一步要结合合闸初相角的扫描一起做,确保在最严重涌流工况下也能可靠闭锁。
这个排查链路的核心逻辑是:从稳态到暂态,从区外到区内,从故障到涌流,每一步都建立在前一步正确的基础上。如果跳过前面的步骤直接查后面的问题,很容易被前面的错误误导,浪费大量时间。
7. 几个让我印象深刻的踩坑实录
理论讲完了,分享几个实际踩过的坑,都是那种“当时觉得没问题、后来发现错得离谱”的类型,希望能帮你提前避开。
7.1 那个让差动电流“凭空多出来”的接地方式
有一次搭一个Y/Δ接线的变压器差动保护模型,稳态工况下差动电流怎么都调不到零,总有一个不大不小的稳态分量。查了相位折算、幅值折算、平衡系数,都没问题。最后发现是变压器中性点接地方式设错了——模型里Y侧中性点直接接地,但实际系统的接地方式是通过一个接地变压器或者消弧线圈接地,中性点电流的路径不一样,导致零序分量在差动回路里产生了额外的差流。
这个坑的教训是:变压器中性点接地方式,在差动保护仿真里不是一个小细节,它直接影响零序电流的流通路径,进而影响差动电流。建模时一定要和实际系统的接地方式对齐,不能想当然地直接接地。
7.2 傅里叶窗长差一个点,结果差出十万八千里
还有一次,仿真里差动电流的计算结果总是有一个小幅度的周期性波动,幅值不大但一直存在,导致保护在区外故障时偶尔误动。查了很久,最后发现是傅里叶滤波器的窗长设成了63个采样点,而采样率是每周期64点。就差这一个点,滤波器的正交性被破坏,基波提取出现了泄漏,差动电流里混进了一个与采样率相关的纹波。
这个坑特别隐蔽,因为差一个点不会导致结果完全错误,只会引入一个小偏差,很容易被当成“正常的仿真误差”忽略掉。但差动保护的定值往往很敏感,这个小偏差在边界工况下就可能导致误判。所以傅里叶窗长和采样率的匹配,一定要严格检查。
7.3 求解器步长设太大,饱和波形全是毛刺
CT饱和的仿真,对求解器步长很敏感。有一次为了加快仿真速度,把步长从微秒级放大到几十微秒,结果CT饱和后的二次电流波形全是毛刺,傅里叶滤波后的差动电流也出现了大幅波动,保护误动。一开始以为是CT模型参数问题,调了半天没效果,后来把步长改回微秒级,波形立刻干净了。
这个坑说明,含非线性元件的仿真,步长不能随便放大。CT饱和过程涉及磁通的快速变化,步长太大,积分误差累积,波形就会失真。建议含CT饱和模型的仿真,步长设在1到10微秒之间,用变步长求解器,让求解器在饱和时刻自动调整。
7.4 二次谐波含量算对了,但闭锁逻辑的“与门”接错了
最后一个坑比较“低级”但很典型。涌流工况下,二次谐波含量算出来是对的,确实超过了制动阈值,但保护还是动作了。查了半天,发现闭锁逻辑里,二次谐波制动信号和差动动作信号之间的逻辑门接错了——本该是“制动信号有效时闭锁动作”,结果接成了“制动信号有效时允许动作”。这种逻辑连接错误在图形化建模里特别容易发生,因为模块一多,连线就看花眼了。
教训是:保护逻辑的每个逻辑门,接完之后要单独做真值表验证,输入各种组合看输出对不对。别嫌麻烦,这一步比事后排查省时间得多。
8. 写在最后:仿真做对了,现场才敢信
差动保护仿真这件事,说到底是在用模型复现一个物理比较过程。模型里每一个环节的误差,都会累积到最后的差动电流上,而差动保护的定值又偏偏很敏感,所以对建模精度的要求比一般仿真高得多。我这些年做下来最大的体会是:仿真结果“看起来对”和“真的对”之间,差的就是那些不起眼的细节——一个折算系数、一个窗长、一个步长、一个逻辑门。这些细节在教材里往往一笔带过,但在实操里每一个都能让你卡上半天。
所以我的建议是,每搭完一个环节,就单独验证这个环节的输入输出关系,别等整个模型搭完再一起调。分环节验证虽然前期慢一点,但后期排查的时候,你能快速定位问题出在哪个环节,而不是面对一个黑箱模型瞎猜。另外,仿真参数和实际装置的定值一定要对齐,别拿一个“差不多”的参数去跑,差动保护这个场景里,“差不多”往往就是“差很多”。
如果你也在做变压器差动保护的仿真,遇到什么奇怪的现象,欢迎一起交流。很多时候问题就出在某个你根本没注意到的建模细节上,说出来大家一起看看,往往比自己闷头查快得多。