Multi-source fusion data-based urban industrial pollution site grading management and control method

Abstract

The invention discloses an urban industrial pollution site grading management and control method based on multi-source fusion data, relates to the technical field of pollution site management and control, and solves the technical problems of inaccurate analysis caused by single-dimensional analysis and low management and control efficiency caused by one-step management and control measures. According to the method, the pollution characteristics are described from four dimensions of pollution complexity, extreme concentration, overall strength and space scale by integrating four core indexes of pollutant standard exceeding quantity, maximum standard exceeding multiple, comprehensive standard exceeding multiple and pollution volume, complex scenes of multi-pollutant combined pollution, low-concentration large-range pollution and the like can be covered, the risk level misjudgment rate is reduced, and the risk level accuracy is improved. Quantitative assignment standards of the four indexes are defined, a weight interval is set, subjective experience interference is reduced, large-scale management requirements are met, differential schemes are made for low, medium, high and extremely high risk levels, and compared with existing one-step management and control, the low-risk site treatment cost can be reduced.

Description

Technical Field

[0001] This invention relates to the field of contaminated site management technology, specifically a hierarchical management method for urban industrial contaminated sites based on multi-source fusion data. Background Technology

[0003] According to patent application number 202210269818.8, a risk management method for contaminated sites and an in-situ risk management domain for contaminated sites are disclosed. The in-situ risk management domain for contaminated sites includes a defined in-situ reaction zone and injection sites and / or extraction wells set within the in-situ reaction zone. The injection sites are either injection sites or direct injection sites. The method for determining the number of injection sites and / or extraction wells is as follows: the number of injection sites and / or extraction wells is obtained based on the basic parameters of the in-situ reaction zone, the design length of the in-situ reaction zone, the influence radius of the injection sites and / or extraction wells, and the layout method.

[0004] However, in the existing technology, the industry generally adopts the basic logic of pollutant concentration detection + risk level determination, but there are obvious limitations in practice. The selection of pollution indicators often focuses on a single dimension, ignoring the comprehensive impact of other indicators. Moreover, the control measures after classification are mostly one-size-fits-all, without combining risk level to formulate differentiated and implementable management plans, resulting in insufficient control of high-risk sites or over-treatment of low-risk sites. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this invention provides a hierarchical management method for urban industrial contaminated sites based on multi-source fusion data. This method solves the problems of inaccurate analysis due to single-dimensional analysis and low management efficiency due to one-size-fits-all management measures.

[0006] To achieve the above objectives, the present invention provides the following technical solution: a method for graded management and control of urban industrial contaminated sites based on multi-source fusion data, which specifically includes the following steps: Multi-source data of urban industrial contaminated sites are collected, processed and cleaned, and the cleaned multi-source data are merged. At the same time, four pollution indicators are extracted: the number of pollutants exceeding the standard, the maximum multiple of exceeding the standard, the comprehensive multiple of exceeding the standard, and the volume of pollution. The number of pollutants exceeding the standard and the maximum multiple of exceeding the standard are assigned values ​​respectively. Different levels of values ​​are assigned based on the number of pollutants exceeding the standard, and values ​​are assigned based on the multiple relationship between the maximum multiple of exceeding the standard and the corresponding evaluation standard value. Assign values ​​to the comprehensive exceedance multiple and the pollution volume, calculate and assign values ​​to the comprehensive exceedance multiple of all pollutants exceeding the standard, and calculate and assign values ​​to the total pollution volume; The weighted sum of the values ​​assigned to different pollution indicators is used to obtain a comprehensive score, which is then matched with the evaluation criteria to determine the corresponding pollution level. Different control measures are implemented based on the pollution levels of low risk, medium risk, high risk, and extremely high risk, generating pollution control information.

[0007] As a further aspect of the present invention, the multi-source data includes soil pollution data and groundwater pollution data, and the soil pollution data includes the pollutant concentration of the soil sample, the sampling depth, and the sampling location coordinates. Groundwater pollution data includes pollutant concentrations, pH values, water level depth, water flow direction, and aquifer type in groundwater samples.

[0008] As a further aspect of the present invention, the method for assigning values ​​to the number of pollutants exceeding the standard is as follows: count the number of pollutant types whose concentrations exceed the risk screening value for the application site type; if the number of types is 1, it is judged as mild and assigned a value of 1; if the number of types is 2-3, it is judged as moderate and assigned a value of 2; if the number of types is ≥4, it is judged as severe and assigned a value of 3. The method for assigning the maximum exceedance multiple is as follows: calculate the individual exceedance multiple of each pollutant according to the formula (measured concentration - evaluation standard value) / evaluation standard value and take the maximum value. If the measured concentration corresponding to the maximum value is ≤2×evaluation standard value, it is judged as mild and assigned a value of 1; if 2×evaluation standard value < measured concentration ≤4×evaluation standard value, it is judged as moderate and assigned a value of 2; if 4×evaluation standard value < measured concentration ≤11×evaluation standard value, it is judged as severe and assigned a value of 3; if the measured concentration >11×evaluation standard value, it is judged as extreme and assigned a value of 4.

[0009] As a further aspect of the present invention, the method for assigning the comprehensive exceedance multiple is as follows: calculate the average exceedance multiple of each pollutant and the proportion of exceedance points (proportion of exceedance points = number of exceedance sampling points for that pollutant / total number of sampling points for that pollutant), and calculate the comprehensive exceedance multiple according to the formula: comprehensive exceedance multiple = Σ(average exceedance multiple of a pollutant × proportion of exceedance points for that pollutant) / number of pollutant types; set an exceedance threshold, if the comprehensive exceedance multiple ≤ 0.5 × exceedance threshold, it is judged as mild and assigned a value of 1; if 0.5 × exceedance threshold < comprehensive exceedance multiple ≤ 1 × exceedance threshold, it is judged as moderate and assigned a value of 2; if 1 × exceedance threshold < comprehensive exceedance multiple ≤ 3 × exceedance threshold, it is judged as severe and assigned a value of 3; if the comprehensive exceedance multiple ≥ 3 × exceedance threshold, it is judged as extreme and assigned a value of 4.

[0010] As a further aspect of the present invention, the method for assigning the pollution volume is as follows: k vertical cross-sections are arranged at intervals L along the pollution diffusion direction; the pollution depth range of each cross-section is determined by plotting a depth-concentration curve based on borehole data; and the pollution area S of each cross-section is calculated using a geometric formula in conjunction with the horizontal width of the cross-section; the total pollution volume is calculated according to the formula = (S1×L+S2×L+…+S) k-1The total pollution volume is calculated using ×L, where k-1 represents the number of corresponding cross sections. If the total pollution volume is ≤1000 m³, ... 3 If it is determined to be mild and assigned a value of 1; if 1000m 3 Total pollution volume ≤ 10000 m³ 3 It is classified as moderate and assigned a value of 2; if 10000m 3 Total pollution volume ≤ 50,000 m³ 3 If the total pollution volume is ≥50000 m³, it is classified as severe and assigned a value of 3; 3 It is determined to be extreme and assigned the value 4.

[0011] As a further aspect of the present invention, the method for determining the corresponding pollution level is as follows: The comprehensive score is calculated using the formula: Comprehensive Score = Number of Pollutants Exceeding Standards × First Weight + Maximum Exceedance Multiple × Second Weight + Comprehensive Exceedance Multiple × Third Weight + Pollution Volume × Fourth Weight. The sum of the first weight, second weight, third weight, and fourth weight is 1. Matching the comprehensive score with the pollution level: If the comprehensive score is ≤1.5, it is judged as low risk; if 1.5 < comprehensive score ≤2.5, it is judged as medium risk; if 2.5 < comprehensive score ≤3.5, it is judged as high risk; if the comprehensive score is >3.5, it is judged as extremely high risk.

[0012] As a further aspect of the present invention, the method for generating pollution control information is as follows: For low-risk sites: implement routine supervision, include them in the urban contaminated site database and update site use and monitoring data regularly; if monitoring finds that pollutant concentrations are close to the standard limits, take simple measures such as covering the topsoil with green vegetation. For medium-risk sites: Whether to remediate depends on development needs; undeveloped sites are included in the urban contaminated site database for regular monitoring; if monitoring detects an increase in pollutant concentration, emergency control measures are immediately initiated and the risk level is reassessed. For high-risk sites: restrict the use of the site, set a repair period, and set up temporary fencing during the repair period; For extremely high-risk sites: Initiate emergency control measures to cut off contamination exposure pathways, prevent the spread of pollutants, and prioritize remediation.

