砂岩型铀矿航磁微弱异常提取方法

杨玉勤,张翔,石连成,邓德伟

(1.核工业航测遥感中心，河北 石家庄 050002;2.河北省航空探测与遥感技术重点实验室，河北 石家庄 050002)

0 引言

砂岩型铀矿的航磁微弱信息研究一直是航空物探领域重要课题之一，前人进行过多项研究工作[1-6]，认为多数砂岩型铀矿上方有微弱航磁异常。但在中—新生代盆地区，受盆地盖层影响，砂岩型铀矿的航磁异常不明显，直接影响了利用航磁资料在中—新生代盆内寻找砂岩型铀矿的应用效果，需寻找总结有效提取航磁微弱异常信息的方法。在我国北方中—新生代盆地，已开展了多年以寻找砂岩型铀矿为目的的高精度航磁、航放综合测量工作，依据最新的航磁资料，进一步开展砂岩型铀矿的航磁微弱信息研究，对发挥高精度航磁资料在砂岩型铀矿勘查中的作用具有非常重要的意义。

1 研究依据

许多研究者曾开展过利用航磁微(弱)磁异常寻找油气和金属矿产等工作，进行了多种提取与油气田、金属矿及构造有关的微(弱)磁异常的方法研究[7-13]。

在铀矿勘探领域，某些铀矿(砂岩型和热液型)在形成与富集过程中,可使其周围铁磁性矿物遭受破坏或交代,使矿化段呈现弱磁性的规律。为此，国内外很多研究者在砂岩型铀矿的航磁微弱磁异常的形成机理和提取方法等方面进行了多项研究[1-6,14-17]。俄罗斯研究者于1995～1998年，在古河床内与透水层中的层间氧化带尖灭界线有关的巴尔克夫铀矿床上研究发现，高精度航磁ΔT出现幅值为±3～±12 nT的高频异常[18]。国内学者在研究伊犁盆地砂岩型铀矿航磁特征时，发现在3个已知铀矿床上方均存在清晰的航磁弱异常，这些弱磁异常叠加在宽缓负异常背景场上，幅值为5～25 nT[1-2]。通过对比可地浸砂岩型铀矿(伊犁某矿)氧化带、氧化还原过渡带和还原带土壤铁磁性矿物含量的变化(图1)，对砂岩型铀矿上方微磁异常形成的机理进行了探讨，认为氧化还原过渡带地表土壤铁磁性矿物含量明显增高，使相应地段磁性增强，从而在砂岩型铀矿的氧化还原过渡带(铀矿体)范围内产生微弱磁异常。

图1 伊犁某矿床铀矿体地球化学分带与地表土壤含铁矿物分布[1]Fig.1 The map of geochemical zoning of ore bodies and distribution of iron minerals in surface soil of Yili uranium deposit[1]

2 航磁微弱磁异常的界定

虽然在利用航磁弱异常区寻找油气和砂岩型铀矿方面已开展了很多研究工作，但现有的规范和资料并没有对航磁微弱磁异常进行严格定义和量化。现行航空磁测技术规范(DZ/T 0142-2010)中在规定航磁局部异常的确定原则时，规定局部磁场变化大于或等于测量总精度的2倍，并且与一定的地质现象联系时，即为局部异常；同时还要求对小于测量总精度的2倍，但与地质任务关系密切的“弱异常”也应加以分析。航磁测量规范还要求，航空磁测磁力仪动态噪声水平(Si)最低要求要≤0.2 nT(三级数据)，把3倍的Si作为剖面提取异常的依据。所以认为对于合格的航磁测量数据，出现>0.6 nT的异常都是可信的。由此可以将≥0.6 nT且≤2倍测量精度的航磁异常定义为微弱异常。根据近年来多个测区1∶5万高精度航磁测量的测量精度统计结果来看，全测区的航磁测量总精度一般在1.5～2.9 nT，多数在1.8 nT左右。如果除去基岩出露区数据，只使用航磁数据变化较小的盆地区数据计算航磁测量精度，其精度会更小。所以2倍测量精度最高可界定在6 nT。

通过以上分析认为，对合格的航磁数据而言，可以把≥3Si(即≥0.6 nT)且小于2倍测量总精度(约6 nT)界定为航磁微弱异常的范围，即把0.6～6 nT的异常界定为航磁微弱异常。

