scholarly journals Estimating the Extreme Values in Gridding Data using Gauss Elimination Method: A New Approach in Potential Field Processing

2021 ◽  
Vol 54 (2E) ◽  
pp. 150-163
Author(s):  
Nguyen Kim Dung
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Extreme Points ◽  
Quadratic Function ◽  
Real Data ◽  
Maximum Point ◽  
Geological Structures ◽  
Potential Field Data ◽  
Straight Line ◽  
Model Structures

The position of a maximum point of a function depends on its coefficients and order. The maximum horizontal gradient method is a popular method that greatly contributes to the detection of maximum points and approximation of geological structures edges. By adopting a mathematical logic, Blakely and Simpson established a quadratic function based on the characteristic of three points of a straight line in the fundamental directions. However, for potential field data like gravity and magnetic data, the coefficients of a quadratic function in each direction are not only dependent on the values of three points on a straight line, but also, they depend on the values of the surrounding points. This article proposes an algorithm which can detect maximum points more effectively in order to delineate geological structures boundaries from potential field data. The proposed algorithm uses a 3×3 neighborhood data grid to establish a two-variables function and to determine its coefficients by applying the Gaussian elimination method. After the two-variables function has been established, the algorithm estimates any extreme points and their positions from a set of four single-variable functions which correspond to the horizontal, vertical and the two diagonal directions by the cases x = 0, y = 0, y = -x and y = x of the main function. Finally, the conditions to detect the maximum point from the extreme points are defined. The validity of the algorithm was demonstrated on synthetic datasets generated by two different model structures. A real data application of the method has also been realized by estimating the geological boundaries by gravity data in the Vietnam’s continental shelf. The results obtained from the synthetic applications of the algorithm proved that it can determine more maximum points as compared to Blakely and Simpson method, and as a result, in all the test cases, it has drawn the real boundaries of the model structures more accurately. The application results of the method on real data revealed new boundary delineations in the research area, interpreted to be faults or fractures which lies between deep trench in the East Vietnam Sea.

Discrete linear transformations of potential field data

Geophysics ◽  
1989 ◽  
Vol 54 (4) ◽  
pp. 497-507 ◽  
Author(s):  
Jorge W. D. Leão ◽  
João B. C. Silva
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Amazon Basin ◽  
Real Data ◽  
Damping Factor ◽  
Linear Transformations ◽  
Potential Field Data ◽  
Data Window ◽  
Green’S Matrix ◽  
Equivalent Layer

We present a new approach to perform any linear transformation of gridded potential field data using the equivalent‐layer principle. It is particularly efficient for processing areas with a large amount of data. An N × N data window is inverted using an M × M equivalent layer, with M greater than N so that the equivalent sources extend beyond the data window. Only the transformed field at the center of the data window is computed by premultiplying the equivalent source matrix by the row of the Green’s matrix (associated with the desired transformation) corresponding to the center of the data window. Since the inversion and the multiplication by the Green’s matrix are independent of the data, they are performed beforehand and just once for given values of N, M, and the depth of the equivalent layer. As a result, a grid operator for the desired transformation is obtained which is applied to the data by a procedure similar to discrete convolution. The application of this procedure in reducing synthetic anomalies to the pole and computing magnetization intensity maps shows that grid operators with N = 7 and M = 15 are sufficient to process large areas containing several interfering sources. The use of a damping factor allows the computation of meaningful maps even for unstable transformations in the presence of noise. Also, an equivalent layer larger than the data window takes into account part of the interfering sources so that a smaller damping factor is employed as compared with other damped inversion methods. Transformations of real data from Xingú River Basin and Amazon Basin, Brazil, demonstrate the contribution of this procedure for improvement of a preliminary geologic interpretation with minimum a priori information.

A method for inversion of 1D vertical soundings of gravity anomalies

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. G15-G23
Author(s):  
Andrea Vitale ◽  
Domenico Di Massa ◽  
Maurizio Fedi ◽  
Giovanni Florio
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Gravity Data ◽  
Finite Size ◽  
Gravity Anomalies ◽  
Real Data ◽  
Gravity Inversion ◽  
Chain Process ◽  
Potential Field Data ◽  
Gamma Density

