非均匀积污下高压绝缘子电位及电场分布特性
1)
苏州大学机电工程学院, 江苏苏州 215021
DISTRIBUTION CHARACTERISTICS OF ELECTRIC POTENTIAL AND ELECTRIC FIELD ALONG HV INSULATORS UNDER NON-UNIFORM CONTAMINATION
1)
School of Mechanical and Electric Engineering, Soochow University, Suzhou 215021, Jiangsu, China
通讯作者:2)高强, 教授, 主要研究方向: 环境流体力学及应用力学. E-mail:gaoqiang@suda.edu.cn
收稿日期:2020-01-10接受日期:2020-04-21网络出版日期:2020-05-18
基金资助: |
|
Received:2020-01-10Accepted:2020-04-21Online:2020-05-18
作者简介 About authors
高压绝缘子在输变电工程中承担着绝缘的重要角色. 由于工业化的飞速发展, 雾霾和沙尘等天气的频繁出现, 污秽颗粒在大气流场的作用下在绝缘子表面沉积, 容易发生污闪事故, 严重影响着电网的安全运行. 将高压绝缘子污闪发生的主要过程-绝缘子的动态积污和表面电场畸变过程进行耦合分析, 探究污闪发生的机理. 建立绝缘子周围的欧拉-欧拉气固两相流模型, 对绝缘子动态积污过程进行数值模拟, 获得绝缘子表面的非均匀积污层分布; 创建其表面导电层单元, 结合有限元方法建立积污绝缘子电场模型, 对非均匀积污层分布的绝缘子串进行电位及电场分析. 研究结果表明: 在积污层非均匀分布情况下, 绝缘子表面积污层周向位置对其电位分布的趋势影响较小, 但对绝缘子表面电场畸变的发生位置影响较大. 随着周向位置的变化, 电场畸变位置的偏移也不同, 且侧风面电场畸变的偏移最大. 随着积污层电导率的增大, 电场强度大于清洁绝缘子表面的场强, 且背风侧的场强畸变最大.
关键词:
HV insulators take on an important role in power transmission engineering. Because of the rapid development of industrialization, sandstorms and smog often happen. Contamination particles deposit on insulator surface under air flow field, pollution high voltage flashover accidents are prone occur, which seriously affects the safe operation of the power grid. In this paper, to investigate the mechanism of pollution high voltage flashover, a coupling analysis of dynamic accumulation and electric field distortion on insulator surface is explored. A model of gas-solid two phase flow around insulator is established by Euler-Euler approach to simulate the dynamic pollution accumulation process. Thus the distribution of non-uniform contamination layer on the insulator surface is obtained. The surface conductive layer element of the contamination layer is created, and it is used in finite element model of electric field on the polluted insulator, and the electric potential and electric field of the insulator string with non-uniform contamination are analyzed. The results show that: in the case of non-uniform contamination layer, the circumferential position of the polluted layer on the insulator surface has less influence on the trend of its potential distribution, but has greater influence on the position of electric field distortion. With the variation of circumferential position, the deviation of electric field distortion position is different, and maximum deviation occurs on crosswind. The electric field intensity is larger than clean insulator surface with increase of contamination layer conductivity. Furthermore, the field intensity distortion on leeward is the largest.
Keywords:
本文引用格式
毛东, 唐秋明, 高强.
Mao Dong, Tang Qiuming, Gao Qiang.
引言
近年来, 随着我国工业化进程的快速推进, 高压和特高压等电网的安全运行已成为国民经济平稳发展不可或缺的基础保障[1-2]. 由于工业化的发展导致了一系列环境问题, 其中包括沙尘及雾霾等天气频繁出现. 当这些极端天气出现时, 空气中的污秽颗粒在高压绝缘子表面沉积后形成积污层, 导致其绝缘性能降低, 易发生污闪事故[3-5], 导致电网瘫痪, 严重影响工业的生产和民众的生活[6-9]. 通过对绝缘子的积污形成机理以及积污层在高压电场下的电位及电场分布的分析, 研究绝缘子污闪发生的机理, 有望为超、特高压输电事业的顺利发展和国家电网系统的安全运行提供有效的保障和设计依据, 具有重要的研究意义和应用价值.
综上, 未见有从绝缘子的非均匀积污层的形成到其非均匀电场分析的这两大导致污闪发生的主要过程进行耦合研究的文献报道. 因为首先无论是从人工模拟还是野外观测等实验的角度等均难以实现积污层的量化和相应电场的测试, 其次在进行CFD模拟动态积污时, 通常采用的是欧拉-拉格朗日方程, 由于颗粒为离散相, 无法获得有效的积污层空间分布. 本研究采用欧拉-欧拉方程, 运用有限体积法, 有效地获取污秽颗粒在绝缘子表面形成的非均匀积污层; 运用表面导电层单元法创建积污层单元并作为边界条件, 在此基础上建立积污绝缘子的电场的有限元模型, 求解污秽绝缘子表面的电准静态场, 对非均匀积污层分布情况下的绝缘子表面电位及电场分布进行研究, 探讨污闪发生的机理.
1 数学模型
1.1 欧拉-欧拉气固两相流模型
式中, $V_{c}$是控制体单元的体积, $V_{s}$是污秽颗粒相的体积, $V_{g}$是空气相的体积, $\alpha_{s}$是指污秽颗粒相的体积分数.
空气相和污秽颗粒相均需满足的质量守恒方程如下
式中, $v_{g}$和$v_{s}$分别是空气相和污秽颗粒相的速度, $\rho_{g}$和$\rho_{s}$分别是空气相和污秽颗粒相的密度.
空气相和污秽颗粒相同时需满足动量守恒方程
式中, $K_{gs}=K_{sg}$是空气相与污秽颗粒相的动量交换系数, $p_{s}$为污秽颗粒相固体压力, $\tau_{s}$和$\tau_{g}$分别是污秽颗粒相与空气相的应力张量, $p$为相间的共享压力, $F_{qs}$是作用在污秽颗粒上的外力(例如电场力等), $g$为重力加速度.
图1
图2
积污计算中将数值风洞的入口边界条件设为速度进口, 出口边界条件设为充分发展出流, 绝缘子表面为无滑移壁面边界, 污秽颗粒一旦与绝缘子表面发生碰撞就认为污秽颗粒将直接吸附在壁面上. 为减少数值风洞壁面对内部流场的影响, 将数值风洞外壁面设为对称边界. 相关研究表明绝缘子表面的污秽颗粒的粒径多集中在0 $\sim$ 50 $\mu$m[24], 大气污染物的质量主要集中在PM10[25-26], 因此计算中将污秽颗粒近似地等效为球体, 粒径设为10 $\mu$m, 密度为2200 kg/m$^{3}$. 因空气中污秽颗粒的体积分数较小, 设为稀相, 建立 Syamlal-Obrien方程描述污秽颗粒间的颗粒黏度. 采用Wen-Yu曳力模型来描述稀疏污秽颗粒相与空气相间的动量交换系数.
1.2 积污绝缘子电场模型
式中, $E$为电场强度, $H$为磁场强度, $D$为电位移矢量或电通密度, $J$为电流密度, $\rho_{e}$为电荷密度, $B$为磁感应强度.
电网运行电压下绝缘子的电位分布应满足泊松方程
式中, $\varphi$为电势, $\rho_{c}$为自由电荷的空间密度, $\varepsilon_{c}$为介电常数. 当$\rho_{c}=0$时, 泊松方程变为关于电位的拉普拉斯方程
则电场强度$E$与电势$\varphi $的关系为
假定绝缘子表面的导电积污层为各向同性介质, 则在谐波状态下导电积污层的控制方程可表示为
式中, $J_{V}$为积污层中的体电流密度, $\omega$为运行电压的角频率.
非均匀积污层的任一区域$\varOmega$的厚度$d(\varOmega$)的面电导率$\sigma_{S}$可以用相应的体积电导率$\sigma$ (即通称的电导率)折算得到
式中, $D_{s}=d(\varOmega )D$为面电通密度, $J_{S}= d(\varOmega )J_{V}$为面电流密度. $\rho_{S}=\rho_{e}d(\varOmega)$为面电荷密度, 其中$\rho_{e}$为体电荷密度. 用电位$\varphi$可将式(12)表示为
式中, $\varepsilon$为介电常数, $\varepsilon_{S}= d(\varOmega)\varepsilon$为面介电常数, 将$\sigma $和$\varepsilon$代入式(12)中可以得到下式
2 结果与分析
2.1 非均匀积污层分布
为分析绝缘子表面污秽颗粒的沉积分布情况, 利用欧拉双流体模型进行积污计算. 设置颗粒相入口处体积分数为0.04[25], 风速以6 m/s水平吹向绝缘子, 污秽颗粒的粒径设为10 $\mu $m. 气流携带体积分数为$\alpha_{in}=0.04$的污秽颗粒流向绝缘子, 在与绝缘子发生碰撞, 污秽颗粒因增加附作用沉积形成积污层. 由于污秽颗粒在绝缘子表面沉积导致壁面相邻单元的体积分数高于来流体积分数. 因此假设当绝缘子表面的体积分数大于来流体积分数时, 形成积污层, 体积分数小于来流体积分数时, 绝缘子表面不会形成积污层. 即污秽颗粒的体积分数可直接表征绝缘子表面污秽颗粒的沉积分布状况. 在风速为6 m/s的情况下, 污秽颗粒的体积分数在绝缘子串表面上的分布情况如图3所示. 绝缘子串中各个绝缘子下伞面污秽颗粒体积分数如图4所示, 绝缘子串中绝缘子从上到下依次编号为1, 2和3.
图3
图4
图4各片绝缘子下表面颗粒相体积分数
Fig.4Particle volume fraction on lower surface of each insulator
从图3中可以看出, 在迎风面上伞面的边缘区域和被绝缘子钢帽遮挡形成的背风区域, 污秽颗粒的体积分数较大, 积污较为严重. 在绝缘子下表面的背风面污秽颗粒的体积分数也较大. 第1片绝缘子上表面的体积分数小于其余2片绝缘子上表面的体积分数. 但其下表面的体积分数与第2片绝缘子的体积分数并无较大差异, 与第3片绝缘子下表面的体积分数差异最大. 第1片绝缘子上表面背风区只受到其自身钢帽的遮挡, 该区域较难形成强烈的涡流, 并无明显的颗粒沉积. 而第2片和第3片绝缘子上表面背风区因自身钢帽及上一片绝缘子下表面的共同影响, 在该区域形成了强烈的涡流, 污秽颗粒沉积明显, 造成第2片与第3片绝缘子的上表面积污大于第1片绝缘子. 同样第3片绝缘子下表面的体积分数与第1片和第2片绝缘子下表面也有较大差异, 这是因为在第2片绝缘子下表面由于存在第3片绝缘子钢帽的遮挡, 下表面背风区棱与棱间形成了涡流, 增加了积污. 而第3片绝缘子下表面仅在棱边缘上有部分积污, 其积污量相对较少. 综上, 在绝缘子上表面迎风区域和下表面正对着来流方向的区域, 污秽颗粒的碰撞对积污起主要作用, 即污秽颗粒与伞裙表面发生碰撞或与已沉积的污秽颗粒发生碰撞. 其余部位的积污主要是由于涡流携带着污秽颗粒以较小的速度撞击绝缘子表面从而形成污秽颗粒沉积.
通过污秽颗粒的体积分数, 在FLUENT 15.0中使用UDM自定义积污层厚度变量, 并使用C语言编写UDF函数, 实现绝缘子表面积污层厚度$d(\varOmega)$的计算
式中, $\rho _{SL} $为积污层的密度, $\rho _{SL} =(1-\varepsilon _{V})\rho _{S} $, 积污层中污秽颗粒的堆积密度$\varepsilon _{V}=0.45$[20], $\rho_{s}$为污秽颗粒的材料密度, $S_{L}$为$\varOmega$的积污面积. $m_{s}$为沉积颗粒的质量, $m_{s}=(\alpha_{c}-\alpha_{in})V_{c}\rho_{s}$, $S_{L}$为绝缘子表面微元的面积, $\alpha_{c}$是绝缘子表面污秽颗粒的体积分数. 当$\alpha_{c}<\alpha_{in}$时, $d(\varOmega)=0$. 计算得到的积污层厚度分布如图5所示.
图5
图5绝缘子表面污秽沉积厚度分布(单位: mm)
Fig.5Thickness of contamination layer on insulator surface (unit: mm)
人工风洞积污试验[25]显示积污在绝缘子表面迎风面和背风面分布较多, 而侧风面几乎没有积污存在, 计算的积污层的分布情况与其基本吻合, 说明计算所得的绝缘子表面污层厚度分布趋势比较正确地反映了试验中绝缘子表面的积污分布规律, 验证了积污模型具有一定的合理性, 计算结果与文献的试验结果在积污分布区域的大小上的一些差异可以认为是绝缘子和试验条件的略有不同导致的, 特别是上伞裙倾角、下伞裙结构形状和颗粒的粒径分布.
2.2 积污绝缘子电场分析
建立绝缘子表面非均匀积污层电场分布模型, 对其电位及电场分布情况进行数值计算分析. 以积污层电导率$\sigma =5.0\times 10^{-4}$为例[11], 对积污后的电场分布等进行阐述.
当均匀积污层覆盖绝缘子表面时, 其表面电位分布如图6(a)所示. 由于绝缘子表面迎风面、侧风面和背风面积污层的厚度各不相同, 导致绝缘子表面不同周向位置处的电位降存在差异, 但整体而言趋势影响并不是很大, 即周向位置的差异对其电位分布的影响不是很大. 如图6(b)所示, 当绝缘子表面积污层分布不均匀时, 不同周向位置处的沿面电场也发生了较大的畸变, 电场分布变得极不均匀. 迎风面和背风面沿面最大场强都出现在高压侧, 最大值分别为51.37 kV/cm和24.56 kV/cm. 绝缘子表面积污层的不均匀分布使得绝缘子串中每片绝缘子沿面电场强度的极值出现的位置与其周向位置有关, 侧风面相对于迎风面与背风面发生了明显的偏移, 据此可以推断出绝缘子表面放电的位置也随之偏移.
图6
图6非均匀积污层分布下沿面电位和电场分布($\sigma =5.0\times 10^{-4}$ S/m)
Fig.6Electric field strength and potential distributions for non-uniform ($\sigma =5.0\times 10^{-4}$ S/m)
图7(a)是非均匀积污层分布条件下绝缘子串表面场强分布云图. 由图7(b)和图7(c)可知, 在清洁和均匀积污条件下, 绝缘子串的沿面电场沿周向均匀一致分布. 由图8可知, 在清洁条件下, 场强最大值发生在钢脚处; 在积污层均匀覆盖条件下, 场强值在一定范围内波动. 而在非均匀分布条件下图7(a)所示, 绝缘子串沿面电场分布极不均匀, 背风面及侧风面处的场强较大, 而迎风面的场强较小. 同一绝缘子表面在背风面与侧风面积污层交界处场强较大, 这是由于导电积污层的存在, 使得两侧的电压降相差较大, 从而在该处具有较大的场强. 非均匀积污层分布条件下, 绝缘子串中不同片绝缘子的电场强度分布也不相同.
图7
图7不同积污层分布条件下绝缘子串表面电场分布(单位: V/m)
Fig.7Electric field strength distributions for different condition (unit: V/m)
图8
图8绝缘子串沿面场强分布(单位: V/m)
Fig.8Electric field strength distributions along an insulator string (unit: V/m)
当积污层的湿润状态不同时, 其电导率也不相同. 不同积污层电导率下绝缘子不同位置处的沿面电场分布情况如图9所示. 从图中可以看出, 随着积污层电导率的增大, 迎风面、背风面和侧风面的平均场强均大于清洁绝缘子表面的平均场强, 且背风侧的场强畸变幅度最大. 同时可以看出随着积污层电导率的变化, 绝缘子表面场强畸变的位置也随之迁移变化, 其中背风面, 电场畸变位置迁移距离最大. 当积污层的电导率为$5.0\times 10^{-3}$ S/m时, 背风侧的场强最大位置由钢脚处偏移到上表面伞裙边缘, 且最大场强出现在高压端处的绝缘子上表面伞裙边缘. 这是因为在背风面伞裙边缘沉积的污秽较少, 使得该区域积污层的电导率几乎为0.
图9
图9不同积污层电导率下绝缘子串不同位置处的沿面电场分布
Fig.9Electric field strength distributions along the insulator string under different conductivity
通过对非均匀积污条件下绝缘子串沿面电位和电场分布的计算, 发现绝缘子串表面由于非均匀积污层分布使得其沿面电场电位分布发生了极大的变化, 特别是背风面的电场畸变幅度最大. 根据以上研究, 可以发现随着积污层对绝缘子表面电位及电场的影响, 使得绝缘子表面电场发生严重畸变, 当电场强度超过空气的击穿场强(33 Kv/cm)时, 绝缘子表面就会形成放电通道, 从而引发绝缘子污闪.
3 结论
将引起绝缘子污闪事故的两大过程--绝缘子积污和表面电场畸变, 进行耦合分析污闪发生的机理. 通过建立欧拉-欧拉气固两相流数值模型, 对绝缘子的积污过程进行仿真模拟, 得到非均匀积污层后, 创建积污层的表面导电层单元, 解决了积污层太薄而导致的网格畸变问题, 从而建立了有效的积污绝缘子电场的有限元分析模型, 对非均匀积污层分布的绝缘子串进行电位及电场分析, 探讨了绝缘子表面积污层周向位置对电位及电场分布的影响以及电场畸变引发绝缘子污闪的情况. 结论如下:
(1) 在绝缘子上表面迎风区域和下表面正对着气流来流方向的区域, 污秽颗粒的碰撞对积污起主要作用. 其余部位的积污主要是由于涡流携带着污秽颗粒以较小的速度撞击绝缘子表面从而形成污秽颗粒沉积.
(2) 在积污层非均匀分布情况下, 绝缘子表面积污层周向位置对其电位分布的趋势影响较小, 但对绝缘子表面电场畸变的发生位置影响较大. 周向位置的不同, 电场畸变位置的偏移也不同, 且背风面的电场畸变的偏移位置最大.
(3) 随着积污层电导率的增大, 迎风面、背风面和侧风面的平均场强均大于清洁条件下绝缘子表面的平均场强, 且背风侧的场强畸变幅度最大. 同时可以看出随着积污层电导率的变化, 绝缘子表面场强畸变的位置也随之迁移变化, 其中背风面电场畸变位置迁移位移最大.