优选向上延拓技术在鞍山—本溪示范区航磁应用
张建兵1, 王明2, 赵百民2, 刘前坤2
1.山西省地质调查院, 山西 太原 030006
2.中国国土资源航空物探遥感中心, 北京 100083
通讯作者:王明(1981-),男,山东人,高级工程师,硕士研究生,主要从事航空物探方法技术研究。Email:wangmingS060880@yeah.net

作者简介: 张建兵(1964-),男,山西人,高级工程师,大学学历,主要从事电磁勘查工作。Email:776587682@qq.com

摘要

优选向上延拓技术可以保持深源低频信号基本不变的同时,使浅部高频信号衰减,能够较为准确地提取区域异常,进而分离出局部异常;避免了传统向上延拓浅源高频信号与深源低频信号同时衰减。笔者利用优选向上延拓技术处理鞍山—本溪示范区航磁资料,实现了深部区域异常的增强和分离,获得了较好的效果;并进一步对西鞍山铁矿床航磁异常进行三维反演,初步推断了西鞍山铁矿体空间形态和规模,为后期资料的成果解释提供了重要的参考依据。

关键词: 优选向上延拓; 格林等效层; 对数功率谱; 鞍山—; 本溪示范区; 航磁
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2017)05-0951-07
The application of preferential upward continuation to aeromagnetic survey in Anshan-Benxi area
ZHANG Jian-Bing1, WANG Ming2, ZHAO Bai-Min2, LIU Qian-Kun2
1.Shanxi Institute of Geological Survey,Taiyuan 030006,China
2.China Aero Geophysical Survey and Remote Sensing Center for Land and Resources,Beijing 100083,China
Abstract

The advantage of preferential upward continuation method lies in the fact that, when it is used for upward continuation, the low frequent signal is not attenuated excessively and it can retain as much information on geological features as possible, so it can accurately separate regional anomaly and local anomaly. However, in conventional methods both high and low frequent signals are attenuated simultaneously for upward continuation. The preferential upward continuation breaks up the restrictions of the conventional upward continuation. So the preferential upward continuation is a method that has the capability for self-adaptive qualitative enhancing and separating the aeromagnetic anomalies. In this paper, the aeromagnetic anomaly in Anshan-Benxi area was processed using the preferential upward continuation to realize the enhancement and separation of the anomaly in the deep source field, with a good effect obtained. Furthermore, the authors studied 3D inversion of typical Xianshan iron ore deposit and inferred the distribution of underground lithologic boundaries and the spatial morphology of the Xianshan iron ore deposit. The results provide an important reference for the qualitative interpretation of the aeromagnetic anomaly.

Keyword: preferential upward continuation; Green's equivalent layer; radial average logarithmic power spectrum; Anshan-Benxi area; aeromagnetic survey
0 引言

鞍山— 本溪示范区位于华北克拉通北缘东段的辽东地区, 是我国重要的铁矿成矿区, 分布了诸多大型、超大型铁矿床, 国内众多学者作了大量的研究工作。目前, 鞍山— 本溪示范区已属于比较成熟的老矿区, 露天矿和浅部矿床基本已经开发或掌握储量, 现阶段找矿模式已从浅部矿产资源勘探开发向深部地质找矿转变, 特别是老矿区的深部、覆盖区等研究程度不高的区域; 同时, 国内外不断涌现的新的成矿理论对深部找矿突破提供了理论指导作用。

1971~1972年, 鞍钢地质勘探公司在西鞍山矿区及其外围先后开展了大量的重磁详查工作和钻探验证, 证实了西鞍山铁矿床的磁异常并非由地表矿体引起, 而是由深部矿体引起, 扩大了铁矿床的规模和储量。2013~2015年中国国土资源航空物探遥感中心在鞍山— 本溪示范区开展了1∶ 1万大比例尺硬架式直升机航磁测量, 获得了高精度的航磁数据, 为进一步攻深探盲、圈定鞍山— 本溪示范区重点找矿靶区提供了重要的数据资料。

在磁法勘探的数据处理中, 深源场和浅源场的分离是航磁数据处理的主要内容之一。根据分离出来航磁异常信息, 对研究或突出某一特定的地质问题有重要意义。利用深源场研究区域构造, 预测找矿远景靶区, 利用浅源场研究局部构造, 间接或直接找矿。目前, 重磁场分离的方法众多[1,15], 无论空间域还是频率域中, 理论研究与应用较为成熟, 均有各自的优势和局限性。对于大面积的航磁资料而言, 通常采用向上延拓方法突出深源低频信息, 消除浅部高频干扰, 进而反映深部构造、初步推断具有深部找矿潜力的异常。

