The problem of hollow space correction in subsurface gravity surveys

Geophysics ◽  
1983 ◽  
Vol 48 (10) ◽  
pp. 1406-1408
Author(s):  
B. V. Satyanarayana Murty ◽  
P. Chandra Reddy

In spite of the extensive computational work involved in the reduction of data especially in estimating the effects of terrain, the gravity method has earned its own esteemed place in mineral exploration. Measurement of gravity in subsurface openings such as mine shafts, drives, and adits is much more valuable in mineral exploration and also useful in the planning, design, and maintenance of mines. However, because the anomalies being sought in this context are usually of very small magnitudes, accuracy at every stage of the gravity investigation is essential. Besides those corrections which are known for the surface gravity measurements, an additional correction for the subsurface gravity data is the result of hollow spaces.

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 925-937 ◽  
Author(s):  
Michael S. Zhdanov ◽  
Robert Ellis ◽  
Souvik Mukherjee

We develop a new method for interpretation of tensor gravity field component data, based on regularized focusing inversion. The focusing inversion makes its possible to reconstruct a sharper image of the geological target than conventional maximum smoothness inversion. This new technique can be efficiently applied for the interpretation of gravity gradiometer data, which are sensitive to local density anomalies. The numerical modeling and inversion results show that the resolution of the gravity method can be improved significantly if we use tensor gravity data for interpretation. We also apply our method for inversion of the gradient gravity data collected by BHP Billiton over the Cannington Ag‐Pb‐Zn orebody in Queensland, Australia. The comparison with the drilling results demonstrates a remarkable correlation between the density anomaly reconstructed by the gravity gradient data and the true structure of the orebody. This result indicates that the emerging new geophysical technology of the airborne gravity gradient observations can improve significantly the practical effectiveness of the gravity method in mineral exploration.


2021 ◽  
Vol 95 (2) ◽  
Author(s):  
Mirjam Bilker-Koivula ◽  
Jaakko Mäkinen ◽  
Hannu Ruotsalainen ◽  
Jyri Näränen ◽  
Timo Saari

AbstractPostglacial rebound in Fennoscandia causes striking trends in gravity measurements of the area. We present time series of absolute gravity data collected between 1976 and 2019 on 12 stations in Finland with different types of instruments. First, we determine the trends at each station and analyse the effect of the instrument types. We estimate, for example, an offset of 6.8 μgal for the JILAg-5 instrument with respect to the FG5-type instruments. Applying the offsets in the trend analysis strengthens the trends being in good agreement with the NKG2016LU_gdot model of gravity change. Trends of seven stations were found robust and were used to analyse the stabilization of the trends in time and to determine the relationship between gravity change rates and land uplift rates as measured with global navigation satellite systems (GNSS) as well as from the NKG2016LU_abs land uplift model. Trends calculated from combined and offset-corrected measurements of JILAg-5- and FG5-type instruments stabilized in 15 to 20 years and at some stations even faster. The trends of FG5-type instrument data alone stabilized generally within 10 years. The ratio between gravity change rates and vertical rates from different data sets yields values between − 0.206 ± 0.017 and − 0.227 ± 0.024 µGal/mm and axis intercept values between 0.248 ± 0.089 and 0.335 ± 0.136 µGal/yr. These values are larger than previous estimates for Fennoscandia.


2017 ◽  
Vol 5 (2) ◽  
pp. T243-T257 ◽  
Author(s):  
Martin Landrø ◽  
Mark Zumberge

We have developed a calibrated, simple time-lapse seismic method for estimating saturation changes from the [Formula: see text]-storage project at Sleipner offshore Norway. This seismic method works well to map changes when [Formula: see text] is migrating laterally away from the injection point. However, it is challenging to detect changes occurring below [Formula: see text] layers that have already been charged by some [Formula: see text]. Not only is this partly caused by the seismic shadow effects, but also by the fact that the velocity sensitivity for [Formula: see text] change in saturation from 0.3 to 1.0 is significantly less than saturation changes from zero to 0.3. To circumvent the seismic shadow zone problem, we combine the time-lapse seismic method with time-lapse gravity measurements. This is done by a simple forward modeling of gravity changes based on the seismically derived saturation changes, letting these saturation changes be scaled by an arbitrary constant and then by minimizing the least-squares error to obtain the best fit between the scaled saturation changes and the measured time-lapse gravity data. In this way, we are able to exploit the complementary properties of time-lapse seismic and gravity data.


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
K. Breili ◽  
R. Hougen ◽  
D. I. Lysaker ◽  
O. C. D. Omang ◽  
B. Tangen

