Method and system for constructing industrial site soil pollutant emission inventory

Abstract

The application provides a kind of industrial site soil pollutant discharge list construction method and system, comprising: according to site information, the characteristic pollutant of site is determined;Identify the main process source of pollution of site, main discharge path, pollution area, establish the soil pollution account of the site;In time scale, spatial scale, the pollution discharge condition recorded in soil pollution account is determined analysis scale;In corresponding time scale, spatial scale, the emission of characteristic pollutant is calculated, and a calculation table is obtained;Draw time discharge list, space discharge list;Carry out data analysis, obtain the soil pollutant discharge list of site, soil pollutant discharge list includes: the characteristic pollutant, the main discharge path, the pollution area, the pollution contribution under multiple dimensions;Multiple dimensions include: time scale, spatial scale, main process source of pollution, main discharge path.The application can completely reflect the whole production process.

Description

Technical Field

[0001] This invention relates to the field of soil pollutant emission inventory, and more particularly to a method and system for constructing an industrial site soil pollutant emission inventory. Background Technology

[0002] Historically, extensive industrial development has resulted in significant site pollution, necessitating source tracing and management of both site and groundwater contamination. This necessitates research into soil pollution emissions from these sites. However, obtaining actual pollution levels and data for historical sites is often challenging, particularly for historical industrial sites where substantial data gaps exist (e.g., enterprise wastewater discharge data, industrial activity information, etc.).

[0003] Previous investigations of historical sites typically relied on on-site exploration, data collection, and multiple rounds of on-site sampling (preliminary site investigation, detailed site investigation, etc.), followed by testing of pollutants in the sampled soil. This approach neglected process information, production information, site process distribution maps, and the intensity and location of emission sources. It failed to establish a quantitative relationship between soil pollution emissions and multi-source information from enterprises and actual site survey data. Furthermore, the results obtained only reflected the pollution situation at the sampling point. Pollutants undergo migration and transformation over a certain period, resulting in a single-point pollution situation. This approach cannot reflect changes in soil pollution emissions over time, space, and process scales, and therefore cannot provide a complete picture of the entire production process. Summary of the Invention

[0004] This invention provides a method and system for constructing an inventory of soil pollutants at industrial sites, in order to solve the problem that existing survey methods ignore a lot of information and cannot fully reflect the entire production process.

[0005] To solve the above-mentioned technical problems, the present invention is achieved through the following technical solution:

[0006] According to a first aspect of the present invention, a method for constructing an industrial site soil pollutant emission inventory is provided, comprising:

[0007] Identify characteristic pollutants: Based on the site information, identify the characteristic pollutants of the site;

[0008] Establish a soil pollution account: Analyze the site information to identify the main technological sources, main emission pathways, and polluted areas of the site; establish a soil pollution account for the site based on the characteristic pollutants, the main technological sources, the main emission pathways, and the polluted areas;

[0009] Obtain the calculation table: Determine the analysis scale for the pollution discharge recorded in the soil pollution account under the time and spatial scales; calculate the discharge of the characteristic pollutants under the corresponding time and spatial scales to obtain the calculation table;

[0010] Emissions Inventory: Based on the polluted area and the analytical scale, generate temporal and spatial emission inventories;

[0011] Obtaining a soil pollutant emission inventory: By analyzing the calculated table, the temporal emission inventory, and the spatial emission inventory, a soil pollutant emission inventory for the site is obtained. The soil pollutant emission inventory includes: the characteristic pollutants, the main emission pathways, the polluted areas, and the pollution contribution under multiple dimensions; the multiple dimensions include: temporal scale, spatial scale, main pollution process sources, and main emission pathways.

[0012] Preferably, the establishment of the soil pollution account further includes: determining the distribution parameters of the characteristic pollutants in the polluted area.

[0013] Preferably, the allocation parameters include: range allocation parameters or pollution amount allocation parameters;

[0014] The allocation parameters within the specified range are allocated according to the spatial size within the contaminated area.

[0015] The pollution allocation parameters are as follows: within the pollution area, the pollution is allocated according to the emission intensity.

[0016] Preferably, the calculation of the emission amount of the characteristic pollutant in the calculation table specifically involves: calculating the emission amount of the characteristic pollutant using an emission pathway function;

[0017] The emission pathway functions include: atmospheric emission function, wastewater emission function, solid waste emission function, and leakage / spillage function.

[0018] Preferably, the atmospheric emission function is specifically:

[0019] D i，atmo =P i，atmo ×k i，atmo ×(1-γ);

[0020] Where D i，atmo The pollutant flux emitted from the industrial section of product i to the site via atmospheric deposition; k i，atom γ is the site settlement ratio coefficient, and γ is the end-of-pipe treatment efficiency.

[0021] Preferably, the wastewater discharge function is specifically:

[0022]

[0023] Where D i，wastewater(pipe) With D i，watewater(pond) These represent the amounts of pollutants discharged into the soil via pipeline leaks and treatment pond leaks, respectively, k i，wastewater(pipe) With k i，wastewater(pond) L represents the leakage coefficient under normal operating conditions for the pipeline and the treatment tank, respectively. pipe S is the length of the pipe. pond c is the area of ​​the pool. i，wastewater The concentration of pollutants in the treatment tank is denoted by T, and the operating time is denoted by T.

[0024] Preferably, the solid waste emission function is specifically:

[0025]

[0026] Where D i，waste c represents the amount of pollutants discharged into the soil via leachate from solid waste stockpiles. i，waste The concentration is obtained from the pollutant leaching test, I leaching η is the ratio of the amount of contaminants leached by dynamic leaching to the amount leached by static leaching. waste This is an open management factor obtained by the enterprise.

[0027] Preferably, the leakage function is specifically:

[0028]

[0029] D i，tank Q is the amount of pollutants discharged into the soil through leakage from storage tanks. L Where is the fluid flow rate, T is the duration of a single leak, Y is the operating life of the storage tank, N is the number of storage tanks, and η is the number of tanks. tank It is the leakage frequency, c i，j It refers to the concentration of this pollutant in the liquid inside the storage tank, C. d These are the parameters for liquid leakage: A is the area of ​​the rupture, ρ is the density of the leaking liquid, P is the pressure of the medium in the storage tank, P0 is the air pressure, g is the acceleration due to gravity, and h is the liquid level above the rupture point.

