Polluted Soil Mapping Simulation Method, System, Device and Storage Medium

Abstract

The present invention discloses a contaminated soil mapping simulation method, system, device, and storage medium. The method comprises establishing an excess standard analysis table and a layer number-thickness comparison table based on soil pollution analysis results; importing the excess standard analysis table and the layer number-thickness comparison table; filtering the excess standard analysis data for the corresponding stratum according to the layer number; eliminating sampling points for each pollution category that were not submitted for inspection; interpolating the excess standard analysis data for each pollution category to obtain patches with different excess standard attributes for the corresponding pollution category; cropping the patches with different excess standard attributes for each pollution category and fusing patches with the same excess standard attributes; eliminating patches that do not exceed the standard to obtain a pollution distribution map for the corresponding pollution category; and superimposing the pollution distribution maps for all pollution categories in the stratum to form a pollution distribution map for the corresponding stratum. The present invention greatly reduces manual intervention and significantly shortens mapping time.

Description

Technical Field

[0001] The present invention belongs to the technical field of environmental protection data processing, and in particular relates to a GIS-based contaminated soil mapping simulation method, system, equipment and storage medium. Background Art

[0002] According to relevant national policies, before construction land is transferred, developed or reused, it is necessary to conduct a soil pollution survey, risk assessment and remediation / control plan formulation. The purpose includes obtaining existing and potential pollution factors and their impact range, and formulating remediation processes based on the distribution of different pollution factors.

[0003] After the soil pollution survey results are completed, if the risk assessment and remediation plan stage is entered, the pollution distribution needs to be simulated and mapped. At the same time, considering the subsequent project implementation costs, multiple processes are required to conduct multi-faceted comparisons from the perspectives of plane layout, implementation cycle, project volume and cost. Often, the same batch of survey data needs to be simulated and mapped multiple times. With the development of smart construction sites, considering the needs of the subsequent construction stage, the results need to be able to be quickly converted into BIM and CAD to facilitate project implementation management and quickly enter the design stage.

[0004] Currently, pollution distribution mapping simulation typically uses a manual single-factor mapping synthesis method. The specific operations include: interpolation simulation of a single pollution factor (or pollution index), clipping of the plot redline boundary, and screening of the boundary line exceeding the standard. If there are n layers from the ground to the depth of the exceedance, a single pollution factor needs to be processed n times. If each pollution category contains m pollution factors, n×m processing times are required. This makes the mapping simulation operation cumbersome and inefficient, often requiring a large amount of manpower and time, and unable to quickly implement batch processing, multiple scheme comparison, and conversion between different software. At the same time, for each pollution category, the mapping of all pollution factors of that pollution category needs to be superimposed according to the stratigraphic layer to obtain the pollution distribution map of the corresponding pollution category. The pollution distribution maps of all pollution categories are then merged to obtain a comprehensive pollution distribution map. Based on the comprehensive pollution distribution map, the contaminated area of each pollution category is manually calculated using the fill function of software (such as ArcGIS software, CAD software, or SUFER software), which is inefficient (especially CAD software or SUFER software). Summary of the Invention

[0005] The purpose of the present invention is to provide a contaminated soil mapping simulation method, system, equipment and storage medium to solve the problem that the traditional manual single-factor mapping method leads to cumbersome mapping simulation operations and low efficiency.

[0006] The present invention solves the above technical problems through the following technical solutions: a contaminated soil mapping simulation method, the method comprising the following steps:

[0007] Step 1: Obtain soil pollution analysis results, and create a raw data table, an indicator threshold table, an exceeding standard analysis table, and a layer number-layer thickness comparison table based on the soil pollution analysis results; wherein the exceeding standard analysis table includes the sampling point number, sampling point coordinates, the layer number to which the sampling point belongs, the exceeding standard multiple of each pollution indicator, and the exceeding standard information of each pollution category;

[0008] Step 2: Import the above-mentioned excess analysis table and layer number-layer thickness comparison table into ArcGIS software;

[0009] Step 3: Based on the exceeding standard analysis table, filter out the exceeding standard analysis data of the corresponding stratum according to the layer number;

[0010] Step 4: Based on the excessive analysis data of the strata screened out in step 3, the sampling points that have not been inspected in each pollution category are eliminated;

[0011] Step 5: performing interpolation calculation on the exceeding standard analysis data of each pollution category processed in step 4 to obtain different exceeding standard attribute patches corresponding to the pollution category;

[0012] Step 6: Crop the patches with different exceeding attributes for each pollution category and fuse the patches with the same exceeding attributes;

[0013] Step 7: Eliminate the patches that do not exceed the standard from the patches fused in step 6 to obtain a pollution distribution map corresponding to the pollution category;

[0014] Step 8: superimposing the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum;

[0015] Step 9: Repeat steps 3 to 8 to obtain the pollution distribution map of each stratum;

[0016] Step 10: For the pollution distribution map of each stratum, add the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, add the pollution combination type, calculate the pollution volume of each pollution category, and add the pollution patch number and inflection point.

[0017] Furthermore, in step 1, obtaining soil pollution analysis results specifically includes:

[0018] Obtain basic data during the implementation of the basic sampling plan, wherein the basic data includes the initial sampling point number, the inspection depth, the actual sampling location, the test results of each pollution index in the sample, and the maximum sampling depth of each initial sampling point;

[0019] Conduct an index exceeding standard analysis based on the test results of each pollution index, and mark the untested pollution index, undetected pollution index and pollution index exceeding the standard according to the index exceeding standard analysis results;

[0020] Perform an indicator bottom-up analysis based on the indicator exceeding the standard analysis results. If the bottom-up analysis is not successful, determine the deepening plan for the initial sampling point based on the indicator bottom-up analysis results; repeat the steps of obtaining basic data, indicator exceeding the standard analysis, and indicator bottom-up analysis during the implementation of the deepening plan until the bottom-up analysis is successful, and obtain the analysis results of all indicator exceeding the standard;

[0021] Conduct perimeter enclosure analysis on the encrypted sampling points in the basic layout plan to determine the encrypted sampling point layout plan, wherein the encrypted sampling point layout plan includes the drilling depth, inspection depth, pollution index and index detection method of each encrypted sampling point; repeat the steps of obtaining basic data, index exceeding standard analysis and index guarantee analysis during the implementation of the encrypted sampling point layout plan.

[0022] Furthermore, the basic point distribution plan is formulated based on the pollution type of the land to be investigated or analyzed and relevant policies. The basic point distribution plan includes an initial sampling point distribution plan and an encrypted sampling point distribution plan. The initial sampling point distribution plan includes the distribution location and density of the initial sampling points, the planned drilling depth, the point number, the inspection depth, the pollution index and the index detection method; the encrypted sampling point distribution plan includes the distribution location, density and point number of the encrypted sampling points.

[0023] Furthermore, the specific implementation process of the indicator exceeding standard analysis is as follows:

[0024] The detection result of the pollution index is compared with the corresponding index threshold. If it exceeds the index threshold, the pollution index exceeds the standard; if it does not exceed the index threshold, the pollution index does not exceed the standard.

[0025] Furthermore, the specific implementation process of the indicator bottom-up analysis is as follows:

[0026] Determine the maximum exceeding depth and the maximum inspection depth according to the exceeding analysis results of the indicators;

[0027] If the maximum exceeding depth is equal to the maximum inspection depth, it indicates that the bottom is not covered. If the bottom is not covered and no samples are retained, the initial sampling point corresponding to the non-coverage factor will be drilled deeper. If the bottom is not covered and samples are retained, there is no need to drill deeper and the retained samples will be sent for inspection.

[0028] If the maximum exceeding depth is less than the maximum inspection depth, it indicates a safety net.

[0029] Furthermore, the specific implementation process of the perimeter enclosure analysis is as follows:

[0030] Taking any encrypted sampling point as the target point, taking the initial sampling points around the target point as the pending points, entering the point code of the pending point, obtaining and updating the excess analysis results of the pending point according to the point code, and then obtaining and updating the excess analysis results of the pending point when different encrypted sampling points are used as the target points;

[0031] Based on the exceeding standard analysis results of all pending points at each target point, the drilling depth, inspection depth, pollution index and index detection method of the target point are determined.

[0032] Furthermore, the inspection depth of the target point is equal to the union of the inspection depths of all the points to be determined that exceed the standard, the pollution index of the target point is equal to the union of the exceeding standard indexes of all the points to be determined, and the drilling depth of the target point is greater than the maximum inspection depth of the target point.

[0033] Furthermore, in step 5, Thiessen polygons are used to perform interpolation calculation on the exceeding standard analysis data of each pollution category.

[0034] Furthermore, in step 6, when the exceeding attribute includes exceeding and not exceeding, the patches corresponding to all exceeding sampling points are fused, and the patches corresponding to all not exceeding sampling points are fused to obtain exceeding patches and not exceeding patches for each pollution category;

[0035] When the exceeding standard attribute includes different exceeding standard multiples, the patches corresponding to the sampling points with the same exceeding standard multiples are merged to obtain the exceeding standard multiples patches for each pollution category; the exceeding standard multiple corresponding to no exceeding standard is 0.

[0036] Furthermore, before step 2, the method further includes using ArcGIS software to create a folder for storing process files under the operation directory.

[0037] Based on the same concept, the present invention also provides a contaminated soil mapping simulation system, which is obtained by secondary development of ArcGIS software and includes:

[0038] An import module is used to import an exceeding standard analysis table and a layer number-layer thickness comparison table. The exceeding standard analysis table and the layer number-layer thickness comparison table are established based on the soil pollution analysis results. The exceeding standard analysis table includes the sampling point number, the sampling point coordinates, the layer number to which the sampling point belongs, the exceeding standard multiple of each pollution indicator, and the exceeding standard information of each pollution category;

[0039] A screening module, configured to screen out the exceeding standard analysis data of the corresponding stratum according to the layer number based on the exceeding standard analysis table;

[0040] A first elimination module is used to eliminate the sampling points that have not been inspected in each pollution category based on the excessive analysis data of the stratum screened by the screening module;

[0041] An interpolation calculation module, configured to perform interpolation calculation on the excessive-standard analysis data of each pollution category processed by the first elimination module to obtain different excessive-standard attribute patches corresponding to the pollution category;

[0042] The cropping and fusion module is used to crop patches with different excessive attributes for each pollution category and to fuse patches with the same excessive attributes;

[0043] A second elimination module is used to eliminate the patches that do not exceed the standard from the patches processed by the clipping and fusion module to obtain a pollution distribution map corresponding to the pollution category;

[0044] an overlay module, configured to overlay the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum;

[0045] The field adding and calculating module is used for the pollution distribution map of each stratum, adding the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, adding the pollution combination type, calculating the pollution volume of each pollution category, and adding the pollution patch number and inflection point.

[0046] Based on the same concept, the present invention further provides an electronic device, comprising:

[0047] Memory for storing computer programs;

[0048] A processor is configured to implement the contaminated soil mapping simulation method as described above when executing the computer program.

[0049] Based on the same concept, the present invention also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the contaminated soil mapping simulation method described above is implemented.

[0050] Beneficial effects

[0051] Compared with the prior art, the advantages of the present invention are:

[0052] The present invention conducts secondary development of ArcGIS software. It only needs to import the exceedance analysis table and the layer number-layer thickness table, select the interpolation method and the clipping range, and then automatically screen the exceedance analysis data of each stratum, eliminate the sampling points of each pollution category in the stratum that have not been inspected, perform interpolation calculations, clip and merge, and overlay to obtain the pollution distribution map of each stratum. This greatly reduces manual participation and significantly shortens the mapping time. It can quickly complete the comparison of the engineering quantities and costs of different plans, and facilitates the determination of the optimal remediation plan. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] In order to more clearly illustrate the technical solution of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only one embodiment of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0054] Figure 1 This is a flow chart of a contaminated soil mapping simulation method according to an embodiment of the present invention;

[0055] Figure 2 This is a schematic diagram of the excess analysis in the embodiment of the present invention;

[0056] Figure 3 is a pollution index diagram included in different pollution categories in an embodiment of the present invention;

[0057] Figure 4 This is a layer number-layer thickness comparison table in an embodiment of the present invention;

[0058] Figure 5 is a flow chart for obtaining soil pollution analysis results in an embodiment of the present invention;

[0059] Figure 6 This is the interface for importing the over-standard analysis table and the layer number-layer thickness comparison table in the embodiment of the present invention;

[0060] Figure 7 This is the pollution distribution map of the first stratum (0-0.5 m) in the embodiment of the present invention;

[0061] Figure 8 This is a structural block diagram of the contaminated soil mapping simulation system in an embodiment of the present invention. DETAILED DESCRIPTION

[0062] The following is a clear and complete description of the technical solutions of the present invention in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts are within the scope of protection of the present invention.

