On the use of the profiled singular-function expansion in gravity gradiometry

Abstract

We demonstrate a method of inverting gravity data, based on profiled singular-function expansions. It is well known that the inverse problem of determining underground density variations from gravity data is severely ill-posed and prior knowledge is needed to restrict the range of possible solutions. Viewed as a linear inverse problem, various standard methods, all of which produce solutions approximating the generalised solution, tend to give density variations concentrated near the surface. To overcome this potentially undesirable trait, various authors have introduced depth weighting. In this paper we carry this idea a step further and introduce a method based on a profiled singular-function expansion which uses the prior knowledge that the underground object is centred at a particular depth. The use of an appropriate depth-weighting profile leads to a solution at the correct depth. Furthermore, we show that if the depth of the object is unknown, a range of solutions at different depths can be produced, allowing other prior knowledge, such as object size or density, to be introduced to determine which solution is the most plausible. The truncation point of the profiled singular-function expansion is determined by the level of noise on the data. We examine how the achievable horizontal resolution varies with this truncation point. A notable property of our approach is that when the centre of the profile corresponds to the true object depth, the solution appears, in a certain sense, to be the most focussed one. Finally we consider a gravimetry example using real data.

Quantum sensing for gravity cartography

“…3c ). The observed spread of the POE is expected, due to measurement uncertainty and the ambiguity that exists between model parameters specifying depth, area and density, typical of inference from potential field data 37 . A signal-to-noise ratio of 8 for the detection is estimated from the data, finding the deduced horizontal position of the tunnel centre at (0.19 ± 0.19) m along the survey line and a depth to the centre of (1.89 −0.59/+2.3) m (see Methods ).…”

Quantum sensing for gravitational cartography

“…3c). The observed spread of the POE is expected, due to measurement uncertainty and the ambiguity that exists between model parameters specifying depth, area, and density, typical of inference from potential eld data 37 . A signal to noise ratio of 8 for the detection is estimated from the data, nding the deduced horizontal position of the tunnel centre at (0.19 ± 0.19) m along the survey line and a depth to the centre of (1.89 -0.59/+2.3) m (see Methods).…”

Advanced Potential Field Analysis Applied to Mise‐à‐la‐Masse Surveys for Leakage Detection

“…A number of different approaches have been proposed and used to this end, in particular in search for the depth of localized anomaly sources (Binley et al., 1999; Ling et al., 2019; Mary et al., 2020; Peruzzo et al., 2020; Shao et al., 2018; Wondimu et al., 2018). A priori information is always needed to guide the result (de Villiers et al., 2019) and a correction for the influence of the resistivity on the current density distribution is needed.…”

“…is always needed to guide the result (de Villiers et al, 2019) and a correction for the influence of the resistivity on the current density distribution is needed.…”