小波多尺度分析在磁测数据处理中的应用.pdfVIP

小波多尺度分析在磁测数据处理中的应用

喻翔;罗照华;梁涛;陈聪;崔广贺;邓俊峰

【摘要】河南栾川罗村钼多金属矿床是与罗村岩体有关的钼多金属隐伏矿床,但隐

伏岩体的形态、位置、规模、边界及地下场源的异常特征等还没有明确的认识.为

了查明这些问题,在前人工作的基础上,运用小波多尺度分析方法对磁测数据进行分

解处理,同时对比了剩余异常和小波多尺度分析在资料处理中的应用效果,为今后钼

多金属矿的研究提供新的思路.

【期刊名称】《世界核地质科学》

【年(卷),期】2014(031)004

【总页数】7页(P607-613)

【关键词】小波多尺度分析;磁异常特征;钼多金属矿;隐伏岩体

【作者】喻翔;罗照华;梁涛;陈聪;崔广贺;邓俊峰

【作者单位】中国地质大学,北京100083;核工业北京地质研究院,中核集团铀资源

勘查与评价技术重点实验室,北京100029;中国地质大学,北京100083;河南省有色

金属地质勘查总院,郑州450052;核工业北京地质研究院,中核集团铀资源勘查与评

价技术重点实验室,北京100029;核工业北京地质研究院,中核集团铀资源勘查与评

价技术重点实验室,北京100029;中国地质大学,北京100083

【正文语种】中文

【中图分类】P631.2

近20年来，小波分析已经有了比较完善的理论与计算方法，并且广泛应用于信号

处理、图像处理、故障检测、CT成像、语言识别、通信和电子系统等众多学科与

相关领域，其应用在地球物理重磁勘查领域中也取得较好的效果［1］。侯遵泽、

杨文采等［2］对中国布格重力异常进行小波变换与多尺度分解，得到了中国大陆

地壳密度的不均匀分布情况；刘天佑等［3］将小波分析应用于大冶铁矿危机矿山

的深部找矿并且获得了明显的地质效果；李才明等［4］基于小波能谱分析对岩溶

区探地雷达目标体进行了识别；刘宏等［5］将小波分析运用到MT的资料处理；

高德章等［6］采用小波多尺度分析得到了沉积基底面和莫霍面在东海及邻区产生

的重力异常。笔者根据实测的磁异常特征和急待解决的地质问题，对罗村钼矿区磁

测资料做了化极、向上延拓、求取剩余异常和小波多尺度分解，对比分析了剩余异

常和小波分析在磁测资料处理中的应用效果，并通过反演计算推断深部存在的隐伏

矿体位置。

1方法原理［7］

1.1小波变换定义

设函数f（x）∈L2（R），定义连续小波变换（ContinuousWaveletTransform）

为：

式中：ψ（x）∈L2（R）—小波函数；f（x）—能量有限的信号的共轭函数；

a—尺度因子；b—位移因子。

1.2小波多尺度分析原理

基于多尺度分析的理论，Mallat提出了一个塔式分解算法，设{Vj}是一给定的多尺

度分析，ψ和φ分别是相应的小波函数和尺度函数。

假设{Vj2}是一个二维多尺度分析，其中Vj2＝VjVj，它们的尺度函数定义为：

小波函数为：

设二维函数小波多尺度分解方法原理：

式中：ΔT（x，y）—二维异常。

式（3）可简化为：ΔT（x，y）＝A4f＋D1f＋D2f＋D3f＋D4f

式（3）说明一个二维磁异常可以由一个4阶逼近A4f及一阶、二阶、三阶和四阶

4个细节（即D1f、D2f、D3f、D4f）所构成，这也就是磁异常多重分解［8］。

通过理论还证明了离散的二维小波变换产生的低阶小波细节具有尺度不变的特征。

它们不随小波变换的总阶数改变而改变，总阶数的增加仅仅是增加了高阶小波细节

的个数和改变了最后一个高阶的逼近。因此笔者可以根据地质目标来组合小波细节，

选择合适的高阶逼近，实现地质意义的分解。

2实例应用

2.1区域地质概况

河南栾川罗村钼多金属矿床地处华北地台华熊台缘坳陷黑沟—栾川断裂带与马超

营断裂带之间的卢氏—栾川陷褶束内，属于秦岭多金属成矿带东秦岭钼成矿亚带

［9］。区域地层、构造、变质带及岩浆岩总体呈NWW、NW向展布，区内前古

生界广泛发育，构造较为复杂，岩浆活动具明显的多旋回、多期性特征，变质作用

类型及其相应岩石种类繁多，成矿地质条件十分有利［10］。矿区南部岩浆岩分

布广泛，主要为加里东期龙王幢富铁钠闪长花岗岩体、辉长岩脉和燕山期罗村花岗

斑岩体（图1）。多金属矿床的形成与燕山期花岗斑岩小岩体的侵入密切相关，含

钼斑岩体赋存于龙王幢花岗岩体与太华群接触带附近并且矿化具明显的高温-中低

温水平分带现象。