CN114325057A - 基于逆黑盒及逆电磁对偶模型的pt一次电压重构方法 - Google Patents

Abstract

本发明涉及一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，包括通过PT采集电力系统的二次电压信号并将其分为低频电压分量和高频电压分量；对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；将一次电压高频电压分量和一次电压低频电压分量整合，得到电力系统的一次电压。该方法通过对PT采集二次电压信号划分的高频电压分量和低频电压分量分别采用逆黑盒模型和逆电磁对偶模型处理，得到一次电压高频电压分量和一次电压低频电压分量后将其相加得到一次电压，得到的一次电压不受PT采集过程中失真的影响，数据准确率高。

Description

基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法

技术领域

本发明涉及测量技术领域，尤其涉及一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法、装置及设备。

背景技术

电压测量及其在线监测是电力系统中计量、故障诊断和故障保护可靠运行的关键。实测电压波形是电力系统中包含信息最多、最具说服力的波形之一。在35kV及以下的配网中，电压系统的电压往往通过电磁式电压互感器(Potential transformer，PT)进行测量。PT是一种仪用变压器，PT的一次侧绕组与电网直接相连，PT的二次侧绕组与计量仪表相连，PT的一次绕组与PT的二次绕组之间没有直接的电路连接而是通过磁场进行耦合测量。因此，PT能够通过磁耦合实现与一次电力系统的电磁隔离，并且成本较低、测量准确、安全可靠。在电力系统中，与电压相关的故障诊断和故障保护依赖于PT二次侧输出的准确的电压信号。

PT工作在其额定频率(50/60Hz)和额定电压范围内时，可以提供准确稳定的测量结果，其电压传递特性恒定，一次电压与二次电压之间几乎没有相位差，且幅值之比为匝数比。然而，当PT一次侧(一次绕组)被高频暂态电压或低频过电压激励时，PT的二次侧信号可能会失真，与原始一次侧电压呈现显著差异，这意味着PT在这些暂态电压激励下提供的暂态电压测量结果非常不准确。失真的PT二次信号对基于电压信号的故障诊断和保护等操作造成潜在隐患。同时，失真电压信号会严重误导事故后的分析及复盘。

发明内容

本发明实施例提供了一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法、装置及设备，用于解决现有电力系统采用PT测量电压，在测量过程中PT存在失真情况，导致测量数据不准确的技术问题。

为了实现上述目的，本发明实施例提供如下技术方案：

一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，包括以下步骤：

通过PT采集电力系统的二次电压信号，并对所述二次电压信号分为低频电压分量和高频电压分量；

对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；

对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；

将所述一次电压高频电压分量和所述一次电压低频电压分量整合，得到电力系统的一次电压。

优选地，对所述二次电压信号分为低频电压分量和高频电压分量包括：

对所述二次电压信号采用傅里叶变换处理，得到二次电压信号频域形式；

对所述二次电压信号频域形式中的频率是否大于过渡频率划分为二次电压信号低频频域和二次电压信号高频频域；

分别对所述二次电压信号低频频域和所述二次电压信号高频频域采用傅里叶逆变换，得到对应的低频电压分量和高频电压分量。

优选地，对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量的步骤包括：

将所述高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

其中，所述传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数。

优选地，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换包括：

对所述传递函数进行拟合转换，得到传递函数的状态方程；

引入变量x和中心差分法对所述状态方程进行转换，得到离散电压重构函数；

采用迭代方式对所述离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，所述离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

优选地，对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量的步骤包括：

将所述低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对所述低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量；

其中，所述基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

优选地，该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法包括：将所述一次电压高频电压分量与所述一次电压低频电压分量相加整合，得到电力系统的一次电压。

本发明还提供一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置，包括：频率分量提取模块、高频反算模块、低频反算模块和整合模块；

所述频率分量提取模块，用于通过PT采集电力系统的二次电压信号，并对所述二次电压信号分为低频电压分量和高频电压分量；

所述高频反算模块，用于对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；

所述低频反算模块，用于对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；

所述整合模块，用于将所述一次电压高频电压分量和所述一次电压低频电压分量整合，得到电力系统的一次电压。

优选地，所述高频反算模块还用于将所述高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

所述传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数；

其中，所述高频反算模块包括转换子模块和计算子模块；

所述转换子模块，用于对所述传递函数进行拟合转换，得到传递函数的状态方程；并引入变量x和中心差分法对所述状态方程进行转换，得到离散电压重构函数；

