Allocation of Potentially Environmentally Hazardous Sections on Pipelines

: The paper presents an analysis of the method of recording the magnetic component of the Earth’s natural pulse electromagnetic ﬁeld in an urban environment. This method of recording has already proved itself to be a method that allocates the stressed sections of rock mass at mining and, therefore, the authors consider it effective for allocating active tectonic disturbances and forecasting accidents at underground utilities, which will help reduce the potential environmental hazard of these objects.


Introduction
Pipelines play a significant role in the life of a modern person. They are used for the transport of oil and gas, drain sewage, supplying water and other purposes. Nowadays, this type of transportation is actively developing, and its total extensiveness is growing. However, pipelines are also a cause of environmental hazards. The consequences of pipeline accidents are diverse, but the main ones are pollution of the environmental components and decrease of biological productivity.
In cities, the largest number of pipeline accidents occur in sewerage and heating systems. These accidents directly affect the life and health of the population. Accidents in sewer networks can lead to the ingress of wastewater into the groundwater, contamination of the area adjacent to the accident site, favorable conditions for the development of pathogenic bacteria, and, as a result, epidemics.
In St. Petersburg, the situation is complicated by the proximity of the Gulf of Finland and the need to ensure a positive environmental situation in the Baltic Sea in accordance with the requirements of the Helsinki Commission. One of the causes of pipeline accidents is the stress-strain state of rock mass, wherein the generation of pulsed electromagnetic radiation occurs. The immediate causes of accidents are movements along faults edges. The possibility of such movements is described by the "active" characteristic. The direction of movement of the fault edges can be either vertical, horizontal, or combined. Vertical movements break the pipe, and horizontal movements stretch or compress it. Both scenarios lead to the destruction of the pipeline material.
The method considered in this article is based on the fact that movement along the edges of faults creates a mechanoelectrical effect of accumulation of electric charges, and then relaxation of charges in the form of pulsed spark discharges that create pulsed electromagnetic radiation in the rock mass [1].
A geodynamic fault, in the general case, is represented in the form of tectonically stressed subsoil and land surface areas, having a tendency towards dynamics [2]. Within these areas, Geosciences 2021, 11, 3 2 of 11 the risk of natural or man-made transportation and environmental accidents increases [3][4][5][6]. Small-amplitude movements along faults can result in discontinuity in residential engineering structures, accidents in heat distribution networks and water discharge systems, gas and oil pipelines, as well as in tunnels [7]. Currently, the issue of geodynamic activity at mining facilities and related accidents at mine workings is a subject of intensive study [8][9][10][11].
Various scientific schools are engaged in researching active faults. Particular attention in this area is paid to predicting dangerous geodynamic processes in the mining industry and geodynamic risks. A group of Russian scientists have achieved significant results in this area [12][13][14][15], as well as members of the scientific school of China [16,17] and others. The obtained results can also be useful for engineering facilities that are not included in the field of mining [18,19].
In urban environments, underground utilities are particularly prone to being affected by active faults as they are linear-type facilities of significant length, with high chances of crossing underground water in the city territory fault areas. These facilities are designed with a lower margin of safety than residential buildings, industrial structures and underground tunnels. Accidents involving utilities result in significant costs of repair, restoration of the territory, modernization of the infrastructure, and, inter alia, cause significant inconvenience to the population. Hot water outbursts are hazardous to human life and health. Finally, they may negatively affect the operation of city transportation [20][21][22].
In addition to timely diagnostics of pipeline faults, it is necessary to conduct timely hydrogeological monitoring of underground water in cities since, if there are breaches in urban communications, ground water is subject to negative impacts. A study of the geoecological state of underground water in the city territory was carried out by a team of young scientists, whose research results were presented in a series of works [23][24][25][26]. The impact of microbial activity in the underground space including metro tunnels, utility networks, deep mines, tailing dumps, infrastructure with shallow and deep foundations was studied [27].
In 2018, more than 5000 accidents involving urban heat distribution networks occurred in St. Petersburg. In 2019, in St. Petersburg, there were five major accidents in the main central heating system, and, during the winter heating season, accidents on pipelines occur on a regular basis. In February 2019, hot water pipes broke in the Kalininsky and Frunzensky districts of the city. Over the repair period, which lasted more than a day, 677 and 32 buildings, in the two districts, including schools and healthcare facilities, were disconnected from the heating system. Healthcare facilities were heated via mobile heating systems, some houses were switched to a reserve heat supply network, some houses were left without hot water supply until the elimination of the consequences of the accident. In March 2019, 280 buildings were disconnected in the Kolpinsky district, which makes nearly a half of all residential buildings in the area. Moreover, as a rule, such accidents may have tragic consequences. In the spring of 2019, a hot water pipe broke in the Krasnoselsky district of the city, as a result of which one man received thermal burns of moderate severity, and one woman received lethal thermal burns. Moreover, more than 20,000 people were left without hot water supply. At the end of September 2018, in the same district, there was a rupture of the main pipeline with the formation of a failure. As a result of the accident, two people received lethal thermal burns, and damage was caused to the property of some citizens.
To reduce the number of accidents, municipal heat supply services perform replacement of old utilities and conduct timely diagnostics of highways [28][29][30][31]. The authors of this paper assume that the application of the method of measuring the magnetic component of the Earth's natural pulse electromagnetic field is expedient for allocating emergency sections induced by the stress-strain state of the soils surrounding the pipeline.
The study of tectonic disturbance activity began in the 1950s, but the concept of "active fault" was introduced in the 1970s [32][33][34][35][36]. At the same time, despite the similarity of the definitions of "active fault" proposed by scientists of multiple countries, a single definition still does not exist. According to [32], which made an attempt to generalize all available definitions, an active fault is a disjunctive tectonic disturbance of geological bodies in the relief or in subsoil, characterized by evidence of directional movement of blocks being separated by tectonic disturbance during the past hundreds of thousands of years at an average speed of at least hundredths of a mm per year. However, the disturbance activity duration is not a category that is convenient to apply practically, therefore it is expedient to proceed with definite quantitative values that can be recorded. The magnetic component (MC) of the Earth's natural pulse electromagnetic field (ENPEMF) can be considered to be one of such values. It is known that the Earth's electromagnetic field is a superposition of electromagnetic fields of different nature [37]. The resultant electromagnetic field is alternating and has a pulsed nature. The ENPEMF has the background values of the concentration of charges in rock mass, the magnitudes of strength of the electric E and magnetic H fields, magnetic induction B and magnetic flux Φ. The background value depends on the properties of the surrounding space-the structure of rock mass, dynamoelectric phenomena occurring in rock mass, microseismic motion, manmade load, etc. Further, the ENPEMF is characterized by the presence of pulsed signals from remote sources, such as atmospherics and deep lithospheric processes [10,38]. Such pulsed signals are characterized by equal peak value of the signal at each point; therefore, they are not recorded during measurements, but increase the background value. In addition, it is worth noting that the background value is also affected by various manmade factors. However, since the operating frequency of manmade sources (cable power lines, operating electromechanical installations, etc.) is 50 Hz and the sources generate a signal harmonically varying according to the sine law, this component can be easily excluded. The method of recording the ENPEMF magnetic component consists in measuring the magnitude of the EMF induced in the receiving antenna and falling within an alternating electromagnetic field.
Earthquake prediction was one of the first fields of application for the method of recording the ENPEMF MC. This is due to the fact that a high concentration of charges is not typical for a stable equilibrium of the system; therefore, it is an earthquake precursor [39,40]. Later this method was transformed and developed in this area. For instance, the VAN method is known. The method is based on the measurement of low frequency electrical signals (SES). According to this method, SES are recorded and evaluated in different stations [41,42].
To adapt the EMF method to the conditions of underground mining, to determine the sufficient sensitivity of the method for recording seismic events of much lower energy than an earthquake, a joint research with other noncontact and contact geophysical methods was implemented. Along with others, the core disking method was widely used as the most studied one [43]. Further, specialists of the Mining University conducted work on the comparison of EMF data and acoustic emission [44,45].
According to the results of multiyear research aimed at predicting a rock-bump hazard at ore and nonore deposits, it was revealed that the indicator of the modern geodynamic activity of faults is the variability of the ENPEMF characteristics in space and time compared to the adjacent areas. Therefore, the team of authors proposed to adapt the method of measuring the ENPEMF MC for the urban environment. For this purpose, a number of streets have been selected along which underground utilities are laid and where accidents have been recently registered, namely, pipeline breaks, ground collapses, and associated liquid spills.

