小波分析和其在地球物理学中的应用.pdfVIP

小波分析及其在地球物理学中的应用 一、引言 曾经听过刘光鼎先生说过，地球物理学是地球科学中高科技，上天入地无所不能。当时 我就在想，这高科技到底体现在哪呢？难道只能说是其本身专业的优势吗，应该不止如此吧。 之所以能称之为高科技，就是相当多的高深的数学与物理的知识，还有计算机的引入。记得 近代物理课的老师曾经说过，我们学物理课的最初动机，就是能把最先进的物理知识应用到 专业课的学习中，数学的学习也应该是这样的吧。 地球物理方法是研究地壳深部构造有力且可靠的工具, 其中重磁方法所取得的位场资 料是研究地质构造、岩石圈结构和矿产勘查的重要基础性资料川。多年来, 位场资料处理的 主要工具仍是 Fourier 分析, 而面对非平稳的位场资料时,Fourier 分析已突显其局限性。近年 来引人位场的小波分析是源于 Fourier 分析而又胜于 Fourier 分析, 有“ 数学显微镜” 之称。 小波分析因其具有自适应性和多尺度分析特性成为位场数据处理的有力工具。位场领域的小 波方法研究经历十余年, 且取得了诸多研究成果,主要应用于位场多尺度分解、重磁异常奇异 性检测、位场反演等方面[1]。 本文先介绍最近 20 年才发展起来的数学分支——小波分析，接着介绍小波分析在地球 物理学中的应用，把最先进的数学方法引入地球物理学中，这才是地球物理学的高科技之处。 二、小波简介 1、小波发展历史 众所周知, 傅里叶分析从 19 世纪初、中期到今, 在数理科学一直享有盛誉, 它为信号 分析、量子力学、声学、电学、数理方程、图像处理、地球物理等许多学科领域提供了科学 的分析方法和强有力的工具。然而, 在应用中人们不断地发现了它的不足之处, 它不能刻划 函数所在的空间L, 即e 并不是L 的无条件基，它只能获得 f 的整体频谱, 而不能获得信 号的局部特性, 即它不能用于局部分[1]。近几年来, 一种被称为小波变换的数学理论和方法 正在科学技术界引起一场轩然大波，小波变换综合了三角函数系与 Haar 系两者的优点, 用 小波基来分解任意函数, 它具有优良的 “变焦”性能; 小波分析是一个新的现代分析学如泛 函分析、数值分析、Fourier 分析、样条分析等的完善结晶, 被誉称为 “数学显微镜”, 它是 一种窗口大小不变但形状可变的时频局部化分析方法。小波变换克服了上述 Fourier 变换存 在的问题，能将各种交织在一起的不同频率的混合信号分解成不同频率的块信号, 因而能有 效地用于如信噪分离、编码解码、奇异性检测、压缩数据、识别模式以及将非线性问题线性 化、非平稳过程平稳化等问[2]，这已成为应用数学的新趋势。 小波变换的思想来源于伸缩与平移方法。小波分析方法的提出, 最早应属 1910 年 Haar 提出的规范正交基, 但当时并没有出现 “ 小波”这个词。1936 年 Littlewood 和 Paley 对傅立 叶级数建立了二进制频率分量分组理论, 对频率按二进制进行划分, 其傅立叶变换的相位变 化并不影响函数的大小, 这是多尺度分析思想的最早来源。1946 年 Gabor 提出的加窗傅立叶 变换或称短时傅立叶变换对弥补傅立叶变换的不足起到了一定的作用。后来,Calderon 、 Zgymund 、Stem 等将 L 一 P 理论推广到高维, 并建立了奇异积分算子理论。1965 年 Calderon 发现了再生核公式, 它的离散形式已接近小波展开, 只是还无法得到一个正交系的结论。 1981 年, Stormberg 对 haar 系进行了改进, 证明了小波函数的存在性。1982 年 Battle 在构造 量子场论中采用了再生核公式的展开形式。 1984 年, 法国地球物理学家 J.