发明名称 一种复合绝缘子老化状态分级及判别方法
摘要 本发明公开了属于高电压与绝缘技术领域的一种复合绝缘子老化状态分级及判别方法,该方法采用聚类分析方法,以选取憎水性丧失试验的丧失静态接触角均值X<sub>1</sub>、丧失静态接触角最小值X<sub>2</sub>和憎水性恢复试验的恢复静态接触角均值X<sub>3</sub>以及热刺激电流(TSC)的陷阱电荷量X<sub>4</sub>、陷阱能级X<sub>5</sub>这5个变量对18支试样进行老化状态量化分级;并通过判别分析方法对3支未知老化状态试样进行判别,得到各个老化等级判别函数;其中,判别标准根据陷阱电荷量大小,由聚类分析结果给出,老化状态分级结果通过判别分析方法来验证精确性,本发明可以通过测量陷阱的电荷量来判断出聚合物材料的老化程度,避免事故发生造成各种损失。
申请公布号 CN105403777A 申请公布日期 2016.03.16
申请号 CN201510570741.8 申请日期 2015.09.09
申请人 华北电力大学 发明人 屠幼萍;王璁;李天福;龚博;张辉;佟宇梁
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 张文宝
主权项 一种复合绝缘子老化状态分级及判别方法,其特征在于,采用聚类分析方法,以选取憎水性丧失试验的丧失静态接触角均值X1、丧失静态接触角最小值X2和憎水性恢复试验的恢复静态接触角均值X3以及热刺激电流(TSC)的陷阱电荷量X4、陷阱能级X5这5个变量对18支试样进行老化状态量化分级;并通过判别分析方法对3支未知老化状态试样进行判别,得到各个老化等级判别函数;其中,判别标准根据陷阱电荷量大小,由聚类分析结果给出,老化状态分级结果通过判别分析方法来验证精确性,具体包括:1)对复合绝缘子参数包括静态接触角、陷阱电荷和年限进行检测,得到如表1所示的原始数据,如表1所示的数据具有不同的量纲、而且数值呈数量级之差的数据要在一起进行比较,必须进行数据变换处理;然后采用聚类分析,获取老化状态试样判别标准的基本信息;表1试样基本信息<img file="FDA0000799033960000011.GIF" wi="1712" he="1519" /><img file="FDA0000799033960000021.GIF" wi="1710" he="1295" />其中,1~18号试样用于聚类分析标准,17~19号试样为新的复合绝缘子试样,19~21号试样用于判别分析;i为第几个样品,i=1、2、…n;,n为样品个数=21;2)复合绝子老化状态分级1.1采用聚类分析对样品进行分并类,采用组间连接法将n个样品根据陷阱电荷量大小分组,一个组;先计算n(n-1)/2个相似性测度即相似度大小,并且把具有最小测度的两个样品合并成成为一类,然后按照系统聚类分析方法计算复合绝缘子每一种老化状态等级和其它n-2个样品之间的距离进行并类;在并类过程中,每一步所做的并类(样品与样品、样品与类、类与类)都要使相似性测度在系统中保持最小,这样每次减少一类,直至所有样品都归类为止;一个组包含的样品个数,计算该组的陷阱电荷量平均值,一个平均值为一类,对表1中1~18号复合绝缘子样品聚类分析结果如表2所示;表1所示,1‑18个样品的聚类分析结果为4个类,表2复合绝缘子聚类分析结果<img file="FDA0000799033960000022.GIF" wi="1989" he="145" /><img file="FDA0000799033960000031.GIF" wi="1990" he="719" />1.2数据变换处理采用标准化变换对变量的属性进行变换处理,首先对列进行中心化,然后用标准差给予标准化,即<img file="FDA0000799033960000032.GIF" wi="399" he="134" />其中<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799033960000033.GIF" wi="967" he="209" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>S</mi><mi>j</mi></msub><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000799033960000034.GIF" wi="997" he="284" /></maths>上述式(1)、式(2)中,i代表第i个样本,j代表样本第j个变量;通过变换处理后,每列数据的平均值为0,方差为1,使用标准差作标准化处理后,计算距离系数,在抽样样本改变时仍保持相对稳定性;1.3计算距离系数,包括:欧式距离:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>X</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>X</mi><mrow><mi>k</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FDA0000799033960000035.GIF" wi="562" he="154" /></maths>Minkowski距离:<img file="FDA0000799033960000036.GIF" wi="549" he="135" />其中,d<sub>ik</sub>为第i个样本和第k个样本的距离大小;根据表1所示数据计算,得到表3所示的聚类过程中不同老化等级之间平均中心(指什么地方的:就是指类)之间距离;表3聚类过程中不同类之间距离<img file="FDA0000799033960000037.GIF" wi="1749" he="144" /><img file="FDA0000799033960000041.GIF" wi="1750" he="1350" />3)老化状态分级,按照步骤1.1分组中的一个组为一类,表1中1‑18个样品的分为4个类,与表2对照,复合绝缘子老化状态分级结果如下面表4.表4复合绝缘子老化状态分级结果<img file="FDA0000799033960000042.GIF" wi="1797" he="631" />4)复合绝缘子老化状态分级结果验证(判别分析)4.1判别分析判别分析法通常从各训练样本中提取已有的各总体信息,构造一定的距离判别准则,来判断新样本属于那个总体,其中距离判别的基本思想是:比较样本和每个总体的马氏距离,并将其判定属于马氏距离最近的那个总体,最后通过Fish准则得到Fish判别函数;聚类分析将表1中1‑18支复合绝缘子试样老化状态分类(分成4),从分类结果来看,微观特性与宏观特性基本相符合,但是需要通过另一些试样进行判别验证;所述马氏距离为设总体G={X<sub>1</sub>,X<sub>2</sub>,X<sub>3</sub>,X<sub>4</sub>,X<sub>5</sub>}<sup>T</sup>为m元总体(这里考虑m个指标,m=5),令μ<sub>i</sub>=E(X<sub>i</sub>)(i=1,2,···m),则总体均值向量为μ={μ<sub>1</sub>,μ<sub>2</sub>,···,μ<sub>m</sub>}<sup>T</sup>,总体G的协方差矩阵为:Σ=cov(G)=E[(G‑μ)(G‑μ)<sup>T</sup>]   (5)定义样本x与总体G的马氏距离为:d<sup>2</sup>(x,G)=(x‑μ)<sup>T</sup>Σ<sup>‑1</sup>(x‑μ)   (6)设有4个m元总体:G<sub>1</sub>,G<sub>2</sub>,G<sub>3</sub>,G<sub>4</sub>,其各自的均值向量和协差矩阵分别为μ<sub>i</sub>,Σ<sub>i</sub>(i=1,2,3,4),对任意给定的m元样本x={x<sub>1</sub>,x<sub>2</sub>,···,xm},判断其属于哪个总体,可以按马氏距离最近准则对x进行判别分类,首先计算样本x到k个总体的马氏距离d<sup>2</sup>(x,G<sub>i</sub>)(i=1,2,3,4),然后进行比较,把样本x判归距离最小的那个总体,设i=l时,若d<sup>2</sup>(x,G<sub>l</sub>)=min<sub>i=1,2,3,4</sub>{d<sup>2</sup>(x,G<sub>t</sub>)}   (7)则样本x∈G<sub>l</sub>;其中,X<sub>1</sub>,X<sub>2</sub>,X<sub>3</sub>,X<sub>4</sub>,X<sub>5</sub>与前面“丧失静态接触角均值X1、丧失静态接触角最小值X2和憎水性恢复试验的恢复静态接触角均值X3以及热刺激电流(TSC)的陷阱电荷量X4、陷阱能级X5这5个变量的含义是一致的;4.2Fisher判别及结果Fisher判别是一种依据方差分析原理建立的判别方法,对马氏距离判别的老化状态分级结果进行重新训练和验证,对p维空间中的点x<sub>i</sub>=(x<sub>i1</sub>,x<sub>i2</sub>,···,x<sub>ip</sub>),i=1,2,···,n,找一组线性函数y<sub>m</sub>(x<sub>i</sub>)=Σc<sub>f</sub>,m=1,2,···,m,一般m<p,用它们把p维空间中的观测点都转换为m维的,再在m维空间中对观测集进行分类,降维后的数据应最大限度地缩小同类中观测之间的差异,并最大限度扩大不同类别间的差异,如此才能获得较高的判别效率,在此采用方差分析的思想,依据使组间均方差与组内均方差之比最大的原则,选择最优的线性函数,以对新抽检的试样的信息进行判别,然后对抽检试样进行老化状态等级归类。
地址 102206 北京市昌平区朱辛庄北农路2号