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.

scholarly journals Measuring change at Earth’s surface: On-demand vertical and three- dimensional topographic differencing implemented in OpenTopography

Geosphere ◽  
2021 ◽  
Author(s):  
Chelsea Scott ◽  
Minh Phan ◽  
Viswanath Nandigam ◽  
Christopher Crosby ◽  
J Ramon Arrowsmith
Keyword(s):  
Three Dimensional ◽  
Point Clouds ◽  
Data Sets ◽  
Data Set ◽  
Landscape Characteristics ◽  
On Demand ◽  
Digital Elevation ◽  
Rigid Deformation ◽  
Optimal Resolution ◽  
Topography Data

Topographic differencing measures landscape change by comparing multitemporal high-resolution topography data sets. Here, we focused on two types of topographic differencing: (1) Vertical differencing is the subtraction of digital elevation models (DEMs) that span an event of interest. (2) Three-dimensional (3-D) differencing measures surface change by registering point clouds with a rigid deformation. We recently released topo­graphic differencing in OpenTopography where users perform on-demand vertical and 3-D differencing via an online interface. OpenTopography is a U.S. National Science Foundation–funded facility that provides access to topographic data and processing tools. While topographic differencing has been applied in numerous research studies, the lack of standardization, particularly of 3-D differencing, requires the customization of processing for individ­ual data sets and hinders the community’s ability to efficiently perform differencing on the growing archive of topography data. Our paper focuses on streamlined techniques with which to efficiently difference data sets with varying spatial resolution and sensor type (i.e., optical vs. light detection and ranging [lidar]) and over variable landscapes. To optimize on-demand differencing, we considered algorithm choice and displacement resolution. The optimal resolution is controlled by point density, landscape characteristics (e.g., leaf-on vs. leaf-off), and data set quality. We provide processing options derived from metadata that allow users to produce optimal high-quality results, while experienced users can fine tune the parameters to suit their needs. We anticipate that the differencing tool will expand access to this state-of-the-art technology, will be a valuable educational tool, and will serve as a template for differencing the growing number of multitemporal topography data sets.

scholarly journals IMPLEMENTATION OF SURVEY AND THREE-DIMENSIONAL MONITORING OF ARCHAEOLOGICAL EXCAVATIONS OF THE KHIRBAT AL-DUSAQ SITE, JORDAN

2019 ◽  
Vol XLII-2/W15 ◽  
pp. 41-47 ◽  
Author(s):  
E. Alby ◽  
E. Vigouroux ◽  
R. Elter
Keyword(s):  
3D Reconstruction ◽  
Three Dimensional ◽  
Archaeological Site ◽  
Archaeological Sites ◽  
Data Sets ◽  
Data Set ◽  
Entire Complex ◽  
Future Data ◽  
Acquisition Processes ◽  
Building Walls

<p><strong>Abstract.</strong> In this paper will be presented the use of photogrammetry integrated to the process of representation of an archaeological site. The Khirbat al-Dūsaq site, Jordan, is an architectural complex composed by three remaining buildings with different shapes and functions. The first one is a reception building name īwān. The second one is vaulted and its function has not been determined yet. The third is a bath with all the complexity that are required for such a function (multiple rooms and sequence of spaces). The site is being excavated and there remains unknown information archaeologists want to discover and represent. This project takes places after several years of collaboration on different other archaeological sites. During these different projects, methods of acquisition, processes and drawings at different places and stages have been developed and work methods that includes the use of photogrammetry are now integrated to the archaeological practices. There is now a need by archaeologists for ortho-photos to draw precise plans. The integration of photogrammetry into the practice of archaeology on site helps also to reduce the time consumption to survey and to represent excavation activity. The data sets obtained year after year can also be used as a support for 3D reconstruction. The 3D modelling stage begins by integrating the context represented here by 3D textured mesh produced during the process of ortho-photos. The integration of photogrammetry started in 2015 by acquiring pictures from bath building. This work had to be extended to the entire complex so that it has been decided to manage it, in a proper way. In 2016, a survey network has been implemented, and complete photogrammetric data set have been produced. At this time there was a photogrammetric survey reference for all the data sets of the site. Several years of survey means that the project has to adapt to its specific context. The site life during 11 months without archaeological preoccupations signifies that it is evolving, so that in 2017, ground points had disappeared. The possibility to geo-reference future data sets imposes to integrate targets on pictures from 2016 data set. The remaining building walls on site keep their shape enough to be integrated as constant structures over the years. At first it has been decided to integrate photogrammetry technic to the representation process of the Khirbat al-Dūsaq site. It has proved, by the precision and flexibility of processes that good quality representations could be produced and the 3D documentation could be used as a support of 3D reconstruction stage also. Photogrammetric documentation, as soon as it is properly managed over the years can thus be integrated in archaeologic practices and can help to reduce time consuming stages and propose other activity support as 3D reconstruction.</p>

Sign in / Sign up

Export Citation Format

Share Document