一种基于校验的同步相量自适应计算方法与流程

本发明涉及电力系统自动测量
技术领域：
，是一种基于校验的同步相量自适应计算方法，特别涉及一种高精度和快速响应的同步相量测量方法。
背景技术：
：近年来，以同步相量测量装置(phasormeasurementunit，PMU)为基础的广域测量系统在电力系统动态过程监视、在线辨识、安全稳定分析以及广域控制等领域中得到广泛的应用。随着广域测量系统应用研究的不断深入，PMU装置对同步相量测量的要求越来越高，其相量算法的准确性和快速性将直接影响到相关应用功能的可靠性。已有基于DFT(DiscreteFourierTransform，离散傅里叶变换)的修正算法，相对于传统DFT算法，较大的提高了计算精度，但对突变等动态过程响应速度有限，且在幅值时刻变化时难以满足精度要求。而已有基于频域动态模型的算法，利用同一数据窗不同频点滤波器的响应来修正DFT的估计结果，提高了对突变等动态过程的响应速度，但此算法存在噪声和谐波抑制能力弱，且在故障过程中精度低的问题。技术实现要素：本发明的目的在于提供一种基于校验的同步相量自适应计算方法，该方法对稳态和动态情况，分别设计了稳态算法和动态算法两种不同计算模式，通过反推校验实现二者间的自适应切换，以同时兼顾稳态情况的计算精度及暂态情况的动态性能要求。本发明方法可在不同情况下，通过校验实现计算模式自适应切换，可以兼顾快速响应和高精度的要求。本发明采用的技术方案是：本发明的基于校验的同步相量自适应计算方法，包括以下步骤：(1)初始化，确定每周波采样点数N，稳态和动态算法的电力信号模型中幅值和相角的阶数K，稳态算法中两数据窗的间隔l点；(2)信号建模，该方法采用复信号P(t)表示电力信号的动态相量为P(t)＝a(t)ejθ(t)，电力信号x(t)表示为：式中：a(t)和θ(t)分别表示电力信号幅值和相角的多项式；f0为额定频率；稳态算法模型假设计算周期内幅值恒定，频率偏差恒定，即a(t)＝a，θ(t)＝θ0+θ1t；动态算法中为更好反映信号的动态特征，假设相量模型中幅值和相角均为K阶模型，即(3)通过两个数据窗数据的DFT(DiscreteFourierTransform，离散傅里叶变换)结果，求解电力信号的同步相量，模型参数：校正后相角：校正后幅值：式中：θc和θcl为校正前两数据窗分别计算角度，ac为校正前DFT计算的幅值，(4)通过比较反推值和实测值，校验计算结果是否正确，如果校验通过则结束，若未通过则说明当前可能处于突变等动态过程，进入步骤(5)；(5)假设电力信号为K阶模型，采用已有动态算法求解同步相量，如通过基于同一数据窗各谐波含量计算同步相量的方法；(6)同步骤(4)，校验动态算法计算结果是否正确，若校验通过则使用动态算法的结果作为最终结果，如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况，为避免波动较大仍使用稳态算法结果。上述步骤(4)的具体地校验过程为：通过计算得到的相量，反推各采样点的理论计算值，并与实际测得的点进行比较，若总偏差较小，则校验通过。本发明针对稳态算法，由校正后θ0、θ1和a，反推各抽取点理论计算值，抽取点的密度根据实际情况确定。本发明抽取点的理论计算值为：xc(ni)=acos(θ0+niθ1Nf0+2πniN)]]>式中：ni(i＝1,2,···,L)为所抽取点，且1≤ni≤N；xc(ni)为所抽取点的理论计算值；理论计算值与实测值之间的总偏差为：ϵ=Σi=1L|xc(ni)-xm(ni)|]]>式中：xm(ni)对应采样点实测值。本发明由于采取以上技术方案，与现有技术相比，具有以下优点：本发明在同时考虑PMU算法快速性和准确性要求下，提出一种基于校验的同步相量自适应计算方法，在稳态和动态情况下自适应采用不同的计算模式。本发明测量方法对稳态和动态情况，分别设计了稳态算法和动态算法两种不同计算模式，通过反推校验实现二者间的自适应切换，从而同时兼顾稳态情况的计算精度及暂态情况的动态性能要求。本发明针对已有算法无法同时兼顾稳态情况的计算精度及暂态情况的动态性能问题，本发明方法可在不同情况下，通过校验实现计算模式自适应切换，可以兼顾快速响应和高精度的要求。附图说明图1为本发明方法的流程示意图。图2为本发明方法实施例中的算法流程示意图。图3为本发明方法的仿真测试中阶跃响应示意图。具体实施方式以下结合附图对本发明进行详细的描述。本发明提出来的基于校验的同步相量自适应计算方法，可以采用多种硬件方案来实现，在此不再赘述。本发明所提出的测量算法流程如图1所示，其中PMU算法的流程如图2所示。测量方法包括以下步骤：(1)初始化，确定每周波采样点数N，稳态和动态算法的电力信号模型中幅值和相角的阶数K，稳态算法中两数据窗的间隔l点。(2)信号建模，该方法采用复信号P(t)表示电力信号的动态相量为P(t)＝a(t)ejθ(t)，电力信号x(t)可以表示为：式中：a(t)和θ(t)分别表示电力信号幅值和相角的多项式；f0为额定频率。稳态算法模型假设计算周期内幅值恒定，频率偏差恒定，即a(t)＝a，θ(t)＝θ0+θ1t。动态算法中为更好反映信号的动态特征，假设相量模型中幅值和相角均为K阶模型，即(3)通过两个数据窗数据的DFT结果，求解电力信号的同步相量。模型参数：校正后相角：校正后幅值：式中：θc和θcl为校正前两数据窗分别计算角度，ac为校正前DFT计算的幅值，(4)通过比较反推值和实测值，校验计算结果是否正确，如果校验通过则结束，若未通过则说明当前可能处于突变等动态过程，进入(5)；具体地校验过程为：通过计算得到的相量，反推各采样点的理论计算值，并与实际测得的点进行比较，若总偏差较小，则校验通过。针对稳态算法，由校正后θ0、θ1和a，反推各抽取点理论计算值，抽取点的密度根据实际情况确定。抽取点的理论计算值为：xc(ni)=acos(θ0+niθ1Nf0+2πniN)]]>式中：ni(i＝1,2,···,L)为所抽取点，且1≤ni≤N；xc(ni)为所抽取点的理论计算值。理论计算值与实测值之间的总偏差为：ϵ=Σi=1L|xc(ni)-xm(ni)|]]>式中：xm(ni)对应采样点实测值。(5)假设电力信号为K阶模型，采用已有动态算法求解同步相量，如通过基于同一数据窗各谐波含量计算同步相量的方法。(6)同步骤(4)，校验动态算法计算结果是否正确，若校验通过则使用动态算法的结果作为最终结果。如果校验未通过则可能故障发生在本数据窗内等动态算法同样无法处理的情况，为避免波动较大仍使用稳态算法结果。为进一步说明本发明方法，下面对本发明方法进行仿真测试。仿真过程中算法的采样率为6400Hz等间隔采样。1、频率偏差测试为了测试相量测量算法在频率偏离额定频率时的性能，国标《电力系统同步相量测量装置检测规范》规定频率测量范围是45Hz至55Hz，在基波频率偏离额定值5Hz时，电压、电流幅值测量误差改变量应小于额定频率时测量误差极限值的100％，相角测量误差改变量应不大于1°。表1为本发明方法在频率偏离额定频率5Hz时的测试结果。可以看出，本发明方法的量测精度高于标准要求。表1频率偏差为5Hz时相量测试结果最大误差平均绝对误差均方根误差角度误差-0.02°0.0174°0.0175°幅值误差0.014％0.0092％0.0097％2、为测试算法对突变等动态过程的响应性能，参考国标《电力系统同步相量测量装置检测规范》，施加90°相角阶跃信号：x(t)=acos(2πf0t+π/6),t<40msacos(2πf0t+π/6+π/2),t≥40ms]]>本发明方法测试结果如图3所示，可以看出，本发明方法的阶跃响应时间为20ms明显快于国标规定的30ms。当前第1页1 2 3