A solution to unexploded ordnance detection problem from its magnetic anomaly using Kaczmarz regularization

2016 ◽  
Vol 4 (3) ◽  
pp. SH61-SH69 ◽  
Author(s):  
Maha Abdelazeem ◽  
Mohamed M. Gobashy

Old military events pose a critical and severe problem for many countries, including the Egyptian northwestern coast. These result in extensive areas of surface/subsurface landmines that affect the economic development plans of many countries. Detection of these landmines becomes a target for many geophysical research teams. Currently, unconventional near-surface flight technologies, such as quad-hexacopters instead of regular land surveys, are used for safety reasons in the acquisition phase. We have introduced a new processing and modeling technique of magnetic data conducted over mines or near-surface geophysical targets for accurate and precise determination of location and depth. The technique is based on the application of the Kaczmarz regularization method to the ill-posed magnetic inverse problem. The advantage of this method is the optimum transformation of regularized normal equations to an equivalent augmented regularized normal system of equations. The condition number of the updated system, which determines the degree of ill posedness, is greatly lower than the original one; this improves and guarantees a good solution to the system. The method is applied to an unexploded ordnance (UXO) test site in the United Kingdom. Our results have determined that the technique is appropriate and promising in efficiently addressing a wide number of problems that are important to near-surface geophysicists, including UXO detection.

2021 ◽  
Vol 51 (3) ◽  
pp. 277-294
Author(s):  
Roman PAŠTEKA ◽  
Miroslav HAJACH ◽  
Bibiana BRIXOVÁ ◽  
Ján MIKUŠKA ◽  
John STANLEY

In this contribution we present results from a case-study, which was performed in collaboration between geophysicists and explosive ordnance disposal technicians at the Rohožník military training range in SW Slovakia. The aim of this study was to locate a deep-penetrated unexploded Mk-82 aerial bomb using high-definition digital magnetometry. The location where this bomb had entered the ground was known but its final position needed to be determined so that a safe excavation and disposal could be conducted. However, the detection of this unexploded ordnance object was complicated by the presence of intense magnetic interference from a number of near surface ferrous items including non-explosive test bombs, fragmentation and other iron junk. These items contributed a localised, high amplitude of magnetic clutter masking any deeper source. Our strategy was to approach the problem in three stages. First, we used magnetic data to locate the near surface items. After the detection and before the excavation of the searched objects, two quantitative interpretation methods were used. These involved an optimised modelling of source bodies and the application of a 3D Euler deconvolution. Both methods yielded acceptable results, but the former was found to be more accurate. After the interpretation phase, many of the items were then safely excavated and removed individually. A second magnetic mapping was then performed and from this data which was now significantly less cluttered, we were able to identify but not quantify, two deep source items and to confirm that all remaining near surface items were significantly smaller in size than a Mk-82 bomb. As the remaining near surface sources were interpreted as being contained within the surface one metre of soil and being small they could be assured to be non-explosive, it was considered most practical to mechanically excavate and remove this soil and the remaining objects contained.


2014 ◽  
Vol 644-650 ◽  
pp. 2670-2673
Author(s):  
Jun Wang ◽  
Xiao Hong Meng ◽  
Fang Li ◽  
Jun Jie Zhou

With the continuing growth in influence of near surface geophysics, the research of the subsurface structure is of great significance. Geophysical imaging is one of the efficient computer tools that can be applied. This paper utilize the inversion of potential field data to do the subsurface imaging. Here, gravity data and magnetic data are inverted together with structural coupled inversion algorithm. The subspace (model space) is divided into a set of rectangular cells by an orthogonal 2D mesh and assume a constant property (density and magnetic susceptibility) value within each cell. The inversion matrix equation is solved as an unconstrained optimization problem with conjugate gradient method (CG). This imaging method is applied to synthetic data for typical models of gravity and magnetic anomalies and is tested on field data.


1969 ◽  
Vol 59 (6) ◽  
pp. 2271-2281
Author(s):  
R. M. Hamilton ◽  
J. H. Healy

abstract The Benham nuclear explosion, a 1.1 megaton test 1.4 km beneath Pahute Mesa at the Nevada Test Site, initiated a sequence of earthquakes lasting several months. The epicenters of these shocks were located within 13 km of ground zero in several linear zones that parallel the regional fault trends. Focal depths range from near surface to 6 km. The earthquakes are not located in the zone of the major ground breakage. The earthquake distribution and fault plane solutions together indicate that both right-lateral strike-slip fault movement and dip-slip fault movement occurred. The explosion apparently caused the release of natural tectonic strain.


Sensors ◽  
2019 ◽  
Vol 19 (1) ◽  
pp. 170 ◽  
Author(s):  
Xianyang Gao ◽  
Frank J. W. Podd ◽  
Wouter Van Verre ◽  
David J. Daniels ◽  
Anthony J. Peyton

