缆-线混合配电网的单相接地故障选线是个特别容易翻车的问题。我最早被它缠上是接手一批10千伏混合线路的录波数据用零序幅值比较和五次谐波两种成熟判据去选结果结论互相矛盾现场人员按其中一个结果排查了半天故障点根本不在这条线。后来把思路转到故障后的短窗暂态数据上用S变换提取时频能量特征才把线路特征真正剥出来整套方法算是在PSCAD里反复仿真验证后落了地。它解决的是小电流接地系统里单相接地后怎么从多条出线中准确挑出故障线路的问题。这篇就把我的思路、建模过程、算法实现和踩坑记录完整写出来给正在做配网故障选线的研究生和工程人员一条能直接复现的路径。1. 问题拆解缆-线混合网络为什么难选线1.1 缆-线混合结构的电气差异配电网越建越密城市核心区电缆化率很高到了城郊结合部又接回了架空线这就形成了典型的缆-线混合馈线。电缆和架空线的单位参数差着数量级三芯交联电缆每相对地电容普遍在0.2到0.3微法每公里架空线的对地电容只有0.008到0.012微法每公里波阻抗上电缆大概二三十欧姆架空线三四百欧姆。这意味着同样一次单相接地暂态分量在电缆段和架空线段里的传播速度、衰减程度、反射折射特性都不一样零序电流波形在母线上看过去是两种线路特征叠加之后的混合体想用一个简单的幅值倍率去标定它就很难。我还记得第一次在PSCAD里搭出这种混合结构的时候把故障设置在缆-线连接点附近录出来的三条非故障出线零序电流居然有一个波形高阶分量高度相似单凭肉眼完全看不出差异。这个现象背后的物理原因就是行波在波阻抗不连续处的折反射故障产生的暂态行波每经过一个电缆-架空线接头就分出一部分能量继续传播其余部分反射回来叠加母线上看到的必然是各线路、各接线段折反射波的叠加。所以选择选线特征量必须对暂态过程有足够细致的刻画能力能把这个叠加结果按频带和时间切分开来看这也是我后来坚持用S变换而不是普通傅里叶变换的直接原因。1.2 常规判据在混合网络中的失效机理配电网中性点不接地或经消弧线圈接地时单相接地故障电流很小稳态量本身就没有明显的幅值优势。经典的零序电流幅值比较法在这种系统里原本就表现一般到了缆-线混合结构里更不可靠——电缆段较长的非故障线路其零序电容电流幅值完全可能超过电缆段短而架空线段长的故障线路。五次谐波法同理谐波电流受消弧线圈补偿度和非线性负载影响混合线路里各段谐波分布又不同实际用起来经常误判。首半波法依赖故障初始相角如果刚好在相角过零附近接地暂态首半波的极性特征会被淹没高阻接地时更明显。我做过一个最简单的对比仿真中性点经消弧线圈接地系统里L2线路电缆段发生200欧姆单相接地四条出线零序电流稳态幅值之比大概是1比0.9比1.1比1.2故障线路L2反而排不进前两名。这说明在混合网络中基于稳态幅值或单一谐波量的方法本质上是把时间窗口拉得太长把最有辨识度的暂态差异平均掉了。后来我把注意力放到故障后几毫秒的短窗暂态信号上情况立刻不一样——这也是整套方案选择短窗数据的出发点。配网里做保护配合时大家还会经常和反时限过流模型打交道但那种稳态延时逻辑在选线问题上是帮不上忙的选线必须走暂态信息这条路。2. 核心方案S变换与短窗能量判据的设计逻辑2.1 S变换的时频分析优势要用短窗数据捕捉暂态能量分布手里的工具得有足够的时频分辨率。短时傅里叶变换的窗长固定频率高了时间分辨率必然变差小波变换能自适应但选基小波、分解层数和频带划分都影响结果物理量纲也不够直观。S变换可以看作带相位信息的连续小波变换它的窗函数是高斯窗窗宽随频率自动调整高频部分窗窄、时间定位准低频部分窗宽、频率分辨好时频矩阵直接反映信号在各频率和各时刻的幅值和相位。这个特性跟故障暂态信号很搭故障初瞬间的高频分量本身持续时间极短需要高时间分辨率才能定位而后面的低频衰减分量频率低需要足够的频率分辨率才能区分。公式层面S变换的定义是信号先乘一个随频率收缩的高斯窗再做傅里叶变换离散实现时等价于把信号频谱平移后加高斯窗再反变换。它保留了信号的绝对相位信息而且不像Wigner-Ville分布那样有交叉项干扰对非平稳故障暂态信号来说非常合适。做配网故障选线的文献里用S变换提取零序电流特征的成功案例不少这里就不再重复推导直接讲工程上怎么用更实际。2.2 短窗数据的选择逻辑所谓短窗指的是取故障发生后很短一段时间的采样数据一般取5毫秒到20毫秒也就是工频周期的四分之一到一个周期。这个窗口选择有讲究故障暂态能量集中在前几个毫秒取短窗可以把稳态工频分量和后续衰减振荡剔除掉大部分但窗口太短又会丢失低频衰减成分的能量同时影响S变换的频率分辨率。从我的仿真经验看10毫秒窗口基本够用既拿到2到3个完整的振荡周期数据用于低频段分离又不会让工频稳态分量占据能量主导。采样频率方面故障暂态的有用频率成分一般到3到6千赫兹为了把高频段信息保留下来PSCAD里的仿真步长我习惯设成10微秒相当于采样率100千赫兹。100千赫兹意味着S变换能做有效分析的频率上限到50千赫兹实际故障暂态频率根本用不到那么高取这个步长主要是为了保证波形细节不丢。数据窗长与频率分辨率是一对矛盾窗口短则每个频率点的带宽粗做S变换时相邻频带可能粘连。折中办法是加采样点数或做尾部补零但补零会增加计算量我在后面实现部分会给出一个相对均衡的配置。2.3 能量判据的构造与增强设计拿到各条出线的零序电流短窗数据之后做S变换得到时频复矩阵取模平方就是信号能量在时频平面上的分布。构造判据的第一步是确定特征频带把故障后10毫秒内的零序电流做频谱分析挑出能量集中的连续频段常见的是300到3000赫兹这个区间具体边界根据仿真和实测数据微调。第二步是计算每条线路在该频带内的能量用一个矩阵切片求和就能得到。第三步是算相对指标用最大能量与次大能量之比作为选线置信度同时配合能量占比判断防止某条非故障线路因为离故障点近而能量接近故障线路。针对缆-线混合结构我额外加了一个增强措施把特征频带细分为高频段和低频段分别计算能量权重后再合成。原因是电缆段对高频分量的传输衰减小、耦合强架空线段对低频分量占比贡献大如果只用一个宽频带求和某些接线方式下非故障线路的宽频能量会把故障线路顶下去。多频带加权的好处是让两个频段的能量互为佐证故障线路总是至少在一个频段上明显占优。这个现象在大量仿真组合里反复出现后来的判据就按这个思路固化了。3. PSCAD仿真建模与实操步骤3.1 搭建10kV缆-线混合配电网模型PSCAD里做这套仿真关键是把线路模型选对。架空线用贝瑞隆模型或者频率相关模型都行但电缆一定要用频率相关相域模型因为电缆的频变特性对暂态波形的准确性影响非常大。我在Master Library里取的是Frequency Dependent (Phase) Model一个三相电缆段、一个架空线段中间用母线节点连接模拟缆-线接头。电缆的几何参数按10千伏三芯交联聚乙烯电缆的典型值填写我用到的核心参数整理成了下面这张表。参数项电缆段取值架空线段取值导体类型铜芯外半径4.5毫米钢芯铝绞线LGJ-120绝缘/半径绝缘厚5.5毫米介电常数2.3导线外径15.2毫米护套/弧垂护套外半径10毫米埋深0.8米电路最低点高度10米土壤电阻率100欧姆米100欧姆米典型单位电容约0.25微法每公里约0.01微法每公里典型波阻抗约30欧姆约350欧姆系统侧用一个110千伏等值电源接主变压器容量20兆伏安接线组别YNd11变压器低压侧10.5千伏中性点引出接消弧线圈。消弧线圈电感要按系统总对地电容电流来配我这里四条馈线加母线设备的总对地电容折算后约为每相1.1微法按过补偿5%计算电感取3.0亨左右。母线侧挂四条馈线馈线1是重点研究的缆-线混合线架空段4公里接电缆段2公里馈线2是纯架空线6公里馈线3是纯电缆线2公里馈线4是另一条混合线架空3公里加电缆1.5公里。四条线有差异选线算法的检验才有意义。搭模型时如果PSCAD版本较新务必确认Master Library路径挂载正确不然Cable和Fault这类元件经常调不出来。3.2 故障模块与测量配置故障元件直接选PSCAD的三相故障模块故障类型设成单相接地接地电阻我设成可切换的5欧姆、50欧姆、200欧姆、500欧姆、1000欧姆都跑一遍。故障起始时间设0.2秒持续0.1秒故障位置分三组电缆段中点、架空线段中点、缆-线连接点临近点。为了模拟工程上最常遇到的边界条件我还把故障设置在母线和馈线1的连接处、以及故障线路末端每个位置都配合不同过渡电阻组合成工况矩阵一共跑了三十多组仿真。零序电流的测量不能直接在PSCAD图形界面上读三相电流我是在每条馈线出口加一个电流测量模块输出三相瞬时电流再用加法器把三相相加得到三倍零序电流。这个信号就是后面S变换的输入。需要说明的是PSCAD里的加法器延迟和采样步长会引入微小误差步长10微秒完全可以忽略但别把数据记录间隔设得比仿真步长还大否则高频细节会被抽掉。故障模块里设置接地电阻时PSCAD对纯电阻和带弧光电阻的表现也不同想模拟真实电弧接地可以把电阻换成随时间变化的元件但作为算法验证恒定电阻已经足够。3.3 仿真参数与工况矩阵设计仿真步长设10微秒时长0.5秒故障前后预留0.1秒的稳态数据。PSCAD的数据记录我习惯同时开两类一类是运行时Plot窗口里的快速观察通道看波形形态和趋势另一类是输出到.out文件的数据通道供后续MATLAB离线分析。不同PSCAD版本读.out文件的方式略有差异新版本可以用Python接口直接读旧版本要么用软件自带的转换工具要么手写二进制解析我在实施的时候用MATLAB写了个读取函数一次性把所有通道读进工作区省得每次手工导出。工况矩阵的设计建议按故障周期来做固定系统参数和故障类型先扫故障位置再扫过渡电阻最后扫故障初相角。初相角这个变量很容易被忽略但它对暂态首半波能量影响很大我单独跑了一个0度、45度、90度初相角的对比组结果在判据稳定性分析里很有价值。每个工况跑完把四条馈线的零序电流通道数据导出按故障时刻对齐截窗保存成统一的mat文件后续算法处理全部自动化。这里提醒一句导出通道命名要规范我一开始命名成I0_1、I0_2这种好在还有注释文档否则几十个文件混在一起根本分不清哪条是馈线1、哪条是馈线3。4. 选线判据实现与仿真结果验证4.1 S变换算法的MATLAB实现PSCAD负责出数据S变换和能量计算在MATLAB里做。S变换的离散实现并不复杂核心代码不到二十行我在这里给出一个教学版本function ST s_transform_energy(sig, fs) % sig输入零序电流序列列向量 % fs采样频率 N length(sig); H fft(sig); halfN floor(N/2); ST zeros(N, halfN - 1); for m 2:halfN % 频域平移索引 km mod((0:N-1) (m-1), N) 1; % 高斯窗乘以平移后的频谱再反变换 G H(km) .* exp(-2*pi^2 * (0:N-1).^2 / (m-1)^2); ST(:, m-1) ifft(G); end f_axis (2:halfN) * fs / N; end这段代码去掉了直流项m从2开始避免了分母为零的问题。实际使用中我会先把零序电流减掉其均值再做一次带通滤波例如二阶巴特沃斯带通范围200到5000赫兹这样低频漂移和高频噪声都不会污染时频矩阵。S变换的结果是复数矩阵每一行对应一个时间平移每一列对应一个频率点取模平方后就是该信号在时频平面上的能量分布。要注意的是这个循环版实现速度偏慢数据量大时建议用向量化或直接调用开源工具包。4.2 特征频带能量计算与选线判据能量计算按前面说的流程来先截取故障后10毫秒的数据调用S变换函数得到时频矩阵再按特征频带范围做频率索引切片把该频带内所有时间和频率点上的能量求和。多频带增强方式就是把频带切成300到1500赫兹与1500到3000赫兹两段分别求能量后加权相加我用的权重系数是低频0.6、高频0.4这个比例是在大量仿真组合里收敛出来的不同系统可以微调。选线判据具体操作分三步。第一步计算四条馈线的综合能量值第二步找出最大能量值和次大能量值求比值第三步把最大能量线路判为故障线路同时输出置信度指标。我设的置信度判据是最大次大能量比大于1.8判定有效低于这个值就报“疑似故障但需进一步分析”。这个阈值不能定太死因为高阻接地故障时暂态能量整体衰减比值会收窄定太大容易漏判定太小又容易受到非故障线路接近故障点时的干扰。实际工程里这个阈值还可以跟消弧线圈补偿度校验结果联动我现在的版本是固定阈值加多频带权重的组合。4.3 典型工况仿真结果解读下面这组数据取自我在不同工况下跑出来的典型结果能量值做了归一化处理这里以标幺值给出。表格里的最大次大比是选线置信度的核心判据。工况故障位置过渡电阻(Ω)L1能量L2能量L3能量L4能量最大/次大比1L2电缆段中点50.316.520.450.2214.52L2电缆段中点5000.291.830.360.195.13L2缆-线连接点边10000.330.910.400.262.34L3架空线段中点2000.380.271.740.314.6工况1和工况2分别代表低阻和高阻情况故障线路L2的能量优势非常明显低阻时最大次大能量比达到14.5几乎不会误判。工况3是缆-线连接点附近加1000欧高阻接地这种工况最考验判据比值降到2.3虽然仍然能正确选到L2但裕度已经压到警戒线。工况4换到L3故障架空线段的暂态能量特征和电缆段不一样L3能量依然是其他线路的4倍以上。这个趋势告诉我缆-线混合网络里故障线路在宽频带综合能量上始终占据主导但高阻和弱点附近的极端工况必须靠多频带加权和阈值配合才能守住边界。整个仿真过程中三十多组工况全部选线正确说明采用短窗数据加S变换能量的路线在缆-线混合网络里是站得住脚的。5. 常见问题排查与仿真避坑技巧5.1 模型搭建相关坑点PSCAD里搭混合配电网最容易踩的坑是电缆模型选错。用Bergeron模型模拟电缆频率相关特性全丢仿真出来的暂态波形在几百赫兹以上衰减速度明显不对直接导致后续S变换的能量分布和集中式线路模型算出来的对不上。另一个坑是电缆和架空线连接处的节点处理PSCAD里如果直接拿普通导线连接波阻抗不连续处的折反射条件会失真应该在连接处用内部节点或适当长度的母线段连接。变压器中性点引出消弧线圈时也别忘了在消弧线圈上并联一个阻尼电阻否则补偿度附近的谐振效应对零序暂态影响很大那一部分能量会被完全放大非故障线路的能量甚至会出现虚假峰值。此外PSCAD里的电源模型用理想电压源加内阻抗作为系统等值时内阻抗参数一定要按短路容量折算不能随便填。我第一个版本电源内阻填的默认值结果零序回路阻抗严重偏移仿真出来的单相接地电流比实际系统大出一倍多选线结果倒是没翻车但能量值和标幺化后的分布比例已经很不可信。这些问题排查下来其实都在模型合理性上仿真工具的故障往往不是软件报错而是物理过程失真。5.2 S变换实现相关坑点MATLAB实现S变换时最隐蔽的问题是循环里频域索引写错。我曾经把频域平移的取模索引写错一个位置结果全频带能量都偏离真实值选线结果时对时错排查了整整两天。建议先在简单正弦信号上验证S变换正确性再接入真实数据。另一个常见问题是直流项。零序电流在故障前后往往有直流偏移如果不做去均值或高通滤波直流分量会在时频矩阵的低频列里形成很宽的能量带把特征频带的能量算到失真。我的做法是进S变换之前先减均值再做200赫兹高通滤波效果稳定。数据窗截取也要注意边界效应。S变换在时间窗两端存在一定的边缘失真截的窗越短越明显。解决方式是把数据窗前后各多取一段数据做S变换之后只保留中间感兴趣的区间用于能量积分边缘效应就能隔离在统计范围之外。这个细节不处理短窗数据的分辨率优势就体现不出来。还有一个容易被忽略的点是零序电流三条相线相加后的比例系数有的仿真记录里记的是3I0有的记的是I0做能量计算前必须确认单位一致否则四条馈线之间差三倍判据直接就乱了。提示S变换计算前务必先对零序电流做去均值和带通滤波这是能量判据稳定性的基础。跳过了这一步后面所有结果都不可信。5.3 选线判据边界条件验证心得整套判据在普通工况下表现很稳真正花时间的都在极端边界条件上。我总结了几条比较实用的经验。第一过渡电阻从低到高扫的时候能量值下降是非线性的从几十欧到几百欧变化最剧烈1000欧以上反而趋于平缓所以不能单纯按比例缩放阈值。第二故障初相角接近0度时首半波能量很低这时候能量判据依赖的是后续衰减振荡成分S变换的短窗优势仍然能保持辨识度但初相角过零附近的信号一定不能截得太短至少保留10到15毫秒。第三缆-线连接点附近故障时折反射波会让故障线路与非故障线路的波形时域上高度相似但频带能量分布仍然有差异多频带加权判据比单频带判据的裕度要高出一截这个在选线装置的工程实现中非常关键。做这些边界工况仿真时我还养成了一个习惯每组参数固定跑三次小扰动比如把消弧线圈抽头电压调高百分之二、或者把故障初相角微调几度观察选线结论会不会跳变。这个方法能快速暴露判据的脆弱面。整套判据后续要落地到实际装置还需要跟现有的反时限过流保护模型做配合用选线结果驱动出口逻辑而不是单纯停留在仿真分析层。5.4 高频问题与排查速查表这段时间反复仿真下来遇到的高频问题基本集中在几个固定环节整理成速查表放在下面方便有需要的朋友直接对照。现象可能原因处理办法零序电流波形毛刺多仿真步长过大或记录间隔过大步长降到10微秒记录间隔与步长一致S变换能量集中在极低频率未去直流或工频分量未滤除减均值加200赫兹高通滤波高阻接地时判据翻车特征频带选得太窄放宽到300至3000赫兹并启用多频带加权非故障线路能量接近故障线路消弧线圈谐振效应或模型缺阻尼电阻消弧线圈并联阻尼电阻检查补偿度设置缆-线接头附近误判短时间内窗截断导致边缘失真前后扩窗再截取只保留中间有效区间不同PSCAD版本读不出.out通道格式或字节序不同确认版本用Python接口或自写解析读取判据阈值不稳定系统对地电容和消弧线圈参数未归算按实际系统短路容量和电容电流重新整定仿真做到最后我最大的体会是故障选线这类问题算出来只是第一步真正值钱的是在各种边界工况下都能稳定选对线。S变换在缆-线混合网络里最大的价值是把暂态信息按频带拆开让故障特征从电缆和架空线结构差异造成的信号混叠里暴露出来。我个人现在做新课题拿到的第一批数据一定是先过一遍这个流程确认频谱和时频能量形态再谈选线。下一步我打算把这套判据移植到硬件在环平台做实时验证有条件的话再跟现场录波数据做一轮盲测。先把这条技术路线放出来给同样在这个方向上努力的同行当块铺路石。