Method for determining pollution range and pollution control priority of target enterprise

Abstract

The invention relates to a method for determining the pollution range and pollution control priority of a target enterprise, and the method comprises the four steps of database establishment, regional content distribution diagram drawing, pollution region determination and pollution mode screening, integrates spatial interpolation, a random walk model, a skeleton algorithm, an IOU algorithm and other technical means, achieves the multi-level and multi-method linkage, and achieves the purpose of determining the pollution range and the pollution control priority of the target enterprise. And the spatial distribution information and the pollution mode of soil heavy metal pollution caused by enterprises are clear in real time, quickly, efficiently and clearly. Scientific algorithms and actual requirements are considered, and later soil pollution prevention and agricultural product quality safety guarantee are facilitated.

Description

Technical Field

[0001] The invention belongs to the field of agricultural environmental science, and specifically relates to a method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise. Background Art

[0002] Heavy metal pollution in soil will directly affect the quality and safety of agricultural products and human health (Qian Junjun. Research on the problem and control measures of heavy metal pollution in soil [J]. Cleaning the World. 2024, 40(05): 78-80). Comprehensive, real-time and accurate soil environmental monitoring has important practical significance for the quality of agricultural products and soil. Accurately obtaining the scope and mode of heavy metal pollution caused by target enterprises is the basis for carrying out soil pollution prevention and control work. At present, the methods for obtaining enterprise heavy metal monitoring data are mostly fixed-point sampling and point source observation, and the research is mostly carried out on the type and degree of heavy metal pollution (Li Chao, Xu Huanhuan, Gao Fangying. Investigation and research on the status of heavy metal pollution in soil of renewable resource recycling enterprises [J]. Environmental Engineering 2023, 41(S2): 1190-1195; Han Zhouyang. Investigation, evaluation and analysis of heavy metal industrial pollution sources in Suqian City [J]. Green Technology, 2017(24): 23-25, 29.).

[0003] The main problems with the above-mentioned technical methods are: (1) soil monitoring technology is single, data acquisition is difficult, and time and financial investment are high; (2) there is little research on the scope and mode of soil heavy metal pollution caused by enterprises; (3) the implementation of prevention and control measures for soil heavy metal pollution caused by enterprises is blind. Summary of the invention

[0004] The present invention uses four major steps, including database establishment, regional content distribution map drawing, pollution area clarification, and pollution mode screening, and integrates technical means such as spatial interpolation, random walk model, skeleton algorithm, and IOU algorithm, to link multiple levels and methods, and to clarify the spatial distribution information and pollution mode of soil heavy metal pollution caused by enterprises in real time, quickly, and efficiently. The present invention takes into account both scientific algorithms and actual needs, and is conducive to the prevention and control of soil pollution and the quality and safety of agricultural products in the later stage.

[0005] In order to solve the above technical problems, the present invention discloses a method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise, comprising the following steps:

[0006] (1) Establishing a database

[0007] The database includes target enterprise information, environmental information, other enterprise information and soil monitoring information of the research area; wherein the target enterprise refers to one or more polluting enterprises of known heavy metal pollution types selected in the research area, and the polluting enterprises are set as the research subjects, and the target enterprise information includes one or more of enterprise type, geographical location, business information, polluting heavy metals, and historical pollution events; obtain environmental factor information on the spread of soil heavy metal pollution caused by the target enterprise, including water source information, wind direction information, and traffic information; the other enterprises refer to all enterprises in the research area except the target enterprise;

[0008] (2) Drawing of regional content distribution map

[0009] (2.1) Obtaining soil heavy metal content

[0010] (2.1.1) Determine the number of points

[0011] (2.1.2) Divide the study area into grids according to the number of points. The center of the grid is the monitoring point. Obtain the longitude and latitude information of the monitoring point and upload it to the database;

[0012] When the monitoring point exceeds the scope of the study area and the area of ​​the study area occupies more than 25% of the area of ​​the grid, the grid should select the point closest to the center point of the grid in a straight line and located in the study area as the monitoring point;

[0013] (2.2) Draw a regional content distribution map

[0014] (2.2.1) Obtaining the soil at the monitoring point, determining the content of each heavy metal in the soil, and uploading the metal detection data to the database;

