System and method for analyzing soil pollutants in land resources reutilization of abandoned mines

Abstract

The application provides a waste mine land resource recycling soil pollutant analysis system and method, identifies and circumscribes a plurality of apparent resistivity anomaly zones caused by pollutants as suspected pollution path zones, extracts apparent resistivity values in the spatial range of each suspected pollution path zone and characteristic pollutant concentration values of corresponding positions and corresponding depth intervals, and then obtains a physical-chemical response synergy degree of each suspected pollution path zone representing the spatial coupling degree of physical-chemical anomalies, determines the dominant channel index of each suspected pollution path zone according to the physical-chemical response synergy degree of each suspected pollution path zone and the extension of the suspected pollution path zone along the direction of potential underground water flow on the profile, judges the suspected pollution path zone with the highest dominant channel index as the dominant migration channel of the pollutants, and labels in combination with the geological profile. The scheme of the application can analyze the dominant migration path of underground pollutants from the perspective of spatial coupling of physical fields and chemical fields.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of soil pollutant analysis, in particular, the present application relates to a soil pollutant analysis system and method for abandoned mine land resource recycling. BACKGROUND

[0002] With the increasing emphasis on ecological environment protection in China, the pollution treatment and land resource safe recycling of a large number of historical abandoned mines have become an urgent task. Among them, accurately identifying the dominant migration path of underground pollutants such as heavy metals, acid wastewater and organic agents from the pollution source to the downstream soil and groundwater environment is the premise and key to implementing efficient and accurate monitoring and repair engineering.

[0003] Currently, the exploration methods for such problems mainly rely on single technical means or simple technical combination. Common methods include geophysical exploration methods and geochemical exploration methods. In recent years, the industry has also tried to combine geophysical and geochemical methods, but usually stays at the level of "working first and then manually comparing", that is, a small number of verification drill holes are laid out in geophysical anomaly areas. This combination lacks systematic data fusion and quantitative analysis logic, and has obvious defects, that is: first, the selection of verification points is subjective and may miss key areas; second, it cannot uniformly and quantitatively prioritize multiple anomaly areas; third, it is difficult to scientifically reveal the dominant migration rule of pollutants, resulting in insufficient targeting of subsequent treatment engineering such as monitoring well layout and barrier wall construction, causing resource waste and ineffective treatment. Therefore, how to analyze the dominant migration path of underground pollutants from the perspective of spatial coupling of physical and chemical fields, and thus provide a targeted position for setting up groundwater monitoring wells or implementing in-situ barrier repair engineering has become a difficult problem in the industry. SUMMARY

[0004] Based on this, the present application provides a soil pollutant analysis system and method for abandoned mine land resource recycling for analyzing the dominant migration path of underground pollutants from the perspective of spatial coupling of physical and chemical fields.

[0005] In a first aspect, the present application provides a soil pollutant analysis method for abandoned mine land resource recycling for analyzing the dominant migration path of pollutants in an abandoned mine. The method comprises the following steps:

[0006] According to the stratum lithology data of the abandoned mine, a plurality of resistivity anomaly areas caused by pollutants are identified and circled as suspected pollution path areas;

[0007] extracting apparent resistivity values within a spatial range of each suspected pollution path area and concentration values of a characteristic pollutant corresponding to a position and a corresponding depth interval, and then analyzing a correlation of spatial distribution trends of both to obtain a degree of coupling of a geochemical response of each suspected pollution path area;

[0008] determining an advantage channel index of each suspected pollution path area according to the degree of coupling of the geochemical response of each suspected pollution path area and an extension of each suspected pollution path area along a direction of a potential groundwater flow on a profile;

[0009] determining a pollutant advantage migration channel as the suspected pollution path area with the highest advantage channel index, and labeling the suspected pollution path area in combination with a geological profile to provide a target position for setting a groundwater monitoring well or implementing an in-situ barrier repair project.

[0010] In some embodiments, the identifying and delineating a plurality of apparent resistivity anomaly areas caused by pollutants as suspected pollution path areas according to stratum lithology data of the abandoned mine specifically includes:

[0011] obtaining apparent resistivity profile data and stratum lithology data of the abandoned mine;

[0012] performing anomaly detection on the apparent resistivity profile data to extract apparent resistivity anomaly areas;

[0013] performing cause filtering on the extracted apparent resistivity anomaly areas according to the stratum lithology data, and then delineating anomaly areas unrelated to non-pollution geological bodies as suspected pollution path areas.

[0014] In some embodiments, the performing anomaly detection on the apparent resistivity profile data to extract apparent resistivity anomaly areas specifically includes:

[0015] determining a background resistivity statistical characteristic value of the apparent resistivity profile data;

[0016] setting anomaly determination threshold values of relatively high resistance and relatively low resistance according to the background resistivity statistical characteristic value;

[0017] performing image segmentation on the apparent resistivity profile data based on the anomaly determination threshold values to extract and delineate independent apparent resistivity anomaly areas.

[0018] In some embodiments, the extracting apparent resistivity values within a spatial range of each suspected pollution path area and concentration values of a characteristic pollutant corresponding to a position and a corresponding depth interval specifically includes:

[0019] For each suspected pollution path area, performing geochemical sampling and detection in the suspected pollution path area and adjacent areas to obtain concentration data of a characteristic pollutant in the suspected pollution path area;

[0020] aligning the concentration data of the characteristic contaminant and the apparent resistivity profile data in spatial position and depth interval, to construct a spatially aligned chemical anomaly concentration profile;

[0021] extracting the aligned apparent resistivity values and corresponding concentration values of the characteristic contaminant from the spatial boundary range of the suspected pollution path area, to form a data pair set for subsequent analysis of correlation.

[0022] In some embodiments, aligning the concentration data of the characteristic contaminant and the apparent resistivity profile data in spatial position and depth interval, to construct a spatially aligned chemical anomaly concentration profile specifically comprises:

[0023] unifying the spatial coordinate system and elevation datum adopted by the concentration data of the characteristic contaminant and the apparent resistivity profile data;

[0024] based on the concentration data of the characteristic contaminant in the unified spatial coordinate system, reconstructing a continuous two-dimensional concentration distribution profile using a spatial interpolation algorithm;

[0025] aligning the reconstructed concentration distribution profile and the apparent resistivity profile at the grid nodes to generate the spatially aligned chemical anomaly concentration profile.

[0026] In some embodiments, analyzing the correlation between the spatial distribution trends to obtain the physico-chemical response synergy degree of each suspected pollution path area representing the degree of spatial coupling of the physico-chemical anomalies specifically comprises:

[0027] for each suspected pollution path area, obtaining a set of apparent resistivity-concentration data pairs of the suspected pollution path area;

[0028] performing rank conversion on the data in the set of apparent resistivity-concentration data pairs to eliminate dimensional differences and convert them into sequences suitable for non-parametric statistical analysis;

[0029] based on the converted sequences, determining a non-parametric statistical quantity that can represent the monotonic correlation between the apparent resistivity values and the concentration values of the characteristic contaminant;

[0030] taking the non-parametric statistical quantity as the physico-chemical response synergy degree of the suspected pollution path area.

[0031] In some embodiments, determining the dominant channel index of each suspected pollution path area according to the physico-chemical response synergy degree of the suspected pollution path area and its extension along the direction of the potential underground water flow on the profile specifically comprises:

[0032] determining a unified potential underground water flow direction in the abandoned mine based on topographic and geological structure data;

[0033] For each suspected pollution path area, a projection length of the suspected pollution path area along a direction of the potential underground water flow on a profile is determined as a quantitative indicator of the extension of the suspected pollution path area;

[0034] The physicochemical response synergy degree and the projection length are normalized respectively to obtain corresponding standardized evaluation values;

