Constrained simultaneous recovery of the depth to basement and lateral susceptibility variation

The basin environment is a widely studied subject in both geology and geophysics for its economic significance in energy and mineral explorations. However, the estimation of the basement depth is often a challenging task given the complexity of the basement relief and lateral physical property change. Previous works simplify the problem by only inverting for the depth to the basement, and more recent studies have suggested the need to incorporate the variation of physical properties to improve basement structure imaging. In this study, we develop an inversion method with the associated workflow to simultaneously recover both the depth to a magnetic basement and a laterally varying magnetic susceptibility in the basement rock. To achieve this, we employ a set of constraints on the inverse problem. Particularly, both the recovered susceptibility and basement depth models are bounded below a possible maximum value, and the depth model is guided by a few depth points obtained from the resistivity models that are obtained from the one-dimensional blocky inversions of magnetotelluric (MT) data. In addition, we apply the fuzzy C-means (FCM) clustering to the susceptibility model during the inversion and use the inverted cluster centers to differentiate for different geological units in the basement. To show the effectiveness of our work, we compare the existing approaches and our method using two test inversions on one synthetic model resembling the basin-basement environment before demonstrating our method on a field data example with magnetic data collected by the U.S. Geological Survey (USGS) over the Illinois Basin. Our results show improved recovery in both basement relief and susceptibility in the basement rock, and inversion with field data is able to identify three different susceptibility zones in basement rock below the Illinois Basin.