[0013] As a further aspect of the present invention, when calculating the overall exceedance multiple, the average exceedance multiple of a single item is the arithmetic mean of the exceedance multiples of the single items at the sampling points for that pollutant. The geometric formulas include the area formula of a rectangle and the area formula of a trapezoid. When the pollution depth is uneven, the cross-sectional pollution area is calculated using the integral method or the segmentation method.

[0014] As a further aspect of the present invention, multi-source data is processed and cleaned. Cleaning includes standardization, data quality verification, and anomaly handling. Standardization includes format unification, unit unification, and completion of traceability information. Data quality verification and anomaly handling include missing value processing, outlier identification and processing, and removal of duplicate and invalid values.

[0015] As a further aspect of the present invention, the missing value handling is specifically as follows: when the missing rate of key data is ≤5%, spatial interpolation or the mean of adjacent samples is used to supplement it; when the missing rate is >5%, resampling is performed; when non-key data is missing, it is marked as missing; the outlier identification and handling are specifically as follows: outliers are identified using box plot method, and after checking the original records, data errors are corrected, and real outliers are retained and marked.

[0016] This invention provides a method for the graded management and control of urban industrial contaminated sites based on multi-source fusion data. Compared with existing technologies, it has the following advantages: This invention integrates four core indicators—the number of pollutants exceeding standards, the maximum multiple of exceeding standards, the comprehensive multiple of exceeding standards, and the volume of pollution—to characterize pollution features from four dimensions: pollution complexity, extreme concentration, overall intensity, and spatial scale. Compared with existing single-indicator assessment methods, it can cover complex scenarios such as multi-pollutant compound pollution and low-concentration large-scale pollution, reducing the misjudgment rate of risk levels. It clarifies the quantitative assignment standards for the four indicators and sets weight ranges to reduce subjective experience interference and meet the needs of large-scale management. Differentiated solutions are formulated for the four risk levels of low, medium, high, and extremely high. Compared with the existing one-size-fits-all control, it can reduce the cost of remediation of low-risk sites. Attached Figure Description

[0017] Figure 1 This is a diagram illustrating the steps and methods of the present invention. Detailed Implementation

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

[0019] Please see Figure 1 This application provides a method for graded management of urban industrial contaminated sites based on multi-source fusion data, which specifically includes the following steps: Step 1: Collect multi-source data corresponding to urban industrial contaminated sites, including soil pollution data and groundwater pollution data. Specific soil pollution data includes pollutant concentration, sampling depth, and sampling location coordinates of soil samples. Groundwater pollution data includes pollutant concentration, pH value, water level depth, flow direction, and aquifer type of groundwater samples. Simultaneously, the obtained multi-source data is processed and cleaned. Cleaning specifically involves standardizing the data, verifying data quality, and handling anomalies. Data standardization includes format unification, unit unification, and completion of source information. Data quality verification and anomaly handling include handling missing values, identifying and processing outliers, and removing duplicate and invalid values. The cleaned multi-source data is then fused, specifically by integrating the data spatially, by attribute, and by dimension, to extract pollution indicators corresponding to the contaminated sites. These indicators include the number of pollutants exceeding standards, the maximum exceeding multiple, the comprehensive exceeding multiple, and the pollution volume. Step 2: Assign values ​​to the quantity and maximum exceedance multiple of pollutants in the pollution indicators, and the specific processing method is as follows: For the number of pollutants exceeding the standard, the number of pollutant types whose concentrations exceed the risk screening value for the application site type is counted, and corresponding values ​​are assigned according to their corresponding quantities. Specifically, if there is only one type, it is indicated as mild and assigned a value of 1; if there are no more than three types, it is indicated as moderate and assigned a value of 2; if there are more than three types, it is indicated as severe and assigned a value of 3. For the maximum exceedance multiple, all pollutants exceeding the standard are collected, and the exceedance multiple of a single pollutant is calculated according to the formula (measured concentration - evaluation standard value) / evaluation standard value. The maximum value is obtained and compared with the evaluation standard value. If the measured concentration ≤ 2 × evaluation standard value, it is indicated as mild and assigned a value of 1. If 2 × evaluation standard value < measured concentration ≤ 4 × evaluation standard value, it is indicated as moderate and assigned a value of 2. If 4 × evaluation standard value < measured concentration ≤ 11 × evaluation standard value, it is indicated as severe and assigned a value of 3. If the measured concentration > 11 × evaluation standard value, it is indicated as extreme and assigned a value of 4.

