一种典型场景下谐波电压电流测量与特征分析方法与流程

1.本发明涉及电网技术领域，尤其涉及一种典型场景下谐波电压电流测量与特征分析方法。


背景技术：

2.快速准确地检测电力谐波分量是实现谐波治理的前提条件。电力平稳谐波信号分析主要采用快速傅里叶变换(fastfouriertransform，fft)法、加窗法、插值法、频谱校正等；对于非平稳信号的谐波分析有小波变换、希尔伯特黄变换(hilbert-huangtransform，hht)、神经网络、prony法、局域均值分解(localmeandecomposition，lmd)算法等。为了更好地指导工程实践，ieee和iec分别制订了相关国际标准，即ieeestd.519-1992和国际标准iec 61000-4-7。在这些标准中，傅里叶分析作为建议的标准信号处理工具，用于分析电压和电流波形的谐波含量。
3.但快速傅里叶变换方法存在以下问题：
4.(1)由于非同步采样和时域截断，常常会产生频谱泄露和栅栏效应，造成谐波检测精度达不到相关标准；
5.(2)采用快速傅里叶变换方法进行谐波分析时，还需考察其对电网频率波动、不同强度噪声以及含间谐波的情况的适应性，以及考虑应用场景对计算速度的要求。


技术实现要素：

