CN112464482A - 一种sf6电气设备中特征组份气体扩散效应分析方法 - Google Patents

一种sf6电气设备中特征组份气体扩散效应分析方法 Download PDF

Info

Publication number
CN112464482A
CN112464482A CN202011403631.XA CN202011403631A CN112464482A CN 112464482 A CN112464482 A CN 112464482A CN 202011403631 A CN202011403631 A CN 202011403631A CN 112464482 A CN112464482 A CN 112464482A
Authority
CN
China
Prior art keywords
gas
diffusion
concentration
equipment
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011403631.XA
Other languages
English (en)
Other versions
CN112464482B (zh
Inventor
张英
贺毅
张靖
王明伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guizhou Power Grid Co Ltd
Original Assignee
Guizhou Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guizhou Power Grid Co Ltd filed Critical Guizhou Power Grid Co Ltd
Priority to CN202011403631.XA priority Critical patent/CN112464482B/zh
Publication of CN112464482A publication Critical patent/CN112464482A/zh
Application granted granted Critical
Publication of CN112464482B publication Critical patent/CN112464482B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Sampling And Sample Adjustment (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种SF6电气设备中特征组份气体扩散效应分析方法,它包括:步骤1、建立气体绝缘全封闭组合电器室几何模型;步骤2、建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程;步骤3、设定初始条件;步骤4、分析SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散情况下的气体扩散效应;步骤5、对气体扩散效应的结果进行分析得出结论;解决了现有技术针对电气设备SF6气体进行分析,由于不了解气体组分在内部的扩散特性,导致设备故障诊断方面尚存在不足,严重阻碍了设备故障诊断技术的发展等技术问题。

Description

一种SF6电气设备中特征组份气体扩散效应分析方法
技术领域
本发明属于六氟化硫气体监测技术,尤其涉及一种SF6电气设备中特征组份气体扩散效应分析方法。
背景技术
六氟化硫(SF6)气体由于其优良的绝缘以及灭弧性能而被广泛应用于各类高压电气设备中。但是由于在设计、制造、安装和运行维护等诸方面存在的一些缺陷,设备故障时有发生。充装SF6的全封闭组合电器(GIS)在运行过程中曾多次发生故障,当设备内部发生放电或过热等故障时,SF6气体会分解产生SO2、HF、H2S等表征设备故障的特征组分气体。关于SF6特征组分气体的检测方法有两种,一种是在线监测,另一种是离线检测。在线监测是指安装在设备取气口处的特征组分气体在线测量装置对封闭体系下的特征组份气体进行在线监测,在线监测方法有红外光谱法、电化学传感器法等。离线检测则是将设备中的气体引出后再利用分解产物仪进行测量,气体组分检测技术曾发现多起设备故障,在保证设备安全运行方面起到重要的作用。
然而,由于不了解气体组分在内部的扩散特性,导致设备故障诊断方面尚存在不足,严重阻碍了设备故障诊断技术的发展。
发明内容
本发明要解决的技术问题是:提供一种SF6电气设备中特征组份气体扩散效应分析方法,以解决现有技术针对电气设备SF6气体进行分析,由于不了解气体组分在内部的扩散特性,导致设备故障诊断方面尚存在不足,严重阻碍了设备故障诊断技术的发展等技术问题。
本发明的技术方案是:
一种SF6电气设备中特征组份气体扩散效应分析方法,它包括:
步骤1、建立气体绝缘全封闭组合电器室几何模型;
步骤2、建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程;
步骤3、设定初始条件;
步骤4、分析SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散情况下的气体扩散效应;
步骤5、对气体扩散效应的结果进行分析,得出结论。
步骤1所述建立气体绝缘全封闭组合电器室几何模型为:模型分为两个部分,其中大圆柱为气体绝缘全封闭组合电器室GIS,圆柱直径为650mm,长为2000mm,上方有小圆柱取气口,直径为10mm,长为150mm;采用Ansys Fluent Meshing工具划分网格,网格数为12843个,网格质量均在0.8以上。
步骤2所述建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程的公式为:
质量守恒方程:
Figure BDA0002817872150000021
式中,ρ为气体密度,μj为湍流粘性系数的径向分量,t为时间,xj表示特征组分气体气体沿腔体径向的位移,Sm为源项;
动量守恒方程:
X方向:
Figure BDA0002817872150000022
Y方向:
Figure BDA0002817872150000031
Z方向:
Figure BDA0002817872150000032
u、v、w为速度矢量在x、y、z方向的分量,Su、Sv和Sw和是动量方程的广义源项;μ为湍流黏性系数;
能量守恒方程:
Figure BDA0002817872150000033
T为流体温度;kt为流体的湍流导热系数;cp为混合气体的定压比热;cpr为SO2定压比热;cpa为SF6的定压比热;ω为重气质量比例;为常数,常取1;μt为流体的湍流粘性系数;
组分输运方程:
SF6分解产生特征组分气体,需要利用组分输运方程计算局部每相的质量分数。由组分质量守恒定律可以得到气体扩散组分方程:
Figure BDA0002817872150000034
式中:ω为每一组分的质量分数,Dm为湍流扩散系数。
所述初始条件包括:气室内压力为0.4Mpa,室内的温度为300K,采用非耦合瞬态的方法进行求解,整个模型是封闭体系即取气口不打开,气体的扩散通过浓度差以及分子热运动进行;忽略设备内部母线或机构在空间上对扩散过程的影响,定义设备内部特征组份气体最大和最小浓度差小于等于1μL·L-1时,认为特征组份气体扩散达到均匀。
所述分析HF气体单独扩散情况下气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SF6气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的HF,使用FLUENT软件模拟HF气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,经过1000s的扩散,设备内部的HF的最大浓度由原来的5000μL·L-1下降到44.6μL·L-1,设备内部HF的浓度呈现一个从左上到右下、从大到小的一个浓度梯度,最大浓度在设备的左上方,浓度大小为44.6μL·L-1,最小浓度在设备的底部,浓度大小为24.4μL·L-1;经过3000s的扩散,设备内部的最大浓度继续下降,下降到36.5μL·L-1,最大浓度出现在取气口位置,最小浓度位于设备右下方,此时最大浓度为36.5μL·L-1,最小浓度为32.4μL·L-1,经过6420s的扩散,设备内部最大浓度为34.8μL·L-1,最小浓度为33.8μL·L-1,浓度差为1μL·L-1,认为扩散已经达到均匀;
当设备内部HF初始浓度分别为5000、2000、1000、500、200和100μL·L-1时,使用FLUENT软件模拟HF气体从开始扩散到扩散均匀所需要的时间,如下表1所示;得出在封闭体系中,HF气体仅通过浓度差以及分子热运动驱动扩散的情况下,达到扩散均匀所需要的时间分别为:6420s、5890s、5190s、4760s、2580s、200s。
表1不同HF浓度下扩散均匀的时间及最大最小浓度
Figure BDA0002817872150000041
Figure BDA0002817872150000051
所述SO2气体单独扩散情况下气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SO2气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的SO2,使用FLUENT软件模拟HF气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,当扩散时间为1000s时,SO2浓度由原来的5000μL·L-1下降到72.2μL·L-1,SO2气体呈现一个由下至上、由大到小的一个浓度梯度,最大浓度在设备底部,浓度为72.2μL·L-1,最小浓度出现在取气口处,浓度为0.121μL·L-1;当扩散到达中期,扩散时间为4800s,设备本体内最大浓度已经下降到35.4μL·L-1,设备底部有最大浓度,取气口顶部存在最小浓度,最小浓度为25.2μL·L-1;当扩散时间为9260s时,最大浓度为34.4μL·L-1,最小浓度为33.4μL·L-1,此时扩散达到均匀;
当设备内部SO2浓度分别为5000、2000、1000、500、200和100μL·L-1时,达到扩散均匀所需要的时间分别为9260s、7310s、6130s、4850s、3560s、2760s;
表2不同SO2浓度下扩散均匀的时间及最大最小浓度
Figure BDA0002817872150000052
Figure BDA0002817872150000061
分析H2S气体单独扩散情况下的气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SF6气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的H2S,使用FLUENT软件模拟H2S气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,使用FLUENT软件模拟H2S气体从指定位置扩散到整个设备内部的过程,H2S气体从扩散初期到中期以及到扩散均匀后的结果;
当扩散时间为1000s时,H2S气体的浓度由最开始的5000μL·L-1下降到83.5μL·L-1,此时设备内部呈现的浓度规律为设备右下方浓度最高,最高浓度为83.5μL·L-1,设备左上方浓度最低,最低浓度为0.0557μL·L-1,呈现一个从右下到左上、从大到小的一个浓度梯度;当扩散时间达到4800s时,此时的最大浓度下降到35.8μL·L-1,最小浓度为25.7μL·L-1。此时设备内部的H2S气体浓度呈现一个设备中部浓度高,设备左右侧浓度低的一个规律,浓度最低出现在取气口;当扩散时间为8880s,此时设备内部的最大浓度为34.6μL·L-1,最小浓度为最小浓度为33.6μL·L-1,扩散达到均匀;
当设备内部SO2浓度分别为5000、2000、1000、500、200、100μL·L-1时,其他条件不变,达到扩散均匀所需要的时间分别为8880s、7140s、6850s、6310s、4960s、3430s;
表3不同H2S浓度下扩散均匀的时间及最大最小浓度
Figure BDA0002817872150000071
步骤5所述对气体扩散效应的结果进行分析,得出结论的方法包括:
SO2的扩散初期由设备左边缘位置向右扩散,形成设备底部浓度高、顶部浓度低的浓度梯度。扩散达到中期,SO2气体依靠浓度差以及分子热运动由设备底部向顶部扩散,最终扩散均匀;
H2S气体与SO2气体的扩散相似,扩散初期由设备左边缘向右扩散,形成了设备右下方浓度高、左上方浓度低的一个浓度梯度;扩散中期时,H2S气体呈现设备中部浓度高,设备两边浓度低的浓度分布,最终扩散均匀;
HF气体的扩散规律则与H2S、SO2不同,扩散初期HF气体先在设备边缘向上扩散,形成了左上方浓度高、右下方浓度低的一个浓度梯度,当扩散来到中期,取气口顶部浓度最高,设备本体浓度低,最终扩散均匀;
在相同浓度的扩散情况下,HF气体扩散达到均匀的时间最短;气体浓度越大,扩散达到均匀的时间越长,反之时间越短;当设备内部故障产生的特征组分气体浓度小于等于100μL·L-1时,扩散到取气口浓度低于1μL·L-1,在现有的检测技术条件下,检测不到特征气体,从而造成设备故障漏判。
本发明有益效果:
本发明基于FLUENT仿真得到特征组分气体在SF6电气设备内的扩散规律、扩散达到均匀的时间、浓度及特征组份气体扩散至检测装置能够检测到气体的时间,为SF6电气设备故障判断、故障发生时间预测等;解决了现有技术针对电气设备SF6气体进行分析,由于不了解气体组分在内部的扩散特性,导致设备故障诊断方面尚存在不足,严重阻碍了设备故障诊断技术的发展等技术问题。
具体实施方式
一种SF6电气设备中特征组份气体扩散效应分析方法,它包括:
建立气体绝缘全封闭组合电器室几何模型;
建立几何模型及网格划分
GIS模型等比例建模,该物理模型分为两个部分,其中大圆柱为气体绝缘全封闭组合电器室(GIS),圆柱直径为650(mm),长为2000(mm),上方小圆柱为取气口,直径为10(mm),长为150(mm)。采用Ansys Fluent Meshing工具划分网格,网格数为12843个,网格质量均在0.8以上,满足要求。
步骤2、建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程;
质量守恒方程
Figure BDA0002817872150000081
式中,ρ为气体密度,μj为湍流粘性系数的径向分量,t为时间,xj表示特征组分气体气体沿腔体径向的位移,Sm为源项。
动量守恒方程
X方向:
Figure BDA0002817872150000091
Y方向:
Figure BDA0002817872150000092
Z方向:
Figure BDA0002817872150000093
u、v、w为速度矢量在x、y、z方向的分量,Su、Sv和Sw和是动量方程的广义源项;μ为湍流黏性系数。
能量守恒方程
Figure BDA0002817872150000094
T为流体温度;kt为流体的湍流导热系数;cp为混合气体的定压比热;cpr为SO2定压比热;cpa为SF6的定压比热;ω为重气质量比例;为常数,常取1;μt为流体的湍流粘性系数。
组分输运方程
SF6分解产生特征组分气体,需要利用组分输运方程计算局部每相的质量分数。由组分质量守恒定律可以得到气体扩散组分方程:
Figure BDA0002817872150000095
式中:ω为每一组分的质量分数,Dm为湍流扩散系数。
步骤3、设定初始条件;
气室内压力为0.4Mpa,室内的温度为300K,采用非耦合瞬态的方法进行求解,整个模型是封闭体系即取气口不打开,气体的扩散仅通过浓度差以及分子热运动进行。忽略设备内部母线或其他机构在空间上对扩散过程的影响,定义设备内部特征组份气体最大和最小浓度差小于等于1μL·L-1时,认为特征组份气体扩散达到均匀。研究如下几种不同情况:SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散,每种情况分别模拟六种不同的浓度,最后得出结论.
步骤4、分析SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散情况下的气体扩散效应;
步骤5、对气体扩散效应的结果进行分析,得出结论。
HF扩散模拟
为模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SF6气体。在设备左下位置(0.1,-0.225,0)处由于局部放电故障导致产生了体积分数为5000μL·L-1的HF。使用FLUENT软件模拟HF气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果。
封闭体系下,经过1000s的扩散,设备内部的HF的最大浓度由原来的5000μL·L-1下降到只有44.6μL·L-1,设备内部HF的浓度呈现一个从左上到右下、从大到小的一个浓度梯度,最大浓度在设备的左上方,浓度大小为44.6μL·L-1,最小浓度在设备的底部,浓度大小为24.4μL·L-1。经过3000s的扩散,设备内部的最大浓度继续下降,下降到36.5μL·L-1,最大浓度出现在取气口位置,最小浓度位于设备右下方,此时最大浓度为36.5μL·L-1,最小浓度为32.4μL·L-1,相比于1000s的扩散情况,最大与最小浓度差有明显的减小。经过6420s的扩散,设备内部最大浓度为34.8μL·L-1
最小浓度为33.8μL·L-1,浓度差为1μL·L-1,认为扩散已经达到均匀。
当设备内部HF浓度分别为5000、2000、1000、500、200、100μL·L-1时,其他条件不变,使用FLUENT软件模拟HF气体从开始扩散到扩散均匀所需要的时间如下表1所示。从表中可以发现,封闭体系中,HF气体仅通过浓度差以及分子热运动驱动扩散的情况下,达到扩散均匀所需要的时间分别为:6420s、5890s、5190s、4760s、2580s、200s。
表1不同HF浓度下扩散均匀的时间及最大最小浓度
1 The time required for uniform diffusion and the maximum and minimumconcentration under different HF concentrations
Figure BDA0002817872150000111
SO2扩散模拟
使用FLUENT软件模拟SO2气体从指定位置扩散到整个设备内部的过程,封闭体系下SO2气体从扩散初期到中期以及扩散均匀后的结果,结果所示在封闭体系下,当扩散时间为1000s时,SO2浓度由原来的5000μL·L-1下降到72.2μL·L-1,与HF气体的扩散情况不同,SO2气体呈现一个由下至上、由大到小的一个浓度梯度,最大浓度在设备底部,浓度为72.2μL·L-1,最小浓度出现在取气口处,浓度为0.121μL·L-1。当扩散到达中期,扩散时间为4800s,设备本体内最大浓度已经下降到35.4μL·L-1,设备底部有最大浓度,取气口顶部存在最小浓度,最小浓度为25.2μL·L-1。当扩散时间为9260s时,最大浓度为34.4μL·L-1,最小浓度为33.4μL·L-1,此时扩散达到均匀。
表2不同SO2浓度下扩散均匀的时间及最大最小浓度
2 The time required for uniform diffusion and the maximum and minimumconcentration under different SO2 concentrations
Figure BDA0002817872150000121
H2S扩散模拟
封闭体系下,使用FLUENT软件模拟H2S气体从指定位置扩散到整个设备内部的过程,H2S气体从扩散初期到中期以及到扩散均匀后的结果。
当扩散时间为1000s时,H2S气体的浓度由最开始的5000μL·L-1下降到83.5μL·L-1,此时设备内部呈现的浓度规律为设备右下方浓度最高,最高浓度为83.5μL·L-1,设备左上方浓度最低,最低浓度为0.0557μL·L-1,呈现一个从右下到左上、从大到小的一个浓度梯度。当扩散时间达到4800s时,此时的最大浓度下降到35.8μL·L-1,最小浓度为25.7μL·L-1。此时设备内部的H2S气体浓度呈现一个设备中部浓度高,设备左右侧浓度低的一个规律,浓度最低出现在取气口。当扩散时间为8880s,此时设备内部的最大浓度为34.6μL·L-1,最小浓度为最小浓度为33.6μL·L-1,扩散达到均匀。
当设备内部SO2浓度分别为5000、2000、1000、500、200、100μL·L-1时,其他条件不变,达到扩散均匀所需要的时间分别为8880s、7140s、6850s、6310s、4960s、3430s。
表3不同H2S浓度下扩散均匀的时间及最大最小浓度
3 The time required for uniform diffusion and the maximum and minimumconcentration under different H2S concentrations
Figure BDA0002817872150000131
为了与工程实际相符合,同时进一步研究特征组份气体扩散到取气口处的浓度随时间的变化规律,在取气口位置,坐标为(1,0.475,0)处设置监测点,对于特征组分气体在线监测装置来说,设备精度一般都是1μL·L-1,根据监测点的数据,得到不同浓度、不同特征组份气体在取气口处达到1μL·L-1的时间,结果如下表4。
在封闭体系下,在扩散初期的2500s内,取气口HF气体的浓度有一个快速增长的过程,由0增长到35.5μL·L-1,此时SO2和H2S气体浓度较小。在随后的扩散时间里,取气口HF气体的浓度缓慢下降,SO2和H2S气体的浓度快速上升。当扩散时间达到10000s时,三种气体的浓度变化已经到达了平衡状态,基本上维持在34.5μL·L-1水平,证明设备内部的扩散已经达到均匀。
表4取气口处浓度达到1μL·L-1的时间
Figure BDA0002817872150000141
步骤4、分析SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散情况下的气体扩散效应;
步骤5、对气体扩散效应的结果进行分析,得出结论。
SO2的扩散初期由设备左边缘位置向右扩散,形成设备底部浓度高、顶部浓度低的浓度梯度。扩散达到中期,SO2气体依靠浓度差以及分子热运动由设备底部向顶部扩散,最终扩散均匀。
H2S气体与SO2气体的扩散相似,扩散初期由设备左边缘向右扩散,形成了设备右下方浓度高、左上方浓度低的一个浓度梯度。扩散中期时,H2S气体呈现设备中部浓度高,设备两边浓度低的浓度分布,最终扩散均匀。
HF气体的扩散规律则与H2S、SO2不同,扩散初期HF气体先在设备边缘向上扩散,形成了左上方浓度高、右下方浓度低的一个浓度梯度,当扩散来到中期,取气口顶部浓度最高,设备本体浓度低,最终扩散均匀。
在相同浓度的扩散情况下,HF气体扩散达到均匀的时间最短。气体浓度越大,扩散达到均匀的时间越长,反之时间越短。预测当设备内部故障较轻,产生的特征组分气体浓度较小,可能扩散到取气口浓度低于1μL·L-1,在现有的检测技术条件下,检测不到特征气体,从而造成设备故障漏判。
为了研究气体绝缘全封闭组合电器(GIS)设备内部由于局放或过热等故障产生的不同特征组份气体在设备内部的扩散规律,文章通过CFD技术模拟设备内部不同特征组份气体的扩散情况,针对HF、H2S和SO2三种特征组分气体进行了定性和定量分析,基于FLUENT仿真得出三种不同特征组份气体在GIS设备内部的扩散规律、扩散达到均匀的浓度、时间及特征组份气体扩散至检测装置能够检测到的时间。本发明仅研究了简单模型下的GIS设备中三种气体组分扩散时间、浓度的一定规律,与真实设备存在一定差异,后期可持续研究复杂模型下的GIS设备中特征气体组分的扩散特性。在后续工作中,本发明的研究可为GIS设备故障诊断、故障发生时间预测等研究提供理论依据。

Claims (8)

1.一种SF6电气设备中特征组份气体扩散效应分析方法,它包括:
步骤1、建立气体绝缘全封闭组合电器室几何模型;
步骤2、建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程;
步骤3、设定初始条件;
步骤4、分析SO2气体单独扩散、HF气体单独扩散以及H2S气体单独扩散情况下的气体扩散效应;
步骤5、对气体扩散效应的结果进行分析,得出结论。
2.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:步骤1所述建立气体绝缘全封闭组合电器室几何模型为:模型分为两个部分,其中大圆柱为气体绝缘全封闭组合电器室GIS,圆柱直径为650mm,长为2000mm,上方有小圆柱取气口,直径为10mm,长为150mm;采用Ansys Fluent Meshing工具划分网格,网格数为12843个,网格质量均在0.8以上。
3.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:步骤2所述建立数学模型;所述数学模型包括质量守恒方程、动量守恒方程、能量守恒方程和组分运输方程的公式为:
质量守恒方程:
Figure FDA0002817872140000011
式中,ρ为气体密度,μj为湍流粘性系数的径向分量,t为时间,xj表示特征组分气体气体沿腔体径向的位移,Sm为源项;
动量守恒方程:
X方向:
Figure FDA0002817872140000012
Y方向:
Figure FDA0002817872140000021
Z方向:
Figure FDA0002817872140000022
u、v、w为速度矢量在x、y、z方向的分量,Su、Sv和Sw和是动量方程的广义源项;μ为湍流黏性系数;
能量守恒方程:
Figure FDA0002817872140000023
T为流体温度;kt为流体的湍流导热系数;cp为混合气体的定压比热;cpr为SO2定压比热;cpa为SF6的定压比热;ω为重气质量比例;为常数,常取1;μt为流体的湍流粘性系数;
组分输运方程:
SF6分解产生特征组分气体,需要利用组分输运方程计算局部每相的质量分数。由组分质量守恒定律可以得到气体扩散组分方程:
Figure FDA0002817872140000024
式中:ω为每一组分的质量分数,Dm为湍流扩散系数。
4.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:所述初始条件包括:气室内压力为0.4Mpa,室内的温度为300K,采用非耦合瞬态的方法进行求解,整个模型是封闭体系即取气口不打开,气体的扩散通过浓度差以及分子热运动进行;忽略设备内部母线或机构在空间上对扩散过程的影响,定义设备内部特征组份气体最大和最小浓度差小于等于1μL·L-1时,认为特征组份气体扩散达到均匀。
5.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:所述分析HF气体单独扩散情况下气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SF6气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的HF,使用FLUENT软件模拟HF气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,经过1000s的扩散,设备内部的HF的最大浓度由原来的5000μL·L-1下降到44.6μL·L-1,设备内部HF的浓度呈现一个从左上到右下、从大到小的一个浓度梯度,最大浓度在设备的左上方,浓度大小为44.6μL·L-1,最小浓度在设备的底部,浓度大小为24.4μL·L-1;经过3000s的扩散,设备内部的最大浓度继续下降,下降到36.5μL·L-1,最大浓度出现在取气口位置,最小浓度位于设备右下方,此时最大浓度为36.5μL·L-1,最小浓度为32.4μL·L-1,经过6420s的扩散,设备内部最大浓度为34.8μL·L-1,最小浓度为33.8μL·L-1,浓度差为1μL·L-1,认为扩散已经达到均匀;
当设备内部HF初始浓度分别为5000、2000、1000、500、200和100μL·L-1时,使用FLUENT软件模拟HF气体从开始扩散到扩散均匀所需要的时间,如下表1所示;得出在封闭体系中,HF气体仅通过浓度差以及分子热运动驱动扩散的情况下,达到扩散均匀所需要的时间分别为:6420s、5890s、5190s、4760s、2580s、200s;
表1 不同HF浓度下扩散均匀的时间及最大最小浓度
Figure FDA0002817872140000031
Figure FDA0002817872140000041
6.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:所述SO2气体单独扩散情况下气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SO2气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的SO2,使用FLUENT软件模拟HF气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,当扩散时间为1000s时,SO2浓度由原来的5000μL·L-1下降到72.2μL·L-1,SO2气体呈现一个由下至上、由大到小的一个浓度梯度,最大浓度在设备底部,浓度为72.2μL·L-1,最小浓度出现在取气口处,浓度为0.121μL·L-1;当扩散到达中期,扩散时间为4800s,设备本体内最大浓度已经下降到35.4μL·L-1,设备底部有最大浓度,取气口顶部存在最小浓度,最小浓度为25.2μL·L-1;当扩散时间为9260s时,最大浓度为34.4μL·L-1,最小浓度为33.4μL·L-1,此时扩散达到均匀;
当设备内部SO2浓度分别为5000、2000、1000、500、200和100μL·L-1时,达到扩散均匀所需要的时间分别为9260s、7310s、6130s、4850s、3560s、2760s;
表2 不同SO2浓度下扩散均匀的时间及最大最小浓度
Figure FDA0002817872140000042
Figure FDA0002817872140000051
7.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:分析H2S气体单独扩散情况下的气体扩散效应的方法为:模拟设备内部出现局部放电故障产生特征组分气体的扩散情况,假设设备本体内充满SF6气体,在设备左下位置处由于局部放电故障导致产生了体积分数为5000μL·L-1的H2S,使用FLUENT软件模拟H2S气体从指定位置扩散到整个设备内部的过程,封闭体系下从扩散初期到中期最后扩散均匀后的结果;它具体包括:
在封闭体系下,使用FLUENT软件模拟H2S气体从指定位置扩散到整个设备内部的过程,H2S气体从扩散初期到中期以及到扩散均匀后的结果;
当扩散时间为1000s时,H2S气体的浓度由最开始的5000μL·L-1下降到83.5μL·L-1,此时设备内部呈现的浓度规律为设备右下方浓度最高,最高浓度为83.5μL·L-1,设备左上方浓度最低,最低浓度为0.0557μL·L-1,呈现一个从右下到左上、从大到小的一个浓度梯度;当扩散时间达到4800s时,此时的最大浓度下降到35.8μL·L-1,最小浓度为25.7μL·L-1。此时设备内部的H2S气体浓度呈现一个设备中部浓度高,设备左右侧浓度低的一个规律,浓度最低出现在取气口;当扩散时间为8880s,此时设备内部的最大浓度为34.6μL·L-1,最小浓度为最小浓度为33.6μL·L-1,扩散达到均匀;
当设备内部SO2浓度分别为5000、2000、1000、500、200、100μL·L-1时,其他条件不变,达到扩散均匀所需要的时间分别为8880s、7140s、6850s、6310s、4960s、3430s;
表3 不同H2S浓度下扩散均匀的时间及最大最小浓度
Figure FDA0002817872140000061
8.根据权利要求1所述的一种SF6电气设备中特征组份气体扩散效应分析方法,其特征在于:步骤5所述对气体扩散效应的结果进行分析,得出结论的方法包括:
SO2的扩散初期由设备左边缘位置向右扩散,形成设备底部浓度高、顶部浓度低的浓度梯度。扩散达到中期,SO2气体依靠浓度差以及分子热运动由设备底部向顶部扩散,最终扩散均匀;
H2S气体与SO2气体的扩散相似,扩散初期由设备左边缘向右扩散,形成了设备右下方浓度高、左上方浓度低的一个浓度梯度;扩散中期时,H2S气体呈现设备中部浓度高,设备两边浓度低的浓度分布,最终扩散均匀;
HF气体的扩散规律则与H2S、SO2不同,扩散初期HF气体先在设备边缘向上扩散,形成了左上方浓度高、右下方浓度低的一个浓度梯度,当扩散来到中期,取气口顶部浓度最高,设备本体浓度低,最终扩散均匀;
在相同浓度的扩散情况下,HF气体扩散达到均匀的时间最短;气体浓度越大,扩散达到均匀的时间越长,反之时间越短;当设备内部故障产生的特征组分气体浓度小于等于100μL·L-1时,扩散到取气口浓度低于1μL·L-1,在现有的检测技术条件下,检测不到特征气体,从而造成设备故障漏判。
CN202011403631.XA 2020-12-04 2020-12-04 一种sf6电气设备中特征组份气体扩散效应分析方法 Active CN112464482B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011403631.XA CN112464482B (zh) 2020-12-04 2020-12-04 一种sf6电气设备中特征组份气体扩散效应分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011403631.XA CN112464482B (zh) 2020-12-04 2020-12-04 一种sf6电气设备中特征组份气体扩散效应分析方法

Publications (2)

Publication Number Publication Date
CN112464482A true CN112464482A (zh) 2021-03-09
CN112464482B CN112464482B (zh) 2023-05-12

Family

ID=74805389

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011403631.XA Active CN112464482B (zh) 2020-12-04 2020-12-04 一种sf6电气设备中特征组份气体扩散效应分析方法

Country Status (1)

Country Link
CN (1) CN112464482B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640633A (zh) * 2021-08-12 2021-11-12 贵州大学 一种气体绝缘组合电器设备故障定位方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101059485A (zh) * 2007-05-23 2007-10-24 重庆大学 六氟化硫放电分解气体组分分析系统及其使用方法
CN103884646A (zh) * 2014-03-04 2014-06-25 中国空间技术研究院 一种基于光纤传感的喷口电弧烧蚀金属蒸汽浓度测量系统
CN109580431A (zh) * 2019-01-28 2019-04-05 贵州电网有限责任公司 一种sf6/n2混合气体扩散特性确定方法
CN110793888A (zh) * 2019-11-20 2020-02-14 国网电力科学研究院武汉南瑞有限责任公司 Sf6/n2混合气体扩散特性实验方法
WO2020133401A1 (zh) * 2018-12-27 2020-07-02 中国科学院光电技术研究所 超大面积扫描式反应离子刻蚀机及刻蚀方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101059485A (zh) * 2007-05-23 2007-10-24 重庆大学 六氟化硫放电分解气体组分分析系统及其使用方法
CN103884646A (zh) * 2014-03-04 2014-06-25 中国空间技术研究院 一种基于光纤传感的喷口电弧烧蚀金属蒸汽浓度测量系统
WO2020133401A1 (zh) * 2018-12-27 2020-07-02 中国科学院光电技术研究所 超大面积扫描式反应离子刻蚀机及刻蚀方法
CN109580431A (zh) * 2019-01-28 2019-04-05 贵州电网有限责任公司 一种sf6/n2混合气体扩散特性确定方法
CN110793888A (zh) * 2019-11-20 2020-02-14 国网电力科学研究院武汉南瑞有限责任公司 Sf6/n2混合气体扩散特性实验方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
DACHANG CHEN 等: "High selectivity n-type InSe monolayer toward decomposition products of sulfur hexafluoride:A density functional theory study", 《APPLIED SURFACE SCIENCE》 *
丁慧心等: "地下水电站出线竖井中SF_6泄露扩散的数值模拟分析", 《制冷与空调(四川)》 *
王东东等: "FLUENT在公路隧道有毒气体事故泄漏扩散研究中的应用", 《安全与环境学报》 *
钱进 等: "SF6电气设备中分解物SO2的在线监测技术及应用", 《水电能源科学》 *
陈邦国等: "GIL廊道内SF_6气体泄漏扩散及风机安装优化研究", 《水电能源科学》 *
马尚权等: "基于FLUENT的SF_6气体扩散数值模拟", 《内蒙古煤炭经济》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640633A (zh) * 2021-08-12 2021-11-12 贵州大学 一种气体绝缘组合电器设备故障定位方法
CN113640633B (zh) * 2021-08-12 2024-04-09 贵州大学 一种气体绝缘组合电器设备故障定位方法

Also Published As

Publication number Publication date
CN112464482B (zh) 2023-05-12

Similar Documents

Publication Publication Date Title
CN102654446B (zh) 一种高温高压多相流腐蚀实验方法及装置
CN103412014B (zh) 一种氟化氢气体在线检测仪及方法
David et al. Dynamic modelling of alkaline self-pressurized electrolyzers: a phenomenological-based semiphysical approach
CN112464482A (zh) 一种sf6电气设备中特征组份气体扩散效应分析方法
CN106198655A (zh) 一种电化学传感器检测sf6分解气体的矩阵计算方法
Han et al. Coordinated analysis of multiple factors of argon blowing parameters on the effect of circulation flow rate in RH vacuum refining process
CN110648724B (zh) 一种电渗析传质过程的数值模拟方法
Pao et al. Numerical validation of gas-liquid slug flow inside horizontal pipe
Gevorkov et al. PLC-based hardware-in-the-loop simulator of a centrifugal pump
CN108090274A (zh) 一种飞机油箱机载惰化系统的工程设计方法
CN107703215A (zh) 一种gis设备放电故障化学诊断方法
CN112924332B (zh) 基于cfd技术的gis设备中so2扩散特性确定方法
CN108052706A (zh) 空气放电组分在空气开关柜下的气体扩算仿真方法及系统
CN114415004A (zh) 一种模拟真型gis设备故障特征气体扩散方法
CN211347828U (zh) 一种sf6混合绝缘气体试验装置
CN105424819B (zh) 混合气水合反应过程中的离子浓度多层位监测装置及监测方法
CN104111301B (zh) 一种油色谱在线监测数据与离线数据智能校验方法
CN113720973A (zh) Gis设备内部sf6气体监测系统以及监测方法
Riedmann et al. Gas Losses in Transformers-Influences and Consideration
Heiselberg et al. Buoyancy driven natural ventilation through horizontal openings
Matas et al. Numerical simulation of flow and determination of aerodynamic forces in the balanced control valve
Chen et al. The effect of quantity of salt on the drying characteristics of fresh noodles
Bagul et al. Computational study of behavior of gas absorption in data center equipment and its effects on the rate of corrosion/contamination
Choi et al. Analysis of Diffusion and Convection Characteristics of Humid Air in a Transportation Pipe
Pizano-Moreno et al. Leak isolation with temperature compensation in pipelines

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant