A method and system for rapid positioning of leakage channels in dam core walls over long distances

Abstract

The present invention relates to the technical field of locating leakage channels in the core wall of a dam, and specifically to a method and system for rapid positioning of leakage channels in the core wall of a dam by long-distance detection. The present invention sets up a physical detection line in the core wall of the dam, and regularly collects detection parameters of the physical detection line; when the detection parameters are abnormal, the initial abnormal time t0 and the initial leakage infiltration area are recorded; secondly, the present invention obtains the corresponding detected leakage infiltration area according to the detection parameters of the physical detection line collected after t0; thirdly, the position of the center point of the leakage area is inferred according to the change of the leakage infiltration area, and the aperture of the leakage area is obtained according to the leakage rate; finally, the present invention obtains the leakage level according to the position of the center point of the leakage area and the aperture of the leakage area, and issues an early warning. The present invention can accurately obtain the position of the center point of the leakage area and the aperture of the leakage area through this method, and then accurately evaluate the leakage level.

Description

Technical Field

[0001] The invention relates to the technical field of locating leakage channels of dam core walls, and in particular to a method and system for quickly locating leakage channels of dam core walls by long-range detection. Background Art

[0002] A dam is usually composed of a multi-layer structure, including an impermeable body and a core wall on the upstream and downstream slopes. The core wall of a dam refers to a main impermeable body located inside the main structure of a water conservancy project, usually in the center of the entire dam cross section. The core wall of a dam is a usually thicker concrete or impermeable soil engineering layer, which is used to prevent water from flowing through the main structure of the dam, thereby ensuring the impermeability and structural stability of the dam.

[0003] Damage to the core wall structure of the dam causes water to flow through the core wall, forming a leakage channel, which affects the anti-seepage performance and stability of the dam. Therefore, it is necessary to locate the leakage channel of the core wall of the dam in time and deal with the leakage problem.

[0004] Most of the existing studies on core wall leakage evaluate the seepage safety of the core wall from the overall state of the core wall. For example, the Chinese patent with publication number CN115935710B proposes a method for calculating the permeability coefficient of the core wall material of the core wall rockfill dam and evaluating the seepage safety. This method conducts on-site sampling and random inspection of the core wall rockfill dam to be evaluated, and measures the saturated permeability coefficient of a small number of random inspection point samples, and calculates the mixing coefficient; then, based on the mixing coefficient and the calculation method of the permeability coefficient of the gravel clay core wall material, the permeability coefficient of other random inspection point samples whose saturated permeability coefficient has not been measured is calculated, and then the seepage safety of the core wall is evaluated. For example, the literature ("Research on Leakage Detection Method of Core Wall Materials of Different Densities in Earth-Rock Dams", "Water Conservancy Technical Supervision", Li Zhou et al., 2022-11-08) conducted leakage detection tests on core wall materials of different densities based on thermal pulse technology and time domain reflection technology to improve the operating safety and stability of earth-rock dams, and determined the leakage performance of sandy loam with different densities. The above method detects the overall state of the core wall and cannot accurately detect the leakage channel of the core wall. In the prior art, there is also a method of detecting the distribution of the leakage area of ​​the core wall by setting up a physical detection line. For example, the Chinese patent with the publication number CN116907753B proposes a rapid positioning method and system for detecting the leakage channel of the core wall of an earth-rock dam. This method sets a corresponding physical detection line layout scheme according to the structural characteristics of the dam type of different types of core wall earth-rock dams; establishes simulation models of different types of core wall earth-rock dams with preset leakage channels, and obtains the distribution form and characteristics of the model leakage infiltration area through numerical simulation; measures the physical property parameters of the dam material of the core wall of the earth-rock dam to be located, and identifies the distribution form and characteristics of the leakage infiltration area through inversion of the measured results of the dam material physical property parameters; compares the measured distribution form and characteristics with the corresponding dam type numerical simulation results for the degree of consistency, and the preset leakage channel position corresponding to the sample data with the consistent comparison is the actual engineering leakage point. However, when the area of ​​the leakage infiltration area is too large, the leakage center point cannot be accurately located, thereby reducing the accuracy of the evaluation of the leakage channel of the dam core wall.

[0005] Therefore, a method and system for long-distance detection and rapid positioning of leakage channels in dam core walls are proposed. Summary of the invention

[0006] The purpose of the present invention is to provide a method and system for rapid positioning of leakage channels in the core wall of a dam by setting up a physical detection line in the core wall of the dam and regularly collecting detection parameters of the physical detection line; when the detection parameters are abnormal, the initial abnormal time is recorded. and the initial leakage infiltration area; secondly, the present invention obtains the corresponding detected leakage infiltration area according to the detection parameters of the object detection line collected each time; thirdly, the position of the center point of the leakage area is inferred according to the change of the leakage infiltration area, and the aperture of the leakage area is obtained according to the leakage rate; finally, the present invention obtains the leakage level according to the center point position of the leakage area and the aperture of the leakage area, and issues an early warning.

[0007] To achieve the above object, the present invention provides the following technical solutions:

[0008] A method for rapid positioning of leakage channels in a dam core wall by long-distance detection, comprising:

[0009] Establish a physical detection line on the core wall of the dam and regularly collect detection parameters of the physical detection line; the detection parameters include resistivity and magnetic flux; when the detection parameters are abnormal, record the initial abnormal time and initial leakage infiltration area;

