CN113625067B - 一种离子流场中悬浮液滴荷电特性测量装置及方法 - Google Patents

一种离子流场中悬浮液滴荷电特性测量装置及方法 Download PDF

Info

Publication number
CN113625067B
CN113625067B CN202110922336.3A CN202110922336A CN113625067B CN 113625067 B CN113625067 B CN 113625067B CN 202110922336 A CN202110922336 A CN 202110922336A CN 113625067 B CN113625067 B CN 113625067B
Authority
CN
China
Prior art keywords
charge density
value
charge
electric field
suspended liquid
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
Application number
CN202110922336.3A
Other languages
English (en)
Other versions
CN113625067A (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.)
Huaneng Clean Energy Research Institute
North China Electric Power University
Huaneng Group Technology Innovation Center Co Ltd
Original Assignee
Huaneng Clean Energy Research Institute
North China Electric Power University
Huaneng Group Technology Innovation Center 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 Huaneng Clean Energy Research Institute, North China Electric Power University, Huaneng Group Technology Innovation Center Co Ltd filed Critical Huaneng Clean Energy Research Institute
Priority to CN202110922336.3A priority Critical patent/CN113625067B/zh
Publication of CN113625067A publication Critical patent/CN113625067A/zh
Application granted granted Critical
Publication of CN113625067B publication Critical patent/CN113625067B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/24Arrangements for measuring quantities of charge

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明涉及一种离子流场中悬浮液滴荷电特性测量装置及方法。该方法包括:确定计算区域并对计算区域进行剖分;给定各节点电荷密度初值和荷电因子初值;计算悬浮液滴空间电荷密度和各节点总空间电荷密度;使用有限元法求解泊松方程,得到节点电位分布,计算空间电场;使用上流元法求解电流连续性方程,计算各节点空间电荷密度;若满足误差限定条件,则计算误差函数值,若不满足误差限定条件,则更新表面电荷密度初值,重复上述步骤,直到满足误差限定条件为止。本发明能够准确计算得到离子流场中悬浮液滴的荷电因子。

Description

