一种基于正序瞬时功率算法的宽频振荡判别系统及方法与流程

1.本发明涉及一种基于正序瞬时功率算法的宽频振荡判别系统及方法，属于电力系统自动化测量技术领域。


背景技术：

2.同步相量测量装置(pmu)利用卫星同步时钟系统为广域范围内的全网同步采样提供统一的采样脉冲和标准时间，使得各个站点之间有了相同的时间基准点和采样参考基准点，在同步采样和计算之后所得到的同步相量能准确描述实际系统的动态过程，为电力系统新型保护、测控、安全稳定控制提供了新的数据源。
3.传统pmu装置一般测量基频分量，而新型“双高”电力系统中，谐波和间谐波含量较多，为了更好地描述电力系统的动态过程，需要同步测量谐波和间谐波，并在电力系统出现宽频振荡时发出宽频振荡告警。
4.依据相关标准，宽频振荡的瞬时功率频率范围是100-300hz。但是，当电力系统中存在三相不平衡时，瞬时功率中将出现二倍工频的瞬时功率分量，该分量恰好位于宽频振荡的频率判别区间。现有的宽频测量装置在用功率判别宽频振荡时，无法对三相不平衡和宽频振荡导致的瞬时功率变化进行区分，从而会导致当电力系统中存在三相不平衡时，装置发生宽频振荡误告警。
5.同时，cpu只能进行串行计算，在同时测量多元件的谐波和间谐波以及计算多元件的正序瞬时功率时，运算量大，存在运算速度慢、cpu负荷过高的问题。发明专利《一种基于fpga硬件dft递推的同步相量计算方法》(专利号：201210310767.5)和《一种快速同步相量修正方法》(专利号：201410529646.9)中均阐述了利用fpga进行多元件同步相量测量并行计算的方法，但这些方法并未解决多元件的电压、电流正序分量和正序瞬时功率并行计算问题。
6.公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域普通技术人员所公知的现有技术。


技术实现要素：