在实际应用中, 传统向上延拓在压制浅源高频信息的同时, 深源低频信息也受到了压制; 针对上述存在的问题, Pawlowski根据温纳滤波和格林等效层原理, 提出了优选向上延拓算子, 用于突出或增强深源低频信息、有效压制高频噪声干扰的异常分离方法[9]。优选向上延拓算子的基本思想是在压制浅源局部信号的同时, 深源信号衰减最小。而这样的结果正是位场定性划分所要求的[9, 11, 16, 17]

1 基本原理

优选向上延拓算子在假定深源、浅源位场无关的前提下, 在压制浅源局部信号同时, 保持深源位场信号衰减最小。优选向上延拓算子并不是一种基于空间变换为目的的单纯延拓方法, 而是一种用于突出区域场的异常分离方法, 有助于保留位场空间换算有用信息含量。

优选向上延拓的算子为:

W=ekhup1-PdeepP(1-e-khup), (1)

其中, k是角频率, hup是向上延拓高度, Pdeep是深源场功率谱密度函数, P是位场功率谱密度函数; 式(1)括号外的指数项是传统的(理论)向上延拓的算子。

从优选向上延拓算子与传统延拓算子对比图(图1)可以明显看出, 优选向上延拓算子相对于传统向上延拓算子优势在于可以保证深源低频成分尽可能衰减小, 能够较为精确地提取区域异常; 随着不同上延高度的增加, 优选上延在保持深部低频信号不变同时, 浅部高频信号逐步衰减。而传统向上延拓算子在压制浅源信号的同时, 深源长波信号也被压制, 损失了有用信息。

图1 优选向上延拓算子与传统向上延拓算子对比

根据估算出位场的功率谱密度函数和对应的深源场、浅源场功率谱密度函数, 实现优选向上延拓计算。利用格林等效层概念建立位场的傅里叶功率谱模型, 通过对实测位场数据的径向平均对数功率谱曲线进行线性拟合(图2), 计算各等效源层频段范

图2 径向平均对数功率谱曲线(红色)和分频段线性拟合(蓝线)示意

围的线性拟合斜率和截距, 求出等效层深度和功率谱密度函数参数。

傅里叶功率谱模型是假设位场的各等效源薄层互不相关, 由随深度变化数目适当的m+n个等效源薄层叠加而成的, 其功率谱密度函数的理论模型可表示为:

P=Pd1+Pd2++Pdm+Ps1+Ps2++Psn+E, (2)

其中, 深源场成分用m个格林等效源薄层建立模型; 浅源场成分用n个格林等效源薄层建立模型; 误差E来源于对实测位场数据的径向平均对数功率谱的不完全拟合和实测数据中的噪声。

按照Naidu[18]等对于实测位场数据中任意深度处的格林等效源层的径向功率谱密度函数为:

P=se-2ωh, (3)

其中, s是与期望的等效源薄层强度成正比的常数, ω 是角频率。

2 航磁资料处理
2.1 研究区地质概况

鞍山— 本溪示范区地势东南高、西北低, 东南海拔一般在200~500 m, 属于辽东山地丘陵区, 西部属于辽河平原区, 海拔一般在10~20 m左右, 仅个别山头海拔在50~200 m。

示范区出露的地层主要有:太古界鞍山群, 元古界辽河群, 两者构成本区的结晶基底, 上覆青白口系、震旦系、寒武系、白垩系和第四系等。区内太古宙结晶基底鞍山群变质岩系发育, 分布广泛, 为鞍山式铁矿的赋存提供了良好的地层条件。含矿层以变质表壳岩的形式呈类包体状态赋存于太古宙变质花岗岩之中, 多与花岗岩呈构造或韧性剪切构造接触。

2.2 航磁异常特征

图3所示为鞍山— 本溪示范区高精度航磁异常, 异常值范围为-1 315.1~16 835.8 nT; 从图中可以看到丰富的磁场信息, 磁场面貌清晰, 特征明显, 展现了多个不同磁场背景与磁异常特征区域; 主体异常带走向以东西向和北东向为主, 前者反映了基底构造特征, 后者显示印支— 燕山运动对基底的构造改造的特点, 局部出现的特强磁场预示着太古界硅铁建造的存在, 是寻找铁矿的直接找矿标志。

