A new approach to depth determination from magnetic anomalies

Geophysics ◽  
2002 ◽  
Vol 67 (5) ◽  
pp. 1524-1531 ◽  
Author(s):  
El‐Sayed M. Abdelrahman ◽  
Hesham M. El‐Araby ◽  
Tarek M. El‐Araby ◽  
Khalid S. Essa

We have developed a semiautomatic method to determine the depth to shallow and deep‐seated structures from a magnetic anomaly profile. It involves using a relationship between the depths to two coaxial sources obtained by combining observations at symmetric points with respect to the coordinate of the sources center. For five established, fixed data points, the depth to the shallow structure is determined for each preassigned depth of the deep‐seated structure. The computed depths to the shallow structure are plotted against the computed depths to the deep‐seated structure, yielding a continuous, monotonically increasing depth curve. The spacing between the observations is then modified, producing several curves. The accepted estimates for the depths to both structures are read at the common intersection of these curves. The effective intensity and the angle of magnetization of both structures are also estimated. The proposed method was tested both on noisy synthetic and real magnetic data. In the case of synthetic data, the depth curves method determined the correct depths of both coaxial and laterally offset sources. In the case of practical data (vertical component anomaly over a chromite body in the Guleman concession, Turkey), the method suggested the shape of the buried shallow structure resembles a horizontal cylinder model buried at a depth of 31 m and the shape of the buried deep seated structure resembles a dike model buried at a depth of 62 m. The estimated shape and depth of the shallow structure are in very good agreement with the results obtained from drilling and surface geology. The area appears to still hold promise for chromite exploration from the deeper structure.

2014 ◽  
Vol 644-650 ◽  
pp. 3793-3796
Author(s):  
Liang Hui Guo ◽  
Rui Gao ◽  
Guo Li Zhang

Under the effects of remanent magnetization, total magnetization direction is different from geomagnetic field direction, which makes magnetic data processing and interpretation complexity. In this paper, we present a new approach for estimating the total magnetization direction of sources via cross-correlation between the reduced-to-pole anomaly and the normalized source strength (who is less sensitive to remanent magnetization). The geomagnetic field direction is used to calculated the normalized source strength, while various assumed total magnetization directions are used to calculated the RTP anomalies. The maximum correlation between the RTP anomalies and the normalized corresponds to the estimated total magnetization direction. Test on synthetic data showed that the new approach is simple and effective.


Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 205-212 ◽  
Author(s):  
E. M. Abdelrahman ◽  
E. R. Abo‐Ezz

This paper presents a new approach for determining the depth of a buried structure from numerical second‐, third‐, and fourth‐horizontal‐derivative anomalies obtained from 2-D magnetic data using filters of successive graticule spacings. The problem of depth determination has been transformed into the problem of finding a solution to a nonlinear equation of the form z = f(z). Formulas have been derived for a horizontal cylinder and a dike. The depths obtained from the second‐, third‐, and fourth‐derivative anomaly values can be used to determine simultaneously the actual depth to the buried structure and the optimum order of the regional magnetic field along the profile. This powerful technique can solve two major potential field problems: regional residual separation and depth determination. The method is applied to theoretical data with and without random errors and is tested on a field example from Arizona.


Geophysics ◽  
2000 ◽  
Vol 65 (1) ◽  
pp. 126-131 ◽  
Author(s):  
El‐Sayed M. Abdelrahman ◽  
Hamdy I. Hassanein

We have developed a simple method to simultaneously determine the shape (shape factor) and the depth of a buried structure from magnetic data. The method is similar to Euler deconvolution, but it solves for shape and depth independently. The method involves using a relationship between the shape factor, the depth to the source, and a combination of observations at symmetric points with respect to the coordinate of the source center with a free parameter (graticule spacing). The relationship represents a parametric family of curves. For a fixed free parameter, the depth is determined for each shape factor. The computed depths are plotted against the shape factors representing a continuous monotonically increasing curve. The solution for the shape and depth of the buried structure is read at the common intersection of the parametric curves. The parametric curves method is applied to two theoretical magnetic anomaly profiles due to a thin dike and a horizontal cylinder. The method is also tested on a field example from Brazil. In all cases, the shape and depth solutions obtained are in good agreement with the actual ones.


Geofizika ◽  
2020 ◽  
Vol 37 (2) ◽  
pp. 213-235
Author(s):  
Ata Eshaghzadeh ◽  
Sanaz Seyedi Sahebari