[0030] Preferably, the determination of the characteristic pollutant specifically includes:

[0031] Obtaining company background information;

[0032] Acquisition of industrial activity information;

[0033] Information obtained through on-site surveys;

[0034] The enterprise background information, industrial activity information, and on-site survey information are integrated and analyzed to determine the characteristic pollutants of the site.

[0035] Preferably, the process of creating the emissions inventory specifically includes:

[0036] Based on the polluted area and the analytical scale, the time-based emission inventory is plotted in tabular form.

[0037] Based on the polluted area and the analysis scale, a certain time point and spatial scale are determined, and corresponding table data are calculated and distributed to each polluted area in a preset allocation method.

[0038] Preferably, after obtaining the soil pollutant emission inventory, the method further includes: verifying the pollutant emission inventory.

[0039] According to a second aspect of the present invention, a system for constructing an industrial site soil pollutant emission inventory is provided, comprising:

[0040] A characteristic pollutant determination unit is used to determine the characteristic pollutants of the site based on site information;

[0041] The soil pollution account establishment unit is used to analyze the site information, identify the main sources of pollution, main emission pathways, and polluted areas of the site, and establish the soil pollution account of the site based on the characteristic pollutants, the main sources of pollution, the main emission pathways, and the polluted areas.

[0042] The calculation table acquisition unit is used to determine the analysis scale for the pollution discharge situation recorded in the soil pollution account at both time and spatial scales; and to calculate the discharge amount of the characteristic pollutants at the corresponding time and spatial scales to obtain the calculation table.

[0043] The emission inventory drawing unit is used to draw time-based emission inventories and spatial emission inventories based on the pollution area and the analysis scale.

[0044] The soil pollutant emission inventory acquisition unit is used to obtain the soil pollutant emission inventory of the site by performing data analysis on the calculation table, the time emission inventory, and the spatial emission inventory. The soil pollutant emission inventory includes: the characteristic pollutants, the main emission pathways, the polluted areas, and the pollution contribution under multiple dimensions; the multiple dimensions include: time scale, spatial scale, main pollution process sources, and main emission pathways.

[0045] The present invention provides a method and system for constructing industrial site soil pollutant emission inventories, which establishes a quantitative relationship between soil pollution emissions and multi-source information of enterprises and actual site surveys, enabling a certain degree of analysis of historical dynamic soil pollution at the site; and allowing the source tracing results to be refined to specific processes or pollution pathways; describing the change process of soil pollution emissions across multiple dimensions, including time scale, spatial scale, scale of major pollution process sources, and scale of major emission pathways, constructing time-based and spatial emission inventories that can comprehensively reflect the entire production process; and simultaneously, conducting analysis at a relatively low cost, forming a widely applicable soil pollution analysis method.

[0046] In one optional embodiment of the present invention, range allocation parameters and pollution quantity allocation parameters are proposed, making the pollution quantity classification of the polluted area more reasonable.

[0047] In one alternative embodiment of the present invention, the pollution emissions under this method are calculated using an innovative atmospheric emission function, and the estimation is performed in a relatively simple way, making it easy to calculate the maximum value formed by atmospheric deposition.

[0048] In one alternative embodiment of the present invention, the pollution discharge amount under this method is calculated using an innovative wastewater discharge function, which facilitates the accurate calculation of the pollution discharge amount under this method.

[0049] In one alternative embodiment of the present invention, an innovative solid waste emission function is used to calculate the pollution emissions under this method. The function is set to calculate how much pollutant can be leached from the same amount of dry matter, and the concept of dynamic leaching is incorporated. The function is set based on the assumption that the amount of pollutants leached by dynamic leaching is less than or equal to that of complete leaching, which facilitates the accurate calculation of pollution emissions under this method.

[0050] In one alternative embodiment of the present invention, the pollution emissions under this method are calculated using an innovative run-through / leakage function, which facilitates the accurate calculation of pollution emissions under this method. Attached Figure Description

[0051] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0052] Figure 1 A flowchart illustrating a method for constructing an industrial site soil pollutant emission inventory according to an embodiment of the present invention;

[0053] Figure 2This is a flowchart illustrating the determination of a specific pollutant according to an embodiment of the present invention;

[0054] Figure 3a This is a schematic diagram illustrating the method for verifying a soil discharge inventory according to an embodiment of the present invention, in which the site is divided into a 1m × 1m grid.

[0055] Figure 3b This is a schematic diagram illustrating the estimation of pollutant content in each soil layer in a soil emission inventory verification method according to an embodiment of the present invention.

[0056] Figure 3c This is a schematic diagram illustrating the estimation of pollutant content in a single soil column using a soil emission inventory verification method according to an embodiment of the present invention.

[0057] Figure 3d This is a schematic diagram illustrating the estimation of pollutant content across the entire site using a soil discharge inventory verification method according to an embodiment of the present invention.

[0058] Figure 4 This is a schematic diagram of a system for constructing an industrial site soil pollutant emission inventory according to an embodiment of the present invention;

[0059] Figure 5 This is a schematic diagram illustrating the content distribution of five heavy metals in 12 POIs, as an example of the present invention.

[0060] Figure 6 This is a schematic diagram of a contaminated area as an example of the present invention;

[0061] Figure 7 This is a schematic diagram of a time-based emissions inventory (four periods from 1972 to 2010) as an example of the present invention;

[0062] Figure 8 This is a schematic diagram of a 10m × 10m (grid) spatial emission inventory (total arsenic flux between 1972 and 2010) as an example of the present invention;

[0063] Figure 9a This is an example of the As pollution distribution map in the first layer (0.0m-0.5m) of the site according to the present invention;

[0064] Figure 9b This is an example of the As contamination distribution map in the second layer (0.5m-1.0m) of the site according to the present invention;

[0065] Figure 9c This is an example of the As contamination distribution map in the third layer (1.0m-1.5m) of the site according to the present invention;

[0066] Figure 9d This is an example of the As contamination distribution map in the fourth layer (1.5m-2.0m) of the site according to the present invention;

[0067] Figure 9e This is an example of the As contamination distribution map in the fifth layer (2.0m-4.0m) of the site according to the present invention;

[0068] Figure 9f This is an example of the As contamination distribution map in the sixth (4.0m-6.0m) layer of the site, which is an example of the present invention.

[0069] Figure 9g This is an example of the As contamination distribution map in the seventh layer (6.0m-8.0m) of the site, which is an example of the present invention.