所述计算子模块，用于采用迭代方式对所述离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，所述离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

优选地，所述低频反算模块还用于将所述低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对所述低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量：

其中，所述基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

本发明还提供一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构设备，包括处理器以及存储器；

所述存储器，用于存储程序代码，并将所述程序代码传输给所述处理器；

所述处理器，用于根据所述程序代码中的指令执行上述所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法。

从以上技术方案可以看出，本发明实施例具有以下优点：该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法、装置及设备，该方法步骤包括：通过PT采集电力系统的二次电压信号，并对二次电压信号分为低频电压分量和高频电压分量；对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；将一次电压高频电压分量和一次电压低频电压分量整合，得到电力系统的一次电压。该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法通过对PT采集二次电压信号划分的高频电压分量和低频电压分量分别采用逆黑盒模型和逆电磁对偶模型处理，得到一次电压高频电压分量和一次电压低频电压分量后将其相加得到一次电压，得到的一次电压不受PT采集过程中失真的影响，数据准确率高；解决了现有电力系统采用PT测量电压，在测量过程中PT存在失真情况，导致测量数据不准确的技术问题。

附图说明

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

图1为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法的步骤流程图；

图2为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法信号划分的步骤流程图；

图3为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法逆电磁对偶模型的示意图；

图4为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法与失真信号的对比图；

图5为本发明另一实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法与失真信号的对比图；

图6为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置的框架图。

具体实施方式

为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

本申请的术语解释如下：

电磁式电压互感器是通过变压器实现电磁隔离的电压测量装置，也是一种仪用变压器，基本原理与变压器完全相同。

在线监测指的是在被测设备处于运行的条件下，对设备的状况进行连续或定时的监测，通常是自动进行的。

PT一次侧指的是PT一次侧绕组(高压绕组)与电网直接相连。

PT一次电压指的是PT一次绕组两端电压。

PT二次侧指的是PT二次侧绕组(低压绕组)与计量装置等直接相连。

PT二次信号指的是PT二次绕组两端电压，也是被测量的信号。

黑盒模型指的是一种端口等效模型，不具备物理意义，仅能实现对端口特性与被建模设备一致。PT的黑盒模型在高频特性模拟方面具备显著的准确度，但是在低频特性上由于测量等因素，误差较大。

电磁对偶模型是基于电磁对偶原理推导而来的模型，通过电量及磁量的对偶关系，将设备的磁路模型转换成电路进行表征，具备物理意义。电磁对偶模型依据适用的频率范围可用不同精细度的模型，然而，该电磁对偶模型在高频特性模拟上需要非常复杂的模型拓扑，对参数的准确度的要求十分高，难以通过试验测量得到全部参数，需要设备的详细设计参数。应用于中低频的电磁对偶模型具备很高的精度。

逆模型是与正模型相对，以PT为例，PT的输入是一次电压，输出是二次信号，基于一次电压得到二次信号的模型为正模型，而逆模型的输入是二次信号，输出是实际的一次电压，即基于二次信号输出一次电压的模型及为逆模型。

本申请实施例提供了一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法、装置及设备，应用于配电网的小电流接地系统上，用于解决了现有电力系统采用PT测量电压，在测量过程中PT存在失真情况，导致测量数据不准确的技术问题。

实施例一：

图1为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法的步骤流程图。

如图1所示，本发明实施例提供了一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，方法包括以下步骤：

S1.通过PT采集电力系统的二次电压信号，并对二次电压信号分为低频电压分量和高频电压分量。

需要说明的是，在步骤S1中主要是将PT采集电力系统的二次电压信号划分为低频电压分量和高频电压分量，便于后续步骤对二次电压信号处理转换为一次电压。

S2.对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量。

需要说明的是，在步骤S2中主要是根据步骤S1获得的高频电压分量，再通过逆黑盒模型分析处理，得到一次电压高频电压分量。该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法通过逆黑盒模型实现基于离散二次电压信号数据的一次电压信号的重建，便于得到一次电压高频电压分量，该方式稳定性好。

S3.对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量。

需要说明的是，在步骤S3中主要是根据步骤S1获得的低频电压分量，再通过逆电磁对偶模型分析处理，得到一次电压低频电压分量。该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法通过逆电磁对偶模型实现考虑深度饱和的基于二次电压信号低频分量的一次电压重构，且该逆电磁对偶模型参数获取方法成熟，无需大量现场实测数据进行建模前的训练，计算简单化。