We have developed a method to interpret potential fields, which obtains 1D models by inverting vertical soundings of potential field data. The vertical soundings are built through upward continuation of potential field data, measured on either a profile or a surface. The method assumes a forward problem consisting of a volume partitioned in layers, each of them homogeneous and horizontally finite, but with the density changing versus depth. The continuation errors, increasing with the altitude, are automatically handled by determining the coefficients of a third-order polynomial function of the altitude. Due to the finite size of the source volume, we need a priori information about the total horizontal extent of the volume, which is estimated by boundary analysis and optimized by a Markov chain process. For each sounding, a 1D inverse problem is independently solved by a nonnegative least-squares algorithm. Merging of the several inverted models finally yields approximate 2D or 3D models that are, however, shown to generate a good fit to the measured data. The method is applied to synthetic models, producing good results for either perfect or continued data. Even for real data, i.e., the gravity data of a sedimentary basin in Nevada, the results are interesting, and they are consistent with previous interpretation, based on 3D gravity inversion constrained by two gamma-gamma density logs.

scholarly journals An edge-assisted smooth method for potential field data

2021 ◽  
Vol 18 (1) ◽  
pp. 113-123
Author(s):  
Shijing Zheng ◽  
Xiaohong Meng ◽  
Jun Wang
Keyword(s):  
Edge Detection ◽  
Potential Field ◽  
Field Data ◽  
Real Data ◽  
Linear Functions ◽  
Detection Methods ◽  
Edge Preserving ◽  
Potential Field Data ◽  
Tectonic Boundaries ◽  
Derivatives Of

Abstract Edge detection is one of the most commonly used methods for the interpretation of potential field data, because it can highlight the horizontal inhomogeneous of underground geological bodies (faults, tectonic boundaries, etc.). A variety of edge detection methods have been reported in the literature, most of which are based on the combined transformation results of horizontal and vertical derivatives of the observations. Consequently, these edge detection methods are sensitive to noise. Therefore, noise reduction is desirable ahead of applying edge detection methods. However, the application of conventional filters smears discontinuities in the data to a certain extent, which would inevitably induce unfavourable influence on subsequent edge detection. To solve this problem, a novel edge-preserving smooth method for potential field data is proposed, which is based on the concept of guided filter developed for image processing. The new method substitutes each data point by a combination of a series of coefficients of linear functions. It was tested on synthetic model and real data, and the results showed that it can effectively smooth potential field data while preserving major structural and stratigraphic discontinuities. The obtained data from the new filter contain more obvious features of existing faults, which brings advantageous to further geological interpretations.

Sparsity constraints using a Laplacian kernel to get geological structures from potential field data

2017 ◽  
Vol 16 (1) ◽  
pp. 39-52 ◽  
Author(s):  
Zhaohai Meng ◽  
Fengting Li ◽  
XueChun Xu ◽  
Wenyue Zhou ◽  
Danian Hang
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Geological Structures ◽  
Potential Field Data ◽  
Sparsity Constraints

An improved approach for detecting the locations of the maxima in interpreting potential field data

2021 ◽  
Vol 14 (1) ◽  
Author(s):  
Luan Thanh Pham ◽  
Ozkan Kafadar ◽  
Erdinc Oksum ◽  
Ahmed M. Eldosouky
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Potential Field Data

Adaptive sampling of potential-field data: A direct approach to compressive inversion

Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. IM1-IM9 ◽  
Author(s):  
Nathan Leon Foks ◽  
Richard Krahenbuhl ◽  
Yaoguo Li
Keyword(s):  
Potential Field ◽  
Adaptive Sampling ◽  
Field Data ◽  
Computational Cost ◽  
Large Field ◽  
Magnetic Data ◽  
Data Sampling ◽  
Data Types ◽  
Data Set ◽  
Potential Field Data

Compressive inversion uses computational algorithms that decrease the time and storage needs of a traditional inverse problem. Most compression approaches focus on the model domain, and very few, other than traditional downsampling focus on the data domain for potential-field applications. To further the compression in the data domain, a direct and practical approach to the adaptive downsampling of potential-field data for large inversion problems has been developed. The approach is formulated to significantly reduce the quantity of data in relatively smooth or quiet regions of the data set, while preserving the signal anomalies that contain the relevant target information. Two major benefits arise from this form of compressive inversion. First, because the approach compresses the problem in the data domain, it can be applied immediately without the addition of, or modification to, existing inversion software. Second, as most industry software use some form of model or sensitivity compression, the addition of this adaptive data sampling creates a complete compressive inversion methodology whereby the reduction of computational cost is achieved simultaneously in the model and data domains. We applied the method to a synthetic magnetic data set and two large field magnetic data sets; however, the method is also applicable to other data types. Our results showed that the relevant model information is maintained after inversion despite using 1%–5% of the data.

Sign in / Sign up

Export Citation Format

Share Document