[0070] Figure 9h This is a site total As pollution distribution map as an example of the present invention;

[0071] Explanation of reference numerals in the attached figures:

[0072] 1-Characteristic pollutant identification unit,

[0073] 2- Soil pollution accounting unit,

[0074] 3-Calculate the cells obtained from the table.

[0075] 4-Emissions Inventory Drawing Unit

[0076] 5- Soil pollutant emission inventory acquisition unit. Detailed Implementation

[0077] 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.

[0078] In the description of this invention, it should be understood that the terms "upper part", "lower part", "upper end", "lower end", "lower surface", "upper surface", etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the invention.

[0079] In the description of this invention, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature.

[0080] In the description of this invention, "a plurality of" means multiple, such as two, three, four, etc., unless otherwise explicitly specified.

[0081] In the description of this invention, unless otherwise explicitly specified and limited, the term "connection" and other such terms should be interpreted broadly. For example, it can refer to a fixed connection, a detachable connection, or an integral connection; it can refer to a mechanical connection, an electrical connection, or a connection that allows communication between the components; it can refer to a direct connection or an indirect connection through an intermediate medium; it can refer to the internal communication between two components or the interaction between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances.

[0082] The technical solution of the present invention will be described in detail below with reference to specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments.

[0083] In one embodiment, a method for constructing an industrial site soil pollutant emission inventory is provided; please refer to [reference needed]. Figure 1 It includes:

[0084] S1: Identify characteristic pollutants: Based on site information, identify the characteristic pollutants of the site;

[0085] S2: Establish a soil pollution account: Analyze site information to identify the main sources of pollution, main emission pathways, and contaminated areas; establish a soil pollution account for the site based on characteristic pollutants, main sources of pollution, main emission pathways, and contaminated areas;

[0086] S3: Obtain the calculation table: Determine the analysis scale for the pollution discharge recorded in the soil pollution account at the time and spatial scales; calculate the discharge of characteristic pollutants at the corresponding time and spatial scales to obtain the calculation table;

[0087] S4: Draw emission inventories: Based on the pollution area and analysis scale, draw time-based and spatial emission inventories;

[0088] S5: Obtain a soil pollutant emission inventory: By analyzing the calculated tables, time emission inventory, and spatial emission inventory, a soil pollutant emission inventory of the site is obtained. The soil pollutant emission inventory includes: characteristic pollutants, main emission pathways, polluted areas, and pollution contribution under multiple dimensions; multiple dimensions include: time scale, spatial scale, main pollution process sources, and main emission pathways.

[0089] In one embodiment, please refer to Figure 2 S1 specifically includes:

[0090] S11: Obtaining company background information;

[0091] Referring to the "Technical Regulations for Information Collection on Land Use Survey of Key Industry Enterprises (Trial Implementation)," the information covers the basic enterprise information, pollution source information, migration pathway information, sensitive receptor information, and existing environmental survey and monitoring information described in the regulations. Please refer to Table 1.

[0092] Table 1 Company Background Information

[0093]

[0094] S12: Industrial Activity Information;

[0095] By collecting production and emissions data from the site;

[0096] S13: Obtaining information from actual site surveys;

[0097] Referring to the "Technical Guidelines for Soil Pollution Status Investigation of Construction Land" (HJ 25.1—2019) and the "Technical Guidelines for Risk Control and Remediation Monitoring of Soil Pollution in Construction Land" (HJ 25.2—2019), the main work includes the work in the guidelines other than sampling. At the same time, in order to determine the sampling of characteristic pollutants in the soil, the preliminary sampling and monitoring methods in the guidelines should be referred to. The number and density of soil monitoring points can be reduced or based on the preliminary sampling and monitoring methods according to the actual situation.

[0098] S14: Integrate and analyze enterprise background information, industrial activity information, and on-site survey information to determine the characteristic pollutants of the site;

[0099] By combining enterprise background information (such as industry sub-category, specific processes and pollution-generating processes), industrial activity information, and actual site survey information, and by combining the data obtained from sampling and analysis (to conduct descriptive statistics and compare with standards to determine soil pollutants exceeding the standards).

[0100] In one embodiment, S2 specifically includes: To meet the need for estimating industrial site pollution, considering the randomness, unpredictability, and complexity of soil pollution distribution, pollution patterns are divided into four main channels (atmosphere A, wastewater B, solid waste C, and leakage D). The main identification methods are derived from the "Handbook of Production and Discharge Coefficients of Industrial Pollution Sources in the First National Pollution Source Census," the "Handbook of Production and Discharge Accounting Methods and Coefficients for Emission Source Statistical Survey," and soil pollution pathways discovered through site records or actual site investigations. Based on existing industrial soil pollution investigation standards, a data collection and estimation process using standardized soil pollution accounting is implemented to lay a foundation for achieving refined soil media emission flux. S2 also includes: confirming the allocation parameters of the polluted area (i.e., the proportion of pollution allocated to the polluted area); for soil pollution with unclear geographical boundaries, it can also be recorded, and suitable methods can be found for allocation. Therefore, the concepts of range allocation parameters and pollution allocation parameters are proposed, and the specific parameter tables and parameter explanations are presented in Table 2.

[0101] Table 2 Soil Pollution Account

[0102]

[0103]

[0104] Please refer to Table 2. The allocation parameters within the scope are spatial allocation weights. Each soil pollution account describes the polluted area and allocates it according to spatial conditions, such as allocation within a grid based on factors like area and length in the emission inventory. Specifically, for a pollution pathway recorded in a soil pollution account, its pollution range (points, lines, areas) is defined, and the emission volume is allocated to that polluted area in a certain way (either averaged or non-averaged).

[0105] Because soil pollution emission inventories are very small in scale, the differences in emission amounts caused by different emission intensities within a small area cannot be ignored. To address this issue, pollution allocation parameters have been proposed. For example, a pipeline may have n orifices with different emission intensities. If sufficient data is available, the pollution amount can be allocated to these points separately, or according to the orifice size or other parameters related to emission amounts. This is what pollution allocation parameters are.

[0106] The range allocation parameter can be used when allocating emissions at a fine spatial scale; the pollution amount allocation parameter can be used when there are non-uniform intensity emissions and careful consideration is required. If the emissions are uniform, they can be averaged.

