A method for evaluating the effect of air pollution control at monitoring sites

Abstract

The present invention discloses an air pollution control effect evaluation method for monitoring sites, which belongs to the field of air pollution prevention and control technology, including: obtaining pollutant concentration data and meteorological data of a target site in a target period and a reference period, and dividing wind direction intervals in the two periods respectively; based on the obtained pollutant concentration data and the divided wind direction intervals, calculating the moment-by-moment pollutant contribution concentration and the mean of the pollutant local contribution concentration in each wind direction interval of the two periods; based on the mean of the pollutant local contribution concentration, calculating the difference between the mean of the pollutant local contribution concentration in the reference period and the target period in each wind direction interval; calculating the ratio of the difference in each wind direction interval to the mean of the pollutant local contribution concentration in the reference period, defining the ratio as a relative rate of change, and characterizing the pollution control effect of different positions of the site. The present invention can realize the rapid evaluation of the air pollution control effect in the city, which is of great significance to the improvement of the air pollution situation.

Description

Technical Field

[0001] The present invention belongs to the technical field of air pollution prevention and control, and specifically relates to an air pollution control effect evaluation method for monitoring stations. Background Art

[0002] With the deepening understanding of air pollution problems, small-scale refined control at the site level has become the focus of current air pollution control. Evaluating the effect of air pollution control at monitoring sites is an important part of refined control. Simple and rapid evaluation can achieve rapid response to air pollution control, obtain control effects in a timely manner and make corresponding improvements, and promote further improvement of ambient air quality. Therefore, when implementing refined control of air pollution at monitoring sites, it is particularly necessary to establish a set of air pollution control effect evaluation methods for monitoring sites.

[0003] At present, the air pollution control effect evaluation method mainly relies on air quality models, but its calculation process is complicated and it takes a long time to obtain the results. In addition, the accuracy of the evaluation results may be affected by some unknown pollution mechanisms and the accuracy of the pollution source emission inventory. At the same time, the requirements for resource allocation are also high, and the practical application is not convenient and fast enough. On the other hand, it is difficult to obtain pollutant emission inventories within a small scale, which makes the model-based evaluation method difficult to implement within a small scale. However, by building atmospheric environment monitoring stations, it is possible to visualize the factors that affect air quality, such as pollutant concentrations and meteorological conditions. This measure makes the evaluation of air pollution control effects within a small scale more convenient and efficient. Summary of the invention

[0004] In view of the shortcomings in the prior art, the present invention provides an air pollution control effect evaluation method for monitoring stations, with the aim of realizing a rapid evaluation of the urban air pollution control effect, thereby providing feedback on existing problems to relevant departments, providing support for improving and perfecting control measures, and further strengthening the control of polluted weather, which is of great significance to the improvement of air pollution conditions.

[0005] The present invention discloses an air pollution control effect evaluation method for monitoring stations to achieve rapid evaluation of air pollution control effects, including:

[0006] Step 1: Obtain pollutant concentration data and meteorological data of the target site during the target period and the reference period, and divide the wind direction intervals in the two periods based on the meteorological data; wherein the pollutant concentration data includes PM 2.5 、PM 10 、SO 2 、NO 2 , CO and O 3concentration data of the meteorological data, wherein the meteorological data includes wind direction data;

[0007] Step 2: Based on the pollutant concentration data obtained in step 1 and the divided wind direction intervals, calculate the moment-by-moment local contribution concentration of pollutants and the average local contribution concentration of pollutants in each wind direction interval in the two time periods;

[0008] Step 3: Based on the mean values ​​of local contribution concentrations of pollutants in each wind direction interval of the two time periods obtained in step 2, the difference between the mean values ​​of local contribution concentrations of pollutants in the reference time period and the target time period is calculated in each wind direction interval;

[0009] Step 4: Based on the difference obtained in step 3, calculate the ratio of the difference in each wind direction interval to the average local contribution concentration of pollutants in the reference period, and define the ratio as the relative rate of change to characterize the pollution control effect at different locations of the site.

[0010] As a further improvement of the present invention, in step 1,

[0011] The method for selecting the benchmark period is: using statistical similarity analysis to analyze the similarity of meteorological conditions between the target period and the historical period, and selecting the historical period that meets the requirements as the benchmark period;

[0012] Pollutant concentration data and meteorological data of target sites during the target period and the baseline period can be obtained through data sharing platforms or the corresponding websites of the Meteorological Bureau and the Environmental Protection Bureau;

[0013] The method for dividing wind direction intervals is: based on the wind direction division standards, the wind direction intervals are divided in combination with the actual conditions around the target site. It can be divided into eight wind direction intervals: north, northeast, east, southeast, south, southwest, west, and northwest.

[0014] As a further improvement of the present invention, in step 2, the calculation method of the moment-by-moment local contribution concentration of pollutants and the mean local contribution concentration of pollutants in each wind direction interval in two time periods is:

[0015] According to the calculation method of local and foreign contribution of pollutants, the local contribution ratio of pollutants at each moment in each wind direction interval in two periods of the target station is calculated respectively (C i ), and multiply it by the pollutant concentration (μ) at the same time in the corresponding interval to obtain the local contribution concentration (μ) of pollutants in each interval of the two periods. t ), respectively count the number of local contribution concentrations of pollutants in each wind direction interval in the two periods and calculate the sum (n), and obtain the mean local contribution concentration of pollutants in each wind direction in the two periods The calculation formula is as follows:

[0016] μ t =μ·C i

[0017]

[0018] Where:

[0019] μ t —Local contribution concentration of pollutants in each wind direction interval at each moment;

[0020] μ—the monitoring concentration value of pollutants at each wind direction interval at each moment;

[0021] C i —The local contribution of pollutants in each wind direction interval at each moment;

[0022] —The average local contribution concentration of pollutants in each wind direction interval;

[0023] n—the number of local contribution concentrations of pollutants in each wind direction range.

[0024] As a further improvement of the present invention, in step 3, the difference is calculated as follows:

[0025] The average local contribution concentration of pollutants in each wind direction interval during the benchmark period obtained in step 2 And the mean local contribution concentration of pollutants in each wind direction interval during the target period The difference between the mean of the reference period and the target period in the same wind direction interval is calculated to obtain the difference Δμ, which is calculated as follows:

[0026]

[0027] Δμ—the difference between the mean local contribution concentration of pollutants in the reference period and the target period in each wind direction interval;

[0028] —The average local contribution concentration of pollutants during the reference period;

[0029] —Average local contribution concentration of pollutants during the target period.

[0030] As a further improvement of the present invention, in step 4, the calculation method of the pollution control effect is:

[0031] Based on the difference Δμ between the mean local contribution concentration of pollutants in two periods of each wind direction interval obtained in step 3, calculate the difference Δμ and the mean local contribution concentration of pollutants in the corresponding wind direction interval benchmark period. The ratio η is defined as the relative change rate, which is used to characterize the pollution control effect of each direction of the station during the target period. The calculation formula is as follows:

[0032]

[0033] Where:

[0034] η—Pollution control effect of each direction of the station during the target period;

[0035] Δμ—the difference between the mean local contribution concentration of pollutants in the reference period and the target period of each wind direction;

[0036] —The average value of local contribution concentration during the reference period of each wind direction.

[0037] Compared with the prior art, the present invention has the following beneficial effects:

[0038] The present invention takes into account meteorological conditions and regional transmission, analyzes the differences in site monitoring data at different time periods, and obtains the air pollution control effect simply and quickly, providing theoretical support for environmental-related departments to formulate pollution control measures, which is of great significance to the prevention and control of urban air pollution. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 A flowchart of an air pollution control effect evaluation method for monitoring sites disclosed in an embodiment of the present invention;

[0040] Figure 2 It is the NO between June 2022 and June 2021 in Shilibao, Linyi City, Shandong Province 2 Radar chart of pollution control effectiveness. DETAILED DESCRIPTION

[0041] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0042] The present invention is further described in detail below in conjunction with the accompanying drawings:

[0043] like Figure 1 As shown, the present invention provides an air pollution control effect evaluation method for monitoring stations, comprising:

[0044] S1. Obtain pollutant concentration data and meteorological data of the target site during the target period and the reference period, and divide the wind direction intervals in the two periods based on the meteorological data; wherein,

[0045] The pollutant concentration data include PM 2.5 、PM 10 、SO 2 、NO2 , CO and O 3 concentration data of the meteorological data, wherein the meteorological data includes wind direction data;

[0046] The method for selecting the reference period is: using the statistical similarity analysis method to analyze the similarity of the meteorological conditions of the target period and the historical period, and selecting the historical period that meets the requirements as the reference period. For example, the same month in two years can be considered as two months with similar meteorological conditions; the period is preferably the month, and the time is preferably the hour;

[0047] Pollutant concentration data and meteorological data of target sites during the target period and the baseline period can be obtained through data sharing platforms or the corresponding websites of the Meteorological Bureau and the Environmental Protection Bureau;

[0048] The method for dividing wind direction intervals is: based on the wind direction division standards, the wind direction intervals are divided in combination with the actual conditions around the target site. It can be divided into eight wind direction intervals: north, northeast, east, southeast, south, southwest, west, and northwest.

[0049] S2. Based on the pollutant concentration data obtained in S1 and the divided wind direction intervals, calculate the local contribution concentration of pollutants and the average of the local contribution concentration of pollutants at each moment in each wind direction interval of the two time periods; wherein, the calculation method of the local contribution concentration of pollutants and the average of the local contribution concentration of pollutants at each moment in each wind direction interval of the two time periods is:

