(12)发明专利申请
(10)申请公布号 CN 109902735 A(43)申请公布日 2019.06.18
(21)申请号 201910132876.4(22)申请日 2019.02.22
(71)申请人 武汉格蓝若智能技术有限公司
地址 430000 湖北省武汉市东湖新技术开
发区光谷大道303号光谷.芯中心2-07栋1803-1805室(72)发明人 窦峭奇
(74)专利代理机构 北京开林佰兴专利代理事务
所(普通合伙) 11692
代理人 刘帅帅(51)Int.Cl.
G06K 9/62(2006.01)G06F 17/16(2006.01)
权利要求书2页 说明书5页 附图1页
(54)发明名称
一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法
(57)摘要
本发明公开了一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，通过以下步骤实现：（1）采集校准后刚投运的三相电压互感器测量信号，提取其基波幅值构成测量矩阵；（2）利用独立成分分析法对上述测量矩阵进行分解，获取独立分量，以及对应的解混矩阵；（3）利用二范数法筛选出残余分量，同时得到计算残余分量的参数矩阵，并计算其控制限；（4）采集三相电压互感器输出电压幅值，依照（2）与（3）计算其对应统计量，并与控制限进行比较，若超出控制限，则判定有电压互感器超差。优点：本方法解决了无标准器条件下的三相电压互感器在线评估问题，评估过程无需停电、不依赖高精度标准器，评估流程简便、效率更高。
CN 109902735 ACN 109902735 A
权 利 要 求 书
1/2页
1.一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，其特征在于，包括以下步骤：
(1)采集校准后的三相电压互感器测量信号，提取其基波幅值构成测量矩阵
(2)利用独立成分分析法对上述测量矩阵进行分解，获取独立分量S1，S2和S3，以及对应的解混矩阵W；
(3)利用二范数法筛选出残余分量Se，同时得到计算残余分量的参数矩阵We，建立反映残余分量变化的统计量SPE，并计算其控制限；
(4)采集三相电压互感器输出电压幅值，依照(2)与(3)计算其对应统计量，并与控制限进行比较，若超出控制限，则判定三相电压互感器中有电压互感器超差。
2.一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，其特征在于：其步骤(2)中基于独立成分分析方法的测量矩阵分解过程分为2.1、2.2和2.3，具体步骤如下：
2.1)对测量矩阵X白化，获得白化矩阵Z，Z中各列向量特征为：零均值、单位方差且互不相关；
2.1.1)对测量矩阵X中心化，使各变量均值为零，即元素全为1，
2.1.2)获得矩阵
的协方差矩阵
为X各行均值；
其中，I中
2.1.3)对协方差矩阵进行特征值分解，即RX＝EDET，其中E为特征向量，D为特征值对角矩阵；
2.1.4)白化矩阵
其中Q＝D1/2UT；
2.2)估计解混矩阵W；
2.2.1)确定待估计的独立成分个数为3，置计数变量p＝1；2.2.2)生成单位范数的随机变量wp；2.2.3)
对w′得到p进行正交化，
其中E{·}为求期望，g(x)＝tanh x，再对w″得到w″p进行标准化，p＝w″p/
||w″p||；
2.2.4)若w″则继续执行步骤(2.2.3)。若收敛，则置p＝p+1；p未收敛，2.2.5)若p≤3，则返回步骤(2.2.2)，若p＞4，则W＝[w1 w2 w3]T；
2.3)计算独立分量
3.一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，其特征在于：步骤(3)的具体步骤如下：
2
CN 109902735 A
权 利 要 求 书
2/2页
3.1)计算解混矩阵W中每行向量wi对应的二范数占比根据
占比排序，选择占比最低的行向量为残余分量对应行向量作为解混矩阵We，该行向量对应独立成分作为残余分量Se；
3.2)构建反映残余分量的统计量SPE＝Q-1WeSe；3.3)采用核密度估计方法，计算统计量SPE对应控制限；3.3.1)记f(x)为统计量SPE的频次统计函数，其自变量x为统计量SPE的幅值，因变量f(x)为对应频次，采用核密度估计方法估计SPE的概率密度函数
估计方法为
其中，n为样本容量，h为带宽，优选的，h＝1.0592n-0.2，K为核函
数，优选的，采用标准高斯函数；
3.3.2)记
为统计量SPE的概率分布函数，α为置信度，优选的，α＝
0.99，则统计量SPE在置信水平α下的控制限为xα＝F-1(α)；进一步、统计量控制限计算方法。
4.一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，其特征在于，步骤(4)的具体步骤如下：
4.1)在线采集三相电压互感器的输出电压幅值，记t时刻的三相电压幅值为V(t)＝[Va Vb Vc]T，对其进行中心化，得到测量数据
已在步骤(2.1.1)定义。4.2)计算其对应SPE统计量为
4.3)比较SPE(t)与xα大小关系，若SPE(t)＞xα，则发出告警信号，认为三相电压互感器已出现了超差故障，需要检修。若SPE(t)＜xα，则认为三相电压互感器并未发生超差故障，可继续运行。
3
CN 109902735 A
说 明 书
1/5页
一种实现变电站内三相电压互感器运行误差评估的二范数独
立成分分析方法
技术领域
[0001]本发明属于电压互感器技术领域，具体涉及一种变电站内三相电压互感器运行误差的评估方法。
背景技术
[0002]电压互感器是一种广泛应用于电力系统的高压信号传感设备，是电力系统测控、保护和计量等设备所用电压信号的来源，准确的电压测量是电力系统安稳运行的基础，因此，评估电压互感器的测量误差并及时更换超差的不合格电压互感器非常重要。不仅如此，随着智能电网的建设，在运电压互感器的数量快速增长，对评估效率也提出了更高的要求。[0003]目前广泛采用的电压互感器误差评估方法被称之为周检，即在变电站现场周期性开展电压互感器误差的离线检测。周检有以下缺点：1)离线检测需要变电站配合停电，损害了供电可靠性；2)周期性检测无法第一时间发现超差故障，导致超差互感器的更换滞后；3)周检没有针对性，评估效率低。[0004]综上所述，传统的电压互感器周检技术已无法适应智能电网快速发展下的电压互感器运行误差评估需求，为了实现对电压互感器运行误差的在线、长期以及高效评估，本发明提出了一种基于二范数独立成分分析的变电站内三相电压互感器运行误差评估方法。发明内容
[0005]针对现有周检的技术缺陷及改进需求，本发明提供了一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，目的在于脱离标准电压互感器条件下的电压互感器运行误差评估。[0006]为实现上述目的，本发明提供的一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，包括以下步骤：[0007](1)采集校准后刚投运的三相电压互感器测量信号，提取其基波幅值构成测量矩阵
其中，ua是指A相电压互感器输出信号的基波幅值，ub是指B相电压互感器输出信号的基波幅值，uc是指C相电压互感器输出信号的基波幅值，n为采集样本数量；[0009](2)利用独立成分分析法(Independent Component Analysis,ICA)对上述测量矩阵进行分解，获取独立分量S1，S2和S3，以及对应的解混矩阵W；[0010](3)利用二范数法筛选出残余分量Se，同时得到计算残余分量的参数矩阵We，建立反映残余分量变化的统计量SPE，并计算其控制限；[0011](4)采集三相电压互感器输出电压幅值，依照(2)与(3)计算其对应统计量，并与控制限进行比较，若超出控制限，则判定三相电压互感器中有电压互感器超差，应组织检修；
[0008]
4
CN 109902735 A[0012]
说 明 书
2/5页
进一步、上述电压互感器误差状态评估方法，其步骤(2)中基于独立成分分析方法
的测量矩阵分解过程如下：
[0013]2.1)对测量矩阵X白化，获得白化矩阵Z，Z中各列向量特征为：零均值、单位方差且互不相关；
[0014]2.2)估计解混矩阵W；
[0015][0016][0017]
2.3)计算独立分量
进一步、上述电压互感器误差状态评估方法，其步骤(2.1)包括以下子步骤：2.1.1)对测量矩阵X中心化，使各变量均值为零，即
为X各行均值。
的协方差矩阵
其
中，I中元素全为1，
[0018][0019]
2.1.2)获得矩阵
2.1.3)对协方差矩阵进行特征值分解，即RX＝EDET，其中E为特征向量，D为特征值
对角矩阵；
2.1.4)白化矩阵
其中Q＝D1/2UT；
[0020][0021][0022][0023][0024]
进一步、上述电压互感器误差状态评估方法，其步骤(2.2)包括以下子步骤：2.2.1)确定待估计的独立成分个数为3，置计数变量p＝1；2.2.2)生成单位范数的随机变量wp；2.2.3)
其中E{·}为求期望， g(x)＝
再对w′得到p进行标准化，
tanhx，对w′得到p进行正交化，
w″p＝w″p/||w′p||；
[0025]2.2.4)若w″则继续执行步骤(2.2.3)。若收敛，则置 p＝p+1；p未收敛，[0026]2.2.5)若p≤3，则返回步骤(2.2.2)，若p＞4，则 W＝[w1 w2 w3]T。[0027]进一步、上述电压互感器误差状态评估方法，其步骤(3)包括以下子步骤：
[0028]
3.1)计算解混矩阵W中每行向量wi对应的二范数占比
根据占比排序，选择占比最低的行向量为残余分量对应行向量作为解混矩阵We，该行向量对应独立成分作为残余分量 Se；
[0029]3.2)构建反映残余分量的统计量SPE＝Q-1WeSe；[0030]3.3)采用核密度估计方法，计算统计量SPE对应控制限。[0031]进一步、上述电压互感器误差状态评估方法，其步骤(3.3)包括以下子步骤：[0032]3.3.1)记f(x)为统计量SPE的频次统计函数，其自变量x为统计量SPE的幅值，因变量f(x)为对应频次。采用核密度估计方法估计 SPE的概率密度函数
估计方法为
其中，n为样本容量，h为带宽，优选的，h＝1.0592n-0.2，K为核函
数，优选的，采用标准高斯函数。
5
CN 109902735 A[0033]
说 明 书
3/5页
3.3.2)记为统计量SPE的概率分布函数，α为置信度，优选
的，α＝0.99，则统计量SPE在置信水平α下的控制限为xα＝F-1(α)；
[0034]进一步、统计量控制限计算方法：[0035]进一步、上述电压互感器误差状态评估方法，其步骤(4)包括以下子步骤：[0036]4.1)在线采集三相电压互感器的输出电压幅值，记t时刻的三相电压幅值为V(t)＝[Va Vb Vc]T，对其进行中心化，得到测量数据
已在步骤(2.1.1)定义；
[0037][0038]
4.2)计算其对应SPE统计量为
4.3)比较SPE(t)与xα大小关系，若SPE(t)＞xα，则发出告警信号，认为三相电压互感器已出现了超差故障，需要检修。若SPE(t)＜xα，则认为三相电压互感器并未发生超差故障，可继续运行，优选的，在线采集数据之间的采样间隔为15min。[0039]本发明的有益效果是：
[0040]本发明首先利用校准后刚投运的三相电压互感器采集数据建立了误差正常状态下的特征量模型，并估计其概率分布，将高置信度下的特征量数值作为控制限。其次，采集该组三相电压互感器在线运行过程中的输出电压幅值，并采用相同的步骤计算其特征量，最后，通过比较在线采集的特征量数值与控制限大小关系，判断三相电压互感器是否出现了超差故障，从而实现了脱离高精度标准器条件下的三相电压互感器运行误差的实时评估。
[0041]本发明采用了基于二范数的独立成分分析方法，由于该方法对输入变量的概率分布并无特殊要求，对电力系统一次电压幅值的概率分布并不提出特殊要求，因此适用面广。[0042]本发明提出的三相电压互感器运行误差评估方法同传统方法相比，减少了现场工作量，不需要停电并且实现了实时评估，在此方法的辅助下，检修人员能够有针对性的对超差互感器进行检修，减少了大量人力物力的投入。附图说明
[0043]为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体的实施方式或现有技术描述中所需要使用的附图作简单地介绍，附图中，各元件或部分并不一定按照实际的比例绘制。
[0044]图1为本发明的分析流程示意图。
具体实施方式
[0045]下面将结合附图对本发明技术方案的实施例进行详细的描述，以下实施例仅用于更加清楚地说明本发明的技术方案，因此只作为示例，而不能以此来限制本发明的保护范围。
[0046]本发明的一种实现变电站内三相电压互感器运行误差评估的二范数独立成分分析方法，流程如图1所示，包括以下步骤：[0047](1)采集校准后刚投运的三相电压互感器测量信号，提取其基波幅值构成测量矩
6
CN 109902735 A
说 明 书
4/5页
阵
其中，ua是指A相电压互感器输出信号的基波幅值，ub是指B相电压互感器输出信号
的基波幅值，uc是指C相电压互感器输出信号的基波幅值，n为采集样本数量；[0049](2)利用独立成分分析法(Independent Component Analysis，ICA)对上述测量矩阵进行分解，获取独立分量S1，S2和S3，以及对应的解混矩阵W；[0050](3)利用二范数法筛选出残余分量Se，同时得到计算残余分量的参数矩阵We，建立反映残余分量变化的统计量SPE，并计算其控制限；[0051](4)采集三相电压互感器输出电压幅值，依照(2)与(3)计算其对应统计量，并与控制限进行比较，若超出控制限，则判定三相电压互感器中有电压互感器超差，应组织检修；[0052]本实施例中、上述电压互感器误差状态评估方法，其步骤(2)中基于独立成分分析方法的测量矩阵分解过程如下：[0053]2.1)对测量矩阵X白化，获得白化矩阵Z，Z中各列向量特征为：零均值、单位方差且互不相关；
[0054]2.2)估计解混矩阵W；
[0055][0056][0057][0048]
2.3)计算独立分量
本实施例中、上述电压互感器误差状态评估方法，其步骤(2.1)包括以下子步骤：2.1.1)对测量矩阵X中心化，使各变量均值为零，即
为X各行均值；
其中，
I中元素全为1，
[0058][0059]
2.1.2)获得矩阵的协方差矩阵
2.1.3)对协方差矩阵进行特征值分解，即RX＝EDET，其中E为特征向量，D为特征值
对角矩阵；
2.1.4)白化矩阵
其中Q＝D1/2UT。
本实施例中、上述电压互感器误差状态评估方法，其步骤(2.2)包括以下子步骤：2.2.1)确定待估计的独立成分个数为3，置计数变量p＝1；2.2.2)生成单位范数的随机变量wp；2.2.3)
g(x)＝tanh x，对w′得到p进行正交化，
其中E{·}为求期望，
再对w′p进
[0060][0061][0062][0063][0064][0065]
行标准化，得到w″p＝w″p/||w′p||；
[0066]2.2.4)若w″则继续执行步骤(2.2.3)。若收敛，则置 p＝p+1；p未收敛，[0067]2.2.5)若p≤3，则返回步骤(2.2.2)，若p＞4，则 W＝[w1 w2 w3]T。[0068]本实施例中、上述电压互感器误差状态评估方法，其步骤(3)包括以下子步骤：
[0069]
3.1)计算解混矩阵W中每行向量wi对应的二范数占比
根据占比排序，选择占比最低的行向量为残余分量对应行向量作为解混矩阵We，该行向量对应独立成分作为残余分量 Se；
7
CN 109902735 A[0070][0071][0072][0073]
说 明 书
5/5页
3.2)构建反映残余分量的统计量SPE＝Q-1WeSe；3.3)采用核密度估计方法，计算统计量SPE对应控制限。本实施例中、上述电压互感器误差状态评估方法，其步骤(3.3)包括以下子步骤：3.3.1)记f(x)为统计量SPE的频次统计函数，其自变量x为统计量SPE的幅值，因变
估计方法为
量f(x)为对应频次。采用核密度估计方法估计 SPE的概率密度函数
其中，n为样本容量，h为带宽，优选的，h＝1.0592n-0.2，K为核函
数，优选的，采用标准高斯函数；
[0074]
3.3.2)记为统计量SPE的概率分布函数，α为置信度，优选
的，α＝0.99，则统计量SPE在置信水平α下的控制限为xα＝F-1(α)；
[0075]本实施例中、上述电压互感器误差状态评估方法，其步骤(4)包括以下子步骤：[0076]4.1)在线采集三相电压互感器的输出电压幅值，记t时刻的三相电压幅值为V(t)＝[Va Vb Vc]T，对其进行中心化，得到测量数据
已在步骤(2.1.1)定义；
[0077][0078]
4.2)计算其对应SPE统计量为
4.3)比较SPE(t)与xα大小关系，若SPE(t)＞xα，则发出告警信号，认为三相电压互感器已出现了超差故障，需要检修。若SPE(t)＜xα，则认为三相电压互感器并未发生超差故障，可继续运行。优选的，在线采集数据之间的采样间隔为15min。[0079]最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。
8
CN 109902735 A
说 明 书 附 图
1/1页
图1
9
因篇幅问题不能全部显示，请点此查看更多更全内容