Industrial park maximum allowable emission balance calculation method based on emission increase and decrease balance

Abstract

The invention discloses an industrial park maximum allowable emission balance calculation method based on emission increase and decrease balance. The method comprises the following steps: S1, quantifying environmental concentration change caused by emission reduction of atmospheric pollution sources around an industrial park; s2, simulating environment concentration change caused by the newly added emission source response sequence of the industrial park; s3, determining an atmospheric pollutant environment concentration increment space of the industrial park site; s4, establishing a composite response relationship of the industrial park; s5, solving the maximum allowable discharge under the constraint of the atmospheric pollutant environmental concentration increment of the industrial park site; and S6, checking the influence of the evaluation site. According to the method, the dynamic balance relation of increase and emission reduction of pollution sources inside and outside the park is constructed, accurate control over the total pollutant emission amount of the park is achieved through scientific calculation and dynamic regulation, and technical support is provided for green development of the park.

Description

Technical Field

[0001] The present invention relates to the field of air quality control, and in particular to a method for calculating the maximum allowable emissions of an industrial park based on an emission increase and decrease balance. Background Art

[0002] As my country's industrialization continues to advance, industrial parks, as key vehicles for industrial cluster development, are driving regional economic growth while also facing severe challenges from air pollution. The challenge of ensuring sustained improvement in air quality within these parks while ensuring economic growth has become a pressing issue.

[0003] Traditional methods for controlling total pollutant emissions are often based on the theory of environmental capacity, using source inventories and atmospheric diffusion models to calculate the atmospheric environmental capacity of a park. However, these methods primarily focus on reducing emissions from pollution sources within the park, using a "one-size-fits-all" approach that makes it difficult to account for the differences and development needs of different regions and industries. This can easily lead to a situation where "emission reductions do not increase efficiency," making it difficult to effectively resolve the contradiction between park development and environmental protection, and thus exhibiting certain limitations. Therefore, there is an urgent need to explore a more scientific, precise, and dynamic method for controlling pollutant emissions to achieve a synergistic and win-win situation between economic development and environmental protection. Summary of the Invention

[0004] The purpose of the present invention is to overcome the shortcomings of the existing technology and provide a method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases.

[0005] The object of the present invention is achieved through the following technical solutions: In a first aspect, the present invention provides a method for calculating the maximum allowable emissions of an industrial park based on an emission increase and decrease balance, comprising the following steps: S1: Quantify the changes in ambient concentrations caused by emission reductions from air pollution sources around industrial parks, including: establishing a buffer zone around the industrial park in the simulation grid by setting a mask area; mapping the emission reduction coefficients corresponding to different pollution sources to the emission inventory adjustment coefficients within the mask area; using an air quality model to simulate the emission inventory of air pollution sources before and after emission reductions, and calculating the reduction in ambient concentrations at different monitoring stations caused by emission reductions from air pollution sources; S2: Simulate the changes in environmental concentration caused by the response sequence of new emission sources in industrial parks, including: taking industrial parks as units, gridding the new emissions in industrial parks according to the types of planned land in the industrial parks, and establishing a serialized grid emission inventory that can be applied to air quality model simulation; for serialized grid emission inventories with different emission intensities per unit area, based on the urban air pollution source emission inventory, using the serialized grid emission inventory as the emission increment, and using the air quality model to complete the simulation of the corresponding serialized grid emission inventory, respectively. Calculate the difference between the simulation results of the serialized grid emission inventory for different emission amounts and the simulation results of the urban air pollution source emission inventory, and obtain the corresponding environmental concentration changes under the emission inventory increments of all industrial parks; S3: Determining the incremental space for ambient concentrations of atmospheric pollutants at the industrial park site, including: combining the ambient concentration reductions from the quantitative assessment of emission reduction effectiveness of atmospheric pollution sources surrounding the industrial park calculated in step S1, the current status of the industrial park monitoring sites, and the target atmospheric pollution concentration values ​​for a specific year, to determine the upper limit of concentration increases caused by increased pollution emissions from the industrial park as a basis for calculating the maximum allowable emissions; S4: Establishing a composite response relationship for the industrial park, including: combining the emission intensity per unit area in step S2 with the corresponding environmental concentration changes under the corresponding industrial park emission inventory increment, calculating the response function relationship between the emission increment and environmental concentration of each monitoring site and each industrial park, and the corresponding different pollutants; and integrating the response functions of different industrial parks to obtain a composite response function of the emissions of different industrial parks to the pollutant concentration of a specific monitoring site; S5: solving the maximum allowable emission of atmospheric pollutants at the industrial park site under the constraint of the increase in ambient concentration of atmospheric pollutants, including: using the composite response function integrated in step S4, substituting the emission of different pollutants from different industrial parks under the condition of the upper limit of the concentration increase caused by the increase in pollution emissions from the industrial park in step S3, and solving the concentration contribution of different industrial parks to the specific environmental monitoring point; S6: Verification of the impact of the evaluation site, including: using the composite response function, substituting the maximum allowable emission data of different parks obtained, solving the change of pollutant concentration at the evaluation environment evaluation site, and verifying whether it meets the evaluation environment evaluation site concentration limit pre-set in step S3.

[0006] Furthermore, in step S1, the quantification of the environmental concentration change caused by the emission reduction of air pollution sources around the industrial park includes: S11: Use Geographic Information System (GIS) software to establish an industrial park regional buffer zone based on the industrial park boundary, and set the required buffer zone radius based on the local emission reduction potential and concentration reduction needs; S12: Based on the relationship between the latitude and longitude of the center of the air quality simulation grid and the industrial park buffer zone, the simulation grids within the buffer zone are screened, and an emission reduction mask file for the industrial park buffer zone is established based on the industrial park as a unit; S13: Prepare emission inventory adjustment files based on the emission reduction ratios of different industrial parks. Use the emission model combined with the emission inventory reduction ratios to adjust the air quality model emission inventory data. Align the emission reduction coefficients corresponding to different pollution sources with the emission inventory adjustment coefficients within the mask area. S14: Use the emission inventory before and after emission reduction to carry out air quality simulation and obtain the p Industrial parks e Pollutant concentration before emission reduction C pe0 and concentration after emission reduction C pe1 ; Combined with the concentration of environmental monitoring stations in industrial parks C peo , correct the pollutant concentration simulation results, take the relative concentration change, and calculate the pollutant concentration reduction caused by emission reduction according to the following formula Δ C pe : Evaluation sites for industrial parks E Emission reduction concentration, using the concentration before emission reduction at the evaluation site C Ee0 and concentration after emission reduction C Ee1 , monitoring concentration at evaluation sites C Eeo , and calculate after correction: .