S4.将一次电压高频电压分量和一次电压低频电压分量整合，得到电力系统的一次电压。

需要说明的是，主要是根据步骤S2的一次电压高频电压和步骤S3的一次电压低频电压分量相加整合，得到电力系统的一次电压。该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法对PT采集的二次电压信号进行处理，得到一次电压，避免了PT测量电力系统过程中失真导致测量数据不准确的问题，该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法实现重构PT一次侧的高频暂态，也能重构PT一次侧的低频暂态过电压。

本发明提供的一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，步骤包括：通过PT采集电力系统的二次电压信号，并对二次电压信号分为低频电压分量和高频电压分量；对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；将一次电压高频电压分量和一次电压低频电压分量整合，得到电力系统的一次电压。该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法通过对PT采集二次电压信号划分的高频电压分量和低频电压分量分别采用逆黑盒模型和逆电磁对偶模型处理，得到一次电压高频电压分量和一次电压低频电压分量后将其相加得到一次电压，得到的一次电压不受PT采集过程中失真的影响，数据准确率高；解决了现有电力系统采用PT测量电压，在测量过程中PT存在失真情况，导致测量数据不准确的技术问题。

图2为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法信号划分的步骤流程图。

如图2所示，在本发明的一个实施例中，对二次电压信号分为低频电压分量和高频电压分量包括：

对二次电压信号采用傅里叶变换处理，得到二次电压信号频域形式；

对二次电压信号频域形式中的频率是否大于过渡频率划分为二次电压信号低频频域和二次电压信号高频频域；

分别对二次电压信号低频频域和二次电压信号高频频域采用傅里叶逆变换，得到对应的低频电压分量和高频电压分量。

需要说明的是，主要是将二次电压信号分为低频电压分量及高频电压分量。具体为：FFT将二次电压信号从时域信号转换为频域信号，然后利用过渡频率f_s进行选择，若频率大于过渡频率f_s的为二次电压信号高频频域，小于等于过渡频率f_s的为二次电压信号低频频域量；将频率大于过渡频率f_s的分量集合起来再进行傅里叶逆变换(iFFT)即可得到与二次电压信号对应时域的高频电压分量。同理，将小于等于过渡频率f_s的分量集合起来进行iFFT即可得到与二次电压信号对应时域的低频电压分量。低频电压分量及高频电压分量采用反算过渡频率f_s进行区分。过渡频率f_s取决于PT的电压传递特性的散射参数测量结果，一般过渡频率f_s远小于第一个PT随频率变化的电压传递特性的第一个谐振点，且一般需要小于0.1倍的该频率点频率。由此，过渡频率f_s可以根据需求限定，此处不做限制。

在本发明的一个实施例中，对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量的步骤包括：

将高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

其中，传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数。

在本发明实施例中，通过逆黑盒模型的传递函数对高频电压分量进行重构变换包括：

对传递函数进行拟合转换，得到传递函数的状态方程；

引入变量x和中心差分法对状态方程进行转换，得到离散电压重构函数；

采用迭代方式对离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

需要说明的是，逆黑盒模型主要是将基于二次信号的高频电压分量重构一次电压的高频分量，具体为：通过逆黑盒模型将高频电压分量中得到的数据采用散射矩阵进行化简，即可得到PT一次电压与二次电压信号的电压传递特性，即H_m(s)：

v_sh(s)为二次的高频高压分量，v_ph(s)为一次电压的高频分量，S₁₁、S₁₂、S₂₁、S₂₂均为散射矩阵S的矩阵元素。由此，可得到：v_sh(s)＝H_m(s)v_ph(s)，由于逆黑盒模型是通过二次电压信号重构一次电压的，因此，可改写为：

由此可知，H_m ^-1(s)就是逆黑盒模型的传递函数，通过矢量匹配法拟合H_m ^-1(s)即可得到有理分式形式的H_i(s)：

式中，d为常数项，e为线性项系数，r_k和p_k为频域响应H_i(s)的零点和极点，N为拟合阶数。由有理分式形式的H_i(s)转化为传递函数的状态方程，即是：H_i(s)＝C(sI-A)^-1B+D+Es，式中，I为N*N的标准单位阵，对角元素全为1，其余为0；D和E分别对应于d和e；在散射矩阵中，E通常等于0。得到逆电压传递函数H_i(s)后，可定义一个新的变量x将一次电压的高频分量v_ph及二次信号的高频电压分量v_sh引入。x的定义为：x＝(sI-A)^-1Bv_sh，得到

