Geoelectric Model Research Articles

Research and application of joint-constrained inversion of transient electromagnetic multivariate parameter

Due to the phenomena of stratigraphic inclination, complex structure, and lateral discontinuity of resistivity or layer thickness in most of the coal seams, the traditional one-dimensional transient electromagnetic inversion method has limitations in interpretation accuracy. In addition, two- and three-dimensional inversion and artificial intelligence inversion have problems of large computation and large sample size, respectively, which limit their application in small- and medium-sized engineering exploration. To improve the inversion effect, this study proposes a method of joint-constrained inversion of transient electromagnetic multivariate parameters. This method achieves the joint constraint inversion of the transient electromagnetic multi-parameter by making full use of the geological data and a priori information to construct the initial model and adding the constraints such as the resistivity, the thickness, and the layer interface of each layer in the inversion objective function, and at the same time, taking into account the spatial correlation of the stratigraphic structure between the neighboring measurement points, as well as the transverse and vertical constraints between the measurement points along the direction of the survey line and perpendicular to the survey line. First, a series of typical geoelectric models are established and numerically simulated, and the results are compared with those of the traditional inversion method to verify the applicability and effectiveness of the method. Then, the constrained inversion is carried out on the physical simulation and measured data, and the results are in good agreement with the actual geological conditions. The numerical simulation, physical simulation and measured data inversion results consistently prove that this method can effectively reduce the uncertainty of the inversion at the isolated measuring points, improve the spatial continuity of the formation boundary, and better reflect the actual geoelectric characteristics of the formation.

Evaluation of the hydrogeological response to neotectonics of the Quequén Grande River Watershed using geoelectrical methods

The Quequén Grande River Watershed (Southeast of Buenos Aires Province, Argentina) is a strong asymmetrical watershed in the Pampean Plain. It is marked by the stream capture of the Quequén Grande River eastward. The topographic slopes and hydraulic gradients are typically low in this area, with streams and rivers primarily effluent along their courses. This behavior is attributable to the influence of the unconfined and shallow Pampean aquifer, which serves as the region's main source of water supply. In the southwestern limit of this drainage basin, several shallow lakes, disconnected from the surface drainage network and related to an old deflation/accumulation relief, occur over hillock formations. In this area, Paleozoic rocks crop out and constitute the hydrological basement of the aquifer. The bedrock architecture of this formation could explain the stream piracy of the Quequén Grande River eastward and the local geomorphology, where structural lineaments related to neotectonic movements have been reported. However, there is a lack of information related to the hydrological basement in depth, which is needed to understand the surface water dynamics of this drainage basin, linked to the topographic control of the landscape. To fill this gap, a geoelectric model of the hydrological basement is proposed from 37 Vertical Electrical Soundings and 6 Electric Resistivity Tomographies, aiming to contribute to understanding this portion of the Earth's crust. Geoelectrical surveys were performed over 3000 km2 of the southwestern limit of the Quequén Grande River Watershed. Results obtained and field observations suggest the deviation of the Quequén Grande River catchment eastward is caused by the hydrogeological basement structure in depth, in response to neotectonic movements inferred by structural lineaments.

Study on ground-penetrating radar wave field characteristics for earth dam disease considering the medium randomness

Ground-Penetrating Radar (GPR) has been widely used for non-destructive testing of earth dam disease. However, the forward simulation of GPR for earth dam disease often employs layered homogeneous models, neglecting the influence of medium randomness on its wave field characteristics. Therefore, considering the randomness of the medium, a geoelectrical model for earth dam disease is established, which is based on the mixed-type autocorrelation function and the finite element time-domain method. The influence of random medium model parameters on the single-channel wave of GPR is analyzed. The electromagnetic wave propagation characteristics under different medium models are explored. The forward simulation of GPR for earth dam disease such as panel voiding, concentrated seepage, and loosening are performed. The differences in propagation characteristics for earth dam disease between uniform medium model and random medium model are compared. Compared to the calculation results of the uniform medium model, the propagation speed and amplitude of electromagnetic waves in the random medium model changes, and a number of diffraction waves are present. When performing forward simulation of GPR for earth dam disease, considering medium randomness can deepen the understanding of the GPR section view and help improve the accuracy of image interpretation.

A Deep Learning Inversion Method for Airborne Time-Domain Electromagnetic Data Using Convolutional Neural Network