Antennas are an important component in ground penetrating radar (GPR) systems. Although there has been much research reported on the design of individual antennas, there is less research reported on the design of the geometry of bi-static antennas. This paper considers the effects of key parameters in the setup of a GPR head consisting of a bi-static bow-tie pair to show the effect of these parameters on the GPR performance. The parameters investigated are the antenna separation, antenna height above the soil, and antenna input impedance. The investigation of the parameters was performed by simulation and measurements. It was found when the bi-static antennas were separated by 7 cm to 9 cm and were operated close to the soil (2 cm to 4 cm), the reflected signal from a near-surface object is relatively unaffected by height variation and object depth. An antenna input impedance of 250 Ω was chosen to feed the antennas to reduce the late-time ringing. Using these results, a new GPR system was designed and then evaluated at a test site near Benkovac, Croatia.


Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. A25-A28 ◽  
Author(s):  
Elena Pettinelli ◽  
Giuliano Vannaroni ◽  
Barbara Di Pasquo ◽  
Elisabetta Mattei ◽  
Andrea Di Matteo ◽  
...  

We explore a new approach to evaluate the effect of soil electromagnetic parameters on early-time ground-penetrating radar (GPR) signals. The analysis is performed in a time interval which contains the direct airwaves and ground waves, propagating between transmitting and receiving antennas. To perform the measurements we have selected a natural test site characterized by very strong lateral gradient of the soil electrical properties. To evaluate the effect of the subsoil permittivity and conductivity on the radar response we compare the envelope amplitude of the GPR signals received in the first [Formula: see text] within [Formula: see text]-wide windows, with the electrical properties ([Formula: see text] and [Formula: see text]) determined using time-domain reflectometry (TDR). The results show that the constitutive soil parameters strongly influence early-time signals, suggesting a novel approach for estimating the spatial variability of water content with GPR.


Geophysics ◽  
1993 ◽  
Vol 58 (5) ◽  
pp. 645-652 ◽  
Author(s):  
Derek J. Woodward

Although draped magnetic surveys contain more information about the magnetization of the rocks near the surface of the earth than surveys at constant elevation, allowance for the effects of the terrain is critical for their correct interpretation. A new method for calculating the magnetic effect of the topography from a digital terrain model by integrating analytically in the vertical direction and then numerically in the horizontal plane is presented. This method lends itself to the calculation of anomalies when the magnetization of the rocks varies with position and thus is well suited to the inversion of draped aeromagnetic surveys to obtain the apparent magnetization of the surficial rocks. This inversion is achieved by repeated use of an approximate inverse function in the form of a two‐dimensional (2-D) filter that is applied to gridded data. An example, using draped magnetic data collected over White Island, an active volcanic island of high relief, shows that although the anomaly pattern is dominated by topographic effects, the distribution of near‐surface magnetic bodies can be determined by a joint inversion of the data and the topography. One of the highly magnetized areas of White Island is interestingly in the vicinity of the active crater, with another near the inner wall of the caldera where there are numerous fumaroles. It may be expected that the higher temperatures in these areas would reduce the magnetization. However, it appears that an explanation for the higher magnetization can be found in the stability field of the mineral magnetite.


Geophysics ◽  
2020 ◽  
pp. 1-41 ◽  
Author(s):  
Jens Tronicke ◽  
Niklas Allroggen ◽  
Felix Biermann ◽  
Florian Fanselow ◽  
Julien Guillemoteau ◽  
...  

In near-surface geophysics, ground-based mapping surveys are routinely employed in a variety of applications including those from archaeology, civil engineering, hydrology, and soil science. The resulting geophysical anomaly maps of, for example, magnetic or electrical parameters are usually interpreted to laterally delineate subsurface structures such as those related to the remains of past human activities, subsurface utilities and other installations, hydrological properties, or different soil types. To ease the interpretation of such data sets, we propose a multi-scale processing, analysis, and visualization strategy. Our approach relies on a discrete redundant wavelet transform (RWT) implemented using cubic-spline filters and the à trous algorithm, which allows to efficiently compute a multi-scale decomposition of 2D data using a series of 1D convolutions. The basic idea of the approach is presented using a synthetic test image, while our archaeo-geophysical case study from North-East Germany demonstrates its potential to analyze and process rather typical geophysical anomaly maps including magnetic and topographic data. Our vertical-gradient magnetic data show amplitude variations over several orders of magnitude, complex anomaly patterns at various spatial scales, and typical noise patterns, while our topographic data show a distinct hill structure superimposed by a microtopographic stripe pattern and random noise. Our results demonstrate that the RWT approach is capable to successfully separate these components and that selected wavelet planes can be scaled and combined so that the reconstructed images allow for a detailed, multi-scale structural interpretation also using integrated visualizations of magnetic and topographic data. Because our analysis approach is straightforward to implement without laborious parameter testing and tuning, computationally efficient, and easily adaptable to other geophysical data sets, we believe that it can help to rapidly analyze and interpret different geophysical mapping data collected to address a variety of near-surface applications from engineering practice and research.


Sign in / Sign up

Export Citation Format

Share Document