AbstractThe Norwegian Mapping Authority (NMA) has recently established a new gravity laboratory in Ny-Ålesund at Svalbard, Norway. The laboratory consists of three independent pillars and is part of the geodetic core station that is presently under construction at Brandal, approximately 1.5 km north of NMA’s old station. In anticipation of future use of the new gravity laboratory, we present benchmark gravity values, gravity gradients, and final coordinates of all new pillars. Test measurements indicate a higher noise level at Brandal compared to the old station. The increased noise level is attributed to higher sensitivity to wind.We have also investigated possible consequences of moving to Brandal when it comes to the gravitational signal of present-day ice mass changes and ocean tide loading. Plausible models representing ice mass changes at the Svalbard archipelago indicate that the gravitational signal at Brandal may differ from that at the old site with a size detectable with modern gravimeters. Users of gravity data from Ny-Ålesund should, therefore, be cautious if future observations from the new observatory are used to extend the existing gravity record. Due to its lower elevation, Brandal is significantly less sensitive to gravitational ocean tide loading. In the future, Brandal will be the prime site for gravimetry in Ny-Ålesund. This ensures gravity measurements collocated with space geodetic techniques like VLBI, SLR, and GNSS.


2020 ◽  
Vol 55 (3) ◽  
pp. 100-117
Author(s):  
Viktor Szabó ◽  
Dorota Marjańska

AbstractGlobal satellite gravity measurements provide unique information regarding gravity field distribution and its variability on the Earth. The main cause of gravity changes is the mass transportation within the Earth, appearing as, e.g. dynamic fluctuations in hydrology, glaciology, oceanology, meteorology and the lithosphere. This phenomenon has become more comprehensible thanks to the dedicated gravimetric missions such as Gravity Recovery and Climate Experiment (GRACE), Challenging Minisatellite Payload (CHAMP) and Gravity Field and Steady-State Ocean Circulation Explorer (GOCE). From among these missions, GRACE seems to be the most dominating source of gravity data, sharing a unique set of observations from over 15 years. The results of this experiment are often of interest to geodesists and geophysicists due to its high compatibility with the other methods of gravity measurements, especially absolute gravimetry. Direct validation of gravity field solutions is crucial as it can provide conclusions concerning forecasts of subsurface water changes. The aim of this work is to present the issue of selection of filtration parameters for monthly gravity field solutions in RL06 and RL05 releases and then to compare them to a time series of absolute gravimetric data conducted in quasi-monthly measurements in Astro-Geodetic Observatory in Józefosław (Poland). The other purpose of this study is to estimate the accuracy of GRACE temporal solutions in comparison with absolute terrestrial gravimetry data and making an attempt to indicate the significance of differences between solutions using various types of filtration (DDK, Gaussian) from selected research centres.


Author(s):  
K. Jourde ◽  
D. Gibert ◽  
J. Marteau

Abstract. This paper examines how the resolution of small-scale geological density models is improved through the fusion of information provided by gravity measurements and density muon radiographies. Muon radiography aims at determining the density of geological bodies by measuring their screening effect on the natural flux of cosmic muons. Muon radiography essentially works like medical X-ray scan and integrates density information along elongated narrow conical volumes. Gravity measurements are linked to density by a 3-D integration encompassing the whole studied domain. We establish the mathematical expressions of these integration formulas – called acquisition kernels – and derive the resolving kernels that are spatial filters relating the true unknown density structure to the density distribution actually recovered from the available data. The resolving kernels approach allows to quantitatively describe the improvement of the resolution of the density models achieved by merging gravity data and muon radiographies. The method developed in this paper may be used to optimally design the geometry of the field measurements to perform in order to obtain a given spatial resolution pattern of the density model to construct. The resolving kernels derived in the joined muon/gravimetry case indicate that gravity data are almost useless to constrain the density structure in regions sampled by more than two muon tomography acquisitions. Interestingly the resolution in deeper regions not sampled by muon tomography is significantly improved by joining the two techniques. The method is illustrated with examples for La Soufrière of Guadeloupe volcano.


Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 807-816 ◽  
Author(s):  
Jérôme Verdun ◽  
Roger Bayer ◽  
Emile E. Klingelé ◽  
Marc Cocard ◽  
Alain Geiger ◽  
...  

