adaline神经网络和FFT补偿的闪变值测量方法与流程

adaline神经网络和fft补偿的闪变值测量方法
技术领域
1.本技术涉及电力系统领域，尤其涉及一种基于adaline神经网络和fft补偿的闪变值测量方法。


背景技术：

2.随着当前世界形势的不断复杂变化以及我国国民经济的高速健康发展，现代社会对电力供应的稳定可靠的性能要求也越来越高，即是对电能质量提出新要求，电力的电能质量作为工业的主导地位比以往任何时候都更加重要。目前随着工业技术的发展，大量非线性波动性负荷广泛应用，同时，多种形式的具有不规律特性的新能源逐渐并网。这些都会向电网注入多种波动，波动会引起白炽灯工作异常也即是电压闪变。闪变情况会导致精密器件无法正常工作以至于造成重大经济损失。采用治理措施的基本依据是闪变的准确测量，因此需要对电压闪变进行精准测量。
3.现有的闪变测量算法主要有iec推荐的平方检波法和傅里叶变换fft。平方检波法具体实现是将原始信号平方运算之后通过滤波器得到调幅波，这种方法计算量大，实现复杂。fft由于傅里叶变换存在频谱泄露从而导致栅栏效应，最终导致计算结果不准确，现有的工作有的提出了一些fft改进方法，使用加窗插值来减小频谱泄露，但是栅栏效应依然存在，而且结果十分依赖窗函数大小，因此需要很高的计算成本结果依然存在误差。


技术实现要素：

