摘要
文章基于谐振原理设计了一款手持式智能磁化率仪,推导了谐振原理法测量岩石磁化率计算公式,仿真探究谐振法测量岩石磁化率的影响因素,结合仿真结果与实际探测需要对仪器进行整体设计。基于LDC1614搭建外围电路,设计了4层螺旋线圈等距并联的电感传感器、信号提取程序和上位机操作界面,完成了磁化率仪样机制作。该仪器轻便易携(重量仅0.2 kg)、操作简单,测量结果与国外同类仪器比较误差低于5%,仪器精度可达1
磁化率是表征物质被磁化难易程度的物理
岩石磁化率测量方法主要有交流互感电桥法和自感
笔者基于LC谐振原理,设计了智能化手持式岩石磁化率仪。从测量原理出发分析了基于LC谐振原理的岩石磁化率计算方法,通过仿真探究影响岩石磁化率测量的因素,通过设计等距并联的4层螺旋电感线圈以减小仪器重量,基于LDC1614搭建硬件电路,设计了信号提取程序和上位机操作界面,保证了岩石视磁化率测量的准确性和一致性,进一步实现了仪器智能化和小型化,为后续设计提供了新思路。
在真空中磁感应强度与磁场强度的关系为
, | (1) |
式中:为真空磁导率;I为磁化电流;对于任意电感线圈,未放置磁介质时线圈电感为
。 | (2) |
对于LC并联传感器,其谐振频率[
, | (3) |
式中:为传感器电容。传感器电容将磁介质置于振荡状态的LC传感器内,磁介质会因为磁场的磁化作用而产生磁化强度,对于各向同性的物质,磁化强度与原磁场强度之比即为介质的磁化率,在国际单位制中(SI)三者关
。 | (4) |
此时受磁介质的影响,磁感应强度为:
, | (5) |
。 | (6) |
放置磁介质时线圈电感为
。 | (7) |
线圈电感会受填充系数影响,填充系数取决于线圈磁场与磁介质磁场耦合程度,考虑填充系数放置磁介质时电感量为
。 | (8) |
放置磁介质后电感传感器谐振频率为
。 | (9) |
由于岩石磁化率值比较小,对
。 | (10) |
岩石引起的线圈频率变化量为
。 | (11) |
磁介质磁化率为
。 | (12) |
根据
本研究的关键之处是岩石磁化率与线圈自感量的函数关系。利用COMSOL Mutiphysics仿真软件搭建电感磁化岩石模型,验证测量方法的可行性,并探究影响谐振法测量岩石磁化率的关键因素,根据仿真结果优化实验仪器设计。如

图1 岩石磁化模型示意图
Fig. 1 Schematic diagram of rock magnetization model
提离距离是线圈表面与岩样的垂直距离,它主要影响磁化场与待测样品的耦合程度。自然界中绝大部分岩石的相对磁导率为1~2,因此,仿真探究了线圈提离距离为1~10 mm时,岩石模型相对磁导率(1~2)与电感线圈电感量的关系。岩石模型电导率设置为1 S/m,激励频率设置为1 kHz,线圈与岩石模型比例设置为1∶1。仿真结果如

图2 提离距离1~10 mm时相对磁导率-线圈电感关系图
Fig. 2 Relationship between relative permeability and coil inductance at lifting distance of 1~10 mm
从
在交流电的磁化作用下,导电物体内会产生涡流效应,涡流效应会产生反磁场从而抵消岩石磁化场的作用,导致测量到的电感变化量不能真实反映岩石的磁化率。仿真分析岩石电导率(1

图3 不同电导率下相对磁导率-线圈电感关系图
Fig. 3 Relative permeability and coil inductance at different conductivities
由
在实际探测中,岩石的形状大小往往是不确定的。对于体积较小的岩石,要想准确评估岩石整体磁性,应将岩样完全置于磁场内。为便于分析,将岩石设置为正方体,探究边长为2~7 cm时相对磁导率与电感量的关系如

图4 相对磁导率-线圈电感关系
Fig. 4 Relative permeability and coil inductance
由
在测量形状较大的岩石或者巨型岩体某一截面的磁性时,由于岩体的面积远大于线圈的截面积,测量时只能将传感器贴近岩样。仿真研究时,岩石设置为正方体,将岩样放置在线圈外部,线圈半径设置为5 cm,探究岩样厚度为2 cm,边长分别为30、40、60、90 cm时(岩样面积远大于线圈面积)岩石磁导率与线圈电感的关系,岩石模型电导率设置为1 S/m,激励频率设置为1 kHz,提离距离设定为1 mm。仿真结果如

图5 岩体体积不同时相对磁导率-线圈电感关系图
Fig. 5 Coil inductance and relative permeability of rock mass of different sizes
由
基于FPGA设计了手持式磁化率仪测量系统,如

图6 仪器系统工作原理图
Fig. 6 Working principle of instrument system
硬件电路基于数字电感传感器LDC1614进行设计,此传感器工作原理如

图7 LDC1614工作原理图
Fig. 7 Working principle of LDC1614
芯片外接由探测线圈与谐振电容并联而成的LC振子,在测量操作指令下实时获取LC电路的谐振频率并以28位数字量反馈至上位机,电路中线圈电感量转换为数字量。在实际测量中不同种类的岩石磁性有差异,LC传感器的振荡频率随岩样磁性不同而变化,LDC1614输出值也因此不同,测量不同磁性岩样时传感器的振荡频率如

图8 不同磁性岩样磁化时振荡波形图
Fig. 8 Oscillations of magnetized rock samples with different magnetic properties
由
通过仿真可知,线圈电感与相对磁导率近似呈线性关系,要通过电感变化量反映岩石磁化率大小,应保证磁化线圈满足以下特征:
1)线圈能提供一定范围内的均匀磁场,磁场均匀度越好,测量结果越精确;
2)根据磁化率的测量原理,对同一样品,电感变化量越大,仪器的分辨率就越高,故应保证线圈电感变化量尽可能大;
3)为方便野外探测,该测量仪器应方便易携,传感器体积和重量应尽可能小。
基于以上要求,文中以多层平面电感线圈作为磁化线圈,平面线圈拥有体积小、质量轻,且易携带等优点,与电容配合时拥有较高的品质因素,LDC1614使用手册中也推荐使用平面电感。平面电感线圈设计借助于Webench Inductor Design Tool。根据实际需要设计线圈的层数、匝数、直径、线径宽度厚度和谐振电容大小等信息,软件将自动生成线圈的PCB原理图,给出线圈的电感量、谐振频率和电阻等信息。
为使线圈与谐振电容良好地配合,平面线圈可通过调节线圈匝数、层数得到所需的电感值。适当增大线圈层数,可在不改变线圈体积的情况下增大线圈自身电感量、提升线圈磁场强度。综合考虑,电感传感器采用4层平面线圈等距并联的方式,电感线圈绕线模型如