[0035] The standardized evaluation value of the physicochemical response synergy degree and the standardized evaluation value of the projection length are weighted and summed to obtain an advantage channel index of the suspected pollution path area.

[0036] In a second aspect, the present application provides a soil pollutant analysis system for reusing land resources of abandoned mines, which comprises:

[0037] An identification module is configured to identify and delineate a plurality of apparent resistivity anomaly areas caused by pollutants as suspected pollution path areas according to stratum lithology data of the abandoned mine;

[0038] A processing module is configured to extract apparent resistivity values in a spatial range of each suspected pollution path area and characteristic pollutant concentration values of corresponding positions and corresponding depth intervals, and then analyze the correlation of spatial distribution trends of the apparent resistivity values and the characteristic pollutant concentration values to obtain a physicochemical response synergy degree of each suspected pollution path area representing a spatial coupling degree of physicochemical anomalies;

[0039] The processing module is further configured to determine an advantage channel index of each suspected pollution path area according to the physicochemical response synergy degree of the suspected pollution path area and the extension of the suspected pollution path area along a direction of the potential underground water flow on a profile;

[0040] An execution module is configured to determine a suspected pollution path area with the highest advantage channel index as a pollutant advantage migration channel, and mark the suspected pollution path area in combination with a geological profile to provide a target position for setting a groundwater monitoring well or implementing an in-situ barrier repair project.

[0041] In a third aspect, the present application provides a computer device comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements steps of the soil pollutant analysis method for reusing land resources of abandoned mines when executing the computer program.

[0042] In a fourth aspect, the present application provides a computer readable storage medium storing a computer program, wherein the computer program is executed by a processor to implement steps of the soil pollutant analysis method for reusing land resources of abandoned mines.

[0043] The technical scheme provided by the embodiments of the present application has the following beneficial effects:

[0044] The abandoned mine land resource recycling soil pollutant analysis system and method provided by the application, first, according to the stratum lithology data of the abandoned mine, a plurality of apparent resistivity anomaly areas caused by pollutants are identified and circled as suspected pollution path areas, this step can automatically and objectively screen out the key target area related to pollution from the wide area geophysical scanning, and filter through geological knowledge to significantly reduce the "false anomaly" interference caused by pure water-bearing layer, clay layer and other non-pollution geological bodies, thereby improving the pertinence and starting point reliability of the subsequent exploration work, and avoiding invalid investment; secondly, the apparent resistivity value in the spatial range of each suspected pollution path area and the characteristic pollutant concentration value of the corresponding position and corresponding depth interval are extracted, and then the correlation of the spatial distribution trend of the two is analyzed, and the materialization response synergy degree of each suspected pollution path area representing the spatial coupling degree of the materialization anomaly is obtained, this step can realize the accurate fusion and quantitative correlation analysis of the geophysical field and the geochemical field in space, through the calculation of the "materialization response synergy degree" which is independent of the data distribution, the spatial trend consistency of the two kinds of data is converted into comparable quantitative evidence, thereby improving the objectivity and scientificity of the pollution cause identification, and overcoming the defect of strong multiple solutions of single data; then, the advantage channel index of each suspected pollution path area is determined according to the materialization response synergy degree of each suspected pollution path area and the extension along the potential groundwater flow direction on the profile, this step can comprehensively evaluate the potential channel from the two dimensions of "abnormal coupling strength" and "hydrogeological function", and combine the materialization response synergy degree and the extension along the groundwater flow direction into a unified "advantage channel index", thereby improving the logical rationality of the priority ranking of the plurality of suspected areas and the decision support intensity, and making the channel identification more in line with the physical mechanism of pollutant migration; finally, the suspected pollution path area with the highest advantage channel index is determined as the pollutant advantage migration channel, and is marked in combination with the geological profile, thereby providing a target position for setting an underground water monitoring well or implementing an in-situ barrier repair project, this step can directly convert the quantitative analysis result into an intuitive and operable engineering guidance drawing, and through the accurate marking of the advantage migration channel and the recommended engineering target position on the geological profile, thereby improving the targeting, one-time success rate and treatment benefit of the subsequent underground water monitoring well layout or in-situ barrier repair project, and realizing the close connection from technical detection to engineering implementation; in summary, the scheme of the application can analyze the advantage migration path of the underground pollutant from the perspective of the spatial coupling of the physical field and the chemical field. BRIEF DESCRIPTION OF DRAWINGS

[0045] Figure 1 is an example flow chart of the abandoned mine land resource recycling soil pollutant analysis method according to some embodiments of the application;

[0046] Figure 2 is an application scenario schematic diagram of the pollutant advantage migration path analysis data processing system according to some embodiments of the application;

[0047] Figure 3 This is a flowchart illustrating the determination of the degree of synergy of physical response according to some embodiments of this application;

[0048] Figure 4 This is a schematic diagram of the structure of a soil pollutant analysis system for the reuse of abandoned mine land resources, as shown in some embodiments of this application;

[0049] Figure 5 This is a schematic diagram of the structure of a computer device for implementing a method for analyzing soil pollutants in the reuse of abandoned mine land resources, according to some embodiments of this application. Detailed Implementation

[0050] To better understand the above technical solutions, the following will provide a detailed explanation of the technical solutions in conjunction with the accompanying drawings and specific implementation methods.

[0051] refer to Figure 1 The figure is an exemplary flowchart of a method for analyzing soil pollutants in the reuse of abandoned mine land resources according to some embodiments of this application. The method for analyzing soil pollutants in the reuse of abandoned mine land resources mainly includes the following steps:

[0052] In step 101, based on the stratigraphic lithology data of the abandoned mine, multiple areas of apparent resistivity anomalies caused by pollutants are identified and delineated as suspected pollution pathway areas.

[0053] In some embodiments, identifying and delineating multiple areas of apparent resistivity anomalies caused by pollutants as suspected pollution pathways based on stratigraphic lithology data of abandoned mines can be achieved through the following steps:

[0054] Obtain apparent resistivity profile data and stratigraphic lithology data of abandoned mines;

[0055] Anomaly detection is performed on the apparent resistivity profile data to extract the apparent resistivity anomaly areas;

[0056] Based on the stratigraphic lithology data, the extracted apparent resistivity anomaly areas are filtered for their origins, and then anomaly areas unrelated to non-contaminated geological bodies are delineated as suspected contamination pathway areas.

[0057] It should be noted that the apparent resistivity profile data in the present application refers to the image data set obtained by high-density resistivity exploration and inversion calculation, which systematically reflects the spatial continuous change of the underground medium conductivity in the form of two-dimensional profile. The role is to provide the basic geophysical field information covering the investigation area and revealing the heterogeneous underground electrical structure for the whole method; the stratum lithology data is the information describing the spatial distribution and lithology characteristics of different stratum units in the investigation area, which provides the key geological background knowledge and constraint conditions for distinguishing the causes of geophysical anomalies, and is an important basis for judging whether the anomaly is caused by non-polluted geological body, thus for preliminary screening and true existence filtering; the apparent resistivity anomaly area refers to the continuous spatial area whose resistivity value is significantly different from the surrounding background medium identified from the apparent resistivity profile data; the suspected pollution path area refers to the underground area where the pollution may exist or migrate, which is screened out by preliminary geophysical scanning and geological common sense filtering and needs to be further verified by geochemical means.

[0058] In specific implementation, the apparent resistivity profile data and stratum lithology data of the abandoned mine can be realized in the following manner, that is, based on the topography and hydrogeological conditions of the abandoned mine, high-density resistivity method measuring lines are laid out, multi-electrode measurement is used to collect the apparent resistivity raw data of the underground medium, the apparent resistivity raw data is preprocessed to eliminate distortion points and terrain correction is performed, then a resistivity inversion algorithm such as smooth constrained least squares inversion is used for calculation to obtain two-dimensional apparent resistivity profile data reflecting the real electrical structure of the underground; at the same time, the existing geological exploration report, drilling columnar graph and hydrogeological map of the region are collected, from which the stratum lithology distribution map is extracted and digitized as the stratum lithology data. In other embodiments, other methods can also be used, which are not limited in the present application.

