Comparison of grid Euler deconvolution with and without 2D constraints using a realistic 3D magnetic basement model

Geophysics ◽  
2005 ◽  
Vol 70 (3) ◽  
pp. L13-L21 ◽  
Author(s):  
Simon E. Williams ◽  
J. Derek Fairhead ◽  
Guy Flanagan
Keyword(s):  
Three Dimensional ◽  
Euler Method ◽  
Magnetic Data ◽  
Data Sets ◽  
Euler Deconvolution ◽  
Data Set ◽  
Distinct Source ◽  
Magnetic Basement ◽  
Modeling Strategy ◽  
Source Geometry

We describe the application of a 2D-constrained grid Euler deconvolution method which is able to determine for each solution window whether the source structure is two dimensional, three dimensional, or poorly defined and to estimate the source location and depth. In each solution window, eigenvalues and eigenvectors are derived from the Euler equations and compared to threshold levels. A single eigenvalue below the given threshold and lying in the x–y-plane is shown to indicate a 2D source, while the absence of such an eigenvalue indicates a 3D source geometry. Two small eigenvalues indicate the field in the window has no distinct source. Applying these criteria to each solution window allows us to generate a map of source-geometry distribution. We evaluate the effectiveness of 2D-constrained grid Euler deconvolution using synthetic magnetic data generated from a 3D basement model based on real topography from an area with surface-exposed faulting. This modeling strategy provides a complex, nonidealized data set that compares Euler depth estimates directly to the known basement surface depth. Our results indicate that noninteger structural indices can be the most appropriate choice for some data sets, and the 2D-constrained grid Euler method images magnetic basement structure more clearly and unambiguously than the conventional grid Euler method.

scholarly journals An improved space domain algorithm for determining the 3-D structure of the magnetic basement

2019 ◽  
Vol 41 (1) ◽  
pp. 69-80
Author(s):  
Nguyen Thi Thu Hang ◽  
Erdinc Oksum ◽  
Le Huy Minh ◽  
Do Duc Thanh
Keyword(s):  
Magnetic Anomaly ◽  
Science And Technology ◽  
Three Dimensional ◽  
Magnetic Anomalies ◽  
Magnetic Data ◽  
Marine Science ◽  
Euler Deconvolution ◽  
The Earth ◽  
Magnetic Basement ◽  
Basement Structures

