本发明涉及关口计量设备状态评估
技术领域：
，尤其是一种用于电能误差分解的时频分析方法。
背景技术：
：随着电力行业市场化改革的深入进行，发电上网、跨区输电、跨省输电及省级供电等关口电量交换日益增加，电力企业开始更加关注自身经济利益的维护。如何确保电能计量装置的准确性，维持电力市场公平有序的运转是电力市场当前研究的一个重要课题。关口电能计量装置运行管理主要包括现场检验及周期轮换，DL/T448-2000《电能计量装置技术管理规程》规定I类、II类和III类电能表的现场检验周期分别为至少3个月、6个月和1年。目前，这种现场检验方式已难以适应电能计量技术和电力公司精细化管理的要求，主要表现为：1)电网规模正在不断扩大，交易电量和电能计量装置也越来越多，在有限的人力条件下实现规范化的技术管理非常困难；2)传统人工现场检验的工作效率低，不能对装置进行实时监测和故障预警，尤其难以及时发现两次现场检验之间出现的计量问题，对电量追补工作带来了很大的困难。因此，须亟需在线监测电能计量装置的运行状态，对其误差水平进行实时测量和趋势预测，从而克服计量装置现场检验模式的弊端，提高其运行管理水平。特别是电能计量装置误差的变化趋势预测是其运行状态评估的重要环节，对于实现电能计量装置超差预警及优化检验周期具有重要的实用价值。电能计量装置的误差水平由电能表、电压和电流互感器、二次回路4个部分的误差综合而成，由于各部分误差来源的因素多样性和内在关联性，使得电能计量误差时间序列呈现显著的非平稳变化特征。现有研究仅对计量装置的误差进行在线监测，未给出计量误差趋势预测的具体方法。计量误差预测是其运行状态预报和故障预警的基础，考虑计量误差时间序列的随机性和不确定性，可采用目前普遍用于新能源发电预测的持续时间法、自回归滑动平均(AutoRegressiveMovingAverage，ARMA)法、卡尔曼滤波法及其他人工智能方法进行预测，也有研究提出将部分模型组合的预测方法，但是前述方法对于计量误差预测的适应性还需进一步研究。电能计量误差包含多种时间和频率尺度的分量，直接对该误差进行预测难以准确获得未来误差趋势的估计值。因此，须采用适当的时频分析法将电能误差分解成局部时间和频率尺度的序列分量，进而对各分量分别进行预测，以提高电能计量误差水平预测的准确性。目前主要采用小波变换、经验模式分解(EmpiricalModeDecomposition，EMD)和局部均值分解，其中小波变换利用固定的基函数分解时间序列，对其局部特征的表征能力有限，也无法进行自适应分解；而经验模式分解法则能把非平稳时间序列自适应地分解为多个本征模式分量，但该方法仍存在端点效应、模式混叠、过/欠包络等不足；为了解决过/欠包络问题和减缓端点效应，局部均值分解法将复杂时间序列分解成多个由调幅和调频信号相乘的乘积函数，但该方法的最优平滑步长选择困难，且计算量大无法用于在线分析。技术实现要素：本发明所要解决的技术问题是提供一种能够提高在线分解效率且能够减小计算耗时，同时也能够很好地保留原始信号的在时间尺度上的信息的用于电能误差分解的时频分析方法。本发明解决其技术问题所采用的技术方案为：该用于电能误差分解的时频分析方法，包括以下步骤：a、获取从0至t时刻的计量装置误差时间序列x(t)，设时间序列x(t)包含M个极值点Xi，对应的出现时刻为Ti(i＝1,2,…,M)，通过如下公式提取Ti至Ti+1时间段内基线信号L(t)，L(t)=Li+(Li+1-LiXi+1-Xi)[x(t)-Xi],t∈[Ti,Ti+1],Li+1=α[Xi+(Ti+1-TiTi+2-Ti)(Xi+2-Xi)]+(1-α)Xi+1,]]>式中：Li+1为Ti至Ti+1时间段内的基线提取算子，α∈[0,1]为常系数，若基线信号提取的极值点出现时刻T0＝0且仅提取至Ti-2至Ti-1时间段；b、计算Ti至Ti+1时间段内旋转分量提取算子Hi+1，Hi+1＝(1-Li+1)；c、利用Ti至Ti+1时间段内基线提取算子Li+1和旋转分量提取算子Hi+1将Ti至Ti+1时间段内时间序列x(t)依次分解为：x(t)=Li+1x(t)+Hi+1x(t)=Li+1(Li+1+Hi+1)x(t)+Hi+1x(t)=[Hi+1(1+Li+1)+Li+12]x(t)=[Hi+1(1+Li+1+Li+12)+Li+13]x(t)=(Hi+1Σj=0n-1Li+1j+Li+1n)x(t)]]>Ti至Ti+1时间段内时间序列x(t)被分解为n个旋转分量和1个趋势分量即为：xRn(t)＝Hi+1Li+1n-1x(t),xT(t)＝Li+1nx(t)当趋势分量xT(t)为单调或者常函数时，则停止分解输出Ti至Ti+1时间段内时间序列x(t)的各个分量。d、重复步骤A、B、C逐步分解出x(t)的旋转分量xRn(t)和趋势分量xT(t)。进一步的是，在步骤b中，所述α取值为0.5。本发明的有益效果：该用于电能误差分解的时频分析方法采用相邻极值点间分段线性化的方式有效提高了在线分解效率，计算耗时明显减少，同时也很好地保留了原始信号的在时间尺度上的信息。附图说明图1表示采用本发明所述ITD以及与现有的EMD和LMD三种方法分解时前3个(高频)分量与原始时间序列的相关系数和计算耗时(分解10次的平均耗时)。具体实施方式本发明所述的用于电能误差分解的时频分析方法，包括以下步骤：a、获取从0至t时刻的计量装置误差时间序列x(t)，设时间序列x(t)包含M个极值点Xi，对应的出现时刻为Ti(i＝1,2,…,M)，通过如下公式提取Ti至Ti+1时间段内基线信号L(t)，L(t)=Li+(Li+1-LiXi+1-Xi)[x(t)-Xi],t∈[Ti,Ti+1],Li+1=α[Xi+(Ti+1-TiTi+2-Ti)(Xi+2-Xi)]+(1-α)Xi+1,]]>式中：Li+1为Ti至Ti+1时间段内的基线提取算子，α∈[0,1]为常系数，若基线信号提取的极值点出现时刻T0＝0且仅提取至Ti-2至Ti-1时间段；b、计算Ti至Ti+1时间段内旋转分量提取算子Hi+1，Hi+1＝(1-Li+1)；c、利用Ti至Ti+1时间段内基线提取算子Li+1和旋转分量提取算子Hi+1将Ti至Ti+1时间段内时间序列x(t)依次分解为：x(t)=Li+1x(t)+Hi+1x(t)=Li+1(Li+1+Hi+1)x(t)+Hi+1x(t)=[Hi+1(1+Li+1)+Li+12]x(t)=[Hi+1(1+Li+1+Li+12)+Li+13]x(t)=(Hi+1Σj=0n-1Li+1j+Li+1n)x(t)]]>Ti至Ti+1时间段内时间序列x(t)被分解为n个旋转分量和1个趋势分量即为：xRn(t)＝Hi+1Li+1n-1x(t),xT(t)＝Li+1nx(t)当趋势分量xT(t)为单调或者常函数时，则停止分解输出Ti至Ti+1时间段内时间序列x(t)的各个分量。d、重复步骤A、B、C逐步分解出x(t)的旋转分量xRn(t)和趋势分量xT(t)。该分解方法采用相邻极值点间分段线性化的方式有效提高了在线分解效率，计算耗时明显减少，同时也很好地保留了原始信号的在时间尺度上的信息。进一步的是，在步骤b中，所述α取值优选为0.5。图1表示采用本发明所述ITD以及与现有的EMD和LMD三种方法分解时前3个(高频)分量与原始时间序列的相关系数和计算耗时(分解10次的平均耗时)，其中第1至第3个分量相关系数依次减小，前两个分量与原始信号的相关性较强，采用本发明所述ITD得到的第1个相关系数和其余两种方法基本一致，但其第2个相关系数高于EMD和LMD法。此外，ITD法的计算耗时明显低于其他方法，可见ITD在保留原始时间序列信息和在线分解方面均具有一定优势。当前第1页1 2 3