Pollutant source analysis method and system based on air quality model

Abstract

The invention discloses a pollutant source analysis method and system based on an air quality model, and belongs to the technical field of environment monitoring, and the method comprises the following specific steps: simulating an output result file based on a three-dimensional air quality numerical model CAMx, extracting and calculating a pollution source concentration contribution information unit in a simulation time period, and then displaying a result. Meanwhile, the invention discloses a system based on the method, and by the adoption of the pollutant source analysis method and system based on the air quality model, the model simulation resolution in a receptor area is improved in a double-layer nesting mode; the boundary conveying contribution, the contribution of the pollution source area outside the acceptor area and the mutual contribution of different acceptor points in the acceptor area are comprehensively considered, and a reliable data basis is provided for subsequent pollution source analysis.

Description

Technical Field

[0001] The present invention relates to the technical field of environmental monitoring, and in particular to a pollutant source analysis method based on an air quality model and a system thereof. Background Art

[0002] With the acceleration of urbanization, industrialization and regional economic integration, air pollution is changing from a single type of air pollution (such as coal smoke pollution, motor vehicle pollution, and petrochemical pollution) to complex air pollution, which is mainly manifested in the increase in the concentration of atmospheric oxidative pollutants (such as PM2.5 and ozone), the significant decrease in atmospheric visibility, and the deterioration of the atmospheric environment spreading to the entire region. In recent years, many regions have experienced severe air pollution on a regional scale. For example, the Beijing-Tianjin-Hebei region, the Yangtze River Delta, the Pearl River Delta, and other urban agglomerations have shown obvious regional atmospheric complex pollution characteristics, which seriously restricts the sustainable development of regional social economy and harms the health of the public. However, the relationship between atmospheric oxidative pollutants and precursors is nonlinear, making their prevention and control very complicated and difficult. More importantly, most of the existing analytical methods are for the analysis of a single receptor point, lacking in-depth and precise analysis of the impact of pollution source areas on multiple receptor points and the mutual influence between multiple receptor points in the receptor area. Summary of the invention

[0003] The purpose of the present invention is to provide a pollutant source analysis method and system based on an air quality model to solve the above technical problems.

[0004] To achieve the above object, the present invention provides a method for pollutant source analysis based on an air quality model, and the specific steps are as follows:

[0005] Step S1: based on the simulation output results of the air quality model CAMx, the pollution source concentration information unit of the required receptor point is extracted and stored in the original data format;

[0006] Step S2: Read the data storage file generated in step S1, and reprocess and analyze the pollution source concentration information unit of the data storage file; the specific process is as follows:

[0007] Step S21: extracting the pollution source concentration contribution data of all source areas to the set receptor point, the pollution source concentration contribution data of all industries of other receptor points in the receptor area to the set receptor point, the boundary transport concentration and the background pollution source concentration from all pollution source concentration contribution information units;

[0008] Step S22: The source area is the area within the boundary grid except the receptor area. The source area includes the main source area and the secondary source area. The concentration contribution data of pollution sources in various industries in the main source area to the set receptor point are summed to obtain the total main source area transport volume. The concentration contribution data of pollution sources in corresponding industries in the secondary source area to the set receptor point are summed to obtain the total secondary source area transport volume.

[0009] Step S23: summing up the concentration contribution data of the corresponding industry pollution sources of the receptor points in the receptor area to obtain the local generation amount;

[0010] Step S24: the boundary transport concentration is the transport amount of the area outside the boundary grid;

[0011] Step S3: Display the processed pollution source concentration information in a graph.

[0012] Preferably, in step S1, the simulation result file output by the air quality model CAMx is read through a character stream editor, and the pollution source concentration contribution information unit includes the simulation time, pollution source area, pollutant type, and concentration contributed by the pollution source industry.

[0013] Preferably, the main source area meets the following conditions:

[0014] The distance between the center point of the main source area and the center point of the receptor point is less than or equal to the set influence distance and the pollution degree of the main source area is greater than or equal to the set pollution degree; the set influence distance is related to the wind field at the simulation time; the secondary source area is the area within the boundary grid except the main source area.

[0015] Preferably, the results are displayed in a bar graph or a pie chart and a text file is output simultaneously.

[0016] A system based on the above-mentioned pollutant source analysis method based on an air quality model includes:

[0017] Air quality model CAMx, used to simulate pollution sources based on relevant gas concentrations and meteorological data and output simulation result files;

[0018] An extraction module, used to extract pollution source concentration contribution information units;

[0019] An analysis and processing module is used to analyze the pollution source concentration contribution information unit, classify and count the data in the pollution source concentration contribution information unit, and generate a text file;

[0020] The display module is used to display classification and statistical results as well as text files.