[0050] According to the calculation method of local and foreign contribution of pollutants, the local contribution ratio of pollutants at each moment in each wind direction interval in two periods of the target station is calculated respectively (C i ), and multiply it by the pollutant concentration (μ) at the same time in the corresponding interval to obtain the local contribution concentration (μ) of pollutants in each interval of the two periods. t ), respectively count the number of local contribution concentrations of pollutants in each wind direction interval in the two periods and calculate the sum (n), and obtain the mean local contribution concentration of pollutants in each wind direction in the two periods The calculation formula is as follows:

[0051] μ t =μ·C i

[0052]

[0053] Where:

[0054] μ t —Local contribution concentration of pollutants in each wind direction interval at each moment;

[0055] μ—the monitoring concentration value of pollutants at each wind direction interval at each moment;

[0056] C i—The local contribution of pollutants in each wind direction interval at each moment;

[0057] —The average local contribution concentration of pollutants in each wind direction interval;

[0058] n—the number of local contribution concentrations of pollutants in each wind direction range.

[0059] S3, the average local contribution concentration of pollutants in each wind direction interval during the benchmark period obtained based on S2 And the mean local contribution concentration of pollutants in each wind direction interval during the target period The difference between the mean of the reference period and the target period in the same wind direction interval is calculated to obtain the difference Δμ, which is calculated as follows:

[0060]

[0061] Δμ—the difference between the mean local contribution concentration of pollutants in the reference period and the target period in each wind direction interval;

[0062] —The average local contribution concentration of pollutants during the reference period;

[0063] —Average local contribution concentration of pollutants during the target period.

[0064] S4, based on the difference Δμ between the mean local contribution concentration of pollutants in two time periods in each wind direction interval obtained in S3, calculate the difference Δμ and the mean local contribution concentration of pollutants in the corresponding wind direction interval benchmark period The ratio η is defined as the relative change rate, which is used to characterize the pollution control effect of each direction of the station during the target period. The calculation formula is as follows:

[0065]

[0066] Where:

[0067] η—Pollution control effect of each direction of the station during the target period;

[0068] Δμ—the difference between the mean local contribution concentration of pollutants in the reference period and the target period of each wind direction;

[0069] —The average value of local contribution concentration during the reference period of each wind direction.

[0070] Example:

[0071] The present invention provides an air pollution control effect evaluation method for monitoring sites, with Shilibao, Linyi City, Shandong Province as the research site, and the research time is June 2021 and June 2022. June 2021 is used as the reference period, June 2022 is used as the target period, and NO2 As a research pollutant, the pollution control measures of the embodiment are evaluated, including:

[0072] S1. Obtain NO in June 2021 and June 2022 through the national control station in Luozhuang District, Linyi City, Shandong Province. 2 Concentration data and wind direction data; Based on the obtained wind direction data and the wind direction division standards in Table 1, the two time periods are divided into eight wind direction intervals: north, northeast, east, southeast, south, southwest, west and northwest.

[0073] Table 1: Wind direction reference standards

[0074] wind direction Standard(°) North Wind 0-22.5&337.5-360 Northeast Wind 22.5-67.5 Dongfeng 67.5-112.5 Southeast Wind 112.5-157.5 south wind 157.5-202.5 Southwest Wind 202.5-247.5 West Wind 247.5-292.5 Northwest Wind 292.5-337.5

[0075] S2, NO of two time periods obtained based on S1 2 Hourly station concentration data, calculate the hourly NO in each wind direction range in Shilibao Residential Committee, Linyi City in 2021 and 2022 2 Local emissions share (C i ), and the NO at the corresponding time in the corresponding wind direction interval 2 The hourly NO concentrations in each wind direction range in Shilibao, Linyi City in June 2021 and June 2022 were obtained by multiplying the monitoring concentration values ​​(μ) at the site. 2 Local contribution concentration (μ t ), and calculate NO in eight wind directions at the same time 2 Local contribution concentration mean The calculation formula is as follows, and the specific data is shown in Table 2;

[0076] μ t =μ·C i

[0077]

[0078] Where:

[0079] μ t —Hour by hour NO 2 Local contribution concentration;

[0080] μ—NO per hour 2 Site monitoring concentration values;

[0081] C i —Hour by hour NO 2 Percentage of local contribution;

[0082] —Two time periods and eight wind direction intervals NO 2 Mean local contribution concentration;

[0083] n—NO in each wind direction interval2 The amount of local contribution concentration;

[0084] S3, hourly NO in each wind direction range in Shilibao Residential Committee, Linyi City, in June 2021 and June 2022, calculated based on S2 2 The local contribution concentration average is calculated for the corresponding wind direction interval reference period and target period NO 2 The difference between the mean values ​​of local contribution concentrations is given by the following formula, and the detailed data is shown in Table 2;