示范区内弱磁异常带主要位于北部边缘的辽阳北弱磁异常带和东部边缘的辽阳东弱磁异常带。辽阳北弱磁异常带以负磁异常为特征, 呈北东东向带状展布, 向西、南逐步升高, 异常磁场最小值近-650 nT; 在区域磁场图上, 异常宏观上表现为一规模宏大的负磁异常带, 呈近东西向的“ J” 型展布。辽阳东弱磁异常带, 与辽阳北弱磁异常相邻, 以负磁异常为特征, 呈北北西向带状展布; 向西、南逐步升高, 异常磁场最小值达-750 nT。

图3 鞍山— 本溪示范区化极航磁异常

强磁异常带主要是中东部的辽阳强磁异常带和西南部东西鞍山强磁异常带。辽阳强磁异常带, 呈北北西向带状展布, 磁场值高达16 000 nT, 为示范区内最强的磁异常。东西鞍山强磁异常带, 呈北西向收敛的扇形展布, 与鞍山平缓升高磁异常相邻, 北部磁场较强, 磁场值高达12 200 nT, 向南逐步减弱。

示范区内的化极航磁异常图普遍存在着明显的高频噪声干扰, 如图3所示, 尤其是示范区中部及鞍山市区附近的磁异常区, 对后期航磁资料处理和解释产生了很大的影响。为了更好地分离出区域异常, 突出或增强深源低频信息, 有效压制高频噪声干扰, 尽可能的减少有效磁场信息的损失, 我们利用优选向上延拓技术处理示范区航磁数据, 并与传统向上延拓效果进行对比分析。

图4所示为鞍山— 本溪示范区化极航磁异常的径向平均对数功率谱曲线, 根据径向平均对数功率谱曲线变化特征, 自上而下划分为6个等效源层并进行线性拟合, 各等效源层频段范围分别是:0~0.012 695(1/m)(频段1)、0.012 695~0.600 59(1/m) (频段2)、0.600 59~0.119 629(1/m) (频段3)、0.119 629~0.358 887(1/m) (频段4)、0.358 887~0.410 156(1/m) (频段5)、0.410 156~0.5(1/m) (频段6)。结合已有的物化探和地质等资料, 分析认为频段5~6对应着高频噪声干扰, 频段1~4对应着深源场和浅源场信息。为了在后期资料解释中进一步研究示范区域航磁深部构造特征、区域构造单元、确定断裂带和地质构造走向等方面提供更多的有效信息, 我们对频段1~2进行优选向上延拓1.5 km, 分离出示范区的区域航磁异常, 并与传统向上延拓1.5 km进行对比分析, 计算结果见图5和图6。

图4 径向平均对数功率谱曲线(红色)和分频段线性拟合(蓝色)示意

从优选向上延拓1.5 km与传统向上延拓1.5 km的结果对比中, 可以明显看出, 优选向上延拓在压制浅源高频信号的同时, 尽可能更多的保留深部异常信息; 异常形态更加清晰突出, 不仅浅部高频噪声得到了很好的压制, 而且异常的形态和走向与原化极航磁异常保持着很多的对应关系, 保留了异常梯度带的陡然突变特点。而传统向上延拓在压制浅部高频信号的同时深部低频信号也被衰减, 异常变得比较宽大平缓, 损失了部分细节信息, 特别是示范区西南部的东西鞍山强磁异常和小岭子正磁异常被过分的压制而难以准确的反映出来。

图5 优选向上延拓1.5 km(左)和传统向上延拓1.5 km(右)

图6 优选向上延拓剩余异常(左)和传统向上延拓剩余异常(右)

从延拓后计算得到剩余异常结果进一步看出, 优选向上延拓剩余异常相比于传统向上延拓剩余异常而言, 包含少量的深部低频信息, 大部分为浅部高频干扰信息, 尤其是位于示范区中部的平缓磁异常带和位于示范区中西部的高磁异常带, 浅部高频噪声得到了很好的压制, 为后期资料解释等提供了更丰富细节信息; 而传统向上延拓剩余异常包含大量的深部低频信息。

因此, 优选向上延拓与传统向上延拓最大不同是可以保持深源低频信号不变, 同时使浅部高频信号得到衰减, 能够较为精确地提取深源(区域)异常, 进而分离出浅源(局部)异常; 避免了传统向上延拓浅源高频信号与深源低频信号同时衰减。