[0107] In one embodiment, taking As as a characteristic pollutant as an example, the calculation of the emission amount of the characteristic pollutant in S3 specifically includes: calculating the source strength of the industrial emission of the enterprise.

[0108] Methods for calculating the source strength of industrial emissions from enterprises:

[0109]

[0110] Where P i，j For Class i industrial products (industrial pollution source), the amount (t) of pollutants (such as As) emitted through Class j emission modes (Atmospheric emission; Wastewater emission; Solid waste emission; Tank leakage) is expressed as a quantity of pollutants. For the production volume of industrial product type i, t; F i，j The coefficient for the amount of pollutant carrier (e.g., wastewater, dust) of type j, generated from type i industrial products, is derived from the "Handbook of Industrial Pollution Source Production and Discharge Coefficients from the First National Pollution Source Census" and the "Handbook of Production and Discharge Accounting Methods and Coefficients from the Statistical Survey of Emission Sources," where t is the pollutant carrier and t is the output, and c is the output. i，j The specific content parameters of pollutants (such as As) obtained from the literature are t pollutant / t pollutant carrier.

[0111] Note: The function derivation method is based on the production volume of a certain industrial product. Multiply by F i，j The amount of sewage medium can be obtained, then multiplied by c. i，j That is, the amount of As discharged.

[0112] Preferably, the emission amounts of characteristic pollutants can be calculated using a pollution pathway correspondence model. For example, four types of emission pathway functions can be used: atmospheric emission function, wastewater emission function, solid waste emission function, and leakage / spillage function. For sites with sufficient industrial activity data on pollution pathways, the source strength P of the pollution emission can be directly obtained. i，j P obtained directly i，j Replace P in the atmospheric emission function i，Atmo Solid waste emission function The four emission pathway functions are described in detail below.

[0113] (1) Atmospheric emission function:

[0114] D i，atmo =P i，atmo ×k i，atmo ×(1-γ);

[0115] Where D i，atmo The flux of pollutants (such as As) emitted from the industrial section of product i into the site via atmospheric deposition is expressed in t. i，atomγ is the site settlement ratio coefficient (dimensionless, ranging from 0 to 1; the approximate estimate or model estimate such as AERMOD is used depending on the amount of emissions from atmospheric emission pathways); γ is the end-of-pipe treatment efficiency (dimensionless); since atmospheric emission pathways account for 1% in this example and atmospheric model parameters are lacking, k is... i，atom The values ​​are approximate estimates and are not finely segmented.

[0116] Note: The function derivation method is as follows, P i，atmo Multiply by (1-γ) to get the amount of As emitted through chimneys and other pathways, and then multiply by k. i，atmo This yields the total amount that can settle on the site. And k i，atmo The maximum value is 1, meaning all emissions fall within the site. The specific proportion can be estimated using an atmospheric model. Since atmospheric deposition is generally small, a maximum estimate of 1 can be used. When using the maximum value estimation function, uncertainty arises from overestimating emissions and the amount entering the soil.

[0117] Existing methods typically use models such as AERMOD to simulate pollutant deposition (requiring a large amount of meteorological and emission data), which can accurately estimate the extent and flux of deposition. However, due to the large amount of data required and the high cost, and the fact that the amount of soil pollution emitted by atmospheric deposition is usually relatively small, the atmospheric emission function proposed in the above embodiments of this invention estimates the maximum value of atmospheric deposition in a relatively simple way, which is convenient for calculation.

[0118] (2) Wastewater discharge function:

[0119] Wastewater discharge pathways are categorized into pipeline leaks and treatment pool leaks. Based on the allowable leakage under normal operating conditions, an estimate is made, which serves as the basis for obtaining a wastewater leakage model.

[0120] D i，wastewater(pipe) =T×k i，wastewater(pipe) ×L pipe ×c i，wastewater ;

[0121] D i，wastewater(pond) =S pond ×k i，wastewater(pond) ×T×c i，wastewater ;

[0122] Where D i，wastewater(pipe) With D i，wastewater(pond) These represent the amounts of pollutants (e.g., As) discharged into the soil via leaks from pipelines and treatment ponds, in tons (t); k. i，wastewater(pipe) With k i，wastewater(pond) Leakage coefficients under normal operating conditions for the pipeline and treatment tank are L / (h·km) and L / (d·m), respectively. 2 );L pipe S represents the pipe length, in meters (m); pondThe area of ​​the pool is in meters (m); c i，wastewater The concentration of As in the pool is mg / L; T is the operating time, d.

[0123] D i，wastewater(pipe) The allowable leakage rate is derived from the allowable leakage rate in the "Standard for Construction and Acceptance of Water Supply and Drainage Pipeline Engineering" (GB50268-2008), and is assumed to be within the allowable range.

[0124] D i，wastewater(pond) The allowable leakage rate in the "Standard for Construction and Acceptance of Water Supply and Drainage Structures" (GB50141-2008) is derived from the assumption that it is within the allowable range.

[0125] Existing research on pipeline leakage relies heavily on precise calculations using simulation software like FLUENT, with limited attention paid to historical site pipeline calculations and even fewer studies on specific calculations of water tank leakage. The wastewater discharge function proposed in the embodiments of this invention estimates both pipeline network leakage and treatment tank leakage.

[0126] (3) Solid waste emission function:

[0127]

[0128] Where D i，waste The amount of pollutants (such as As) discharged into the soil medium via solid waste leaching, in tons; c i，waste The concentration obtained from the As leaching test, mg / kg; I leaching This is a parameter specific to this embodiment, used to calculate the percentage of statically leached arsenic released under dynamic leaching conditions. It is the ratio of the amount of As leached statically to the dynamic leaching ratio (obtained from the total dynamic leaching / static leaching amount), expressed as mg of arsenic leached dynamically / mg of arsenic leached statically; η waste The management factor obtained by the enterprise is derived from the open storage rate and is dimensionless.

[0129] Note: The derivation process is as follows: The amount of solid waste or raw materials discharged is obtained, multiplied by c. i，waste The leaching concentration under static leaching experiments yields the maximum amount of all solid waste that can be leached, multiplied by I. leaching The dynamic leaching ratio compared to the static leaching ratio gives the amount of dynamic leaching, which can then be multiplied by η. waste By opening the management factors, the amount of solid waste leaching under dynamic leaching under open-air conditions can be obtained.

