CN111624389B - 一种用于倾斜多导体系统的电流测量系统及方法 - Google Patents

一种用于倾斜多导体系统的电流测量系统及方法 Download PDF

Info

Publication number
CN111624389B
CN111624389B CN202010501833.1A CN202010501833A CN111624389B CN 111624389 B CN111624389 B CN 111624389B CN 202010501833 A CN202010501833 A CN 202010501833A CN 111624389 B CN111624389 B CN 111624389B
Authority
CN
China
Prior art keywords
magnetic field
conductor
current
formula
sensor
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
CN202010501833.1A
Other languages
English (en)
Other versions
CN111624389A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202010501833.1A priority Critical patent/CN111624389B/zh
Publication of CN111624389A publication Critical patent/CN111624389A/zh
Application granted granted Critical
Publication of CN111624389B publication Critical patent/CN111624389B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0092Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring current only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/005Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Abstract

本发明涉及一种用于倾斜多导体系统的电流测量系统及方法,属于电气技术领域。该方法包括以下步骤:步骤S1,建立空间磁场与待测电流的逆问题模型;步骤S2,测得空间6个点的三维磁场;步骤S3,将磁场测量值代入逆问题模型,采用进化算法逆推得到各导体的坐标位置,倾斜角以及各相电流;步骤S4,利用所计算得到的坐标位置和倾斜角,计算系数矩阵,完成电流传感器的校准过程;步骤S5,电流传感器的实际应用过程中,将所测磁场值所构成的列向量乘以系数矩阵的伪逆即可得到各相导体的电流。最后,通过对比仿真,给出了仅利用二维平面的磁场测量值(即未考虑导体倾斜角)和本发明提出的方法所计算得到的各相电流误差,进而证明本专利提出的方法有效提升了电流的测量精度,具有实际的应用价值。

Description

