A method for monitoring air quality in a park

Abstract

The present invention discloses a method for monitoring the air quality of a park, which belongs to the field of data analysis technology, and includes the following steps: S1, constructing an air index set and a state index set for the park; S2, determining the flow impact value of the state index set on the air index set; S3, determining whether the air quality of the park is qualified according to the flow impact value of the state index set on the air index set. The present invention takes into account that wind speed is one of the important factors affecting air quality. By monitoring the wind speed within a certain period of time, the air flow impact value in the park is determined, which is helpful for analyzing the diffusion degree of carbon monoxide concentration caused by wind speed; the present invention also performs grid splitting on the park, and performs data analysis technology on the carbon monoxide concentration in the grid, so as to realize comprehensive monitoring of the air quality in the park and improve the accuracy and reliability of the monitoring results.

Description

Technical Field

[0001] The invention belongs to the technical field of data analysis, and in particular relates to an air quality monitoring method for a park. Background Art

[0002] With the acceleration of industrialization, the air quality of the industrial park has become an important factor affecting public health and environmental sustainable development. Due to its special production characteristics, chemical parks emit various pollutants, which have a significant impact on the air quality of the surrounding environment. With the strengthening of environmental protection laws and regulations and the improvement of public environmental awareness, accurately predicting the air quality of chemical parks has become an important task for environmental supervision and corporate environmental management. The air quality of chemical parks is affected by many factors, but traditional air quality monitoring methods have problems such as inaccurate monitoring and failure to consider meteorological factors. Summary of the invention

[0003] In order to solve the above problems, the present invention proposes a method for monitoring the air quality of a park.

[0004] The technical solution of the present invention is: a method for monitoring air quality in a park comprises the following steps:

[0005] S1, construct an air index set and a state index set for the park;

[0006] S2. Determine the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid in the park;

[0007] S3. Determine whether the air quality of the park is qualified based on the flow impact value of the state index set on the air index set.

[0008] Furthermore, in S1, the method for constructing the air index set is: dividing the park into a number of Euler grids, and determining the carbon monoxide concentration corresponding to each Euler grid at the current moment, to generate the air index set;

[0009] In S1, the method for constructing the state index set is: determining the wind speed of each Euler grid in the historical time period and the current moment, and generating the state index set.

[0010] Euler grid is a numerical simulation method widely used in fields such as fluid dynamics, atmospheric science and environmental science. It can divide the calculation area into a series of fixed grid cells.

[0011] Furthermore, S2 includes the following sub-steps:

[0012] S21, taking the wind speed of each Euler grid in the historical time period as the wind speed sequence of each Euler grid, and taking the autocorrelation of the wind speed sequence as the change value of the corresponding Euler grid;

[0013] S22, determining the flow coefficient of each Euler grid according to the safe wind speed of the park and the wind speed of each Euler grid at the current moment;

[0014] S23. Determine the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid.

[0015] The beneficial effect of the above further scheme is: in the present invention, the carbon monoxide concentration is a dynamic variable, and the carbon monoxide concentration in each area of ​​the entire park is not necessarily exactly the same, so it is necessary to grid the park and perform zoning analysis. The wind speed in the historical time period and the current moment in the park will affect the diffusion speed of the carbon monoxide concentration in the park at the current moment. The greater the wind speed, the faster the diffusion speed of carbon monoxide in each grid. Therefore, it is necessary to effectively analyze the impact of wind speed on carbon monoxide concentration at different times, so that the air quality monitoring results are more accurate.

[0016] Safe wind speed refers to the maximum wind speed at which, under specific circumstances, the wind will not cause damage to personnel, equipment or structures.

[0017] Furthermore, in S22, the flow coefficient F of the i-th Euler grid is i The calculation formula is: ; Where V safe Indicates the safe wind speed in the park, v i represents the wind speed of the ith Euler grid at the current moment, M represents the number of Euler grids adjacent to the ith Euler grid, v m represents the wind speed of the mth Euler grid adjacent to the ith Euler grid at the current moment, e represents the index, v1 represents the wind speed of the first Euler grid at the current moment, and v I represents the wind speed of the Ith Euler grid at the current moment, max(·) represents the maximum value function, and I represents the number of Euler grids in the park.

[0018] Further, S23 includes the following sub-steps:

[0019] S231, extracting the maximum value between the standard deviation of all Euler grid change values ​​and the standard deviation of all Euler grid flow coefficients as the traversal initial value;

[0020] S232, using the traversal initial value to traverse each Euler grid in turn, to obtain the current flow traversal value of each Euler grid;

[0021] S233. Determine the flow influence value of the state index set on the air index set according to the current flow traversal values ​​of all Euler grids.

[0022] The beneficial effect of the above further scheme is: in the present invention, the change value and flow coefficient corresponding to the Euler grid can reflect the change of wind speed in the grid within a certain period of time. Therefore, the value with the largest standard deviation is selected as the traversal initial value, and the current flow traversal value is calculated for each grid, thereby obtaining the influence of the entire wind speed value on the entire carbon monoxide concentration value, that is, the flow influence value.

[0023] Further, in S232, the current flow traversal value B of the i-th Euler grid i The calculation formula is: ; In the formula, Start represents the initial value of the traversal, B i-1 represents the current flow traversal value of the i-1th Euler grid, F i represents the flow coefficient of the ith Euler grid, C i represents the change value of the i-th Euler grid, max(·) represents the maximum function, min(·) represents the minimum function, and ω represents the learning rate.

[0024] Furthermore, in S233, the method for determining the flow influence value of the state index set on the air index set is: taking the mean of the current flow traversal values ​​of all Euler grids and rounding them down as the flow influence value of the state index set on the air index set.

[0025] Furthermore, S3 includes the following sub-steps:

[0026] S31, extracting the row and column of the Euler grid corresponding to the highest carbon monoxide concentration as the first reference row and the first reference column;

[0027] S32, extracting the row and column of the Euler grid corresponding to the lowest carbon monoxide concentration as the second reference row and the second reference column;

[0028] S33, determining an air state threshold value according to the first reference row, the first reference column, the second reference row, the second reference column, and the flow influence value of the state index set on the air index set;

[0029] S34. Determine whether the air quality in the park is qualified based on the air status threshold.

[0030] The beneficial effect of the above further scheme is: in the present invention, the rows and columns where the highest and lowest carbon monoxide concentrations are located are extracted, which are more representative to generate thresholds. When the carbon monoxide concentration in the Euler grid is lower than the air state threshold, it means that the carbon monoxide concentration at that location is abnormal and the air quality is unqualified.

[0031] Furthermore, in S33, the calculation formula of the air state threshold η is: ; In the formula, O 11represents the mean carbon monoxide concentration of all Eulerian grids in the first reference row, O 12 represents the mean carbon monoxide concentration of all Eulerian grids in the second reference row, O 21 represents the mean carbon monoxide concentration of all Eulerian grids in the first reference column, O 22 represents the mean carbon monoxide concentration of all Euler grids in the second reference column, and I represents the flow impact value of the state index set on the air index set.

[0032] The beneficial effects of the present invention are as follows: the present invention takes into account that wind speed is one of the important factors affecting air quality, and by monitoring the wind speed within a certain period of time, the air flow impact value in the park is determined, which is helpful for analyzing the degree of diffusion of carbon monoxide concentration caused by wind speed; the present invention also divides the park into grids, and by performing data analysis technology on the carbon monoxide concentration in the grid, comprehensive monitoring of the air quality in the park is achieved, thereby improving the accuracy and reliability of the monitoring results. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] Figure 1 Flowchart of the park's air quality monitoring method. DETAILED DESCRIPTION

[0034] The embodiments of the present invention will be further described below in conjunction with the accompanying drawings.

[0035] like Figure 1 As shown, the present invention provides a method for monitoring air quality in a park, comprising the following steps:

[0036] S1, construct an air index set and a state index set for the park;

[0037] S2. Determine the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid in the park;

[0038] S3. Determine whether the air quality of the park is qualified based on the flow impact value of the state index set on the air index set.

[0039] In the embodiment of the present invention, in S1, the method for constructing the air index set is: dividing the park into a plurality of Euler grids, and determining the carbon monoxide concentration corresponding to each Euler grid at the current moment, and generating the air index set;

[0040] In S1, the method for constructing the state index set is: determining the wind speed of each Euler grid in the historical time period and the current moment, and generating the state index set.

[0041] Euler grid is a numerical simulation method widely used in fields such as fluid dynamics, atmospheric science and environmental science. It can divide the calculation area into a series of fixed grid cells.