[0063] The following specific embodiments are used to describe the technical solution of the present application in detail. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.

[0064] The present invention conducts secondary development of ArcGIS software and uses the secondary developed ArcGIS software to simulate contaminated soil mapping, which greatly reduces the degree of manual participation and improves the efficiency of mapping simulation. Figure 1 As shown, a contaminated soil mapping simulation method provided by an embodiment of the present invention includes the following steps:

[0065] Step 1: Obtain soil pollution analysis results, and create a raw data table, an indicator threshold table, an over-standard analysis table, and a layer number-layer thickness comparison table based on the soil pollution analysis results;

[0066] Step 2: Import the exceedance analysis table and layer number-layer thickness comparison table into ArcGIS software;

[0067] Step 3: Based on the exceedance analysis table, filter out the exceedance analysis data of the corresponding stratum according to the layer number;

[0068] Step 4: Based on the excessive analysis data of the strata screened in step 3, the sampling points that have not been inspected in each pollution category are eliminated;

[0069] Step 5: Perform interpolation calculation on the exceeding standard analysis data of each pollution category processed in step 4 to obtain different exceeding standard attribute patches of the corresponding pollution category;

[0070] Step 6: Crop the patches with different exceeding attributes for each pollution category, remove the areas outside the red line boundary during the interpolation calculation, and merge the patches with the same exceeding attributes;

[0071] Step 7: Eliminate the patches that do not exceed the standard from the patches fused in step 6 to obtain the pollution distribution map of the corresponding pollution category;

[0072] Step 8: Overlay the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum;

[0073] Step 9: Repeat steps 3 to 8 to obtain the pollution distribution map of each stratum;

[0074] Step 10: For the pollution distribution map of each stratum, add the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, add the pollution combination type, calculate the pollution volume of each pollution category, and add the pollution patch number and inflection point.

[0075] In step 1, after obtaining the soil pollution analysis results, an original data table, an indicator threshold table, an exceeding standard analysis table, and a layer number-layer thickness comparison table are created in an Excel spreadsheet based on the soil pollution analysis results. The indicator threshold table contains the standard values or thresholds for each pollution indicator. These standard values or thresholds are determined based on the planned use of the soil, surrounding sensitivity, and economic factors in the risk assessment report. After determining the standard values or thresholds for each pollution indicator, an exceeding standard analysis table is created. When the standard values or thresholds of pollution indicators change according to the risk assessment report, the exceeding standard analysis table can be automatically updated to prevent the engineering quantities obtained by setting the standard values or thresholds during the risk assessment process from failing to meet economic requirements.

[0076] Figure 2 Shows part of the exceeding standard analysis table, Figure 2 It can be seen that the exceedance analysis table contains the sampling point number, sampling point coordinates, the layer number to which the sampling point belongs, the exceedance multiples of each pollution index, and the exceedance information of each pollution category. The exceedance information of each pollution category includes not submitted for inspection, not exceeded, exceedance category (that is, which pollution category exceeds the standard, for example Figure 2 polycyclic aromatic hydrocarbons in), maximum multiples exceeding the standard.

[0077] Pollution categories are the classification of pollution indicators based on possible remediation processes. Each pollution category contains at least one pollution indicator. For example, the simulation of contaminated soil mapping for a certain plot involves 7 pollution categories, namely petroleum hydrocarbons, benzene series, halogenated aliphatic hydrocarbons, polycyclic aromatic hydrocarbons, heavy metals, other VOCs, and other SVOCs. Petroleum hydrocarbons contain 1 pollution indicator, benzene series contain 11 pollution indicators, halogenated aliphatic hydrocarbons contain 21 pollution indicators, polycyclic aromatic hydrocarbons contain 9 pollution indicators, heavy metals contain 1 pollution indicator, other VOCs contain 3 pollution indicators, and other SVOCs contain 4 pollution indicators. Figure 3 shown.

[0078] When investigating soil pollution, stratification is performed from the ground to the depth where the soil exceeds the standard. This allows us to know the layer number to which each sampling point belongs and the thickness of the corresponding layer, such as Figure 4 As shown, the layer thickness is in meters. Figure 4 It can be seen that the pollution depth involved in this embodiment is 38 meters and the pollution layer is 18 layers.

[0079] Because soil and groundwater pollution problems within the plots are hidden, investigation plans formulated based on theories often have the following limitations during implementation:

[0080] (1) Investigation includes internal work (which generally refers to the work of project data and archive management, planning, statistical management, and internal secretarial management during project construction) and field work. Internal staff will revise the plan based on the results of the on-site investigation and the test data. However, the test data results have a certain time cycle. If the internal staff cannot revise the results in time, the field staff will still have to implement the plan with deviations or wait;

[0081] (2) The distribution of multiple types of soil pollutants is often spatially diverse. The investigation process needs to consider both the vertical (i.e., the drilling depth of each sampling point) and the horizontal (i.e., the distribution location of the sampling points) aspects simultaneously to ensure the accuracy of the pollution distribution data. However, due to the complexity of data processing and the implementation cycle, the actual investigation process cannot take into account both vertical and horizontal distribution as a whole, resulting in low accuracy of the pollution distribution range or an excessively large deep pollution distribution area, which in turn leads to excessively high subsequent remediation costs.

[0082] (3) During the on-site investigation, in order to avoid repeated drilling, drilling is often done by taking deeper samples, and whether to send the deeper samples for inspection is determined based on the test data. However, some detection methods are time-sensitive. If the internal work cannot respond in time, re-drilling may be required, resulting in a waste of manpower and costs.

[0083] (4) Data processing and analysis are generally done manually. This processing method requires manual real-time monitoring and calculation of project progress. The labor intensity is high and manual processing is prone to errors, which affects the construction progress.

[0084] In order to solve the problems in the above-mentioned traditional technologies that cannot take into account the vertical and horizontal distribution as a whole, resulting in low accuracy of pollution distribution range or too large a deep pollution distribution area, which in turn leads to excessively high subsequent repair costs, the industry's work cannot be handled in a timely manner, resulting in deviations in field work or idle work and re-drilling problems, as well as manual monitoring and calculation of project progress resulting in low work efficiency, prone to errors and high labor intensity, such as Figure 5 As shown, the present embodiment obtains the soil pollution analysis result specifically including the following steps:

[0085] Step 1.1: Obtain basic data.

[0086] Before analyzing the pollution distribution (both horizontally and vertically) and pollution indicators for a plot, a basic pollution distribution plan should be developed based on the pollution type and relevant national policies (such as GB_36600 and HJ·25.1 standards). Different pollution types correspond to different pollution indicators and require different distribution density.

[0087] The basic sampling point layout plan includes the initial sampling point layout plan and the encrypted sampling point layout plan. The initial sampling point layout plan includes the layout location and density of the initial sampling points, the planned drilling depth, the point number, the inspection depth (i.e., the sampling depth of the inspection sample), the pollution index (i.e., the index that the inspection sample needs to be tested) and the index detection method. For example, the density of the initial sampling points is a 40m*40m grid, and the boreholes can be layered from top to bottom, for example, 0~0.5m is the first layer, 0.5m~1m is the second layer, 1m~2m is the third layer, ..., different inspection depths correspond to different layers.

[0088] The intensified sampling point layout plan includes the location, density, and numbering of intensified sampling points. The drilling depth, inspection depth, pollution indicators, and indicator detection methods of the intensified sampling points are determined through perimeter enclosure analysis of the intensified sampling points. According to relevant national policies on the density of intensified sampling points for contaminated land, intensified sampling points include primary intensified sampling points, secondary intensified sampling points, and so on, and N-th intensified sampling points, where N is determined by the density requirements. For example, if the point density requirement is a 20m*20m grid, when the density of the initial sampling points is a 40m*40m grid and the density of the primary encrypted sampling points is a 20m*20m grid, the density of the primary encrypted sampling points already meets the point density requirement, and the encrypted sampling points only include the primary encrypted sampling points; if the point density requirement is a 10m*10m grid, when the density of the initial sampling points is a 40m*40m grid, the density of the primary encrypted sampling points is a 20m*20m grid, and the density of the secondary encrypted sampling points is a 10m*10m grid, the density of the secondary encrypted sampling points meets the point density requirement, and the encrypted sampling points include the primary encrypted sampling points and the secondary encrypted sampling points. This embodiment takes the case where the encrypted sampling points only include the primary encrypted sampling points as an example, that is, the encrypted sampling points are the primary encrypted sampling points.

[0089] When the plot to be investigated or analyzed is large, the plot is divided into regions, and the contaminated soil analysis method of the present invention is performed on each region. After the contaminated soil analysis method of the present invention is completed in one region, the contaminated soil analysis method of the present invention is performed on the next region, and the contaminated soil analysis method of the present invention is performed on all regions in the plot, thereby obtaining complete analysis data for the plot. When the plot to be investigated or analyzed is small, the contaminated soil analysis method of the present invention can also be performed on the entire plot as a unit.

[0090] Basic data will be obtained during the implementation of different plans:

[0091] During the implementation of the basic site layout plan in the field, soil samples are taken at the initial sampling points, unified rules are formulated for the sample names (such as point number-depth-others), and the samples are sent to the testing unit for testing. The testing unit gradually feeds back the test results of each pollution indicator, that is, the basic data in the implementation of the basic site layout plan is obtained (for distinction, referred to as the first basic data). The first basic data includes the initial sampling point number, the inspection depth, the actual sampling position (that is, the actual sampling depth of the sample submitted for inspection), the test results of each pollution indicator in the sample, and the maximum sampling depth of each initial sampling point.

[0092] During the implementation of the deepening plan for the initial sampling points, the initial sampling points are sampled more deeply, and the samples are sent to the testing unit for testing. The testing unit gradually feeds back the test results of each pollution index, that is, the basic data (referred to as the second basic data for the sake of distinction) during the implementation of the deepening plan for the initial sampling points are obtained, where the second basic data include the initial sampling point number, the inspection depth, the actual sampling position (that is, the actual sampling depth of the sample submitted for inspection), the test results of each pollution index in the sample, and the maximum sampling depth of each initial sampling point.

[0093] When conducting a fallback analysis based on the results of the second basic data's indicator exceeding standard analysis, if the fallback is not achieved, it is necessary to determine a further deepening plan for the initial sampling point, and then obtain the basic data, exceeding standard analysis, and fallback analysis during the implementation of the further deepening plan until the fallback is achieved.

[0094] During the implementation of the encrypted sampling point layout plan, samples are taken at the encrypted sampling points and the samples are sent to the testing units for testing. The testing units gradually feedback the test results of various pollution indicators, that is, the basic data during the implementation of the encrypted sampling point layout plan is obtained (for the sake of distinction, referred to as the third basic data). The third basic data includes the encrypted sampling point number, the inspection depth, the actual sampling position (that is, the actual sampling depth of the sample submitted for inspection), the test results of each pollution indicator in the sample, and the maximum sampling depth of each encrypted sampling point.

[0095] When conducting a fallback analysis based on the results of the third basic data's indicator exceeding the standard analysis, if the fallback is not achieved, it is necessary to determine a deepening plan for the encrypted sampling points, and then obtain the basic data during the implementation of the deepening plan for the encrypted sampling points. Similar to the deepening implementation plan for the initial sampling points, repeat the steps of exceeding the standard analysis and fallback analysis until the fallback is achieved.

[0096] The encrypted sampling points in this embodiment only include one-time encrypted sampling points. If the encrypted sampling points include multiple-time encrypted sampling points, each encrypted sampling point will execute the steps of determining the corresponding secondary encrypted sampling point layout plan by perimeter enclosure analysis and determining the encryption plan of the corresponding secondary encrypted sampling point by bottom-up analysis.

[0097] Step 1.2: Analysis of pollution index exceeding the standard.

[0098] Corresponding to step 1.1, the pollution index exceeding standard analysis includes performing the index exceeding standard analysis based on the test results of each pollution index in the basic data of different schemes:

[0099] After obtaining the first basic data, the detection results of each pollution indicator in the first basic data are analyzed for index exceeding the standard to obtain the index exceeding the standard analysis result (for the sake of distinction, referred to as the first exceeding the standard analysis result): the detection result of the pollution indicator is compared with the corresponding indicator threshold. If it exceeds the indicator threshold, the pollution indicator exceeds the standard; if it does not exceed the indicator threshold, the pollution indicator does not exceed the standard.

[0100] After obtaining the second basic data, an indicator exceedance analysis is performed on the test results of each pollution indicator in the second basic data to obtain an indicator exceedance analysis result (for the sake of distinction, simply referred to as the second exceedance analysis result). If a re-deepening plan exists for the initial sampling point, an indicator exceedance analysis is also performed on the test results of each pollution indicator in the basic data during the implementation of the re-deepening plan.

