第二代小波变换及在不规则测点重磁资料处理中的应用.pdfVIP

第二代小波变换及在不规则测点重磁 资料处理中的应用1 刘天佑，史辉，吴小羊 中国地质大学，湖北武汉（430074 ） E-mail：liuty4508@126.com 摘 要：1994年swelden提出了基于提升算法的第二代小波，它继承了第一代小波的多分辨 特性；不依赖傅立叶变换，小波变换后的系数是整数，运算速度快；并且可以实现不规则测 网数据的小波分析。本文实现了基于提升算法的第二代小波变换，并把它应用于不规则测点 的重力数据的处理，该方法比预先将不规则测点的重力数据经过二次插值网格化，再进行第 一代小波分析的方法不仅精度高，而且失真小。它可用于1：5万～1：20万石油高精度重磁 勘探中对不规则测网数据的处理。最后利用第二代小波变换处理了江南古陆CHAMP卫星磁 测数据。 关键词：第二代小波 提升算法 高精度重磁勘探 不规则测网 江南古陆 CHAMP 卫星磁测 中图分类号：P3 1．引言 重磁勘探是方法理论成熟，覆盖面积广，应用领域十分广泛的两种地球物理方法。在1： 20 万或更小的比例尺重磁勘探的数据采集中，通常采用不规则测网。近年，随着人们绿色 与环保理念的增强，为了在施工中不砍伐树木、破坏生态环境，在 1：10 万，1：5 万比例 尺的重磁数据采集中也常常采用不规则测网。在石油重磁勘探中，由于被探测的目标埋深大 （通常密度界面、磁性界面深度在3～10km），产生的重磁异常弱，为了探测深部构造，近 年国内已开始采用“高精度三维重磁采集方法”，其做法是沿测点号观测一次，再沿测线号观 测一次，通过多次观测来提高观测精度。例如在我国南方复杂地形的石油重磁勘探，1：5 万重力设计精度为 0.09×10-5m/s2，而实际可达到 0.065×10-5m/s2，在野外采集这一环节， 目前国内已经可以达到相当高的精度。重磁资料数据处理，如利用傅立叶变换的频率域位场 转换，小波分析等，都要求观测数据是等间距的，即规则测网数据。对于实测不规则测点数 据，通常要先做网格化处理变为网格数据，由于对不规则测点重磁数据做了网格化（如采用 距离平方反比、克吕金法等等），原本野外采集的数据其高精度将由于网格化过程而丧失。 因此，寻找一种能够保持原始重磁观测数据高精度的处理方法具有十分重要的意义，它是实 现野外采集与室内资料处理同时高精度的重要途径。本文介绍的第二代小波变换是一种能够 直接对不规则测点重磁资料进行小波分析的新方法，由此可以对不规则测点重磁资料进行去 噪，位场分离等等处理与解释。1994 年，W.Sweldens 等人针对第一代小波的局限性，提出 [1] [2] 了一种不依赖傅立叶变换的新的小波构造算法－提升算法（Lifting scheme ） ，称之为 第二代小波变换，其主要特点有：继承了第一代小波的多分辨特性；不依赖傅立叶变换；小 波变换后的系数是整数[3] ，；基于多项式内插的思想，所有运算都在空间域进行，从而摆脱 了对频域的依赖[4] 。由于无需傅立叶分析，运算速度大大加快，且逆变换也容易实现，它还 简化了小波函数的构造（将小波构造转化为选用合适的插值算法）。对于重磁数据处理，第 二代小波变换还有一个重要应用就是可以实现不规则点数据的小波分析。 2 ．第二代小波变换的基本原理 1本课题得到高等学校博士学科点专项科研基金（项目编号：20050491504 ）的资助。 -1- 第二代小波是一种基于提升算法（Lifting scheme ）的小波，它主要包括分裂（Split ）、 预测（Predict ）、更新（Update ）三个步骤。 2.1 分裂（Split，又称为惰性小波变换） 此步骤是将原始的离散数据λ 分割为两部分r 和λ 。通常是将数列进行