Due to the high detection efficiency of the airborne time-domain electromagnetic method, it can quickly collect electromagnetic response data for large area-wide regions, but it also brings great challenges to the inversion interpretation of the data because there are numerous survey data that need to be inverted. Conventional optimal inversion and fast imaging methods still take a long time to obtain conductivity and depth information, which affect the efficiency of real-time data interpretation. In this paper, we present a deep learning inversion method that can be used to solve the fast inversion problem of airborne time-domain electromagnetic data; the method uses a one-dimensional convolutional neural network. The network structure consists of two parts containing different numbers of convolutional and pooling layers. The training sample dataset was generated via two ways of constructing geoelectric models through forward modelling. To check the effectiveness of our deep learning inversion strategy, we tested it on synthetic data and two types of survey data. The experimental results show that this inversion method is effective and that it can be applied to airborne time-domain electromagnetic data collected using different observation systems. The proposed inversion method can obtain better inversion results for both simple and complex stratigraphic structures and requires significantly less computation time compared to conventional optimal inversion methods.

Transient Electromagnetic Monitoring of the Cryolithozone: Mathematical Modeling

Abstract ––The paper presents theoretical results on mathematical modeling and numerical inversion of transient electromagnetic data for the monitoring of the state of the cryolithozone. Solutions to the direct and inverse problems are considered based on the vector finite element method, the Sumudu integral transform and the apparatus of artificial neural networks. We show the capability of spatial localization of thaw zones (taliks) with transient cross-borehole exploration systems in a geoelectric model that takes into account the dispersive properties of frozen rocks. The performance and accuracy of the developed algorithms are assessed.

A potential deep geothermal reservoir in eastern Ladakh as inferred from the upper crustal geoelectric structure of the region

The Ladakh Himalaya is a repository of the evolutionary history of the India-Eurasia convergent plate margin. We present the results of a magnetotelluric (MT) study carried out in Eastern Ladakh along a 40-km-long profile traversing across various tectonic domains of this plate margin. We modeled the MT dataset by two-dimensional (2-D) and three-dimensional (3-D) inversion algorithms. The upper crustal geoelectric model obtained by 2-D inversion of distortion corrected and decomposed data of 18 sites shows an excellent match with the surface geology. In this model, the Indus Suture Zone (ISZ) appears as a steeply dipping contact and the Ladakh batholith is about 10 to 12 km thick. The model yields a prominent electrical conductor of ∼6 km diameter at 4 km depth beneath the Tso Morari Crystallines with an offshoot of a small conductor that rises upward along the ISZ. The geoelectrical model obtained by 3-D inversion is broadly consistent with the 2-D model but also reveals off-profile features and a connectivity of the shallow conductors with a widespread conductive zone at deeper level of the upper crust. In view of the presence of Puga and Chumathang hot springs towards west and Demchok hot springs towards east of our study region, we infer these conductors to be representing a potential major geothermal reservoir system connected to the widespread deeper conductive zone of partial melts, which extends from the southern Tibet to the eastern Ladakh region.

A study on geoelectrical recognition model of seawater/freshwater interface based on convolutional neural network: an application in sand tank experiments

A study on geoelectrical recognition model of seawater/freshwater interface based on convolutional neural network: an application in sand tank experiments

3-D parallel anisotropic inversion of controlled-source electromagnetic data using nested tetrahedral grids

SUMMARY Geophysicists today face the challenge of quickly and reliably interpreting extensive controlled-source electromagnetic (CSEM) data sets to map subsurface conductivity structures within realistic geological environments. An ideal 3-D CSEM inversion algorithm using tetrahedral grids should be capable of distinguishing different resolution requirements between forward modelling and inversion grids, have an optimal parallel strategy that fully exploits the inherent independence of CSEM data sets while also possessing the capability to handle large-scale geo-electrical models, and incorporate conductivity anisotropy which should be a common characteristic in realistic subsurface environments. However, existing tools in the geo-electromagnetic community often fall short of these three demands. Addressing this gap, our study introduces a scalable and parallel anisotropic inversion technique for CSEM data, capitalizing on the potential of unstructured tetrahedral grids. We first apply the tetrahedral longest-edge bisection method to create a refined dense, heterogeneous forward modelling grid from a coarse inversion grid. This refinement, focused on areas around transmitters and receivers, is seamlessly integrated within the coarser inversion grid’s topology, enabling precise conductivity mapping and preserving electromagnetic response accuracy during model updates. We further innovate with a source-mesh double-level parallel strategy, utilizing the message passing interface technique for parallel handling of independent CSEM data sets and large-scale geo-electrical models. Externally, we dedicate a processor for inversion model updates employing the Limited-memory Broyden–Fletcher–Goldfarb–Shanno optimization algorithm and divide other processors into groups, each associated with specific transmitting sources and frequencies. Internally, in each group, we employ a domain-decomposition-based scalable and robust iterative solvers using the Auxiliary-Space Maxwell pre-conditioner to parallel quickly calculate the electromagnetic responses from its assigned source-frequency set. Additionally, recognizing the potential for electrical conductivity anisotropy in field data, we incorporate the case of vertical transverse isotropy. We validate the effectiveness of our method through examples, including an isotropic land model with undulating topography, an anisotropic marine model and a real-field data case. Results from both synthetic and field data inversions underscore our method’s significant advancements in efficiency and practicality, particularly in addressing large-scale 3-D CSEM data sets inversion challenges in realistic geological environments.