[0042] In this embodiment of the present invention, S2 includes the following sub-steps:

[0043] S21, taking the wind speed of each Euler grid in the historical time period as the wind speed sequence of each Euler grid, and taking the autocorrelation of the wind speed sequence as the change value of the corresponding Euler grid;

[0044] S22, determining the flow coefficient of each Euler grid according to the safe wind speed of the park and the wind speed of each Euler grid at the current moment;

[0045] S23. Determine the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid.

[0046] In the present invention, the carbon monoxide concentration is a dynamic variable. The carbon monoxide concentration in each area of ​​the entire park is not necessarily exactly the same. Therefore, the park needs to be gridded and analyzed in different areas. The wind speed in the historical time period and the current moment in the park will affect the diffusion speed of the carbon monoxide concentration in the park at the current moment. The greater the wind speed, the faster the diffusion speed of carbon monoxide in each grid. Therefore, it is necessary to effectively analyze the impact of wind speed on carbon monoxide concentration at different times, so that the air quality monitoring results are more accurate.

[0047] Safe wind speed refers to the maximum wind speed at which, under specific circumstances, the wind will not cause damage to personnel, equipment or structures.

[0048] In the embodiment of the present invention, in S22, the flow coefficient F of the i-th Euler grid i The calculation formula is: ; Where V safe Indicates the safe wind speed in the park, v i represents the wind speed of the ith Euler grid at the current moment, M represents the number of Euler grids adjacent to the ith Euler grid, v m represents the wind speed of the mth Euler grid adjacent to the ith Euler grid at the current moment, e represents the index, v1 represents the wind speed of the first Euler grid at the current moment, and v I represents the wind speed of the Ith Euler grid at the current moment, max(·) represents the maximum value function, and I represents the number of Euler grids in the park.

[0049] In this embodiment of the present invention, S23 includes the following sub-steps:

[0050] Extract the maximum value between the standard deviation of all Euler grid change values ​​and the standard deviation of all Euler grid flow coefficients as the traversal initial value;

[0051] S232, using the traversal S231, extracting the maximum value between the standard deviation of all Euler grid change values ​​and the standard deviation of all Euler grid flow coefficients as the traversal initial value;

[0052] S232, using the traversal initial value to traverse each Euler grid in turn, to obtain the current flow traversal value of each Euler grid;

[0053] S233. Determine the flow influence value of the state index set on the air index set according to the current flow traversal values ​​of all Euler grids.

[0054] In the present invention, the change value and flow coefficient corresponding to the Euler grid can reflect the change of wind speed in the grid within a certain period of time. Therefore, the value with the largest standard deviation is selected as the traversal initial value, and the current flow traversal value is calculated for each grid, thereby obtaining the influence of the entire wind speed value on the entire carbon monoxide concentration value, that is, the flow influence value.

[0055] In the embodiment of the present invention, in S232, the current flow traversal value B of the i-th Euler grid i The calculation formula is: ; In the formula, Start represents the initial value of the traversal, B i-1 represents the current flow traversal value of the i-1th Euler grid, F i represents the flow coefficient of the ith Euler grid, C i represents the change value of the i-th Euler grid, max(·) represents the maximum function, min(·) represents the minimum function, and ω represents the learning rate.

[0056] In an embodiment of the present invention, in S233, the method for determining the flow influence value of the state index set on the air index set is: taking the mean of the current flow traversal values ​​of all Euler grids and then rounding them down as the flow influence value of the state index set on the air index set.

[0057] In this embodiment of the present invention, S3 includes the following sub-steps:

[0058] S31, extracting the row and column of the Euler grid corresponding to the highest carbon monoxide concentration as the first reference row and the first reference column;

[0059] S32, extracting the row and column of the Euler grid corresponding to the lowest carbon monoxide concentration as the second reference row and the second reference column;

[0060] S33, determining an air state threshold value according to the first reference row, the first reference column, the second reference row, the second reference column, and the flow influence value of the state index set on the air index set;

[0061] S34. Determine whether the air quality in the park is qualified based on the air status threshold.

[0062] In the present invention, the rows and columns where the highest and lowest carbon monoxide concentrations are located are extracted, which are more representative, to generate thresholds. When the carbon monoxide concentration in the Euler grid is lower than the air state threshold, it means that the carbon monoxide concentration at that location is abnormal and the air quality is unqualified.