Materials and Methods
To analyze the efficiency of the ENPEMF MC measurements, carried out to identify active faults and detect hazardous areas, the following actions were carried out: However, it is worth noting that not all faults plotted on the map [46] are active and pose a threat to engineering facilities, subsequently, a simple comparison between accidents and the map will not ensure a reliable result. Moreover, information on accidents at underground utilities is not structured to date. Therefore, this study involves the analysis of both scientific and technical sources as well as mass media news feeds.
To conduct the study, a receiving antenna in the form of an inductance coil with a constant ferromagnetic core was used. Ferromagnetic materials were used to expand the measurement limits of the recorded signals since the induced EMF ranges from several µV to tens of mV. An alternating voltage voltmeter served as a recording device.
The method of recording the ENPEMF MC included recording the magnitude of the EMF induced in the receiving antenna falling within the alternating electromagnetic field in accordance with the law of electromagnetic induction. A rod of the antenna had the following characteristics: core material-M400HH, according to the Russian material classification; initial magnetic permeability-400 H/m; rod diameter-10 mm; rod length-200 mm; winding-180 mm; the diameter of the copper wire was 0.1 mm; the number of turns was about 3000. Constant sensitivity was confirmed in the frequency range from 1 to 18 kHz. The APPA-109N device was used as a recorder. Measurements in three directions of the antenna were carried out along the streets with a 5 m measurement pitch. The holding time at each pitch was 10 s.
Subsequently, the studied areas were analyzed geologically. The upper layers of a typical geological record in this area of the city are Quaternary deposits, which are interbedded sandy-clay soils. The thickness of the layer of Quaternary deposits is not consistent and varies from 30 to 40 m in the studied area. There are peat interlayers in the upper part of the Quaternary deposits. Quaternary sediments are everywhere covered by technogenic soils from above, the thickness of which is less than 5 m. The depth of the groundwater level is 2 m. The bedding rocks are the Upper Vendian Verkhnekotlinsky clays. They are dense, dry, hard clays, with interbedded sandstone. The upper part is more humid; the lower part is drier and where the sandstones may appear. The crystalline basement lies at a depth of about 200 m and consists of granites, quartzites, gneisses [47,48].
Tectonic faults occur in the sedimentary cover rocks and are prolonged by basement faults. The average width of the investigated faults is 80-100 m. Tectonic faults have no outlet in the series of the Quaternary deposits and are traced only by indirect indicators. Thus, for this study, it is necessary to take into account that the signal is generated in bedrocks and basement rocks and passes through the series of loose sediments.
To establish the validity of the study, it is necessary to determine the possibility to register on the surface a signal from a source located at the average depth of the base of the Quaternary sediments (about 35 m). To that end, it is necessary to take into account the values of the electrical and magnetic parameters of soils. Table 1 shows the averaged values of the electrical and magnetic parameters of soils typical for the investigated area for a frequency of 1 kHz.
As can be seen from the table, soil parameters cannot obstruct the passage of a pulsed electromagnetic field generated by active tectonic faults. The pulse frequency varies from 10 3 to 10 5 Hz, so the length of the electromagnetic wave will be hundreds of times longer than the distance to the receiving device.
According to Faraday's law, the magnitude of the EMF occurring in a coil with a constant core is equal to the rate of change in magnetic flux through a surface limited by the coil loop: The magnetic flux is proportional to the magnitude of the magnetic induction vector B, the area of the limited loop S and the cosine of the angle between the magnetic induction vector B and the normal to the loop plane: As is known, the magnetic induction vector B is directly proportional to the magneticfield vector H: where µ 0 is the absolute magnetic permeability of the vacuum, µ is the relative magnetic permeability of the medium. According to the dependences depicted above, the magnitude of the EMF E, typical for the studied section with a length of several hundred meters, will also have a background value since the nature of the change in magnetic flux will be constant in the absence of geodynamic activity. Nonetheless, if stress-strain states occur in rock mass, the recorded signal acquires a step-like variable shape. The amplitudes of such anomalous values of E can exceed the background value by many times and have both positive and negative values.
The occurrence of anomalous values is due to the occurrence of an elastic wave formed under the action of deforming forces that arise with the change in mechanical and hydrogeological properties, the stress state of the soil massif, which also comes with increased intensity of corrosion activity. The deforming forces cause the redistribution of charges at the edges of geological faults, and, as is known, the redistribution of charges is accompanied by their movement, which results in induction of short-term electromagnetic fields of greater intensity.
The value of EMF before the occurrence of stress-strain state is a variable magnitude of a pulsed nature, which has a background value of E F (t). The occurrence of stress-strain state in massif results in a change in E F by the value of ∆E. Therefore, the induced EMF E F in the disturbed massif will be equal to: The value of ∆E is the anomaly amplitude and can take both positive and negative values. Therefore, it can increase or decrease the background value of E F . This is explained by the principle of superposition. The increase in the background value is due to the enhancement of electromagnetic emission, and, thereby, the electromagnetic field enhancement, i.e., increase in vectors of electric and magnetic fields, takes place. The enhancement is due to the accordant direction of the background field and the field induced by stress-strain state.
At a negative value of the amplitude, the pulse electromagnetic field, occurring under deformation, is directed opposite to the background field.
Thus, according to this method, not only the excess over the background value, but also the minimum values of the EMF, should be taken into account.

