Marquardt inverse modeling of gravity anomalies due to anticlinal and synclinal structures based on parabolic density contrast. (A case study: Kerend region)

Authors

1 Ph. D. student. Dept., of Geology, Faculty of sciences, University of Isfahan, Isfahan

2 Dept., of Civil Engineering, Nabi Akram Higher education Institute, Tabriz

Abstract

In this paper, a method based on the 2D Marquardt inverse modeling of anticlinal and synclinal structures using the gravity field data is introduced and the density contrast of the difference depths of the earth is computed based on a parabolic density function (PDF). The normal and inverted isosceles triangular models are generally used to describe the geometries of these structures in analyzing gravity anomalies. In the absence of known geology, it may not be possible to isolate completely the gravity signature due to a geological structure from the regional gravity background. In this paper a computer program solves three different parameters of the structure, z1, z2 and i, in addition to estimating two coefficients of regional gravity anomaly. The modeling process begins with computing the theoretical gravity anomaly of an anticline or syncline prototype, defined by approximate shape parameters in each case that can be attained from known geology. The program reduces the error between the observed gravity anomaly and the estimated gravity anomaly by improving the initial parameters of the model during repetition until the misfit function falls below a predefined allowable error or the damping factor acquires a large value or the specified number of iterations is completed. The efficiency of the algorithm is illustrated with a set of synthetic gravity anomalies over an anticlinal and a synclinal structure both with and without regional background, further, the code is exemplified with the gravity data from Kerend region, Iran. The target of gravity survey studies is determination the limits and depth values of anticlinal structures as a probable hydrocarbon traps. The results show the top depth of the anticlinal about 2800 m, the bottom depth of the anticlinal about 4200 m and the bottom depth of the synclinal about 5600 m.

Keywords

References

آقانباتی، ع (1383) زمین­شناسی ایران، انتشارات سازمان زمین­شناسی و اکتشافات معدنی کشور.
افشار، الف.، سهیلی، م.، والی، ن.، نقشه زمین­شناسی 1:250000 کوه کورخود. سازمان زمین­شناسی و اکتشافات معدنی کشور.
Alessandrello, E., Bichara, M. and Lakshmanan, J (1983) Automatic three layer, three-dimensional deconvolution of the Pays De Bray anticline. Geophysical Prospecting, 31: 608–626.
Chai, Y. and Hinze, W. J (1988) Gravity inversion of an interface above which the density contrast varies exponentially with depth. Geophysics, 53: 837 -845.
Chakravarthi, V (2003) US patent, 6: 615, 139.
Chakravarthi, V. and Sundararajan, N (2004) Ridge regression algorithm for gravity inversion of fault structures with variable density. Geophysics, 69: 1394–1404.
Chakravarthi, V. and Sundararajan, N (2005) Gravity modeling of 21/2-D sedimentary basins—a case of variable density contrast. Computers & Geosciences, 31: 820–827.
Chakravarthi, V. and Sundararajan, N (2007) Marquardt optimization of gravity anomalies of anticlinal and synclinal structures with prescribed depth dependent density. Geophysical Prospecting, 55: 571–587.
Chakravarthi, V. and Sundararajan, N (2008) TODGINV—A code for optimization of gravity anomalies due to anticlinal and synclinal structures with parabolic density contrast. Computers & Geosciences, 34: 955–966
Eshaghzadeh, A. and Kalantary, R. A (2017) Anticlinal Structure Modeling with Feed Forward Neural Networks for Residual Gravity Anomaly Profil, 8th congress of the Balkan Geophysical Society, DOI: 10.3997/2214-4609.201414210, 2015.
Eshaghzadeh, A., Dehghanpour, A. and Kalantari, R. A (2019) Marquardt inverse modeling of the residual gravity anomalies due to simple geometric structures: A cast study of Chromite deposit. Contributions to Geophysics and Geodesy, 49: 153-180.
Eshaghzadeh, A., Sahebari S. S. and Kalantari, R. A (2020) Anticlinal structure modeling using Feed-forward neural network (FNN) inversion for 2D gravity field data. Accepted in journal of the earth and space physics.
Eshaghzadeh, A., Sahebari S. S. and Kalantari, R. A (2020) Determination of sheet-like geological structures parameters using Marquardt inversion of the magnetic data. Indian Journal of Geo Marine Sciences, 49: 450-457.
Heiland, C. A (1968) Geophysical Exploration, 2nd ed. Hafner Publishing Co., New York, 1013pp.
Jianghai, X. and Sprowl, D. R (1995) Moho depths in Kansas from gravity inversion assuming exponential density contrast. Computers & Geosciences, 21: 237–244.
Lyons, P. L (1956) Geophysical Case Histories, Society of Exploration Geophysicists, 237–518.
Marquardt, D. W (1963) An algorithm for least squares estimation of nonlinear parameters. Journal of the Society for Industrial and Applied Mathematics, 11: 431–441.
Rao, D. B (1990) Analysis of gravity anomalies of sedimentary basins by an asymmetrical trapezoidal model with quadratic density function. Geophysics, 55: 226–231.
Rao, B. S. R. and Murty, I. V. R (1978) Gravity and Magnetic Methods of Prospecting. Arnold-Heinemann Publishers, New Delhi, India, 390pp.  
Jianghai, X. and Sprowl, D. R (1995) Moho depths in Kansas from gravity inversion assuming exponential density contrast. Computers & Geosciences, 21: 237–244.
Zhang, J., Zhong, B., Zhou, X. and Dai, Y (2001) Gravity anomalies of 2-D bodies with variable density contrast. Geophysics, 66: 809–813.

Files

History

  • Receive Date: 20 November 2020
  • Revise Date: 20 December 2020
  • Accept Date: 20 January 2021

Share

How to cite

Statistics