CN108802491B - 基于谐波源独立性的谐波发射水平估计求解新方法 - Google Patents

Abstract

本发明公开了基于谐波源独立性的谐波发射水平估计求解新方法，所述方法包括以下步骤：步骤(A)：对一定时间段内的公共耦合点电压、电流值进行采样，使用采用的量测数据建立量测数据矩阵；步骤(B)：对步骤(A)中的量测数据矩阵，计算其协方差矩阵；同时根据已知的电路方程，建立其协方差矩阵的另一个表达式，根据两个协方差矩阵表达式的对应关系，建立阻抗计算方程；步骤(C)：基于PCC点两侧谐波源独立的特性，对步骤(B)中的协方差矩阵的另一个表达式进行简化，由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗。

Description

基于谐波源独立性的谐波发射水平估计求解新方法

技术领域

本发明涉及谐波源定位中的谐波发射水平研究领域，具体地，涉及一种基于谐波源独立性的谐波发射水平估计求解新方法。

背景技术

随着电网的不断发展，越来越多的电力电子元件接入了电网中，从而引发的谐波问题越来越引起重视，如何定量划分接入***中各谐波源的谐波责任是实现电能质量智能管理、控制的前提。

通常采用等效电路的方法。即选定一个公共连接点(PCC),将PCC两端的***侧电路和用户侧电路等效成谐波电压源和谐波电流源，并定义***侧和用户侧谐波参考阻抗，将谐波责任表示成PCC处量测电压、电流和参考阻抗，并划分两侧的谐波责任。

目前进行谐波责任划分的方法主要有波动量法、线性回归方法、独立分量法等。

综上所述，本申请发明人在实现本申请发明技术方案的过程中，发现上述技术至少存在如下技术问题：

在现有技术中，现有的波动量方法、线性回归方法存在估计误差随着两侧谐波源发射水平逐渐接近，误差增长速度较快的技术问题。波动量法根据波动量比值特征符号进行估计，在某些场景下估计精度不能保证；线性回归方法适用于***侧和用户侧平稳的条件，在波动量较大时拟合效果不好。

发明内容

本发明的目的在于提供一种根据谐波源独立性计算谐波阻抗的计算方法，依据谐波源之间的独立性，通过PCC处测量得到的电压、电流波动量，建立量测数据矩阵并计算其协方差矩阵；同时通过等效电路建立PCC处测量电压、电流与谐波源电流波动量、参考阻抗的方程；由***侧谐波源电流波动量远小于用户侧谐波源电流波动量的特性，可对方程进行简化，并求解谐波阻抗和谐波发射水平。

为实现上述发明目的，本申请提供了基于谐波源独立性的谐波发射水平估计求解新方法，包括以下步骤：

步骤(A)：对一定时间段内的公共耦合点电压、电流值进行采样，使用采用的量测数据建立量测数据矩阵；

步骤(B)：对步骤(A)中的量测数据矩阵，计算其协方差矩阵；同时根据已知的电路方程，建立其协方差矩阵的另一个表达式，根据两个协方差矩阵表达式的对应关系，建立阻抗计算方程；

步骤(C)：基于PCC点两侧谐波源独立的特性，对步骤(B)中的协方差矩阵的另一个表达式进行简化，由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗。

进一步的，所述步骤(A)中的量测数据矩阵建立按照以下步骤进行：

步骤(A1)：假设采样电压、电流为V_pcc、I_pcc，采样点数为N；对采样值去均值后，建立量测数据矩阵A；

上式中ΔV_pcc为PCC点电压波动量，ΔI_pcc为PCC点电流波动量；

步骤(A2)：上式中的ΔV_pcc、ΔI_pcc可由下式计算得到：

ΔV_pcc＝V_pcc(n+1)-V_pcc(n)n＝1，2，...，N-1；

ΔI_pcc＝I_pcc(n+1)-I_pcc(n)n＝1，2，...，N-1。

进一步的，所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行：

步骤(B1)：量测数据矩阵A的协方差矩阵(1)建立如下：

上式中上注符^H表示共轭转置；上注符^*表示共轭；

步骤(B2)：根据PCC诺顿等效电路可建立方程：

A＝ZI

上式中，Z_u、Z_c分别为PCC诺顿等效电路***侧和用户侧的等效阻抗；ΔI_u和ΔI_c分别为***侧和用户侧电流波动量，Z为等效阻抗矩阵，I为PCC两侧电流波动量矩阵；

步骤(B3)：由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下：

AA^H＝(ZI)(ZI)^H＝Z(II^H)Z^H (2)。

进一步的，所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行：

步骤(C1)：II^H可简化如下：

其中，λ₁、λ₂为PCC***侧和用户侧的电流波动量幅值平方和，可由下式表示：

步骤(C2)：由(B2)计算协方差矩阵AA^H如下：

步骤(C3)：联立(B1)和(C2)中协方差矩阵AA^H计算方程可建立方程组：

步骤(C4)：(C3)中方程可简化如下，求解简化后方程即可求解谐波阻抗并计算谐波发射水平：

本发明与现有技术相比，具有如下的优点和有益效果：

本发明基于谐波源独立性的谐波发射水平估计求解新方法，在建立量测数据协方差矩阵时可忽略谐波源电流波动量之间交叉相乘项，减少了PCC两侧谐波源的相互影响对计算结果的影响；2本发明基于谐波源独立性的谐波发射水平估计求解新方法，可减少谐波责任计算的计算量，加快计算速度；3本发明基于谐波源独立性的谐波发射水平估计求解新方法，算法计算结果较为准确快速，误差较小，具有良好的实用性。

附图说明

此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定；

图1为本发明流程示意图；

图2为实施例一的PCC点等效电路模型图。

具体实施方式

