一种无功功率自适应计算方法
【技术领域】
[0001 ]本发明属于电网电能计量技术领域，特别涉及一种无功功率自适应计算方法。
【背景技术】
[0002] 随着各种节能电器以及类似快速充电站这样的脉冲电能负载的产生，电力系统参 数越来越复杂，给电能的准确计量带来一定困扰。在保证电力系统稳定性方面，通过调节有 功功率和无功功率得以实现控制策略，因此，研究系统参数随机可调情况下的功率测量方 法对电力系统供电可靠性和经济性具均有指导作用。
[0003] 在实际电力系统电能计量时，通常需要把电压信号或电流信号实时分解成基波分 量和多个频率已知的谐波分量，并计算基波和各次谐波的幅值、相位，虽然取得了一些效 果，但同时也存在一定的不足之处，例如：
[0004] 基于傅里叶变换的分析方法，离散傅里叶变换能分析出非整数次谐波，但由于计 算量过大，对硬件要求较高，快速傅里叶变换虽然降低了离散傅里叶变换的计算量，但是由 于存在栅栏效应无法从混叠的频谱中提取间谐波;基于小波变换的分析方法，由于不同小 波滤波器的分频功能以及能量泄露差别等问题都尚不清楚，因此有必要提出一种方案来实 现在系统参数随机变化下，精确测量出电力系统的交流负载功率，为电力系统稳定运行提 供有效的控制策略。
【发明内容】
[0005] 本发明为了克服上述现有技术的不足，提供了一种无功功率自适应计算方法，实 现在系统参数随机变化下，精确测量出电力系统的交流负载功率，为电力系统稳定运行提 供有效的控制策略。
[0006] 为实现上述目的，本发明采用了以下技术措施：
[0007] -种无功功率自适应计算方法，包括以下步骤：
[0008] S1、设定采样周期为T，定时对电网负载端的电压信号和电流信号进行同步采样， 得到同一时刻的电压采样数据u(t)和电流采样数据i(t);
[0009] S2、对电压采样数据u(t)或电流采样数据i(t)进行工况匹配分析，获得交流分量 个数N;
[0010] S3、对所述电压采样数据u(t)进行希尔伯特变换；
[0011] S4、采用改进扩展Prony方法拟合出N组交流参数分别为交流电压ui、U2、…、un、希 尔伯特变换后的电压UW2、…、U、、交流电流h、i2、…、iN、以及频率fl、f2、-_、f N;
[0012] S5、计算出每一个交流分量的无功功率IQ2、…、Qn。
[0013] 所述步骤S2具体包括：
[0014] a2:采用频谱分析法处理所述电压采样数据u(t)或电流采样数据i(t)，获得电压 信号或电流信号的频谱信息，从而得出所述电压信号或电流信号的交流分量个数N;
[0015] b2:依据所述交流分量个数N，设定改进扩展Prony检测模型的阶数P，P = N。
[0016] 所述步骤S3具体包括：
[0017] a3:对于电压采样数据u(t)的希尔伯特变换形式为：
[0018]
，其中，C(t)=u(t)，H[C(t)]为对函数 C(t) 进行希尔伯特变换，H(t)为希尔伯特变换结果，λ为积分变量，π为圆周率；
[0019] b3:构造一个解析信号Z(t)，令Z(t) = C(t) + iH(t) = A(t)e10(t)，其中，i为单位虚 数，i2 = _l，并且，
[0020] 幅值函$ .......单位:V
[0021 ] 幅角函数：
单位:rad
[0022] 瞬时频率：
单位:Hz。
[0023]所述步骤S4具体包括：
[0024] a4:构造扩展阶矩阵
[0025] 将所述电压采样数据u(t)或电流采样数据i(t)作为观测数据序列x(n)和样本函 数公式r (i，j)，构造出扩展阶矩阵歹
[0026]
[0027] 其中，
>p，pe和P分别为所述扩展矩阵瓦的阶 数、以及有效秩，xT(n)表示求一维向量x(n)的转置，N为交流分量的个数，并且0<i，j<N-1;
[0028] b4:对扩展阶矩阵互求解特征多项式系数
[0029] 采用总体最小二乘估计方法求解法方程：
[0030]
[0031] 定义最小误差能量εΡ为：
[0032]其中，aj为特征多项式系数，0<j<p;
[0033] 〇4:计算21，22，."，乙[
[0034] 首先，利用b4中得到的系数ai，a2，…，aP，求解特征多项式
的根，即可得到一组复指数Z1，Z2,…，Ζρ，所述Z1，Z2，···，Ζρ为特征多项式的Φ的根，令表 示X(n)的估计值，再利用估计值i(n)的递推差分方程式
塞个计 算出i.(ra)，其中，n = 0,…，N-1,i(0) = -t(o);
[0035] d4:计算bib,…，bP
[0036] 构造方程组并代换成矩阵形式，可得：?(") (1)
[0037] 其中
C2)
[0038] b=[bi,b2,---,bP]T
[0039]
[0040] 公式(2)所示矩阵Z是一个N维的范德蒙矩阵，由于矩阵Z中的各个列均不相同，根 据范德蒙矩阵的性质可知矩阵Z的各列线性独立，即矩阵Z是一个列满秩的矩阵，因此公式 (1)采用最小二乘法求解的结果为? = 〃），其中，ZH为矩阵z的共辄转置矩阵；
[0041] e4:计算电压采样数据u(t)或电流采样数据i(t)的频域信息
[0042] 在ζι，Ζ2,…，zp和bi，b2,…，bp的计算结果基础上，利用公式⑶可以依次计算出所 述观测数据序列X (η)的各特征量包括幅值Am、相角3m、角频率fm，
[0043]
[0044] 其中：T表示采样时间间隔，π为圆周率，Re表示取复数的实部，Im表示取复数的虚 部；
[0045] 所述电流采样数据i(t)通过公式(3)处理后，输出N组电流信号谐波分量的电流幅 值Ik，单位为A、电流相位誕，单位为rad、以及电流频率fk，单位为Hz，所述电流幅值Ik、电流 相位也、以及电流频率fk满足公式(4)
[0046]
(4)
[0047] 所述电压采样数据u(t)通过公式(3)处理后，输出N组电压信号谐波分量的电压幅 值Uk，单位为V、电压相位約:，单位为rad、以及电压频率f k，单位为Hz，所述电压幅值Uk、电压 相位涔、以及电压频率fk满足公式(5)
[0048]
(5)
[0049]所述电压采样数据u(t)希尔伯特变换得到Y (thY (t)通过公式(3)处理后，输出 N组电压信号谐波分量的电压幅值lTk，单位为V、电压相位单位为rad、以及电压频率fk， 单位为Hz，所述电压幅值l/k、电压相位<、以及电压频率f k满足公式(6)
[0050]
(6；)
[0051] 其中Ui^U'k，式=% +寻 6 2
[0052]所述步骤S5具体包括：
[0053] a5:依据IEEE标准无功功率Q的功率计算公式为
计算出各个交流分量的无功功率Qi、Q2、…、Qn，其中，9k为电压相位和电流相位差，且 单位为= 单位为 var，Qi 为基波无功功率；
,单位为var，Qh为所有谐波无功功率总和。
[0054] 本发明的有益效果在于：
[0055] 1)、本发明公开了一种无功功率自适应计算方法，该方法包括首先设定采样周期 为T，定时对电网负载端的电压信号和电流信号进行同步采样，得到同一时刻负载的电压采 样数据u(t)和电流采样数据i(t);然后对所述电压采样数据u(t)或电流采样数据i(t)进行 工况匹配分析，获得交流分量个数N;再对所述电压采样数据u(t)进行希尔伯特变换;最后 采用改进扩展Prony方法拟合出交流电压ui、U2、…、un和希尔伯特变换后的电压ι/ 1、ι/ 2、…、 u%，交流电流ii、i2、…、iN，频率fi、f2、…、fN;最后计算出每一个交流分量的无功功率&、 Q2、…、Qn ;本发明实现在系统参数随机变化下，精确测量出电力系统的交流负载功率，为电 力系统稳定运行提供有效的控制策略，而且本方法执行简单、精确度高、易于实现，可以广 泛应用于电网电能功率计量中。
【附图说明】
[0056] 图1为本发明的流程图；
[0057] 图2为本发明的负载二阶电路图。
【具体实施方式】
[0058]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0059] 设定一个随机可变的电流信号：
[0060]
[0061] 如图1所示，一种负载随机可调下的无功功率自适应计算方法，包括以下步骤：
[0062] S1、设定采样周期T = 0.001s，定时对图2中输入端的电压信号和电流信号进行同 步采样，得到同一时刻负载的电压采样数据u(t)，单位为V，电流采样数据i(t)，单位为A;
[0063] S2、采用频谱分析方法处理电压采样数据u(t)或电流采样数据i(t)，获得电压信 号或电流信号的频谱信息，从而得到电压信号或电流信号的交流分量个数为2,依据交流分 量个数设定扩展Prony检测模型的阶数P = 2。
[0064] S3、对所述电压采样数据u(t)进行希尔伯特变换的具体步骤包括：
[0065] a3:对于电压采样数据u(t)的希尔伯特变换形式为：
[0066]
其中，C(t) = u(t)，H[C(t)]为对函数C(t) 进行希尔伯特变换，H(t)为希尔伯特变换结果，λ为积分变量，π为圆周率；
[0067] b3:构造一个解析信号Z(t)，令Z(t) = C(t) + iH(t) = A(t)e10(t)，其中，i为单位虚 数，i2 = _l，并且，
[0068] 幅值函数： _
丨，单位:V
[0069] 幅角函数：
Μ立:rad
[0070] 瞬时频率：
，单位:Hz。
[0071] S4、采用改进扩展Prony方法拟合出N组交流参数分别为交流电压m、u2,单位为V， 希尔伯特变换后的电压11 /1、11/2，单位为￥，交流电流11、12，单位为4，频率心、5，单位为他的 具体操作步骤为：
[0072] a4:构造扩展阶矩阵
[0073] 将所述电压采样数据u(t)或电流采样数据i(t)作为观测数据序列x(n)和样本函 数公式r (i，j)，构造出扩展阶矩阵瓦:
[0074]
[0075]其中
hPe和P分别为所述扩展矩阵1的阶 数、以及有效秩，xT(n)表示求一维向量x(n)的转置，N为交流分量的个数，并且0<i，j<N-1;
[0076 ] b4:对扩展阶矩阵瓦求解特征多项式系数
[0077] 采用总体最小二乘估计方法求解法方程：
[0078]
[0079]定义最小误差能量ερ为：
[0080] 其中，aj为特征多项式系数，0<j<p;
[0081] 〇4:计算21，22，."，乙[)
[0082] 首先，利用b4中得到的系数&1，&2，"_，如求解特征多项岁
" 的根，即可得到一组复指数Ζ1，Ζ2,…，Ζρ，所述Ζ1，Ζ2，···，Ζρ为特征多项式的Φ的根，令对〃)表 示X(η)的估计值，再利用估计值i(?)的递推差分方程式
?个计 算出场?