用户名: 密码: 验证码:
基于小波变换的位场场源识别与异常分离方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
重磁勘探是重要的物探方法,场源识别和位场分离是重磁数据处理中两个重要的问题。经典的频率域位场分离方法能有效分离不同深度、不同规模的异常,但由于傅立叶变换没有空间信息,存在分离不彻底和不能空间定位分离的缺点。小波变换具有多尺度分析、空间--尺度(频率)域定位分析等优点,有望解决上述位场分离和场源识别问题。
     本文首先在阐述小波分析基本理论的基础上,开展基于小波变换的位场场源识别和位场分离方面的研究,主要成果如下:
     1.在详细推导连续小波变换识别位场场源理论的基础上,设计各种类型场源模型,详细研究连续小波变换识别位场场源方法适用条件、母小波选择原则、噪声影响和尺度因子选取方法。
     2.提出选择最佳母小波的定量指标稀疏指数和带空间-尺度定位的离散小波变换(DWT)位场异常分离方法。该方法对位场数据进行小波时频分析,确定出局部异常/区域异常对应的小波系数在空间-尺度域的范围,用特定范围内的小波系数重构出局部异常/区域异常,实现带空间-尺度域定位的位场异常分离。实验结果表明,先利用稀疏指数指标优选母小波,再采用带空间-尺度域定位的DWT位场异常分离方法,能提高位场异常分离的准确性,有效地从总异常中分离出区域异常和局部异常。
     3.提出了基于格林等效层和维纳滤波的小波域优化位变滤波方法。该方法具有空间变化滤波能力,从而可解决局部频谱与全局频谱存在较大差异的位场异常分离问题。理论模型试验结果表明,在局部频谱与全局频谱差异较小的情况下,该方法的分离效果与Butterworth滤波法和优化滤波法相当;在局部频谱与全局频谱差异较大的情况下,该方法相对于Butterworth滤波法和优化滤波法具有一定优势。
     4.对理论模型合成数据进行先位场分离再场源识别和直接场源识别进行对比实验,实验结果表明:采用带空间-尺度定位的DWT的分离方法,先分离出各个孤立局部异常再分别进行场源识别,能提高局部异常和区域异常场源识别的准确性。最后,通过实测资料处理检验上述结论的正确性。
The gravity and magnetic survey are the important geophysical methods. Theseparation and source identification of potential field are two key problems of gravityand magnetic field data processing. The classical separation method in frequencydomain can separate anomaly caused by sources of different depths and scaleseffectively. But they have some disadvantages, such as incomplete separation and cannot separate potential field with spatial localization, due to the lack of spatialinformation of Fourier transform. The wavelet transform, which has the ability ofmulti-resolution analysis and spatial localization, has been expected to solve theproblems presented above.
     In this paper, after the basic theory of wavelet analysis has been introduced, thecausing sources identification and separation of potential field have been studied. Themain results are listed as follow:
     1. The applicability of the identification of causing sources of potential field withContinuous Wavelet Transform(CWT) method, the principal of the choice of motherwavelet, the inference of noise and the choice of scaling factor have been studied,using several designed typical models, after the basic theory is re-derived in detail.
     2. The quantitative indictor, Sparse Index(SI), for choosing the preferentialmother wavelet, and the potential field separation method based on Discrete WaveletTransform (DWT) with spatial-scale localization have been proposed. The range ofwavelet coefficients in the spatial-scale domain corresponding to the regional andlocal anomaly respectively is determined by means of wavelet time-frequencyanalysis to the potential field data. The regional and local anomaly can bereconstructed using the wavelet coefficients within the range determined aboverespectively, and then the potential field separation with spatial-scale localization isaccomplished. The numerical results show that the total field anomaly can beseparated into regional anomaly and several isolated local ones, and the accuracy ofthe separation result of potential field can be improved obviously by the method ofpotential field separation with spatial-scale localization with the preferable mother wavelet chosen by the indicator SI.
     3. The preferable spatial-varying filtering method in wavelet domain has beenproposed based on the theories of the scale filtering in wavelet domain, Wienerfiltering and Green equivalent layer. The proposed method can solve the potentialfield separation problem whose local spectrum is different to global one obviouslybecause the spatial-varying filter parameter vary with position. The numerical resultsshow that the proposed method is as better as Butterworth filter and preferential filtermethod in the case of local spectrum is similar to global one, while the proposedmethod is better than the Butterworth filter and preferential filter method in the caseof local spectrum is much different to global one.
     4. The comparing experiments between identification the causing source afterseparation of potential field, and identification the causing source directly has beenmade. The conclusion can be drawn that separating the potential field into isolatedlocal anomaly before identifying using separation method based on DWT withspatial-scale localization can improve the accuracy of identification result of bothlocal and regional sources which interferes each other before separation. Theconclusion is verified by the field data processing result.