[0059] In some embodiments, the apparent resistivity profile data is subjected to anomaly detection, and the apparent resistivity anomaly area is extracted by the following steps:

[0060] The statistical characteristic value of the background resistivity of the apparent resistivity profile data is determined;

[0061] According to the statistical characteristic value of the background resistivity, the anomaly judgment threshold of relative high resistance and relative low resistance is set;

[0062] Based on the anomaly judgment threshold, the apparent resistivity profile data is subjected to image segmentation, and the independent apparent resistivity anomaly area is extracted and circled.

[0063] It should be noted that the statistical characteristic value of the background resistivity in the present application refers to the statistical quantity used to represent the general electrical property level of the non-polluted medium in the abandoned mine.

[0064] In a specific implementation, the background resistivity statistical characteristic value of the apparent resistivity profile data can be determined by the following method: taking the resistivity values of all data points in the two-dimensional apparent resistivity profile data as the overall sample, first, using an iterative statistical method to eliminate obvious outliers, specifically, calculating the initial average value and standard deviation of all resistivity values, and eliminating data points with resistivity values exceeding the initial average value plus or minus three times the standard deviation as outliers; then, based on the remaining data points, the average value and standard deviation of the resistivity values are recalculated, and the recalculated average value is taken as the background resistivity average value in the background resistivity statistical characteristic value, and the recalculated standard deviation is taken as the background resistivity standard deviation in the background resistivity statistical characteristic value; this process aims to obtain a stable statistic that can represent the general electrical characteristics of the uncontaminated medium in the abandoned mine, and other methods can also be used in other embodiments, which are not limited in the present application.

[0065] In a specific implementation, according to the background resistivity statistical characteristic value, the abnormal judgment threshold of relative high resistance and relative low resistance can be set by the following method: taking the background resistivity average value in the background resistivity statistical characteristic value as the center reference, defining the numerical range higher than the background resistivity average value plus twice the background resistivity standard deviation as the judgment threshold interval of the relative high resistance anomaly; at the same time, defining the numerical range lower than the background resistivity average value minus twice the background resistivity standard deviation as the judgment threshold interval of the relative low resistance anomaly; the multiple of "twice" can be adaptively adjusted according to the actual data variation degree and the survey accuracy requirement, for example, between one times and three times the standard deviation, and other methods can also be used in other embodiments, which are not limited in the present application.

[0066] In a specific implementation, based on the abnormality determination threshold, the apparent resistivity profile data is subjected to image segmentation, and the independent apparent resistivity abnormal area is extracted and circled. The following method can be used: the apparent resistivity profile data is regarded as a two-dimensional gray image, and the pixel value is the resistivity value. First, according to the set relative high-resistance and relative low-resistance abnormality determination threshold, the image is subjected to binaryzation processing, the pixel with the resistivity value falling into the high-resistance threshold interval is marked as a high-resistance potential abnormal pixel, and the pixel with the resistivity value falling into the low-resistance threshold interval is marked as a low-resistance potential abnormal pixel. Then, a connected component analysis algorithm in image processing is used to scan the pixel sets marked as high-resistance and low-resistance respectively, and the same kind of pixels adjacent in space (using four-connected or eight-connected rule) are aggregated together to form a connected pixel block. If the area (i.e. the number of pixels contained) of each connected pixel block is greater than a preset minimum area threshold, for example, the actual area is greater than 10 square meters, the connected pixel block is identified as an independent abnormal area with a clear spatial boundary, i.e. an extracted apparent resistivity abnormal area. Finally, the spatial position coordinates and boundary range of all independent apparent resistivity abnormal areas meeting the area condition are output. In other embodiments, other methods can also be used, which are not limited in the present application.

[0067] It should be noted that the setting of the minimum area threshold in the present application needs to be combined with the accuracy of the geophysical exploration and the minimum resolvable scale of the target geological body. The following method can be used to determine the minimum area threshold: first, according to the electrode spacing and the measurement depth of the high-density resistivity method, the spatial resolution of the apparent resistivity profile data is estimated , which is usually 1-1.5 times the electrode spacing; then, the minimum area threshold is set to the minimum area of a geological abnormal body with a certain continuity and scale that can constitute a pollution migration path on the profile, such as a fracture zone and a pollution plume; a practical method is to set the threshold value to , where is an empirical coefficient, usually taking a value between 4 and 10, to ensure that isolated small anomalies caused by random noise or local heterogeneity are filtered out, while the abnormal area with analysis value is retained; the threshold value can be adjusted as a adjustable parameter in the data processing software according to the electrical background noise level of the actual exploration area and the exploration target; in other embodiments, small-area anomalies can also be removed by automatic methods such as image morphological opening operation, which are not limited in the present application

[0068] In a specific implementation, the extracted apparent resistivity abnormal area is subjected to genetic filtering according to the stratum lithology data, and then the abnormal area irrelevant to the non-pollution geological body is circled as a suspected pollution path area. The following method can be used: the extracted apparent resistivity abnormal area is subjected to genetic filtering according to the stratum lithology data, which can be judged according to the following rules:

[0069] a. If the apparent resistivity anomaly zone is located entirely within the thick, continuously distributed clay or silty clay layer region identified by the stratigraphic lithology data, and its resistivity exhibits a uniform low resistivity characteristic, then the anomaly is determined to be caused by a low-resistivity clay layer and is filtered out. The thickness of the thick layer is greater than a preset value, such as 2 meters, and the low resistivity characteristic is that the resistivity value is lower than the average background resistivity value minus one standard deviation.

[0070] b. If the apparent resistivity anomaly zone is located entirely within the intact bedrock area and exhibits isolated high resistivity characteristics, and there are no known pollution sources or tectonic fracture zones in the surrounding area, then the anomaly is determined to be caused by a high resistivity bedrock mass and is filtered out. The intact bedrock is such as unweathered granite or limestone, and the high resistivity characteristic is that the resistivity value is higher than the average background resistivity plus one standard deviation.

[0071] c. If the spatial location of the apparent resistivity anomaly zone highly overlaps with the known distribution area of ​​pure, saturated sand and gravel aquifers, and its resistivity value is within the typical low to medium resistivity range of that aquifer, then the anomaly is determined to be caused by a clean aquifer and should be filtered out.

[0072] d. For apparent resistivity anomaly areas that do not meet the above rules a, b, and c, especially those located in lithologically complex areas, known goaf areas, waste rock piles, or downstream of tailings ponds, or whose resistivity anomaly morphology cannot be reasonably explained by known clean lithology, they are retained and delineated as suspected contamination pathway areas. The resistivity anomaly morphology is such as strip-shaped or funnel-shaped. In other embodiments, a more detailed electrical-lithological knowledge base can be established or machine learning models can be used for automatic genetic classification. This application does not limit this.

[0073] It should be noted that the above steps automatically and objectively screen out key target areas related to pollution from the wide-area geophysical scan, and significantly reduce the interference of "false anomalies" caused by non-polluted geological bodies such as pure aquifers and clay layers through geological knowledge filtering, thereby improving the pertinence and reliability of the starting point of subsequent exploration work and avoiding ineffective investment.