This paper presents a nature-based algorithm, titled multivariable teaching-learning-based optimization (MTLBO) algorithm. MTLBO algorithm during an iterative process can estimates the best values of the buried structure (model) parameters in a multi-objective problem. The algorithm works in two computational phases: the teacher phase and the learner phase. The major purpose of the MTLBO algorithm is to modify the value of the learners and thus, improving the value of the model parameters which leads to the optimal solution. The variables of each learner (model) are the depth (z), amplitude coefficient (k), shape factor (q), angle of effective magnetization (θ) and axis location (x0) parameters. We employ MTLBO method for the magnetic anomalies caused by the buried structures with a simple geometric shape such as sphere and horizontal cylinder. The efficiency of the MTLBO is also studied by noise corruption synthetic data, as the acceptable results were obtained. We have applied the MTLBO for the interpretation of the four magnetic anomaly profiles from Iran, Brazil and India.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Manzar Fawad ◽  
Nazmul Haque Mondol

AbstractGeological CO2 storage can be employed to reduce greenhouse gas emissions to the atmosphere. Depleted oil and gas reservoirs, deep saline aquifers, and coal beds are considered to be viable subsurface CO2 storage options. Remote monitoring is essential for observing CO2 plume migration and potential leak detection during and after injection. Leak detection is probably the main risk, though overall monitoring for the plume boundaries and verification of stored volumes are also necessary. There are many effective remote CO2 monitoring techniques with various benefits and limitations. We suggest a new approach using a combination of repeated seismic and electromagnetic surveys to delineate CO2 plume and estimate the gas saturation in a saline reservoir during the lifetime of a storage site. This study deals with the CO2 plume delineation and saturation estimation using a combination of seismic and electromagnetic or controlled-source electromagnetic (EM/CSEM) synthetic data. We assumed two scenarios over a period of 40 years; Case 1 was modeled assuming both seismic and EM repeated surveys were acquired, whereas, in Case 2, repeated EM surveys were taken with only before injection (baseline) 3D seismic data available. Our results show that monitoring the CO2 plume in terms of extent and saturation is possible both by (i) using a repeated seismic and electromagnetic, and (ii) using a baseline seismic in combination with repeated electromagnetic data. Due to the nature of the seismic and EM techniques, spatial coverage from the reservoir's base to the surface makes it possible to detect the CO2 plume’s lateral and vertical migration. However, the CSEM low resolution and depth uncertainties are some limitations that need consideration. These results also have implications for monitoring oil production—especially with water flooding, hydrocarbon exploration, and freshwater aquifer identification.


Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. F239-F250 ◽  
Author(s):  
Fernando A. Monteiro Santos ◽  
Hesham M. El-Kaliouby

Joint or sequential inversion of direct current resistivity (DCR) and time-domain electromagnetic (TDEM) data commonly are performed for individual soundings assuming layered earth models. DCR and TDEM have different and complementary sensitivity to resistive and conductive structures, making them suitable methods for the application of joint inversion techniques. This potential joint inversion of DCR and TDEM methods has been used by several authors to reduce the ambiguities of the models calculated from each method separately. A new approach for joint inversion of these data sets, based on a laterally constrained algorithm, was found. The method was developed for the interpretation of soundings collected along a line over a 1D or 2D geology. The inversion algorithm was tested on two synthetic data sets, as well as on field data from Saudi Arabia. The results show that the algorithm is efficient and stable in producing quasi-2D models from DCR and TDEM data acquired in relatively complex environments.


2012 ◽  
Vol 2012 ◽  
pp. 1-10 ◽  
Author(s):  
Lev V. Utkin

A fuzzy classification model is studied in the paper. It is based on the contaminated (robust) model which produces fuzzy expected risk measures characterizing classification errors. Optimal classification parameters of the models are derived by minimizing the fuzzy expected risk. It is shown that an algorithm for computing the classification parameters is reduced to a set of standard support vector machine tasks with weighted data points. Experimental results with synthetic data illustrate the proposed fuzzy model.


Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. V41-V59 ◽  
Author(s):  
Olena Tiapkina ◽  
Martin Landrø ◽  
Yuriy Tyapkin ◽  
Brian Link

