摘要
重磁勘探是重要的物探方法,场源识别和位场分离是重磁数据处理中两个重要的问题。经典的频率域位场分离方法能有效分离不同深度、不同规模的异常,但由于傅立叶变换没有空间信息,存在分离不彻底和不能空间定位分离的缺点。小波变换具有多尺度分析、空间--尺度(频率)域定位分析等优点,有望解决上述位场分离和场源识别问题。
本文首先在阐述小波分析基本理论的基础上,开展基于小波变换的位场场源识别和位场分离方面的研究,主要成果如下:
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