[0007] Furthermore, in step S2, the simulation of the environmental concentration change caused by the response sequence of the newly added emission sources in the industrial park includes: S21: Using the industrial park as the boundary, establish emission inventory data in combination with the land use type structure within the industrial park; S22: Based on historical data, establish a time distribution series for industrial source emissions, including monthly distribution coefficients Fac_m , weekly distribution coefficient Fac_w and hourly distribution coefficient Fac_h ; Establish species partition coefficient Fac_s , to achieve the distribution of complex pollutants to air quality simulation components; S23: A model for processing incremental emissions from industrial parks, which is used to convert annual emission data from industrial parks under a specific emission intensity into hourly emission data with time variations, and to incorporate incremental emissions from industrial parks into the basic emission inventory data of the air quality model; p The emission intensity per unit area of ​​industrial parks EMIS pa Down, i Month, week j 、 k Hours, components l The emission data is calculated as follows: S24: For different industrial parks, different unit area emissions EMIS pa Establish an emission sequence and use the industrial park emission increase model to operate the basic inventory of the air quality model. Add the hourly emission data of each pollutant component newly emitted by the industrial park to the corresponding pollutant component in the basic emission inventory of the air quality model to obtain the scenario emission inventory data. Call the air quality model to complete the simulation of the basic emission inventory and different scenario emission inventory data, according to the following formula, using c Emission inventory simulated concentrations for scenarios C pec Simulated concentrations compared to the base case C peb , Computing Industrial Park p New emission sources e In the industrial park p Pollutant response sequence concentration ΔC pea : Evaluation Site E , extract review sites E Corresponding simulation grid scenario emission inventory concentration C Eec , basic scenario simulation concentration C Eeb , the corresponding pollutant response sequence concentration can be calculated ΔC Eea , the method is as follows: .

[0008] Furthermore, in step S21, in the absence of detailed planning, the average distribution of industrial land is considered and the emission intensity is established as 1t / km 2 ·a unit emission intensity emission inventory data to form a baseline emission inventory.

[0009] Furthermore, in step S3, determining the incremental space of atmospheric pollutant environmental concentration at the industrial park site includes: Combined with the pollutant concentration reduction obtained from the quantitative assessment of the emission reduction effectiveness of air pollution sources around the industrial park in step S1 , using the current concentration of site p in the industrial park C peo0 and target atmospheric pollution concentration values ​​for specific years C peoy Determine the concentration increase limit in the park caused by the increase in pollution emissions in the park , as the basis for calculating the maximum allowable emissions, the calculation formula is as follows: Among them, the current concentration of site p in the industrial park C peo0 Corresponding to the concentration of the industrial park environmental monitoring station in step S1 C peo .

[0010] Furthermore, in step S4, the establishment of the industrial park composite response relationship includes: S41: Combining different EMIS pa Simulation results of the intensity and corresponding industrial park emissions’ contribution to the concentration of the park’s monitoring stations Using cubic natural spline interpolation, we can obtain the response function relationship between emission increment and environmental concentration of different environmental monitoring stations and industrial parks, and different pollutants. C pe =f P (EMIS pa(x) ) , specifically including: The response function between new emissions from industrial parks and environmental concentrations at park sites f p (EMIS pax ) The establishment of the interpolation process needs to be based on the emission intensity of different unit areas. EMIS pa The corresponding response sequence concentration of the park site obtained from the air quality simulation ΔC pea Divide the simulation into segments and calculate the interpolation point spacing of the nth segment in sequence. For example, the first segment needs to use 0 and EMIS pa(1) 、 Δ C pea(1) , establish a piecewise fitting function: The key to establishing the response function is to solve the a 、 b 、 c 、 d The coefficient of the fitting function is continuous; since the solution is between the value of each segment point and the breakpoint The values ​​are equal, and 2n equations are established, among which EMIS pa(n) are all known, namely: … Since the fitting function is continuous and differentiable, the first-order derivatives at each segment point are equal, and the following n-1 equations are established: … At the same time, the second-order derivatives of the functions at each segment point are equal, and the following n-1 equations are established: … Finally, the second-order derivative of the first and last points is 0, which serves as the boundary condition, namely: Since in the above equation, EMIS pa(1) and ΔC pea(1) The unit area park emission intensity input into the model and the corresponding park site concentration increment calculated by the model are both known quantities. Therefore, after combining the above equations into a system of equations, the linear equations can be solved to obtain the corresponding values ​​for different segments. a 、 b 、 c 、 d coefficient, we get f p (EMIS pa(x) ) , and for the functional relationship between industrial park emissions and evaluation sites, it is only necessary to replace Replace with That is, the corresponding functional relationship is obtained fE (EMIS pa(x) ) ; S42: Integrate the response functions of different industrial parks to obtain the composite response function of different industrial park emissions to the pollutant concentration of a specific environmental monitoring station; the contribution concentration of P different industrial parks to the pollutant e at a certain industrial park monitoring station C e Calculate according to the following formula and substitute the unit area emissions of different parks EMIS pa(x) The concentration increment at the corresponding station can be solved C e : .