v_ph＝Cx+Dv_sh。由于采集的二次电压信号是离散信号而非连续数据，而适用于连续数据，因此，采用中心差分法对状态方程进行转换，使其能够适用于离散数据得电压重构的离散电压重构函数，即是：

v_ph(k)＝Cx_k+Dv_sh(k)

对离散电压重构函数进行简化得到：x_k＝αx_k-1+λBv_sh(k)+μBv_sh(k-1)，v_ph(k)＝Cx_k+Dv_sh(k)；其中，

式中，λ、α、μ均为引入的变量，无意义，便于阅读。简化后的状态变量x_k与同一时间点得的输入变量v_sh(k)相关。因此，需要明引入了一个新的状态变量x_k'用以规避迭代计算得矛盾，状态变量x_k'为：x′_k＝x_k-λBv_sh(k)，对应的离散电压重构函数转换为离散状态空间方程(也称逆黑盒模型)，该离散状态空间方程为：

v_ph(k)＝Cx′_k+Gv_sh(k)，

G＝D+CλB，式中的G、

均为引入的变量，无意义，便于阅读；v_ph(k)为重构后输出的一次电压高频电压分量。逆黑盒模型的输入变量为二次信号的高频高压分量，逆黑盒模型的输出变量重构得到得一次电压。

在本发明的一个实施例中，对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量的步骤包括：

将低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量；

其中，基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

图3为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法逆电磁对偶模型的示意图。

需要说明的是，逆电磁对偶模型主要是用于重构一次电压低频电压分量。其由PT的正向电磁对偶模型推导得到，如图3所示的PT低频电磁对偶模型，其中，R_s1和R_s2分别为逆电磁对偶模型一次绕组的电阻和逆电磁对偶模型二次绕组的电阻；L_s为逆电磁对偶模型的漏感。漏感和绕组的电阻是恒定的，而两个磁化电感分别为图3中的L_m1和L_m2且是高度非线性的。L_m1和L_m2与铁芯不同部位的磁导有关。逆电磁对偶模型的分流电阻R_m1和R_m2代表PT的磁芯损耗，它们远大于磁化阻抗。N₀、N₁和N₂分别为参考匝数、PT一次绕组匝数和二次绕组匝数。v_pl、v_sl和i_pl分别为测量的两个绕组两端的端电压和一次电流；R_L和i_s分别为负载和负载电流；i_Ls为流经漏感上电流；v_m1、v_m2、i_m1、i_m2分别为励磁支路1、2的电压和电流；i_L1、i_R1、i_L2、i_R2分别为流经L_m1、R_m1、L_m2、R_m2的电流；i_m1和i_m2分别是流经磁化支路1和2的电流。其中，假设N₀＝N₂，匝数比n＝N₁/N₂。

如图3所示，逆电磁对偶模型通过正向电磁对偶模型中电压电流的关系推导而来，其输入变量为二次信号的低频电压分量，而输出变量为一次电压低频电压分量，其一次电压低频电压分量重构具体为：

负载电流i_L是通过仪器测量，也可以通过负载阻抗和负载两端电压计算得到。因此，磁化支路2两端的电压的计算公式为：v_m2＝v_sl+i_sR_s2，磁链通过对电压积分获得。因此，磁化支路2(λ_m2)上的磁链计算公式为：

式中，t₁为积分的持续时间，λ(0)是λ的初始值，磁化电感的λ-i曲线可以通过空载测试和饱和测试获得。那么流过磁化电阻的电流用式计算，然后流过磁化支路2的电流用式计算。同样，i_L1、i_R1和i_m1可以采用相同的方式计算得到。其中，

i_m2＝i_L2+i_R2

在逆电磁对偶模型中，漏感是常数。流过漏感的电流等于i_m2。因此，漏感两端的磁链和电压(λ_Ls和v_Ls)可以分别通过如下梁式计算，公式为：

i_m2＝i_L2+i_R2，i_Ls＝i_m2+i_s，λ_Ls＝i_LsL_s；

