本发明涉及交直流大电网广域状态监测技术领域，具体涉及到一种基于中心距和样本特征的换流变压器状态评估方法。
背景技术：
换流变压器是直流输电系统中向换流器提供适当电压等级的三相交流电压源设备，具有参与交直流变换、抑制直流故障电流、削弱交流系统入侵直流系统的过电压、隔离交直流系统、减少交流系统谐波的作用，其安全可靠运行直接关系到直流输电系统的运行可靠性。
随着直流输电工程运行年限的增长，已投运的换流变压器的绝缘问题也日益突出。葛上工程中第一批次换流变压器主绝缘试验中多次出现阀侧绕组、引线、分接开关等多个部位的绝缘问题，在西门子、abb各自生产的首台±800kv换流变试验过程中均出现了绝缘故障。因此对运行中期的换流变压器进行状态评估，并根据评估结果制定相应的检修方案，对提高运行中期的换流变压器的安全可靠性和保障直流输电工程的安全稳定运行具有十分重要的意义。
目前基于状态检修的要求，国内外已经涌现了大量的电力变压器状态评估的方法，主要分为两大类：一类是基于现场运维经验和专家知识的主观评价方法，这种方法的优点是实用性强，但是评价结果精度容易受到评估人员主观因素的影响；另一类是基于运维数据和人工智能算法的客观评价方法，这种方法具有全面客观的优点，但是像神经网络、贝叶斯网络等算法需要大量的故障样本数据来驱动，换流变压器因投运的时间较短可用于驱动人工智能算法的数据样本规模小，不能满足评价精度要求。另外，换流变压器与普通电力变压器在短路阻抗、阀侧绕组绝缘、阀侧套管积污、内部电磁环境等方面存在较大差异，且在运行负荷和外部环境的影响下，换流变压器的运行特征呈现出较强的时间差异性，为了对其运行状态进行精确评估，需要建立能体现时间差异性的换流变压器状态评估方法。
虽然对电力变压器运行状态评价的方法较多，但目前还没有针对换流变压器的运行状态评价方法，所以提出一种基于换流变压器运行特征的状态评价方法具有重要的工程实用价值。
技术实现要素：
本发明的目的在于克服上述现有技术的不足，提供一种基于中心距和样本特征的换流变压器状态评估方法，以实现对换流变压器状态评进行精确评估。
为实现上述目的，本发明的技术方案是：
一种基于中心距和样本特征的换流变压器状态评估方法，所述方法包括如下步骤：
步骤一：根据换流变压器的结构及运行特点，划分换流变压器评估层次，确定故障层的故障类型及指标层的状态量；
步骤二：统计换流变压器在检修周期内不同时间段p个的运维数据样本cn，同时采集待评价换流变压器指标层状态量数据，并进行归一化处理，得到换流变压器实测指标量的状态量的归一化评分gi；
步骤三：根据其中的一运维数据样本来确定其故障类型的相关状态量，列写模糊判断矩阵，求得指标量权重ωi和故障类型评分xn；
步骤四：利用各故障类型评分，确定故障类型模糊判断矩阵，求得故障类型权重向量ωi′，以及基于该样本的换流变状态评估值s1；
步骤五：按步骤三、四的流程依次对p个运维数据样本进行换流变状态评估，得到基于样本的换流变状态值
步骤六：依据故障类型权重集合求取样本中心距dl及样本融合权重ωl，以求得换流变最终的状态评估值w，根据w的置信区间确定检修方案。
本发明与现有技术相比，其有益效果在于：
本发明提供了一种适用于大量样本数据特征的换流变压器状态评估方法。本发明主要包括样本特征提取和结果评价两个阶段。样本特征提取阶段包括了换流变压器特征的状态评估指标体系的建立、特征指标量及故障类型权重计算、样本特征矩阵的求取。结果评价阶段给出基于不同样本不同特征的换流变压器状态评分，依据中心距法给出换流变压器状态的最终评估结果。
本方法具有如下技术优势：
(1)首次将模糊评法运用于换流变压器状态评价中，实现了对换流变压器健康状态的检测，为换流变压器的运行与维护提供有力依据。
(2)首次建立较全面的换流变压器状态评估指标体系，综合考察24项换流变压器运维过程中指标量和8种故障类型，基本覆盖了换流变压器常见的运维故障和缺陷，充分体现了换流变结构特点和运行特性。
(3)本方法以实测数据和运维统计数据为基础，确定换流变压器模糊层次分析过程中的权重值，避免了模糊分析过程中主观因素的影响，使评价结果更具客观性。
(4)本方法以大量换流变运维数据样本为基础，建立了多个带有样本特征的换流变压器模糊状态评估模型，通过中心距融合法，提高了换流变压器综合模糊状态评估模型的精确度。
附图说明
图1为换流变压器阀侧绕组工作原理图；
图2为本发明的系统图；
图3为本发明实施例的样本训练流程图。
具体实施方式
下面结合附图和具体实施方式对本发明的内容做进一步详细说明。
实施例：
一种基于中心距和样本特征的换流变压器状态评估方法，具体步骤为：
步骤一：根据换流变压器的结构及运行特点，划分换流变压器评估层次，确定故障层的故障类型及指标层的状态量；
步骤二：统计换流变压器在检修周期内不同时间段p个的运维数据样本cn，同时采集待评价换流变压器指标层状态量数据，并进行归一化处理，得到换流变压器实测指标量的状态量的归一化评分gi；
步骤三：根据其中的一运维数据样本来确定其故障类型的相关状态量，列写模糊判断矩阵，求得指标量权重ωi和故障类型评分xn；
步骤四：利用各故障类型评分，确定故障类型模糊判断矩阵，求得故障类型权重向量ωi′，以及基于该样本的换流变状态评估值s1；
步骤五：按步骤三、四的流程依次对p个样本进行换流变状态评估，得到基于样本的换流变状态值
步骤六：依据故障类型权重集合求取样本中心距dl及样本融合权重ωl，以求得换流变最终的状态评估值w。
其中，在步骤一中，换流变压器的结构特点如下：
1)侧绕组及其引线的绝缘复杂。由图1来描述其阀侧绕组工作原理图，相较于普通电力变压器，其阀侧绕组在一个工频周期内要承受交直流瞬变冲击电压和极性反转电压的应力作用，而且在直流电场下，阀侧绕组的绝缘纸将承担大部分的绝缘压力，这使换流变阀侧绕组及其引线的绝缘复杂多变，成为运维过程中绝缘监测的重点。
2)阀侧套管工作于交直流电场。阀侧套管连接阀侧绕组与换流阀，在一个工频周期内承受交流电场与直流电场共同作用。直流电场作用下，套管表面更易吸引空气中的带电颗粒形成污秽，从而加大发生污闪及不均匀受潮放电的几率，套管的运行状态受空气污染的影响而呈现季节性差异。
3)内部电磁环境更恶劣。为降低阀侧绕组端部对地高电压的接地风险，在设计时会增加绕组端部与铁芯轭的距离，这不可避免地导致换流变压器内部漏磁和局部损耗的增加。随着负荷的增加，换流变压器内部电磁环境更加恶劣。
步骤一中，换流变压器的运行特点为：故障频发的主要部位的次序大致为绕组、铁芯、分接开关、套管、绝缘油。
步骤一中，换流变压器评估层次分为系统层、故障层和指标层，故障层中的故障类型有绕组故障、铁心故障、内部过热、绝缘受潮、固体绝缘劣化、电弧放电、局部放电、阀侧套管故障；指标层状态量如表1所示。
表1指标层状态量
步骤二中，一个检修周期为一年，不同时间段可按季节划分，即每个季节为一个时间段，有4个不同时间段的运维数据样本。
步骤二中的运维数据包括故障层中故障发生的次数，以及某种类型的故障发生时相关指标状态量超标的次数。
步骤二中归一化处理方法为；定义实际值在正常范围内相对于正常状态阈值的裕度为相对完好度，则对于测量值越大状态越优的指标，如极化指数、绝缘电阻，其相对完好度可按公式(1)确定：
对测量值越小状态越优的指标，如乙炔含量、局部放电量，其相对完好度按公式(2)确定：
其中，g(xi)为指标xi的相对完好度；xi为该指标测量值，i表示指标层编号；xip表示该指标量的初值；xi0为该指标量的警示值，其值均按《油浸式变压器(电抗器)状态评价导则》要求确定。各指标量的警示值如表2所示。
步骤三中某一故障类型的相关状态量为该故障发生情况下，测量值出现过异常的所有指标量。
步骤三中对某故障类型状态评价时需要建立其相关状态量模糊判断矩阵ra，模糊判断矩阵ra可按公式(3)确定：
其中矩阵元素rij可按公式(4)确定：
公式(4)中pi(i＝1,2,…,n)表示在同一统计样本下某故障类型的相关状态量的频数。
步骤三中判断矩阵的一致性检验可按如下方法进行验证：
模糊层次分析法判断矩阵r＝(rij)n×n的一致性应满足：
rij＝rik﹣rik+0.5公式(5)
或满足顺序一致性：对模糊判断矩阵r，
若
所述步骤三中指标量权重可按公式(6)确定：
公式(6)中β为权重分辨系数，通过增大β值可以提高各元素权重的分辨率，为充分体现各相关状态量之间的差异性本发明中β取ee(e为自然对数)。
公式(6)中rij为公式(3)中矩阵元素。
步骤三中故障类型状态评分可按公式(7)确定：
步骤四中故障类型模糊判断矩阵按公式(8)确定：
其中，矩阵元素rij可按公式(9)确定：
公式(9)中xi为按公式(7)确定的故障类型状态评分。
步骤四中故障类型的权重ωi′可按照公式(10)确定：
公式(10)中权重分辨系数β取ee(e为自然对数)。
步骤四中，换流变压器状态值s按公式(11)确定：
公式(11)中s为换流变压器状态值，ω′i为故障类型权重，xi为故障类型评分。
步骤五中，对p个样本进行换流变压器状态按照步骤三和步骤四所述内容进行评估，p个换流变压器样本状态值sp，组成状态向量s＝[s1，s2，…，sp]。
步骤五中样本特征矩阵为：
其中vp为样本p的特征向量，vp＝[v1pv2p…vnp]，vnp为对第p号样本进行评估时第n类故障的权重。
步骤六中，样本特征向量中心距dj，可按照公式(12)确定：
公式(12)中vij表示故障类型fi的权重值，j为样本编号，j＝1,2，…，p；i为故障类型编号，i＝1,2，...，n。
公式(12)中
步骤六中，样本融合权重wj具有如下性质：
1.样本融合权重表征某样本评价结果与所有样本初评结果的一致性程度；
2.样本融合权重与中心距负相关，即融合权重越大中心距越小，融合权重越小中心距越大。
步骤六中，样本融合权重wj可以按照公式(14)确定：
公式(14)中dj为样本个体的特征向量中心距。
公式(14)中a为与样本数p有关的常数，可以按照公式(15)确定。
公式(15)中k为差异系数，k值越大对样本特征的区分度就越大，一般取100。用matlab对公式(15)求解，得到a与样本数p之间的对应关系，如表3所示。
表3常数a与样本个数p的关系
步骤六中，换流变最终的状态评估值w，可以按照公式(16)确定：
w＝w·s公式(16)
公式(16)中w为样本融合权重w＝[w1，w2，···，wp]；
公式(16)中s为依据样本的换流变压器模糊状态评估结果s＝[s1，s2，···，sp]。
由此可知，本方法提供了一种适用于大量样本数据特征的换流变压器状态评估方法。本发明主要包括样本特征提取和结果评价两个阶段。样本特征提取阶段包括了换流变压器特征的状态评估指标体系的建立、特征指标量及故障类型权重计算、样本特征矩阵的求取。结果评价阶段给出基于不同样本不同特征的换流变压器状态评分，依据中心距法给出换流变压器状态的最终评估结果。
下面以收集整理到的中国南方电网公司所辖换流变压器故障和缺陷信息为统计数样本，以中国南方电网公司某换流站极1c相换流变压器为评估对象，对该换流变压器运行状态进行评估，本实例整体框架如图2所示，主要包括样本特征提取(即步骤一至步骤四)和中心距融合评价(即步骤五、步骤六)两个阶段。用于提取样本特征的数据为收集整理到的中国南方电网公司换流变压器故障和缺陷统计信息，如表4所示。待测评换流变压器的24个指标量实测值、规定注意值及状态评分如表5所示。样本特征提取阶段包括了换流变压器特征的状态评估指标体系的建立、特征指标量及故障类型权重计算、个体样本特征矩阵的求取。结果评价阶段给出基于不同样本不同特征的换流变压器状态评分，依据中心距法给出换流变压器状态的最终评估结果。
表4换流变压器故障信息样本统计数据
表5换流变压器指标量实测数据及其状态评价值
作为本实施例的一种优选方案，上述的样本特征提取部分的实施步骤如下：
(1)统计换流变压器发生故障或存在缺陷时所选的24个指标量状态值异常的次数，如表4所示。
(2)以表中数据为基础，根据公式(3)和公式(4)求出状态量评价矩阵r＝(rij)n×n，并根据公式(5)或顺序一致性原则检验状态量评价矩阵r＝(rij)n×n的一致性。
(3)在状态量评价矩阵r＝(rij)n×n满足一致性检测的前提下，根据公式(6)计算换流变压器指标量权重ωi，根据公式(7)计算各故障类型状态评分xn。
(4)根据公式(8)和公式(9)确定故障类型模糊判断矩阵rf，根据公式(5)或顺序一致性原则检验验证其一致性。
(5)在故障类型模糊判断矩阵rf满足一致性检测的前提下，根据公式(10)确定故障类型的权重ω′i，按公式(11)确定换流变压器模糊状态评分s。
作为本实施例的的一种优选方案，上述的中心距融合法部分的实施步骤如下：
(1)根据样本特征提取部分故障类型的权重ω′i和故障类型评分s，得到样本p的特征向量vp＝[v1pv2p…vnp]。
(2)求样本特征向量中心距dj，先根据公式(13)求出样本中故障类型fi权重的平均值
(3)求样本融合权重wj，根据公式(14)确定样本融合权重wj，根据公式(16)确定换流变最终的状态评估值w。
作为本发明的一种优选方案，换流变压器状态评估值w的置信区间如下：
当换流变压器刚投运时，其状态评估值为1；当发生重大故障而被迫退出运行时，状态评估值为0。
当w∈(0.8,1]时，将其状态等级评定为良好，表明换流变压器运行状态良好能继续稳定运行可适当推迟检修计划。
当w∈(0.6,0.8]时，将其状态等级评定为注意，表明换流变压器运行状态一般，能继续运行，但其内部可能存在微小缺陷，应当按计划检修。
当w∈(0.2,0.6]时，将其状态等级评定为异常，表明换流变压器运行状态较差，能勉强运行，但其内部可能存在较大缺陷，应当尽快安排检修。
当w∈(0,0.2]时，将其状态等级评定为严重，表明换流变压器运行状态很差，不能继续运行，但其内部可能存在影响电网运行的缺陷，应当立即退出运行安排检修。
本实施例最终的评定结果，w＝0.698其状态等级评定为注意，应当对其进行按期检修。但该换流变压器近期经历过短时间短路电流冲击，在后续运行过程中该换流变压器运行状态持续下降，电网公司决定对该换流变压器进行解体大修。解体大修发现了该换流变压器的确存在绕组缺陷，证明了本评估方法的可行性与准确性。
上述实施例只是为了说明本发明的技术构思及特点，其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所做出的等效的变化或修饰，都应涵盖在本发明的保护范围内。