CN102551711B - 用超导量子干涉器测量磁场并定位磁场电流源的方法 - Google Patents
用超导量子干涉器测量磁场并定位磁场电流源的方法 Download PDFInfo
- Publication number
- CN102551711B CN102551711B CN 201010580178 CN201010580178A CN102551711B CN 102551711 B CN102551711 B CN 102551711B CN 201010580178 CN201010580178 CN 201010580178 CN 201010580178 A CN201010580178 A CN 201010580178A CN 102551711 B CN102551711 B CN 102551711B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- field
- maximum
- minimum
- interference device
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明涉及用超导量子干涉器测量磁场并定位磁场电流源的方法,该方法包括以下步骤:1)利用超导量子干涉器采集少量测点上的磁场信号;2)根据测量到的磁场信号求出每个时刻的磁场极大值和磁场极小值;3)计算机对步骤2)中求得的每个时刻的磁场极大值和磁场极小值处理后求出磁场源的位置和强度参数。与现有技术相比,本发明具有方法简单、有效,且可以大大节省硬件成本等优点。
Description
技术领域
本发明涉及磁场电流源定位的方法,尤其是涉及用超导量子干涉器测量磁场并定位磁场电流源的方法。
背景技术
1970-1976年D.Cohen等首次提出了在实验室里用超导量子干涉器(superconducting quantum interference device,SQUID)测量心脏磁场信号,以及定位一个等效电流偶极子源的方法。近40年来,磁场检测及源定位的研究引起各国高度的重视。应用范围不断扩大。为了提高源定位的精度,用SQUID测量磁场的通道数不断增加,这也导致设备成本增加。一般商用心脏磁场测量***需要采用6X 6排列的36-64个通道SQUID同步测量磁场信号,不仅硬件***复杂,而且磁场源定位逆问题的计算方法和计算软件也面临诸多问题。因此,一种可通过少量磁场测点进行磁场电流源准确定位的方法属于新的发明。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供方法简单、高效,且可以节省磁场测量***的硬件成本,简化相关算法和软件,得到准确的源定位结果的用超导量子干涉器测量磁场并定位磁场电流源的方法。
本发明的目的可以通过以下技术方案实现:用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,该方法包括以下步骤:1)利用超导量子干涉器采集少量测点上的磁场信号;2)根据测量到的磁场信号求出每个时刻的磁场极大值和磁场极小值;3)计算机对步骤2)中求得的每个时刻的磁场极大值和磁场极小值进行处理后求出磁场电流源的位置和强度参数。
所述的步骤1)中磁场信号需要测量的点数可以根据磁场电流源位置固定与否来确定,对于位置固定的磁场源,通常只需要两个测量点;如果源位置周期性的变化,则需要采集几个相邻检测点上的磁场信号,通过插值求出测量磁场的极大值和极小值。
所述的步骤1)、步骤2)中的磁场信号包括每个时刻的磁场值和该磁场值对应的X-Y坐标。
所述的步骤3)具体为:
31)首先利用公式
求出磁场电流源的深度d;
33)利用公式
求出磁场强度Q和磁场电流源的坐标(x0,y0),其中,Qx和Qy是该磁场电流源在X和Y坐标的强度分量。
所述的步骤2)中,当根据磁场信号可以求出磁场极大值和磁场最小值时,就用该磁场最小值近似替代磁场的极小值。
所述的步骤2)中,当根据磁场信号只能求得磁场最大值和磁场极小值时,就用该磁场的最大值近似替代磁场极大值。
因为本发明采用以上技术方案,即利用超导量子干涉器采集2~8个测量点上的磁场信号进行磁场电流源定位,所以,定位方法简单、有效,与现有技术采用36或64多通道的SQUID阵列检测方法相比,可以大大节省硬件成本。
附图说明
图1为利用两个SQUID测量点上的磁场信号定位一个磁场电流源的示意图;
图2为测量平面上的测点分布图。
图1中1表示一个电流偶极子源。
具体实施方式
下面结合附图1和具体实施例对本发明进行详细说明。
实施例
图1给定坐标系中有A和B两个测量点。A点表示磁场的极大值所在的位置,B点表示磁场极小值的位置。X-Y坐标组成的测量平面的下方,位于深度d处有一个产生被测量磁场的电流偶极子源。测量平面上测量点A与该电流偶极子的连线的夹角为θ,测量平面上A、B测量点的连线与X轴的夹角为
用该方法进行测量磁场的源定位时,基本原理如图1所示。
在Maxwell方程准静态近似条件下,根据Biot-Savart定律,空间中一个电流偶极子产生的磁感应强度为:
其中,Qx和Qy是磁场强度Q在X和Y方向的分量。产生被测磁场的电流偶极子源的位置坐标是(x0,y0,z0),电流偶极子源的深度d=-z0。
当测量磁场的极大值和极小值相等时,产生磁场的电流偶极子源的深度d可以根据图1用公式计算。其中,D是测量平面上坐标(x1,y1)处磁场极大值Bzmax和坐标(x,y)处磁场极小值Bzmin之间的距离。该电流源的位置(x0,y0)可以用中点法计算。但是,实际测量到的磁场往往不是关于测量平面对称分布的,而且磁场与产生磁场的电流源的参数之间呈非线性关系,在求解此类逆问题时解不唯一。为此,我们提出了一种新的源定位方法。
图1给出了检测平面上两个测量点上的磁场极大值和极小值与电流源位置的几何关系:
其中,
磁场极大值Bzmax和磁场极小值Bzmin可以通过测量和计算得到。测量平面上,测点(x1,y1)处检测到的Z方向的磁场极大值:
位于(x,y)的磁场极小值:
当已知测量磁场的极大值Bzmax和极小值Bzmin时,公式(8)和(9)中的Q和d可以由这两个方程解出。从而由(5)、(6)和(7)式可以求出电流偶极子在X和Y方向上的电流强度分量,并由(2)和(3)式计算出电流源的平面坐标。当磁场的极大值和极小值没有被检测到时,可以用检测到的磁场最大值和最小值近似替代。
实施步骤如下:
1.将两个SQUID探头分别对准可以检测到磁场极大值Bzmax和极小值Bzmin的测量位置,并确定这两个测点在检测平面上的坐标,即(x1,y1)处是磁场极大值Bzmax,(x,y)是测量到的磁场极小值Bzmin的位置。
4.用公式(2),(3)计算该电流偶极子的X-Y坐标(x0,y0)。
5.根据公式(8)计算出Q。
6.根据公式(7)求出Qx和Qy分量。
如果没有磁场峰值和谷值随时间变化的情况,通常只需要两个测点就可以按照上述方法完成源定位的计算。如果测量磁场的峰值和谷值有变化,可以通过相邻测点上的信号取值,再通过插值求出磁场信号的极大值和极小值,然后进行电流源参数的计算。每个测点的间距可以是3.5-4cm。
采集磁场信号时,如图2所示,先根据经验确定磁场极大值Bzmax所在的位置,磁场极小值Bzmin的位置通常是在磁场极大值的右侧,需要根据经验搜索。磁场的峰值和谷值的测量位置确定后,在这两个测点上同步检测心脏磁场的信号。
Claims (5)
1.用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,该方法包括以下步骤:
1)利用超导量子干涉器采集少量测点上的磁场信号;
2)根据测量到的磁场信号求出每个时刻的磁场极大值和磁场极小值;
3)计算机对步骤2)中求得的每个时刻的磁场极大值和磁场极小值处理后,求出磁场电流源的位置和强度参数;
所述步骤3)的具体为:
31)首先利用公式
求出磁场电流源的深度d;
33)利用公式
求出磁场强度Q和磁场电流源的坐标(x0,y0),其中,Qx和Qy是该磁场电流源在X和Y坐标的强度分量。
2.根据权利要求1所述的用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,所述的步骤1)中磁场信号需要测量的点数根据磁场电流源的位置固定与否来确定,对于位置固定的磁场源,只需要两个测量点;如果源位置周期性的变化,则需要采集几个相邻测量点上的磁场信号,通过插值求出测量磁场的极大值和极小值。
3.根据权利要求1所述的用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,所述的步骤1)和步骤2)中的磁场信号包括每个时刻的磁场值和该磁场值对应的X-Y坐标。
4.根据权利要求1所述的用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,所述的步骤2)中,当根据磁场信号可以得到磁场极大值和磁场最小值时,就用该磁场最小值近似替代磁场的极小值。
5.根据权利要求1所述的用超导量子干涉器测量磁场并定位磁场电流源的方法,其特征在于,所述的步骤2)中,当根据磁场信号只能求得磁场最大值和磁场极小值时,就用该磁场的最大值近似替代磁场极大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010580178 CN102551711B (zh) | 2010-12-09 | 2010-12-09 | 用超导量子干涉器测量磁场并定位磁场电流源的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010580178 CN102551711B (zh) | 2010-12-09 | 2010-12-09 | 用超导量子干涉器测量磁场并定位磁场电流源的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102551711A CN102551711A (zh) | 2012-07-11 |
CN102551711B true CN102551711B (zh) | 2013-10-16 |
Family
ID=46399289
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010580178 Expired - Fee Related CN102551711B (zh) | 2010-12-09 | 2010-12-09 | 用超导量子干涉器测量磁场并定位磁场电流源的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102551711B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092033B (zh) * | 2017-04-19 | 2019-06-11 | 中南大学 | 一种用于调控测量信号强度的方法 |
CN109375368B (zh) * | 2018-10-24 | 2021-04-30 | 泉州师范学院 | 一种基于空间偶极子阵的三维多焦斑阵列的产生方法 |
CN110850483B (zh) * | 2019-11-12 | 2021-09-14 | 中国船舶重工集团公司第七一九研究所 | 一种基于电场电极阵列布置的水中目标探测和定位方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1552285A (zh) * | 2003-05-26 | 2004-12-08 | 中国科学院物理研究所 | 一种心磁波组图及其分析方法 |
US7521708B1 (en) * | 2004-12-29 | 2009-04-21 | The United States Of America As Represented By The Secretary Of The Navy | High sensitivity ring-SQUID magnetic sensor |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4000343B2 (ja) * | 1999-05-20 | 2007-10-31 | 株式会社島津製作所 | 生体活動電流源推定装置 |
-
2010
- 2010-12-09 CN CN 201010580178 patent/CN102551711B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1552285A (zh) * | 2003-05-26 | 2004-12-08 | 中国科学院物理研究所 | 一种心磁波组图及其分析方法 |
US7521708B1 (en) * | 2004-12-29 | 2009-04-21 | The United States Of America As Represented By The Secretary Of The Navy | High sensitivity ring-SQUID magnetic sensor |
Non-Patent Citations (5)
Title |
---|
JP特开2000-325323A 2000.11.28 |
基于单电流偶极子的心脏磁场反演计算;张磊等;《现代医学仪器》;20051231;第2005卷(第4期);31-34 * |
张磊等.基于单电流偶极子的心脏磁场反演计算.《现代医学仪器》.2005,第2005卷(第4期),31-34. |
浅谈心脏磁场的检测及心脏源模型建模;石明伟等;《硅谷》;20111231;44-45 * |
石明伟等.浅谈心脏磁场的检测及心脏源模型建模.《硅谷》.2011,44-45. |
Also Published As
Publication number | Publication date |
---|---|
CN102551711A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597555B (zh) | 一种瞬变电磁法的接地网腐蚀程度评价方法 | |
CN107328849A (zh) | 薄化检测***和薄化检测方法 | |
CN106772639B (zh) | 地下铁质管线埋深磁偶极子构造法优化反演 | |
CN112525201B (zh) | 一种基于电磁场特征多信息融合的水下目标跟踪方法 | |
CN105043381A (zh) | 一种基于磁钉的定位方法 | |
CN108151766A (zh) | 磁钉的定位方法、磁钉定位导航误差修正方法及定位装置 | |
CN108152687B (zh) | 一种利用超声反演的电力变压器局部放电定位方法 | |
CN102621582A (zh) | 空间矢量法地下金属管线位置探测方法及装置 | |
CN102551711B (zh) | 用超导量子干涉器测量磁场并定位磁场电流源的方法 | |
CN114235962B (zh) | 一种面向各向异性结构的超声导波成像方法及*** | |
CN110161572A (zh) | 一种接地网水平拓扑结构快速无损探测***与方法 | |
CN103292773A (zh) | 一种基于最小区域的对称度误差评定方法 | |
CN106483583B (zh) | 电子标签的定位***、方法及电子标签探测仪 | |
CN102809687A (zh) | 一种交流电频率的数字化测量方法 | |
CN111398752A (zh) | 一种基于多探测器信息融合的电力变压器局部放电定位装置及方法 | |
CN107238364B (zh) | 接触式测量测针球头半径精密补偿方法 | |
CN112344925B (zh) | 基于归一化磁矩的水下目标相关性检测定位方法 | |
CN102928666B (zh) | 一种交流电的相位差的数字化测量方法 | |
CN103575981A (zh) | 一种交流电频率的精确测量方法 | |
CN204807074U (zh) | 一种大型盖类产品的平面测量装置 | |
CN114236620A (zh) | 一种金属管线定位方法及装置 | |
CN203259109U (zh) | 一种柔性汽车玻璃检具 | |
CN110687609B (zh) | 一种埋地管线交直流杂散电流检测中平行管线干扰消除方法 | |
CN105372710B (zh) | 一种导线探测方法及*** | |
CN208704763U (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131016 Termination date: 20161209 |