CN117272213B - Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters - Google Patents
Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters Download PDFInfo
- Publication number
- CN117272213B CN117272213B CN202311551696.2A CN202311551696A CN117272213B CN 117272213 B CN117272213 B CN 117272213B CN 202311551696 A CN202311551696 A CN 202311551696A CN 117272213 B CN117272213 B CN 117272213B
- Authority
- CN
- China
- Prior art keywords
- parameter
- normalization
- parameters
- scanning
- correlation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000010408 sweeping Methods 0.000 title claims abstract description 23
- 239000003344 environmental pollutant Substances 0.000 title abstract description 24
- 231100000719 pollutant Toxicity 0.000 title abstract description 24
- 238000010606 normalization Methods 0.000 claims abstract description 170
- 230000002159 abnormal effect Effects 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims description 38
- 239000000356 contaminant Substances 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 14
- 238000010219 correlation analysis Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 230000005856 abnormality Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 241000282414 Homo sapiens Species 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000010287 polarization Effects 0.000 description 19
- 230000006870 function Effects 0.000 description 9
- 230000005284 excitation Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Abstract
The invention provides a method, a device, equipment and a medium for sweeping underground pollutant by using ground-based comprehensive parameters, wherein the method for sweeping underground pollutant by using ground-based comprehensive parameters comprises the following steps: acquiring various scanning data of each measuring point of a target area; normalizing the first scanning surface characteristics and the second scanning surface characteristics of all the measuring points; calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result; according to the correlation result, calculating the comprehensive parameters of each measuring point according to the first normalization parameter and the second normalization parameter; and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of the target area. The beneficial effects of the invention are as follows: the method improves the abnormal recognition precision of the polluted site, reduces the probability of human misjudgment, provides richer information support for the treatment of the polluted site, and improves the exploration effect and precision of underground pollutants.
Description
Technical Field
The invention relates to the technical field of the earth science, in particular to a method, a device, equipment and a medium for sweeping the ground physical comprehensive parameters of underground pollutants.
Background
Along with the increase of human pollutants, the influence on the living environment of human beings is larger and larger, how to better judge the distribution situation of the pollutants, and provide information with higher precision and higher accuracy for subsequent treatment, and further improvement on the existing methods of geology, geophysics and geochemistry is necessary, so that the exploration precision and accuracy of the existing methods are improved.
The existing geological, geophysical and geochemical scanning achievements generally adopt scanning data to carry out single-parameter mapping interpretation according to coordinates of measuring points, and though the abnormal condition of the single parameters can be well analyzed, the interrelationship among different parameters is difficult to judge; even if judgment among multiple parameters is carried out, the judgment is carried out manually through personal experience; because different individuals have different individual differences, interpretation and judgment results of the corresponding methods are often different from person to person, so that the results have stronger difference, randomness and individuation, comprehensive interpretation of the results is inconvenient, and misjudgment of judgment results is easy to cause.
Disclosure of Invention
The invention mainly aims to provide a method, a device, equipment and a medium for sweeping the ground physical comprehensive parameters of underground pollutants, which improve the abnormal recognition precision of a polluted site, reduce the probability of human misjudgment and provide richer information support for the treatment of the polluted site.
In one aspect, the invention provides a method for cleaning underground pollutants by using ground-based comprehensive parameters, comprising the following steps:
according to the parameter scanning request, scanning data of each measuring point in a target area are obtained, wherein the scanning data comprise a first scanning characteristic and a second scanning characteristic;
normalizing the first scanning surface features and the second scanning surface features of all measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface features, and the second normalization parameters are used for representing normalization processing results of the second scanning surface features;
calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, wherein the correlation result comprises one of positive correlation, negative correlation and uncorrelation;
according to the correlation result, calculating comprehensive parameters of each measuring point according to the first normalization parameters and the second normalization parameters, wherein the comprehensive parameters are used for representing whether the scan data of the measuring point is abnormal or not;
and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of a target area.
According to the method for scanning the underground pollutant by the ground physical comprehensive parameters, scanning data of each measuring point of a target area are acquired according to a parameter scanning request, and the method comprises the following steps:
and performing scanning processing on the target area by adopting a scanning method according to preset scanning grid parameters to obtain scanning data of each measuring point, wherein the scanning grid parameters comprise transverse adjacent point distances and longitudinal adjacent point distances, and the scanning method comprises at least one of geophysical scanning, geological scanning and chemical detection scanning.
The method for scanning the surface physical and chemical comprehensive parameters of the underground pollutant comprises the following steps of:
performing the following processing according to the relation between the number of first scanning surface features or the number of second scanning surface features in the scanning surface data of each measuring point and the total number of the measuring points;
when the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, adopting a normalization formula
Performing normalization processing on the first scan feature or the second scan feature, whereinRepresenting one of said first normalization parameter and said second normalization parameter,/or->Representing one of said first and said second glance characteristics, wherein ∈>Representing the plane coordinates of the measuring point, wherein ∈>Numbering the measuring points of the scanning surface, and，/>the total number of measuring points for the scanning surface, wherein +.>Is a natural number not less than 1, wherein +.>A natural number of not less than 2;
and when the number of the first scanning surface features or the number of the second scanning surface features is inconsistent with the total number of the measuring points, performing gridding processing on coordinates and scanning surface data according to the measuring points so that the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, and performing normalization processing on the first scanning surface features or the second scanning surface features by adopting a normalization formula, wherein grid parameters obtained by the gridding processing comprise an ordinate range, an abscissa range, an ordinate gridding interval and an abscissa gridding interval.
According to the method for sweeping the ground physical synthesis parameters of the underground pollutants, the correlation between the first normalization parameters and the second normalization parameters is calculated to obtain a correlation result, and the method comprises the following steps:
representing the said by means of a correlation parameterCorrelation of a first normalized parameter and said second normalized parameter, wherein the correlation parameterIs that
Wherein,indicating positive correlation +.>Indicating a negative correlation +.>Representing uncorrelation;
wherein,for the total number of measuring points participating in normalized parameter correlation analysis, +.>，For measuring point plane coordinates +.>One of said first normalized parameter and said second normalized parameter participating in a correlation analysis,/i->For measuring point plane coordinates +.>Is involved in correlation analysis and the second normalizationAnother one of the parameters->Parameter number for scan, and +.>，/>Parameter number for scan, and +.>，/> ，/>For the total number of parameters for acquiring scan data at scan points, wherein +.>A natural number of not less than 2, wherein
。
The method for sweeping the ground physical comprehensive parameters of the underground pollutant comprises the following steps:
wherein the method comprises the steps ofAnd->Are natural numbers of not less than 5.
According to the method for sweeping the ground physical feature comprehensive parameters of the underground pollutants, according to the correlation result, the comprehensive parameters of each measuring point are calculated according to the first normalized parameters and the second normalized parameters, and the method comprises the following steps:
the correlation result of the first normalization parameter and the second normalization parameter is positive correlation, and the method adopts
Calculating a composite parameter, whereinRepresenting the comprehensive parameters;
the correlation result of the first normalization parameter and the second normalization parameter is negative correlation, and the method adopts
Calculating comprehensive parameters;
and if the correlation result of the first normalization parameter and the second normalization parameter is uncorrelated, not executing the calculation of the comprehensive parameter.
The method for scanning the ground physical comprehensive parameters of the underground pollutant, wherein at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters is interpreted to obtain the parameter scanning of the target area, comprises the following steps:
and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters of all the measuring points, and displaying an interpretation result through a contour map.
The invention also discloses a ground physical comprehensive parameter scanning device of the underground pollutant, which comprises:
the first module is used for acquiring scanning data of each measuring point of the target area according to the parameter scanning request, wherein the scanning data comprises a first scanning characteristic and a second scanning characteristic;
the second module is used for carrying out normalization processing on the first scanning surface characteristics and the second scanning surface characteristics of all measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface characteristics, and the second normalization parameters are used for representing normalization processing results of the second scanning surface characteristics;
the third module is used for calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, and the correlation result comprises one of positive correlation, negative correlation and uncorrelation;
a fourth module, configured to calculate, according to the correlation result, a comprehensive parameter of each measurement point according to the first normalization parameter and the second normalization parameter, where the comprehensive parameter is used to characterize whether the scan data of the measurement point has an abnormality;
and a fifth module, configured to interpret at least one of the integrated parameter, the first normalized parameter and the second normalized parameter to obtain a parameter scan of the target area.
In another aspect, the invention provides an electronic device comprising a processor and a memory;
the memory is used for storing programs;
the processor executes the program to implement the method as described above.
The invention also discloses a computer program product or a computer program comprising computer instructions stored in a computer readable storage medium. The computer instructions may be read from a computer-readable storage medium by a processor of a computer device, and executed by the processor, cause the computer device to perform the method described previously.
The beneficial effects of the invention are as follows: parameters acquired based on a ground surface feature scanning mode are acquired by adopting different measuring points, and abnormal points are easier to identify through normalization and comprehensive association processing of different parameters, so that the abnormal identification precision of a polluted site can be improved, the probability of artificial misjudgment is reduced, richer information support is provided for the treatment of the polluted site, and the exploration effect and precision of underground pollutants are improved.
Drawings
The foregoing and/or additional aspects and advantages of the invention will become apparent and may be better understood from the following description of embodiments taken in conjunction with the accompanying drawings in which:
FIG. 1 is a schematic diagram of a geochemical integrated parameter sweep system for subsurface contaminants in accordance with an embodiment of the invention.
FIG. 2 is a flow chart of a method for surface characterization of subsurface contaminants, according to an embodiment of the invention.
FIG. 3 is a plan view of a working arrangement of a geochemical sweep of subsurface contaminants in accordance with an embodiment of the invention.
FIG. 4 is a contour plot of the polarization ratio of a laser plane obtained based on the station arrangement of FIG. 3, in accordance with an embodiment of the present invention.
FIG. 5 is a resistivity contour plot of a laser plane obtained based on the FIG. 3 station arrangement in accordance with an embodiment of the invention.
FIG. 6 is a contour plot of normalized parameters of the polarization of a laser plane obtained based on the FIG. 3 station arrangement in accordance with an embodiment of the present invention.
FIG. 7 is a contour plot of normalized parameters of resistivity of a laser plane obtained based on the FIG. 3 station arrangement, in accordance with an embodiment of the present invention.
FIG. 8 is a contour plot of the overall parameters of a laser scan obtained based on the FIG. 3 station arrangement in accordance with an embodiment of the present invention.
FIG. 9 is a schematic diagram of an apparatus for the geochemical synthesis parameter sweep of subsurface contaminants in accordance with an embodiment of the invention.
Detailed Description
Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. In the following description, suffixes such as "module", "part" or "unit" for representing elements are used only for facilitating the description of the present invention, and have no particular meaning in themselves. Thus, "module," "component," or "unit" may be used in combination. "first", "second", etc. are used for the purpose of distinguishing between technical features only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated. In the following description, the continuous reference numerals of the method steps are used for facilitating examination and understanding, and the technical effects achieved by the technical scheme of the invention are not affected by adjusting the implementation sequence among the steps in combination with the overall technical scheme of the invention and the logic relations among the steps. The embodiments described below by referring to the drawings are illustrative only and are not to be construed as limiting the invention.
Referring to fig. 1, fig. 1 is a schematic diagram of a geodesic integrated parameter scanning system of an underground pollutant according to an embodiment of the present invention, which includes a scanning device 100 and a client 200, wherein the scanning device 100 mainly includes a measuring point 110 and a scanning grid 120, and a user of the client 200 obtains scanning data of each measuring point in a target area, and the scanning data includes a first scanning feature and a second scanning feature; normalizing the first scanning surface characteristics and the second scanning surface characteristics of all the measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface characteristics, and the second normalization parameters are used for representing normalization processing results of the second scanning surface characteristics; calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, wherein the correlation result comprises one of positive correlation, negative correlation and uncorrelation; according to the correlation result, calculating the comprehensive parameters of each measuring point according to the first normalization parameter and the second normalization parameter, wherein the comprehensive parameters are used for representing whether the scan data of the measuring point is abnormal or not; and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of the target area, and displaying the parameter scan of the target area in a contour map mode. Referring to fig. 2, fig. 2 is a flowchart of a method for cleaning the surface of a subsurface contaminant according to an embodiment of the present invention, which includes, but is not limited to, steps S100 to S500:
s100, acquiring scanning data of each measuring point of a target area according to a parameter scanning request, wherein the scanning data comprises a first scanning feature and a second scanning feature.
In some embodiments, the first and second scan features are scan features obtained by different scan methods, such as polarization and resistivity.
In some embodiments, referring to the plan view of the working layout of the geochemical sweeping surface of the underground pollutant shown in fig. 3, the technical scheme of the embodiment of the invention obtains the plane coordinates of each measuring point by expanding the area of the geophysical excitation sweeping surface in the target area and designing the grid parameters of the sweeping surface according to the requirements of the working of the sweeping surface (as shown in fig. 3, 324 measuring points are arranged according to the grid parameters of 10 meters between the longitudinal adjacent points and 50 meters between the transverse adjacent points)(as shown in the vertical and horizontal coordinates in fig. 3); performing excitation sweeping work to obtain parameters of sweeping of each measuring point>(As shown in FIG. 3, table 1 and Table 2, the geophysical induced-polarization scanning was performed in this embodiment, and 2 scanning parameters were obtained, namely, the polarization ratio and the resistivity, respectively; and were setFor the susceptibility parameter of the scan, +.>Resistivity parameters for the sweep); wherein->Numbering measuring points for scanning, and；/>to develop the total number of stations of the scan, wherein +.>Is not less than 1Natural number (n=324 in this embodiment as shown in fig. 3); />Parameter number for scan, and +.>；/>For the total number of parameters acquired at a certain measuring point during scanning, wherein +.>For a natural number not less than 2 (as shown in fig. 3 and tables 1 and 2, m=2 in this embodiment, i.e. the present scan, 2 parameters are obtained at each measurement point, namely, a polarization rate parameter and a resistivity parameter); wherein parameter->Not all 0 (as shown in tables 1 and 2, in this embodiment, the polarization and resistivity parameters of the geophysical excitation plane are not 0).
In some embodiments, where geophysical electro-sweeping is only one of the sweeping methods, the sweeping method may also be a geological sweeping and a chemical detection sweeping.
And S200, carrying out normalization processing on the polarizability and the resistivity of all the measuring points to obtain a first normalization parameter and a second normalization parameter, wherein the first normalization parameter is used for representing the normalization processing result of the polarizability, and the second normalization parameter is used for representing the normalization processing result of the resistivity.
In some embodiments, the following process is performed according to the number of polarizability or resistivity in the scan data for each station versus the total number of stations:
when the number of the polarizability or the number of the resistivity is consistent with the total number of the measuring points, calculating by adopting a normalization formula; and when the number of the polarizability or the resistivity is inconsistent with the total number of the measuring points, performing gridding processing on the coordinate and the scan data according to the measuring points so that the number of the polarizability or the resistivity is consistent with the total number of the measuring points, and performing normalization processing on the polarizability or the resistivity by adopting a normalization formula, wherein grid parameters obtained by the gridding processing comprise an ordinate range, an abscissa range, an ordinate gridding interval and an abscissa gridding interval.
In some embodiments, the polarization rate and resistivity parameters of all the measuring points (such as 324 measuring points shown in fig. 3) obtained by the technical scheme of the corresponding embodiment in fig. 3 are normalized to obtain normalized parameters of the parameters(As shown in tables 1 and 2, the polarization ratio normalization parameter +.>And resistivity normalization parameter->) The method comprises the steps of carrying out a first treatment on the surface of the Wherein->Is the total number of certain parameters, wherein +.>A natural number of not less than 2; since the polarization and resistivity parameters are obtained at all the measurement points in the present embodiment, the polarization and resistivity normalization parameter calculation is performed by directly adopting the normalization formula (1) for the measured parameter values (the calculation results are shown in tables 1 and 2, the polarization normalization parameter +.>And resistivity normalization parameter->) In this case +.>(as shown in figures 3 and tables 1 and 2,）；
;
wherein, table 1 and table 2 respectively show the polarization rate, the resistivity and the polarization rate normalization parameters of the geophysical excitation plane on which the measuring points 1-10 and the measuring points 310-322 are arranged,Resistivity normalization parameters->Correlation coefficient, < >>Comprehensive parameters,A kind of electronic device
TABLE 1 parameter Table for measuring points 1-10 geophysical excitation scan
TABLE 2 parameter Table for measurement points 313-322 geophysical excitation scan
S300, calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, wherein the correlation result comprises one of positive correlation, negative correlation and uncorrelation.
In some embodiments, the correlation coefficient is obtained by adopting the formula (2) to analyze the correlation between one normalized parameter and another normalized parameter(correlation coefficient shown in Table 1 +.>A correlation coefficient of the polarization ratio normalization parameter and the resistivity normalization parameter); if the correlation coefficient->Then there is a positive correlation between the two normalized parameters; if the correlation coefficient->There is a negative correlation between the two normalized parameters (as shown in table 1, the correlation coefficient r=0.11 in the present embodiment, that is, the relationship between the polarization rate parameter and the resistivity parameter is a negative correlation); if the correlation coefficient->Then there is no correlation between the two normalized parameters; />Total number of measurement points (n=324 in this embodiment) for participation in normalized parameter correlation analysis;(n=324 in this embodiment, as shown in fig. 3); />For measuring point plane coordinates +.>One of the normalization parameters involved in the correlation analysis (as shown in Table 1 +.in this embodiment)>=1, i.e. the polarizability parameter); />For measuring point plane coordinates +.>Another normalized parameter involved in correlation analysis (as shown in Table 1, in this embodiment +.>=2, i.e. resistivity parameter); />Parameter number for scan, and +.>；/>Parameter number for scan, and +.>；/>；/>For the total number of parameters acquired at a certain measuring point during scanning, wherein +.>A natural number of not less than 2 (in this embodiment, m=2 as shown in fig. 3 and table 1); wherein the method comprises the steps of，；
S400, according to the correlation result, calculating the comprehensive parameters of each measuring point according to the first normalization parameter and the second normalization parameter, wherein the comprehensive parameters are used for representing whether the scan data of the measuring point is abnormal or not.
In some embodiments, if there is a positive correlation between the two normalized parameters, then calculate the combined parameters of the two normalized parameters for all points using equation (3)The method comprises the steps of carrying out a first treatment on the surface of the If the two normalization parameters have negative correlation, calculating the comprehensive parameters of the two normalization parameters of the same measuring point by adopting a formula (4)>If a certain measuring point in the formula (4)The comprehensive parameters of the measuring point are not calculated; if the two normalization parameters are not related, not performing comprehensive parameter calculation; (according to the table 1, the correlation coefficient of the polarization ratio normalization parameter and the resistivity normalization parameter is-0.11, and the resistivity normalization parameter of each measuring point +.>All are not 0, so in this embodiment, formula (4) is selected to perform the overall parameter +.>The results obtained are shown as "comprehensive parameters +.>"column; )
S500, interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of the target area.
Comprehensive parameters of all measuring pointsAnd/or normalizing parameters to obtain all the parameters in the regionA plurality of glance interpretation results; if the comprehensive parameters cannot be calculated at all the measuring points, analyzing and interpreting the normalized parameters of all the measuring points to obtain the scan interpretation result in the measuring area. (in this embodiment, all the measuring points can calculate the comprehensive parameters and the normalization parameters, so that the interpretation of the excitation scan results can be performed according to the comprehensive parameters and the normalization parameters in this embodiment; as shown in fig. 4, 5, 6, 7 and 8); the contour diagram of the comprehensive parameters obtained in fig. 8 can better embody the distribution condition of the detection targets because the polarizability and the resistivity show negative correlation characteristics, and the result (fig. 8) of the comprehensive parameters is obtained; fig. 7 better shows the strong polarization and low resistivity features of the present embodiment, and thus better finds anomalies of the strong polarization and low resistivity features.
FIG. 9 is a schematic diagram of an apparatus for the geochemical synthesis parameter sweep of subsurface contaminants in accordance with an embodiment of the invention. The apparatus includes a first module 910, a second module 920, a third module 930, a fourth module 940, and a fifth module 950.
The first module 910 is configured to obtain, according to a parameter scan request, scan data of each measurement point in the target area, where the scan data includes a first scan feature and a second scan feature; the second module 920 is configured to normalize the first scan characteristic and the second scan characteristic of all the measurement points to obtain a first normalization parameter and a second normalization parameter, where the first normalization parameter is used to represent a normalization result of the first scan characteristic, and the second normalization parameter is used to represent a normalization result of the second scan characteristic; a third module 930, configured to calculate a correlation between the first normalized parameter and the second normalized parameter to obtain a correlation result, where the correlation result includes one of a positive correlation, a negative correlation, and an uncorrelation; a fourth module 940, configured to calculate, according to the correlation result, a comprehensive parameter of each measurement point according to the first normalization parameter and the second normalization parameter, where the comprehensive parameter is used to characterize whether the scan data of the measurement point has an anomaly; a fifth module 950 is configured to interpret at least one of the integrated parameter, the first normalized parameter, and the second normalized parameter to obtain a parameter scan of the target area.
For example, in cooperation with the first module 910, the second module 920, the third module 930, the fourth module 940, and the fifth module 950 in the apparatus, the apparatus of the embodiment may implement any of the foregoing ground-feature integrated parameter scanning methods for underground pollutants, that is, obtain scanning data of each measuring point of the target area according to a parameter scanning request, where the scanning data includes a first scanning feature and a second scanning feature; normalizing the first scanning surface characteristics and the second scanning surface characteristics of all the measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface characteristics, and the second normalization parameters are used for representing normalization processing results of the second scanning surface characteristics; calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, wherein the correlation result comprises one of positive correlation, negative correlation and uncorrelation; according to the correlation result, calculating the comprehensive parameters of each measuring point according to the first normalization parameter and the second normalization parameter, wherein the comprehensive parameters are used for representing whether the scan data of the measuring point is abnormal or not; and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of the target area. The beneficial effects of the invention are as follows: the first scanning surface features and the second scanning surface features which are acquired based on the geophysical induced polarization scanning surface mode are acquired for different measuring points, and the abnormal points are easier to identify through normalization and comprehensive association processing of the first scanning surface features and the second scanning surface features, so that the abnormal identification precision of a polluted site can be improved, the probability of human misjudgment is reduced, richer information support is provided for the treatment of the polluted site, and the exploration effect and precision of underground pollutants are improved.
The embodiment of the invention also provides electronic equipment, which comprises a processor and a memory;
the memory stores a program;
the processor executes a program to execute the above-mentioned ground physical comprehensive parameter sweeping method of the underground pollutant; the electronic equipment has the functions of a software system for carrying and running the ground-based comprehensive parameter scanning of the underground pollutant, such as a computer, a speed regulator, a speed regulating device and the like.
Embodiments of the present invention also provide a computer readable storage medium storing a program for execution by a processor to implement a geochemical synthetic parameter sweeping method for subsurface contaminants as described above.
In some alternative embodiments, the functions/acts noted in the block diagrams may occur out of the order noted in the operational illustrations. For example, two blocks shown in succession may in fact be executed substantially concurrently or the blocks may sometimes be executed in the reverse order, depending upon the functionality/acts involved. Furthermore, the embodiments presented and described in the flowcharts of the present invention are provided by way of example in order to provide a more thorough understanding of the technology. The disclosed methods are not limited to the operations and logic flows presented herein. Alternative embodiments are contemplated in which the order of various operations is changed, and in which sub-operations described as part of a larger operation are performed independently.
Embodiments of the present invention also disclose a computer program product or computer program comprising computer instructions stored in a computer readable storage medium. The computer instructions may be read from a computer-readable storage medium by a processor of a computer device, and executed by the processor, to cause the computer device to perform the above-described geochemical synthetic parameter scanning method of subsurface contaminants.
Furthermore, while the invention is described in the context of functional modules, it should be appreciated that, unless otherwise indicated, one or more of the described functions and/or features may be integrated in a single physical device and/or software module or one or more functions and/or features may be implemented in separate physical devices or software modules. It will also be appreciated that a detailed discussion of the actual implementation of each module is not necessary to an understanding of the present invention. Rather, the actual implementation of the various functional modules in the apparatus disclosed herein will be apparent to those skilled in the art from consideration of their attributes, functions and internal relationships. Accordingly, one of ordinary skill in the art can implement the invention as set forth in the claims without undue experimentation. It is also to be understood that the specific concepts disclosed are merely illustrative and are not intended to be limiting upon the scope of the invention, which is to be defined in the appended claims and their full scope of equivalents.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Logic and/or steps represented in the flowcharts or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer readable medium may even be paper or other suitable medium on which the program is printed, as the program may be electronically captured, via, for instance, optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present invention have been shown and described, it will be understood by those of ordinary skill in the art that: many changes, modifications, substitutions and variations may be made to the embodiments without departing from the spirit and principles of the invention, the scope of which is defined by the claims and their equivalents.
While the preferred embodiment of the present invention has been described in detail, the present invention is not limited to the embodiments described above, and those skilled in the art can make various equivalent modifications or substitutions without departing from the spirit of the present invention, and these equivalent modifications or substitutions are included in the scope of the present invention as defined in the appended claims.
Claims (7)
1. A method for cleaning a surface of a subsurface contaminant with a ground-based synthetic parameter, comprising:
according to the parameter scanning request, scanning data of each measuring point in a target area are obtained, wherein the scanning data comprise a first scanning characteristic and a second scanning characteristic;
normalizing the first scanning surface features and the second scanning surface features of all measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface features, and the second normalization parameters are used for representing normalization processing results of the second scanning surface features;
the normalizing processing is performed on the first scanning surface characteristics and the second scanning surface characteristics of all measuring points, and the normalizing processing comprises the following steps:
performing the following processing according to the relation between the number of first scanning surface features or the number of second scanning surface features in the scanning surface data of each measuring point and the total number of the measuring points;
when the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, adopting a normalization formula
Performing on the first or second scan featureNormalization process in whichRepresenting one of said first normalization parameter and said second normalization parameter,/or->Representing one of said first and said second glance characteristics, wherein ∈>Representing the plane coordinates of the measuring point, wherein ∈>Parameter number for scan, and +.>Wherein->The scanning points are numbered, and +.>，/>The total number of measuring points for the scanning surface, wherein +.>A natural number of not less than 5;
when the number of the first scanning surface features or the number of the second scanning surface features is inconsistent with the total number of the measuring points, performing gridding processing on coordinates and scanning surface data according to the measuring points so that the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, and performing normalization processing on the first scanning surface features or the second scanning surface features by adopting a normalization formula, wherein gridding parameters comprise an ordinate range, an abscissa range, an ordinate gridding interval and an abscissa gridding interval;
calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, wherein the correlation result comprises one of positive correlation, negative correlation and uncorrelation;
the calculating the correlation between the first normalization parameter and the second normalization parameter to obtain a correlation result includes: representing the correlation of the first normalized parameter and the second normalized parameter by a correlation parameter, wherein the correlation parameterIs that
Wherein,indicating positive correlation +.>Indicating a negative correlation +.>Representing uncorrelation;
wherein,for the total number of measuring points participating in normalized parameter correlation analysis, +.>，Is the measuring pointThe plane coordinates are +.>One of said first normalized parameter and said second normalized parameter participating in a correlation analysis,/i->For measuring point plane coordinates +.>Is the other of said first normalized parameter and said second normalized parameter participating in the correlation analysis,/->Numbering parameters of the scan, and，/>parameter number for scan, and +.>，/>，/>For the total number of parameters for acquiring scan data at scan points, wherein +.>A natural number of not less than 2, wherein
According to the correlation result, calculating comprehensive parameters of each measuring point according to the first normalization parameters and the second normalization parameters, wherein the comprehensive parameters are used for representing whether the scan data of the measuring point is abnormal or not;
the calculating, according to the correlation result and the first normalization parameter and the second normalization parameter, the comprehensive parameter of each measuring point includes:
the correlation result of the first normalization parameter and the second normalization parameter is positive correlation, and the method adopts
Calculating a composite parameter, whereinRepresenting the comprehensive parameters;
the correlation result of the first normalization parameter and the second normalization parameter is negative correlation, and the method adopts
Calculating comprehensive parameters;
if the correlation results of the first normalization parameter and the second normalization parameter are uncorrelated, not executing the calculation of the comprehensive parameter;
and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters to obtain a parameter scan of a target area.
2. The method for geodesic integrated parameter scanning of subsurface contaminants according to claim 1, wherein said obtaining scan data for each measurement point of the target zone based on the parameter scan request comprises:
and performing scanning processing on the target area by adopting a scanning method according to preset scanning grid parameters to obtain scanning data of each measuring point, wherein the scanning grid parameters comprise transverse adjacent point distances and longitudinal adjacent point distances, and the scanning method comprises at least one of geophysical scanning, geological scanning and chemical detection scanning.
3. The method for geochemical integrated parameter scanning of subsurface contaminants according to claim 1, further comprising:
wherein N and N are natural numbers of not less than 5.
4. The method of claim 1, wherein interpreting at least one of the synthetic parameters, the first normalized parameters, and the second normalized parameters to obtain a parameter scan of a target zone comprises:
and interpreting at least one of the comprehensive parameters, the first normalization parameters and the second normalization parameters of all the measuring points, and displaying an interpretation result through a contour map.
5. A geodesic integrated parameter scanning device for subsurface contaminants, comprising:
the first module is used for acquiring scanning data of each measuring point of the target area according to the parameter scanning request, wherein the scanning data comprises a first scanning characteristic and a second scanning characteristic;
the second module is used for carrying out normalization processing on the first scanning surface characteristics and the second scanning surface characteristics of all measuring points to obtain first normalization parameters and second normalization parameters, wherein the first normalization parameters are used for representing normalization processing results of the first scanning surface characteristics, and the second normalization parameters are used for representing normalization processing results of the second scanning surface characteristics;
the second module is further used for executing the following processing according to the relation between the number of first scanning surface features or the number of second scanning surface features and the total number of the measuring points in the scanning surface data of each measuring point;
when the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, adopting a normalization formula
Performing normalization processing on the first scan feature or the second scan feature, whereinRepresenting one of said first normalization parameter and said second normalization parameter,/or->Representing one of said first and said second glance characteristics, wherein ∈>Representing the plane coordinates of the measuring point, wherein ∈>Numbering parameters of the scan, andwherein->The scanning points are numbered, and +.>，/>The total number of measuring points for the scanning surface, wherein +.>A natural number of not less than 5;
when the number of the first scanning surface features or the number of the second scanning surface features is inconsistent with the total number of the measuring points, performing gridding processing on coordinates and scanning surface data according to the measuring points so that the number of the first scanning surface features or the number of the second scanning surface features is consistent with the total number of the measuring points, and performing normalization processing on the first scanning surface features or the second scanning surface features by adopting a normalization formula, wherein gridding parameters comprise an ordinate range, an abscissa range, an ordinate gridding interval and an abscissa gridding interval;
the third module is used for calculating the correlation of the first normalization parameter and the second normalization parameter to obtain a correlation result, and the correlation result comprises one of positive correlation, negative correlation and uncorrelation;
the third module is further configured to calculate a correlation between the first normalized parameter and the second normalized parameter, to obtain a correlation result, where the calculating includes: representing the correlation of the first normalized parameter and the second normalized parameter by a correlation parameter, wherein the correlation parameterIs that
Wherein,indicating positive correlation +.>Indicating a negative correlation +.>Representing uncorrelation;
wherein,for the total number of measuring points participating in normalized parameter correlation analysis, +.>，For measuring point plane coordinates +.>One of said first normalized parameter and said second normalized parameter participating in a correlation analysis,/i->For measuring point plane coordinates +.>Is the other of said first normalized parameter and said second normalized parameter participating in the correlation analysis,/->Numbering parameters of the scan, and，/>parameter number for scan, and +.>，/>，/>For the total number of parameters for acquiring scan data at scan points, wherein +.>A natural number of not less than 2, wherein
A fourth module, configured to calculate, according to the correlation result, a comprehensive parameter of each measurement point according to the first normalization parameter and the second normalization parameter, where the comprehensive parameter is used to characterize whether the scan data of the measurement point has an abnormality;
the fourth module is further configured to calculate, according to the correlation result and the first normalization parameter and the second normalization parameter, a comprehensive parameter of each measurement point, where the calculation includes:
the correlation result of the first normalization parameter and the second normalization parameter is positive correlation, and the method adopts
Calculating a composite parameter, whereinRepresenting the comprehensive parameters;
the correlation result of the first normalization parameter and the second normalization parameter is negative correlation, and the method adopts
Calculating comprehensive parameters;
if the correlation results of the first normalization parameter and the second normalization parameter are uncorrelated, not executing the calculation of the comprehensive parameter;
and a fifth module, configured to interpret at least one of the integrated parameter, the first normalized parameter and the second normalized parameter to obtain a parameter scan of the target area.
6. An electronic device comprising a processor and a memory;
the memory is used for storing programs;
the processor executing the program to implement the geochemical integrated parameter sweeping method of subsurface contaminants according to any one of claims 1-4.
7. A computer readable storage medium, wherein the storage medium stores a program for execution by a processor to implement the geochemical synthesis parameter scanning method of subsurface contaminants according to any one of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311551696.2A CN117272213B (en) | 2023-11-21 | 2023-11-21 | Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311551696.2A CN117272213B (en) | 2023-11-21 | 2023-11-21 | Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117272213A CN117272213A (en) | 2023-12-22 |
| CN117272213B true CN117272213B (en) | 2024-02-02 |
Family
ID=89210864
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311551696.2A Active CN117272213B (en) | 2023-11-21 | 2023-11-21 | Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117272213B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1701336A (en) * | 2002-09-23 | 2005-11-23 | 哥伦比亚技术公司 | Smart data subsurface data repository system, method and computer program product |
| CN108323295A (en) * | 2017-12-05 | 2018-07-27 | 江苏大学 | A kind of seedling stage crop liquid manure based on multiple dimensioned habitat information detects and controls method and device |
| CN108445156A (en) * | 2018-03-19 | 2018-08-24 | 攀钢集团攀枝花钢铁研究院有限公司 | The identification method of high-titanium blast furnace slag object phase and application |
| CN108845361A (en) * | 2018-06-06 | 2018-11-20 | 中南大学 | A kind of sharp electric surface sweeping method that multiple supply terminals are powered simultaneously |
| CN112965141A (en) * | 2021-02-06 | 2021-06-15 | 核工业北京地质研究院 | Delineation method for favorable section of uranium polymetallic ore |
| CN113156516A (en) * | 2021-04-21 | 2021-07-23 | 中南大学 | New induced polarization method adopting accurate polarization constant |
| CN115640995A (en) * | 2022-09-30 | 2023-01-24 | 长江地球物理探测（武汉）有限公司 | Tunnel advanced geological forecast risk evaluation method and device |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100164513A1 (en) * | 2006-08-21 | 2010-07-01 | Uri Rapoport | non-destructive on-line method for measuring predetermined physical, electrochemical, chemical or biological state transformation of a substance and a system thereof |
| US9286419B2 (en) * | 2010-07-19 | 2016-03-15 | Terje Graham Vold | Computer simulation of electromagnetic fields |
| US20190018041A1 (en) * | 2013-11-06 | 2019-01-17 | Yeda Research And Development Co. Ltd. | Nanotube based transistor structure, method of fabrication and uses thereof |
| WO2021129965A1 (en) * | 2019-12-24 | 2021-07-01 | Navix International Limited | System, method and accesories for dielectric-mapping |
| CN111812724B (en) * | 2020-06-29 | 2021-08-27 | 山东大学 | Hidden karst pipeline detection method and system |
| CN114384593A (en) * | 2022-01-18 | 2022-04-22 | 中国地质科学院 | Distributed three-dimensional induced polarization data acquisition and processing device and method |
-
2023
- 2023-11-21 CN CN202311551696.2A patent/CN117272213B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1701336A (en) * | 2002-09-23 | 2005-11-23 | 哥伦比亚技术公司 | Smart data subsurface data repository system, method and computer program product |
| CN108323295A (en) * | 2017-12-05 | 2018-07-27 | 江苏大学 | A kind of seedling stage crop liquid manure based on multiple dimensioned habitat information detects and controls method and device |
| CN108445156A (en) * | 2018-03-19 | 2018-08-24 | 攀钢集团攀枝花钢铁研究院有限公司 | The identification method of high-titanium blast furnace slag object phase and application |
| CN108845361A (en) * | 2018-06-06 | 2018-11-20 | 中南大学 | A kind of sharp electric surface sweeping method that multiple supply terminals are powered simultaneously |
| CN112965141A (en) * | 2021-02-06 | 2021-06-15 | 核工业北京地质研究院 | Delineation method for favorable section of uranium polymetallic ore |
| CN113156516A (en) * | 2021-04-21 | 2021-07-23 | 中南大学 | New induced polarization method adopting accurate polarization constant |
| CN115640995A (en) * | 2022-09-30 | 2023-01-24 | 长江地球物理探测（武汉）有限公司 | Tunnel advanced geological forecast risk evaluation method and device |
Non-Patent Citations (4)
| Title |
|---|
| Identification of karst sinkholes in a forested karst landscape using airborne laser scanning data and water flow analysis;Jaroslav Hofierka等;《Geomorphology》;265-277 * |
| 大功率激电法在某地区金属矿普查中的应用;蔡伟涛;;世界有色金属(第02期);135-136 * |
| 激发极化法在有色金属矿产勘查中的研究进展;刘春明等;《中国有色金属学报》;203-221 * |
| 激电中梯装置测量长方异常体的三维正演模拟及启示;曹创华;邓专;康方平;柳建新;童孝忠;;物探化探计算技术(第04期);51-59 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117272213A (en) | 2023-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105221133B (en) | A kind of method and apparatus that content of organic carbon of hydrocarbon source rock is determined based on well logging multi-parameter | |
| Wulder et al. | Local spatial autocorrelation characteristics of remotely sensed imagery assessed with the Getis statistic | |
| Heil et al. | Comparison of the EM38 and EM38-MK2 electromagnetic induction-based sensors for spatial soil analysis at field scale | |
| MXPA04001369A (en) | Interpretation methods for nmr diffusion-t2 maps. | |
| CN105093299A (en) | Observation system optimization method based on offset vector tile technology and apparatus thereof | |
| NO322790B1 (en) | Method and apparatus for interpreting formation interfaces in geological formations using source log curves | |
| Schmidt et al. | Dreaming of perfect data: Characterizing noise in archaeo-geophysical measurements | |
| CN206378448U (en) | Polluted Soil detection arrangement of measuring-line structure based on comprehensive geophysical prospecting methods | |
| Pearce et al. | Numerical analysis of fold curvature using data acquired by high-precision GPS | |
| Sokolov et al. | Prospecting and evaluation of underground massive ice by ground-penetrating radar | |
| Karatson et al. | Preferred clast orientation in volcaniclastic mass-flow deposits: application of a new photo-statistical method | |
| CN112946751B (en) | Method and system for predicting plane distribution of different reservoirs | |
| CN107085240A (en) | A kind of side slope magnetic fluid detection system and method | |
| Abdelrahman et al. | A least-squares approach to depth determination from numerical horizontal self-potential gradients | |
| CN117272213B (en) | Method, device, equipment and medium for sweeping underground pollutant by using ground physical comprehensive parameters | |
| Green | Magnetic profile analysis | |
| CN110231661B (en) | Rock specimen and corresponding method for field exploration measurement of induced polarization parameters | |
| Pereira et al. | Assessing dataset equivalence and leveling data in geochemical mapping | |
| CN115420778A (en) | Method and system for quickly investigating and positioning shallow soil pollution | |
| CN206348292U (en) | Polluted Soil detection arrangement of measuring-line structure based on three-dimensional high-density resistivity method | |
| CN111221043B (en) | Working parameter optimization method of passive source electric field method | |
| CN113216945A (en) | Permeability quantitative evaluation method for tight sandstone reservoir | |
| CN109633760B (en) | Underground fluid monitoring method based on natural potential equivalent fusion imaging | |
| Gobashy et al. | High resolution ground magnetic survey (HRGM) for determining the optimum location of subsurface dam in Wadi Nu’man, Makkah Al Mukarammah, KSA | |
| US20090265109A1 (en) | Method for acquiring and processing magnetometric data by local updates in real time |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |