抗差输电线路正序参数辨识方法与流程

本发明涉及电力系统运行控制技术领域，尤其涉及一种抗差输电线路正序参数辨识方法。

背景技术：

随着电力系统的不断发展，输电线路也越来越多，系统也就变的越来越复杂，因此我们对输电线路的参数准确度的要求也就比较高。准确的线路参数是正确进行潮流计算，故障分析，继电保护整定，网损计算等各种电力系统计算的基础。线路参数不准确时，对短路电流的计算，稳定计算，潮流计算的正确性影响很大。所以输电线路参数的准确性也就影响着电力系统的安全，稳定和可靠运行。

输电线路可以分为架空输电线路和电缆线路，是电力输送的载体，其电压等级远高于配电网，且主网的大部分参数是已知的，与配电网有很大的不同。

线路参数测量方法有很多，有理论计算，离线停电测量，在线带电测量。传统线路理论计算以carson为基础，利用线路几何均距，材料结构等物理参数，结合气温，地理位置等根据公式计算出电抗，电阻，电纳。或者根据手册公式，直接计算线路工频参数，这种方法忽略了架空地线的影响。而且理论计算只考虑了完全对称的情况，没有考虑到实时温度，弧垂等实际问题，因此该方法的结果和实际参数存在很大的差异。离线停电测量是将被测线路停电，并脱离电网，采用外加电源，利用电压表、电流表等各种表计测量线路数据，经人工读取表值并结合相应的公式计算出各个参数。该方法存在需要将被测线路停电，多回运行的平行线路无法测量互感等问题。

技术实现要素：

本发明的目的是提供一种抗差输电线路正序参数辨识方法，计算简单，且不受外界环境，负荷大小影响，可以准确辨识线路参数，实用性好，具有很好的抗差能力。

本发明的目的是通过以下技术方案实现的：

一种抗差输电线路正序参数辨识方法，包括：

通过单回输电线路两端装设的相量量测装置测量得到线路两端的三相电流与电压相量；

根据单回输电线路三相等值模型并结合线路两端的三相电流与电压相量，建立线路阻抗参数和导纳参数的计算方程；

构建具有抗差能力的huber估计目标函数作为线路阻抗参数和导纳参数辨识的目标函数，从而求解得到线路阻抗矩阵和导纳矩阵参数；

采用对称分量法将线路阻抗矩阵和导纳矩阵解耦，得到正序电阻、电抗和对地导纳参数。

所述根据单回输电线路三相等值模型并结合线路两端的三相电流与电压相量，建立线路阻抗参数和导纳参数的计算方程包括：

单回输电线路两端的三相电流与电压满足如下方程：

其中，分别为m端与n端的a、b、c三相电流相量；分别为m端与n端的a、b、c三相电压相量；z、yc分别为待辨识的线路阻抗参数、导纳参数；

将上述方程转换为下述形式：

上式中，

则当前时刻，待辨识的参数x为：

x＝[x11x12...x19x21x22...x29]^t；

满足如下线性方程：

ax＝b+v；

其中，a为当前时刻三相电压与电流相量构成的系数矩阵，b为当前时刻两端电压降相量组成的常数项，v为方程残差相量；

辨识的参数x中包含18个未知量，则方程的个数要大于等于18个，取多组不同时刻的数据求解，其中，a的维度是6n×18，n≥3，b的维度是6n×1，具体目标方程如下：

所述构建具有抗差能力的huber估计目标函数作为线路阻抗参数和导纳参数辨识的目标函数，从而求解得到线路阻抗矩阵和导纳矩阵包括：

huber抗差估计理论是假设实际观测数据服从huber分布，huber分布是污染分布的一种，其主体是正态分布，干扰部分服从laplace分布；

huber分布的概率密度为：

其中，ε是污染率,它表示污染部分的数据在整个数据中所占的比例；λ是观测值；

为标准正态分布密度，在区间-c≤λ≤c内，观测值服从正态分布；在λ＞c时，观测值服从laplace分布；c的取值在1.0～2.0之间；