Results
The averaged value of amplitudes that exceed significance levels is deemed to be an informative indicator of electromagnetic radiation of a geological nature. In this study, the criterion for allocating abnormal areas is a change in the indicator by more than 20%, regardless of its sign.
Profiling was performed under the mode of discrete measurements at a pitch of 20 m along the surface. The exposure time was 10 s at a point. Signals were recorded at three orientations of the antenna magnetic moment: one vertical and two horizontal-forward along the profile and across the profile. Figure 1 illustrates some of the results of the study. It can be seen from the figures that discontinuous faults are clearly recorded on the profiles of the ENPEMF MC, wherein most accidents are confined to the areas of these faults. In this figure, the diagram of the vertical component A Z is taken as representative. As can be seen from the diagram, in the measurement areas from 0 to 30 m, from 460 to 500 m, from 790 to 820 m and from 1080 to 1110 m, the parameters of ∆E exceed the background value of E F by more than 20%. These areas coincide with the sites of accidents at underground utilities.
Geosciences 2021, 11, x FOR PEER REVIEW 6 of 11 the enhancement of electromagnetic emission, and, thereby, the electromagnetic field enhancement, i.e., increase in vectors of electric and magnetic fields, takes place. The enhancement is due to the accordant direction of the background field and the field induced by stress-strain state. At a negative value of the amplitude, the pulse electromagnetic field, occurring under deformation, is directed opposite to the background field. Thus, according to this method, not only the excess over the background value, but also the minimum values of the EMF, should be taken into account.

