Correct structural index defined by base level estimates in Euler deconvolution

2017 ◽  
Author(s):  
Felipe Ferreira de Melo ◽  
Valeria Cristina Ferreira Barbosa
Keyword(s):  
Euler Deconvolution ◽  
Structural Index ◽  
Base Level

Estimating the nature and the horizontal and vertical positions of 3D magnetic sources using Euler deconvolution

Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. J87-J98 ◽  
Author(s):  
Felipe F. Melo ◽  
Valeria C. F. Barbosa ◽  
Leonardo Uieda ◽  
Vanderlei C. Oliveira Jr. ◽  
João B. C. Silva
Keyword(s):  
Synthetic Data ◽  
Real Data ◽  
Angular Coefficient ◽  
Coefficient Estimates ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Vertical Coordinate ◽  
Base Level ◽  
The One ◽  
Inclined Planes

We have developed a new method that drastically reduces the number of the source location estimates in Euler deconvolution to only one per anomaly. Our method employs the analytical estimators of the base level and of the horizontal and vertical source positions in Euler deconvolution as a function of the [Formula: see text]- and [Formula: see text]-coordinates of the observations. By assuming any tentative structural index (defining the geometry of the sources), our method automatically locates plateaus, on the maps of the horizontal coordinate estimates, indicating consistent estimates that are very close to the true corresponding coordinates. These plateaus are located in the neighborhood of the highest values of the anomaly and show a contrasting behavior with those estimates that form inclined planes at the anomaly borders. The plateaus are automatically located on the maps of the horizontal coordinate estimates by fitting a first-degree polynomial to these estimates in a moving-window scheme spanning all estimates. The positions where the angular coefficient estimates are closest to zero identify the plateaus of the horizontal coordinate estimates. The sample means of these horizontal coordinate estimates are the best horizontal location estimates. After mapping each plateau, our method takes as the best structural index the one that yields the minimum correlation between the total-field anomaly and the estimated base level over each plateau. By using the estimated structural index for each plateau, our approach extracts the vertical coordinate estimates over the corresponding plateau. The sample means of these estimates are the best depth location estimates in our method. When applied to synthetic data, our method yielded good results if the bodies produce weak- and mid-interfering anomalies. A test on real data over intrusions in the Goiás Alkaline Province, Brazil, retrieved sphere-like sources suggesting 3D bodies.

Stability analysis and improvement of structural index estimation in Euler deconvolution

Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 48-60 ◽  
Author(s):  
Valéria C. F. Barbosa ◽  
João B. C. Silva ◽  
Walter E. Medeiros
Keyword(s):  
Stability Analysis ◽  
Vertical Position ◽  
Horizontal Position ◽  
Total Field ◽  
Source Position ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Data Window ◽  
Base Level ◽  
New Criterion

Euler deconvolution has been widely used in automatic aeromagnetic interpretations because it requires no prior knowledge of the source magnetization direction and assumes no particular interpretation model, provided the structural index defining the anomaly falloff rate related to the nature of the magnetic source, is determined in advance. Estimating the correct structural index and electing optimum criteria for selecting candidate solutions are two fundamental requirements for a successful application of this method. We present a new criterion for determining the structural index. This criterion is based on the correlation between the total‐field anomaly and the estimates of an unknown base level. These estimates are obtained for each position of a moving data window along the observed profile and for several tentative values for the structural index. The tentative value for the structural index producing the smallest correlation is the best estimate of the correct structural index. We also propose a new criterion to select the best solutions from a set of previously computed candidate solutions, each one associated with a particular position of the moving data window. A current criterion is to select only those candidates producing a standard deviation for the vertical position of the source smaller than a threshold value. We propose that in addition to this criterion, only those candidates producing the best fit to the known quantities (combinations of anomaly and its gradients) be selected. The proposed modifications to Euler deconvolution can be implemented easily in an automated algorithm for locating the source position. The above results are grounded on a theoretical uniqueness and stability analysis, also presented in this paper, for the joint estimation of the source position, the base level, and the structural index in Euler deconvolution. This analysis also reveals that the vertical position and the structural index of the source cannot be estimated simultaneously because they are linearly dependent; the horizontal position and the structural index, on the other hand, are linearly independent. For a known structural index, estimates of both horizontal and vertical positions are unique and stable regardless of the value of the structural index. If this value is not too small, estimates of the base level for the total field are stable as well. The proposed modifications to Euler deconvolution were tested both on synthetic and real magnetic data. In the case of synthetic data, the proposed criterion always detected the correct structural index and good estimates of the source position were obtained, suggesting the present theoretical analysis may lead to a substantial enhancement in practical applications of Euler deconvolution. In the case of practical data (vertical component anomaly over an iron deposit in the Kursk district, Russia), the estimated structural index (corresponding to a vertical prism) was in accordance with the known geology of the deposit, and the estimates of the depth and horizontal position of the source compared favorably with results reported in the literature.