[0015] (2.2.2) Combined with the information of the research area, the spatial interpolation method is used to obtain the regional content distribution map of the polluting heavy metals. Each polluting heavy metal corresponds to a regional content distribution map;

[0016] (3) The pollution area is clear

[0017] (3.1) Determination of the main pollution area

[0018] (3.1.1) Based on the regional content distribution map of heavy metal pollution, obtain the content of heavy metal pollution in each grid in the study area; the grid is the minimum area unit for spatial interpolation;

[0019] The number of grids per square meter in the study area should be greater than 100;

[0020] (3.1.2) Using a random walk model, setting the walk model parameters to obtain a closed area; the closed area containing the target enterprise is the main pollution area; the parameters include kernel function, step size, switch parameter, boundary control function, stop condition, and path; the path parameter is set as a closed loop;

[0021] The closed area is based on the heavy metal content of each grid, the difference in heavy metal content in adjacent grids is calculated, the grid is probabilized using a kernel function, a grid is randomly selected as the starting point, and a random walk (8 directions) is performed according to the probability of adjacent grids. The step size is set to k, and a judgment is made after walking k steps, and the walk continues in the direction with the maximum probability until it stops at the starting point. The area contained in the irregular figure finally formed is the closed area, and the closed area is one or more;

[0022] (3.2) Determination of contaminated radiation areas

[0023] (3.2.1) Based on the target enterprise information in the database, identify all product categories of the target enterprise;

[0024] (3.2.2) For all closed areas, filter out the closed areas that do not contain the target enterprise, and obtain other enterprise information of all other enterprises in the corresponding closed areas based on the database;

[0025] (3.2.3) Based on the information obtained from other enterprises, determine the types of raw materials of all other enterprises;

[0026] (3.2.4) Associating and matching the target enterprise's product categories with the raw material categories of other enterprises. If the target enterprise's product categories are consistent with the raw material categories of other enterprises, then the closed area where the other enterprises are located is the contaminated radiation area; the contaminated radiation area may be 0, 1, or more;

[0027] (3.3) Determination of contaminated areas

[0028] (3.3.1) If there is no contaminated radiation area in the study area, the main contaminated area is the contaminated area;

[0029] (3.3.2) If there is a contaminated radiation area in the study area, the contaminated area is the sum of the main contaminated area and the contaminated radiation area;

[0030] (4) Pollution mode screening: Pollution mode screening is only carried out in the main pollution area

[0031] (4.1) Calculation of contaminated main ridge

[0032] Based on the content of polluted heavy metals in each grid of the main pollution area and the boundary of the closed area, the main pollution ridge is extracted by using a skeleton algorithm; the skeleton algorithm uses a thinning algorithm, that is, gradually eroding from the boundary of the main pollution area into the area until a main pollution ridge is retained, thereby ensuring the topological structure and geometric structure of the main pollution area; there is only one main pollution ridge in the main pollution area;

[0033] (4.2) Surface element calculation

[0034] (4.2.1) Using the professional software, based on the area of ​​the main pollution area, draw a circumscribed ellipse of the main pollution area and obtain the short radius of the circumscribed ellipse; the circumscribed ellipse is the smallest unit that includes all the main pollution areas;

[0035] (4.2.1) Using 5-95% of the short radius of the circumscribed ellipse as the width unit, widen the pollution main ridge, convert the line element into a surface element, and obtain the pollution main ridge surface element;

[0036] (4.2.3) Read the environmental information in the database, obtain the environmental factor information, and convert all environmental factors into line elements; the line elements include water source line elements, traffic line elements and wind direction line elements. The conversion method of the line elements is to extract the center lines of rivers, irrigation canals, railways and highways as water source or traffic line elements; the wind direction line elements are drawn along the wind direction with the boundary of the research area as the starting point and the end point and the target enterprise as the passing point;

[0037] (4.2.4) Using 5-95% of the short radius of the circumscribed ellipse as the width unit, widen the line elements of all environmental factors, convert them from line elements to surface elements, and obtain the surface elements of environmental factors;

[0038] (4.3) Pollution prevention and control priority determination

