Assembly, method and device for improving the precision of dam leakage double-mode parallel electrical method inversion

Abstract

The application discloses a kind of components, methods and devices for improving dam leakage double-mode parallel electric method inversion precision, belong to water conservancy engineering safety detection technical field, the component includes: cable, it is arranged in dam top longitudinal axis direction in parallel;Sensing element includes power supply electrode and measuring electrode, power supply electrode and measuring electrode are isolated from each other, and power supply electrode and measuring electrode are arranged on the two sides of cable, and are arranged at equal intervals, the present scheme can solve the current parallel electric method massive data body utilization rate is lower, power supply electrode is switched into measuring electrode There is polarization potential that has not all dissipated phenomenon, different data bodies are not fully embodied in the characteristics of electric method device in joint inversion process and different geotechnical area inversion resistivity depth exists difference and other problems.

Description

Technical Field

[0001] This invention belongs to the field of water conservancy engineering safety detection technology, specifically relating to a component, method and device for improving the accuracy of dual-mode parallel electrical resistivity tomography inversion of dam leakage. Background Technology

[0002] Parallel electrical resistivity tomography (EPM) employs a quasi-seismic data acquisition mode, achieving synchronous and transient recording of the entire electric field data. This significantly improves work efficiency and avoids the inconsistency in the time domain of steady electric fields formed at different times. It has been widely applied in engineering geology, hydrogeology, environmental surveys, disaster assessment, and mineral exploration. Dam leakage has always been a hot and challenging issue in the water conservancy industry. To identify the causes and lesions of leakage, parallel EPM technology, as a rapid detection method, has been introduced into dam detection, playing a crucial role in ensuring the safe operation of dams.

[0003] However, due to limitations in forward and inverse algorithms, the massive amounts of data collected by parallel resistivity methods are not fully utilized. They are typically converted into data volumes from conventional, standardized devices for processing, resulting in a significant waste of data and hindering the improvement of accuracy in detecting dam leakage hazards. Furthermore, current resistivity inversion relies excessively on initial structured model information, assigning the same weight to data involved in the joint inversion of different data volumes. This fails to highlight the characteristics of a single resistivity device. Moreover, the depth coefficients involved in the constructed forward grid model are essentially uniform, lacking prior information that considers the characteristics of different soil and rock masses in different dam sections and sets different depths. Especially in areas where the dam body and rock mass differ significantly, directly using resistivity inversion results can introduce errors in determining the depth of leakage hazards, negatively impacting dam reinforcement design and cost estimation.

[0004] Invention patent CN 110702587 B combines data from arrangements of the Wenner quadrupole device, Wenner dipole device, and Wenner differential device for joint inversion, which improves the efficiency and accuracy of seepage diagnosis for earth-rock dams to some extent, but still does not consider the data characteristics of different devices. Furthermore, engineering practice has shown that the three-electrode electrical resistivity method is more sensitive to dam seepage anomalies. However, in the parallel electrical resistivity method, the same electrode is used for both the power supply and measurement electrodes. Due to the short measurement time of the parallel electrical resistivity method, when switching from the power supply electrode to the measurement electrode, the electrode polarization effect generated by the metal electrode itself is not completely eliminated, thus introducing system measurement errors into the MNB data acquired by the AM method. Summary of the Invention

[0005] The purpose of this invention is to provide a component, method, and apparatus for improving the accuracy of dual-mode parallel electrical resistivity inversion of dam seepage. This invention can solve problems such as the low utilization rate of massive parallel electrical resistivity data, the phenomenon that polarization potential is not completely dissipated when switching from power supply electrode to measurement electrode, the failure of different data volumes to fully reflect the characteristics of the electrical resistivity device during joint inversion, and the difference in resistivity depth during inversion in different soil and rock regions.

[0006] To solve the above-mentioned technical problems, the present invention is implemented as follows:

[0007] In a first aspect, embodiments of the present invention provide a component for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, comprising:

[0008] The cable is arranged parallel to the longitudinal axis of the dam crest.

[0009] A sensing element, comprising a power supply electrode and a measuring electrode, wherein the power supply electrode and the measuring electrode are isolated from each other and are disposed on both sides of the cable and arranged at equal intervals.

[0010] Optionally, the power supply electrode is a metal electrode used to supply power to the underground to form a stable electric field within the geological body; the measuring electrode is a non-polarized electrode used to measure the potential difference formed by the stable electric field within the geological body.