[0130] Existing technologies for solid waste leaching typically rely on experimental research papers, rarely calculating emissions from real-world sites. The solid waste emission function proposed in the embodiments of this invention can accurately calculate how much pollutant can be leached from the same amount of dry matter. Furthermore, it incorporates the concept of dynamic leaching, creatively setting a function based on the assumption that the amount of pollutants leached through dynamic leaching is less than or equal to that leached completely.

[0131] (4) Leakage function:

[0132]

[0133] D i，tank Q is the amount (As) discharged into the soil via leakage from the storage tank, expressed in tons (t). L It is the fluid flow rate, kg / s (10 -3 t / s); T is the duration of a single leak, s (obtained from the guidelines); Y is the service life of the storage tank, a; N is the number of storage tanks; η tank This is the leakage frequency, / a (obtained from the guidelines). The guidelines provide a relatively small probability. If the company has a known history of leakage, it can be calculated using the company's annual leakage frequency, or estimated based on once a year; C d This is the liquid leakage parameter (dimensionless) (selected by the guidelines, maximum value); A is the area of ​​the damage, m. 2 ρ is the density of the leaked liquid, kg / m³ 3 P is the pressure of the medium in the storage tank, in Pa; P0 is the air pressure, in Pa; g is the acceleration due to gravity, 9.81 m / s². 2 h is the liquid level above the break, m.

[0134] Note: This formula is derived based on the "Technical Guidelines for Environmental Risk Assessment of Construction Projects HJ 169-2018", Q L The formula is derived from the guidelines and will not be elaborated further. Multiply by the single leakage time to get the single leakage volume, then multiply by c. i，j That is, the amount of As leaked in a single incident, multiplied by YN, gives the situation for multiple storage tanks over many years, multiplied by η tank The total As amount leaked from multiple storage tanks over many years was obtained, where η tank The value can be selected according to the factory's situation. The value in the guide is a probability over a wide range, which may be relatively small for a single site.

[0135] In one embodiment, S4 specifically includes: based on the calculation table, and combined with the polluted areas recorded in the soil pollution ledger and the selected analysis scale, drawing a temporal emission inventory and a spatial emission inventory. The temporal emission inventory can present the data from the calculation table (three dimensions: time period, emission pathway and process, and soil pollution emission amount) using Excel, Origin, or similar methods. The spatial emission inventory requires determining a specific point in time and spatial scale, corresponding to the data in the calculation table, and allocating it to each polluted area using methods such as average or non-average distribution. It is then presented using geographic information mapping tools such as ArcMAP, for example, displaying the soil pollution emission amount in each 10m × 10m grid of the site.

[0136] In one embodiment, S5 is followed by;

[0137] S6: Verify the obtained soil pollutant emission inventory;

[0138] Specifically, (1) for soil pollutants with weak migration and transformation capabilities, such as arsenic, the principle that the emission amount calculated from the soil pollutant emission inventory is approximately equal to the amount of pollutants calculated from actual sampling can be used for verification; whereby the method for calculating the amount of pollutants from actual sampling is referred to Figures 3a-3d , Figure 3a To divide the site into a 1m x 1m grid, Figure 3b To estimate the pollutant content of each soil layer, Figure 3c To estimate the contaminant content of a single soil column; Figure 3d To estimate the pollutant content of the entire site. (2) For soil pollutants with strong migration and transformation capabilities, the migration and transformation can be simulated by combining the Hydrus and fugacity models, and verified by combining the verification method in (1).

[0139] In one embodiment, a system for constructing an industrial site soil pollutant emission inventory is also provided; please refer to [reference needed]. Figure 4 It includes: Characteristic pollutant identification unit 1, soil pollution account establishment unit 2, calculation table acquisition unit 3, emission inventory drawing unit 4, and soil pollutant emission inventory acquisition unit 5.

[0140] Among them, the characteristic pollutant determination unit 1 is used to determine the characteristic pollutants of the site based on the site information;

[0141] The Soil Pollution Account Establishment Unit 2 is used to analyze site information, identify the main sources of pollution, main emission pathways, and polluted areas of the site, and establish a soil pollution account for the site based on characteristic pollutants, main sources of pollution, main emission pathways, and polluted areas.

[0142] Unit 3, which obtains the calculation table, is used to determine the analysis scale for the pollution discharge recorded in the soil pollution account at both time and spatial scales; and to calculate the discharge of characteristic pollutants at the corresponding time and spatial scales to obtain the calculation table.

[0143] Emission inventory drawing unit 4 is used to draw time-based and spatial emission inventories based on pollution areas and analysis scales;

[0144] Unit 5, which obtains the soil pollutant emission inventory, is used to analyze data from the calculation tables, time emission inventory, and spatial emission inventory to obtain the soil pollutant emission inventory of the site. The soil pollutant emission inventory includes: characteristic pollutants, main emission pathways, polluted areas, and pollution contribution under multiple dimensions. The multiple dimensions include: time scale, spatial scale, main pollution process sources, and main emission pathways.

[0145] The feasibility and reliability of constructing and implementing the industrial site soil pollutant emission inventory system described in the above embodiments are verified below with specific examples.

[0146] In this example, the site is a phosphate fertilizer-sulfuric acid production site. Its overview is that it mainly produces four products and processes, covering superphosphate production, compound fertilizer production, sulfuric acid from pyrite, and sulfuric acid from sulfur. It also involves some organic matter production processes such as aniline. The site area is approximately 80,000 square meters. The specific process is as follows:

[0147] S1: Analyze enterprise background information, industrial activity information, and on-site survey information to determine the characteristic pollutants of the site;

[0148] The company background information for the site was obtained according to the required forms, and certain industrial activity information regarding production volume was obtained through previous site investigation reports (Table 3). An on-site investigation of the contaminated areas was also conducted. A total of 73 sampling points (506 samples) were collected, including 12 points of interest (POIs: D-1 to D-12) and 61 validation sampling points (D-13 to D-73). Each point was sampled at seven layers: 0.0–0.5m, 0.5–1.0m, 1.0–1.5m, 2.0–4.0m, 4.0–6.0m, and 6.0–8.0m. Five heavy metals and 31 volatile organic compounds (VOCs) were detected in the POIs. Based on the process characteristics, 12 POIs (12 × 7 layers = 84 samples) were selected and analyzed for V, Ni, Pb, As, and Hg using OriginPro 2022. The data were analyzed, and descriptive statistical results were obtained (Table 4). Figure 5 The results showed that arsenic was the main pollutant distributed across all soil layers. Therefore, for the remaining 61 samples, the As content was primarily determined to obtain detailed As pollution levels and distribution.

