Standard

Treatment of singularity in the method of boundary integral equations for 2.5D electromagnetic modeling. / Dyatlov, Gleb; Kushnir, Dmitry; Dashevsky, Yuliy.

In: Geophysics, Vol. 82, No. 2, 01.03.2017, p. E57-E75.

Research output: Contribution to journalArticlepeer-review

Harvard

APA

Vancouver

Dyatlov G, Kushnir D, Dashevsky Y. Treatment of singularity in the method of boundary integral equations for 2.5D electromagnetic modeling. Geophysics. 2017 Mar 1;82(2):E57-E75. doi: 10.1190/GEO2015-0645.1

Author

Dyatlov, Gleb ; Kushnir, Dmitry ; Dashevsky, Yuliy. / Treatment of singularity in the method of boundary integral equations for 2.5D electromagnetic modeling. In: Geophysics. 2017 ; Vol. 82, No. 2. pp. E57-E75.

BibTeX

@article{dfaccae6352145f1b9fab7c418ed9abb,
title = "Treatment of singularity in the method of boundary integral equations for 2.5D electromagnetic modeling",
abstract = "Fast and accurate simulation of responses of electromagnetic logging-while-drilling tools in 2D and 3D formations is important for further development of the proactive geosteering service,which is currently based on the 1D layered-model inversion. In search of an optimal method of forward simulation, we have developed boundary integral equations in a 2D formation model with planar boundaries. A peculiarity of the method is that the number of grid nodes required for reaching a certain accuracy,and hence the computation time increases as the transmitter approaches the boundary. To retain efficiency in the case of transmitter-boundary proximity, we use the following modification of the method. The main term of the calculated electromagnetic field, namely, the solution for the two-layered model corresponding to the closest boundary, is calculated by the well-known explicit formulas, whereas the anomalous field is found from the system of integral equations. The modified method is implemented and tested on complex realistic 2D formation models as well as basic models to be used in 2D inversion. Numerical simulations of real propagation and azimuthal resistivity tools show that software based on the modified method can be successfully used for synthetic log simulation in complex formations.However, the proposed method has two shortcomings: poor efficiency when the transmitter is close to the boundary endpoints and limitation to isotropic formations. Both are topics for further research. Currently, the singularity near the endpoints is treated by simply refining the grid, which deteriorates efficiency.Generalization to the anisotropy case would include modification of the basic method and derivation of the closed-form solution for the layered model without the assumption of transverse isotropy.",
author = "Gleb Dyatlov and Dmitry Kushnir and Yuliy Dashevsky",
note = "Publisher Copyright: {\textcopyright} 2017 Society of Exploration Geophysicists. Copyright: Copyright 2020 Elsevier B.V., All rights reserved",
year = "2017",
month = mar,
day = "1",
doi = "10.1190/GEO2015-0645.1",
language = "English",
volume = "82",
pages = "E57--E75",
journal = "Geophysics",
issn = "0016-8033",
publisher = "SOC EXPLORATION GEOPHYSICISTS",
number = "2",

}

RIS

TY - JOUR

T1 - Treatment of singularity in the method of boundary integral equations for 2.5D electromagnetic modeling

AU - Dyatlov, Gleb

AU - Kushnir, Dmitry

AU - Dashevsky, Yuliy

N1 - Publisher Copyright: © 2017 Society of Exploration Geophysicists. Copyright: Copyright 2020 Elsevier B.V., All rights reserved

PY - 2017/3/1

Y1 - 2017/3/1

N2 - Fast and accurate simulation of responses of electromagnetic logging-while-drilling tools in 2D and 3D formations is important for further development of the proactive geosteering service,which is currently based on the 1D layered-model inversion. In search of an optimal method of forward simulation, we have developed boundary integral equations in a 2D formation model with planar boundaries. A peculiarity of the method is that the number of grid nodes required for reaching a certain accuracy,and hence the computation time increases as the transmitter approaches the boundary. To retain efficiency in the case of transmitter-boundary proximity, we use the following modification of the method. The main term of the calculated electromagnetic field, namely, the solution for the two-layered model corresponding to the closest boundary, is calculated by the well-known explicit formulas, whereas the anomalous field is found from the system of integral equations. The modified method is implemented and tested on complex realistic 2D formation models as well as basic models to be used in 2D inversion. Numerical simulations of real propagation and azimuthal resistivity tools show that software based on the modified method can be successfully used for synthetic log simulation in complex formations.However, the proposed method has two shortcomings: poor efficiency when the transmitter is close to the boundary endpoints and limitation to isotropic formations. Both are topics for further research. Currently, the singularity near the endpoints is treated by simply refining the grid, which deteriorates efficiency.Generalization to the anisotropy case would include modification of the basic method and derivation of the closed-form solution for the layered model without the assumption of transverse isotropy.

AB - Fast and accurate simulation of responses of electromagnetic logging-while-drilling tools in 2D and 3D formations is important for further development of the proactive geosteering service,which is currently based on the 1D layered-model inversion. In search of an optimal method of forward simulation, we have developed boundary integral equations in a 2D formation model with planar boundaries. A peculiarity of the method is that the number of grid nodes required for reaching a certain accuracy,and hence the computation time increases as the transmitter approaches the boundary. To retain efficiency in the case of transmitter-boundary proximity, we use the following modification of the method. The main term of the calculated electromagnetic field, namely, the solution for the two-layered model corresponding to the closest boundary, is calculated by the well-known explicit formulas, whereas the anomalous field is found from the system of integral equations. The modified method is implemented and tested on complex realistic 2D formation models as well as basic models to be used in 2D inversion. Numerical simulations of real propagation and azimuthal resistivity tools show that software based on the modified method can be successfully used for synthetic log simulation in complex formations.However, the proposed method has two shortcomings: poor efficiency when the transmitter is close to the boundary endpoints and limitation to isotropic formations. Both are topics for further research. Currently, the singularity near the endpoints is treated by simply refining the grid, which deteriorates efficiency.Generalization to the anisotropy case would include modification of the basic method and derivation of the closed-form solution for the layered model without the assumption of transverse isotropy.

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

U2 - 10.1190/GEO2015-0645.1

DO - 10.1190/GEO2015-0645.1

M3 - Article

AN - SCOPUS:85013820660

VL - 82

SP - E57-E75

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 2

ER -

ID: 26146606