一种用于倾斜多导体系统的电流测量系统及方法
技术领域
本发明属于电气技术领域,涉及一种用于倾斜多导体系统的电流测量系统及方法。
背景技术
电流是诸多应用领域中不可或缺的电气参数,如潮流计算,继电保护和故障诊断。惯用的电流测量方法基于电磁感应定律的电磁式测量,即电流互感器(CurrentTransformer,CT)。然而,由于铁芯的存在,磁芯饱和现象会导致二次检测信号畸变,甚至导致电力系统铁磁谐振,造成系统不稳定。近年来,基于磁场传感器阵列的电流测量方法备受学者的关注,并研究出了针对单根导体以及多根导体的电流测量方法。然而,所报道的方法中,大多假设导体与磁场传感器阵列所在平面垂直,鲜有考虑当不垂直时,电流测量误差是否远远超过可接受的误差范围。基于上述背景,本发明提出了基于磁场传感器阵列的倾斜多导体系统的电流测量方法,从电流传感器的结构到电流计算算法再到误差分析等角度,来说明本发明提出的电流测量方法用于电流测量的有效性。
发明内容
有鉴于此,本发明的目的在于提供一种用于倾斜多导体系统的电流测量系统及方法。
为达到上述目的,本发明提供如下技术方案:
一种用于倾斜多导体系统的电流测量系统,该系统包括电流传感器,数据采集卡和上位机,上述电流传感器包括磁场检测组件、差分放大器,且都焊接在PCB板上;
所述电流传感器呈圆环形结构,多根导体由中空圆孔穿过,导体电流在空间产生的磁场由磁场检测组件测量得到;
所述磁场检测组件包括6个TMR三维磁场传感器,呈均匀的圆周分布;以电流传感器的中心作为坐标原点,建立三维的笛卡尔坐标系;TMR三维磁场传感器所测三个磁场分量分别朝向x,y,z轴方向;
所述TMR三维磁场传感器的输出信号以差动模式输出,为抑制共模干扰,后续信号处理电路将传感器输出信号差分放大转化为单端信号后,再以差分方式输出,保证信号都以差分方式传输,从而提升电路的噪声抑制能力;
所述差动放大器包括2个运算放大器,分别记为放大器1和放大器2,磁场传感器的输出信号为差分模式,放大器1以差分方式接收上述输出信号,并以一定增益放大,从而将双端信号转化为单端信号,随后,放大器2上述单端信号再转化为双端信号,即以差分模式输出;
所述数据采集卡为32路单端或16路差分模拟输入,采样率为100kS/s,采用16路差分模拟输入模式,采用2个上述采集卡模块来采集6个TMR磁场传感器的18路输出信号,并将信号显示在上位机系统;
所述上位机系统与数据采集卡进行数据传输,后续利用所采集的磁场数据进行逆问题计算来求解流经各相导体的电流,同时将数据可视化;
测量时,多导体穿过所述中空圆孔。
可选的,所述导体电流与其所产生磁场的关系以矩阵形式表示为
Figure BDA0002525030270000021
当系数矩阵确定后,电流[I]通过磁场[B]乘以系数矩阵[M]的伪逆得到;电流传感器的校准过程即求解系数矩阵[M]。
可选的,所述校准过程包括以下步骤:
步骤S1:根据无限长导体在空间产生的磁场与电流的关系,当单根导体与磁场传感器阵列倾斜时,建立流经导体的电流与空间磁场的数学模型;
步骤S2:由步骤S1所得结果,利用叠加原理,得到三相倾斜导体与空间磁场的数学模型;
步骤S3:由所述磁场检测组件测量空间6个点的三维磁场分量;
步骤S4:利用步骤S3所测磁场值,并结合步骤S2的逆问题模型,采用进化算法逆推各导体的坐标位置,倾斜角和流经导体的电流;
步骤S5:由步骤S4得到的导体坐标位置和倾斜角计算系数矩阵[M],最终完成电流传感器的校准过程。
可选的,所述步骤S1中,计算公式为:
Figure BDA0002525030270000022
式中,Bs,Bsx,Bsy,Bsz分别为空间某一点处的合成磁场,x,y,z方向的磁场分量;(xs,ys)为磁场传感器阵列所在平面的坐标位置,(xc,yc)为导体与磁场传感器阵列所在平面的交点坐标,I为流经导体的电流,m=sin(α)cos(β),n=sin(α)sin(β),p=cos(α),α,β分别为倾斜导体与z轴方向的夹角,倾斜导体在xoy平面的投影与x轴方向的夹角。
可选的,所述步骤S2中,计算公式为:
Figure BDA0002525030270000031
由于每个磁场传感器测量得到的是x,y,z轴三个方向的磁场,因此分别对三个方向的磁场进行叠加,从而得到公式(2)所示的三种形式的以X为变量的方程;公式(2)中,αi,βi(i=1,2,3)为三相导体的倾斜角,xci,yci(i=1,2,3)分别为三相导体的坐标位置,Ii(i=1,2,3)分别为流经三相导体的电流;Fi(X),Gi(X),Hi(X)的表达式如公式(3);式(3)中的f(i,j)表达式由式(4)表示;与倾斜角相关的参数如式(5)所示;
Figure BDA0002525030270000032
f(i,j)=pj 2(ysi-ycj)2+pj 2(xsi-xcj)2+[mj(ysi-ycj)-nj(xsi-xcj)]2 (4)
Figure BDA0002525030270000033
可选的,在所述步骤S1和S2中,得到倾斜三相导体电流与其在空间中所产生的磁场关系的数学模型;由公式(2),含有的未知变量数为15,分别为导体倾斜角αii(i=1,2,3),导体的坐标位置xci,yci以及各相导体的电流Ii(i=1,2,3);采用6个三维磁场传感器,测量得到18个磁场分量,由叠加定理确定18个方程;采用进化算法进行多目标寻优时,非线性问题所满足的约束条件为
|I1|,|I2|,|I3|<MAXcurrent (6)
Figure BDA0002525030270000041
Figure BDA0002525030270000042
上式中,MAXcurrent,MAXradius分别为导体允许流过的电流,导体所在限制区域的最大半径;
涉及的目标函数为
Figure BDA0002525030270000043
上式中,Bix,Biy,Biz为磁场测量值,
Figure BDA0002525030270000044
为由式(2)计算所得磁场值,要求所计算的电流不能超过允许值,因此添加了一个惩罚项,并将惩罚因子设置为1000;进化算法的目的就是通过群体的进化迭代得到最优个体,使得式(9)的值尽量小;
可选的,所述步骤S4中,确定几何参数α11,xc1,yc122,xc2,yc233,xc3,yc3后,公式(2)写成矩阵形式,即
Figure BDA0002525030270000045
上式中,[I]为待求的导体电流,[B]为磁场传感器的测量值,[M]为系数矩阵。
电流通过下式得到
[I]3×1=pseudo([M]18×3)[B]18×1
上式中,[B]18×1为磁场测得的磁通密度,[I]3×1为流经导体的电流,[M]18×3为系数矩阵。
本发明的有益效果在于:
1、不含铁芯,有效避免了磁芯饱和带来的信号畸变的问题;
2、结构简单,体积小,便于安装;
3、属于非接触的测量方法,无需断电后安装;
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为单根倾斜载流导体在三维磁场传感器S产生的磁场示意图;
图2为三相倾斜导体在三维磁场传感器S产生的磁场示意图;
图3为六个三维磁场传感器阵列均匀分布在一个圆环区域;
图4为三维磁场传感器的引脚图;
图5为磁场传感器差分输出信号的后续处理电路;
图6为本发明所提出的电流测量系统;
图7为本发明提出的方法所计算的三相电流幅值误差;
图8为采用二维磁场传感器阵列测量三相倾斜导体电流的幅值误差。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本发明的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
本发明实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本发明的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
本发明提出的方法利用磁场传感器阵列来测量倾斜多导体系统的电流。先分析单根倾斜导体的测量方法。
图1为单根导体穿过磁场传感器阵列的中空区域并在某个三维磁场传感器产生的磁场示意图。图1中,A(xc,yc,0)为倾斜导体与阵列的交点坐标,α,β分别为倾斜导体与z轴的夹角,导体在磁场传感器阵列的投影与x轴方向的夹角。S(xs,ys,0)为三维磁场传感器的坐标。Bs,Bsx,Bsy,Bsz分别为点S处的合成磁场,x,y,z方向的磁场分量。I为流经导体的电流。根据电磁理论,可以推导得到电流与磁场的关系,如式(1)所示。当导体与磁场传感器阵列有一定的倾斜角时,以阵列方式测量得到的磁场含有Z轴方向的分量,如果此时仍用二维的磁场传感器阵列测量导体电流,所得到的误差将远远超过可接受范围,本发明后续会进行详细的误差分析。由公式(1)可知,未知参数有α,β,xc,yc和I,为了求解上述五个参数,必须建立五个独立方程,为了避免多解问题的出现,可以采用两个三维磁场传感器,从而得到由6个方程构成的超定方程组。
图2为三相导体穿过磁场传感器阵列的中空区域,并且在某个三维磁场传感器产生的磁场示意图。根据叠加定理,此时传感器在x,y,z方向的磁场分别为三根导体在对应方向的叠加。由于x,y,z轴方向的计算磁场的表达式不同,因此,由叠加定理可以得到三种类型的方程,如式(2)所示。公式(2)所表述的方程中,含有的未知参数有αii,xci,yci,Ii(i=1,2,3)共15个未知参数,需要15个独立方程来确定上述参数。考虑多解问题的出现,本专利采用6个三维磁场传感器来测量磁场,从而得到18个方程构成的超定方程组。
本文研究的倾斜多导体系统的测量方法,实际是一个非线性优化问题。约束条件包括等式约束条件和不等式约束条件。公式(2)为表征磁场与电流关系的等式约束条件,不等式约束条件由式(6)~(8)决定。上述非线性优化问题采用进化算法来求解,经过种群的进化迭代,所得最优个体使目标函数(式(9)所示)的值最小。待优化的参数有α11,xc1,yc1,I1,Real,I1,Imag22,xc2,yc2,I2,Real,I2,Imag33,xc3,yc3,I3,Real,I3,Imag共18个。其中I1,Real,I1,Imag,I2,Real,I2,Imag,I3,Real,I3,Imag分别为流经导体1,导体2和导体3的电流的实部和虚部。优化问题的维度为18,然而,由公式(2)可知,当磁场经测量确定,并且假设αii,xci,yci(i=1,2,3)已知时,各相导体的电流可直接计算得到。从而,优化问题的维度由18降低到12,即αii,xci,yci(i=1,2,3),从而大大提升进化算法的计算效率。
目前,已有的进化算法很多。为了研究进化算法的类型是否影响此非线性优化问题的求解效率和准确度,本发明对比研究了四种进化算法,即改进粒子群算法,蛙跳算法,帝国进化算法和差分进化算法。由于进化算法的随机特性,每一次进化迭代寻优,未必能得到准确的结果。因此,针对三相导体穿过电流传感器的13种不同情况,分别采用上述四种进化算法来求解非线性优化问题,所得结果如表1所示。
表1四种进化算法求解非线性优化问题结果表
Figure BDA0002525030270000061
Figure BDA0002525030270000071
表1的第一行为三相导体穿过电流传感器不同情况的索引。针对每一种情况,进化算法均运行了50次,并将这50次运行所对应的目标函数值的最大值,最小值和平均值列于表中。同时将上述三种类型的目标函数值的平均值列于表1的最后一列。经仿真研究发现,当目标函数值小于5e-5时,三相电流的误差小于2%,而目标函数值大于5e-5,电流的计算精度较差,因此,目标函数值的大小可以作为衡量算法是否收敛较高精度的依据。引入参数VOF为收敛到较高精度结果的个数,用于表征各个算法的收敛率。参数VOF也列于表1中。仿真使用的笔记本电脑的配置为2.6GHz的因特尔核,6700HQ的CPU以及8GB的内存。需要说明的是,表1中的目标函数值的数量级为1e-5,即实际的目标函数值的大小为表中值乘以1e-5。各个算法的平均运行时间列于表1第一列。由表1可知,帝国进化算法的计算时间较其它三种稍长,但是拥有最小的目标函数值和最大的VOF,因此实际应用中,可以采用帝国竞争算法进行非线性优化问题的求解。
介绍电流测量原理后,接下来详细介绍用于电流测量的电流传感器。详见图3-6所示。
图3中,6个三维磁场传感器焊接在圆环形PCB板上,并且呈圆周均匀分布。Rc,Rs分别为圆环形PCB板的内径和外径。载流导体穿过半径为Rc的中空区域,导体可与PCB垂直,或具有一定的倾斜角。
本发明所采用的是基于穿隧磁阻效应的三维磁场传感器,选用的是多维科技公司生产的TMR2301磁场传感器。其引脚如图4所示,可以同时测量x,y,z轴方向的磁场分量。
本发明所采用的三维磁场传感器为差动输出方式,为了抑制共模干扰信号,提高磁场测量精度,本发明采用差动放大方式,图5为磁场传感器后续的放大电路。运算放大器ad8220采用差分输入方式,将磁场传感器的差分输出信号进行放大,转换为单端信号后,其与运算放大器ad8513组合,又将单端信号转化为差分信号,然后由DAQ以差分方式采集,从而得到磁场信号。
图6为本发明提出的电流测量系统。6个三维磁场传感器均匀分布在圆环形PCB板上,三相导体穿过中空区域,导体电流在周围空间产生的磁场由磁场传感器测量,信号经DAQ采集,传输到上位机,经过逆推算法计算得到各相电流,并将数据可视化。
为了证明本发明的有效性和具有潜在的应用价值,分别采用本发明提出的电流测量方法和基于二维磁场传感器阵列的电流测量方法来计算导体与电流传感器倾斜情况下的电流,并计算各相电流的误差,通过误差的对比分析来说明本发明所提方法的优越性。
仿真研究中的给定参数:
导体1与电流传感器的交点为(-0.2,0.2),α11分别在0~π/4,0~π/3,且变化的步长为π/180;导体2与电流传感器的交点为(0.6,0.4),α22分别为π/6,π/3;
导体3与电流传感器的交点为(0.4,-0.5),α33分别为π/7,π/9;
导体1,2,3所载电流分别为15A,15∠-2π/3A,15∠2π/3A。
采用本发明提出的电流测量方法来计算电流误差的步骤为:
步骤S1:由给定参数计算磁场传感器所在位置的三维磁场,考虑磁场传感器在实际测量中,易受电压偏移和环境噪声磁场的影响,在计算磁场中施加40dB的高斯白噪声,以此来模拟实际的磁场测量值;
步骤S2:将上述磁场值代入逆推算法,从而计算各导体的电流,与电流传感器的交点以及倾斜角;
步骤S3:改变α11的值,重复步骤S1,S2;
步骤S4:以α11为x,y坐标,以电流误差作为z坐标绘制各相电流的误差随α11变化的三维图,如图7所示;
基于二维磁场传感器阵列的电流测量方法实现电流准确测量的前提为导体必须与电流传感器垂直。然而,实际应用中有时难以保证上述的垂直关系。而本文提出的方法很好地解决了导体倾角给电流测量带来较大误差的问题。基于二维磁场传感器阵列的电流测量方法与本发明提出方法的本质区别在于所求解的参数为xci,yci,Ii(i=1,2,3),因此仅需3个二维的磁场传感器测量9个磁场值,构建9个方程即可求解得到上述参数。令本发明提出的逆问题模型中的α,β为零即可得到相应的逆问题算法的约束条件和目标函数。
基于二维磁场传感器阵列的电流测量方法的电流计算步骤为:
步骤S1:由给定参数,利用倾斜导体电流与空间磁场的数学模型计算二维磁场传感器所在位置的x,y方向的磁场,不考虑z方向的磁场,忽略高斯白噪声,直接以计算值来模拟实际的磁场测量值;
步骤S2:将上述磁场值代入基于二维磁场传感器阵列的电流测量方法中的逆推算法,从而计算各导体的电流以及导体与电流传感器的交点;
步骤S3:改变α11的值,重复步骤S1,S2;
步骤S4:以α11为变量,以电流误差作为因变量绘制各相电流的误差随α11变化的三维图,如图8所示;
由图7可知,采用本发明提出的电流测量方法所计算的各相电流误差均在2%以下,即其测量误差基本不受导体倾斜的影响,具有稳定的测量精度;然而,当采用基于二维磁场传感器阵列的电流计算方法时,三相导体电流的误差具有不同的变化趋势,同时误差很大。如图8所示的a相电流的误差趋势图,电流误差可以达到40%,远远超过可接受的误差范围,同样的,b,c相的电流误差同样较大,变化明显。从而证明了本文提出的电流测量方法的有效性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (2)

