电力谐波分析方法技术

本发明专利技术公开了一种电力谐波分析方法，包括获取输入电压实际周期采样点数；计算每个采样点所需平移量；计算该插值点所需总平移量；将总平移量的整数部分作为偏移量获取对应的实际采样点值，根据得到的三个采样值和总平移量的小数部分进行二次插值得到该插值点对应的插值数据；重复上述步骤直至已处理点数达到总点数，从而得到全新的插值序列；再存储的电压和电流的插值序列进行分析和计算电力谐波参数。本发明专利技术方法有效抑制了非同步采样时FFT的频谱泄露引起的谐波分析误差；省去了硬件同步采样中复杂的锁相环电路，也无需对ADC采样间隔进行微调，因而适用面极广；有效减小插值对各次谐波幅值的衰减，大大提高谐波分析准确度和运算效率。

【技术实现步骤摘要】
电力谐波分析方法
本专利技术属于电气自动化领域，具体涉及一种电力谐波分析方法。
技术介绍
随着电力电子技术的飞速发展，如光伏电源、风电、电弧炉、电气化铁路和轧钢机等负荷的广泛使用，非线性负荷在电网的比重逐渐增加。目前非线性负荷在电网中造成的谐波污染问题日益凸显，电力管理部门对谐波表计量准确性的关注度也日益提高。国内于2014年发布实施了谐波表标准《GB/T17215.302-2013静止式谐波有功电能表》，对谐波表进行了规范统一，以使其在电网中更有效、更合理地发挥作用。目前应用最广泛的谐波分析方法是FFT算法，但众所周知的是欲获得准确可靠的分析结果，必须解决非同步采样引起的FFT频谱泄漏问题，非同步采样表现为采样频率与电网基波频率的不同步。同时FFT算法要求信号的一周期采样点数为2的整数幂。一般情况下，电网基波频率波动范围很小，通常在±0.5Hz以内；但在某些大功率的非线性负荷下，电网频率波动范围较大，可达几赫兹。在谐波表国家标准GB/T17215.302-2013中，要求频率改变量在±2％以内，高精度的1级谐波表的谐波电能计量误差改变量小于±0.5％，所以如果采用FFT算法实现谐波电能计量必须考虑频率波动带来的影响。目前减小非同步采样误差的方法主要有加窗插值算法和同步采样技术两大类。在时域加余弦窗可以有效地减少频谱泄漏，在频域对FFT结果进行频谱插值可以减小栅栏效应引起的误差，但加窗和频谱插值处理的数据运算量大，并且包括解高次线性方程、除法运算、谱线峰值搜索等处理，对CPU性能要求高，不适用于电表平台。同步采样技术有硬件同步采样和软件同步采样两种。硬件同步采样是利用锁相环实现采样频率对基波频率的跟踪，但需要对输入信号进行低通滤波滤除工频以外的直流和谐波成分，采样频率较高时锁相环分频数较大，锁相倍频电路设计难度较大，而且该技术只能用在基于逐次逼近(SAR)并且可由硬件触发采样的ADC上，某些其他类型的ADC如Σ-ΔADC无法使用该技术，其采样间隔无法微调。软件同步采样是根据测量得到的输入信号频率，通过定时器或其他可编程延时模块对ADC采样间隔进行微调，从而实现频率跟踪。该算法也只能用在SARADC上，高分辨率、高集成度且低成本的Σ-ΔADC启动时间过长，无法使用该技术，而高精度电能计量系统中的模拟采样通常需要采用Σ-ΔADC。因此，Σ-ΔADC的非同步采样及其在电能计量应用中采用FFT进行谐波分析是目前急需解决的问题。
技术实现思路
本专利技术的目的在于提供一种高精度的电力谐波分析方法。本专利技术提供的这种电力谐波分析方法，包括如下步骤：S1.获取输入电压实际的周期采样点数；S2.根据步骤S1中得到的实际的周期采样点数和事先设定的分析点数计算每个采样点所需的平移量；S3.根据插值点的序号和每个点的平移量计算该插值点所需的总平移量，同时获取该总平移量的整数部分和小数部分；S4.将步骤S3中得到的总平移量的整数部分作为偏移量获取对应的实际采样点值，并根据得到的三个采样值和步骤S3中得到的总平移量的小数部分进行二次插值，得到该插值点对应的插值数据；S5.重复步骤S3～S4直至已处理点数达到总点数，从而得到全新的插值序列；再利用步骤S4中存储的电压和电流的插值序列进行分析，根据分析结果计算电力谐波参数。步骤S2中所述的计算每个采样点所需要的平移量，具体为采用如下公式进行计算：Δ＝(M-N)/N式中，Δ为每个采样点的平移量；M为实际的周期采样点数；N为事先设定的分析点数。步骤S3中所述的计算当前插值点所需的总平移量并获取其整数部分和小数部分，具体为采用如下公式进行计算：Δn＝n·Δp＝[Δn]＝[n·Δ]q＝Δn-p＝n·Δ-p式中，n为插值点的序号，取值范围为0～N-1；Δn为当前插值点所需的总平移量；p为总平移量的整数部分，[Δn]表示对Δn取整；q为总平移量的小数部分。步骤S4中所述的将步骤S3中得到的总平移量的整数部分作为偏移量获取对应的实际采样点值，具体为根据当前插值点的序号n和步骤S3中的总平移量整数部分p从电压和电流采样值序列中查找输入数据块起始点之后n+p、n+p+1和n+p+2点的采样值。步骤S4中所述的二次插值为Lagrange二次插值。步骤S4中所述的根据得到的三个采样值和步骤S3中得到的总平移量的小数部分进行二次插值得到该插值点对应的插值数据，具体为采用如下算式进行计算：u'n＝(q-1)(q-2)/2·un+p-q(q-2)·un+p+1+q(q-1)/2·un+p+2i'n＝(q-1)(q-2)/2·in+p-q(q-2)·in+p+1+q(q-1)/2·in+p+2式中，n为插值点的序号，取值范围为0～N-1；p为总平移量的整数部分；q为总平移量的小数部分；un+p、un+p+1和un+p+2及in+p、in+p+1和in+p+2分别为输入的电压和电流数据块起始点之后第n+p、n+p+1和n+p+2点的采样点值；u′n和i′n为插值后的第n点数据。步骤S5中所述的对存储的电压和电流的插值序列进行分析，具体为对存储的电压和电流的插值序列进行FFT分析。步骤S5中所述的电力谐波参数包括各次谐波电压、各次谐波电流和总谐波有功功率。本专利技术提供的这种电力谐波分析方法，通过使用重采样技术，将一周期采样点数并非2的幂数的原始采样序列转换为FFT算法所需的2的幂数个点，即使在非同步采样下也可得到完整周期并且适用于FFT的数据序列，这有效抑制了非同步采样时FFT的频谱泄露引起的谐波分析误差；本专利技术方法采用纯软件算法实现，省去了硬件同步采样中复杂的锁相环电路，也无需对ADC采样间隔进行微调，因而不仅适用于SARADC也便于在难以修改ADC周期性采样间隔的Σ-ΔADC的电能计量应用中正常使用FFT算法；此外，由于采用Lagrange二次插值法实现重采样，与线性插值法相比，有效减小插值对各次谐波幅值的衰减，大大提高谐波分析准确度，而且与加窗插值FFT算法相比具有更高的运算效率。附图说明图1为本专利技术方法的基本原理框图。图2为本专利技术的方法流程图。图3为本专利技术方法中的Lagrange二次插值示例示意图。图4为本专利技术方法中采用线性插值时的精度结果示意图。图5为本专利技术方法中采用Lagrange二次插值时的精度结果示意图。具体实施方式如图1所示为本专利技术的基本原理框图：本专利技术主要基于Lagrange插值和FFT重采样技术。首先利用Σ-ΔADC采集电压和电流，得到离散的电压采样值和电流采样值序列，并计算电压的周期值；然后根据电压的周期对电压和电流进行插值重采样；接着通过FFT算法对插值后的电压和电流数据序列进行谐波分析；最后根据FFT分析结果计算谐波电压、谐波电流、总谐波有功等谐波计量指标。如图2所示为本专利技术的方法流程图：本专利技术提供的这种电力谐波分析方法，包括如下步骤：S1.获取输入电压实际的周期采样点数；S2.根据步骤S1中得到的实际的周期采样点数和事先设定的分析点数计算每个采样点所需的平移量，具体为采用如下公式进行计算：Δ＝(M-N)/N式中，Δ为每个采样点的平移量；M为实际的周期采样点数；N为事先设定的分析点数；S3.根据插值点的序号和每个点的平移量计算该插值点所需的总平移量，同时获取该总平移量的整数部分和小数部分，本文档来自技高网...

