CN105588982A - 一种电力***谐波阻抗的测量计算方法 - Google Patents

Abstract

本发明公开了一种电力***谐波阻抗的测量计算方法，将复最小二乘法引入电力***谐波阻抗计算，将测量电压、测量电流等相量实部和虚部作为一个有机的整体进行计算，从而得到最优最小二乘解。本发明引入数据处理过程，采用基于相对复残差排序的方法对原始数据进行筛选，有效排除异常点对回归系数的影响。在此原始数据筛选的基础上，采用复数域反复加权的方法进行循环回归计算，有效减小强影响点对回归系数的影响，增加回归的稳健性。

Description

一种电力***谐波阻抗的测量计算方法

技术领域

本发明涉及谐波阻抗计算方法设计领域，特别是一种电力***谐波阻抗的测量计算方法。

背景技术

电力***谐波阻抗对电力***状态评估及接入新电力设备的设计十分重要，在电力***建模、仿真以及线路保护等电力装置的配置中起着不可或缺的作用。同时，随着国内外电力电子技术的发展，大量由电力电子开关构成的，具有非线性特性的用电设备广泛应用于冶金、钢铁、交通、化工等需求领域，如电解装置、电气机车、轧制机械、高频炉等，导致电力***中谐波问题日益严重，对谐波相关问题，如谐波潮流计算、谐波责任划分、谐波发射水平、谐波源识别等的研究一直是电力行业关注的问题。而***谐波阻抗的准确测量是研究谐波相关问题基础。

目前，谐波阻抗测量估计方法主要分为“干预式”和“非干预式”两大类。“干预式”方法主要是通过人为方式产生扰动，如向***注入谐波电流或是间谐波电流，或是通过开断***某一支路来进行***侧谐波阻抗的测量。这种扰动可能会对电力***正常运行产生不利的影响，因此该类方法不能得到广泛应用。“非干预式”方法则是利用***或谐波源负荷本身的扰动，通过可测量参数来进行相关计算，不会影响到***的安全、稳定运行，因而成为当前相关课题研究的主流方法。典型的“非干预式”方法有波动量法和线性回归法。线性回归法在谐波发射水平评估方面占有重要地位，是“非干预式”方法的核心。传统的线性回归方法包括二元回归法、稳健回归法和偏最小二乘法，这些方法均没有数据的预处理过程，其本质都是将相量实部、虚部分割开来，在实数域内利用最小二乘法分别进行回归系数的求解。这种处理方法忽视了相量是作为一个有特定物理意义的整体而存在的，其求解过程并不是实部、虚部结果的简单相加和组合。

发明内容

本发明的目的是要提供一种电力***谐波阻抗的测量计算方法。

为达到上述目的，本发明是按照以下技术方案实施的：

一种电力***谐波阻抗的测量计算方法，包括下列步骤:

步骤1：采集公共连接点的母线电压瞬时值和用户注入***的电流瞬时值，通过傅里叶变换得到谐波电压和谐波电流相量数据序列，并建立谐波电压和谐波电流相量的复回归方程；

步骤2：对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，排除异常点对回归系数的影响；

步骤3：在初步筛选基础上，求解***谐波阻抗回归系数；求出最优的赋权回归系数，并得到谐波阻抗的值。

进一步的，对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，所述初步筛选采用基于相对复残差排序的方法，定义复残差列向量，e=[e₁,e₂,e₃…e_n]＇，其中，指没有经过数据预处理过程回归参数的估计值，Y为谐波电压相量数据序列组成的列向量，X为谐波电流相量数据序列组成的增广矩阵，为X的共轭转置矩阵；定义相对复残差，（i=1,2，…，n），e_j为第i个复残差值，Y_j为第i个谐波电压相量，n为列向量的个数；将e_ri由小到大排序，得到相对复残差序列{e_ri}，删除处于尾部5%相对复参差e_ri所对应的谐波电压和谐波电流数据，并将余下谐波电压和谐波电流相量记为Y_95%、X_95%。

具体的，所述求解***谐波阻抗回归系数，采用Ramsay影响函数赋权的反复加权复最小二乘的稳健回归方法；具体步骤为：

a）选择剩余的谐波电压X_95%和谐波电流Y_95%，根据式，计算回归参数；为X_95%的共轭转置矩阵；

b）由式，计算复残差f=[f₁,f₂,f₃,…,f_m]＇,即f={f_j}，（i=1,2,3…，m，m=0.95n）；

c）由式计算Ramsay影响函数自变量，其中，=1.5med（）；med（）是的中位数；由式，计算基于Ramsay影响函数的赋权函数，其中a为调整Ramsay影响函数形状的调整因子；