1.一种用于倾斜多导体系统的电流测量系统,其特征在于:该系统包括电流传感器,数据采集卡和上位机,上述电流传感器包括磁场检测组件、差分放大器,且均集成在圆环形PCB板上;
所述电流传感器呈圆环形结构,多根导体由中空圆孔穿过,导体电流在空间产生的磁场由磁场检测组件测量得到;
所述磁场检测组件包括5个隧穿磁阻效应TMR三维磁场传感器,呈均匀的圆周分布;以电流传感器的中心作为坐标原点,建立三维的笛卡尔坐标系;TMR三维磁场传感器x,y,z轴方向的测量分量分别朝向上述笛卡尔坐标系的三个坐标方向;
所述TMR三维磁场传感器的输出信号为差动输出模式,为抑制共模干扰,后续信号处理电路将传感器输出信号差分放大转化为单端信号后,再以差分方式输出,保证信号都已差分方式传输,从而提升电路的噪声抑制能力;
所述差分放大器包括2个运算放大器,分别记为放大器1和放大器2,磁场传感器的输出信号为差分模式,放大器1以差分方式接收上述输出信号,并以一定增益放大,从而将双端信号转化为单端信号输出,随后,放大器2与放大器1的输出构成差分的输出结构,进而将上述单端信号再转化为双端信号,即差分信号输出;
所述数据采集卡为32路单端或16路差分模拟输入,采样率为100kS/s,采用16路差分模拟输入模式,需2块上述数据采集卡同步采样6个三维TMR磁场传感器的18路差分信号,并将信号显示在上位机系统;
所述上位机系统与数据采集卡进行数据传输,后续利用所采集的磁场数据进行逆问题计算来求解流经各相导体的电流,同时将数据可视化;
测量时,多导体穿过所述中空圆孔;
校准过程包括以下步骤:
步骤S1:根据无限长导体在空间产生的磁场与电流的关系,当单根导体与磁场传感器阵列倾斜时,建立流经导体的电流与空间磁场的数学模型;
步骤S2:由步骤S1所得结果,利用叠加原理,得到多根倾斜导体与空间磁场的数学模型;
步骤S3:由所述磁场检测组件测量空间6个点的三维磁场分量;
步骤S4:利用步骤S3所测磁场值,并结合步骤S2的逆问题模型,采用进化算法逆推各导体的坐标位置,倾斜角和流经导体的电流;
步骤S5:由步骤S4得到的导体坐标位置和倾斜角计算系数矩阵[M],最终完成电流传感器的校准过程;
步骤S6:电流传感器投入实际应用时,将所测磁场值构成的列向量乘以由步骤S5得到的系数矩阵[M]的伪逆,从而得到待测电流值;
所述步骤S1中,计算公式为:
Figure FDA0003725893270000021
式中,Bs,Bsx,Bsy,Bsz分别为空间某一点处的合成磁场,x,y,z方向的磁场分量;(xs,ys)为各个三维磁场传感器的坐标位置,(xc,yc)为导体与磁场传感器阵列所在平面的交点坐标,I为流经导体的电流,m=sin(α)cos(β),n=sin(α)sin(β),p=cos(α),α,β分别为倾斜导体与z轴方向的夹角,倾斜导体在xoy平面的投影与x轴方向的夹角;
所述步骤S2中,计算公式为:
Figure FDA0003725893270000022
由于每个磁场传感器测量得到的是x,y,z轴三个方向的磁场,因此分别对三个方向的磁场进行叠加,从而得到公式(2)所示的三种形式的以X为变量的方程;公式(2)中,αi,βi(i=1,2,3)为三相导体的倾斜角,xci,yci(i=1,2,3)为三相导体的坐标位置,Ii(i=1,2,3)为流经三相导体的电流;Fi(X),Gi(X),Hi(X)的表达式如公式(3);式(3)中的f(i,j)表达式由式(4)表示;与倾斜角相关的参数如式(5)所示;
Figure FDA0003725893270000023
f(i,j)=pj 2(ysi-ycj)2+pj 2(xsi-xcj)2+[mj(ysi-ycj)-nj(xsi-xcj)]2 (4)
Figure FDA0003725893270000031
在所述步骤S1和S2中,得到空间磁场与倾斜多导体电流的数学模型;由公式(2),含有的未知变量数为15,分别为导体倾斜角αi,βi,导体的坐标位置xci,yci以及各相导体的电流Ii(i=1,2,3);采用6个三维磁场传感器,测量得到18个x,y,z轴方向的磁场分量,由叠加定理确定18个方程;采用进化算法进行多目标寻优时,非线性问题所满足的约束条件为
|I1|,|I2|,|I3|<MAXcurrent (6)
Figure FDA0003725893270000032
Figure FDA0003725893270000033
上式中,MAXcurrent,MAXradius分别为导体允许流过的电流,导体所在限制区域的最大半径;
涉及的目标函数为
Figure FDA0003725893270000034
上式中,Bix,Biy,Biz为磁场测量值,
Figure FDA0003725893270000035
为由式(2)计算所得磁场值,要求所计算的电流不能超过允许值,因此添加了一个惩罚项,并将惩罚因子设置为1000;进化算法的目的就是通过群体的进化迭代得到最优个体,使得式(6)的值尽量小。
2.根据权利要求1所述的一种用于倾斜多导体系统的电流测量系统 ,其特征在于:所述步骤S4中,确定几何参数α11,xc1,yc122,xc2,yc233,xc3,yc3后,公式(2)写成矩阵形式,即
Figure FDA0003725893270000036
上式中,[I]为待求的导体电流,[B]为磁场传感器的测量值,[M]为系数矩阵;
电流由下式计算
[I]3×1=pseudo([M]18×3)[B]18×1
上式中,[B]18×1为磁场测得的磁通密度,[I]3×1为流经导体的电流,[M]18×3为系数矩阵。
CN202010501833.1A 2020-06-04 2020-06-04 一种用于倾斜多导体系统的电流测量系统及方法 Active CN111624389B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010501833.1A CN111624389B (zh) 2020-06-04 2020-06-04 一种用于倾斜多导体系统的电流测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010501833.1A CN111624389B (zh) 2020-06-04 2020-06-04 一种用于倾斜多导体系统的电流测量系统及方法