[0074] In some embodiments, reference Figure 2 As shown in the figure, this figure is a schematic diagram of the application scenario of the pollutant dominant migration path analysis data processing system shown in some embodiments of this application. The figure includes three main components: acquisition equipment, server and data storage equipment. The acquisition equipment is responsible for collecting stratigraphic lithology data of abandoned mines and sending the collected stratigraphic lithology data of abandoned mines to the server through a communication network. The pollutant dominant migration path analysis data processing system runs on the server. The server stores the processing results in the data storage equipment and visualizes them.

[0075] In step 102, the apparent resistivity value and the characteristic pollutant concentration value of the corresponding location and depth range within the spatial range of each suspected pollution path area are extracted, and then the correlation between the spatial distribution trends of the two is analyzed to obtain the physicochemical response synergy degree, which characterizes the spatial coupling degree of physicochemical anomalies in each suspected pollution path area.

[0076] In some embodiments, extracting the apparent resistivity values ​​and corresponding characteristic pollutant concentration values ​​at corresponding locations and depths within the spatial range of each suspected contamination path area can be achieved using the following steps:

[0077] For each suspected contamination pathway area, geochemical sampling and detection are carried out in the suspected contamination pathway area and its adjacent areas to obtain the concentration data of characteristic pollutants in the suspected contamination pathway area;

[0078] The concentration data of the characteristic pollutants are registered with the apparent resistivity profile data in terms of spatial location and depth range to construct a spatially aligned chemical anomaly concentration profile.

[0079] Within the spatial boundary of the suspected contamination pathway area, the registered apparent resistivity values ​​and corresponding characteristic pollutant concentration values ​​are extracted to form a set of data pairs for subsequent correlation analysis.

[0080] It should be noted that the chemical anomaly concentration profile in this application refers to a two-dimensional continuous concentration distribution image in which the pollutant concentration data collected at discrete points are strictly aligned with the apparent resistivity profile in spatial location. Its function is to realize the spatial visualization overlay and point-by-point comparison of the geochemical field and the geophysical field. The apparent resistivity value is a physical quantity that characterizes the conductivity of the medium at a certain point underground. It is the basic data reflecting the change in electrical properties of the underground medium caused by the presence of pollutants or changes in geological structure. The characteristic pollutant concentration value is a chemical quantity that characterizes the level of a specific pollutant in the soil or groundwater.

[0081] In a specific implementation, the geochemical sampling and detection are performed in the suspected pollution path area and its adjacent area, and the concentration data of the characteristic pollutants in the suspected pollution path area can be obtained by the following method: taking the spatial boundary range of each suspected pollution path area defined in the previous step as the core, arranging soil or groundwater sampling points in the form of an encrypted grid along the direction of the corresponding geophysical detection line, and covering the sampling points in the abnormal area and a certain range of background control area outside the abnormal area; using a light drilling equipment to collect undisturbed soil column samples or pore water samples at different depths at each sampling point, wherein the sampling depth interval should at least cover the shallow unsaturated zone, the phreatic surface fluctuation zone and the shallow aquifer, so as to ensure the capture of the vertical distribution characteristics of the pollutants; then, using a portable rapid detection equipment to perform instant analysis on the samples, using a portable X-ray fluorescence spectrometer to measure the concentration of heavy metal elements such as arsenic, lead, cadmium and zinc, and using a portable gas chromatograph-mass spectrometer or a photoionization detector to measure the concentration of volatile organic pollutants, so as to obtain the concentration data of the characteristic pollutants at different depths at each sampling point. In other embodiments, laboratory analysis or other field detection methods can also be used, which are not limited in the present application.

[0082] In some embodiments, the concentration data of the characteristic pollutants is matched with the apparent resistivity profile data in terms of spatial position and depth interval, and the construction of the spatially aligned chemical anomaly concentration profile can be achieved by the following steps:

[0083] The spatial coordinate system and the elevation reference used by the concentration data of the characteristic pollutants and the apparent resistivity profile data are unified;

[0084] Based on the concentration data of the characteristic pollutants in the unified spatial coordinate system, a continuous two-dimensional concentration distribution profile is reconstructed by using a spatial interpolation algorithm;

[0085] The reconstructed concentration distribution profile is aligned with the apparent resistivity profile at the grid nodes to generate the spatially aligned chemical anomaly concentration profile.

[0086] It should be noted that the concentration distribution profile in the present application is a two-dimensional image reflecting the spatial continuous variation trend of the concentration of pollutants, which is preliminarily reconstructed by spatial interpolation and is an intermediate product generated in the process of constructing the chemical anomaly concentration profile.

[0087] In a specific implementation, the spatial coordinate system and the height reference used for unifying the concentration data of the characteristic contaminant and the apparent resistivity profile data can be implemented in the following manner: first, determine the plane coordinate system and the height reference used in geophysical surveying, such as the specific projection zone coordinates of the national geodetic coordinate system CGCS2000 and the 1985 national height reference; subsequently, in the implementation of geochemical sampling, use a global satellite navigation system receiver of the same model or matching accuracy as that used in geophysical surveying to perform real-time differential positioning on each sampling point, directly obtain the coordinates of the sampling point in the same plane coordinate system, and use a level or a global satellite navigation system height measurement function to obtain the height of the sampling point in the same height reference; for all sampling points, unify the position information, including the plane coordinates and the height, to the same coordinate system and height reference as the apparent resistivity profile data, so as to ensure that the concentration data of the characteristic contaminant and the apparent resistivity profile data have a consistent spatial reference framework. In other embodiments, other methods can also be used, which are not limited in the present application.

[0088] In a specific implementation, based on the concentration data of the characteristic contaminant in the unified spatial coordinate system, the continuous two-dimensional concentration distribution profile can be reconstructed by using a spatial interpolation algorithm in the following manner: first, establish a one-to-one correspondence between the concentration value of the characteristic contaminant obtained at different depths of each sampling point and the unified three-dimensional spatial coordinates of the sampling point, that is, the horizontal distance along the survey line and the depth, to form a three-dimensional discrete data set; then, for a target depth interval, select a two-dimensional profile formed by the direction along the survey line and the vertical depth direction as an interpolation region, and estimate the concentration in the interpolation region by using a Kriging spatial interpolation algorithm. Specifically, analyze the variation characteristics of the discrete concentration data in the two-dimensional space, fit an optimal semi-variation function model, such as a spherical model or an exponential model, and determine the corresponding range, nugget value and sill value parameters. Subsequently, based on the model, calculate the optimal unbiased estimate value of the concentration of the characteristic contaminant at each regular grid node on the profile by using ordinary Kriging, and finally generate a two-dimensional raster image composed of continuous grid node values, which reflects the spatial distribution trend of the concentration of the contaminant, that is, reconstruct the continuous two-dimensional concentration distribution profile. In other embodiments, other methods can also be used, which are not limited in the present application.

[0089] In specific implementation, aligning the reconstructed concentration distribution profile with the apparent resistivity profile on the grid nodes to generate the spatially aligned chemical anomaly concentration profile can be achieved as follows: Obtain the grid definition parameters of the apparent resistivity profile data, including the profile start point coordinates, profile azimuth angle, horizontal grid spacing, and vertical grid layering information; using these as a standard, perform grid resampling processing on the aforementioned reconstructed concentration distribution profile. Specifically, for the center point coordinates of each grid node of the apparent resistivity profile, calculate its corresponding position in the reconstructed concentration distribution profile raster image, and use bilinear... The concentration estimate at the corresponding node position is extracted from the concentration distribution profile raster image by interpolation. After traversing and interpolating all apparent resistivity profile grid nodes, a brand new concentration profile is obtained. This brand new concentration profile has the same number of grid nodes, node spatial position, and depth layering as the apparent resistivity profile, that is, the two are spatially aligned node by node. This final concentration profile, which strictly matches the apparent resistivity profile grid, is the spatially aligned chemical anomaly concentration profile. Other methods can be used in other embodiments, and this application does not limit them.