[0020] For example, if the measured benzene concentration at a site is 12 mg / kg and the benzene screening value for the contaminated site is 4 mg / kg, then the benzene exceedance multiple is (12-4) / 4 = 2. If the measured mercury concentration is 6 mg / kg and the mercury screening value is 1.8 mg / kg, then the mercury exceedance multiple is (6-1.8) / 1.8 ≈ 2.33. Further screening yields a maximum exceedance multiple of 2.33. Comparing this with the screening value, we can assign a value of 1.

[0021] Step 3: Assign values ​​to the comprehensive exceedance multiple and pollution volume of the pollution indicators, and the specific processing method is as follows: To determine the overall exceedance multiple, we obtain the pollutants exceeding the standard at the contaminated site and the corresponding number of sampling points. We also obtain the number of exceeding sampling points and their exceedance multiples. Here, the exceedance multiple represents the average exceedance multiple across all exceeding sampling points. First, we calculate the exceedance multiple for each individual exceeding sampling point. Then, we calculate the overall average multiple and the corresponding percentage of exceeding sampling points. The percentage of exceeding sampling points = number of exceeding sampling points for that pollutant / total number of sampling points for that pollutant. Finally, we apply the formula: Overall Exceedance Multiple = (The multiple of a pollutant exceeding the standard × the percentage of locations exceeding the standard for that pollutant) / the number of types of pollutants exceeding the standard. The overall multiple of exceeding the standard is compared with the threshold value, and the specific value of the threshold value is set by the operator. If the overall multiple of exceeding the standard is ≤ 0.5 times the threshold value, it is indicated as mild and assigned a value of 1. If the overall multiple of exceeding the standard is 0.5-1 times the threshold value, it is indicated as moderate and assigned a value of 2. If the overall multiple of exceeding the standard is 1-3 times the threshold value, it is indicated as severe and assigned a value of 3. If the overall multiple of exceeding the standard is ≥ 3 times the threshold value, it is indicated as extreme and assigned a value of 4. To determine the pollution volume, the direction of pollution diffusion is obtained, and vertical cross-sections are laid out at intervals L. The number of vertical cross-sections, k, is also obtained. Then, a depth-concentration curve is plotted based on the borehole data to determine the pollution depth range of each cross-section. The pollution area S of a single cross-section is calculated using a geometric formula, which can be applied to rectangles or trapezoids. For uneven pollution depths, the area is calculated using integration or segmentation methods. The total pollution volume is calculated using the formula: (S1×L+S2×L+…+S) k-1 The total pollution volume is calculated using the formula (k-1, S1, S2, S...). k-1 This represents the contaminated area of ​​a corresponding cross-section. For example, in a narrow contaminated zone formed by a pipeline leak, five cross-sections are laid out with a spacing of 10m. The contaminated area of ​​each cross-section is: S1 = 30m². 2 S2=45m 2 S3=50m 2 S4=35m 2 S5=20m 2 The total volume = ((30+45) / 2+(45+50) / 2+(50+35) / 2+(35+20) / 2)×10 = (37.5+47.5+42.5+27.5)×10 = 1550m³ 3 ; The total pollution volume is then assigned a value. If the total pollution volume is ≤1000 m³, the value is adjusted accordingly. 3 If the total pollution volume is 1000-10000 m³, it is considered mild and assigned a value of 1. 3If the total pollution volume is 10,000-50,000 m³, it is considered moderate and assigned a value of 2. 3 If the total pollution volume is ≥50000 m³, it indicates severe pollution and is assigned a value of 3. 3 If , it indicates extreme, and is assigned the value 4.

[0022] Step 4: The obtained pollution indicators, including the number of pollutants exceeding the standard, the maximum exceeding multiple, the comprehensive exceeding multiple, and the pollution volume, are weighted and summed. The comprehensive score is calculated according to the formula: Comprehensive Score = Number of Pollutants Exceeding the Standard × First Weight + Maximum Exceeding Multiple × Second Weight + Comprehensive Exceeding Multiple × Third Weight + Pollution Volume × Fourth Weight. The specific values ​​of the weights are set by the operator. At the same time, the obtained comprehensive score is matched with the corresponding evaluation criteria to generate the corresponding pollution level. If the overall score is ≤1.5, it indicates low risk; if 1.5 < overall score ≤2.5, it indicates medium risk; if 2.5 < overall score ≤3.5, it indicates high risk; and if the overall score >3.5, it indicates extremely high risk.