The paper presents an improved algorithm based on Bhaskara Rao and Ramesh Babu’s algorithm to invert magnetic anomalies of three-dimensional basement structures. The magnetic basement is approximated by an ensemble of juxtaposed vertical prisms whose bottom surface coincides with Curie surface with the known depth. The computer program operating with the proposed algorithm is built in Matlab environment. Test applications show that the proposed method can perform computations with fast and stable convergence rate where the results also coincide well with the actual model structure. The effectiveness of the method is demonstrated by inverting magnetic anomalies of the southeast part of Vietnam continental shelf. The calculated magnetic basement relief of the study area provides useful additional information for studies in the aim of dealing with the geological structure of the area.References Beiki M., 2010. Analytic signals of gravity gradient tensor and their application to estimate source location, Geophysics, 75(6), i59–i74.Bui C.Q. (chief author), Le T., Tran T. D., Nguyen T. H., Phi T.T., 2007. Map of deep structure of the Earth’s crust, Atlas of the characteristics of natural conditions and environment in Vietnam’s waters and adjacent region. Publisher of Science and Technology, Ha Noi. Do D.T., Nguyen T.T.H., 2011. Atempt the improvement of inversion of magnetic anomalies of two dimensional polygonal cross sections to determine the depth of magnetic basement in some data profile of middle off shelf of Vietnam. Journal of Science and Technology, Vietnam Academy of Science and Technology, 49(2), 125–132.Do D.T., 2013. Study for application of 3D magnetic and gravity method to determine density contribution of basement rock and depth of magnetic basement on Vietnam’s shelf for oil research and prospecting Vietnam National University, Hanoi, Project code QG-11-04. Keating P. and Pilkington M., 2000, Euler deconvolution of the analytic signal, 62nd Annual International Meeting, EAGE, Session P0193.Keating P., Zerbo L., 1996. An improved technique for reduction to the pole at low latitudes, Geophysics, 61, 131–137.Le H.M., Luu V.H., 2003. Preliminary interpretation of the magnetic anomalies of the Eastern Vietnam sea and adiacent regions. J.  Sci. of the Earth, 25(2), 173–181. Mai T.T., Pham V.T., Dang V.B., Le D.B., Nguyen B., Le V.D., 2011. Characteristics of Pliocene - Quaternary geology and Geoengineering in the Center and Southeast parts of Continental Shelf of Vietnam. J.  Sci.  of the Earth, 33(2), 109-118.Mushayandebvu M.F., Lesur V., Reid A.B., Fairhead J.D., 2004. Grid Euler deconvolution with constraints for 2D structures, Geophysics, 69, 489–496.Nguyen N.T., Bui V.N., Nguyen T.T.H., Than D.L., 2014a. Application of power density spectrum of magnetic anomaly to estimate the structure of magnetic layer of the earth crust in the Bac Bo gulf. Journal of Marine Science and Technology, 14(4A), 137–148.Nguyen N.T., Bui V.N., Nguyen T.T.H., 2014b. Determining the depth to the magnetic basementand fault systems in Tu Chinh - Vung May area  by magnetic data interpretation. Journal of Marine Science and Technology, 14(4A), 16–25.Nguyen T.T.H., Pham T.L., Do D.T., Le H.M., 2018. Improving algorithm of determining the coordinates of the vertices of the polygon to invert magnetic anomalies of two-dimensional basement structures in space domain, Journal of Marine Science and Technology (preparing to print).Parker R.L., 1973. The rapid calculation of potential anomalies, Geophys. J. Roy. Astron. Soc, 31, 447–455. Pilkington M., Gregotski M.E., Todoeschuck J.P., 1994. Using fractal crustal magnetization models in magnetic interpretation, Geophysical Prospecting, 42, 677–692.Pilkington M., 2006. Joint inversion of gravity and magnetic data for two-layer models, Geophysics, 71, L35–L42.Rao D.B., Babu N.R., 1993. A fortran 77 computer program for three dimensional inversion of magnetic anomalies resulting from multiple prismatic bodies, Computer & Geosciences, 19(8), 781–801.Tanaka A., Okubo Y., Matsubayashi O., 1999. Curie point depth based on spectrum analysis of the magnetic anomaly data in East and Southeast Asia, Tectonic Pphysics, 306, 461–470.Thompson D.T., 1982. EULDTH – A new technique for marking computer-assisted depth estimates from magnetic data, Geophysics, 47, 31–37.Vo T.S., Le H.M., Luu V.H., 2005. Determining the horizontal position and depth of the density discontinuties in Red River Delta by using the vertical derivative and Euler deconvolution for the gravity anomaly data, Vietnam. Journal of Geology, Series A, 287(3–4), 39–52.  Werner S., 1955. Interpretation of magnetic anomalies of sheet-like bodies, Sveriges Geologiska Undersokning, Series C, Arsbok, 43, 6.Xu S.Z., 2006. The integral-iteration method for continuation of potential fields, Chinese journal of geophysics (in Chinese), 49(4), 1176–1182.Zhang C., Huang D.N., Zhang K., Pu Y.T., Yu P., 2016. Magnetic interface forward and inversion method based on Padé approximation, Applied Geophysics, 13(4), 712–720.CCOP, 1996. Magnetic anomaly map of East Asia, scale 1:4.000.000, Geological survey of Japan and Committee for co-ordination of joint prospecting for mineral resources in asian offshore areas.

scholarly journals A fast methodology for large-scale focusing inversion of gravity and magnetic data using the structured model matrix and the 2-D fast Fourier transform

2020 ◽  
Vol 223 (2) ◽  
pp. 1378-1397
Author(s):  
Rosemary A Renaut ◽  
Jarom D Hogue ◽  
Saeed Vatankhah ◽  
Shuang Liu
Keyword(s):  
Fourier Transform ◽  
Fast Fourier Transform ◽  
Linear Systems ◽  
Large Scale ◽  
Surface Measurement ◽  
Magnetic Data ◽  
Uniform Grid ◽  
Data Sets ◽  
Inversion Algorithm ◽  
Data Set

