简介一份基于MATLAB与PSCAD4.2实现电力系统距离保护仿真的技术资料适合继电保护方向学生、电力系统仿真初学者及从事输电线路保护研究的工程人员。内容从距离保护三段式原理与整定计算入手逐步给出接地距离保护、相间距离保护的主电路模型和控制电路模型搭建方法并针对单相接地、相间短路等典型故障进行仿真结合电压、电流、功率波形分析保护动作过程便于理解阻抗继电器与断路器联动逻辑。资源为1个doc文档压缩包约152KB文字与图表完整可配合PSCAD4.2实际操作参考。该文档已被422人学习适合用来入门距离保护仿真、撰写课程设计或作为建模排错的参考。1. 输电线路距离保护的核心用测量阻抗替代过流判据距离保护的名字里带着距离实际测的却是阻抗这个反差点正是它在复杂电网中比电流保护更可靠的原因。在PSCAD4.2里复现这个逻辑并不复杂搭一个单端电源两级线路模型装上三段式距离保护再分别设置A相接地和AB相间短路观察断路器在故障发生后的动作时序。我按自己拆这个模型时的顺序写这篇博文从整定计算、元件参数到波形判读每一步尽量给出能直接抄的参数和常见坑。适合做继电保护课程设计、刚接触PSCAD/EMTDC的电气专业学生以及想快速验证保护逻辑的调试工程师。2. 三段式距离保护整定计算与接地/相间测量阻抗的差异2.1 为什么是阻抗而不是电流线路发生短路时流过保护的短路电流大小受系统运行方式影响很大系统阻抗小的时候短路电流大阻抗大的时候电流小仅靠电流幅值很难划定稳定的保护边界。而保护安装处的电压与电流之比在故障点电气距离内近似正比于线路长度这就是距离保护不受运行方式影响的由来。PSCAD4.2里测量阻抗的计算也是这个思路正常运行时测量阻抗落在负荷区故障后迅速进入保护动作区断路器在几十毫秒内切开故障。2.2 三段式整定边界与动作时间配合距离保护要做到选择性、速动性和灵敏性兼顾靠的是三段配合。I段不能伸出本线路末端否则下一级线路出口故障时本段会误动II段要覆盖本线路全长并延伸到相邻线路的一部分作为I段末端死区的补充III段则要躲开正常负荷阻抗作为I、II段的后备。以下是常见整定原则PSCAD4.2仿真模型里直接按这个思路填阻抗圆参数段别保护范围整定原则动作时间I段本线路全长的80%~90%Z_1set K_rel1 * Z_ABK_rel1取0.8~0.9瞬时II段本线路全长 相邻线路一部分Z_2set K_rel2 * (Z_AB K_rel2 * Z_BC)0.3~0.5sIII段相邻线路及更远后备躲最小负荷阻抗折算自启动系数和返回系数阶梯配合I段的0.8~0.9可靠系数就是为了躲开本线路末端经过渡电阻短路和测量误差的影响。II段的整定里括号内的第二项是相邻线路的I段保护范围两个K_rel2相乘表示两个不确定环节的叠加。III段不按线路阻抗算而是按最小负荷阻抗打折原因很简单正常运行时负荷也能让测量阻抗变小如果III段定值不够保守重负荷时保护就会误动。2.3 接地与相间测量阻抗的公式差异距离保护按测量对象分成两类。接地距离保护测量的是相电压与带零序补偿的相电流之比因为单相接地时故障电流里含有零序分量直接用相电流会算出一个偏大的测量阻抗导致保护范围缩短。工程中常见的补偿形式为Z_A U_A / (I_A 3 * K0 * I_0)其中K0是零序补偿系数理论上与线路零序阻抗和正序阻抗有关PSCAD4.2模型里通常按实数近似输入。相间距离保护则测量线电压与线电流之比不涉及零序补偿Z_AB (U_A - U_B) / (I_A - I_B)两种公式的物理含义不同对应在控制回路里就要分别搭线对地阻抗和相间阻抗两个测量模块。仿真时如果只配置了接地保护AB相间短路故障时测量阻抗不会被正确计算断路器不会动作这是新手最容易漏掉的一点。2.4 用MATLAB先把整定值算清楚进PSCAD之前我会先用MATLAB把三段整定值算出来这样后面填阻抗圆参数时心里有数也方便把整定结果和仿真波形对照。下面这个脚本是常见做法Z_AB和Z_BC的数值需要按你模型里的线路参数替换% distance_relay_setting.m % 距离保护三段整定阻抗计算, 单位统一为欧姆 Z_AB 10 1j*35; % 线路AB正序阻抗 Z_BC 8 1j*28; % 相邻线路BC正序阻抗 K_rel1 0.85; % I段可靠系数, 0.8~0.9 K_rel2 0.8; % II段可靠系数, 考虑分支影响 K_ss 1.5; % 电动机自启动系数 K_re 1.2; % 继电器返回系数 U_ph 230e3/sqrt(3); % 最低运行相电压(V) I_max 600; % 最大负荷电流(A) % I段: 只保护本线路的一部分 Z_set1 K_rel1 * Z_AB; % II段: 本线路全长 相邻线路I段范围 Z_set2 K_rel2 * (Z_AB K_rel2 * Z_BC); % III段: 按躲最小负荷阻抗整定, 0.9为电压系数 Z_load_min U_ph / I_max; Z_set3 0.9 * Z_load_min / (K_ss * K_re); fprintf(I段整定阻抗 : %.2f j%.2f ohm\n, real(Z_set1), imag(Z_set1)); fprintf(II段整定阻抗: %.2f j%.2f ohm\n, real(Z_set2), imag(Z_set2)); fprintf(III段整定阻抗: %.2f j%.2f ohm\n, real(Z_set3), imag(Z_set3));K_rel1取0.85意味着I段保护范围是线路全长的85%末端15%留给了II段避免本线路末端经过渡电阻短路时I段超越误动。K_ss和K_re只影响III段III段的目的是躲负荷而不是保护某条具体线路所以数值上通常比I、II段大得多。算完这三个值Z_set2会直接用在第四章画阻抗圆和R-X轨迹验证里。3. PSCAD4.2 主电路建模单端电源两级线路的元件参数3.1 主电路拓扑与断路器位置主电路采用单端电源两级线路结构三相电源经变压器升压后接线路1线路1末端是断路器B1B1之后是线路2线路2末端接负荷。B1装在线路1与变压器之间对线路1是主保护对线路2是后备保护。单端电源的好处是故障电流方向固定距离继电器不会出现方向判别上的歧义适合先把三段式配合逻辑验证清楚。如果换成双端电源还要考虑功率倒向问题属于后话了。3.2 元件参数与选型表PSCAD4.2元件库里的模型参数按以下表格设置基本对应原模型里其余默认的部分。电源电压和负荷功率是固定的线路阻抗则要和第二章MATLAB算整定值时用到的数值保持一致否则整定与仿真就对不上。元件使用的模型关键参数三相电源Three-Phase Source额定电压230kV频率50Hz单线视图内阻抗默认三相断路器Three-Phase Breaker高压模式单线视图初始闭合由控制信号驱动开断输电线路输电线路模型正序阻抗按Z_AB、Z_BC设定长度对应整定计算假设三相负荷三相负荷模型每相有功100MW无功25Mvar频率50Hz三相故障源Three-Phase Fault内部控制方式中性点接地可配置接地和相间故障电源内阻抗默认即可仿真重点是保护动作逻辑而不是电源的暂态特性。断路器选高压模式是为了让开断电弧的显示更直观实际上PSCAD里这个选项主要影响可视化效果不影响控制逻辑。负荷每相100MW、25Mvar对应的是一个重负荷工况III段整定时要确保最小负荷阻抗不落在动作圆内。3.3 故障源内部控制与故障时序配置故障源采用内部控制方式即由故障源内部的时控逻辑触发不需要外部输入信号。仿真中把故障发生时刻设在0.2s持续时间0.2s也就是说0到0.2s系统正常运行0.2s到0.4s是故障窗口0.4s之后依赖断路器动作切除故障。设置步骤在元件库拖入Three-Phase Fault控制方式选内部控制中性点接地勾选。故障发生时刻填0.2s持续时长填0.2s仿真总时长设为0.5s。做A相接地故障时只闭合A相与地之间的故障开关做AB相间故障时短接A相和B相并且不接地。相间故障设置里最容易踩的坑是中性点接地选项忘关。A相接地要求故障点对地形成通路所以必须接地AB相间故障如果还留着接地就成了AB两相接地短路故障电流的零序分量会明显改变波形特征和论文里对照的相间短路就不是一回事了。4. 控制回路从FFT谐波提取到阻抗圆判据再到断路器跳闸4.1 控制回路的信号流结构主电路负责提供故障电流路径控制回路负责在故障发生后几十毫秒内完成判断。整体信号流是母线电压互感器和电流互感器采集电压电流送入FFT组件提取基频幅值相角再送到线对地阻抗和相间阻抗模块计算R和X阻抗值进入阻抗圆组件得到0/1逻辑逻辑信号经过时间元件和或门综合最终驱动断路器B1。这个链路里任何一个环节的输入输出对不上波形就会表现成故障发生了但断路器不动。4.2 FFT组件与序电流处理FFT组件的作用是把畸变的故障波形拆成基频、谐波和直流分量距离保护只关心基频所以谐波和衰减直流都要在这里滤掉。PSCAD4.2的FFT组件可以同时输出三相电压电流的基频幅值和相角还能输出零序、正序、负序分量。接地距离保护需要零序电流来合成补偿电流这一步就在FFT组件内部完成不需要额外搭建序分量滤波器。接地故障时短路电流里含有明显的衰减直流分量如果不做FFT直接拿瞬时值算阻抗测量阻抗会在故障初期剧烈摆动阻抗圆判断容易误动。FFT之后取基波相量阻抗轨迹就稳定了。相间距离保护用的线电压和线电流同样来自FFT输出只是组合方式不同。4.3 线对地阻抗、相间阻抗与阻抗圆配置线对地阻抗模块接收A相电压的幅值相角、带零序补偿的A相电流幅值相角输出直角坐标下的R和X。相间阻抗模块接收AB线电压和AB线电流同样输出R和X。两个模块并联接入保证接地故障和相间故障都覆盖到。阻抗圆组件按以下参数设置参数取值说明圆心Z_set / 2全阻抗圆特性动作区为经过原点的圆半径|Z_set| / 2按第2章计算出的II段整定值输出逻辑区内为1区外为0下游接时间元件和逻辑门这里用的是全阻抗圆特性整定阻抗作为圆的直径。实际工程里更常见的是姆欧圆偏移圆圆心在Z_set/2处且圆不过原点目的是避免故障点靠近保护安装处时电压过低形成死区。仿真阶段先用全阻抗圆把三段逻辑跑通再改成偏移圆对比动作边界的变化是加深理解的一个好办法。阻抗圆输出的0/1逻辑还要配合时间元件I段直接输出II段加0.3到0.5s延时III段按阶梯时间配合。三个时间元件的输出进同一个或门或门输出就是控制信号TS。TS接到断路器控制端后具体通过一个选择器实现TS为高电平时选择器把常数1送到断路器开断控制端断路器随即跳闸。原文波形里TS和B1完全重合就是这个选择器直通逻辑的结果。4.4 用MATLAB还原R-X阻抗轨迹验证动作边界PSCAD自带的波形图能看电压电流但要看测量阻抗在R-X平面上的运动轨迹还是导出数据用MATLAB画更方便。常见做法是把FFT组件输出的基频幅值和相角通道导出为CSV再用下面的脚本处理% plot_rx_trajectory.m % 读取PSCAD导出的基频幅值/相角数据, 在R-X平面绘制阻抗轨迹 % 列顺序: 时间, U幅值, U相角(度), I幅值, I相角(度) data readmatrix(phase_a_ground.csv); t data(:,1); Um data(:,2); Up data(:,3); % A相电压基波幅值与相角 Im data(:,4); Ip data(:,5); % A相电流基波幅值与相角 U Um .* exp(1j * deg2rad(Up)); % 电压相量 I Im .* exp(1j * deg2rad(Ip)); % 电流相量 Z U ./ I; % 测量阻抗 % 取故障中段稳态窗口, 避开0.2s故障投入和0.4s跳闸暂态 idx (t 0.22) (t 0.38); figure; plot(real(Z(idx)), imag(Z(idx)), LineWidth, 1.5); hold on; % 画II段整定阻抗圆: 圆心Z_set2/2, 半径|Z_set2|/2 Cx real(Z_set2) / 2; Cy imag(Z_set2) / 2; r abs(Z_set2) / 2; theta linspace(0, 2*pi, 200); plot(Cx r*cos(theta), Cy r*sin(theta), k--); xlabel(R / ohm); ylabel(X / ohm); axis equal; grid on;readmatrix读取的列顺序必须和PSCAD导出时的通道顺序一致建议导出前把无关通道删掉只留时间、电压幅值、电压相角、电流幅值、电流相角五列。deg2rad是把PSCAD输出的角度单位转成弧度这一步漏了的话相量角度整体偏移阻抗实部和虚部都会错。取0.22到0.38s窗口是为了避开故障投入和断路器跳闸这两个暂态点暂态窗口里的阻抗点会乱飞不适合判断稳态动作边界。轨迹进入虚线圆内的时刻就是距离保护判出区内故障的时刻。5. 故障波形判读断路器动作时刻与相角突变特征5.1 A相接地故障的三个关键时段A相接地故障波形按时间轴分成三段。0到0.2s是正常运行段三相电压电流对称正弦有功稳定在100MW无功稳定在25Mvar。0.2s故障投入后A相电压明显跌落电流急剧增大三相电压电流开始不对称。紧接着控制系统输出高电平断路器B1动作电流迅速回落到接近0电压恢复三相对称状态。判断保护是否正常动作就看故障发生后到电流归零之间是否有约一个周波的延时这个延时主要来自FFT计算和断路器动作时间是正常的。5.2 AB相间短路与接地故障的波形差异把故障源类型切成AB相间短路、其他参数不动得到的波形总体趋势和单相接地一致但故障期间的波形形态不同。单相接地时非故障相电压会抬升这是零序电压叠加的结果AB相间短路时C相基本维持正常电压故障相间电压降低短路电流的幅值也和接地故障不同原因是接地回路电阻和零序阻抗的参与程度不一样。看波形时先认故障相再认有无零序分量基本就能判断设置的是接地还是相间故障。5.3 用整定边界反向验证模型波形正确不等于整定正确我一般会用两个小实验反向验证。第一个是把故障发生时刻从0.2s改到0.35s断路器动作时序应该整体后移且故障发生后到跳闸的延时不变这验证了控制逻辑是纯故障触发的不依赖仿真起始时刻。第二个是调整故障点在线路上的位置或者修改Z_AB的数值让测量阻抗刚好落在II段阻抗圆边界附近观察断路器是否在I段快动和II段延时动作之间切换。把这个边界找出来和第二章MATLAB算出的Z_set2对照偏差在10%以内说明模型参数前后一致如果偏差很大优先检查线路阻抗单位是欧姆还是标幺值这是PSCAD里最隐蔽的单位坑。本文还有配套的精品资源点击获取