[0023] Step 5: Different management methods are implemented based on the pollution level. For low-risk sites, routine supervision is carried out, and the sites are included in the urban contaminated site database. The site use and monitoring data are updated regularly. If the monitoring finds that the pollutant concentration is close to the standard limit, simple measures can be taken, such as covering the top soil with green vegetation to reduce dust exposure, without the need to start formal remediation. For medium-risk sites, the decision to remediate depends on the site's development needs. Undeveloped sites may not require remediation for the time being, but they will be included in the city's contaminated site database and tracked regularly. If monitoring reveals an increase in pollutant concentration, emergency control measures must be initiated immediately, and the risk level must be reassessed. For high-risk sites, the use of the site is restricted, and a corresponding repair period is set. Temporary fencing is also set up during the repair period. In response to the extremely high risk, emergency control measures were initiated to immediately cut off the pollution exposure routes, prevent the spread of pollutants, and prioritize recovery efforts.

[0024] Some of the data in the above formulas are numerical calculations with dimensions removed, and the contents not described in detail in this specification are all prior art known to those skilled in the art.

[0025] The above embodiments are only used to illustrate the technical methods of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical methods of the present invention without departing from the spirit and scope of the technical methods of the present invention.

Claims

1. A method for hierarchical management and control of urban industrial contaminated sites based on multi-source fusion data, characterized in that, The method specifically comprises the following steps: Collecting multi-source data of urban industrial contaminated sites, and performing arrangement and cleaning, fusing the cleaned multi-source data, and extracting four pollution indexes of pollution exceeding number, maximum exceeding multiple, comprehensive exceeding multiple, and pollution square amount; Assigning values to the pollution exceeding number and the maximum exceeding multiple, respectively, assigning values based on the pollution exceeding number in different grades, and assigning values based on the multiple relationship between the maximum exceeding multiple and the corresponding evaluation standard value; Assigning values to the comprehensive exceeding multiple and the pollution square amount, respectively, calculating the comprehensive exceeding multiple of all exceeding pollutants and assigning values, and calculating the total pollution square amount and assigning values; Weighted summing of the values of different pollution indexes to obtain a comprehensive score, and matching with the judgment standard to determine the corresponding pollution grade; According to the low risk, medium risk, high risk and extremely high risk in the pollution grade, different management and control treatments are performed to generate pollution management and control information.

2. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 1, characterized in that, The multi-source data includes soil pollution data and groundwater pollution data, the soil pollution data includes the concentration of pollutants of soil samples, sampling depth and sampling position coordinates; The groundwater pollution data includes the concentration of pollutants of groundwater samples, pH value, water level depth, water flow direction and aquifer type.

3. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 1, characterized in that, The assignment method of the pollution exceeding number is: counting the number of pollutant types whose concentration exceeds the risk screening value of the corresponding land type, if the number is 1, it is determined as mild and assigned a value of 1; if the number is 2-3, it is determined as moderate and assigned a value of 2; if the number is ≥4, it is determined as severe and assigned a value of 3; The assignment method of the maximum exceeding multiple is: calculating the single exceeding multiple of each exceeding pollutant according to the formula (measured concentration - evaluation standard value) / evaluation standard value and taking the maximum value, if the measured concentration corresponding to the maximum value is ≤2×evaluation standard value, it is determined as mild and assigned a value of 1; if 2×evaluation standard value < measured concentration ≤4×evaluation standard value, it is determined as moderate and assigned a value of 2; if 4×evaluation standard value < measured concentration ≤11×evaluation standard value, it is determined as severe and assigned a value of 3; if measured concentration > 11×evaluation standard value, it is determined as extremely severe and assigned a value of 4.

4. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 1, characterized in that, The assignment method of the comprehensive exceeding multiple is: calculating the single exceeding multiple mean value of each exceeding pollutant and the exceeding point position ratio (exceeding point position ratio = the number of exceeding sampling points of the pollutant / the total sampling points of the pollutant), and calculating the comprehensive exceeding multiple according to the formula comprehensive exceeding multiple = Σ (the multiple mean value of a certain exceeding pollutant × the exceeding point position ratio of the pollutant) / the number of exceeding pollutant types; Setting an exceeding threshold, if the comprehensive exceeding multiple is ≤0.5×exceeding threshold, it is determined as mild and assigned a value of 1; if 0.5×exceeding threshold < comprehensive exceeding multiple ≤1×exceeding threshold, it is determined as moderate and assigned a value of 2; if 1×exceeding threshold < comprehensive exceeding multiple ≤3×exceeding threshold, it is determined as severe and assigned a value of 3; if comprehensive exceeding multiple ≥3×exceeding threshold, it is determined as extremely severe and assigned a value of 4.

5. The urban industrial contaminated site hierarchical management and control method based on multi-source fusion data according to claim 1, characterized in that, The assignment method of the pollution square amount is: arranging k vertical sections with an interval L along the pollution diffusion direction, determining the pollution depth range of each section according to the depth-concentration curve of the drilling data, and calculating the pollution area S of each section by using the geometric formula combined with the horizontal width of the section. The total pollution square amount is calculated according to the formula total pollution square amount = (S1×L+S2×L+…+Sk-1×L+Sk×L) wherein k-1 represents the number of corresponding sections. k-1 wherein k-1 represents the number of corresponding sections. If the total pollution square amount ≤ 1000 m 3 , it is determined to be light and assigned a value of 1; if 1000 m 3 < total pollution square amount ≤ 10000 m 3 , it is determined to be moderate and assigned a value of 2; if 10000 m 3 < total pollution square amount ≤ 50000 m 3 , it is determined to be severe and assigned a value of 3; if the total pollution square amount ≥ 50000 m 3 , it is determined to be extreme and assigned a value of 4.

6. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 1, characterized in that, The way to determine the corresponding pollution grade is: The comprehensive score is calculated according to the formula: comprehensive score = (number of pollutants exceeding the standard) * first weight + (maximum exceeding multiple) * second weight + (comprehensive exceeding multiple) * third weight + (pollution amount) * fourth weight, wherein the sum of the first weight, the second weight, the third weight and the fourth weight is 1; The comprehensive score is matched with the pollution level: if the comprehensive score is less than or equal to 1.5, it is determined as low risk; if 1.5 < comprehensive score < 2.5, it is determined as medium risk; if 2.5 < comprehensive score < 3.5, it is determined as high risk; and if the comprehensive score is greater than 3.5, it is determined as extremely high risk.

7. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 1, characterized in that, The way of generating pollution control information is: For low-risk sites, routine supervision is implemented, and the sites are included in the urban pollution site database and the site use and monitoring data are updated regularly. If the monitoring finds that the concentration of pollutants is close to the standard limit, surface soil is covered with green vegetation; For medium-risk sites, it is decided whether to repair according to the development needs. The sites that are not developed are included in the urban pollution site database for regular tracking. If the monitoring finds that the concentration of pollutants is rising, emergency control is started immediately and the risk level is re-evaluated; For high-risk sites, the site use is limited, and a repair period is set. Temporary fences are set up during the repair period; For extremely high-risk sites, emergency control is started to cut off the exposure path of pollution and prevent the spread of pollutants, and repair is started preferentially.

8. The urban industrial contaminated site grading management method based on multi-source fusion data according to claim 4, characterized in that, When calculating the comprehensive exceeding multiple, the single exceeding multiple mean value is the arithmetic mean of the single exceeding multiple of the sampling points of the pollutant. The geometric formula includes the rectangular area formula and the trapezoidal area formula. When the pollution depth is uneven, the integral method or the segmentation method is used to calculate the cross-sectional pollution area. 9.The method of claim 1, wherein, The multi-source data are arranged and cleaned. The cleaning includes standardized arrangement, data quality verification and abnormality processing. The standardized arrangement includes format unification, unit unification and traceability information completion. The data quality verification and abnormality processing include missing value processing, abnormal value identification and processing, and repeated value and invalid value elimination.

10. The method according to claim 1, wherein, The missing value processing is specifically: when the missing rate of key data is less than or equal to 5%, the spatial interpolation method or the adjacent sample mean method is used for supplement; when the missing rate is greater than 5%, resampling is performed; when the missing rate of non-key data is less than or equal to 5%, the missing value is marked; the abnormal value identification and processing is specifically: the box plot method is used to identify abnormal values, and after checking the original records, the data errors are corrected, and the true abnormal values are retained and marked.

Images