7.本发明的目的在于克服现有技术中的不足，提供一种基于正序瞬时功率算法的宽频振荡判别系统及方法，采用正序瞬时功率判别宽频振荡，可以滤除因三相不平衡导致的二倍工频谐波成分，同时还可以滤除该成分叠加至宽频振荡瞬时功率中后引发的新的频率成分。
8.为达到上述目的，本发明是采用下述技术方案实现的：
9.一方面，本发明公开了一种基于正序瞬时功率算法的宽频振荡判别系统，包括，
10.a/d采样模块，用于获取模拟信号，根据所述模拟信号，得到同步采样数字信号；
11.fpga模块，用于根据同步采样数字信号，计算正序瞬时功率，得到第二运算结果；
12.cpu模块，用于根据第二运算结果，得到宽频振荡的判别结果；
13.其中，所述fpga模块分别与a/d采样模块和cpu模块之间采用并行数据总线连接。
14.进一步的，所述fpga模块包括，
15.同步单元，用于获取b码时间信号，得到1pps信号；
16.第一fft单元，用于根据同步采样数字信号和1pps信号，基于fft运算，得到第一运算结果；
17.电压电流缓存单元，用于转存第一运算结果；
18.第二fft单元，用于根据第一运算结果，基于fft运算，得到第二计算结果；
19.正序瞬时功率缓存单元，用于转存第二运算结果。
20.进一步的，所述第一运算结果，包括打上时标的三相电压的谐波、间谐波幅值信息和相位信息，以及三相电流的谐波、间谐波幅值信息和相位信息。
21.进一步的，所述第二运算结果，包括打上时标的正序瞬时功率的幅值信息。
22.进一步的，所述fpga模块还包括，
23.第一补偿单元，设于第一fft单元与电压电流缓存单元之间，用于修正第一运算结果。
24.进一步的，所述fpga模块还包括，
25.第二补偿单元，设于第二fft单元与正序瞬时功率缓存单元之间，用于修正第二运算结果。
26.进一步的，所述fpga模块还包括，
27.dft单元，用于根据同步采样数字信号和1pps信号，基于dft运算，得到相量数据；
28.相量存储单元，用于存储相量数据。
29.另一方面，本发明公开了基于正序瞬时功率算法的宽频振荡判别系统的判别方法，包括如下步骤：
30.获取模拟信号，根据所述模拟信号，得到同步采样数字信号；
31.根据同步采样数字信号，计算正序瞬时功率，得到第二运算结果；
32.根据第二运算结果，得到宽频振荡的判别结果。
33.与现有技术相比，本发明所达到的有益效果：
34.本发明采用正序瞬时功率判别宽频振荡，可以滤除因三相不平衡导致的二倍工频谐波成分，同时还可以滤除该成分叠加至宽频振荡瞬时功率中后引发的新的频率成分。滤除三相不平衡导致的频率成分后，可以避免宽频振荡告警的误报，使得宽频振荡告警功能能够准确地在电力系统安全稳定运行的过程中发挥作用，更好地满足用户的需求。
35.本发明利用fpga的同步并发计算能力，可以实现多个元件本时刻电压、电流的谐波和间谐波与上一时刻电压、电流的正序分量以及正序瞬时功率的谐波和间谐波的快速计算。上述测量过程在卫星同步时钟同步后，所计算的广域范围内不同电气元件的谐波、间谐波信号可以直接进行比较，方便定位低频振荡、次超同步振荡、宽频振荡的振荡源。
附图说明
36.图1是一种基于正序瞬时功率算法的宽频振荡判别系统的模块原理图；
37.图2是fpga模块的原理图；
38.图3是一种基于正序瞬时功率算法的宽频振荡判别系统的方法流程图；
39.图4是使用本发明判别正常状态的结果；
40.图5是使用本发明判别三相不平衡状态的结果；
41.图6是使用本发明判别宽频振荡状态的结果。
具体实施方式
42.下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
43.实施例
44.本实施例公开了一种基于正序瞬时功率算法的宽频振荡判别系统，如图1-3所示，包括，
45.a/d采样模块，用于获取模拟信号，根据模拟信号，得到同步采样数字信号；
46.fpga模块，用于根据同步采样数字信号，计算正序瞬时功率，得到第二运算结果；
47.cpu模块：用于根据第二运算结果，得到宽频振荡的判别结果；
48.其中，fpga模块分别与a/d采样模块和cpu模块之间采用并行数据总线连接。
49.需要说明的是，a/d采样模块中获取的模拟信号是通过ac模块变换整形得到的。ac模块用于将原始采样得到的电压、电流的模拟信号进行变换整形，使其转变为满足a/d采样要求的
±
5v弱电电压信号后，再输入至a/d采样模块。
50.本实施例中的ac模块的额定相电压是57.74v，额定相电流是5a/1a。
51.具体的，fpga模块内设有同步单元并使之驯服于卫星同步时钟，用于获取b码时间信号，即卫星同步时钟信号，得到1pps信号。
52.需要说明的是，此处的1pps信号还需要进行倍频处理，得到同步采样脉冲和同步计算脉冲。其中，同步采样脉冲用于控制a/d采样模块完成广域同步采样过程，同步计算脉冲用于后续的fft运算。
53.a/d采样模块，将获取模拟信号转换为数字信号，转换过程在fpga和外部卫星同步时钟控制下进行，得到同步采样数字信号。
54.本实施例中，pmu需接入b码对时，装置采集指定电压、电流的模拟信号时，接入b码对时，实现同步采样，1pps时刻和采样时刻重合，采样频率为1024hz。
55.转换后，对应的同步采样数字信号x(k)表示为：
[0056][0057]
其中，分别为第i个分量的幅值、频率、初相位，k为当前采样点序号，δt为采样间隔。
[0058]
进一步的，fpga模块包括，
[0059]
第一fft单元，用于根据同步采样数字信号和1pps信号，基于fft运算，得到第一运算结果；
[0060]
电压电流缓存单元，用于转存第一运算结果；
[0061]
第二fft单元，用于根据第一运算结果，基于fft运算，得到第二计算结果；
[0062]
正序瞬时功率缓存单元，用于转存第二运算结果。
[0063]
进一步的，fpga模块还包括，
[0064]
第一补偿单元，设于第一fft单元与电压电流缓存单元之间，用于修正第一运算结果。
[0065]
第二补偿单元，设于第二fft单元与正序瞬时功率缓存单元之间，用于修正第二运算结果。
[0066]
dft单元，用于根据同步采样数字信号和1pps信号，基于dft运算，得到相量数据；
[0067]
相量存储单元，用于存储相量数据。
[0068]
具体的，第一fft单元在上述1pps信号同步下，对所得到的同步采样数字信号进行fft运算，得到打上高精度时标的第一运算结果。通过第一补偿单元对第一运算结果进行修正再存入电压电流数据缓存单元。
[0069]
第一运算结果即工频分量和各谐波、间谐波分量的频率、幅值和相位信息，包括打上时标的三相电压的谐波、间谐波幅值信息和相位信息，以及三相电流的谐波、间谐波幅值信息和相位信息。
[0070]
第一fft单元的工作原理为：对同步采样数字信号进行fft运算，窗内采样点数为n1，得出三相电压谐波、间谐波幅值的频谱序列，得出三相电压谐波、间谐波幅值的频谱序列和相位的频谱序列和相位的频谱序列以及三相电流谐波、间谐波幅值的频谱序列和相位的频谱序列使用采样中间点的时间作为该频谱序列的绝对时标。
[0071]
第二fft单元在上述1pps信号同步下，读取电压电流数据缓存单元的第一运算结果，进行fft运算，得到打上高精度时标的第二运算结果。通过第二补偿单元对第二运算结果进行修正再存入正序瞬时功率缓存单元。
[0072]
第二运算结果即正序瞬时功率中各谐波和间谐波分量的频率和幅值信息，包括打上时标的正序瞬时功率的幅值信息。
[0073]
第二fft单元的工作原理为：基于第一运算结果中的三相电压、电流频谱序列计算电压、电流正序分量序列，进而计算正序瞬时功率。
[0074]
第k个正序电压分量和正序电流分量的计算方法为：
[0075][0076]
[0077][0078][0079][0080][0081]
其中，e表示自然常数，j表示复数的虚部，分别表示a、b、c三相电压和电流中第k个分量的复频域值，满足：
[0082][0083][0084][0085][0086][0087][0088]
因此，第k个正序分量可以表示为：
[0089]ua1k
(t)＝u
a1k
cos(2πfkt α
a1k
)
[0090]ub1k
(t)＝u
b1k
cos(2πfkt α
b1k
)
[0091]uc1k
(t)＝u
c1k
cos(2πfkt α
c1k
)
[0092]ia1k
(t)＝i
a1k
cos(2πfkt β
a1k
)
[0093]ib1k
(t)＝i
b1k
cos(2πfkt β
b1k
)
[0094]ic1k
(t)＝i
c1k
cos(2πfkt β
c1k
)
[0095]
其中，u
a1k
、α
a1k
分别是幅值和相位，u
b1k
、α
b1k
分别是的幅值和相位，u
c1k
、α
c1k
分别是的幅值和相位，i
a1k
、β
a1k
分别是的幅值和相位，i
b1k
、β
b1k
分别是的幅值和相位，i
c1k
、β
c1k
分别是的幅值和相位，t是时间。
[0096]
因此，上一数据窗内的时刻对应的瞬时功率为
[0097][0098]
进一步地，fpga在1pps信号同步下，对得到的正序瞬时功率进行fft运算，窗内采样点数为n2，得出正序瞬时功率幅值的频谱序列为使用采样中间点的时间
作为该频谱序列的绝对时标。
[0099]
cpu模块的工作原理如下：
[0100]
取100-300hz内频谱序列分量的最大值，假设为pk，若满足下式，则判定为发生宽频振荡，装置发出宽频振荡告警。
[0101][0102]
式中，pm是宽频振荡功率一次侧门槛值，ku为一二次测的电压变比，ki为一二次测的电流变比。
[0103]
在本实施例中，正常状态下的电压、电流只含有工频分量，其幅值大小等于额定值，相位互差120
°
；三相不平衡状态下的电压、电流仍只含有工频分量，和正常状态不同的是，a相电压幅值变为40v；宽频振荡状态下，电压只含有工频分量，电流除了工频分量外，还含有175hz、204hz和220hz的谐波分量。通过探讨正常状态、三相不平衡状态以及宽频振荡状态三种情形下的测量，得出的正常状态的振荡检测结果如图4所示，三相不平衡状态的振荡检测结果如图5所示，宽频振荡状态的振荡检测结果如图6所示。
[0104]
由图4-6所示的三种情形下的正序瞬时功率频谱图可知，本发明的判别方法能够正确判别宽频振荡的发生并发出告警，且不会将三相不平衡状态误判为宽频振荡状态。
[0105]
本发明着重阐述基于fpga正序瞬时功率算法的宽频振荡判别方法在pmu中的应用过程，实际中，该方法还可以在宽频测量装置、电能质量检测装置等场景中应用。
[0106]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0107]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0108]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0109]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0110]
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应
视为本发明的保护范围。