Publications (2)

Publication Number Publication Date
CN111624389A CN111624389A (zh) 2020-09-04
CN111624389B true CN111624389B (zh) 2022-08-09

Family

ID=72259563

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010501833.1A Active CN111624389B (zh) 2020-06-04 2020-06-04 一种用于倾斜多导体系统的电流测量系统及方法

Country Status (1)

Country Link
CN (1) CN111624389B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112595873B (zh) * 2020-09-18 2022-05-17 国网江苏省电力有限公司徐州供电分公司 基于三轴隧道磁阻阵列的电流传感器及其测量方法
CN112213679B (zh) * 2020-10-22 2022-11-08 国网福建省电力有限公司 基于位置信息的磁敏电流互感器估值方法
CN112834805A (zh) * 2021-01-04 2021-05-25 中国电力科学研究院有限公司 具有位置误差校准功能的隧穿磁阻电流传感器及校准方法
CN113049874B (zh) * 2021-03-17 2022-03-08 南方电网数字电网研究院有限公司 电流传感器、电流测量设备、系统、装置和存储介质
CN113325343B (zh) * 2021-05-24 2022-05-03 合肥工业大学 一种单轴隧道磁阻电流测量方法
CN113848367A (zh) * 2021-09-22 2021-12-28 温州大学乐清工业研究院 一种非接触式电压测量的自适应动态补偿方法及装置
CN113820532B (zh) * 2021-09-23 2022-04-15 南方电网数字电网研究院有限公司 非接触式双芯电缆电流测量方法和装置
CN114441833B (zh) * 2022-01-24 2023-03-31 南方电网数字电网研究院有限公司 电流测量方法、装置、计算机设备、存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10158836A1 (de) * 2001-11-30 2003-06-18 Infineon Technologies Ag Verfahren und Vorrichtung zur Kalibrierung eines Sonsorsystems
CN102012446A (zh) * 2009-09-04 2011-04-13 罗伯特·博世有限公司 电流传感器、电流测量模块以及用于电流测量的方法
CN108594151A (zh) * 2018-06-02 2018-09-28 福州大学 无磁芯电流传感器的位置误差补偿法
CN110546519A (zh) * 2017-04-07 2019-12-06 西门子股份公司 电流测量方法和电流测量装置
CN111157777A (zh) * 2020-01-14 2020-05-15 清华大学 一种双磁芯测量差分泄漏电流传感器设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10158836A1 (de) * 2001-11-30 2003-06-18 Infineon Technologies Ag Verfahren und Vorrichtung zur Kalibrierung eines Sonsorsystems
CN102012446A (zh) * 2009-09-04 2011-04-13 罗伯特·博世有限公司 电流传感器、电流测量模块以及用于电流测量的方法
CN110546519A (zh) * 2017-04-07 2019-12-06 西门子股份公司 电流测量方法和电流测量装置
CN108594151A (zh) * 2018-06-02 2018-09-28 福州大学 无磁芯电流传感器的位置误差补偿法
CN111157777A (zh) * 2020-01-14 2020-05-15 清华大学 一种双磁芯测量差分泄漏电流传感器设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维巨磁阻磁场传感器的电路设计和实验标定;林楚 等;《机电技术》;20161231;第62-65页 *
利用正交型锁相放大器实现三维磁场微弱信号检测;康中尉 等;《传感器技术》;20041231;第23卷(第12期);第69-72页 *