[0011] Furthermore, in step S5, the maximum allowable emission of atmospheric pollutants at the industrial park site under the constraint of the increase in ambient concentration of atmospheric pollutants is solved, including: Based on the functional relationship established in step S4, the emission intensity data of different industrial parks are introduced EMIS pa(x) To solve, due to the pollution emission intensity adjustment involving multiple industrial parks, the same pollutant concentration increment C e It often corresponds to multiple park emission combinations. In order to cover the emission increment space of different parks, the exhaustive method with maximum value constraints is used to solve it. Using the response function establishment method in step S4, the dependent variable and the independent variable are exchanged to obtain the functional relationship between the emission intensity change and the concentration increment, that is, , using this function and the corresponding concentration increase limit in step S3 That is, calculate the maximum emission when only this park emits EMIS paxm , as the upper limit of the park emission increment, and then set the number of times the response equation is solved S , calculate the incremental step size as follows: From 0 to EMIS paxm By step length Establish an emission increment sequence, orthogonalize it by industrial park, form an input data set, and bring it into the function relationship established in step S4 for solution calculation. The number of calculations required is related to the number of solutions S and the number of parks p, which is Second-rate.

[0012] Furthermore, in step S6, the evaluation site impact verification includes: The solution obtained in step S5 is analyzed and screened. The result of the solution is taken as the solution that meets the conditions, and the value trend of different parks in the actual park development plan is combined to select the final unit area emission increment of different parks. , using the method described in step S4 to establish a functional relationship between park emissions and evaluation sites f E (EMIS pax ) , combined with the current concentration of the evaluation site C Eeo0 , target concentration C Eeoy , reduce emissions and concentrations Carry out the concentration impact verification of the evaluation site. It is considered that the increased emissions from the park have not caused the deterioration of the air quality at the evaluation site. Otherwise, the emissions from the park should be reduced to ensure the contribution to the comprehensive concentration of the evaluation site. Below 0: .

[0013] The beneficial effects of the present invention are: In an exemplary embodiment of the present invention, full use is made of the space for pollutant concentration reduction brought about by emission reduction of pollution sources around the industrial park, and the increase in pollutant concentration caused by the new pollution emissions of the industrial park is offset, so as to develop the industrial park without causing the deterioration of the current atmospheric environment. The core of this method is to establish a dynamic balance relationship between the increase and decrease of pollution sources inside and outside the park, and through scientific measurement and dynamic regulation, to achieve precise control of the total amount of pollutant emissions in the park, and provide technical support for the green development of the park. Subsequently, while ensuring the quality of the environment, by vacating high-pollution, low-output emission sources, the maximum value of the pollution emission increase required for the development of enterprises in the park can be calculated, thereby facilitating green and high-quality development. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] Figure 1 A flowchart of a method for calculating the maximum allowable emissions of an industrial park based on an emission increase and decrease balance provided in an exemplary embodiment of the present invention; Figure 2 A schematic diagram of detailed steps of a method for calculating the maximum allowable emissions of an industrial park based on an emission increase and decrease balance provided in an exemplary embodiment of the present invention; Figure 3 A schematic diagram of delineating an industrial park buffer zone provided in an exemplary embodiment of the present invention; Figure 4 A schematic diagram of an emission reduction mask file for an industrial park buffer zone provided in an exemplary embodiment of the present invention; Figure 5 It is a schematic diagram of function fitting provided in an exemplary embodiment of the present invention. DETAILED DESCRIPTION

[0015] The technical solution of the present invention is described clearly and completely below with reference to the accompanying drawings. It is apparent that the embodiments described are only a portion of the embodiments of the present invention, not all of them. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without inventive effort are intended to fall within the scope of protection of the present invention.

[0016] In the description of the present invention, it should be noted that the directions or positional relationships indicated by terms such as "center," "up," "down," "left," "right," "vertical," "horizontal," "inside," and "outside" are based on the directions or positional relationships described in the accompanying drawings and are intended solely to facilitate and simplify the description of the present invention. They are not intended to indicate or imply that the devices or components referred to must have a specific orientation, be constructed, or operate in a specific orientation. Therefore, they should not be construed as limitations on the present invention. Furthermore, terms such as "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0017] In the description of the present invention, it should be noted that, unless otherwise expressly specified or limited, the terms "mounted," "connected," and "connected" should be understood in a broad sense. For example, they may refer to fixed connections, detachable connections, or integral connections; mechanical connections or electrical connections; direct connections or indirect connections through an intermediate medium; and internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the present invention.

[0018] In addition, the technical features involved in different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0019] See also Figure 1 , Figure 1 A flowchart of a method for calculating the maximum allowable emissions of an industrial park based on an emission increase and decrease balance provided in an exemplary embodiment of the present invention is shown, comprising the following steps: S1: Quantify the changes in ambient concentrations caused by emission reductions from air pollution sources around industrial parks, including: establishing a buffer zone around the industrial park in the simulation grid by setting a mask area; mapping the emission reduction coefficients corresponding to different pollution sources to the emission inventory adjustment coefficients within the mask area; using an air quality model to simulate the emission inventory of air pollution sources before and after emission reductions, and calculating the reduction in ambient concentrations at different monitoring stations caused by emission reductions from air pollution sources; S2: Simulate the changes in environmental concentration caused by the response sequence of new emission sources in industrial parks, including: taking industrial parks as units, gridding the new emissions in industrial parks according to the types of planned land in the industrial parks, and establishing a serialized grid emission inventory that can be applied to air quality model simulation; for serialized grid emission inventories with different emission intensities per unit area, based on the urban air pollution source emission inventory, using the serialized grid emission inventory as the emission increment, and using the air quality model to complete the simulation of the corresponding serialized grid emission inventory, respectively. Calculate the difference between the simulation results of the serialized grid emission inventory for different emission amounts and the simulation results of the urban air pollution source emission inventory, and obtain the corresponding environmental concentration changes under the emission inventory increments of all industrial parks; S3: Determining the incremental space for ambient concentrations of atmospheric pollutants at the industrial park site, including: combining the ambient concentration reductions from the quantitative assessment of emission reduction effectiveness of atmospheric pollution sources surrounding the industrial park calculated in step S1, the current status of the industrial park monitoring sites, and the target atmospheric pollution concentration values ​​for a specific year, to determine the upper limit of concentration increases caused by increased pollution emissions from the industrial park as a basis for calculating the maximum allowable emissions; S4: Establishing a composite response relationship for the industrial park, including: combining the emission intensity per unit area in step S2 with the corresponding environmental concentration changes under the corresponding industrial park emission inventory increment, calculating the response function relationship between the emission increment and environmental concentration of each monitoring site and each industrial park, and the corresponding different pollutants; and integrating the response functions of different industrial parks to obtain a composite response function of the emissions of different industrial parks to the pollutant concentration of a specific monitoring site; S5: solving the maximum allowable emission of atmospheric pollutants at the industrial park site under the constraint of the increase in ambient concentration of atmospheric pollutants, including: using the composite response function integrated in step S4, substituting the emission of different pollutants from different industrial parks under the condition of the upper limit of the concentration increase caused by the increase in pollution emissions from the industrial park in step S3, and solving the concentration contribution of different industrial parks to the specific environmental monitoring point; S6: Verification of the impact of the evaluation site, including: using the composite response function, substituting the maximum allowable emission data of different parks obtained, solving the change of pollutant concentration at the evaluation environment evaluation site, and verifying whether it meets the evaluation environment evaluation site concentration limit pre-set in step S3.