This paper introduces a new approach to airborne gravity data reduction well‐suited for surveys flown at high altitude with respect to gravity sources (mountainous areas). Classical technique is reviewed and illustrated in taking advantage of airborne gravity measurements performed over the western French Alps by using a LaCoste & Romberg air‐sea gravity meter. The part of nongravitational vertical accelerations correlated with gravity meter measurements are investigated with the help of coherence spectra. Beam velocity has proved to be strikingly correlated with vertical acceleration of the aircraft. This finding is theoretically argued by solving the equation of the gravimetric system (gravity meter and stabilized platform). The transfer function of the system is derived, and a new formulation of airborne gravity data reduction, which takes care of the sensitive response of spring tension to observable gravity field wavelengths, is given. The resulting gravity signal exhibits a residual noise caused by electronic devices and short‐wavelength Eötvös effects. The use of dedicated exponential filters gives us a way to eliminate these high‐frequency effects. Examples of the resulting free‐air anomaly at 5100‐m altitude along one particular profile are given and compared with free‐air anomaly deduced from the classical method for processing airborne gravity data, and with upward‐continued ground gravity data. The well‐known trade‐off between accuracy and resolution is discussed in the context of a mountainous area.


2012 ◽  
Vol 19 (2) ◽  
pp. 291-296 ◽  
Author(s):  
M. Pilkington ◽  
P. Keating

Abstract. Most interpretive methods for potential field (magnetic and gravity) measurements require data in a gridded format. Many are also based on using fast Fourier transforms to improve their computational efficiency. As such, grids need to be full (no undefined values), rectangular and periodic. Since potential field surveys do not usually provide data sets in this form, grids must first be prepared to satisfy these three requirements before any interpretive method can be used. Here, we use a method for grid preparation based on a fractal model for predicting field values where necessary. Using fractal field values ensures that the statistical and spectral character of the measured data is preserved, and that unwanted discontinuities at survey boundaries are minimized. The fractal method compares well with standard extrapolation methods using gridding and maximum entropy filtering. The procedure is demonstrated on a portion of a recently flown aeromagnetic survey over a volcanic terrane in southern British Columbia, Canada.


1999 ◽  
Vol 36 (1) ◽  
pp. 75-89 ◽  
Author(s):  
Hamid Telmat ◽  
Jean-Claude Mareschal ◽  
Clément Gariépy

Gravity data were obtained along two transects on the southern coast of Ungava Bay, which provide continuous gravity coverage between Leaf Bay and George River. The transects and the derived gravity profiles extend from the Superior craton to the Rae Province across the New Quebec Orogen (NQO). Interpretation of the transect along the southwestern coast of Ungava Bay suggests crustal thickening beneath the NQO and crustal thinning beneath the Kuujjuaq Terrane, east of the NQO. Two alternative interpretations are proposed for the transect along the southeastern coast of the bay. The first model shows crustal thickening beneath the George River Shear Zone (GRSZ) and two shallow bodies correlated with the northern extensions of the GRSZ and the De Pas batholith. The second model shows constant crustal thickness and bodies more deeply rooted than in the first model. The gravity models are consistent with the easterly dipping reflections imaged along a Lithoprobe seismic line crossing Ungava Bay and suggest westward thrusting of the Rae Province over the NQO. Because no gravity data have been collected in Ungava Bay, satellite altimetry data have been used as a means to fill the gap in data collected at sea. The satellite-derived gravity data and standard Bouguer gravity data were combined in a composite map for the Ungava Bay region. The new land-based gravity measurements were used to verify and calibrate the satellite data and to ensure that offshore gravity anomalies merge with those determined by the land surveys in a reasonable fashion. Three parallel east-west gravity profiles were extracted: across Ungava Bay (59.9°N), on the southern shore of the bay (58.5°N), and onshore ~200 km south of Ungava Bay (57.1°N). The gravity signature of some major structures, such as the GRSZ, can be identified on each profile.


Sign in / Sign up

Export Citation Format

Share Document