[0039] (4.3.1) Calculate the area of ​​the main pollution ridge surface element and the area of ​​all environmental factor surface elements respectively;

[0040] (4.3.2) The IOU is calculated by comparing the area of ​​the main pollution ridge element with the area of ​​all environmental factor elements, and the IOU results are sorted. The larger the IOU, the greater the impact of the corresponding environmental factor on the pollution of the target enterprise. The specific IOU calculation formula is as follows:

[0041]

[0042] Where R is the area of ​​the main pollution ridge element, E is the area of ​​the environmental factor surface element; the IOU range is 0-1;

[0043] (4.4) Selection of prevention and control measures

[0044] Based on the IOU calculated for each environmental factor, the IOU is sorted from large to small. During pollution prevention and control, the prevention and control measures corresponding to the environmental factors can be selected in sequence to ensure soil quality safety.

[0045] Furthermore, the types of enterprises include animal husbandry, electronic and electrical industry, mechanical and electrical industry, food processing industry, metallurgical and mining industry, construction and building materials industry, petrochemical industry, etc.;

[0046] Furthermore, the polluting heavy metals include one or more of Cd, Hg, As, Pb, Cr, Cu, Zn, Ni, Mg, Fe, Mn, Ca, Si, and Se;

[0047] Furthermore, the business information includes one or more of the following: enterprise name, business scope, raw material type, product type, output volume, enterprise scale, and business years;

[0048] Further, environmental information of the study area is obtained, wherein the environmental information includes one or more of the study area area, soil type, topography, landforms, landforms, planting systems, crop types, and soil historical monitoring information;

[0049] Furthermore, the pollution mode screening is only carried out on the main pollution area;

[0050] Furthermore, the other enterprise information is consistent with the target enterprise information, including enterprise type, geographical location, operating information, and historical pollution events;

[0051] Furthermore, the heavy metal content of the soil is measured by using an X-fluorescence detector;

[0052] Furthermore, the spatial interpolation methods include Kriging interpolation, inverse distance weighted average interpolation (IDW), trend surface method, spline function method, etc.; the professional software includes ArcGIS, R, python, matlab, CAD, etc.;

[0053] Furthermore, based on the historical soil monitoring information in the database, the minimum number of monitoring points in the study area is calculated using the following formula:

[0054]

[0055] Wherein, N is the minimum number of monitoring points, t is the passing threshold of the T test, m is the tolerance error, Cv is the coefficient of variation, Cv=σ / μ, σ and μ are the standard deviation and average value of the heavy metal pollution corresponding to all historical soil monitoring points, respectively, and the number of monitoring points is more than the minimum number of monitoring points;

[0056] (2.2) Content distribution diagram of receipt area

[0057] (2.2.1) Obtaining soil from the monitoring point, measuring the content of each heavy metal in the soil, and uploading the heavy metal detection data consistent with the type of heavy metal pollution of the target enterprise to the database;

[0058] (2.2.2) Combined with the information of the research area, the spatial interpolation method is used to obtain the regional content distribution map of the polluting heavy metals. Each polluting heavy metal corresponds to a regional content distribution map;

[0059] The method of the present invention for clarifying the pollution scope and pollution prevention and control priority of a target enterprise has the following advantages:

[0060] 1. The present invention adopts the method of global point distribution combined with X-fluorescence rapid detection, which can quickly complete the detection work on the basis of meeting the data volume requirements, provide a data basis for the subsequent determination of the pollution scope of the target enterprise, and ensure the accuracy of the pollution area division results;

[0061] 2. The random walk combined with the skeleton algorithm adopted in the present invention realizes the linkage innovation of technology, ensures the accuracy of the pollution area division and the representativeness of the pollution main ridge, and provides technical support for the selection of subsequent prevention and control measures;

[0062] 3. The present invention adopts the IOU algorithm and combines regional environmental factors to analyze the pollution mode, and then selects appropriate prevention and control measures to target the farmland pollution caused by target enterprises, reducing about 20-35% of the financial and time cost losses caused by the blind implementation of measures.