2.3 典型铁矿床反演

鞍山— 本溪示范区是我国重要的条带状铁矿资源分布区, 已经发现和开采了众多大型、超大型铁矿床, 如东西鞍山、齐大山、胡家庙子和关门山等铁矿。结合实际地质情况和矿产预测需要, 依据综合信息地质— 地球物理特征, 充分利用已知重要铁矿床所在地区地质、物探资料, 建立地质— 地球物理模型, 对优选向上延拓技术分离出来的典型铁矿床航磁异常进行三维反演, 获取地下磁性体分布的定量信息, 以西鞍山航磁异常为例。

西鞍山铁矿床位于鞍山市区的南部, 是海相火山沉积变质型矿床, 火山喷发为鞍山式铁矿提供了成矿的物质来源, 沉积作用使铁质矿源集中沉积富集成矿; 铁矿石中全铁(TFe)含量一般在28%~42%, 全矿床平均为34.22%, 铁矿床矿石类型主要为磁铁石英岩和假象赤铁石英岩、磁铁假象赤铁石英岩, 平均品位约为33%, 富铁矿石多呈脉状产出; 成矿构造背景为古华北板块的东北部鞍山— 本溪太古宙绿岩盆地, 区内出露的地层有中太古界鞍山群樱桃园组, 元古界震旦系地层主要出露在矿区的北部(图7), 它不整合覆盖在太古界鞍山群之上, 可分为3个组:钓鱼台组、南芬组和桥头组; 第四系广泛分布在矿区北部、西部和南部的平原地带。

图7 西鞍山矿区地质

鞍山式铁矿的成矿层位为中— 新太古界鞍山群樱桃园组, 该地层主要由绿泥千枚岩、绢云千枚岩、二云变粒岩、碳质千枚岩、绢云石英片岩及铁矿层所组成, 自上而下可分为3层, 即下部千枚岩层、铁矿层和上部千枚岩层。鞍山群地层与上覆青白口系地层为不整合接触关系, 普遍接受绿片岩相的区域变质作用。下部千枚岩层主要出露在矿区南部和西部, 其下部为新太古代二长花岗片麻岩, 该千枚岩厚度变化较大, 呈现东薄西厚的趋势; 上部千枚岩在矿区内零星出露。

研究成果表明[19], 矿体为巨厚层状的单一矿体, 受断裂错断影响, 分为东西两个矿体, 东部矿体走向近东西向, 而西部矿体呈北西向, 东西延长达 4 593 m, 向东与东鞍山矿层相连接, 向西隐没于冲积平原之下。东部矿体倾斜延深达1 000 m, 矿体厚度355 m; 西部矿体勘查深度达659 m, 矿体厚度为230 m; 整个矿层一般倾斜延深为800 m左右, 矿层的厚度由上向下不仅没有减薄反而有加厚的趋势, 可以推测矿层向下延深可达1 000 m以下, 倾斜延深将超过1 500 m。

在磁场图上, 西鞍山磁异常位于鞍山大型磁异常带的西南部, 航磁异常形态为北西向条带状异常, 异常体规模较大, 异常峰值高达4 397 nT, 走向长度为3 400 m, 宽度900 m, 与东鞍山异常连成一体。

图8 西鞍山铁矿体空间分布

利用磁法UBC三维反演技术对西鞍山进行地质— 地球物理联合反演(图8), 从反演结果可以看出, 西鞍山矿体为巨厚层状, 向下延伸较大, 主矿体地下尖灭深度为2 600 m, 富矿体地下延伸至2 400 m左右终止, 东西向最大延伸约900 m, 南北向最大宽度可达2 650 m, 说明已探明矿体只占整个西鞍山矿体很小一部分, 深部还有规模较大的、未探明铁矿体。三维反演结果在浅部与铁矿区16线勘探剖面对应较好, 且与已知矿体的边界、倾向、走向基本吻合, 为后期准确预测西鞍山式铁矿体形态、规模和潜在的资源量提供了数据支持和参考依据。

3 结论

优选向上延拓在假设格林等效源层的各层场源信息互不相关的前提下, 保持深源低频信息不变的同时, 压制浅源高频信息, 区别于传统向上延拓深源低频信息与浅源高频信息同时衰减。对于大面积航磁资料的常规处理与解释而言, 是一种突出深源区域异常, 消除浅部局部高频干扰的有效方法, 为定性研究区域地质构造提供了丰富有价值的细节信息。

