Skip to main content
Log in

1-D DC Resistivity Modeling and Interpretation in Anisotropic Media Using Particle Swarm Optimization

  • Published:
Pure and Applied Geophysics Aims and scope Submit manuscript

Abstract

We examine the one-dimensional direct current method in anisotropic earth formation. We derive an analytic expression of a simple, two-layered anisotropic earth model. Further, we also consider a horizontally layered anisotropic earth response with respect to the digital filter method, which yields a quasi-analytic solution over anisotropic media. These analytic and quasi-analytic solutions are useful tests for numerical codes. A two-dimensional finite difference earth model in anisotropic media is presented in order to generate a synthetic data set for a simple one-dimensional earth. Further, we propose a particle swarm optimization method for estimating the model parameters of a layered anisotropic earth model such as horizontal and vertical resistivities, and thickness. The particle swarm optimization is a naturally inspired meta-heuristic algorithm. The proposed method finds model parameters quite successfully based on synthetic and field data. However, adding 5 % Gaussian noise to the synthetic data increases the ambiguity of the value of the model parameters. For this reason, the results should be controlled by a number of statistical tests. In this study, we use probability density function within 95 % confidence interval, parameter variation of each iteration and frequency distribution of the model parameters to reduce the ambiguity. The result is promising and the proposed method can be used for evaluating one-dimensional direct current data in anisotropic media.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12
Fig. 13
Fig. 14

Similar content being viewed by others

