CN113672863A - 一种发电机启机保护的简化相量计算方法及系统 - Google Patents
一种发电机启机保护的简化相量计算方法及系统 Download PDFInfo
- Publication number
- CN113672863A CN113672863A CN202110773937.2A CN202110773937A CN113672863A CN 113672863 A CN113672863 A CN 113672863A CN 202110773937 A CN202110773937 A CN 202110773937A CN 113672863 A CN113672863 A CN 113672863A
- Authority
- CN
- China
- Prior art keywords
- generator
- frequency
- phasor
- sampling
- analog quantity
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
本发明公开了一种发电机启机保护的简化相量计算方法及系统,本发明对发电机进行测频,根据发电机频率和采样频率调整数据窗,在调整后数据窗的基础上将傅氏算法推导成递归表达式,采用递归表达式计算模拟量离散化后的谐波相量,减少计算量,在递归计算过程中,每隔预设时间采用全波傅氏算法进行一次非递归校正,确保计算精度,实现了精度高、计算量少的谐波相量计算方法。便于滤除静止变频器启动时电压或电流中的整次谐波。
Description
技术领域
本发明涉及一种发电机启机保护的简化相量计算方法及系统,具体涉及一种适用于大型燃气发电机启机保护的简化相量计算方法及系统,属于发变组启机保护技术领域。
背景技术
大型燃气发电机常采用静止变频器(Static Frequency Convertor,SFC)启动。燃气发电机在启动过程中,输入到继电保护装置的二次电流、二次电压的频率变化范围较大,变频状态下的基波相量与各次谐波相量的计算较为困难。此外,SFC工作时输出电流谐波较大,影响大型燃气发电机启机保护的计算准确性。
为了能够在变频环境中滤除发电机或SFC产生的整次谐波电流,并对变频过程中的发电机电气量的基波进行精确测量,需要进行谐波相量计算。目前常用的方法为傅氏算法,该方法需要对一周波内所有有效数据点进行傅氏计算,虽然该种方法精度较高,但计算量偏大。
发明内容
本发明提供了一种发电机启机保护的简化相量计算方法及系统,解决了现有方法计算量偏大的问题。
为了解决上述技术问题,本发明所采用的技术方案是:
一种发电机启机保护的简化相量计算方法,包括:
利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
根据采样频率和发电机工作频率,计算每周波的采样点数;
对采样点数进行抽点处理,获得新的数据窗;
基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间或发电机频率变化超过阈值采用全波傅氏算法进行一次非递归校正。
利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率,具体过程为:
对发电机机端相间电压或中性点电流采样值数据窗进行低通滤波,获得测频数据窗;
根据测频数据窗进行过零点位置判别,获取相邻两过零点的时间差;
根据相邻两过零点的时间差,获得发电机工作频率。
优先利用发电机机端相间电压采样值数据窗进行测频,若发生三相短路故障、且发电机机端相间电压低于测频门槛值,利用中性点电流采样值数据窗进行测频。
采样点数抽点处理的公式为:
其中,Nf为抽点处理后每周波的采样点数,N为Nf的整数部分,ΔN为Nf的小数部分;
递归表达式为:
其中,时刻模拟量X的k次谐波相量,f0为发电机工作频率,Rk(t1)为t1时刻模拟量X的k次谐波相量的实部,Ik(t1)为t1时刻模拟量X的k次谐波相量的虚部,Rk(t1-Kn/fs)为t1-Kn/fs时刻模拟量X的k次谐波相量的实部,Ik(t1-Kn/fs)为t1-Kn/fs时刻模拟量X的k次谐波相量的虚部,T为待测模拟量X(t)的基波周期。
模拟量离散化后的谐波相量的公式为:
其中,为t1时刻模拟量X的k次谐波相量,模拟量X的采样点序列依次为x(n)、x(n-1)、x(n-2)、…,x(n)为当前最新采样点,x(n-i)表示x(n)之前i/fs时间的采样点,t1时刻位于x(n)、x(n-1)的采样时刻之间,rk(n)、ik(n)分别为的实部计算值和虚部计算值,Kn为抽点间隔,Nf为抽点处理后每周波的采样点数,N为Nf的整数部分,f0为发电机工作频率,fs为采样频率;
非递归校正的公式为:
一种发电机启机保护的简化相量计算系统,包括:
工作频率获取模块:利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
采样点数计算模块:根据采样频率和发电机工作频率,计算每周波的采样点数
抽点模块:对采样点数进行抽点处理,获得新的数据窗;
递归模块:基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间或发电机频率变化超过阈值采用全波傅氏算法进行一次非递归校正。
一种存储一个或多个程序的计算机可读存储介质,所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行发电机启机保护的简化相量计算方法。
一种计算设备,包括一个或多个处理器、一个或多个存储器以及一个或多个程序,其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行,所述一个或多个程序包括用于执行发电机启机保护的简化相量计算方法的指令。
本发明所达到的有益效果:本发明对发电机进行测频,根据发电机频率和采样频率调整数据窗,在调整后数据窗的基础上将傅氏算法推导成递归表达式,采用递归表达式计算模拟量离散化后的谐波相量,减少计算量,在递归计算过程中,每隔预设时间采用全波傅氏算法进行一次非递归校正,确保计算精度,实现了精度高、计算量少的谐波相量计算方法。便于滤除静止变频器启动时电压或电流中的整次谐波。
附图说明
图1为本发明方法的流程图;
图2为低通滤波器幅频特性示意图;
图3为相量实部的处理示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,一种发电机启机保护的简化相量计算方法,包括以下步骤:
步骤1,利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
步骤2,根据采样频率和发电机工作频率,计算每周波的采样点数;
步骤3,对采样点数进行抽点处理,获得新的数据窗;
步骤4,基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间或发电机频率变化超过阈值采用全波傅氏算法进行一次非递归校正。
上述方法运用在继电保护装置,对发电机进行测频,根据发电机频率和采样频率调整数据窗,在调整后数据窗的基础上将傅氏算法推导成递归表达式,采用递归表达式计算模拟量离散化后的谐波相量,减少计算量,在递归计算过程中,每隔预设时间采用全波傅氏算法进行一次非递归校正,确保计算精度,实现了精度高、计算量少的谐波相量计算方法,便于滤除静止变频器启动时电压或电流中的整次谐波。
继电保护装置先对发电机机端相间电压或中性点电流、机端零序电压进行采样,得到相间电压、机端零序电压的采样瞬时值,然后基于采样进行测频。
上述测频,优先利用发电机机端相间电压采样值数据窗,若发生三相短路故障、且发电机机端相间电压低于测频门槛值,利用中性点电流采样值数据窗;具体过程如下:
11)对发电机机端相间电压或中性点电流采样值数据窗进行低通滤波,获得测频数据窗;
12)根据测频数据窗进行过零点位置判别,获取相邻两过零点的时间差;
13)根据相邻两过零点的时间差,获得发电机工作频率。
以发电机机端相间电压采样值数据窗为例,机端AB相间电压采样值数据窗经过如图1所示的低通滤波器形成测频数据窗,然后对测频数据窗进行过零点位置判别,当找到过零点位置时,根据过零点两侧的采样值进行插值,求得过零点的精确时间,同样的方法查找下一个过零点时间,根据相邻两过零点的时间差,得到发电机工作频率f0;其中,低通滤波器要求0-60Hz信号不衰减,100Hz以上能够衰减至0.1倍以下。
根据发电机工作频率f0和采样频率fs(继电保护装置固有采样频率),可计算出每周波的采样点数,具体公式如下:
其中,Nall为抽点处理前每周波的采样点数;
对数据窗中的Nall进行抽点处理,抽点间隔为:
其中,m为整数,表示一周波采样点是否进行抽点处理的门槛值,其取值范围为12≤m≤24;
因此抽点处理后的采样点数为:
其中,Nf为抽点处理后每周波的采样点数,N为Nf的整数部分,ΔN为Nf的小数部分。
在时间不变的情况下,减小采样点,从而构成降低计算量的新数据窗。
在新数据窗的基础上,将全波傅氏算法推导成递归表达式,以t1时刻模拟量X的k次谐波相量为例:
推导得到递归表达式为:
其中,为t1时刻模拟量X的k次谐波相量,f0为发电机工作频率,Rk(t1)为t1时刻模拟量X的k次谐波相量的实部,Ik(t1)为t1时刻模拟量X的k次谐波相量的虚部,Rk(t1-Kn/fs)为t1-Kn/fs时刻模拟量X的k次谐波相量的实部,Ik(t1-Kn/fs)为t1-Kn/fs时刻模拟量X的k次谐波相量的虚部,T为待测模拟量X(t)的基波周期。
模拟量X的采样点序列依次为x(n)、x(n-1)、x(n-2)、…,x(n)为当前最新采样点,x(n-i)表示x(n)之前i/fs时间的采样点,假设t1时刻位于x(n)、x(n-1)的采样时刻之间,将该递归表达式按照等时间间隔离散化处理,并按照梯形法则对函数X(t)cos(2kπf0t)在时间轴上进行积分,求得f1,如图2所示,f1可用图2中两块阴影部分的面积之差求得;同理对函数X(t)sin(2kπf0t)在时间轴上进行积分,求得f2,则离散化后模拟量X的k次谐波相量的实部与虚部计算公式如下:
根据上式,模拟量的谐波相量的公式为:
为保证计算精度,每隔预设时间(40ms)或发电机频率变化超过阈值(0.2Hz)采用全波傅氏算法进行一次非递归校正,具体非递归校正的公式为:
上述算法中使用大量的sin和cos,为了节约计算量,采用拟合函数fsin和fcos替换,具体如下:
上述方法在变频环境下精度较高,响应速度较快,能够滤除整次谐波,且计算效率较高。
一种发电机启机保护的简化相量计算系统,包括:
工作频率获取模块:利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
采样点数计算模块:根据采样频率和发电机工作频率,计算每周波的采样点数
抽点模块:对采样点数进行抽点处理,获得新的数据窗;
递归模块:基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间采用全波傅氏算法进行一次非递归校正。
一种存储一个或多个程序的计算机可读存储介质,所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行发电机启机保护的简化相量计算方法。
一种计算设备,包括一个或多个处理器、一个或多个存储器以及一个或多个程序,其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行,所述一个或多个程序包括用于执行发电机启机保护的简化相量计算方法的指令。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。
Claims (10)
1.一种发电机启机保护的简化相量计算方法,其特征在于,包括:
利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
根据采样频率和发电机工作频率,计算每周波的采样点数;
对采样点数进行抽点处理,获得新的数据窗;
基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间或发电机频率变化超过阈值采用全波傅氏算法进行一次非递归校正。
2.根据权利要求1所述的一种发电机启机保护的简化相量计算方法,其特征在于,利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率,具体过程为:
对发电机机端相间电压或中性点电流采样值数据窗进行低通滤波,获得测频数据窗;
根据测频数据窗进行过零点位置判别,获取相邻两过零点的时间差;
根据相邻两过零点的时间差,获得发电机工作频率。
3.根据权利要求1或2所述的一种发电机启机保护的简化相量计算方法,其特征在于,优先利用发电机机端相间电压采样值数据窗进行测频,若发生三相短路故障、且发电机机端相间电压低于测频门槛值,利用中性点电流采样值数据窗进行测频。
6.根据权利要求1所述的一种发电机启机保护的简化相量计算方法,其特征在于,模拟量离散化后的谐波相量的公式为:
其中,为t1时刻模拟量X的k次谐波相量,模拟量X的采样点序列依次为x(n)、x(n-1)、x(n-2)、…,x(n)为当前最新采样点,x(n-i)表示x(n)之前i/fs时间的采样点,t1时刻位于x(n)、x(n-1)的采样时刻之间,rk(n)、ik(n)分别为的实部计算值和虚部计算值,Kn为抽点间隔,Nf为抽点处理后每周波的采样点数,N为Nf的整数部分,f0为发电机工作频率,fs为采样频率;
8.一种发电机启机保护的简化相量计算系统,其特征在于,包括:
工作频率获取模块:利用发电机机端相间电压或中性点电流采样值数据窗进行测频,获得发电机工作频率;
采样点数计算模块:根据采样频率和发电机工作频率,计算每周波的采样点数
抽点模块:对采样点数进行抽点处理,获得新的数据窗;
递归模块:基于新的数据窗,采用递归表达式计算模拟量离散化后的谐波相量;其中,递归表达式为全波傅氏算法推导而成的表达式,并且在递归计算过程中,每隔预设时间或发电机频率变化超过阈值采用全波傅氏算法进行一次非递归校正。
9.一种存储一个或多个程序的计算机可读存储介质,其特征在于:所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行根据权利要求1至7所述的方法中的任一方法。
10.一种计算设备,其特征在于,包括:
一个或多个处理器、一个或多个存储器以及一个或多个程序,其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行,所述一个或多个程序包括用于执行根据权利要求1至7所述的方法中的任一方法的指令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110773937.2A CN113672863A (zh) | 2021-07-08 | 2021-07-08 | 一种发电机启机保护的简化相量计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110773937.2A CN113672863A (zh) | 2021-07-08 | 2021-07-08 | 一种发电机启机保护的简化相量计算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113672863A true CN113672863A (zh) | 2021-11-19 |
Family
ID=78538673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110773937.2A Pending CN113672863A (zh) | 2021-07-08 | 2021-07-08 | 一种发电机启机保护的简化相量计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113672863A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT525717A1 (de) * | 2021-11-24 | 2023-06-15 | Sprecher Automation Gmbh | Verfahren zur Phasorbestimmung eines Strom- oder Spannungssignals |
-
2021
- 2021-07-08 CN CN202110773937.2A patent/CN113672863A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT525717A1 (de) * | 2021-11-24 | 2023-06-15 | Sprecher Automation Gmbh | Verfahren zur Phasorbestimmung eines Strom- oder Spannungssignals |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8108165B2 (en) | Acquiring phasors outside the frequency tracking range for power protective relays | |
Babu et al. | Analysis of SDFT based phase detection system for grid synchronization of distributed generation systems | |
CN107589299B (zh) | 基于多频率相量模型的电力信号同步相量测量方法 | |
CN111398731A (zh) | 基于多频率-泰勒模型滤除衰减直流的动态相量测量方法 | |
CN108896820A (zh) | 一种适用于静止变频器启动的调相机启机保护相量计算方法 | |
Luna et al. | Advanced PLL structures for grid synchronization in distributed generation | |
CN112505407A (zh) | 电网宽频振荡监测方法、系统、设备及可读存储介质 | |
CN113672863A (zh) | 一种发电机启机保护的简化相量计算方法及系统 | |
CN111190042B (zh) | 一种智能传感终端及对电网信号全频段测量的方法 | |
CN110412415B (zh) | 一种基于dft、多阶滤波和突变判据的同步相量计算方法 | |
CN111800055B (zh) | 一种双凸极电机平均转矩确定方法和装置 | |
CN110596455B (zh) | 一种工频电参数提取方法、系统及计算机可读存储介质 | |
CN115586345B (zh) | 一种基于定子电流的电机自适应转速估计方法及系统 | |
Johns et al. | Fundamental digital approach to the distance protection of EHV transmission lines | |
Nakano et al. | Frequency detection method based on recursive DFT algorithm | |
CN112595896B (zh) | 一种对地电容参数检测方法 | |
Yang et al. | A novel phase locked loop for grid-connected converters under non-ideal grid conditions | |
CN110535091B (zh) | 一种适应电气设备变频特性的微机保护方法和装置 | |
CN112462138A (zh) | 一种谐波测量方法及系统 | |
CN108037351A (zh) | 一种变频电源功率因素测量方法 | |
CN103176030A (zh) | 配电系统的间谐波检测方法 | |
CN112255457A (zh) | 适用于自动准同期装置的相角差测量方法 | |
CN104820129A (zh) | 一种基波正序有功电流的检测方法 | |
Kunac et al. | Grid Voltage Amplitude and Frequency Real-Time Estimation using Linear Kalman Filter | |
Bertocco et al. | Numerical algorithms for power measurements |
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 |