[0090] In specific implementation, extracting the registered apparent resistivity values ​​and corresponding characteristic pollutant concentration values ​​from the spatial boundary of the suspected contamination path area to form a data pair set for subsequent correlation analysis can be achieved in the following way: For each suspected contamination path area, its delineated spatial boundary range is used as a spatial mask and applied to the spatially registered apparent resistivity profile and chemical anomaly concentration profile; at each node with the same horizontal coordinate and the same depth interval within the spatial boundary range, the data values ​​corresponding to the two profiles are read simultaneously: one value is the apparent resistivity value at that location, and the other value is the characteristic pollutant concentration value at that location; such a pair of data with a clear spatial coordinate correlation is recorded as a record, and all valid spatial nodes within the boundary of the anomaly area are traversed to finally form a data pair set of apparent resistivity-concentration specific to the suspected contamination path area; this set is the direct input data for the next step of analyzing the correlation between the spatial distribution trends of the two to calculate the physicochemical response synergy. In other embodiments, batch data export or automatic extraction by programming can also be used, which is not limited in this application.

[0091] In some embodiments, reference Figure 3 As shown in the figure, this is a flowchart illustrating the determination of physicochemical response synergy in some embodiments of this application. In this embodiment, the correlation between the spatial distribution trends of the two is analyzed, and the physicochemical response synergy, which characterizes the degree of spatial coupling of physicochemical anomalies in each suspected contamination path area, can be obtained by the following steps:

[0092] In step 1031, for each suspected pollution path area, a set of apparent resistivity-concentration data pairs of the suspected pollution path area is obtained;

[0093] In step 1032, the data in the set of apparent resistivity-concentration data pairs is rank converted to eliminate dimensional differences and converted into a sequence suitable for non-parametric statistical analysis;

[0094] In step 1033, a non-parametric statistical quantity capable of characterizing the monotonic correlation between the apparent resistivity value and the characteristic pollutant concentration value is determined based on the converted sequence;

[0095] In step 1034, the non-parametric statistical quantity is directly used as the geochemical response synergy degree of the suspected pollution path area.

[0096] It should be noted that the non-parametric statistical quantity in the present application is a mathematical index for measuring the strength and direction of the monotonic correlation between two variable sequences, which is used to robustly quantify the statistical correlation of the spatial trend of geophysical and geochemical data; the geochemical response synergy degree is a comprehensive index for quantifying the consistency of the distribution trend of the apparent resistivity anomaly and the pollutant concentration anomaly in the same spatial range, and its absolute value size represents the coupling strength, and its sign represents the coupling direction (positive / negative correlation), which is a key fusion evidence for judging whether the geophysical anomaly is caused by pollution.

[0097] In specific implementation, the rank conversion of the data in the set of apparent resistivity-concentration data pairs to eliminate dimensional differences and convert into a sequence suitable for non-parametric statistical analysis can be implemented in the following manner, that is, the independent sorting processing is respectively performed on all apparent resistivity values and all characteristic pollutant concentration values in the data pair set; for the apparent resistivity value sequence, all values are arranged in ascending order, and the sorting position of each value, i.e., the rank, is assigned to it, and if the same values are encountered, the position occupied by these same values is averaged as their common rank; the same independent sorting and rank assignment rules are used to process the characteristic pollutant concentration value sequence; after this process, the original apparent resistivity values and concentration values with different physical units and magnitudes are respectively converted into rank sequences representing only their relative size positions in the sequence, thereby eliminating the dimensional differences and non-normal distribution of the original data, and obtaining two columns of rank data with equal length and one-to-one correspondence, i.e., the converted sequence suitable for non-parametric statistical analysis. In other embodiments, other methods can also be used to implement it, which is not limited in the present application.

[0098] In specific implementation, determining the nonparametric statistic that characterizes the monotonic correlation between apparent resistivity values ​​and characteristic pollutant concentration values ​​based on the transformed sequence can be achieved as follows: The Spearman rank correlation coefficient method is used to determine the nonparametric statistic. Specifically, first, the difference between each pair of apparent resistivity ranks and concentration value ranks in the transformed sequence is calculated, and then each difference is squared. Next, the squares of all differences are summed to obtain the sum of squares. Subsequently, according to the Spearman rank correlation coefficient calculation formula, the sum of squares and the total number of data pairs are used to calculate the final result. A value between -1 and 1 is the Spearman rank correlation coefficient. This Spearman rank correlation coefficient is the nonparametric statistic that can characterize monotonic correlation. Its absolute value indicates the strength of the monotonic correlation between apparent resistivity and pollutant concentration in spatial distribution trends. Its positive or negative sign indicates the direction of this correlation, i.e., positive or negative correlation. This nonparametric statistic does not require the original data to meet the assumption of normal distribution or linear relationship. It is suitable for analyzing the spatial trend correlation between geophysical and geochemical data. Other methods can be used in other embodiments, and this application does not limit them.

[0099] It should be noted that, in order to ensure the statistical reliability of the physicochemical response synergy, this application includes a step of performing a statistical significance test on the correlation coefficient after calculating the nonparametric statistics. This can be implemented by: determining the number of data pairs in the apparent resistivity-concentration data set. and the preset significance level ,For example Query the Spearman rank correlation coefficient critical value table to obtain the value at the significance level. The critical value below If the calculated If the correlation coefficient is statistically significant, then it will be determined that the correlation coefficient is statistically significant. The value is directly used as the physicochemical response synergy of the suspected contamination pathway area; if If the correlation coefficient is not statistically significant, the physical-chemical response synergy of the suspected contamination path area is assigned to zero or a preset value that represents no significant correlation. By introducing a statistical significance test, weak or false correlations caused by random data fluctuations can be effectively distinguished from true physical-chemical spatial coupling, thereby improving the credibility of the physical-chemical response synergy as evaluation evidence. In other embodiments, the t-test or other non-parametric test methods can also be used for significance judgment, which is not limited in this application.

[0100] In addition, it should be noted that the above steps achieve accurate fusion and quantitative correlation analysis of the geophysical field and the geochemical field in space, and through calculation of the "geochemical response synergy degree" independent of data distribution, the spatial trend consistency of the two kinds of data is converted into quantifiable evidence for comparison, thereby improving the objectivity and scientificity of the pollution cause identification and overcoming the defect of strong multi-solution of single data.

[0101] In step 103, the dominant channel index of each suspected pollution path area is determined according to the geochemical response synergy degree of each suspected pollution path area and the extension of the suspected pollution path area along the potential groundwater flow direction on the profile.

[0102] In some embodiments, the determination of the dominant channel index of each suspected pollution path area according to the geochemical response synergy degree of each suspected pollution path area and the extension of the suspected pollution path area along the potential groundwater flow direction on the profile can be achieved by the following steps, that is:

[0103] determining a unified potential groundwater flow direction in the abandoned mine based on terrain and geological structure data;

[0104] For each suspected pollution path area, determining the projection length of the suspected pollution path area along the potential groundwater flow direction on the profile as a quantitative indicator of the extension;

[0105] respectively normalizing the geochemical response synergy degree and the projection length to obtain corresponding standardized evaluation values;

[0106] performing weighted summation on the standardized evaluation value of the geochemical response synergy degree and the standardized evaluation value of the projection length to obtain the dominant channel index of the suspected pollution path area.

[0107] It should be noted that the dominant channel index in the present application is a composite evaluation index that combines the geochemical response synergy degree and the extension along the groundwater flow direction, and its role is to prioritize all suspected areas to identify the channel that is most likely to become the main migration path of the pollutants.

[0108] In some embodiments, the determination of a unified potential groundwater flow direction in the abandoned mine based on terrain and geological structure data can be achieved by the following steps, that is:

[0109] obtaining digital elevation model data and geological structure maps of the abandoned mine;

[0110] analyzing the digital elevation model data to extract the convergence trend of the surface water system and the terrain slope direction of the abandoned mine;

[0111] According to the geological structure maps, the occurrence of the rock stratum of the main aquifer in the abandoned mine and the trend of the dominant structure line are identified;

[0112] According to the flow convergence tendency, the terrain slope direction, the rock stratum occurrence and the dominant structural line direction, a unified potential underground water flow direction in the abandoned mine is determined.

[0113] It should be noted that the flow convergence tendency in the present application is tendency information representing the overall convergence and flow direction of the surface water body, and its role is to provide macroscopic surface hydrology basis for inferring the regional potential underground water flow field; the terrain slope direction refers to the main direction of the overall inclination of a region, and its role is to serve as an important topographic factor indicating the potential direction of underground water flow, which is usually related to the direction of underground water potential energy reduction; the rock stratum occurrence refers to the three-dimensional spatial orientation information of the underground rock stratum in space, and its role is to reveal the advantage direction of the rock stratum itself that may control or affect the underground water flow from the geological structure level; the dominant structural line direction refers to the extension direction of the largest and longest linear geological structure such as fault and fracture zone in the investigation area, and its role is to identify the direction of the key geological structure that may become a high-speed migration channel or control the boundary of water flow.

[0114] In a specific implementation, the digital elevation model data and the geological structure map of the abandoned mine can be obtained in the following manner: the digital elevation model data of the abandoned mine can be obtained from the national or provincial basic geographic information center or the relevant surveying and mapping department, and the digital elevation model data is usually stored in a regular grid form, and the precision can be selected according to the investigation range, such as 5-meter or 10-meter grid resolution; at the same time, the regional geological map, hydrogeological map and existing geological exploration report prepared by the geological and mineral resources department are collected, and the geological structure elements are extracted and digitized to form a digital geological structure map containing information such as stratum boundary, fault line and rock stratum occurrence symbol. In other embodiments, other methods can also be used to achieve this, and the present application does not limit this.

[0115] In a specific implementation, the analysis of the digital elevation model data to extract the runoff trend of the surface water system of the abandoned mine and the terrain slope direction can be achieved by using the following method: using the hydrological analysis tool of the geographic information system software to process the digital elevation model data, first, performing depression filling to eliminate data depressions; second, using a standard water flow direction algorithm, such as the D8 single flow direction algorithm, to calculate the water flow direction of each grid cell to obtain water flow direction data, the principle of the D8 single flow direction algorithm is: for the center grid, comparing the elevation difference with the eight surrounding grids, and assigning the water flow direction to the direction of the grid with the largest elevation difference; then, based on the obtained water flow direction data, using the flow accumulation algorithm to calculate, that is, starting from all grid cells, following the water flow direction, accumulating the number of all upstream grids (including itself) received by each grid cell, and the accumulation value is the flow accumulation of each grid cell, which simulates the upstream catchment area passing through the point; then, by setting a flow accumulation threshold, for example, greater than 500 grid cells, connecting the grid cells with a flow accumulation greater than the flow accumulation threshold, which can identify and extract the main surface water system network; finally, analyzing the overall structure of the surface water system network, especially the extension direction of the main stream, as the runoff trend of the surface water system, at the same time, using the digital elevation model data to calculate the surface slope and slope direction of the entire region, and determining the terrain slope direction by statistical analysis (such as slope rose diagram), that is, the main direction of the overall terrain inclination, in other embodiments, other methods can also be used to achieve this, which is not limited in this application.

[0116] In a specific implementation, according to the geological structure map, the occurrence of the main aquifer rock layer of the abandoned mine and the trend of the dominant structure line can be achieved by using the following method: on the digitized geological structure map, positioning and identifying the distribution range of the main water-bearing rock layer in the abandoned mine, such as sand and gravel layer, fractured bedrock, etc.; reading the rock layer occurrence symbol marked on the geological structure map of the main water-bearing rock layer, usually represented by trend, inclination, and dip angle, to determine the occurrence of the rock layer; at the same time, identifying the largest scale and longest extension fault or large fracture zone on the geological structure map, defining it as the dominant structure line, and measuring its extension direction on the geological map, that is, the trend of the dominant structure line, in other embodiments, other methods can also be used to achieve this, which is not limited in this application.

[0117] In a specific implementation, the unified potential underground water flow direction in the abandoned mine can be determined according to the convergence trend, the topographic slope direction, the rock stratum occurrence and the dominant structure line direction in the following manner: the convergence trend of the surface water system and the topographic slope direction are used as important surface references for indicating the general flow direction of the shallow groundwater; meanwhile, the rock stratum occurrence (especially the rock stratum tendency) of the main aquifer and the dominant structure line direction are used as key underground geological factors for controlling the deep groundwater flow; the four direction information are comprehensively compared, if they show high consistency, for example, the topographic slope direction, the water system convergence direction, the rock stratum tendency and the structure line direction all roughly point to the same direction, the direction is determined as the unified potential underground water flow direction in the abandoned mine; if there is a difference, that is, when the hydrogeological structure of the abandoned mine is complex, there are multiple independent hydrogeological units or the underground water flow field direction changes significantly, the rock stratum occurrence and the dominant structure line direction reflecting the control of the underground structure are given greater weight according to the principle of hydrogeology, and the topographic information is combined for judgment, and finally a representative water flow direction for the whole area evaluation is determined as the unified potential underground water flow direction in the form of angle value or direction description. In other embodiments, other methods can also be used to achieve this, which are not limited in the present application.

[0118] It should be noted that in some implementations, when the hydrogeological structure of the abandoned mine is complex, there are multiple independent hydrogeological units or the underground water flow field direction changes significantly, the potential underground water flow direction can be determined in a regional or profile manner, which can be achieved in the following steps: according to the digital elevation model, the geological structure zoning and the drilling hydrogeological data, the entire abandoned mine research area is divided into several sub-areas with relatively uniform hydrogeological characteristics; for each sub-area, independently perform: analyze the topographic slope direction, the surface water system convergence trend, the rock stratum occurrence of the main aquifer and the structure line direction in the sub-area, and comprehensively determine the potential underground water flow direction in the sub-area; accordingly, when calculating the projection length of each suspected pollution path area, the potential underground water flow direction determined by the sub-area where the suspected pollution path area is located is used for calculation; if the suspected pollution path area crosses multiple sub-areas, the direction of the sub-area where the center point belongs or the weighted value of the projection length in different sub-area directions can be taken; this way enhances the adaptability of the method to complex hydrogeological conditions, making the evaluation of the dominant channel index more consistent with the actual underground water dynamic field characteristics; in other embodiments, a spatially continuous water flow direction field can also be obtained through numerical simulation of the underground water flow, which is not limited in the present application.

[0119] In a specific implementation, the projection length of the suspected pollution path area along the direction of the underground water flow on the profile can be quantified by the following method: on a two-dimensional profile containing all suspected pollution path areas, the two-dimensional profile contains horizontal distance and depth coordinate axes, and an arrow indicating line representing the unified direction of the underground water flow is drawn; for each suspected pollution path area, the boundary of the abnormal area circled on the two-dimensional profile is projected along the direction perpendicular to the arrow line of the water flow direction onto a straight line parallel to the arrow line of the water flow direction, and the length of the line segment covered by the projection of the abnormal area on the straight line is measured, which is defined as the projection length of the corresponding suspected pollution path area along the direction of the underground water flow on the profile. The projection length quantifies the spatial distribution scale of the corresponding abnormal area along the possible hydraulic migration direction. In other embodiments, other methods can also be used, which are not limited in the present application.

