16 Jun 2021

16 Jun 2021

Review status: this preprint is currently under review for the journal SE.

Constraining 3D geometric gravity inversion with 2D reflection seismic profile using a generalized level-set approach: application to Eastern Yilgarn craton

Mahtab Rashidifard1,3, Jérémie Giraud1,3, Mark Lindsay1,3, Mark Jessell1,3, and Vitaliy Ogarko1,2 Mahtab Rashidifard et al.
  • 1Centre of Exploration Targeting (School of Earth Sciences), University of Western Australia, 35 Stirling Highway, Crawley, WA, 6009, Australia
  • 2International Centre for Radio Astronomy Research (ICRAR), University of Western Australia, 7 Fairway, Crawley, WA 6009, Australia
  • 3Mineral Exploration Cooperative Research Centre, School of Earth Sciences, University of Western Australia, 35 Stirling Highway, WA Crawley 6009, Australia

Abstract. One of the main tasks in 3D geological modelling is the boundary parametrization of the subsurface from geological observations and geophysical inversions. Several approaches have been developed for geometric inversion and joint inversion of geophysical datasets. However, the robust, quantitative integration of models and datasets with different spatial coverage, resolution, and levels of sparsity remains challenging. One promising approach for recovering the boundary of the geological units is the utilization of a level-set inversion method with potential field data. We focus on constraining 3D geometric gravity inversion with sparse lower-uncertainty information from a 2D seismic section.

We use a level-set approach to recover the geometry of geological bodies using two synthetic examples and data from the geologically complex Yamarna terrane (Yilgarn craton, Western Australia). In this study, a 2D seismic section has been used for constraining the location of rock unit boundaries being solved during the 3D gravity geometric inversion. The proposed work is the first we know of that automates the process of adding spatially distributed constraints to the 3D level-set inversion. In many hard-rock geoscientific investigations, seismic data is sparse and our results indicate that unit boundaries from gravity inversion can be much better constrained with seismic information even though they are sparsely distributed within the model. Thus, we conclude that it has the potential to bring the state of the art a step further towards building a 3D geological model incorporating several sources of information in similar regions of investigation.

Mahtab Rashidifard et al.

Status: open (until 06 Aug 2021)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Mahtab Rashidifard et al.

Data sets

Constrained 3D geometric gravity inversion Mahtab Rashidifard; Jérémie Giraud; Mark Lindsay; Mark Jessell; Vitaliy Ogarko

Mahtab Rashidifard et al.


Total article views: 365 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
293 68 4 365 2 4
  • HTML: 293
  • PDF: 68
  • XML: 4
  • Total: 365
  • BibTeX: 2
  • EndNote: 4
Views and downloads (calculated since 16 Jun 2021)
Cumulative views and downloads (calculated since 16 Jun 2021)

Viewed (geographical distribution)

Total article views: 328 (including HTML, PDF, and XML) Thereof 328 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 28 Jul 2021
Short summary
One motivation for this study is to develop a workflow that enables the integration of geophysical datasets with different coverage which are quite common in exploration geophysics. We have utilized a level-set approach to achieve this goal. The utilized technique parametrizes the subsurface in the same fashion as geological models. Our results indicate that the approach is capable of integrating the information from the seismic in 2D to guide the 3D inversion results of the gravity data.