Estimation of regularization parameter by active constraint balancing for 2D inversion of gravity data

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.

MOGHADASI M.; NEJATI KALATEH A.; REZAIE M.

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

Journal Paper

Paper Information

Journal: Journal of the Earth and Space Physics Year:2018 | Volume:44 | Issue:3 Page(s): 575-583

Download Full-Text

Persian Verion

View:

749

Download:

Cites:

Information Journal Paper

Title

Estimation of regularization parameter by active constraint balancing for 2D inversion of gravity data

Author(s)

MOGHADASI M. | NEJATI KALATEH A. | REZAIE M. | Issue Writer Certificate

Keywords

InversionQ2

Gravity dataQ2

Regularization parameterQ1

ACB methodQ1

Abstract

Inversion method is very common in the interpretation of practical Gravity data. The goal of 3D Inversion is to estimate density distribution of an unknown subsurface model from a set of known gravity observations measured on the surface. The Regularization parameter is one of the effective parameters for obtaining optimal model in Inversion of the Gravity data for similar Inversion of other geophysical data. For estimation of the optimum Regularization parameter the statistical criterion of Akaike’ s Bayesian Information Criterion (ABIC) usually used. This parameter is experimentally estimated in most Inversion methods. The choice of the Regularization parameter, which balances the minimization of the data misfit and model roughness, may be a critical procedure to achieve both resolution and stability. In this paper the Active Constraint Balancing (ACB) as a new method is used for estimating the Regularization parameter in two-dimensional (2-D) Inversion of Gravity data. This technique is supported by smoothness-constrained least-squares Inversion. We call this procedure “ active constraint balancing” (ACB). Introducing the Lagrangian multiplier as a spatially-dependent variable in the regularization term, we can balance the regularizations used in the Inversion. Spatially varying Lagrangian multipliers (Regularization parameters) are obtained by a parameter resolution matrix and Backus-Gilbert spread function analysis. For estimation of Regularization parameter by ACB method use must computed the resolution matrix R. The parameter resolution matrix R can be obtained in the Inversion process with pseudo-inverse G􀬾 multiplied by the kernel G. 􀜴 = 􀜩 􀬾 􀜩 (1) The spread function, which accounts for the inherent degree of how much the ith model parameter is not resolvable, defined as: 􀜵 􀜲 􀯜 = Σ 􀵫 􀝓 􀯜 􀯝 􀵫 1 − 􀝏 􀯜 􀯝 􀵯 􀜴 􀯜 􀯝 􀯇 􀵯 􀬶 􀯝 􀭀 􀬵 (2) where M is the total number of Inversion parameters, 􀝓 􀯜 􀯝 is a weighting factor defined by the spatial distance between the ith and jth model parameters, and 􀝏 􀯜 􀯝 is a factor which accounts for whether the constraint or regularization is imposed on the ith parameter and its neighboring parameters. In other words, the spread function defined here is the sum of the squared spatially weighted spread of the ith model parameter with respect to all of the model parameters excluding ones upon which a smoothness constraint is imposed. In this approach, the Regularization parameter λ (x, z) is set by a value from log-linear interpolation. log􁈺 λ 􀭧 􁈻 = log􁈺 λ 􀭫 􀭧 􀭬 􁈻 + 􀭪 􀭭 􀭥 􁈺 􀮛 􀱣 􀱗 􀱮 􁈻 􀬿 􀭪 􀭭 􀭥 􁈺 􀮛 􀱣 􀱟 􀱤 􁈻 􀭪 􀭭 􀭥 􁈺 􀭗 􀭔 􀱣 􀱗 􀱮 􁈻 􀬿 􀭪 􀭭 􀭥 􁈺 􀭗 􀭔 􀱣 􀱟 􀱤 􁈻 × 􁈼 log􁈺 SP􀭧 􁈻 − log􁈺 SP􀭫 􀭧 􀭬 􁈻 􁈽 (3) where 􀜵 􀜲 􀯠 􀯜 􀯡 and 􀜵 􀜲 􀯠 􀯔 􀯫 are the minimum and maximum values of spread function 􀜵 􀜲 􀯜 , respectively, and the 􀟣 􀯠 􀯜 􀯡 and 􀟣 􀯠 􀯔 􀯫 are minimum and maximum values of the Regularization parameter λ (x, z), which must be provided by the user. With this method, we can automatically set a smaller value λ (x, z) of the Regularization parameter to the highly resolvable model parameter, which corresponds to a smaller value of the spread function SP􀭧 in the Inversion process and vice versa. Users can choose these minimum and maximum Regularization parameters by setting variables LambdaMin and LambdaMax. For getting the target an algorithm is developed that estimates this parameter. The validity of the proposed algorithm has been evaluated by Gravity data acquired from a synthetic model. Then the algorithm used for Inversion of real Gravity data from Matanzas Cr deposit. The result obtained from 2D Inversion of Gravity data from this mine shows that this algorithm can provide good estimates of density anomalous structures within the subsurface.

Cites

No record.

References

No record.

Cite

APA: Copy

MOGHADASI, M., NEJATI KALATEH, A., & REZAIE, M.. (2018). Estimation of regularization parameter by active constraint balancing for 2D inversion of gravity data. JOURNAL OF THE EARTH AND SPACE PHYSICS, 44(3 ), 575-583. SID. https://sid.ir/paper/80406/en

Vancouver: Copy

MOGHADASI M., NEJATI KALATEH A., REZAIE M.. Estimation of regularization parameter by active constraint balancing for 2D inversion of gravity data. JOURNAL OF THE EARTH AND SPACE PHYSICS[Internet]. 2018;44(3 ):575-583. Available from: https://sid.ir/paper/80406/en

IEEE: Copy

M. MOGHADASI, A. NEJATI KALATEH, and M. REZAIE, “Estimation of regularization parameter by active constraint balancing for 2D inversion of gravity data,” JOURNAL OF THE EARTH AND SPACE PHYSICS, vol. 44, no. 3 , pp. 575–583, 2018, [Online]. Available: https://sid.ir/paper/80406/en

Related Journal Papers