基于S变换的配电网单相接地故障精准定位与诊断研究

基于S变换的配电网单相接地故障精准定位与诊断研究一、引言1.1研究背景与意义随着城市化进程的加速和经济的飞速发展，配电网作为电力系统与用户连接的关键环节，其规模持续扩张，结构也愈发复杂。同时，分布式电源、电动汽车充电桩、储能装置等新型设备的广泛接入，使配电网从传统的无源单向辐射网络逐步向有源双向交互系统转变，这在提升能源利用效率和促进可再生能源消纳的同时，也给配电网的安全稳定运行带来了新的挑战。单相接地故障是配电网中最为常见的故障类型，在中压配电网中，尤其是采用小电流接地方式（如中性点不接地、经消弧线圈接地或经高电阻接地）的系统中，单相接地故障发生概率较高。据相关统计资料显示，在我国中压配电网故障中，单相接地故障占比可达70%-80%。当发生单相接地故障时，虽然故障电流相对较小，系统线电压仍能保持对称，在一定时间内可维持对用户的供电，一般允许继续运行1-2小时，但这并不意味着可以忽视此类故障。长时间的单相接地运行会对变电设备、配电设备、区域电网以及人畜安全等造成严重危害。对变电设备而言，故障发生后，变电站母线上的电压互感器可能检测到零序电流，在开口三角形上产生零序电压，导致电压互感器铁芯饱和，励磁电流急剧增加，若长时间运行，极易烧毁电压互感器。例如，在某地区配电网中，曾因单相接地故障持续时间过长，造成多台电压互感器损坏，进而引发大面积停电事故。此外，单相接地故障还可能产生谐振过电压，其幅值可达正常电压的数倍，这对变电设备的绝缘构成极大威胁，严重时会使设备绝缘击穿，引发更严重的事故。对于配电设备，单相接地故障可能引发间歇性弧光接地，产生的谐振过电压会进一步使线路上的绝缘子击穿，造成严重的短路事故，甚至可能烧毁部分配电变压器，使线路上的避雷器、熔断器绝缘击穿、烧毁，还可能引发电气火灾事故。在一些老旧城区的配电网中，由于线路老化、绝缘性能下降，单相接地故障引发的短路和设备烧毁事故时有发生，给居民生活和企业生产带来了极大不便。严重的单相接地故障还可能破坏区域电网的稳定，导致电网电压波动、频率变化，甚至引发连锁反应，造成更大范围的停电事故。在某些极端情况下，单相接地故障可能引发电力系统的振荡，威胁整个电网的安全运行。同时，对于导线落地等单相接地故障，如果配电线路未及时停运，行人和线路巡视人员（特别是夜间）可能发生跨步电压引起的人身电击事故，牲畜也可能遭受电击伤亡。单相接地故障对供电可靠性的影响也不容忽视。发生故障后，一方面需要进行人工选线，这期间未发生故障的配电线路也可能会停电，中断正常供电；另一方面，查找和消除故障需要停运故障线路，导致用户停电时间延长。特别是在庄稼生长期、恶劣气候条件下，以及山区、林区等复杂地区和夜间，故障查找和修复难度加大，将造成长时间、大面积停电，严重影响供电可靠性，给社会经济带来巨大损失。据估算，每年因配电网单相接地故障导致的少供电量可达数亿千瓦时，直接影响供电企业的经济效益和社会效益。因此，快速准确地进行配电网单相接地故障选线与故障定段，对于保障配电网的安全可靠运行、提高供电可靠性、减少停电损失以及保护人身和设备安全具有至关重要的意义。它不仅能及时隔离故障线路，避免故障扩大，还能缩短停电时间，提高电力系统的运行效率和稳定性，满足用户对高质量电能的需求。1.2国内外研究现状在国外，对配电网单相接地故障选线与故障定段的研究起步较早，且取得了一系列成果。早期主要采用基于稳态量的故障选线方法，如零序电流幅值比较法、零序电流相位比较法等。这些方法原理简单，但受系统运行方式、线路参数等因素影响较大，在实际应用中存在一定局限性。随着信号处理技术和人工智能技术的发展，基于暂态量的故障选线与定段方法逐渐成为研究热点。例如，利用小波变换对故障信号进行分析，提取故障特征量，实现故障选线和定位。小波变换具有良好的时频局部化特性，能够有效地捕捉故障暂态信号的突变特征，但小波基函数的选择和分解层数的确定较为困难，且对噪声敏感。近年来，S变换作为一种新的时频分析方法，在配电网故障诊断领域得到了广泛关注。S变换结合了短时傅里叶变换和小波变换的优点，能够在不同频率下具有自适应的时间分辨率，更准确地反映信号的时频特性。国外一些学者利用S变换对配电网单相接地故障零序电流进行分析，通过提取特征频率下的能量等特征量，构建故障选线判据，取得了较好的选线效果。同时，在故障定段方面，将S变换与行波理论相结合，通过分析故障行波的时频特性，实现故障点的精确定位。然而，实际配电网中存在大量的噪声干扰和复杂的电磁环境，如何提高S变换在复杂情况下的抗干扰能力和故障诊断准确性，仍是需要进一步研究的问题。在国内，配电网单相接地故障选线与故障定段的研究也得到了众多学者和科研机构的重视。早期主要借鉴国外的研究成果，并结合国内配电网的实际特点进行改进和完善。随着国内电力技术的快速发展，对该领域的研究逐渐深入，提出了许多具有创新性的方法和技术。除了传统的稳态量和暂态量分析方法外，还涌现出了一些基于人工智能、大数据、物联网等新技术的故障诊断方法。例如，利用人工神经网络对故障数据进行学习和训练，实现故障线路和故障点的智能识别；通过建立故障数据库，运用大数据分析技术挖掘故障特征和规律，提高故障诊断的准确性和可靠性。在基于S变换的研究方面，国内学者也做了大量工作。通过对S变换算法的优化和改进，提高其计算效率和精度；结合配电网的拓扑结构和运行参数，提出了多种基于S变换的故障选线和定段方案，并通过仿真和实验进行验证。然而，目前国内的研究仍存在一些问题和不足。一方面，部分研究成果在实际应用中还存在一定的局限性，如对硬件设备要求较高、算法复杂度较大等，导致难以大规模推广应用；另一方面，对于不同接地方式和复杂运行工况下的配电网，基于S变换的故障诊断方法的适应性和通用性还有待进一步提高。总体而言，国内外在基于S变换的配电网单相接地故障选线与故障定段研究方面取得了一定的进展，但仍存在一些问题需要解决。如在复杂电磁环境下的抗干扰能力、不同接地方式和运行工况下的适应性、算法的实时性和准确性等方面，都还有待进一步深入研究和改进。1.3研究内容与方法本研究基于S变换，深入开展配电网单相接地故障选线与故障定段的研究，主要内容包括以下几个方面：S变换原理及特性分析：详细阐述S变换的基本原理、数学表达式以及时频特性，分析S变换在处理配电网故障信号时的优势和特点，为后续的故障选线与定段研究奠定理论基础。配电网单相接地故障电气特征分析：深入研究不同中性点接地方式下（中性点不接地、经消弧线圈接地、经高电阻接地等）配电网发生单相接地故障时的稳态和暂态电气特征，包括零序电流、零序电压、故障行波等的变化规律，为提取有效的故障特征量提供依据。基于S变换的故障选线方法研究：利用S变换对故障零序电流信号进行时频分析，提取能够表征故障线路的特征量，如特征频率下的能量、相位等。通过对这些特征量的分析和比较，构建合理的故障选线判据，实现准确的故障选线。同时，研究如何提高选线方法在不同运行工况和干扰条件下的可靠性和适应性。基于S变换的故障定段方法研究：结合S变换和故障行波理论，分析故障行波在配电网中的传播特性，通过对故障行波信号的S变换时频分析，确定故障行波到达不同检测点的时间差，进而实现故障点的精确定段。研究如何克服行波传播过程中的衰减、畸变以及噪声干扰等因素对故障定段精度的影响。仿真实验与实际案例验证：利用MATLAB/Simulink等仿真软件搭建配电网单相接地故障仿真模型，对所提出的基于S变换的故障选线与定段方法进行仿真实验，验证其有效性和准确性。同时，收集实际配电网中的故障数据，进行实际案例分析，进一步检验方法在实际工程中的应用效果，并根据实际情况对方法进行优化和改进。在研究方法上，本研究采用理论分析、仿真实验和实际案例验证相结合的方式。首先，通过对S变换原理、配电网故障电气特征以及故障选线与定段方法的理论分析，构建完整的研究体系；然后，利用仿真软件进行大量的仿真实验，模拟不同的故障场景和运行工况，对所提出的方法进行验证和优化；最后，结合实际配电网故障数据，进行实际案例分析，将研究成果应用于实际工程中，检验其可行性和实用性。通过这种多方法结合的研究方式，确保研究成果既具有坚实的理论基础，又能满足实际工程的需求，为提高配电网的安全可靠运行水平提供有效的技术支持。二、配电网单相接地故障特性分析2.1配电网中性点接地方式2.1.1中性点不接地系统中性点不接地系统是指除保护或测量用途的高阻抗接地以外，中性点没有人工接地的非有效接地系统。在正常运行状态下，系统三相电压对称，各相对地电容电流的数值相等而相位相差120^{\circ}，其向量和等于零，地中没有电容电流通过，中性点对地电位为零，即中性点与地电位一致。此时，中性点接地与否对各相对地电压没有任何影响。然而，当中性点不接地系统的各相对地电容不相等时，即使在正常运行状态下，中性点的对地电位便不再是零，通常称此情况为中性点位移，即中性点不再是地电位。这种现象多是由于架空线路排列不对称而又换位不完全所导致。当系统发生单相接地故障时，故障点会流过电容电流，该电流大小为系统正常运行时单相对地电容电流的3倍。由于故障电流相对较小，接地电弧通常能自动熄灭，系统线电压仍能保持对称，三相系统的平衡未遭到破坏，因此可继续运行一段时间，这是中性点不接地系统的最大优点，能在一定程度上提高供电可靠性。例如，在某些对供电连续性要求较高的场合，如医院、数据中心等，中性点不接地系统可以在发生单相接地故障时，维持短时间的供电，为故障排查和处理争取时间。但中性点不接地系统也存在明显的缺点。一方面，未接地两相的对地电压会升高到线电压，这对系统的绝缘水平提出了更高要求，增加了设备成本。若一相接地运行时间过长，可能会造成两相短路，扩大故障范围。另一方面，接地点的电容电流可能引发间歇性弧光接地，产生的谐振过电压幅值可达正常电压的数倍，对设备绝缘构成极大威胁，严重时会导致设备损坏。同时，由于故障电流较小，传统的保护装置难以准确检测和定位故障线路，故障定位难度较大，可能导致故障排查时间延长，影响供电可靠性。2.1.2中性点经消弧线圈接地系统中性点经消弧线圈接地系统是在电网中性点接入一个具有感性的消弧线圈。当系统发生单相接地故障时，故障点的电容电流会使接地点产生电弧，而消弧线圈产生的感性电流可以抵消电弧电流，使通过接地点的电流减小到能够自行熄弧的范围，从而达到熄灭电弧、防止电网故障扩大的目的。消弧线圈的补偿方式主要有欠补偿、全补偿和过补偿三种。欠补偿是指消弧线圈的电感电流小于接地电容电流，此时接地点仍有剩余的电容电流，可能会导致电弧无法完全熄灭，并且在系统运行方式发生变化或切除部分线路时，可能会使系统接近或达到全补偿状态，引发串联谐振过电压，因此欠补偿方式一般较少采用。全补偿是指消弧线圈的电感电流等于接地电容电流，从理论上讲，此时接地点电流为零，电弧可以迅速熄灭，但在实际运行中，由于电网参数的变化和测量误差等因素，很难精确实现全补偿，一旦出现偏差，也容易引发谐振过电压，所以全补偿方式也很少应用。过补偿是指消弧线圈的电感电流大于接地电容电流，接地点存在一定的电感电流，这种方式可以有效避免谐振过电压的产生，并且在系统运行方式变化时，仍能保证可靠的补偿效果，因此在实际工程中应用最为广泛。在单相接地故障发生瞬间，暂态过程较为复杂。暂态零序电流由电容电流的自由分量和强制分量组成，其中自由分量衰减较快，强制分量衰减较慢。暂态电容电流中包括自由分量和强制分量，在相电压接近最大值瞬间单相接地过程中，暂态电容电流比流过消弧线圈的暂态电感电流大很多，暂态电感电流可忽略不计，因此在同一电网中，即使中性点经消弧线圈接地，其过渡过程与中性点不接地情况下近似相同。但随着故障电阻的变化，系统暂态特性也会发生改变。当故障电阻较小时，系统处于欠阻尼状态；随着故障电阻的增大，系统进入过阻尼状态；最后又过渡到另一种等值电路下的欠阻尼状态。此外，电缆线路的引入会使暂态电流频谱主成分向低频段移动，故障电阻对含电缆线路的谐振接地系统单相接地故障暂态电流频谱特征影响很大，电阻增大，频谱主成分向低频段移动，当故障电阻增大到一定数值时，频谱主成分趋于工频。2.1.3中性点经电阻接地系统中性点经电阻接地系统是将变压器或发电机的中性点通过一个电阻接地。根据电阻阻值的大小，可分为高电阻接地和低电阻接地。在高电阻接地系统中，电阻阻值较大，一般为几百欧姆到数千欧姆，其主要目的是限制故障电流，同时利用接地电流产生的零序电压和零序电流信号，实现故障检测和保护。高电阻接地系统能有效降低故障时的对地过电压，对设备绝缘要求相对较低，并且可以避免谐振过电压的产生。在一些对供电可靠性要求较高，同时又希望限制故障电流的场合，如城市配电网的部分区域，高电阻接地系统得到了应用。低电阻接地系统的电阻阻值较小，一般为几欧姆到几十欧姆，主要作用是在发生单相接地故障时，使故障电流足够大，以便快速切除故障线路，提高系统的安全性。当系统发生单相接地故障时，故障电流通过小电阻流入大地，断路器迅速跳闸，从而快速切除故障。这种方式能够有效抑制谐振过电压，并且设备绝缘水平要求较低，但由于故障时会产生较大的故障电流，可能会对通信设备产生电磁干扰，影响通信设备的正常运行。在一些工业企业的配电网中，由于对供电连续性要求相对较低，而更注重系统的安全性和故障切除的快速性，低电阻接地系统得到了广泛应用。2.2单相接地故障电气量特征分析2.2.1稳态电气量特征在配电网发生单相接地故障时，稳态电气量会发生明显变化。以中性点不接地系统为例，当发生单相接地故障时，故障相电压降为零，非故障相电压升高到线电压，此时母线处的零序电压大小等于电网的相电压。各线路的零序电流大小等于本线路的对地电容电流，方向为从母线流向线路。故障线路的零序电流等于所有非故障线路零序电流之和，方向与非故障线路相反，从线路流向母线。因此，通过比较各线路零序电流的大小和方向，可以初步判断故障线路。但这种方法在实际应用中存在一定局限性，当系统运行方式发生变化，如线路投切、负荷变化等，会导致线路对地电容发生改变，从而影响零序电流的大小和分布，降低选线的准确性。在中性点经消弧线圈接地系统中，由于消弧线圈的补偿作用，故障线路和非故障线路的零序电流幅值差异减小，且方向可能不再具有明显的特征。在过补偿情况下，故障线路和非故障线路的零序电流方向都为从母线流向线路，仅通过零序电流的大小和方向难以准确选线。不过，此时零序功率方向会发生变化，故障线路的零序功率方向是由线路指向母线，而非故障线路的零序功率方向是由母线指向线路，可利用零序功率方向来进行故障选线，但该方法也容易受到过渡电阻、消弧线圈调谐不准确等因素的影响。2.2.2暂态电气量特征配电网单相接地故障暂态过程持续时间较短，一般为几个工频周期，但其中包含了丰富的故障信息。在暂态过程中，暂态电流和暂态电压的变化特点和频谱特性对于故障选线和定段具有重要意义。当发生单相接地故障瞬间，系统的电容和电感元件会发生能量的重新分配，产生暂态电流和暂态电压。暂态电流主要由电容电流的自由分量和强制分量组成，其中自由分量衰减较快，其频率较高，一般在数千赫兹以上；强制分量衰减较慢，频率接近工频。暂态电压也包含高频分量和低频分量，高频分量主要反映故障瞬间的突变信息，低频分量则与系统的稳态运行状态相关。通过对暂态电流和暂态电压进行频谱分析，可以发现其频谱特性具有明显的特征。在故障暂态初期，高频分量的幅值较大，随着时间的推移，高频分量逐渐衰减，低频分量逐渐占据主导。不同线路的暂态电流和暂态电压在频谱特性上存在差异，故障线路的暂态电流和暂态电压在某些特征频率下的幅值和相位与非故障线路不同。例如，在中性点不接地系统中，故障线路暂态零序电流的高频分量幅值通常大于非故障线路；在中性点经消弧线圈接地系统中，虽然消弧线圈对暂态电流有一定的补偿作用，但故障线路和非故障线路在暂态电流的频谱特性上仍存在可区分的特征。利用这些特征，可以提取有效的故障特征量，用于故障选线和定段。三、S变换原理及其在故障信号分析中的应用3.1S变换基本原理3.1.1S变换的定义与公式推导S变换（S-transform）是一种时频分析方法，由R.G.Stockwell于1996年提出，它结合了短时傅里叶变换（STFT）和小波变换的优点，能够在时频域对信号进行有效的分析。S变换的定义如下：设x(t)为时间t的函数，其S变换S(\tau,f)定义为：S(\tau,f)=\int_{-\infty}^{\infty}x(t)\frac{|f|}{\sqrt{2\pi}}e^{-\frac{(t-\tau)^2f^2}{2}}e^{-j2\pif(t-\tau)}dt其中，\tau为时移因子，控制高斯窗函数在时间轴上的位置；f为频率；x(t)为分析信号；e^{-\frac{(t-\tau)^2f^2}{2}}是高斯窗函数，其标准差为\frac{1}{|f|}，这使得窗函数的宽度会随着频率f的变化而自适应调整，在低频时窗宽较大，具有较高的频率分辨率；在高频时窗宽较小，具有较高的时间分辨率。下面对S变换公式进行详细推导。首先回顾傅里叶变换，对于信号x(t)，其傅里叶变换X(f)为：X(f)=\int_{-\infty}^{\infty}x(t)e^{-j2\pift}dt为了实现对信号的时频局部化分析，在傅里叶变换的基础上引入窗函数w(t-\tau)，得到短时傅里叶变换（STFT）：STFT(\tau,f)=\int_{-\infty}^{\infty}x(t)w(t-\tau)e^{-j2\pift}dtSTFT的窗函数w(t-\tau)是固定的，无法根据信号频率特性自适应调整，而S变换采用了高斯窗函数w(t-\tau)=\frac{|f|}{\sqrt{2\pi}}e^{-\frac{(t-\tau)^2f^2}{2}}，将其代入STFT公式，即可得到S变换的表达式：S(\tau,f)=\int_{-\infty}^{\infty}x(t)\frac{|f|}{\sqrt{2\pi}}e^{-\frac{(t-\tau)^2f^2}{2}}e^{-j2\pift}dt进一步变形可得：S(\tau,f)=\int_{-\infty}^{\infty}x(t)\frac{|f|}{\sqrt{2\pi}}e^{-\frac{(t-\tau)^2f^2}{2}}e^{-j2\pif(t-\tau)}e^{-j2\pif\tau}dt=e^{-j2\pif\tau}\int_{-\infty}^{\infty}x(t)\frac{|f|}{\sqrt{2\pi}}e^{-\frac{(t-\tau)^2f^2}{2}}e^{-j2\pif(t-\tau)}dt对于离散信号x[n]，其离散傅里叶变换（DFT）为X[k]，采样时间间隔为T，则离散S变换可表示为：S[n,k]=\sum_{m=0}^{N-1}x[m]\frac{|k|}{\sqrt{2\pi}}e^{-\frac{(m-n)^2k^2}{2N^2}}e^{-j2\pi\frac{k}{N}(m-n)}其中，N为离散信号的长度，n=0,1,\cdots,N-1，k=0,1,\cdots,N-1。通过上述推导过程，明确了S变换从时域到频域的变换过程，为后续利用S变换分析信号的时频特性奠定了基础。3.1.2S变换的时频特性分析S变换的时频特性主要体现在其对不同频率信号具有自适应的时间分辨率和频率分辨率。从S变换的定义可知，其高斯窗函数的宽度与频率成反比，即\sigma=\frac{1}{|f|}。这意味着在低频段，高斯窗函数较宽，对信号的时间分辨率较低，但频率分辨率较高，能够更准确地分辨低频信号的频率成分；在高频段，高斯窗函数较窄，对信号的时间分辨率较高，能够捕捉到高频信号的快速变化，但频率分辨率相对较低。以一个包含低频和高频成分的合成信号为例，对其进行S变换分析。合成信号x(t)由低频正弦信号x_1(t)=A_1\sin(2\pif_1t)和高频正弦信号x_2(t)=A_2\sin(2\pif_2t)组成，即x(t)=x_1(t)+x_2(t)，其中A_1和A_2分别为低频和高频信号的幅值，f_1=50Hz，f_2=500Hz。对x(t)进行S变换后，得到时频图。在时频图中，可以清晰地看到低频信号x_1(t)的能量主要集中在频率f_1附近，且在时间轴上的分布较为平滑，这是因为低频时高斯窗函数较宽，对低频信号的时间分辨率较低，但能较好地分辨其频率；高频信号x_2(t)的能量主要集中在频率f_2附近，在时间轴上的分布较为集中，能够准确地定位高频信号出现的时间，这是由于高频时高斯窗函数较窄，时间分辨率较高。与短时傅里叶变换相比，短时傅里叶变换使用固定的窗函数，无法根据信号频率自适应调整窗宽，因此在时频分辨率上存在局限性，对于不同频率的信号不能同时兼顾良好的时间分辨率和频率分辨率。而S变换通过自适应调整窗宽，能够在不同频率下都能较好地反映信号的时频特性，在时频分辨率上具有明显优势。S变换在时频域的分辨率特性使其非常适合分析非平稳信号，如配电网中的故障信号。在配电网发生单相接地故障时，故障信号包含丰富的暂态信息，既有低频的稳态分量，也有高频的暂态分量，S变换能够有效地分析这些不同频率成分的时频特性，为故障特征提取和故障诊断提供有力支持。3.2S变换在配电网故障信号分析中的优势3.2.1与傅里叶变换、小波变换的对比傅里叶变换（FT）是一种将时域信号转换为频域信号的分析方法，其核心思想是将任何周期函数或非周期函数分解为一系列不同频率的正弦和余弦函数的叠加。傅里叶变换在分析平稳信号时具有重要作用，能够清晰地展示信号的频率组成。然而，傅里叶变换缺乏时间信息，它将信号整体进行变换，无法反映信号在不同时刻的频率变化情况。对于配电网中的故障信号，其具有明显的非平稳特性，故障发生瞬间信号会发生突变，傅里叶变换难以准确捕捉这些瞬态变化信息，因此在处理配电网故障信号时存在较大局限性。小波变换（WT）是一种时频分析方法，它通过将信号分解为一系列小波函数的线性组合，实现对信号的多分辨率分析。小波变换的窗函数可以根据尺度因子和平移因子进行调整，具有良好的时频局部化特性，能够有效地处理非平稳信号。但是，小波变换在选择小波基函数时存在一定困难，不同的小波基函数对信号分析的结果可能有较大差异，且小波变换的分解层数需要根据具体问题进行选择，缺乏明确的理论指导。S变换结合了傅里叶变换和小波变换的优点。与傅里叶变换相比，S变换能够同时提供信号的时间和频率信息，通过高斯窗函数对信号进行局部化分析，克服了傅里叶变换无法处理非平稳信号的缺陷。与小波变换相比，S变换无需选择小波基函数，其高斯窗函数具有自适应性，能够根据信号频率自动调整窗宽，在不同频率下都能获得较好的时频分辨率。在处理配电网单相接地故障信号时，S变换能够更准确地分析故障信号的时频特性。故障信号中的暂态分量包含丰富的故障信息，S变换可以清晰地展示暂态分量在不同时刻的频率变化，为故障特征提取提供更全面的信息。而傅里叶变换无法有效分析暂态分量，小波变换在选择合适的小波基函数和分解层数时可能面临困难，影响对故障信号的分析效果。3.2.2对故障信号特征提取的有效性为了验证S变换对故障信号特征提取的有效性，以实际配电网单相接地故障信号为例进行分析。通过在MATLAB/Simulink中搭建配电网仿真模型，模拟不同类型的单相接地故障，获取故障零序电流信号。对故障零序电流信号进行S变换，得到时频图。在时频图中，可以观察到故障信号在不同频率和时间上的能量分布情况。故障发生瞬间，时频图上会出现明显的能量集中区域，且不同故障类型和故障位置对应的能量分布特征存在差异。例如，在中性点不接地系统中发生金属性单相接地故障时，故障零序电流的高频分量在时频图上表现为在故障发生时刻出现高频能量尖峰，且随着时间推移，高频能量逐渐衰减；而在经消弧线圈接地系统中发生单相接地故障时，由于消弧线圈的补偿作用，故障零序电流的高频分量相对较弱，但在时频图上仍能观察到与故障相关的特征频率成分。通过对时频图的分析，可以提取出能够表征故障的特征量，如特征频率下的能量、相位等。将这些特征量作为故障诊断的依据，能够有效地识别故障线路和故障类型。与传统的基于时域或频域分析的故障特征提取方法相比，基于S变换的特征提取方法能够充分利用故障信号的时频信息，提取的特征量更加全面、准确，从而提高故障诊断的准确性和可靠性。四、基于S变换的配电网单相接地故障选线方法研究4.1基于S变换的故障选线判据构建4.1.1零序电流S变换能量特征分析在配电网发生单相接地故障时，零序电流包含了丰富的故障信息。对故障线路和非故障线路的零序电流进行S变换，能够将时域的零序电流信号转换为时频域信号，从而更清晰地展现其能量在不同频率和时间上的分布特征。以中性点不接地系统为例，当发生单相接地故障时，故障线路的零序电流等于所有非故障线路对地电容电流之和。对故障线路和非故障线路的零序电流进行S变换后，在时频图上可以观察到，故障线路零序电流在某些特征频率下的能量明显高于非故障线路。这是因为故障线路零序电流的幅值相对较大，经过S变换后，其在特征频率处的能量也相应增强。在中性点经消弧线圈接地系统中，虽然消弧线圈的补偿作用会使故障线路和非故障线路的零序电流幅值差异减小，但通过S变换分析仍能发现两者在能量分布上的不同。故障线路零序电流的暂态过程中，包含了自由分量和强制分量，自由分量的频率较高，在S变换时频图上，这些高频分量对应的能量分布与非故障线路存在差异。例如，在故障发生后的短时间内，故障线路零序电流的高频能量可能会出现一个尖峰，而后逐渐衰减，而非故障线路的高频能量变化相对较为平缓。进一步分析不同故障位置和过渡电阻对零序电流S变换能量特征的影响。当故障位置靠近母线时，故障线路零序电流的幅值相对较大，S变换后在特征频率下的能量也较大；随着故障位置远离母线，零序电流幅值逐渐减小，能量也相应降低。过渡电阻的存在会使故障电流减小，从而影响零序电流S变换的能量分布。当过渡电阻增大时，故障线路零序电流在特征频率下的能量会减弱，但即使在高阻接地故障情况下，通过S变换仍能在时频图上观察到故障线路与非故障线路能量分布的差异。4.1.2故障选线判据的建立根据上述对零序电流S变换能量特征的分析，建立以能量最大或能量变化率为判据的故障选线方法。以能量最大为判据时，定义各线路零序电流在S变换后某一特征频带内的能量为：E_i=\sum_{f=f_1}^{f_2}\sum_{t=t_1}^{t_2}|S_i(t,f)|^2其中，E_i为第i条线路零序电流在特征频带内的能量，S_i(t,f)为第i条线路零序电流的S变换结果，f_1和f_2为特征频带的下限和上限频率，t_1和t_2为分析的时间区间。在实际应用中，通过比较各线路的能量E_i大小，能量最大的线路判定为故障线路。例如，在一个包含多条线路的配电网中，对各线路零序电流进行S变换并计算能量，若线路L_3的能量E_3在所有线路中最大，则判定线路L_3为故障线路。以能量变化率为判据时，定义各线路零序电流在S变换后某一特征频带内的能量变化率为：\DeltaE_i=\frac{E_{i2}-E_{i1}}{E_{i1}}其中，\DeltaE_i为第i条线路零序电流在特征频带内的能量变化率，E_{i1}和E_{i2}分别为故障前和故障后第i条线路零序电流在特征频带内的能量。在故障发生后，故障线路的能量变化率通常会大于非故障线路。通过比较各线路的能量变化率\DeltaE_i，能量变化率最大的线路判定为故障线路。例如，在某一配电网中，故障发生后，线路L_2的能量变化率\DeltaE_2显著大于其他线路，据此可判定线路L_2为故障线路。通过理论分析和仿真验证，这两种判据在不同中性点接地方式下的配电网中都具有较好的故障选线性能。能量最大判据在故障信号较强、干扰较小时，能够准确地选出故障线路；能量变化率判据则对故障的快速响应能力较强，在故障初期就能有效地识别故障线路。同时，为了提高选线的可靠性和准确性，可以综合考虑这两种判据，当两种判据的结果一致时，进一步增强对故障线路判定的可信度；当两种判据结果不一致时，可以结合其他辅助信息，如零序电流的相位信息等，进行综合判断。4.2故障选线算法实现与流程设计4.2.1算法步骤详细描述基于S变换的故障选线算法具体步骤如下：数据采集：利用安装在配电网各线路上的电流互感器等测量装置，实时采集各线路的零序电流数据。采样频率应根据配电网的实际情况和故障信号的频率特性进行合理选择，一般要求采样频率足够高，以准确捕捉故障信号的变化。例如，对于50Hz工频的配电网，采样频率可选择为6400Hz或更高，以确保能够完整地获取故障暂态过程中的零序电流信息。S变换计算：将采集到的各线路零序电流数据进行S变换。根据S变换的离散形式公式，对每一条线路的零序电流序列进行计算，得到各线路零序电流在不同时间和频率下的S变换结果S_i(t,f)。在计算过程中，可利用快速傅里叶变换（FFT）算法来提高计算效率，减少计算时间。例如，在MATLAB编程实现中，可以调用FFT函数进行快速计算。判据计算：根据前面建立的故障选线判据，计算各线路的能量E_i或能量变化率\DeltaE_i。若采用能量最大判据，按照公式E_i=\sum_{f=f_1}^{f_2}\sum_{t=t_1}^{t_2}|S_i(t,f)|^2计算各线路在特征频带内的能量；若采用能量变化率判据，则先计算故障前和故障后各线路在特征频带内的能量E_{i1}和E_{i2}，再根据公式\DeltaE_i=\frac{E_{i2}-E_{i1}}{E_{i1}}计算能量变化率。故障线路判定：比较各线路的能量E_i或能量变化率\DeltaE_i大小，能量最大或能量变化率最大的线路判定为故障线路。在判定过程中，可设置一定的阈值，以避免因干扰等因素导致的误判。例如，当某线路的能量或能量变化率与其他线路相比超过一定阈值时，才判定为故障线路；若各线路的能量或能量变化率差异不明显，可进一步分析其他特征或重新进行数据采集和计算。4.2.2算法流程图绘制基于S变换的故障选线算法流程图如图1所示：@startumlstart:数据采集;:进行S变换计算;:计算能量E或能量变化率ΔE;if(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@endumlstart:数据采集;:进行S变换计算;:计算能量E或能量变化率ΔE;if(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:数据采集;:进行S变换计算;:计算能量E或能量变化率ΔE;if(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:进行S变换计算;:计算能量E或能量变化率ΔE;if(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:计算能量E或能量变化率ΔE;if(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@endumlif(采用能量最大判据?)then(是):比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:比较各线路能量E大小;:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:能量最大的线路判定为故障线路;else(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@endumlelse(否):比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:比较各线路能量变化率ΔE大小;:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@enduml:能量变化率最大的线路判定为故障线路;endif:输出故障线路结果;stop@endumlendif:输出故障线路结果;stop@enduml:输出故障线路结果;stop@endumlstop@enduml@enduml图1中，首先进行数据采集，获取各线路零序电流数据；然后对数据进行S变换计算，得到S变换结果；接着根据选择的判据（能量最大判据或能量变化率判据）计算相应的特征量（能量E或能量变化率\DeltaE）；通过比较各线路的特征量大小，判定故障线路；最后输出故障线路结果。该流程图清晰地展示了故障选线算法的执行逻辑和流程，便于理解和实现。4.3仿真分析与验证4.3.1仿真模型搭建利用MATLAB/Simulink软件搭建配电网仿真模型，以验证基于S变换的故障选线方法的准确性和可靠性。仿真模型的结构如图2所示：@startumlpackage"电源"assource{component"三相交流电源"aspower}package"中性点接地方式"asgrounding{component"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlpackage"电源"assource{component"三相交流电源"aspower}package"中性点接地方式"asgrounding{component"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"三相交流电源"aspower}package"中性点接地方式"asgrounding{component"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@enduml}package"中性点接地方式"asgrounding{component"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlpackage"中性点接地方式"asgrounding{component"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"中性点不接地"asungroundedcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"中性点经消弧线圈接地"asarc_coil_groundedcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"中性点经电阻接地"asresistance_grounded}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@enduml}package"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlpackage"线路"aslines{component"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"线路L1"asline1component"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"线路L2"asline2component"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"线路L3"asline3component"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"线路L4"asline4}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@enduml}package"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlpackage"负荷"asloads{component"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"负荷1"asload1component"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_point-->lineslines-->ctct-->s_transforms_transform-->fault_selection@endumlcomponent"负荷2"asload2component"负荷3"asload3component"负荷4"asload4}package"故障设置"asfault{component"故障点"asfault_point}package"测量与分析"asmeasurement{component"电流互感器"asctcomponent"S变换模块"ass_transformcomponent"故障选线模块"asfault_selection}power-->groundinggrounding-->lineslines-->loadslines-->fault_pointfault_poin