GRAVITY DATA INVERSION USING L1-NORM STABILIZER

Q: How can I download an article?

To download an article from SID, first log in to the site, search for the article title, and click on the 'Download Article' option.

Q: How can I download an ISI article?

To download an ISI article on SID, enter the keyword or article title in the search bar, view the relevant results, click on the desired article, and select the 'Download Article' option.

Q: How can I access the SID database?

To access the SID database, visit SID.ir, create an account, and log in to access scientific resources.

Q: Is downloading articles from SID free?

Some articles on SID are available for free, while others require payment. Details are specified on the article's page.

VATANKHAH SAEED

مرکز اطلاعات علمی Scientific Information Database (SID) - Trusted Source for Research and Academic Resources

Journal Paper

Paper Information

Journal: Journal of the Earth and Space Physics Year:2016 | Volume:42 | Issue:2 Page(s): 337-348

Download Full-Text

Persian Verion

View:

1,089

Download:

Cites:

Information Journal Paper

Title

GRAVITY DATA INVERSION USING L1-NORM STABILIZER

Author(s)

VATANKHAH SAEED | | Issue Writer Certificate

Keywords

GRAVIMETRYQ2

INVERSIONQ2

REGULARIZATIONQ1

L1Q1

NORMQ1

REGULARIZATION PARAMETERQ1

Abstract

In this paper the INVERSION of gravity data using L1–norm stabilizer is considered. The INVERSION is an important step in the interpretation of data. In gravity data INVERSION, the goal is to estimate density and geometry of the unknown subsurface model from a set of known observation measured on the surface. Commonly, rectangular prisms are used to model the subsurface under the survey area. The unknown density contrasts within each prism are the parameters which should be estimated. The INVERSION of gravity data is an example of underdetermined and ill-posed problem, i.e. the solution can be non-unique and unstable. Thus, in order to find an acceptable solution REGULARIZATION should be imposed. Solution is usually obtained by minimizing a global objective function consisting of two terms, data misfit and the REGULARIZATION term. Data misfit measures how well an obtained model can reproduce the observed data. Usually, it is assumed noise in gravity data is Gaussian, therefore a L2–norm measure of the error between observed and predicted data is well suited for data misfit. There are several choices for a stabilizer, depends on type of features one wants to see from inverted model. A typical choice is a L2 –norm of a low-order differential operator applied to the model, which also a priori information and depth weighting can be incorporated (Li and Oldenburg, 1996). In this case the objective function is quadratic, then minimization of the function results a linear system to be solved. However, the models recovered in this way are characterized by smooth feature which are not always consistent with the real geological structures. There are situations in which the sources are localized and separated by sharp, distinct interfaces. To deal with this problem, during last decades, researchers have proposed a few types of stabilizer. Last and Kubik (1983) presented a compactness criterion for gravity INVERSION that seeks to minimize the area (or volume in 3D) of the causative body. Portniaguine and Zhdanov (1999) based on this stabilizer, who named the minimum support (MS), developed the minimum gradient support (MGS) stabilizer. For both constraint, the REGULARIZATION term can be written as the weighted L2–type NORM of the model. Therefore, the problem of the minimization of the objective function can be treated same as conventional Tikhonov functional. The only difference is that a priori variable weighting matrix for model parameters incorporated in the REGULARIZATION term. Thus the Iteratively Reweighted Least Square (IRLS) algorithm is required to solve the problem. Other possibility for stabilizer is the minimization of the L1-norm of model or gradient of model, the latter indicates total variation REGULARIZATION. The L1–norm stabilizer allows occurrence of large elements in the inverted model among mostly small values. Therefore, it can be used to obtain sharp boundaries and blocky features. Although the L1–norm stabilizer has favorable properties, in reconstruction of sparse models, its numerical implementation in a minimization problem can be difficult because its derivatives with respect to an element is not defined at zero. To overcome this difficulty, in this paper, the L1–norm stabilizer is approximated by a reweighted L2 –norm term. The algorithm is extended to gravity inverse problem, which needs depth weighting and other priori information to be included in the objective function. For estimating the REGULARIZATION parameter, which balances between two terms of objective function, the Unbiased Predictive Risk Estimator (UPRE) method is used. The solution of the resulting objective functional is found using Generalized Singular Value Decomposition (GSVD), also provides for efficient determination of the REGULARIZATION parameter at each iteration. Simulation using synthetic data of a dipping dike demonstrates that the method is capable to reconstruct focused image, boundaries and slop of the reconstructed model are close to those of the original model. The method is applied on gravity data acquired over the Gotvand dam site, in the south-west of Iran. The results show rather good agreement with those obtained from the boreholes.

Cites

No record.

References

No record.

Cite

APA: Copy

VATANKHAH, SAEED. (2016). GRAVITY DATA INVERSION USING L1-NORM STABILIZER. JOURNAL OF THE EARTH AND SPACE PHYSICS, 42(2), 337-348. SID. https://sid.ir/paper/80610/en

Vancouver: Copy

VATANKHAH SAEED. GRAVITY DATA INVERSION USING L1-NORM STABILIZER. JOURNAL OF THE EARTH AND SPACE PHYSICS[Internet]. 2016;42(2):337-348. Available from: https://sid.ir/paper/80610/en

IEEE: Copy

SAEED VATANKHAH, “GRAVITY DATA INVERSION USING L1-NORM STABILIZER,” JOURNAL OF THE EARTH AND SPACE PHYSICS, vol. 42, no. 2, pp. 337–348, 2016, [Online]. Available: https://sid.ir/paper/80610/en

Related Journal Papers