[0149] Table 3. Production Capacity of Site Process in Examples

[0150] Process (output, t) 1972–1989 1990–2000 2000–2005 2006–2010 Superphosphate (dilute acid method) 900000 550000 250000 250000 Sulfuric acid production from pyrite 0 440000 200000 200000 sulfuric acid production 0 0 150000 150000 Compound fertilizer - slurry method 0 0 0 1000000

[0151] Box plots for five heavy metals were generated from the analysis of 12 POIs samples. Please refer to the provided text. Figure 5 The table shows the content distribution of five heavy metals in 12 POIs. According to the primary screening values ​​(GB 36600-2018) (Table 4), As had the highest exceedance rate, and As was selected as the characteristic pollutant of the soil.

[0152] Table 4 shows the descriptive statistics of the five heavy metals in 12 POIs in the examples.

[0153] element Minimum value Maximum value average value variance coefficient of variation First type of screening value mg / kg mg / kg mg / kg mg / kg mg / kg V 12.00 182.00 78.54 26.63 0.34 165 Ni 3.00 189.00 43.71 35.66 0.82 150 Pb 15.00 1150.00 107.9 207.8 1.93 400 As 5.70 1100.00 70.19 171.1 2.44 20 Hg 0.00 11.00 0.34 1.27 3.75 8

[0154] S2: Identify the main pollution pathways at the site and complete the soil pollution record.

[0155] Based on the investigation of the "Handbook of Production and Discharge Coefficients of Industrial Pollution Sources in the First National Pollution Source Census" and the "Handbook of Production and Discharge Accounting Methods and Coefficients for Statistical Survey of Emission Sources," combined with the actual site survey, soil pollution records were filled out, and polluted areas were recorded using GIS (please refer to...). Figure 6 ).

[0156] This site in the example covers nine main pollution modes, which are mainly divided into air pollution, wastewater pollution, solid waste leaching pollution, and spillage pollution. For the specific pollution modes of the phosphate fertilizer site, please refer to Tables 5-1 and 5-2. Due to the page size constraints, Table 5 has been divided into two. If the page display is unrestricted, the two tables can be merged.

[0157] Table 5-1 Example of Site Contamination Investigation Results (Based on Soil Contamination Account)

[0158]

[0159]

[0160] Table 5-2 Example Site Contamination Investigation Results (Based on Soil Contamination Account)

[0161]

[0162]

[0163] S3: Determine the analysis scale and perform calculations using four emission pathway functions to obtain the calculation table.

[0164] Based on the requirements of this case, a grid scale of 10m×10m was divided, and calculations were performed using the four emission pathway functions in S3 through five tables: the atmospheric emission calculation table in Table 6, the wastewater emission calculation table in Table 7, the solid waste emission calculation table in Table 8, the run-in-leak emission calculation table in Table 9, and the parameter selection and source of the calculation table in Table 10.

[0165] (1) Atmospheric emissions

[0166] Table 6 shows the atmospheric emission calculation table in the embodiments.

[0167]

[0168]

[0169] Note: Letters such as "a", "b", and "c" are reference labels, see Table 10.

[0170] (2) Wastewater discharge

[0171] Table 7 Wastewater Discharge Calculation Table in Examples

[0172]

[0173]

[0174] Wastewater discharge calculation table (including parameter selection and reference). Letters "a", "b", "c", etc. are reference labels, see Table 10.

[0175] __________________

[0176] (3) Solid waste discharge

[0177] Table 8 shows the solid waste emission calculation table in the embodiments.

[0178]

[0179] Solid waste emission calculation table (including parameter selection and reference). Letters "a", "b", "c", etc. are reference labels, see Table 10.

[0180] (4) Leakage

[0181] Table 9 shows the calculation table for leakage and spillage in the embodiments.

[0182]

[0183]

[0184]

[0185] Note: If ηtank is set to 1, the total amount is only 0.002t, so the discharge from the storage tank is very small.

[0186] Tank Leakage Calculation Table (including parameter selection and reference). Note: Letters "a", "b", "c", etc. are reference labels, see Table 10.

[0187] (5) Parameter selection and sources in the calculation of emission pathway functions for the four categories

[0188] Table 10 shows the selection and source of calculation table parameters in the embodiments.

[0189]

[0190]

[0191]

[0192] S4: Based on the calculation tables, and combined with the polluted areas recorded in the soil pollution ledger and the selected analysis scale, construct the temporal and spatial emission inventories. The temporal emission inventory uses the data from the calculation tables presented in Origin format; please refer to [reference needed]. Figure 7 The spatial emissions inventory needs to be calculated using the current time and with a spatial scale of 10m × 10m. Correspondingly, the data in the calculation table is allocated to each contaminated area in this embodiment using an average distribution method and presented using ArcMAP. Please refer to [the relevant documentation / reference]. Figure 8 .

[0193] S5: By performing data analysis (mainly descriptive statistics in this embodiment) and comprehensive information analysis on the calculated spatial emission inventory, temporal emission inventory, and calculation tables, combined with enterprise background information and actual site survey information, it can be determined that the main source of pollution at the site is sulfuric acid production from pyrite, and the main pollution pathway is solid waste dumping, providing a basis and support for subsequent site management.

[0194] This example study investigated four pollution pathways: atmospheric emissions, wastewater discharge, solid waste discharge, and leaks / leaks. Since the study primarily focused on historical sites no longer in operation, usable production data was difficult to obtain. Therefore, we mainly used emission coefficients listed in the "Handbook of Production and Discharge Coefficients of Industrial Pollution Sources from the First National Pollution Source Census" and the "Handbook of Production and Discharge Accounting Methods and Coefficients from the Statistical Survey of Emission Sources," which approximate the actual node emissions. The site's atmospheric emissions were 0.032 t, wastewater discharge was 0.884 t, and solid waste leaching was 23.682 t (when I... leaching (55.48%), and the leakage from the storage tank was 0.02t. (I)leaching When the emissions are 40%, 50%, 60%, and 100%, the solid waste emissions are 17.07t, 21.34t, 25.61t, and 42.69t, respectively. (i) Atmospheric emissions: The calculated emissions are negligible compared to other pollution pathways. Therefore, we can simplify the allocation of atmospheric emissions. In this example, atmospheric emissions are evenly distributed across the entire site. (ii) Wastewater discharge: Since there are no records of pipe leaks, this example assumes that pipe leaks are evenly distributed along the entire line and comply with standards. There were no previous leaks in the wastewater tank, so this example assumes that wastewater tank leaks are permissible and comply with standards (MOHURD, 2008). (iii) Solid waste: In this example, the waste leaching area is divided into three zones (C1, C2, C3), and the corresponding pollution emissions are allocated to these zones according to the calculation tables. (iv) Tank leaks: In this example, pollution emissions are evenly distributed to zones (D1 and D2).

