To the solution of geoelectric problem on a direct current based on the global decomposition

The solution of a three-dimensional problem of geoelectric on a direct current is discussed. The solution of this task is constructed on the basis of global decomposition (GDAM). The algorithm consists in iterative reconstruction of the decision from two basic elements — field of a point source in the layered (background) environment and solution in a local insertion. As a background (normal) model, the horizontally layered environment is used, in one of layers of which the point source of direct current is located. The background environment contains limited inclusion with enough arbitrary distribution of specific electrical resistance. The problem consists in studying of influence of an insertion on the normal field. The work of the algorithm of the solution of a task is illustrated on the example of a homogeneous half-space (the three-layer horizontally layered environment with identical conductivity of layers). The point source is located in the first layer. The second layer contains a homogeneous local body (a rectangular parallelepiped or a finite cylinder). Solutions of tasks of Dirikhle for Laplace’s equation in these bodies are known. For the improving of the accuracy of calculations in a parallelepiped, the little change of an algorithm of calculation of potential was needed. Acceleration of calculations was provided by the calculations of Fourier coefficients with the use of direct discrete sine transformation.On the basis of the algorithms given in the work, the programs of calculation of potential in local insertions are developed and the numerical experiments which allowed presenting graphically calculation results and making assessments of accuracy of calculations, are executed.

1. Budak B.M., Samarsky A.A., Tihonov A.N. Collection of problems in mathematical physics. M., Fizmatlit Publ., 2003, 688 p. (In Russian).

2. Van’yan L.L., Debabov A.S., Yudin M.N. Interpretation of magnetotelluric soundings of inhomogeneous media. M., Nedra Publ., I984, 197 p. (In Russian).

3. Zaborovsky A.I. Electrical prospecting. M., Gostoptekhizdat Publ., 1963, 423 p. (In Russian).

4. Serebrennikova N.N. Numerical simulation of the electric field of a point source in a layered medium with local inclusions. Ph.D. thesis, Perm, 1987, 127 p. (In Russian).

5. Wait J.R. Geoelectromagnetism. M., Nedra Publ., 1987, 235 p. (In Russian).

6. Yudin M.N., Spassky B.A On the assessment of the modified Schwarz‘s algorithm convergence in VES direct problem. Geophysics, 2016, no 5, pp. 23—31. (In Russian).

7. Yudin M.N., Daev D.S., Karinsky A.D., Sevost’yanov N.A. The Schwartz algorithm. Calculation the anomalous potential Fourier-spectr in a layered medium containing local bodies. Collection of scientific papers «Theory and practice of exploration and petroleum geophysics», Perm, 2018, pp. 312—317. (In Russian).

8. Yudin M.N., Rumyanceva A.A. Multiscale approach to modeling electromagnetic fields in geoelectrics. Proceedings of the Inter-national Scientific and Practical Conference «New Ideas in Earth Sciences» MGRI-RGGRU, April 9—12, 2013, Moscow, 2013, Vol. 3, pp. 261—264. (In Russian).