一种加窗谱线插值的电力系统频率计算方法
【专利摘要】本发明提供了一种加窗谱线插值的电力系统频率计算方法，包括以下步骤：1)通过电压测量电路获得公共连接点的10个连续采样周期的电压采样数据；2)对电压采样数据进行时域加窗函数处理；3)对加窗后的数据进行局部离散傅里叶变换，获得特定谐波次数的电压幅值；4)对计算得到的特定电压幅值，在相邻两个最大幅值之间进行插值；5)通过多项式曲线拟合方法计算得到频率系数，从而确定电网的基波频率；本发明方法，对N个采样点加布莱克曼窗函数，并且只对指定的谐波频率进行局部离散傅里叶(DFT)变换，经试验证明，本发明方法抑制频谱泄漏效果好，计算速度快，计算过程简单、精度高、通用性强，可用于高精度电能计量和电能质量分析。
【专利说明】
-种加窗谱线插值的电力系统频率计算方法
技术领域
[0001] 本发明设及电能质量分析与控制领域，尤其是设及一种加窗谱线插值的电力系统 频率计算方法。
【背景技术】
[0002] 频率是电能生产、消费必须具备的质量检验指标，也是衡量电力系统运行状态的 重要参数，是电力系统最主要的特征量之一，因此，频率测量是电力系统测量中十分重要的 环节，若电网频率测量精度达不到要求将导致与频率相关的电气参数测量或计算严重偏离 实际值。
[0003] 电网频率测量技术研究的意义首先在于电网频率是电能质量体系中的一个重要 指标，其次电网频率测量不准将导致与频率相关的电气参数测量或计算出现较大误差；另 外在电力系统中，一些高新技术的应用需要高精度频率测量作保证，同时它还是后续谐波 分析的基础;频率测量主要基于硬件测频技术和软件测频技术;常用的软件测频算法有周 期法(或称交零法）、解析法、误差最小化原理类算法、傅里叶算法、正交去调制法等;周期法 在实时测量仪器中应用较多；实践证明周期法虽有原理简单、易于实现的优点，但其对噪 声、谐波及间谐波的抗干扰能力较弱；
[0004] 快速傅里叶变换(FFT)易于在嵌入式数字信号系统中实现，是频率软件测量的常 用方法，也是谐波分析的主要方法，但由于很难做到同步采样和整周期截断，存在频谱泄漏 和栅栏效应现象，影响频率测量精度;基于W上原因，现有发明[CN101852826A]中，采用基 于Blackman-化rris窗的四阶插值，对基波频率进行S次逼近;虽然精度达到要求，但该算 法的计算复杂程度大大增加了，运对于实时性要求很高的电力系统检测装置而言是不允许 的，且精度无法根据实际情况进行调整；同时，插值函数的阶次越高，计算精度有时未必就 越高，使用阶次较高的Blackman-化rris窗插值得到的电力系统频率与真实频率相比，有时 误差不小;而现有发明[CN104849545A]中，首先对电流信号采样，得到n时刻和n-1时刻的采 样矩阵，最终通过求取电压电流矩阵的特征根来确定相位及其系统频率;实际系统中，由于 电流信号初相角度变化通常较大，因此仅基于电流采样信号计算系统频率误差也较大，而 且该方法需要求取电流矩阵的广义逆矩阵，明显加大了计算量，运行时间加长。
【发明内容】
[0005] 针对现有技术中对电网系统频率的计算速度慢、精度低，本发明提供了一种加窗 谱线插值的电力系统频率计算方法，可对电网系统频率实现精确快速计算。
[0006] -种加窗谱线插值的电力系统频率计算方法，其中，包括W下步骤：
[0007] 1)采集获得公共连接点的10个连续工频周期的电压采样数据y(n);
[000引式中n = 0，l，…..，N-1，其中，n为整数，表示当前采样点;N为正整数，表示所采样 数据的总长度；
[0009] 之后，进入步骤2);
[0010] 2)对步骤I)得到的电压采样数据y(n)进行时域加窗函数处理，采用=项系数二阶 升余弦窗即布莱克曼窗函数对电压采样数据进行加窗处理，布莱克曼窗函数w(n)的时域表 达式为：
[0011]
[0012] 将步骤1)得到的电压采样数据y(n)和布莱克曼窗函数的离散序列相乘，即可得到 加窗处理后的数据x(n)，具体如下式所示：
[0013] x(n) =y(n)w(n)
[0014] 3)对从步骤2)得到的加窗处理后的数据x(n)进行局部离散傅里叶变换，计算第9、 第10和第11根谱线分量，局部离散傅里叶变换的计算公式为：
[0015]
[0016]
[0017]
[0018] 其中，
，k为整数，此处4 = 9、4=10、4=11，即计算第9、第10和第11根 谱线分量；j为虚数单位，e为自然指数函数的底数；丈(9)、方(iO巧ILY(Ii):分别是第9、10和 11根谱线分量，为复数形式，分别取绝对值即可得到第9、10和11谱线的幅值：
[0019]
[0020]
[0021]
[0022] 4)由于理想电力基波频率为50化，实际电力频率在50化周围变化，步骤3)计算得 到的第9、10和11根谱线对应幅值X(9)、X(10)和X(Il)中，X(IO)是最大的，作为最大值Xmaxl =X( 10 )，令常数kmaxl = 10 ;再选出X ( 9 )和X ( 11 )中的较大值，作为次大值Xmax2，则：
[0023] 当 X(9)〉X(11)时，Xmax2 = X(9)，令常数 kmax2 = 9;
[0024] 当 X(9)<X(11)时，Xmax2 = X(ll)，令常数 kmax2=ll;
[002引设km是电力系统实际频率所对应的谱线系数，令a = km-kmaxl-0 . 5，在最大值Xmaxl和 次大值Xmax2之间进行插值，插值的表达式为：
[0026]
[0027] 其中e是比例参数，W(2时)是布莱克曼窗函数。由于X化maxi)和X化max2)是已知的，可 W计算出e值；
[002引 5)对步骤4)得到的0表达式，计算电力系统实际频率所对应的谱线系数km = a+kmaxl +0.5,其中：
[00 巧]a = 1.9600+0.153 炉+0.074护
[0030] 则电网电压的实际频率为fm = 5*km。
[0031 ]与现有技术相比，本发明具有W下优点：
[0032] -、数据长度小，降低了对于硬件内存的要求:所提方法仅需要10个连续采样周期 的电压采样数据，数据长度小，降低了对于硬件内存的要求；
[0033] 二、计算过程更加简便，计算量较小，节省运算时间和物理内存:进行局部离散傅 里叶变换，仅计算第9、10和11根谐波谱线幅值，无需计算其他谱线幅值，计算量较小，节省 运算时间和物理内存；
[0034] =、计算精度较高，实际运行误差很小:使用基于布莱克曼窗插值的离散傅里叶变 换在一定程度上消除了频谱泄漏和栅栏效应的影响，大大提高了计算精度。
【附图说明】
[0035] 图1为本发明提供的一种加窗谱线插值的电力系统频率计算方法流程图；
[0036] 图2为本发明加窗谱线插值的电力系统频率计算方法步骤4)的谐波谱线图。
【具体实施方式】
[0037] 下面结合流程图和具体实施例对本发明进行详细说明；
[0038] 本发明采用基于加窗谱线插值的频率计算方法，见图1所示流程图：
[0039] 1)获取电压数据
[0040] 根据IEC 61000-4-7标准，获得公共连接点的10个连续采样周期的电压数据；
[0041 ] 2)加窗函数
[0042] 对时域样本加窗函数处理;窗函数通常用来做频率分析;布莱克曼窗实质上是一 种=项系数二阶升余弦窗，其时域表达式如下：
[0043]
(1)
[0044] 式中n = 0，l，.... .，N-1;
[0045] 在本发明应用中，较为关注的是谐波幅值和相角大小的精确度，因此，选择具有最 大的主瓣宽度窗函数更适用于本发明，综合考虑后，本发明采用布莱克曼窗对电压采样数 据进行加窗处理；
[0046] 3)执行局部离散傅里叶变换
[0047] 该步骤是对采样数据进行局部离散傅里叶变换;对10个连续周期的N个采样点进 行离散傅里叶变换，表达式可W写成下式：
[004引
（2)
[0049] 式中，Wf指第k个DFT系数；
[00加]
㈱
[0051] N个点序列的DFT变换也可W用矩阵来计算，（相应采样周期的各采样数据按照时 间顺序组成点序列）
[0化2] 如；
[0化 3
(4)
[0054] 其中，X(k)是第k次谐波谱线的幅值，x(n)是输入数据(指的是相应采样周期的采 样数据）；
[0055] 离散傅里叶变换的频谱分辨率是：
[0056]
(、）
[0057] 式中，f S是电网理想基波频率50化，其中，0.02s是工频周期(理想基波频率50化的 倒数），10指的是上文所述10个采样周期，A指代离散傅里叶变换的频谱分辨率；
[005引工程上，电网频率一般在45化至55化之间，因此基波频率fm所对应的谱线一般落 在第9、第10根谱线之间或者第10、第11根谱线之间，如图2所示，因此，确定主频时只需对第 9、10和11根谱线分量进行计算，并不需要对所有的谱线分量进行计算;局部离散傅里叶变 换的表达式为：
[0化9] (6)
[0060]
[0061]
[0062] 其中，
，x(n)是加窗数据样本值，X(9)、X(IO)和乂(川分别是第9、10 和11根谱线分量，为复数形式，分别取绝对值即可得到第9、10和11谱线的幅值：
[0063]
[0064]
[00 化]
[0066] 4)确定基波频率
[0067] 该步骤是确定基波频率；从步骤3)获得第9、10和11根谱线幅值，分别为X(9)、X (10)和X(11)。由于理想电力基波频率为50Hz，实际电力频率在50化周围变化。步骤3)计算 得到的第9、10和11根谱线对应幅值X(9)、X(10)和X(Il)中，X(IO)应是最大的，作为最大值 Xmaxl = X( 10 )，令谱线系数kmaxl = 10 ;再选出X(9 )和X( 11 )中的较大值，作为次大值Xmax2，则： [006引当 X(9)〉X(11)时，Xmax2 = X(9)，令谱线系数 kmax2 = 9;
[0069] 当 X(9)<X(11)时，Xmax2 = X(ll)，令谱线系数 kmax2=ll;
[0070] 由工程经验可知，电力系统的基波频率通常是波动的，令电力系统的基波频率为 fm，则：
[0071]
[0072] 式中，km是电力系统基波频率对应的谱线系数，通常不是整数，应介于kl和k2之间。
[0073] 本发明中，提出一种方法来确定系统主频，即采用在最大值Xmaxl和次大值Xmax2之间 进行插值，插值的表达式为：
[0074]
巧）
[0075] 其中0是比例参数，W(2村）是布莱克曼窗函数，f表示谐波频率；由于X(kmaxl)和X 化max2)是已知的，可W计算出0值：
[0076] 令
[0077] a = km-ki〇-〇.5 (9)
[0078] 计算电力系统实际频率所对应的谱线系数km = a+kmaxl+0.5,其中：
[00 巧]a = 1.9600+0.153 炉+0.074护（10)
[0080] 则电网电压的实际频率为fm = 5*km。
[0081] 图2为本发明加窗谱线插值的电力系统频率计算方法步骤4)的谐波谱线图。
[0082] 用计算机模拟验证本发明提出的确定基波频率方法的有效性，根据表1所示的 IEC-4-30标准来建立模拟系统：
[0083] 表1模拟参数
[0084]
[0085] 测试的输入信号可表示为：
[0086]
[0087] 式中u(t)表示时变电压信号，CO是可W进行设置的基波频率(角速度），U是基波电 压的幅值，0为谐波信号的初始相角(对应上表中的目），noise signal为噪声信号(对应上表 中的60地）；
[0088] 用本方法确定输入信号的基波频率，计算结果的误差如表2所示：
[0089] 表2仿真条件 r00901
[0091]理论上，根据标准IEC61000-4-30的模拟条件，本技术方案的方法具有很高的精确 度。
【主权项】
1. 一种加窗谱线插值的电力系统频率计算方法，其特征在于，包括以下步骤： 1) 采集获得公共连接点的10个连续工频周期的电压采样数据y(n); 式中n = 0,1，…..，N-1，其中，η为整数，表示当前采样点;N为正整数，表示所采样数据 的总长度； 之后，进入步骤2); 2) 对步骤1)得到的电压采样数据y(n)进行时域加窗函数处理，采用三项系数二阶升余 弦窗即布莱克曼窗函数对电压采样数据进行加窗处理，布莱克曼窗函数w(n)的时域表达式 为：将步骤1)得到的电压采样数据y(n)和布莱克曼窗函数的离散序列相乘，即可得到加窗 处理后的数据χ(η)，具体如下式所示： x(n)=y(n)w(n) 3) 对从步骤2)得到的加窗处理后的数据x(n)进行局部离散傅里叶变换，计算第9、第10 和第11根谱线分量，局部离散傅里叶变换的计算公式为：其中e为整数，此处k = 9、k= 10、k= 11，即计算第9、第10和第11根谱线分 量;j为虚数单位，e为自然指数函数的底数；义(9)、1卩0；)和1(：11：)分别是第9、10和11根谱线 分量，为复数形式，分别取绝对值即可得到第9、10和11谱线的幅值： X(9) = |Χ(9) JT(IO) = P(IO)I ^(ii) = |x(H) 4) 由于理想电力基波频率为50Hz，实际电力频率在50Hz周围变化，步骤3)计算得到的 第9、10和11根谱线对应幅值X(9)、X(10)和X(Il)中，X(IO)是最大的，作为最大值X maxl = X (10)，令常数kmaxl = 10;再选出X(9)和X(11)中的较大值，作为次大值Xmax2，则： 当 X(9)>X(11)时，Xmax2 = X(9)，令常数 kmax2 = 9; 当X(9)〈X(11)时，Xmax2 = X(Il)，令常数kmax2=ll; 设1^是电力系统实际频率所对应的谱线系数，令a = km-kmaxl-0.5，在最大值Xmaxl和次大 值Xmax2之间进行插值，插值的表达式为：其中β是比例参数，W ( 2lTf )是布莱克曼窗函数。由于X( kmaxl )和X ( kmax2 )是已知的，可以计 算出β值； 5)对步骤4)得到的β表达式，计算电力系统实际频率所对应的谱线系数1^ = 〇+1^+ 0.5,其中： a = 1.960β+0.153β3+0.074β5 则电网电压的实际频率为f m=5*km。
【文档编号】G01R23/16GK106018956SQ201610653629
【公开日】2016年10月12日
【申请日】2016年8月10日
【发明人】林顺富, 胡飞, 顾春艳, 高健飞, 陈中奎, 许亮峰
【申请人】北京妙微科技有限公司