The advent of single receiver point, multi-component geophones has necessitated that ground roll be removed in the processing flow rather than through acquisition design. A wide class of processing methods for ground-roll elimination is polarization filtering. A number of these methods use singular value decomposition (SVD) or some related transformations. We focus on a single-station SVD-based polarization filter that we consider to be one of the best in the industry. The method is comprised of two stages: (1) ground-roll detection and (2) ground-roll estimation and filtering. To detect the ground roll, a special attribute dependent on the singular values of a three-column matrix formed by a sliding time window is used. The ground roll is approximated and subtracted using the first two eigenimages of this matrix. To limit the possible damage to the signal, the filter operates within the record intervals where the ground roll is detected and within the ground-roll frequency bandwidth only. We improve the ground-roll detector to make it theoretically insensitive to ambient noise and more sensitive to the presence of ground roll. The advantage of the new detector is demonstrated on synthetic and field data sets. We estimate theoretically and with synthetic data the attenuation of the underlying reflections that can be caused by the polarization filter. We show that the underlying signal always loses almost all the energy on the vertical component and on the horizontal component in the ground-roll propagation plane and within the ground-roll frequency bandwidth. The only signal component, if it exists, that can retain a significant part of its energy is the horizontal component orthogonal to the above plane. When 2D 3C field operations are conducted, the signal particle motion can deviate from the ground-roll propagation plane and can therefore retain some of its energy due to a set of offline reflections. In the case of 3D 3C seismic surveys, the reflected signal always deviates from the ground-roll propagation plane on the receiver lines that do not contain the source. This is confirmed with a 2.5D 3C synthetic data set. We discuss when the ability of the filter to effectively subtract the ground roll may, or may not, allow us to ignore the inevitable harm that is done to the underlying reflected waves.


2019 ◽  
Vol 486 (1) ◽  
pp. L105-L109 ◽  
Author(s):  
David W Pearson ◽  
Lado Samushia

ABSTRACT As we move towards future galaxy surveys, the three-point statistics will be increasingly leveraged to enhance the constraining power of the data on cosmological parameters. An essential part of the three-point function estimation is performing triplet counts of synthetic data points in random catalogues. Since triplet counting algorithms scale at best as $\mathcal {O}(N^2\log N)$ with the number of particles and the random catalogues are typically at least 50 times denser than the data; this tends to be by far the most time-consuming part of the measurements. Here, we present a simple method of computing the necessary triplet counts involving uniform random distributions through simple one-dimensional integrals. The method speeds up the computation of the three-point function by orders of magnitude, eliminating the need for random catalogues, with the simultaneous pair and triplet counting of the data points alone being sufficient.


2020 ◽  
Vol 110 (5) ◽  
pp. 2530-2540 ◽  
Author(s):  
Adam T. Ringler ◽  
Robert E. Anthony ◽  
David C. Wilson ◽  
Abram C. Claycomb ◽  
John Spritzer

ABSTRACT Seismometers are highly sensitive instruments to not only ground motion but also many other nonseismic noise sources (e.g., temperature, pressure, and magnetic field variations). We show that the Alaska component of the Transportable Array is particularly susceptible to recording magnetic storms and other space weather events because the sensors used in this network are unshielded and magnetic flux variations are stronger at higher latitudes. We also show that vertical-component seismic records across Alaska are directly recording magnetic field variations between 40 and 800 s period as opposed to actual ground motion during geomagnetic events with sensitivities ranging from 0.004 to 0.48  (m/s2)/T. These sensitivities were found on a day where the root mean square variation in the magnetic field was 225 nT. Using a method developed by Forbriger (2007, his section 3.1), we show that improving vertical seismic resolution of an unshielded sensor by as much as 10 dB in the 100–400 s period band using magnetic data from a collocated three-component magnetometer is possible. However, due to large spatial variations in Earth’s magnetic field, this methodology becomes increasingly ineffective as the distance between the seismometer and magnetometer increases (no more than 200 km separation). A potential solution to this issue may be to incorporate relatively low-cost magnetometers as an additional environmental data stream at high-latitude seismic stations. We demonstrate that the Bartington Mag-690 sensors currently deployed at Global Seismographic Network sites are not only acceptable for performing corrections to seismic data, but are also capable of recording many magnetic field signals with similar signal-to-noise ratios, in the 20–1000 s period band, as the observatory grade magnetometers operated by the U.S. Geological Survey Geomagnetism Program. This approach would densify magnetic field observations and could also contribute to space weather monitoring by supplementing highly calibrated magnetometers with additional sensors.


Sign in / Sign up

Export Citation Format

Share Document