[0010] from Start, every The detection parameters are collected once at a time, and the detection parameters are collected N times in total;

[0011] Acquire the corresponding leakage infiltration area according to the detection parameters of the object detection line collected each time;

[0012] The central point position of the leakage area and the aperture of the leakage area of ​​the dam core wall are obtained according to the N detected leakage infiltration areas and the initial leakage infiltration areas. The step of obtaining the aperture of the leakage area includes:

[0013] Obtaining the leakage parameters of the core wall of the dam; the leakage parameters include the leakage value of the core wall material, the leakage rate and the leakage pressure difference; the leakage value of the core wall material is the permeability coefficient of the core wall; the leakage pressure difference is expressed as the absolute value of the water pressure difference on both sides of the center point of the leakage area;

[0014] The calculation formula of the pore size of the leakage area is:

[0015] ;

[0016] Among them, SLCC represents the pore size of the leakage area; Expressed as leakage rate; It is expressed as the leakage value of the core wall material; Expressed as leakage pressure difference; Expressed as leakage rate coefficient; Expressed as the core wall material influence coefficient; Expressed as leakage pressure difference influence coefficient; It is expressed as an exponential function with a natural constant as base;

[0017] The leakage level is obtained according to the central point position of the leakage area and the aperture of the leakage area, and an early warning is issued.

[0018] Furthermore, the process of determining whether the detection parameters are abnormal includes: performing spectral analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal.

[0019] Furthermore, the step of obtaining the central point position of the leakage area includes:

[0020] A wetting and spreading coordinate system is constructed based on the initial leakage and infiltration area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and infiltration area; the plane formed by the X-axis and the Y-axis of the wetting and spreading coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the wetting and spreading coordinate system are perpendicular to each other;

[0021] Selecting M infiltration and propagation directions, obtaining the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generating M infiltration boundary point coordinate sets;

[0022] Each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration propagation direction;

[0023] The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; by analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the initial leakage infiltration area at time;

[0024] Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period;

[0025] Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period;

[0026] According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period;

[0027] Obtain the leakage infiltration area to be tested at K time points according to the coordinates of the infiltration boundary points in the M infiltration propagation directions at each time point; when the overlap between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap, the leakage infiltration area to be tested is eliminated;

[0028] The areas of the remaining leakage infiltration regions to be detected are detected, and the center point position of the leakage infiltration region to be detected with the smallest area is taken as the center point position of the leakage region.

[0029] Furthermore, the step of obtaining the infiltration boundary point is:

[0030] Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area;

[0031] The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

[0032] Furthermore, the leakage rate acquisition process includes:

[0033] The leakage infiltration area is obtained from the initial abnormal time according to the M infiltration boundary point coordinate sets. arrive In the time, every The unit leakage rate of time; the calculation formula of the unit leakage rate is:

[0034] ;

[0035] in, Expressed as The unit leakage rate of the leakage infiltration area at the moment; the value range of i is 1 to N; Expressed as The distance between the infiltration boundary point and the origin corresponding to the j-th infiltration propagation direction in the detected leakage infiltration area at the moment; the leakage rate is the average value of N unit leakage rates.

[0036] Furthermore, the process of obtaining the leakage level includes:

[0037] Obtain the straight-line distance L between the center point of the leakage area and the center point of the bottom of the dam core wall; at the same time, detect the angle between the line connecting the center point of the leakage area and the center point of the bottom of the dam core wall and the ground plane ; The angle is less than 90 degrees;

[0038] Calculate the leakage risk value using the following formula:

[0039] ;

[0040] in, Expressed as leakage risk value; Expressed as core wall width; It is expressed as the core wall height; Expressed as leakage aperture parameter; It is expressed as the leakage lateral parameter; Expressed as leakage height parameter; It is expressed as an exponential function with a natural constant as the base; the leakage risk level is obtained according to the leakage risk value, and an early warning is issued.

[0041] A system for quickly locating leakage channels in dam core walls over long distances, comprising:

[0042] The first detection parameter acquisition unit is used to regularly collect the detection parameters of the object detection line on the core wall of the dam;

[0043] Abnormality judgment unit: used to judge whether the detection parameters are abnormal, and record the initial abnormal time when the detection parameters are abnormal and initial leakage infiltration area;

[0044] The second detection parameter acquisition unit is used to Start, every The detection parameters are collected once at a time, and the detection parameters are collected N times in total;

[0045] A detection leakage infiltration area acquisition unit is used to acquire a corresponding detection leakage infiltration area according to the detection parameters of the object detection line collected each time;

[0046] A leakage information acquisition unit: used for acquiring the central point position of the leakage area and the aperture of the leakage area of ​​the dam core wall according to the N detected leakage infiltration areas and the initial leakage infiltration area;

[0047] Leakage level acquisition unit: used to acquire the leakage level according to the center point position of the leakage area and the aperture of the leakage area, and to issue an early warning.

[0048] Furthermore, the process of determining whether the detection parameters are abnormal includes: performing spectral analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal.

[0049] Furthermore, the step of obtaining the central point position of the leakage area includes:

[0050] A wetting and spreading coordinate system is constructed based on the initial leakage and wetting area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and wetting area;

[0051] The plane formed by the X-axis and the Y-axis of the infiltration and propagation coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the infiltration and propagation coordinate system are perpendicular to each other;

