在配电网单相接地故障选线这件事上缆-线混合网络是我做过最“磨人”的一类对象。电缆段和架空线段交错之后整条馈线的波阻抗、对地电容和暂态衰减特性都变得不一致故障线路的零序电压电流特征会被不同段落的参数“带偏”传统基于稳态工频量或简单谐波幅值比较的选线逻辑经常在仿真阶段就被打回原形。这篇内容我把一套基于短窗数据S变换能量的选线方法从问题拆解、特征提取、判据设计到PSCAD仿真建模和工况验证完整梳理一遍重点讲清楚每一步为什么这么做以及实操中哪些环节最容易坑人。适合正在做小电流接地系统选线研究、或者被缆-线混合配电网仿真结果折磨的工程师和研究生参考。1. 缆-线混合馈线单相接地传统选线方法为什么频频翻车1.1 单相接地后的零序暂态才是选线的“富矿”中性点不接地或经消弧线圈接地的配电系统单相接地故障最麻烦的点在于故障电流小。稳态情况下故障点流过的电流由系统对地电容决定往往只有几安到几十安比负荷电流小一个数量级。靠稳态基波电流幅值做选线在现场基本没有裕度更何况消弧线圈投入后稳态零序电流还要再打一个折扣。但故障发生的初始瞬间完全是另一回事。母线电压在故障相上突然塌陷所有健全相通过对地电容向故障点充电故障线路对地电容则通过故障点放电这个暂态过程会产生幅值远大于稳态电流的高频分量有时候能达到稳态电流的十几倍甚至几十倍。更重要的是故障线路与健全线路在这个暂态里的零序电流极性相反幅值差异明显。所以现在主流的选线思路几乎都围绕暂态量做文章关键是怎么把这段转瞬即逝的信号稳定地换算成一个能比较的特征量。1.2 缆-线混合网络的参数离散让信号变得“非典型”架空线和电缆在电气参数上差别非常大。电缆的单位对地电容通常是架空线的十几倍到几十倍波阻抗却只有架空线的五分之一到六分之一。一条馈线里既有电缆又有架空线相当于一条链路被分成了好几段电气特性完全不同的“零件”故障暂态在这条链路里传播时会产生折反射零序电流的波形不仅取决于故障点在哪一段还取决于该段前后相邻线路的参数匹配情况。我做过一个对比仿真同样在距离母线5公里的位置发生单相接地全架空线路的零序暂态主频分量能到几千赫兹衰减也快而如果是电缆段故障信号频率整体被拉低波形明显更“钝”。如果用固定频点的幅值比较法比如只看5次谐波或某个特征频带故障点所在段一变这个频点上的幅值排序就可能跟着变选线结果当然不稳定。这也是为什么简单傅里叶方法在混合网络上表现不佳的根本原因。1.3 消弧线圈补偿让稳态方向判断更难如果系统中性点还接了消弧线圈稳态问题会被进一步放大。消弧线圈按过补偿方式运行稳态零序电流不再是容性电流占主导而是补偿后的残余电流本身数值很小相位关系也被改变。很多基于稳态分量或功率方向原理的选线装置在这种系统里基本失去方向分辨能力。不过消弧线圈有一个特点它的电感电流建立需要时间。故障后几个毫秒到十几毫秒内电感支路还来不及充分参与暂态过程故障初期的零序暂态特征反而还保留着原始的电容放电形态。这意味着只要数据窗取得足够短、足够靠前就能绕开消弧线圈对稳态分量的“污染”拿到可靠的故障特征。短窗数据方案的核心优势就在这里。2. 从零序暂态里提取可判特征短窗数据与S变换2.1 窗长取多少这是一个先有靶子再射箭的问题短窗数据到底多短我试过从半个周波到几个周波不等的窗口长度最后定在故障后40ms左右也就是工频两个周波。理由有三个第一消弧线圈的电感电流通常在故障后10到20ms内开始明显建立40ms窗口里后半段虽然已经受一点影响但前半段暂态的主导特征还在。如果窗长拉到100ms以上稳态补偿量就会严重稀释暂态特征。第二S变换在低频段需要足够的时域长度才能分辨窗太短了低频段能量算不准不利于处理高阻接地这种能量集中在低频的故障。第三采样率20kHz时40ms就是800个点做离散S变换的计算量完全可以接受。实际操作时我习惯在故障前再留3到5ms的“预触发数据”一方面用来和故障标志位对齐另一方面可以用故障前信号做差分或者基线校正。这样最终截出来的数据窗大概是故障前5ms到故障后40ms共45ms。2.2 S变换和STFT、小波到底有什么区别S变换可以理解成短时傅里叶变换STFT的改良版也可以理解成小波变换的一种特殊形式。STFT用固定宽度的高斯窗窗宽定了频率分辨率就定了高频的时间分辨率和低频的频率分辨率没法兼顾。小波变换能自适应调整时频分辨率但小波基函数的选择很依赖经验而且一般小波变换的相位信息不直观。S变换的核心是让高斯窗的宽度随频率反向变化频率低的时候窗宽能看到更窄的低频频带频率高的时候窗窄能精确定位突变时刻。公式上可以写成S(τ, f) ∫ h(t) · (|f| / √(2π)) · exp(-f²(τ − t)² / 2) · exp(−j2πft) dt从使用者的角度看不需要死磕公式只需要把它当成一个能输出二维时频矩阵的工具。矩阵的行对应频率列对应时间每个元素是一个复数模值就代表该时刻该频率分量的强度。和固定滤波器组相比S变换不用费劲设计带通滤波器组改频带范围只需在结果矩阵里切行这对前期研究和参数扫描特别友好。2.3 为什么用“能量”而不是幅值或相位非平稳故障信号里单看某一时刻的瞬时幅值和相位都很不稳定。故障初相角不同第一个波峰的极性甚至幅值都会变而相位信息在故障点和母线之间测量时还受线路反射波干扰严格来说每条线路测到的相位也不是理想的反相180度关系。能量是对整个短窗内所有时刻、所有频点的模值平方做累积相当于把整个暂态过程“积分”了一遍对局部噪声、相位抖动、采样起始点微小偏移的鲁棒性都更好。在S变换结果矩阵上就是取模后逐点平方再在时间和频带上求和。这样做还有一个好处故障线路因为故障点直接泄放对地电容储能它的暂态能量天然远大于健全线路而健全线路之间能量差异则主要由线路长度和电容决定规律性更强方便后续归一化处理。3. 选线判据搭建短窗能量怎么算、怎么比较3.1 从零序电流到特征能量的完整计算流程假设系统有N条馈线先采集每条馈线的零序电流 i_k(t)k 1,2,...,N。对每条电流按前面说的短窗截取做S变换得到复时频矩阵 M_k[t][f]。特征能量的计算分三步对矩阵每个元素取模并平方得到时频能量谱 P_k[t][f] |M_k[t][f]|²。在选定的频率范围内累加通常取下限 f_min 50Hz上限 f_max 2kHz得到随时间变化的能量轨迹 E_k(t)。再对短窗内的时间轴累加得到每条馈线在这个数据窗内的总特征能量 E_k。最后得到的 E_k 就是选线判据的基础数据。在编程实现时S变换矩阵可以按列优先存储能量累加用矩阵切片操作速度非常快。3.2 线路长度和电容差异一定要归一化这里有一个新手很容易忽略的坑直接比较 E_k 往往得到的是“长线路优先”。因为长线路对地电容大健全时自身充放电能量本来也大一条20km的健全线路能量可能比一条5km的故障线路还高。我一般用单位长度平均能量来处理也就是 E_k E_k / L_k其中 L_k 是该馈线的总长度。如果知道各条线路的对地电容 C_k用 E_k E_k / C_k 更准确。混合网络中电缆段占比不同的两条线路总长一样但电容差几倍这时候按总长度归一化仍然不够建议优先按对地电容归一化。归一化之后故障线路只要不是极端高阻接地它的 E_k 通常还是会明显高于其他线路。如果两条线路的 E_k 太接近下一步可以用极性校验兜底。3.3 相对能量系数与极性校验为了把数值变成一个直观、可设阈值的判据我习惯再定义一个相对能量系数p_k E_k / Σ_{i1..N} E_ip_k 表示第 k 条线路在全部馈线总特征能量里占的比例。故障线路的 p_k 往往超过0.5健全线路通常在0.1到0.3之间。在输出结果时直接把每条线路的 p_k 列出来比看原始能量值直观很多。极性校验的做法是在S变换矩阵中选取主频带通常是能量最大的那几个频率点看短窗内各线路S变换复数实部的初始方向。故障线路的初始暂态极性与健全线路相反健全线路之间基本一致。如果母线故障则所有线路的极性都一致且 p_k 分布相对平均没有特别突出的最大值。极性校验不需要很高精度只需要判断正负方向所以受高阻接地导致的幅值衰减影响较小。3.4 初始阈值怎么设在没有做扫描之前可以先按下面的经验值启动后面再用仿真扫描结果来标定存在一条线路满足 p_max 0.5且它的极性与至少两条健全线路相反判这条线路故障。所有线路 p_k 都小于0.4且极性一致判母线故障。如果最大 p_k 在0.4到0.5之间需要结合极性校验和原始波形二次确认不直接给出选线结果宁可不选也不要误选。这个初值不保证在任何系统都适用但作为起点去扫描各种故障条件能节省大量调参时间。4. PSCAD里搭出可用的缆-线混合网络仿真模型4.1 主接线结构与中性点接地方式模型不需要太复杂但必须覆盖缆-线混合的典型特征。我搭的是一套10kV配电网仿真算例主变压器二次侧星形中性点通过开关连接消弧线圈。消弧线圈的补偿度按10%过补偿设置这是国内中性点经消弧线圈接地系统最常见的运行方式。馈线设计上建议至少4回包含三类纯电缆馈线、纯架空馈线、缆-线混合馈线。混合馈线如果条件允许可以做成“架空电缆架空”和“电缆架空电缆”两种结构。这样在验证判据时能同时检验故障点位于电缆段和架空段两种情况下的表现。每回馈线长度不要太悬殊我用了6km到14km这个范围算出总对地电容后反推消弧线圈电感值。4.2 电缆和架空线的模型选择别用集中参数糊弄PSCAD里线路模型有杆塔模型、贝杰龙模型、频率依赖模态模型、频率依赖相域模型。做单相接地选线必须考虑零序暂态过程电缆和架空线我都建议用频率相关相域模型。贝杰龙模型对参数随频率变化的模拟很弱高阻接地时误差会明显影响高频分量的衰减规律而选线判据恰恰依赖暂态高频特性。电缆模型参数可以参考10kV交联聚乙烯电缆的典型结构导体截面120平方毫米级别内半导电层、绝缘层、金属屏蔽层、外护套分层设置埋深和土壤电阻率也要填。架空线用典型10kV配电线路参数导线半径8mm左右相间距离、导线弧垂、地线参数按常规杆塔布置填入。参数不需要极度精确但数量级不能错否则频率依赖关系会失真。4.3 故障模块和短窗触发条件设置PSCAD的故障元件直接选单相接地故障故障相、接地电阻、故障起始时刻都是可以自由设置的。故障持续时间一般设置120ms左右足够覆盖整个短窗。需要注意的是PSCAD的故障元件只在仿真运行开始时启动一次如果要做多工况扫描可以配合Multi-Run或者用外部脚本循环改写故障参数。采样步长设置成20微秒通道步长也保持一致对应采样频率50kHz实际存储时再降采样到20kHz。这个设置比10kHz采样更稳因为PSCAD里高频暂态如果有5kHz以上的分量10kHz采样会混叠。另外一定要输出一个“故障启动”标志信号可以用逻辑表达式在故障开始时置1后面在数据处理时用它精确定位故障发生的仿真波形点。4.4 数据分析放在PSCAD外做S变换不是PSCAD内置的模块硬要在PSCAD里用Fortran自定义组件实现S变换开发成本很高而且调试麻烦。我采用的方式是把仿真结果导出再用外部数值分析软件做S变换和能量计算。PSCAD可以把波形数据保存成.out文件或通过导出功能生成CSV。读取后要做的第一件事不是直接截窗而是检查时间轴是否等间隔。PSCAD在某些输出配置下会按仿真步长输出但记录的数据点数可能和期望不符存在微小抖动需要用线性插值统一到20kHz的等间隔时间轴然后根据故障标志信号把故障前5ms到故障后40ms的数据窗口截出来。5. 故障工况扫描与阈值、误判边界5.1 典型算例结果一条电缆段故障的判据表现以4回馈线的仿真算例为例第1回是混合馈线其中电缆段4km加架空段6km故障点设置在电缆段中点金属性接地。经过短窗S变换能量计算和长度归一化后各线路的相对能量系数如下表馈线编号线路结构长度(km)归一化特征能量相对能量系数 p_kL1电缆4架空6100.830.74L2全架空120.090.08L3全电缆80.120.11L4架空6电缆280.080.07故障线路L1的p_k为0.74远超其他线路。注意这里如果不做长度归一化L2这条12km全架空线路的原始能量可能会比L1更高造成漏判。这个算例也再次说明归一化步骤的必要性。5.2 不同故障条件扫描结果我系统扫了故障电阻0Ω、10Ω、100Ω、500Ω、1000Ω、故障初相角0度、30度、90度、故障位置电缆段、架空段、电缆与架空连接头处等组合。整体规律如下工况特征最突出线路p_k范围第二突出线路p_k范围选线可靠程度金属接地、90度故障角0.65~0.850.06~0.15可靠100Ω、任意故障角0.55~0.700.08~0.18可靠500Ω、故障角接近0度0.42~0.520.18~0.25需极性校验兜底1000Ω、故障角接近0度0.30~0.450.20~0.30存在误选风险母线单相接地0.15~0.400.15~0.35按母线故障逻辑处理可以看到故障电阻增大和故障初相角接近零都会让故障暂态能量变弱故障线路和健全线路的能量差距缩小。这也是这种短窗能量类方法的固有边界不能指望在超高频、超低信噪比工况下还能保持绝对优势。5.3 针对自己系统的阈值标定阈值不能照抄论文也不能照抄我这里的数字。每个网络拓扑、线路长度、电缆占比和补偿度不同p_k的分布范围会整体平移。我的标定方法是把扫描结果中故障线路的最小p_k记为 p_min_fault把健全线路的最大p_k记为 p_max_healthy阈值设在这两个值的中值threshold (p_min_fault p_max_healthy) / 2如果 threshold 低于0.45说明判据裕度不足可以考虑收窄S变换频带到主能量频带或者把归一化方式从按长度改为按对地电容。如果 threshold 高于0.6说明这个网络条件下选线非常容易可以在一定程度上放宽裕度来提高灵敏性。5.4 母线故障怎么和线路故障分开母线发生单相接地时所有出线都流过从母线向故障点放电的暂态电流它们的极性一致能量分布主要取决于各线路对地电容。经过归一化后各线路p_k不会出现一个特别突出的值通常都在0.2到0.4之间。判据上我加了两个条件一是最大p_k和最小p_k的差值不超过0.3二是所有线路初始极性一致。同时满足就判母线故障。实际仿真中母线高阻接地时各线路能量差异会更小反而更容易区分反而是在母线附近线路侧发生故障时由于故障距离短、暂态能量大可能被误判为母线故障这种情况要靠波形起始突变程度辅助判断。6. 做这套仿真实战最容易踩的坑6.1 数据不同步S变换最怕起始点错位这是仿真里最容易出问题的一步。PSCAD导出数据后如果直接按“仿真时间大于故障起始时间”作为截窗条件会因为数字信号处理和故障元件内部逻辑延时导致实际截取的窗口偏移几个采样点。S变换对窗口起点非常敏感几个点的偏移就可能让能量计算结果出现明显波动。我的做法是单独加一个故障标志信号通道故障启动瞬间从0跳到1在数据后处理时用上升沿精确定位而不是用仿真时间判断。一次仿真结束后先检查标志位上升沿对应的数据点索引再结合采样率反推实际故障时刻确保窗口起点误差不超过一个采样周期。6.2 通道步长和重采样问题PSCAD里仿真步长和输出通道步长如果设置不一致导出的波形时间轴并不是严格等间隔的。直接用这种数据做S变换离散傅里叶部分会出现频谱泄漏高频段能量分布会被污染。建议先在PSCAD里把通道步长设成与仿真步长一致比如都是20微秒再在导出后统一用重采样函数插值到20kHz等间隔时间轴。如果嫌数据量大也可以等间隔降采样到20kHz但降采样之前必须先做抗混叠滤波。我吃过亏之后就再也没跳过这步。6.3 高频反射分量被误当成故障特征缆-线混合网络里暂态行波会在电缆段和架空段的连接点反复折反射。S变换能量谱上会出现多个能量峰值时间上对应反射波到达母线的时刻。如果分析频带上限开得过高比如超过5kHz这些反射分量会显著抬高某些健全线路的“能量”造成误判。我根据故障距离和电缆波速估算反射时间粗略算下来10km距离上的主要反射频率在7kHz级别甚至更高。所以把S变换的分析频带上限限制在2kHz附近既能保留主要暂态能量又能避开大部分一两次反射峰。对更长距离的线路这个上限还可以进一步放宽。6.4 高阻接地能量太弱时的对策500Ω以上高阻接地是短窗暂态能量法的软肋故障线路能量会变得很弱。我试过的有效对策有三个一是把数据窗适当加长到60ms左右多积累一些低频暂态能量二是把S变换的分析频带下限从100Hz降到50Hz让工频附近的暂态能量也能进入统计范围三是改用对地电容归一化之后再做相对能量系数通常能把故障线路p_k往上抬几个百分点。不过说句实在话到1000Ω以上结合小故障角时任何暂态量方法都很难保证绝对可靠工程上还需要其他判据冗余。6.5 阈值标定要覆盖“最不利组合”不要只看典型工况就定阈值。我在仿真里发现单看100Ω接地故障阈值设0.4很安全但一旦故障初相角变成15度健全线路的最大p_k能到0.23故障线路的最小p_k降到0.48阈值虽然还能放中间但裕度明显变小了。如果继续把故障角压到接近0度同时故障电阻升到500Ωp_max_healthy和p_min_fault几乎要贴在一起。所以最终标定阈值时一定要扫到故障电阻、故障角、故障位置三个维度交叉的最不利组合而不是取每个维度的典型值。这一步多花的时间价值很大。6.6 消弧线圈补偿度不能只设一种既然标题里是缆-线混合配电网络消弧线圈的作用不能忽略。过补偿10%是常见设置但我建议把补偿度从5%到20%之间各跑一遍。补偿度不同时暂态低频分量的衰减速度和稳态残余电流方向都会有变化从而影响短窗后面半段的能量分布。如果时间紧张至少跑三种欠补偿5%、过补偿10%、过补偿20%。我遇到过一种情况过补偿10%时判据很干净欠补偿时有一回健全线路的p_k达到了0.35虽然没有越过阈值但和第二名之间的差距非常小说明系统对补偿度是敏感的。这套方法从原理到仿真验证真正需要花心思的地方不在S变换本身而在于把故障特征、数据窗口和网络参数这三者之间的关系摸清楚。我自己在调试过程中最大的体会是短窗数据给S变换提供了一块“干净”的信号切片S变换再给出一个同时包含时域和频域信息的能量地图而判据本身反而可以设计得很简单。每次遇到选线结果不稳定我都会先回到PSCAD模型里去检查线路参数和故障设置是否合理而不是急着调阈值——因为仿真模型失真的话再精细的算法也算不出正确结论。