Correct structural index in Euler deconvolution via base-level estimates

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. J87-J98 ◽  
Author(s):  
Felipe F. Melo ◽  
Valéria C. F. Barbosa
Keyword(s):  
Standard Deviation ◽  
Field Measurements ◽  
Real Data ◽  
Magnetic Anomalies ◽  
Magnetic Data ◽  
Total Field ◽  
Minimum Standard ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Base Level

In most applications, the Euler deconvolution aims to define the nature (type) of the geologic source (i.e., the structural index [SI]) and its depth position. However, Euler deconvolution also estimates the horizontal positions of the sources and the base level of the magnetic anomaly. To determine the correct SI, most authors take advantage of the clustering of depth estimates. We have analyzed Euler’s equation to indicate that random variables contaminating the magnetic observations and its gradients affect the base-level estimates if, and only if, the SI is not assumed correctly. Grounded on this theoretical analysis and assuming a set of tentative SIs, we have developed a new criterion for determining the correct SI by means of the minimum standard deviation of base-level estimates. We performed synthetic tests simulating multiple magnetic sources with different SIs. To produce mid and strongly interfering synthetic magnetic anomalies, we added constant and nonlinear backgrounds to the anomalies and approximated the simulated sources laterally. If the magnetic anomalies are weakly interfering, the minima standard deviations either of the depth or base-level estimates can be used to determine the correct SI. However, if the magnetic anomalies are strongly interfering, only the minimum standard deviation of the base-level estimates can determine the SI correctly. These tests also show that Euler deconvolution does not require that the magnetic data be corrected for the regional fields (e.g., International Geomagnetic Reference Field [IGRF]). Tests on real data from part of the Goiás Alkaline Province, Brazil, confirm the potential of the minimum standard deviation of base-level estimates in determining the SIs of the sources by applying Euler deconvolution either to total-field measurements or to total-field anomaly (corrected for IGRF). Our result suggests three plug intrusions giving rise to the Diorama anomaly and dipole-like sources yielding Arenópolis and Montes Claros de Goiás anomalies.

Effects of low-pass filtering on the calculated structural index from magnetic data

Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. L23-L28 ◽  
Author(s):  
Kristofer Davis ◽  
Yaoguo Li ◽  
Misac N. Nabighian
Keyword(s):  
Magnetic Data ◽  
Cutoff Wavelength ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Pass Filter ◽  
Low Pass Filter ◽  
Source Type ◽  
A Value ◽  
Low Pass ◽  
Unexploded Ordnance Detection

Euler and extended Euler deconvolution applications use an assumed structural index (SI) or calculate the SI, respectively, for magnetic anomaly data within a specified window. The structural index depends on the source type: specifically, the rate at which the field produced by the source decays. We have examined the effects that the application of low-pass filtering to magnetic data has on estimating the SI. Using a simple low-pass filter, we derived the SI for filtered-field solutions directly over, and away from, a target based on the magnetic potential of a vertical dipole [Formula: see text]. We validated this approach by applying extended Euler deconvolution to synthetic and field examples. In general, filtered magnetic data will decrease the numerically determined SI to a value lower than the theoretical one. The slope and cutoff wavelength of the filter directly affect the estimated SI solutions. The results prove that one must take into account filtering for the application of Euler deconvolution to locate dipole anomalies for unexploded ordnance detection.

Scattering, symmetry, and bias analysis of source‐position estimates in Euler deconvolution and its practical implications

Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 1149-1156 ◽  
Author(s):  
J. B. C. Silva ◽  
V. C. F. Barbosa ◽  
W. E. Medeiros
Keyword(s):  
Vertical Position ◽  
Lower Sensitivity ◽  
Horizontal Position ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Bias Analysis ◽  
Data Noise ◽  
Data Window ◽  
Symmetry Properties ◽  
Window Position

