1887
  1. Home
  2. A-Z Publications
  3. Near Surface Geophysics
  4. Volume 12, Issue 6
  5. Article
Volume 12 Number 6
  • ISSN: 1569-4445
  • E-ISSN: 1873-0604

Abstract

This paper describes an efficient wavelet‐based deconvolution method for the downward continuation of airborne potential (magnetic and gravity) field data. We formulate the downward continuation process as a linear ill‐posed deconvolution problem. To obtain a reasonable downward continued field data, it is stabilized in a wavelet domain by minimizing the L1‐norm of the coefficients subject to the constraint that is the agreement of the upward response with the original observed data up to a white Gaussian noise level. The resulting convex constrained problem is then solved very fast and efficiently via the split Bregman iterations method. The generalized cross‐validation (GCV) criterion as a plotted curve versus the iteration count with new formulation is used to determine the optimum number of Bregman iterations without prior knowledge about noise properties. A synthetic magnetic field data embedded in a Gaussian noise is constructed from isolated multi‐source anomalies, whose results show that the airborne magnetic data can be continued stably downward by the proposed automatic sparse deconvolution method in a few iterations. When compared with conventional methods such as Tikhonov and Edge‐preserving with the proposed method, similar results were obtained. To test the performance on real data, we chose to use the airborne magnetic data of the central Iranian volcanic‐sedimentary belt. The enhanced downward continued data to a depth level of 200 m beneath the ground surface indicates adequate matching of high potential zones with previous working mines and copper deposits. Obtained results from both synthetic and real data confirm that the performance of the proposed technique is as good as the classical deconvolution ones.

© European Association of Geoscientists and Engineers © 2014 European Association of Geoscientists and Engineers
Loading

Article metrics loading...

