Algorithm for determining inclusion parameters in solving inverse problems of geoelectrical exploration using the profiling method

Received: March 30, 2024
Lviv Polytechnic National University, Software Engineering Department
Carpathian Branch of Subbotin Institute of Geophysics of the National Academy of Sciences of Ukraine

The paper aims to develop an algorithm for recognizing the physical and geometric parameters of inclusion, using indirect methods of boundary, near-boundary, and partially-boundary elements based on the data of the potential field. Methodology. The direct and inverse two-dimensional problems of the potential theory concerning geophysics were solved when modeling the earth's crust with a piecewise-homogeneous half-plane composed of a containing medium and inclusion that are an ideal contact. To construct the integral representation of the solution of the direct problem, a special fundamental solution for the half-plane (Green's function) of Laplace's equation, which automatically satisfies the zero-boundary condition of the second kind on the day surface, and a fundamental solution for inclusion were used. To find the intensities of unknown sources introduced in boundary, near-boundary, or partially-boundary elements, the collocation technique was used, i.e. the conditions of ideal contact are satisfied in the middle of each boundary element. After solving the obtained SLAE, the unknown potential in the medium and inclusion and the flow through their boundaries are found, considering that the medium and inclusion are considered as completely independent domains. Results. The computational experiment for the task of electric prospecting with a constant artificial field using the resistance method, in particular, electrical profiling, was carried out, while focusing on the physical and geometric interpretation of the data. Initial approximations for the electrical conductivity of the inclusion, its center of mass, orientation and dimensions are determined by the nature of the change in apparent resistivity. To solve the inverse problem two cascades of iterations are organized: the first one is to specify the location of the local heterogeneity and its approximate dimensions, the second one is to specify its shape and orientation in space. At the same time, the minimization of the functional considered on the section of the boundary, where an excess of boundary conditions is set, is carried out. Originality. The method of boundary integral equations is shown to be effective for constructing numerical solutions of direct and inverse problems of potential theory in a piecewise homogeneous half-plane, using indirect methods of boundary, near-boundary, and partial-boundary elements as variants. Practical significance. The proposed approach for solving the inverse problem of electrical exploration with direct current is effective because it allows fora step-by-step, "cascade" recognition of the shape, size, orientation, and electrical conductivity of the inclusion. We follow the principle of not using all the details of the model and not attempting to recognize parameters with little effect on the result, especially with imprecise initial approximations.

  1. Brebbia, C. A., Telles, J. C. F., & Wrobel, L. C. (2012). Boundary element techniques: theory and applications in engineering. Springer Science & Business Media.
  2. Foks, N. L., & Krahenbuhl, R., Li, Y. (2014). Adaptive sampling of potential-field data: A direct approach to compressive inversion. Geophysics, 79(1), IM1–IM9
  3. Li, Z., Ma, G., Meng, Q., Wang, T., & Li, L. (2023). Gravity and magnetic fast inversion method with cross-gradient based on function fitting. Geophysical Journal International, 232(2), 1209-1218.
  4. Lv, M., Zhang, Y., & Liu, S. (2023). Fast forward approximation and multitask inversion of gravity anomaly based on UNet3+. Geophysical Journal International, 234(2), 972-984.
  5. Mikheeva T., Lapina O., Kyshman-Lavanova T., Prychepiy T. Geophysical data interpretation technologies in the study and exploration of oil-and-gas deposits. Geofizicheskiy Zhurnal, 2023. 44(5),106–122. (in Ukrainian).
  6. Milson J., Eriksen A. (2011). Field Geophysics: 4th ed. John Wiley & Sons, 287 p.
  7. Mukanova, B., & Modin, I. (2018). The boundary element method in geophysical survey. Springer International Publishing.
  8. Nikolaev, I. Y., Kushnir, A. M., Ilyenko, V. A., & Nikolaev, Y. I. (2019). Electromagnetic studies of the western part of the Ukrainian Shield. Geofizičeskij žurnal41(3), 120–133. (in Ukrainian).
  9. Pierre van Baal. (2014). A Course in Field Theory, 1st Edition, CRC Press,
  10. Qu, W., Chen, W., & Fu, Z. (2015). Solutions of 2D and 3D non-homogeneous potential problems by using a boundary element-collocation method. Engineering Analysis with Boundary Elements, 60, 2-9.
  11. Zhang, Y., Qu, W., & Chen, J. (2013). A new regularized BEM for 3D potential problems. Scientia Sinica Physica, Mechanica & Astronomica, 43(3), 297.
  12. Zhdanov, M. S. (2009). Geophysical electromagnetic theory and methods. Methods in Geochemistry and Geophysics, (vol 43), Elsevier.
  13. Zhou, S., Jia, H., Lin, T., Zeng, Z., Yu, P., & Jiao, J. (2023). An Accelerated Algorithm for 3D Inversion of Gravity Data Based on Improved Conjugate Gradient Method. Applied Sciences, 13(18), 10265. 
  14. Zhuravchak, L. M., & Zabrodska, N. V. (2021). Using of partly-boundary elements as a version of the indirect near-boundary element method for potential field modeling. Mathematical modeling and computing, 8(1), 1-10.
  15. Zhuravchak, L. M., & Zabrodska, N. V. (2021). Indentification of inhomogeneities according to potential field data using the near-boundary element method. 16th International Conference on Computer Sciences and Information Technologies (CSIT),
  16. Zhuravchak, L. M. (2019). Computation of pressure change in piecewise-homogeneous reservoir for elastic regime by indirect near-boundary element method. 14th International Conference on Computer Science and Information Technologies” (CSIT),  141-144.