然后分别用和计算λ_m1和v_m1。采用逆对偶导出模型中的磁链守恒定律为：λ_m1＝λ_Ls+λ_m2；v_m1＝v_Ls+v_m2，那么流经一次绕组的电流由i_pl＝(i_m1+i_Ls)/n获得，然后采用v_pl＝nv_m1+R_s1i_pl计算一次电压。式中，λ_m1为逆电磁对偶模型中磁化支路1的磁链，λ_m2为逆电磁对偶模型中磁化支路2的磁链，λ_Ls为逆电磁对偶模型中漏感的磁链。PT的电阻、电感、匝数比通过测试测量或由设备出厂商提供。二次绕组两端的电压由现场的PT测量。因此，可以直接用于重构PT的一次电压低频电压分量。通过使用梯形积分和中心差分方程来转换为和，用于基于离散二次信号的低频电压分量的一次电压低频电压分量重构。

其中Λ(k)、V(k)、V_Ls(k)和I_m2(k)分别是λ(t)、v(t)、v_Ls(t)和i_m2(t)的离散形式。k＝1，2，3，……。

图4为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法与失真信号的对比图，图5为本发明另一实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法与失真信号的对比图。

在本发明实施例中，该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法对低频暂态铁磁谐振过电压工况和雷电冲击工况这两种失真信号的二次电压信号进行重构处理，具体如：在低频暂态铁磁谐振过电压工况中，图4给出了实际一次电压，二次失真信号及重构后一次电压的对比，证明了通过该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法测量的一次电压的低频精度，测量精度满足现场要求。由图4可知，在铁磁谐振下，PT饱和，二次侧信号显著失真，然而，通过该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法能够重构得到一次侧电压，且结果几乎和实际一次侧电压重合。在雷电冲击工况中，图5给出了实际一次电压，二次失真信号及重构后一次电压的对比，证明了该于逆黑盒及逆电磁对偶模型的PT一次电压重构方法测量的一次电压的高频准确性，准确测量精度满足现场要求；由图5可知，在雷电冲击下，PT的频率依赖性使其二次侧信号发生显著失真，通过该基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法能够重构得到准确的一次侧雷电电压，与施加的雷电冲击几乎重合。

实施例二：

图6为本发明实施例所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置的框架图。

如图6所示，本发明实施例还提供一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置，包括：频率分量提取模块10、高频反算模块20、低频反算模块30和整合模块40；

频率分量提取模块10，用于通过PT采集电力系统的，并对二次电压信号分为低频电压分量和高频电压分量；

高频反算模块20，用于对高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；

低频反算模块30，用于对低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；

整合模块40，用于将一次电压高频电压分量和一次电压低频电压分量整合，得到电力系统的一次电压。

在本发明实施例中，高频反算模块20还用于将高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数；

其中，高频反算模块20包括转换子模块和计算子模块；

转换子模块，用于对传递函数进行拟合转换，得到传递函数的状态方程；并引入变量x和中心差分法对状态方程进行转换，得到离散电压重构函数；

计算子模块，用于采用迭代方式对离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

在本发明实施例中，低频反算模块30还用于将低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量；

其中，磁链守恒以及基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

需要说明的是，实施例二装置中的模块对应于实施例一方法中的步骤，实施例一方法中的步骤已在实施例一中详细阐述了，在此实施例二中不再对装置中的模块内容进行详细阐述。

实施例三：

本发明实施例提供了一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构设备，包括处理器以及存储器；

存储器，用于存储程序代码，并将程序代码传输给处理器；

处理器，用于根据程序代码中的指令执行上述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法。

需要说明的是，处理器用于根据所程序代码中的指令执行上述的一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法实施例中的步骤。或者，处理器执行计算机程序时实现上述各系统/装置实施例中各模块/单元的功能。

示例性的，计算机程序可以被分割成一个或多个模块/单元，一个或者多个模块/单元被存储在存储器中，并由处理器执行，以完成本申请。一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述计算机程序在终端设备中的执行过程。

终端设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。终端设备可包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，并不构成对终端设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如终端设备还可以包括输入输出设备、网络接入设备、总线等。

所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