[0101] After obtaining the third basic data, an index exceedance analysis is performed on the test results of each pollution indicator in the third basic data to obtain the index exceedance analysis results (for the sake of distinction, simply referred to as the third exceedance analysis results). If a deepening plan or a further deepening plan with denser sampling points exists, an index exceedance analysis is also performed on the test results of each pollution indicator in the basic data during the implementation of the deepening plan or even the further deepening plan.

[0102] If there are multiple encrypted sampling points, similarly, the detection results of various pollution indicators in the basic data of the implementation process of each encrypted sampling point layout plan, deepening plan, and re-deepening plan will be analyzed for indicators exceeding the standard.

[0103] The test results of pollution indicators fed back by the testing units are usually displayed in the form of Excel tables or Word, and the test results are usually expressed in text format. In order to facilitate identification in the subsequent processing process, the undetected pollution indicators, undetected pollution indicators and exceeded pollution indicators are digitally identified according to the results of the indicator exceeding the standard analysis, that is, the test results in text format are automatically changed to digital identification. For example, if the pollution indicator is an undetected indicator (that is, the indicator has not been tested), it is automatically identified as "-2" (corresponding to the Excel table, it may be displayed as " / " or "-"); if the pollution indicator is an undetected indicator (that is, the content of the indicator is not detected in the sample or the indicator is not present), it is automatically identified as "0"; if the pollution indicator is an exceeded indicator, it is automatically identified as the multiple of the exceeded standard.

[0104] When marking the pollution index, the sampling point number and the inspection depth must also be marked. For example, the sampling point number-inspection depth.

[0105] Step 1.3: Comprehensive analysis of pollution indicators.

[0106] The bottom-line analysis is used to determine whether it is necessary to conduct a deeper test (i.e., a deepening plan or a re-deepening plan) based on the sampling point layout plan. The specific implementation process of the bottom-line analysis is: determine the maximum exceeding depth and the maximum inspection depth based on the indicator exceeding analysis results;

[0107] If the maximum exceeding depth is equal to the maximum inspection depth, it indicates that the bottom is not covered. If the bottom is not covered and no samples are retained, the initial sampling point corresponding to the non-covered factor shall be drilled deeper; if the bottom is not covered and samples are retained, there is no need to deepen the drilling and the retained samples shall be sent for inspection; if the maximum exceeding depth is less than the maximum inspection depth, it indicates that the bottom is covered.

[0108] The maximum exceeding depth and maximum inspection depth of each sampling point can be automatically obtained through the point number. The exceeding depth refers to the inspection depth corresponding to the exceeding pollution index.

[0109] If the maximum exceeding depth is less than the maximum inspection depth, it indicates that there is no pollution at the maximum inspection depth, which is a guarantee and does not require further deepening for pollution analysis; if the maximum exceeding depth is equal to the maximum inspection depth, it indicates that there is still pollution at the maximum inspection depth (i.e., there are pollution indicators exceeding the standard), which is not a guarantee and requires further drilling to further determine the pollution depth. The pollution indicators exceeding the standard are the unguaranteed factors. When there is no guarantee, it is necessary to deepen the test, that is, to determine the deepening plan or deepen the plan again. During the implementation of the plan, if the maximum sampling depth is greater than the maximum inspection depth, it indicates that there are retained samples (i.e., samples have been taken but not sent for inspection), and there is no need to deepen the drilling of the corresponding sampling points (i.e., the sampling points corresponding to the pollution indicators exceeding the standard). The retained samples can be directly sent for inspection; if the maximum sampling depth is equal to the maximum inspection depth, it indicates that there are no retained samples and it is necessary to deepen the drilling of the corresponding sampling points, take samples, and then send the samples for inspection.

[0110] Corresponding to step 1.2, the pollution index bottom-up analysis includes the following:

[0111] Based on the first exceeding standard analysis result, an indicator bottom-up analysis is performed to obtain the indicator bottom-up analysis result (for distinction, referred to as the first bottom-up analysis result); when there is no bottom-up, the deepening plan of the initial sampling point is determined according to the first bottom-up analysis result; when the bottom-up is achieved, all exceeding standard analysis results of the initial sampling point are obtained.

[0112] Based on the second exceedance analysis results, a fallback analysis is performed to obtain the fallback analysis results (referred to as the second fallback analysis results for clarity). If the fallback is not achieved, a further deepening plan for the initial sampling point is determined based on the second fallback analysis results. If the fallback is achieved, all exceedance analysis results for the initial sampling point are obtained. Repeat the acquisition of basic data, the exceedance analysis, and the fallback analysis for the initial sampling point for the further deepening plan until the fallback is achieved.

[0113] Based on the third exceedance analysis results, perform an indicator bottom-line analysis to obtain the indicator bottom-line analysis results (for distinction, referred to as the third bottom-line analysis results). If the bottom line is not reached, determine the deepening plan for the intensified sampling points based on the third bottom-line analysis results. If the bottom line is reached, obtain all exceedance analysis results for the intensified sampling points. Repeat the acquisition of basic data, indicator exceedance analysis, and indicator bottom-line analysis for the deepening plan for the intensified sampling points until the bottom line is reached.

[0114] If the encrypted sampling points include single-time encrypted sampling points, secondary encrypted sampling points, ..., N-times encrypted sampling points, an indicator backup analysis shall be performed on the exceeding standard analysis results corresponding to each encrypted sampling point layout plan, and the acquisition of basic data, indicator exceeding standard analysis, and indicator backup analysis shall be repeated until the backup is achieved.

[0115] Step 1.4: Conduct perimeter enclosure analysis of the encrypted sampling points.

[0116] With the implementation of the initial sampling point layout plan, deepening plan, and even further deepening plan, the drilling depths of different sampling points are different, which leads to deeper or larger sampling layers and a higher density of sampling points (that is, the deeper the sampling, the fewer sampling points in the corresponding layer), resulting in a large amount of restoration work, so a deep perimeter enclosure analysis is needed.

[0117] The basic point distribution plan specifies the location, density and point number of the encrypted sampling points, but does not specify the drilling depth, inspection depth, pollution index and index detection method of the encrypted sampling points. The encrypted sampling point distribution plan is determined by conducting perimeter enclosure analysis on the encrypted sampling points, that is, determining the drilling depth, inspection depth, pollution index and index detection method of each encrypted sampling point.

[0118] If the encrypted sampling points include primary encrypted sampling points, secondary encrypted sampling points, ..., N secondary encrypted sampling points, a perimeter shielding analysis is performed on each encrypted sampling point in sequence to determine the corresponding encrypted sampling point layout plan. For example, a perimeter shielding analysis is performed on each primary encrypted sampling point to determine the primary encrypted sampling point layout plan; a perimeter shielding analysis is performed on each secondary encrypted sampling point to determine the secondary encrypted sampling point layout plan, and so on.