[0120] In a specific implementation, the normalization of the materialization response synergy degree and the projection length can be realized by the following method: first, collect the materialization response synergy degree values of all suspected pollution path areas, and find the maximum and minimum values; then, for any suspected pollution path area, subtract the minimum value of the overall synergy degree from the materialization response synergy degree value, and divide by the difference between the maximum and minimum values of the overall synergy degree, to calculate the normalized evaluation value of the materialization response synergy degree of the suspected pollution path area, which will fall within the interval of 0 to 1; in the same way, based on the projection length data of all suspected pollution path areas, the maximum and minimum values are found, and the normalized evaluation value of the projection length of each suspected pollution path area is calculated. In other embodiments, other methods can also be used, which are not limited in the present application.

[0121] In a specific implementation, the weighted sum of the normalized evaluation value of the materialization response synergy degree and the normalized evaluation value of the projection length can be realized by the following method: according to the formation mechanism of the pollution dominant migration channel, the normalized evaluation value of the materialization response synergy degree is given a weight , and the normalized evaluation value of the projection length is given a weight , where ; the specific weight can be determined by one of the following methods: method one, empirical value, based on domain knowledge, since the materialization response synergy degree directly represents the evidence strength of pollution, it is usually given a higher weight, for example, set between 0.6 and 0.7, and the normalized evaluation value of the projection length between 0.3 and 0.4, the second mode objective weighting is based on historical case data or an evaluation matrix of multiple suspected pollution path areas in the same site, the weight is determined by using an entropy weight method, specifically including, constructing an evaluation matrix composed of all suspected pollution path areas and composing an evaluation matrix composed of all suspected pollution path areas and ; for the suspected pollution path area, the advantage channel index is calculated according to the formula , in other embodiments, the weight can also be dynamically adjusted according to a specific remediation engineering target, such as priority blocking or priority monitoring, by using an analytic hierarchy process, a coefficient of variation method or the like, which is not limited in the present application.

[0122] It should be noted that the above steps comprehensively evaluate the potential channel from two dimensions of "abnormal coupling strength" and "hydrogeological function", combine the physical and chemical response synergy degree and the extension along the groundwater flow direction into a unified "advantage channel index", thereby improving the logical rationality of prioritizing multiple suspected areas and the decision support intensity, and making the channel identification more in line with the physical mechanism of pollutant migration.

[0123] In step 104, the suspected pollution path area with the highest advantage channel index is determined as the dominant migration channel of the pollutant, and is labeled in combination with the geological profile to provide a targeted position for setting a groundwater monitoring well or implementing an in-situ barrier remediation engineering.

[0124] In some embodiments, the suspected pollution path area with the highest advantage channel index is determined as the dominant migration channel of the pollutant, and is labeled in combination with the geological profile to provide a targeted position for setting a groundwater monitoring well or implementing an in-situ barrier remediation engineering, which can be realized by the following steps, namely:

[0125] sorting and screening the advantage channel indexes of all suspected pollution path areas to determine the suspected pollution path area with the highest advantage channel index;

[0126] determining the suspected pollution path area with the highest advantage channel index as the dominant migration channel of the pollutant, and extracting spatial feature information of the dominant migration channel of the pollutant;

[0127] superimposing and labeling the spatial feature information on the geological profile map of the abandoned mine to form a comprehensive result map for guiding engineering implementation, so as to provide a targeted position for setting a groundwater monitoring well or implementing an in-situ barrier remediation engineering.

[0128] In specific implementation, the dominant channel index of all suspected contamination path areas is sorted and screened to determine the suspected contamination path area with the highest dominant channel index. This can be achieved by sorting the dominant channel indices of all suspected contamination path areas in descending order of value to determine the suspected contamination path area with the highest dominant channel index. Other methods can also be used in other embodiments, and this application does not limit them.

[0129] In specific implementation, the selected suspected pollution path areas are identified as dominant pollutant migration channels, and the spatial feature information of the dominant pollutant migration channels is extracted. This can be achieved in the following way: the selected suspected pollution path areas are formally assigned the attribute label of dominant pollutant migration channels; then, for the suspected pollution path areas corresponding to the dominant pollutant migration channels, the planar distribution range and spatial location are extracted from the abnormal area boundary data delineated in step 101, the corresponding potential groundwater flow direction is extracted from the water flow direction data determined in step 104, and the specific value is extracted from the dominant channel index data calculated in step 104. These data, namely spatial location, boundary range, dominant channel index, and potential groundwater flow direction, are integrated into the spatial feature information of the dominant pollutant migration channel. In other embodiments, additional information such as the center point coordinates, area, and major axis direction can also be extracted, but this application does not limit this.

[0130] In practice, the spatial feature information is overlaid and annotated onto the geological profile of the abandoned mine to form a comprehensive result map guiding the project implementation. This map provides target locations for setting up groundwater monitoring wells or implementing in-situ barrier remediation projects. This can be achieved in the following way: First, obtain a base map of the geological profile of the abandoned mine, which includes stratigraphic boundaries, lithological symbols, and necessary scale and legend. Then, using a geographic information system or professional mapping software, the spatial feature information of the extracted dominant pollutant migration channels is overlaid onto the base map as graphic elements. Specifically, this includes: clearly outlining the boundary of each channel using closed polygons, marking its dominant channel index within or beside the polygons, and... Arrow symbols are used to indicate the potential direction of groundwater flow. Simultaneously, based on a comprehensive consideration of the boundary range, flow direction, and the location of downstream sensitive targets, prominent markings, such as asterisks or drilling rig symbols, are used to clearly indicate the preferred target locations for deploying groundwater monitoring wells or implementing in-situ vertical barrier walls within the channel boundary or at appropriate locations adjacent to the downstream. Finally, the map is refined by adding map titles, descriptions, etc., to generate a comprehensive result map that intuitively displays the spatial distribution of dominant pollutant migration channels and engineering layout recommendations. This comprehensive result map can be directly used to guide the design and construction of subsequent environmental remediation projects. In other embodiments, it can also be presented in the form of three-dimensional visualization or thematic maps; this application does not limit this.

[0131] It should be noted that the above steps directly convert the quantitative analysis results into intuitive and operable engineering guidance drawings. By accurately marking the identified advantage migration channel and the recommended engineering target position on the geological profile, the targeting, one-time success rate and treatment efficiency of subsequent groundwater monitoring well layout or in-situ barrier repair engineering are improved, and the close connection from technical detection to engineering implementation is realized.

[0132] In addition, another aspect of the present application, in some embodiments, the present application provides a waste mine land resource recycling soil pollutant analysis system, referring to Figure 4 The figure is a structural schematic diagram of a waste mine land resource recycling soil pollutant analysis system according to some embodiments of the present application, which includes an identification module 401, a processing module 402 and an execution module 403, which are described as follows:

[0133] The identification module 401 is mainly used for identifying and delineating a plurality of apparent resistivity anomaly zones caused by pollutants as suspected pollution path zones according to the stratum lithology data of the waste mine in the present application;

[0134] The processing module 402 is mainly used for extracting the apparent resistivity value within the spatial range of each suspected pollution path zone and the characteristic pollutant concentration value of the corresponding position and corresponding depth interval, and then analyzing the correlation of the spatial distribution trend of the two to obtain the degree of coupling of the physical and chemical response of each suspected pollution path zone;

[0135] The processing module 402 in the present application is also used for determining the advantage channel index of each suspected pollution path zone according to the degree of coupling of the physical and chemical response of each suspected pollution path zone and its extension along the direction of potential groundwater flow on the profile;

[0136] The execution module 403 is mainly used for determining the suspected pollution path zone with the highest advantage channel index as the pollutant advantage migration channel, and marking it in combination with the geological profile to provide a target position for setting a groundwater monitoring well or implementing an in-situ barrier repair engineering.

[0137] Each module in the above waste mine land resource recycling soil pollutant analysis system can be realized by software, hardware and their combinations. The above modules can be embedded in or independent of the processor in the computer device in hardware form, or can be stored in the memory in the computer device in software form, so as to be called and executed by the processor to perform the operations corresponding to the above modules.

[0138] In addition, in one embodiment, the present application provides a computer device, which can be a server, and its internal structure diagram can be as follows:Figure 5 As shown in the figure. The computer device includes a processor, a memory and a network interface connected by a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium, an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store the abandoned mine land resource recycling soil pollution analysis data. The network interface of the computer device is used to communicate with the external terminal through the network connection. The computer program is executed by the processor to realize an abandoned mine land resource recycling soil pollution analysis method.

[0139] Those skilled in the art can understand that, Figure 5 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0140] In one embodiment, a computer device is also provided, including a memory and a processor, the memory storing a computer program, and the processor executing the computer program to realize the steps in the above-mentioned abandoned mine land resource recycling soil pollution analysis method embodiments.

[0141] In one embodiment, a computer readable storage medium is provided, storing a computer program, which is executed by a processor to realize the steps in the above-mentioned abandoned mine land resource recycling soil pollution analysis method embodiments.

[0142] In one embodiment, a computer program product or computer program is provided, which includes computer instructions stored in a computer readable storage medium. The processor of the computer device reads the computer instructions from the computer readable storage medium, and the processor executes the computer instructions, so that the computer device executes the steps in the above-mentioned abandoned mine land resource recycling soil pollution analysis method embodiments.

[0143] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when executed, can include the processes of the above-mentioned embodiment methods. Any reference to memory, storage, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0144] Any combination of the technical features of the above embodiments can be made. In order to make the description simple, all possible combinations of the technical features in the above embodiments are not described, but as long as the combination of the technical features does not exist, it should be considered as the scope of the present application.

[0145] The above embodiments only express several implementation manners of the present application, and the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent. It should be pointed out that for those skilled in the art, without departing from the concept of the present application, some modifications and improvements can be made, which are all within the protection scope of the present application. Therefore, the protection scope of the patent of the present application should be subject to the appended claims.

Claims

1. A method for analyzing soil pollutants in a land resource recycling of an abandoned mine, for analyzing a dominant migration path of pollutants in an abandoned mine, characterized by, The method comprises the following steps: According to the stratum lithology data of the abandoned mine, a plurality of apparent resistivity anomaly zones caused by pollutants are identified and circled as suspected pollution path zones; Extract the apparent resistivity values within the spatial range of each suspected pollution path zone and the characteristic pollutant concentration values corresponding to the position and depth interval, and then analyze the correlation of the spatial distribution trends of the two to obtain the physical and chemical response synergy degree of each suspected pollution path zone representing the spatial coupling degree of the physical and chemical anomalies; According to the physical and chemical response synergy degree of each suspected pollution path zone and its extension along the potential groundwater flow direction on the profile, determine the dominant channel index of each suspected pollution path zone; Determine the suspected pollution path zone with the highest dominant channel index as the dominant migration channel of the pollutants, and mark it in combination with the geological profile to provide a targeted position for setting up a groundwater monitoring well or implementing an in-situ barrier repair project; The method comprises the following steps: Determine the unified potential groundwater flow direction in the abandoned mine based on the terrain and geological structure data; For each suspected pollution path zone, determine the projection length of the suspected pollution path zone along the potential groundwater flow direction on the profile as a quantitative indicator of its extension; Respectively normalize the physical and chemical response synergy degree and the projection length to obtain the corresponding standardized evaluation values; Weighted sum the standardized evaluation values of the physical and chemical response synergy degree and the projection length to obtain the dominant channel index of the suspected pollution path zone.

2. The method of claim 1, wherein, According to the stratum lithology data of the abandoned mine, a plurality of apparent resistivity anomaly zones caused by pollutants are identified and circled as suspected pollution path zones, which specifically comprises: Obtain the apparent resistivity profile data and stratum lithology data of the abandoned mine; Perform anomaly detection on the apparent resistivity profile data to extract apparent resistivity anomaly zones; According to the stratum lithology data, filter the extracted apparent resistivity anomaly zones to further circulate the anomaly zones unrelated to non-pollution geological bodies as suspected pollution path zones.

3. The method of claim 2, wherein, The method comprises the following steps: Determine the background resistivity statistical characteristic value of the apparent resistivity profile data; According to the background resistivity statistical characteristic value, set the anomaly judgment threshold of the relative high resistance and the relative low resistance; Based on the anomaly judgment threshold, perform image segmentation on the apparent resistivity profile data to extract and circulate independent apparent resistivity anomaly zones.

4. The method of claim 1, wherein, The method comprises the following steps: For each suspected pollution path zone, perform geochemical sampling and detection in the suspected pollution path zone and its adjacent area to obtain the concentration data of the characteristic pollutants in the suspected pollution path zone; Align the concentration data of the characteristic pollutants with the apparent resistivity profile data in space position and depth interval to construct a spatially aligned chemical anomaly concentration profile; From the spatial boundary range of the suspected pollution path area, the registered apparent resistivity value and the corresponding characteristic pollutant concentration value are extracted to form a data pair set for subsequent analysis of correlation.

5. The method of claim 4, wherein, The concentration data of the characteristic pollutant and the apparent resistivity profile data are registered in terms of spatial position and depth interval to construct a spatially aligned chemical anomaly concentration profile, which specifically includes: Unifying the spatial coordinate system and the elevation reference system used by the concentration data of the characteristic pollutant and the apparent resistivity profile data; Based on the concentration data of the characteristic pollutant in the unified spatial coordinate system, a continuous two-dimensional concentration distribution profile is reconstructed using a spatial interpolation algorithm; Aligning the reconstructed concentration distribution profile and the apparent resistivity profile at the grid nodes to generate the spatially aligned chemical anomaly concentration profile.

6. The method of claim 1, wherein, Analyzing the correlation of the spatial distribution trends to obtain the degree of coupling of the physical and chemical anomalies in each suspected pollution path area, which specifically includes: For each suspected pollution path area, a set of apparent resistivity-concentration data pairs of the suspected pollution path area is obtained; The data in the set of apparent resistivity-concentration data pairs is rank-transformed to eliminate dimensional differences and converted into a sequence suitable for non-parametric statistical analysis; Based on the converted sequence, a non-parametric statistical quantity is determined, which can represent the monotonic correlation between the apparent resistivity value and the concentration value of the characteristic pollutant; The non-parametric statistical quantity is directly used as the physical and chemical response synergy degree of the suspected pollution path area.

7. A system for analyzing soil pollutants in land resources for reclamation of abandoned mines, which analyzes soil pollutants in land resources for reclamation of abandoned mines using the method according to any one of claims 1 to 6, characterized by, The system includes: An identification module for identifying and delineating a plurality of apparent resistivity anomaly areas caused by pollutants as suspected pollution path areas according to the stratum lithology data of the abandoned mine; A processing module for extracting apparent resistivity values within the spatial range of each suspected pollution path area and corresponding position and depth interval characteristic pollutant concentration values, and then analyzing the correlation of the spatial distribution trends to obtain the degree of coupling of the physical and chemical anomalies in each suspected pollution path area, which specifically includes: The processing module is also used to determine the dominant channel index of each suspected pollution path area according to the physical and chemical response synergy degree of each suspected pollution path area and its extension along the direction of the potential groundwater flow on the profile; An execution module for determining the suspected pollution path area with the highest dominant channel index as the dominant pollutant migration channel, and labeling it in combination with the geological profile to provide a targeted position for setting up a groundwater monitoring well or implementing an in-situ remediation project.

Images