[0063] The above description is only an overview of the technical solution of the present application. In order to more clearly understand the technical means of the present application so that it can be implemented in accordance with the contents of the specification, and to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the following is a detailed description of the preferred embodiments of the present application in conjunction with the accompanying drawings as follows.

[0064] Based on the detailed description of the specific embodiments of the present application in combination with the accompanying drawings below, those skilled in the art will become more aware of the above and other objects, advantages and features of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0065] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can also be obtained based on these drawings without creative work. In all drawings, similar elements or parts are generally identified by similar reference numerals. In the drawings, each element or part is not necessarily drawn according to the actual scale.

[0066] Figure 1 It is a clear method and technical flow chart for the pollution scope and pollution prevention and control priorities of the target enterprises;

[0067] Figure 2 This is the distribution map of monitoring points in T County;

[0068] Figure 3 It is the distribution map of the main pollution ridge and closed area;

[0069] Figure 4 It is the surface element distribution map; Specific embodiments

[0070] The present invention is further described in detail below through examples so that those skilled in the art can implement the invention with reference to the description.

[0071] It should be understood that the terms such as “having”, “including” and “comprising” used herein do not exclude the existence or addition of one or more other elements or combinations thereof.

[0072] Example 1

[0073] 1. Database establishment

[0074] (1) Target enterprise information

[0075] Taking enterprise A in T County as the research object, through consulting literature and reports, it is determined that the main polluting heavy metal of target enterprise A is Cu, and two pollution incidents have occurred in the past two years. The main finished products are steel, copper materials, and building supplies;

[0076] (2) Environmental information

[0077] By consulting yearbook reports and documents, it was determined that the environmental factors affecting pollution in T County include water sources and wind direction, and that the river runs across the county and the wind direction is northwesterly. In addition, it was determined that T County is a hilly area with little rainfall and mainly grows wheat. More than one thousand pieces of historical soil monitoring information were obtained.

[0078] (3) Other company information

[0079] Through online searches and field visits, it was determined that in addition to the target enterprise A, there were 267 other enterprises (in operation) in T County, with a total of two industrial parks, and most of the other enterprises were scattered; and information on other enterprises was collected, of which 78 other enterprises had incomplete information;

[0080] 2. Regional content distribution map drawing

[0081] (1) Using the minimum monitoring point calculation formula, it is calculated that the minimum number of monitoring points in County T is 118;

[0082] (2) With a grid size of 5*5 km, the study area was divided into 141 grids and 129 monitoring points were selected, which is more than the minimum number of monitoring points calculated. Figure 2 , which includes 7 adjusted monitoring points, uses GPS to obtain the latitude and longitude information of the monitoring points and uploads them to the database;

[0083] (3) Using an X-ray fluorescence detector, collect the surface soil at the monitoring point in situ, test the Cu content, and upload the test results to the database;

[0084] (4) The terrain of T County is relatively flat, so the Kriging interpolation method was selected and ArcGIS software was used to perform spatial interpolation on the Cu content obtained in the previous step to obtain the Cu regional content distribution map;

[0085] 3. Polluted areas are clear

[0086] (1) Based on the Cu regional content distribution map, the random walk model was selected, and the parameters were debugged using Python software to obtain three closed areas. Among them, closed area 1 contained the target enterprise A, which was the main pollution area;

[0087] (2) Check the database information and confirm that there are 67 other enterprises in closed area 2 and closed area 3. The other enterprises with incomplete enterprise information have completed the enterprise information;

[0088] (3) The product types of target enterprise A were correlated and matched with the raw material types of 67 enterprises. The results showed that there were other enterprises in both areas with raw material types that were consistent with the product types of target enterprise A. Therefore, closed area 2 and closed area 3 were both pollution radiation areas of target enterprise A.

[0089] (4) The pollution area caused by target enterprise A in the research area includes 1 main pollution area and 2 pollution radiation areas, a total of 3 closed areas, see Figure 3 ;

[0090] 4. Pollution screening

[0091] (1) Using Python software, call the skeleton function of skimage to obtain one main contamination ridge in the main contamination area;

[0092] (2) Using CAD software, draw a circumscribed ellipse of the main pollution area. The short radius of the circumscribed ellipse is about 7 km. 20% of the radius is selected as the width unit, i.e. 1.4 km, to widen the main pollution ridge;

