一种基于录波文件基准通道的二次设备时标对齐方法技术

本发明专利技术公开了一种基于录波文件基准通道的二次设备时标对齐方法，选取不同设备间的同电气节点处的录波通道，设定参考录波信号和待对齐信号，利用总谐波畸变率分布特点提取待对齐信号的特征信号波段，建立短时归一化相关系数模型，定义时间识别指标，计算并修正时标，最后完成不同二次设备录波波形及SOE时标对齐，准确还原电网事故关键事件状态序列，为事故回放、分析和处理提供参考依据。本发明专利技术所述方法不受不同设备录波极性影响，时间偏差计算准确度高，具有良好的抗噪声能力。

A time scale alignment method of secondary equipment based on the reference channel of recording file

【技术实现步骤摘要】
一种基于录波文件基准通道的二次设备时标对齐方法
本专利技术涉及电力系统继电保护领域，尤其涉及一种基于录波文件基准通道的二次设备时标对齐方法。
技术介绍
现代电力系统的覆盖范围非常广泛，电网运行瞬息万变，为全面地、实时地、准确地监控电力系统的运行状态，以便分析事故发展的过程与原因，电力系统各系统之间需采用一个统一的时间基准。外部同步时钟信号中断或异常时，变电站自动化装置仅依靠内部自带时钟，其固有误差难以避免，而且误差会随着运行时间的增加而累积，从而失去准确的时间统一作用。难以准确描述电力系统的事件顺序和发展过程，无法给电网事故分析提供有效的分析依据。因此，统一精确的时间是保证电力系统安全运行，提高运行水平的一个重要措施。一旦外部同步时钟信号异常或缺失，会造成电力系统二次设备缺少统一的时间基准，其提供的SOE记录和录波文件不可避免地会发生时间顺序错位，仅靠人工比对时间序列，不仅耗时耗力，而且修正精度低。因此迫切需要实现对于不同装置故障录波文件的时标对齐功能，并利用计算出的二次设备间的时间误差修正相关SOE时标，准确还原电网事故关键事件状态序列，为事故回放、分析和处理提供参考依据，为此提出一种基于录波文件基准通道的二次设备时标对齐方法。
技术实现思路
本专利技术所要解决的技术问题是提供一种基于录波文件基准通道的二次设备时标对齐方法，利用短时归一化相关系数模型，准确地统一不同设备间故障录波和SOE事件时标。为解决上述技术问题，本专利技术采用的技术方案是：一种基于录波文件基准通道的二次设备时标对齐方法，其特征在于，包括以下步骤：步骤1：考虑变电站二次设备录波数据各通道波形为能量有限的离散信号，而且每个录波文件内部通道时标一致。指定参考录波文件为A，选取参考录波信号，包括；步骤1.1：在A文件中选取某一录波通道，其信号可以表示为：其中，Dt为采样间隔，Dt=1/fs，fs是采样频率；n=0,1,2,…,N-1，N为采样点数；步骤1.2：将x[n]定义为参考录波信号，即取此信号采样时刻为基准时刻，求解其他待对齐录波信号的相对时间偏差。步骤2：设待对齐录波文件为B，在B文件中选取待对齐信号通道，定义待对齐信号为y[n]，可以表示为：步骤3：待对齐信号的选取原则包括：步骤3.1：录波文件B选取的通道信号y[n]与录波文件A选取的参考信号x[n]同为电压或电流信号；步骤3.2：信号采集于电力系统的同一电气节点处；步骤4：若两个信号采样频率不一致，分别为fs1和fs2，如果fs1≠fs2，则需要对高频采样信号进行重采样，重采样频率为min(fs1,fs2)；步骤5：提取待对齐信号y[n]的特征信号波段，包括：步骤5.1：在录波启动时刻附近△N范围内计算总谐波畸变率，表示为：其中，Gh为第h次谐波分量有效值，G1为基波分量有效值，计算窗长为采样点n的前一个周波采样点数，H为需要计算的最高谐波次数；步骤5.2：计算总谐波畸变率最大值，表示为：步骤5.3：采样点np处的总谐波畸变率最高，可计算出k为：步骤5.4：待对齐信号y[n]的特征信号波段可表示为：其中，k为特征信号区间的起始时刻，h[n]为矩形窗，所示为：其中，M为矩形窗窗长，也是特征信号长度；步骤6：定义x[n]和ys[n]的短时归一化相关系数可以表示为：其中，l为矩形窗移动的采样点数，且；步骤7：考虑到不同设备录波极性可能不一致，定义偏差时间识别指标，表示为：步骤8：计算偏差时间识别指标最大值位置，表示为：其中，lp即为计算出的偏差时间识别指标最大值位置；步骤9：得到两组信号时间之差为：步骤10：y[n]信号修正时标后表示为：步骤11：定义tSWx和tSWy分别为参考信号和待对齐信号所属二次设备录波波形及SOE事件的时间，计算不同二次设备录波波形及SOE时标对齐结果为：其中，为待对齐信号所属二次设备录波波形及SOE事件的校正后的时间。本专利技术的有益效果：采用本专利技术所述方法，对于不同装置的故障录波文件，提出了信号选取原则，利用总谐波畸变率分布特点提取待对齐信号的特征信号波段，建立了短时归一化相关系模型，实现了不同录波文件的时标校正，计算结果准确性高，不受CT极性接线方式的影响，且受噪声影响较小，有良好的抗干扰能力，提高了同一时间坐标下的故障时序分析的准确度。附图说明图1为本专利技术的实施例的参考信号和待对齐信号；图2为本专利技术的实施例提取待对齐信号的特征信号结果；图3为本专利技术的实施例偏差时间识别指标最大值位置计算结果。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下根据附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施仅仅用以解释本专利技术，并不限定本专利技术。下面结合附图和具体实施方式，对本专利技术做进一步的说明。步骤包括：步骤1：考虑变电站二次设备录波数据各通道波形为能量有限的离散信号，而且每个录波文件内部通道时标一致。指定参考录波文件为A，选取参考录波信号，包括；步骤1.1：在A文件中选取某一录波通道，其信号可以表示为：其中，Dt为采样间隔，Dt=1/fs，fs是采样频率；n=0,1,2,…,N-1，N为采样点数；步骤1.2：将x[n]定义为参考录波信号，即取此信号采样时刻为基准时刻，求解其他待对齐录波信号的相对时间偏差。步骤2：设待对齐录波文件为B，在B文件中选取待对齐信号通道，定义待对齐信号为y[n]，可以表示为：步骤3：待对齐信号的选取原则包括：步骤3.1：录波文件B选取的通道信号y[n]与录波文件A选取的参考信号x[n]同为电压或电流信号；步骤3.2：信号采集于电力系统的同一电气节点处；步骤4：若两个信号采样频率不一致，分别为fs1和fs2，如果fs1≠fs2，则需要对高频采样信号进行重采样，重采样频率为min(fs1,fs2)；步骤5：提取待对齐信号y[n]的特征信号波段，包括：步骤5.1：在录波启动时刻附近△N范围内计算总谐波畸变率，表示为：其中，Gh为第h次谐波分量有效值，G1为基波分量有效值，计算窗长为采样点n的前一个周波采样点数，H为需要计算的最高谐波次数；步骤5.2：计算总谐波畸变率最大值，表示为：步骤5.3：采样点np处的总谐波畸变率最高，可计算出k为：步骤5.4：待对齐信号y[n]的特征信号波段可表示为：其中，k为特征信号区间的起始时刻，h[n]为矩形窗，所示为：其中，M为矩形窗窗长，也是特征信号长度；步骤6：定义x[n]和ys[n]的短时归本文档来自技高网...

