ERT, dam sounding, integral equations, triangulation, adaptive grid


To numerically solve the system of integral equations, it is customary to establish a discrete grid within each integration area. In the context of 3D modeling, these areas correspond to surfaces situated in space. The standard discretization technique employed for the computational domain is triangulation. This study addresses the integral equation system pertinent to the electrical tomography of dams. The structural model encompasses an embankment dam, the upstream and downstream water bodies, the dam's base, and a potential leakage region on the upstream side. An alternative configuration may be encountered in specific scenarios with no water downstream. Consequently, the model may incorporate up to nine distinct contact boundaries. Accordingly, the system of integral equations comprises an equivalent number of equations. Effectively resolving this system through numerical methods necessitates applying triangulation techniques to these diverse surfaces. While mathematical packages like Matlab offer triangulation functions, they may not fully address the specific demands of the problem. Additionally, the grid resolution should be heightened in proximity to key elements such as the sounding line, the supply electrode, and the various contact lines within the medium. These considerations transform the triangulation task into a distinct subtask within the numerical simulation of the resistivity tomography problem. In this paper, we provide our specific approach to this problem. The simplification of the triangulation algorithm is rooted in the predominant utilization of the two-dimensional geometric properties inherent to the object under study. For most contact boundaries, the triangulation is constructed layer by layer with a gradual modulation in triangle dimensions as one progresses from one layer to the next, orthogonal to the axis of the dam. Concerning the surface corresponding to the leakage area, cylindrical coordinates are used for surface parameterization. This approach enables partitioning the surface into discrete strata, facilitating a systematic, layer-by-layer grid construction. Additionally, points at the intersections of contact boundaries are integrated into the pre-existing triangulation by applying a standard function within the Matlab package. In the future, the mathematical modeling based on the Integral Equation Method with adaptive discretization will help incorporate real-time computations into information systems related to monitoring hydraulic structures.


Barker, R.D. (1981). The offset system of electrical resistivity sounding and its use with a multicore cable. Geophysical Prospecting, 29(1), 128-143.

Barker, R.D. (1993). A simple algorithm for electrical imaging of the subsurface. First Break, 10(2), 53-62.

Dahlin, T. (1993). On the automation of 2D resistivity surveying for engineering and environmental applications [PhD thesis]. Lund University.

Dahlin, T. (1996). 2D resistivity surveying for environmental and engineering applications. First Break, 14, 275–283.

Dahlin, T. (2001). The development of DC resistivity imaging techniques. Computers &Geosciences, 27, 1019–1029.

Dahlin, T. & B. Zhou, B. (2004). A numerical comparison of 2D resistivity imaging with 10 electrode arrays. Geophysical Prospecting, 52, 379-398.

Gunther, T. & Rucker, C. (2013). Boundless Electrical Resistivity Tomography. BERT 2 - the user tutorial, Ver.2.0.

Loke, M.H. & Barker, R.D. (1996). Rapid least-squares inversion of apparent resistivity pseudosections using a quasi-Newton method. Geophysical Prospecting, 44, 131-152.

Loke, M.H. (2000). Topographic modeling in electrical imaging inversion. 62nd Conference and Technical Exhibition, EAGE, Extended Abstracts, D-2.

Turarova, M., Mirgalikyzy, T., Mukanova, B., Modin, I., & Kaznacheev, P. (2022). Evaluation of the 3D Topographic Effect of Homogeneous and Inhomogeneous Media on the Results of 2D Inversion of Electrical Resistivity Tomography Data. Modelling and Simulation in Engineering.

Rakisheva, D., & Mukanova, B. (2021). Fourier Transformation method for solving integral equation in the 2.5 D problem of electric sounding. Journal of Physics: Conference Series, 2092(1), 012018.

Rakisheva, D., Mukanova, B., & Modin, I. (2020). Simulation Of Electrical Monitoring Of Dams With Leakage With A Transverse Placement Of The Measuring Installation. Eurasian Journal of Math. and Computer Appl., 8(4), 69-82.

Mukanova, B., & Modin, I. (2018). The Boundary Element Method in Geophysical Survey. Springer. ISBN 978-3-319-72908-4.

Shi-zhe, Xu. (2001). The Boundary Element Method in Geophysics. Geophysical Monograph Series, (9), SEG Books, 217.

Zhdanov, M., & Michael, S. (2015). Geophysical Inverse Theory and Applications. Elsevier, 2,730.

Turarova, M., Mirgalikyzy, T., Mukanova, B., & Modin, I. (2022). Elimination of the ground surface topographic effect in the 2d inversion results of electrical resistivity tomography data. Eurasian Journal of Math. and Computer Appl., 10(3), 84–104.

Frey, P.J. & George, P. (2008). Mesh Generation: Application to Finite Elements. ISTE Ltd,

Cuillière, JC. (1998). An adaptive method for the automatic triangulation of 3D parametric surfaces. Computer-Aided Design, 30(2), 139–149.

Lo, S. H. (2002). Finite element mesh generation and adaptive meshing. Progress in Structural Engineering and Materials, 4(4), 381–399.

Feng, L., Alliez, P., Busé, L., Delingette, H., & Desbrun, M. (2018). Curved optimal delaunay triangulation. ACM Transactions on Graphics, 37(4), 1-16.

Schmidt, P.P., & Kobbelt, L. (2023). Surface Maps via Adaptive Triangulations. Computer Graphics Forum, 42(2), 103 – 117.

Shewchuk, J. R. (2000). Delaunay refinement algorithms for triangular mesh generation. Computational Geometry: Theory and Applications, 22, (1-3), 21–74.

Meng, X., & Mei, G. (2013). Robust delaunay triangulation for domain with acute angles. Proceedings of SPIE - The International Society for Optical Engineering, 88782013(88783D) 5th International Conference on Digital Image Processing, ICDIP 2013.

Schier, A., & Klein, R. (2023). Discrete exterior calculus for meshes with concyclic polygons. Computer Aided Geometric Design, 101(C), 102170.

Shewchuk, J.R. (2012). Lecture Notes on Delaunay Mesh Generation. University of California at Berkeley, Berkeley. delnotes.dvi (

Engwirda, D. (2014). Locally-optimal Delaunay-refinement and optimisation-based mesh generation [Ph.D. Thesis, School of Mathematics and Statistics]. The University of Sydney.

Engwirda, D. (2005). Unstructured mesh methods for the Navier-Stokes equations. [Honours Thesis, School of Aerospace, Mechanical and Mechatronic Engineering]. The University of Sydney.

Mirgalikyzy, T., Mukanova, B., & Modin, I. (2015). Method of Integral Equations for the Problem of Electrical Tomography in a Medium with Ground Surface Relief. Journal of Applied Mathematics. DOI:

George, B.A., Hans, J.W., & Frank, E.H. (2012). Mathematical Methods for Physicists. Elsevier, 7. DOI:




How to Cite

Mukanova, B. (2023). SIMPLIFIED ADAPTIVE TRIANGULATION OF THE CONTACT BOUNDARIES OF THE DAM MODEL. Scientific Journal of Astana IT University, 15(15), 16–26.