Results
The averaged value of amplitudes that exceed significance levels is deemed to be an informative indicator of electromagnetic radiation of a geological nature. In this study, the criterion for allocating abnormal areas is a change in the indicator by more than 20%, regardless of its sign.
Profiling was performed under the mode of discrete measurements at a pitch of 20 m along the surface. The exposure time was 10 s at a point. Signals were recorded at three orientations of the antenna magnetic moment: one vertical and two horizontal-forward along the profile and across the profile. Figure 1 illustrates some of the results of the study. It can be seen from the figures that discontinuous faults are clearly recorded on the profiles of the ENPEMF MC, wherein most accidents are confined to the areas of these faults. In this figure, the diagram of the vertical component AZ is taken as representative. As can be seen from the diagram, in the measurement areas from 0 to 30 m, from 460 to 500 m, from 790 to 820 m and from 1080 to 1110 m, the parameters of Δ exceed the background value of by more than 20%. These areas coincide with the sites of accidents at underground utilities. The results of the study, presented in Figure 2, also demonstrate the convergence of the location of the anomalies by AZ and the sites of accidents on the pipelines. However, at 450 m, a fluctuation in the AEMF MС curve, which coincides with the place of the pipeline break but does not coincide with any disturbance shown on the map, is observed [30]. This is probably due to the fact that either the disturbance actually exists but is not plotted on the map, and, therefore, adjustments are required, or this disturbance was not active during the map compilation period but during the research period the movement began again. It is also possible that there is a low-ranked disturbance in this section, wherein The results of the study, presented in Figure 2, also demonstrate the convergence of the location of the anomalies by AZ and the sites of accidents on the pipelines. However, at 450 m, a fluctuation in the AEMF MC curve, which coincides with the place of the pipeline break but does not coincide with any disturbance shown on the map, is observed [30]. This is probably due to the fact that either the disturbance actually exists but is not plotted on the map, and, therefore, adjustments are required, or this disturbance was not active during the map compilation period but during the research period the movement began again. It is also possible that there is a low-ranked disturbance in this section, wherein high voltage is concentrated, and at the time of the massif unloading (discontinuous fault expansion), an accident occurred.        Designations, used in Figures 1-3, are given in Table 2.  Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamica amounted to about 20% of the profile total length. Meanwhile, t 51 known accidents coincided with the zones affected by GAF, w ratio of the number of accidents per running meter within zones accidents beyond their bounds made 17.3. The data obtained ind accidents at underground utilities is much greater within the zon beyond their bounds, however, to conclude more precisely, it is n ther collection of statistical data on accidents. Table 3 summarize Table 3. Designations to figures.