huber分布的极大似然估计为huber估计，其目标函数为：

其中，k为抗差阈值调制系数，与污染率ε有关；vi为量测残差，vi＝(yci-ymi)/s，yci和ymi分别表示输出矢量的计算值和测量值；s在|vi|≤k区间，取量测误差的标准差σi，s在|vi|＞k区间，取kmad，kmad＝med|yci-ymi|；根据目标函数，当|vi|＞k时，降低大偏差量测影响；

要建立的huber估计目标函数是方程残差的函数，可定义为：

pi为第i个量测量的权重，权值为1/σi²，ρ(vi)为huber法对应的目标函数，vi为量测残差；

考虑pmu量测误差一致，建立huber估计目标函数为：

具体的，根据线性方程ax＝b+v，此估计目标函数可以进一步表示为：

其中，aixi-bi表示采样时刻为i时的残差，最后便是求取满足此目标函数为极小值时的相量x；n为采样总数；该huber估计目标函数在观测误差小于阈值时即为最小二乘方法的目标函数，在观测误差大于阈值时，能够削弱其对目标函数的不利影响，自动剔除主网相量量测装置数据中的坏数据；

再利用优化方法求解可得线路阻抗矩阵和导纳矩阵。

采用对称分量法将线路阻抗矩阵和导纳矩阵解耦，得到正序电阻、电抗和对地导纳参数包括：

构造矩阵t：

此处构造的矩阵t是一种对称分量变换，其作用是把a,b,c三相的相分量转化为正序，负序和零序分量；其中，a为旋转因子，

利用线路相分量模型建立三相电压降与三相电流关系的矩阵方程：

上式中，表示三相电压降，表示三相电流；

将三相电压降和三相电流替换为序分量，得到：

表示电压降的各序分量，即正、负、零序电压降；表示电流的各序分量，即正、负、零序电流；；

进一步表示为：

zp即为序分量的阻抗矩阵，即：

其中的z(1)为线路正序阻抗，z(2)表示线路负序阻抗，z(0)表示线路零序阻抗；

同理，对于导纳矩阵，同样推出：

其中的yc(1)为线路正序对地导纳，yc(2)表示线路负序对地导纳，yc(0)表示线路零序对地导纳。

由上述本发明提供的技术方案可以看出，仅需采集线路两端的三相电流、电压相量信息，可操作性强，易于实施，且辨识得到的正序阻抗参数精度高，具有很好的抗差能力，辨识结果更加可信。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。

图1为本发明实施例提供的一种抗差输电线路正序参数辨识方法的流程图；

图2为本发明实施例提供的单回输电线路的三相模型示意图；

图3为本发明实施例提供的500kv仿真系统示意图。

具体实施方式

下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。

本发明实施例提供一种抗差输电线路正序参数辨识方法，该方法属于在线辨识方法。在线辨识方法，可以不用对测量线路停电，即可对两端数据，如电压，电流，相角等进行同步采集，从而实现在线线路参数的辨识，而且其精度，效率更高。如图1所示，为本发明实施例提供的一种抗差输电线路正序参数辨识方法的流程图，其主要包括如下步骤：

步骤11、通过单回输电线路两端装设的相量量测装置测量得到线路两端的三相电流与电压相量。

本发明实施例中，预先在所需要量测的单回输电线路两端装设输电线路相量量测装置，测量单回输电线路两端的三相电压、电流相量。

示例性的，单回输电线路的三相模型可以参见图2，图2中各个参数的含义将在后文逐一进行说明。

步骤12、根据单回输电线路三相等值模型并结合线路两端的三相电流与电压相量，建立线路阻抗参数和导纳参数的计算方程。

本发明实施例中，单回输电线路两端的三相电流与电压满足如下方程：

其中，分别为m端与n端的a、b、c三相电流相量；分别为m端与n端的a、b、c三相电压相量；z、yc分别为待辨识的线路阻抗参数、导纳参数；