/content/journals/10.3997/1873-0604.2014027
2014-03-01
2024-04-24

  • From This Site

    /content/journals/10.3997/1873-0604.2014027
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. AbediM., GholamiA. and NorouziG.H.2013a. A stable downward continuation of airborne magnetic data: A case study for mineral prospectivity mapping in Central Iran. Computers & Geosciences52, 269–280.
     [Google Scholar]
  2. AbediM., GholamiA., NorouziG.H. and FathianpourN.2013b. Fast inversion of magnetic data using Lanczos bidiagonalization method. Journal of Applied Geophysics90, 126–137.
     [Google Scholar]
  3. AminiA., Fazel BidgoliY. and Rokn AbadiH.H.H.2010. Linear filter for upward‐ down ward continuation by using wavelet transform and its application in magnetic data processing. Journal of the Earth and Space Physics36, 25–39.
     [Google Scholar]
  4. AsterR.C., BorchersB. and ThurberC.2003. Parameter Estimation and Inverse Problems.Academic Press, New York, NY.
     [Google Scholar]
  5. AtapourH. and AftabiA.2007. The geochemistry of gossans associated with Sarcheshmeh porphyry copper deposit, Rafsanjan, Kerman, Iran: Implications for exploration and the environment. Journal of Geochemical Exploration93, 47–65.
     [Google Scholar]
  6. Bertero‐AguirreH., CherkaevE. and OristaglioM.2002. Non‐smooth gravity problem with total variation penalization functional. Geophysical Journal International149, 499–507.
     [Google Scholar]
  7. BlakelyR.J.1995. Potential Theory in Gravity and Magnetic Applications.Cambridge Univ. Press.
     [Google Scholar]
  8. BoschettiF., HornbyP. and HorowitzF.2001. Wavelet based inversion of gravity data. Exploration Geophysics32, 48–55.
     [Google Scholar]
  9. CandesE.J. and GuoF.2002. New multiscale transforms, minimum total variation synthesis: applications to edge‐preserving image reconstruction. Signal Processing82, 1519–1543.
     [Google Scholar]
  10. CharbonnierP., Blanc‐FeraudL., AubertG. and BarlaudM.1997. Deterministic edge‐reserving regularization in computed imaging. IEEE Trans. Image Processing6, 298–310.
     [Google Scholar]
  11. ClarkD.A.1999. Magnetic petrology of igneous intrusions – Implications for exploration and magnetic interpretation. Exploration Geophysics20, 5–26.
     [Google Scholar]
  12. CooperG.2004. The stable downward continuation of potential field data. Exploration Geophysics35, 260–265.
     [Google Scholar]
  13. CooperG.2006. Interpreting potential field data using continuous wavelet transforms of their horizontal derivatives. Computers & Geosciences32, 984–992.
     [Google Scholar]
  14. DaneshfarB.1997. An evaluation of indicators of prospectivity and potential mapping of porphyry deposits in middle and southern British Columbia by a GIS study of regional geochemical and other geoscientific data. Ph.D. thesis, University of Ottawa, Canada.
     [Google Scholar]
  15. DaubechiesI.1988. Orthonormal bases of compactly supported wavelets. Communications on Pure & Applied Mathematics41, 909–996.
     [Google Scholar]
  16. DaubechiesI., DefriseM. and De MolC.2004. An iterative thresholding algorithm for linear inverse problems with a sparsity constraint. Communications on Pure & Applied Mathematics57, 1413–1457.
     [Google Scholar]
  17. DavisK. and LiY.2011. Fast solution of geophysical inversion using adaptative mesh, space‐filling curves and wavelet compression. Geophysical Journal International185, 157–166.
     [Google Scholar]
  18. de Oliveira LyrioJ.C.S., TenorioL. and LiY.2004. Efficient automatic denoising of gravity gradiometry data. Geophysics69, 772–782.
     [Google Scholar]
  19. FediM. and QuartaT.1998. Wavelet analysis for the regional‐residual and local separation of potential field anomalies. Geophysical Prospecting46, 507–525.
     [Google Scholar]
  20. FediM. and FlorioG.2002. A stable downward continuation by using ISVD method. Geophysical Journal International151, 146–156.
     [Google Scholar]
  21. GholamiA. and SiahkoohiH.R.2009. Simultaneous constraining of model and data smoothness for regularization of geophysical inverse problems. Geophysical Journal International176, 151–163.
     [Google Scholar]
  22. GholamiA. and SacchiM.D.2012. A fast and automatic sparse decon‐volution in the presence of outliers. IEEE Transactions on Geoscience and Remote Sensing50, 4105–4116.
     [Google Scholar]
  23. GuoqingM., LiuC., HuangD. and LiL.2013. A stable iterative downward continuation of potential field data. Journal of Applied Geophysics, (in press).
     [Google Scholar]
  24. Jin‐ShengN., Hai‐HongW. and Zhi‐CaiL.2005. Downward Continuation of Gravity Signal Based on The Multiscale Edge Constraint. Chinese Journal of Geophysics48, 74–80.
     [Google Scholar]
  25. JohnD.A., AyusoR.A., BartonM.D., BlakelyR.J., BodnarR.J., DillesJ.H.et al.2010. Porphyry copper deposit model, chap. B of Mineral deposit models for resource assessment. U.S. Geological Survey Scientific Investigations Report 2010–5070–B, 169 p.
     [Google Scholar]
  26. LeblancG.E. and MorrisW.A.2001. Denoising of aeromagnetic data via the wavelet transform. Geophysics66, 1793–1804.
     [Google Scholar]
  27. LiY., DevrieseS.G.R., KrahenbuhlR. and DavisK.2013. Enhancement of Magnetic Data by Stable Downward Continuation for UXO Application. IEEE Transactions on Geoscience and Remote Sensing51, 3605–3614.
     [Google Scholar]
  28. LiY. and OldenburgD.W.2003. Fast inversion of large‐scale magnetic data using wavelet transforms and a logarithmic barrier method. Geophysical Journal International152, 251–265.
     [Google Scholar]
  29. MoreauF., GilbertD. and SaraccoG.1996. Filtering non stationary in geophysical data with orthogonal wavelets. Geophysical Research Letters23, 407–410.
     [Google Scholar]
  30. MoreauF., GibertD., HolschneiderM. and SaraccoG.1997. Wavelet analysis of potential fields. Inverse Problems13, 165–178.
     [Google Scholar]
  31. MoreauF., GilbertD., HolschneiderM. and SaraccoG.1999. Identification of potential fields with the continuous wavelet transform: basic theory. Journal of Geophysical Research104, 5003–5013.
     [Google Scholar]
  32. NamakiL., GholamiA. and HafiziM.A.2011. Edge‐preserved 2‐D inversion of magnetic data: an application to the Makran arc‐trench complex. Geophysical Journal International184, 1058–1068.
     [Google Scholar]
  33. PaolettiV., FediM., FlorioG. and RapollaA.2007. Localized cultural denoising of high‐resolution aeromagnetic data. Geophysical Prospecting55, 421–432.
     [Google Scholar]
  34. PaštekaR., KarcolR., KušnirákD. and MojzešA.2012. REGCONT: A Matlab based program for stable downward continuation of geophysical potential fields using Tikhonov regularization. Computers & Geosciences49, 278–289.
     [Google Scholar]
  35. PazandK., HezarkhaniA. and AtaeiM.2012. Using TOPSIS approaches for predictive porphyry Cu potential mapping: A case study in Ahar‐Arasbaran (NW‐Iran). Computers & Geosciences49, 62–71.
     [Google Scholar]
  36. RanjbarH. and HonarmandM.2004. Integration and analysis of airborne geophysical and ETM+ data for exploration of porphyry type deposits in the Central Iranian Volcanic Belt using fuzzy classification. International Journal of Remote Sensing25, 4729–4741.
     [Google Scholar]
  37. Ridsdill‐SmithT.A. and DentithM.C.1999. The wavelet transform in aeromagnetic processing. Geophysics64, 1003–1013.
     [Google Scholar]
  38. Ridsdill‐SmithT.A. and DentithM.C.2000. Drape corrections of aeromagnetic data using wavelets. Exploration Geophysics31, 39–46.
     [Google Scholar]
  39. Ridsdill‐SmithT.A.2000. The Application of the Wavelet Transform to the Processing of Aeromagnetic Data. PhD Thesis, The University of Western Australia, 179 p.
     [Google Scholar]
  40. SaraccoG., MoreauF., MathéP.‐E., HermitteD. and MichelJ.‐M.2007. Multiscale tomography of buried magnetic structures: Its use in the localization and characterization of archaeological structures. Geophysical Journal International171, 87–103.
     [Google Scholar]
  41. ThomanM.W., ZongeK.L. and LiuD.2000. Geophysical case history of North Silver Bell, Pima County, Arizona – A supergene‐enriched porphyry copper deposit. In: EllisR.B. , IrvineR. and FritzF. (eds), Northwest Mining Association 1998 Practical Geophysics Short Course Selected Papers on CD‐ROM.Spokane, Washington. Northwest Mining Association (paper 4, 42 p.).
     [Google Scholar]
  42. TikhonovA.N. and ArseninV.Y.1977. Solutions of Ill‐Posed Problems,Winston, Washington, D.C.
     [Google Scholar]
  43. TrompatH., BoschettiF. and HornbyP.2003. Improved downward continuation of potential field data. Exploration Geophysics34, 249–256.
     [Google Scholar]
  44. ValeeM.A., KeatingP., SmithR.S. and St‐HilaireC.2004. Estimating depth and model type using the continuous wavelet transform of magnetic data. Geophysics69, 191–199.
     [Google Scholar]
  45. WahbaG.1977. Practical approximate solution of linear operator equations when the data are noisy. SIAM Journal of Numerical Analysis14, 651–667.
     [Google Scholar]
  46. XuS., YangJ., YangC., XiaoP., ChenS. and GuoZ.2007. The iteration method for downward continuation of a potential field from a horizontal plane. Geophysical Prospecting55, 883–889.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.3997/1873-0604.2014027
Loading
A new stable downward continuation of airborne magnetic data based on Wavelet deconvolution
Near Surface Geophysics 12, 751 (2014); https://doi.org/10.3997/1873-0604.2014027
/content/journals/10.3997/1873-0604.2014027
/content/journals/10.3997/1873-0604.2014027
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error