TY - JOUR
T1 - A numerical 1.5D method for the rapid simulation of geophysical resistivity measurements
AU - Shahriari, Mostafa
AU - Rojas, Sergio
AU - Pardo, David
AU - Rodríguez-Rozas, Angel
AU - Bakr, Shaaban A.
AU - Calo, Victor M.
AU - Muga, Ignacio
N1 - Funding Information:
Acknowledgments: This work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 777778, the Projects of the Spanish Ministry of Economy and Competitiveness with reference MTM2016-76329-R (AEI/FEDER, EU), and MTM2016-81697-ERC/AEI, the BCAM “Severo Ochoa” accreditation of excellence SEV-2013-0323, and the Basque Government through the BERC 2014–2017 program, the Consolidated Research Group Grant IT649-13 on “Mathematical Modeling, Simulation, and Industrial Applications (M2SI)”. This publication was also made possible in part by the CSIRO Professorial Chair in Computational Geoscience at Curtin University and the Deep Earth Imaging Enterprise Future Science Platforms of the Commonwealth Scientific Industrial Research Organisation, CSIRO, of Australia. Additional support was provided, the Mega-grant of the Russian Federation Government (N 14.Y26.31.0013), and The Institute for Geoscience Research (TIGeR) and the Curtin Institute for Computation at Curtin University, and the Chilean FONDECYT research project No 1160774.
Funding Information:
This work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 777778, the Projects of the Spanish Ministry of Economy and Competitiveness with reference MTM2016-76329-R (AEI/FEDER, EU), and MTM2016-81697-ERC/AEI, the BCAM “Severo Ochoa” accreditation of excellence SEV-2013-0323, and the Basque Government through the BERC 2014–2017 program, the Consolidated Research Group Grant IT649-13 on “Mathematical Modeling, Simulation, and Industrial Applications (M2SI)”. This publication was also made possible in part by the CSIRO Professorial Chair in Computational Geoscience at Curtin University and the Deep Earth Imaging Enterprise Future Science Platforms of the Commonwealth Scientific Industrial Research Organisation, CSIRO, of Australia. Additional support was provided, the Mega-grant of the Russian Federation Government (N 14.Y26.31.0013), and The Institute for Geoscience Research (TIGeR) and the Curtin Institute for Computation at Curtin University, and the Chilean FONDECYT research project No 1160774.
Publisher Copyright:
© 2018 by the authors. Licensee MDPI, Basel, Switzerland.
PY - 2018/6
Y1 - 2018/6
N2 - In some geological formations, borehole resistivity measurements can be simulated using a sequence of 1D models. By considering a 1D layered media, we can reduce the dimensionality of the problem from 3D to 1.5D via a Hankel transform. The resulting formulation is often solved via a semi-analytic method, mainly due to its high performance. However, semi-analytic methods have important limitations such as, for example, their inability to model piecewise linear variations on the resistivity. Herein, we develop a multi-scale finite element method (FEM) to solve the secondary field formulation. This numerical scheme overcomes the limitations of semi-analytic methods while still delivering high performance. We illustrate the performance of the method with numerical synthetic examples based on two symmetric logging-while-drilling (LWD) induction devices operating at 2 MHz and 500 KHz, respectively.
AB - In some geological formations, borehole resistivity measurements can be simulated using a sequence of 1D models. By considering a 1D layered media, we can reduce the dimensionality of the problem from 3D to 1.5D via a Hankel transform. The resulting formulation is often solved via a semi-analytic method, mainly due to its high performance. However, semi-analytic methods have important limitations such as, for example, their inability to model piecewise linear variations on the resistivity. Herein, we develop a multi-scale finite element method (FEM) to solve the secondary field formulation. This numerical scheme overcomes the limitations of semi-analytic methods while still delivering high performance. We illustrate the performance of the method with numerical synthetic examples based on two symmetric logging-while-drilling (LWD) induction devices operating at 2 MHz and 500 KHz, respectively.
KW - Finite element method
KW - Hankel transform
KW - Logging-while-drilling (LWD)
KW - Multi-scale method
KW - Resistivity measurements
KW - Secondary field
UR - http://www.scopus.com/inward/record.url?scp=85048999972&partnerID=8YFLogxK
U2 - 10.3390/geosciences8060225
DO - 10.3390/geosciences8060225
M3 - Article
AN - SCOPUS:85048999972
VL - 8
JO - Geosciences (Switzerland)
JF - Geosciences (Switzerland)
SN - 2076-3263
IS - 6
M1 - 225
ER -