6.本发明的目的在于克服上述技术的不足，提出一种典型场景下谐波电压电流测量与特征分析方法，该方法针对若干典型场景，提出充分考虑其负荷时空动态特性的谐波测量方案，并根据该方案可提出适用于各场景谐波特性的改进谐波分析方法，准确提取谐波分量。
7.为解决以上技术问题，本发明所采用的技术方案是：
8.一种典型场景下谐波电压电流测量与特征分析方法，包括一种典型场景下谐波电压电流测量与特征分析方法，包括以下步骤：
9.s1：选择若干典型场景，研究确定谐波流进行现场测量，基于组合余弦窗函数参数优化的电力谐波分析方法提取pcc处的电压电流谐波分量，利测量方案；
10.s2：对电网公共连接点谐波电压电用窗函数对电压电流谐波分量进行优化，并采用粒子群算法对窗函数进行优化，以旁瓣峰值电平和旁瓣渐近衰减速率作为窗函数参数优化的目标；
11.s3：分析各频次谐波电压电流的主要特征及演变规律。
12.进一步地，所述s2中采用粒子群算法对窗函数进行优化的具体步骤如下：
13.s2-1：对窗函数的各项数参数进行随机初始化，并初始化各项系数对应的速度和位置；
14.s2-2：以旁瓣峰值电平和旁瓣渐近衰减速率为窗函数参数优化的目标，根据公式
(1)计算适应函数值f
fit
，并且得到粒子的历史最优位置和群体的全局最优位置；
15.其中，φ为旁瓣峰值电平，λ为旁瓣渐近衰减速率，α为旁瓣峰值电平权重，β为旁瓣渐近衰减速率权重，θ为惩罚因子，α β＝1。惩罚因子一般情况下为1，在寻优过程中出现λ《0、φ《0等异常情况时为0；
16.s2-3：粒子通过跟踪两个极值pbest和gbest来更新，根据粒子自身的历史最优位置和全局的最优位置，运用公式(2)(3)更新每个粒子的速度和位置；vi＝ω
×
vi c1×
rand()
×
(pbest
i-xi) c2×
rand()
×
(gbest
i-xi)(3)x'i＝xi vi(4)
17.其中，vi为第i个粒子的速度；rand()为介于(0,1)之间的随机数；xi为粒子的当前位置；x'i为粒子更新后的位置，c1，c2分别为两个极值的权重值，ω为惯性因子；pbesti为自身的历史最优位置，gbesti为全局的最优位置；
18.s2-4：运用公式(1)对更新速度和位置后的粒子计算适应度，进而更新粒子的历史最优位置和全局的最优位置；
19.s2-5：检查粒子群是否满足如式(5)所示的约束条件，若满足约束条件，则输出全局最优结果并结束程序，若不满足约束条件，则转向步骤s2-3继续执行；
20.其中，m为窗函数；m为窗函数项数，bm为初始粒子群为窗函数各项系数，每个系数包括对应的位置和速度。
21.进一步地，所述s2-3中c1＝0.7，c2＝0.3。
22.进一步地，所述s2-3中惯性因子的计算公式为ω
(t)
＝(ω
ini-ω
end
)(g
k-g)/gk ω
end
(5)
23.进一步地，其中g为迭代次数，gk为最大迭代次数，ω
ini
为初始惯性权值，ω
end
为迭代至最大进化代数时的惯性权值；进一步地，ω
ini
＝0.9，ω
end
＝0.4。
24.进一步地，所述s1的具体操作方法为：
25.s1-1：充分调研并选择适合配置谐波测量设备的典型场景的110/220kv接入处，结合其二次接线图确定测点；
26.s1-2：考虑并网点处负荷动态行为，研究制定细致的谐波测量方案。
27.进一步地，所述s3的具体操作方法为：
28.s3-1：利用概率统计方法提取各频次谐波电压电流的时频域分布特征；
29.s3-2：对负荷动态特性划分；
30.s3-3：研究公共连接点处谐波电压电流的主要特征随负荷动态特性的演变规律，基于s3-1和s3-2的结果研究公共连接点处谐波电压电流的主要特征在24小时周期内的演变规律。
31.进一步地，所述典型场景包括钢铁厂10kv/35kv系统、光伏电站并网点和高速铁路牵引供电系统。
32.本发明的有益效果为：
33.1、本发明考虑负荷动态特性进行工况划分制定细致谐波测量方案的策略，提高电力系统谐波问题分析能力，典型场景下的谐波特性分析及演变规律为谐波源模型建立提供了数据基础，为提高电网谐波污染管理水平奠定坚实基础。
34.2、基于组合余弦窗函数参数优化实现的电力谐波分析方法，通过半实物仿真平台验证其效果，充分考虑区域电网负荷动态特性，进一步优化参数提供谐波分析方法的准确性和计算速度；以旁瓣峰值电平和旁瓣渐近衰减速率为窗函数参数优化的目标，旁瓣峰值电平小且旁瓣渐近衰减速率大的窗函数能有效地抑制频谱泄漏。
附图说明
35.图1为本发明实施例中典型场景下谐波电压电流测量与特征分析方法的整体技术线路示意图。
36.图2为本发明实施例中粒子群算法优化的流程示意图。
具体实施方式
37.下面结合附图和实施例对本发明的实施方式作进一步描述。需要说明的是，实施例并不对本发明要求保护的范围构成限制。
38.实施例1
39.如附图1所示，一种典型场景下谐波电压电流测量与特征分析方法，包括以下步骤：
40.s1：选择若干典型场景，研究确定谐波测量方案；
41.本发明的典型场景拟选择钢铁厂10kv/35kv系统、光伏电站并网点或高速铁路牵引系统 110/220kv并网点。考虑负荷动态行为，研究制定细致的谐波测量方案；研究根据钢铁厂日生产计划划分典型工况，确定测量方案；研究高速铁路牵引负荷运行周期性规律，充分考虑牵引负荷处于牵引、再生制动、惰行等工况，确定测量方案。实施步骤如下：
42.s1-1：充分调研并选择适合配置谐波测量设备的钢铁厂10kv/35kv系统、光伏电站并网点、高速铁路牵引供电系统110/220kv接入处，结合其二次接线图确定测点；
43.s1-2：考虑并网点处负荷动态行为，研究制定细致的谐波测量方案，研究根据钢铁厂日生产计划划分典型工况，确定测量方案；研究高速铁路牵引负荷运行周期性规律，充分考虑牵引负荷处于牵引、再生制动、惰行等工况，确定测量方案。
44.s2：对电网公共连接点谐波电压电流进行现场测量；
45.研究基于组合余弦窗函数参数优化实现的电力谐波分析方法，通过半实物仿真平台验证其效果，充分考虑区域电网负荷动态特性，进一步优化参数提供谐波分析方法的准确性和计算速度。其中，对于如何通过选择合适的窗函数及其参数来提高窗函数频谱特性，进而提高检测精度，窗函数的频域特性主要包括主瓣宽度、旁瓣峰值、旁瓣衰减速率以及主瓣旁瓣能量比等，拟采用的解决方案是针对组合余弦窗函数等常用窗函数，以旁瓣峰值电平和旁瓣渐近衰减速率作为窗函数参数优化的目标，对应合适的约束条件，采用粒子群优化算法进行求解，对常用窗函数参数进行优化，然后对比分析其谐波分析精度，进而得到频谱特性最佳的窗函数及其优化参数。操作步骤如下：
46.首先，对于选择的钢铁厂10kv/35kv系统、光伏电站并网点、高速铁路牵引供电系
统110/220kv接入处，结合其二次接线图确定具体量测点；测试前对测试仪器和探头进行可靠性确认并正确设置，测试前及测试中定期对不同测点的测试系统进行时钟校准；基于组合余弦窗参数优化的电力谐波分析方法提取pcc处的电压电流谐波分量。
47.采用粒子群算法优化如下：
48.s2-1：随机初始化。对常用组合余弦窗函数的各项数参数进行随机初始化，并初始化各项数对应的速度和位置；
49.s2-2：根据公式(2)计算适应函数值，并且得到粒子的历史最优位置和群体的全局最优位置。
50.s2-3：粒子的速度和位置更新。根据粒子自身的历史最优位置和全局的最优位置，运用公式(3)(4)更新每个粒子的速度和位置。
51.s2-4：运用公式(2)对更新速度和位置后的粒子计算适应度，进而更新粒子的历史最优位置和全局的最优位置。
52.s2-5：如果满足约束条件，则输出全局最优结果并结束程序，否则转向步骤2继续执行。
53.(1)对于组合余弦窗函数参数优化，粒子群算法的初始粒子群为窗函数各项系数bm,每个系数包括对应的位置和速度，因此优化的约束条件为：
[0054][0055]
其中，m为窗函数；m为窗函数项数
[0056]
(2)确定微粒的适应度。本发明以旁瓣峰值电平和旁瓣渐近衰减速率为窗函数参数优化的目标，旁瓣峰值电平小且旁瓣渐近衰减速率大的窗函数能有效地抑制频谱泄漏。因此适应度函数f
fit
可构造为：
[0057][0058]
公式中：φ为旁瓣峰值电平，λ为旁瓣渐近衰减速率，α为旁瓣峰值电平权重，β为旁瓣渐近衰减速率权重，θ为惩罚因子，α β＝1。惩罚因子一般情况下为1，在寻优过程中出现λ《0、α《0等异常情况时为0。
[0059]
(3)粒子通过跟踪两个极值pbest和gbest来更新，粒子的速度和位置更新公式如下：
[0060]
vi＝ω
×
vi c1×
rand()
×
(pbest
i-xi) c2×
rand()
×
(gbest
i-xi)(3)
[0061]
x'i＝xi vi(4)
[0062]
其中vi为粒子的速度；rand()为介于(0,1)之间的随机数；xi为粒子的当前位置；x'i为粒子更新后的位置，c1，c2分别为两个极值的权重值，ω为惯性因子；c1＝0.7，c2＝0.3。
[0063]
惯性因子ω的计算公式如下：
[0064]
ω
(t)
＝(ω
ini-ω
end
)(g
k-g)/gk ω
end
(5)
[0065]
其中g为迭代次数，gk为最大迭代次数，ω
ini
为初始惯性权值，ω
end
为迭代至最大进化代数时的惯性权值；进一步地，ω
ini
＝0.9，ω
end
＝0.4。
[0066]
根据上述优化的步骤通过对粒子群算法的参数不断调试，可对任意项数的组合余
弦窗函数进行优化。
[0067]
s3：分析各频次谐波电压电流的主要特征及演变规律。
[0068]
研究各频次谐波电压电流的时频域分布特征，利用概率统计方法提出电网谐波的统计性描述参数，研究典型场景下负荷动态特性划分，如高速铁路牵引供电系统，以行车运行图为基础研究牵引变电所全天 24小时行车运行时空分布，充分考虑牵引负荷处于牵引、再生制动、惰行等工况及空载情况。研究公共连接点处谐波电压电流的主要特征在24小时周期内随行车运行及其工况变化的演变规律。步骤如下
[0069]
s3-1：利用概率统计方法提取各频次谐波电压电流的时频域分布特征；
[0070]
s3-2：对负荷动态特性划分，如高速铁路牵引供电系统，以行车运行图为基础，充分考虑牵引负荷处于牵引、再生制动、惰行等工况以及空载情况，对牵引变电所全天24小时行车运行工况进行划分；
[0071]
s3-3：研究公共连接点处谐波电压电流的主要特征随负荷动态特性的演变规律如针对高速铁路牵引供电系统，基于步骤一和步骤二结果研究公共连接点处谐波电压电流的主要特征在24小时周期内随行车运行及其工况变化的演变规律。
[0072]
按照如上步骤针对若干典型场景，提出充分考虑其负荷时空动态特性的谐波测量方案，并针对性地提出适用于各场景谐波特性的改进快速傅里叶变换分析方法，准确提取谐波分量，提取各频次谐波电压电流的时频域统计分布特性，揭示公共连接点处各频次谐波电压电流的主要特征及演变规律。
[0073]
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。