具体来说：

式中，zkk＝rkk+jxkk，zkk为k(k＝a,b,c)相的线路自阻抗，相应的，rkk表示k相的线路串联电阻，xkk表示k相的线路串联电抗，j表示虚数单位；zmn＝rmn+jxmn，zmn(rmn，xmn)，为m(m＝a,b,c)相和n(n＝a,b,c)相之间的互阻抗(互电阻、互电抗)(m≠n)，相应的，rmn表示m相和n相之间的互电阻，xmn表示m相和n相之间的互电抗；ykk为k(k＝a,b,c)相的线路对地导纳，ymn为m(m＝a,b,c)相和n(n＝a,b,c)相之间的互导纳。zmm(rmm，xmm)，为m(m＝a,b,c)相的自阻抗(自电阻，自电抗)。

将上述方程转换为下述形式：

上式中，

则当前时刻，待辨识的参数x为：

x＝[x11x12...x19x21x22...x29]^t(7)

满足如下线性方程：

ax＝b+v(8)；

其中，a为当前时刻三相电压与电流相量构成的系数矩阵，b为当前时刻两端电压降相量组成的常数项，v为方程残差相量；

辨识的参数x中包含18个未知量，则方程的个数要大于等于18个。因此，取多组不同时刻的数据求解，其中，a的维度是6n×18，n≥3，b的维度是6n×1，具体目标方程如下：

上述矩阵方程(9)由式(4)的单回输电线路三相等值模型根据线路两侧多组量测相量获得的。由式(9)可计算获得各相电阻、电抗、电纳，即z和yc两个矩阵中的未知参数。

式(9)将两个待辨识的参数矩阵(3×3)写为了一个线性向量x(18×1)的形式，从而方便求解。

进一步，再通过对称分量法，将相参数转化为序参数，便得到了线路的正序电阻、电抗及对地电纳。后文结合huber估计的抗差方法形成目标函数进行抗差辨识的过程中需要用到上述矩阵方程，并且在实例中，运用到结合上述矩阵方程的huber估计抗差辨识算法来辨识线路的正序电阻，电抗和电纳参数。

步骤13、构建具有抗差能力的huber估计目标函数作为线路阻抗参数和导纳参数辨识的目标函数，从而求解得到线路阻抗矩阵和导纳矩阵参数。

在传统最小二乘(ls)的方法中，等权的看待每一组数据，当出现个别数据偏差较大的情形时，最终的辨识结果将大幅度偏离真实值，使得算法的抗差性能较差。

为使辨识结果具有抗差性能，可引入具有抗差能力的huber估计目标函数，huber抗差估计理论是假设实际观测数据服从huber分布，huber分布是污染分布的一种，其主体是正态分布，干扰部分服从laplace分布；

huber分布的概率密度为：

其中，ε是污染率,它表示污染部分的数据在整个数据中所占的比例；λ是观测值。

为标准正态分布密度，在区间-c≤λ≤c内，观测值服从正态分布；在λ＞c时，观测值服从laplace分布；c的取值在1.0～2.0之间，示例性的，可将c的取值设置为1.5。

huber分布的极大似然估计为huber估计，其目标函数为：

其中，k为抗差阈值调制系数，与污染率ε有关；vi为量测残差，vi＝(yci-ymi)/s，yci和ymi分别表示输出矢量的计算值和测量值；s在|vi|≤k区间，取量测误差的标准差σi，s在|vi|＞k区间，取kmad，kmad＝med|yci-ymi|；根据目标函数，当|vi|＞k时，降低大偏差量测影响。

要建立的huber估计目标函数是方程残差的函数，可定义为：

pi为第i个量测量的权重，权值为1/σi²，ρ(vi)为huber法对应的目标函数，vi为量测残差。

考虑pmu量测误差一致，建立huber估计目标函数为：

具体的，根据线性方程ax＝b+v，此估计目标函数可以进一步表示为：

aixi-bi表示采样时刻为i时的残差，最后便是求取满足此目标函数为极小值时的相量x即可。n为采样总数；该huber估计目标函数在观测误差小于阈值时即为最小二乘方法的目标函数，在观测误差大于阈值时，能够削弱其对目标函数的不利影响，自动剔除主网相量量测装置数据中的坏数据；

再利用优化方法求解可得线路阻抗矩阵和导纳矩阵，示例性的，可以通过遗传算法来计算得到线路阻抗矩阵和导纳矩阵。

步骤14、采用对称分量法将线路阻抗矩阵和导纳矩阵解耦，得到正序电阻、电抗和对地导纳参数。

通过遗传算法得到的阻抗矩阵和导纳矩阵后，可以解耦得到正序分量；具体过程如下：

构造矩阵t：

此处构造的矩阵t是一种对称分量变换，其作用是把a,b,c三相的相分量转化为正序，负序和零序分量。其中，a为旋转因子，

此式是由线路相分量模型建立的三相电压降与三相电流关系的矩阵方程。式中，表示三相电压降，表示三相电流。

将三相电压降和三相电流替换为序分量，得到：

表示电压降的各序分量，即正、负、零序电压降；表示电流的各序分量，即正、负、零序电流；阻抗矩阵z就是前文提到的相分量阻抗矩阵。

进一步表示为：

zp即为序分量的阻抗矩阵，即：

其中的z(1)为线路正序阻抗(即，正序电阻、电抗)，z(2)表示线路负序阻抗，z(0)表示线路零序阻抗。

同理，对于导纳矩阵，同样推出：

其中的yc(1)为线路正序对地导纳，yc(2)表示线路负序对地导纳，yc(0)表示线路零序对地导纳。

下面以具体的示例来对本发明提供的抗差输电线路正序参数辨识方法(以下简称，本发明方法)进行论证：

本示例利用pscad搭建500kv仿真系统，如图3所示，对单回配电线路l1的正序参数进行辨识。500kv线路l1为单回线，线路长度为200km；正序参数设计值为：电阻r1＝3.1524ω、电抗xl1＝55.2061ω、电纳b1＝7.0564×10^-4s。图3中，主要辨识的是输电线路两端(m、n)之间的正序电阻、电抗和对地电纳参数，图中p+jq表示线路末端的负荷功率，p为负荷的有功，q为负荷的无功。

假设线路l1两端装设了相量量测装置，测量到了正常运行时线路l1两端的三相不对称电压、电流相量(此处的三相电压、电流相量可以是不对称的，故该方法只需要采用线路正常运行时多个稳态时刻的两端三相不对称的量测数据即可辨识出该线路的正序参数)，采样间隔为10ms。采样时间为10s，取稳态后的900组数据并根据实施方式辨识线路l1的正序阻抗，并设置如下实验，以表明本发明所述方法的有效性。

实验一：直接利用仿真数据，不进行处理；

实验二：在理想仿真数据中叠加随机高斯噪声，模拟真实的输电线路相量量测装置数据；其中，电压电流幅值的量测误差标准差为0.1％，相角误差为0.1°；

实验三：在实验一的基础上随机将5组电流幅值量测置零，以模拟输电线路相量量测装置量测出现坏数据。

两种实验方案下，本发明方法得到的正序参数辨识结果如表1所示

表1正序参数辨识结果

表1表明，在理想仿真数据下本发明辨识结果与设计值近乎一致，表明本发明所述方法是基本可行的；实际运行中，当输电线相量量测装置获得的数据含一定量测噪声时，所述方法略优于最小二乘法。但若在输电线相量量测装置量测数据存在坏数据时，所述方法要明显优于最小二乘法，辨识结果的可信度高。

对比实验有力地表明本发明方法更适用于含量测噪声甚至坏数据的实际输电线路相量量测装置的量测数据，能有效降低坏数据对正序阻抗参数辨识的不利影响，所得正序阻抗参数值更可信。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是cd-rom，u盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。