Full text loading...
-
Fast 3D Inversion of the Magnetic Vector Amplitude Produced by a Basement Relief Using Gauss-Legendre Quadrature
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 80th EAGE Conference and Exhibition 2018, Jun 2018, Volume 2018, p.1 - 5
Abstract
We present a fast 3D magnetic inversion algorithm for depth-to-basement estimate based on an efficient way to compute the amplitude of the magnetic vector components produced by an arbitrary interface separating nonmagnetic sedimentary rocks from a magnetic basement. We parametrized the basement layer as a grid of juxtaposed right rectangular prisms with a known homogeneous magnetization intensity. Because the amplitude data inversion is weakly dependent on the magnetization vector direction, our method does not require the knowledge about this direction. To produce a new forward modelling which is computationally faster than the one using a prism-based forward modelling, the x-, y-, and z-components of the magnetic vector of 3D prism are calculated numerically using the Gauss-Legendre quadrature (GLQ) produced by dipoles located along the vertical axis passing through the prism center. With this new forward modeling, we define a new sensitivity matrix which is simpler and computationally faster to compute than the equivalent matrix using a prism-based forward modelling. Test on synthetic amplitude data yields an acceptable data fit and a reasonable basement relief estimate.