[0011] Optionally, the cable includes at least two conductors, the first ends of which are respectively connected to the power supply electrode and the measuring electrode, and the second ends of which are both connected to the dual-mode parallel electrical resistivity meter. Parameter data can be transmitted to the dual-mode parallel electrical resistivity meter through the power supply electrode and the measuring electrode.

[0012] Secondly, embodiments of the present invention provide a method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EPM) inversion of dam leakage. The power supply electrodes include a first power supply electrode and a second power supply electrode, and the measuring electrodes include a first measuring electrode and a second measuring electrode. The distance between the first power supply electrode and the first measuring electrode, the distance between the first measuring electrode and the second measuring electrode, and the distance between the second measuring electrode and the second power supply electrode are all equal. The Winner AMN represents the data between the first power supply electrode, the first measuring electrode, and the second measuring electrode, and the Winner MNB represents the data between the first measuring electrode, the second measuring electrode, and the second power supply electrode. The method for improving the accuracy of dual-mode parallel EPM inversion of dam leakage includes:

[0013] S101: Obtain the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at the first, second and third depths, wherein the second depth is greater than the first depth and the third depth is greater than the second depth;

[0014] S102: Perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth.

[0015] S103: Based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, obtain the initial model of the dam structure, and obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB.

[0016] S104: Based on the combined data inversion resistivity map of the Winner AMN and the Winner MNB, obtain the inversion resistivity of the Winner AMN and the Winner MNB at different locations at the same depth, and perform depth correction to obtain the spatial location of the dam leakage.

[0017] Optionally, the weighted correction uses the formula:

[0018]

[0019]

[0020] Where, α n *β n =1; O is the lateral center point of the dam; n is the different depths at the same location on the dam; The resistivity of the Wenner AMN at different depths at the same location on the dam; The resistivity of the Wenner MNB at different depths at the same location on the dam.

[0021] Optionally, after performing weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB at the first depth, the second depth, and the third depth, the method further includes:

[0022] The resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first, second, and third depths after weighted correction are corrected to obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first, second, and third depths.

[0023] Optionally, the correction process uses the following formula:

[0024] When the resistivity of the Winner AMN is located to the left of the center point, the weighted corrected resistivity of the Winner AMN is,

[0025]

[0026] When the resistivity of the Winner AMN is located to the right of the center point, the weighted corrected resistivity of the Winner AMN is,

[0027]

[0028] When the resistivity of the Winner MNB is located to the left of the center point, the weighted corrected resistivity of the Winner MNB is,

[0029]

[0030] When the resistivity of the Winner MNB is located to the right of the center point, the weighted corrected resistivity of the Winner MNB is,

[0031]

[0032] The corrected Winner AMN resistivity and the corrected Winner MNB resistivity are,

[0033]

[0034]

[0035] Where n represents different depths at the same location of the dam, L represents the location to the left of the center point, and R represents the location to the right of the center point.

[0036] Optionally, S104 specifically includes:

[0037] Obtain the ratio of dam body depth to dam body resistivity at different locations within the same dam layer, and the ratio of dam shoulder depth to dam shoulder resistivity.

[0038] Based on the ratio of the dam body depth to the dam body resistivity and the ratio of the dam shoulder depth to the dam shoulder resistivity, the dam body depth correction coefficient and the dam shoulder depth correction coefficient are obtained.

[0039] The spatial location of the dam seepage can be obtained based on the dam body depth correction coefficient, the dam shoulder depth correction coefficient, and different depths.

[0040] Optionally, the ratio of dam body depth to dam body resistivity and the ratio of dam abutment depth to dam abutment resistivity are expressed by the following formulas:

[0041]

[0042]

[0043] Wherein, η is the dam body depth correction coefficient, and ε is the dam abutment depth correction coefficient.

[0044] Thirdly, embodiments of the present invention provide an apparatus for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, comprising:

[0045] The first acquisition module is used to acquire the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at a first depth, a second depth, and a third depth, wherein the second depth is greater than the first depth, and the third depth is greater than the second depth.

[0046] The second acquisition module is used to perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth.

[0047] The third acquisition module is used to obtain an initial model of the dam structure based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, and to obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB.

[0048] The correction module is used to obtain the inversion resistivity of the Wenner AMN and the Wenner MNB at different locations at the same depth based on the joint data inversion resistivity map of the Wenner AMN and the joint data inversion resistivity map of the Wenner MNB, and to perform depth correction, thereby obtaining the spatial location of the dam leakage.

[0049] Fourthly, embodiments of the present invention provide a device, comprising:

[0050] processor;

[0051] Memory used to store processor-executable instructions;

[0052] In this method, the processor is configured to invoke instructions stored in memory to execute the second aspect of the method.

[0053] Fifthly, embodiments of the present invention provide a computer-readable storage medium having computer program instructions stored thereon, which, when executed by a processor, implement the method of the second aspect.

[0054] This scheme extracts the joint resistivity values ​​of the Wenner AMN and Wenner MNB at equal intervals from the AM method data volume acquired by the dual-mode parallel electrical resistivity meter. The apparent resistivity at the center point of a single probe profile and at different depths is used as a normalization parameter for the asymmetry of the apparent resistivity data. The Wenner AMN and Wenner MNB data volumes are normalized separately. All data are then used for forward and inverse calculations in a priori model of the existing dam structure. The inversion results are then processed at different depths, thereby improving the inversion accuracy while reducing the error of the probe depth. Attached Figure Description

[0055] Figure 1 This is a flowchart illustrating a method for improving the accuracy of dual-mode parallel electrical resistivity inversion of dam seepage, provided by an embodiment of the present invention.

[0056] Figure 2 This is a schematic diagram of the structure of a device for improving the accuracy of dual-mode parallel electrical resistivity tomography inversion of dam seepage, provided in an embodiment of the present invention.

[0057] The realization of the objective, functional characteristics and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0058] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of the embodiments of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. Based on the embodiments of this invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this invention.

[0059] The terms "first," "second," etc., used in the specification and claims of this invention are used to distinguish similar objects and are not used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of the invention can be implemented in orders other than those illustrated or described herein, and that the objects distinguished by "first," "second," etc., are generally of the same class and do not limit the number of objects; for example, a first object can be one or more. It should be understood that in the various embodiments of this disclosure, the sequence number of each process does not imply the order of execution; the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this disclosure.

[0060] It should be understood that in this disclosure, "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion, for example, a process, method, system, product, or device that includes a series of steps or units is not necessarily limited to those steps or units that are explicitly listed, but may include other steps or units that are not explicitly listed or that are inherent to such process, method, product, or device.

[0061] It should be understood that in this disclosure, "multiple" refers to two or more. "And / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, "and / or B" can represent: A existing alone, A and B existing simultaneously, or B existing alone. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "Contains A, B, and C", "Contains A, B, and C" means that all three A, B, and C are contained; "Contains A, B, or C" means that one of A, B, and C is contained; "Contains A, B, and / or C" means that any one, two, or three of A, B, and C are contained.

[0062] It should be understood that in this disclosure, "B corresponding to A", "B corresponding to A", "A corresponds to B", or "B corresponds to A" means that B is associated with A, and B can be determined based on A. Determining B based on A does not mean determining B solely based on A; B can also be determined based on A and / or other information. Matching A and B is defined as a similarity between A and B that is greater than or equal to a preset threshold.

[0063] Depending on the context, "if" as used here can be interpreted as "when," "when," "in response to determination," or "in response to detection."

[0064] The following description, in conjunction with the accompanying drawings, details a component, method, and apparatus for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, provided by the present invention, through specific embodiments and application scenarios.

[0065] Example 1

[0066] This invention provides a component for improving the accuracy of dual-mode parallel electrical resistivity tomography (EPR) inversion of dam leakage, comprising:

[0067] The cable is arranged parallel to the longitudinal axis of the dam crest.

[0068] A sensing element, comprising a power supply electrode and a measuring electrode, wherein the power supply electrode and the measuring electrode are isolated from each other and are disposed on both sides of the cable and arranged at equal intervals.

[0069] Optionally, the power supply electrode is a metal electrode used to supply power to the underground to form a stable electric field within the geological body; the measuring electrode is a non-polarized electrode used to measure the potential difference formed by the stable electric field within the geological body.

[0070] Optionally, the cable includes at least two conductors, the first ends of which are connected to the power supply electrode and the measuring electrode, respectively, and the second ends of which are both connected to a dual-mode parallel electrical resistivity meter. Parameter data can be transmitted to the dual-mode parallel electrical resistivity meter through the power supply electrode and the measuring electrode. The parameters include grounding resistance, power supply time, acquisition interval, input waveform, starting electrode, number of electrodes, and acquisition mode. Current and voltage data are acquired in AM method acquisition mode.

[0071] Example 2

[0072] Reference Figure 1 The diagram shows a flowchart of a method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, provided by an embodiment of the present invention.

[0073] This invention provides a method for improving the accuracy of dual-mode parallel electrical resistivity tomography (DEM) inversion of dam leakage, applied to the component for improving the accuracy of dual-mode parallel electrical resistivity tomography of dam leakage as described in the first aspect. The power supply electrodes include a first power supply electrode and a second power supply electrode, and the measuring electrodes include a first measuring electrode and a second measuring electrode. The distance between the first power supply electrode and the first measuring electrode is equal to the distance between the first measuring electrode and the second measuring electrode, and the distance between the second measuring electrode and the second power supply electrode. The Winner AMN represents the data between the first power supply electrode, the first measuring electrode, and the second measuring electrode, and the Winner MNB represents the data between the first measuring electrode, the second measuring electrode, and the second power supply electrode. The method for improving the accuracy of dual-mode parallel electrical resistivity tomography of dam leakage includes:

[0074] S101: Obtain the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at a first depth, a second depth, and a third depth, wherein the second depth is greater than the first depth, and the third depth is greater than the second depth.

[0075] Based on the current and voltage data obtained in the AM method acquisition mode, the resistivity values ​​of Wenner AMN and Wenner MNB can be obtained. Subsequently, the format can be converted according to the horizontal position of the recording point, the isolation coefficient, 1, the resistivity value, etc.

[0076] In some embodiments, the power supply electrode and the measuring electrode are located on opposite sides of the cable, and the distance between them and the cable is less than 0.2m.

[0077] S102: Perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth.

[0078] Optionally, the transverse center point of the detected apparent resistivity profile is designated as point O, and the resistivity of the Wenner AMN and Wenner MNB at different layers at point O is extracted. For example, the apparent resistivity of the Wenner AMN at point O in the first layer is denoted as... The apparent resistivity of the Wenner MNB at point O in the first layer is denoted as... The apparent resistivity of the Wenner AMN at point O in the second layer is denoted as... The apparent resistivity of the Wenner MNB at point O in the second layer is denoted as... The apparent resistivity of the Winner AMN at point O in the third layer is denoted as... The apparent resistivity of the Wenner MNB at point O in the third layer is denoted as... The resistivity representations for different layers, from the shallow layer to the deep layer at point O, are derived sequentially.

[0079] Optionally, at the transverse center point O of the apparent resistivity profile, the apparent resistivity of different layers is less affected by the asymmetry of the triode electrical resistivity device. Therefore, the apparent resistivity values ​​at other transverse locations can be corrected based on the difference in apparent resistivity between the Wenner AMN and Wenner MNB devices at the transverse center point O, thereby improving the accuracy of the inversion.

[0080] The apparent resistivity of Wenner AMN and Wenner MNB at the transverse center point O is extracted, and the apparent resistivity weighting correction coefficient is calculated based on the resistivity difference at the corresponding layer depth.

[0081] The weighted correction uses the following formula:

[0082]

[0083]

[0084] Where, α n *β n =1; O is the lateral center point of the dam; n is the different depths at the same location on the dam; The resistivity of the Wenner AMN at different depths at the same location on the dam; The resistivity of the Wenner MNB at different depths at the same location on the dam.

[0085] Optionally, after weighting and correcting the resistivity of the Winner AMN and the Winner MNB at the first depth, the second depth, and the third depth, the method further includes:

[0086] The resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first, second, and third depths after weighted correction are corrected to obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first, second, and third depths.

[0087] Optionally, the correction process uses the following formula:

[0088] When the resistivity of the Winner AMN is located to the left of the center point, the weighted corrected resistivity of the Winner AMN is,

[0089]

[0090] When the resistivity of the Winner AMN is located to the right of the center point, the weighted corrected resistivity of the Winner AMN is,

[0091]

[0092] When the resistivity of the Winner MNB is located to the left of the center point, the weighted corrected resistivity of the Winner MNB is,

[0093]

[0094] When the resistivity of the Winner MNB is located to the right of the center point, the weighted corrected resistivity of the Winner MNB is,

[0095]

[0096] The corrected Winner AMN resistivity and the corrected Winner MNB resistivity are,

[0097]

[0098]

[0099] Where n represents different depths at the same location of the dam, L represents the location to the left of the center point, and R represents the location to the right of the center point.

[0100] S103: Based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, obtain the initial model of the dam structure, and obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB.

[0101] Optionally, the resistivity values ​​of the dam body and foundation are measured in situ to obtain the statistical resistivity values ​​of the dam rock and soil at different locations. Based on the geological data of the reservoir, a geological structure model of the dam is constructed, and resistivity values ​​are assigned to the dam body and foundation, thus establishing an initial geological-resistivity dam structure model.

[0102] The weighted and corrected apparent resistivity of the Winner AMN device and the Winner MNB device were used to replace the original data, and the data from both devices were arranged and combined sequentially according to the format of the horizontal position of the recording point, the isolation coefficient, 1, and the resistivity value, thereby establishing a [database structure]. and Joint inversion documents. In the context of... and When performing the inversion process, the data volumes of the two devices should be inverted separately in sequence.

[0103] First, the resistivity of the Wenner AMN data is inverted. Based on the characteristics of the Wenner AMN resistivity data, the initial geological-resistivity dam structure model is automatically meshed. The electric field distribution in the initial geological-resistivity dam structure model is calculated using the finite element method. The calculated apparent resistivity at each recording point is then compared with... The root mean square error (RMSE) is calculated and used as the termination condition for the inversion. When the RMSE exceeds a preset value, the initial geological-resistivity dam structure model is modified and used as the reference model for the next inversion. Electric field calculations are then performed on the reference model, and the apparent resistivity at each recorded point is compared with... If the root mean square error between the two is greater than the preset value, the previous reference model is modified again and used as the reference model for the next inversion. The reference model is iterated and calculated continuously until the inversion model of the Wenner AMN device within the allowable error range is finally obtained.

[0104] Secondly, the resistivity of the Winner MNB was inverted. Based on the characteristics of the Winner MNB resistivity data, the Winner AMN inversion model was automatically meshed. The electric field distribution in the Winner AMN inversion model was calculated using the finite element method. The calculated apparent resistivity of each recording point was then compared with... The root mean square error (RMSE) is calculated and used as the termination condition for the inversion. If the RMSE exceeds a preset value, the initial dam structure model of the Wenner-AMN inversion model is modified and used as the reference model for the next inversion. Electric field calculations are then performed on the reference model, and the apparent resistivity at each recorded point is compared with... If the root mean square error between the two is greater than the preset value, the previous reference model is modified again and used as the reference model for the next inversion. The reference model is iterated continuously to obtain the inversion model of the Winner MNB device within the allowable error range. The inversion model of Winner MNB is used as the joint data inversion resistivity map of Winner AMN and Winner MNB.

[0105] S104: Based on the combined data inversion resistivity map of the Winner AMN and the Winner MNB, obtain the inversion resistivity of the Winner AMN and the Winner MNB at different locations at the same depth, and perform depth correction to obtain the spatial location of the dam leakage.

[0106] Optionally, S104 specifically includes:

[0107] Obtain the ratio of dam body depth to dam body resistivity at different locations within the same dam layer, and the ratio of dam shoulder depth to dam shoulder resistivity.

[0108] Based on the ratio of the dam body depth to the dam body resistivity and the ratio of the dam shoulder depth to the dam shoulder resistivity, the dam body depth correction coefficient and the dam shoulder depth correction coefficient are obtained.

[0109] The spatial location of the dam seepage can be obtained based on the dam body depth correction coefficient, the dam shoulder depth correction coefficient, and different depths.

[0110] Optionally, the ratio of the dam body depth to the dam body resistivity and the ratio of the dam abutment depth to the dam abutment resistivity are expressed by the following formula:

[0111]

[0112]

[0113] Wherein, η is the dam body depth correction coefficient, and ε is the dam abutment depth correction coefficient.

[0114] Based on the dam structure, the transverse resistivity profile of the dam is divided into two sections: the dam body and the dam abutment. The depth of the measuring points in each section is multiplied by η and ε to obtain a correction, thus obtaining the corrected resistivity profile. Based on the characteristics and depth of the resistivity profile, the spatial location of dam leakage can be inferred.

[0115] At least three boreholes shall be arranged, and at least one borehole shall be located in the dam body and abutment. The depth of the boreholes shall extend into the dam foundation.

[0116] In one specific embodiment, the reservoir's catchment area is approximately 0.12 km². 2The main length is approximately 0.3 km, and the total reservoir capacity is 250,000 m³. 3 Normal storage capacity is 20.0m³. 3 It is a small reservoir primarily used for irrigation and aquaculture, irrigating an area of ​​approximately 600 mu (about 40 hectares). The dam is 9 meters high, 50 meters long at the crest, and 3.0 meters wide at the crest. The upstream slope ratio is 1:2, with hexagonal precast concrete blocks forming the revetment. The downstream slope ratio is 1:2.2, with frame-beam turf forming the revetment. A total of 57 electrodes were laid out on the dam crest using electrical resistivity tomography (EMG), with an electrode spacing of 1 meter, running from the right bank to the left bank. AM data was collected on-site using a dual-mode parallel ERG, with a sampling interval of 0.05 seconds and a single power supply duration of 0.5 seconds. Boreholes ZK1, ZK2, and ZK3 correspond to depths of 44m, 25m, and 12.5m on the electrical resistivity survey line, respectively. The bedrock depth of borehole ZK1 is 6.5m, that of borehole ZK2 is 11.2m, and that of borehole ZK3 is 7.8m. Based on the borehole bedrock depths, the inverted resistivity profile depth of the right dam bedrock (0–15m) matches the borehole depth, therefore no correction is needed for this section. In the riverbed dam body (15–35m), the bedrock depth reflected in the resistivity profile at borehole ZK2 is 13.3m, resulting in a correction factor of 0.842 relative to the actual depth. In the left dam bedrock (35–55m), the bedrock depth reflected in the resistivity profile at borehole ZK3 is 5.5m, resulting in a correction factor of 1.418 relative to the actual depth.

[0117] Based on resistivity inversion, there are weak leakage zones on both sides of the dam.

[0118] This solution has the following prominent advantages:

[0119] (1) The present invention uses a dual-mode parallel electrical method to obtain the data volume of the entire electric field in one detection. Compared with the conventional high-density electrical method, it has high detection efficiency, rich data volume, and effectively suppresses the polarization problem of the power supply electrode.

[0120] (2) By optimizing and weighting the data volumes of the Winner AMN and Winner MNB devices, the asymmetry problem of conventional three-pole arrangement was avoided, and the optimization and inversion of a large amount of data improved the detection accuracy and resolution.

[0121] (3) By adopting the lateral correction method of the dam resistivity profile, the low resistivity anomalies in the dam body section and dam shoulder section were located in the depth direction, which improved the ability to distinguish the depth of the rock and soil layer and increased the accuracy of the target area for handling leakage hazards.

[0122] Example 3

[0123] Reference Figure 2The diagram shows a schematic of the structure of a device 30 for improving the accuracy of dual-mode parallel electrical resistivity tomography inversion of dam seepage, provided by an embodiment of the present invention.

[0124] The apparatus 30 provided in this embodiment of the invention for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage includes:

[0125] The first acquisition module 301 is used to acquire the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at a first depth, a second depth and a third depth, wherein the second depth is greater than the first depth and the third depth is greater than the second depth.

[0126] The second acquisition module 302 is used to perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth.

[0127] The third acquisition module 303 is used to obtain an initial model of the dam structure based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, and to obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB.

[0128] The correction module 304 is used to obtain the inversion resistivity of the Wenner AMN and the Wenner MNB at different locations at the same depth based on the joint data inversion resistivity map of the Wenner AMN and the joint data inversion resistivity map of the Wenner MNB, and to perform depth correction, thereby obtaining the spatial location of the dam leakage.

[0129] This invention extracts the joint resistivity values ​​of the Winner AMN and Winner MNB at equal intervals from the AM method data volume acquired by a dual-mode parallel resistivity transducer. Using the apparent resistivity at the center point of a single probe profile at different depths as a normalization parameter for the asymmetry of the apparent resistivity data, the Winner AMN and Winner MNB data volumes are normalized separately. All data are then used for forward and inverse calculations within a priori model of the existing dam structure. The inversion results are then processed at different depths, thereby improving inversion accuracy while reducing the error in probe depth.

[0130] The virtual system in this embodiment of the invention may be a device, or a component, integrated circuit, or chip in a terminal.

[0131] Furthermore, it should be noted that the device embodiments described above are merely illustrative and do not constitute a limitation on the scope of protection of the present invention. In practical applications, those skilled in the art can select some or all of the modules to achieve the purpose of the embodiment according to actual needs, and no restrictions are imposed here.

[0132] Example 4

[0133] This invention provides a device, a processor;

[0134] Memory used to store processor-executable instructions;

[0135] The processor is configured to invoke instructions stored in the memory to execute the method described in Embodiment 2.

[0136] This invention extracts the joint resistivity values ​​of the Winner AMN and Winner MNB at equal intervals from the AM method data volume acquired by a dual-mode parallel resistivity transducer. Using the apparent resistivity at the center point of a single probe profile at different depths as a normalization parameter for the asymmetry of the apparent resistivity data, the Winner AMN and Winner MNB data volumes are normalized separately. All data are then used for forward and inverse calculations within a priori model of the existing dam structure. The inversion results are then processed at different depths, thereby improving inversion accuracy while reducing the error in probe depth.

[0137] Example 5

[0138] This invention provides a computer-readable storage medium storing computer program instructions thereon, which, when executed by a processor, implement the method described in Embodiment 2.

[0139] This invention extracts the joint resistivity values ​​of the Winner AMN and Winner MNB at equal intervals from the AM method data volume acquired by a dual-mode parallel resistivity transducer. Using the apparent resistivity at the center point of a single probe profile at different depths as a normalization parameter for the asymmetry of the apparent resistivity data, the Winner AMN and Winner MNB data volumes are normalized separately. All data are then used for forward and inverse calculations within a priori model of the existing dam structure. The inversion results are then processed at different depths, thereby improving inversion accuracy while reducing the error in probe depth.

[0140] Computer-readable storage media can be tangible devices capable of holding and storing instructions for use by an instruction execution device. Computer-readable storage media can be, for example—but not limited to—electrical storage devices, magnetic storage devices, optical storage devices, electromagnetic storage devices, semiconductor storage devices, or any suitable combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), static random access memory (SRAM), portable compact disc read-only memory (CD-ROM), digital multifunction disc (DVD), memory sticks, floppy disks, mechanical encoding devices, such as punch cards or recessed protrusions storing instructions thereon, and any suitable combination thereof. The computer-readable storage media used herein are not to be construed as transient signals themselves, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through waveguides or other transmission media (e.g., light pulses through fiber optic cables), or electrical signals transmitted through wires.

[0141] The computer-readable program instructions described herein can be downloaded from computer-readable storage media to various computing / processing devices, or downloaded via a network, such as the Internet, local area network, wide area network, and / or wireless network, to an external computer or external storage device. The network may include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers, and / or edge servers. A network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards them to the computer-readable storage media in the respective computing / processing device.

[0142] The computer program instructions used to perform the operations of this invention may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-dependent instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Smalltalk, C++, etc., and conventional procedural programming languages ​​such as the "C" language or similar programming languages. The computer-readable program instructions may be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or may be connected to an external computer (e.g., via the Internet using an Internet service provider). In some embodiments, electronic circuitry, such as programmable logic circuitry, field-programmable gate arrays (FPGAs), or programmable logic arrays (PLAs), is personalized by utilizing state information from the computer-readable program instructions. This electronic circuitry can execute the computer-readable program instructions to implement various aspects of the invention.

[0143] Various aspects of the present invention are described herein with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer-readable program instructions.

[0144] These computer-readable program instructions can be provided to a processing unit of a general-purpose computer, a special-purpose computer, or other programmable data processing apparatus to produce a machine such that, when executed by the processing unit of the computer or other programmable data processing apparatus, they create means for implementing the functions / actions specified in one or more blocks of the flowchart and / or block diagram. These computer-readable program instructions can also be stored in a computer-readable storage medium that causes a computer, programmable data processing apparatus, and / or other device to operate in a particular manner. Thus, the computer-readable medium storing the instructions comprises an article of manufacture that includes instructions for implementing aspects of the functions / actions specified in one or more blocks of the flowchart and / or block diagram.

[0145] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions that execute on the computer, other programmable data processing apparatus, or other device to perform the functions / actions specified in one or more boxes of a flowchart and / or block diagram.

[0146] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of an instruction, which contains one or more executable instructions for implementing a specified logical function. In some alternative implementations, the functions marked in the blocks may occur in a different order than those shown in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, may be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.

[0147] Note that, unless otherwise explicitly stated, all features disclosed in this specification (including any appended claims, abstract, and drawings) may be replaced by alternative features for achieving the same, equivalent, or similar purpose. Therefore, unless explicitly stated otherwise, each disclosed feature is merely one example of a set of equivalent or similar features. Where used, "further," "preferably," "even further," and "more preferably" are simple starting points for describing another embodiment based on the foregoing embodiments, the combination of which with the foregoing embodiments constitutes the complete configuration of another embodiment. Any combination of several "further," "preferably," "even further," or "more preferably" settings following the same embodiment constitutes yet another embodiment.

[0148] Those skilled in the art should understand that the embodiments of the present invention described above and shown in the accompanying drawings are merely examples and do not limit the present invention. The objectives of the present invention have been fully and effectively achieved. The functions and structural principles of the present invention have been shown and explained in the embodiments, and any modifications or variations of the embodiments of the present invention may be made without departing from the stated principles.

[0149] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this disclosure, and are not intended to limit them. Although this disclosure has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of this disclosure.

Claims

1. A method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, characterized in that, The power supply electrodes include a first power supply electrode and a second power supply electrode, and the measuring electrodes include a first measuring electrode and a second measuring electrode; the distance between the first power supply electrode and the first measuring electrode, the distance between the first measuring electrode and the second measuring electrode, and the distance between the second measuring electrode and the second power supply electrode are all equal; Winner AMN represents the data between the first power supply electrode, the first measuring electrode, and the second measuring electrode; Winner MNB represents the data between the first measuring electrode, the second measuring electrode, and the second power supply electrode; the method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EPM) inversion of dam seepage includes: S101: Obtain the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at the first, second, and third depths, wherein the second depth is greater than the first depth and the third depth is greater than the second depth; S102: Perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth. S103: Based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, obtain the initial model of the dam structure, and obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB. S104: Based on the combined data inversion resistivity map of the Winner AMN and the combined data inversion resistivity map of the Winner MNB, obtain the inversion resistivity of the Winner AMN and the Winner MNB at different locations at the same depth, and perform depth correction to obtain the spatial location of the dam leakage. The weighted correction uses a formula. ； ； in, O represents the lateral center point of the dam; n represents different depths at the same location on the dam. The resistivity of the Wenner AMN at different depths at the same location on the dam; The resistivity of the Wenner MNB at different depths at the same location on the dam; When the resistivity of the Winner AMN is located to the left of the center point, the weighted corrected resistivity of the Winner AMN is, ； When the resistivity of the Winner AMN is located to the right of the center point, the weighted corrected resistivity of the Winner AMN is, ； When the resistivity of the Winner MNB is located to the left of the center point, the weighted corrected resistivity of the Winner MNB is, ； When the resistivity of the Winner MNB is located to the right of the center point, the weighted corrected resistivity of the Winner MNB is, ； The corrected Winner AMN resistivity and the corrected Winner MNB resistivity are, ； ； Where n represents different depths at the same location of the dam, L represents the location to the left of the center point, and R represents the location to the right of the center point.

2. The method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage according to claim 1, characterized in that, S104 specifically includes: Obtain the ratio of dam body depth to dam body resistivity at different locations within the same dam layer, and the ratio of dam shoulder depth to dam shoulder resistivity. Based on the ratio of the dam body depth to the dam body resistivity and the ratio of the dam shoulder depth to the dam shoulder resistivity, the dam body depth correction coefficient and the dam shoulder depth correction coefficient are obtained. The spatial location of the dam seepage can be obtained based on the dam body depth correction coefficient, the dam shoulder depth correction coefficient, and different depths.

3. The method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage according to claim 2, characterized in that, include: The ratio of the dam body depth to the dam body resistivity, and the ratio of the dam abutment depth to the dam abutment resistivity, are expressed by the following formula: ； Wherein, η is the dam body depth correction coefficient, and ε is the dam abutment depth correction coefficient.

4. A device for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage, characterized in that, Specifically, the method for improving the accuracy of dual-mode parallel electrical resistivity tomography (EDT) inversion of dam seepage as described in any one of claims 1 to 3 includes: The first acquisition module is used to acquire the resistivity of the Wenner AMN and the resistivity of the Wenner MNB at the same location on the dam at a first depth, a second depth, and a third depth, wherein the second depth is greater than the first depth, and the third depth is greater than the second depth. The second acquisition module is used to perform weighted correction on the resistivity of the Winner AMN and the resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth, and obtain the corrected resistivity of the Winner AMN and the corrected resistivity of the Winner MNB in ​​the first depth, the second depth and the third depth. The third acquisition module is used to obtain an initial model of the dam structure based on the corrected Winner AMN resistivity and the corrected Winner MNB resistivity, and to obtain the joint data inversion resistivity map of the Winner AMN and the joint data inversion resistivity map of the Winner MNB. The correction module is used to obtain the inversion resistivity of the Wenner AMN and the Wenner MNB at different locations at the same depth based on the joint data inversion resistivity map of the Wenner AMN and the joint data inversion resistivity map of the Wenner MNB, and to perform depth correction, thereby obtaining the spatial location of the dam leakage.

Images