图9 4层螺旋线圈等距并联的平面电感线圈图
Fig. 9 Planar inductor coil with 4 helical coils equidistant in parallel
为探究4层平面电感线圈磁场特性及磁导率-线圈电感关系,采用有限元仿真软件ANSYS Maxwell建立4层螺旋线圈模型,模拟出4层电感线圈的磁场分布特性如

图10 4层平面电感线圈磁场特性图
Fig. 10 Load layout, section form and measuring points
由
基于FPGA设计了控制系统,控制系统的主要作用是控制和协调整个仪器的工作。工作流程如

图11 化率仪工作流程图
Fig. 11 Flowchart of the magnetic susceptibility meter
LDC1614芯片外接由探测线圈与谐振电容并联而成的LC振荡器,实时获取LC电路的谐振频率信息并以数字的形式存放在寄存器
。 | (13) |
式中:为寄存器中的数值;为分频系数。通过配置相应的寄存器可将芯片内部基准频率40 MHz进行分频,分频后参考频率
。 | (14) |
传感器初始谐振频率应结合线圈电感、寄生电容和谐振电容来确定。利用高精度谐振频率测试仪测试得到其准确值。传感器初始谐振频率还应保证当测试量程内最大磁化率时寄存器中数据不会溢出。根据实验测试结果,设置芯片时钟分频寄存器的分频模式。
传感器数据转换所需要的时间由名为RCOUNT的寄存器决定,令其中数值为,则为
。 | (15) |
传感器测量过程激励所需要的时间由名为SETTLECOUNT的寄存器决定,令其中数值为g,则为
。 | (16) |
为保证数据稳定,将寄存器RCOUNT中值设置为0xFFFF,SETTLECOUNT中的值设置为0x0100。
寄存器DRIVE_CURRENT中数值决定线圈驱动电流大小,LDC1614的激励电流范围为16 µA~1.6 mA,电流越大磁场越强,信号分辨率就越好,但电流过大易造成磁饱和。在设置驱动电流时应对仪器量程内的最大磁化率与最小磁化率样品进行实验,保证在最恶劣的情况下传感器仍能维持稳定的振荡。
电感变化量表征的是岩样的整体磁性,为保证测量的准确性,针对不同形状、大小的岩石建议更换不同规格的线圈,保证线圈与岩样完全贴合。对于磁化线圈,选用相对磁导率已知的样品对仪器进行标定,绘制出岩样磁化率与谐振频率之间的关系曲线,并将关系曲线储存在测量程序中,下次使用时可根据之前储存的关系曲线,推测并显示未知样品的磁化率。
设计了磁化率仪上位机,操作程序可支持电脑和平板使用。如

