自适应滤波的磁共振信号抗饱和消噪装置及消噪方法
【专利摘要】本发明涉及一种自适应滤波的磁共振信号抗饱和消噪装置及消噪方法，采用三个多匝小线圈作为参考通道，前端处理模块经高精度AD采集模块连接至FPGA模块，再经数据传输模块连接至PC模块，避免了核磁信号的感应。通过FPGA中抗饱和模块使信号被尽可能放大且不饱和；再用自适应消噪算法对信号进行实时参考消噪，利用PC模块实时修改自适应滤波器参数，能适应各种复杂的探测环境。本发明的消噪装置和消噪方法，经在强电磁干扰环境中实验，不仅保证核磁信号不失真，提高信噪比，抑制工频谐波干扰，还起到抗饱和的作用，并能在低信噪比且易饱和的探测现场实时检测到核磁信号，为后续核磁信号的数据反演奠定可靠基础。
【专利说明】
自适应滤波的磁共振信号抗饱和消噪装置及消噪方法
技术领域
[0001] 本发明设及一种数据采集装置及数字信号处理方法，尤其适用于核磁共振 (Nuclear Magnetic Resonance,NMR)探水信号的采集并进行抗饱和处理和自适应实时参 考消噪。
【背景技术】
[0002] 利用核磁共振原理找水时，由于探测线圈接收到的核磁信号是非常微弱（纳伏级） 的，很容易受尖峰噪声，工频谐波和随机噪声的影响。尤其是在矿井隧道、城市中屯、运种强 电磁干扰环境中，核磁信号基本处于被工频谐波淹没且饱和的状态，常规的滤波方法不仅 无法有效提取核磁信号，还可能使信号饱和失真。
[0003] CN201010537465.2公开了一种带有参考线圈的核磁共振地下水探测系统及探测 方法，通过多路模数转换单元同步采集主通道和参考通道的全波形数据，采用的变步长自 适应算法来滤除工频谐波噪声，提高仪器抗干扰能力。但是，在采集完全波数据后，再进行 自适应消噪，没有现场实时处理，信号还是会被淹没在噪声中，不能保证数据的可靠性;另 夕h其变步长算法的步长取决于被测区信噪比。在进行自适应滤波前需要计算被测区信噪 比，过程相当繁琐;除此之外，还要使用几十阶变步长自适应滤波器，才能取得一定效果。
[0004] CN201410142630.2公开了一种多个近端参考线圈的核磁共振信号实时噪声抵消 装置，由主接收线圈与一个W上参考收线圈分别经地面核磁共振仪器的多通道信号调理电 路后与数模转换器连接，转换成数字量进入数字信号处理化SP)忍片，数模转换器与主自适 应滤波器相连，参考信号通道数模转换器分别于自适应滤波器连接构成。其利用数字信号 处理忍片来做自适应算法，速度较慢，难W实时实现;采用的最小均方(LMS)算法、归一化最 小均方(NLMS)算法，收敛速度慢，稳态误差较大;而最小二乘算法(RLS)计算量过大，也很难 实时实现。另外，该专利也不能解决强电磁干扰环境下信号饱和的问题。
【发明内容】
：
[0005] 本发明的目的就是针对上述技术存在的问题，提供一种适合于在矿井隧道、城市 中屯、运种强电磁干扰环境中探测水源的抗饱和实时消噪装置
[0006] 本发明的另一目的是提供一种自适应滤波的磁共振信号抗饱和消噪装置的消噪 方法。
[0007] 本发明的目的是通过W下技术方案实现的：
[000引自适应滤波的磁共振信号抗饱和消噪装置，是由参考线圈曰1经前端处理模块曰4与 AD采集模块8连接，探测线圈0经前端处理模块b5与AD采集模块8连接，参考线圈b2经前端处 理模块c6与AD采集模块8连接，参考线圈c3经前端处理模块d7与AD采集模块8连接，前端处 理模块a4经前端处理模块b5、前端处理模块c6、前端处理模块d7和FPGA模块9分别与AD采集 模块8和前端处理模块a4连接，AD采集模块8经FPGA模块9、数据传输模块10和PC模块11与 FPGA模块9连接构成。
[0009] 前端处理模块a4由继电器控制电路al2经低噪声前放电路al3和带通滤波器al4与 程控放大器al5连接构成；前端处理模块b5由继电器控制电路bl6经低噪声前放电路bl7和 带通滤波器bl8与程控放大器bl9连接构成;前端处理模块c6由继电器控制电路c20经低噪 声前放电路c21和带通滤波器c22与程控放大器c23连接构成；前端处理模块d7由继电器控 制电路d24经低噪声前放电路d25和带通滤波器d26与程控放大器d27连接构成。
[0010] FPGA模块9由AD控制模块28，抗饱和模块29和自适应滤波模块30组成。其中AD控制 模块28由PC11控制产生同步信号;抗饱和模块29自动调节程控放大器al 5、程控放大器b 19、 程控放大器c23和程控放大器d27的放大倍数。
[0011] 自适应滤波的磁共振信号抗饱和消噪装置的消噪方法，包括W下步骤：
[0012] a.铺设一个方形探测线圈0和Ξ个多应参考线圈al、参考线圈b2、参考线圈c3;参 考线圈尽可能靠近干扰源铺设，参考线圈尽可能靠近干扰源铺设，参考线圈中屯、到探测线 圈中屯、的距离大于正方形探测线圈的边长。
[001引b.系统上电，通过P对莫块11初始化各模块；
[0014] C.在采集开始时，PC模块11发出采集控制指令，经数据传输模块10送到FPGA模块9 中，由AD控制模块28产生同步信号，驱动继电器控制电路al2、继电器控制电路bl6、继电器 控制电路c20和继电器控制电路d24与高精度AD采集模块8同步工作；
[0015] d.继电器打开后，信号经继电器控制电路a 12、继电器控制电路b 16、继电器控制电 路c20和继电器控制电路d24到低噪声前放电路曰13、低噪声前放电路bl7、低噪声前放电路 c21和低噪声前放电路d25放大，之后经宽频带带通滤波器al4、宽频带带通滤波器bl8、宽频 带带通滤波器c22和宽频带带通滤波器d26滤波放大，再通过程控放大器al5、程控放大器 bl9、程控放大器c23和程控放大器d27被再次放大，模拟信号dO(t)、Xl(t)、X2(t)和X3(t)通 过高精度AD采集模块8转换成数字信号d0(n)、X1 (η)、X2(n)和X3(n);
[0016] e.FPGA模块9产生采集时钟频率并同步地将四路数字信号读入FPGA忍片中，经抗 饱和模块29判断信号是否饱和，饱和时，FPGA模块会自适应的调节程控放大器al5、程控放 大器bl9、程控放大器c23和程控放大器d27的倍数直到信号处于临界饱和状态；
[0017] f.抗饱和模块29的四个不饱和信号同时进入自适应滤波模块30,进行自适应实时 消噪处理；
[0018] g.自适应滤波模块30的输出及自适应滤波器的权系数，经数据传输模块10打包后 送到P对莫块11中，并显示波形图；
[0019] h.观察自适应滤波器的权系数，判断权系数是否快速稳定收敛，对不收敛、慢收敛 或稳态误差明显的，通过PC模块11重新设置自适应滤波器参数，PC模块11直接将运些参数 送到FPGA模块9中的自适应滤波模块30，进行滤波器参数的修改；
[0020] i.修改好滤波器参数后，从步骤b到步骤h重复进行，直到权系数快速稳定的收敛 为止，固定滤波器参数，进行不同深度水体的探测，存储数据，并作核磁信号的数据反演。
[0021] 有益效果：采用Ξ个多应小线圈作为参考通道，Ξ个多应小线圈与探测线圈的距 离大于探测线圈的边长，可避免了核磁信号的感应。通过FPGA中抗饱和模块使信号被尽可 能放大且不饱和;再利用自适应消噪算法对信号进行实时参考消噪。利用PC模块实时修改 自适应滤波器参数，能适应各种复杂的探测环境。在强电磁干扰环境中，本发明的消噪装置 和消噪方法，不仅保证核磁信号不失真，还起到抗饱和、抑制工频谐波干扰的作用，并能在 低信噪比且易饱和的探测现场实时检测到核磁信号，为后续核磁信号的数据反演奠定可靠 基础。
【附图说明】
[0022] 图1为自适应滤波的磁共振信号抗饱和消噪装置结构框图
[0023] 图2为附图1中前端处理模块结构框图
[0024] 图3为附图1中FPGA模块9的结构框图
[0025] 图4为附图3中自适应滤波模块30算法原理图
[0026] 图5为Ξ个参考信号的5阶变步长梯度自适应格型联合算法原理图
[0027] 图6为参考线圈与探测线圈野外铺设图
[002引 0探测线圈，1参考线圈a，2参考线圈b，3参考线圈C，4前端处理模块a，5前端处理模 块b，6前端处理模块C，7前端处理模块d，8高精度AD采集模块，9FPGA模块，10数据传输模块， 11PC模块，12继电器控制电路a，13低噪声前放电路a，14带通滤波器a，15程控放大器a，16继 电器控制电路b，17低噪声前放电路b，18带通滤波器b，19程控放大器b，20继电器控制电路 c，21低噪声前放电路c，22带通滤波器c，23程控放大器c，24继电器控制电路d，25低噪声前 放电路d，26带通滤波器d，27程控放大器d，28AD控制模块，29抗饱和模块，30自适应滤波模 块。
【具体实施方式】
[0029] 下面结合附图与实施例对本发明作进一步的详细说明：
[0030] 自适应滤波的磁共振信号抗饱和消噪装置，是由参考线圈al经前端处理模块a4与 AD采集模块8连接，探测线圈0经前端处理模块b5与AD采集模块8连接，参考线圈b2经前端处 理模块c6与AD采集模块8连接，参考线圈c3经前端处理模块d7与AD采集模块8连接，前端处 理模块a4经前端处理模块b5、前端处理模块c6、前端处理模块d7和FPGA模块9分别与AD采集 模块8和前端处理模块a4连接，AD采集模块8经FPGA模块9、数据传输模块10和PC模块11与 FPGA模块9连接构成。
[0031] 前端处理模块a4由继电器控制电路al2经低噪声前放电路al3和带通滤波器al4与 程控放大器al5连接构成；前端处理模块b5由继电器控制电路bl6经低噪声前放电路bl7和 带通滤波器bl8与程控放大器bl9连接构成;前端处理模块c6由继电器控制电路c20经低噪 声前放电路c21和带通滤波器c22与程控放大器c23连接构成；前端处理模块d7由继电器控 制电路d24经低噪声前放电路d25和带通滤波器d26与程控放大器d27连接构成。
[0032] FPGA模块9由AD控制模块28，抗饱和模块29和自适应滤波模块30组成。其中AD控制 模块28由PC11控制产生同步信号;抗饱和模块29自动调节程控放大器al 5、程控放大器b 19、 程控放大器c23和程控放大器d27的放大倍数。
[0033] 自适应滤波的磁共振信号抗饱和消噪装置的消噪方法，包括W下步骤：
[0034] a.铺设一个方形探测线圈0和Ξ个多应参考线圈al、参考线圈b2、参考线圈c3;参 考线圈尽可能靠近干扰源铺设，参考线圈尽可能靠近干扰源铺设，参考线圈中屯、到探测线 圈中屯、的距离大于正方形探测线圈的边长。
[0035] b.系统上电，通过P对莫块11初始化各模块；
[0036] c.在采集开始时，PC模块11发出采集控制指令，经数据传输模块10送到FPGA模块9 中，由AD控制模块28产生同步信号，驱动继电器控制电路al2、继电器控制电路bl6、继电器 控制电路c20和继电器控制电路d24与高精度AD采集模块8同步工作；
[0037] d.继电器打开后，信号经继电器控制电路a 12、继电器控制电路b 16、继电器控制电 路c20和继电器控制电路d24到低噪声前放电路曰13、低噪声前放电路bl7、低噪声前放电路 c21和低噪声前放电路d25放大，之后经宽频带带通滤波器al4、宽频带带通滤波器bl8、宽频 带带通滤波器c22和宽频带带通滤波器d26滤波放大，再通过程控放大器al5、程控放大器 bl9、程控放大器c23和程控放大器d27被再次放大，模拟信号dO(t)、Xl(t)、X2(t)和X3(t)通 过高精度AD采集模块8转换成数字信号d0(n)、X1 (η)、X2(n)和X3(n);
[003引e.FPGA模块9产生采集时钟频率并同步地将四路数字信号读入FPGA忍片中，经抗 饱和模块29判断信号是否饱和，饱和时，FPGA模块会自适应的调节程控放大器al5、程控放 大器bl9、程控放大器c23和程控放大器d27的倍数直到信号处于临界饱和状态；
[0039] f.抗饱和模块29的四个不饱和信号同时进入自适应滤波模块30,进行自适应实时 消噪处理；
[0040] g.自适应滤波模块30的输出及自适应滤波器的权系数，经数据传输模块10打包后 送到P对莫块11中，并显示波形图；
[0041] h.观察自适应滤波器的权系数，判断权系数是否快速稳定收敛，对不收敛、慢收敛 或稳态误差明显的，通过PC模块11重新设置自适应滤波器参数，PC模块11直接将运些参数 送到FPGA模块9中的自适应滤波模块30，进行滤波器参数的修改；
[0042] i.修改好滤波器参数后，从步骤b到步骤h重复进行多次，直到权系数快速稳定的 收敛为止，固定滤波器参数，进行不同深度水体的探测，存储数据，并作核磁信号的数据反 演。
[0043] 附图1是自适应滤波的磁共振信号抗饱和消噪装置结构框图，是由探测线圈0与参 考线圈1、2、3分别经前端处理模块5、4、6、7后与AD采集模块8相连，转换成数字量后进入 FPGA模块9中，经数据传输模块10连接至PC11。
[0044] 附图2是前端处理模块内部结构框图，前端处理模块4由继电器控制电路12,低噪 声前放电路13,带通滤波器14,程控放大器15连接构成;前端处理模块5由继电器控制电路 16,低噪声前放电路17,带通滤波器18,程控放大器19连接构成;前端处理模块6由继电器控 制电路20,低噪声前放电路21，带通滤波器22,程控放大器23连接构成;前端处理模块7由继 电器控制电路24,低噪声前放电路25,带通滤波器26,程控放大器27连接构成。
[0045] 附图3是FPGA模块内部结构框图，FPGA模块9主要由AD控制模块28,抗饱和模块29 和自适应滤波模块30组成。其中AD控制模块28由PC11控制产生同步信号;抗饱和模块29自 动调节程控放大器15、19、23、27的放大倍数。
[0046] 附图4是自适应滤波模块算法的原理图，由多级格型预测器和VSLMS滤波器构成。
[0047] PC模块11控制自适应滤波模块30的滤波器参数。
[0048] 附图5是变步长梯度自适应格型联合算法原理图，是由格型预测器和变步长LMS期 望响应估计器(VSLMS滤波器)组成，其中格型预测器迭代公式为：
[0049] f〇(n) =bo(n) =X(n)
[0050] fm(n) = fm-i(n)+kmbm-i(n-l)
[0060] 其中VSLMS滤波器的时变步长参数是由记忆因子U(n-l)和指数因子s-Sw+Λ。》组 成。算法的稳态误差主要由常参量α决定，收敛速度主要由常参量β决定。通过修改运两个参 数，可W使算法拥有优越的跟踪性能与较好的稳态误差。同一参考信号的所有反射系数km (n)(m=l，2……Μ)的步长U都设置一样；同一参考信号的所有记忆因子的常参量α和指数因 子的常参量β也都分别设置一样，W便现场实时修改，尽快达到最好的滤波性能。
[0061] 附图6是线圈铺设的平面示意图，一个大的方形探测线圈0,3个多应小型参考线圈 1，参考线圈2,参考线圈3;探测线圈0与参考线圈1，参考线圈2,参考线圈3的中屯、距离要大 于方形探测线圈0的边长，参考线圈尽可能靠近干扰源。
[0062] 实施例1
[0063] 在长春市烧锅镇做抗饱和自适应实时参考消噪实验。具体实施步骤如下：
[0064] a.铺设100m xlOOm的方形线圈作为发射线圈，W发射线圈的中屯、为基准，铺设25m x25m的探测线圈0,在距离基准30m左右的3个参考点铺设20应50cm x50cm的参考线圈1、2、 3。
[0065] b.对整个系统上电，通过PC初始化各模块(程控放大器15、19、23、27的倍数设为最 大的16倍）。
[0066] C.在采集开始时，PC11发出采集的相关控制指令，经数据传输模块10送到FPGA模 块9中，由AD控制模块28产生同步信号，驱动继电器控制电路12、16、20、24与高精度40采集 模块8同步开始工作。
[0067] d.在继电器打开后，信号经过继电器控制电路12、16、20、24，到低噪声前放电路 13、17、21、25被放大，之后被宽频带带通滤波器14、18、22、26滤波放大，再通过程控放大器 15、19、23、27被再次放大。模拟信号(10(〇心1(〇心2(〇心3(〇通过高精度40采集模块8转 换成数字信号 dO(n)、Xl(n)、X2(n)、X3(n)。
[006引e .FPGA模块9控制采集时钟频率并同步地将数字信号读入FPGA忍片中。经抗饱和 模块29判断信号不饱和。
[0069] f.4个通道的不饱和信号同时进入自适应滤波模块30,采用图4中的变步长梯度自 适应格型联合算法对其进行自适应实时消噪处理。
[0070] g.自适应滤波模块30的输出及自适应滤波器的系数，在数据传输模块10打包好之 后会被送到PC11中，并在波形图中显示。
[0071 ] h.观察自适应滤波器的系数，发现系数收敛较慢。通过PC11重新设置自适应滤波 器参数。PCI 1直接将运些参数送到FPGA模块9中的自适应滤波模块30，进行滤波器参数的修 改。
[0072] i.从步骤b开始到步骤h，重复进行多次实验，直到自适应滤波器的系数快速稳定 的收敛为止。固定滤波器参数进行不同深度水体探测实验，存储数据，用来做核磁信号的数 据反演。由实验数据算得核磁信号的信噪比提高3(WbW上。
[0073] 实施例2:
[0074] 在吉林大学中屯、校区做抗饱和自适应实时参考消噪实验。具体实施步骤如下：
[0075] a.铺设100m xlOOm的方形线圈作为发射线圈，W发射线圈的中屯、为基准，铺设25m x25m的探测线圈0,在距离基准35m左右的3个参考点铺设20应Im xlm的参考线圈1、2、3。
[0076] b.对整个系统上电，通过PC初始化各模块(程控放大器15、19、23、27的倍数设为最 大的16倍）。
[0077] C.在采集开始时，PC11发出采集的相关控制指令，经数据传输模块10送到FPGA模 块9中，由AD控制模块28产生同步信号，驱动继电器控制电路12、16、20、24与高精度40采集 模块8同步开始工作。
[0078] d.在继电器打开后，信号经过继电器控制电路12、16、20、24，到低噪声前放电路 13、17、21、25被放大，之后被宽频带带通滤波器14、18、22、26滤波放大，再通过程控放大器 15、19、23、27被再次放大。模拟信号(10(〇心1(〇心2(〇心3(〇通过高精度40采集模块8转 换成数字信号 dO(n)、Xl(n)、X2(n)、X3(n)。
[0079] e.FPGA模块9控制采集时钟频率并同步地将数字信号读入FPGA忍片中。经抗饱和 模块29判断信号饱和。抗饱和模块29自动调整程控放大倍数，到2倍时信号不饱和。
[0080] f.4个通道的不饱和信号同时进入自适应滤波模块30。采用图4中的变步长梯度自 适应格型联合算法对其进行自适应实时参考消噪。
[0081] g.自适应滤波模块30的输出及自适应滤波器的系数，在数据传输模块10打包好之 后会被送到PC11中，并在波形图中显示。
[0082] h.观察自适应滤波器的系数，发现系数不收敛。通过PC11重新设置自适应滤波器 参数。PC11直接将运些参数送到FPGA模块9中的自适应滤波模块30,进行滤波器参数的修 改。
[0083] i.从步骤b开始到步骤h，重复进行多次实验，直到自适应滤波器的系数快速稳定 的收敛为止。固定滤波器参数进行不同深度水体探测实验，存储数据，用来做核磁信号的数 据反演。由实验数据算得核磁信号信噪比提高35化W上。
【主权项】
1. 一种自适应滤波的磁共振信号抗饱和消噪装置，其特征在于，是由参考线圈a(l)经 前端处理模块a (4)与AD采集模块(8)连接，探测线圈（0)经前端处理模块b (5)与AD采集模块 (8) 连接，参考线圈b(2)经前端处理模块c(6)与AD采集模块(8)连接，参考线圈c(3)经前端 处理模块d(7)与AD采集模块(8)连接，前端处理模块a(4)经前端处理模块b(5)、前端处理模 块c (6)、前端处理模块d(7)和FPGA模块(9)分别与AD采集模块(8)和前端处理模块a(4)连 接，AD采集模块(8)经FPGA模块(9)、数据传输模块(10)和PC模块(11)与FPGA模块(9)连接构 成。2. 按照权利要求1所述的自适应滤波的磁共振信号抗饱和消噪装置，其特征在于，前端 处理模块a(4)由继电器控制电路a(12)经低噪声前放电路a(13)和带通滤波器a(14)与程控 放大器a(15)连接构成；前端处理模块b(5)由继电器控制电路b(16)经低噪声前放电路b (17)和带通滤波器b(18)与程控放大器b(19)连接构成;前端处理模块c(6)由继电器控制电 路c(20)经低噪声前放电路c(21)和带通滤波器c(22)与程控放大器c(23)连接构成;前端处 理模块d(7)由继电器控制电路d(24)经低噪声前放电路d(25)和带通滤波器d(26)与程控放 大器d(27)连接构成。3. 按照权利要求1所述的自适应滤波的磁共振信号抗饱和消噪装置，其特征在于，FPGA 模块(9)由AD控制模块(28)，抗饱和模块(29)和自适应滤波模块(30)组成。其中AD控制模块 (28)由PC(ll)控制产生同步信号;抗饱和模块(29)自动调节程控放大器a(15)、程控放大器 b (19)、程控放大器c (23)和程控放大器d(27)的放大倍数。4. 按照权利要求1所述的自适应滤波的磁共振信号抗饱和消噪装置的消噪方法，其特 征在于，包括以下步骤： a. 铺设一个正方形探测线圈（0)和三个多匝参考线圈a(l)、参考线圈b(2)、参考线圈c (3)； b. 系统上电，通过PC模块(11)初始化各模块； c. 在采集开始时，PC模块（11)发出采集控制指令，经数据传输模块（10)送到FPGA模块 (9) 中，由AD控制模块（28)产生同步信号，驱动继电器控制电路a(12)、继电器控制电路b (16 )、继电器控制电路c (20)和继电器控制电路d (24)与高精度AD采集模块(8)同步工作； d. 继电器打开后，信号经继电器控制电路a(12)、继电器控制电路b(16)、继电器控制电 路c(20)和继电器控制电路d(24)到低噪声前放电路a(13)、低噪声前放电路b(17)、低噪声 前放电路c(21)和低噪声前放电路d(25)放大，之后经宽频带带通滤波器a(14)、宽频带带通 滤波器b(18)、宽频带带通滤波器c(22)和宽频带带通滤波器d(26)滤波放大，再通过程控放 大器a(15)、程控放大器b (19)、程控放大器c (23)和程控放大器d(27)被再次放大，模拟信号 (10(〇41(〇42(〇和乂3(〇通过高精度40采集模块(8)转换成数字信号(10(11)41(11)、父2 (η)和X3(n); e . FPGA模块(9)产生采集时钟频率并同步地将数字信号读入FPGA芯片中，经抗饱和模 块(29)判断信号是否饱和，饱和时，FPGA模块会自适应的调节程控放大器a(15)、程控放大 器b(19)、程控放大器c(23)和程控放大器d(27)的倍数直到信号处于临界饱和状态； f. 抗饱和模块(29)的四个不饱和信号同时进入自适应滤波模块(30)，进行自适应实时 消噪处理； g. 自适应滤波模块(30)的输出及自适应滤波器的权系数，经数据传输模块（10)打包后 送到PC模块(11)中，并显示波形图； h.观察自适应滤波器的权系数，判断权系数是否快速稳定收敛，对不收敛、慢收敛或稳 态误差明显的，通过PC模块（11)重新设置自适应滤波器参数，PC模块（11)直接将这些参数 送到FPGA模块(9)中的自适应滤波模块(30)，进行滤波器参数的修改； i .修改好滤波器参数后，从步骤b到步骤h重复进行，直到权系数快速稳定的收敛为止， 固定滤波器参数，进行不同深度水体的探测，存储数据，并作核磁信号的数据反演。5.按照权利要求4所述的自适应滤波的磁共振信号抗饱和消噪装置的消噪方法，其特 征在于，参考线圈尽可能靠近干扰源铺设，参考线圈尽可能靠近干扰源铺设，参考线圈中心 到探测线圈中心的距离大于正方形探测线圈的边长。
【文档编号】G01V3/14GK105824053SQ201610344541
【公开日】2016年8月3日
【申请日】2016年5月23日
【发明人】林君, 舒旭, 张扬, 张洋, 林婷婷, 赵静
【申请人】吉林大学