重磁异常分离Nuady非线性滤波的研究毕业论文

重磁异常分离：Nuady非线性滤波的研究

Gravity and magnetic anomaly separation : study of the

Nuady nonlinear filtering

目录

绪论 (1)

第一章研究背景 (3)

第二章方法原理 (5)

第一节Nuady非线性滤波法 (5)

第二节插值切割法 (6)

第三章理论模型 (8)

第一节参数设置 (8)

第二节matlab编程成图解释分析对比 (9)

第三节理论模型的误差分析表格 (20)

第四节理论模型结论 (23)

第四章实际资料处理 (23)

第一节资料获取工区概况 (23)

第二节 matlab编程成图解释分析对比 (24)

结论和建议 (31)

参考文献 (32)

致 (34)

绪论

近些年，随着我国经济的飞速发展，矿产资源的短缺的问题已成为限制国民经济持续快速发展的重大问题。所以，资源勘探成了重中之重，地球物理勘探方法有了用武之地。资源勘探中，重磁勘探是最重要勘探方法之一，重磁测量及其数据解释方法具有广阔的前景。尤其是在寻找磁铁矿及其共生磁性矿藏的工作里，重磁勘探有着不可代替的作用。重磁勘探具有仪器轻便、采集数据快、经济等特点，可以在地面、海洋、空中以及卫星上获取有用的观测数据。对于沙漠、森林等环境非常恶劣的地域，资源勘探仍然存在很多盲区，重磁航测的发展对这些盲区有非常重要的意义。这些都充分彰显重磁方法在地球物理勘探里的先导和直接寻找资源的主导作用。

重磁勘探中，数据资料处理是最重要的环节之一，同时也是重磁资料解释的前提。在生产过程中证实:实际的勘探工作里采集到的重磁资料，通常是由于各种不同的因素产生的叠加异常，其包含了许多有用的宝贵信息，是不同的规模和埋藏深度的地质体的综合反馈。若是能从采集到的重磁数据中将不同地质体产生的各种异常分离出来，就可以分别研究它的地质意义。该怎样将重磁数据中的有用的异常很好的提取出来是现在重磁勘探中需要解决的主要问题，所以一直受到国外的众多学者的重视。一般我们可以将叠加异常看成为区域异常与局部异常的叠加组成，将从重磁数据中分离出的有用异常的问题，转变为区域异常与局部异常的分离。这种分离是多种处理方法的最基本的，准确的分离区域异常与局部异常，可以为重磁数据的反演解释打下坚实基础，可以让资料解释有较好的地质效果。

在最开始的时候, 人们把区域场视为恒定的场, 采取徒手圆滑的技术, 再后来渐变为数字处理的方法, 例如移动平均法等。50 年代开始有了数字滤波概念而且用于异常分离。再后来后来还有人对异常分离滤波的误差性质做了比较深入探讨, 还提出提高最佳线性滤波的技术等, 为了克服摺积滤波对波带较宽位场波谱的不适应性, 同时也提出用非线性滤波方法对叠加场进行的异常分离, 这种方法仅仅限于航磁剖面数据的处理。

区域场概念方面的多项式拟合曲面,最初做了很多探索正交多项式类型但效果并不好,如重力区域领域多元多项式系数多项式拟合的地区区域异常效果,导致获得的数学的表面可能不是一个真正的区域场近似。因此,我们建议根据一些数据的一部分没有区域场效应,使用迭代方法的逐次逼近的方法或试图消除区域异常效应的数据,然后寻求表面多项式系数,这样多项式表面提出法律障碍要克服上述需要适当的同时,大大增强了活力的方法。60年代初,美国地球物理气的哥哥亚历克斯,一个多项式拟合曲面区域的领域,后来开发的更先进的表面来取代并接近。50年代后期,我们建议常规多项式描述区域场,为了提高计算速度和准确度多项式拟合。在1970年代联使用双正交函数转换公式来转换多项式近似的引力场,从而达到目的的分离区域场异常。孤立的异常通常有很多的歪曲,很难定量解释。一些人认为的标题字段将导致混乱的地区在感觉仿佛被重新命名为“领域的趋势”是适当的。此外,联还可以根据地质、地球物理场源的先验信息给模型(或接口深度参数)是该算法分离的例外不同类型的字段来源(或相关性分析来分离不同的接口),它将会分离不同寻常的主题和建立的参数区域背景场源。到目前为止,外国地球物理学家相信势场异常分离理论虽然它是一个反问题,但是在确定区域领域和地区分布有着本质的联系,所以,作为解决方案的反问题,分离的区域场也有多个解决严格数学的观点,不可能获得的独特的解决方案。

naudy非线性滤波法与插值切割法在异常分离中产生的效果是非常好的。此次论文就是以此两种方法作为一次对比，对两种方法产生的效果进行透彻的分析。

第一章研究背景

重磁异常分离方法，国外多年以来进行了许多的研究与开发，发表了许多不同有效的新观点技术，近来还有人在发表这方面的研究成果。分离区域场与局部场的方法技术大体上分两类:一类是空间域滤波;另一类是频率域滤波。在分离空间域中有徒手圆滑法、最小二乘法、插值法等;在分离频率域中有维纳滤波法、小波逼近法、最佳线性滤波法等。有关异常分离方法的分类如下图：

非线性滤波的使用是Naudy（1968）首次推荐使用的，非线性滤波在磁航剖面图去噪音过程中取得了很好的效果。从那时起，非线性滤波用来去除各种地球物理数据异常，基本的原则是异常完全被移除并且这些剩余异常应该保持不被改变。非线性滤波在地球物理的应用中不断的被完善，例如在2001年，Pierre Keating ,和 Nicolas Pinet发表的非线性滤波在分离局部场和区域场中的应用。并且在2003年，由、熊盛青、郭志红等人发表的改进的非线性滤波在高山地区的应用都能能说明非线性滤波在地球物理异常分离中取得了较好的效果。

非线性滤波中的插值切割法是文百红和程方道在1990年提出的一种新的异常分离方法技术，该方法以当前计算点场值与四点圆周平均值的插值运算为切割算子, 通过连续切割, 得到磁异常的切割区域场, 而从磁异常中减去后, 就得

到切割局部场。在1991年的2月，文百红更加详细的介绍了插值切割在消除磁异常在干扰中的应用，在经过前人的不断研究完善，插值切割法已经对异常分离处理有了非常明显的效果。