SUMMARY We discuss the focusing inversion of potential field data for the recovery of sparse subsurface structures from surface measurement data on a uniform grid. For the uniform grid, the model sensitivity matrices have a block Toeplitz Toeplitz block structure for each block of columns related to a fixed depth layer of the subsurface. Then, all forward operations with the sensitivity matrix, or its transpose, are performed using the 2-D fast Fourier transform. Simulations are provided to show that the implementation of the focusing inversion algorithm using the fast Fourier transform is efficient, and that the algorithm can be realized on standard desktop computers with sufficient memory for storage of volumes up to size n ≈ 106. The linear systems of equations arising in the focusing inversion algorithm are solved using either Golub–Kahan bidiagonalization or randomized singular value decomposition algorithms. These two algorithms are contrasted for their efficiency when used to solve large-scale problems with respect to the sizes of the projected subspaces adopted for the solutions of the linear systems. The results confirm earlier studies that the randomized algorithms are to be preferred for the inversion of gravity data, and for data sets of size m it is sufficient to use projected spaces of size approximately m/8. For the inversion of magnetic data sets, we show that it is more efficient to use the Golub–Kahan bidiagonalization, and that it is again sufficient to use projected spaces of size approximately m/8. Simulations support the presented conclusions and are verified for the inversion of a magnetic data set obtained over the Wuskwatim Lake region in Manitoba, Canada.

Rapid multi-scale analysis of near-surface geophysical anomaly maps: Application to an archaeo-geophysical data set

Geophysics ◽  
2020 ◽  
pp. 1-41 ◽  
Author(s):  
Jens Tronicke ◽  
Niklas Allroggen ◽  
Felix Biermann ◽  
Florian Fanselow ◽  
Julien Guillemoteau ◽  
...  
Keyword(s):  
Random Noise ◽  
Geophysical Data ◽  
Magnetic Data ◽  
Engineering Practice ◽  
Data Sets ◽  
Stripe Pattern ◽  
Data Set ◽  
Near Surface ◽  
Multi Scale ◽  
Topographic Data

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.

Limitations on digital filtering of the DNAG magnetic data set for the conterminous U.S.

Geophysics ◽  
1993 ◽  
Vol 58 (9) ◽  
pp. 1281-1296 ◽  
Author(s):  
V. J. S. Grauch
Keyword(s):  
Spatial Frequency ◽  
North American ◽  
Digital Filtering ◽  
Magnetic Data ◽  
Data Sets ◽  
Frequency Content ◽  
Long Distance ◽  
High Quality ◽  
Data Set ◽  
The U.S

The magnetic data set compiled for the Decade of North American Geology (DNAG) project presents an important digital data base that can be used to examine the North American crust. The data represent a patchwork from many individual airborne and marine magnetic surveys. However, the portion of data for the conterminous U.S. has problems that limit the resolution and use of the data. Now that the data are available in digital form, it is important to describe the data limitations more specifically than before. The primary problem is caused by datum shifts between individual survey boundaries. In the western U.S., the DNAG data are generally shifted less than 100 nT. In the eastern U.S., the DNAG data may be shifted by as much as 300 nT and contain regionally shifted areas with wavelengths on the order of 800 to 1400 km. The worst case is the artificial low centered over Kentucky and Tennessee produced by a series of datum shifts. A second significant problem is lack of anomaly resolution that arises primarily from using survey data that is too widely spaced compared to the flight heights above magnetic sources. Unfortunately, these are the only data available for much of the U.S. Another problem is produced by the lack of common observation surface between individual pieces of the U.S. DNAG data. The height disparities introduce variations in spatial frequency content that are unrelated to the magnetization of rocks. The spectral effects of datum shifts and the variation of spatial frequency content due to height disparities were estimated for the DNAG data for the conterminous U.S. As a general guideline for digital filtering, the most reliable features in the U.S. DNAG data have wavelengths roughly between 170 and 500 km, or anomaly half‐widths between 85 and 250 km. High‐quality, large‐region magnetic data sets have become increasingly important to meet exploration and scientific objectives. The acquisition of a new national magnetic data set with higher quality at a greater range of wavelengths is clearly in order. The best approach is to refly much of the U.S. with common specifications and reduction procedures. At the very least, magnetic data sets should be remerged digitally using available or newly flown long‐distance flight‐line data to adjust survey levels. In any case, national coordination is required to produce a consistent, high‐quality national magnetic map.

Connecting onshore-offshore Campos Basin structures: Interpretation of high-resolution airborne magnetic data

2014 ◽  
Vol 2 (4) ◽  
pp. SJ35-SJ45 ◽  
Author(s):  
Juarez Lourenço ◽  
Paulo T. L. Menezes ◽  
Valeria C. F. Barbosa
Keyword(s):  
Source Rocks ◽  
Magnetic Data ◽  
Data Sets ◽  
Hydrocarbon Migration ◽  
Data Set ◽  
Campos Basin ◽  
Turbidite Sedimentation ◽  
Edge Detection Method ◽  
Transfer Faults ◽  
Airborne Magnetic