[0020] Specifically, in this exemplary embodiment, the problem of balancing incremental pollution emissions required for the development of industrial parks is solved through the dynamic balance of emission increases and reductions, and a composite response model of the mutual influence of emissions from multiple parks can be established to realize the analysis of the mutual influence between industrial parks.

[0021] Specifically, in steps S1 and S3, an air quality model is used to assess the reduction in environmental concentrations resulting from pollutant emission reductions around the industrial park. The effectiveness of pollution source emission reductions is quantified, and applicable environmental pollutant concentration limits are determined based on the current and target concentrations of environmental pollutants at industrial park sites and control sites. This determines the upper limit for the increase in environmental concentrations caused by increased emissions from the industrial park. Furthermore, in steps S2 and S4, a virtual atmospheric pollution source emission inventory is established for newly added industrial emission plots in the industrial park. Emission sequences are established, and the environmental pollutant concentrations corresponding to different emission sequences are simulated to obtain response functions from emissions to environmental concentrations in different industrial parks. Finally, in steps S5 and S6, the response functions of different industrial parks for different evaluation sites are integrated to obtain a composite response function. The maximum allowable emissions of the industrial park are calculated using the concentration limits of the industrial park monitoring sites, and the maximum allowable emissions of the industrial park are verified and adjusted using the concentration limits of the evaluation sites.

[0022] In summary, this exemplary embodiment proposes a method for calculating the maximum allowable emissions based on the dynamic composite balance of emission increases and decreases in industrial parks, aiming to make full use of the space for pollutant concentration reductions brought about by emission reductions of pollution sources around the industrial park, offset the increase in pollutant concentrations caused by new pollution emissions in the industrial park, and develop the park industry without causing a deterioration in the current atmospheric environment. The core of this method is to establish a dynamic balance relationship between emission increases and decreases of pollution sources inside and outside the park, and through scientific measurement and dynamic regulation, to achieve precise control of the total amount of pollutant emissions in the park, and provide technical support for the green development of the park. Subsequently, while ensuring environmental quality, by vacating high-pollution, low-output emission sources, the maximum value of the pollution emission increment required for the development of enterprises in the park can be calculated, thereby contributing to green and high-quality development.

[0023] The following content will explain each step in detail. For detailed step-by-step process, please refer to Figure 2 : More preferably, in an exemplary embodiment, in step S1, a quantitative assessment of the effectiveness of emission reduction of air pollution sources around the industrial park is achieved. By using regional masks, the buffer zone range around the industrial park in the simulation grid is marked, and the emission reduction coefficients corresponding to different pollution sources are used to reduce pollutant emissions around the industrial park. The air quality model is then used to simulate the emission inventory of air pollution sources before and after emission reduction, and the changes in environmental concentration at different observation points caused by emission reduction of air pollution sources are calculated. Figure 2 “Emission Reduction Effectiveness Evaluation” section in the .