[0093] (3) Consult the database to obtain the environmental information of County T, extract the water source line elements and wind direction line elements, and widen the water source line elements and wind direction line elements with a width unit of 1.4 km to obtain the water source surface elements and wind direction surface elements. Figure 4 ;

[0094] (4) Calculate the surface area of ​​the main pollution ridge, the surface area of ​​the water source, and the surface area of ​​the wind direction respectively, and use the IOU algorithm to obtain the IOU water source and IOU wind direction. The results show that IOU wind direction > IOU water source;

[0095] (5) Based on the calculation results of the previous step, the prevention and control measures for target enterprise A give priority to technical means of blocking wind propagation, such as reducing the smoke emissions of target enterprise A and planting trees near target enterprise A.

[0096] Although the embodiments of the present invention have been disclosed as above, they are not limited to the applications listed in the specification and the implementation modes. They can be fully applied to various fields suitable for the present invention. For those familiar with the art, additional modifications can be easily implemented. Therefore, without departing from the general concept defined by the claims and the scope of equivalents, the present invention is not limited to the specific details and the embodiments shown and described herein.

Claims

1. A method for clarifying the pollution scope and pollution prevention and control priorities of target enterprises, characterized by: The following steps are involved: (1) Establishing a database The database includes target enterprise information, environmental information, other enterprise information and soil monitoring information of the research area; (2) Drawing of regional content distribution map (2.1) Obtaining soil heavy metal content (2.1.1) Determine the number of points (2.1.2) Divide the study area into grids according to the number of points. The center of the grid is the monitoring point. Obtain the longitude and latitude information of the monitoring point and upload it to the database; When the monitoring point exceeds the scope of the study area and the area of ​​the study area occupies more than 25% of the area of ​​the grid, the grid should select the point closest to the center point of the grid in a straight line and located in the study area as the monitoring point; (2.2) Draw a regional content distribution map (2.2.1) Obtaining the soil at the monitoring point, determining the content of each heavy metal in the soil, and uploading the metal detection data to the database; (2.2.2) Combined with the information of the research area, the spatial interpolation method is used to obtain the regional content distribution map of the polluting heavy metals. Each polluting heavy metal corresponds to a regional content distribution map; (3) The pollution area is clear (3.1) Determination of the main pollution area (3.1.1) Based on the regional content distribution map of heavy metal pollution, obtain the content of heavy metal pollution in each grid in the study area; the grid is the minimum area unit for spatial interpolation; (3.1.2) Using a random walk model, setting the walk model parameters to obtain a closed area; the closed area containing the target enterprise is the main pollution area; the parameters include kernel function, step size, switch parameter, boundary control function, stop condition, and path; the path parameter is set as a closed loop; The closed area is based on the heavy metal content of each grid, calculates the difference in heavy metal content in adjacent grids, uses kernel function to probabilize the grid, randomly selects a grid as the starting point, and randomly walks in 8 directions according to the probability of adjacent grids. The step size is set to k, and after walking k steps, a judgment is made, and the walk continues in the direction with the maximum probability until it stops at the starting point. The area contained in the irregular figure finally formed is the closed area, and the closed area is one or more; (3.2) Determination of contaminated radiation areas (3.2.1) Based on the target enterprise information in the database, identify all product categories of the target enterprise; (3.2.2) For all closed areas, filter out the closed areas that do not contain the target enterprise, and obtain other enterprise information of all other enterprises in the corresponding closed areas based on the database; (3.2.3) Based on the information obtained from other enterprises, determine the types of raw materials of all other enterprises; (3.2.4) Associating and matching the target enterprise's product categories with the raw material categories of other enterprises. If the target enterprise's product categories are consistent with the raw material categories of other enterprises, then the closed area where the other enterprises are located is the contaminated radiation area; the contaminated radiation area may be 0, 1, or more; (3.3) Determination of contaminated areas (3.3.1) If there is no contaminated radiation area in the study area, the main contaminated area is the contaminated area; (3.3.2) If there is a contaminated radiation area in the study area, the contaminated area is the sum of the main contaminated area and the contaminated radiation area; (4) Pollution screening (4.1) Calculation of contaminated main ridge Based on the content of polluted heavy metals in each grid of the main pollution area and the boundary of the closed area, the main pollution ridge is extracted by using a skeleton algorithm; the skeleton algorithm uses a thinning algorithm, that is, gradually eroding from the boundary of the main pollution area into the area until a main pollution ridge is retained, thereby ensuring the topological structure and geometric structure of the main pollution area; there is only one main pollution ridge in the main pollution area; (4.2) Area element calculation (4.2.1) Using the professional software, based on the area of ​​the main pollution area, draw a circumscribed ellipse of the main pollution area and obtain the short radius of the circumscribed ellipse; the circumscribed ellipse is the smallest unit that includes all the main pollution areas; (4.2.1) Using 5-95% of the short radius of the circumscribed ellipse as the width unit, widen the pollution main ridge, convert the line element into a surface element, and obtain the pollution main ridge surface element; (4.2.3) Read the environmental information in the database, obtain the environmental factor information, and convert all environmental factors into line elements; the line elements include water source line elements, traffic line elements and wind direction line elements. The conversion method of the line elements is to extract the center lines of rivers, irrigation canals, railways and highways as water source or traffic line elements; the wind direction line elements are drawn along the wind direction with the boundary of the research area as the starting point and the end point and the target enterprise as the passing point; (4.2.4) Using 5-95% of the short radius of the circumscribed ellipse as the width unit, widen the line elements of all environmental factors, convert them from line elements to surface elements, and obtain the environmental factor surface elements; (4.3) Pollution prevention and control priority determination (4.3.1) Calculate the area of ​​the main pollution ridge surface element and the area of ​​all environmental factor surface elements respectively; (4.3.2) The IOU is calculated by comparing the area of ​​the main pollution ridge element with the area of ​​all environmental factor elements, and the IOU results are sorted. The larger the IOU, the greater the impact of the corresponding environmental factor on the pollution of the target enterprise. The specific IOU calculation formula is as follows: Where R is the area of ​​the main pollution ridge element, E is the area of ​​the environmental factor surface element; the IOU range is 0-1; (4.4) Selection of prevention and control measures Based on the IOU calculated for each environmental factor, the IOU is sorted from large to small. During pollution prevention and control, the prevention and control measures corresponding to the environmental factors can be selected in sequence to ensure soil quality safety.