We interpreted northwest-trending transfer faults whose extensions are not entirely mapped in the Precambrian basement of the onshore and offshore Campos Basin. To enhance the subtle northwest–southeast lineaments not clearly seen in the total-field data, we reprocessed and merged two airborne magnetic data sets aiming at producing a single merged magnetic data set. Next, we applied a directional filter to these integrated magnetic data. Finally, we applied a multiscale edge detection method to these filtered data. This combination allowed the detection of edges and ridges that are used to produce several northwest–southeast lineations. We interpreted these northwest-trending lineations as magnetic expressions of transfer faults that cut across the onshore adjacent basement of the Campos Basin to the shallow and deep Campos Basin waters. These interpreted northwest-trending faults suggested the continuity of the known northwest-striking transfer faults in the deep Campos Basin waters toward the shallow Campos Basin waters and the adjacent continent. Moreover, our interpreted northwest-trending faults revealed the control of several known oilfields in the Campos Basin. This result supported the hypothesis of the influence of the northwest–southeast-trending transfer faults on the petroleum system of Campos Basin, which were reactivated in the Tertiary providing a pathway for the turbidite sedimentation, reworking, and redistribution of several deepwater reservoirs. In addition, it was hypothesized that this faulting system controlled the hydrocarbon migration paths from the presalt source rocks through salt windows into basal suprasalt layers.

scholarly journals Accuracy of Low-Cost Alternative Facial Scanners: A Prospective Cohort Study.

2021 ◽  
Author(s):  
Alexander K. Bartella ◽  
Josefine Laser ◽  
Mohammad Kamal ◽  
Dirk Halama ◽  
Michael Neuhaus ◽  
...  
Keyword(s):  
Low Cost ◽  
Three Dimensional ◽  
Neck Surgery ◽  
Clinical Analysis ◽  
Smart Device ◽  
Data Sets ◽  
Average Deviation ◽  
Technical Improvement ◽  
Data Set ◽  
Significant Difference

Abstract Introduction: Three-dimensional facial scan images have been showing an increasingly important role in peri-therapeutic management of oral and maxillofacial and head and neck surgery cases. Face scan images can be open using optical facial scanners utilizing line-laser, stereophotography, structured light modality, or from volumetric data obtained from cone beam computed tomography (CBCT). The aim of this study is to evaluate, if two low-cost procedures for creating a three-dimensional face scan images are able to produce a sufficient data set for clinical analysis. Materials and methods: 50 healthy volunteers were included in the study. Two test objects with defined dimensions were attached to the forehead and the left cheek. Anthropometric values were first measured manually, and consecutively, face scans were performed with a smart device and manual photogrammetry and compared to the manually measured data sets.Results: Anthropometric distances on average deviated 2.17 mm from the manual measurement (smart device scanning 3.01 mm vs. photogrammetry 1.34 mm), with 7 out of 8 deviations were statistically significant. Of a total of 32 angles, 19 values showed a significant difference to the original 90° angles. The average deviation was 6.5° (smart device scanning 10.1° vs. photogrammetry 2.8°).Conclusion: Manual photogrammetry with a regular photo-camera shows higher accuracy than scanning with smart device. However, the smart device was more intuitive in handling and further technical improvement of the cameras used should be watched carefully.

scholarly journals Petrophysically and geologically guided multi-physics inversion using a dynamic Gaussian mixture model

2020 ◽  
Vol 224 (1) ◽  
pp. 40-68 ◽  
Author(s):  
Thibaut Astic ◽  
Lindsey J Heagy ◽  
Douglas W Oldenburg
Keyword(s):  
Gaussian Mixture Model ◽  
Mixture Model ◽  
Gaussian Mixture ◽  
Quantitative Information ◽  
Geophysical Data ◽  
Magnetic Data ◽  
Data Sets ◽  
Data Set ◽  
Geophysical Surveys ◽  
Geological Information