Also Published As

Publication number Publication date
CN111624389A (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
CN111624389B (zh) 一种用于倾斜多导体系统的电流测量系统及方法
CN110007350B (zh) 一种磁探测方法盲区的分析方法
US20200326389A1 (en) Method for Making a Magnetic Gradiometer with High Detection Accuracy and Success Rate
CN107817457B (zh) 一种地磁梯度张量测量阵列的设计方法
Bazzocchi et al. Interference rejection algorithm for current measurement using magnetic sensor arrays
CN110146839A (zh) 一种移动平台磁梯度张量系统校正方法
CN109725360A (zh) 基于磁梯度张量不变量的单点定位方法
CN111522835B (zh) 一种基于数据库特征匹配的多磁性目标位置探测方法
CN111399066B (zh) 一种基于正交基函数处理标量磁异常梯度信号的方法
CN110852025A (zh) 一种基于超收敛插值逼近的三维电磁慢扩散数值模拟方法
CN116539940A (zh) 一种电流测量装置和方法
CN109725361B (zh) 基于磁梯度张量不变量的一种磁性目标定位方法
CN113609749B (zh) 适用于多场景的基于磁场信号的电流计算方法
CN109633540B (zh) 一种磁源的实时定位系统及实时定位方法
Yousefi et al. A fast and robust magnetic localization technique based on elimination of the orientation variables from the optimization
CN115979245B (zh) 一种自校准估计的磁感应网络定位方法
CN110646707B (zh) 一种基于微量处理法与粒子群算法的接地网故障诊断方法
Primin et al. Inverse problem solution algorithms in magnetocardiography: new analytical approaches and some results
Laitinen et al. A first/third-order probe correction technique for spherical near-field antenna measurements using three probe orientations
CN111596119A (zh) 一种基于导线磁场的检测方法及装置
CN109188320B (zh) 一种基于磁阻效应的流场成像系统及成像方法
CN116184512A (zh) 一种基于磁矩不变量的矢量传感器网络实时定位方法及系统
CN114660511A (zh) 一种磁传感器阵列的标定方法
JPH04319334A (ja) 生体磁気イメージング装置
Skvortsov et al. Estimation of localization of point sources from a printed circuit board in the near field

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