本发明涉及继电保护领域，具体涉及一种继电保护故障回放方法。
背景技术：
目前，我国电力系统故障录波发展迅猛，它可以自动准确地记录故障前后过程中各种电气量的变化情况，记录系统频率的全过程变化现象，以及记载暂态信息，通过对这些信息的记录对电力故障发生原因提供一定的参考价值，录波器是独立于保护之外的记录设备，通过记载模拟量和开关量来真实反应事故情况，是判断事故发展过程的最有效的工具。
故障录波在电力系统发生故障时，自动地、准确地记录电力系统故障前、后过程的各种电气量(主要数字量，比如开关状态变化，模拟量，主要是电压、电流数值)的变化情况，通过这些电气量的分析、比较，对分析处理事故、判断保护是否正确动作、提高电力系统安全运行水平的作用。
故障信息系统中含有海量的故障录波数据，其往往具有时序性、时钟不同步性、不一致性、不完整性、冗余性等特征，同一时标下的输电线路两端录波数据并不一定是匹配的。把输电线路两端双端故障录波数据匹配，将其应用于双端故障测距，保护行为分析，故障回放、事故状态下等值校验等，将能更好的发挥数据的价值，对故障分析与故障恢复具有重要的意义。
故障录波数据是高速采样的，对不同装置产生的录波数据进行综合分析便对时钟精度提出了很高的要求，对进行深层次的分析造成了障碍。
技术实现要素：
本发明提供一种继电保护故障回放方法，本发明采用数据压缩模块，基于波形周期性特性的前后波形差分压缩算法，具有低功耗、计算量小、压缩比例高的特点，有较好的波形数据压缩效果，较大幅度的降低继电保护故障录制以及回放系统的制造成本；本发明在故障录波数据的基础上，对电力系统一次设备进行建模，将一次设备模型与故障录波数据文件中的录波通道进行关联，得到录波分析模型，在故障发生时，能够根据电压、电流的突变特性，进行多录波数据的波形对齐；本发明分别从各故障录波文件中提取故障时刻前后一段时间的故障录波数据，以一次设备元件为单位，进行各一次设备元件的故障特征量的计算，形成中间数据，克服了故障录波文件本身的限制；本发明从多个角度分析录波数据，将录波采样信息用不同维度的特征量进行了数字化和量化处理，得到了多个维度的故障特征量，能够保证可靠性和快速性。
为实现上述目的，本发明的一种继电保护故障回放方法，所述方法包括以下步骤：
s1.实时获取电力设备波形，并转化成波形数据；
s2.压缩所述波形数据；
s3.分析所述波形数据，得到故障波形及故障数据；
s4.回放所述故障数据和波形。
优选的，所述步骤s2具体包括如下步骤：
s21.提取所压缩的电流、电压信号具备周期性的特征，允许包含干扰信号；所述的电流、电压信号在压缩过程中会产生前后周期的差值；
s22.将步骤s21的电流、电压信号在整个时间轴上的波形数据按照周期性封装为多个独立的数据单元，每个数据单元配置数据头；
s23.压缩过程采用多种子算法依次实验的方式，获得最优的压缩比例。
优选的，所述步骤s3具体包括如下步骤：
s31.基于录波数据，建立电力系统一次设备模型，将一次设备模型与故障录波文件中的录波通道进行关联，得到录波分析模型；
s32.提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件的故障特征量；
s33.将录波分析模型和故障特征量存入数据库，供故障诊断和事故分析使用。
优选的，所述步骤s31中，基于故障录波数据，建立电力系统一次设备模型，将一次设备模型与智能录波终端的录波通道进行关联，得到录波分析模型，包括：
s311.基于各故障的录波数据，建立电力系统一次设备模型；
s322.在一次设备模型中增加各录波数据中对应的多个模拟量通道号和多个开关量通道的通道号，形成录波分析模型。
优选的，所述录波分析模型包括母线模型、输电线路模型、变压器模型、断路器模型和保护装置模型；
其中，所述母线模型包括母线名称、母线编号、电压等级、相关的保护装置编号、相关的断路器设备编号和对应录波数据中的电压通道号；
所述输电线路模型包括线路名称、线路编号、电压等级、线路阻抗参数、线路长度、相关的保护装置编号、相关的断路器设备编号、对应录波数据中的电流通道号；
所述变压器模型包括变压器名称、变压器编号、二或三侧变压器卷、相关的保护装置编号；
所述断路器模型包括对应故障录波数据中开关量通道号和开关量所属信号类型；
所述保护装置模型包括对应故障录波数据中开关量通道号和开关量所属信号类型。
优选的，得到录波分析模型后，在故障发生时，通过遍历波形数据得到电流突变点，根据电流突变点进行多故障录波数据的波形对齐。
优选的，所述故障特征量包括故障前后电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量、断路器动作特征量和波形数据；一次设备元件包括母线、输电线路、变压器和断路器。
优选的，提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件故障前后电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量，包括：
从各故障录波文件中提取故障时刻前后一段时间的故障录波数据；
在故障前到故障后一段时间内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量；
在重合闸时刻到重合后时间范围内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量。
优选的，提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件故障前后波形数据，包括：
分别从各故障录波文件中提取母线对应的三相电压的波形数据、线路对应的三相电压和三相电流的波形数据、变压器各侧及中性点对应的三相电压和三相电流的波形数据以及断路器对应的三相电流的波形数据；
将提取到的波形数据进行统一归一化处理，并计算各一次设备元件故障前3周波波形数据、故障后10周波波形数据和重合闸后10周波波形数据一次值；
按照comtrade文件格式存储各一次设备元件故障前3周波波形数据、故障后10周波波形数据和重合闸后10周波波形数据一次值。
本发明的技术方案具有如下优点：
(1)本发明采用数据压缩模块，基于波形周期性特性的前后波形差分压缩算法，具有低功耗、计算量小、压缩比例高的特点，有较好的波形数据压缩效果，较大幅度的降低继电保护故障录制以及回放系统的制造成本。
(2)本发明在故障录波数据的基础上，对电力系统一次设备进行建模，将一次设备模型与故障录波数据文件中的录波通道进行关联，得到录波分析模型，在故障发生时，能够根据电压、电流的突变特性，进行多录波数据的波形对齐；本发明分别从各故障录波文件中提取故障时刻前后一段时间的故障录波数据，以一次设备元件为单位，进行各一次设备元件的故障特征量的计算，形成中间数据，克服了故障录波文件本身的限制。
(3)本发明从多个角度分析录波数据，将录波采样信息用不同维度的特征量进行了数字化和量化处理，得到了多个维度的故障特征量，能够保证可靠性和快速性。
附图说明
图1示出了本发明的一种继电保护故障回放系统的框图；
图2示出了本发明的一种继电保护故障回放方法。
具体实施方式
图1示出了本发明的一种继电保护故障回放系统，该系统包括智能录波终端1和录波分析及回放装置2；
该智能录波终端1包括：
波形获取模块11，用于实时获取电力设备的波形信息；
数据转换模块12，用于实时将波形信息转换成波形数据；
数据压缩模块13，用于将波形数据进行压缩；
录波数据发送模块15，用于通过无线通信网络将压缩后的波形数据发送给录波分析及回放装置；
终端控制模块14，并用于智能录波终端中的各模块工作；
所述录波分析及回放装置2包括：
录波数据接收模块21，用于接收智能终端发送过来的录波数据；
录波分析模块22，用于分析所述录波数据，得到故障波形及数据；
录波回放模块23，用于回放故障波形及数据；
中控模块24，用于协调控制整个继电保护故障回放系统。
数据压缩模块13包括：
周期性特征提取单元，用于提取所压缩的电流、电压信号具备周期性的特征，允许包含干扰信号，所述的电流、电压信号在压缩过程中会产生前后周期的差值；
数据封装单元，将所述电流、电压信号在整个时间轴上的波形数据按照周期特征封装为多个独立的数据单元，每个数据单元配置数据头。
优选的，数据压缩模块13在压缩过程采用多种子算法依次实验的方式，获得最优的压缩比例。
具体的，数据压缩模块13采用以下算法进行数据压缩：
对于故障指示器录取的波形，采用前后波形相减的方式，将后一个波形中的某一个数据点减去前一个波形的相应的同一个数据点，以获得更小的差分数据，此差分数据通过最高位是否为1来区分正负。
将时间轴上的波形数据按照其周期时间分为多个子数据单元，每个子数据单元先按照前后波形差分的方式，得到当前波形与前一个波形的差值，并将差值进行位宽压缩。
在进行前后波形差分的过程中，每个数据单元被简化为相比于前一个波形的差值表示；搜索波形差值的最大偏差量，记录最大偏差量位数作为基准，进行位宽压缩，其差分最大位数包括表示正负的符号位。
对于首个波形作为后续波形的差分参照波形，采用两种压缩算法，一种为传统的位宽压缩算法，一种为带符号的位宽压缩算法，两种算法各统计压缩完后的数据字节数，选取数据字节数最少的作为首个波形的压缩方式。
优选的，采用多个已存在的子算法的组合，通过对子算法的最优选择，达到最大压缩比例。
所述录波分析模块22包括：
一次设备模型建立单元，用于基于故障录波数据，建立一次设备模型；
录波分析模型建立单元，用于将一次设备模型与故障录波文件中的录波通道进行关联，得到录波分析模型；
数据提取单元，用于从录波数据中提取故障时刻前后一段时间的故障录波数据；
故障特征量计算单元，用于根据所提取的故障时刻前后一段时间内的故障录波数据，计算各一次设备元件的故障特征量；
存储单元，用于存储录波分析模型和故障特征量，供故障诊断和事故分析使用。
优选的，所述录波分析模块22还包括：
波形对齐单元，用于在故障发生时，通过遍历波形数据找到该突变点，进行多故障录波数据的波形对齐。
优选的，故障特征量计算单元，具体用于：
在故障前到故障后一段时间内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点的故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量；
接着在重合闸时刻到重合后时间范围内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点的故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量。
图2示出了本发明的一种继电保护故障回放方法，所述方法包括以下步骤：
s1.实时获取电力设备波形，并转化成波形数据；
s2.压缩所述波形数据；
s3.分析所述波形数据，得到故障波形及故障数据；
s4.回放所述故障数据和波形。
优选的，所述步骤s2具体包括如下步骤：
s21.提取所压缩的电流、电压信号具备周期性的特征，允许包含干扰信号；所述的电流、电压信号在压缩过程中会产生前后周期的差值；
s22.将步骤s21的电流、电压信号在整个时间轴上的波形数据按照周期性封装为多个独立的数据单元，每个数据单元配置数据头；
s23.压缩过程采用多种子算法依次实验的方式，获得最优的压缩比例。
具体的，数据压缩模块13采用以下算法进行数据压缩：
对于故障指示器录取的波形，采用前后波形相减的方式，将后一个波形中的某一个数据点减去前一个波形的相应的同一个数据点，以获得更小的差分数据，此差分数据通过最高位是否为1来区分正负。
将时间轴上的波形数据按照其周期时间分为多个子数据单元，每个子数据单元先按照前后波形差分的方式，得到当前波形与前一个波形的差值，并将差值进行位宽压缩。
在进行前后波形差分的过程中，每个数据单元被简化为相比于前一个波形的差值表示；搜索波形差值的最大偏差量，记录最大偏差量位数作为基准，进行位宽压缩，其差分最大位数包括表示正负的符号位。
对于首个波形作为后续波形的差分参照波形，采用两种压缩算法，一种为传统的位宽压缩算法，一种为带符号的位宽压缩算法，两种算法各统计压缩完后的数据字节数，选取数据字节数最少的作为首个波形的压缩方式。
优选的，采用多个已存在的子算法的组合，通过对子算法的最优选择，达到最大压缩比例。
优选的，对电流电压的采样精度按照12位ad采样，每组波形包含960个采样点，共12个波形，第一步是将12个波形分为12个数据单元，每个数据单元增加数据头，每个数据单元内容有压缩类型、压缩后单点数据位宽和所压缩数据采样点数；
第二步是对12个波形中的首个波形进行压缩：在压缩时选用两种压缩方式，位宽压缩与前后点差分压缩，在进行位宽压缩时，采用12bit的压缩方式，压缩比例为1/4；在进行前后采样点差分压缩时，将差分值最大值位宽加上符号位作为参照位宽；比较两种压缩方式，选用压缩后位宽最小的压缩算法作为首个波形的压缩方式。
第三步是取后续11个波形中的第n个波形，设波形中采样数据点表示为pn，i，n＝[2，3，..，12]，i＝[1，2，..，80]，则对于波形wn可表达为wn＝{pn，i}，i＝[1，2，..，80]；根据
其中n＝[2，3，..，11]，i＝[1，2，..，80]；
则全波形数据可做如下等价：
压缩后的波形数据单元位宽若小于等于12位，则执行n＝n+1，回到重复第三步执行；若压缩后的位宽大于12位，执行下述的第四步；
第四步，对于前后波形差分压缩算法非最优的情况，则使用前后采样点差分压缩算法，对于波形单元wn，n＝[2，3，..1，2]，pn可做如下等效
当压缩后的波形数据单元位宽小于等于12位时，则执行n＝n+1，重复第三步；当压缩后的位宽大于12位，则执行下述的第五步；
第五步，对于前后采样点差分压缩算法非最优的情况，则使用12bit位宽压缩方式，执行n＝n+1，回到重复第三步执行。
优选的，所述步骤s3具体包括如下步骤：
s31.基于录波数据，建立电力系统一次设备模型，将一次设备模型与故障录波文件中的录波通道进行关联，得到录波分析模型；
s32.提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件的故障特征量；
s33.将录波分析模型和故障特征量存入数据库，供故障诊断和事故分析使用。
优选的，所述步骤s31中，基于故障录波数据，建立电力系统一次设备模型，将一次设备模型与智能录波终端的录波通道进行关联，得到录波分析模型，包括：
s311.基于各故障的录波数据，建立电力系统一次设备模型；
s322.在一次设备模型中增加各录波数据中对应的多个模拟量通道号和多个开关量通道的通道号，形成录波分析模型。
优选的，所述录波分析模型包括母线模型、输电线路模型、变压器模型、断路器模型和保护装置模型；
其中，所述母线模型包括母线名称、母线编号、电压等级、相关的保护装置编号、相关的断路器设备编号和对应录波数据中的电压通道号；
所述输电线路模型包括线路名称、线路编号、电压等级、线路阻抗参数、线路长度、相关的保护装置编号、相关的断路器设备编号、对应录波数据中的电流通道号；
所述变压器模型包括变压器名称、变压器编号、二或三侧变压器卷、相关的保护装置编号；
所述断路器模型包括对应故障录波数据中开关量通道号和开关量所属信号类型；
所述保护装置模型包括对应故障录波数据中开关量通道号和开关量所属信号类型。
优选的，得到录波分析模型后，在故障发生时，通过遍历波形数据得到电流突变点，根据电流突变点进行多故障录波数据的波形对齐。
优选的，所述故障特征量包括故障前后电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量、断路器动作特征量和波形数据；一次设备元件包括母线、输电线路、变压器和断路器。
优选的，提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件故障前后电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量，包括：
从各故障录波文件中提取故障时刻前后一段时间的故障录波数据；
在故障前到故障后一段时间内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量；
在重合闸时刻到重合后时间范围内，每间隔10ms取一个点，根据该点的故障录波数据，计算该点故障前后各一次设备对应的电压电流有效值、5-15次谐波值、直流分量及时间衰减常数、序分量、差流值、线路测量阻抗、频率特征量、故障测距特征量、保护动作特征量和断路器动作特征量。
优选的，提取故障时刻前后一段时间内的故障录波数据，计算各一次设备元件故障前后波形数据，包括：
分别从各故障录波文件中提取母线对应的三相电压的波形数据、线路对应的三相电压和三相电流的波形数据、变压器各侧及中性点对应的三相电压和三相电流的波形数据以及断路器对应的三相电流的波形数据；
将提取到的波形数据进行统一归一化处理，并计算各一次设备元件故障前3周波波形数据、故障后10周波波形数据和重合闸后10周波波形数据一次值；
按照comtrade文件格式存储各一次设备元件故障前3周波波形数据、故障后10周波波形数据和重合闸后10周波波形数据一次值。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，做出若干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明的保护范围。