4.本技术提供了一种基于adaline神经网络和fft补偿的闪变测量方法，解决了现有技术存在的计算成本高计算误差大的问题。
5.为解决上述技术问题，本技术提供了一种基于adaline神经网络和fft补偿的初始闪变值测量方法，包括：
6.根据adaline神经网络和预先建立的实时电网电压闪变模型提取电压闪变的包络信号；
7.依据fft变换对所述包络信号进行计算以得初始幅值和初始频率，并用pso算法对所述初始幅值和所述初始频率进行修正以得所述包络信号的估计幅值和估计频率；
8.依据所述估计幅值和所述估计频率计算得出闪变包络信号，并对所述闪变包络信号采用国标闪变值计算方法进行计算，并进行误差补偿得闪变值。
9.优选地，所述实时电网电压闪变模型为：
10.u(t)＝a[1+v(t)]cos(2πf0t+θ0)
[0011]
式中，a为电网电压幅值；f0为工频载波电压的频率；θ0为电压的初始相位；v(t)为u(t)假想零线的电压瞬时值的包络波相对于a的大小，考虑正弦波电压波动，v(t)表达式为：
[0012][0013]
其中，n为电压闪变分量的数量，mn为电网电压闪变的第n个分量的幅值，fn为电网
电压闪变的第n个分量的频率，θn为第n个分量相对于工频电压的角度。
[0014]
优选地，所述实时电网电压闪变模型为：
[0015]
u(t)＝w1x1+w2x2[0016]
式中，自适应滤波器权值w1和w2分别为(a
en cos(θ0))和(-a
en sin(θ0))，输入信号x1和x2的值分别为cos(2πf0t)和sin(2πf0t)。
[0017]
优选地，所述adaline神经网络的系统函数为：
[0018][0019]
式中，为adaline系统的输出，为w1的估计值，为w2的估计值；
[0020]
t时刻的预测误差e(t)为：
[0021][0022]
式中，u(t)为电力系统实测电压；
[0023]
对应地，根据预测误差和自适应规则调整加权系数，直到e(t)为零时，得出所述电压闪变的包络信号a
en
为：
[0024][0025]
优选地，根据最小均方lms算法作为自适应规则对所述adaline神经网络的系统函数中的和进行更新，公式如下：
[0026][0027]
式中，x(k)为输入向量的k个样本，e(k)为k个样本的误差信号，α为学习速率。
[0028]
优选地，所述依据fft变换对所述包络信号进行计算以得和初始估计幅值m
range
和初始频率f
range
的范围为：
[0029][0030]mm
≤m
range
≤mm+γmm[0031]
式中，tw为时间窗口，γ为小于1.0的常数。
[0032]
优选地，所述用pso算法对所述初始幅值和所述初始频率进行修正为：
[0033]
对所述pso算法中的速度和位置进行更新，公式为：
[0034]vi,k+1
＝w
pvik
+c1r1(p
i,k-x
i,k
)+c2r2(p
g,k-x
i,k
)
[0035]
x
i,k+1
＝x
i,k
+v
i,k+1
[0036]
式中，v
i,k
为第k代第i个粒子的速度矢量，x
i,k
为位置矢量，c1和c2是影响粒子速度的学习因子，r1和r2是介于0和1之间的随机数，系数w
p
减小公式如下：
[0037][0038]
式中，w
max
为w
p
上限值，w
min
为w
p
下限值，iter是迭代次数，iter
max
为最大迭代次数。
[0039]
优选地，所述依据所述估计幅值和所述估计频率计算得出闪变包络信号为：
[0040]
[0041]
式中，为mn的估计值，为fn的估计值。
[0042]
优选地，所述对所述闪变包络信号采用国标闪变值计算方法进行计算，并进行误差补偿得闪变值为：
[0043][0044]
其中，a为工频电压的幅值，为pso算法估计得出的闪变幅值，通过pso算法估计的频率查表或者带入多项式进行计算获得补偿值c和标准波动值d
lim
。
[0045]
由以上技术方案可知，本技术所提供的一种基于adaline神经网络和fft补偿的闪变值测量方法，通过使用adaline神经网络提取电网电压中的波动量；通过fft变换与pso算法相结合的方式，可以准确估计出波动量的幅值与频率；最后考虑了相应频率的闪变值补偿误差对得出的闪变包络信号进行补偿，得出最终的闪变值，在降低计算成本的时提高了计算精度。
附图说明
[0046]
为了更清楚的说明本技术的技术方案，下面将对实施例中所需要使用的附图作简要的介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性的前提下，还可以根据这些附图获得其他的附图。
[0047]
图1为本发明实施例所提供的一种基于adaline神经网络和fft补偿的闪变值测量方法流程图；
[0048]
图2为本发明实施例所提供的一种adaline神经网络的结构图；
[0049]
图3为本发明实施例所提供的一种pso算法的流程图。
具体实施方式
[0050]
为了使本技术领域的人员更好地理解本技术中的技术方案，下面将结合附图，对本技术实施例中的技术方案进行清楚完整的描述。
[0051]
本技术的核心是提供一种基于adaline神经网络和fft补偿的闪变值测量方法，可以解决现有技术存在的计算成本高计算误差大的问题。
[0052]
图1为本发明实施例所提供的一种基于adaline神经网络和fft补偿的闪变值测量方法流程图，图2为本发明实施例所提供的一种adaline神经网络的结构图；图3为本发明实施例所提供的一种pso算法的流程图。
[0053]
如图1至图3所示。
[0054]
在执行之前，需要构建实时电网电压闪变模型。闪变是电网中电压波动作用的结果，因此，考虑波动分量，可以建立如下的实时电网电压闪变模型：
[0055]
u(t)＝a[1+v(t)]cos(2πf0t+θ0)
ꢀꢀꢀ
(1)
[0056]
式中，a为电网电压幅值；f0为工频载波电压的频率；θ0为电压的初始相位；v(t)为u(t)假想零线的电压瞬时值的包络波相对于a的大小，考虑正弦波电压波动，v(t)表达式为：
[0057]
[0058]
其中，n为电压闪变分量的数量，mn为电网电压闪变的第n个分量的幅值，fn为电网电压闪变的第n个分量的频率，θn为第n个分量相对于工频电压的角度。
[0059]
电压闪变信号实质是叠加了无数个正弦信号的结果，式1电网电压闪变模型可以变形如下，a
en
为电压闪变的包络信号：
[0060][0061][0062]
做完上述工作之后，一种基于adaline神经网络和fft补偿的闪变值测量方法，该方法包括以下步骤：
[0063]
s101：根据adaline神经网络和预先建立的实时电网电压闪变模型提取电压闪变的包络信号。
[0064]
adaline是一个具有多个节点、一个输出和多个可调权值的单层神经网络。如果对一个信号直接使用fft时，有必要以远高于信号频率的采样频率对信号进行采样，并且在时间窗内使用的内存需要足够大。因此使用adaline系统提取电力系统电压的包络信号，从而有效地降低fft计算的计算成本和内存使用。具体操作为，调整adaline权值，使其跟踪电力系统电压的包络线。
[0065]
为了使用adaline，将式(2)变为式(3)，将式(1)变为式(5)
[0066]
u(t)＝w1x1+w2x2ꢀꢀ
(5)
[0067]
式中，自适应滤波器权值w1和w2分别为(a
en cos(θ0))和(-a
en sin(θ0))，输入信号x1和x2的值分别为cos(2πf0t)和sin(2πf0t)。
[0068]
adaline神经网络的系统函数为：
[0069][0070]
式中，为adaline系统的输出，为w1的估计值，为w2的估计值，e(t)为t时刻的预测误差；
[0071][0072]
式中，u(t)为电力系统实测电压；
[0073]
对应地，根据预测误差和自适应规则调整加权系数，直到e(t)为零时，得出包络信号a
en
为：
[0074][0075]
在adaline神经网络中，使用最小均方(lms)算法，即widro-hoff delta规则作为自适应规则。对于第k个样本，权值更新公式如下：
[0076][0077]
式中，x(k)为输入向量的k个样本，e(k)为kk个样本的误差信号，α为学习速率。
[0078]
s102：依据fft变换对包络信号进行计算，此时是粗略估计，得出初始幅值和初始
频率，并用pso算法对初始幅值和初始频率进行修正以得包络信号的估计幅值和估计频率。
[0079]
如果电网电压没有闪变，则a
en
是常数，a
en
的fft在零频率处只有一个非零分量。反之，电压闪变的频率范围为0.5-35hz，对于其中的某一频率估计时，由于fft自身存在的频谱泄露问题，会导致栅栏效应。求出的值并不准确，设最大幅值为mm，对应频率为ff。真实闪变的频率范围和相应的幅值(frange，mrange)为：
[0080][0081]mm
≤m
range
≤mm+γmmꢀꢀ
(11)
[0082]
式中，tw为时间窗口，γ为小于1.0的常数。
[0083]
为了准确估计幅值和频率，采用pso算法粒子群优化算法属于进化算法的一类，由kennedy和eberhart在1995年提出。近年来，这类优化算法在电力系统优化问题中得到了越来越多的应用。粒子群优化算法是一组在搜索空间中扩散的粒子(优化变量)。每个粒子都可能是一个潜在的解。很明显，有些粒子比其他粒子有更好的位置，而其余的粒子则试图将自己的位置提升到优势粒子的位置。与此同时，优越粒子的位置也发生了变化。每个粒子位置的变化是基于粒子在先前运动中的经验和相邻粒子的经验。事实上，每个粒子都知道自己对邻近的粒子以及对整个粒子群的优越性和非优越性。
[0084]
s103：依据估计幅值和估计频率计算得出闪变包络信号，并对闪变包络信号采用国标闪变值计算方法进行计算，并进行误差补偿得闪变值。
[0085]
依据估计幅值和估计频率计算得出闪变包络信号为：
[0086][0087]
式中，为mn的估计值，为fn的估计值。
[0088]
在第k个样本的每个时刻，将真实闪烁分量与估计分量的差值提供给算法，并根据fft算法获得的目标函数和优化变量的范围，计算出最佳的估计。
[0089]
在粒子群算法中，粒子寻找最佳位置，直到遇到相对不变的状态或超过计算限制。粒子群算法的关键步骤是在每一个时间步长的时候，改变每个粒子朝向其个体历史最优位置pbest和群体历史最优gbest位置的速度。随机生成一个加速度，对于朝向pbest和gbest位置的速度产生单独的随机数。pso的过程可以描述如下：
[0090]
1.在目标空间的d维空间中初始化一个随机位置和速度的粒子群(数组)，第i个粒子的位置向量为速度向量为
[0091]
2.对于每个粒子，在d个变量中计算适应度。
[0092]
3.比较粒子的适应值也即是目标函数和它迄今为止搜索到的最优位置pbest。如果当前粒子所在的位置计算的目标函数值优于pbest，则将pbest位置设置为d维空间中的当前位置。
[0093]
4.将适应值与整个粒子群迄今为止搜索到的总体最佳值进行比较。如果当前值比gbest好，则重置gbest为当前粒子的位置。
[0094]
5.用pso算法对初始幅值和初始估计频率进行修正为：
[0095]
根据式(8)对pso算法中的速度和位置进行更新，公式为：
[0096]vi,k+1
＝w
pvik
+c1r1(p
i,k-x
i,k
)+c2r2(p
g,k-x
i,k
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0097]
x
i,k+1
＝x
i,k
+v
i,k+1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0098]
式中，v
i,k
为第k代第i个粒子的速度矢量，x
i,k
为位置矢量，c1和c2是影响粒子速度的学习因子，r1和r2是介于0和1之间的随机数。
[0099]
在实际使用中，若w
max
设置过大可能使粒子越过最小值点，相反的w
max
设置过小可能使粒子绕其位置旋转，从而无法搜索目标空间。w
max
的数量通常选择在变量范围的10％到20％之间。c1和c2选择范围为c1+c2＜4。此外，适当选择w的结果可以使得算法以较小的重复次数，达到最优点。在典型的pso算法实现过程中，系数w从0.9减小到0.4，其减小公式如下：
[0100][0101]
式中，w
max
为上限值，w
min
为是下限值，iter是迭代次数，iter
max
为最大迭代次数。
[0102]
6.返回到式2，直到满足一个标准，通常是一个足够好的适合度或最大的迭代数量，最后输出
[0103]
根据《gbt 12326-2008电能质量电压波动和闪变》中规定的计算方法
[0104][0105]
便可以根据估计所得的频率查表1对应的d
lim
，结合估计的幅值即可求出闪变值。由于实际操作中存在衰减，故增加了闪变值补偿环节。求出的闪变值与补偿值相乘就能够求出准确的闪变值。当求出的频率不在表1范围内时，可用表2的多项式系数求得对应的值，再进行运算。
[0106][0107]
最后可求得闪变值为：
[0108][0109]
其中，a为工频电压的幅值，为pso算法估计得出的闪变幅值，通过pso算法估计的频率查表或者带入多项式进行计算获得补偿值c和标准波动值d
lim
。
[0110]
表1不同频率下单位闪变波动值和补偿值
[0111]
f/hzd
lim
/％补偿值cf/hzd
lim
/％补偿值c0.52.341.025641026100.261.08342361911.4321.01522842610.50.271.0869565221.51.080.9765625110.2821.09051254120.8820.98619329411.50.2961.0940919042.50.7540.987166831120.3121.09529025230.6540.988142292130.3481.0989010993.50.5681.006036217140.3881.12739571640.51.016260163150.4321.148105626
4.50.4461.018329939160.481.15740740750.3981.031991744170.531.1723329435.50.361.028806584180.5841.20048019260.3281.052631579190.641.2330456236.50.31.052631579200.71.25156445670.281.048218029210.761.2755102047.50.2661.025641026220.8241.30890052480.2561.048218029230.891.3513513518.80.251.058201058240.9621.3869625529.50.2541.067235859251.0421.416430595
[0112]
表2多项式拟合系数
[0113][0114]
本领域技术人员在考虑说明书及实践这里公开的申请后，将容易想到本技术的其他实施方案。本技术旨在涵盖本技术的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本技术的一般性原理并包含本技术公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为实例性的，本技术的真正范围由权利要求指出。
[0115]
应当理解的是，本技术并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。以上所述的本技术实施方式并不构成对本技术保护范围的限定。