Property analysis and application of multi-scale wavelet decomposition of gravity potential field
MENG Qing-Kui1,2(), ZHANG Wen-Zhi2, GAO Wei2, SHU Qing1,2, LI Rui2, XU Guang-Jing2, ZHANG Kai-Song2
1. Key Laboratory of Airborne Geophysics and Remote Sensing Geology, Ministry of Natural Resources, Beijing 100083, China 2. China Aero Geophysical Survey and Remote Sensing Center for Natural Resources, Beijing 100083, China
Multi-scale wavelet decomposition is one of the common methods for gravity potential field separation. The biggest advantage is that it breaks through the traditional concept of dichotomy gravity anomalies and achieves the multiple decomposition of gravity anomalies. However, no systematic study has been carried out on the limitations and properties of the multi-scale wavelet decomposition. To systematically investigate the multi-scale wavelet decomposition and provide guidance for its practical application, this study, starting from the definition of multi-scale wavelet decomposition of gravity potential field based on profiles and grids, expounded three important properties such as the criterion that low-order wavelet keeps details invariant. Then, it analyzed the definition and properties of multi-scale wavelet by designing simple and complex theoretical models. Afterward, this study compared the multi-scale wavelet decomposition with the interpolation cutting method using field data. The results show that multi-scale wavelet decomposition can achieve multi-layer separation of gravity potential field and estimate the burial depths of source bodies. In addition, the multi-scale wavelet decomposition can provide some ideas for solving the limitations pointed out in this study, such as abnormal scale aliasing and difficulty with the determination of scale coefficients. The above basic research can provide a certain degree of references for the processing and interpretation of gravitational potential field data.
Nettleton L L. Regionals,residuals,and structures[J]. Geophysics, 1954, 19(1):1-22.
doi: 10.1190/1.1437966
[2]
Agocs W B. Least-squares residual anomaly determiration[J]. Geoghysics, 1951, 16(4):686-696.
[3]
Oldham C H G, Sutherland D B. Orthogonal polynomials: Their use in estimating the regional effect[J]. Geophysics, 1955, 20(2):295-306.
doi: 10.1190/1.1438143
Liu J Z, Zhu L T, Zhang S X, et al. Gravity anomalies vertical separation at different depths by an improved 2D polynomial fitting and weighting iterative algorithm[J]. Journal of Geodesy and Geodynamics, 2018, 38(9): 897-902,907.
[5]
Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis[J]. Proceedings of the Royal Society, 1998, 454(1971):903-995.
doi: 10.1098/rspa.1998.0193
[6]
Mickus K L, Aiken C L V, Kennedy W D. Regional-residual gravity anomaly separation using the minimum-curvature technique[J]. Geophysics, 1991, 56(2):279-283.
doi: 10.1190/1.1443041
Ji X L, Wang W Y, Qiu Z Y. The research to the minimum curvature technique for potential field data separation[J]. Chinese Journal of Geophysics, 2015, 58(3):1042-1058.
Ji X L, Wang W Y, Qiu Z Y. The parameter choose experimental research to the minimum curvature technique potential field data separation method[J]. Progress in Geophysics, 2019, 34(4):1441-1452.
Cheng F D, Liu D J, Yao R X. A study on the identification of regional and local gravity fields[J]. Computing Techniques for Geophysical and Geochemical Exploration, 1987, 9(1):1-9.
Wen B H, Cheng F D. A new interpolating cut method for identifying regional and local fields of magnetic anomaly[J]. J.Cent-South Inst.Min.Metall, 1990, 21(3):229-235.
Liu D J, Chen F D. A multiple-cut method for identification of regional and local gravity fields[J]. Computing Techniques for Geophysical and Geochemical Exploration, 1997, 19(1):32-36.
Xu S Z, Zhang Y, Wen B H, et al. Application of cut method in magnetic anomaly interpretation in LuDong area[J]. Geophysical Prospecting for Petroleum, 2006, 45(3):316-318.
Ge C, Ren S L, Li Y D, et al. Improvement and application of the layered separation method for gravity anomalies: An example of the wuhe area,Anhui Province[J]. Chinese Journal of Geophysics, 2017, 60(12):4826-4839.
Yu C C, Xiong S Q, Guo Z H, et al. The improved nonlinear filtering mthod and its application in middle and high mountain areas[J]. Geophysical and Geochemical Exploration, 2003, 27(1):39-42,62.
Guo Z H, Liu H J, Xiong S Q. The nonlinear curvature filtering technique of plane potential grid data in space domain[J]. Progress in Geophysics, 2003, 18(1):134-137.
[16]
Kaftan I, Salk M, Suri C. Application of the finite element method to gravity data case study: Western Turkey[J]. Journal of Geodynamics, 2005, 39(5):431-443.
doi: 10.1016/j.jog.2005.04.003
[17]
Kannan S, Mallick K. Accurate regional-residual separation by finite element approach, Bouguer gravity of Precambrian mineral prospect in northwestern Ontario[J]. First Break, 2003, 21:39-42.
[18]
Mallick K, Sharma K K. A finite element method for computation of the regional gravity anomaly for computation of the regional gravity anomaly[J]. Geophysics, 1999, 64(2):461-469.
doi: 10.1190/1.1444551
[19]
Mallick K, Sarma G S. Finite element formulation for seismic wave propagation in oil-bearing geological formations[C]// International Symposium on Mathematical Modeling and Scientific Computation, 1992.
[20]
Sarma G S, Gadhinglajkar V R, Mallick K. Finite element simulation of bright-spot structures[J]. Journal of Association Exploration Geophysics, 1993, 14(2):43-47.
[21]
Agarwal B N P, Shalivahan S. A fortran program to implement the method of finite elements to compute regional and residual anomalies from gravity data[J]. Computers & Geosciences, 2010, 36(7):848-852.
doi: 10.1016/j.cageo.2009.11.002
Xiao F, Meng L S. Use finite element method for computation of the regional gravity anomaly[J]. Journal of Jilin University:Earth Science Edition, 2006, 36(S):5-8.
[23]
Yu Y Y, Huang G X. Application of kalman filtering theory to regional-residual separation of gravity or magnetic anomalies[J]. Computing Techniques for Geophysical and Geochemical Exploration, 1991, 13(3):220-228.
Xu J S, Meng L S, Liu Y P. Application of kalman filtering in the separation of potential field[J]. Science Technology and Engineering, 2010, 10(2):394-399.
[25]
Pawlowski R S, Hansen R O. Gravity anomaly separation by Wiener filtering[J]. Geophysics, 1990, 55 (5):539-548.
doi: 10.1190/1.1442865
[26]
Spector A, Grant F S. Statistical model for interpreting aeromagnetic data[J]. Geophysics, 1970, 35(2):293-302.
doi: 10.1190/1.1440092
[27]
Syberg F J R. A Fourier method for the regional residual problem of potential fields[J]. Geophysical Prospecting, 1972, 20(1):47-75.
doi: 10.1111/j.1365-2478.1972.tb00619.x
[28]
Charaborty L B. Mapping of crustal discontinuities by wavelength filtering of gravity field[J]. Geophysical Prospecting, 1992, 40(7):801-822.
doi: 10.1111/j.1365-2478.1992.tb00553.x
Wang M, Wang L F, He H. The application of the matched filtering technology to the separation of gravity field sources[J]. Geophysical and Geochemical Exploration, 2015, 39(S1) : 126-132.
Yang W C, Shi Z Q, Hou Z Z, et al. Discrete wavelet transform for multiple decomposition of gravity anomalies[J]. Chinese Journal of Geophysics, 2001, 44(4):534-541.
Gao D Z, Hou Z Z, Tang J. Multiscale analysis of gravity anomalies one east China sea and adjacent regions[J]. Chinese Journal of Geophysics, 2000, 43(6):842-849.
Yang W C, Sun Y Y, Hou Z Z, et al. An multi-scale scratch analysis method for quantitative interpretation of regional gravity fields[J]. Chinese Journal of Geophysics, 2015, 58(2) : 520-531.
doi: 10.1002/cjg2.20192
Yang C B, Liu J Y, Wu Y G, et al. Two-dimensional multi-scale discrete wavelet/wavelet packet analysis for models of regional and residual gravity anomaly separation[J]. Progress in Geophysics, 2010, 25(3):1007-1014.
[36]
Banerjee B, Das Gupta S P. Gravitational attraction of a rectangular parallelepiped. Geophysics, 1977, 42(5):1053-1055.
doi: 10.1190/1.1440766
Wang S R, Zhang H D, Chen J. Application of interpolation cut method to gravity data processing in South China Sea[J]. Marine Geology & Quaternary Geology, 2017, 39(1):173-183.