磁铁矿矿致磁异常测定方法及其应用
A novel method for determining magnetite ore-induced magnetic anomalies and its application
通讯作者: 何敬梓(1990-),男,博士,高级工程师,主要从事综合地球物理勘查研究工作。Email:749079407@qq.com
第一作者:
收稿日期: 2025-09-16 修回日期: 2025-10-14
| 基金资助: |
|
Received: 2025-09-16 Revised: 2025-10-14
准确、快速地确定磁异常源的地质属性是磁法勘探的关键技术之一,直接影响到利用磁测数据开展地质解释的准确性。感应磁化强度随时间变化,剩余磁化强度一般不随时间变化,因此,理论上可以通过探测磁异常强度随时间变化特征来研究磁异常场源体性质。虽然前人开展过相关研究,但目前尚未形成实用的技术方法。本文提出利用地磁日变观测数据,构建磁异常强度变化量A参数、磁异常强度变化率η参数和归一化磁异常强度变化量F参数来评估强磁异常场源体性质,确定磁异常是否为磁铁矿引起的新方法,该方法对磁铁矿勘探具有重要意义。
关键词:
Accurately and quickly determining the geological properties related to magnetic anomaly sources is a key technical challenge in magnetic prospecting, directly influencing the accuracy of geological interpretation using magnetic survey data. The induced magnetization varies with time, whereas the remanent magnetization typically remains constant over time. Therefore, there exists a theoretical basis for investigating the properties of magnetic anomaly source bodies by detecting the time variations of magnetic anomaly intensity. Despite geophysicists' relevant research in this field, practical technical methods have not been established. Hence, this study proposed a novel method for determining magnetite ore-induced magnetic anomalies. Based on the observational data of geomagnetic diurnal variations, the proposed method constructed parameters, including the variations (A), variation rate (η), and normalized variations (F) of magnetic anomaly intensity, to evaluate the properties of strong magnetic anomaly source bodies. Accordingly, the proposed method determined the possibility of magnetite ore-induced magnetic anomalies, showing critical significance for magnetite exploration.
Keywords:
本文引用格式
范正国, 杨海, 葛藤菲, 何敬梓, 贾志业, 范振宇, 刘前坤, 杨雪.
FAN Zheng-Guo, YANG Hai, GE Teng-Fei, HE Jing-Zi, JIA Zhi-Ye, FAN Zhen-Yu, LIU Qian-Kun, YANG Xue.
0 引言
截至2021年,我国利用航磁测量结果筛选的磁异常有5万多处,仅对其中1万余处磁异常开展过查证工作,发现各类矿产1 300多处(大部分为铁矿),其中大型及以上规模的矿床273处[8-9],受限于地质认识和磁异常查证手段的不足,目前尚存大量的磁异常性质不明确。已发现的这些磁异常是否具有找铁矿的潜力,需要利用一些经济、高效的磁异常查证新方法,对其进行快速评价。另外,地表矿已基本发现殆罄,找矿工作正向隐伏区、甚至是深部进军,一般的磁异常查证方法难以直接确定磁异常源的地质属性[3,10-
在前人研究基础上,笔者通过理论研究及在鲁西邹平火山岩区和金岭磁铁矿区的试验,发明了快速确定磁异常源地质属性的简便方法,对磁异常的快速评价和磁铁矿找矿突破具有重要意义。本文简要介绍该项方法技术及其试验结果,供磁法勘探工作者参考。
1 基本原理
图1
图1
山东省平阴地区地磁场随时间变化曲线(据参考文献[25])
Fig.1
Time variations of geomagnetic field intensity in the Pingyin area, Shandong Province
磁异常强度变化量A定义为磁异常区某处的异常值在两个时间点t0、t的强度差值,计算公式为
背景磁场变化幅度δTb和异常磁场变化幅度δTa为背景站磁场强度和异常站磁场强度在两个时间点t0、t的差值,计算公式为
异常站磁场变化幅度大于背景站磁场变化幅度的百分比,简称为磁异常强度变化率η,计算公式为
由式(1)~式(5)可知,磁异常强度变化量A和磁异常强度变化率η不仅与磁异常源的磁性(Mi、Mr)有关,还与地磁场强度及其随时间变化幅度、异常值等因素有关,无法根据其大小直接判断磁异常源的地质属性。为此,引入一个新参数F,定义为t时刻的归一化磁异常强度变化量[28]。通过归一化处理,F参数成为一个仅与磁异常源的磁性(Mi、Mr)有关、而与地磁场强度及其随时间变化幅度和异常值等无关的量,具体关系为:Mi≪Mr时F≈0,Mi≈Mr时F≈0.5,Mi≫Mr时F≈1。因此,根据F参数的数值大小,可以确定Mi和Mr的相对大小,即确定磁异常源的地质属性。
2 地质与地球物理概况
2.1 区域地质背景
试验区位于华北克拉通东部的鲁西地块东北缘、齐河—广饶断裂南侧,存在多期构造—热液事件的叠加,是我国重要的矽卡岩型富铁矿成矿带[29-
图2
试验区矿产主要为磁铁矿,分布于金岭杂岩体周边,是山东省重要的矽卡岩型富铁矿集中发育区之一,累计探明富铁矿资源储量超过2亿t[38],2 000 m 以浅预测资源超过5亿t[39]。金岭杂岩体主要岩性包括辉长闪长岩、角闪石闪长岩和二长岩,围岩主要为中奥陶统的灰岩和石炭系—二叠系的砂岩、页岩。矽卡岩型铁矿床(点)在空间上与金岭杂岩体紧密伴生[35],分布于金岭杂岩体与碳酸盐岩地层接触带附近,已发现大、中、小矿床十余处,其中王旺庄磁异常反映的铁矿,目前已控制铁矿储量超过1亿t;侯家庄磁异常反映的铁矿,目前已控制铁矿储量超过2 000万t;北金召磁异常反映的铁矿,目前已控制铁矿储量约2 400万t[40-41]。
2.2 地球物理特征
2.2.1 岩(矿)石磁性特征
表1 邹平磁异常区和金岭磁异常区磁性特征及异常强度
Table 1
2.2.2 航磁异常特征
图3
图3
山东省邹平—金岭地区航磁ΔT等值线及磁日变站分布示意
Fig.3
Aeromagnetic anomaly map and distribution of magnetic diurnal variation stations in the Zouping-Jinling area, Shandong Province
金岭磁异常区为NNE向展布的块状异常,其中叠加了多处孤立的等轴状局部异常;王旺庄磁异常位于金岭异常块东北端,侯家庄、北金召磁异常位于金岭异常块北侧边缘(图3)。金岭磁异常区于1959年1∶10万航空物探(磁及放射性)测量时首次发现,表现为负背景场中的高值异常,整体呈轴向NE的椭园状,内部升高次级异常较多,部分为强度可达+1 400 nT之尖峰状次级异常[36]。1979年,通过1∶5万航空磁力和放射性测量,从金岭异常中分解出侯家庄异常(60-4-1)、张庄异常(60-4-2)、北金召异常(60-4-4)、王旺庄异常(60-4-5)等,钻探结果证实这些局部叠加异常都是由矽卡岩型铁矿所引起[42]。2020年,在金岭磁异常区东北部开展了1∶1万航空磁力测量,涉及了侯家庄异常(60-4-1)、北金召异常(60-4-4)和王旺庄异常(60-4-5),更加精细地刻画了局部异常特征,其中王旺庄异常峰值达1 520 nT、北侧伴生-1 075 nT异常,侯家庄异常峰值达880 nT、北侧伴生-810 nT异常,北金召异常峰值达1 518 nT、北侧伴生-750 nT异常。
3 应用试验
2024年10月31日~11月1日,在山东省邹平磁异常区和金岭磁异常区开展了本次新方法技术的应用试验,背景站位于金岭磁异常区西侧的平静磁场区,异常站分别位于邹平磁异常区内的孟白庄局部磁异常的极大值附近,以及金岭磁异常区内的王旺庄和北金召局部磁异常的极大值附近、侯家庄局部磁异常的极小值附近(图3)。
本次野外磁场测量使用了2台Overhauser型 GSM-19磁力仪和2台Overhauser型GSM-19T磁力仪,其中GSM-19磁力仪的灵敏度为0.015 nT,分辨率为0.01 nT,绝对精度为±0.1 nT;GSM-19T磁力仪的灵敏度为0.05 nT,分辨率为0.01 nT,绝对精度为±0.2 nT。在野外施工前,对使用的磁力仪等设备进行了检测,确保它们均处于正常工作状态,而且不同设备都进行了测量时间同步。磁场测量的采样率为1次/2 min,测量数据处理结果列于表2。
表2 孟白庄、王旺庄、侯家庄和北金召4处磁异常的测量数据处理结果
Table 2
| 异常名称 | 地质属性 | 地磁场强度/nT | 观测点异常值/nT | 背景站磁场 变化幅度/nT | 异常站磁场 变化幅度/nT | A/nT | η/% | F |
|---|---|---|---|---|---|---|---|---|
| 王旺庄 | 磁铁矿 | 52850 | 824.8 | 35.62 | 36.11 | 0.49 | 1.38 | 0.88 |
| 侯家庄 | 磁铁矿 | 52850 | -184.3 | 44.43 | 44.29 | -0.14 | -0.32 | 0.91 |
| 北金召 | 磁铁矿 | 52850 | 481.7 | 34.91 | 35.12 | 0.21 | 0.60 | 0.66 |
| 孟白庄 | 火山岩 | 52850 | 979.0 | 38.28 | 38.60 | 0.32 | 0.85 | 0.45 |
下面详细介绍具有代表性的孟白庄、王旺庄和侯家庄3处磁异常的测量数据处理过程,其中孟白庄地区为火山岩区,王旺庄和侯家庄为磁铁矿区;孟白庄和王旺庄异常站处的磁异常值为正,侯家庄异常站处的磁异常值为负。
3.1 孟白庄磁异常应用试验
2024年10月31日开展了磁异常背景站和孟白庄异常站的磁场测量,结果(图4)显示,时间从11:54:02到16:24:02,地磁场强度持续增加,背景站的磁场值从52 834.24 nT增加到了52 872.52 nT,增幅达38.28 nT。背景站和异常站的磁场曲线较光滑且变化趋势一致,表明背景站和异常站的磁场变化规律基本一致,数据中噪声较小。
图4
图4
孟白庄磁异常背景日变站(Tb)和异常日变站(Ta)磁日变曲线
Fig.4
Curves showing the magnetic diurnal variations at the background(Tb) and anomaly stations of the Mengbaizhuang magnetic anomaly(Ta)
根据F参数值,孟白庄磁异常源的磁性特征为剩余磁化强度大于感应磁化强度,与孟白庄磁异常区岩矿标本的磁性特征一致。
3.2 王旺庄异常应用试验
2024年10月31日开展了磁异常背景站和王旺庄异常站的磁场测量,结果(图5)显示,时间从11:54:02到16:24:02,地磁场强度持续增加,背景站和异常站的磁场变化趋势基本一致,但背景站的磁场曲线较光滑,而异常站的磁场曲线波动较明显,表明背景站的磁场数据中噪声较小,异常站的磁场数值中噪声较大。
图5
图5
王旺庄磁异常背景日变站(Tb)和异常日变站(Ta)磁日变曲线及总噪声(S)
Fig.5
Curves showing the magnetic diurnal variations at the background(Tb) and anomaly stations of the Wangwangzhuang magnetic anomaly(Ta)
为了进一步量化噪声分布特征,将计算出的背景站和异常站磁场日变数据四阶差分的平方根S作为总噪声,结果表明大部分数据中的总噪声S>80 pT(图5中仅显示了S≤80 pT部分,S>80 pT时被限幅)。
根据F参数值,王旺庄磁异常源的磁性特征为剩余磁化强度小于感应磁化强度,与金岭磁异常区铁矿石标本的磁性特征一致。
3.3 侯家庄异常应用试验
2024年11月1日开展了磁异常背景站和侯家庄异常站的磁场测量,结果(图6)显示,时间从12:10:02到16:08:02,地磁场强度持续增加,背景站的磁场值从52 828.33 nT增加到了52 872.76 nT,增幅达44.43 nT。背景站和异常站的磁场曲线较光滑、变化趋势一致,表明背景站和异常站的磁场变化规律基本一致,数据中噪声较小。
图6
图6
侯家庄磁异常背景日变站(Tb)和异常日变站(Ta)磁日变曲线
Fig.6
Curves showing the magnetic diurnal variations at the background(Tb) and anomaly stations of the Houjiazhuang magnetic anomaly(Ta)
根据F参数值,侯家庄磁异常源的磁性特征为剩余磁化强度小于感应磁化强度,与金岭磁异常区铁矿石标本的磁性特征一致。
4 讨论
4.1 用磁异常强度变化量及其变化率判断磁异常源地质属性的效果
表3 磁异常强度变化量A值及磁异常强度变化率η
Table 3
以往学者们开展过零星相关研究,但均未形成系统性成果。苏联地磁学家Б·М·扬诺夫1938年在含铁石英岩(库尔斯克磁异常)附近观测到地磁场日变和其他长周期变化幅度非常之大[17]。美国地球物理学家Ward S.H.和研制仪器厂家于20世纪60年代初在一个小型致密块状磁铁矿体露头附近观测磁场日变,发现周期为20~30 s的地磁场短周期微扰内,磁异常强度变化率η>30%[13,17]。随后, Goldstein N. E.和Ward S. H.在美国西部的2个磁性铁矿体(磁化率分别为75 398×10-5 SI和125 664×10-5 SI,Q为0.06~0.7)附近和2处有火山岩(一处为火山玻璃,其剩余磁化强度大约为200 A/m, 另一处为熔结凝灰岩)分布区观测磁场日变,发现Pc3微扰内磁性铁矿体磁异常强度变化量A分别达到0.035 nT和0.2 nT,火山岩区磁场随时间变化曲线与正常场随时间变化曲线完全一致,Pc3微扰的周期范围为10~45 s,地磁场变化幅度约0.5~1.0 nT[14,17]。20世纪70年代,冯剑萍在小兴安岭一个性质未知的磁异常区和背景磁场区进行了太阳静日变化(Sq)磁场日变观测,未发现明显差异,认为该磁异常不是由主要为感应磁化的磁铁矿体引起[17]。20世纪80~90年代,澳大利亚地球物理学家Clark D. A.和Schmidt P. W.等也进行了区分磁性体的感应磁化与剩余磁化的试验工作,认为原地确定磁性体磁化特征的方法正确有效[15,17]。这些成果表明,磁异常强度变化量A及其变化率η在磁铁矿区的值较大,在火山岩区的值较小。
本次应用试验中,王旺庄、侯家庄、北金召3个磁铁矿区的磁异常强度变化量A分别为0.49 nT、-0.14 nT和0.21 nT,磁异常强度变化率η分别为1.38%、-0.32%和0.6%;孟白庄火山岩区的A为0.32 nT,η为0.85%。对比可知,这4处试验区中,1处磁铁矿区的A、η值较大,2处磁铁矿区的A、η值较小,火山岩区居中。
精度应分别达到或优于0.05 nT、0.01 nT和±0.2 nT,应具有测量精度(四阶差分)优于±80 pT的大量数据。
4.2 用F参数判断磁异常源地质属性的效果
F参数,即归一化磁异常强度变化量,是笔者首次提出的一种表征磁异常强度变化特征的参数[28]。
本次应用试验过程中,在王旺庄、侯家庄及北金召磁铁矿区和孟白庄火山岩区进行了磁日变观测。计算结果:磁铁矿区的F参数分别为0.88、0.91、0.66,火山岩区F=0.45。前面已提到,Mi≪Mr时F≈0,Mi≈Mr时F≈0.5,Mi≫Mr时F≈1,因此,王旺庄、侯家庄、北金召3个磁铁矿区以感应磁化强度Mi为主,孟白庄火山岩区以剩余磁化强度Mr为主。
据岩(矿)石标本磁性测量结果(表1),王旺庄、侯家庄和北金召磁铁矿区的Q值为0.01~0.52,表明Mi>Mr;孟白庄火山岩区的Q值为35.7~178,表明Mi<Mr。所以,用磁日变数据的计算结果与岩矿石标本的测量结果一致。
综上分析,用F参数可以准确判断磁异常源是感应磁化强度Mi为主还是剩余磁化强度Mr为主。
4.3 用F参数判断铁矿矿致异常的可能性
前文已提到,用F参数可以准确判断磁异常源的磁化强度性质。实践表明,磁铁矿,尤其是埋藏较浅的磁铁矿,一般可以引起强度大、梯度陡的局部磁异常。因此,根据磁异常特征,结合F参数,可以辅助判断磁异常是否属于磁铁矿矿致异常。
5 结论与建议
通过观测磁异常区和背景地区的磁场随时间变化来确定磁异常源磁性特征的研究已持续了近90年,并取得了一定结果。在本研究中,笔者首次提出了F参数,并在山东省邹平磁异常区和金岭磁异常区进行了应用试验。对磁异常强度变化量A及其变化率η和F参数试验结果进行比较,得出了以下结论:
1)本技术既可利用太阳静日变化(Sq)的磁日变,也可利用磁暴、微扰等的磁日变。
2)虽然在大多数情况下能用磁异常强度变化量A及其变化率η大致判断磁异常源是以感应磁化强度Mi为主还是以剩余磁化强度Mr为主,但没有量化的判断标准,仅为定性分析。
3)可以用F参数准确判断磁异常源是以感应磁化强度Mi为主还是以剩余磁化强度Mr为主。
4)可以用F参数辅助判断磁异常是否属于磁铁矿矿致异常。
此外,本文所提方法技术的实际应用尚不多,其有效性和适应性尚需在后续工作中进一步验证。未来将持续研究、完善该方法的技术及其工作流程,促使其成为一种经济、高效的磁异常快速评价和识别铁矿矿致磁异常的新方法技术。
参考文献
解决物探异常解释多解性的一次尝试
[J].
An attempt to reduce ambiguity in geophysical interpretation
[J].
航空地球物理在战略性矿产勘查中的应用前景
[J].
Application prospect of aerogeophysics in strategic mineral exploration
[J].
New insights into the crustal magnetization of the eastern Indian shield inferred from the aeromagnetic data-based depth to the bottom of the magnetic sources (DBMS)
[J].
Crustal structure across the central part of Narmada-Son Lineament,India based on the interpretation of aeromagnetic and gravity data:Geological implications
[J].
秦岭及天山等重点成矿区带航空物探调查成果报告
[R].
Aerial geophysical survey of key mineralization belts such as Qinling and Tianshan Mountains
[R].
邯邢式铁矿深部探测技术及综合找矿模式研究——以河北省沙河市白涧铁矿床为例
[J].
Application of deep detecting technology in Hanxing subtype iron deposits and an integrated prospecting model:A case study of Baijian iron deposit
[J].
厚覆盖区磁铁矿空—地—井协同勘查技术体系
[J].
Airborne-surface-borehole cooperative exploration technical system for magnetite exploration in areas with thick overburden
[J].
基于反射波地震勘探技术的深覆盖区矽卡岩型富铁矿找矿方法探讨
[J/OL].
Discussion on prospecting methods for rich iron ore in deep coverage areas based on reflected wave seismic exploration technology
[J/OL].
A field experiment with a rubidium-vapor magnetometer
[J].
The separation of remanent from induced magnetism in situ
[J].
Remote determination of magnetic properties and improved drill targeting of magnetic anomaly sources by differential vector magnetometry (DVM)
[J].
关于铁矿物探工作的几点建议
[EB/OL].(
Some suggestions on geophysical exploration of iron minerals
[EB/OL].(
磁异常解释中的剩余磁化问题
[J].
Remanent magnetization in the interpretation of magnetic anomaly
[J].
提高磁测日变改正精度的方法
[J].
The method of increasing precision of diurnalcorrection in high-precision magnetic survey
[J].
两种典型磁扰对航空磁测的影响
[J].
The influence of two typical kinds of magnetic disturbance on the airborne magnetic survey
[J].
地球磁场的变化及磁暴现象的原因分析
[J].
Changes of earth’s magnetic field and cause analysis of magnetic storm phenomenon
[J].
山东齐河—禹城地区1∶5万航空物探(磁)测量成果报告
[R].
1∶50000 aerial geophysical (gravimetric and magnetic) survey in Qihe-Yucheng area,Shandong Province
[R].
鲁西金岭地区含矿闪长岩体成因:来自锆石U-Pb年代学和地球化学证据
[J].
Petrogenesis of ore-bearing dioritic pluton in jinling area in western Shandong:Evidence from zircon U-Pb chronology and petro-geochemistry
[J].
Petrogenesis of the zhangmatun gabbro in the ji’nan complex,North China Craton:Implications for skarn-type iron mineralization
[J].
华北克拉通矽卡岩型富铁矿成矿关键控制因素:来自地层学的制约
[J].
Key ore-controlling factors of the skarn type high-grade iron deposits in North China Craton:Constraints from stratigraphy
[J].
鲁西地区中生代高镁闪长岩角闪石成分变化特征及其对矽卡岩型铁矿成矿作用的启示
[J].
Compositional variations of amphiboles from Mesozoic high-Mg diorite in western Shandong:Implications for the mineralization of skarn-type iron deposits
[J].
Crust-mantle interaction during the tectono-thermal reactivation of the North China Craton:Constraints from SHRIMP zircon U-Pb chronology and geochemistry of Mesozoic plutons from western Shandong
[J].
Two tales of the continental lithospheric mantle prior to the destruction of the North China Craton:Insights from Early Cretaceous mafic intrusions in western Shandong,East China
[J].
Genetic relationship of high-Mg dioritic pluton to iron mineralization:A case study from the Jinling skarn-type iron deposit in the North China Craton
[J].
山东中部及东部地区航空物探(磁及放射性)工作结果报告
[R].
Aerogeophysical (magnetic and radiometric) survey results in central and eastern Shandong Province
[R].
鲁西中生代高Mg闪长岩的成因:年代学与岩石地球化学证据
[J].
Petrogenesis of the Mesozoic high-Mg diorites in west Shandong:Evidence from chronology and petro-geochemistry
[J].
山东淄博金岭地区辛庄铁矿地质地球物理特征及深部找矿预测
[J].
Geological and geophysical characteristics and deep prospecting prediction of Xinzhuang iron deposit in Jinling area,Zibo,Shandong
[J].
山东淄博金岭杂岩体分布区深部富铁矿资源潜力及找矿方向
[J].
Iron-rich ore exploration potential and favorable sites for deep horizon of Jinling magmatic complex distribution in Zibo,Shandong Province
[J].
鲁中隆起区中北部矽卡岩型铁矿成矿预测
[J].
Metallogenic prediction of skarn iron deposits in the central and northern Luzhong uplift,East China
[J].
鲁中地区航空磁力和放射性测量成果报告
[R].
Aeromagnetic and radiometric survey results in central Shandong region
[R].
大兴安岭地区的强磁性火山岩
[J].
On the magnetic property of volcanic rocks occurring in the great shingan mountain region,Northeast China
[J].
/
| 〈 |
|
〉 |