图12 磁化率仪工作界面图
Fig. 12 Working interface of the susceptibility meter
为验证磁化率仪的测量效果,选用不同种类的岩石标本,包含火成岩、变质岩、沉积岩共计216种,将测量结果与实际工程中应用较为广泛的SM-30、KM-7磁化率仪进行比较,部分岩样测量结果如
编号 | 岩石 种类 | SM-30磁化率仪结果 | KM-7磁化率仪结果 | 文中方法测量结果 | |||
---|---|---|---|---|---|---|---|
一次测量 | 二次测量 | 一次测量 | 二次测量 | 一次测量 | 二次测量 | ||
1 | 磁铁矿 | 801.000 | 794.000 | 794.000 | 799.000 | 803.800 | 805.600 |
2 | 钛铁矿 | 424.000 | 426.000 | 414.000 | 412.000 | 425.600 | 424.100 |
3 | 绿帘石 | 130.000 | 125.000 | 128.000 | 126.000 | 128.700 | 127.100 |
4 | 辉绿岩 | 47.200 | 45.800 | 46.800 | 45.500 | 45.200 | 46.800 |
5 | 辉岩 | 10.100 | 10.600 | 9.800 | 9.700 | 9.500 | 9.900 |
6 | 玄武岩 | 3.180 | 3.500 | 3.400 | 3.800 | 2.700 | 2.800 |
7 | 珍珠岩 | 0.780 | 0.850 | 0.590 | 0.650 | 0.860 | 0.890 |
8 | 明矾石 | 0.070 | 0.050 | 0.080 | 0.060 | 0.050 | 0.060 |
9 | 页岩 | 0.020 | 0.010 | 0.030 | 0.020 | 0.010 | 0.010 |
10 | 砾岩 | 0.001 | 0.002 | 0.003 | 0.004 | 0.002 | 0.003 |
注: 表中数据需乘以,单位为国际单位SI。
1)根据岩石磁化率测量要求,提出了基于LC谐振原理的岩石磁化率测量方法,分析并推导了运用此方法计算岩石磁化率的表达式。
2)搭建仿真模型,验证了测量方法的正确性,并分析了提离距离、岩石电导率、岩样大小等因素对测量的影响,明确了仪器设计的优化条件和实际测量时的注意事项,为仪器具体设计提供了支持。
3)基于FPGA设计了仪器整体测量系统,选用高精度数字电感传感器LDC1614搭建了测量的硬件电路,提出4层电感线圈等距并联的传感器设计方案,极大地提高了测量精度(可达1
4)设计仪器整体工作流程,即空气-岩石-空气3次测量方式,极大地减小了温度等不确定因素造成的误差,完成了仪器控制系统、软件部分、上位机设计及样机制作。实测结果表明本仪器测量误差小、精度高、测量一致性好,且仪器操作简单,适合野外作业。
参考文献
Mitchell J. Monitoring lithology variations in drilled rock formations using NMR apparent magnetic susceptibility contrast[J]. Applied Magnetic Resonance, 2020, 51(3): 205-219. [百度学术]
Nyam G G, Adeeko T O, Umar M, et al. Assessment of magnetic susceptibility of some selected rock samples from Karu area, northcentral Nigeria[J]. Asian Journal of Advanced Research and Reports, 2019: 1-6. [百度学术]
Yan J Y, Lü Q T, Luo F, et al. A gravity and magnetic study of lithospheric architecture and structures of South China with implications for the distribution of plutons and mineral systems of the main metallogenic belts[J]. Journal of Asian Earth Sciences, 2021, 221: 104938. [百度学术]
Tan M J, Fan L J, Mao K Y, et al. Influential factors analysis and porosity correction method of nuclear magnetic resonance measurement in igneous rocks[J]. Journal of Applied Geophysics, 2019, 161: 153-166. [百度学术]
Roy S. A simple bridge to operate an AC susceptometer[J]. European Journal of Physics, 2020, 41(5): 055501. [百度学术]
Li T, Shang H P, Zhang Q Y, et al. Optimizing nonlinearity of piezoresistive pressure sensors by an asymmetric Wheatstone bridge[J]. Microsystem Technologies, 2021, 27(8): 2989-2994. [百度学术]
詹文山. 交流磁化率的测量[J]. 物理, 1982(12):732-736. [百度学术]
Zhan W S. The measurement of AC susceptibility[J]. Physics, 1982(12):732-736. (in Chinese) [百度学术]
裴朝, 崔益民, 李华, 等. 提高锁相放大器测量交流磁化率精度的方法[J]. 物理实验, 2010, 30(11): 38-40. [百度学术]
Pei Z, Cui Y M, Li H, et al. Improving accuracy of AC susceptibility measurement with lock-in amplifier[J]. Physics Experimentation, 2010, 30(11): 38-40. (in Chinese) [百度学术]
谷静, 杨雁, 陆青, 等. 基于数字比例技术的高精度交流电桥研究[J]. 仪器仪表学报, 2020, 41(7): 29-37. [百度学术]
Gu J, Yang Y, Lu Q, et al. Research on high precision AC bridge based on digital ratio technique[J]. Chinese Journal of Scientific Instrument, 2020, 41(7): 29-37. (in Chinese) [百度学术]
王力波. HKB-1型磁化率仪研制及计算机联机的实现[J]. 火山地质与矿产, 1991, 12(4): 40-52. [百度学术]
Wang L B. Model HKB-1 magnetic susceptibility meter and connected with computer[J]. Volcanology & Mineral Resources, 1991, 12(4): 40-52. (in Chinese) [百度学术]
Talebi A, Hosseini S V, Parvaz H, et al. Design and fabrication of an online inductive sensor for identification of ferrous wear particles in engine oil[J]. Industrial Lubrication and Tribology, 2021, 73(4): 666-675. [百度学术]
查忆秋, 肖锋, 黄为俊, 等. SM-30磁化率仪的工作原理探讨[J]. 工程地球物理学报, 2017, 14(1): 45-49. [百度学术]
Zha Y Q, Xiao F, Huang W J, et al. The operating principle of SM-30 magnetic susceptibility instrument[J]. Chinese Journal of Engineering Geophysics, 2017, 14(1): 45-49. (in Chinese) [百度学术]
叶春, 陈晔, 陈仕涛. 智能化磁化率仪原理及研制[J]. 南京师范大学学报(工程技术版), 2001, 1(1): 55-59, 63. [百度学术]
Ye C, Chen Y, Chen S T. Study of intelligent magnetic susceptibility instrument[J]. Journal of Nanjing Normal University (Engineering and Technology), 2001, 1(1): 55-59, 63. (in Chinese) [百度学术]
周锡华, 孟祥聪, 王金龙, 等. 智能磁化率仪的研制[J]. 物探与化探, 2005, 29(5): 428-430. [百度学术]
Zhou X H, Meng X C, Wang J L, et al. The developing of intelligent magnetic susceptibility meter[J]. Geophysical and Geochemical Exploration, 2005, 29(5): 428-430. (in Chinese) [百度学术]
王振宇, 成立. 模拟电子技术基础[M]. 3版. 南京: 东南大学出版社, 2019. [百度学术]
Wang Z Y, Cheng L. Fundamentals of analog electronic technology[M]. 3rd ed. Nanjing: Southeast University Press, 2019. (in Chinese) [百度学术]
刘秉裕. 矿石矿物磁选[M]. 北京: 冶金工业出版社, 2020. [百度学术]
Liu B Y. Ore mineral magnetic separation[M]. Beijing: Metallurgical Industry Press, 2020. (in Chinese) [百度学术]
Wang X, Lv D, Sun L, et al. Magnetic behaviors of a ferrimagnetic decorated kagome-like lattice under an external magnetic field[J]. Journal of Magnetism and Magnetic Materials, 2021,538(5): 168259. [百度学术]
胡瑞霞, 马卫华, 史天成, 等. 基于LDC1614的磁浮电涡流位移传感器设计[J]. 传感器与微系统, 2019, 38(7): 67-69. [百度学术]
Hu R X, Ma W H, Shi T C, et al. Design of maglev eddy current displacement sensor based on LDC1614[J]. Transducer and Microsystem Technologies, 2019, 38(7): 67-69. (in Chinese) [百度学术]
Li C, Zhang W P, Zhang Q, et al. Research on the detection of curved ferromagnetic materials by flexible eddy current array sensor based on LDC1614[J]. Journal of Physics: Conference Series, 2020, 1544(1): 012135. [百度学术]