3-D probability density imaging of Euler solutions using gravity data: a case study of Mount Milligan, Canada

• Published in: Acta geophysica Vol. 72; no. 5; pp. 3371 - 3391
• Main Authors: Cao, Shujin, Deng, Yihuai, Yang, Bo, Lu, GuangYin, Zhu, Ziqiang, Chen, Peng, Xie, Jingtao, Chen, Xinyue
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 2024; Springer Nature B.V
• Subjects: Algorithms; Deconvolution; Density distribution; Earth and Environmental Science; Earth Sciences; Geophysics/Geodesy; Geotechnical Engineering & Applied Earth Sciences; Multilayers; Parameter estimation; Potential fields; Research Article - Applied Geophysics; Structural Geology; Tensors; Canada; B-spline density estimation; Euler deconvolution; Probability density imaging; Full tensor gravity gradient
• ISSN: 1895-7455; 1895-6572; 1895-7455
• DOI: 10.1007/s11600-023-01279-y

Euler deconvolution is a widely used automatic or semi-automatic method for potential field data. However, it yields many spurious solutions that complicate interpretation and must be reduced, eliminated, recognized, or ignored during interpretation. This study proposes a post-processing algorithm that converts Euler solutions produced by tensor Euler deconvolution of gravity data with an unprescribed structural index into probability values ( p values) using the B-spline series density estimation (BSS) method. The p values of the Euler solution set form a probability density distribution on the estimation grid. The BSS method relies on the fact that while spurious solutions are sparse and ubiquitous, Euler deconvolution yields many similar or duplicate solutions, which may tightly cluster near real sources. The p values of the Euler solution clusters form multi-layered isosurfaces that can be used to discriminate neighboring target sources because the p values of spurious solutions are vanishingly small, making it simple to remove their interference from the probability density distribution. In all synthetic cases, the geometric outlines of anomaly sources are estimated from probability density isosurfaces approximating synthetic model parameters. The BSS method was then applied to airborne gravity data from Mount Milligan, British Columbia, Canada. Subsequently, results from synthetic models and field data show that the proposed method can successfully localize meaningful geological targets.