[0195] Atmospheric deposition mainly originates from industrial dust emissions generated by the superphosphate process (A1); ore crushing in the compound fertilizer process also releases As (A2). Wastewater from pipelines and wastewater treatment ponds (B1, B2) originates from multiple processes (sulfuric acid production from pyrite, sulfuric acid production from sulfur, and compound fertilizer production). These processes (sulfuric acid production from pyrite, sulfuric acid production from sulfur, and compound fertilizer production) (C1, C2, C3) primarily generate tailings, phosphogypsum, and other solid wastes, which pollute the soil through stockpiling and leaching.

[0196] Site investigation and soil pollution identification indicate that the site is mainly contaminated by C1 and C3 pollution areas.

[0197] We input pollution survey and industrial activity data (Table S3 and spreadsheets) into our method and then referenced literature to obtain industrial design parameters. Production changes over four periods were used as time parameters. Detailed calculations were performed for four emission pathways over the four time periods: 1972–1989, 1990–2000, 2001–2005, and 2006–2010. The calculation process and parameters are shown in Tables 6 to 10. Figure 7 For time-based emission inventory.

[0198] The amount of As pollution emitted varies considerably from region to region. Figure 8 The main contributors (C1 and C3) are sulfuric acid production from pyrite and compound fertilizer production. These processes generate large amounts of tailings (pyrite slag) and phosphogypsum. In the past few decades, many plants had lax environmental management requirements. The large-scale generation of solid waste (high pollutant concentrations) became a major source of soil pollution through leaching. From 1990 to 2000, C1 contributed 42.24% of total arsenic, representing the highest period of emissions. Figure 7Over four time periods, C1 continuously emitted As into the soil, accounting for 80.59% of total As emissions. The second largest contributor was C3, with compound fertilizer production reaching 1 million tons (2006-2010), and its byproduct, lumpy phosphogypsum, accounting for 15.23% of As emissions. The remaining three emission pathways (atmospheric emissions, wastewater discharge, and leaks) contributed only 3.81% of As pollution emissions. Wastewater discharge was the largest contributor, while atmospheric deposition and leaks accounted for a very small proportion (approximately 0.1%). Therefore, the primary contributor to As emissions was sulfuric acid production from pyrite, followed by compound fertilizer production; the remaining pollution pathways contributed to the remaining As emissions.

[0199] Pollution hotspots are located in C1, C3 and "other organic production areas" ( Figure 9a This is not entirely consistent with the on-site investigation. This is because there are historical records of solid waste landfilling in "other organic production areas." Main hotspots ( Figure 9b Randomly distributed in areas such as Zone C1, "Other Organic Production Areas," raw material storage, and compound fertilizer production areas. Solid waste landfill can also explain this. (Hotspot) Figure 9c The area is clearly located in the pyrite crushing zone and the pyrite-to-sulfuric acid production zone (C1), consistent with the results of the actual site survey. It also verifies the findings of Mao et al. (2022) that the arsenic concentration was high in the raw material area, waste residue area, and sulfuric acid production area. (Contaminated area) Figure 9d The site is located around C1, C2, and B2, but the As content in layer 4 is low, suggesting the possible presence of small-scale random solid waste landfilling. Hotspots in layers 5–7 are located at C1, “Waste Slag Disposal Site 2,” and B2, consistent with the site investigation, and can be interpreted as sulfuric acid production primarily from pyrite. As migration is relatively weak; high concentrations of arsenic in soil are typically associated with industrial activities, consistent with previous studies (Li et al., 2024; Yang et al., 2018). Simultaneously, groundwater contamination is possible at this site, with non-confined aquifers potentially more susceptible to As contamination (Wang et al., 2018). This validates the groundwater contamination record from the site investigation. Repeated disposal of pyrite slag, natural soil, and other waste mixtures can lead to widespread heavy metal contamination. While developed countries have implemented solid waste management measures, pyrite slag dumping remains a common problem in the development process. It is essential to consider developing appropriate solid waste treatment standards and measures.

[0200] S6: This step is optional. The embodiment verifies the soil emission inventory obtained from the site. For soil pollutants with weak migration and transformation capabilities, such as arsenic, the verification method uses the principle that the emission amount calculated from the soil emission inventory proposed in this patent is approximately the same as the pollutant amount calculated from actual sampling. Figures 3a-3d ).

[0201] By calculating 7 layers ( Figures 9a-9gThe total As is calculated and IDW interpolation (1m×1m grid) is performed to estimate the As amount for the entire site. Figure 9h The majority of As pollution emissions are still in C1, which proves that sulfuric acid production using pyrite as a raw material is the decisive pollution pathway.

[0202] It can be observed that C1 is the main hotspot in both maps. Figure 8 and Figure 9h The most severe problems are found in sulfuric acid production areas dominated by pyrite and in pyrite crushing areas. (In the spatial emissions inventory...) Figure 8 In ), hotspot areas are more Figure 9h The large scale is due to a lack of precise information on the depositional area (geological parameters). Other areas within Functional Zone 1 also contain high-intensity pollution emissions, which are difficult to interpret with existing field investigation data. These can be investigated through limited, targeted sampling and analysis.

[0203] The total arsenic value was calculated using IDW interpolation based on the overall As contamination distribution (73 points) across the site (1m×1m).

[0204] Flux mean ×S site = 0.28 kg·(1m×1m grid) -1 ×75264m 2 =21.17t

[0205] Flux mean It is the average value of each soil column in a 1m×1m grid that receives pollution emissions. Figure 9h The flux was obtained from ArcMAP, kg·(1m×1m grid). -1 S site This is the area of ​​the site, in meters. 2 .

[0206] Based on the previous calculations of plant emissions, the calculated arsenic content is 24.62 tons, which differs from the total arsenic content obtained by IDW interpolation at the actual site by +16%. The difference is relatively small. Therefore, the above embodiments of the present invention have good reliability.

[0207] In the description of this specification, the references to terms such as "an embodiment," "an example," "a specific implementation process," and "an example" indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples.

[0208] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for constructing a soil pollutant emission inventory for industrial sites, characterized in that, include: Identify characteristic pollutants: Based on the site information, identify the characteristic pollutants of the site; Establish a soil pollution account: Analyze the site information to identify the main technological sources, main emission pathways, and polluted areas of the site; establish a soil pollution account for the site based on the characteristic pollutants, the main technological sources, the main emission pathways, and the polluted areas; Obtaining the calculation table: For the pollution emissions recorded in the soil pollution ledger, determine the analysis scale based on the time and spatial scales; calculate the emissions of the characteristic pollutants at the corresponding time and spatial scales to obtain the calculation table; specifically, the calculation of the emissions of the characteristic pollutants in the obtained calculation table involves calculating the emissions of the characteristic pollutants using an emission pathway function; the emission pathway function includes a solid waste emission function, which is specifically: in This refers to the amount of pollutants discharged into the soil through leachate from solid waste dumps. The concentration is obtained from pollutant leaching tests. This is the ratio obtained by leaching pollutants in dynamic leaching to the amount leached in static leaching. Open management factors obtained by enterprises; Emissions Inventory: Based on the polluted area and the analytical scale, generate temporal and spatial emission inventories; Obtaining a soil pollutant emission inventory: By analyzing the calculated table, the temporal emission inventory, and the spatial emission inventory, a soil pollutant emission inventory for the site is obtained. The soil pollutant emission inventory includes: the characteristic pollutants, the main emission pathways, the polluted areas, and the pollution contribution under multiple dimensions; the multiple dimensions include: temporal scale, spatial scale, main pollution process sources, and main emission pathways.

2. The method for constructing an industrial site soil pollutant emission inventory according to claim 1, characterized in that, The establishment of the soil pollution account also includes: determining the distribution parameters of the characteristic pollutants in the polluted area.

3. The method for constructing an industrial site soil pollutant emission inventory according to claim 2, characterized in that, The allocation parameters include: range allocation parameters or pollution amount allocation parameters; The allocation parameters within the specified range are allocated according to the spatial size within the contaminated area. The pollution allocation parameters are as follows: within the pollution area, the pollution is allocated according to the emission intensity.

4. The method for constructing an industrial site soil pollutant emission inventory according to claim 1, characterized in that, The emission pathway function also includes: atmospheric emission function, wastewater emission function, and leakage function.

5. The method for constructing an industrial site soil pollutant emission inventory according to claim 4, characterized in that, The atmospheric emission function is specifically as follows: in The amount of pollutants emitted from the industrial section of product i to the site via atmospheric deposition. It is the settlement ratio coefficient within the site. This refers to the efficiency of end-of-pipe emission treatment.

6. The method for constructing an industrial site soil pollutant emission inventory according to claim 4, characterized in that, The wastewater discharge function is specifically as follows: in and These represent the amounts of pollutants discharged into the soil via pipeline leaks and treatment pond leaks, respectively. and These are the leakage coefficients under normal operating conditions for the pipeline and the treatment tank, respectively. For the length of the pipe, The area of ​​the pool, The concentration of pollutants in the treatment tank. This refers to the runtime.

7. The method for constructing an industrial site soil pollutant emission inventory according to claim 4, characterized in that, The aforementioned leakage function is specifically: This refers to the amount of pollutants discharged into the soil through leakage from storage tanks. It is the fluid flow rate. It is the duration of a single leak. It refers to the service life of the storage tank. It refers to the number of storage tanks. It's the leakage frequency. It refers to the concentration of this contaminant in the liquid inside the storage tank. These are parameters related to liquid leakage. It refers to the area of ​​damage. It is the density of the leaked liquid. It refers to the pressure of the medium in the storage tank. It is the air pressure. It is gravitational acceleration. It refers to the liquid level above the damaged area.

8. The method for constructing an industrial site soil pollutant emission inventory according to claim 1, characterized in that, The specific pollutants identified include: Obtaining company background information; Acquisition of industrial activity information; Information obtained through on-site surveys; The enterprise background information, industrial activity information, and actual site survey information are integrated and analyzed to determine the characteristic pollutants of the site.

9. The method for constructing an industrial site soil pollutant emission inventory according to claim 1, characterized in that, The creation of the emissions inventory specifically includes: Based on the polluted area and the analytical scale, the time-based emission inventory is plotted in tabular form. Based on the polluted area and the analysis scale, a certain time point and spatial scale are determined, and corresponding table data are calculated and distributed to each polluted area in a preset allocation method.

10. The method for constructing an industrial site soil pollutant emission inventory according to claim 1, characterized in that, The process of obtaining the soil pollutant emission inventory also includes: verifying the pollutant emission inventory.

11. A system for constructing an industrial site soil pollutant emission inventory, characterized in that, A construction system for realizing the method for constructing an industrial site soil pollutant emission inventory as described in any one of claims 1 to 10; The construction system includes: A characteristic pollutant determination unit is used to determine the characteristic pollutants of the site based on site information; The soil pollution account establishment unit is used to analyze the site information, identify the main sources of pollution, main emission pathways, and polluted areas of the site, and establish the soil pollution account of the site based on the characteristic pollutants, the main sources of pollution, the main emission pathways, and the polluted areas. The calculation table acquisition unit is used to determine the analysis scale for the pollution discharge situation recorded in the soil pollution account at both time and spatial scales; and to calculate the discharge amount of the characteristic pollutants at the corresponding time and spatial scales to obtain the calculation table. The emission inventory drawing unit is used to draw time-based emission inventories and spatial emission inventories based on the pollution area and the analysis scale. The soil pollutant emission inventory acquisition unit is used to obtain the soil pollutant emission inventory of the site by performing data analysis on the calculation table, the time emission inventory, and the spatial emission inventory. The soil pollutant emission inventory includes: the characteristic pollutants, the main emission pathways, the polluted areas, and the pollution contribution under multiple dimensions; the multiple dimensions include: time scale, spatial scale, main pollution process sources, and main emission pathways.

Images