[0052] Selecting M infiltration and propagation directions, obtaining the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generating M infiltration boundary point coordinate sets;

[0053] Each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration propagation direction;

[0054] The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; by analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the initial leakage infiltration area at time;

[0055] Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period;

[0056] According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period;

[0057] Obtain the leakage infiltration area to be tested at K time points according to the coordinates of the infiltration boundary points in the M infiltration propagation directions at each time point; when the overlap between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap, the leakage infiltration area to be tested is eliminated;

[0058] The areas of the remaining leakage infiltration regions to be detected are detected, and the center point position of the leakage infiltration region to be detected with the smallest area is taken as the center point position of the leakage region.

[0059] Furthermore, the step of obtaining the infiltration boundary point is:

[0060] Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area;

[0061] The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

[0062] Compared with the prior art, the present invention has the following beneficial effects:

[0063] 1. When the present invention detects that there is an abnormality in the detection parameters, the initial abnormal time and the initial leakage infiltration area are recorded, and the detected leakage infiltration area is obtained at a set time interval. The present invention reversely deduces the leakage infiltration area to be tested before the initial abnormal time based on the changes in the boundary points of the detected leakage infiltration area at different times, and selects the center point position of the leakage infiltration area to be tested with the smallest area as the center point position of the leakage area. This method can reduce the area of ​​the leakage infiltration area, more accurately obtain the center point position of the leakage area, and thus improve the accuracy of the evaluation of the leakage channel of the dam core wall.

[0064] 2. The present invention obtains the aperture of the leakage area based on the leakage rate, core wall material and the water pressure difference on both sides of the core wall to detect the aperture of the leakage channel. This method can accurately calculate the aperture of the leakage channel according to the intensity and extension of the leakage, thereby improving the accuracy of the leakage channel assessment of the dam core wall.

[0065] 3. The present invention obtains the leakage risk value based on the distance between the center point of the leakage area and the center point of the bottom of the core wall and the leakage aperture. This method combines the location and size of the leakage point. The location information reflects the specific location of the leakage on the core wall, while the size reflects the scale and extension of the leakage. Taking these two factors into consideration, the actual risk of leakage to the core wall of the dam can be more accurately assessed, so as to formulate targeted management and response strategies. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] Figure 1 Schematic diagram of dam structure provided by prior art;

[0067] Figure 2 A flow chart of a method for rapid positioning of a leakage channel in a dam core wall over a long distance provided by an embodiment of the present invention;

[0068] Figure 3 A structural diagram of a system for long-distance detection and rapid positioning of leakage channels in a dam core wall provided by an embodiment of the present invention;

[0069] Figure 4 A schematic diagram of the initial leakage infiltration area of ​​the dam core wall numbered 2 provided in the first embodiment of the present invention;

[0070] Figure 5 A schematic diagram of the initial leakage infiltration area of ​​the dam core wall numbered 3 provided in Example 1 of the present invention.

[0071] In the figure: 1. Downstream dam body; 2. Dam core wall; 3. Upstream dam body. DETAILED DESCRIPTION

[0072] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0073] A dam usually consists of a multi-layer structure, including the upstream and downstream anti-seepage bodies, the core wall, and possible auxiliary anti-seepage measures. The core wall of a dam refers to a main anti-seepage body located inside the main structure of a water conservancy project. Figure 1 , Figure 1 It includes the downstream dam body 1, the dam core wall 2 and the upstream dam body 3; the dam core wall 2 is usually located in the center of the entire dam cross section. The dam core wall 2 is a usually thick concrete or impermeable soil engineering layer, which is used to prevent water from flowing through the main structure of the dam, thereby ensuring the impermeability and structural stability of the dam.

[0074] Damage to the structure of the dam core wall 2 causes water to flow through the core wall, forming a leakage channel, thereby affecting the anti-seepage performance and stability of the dam. Therefore, it is necessary to locate the leakage channel in the dam core wall 2 in a timely manner and deal with the leakage problem.

[0075] The present invention provides a method for rapid positioning of leakage channels in the core wall of a dam over long distance detection. The method is applied to a rapid positioning system for rapid positioning of leakage channels in the core wall of a dam. For a specific method flow chart and system structure chart, refer to Figure 2 and Figure 3 .

[0076] Embodiment 1

[0077] Reference Figure 2 S10 in the figure is applied to the first detection parameter acquisition unit of a long-distance detection and rapid positioning system for a leakage channel in a dam core wall.

[0078] Furthermore, this embodiment monitors the core wall of the dam in area A; the detection parameters include resistivity and magnetic flux. There are 4 dams in area A, and the core wall information of the 4 dams is shown in Table 1.

[0079] Table 1 Information of the core wall of the dam in area A

[0080]

[0081] Reference Figure 2 S20 in the figure is applied to an abnormality judgment unit of a long-distance detection and rapid positioning system for a leakage channel in a dam core wall.

[0082] Furthermore, the process of determining whether the detection parameters are abnormal includes: performing spectral analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal. Based on the detection parameter inversion, the horizontal mutation starting point, horizontal mutation end point, mutation center point, vertical mutation starting point and vertical mutation end point on the survey line section are identified, and an initial leakage infiltration area is generated. This embodiment detects that the core walls of the dam core wall numbers 2 and 3 are abnormal, and constructs a core wall coordinate system. The X-axis of the core wall coordinate system is parallel to the ground plane, the Y-axis is perpendicular to the ground plane, and the Y-axis is located on the leftmost side of the dam. Figure 4 and Figure 5 , generating the initial seepage and infiltration areas of the two dam core walls, Figure 4 B1, B2, B3, B4 and B5 in the figure represent the vertical mutation starting point, horizontal mutation end point, vertical mutation end point, horizontal mutation starting point and mutation center point on the survey line section with serial number 2 of the dam core wall, respectively; Figure 5 A1, A2, A3, A4 and A5 in the figure represent the vertical mutation starting point, lateral mutation end point, vertical mutation end point, lateral mutation starting point and mutation center point on the survey line section with serial number 3 of the dam core wall.

[0083] In this step of the present embodiment, the resistivity and magnetic flux of the physical detection line are collected, and the mutation point is detected according to the spectrum analysis, and the detection parameter is judged to be abnormal according to whether the mutation point exists. When the detection parameter is abnormal, the initial leakage infiltration area is obtained according to the mutation point. This method can accurately identify the leakage infiltration area of ​​the dam core wall, providing a basis for the subsequent leakage location.

[0084] Reference Figure 2 S30 and S40 are applied to a detection parameter second acquisition unit and a detection leakage infiltration area acquisition unit of a long-distance detection and rapid positioning system for a leakage channel of a dam core wall.

[0085] Furthermore, T is a fixed value, and T can be 1 minute, 5 minutes, 10 minutes, etc.

[0086] Reference Figure 2 S50 is applied to a leakage information acquisition unit of a long-distance detection and rapid positioning system for leakage channels in a dam core wall.

[0087] Furthermore, the step of obtaining the central point position of the leakage area includes:

[0088] A wetting and spreading coordinate system is constructed based on the initial leakage and infiltration area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and infiltration area; the plane formed by the X-axis and the Y-axis of the wetting and spreading coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the wetting and spreading coordinate system are perpendicular to each other;

[0089] Select M infiltration and propagation directions, obtain the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generate M infiltration boundary point coordinate sets; in this embodiment, 10 infiltration and propagation directions are selected, and the angles between two adjacent infiltration and propagation directions are the same; the steps for obtaining the infiltration boundary points are:

[0090] Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area;

[0091] The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

[0092] In this step of the present embodiment, multiple infiltration spreading directions are selected, and the infiltration boundary points of the leakage infiltration area at different times are detected. This method can obtain the spreading status of the initial leakage infiltration area at different times, and provide an accurate data basis for the subsequent leakage rate calculation.

[0093] Furthermore, each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration spreading direction;

[0094] The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; by analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the initial leakage infiltration area at time;

[0095] Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period;

[0096] According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period;

[0097] According to the coordinates of the infiltration boundary points of M infiltration propagation directions at each time point, the leakage infiltration area to be tested at K time points is obtained; when the overlap degree between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap degree, the leakage infiltration area to be tested is eliminated; the area of ​​the remaining leakage infiltration area to be tested is detected, and the center point position of the leakage infiltration area to be tested with the smallest area is taken as the center point position of the leakage area.

[0098] Furthermore, the calculation formula of the overlap degree is:

[0099] ;

[0100] in, It is expressed as the overlap between the kth leakage infiltration area to be tested and the initial leakage infiltration area; It is represented as the overlapping area of ​​the kth leakage infiltration area to be tested and the initial leakage infiltration area; It is expressed as the area of ​​the kth leakage infiltration area to be tested.

[0101] In order to verify the effectiveness of the method for obtaining the center point position of the leakage area in this embodiment, this embodiment also calculates the center point position of the leakage area obtained by the initial leakage infiltration area, and compares the center point position of the leakage area obtained by the two methods with the actual center point position according to Figure 4 and Figure 5 The center point of the leakage area obtained by the two methods is represented in the core wall coordinate system shown in FIG. 1 , and the position deviation between the center point of the leakage area obtained by the two methods and the actual center point is calculated. The calculation results are shown in Table 2.

[0102] Table 2 Example 1 Leakage area center point positioning verification table

[0103]

[0104] Method 1 in Table 2 is the method proposed in this embodiment; Method 2 is to use the initial leakage infiltration area to calculate the center point of the initial leakage infiltration area and obtain the position of the center point of the leakage area. It can be seen from Table 2 that the position deviation between the center point position of the leakage area detected on the core walls of the two dams and the actual center point position provided by the method for obtaining the center point of the leakage area provided in this embodiment is better than that of Method 2, so the method proposed in this embodiment can more accurately locate the center point of the leakage area.

[0105] This step of the present embodiment records the initial abnormal time and the initial leakage infiltration area when the detection parameters are detected to be abnormal, and obtains the detected leakage infiltration area at the set time interval. This embodiment infers the leakage infiltration area to be tested before the initial abnormal time based on the changes of the boundary points of the detected leakage infiltration area at different times, and selects the center point position of the leakage infiltration area to be tested with the smallest area as the center point position of the leakage area. This method can reduce the area of ​​the leakage infiltration area, obtain the center point position of the leakage area more accurately, and thus improve the accuracy of the evaluation of the leakage channel of the dam core wall.

[0106] Furthermore, the step of obtaining the pore size of the leakage area includes:

[0107] Obtaining the leakage parameters of the core wall of the dam; the leakage parameters include the leakage value of the core wall material, the leakage rate and the leakage pressure difference; the leakage value of the core wall material is the permeability coefficient of the core wall; the leakage pressure difference is expressed as the absolute value of the water pressure difference on both sides of the center point of the leakage area;

[0108] The calculation formula of the pore size of the leakage area is:

[0109] ;

[0110] Among them, SLCC represents the pore size of the leakage area; Expressed as leakage rate; It is expressed as the leakage value of the core wall material; Expressed as leakage pressure difference; Expressed as leakage rate coefficient; Expressed as the core wall material influence coefficient; It is expressed as leakage pressure difference influence coefficient; the leakage rate coefficient, the core wall material influence coefficient and the leakage pressure difference influence coefficient are obtained through training with a large amount of historical data; Expressed as an exponential function with a natural constant as base.

[0111] In this embodiment, the step obtains the aperture of the leakage area based on the leakage rate, the core wall material and the water pressure difference on both sides of the core wall to detect the aperture of the leakage channel. This method can accurately calculate the aperture of the leakage channel according to the intensity and extension of the leakage, thereby improving the accuracy of the leakage channel assessment of the dam core wall.

[0112] Furthermore, the leakage rate acquisition process includes:

[0113] The leakage infiltration area is obtained from the initial abnormal time according to the M infiltration boundary point coordinate sets. arrive In the time, every The unit leakage rate of time; the calculation formula of the unit leakage rate is:

[0114] ;

[0115] in, Expressed as The unit leakage rate of the leakage infiltration area at the moment; the value range of i is 1 to N; Expressed as The distance between the infiltration boundary point and the origin corresponding to the j-th infiltration spreading direction in the detected leakage infiltration area at the moment; when i=1, It is expressed as the distance between the infiltration boundary point and the origin corresponding to the j-th infiltration propagation direction in the initial leakage infiltration area; the leakage rate is the average value of N unit leakage rates.

[0116] In this step of the present embodiment, according to the changes in the infiltration boundary points at different times in different infiltration propagation directions of the leakage infiltration area, the unit leakage rate in different time periods is obtained, and the unit leakage rate in different time periods is averaged to obtain the leakage rate of the leakage infiltration area. This method can accurately obtain the leakage rate of the leakage infiltration area and provide a data basis for calculating the aperture of the leakage area.

[0117] Reference Figure 2 The S60 in the embodiment is applied to a leakage grade acquisition unit of a long-distance detection and rapid positioning system for leakage channels in the core wall of a dam.

[0118] Furthermore, the process of obtaining the leakage level includes:

[0119] Obtain the straight-line distance L between the center point of the leakage area and the center point of the bottom of the dam core wall; at the same time, detect the angle between the line connecting the center point of the leakage area and the center point of the bottom of the dam core wall and the ground plane ; The angle is less than 90 degrees;

[0120] Calculate the leakage risk value using the following formula:

[0121] ;

[0122] in, Expressed as leakage risk value; Expressed as core wall width; It is expressed as the core wall height; Expressed as leakage aperture parameter; It is expressed as the leakage lateral parameter; It is expressed as a leakage height parameter; the leakage aperture parameter, the leakage lateral parameter and the leakage height parameter are obtained through training with a large amount of historical data; It is expressed as an exponential function with a natural constant as the base; the leakage risk level is obtained according to the leakage risk value, and an early warning is issued. This embodiment obtains four leakage risk level thresholds, namely x1, x2, x3 and x4; five leakage risk levels are set according to the four leakage risk level thresholds, and the lower the leakage risk level, the lower the threat of leakage; when the leakage risk value is less than x1, the corresponding leakage risk level is level one; when the leakage risk value is within the interval When the leakage risk value is within the range When the leakage risk value is within the range When x4 is greater than or equal to x4, the corresponding leakage risk level is level 4; when the leakage risk value is greater than or equal to x4, the corresponding leakage risk level is level 5. After calculation, it is found that the leakage risk level of the dam core wall with serial number 2 is level 4; the leakage risk level of the dam core wall with serial number 3 is level 3.

[0123] This step of the present embodiment obtains the leakage risk value based on the distance between the center point of the leakage area and the center point of the bottom of the core wall and the leakage aperture. This method combines the location and size of the leakage point. The location information reflects the specific location of the leakage on the core wall, while the size reflects the scale and extension of the leakage. Taking these two factors into consideration, the actual risk of leakage to the core wall of the dam can be more accurately assessed, so as to formulate management and response strategies in a targeted manner.

[0124] Embodiment 2

[0125] Reference Figure 2 S10 in the figure is applied to the first detection parameter acquisition unit of a long-distance detection and rapid positioning system for a leakage channel in a dam core wall.

[0126] Furthermore, this embodiment monitors the core wall of the dam in area B; the detection parameters include resistivity and magnetic flux. There are 6 dams in area B, and the core wall information of the 6 dams is shown in Table 3.

[0127] Table 3 Information of the core wall of the dam in area B

[0128]

[0129] Reference Figure 2 S20 in the figure is applied to an abnormality judgment unit of a long-distance detection and rapid positioning system for a leakage channel in a dam core wall.

[0130] Furthermore, the process of determining whether the detection parameters are abnormal includes: performing spectrum analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal. This embodiment detects that the core walls of the dam with serial numbers 4 and 5 have abnormalities. Based on the detection parameters, the horizontal mutation starting point, horizontal mutation end point, mutation center point, vertical mutation starting point and vertical mutation end point on the survey line section are inverted to identify and generate the initial leakage infiltration area.

[0131] Reference Figure 2 S30 and S40 are applied to a detection parameter second acquisition unit and a detection leakage infiltration area acquisition unit of a long-distance detection and rapid positioning system for a leakage channel of a dam core wall.

[0132] Furthermore, T is a fixed value, and T can be 1 minute, 5 minutes, 10 minutes, etc.

[0133] Reference Figure 2 S50 is applied to a leakage information acquisition unit of a long-distance detection and rapid positioning system for leakage channels in a dam core wall.

[0134] Furthermore, the step of obtaining the central point position of the leakage area includes:

[0135] A wetting and spreading coordinate system is constructed based on the initial leakage and infiltration area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and infiltration area; the plane formed by the X-axis and the Y-axis of the wetting and spreading coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the wetting and spreading coordinate system are perpendicular to each other;

[0136] Select M infiltration and propagation directions, obtain the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generate M infiltration boundary point coordinate sets; in this embodiment, 10 infiltration and propagation directions are selected, and the angles between two adjacent infiltration and propagation directions are the same;

[0137] The steps for obtaining the infiltration boundary point are:

[0138] Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area;

[0139] The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

[0140] Furthermore, each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration spreading direction;

[0141] The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time;

[0142] By analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the initial leakage infiltration area at time;

[0143] Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period;

[0144] According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period;

[0145] According to the coordinates of the infiltration boundary points of M infiltration propagation directions at each time point, the leakage infiltration area to be tested at K time points is obtained; when the infiltration overlap between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap, the leakage infiltration area to be tested is eliminated; the area of ​​the remaining leakage infiltration area to be tested is detected, and the center point position of the leakage infiltration area to be tested with the smallest area is taken as the center point position of the leakage area.

[0146] In order to verify the effectiveness of the method for obtaining the center point position of the leakage area in this embodiment, this embodiment also calculates the center point position of the leakage area obtained from the initial leakage infiltration area, and compares the center point position of the leakage area obtained by the two methods with the actual center point position. The comparison results are shown in Table 4.

[0147] Table 4 Example 2 Leakage area center point positioning verification table

[0148]

[0149] Method 1 in Table 4 is the method proposed in this embodiment; Method 2 is to use the initial leakage infiltration area to calculate the center point of the initial leakage infiltration area and obtain the position of the center point of the leakage area. It can be seen from Table 4 that the position deviation between the center point position of the leakage area detected on the core walls of the two dams and the actual center point position provided by the method for obtaining the center point of the leakage area provided in this embodiment is better than that of Method 2, so the method proposed in this embodiment can more accurately locate the center point of the leakage area.

[0150] Furthermore, the step of obtaining the pore size of the leakage area includes:

[0151] The leakage parameters of the core wall of the dam are obtained; the leakage parameters include the leakage value of the core wall material, the leakage rate and the leakage pressure difference; the leakage value of the core wall material is the permeability coefficient of the core wall; the leakage pressure difference is expressed as the absolute value of the water pressure difference on both sides of the center point of the leakage area; the calculation formula of the aperture of the leakage area is:

[0152] ;

[0153] Among them, SLCC represents the pore size of the leakage area; Expressed as leakage rate; It is expressed as the leakage value of the core wall material; Expressed as leakage pressure difference; Expressed as leakage rate coefficient; Expressed as the core wall material influence coefficient; Expressed as leakage pressure difference influence coefficient; Expressed as an exponential function with a natural constant as base.

[0154] Furthermore, the leakage rate acquisition process includes:

[0155] The leakage infiltration area is obtained from the initial abnormal time according to the M infiltration boundary point coordinate sets. arrive In the time, every The unit leakage rate of time; the calculation formula of the unit leakage rate is:

[0156] ;

[0157] in, Expressed as The unit leakage rate of the leakage infiltration area at the moment; the value range of i is 1 to N; Expressed as The distance between the infiltration boundary point and the origin corresponding to the j-th infiltration propagation direction in the detected leakage infiltration area at the moment; the leakage rate is the average value of N unit leakage rates.

[0158] Reference Figure 2 The S60 in the embodiment is applied to a leakage grade acquisition unit of a long-distance detection and rapid positioning system for leakage channels in the core wall of a dam.

[0159] Furthermore, the process of obtaining the leakage level includes:

[0160] Obtain the straight-line distance L between the center point of the leakage area and the center point of the bottom of the dam core wall; at the same time, detect the angle between the line connecting the center point of the leakage area and the center point of the bottom of the dam core wall and the ground plane ; The angle is less than 90 degrees;

[0161] Calculate the leakage risk value using the following formula:

[0162] ;

[0163] in, Expressed as leakage risk value; Expressed as core wall width; It is expressed as the core wall height; Expressed as leakage aperture parameter; It is expressed as the leakage lateral parameter; Expressed as leakage height parameter; It is expressed as an exponential function with a natural constant as the base; the leakage risk level is obtained according to the leakage risk value, and an early warning is issued.

[0164] This embodiment obtains four leakage risk level thresholds, namely x1, x2, x3 and x4; five leakage risk levels are set according to the four leakage risk level thresholds, and the lower the leakage risk level, the lower the threat of leakage; when the leakage risk value is less than x1, the corresponding leakage risk level is level one; when the leakage risk value is within the interval When the leakage risk value is within the range When the leakage risk value is within the range When x4 is greater than or equal to x4, the corresponding leakage risk level is level 4; when the leakage risk value is greater than or equal to x4, the corresponding leakage risk level is level 5. After calculation, it is found that the leakage risk level of the dam core wall with serial number 4 is level 3; the leakage risk level of the dam core wall with serial number 5 is level 2.

[0165] When the detection parameters are detected to be abnormal, this embodiment records the initial abnormal time and the initial leakage infiltration area, and obtains the detected leakage infiltration area at the set time interval. This embodiment reversely deduces the leakage infiltration area to be detected before the initial abnormal time based on the changes in the boundary points of the detected leakage infiltration area at different times, and selects the center point position of the leakage infiltration area to be detected with the smallest area as the center point position of the leakage area. This method can reduce the area of ​​the leakage infiltration area, more accurately obtain the center point position of the leakage area, and thus improve the accuracy of the leakage channel assessment of the dam core wall. This embodiment obtains the leakage area aperture based on the leakage rate, the core wall material and the water pressure difference on both sides of the core wall to detect the leakage channel aperture. This method can accurately calculate the aperture of the leakage channel based on the intensity and extension of the leakage, thereby improving the accuracy of the leakage channel assessment of the dam core wall. This embodiment obtains the leakage risk value based on the distance between the center point of the leakage area and the center point of the bottom of the core wall and the leakage aperture. This method combines the location and size of the leakage point. The location information reflects the specific location of the leakage on the core wall, while the size reflects the scale and extension of the leakage. Taking these two factors into consideration, the actual risk of leakage to the core wall of the dam can be more accurately assessed, so that targeted management and response strategies can be formulated.

[0166] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. A method for rapid positioning of leakage channels in dam core walls over long distances, characterized in that: include: Establish a geophysical detection line on the core wall of the dam and regularly collect detection parameters of the geophysical detection line; The detection parameters include resistivity and magnetic flux; when the detection parameters are abnormal, the initial abnormal time is recorded and initial leakage infiltration area; from Start, every The detection parameters are collected once at a time, and the detection parameters are collected N times in total; Obtaining the corresponding detection leakage 1 infiltration area according to the detection parameters of the object detection line collected each time; Obtaining the central point position and the aperture of the leakage area of ​​the dam core wall according to the N detected leakage infiltration areas and the initial leakage infiltration area; The step of obtaining the central point position of the leakage area includes: A wetting and spreading coordinate system is constructed based on the initial leakage and infiltration area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and infiltration area; the plane formed by the X-axis and the Y-axis of the wetting and spreading coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the wetting and spreading coordinate system are perpendicular to each other; Selecting M infiltration and propagation directions, obtaining the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generating M infiltration boundary point coordinate sets; Each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration propagation direction; The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; by analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the initial leakage infiltration area at time; Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period; According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period; According to the coordinates of the infiltration boundary points of M infiltration propagation directions at each time point, the leakage infiltration area to be tested at K time points is obtained; when the overlap between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap, the leakage infiltration area to be tested is eliminated; the area of ​​the remaining leakage infiltration area to be tested is detected, and the center point position of the leakage infiltration area to be tested with the smallest area is taken as the center point position of the leakage area; The step of obtaining the pore size of the leakage area comprises: Obtaining the leakage parameters of the core wall of the dam; the leakage parameters include the leakage value of the core wall material, the leakage rate and the leakage pressure difference; the leakage value of the core wall material is the permeability coefficient of the core wall; the leakage pressure difference is expressed as the absolute value of the water pressure difference on both sides of the center point of the leakage area; The calculation formula of the pore size of the leakage area is: ； Among them, SLCC represents the pore size of the leakage area; Expressed as leakage rate; It is expressed as the leakage value of the core wall material; Expressed as leakage pressure difference; Expressed as leakage rate coefficient; Expressed as the core wall material influence coefficient; Expressed as leakage pressure difference influence coefficient; It is expressed as an exponential function with a natural constant as base; The leakage level is obtained according to the central point position of the leakage area and the aperture of the leakage area, and an early warning is issued.

2. The method for rapid positioning of leakage channels in dam core walls over long distances according to claim 1, characterized in that: The process of determining whether the detection parameters are abnormal includes: performing spectrum analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal.

3. The method for rapid positioning of leakage channels in dam core walls over long distances according to claim 1, characterized in that: The steps for obtaining the infiltration boundary point are: Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area; The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

4. The method for rapid positioning of leakage channels in dam core walls over long distances according to claim 1, characterized in that: The leakage rate acquisition process includes: The leakage infiltration area is obtained from the initial abnormal time according to the M infiltration boundary point coordinate sets. arrive In the time, every The unit leakage rate of time; the calculation formula of the unit leakage rate is: ； in, Expressed as The unit leakage rate of the leakage infiltration area at the moment; the value range of i is 1 to N; Expressed as The distance between the infiltration boundary point and the origin corresponding to the j-th infiltration propagation direction in the detected leakage infiltration area at the moment; the leakage rate is the average value of N unit leakage rates.

5. The method for rapid positioning of leakage channels in dam core walls over long distances according to claim 1, characterized in that: The process of obtaining the leakage level includes: Obtain the straight-line distance L between the center point of the leakage area and the center point of the bottom of the dam core wall; at the same time, detect the angle between the line connecting the center point of the leakage area and the center point of the bottom of the dam core wall and the ground plane ; The angle is less than 90 degrees; Calculate the leakage risk value using the following formula: ； in, Expressed as leakage risk value; Expressed as core wall width; It is expressed as the core wall height; Expressed as leakage aperture parameter; It is expressed as the leakage lateral parameter; Expressed as leakage height parameter; It is expressed as an exponential function with a natural constant as the base; the leakage risk level is obtained according to the leakage risk value, and an early warning is issued.

6. A long-distance detection and rapid positioning system for leakage channels in dam core walls, characterized in that: include: The first detection parameter acquisition unit is used to regularly collect the detection parameters of the object detection line on the core wall of the dam; Abnormality judgment unit: used to judge whether the detection parameters are abnormal, and record the initial abnormal time when the detection parameters are abnormal and initial leakage infiltration area; The second detection parameter acquisition unit is used to Start, every The detection parameters are collected once at a time, and the detection parameters are collected N times in total; A detection leakage infiltration area acquisition unit is used to acquire a corresponding detection leakage infiltration area according to the detection parameters of the object detection line collected each time; A leakage information acquisition unit: used for acquiring the central point position of the leakage area and the aperture of the leakage area of ​​the dam core wall according to the N detected leakage infiltration areas and the initial leakage infiltration area; The step of obtaining the central point position of the leakage area includes: A wetting and spreading coordinate system is constructed based on the initial leakage and infiltration area; the origin of the wetting and spreading coordinate system is the center point of the initial leakage and infiltration area; the plane formed by the X-axis and the Y-axis of the wetting and spreading coordinate system is perpendicular to the ground plane; the X-axis and the Y-axis of the wetting and spreading coordinate system are perpendicular to each other; Selecting M infiltration and propagation directions, obtaining the infiltration boundary point coordinates of the initial leakage infiltration area and the N detected leakage infiltration areas in the M infiltration and propagation directions, and generating M infiltration boundary point coordinate sets; Each of the infiltration boundary point coordinate sets includes the infiltration boundary point coordinates of the initial leakage infiltration area and N detected leakage infiltration areas in the same infiltration propagation direction; The first data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the second data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; by analogy, the Nth data in the infiltration boundary point coordinate set is The coordinates of the infiltration boundary points of the leakage infiltration area detected at the time; the N+1th data in the infiltration boundary point coordinate set is The coordinates of the initial infiltration boundary point of the leakage infiltration area at time; Input each of the infiltration boundary point coordinate sets into the time series model to obtain The predicted values ​​of the infiltration boundary point coordinates at K time points in the infiltration spreading direction corresponding to each infiltration boundary point coordinate set within the time period; According to the predicted values ​​of the infiltration boundary point coordinates at K time points corresponding to the infiltration spreading direction of each infiltration boundary point coordinate set, obtain The coordinates of the infiltration boundary points in M ​​infiltration spreading directions at K time points within the time period; Obtain the leakage infiltration area to be tested at K time points according to the coordinates of the infiltration boundary points in the M infiltration propagation directions at each time point; when the overlap between the leakage infiltration area to be tested and the initial leakage infiltration area is less than the set standard overlap, remove the leakage infiltration area to be tested; detect the area of ​​the remaining leakage infiltration area to be tested, and take the center point position of the leakage infiltration area to be tested with the smallest area as the center point position of the leakage area; The step of obtaining the pore size of the leakage area comprises: Obtaining the leakage parameters of the core wall of the dam; the leakage parameters include the leakage value of the core wall material, the leakage rate and the leakage pressure difference; the leakage value of the core wall material is the permeability coefficient of the core wall; the leakage pressure difference is expressed as the absolute value of the water pressure difference on both sides of the center point of the leakage area; The calculation formula of the pore size of the leakage area is: ； Among them, SLCC represents the pore size of the leakage area; Expressed as leakage rate; It is expressed as the leakage value of the core wall material; Expressed as leakage pressure difference; Expressed as leakage rate coefficient; Expressed as the core wall material influence coefficient; Expressed as leakage pressure difference influence coefficient; It is expressed as an exponential function with a natural constant as base; Leakage level acquisition unit: used to acquire the leakage level according to the center point position of the leakage area and the aperture of the leakage area, and to issue an early warning.

7. The long-distance detection and rapid positioning system for leakage channels in dam core walls according to claim 6 is characterized in that: The process of determining whether the detection parameters are abnormal includes: performing spectrum analysis on the detection parameters collected in the current period to obtain mutation points on the detection parameters; when the mutation point does not exist, it indicates that the detection parameters are not abnormal; when the mutation point exists, it indicates that the detection parameters are abnormal.

8. The long-distance detection and rapid positioning system for leakage channels in dam core walls according to claim 6 is characterized in that: The steps for obtaining the infiltration boundary point are: Make M rays from the origin of the infiltration and propagation coordinate system in the M infiltration and propagation directions; take the intersection of the M rays and the boundary point of the initial leakage and infiltration area as the infiltration boundary point corresponding to the initial leakage and infiltration area; The intersection point of the M rays and each boundary point of the detected leakage infiltration area is taken as the infiltration boundary point corresponding to the detected leakage infiltration area.

Images