SUMMARY In a previous paper, we introduced a framework for carrying out petrophysically and geologically guided geophysical inversions. In that framework, petrophysical and geological information is modelled with a Gaussian mixture model (GMM). In the inversion, the GMM serves as a prior for the geophysical model. The formulation and applications were confined to problems in which a single physical property model was sought, and a single geophysical data set was available. In this paper, we extend that framework to jointly invert multiple geophysical data sets that depend on multiple physical properties. The petrophysical and geological information is used to couple geophysical surveys that, otherwise, rely on independent physics. This requires advancements in two areas. First, an extension from a univariate to a multivariate analysis of the petrophysical data, and their inclusion within the inverse problem, is necessary. Secondly, we address the practical issues of simultaneously inverting data from multiple surveys and finding a solution that acceptably reproduces each one, along with the petrophysical and geological information. To illustrate the efficacy of our approach and the advantages of carrying out multi-physics inversions coupled with petrophysical and geological information, we invert synthetic gravity and magnetic data associated with a kimberlite deposit. The kimberlite pipe contains two distinct facies embedded in a host rock. Inverting the data sets individually, even with petrophysical information, leads to a binary geological model: background or undetermined kimberlite. A multi-physics inversion, with petrophysical information, differentiates between the two main kimberlite facies of the pipe. Through this example, we also highlight the capabilities of our framework to work with interpretive geological assumptions when minimal quantitative information is available. In those cases, the dynamic updates of the GMM allow us to perform multi-physics inversions by learning a petrophysical model.

Semi-Automated Roller Parameters Extraction from Terrestrial Lidar

2021 ◽  
Vol 87 (12) ◽  
pp. 879-890
Author(s):  
Sagar S. Deshpande ◽  
Mike Falk ◽  
Nathan Plooster
Keyword(s):  
Rotation Axis ◽  
Moving Average ◽  
Three Dimensional ◽  
Data Sets ◽  
Steel Mill ◽  
Data Set ◽  
Terrestrial Lidar ◽  
Two Dimensional Image ◽  
Average Rotation

Rollers are an integral part of a hot-rolling steel mill. They transport hot metal from one end of the mill to another. The quality of the steel highly depends on the surface quality of the rollers. This paper presents semi-automated methodologies to extract roller parameters from terrestrial lidar points. The procedure was divided into two steps. First, the three-dimensional points were converted to a two-dimensional image to detect the extents of the rollers using fast Fourier transform image matching. Lidar points for every roller were iteratively fitted to a circle. The radius and center of the fitted circle were considered as the average radius and average rotation axis of the roller, respectively. These parameters were also extracted manually and were compared to the measured parameters for accuracy analysis. The proposed methodology was able to extract roller parameters at millimeter level. Erroneously identified rollers were identified by moving average filters. In the second step, roller parameters were determined using the filtered roller points. Two data sets were used to validate the proposed methodologies. In the first data set, 366 out of 372 rollers (97.3%) were identified and modeled. The second, smaller data set consisted of 18 rollers which were identified and modelled accurately.

Geostatistical Joint Interpretation of Gravity and Magnetotelluric Data of the US Cordillera

2021 ◽  
Author(s):  
Mohammad Shehata ◽  
Hideki Mizunaga
Keyword(s):  
Gravity Data ◽  
Three Dimensional ◽  
Data Sets ◽  
Spatial Expansion ◽  
Magnetotelluric Data ◽  
Data Set ◽  
Entire Area ◽  
Mean Values ◽  
Local Means ◽  
The Us

<p>Long-period magnetotelluric and gravity data were acquired to investigate the US cordillera's crustal structure. The magnetotelluric data are being acquired across the continental USA on a quasi-regular grid of ∼70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. International Gravimetreique Bureau compiled gravity Data at high spatial resolution. Due to the difference in data coverage density, the geostatistical joint integration was utilized to map the subsurface structures with adequate resolution. First, a three-dimensional inversion of each data set was applied separately.</p><p>The inversion results of both data sets show a similarity of structure for data structuralizing. The individual result of both data sets is resampled at the same locations using the kriging method by considering each inversion model to estimate the coefficient. Then, the Layer Density Correction (LDC) process's enhanced density distribution was applied to MT data's spatial expansion process. Simple Kriging with varying Local Means (SKLM) was applied to the residual analysis and integration. For this purpose, the varying local means of the resistivity were estimated using the corrected gravity data by the Non-Linear Indicator Transform (NLIT), taking into account the spatial correlation. After that, the spatial expansion analysis of MT data obtained sparsely was attempted using the estimated local mean values and SKLM method at the sections where the MT survey was carried out and for the entire area where density distributions exist. This research presents the integration results and the stand-alone inversion results of three-dimensional gravity and magnetotelluric data.</p>

Sign in / Sign up

Export Citation Format

Share Document