在重磁弱异常可信度的评价方面，刘云祥经过统计研究发现，随着弱异常规模增大和规律性增加，即使幅度在1～2.5倍异常总精度，其可信度也会大大增强[18]。当前使用的高灵敏度氦光泵磁力仪进行的高精度航磁测量，完全可以探测到这些微弱磁异常。所以利用航磁测量数据进行砂岩型铀矿成矿预测的关键是如何将高频低幅0.6～6 nT的微弱磁异常提取出来。

3 航磁微弱异常特征及提取方法

为了对砂岩型铀矿的航磁微弱异常进行分析研究，本次以最新的1∶5万高精度航磁数据为基础，内蒙古中西部已发现的砂岩型铀矿为研究对象，对提取砂岩型铀矿上方微弱磁异常的方法进行了探索实验。

3.1 内蒙古西部某铀矿航磁ΔT异常特征

内蒙古西部某铀矿床为层间氧化带型铀矿，据普查资料显示，该区域巴音戈壁组上段红色和黄色层间氧化带发育，呈不规则的带状，近EW向展布，前锋线为不规则蛇曲状。矿区自北向南为完全氧化带—氧化还原过渡带—还原带。钻探查明的矿体和矿化范围集中分布在氧化还原过渡带内[19](图2)。

图2 内蒙古西部某铀矿航磁剖面平面Fig.2 The aeromagnetic map of uranium deposit in the west of Inner Mongolia

从矿区航磁ΔT剖面平面图上可以看出，由于矿床所在盆地中—新生界盖层均为无磁或弱磁性，基底为古生界弱磁性，矿床区域位于航磁ΔT负磁异常中，航磁异常仅表现为淹没在负背景场中向零值线方向的缓状突起，航磁基础图件上异常表现不十分明显。

将航磁剖面平面图纵向比例尺放大后，矿床所在的氧化还原过渡带能够在连续数条剖面上显示有弱磁异常存在(图3)。为突出显示这些微弱磁异常信息，需要寻找一种提取这类异常信息的有效方法。

图3 内蒙古西部某铀矿航磁剖面平面(纵向比例放大后)Fig.3 The aeromagnetic map of uranium deposit in the west of Inner Mongolia (large vertical scale)

3.2 航磁弱磁异常提取方法

砂岩型铀矿产于中—新生代盆地，由于盆地盖层一般为无磁(或弱磁)，区域磁场通常表现为平缓的负背景场(或正背景场)。砂岩型铀矿引起的航磁弱异常一般为叠加在平缓磁场背景上的低幅高频异常，航磁基础图上一般没有明显的异常显示。提取低幅高频异常较为有效的一般方法是选用合适的滤波器和滤波窗口对原始磁异常数据进行高通滤波，能够滤掉低频、突出高频的信息。

3.2.1 滤波器选择

首先建立已知铀矿及其外围航磁数据库，数据范围覆盖已知的氧化带、氧化还原过渡带(见矿区域)及还原带。选取穿过矿区的L3790测线航磁数据进行滤波器选择对比实验。共选择了高通、巴特沃斯滤波、余弦滤波、高斯滤波4种滤波器对磁数据方法进行高通滤波处理。4种滤波器进行高通滤波处理后，在氧化还原过渡带均能提取出弱磁异常，利用巴特沃斯及余弦滤波器提取的弱磁异常局部异常较为突出，与氧化还原环境的分带性对应较好，本次选择使用了巴特沃斯滤波器。

3.2.2 滤波器窗口宽度选择

选定滤波器后进行了不同尺度滤波窗口对比实验，窗口尺度分别选择100、250、500、1 000、1 500、2 000、2 500、5 000、7 500 m，处理结果见图4。

图4 L3790线航磁原始数据及不同滤波窗口提取微弱磁异常剖面对比Fig.4 Aeromagnetic ΔT data and different filtering window to filtering results extracted weak magnetic anomaly profile comparison chart of L3790

经对比分析，如果滤波窗口选择太窄(100～250 m)，获得的高频信号振幅太小，多数小于0.6 nT，其中多为地面噪声信号，如输电线塔等；滤波窗口选择在1 500～2 500 m就可以提取出小于2倍测量总精度(约6 nT)的微弱异常；如果滤波窗口选择太宽(5 000 m以上)，高频信息较弱，异常宽度变大，局部异常信息减弱。为此，在滤波器的组合选择时，首先选择了巴特沃斯高通滤波器(窗口宽度为250 m)以滤除地面高频干扰异常，再选择巴特沃斯高通滤波器(窗口宽度为2 500 m)滤除埋藏较深的磁源体引起的低频异常；实际应用过程可采用航磁网格数据进行巴特沃斯带通滤波(滤波窗口250～2 500 m)。处理结果显示，在氧化还原过渡带(见矿范围)都有明显的弱磁异常(图4)；矿床上方的航磁弱磁(微磁)异常在1～7 nT，异常靠近还原带一侧；完全氧化带弱磁异常信息则小于0.6 nT，多为噪声信息。

