Determination of Optical Constants by Double Thickness Transmittance Model Based on Polynomial Root
-
摘要: 为解决光谱反演法确定物质光学常数的一些问题,基于传统的双厚度透射率模型,建立厚度分别为L和2L的光谱透射率方程,通过代数运算获得与衰减系数有关的八次多项式方程,求解并选择其大于0小于1的实数根来计算衰减系数和消光系数;再求解关于界面反射率的一元二次方程,选择其大于0小于1的根来计算折射率。在确定光学常数的过程中,新方法没有反演误差和迭代计算耗时问题。利用已知文献中庚烷的光学常数验证新方法的可靠性,并分析了双厚度不满足2倍关系时对计算结果的影响,结论是第二厚度2L的相对误差不超过1%时,消光系数的计算误差不超过2.03%,不考虑3个强吸收点时,折射率的计算误差不超过1%。Abstract: To solve some problems in the determination of substance optical constants by spectral inversion, based on the traditional double thickness transmittance model, spectral transmittance equations with thicknesses of L and 2L are established. The eighth-order polynomial equation related to the attenuation coefficient is obtained through an algebraic operation. Real roots greater than 0 and less than 1 are solved to calculate the attenuation coefficient and extinction coefficient. Then, the quadratic equation with one unknown quantity about the interface reflectivity is solved, and the roots greater than 0 and less than 1 are selected to calculate the refractive index. In the process of determining the optical constants, the new method has no inversion error or iterative calculation time-consuming problems. The reliability of the new method is verified using the known optical constants of heptane from the literature, and the influence on the calculation results was analyzed when the double thickness does not satisfy the double relationship. In conclusion, when the relative error of the second thickness 2L was no more than 1%, the calculation error of the extinction coefficient was no more than 2.03%, and when three strong absorption points were not considered, the calculation error of the refractive index was no more than 1%.
-
0. 引言
线路绝缘子主要承担着连接导体和电气绝缘的功能,是整个电力系统的重要的组成部分[1-3]。由于绝缘子在生产过程中会造成一定的缺陷,以及受到自然环境的影响,绝缘子会自然劣化,绝缘性能不断减弱,绝缘子串的闪络概率将增大,最终造成电网运行的不稳定,而电力事故发生的概率也将增大,给生产生活带来不利的影响。因此绝缘子的定期检测与及时维修对于维护保障电网的安全至关重要[4-5]。
目前主流的低零值(根据DL/T596-1996《电力设备预防性试验规程》中的要求,每片悬式绝缘子的绝缘电阻不应低于300 MΩ,500 kV悬式绝缘子不低于500 MΩ。低于上述水平的,一般就认为是低值或零值绝缘子。)绝缘子检测方法有:光谱法、紫外脉冲法、径向温度法、超声波检测法等。但分析文献,发现不少方法存在危险性高、算法复杂等问题,均需进一步深入研究。由于低零值绝缘子在线路中温度变化明显,目前许多电力公司逐步采用红外成像技术对低零值绝缘子进行检测。然而目前图像处理技术对电力设备进行在线检测的研究面临着图像特征提取的困难,现有算法无法有效解决绝缘子状态检测的多分类问题,且面临处理海量数据检测耗时,检测正确率低的问题。因此,面向大数据的低零值绝缘子检测方法是今后研究的重点。
随着神经网络算法的不断改进,以深度学习为代表的人工智能理论与应用研究越来越多的被应用到故障检测识别中。目前BP(back propagation)神经网络、遗传算法[6]、Petri网络及决策树等不少数据挖掘的方法被成功应用到劣化绝缘子的诊断识别中。支持向量机同样被应用在电力系统的故障诊断领域,然而,直接采用支持向量机模型对绝缘子样本进行检测效果不尽理想。如何优化原有的支持向量机模型,解决大数据环境下绝缘子检测问题,是当下需要着重研究的方向。目前网格搜索法[7]、布谷鸟搜索算法[8]、粒子群算法[9]等等都成功被应用到优化向量机的参数寻优中。灰狼算法与支持向量机相结合应用在诸多领域中,如医学信号识别、植物种类识别、医学图像识别,其实验结果都有所改善,但很少被应用到电气设备故障诊断。本文提出的灰狼优化算法与支持向量机结合针对低零值绝缘子检测识别的应用尚属空白。
本文通过对绝缘子红外图像进行处理,对绝缘子红外图像样本进行多层次深度特征提取用于支持向量机分类识别,并采用灰狼算法实现对支持向量机参数的优化,实现对低零值绝缘子检测识别。
1. 绝缘子串红外图像预处理
1.1 信号绝缘子串红外图像增强
灰度变换增强可以增强红外图像中的目标与背景的对比度,提高图像的质量。灰度变换作为一种应用广泛的图像增强技术可使图像清晰、特征明显[10-12]。
设原图像为f(x, y),其灰度范围为[a, b];变换后的图像g(x, y),其灰度范围线性的扩展至[c, d]。
$$g\left( {x, y} \right) = \frac{{d - c}}{{b - a}}\left[ {f\left( {x, y} \right) - a} \right] + c$$ (1) 图像中大部分灰度级分布在区间[a, b]内,有少许部分在此区间之外,为了改善增强效果,可以令:
$$g\left( {x, y} \right) = \left\{ \begin{array}{l} c \\ \frac{{d - c}}{{b - a}}\left[ {f\left( {x, y} \right) - a} \right] + c \\ d \\ \end{array} \right.$$ (2) 直方图均衡化作为一种应用广泛的图像增强方法,可使绝缘子串同背景对比度增大,方便后期提取绝缘子串。
如图 1所示,直方图均衡化增大了绝缘子串与背景的灰度级,图像的对比度也增强了,这样有利于后期正确分割出绝缘子串和背景。
1.2 基于Ostu算法的图像自适应阈值分割
首先利用最大类间方差法(Ostu)对增强后的绝缘子红外图像进行分割[13],如图 2所示。该方法可自动选取阈值,分割效果好、速度快。
1.3 绝缘子分割
对Ostu分割得到的二值图像进行切割,提取图像中完整的绝缘子串,如图 3所示,为绝缘子缺陷检测智能认知工作做好准备。
1.4 基于Randon变换的图像倾斜校正
运用Ostu算法将绝缘子串的候选区域分割提取后,用Randon变换[14]的图像倾斜校正算法进行绝缘子角度校正,如图 4所示。
2. 灰狼算法支持向量机模型
2.1 支持向量机模型
支持向量机(support vector machine, SVM)是Bell实验室以V. Vapnik教授为首的研究小组针对小样本机器学习方法提出的一种新型模式识别方法[15]。
对于线性不可分的特征向量,需采用核函数将向量投放到高维空间中达到可以分类的效果。高斯径向基核函数作为应用最广泛的核函数,在缺乏样本数据的先验知识时,可通过调整参数取得较好的学习效果。本文采用高斯径向基核函数:
$$K\left( {{x_i}, {x_j}} \right) = {\rm{exp}}\left( { - \frac{{\left\| {{x_i} - {x_j}} \right\|}}{{2{\delta ^2}}}} \right)$$ (3) 式中:δ>0为高斯核的带宽。
设h维的空间上,针对线性可分问题,所有样本均满足约束件:
$${y_i}\left( {{\mathit{\boldsymbol{\omega }}^{\rm{T}}}{X_i} + {\omega _0}} \right) - 1≥0, i = 1, 2, {\rm{L}} $$ (4) 求解支持向量机可转化为分类间隔问题:
$$\mathop {{\rm{min}}}\limits_{\omega , b} \frac{1}{2}{\left\| \mathit{\boldsymbol{\omega }} \right\|^2}$$ (5) 在实际机器学习时,为了允许机器出现一些错分的点,通常在约束条件中加入松弛变量ζ>0,增加一个常数C作为惩罚因子:
$$\mathop {{\rm{min}}}\limits_{\omega , b, \zeta } \left( {\frac{1}{2}{{\left\| \mathit{\boldsymbol{\omega }} \right\|}^2} + C\sum\limits_{i = 1}^h {{\zeta _i}} } \right)$$ (6) 以上问题的求解,可以得到SVM回归表达式:
$$y\left( x \right) = {\rm{sgn}}\left( {\sum\limits_{j = 1}^h {{\alpha _j}{y_j}K\left( {{x_j}, x} \right) - b} } \right)$$ (7) 式中:αj为拉格朗日乘积因子。
2.2 灰狼优化算法
灰狼优化算法(grey wolf optimizer,GWO),是一种通过模拟灰狼捕猎过程中的狩猎和搜索行为建立的全局随机搜索算法。由澳大利亚学者Seyedali. Mirjalili等人在2014年提出的新型算法[16]。GWO算法与粒子群优化算法(particle swarm optimization,PSO)类似,是一个从随机解出最优解的过程。该方法相较于PSO、网格搜索算法(GS)等算法参数少,结构简单,同时又有较强的收敛性,已成功应用于图像处理等领域中。
将最优解设为α,第二个和第三个最佳解分别命名为β和δ,而其余的解均设为ω。狼群通过3只个体狼α、β和δ为初始解带领狼群ω在空间中向猎物(最优解)逼近,经过图 5所示的狼群移动方式,不断迭代,引导狼群不断靠近全局最优解。搜索过程狼群捕食位置更新:
$$D = \left| {C \cdot {X_{\rm{p}}}\left( t \right) - X\left( t \right)} \right|$$ (8) $$X\left( {t + 1} \right) = {X_{\rm{p}}}\left( t \right) - A \cdot D$$ (9) 式中:D为当前灰狼距猎物距离;A和C为系数向量;Xp是猎物的位置向量。
$$A = 2 \cdot a \cdot {r_1} - a$$ (10) $$C = 2 \cdot {r_2}$$ (11) 式中:a随迭代次数从2~0递减;r1,r2是[0, 1]内的随机向量。
为了模拟狩猎行为,假设α,β和δ对猎物的潜在位置有更好的了解,在每次迭代过程中,保留当前最优的α,β和δ解。
$$\left\{ \begin{array}{l} {D_\alpha } = \left| {{C_1} \cdot {X_\alpha }\left( t \right) - X\left( t \right)} \right| \\ {D_\beta } = \left| {{C_2} \cdot {X_\beta }\left( t \right) - X\left( t \right)} \right| \\ {D_\delta } = \left| {{C_3} \cdot {X_\delta }\left( t \right) - X\left( t \right)} \right| \\ \end{array} \right.$$ (12) $$\left\{ \begin{array}{l} {X_1} = {X_\alpha } - {A_1} \cdot {D_\alpha } \\ {X_2} = {X_\beta } - {A_2} \cdot {D_\beta } \\ {X_3} = {X_\delta } - {A_3} \cdot {D_\delta } \\ \end{array} \right.$$ (13) $${X_{\rm{p}}}\left( {t + 1} \right) = \frac{{{X_1} + {X_2} + {X_3}}}{3}$$ (14) 式中:Xα,Xβ和Xδ分别代表α狼,β狼和δ狼当前位置;Dα,Dβ和Dδ分别代表当前狼位置和3只头狼的位置间的距离;A1,A2和A3为随机系数向量;t表示迭代次数。
2.3 GWO-SVM
采用GWO优化算法对绝缘子红外图谱识别的SVM网络核参数惩罚因子C与核宽度δ进行参数优化,以达到图谱分类识别的准确性和泛化能力。
① 输入绝缘子图谱的特征量,选取部分作为SVM的训练集,并将剩余的特征向量集作为测试集,以验证SVM识别的准确率。
② 初始化狼群数量、迭代次数,设置惩罚因子C与核宽度δ的范围。
③ SVM根据初始参数C与δ进行训练和测试,并以错误率最小化为目标。
④ GWO以C与δ为猎物进行优化,达到最大迭代次数时输出GWO全局最优值。
⑤ 将处理后的绝缘子图谱样本分别作为SVM的训练集与测试集。采用最佳参数C与δ建立识别模型,并对测试样本进行预测、分析。
3. 实验结果与分析
3.1 数据描述
为了验证所提出的基于GWO-SVM劣化绝缘子状态检测的可行性,由于目前未建立绝缘子红外图像数据库,我们选取200幅绝缘子图像作为绝缘子样本库。采取随机抽样的方法选取两类样本,其中训练样本120幅,测试样本80幅。所有样本由多位人工分拣专家投票分为完好和低零值两类(如图 6所示)。
3.2 绝缘子检测结果
本文利用灰狼优化算法、粒子群优化算法(random-search)和网格搜索算法(grid-search)对支持向量机参数进行优化。我们对比了3种算法的寻优时间、寻参效率和训练准确率。如表 1所示,灰狼优化算法的各项性能都要比另外两种算法好,其准确率及寻参效率都高于其余两种优化算法。
表 1 参数寻优方法对比Table 1. Comparison of parameter optimization methodsParameter optimization method Accuracy/% Optimization time/s Seeking efficiency/(s/time) Grid-search 91.523 12.693 0.2487 Random-search 92.267 8.159 0.3156 Grey wolf optimizer 95.246 6.251 0.1145 从3种算法优化支持向量机的结果看,网格搜索耗时长且识别准确率低,且寻优时存在复杂度高,运算量大等不足。粒子群优化算法收敛速度快,算法简单,但也存在很明显的缺点,它对于有多个局部极值点的函数,容易陷入到局部极值点中,得不到正确的结果,因此其优化向量机识别的正确率不高。而灰狼优化算法识别准确率可达到95.246%,寻优时间最少且寻参效率高。灰狼算法充分利用先验知识,避免由于惩罚参数过大而导致算法陷入局部最优的风险。因此灰狼优化相比于粒子群搜索算法和网格搜索算法能高效的对低零值绝缘子进行识别。图 7为SVM参数寻优过程。
本文通过网格优化、粒子群优化和灰狼优化这3种优化算法对支持向量机的参数C和δ进行优化。图 8为不同优化算法的故障识别对比。
通过图 8不同优化算法的识别对比,可以看出GWO-SVM诊断方式相比于GS-SVM和PSO-SVM识别正确率更高。对于图 8分类的结果,结合表 1,GWO-SVM错误识别的绝缘子仅有一个,且GWO寻优时间及寻参效率明显优于GS与PSO,满足预设要求。整个绝缘子串检测系统可以实现有效地对低零值绝缘子进行故障诊断,具有工程实际意义。
4. 结束语
1)本文绝缘子红外图像样本进行了多层次深度特征提取,相比于现有的深度模型提取的特征具有更强的鉴别能力。
2)本文提出红外图像和灰狼算法优化支持向量机相结合的方法实现对低零值绝缘子的检测识别,能够在大数据层面准确地识别低零值绝缘子,减少人力,物力以及财力。
3)本文采用灰狼算法优化支持向量机参数,并采用高斯径向基核函数,得到的识别模型识别效果好。
-
-
[1] El-Zaiat S Y, El-Den M B, El-Kameesy S U, et al. Spectral dispersion of linear optical properties for Sm2O3 doped B2O3-PbO-Al2O3 glasses[J]. Opt. Laser Technol. , 2012, 44: 1270-1276. DOI: 10.1016/j.optlastec.2011.12.051
[2] Bridou F, Cuniot-Ponsard M, Jean-Michel D. Experimental determination of optical constants in the vacuum ultra violet wavelength region between 80 and 140 nm: a reflectance versus thickness method and its application to ZnSe [J]. Opt. Commun. , 2007, 271: 353-360. DOI: 10.1016/j.optcom.2006.10.042
[3] Dhanasekaran V, Mahalingam T, Rhee J K, et al. Structural and optical properties of electrosynthesized ZnSe thin films[J]. Optik, 2013, 124(3): 255-260. DOI: 10.1016/j.ijleo.2011.11.063
[4] 王贵全, 张锦荣, 邵毅, 等. 基于透射光谱的类金刚石膜光学参数反演[J]. 红外技术, 2021, 43(5): 473-477. http://hwjs.nvir.cn/article/id/916a8eda-7b0a-4936-a8e2-e659b75abb96 WANG G, ZHANG J, SHAO Y, et al. Calculation of Optical Parameters of Diamond-like Carbon Film Based on Transmission Spectrum [J]. Infrared Technology, 2021, 43 (5): 473-477. http://hwjs.nvir.cn/article/id/916a8eda-7b0a-4936-a8e2-e659b75abb96
[5] 侯典心, 路远, 宋福印. 基于全光谱拟合法的VO2薄膜光学常数计算[J]. 红外技术, 2017, 39(3): 243-249. http://hwjs.nvir.cn/article/id/hwjs201703008 HOU D, LU Y, SONG F. Optical constants of VO2 thin films based on whole optical spectrum fitting. Infrared Technology, 2017, 39 (3): 243-249. http://hwjs.nvir.cn/article/id/hwjs201703008
[6] Tuntomo A, Tien C L, Park S H. Optical constants of liquid hydrocarbon fuels[J]. Combustion Science and Technology, 1992, 84: 1-6, 133-140. DOI: 10.1080/00102209208951841
[7] 李全葆, 宋炳文, 魏天衢. Hg(1-x)CdxTe光学常数测量[J]. 红外技术, 1991, 13 (5): 17-20. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS199105004.htm LI Q, SONG B, WEI T. The measurement of optical constants of Hg(1-x)CdxTe [J]. Infrared Technology, 1991, 13 (5): 17-20. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS199105004.htm
[8] Bohren C F, Huffman D R. Absorption and Scattering of Light by Small Particles[M]. New York: John Wiley & Sons, 1983: 36-41.
[9] 苏星, 李正芬, 刘成赞, 等. 一种红外硒化物玻璃的光学常数及其增透膜[J]. 红外技术, 1996, 18 (5): 15-18. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS605.004.htm SU X, LI Z, LIU C, et al. Optical constants of a Selenide glass and its AR coatings. Infrared Technology, 1996, 18 (5): 15-18. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS605.004.htm
[10] 李栋, 齐晗兵, 吴国忠. 柴油的透射光谱测量和热辐射物性参数反演[J]. 光谱学与光谱分析, 2015, 35(3): 719-723. DOI: 10.3964/j.issn.1000-0593(2015)03-0719-05 LI D, QI H, WU G. Transmittance spectra measurement and thermal radiative physical parameters inversion of diesel fuel[J]. Spectroscopy and Spectral Analysis, 2015, 35(3): 719-723. DOI: 10.3964/j.issn.1000-0593(2015)03-0719-05
[11] 李栋, 艾青, 夏新林. 液态碳氢燃料热辐射物性参数反演方法[J]. 航空动力学报, 2012, 27(8): 1712-1717. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201208006.htm LI D, AI Q, XIA X. Inverse method investigation of thermal radiation property of liquid hydrocarbon fuel[J]. Journal of Aerospace Power, 2012, 27(8): 1712-1717. https://www.cnki.com.cn/Article/CJFDTOTAL-HKDI201208006.htm
[12] 李栋, 艾青, 夏新林. 透射法测量半透明液体热辐射物性的双厚度模型[J]. 化工学报, 2012, 63(S1): 23-129. https://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ2012S1022.htm LI D, AI Q, XIA X. Double-thickness model of thermal radiation physical property measurement of semi-transparent liquid with transmission method[J]. CIESC Journal, 2012, 63 (S1): 123-129. https://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ2012S1022.htm
[13] 王程超, 李兴灿, 谭建宇, 等. 生物柴油光学常数的双光程法实验测量[J]. 激光与光电子学进展, 2015, 52(5): 051206. https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201505019.htm WANG C, LI X, TAN J, et al. Experimental measurement of optical constant of biodiesel by double optical pathlength transmission method[J]. Laser & Optoelectronics Progress, 2015, 52(5): 051206. https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201505019.htm
[14] 李栋, 艾青, 夏新林. 利用透射光谱反演硒化锌的光学常数[J]. 光谱学与光谱分析, 2013, 33(4): 930-934. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201304019.htm LI D, AI Q, XIA X. optical constants determination of Zinc Selenide by Inversing transmittance spectrogram. Spectroscopy and Spectral Analysis, 2013, 33 (4): 930-934. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201304019.htm
[15] 李栋, 夏新林, 艾青. 两种反演半透明液体光学常数的方法对比[J]. 哈尔滨工业大学学报, 2012, 44(9): 73-77. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX201209015.htm LI D, XIA X, AI Q. Comparison of two inversion methods on optical constants of semitransparent liquid. Journal of Harbin Institute of Technology, 2012, 44 (9): 73-77. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX201209015.htm
[16] 吴国忠, 胡鑫浩, 齐晗兵, 等. 基于粒子群算法反演液体光学常数新方法[J]. 光学与光电技术, 2018, 16(4): 23-29. https://www.cnki.com.cn/Article/CJFDTOTAL-GXGD201804005.htm WU G, HU X, QI H, et al. A new method based on particle swarm optimization to invert liquid optical constants [J]. Optics & Optoelectronic Technology, 2018, 16 (4): 23-29. https://www.cnki.com.cn/Article/CJFDTOTAL-GXGD201804005.htm
-
期刊类型引用(2)
1. 付沛,崔岚,李硕. 基于高光谱成像的光敏印油种类区分实验. 中国无机分析化学. 2024(06): 836-841 . 百度学术
2. 李硕,崔岚,付沛. 基于高光谱成像结合分光光度技术的喷墨打印墨水种类鉴别方法. 中国无机分析化学. 2024(06): 826-835 . 百度学术
其他类型引用(0)