【技术保护点】
一种电力谐波分析方法，包括如下步骤：S1.获取输入电压实际的周期采样点数；S2.根据步骤S1中得到的实际的周期采样点数和事先设定的分析点数计算每个采样点所需的平移量；S3.根据插值点的序号和每个点的平移量计算该插值点所需的总平移量，同时获取该总平移量的整数部分和小数部分；S4.将步骤S3中得到的总平移量的整数部分作为偏移量获取对应的实际采样点值，并根据得到的三个采样值和步骤S3中得到的总平移量的小数部分进行二次插值，得到该插值点对应的插值数据；S5.重复步骤S3～S4直至已处理点数达到总点数，从而得到全新的插值序列；再利用步骤S4中存储的电压和电流的插值序列进行分析，根据分析结果计算电力谐波参数。

【技术特征摘要】
1.一种电力谐波分析方法，包括如下步骤：S1.获取输入电压实际的周期采样点数；S2.根据步骤S1中得到的实际的周期采样点数和事先设定的分析点数计算每个采样点所需的平移量；S3.根据插值点的序号和每个点的平移量计算该插值点所需的总平移量，同时获取该总平移量的整数部分和小数部分；S4.将步骤S3中得到的总平移量的整数部分作为偏移量获取对应的实际采样点值，并根据得到的三个采样值和步骤S3中得到的总平移量的小数部分进行二次插值，得到该插值点对应的插值数据；S5.重复步骤S3～S4直至已处理点数达到总点数，从而得到全新的插值序列；再利用步骤S4中存储的电压和电流的插值序列进行分析，根据分析结果计算电力谐波参数。2.根据权利要求1所述的电力谐波分析方法，其特征在于步骤S2中所述的计算每个采样点所需要的平移量，具体为采用如下公式进行计算：Δ＝(M-N)/N式中，Δ为每个采样点的平移量；M为实际的周期采样点数；N为事先设定的分析点数。3.根据权利要求1所述的电力谐波分析方法，其特征在于步骤S3中所述的计算当前插值点所需的总平移量并获取其整数部分和小数部分，具体为采用如下公式进行计算：Δn＝n·Δp＝[Δn]＝[n·Δ]q＝Δn-p＝n·Δ-p式中，n为插值点的序号，取值范围为0～N-1；Δn为当前插值点所需的总平移量；p为总平移量的整数部分，[Δn]表示对Δn取整；q为总平移量的小数部分。4.根据权利要求1所述的电力谐波分析方法，其特征在于步骤S4中所述的将步骤S3中得到的总...

【专利技术属性】
技术研发人员：汪龙峰，周杰文，任智仁，
申请(专利权)人：威胜集团有限公司，
类型：发明
国别省市：湖南,43