References

  • Anderson W. L. (1979), Numerical integration of related Hankel transforms of orders 0 and 1 by adaptive digital filtering, Geophysics 44, 1287–1305.

  • Bhattacharya P. K. and Patra H.P. (1968), Direct current geoelectric sounding: Principles and applications (Elsevier/Amsterdam).

  • Bhattacharya B. B. and Sen B. K. (1981), Depth of investigation of collinear electrode arrays over homogeneous anisotropic half-space in direct current methods, Geophysics 46, 768–780.

  • Christensen N. B. (1990), Optimized fast Hankel transform filters, Geophysical Prospecting 38, 545–558.

  • Dey A. and Morrison H. F. (1979), Resistivity modeling for arbitrarily shaped two-dimensional structures, Geophysical Prospecting 27, 106–136.

  • Fernandez–Martinez J. L. and Garcia–Gonzalo E. (2008), The generalized PSO: A new door to PSO, Journal of Artificial Evolution and Applications 2008, Article ID 861275.

  • Fernandez-Martinez J. L. and Garcia-Gonzalo E. (2009), The PSO family: deduction, stochastic analysis and comparison, Swarm Intelligence 3, 245–273.

  • Fernandez-Martinez J. L. and Garcia-Gonzalo E. (2011), Stochastic stability analysis of the linear continues and discrete PSO models, IEEE Transactions on Evolutionary Computations 15, 405–423.

    Google Scholar 

  • Fernandez-Martinez J. L., Garcia-gonzalo E., Fernandez-alvarez J. P., Kuzma H. A., and Menendez-perez C. O. (2010a), PSO: A powerful algorithm to solve geophysical inverse problems application to a 1D-DC case, Journal of Applied Geophysics 71, 13–25.

  • Fernandez-martinez J. L., Garcia-Gonzalo E., and Naudet V. (2010b), Particle swarm optimization applied to solving and appraising the streaming-potential inverse problem, Geophysics 75, WA3–WA15.

  • Fernandez-Martinez J. L., Fernandez-Muniz M. Z., and Tompkins M. J. (2012), On the topography of the cost functional in linear and nonlinear inverse problems, Geophysics 77, W1–W15.

  • Ghosh D. P. (1971a) The application of linear filter theory to the direct interpretation of geoelectrical resistivity measurements, Geophysical Prospecting 19, 192–217.

  • Ghosh D. P. (1971b), Inverse filter coefficients for the computation of apparent resistivity standard curves for a horizontally stratified earth, Geophysical Prospecting 19, 769–775.

  • Grant F. S. and West G. F. (1965) Interpretation theory in applied geophysics (New York:McGraw-Hill).

  • Greenhalgh S. A., Marescot L., Zhou B., Greenhalgh, M., and Wiese T. (2009), Electric potential and Frechet derivatives for a uniform anisotropic medium with a tilted axis of symmetry, Pure and Applied Geophysics 166, 673–699.

  • Guptasarma D. (1982), Optimization of short digital linear filters for increased accuracy, Geophysical Prospecting 30, 501–514.

  • Herwanger J. V., Pain C. C., Binley A., De Olivereira, C. R. E., and Wortington M. H. (2004), Anisotropic resistivity tomography, Geophysical Journal International 158, 409–425.

  • Johansen H. K. and Sorensen K. I. (1979), Fast Hankel transform, Geophysical Prospecting 27, 876–901.

  • Jupp D. L. B. and Vozoff K. (1977), Resolving anisotropy in layered media by joint inversion, Geophysical Prospecting 25, 460–470.

  • Kara I., Durdu M., and Arigün Z. (2009) Konya-Tuzlukçu (KT-1) Jeotermal Sondajı Kuyu Bitirme Raporu (MTA-Report [in Turkish]).

  • Keller G. V. and Frischknecht F. C. (1966), Electrical methods in geophysical prospecting (NewYork: Permagon).

  • Kennedy J. and Eberhart R. C. (1995), Particle swarm optimization, IEEE International Conference on Neural Networks 4, 1942–1948.

  • Key K. (2012), Is the fast Hankel transform faster than quadrature? Geophysics 77, F21–F30.

  • Kim J., Yi M., Cho S., Son J., and Song, W. (2006), Anisotropic cross hole resistivity tomography for ground safety analysis of a high-storied building over an abandoned mine, Journal of Environmental Engineering Geophysics 11, 225–235.

  • Koefoed O. (1979), Geosounding principles resistivity sounding measurements Vol. 14A (Amsterdam: Elsevier).

  • Kreyszig E. (1993), Advanced Engineering Mathematics (New York: Wiley).

  • Li Y. and Spitzer K. (2005), Finite element resistivity modelling for three-dimensional structures with arbitrary anisotropy, Physics of the Earth and Planetary Interiors 150, 15–27.

  • Maillet R. (1947), The fundamental equations of electrical prospecting, Geophysics 12, 529–556.

  • Matias M. J. S. and Habberjam G. M. (1986), The effect of structure and anisotropy on resistivity measurements, Geophysics 51, 964–971.

  • Pain C., Herwanger J., Saunders J., Wortington M., and Oliveira C. (2003), Anisotropic resistivity inversion, Inverse Problems 19, 1081–1111.

  • Pekşen E., Yas T., Kayman A. Y., and Özkan C. (2011), Application of particle swarm optimization on self-potential data, Journal of Applied Geophysics 75, 305–318.

  • Pervago E., Mousatov A., and Shevnin V. (2006), Analytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order, Geophysical Prospecting 54, 651–661.

  • Pidlisecky A. and Knight R. (2008), FW2_5D: A MATLAB 2.5-D electrical resistivity modeling Code, Computers & Geosciences 34, 1645–1654.

  • Sen M. and Stoffa P. L. (1995), Global optimization methods in geophysical inversion, (Amsterdam: Elsevier).

  • Shaw R. and Srivastava S. (2007), Particle swarm optimization: A new tool to invert geophysical data, Geophysics 72, F75–F83.

  • Shi Y. and Eberhart R. C. (1998), A modified particle swarm optimizer, IEEE International Conference on Evolutionary Computation (Anchorage, Alaska), 69–73.

  • Tarantola A. (1997), Inverse problem theory of data fitting and model parameter estimation (Amsterdam: Elsevier).

  • Telford W. M., Geldart L. P., and Sheriff R. E. (1990), Applied Geophysics (Cambridge: Cambridge).

  • Watson K. A. and Barker R. D. (1999), Differentiating anisotropy and lateral effects using azimuthal resistivity offset Wenner soundings, Geophysics 64, 739–745.

  • Watson K. A. and Barker R. D. (2010), Tank modeling of azimuthal resistivity surveys over anisotropic bedrock with dipping overburden, Near Surface Geophysics 8, 297–309.

  • Wiese T., Greenhalgh S. A. and Marescot L. (2009), DC sensitivity for surface patterns for tilted transversely isotropic media, Near Surface Geophysics 7, 125–139.

  • Xu S., Duan B., and Zhang D. (2000), Selection of the wavenumbers k using an optimization method for the inverse Fourier transform in 2.5D electrical modeling, Geophysical Prospecting 48, 789–796.

  • Yang X. S. (2010), Engineering optimization: An introduction with metaheuristic applications (New Jersey: Wiley).

  • Yin C. and Wiedelt P. (1999), Geoelectrical fields in layered earth with arbitrary anisotropy, Geophysics 64, 426–434.

  • Zhou B., Greenhalgh M., and Greenhalgh S. A. (2009), 2.5/3-D resistivity modelling in anisotropic media using Gaussian quadrature grids, Geophysical Journal International 176, 63–80.

Download references

Acknowledgments

This work was supported by the Scientific and Technical Research Council of Turkey (TÜBİTAK) under Grant 110Y354. We would like to thank the General Directorate of Mineral Research and Exploration (MTA) for providing us with the VES data and the well information. We would like to thank anonymous reviewers for constructive comments which helped improve the quality of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ertan Pekşen.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Pekşen, E., Yas, T. & Kıyak, A. 1-D DC Resistivity Modeling and Interpretation in Anisotropic Media Using Particle Swarm Optimization. Pure Appl. Geophys. 171, 2371–2389 (2014). https://doi.org/10.1007/s00024-014-0802-2

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00024-014-0802-2

Keywords

Navigation