Despite being widely used, Euler deconvolution has received little attention in its theoretical aspects. We obtain analytical estimators [Formula: see text] and [Formula: see text] for the horizontal [Formula: see text] and the vertical [Formula: see text] source positions, which are the parameters to be estimated when the structural index is assigned a tentative value. By analyzing the estimators’ properties, we show two things. (1) Scattering of the alternative solutions associated with different data window positions has two components: data noise and the use of a wrong tentative value for the structural index, (2) The differences [Formula: see text] and [Formula: see text], as a function of the data window position, exhibit antisymmetric and symmetric behaviors, respectively, about [Formula: see text] when the magnetization and the geomagnetic field inclinations are 90°, 0°, or 45°. For intermediate inclinations there is a moderate departure from the symmetry properties. This analysis shows that the criterion for determining the structural index as the tentative value producing the smallest solution scattering is theoretically sound but occasionally fails in practice because data noise also contributes to solution scattering. In addition, the lower sensitivity of the averages of estimates of the horizontal position (as compared with the vertical position) occurs because [Formula: see text] displays an antisymmetric behavior with respect to the true horizontal position while the estimates of the vertical position display a symmetric behavior with respect to the same point.

Euler deconvolution of gravity anomalies from thick contact/fault structures with extended negative structural index

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. I51-I58 ◽  
Author(s):  
Petar Stavrev ◽  
Alan Reid
Keyword(s):  
Field Data ◽  
Gravity Anomalies ◽  
Potential Fields ◽  
Contact Structures ◽  
Euler Deconvolution ◽  
Structural Index ◽  
Infinite Slab ◽  
Fault Structures ◽  
Structural Indices ◽  
Analytical Expressions

The concept of extended Euler homogeneity of potential fields is examined with respect to all variables of length dimension in their analytical expressions. This reveals the possible existence of positive degrees of homogeneity or corresponding negative structural indices considered as extensions of the Thompson’s structural indices in Euler deconvolution. This approach is implemented for a contact gravity model, represented by a 2D semi-infinite slab with large thickness relative to its depth. Applying Euler deconvolution on synthetic and field data indicates that the positive degree of homogeneity, i.e., the extended negative structural index, is the appropriate one for the inversion of gravity anomalies from contact structures.

scholarly journals Edge Detection and Depth to Magnetic Source Estimation in Part of Central Nigeria

2020 ◽  
pp. 54-67
Author(s):  
Mam D. Tawey ◽  
Abbass A. Adetona ◽  
Usman D. Alhassan ◽  
Abdulwaheed A. Rafiu ◽  
Kazeem A. Salako ◽  
...  
Keyword(s):  
Edge Detection ◽  
Trend Analysis ◽  
Source Parameter ◽  
Fault System ◽  
Euler Deconvolution ◽  
Structural Index ◽  
The South ◽  
Source Estimation ◽  
Magnetic Source ◽  
3D Euler Deconvolution

High-resolution aeromagnetic data covering an area of 24, 200 km2 in north central Nigeria has been acquired and analyzed with the aim of carrying out trend analysis, edge detection (structural delineation) and depth to magnetic source estimation using reduce to the pole (RTP), horizontal gradient magnitude (HGM), center for exploration targeting plug-in (CET), 3D Euler deconvolution and source parameter imaging (SPI) techniques. Trend analysis was applied to the RTP data to delineate structures that have dissected the area. The 3D Euler deconvolution and HGM were correlated by plotting the estimated Euler solutions for a structural index of one (SI=1) on HGM map and the resulting map produced have shown that both methods can contribute in the interpretation of the general structural framework of the study area. The structural delineation based on HGM and CET maps showed that two predominant trends (ENE-WSW) and (WNW-ENE) have affected the area. The trend/depth/contacts of these faults were classified into four groups: Faults <150 m, 150 m - 300 m, 300m - 450 m which are the most predominant fault system based on Euler solutions with a structural index of one (SI=1) and those deeper than 450 m while the result of source parameter imaging (SPI) revealed a depth to source varying from 58 m specifically for areas with shallow depth to the magnetic source to those from deeper source occurring at 588.153m depth especially the south-central portion and the south-eastern portion of the study area.

Sign in / Sign up

Export Citation Format

Share Document