【技术保护点】
1.一种基于录波文件基准通道的二次设备时标对齐方法，其特征在于，包括以下步骤：/n步骤一：指定参考录波文件为A，选取参考录波信号，包括；/n步骤1.1：在A文件中选取某一录波通道，其信号表示为：/n

【技术特征摘要】
1.一种基于录波文件基准通道的二次设备时标对齐方法，其特征在于，包括以下步骤：
步骤一：指定参考录波文件为A，选取参考录波信号，包括；
步骤1.1：在A文件中选取某一录波通道，其信号表示为：



其中，Dt为采样间隔，Dt=1/fs，fs是采样频率；n=0,1,2,…,N-1，N为采样点数；
步骤1.2：将x[n]定义为参考录波信号，求解其他待对齐录波信号的相对时间偏差；
步骤二：设待对齐录波文件为B，在B文件中选取待对齐信号通道，定义待对齐信号为y[n]：

；
步骤三：选取待对齐的信号；
步骤四：若两个信号采样频率不一致，分别为fs1和fs2，fs1≠fs2，则需要对高频采样信号进行重采样，重采样频率为min(fs1,fs2)；
步骤五：提取待对齐信号y[n]的特征信号波段；
步骤六：定义x[n]和ys[n]的短时归一化相关系数：



其中，l为矩形窗移动的采样点数，且，h[n-l]为矩形窗h[n]右移l个采样点数，ys[n-l]为待对齐信号y[n]右移l个采样点数；
步骤七：设定偏差时间识别指标，表示为：

；
步骤八：计算偏差时间识别指标最大值位置：

；
其中，lp即为计算出的偏差时间识别指标最大值位置；
步骤九：得到两组信号时间之差为：

；
其中，k为特征信号区间的起始时刻；
步骤十：y[n]信号修正...

【专利技术属性】
技术研发人员：陈军，孙仲民，赵子涵，庞晓艳，朱小红，向博，文继锋，吕航，叶翔，
申请(专利权)人：南京南瑞继保电气有限公司，国网四川省电力公司，
类型：发明
国别省市：江苏;32