DOI QR코드

DOI QR Code

Evaluation of Inverse Fourier Integral Considering the Distances from the Source Point in 2D Resistivity Modeling

전기비저항탐사 2차원 모델링에서 송수신 간격을 고려한 푸리에 역변환

  • Cho, In-Ky (Division of Geology and Geophysics, Kangwon National University) ;
  • Jeong, Da-Bhin (Division of Geology and Geophysics, Kangwon National University)
  • 조인기 (강원대학교 지질.지구물리학부) ;
  • 정다빈 (강원대학교 지질.지구물리학부)
  • Received : 2017.10.11
  • Accepted : 2017.12.12
  • Published : 2018.02.28

Abstract

In the two-dimensional (2D) modeling of electrical method, the potential in the space domain is reconstructed with the calculated potentials in the wavenumber domain using inverse Fourier transform. The inverse Fourier integral is numerically evaluated using the transformed potential at different wavenumbers. In order to improve the precision of the integration, either the logarithmic or exponential approximation has been used depending on the size of wavenumber. Two numerical methods have been generally used to evaluate the integral; interval integration and Gaussian quadrature. However, both methods do not consider the distance from the current source. Thus the resulting potential in the space domain shows some error. Especially when the distance from the current source is very small or large, the error increases abruptly and the evaluated potential becomes extremely unstable. In this study, we developed a new method to calculate the integral accurately by introducing the distance from the current source to the rescaled Gauss abscissa and weight. The numerical tests for homogeneous half-space model show that the developed method can yield the error level lower than 0.4 percent over the various distances from the current source.

전기탐사 2차원 모델링에서는 다수의 파수영역 전위를 계산하고 이를 푸리에 역변환하여 공간영역 전위를 계산한다. 푸리에 역변환은 여러 개의 서로 다른 파수에서의 파수영역 전위를 사용하여 수치적으로 얻어진다. 적분의 정확도를 향상시키기 위하여 파수의 크기에 따라 적분 구간을 지수 근사와 대수 근사 구간으로 분할하는 방법이 널리 사용되고 있다. 푸리에 역변환에는 크게 구간 적분법과 가우스 적분법이 사용되고 있다. 그러나 이들 방법은 송수신 간격을 고려하지 못하므로 송수신 간격에 따른 오차를 피할 수 없다. 특히 송수신 간격이 매우 작거나 클 경우 오차가 급격하게 증가하는 문제점을 가지고 있다. 이 연구에서는 송수신 간격을 고려하여 가우스 좌표값 및 가중값을 적용하는 새로운 수치 적분법을 개발하였다. 반무한 공간에 대한 수치 실험 결과, 개발된 수치 적분법은 송수신 간격에 관계없이 0.4% 이하의 정밀도를 나타내었다.

Keywords

References

  1. Abramowitz, M., and Stegun, I. A., 1970, Handbook of mathematical functions, Dover publications.
  2. Dey, A., and Morrison, H. F., 1979, Resistivity modeling for arbitrarily shaped two-dimensional structures, Geophys. Prosp., 27, 106-136. https://doi.org/10.1111/j.1365-2478.1979.tb00961.x
  3. Kemna, A., Binley, A., and Slater, L., 2004. Crosshole IP imaging for engineering and environmental applications, Geophysics, 69, 97-107. https://doi.org/10.1190/1.1649379
  4. Kemna, A., 2000. Tomographic inversion of complex resistivity - Theory and application. PhD thesis, Ruhr-University of Bochum.
  5. Kim, H. J., 1992. Numerical quadrature techniques for inverse Fourier transform in two-dimensional resistivity modeling, J. Korea Inst. Mining Geol., 25, 73-77.
  6. LaBrecque, D. J., Miletto, M., Daily, W. D., Ramirez, A. L., and Owen, E., 1996, The effects of noise on Occam's inversion of resistivity tomography data, Geophysics, 61, 538-548. https://doi.org/10.1190/1.1443980
  7. Loke, M. H., Chambers, J. E., and Ogilvy, R. D., 2006, Inversion of 2D spectral induced polarization imaging data, Geophys. Prospect., 54, 287-301. https://doi.org/10.1111/j.1365-2478.2006.00537.x
  8. Queralt, P., Pous, J., and Marcuello, A., 1991, 2D resistivity modeling: An approach to arrays parallel to the strike direction, Geophysics, 56, 941-950. https://doi.org/10.1190/1.1443127
  9. Tripp, A. C., Hohmann, G. W., and Swift, C. M., Jr, 1984. Twodimensional resistivity inversion, Geophysics, 49, 1708-1717. https://doi.org/10.1190/1.1441578
  10. Weller, M., Seichter, M., and Kampke, A., 1996, Inducedpolarization modeling using complex electrical conductivities, Geophys. J. Int., 127, 387-398. https://doi.org/10.1111/j.1365-246X.1996.tb04728.x