本发明专利技术涉及一种声脉冲激励下确定空间电荷密度分布的方法，首先对待测空间进行网格剖分形成单元，然后计算单元内空间电荷振动产生的电场信号，最后确定连续声波调制下待测空间内空间电荷振动产生的电场信号。该电场信号解耦算法，用电场信号、空间电荷与声脉冲之间的数值关系代替三者之间的函数关系，一方面避免了上述问题，另一方面，有利于后期三维空间电荷测量系统的设计。
【技术实现步骤摘要】
一种声脉冲激励下确定空间电荷密度分布的方法
本专利技术涉及空间电荷测量
，具体指一种声脉冲激励下确定空间电荷密度分布的方法。
技术介绍
伴随特高压直流输电的迅猛发展，特高压直流输电电磁环境等相关问题受到了越来越多人的关注。电晕放电是电磁环境的一个重要组成部分。正常工作电压下的输电线路通常都会发生电晕效应。由于输电线电压极性固定，直流电晕产生的离子在背离导线的方向上运动，充斥着整个输电线空间形成离子流。输电线周围空间电场是输电导线产生标称场强和离子流场共同作用的结果。电晕放电离子流场数值计算的方法已经相当完备，但少有学者从实验测量上对计算结果进行验证。而关于声脉冲激励下对电场信号及空间电荷的分布方面，Hazmi推导了与声脉冲激励下电场信号与声脉冲及空间电荷之间的函数关系，但是存在如下问题：(1)声脉冲传播防线的假设过于理想化，声脉冲激发产生的声场很难按照理想圆锥体形状分布；(2)对于不均匀结构的电极，电场信号与声脉冲及空间电荷之间的函数关系式不再适用；(3)多数情况下，函数表达式很难描述空间电荷分布。
技术实现思路
针对现有技术存在的上述问题，本专利技术的目的是提供一种能使用于多种形状电极的，声脉冲激励下确定空间电荷密度分布的方法。为实现上述目的，本专利技术采用如下技术方案：一种声脉冲激励下确定空间电荷密度分布的方法，包括如下步骤：S1：采用有限差分法对待测空间进行网格剖分，将待测空间分成多个单元Ni；S2：对步骤S1形成的单元Ni进行剖分，形成多个亚单元Ni-abc；S3：计算电场信号：S3a：计算第i个单元Ni内空间电荷振动产生的电场信号：每个单元Ni内空间电荷产生的电场信号等于构成该单元的所有亚单元内空间电荷产生的电场信号的叠加，具体如下：电场信号在声信号作用下的变化表示为：其中，ρ表示电场的电荷密度，ρi表示第i个单元内的空间电荷密度，A(x)是声波幅值，α是衰减参数，为常数，γ表示定压比热容，j为虚数单位，没有实际含义，ω是声波的角频率，K是波系数，为常数，r表示场点距源点的距离；其中，Ei-abc表示亚单元内空间电荷产生的电场信号，eR表示场点到源点的向量，qi-abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离；Eix-abc，Eiy-abc和Eiz-abc分别是Ei-abc在X，Y和Z轴上的投影，ρi表示第i个单元内的空间电荷密度，A(x)表示声压幅值，P0是标准大气压，ε0表示空气的介电常数，ti表示第i个单元开始振动时间，eR表示场点到源点的向量，qi-abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离，ΔX、ΔY、ΔZ表示R在X，Y和Z轴上的投影，ΔV表示组成第i个单元的一个亚单元的体积；令Ci,Px,Py和Pz是为了方便表示而引入的参量，不具有实际含义；第i个单元内的电场信号是组成其的所有亚单元内空间电荷振动产生的电压信号的叠加：Ei-x＝∑Eix-abc(7)；Ei-y＝∑Eiy-abc(8)；Ei-z＝∑Eiz-abc(9)；其中Ei-x，Ei-y和Ei-z分别表示第i个单元内空间电荷振动产生的电场信号在X，Y和Z轴上的投影，待测空间中第i个单元内电场信号表示为：：令Di-x＝∑Px(11)；Ei＝CiρiDi(13)；其中，Ei表示第i个单元内空间电荷产生的电场信号，Di是为了方便表示而引入的参量，不具有实际含义，Di-x，Di-y和Di-z分别表示Di在X，Y和Z轴上的投影；S3b：连续声波调制下待测空间内空间电荷振动产生的电场信号表示为：其中，ki＝CiDi，E'(ti)表示ti时刻待测空间内的电场信号，T和f分别表示声波的周期和频率；S4：求解公式(14)，得到在连续声波调制下待测空间内空间电荷的密度分布。相对于现有技术，本专利技术具有如下优点：本专利技术提出了改进后的电场信号解耦算法，用电场信号、空间电荷与声脉冲之间的数值关系代替三者之间的函数关系，一方面避免了上述问题，另一方面，有利于后期三维空间电荷测量系统的设计。本专利技术首先引入有限差分法对待测区域进行了网格剖分，避免了由于待测区域形状多样性给积分算法带来的麻烦；其次剖分后，待测区域每一点的电荷密度大小更方便的获取和表示；最后幅值随距离的变化通过实验获取，更加精确，同时声场的分布特性也更加明确。具体实施方式下面对本专利技术作进一步详细说明。一种声脉冲激励下确定空间电荷密度分布的方法，具体包括如下步骤：S1：采用有限差分法对待测空间进行网格剖分，形成多个单元Ni。S2：对步骤S1形成的单元Ni进行剖分，形成多个亚单元Ni-abc；对单元Ni进行剖分可以有效提高求解函数的精度，该步骤中也采用有限差分法对单元Ni进行剖分；S3：计算电场信号：S3a：计算第i个单元Ni内空间电荷振动产生的电场信号：每个单元Ni内空间电荷产生的电场信号等于构成该单元的所有亚单元内空间电荷产生的电场信号的叠加，具体如下：电场信号在声信号作用下的变化表示为：其中，ρ表示电场的电荷密度，ρi表示第i个单元内的空间电荷密度，即第i个单元内电场稳态时的电荷密度(组成每个单元的所有亚单元中电荷密度相同，即第i个单元的空间电荷密度即是组成该单元的亚单元的空间电荷密度)，A(x)是声波幅值(声波幅值的大小随声源到场点距离的增加而减少)，α是衰减参数，为常数，γ表示定压比热容，j为虚数单位，没有实际含义，ω是声波的角频率，K是波系数，为常数，K＝ω/c，c表示声速；r表示场点距源点的距离，其中，场点是指待测空间中的某一点，源点指发出的声波的点，即声波发生装置所在的位置；其中，Ei-abc表示亚单元内空间电荷产生的电场信号，eR表示场点到源点的向量，qi-abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离；Eix-abc，Eiy-abc和Eiz-abc分别是Ei-abc在X，Y和Z轴上的投影，ρi表示第i个单元内的空间电荷密度，A(x)表示声压幅值，P0是标准大气压，ω是声波的角频率，ti表示第i个单元开始振动时间，ε0表示空气的介电常数，ΔX、ΔY、ΔZ表示R在X，Y和Z轴上的投影；ΔV表示组成第i个单元的一个亚单元的体积，组成第i个单元的所有亚单元的体积均相等；令Ci,Px,Py和Pz是为了方便表示而引入的参量，不具有实际含义；第i个单元内的电场信号是组成其的所有亚单元内空间电荷振动产生的电压信号的叠加：Ei-x＝∑Eix-abc(22)；Ei-y＝∑Eiy-abc(23)；Ei-z＝ΣEiz-abc(24)；其中Ei-x，Ei-y和Ei-z分别表示第i个单元内空间电荷振动产生的电场信号在X，Y和Z轴上的投影，待测空间中第i个单元内电场信号表示为：令Di-x＝∑Px(26)；Ei＝CiρiDi(28)；其中，Ei表示第i个单元内空间电荷振动产生的电场信号，Di是为了方便表示而引入的参量，不具有实际含义，Di-x，Di-y和Di-z分别表示Di在X，Y和Z轴上的投影；S3b：连续声波调制下待测空间内空间电荷振动产生的电场信号表示为：其中，ki＝CiDi本文档来自技高网...
【技术保护点】
一种声脉冲激励下确定空间电荷密度分布的方法，其特征在于，包括如下步骤：S1：采用有限差分法对待测空间进行网格剖分，将待测空间分成多个单元Ni；S2：对步骤S1形成的单元Ni进行剖分，形成多个亚单元Ni‑abc；S3：计算电场信号：S3a：计算第i个单元Ni内空间电荷振动产生的电场信号：每个单元Ni内空间电荷产生的电场信号等于构成该单元的所有亚单元内空间电荷产生的电场信号的叠加，具体如下：电场信号在声信号作用下的变化表示为：ρ=ρi[1+A(x)γP0e-αγ·ej(ωt-Kr)],(i=1,2,3,4,...n)---(1);]]>其中，ρ表示电场的电荷密度，ρi表示第i个单元内的空间电荷密度，A(x)是声波幅值，α是衰减参数，为常数，γ表示定压比热容，j为虚数单位，没有实际含义，ω是声波的角频率，K是波系数，为常数，r表示场点距源点的距离；Ei-abc=qi-abc4πϵ0R2eR---(2);]]>Eix-abc=ρiA(x)γP0cos[ω(t-ti)]·ΔV·ΔX4πϵ0R3---(3);]]>Eiy-abc=ρiA(x)γP0cos[ω(t-ti)]·ΔV·ΔY4πϵ0R3---(4);]]>Eiz-abc=ρiA(x)γP0cos[ω(t-ti)]·ΔV·ΔZ4πϵ0R3---(5);]]>其中，Ei‑abc表示亚单元内空间电荷产生的电场信号，eR表示场点到源点的向量，qi‑abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离；Eix‑abc，Eiy‑abc和Eiz‑abc分别是Ei‑abc在X，Y和Z轴上的投影，ρi表示第i个单元内的空间电荷密度，A(x)表示声压幅值，P0是标准大气压，ε0表示空气的介电常数，ti表示第i个单元开始振动时间，eR表示场点到源点的向量，qi‑abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离，ΔX、ΔY、ΔZ表示R在X，Y和Z轴上的投影，ΔV表示组成第i个单元的一个亚单元的体积；令Ci=A(x)·ΔV4πϵ0γP0,]]>Px=ΔXR3cos[ω(t-ti)]Py=ΔYR3cos[ω(t-ti)]Pz=ΔZR3cos[ω(t-ti)];]]>Ci,Px,Py和Pz是为了方便表示而引入的参量，不具有实际含义；Eix-abc=Ci·Px·ρiEiy-abc=Ci·Py·ρiEiz-abc=Ci·Pz·ρi---(6);]]>第i个单元内的电场信号是组成其的所有亚单元内空间电荷振动产生的电压信号的叠加：Ei‑x＝ΣEix‑abc (7)；Ei‑y＝ΣEiy‑abc (8)；Ei‑z＝ΣEiz‑abc (9)；其中Ei‑x，Ei‑y和Ei‑z分别表示第i个单元内空间电荷振动产生的电场信号在X，Y和Z轴上的投影，待测空间中第i个单元内电场信号表示为：Ei=Ei-x2+Ei-y2+Ei-z2=(ΣEix-abc)2+(ΣEiy-abc)2+(ΣEiz-abc)2=(ΣPx)2+(ΣPy)2+(ΣPz)2·Ci·ρi---(10);]]>令Di‑x＝ΣPx (11)；Di=Di-x2+Di-y2+Di-z2---(12);]]>Ei＝CiρiDi (13)；其中，Ei表示第i个单元内空间电荷产生的电场信号，Di是为了方便表示而引入的参量，不具有实际含义，Di‑x，Di‑y和Di‑z分别表示Di在X，Y和Z轴上的投影；S3b：连续声波调制下待测空间内空间电荷振动产生的电场信号表示为：ti=i&Center...
【技术特征摘要】
1.一种声脉冲激励下确定空间电荷密度分布的方法，其特征在于，包括如下步骤：S1：采用有限差分法对待测空间进行网格剖分，将待测空间分成多个单元Ni；S2：对步骤S1形成的单元Ni进行剖分，形成多个亚单元Ni-abc；S3：计算电场信号：S3a：计算第i个单元Ni内空间电荷振动产生的电场信号：每个单元Ni内空间电荷产生的电场信号等于构成该单元的所有亚单元内空间电荷产生的电场信号的叠加，具体如下：电场信号在声信号作用下的变化表示为：其中，ρ表示电场的电荷密度，ρi表示第i个单元内的空间电荷密度，A(x)是声波幅值，α是衰减参数，为常数，γ表示定压比热容，j为虚数单位，没有实际含义，ω是声波的角频率，K是波系数，为常数，r表示场点距源点的距离；其中，Ei-abc表示亚单元内空间电荷产生的电场信号，eR表示场点到源点的向量，qi-abc表示组成第i个单元的一个亚单元内空间电荷的电量，R表示亚单元与接受待测空间内电荷振动形成电场信号的电场天线之间的距离；Eix-abc，Eiy-abc和Eiz-abc分别是Ei-abc在X，Y和Z轴上的投影，ρi表示第i个单元内的空间电荷密度，A(x)表示声压幅...
【专利技术属性】
技术研发人员：廖瑞金，赵学童，王平，杨帆，刘康淋，伍飞飞，刘宏波，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;85
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。