[0021] Therefore, the present invention adopts the above-mentioned pollutant source analysis method and system based on air quality model, which has the beneficial effects of: using air quality model CAMx to obtain simulation results, using double-layer nesting to improve the resolution in the analysis area, and improving simulation accuracy. The analysis and processing module is used to classify and count the data in the pollution source concentration contribution information unit, not only considering the impact of boundary transport, but also considering the impact of the external source area of ​​the study area and the mutual contribution rate of different receptor points in the study area, and dividing the source area into the main source area and the secondary source area, providing reliable pollution source analysis data for subsequent pollution source research.

[0022] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 This is a flow chart of a pollutant source analysis method based on an air quality model of the present invention;

[0024] Figure 2 A schematic diagram of a simulation area according to an embodiment of the present invention;

[0025] Figure 3 This is a schematic diagram of the receptor region of an embodiment of the present invention;

[0026] Figure 4 This is a pie chart of the average component contribution rate of O3 concentration in Liaoning region from May to August 2023 of the present invention;

[0027] Figure 5 This is a bar graph of the average component contribution rate of O3 concentration in Liaoning region from May to August 2023 of the present invention;

[0028] Figure 6 This is a bar graph of the average component contribution rate of O3 concentration in the source area from May to August 2023 of the present invention;

[0029] Figure 7 This is a bar graph of the average component contribution rate of O3 concentration in the receptor area from May to August 2023 of the present invention. DETAILED DESCRIPTION

[0030] In the description of the present invention, it should be noted that the terms "upper", "lower", "inside", "outside", etc. indicate positions or positional relationships based on the positions or positional relationships shown in the accompanying drawings, or the positions or positional relationships in which the invented product is usually placed when in use. They are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In the description of the present invention, it should also be noted that, unless otherwise clearly specified and limited, the terms "setting", "installation", and "connection" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be an indirect connection through an intermediate medium, or it can be a connection between the two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0031] The embodiments of the present invention are described in detail below in conjunction with the accompanying drawings.

[0032] like Figure 2-Figure 3 As shown in the figure, Liaoning is taken as the research object (receptor area). The meteorological field driving of CAMx is provided by WRF (version4.1.2). The model simulation area is double nested, and the center point is (116.5°E, 39.5°N). The horizontal resolution of the first layer is 27km, the number of grids is 119 (east-west) × 119 (north-south), covering most of the central and eastern regions of my country (source area); the horizontal resolution of the second layer is 9km, the number of grids is 136 (east-west) × 103 (north-south), covering Liaoning, Hebei, Shandong and southeastern Inner Mongolia, where the points are the main source areas about Liaoning. The vertical direction is divided into 27 layers, and Lambert projection is selected. The parameterization schemes are: Lin microphysics scheme, RRTM longwave radiation scheme, Goddard shortwave radiation, Grell3D cumulus convection scheme, Noah land surface process scheme, and YSU boundary layer scheme. NCEP GFS data are used as the meteorological initial conditions and boundary conditions of the model. The boundary conditions are updated every 6 hours, and the horizontal resolution of the data is 1°×1°.

[0033] The receptor areas include 14 prefecture-level cities in Liaoning (Shenyang, Dalian, Anshan, Fushun, Benxi, Dandong, Jinzhou, Yingkou, Fuxin, Liaoyang, Panjin, Tieling, Chaoyang, and Huludao), eight major source areas (Jilin, Heilongjiang, Inner Mongolia, Shandong, Hebei, Beijing, Tianjin, and Jiangsu), and “other areas” other than Liaoning and eight surrounding provinces as secondary source areas.

[0034] NO based on the 2019 multi-scale emission inventory model MEIC (Multi-resolution Emission Inventory for China) x , CO, NMVOC, PM 2.5 The emission inventory of CAMx model operation is prepared using ISAT V2.0 tool, and the inventory of various pollutants and air pollutants (resolution of 0.1°) of five sectors, including power, industry, civil, transportation and agriculture. The model simulation period is from 00:00 on April 24, 2023 to 18:00 on August 31, 2023 (UTC), of which the period before 00:00 on May 1 is the model start adjustment period, and the model output frequency is 1h.

[0035] like Figure 1 As shown, a pollutant source analysis method based on an air quality model is characterized by the following specific steps:

[0036] Step S1: Based on the simulation output results of the air quality model CAMx, the pollution source concentration information unit of the required receptor point is extracted and stored in the original data format.

[0037] Step S2: Read the data storage file generated in step S1, and read the simulation result file output by the air quality model WRF-CAMx through the character stream editor. The pollution source concentration contribution information unit includes the simulation time, pollution source area, pollutant type, and the concentration contributed by the pollution source industry. Reprocess and analyze the pollution source concentration information unit of the data storage file; the specific analysis process is as follows:

[0038] Step S21: extracting the pollution source concentration contribution data of all source areas to the set receptor point, the pollution source concentration contribution data of all industries of other receptor points in the receptor area to the set receptor point, the boundary transport concentration and the background pollution source concentration from all pollution source concentration contribution information units;