3.3 铀矿床航磁弱磁(微磁)信息提取结果

从图2看出，矿区范围内的航磁ΔT异常很不明显，基本上为平缓的负磁场。但航磁ΔT数据经过带通滤波处理后，突出了磁数据的低幅高频信息，处理结果成图后(图5a、图5b)可以看出氧化还原过渡带(见矿范围)基本上位于弱磁异常的正负变化梯度带位置，靠近还原带微弱磁异常明显，异常幅值在1～7 nT，而且在连续相邻测线上都有异常反映，异常分带特征明显。由此看出本次研究选用的提取航磁微弱磁异常处理方法是较为合理有效的。

图5 内蒙古西部某铀矿床航磁微弱异常信息Fig.5 The aeromagnetic weak information maps of a uranium deposit in lnner Mongolia

3.4 验证效果

内蒙古二连盆地某铀成矿带已发现多个砂岩型铀矿为古河谷层间氧化带砂岩型铀矿床，是该铀成矿带重点找矿类型之一[20]。为了验证砂岩型铀矿上方是否同样存在航磁微弱异常，并确定航磁微弱异常提取方法的有效性，在该区域开展了航磁微弱异常提取实验。

内蒙古二连盆地大部分区域已于2019年完成1∶5万高精度航磁测量，航磁数据覆盖3个已发现的古河道砂岩型铀矿区，分布在工作区的东部(A矿)、中部(B矿)和南部(C矿)[20-21],这为在该区域进行砂岩型铀矿的航磁微弱异常研究创造了良好的前提条件。

首先对3个已知铀矿的航磁ΔT异常特征进行分析。在3个铀矿上方的航磁ΔT剖面平面图上均未见明显的航磁异常，利用巴特沃斯带通滤波(250～2 500 m)方法处理后，在3个铀矿床上方均提取出明显的航磁微弱异常。

3.4.1 东部(A铀矿区)

二连盆地东部A铀矿床控矿成因为辫状河—层间氧化型，控矿沉积体系为辫状河沉积体系，矿体受水退体系域内辫状河控矿沉积体系控制。地球化学结构面主要呈NE-SW分布，呈带状展布，矿体明显受氧化带前锋线控制[22]。

据刘波等在该铀矿所做的铀矿床环境参数中与铁磁性矿物有关的的参量(TFe/%)统计结果显示，不同的地球化学分带其TFe含量有明显不同，其中黄色氧化带为1.42%，黄绿色氧化带为2.20%，氧化还原过渡带为2.98%，还原带为3.03%[22]，说明在氧化还原过渡带(铀成矿有利地段)和还原带其TFe含量明显高于氧化带。这些铁磁性矿物含量的变化可以在铀矿床上方引起航磁微弱的磁场变化。另外，核工业208队在A矿床上方的地磁测量剖面上发现存在局部跳跃变化的现象(图6)[23]。而过渡带则位于前锋线附近，特别是偏向蚀源区方向。总体上看，氧化还原过渡带引起的磁异常变化值不大，一般在10～30 nT。

图6 A砂岩型铀矿地磁测量剖面[25]Fig.6 Geomagnetic survey map of of A uranium deposit[25]

现在需要解决的问题是如何在高精度航磁数据中提取这些微弱磁异常。其中A矿区域为20 nT平缓的正磁背景，未见明显的局部异常，推测为深部磁性基底的磁场特征，其南侧见由隐伏岩体引起的宽缓的磁异常(图7a)。航磁ΔT数据进行微弱磁异常提取处理后，在A矿区范围内提取出1～8.0 nT左右的微弱磁异常，多条测线均见微弱异常反映，大致呈NE向带状展布(图7b、图7c)。西南段的两个已知矿体位于NE向展布的0.5～6 nT微弱磁异常南侧由正负变化较陡的梯度带上，矿体走向与异常展布方向一致。中段的已知矿位于片状微弱异常的西北部。东北部主矿体位于NE向展布的0.5～8 nT微弱磁异常的北侧，矿体走向与异常展布方向一致。部分矿化孔及异常孔与微弱磁异常对应。

图7 A铀矿区航磁ΔT、微弱磁异常Fig.7 The aeromagnetic ΔT and weak magnetic information maps of uranium deposit area A

