Open AccessEstimating the Extreme Values in Gridding Data using Gauss Elimination Method: A New Approach in Potential Field Processing
Author(s) -
Nguyen Kim Dung,
David Gómez-Ortíz,
Zaidoon Taha Abdulrazzaq,
Saada Ahmed Saada,
Tran Tuan Dung
Publication year - 2021
Publication title -
iraqi geological journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.2
H-Index - 5
eISSN - 2663-8754
pISSN - 2414-6064
DOI - 10.46717/igj.54.2e.10ms-2021-11-26
Subject(s) - diagonal , function (biology) , mathematics , extreme point , field (mathematics) , algorithm , data point , gaussian , quadratic equation , position (finance) , quadratic function , point (geometry) , line (geometry) , gauss , gaussian function , quadratic gauss sum , variable (mathematics) , mathematical analysis , geometry , combinatorics , physics , finance , quantum mechanics , evolutionary biology , pure mathematics , economics , biology
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.