TY - JOUR

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

AU - Rojas, Sergio

AU - Muga, Ignacio

AU - Pardo, David

N1 - Publisher Copyright:
© 2016, Springer International Publishing Switzerland.

PY - 2016/12/1

Y1 - 2016/12/1

N2 - 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.

AB - 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.

KW - Borehole measurements

KW - Hankel transform

KW - Planarly layered media

KW - Quadrature-free method

UR - http://www.scopus.com/inward/record.url?scp=84989172546&partnerID=8YFLogxK

U2 - 10.1007/s10596-016-9592-1

DO - 10.1007/s10596-016-9592-1

M3 - Article

AN - SCOPUS:84989172546

SN - 1420-0597

VL - 20

SP - 1301

EP - 1318

JO - Computational Geosciences

JF - Computational Geosciences

IS - 6

ER -