笔者根据鞍山— 本溪示范区航磁异常的径向平均对数功率谱变化特征, 采用多层格林等效源层分频段精细拟合, 通过优选向上延拓, 实现了深部区域异常的增强和分离; 在此基础上, 并进一步对西鞍山铁矿床航磁异常进行三维反演, 初步推断了西鞍山铁矿体空间形态和规模, 为后期资料的定性解释及寻找铁矿靶区提供了准确而翔实的基础数据和依据。

The authors have declared that no competing interests exist.

参考文献
[1] Spector A, Grant F S. Statistical models for interpreting aeromagnetic data[J]. Geophysics, 1970, 35(2): 293-302. [本文引用:1]
[2] Archibald N, Gow P, Boschetti F. Multiscal edge analysis of potential field data[J]. Exploration Geophysics, 1999, 30(2): 38-44. [本文引用:1]
[3] Holden D J, Archibald N J, Boschetti F, et al. Infering geological structures using wavelet-based multiscale edge analysis and forward models[J]. Exploration geophysics, 2001, 31(4): 67-71. [本文引用:1]
[4] Marlet G, Sailhac P, Moreau F, et al. Characterization of geological boundaries using 1-D wavelet transform on gravity data: Theory and application to the Himalayas[J]. Geophysics, 66(4): 1116-1129. [本文引用:1]
[5] 侯重初. 补偿圆滑滤波方法[J]. 石油物探, 1981(2): 22-29. [本文引用:1]
[6] 安玉林, 管志宁. 滤波高频干扰的正则化稳定因子[J]. 物探化探计算技术, 1985, 7(1): 13-23. [本文引用:1]
[7] Pawlowski R S, Hansen R O. Gravity anomaly separation by Wiener filtering[J]. Geophysics, 1990, 55(5): 539-548. [本文引用:1]
[8] Pawlowski R S. Green s equivalent-layer concept in gravity band -pass filter design[J]. Geophysics, 1994, 55(5): 69-76. [本文引用:1]
[9] Pawlowski R S. Preferential continuation for potential-field anomaly enhancement[J]. Geophysics, 1995, 60(2): 390-398. [本文引用:2]
[10] Zeng H, Xu D, Tan H. A model study for estimating optimum upward-continuation height for gravity separation with application to a Bouguer gravity anomaly over a mineral deposit, Jilin province, northeast China[J]. Geophysics, 2008, 72(4): 145-150. [本文引用:1]
[11] Meng X H, Guo L H, Chen Z X, et al. A method for gravity anomaly separation based on preferential continuation and its application[J]. Applied Geophysics, 2009, 6(3): 217-225. [本文引用:1]
[12] 侯尊泽, 杨文采. 中国重力异常的小波变换与多尺度分析[J]. 地球物理学报, 1997, 40(1): 85-95. [本文引用:1]
[13] Fedi M, Quarta T. Wavelet analysis for the regional-residual and local separation of potential field anomalies[J]. Geophysical Prospecting, 1998, 46(5): 507-525. [本文引用:1]
[14] Keating P, Pinet N. Use of non-linear filtering for the regional-residual separation of potential field data[J]. Journal of Applied Geophysics, 2011, 73(4): 315-322. [本文引用:1]
[15] Nikitin A A, Vasov O K, Belov A P, et al. Vozmozhnosti kompleksnoy geofizicheskoy interpretatsii na baze entropiynogo filtra[J]. Izvestiya Akademii Nauk Turkmenskoy SSR. Seriya Fiziko-Tekhnicheskikh, Khimicheskikhi Geologicheskikh Nauk (in Russia), 1984, 2: 79-82. [本文引用:1]
[16] 许德树, 曾华霖. 优选延拓技术及其在中国布格重力异常图处理上的应用[J]. 现代地质, 2000, 14(2): 215-222. [本文引用:1]
[17] 郭良辉, 孟小红, 石磊, . 优化滤波方法及其在中国大陆布格重力异常数据处理中的应用[J]. 地球物理学报, 2012, 55(12): 4078-4088. [本文引用:1]
[18] Naidu P. Spectrum of the potential field due to rand omly distributed sources[J]Geophysics, 1968, 33: 337-345. [本文引用:1]
[19] 崔培龙. 鞍山—本溪地区铁建造型铁矿成矿构造环境与成矿、找矿模式研究[D]. 长春: 吉林大学, 2014. [本文引用:1]