Estimation method for differential magnetic gradient tensor measurement limits

差分磁梯度张量测量极限估计 Estimation method for differential magnetic gradient tensor measurement limits first-author 李青竹（1993-），男，四川绵阳人，博士研究生，2016年于西南交通大学获得学士学位，2018年于陆军工程大学获得硕士学位，主要从事磁异常探测，弱磁测试信号处理等方面的研究。E-mail： laznlqz666@qq.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777817&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777830&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777834&type= 22.01333237 32.00399017 李青竹 （1993-），男，四川绵阳人，博士研究生，2016年于西南交通大学获得学士学位，2018年于陆军工程大学获得硕士学位，主要从事磁异常探测，弱磁测试信号处理等方面的研究。E-mail： laznlqz666@qq.com laznlqz666@qq.com 石志勇（1965-），男，四川眉山人，教授，博士生导师，1997年、2000年于军械工程学院分别获得硕士和博士学位，主要从事弱磁测试和惯性导航技术的研究。E-mail： szy651010@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777836&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777849&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777839&type= 22.01332855 32.00399017 石志勇 （1965-），男，四川眉山人，教授，博士生导师，1997年、2000年于军械工程学院分别获得硕士和博士学位，主要从事弱磁测试和惯性导航技术的研究。E-mail： szy651010@163.com szy651010@163.com corresp E-mail： lgdsxq@163.com E-mail： lgdsxq@163.com 李志宁（1972-），男，河北石家庄人，副教授，博士生导师，1999年于军械工程学院获得硕士学位，2007年于清华大学获得博士学位，主要从事测试技术与信号处理的研究。E-mail： lgdsxq@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777852&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777856&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777854&type= 22.01334000 32.00399780 李志宁 （1972-），男，河北石家庄人，副教授，博士生导师，1999年于军械工程学院获得硕士学位，2007年于清华大学获得博士学位，主要从事测试技术与信号处理的研究。E-mail： lgdsxq@163.com lgdsxq@163.com 陆军工程大学 车辆与电气工程系，河北 石家庄 050003 陆军工程大学 车辆与电气工程系，河北 石家庄 050003 Department of Vehicle and Electrical Engineering， Army Engineering University of PLA， Shijiazhuang 050003， China Department of Vehicle and Electrical Engineering， Army Engineering University of PLA， Shijiazhuang 050003， China 为了研究磁梯度张量系统（Magnetic Gradient Tensor System， MGTS）的理论探测极限，利用磁偶极子正演公式、张量矩阵特征方程和张量不变量推导了差分磁梯度张量测量范围公式，提出了MGTS理论探测极限估计方法。利用磁梯度张量空间衍生不变关系定位方法估计目标体磁矩，并利用差分计算原理在背景场静态采样中估计MGTS的张量分量测量准确度。利用提出的方法估计现场MGTS的理论探测极限。滑动磁性物体远离MGTS并进行连续采样，记录张量信号实际可观测范围并对估计结果进行验证。结果表明，MGTS的理论探测极限与基线距离 d 、测量准确度 q 、目标磁矩强度 M 、观测点与目标磁矩间夹角 θ 等有关； d 越长， q 越高， M 越大，MGTS理论可探测距离越远； d 越长，观测距离越近，MGTS测量的理论误差越大；探测距离在MGTS平行于磁矩方向时最大，在垂直于磁矩方向时急剧减小。实验表明，平面十字形MGTS针对4个典型磁铁的探测极限估计准确度为±0.4 m。 To study the theoretical detection limit of the magnetic gradient tensor system （MGTS）， the differential tensor measurement range formula was derived with the magnetic dipole forward equation， the tensor matrix characteristic equation and the tensor invariant， and the MGTS theoretical detection limit estimation method was proposed. First， the target magnetic moment was estimated with the magnetic gradient tensor derivative invariant relations， and the difference calculation principle was used to estimate the measurement accuracy of the tensor component of the MGTS in the static sampling of the background field； then， the theoretical detection limit of on-site MGTS was estimated with the proposed method； and finally， the magnetic object was slided away from the MGTS and continuous sampling was performed， recording the actual tensor signal observable range， and verifying the estimation result of detection limit. The results show that the theoretical detection limit of MGTS is related to the baseline distance d ， the measurement accuracy q ， the target magnetic moment intensity M ， the angle θ between the target magnetic moment and the observation point， etc.； The longer d ， the higher q ， and the larger M ， the farther detect distance； The longer d and the closer observation， the greater the theoretical differential measurement error； The detection range of MGTS reaches the maximum in the direction parallel to the magnetic moment， but decreases sharply in the direction perpendicular. The experiments show that the estimation accuracy of detection limit of the planar-cross MGTS for four typical magnets is ± 0.5 m. 磁梯度张量 差分测量系统 探测极限 目标定位 magnetic gradient tensor differential measurement system detection limits target positioning 1 引　言 1 磁异常探测技术可应用于地下小尺度磁目标的定位与识别 ［ 1 1 - 2 2 1-2 ］ 。与磁场矢量和总磁场强度（Total Magnetic Intensity， TMI）相比，磁梯度张量（Magnetic Gradient Tensor， MGT）可以提供更丰富的目标体方位信息，抗干扰能力强，能够更好地适应复杂的测量环境 ［ 3 3 ］ 。MGT探测具有广泛的应用前景，如航空磁测、矿产勘探、未爆弹药搜索和排雷等 ［ 4 4 - 5 5 4-5 ］ 。 直接测量磁场梯度本质上很困难，但可利用矢量传感器在基线两端的读数差异近似估计磁标势的二次偏微分 ［ 6 6 ］ 。目前，常以磁传感器阵列形式构建磁梯度张量系统（Magnetic Gradient Tensor System， MGTS），并实现张量差分测量。MGTS主要分为两类：（1）基于超导效应 ［ 7 7 ］ ，这类系统由具有极高灵敏度和较小量程的超导量子干涉装置（Superconducting Quantum Interference Devices， SQUID）组成，其制造成本高，对测量环境要求严格，适用于生物磁检测、金属无损检测、航空磁测等灵敏度要求高但磁异常较弱的工况；（2）基于磁通门法 ［ 8 8 - 9 9 8-9 ］ ，这类系统由多个磁通门传感器组成，利于批量生产和制造，成本较低，安装要求较简单。目前，最先进的磁通门探头的灵敏度噪声可以达到6 pTrms/ <math id="M1"><mroot><mrow><mi mathvariant="normal">H</mi><mi mathvariant="normal">z</mi></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777860&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777858&type= 7.70466709 4.06400013 @1 Hz ［ 10 10 ］ ，量程是SQUID的数千倍 ［ 11 11 - 12 12 11-12 ］ 。 近年来，国内外研究团队搭建了各类磁通门法MGTS，包括直角四面体、正四面体、正方形、十字形、三角形等结构 ［ 13 13 - 16 16 13-16 ］ ，但针对特定结构MGTS探测极限的研究却鲜有报道。实测经验表明，差分方法测量磁梯度张量时系统的探测极限常受到结构误差、噪声、磁源强弱和探测方位等因素的影响 ［ 17 17 ］ 。为了定性且定量研究MGTS的理论探测极限，本文利用磁偶极子正演方程、张量矩阵特征方程和张量不变量联合推导出在系统基线距离、传感器测量精度、系统观测方位和磁源强度等相关参数约束下MGTS空间理论探测范围的计算公式。针对单一磁源，张量衍生不变关系定位方法 ［ 18 18 ］ 利用MGTS实现单点精确定位并准确估计出磁源的磁矩强度，使实测MGTS的理论探测极限成为可能。 2 MGT测量理论和张量不变量 1 2.1　差分MGT测量原理 2 MGT为3个正交方向的磁场矢量空间变化率 ［ 3 3 ］ ，共9个元素，可表示为： <math display="block" id="M2"><mi mathvariant="bold-italic">G</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>y</mi><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>y</mi><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>z</mi><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>z</mi><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>φ</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced><mtext> </mtext><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>y</mi></mrow></mfrac></mtd><mtd><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>z</mi></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>x</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>x</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>y</mi></mrow></msub></mtd><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777867&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777864&type= 140.54666138 32.51200104 ， 其中： B 为磁场强度矢量， G 表示MGT矩阵， φ m 是磁标势， B m （ m = x ， y ， z ）表示 B 的正交分量， B ij （ i ， j = x ， y ， z ）表示MGT分量。在没有电流的静态磁环境中，麦克斯韦方程约束下磁场的散度和旋度为0，有 <math id="M3"><mo>∇</mo><mo>·</mo><mi mathvariant="bold-italic">B</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777844&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777870&type= 13.12333298 3.13266683 ， <math id="M4"><mo>∇</mo><mo>×</mo><mi mathvariant="bold-italic">B</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777879&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777847&type= 16.93333435 3.13266683 ，因此 G 是对称且无迹的。故 G 的9个元素可以用5个独立的分量来表示，即 B xx ， B xy ， B xz ， B yy 和 B yz 。 MGT是磁标量势 φ m 的二阶偏微分，也是磁场矢量 B 的偏微分，难以直接测量。事实上，常利用跨测量基线的分量读数差异来代替磁矢量偏微分，以近似MGT分量，如： <math display="block" id="M5"><msub><mrow><mi>B</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>j</mi></mrow></mfrac><mo>≈</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><msub><mrow><mi>B</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">Δ</mi><msub><mrow><mi>d</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfrac><mtext>  </mtext><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777881&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777873&type= 48.17533112 8.89000034 ， 其中： <math id="M6"><mi mathvariant="normal">Δ</mi><msub><mrow><mi>B</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777877&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777875&type= 5.75733376 3.72533321 是两个相邻磁传感器测量的 i 分量读数差异， <math id="M7"><mi mathvariant="normal">Δ</mi><msub><mrow><mi>d</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777888&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777886&type= 5.33400011 3.72533321 是两个磁传感器在 j 方向上的距离，定义为基线距离。 2.2　MGT不变量及衍生不变关系 2 MGT不变量是当观测点固定不随测量系统不同方向而改变的量或对应关系。一些基本的MGT 不变量包括矩阵迹、特征值、特征方程系数和 Frobenius范数等 ［ 19 19 ］ 。MGT矩阵是实对称矩阵，可对角化。令 λ 1 ， λ 2 和 λ 3 表示 G 的特征值，满足特征方程 <math id="M8"><msup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>-</mo><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>λ</mi><mo>-</mo><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777897&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777895&type= 36.23733139 4.06400013 ，其中： <math display="block" id="M9"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mi>t</mi><mi>r</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">G</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub><mo>=</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub><mo>+</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub><mo>-</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>+</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mi mathvariant="normal">d</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">t</mi><mi mathvariant="bold-italic">G</mi><mo>=</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow></msub><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow></msub><mo>-</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msub><mo>=</mo><msub><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mi mathvariant="bold-italic">G</mi></mrow></mfenced></mrow><mrow><mi>F</mi></mrow></msub><mo>=</mo><mroot><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>z</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><mn mathvariant="normal">2</mn><mfenced separators="|"><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfenced></mrow><mrow/></mroot></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777908&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777898&type= 118.02532196 24.80733490 ， 其中： I 0 是 G 的迹， I 1 是关于 G 对称且无迹的量， I 2 是矩阵 G 的行列式， C T 是张量收缩。给定一个磁偶极子，磁矩为 m =（ m x ， m y ， m z ） ［ 20 20 ］ ，偶极子到观测点的位置矢量为 r =（ x ， y ， z ）。设 M =|| m ||， r =|| r ||， r 为观测距离，则磁场矢量和5个独立MGT分量的正演方程为 ［ 19 19 - 20 20 19-20 ］ ： <math display="block" id="M10"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfrac><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi></mrow></mfrac><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">3</mn><mfenced separators="|"><mrow><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow></mfenced><mi mathvariant="bold-italic">r</mi></mrow><mrow><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mi mathvariant="bold-italic">m</mi></mrow><mrow><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac></mrow></mfenced><mo>=</mo><mfrac><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">3</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd><mtd><mn mathvariant="normal">3</mn><mi>x</mi><mi>y</mi></mtd><mtd><mn mathvariant="normal">3</mn><mi>x</mi><mi>z</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">3</mn><mi>x</mi><mi>y</mi></mtd><mtd><mn mathvariant="normal">3</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd><mtd><mn mathvariant="normal">3</mn><mi>y</mi><mi>z</mi></mtd></mtr><mtr><mtd><mn mathvariant="normal">3</mn><mi>x</mi><mi>z</mi></mtd><mtd><mn mathvariant="normal">3</mn><mi>y</mi><mi>z</mi></mtd><mtd><mn mathvariant="normal">3</mn><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd></mtr></mtable></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777904&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777909&type= 112.86066437 16.17133331 ， <math display="block" id="M11"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>z</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><mi>μ</mi></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">7</mn></mrow></msup></mrow></mfrac><mo>×</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">3</mn><mi>x</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mtd><mtd><mi>y</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>y</mi></mtd><mtd><mi>z</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>z</mi></mtd></mtr><mtr><mtd><mi>y</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>y</mi></mtd><mtd><mi>x</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd><mtd><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><mi>y</mi><mi>z</mi></mtd></mtr><mtr><mtd><mi>z</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>z</mi></mtd><mtd><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><mi>y</mi><mi>z</mi></mtd><mtd><mi>x</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd></mtr><mtr><mtd><mi>x</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd><mtd><mn mathvariant="normal">3</mn><mi>y</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mtd><mtd><mi>z</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>z</mi></mtd></mtr><mtr><mtd><mo>-</mo><mn mathvariant="normal">5</mn><mi>x</mi><mi>y</mi><mi>z</mi></mtd><mtd><mi>z</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>z</mi></mtd><mtd><mi>y</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>-</mo><mn mathvariant="normal">5</mn><mi>y</mi><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd></mtr></mtable></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>x</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>m</mi></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>⇒</mo><mi mathvariant="bold-italic">I</mi><mo>=</mo><mi>C</mi><mi>H</mi><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777923&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777911&type= 112.52200317 27.94000053 ， 其中： μ 为介质的磁导率，空气中 μ ≈ μ 0 =4π×10 -7 N•A -2 ， μ 0 为真空磁导率。 空间衍生不变关系如下 ［ 18 18 ］ ： <math display="block" id="M12"><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi><mo>=</mo><mover accent="true"><mi mathvariant="bold-italic">m</mi><mo>¯</mo></mover><mo>⋅</mo><mover accent="true"><mi mathvariant="bold-italic">r</mi><mo>¯</mo></mover><mo>=</mo><mfrac><mrow><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mi mathvariant="bold-italic">m</mi></mrow></mfenced><mo>⋅</mo><mfenced open="‖" close="‖" separators="|"><mrow><mi mathvariant="bold-italic">r</mi></mrow></mfenced></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mroot><mrow><mo>-</mo><msubsup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>-</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow/></mroot></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><mi>u</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777928&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777914&type= 73.23667145 10.32933331 ， （6） 其中： θ 是 r 和 m 之间的角度， <math id="M13"><mn mathvariant="normal">0</mn><mo>≤</mo><mi>θ</mi><mo>≤</mo><mn mathvariant="normal">180</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777934&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777931&type= 19.98133469 3.30200005 ， u 是归一化源强度（Normalized Source Strength， NSS） ［ 19 19 ］ 。 由 式（3） 式（3） 、 式（4） 式（4） 、式（6）分别得到MGT矩阵的特征值 <math id="M14"><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777938&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777916&type= 2.87866688 3.72533321 ， <math id="M15"><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777917&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777941&type= 2.87866688 3.72533321 和 <math id="M16"><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777947&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777918&type= 2.87866688 3.72533321 ，即有： <math display="block" id="M17"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow><mrow><mn mathvariant="normal">8</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mfenced separators="|"><mrow><mo>-</mo><mroot><mrow><mn mathvariant="normal">5</mn><mo>+</mo><mn mathvariant="normal">4</mn><mfrac><mrow><msup><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow><mrow/></mroot><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow><mrow><mn mathvariant="normal">8</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mfenced separators="|"><mrow><mroot><mrow><mn mathvariant="normal">5</mn><mo>+</mo><mn mathvariant="normal">4</mn><mfrac><mrow><msup><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mfenced separators="|"><mrow><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow><mrow/></mroot><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mtd></mtr><mtr><mtd><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mi mathvariant="bold-italic">m</mi><mo>⋅</mo><mi mathvariant="bold-italic">r</mi><mtext> </mtext></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777951&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777949&type= 61.38333893 37.16867065 . 则NSS可以表示为： <math display="block" id="M18"><mi>u</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>M</mi></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777955&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777952&type= 16.34066582 8.21266651 ， 其中： <math id="M19"><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced><mo>≥</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777971&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777956&type= 15.15533447 4.82600021 ， <math id="M20"><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced><mo>≥</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777962&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777972&type= 15.15533447 4.82600021 和 <math id="M21"><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>≥</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>≥</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777966&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777974&type= 18.37266731 4.14866638 。 式（8） 式（8） 可为MGTS的理论探测极限估计提供条件约束。 3 MGTS的理论探测极限 1 相关研究表明，相较于正四面体、直角四面体、三角形和正方形等结构，平面十字形MGTS不仅结构简单、传感器安装方便、对准精度高，其结构误差较小且同基线距离下的最近观测范围最大 ［ 17 17 ］ 。 图1 图1 为平面十字形MGTS的结构设计，由4个磁通门传感器和十字架构成。设 d 为基线距离， <math id="M22"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>b</mi></mrow><mrow><mi>a</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777977&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777983&type= 3.64066648 3.97933316 （ a =1，2，3，4； b = x ， y ， z ）表示传感器 a 在 b 轴方向的读数。结合 式（2） 式（2） ，平面十字形MGTS的张量测量矩阵 G m 为： <math display="block" id="M23"><msub><mrow><mi mathvariant="bold-italic">G</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>=</mo><mtext>  </mtext><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>d</mi></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><mi>z</mi></mrow></msub></mtd><mtd><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><mo stretchy="false">(</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>+</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778001&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777990&type= 88.64600372 14.30866718 ， 式中： B yx 和 B xy 的测量方式存在差异（事实上二者实测值几乎相同），故 G m 共有6个独立分量。 令 n 为一单位空间向量， <math id="M24"><mrow><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">B</mi><mo stretchy="false">)</mo></mrow><mo>/</mo><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">n</mi><mo stretchy="false">)</mo></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778004&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777994&type= 16.67933464 3.97933316 为 B 对 n 的方向导数。用MGT矩阵 G 的特征向量 v 1 ， v 2 和 v 3 表示 n ，有 n = a v 1 + b v 2 + c v 3 且 a 2 + b 2 + c 2 =1，得到： <math display="block" id="M25"><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">B</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">n</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><mi>G</mi><mi mathvariant="bold-italic">n</mi></mrow></mfenced><mo>=</mo><mo stretchy="false">(</mo><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msubsup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msubsup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msup><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msubsup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777996&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778005&type= 59.68999863 9.82133293 . 定义 q 为传感器的测量准确度（单位：±nT）， Q 为MGTS的测量准确度（单位：±nT/m），结合 式（9） 式（9） 有 q = Qd /2。注意， Q 和 q 不仅反映了传感器本底噪声引起的测量误差，还包括传感器系统误差、磁干扰误差和未对准误差引起的不确定性 ［ 8 8 ］ 。将考虑全部实测误差因素时的最大测量偏差作为准确度 Q 和 q 的取值范围。 在观测点 r 0 =（ x 0 ， y 0 ， z 0 ）处捕捉目标需要 q 满足可靠的差分计算。因此， q 要小于磁场强度在基线距离之间的实际变化值。令 r 为MGTS的观测距离，结合 式（8） 式（8） 和 式（10） 式（10） ，传感器测量准确度 q 需满足： <math display="block" id="M26"><mtable><mtr><mtd><mn mathvariant="normal">2</mn><mfenced open="|" close="|" separators="|"><mrow><mi>q</mi></mrow></mfenced><mo>≤</mo><mi>d</mi><mo>⋅</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><msub><mrow><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">B</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">n</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced></mrow><mrow><mtable columnalign="left"><mtr><mtd><mi>x</mi><mo>=</mo><msub><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mtd></mtr><mtr><mtd><mi>y</mi><mo>=</mo><msub><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mtd></mtr><mtr><mtd><mi>z</mi><mo>=</mo><msub><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mtd></mtr></mtable></mrow></msub><mo>=</mo><mi>d</mi><mo>⋅</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>M</mi><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi></mrow></mfenced><mo>.</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29777997&type= 57.31933212 23.36800003 使得在基线距离 d 、该处磁异常强度 B 条件下的差分计算可靠。 差分测量可靠的必要条件如下： <math display="block" id="M27"><mi>r</mi><mo>≤</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>μ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>M</mi><mi>d</mi></mrow><mrow><mn mathvariant="normal">8</mn><mi mathvariant="normal">π</mi><mfenced open="|" close="|" separators="|"><mrow><mi>q</mi></mrow></mfenced></mrow></mfrac><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>θ</mi></mrow></mfenced></mrow></mfenced></mrow><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">4</mn></mrow></mfrac></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778008&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778015&type= 33.35866928 14.13933277 . 式（12） 式（12） 即为差分磁梯度张量测量范围的计算公式。由此可知，某一MGTS系统的理论探测极限不仅与传感器的测量准确度 q 和系统基线距离 d 有关，还与磁源目标的磁矩大小 M 和观测角度 θ 有关。这表明MGTS对同一目标在不同方向上具有不同的探测极限。 由于 式（12） 式（12） 中与磁源相关的磁矩 M 和观测角度 θ 为单一磁源目标参数，当测区内存在多个磁异常目标时， M 应为空间某处的磁偶极子等效源模型磁矩。该等效源是由各个磁异常目标在观测点处的磁场叠加后再反演得到的磁偶极子， θ 应为该等效源磁偶极子的磁矩与观测点间的夹角。 在MATLAB仿真环境中，建立空间笛卡尔坐标系，令 x 轴朝东， y 轴朝北， z 轴朝上。将磁偶极子放置在 r 1 =（10 m，0， - 20 m）处，磁矩 M =5 000 A·m 2 ，磁偏角为20°（东偏），磁倾角为 - 60°（向上）。设置 MGTS的初始参数为基线距离 d =0.5 m，传感器的测量准确度 q =±0.1 nT。 图2 图2 （a）中绘制了MGTS在三维空间中的理论探测极限。 图2 图2 （b）显示了在 z = - 20 m平面内不同 d 和 q 时MGTS的理论探测极限。 图3 图3 （a）显示了空间中不同 θ 和磁矩强度 M 时MGTS的理论探测极限距离。 图3 图3 （b）显示MGTS最大探测极限距离同时随 d 和 q 的变化而变化。 1;2; 1;2; 显然，在平行于磁矩矢量的方向（即 θ =180°或0°），MGTS的探测距离最大。 q 越小， d 越长， M 越大，理论探测极限距离越远。 然而实际情况表明，当探测距离不变而 d 值增长到一定程度时，MGT测量值会失真。由于传感器差分计算得到的张量值与真实值之间存在理论误差，以分量 B xx 为例说明 式（3） 式（3） 计算得到的MGT分量值沿 x 轴梯度测量方向的逼近过程。间隔 d 的两传感器 x 轴测量值 B 1 x 和 B 2 x 分别自 <math id="M28"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778063&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778061&type= 3.72533321 3.97933316 展开为泰勒级数，其中 <math id="M29"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778057&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778056&type= 3.72533321 3.97933316 为观测点处磁场矢量的 x 轴真实分量。令 B xx 测量值为 <math id="M30"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">m</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mrow><mrow><mfenced separators="|"><mrow><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub></mrow></mfenced></mrow><mo>/</mo><mrow><mi>d</mi></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778065&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778074&type= 33.02000046 5.33400011 ，真实值为 <math id="M31"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">t</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mrow><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mo>/</mo><mrow><mo>∂</mo><mi>x</mi></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778068&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778066&type= 20.40466690 4.57200003 ，则 <math id="M32"><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">m</mi><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778072&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778070&type= 5.92666674 3.97933316 可由泰勒级数表示： <math display="block" id="M33"><mtable><mtr><mtd><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">m</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mfrac><mrow><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><msup><mrow><mi>B</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mi>d</mi></mrow></mfrac><mo>=</mo><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mo>∂</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>O</mi></mrow></msub><mo>+</mo><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">3</mn><mo>!</mo></mrow></mfrac><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>O</mi></mrow></msub><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac><mo stretchy="false">(</mo><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">5</mn><mo>!</mo></mrow></mfrac><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>O</mi></mrow></msub><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac><mo stretchy="false">(</mo><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mo>⋯</mo><mo>=</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi><mi>x</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">t</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>+</mo><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">24</mn></mrow></mfrac><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>O</mi></mrow></msub><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mfenced open="" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">1920</mn></mrow></mfrac><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow><mrow><mo>∂</mo><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msup></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>O</mi></mrow></msub><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup><mo>+</mo><mo>⋯</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778076&type= 117.77133942 23.28333473 . 式（13） 式（13） 表明，在差分过程中，由于忽略了泰勒级数的三阶和高阶奇数项，测量值与基线距离 d 存在一定的正相关偏差。当 d 值增大且观测距离很近时，磁场高阶偏微分也变得不可忽略。此时， 式（13） 式（13） 中第二项和后续项与首项相比不能忽略，差分测量结果失真愈发明显。 图4 图4 显示了系统由窗口（10 m，0， - 20 m）滑动测量至（100 m，0， - 20 m），即测点由磁偶极子中心逐渐远离的过程中 B xx 和 B xy 的理论测量误差。因此，随着观测距离变短，由差分计算的MGT测量值可靠性会变低，且 d 越长这种趋势越明显。从这些结果来看，MGTS的理论探测极限与基线距离、测量准确度、磁矩、观测点与磁矩矢量的夹角等有关；基线距离越长，传感器测量准确度越高，MGTS的理论可探测距离越远；基线距离越长，观测距离越近，张量差分测量理论误差越大；探测距离在平行于磁矩矢量的方向上达到最大，而在垂直于磁矩矢量的方向上急剧减小。 在工程实际中，基线距离的设置存在一个较优解，需综合考虑精度需求、距离需求及仪器尺寸等。 4 实　验 1 本文构建了一个实际的平面十字形MGTS，基线距离为0.5 m，其中包含4个Barrington公司生产的Mag-03磁通门传感器，以及一个由非磁性塑性树脂材料制成的十字架，见 图5 图5 。 将 图6 图6 中4种典型的磁体标记为 m 1 ， m 2 ， m 3 和 m 4 ，其中 m 1 和 m 2 为圆柱体， m 3 和 m 4 为长方体。这些磁体的等效偶极矩未知。 实验过程如 图7 图7 所示，以实现典型磁体的MGTS探测极限估计和验证。磁异常测量经验法表明，在超过物体长度约2.5倍的距离处，偶极矩会占主导作用 ［ 21 21 ］ ，磁目标可近似于磁偶极子。将4个磁铁视为远离目标位置的磁偶极子，并试图找到它们相对于几何形状的磁矩矢量。 图7 图7 中，环境测量用以估算传感器和MGTS的实际测量准确度， Q 为MGTS的测量准确度（单位：±nT/m），有 q = Qd /2。采用张量衍生不变关系定位方法 ［ 18 18 ］ 来估计磁铁位置，若位置处于磁铁的物理尺寸内，则表明磁矩估计是有效的。然后，利用放置磁铁前后的MGT测量读数和准确度 Q 判断MGTS是否已到达其探测极限。在测量中，这里采用零相位低通滤波来消除信号干扰。 4.1　传感器和MGTS的测量准确度估计 2 Mag-03传感器的出厂本底噪声幅度在±0.01~0.02 nT内 ［ 10 10 ］ ，然而传感器的真实读数几乎不可能达到此量级。实际上，实测读数不可避免地受到测量设备的电磁干扰，地磁场的本底噪声以及环境中其他未知涡流磁信号的干扰。通过信号处理的方法能够过滤掉一些有规律的高频信号，尽可能保留真实的磁场信息。 实验地点为中国石家庄某空旷野外，划定10 m×10 m且磁场较稳定的测区，标量质子磁强计测得该区域的磁总场强度均值为53 162 nT （±20 nT）。首先进行了环境测量，以确定滤波和降噪的截止频率。采样频率为500 Hz。 图8 图8 （a）、8（b）显示了经过约54 s静态采样后传感器1的 x 轴分量和频谱，并突出了前后各一秒内的读数变化细节。静磁场信号是一种超低频信号，而设备电流采集卡的噪声频段主要集中在50，235 Hz等频段。利用零相位低通滤波器来降低地磁信号噪声，并将截止频率保守地设置为 5 Hz。由于地磁日变的影响以及滤波器中遗漏了部分低频噪声分量，实际磁场强度在不断变化，传感器测量准确度 q 难以直接获得。然而，在MGT数据的差分计算过程中，可以等效地减去该区域地磁场的协同变化。故若已知MGTS的测量准确度 Q ，可以用差分原理估计 q ，即 q = Qd /2。 1;2;3;4; 图8 图8 （c）显示了开始采样后一段时间4个传感器的磁场强度分量读数，各轴读数明显地随地磁场波动而变化。由于采样时间一分钟内温度相对稳定，并且传感器的出厂偏移温度系数（nT/°C）小于±0.1 nT ［ 10 10 ］ ，因此 图8 图8 （c）中的磁场读数漂移几乎全部来自地磁在短时间内的随机日静变化，此时直接估计 q 是不可靠的。然而，在静态采样期间测量的MGT应该是恒定的，因为它代表了磁场强度的空间变化率，不受地磁日变化的影响。MGT读数随采样时间的绝对误差如 图8 图8 （d）所示。Mag-03传感器的出厂测量准确度已达到± 0.02 nT，确保了此时MGT读数的可靠性。 表1 表1 列出了MGTS读数的精确度（均方根误差（Root Mean Square Error， RMSE）、平均误差（Mean Error， ME）和准确度 Q 。 指　标 B xx B xy B xz B yx B yy B yz 估计值 RMSE 3.219 0.916 1.398 1.577 3.126 2.676 1.844 ME 2.647 0.720 1.127 1.228 2.529 2.135 1.484 Q ±8.285 ±3.804 ±4.392 ±5.152 ±8.177 ±6.903 ±8.285 在约54 s的静态采样阶段内，测得MGT分量读数的最大绝对误差收敛到±8.285 nT/m，作为MGTS的估计张量测量准确度 Q 。重复采样后得到的MGT值也在此范围内波动，表明局部磁场相对稳定且 Q 值估计有效。由于 q = Qd /2，给定 d =0.5 m，故 q 真值应收敛于±2.121 nT附近，作为现有仪器和降噪条件下估计的传感器测量准确度。 4.2　MGTS的探测极限估计 2 一旦获得MGT分量和位置矢量，可直接利用 式（5） 式（5） 磁偶极子正演方程来反演磁矩矢量 m ，即： <math display="block" id="M34"><mi mathvariant="bold-italic">m</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>C</mi></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mo>+</mo></mrow></msup><mi mathvariant="bold-italic">I</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>C</mi></mrow></mfrac><msup><mrow><mfenced separators="|"><mrow><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">H</mi></mrow></mfenced></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><msup><mrow><mi mathvariant="bold-italic">H</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778148&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778131&type= 50.96933365 7.87400007 ， 其中： H 仅与位置矢量有关， I 是具有5个独立MGT分量的列向量。请注意， 式（13） 式（13） 中的 H + 是Moore-Pennrose逆。然而，现有情况不能直接判断估计的 m 是否有效。环境中可能存在其他磁异常，使得由反演公式得到的叠加磁场等效磁矩代替了目标的真实磁矩。此外，前文结论表明，距离过近同样会使张量差分测量值变得不可靠。 为了确保磁矩估计是有效的，利用张量衍生不变关系法定位磁体 ［ 18 18 ］ 。该方法从以下方程中提供4个可能的坐标解，其中一个是正确的。 <math display="block" id="M35"><mtable><mtr><mtd><mi mathvariant="bold-italic">r</mi><mo>=</mo><mo>±</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">M</mi></mrow></msub><mroot><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>θ</mi><mo>+</mo><mn mathvariant="normal">2</mn></mrow><mrow/></mroot></mrow><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msub><mroot><mrow><mn mathvariant="normal">3</mn><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>θ</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow><mrow/></mroot></mrow></mfrac><mroot><mrow><mfrac><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>-</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot><mfrac><mrow><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced></mrow></mfrac><mo>±</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">M</mi></mrow></msub><mroot><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>θ</mi><mo>+</mo><mn mathvariant="normal">2</mn></mrow><mrow/></mroot></mrow><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msub><mroot><mrow><mn mathvariant="normal">3</mn><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>θ</mi><mo>+</mo><mn mathvariant="normal">1</mn></mrow><mrow/></mroot></mrow></mfrac><mroot><mrow><mfrac><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub></mrow><mrow><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfrac></mrow><mrow/></mroot><mfrac><mrow><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced></mrow></mfrac><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778136&type= 64.09266663 24.80733490 式中 v 1 和 v 2 分别对应于特征值 λ 1 和 λ 2 的特征向量， C T 是张量缩并， <math id="M36"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">M</mi></mrow></msub><mo>=</mo><mroot><mrow><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=29778138&type= 33.52799988 5.24933338 是仅由磁偶极子产生的磁总场强度TMI。因此，在计算 I TM 时应提前测量地磁场进行补偿。 对坐标已知的磁体进行定位，一旦估计坐标与目标真实位置间的测量偏差控制在磁体尺寸范围内，则反演得到的磁偶极矩是有效的。MGTS置于地面高50 cm处，以观测点为原点，将4块磁铁分别放置于位置1（50 cm，50 cm， - 50 cm）、位置2（100 cm，100 cm， - 50 cm），如 图9 图9 所示。 对4块磁铁共进行了4×2次定位实验。采样频率为500 Hz，单次采样时间为10 s左右。定位结果为每单次采样时间内计算的平均值。设置相同的滤波条件，截止频率均为5 Hz。定位结果列于 表2 表2 。显然，位置1处的定位结果均有偏差，而位置2处的定位结果均控制在4个磁体的几何空间范围内。位置1处由于探测距离过近而差分测量失真，而位置2的结果都是有效的。4个磁铁估计的磁矩矢量如 图10 图10 所示。 磁铁 尺寸 /cm 真实位置 /cm 估计位置 /cm 估计 m / （A·m 2 ） 估计 M / （A·m 2 ） 是否有效 2 m 1 2 ϕ2.0× 2.0 （50，50，-50） （66.5，67.0，-52.8） （1.68， 2.41， -6.61） 7.24 否 （100，100，-50） （100.0，100.3，-49.9） （0.04，0.08，-4.13） 4.13 是 2 m 2 2 ϕ6.0× 1.3 （50，50，-50） （66.9，65.1，53.0） （4.98， 2.76， -11.82） 13.12 否 （100，100，-50） （100.0，99.6，-50.7） （0.45， 0.36， -6.67） 6.69 是 2 m 3 2 6.0× 4.0×1.0 （50，50，-50） （69.8，62.7，-55.7） （7.74， 2.44， -15.75） 17.72 否 （100，100，-50） （102.7，103.6，-50.4） （-0.63， -0.81， -11.46） 11.51 是 2 m 4 2 10.0× 3.0×1.0 （50，50，-50） （64.6，67.4，-56.4） （-3.44， 10.05， -32.57） 34.26 否 （100，100，-50） （101.3，104.5，-50.7） （-2.84， -0.27， -22.77） 22.9467 是 一旦确定了磁矩的方向和大小、传感器读数准确度以及MGTS的基线距离，就可以得到该MGTS对4个磁体的空间理论探测极限范围，如 图11 图11 （a）和11（b）所示。此外， 图11 图11 （c）显示了4个磁铁的探测极限随角度 θ 的变化。MGTS相关参数和最大探测极限距离（ r max ）的估计结果列于 表3 表3 。 1;2;3; d /m q /nT θ /（°） 磁铁 m /（A·m 2 ） M /（A·m 2 ） 估计 r max /m 实测 r max /m 4 0.5 4 ±2.1 4 0 or 180 m 1 （0，0.1，-4.1） 4.1 6.2 5.8 m 2 （0.4，0.4，-6.7） 6.7 7.0 6.6 m 3 （-0.6，-0.8，-11.5） 11.5 8.0 7.7 m 4 （-2.8，-0.3，-22.8） 23.0 9.5 9.2 为了验证估计结果，将4块磁铁沿MGTS的 x 轴滑动并连续采样。采样率为500 Hz，滑动窗口为（1 m，0，0）到（10 m，0，0）。规定每0.01 m采集不少于20个点时，实测MGT的准确度可信，故滑动速度应小于0.25 m/s。设置滑动速度为0.2 m/s，速度偏差为±0.05 m/s，即可满足要求。滑动过程中始终保持角度 θ 为0°，即磁矩 m 的方向始终指向观察点。 通过观察是否可以区分不同MGT分量读数来判断MGTS是否已达到对4种磁铁的真实探测极限边界，并估计出实测 r max 。为了消除环境中其他未知磁异常对结果的干扰，利用背景磁场的空侧数据对滑动采样阶段测量的磁场信号进行补偿。一旦磁体滑动时测得的MGT分量读数之差不再超过MGTS的测量准确度 Q （±8.285 nT/m），判定磁体已脱离了系统的有效探测范围。测量结果如 图12 图12 所示，图中标出了每个磁体的实测探测极限 r max ，并记录在 表3 表3 中，其中测量值和估计值最大偏差为0.4 m。 上述结果表明，MGTS实测探测极限与估计值吻合较好，探测极限估计准确度为±0.4 m，从而验证了所提差分磁梯度张量测量极限估计方法的有效性。 5 结　论 1 本文提出了差分MGT测量极限估计方法。根据差分MGT测量范围公式，MGTS的理论探测极限与基线距离、传感器测量准确度、目标磁偶极矩、观测点与磁矩矢量间的夹角有关。基线距离越长，传感器测量准确度越高，MGTS的有效探测距离越远。探测距离在平行于磁矩矢量的方向上达到最大，而在垂直于磁矩矢量的方向上急剧减小。在实际测量中，成功估计并验证了搭建的平面十字形MGTS的理论探测极限，该系统针对4块典型磁铁的探测极限估计准确度为±0.4 m。然而，本文尚未深入考虑地磁日变化、涡流磁场干扰和磁传感器非线性误差等因素对实测中系统探测极限的影响。 参考文献 周子童 ， 闫韶华 ， 赵巍胜 ， 等 . 隧穿磁阻传感器研究进展 ［J］. 物理学报 ， 2022 ， 71 （ 5 ）： 058504 . doi: 10.7498/aps.71.20211883 http://dx.doi.org/10.7498/aps.71.20211883 隧穿磁阻传感器研究进展 ZHOU Z T ， YAN SH H ， ZHAO W SH ， et al . Research progress of tunneling magnetoresistance sensor ［J］. Acta Physica Sinica ， 2022 ， 71 （ 5 ）： 058504 . （in Chinese） . doi: 10.7498/aps.71.20211883 http://dx.doi.org/10.7498/aps.71.20211883 Research progress of tunneling magnetoresistance sensor 李金朋 ， 任国全 ， 张英堂 ， 等 . 改进二维变分模态分解的磁源分离 ［J］. 光学 精密工程 ， 2020 ， 28 （ 5 ）： 1200 - 1211 . doi: 10.3788/OPE.20202805.1200 http://dx.doi.org/10.3788/OPE.20202805.1200 改进二维变分模态分解的磁源分离 LI J P ， REN G Q ， ZHANG Y T ， et al . Source separation based on improved two-dimensional variational mode decomposition ［J］. Optics and Precision Engineering ， 2020 ， 28 （ 5 ）： 1200 - 1211 . （in Chinese） . doi: 10.3788/OPE.20202805.1200 http://dx.doi.org/10.3788/OPE.20202805.1200 Source separation based on improved two-dimensional variational mode decomposition 刘世斌 . 磁梯度张量测量系统野外现场校正与标定方法研究 ［D］. 长春 ： 吉林大学 ， 2021 . doi: 10.23977/agrfem.2021.040105 http://dx.doi.org/10.23977/agrfem.2021.040105 磁梯度张量测量系统野外现场校正与标定方法研究 LIU SH B . Research on Correction and Calibration Methods for Magnetic Gradient Tensor Measurement Systems in the Field Surveys ［D］. Changchun ： Jilin University ， 2021 . （in Chinese） . doi: 10.23977/agrfem.2021.040105 http://dx.doi.org/10.23977/agrfem.2021.040105 Research on Correction and Calibration Methods for Magnetic Gradient Tensor Measurement Systems in the Field Surveys ZHENG C ， ZHU K ， CARDOSO DE FREITAS S ， et al . Magnetoresistive sensor development roadmap （non-recording applications） ［J］. IEEE Transactions on Magnetics ， 2019 ， 55 （ 4 ）： 1 - 30 . doi: 10.1109/tmag.2019.2896036 http://dx.doi.org/10.1109/tmag.2019.2896036 Magnetoresistive sensor development roadmap （non-recording applications） 张云洲 ， 吴峻 . Prewitt霍尔磁梯度张量系统结构设计 ［J］. 仪器仪表学报 ， 2021 ， 42 （ 9 ）： 115 - 123 . Prewitt霍尔磁梯度张量系统结构设计 ZHANG Y ZH ， WU J . Structure design of Prewitt hall magnetic gradient tensor system ［J］. Chinese Journal of Scientific Instrument ， 2021 ， 42 （ 9 ）： 115 - 123 . （in Chinese） Structure design of Prewitt hall magnetic gradient tensor system WANG C ， QU X D ， ZHANG X J ， et al . A fast calibration method for magnetometer array and the application of ferromagnetic target localization ［J］. IEEE Transactions on Instrumentation and Measurement ， 2017 ， 66 （ 7 ）： 1743 - 1750 . doi: 10.1109/tim.2017.2668558 http://dx.doi.org/10.1109/tim.2017.2668558 A fast calibration method for magnetometer array and the application of ferromagnetic target localization STOLZ R ， ZAKOSARENKO V ， SCHULZ M ， et al . Magnetic full-tensor SQUID gradiometer system for geophysical applications ［J］. The Leading Edge ， 2006 ， 25 （ 2 ）： 178 - 180 . doi: 10.1190/1.2172308 http://dx.doi.org/10.1190/1.2172308 Magnetic full-tensor SQUID gradiometer system for geophysical applications LI Q Z ， LI Z N ， ZHANG Y T ， et al . Integrated compensation and rotation alignment for three-axis magnetic sensors array ［J］. IEEE Transactions on Magnetics ， 2018 ， 54 （ 10 ）： 1 - 11 . doi: 10.1109/tmag.2018.2860950 http://dx.doi.org/10.1109/tmag.2018.2860950 Integrated compensation and rotation alignment for three-axis magnetic sensors array WANG B ， REN G Q ， LI Z N ， et al . The stability optimization algorithm of second-order magnetic gradient tensor ［J］. AIP Advances ， 2021 ， 11 （ 7 ）： 075322 . doi: 10.1063/5.0056361 http://dx.doi.org/10.1063/5.0056361 The stability optimization algorithm of second-order magnetic gradient tensor 周鹤 . 基于μCOS内核的磁通门磁力仪研制 ［D］. 长春 ： 吉林大学 ， 2016 . 基于μCOS内核的磁通门磁力仪研制 ZHOU H . The Development of Fluxgate Magnetometer Based on an μ COS System ［D］. Changchun ： Jilin University ， 2016 . （in Chinese） The Development of Fluxgate Magnetometer Based on an μCOS System DUMAS R K ， HOGAN T . Recent Advances in SQUID Magnetometry ［J］. Magnetic Measurement Techniques for Materials Characterization ， 2021 ： 39 - 62 . doi: 10.1007/978-3-030-70443-8_3 http://dx.doi.org/10.1007/978-3-030-70443-8_3 Recent Advances in SQUID Magnetometry 李青竹 ， 李志宁 ， 张英堂 ， 等 . 磁梯度张量系统传感器阵列的快速旋转校准 ［J］. 光学 精密工程 ， 2018 ， 26 （ 7 ）： 1813 - 1826 . doi: 10.3788/OPE.20182607.1813 http://dx.doi.org/10.3788/OPE.20182607.1813 磁梯度张量系统传感器阵列的快速旋转校准 LI Q ZH ， LI ZH N ， ZHANG Y T ， et al . Fast rotation calibration of sensor array in magnetic gradient tensor system ［J］. Optics and Precision Engineering ， 2018 ， 26 （ 7 ）： 1813 - 1826 . （in Chinese） . doi: 10.3788/OPE.20182607.1813 http://dx.doi.org/10.3788/OPE.20182607.1813 Fast rotation calibration of sensor array in magnetic gradient tensor system LI Q Z ， SHI Z Y ， LI Z N ， et al . Magnetic object positioning based on second-order magnetic gradient tensor system ［J］. IEEE Sensors Journal ， 2021 ， 21 （ 16 ）： 18237 - 18248 . doi: 10.1109/jsen.2021.3085573 http://dx.doi.org/10.1109/jsen.2021.3085573 Magnetic object positioning based on second-order magnetic gradient tensor system PEI Y H ， YEO H G . UXO survey using vector magnetic gradiometer on autonomous underwater vehicle ［C］. OCEANS 2009 . October 26-29， 2009 ， Biloxi， MS， USA . IEEE ， 2009 ： 1 - 8 . doi: 10.23919/oceans.2009.5422446 http://dx.doi.org/10.23919/oceans.2009.5422446 UXO survey using vector magnetic gradiometer on autonomous underwater vehicle 李青竹 ， 李志宁 ， 张英堂 ， 等 . 平面十字磁梯度张量系统的两步线性校正 ［J］. 仪器仪表学报 ， 2017 ， 38 （ 9 ）： 2232 - 2242 . doi: 10.3969/j.issn.0254-3087.2017.09.018 http://dx.doi.org/10.3969/j.issn.0254-3087.2017.09.018 平面十字磁梯度张量系统的两步线性校正 LI Q ZH ， LI ZH N ， ZHANG Y T ， et al . Two-step linear calibration of planar cross magnetic gradient tensor system ［J］. Chinese Journal of Scientific Instrument ， 2017 ， 38 （ 9 ）： 2232 - 2242 . （in Chinese） . doi: 10.3969/j.issn.0254-3087.2017.09.018 http://dx.doi.org/10.3969/j.issn.0254-3087.2017.09.018 Two-step linear calibration of planar cross magnetic gradient tensor system COCCHI L ， CARMISCIANO C ， PALANGIO P ， et al . S3MAG-Low magnetic noise AUV for multipurpose investigations ［C］. OCEANS 2015 - Genova . 1821，2015 ， Genova ， Italy . IEEE ， 2015 ： 1 - 4 . doi: 10.1109/oceans-genova.2015.7271480 http://dx.doi.org/10.1109/oceans-genova.2015.7271480 S3MAG-Low magnetic noise AUV for multipurpose investigations LI Q Z ， SHI Z Y ， LI Z N ， et al . Preferred configuration and detection limits estimation of magnetic gradient tensor system ［J］. IEEE Transactions on Instrumentation and Measurement ， 2021 ， 70 ： 1 - 14 . doi: 10.1109/tim.2021.3109741 http://dx.doi.org/10.1109/tim.2021.3109741 Preferred configuration and detection limits estimation of magnetic gradient tensor system 李青竹 ， 李志宁 ， 张英堂 ， 等 . 张量衍生不变关系下的磁源单点定位 ［J］. 光学 精密工程 ， 2019 ， 27 （ 8 ）： 1880 - 1893 . doi: 10.13810/j.cnki.issn.1000-7210.2019.04.024 http://dx.doi.org/10.13810/j.cnki.issn.1000-7210.2019.04.024 张量衍生不变关系下的磁源单点定位 LI Q ZH ， LI ZH N ， ZHANG Y T ， et al . Magnetic source single-point positioning by tensor derivative invariant relations ［J］. Optics and Precision Engineering ， 2019 ， 27 （ 8 ）： 1880 - 1893 . （in Chinese） . doi: 10.13810/j.cnki.issn.1000-7210.2019.04.024 http://dx.doi.org/10.13810/j.cnki.issn.1000-7210.2019.04.024 Magnetic source single-point positioning by tensor derivative invariant relations CLARK D A . New methods for interpretation of magnetic vector and gradient tensor data I： eigenvector analysis and the normalised source strength ［J］. Exploration Geophysics ， 2012 ， 43 （ 4 ）： 267 - 282 . doi: 10.1071/eg12020 http://dx.doi.org/10.1071/eg12020 New methods for interpretation of magnetic vector and gradient tensor data I： eigenvector analysis and the normalised source strength MOSCHITTA A ， DE ANGELIS A ， SANTONI F ， et al . Estimation of the magnetic dipole moment of a coil using AC voltage measurements ［J］. IEEE Transactions on Instrumentation and Measurement ， 2018 ， 67 （ 10 ）： 2495 - 2503 . doi: 10.1109/tim.2018.2828700 http://dx.doi.org/10.1109/tim.2018.2828700 Estimation of the magnetic dipole moment of a coil using AC voltage measurements 张朝阳 ， 肖昌汉 ， 高俊吉 ， 等 . 磁性物体磁偶极子模型适用性的试验研究 ［J］. 应用基础与工程科学学报 ， 2010 ， 18 （ 5 ）： 862 - 868 . doi: 10.3969/j.issn.1005-0930.2010.05.016 http://dx.doi.org/10.3969/j.issn.1005-0930.2010.05.016 磁性物体磁偶极子模型适用性的试验研究 ZHANG ZH Y ， XIAO CH H ， GAO J J ， et al . Experiment research of magnetic dipole model applicability for a magnetic object ［J］. Journal of Basic Science and Engineering ， 2010 ， 18 （ 5 ）： 862 - 868 . （in Chinese） . doi: 10.3969/j.issn.1005-0930.2010.05.016 http://dx.doi.org/10.3969/j.issn.1005-0930.2010.05.016 Experiment research of magnetic dipole model applicability for a magnetic object 10.37188/OPE.20223011.1325 TH763 TP212.9 A 1004-924X（2022）11-1325-12 李青竹,石志勇,李志宁等.差分磁梯度张量测量极限估计[J].光学精密工程,2022,30(11):1325-1336. 李青竹,石志勇,李志宁等.差分磁梯度张量测量极限估计[J].光学精密工程,2022,30(11):1325-1336. LI Qingzhu,SHI Zhiyong,LI Zhining,et al.Estimation method for differential magnetic gradient tensor measurement limits[J].Optics and Precision Engineering,2022,30(11):1325-1336. LI Qingzhu,SHI Zhiyong,LI Zhining,et al.Estimation method for differential magnetic gradient tensor measurement limits[J].Optics and Precision Engineering,2022,30(11):1325-1336. 水下测控技术重点实验室基金资助项目（No.6142407190307）；军内科研资助项目 2022-01-11 2022-02-17 2022-06-10 2022-06-21 光学精密工程 11 30