[0063] In the embodiment of the present invention, in S33, the calculation formula of the air state threshold η is: ; In the formula, O 11 represents the mean carbon monoxide concentration of all Eulerian grids in the first reference row, O 12 represents the mean carbon monoxide concentration of all Eulerian grids in the second reference row, O 21 represents the mean carbon monoxide concentration of all Eulerian grids in the first reference column, O 22 represents the mean carbon monoxide concentration of all Euler grids in the second reference column, and I represents the flow impact value of the state index set on the air index set.

[0064] Those skilled in the art will appreciate that the embodiments described herein are intended to help readers understand the principles of the present invention, and should be understood that the protection scope of the present invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific variations and combinations that do not deviate from the essence of the present invention based on the technical revelations disclosed by the present invention, and these variations and combinations are still within the protection scope of the present invention.

Claims

1. A method for monitoring air quality in a park, characterized in that: The following steps are involved: S1, construct an air index set and a state index set for the park; S2. Determine the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid in the park; S3. Determine whether the air quality of the park is qualified according to the flow impact value of the state index set on the air index set; In S1, the method for constructing the air index set is: dividing the park into a number of Euler grids, and determining the carbon monoxide concentration corresponding to each Euler grid at the current moment, to generate the air index set; In S1, the method for constructing the state index set is: determining the wind speed of each Euler grid in the historical time period and the current moment, and generating the state index set; The S2 comprises the following sub-steps: S21, taking the wind speed of each Euler grid in the historical time period as the wind speed sequence of each Euler grid, and taking the autocorrelation of the wind speed sequence as the change value of the corresponding Euler grid; S22, determining the flow coefficient of each Euler grid according to the safe wind speed of the park and the wind speed of each Euler grid at the current moment; S23, determining the flow influence value of the state index set on the air index set according to the change value and flow coefficient of each Euler grid; In S22, the flow coefficient F of the i-th Euler grid i The calculation formula is: ; Where V safe Indicates the safe wind speed in the park, v i represents the wind speed of the ith Euler grid at the current moment, M represents the number of Euler grids adjacent to the ith Euler grid, v m represents the wind speed of the mth Euler grid adjacent to the ith Euler grid at the current moment, e represents the index, v1 represents the wind speed of the first Euler grid at the current moment, and v I represents the wind speed of the Ith Euler grid at the current moment, max(·) represents the maximum value function, and I represents the number of Euler grids in the park; The S23 comprises the following sub-steps: S231, extracting the maximum value between the standard deviation of all Euler grid change values ​​and the standard deviation of all Euler grid flow coefficients as the traversal initial value; S232, using the traversal initial value to traverse each Euler grid in turn, to obtain the current flow traversal value of each Euler grid; S233, determining the flow influence value of the state index set on the air index set according to the current flow traversal values ​​of all Euler grids; In S232, the current flow traversal value B of the i-th Euler grid i The calculation formula is: ; In the formula, Start represents the initial value of the traversal, B i-1 represents the current flow traversal value of the i-1th Euler grid, F i represents the flow coefficient of the ith Euler grid, C i represents the change value of the i-th Euler grid, max(·) represents the maximum value function, min(·) represents the minimum value function, and ω represents the learning rate; In the S233, the method for determining the flow influence value of the state index set on the air index set is: taking the mean of the current flow traversal values ​​of all Euler grids and performing a rounding operation downward as the flow influence value of the state index set on the air index set; The S3 comprises the following sub-steps: S31, extracting the row and column of the Euler grid corresponding to the highest carbon monoxide concentration as the first reference row and the first reference column; S32, extracting the row and column of the Euler grid corresponding to the lowest carbon monoxide concentration as the second reference row and the second reference column; S33, determining an air state threshold value according to the first reference row, the first reference column, the second reference row, the second reference column, and the flow influence value of the state index set on the air index set; S34. Judging whether the air quality in the park is qualified according to the air state threshold; In S33, the calculation formula of the air state threshold η is: ; In the formula, O 11 represents the mean carbon monoxide concentration of all Eulerian grids in the first reference row, O 12 represents the mean carbon monoxide concentration of all Eulerian grids in the second reference row, O 21 represents the mean carbon monoxide concentration of all Eulerian grids in the first reference column, O 22 represents the mean carbon monoxide concentration of all Euler grids in the second reference column, and I represents the flow impact value of the state index set on the air index set.

Images