[0119] In this embodiment, the specific implementation process of the perimeter enclosure analysis is as follows:

[0120] Take any encrypted sampling point as the target point, and the initial sampling points around the target point as the pending points, enter the point code of the pending point, obtain and update the exceeding standard analysis results of the pending point according to the point code, and then obtain and update the exceeding standard analysis results of the pending point when different encrypted sampling points are used as the target point; determine the drilling depth, inspection depth, pollution index and index detection method of the target point based on the exceeding standard analysis results of all the pending points of each target point.

[0121] Since both the initial sampling points and the intensified sampling points are arranged in a grid pattern, the sampling points surrounding the target point refer to the four initial sampling points surrounding the target point. If the intensified sampling points include primary intensified sampling points, secondary intensified sampling points, ..., N-th intensified sampling points, when the target point is an intensified sampling point in the i-th intensified sampling point set, the target point's pending points are the four intensified sampling points in the i-1th intensified sampling point set, where i = 1, 2, ..., N. When i = 1, the target point's pending points are the four initial sampling points surrounding the target point.

[0122] Since the implementation of the plan takes a long time, in order to narrow the cycle time difference, after gradually obtaining the exceedance analysis results according to the point number, the drilling depth, inspection depth, pollution index, and index detection method of the target point can be determined based on the exceedance analysis results of the pending point. That is, when the first exceedance analysis results are gradually obtained according to the point number (without waiting for the second exceedance analysis results to be updated), the first intensified sampling point layout plan is determined based on the first exceedance analysis results; when the exceedance analysis results corresponding to the first intensified sampling point layout plan are gradually obtained according to the point number, the second intensified sampling point layout plan is determined based on the exceedance analysis results, and so on, without waiting for the exceedance analysis results corresponding to the deepening plan or even the second deepening plan of the sampling point to be updated.

[0123] After obtaining the exceeding standard analysis results step by step according to the point number, the status of each pending point is converted into a summary of information such as the sampling point number, pollution index and exceeding standard analysis results. By classifying according to the pollution index, it can be clearly shown which inspection depth positions of each pollution index around the target point exceed the standard, thereby ensuring that the inspection depth (or inspection level) of the pollution index of the target point includes all exceeding standard inspection depths of the surrounding sampling points. Finally, the inspection depths of each pollution index of the target point are summarized, and the maximum sampling depth of the target point is obtained, and an encrypted sampling point layout plan is automatically formed.

[0124] The inspection depth of the target point is equal to the union of the inspection depths of all the pending points with exceeded standards, the pollution index of the target point is equal to the union of the exceeded standards of all the pending points, and the drilling depth of the target point is greater than the maximum inspection depth of the target point.

[0125] For example, the pollution indexes exceeding the standard and the exceeding depths of the four pending points of the target point are as follows: pollution index A exceeds the standard and the exceeding depth is L1 (pending point 1), pollution index B exceeds the standard and the exceeding depth is L2 (pending point 2), pollution index C exceeds the standard and the exceeding depth is L3 (pending point 3), pollution index D exceeds the standard and the exceeding depth is L4 (pending point 4), then the inspection depths of the target point are L1, L2, L3 and L4, the pollution indexes (or inspection indexes) of the target point are pollution index A, pollution index B, pollution index C and pollution index D, the target point The drilling depth of the point is greater than the maximum value of the inspection depths L1, L2, L3 and L4, that is, sampling is taken at the depth L1 of the target point, and the pollution index A, pollution index B, pollution index C and pollution index D are tested; sampling is taken at the depth L2 of the target point, and the pollution index A, pollution index B, pollution index C and pollution index D are tested; sampling is taken at the depth L3 of the target point, and the pollution index A, pollution index B, pollution index C and pollution index D are tested; sampling is taken at the depth L4 of the target point, and the pollution index A, pollution index B, pollution index C and pollution index D are tested.

[0126] In the basic point distribution plan, the number of encrypted sampling times N included in the encrypted sampling points has been determined according to the point distribution density requirements, that is, the encrypted sampling points include single-time encrypted sampling points, secondary-time encrypted sampling points, ..., N-times encrypted sampling points. After the N-times encrypted sampling points are subjected to perimeter enclosure analysis to determine the N-times encrypted sampling point distribution plan, the steps of obtaining basic data, exceeding-standard index analysis, and bottom-line index analysis are repeated until the bottom line is reached. Then, the point distribution density of the area or plot to be analyzed meets the point distribution density requirements stipulated by relevant national policies.

[0127] The present invention dynamically obtains basic data during the implementation of each plan, conducts an over-standard analysis on the indicator detection results in the basic data, conducts a bottom-up analysis and a peripheral enclosure analysis based on the over-standard analysis results, and dynamically corrects the layout plan of the initial sampling points and the encrypted sampling points through the bottom-up analysis, that is, deepens the plan, and automatically forms an encrypted sampling point layout plan through the peripheral enclosure analysis, so that the layout plan is more accurate and reasonable, the on-site management is smoother, the data analysis time is greatly shortened, the cycle time difference between internal and external work is reduced, and the problem of deviation in external work or idle waiting and re-drilling caused by the inability to handle internal work in a timely manner is improved, thereby reducing manpower and material costs; the present invention greatly improves the accuracy of data processing, reduces the probability of errors in human analysis, and makes the final data results more scientific and objective.

[0128] The basic data of the present invention includes the location of each sampling point and the inspection depth. At the same time, considering the horizontal distribution and vertical depth, the accuracy of the soil pollution distribution data is ensured, and thus the accuracy of the pollution distribution range analysis is ensured.

[0129] In order to facilitate the preservation of the process files of steps 4 to 9, use ArcGIS software to create a folder for saving process files under the operation directory. The process files of steps 4 to 9 will be saved in this folder.

[0130] In order to avoid errors when ArcGIS software recognizes the layer number-layer thickness comparison table, ArcGIS software is used to convert the layer number-layer thickness comparison table into bbf format.

[0131] In step 2, the exceedance analysis table and the layer number-layer thickness comparison table are imported into ArcGIS software, and the interpolation calculation method and clipping elements (i.e., the red line range input by the clipping module) are set to automatically run steps 2 to 9 to obtain the pollution distribution map of different strata, such as Figure 6 As shown. Figure 6 As can be seen, this embodiment uses Thiessen polygons to interpolate the excess pollution analysis data for each pollution category. Each sampling point corresponds to a Thiessen polygon. After interpolation, patches with different excess pollution attributes are obtained for each pollution category. The patches are clipped according to the red line range and fused together to form a single patch instead of separate patches.

[0132] When the exceeding attribute includes exceeding and not exceeding the standard, the interpolation calculation results in the exceeding small patches (i.e., patches corresponding to the sampling points that exceed the standard) and not exceeding small patches (i.e., patches corresponding to the sampling points that do not exceed the standard) for each pollution category. All the exceeding small patches are fused, and all the not exceeding small patches are fused to obtain the exceeding patches and not exceeding patches for each pollution category.

[0133] When the exceedance attribute includes different exceedance multiples, interpolation calculations yield small patches with different exceedance multiples for each pollution category. Small patches with the same exceedance multiple are fused to obtain patches with each exceedance multiple for each pollution category. For example, all patches with an exceedance multiple of 0.5 are fused, and all patches with an exceedance multiple of 0 are fused, where the exceedance multiple corresponding to no exceedance is 0. This fusion process can reduce the number of boundary lines drawn at sampling points with the same exceedance attribute, reducing inflection points.

[0134] In step 3, the exceeding standard analysis data of different strata include the sampling point number of the stratum, the sampling point coordinates, the exceeding standard multiples of each pollution index and the exceeding standard information of each pollution category.

[0135] Repeating steps 4 to 7 can obtain the pollution distribution map of different pollution categories in a certain stratum. Repeating steps 3 to 8 can obtain the pollution distribution map of strata corresponding to different layer numbers, such as Figure 7 The pollution distribution map of the first layer is shown. In step 8, the layer number is lost during overlay, and a layer number value is read at the same time. Therefore, it is necessary to add a layer number to the pollution distribution map of each stratum. The layer number can be added by reading the layer number value. The corresponding layer thickness can be added according to the layer number-layer thickness comparison table and the layer number; the pollution combination type (i.e., pollution category) can be added according to the exceedance analysis table.

[0136] like Figure 8 As shown, an embodiment of the present invention further provides a contaminated soil mapping simulation system, comprising:

[0137] An import module is used to import an exceeding standard analysis table and a layer number-layer thickness comparison table. The exceeding standard analysis table and the layer number-layer thickness comparison table are established based on the soil pollution analysis results. The exceeding standard analysis table includes the sampling point number, the sampling point coordinates, the layer number to which the sampling point belongs, the exceeding standard multiple of each pollution indicator, and the exceeding standard information of each pollution category;

[0138] A screening module, configured to screen out the exceeding standard analysis data of the corresponding stratum according to the layer number based on the exceeding standard analysis table;

[0139] A first elimination module is used to eliminate the sampling points that have not been inspected in each pollution category based on the excessive analysis data of the stratum screened by the screening module;

[0140] An interpolation calculation module, configured to perform interpolation calculation on the excessive-standard analysis data of each pollution category processed by the first elimination module to obtain different excessive-standard attribute patches corresponding to the pollution category;

[0141] The clipping and fusion module is used to clip patches with different excessive attributes for each pollution category and fuse patches with the same excessive attributes. The clipping and fusion module needs to cancel the creation of multi-component features.

[0142] A second elimination module is used to eliminate the patches that do not exceed the standard from the patches processed by the clipping and fusion module to obtain a pollution distribution map corresponding to the pollution category;

[0143] an overlay module, configured to overlay the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum;

[0144] The field adding and calculating module is used for the pollution distribution map of each stratum, adding the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, adding the pollution combination type, calculating the pollution volume of each pollution category, and adding the pollution patch number and inflection point.

[0145] The system also includes a process design module for determining a remediation process based on a pollution distribution map of each stratum, the number of times the pollution exceeds the standard, the type of pollution combination, and the stratum number.

[0146] The number of sampling points exceeded 1,000, the depth exceeding the standard was 27m, the number of pollution combination types was 31, the number of soil layers was 15, and 500,000 data were involved. After setting the indicator thresholds, the exceedance analysis table and layer number-layer thickness comparison table were imported, and the interpolation calculation method and cropping parameters were set. The system only took 30 minutes to run and realize the mapping simulation of 15 layers.

[0147] An embodiment of the present invention further provides an electronic device, comprising: a processor and a memory storing a computer program, wherein the processor is configured to implement the above-mentioned contaminated soil mapping simulation method when executing the computer program.

[0148] Although not shown, the electronic device includes a processor that can perform various appropriate operations and processes based on programs and / or data stored in a read-only memory (ROM) or programs and / or data loaded from a storage portion into a random access memory (RAM). The processor can be a multi-core processor or can include multiple processors. In some embodiments, the processor can include a general-purpose main processor and one or more special coprocessors, such as a central processing unit (CPU), a graphics processing unit (GPU), a neural network processor (NPU), a digital signal processor (DSP), etc. Various programs and data required for the operation of the electronic device are also stored in the RAM. The processor, ROM, and RAM are connected to each other via a bus. An input / output (I / O) interface is also connected to the bus.

[0149] The processor and memory are used together to execute the program stored in the memory. When the program is executed by the computer, the methods, steps or functions described in the above embodiments can be implemented.

[0150] Although not shown, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the contaminated soil mapping simulation method as described above is implemented.

[0151] Storage media in embodiments of the present invention include permanent and non-permanent, removable and non-removable items that can be used to store information using any method or technology. Examples of storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device.

[0152] The above disclosure is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any technician familiar with this technical field can easily think of changes or modifications within the technical scope disclosed in the present invention, and they should all be covered by the scope of protection of the present invention.

Claims