d）取赋权函数的值=（i=1,2,3…，m），记权重对角阵为W=diag（），（i=1,2,3…，m，m=0.95）；

e）由加权复最小二乘计算公式：；计算出赋权的回归参数，即；

f）计算，即计算相邻两次回归参数的差值；其中j代表回归计算循环次数，其初始值为0，当j=0时，即是步骤a求出的回归参数，而是步骤e求出的回归参数；和都是复数；

g）计算和模值和，并判断max{，}<是否成立，若成立，则执行步骤h；否则，令=，返回步骤b；

h）由即可求出谐波阻抗。

与现有技术相比，本发明的有益效果：

将复最小二乘法引入电力***谐波阻抗计算课题，将测量电压、测量电流等相量实部和虚部作为一个有机的整体进行计算，从而得到最优最小二乘解。同时，本发明引入数据处理过程，采用基于相对复残差排序的方法对原始谐波电压电流数据进行筛选，有效排除异常点对回归系数的影响。在此原始谐波电压和谐波电流数据筛选的基础上，采用复数域反复加权的方法进行循环回归计算，有效减小强影响点对回归系数的影响，增加回归的稳健性。

附图说明

图1为本发明提供的一种电力***谐波阻抗的测量计算方法的***和用户等值电路示意图；

图2为本发明提供的一种电力***谐波阻抗的测量计算方法的电力***谐波阻抗计算流程图；

图3为本发明提供的一种电力***谐波阻抗的测量计算方法的IEEE14节点标准测试***；

图4为本发明提供的一种电力***谐波阻抗的测量计算方法的谐波源HL1注入的5次波电流；其中，a是5次谐波电流实部；b是5次电流虚部；

图5为本发明提供的一种电力***谐波阻抗的测量计算方法的母线11处5次谐波电压有效值曲线；

图6为本发明提供的一种电力***谐波阻抗的测量计算方法的谐波源HL1注入母线11的5次谐波电流有效值曲线。

具体实施方式

下面结合具体实施例对本发明作进一步描述，在此发明的示意性实施例以及说明用来解释本发明，但并不作为对本发明的限定。

如图1、图2所示的一种电力***谐波阻抗的测量计算方法，包括下列步骤:

步骤1：采集公共连接点的母线电压瞬时值和用户注入***的电流瞬时值，通过傅里叶变换得到谐波电压和谐波电流相量数据序列，并建立谐波电压和谐波电流相量的复回归方程；

步骤2：对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，排除异常点对回归系数的影响；

步骤3：在初步筛选基础上，求解***谐波阻抗回归系数；求出最优的赋权回归系数，并得到谐波阻抗的值。

对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，所述初步筛选采用基于相对复残差排序的方法，定义复残差列向量，e=[e₁,e₂,e₃…e_n]＇，其中，指没有经过数据预处理过程回归参数的估计值，Y为谐波电压相量数据序列组成的列向量，X为谐波电流相量数据序列组成的增广矩阵，为X的共轭转置矩阵；定义相对复残差，（i=1,2，…，n），e_j为第i个复残差值，Y_j为第i个谐波电压相量，n为列向量的个数；将e_ri由小到大排序，得到相对复残差序列{e_ri}，删除处于尾部5%相对复参差e_ri所对应的谐波电压和谐波电流数据，并将余下谐波电压和谐波电流相量记为Y_95%、X_95%。

求解***谐波阻抗回归系数，采用Ramsay影响函数赋权的反复加权复最小二乘的稳健回归方法；具体步骤为：

a）选择剩余的谐波电压X_95%和谐波电流Y_95%，根据式，计算回归参数；为X_95%的共轭转置矩阵；

b）由式，计算复残差f=[f₁,f₂,f₃,…,f_m]＇,即f={f_j}，（i=1,2,3…，m，m=0.95n）；

c）由式计算Ramsay影响函数自变量，其中，=1.5med（）；med（）是的中位数；由式，计算基于Ramsay影响函数的赋权函数，其中a为调整Ramsay影响函数形状的调整因子；

d）取赋权函数的值=（i=1,2,3…，m），记权重对角阵为W=diag（），（i=1,2,3…，m，m=0.95）；

e）由加权复最小二乘计算公式：；计算出赋权的回归参数，即；

f）计算，即计算相邻两次回归参数的差值；其中j代表回归计算循环次数，其初始值为0，当j=0时，即是步骤a求出的回归参数，而是步骤e求出的回归参数；和都是复数；

g）计算和模值和，并判断max{，}<是否成立，若成立，则执行步骤h；否则，令=，返回步骤b；

h）由即可求出谐波阻抗。

图3是本发明提供的一种电力***谐波阻抗的测量计算方法的IEEE14节点标准测试***。图3中，G是发电机，C是调相机，是双绕组变压器，利用IEEE14节点标准测试***，对本发明所提出的算法进行仿真验证，该测试***由2台发电机、3台同步调相机、14条母线、15条输电线路和3台变压器组成。选择母线11为关注母线，HL1及L2为母线11所接负荷，其中HL1为非线性负荷，代表谐波源用户，L2为线性负荷。以母线11为公共连接处母线，将负荷HL1单独视为用户侧，将除去负荷HL1的网络其余部分视为该仿真算例的***侧，以***侧正常工作时产生的谐波作为分析问题时的***谐波，这与图1所示等值电路吻合。本发明以5次谐波为例进行仿真验证。

谐波源HL1的5次谐波电流曲线中一分钟1个样本点，一天对应1440个点，如图4所示。测量母线11处谐波电压有效值曲线和谐波源HL1注入母线11谐波电流有效值曲线分别如5和图6所示。

采用本发明提出的基于稳健复最小二乘法计算***等值谐波阻抗，计算结果与真实值对比如下表所示，符号“-”表示此值不需计算。

表1***等值谐波阻抗计算结果对比

从表1可以看出，稳健复最小二乘法可以精确计算***等值谐波阻抗。

上述实施的仿真实验表明，基于相对复残差排序的原始谐波电压和谐波电流数据序列筛选和复数域反复加权的稳健复最小二乘法，可以有效消除或减小异常点及强影响点对回归系数的影响，提高***谐波阻抗求解精度。

本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案做出的技术变形，均落入本发明的保护范围之内。

Claims (3)

1.一种电力***谐波阻抗的测量计算方法，其特征在于，包括下列步骤:

步骤1：采集公共连接点的母线电压瞬时值和用户注入***的电流瞬时值，通过傅里叶变换得到谐波电压和谐波电流相量数据序列，并建立谐波电压和谐波电流相量的复回归方程；

步骤2：对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，排除异常点对回归系数的影响；

步骤3：在初步筛选基础上，求解***谐波阻抗回归系数；求出最优的赋权回归系数，并得到谐波阻抗的值。

2.根据权利要求1所述的电力***谐波阻抗计算方法，其特征在于，对步骤1得到的谐波电压和谐波电流相量数据序列进行初步筛选，所述初步筛选采用基于相对复残差排序的方法，定义复残差列向量，e=[e₁,e₂,e₃…e_n]＇，其中，指没有经过数据预处理过程回归参数的估计值，Y为谐波电压相量数据序列组成的列向量，X为谐波电流相量数据序列组成的增广矩阵，为X的共轭转置矩阵；定义相对复残差，（i=1,2，…，n），e_j为第i个复残差值，Y_j为第i个谐波电压相量，n为列向量的个数；将e_ri由小到大排序，得到相对复残差序列{e_ri}，删除处于尾部5%相对复参差e_ri所对应的谐波电压和谐波电流数据，并将余下谐波电压和谐波电流相量记为Y_95%、X_95%。

3.根据权利要求1所述的电力***谐波阻抗的测量计算方法，其特征在于，所述求解***谐波阻抗回归系数，采用Ramsay影响函数赋权的反复加权复最小二乘的稳健回归方法；具体步骤为：

a）选择剩余的谐波电压X_95%和谐波电流Y_95%，根据式，计算回归参数；为X_95%的共轭转置矩阵；

b）由式，计算复残差f=[f₁,f₂,f₃,…,f_m]＇,即f={f_j}，（i=1,2,3…，m，m=0.95n）；

c）由式计算Ramsay影响函数自变量，其中，=1.5med（）；med（）是的中位数；由式，计算基于Ramsay影响函数的赋权函数，其中a为调整Ramsay影响函数形状的调整因子；

d）取赋权函数的值=（i=1,2,3…，m），记权重对角阵为W=diag（），（i=1,2,3…，m，m=0.95）；

e）由加权复最小二乘计算公式：；计算出赋权的回归参数，即；

f）计算，即计算相邻两次回归参数的差值；其中j代表回归计算循环次数，其初始值为0，当j=0时，即是步骤a求出的回归参数，而是步骤e求出的回归参数；和都是复数；

g）计算和模值和，并判断max{，}<是否成立，若成立，则执行步骤h；否则，令=，返回步骤b；

h）由即可求出谐波阻抗。