[0024] Specifically, the changes in environmental concentrations caused by emission reductions of air pollution sources around the quantified industrial park include: S11: Use geographic information system (GIS) software to establish industrial park regional buffer zones based on the boundaries of industrial parks, and calculate the local emission reduction potential ( Figure 2The buffer zone radius required for setting the concentration reduction demand is based on the emission reduction potential analysis in the “Emission Reduction Effectiveness Evaluation” section of the Figure 2 Considering that pollution emissions mainly affect nearby areas, the radius of the buffer zone can be set to 2-5 km. One specific implementation method is as follows: Figure 3 As shown; S12: Figure 4 As shown, based on the relationship between the latitude and longitude of the center of the air quality simulation grid and the industrial park buffer zone, the simulation grids in the buffer zone are screened, and the emission reduction mask file of the industrial park buffer zone is established with the industrial park as the unit; S13: Prepare emission inventory adjustment documents based on emission reduction ratios of different industrial parks ( Figure 2 (Compilation of regional emission reduction scenario documents in the "Emission Reduction Effectiveness Evaluation" section of the IEEE International Conference on Climate Change and Climate Change, 2016) Use the emission model combined with the emission inventory reduction ratio to adjust the air quality model emission inventory data ( Figure 2 Emission reduction scenario list in the "Emission Reduction Effectiveness Assessment" section of the ); where the emission reduction coefficients corresponding to different pollution sources are mapped to the emission inventory adjustment coefficients within the mask area (those outside the mask are set to 1); S14: Use the emission reduction before ( Figure 2 The basic emission inventory in the “Emission Reduction Effectiveness Assessment” section of the Figure 2 Air quality simulation (based on the emission reduction inventory in the “Emission Reduction Effectiveness Evaluation” section of the Figure 2 The CMAQ / CAMx and other air quality models in the “Emission Reduction Effectiveness Evaluation” section of the paper were used to obtain the p Industrial parks e Pollutant concentration before emission reduction C pe0 and concentration after emission reduction C pe1 ; Combined with the concentration of environmental monitoring stations in industrial parks C peo , correct the pollutant concentration simulation results, take the relative concentration change, and calculate the pollutant concentration reduction caused by emission reduction according to the following formula ΔC pe ( Figure 2 The emission reduction concentration reduction of the park site in the "Emission Reduction Effectiveness Evaluation" section ): Evaluation sites for industrial parks E Emission reduction concentration, using the concentration before emission reduction at the evaluation site C Ee0 and concentration after emission reduction C Ee1 , monitoring concentration at evaluation sites C Eeo, and then calculate after correction ( Figure 2 Emission reduction and concentration reduction of the evaluation site in the "Emission Reduction Effectiveness Evaluation" section ): .

[0025] Among them, it should be noted that: (1) The mask file is spatially aligned with the air quality emission inventory, and the emission amount of the inventory within the mask range is adjusted according to the emission reduction ratio, thereby simulating the change in environmental concentration caused by emission reduction of atmospheric pollution sources and quantifying the emission reduction amount; (2) The main purpose of this exemplary embodiment is to utilize the space for pollutant concentration reduction brought about by emission reduction of pollution sources around the park to offset the increase in pollutant concentration caused by new pollution emissions from the park, so as to develop the park industry without causing the deterioration of the current atmospheric environment, and use the emission reduction in the surrounding area to provide the reduction amount for the park development. Therefore, a mask file is used for emission reduction, and the increase in emission is mainly caused by industry, which has a specific park boundary. Therefore, the subsequent increase in emission does not require a mask file.

[0026] (3) Figure 3 is the vector diagram of the emission reduction area mask file in the model, which is only the mask range. In steps S13 and 14, the emission inventory within the range is screened according to this range, and the emission inventory within the range is adjusted according to the emission reduction ratio. After adjustment, the emission inventory is input into the model for simulation. The concentration difference before and after emission reduction is obtained by subtracting the simulation results of the inventory before and after adjustment, that is, C peo -C pel .

[0027] (4) The simulated concentration in step S14 is the simulation result, which is a grid result. The simulation result of the grid point is extracted using the longitude and latitude of the station, which is the corresponding station simulation concentration.

[0028] More preferably, in an exemplary embodiment, in step S2, the industrial park newly added emission source response sequence simulation is carried out. In order to realize the serialized simulation of the newly added emission response of the industrial park, it is necessary to establish an emission inventory and complete the air quality simulation through this step. Figure 2 The “Emission Impact Analysis” section of the .

[0029] Specifically, the changes in environmental concentrations caused by the response sequence of newly added emission sources in the simulated industrial park include: S21: Taking the industrial park as the boundary, establish emission inventory data in combination with the land use type structure within the industrial park ( Figure 2In the "Emission Impact Analysis" section of the industrial park boundary gridding and unit emission intensity inventory preparation); in a specific exemplary embodiment, in the absence of detailed planning, considering the average distribution within the industrial land area, an emission intensity of 1t / km is established. 2 ·a unit emission intensity emission inventory data to form a baseline emission inventory.

[0030] S22: Based on historical data, establish a time distribution series for industrial source emissions, including monthly distribution coefficients Fac_m , weekly distribution coefficient Fac_w and hourly distribution coefficient Fac_h ; Establish species partition coefficient Fac_s , to achieve the distribution of complex pollutants to air quality simulation components ( Figure 2 Emission allocation coefficients for industrial parks in the “Emission Impact Analysis” section of the report); S23: A model for processing incremental emissions from industrial parks, used to convert annual emission data from industrial parks under a specific emission intensity into hourly emission data with time variations, and to incorporate incremental emissions from industrial parks into the basic emission inventory data of the air quality model (corresponding to Figure 2 Emission sequence per unit area of ​​industrial parks in the “Emission Impact Analysis” section of the 2nd edition); p The emission intensity per unit area of ​​industrial parks EMIS pa Down, i Month, week j 、 k Hours, components l The emission data is calculated as follows: S24: For different industrial parks, different unit area emissions EMIS pa Establish an emission sequence, use the industrial park emission increase model to operate the basic inventory of the air quality model, add the hourly emission data of each pollutant component newly emitted by the industrial park to the corresponding pollutant component in the basic emission inventory of the air quality model, and obtain the scenario emission inventory data (corresponding to Figure 2 List of sequential scenarios in the “Emission Impact Analysis” section of the Implementation Plan); Call the air quality model to complete the simulation of the basic emission inventory and different scenario emission inventory data (corresponding to Figure 2 The CMAQ / CAMx air quality model in the “Emission Impact Analysis” section of the IEEE Conference on Air Quality (hereinafter referred to as the “Impact of Emissions”) is used as follows: c Emission inventory simulated concentrations for scenarios C pec Simulated concentrations compared to the base case C peb , Computing Industrial Park pNew emission sources e In the industrial park p Pollutant response sequence concentration ΔC pea (correspond Figure 2 Simulation results of the park site in the "Emission Impact Analysis" section ): Evaluation Site E , extract review sites E Corresponding simulation grid scenario emission inventory concentration C Eec , basic scenario simulation concentration C Eeb , the corresponding pollutant response sequence concentration can be calculated ΔC Eea , the method is as follows (corresponding to Figure 2 Simulation results of the park site in the "Emission Impact Analysis" section ): .

[0031] More preferably, in an exemplary embodiment, in step S3, the incremental space of the atmospheric pollutant environmental concentration at the industrial park site is determined to achieve the balance basis of the maximum allowable emission amount. Figure 2 Refer to the “Calculation and Determination of Maximum Allowable Emissions” section.

[0032] More specifically, the determination of the incremental space of atmospheric pollutant environmental concentrations at industrial park sites includes: Combined with the pollutant concentration reduction obtained from the quantitative assessment of the emission reduction effectiveness of air pollution sources around the industrial park in step S1 , using the current concentration of site p in the industrial park C peo0 and target atmospheric pollution concentration values ​​for specific years C peoy Determine the concentration increase limit in the park caused by the increase in pollution emissions in the park , as the basis for calculating the maximum allowable emissions (corresponding to Figure 2 The concentration increase limit of the park in the "Calculation and determination of maximum allowable emissions" section ), calculated as follows: Among them, the current concentration of site p in the industrial park C peo0 Corresponding to the concentration of the industrial park environmental monitoring station in step S1 C peo .

[0033] More preferably, in an exemplary embodiment, the establishment of the industrial park composite response relationship in step S4 is used to fit the simulation results of the serialized grid emission inventories of different industrial parks to obtain the response function relationship between the emission increments and environmental concentrations of different environmental monitoring sites and industrial parks, and different pollutants. At the same time, the response functions of different parks are integrated to obtain the composite response function of the emissions of different industrial parks to the pollutant concentrations of specific environmental monitoring sites. Figure 2 The “Emission Impact Analysis” section of the .

[0034] More specifically, in step S4, establishing the industrial park composite response relationship includes: S41: Combining different EMIS pa Simulation results of the intensity and corresponding industrial park emissions’ contribution to the concentration of the park’s monitoring stations Using cubic natural spline interpolation, we can obtain the response function relationship between emission increment and environmental concentration of different environmental monitoring stations and industrial parks, and different pollutants. C pe =f P (EMIS pa(x) ) , the function fitting is as follows Figure 5 As shown, specifically including: The response function between new emissions from industrial parks and environmental concentrations at park sites f p (EMIS pax ) The establishment of the interpolation process needs to be based on the emission intensity of different unit areas. EMIS pa The corresponding response sequence concentration of the park site obtained from the air quality simulation ΔC pea Divide the simulation into segments and calculate the interpolation point spacing of the nth segment in sequence. For example, the first segment needs to use 0 and EMIS pa(1) 、 Δ C pea(1) , establish a piecewise fitting function: The key to establishing the response function is to solve the a 、 b 、 c 、 d The coefficient of the fitting function is continuous; since the solution is between the value of each segment point and the breakpoint The values ​​are equal, and 2n equations are established, among which EMIS pa(n) are all known, namely: … Since the fitting function is continuous and differentiable, the first-order derivatives at each segment point are equal, and the following n-1 equations are established: … At the same time, the second-order derivatives of the functions at each segment point are equal, and the following n-1 equations are established: … Finally, the second-order derivative of the first and last points is 0, which serves as the boundary condition, namely: think 、 , since in the above equation, EMIS pa(1) and ΔC pea(1) The unit area park emission intensity input into the model and the corresponding park site concentration increment calculated by the model are both known quantities. Therefore, after combining the above equations into a system of equations, the linear equations can be solved to obtain the corresponding values ​​for different segments. a 、 b 、 c 、 d coefficient, we get f p (EMIS pa(x) ) (correspond Figure 2 The park site response function in the "Emission Impact Analysis" section f p (EMIS pa(x) ), and for the functional relationship between industrial park emissions and evaluation sites, it is only necessary to replace Replace with That is, the corresponding functional relationship is obtained f E (EMIS pa(x)) (correspond Figure 2 Evaluation site response function in the "Emission Impact Analysis" section f E (EMIS pa(x) ); S42: Integrate the response functions of different industrial parks to obtain the composite response function of different industrial park emissions to the pollutant concentration of a specific environmental monitoring station; the contribution concentration of P different industrial parks to the pollutant e at a certain industrial park monitoring station C e Calculate according to the following formula and substitute the unit area emissions of different parks EMIS pa(x) The concentration increment at the corresponding station can be solved C e : .

[0035] More preferably, in an exemplary embodiment, the maximum allowable emission under the constraint of the environmental concentration increment of atmospheric pollutants at the industrial park site in step S5 is solved based on the functional relationship established in step S4, and the emission intensity data of different industrial parks can be brought into the solution. EMIS pa(x) To solve, since the pollution emission intensity adjustment involves multiple parks, the concentration increment of the same pollutant C e Often corresponding to multiple park emission combinations, in order to cover the emission increment space of different parks, this exemplary embodiment uses the exhaustive method of the maximum value constraint to solve. Figure 2 See the "Solving the Maximum Allowable Emissions" section in the .

[0036] More specifically, the maximum allowable emissions of atmospheric pollutants at industrial park sites under the constraint of incremental ambient concentrations are solved, including: Based on the functional relationship established in step S4, the emission intensity data of different industrial parks are introduced EMIS pa(x) To solve, due to the pollution emission intensity adjustment involving multiple industrial parks, the same pollutant concentration increment C e It often corresponds to multiple park emission combinations. In order to cover the emission increment space of different parks, the exhaustive method with maximum value constraints is used to solve it. Using the response function establishment method in step S4, the dependent variable and the independent variable are exchanged to obtain the functional relationship between the emission intensity change and the concentration increment, that is, , using this function and the corresponding concentration increase limit in step S3 That is, calculate the maximum emission when only this park emits EMIS paxm ( Figure 2The maximum emission increment of the park in the "Maximum allowable emission solution" section EMIS paxm ), as the upper limit of the park emission increment, and then set the number of times the response equation is solved S ( Figure 2 The number of times the response equation is solved in the "Maximum Allowable Emissions Solution" section S Determine), calculate the incremental step size as follows: From 0 to EMIS paxm By step length Establish an emission increment sequence, orthogonalize it by industrial park, form an input data set, and bring it into the functional relationship established in step S4 for solution calculation ( Figure 2 The response equation in the "maximum allowable emission solution" section is vectorized and solved). The required number of calculations is related to the number of solutions S and the number of parks p, which is If there are three parks and each park performs 500 solutions, then 125 million calculations are required. This invention uses vectorized calculations to improve efficiency. Finally, the maximum allowable emissions of the park are calculated. EMIS pa(x) ( Figure 2 The maximum allowable emissions solution for the park in the "Maximum allowable emissions solution" section EMIS pa(x) ).

[0037] More preferably, in an exemplary embodiment, for the impact verification of the evaluation site in step S6, a composite response function is used to bring in the obtained maximum allowable emission data of different parks, solve the change in pollutant concentration of the evaluation environment evaluation site, and verify whether it meets the pre-set concentration limit of the evaluation environment evaluation site. Figure 2 See the "Solving the Maximum Allowable Emissions" section in the .

[0038] More specifically, the evaluation site impact verification ( Figure 2 The assessment site verification in the "Verification and Determination of Maximum Allowable Emissions" section of the evaluation process includes: The solution obtained in step S5 is analyzed and screened. The result of the solution is taken as the solution that meets the conditions, and the value trend of different parks in the actual park development plan is combined to select the final unit area emission increment of different parks. , using the method described in step S4 to establish a functional relationship between park emissions and evaluation sites f E (EMIS pa(x) ) , combined with the current concentration of the evaluation siteC Eeo0 , target concentration C Eeoy , reduce emissions and concentrations Carry out the concentration impact verification of the evaluation site. It is considered that the increased emissions from the park have not caused the deterioration of the air quality at the evaluation site. Otherwise, the emissions from the park should be reduced to ensure the contribution to the comprehensive concentration of the evaluation site. Below 0 ( Figure 2 The maximum allowable emissions balance results of the park are as follows: .

[0039] Obviously, the above embodiments are merely examples for clarity of explanation and are not intended to limit the implementation methods. Those skilled in the art will readily appreciate that other variations or modifications can be made based on the above descriptions. It is not necessary and impossible to enumerate all implementation methods here. Obvious variations or modifications derived therefrom remain within the scope of protection of the present invention.

Claims

1. The method for calculating the maximum allowable emissions of industrial parks based on the balance of emissions increases and decreases is characterized by: Including the following step: S1: Quantify the changes in ambient concentrations caused by emission reductions from air pollution sources around industrial parks, including: establishing a buffer zone around the industrial park in the simulation grid by setting a mask area; mapping the emission reduction coefficients corresponding to different pollution sources to the emission inventory adjustment coefficients within the mask area; using an air quality model to simulate the emission inventory of air pollution sources before and after emission reductions, and calculating the reduction in ambient concentrations at different monitoring stations caused by emission reductions from air pollution sources; S2: Simulate the changes in environmental concentration caused by the response sequence of new emission sources in industrial parks, including: taking industrial parks as units, gridding the new emissions in industrial parks according to the types of planned land in the industrial parks, and establishing a serialized grid emission inventory that can be applied to air quality model simulation; for serialized grid emission inventories with different emission intensities per unit area, based on the urban air pollution source emission inventory, using the serialized grid emission inventory as the emission increment, and using the air quality model to complete the simulation of the corresponding serialized grid emission inventory, respectively. Calculate the difference between the simulation results of the serialized grid emission inventory for different emission amounts and the simulation results of the urban air pollution source emission inventory, and obtain the corresponding environmental concentration changes under the emission inventory increments of all industrial parks; S3: Determining the incremental space for ambient concentrations of atmospheric pollutants at the industrial park site, including: combining the ambient concentration reductions from the quantitative assessment of emission reduction effectiveness of atmospheric pollution sources surrounding the industrial park calculated in step S1, the current status of the industrial park monitoring sites, and the target atmospheric pollution concentration values ​​for a specific year, to determine the upper limit of concentration increases caused by increased pollution emissions from the industrial park as a basis for calculating the maximum allowable emissions; S4: Establishing a composite response relationship for the industrial park, including: combining the emission intensity per unit area in step S2 with the corresponding environmental concentration changes under the corresponding industrial park emission inventory increment, calculating the response function relationship between the emission increment and environmental concentration of each monitoring site and each industrial park, and the corresponding different pollutants; and integrating the response functions of different industrial parks to obtain a composite response function of the emissions of different industrial parks to the pollutant concentration of a specific monitoring site; S5: solving the maximum allowable emission of atmospheric pollutants at the industrial park site under the constraint of the increase in ambient concentration of atmospheric pollutants, including: using the composite response function integrated in step S4, substituting the emission of different pollutants from different industrial parks under the condition of the upper limit of the concentration increase caused by the increase in pollution emissions from the industrial park in step S3, and solving the concentration contribution of different industrial parks to the specific environmental monitoring point; S6: Verification of the impact of the evaluation site, including: using the composite response function, substituting the maximum allowable emission data of different parks obtained, solving the change of pollutant concentration at the evaluation environment evaluation site, and verifying whether it meets the evaluation environment evaluation site concentration limit pre-set in step S3.

2. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 1 is characterized by: In step S1, the quantification of the environmental concentration changes caused by the emission reduction of air pollution sources around the industrial park includes: S11: Use Geographic Information System (GIS) software to establish an industrial park regional buffer zone based on the industrial park boundary, and set the required buffer zone radius based on the local emission reduction potential and concentration reduction needs; S12: Based on the relationship between the latitude and longitude of the center of the air quality simulation grid and the industrial park buffer zone, the simulation grids within the buffer zone are screened, and an emission reduction mask file for the industrial park buffer zone is established based on the industrial park as a unit; S13: Prepare emission inventory adjustment files based on the emission reduction ratios of different industrial parks. Use the emission model combined with the emission inventory reduction ratios to adjust the air quality model emission inventory data. Align the emission reduction coefficients corresponding to different pollution sources with the emission inventory adjustment coefficients within the mask area. S14: Use the emission inventory before and after emission reduction to carry out air quality simulation and obtain the p Industrial parks e Pollutant concentration before emission reduction C pe0 and concentration after emission reduction C pe1 ; Combined with the concentration of environmental monitoring stations in industrial parks C peo , correct the pollutant concentration simulation results, take the relative concentration change, and calculate the pollutant concentration reduction caused by emission reduction according to the following formula ΔC pe : Evaluation sites for industrial parks E Emission reduction concentration, using the concentration before emission reduction at the evaluation site C Ee0 and concentration after emission reduction C Ee1 , monitoring concentration at evaluation sites C Eeo , and calculate after correction: 。 3. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 2 is characterized by: In step S2, the simulation of the environmental concentration changes caused by the response sequence of the newly added emission sources in the industrial park includes: S21: Using the industrial park as the boundary, establish emission inventory data in combination with the land use type structure within the industrial park; S22: Based on historical data, establish a time distribution series for industrial source emissions, including monthly distribution coefficients Fac_m , weekly distribution coefficient Fac_w and hourly distribution coefficient Fac_h ; Establish species partition coefficient Fac_s , to achieve the distribution of complex pollutants to air quality simulation components; S23: A model for processing incremental emissions from industrial parks, which is used to convert annual emission data from industrial parks under a specific emission intensity into hourly emission data with time variations, and to incorporate incremental emissions from industrial parks into the basic emission inventory data of the air quality model; p The emission intensity per unit area of ​​industrial parks EMIS pa Down, i Month, week j 、 k Hours, components l The emission data is calculated as follows: S24: For different industrial parks, different unit area emissions EMIS pa Establish an emission sequence and use the industrial park emission increase model to operate the basic inventory of the air quality model. Add the hourly emission data of each pollutant component newly emitted by the industrial park to the corresponding pollutant component in the basic emission inventory of the air quality model to obtain the scenario emission inventory data. Call the air quality model to complete the simulation of the basic emission inventory and different scenario emission inventory data, according to the following formula, using c Emission inventory simulated concentrations for scenarios C pec Simulated concentrations compared to the base case C peb , Computing Industrial Park p New emission sources e In the industrial park p Pollutant response sequence concentration ΔC pea : Evaluation Site E , extract review sites E Corresponding simulation grid scenario emission inventory concentration C Eec , basic scenario simulation concentration C Eeb , the corresponding pollutant response sequence concentration can be calculated ΔC Eea , the method is as follows: 。 4. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 3 is characterized by: In step S21, in the absence of detailed planning, the average distribution of industrial land is considered and the emission intensity is established as 1t / km 2 ·a unit emission intensity emission inventory data to form a baseline emission inventory.

5. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 3 is characterized by: In step S3, determining the incremental space of atmospheric pollutant environmental concentration at the industrial park site includes: Combined with the pollutant concentration reduction obtained from the quantitative assessment of the emission reduction effectiveness of air pollution sources around the industrial park in step S1 , using the current concentration of site p in the industrial park C peo0 and target atmospheric pollution concentration values ​​for specific years C peoy Determine the concentration increase limit in the park caused by the increase in pollution emissions in the park , as the basis for calculating the maximum allowable emissions, the calculation formula is as follows: Among them, the current concentration of site p in the industrial park C peo0 Corresponding to the concentration of the industrial park environmental monitoring station in step S1 C peo .

6. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 5 is characterized by: In step S4, the establishment of the industrial park composite response relationship includes: S41: Combining different EMIS pa Simulation results of the intensity and corresponding industrial park emissions’ contribution to the concentration of the park’s monitoring stations Using cubic natural spline interpolation, we can obtain the response function relationship between emission increment and environmental concentration of different environmental monitoring stations and industrial parks, and different pollutants. C pe =f P (EMIS pa(x) ) , specifically including: The response function between new emissions from industrial parks and environmental concentrations at park sites f p (EMIS pax ) The establishment of the interpolation process needs to be based on the emission intensity of different unit areas. EMIS pa The corresponding response sequence concentration of the park site obtained from the air quality simulation Δ C pea Divide the simulation into segments and calculate the interpolation point spacing of the nth segment in sequence. For example, the first segment needs to use 0 and EMIS pa(1) 、 Δ C pea(1) , establish a piecewise fitting function: The key to establishing the response function is to solve the a 、 b 、 c 、 d The coefficient of the fitting function is continuous; since the solution is between the value of each segment point and the breakpoint The values ​​are equal, and 2n equations are established, among which EMIS pa(n) are all known, namely: … Since the fitting function is continuous and differentiable, the first-order derivatives at each segment point are equal, and the following n-1 equations are established: … At the same time, the second-order derivatives of the functions at each segment point are equal, and the following n-1 equations are established: … Finally, the second-order derivative of the first and last points is 0, which serves as the boundary condition, namely: Since in the above equation, EMIS pa(1) and ΔC pea(1) The unit area park emission intensity input into the model and the corresponding park site concentration increment calculated by the model are both known quantities. Therefore, after combining the above equations into a system of equations, the linear equations can be solved to obtain the corresponding values ​​for different segments. a 、 b 、 c 、 d coefficient, we get f p (EMIS pa(x) ) , and for the functional relationship between industrial park emissions and evaluation sites, it is only necessary to replace Replace with That is, the corresponding functional relationship is obtained f E (EMIS pa(x) ) ; S42: Integrate the response functions of different industrial parks to obtain the composite response function of different industrial park emissions to the pollutant concentration of a specific environmental monitoring station; the contribution concentration of P different industrial parks to the pollutant e at a certain industrial park monitoring station C e Calculate according to the following formula and substitute the unit area emissions of different parks EMIS pa(x) The concentration increment at the corresponding station can be solved C e : 。 7. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 6 is characterized by: In step S5, the maximum allowable emission of atmospheric pollutants at the industrial park site under the constraint of the ambient concentration increment is solved, including: Based on the functional relationship established in step S4, the emission intensity data of different industrial parks are introduced EMIS pa(x) To solve, due to the pollution emission intensity adjustment involving multiple industrial parks, the same pollutant concentration increment C e It often corresponds to multiple park emission combinations. In order to cover the emission increment space of different parks, the exhaustive method with maximum value constraints is used to solve it. Using the response function establishment method in step S4, the dependent variable and the independent variable are exchanged to obtain the functional relationship between the emission intensity change and the concentration increment, that is, , using this function and the corresponding concentration increase limit in step S3 That is, calculate the maximum emission when only this park emits EMIS paxm , as the upper limit of the park emission increment, and then set the number of times the response equation is solved S , calculate the incremental step size as follows: From 0 to EMIS paxm By step length Establish an emission increment sequence, orthogonalize it by industrial park, form an input data set, and bring it into the function relationship established in step S4 for solution calculation. The number of calculations required is related to the number of solutions S and the number of parks p, which is Second-rate.

8. The method for calculating the maximum allowable emissions of an industrial park based on the balance of emissions increases and decreases according to claim 7 is characterized by: In step S6, the evaluation site impact verification includes: The solution obtained in step S5 is analyzed and screened. The result of the solution is taken as the solution that meets the conditions, and the value trend of different parks in the actual park development plan is combined to select the final unit area emission increment of different parks. , using the method described in step S4 to establish a functional relationship between park emissions and evaluation sites f E (EMIS pax ) , combined with the current concentration of the evaluation site C Eeo0 , target concentration C Eeoy , reduce emissions and concentrations Carry out the concentration impact verification of the evaluation site. It is considered that the increased emissions from the park have not caused the deterioration of the air quality at the evaluation site. Otherwise, the emissions from the park should be reduced to ensure the contribution to the comprehensive concentration of the evaluation site. Below 0: 。

Images