Open AccessMatlab tool REGCONT2: effective source depth estimation by means of Tikhonov’s regularized downwards continuation of potential fields
Author(s) -
Roman Pašteka,
David Kušnirák,
Roland Karcol
Publication year - 2018
Publication title -
contributions to geophysics and geodesy
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.235
H-Index - 18
eISSN - 1338-0540
pISSN - 1335-2806
DOI - 10.2478/congeo-2018-0010
Subject(s) - continuation , tikhonov regularization , extrapolation , interpolation (computer graphics) , unexploded ordnance , regularization (linguistics) , interpretation (philosophy) , sampling (signal processing) , algorithm , geology , mathematics , computer science , mathematical analysis , inverse problem , artificial intelligence , filter (signal processing) , computer vision , motion (physics) , remote sensing , programming language
Transformation based on downward continuation of potential fields is an important tool in their interpretation – depths of shallowest important sources can be determined by means of stable downward continuation algorithms. We analyse here selected properties of one from these algorithms (based on Tikhonov’s regularization approach) from the scope of two most important discretization parameters – dimensions of the areal coverage of the interpreted field and the sampling interval size. Estimation of the source depth is based on the analysis of computed L P -norms for various continuation depths. A typical local minimum of these norms disappears at the source depth. We show on several synthetic bodies (sphere, horizontal cylinder, vertical rod) and also real-world data-sets (results from a magnetic survey for unexploded ordnance detection) that there is a need for relatively large surroundings around the interpreted anomalies. Beside of this also the sampling step plays its important role – grids with finer sampling steps give better interpretation results, when using this downward continuation method. From this point of view, this method is more suitable for the interpretation of objects in near surface and mining geophysics (anomalies from cavities, unexploded ordnance objects and ore bodies). Anomalies should be well developed and separable, and densely sampled. When this is not valid, several algorithms of interpolation and extrapolation (grid padding methods) can improve the interpretation properties of studied downward continuation method.