本发明提供了一种基于双目立体视觉的接触网几何参数检测测量***，解决了现有的基于单目视觉的检测方法存在先天的误差的技术问题，实现了测量***设计合理，几何参数检测测量结果准确的技术效果。

为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在相互不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。

在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述范围内的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。

请参考图1，本申请提供了基于谐波源独立性的谐波发射水平估计求解新方法，所述方法包括以下步骤：

步骤(A)：对一定时间段内的公共耦合点电压、电流值进行采样，使用采用的量测数据建立量测数据矩阵；

步骤(B)：对步骤(A)中的量测数据矩阵，计算其协方差矩阵；同时根据已知的电路方程，建立其协方差矩阵的另一个表达式，根据两个协方差矩阵表达式的对应关系，建立阻抗计算方程；

步骤(C)：基于PCC点两侧谐波源独立的特性，对步骤(B)中的协方差矩阵的另一个表达式进行简化，由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗；

实施例一：在Matlab中搭建如图2所示的PCC点等效电路，各元件参数表1所示，使用本发明对谐波阻抗进行求解。并将计算结果、计算误差与其他方法进行比较如表2-表3所示。

表1

表2

表3

其中，在本申请实施例中，所述步骤(A)中的量测数据矩阵建立按照以下步骤进行：

步骤(A1)：假设采样电压、电流为V_pcc、I_pcc，采样点数为N；对采样值去均值后，建立量测数据矩阵A；

式中ΔV_pcc为PCC点电压波动量，ΔI_pcc为PCC点电流波动量。

步骤(A2)：式(1)中的ΔV_pcc、ΔI_pcc可由下式计算得到：

ΔV_pcc＝V_pcc(n+1)-V_pcc(n)n＝1，2，...，N-1；

ΔI_pcc＝I_pcc(n+1)-I_pcc(n)n＝1，2，...，N-1；

其中，在本申请实施例中，所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行：

步骤(B1)：量测数据矩阵A的协方差矩阵(1)建立如下：

式中上注符^H表示共轭转置；上注符^*表示共轭。

步骤(B2)：根据PCC诺顿等效电路可建立方程：

A＝ZI

式中，Z_u、Z_c分别为PCC诺顿等效电路***侧和用户侧的等效阻抗；ΔI_u和ΔI_c分别为***侧和用户侧电流波动量，Z为等效阻抗矩阵，I为PCC两侧电流波动量矩阵。

步骤(B3)：由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下：

AA^H＝(ZI)(ZI)^H＝Z(II^H)Z^H (2)

其中，在本申请实施例中，所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行：

步骤(C1)：II^H可简化如下：

其中，λ₁、λ₂为PCC***侧和用户侧的电流波动量幅值平方和，可由下式表示：

步骤(C2)：由(B2)计算协方差矩阵AA^H如下：

步骤(C3)：联立(B1)和(C2)中协方差矩阵AA^H计算方程可建立方程组：

步骤(C4)：(C3)中方程可简化如下，求解简化后方程即可求解谐波阻抗并计算谐波发射水平：

尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (1)

1.基于谐波源独立性的谐波发射水平估计求解新方法，其特征在于，所述方法包括以下步骤：

步骤(A)：对一定时间段内的公共耦合点电压、电流值进行采样，使用采用的量测数据建立量测数据矩阵；

步骤(B)：对步骤(A)中的量测数据矩阵，计算其协方差矩阵；同时根据已知的电路方程，建立其协方差矩阵的另一个表达式，根据两个协方差矩阵表达式的对应关系，建立阻抗计算方程；

步骤(C)：基于PCC点两侧谐波源独立的特性，对步骤(B)中的协方差矩阵的另一个表达式进行简化，由两协方差矩阵表达式对应相等建立方程组求解谐波阻抗；

所述步骤(A)中的量测数据矩阵建立按照以下步骤进行：

步骤(A1)：假设采样电压、电流为V_pcc、I_pcc，采样点数为N；对采样值去均值后，建立量测数据矩阵A；

上式中ΔV_pcc为PCC点电压波动量，ΔI_pcc为PCC点电流波动量；

步骤(A2)：上式中的ΔV_pcc、ΔI_pcc可由下式计算得到：

ΔV_pcc＝V_pcc(n+1)-V_pcc(n)n＝1,2,…,N-1；

ΔI_pcc＝I_pcc(n+1)-I_pcc(n)n＝1,2,…,N-1；

所述步骤(B)中的协方差矩阵的两个表达式的建立按照以下步骤进行：

步骤(B1)：量测数据矩阵A的协方差矩阵(1)建立如下：

上式中上注符^H表示共轭转置；上注符^*表示共轭；

步骤(B2)：根据PCC诺顿等效电路可建立方程：

A＝ZI

上式中，Z_u、Z_c分别为PCC诺顿等效电路***侧和用户侧的等效阻抗；ΔI_u和ΔI_c分别为***侧和用户侧电流波动量，Z为等效阻抗矩阵，I为PCC两侧电流波动量矩阵；

步骤(B3)：由电路方程改写量测数据矩阵A的协方差矩阵的另一个表达式(2)如下：

AA^H＝(ZI)(ZI)^H＝Z(II^H)Z^H (2)；

所述步骤(C)中的协方差矩阵简化和谐波阻抗计算按照以下步骤进行：

步骤(C1)：II^H可简化如下：

其中，λ₁、λ₂为PCC***侧和用户侧的电流波动量幅值平方和，可由下式表示：

步骤(C2)：由(B2)计算协方差矩阵AA^H如下：

步骤(C3)：联立(B1)和(C2)中协方差矩阵AA^H计算方程可建立方程组：

步骤(C4)：(C3)中方程可简化如下，求解简化后方程即可求解谐波阻抗并计算谐波发射水平：