浅析小波分析的重磁数据求导方法及应用.docxVIP

第38卷第6期吉林大学学报(地球科学版)V01．38 No．62008年11月JournalofJilinUniversity(EarthScienceEdition)Nov．2008基于小波分析的重磁数据求导方法及应用方东红1，曾昭发1，陈家林21．吉林大学地球探测科学与技术学院，长春1300262．吉林省区域地质矿产调查所。长春130022摘要：在双正交小波分析下，实现了重磁异常数据高精度的导数计算。该算法是在小波分析得到的系数基础上，采用Monte—Carlo方法进行噪声的估计，根据获得参数分别进行去噪处理。理论计算结果表明：算法的计算效率高。计算的结果精度高，抗干扰能力强。实际数据处理结果表明：采用基于小渡变换的方法完成导数处理，结果对地质填图中地质体边界的划分和磁性地层参数分析提供了重要的资料，也为采用基于导数的位场数据反演解释方法提供了优良的基础数据。关键词：重磁异常数据；导数处理；小波分析；Monte—Carlo方法中图分类号：P631．2文献标识码：A文章编号：1671—5888(2008)06—1049一06TheDerivativesCalculationBasedonWaveletAnalysisofG／M DataandItsApplicationFANGDong-hon91，ZENGZhao-fal，CHENJia—lin2t．CollegeofGeoExplorationScienceandTechnology·JilinUniversity·Changchun 130026，China2．SurveyofRegionalGeologyandMineralResourcesofJilinProvince，Changchun 130022，ChinaAbstract：Thewaveletanalysisiswidelyappliedingeophysicaldata processingandinterpretationfor itsmultipleresolutionanalysis，andgotvery goodresults．Thebiorthogonalwaveletanalysiswasused asoriginalbasis ofwaveletforhigh accuracy derivatescalculation，whichservesasimportantinforma— tionofanomalyanalysesandinterpretation，andalsoasthebasic data forfurthercalculation asEulerde— convolution，analyticinversionandsourceparametersinterpretation．Basedon the waveletcoefficients，theMonte-Carlomethodwasappliedtocalculatetheparametersofnoise．Thecasesstudyfromrealdata showsthatthealgorithmhas high accuracy，andits resultscan help thegeologicalmappingandstratastudy．Keywords：gravityandmagneticdata；derivativeprocessing；waveletanalysis；Monte—Carlometh—od常重要的常规数据处理方法。其结果一方面应用于0引言数据解释和异常分析、研究局部异常和分析地质构重磁数据的求导处理是重磁数据处理和解释非造展布方向；另一方面，许多重磁异常解释方法都收稿日期：2008一03—18‘．基金项目：国家自然科学基金项目作者简介：方东红(1963--)，男，安徽贵池人，博士研究生，主要从事地球物理资料的处理与解释方面的应用研究，Tel：0431E-mail：dhfang@jlu．edu．ca通讯联系人：曾昭发(1966--)，男，广西全州人，教授，主要从事应用地球物理教学与研究，Tel：0431Bmail：zengzf@jlu．edu．cn。万方数据吉林大学学报(地球科学版)第38卷基于重磁异常的导数值，例如解析信号的解释口、，(z)=GV2||W，(口，6)似．6(z)—da万db。(3)磁源参数成像L5]，欧拉反褶积【6卅等。莺磁异常导数∥“的计算精度和质昔商接关系到资料解释和反演结果式(2)和(3)构成了小波变换的信号分解和重建的变的精度[9。10]。由于重磁异常导数处理是一个高通滤换对。波处理。其阶数越高，对高频