Global insights on flood risk mitigation in arid regions using geomorphological and geophysical modeling from a local case study

This research provides a comprehensive examination of flood risk mitigation in Saudi Arabia, with a focus on Wadi Al-Laith. It highlights the critical importance of addressing flood risks in arid regions, given their profound impact on communities, infrastructure, and the economy. Analysis of morphometric parameters ((drainage density (Dd), stream frequency (Fs), drainage intensity (Di), and infiltration number (If)) reveals a complex hydrological landscape, indicating elevated flood risk. due to low drainage density, low stream frequency, high bifurcation ratio, and low infiltration number. Effective mitigation strategies are imperative to protect both communities and infrastructure in Wadi Al-Laith. Geophysical investigations, using specialized software, improve the quality of the dataset by addressing irregularities in field data. A multi-layer geoelectric model, derived from vertical electrical sounding (VES) and time domain electromagnetic (TDEM) surveys, provides precise information about the geoelectric strata parameters such as electrical resistivity, layer thicknesses, and depths in the study area. This identifies a well-saturated sedimentary layer and a cracked rocky layer containing water content. The second region, proposed for a new dam, scores significantly higher at 56% in suitability compared to the first region’s 44%. The study advocates for the construction of a supporting dam in the second region with a height between 230 and 280 m and 800 m in length. This new dam can play a crucial role in mitigating flash flood risks, considering various design parameters. This research contributes to flood risk management in Saudi Arabia by offering innovative dam site selection approaches. It provides insights for policymakers, researchers, and practitioners involved in flood risk reduction, water resource management, and sustainable development in arid regions globally.

Mineralization based on CSAMT data: A case study on the Gaoloushan section of Yangshan gold belt

AbstractYangshan gold belt (YSGB) in Gansu Province is a super‐large orogenic gold deposit located in the middle section of the Mianxian‐Lueyang arcuate suture zone (Mian‐Lue suture zone). The YSGB is mainly divided into the Gejiaowan, Anba, and Gaoloushan mine sections. The ore veins are strictly controlled by strong strain zones in the fracture zone and secondary fractures superimposed on the inherited activity and are vein‐like in general and lenticular in localities. For a long time, there have been different understandings regarding its ore‐controlling factors and mineral genesis. By utilizing controlled‐source audio‐frequency magnetotellurics (CSAMT) technology to conduct electrical characteristic studies, the deep electrical structural characteristics of the Gaoloushan section are identified, and a tectonically‐controlled geoelectric model of the deposit is established, which provides geophysical evidence for the view that the Yangshan gold mine is a tectonically‐controlled low‐temperature hydrothermal deposit in altered rock, guiding the arrangement of exploration work and identifying prospective mining areas, and thereby providing a geophysical exploration paradigm for deep exploration of large orogenic gold deposits.

Airborne transient electromagnetic imaging method based on wavelet neural network

The airborne transient electromagnetic method is a crucial technique for near-surface exploration in rugged terrains. In this study, we introduce a multi-input, single-output double-hidden wavelet neural network design for airborne transient electromagnetic pseudo-resistivity imaging. We construct uniform half-space, stratified stratum, and three-dimensional geoelectric models. By evaluating various performance indicators of neural networks that employ different wavelet basis functions as activation functions, we identify the most suitable wavelet basis functions. The quasi-resistivity is computed using both the wavelet neural network and the backpropagation neural network and then juxtaposed with traditional apparent resistivity. Our findings indicate that the wavelet neural network's quasi-resistivity aligns more closely with the resistivity of the real model. It is also more responsive to low resistivity anomalies than the conventional apparent resistivity translation algorithm. The wavelet approach moderates the undershoot or overshoot occurrences during abrupt stratum changes, offering a more accurate representation of subterranean electrical properties. When compared to the backpropagation neural network, the wavelet neural network provides a superior fit for the model, rendering a smoother quasi-resistivity depth curve. Therefore, it stands out as an improved method for pseudo-resistivity imaging. By processing survey data via the trained wavelet neural network, we find that the all-time apparent resistivity and quasi-resistivity align well with real-world situations. The wavelet neural network prominently showcases the low-resistance calculation results, offering a broader resistivity range. This clarity enhances anomaly detection, confirming the wavelet neural network's applicability and practicality for airborne transient electromagnetic imaging.

Анализ эквивалентности трехмерных моделей при 3D-интерпретации данных аэроэлектроразведки

The article is devoted to the problem of the equivalence of three-dimensional geoelectric models obtained as a result of 3D-interpretation of aerial electrical survey data in medium where local conductive bodies are covered by a heterogeneous near-surface layer. To carry out the research, we used data from aerial electrical prospecting carried out in the Creighton area (Canada) when searching for polymetallic ores. A two-stage 3D-inversion was used to build the geoelectric model. The purpose of the first stage of 3D-inversion was to reconstruct near-surface structures based on survey data obtained at the early times. The second stage of 3D-inversion was based on the residual map constructed from the practical signals and the signals obtained after the first stage of 3D-inversion at later times. Using the residual map, the starting positions of local bodies under the overlying layer were specified, and their geometry and properties were obtained during local 3D-inversions. Based on the results of 3D-inversion, a local conductive body was identified and confirmed by drilling data. An equivalence analysis was carried out for a three-dimensional model of this body. Various options for initial models differing in the number of blocks describing the body were considered, as well as options for changing the size of blocks in the model obtained as a result of 3D-inversion. It was found that the main equivalence is related to the large distance between flight lines, which is several times greater than the distance between positions on the profile. Therefore, most equivalent models are characterized by a decrease in the size of blocks in the direction orthogonal to the profiles, due to an increase in their size along the profile and a change in resistivity. At the same time, according to the criterion for drilling wells, the resulting equivalent models were identical, i.e. from a geological point of view, no fundamental difference was found between the resulting equivalent models.

THEORETICAL ANALYSIS OF THE APPLICATION OF PSEUDO-NOISE SIGNALS IN MARINE PULSE ELECTROMAGNETIC PROSPECTING

Using numerical modeling of a useful signal based on a characteristic geoelectric model of underwater permafrost and real recordings of a noise signal, we conducted a comparative analysis of noise suppression during transient sounding using a marine towed dipole-dipole array in accumulation modes with opposite-polar current pulses and pseudo-noise signals (PNS). For a series of current sequences in the form of PNS with different durations and numbers of pulses, as well as for a signal in the accumulation mode, with the superposition of an identical noise signal, transient sounding curves were obtained corresponding to the geoelectric model under consideration at a recording time of the order of 8-13 s per station. Based on the results of comparison of the obtained curves reconstructed from noisy synthetic data in the accumulation and PNS modes, it was established that in the PNS mode with pulses of duration 100 μs and 1 ms, the relative error in the PNS mode on average over the profile turns out to be significantly lower (up to 1.5 times) than a similar one error in the accumulation mode, and remains within an acceptable value (up to several percent) until later times.

RESULTS OF DEEP MAGNETOTELLURIC SURVEYS ALONG THE PUSHKINSKIYE GORY - ANDREAPOL PROFILE

In january-february of 2022 deep magnetotelluric soundings (DMTS) were carried out along the regional (300 km) profile Pushkinskie Gory - Andreapol. The field expedition included professors, graduates and students of the Geophysical department of the Geological Faculty of Lomonosov Moscow State University, as well as scientists from the Laboratory of Magnetotelluric Soundings of the Center of Geoelectromagnetic Research IPE RAS. The article presents an analysis of the results, as well as geoelectric two-dimensional models obtained during the inversion of magnetotelluric data in the ZONDMT2D program.

Frequency-dependent shear wave attenuation across the Central Anatolia region, Türkiye

Abstract. The Central Anatolian Plateau with its volcanic provinces represents a broad transition zone between the compressional deformation in the east and the extensional regime in the west. The Central Anatolian Fault Zone separates the Kırşehir Block in the north and the Anatolide–Tauride Block in the south within the plateau. A proper understanding of physical properties such as seismic attenuation in the crustal volume of this region can provide hints toward the possible source for the geodynamic events in the past and present that likely lead to the observed deformation. In order to model intrinsic and scattering attenuation separately, we perform a nonempirical coda-wave modeling approach in which a fitting process between observed and synthetic coda-wave envelopes is performed for each earthquake in multiple frequency bands. Here, the acoustic radiative transfer theory, assuming multiple isotropic scattering, was utilized for the forward modeling of the synthetic coda-wave envelopes of local earthquakes. Our findings generally highlight the prominent nature of intrinsic attenuation over scattering attenuation, implying the presence of thick volcanic rocks with relatively high attenuation values beneath Central Anatolia. Overall, the spatial distribution of the attenuation at varying frequencies marks the Kırşehir Massif distinctively with its considerable high-attenuating character. Our findings, combined with early seismological and geo-electrical models, suggest a possible partial melt beneath most of the Central Anatolian Volcanic Province, and the resultant zones of elevated fluid-rich content exhibit high and dominant intrinsic attenuation. To the southeast, a gradual decrease in the observed attenuation coincides with the Central Taurus Mountains where high altitude is considered to be evolved following the slab break-off and resulting mantle upwelling.

Geo-electrical prognosis of aquifer protectivity, corrosivity, and vulnerability via index-based models within a major coastal milieu

This work was carried out to investigate the protective capacity, vulnerability, and corrosivity within a major coastal milieu in Southern Nigeria with the use of index-based geo-electrical modeling methods. Vertical electrical soundings were undertaken at twenty locations with the aid of Schlumberger array having a maximum electrode spacing of 400 m. The results indicated that the lithology comprised four subsurface layers having variable values of resistivity and thickness. The Dar-Zarrouk parameter, the Aquifer Vulnerability Index (AVI), and the GOD (Groundwater occurrence G, Overlying lithology O and Depth to aquifer D) models were employed to appraise measures of aquifer protectivity and vulnerability to contamination. The longitudinal conductance values ranged from 0.0071–1.95 mhos with a mean of 0.32 mhos, indicating moderate protectivity. AVI values ranged from 1.73–4.10 with a mean of 3.03, indicating moderate aquifer vulnerability. The GOD indices ranged from 0.35–0.63 with a mean of 0.49, indicating moderate aquifer vulnerability. Corrosivity was also computed based on topsoil resistivity values which ranged from 12.7 to 664.2 Ωm with a mean of 168.17 Ωm, indicating moderate corrosivity, and demonstrating the unsuitability of corrosive locations for laying underground pipes. All the index-based models gave similar interpretations, indicating moderate aquifer protectivity and susceptibility. These results were corroborated by 2D electrical resistivity tomography surveys conducted at four stations. This work has therefore delineated important aquifer geo-hydraulic properties with index-based geo-electrical modeling techniques. The results obtained are critical for effective aquifer management, conservation, and sustainability.

Geoelectric Modelling Based on Bessel Functions Integral and Damped Least-Square Inversion for Layered Earth Models

Apparent resistivity formulations of a geoelectric method for layered Earth have been obtained analytically in the Bessel functions integral. This formulation is obtained through the application of two boundary conditions. First, the current density is zero at the Earth's surface. Second, current density and electric potential are continuous at the boundary between the layers. Geoelectric modelling can be done using several model parameters through this formulation. The modelling results show that apparent resistivity formulations based on Bessel function integrals can simulate apparent resistivity curves for layered earth models. The inversion process uses the apparent resistivity formulation based on Bessel function integrals. In this study, the inversion process uses a dumped Least-squared inversion method. The first step begins by testing the inversion program using synthetic data. After that, the inversion program is used in field data. The inversion results using field data from geoelectric data in Mount Pandan, East Java, show that the program was well done, with a minimal error value of 1.24%.

Application of Electrical Resistivity Sounding Method for Groundwater Exploration in Ugboshi-Afe, Akoko-Edo, Southwestern Nigeria

Schlumberger configuration of Vertical Electrical Sounding (VES) was deplored for groundwater exploration in Ugboshi-Afe study location. A total of nine (9) vertical electrical sounding (VES) locations, spread across the community were occupied to gain insights on the hydrogeological settings within the area. Quantitative computer interpretation of the data yielded three to six geoelectrical layers. The geoelectrical interpretation models were geologically interpreted as: clayey sand/sandy topsoil with resistivity values ranging from 54 Ohm-m to 424 Ohm-m, this is underlain by a relatively high resistivity layer (228-868 Ohm-m) interpreted as lateritic sand, and below this the low resistivity (21.38 Ohm-m – 42.4 Ohm-m) interpreted as clay, then the fractured/weathered geoelectric layer (95Ohm-m to 385Ohm-m), while the fresh basement has resistivity greater than 400 Ohm-m. Several maps such as fractured window thickness map, isoresistivity map of the fractured interval, overburden map were generated to investigate trends that may hold prospects for groundwater exploitation. Furthermore, hydro-resistivity parameters such as Total Transverse Resistance, Total Longitudinal Conductance, Resistivity Reflection Co-efficient and Resistivity Contrast values were calculated for further screening of the study site for groundwater development. The VES1, VES2, VES6 and VES7 locations were concluded as holding prospects for groundwater potentials in the study area.

Three-dimensional inversion for short-offset transient electromagnetic data based on 3D U-Net

Abstract The short-offset transient electromagnetic (SOTEM) method carries out surveys in the near source region; the strong signal makes it suitable for deep detection with high precision. When the underground structure is complex, three-dimensional (3D) inversion of SOTEM data is necessary to meet the need of high-precision detection. Currently, difficulties faced by the conventional 3D inversion methods include high computational complexity, and the influence of the initial model. Deep learning (DL), as a completely nonlinear algorithm, can predict the underground structure from the measured data. DL is completely data-driven and does not use traditional misfit optimization methods. In this study, an efficient method is proposed to conduct 3D inversion for SOTEM data, which trains a 3D U-Net based on massive data to establish a mapping from SOTEM data to geoelectric models. After the training is completed, we input the new SOTEM data into the trained network, and the corresponding geoelectric model can be obtained. Although the training is a time-consuming process, prediction for new data can be completed in seconds. The inversion results for simulated data indicate that the 3D U-Net has good generalization performance and anti-noise ability. The inversion performance of 3D U-Net on the double-anomaly model improves by 51.1% compared to a 3D fully convolutional network (FCN). The inversion results of the 3D U-Net on the field data successfully delineated the aquiferous collapse column. The inversion results for simulated and field data demonstrate that the proposed method can achieve accurate 3D inversion for large volume of data while greatly saving computational time.

GEOCRYOLOGICAL CONDITIONS OF THE FORMATION OF GIANT SPRING AUFEIS AT THE ANMANGYNDA RIVER (MAGADAN REGION) ACCORDING TO GEOPHYSICAL DATA

Giant aufeis fields, common in the Northeast of Russia, are the indicators of water exchange processes in cryosphere. The development of ideas about icing processes is relevant both from the fundamental point of view of studying the permafrost evolution, and from the practical point of view – for the development of aufeis hazard measures. The aufeis in the Anmangynda River basin (aufeis glade area 7 km2) is considered representative of the region, and its studies have been carried out since 1962. In 2022, during the period of maximum thawing of the active layer Electrical Resistivity Tomography (ERT) soundings were carried out at the aufeis glade aiming to identify underchannel taliks and flooded fault zones in bedrock, including local areas of groundwater discharge. It was found that within the main river channels there are underchannel taliks up to 30 m deep. According to the results of 2D inversion, local anomalies of low electrical resistivity mark groundwater filtration channels. In 3D geoelectrical models, pipe-like anomalies of low resistivity are identified in the areas of groundwater discharge, interpreted as filtration channels in the alluvium and the zone of exogenous fracturing in bedrock formed by sandy-clay shales, as well as linear vertical anomalies of low resistivity, interpreted as faults. On vertical sections of 3D resistive models, a connection between faults and filtration channels in alluvium and a layer of exogenous fracturing is traced. In the right bank of the valley, geoelectric signs of taliks in the bedrock, presumably associated with fault tectonics, have been established. It is assumed that the identified faults are the additional transit routes for groundwater in the Anmangynda River valley, along with the alluvial aquifer and the zone of exogenous fracturing of bedrock.