一种行波信号的参数化识别方法
【专利摘要】本发明属于电力系统领域，涉及一种行波信号的参数化识别方法。该方法根据行波信号的统计特性建立行波信号的参数化模型；利用自适应粒子滤波算法和优化方法估计出表示不同故障情况下行波特征的参数信息和时频信息,提取能够全面反映行波信号特点的时域和时频域关联特征量组；联合分析上述关联特征量组，提取行波故障特征，突显故障信息，准确可靠地感知强电磁干扰下微弱故障行波信号和复杂折反射行波信号。本发明方法利用参数化时频分析的方法来识别电网故障行波，提出行波信号的新型识别参数,可有效提高受扰行波信号识别的可靠性和准确性,为构建新型故障行波保护与定位原理提供参数支持。
【专利说明】
一种行波信号的参数化识别方法
技术领域
[0001] 本发明属于电力系统继电保护领域，涉及一种行波信号的参数化识别方法。
【背景技术】
[0002] 目前，行波故障定位法在理论上具有很高的定位精度，使其获得了广泛的研究和 应用。故障行波信号突变点的时刻及其位置都代表了具体的故障信息，因此，准确检测故障 行波信号的突变点就成为行波故障定位的关键。基于小波分析、希尔伯特黄变换等的非参 数时频分析方法，在行波识别中已获得了广泛的应用，并取得了较好的故障定位效果。但这 些非参数化的时频分析方法在实际应用中还存在一定的局限性，如电压过零故障或高阻接 地故障等微弱故障行波信号，以及行波传输过程中的衰减对行波奇异性的削弱；线路雷击 未故障时雷击干扰的识别，以及雷击故障中雷击点和闪络点不一致时，闪络点的精确定位 需要识别和检测雷击侧的第二个行波波头(故障点的反射波）；在单端行波故障定位法中， 故障点反射波与对端母线反射波的准确识别；行波信号在传输和采样过程中，噪声的干扰 使得有用信号不能被有效提取。在这些影响因素下，单纯以时域或时频域特征量提取的行 波识别方法，其可靠性和准确性难以保障，严重制约了行波故障定位的精度和准确度。
【发明内容】
[0003] 为克服已有技术的不足，解决现有行波信号识别方法在噪声干扰和微弱故障情况 下其准确性和可靠性差的技术问题，提出了一种行波信号识别的参数化时频分析方法。本 发明根据行波信号的统计特性，通过参数化建模的方法，利用丰富的参数信息和时频谱信 息提取能够全面反映信号特点的时域和时频域关联特征量组，联合分析上述关联特征量 组，以此提取行波故障特征，准确可靠地感知强电磁干扰下微弱故障行波信号和复杂折反 射行波信号，有效解决了现有行波信号识别可靠性和准确性差等问题，有助于进一步提高 行波故障定位的精确度和实用性。
[0004] 本发明解决上述技术问题的技术方案包括以下步骤： 1) 根据行波信号的脉冲特性和非平稳特性，构建由?稳定分布驱动的时变自回归模 型(?稳定分布-TVAR模型）； 2) 针对模型的非高斯性以及时变参数突变的延迟性，提出以参数噪声统计特性估计 为核心、带时变遗忘因子的自适应粒子滤波算法，用于行波信号建模的参数估计中； 3) 根据估计出的模型系数计算分数低阶时变功率谱密度函数，提取以参数幅值特征 的时域特征量和以分数低阶时变功率谱表征的时频域特征量； 4) 根据估计出的《稳定分布的参数值提取行波信号的时域特征量； 5) 根据时域和时频域的多种关联特征量组，准确提取能够全面反映行波故障信号的 特征量，实现复杂行波信号和微弱行波信号的可靠准确辨识。
[0005] 上述步骤1)中，受扰行波信号的a稳定分布-TVAR模型：
其中，心(1_=;1足-，-:^)是时变的自回归系数，|*代表模型的阶次，1^卷也#:#>是服 从《稳定分布的随机噪声，且〃、:r、#、芦分别指特征指数、分散系数、对称参数、位置参数， 它们的取值能完全控制夂稳定分布特性。
[0006] 上述步骤2)中，为了估计出模型中的参数，需建立如下的状态空间模型：
其中，尤二:读:名；f，且::尤：=[^以产为系统真实的状态变量， € = 为系统扩展的状态变量，以及观测变量&=1^1，系统矩阵
，观测矩阵梦*[1^---,用，峨％^是参数扰动量即人 工噪声，是观测高斯白噪声。
[0007] 基于上述状态空间模型，利用以参数噪声统计特性估计为核心、带时变遗忘因子 的自适应粒子滤波算法，估计出行波信号模型中的参数。
[0008] 上述步骤3)中，根据估计出TVAR模型的自回归系数的奇异值即行波信号的突变点 来初步判断雷击冲击或故障冲击。
[0009] 上述步骤3)中，根据估计出TVAR模型的自回归系数和從稳定分布的参数值计算分 数低阶瞬时功率谱密度函数:
由上式瞬时频率所组成的时频图中，依据行波波头在时频图中表现的高频突变，即可 实现对行波波头的准确检测。
[0010] 上述步骤4)中，提取以a：稳定分布特征指数和分散系数、以及峭度所表示的信号 脉冲水平和故障程度的时域特征量。具体来讲： 特征指数《反映行波信号整体脉冲水平。其值越小，说明行波信号的脉冲性越强，在无 噪声或弱噪声情况下，线路故障越严重;否则线路无故障或故障轻微。该特征量具有较强的 鲁棒性，其值不会受少量干扰脉冲信号的影响，能很好地反映线路故障程度。
[0011] 分散系数7反映行波信号整体幅值的变化，其值越大，说明行波信号整体幅值变 化越大，无论是故障本身引起的变化大，还是雷击干扰所引起的变化大。利用分散系数的这 一特点，可以识别雷击干扰和雷击故障。当线路雷击未造成故障时，雷击如后电流整体幅值 变化不大；当线路雷击造成故障时，雷击前后电流整体幅值变化大。不同雷击情况下电流整 体幅值变化不同的特征可通过分散系数值的变化来表征，即：若分散系数值变化不大，保持 在某一稳定水平，则雷击未造成故障;若分散系数值变化大，则雷击造成故障。当然，在无干 扰情况下，该特征参数与故障程度也是紧密相关的。
[0012] 峭度是描述波形尖峰量的量纲，是一种经典的表征脉冲特性的指标参数。其性能 同特征指数类似，峭度值变大时，则行波信号的脉冲性越强；反之，则行波信号的脉冲性 越弱。但与特征指数相比，峭度对信号脉冲的变化非常敏感，任意单个脉冲的出现都会使峭 度有一个大幅度的增加。在故障信号微弱（如轻微故障）时，这一特征参数也能较好地捕捉 到信号的突变点，即辨识出行波波头。因此，该特征量能准确反映轻微故障的发生。
[0013] 与现有技术相比，本发明所具有的有益效果为:本发明依据行波信号的统计特性， 从时域的角度，利用统计信号处理方法和人工智能技术对行波信号进行分析和建模，从时 域和时频域提取暂态行波信号的参数特征量和瞬时频率特征量，而不仅仅凭借单一故障行 波时域或时频域特征来辨识行波信号。本发明方法是一种典型的参数化时频分析方法，综 合了时域、时频域以及参数化分析方法的优点：丰富的参数信息和较高的时频分辨率。本发 明方法能在保留原始信号特征的同时，有效滤掉混杂在行波信号中的各种噪声;通过综合 分析时域和时频域的多种关联特征量组，实现行波信号的准确辨识。本发明方法的提出，能 有效提高行波故障定位的准确性和可靠性，提高行波保护的实用化，为构建新型故障行波 保护与定位原理提供参数支持。
[0014] 下面结合附图对本发明做出进一步的说明。
【附图说明】
[0015] 图1为本发明行波信号识别的流程图。
【具体实施方式】
[0016] 图1为本发明行波信号识别的流程图，搭建电网EMTP仿真分析模型，设置不同故障 点、故障类型以及故障程度，全面总结不同故障情况下行波信号的传输特性;在学校高压实 验室搭建典型模拟实验系统，模拟测试雷击、开关操作及污闪产生的干扰暂态行波信号的 波形特征，通过统计特性分析，比对不同故障情况下行波信号与干扰暂态行波信号的波形 特征差异，提出辨识故障行波信号的判据;选择株洲liokv线路，采用高速存储示波器和申 请者所在课题组最新研制的10MHZ暂态记录分析仪测试线路弧光故障、高阻故障和金属接 地故障时，不同故障点位置产生的原始行波波形，与线路干扰暂态行波信号进行波形相似 度比对，寻找波形差异规律，进一步改进与完善理论研究结果，为行波信号的有效辨识提供 完备的理论依据。
[0017] 采样的行波信号是被白噪声干扰、具有脉冲特性的非平稳时间序列，它可被描述 为： 式中，(泥为采样长度)是被白噪声干扰的行波序列，^^是纯净的 行波序列，为观测白噪声。
[0018] 根据行波信号的非平稳性和突变性，去噪后纯净的行波序列4可通过?稳定分 布-TVAR模型所描述 式中，气是时变的自回归系数』代表模型的阶次，办，爲#是服 从《稳定分布的随机噪声。H#、芦分别指特征指数、分散系数、对称参数、位置参数，它 们的取值能完全控制稳定分布特性。
[0019] 为了估计出模型中的参数，需建立如下的状态空间模型：
其中，，且:笔为系统真实的状态变量， 見为系统扩展的状态变量，以及观测变量心，系统 矩阵
，观测矩阵葬V叫，瑪氣乂)'是参数扰动量 即人工噪声，：4~網<1^:)是观测高斯白噪声。
[0020] 基于上述状态空间模型，利用两类以参数噪声统计特性估计为核心、带时变遗忘 因子的自适应粒子滤波算法，估计出行波信号模型中的参数。根据估计出的模型参数，提取 以TVAR模型的自回归系数、〃稳定分布中的分布参数a、:：f和峭度所组成的时域特征量组， 以及由瞬时频率所表示的时频特征量。综合分析上述关联特征量组，对比不同影响因素和 故障情况下特征量的差异，准确提取能够全面反映行波故障信息的特征量，实现行波信号 的准确识别。如，线路故障严重时，则特征指数?值较小，否则，线路未故障或轻微故障时， 则特征指数《值较大；当线路雷击未造成故障时，雷击前后电流整体幅值变化不大，则分散 系数:r值变化不大，保持在某一稳定水平，否则，当线路雷击造成故障时，雷击前后电流整 体幅值变化大，则分散系数f值变化大;在故障信号微弱如轻微故障时，则可通过峭度值的 变化来进行判断。
【主权项】
1. 一种行波信号的参数化识别方法，包括如下步骤： 1) 根据行波信号的脉冲特性和非平稳特性，构建由《稳定分布驱动的时变自回归模 型(〃稳定分布-TVAR模型）； 2) 针对模型的非高斯性以及时变参数突变的延迟性，提出以参数噪声统计特性估计 为核心、带时变遗忘因子的自适应粒子滤波算法，用于行波信号建模的参数估计中； 3) 根据估计出的模型系数计算分数低阶时变功率谱，提取以参数幅值特征的时域特 征量和以分数低阶时变功率谱表征的时频域特征量； 4) 根据估计出的α稳定分布参数值提取行波信号的时域特征量； 5) 根据时域和时频域的多种关联特征量组，准确提取能够全面反映行波故障信号的 特征量，实现复杂行波信号和微弱行波信号的可靠准确辨识。2. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤1)中，行波 信号的α稳定分布-TVARM刑先.其中，是时变的自回归系数，Jr代表模型的阶是服 从《稳定分布的随机噪声，且〃、::Γ、及、私分别指特征指数、分散系数、对称参数、位置参数， 它们的取值能完全控制&稳定分布特性； TVAR模型自回归系数吒》的动态变化情况能很好地反映信号的动态特性，即自回归系数 的突变或缓慢变化分别对应行波信号的突变(如雷击冲击或故障冲击)或缓慢变化，由 此可提取以模型自回归系数奇异值表示的时域特征量;而且，根据时变自回归系数％▲的值 可以计算出信号的瞬时频率，由此可提取由瞬时频率所组成的时频特征量； 在α稳定分布中，与行波信号特征相关的两参数《和Γ分别表征信号的脉冲水平和故 障程度;g卩，《值越小，信号的整体脉冲特性越强；7越大，信号的整体幅值变化越大； 这两个参数都与故障程度和故障信息密切相关，因此，可以将它们作为衡量信号脉冲 特性和故障程度的参数特征量。3. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤2)中，模型 参数估计方法如下：- I分布-TVAR模型所有待估计的参数；为了准确估i+m樽型中白冋咁系教以及纯净杆被信号，诜取状杰夺投If， 且 ? 扩 其中，系统矩 是参数扰动量即y基于上述状态空间模型，利用自适应滤波算法即可估计出模型状态和参数。4. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤3)中提取以 模型自回归系数％*奇异值表示的时域特征量，即：自回归系数的奇异值点对应行波信号的 突变，如雷击冲击或故障冲击。5. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤3)中信号的 α谱密度函数的估计如下：由上式瞬时频率所组成的时频图中，行波波头会在该时频图中表现为高频突变，由此 可以通过频率突变点的时间位置确定行波到达时间，实现对行波波头的准确检测。6. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤4)中根据估 计估计的特征指数?来判断线路故障程度，即其值越小，说明行波信号的脉冲性越强，在无 噪声或弱噪声情况下，线路故障越严重，反之亦然。7. 根据权利要求1所述行波信号的参数化识别方法，其特征在于:所述步骤4)中根据估 计的分散系数/来区分雷击干扰和雷击故障，即：若其值变化不大，保持在某一稳定水平， 则雷击未造成故障;若其值变化大，则雷击造成故障。
【文档编号】G06K9/00GK105893976SQ201610246191
【公开日】2016年8月24日
【申请日】2016年4月20日
【发明人】席燕辉, 李泽文, 赵廷, 郭田田
【申请人】长沙理工大学