Joint Inversion of DC Resistivity and Magnetic Data, Constrained by Cross Gradients, Compactness and Depth Weighting
- PDF / 9,051,684 Bytes
- 19 Pages / 547.087 x 737.008 pts Page_size
- 37 Downloads / 174 Views
Pure and Applied Geophysics
Joint Inversion of DC Resistivity and Magnetic Data, Constrained by Cross Gradients, Compactness and Depth Weighting RAMIN VARFINEZHAD,1 BEHROOZ OSKOOI,1 and MAURIZIO FEDI2 Abstract—In this paper we perform a 2-D joint inversion of DC resistivity and magnetic data, constrained by cross-gradients. Inspired by methods developed for potential fields, we introduce into both the separate and joint inversion algorithms also compactness and depth weighting functions, under the form of a model weighting-function. These constraints, usually not considered for DC resistivity inversion, reveal to be decisive for its joint inversion with magnetic data. A linear approximated forward problem of the resistivity is used for the joint inversion so that both the resistivity and magnetic problems are expressed as a linear integral equation under the form of a Fredholm integral of 1st kind. To examine the feasibility of the joint inversion algorithm, we first test the method with two synthetic cases: a thick dyke in a two-layered medium and a cavity located above a conductor. A third synthetic case involves a multisource model. The results are encouraging, revealing that the cross-gradient constraint is an effective tool to improve the separate inversions of DC resistivity and magnetic data. The joint inversion algorithm is also applied to data in the archeological area of the old Pompeii city, nearby Naples. Comparing the results of joint and separate inversions, we obtain a significant improvement in the interpretation of both kind of data in terms of buried walls of an ancient roman villa. In all the studied cases, the cross-gradient constraint appears to be a key-diagnostic tool to assess whether actual coherence is gained among DC resistivity and magnetic susceptibility models. Keywords: Cross-gradient function, joint inversion, magnetic, pompeii, resistivity.
1. Introduction Separate inversion of DC resistivity (e.g., Dimri 1992; Sasaki 1994; Loke and Barker 1995; Gu¨nther et al. 2006; Li and Oldenburg 2007; Hayley et al. 2011; Wiese et al. 2015; Gu¨ndog˘du and Candansayar
1
Department of Geomagnetism, Institute of Geophysics, University of Tehran, Tehran, Iran. E-mail: [email protected] 2 Department of Earth Sciences, Environment and Resources, University of Naples Federico II, Naples, Italy.
2018) and magnetic data (e.g., Dimri 1992; Li and Oldenburg 1996; Fedi and Rapolla 1999; Portniaguine and Zhdanov 1999; Silva et al. 2001; Pilkington 2009; Cella and Fedi 2012; Paoletti et al. 2013; Ganguli et al. 2016) are tools leading to realistic models of the related physical parameters: electrical resistivity and magnetic susceptibility, respectively. In order to build a valid geological model of the subsurface, it could be useful to integrate the information from different geophysical methods. This task may be fulfilled by joint interpretation (f.i., Orlando 2005; Ogaya et al. 2016), cooperative inversion (e.g., Takam Takougang et al. 2015; Le et al. 2016) and joint inversion (e.g., Gallardo and Meju
Data Loading...