一种离子流场中悬浮液滴荷电特性测量装置及方法
技术领域
本发明涉及离子流场中悬浮液滴荷电特性测量领域,特别是涉及一种离子流场中悬浮液滴荷电特性测量装置及方法。
背景技术
为了解决中国能源分布不均的问题,中国已经建成了多条特高压直流输电线路。特高压直流输电线路主要用于远距离输电,跨越区域气候环境复杂,周围环境湿度变化较大。特高压直流输电线路在运行过程中,导线会发生电晕放电现象,电晕放电产生的空间电荷会形成离子流,形成了特有的离子流场问题。随着周围环境温湿度的变化,尤其是在较高湿度环境下,空气中会出现悬浮液滴。悬浮液滴会吸附导线电晕放电产生的空间电荷,形成荷电悬浮液滴,导致周围电场发生畸变,使得离子流场问题变得更加复杂。因此2019年11月中国对GB50790-2013“±800kV直流架空输电线路设计规范”进行了局部修改,修订版进一步明确直流线路设计必须考虑湿度对离子流场的影响。
为了明确湿度对高压直流输电线路离子流场的影响规律,需要计算不同湿度条件下悬浮液滴荷电特性。但是离子流场中悬浮液滴的粒径、数密度等参数难以直接测量,无法使用荷电模型直接计算得到悬浮液滴的荷电特性,导致当前离子流场中悬浮液滴荷电特性尚不明确。
发明内容
本发明的目的是提供一种离子流场中悬浮液滴荷电特性测量装置及方法,以解决当前离子流场中悬浮液滴荷电特性测量准确性差的问题。
为实现上述目的,本发明提供了如下方案:
一种离子流场中悬浮液滴荷电特性测量装置,包括:接地电晕笼、中心电极以及场磨;
所述中心电极横穿所述接地电晕笼,所述接地电晕笼为三段式结构,所述接地电晕笼的中间段为测量段,所述接地电晕笼的两侧为屏蔽段,所述测量段内设有离子流板;所述测量段用于测量高压直流导线离子流场,所述屏蔽段用于削弱端部效应;所述场磨设于所述测量段的底部,所述场磨用于测量接地电晕笼笼壁处的电场强度;利用高压直流电源对所述中心电极施加电压,使得所述中心电极处于电晕放电状态,根据所述接地电晕笼笼壁处的电场强度测量离子流场中悬浮液滴荷电特性。
可选的,所述接地电晕笼为同轴圆柱结构,所述中心电极为裸铜导线,所述中心电极横穿所述接地电晕笼的中心位置。
可选的,还包括:支架以及绝缘子;
所述支架设于两侧所述屏蔽段的底部,且所述场磨设于所述测量段与所述支架之间;所述支架用于支撑所述接地电晕笼;
所述绝缘子设于所述中心电极的两端,所述绝缘子用于支撑所述中心电极。
一种离子流场中悬浮液滴荷电特性测量方法,所述离子流场中悬浮液滴荷电特性测量方法应用于一种离子流场中悬浮液滴荷电特性测量装置,包括:接地电晕笼、中心电极以及场磨;所述中心电极横穿所述接地电晕笼,所述接地电晕笼为三段式结构,所述接地电晕笼的中间段为测量段,所述接地电晕笼的两侧为屏蔽段,所述测量段内设有离子流板;所述测量段用于测量高压直流导线离子流场,所述屏蔽段用于削弱端部效应;所述场磨设于所述测量段的底部,所述场磨用于测量接地电晕笼笼壁处的电场强度;利用高压直流电源对所述中心电极施加电压,使得所述中心电极处于电晕放电状态,根据所述接地电晕笼笼壁处的电场强度测量离子流场中悬浮液滴荷电特性;
所述离子流场中悬浮液滴荷电特性测量方法包括:
设定离子流场内各个节点的电荷密度初值以及荷电因子初值;
根据所述电荷密度初值以及所述荷电因子初值确定悬浮液滴空间荷电密度;
在所述电荷密度初值以及所述荷电因子初值的前提下,获取离子空间电荷密度,并根据所述离子空间电荷密度以及所述悬浮液滴空间荷电密度确定计及荷电悬浮液滴影响的泊松方程;
利用有限元法求解所述计及荷电悬浮液滴影响的泊松方程,确定所述荷电因子初值下各个节点处的电场强度计算值;
获取离子迁移率,并根据所述离子迁移率、所述离子空间电荷密度以及所述电场强度计算值确定电流连续性方程;
利用上流元法求解所述电流连续性方程,计算各个节点的空间电荷密度;
判断所述电场强度计算值以及所述空间电荷密度是否均小于第一误差限定条件;若是,根据所述接地电晕笼笼壁处的所述电场强度计算值与所述接地电晕笼笼壁处的电场强度测量值构建误差函数;若否,更新所述电荷密度初值;
判断所述误差函数是否小于第二误差限定条件;若是,确定所述荷电因子初值为当前湿度条件下所述悬浮液滴的荷电因子;若否,更新所述荷电因子初值。
可选的,所述设定离子流场内各个节点的电荷密度初值以及荷电因子初值,具体包括:
根据所述离子流场中悬浮液滴荷电特性测量装置确定计算区域,并对所述计算区域进行剖分,生成剖分后的网格,确定所述剖分后的网格内节点的编号和坐标;所述节点为所述剖分后的网格的端点;
获取所述中心电极的导线参数,并根据所述导线参数确定所述中心电极的表面上各个节点的电荷密度初值以及荷电因子初值;所述导线参数包括导线起晕电压、导线正下方地面处的标称电场、导线的运行电压、导线半径、导线对地边界高度以及导线表面的起晕场强;
将所述接地电晕笼的空间电荷密度解析解作为剩余节点的电荷密度初值。
可选的,所述第一误差限定条件为:
Figure BDA0003207898770000031
Figure BDA0003207898770000032
其中,σρ为电荷密度迭代误差限值;ρn为第n次迭代时导线表面电荷密度值;ρn-1为第n-1次迭代时导线表面电荷密度值;σE为导线表面电场强度迭代误差限值;Emax为导线表面电场强度最大值;Ec为导线起晕场强。
可选的,所述误差函数为:
Figure BDA0003207898770000041
其中,g(Er)为误差函数;Em为所述接地电晕笼笼壁处的电场强度测量值;Er为所述接地电晕笼笼壁处的电场强度计算值。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明通过引入用于表征悬浮液滴荷电特性的参数—荷电因子,提出了一种离子流场中悬浮液滴荷电特性测量装置及方法,用于测量离子流场中悬浮液滴的荷电特性,不需要测量离子流场中悬浮液滴的粒径、数密度等参数,就能够实现不同环境下悬浮液滴荷电特性的计算和表征,对于固定结构的电晕放电平台,通过测量当前环境下地面合成电场,即可准确计算得到离子流场中悬浮液滴荷电因子,以表征离子流场中悬浮液滴的荷电特性,克服了离子流场中悬浮液滴的粒径、数密度等参数难以直接测量的问题,为分析外界环境因素,尤其是湿度对离子流场的影响机理时,提供理论依据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所提供的离子流场中悬浮液滴荷电特性测量装置结构图;
图2为本发明所提供的离子流场中悬浮液滴荷电特性测量方法流程图;
图3为本发明所提供的另一种离子流场中悬浮液滴荷电特性测量方法的具体计算流程图;
图4为本发明所提供的荷电因子计算结果示意图;
图5为本发明所提供的不同湿度条件下悬浮液滴空间电荷密度。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种离子流场中悬浮液滴荷电特性测量装置及方法,能够准确计算得到离子流场中悬浮液滴的荷电因子。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的离子流场中悬浮液滴荷电特性测量装置结构图,如图1所示,一种离子流场中悬浮液滴荷电特性测量装置,包括:接地电晕笼1、中心电极2以及场磨3;所述中心电极2横穿所述接地电晕笼1,所述接地电晕笼1为三段式结构,所述接地电晕笼1的中间段为测量段4,所述接地电晕笼1的两侧为屏蔽段5,所述测量段4内设有离子流板6;所述测量段4用于测量高压直流导线离子流场,所述屏蔽段5用于削弱端部效应;所述场磨3设于所述测量段4的底部,所述场磨3用于测量接地电晕笼笼壁处的电场强度;利用高压直流电源7对所述中心电极2施加电压,使得所述中心电极2处于电晕放电状态,根据所述接地电晕笼笼壁处的电场强度测量离子流场中悬浮液滴荷电特性。
在实际应用中,所述接地电晕笼1为同轴圆柱结构,所述中心电极2为裸铜导线,所述中心电极2横穿所述接地电晕笼1的中心位置。
本发明还包括:支架8以及绝缘子9;所述支架8设于两侧所述屏蔽段5的底部,且所述场磨3设于所述测量段4与所述支架8之间;所述支架8用于支撑所述接地电晕笼1;所述绝缘子9设于所述中心电极2的两端,所述绝缘子9用于支撑所述中心电极2。
图2为本发明所提供的离子流场中悬浮液滴荷电特性测量方法流程图,如图2所示,一种离子流场中悬浮液滴荷电特性测量方法,该离子流场中悬浮液滴荷电特性测量方法应用于如图1所示的离子流场中悬浮液滴荷电特性测量装置。
该离子流场中悬浮液滴荷电特性测量方法包括:
步骤201:设定离子流场内各个节点的电荷密度初值以及荷电因子初值。
所述步骤201具体包括:根据所述离子流场中悬浮液滴荷电特性测量装置确定计算区域,并对所述计算区域进行剖分,生成剖分后的网格,确定所述剖分后的网格内节点的编号和坐标;所述节点为所述剖分后的网格的端点;获取所述中心电极的导线参数,并根据所述导线参数确定所述中心电极的表面上各个节点的电荷密度初值以及荷电因子初值;所述导线参数包括导线起晕电压、导线正下方地面处的标称电场、导线的运行电压、导线半径、导线对地边界高度以及导线表面的起晕场强;将所述接地电晕笼的空间电荷密度解析解作为剩余节点的电荷密度初值。
在实际应用中,使用高压直流电晕为中心电极施加电压,使得中心电极处于电晕放电状态。
使用场磨测量电晕笼笼壁处的电场强度。
确定并剖分计算区域:根据测量装置的结构确定计算区域后,对计算区域进行剖分,得到剖分后各节点的编号和坐标。
步骤202:根据所述电荷密度初值以及所述荷电因子初值确定悬浮液滴空间荷电密度。
步骤202的理论依据如下:
给定各节点电荷密度初值和荷电因子初值:
导线表面各节点电荷密度初值可由式(1)进行求解。其他节点的电荷密度初值可以采用同轴圆柱电极的空间电荷密度解析解。
Figure BDA0003207898770000061
式中,U0为导线起晕电压;Eg为导线正下方地面处的标称电场;U为导线的运行电压;r为导线半径;H为导线对地边界高度;Eon为导线表面的起晕场强。
在所述电荷密度初值以及所述荷电因子初值的前提下,获取所述离子流场中悬浮液滴的液滴参数,并采用饱和场致荷电模型,根据所述液滴参数计算悬浮液滴荷电量;所述液滴参数包括悬浮液滴所在位置的电场强度、悬浮液滴介电常数、悬浮液滴半径以及真空介电常数。
计算悬浮液滴荷电因子初值,首先需要定义荷电因子的表达式。由于周围环境中悬浮液滴的粒径均大于0.2μm,悬浮液滴的荷电方式以场致荷电为主,离子流场中悬浮液滴荷电量达到饱和的时间为0.5s左右,因此可以采用饱和场致荷电模型计算悬浮液滴荷电量:
Figure BDA0003207898770000071
式中,qs为悬浮液滴饱和场致荷电量;E为悬浮液滴所在位置的电场强度;ε为悬浮液滴介电常数;r为悬浮液滴半径;ε0为真空介电常数。
获取悬浮液滴数密度,并根据所述悬浮液滴数密度以及所述悬浮液滴荷电量确定悬浮液滴空间荷电密度。
根据式(3)可得,悬浮液滴空间电荷密度:
Figure BDA0003207898770000072
式中,N为悬浮液滴数密度。
根据所述悬浮液滴空间荷电密度以及所述悬浮液滴所在位置的电场强度确定悬浮液滴的荷电因子。
为了便于分析和计算,定义悬浮液滴的荷电因子:
Figure BDA0003207898770000073
由式(4)可知,荷电因子只与悬浮液滴本身的粒径、介电常数和数密度有关,表征了悬浮液滴在单位电场强度下吸附周围空间电荷的能力。荷电因子越大表明悬浮液滴在单位电场中可以吸附的空间电荷越多,荷电能力越强。
如式(4),基于荷电因子θ、电场强度E以及空间ρw之间的关系,在设定电荷密度初值的基础上对应一个电场强度,根据该电场强度以及设定的荷电因子初值得到悬浮液滴空间荷电密度,获取该悬浮液滴所在位置的电场强度,以该悬浮液滴空间空间荷电密度以及该悬浮液滴所在位置的电场强度确定悬浮液滴的荷电因子。
由此可得,悬浮液滴数密度与悬浮液滴质量密度的关系为:
Figure BDA0003207898770000081
式中,m为悬浮液滴质量浓度;ρ1为悬浮液滴质量密度
由式(6)可得悬浮液滴荷电因子表达式为:
Figure BDA0003207898770000082
取悬浮液滴质量浓度为饱和水含量、粒径取0.2μm时,根据式(6)即可求得悬浮液滴荷电因子的初值。
步骤203:在所述电荷密度初值以及所述荷电因子初值的前提下,获取离子空间电荷密度,并根据所述离子空间电荷密度以及所述悬浮液滴空间荷电密度确定计及荷电悬浮液滴影响的泊松方程。
步骤204:利用有限元法求解所述计及荷电悬浮液滴影响的泊松方程,确定所述荷电因子初值下各个节点处的电场强度计算值。
根据给定的悬浮液滴荷电因子值,根据式(3)计算悬浮液滴空间电荷密度,计算过程中认为悬浮液滴在空间中均匀分布。各节点总空间电荷密度为离子空间电荷密度和悬浮液滴空间电荷密度之和,计算公式如式(7)所示。
ρ=ρew (7)
式中ρ为总空间电荷密度;ρe为离子空间电荷密度;ρw为悬浮液滴空间电荷密度。
使用有限元法求解泊松方程,得到节点电位分布,计算空间电场。
计及荷电悬浮液滴影响的泊松方程为:
Figure BDA0003207898770000083
式中,φ为标量电位。
步骤205:获取离子迁移率,并根据所述离子迁移率、所述离子空间电荷密度以及所述电场强度计算值确定电流连续性方程。
步骤206:利用上流元法求解所述电流连续性方程,计算各个节点的空间电荷密度。
使用上流元法求解电流连续性方程,计算各节点空间电荷密度;
由于离子在电场中的迁移速率比悬浮液滴大两到三个数量级,因此计算过程中不考虑悬浮液滴的运动特性,此时电流连续性方程为:
Figure BDA0003207898770000091
式中J为离子流密度;k为离子迁移率。
步骤207:判断所述电场强度计算值以及所述空间电荷密度是否均小于第一误差限定条件;若是,执行步骤208;若否,执行209。
若满足第一误差限定条件,则计算误差函数值g(Er),若不满足第一误差限定条件,则更新表面电荷密度值,重复步骤201之前的准备工作,即:“使用高压直流电晕为中心电极施加电压,使得中心电极处于电晕放电状态;使用场磨测量电晕笼笼壁处的电场强度;确定并剖分计算区域:根据测量装置的结构确定计算区域后,对计算区域进行剖分,得到剖分后各节点的编号和坐标;给定各节点电荷密度初值和荷电因子初值”,直到满足第一误差限定条件为止。第一误差限定条件为:
Figure BDA0003207898770000092
Figure BDA0003207898770000093
式中σρ为电荷密度迭代误差限值,ρn和ρn-1分别为第n次和第n-1次迭代时导线表面电荷密度值,σE为导线表面电场强度迭代误差限值,Emax和Ec分别为导线表面电场强度最大值和导线起晕场强。
步骤208:根据所述接地电晕笼笼壁处的所述电场强度计算值与所述接地电晕笼笼壁处的电场强度测量值构建误差函数。
误差函数定义为:
Figure BDA0003207898770000101
式中,Em和Er分别为电晕笼笼壁处电场强度测量值和计算值。
步骤209:更新所述电荷密度初值。
步骤210:判断所述误差函数是否小于第二误差限定条件;若是,执行步骤211;若否,执行步骤212。
步骤211:确定所述荷电因子初值为当前湿度条件下所述悬浮液滴的荷电因子。
步骤212:更新所述荷电因子初值。
若计算的误差函数值g(Er)≤0.01%,则认为该荷电因子为当前湿度条件下悬浮液滴的荷电因子,否则更新荷电因子值,直到满足误差函数值g(Er)≤0.01%为止。图3为本发明所提供的另一种离子流场中悬浮液滴荷电特性测量方法的具体计算流程图,如图3所示。
将具体数值赋值到本发明中时,具体实施方案如下:
具体实施过程中使用的测量装置中心电极为裸铜导线,长度为2.8m,半径为1.1mm。接地电晕笼测量段长2.0m,一侧屏蔽段长0.3m,半径为0.4m。
分别在30%-90%相对湿度条件下,使用高压直流电源对中心电极施加55kV正极性电压,保证中心电极处于电晕放电状态,使用直流场强仪测量电晕笼处的电场强度。
本发明在计算过程中基于以下假设:
(1)计算过程中认为周围空气中的悬浮液滴均匀分布。
(2)在相同的电场中,离子的迁移速率一般高出荷电悬浮液滴迁移速率2-3个数量级,因此计算过程中认为悬浮液滴处于静止状态,不考虑悬浮液滴的动态特性。
(3)由于悬浮液滴在整个空气中所占的质量分数较小,因此认为悬浮液滴并不对空气介电常数产生影响。
在以上假设的基础上,同轴圆柱电极的离子流场控制方程和边界条件为:
Figure BDA0003207898770000111
式中,φ为标量电位;E为电场强度;ρe为离子空间电荷密度;ρw为悬浮液滴空间电荷密度;J为离子流密度;k为离子迁移率;ε0为真空介电常数。
控制方程的边界条件如下:
导线表面电压为运行电压:
Figure BDA0003207898770000112
电晕笼电位为0:
Figure BDA0003207898770000113
导线表面电场强度为起晕场强:
EL=E0 (16)
在以上理论的基础上,悬浮液滴荷电因子的测量方法如下:
(1)确定计算区域并对计算区域进行剖分。
(2)给定各节点电荷密度初值和荷电因子初值。
(3)计算悬浮液滴空间电荷密度和各节点总空间电荷密度。
(4)使用有限元法求解泊松方程,得到节点电位分布,计算空间电场。
(4)使用上流元法求解电流连续性方程,计算各节点空间电荷密度。
(5)若满足误差限定条件,则计算误差函数值g(Er),若不满足误差限定条件,则更新表面电荷密度值,重复步骤2-4,直到满足误差限定条件为止。
误差限定条件为:
Figure BDA0003207898770000121
Figure BDA0003207898770000122
式中,σρ为电荷密度迭代误差限值,ρn和ρn-1分别为第n次和第n-1次迭代时导线表面电荷密度值,σE为导线表面电场强度迭代误差限值,Emax和Ec分别为导线表面电场强度最大值和导线起晕场强。
误差函数定义为:
Figure BDA0003207898770000123
式中,Em和Er分别为电晕笼笼壁处电场强度测量值和计算值。
(6)若计算的误差函数值g(Er)≤0.01%,则认为该荷电因子为当前湿度条件下悬浮液滴的荷电因子,否则更新荷电因子值,重复步骤2-5,直到满足误差函数值g(Er)≤0.01%为止。
根据荷电因子计算方法,分别计算不同湿度条件下,离子流场中悬浮液滴的荷电因子,计算结果如图4所示,并使用线性函数分别对低湿和高湿环境下荷电因子的进行拟合。
由拟合曲线可知,在低湿条件下悬浮液滴荷电因子随湿度变化较小,RH=50%时荷电因子为0.032×10-12C/(V·m2),仅为RH=30%时荷电因子的1.19倍。说明在低湿环境中悬浮液滴含量较少,且随周围环境湿度变化不大。在高湿条件下,悬浮液滴荷电因子随湿度线性增加,RH=90%时荷电因子为13.23×10-12C/(V·m2),为RH=80%时荷电因子的1.47倍。两个拟合曲线的交点横坐标约为60%,说明当周围环境相对湿度大于60%时,空间中会明显出现悬浮液滴,悬浮液滴含量随着湿度的增大迅速增多。
使用计算得到的荷电因子和式(13)分别计算不同湿度条件下,悬浮液滴的空间电荷密度,计算结果如图5所示。由计算结果可知,随着湿度的增大,悬浮液滴空间电荷密度逐渐增大。
通过测量结果可知,本发明公开的测量装置和方法,可以用来测量不同湿度条件下离子流场中悬浮液滴的荷电特性,为明确湿度对离子流场的影响规律提供数据支撑和理论支持。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (4)