1. A contaminated soil mapping simulation method, characterized in that: The method comprises the following steps: Step 1: Obtain soil pollution analysis results, and create a raw data table, an indicator threshold table, an exceeding standard analysis table, and a layer number-layer thickness comparison table based on the soil pollution analysis results; wherein the exceeding standard analysis table includes the sampling point number, sampling point coordinates, the layer number to which the sampling point belongs, the exceeding standard multiple of each pollution indicator, and the exceeding standard information of each pollution category; Step 2: Import the above-mentioned excess analysis table and layer number-layer thickness comparison table into ArcGIS software; Step 3: Based on the exceeding standard analysis table, filter out the exceeding standard analysis data of the corresponding stratum according to the layer number; Step 4: Based on the excessive analysis data of the strata screened out in step 3, the sampling points that have not been inspected in each pollution category are eliminated; Step 5: performing interpolation calculation on the exceeding standard analysis data of each pollution category processed in step 4 to obtain different exceeding standard attribute patches corresponding to the pollution category; Step 6: Crop the patches with different exceeding attributes for each pollution category and fuse the patches with the same exceeding attributes; Step 7: Eliminate the patches that do not exceed the standard from the patches fused in step 6 to obtain a pollution distribution map corresponding to the pollution category; Step 8: superimposing the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum; Step 9: Repeat steps 3 to 8 to obtain the pollution distribution map of each stratum; Step 10: For the pollution distribution map of each stratum, add the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, add the pollution combination type, calculate the pollution volume of each pollution category, and add the pollution patch number and inflection point; Wherein, in said step 6, when the exceeding attribute includes exceeding the standard and not exceeding the standard, the patches corresponding to all the exceeding sampling points are fused, and the patches corresponding to all the not exceeding sampling points are fused to obtain the exceeding patches and not exceeding the standard patches for each pollution category; When the exceeding standard attribute includes different exceeding standard multiples, the patches corresponding to the sampling points with the same exceeding standard multiples are merged to obtain the exceeding standard multiples patches for each pollution category; the exceeding standard multiple corresponding to no exceeding standard is 0.

2. The contaminated soil mapping simulation method according to claim 1, characterized in that: In step 1, obtaining soil pollution analysis results specifically includes: Obtain basic data during the implementation of the basic sampling plan, wherein the basic data includes the initial sampling point number, the inspection depth, the actual sampling location, the test results of each pollution index in the sample, and the maximum sampling depth of each initial sampling point; Conduct an index exceeding standard analysis based on the test results of each pollution index, and mark the untested pollution index, undetected pollution index and pollution index exceeding the standard according to the index exceeding standard analysis results; Perform an indicator bottom-up analysis based on the indicator exceeding the standard analysis results. If the bottom-up analysis is not successful, determine the deepening plan for the initial sampling point based on the indicator bottom-up analysis results; repeat the steps of obtaining basic data, indicator exceeding the standard analysis, and indicator bottom-up analysis during the implementation of the deepening plan until the bottom-up analysis is successful, and obtain the analysis results of all indicator exceeding the standard; Conduct perimeter enclosure analysis on the encrypted sampling points in the basic layout plan to determine the encrypted sampling point layout plan, wherein the encrypted sampling point layout plan includes the drilling depth, inspection depth, pollution index and index detection method of each encrypted sampling point; repeat the steps of obtaining basic data, index exceeding standard analysis and index guarantee analysis during the implementation of the encrypted sampling point layout plan.

3. The contaminated soil mapping simulation method according to claim 2, characterized in that: The specific implementation process of the indicator bottom-up analysis is as follows: Determine the maximum exceeding depth and the maximum inspection depth according to the exceeding analysis results of the indicators; If the maximum exceeding depth is equal to the maximum inspection depth, it indicates that the bottom is not covered. If the bottom is not covered and no samples are retained, the initial sampling point corresponding to the non-coverage factor will be drilled deeper. If the bottom is not covered and samples are retained, there is no need to drill deeper and the retained samples will be sent for inspection. If the maximum exceeding depth is less than the maximum inspection depth, it indicates a safety net.

4. The contaminated soil mapping simulation method according to claim 2, characterized in that: The specific implementation process of the perimeter enclosure analysis is as follows: Taking any encrypted sampling point as the target point, taking the initial sampling points around the target point as the pending points, entering the point code of the pending point, obtaining and updating the excess analysis results of the pending point according to the point code, and then obtaining and updating the excess analysis results of the pending point when different encrypted sampling points are used as the target points; Based on the exceeding standard analysis results of all pending points at each target point, the drilling depth, inspection depth, pollution index and index detection method of the target point are determined.

5. The contaminated soil mapping simulation method according to claim 1, characterized in that: In step 5, Thiessen polygons are used to perform interpolation calculation on the exceeding standard analysis data of each pollution category.

6. The contaminated soil mapping simulation method according to claim 1, characterized in that: Before step 2, the method further includes using ArcGIS software to create a folder for storing process files under the operation directory.

7. A contaminated soil mapping simulation system, characterized in that: The system is obtained by secondary development of ArcGIS software, and the system includes: An import module is used to import an exceeding standard analysis table and a layer number-layer thickness comparison table. The exceeding standard analysis table and the layer number-layer thickness comparison table are established based on the soil pollution analysis results. The exceeding standard analysis table includes the sampling point number, the sampling point coordinates, the layer number to which the sampling point belongs, the exceeding standard multiple of each pollution indicator, and the exceeding standard information of each pollution category; A screening module, configured to screen out the exceeding standard analysis data of the corresponding stratum according to the layer number based on the exceeding standard analysis table; A first elimination module is used to eliminate the sampling points that have not been inspected in each pollution category based on the excessive analysis data of the stratum screened by the screening module; An interpolation calculation module, configured to perform interpolation calculation on the excessive-standard analysis data of each pollution category processed by the first elimination module to obtain different excessive-standard attribute patches corresponding to the pollution category; The cropping and fusion module is used to crop the patches of different exceeding attributes of each pollution category and fuse the patches of the same exceeding attributes. When the exceeding attributes include exceeding and not exceeding, the patches corresponding to all exceeding sampling points are fused, and the patches corresponding to all not exceeding sampling points are fused to obtain the exceeding and not exceeding patches of each pollution category. When the exceeding attributes include different exceeding multiples, the patches corresponding to sampling points with the same exceeding multiple are fused to obtain the patches of each exceeding multiple for each pollution category. The exceeding multiple corresponding to not exceeding is 0. A second elimination module is used to eliminate the patches that do not exceed the standard from the patches processed by the clipping and fusion module to obtain a pollution distribution map corresponding to the pollution category; an overlay module, configured to overlay the pollution distribution maps of all pollution categories in the stratum to form a pollution distribution map of the corresponding stratum; The field adding and calculating module is used for the pollution distribution map of each stratum, adding the layer number and corresponding layer thickness according to the layer number-layer thickness comparison table, adding the pollution combination type, calculating the pollution volume of each pollution category, and adding the pollution patch number and inflection point.

8. An electronic device, characterized in that: The device comprises: Memory for storing computer programs; A processor, configured to implement the contaminated soil mapping simulation method according to any one of claims 1 to 6 when executing the computer program.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the contaminated soil mapping simulation method according to any one of claims 1 to 6 is implemented.

Images