A quadrature-free method for simulation and inversion of 1.5D direct current (DC) borehole measurements

Sergio Rojas, Ignacio Muga, David Pardo

Research output: Contribution to journalArticlepeer-review

4 Scopus citations


Resistivity inverse problems are routinely solved in order to characterize hydrocarbon bearing formations. They often require a large number of forward problems simulations. When considering a one dimensional (1D) planarly layered media, semi-analytical methods can be employed in order to solve a single forward problem in a fraction of a second. However, in some situations, a large number of (over one million) simulations is required, preventing this method to be used as a real time (logging) alternative. In this paper, we propose a novel semi-analytical method that dramatically reduces the total computational time, so it can be employed for real time inversion. In our proposed method, we select an ad hoc basis representation for the spectral solution such that its inverse Hankel transform can be computed analytically. The proposed method requires a pre-process that is expensive when compared with a single evaluation in classical semi-analytical methods. However, subsequent evaluations can be rapidly obtained, decreasing thus the total computational time by orders of magnitude when the number of required forward simulations is large.

Original languageEnglish
Pages (from-to)1301-1318
Number of pages18
JournalComputational Geosciences
Issue number6
StatePublished - 1 Dec 2016


  • Borehole measurements
  • Hankel transform
  • Planarly layered media
  • Quadrature-free method


Dive into the research topics of 'A quadrature-free method for simulation and inversion of 1.5D direct current (DC) borehole measurements'. Together they form a unique fingerprint.

Cite this