存储器可以是终端设备的内部存储单元，例如终端设备的硬盘或内存。存储器也可以是终端设备的外部存储设备，例如终端设备上配备的插接式硬盘，智能存储卡(SmartMedia Card，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。进一步地，存储器还可以既包括终端设备的内部存储单元也包括外部存储设备。存储器用于存储计算机程序以及终端设备所需的其他程序和数据。存储器还可以用于暂时地存储已经输出或者将要输出的数据。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，包括以下步骤：

通过PT采集电力系统的二次电压信号，并对所述二次电压信号分为低频电压分量和高频电压分量；

对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；

对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；

将所述一次电压高频电压分量和所述一次电压低频电压分量整合，得到电力系统的一次电压。

2.根据权利要求1所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，对所述二次电压信号分为低频电压分量和高频电压分量包括：

对所述二次电压信号采用傅里叶变换处理，得到二次电压信号频域形式；

对所述二次电压信号频域形式中的频率是否大于过渡频率划分为二次电压信号低频频域和二次电压信号高频频域；

分别对所述二次电压信号低频频域和所述二次电压信号高频频域采用傅里叶逆变换，得到对应的低频电压分量和高频电压分量。

3.根据权利要求1所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量的步骤包括：

将所述高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

其中，所述传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数。

4.根据权利要求3所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换包括：

对所述传递函数进行拟合转换，得到传递函数的状态方程；

引入变量x和中心差分法对所述状态方程进行转换，得到离散电压重构函数；

采用迭代方式对所述离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，所述离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

5.根据权利要求1所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量的步骤包括：

将所述低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对所述低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量；

其中，所述基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

6.根据权利要求1所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法，其特征在于，包括：将所述一次电压高频电压分量与所述一次电压低频电压分量相加整合，得到电力系统的一次电压。

7.一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置，其特征在于，包括：频率分量提取模块、高频反算模块、低频反算模块和整合模块；

所述频率分量提取模块，用于通过PT采集电力系统的二次电压信号，并对所述二次电压信号分为低频电压分量和高频电压分量；

所述高频反算模块，用于对所述高频电压分量采用逆黑盒模型进行重构一次电压，得到一次电压高频电压分量；

所述低频反算模块，用于对所述低频电压分量采用逆电磁对偶模型进行重构一次电压，得到一次电压低频电压分量；

所述整合模块，用于将所述一次电压高频电压分量和所述一次电压低频电压分量整合，得到电力系统的一次电压。

8.根据权利要求7所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置，其特征在于，所述高频反算模块还用于将所述高频电压分量作为逆黑盒模型的输入，通过逆黑盒模型的传递函数对所述高频电压分量进行重构变换，逆黑盒模型输出一次电压高频电压分量；

所述传递函数为

式中，v_sh(s)为逆黑盒模型的输入变量，v_ph(s)为逆黑盒模型输出变量，H_m ^-1(s)为逆黑盒模型的传递函数；

其中，所述高频反算模块包括转换子模块和计算子模块；

所述转换子模块，用于对所述传递函数进行拟合转换，得到传递函数的状态方程；并引入变量x和中心差分法对所述状态方程进行转换，得到离散电压重构函数；

所述计算子模块，用于采用迭代方式对所述离散电压重构函数进行计算，得到重构的一次电压高频电压分量；

其中，所述离散电压重构函数为：

v_ph(k)＝Cx_k+Dv_sh(k)

式中，x为引用变量符号，k、k-1分别为高频电压分量中第k个、k-1个的时刻点，A为传递函数极点的N×N对角矩阵，B为N×1数组，Δt为第k个高频电压分量与第k-1个高频电压分量之间的时间间隔，C为传递函数零点的1×N的数组，D为常数项。

9.根据权利要求7所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构装置，其特征在于，所述低频反算模块还用于将所述低频电压分量作为逆电磁对偶模型的输入，通过逆电磁对偶模型的磁链守恒以及基尔霍夫电流电压定律对所述低频电压分量进行重构变换，逆电磁对偶模型输出一次电压低频电压分量；

其中，所述基尔霍夫电流电压定律为：

v_pl＝nv_m1+R_s1i_pl；

v_m1＝v_Ls+v_m2；

式中，v_pl为一次电压低频电压分量，n为逆电磁对偶模型的匝数比，v_m1为逆电磁对偶模型中第一励磁支路的电压，v_m2为逆电磁对偶模型中第二励磁支路的电压，v_Ls为逆电磁对偶模型中漏感两端的电压，R_s1为逆电磁对偶模型中一次绕组的电阻，i_pl为逆电磁对偶模型的一次电流，i_m1为逆电磁对偶模型中流过第一励磁支路的电流，i_Ls为流过逆电磁对偶模型中漏感两端的电流。

10.一种基于逆黑盒及逆电磁对偶模型的PT一次电压重构设备，其特征在于，包括处理器以及存储器；

所述存储器，用于存储程序代码，并将所述程序代码传输给所述处理器；

所述处理器，用于根据所述程序代码中的指令执行如权利要求1-6任意一项所述的基于逆黑盒及逆电磁对偶模型的PT一次电压重构方法。

Images