基于ica算法的智能变电站设备状态预警方法

基于ica算法的智能变电站设备状态预警方法
【专利摘要】本发明提供了一种基于ICA算法的智能变电站设备状态预警方法,包括以下步骤：（一） 信息数据采集；（二）信息数据传送；（三）信息数据处理：控制室的控制中心实时读取采集器传送来的数据集，并利用ICA算法对采集的数据集进行清理、数据集成和转换、数据降维，数据离散化等工作；（四）预警诊断：站控制室的控制中心利用贝叶斯算法对步骤（三）中处理好的数据进行分析，通过数据分类、关联和预测的方法，对历史数据训练建模以及集成工作人员现有的经验得到策略，输出相应的预警信号。本发明能够简单、快速地提取有效数据信息，对变电站一次/二次设备的运行状况作出更精准的预警，令对变电站的监测预警更及时、准确。本发明适用于电力系统所有通信设备的风险评估与预警。
【专利说明】
基于ICA算法的智能变电站设备状态预警方法
技术领域
[0001] 本发明属于电力系统的通信领域，具体地说是一种基于ICA算法的智能变电站设 备状态预警方法。
【背景技术】
[0002] 随着变电站智能化的不断发展，无人值守变电站逐渐增多，对于变电站内高压开 关柜、母线接头、电缆接头、输变电线路接头、室外刀闸开关等重要设备在线实时监测的任 务也不断加重，这些在长期运行过程中容易因绝缘老化或接触电阻过大而发热，从而导致 火灾和大面积的停电事故，造成重大的经济损失及严重的社会影响，因此对智能变电站设 备进行温度实时监测是十分必要的。
[0003] 由于变电站设备触头因为接触不良、氧化、老化等原因都会导致电阻变大，出现突 然高热，触点温度升高是设备故障最直观的表现，故障发生前，一般触点温度都是保持一 致，当负载增加就会促使设备温度变化，使其突然升高，因此现有技术中的变电站设备预警 系统主要是对变电站各设备的温度进行预警，而现有的变电站设备预警系统对变电站各种 重要设备的温度监测时，直接采集个设备的温度信息，然后直接对温度信息进行分析处理， 这种处理方式的弊端在于:变电站电磁环境复杂等因素，采集的同一类触点的温度数据其 实为一个数据的集合，其中包括原始温度数据、其他设备电磁噪声干扰数据、旁路触点温度 干扰数据等集合，在预警系统中做数据分析时如果不将这些无效的数据剔除，极易造成错 误的分析结果，而使预警系统做出错误的预判，会极大的影响变电站设备的正常运行。

【发明内容】

