Nothing Special   »   [go: up one dir, main page]

CN109828183B - Waveform synchronous correction method and system suitable for transient recording type fault indicator - Google Patents

Waveform synchronous correction method and system suitable for transient recording type fault indicator Download PDF

Info

Publication number
CN109828183B
CN109828183B CN201811475046.3A CN201811475046A CN109828183B CN 109828183 B CN109828183 B CN 109828183B CN 201811475046 A CN201811475046 A CN 201811475046A CN 109828183 B CN109828183 B CN 109828183B
Authority
CN
China
Prior art keywords
waveform
time
modal
take
transient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811475046.3A
Other languages
Chinese (zh)
Other versions
CN109828183A (en
Inventor
吴楠
史明明
袁宇波
孙健
杨雄
焦昊
祁承磊
周涛
张刘冬
方鑫
曾飞
李娟�
朱卫平
何明
袁晓冬
陈兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811475046.3A priority Critical patent/CN109828183B/en
Publication of CN109828183A publication Critical patent/CN109828183A/en
Application granted granted Critical
Publication of CN109828183B publication Critical patent/CN109828183B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Relating To Insulation (AREA)
  • Locating Faults (AREA)

Abstract

本发明公开了适用于暂态录波型故障指示器的波形同步校正方法及其系统,方法包括获取采集单元的三相电流录波波形;对a相电流波形的最后f/50*2个采样点进行模态分析;读取a相电流波形的触发录波时刻,向前取f/50个采样点,对应时刻t1,t1下一个采样点对应时刻为t2,以此类推,则由tg向后取f/50个采样点;进行模态分析获得误差能量△E1、△E2……△Eg,取其中最小值对应的时间窗起始时刻点tfa;依次完成b、c相分析,得到三相的故障起始时刻点,取其中最小值;将三相电流分别向左平移,得到三相同步录波波形。本发明对暂态录波型故障指示器的录波波形进行校正,保证用于主站分析的录波波形具有高准确性。

Figure 201811475046

The invention discloses a waveform synchronization correction method and a system suitable for a transient recording type fault indicator. The method includes acquiring a three-phase current recording waveform of an acquisition unit; sampling the last f/50*2 of the a-phase current waveform The modal analysis is carried out at the point of t1; read the trigger recording time of the a-phase current waveform, take f/50 sampling points forward, corresponding to time t1, the corresponding time of the next sampling point of t1 is t2, and so on, from tg to Then take f/50 sampling points; perform modal analysis to obtain error energy △E1, △E2...△Eg, and take the time window starting time point tfa corresponding to the minimum value; complete the b and c phase analysis in turn, and obtain three The fault starting time point of the phase, take the minimum value; shift the three-phase currents to the left respectively to obtain the three-phase synchronous recording waveform. The invention corrects the wave recording waveform of the transient wave recording type fault indicator, so as to ensure the high accuracy of the wave recording waveform used for the analysis of the main station.

Figure 201811475046

Description

适用于暂态录波型故障指示器的波形同步校正方法及系统Waveform synchronization correction method and system suitable for transient waveform recording type fault indicator

技术领域technical field

本发明涉及一种适用于暂态录波型故障指示器的波形同步校正方法及系统,属于电力系统配电自动化技术领域。The invention relates to a waveform synchronization correction method and system suitable for a transient recording type fault indicator, belonging to the technical field of power system distribution automation.

背景技术Background technique

目前,电网正积极推进配电自动化建设,着力提升配电自动化应用水平。暂态录波型故障指示器是安装于配电线路上,用于监测线路场强、电流,并在故障时刻进行录波的装置。暂态录波型故障指示器可将告警信号、录波文件传给主站,用于主站故障研判及分析决策。At present, the power grid is actively promoting the construction of distribution automation, focusing on improving the application level of distribution automation. The transient wave recording type fault indicator is a device installed on the distribution line to monitor the field strength and current of the line, and to record the wave at the time of the fault. The transient wave recording type fault indicator can transmit the alarm signal and wave recording file to the main station, which is used for the main station fault judgment and analysis and decision-making.

暂态录波型故障指示器由三个采集单元、一个汇集单元组成,三个采集单元分别装在 a、b、c三相线路上,汇集单元装在电线杆上,采集单元与汇集单元间通过小无线通信,汇集单元与主站间通过无线公网通信。暂态录波型故障指示器的一个难点是,a、b、c三相录波存在不同步的问题,导致零序分量合成时误差较大,影响主站波形分析结果。The transient recording type fault indicator is composed of three acquisition units and one collection unit. The three acquisition units are respectively installed on the a, b, and c three-phase lines. The collection unit is installed on the utility pole. The connection between the collection unit and the collection unit is Through small wireless communication, the collection unit communicates with the master station through the wireless public network. One of the difficulties of the transient recording type fault indicator is that the a, b, and c three-phase recordings are out of synchronization, which leads to a large error in the synthesis of the zero-sequence components, which affects the waveform analysis results of the master station.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的技术问题是克服现有技术的缺陷,提出一种适用于暂态录波型故障指示器的波形同步校正方法及系统,对暂态录波型故障指示器的录波波形进行校正,保证用于主站分析的录波波形具有高准确性。The technical problem to be solved by the present invention is to overcome the defects of the prior art, and propose a waveform synchronization correction method and system suitable for the transient recording type fault indicator. Calibration to ensure high accuracy of recorded waveforms for master station analysis.

为解决上述技术问题,本发明采用的技术方案如下:In order to solve the above-mentioned technical problems, the technical scheme adopted in the present invention is as follows:

本发明的适用于暂态录波型故障指示器的波形同步校正方法,包括以下步骤:The waveform synchronization correction method applicable to the transient recording type fault indicator of the present invention comprises the following steps:

步骤1、故障指示器汇集单元获取采集单元的a、b、c三相电流录波波形;Step 1. The fault indicator collection unit acquires the a, b, and c three-phase current recording waveforms of the acquisition unit;

步骤2、对a相电流波形的最后f/50*2个采样点进行模态分析,其中,f为采集单元采样频率;读取模态中幅值最大的模态参数:幅值A0、频率f0;Step 2. Perform modal analysis on the last f/50*2 sampling points of the a-phase current waveform, where f is the sampling frequency of the acquisition unit; read the modal parameters with the largest amplitude in the modal: amplitude A0, frequency f0;

步骤3、读取a相电流波形的触发录波时刻tg,由tg向前取f/50个采样点,对应时刻t1,t1下一个采样点对应时刻为t2,以此类推,则由tg向后取f/50个采样点,对应时刻 t(1+f/50*2);Step 3. Read the trigger recording time tg of the current waveform of phase a, take f/50 sampling points forward from tg, corresponding to time t1, the corresponding time of the next sampling point of t1 is t2, and so on, from tg to Then take f/50 sampling points, corresponding to time t(1+f/50*2);

步骤4、对[t1,t(f/50*2)]时间区间共f/50*2个采样点进行模态分析,取n个主模态参数:A11、f11;A12、f12……A1n、f1n,分别计算n个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12……ΔE1n,取ΔE11、ΔE12……ΔE1n中的最小值,记为ΔE1;Step 4. Perform modal analysis on a total of f/50*2 sampling points in the [t1, t(f/50*2)] time interval, and take n main modal parameters: A11, f11; A12, f12...A1n , f1n, calculate the error energy ΔE11, ΔE12...ΔE1n of the n main modes and the modes formed by A0, f0 respectively, take the minimum value of ΔE11, ΔE12...ΔE1n, and denote it as ΔE1;

步骤5、将[t1,t(f/50*2)]的时间窗向后平移一个采样点,即对[t2,t(f/50*2+1)]时间区间共f/50*2个采样点进行模态分析,同步骤4,取最小误差能量,记为ΔE2;Step 5. Move the time window of [t1, t(f/50*2)] backward by one sampling point, that is, for the time interval of [t2, t(f/50*2+1)], a total of f/50*2 Perform modal analysis at each sampling point, same as step 4, take the minimum error energy, denoted as ΔE2;

步骤6、依次对[tx,t(x+f/50*2-1)]时间区间共f/50*2个采样点进行模态分析,同步骤 4,取最小误差能量,记为ΔEx,直至完成对tg为起始点的f/50*2个采样点的模态分析;Step 6. Perform modal analysis on a total of f/50*2 sampling points in the [tx, t(x+f/50*2-1)] time interval. Same as step 4, take the minimum error energy, denoted as ΔEx, Until the modal analysis of f/50*2 sampling points with tg as the starting point is completed;

步骤7、获得误差能量ΔE1、ΔE2、…、ΔEx、…、ΔEg,取其中最小值对应的时间窗起始时刻点tfa;Step 7: Obtain the error energies ΔE1, ΔE2, ..., ΔEx, ..., ΔEg, and take the start time point tfa of the time window corresponding to the minimum value;

步骤8、依次完成b相、c相分析,得到三相的故障起始时刻点tfb、tfc,取tfa、tfb、tfc中的最小值,记为故障发生的三相同步时刻点tf;将a相、b相、c相电流分别向左平移tfa-tf、tfb-tf、tfc-tf,得到三相同步录波波形。Step 8. Complete the b-phase and c-phase analysis in turn, and obtain the three-phase fault starting time points tfb and tfc, take the minimum value of tfa, tfb, and tfc, and record it as the three-phase synchronization time point tf when the fault occurs; The phase, b-phase, and c-phase currents are shifted to the left by tfa-tf, tfb-tf, and tfc-tf, respectively, to obtain three-phase synchronous waveforms.

f/50为小数时,向上取整。When f/50 is a decimal, round up.

步骤2中,所述模态分析的方法如下:In step 2, the modal analysis method is as follows:

对波形进行模态分解,将一波形表示成多个模态之和,每个模态表示形式为 Asin(2*π*fx*t+a),其中,A为模态的幅值,fx为模态的频率,a为模态的初相角;The waveform is modal decomposed, and a waveform is expressed as the sum of multiple modes, each mode is represented in the form of Asin(2*π*fx*t+a), where A is the amplitude of the mode, fx is the frequency of the mode, and a is the initial phase angle of the mode;

以下仅对幅值A、频率fx进行应用,即模态可简化表示为Asin(2*π*fx*t),读取模态中幅值最大的模态参数,记为幅值A0、频率f0。The following only applies to the amplitude A and frequency fx, that is, the modal can be simplified as Asin(2*π*fx*t), read the modal parameter with the largest amplitude in the modal, and record it as amplitude A0, frequency f0.

步骤3中,所述触发录波时刻tg位于故障发生后半个工频周波处。In step 3, the trigger wave recording time tg is located at half the power frequency cycle after the fault occurs.

步骤3中,所述时刻t1位于tg前一个工频周波处,且故障时刻落在[t1,tg]间。In step 3, the time t1 is located at the power frequency cycle before tg, and the fault time is between [t1, tg].

步骤4中,取三个主模态参数:A11、f11;A12、f12;A13、f13。In step 4, three main modal parameters are taken: A11, f11; A12, f12; A13, f13.

三个主模态的定义是:三个模态分量相加所得波形W1,与[t1,t(f/50*2)]区间的波形W0误差能量△E最小,两个波形误差能量△E的定义为波形差值的绝对值在时间区间上的积分,此处W1与W0的误差能量表示为:

Figure GDA0002899693510000021
The definition of the three main modes is: the waveform W1 obtained by adding the three modal components, and the waveform W0 in the interval [t1, t(f/50*2)], the error energy △E is the smallest, and the error energy △E of the two waveforms is the smallest. is defined as the integral of the absolute value of the waveform difference over the time interval, where the error energy of W1 and W0 is expressed as:
Figure GDA0002899693510000021

分别计算三个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12、ΔE13,则ΔE11、ΔE12、ΔE13分别表示为:Calculate the error energy ΔE11, ΔE12, ΔE13 of the three main modes and the modes formed by A0, f0, then ΔE11, ΔE12, ΔE13 are expressed as:

Figure GDA0002899693510000022
Figure GDA0002899693510000022

Figure GDA0002899693510000023
Figure GDA0002899693510000023

Figure GDA0002899693510000024
Figure GDA0002899693510000024

取ΔE11、ΔE12、ΔE13中的最小值,记为ΔE1。Take the minimum value of ΔE11, ΔE12, and ΔE13, and denote it as ΔE1.

步骤6中,将f/50*2个采样点的波形分析窗依次向后平移,每次平移一个采样点,直至完成对tg为起始点的f/50*2个采样点的模态分析。In step 6, the waveform analysis window of f/50*2 sampling points is shifted backward in turn, one sampling point at a time, until the modal analysis of f/50*2 sampling points with tg as the starting point is completed.

本发明的适用于暂态录波型故障指示器的波形同步校正系统,所述系统包括网络接口、存储器和处理器;其中,The waveform synchronization correction system applicable to the transient recording type fault indicator of the present invention comprises a network interface, a memory and a processor; wherein,

所述网络接口,用于在与其他外部网元之间进行收发信息过程中,信号的接收和发送;The network interface is used for receiving and sending signals in the process of sending and receiving information with other external network elements;

所述存储器,用于存储能够在所述处理器上运行的计算机程序指令;the memory for storing computer program instructions executable on the processor;

所述处理器,用于在运行所述计算机程序指令时,执行上述适用于暂态录波型故障指示器的波形同步校正方法的步骤。The processor is configured to execute the steps of the above-mentioned waveform synchronization correction method applicable to the transient waveform recording type fault indicator when the computer program instructions are executed.

与现有技术相比,本发明的有益效果在于,本发明的方法对暂态录波型故障指示器的录波波形进行校正,保证了用于主站分析的录波波形具有高准确性,对配电自动化系统的实用化应用具有重要意义。Compared with the prior art, the beneficial effect of the present invention is that the method of the present invention corrects the recorded waveform of the transient recording type fault indicator, which ensures that the recorded waveform used for the main station analysis has high accuracy, It is of great significance to the practical application of distribution automation system.

附图说明Description of drawings

图1是本发明的适用于暂态录波型故障指示器的波形同步校正方法工作流程图。FIG. 1 is a working flow chart of the waveform synchronization correction method applicable to the transient waveform recording type fault indicator of the present invention.

具体实施方式Detailed ways

以下将结合附图对本发明各实施例的技术方案进行清楚、完整的描述,显然,所描述发实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例,都属于本发明所保护的范围。The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present invention.

本发明提供的适用于暂态录波型故障指示器的波形同步校正方法,利用暂态录波型故障指示器的汇集单元对采集单元的录波波形进行校正,如图1所示。具体过程如下:The waveform synchronization correction method applicable to the transient waveform recording type fault indicator provided by the present invention utilizes the collection unit of the transient waveform recording type fault indicator to correct the waveform recording waveform of the acquisition unit, as shown in FIG. 1 . The specific process is as follows:

步骤1、故障指示器汇集单元获取采集单元的a、b、c三相电流录波波形;Step 1. The fault indicator collection unit acquires the a, b, and c three-phase current recording waveforms of the acquisition unit;

步骤2、对a相电流波形的最后f/50*2个采样点进行模态分析,其中f为采集单元采样频率,f/50为小数时,向上取整,以下f/50同理,读取模态中幅值最大的模态参数,即幅值A0、频率f0,模态分析法是指对波形进行模态分解,将一波形表示成多个模态之和,每个模态表示形式为Asin(2*π*fx*t+a),其中,A为模态的幅值,fx为模态的频率,a为模态的初相角,以下仅对幅值A、频率fx进行应用,即模态可简化表示为Asin(2*π*fx*t)。读取模态中幅值最大的模态参数,记为幅值A0、频率f0,该模态理论上接近工频50Hz。;Step 2. Perform modal analysis on the last f/50*2 sampling points of the a-phase current waveform, where f is the sampling frequency of the acquisition unit, and when f/50 is a decimal, round up, the following f/50 is the same, read Take the modal parameter with the largest amplitude in the modal, that is, the amplitude A0 and the frequency f0. The modal analysis method refers to the modal decomposition of the waveform, and a waveform is expressed as the sum of multiple modes. Each mode represents The form is Asin(2*π*fx*t+a), where A is the amplitude of the modal, fx is the frequency of the modal, and a is the initial phase angle of the modal. The following is only for the amplitude A and frequency fx To apply, the mode can be simplified as Asin(2*π*fx*t). Read the modal parameter with the largest amplitude in the modal, and record it as the amplitude A0 and the frequency f0. This modal is theoretically close to the power frequency of 50Hz. ;

步骤3、读取a相电流波形的触发录波时刻tg,由tg向前取f/50个采样点,对应时刻t1,t1下一个采样点对应时刻为t2,以此类推,则由tg向后取f/50个采样点,对应时刻 t(1+f/50*2),触发录波时刻tg应位于故障发生后半个工频周波处,由tg向前取f/50个采样点,对应时刻为t1,故t1位于tg前一个工频周波处,且故障时刻应落在[t1,tg]间。t1下一个采样点对应时刻为t2,以此类推,则由tg向后取f/50个采样点,对应时刻t(1+f/50*2);Step 3. Read the trigger recording time tg of the current waveform of phase a, take f/50 sampling points forward from tg, corresponding to time t1, the corresponding time of the next sampling point of t1 is t2, and so on, from tg to Take f/50 sampling points later, corresponding to time t(1+f/50*2), the trigger recording time tg should be located at half the power frequency cycle after the fault occurs, and take f/50 sampling points forward from tg , the corresponding time is t1, so t1 is located at the power frequency cycle before tg, and the fault time should fall between [t1, tg]. The time corresponding to the next sampling point of t1 is t2, and so on, then take f/50 sampling points backward from tg, corresponding to time t(1+f/50*2);

步骤4、对[t1,t(f/50*2)]时间区间共f/50*2个采样点进行模态分析,取三个主模态参数:A11、f11;A12、f12;A13、f13,分别计算三个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12、ΔE13,取ΔE11、ΔE12、ΔE13中的最小值,记为ΔE1,三个主模态的定义是:三个模态分量相加所得波形W1,与[t1,t(f/50*2)]区间的波形W0误差能量△E 最小,两个波形误差能量△E的定义为波形差值的绝对值在时间区间上的积分,此处W1 与W0的误差能量可示意表示为:

Figure GDA0002899693510000041
取三个主模态参数:A11、f11; A12、f12;A13、f13,分别计算三个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12、ΔE13,则ΔE11、ΔE12、ΔE13可分别表示为:Step 4. Perform modal analysis on a total of f/50*2 sampling points in the [t1, t(f/50*2)] time interval, and take three main modal parameters: A11, f11; A12, f12; A13, f13, calculate the error energy ΔE11, ΔE12, ΔE13 of the three main modes and the modes formed by A0 and f0 respectively, take the minimum value of ΔE11, ΔE12, and ΔE13, and denote it as ΔE1. The definition of the three main modes is: The waveform W1 obtained by adding the three modal components, and the waveform W0 in the interval [t1, t(f/50*2)] have the smallest error energy △E, and the error energy △E of the two waveforms is defined as the absolute value of the waveform difference The integration over the time interval, where the error energy of W1 and W0 can be represented schematically as:
Figure GDA0002899693510000041
Take three main mode parameters: A11, f11; A12, f12; A13, f13, calculate the error energy ΔE11, ΔE12, ΔE13 of the three main modes and the modes formed by A0, f0, then ΔE11, ΔE12, ΔE13 can be expressed as:

Figure GDA0002899693510000042
Figure GDA0002899693510000042

Figure GDA0002899693510000043
Figure GDA0002899693510000043

Figure GDA0002899693510000044
Figure GDA0002899693510000044

取ΔE11、ΔE12、ΔE13中的最小值,记为ΔE1。Take the minimum value of ΔE11, ΔE12, and ΔE13, and denote it as ΔE1.

步骤5、将[t1,t(f/50*2)]的时间窗向后平移一个采样点,即对[t2,t(f/50*2+1)]时间区间共f/50*2个采样点进行模态分析,同上,取最小误差能量,记为ΔE2;Step 5. Move the time window of [t1, t(f/50*2)] backward by one sampling point, that is, for the time interval of [t2, t(f/50*2+1)], a total of f/50*2 1 sampling point for modal analysis, same as above, take the minimum error energy, denoted as ΔE2;

步骤6、依次进行,将f/50*2个采样点的波形分析窗依次向后平移,每次平移一个采样点,对[tx,t(x+f/50*2-1)]时间区间共f/50*2个采样点进行模态分析,同上,取最小误差能量,记为ΔEx,直至完成对tg为起始点的f/50*2个采样点进行模态分析;Step 6. Carry out sequentially, move the waveform analysis window of f/50*2 sampling points backward in turn, and move one sampling point at a time, for the [tx, t(x+f/50*2-1)] time interval A total of f/50*2 sampling points for modal analysis, the same as above, take the minimum error energy and record it as ΔEx, until the modal analysis of f/50*2 sampling points with tg as the starting point is completed;

步骤7、获得误差能量ΔE1、ΔE2、…、ΔEx、…、ΔEg,取其中的最小值对应的时间窗起始时刻点tfa;Step 7: Obtain the error energies ΔE1, ΔE2, . . . , ΔEx, .

步骤8、依次完成b相、c相分析,得到三相的故障起始时刻点tfb、tfc,取tfa、tfb、tfc中的最小值,记为故障发生的三相同步时刻点tf,将a相、b相、c相电流分别向左平移tfa-tf、tfb-tf、tfc-tf,得到三相同步录波波形。Step 8. Complete the b-phase and c-phase analysis in turn, and obtain the three-phase fault starting time points tfb and tfc, take the minimum value of tfa, tfb, and tfc, and record it as the three-phase synchronization time point tf when the fault occurs. The phase, b-phase, and c-phase currents are shifted to the left by tfa-tf, tfb-tf, and tfc-tf, respectively, to obtain three-phase synchronous waveforms.

本发明的适用于暂态录波型故障指示器的波形同步校正系统包括网络接口、存储器和处理器;网络接口用于在与其他外部网元之间进行收发信息过程中,信号的接收和发送;存储器用于存储能够在所述处理器上运行的计算机程序指令;处理器,用于在运行计算机程序指令时,执行上述适用于暂态录波型故障指示器的波形同步校正方法的步骤。The waveform synchronization correction system suitable for the transient waveform recording type fault indicator of the present invention comprises a network interface, a memory and a processor; the network interface is used for receiving and sending signals in the process of sending and receiving information with other external network elements. the memory is used to store computer program instructions that can be run on the processor; the processor is used to execute the steps of the above-mentioned waveform synchronization correction method applicable to the transient waveform recording type fault indicator when the computer program instructions are executed.

以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。The basic principles and main features of the present invention and the advantages of the present invention have been shown and described above. Those skilled in the art should understand that the present invention is not limited by the above-mentioned embodiments, and the descriptions in the above-mentioned embodiments and the description are only to illustrate the principle of the present invention. Without departing from the spirit and scope of the present invention, the present invention will have Various changes and modifications fall within the scope of the claimed invention. The claimed scope of the present invention is defined by the appended claims and their equivalents.

Claims (10)

1.适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,包括以下步骤:1. The waveform synchronization correction method applicable to the transient recording waveform fault indicator is characterized in that, comprising the following steps: 步骤1、获取采集单元的三相电流波形;Step 1. Obtain the three-phase current waveform of the acquisition unit; 步骤2、对任意一相电流波形的最后f/50*2个采样点进行模态分析;f为采集单元采样频率;读取模态中幅值最大的模态参数:幅值A0、频率f0;Step 2. Perform modal analysis on the last f/50*2 sampling points of any phase current waveform; f is the sampling frequency of the acquisition unit; read the modal parameters with the largest amplitude in the modal: amplitude A0, frequency f0 ; 步骤3、读取步骤2选中一相电流波形的触发录波时刻tg,由tg向前取f/50个采样点,对应时刻t1,t1下一个采样点对应时刻为t2,以此类推,则由tg向后取f/50个采样点,对应时刻t(1+f/50*2);Step 3. Read Step 2. Select the trigger recording time tg of the current waveform of one phase, take f/50 sampling points forward from tg, corresponding to time t1, the corresponding time of the next sampling point of t1 is t2, and so on, then Take f/50 sampling points backward from tg, corresponding to time t(1+f/50*2); 步骤4、对[t1,t(f/50*2)]时间区间共f/50*2个采样点进行模态分析,取n个主模态参数:A11、f11、A12、f12……A1n、f1n,分别计算n个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12……ΔE1n,取ΔE11、ΔE12……ΔE1n中的最小值,记为ΔE1;Step 4. Perform modal analysis on a total of f/50*2 sampling points in the [t1, t(f/50*2)] time interval, and take n main modal parameters: A11, f11, A12, f12...A1n , f1n, calculate the error energy ΔE11, ΔE12...ΔE1n of the n main modes and the modes formed by A0, f0 respectively, take the minimum value of ΔE11, ΔE12...ΔE1n, and denote it as ΔE1; 步骤5、将[t1,t(f/50*2)]的时间窗向后平移一个采样点,即对[t2,t(1+f/50*2)]时间区间共f/50*2个采样点进行模态分析,同步骤4,取最小误差能量,记为ΔE2;Step 5. Move the time window of [t1, t(f/50*2)] backward by one sampling point, that is, for the time interval of [t2, t(1+f/50*2)], a total of f/50*2 Perform modal analysis at each sampling point, same as step 4, take the minimum error energy, denoted as ΔE2; 步骤6、依次对[tx,t(x+f/50*2-1)]时间区间共f/50*2个采样点进行模态分析,同步骤4,取最小误差能量,记为ΔEx,直至完成对tg为起始点的f/50*2个采样点的模态分析;Step 6. Perform modal analysis on a total of f/50*2 sampling points in the [tx, t(x+f/50*2-1)] time interval. Same as step 4, take the minimum error energy, denoted as ΔEx, Until the modal analysis of f/50*2 sampling points with tg as the starting point is completed; 步骤7、获得误差能量ΔE1、ΔE2、…、ΔEx、…、ΔEg,取其中最小值对应的时间窗起始时刻点tfa;Step 7: Obtain the error energies ΔE1, ΔE2, ..., ΔEx, ..., ΔEg, and take the start time point tfa of the time window corresponding to the minimum value; 步骤8、依次完成另外两相分析,得到三相的故障起始时刻点tfb、tfc,取tfa、tfb、tfc中的最小值,记为故障发生的三相同步时刻点tf;将三相电流分别向左平移tfa-tf、tfb-tf、tfc-tf,得到三相同步录波波形。Step 8. Complete the other two-phase analysis in turn to obtain the three-phase fault starting time points tfb and tfc, take the minimum value among tfa, tfb, and tfc, and record it as the three-phase synchronization time point tf of the fault; Move tfa-tf, tfb-tf, and tfc-tf to the left, respectively, to obtain three-phase synchronous recording waveforms. 2.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,f/50为小数时,向上取整。2 . The waveform synchronization correction method applicable to the transient waveform recording type fault indicator according to claim 1 , wherein, when f/50 is a decimal, it is rounded up. 3 . 3.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,步骤2中,所述模态分析的方法如下:3. The waveform synchronization correction method applicable to the transient recording type fault indicator according to claim 1, is characterized in that, in step 2, the method of described modal analysis is as follows: 对波形进行模态分解,将一波形表示成多个模态之和,每个模态表示形式为Asin(2*π*fx*t+a),其中,A为模态的幅值,fx为模态的频率,a为模态的初相角;The waveform is modal decomposed, and a waveform is expressed as the sum of multiple modes, each mode is represented in the form of Asin(2*π*fx*t+a), where A is the amplitude of the mode, fx is the frequency of the mode, and a is the initial phase angle of the mode; 以下仅对幅值A、频率fx进行应用,即模态可简化表示为Asin(2*π*fx*t),读取模态中幅值最大的模态参数,记为幅值A0、频率f0。The following only applies to the amplitude A and frequency fx, that is, the modal can be simplified as Asin(2*π*fx*t), read the modal parameter with the largest amplitude in the modal, and record it as amplitude A0, frequency f0. 4.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,步骤3中,所述触发录波时刻tg位于故障发生后半个工频周波处。4. The waveform synchronization correction method suitable for transient wave recording type fault indicator according to claim 1, wherein in step 3, the trigger wave recording time tg is located at the half power frequency cycle after the fault occurs . 5.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,步骤3中,所述时刻t1位于tg前一个工频周波处,且故障时刻落在[t1,tg]间。5 . The waveform synchronization correction method suitable for a transient recording type fault indicator according to claim 1 , wherein, in step 3, the time t1 is located at a power frequency cycle before tg, and the fault time falls. 6 . between [t1, tg]. 6.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,步骤4中,取三个主模态参数:A11、f11;A12、f12;A13、f13。6 . The waveform synchronization correction method suitable for transient wave recording type fault indicator according to claim 1 , wherein, in step 4, three main modal parameters are taken: A11, f11; A12, f12; A13 , f13. 7.根据权利要求6所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,三个主模态的定义是:三个模态分量相加所得波形W1,与[t1,t(f/50*2)]区间的波形W0误差能量△E最小,两个波形误差能量△E的定义为波形差值的绝对值在时间区间上的积分,此处W1与W0的误差能量表示为:
Figure FDA0002899693500000021
7. The waveform synchronization correction method applicable to the transient recording type fault indicator according to claim 6, wherein the definition of the three main modes is: the waveform W1 obtained by adding the three modal components, and the The waveform W0 error energy △E in the interval [t1, t(f/50*2)] is the smallest, and the error energy △E of the two waveforms is defined as the integral of the absolute value of the waveform difference over the time interval, where W1 and W0 The error energy of is expressed as:
Figure FDA0002899693500000021
8.根据权利要求7所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,分别计算三个主模态与A0、f0所构成模态的误差能量ΔE11、ΔE12、ΔE13,则ΔE11、ΔE12、ΔE13分别表示为:8 . The waveform synchronization correction method suitable for transient waveform recording type fault indicator according to claim 7 , wherein the error energies ΔE11 and ΔE12 of the three main modes and the modes formed by A0 and f0 are calculated respectively. 9 . , ΔE13, then ΔE11, ΔE12, ΔE13 are expressed as:
Figure FDA0002899693500000022
Figure FDA0002899693500000022
Figure FDA0002899693500000023
Figure FDA0002899693500000023
Figure FDA0002899693500000024
Figure FDA0002899693500000024
取ΔE11、ΔE12、ΔE13中的最小值,记为ΔE1。Take the minimum value of ΔE11, ΔE12, and ΔE13, and denote it as ΔE1.
9.根据权利要求1所述的适用于暂态录波型故障指示器的波形同步校正方法,其特征在于,步骤6中,将f/50*2个采样点的波形分析窗依次向后平移,每次平移一个采样点,直至完成对tg为起始点的f/50*2个采样点的模态分析。9. The waveform synchronization correction method applicable to the transient recording type fault indicator according to claim 1, wherein in step 6, the waveform analysis window of f/50*2 sampling points is shifted backwards in turn , and translate one sampling point at a time until the modal analysis of f/50*2 sampling points with tg as the starting point is completed. 10.适用于暂态录波型故障指示器的波形同步校正系统,其特征在于,所述系统包括网络接口、存储器和处理器;其中,10. A waveform synchronization correction system suitable for a transient waveform recording type fault indicator, characterized in that the system comprises a network interface, a memory and a processor; wherein, 所述网络接口,用于在与其他外部网元之间进行收发信息过程中,信号的接收和发送;The network interface is used for receiving and sending signals in the process of sending and receiving information with other external network elements; 所述存储器,用于存储能够在所述处理器上运行的计算机程序指令;the memory for storing computer program instructions executable on the processor; 所述处理器,用于在运行所述计算机程序指令时,执行权利要求1至9任一项所述适用于暂态录波型故障指示器的波形同步校正方法的步骤。The processor is configured to execute the steps of the waveform synchronization correction method applicable to the transient waveform recording type fault indicator according to any one of claims 1 to 9 when the computer program instructions are executed.
CN201811475046.3A 2018-12-04 2018-12-04 Waveform synchronous correction method and system suitable for transient recording type fault indicator Active CN109828183B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811475046.3A CN109828183B (en) 2018-12-04 2018-12-04 Waveform synchronous correction method and system suitable for transient recording type fault indicator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811475046.3A CN109828183B (en) 2018-12-04 2018-12-04 Waveform synchronous correction method and system suitable for transient recording type fault indicator

Publications (2)

Publication Number Publication Date
CN109828183A CN109828183A (en) 2019-05-31
CN109828183B true CN109828183B (en) 2021-06-08

Family

ID=66859806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811475046.3A Active CN109828183B (en) 2018-12-04 2018-12-04 Waveform synchronous correction method and system suitable for transient recording type fault indicator

Country Status (1)

Country Link
CN (1) CN109828183B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111220878B (en) * 2019-12-10 2022-03-18 上海思源弘瑞自动化有限公司 Synchronous wave recording method and device
CN111579918B (en) * 2020-04-13 2022-03-04 上海紫通信息科技有限公司 Sampling synchronous correction method for transient recording type fault indicator
CN111366887A (en) * 2020-04-27 2020-07-03 广东电网有限责任公司电力调度控制中心 Mutual inductor transient sampling precision calibration method, device, terminal and medium
CN111679151B (en) * 2020-05-27 2022-06-21 武汉华电顺承科技有限公司 Synchronous alignment point vector calculation method for double-end distance measurement of wave recording data
CN112462172B (en) * 2020-11-11 2023-06-13 国网四川省电力公司电力科学研究院 An Adaptive Sampling Rate Power Quality Event Recording Method

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11352252A (en) * 1998-06-10 1999-12-24 Funai Electric Co Ltd Device and method for setting time based on broadcasting electronic wave and medium with program for setting time correct based on broadcasting electronic wave recorded therein
JP2004045355A (en) * 2002-07-16 2004-02-12 Japan Atom Energy Res Inst Waveform recorder that facilitates monitoring of large-scale facilities that are operating at high frequency
CN103199862B (en) * 2013-02-21 2016-06-22 四川电力科学研究院 A kind of method of data synchronization for modulus mixing record ripple
CN106501656A (en) * 2016-10-12 2017-03-15 国网上海市电力公司 A kind of on-line acquisition system of distribution line failure waveform
CN108226836B (en) * 2016-12-31 2021-06-01 科大智能电气技术有限公司 Calibration method for waveform synchronization of fault indicator
CN106950445B (en) * 2017-03-15 2019-12-20 北京四方继保自动化股份有限公司 Inter-station time difference analysis method based on fault recording data
CN107222280A (en) * 2017-05-24 2017-09-29 云南电力试验研究院(集团)有限公司 A kind of method that time synchronized is realized to different time scales fault recording system
CN107024632B (en) * 2017-06-07 2023-03-24 国网山西省电力公司电力科学研究院 Synchronization device and synchronization method for multi-measuring-point different wave recording equipment
CN107390089B (en) * 2017-07-19 2020-09-04 国网福建省电力有限公司 Method and device for realizing transient zero-sequence current signal synchronization in distributed system
CN107959969B (en) * 2017-12-11 2020-06-09 武汉中原电子信息有限公司 Time synchronization method applied to transient recording type fault indicator

Also Published As

Publication number Publication date
CN109828183A (en) 2019-05-31

Similar Documents

Publication Publication Date Title
CN109828183B (en) Waveform synchronous correction method and system suitable for transient recording type fault indicator
CN108173263B (en) A distribution network topology error identification algorithm based on AMI measurement information
CN108226836B (en) Calibration method for waveform synchronization of fault indicator
CN103743998B (en) One-phase earthing failure in electric distribution network localization method based on cross-correlation coefficient and system
CN108205096B (en) Multisource fault data analysis method based on active time correction
CN110161326A (en) A kind of more network live nuclear-phase methods and system at a distance
CN105044639A (en) Correction method and system for merging unit acquisition channels of power system
CN103884931B (en) A kind of substation bus bar part throttle characteristics test recording equipment
CN111381130A (en) T-connection line fault positioning method and system considering traveling wave velocity
CN104215838A (en) Remote nuclear phase method for intelligent substation
CN104865446B (en) Photovoltaic combining inverter exchanges outlet line impedance online test method
CN105388355A (en) Ground screen shunting vector test system and method with GPS synchronization
CN105261202A (en) Transmission delay measurement method based on B spline interpolation wavelet singularity detection
CN102788901A (en) High accuracy synchronous dynamic phasor measurement method
CN114925973A (en) Power grid running state obtaining method and system based on multi-source data integration
CN112557959B (en) Diagnosis and location method for single open circuit fault of switching tube of direct drive permanent magnet wind power converter
CN112858802B (en) Method, device and system for determining effective value of output electric parameter of switching power converter
CN110988599B (en) Distributed wave recording high-precision synchronization method for power distribution fault indicator
CN109298361B (en) Fault analysis preprocessing method based on transient recording type fault indicator
CN105429629A (en) FPGA-Based Phase-Locked Method and Its Phase-Locked Loop
CN102914378A (en) Method for monitoring temperature of generator rotor coil of nuclear power station based on DCS
CN109142985A (en) A kind of fault oscillograph system for Travelling Wave Fault Location
CN111443260A (en) Power grid phase difference detection method and system
CN114859284A (en) Test method and system for step response of current transformer
CN111682564B (en) Active synchronous control method and system for distributed energy storage fast smooth grid connection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant