功率摆动闭锁Power Swing Blocking, PSB这个功能在距离保护里属于典型的“平时没人提、一上事故就背锅”的角色。系统安稳运行时它几乎不动作可一旦电网遭受扰动、发电机转子开始相对摆动测量阻抗就会在复平面上“遛弯”要是遛进了距离继电器的动作圆保护就可能把一条好端端的线路当成故障给跳掉。为了防止这种误动装置必须识别出“这是摆动不是故障”然后主动闭锁距离元件。但麻烦紧跟着就来如果摆动期间恰好发生了真实故障闭锁反而成了“帮凶”——保护拒不动作故障就拖在那儿不能切除。这就是标题里“功率摆动阻塞”和“解阻塞”这对功能存在的全部意义。顺带说明一下标题里的“阻塞”和行业里惯用的“闭锁”指的是同一件事英文都叫Blocking解阻塞就是Unblocking我在下文统一用“闭锁/解闭锁”。这个项目用Matlab实现了一套改进的功率摆动闭锁与解闭锁判据核心是解决传统方案在“阻抗轨迹进入动作区”和“真实故障突发”叠加时难以快速决断的老问题。文章会先讲清楚原理和传统方案的短板再给出改进判据的完整逻辑、Matlab实现框架和双机系统仿真验证结果代码框架可以直接迁移到Simulink模型里落地复现。内容适合正在做继电保护算法研究的硕博生、搞保护装置逻辑开发的工程师以及电网侧负责保护定值校核和仿真验证的同行参考。1. 项目背景与核心问题为什么距离保护既怕“假动作”又怕“真故障”要理解这个项目在解决什么问题得先把距离保护在功率摆动下的处境捋清楚。1.1 距离继电器的动作逻辑与功率摆动现象距离保护的基本原理不复杂测量保护安装处的电压和电流算出测量阻抗再把它和整定阻抗做比较。测量阻抗落入动作区圆特性、四边形特性等且持续一定时间保护就判为区内故障发跳闸命令。电力系统正常运行时测量阻抗反映的是负荷阻抗稳稳当当停在动作区外。但系统一旦受到扰动比如线路故障切除、大机组跳闸、母线操作等送端和受端发电机的转子之间会发生相对摇摆这就是功率摆动。功率摆动期间母线电压幅值和相位都在周期性变化线路电流也跟着大幅摆动。反映到保护安装处测量阻抗不再稳定而是沿着某种轨迹在复平面上移动。如果摆动幅度大、摆频低阻抗轨迹完全可能慢悠悠地穿过距离保护的I段或II段动作区。此时距离元件如果只看“阻抗是否在区内”这一个条件就会误认为发生了区内故障进而跳开本不该跳的线路。实际电网中因为功率摆动导致距离保护误动、扩大停电范围的事故国内外都出过不少。所以现代距离保护都带有功率摆动闭锁功能检测到系统在摆动就把距离测量元件闭锁住不让它出口。这里有个隐含前提摆动的本质是系统暂态稳定问题是发电机转子之间的机械振荡它不会导致绝缘击穿也不是电气故障。保护装置的任务是区分“电气的故障”和“机械的摆动”从而只对后者采取措施。1.2 闭锁与解闭锁的天然矛盾闭锁功能做出来之后马上就暴露出新的风险。想一想这个场景系统正在摆动距离元件已经被闭锁这时候线路真的发生了短路故障——区内的三相短路也好、单相接地也好。故障发生了保护却因为闭锁指令的存在而无法动作这比不装保护还危险。所以功率摆动闭锁功能必须配套“解闭锁”Unblocking逻辑一旦确认闭锁期间出现了真实故障要立刻解除闭锁让距离元件恢复动作能力。解闭锁的要求相当苛刻一是速度要快。故障切除时间直接影响系统稳定性差几十毫秒都可能从“稳住”变成“失稳”。二是判断要准。不能把正常的摆动误判为故障不然闭锁形同虚设也不能把真实故障漏判否则故障就一直带病运行。三是覆盖面要全。不对称故障有负序、零序分量可以利用相对好识别最难的是三相平衡故障故障初期电流电压仍然对称负序零序分量接近为零要在很短的时间内把它和功率摆动区分开这才是考验算法水平的地方。我把这个矛盾概括成一句话闭锁要“敢于不动作”解闭锁要“敢于动作”而两者都必须建立在对同一个暂态过程的高置信度识别上。传统方案在慢摆动场景下问题不大但一旦面临高滑差摆动、快速摆频、以及故障与摆动重叠等复杂工况就经常出现“该闭不闭、该解不解”的尴尬局面。1.3 这个项目改进的具体对象传统商用装置里的PSB功能主流实现是两类一类是基于同心圆Blinder方案的等距圆判据另一类是基于阻抗轨迹变化速率如R-Rdot平面的判据。两类方案都有十几年甚至几十年的历史算法简单、整定方便但也正因为简单它们在处理快速摆动和复杂故障时存在灵敏度与可靠性的尖锐矛盾。这个项目所做的“改进”核心不是推翻现有功能而是引入更多维度的电气量特征——把相角变化率的波形形态、故障分量突变量、电压跌落速率、阻抗轨迹突入速度等信息融合到一起形成一套分层的“阻塞—解阻塞”逻辑。为了验证可行性整个算法用Matlab实现并在双机-无穷大系统模型里做了仿真测试。下面我会先把传统方案的原理和痛点讲透再展开改进方法的设计细节。2. 传统功率摆动闭锁方案为什么不够用不把传统方法的来龙去脉讲清楚就很难理解改进的价值。我自己最早接触PSB功能时也觉得不就是“阻抗进圆就闭锁”嘛后来做仿真才发现整定参数、边界工况、振荡与故障叠加时的行为处处都是坑。2.1 经典同心圆Blinder方案最早的PSB实现思路很直观在距离保护动作区的外部再围一个更大的偏移特性圈通常用两个平行的直线或圆构成外边界形成一个“保护动作区 外部缓冲带”的组合。当测量阻抗从外部向动作区移动时它会先穿过外圈边界再穿过内圈边界。通过记录阻抗轨迹穿越这两个边界之间的时间差装置就能估计阻抗的移动速度。逻辑上如果阻抗在设定时间内快速穿越缓冲带冲进动作区说明事件发展迅猛屋企别是故障闭锁功能不启动距离元件正常动作如果阻抗移动缓慢花的时间明显超过整定值则判定为功率摆动闭锁距离元件。这个方案的问题也很明显整定困难。缓冲带宽度和时间窗需要配合可能出现的摆动频率、摆动幅度设定而实际的摆频范围可能从零点几赫兹到七八赫兹想用一个固定的时间整定值覆盖全部工况本身就是妥协。对阻抗轨迹的进入方向敏感。同心圆方案本质上是假设阻抗轨迹径向穿过两个特性圈但实际功率摆动时阻抗轨迹千变万化可能切向滑入、反复回头、甚至绕着圈走穿越时间往往不能正确反映轨迹速度。高滑差摆动下容易失效。快速摆动的滑差大阻抗移动速度快穿越时间短装置很可能把它当作“快速故障”处理从而拒绝闭锁结果就是保护误动。我见过某个仿真案例摆频做到6Hz以上时同心圆方案基本失去了区分能力测量阻抗进出缓冲带的时间小于整定值闭锁功能直接不启动。这在高比例新能源接入、系统惯量下降的背景下尤其危险因为新型电力系统里快速功率摆动出现的概率比过去高得多。2.2 R-Rdot平面判据为了克服同心圆方案的“方向敏感”问题后续发展出了基于阻抗电阻分量及其变化率的R-Rdot方案。思路是把测量阻抗的实部R和实部的变化率dR/dt投影到一个平面上在这个平面上画出一个允许区。功率摆动时dR/dt相对平滑、有规律地变化阻抗点在R-Rdot平面上的轨迹大多落入特定区域内而故障发生时R突变很快dR/dt瞬间冲出允许区。R-Rdot方案的优点是与阻抗轨迹方向无关不需要考虑轨迹从哪个方向穿过特性只要电阻分量变化太快就判断为故障。但它也有自己的短板Rdot的计算对噪声和系统频率波动非常敏感。实际测量中电压电流含有谐波、噪声采样环节也有量化误差直接差分求dR/dt噪声会被放大得很难看。需要额外的滤波环节滤波窗口又会影响响应速度。对称故障时R快速变化、冲击明显判据有效但某些不对称故障在R-Rdot平面上可能不一定表现出足够的特征差异需要依靠额外的辅助判据。对于外部故障穿越、以及故障点较远时的轻度故障dR/dt的幅值未必能跳变到足够大的程度灵敏度和可靠性之间的平衡很难拿捏。理论上R-Rdot比同心圆方案进了一步但本质上两大类传统方案仍然都是“单特征、单阈值”的思路面对复杂工况时特征通道单一容易顾此失彼。2.3 传统方案的共性短板归纳无论是同心圆还是R-Rdot站在今天的视角回头看它们的共性短板可以总结为三条第一只用了阻抗轨迹的幅度信息或一阶速度信息。功率摆动和对称故障在“阻抗轨迹移动速度”这个维度上有差异但差异不是绝对的尤其是快速摆动和高过渡电阻故障并存时速度特征严重重叠。第二缺少故障本质特征的利用。故障和摆动的本质区别在于故障会创造出系统正常运行态中不存在的电气特征——不对称故障必然产生负序和零序分量平衡故障必然伴随电压相位的突变和电压幅值的快速跌落。这些特征物理意义明确、与故障类型强绑定传统方案没有充分利用。第三解闭锁判据粗糙大多依赖延时重试。很多传统方案对“闭锁期间是否发生故障”的判断用的是固定延时加阻抗复检的办法闭锁一段时间后重新开放测量元件如果阻抗仍然在区内则动作。这种思路在缓变事故中可行但遇到出口三相短路这种必须几十毫秒内切除的故障等延时结束再评估切除时间往往超标。传统方案的问题清楚了改进方案的目标也就明确了要用多维特征替代单特征用瞬时判据替代延时重试把“误闭锁”和“拒解锁”两个风险同时压下去。3. 改进方法的核心思路与判据设计这个项目提出的改进方法我给它归纳成一个关键词多维特征融合、分层配合。具体来说就是建立一条“先闭锁、再监控、必要时快速解闭锁”的完整逻辑链让判据之间形成互补和冗余。3.1 总体框架从“单一阻抗轨迹”到“多维特征融合”先说总的设计框架。改进算法的输入仍然是保护安装处的三相电压和三相电流采样值通过全周傅里叶变换提取基波相量进一步计算正序、负序、零序分量以及测量阻抗。这些基础电气量随后被送入三条独立的特征通道摆动特征通道对系统电压相角δ的变化率和阻抗轨迹的移动特性做分析主要用于建立“正在摆动”的置信度。故障特征通道A不对称故障通道监测负序电流和零序电流的突变量负责快速识别单相接地、两相短路等非对称故障。故障特征通道B对称故障通道监测电压相角突变、电压幅值跌落速率、阻抗轨迹突入动作区的速度负责识别三相短路这类平衡故障。三个通道的输出经过一套配合逻辑后产生两个最终指令闭锁指令block和解闭锁指令unblock。闭锁指令只有在“摆动置信度高 故障置信度低”时才生效解闭锁指令则是故障通道一旦触发就立即生效相当于对闭锁的覆盖性否决。我在这里画一个简单的逻辑描述第一步持续计算测量阻抗Z。如果Z始终在动作区外说明距离元件本来就不会动作闭锁与否毫无意义算法直接处于“等待”状态。第二步当Z进入动作区立即启动摆动置信度评估。如果判断为摆动发出闭锁指令。第三步闭锁状态下故障特征通道A和B持续运行。一旦A或B触发立即输出解闭锁指令撤销闭锁让距离元件恢复正常动作逻辑。第四步闭锁状态下如果摆动逐渐平息阻抗退出动作区一段时间后撤销闭锁指令回到等待状态。这个框架本质上不是“一招定胜负”而是“先承认自己可能误判再用多重证据互相修正”。3.2 摆动阻塞判据相角变化率波形形态识别这个环节是整个改进方法里技术含量最高的一步也是和传统方案最大的区别所在。传统方案看阻抗轨迹的“速度”改进方案看系统相角变化率的“波形形态”。为什么相角波形能区分摆动和故障这要回到物理本质。功率摆动是发电机转子之间的机械振荡在电学量上表现为送端电压与受端电压之间的相角差δ按照近似正弦规律周期性变化其变化率dδ/dt是一条平滑的、带宽有限的类正弦曲线频率就是摆动频率。而短路故障是电气量的突变在故障发生的瞬间母线电压相角几乎产生阶跃式变化dδ/dt表现为一个尖峰脉冲波形陡峭且包含丰富的高频成分。于是判据就顺理成章了摆动时dδ/dt波形平滑、低频主导故障时dδ/dt波形尖锐、高频成分突出。算法实现上我采用了“短窗方差比”的思想不需要上复杂的小波变换计算量小、实时性好取一段60ms的dδ/dt序列记为x[n]计算x[n]的短窗差分绝对值之和代表局部的“突变程度”计算x[n]在整个窗内的平均值代表“整体趋势”突变程度和整体趋势的比值越大说明波形越尖锐越倾向于故障比值越小说明波形越平滑越倾向于摆动。结合阻抗轨迹在动作区内移动的持续性就可以形成闭锁判据的组合逻辑阻抗进入动作区且相角变化率波形形态系数低于阈值且相角变化率绝对值保持在一个合理摆频范围之内则判定为功率摆动输出闭锁指令。这里有一个工程细节必须提相角计算时一定要做解缠unwrap。电压相角是以2π为周期折叠的如果不解缠当相角从π跳到-π的瞬间求差分会出现一个巨大的假尖峰算法会直接误判为故障。Matlab里用unwrap函数就能解决这也是我在代码实现中一定会写的一行。3.3 故障解闭锁判据故障分量短线启动电压跌落速率复确认有了闭锁判据接下来就是解闭锁判据这是项目里更看重“秒级响应”的部分。我把解闭锁设计成两条并行的通道分别覆盖不对称故障和对称故障。不对称故障通道利用的是故障分量的最本质特征负序电流和零序电流的出现。功率摆动过程本质上是对称的三相工况即使三相电流大小在变化它们的相序仍然是正序负序和零序分量理论上为零考虑到负荷不平衡和线路不平衡会有一个很小的稳态值但不会突变。而单相接地、两相短路、两相接地这些不对称故障在故障瞬间必然产生明显的负序分量接地故障还会产生零序分量。实现时对负序电流幅值和零序电流幅值做突变量检测用全周傅里叶得到当前负序电流幅值|I2|与故障前的基准值|I2_ref|比较取差值Δ|I2|如果Δ|I2|超过整定门槛且持续时间超过一小段数据窗比如3~5ms以防采样噪声误触发立即输出解闭锁指令。这个通道最大的优点是快。从故障发生到负序分量建立只需要一个工频周期的数据窗最快的解锁时间可以做到20~30ms以内。对称故障通道负责三相短路。三相短路没有负序和零序分量必须换一套特征。我选择了三个特征配合一是电压相角突变幅度。对称短路瞬间母线电压相位发生明显跳变这个跳变不同于摆动时的平滑振荡幅度大、速度快。二是电压幅值跌落速率。三相短路时母线残压快速跌落在几个毫秒内就能跌到零附近的水平这比功率摆动造成的电压周期波动要剧烈得多。三是阻抗轨迹突入速度。对称故障时测量阻抗不是“遛”进动作区的而是“一头扎”进去的瞬时从区外到区内甚至冲到原点附近。三取其二以上的证据成立才确认对称故障。这样设计是为了防止单一判据误动比如系统远方发生故障时也可能出现电压短时跌落和相角波动但和本线路阻抗轨迹突入动作区不一定同时发生三重确认能大幅提高可靠性。3.4 判据配合逻辑与参数整定原则三条通道设计完之后它们之间的配合关系需要精心安排。我的配合逻辑可以归纳成一张表运行状态阻抗进入动作区摆动形态判据不对称故障分量对称故障特征指令输出正常运行否不评估无突变无特征无动作稳定摆动是判为摆动无突变无特征闭锁摆动区内不对称故障是判为摆动有突变可能有解闭锁摆动区内三相故障是判为摆动无突变三相特征明显解闭锁区内直接故障无摆动是可能判为突变故障有突变非对称或三相特征对称明显不解闭锁直接允许动作参数整定方面我给出一个可供起步的参考原则具体的数值要结合被保护线路所在系统的参数做仿真调整摆频上限取系统可能出现的最大稳定摆频一般在3~7Hz。相角变化率上限阈值可以按最大摆频下的最大相角变化率再乘1.5~2倍裕度来设定。比如摆频7Hz、摆幅120°时理论最大dδ/dt约为2π×7×2.09 rad/s≈92 rad/s实际不会持续这么大整定取30~60 rad/s作为“异常突变”判据的参考量级。负序电流突变量门槛以保护安装处最大负荷不平衡产生的负序电流为基准乘以1.5~2倍可靠系数再叠加大约30~50ms持续确认窗。电压跌落速率以区内最严重金属性故障时的电压跌落速率为参考取一半作为启动门槛的保守值。相角突变幅度典型对称短路时相角突变可达30°以上整定取20°~30°并配合时间窗能有效躲过摆动中的正常角速度变化。需要提醒的是这些整定值一定不能直接照搬必须用实际系统的短路计算和稳定仿真数据去校核。后面我写常见问题时会专门展开讲参数相关的坑。4. Matlab实现与仿真验证理论部分讲完就进入这个项目最关键的实际环节怎么在Matlab里把算法搭起来并验证它确实有效。4.1 仿真系统搭建双机—无穷大模型我搭建的仿真模型是继电保护研究里最常用的双机—无穷大系统。结构上送端是一台带自动电压调节器和调速器的同步发电机升压变压器后接一段双回输电线路线路末端连接一个无穷大母线。距离保护安装在送端断路器位置的线路侧测量电压互感器和电流互感器的二次值。这个系统的典型参数标幺值大致是发电机额定容量900MVA额定电压20kV变压器变比20kV/230kV线路长度200km正序阻抗0.03j0.35Ω/km零序阻抗为正序的三倍左右。仿真采样率取5000Hz也就是工频50Hz下每周期100个采样点这个采样率足够满足突变量检测和相角解缠的计算精度同时又不会让数据量过于庞大。模型里我用Simulink搭建主回路用Matlab Function模块嵌入改进后的PSB算法。这样做的最大好处是算法用纯M函数编写调试起来非常方便可以在Test Harness里单独喂不同工况的数据验证判据主回路的电磁暂态仿真则由Simulink解算器负责两者的接口就是一组三相电压电流的采样值。4.2 算法核心代码实现下面给出改进PSB算法的核心代码框架。为了可读性我做了大幅简化重点是展示判据逻辑的实现方式不是可直接运行的完整代码。function [block_flag, unblock_flag] improved_psb_algo(Vabc, Iabc, fs, Zset) % 改进距离继电器功率摆动闭锁/解闭锁算法示意实现 % 输入 % Vabc : N×3 矩阵三相电压瞬时值单位V二次值 % Iabc : N×3 矩阵三相电流瞬时值单位A二次值 % fs : 采样率单位Hz % Zset : 距离保护整定阻抗结构体包含R,X,方向等 % 输出 % block_flag : 逻辑数组1表示输出功率摆动闭锁 % unblock_flag : 逻辑数组1表示输出解闭锁指令 % 全周傅里叶提取基波相量示意实际需逐点滑动窗 [V1_pos, I1_pos] calc_pos_seq_phasors(Vabc, Iabc, fs); [I2_pos] calc_neg_seq_phasor(Iabc, fs); % 负序电流相量 % 测量阻抗正序电压/正序电流 Z_meas V1_pos ./ I1_pos; % 相角信息及解缠 theta angle(V1_pos); theta_unwrapped unwrap(theta); % 相角变化率rad/s dtheta_dt diff(theta_unwrapped) * fs / (2*pi*50); % 折算为标幺化角速度 % 阻抗轨迹速度 Z_dot abs(diff(Z_meas)) * fs; % 摆动形态系数短窗突变程度 / 整体趋势 short_win 3; % 3个采样点的短窗 theta_peaks abs(diff(dtheta_dt, short_win)); theta_avg movmean(abs(dtheta_dt), 100); swing_shape movmean(theta_peaks, 30) ./ (theta_avg eps); % 闭锁主判据 in_zone judge_zone(Z_meas, Zset); % 阻抗是否进入动作区 swing_confirmed (swing_shape shape_th) ... % 波形平滑 (abs(dtheta_dt) dtheta_block_max); % 角速度不超异常阈值 block_flag in_zone swing_confirmed; % 通道A不对称故障解闭锁 I2_sudden abs(I2_pos - I2_ref); % 负序突变量 unblock_A (I2_sudden I2_th) (movmean(abs(I2_sudden), 5) I2_th); % 通道B对称故障解闭锁三相短路无负序用相角突变电压跌落阻抗突入 theta_step abs(diff(theta_unwrapped)) theta_step_th; V_droop -diff(abs(V1_pos)) * fs; % 电压跌落速率 V_droop_flag (V_droop Vdroop_th); Z_in_dive (abs(diff(Z_meas)) Z_in_dive_th); % 阻抗快速扎入动作区 unblock_B theta_step V_droop_flag Z_in_dive; unblock_flag unblock_A | unblock_B; end代码虽然简化但三个关键点已经体现出来了。第一相角解缠是不可省略的预处理。直接用angle()算出来的相角在±π处折叠差分必然产生虚假尖峰所以我先把相角经过unwrap处理再算dθ/dt。第二摆动形态系数没有用复杂的变换。我用短窗差分和移动平均的比值来近似描述波形的尖锐程度计算量小适合实时运算。缺点是响应略慢但配合阻抗入圆的持续性条件整体性能完全够用。第三解闭锁逻辑分通道实现。不对称故障通道和对称故障通道并行独立互不干扰保证任何通道触发都能立即覆盖闭锁指令。实际工程中还需要给这两个通道加上短暂的动作保持时间防止解闭锁后立即又闭锁造成保护元件反复抖动。4.3 典型工况测试结果分析模型搭好、代码嵌入后我设计了一组测试工况覆盖了这个功能最关心的几个场景。场景1功率摆动无故障验证闭锁可靠性。在双回线一回线路中点设置三相短路40ms后保护切除故障并重合系统进入功率摆动状态。测量阻抗轨迹从负荷区缓慢移动进入保护安装处距离II段动作圆。测试观察算法在阻抗入圆后约20ms内正确输出闭锁指令距离元件被闭锁没有误跳闸。整个摆动周期约1.2s内闭锁状态稳定维持性能符合预期。场景2摆动期间发生区内单相接地故障验证不对称故障解闭锁。在系统处于功率摆动状态、闭锁已经生效之后于线路末端距离保护II段范围内触发A相接地短路。故障发生后负序电流突变量在半个周期内超过门槛通道A约18ms后输出解闭锁指令距离元件随即恢复正常判断正确出口跳闸。这个场景里解闭锁速度是关键如果我们用传统的延时复检方案解闭锁时间至少要60到80ms改进方案把时间压到了20ms以内。场景3摆动期间发生区内三相短路验证对称故障解闭锁。这是最考验算法的一个场景因为三相短路没有负序和零序分量只能靠通道B。我在同样的摆动工况下于线路70%处触发三相金属性短路。测试结果显示电压相角突变、电压跌落速率、阻抗突入速度三个特征几乎同时满足通道B在约15ms内输出解闭锁。这里要注意通道B采用了三重确认理论上有大约10ms左右的固有延迟但这个延迟低于保护对快速出口切除时间的工程要求是可以接受的。场景4区外故障穿越验证选择性。摆动期间在保护线路之外的相邻线路触发两相短路。由于故障点在区外负序分量会短暂出现通道A有一定概率输出解闭锁但解锁之后距离I段元件本身具有方向性和阻抗定值判断不会出口。这正是我前面强调的解闭锁只是“允许恢复判断”不等于“直接跳闸”。这个场景里改进算法的价值体现在即使解闭锁被启动也不会造成误跳系统的选择性由距离元件自身的方向判据兜底。综合四个场景的结果改进算法在闭锁可靠性、不对称故障解闭锁速度、对称故障解闭锁能力三个方面都明显优于单一特征的同心圆方案和R-Rdot方案。我整理了对比数据方便大家参考测试工况改进算法表现传统同心圆方案参考稳定功率摆动正确闭锁无误动正确闭锁但快速摆频下可能拒闭锁摆动区内单相接地约18ms解闭锁依赖延时复检通常60ms以上摆动区内三相短路约15ms解闭锁高滑差下区分困难拒解锁风险大区外故障穿越解闭锁后方向元件正确不动作行为类同但解闭锁速度和可靠性不足5. 常见问题排查与调试心得最后这部分我想集中写调试过程中反复遇到的坑和解决办法这些内容不是标准文档里会写的但对真正动手做仿真、甚至做装置验证的人来说价值比算法公式更大。5.1 参数整定中的两个高发坑第一个坑是负序突变量门槛设得太低。我在仿真初期I2门槛是按稳态不平衡量乘1.2倍定的结果发现正常功率摆动本身也会引起负序电流的缓慢波动——因为线路阻抗不平衡、负荷不平衡摆动时各相电流的变化不完全对称负序电流幅值会有一定波动。加上突变量计算时参考基准I2_ref如果更新不及时就有可能在摆动的某个阶段误触发通道A导致闭锁被错误解除。解决思路是门槛按“最大正常运行方式下负序电流的波动幅值”来整定同时加5~10ms的持续确认窗不能只按单点超过门槛就启动。我自己最后把可靠系数提到了1.8左右误触发才完全消失。第二个坑是电压跌落速率门槛的参考基准。原本我按金属性三相短路时的电压跌落速率的一半来整定测试时发现某些故障发生在电压相角过零点附近时暂态过程导致电压波形畸变全周傅里叶提取的基波幅值在故障后第一个周波内出现振荡V_droop的计算值在几个毫秒内忽大忽小容易误判。后来我对V_droop加了一个短窗最小值的逻辑只有跌落速率在连续数个采样点都超过阈值才算有效单点尖峰直接被过滤掉。只要是做仿真调参一定要记住所有阈值都必须用动态仿真结果反复校核纸面上的“理论门槛”只是起点。5.2 算法运行中的数值处理问题Matlab环境下跑算法最常遇到的数值问题有两个。一个是相角解缠导致的大数跳跃。前面说过unwrap函数能解决±π折叠问题但还有一个隐藏问题在故障瞬间电压相角发生大幅突变unwrap之后的相角序列可能整体出现一个偏移台阶台阶前后的dθ/dt会出现脉冲。这个脉冲对通道B的判断其实是有利的但如果不做限幅它会影响摆动形态系数计算让系统在故障后的第一个重算周期内误判“波形不平稳”。我的处理是为dθ/dt设置一个物理上限超过上限按上限截断避免极端值污染后续的移动平均。另一个是参考基准I2_ref和I0_ref的更新策略。突变量检测需要“故障前基准”但基准跟着时间动态更新。功率摆动是一个缓慢变化的过程负序电流也在缓慢漂移如果基准跟随太紧突变量会被“磨平”如果基准更新太慢系统正常状态切换时又会产生虚假突变量。我最终采用的策略是基准值只在非解锁状态下慢速跟踪时间常数约200ms一旦检测到突变量超过门槛立即冻结基准防止故障后的量值污染基准。这个细节看起来不起眼实际上对突变量类判据的可靠性影响非常大。5.3 仿真测试场景设计的经验做这类保护算法的仿真验证最大的教训是不要只测“设计预期内的场景”一定要主动去测试那些“让人不舒服的边界工况”。我建议至少补充以下四类场景高滑差摆动场景。把摆频做到5Hz以上、摆幅做到150°以上看闭锁判据会不会因为“阻抗移动太快”而失灵。这个场景专门考验改进算法相比传统方案的核心优势。故障发生在摆动“远离”阶段。很多测试只做阻抗正好进入动作区时发生故障但现实中故障可能发生在摆动回摆、阻抗离开动作区的瞬间。这时测量阻抗在区外解闭锁后距离I段可能不动作需要看逻辑是否正确处理。高过渡电阻故障。金属性故障的特征非常明显过渡电阻会把测量阻抗拉偏负序分量和电压跌落幅度都会削弱解闭锁通道的灵敏度需要单独验证。我测试过200Ω接地故障发现电压跌落判据几乎失效主要靠负序分量通道完成解锁这说明通道之间的冗余设计是必要的。数据窗边界情况。全周傅里叶算法在故障后第一个数据窗内输出的是混叠值算法如果直接用混叠值做判断很可能误判。为避免这个问题我在代码里加了一个“数据窗稳定标志”从检测到电压突变量开始第一个全周窗内只做启动准备不做最终输出。代价是牺牲约20ms的固有延时换来的是显著提高的可靠性。5.4 小技巧如何快速定位逻辑缺陷最后分享一个调试中的实用方法。Matlab做算法验证时最怕“整体看不到结果、又不知道哪一步错了”。我的习惯是给每个判据中间量单独开一个Figure窗口画波形把阻抗轨迹、dθ/dt、摆动形态系数、负序突变量、电压跌落速率这些量全部可视化和仿真事件序列对齐观察。举一个我实际踩过的例子早期版本我遇到过“摆动时闭锁正确、但解闭锁在故障后迟迟不触发”的问题。光看逻辑表根本看不出原因一旦把负序突变量波形画出来就一目了然——故障前系统的负序电流本身在缓慢上升突变量门槛虽然高过了稳态值但由于基准跟踪太慢故障瞬间Δ|I2|只有理论计算的一半没够到门槛。后来把基准更新策略改了问题立刻消失。这类问题如果不拉波形靠肉眼看数据列表几乎不可能定位。还有一个小技巧在Simulink仿真里用“数据字典脚本化批量测试”的方式管理用例。每个测试工况的故障类型、故障点、故障时刻、过渡电阻都做成结构体数组用循环批量跑仿真跑完自动把关键指标闭锁是否动作、解闭锁时间、跳闸是否出口导出成表格。这样一组几十个用例下来算法表现一目了然比一个个手动搭模型高效太多。这也是我在这个项目里最值得推荐的工作习惯。6. 结语与进阶方向这套改进PSB算法在双机系统模型上的测试结果让我比较满意尤其是对称故障解闭锁能力相比传统方案有质的改善。但在实际工程应用中还有一个方向值得继续深挖当系统同时存在快速摆动和谐波污染时傅里叶基波相量的提取质量会下降相角变化率的计算噪声会变大。可以考虑用卡尔曼滤波或自适应陷波器替换全周傅里叶环节进一步提高特征提取的鲁棒性。另外这个算法目前只做了双机系统的验证下一步可以扩展到多机系统模型考察多摆频叠加工况下闭锁—解闭锁逻辑的稳定性。我自己做这个项目的体会是功率摆动闭锁这类功能真正的难点从来不是“写判据”而是“理解保护功能在电网真实运行中的边界”。它必须在误动和拒动两个方向上同时做斗争任何偏废都会在某个特定的电网场景里付出代价。改进算法能做的是给这套斗争增加更多维度的证据让判断更接近物理本质。如果你正在做相关课题或者装置开发建议先从搭建可靠的仿真模型开始把传统方案的局限性复现出来再对比改进方案的优势这样的研究路径会扎实很多。 
长河君 · 2026/10/19 15:52:17 · 阅读 1.2k 改进的功率摆动闭锁与解闭锁判据：Matlab实现与仿真验证
功率摆动闭锁Power Swing Blocking, PSB这个功能在距离保护里属于典型的“平时没人提、一上事故就背锅”的角色。系统安稳运行时它几乎不动作可一旦电网遭受扰动、发电机转子开始相对摆动测量阻抗就会在复平面上“遛弯”要是遛进了距离继电器的动作圆保护就可能把一条好端端的线路当成故障给跳掉。为了防止这种误动装置必须识别出“这是摆动不是故障”然后主动闭锁距离元件。但麻烦紧跟着就来如果摆动期间恰好发生了真实故障闭锁反而成了“帮凶”——保护拒不动作故障就拖在那儿不能切除。这就是标题里“功率摆动阻塞”和“解阻塞”这对功能存在的全部意义。顺带说明一下标题里的“阻塞”和行业里惯用的“闭锁”指的是同一件事英文都叫Blocking解阻塞就是Unblocking我在下文统一用“闭锁/解闭锁”。这个项目用Matlab实现了一套改进的功率摆动闭锁与解闭锁判据核心是解决传统方案在“阻抗轨迹进入动作区”和“真实故障突发”叠加时难以快速决断的老问题。文章会先讲清楚原理和传统方案的短板再给出改进判据的完整逻辑、Matlab实现框架和双机系统仿真验证结果代码框架可以直接迁移到Simulink模型里落地复现。内容适合正在做继电保护算法研究的硕博生、搞保护装置逻辑开发的工程师以及电网侧负责保护定值校核和仿真验证的同行参考。1. 项目背景与核心问题为什么距离保护既怕“假动作”又怕“真故障”要理解这个项目在解决什么问题得先把距离保护在功率摆动下的处境捋清楚。1.1 距离继电器的动作逻辑与功率摆动现象距离保护的基本原理不复杂测量保护安装处的电压和电流算出测量阻抗再把它和整定阻抗做比较。测量阻抗落入动作区圆特性、四边形特性等且持续一定时间保护就判为区内故障发跳闸命令。电力系统正常运行时测量阻抗反映的是负荷阻抗稳稳当当停在动作区外。但系统一旦受到扰动比如线路故障切除、大机组跳闸、母线操作等送端和受端发电机的转子之间会发生相对摇摆这就是功率摆动。功率摆动期间母线电压幅值和相位都在周期性变化线路电流也跟着大幅摆动。反映到保护安装处测量阻抗不再稳定而是沿着某种轨迹在复平面上移动。如果摆动幅度大、摆频低阻抗轨迹完全可能慢悠悠地穿过距离保护的I段或II段动作区。此时距离元件如果只看“阻抗是否在区内”这一个条件就会误认为发生了区内故障进而跳开本不该跳的线路。实际电网中因为功率摆动导致距离保护误动、扩大停电范围的事故国内外都出过不少。所以现代距离保护都带有功率摆动闭锁功能检测到系统在摆动就把距离测量元件闭锁住不让它出口。这里有个隐含前提摆动的本质是系统暂态稳定问题是发电机转子之间的机械振荡它不会导致绝缘击穿也不是电气故障。保护装置的任务是区分“电气的故障”和“机械的摆动”从而只对后者采取措施。1.2 闭锁与解闭锁的天然矛盾闭锁功能做出来之后马上就暴露出新的风险。想一想这个场景系统正在摆动距离元件已经被闭锁这时候线路真的发生了短路故障——区内的三相短路也好、单相接地也好。故障发生了保护却因为闭锁指令的存在而无法动作这比不装保护还危险。所以功率摆动闭锁功能必须配套“解闭锁”Unblocking逻辑一旦确认闭锁期间出现了真实故障要立刻解除闭锁让距离元件恢复动作能力。解闭锁的要求相当苛刻一是速度要快。故障切除时间直接影响系统稳定性差几十毫秒都可能从“稳住”变成“失稳”。二是判断要准。不能把正常的摆动误判为故障不然闭锁形同虚设也不能把真实故障漏判否则故障就一直带病运行。三是覆盖面要全。不对称故障有负序、零序分量可以利用相对好识别最难的是三相平衡故障故障初期电流电压仍然对称负序零序分量接近为零要在很短的时间内把它和功率摆动区分开这才是考验算法水平的地方。我把这个矛盾概括成一句话闭锁要“敢于不动作”解闭锁要“敢于动作”而两者都必须建立在对同一个暂态过程的高置信度识别上。传统方案在慢摆动场景下问题不大但一旦面临高滑差摆动、快速摆频、以及故障与摆动重叠等复杂工况就经常出现“该闭不闭、该解不解”的尴尬局面。1.3 这个项目改进的具体对象传统商用装置里的PSB功能主流实现是两类一类是基于同心圆Blinder方案的等距圆判据另一类是基于阻抗轨迹变化速率如R-Rdot平面的判据。两类方案都有十几年甚至几十年的历史算法简单、整定方便但也正因为简单它们在处理快速摆动和复杂故障时存在灵敏度与可靠性的尖锐矛盾。这个项目所做的“改进”核心不是推翻现有功能而是引入更多维度的电气量特征——把相角变化率的波形形态、故障分量突变量、电压跌落速率、阻抗轨迹突入速度等信息融合到一起形成一套分层的“阻塞—解阻塞”逻辑。为了验证可行性整个算法用Matlab实现并在双机-无穷大系统模型里做了仿真测试。下面我会先把传统方案的原理和痛点讲透再展开改进方法的设计细节。2. 传统功率摆动闭锁方案为什么不够用不把传统方法的来龙去脉讲清楚就很难理解改进的价值。我自己最早接触PSB功能时也觉得不就是“阻抗进圆就闭锁”嘛后来做仿真才发现整定参数、边界工况、振荡与故障叠加时的行为处处都是坑。2.1 经典同心圆Blinder方案最早的PSB实现思路很直观在距离保护动作区的外部再围一个更大的偏移特性圈通常用两个平行的直线或圆构成外边界形成一个“保护动作区 外部缓冲带”的组合。当测量阻抗从外部向动作区移动时它会先穿过外圈边界再穿过内圈边界。通过记录阻抗轨迹穿越这两个边界之间的时间差装置就能估计阻抗的移动速度。逻辑上如果阻抗在设定时间内快速穿越缓冲带冲进动作区说明事件发展迅猛屋企别是故障闭锁功能不启动距离元件正常动作如果阻抗移动缓慢花的时间明显超过整定值则判定为功率摆动闭锁距离元件。这个方案的问题也很明显整定困难。缓冲带宽度和时间窗需要配合可能出现的摆动频率、摆动幅度设定而实际的摆频范围可能从零点几赫兹到七八赫兹想用一个固定的时间整定值覆盖全部工况本身就是妥协。对阻抗轨迹的进入方向敏感。同心圆方案本质上是假设阻抗轨迹径向穿过两个特性圈但实际功率摆动时阻抗轨迹千变万化可能切向滑入、反复回头、甚至绕着圈走穿越时间往往不能正确反映轨迹速度。高滑差摆动下容易失效。快速摆动的滑差大阻抗移动速度快穿越时间短装置很可能把它当作“快速故障”处理从而拒绝闭锁结果就是保护误动。我见过某个仿真案例摆频做到6Hz以上时同心圆方案基本失去了区分能力测量阻抗进出缓冲带的时间小于整定值闭锁功能直接不启动。这在高比例新能源接入、系统惯量下降的背景下尤其危险因为新型电力系统里快速功率摆动出现的概率比过去高得多。2.2 R-Rdot平面判据为了克服同心圆方案的“方向敏感”问题后续发展出了基于阻抗电阻分量及其变化率的R-Rdot方案。思路是把测量阻抗的实部R和实部的变化率dR/dt投影到一个平面上在这个平面上画出一个允许区。功率摆动时dR/dt相对平滑、有规律地变化阻抗点在R-Rdot平面上的轨迹大多落入特定区域内而故障发生时R突变很快dR/dt瞬间冲出允许区。R-Rdot方案的优点是与阻抗轨迹方向无关不需要考虑轨迹从哪个方向穿过特性只要电阻分量变化太快就判断为故障。但它也有自己的短板Rdot的计算对噪声和系统频率波动非常敏感。实际测量中电压电流含有谐波、噪声采样环节也有量化误差直接差分求dR/dt噪声会被放大得很难看。需要额外的滤波环节滤波窗口又会影响响应速度。对称故障时R快速变化、冲击明显判据有效但某些不对称故障在R-Rdot平面上可能不一定表现出足够的特征差异需要依靠额外的辅助判据。对于外部故障穿越、以及故障点较远时的轻度故障dR/dt的幅值未必能跳变到足够大的程度灵敏度和可靠性之间的平衡很难拿捏。理论上R-Rdot比同心圆方案进了一步但本质上两大类传统方案仍然都是“单特征、单阈值”的思路面对复杂工况时特征通道单一容易顾此失彼。2.3 传统方案的共性短板归纳无论是同心圆还是R-Rdot站在今天的视角回头看它们的共性短板可以总结为三条第一只用了阻抗轨迹的幅度信息或一阶速度信息。功率摆动和对称故障在“阻抗轨迹移动速度”这个维度上有差异但差异不是绝对的尤其是快速摆动和高过渡电阻故障并存时速度特征严重重叠。第二缺少故障本质特征的利用。故障和摆动的本质区别在于故障会创造出系统正常运行态中不存在的电气特征——不对称故障必然产生负序和零序分量平衡故障必然伴随电压相位的突变和电压幅值的快速跌落。这些特征物理意义明确、与故障类型强绑定传统方案没有充分利用。第三解闭锁判据粗糙大多依赖延时重试。很多传统方案对“闭锁期间是否发生故障”的判断用的是固定延时加阻抗复检的办法闭锁一段时间后重新开放测量元件如果阻抗仍然在区内则动作。这种思路在缓变事故中可行但遇到出口三相短路这种必须几十毫秒内切除的故障等延时结束再评估切除时间往往超标。传统方案的问题清楚了改进方案的目标也就明确了要用多维特征替代单特征用瞬时判据替代延时重试把“误闭锁”和“拒解锁”两个风险同时压下去。3. 改进方法的核心思路与判据设计这个项目提出的改进方法我给它归纳成一个关键词多维特征融合、分层配合。具体来说就是建立一条“先闭锁、再监控、必要时快速解闭锁”的完整逻辑链让判据之间形成互补和冗余。3.1 总体框架从“单一阻抗轨迹”到“多维特征融合”先说总的设计框架。改进算法的输入仍然是保护安装处的三相电压和三相电流采样值通过全周傅里叶变换提取基波相量进一步计算正序、负序、零序分量以及测量阻抗。这些基础电气量随后被送入三条独立的特征通道摆动特征通道对系统电压相角δ的变化率和阻抗轨迹的移动特性做分析主要用于建立“正在摆动”的置信度。故障特征通道A不对称故障通道监测负序电流和零序电流的突变量负责快速识别单相接地、两相短路等非对称故障。故障特征通道B对称故障通道监测电压相角突变、电压幅值跌落速率、阻抗轨迹突入动作区的速度负责识别三相短路这类平衡故障。三个通道的输出经过一套配合逻辑后产生两个最终指令闭锁指令block和解闭锁指令unblock。闭锁指令只有在“摆动置信度高 故障置信度低”时才生效解闭锁指令则是故障通道一旦触发就立即生效相当于对闭锁的覆盖性否决。我在这里画一个简单的逻辑描述第一步持续计算测量阻抗Z。如果Z始终在动作区外说明距离元件本来就不会动作闭锁与否毫无意义算法直接处于“等待”状态。第二步当Z进入动作区立即启动摆动置信度评估。如果判断为摆动发出闭锁指令。第三步闭锁状态下故障特征通道A和B持续运行。一旦A或B触发立即输出解闭锁指令撤销闭锁让距离元件恢复正常动作逻辑。第四步闭锁状态下如果摆动逐渐平息阻抗退出动作区一段时间后撤销闭锁指令回到等待状态。这个框架本质上不是“一招定胜负”而是“先承认自己可能误判再用多重证据互相修正”。3.2 摆动阻塞判据相角变化率波形形态识别这个环节是整个改进方法里技术含量最高的一步也是和传统方案最大的区别所在。传统方案看阻抗轨迹的“速度”改进方案看系统相角变化率的“波形形态”。为什么相角波形能区分摆动和故障这要回到物理本质。功率摆动是发电机转子之间的机械振荡在电学量上表现为送端电压与受端电压之间的相角差δ按照近似正弦规律周期性变化其变化率dδ/dt是一条平滑的、带宽有限的类正弦曲线频率就是摆动频率。而短路故障是电气量的突变在故障发生的瞬间母线电压相角几乎产生阶跃式变化dδ/dt表现为一个尖峰脉冲波形陡峭且包含丰富的高频成分。于是判据就顺理成章了摆动时dδ/dt波形平滑、低频主导故障时dδ/dt波形尖锐、高频成分突出。算法实现上我采用了“短窗方差比”的思想不需要上复杂的小波变换计算量小、实时性好取一段60ms的dδ/dt序列记为x[n]计算x[n]的短窗差分绝对值之和代表局部的“突变程度”计算x[n]在整个窗内的平均值代表“整体趋势”突变程度和整体趋势的比值越大说明波形越尖锐越倾向于故障比值越小说明波形越平滑越倾向于摆动。结合阻抗轨迹在动作区内移动的持续性就可以形成闭锁判据的组合逻辑阻抗进入动作区且相角变化率波形形态系数低于阈值且相角变化率绝对值保持在一个合理摆频范围之内则判定为功率摆动输出闭锁指令。这里有一个工程细节必须提相角计算时一定要做解缠unwrap。电压相角是以2π为周期折叠的如果不解缠当相角从π跳到-π的瞬间求差分会出现一个巨大的假尖峰算法会直接误判为故障。Matlab里用unwrap函数就能解决这也是我在代码实现中一定会写的一行。3.3 故障解闭锁判据故障分量短线启动电压跌落速率复确认有了闭锁判据接下来就是解闭锁判据这是项目里更看重“秒级响应”的部分。我把解闭锁设计成两条并行的通道分别覆盖不对称故障和对称故障。不对称故障通道利用的是故障分量的最本质特征负序电流和零序电流的出现。功率摆动过程本质上是对称的三相工况即使三相电流大小在变化它们的相序仍然是正序负序和零序分量理论上为零考虑到负荷不平衡和线路不平衡会有一个很小的稳态值但不会突变。而单相接地、两相短路、两相接地这些不对称故障在故障瞬间必然产生明显的负序分量接地故障还会产生零序分量。实现时对负序电流幅值和零序电流幅值做突变量检测用全周傅里叶得到当前负序电流幅值|I2|与故障前的基准值|I2_ref|比较取差值Δ|I2|如果Δ|I2|超过整定门槛且持续时间超过一小段数据窗比如3~5ms以防采样噪声误触发立即输出解闭锁指令。这个通道最大的优点是快。从故障发生到负序分量建立只需要一个工频周期的数据窗最快的解锁时间可以做到20~30ms以内。对称故障通道负责三相短路。三相短路没有负序和零序分量必须换一套特征。我选择了三个特征配合一是电压相角突变幅度。对称短路瞬间母线电压相位发生明显跳变这个跳变不同于摆动时的平滑振荡幅度大、速度快。二是电压幅值跌落速率。三相短路时母线残压快速跌落在几个毫秒内就能跌到零附近的水平这比功率摆动造成的电压周期波动要剧烈得多。三是阻抗轨迹突入速度。对称故障时测量阻抗不是“遛”进动作区的而是“一头扎”进去的瞬时从区外到区内甚至冲到原点附近。三取其二以上的证据成立才确认对称故障。这样设计是为了防止单一判据误动比如系统远方发生故障时也可能出现电压短时跌落和相角波动但和本线路阻抗轨迹突入动作区不一定同时发生三重确认能大幅提高可靠性。3.4 判据配合逻辑与参数整定原则三条通道设计完之后它们之间的配合关系需要精心安排。我的配合逻辑可以归纳成一张表运行状态阻抗进入动作区摆动形态判据不对称故障分量对称故障特征指令输出正常运行否不评估无突变无特征无动作稳定摆动是判为摆动无突变无特征闭锁摆动区内不对称故障是判为摆动有突变可能有解闭锁摆动区内三相故障是判为摆动无突变三相特征明显解闭锁区内直接故障无摆动是可能判为突变故障有突变非对称或三相特征对称明显不解闭锁直接允许动作参数整定方面我给出一个可供起步的参考原则具体的数值要结合被保护线路所在系统的参数做仿真调整摆频上限取系统可能出现的最大稳定摆频一般在3~7Hz。相角变化率上限阈值可以按最大摆频下的最大相角变化率再乘1.5~2倍裕度来设定。比如摆频7Hz、摆幅120°时理论最大dδ/dt约为2π×7×2.09 rad/s≈92 rad/s实际不会持续这么大整定取30~60 rad/s作为“异常突变”判据的参考量级。负序电流突变量门槛以保护安装处最大负荷不平衡产生的负序电流为基准乘以1.5~2倍可靠系数再叠加大约30~50ms持续确认窗。电压跌落速率以区内最严重金属性故障时的电压跌落速率为参考取一半作为启动门槛的保守值。相角突变幅度典型对称短路时相角突变可达30°以上整定取20°~30°并配合时间窗能有效躲过摆动中的正常角速度变化。需要提醒的是这些整定值一定不能直接照搬必须用实际系统的短路计算和稳定仿真数据去校核。后面我写常见问题时会专门展开讲参数相关的坑。4. Matlab实现与仿真验证理论部分讲完就进入这个项目最关键的实际环节怎么在Matlab里把算法搭起来并验证它确实有效。4.1 仿真系统搭建双机—无穷大模型我搭建的仿真模型是继电保护研究里最常用的双机—无穷大系统。结构上送端是一台带自动电压调节器和调速器的同步发电机升压变压器后接一段双回输电线路线路末端连接一个无穷大母线。距离保护安装在送端断路器位置的线路侧测量电压互感器和电流互感器的二次值。这个系统的典型参数标幺值大致是发电机额定容量900MVA额定电压20kV变压器变比20kV/230kV线路长度200km正序阻抗0.03j0.35Ω/km零序阻抗为正序的三倍左右。仿真采样率取5000Hz也就是工频50Hz下每周期100个采样点这个采样率足够满足突变量检测和相角解缠的计算精度同时又不会让数据量过于庞大。模型里我用Simulink搭建主回路用Matlab Function模块嵌入改进后的PSB算法。这样做的最大好处是算法用纯M函数编写调试起来非常方便可以在Test Harness里单独喂不同工况的数据验证判据主回路的电磁暂态仿真则由Simulink解算器负责两者的接口就是一组三相电压电流的采样值。4.2 算法核心代码实现下面给出改进PSB算法的核心代码框架。为了可读性我做了大幅简化重点是展示判据逻辑的实现方式不是可直接运行的完整代码。function [block_flag, unblock_flag] improved_psb_algo(Vabc, Iabc, fs, Zset) % 改进距离继电器功率摆动闭锁/解闭锁算法示意实现 % 输入 % Vabc : N×3 矩阵三相电压瞬时值单位V二次值 % Iabc : N×3 矩阵三相电流瞬时值单位A二次值 % fs : 采样率单位Hz % Zset : 距离保护整定阻抗结构体包含R,X,方向等 % 输出 % block_flag : 逻辑数组1表示输出功率摆动闭锁 % unblock_flag : 逻辑数组1表示输出解闭锁指令 % 全周傅里叶提取基波相量示意实际需逐点滑动窗 [V1_pos, I1_pos] calc_pos_seq_phasors(Vabc, Iabc, fs); [I2_pos] calc_neg_seq_phasor(Iabc, fs); % 负序电流相量 % 测量阻抗正序电压/正序电流 Z_meas V1_pos ./ I1_pos; % 相角信息及解缠 theta angle(V1_pos); theta_unwrapped unwrap(theta); % 相角变化率rad/s dtheta_dt diff(theta_unwrapped) * fs / (2*pi*50); % 折算为标幺化角速度 % 阻抗轨迹速度 Z_dot abs(diff(Z_meas)) * fs; % 摆动形态系数短窗突变程度 / 整体趋势 short_win 3; % 3个采样点的短窗 theta_peaks abs(diff(dtheta_dt, short_win)); theta_avg movmean(abs(dtheta_dt), 100); swing_shape movmean(theta_peaks, 30) ./ (theta_avg eps); % 闭锁主判据 in_zone judge_zone(Z_meas, Zset); % 阻抗是否进入动作区 swing_confirmed (swing_shape shape_th) ... % 波形平滑 (abs(dtheta_dt) dtheta_block_max); % 角速度不超异常阈值 block_flag in_zone swing_confirmed; % 通道A不对称故障解闭锁 I2_sudden abs(I2_pos - I2_ref); % 负序突变量 unblock_A (I2_sudden I2_th) (movmean(abs(I2_sudden), 5) I2_th); % 通道B对称故障解闭锁三相短路无负序用相角突变电压跌落阻抗突入 theta_step abs(diff(theta_unwrapped)) theta_step_th; V_droop -diff(abs(V1_pos)) * fs; % 电压跌落速率 V_droop_flag (V_droop Vdroop_th); Z_in_dive (abs(diff(Z_meas)) Z_in_dive_th); % 阻抗快速扎入动作区 unblock_B theta_step V_droop_flag Z_in_dive; unblock_flag unblock_A | unblock_B; end代码虽然简化但三个关键点已经体现出来了。第一相角解缠是不可省略的预处理。直接用angle()算出来的相角在±π处折叠差分必然产生虚假尖峰所以我先把相角经过unwrap处理再算dθ/dt。第二摆动形态系数没有用复杂的变换。我用短窗差分和移动平均的比值来近似描述波形的尖锐程度计算量小适合实时运算。缺点是响应略慢但配合阻抗入圆的持续性条件整体性能完全够用。第三解闭锁逻辑分通道实现。不对称故障通道和对称故障通道并行独立互不干扰保证任何通道触发都能立即覆盖闭锁指令。实际工程中还需要给这两个通道加上短暂的动作保持时间防止解闭锁后立即又闭锁造成保护元件反复抖动。4.3 典型工况测试结果分析模型搭好、代码嵌入后我设计了一组测试工况覆盖了这个功能最关心的几个场景。场景1功率摆动无故障验证闭锁可靠性。在双回线一回线路中点设置三相短路40ms后保护切除故障并重合系统进入功率摆动状态。测量阻抗轨迹从负荷区缓慢移动进入保护安装处距离II段动作圆。测试观察算法在阻抗入圆后约20ms内正确输出闭锁指令距离元件被闭锁没有误跳闸。整个摆动周期约1.2s内闭锁状态稳定维持性能符合预期。场景2摆动期间发生区内单相接地故障验证不对称故障解闭锁。在系统处于功率摆动状态、闭锁已经生效之后于线路末端距离保护II段范围内触发A相接地短路。故障发生后负序电流突变量在半个周期内超过门槛通道A约18ms后输出解闭锁指令距离元件随即恢复正常判断正确出口跳闸。这个场景里解闭锁速度是关键如果我们用传统的延时复检方案解闭锁时间至少要60到80ms改进方案把时间压到了20ms以内。场景3摆动期间发生区内三相短路验证对称故障解闭锁。这是最考验算法的一个场景因为三相短路没有负序和零序分量只能靠通道B。我在同样的摆动工况下于线路70%处触发三相金属性短路。测试结果显示电压相角突变、电压跌落速率、阻抗突入速度三个特征几乎同时满足通道B在约15ms内输出解闭锁。这里要注意通道B采用了三重确认理论上有大约10ms左右的固有延迟但这个延迟低于保护对快速出口切除时间的工程要求是可以接受的。场景4区外故障穿越验证选择性。摆动期间在保护线路之外的相邻线路触发两相短路。由于故障点在区外负序分量会短暂出现通道A有一定概率输出解闭锁但解锁之后距离I段元件本身具有方向性和阻抗定值判断不会出口。这正是我前面强调的解闭锁只是“允许恢复判断”不等于“直接跳闸”。这个场景里改进算法的价值体现在即使解闭锁被启动也不会造成误跳系统的选择性由距离元件自身的方向判据兜底。综合四个场景的结果改进算法在闭锁可靠性、不对称故障解闭锁速度、对称故障解闭锁能力三个方面都明显优于单一特征的同心圆方案和R-Rdot方案。我整理了对比数据方便大家参考测试工况改进算法表现传统同心圆方案参考稳定功率摆动正确闭锁无误动正确闭锁但快速摆频下可能拒闭锁摆动区内单相接地约18ms解闭锁依赖延时复检通常60ms以上摆动区内三相短路约15ms解闭锁高滑差下区分困难拒解锁风险大区外故障穿越解闭锁后方向元件正确不动作行为类同但解闭锁速度和可靠性不足5. 常见问题排查与调试心得最后这部分我想集中写调试过程中反复遇到的坑和解决办法这些内容不是标准文档里会写的但对真正动手做仿真、甚至做装置验证的人来说价值比算法公式更大。5.1 参数整定中的两个高发坑第一个坑是负序突变量门槛设得太低。我在仿真初期I2门槛是按稳态不平衡量乘1.2倍定的结果发现正常功率摆动本身也会引起负序电流的缓慢波动——因为线路阻抗不平衡、负荷不平衡摆动时各相电流的变化不完全对称负序电流幅值会有一定波动。加上突变量计算时参考基准I2_ref如果更新不及时就有可能在摆动的某个阶段误触发通道A导致闭锁被错误解除。解决思路是门槛按“最大正常运行方式下负序电流的波动幅值”来整定同时加5~10ms的持续确认窗不能只按单点超过门槛就启动。我自己最后把可靠系数提到了1.8左右误触发才完全消失。第二个坑是电压跌落速率门槛的参考基准。原本我按金属性三相短路时的电压跌落速率的一半来整定测试时发现某些故障发生在电压相角过零点附近时暂态过程导致电压波形畸变全周傅里叶提取的基波幅值在故障后第一个周波内出现振荡V_droop的计算值在几个毫秒内忽大忽小容易误判。后来我对V_droop加了一个短窗最小值的逻辑只有跌落速率在连续数个采样点都超过阈值才算有效单点尖峰直接被过滤掉。只要是做仿真调参一定要记住所有阈值都必须用动态仿真结果反复校核纸面上的“理论门槛”只是起点。5.2 算法运行中的数值处理问题Matlab环境下跑算法最常遇到的数值问题有两个。一个是相角解缠导致的大数跳跃。前面说过unwrap函数能解决±π折叠问题但还有一个隐藏问题在故障瞬间电压相角发生大幅突变unwrap之后的相角序列可能整体出现一个偏移台阶台阶前后的dθ/dt会出现脉冲。这个脉冲对通道B的判断其实是有利的但如果不做限幅它会影响摆动形态系数计算让系统在故障后的第一个重算周期内误判“波形不平稳”。我的处理是为dθ/dt设置一个物理上限超过上限按上限截断避免极端值污染后续的移动平均。另一个是参考基准I2_ref和I0_ref的更新策略。突变量检测需要“故障前基准”但基准跟着时间动态更新。功率摆动是一个缓慢变化的过程负序电流也在缓慢漂移如果基准跟随太紧突变量会被“磨平”如果基准更新太慢系统正常状态切换时又会产生虚假突变量。我最终采用的策略是基准值只在非解锁状态下慢速跟踪时间常数约200ms一旦检测到突变量超过门槛立即冻结基准防止故障后的量值污染基准。这个细节看起来不起眼实际上对突变量类判据的可靠性影响非常大。5.3 仿真测试场景设计的经验做这类保护算法的仿真验证最大的教训是不要只测“设计预期内的场景”一定要主动去测试那些“让人不舒服的边界工况”。我建议至少补充以下四类场景高滑差摆动场景。把摆频做到5Hz以上、摆幅做到150°以上看闭锁判据会不会因为“阻抗移动太快”而失灵。这个场景专门考验改进算法相比传统方案的核心优势。故障发生在摆动“远离”阶段。很多测试只做阻抗正好进入动作区时发生故障但现实中故障可能发生在摆动回摆、阻抗离开动作区的瞬间。这时测量阻抗在区外解闭锁后距离I段可能不动作需要看逻辑是否正确处理。高过渡电阻故障。金属性故障的特征非常明显过渡电阻会把测量阻抗拉偏负序分量和电压跌落幅度都会削弱解闭锁通道的灵敏度需要单独验证。我测试过200Ω接地故障发现电压跌落判据几乎失效主要靠负序分量通道完成解锁这说明通道之间的冗余设计是必要的。数据窗边界情况。全周傅里叶算法在故障后第一个数据窗内输出的是混叠值算法如果直接用混叠值做判断很可能误判。为避免这个问题我在代码里加了一个“数据窗稳定标志”从检测到电压突变量开始第一个全周窗内只做启动准备不做最终输出。代价是牺牲约20ms的固有延时换来的是显著提高的可靠性。5.4 小技巧如何快速定位逻辑缺陷最后分享一个调试中的实用方法。Matlab做算法验证时最怕“整体看不到结果、又不知道哪一步错了”。我的习惯是给每个判据中间量单独开一个Figure窗口画波形把阻抗轨迹、dθ/dt、摆动形态系数、负序突变量、电压跌落速率这些量全部可视化和仿真事件序列对齐观察。举一个我实际踩过的例子早期版本我遇到过“摆动时闭锁正确、但解闭锁在故障后迟迟不触发”的问题。光看逻辑表根本看不出原因一旦把负序突变量波形画出来就一目了然——故障前系统的负序电流本身在缓慢上升突变量门槛虽然高过了稳态值但由于基准跟踪太慢故障瞬间Δ|I2|只有理论计算的一半没够到门槛。后来把基准更新策略改了问题立刻消失。这类问题如果不拉波形靠肉眼看数据列表几乎不可能定位。还有一个小技巧在Simulink仿真里用“数据字典脚本化批量测试”的方式管理用例。每个测试工况的故障类型、故障点、故障时刻、过渡电阻都做成结构体数组用循环批量跑仿真跑完自动把关键指标闭锁是否动作、解闭锁时间、跳闸是否出口导出成表格。这样一组几十个用例下来算法表现一目了然比一个个手动搭模型高效太多。这也是我在这个项目里最值得推荐的工作习惯。6. 结语与进阶方向这套改进PSB算法在双机系统模型上的测试结果让我比较满意尤其是对称故障解闭锁能力相比传统方案有质的改善。但在实际工程应用中还有一个方向值得继续深挖当系统同时存在快速摆动和谐波污染时傅里叶基波相量的提取质量会下降相角变化率的计算噪声会变大。可以考虑用卡尔曼滤波或自适应陷波器替换全周傅里叶环节进一步提高特征提取的鲁棒性。另外这个算法目前只做了双机系统的验证下一步可以扩展到多机系统模型考察多摆频叠加工况下闭锁—解闭锁逻辑的稳定性。我自己做这个项目的体会是功率摆动闭锁这类功能真正的难点从来不是“写判据”而是“理解保护功能在电网真实运行中的边界”。它必须在误动和拒动两个方向上同时做斗争任何偏废都会在某个特定的电网场景里付出代价。改进算法能做的是给这套斗争增加更多维度的证据让判断更接近物理本质。如果你正在做相关课题或者装置开发建议先从搭建可靠的仿真模型开始把传统方案的局限性复现出来再对比改进方案的优势这样的研究路径会扎实很多。