1.一种离子流场中悬浮液滴荷电特性测量方法,其特征在于,所述离子流场中悬浮液滴荷电特性测量方法应用于一种离子流场中悬浮液滴荷电特性测量装置,所述装置包括:接地电晕笼、中心电极以及场磨;所述中心电极横穿所述接地电晕笼,所述接地电晕笼为三段式结构,所述接地电晕笼的中间段为测量段,所述接地电晕笼的两侧为屏蔽段,所述测量段内设有离子流板;所述测量段用于测量高压直流导线离子流场,所述屏蔽段用于削弱端部效应;所述场磨设于所述测量段的底部,所述场磨用于测量接地电晕笼笼壁处的电场强度;利用高压直流电源对所述中心电极施加电压,使得所述中心电极处于电晕放电状态,根据所述接地电晕笼笼壁处的电场强度测量离子流场中悬浮液滴荷电特性;
所述离子流场中悬浮液滴荷电特性测量方法包括:
设定离子流场内各个节点的电荷密度初值以及荷电因子初值;
根据所述电荷密度初值以及所述荷电因子初值确定悬浮液滴空间荷电密度;
在所述电荷密度初值以及所述荷电因子初值的前提下,获取离子空间电荷密度,并根据所述离子空间电荷密度以及所述悬浮液滴空间荷电密度确定计及荷电悬浮液滴影响的泊松方程;
利用有限元法求解所述计及荷电悬浮液滴影响的泊松方程,确定所述荷电因子初值下各个节点处的电场强度计算值;
获取离子迁移率,并根据所述离子迁移率、所述离子空间电荷密度以及所述电场强度计算值确定电流连续性方程;
利用上流元法求解所述电流连续性方程,计算各个节点的空间电荷密度;
判断所述电场强度计算值以及所述空间电荷密度是否均小于第一误差限定条件;若是,根据所述接地电晕笼笼壁处的所述电场强度计算值与所述接地电晕笼笼壁处的电场强度测量值构建误差函数;若否,更新所述电荷密度初值;
判断所述误差函数是否小于第二误差限定条件;若是,确定所述荷电因子初值为当前湿度条件下所述悬浮液滴的荷电因子;若否,更新所述荷电因子初值。
2.根据权利要求1所述的离子流场中悬浮液滴荷电特性测量方法,其特征在于,所述设定离子流场内各个节点的电荷密度初值以及荷电因子初值,具体包括:
根据所述离子流场中悬浮液滴荷电特性测量装置确定计算区域,并对所述计算区域进行剖分,生成剖分后的网格,确定所述剖分后的网格内节点的编号和坐标;所述节点为所述剖分后的网格的端点;
获取所述中心电极的导线参数,并根据所述导线参数确定所述中心电极的表面上各个节点的电荷密度初值以及荷电因子初值;所述导线参数包括导线起晕电压、导线正下方地面处的标称电场、导线的运行电压、导线半径、导线对地边界高度以及导线表面的起晕场强;
将所述接地电晕笼的空间电荷密度解析解作为剩余节点的电荷密度初值。
3.根据权利要求1所述的离子流场中悬浮液滴荷电特性测量方法,其特征在于,所述第一误差限定条件为:
Figure FDA0003593338780000021
Figure FDA0003593338780000022
其中,σρ为电荷密度迭代误差限值;ρn为第n次迭代时导线表面电荷密度值;ρn-1为第n-1次迭代时导线表面电荷密度值;σE为导线表面电场强度迭代误差限值;Emax为导线表面电场强度最大值;Ec为导线起晕场强。
4.根据权利要求3所述的离子流场中悬浮液滴荷电特性测量方法,其特征在于,所述误差函数为:
Figure FDA0003593338780000023
其中,g(Er)为误差函数;Em为所述接地电晕笼笼壁处的电场强度测量值;Er为所述接地电晕笼笼壁处的电场强度计算值。
CN202110922336.3A 2021-08-12 2021-08-12 一种离子流场中悬浮液滴荷电特性测量装置及方法 Active CN113625067B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110922336.3A CN113625067B (zh) 2021-08-12 2021-08-12 一种离子流场中悬浮液滴荷电特性测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110922336.3A CN113625067B (zh) 2021-08-12 2021-08-12 一种离子流场中悬浮液滴荷电特性测量装置及方法

Publications (2)

Publication Number Publication Date
CN113625067A CN113625067A (zh) 2021-11-09
CN113625067B true CN113625067B (zh) 2022-06-10

Family

ID=78384730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110922336.3A Active CN113625067B (zh) 2021-08-12 2021-08-12 一种离子流场中悬浮液滴荷电特性测量装置及方法

Country Status (1)

Country Link
CN (1) CN113625067B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114113815B (zh) * 2021-11-19 2022-09-16 华北电力大学 一种颗粒物荷电特性测量方法及系统
CN117074801B (zh) * 2023-10-14 2024-02-13 之江实验室 一种悬浮带电微球测量电场的装置及方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007073920A (ja) * 2005-08-09 2007-03-22 Kansai Electric Power Co Inc:The 避雷器用非線形抵抗素子
CN102609588A (zh) * 2012-02-10 2012-07-25 华北电力大学 一种离子流场的计算方法
CN102608394A (zh) * 2012-04-01 2012-07-25 华北电力大学 一种线路电晕电流测量系统及方法
CN105572510A (zh) * 2016-01-23 2016-05-11 清华大学 一种基于电晕笼的直流输电线路电晕损失测试方法
CN106054049A (zh) * 2016-06-17 2016-10-26 华北电力大学 一种表面积污导线、可控的导线表面积污的方法及系统
CN108846193A (zh) * 2018-06-08 2018-11-20 华北电力大学(保定) 一种考虑实时气候条件的高压直流输电线路合成场强和离子流密度的计算方法
CN112182920A (zh) * 2020-08-24 2021-01-05 中国电力科学研究院有限公司 一种确定直流输电线路合成电场的场强值的迭代方法
CN112444686A (zh) * 2020-08-24 2021-03-05 中国电力科学研究院有限公司 一种确定交叉跨域直流线路的电荷密度初始值的方法
CN112649676A (zh) * 2020-12-18 2021-04-13 国网西藏电力有限公司电力科学研究院 一种考虑风速影响的高海拔地区混合场强的计算方法
CN113033117A (zh) * 2021-03-09 2021-06-25 江苏大学 一种运动荷电液滴诱导电场强度及电场力计算方法和系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012095205A2 (en) * 2011-01-11 2012-07-19 Siemens Aktiengesellschaft Apparatus and method for removal of particulate matter from a gas
US20170322224A1 (en) * 2016-03-12 2017-11-09 Andreas Hieke Rigorous method and apparatuses for the analysis of complex mixtures of organic molecules with an enhanced degree of information extraction
CN107704655A (zh) * 2017-09-05 2018-02-16 国网浙江省电力公司电力科学研究院 一种直流输电导线污秽评估方法及系统
CN108594029B (zh) * 2018-04-26 2020-05-19 大连理工大学 基于图像识别的粉体颗粒荷电量测量装置和方法
CN108984819A (zh) * 2018-05-24 2018-12-11 华北电力大学(保定) 一种计算雾霾条件下高压直流输电线路电晕离子流场的方法
CN109117458A (zh) * 2018-06-28 2019-01-01 浙江省电力有限公司电力科学研究院 一种基于改进无网格法的直流设备积污特性计算方法
CN109344470A (zh) * 2018-09-17 2019-02-15 上海电力学院 一种高压直流离子流场数学模型获取方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007073920A (ja) * 2005-08-09 2007-03-22 Kansai Electric Power Co Inc:The 避雷器用非線形抵抗素子
CN102609588A (zh) * 2012-02-10 2012-07-25 华北电力大学 一种离子流场的计算方法
CN102608394A (zh) * 2012-04-01 2012-07-25 华北电力大学 一种线路电晕电流测量系统及方法
CN105572510A (zh) * 2016-01-23 2016-05-11 清华大学 一种基于电晕笼的直流输电线路电晕损失测试方法
CN106054049A (zh) * 2016-06-17 2016-10-26 华北电力大学 一种表面积污导线、可控的导线表面积污的方法及系统
CN108846193A (zh) * 2018-06-08 2018-11-20 华北电力大学(保定) 一种考虑实时气候条件的高压直流输电线路合成场强和离子流密度的计算方法
CN112182920A (zh) * 2020-08-24 2021-01-05 中国电力科学研究院有限公司 一种确定直流输电线路合成电场的场强值的迭代方法
CN112444686A (zh) * 2020-08-24 2021-03-05 中国电力科学研究院有限公司 一种确定交叉跨域直流线路的电荷密度初始值的方法
CN112649676A (zh) * 2020-12-18 2021-04-13 国网西藏电力有限公司电力科学研究院 一种考虑风速影响的高海拔地区混合场强的计算方法
CN113033117A (zh) * 2021-03-09 2021-06-25 江苏大学 一种运动荷电液滴诱导电场强度及电场力计算方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Finite Volume-Based Approach for the Hybrid;Han Yin, Jinliang He,Bo Zhang, Rong Zeng;《IEEE TRANSACTIONS ON POWER DELIVERY》;20110922;全文 *
基于气象数据的雾天气条件下;申南轩,张远航,徐鹏,卢铁兵;《中国电机工程学报》;20201205;全文 *