引文
Abdelrahman E M, Bayoumi A I, Abdelhady Y E, et al. Gravity interpretation usingcorrelation factors between successive least-squares residual anomalies. Geophysics,1989,54:1614-1621
    Billings S D and Herrmann F J. Automatic detection of position and depth ofpotential UXO using continuous wavelet transforms. AeroSense,2003, InternationalSociety for Optics and Photonics:1012-1022
    Bracewell R N. The Fourier transform and its applications, McGraw-Hill NewYork,1965
    Carmona H., Hwang W-L. and B. Torresani, et al. Practical Time-FrequencyAnalysis, Gabor and Wavelet Transformswith an implementation in S,1998,Academic Press.
    Cohen A, Daubechies I, Feauveau J C. Biorthogonal bases of compactly supportedwavelets. Communications on pure and applied mathematics,1992,45(5):485-560
    Cooper G R J. Interpreting potential field data using continuous wavelet transformsof their horizontal derivatives. Comp. and Geosci,2006,32(7):984-992
    Dampney C. The equivalent source technique. Geophysics,1969,34(1):39-53.
    Daubechies I. Orthonormal bases of compactly supported wavelets.Communications on pure and applied mathematics,1988,41(7):909-996
    Daubechies I. Ten Lectures on Wavelets. Philadelphia, PA., SIAM.,1992
    Fedi M. DEXP: A fast method to determine the depth to the sources of potentialfields. SEG2005Annual Meeting. Houston, U.S.A., Society of ExplorationGeophysicists,2005,72: I1-I11
    Fedi M. DEXP: A fast method to determine the depth and the structural index ofpotential fields sources. Geophysics,2007,72(1): I1-I11.
    Fedi M, Quarta T. Wavelet analysis for the regional-residual and local separation ofpotential field anomalies. Geophysical Prospecting,1998,46(5):507-525.
    Fedi M, Cella F, Quarta T, et al.2D continuous wavelet transform of potential fieldsdue to extended source distributions. Applied and computational harmonic analysis,2010,28(3):320-337
    Fedi M, Quarta T, Santis A D. Inherent power law behavior of magnetic field powerspectra from a Spector and Grant ensenble." Geophysics,1997,64:1143-1150
    Flandrin P, Rilling G, Goncalves P. Empirical mode decomposition as a filter bank.Signal Processing Letters,2004,11(2):112-114
    Florio G, Fedi M, Pasteka R. On the application of Euler deconvolution to theanalytic signal. Geophysics,2006,71(6): L87-L93
    Fuller B D. Two dimensional frequency analysis and design of grid operators.Mining Geophysics,1967,42:658-708
    Gabor D. Theory of communication, Journal of the Institution of ElectricalEngineers,1946,93(26):429-441
    Griffin W R. Residual gravity in theory and practice. Geophysics,1949,14(1):39-56
    Grossmann A, Morlet J. Decomposition of Hardy functions into square integrablewavelets of constant shape. SIAM journal on mathematical analysis,1984,15(4):723-736
    Hammer S. Deep gravity interpretation by stripping. Geophysics,1963,28(3):369-378
    Hood P. Gradient measurements in aeromagnetic surveying. Geophysics,1965,30(5):891-902
    Hornby P F, Boschetti, Horowitz F G. Analysis of potential field data in the waveletdomain, Geophysical Journal International,1999,137:175-196
    Hou Z, Yang W. Multi-scale inversion of density structure from gravity anomaliesin Tarim Basin. Science China, Earth Sciences,2011,54(3):399-409
    Hsu S, Sibuet J, Shyu C. High-resolution detection of geologic boundaries frompotential-field anomalies: An enhanced analytic signal technique. Geophysics,1996,61(2):373-386
    Huang D, Gubbins D, Clark R A, et al. Combined study of Euler's homogeneityequation for gravity and magnetic field.57th EAEG Meeting,1995
    Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and theHilbert spectrum for nonlinear and non-stationary time series analysis. Proceeding ofthe Royal Society,1998,454(1971):903-995
    Jacobsen B H. A case for upward continuation as a standard separation filter forpotential-field maps. Geophysics,1987,52(8):1138-1148
    Kaiser G. Wavelet filtering in the scale domain. SPIE's1996InternationalSymposium on Optical Science, Engineering, and Instrumentation, InternationalSociety for Optics and Photonics,1996
    Keating P, Pilkington M. Euler deconvolution of the analytic signal and itsapplication to magnetic interpretation. Geophysical Prospecting,2004,52(3):165-182
    Keating P, Sailhac P. Use of the analytic signal to identify magnetic anomalies dueto kimberlite pipes. Geophysics,2004,69(1):180-190
    Li X. Understanding3D analytic signal amplitude. Geophysics,2006,71(2):L13-L16
    Li Y, Braitenberg C, Yang Y. Interpretation of gravity data by the continuouswavelet transform: The case of the Chad lineament (North-Central Africa). Journal ofApplied Geophysics,2013,90:62-70
    Mallat S. A wavelet tour of signal processing: the sparse way, Access Online viaElsevier,2008
    Mallat S. A wavelet tour of signal proccessing, Third Edition: The Sparse Way,Elsevier,2009
    Mallat S, Zhong S. Characterization of signals from multiscale edges. IEEETransaction on Pattern Recognition and Machine Intelligence,1992,14(7):710-732
    Mauri G. Multi-scale analysis of multiparameter geophysical and geochemical datafrom active volcanic systems. Ph. D., Simon Fraser University,2009
    Moreau F, Gibert D, Holschneider M, et al. Wavelet analysis of potential fields.Inverse Problems,1997,13(1):165
    Moreau F, Gibert D, Holschneider M, et al. Identification of sources of potentialfields with the continuous wavelet transform: Basic theory. Journal of GeophysicalResearch,1999,104(B3):5003-5013
    Morlet J, Arens G, Fourgeau E, et al. Wave propagation and sampling theory-Part I:Complex signal and scattering in multilayered media. Geophysics,1982,47(2):203-221
    Nabighian M N. The analytic signal of two-dimensional magnetic bodies withpolygonal cross-section: its properties and use for automated anomaly interpretation.Geophysics,1972,37(3):507-517
    Nabighian M N. Additional comments on the analytic signal of two-dimensionalmagnetic bodies with polygonal cross-section. Geophysics,1974,39(1):85-92
    Nabighian M N. Toward a three dimensional auto matic interpretation of potentialfield data via generalized Hilbert transforms-Fundamental relations. Geophysics,1984,49:780-786
    Naidu P. Spectrum of the potential field due to randomly distributed sources.Geophysics,1968,33(2):337-345.
    Nettleton L L. Regionals, residuals, and structures. Geophysics,1954,19(1):1-22
    Newland D E. Random vibrations, spectral and wavelet analysis, Longman,1993
    Oonincx P J, Hermand J P. Empirical mode decomposition of ocean acoustic datawith constraint on the frequency range. Proceedings of the7th European Conferenceon Underwater Acoustics,2004
    Oruc B, Selim H H. Interpretation of magnetic data in the Sinop area of Mid BlackSea, Turkey, using tilt derivative, Euler deconvolution, and discrete wavelet transform.Journal of Applied Geophysics,2011,74(4):194-204
    Pawlowski R S. Green's equivalent-layer concept in gravity band-pass filter design.Geophysics,1994,59(1):69-76
    Pawlowski R S. Preferential continuation for potential-field anomaly enhancement.Geophysics,1995,60(2):390-398
    Pawlowski R S, Hansen R O. Gravity anomaly separation by Wiener filtering.Geophysics,1990,55(5):539-548
    Pilkington M, Keating P. The relationship between local wavenumber and analyticsignal in magnetic interpretation. Geophysics,2005,71(1): L1-L3
    Reid A B, Allsop J M, Granser H, et al. Magnetic interpretation in three dimensionsusing Euler deconvolution. Geophysics,1990,55(1):80-91
    Ridsdill-Smith T A. Separating aeromagnetic anomalies using wavelet matchedfilters.68th Annual International Meeting, SEG, Expanded Abstracts:550-553,1998
    Roest W R, Verhoef J, Pilkington M. Magnetic interpretation using the3-D analyticsignal. Geophysics,1992,57(1):116-125
    Sailhac P, Galdeano A, Gibert D, et al. Identification of sources of potential fieldswith the continuous wavelet transform: Complex wavelets and application toaeromagnetic profiles in French Guiana. Journal of Geophysical Research,2000,105(B8):19455-19475
    Sailhac P, Gibert D, Boukerbout H. The theory of the continuous wavelet transformin the interpretation of potential fields: a review. Geophysical Prospecting,2009,57(4):517-525
    Ridsdill-Smith R S, Thurston J B, Dai T F, et al. iSPI--The improved sourceparameter imaging method. Geophysical Prospecting,1998,46(2):141-151
    Spector A, Grant F S. Statistical models for interpreting aeromagnetic data.Geophysics,1970,35(2):293-302
    Srivastava S, Agarwal B N P. Inversion of the amplitude of the two-dimensionalanalytic signal of the magnetic anomaly by the particle swarm optimization technique.Geophysical Journal International,2010,182:652-662
    Strang G, Nguyen T. Wavelets and Filter Banks. Wellesley, MA,Wellesley-Cambridge Press,1995
    Syberg F J R. A FOURIER METHOD FOR THE REGIONAL-RESIDUALPROBLEM OF POTENTIAL FIELDS. Geophysical Prospecting,1972,20(1):47-75
    Thompson D T. EULDPH: A new technique for making computer-assisted depthestimates from magnetic data. Geophysics,1982,47(1):31-37
    Thurston J B, Ridsdill-Smith R S. Automatic conversion of magnetic data to depth,dip, and susceptibility contrast using the SPI (TM) method. Geophysics,1997,62(3):807-813
    Thurston J B, Ridsdill-Smith R S, Guillon J C. A multimodel method for depthestimation from magnetic data. Geophysics,2002,67(2):555-561
    Tsivouraki-Papafotiou B, Tsokas G N, Tsourlos P I. Regional-residual magneticfield data separation in wavelet domain. Journal of Balkan Geophysical Society,2003,6:209-220
    Tsivouraki-Papafotiou V, Tsokas G N, Hansen R O, et al. Trend removal anddetection of overlapping magnetic field anomalies by wavelet analysis. Proceedingsof SAGEEP2007conference,2007, April:1-5
    Ucan O N, Albora A M, Hisarli Z M. Comments on the gravity and magneticanomalies of Saros Bay using wavelet approach. Marine Geophysical Research,2001,22(4):251-264
    Ucan O N, Seker S, Albora A M, et al. Separation of magnetic fields in geophysicalstudies using a2-D multi-resolution Wavelet analysis approach. Journal of BalkanGeophysical Society,2000,3:53-58
    Vallee M A, Keating P, Ridsdill-Smith R S, et al. Estimating depth and model typeusing the continuous wavelet transform of magnetic data. Geophysics,2004,69(1):191-199
    Vallee M A, Ridsdill-Smith R S, St-Hilaire C. Magnetic depth estimates using thecontinuous wavelet transform. SEG Int'l Exposition and72th Annual Meeting. SaltLake City,Utah,2002
    Xu Y, Hao T, Li Z, et al. Regional gravity anomaly separation using wavelettransform and spectrum analysis. Journal of Geophysics and engineering,2009,6(3):279
    Zeng H, Xu D, Tan H. A model study for estimating optimum upwardcontinuationheight for gravity separation with application to a Bouguer gravity anomaly over amineral deposit, Jilin province, northeast china. Geophysics,2007,72:145-150
    Zhou Z, Adeli H. Time-Frequency Signal Analysis of Earthquake Records UsingMexican Hat Wavelets. Computer-Aided Civil and Infrastructure Engineering,2003,18:379-389
    Zurflueh E G. Applications of two-dimensional linear wavelength filtering.Geophysics,1967,32(6):1015-1035
    曾华霖.重力场与重力勘探.北京,地质出版社,2005
    陈玉东.利用位场连续复小波变换识辨磁场源(上).物探化探计算技术,2003,25(2):113-118
    陈玉东.利用位场连续复小波变换识辨磁场源(下).物探化探计算技术,2003,25(3):220-225
    陈玉东.二维连续复小波变换识别重力场源.物探与化探,2006,30(2):141-147
    刁博,王家林,程顺有.基于重力异常连续复小波变换的断裂分析.石油物探,2009,48(3):303-306
    高德章,侯遵泽,唐建.东海及临区重力异常多尺度分解.地球物理学报,2000,43(6):842-849
    管志宁.地磁场与磁力勘探.北京,地质出版社,2005
    郭良辉,孟小红,石磊,等.优化滤波方法及其在中国大陆布格重力异常数据处理中的应用.地球物理学报,2012,55(12):4078-4088
    侯遵泽,杨文采.中国重力异常的小波变换与多尺度分析.地球物理学报,1997,40(1):85-95
    侯遵泽,杨文采,刘家琦.中国大陆地壳密度差异多尺度反演.地球物理学报,1998,41(5):642-651
    胡广书.现代信号处理教程.北京,清华大学出版社,2005
    李宗杰,杨林,王勤聪.二维小波变换在位场数据处理中的应用实验研究.石油物探,1997,36(3):70-78.
    李宗杰,杨林,王勤聪.小波变换在位场数据处理中的应用.石油物探,1997,36(3):70-78
    马涛,王铁成,王雨.一种改进的网格数据保持梯度滤波方法.石油地球物理勘探,2007,42(2):198-201
    宁津生,王伟,汪海洪,等.应用小波变换确定琉球俯冲带的深部特征.武汉大学学报-信息科学版,2010,35(10):1135-1137
    邱宁,何展翔,昌彦君.分析研究基于小波分析与谱分析提高重力异常的分辨能力.地球物理学进展,2007,22(1):112-120
    汪海洪,王伟,李振,等.基于多尺度边缘的重力异常反演.大地测量与地球动力学,2008,28(5):110-114
    王刚,林江发,张赛珍.胜利油区航磁异常的小波分析.中国图形图像学报,1999,4(9):785-789
    王明,骆遥,罗锋,等.欧拉反褶积在重磁位场中应用与发展.物探与化探,2012,36(5):834-841
    熊光楚.位变滤波及其在重磁异常解释中的应用.物探与化探,1979,3(5):p.43-49
    熊光楚,张福荣.位变滤波器.物探与化探,1981,5(4):p.205-213
    徐亚,郝天珧,周立宏.位场小波变换研究进展.地球物理学进展,2006,21(4):7.
    杨高印.位场数据处理的一项新技术—小子域滤波法.石油地球物理勘探,1995,30(2):240-244.
    杨文采,侯遵泽,程振炎.重力场小波多尺度分析与大别苏鲁造山带岩石圈分区,大地测量与地球动力学进展,2004
    杨文采,施志群,侯遵泽,等.离散小波变换与重力异常多重分解.地球物理学报,2001,44(4):534-542

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700