2. A method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The types of enterprises mentioned include animal husbandry, electronics and electrical industry, mechanical and electrical industry, food processing industry, metallurgy and mining industry, construction and building materials industry, and petrochemical industry.

3. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1 is characterized by: The polluting heavy metals include one or more of Cd, Hg, As, Pb, Cr, Cu, Zn, Ni, Mg, Fe, Mn, Ca, Si, and Se.

4. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The pollution mode screening is only carried out on the main pollution area.

5. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The environmental information includes one or more of the study area area, soil type, topography, landforms and features, planting system, crop type, and soil historical monitoring information; the business information includes one or more of the company name, business scope, raw material types, product types, output, company scale, and years of operation.

6. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: Target enterprises refer to one or more polluting enterprises with known heavy metal pollution types selected in the research area, and the polluting enterprises are set as the research subjects. The target enterprise information includes one or more of enterprise type, geographical location, business information, polluting heavy metals, and historical pollution events; Obtain information on environmental factors that may affect the spread of heavy metal pollution in soil caused by the target enterprise, including water source information, wind direction information, and traffic information; The other enterprises refer to all enterprises in the study area except the target enterprise; the information of the other enterprises is consistent with the information of the target enterprise, including enterprise type, geographical location, operating information, and historical pollution events.

7. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The heavy metal content in the soil was determined by using an X-ray fluorescence detector.

8. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The spatial interpolation methods include Kriging interpolation method, inverse distance weighted average interpolation method, trend surface method and spline function method.

9. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The spatial interpolation method is calculated using professional software, which includes ArcGIS, R, Python, Matlab, and CAD.

10. The method for clarifying the pollution scope and pollution prevention and control priority of a target enterprise according to claim 1, characterized in that: The number of grids per square meter in the study area should be greater than 100.

Images