[0085]

[0086] Where:

[0087] Δμ—Wind directions in June 2021 and June 2022 NO 2 The difference between the mean values ​​of the local contribution concentrations;

[0088] —June 2022NO 2 Mean local contribution concentration;

[0089] —June 2021NO 2 Mean local contribution concentration;

[0090] Table 2: NO in Shilibao, Linyi City, June 2022 compared with June 2021 2 All-round control effect

[0091]

[0092] S4, NO in eight wind direction intervals in Shilibao, Linyi City in June 2021 and June 2022 obtained based on S3 2 The difference between the local contribution concentration mean and the NO in the corresponding wind direction range in June 2021 is calculated. 2 The ratio of the local contribution concentration mean is used to obtain the NO concentration in each wind direction. 2 The relative rate of change is NO 2 The control effect is calculated as follows, and the results are shown in Table 2 and Figure 2 .

[0093]

[0094] Where:

[0095] η—June 2022NO 2 Pollution control effects in all directions;

[0096] Δμ—June 2021 and June 2022NO 2 The difference between the mean values ​​of local contribution concentration in each wind direction interval;

[0097] —June 2021NO 2 Mean local contribution concentration.

[0098] The advantages of the present invention are:

[0099] The present invention takes into account meteorological conditions and regional transmission, analyzes the differences in site monitoring data at different time periods, and obtains the air pollution control effect simply and quickly, providing theoretical support for environmental-related departments to formulate pollution control measures, which is of great significance to the prevention and control of urban air pollution.

[0100] The above are only preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for evaluating the effect of air pollution control at a monitoring station, characterized in that: include: Step 1: Obtain pollutant concentration data and meteorological data of the target site during the target period and the reference period, and divide the wind direction intervals in the two periods based on the meteorological data; wherein the pollutant concentration data includes PM 2.5 、PM 10 , SO2, NO2, CO and O3 concentration data, and the meteorological data includes wind direction data; Step 2: Based on the pollutant concentration data obtained in step 1 and the divided wind direction intervals, calculate the moment-by-moment local contribution concentration of pollutants and the average value of the local contribution concentration of pollutants in each wind direction interval of the two time periods; wherein, in step 2, the calculation method of the moment-by-moment local contribution concentration of pollutants and the average value of the local contribution concentration of pollutants in each wind direction interval of the two time periods is: According to the calculation method of local and foreign contribution of pollutants, the local contribution ratio C of pollutants at each moment in each wind direction interval in two periods of the target station is calculated respectively. i , multiplied by the pollutant station monitoring concentration value μ at the same time in the corresponding wind direction interval, to obtain the local contribution concentration μ of pollutants at each wind direction interval in the two periods t , respectively count the number of local contribution concentrations of pollutants in each wind direction interval in the two periods and calculate the sum n, and obtain the mean local contribution concentration of pollutants in each wind direction interval in the two periods The calculation formula is as follows: m t =μ·C i Step 3: Based on the mean values ​​of the local contribution concentrations of pollutants in each wind direction interval of the two time periods obtained in step 2, the difference between the mean values ​​of the local contribution concentrations of pollutants in the reference time period and the target time period is calculated in each wind direction interval; wherein, in step 3, the difference is calculated as follows: The average local contribution concentration of pollutants in each wind direction interval during the benchmark period obtained in step 2 And the mean local contribution concentration of pollutants in each wind direction interval during the target period The difference between the mean of the reference period and the target period in the same wind direction interval is obtained to obtain the difference Δμ; the calculation formula is: Step 4: Based on the difference obtained in step 3, calculate the ratio of the difference in each wind direction interval to the average local contribution concentration of pollutants in the reference period, and define the ratio as the relative rate of change to characterize the pollution control effect at different locations of the site.

2. The air pollution control effect evaluation method for monitoring stations according to claim 1 is characterized in that: In the step 1, The method for selecting the benchmark period is: using statistical similarity analysis to analyze the similarity of meteorological conditions between the target period and the historical period, and selecting the historical period that meets the requirements as the benchmark period; Channels for obtaining pollutant concentration data and meteorological data of target sites during target and baseline periods include: data sharing platforms and websites of the Meteorological Bureau and Environmental Protection Bureau; The divided wind direction ranges include eight wind direction ranges: north, northeast, east, southeast, south, southwest, west, and northwest.

3. The air pollution control effect evaluation method for monitoring sites according to claim 1 is characterized in that: In step 4, the calculation method of the pollution control effect is: Based on the difference Δμ between the mean local contribution concentration of pollutants in two periods of each wind direction interval obtained in step 3, calculate the difference Δμ and the mean local contribution concentration of pollutants in the corresponding wind direction interval benchmark period. The ratio η is defined as the relative change rate to characterize the pollution control effect of each direction of the station during the target period.

Images