[0039] Step S22: The source area is the area within the boundary grid except the receptor area. The source area includes the main source area and the secondary source area. The concentration contribution data of pollution sources in various industries in the main source area to the set receptor point are summed to obtain the total main source area transport volume. The concentration contribution data of pollution sources in corresponding industries in the secondary source area to the set receptor point are summed to obtain the total secondary source area transport volume.

[0040] Step S23: summing up the concentration contribution data of the corresponding industry pollution sources of the receptor points in the receptor area to obtain the local generation amount;

[0041] Step S24: the boundary transport concentration is the transport amount of the area outside the boundary grid;

[0042] Step S3: Display the processed pollution source concentration information in a chart. The result is displayed in a bar chart or pie chart, and a text file is output at the same time.

[0043] like Figure 4-Figure 7 As shown in the figure, for Liaoning as a whole, the boundary transport contributes the most to the O3 concentration in Liaoning under the average conditions from May to August 2023, reaching 65.20%, followed by the contribution of locally generated O3 concentration (15.96%) and the contribution of transport from 8 main source areas (14.73%); the contribution of secondary source transport is the smallest, at 4.09%. All areas outside the outer grid contribute the most to Liaoning, accounting for a much larger proportion than other contributing components, and the contribution of local generation and transport from surrounding main source areas cannot be ignored.

[0044] A system based on the above-mentioned pollutant source analysis method based on an air quality model includes an air quality model CAMx, an extraction module, an analysis and processing module, and a display module connected in sequence.

[0045] The air quality model CAMx is used to simulate pollution sources based on relevant gas concentrations and meteorological data and output simulation result files.

[0046] The extraction module is used to extract the pollution source concentration contribution information unit.

[0047] The analysis and processing module is used to analyze the pollution source concentration contribution information unit, classify and count the data in the pollution source concentration contribution information unit, and generate text files.

[0048] The display module is used to display classification and statistical results as well as text files.

[0049] Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that they can still modify or replace the technical solution of the present invention with equivalents, and these modifications or equivalent replacements cannot cause the modified technical solution to deviate from the spirit and scope of the technical solution of the present invention.

Claims

1. A pollutant source analysis method based on an air quality model, characterized in that: The specific steps are as follows: Step S1: extracting the pollution source concentration information unit of the required receptor point based on the simulation output result of the air quality model CAMx, and storing it in the original data format; Step S2: Read the data storage file generated in step S1, and reprocess and analyze the pollution source concentration information unit of the data storage file; the specific process is as follows: Step S21: extracting the pollution source concentration contribution data of all source areas to the set receptor point, the pollution source concentration contribution data of all industries of other receptor points in the receptor area to the set receptor point, the boundary transport concentration and the background pollution source concentration from all pollution source concentration contribution information units; Step S22: The source area is the area within the boundary grid except the receptor area. The source area includes the main source area and the secondary source area. The concentration contribution data of pollution sources in various industries in the main source area to the set receptor point are summed to obtain the total main source area transport volume. The concentration contribution data of pollution sources in corresponding industries in the secondary source area to the set receptor point are summed to obtain the total secondary source area transport volume. Step S23: summing up the concentration contribution data of the corresponding industry pollution sources of the receptor points in the receptor area to obtain the local generation amount; Step S24: the boundary transport concentration is the transport volume of the area outside the boundary grid; Step S3: Display the processed pollution source concentration information in a graph.

2. The method for pollutant source analysis based on an air quality model according to claim 1, characterized in that: In step S1, the simulation result file output by the air quality model CAMx is read through the character stream editor, and the pollution source concentration contribution information unit includes the simulation time, pollution source area, pollutant type, and the concentration contributed by the pollution source industry.

3. The method for pollutant source analysis based on an air quality model according to claim 2, characterized in that: The main source areas meet the following conditions: The distance between the center point of the main source area and the center point of the receptor point is less than or equal to the set influence distance and the pollution degree of the main source area is greater than or equal to the set pollution degree; the set influence distance is related to the wind field at the simulation time; the secondary source area is the area within the boundary grid except the main source area.

4. The method for pollutant source analysis based on an air quality model according to claim 3, characterized in that: The results are displayed in a bar graph or pie chart and output to a text file.

5. A system for pollutant source analysis based on an air quality model according to any one of claims 1 to 4, characterized in that: include: Air quality model CAMx, used to simulate pollution sources based on relevant gas concentrations and meteorological data and output simulation result files; An extraction module, used to extract pollution source concentration contribution information units; An analysis and processing module is used to analyze the pollution source concentration contribution information unit, classify and count the data in the pollution source concentration contribution information unit, and generate a text file; The display module is used to display classification and statistical results as well as text files.

Images