Characteristics
The total length of the electrical survey profile, with recording magnetic component, in km The total length of the zones affected by GAF, allocated by the logical data, in running km The total length of internal fault-blocks in running k The total number of accidents at underground utilities of the w system The number of accidents at underground utilities of the water s within the zones affected by GAF The number of accidents at underground utilities of the water s beyond the bounds of GAF The number of accidents at underground utilities of the water s within the zones affected by GAF per running km The number of accidents at the underground utilities of the w system beyond the bounds of GAF per running km The ratio of the number of accidents within the zones affected b number of accidents beyond their bounds  Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamically active faults (GA amounted to about 20% of the profile total length. Meanwhile, the locations of 41 out 51 known accidents coincided with the zones affected by GAF, which also means that ratio of the number of accidents per running meter within zones of GAF to the number accidents beyond their bounds made 17.3. The data obtained indicate that the number accidents at underground utilities is much greater within the zones affected by GAF th beyond their bounds, however, to conclude more precisely, it is necessary to perform f ther collection of statistical data on accidents. Table 3 summarizes the results of the stud Table 3. Designations to figures.

Characteristics
Valu The total length of the electrical survey profile, with recording of the AEMF magnetic component, in km 6.38 The total length of the zones affected by GAF, allocated by the sum of geological data, in running km 1.22 The total length of internal fault-blocks in running km 5.16 The total number of accidents at underground utilities of the water supply system 51 The number of accidents at underground utilities of the water supply system within the zones affected by GAF 41 The number of accidents at underground utilities of the water supply system beyond the bounds of GAF 10 The number of accidents at underground utilities of the water supply system within the zones affected by GAF per running km 33.61 The number of accidents at the underground utilities of the water supply system beyond the bounds of GAF per running km 1.94 The ratio of the number of accidents within the zones affected by GAF to the number of accidents beyond their bounds 17.34

Zones of active faults
Geosciences 2021, 11, x FOR PEER REVIEW Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamica amounted to about 20% of the profile total length. Meanwhile, t 51 known accidents coincided with the zones affected by GAF, w ratio of the number of accidents per running meter within zones accidents beyond their bounds made 17.3. The data obtained ind accidents at underground utilities is much greater within the zon beyond their bounds, however, to conclude more precisely, it is n ther collection of statistical data on accidents. Table 3 summarize Table 3. Designations to figures.

Characteristics
The total length of the electrical survey profile, with recording magnetic component, in km The total length of the zones affected by GAF, allocated by the logical data, in running km The total length of internal fault-blocks in running k The total number of accidents at underground utilities of the w system The number of accidents at underground utilities of the water s within the zones affected by GAF The number of accidents at underground utilities of the water s beyond the bounds of GAF The number of accidents at underground utilities of the water s within the zones affected by GAF per running km The number of accidents at the underground utilities of the w system beyond the bounds of GAF per running km The ratio of the number of accidents within the zones affected b number of accidents beyond their bounds  Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamically active faults (GA amounted to about 20% of the profile total length. Meanwhile, the locations of 41 out 51 known accidents coincided with the zones affected by GAF, which also means that ratio of the number of accidents per running meter within zones of GAF to the number accidents beyond their bounds made 17.3. The data obtained indicate that the number accidents at underground utilities is much greater within the zones affected by GAF th beyond their bounds, however, to conclude more precisely, it is necessary to perform f ther collection of statistical data on accidents. Table 3 summarizes the results of the stud Table 3. Designations to figures.

Characteristics
Valu The total length of the electrical survey profile, with recording of the AEMF magnetic component, in km 6.38 The total length of the zones affected by GAF, allocated by the sum of geological data, in running km 1.22 The total length of internal fault-blocks in running km 5.16 The total number of accidents at underground utilities of the water supply system 51 The number of accidents at underground utilities of the water supply system within the zones affected by GAF 41 The number of accidents at underground utilities of the water supply system beyond the bounds of GAF 10 The number of accidents at underground utilities of the water supply system within the zones affected by GAF per running km 33.61 The number of accidents at the underground utilities of the water supply system beyond the bounds of GAF per running km 1.94 The ratio of the number of accidents within the zones affected by GAF to the number of accidents beyond their bounds 17.34 House numbers Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamica amounted to about 20% of the profile total length. Meanwhile, t 51 known accidents coincided with the zones affected by GAF, w ratio of the number of accidents per running meter within zones accidents beyond their bounds made 17.3. The data obtained ind accidents at underground utilities is much greater within the zon beyond their bounds, however, to conclude more precisely, it is n ther collection of statistical data on accidents. Table 3 summarizes Table 3. Designations to figures.

Characteristics
The total length of the electrical survey profile, with recording magnetic component, in km The total length of the zones affected by GAF, allocated by the logical data, in running km The total length of internal fault-blocks in running k The total number of accidents at underground utilities of the w system The number of accidents at underground utilities of the water s within the zones affected by GAF The number of accidents at underground utilities of the water s beyond the bounds of GAF The number of accidents at underground utilities of the water s within the zones affected by GAF per running km The number of accidents at the underground utilities of the w system beyond the bounds of GAF per running km The ratio of the number of accidents within the zones affected b number of accidents beyond their bounds  Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamically active faults (GA amounted to about 20% of the profile total length. Meanwhile, the locations of 41 out 51 known accidents coincided with the zones affected by GAF, which also means that t ratio of the number of accidents per running meter within zones of GAF to the number accidents beyond their bounds made 17.3. The data obtained indicate that the number accidents at underground utilities is much greater within the zones affected by GAF th beyond their bounds, however, to conclude more precisely, it is necessary to perform f ther collection of statistical data on accidents. Table 3 summarizes the results of the stud Table 3. Designations to figures.

Characteristics
Valu The total length of the electrical survey profile, with recording of the AEMF magnetic component, in km 6.38 The total length of the zones affected by GAF, allocated by the sum of geological data, in running km 1.22 The total length of internal fault-blocks in running km 5.16 The total number of accidents at underground utilities of the water supply system 51 The number of accidents at underground utilities of the water supply system within the zones affected by GAF 41 The number of accidents at underground utilities of the water supply system beyond the bounds of GAF 10 The number of accidents at underground utilities of the water supply system within the zones affected by GAF per running km 33.61 The number of accidents at the underground utilities of the water supply system beyond the bounds of GAF per running km 1.94 The ratio of the number of accidents within the zones affected by GAF to the number of accidents beyond their bounds 17.34 Accident sites at underground utilities Designations, used in Figures 1-3, are given in Table 2. Thus, the length of the zones affected by geodynamically active faults (GA amounted to about 20% of the profile total length. Meanwhile, the locations of 41 ou 51 known accidents coincided with the zones affected by GAF, which also means that ratio of the number of accidents per running meter within zones of GAF to the numbe accidents beyond their bounds made 17.3. The data obtained indicate that the numbe accidents at underground utilities is much greater within the zones affected by GAF th beyond their bounds, however, to conclude more precisely, it is necessary to perform f ther collection of statistical data on accidents. Table 3 summarizes the results of the stu Table 3. Designations to figures.

Characteristics
Valu The total length of the electrical survey profile, with recording of the AEMF magnetic component, in km 6.38 The total length of the zones affected by GAF, allocated by the sum of geological data, in running km 1.22 The total length of internal fault-blocks in running km 5.16 The total number of accidents at underground utilities of the water supply system 51 The number of accidents at underground utilities of the water supply system within the zones affected by GAF 41 The number of accidents at underground utilities of the water supply system beyond the bounds of GAF 10 The number of accidents at underground utilities of the water supply system within the zones affected by GAF per running km 33.6 The number of accidents at the underground utilities of the water supply system beyond the bounds of GAF per running km 1.94 The ratio of the number of accidents within the zones affected by GAF to the number of accidents beyond their bounds 17.3

Anomalous values of Az
Thus, the length of the zones affected by geodynamically active faults (GAF) amounted to about 20% of the profile total length. Meanwhile, the locations of 41 out of 51 known accidents coincided with the zones affected by GAF, which also means that the ratio of the number of accidents per running meter within zones of GAF to the number of accidents beyond their bounds made 17.3. The data obtained indicate that the number of accidents at underground utilities is much greater within the zones affected by GAF than beyond their bounds, however, to conclude more precisely, it is necessary to perform further collection of statistical data on accidents. Table 3 summarizes the results of the study. Table 3. Designations to figures.

Characteristics Value
The total length of the electrical survey profile, with recording of the AEMF magnetic component, in km 6.38 The total length of the zones affected by GAF, allocated by the sum of geological data, in running km 1.22 The total length of internal fault-blocks in running km 5.16 The total number of accidents at underground utilities of the water supply system 51 The number of accidents at underground utilities of the water supply system within the zones affected by GAF 41 The number of accidents at underground utilities of the water supply system beyond the bounds of GAF 10 The number of accidents at underground utilities of the water supply system within the zones affected by GAF per running km 33.61 The number of accidents at the underground utilities of the water supply system beyond the bounds of GAF per running km 1.94 The ratio of the number of accidents within the zones affected by GAF to the number of accidents beyond their bounds 17.34

Discussion
The method of recording the ENPEMF MC in an urban environment has shown good convergence of the recorded anomalies and known active faults, to which accidents at urban utilities, in particular in lines of water supply and water discharge systems, are confined. Nonetheless, not all active faults plotted on the map [46], confirming the geodynamically active pattern according to data of morphometric and lineament analyzes, are recorded within the field, which is due to the high ratio of the thickness of tectonic disturbances and the overlying strata-sediments. In addition, the result obtained indirectly indicates a variable degree of activity, namely, the concentration of natural stress fields within separate geodynamic blocks, bound by the disturbances under consideration. It is worth noting that the complexity associated with allocating low-ranked disturbances, with a high probability, is not substantiated in terms of the method accuracy with respect to the allocation of potential emergency sections on pipelines since, as a rule, in diagrams, such disturbances are presented in the form of a single fluctuation zone, and to implement the most reliable allocation thereof, the existing measurement techniques should be substantially altered.

Conclusions
Hereby, the field studies, which were conducted on more than 6 km of streets in an urban environment, have shown a high efficiency of the method for identifying potentially dangerous areas on underground utilities. The obtained data also demonstrate that the vertical component Az is the most significant one, since it has the highest convergence of data, and it is also the least susceptible to manmade interference.
Eighty percent of accidents which took place in the considered areas occurred in the GAF zones, and all the accidents correlate in the graphs with fluctuations of more than 20%. Two accidents that are not spatially related to the known GAF zones are reflected in the graphs by anomalies of a smaller order.
It is also worth noting that the main advantage of the method in comparison with other geophysical methods used for solving such problems is the ability to perform measurements at any time of the year with few staff and high productivity and also that there is no need for grounding.
The authors consider the possibility of further developing the studied method to implement a monitoring system at underground utilities of water supply and water discharge systems based on the tested instrumental base. The use of this method at the stage of predesign and design surveys in the future, with further consideration of hazardous zones allocation during installation of the utility system facilities, will make it possible to reduce the accident rate at underground utilities of water supply and water discharge systems. Moreover, the use of this method in combination with operating monitoring systems will significantly reduce the time of allocating an emergency section, and allow the determining of the causes of its occurrence. This will reduce the time for performing repair work and allow entering corrections in the design solution with respect to this zone due to applying compensatory measures. As a result, an increase of living standards of the population in the urban environment and a decline if the negative impacts on the environment, the likelihood of injury and the death of people and animals are expected.