[0004] 本发明要解决的技术问题，是提供一种基于ICA算法的智能变电站设备状态预警 方法，通过对大量的采集数据进行精简筛选，获取出有用的信息，令分析出的结果更加准 确。
[0005] 为解决上述技术问题，本发明所采取的技术方案是： 一种基于ICA算法的智能变电站设备状态预警方法， 它包括以下步骤： (一） 信息数据采集:设于变电站内的无线监测系统通过贴于一次/二次设备附近的温 度传感器采集各类设备的温度信息； (二) 信息数据传送:采集的温度信息通过微功率无线通信上传至与变电站最近的采集 器，再由采集器通过集中器网关、经过光纤或以太网传送至站内控制室的控制中心； (三) 信息数据处理:控制室的控制中心实时读取采集器传送来的数据集，并利用ICA算 法对采集的数据集进行清理、数据集成和转换、数据降维，数据离散化等工作； (四） 预警诊断:站控制室的控制中心利用贝叶斯算法对步骤(三）中处理好的数据进行 分析，通过数据分类、关联和预测的方法，对历史数据训练建模，以及集成工作人员现有的 经验得到策略，输出相应的预警信号。
[0006] 作为对步骤(三)的限定:所述步骤(三)包括以下步骤： 1) 数据预处理:控制室的控制中心对接收到的观测数据进行去均值或加上均值向量后 得到正确的源信号波形估计;然后利用白化法去除数据间的相关性； 2) 进入ICA主算法:控制室的控制中心利用傅里叶算法以负熵最大作为搜寻方向对预 处理后的数据进行快速收敛，并通过非高斯性度量来表示分离结果间的相互独立性，当非 尚斯性达到最大时，完成ICA分尚。
[0007] 作为对步骤1)的限定:所述步骤1)对数据预处理后得到的数据信号集必须满足每 个采集点的温度值之间是相互独立的、得到的分离矩阵为满秩矩阵、得到的数据集中至多 只有一个信号为高斯分布函数三个先验条件。
[0008] 作为对步骤2)的限定:所述步骤2)中利用傅里叶算法进行ICA分离包括依次进行 的以下步骤： ① 随机选择初始权值向量； ② 利用公另
更新权值向量； ③ 利用公式Wk+l:Wk+l/| |wk+l| I进行归一化； ④ 如果I Wk+1ik | >雾;，则算法不收敛，返回步骤②，否则利用傅里叶算法的ICA算法估算 出了 一个独立分量，完成了所采集数据集的提取。
[0009] 作为对步骤（四）的限定:所述步骤（四）中利用的贝叶斯算法包括依次进行的以下 步骤： i)建立温度模型； ? )预处理模型参数； m)参数带入模型，进行故障预测； iv )在线采集温度，更新温度； V)得到此时刻温度分布； Vi)结合此时刻温度分布、温度总体分布、温度先验分布，依据观测值预测此时刻后验 分布。
[0010] 作为对步骤m)的限定:所述步骤中利用贝叶斯算法进行预测的过程为:设温度的 初始状态的概率分布和k-l时刻的概率分布分别为：/K'A)& )、bvD，Xk只与 Xk-1有关，得到：
式中不包含yk的先验分布，由系统的状态转移概率分布来计算； 所述步骤vi)更新后验分布参数的过程为：利用最新观测变量yk和先验分布 戶摊导多由贝叶斯公式可得：
由于
：巴上述两式代入贝叶斯公式， 可得：
再根据条件概率分布、联合概率分布和贝叶斯公式
[0011] 由于采用了上述的技术方案，本发明与现有技术相比，所取得的技术进步在于： (1) 本发明采用无线微功率方式进行信息采集，高可靠、低功耗、使用方便灵活的通信 技术，能够在智能变电站复杂环境下实现无线数据的稳定可靠传输； (2) 本发明采用ICA算法进行有用信息的提取，算法收敛速度快，鲁棒性高，具有并行 性、分布式、计算简单、内存消耗少等优点； (3) 本发明采用贝叶斯算法考虑了主观因素和先验数据的关联性，使算法推理出的结 果具有更加精准的预判性； (4) 本发明采用变电站预警指标体系拥有健壮的系统架构，健全的综合指标体系，利用 贝叶斯算法制定风险评估规则，得到较客观的预警信号指示，并给出合理的解决措施方案， 系统具有预判精准，措施合理等优点。
[0012] 综上所述，本发明能够简单、快速地提取有效数据信息，对变电站一次/二次设备 的运行状况作出更精准的预警，令对变电站的监测预警更及时、准确。
[0013] 本发明适用于电力系统所有通信设备的风险评估与预警。
[0014] 本发明下面将结合具体实施例作进一步详细说明。
【附图说明】
[0015] 图1为本发明实施例的流程图； 图2为本发明实施例中ICA算法的分离模型框图； 图3为本发明实施例中贝叶斯算法的故障预测流程图； 图4a为本发明实施例去燥前混合数据的信号波形图； 图4b为图4a的去燥分离后的信号波形图； 图5为本发明实施例的指标评价体系图。
【具体实施方式】
[0016] 实施例基于ICA算法的智能变电站设备状态预警方法 本实施例提供了一种基于ICA算法的智能变电站设备状态预警方法，如图1所示，它包 括以下步骤： (一）信息数据采集:设于变电站内的无线监测系统通过贴于一次/二次设备附近的温 度传感器采集各类设备的温度信息。本实施例中采用现有技术中常用的温度传感器进行温 度采集即可，无需作出更换，节省成本。
[0017] (二)信息数据传送:采集的温度信息通过微功率无线通信上传至与变电站最近的 采集器，再由采集器通过集中器网关、经过光纤或以太网传送至站内控制室的控制中心。本 实施例中的微功率无线通信设备采用现有的无线通信协议，可将采集的数据信息无线传送 至控制中心，通过控制中心进行完整的处理。
[0018] (三)信息数据处理:控制室的控制中心实时读取采集器传送来的数据集，并利用 ICA算法对采集的数据集进行清理、数据集成和转换、数据降维，数据离散化等工作。
[0019]本实施例如图2所示，采用现有技术中ICA算法的基本模型，即X= AS 用矩阵形式表示为：
式中A为mX η阶满混合矩阵;S为η维随机列向量，即为某时亥Ijn个源信号构成的向量;X 为观测信号或混合信号向量，每个混合信号都是一个随机信号，其每个值都是对随机信号 的一次抽样。由矩阵式可以看出，某个时刻每个观测信号都是由各个源信号经过不同的 加权得到的。要想从观测信号X中得到源信号S的估计Υ，需要一个解混矩阵W进行变换对X，W 可以看成是混合矩阵的一个逆估计，具体过程由Y=WX表示。
[0020] ICA方法的估计就是如何得到最优的解混矩阵W，这样才能使分离的信号最接近源 信号。因此需要针对不同的观测信号，选择合适的目标函数和优化算法去训练W，使其达到 极大或极小，这个#'就是所要求的数据。
[0021] 在对信号进行ICA方法处理之前，先对采集的信号进行一些预处理，预处理方法主 要有中心化和白化，中心化的最基本方法就是去均值，其基本思想是从观测信号X中减去信 号的均值向量，使其变为均值为零的变量，定义式如下：
式中E[x]为X的均值。中心化不会影响信号继续分离的过程，但也不一定所有ICA处理 都需要进行中心化，如果X与S是线性相关的，由式子召如卜0!可知S也是零均值的，所以分 离出的信号与源信号在幅度上有偏差，因此分离出的信号之后，需要再加上均值向量，才会 显示正确的源信号波形的估计。
[0022] 白话主要是去除数据间的相关性，达到空域解相关的目的。此时本实施例的任务 是要找到一种线性变换，使得X变为白化向量： 如果一个零均值向量X的协方差矩阵是单位矩阵，即：
称之为为白化信号。本实施例的处理过程就是要找到这样的线性变换中的白化矩阵z， 可得：
式中，如:麵〖％ ::,?提X协方差矩阵分解后特征值从大到小排列组成的;V是这些特征值 对应的特征向量组成的正交矩阵。白化后，混合矩阵A也转化为新的矩阵重，证明过程如下：
白化后就减少了 ICA需要估计的参数，加快收敛的速度，也加强了算法的稳定性。
[0023] 预处理之后就进入ICA主算法，本发明采用ICA演进的一种快速ICA算法，利用傅里 叶算法进行快速收敛，简称FICA。算法是以负熵最大作为一个搜寻方向，源信号比观测信号 具有更强的非高斯性，在分离过程中，需要进行分离的数据集具有三个先验条件:每个采集 点的温度值之间是相互独立的；得到的分离矩阵为满秩矩阵;得到的数据集中至多只有一 个信号为高斯分布函数，并通过分离的非高斯性度量来表示分离结果间的相互独立性，当 非高斯性达到最大时，表示已完成ICA分离。
[0024] 在ICA分离过程中，有线性可逆变换Y=WTX，WTX的负熵最大可近似通过E{g(WTX)} 来进行优化获取。因为W为正交矩阵，所以|li丨~l，E{g(WTX)}的最优值能在满足公式 的点上获得。
[0025] 上述公式中，雜一个恒定值，|，__峰，喂优化后的r值。通过牛顿迭代法来 解决这个等式，上述公式左边的函数用f表示，可以得到它的Jacobian矩阵:沢趵:::E to Yiii·'％ -潭。 由于数据被球化，即?，所以，嗯淚。因而雅可比矩阵变 成了对角阵，并且能比较容易地求逆。因而可以得到下面的近似牛顿迭代公式：
上述公式中，麽5是_的新值，，规格化能提高解的稳定性。简化后就可 以得到FastICA算法的迭代公式：
实际上，FICA中的期望必须由它们的估计值所取代，当然最好的估计是相应的样本均 值。虽然所有的数据应该利用到，但是这并不是一种好方法，因为计算量可能要求太大，所 以通常是用较少的样本估计出平均数，样本的数量在最后的估计上可能有不可忽视的影 响。如果不满足收敛条件，可以增加样本的数量。
[0026] FICA算法的具体步骤描述如下： ① 随机选择初始权值向量w0，k=0; ② 利用3
更新权值矢量Wk+1; ③ 利用公式Wk+l:Wk+l/| |wk+l| I进行归一化； ④ 如果|wk+1-wk|>竄，算法不收敛，返回②，否则利用傅里叶算法的ICA算法估算出了一 个独立分量，完成了所采集数据集的提取。
[0027] 上述公式中釔为w的参数，w为解混矩阵。
[0028] 如图4所示，为本实施例对某一变电站设备温度采集信息的去燥分离图，图4可知， 由于智能变电站内监测设备主要分为一次设备和二次设备，一次设备主要包括主变压器及 其附属设备、GIS设备、开关柜设备、接地变压器、站用变压器、动态无功补偿装置，其他附属 设备如中性点接地隔离开关、中性点避雷器、零序电流互感器等。二次设备有主要包括综合 自动化设备，如线路保护测控柜，主变压器保护测控柜，电能计量屏，频率电压紧急控制装 置，电能质量监测柜，二次安防设备等等，大量温度传感器探头布置于这些设备表面，进行 温度感应。每个传感器里都内嵌了无线通信模块，借助微功率无线通信技术，与距离最近的 采集器通信，不同的采集器之间通过无线mesh拓扑方式和网络层集中器网关相连，网关通 过光纤以太网方式和控制室内的控制中心相连，上传采集的实时温度数据到控制中心。
[0029] 而本实施例中的智能变电站控制室的控制中心由大量数据库和服务器组成，应用 层汇聚信息预处理软件及设备温度监测预警软件。信息预处理软件主要负责完成ICA信号 的主成分分析提取。由于条件限制，本实施例为了验证其算法的可行性，采用计算机模拟产 生的温度信号。温度传感器按一定得时间间隔（500s)在线采样设备触点模拟温度，应用幅 移键控(ASK)来模拟温度信号进行试验。以1 s为采样频率，采样点为500，取两个监测点1，2 数据为原始信号，分离时以监测点1为主成分，监测点2视为干扰信号，试验环境在在 Windows XP下，用matlab7.1.0运行仿真实验，混合数据图形如图4a所示，设备温度的范围 在-10°C - 110°C之间跳变。经过ICA算法后，可以将两个监测点温度还原，在分离之前先进 行预处理，去均值和白话处理，是从观测信号中减去信号的均值向量从而使X成为一个 均值为零的变量，即对它进行中心化。当然，并不是所有的数据都要进行中心化处理，因为X 与y是线性关系，所以y也是零均值的。接下来用FICA算法对信号进行训练，分离后的图形 如图4b所示，很明显，上面的为监测点1的温度信号，下面的为高斯白噪声信号监测点2的信 号。虽然分离出的信号和源信号有一定差别，但是从波形整体形状上来看是大致一致的，说 明ICA算法分离的健壮性很强，使图形吻合度提高 (四）预警诊断:站控制室的控制中心利用贝叶斯算法对步骤(三）中处理好的数据进行 分析，通过数据分类、关联和预测的方法，对历史数据训练建模以及集成工作人员现有的经 验得到策略，输出相应的预警信号。
[0030] 预警诊断方法是根据统计学原理，相邻观测值之间的依赖性，揭示序列的结构与 规律。对采集的温度数据进行预处理，主要完成了数据清理、数据集成和转换、数据降维、数 据离散化等工作，进而通过贝叶斯算法对数据进行分析，通过数据分类、关联和预测的思想 方法，对历史数据训练建模以及集成现有的经验得到策略。该历史数据为已知类型标号、含 有当前运行状态信息的数据，类型标号值等。数据分析对含有当前运行状态信息的实时数 据流进行分析，得到类型标号，如果类型标号显示为某种故障，则发出报警。本实施例采用 贝叶斯算法对处理好的原始数据进行关联预测，所述的贝叶斯算法是基于统计学的一种数 据挖掘算法，其核心观点是:在关于事件A的任何统计推断问题中，除了使用样本信息X所提 供的信息数据除外，还必须对事件A预先预定一个先验分布，它可以是预测者根据自己的经 验来确定，也可以是预测者根据自己的主观认识来确定。在进行统计推断时，先验概率分布 是不可或缺的一个重要要素。图3为贝叶斯算法的故障预测流程图，包括以下步骤： i)建立温度模型。
[0031] ii)预处理模型参数。
[0032] m)参数带入模型，进行故障预测。本步骤中利用贝叶斯公式将主观认识与先验信息 综合，得到后验信息，后验信息不仅体现主观认识还有先验信息。而利用贝叶斯算法进行预测 的过程为:设温度的初始状态的概率分布和k-Ι时刻的概率分布分别为： Xk只与Xk-1有关，得到：
式中不包含yk的先验分布，由系统的状态转移概率分布梦(％|知以来计算； 所述步骤vi)更新后验分布参数的过程为：利用最新观测变量yk和先验分布 M々hn雕导/>(々h.. J，由贝叶斯公式可得：
由于 . ,
、杷上述两式代入贝叶斯公式，可得：
再根据条件概率分布、联合概率分布和贝叶斯公式
代入雜努難)可得至
[0033] iv )在线采集温度，更新温度； V)得到此时刻温度分布； vi)结合此时刻温度分布、温度总体分布、温度先验分布，依据观测值预测此时刻后验 分布。
[0034]图5所示为本实施例的指标评价体系图，即经过本实施例计算处理后各个设备的 温度信息与其相对应的实际风险。根据图5给出的信息，则可以利用本实施例判断变电站各 个设备的运行情况，并采取相应的措施即可。
【主权项】
1. 一种基于ICA算法的智能变电站设备状态预警方法，其特征在于:它包括以下步骤： (一） 信息数据采集:设于变电站内的无线监测系统通过贴于一次/二次设备附近的温 度传感器采集各类设备的温度信息； (二) 信息数据传送:采集的温度信息通过微功率无线通信上传至与变电站最近的采集 器，再由采集器通过集中器网关、经过光纤或以太网传送至站内控制室的控制中心； (三) 信息数据处理:控制室的控制中心实时读取采集器传送来的数据集，并利用ICA算 法对采集的数据集进行清理、数据集成和转换、数据降维，数据离散化等工作； (四） 预警诊断:站控制室的控制中心利用贝叶斯算法对步骤(三）中处理好的数据进行 分析，通过数据分类、关联和预测的方法，对历史数据训练建模，以及集成工作人员现有的 经验得到策略，输出相应的预警信号。2. 根据权利要求1所述的基于ICA算法的智能变电站设备状态预警方法，其特征在于： 所述步骤(三)包括以下步骤： 1) 数据预处理:控制室的控制中心对接收到的观测数据进行去均值或加上均值向量后 得到正确的源信号波形估计;然后利用白化法去除数据间的相关性； 2) 进入ICA主算法:控制室的控制中心利用傅里叶算法，以负熵最大作为搜寻方向对预 处理后的数据进行快速收敛，并通过非高斯性度量来表示分离结果间的相互独立性，当非 尚斯性达到最大时，完成ICA分尚。3. 根据权利要求2所述的基于ICA算法的智能变电站设备状态预警方法，其特征在于： 所述步骤1)对数据预处理后得到的数据信号集必须满足每个采集点的温度值之间是相互 独立的、得到的分离矩阵为满秩矩阵、得到的数据集中至多只有一个信号为高斯分布函数 二个先验条件。4. 根据权利要求2或3所述的基于ICA算法的智能变电站设备状态预警方法，其特征在 于:所述步骤2)中利用傅里叶算法进行ICA分离包括依次进行的以下步骤： ①随机选择初始权值向量；③ 利用公式Wk+l:Wk+l/| |wk+l| I进行归一化； ④ 如果I Wk+1-Wk I >&，则算法不收敛，返回步骤②，否则利用傅里叶算法的ICA算法估算 出了 一个独立分量，完成了所采集数据集的提取。5. 根据权利要求1至3中任意一项所述的基于ICA算法的智能变电站设备状态预警方 法，其特征在于:所述步骤(四）中利用的贝叶斯算法包括依次进行的以下步骤： i)建立温度模型； ? )预处理模型参数； m)参数带入模型，进行故障预测； iv )在线采集温度，更新温度； V )得到此时刻温度分布； Vi)结合此时刻温度分布、温度总体分布、温度先验分布，依据观测值预测此时刻后验 分布。6. 根据权利要求4所述的基于ICA算法的智能变电站设备状态预警方法，其特征在于： 所述步骤(四）中利用的贝叶斯算法包括依次进行的以下步骤： i)建立温度模型； ? )预处理模型参数； m)参数带入模型，进行故障预测； iv )在线采集温度，更新温度； V )得到此时刻温度分布； Vi)结合此时刻温度分布、温度总体分布、温度先验分布，依据观测值预测此时刻后验 分布。7. 根据权利要求5所述的基于ICA算法的智能变电站设备状态预警方法，其特征在于： 所述步骤m)中利用贝叶斯算法进行预测的过程为:设温度的初始状态的概率分布和k-1时 刻的概率分布分别为:.，Xk只与Xk-i有关，得到：式中不包含yk的先验分布，由系统的状态转移概率分布来计算； 所述步骤vi)更新后验分布参数的过程为:利用最新观测变量yk和先验分布沪 推导#，由贝叶斯公式可得：由于龙滅j，把上述两式代入贝叶斯公式， 可得：再根据条件概率分布、联合概率分布和贝叶斯公式8. 根据权利要求6所述的基于ICA算法的智能变电站设备状态预警方法，其特征在于： 所述步骤m)中利用贝叶斯算法进行预测的过程为:设温度的初始状态的概率分布和k-Ι时 刻的概率分布分别为:，Xk只与Xk-i有关，得到：式中不包含yk的先验分布，由系统的状态转移概率分布·来计算； 所述步骤vi)更新后验分布参数的过程为:利用最新观测变量yk和先验分布^ 推导，由贝叶斯公式可得：，把上述两式代入贝叶斯公式， 可得：再根据条件概率分布、联合概率分布和贝叶斯公式
【文档编号】H02J13/00GK105932774SQ201610309884
【公开日】2016年9月7日
【申请日】2016年5月11日
【发明人】孙涛, 刘耀根, 任启, 袁硕, 穆润根, 刘晓丽, 姚学科, 苏子强, 李睿琦, 李鹏飞, 魏英星
【申请人】国网冀北电力有限公司张家口供电公司, 国家电网公司