3.4.2 中部(B铀矿区)

二连盆地中部B矿区，其东部矿体位于宽缓的磁异常边部，西部矿体整体位于负磁背景场中，除矿床的西北和东北侧见隐伏岩体引起的宽缓磁异常外，矿体上方未见明显的航磁ΔT异常(图8a)。

图8 B铀矿区航磁ΔT、微弱磁异常Fig.8 The aeromagnetic ΔT and weak magnetic information maps of uranium deposit area B

航磁ΔT数据进行提取微弱磁异常处理后，在矿区范围提取出1～10 nT左右的微弱磁异常，多条测线均见微弱磁异常反映，大致呈NE向带状展布(图8b、图8c)。西侧主矿体西段对应微弱磁异常，东段对应微弱异常的正负变化梯度带，矿体走向与微弱磁异常展布方向基本一致。矿区内其他见矿孔及矿化孔多数与微弱磁异常对应。

从A、B两个铀矿区提取的航磁微弱异常展布方向与矿体走向的一致性及见矿及矿化钻孔的对应性分析可以推断，这两个区域的铀成矿与航磁微弱异常关系十分密切。分析这些微弱异常的成因可能有两点，一是由矿区的还原带—氧化还原过渡带—氧化带环境的磁性变化引起的异常，另一个原因是由与铀成矿关系密切的局部盖层断裂构造活动引起的磁异常。断裂构造是砂岩型铀矿成矿作用中的重要成矿条件，有利于铀的迁移和富集，它既是还原性物质(炭屑)的储存地又是还原性物质(H2S、烃类)的运输通道。因此，在有利铀成矿地区，断裂构造周围地层往往会形成一个氧化还原过渡带[24]。这些航磁微弱异常的变化对寻找铀成矿有利的成矿环境具有重要的指示意义。

3.4.3 南部(C铀矿区)

二连盆地南部C铀矿床区域，航磁ΔT以平缓的正背景场为主，矿体上方仅见特别低缓的航磁ΔT异常，异常不十分明显(图9a)。航磁ΔT数据进行微弱磁异常提取处理后，在矿区附近提取出1～5 nT左右的微弱磁异常，多条测线均见异常反映，大致呈NE向带状展布(图9b、图9c)。但现查明的已知矿体的走向与提取的微弱异常展布方向有些不一致，矿体走向以近SN向为主，而异常展布方向以NE向为主。矿体的位置与微弱磁异常的关系有待进一步研究。

图9 C铀矿区航磁ΔT、微弱磁异常Fig.9 The aeromagnetic ΔT and weak magnetic information maps of uranium deposit area C

4 结论

通过以上对内蒙古中—新生代盆地内3个砂岩型铀航磁ΔT及航磁微弱异常特征分析，可以看出，中—新生代盆地的砂岩型铀矿一般位于航磁ΔT平静的负磁背景场或平静的正磁背景场中，航磁异常不十分明显。

利用本文研究的方法进行航磁微弱异常提取后发现已知砂岩型铀矿矿体范围、工业孔及矿化孔上方或附近确实存在1～10 nT左右航磁微弱异常存在。这些航磁微弱异常对寻找成矿有利的地球化学环境和构造环境具有很好的指示意义。

1) 通过资料及数据分析对航磁微弱磁异常的量值进行了界定，把0.6～6 nT的航磁异常界定为微弱磁异常。

2) 分析总结了砂岩型铀矿弱磁(微磁)异常的形成机理，在内蒙古西部某铀矿床上方开展航磁微弱异常提取方法实验研究，总结了一套提取航磁微弱异常的方法，提取出明显的弱磁(微磁)异常。

3) 利用研究总结的航磁微弱异常提取方法，采用1∶5万高精度航磁数据，在二连盆地某铀成矿带3个已发现的砂岩型铀矿区域进行了提取航磁微弱磁异常验证实验，结果表明，在A、B、C3个古河道砂岩型铀矿区上方均存在航磁微弱异常，其中A、B铀矿区的航磁微弱异常展布方向与矿体走向一致性较好，C铀矿区较差。本文使用的航磁微弱异常的提取方法具有较好的使用效果。

4) 利用提取的航磁微弱异常可以寻找砂岩型铀矿成矿有利地化环境和构造环境，为铀成矿预测提供更丰富的信息。砂岩型铀矿航磁微弱异常提取方法的研究，为在中—新生代盆地区利用高精度航磁测量数据资料寻找古河道型砂岩型铀矿探索出一条新途径。