刘 刚,徐 庆,石聪聪,张圆明,黄华峰
(国网江苏省电力有限公司宿迁供电公司,江苏 宿迁 223800)
目前我国大部分城区新建配电线路多采用配电全电缆网络,已有配电架空线路也在逐步改造过渡为架空电缆混合网[1]。其中,环网柜因具有维护工作量小、安装方式灵活、尺寸小等优点,在市辖区配电网络中得到了广泛的应用[2]。然而环网柜运行时仓盖封闭,加上其数量大、安装位置繁多,不便于巡视检查,给运维检修工作带来了很大的难度[3]。
近年来,由于巡检不及时且无科学有效的在线监测手段,我国电网辖区发生了多起由环网柜T 型电缆终端接头温度过热、后堵盖脱落、爬电而导致的严重劣化烧蚀故障,从而造成相间短路、单相接地等电网故障[4],严重时甚至引发环网柜灼烧、局部区域停电,给电力系统的安全稳定运行造成不利影响。究其原因,大部分是由绝缘交界面气隙缺陷导致[4-5]。
环网柜空间狭小、不利于散热,T 型终端接头施工工艺复杂且质量难以保证,造成了以下问题:①长期运行在强电场、弱对流的环境下,接头内部热胀冷缩、表面结垢、氧化或腐蚀,导致绝缘交界面松动接触不良、出现气隙,引起发热、漏电,加剧电缆终端绝缘层的劣化速度[5-6];②由于生产和安装的问题,部分环网柜电缆接头处本身就存在一定缺陷,导致接触电阻、弯曲应力过大,在长期的热老化和机械老化作用下,接头根部松动开裂,产生气隙,最终导致内外屏蔽层击穿[7-8]。
研究电缆绝缘及其接头附件在复杂环境下的电场、温度场分布特性,可以揭示其故障缺陷的形成机理和发展过程,为电缆接头故障预防及抑制提供依据,为优化绝缘材料配方、设计安装结构、开展故障检测及诊断提供重要的理论支撑。为此,学者们开展了一定的研究[9-15],文献[12]利用有限元仿真研究了无应力管、有空气隙、有金属微粒附着及钢针扎入4 种缺陷下的电缆头电场分布,结果发现电缆终端电场分布受应力管影响较大,金属微粒会导致明显的局部电场畸变,针尖缺陷的尖端处会因电场急剧增大而导致绝缘放电。文献[13]考虑了空间电荷积累及直流配网拓扑结构,仿真分析了10 kV三芯绝缘电缆截面的温度场和电场分布,发现电缆绝缘层上的电场强度沿半径方向逐渐降低,电场强度的最大值出现在导体表面,而当电缆绝缘温差大于10℃时,电场分布发生反转,电场强度最大值出现在绝缘外表面。文献[14]考虑了直流电压、直流叠加冲击电压的作用,仿真分析了不同温度梯度作用下高压直流电缆接头的电场分布,得到其最大场强随着温度的升高而增大,而高温、高电场时电缆接头会出现局部电场畸变。文献[15]仿真分析了不同导电性硅橡胶下的XLPE 电缆中间接头电场分布特性,提出利用电导率较大的非线性硅橡胶基复合材料作为中间接头绝缘材料,可以改善电缆中间接头的电场分布。
上述研究对掌握高压电缆运行特性及故障机制有参考意义。但是目前针对T型电缆终端接头绝缘交界面缺陷特性的研究较少。为获得T型电缆终端劣化的电场和温度场特性,进而提出合理有效的监测识别手段,本研究建立环网柜用10 kV T 型电缆终端接头有限元模型,研究绝缘交界面典型气隙缺陷下,T型终端接头的电场、温度场幅值变化规律及其在电缆室不同位置处的分布特性,并基于此提出可有效反映T型电缆终端接头劣化的特征量及监测位置。
研究T 型电缆终端接头的电场与温度场分布时,运行工况为50 Hz 工频电压,实际运行中,传导电流密度为位移电流密度的107倍,因此位移电流可以忽略不计。对于T 型电缆终端接头,沿电缆终端接头劣化界面的电场和电位分布受位移电流的影响较大。因此,对于电场的计算,本研究选取准静电场计算模型。
由静电场原理,T 型电缆终端接头整体电位分布满足泊松方程,如式(1)所示。当电场中无自由移动的空间电荷时,ρ(静电场中任意一点的电荷密度)为0,拉普拉斯方程成立,如式(2)所示。
式(1)~(2)中:ε为介质的介电常数;φ为电位;∇2为拉普拉斯算子。
求解T 型电缆终端接头的整体电场分布时,静电场变分公式和计算域单元顶点的电位(Fe(φe))表达式分别为式(3)和式(4)所示。
式(3)~(4)中:e表示有限元计算单元;J表示计算场域;V表示计算曲面包围的空间体积。
Fe(φe)对φe导数为0,则可以得到式(5)。进一步表示为矩阵的形式,如式(6)所示。
式(6)中,[K]为刚度矩阵,通过求解器的迭代计算,最终可以求得T型电缆终端接头的整体电场分布。
为研究T 型电缆终端接头的电热耦合问题,根据麦克斯韦方程,引入矢量磁位Ȧ,对于有激励源存在的电流区的磁矢量位控制方程为式(7),对于无激励源的非电流区的磁矢量位控制方程为式(8)。
式(7)~(8)中:ω为角频率;μ为磁导率;σ为电导率;J̇s为外加电流密度。
T 型电缆终端接头磁场问题,可归纳为式(7)和
式(8)的求解问题,本研究中,给出4 类求解的边界条件,其中,第一类边界条件,其磁力线垂直于边界面,表达式可写为式(9),其中n为单位法向矢量;第二类边界条件,其磁力线平行于边界面,表达式可写为式(10);第三类边界条件,具有面电流密度δ̇s,如式(11)所示,则表达式可写为式(12)。最后一类边界条件,给定矢量磁位具体数量值,如式(13)所示,当此边界离电流源较远时,则边界条件可以近似简化为式(14)。
式(9)~(14)中:Γ1、Γ2、Γ3、Γ4分别代表4类边界条件;Ḣ为磁场强度。
对于T 型电缆终端接头的温度场计算,T 型电缆终端接头属于固体介质,其温度场满足的导热微分方程在直角坐标系中可写为式(15)。
式(15)中:ρ'表示微元体密度;c是比热容;T为温度;t是时间;qv为体积产热率;λ为导热系数。
本研究的T型电缆终端接头的温度场属于稳态温度场,而稳态导热微分方程可以简化为式(16)。
根据传热学理论,通常传热条件有以下3种:①已知边界的温度分布情况,一般采用恒温边界条件,表示为式(17),T0为初始已知边界条件温度;②已知流体温度及对流传热情况的边界条件,控制方程为式(18),Tf为流体温度,h为对流换热系数;③已知周围环境,且产热和散热达到平衡的情况下,其控制方程可以写为式(19),其中Th为环境温度。
本文研究的是T型电缆终端接头在空气中的电场及热场分布,属于以式(19)为控制方程的边界条件,此外,根据斯蒂芬-玻尔兹曼定律,T型电缆终端接头表面的散热边界条件可以表示为式(20)。
式(20)中,σ0、ε1分别为斯蒂芬-玻尔兹曼常数及表面发射率。
此外,电缆接头的电磁场计算和温度场计算是一个双向耦合过程。
仿真计算中,以实际10 kV T 型电缆终端接头为研究对象,设计参数为6/10 kV,其结构剖面图如图1 所示,本研究综合考虑了以下工况:T 型电缆终端接头尾塞、套管与T型接头的绝缘交界面有气隙,以及气隙严重程度对其电场和温度场的影响。
图1 T型电缆终端接头剖面图Fig.1 Section of T-type cable joint
参考文献[12,16-17]的数据,本研究有限元计算的材料物性参数设置如表1所示。
表1 有限元计算的材料物性参数Tab.1 Material physical parameters for finite element calculation
通过几何建模的方式,在交界面建立截面为1 mm×30 mm 和2 mm×50 mm 的空气气隙分别模拟轻度气隙和中度气隙,以轻度气隙为例,模拟结果如图2(a)所示。本研究计算过程中,施加电压为10/kV,根据有限元原理的计算方法,先对T 型电缆终端接头的电场和温度场进行计算。其中,本研究设置空气计算域为环境温度,电缆头与空气接触面为对流传热边界。计算完成后,选取交界面、箱体上表面、法兰底座、上护套外部4 处典型位置,分析T 型电缆终端接头的电场与温度场分布特性,如图2(b)所示。
图2 T型电缆终端接头电场与温度场仿真设置Fig.2 Simulation setting of electric field and temperature field of T-type cable joint
2.1.1 无缺陷整体电场分布
无缺陷T型电缆终端接头整体的电场分布如图3 所示。从图3 可以看出,T 型电缆终端接头在接地眼处电场强度偏大,交界面电场强度次之,但无缺陷的情况下,交界面电场强度没有明显的异常增大情况。
图3 无缺陷整体电场分布Fig.3 Electric field distribution without defect
根据不同路径的电场分布结果,交界面电场、箱体上表面电场、法兰底面电场和上护套外部电场计算结果如图4 所示。从图4 可以看出,交界面电场范围为43.07~162.23 kV/m,箱体上表面电场范围为3.03~30.37 kV/m,法兰底面电场范围为7.00~61.10 kV/m,上护套外部电场范围为57.96~218.65 kV/m。其中,交界面电场和上护套外部电场最大值较大,上护套外部电场最大值可达218.65 kV/m,由此可得,交界面和上护套外部对绝缘要求更高。
图4 无缺陷下不同位置的电场分布仿真结果Fig.4 Simulation results of electric field distribution at different positions without defects
2.1.2 轻度气隙缺陷整体电场分布
T 型电缆终端接头在轻度气隙缺陷下的整体电场分布如图5 所示。从图5 可以看出,整体电场分布情况与无缺陷情况相似。
图5 轻度气隙缺陷下整体电场分布Fig.5 Electric field distribution with slight air gap defect
根据不同路径的电场分布结果,在轻度气隙缺陷下交界面电场、箱体上表面电场、法兰底面电场和上护套外部电场计算结果如图6 所示。从图6 可以看出,交界面电场范围为41.21~197.77 kV/m,箱体上表面电场范围为3.03~30.72 kV/m,法兰底面电场范围为7.00~61.10 kV/m,上护套外部电场范围为56.56~221.09 kV/m。其中,交界面电场和上护套外部电场最大值仍然较大,上护套外部电场最大值可达221.09 kV/m。此外,电场沿仿真路径均不是均匀分布,与无缺陷情况下较为相似。
图6 轻度气隙缺陷不同路径电场分布结果Fig.6 Simulation results of electric field distribution at different positions with slight air gap defect
T 型电缆终端接头在轻度气隙缺陷与无缺陷情况下电场分布的主要区别为:交界面电场最大值从162.23 kV/m 增大至197.77 kV/m,增加了21.91%;上护套外部电场最大值从218.65 kV/m 增大至221.09 kV/m。而箱体上表面电场和法兰底面电场变化较小。
2.1.3 中度气隙缺陷整体电场分布
T 型电缆终端接头在中度气隙缺陷下的整体电场分布如图7 所示。从图7 可以看出,电场分布亦与无缺陷电场分布相似,T 型电缆终端接头在接地点处电场强度仍偏大,交界面电场强度次之。
图7 中度气隙缺陷下整体电场分布Fig.7 Electric field distribution with moderate air gap defect
根据不同路径的电场分布结果,在中度气隙缺陷下交界面电场、箱体上表面电场、法兰底面电场和上护套外部电场计算结果如图8 所示。从图8 可以看出,交界面电场范围为39.72~235.32 kV/m,箱体上表面电场范围为3.03~30.73 kV/m,法兰底面电场值范围为7.00~61.10 kV/m,上护套外部电场范围为55.42~220.10 kV/m。其中,交界面电场和上护套外部电场最大值仍然较大,交界面电场最大值可达235.32 kV/m。此外,电场沿几个仿真路径均不是均匀分布,与无缺陷和轻度缺陷情况下较为相似。
图8 中度气隙缺陷下不同路径电场分布结果Fig.8 Simulation results of electric field distribution at different positions with moderate air gap defect
T 型电缆终端接头在中度气隙缺陷与无缺陷情况下电场分布的主要区别为:交界面电场最大值从162.23 kV/m 增大至235.32 kV/m,增加了45.05%;上护套外部电场最大值从218.65 kV/m 增大至220.10 kV/m。但箱体上表面电场和法兰底面电场变化仍然较小。
2.2.1 无缺陷整体温度场分布
无缺陷T型电缆终端接头整体的温度场分布如图9所示。从图9可以看出,T型电缆终端接头在接地眼处温度高,电缆头本体温度次之,环网柜箱体温度最低。无缺陷的情况下,交界面温度场强度没有明显的异常增大情况。
图9 无缺陷整体温度场分布Fig.9 Temperature field distribution without defect
交界面温度场、箱体上表面温度场、法兰底面温度场和上护套外部温度场计算结果如图10所示。从图10 可以看出,交界面温度范围为30.75~31.12℃、箱体上表面温度范围为25.17~29.60℃,法兰底面温度范围为29.59~30.95℃,上护套外部温度范围为30.61~31.20℃。其中,交界面温度场和上护套外部温度场最大值较大,上护套外部温度场最大值达31.20℃。此外,温度场沿监测路径均不是均匀分布,交界面温度场沿法兰底座方向监测路径温度降低,箱体上表面温度沿法兰底座方向监测路径温度升高,法兰底面温度沿监测路径先升高后降低,上护套外部温度场沿法兰底座方向监测路径温度降低。
图10 无缺陷不同路径温度场分布结果Fig.10 Simulation results of temperature field distribution at different positions without defects
2.2.2 轻度气隙缺陷整体温度场分布
T 型电缆终端接头在轻度气隙缺陷下的整体温度场分布如图11所示。从图11可以看出,轻度气隙缺陷下温度场分布与无缺陷情况下温度场分布相似,T型电缆终端接头在接地眼处温度最高,电缆头本体温度次之,环网柜箱体温度最低。
图11 轻度气隙缺陷下整体温度场分布Fig.11 Temperature field distribution with slight air gap defect
根据不同路径的温度场分布结果,在轻度气隙缺陷下交界面温度场、箱体上表面温度场、法兰底面温度场和上护套外部温度场计算结果如图12 所示。从图12 可以看出,交界面温度范围为31.16~31.53℃,箱体上表面温度范围为25.23~30.01℃,法兰底面温度值范围为29.95~31.16℃,上护套外部温度值范围为31.01~31.60℃。温度场沿监测路径均不是均匀分布,与无缺陷情况下较为相似。
图12 轻度气隙缺陷不同路径温度场分布结果Fig.12 Simulation results of temperature field distribution at different positions with slight air gap defect
T 型电缆终端接头在轻度气隙缺陷下与无缺陷情况下电场分布的主要区别为:4 个分布路径上,最高温度值均有所上升,但幅度较小,温度最大增加值仅为0.41℃。
2.2.3 中度气隙缺陷整体温度场分布
T 型电缆终端接头在中度气隙缺陷下的整体温度场分布如图13 所示。从图13 可以看出,中度气隙缺陷下温度场分布亦与无缺陷情况下温度场分布相似,T 电缆终端接头在接地眼处温度场强度仍偏大,交界面温度场强度次之。
图13 中度气隙缺陷下整体温度场分布Fig.13 Temperature field distribution with moderate air gap defect
根据不同路径的温度场分布结果,在中度气隙缺陷下交界面温度场、箱体上表面温度场、法兰底面温度场和上护套外部温度场计算结果如图14 所示。从图14 可以看出,交界面温度范围为31.83~32.20℃,箱体上表面温度范围为25.69~30.61℃,法兰底面温度范围为30.54~32.03℃,上护套外部温度范围为31.67~32.25℃。温度场沿不同路径分布与轻度气隙缺陷和无缺陷情况下较为相似。
图14 中度气隙缺陷下不同路径温度场分布结果Fig.14 Simulation results of temperature field distribution at different positions with moderate air gap defect
T 型电缆终端接头在中度气隙缺陷下与无缺陷下温度场分布的主要区别为:4 个路径上,最大温度值亦均有所升高,但幅度也较小,温度最大增加值为1.10℃。此外,统计上述数据得到温度的变化幅度在0.68%~3.49%。
上述研究表明气隙状态下的绝缘交界面处最大电场强度不到3 000 kV/m,远达不到空气击穿强度,这一方面是由于10 kV 环网柜T 型接头施加的电压较低、电场强度不够,另一方面是气隙缺陷未严重到使内部发生局部放电。结合文献调研及现场分析,考虑到长期气隙受潮的影响,分析气隙、水膜、水珠并存时T 型电缆接头的多物理场分布特性[18-19]。
模拟气隙受潮下的多物理场分布,针对轻度气隙受潮,设置了23.5 mm 的水膜和4 个接触角为90°的水滴,水滴为半圆形,直径为3.5~4.5 mm,如图15 所示。其中,水的相对介电常数设置为81,导热系数取常温下的值(0.599 W/(m·K)),电导率设置为
图15 T型电缆终端接头气隙受潮模拟Fig.15 Simulation of T-type cable joint with air gap wetting
10 μS/cm。
气隙受潮情况下,不同典型路径下的电场分布仿真结果如图16 所示。从图16 可以看出,气隙受潮情况下,在水滴边缘出现电场强度最大值,超过了3 000 kV/m,达到了空气的临界击穿强度,会造成局部放电、沿面爬电,从而引发电缆头的劣化烧蚀故障。
从图16 还可以看到,交界面电场范围为45.31~3 431.98 kV/m,箱体上表面电场范围为5.81~38.12 kV/m,法兰底面电场范围为14.21~77.32 kV/m,上护套外部电场范围为99.39~261.75 kV/m。其中,交界面电场最大值较大,达到了3 431.98 kV/m。此外,电场沿监测路径分布与气隙缺陷情况下较为相似。
图16 气隙受潮缺陷下不同路径电场分布结果Fig.16 Simulation results of electric field distribution at different positions with air gap wetting defect
在气隙受潮、达到临界击穿的情况下,箱体、法兰、上护套处的最大电场强度增幅明显,其中箱体上表面最大电场增大约26%、法兰底座最大电场增大约27%、上护套外部最大电场增大约20%。由此可见,法兰底座、箱体上表面的电场幅值变化能够表征气隙受潮缺陷的临界情况,且变化幅度较大。
T 型电缆终端接头在气隙受潮缺陷下的温度场分布如图17 所示。从图17 可以看出,气隙受潮缺陷下交界面温度范围为46.58~47.45℃,箱体上表面温度范围为33.12~43.74℃,法兰底面温度范围为43.21~47.06℃,上护套外部温度范围为46.28~47.52℃。温度场沿几个典型路径分布和无缺陷情况下较为相似。
图17 气隙受潮缺陷下不同路径温度场分布结果Fig.17 Simulation results of temperature field distribution at different positions with air gap wetting defect
从图17 还可以看到,在气隙受潮、达到临界击穿的情况下,箱体、法兰、上护套处的最高温度有明显升高的趋势,其中箱体上表面最高温度值升高了约14.1℃、法兰座最高温度值升高了约16.1℃、上护套外部最高温度值升高了约16.3℃。综上,T 型电缆终端接头在气隙受潮缺陷下与无缺陷情况下温度分布的主要区别为:4 个典型路径上,最高温度值均有约15℃的温升。
(1)绝缘交界面处电场分布受气隙缺陷的影响最为显著,随着缺陷加剧,电场最大值呈显著增大趋势,轻度缺陷和中度缺陷较无缺陷时分别增大了21.91%和45.05%
(2)随着气隙增大,4 处典型位置处的最高温度值也呈上升趋势,但是相较于电场分布,温度值的变化幅度较小,升高幅度为0.68%~3.49%。
(3)在气隙受潮的情况下,气隙缺陷处会达到临界击穿,此时箱体上表面最大电场值增大约26%,法兰座最大电场值增大约27%,上护套外部最大电场值增大约20%。
(4)在气隙受潮、达到临界击穿的情况下,箱体、法兰、上护套处的最高温度明显升高,其中箱体上表面最高温度升高约14.1℃,法兰座最高温度升高约16.1℃,上护套外部最高温度升高增大约16.3℃。
(5)考虑到传感器的安装和实施,在实际工程中,可通过监测法兰底座、箱体上表面的电场幅值增量(约20%)以及温升(约15℃),来反映T 型电缆终端接头即将出现局部击穿放电的潜在临界隐患。