CN103076548B - 一种用表面电导率和泄漏电流预测闪络电压的方法 - Google Patents
一种用表面电导率和泄漏电流预测闪络电压的方法 Download PDFInfo
- Publication number
- CN103076548B CN103076548B CN201310040885.3A CN201310040885A CN103076548B CN 103076548 B CN103076548 B CN 103076548B CN 201310040885 A CN201310040885 A CN 201310040885A CN 103076548 B CN103076548 B CN 103076548B
- Authority
- CN
- China
- Prior art keywords
- arc
- insulator
- leakage current
- voltage
- flashover
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 12
- 239000012212 insulator Substances 0.000 claims abstract description 63
- 238000011161 development Methods 0.000 claims abstract description 10
- 238000009422 external insulation Methods 0.000 claims abstract description 9
- 238000009413 insulation Methods 0.000 claims abstract description 6
- 238000010891 electric arc Methods 0.000 claims description 22
- 230000008033 biological extinction Effects 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 8
- 239000000725 suspension Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 6
- 230000005611 electricity Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 3
- 239000011521 glass Substances 0.000 description 2
- 229910052573 porcelain Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Testing Relating To Insulation (AREA)
Abstract
本发明公开了一种用表面电导率和泄漏电流预测闪络电压的方法,测量绝缘子表面电导率,以绝缘子表面电导率作为输入参数,结合绝缘子形状系数参数,建立动态电弧模型,计算泄漏电流发展趋势,得到绝缘子串的最低闪络电压,基于最低闪络电压判断该绝缘子的绝缘裕度及外绝缘水平。本发明根据实测绝缘子串整体电导率,利用动态电弧模型计算泄漏电流的发展趋势,电弧动态模型通过数学方程来描述泄漏电流以及电弧长度随时间的变化规律,更符合真实情况,进而实现绝缘状态判断与污闪预警的目的。
Description
技术领域
本发明属于输电线路外绝缘状态评估领域,以绝缘子表面污层电导率作为输入参数,可仿真得到泄漏电流发展趋势,并预测绝缘子的闪络电压,进而评估输电线路的外绝缘水平,适用于输电线路广泛使用的瓷绝缘子和玻璃绝缘子。
背景技术
污闪是可能导致电力系统出现灾难性事故的重要原因之一。为保证电网的安全运行需要在污闪事故发生前能做出预警。闪络电压是表征绝缘子外绝缘性能的重要参数,通过测量、估算绝缘子串的闪络电压是评价绝缘子串外绝缘性能的有效方式。
实测绝缘子串的污闪电压存在一定的难度。更重要的是在线预警技术需要在闪络发生前即可通过测量绝缘子的电气参数预估其闪络电压。
静态模型对临界污闪状态的分析有很大的作用。但是,对于描述绝缘子电压、电流、电弧的动态变化过程存在明显的不足。
发明内容
本发明提供了一种用电导率和泄漏电流发展趋势预测闪络电压的方法,采用电弧动态模型通过数学方程来描述绝缘子泄漏电流以及电弧长度随时间的变化规律,更符合真实情况。
以绝缘子表面污层电导率作为输入参数,计算得到泄漏电流发展趋势,预测绝缘子的闪络电压,进而评估输电线路的外绝缘水平。
一种用表面电导率和泄漏电流预测闪络电压的方法,其特征在于,包括以下步骤:
(1)测量绝缘子表面电导率,
(2)以绝缘子表面电导率作为输入参数,结合绝缘子形状系数参数,建立动态电弧模型,
(3)计算泄漏电流发展趋势,
(4)得到绝缘子串的最低闪络电压,基于最低闪络电压判断该绝缘子的绝缘裕度及外绝缘水平。
建立动态电弧模型的步骤如下:
设定绝缘子形状系数,给定参数的初始值;
输入测得的绝缘子表面电导率σ;
假设在零时刻,电压u(t)达到波峰,燃弧开始:
u(t)=Um cosωt (1)
其中,Um为电压有效值,ω为角频率,t为时间;
平板型悬式绝缘子剩余污层电阻随着弧根半径和电弧长度的变化公式:
其中,L为泄漏距离总长度,Larc为当前条件下的电弧长度;σe为表面污层电导率;r0为弧根半径,由以下公式求得:
其中,该状态下对应的泄漏电流Iarc由公式(4)求得:
其中,u(t)为此时电压,Rp为绝缘子剩余污层电阻,rarc为单位长度的电弧电阻,Larc为当前条件下的电弧长度;
由以下公式计算电弧场强Earc和剩余污层场强Ep:
Earc=63Iarc -0.5 (5)
电弧延伸条件为:
Earc<Ep (7)
当Ep<Earc时,随着电压的变化,电弧长度变化的表达式如式(8)所示:
当Ep>Earc时,电弧长度变化的表达式如式(9)所示:
电弧长度变化还有附加的延伸参数μ;
建立单位长度电弧电阻变化的动态方程:
其中,N0为单位长电弧的耗散功率;τ为时间常数,由下式计算:
其中,Q0为试验所得耗散功率常数;
定义电弧的熄灭条件为:
Rarc>Rn (12)
其中Rarc为电弧总电阻,Rarc=rarcLarc;Rn为一个设定的电阻值,若此时由于电弧电阻远大于干区的电阻,电流通过干区后再流过剩余污层,电弧熄灭;
电弧熄灭后,说明在此时的表面电导率σ下,Um不足以使得绝缘子发生闪络,故增大Um,直至得到一个最小的Um,使得绝缘子发生闪络,此值即为得到的该绝缘子串的最低闪络电压。
上述公式中参数取值为:N0=60,Q0=0.075,μ=5。
本发明所达到的有益效果:
本发明的用电导率和泄漏电流发展趋势预测闪络电压的方法,以绝缘子表面电导率作为输入参数,通过输入测量所得绝缘子表面电导率,结合系统内置绝缘子形状系数参数,通过动态电弧模型计算泄漏电流发展趋势,进而最终得到绝缘子串的闪络电压,进而预测评估输电线路的外绝缘水平,实现绝缘状态判断与污闪预警的目的。本发明的方法采用电弧动态模型通过数学方程来描述绝缘子泄漏电流以及电弧长度随时间的变化规律,更符合真实情况,适用于输电线路广泛使用的瓷绝缘子和玻璃绝缘子。
附图说明
图1为计算流程图;
图2(a)为未发生闪络时的动态电弧模型仿真泄漏电流波形;泄漏电流波形对比图;
图2(b)为未发生闪络时的实测泄漏电流波形;
图3(a)为发生闪络时的动态电弧模型仿真泄漏电流波形;泄漏电流波形对比图;
图3(b)为发生闪络时的实测泄漏电流波形
图4为污层整体电导率和污闪电压关系图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,计算开始时,设定绝缘子形状系数。
输入试验测得的绝缘子表面电导率σ,给定各个参数的初始值。
假设在零时刻,电压u(t)达到波峰,燃弧开始:
u(t)=Um cosωt (1)
其中,Um为电压有效值,ω为角频率,t为时间。
平板型悬式绝缘子剩余污层电阻Rp随着弧根半径r0和电弧长度Larc的变化公式:
其中,r0为弧根半径;L为泄漏距离总长度,Larc为当前条件下的电弧长度;σe为表面污层电导率;
r0为弧根半径,由以下公式求得:
该状态下对应的泄漏电流Iarc(电弧电流)由公式(4)求得,表达式如下:
其中,u(t)为此时电压,Rp为绝缘子剩余污层电阻,rarc为单位长度的电弧电阻,Larc为当前条件下的电弧长度;
Earc和Ep分别为电弧场强和剩余污层场强。由以下公式计算:
Earc=63Iarc -0.5 (5)
电弧延伸条件为:
Earc<Ep (7)
当Ep<Earc时,随着电压的变化,电弧长度变化的表达式如下:
当Ep>Earc时,电弧长度变化还有附加的延伸参数μ:
单位长度电弧电阻变化的动态方程,如下式所示:
其中,N0为单位长电弧的耗散功率,τ为时间常数,由下式计算:
其中,Q0为试验所得耗散功率常数;
如果电弧的输入能量始终小于散失能量,电弧电阻将逐步增加,最终电弧逐渐熄灭。因此,定义电弧的熄灭条件为:
Rarc>Rn (12)
其中Rarc为电弧总电阻,Rarc=rarcLarc;Rn代表描述电弧熄灭的电阻,远大于常态下的电弧电阻,仿真中可取1MΩ(Rn取值只要足够大,就不会影响电弧发展时的电流和弧长变化趋势。同时该值可根据仿真结果加以优化,以提高计算效率)。当Rarc>Rn时,由于电弧电阻远大于干区的电阻,因此电流通过干区后再流过剩余污层,电弧熄灭。
电弧熄灭后,说明在此时的表面电导率σ下,Um不足以使得绝缘子发生闪络,故增大Um,直至得到一个最小的Um,使得绝缘子发生闪络。此值即为该串绝缘子的闪络电压。
仿真时的电弧参数和实测参数有很大的不同,其中,影响仿真结果的最重要的是污层的耗散功率N0,该参数决定了电弧能量的增减,直接影响电弧电阻的动态变化,进而影响绝缘子的污闪电压和泄漏电流变化规律。因此,本模型中使用的参数较佳地取值为:
N0=60,Q0=0.075,μ=5。
上述参数都是比例系数,单位取1。结果表明模型仿真结果和人工污秽试验结果符合度较好。
对上述建立的动态电弧模型进行仿真,图2(a)、图2(b)和图3(a)、图3(b)分别为未发生污闪和发生污闪时候的泄漏电流仿真波形和实测波形的对比。仿真和实测的对象均为XP-70型绝缘子。由于电弧没有明显的熄灭、重燃现象,无论绝缘子是否发生污闪,模拟电流和实测电流波形非常相似。
在未发生污闪时,泄漏电流先达到最大值,然后在几个周波内迅速减小。发生污闪时,在临闪前,泄漏电流均达到了安培级,泄漏电流逐步增大,最终电弧贯穿绝缘子而形成闪络。仿真波形和实测波形一个显著的差别是:在实测波形中,尽管电弧过零不熄灭,电流波形也有可能出现“零休”现象。例如实测未闪络波形电弧熄灭前的最后一个周波以及闪络波形的前一个周波。
为了进一步验证模型的准确性,对X-4.5型悬式绝缘子在相同污层整体电导率下的污闪电压的理论计算结果和试验结果进行了比较,比较结果如图4所示。可以发现,试验结果和仿真结果的差异很小,证明了该染污放电模型对污闪电压的预测符合实际情况。
本方法在剩余污层电阻、电弧发展速度、电弧熄灭条件、仿真参数等方面建立了与实测波形相似、试验结果吻合的染污放电泄漏电流仿真模型。试验结果和模型仿真结果得到了相互验证。
通过本方法预测得到的最低闪络电压与真实绝缘子闪络电压契合度较高,对于实测XP-70绝缘子,预测值和真实值之间差异<5%。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (2)
1.一种用表面电导率和泄漏电流预测闪络电压的方法,包括以下步骤:
(1)测量绝缘子表面电导率,
(2)以绝缘子表面电导率作为输入参数,结合绝缘子形状系数参数,建立动态电弧模型,
(3)计算泄漏电流发展趋势,
(4)得到绝缘子串的最低闪络电压,基于最低闪络电压判断该绝缘子的绝缘裕度及外绝缘水平;
其特征在于,
建立动态电弧模型的步骤如下:
设定绝缘子形状系数,给定参数的初始值;
输入测得的绝缘子表面电导率σ;
假设在零时刻,电压u(t)达到波峰,燃弧开始:
u(t)=Um cosωt (1)
其中,Um为电压有效值,ω为角频率,t为时间;
平板型悬式绝缘子剩余污层电阻随着弧根半径和电弧长度的变化公式:
其中,L为泄漏距离总长度,Larc为当前条件下的电弧长度;σe为表面污层电导率;r0为弧根半径,由以下公式求得:
其中,该状态下对应的泄漏电流Iarc由公式(4)求得:
其中,u(t)为此时电压,Rp为绝缘子剩余污层电阻,rarc为单位长度的电弧电阻,Larc为当前条件下的电弧长度;
由以下公式计算电弧场强Earc和剩余污层场强Ep:
Earc=63Iarc -0.5 (5)
电弧延伸条件为:
Earc<Ep (7)
当Ep<Earc时,随着电压的变化,电弧长度变化的表达式如式(8)所示:
当Ep>Earc时,电弧长度变化的表达式如式(9)所示:
电弧长度变化还有附加的延伸参数μ;
建立单位长度电弧电阻变化的动态方程:
其中,N0为单位长电弧的耗散功率;τ为时间常数,由下式计算:
其中,Q0为试验所得耗散功率常数;
定义电弧的熄灭条件为:
Rarc>Rn (12)
其中Rarc为电弧总电阻,Rarc=rarcLarc;Rn为一个设定的电阻值,若此时由于电弧电阻远大于干区的电阻,电流通过干区后再流过剩余污层,电弧熄灭;
电弧熄灭后,说明在此时的表面电导率σ下,Um不足以使得绝缘子发生闪络,故增大Um,直至得到一个最小的Um,使得绝缘子发生闪络,此值即为得到的该绝缘子串的最低闪络电压。
2.根据权利要求1所述的用表面电导率和泄漏电流预测闪络电压的方法,其特征在于,
上述公式中参数取值为:N0=60,Q0=0.075,μ=5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310040885.3A CN103076548B (zh) | 2013-02-01 | 2013-02-01 | 一种用表面电导率和泄漏电流预测闪络电压的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310040885.3A CN103076548B (zh) | 2013-02-01 | 2013-02-01 | 一种用表面电导率和泄漏电流预测闪络电压的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103076548A CN103076548A (zh) | 2013-05-01 |
CN103076548B true CN103076548B (zh) | 2015-04-29 |
Family
ID=48153141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310040885.3A Active CN103076548B (zh) | 2013-02-01 | 2013-02-01 | 一种用表面电导率和泄漏电流预测闪络电压的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103076548B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137299A (zh) * | 2015-08-26 | 2015-12-09 | 芜湖市凯鑫避雷器有限责任公司 | 绝缘子闪络电压预测方法 |
CN105785235B (zh) * | 2016-02-29 | 2020-02-18 | 华南理工大学 | 一种评定绝缘子交流闪络电压特性的试验方法 |
CN105868544B (zh) * | 2016-03-25 | 2018-03-02 | 华北电力大学(保定) | 三参数韦伯分布处理闪络电压分析绝缘材料性能的方法 |
CN105810373A (zh) * | 2016-04-29 | 2016-07-27 | 南方电网科学研究院有限责任公司 | 一种断路器灭弧室瓷套及其表面电导率的计算方法 |
CN105719779A (zh) * | 2016-04-29 | 2016-06-29 | 南方电网科学研究院有限责任公司 | 一种支柱复合绝缘子及其表面电导率的计算方法 |
CN109799442B (zh) * | 2019-03-29 | 2021-11-19 | 云南电网有限责任公司电力科学研究院 | 基于机载高光谱的绝缘子污闪预测方法及系统 |
CN111159866B (zh) * | 2019-12-18 | 2023-05-30 | 东南大学 | 一种基于随机漫步理论的复合绝缘子电弧预测方法 |
CN112858851B (zh) * | 2021-01-12 | 2022-11-08 | 国网江苏省电力有限公司电力科学研究院 | 一种适用于特殊工业粉尘地区的复合绝缘子运行性能状态评价方法 |
CN115144790B (zh) * | 2022-09-02 | 2022-11-18 | 山东理工大学 | 一种高压绝缘子泄漏电流在线监测系统及方法 |
CN115203981B (zh) * | 2022-09-14 | 2022-12-06 | 山东理工大学 | 基于严重覆雪绝缘子单电弧发展闪络的等效电路方法 |
CN115877154B (zh) * | 2023-02-27 | 2023-08-15 | 国网江西省电力有限公司电力科学研究院 | 获取积污支柱绝缘子鸟巢材料短接混合放电特性的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102511065B (zh) * | 2011-10-08 | 2013-07-17 | 清华大学深圳研究生院 | 绝缘子及输电线设备 |
CN102645618B (zh) * | 2012-04-16 | 2014-11-12 | 浙江清科电力科技有限公司 | 一种污秽绝缘智能检测方法 |
-
2013
- 2013-02-01 CN CN201310040885.3A patent/CN103076548B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103076548A (zh) | 2013-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103076548B (zh) | 一种用表面电导率和泄漏电流预测闪络电压的方法 | |
Kong et al. | Fault diagnosis and quantitative analysis of micro-short circuits for lithium-ion batteries in battery packs | |
Xie et al. | An improved resistance-based thermal model for a pouch lithium-ion battery considering heat generation of posts | |
CN109307827B (zh) | 一种混合气体介质条件下gil内绝缘状态判断方法 | |
CN103605042B (zh) | 基于自适应粒子群算法的接地网故障诊断方法 | |
CN103400019B (zh) | 基于Matlab/Simulink的交流故障电弧仿真方法 | |
CN106443381A (zh) | 动态电弧模型构建方法与系统 | |
CN107561397B (zh) | Gis触头系统的电接触状态检测方法和系统 | |
CN104217061A (zh) | 低压配电柜的温度场仿真设计方法 | |
CN104361196A (zh) | 一种计及列车速度的弓网离线电弧数学模型 | |
CN103645014A (zh) | 用于gis设备的sf6气体泄漏率检测方法 | |
CN104142455A (zh) | 基于特勒根定理的接地网腐蚀故障诊断方法 | |
CN105205229A (zh) | 一种高压单芯电缆应急负荷下的载流能力动态评估方法 | |
CN110633496B (zh) | 一种基于热-力耦合模型的锂离子电池放电过程中热应力和温度的确定方法 | |
CN105243232B (zh) | 将场分析集成到电网络求解的电磁暂态仿真方法及系统 | |
CN104008288A (zh) | 一种变压器寿命仿真估计方法 | |
CN110543717A (zh) | 一种gis设备温升特性仿真模型的构建方法及装置 | |
Vita et al. | Comparison of ANN and finite element analysis simulation software for the calculation of the electric field around metal oxide surge arresters | |
CN104505813B (zh) | 一种在贝杰龙线路模型下利用瞬时功率进行测后模拟的输电线路纵联保护方法 | |
CN105021871B (zh) | 一种不完全信息下电缆线路泄露电流确定方法 | |
CN110781589B (zh) | 检测气体绝缘金属封闭输电线路搭接处过热故障的方法 | |
CN107122538A (zh) | 一种特高压直流换流站接头端子的温升优化方法及系统 | |
CN103412190A (zh) | 基于参数在线辨识的开关类设备状态评价方法 | |
CN104237673A (zh) | 一种变压器谐波损耗的估算方法 | |
CN103777102A (zh) | 考虑三次电压谐波影响的避雷器绝缘性能检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |