CN-121783453-B - Anti-seepage structure integrity assessment method and system based on resistivity imaging
Abstract
The invention belongs to the technical field of impervious structure evaluation, and in particular relates to a impervious structure integrity evaluation method and system based on resistivity imaging, wherein the method comprises the steps of firstly collecting complex admittance data of monitoring points and extracting phase loss angle characteristics, and stripping interference of humidity fluctuation of an underground environment through phase information; the method comprises the steps of constructing an imaginary part admittance gradient characteristic by combining a structural anisotropy correction factor, enhancing the identification degree of a tiny breakage edge, calculating leakage probability weight by utilizing reference gradient deviation, identifying leakage clusters by using a self-adaptive seed point screening and eight-neighborhood connected domain marking algorithm, obtaining the area of a continuous leakage area, and finally calculating integrity score by introducing a nonlinear model of risk sensitivity index and area penalty weight, so as to realize accurate early warning and positioning. The invention effectively solves the problem of poor leakage identification accuracy under complex background interference, and remarkably improves the robustness and scientificity of the evaluation result.
Inventors
- ZHANG PEIKUN
- ZHANG SENLIN
- WANG RUJUN
- ZHAO ZHONGTAO
- ZHANG LIZHI
- DIAO MINGMING
- SUN PENG
- DONG JUAN
- WANG LEI
- FAN YINGGUANG
- WANG LANSHENG
- WANG GUODONG
- Qiao Erwei
- ZHAO HAOWEI
- CHANG YUAN
- ZHENG YAYA
Assignees
- 山东华新通信科技有限公司
Dates
- Publication Date
- 20260508
- Application Date
- 20260303
Claims (9)
- 1. A resistivity imaging-based barrier structure integrity assessment method, comprising: Arranging a monitoring electrode array in a covering soil layer above the seepage-proofing structure, continuously collecting a response voltage sequence and a response current sequence of each monitoring point, obtaining complex admittance data of each monitoring point through fast Fourier transformation, and extracting phase loss angle characteristics of each monitoring point according to the complex admittance data; Capturing local response jump based on the second partial derivative of the phase loss angle characteristic relative to the space coordinate, calculating the square of the difference between the phase loss angle characteristic and the global average phase loss angle of the characteristic field to represent the deviation degree, carrying out weighted coupling on the local response jump and the deviation degree through a structural anisotropy correction factor, and constructing the imaginary part admittance gradient characteristic of each monitoring point; acquiring leakage probability weights of all monitoring points according to the difference of the imaginary admittance gradient characteristics compared with the reference gradient characteristic values; acquiring the leakage probability weights of all the monitoring points, arranging the leakage probability weights according to the numerical values from large to small, selecting candidate leakage seed points based on abnormal screening quantile coefficients, identifying leakage clusters based on the candidate leakage seed points, and acquiring the areas of continuous leakage areas corresponding to the monitoring points by combining the areas of grid cells; Based on the leakage probability weight, the area of the continuous leakage area and the total monitoring point position number, obtaining an integrity assessment score of the seepage prevention structure, wherein the calculation formula of the integrity assessment score of the seepage prevention structure is as follows: In the formula (I), in the formula (II), Scoring the integrity assessment; The total monitoring point position number is; Is the first Leakage probability weights of the monitoring points; Is a risk sensitive index, and has a value ranging from 1.5 to 3.0; the area penalty weight is the value range of 5.0 to 15.0; Is the first The areas of the continuous seepage areas to which the monitoring points belong; Is the total area of the monitored area; and carrying out leakage early warning and position positioning on the seepage prevention structure based on the integrity evaluation score.
- 2. The resistivity imaging-based anti-seepage structure integrity assessment method according to claim 1, wherein the complex admittance data of each monitoring point position comprises a real part admittance component and an imaginary part admittance component, wherein the current real part component and the voltage real part component of each monitoring point position form a real part admittance component, the current imaginary part component and the voltage imaginary part component of each monitoring point position form an imaginary part admittance component, and the current real part component, the current imaginary part component, the voltage real part component and the voltage imaginary part component of each monitoring point position are obtained by performing fast fourier transformation on the response voltage sequence and the response current sequence.
- 3. The method for evaluating the integrity of an impermeable structure based on resistivity imaging according to claim 2, wherein the phase loss angle characteristics of each monitoring point are calculated as follows: ; in the formula, Monitoring the phase loss angle characteristic of the point location; Is an arctangent function; Is the real component of the current; is the voltage imaginary component; Is the imaginary component of the current; is the real component of the voltage.
- 4. The resistivity imaging-based barrier structure integrity assessment method of claim 1, wherein the imaginary admittance gradient feature of each monitoring point is calculated as: ; in the formula, The imaginary part admittance gradient characteristics of the point positions are monitored; is a three-dimensional space Laplacian; Monitoring the phase loss angle characteristic of the point location; is a structural anisotropy correction factor, and has a value ranging from 1.0 to 1.5; is the global average phase loss angle of the feature field.
- 5. The resistivity imaging-based method for evaluating the integrity of an impermeable structure of a vehicle according to claim 1, wherein the leakage probability weight of each monitoring point is calculated as: ; in the formula, The leakage probability weight of the point location is monitored; Is a natural exponential function; The gain sensitivity coefficient is in the range of 10.0 to 20.0; the imaginary part admittance gradient characteristics of the point positions are monitored; Is the reference gradient eigenvalue.
- 6. The method for evaluating the integrity of an impermeable structure based on resistivity imaging of claim 1, wherein said obtaining the area of the continuous leak area corresponding to each monitoring point comprises: Selecting the leakage probability weight Monitoring points in the range of the first 5% proportion and marking the selected monitoring points as candidate seepage seed points; Performing space growth by adopting an eight-neighborhood connected domain marking algorithm and taking the candidate leakage seed points as cores, and classifying the candidate leakage seed points with adjacent space positions into the same leakage cluster; and counting the number of candidate leakage seed points contained in each leakage cluster and combining the areas of the grid cells to obtain the areas of continuous leakage areas corresponding to the monitoring points.
- 7. The resistivity imaging-based barrier structure integrity assessment method of claim 1, wherein the response voltage sequence and the response current sequence of each monitoring point are acquired while alternating current signals with a frequency range of 10Hz to 1kHz are injected into the monitoring electrode array by using a multi-channel complex impedance analyzer.
- 8. The resistivity imaging-based barrier structure integrity assessment method of claim 1, wherein performing leak warning and location positioning of the barrier structure based on the integrity assessment score comprises: When the integrity evaluation score is continuously lower than a preset alarm threshold three times, judging that the seepage prevention structure has real seepage risk and triggering early warning; and generating a leakage position thermodynamic diagram of the area to be monitored by using the spatial coordinates of each monitoring point and the corresponding imaginary admittance gradient characteristics and adopting a bilinear interpolation algorithm.
- 9. A resistivity imaging based barrier structure integrity evaluation system comprising a processor and a memory, the memory storing computer program instructions that when executed by the processor implement a resistivity imaging based barrier structure integrity evaluation method according to any one of claims 1-8.
Description
Anti-seepage structure integrity assessment method and system based on resistivity imaging Technical Field The invention relates to the technical field of seepage-proofing structure evaluation. More particularly, the present invention relates to a resistivity imaging-based barrier structure integrity assessment method and system. Background In the field of integrity monitoring of hydraulic seepage-proofing engineering, environmental sanitation landfill and chemical storage facilities, a seepage-proofing structure is a core barrier for blocking high-concentration pollutant from overflowing and protecting underground water resources. In order to ensure the reliability of the anti-seepage system in long-term operation, a monitoring electrode array is generally arranged in a covering soil layer above an anti-seepage structure, a main flow apparent resistivity detection technology is used for positioning leakage points by injecting direct current into the soil layer and collecting potential distribution and utilizing local resistivity distortion caused by a conductive path formed at the damaged part of the anti-seepage film, and the monitoring means based on the resistive characteristics has certain effectiveness in controlled environments such as a laboratory and the like, and is a basic scheme for online evaluation of the anti-seepage structure in the engineering community at present. However, under a complex actual engineering environment, the soil resistivity has extremely high sensitivity to the water content, and rainfall, seasonal fluctuation of the ground water level or local water storage of a covering soil layer caused by uneven compaction degree can cause severe uneven deviation of the background resistivity, and the false electrical abnormality caused by fluctuation of the environment humidity is often highly aliased with current distortion caused by real leakage points on the aspects of signal amplitude and spatial distribution, so that the existing monitoring system is difficult to accurately distinguish local conductivity rise caused by the increase of the soil humidity from charge leakage caused by physical breakage of an impermeable film. In addition, the existing evaluation logic is mostly based on simple threshold judgment or linear weighting, complex admittance information of underground media under an alternating electric field is not fully considered, so that the underground media lack sufficient identification dimension in a signal extraction stage, and when tiny breakage or multi-point aggregation leakage is faced, the judgment mode of a fixed threshold is extremely easy to be interfered by environmental noise to cause report missing, or humidity fluctuation with a large area is misjudged to be serious leakage. Due to the signal characteristic ambiguity caused by the environmental fluctuation, the accuracy of leakage identification is reduced, the evaluation system cannot give out a quantitative score matched with the physical damage degree of the seepage-proofing structure, and the reliability and pertinence of the seepage-proofing engineering safety early warning are severely restricted. Disclosure of Invention In order to solve the technical problems that in the prior art, under the conditions of humidity fluctuation and complicated underground background interference, the environmental influence and the real leakage signal are difficult to distinguish, so that the leakage identification accuracy is poor and the structural integrity cannot be quantitatively evaluated, the invention provides a scheme in the following aspects. The invention provides an anti-seepage structure integrity assessment method based on resistivity imaging, which comprises the steps of arranging a monitoring electrode array in a covering soil layer above an anti-seepage structure, continuously collecting response voltage sequences and response current sequences of all monitoring points, obtaining complex admittance data of all monitoring points through fast Fourier transformation, extracting phase loss angle characteristics of all monitoring points according to the complex admittance data, capturing local response jump based on second-order partial derivatives of the phase loss angle characteristics relative to space coordinates, calculating the square of difference between the phase loss angle characteristics and global average phase loss angles of a characteristic field to represent deviation degree, carrying out weighted coupling on the local response jump and the deviation degree through a structural anisotropy correction factor, constructing imaginary part admittance gradient characteristics of all monitoring points, obtaining candidate leakage seed points based on the number of the grid probability weights of all monitoring points and arranged from large to small according to the number, obtaining candidate leakage seed point clusters based on the candidate leakage coefficient, obtaining the total leakage a