Also Published As

Publication number Publication date
CN113625067A (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
CN113625067B (zh) 一种离子流场中悬浮液滴荷电特性测量装置及方法
Malik A review of the charge simulation method and its applications
Bian et al. Corona-generated space charge effects on electric field distribution for an indoor corona cage and a monopolar test line
Abouelatta et al. Measurement and assessment of corona current density for HVDC bundle conductors by FDM integrated with full multigrid technique
Ziedan et al. Onset voltage of corona discharge in wire-duct electrostatic precipitators
CN108828292B (zh) 考虑接地排影响的氧化锌避雷器阻性电流在线监测方法
Talaat Electrostatic field calculation in air gaps with a transverse layer of dielectric barrier
Al-Hamouz et al. Finite-element solution of monopolar corona on bundle conductors
Kontargyri et al. Measurement and simulation of the voltage distribution and the electric field on a glass insulator string
Bian et al. Comparative study on corona-generated audible noise and radio noise of ac long-term operating conductors with two bundle types
Dumitran et al. 2-D corona field computation in configurations with ionising and non-ionising electrodes
Zhang et al. Experimental investigation on humidity effects on the variations of positive DC corona discharge
Baek et al. Experiment and analysis for effect of floating conductor on electric discharge characteristic
Dumitran et al. Computational and experimental study of ionic space charge generated by combined corona–electrostatic electrode systems
CN111753454B (zh) 一种正极性连续电晕放电电流的仿真方法
Ding et al. Discharge simulation of typical air gap considering dynamic boundary and charge accumulation
Kim et al. Influence of multiple insulator rods on potential and electric field distributions at their surface
Amrani et al. Modified Peek formula for calculating positive DC corona inception voltage on polluted insulator
Kontargyri et al. Calculation of the electric field on an insulator string using the finite elements method
Salceanu et al. Upon the influence of charge image on the electric field intensity
Al-Hamouz Corona power loss on bundled conductors: experimental and computational results
CN113376444B (zh) 一种基于离子流密度计算导线起晕场强的方法及系统
Nicolopoulou et al. Electric field and voltage distribution around composite insulators
CN109408868A (zh) 一种变压器绕组空间电场计算方法
Fylladitakis et al. Analytical estimation of the electrostatic field in cylinder-plane and cylinder-cylinder electrode configurations

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