A method and system for online source analysis of atmospheric pollutants

Abstract

The present invention discloses an online source analysis method and system for atmospheric pollutants, and relates to the technical field of pollutant source analysis. Its main technical solution is as follows: inputting measured meteorological data into a concentration loss prediction model to obtain concentration loss; processing the concentration loss and the measured concentration change to obtain the theoretical concentration change; inputting the emission factor and the theoretical concentration change into a chemical mass balance model to obtain the contribution of each pollution source. The influence of meteorological conditions on the concentration loss is taken into account, in order to achieve the purpose of more accurately identifying and quantifying the contribution of different pollution sources. The influence of current meteorological conditions on the concentration loss can be predicted in real time, in order to achieve the purpose of improving the accuracy of assessing the contribution of pollution sources. The loss caused by the generation of secondary pollutants, as well as diffusion dilution, wet deposition during rain / snow and other reasons is taken into account, in order to effectively ensure the accuracy of the theoretical concentration change, and thus improve the accuracy of analyzing the contribution of pollution sources.

Description

Technical Field

[0001] The present invention relates to the technical field of pollutant source analysis, and in particular to an online source analysis method and system for atmospheric pollutants. Background Art

[0002] With the acceleration of industrialization and urbanization, air pollution is becoming increasingly serious and has become one of the major factors affecting human health and environmental quality. To effectively control and manage air pollution, it is crucial to accurately identify and quantify the contribution of different pollution sources to air pollutant concentrations.

[0003] After the primary pollutants (also known as primary pollutants or primary pollutants) produced by pollution sources enter the external environment, they change under the influence of physical, chemical factors or organisms, or react with other substances in the environment to form secondary pollutants with physical and chemical properties different from those of the primary pollutants. The dynamic process of transformation of primary pollutants into secondary pollutants greatly increases the difficulty of accurately assessing the contribution of each pollution source using existing traceability technology. In addition, changes in meteorological conditions further exacerbate the complexity of the transformation process of primary pollutants into secondary pollutants, which may lead to increased errors in the determination of the contribution of pollution sources. Therefore, how to improve the accuracy of the judgment of the contribution of pollution sources. Therefore, how to improve the accuracy of judging the contribution of pollution sources is a technical problem that needs to be solved urgently. Summary of the Invention

[0004] The purpose of the present invention is to provide an online source analysis method and system for atmospheric pollutants, so as to improve the accuracy of determining the contribution of pollution sources.

[0005] In order to solve the above technical problems, the present invention adopts the following technical solutions:

[0006] In a first aspect, a method for online source apportionment of atmospheric pollutants is provided, comprising the following operations:

[0007] Construct a model to predict the concentration loss of primary pollutants as the concentration loss changes with meteorological changes;

[0008] Collect measured meteorological data and measured concentration changes of primary pollutants in the target area;

[0009] Inputting the measured meteorological data into a concentration loss prediction model to obtain the concentration loss of primary pollutants;

[0010] The concentration loss and measured concentration change of the primary pollutant are processed to obtain the theoretical concentration change of the primary pollutant;

[0011] Obtain emission factors of primary pollutants generated by each pollution source in the target area;

[0012] The emission factors and theoretical concentration changes are input into a chemical mass balance model to obtain the contribution of each pollution source in the target area.

[0013] A further solution is that the process of constructing the concentration loss prediction model includes:

[0014] Obtaining historical information of the sample area; the historical information includes historical meteorological data, and historical concentration changes of primary pollutants corresponding to the historical meteorological data and historical concentration generation of primary pollutants generated by each pollution source;

[0015] Processing the historical concentration change and the historical concentration generation to obtain the historical concentration loss of the primary pollutant;

[0016] The historical meteorological data and the historical concentration loss are trained to obtain the concentration loss prediction model.

[0017] A further approach is to perform pollution source cleaning before inputting the emission factors into the chemical mass balance model. The pollution source cleaning process includes:

[0018] Marking the initial pollutants in the target area as a type of pollutant;

[0019] Label the primary pollutants generated by pollution sources in the target area as Class II pollutants;

[0020] Determine whether the secondary pollutants of the pollution source have a corresponding type of pollutant; if so, retain the emission factor of the pollution source; if not, delete the emission factor of the pollution source.

[0021] A further solution is: the process of obtaining the pollution source contribution includes:

[0022] Inputting the emission factors and theoretical concentration changes into a chemical mass balance model to obtain a contribution share combination of pollution sources;

[0023] Determine whether the number of contribution share combinations is greater than 1; if the number of contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of contribution share combinations is greater than 1, clean the contribution share combinations until the number of contribution share combinations is 1.

[0024] A further solution is: the process of cleaning the contribution share portfolio includes:

[0025] Determine whether any contribution share in the contribution share combination is greater than or equal to 0. If so, retain the contribution share combination; if not, delete the contribution share combination.

[0026] A further solution is that the process of cleaning the contribution share combination also includes:

[0027] After deleting the contribution share combinations with contribution shares less than 0, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, determine whether any contribution share in the contribution share combination is an integer, and if so, retain the contribution share combination; if not, delete the contribution share combination.

[0028] A further solution is that the process of cleaning the contribution share combination also includes:

[0029] After deleting the contribution share combinations whose contribution shares are non-integer, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, calculate the pollution coefficient of the contribution share combination, and process the contribution share of each pollution source and the total contribution share of each pollution source in the contribution share combination with the largest pollution coefficient to obtain the contribution degree of the pollution source.

[0030] A further solution is that the calculation process of the pollution coefficient includes:

[0031] Obtain the pollution risk value of each contribution share and each primary pollutant in the contribution share combination;

[0032] The contribution shares and pollution risk values ​​are calculated to obtain the pollution coefficient of each contribution share combination.

[0033] In the second aspect, an online source analysis system for atmospheric pollutants is provided, which is used to implement the method described in the first aspect. The system includes a construction module, an acquisition module, a first processing module, a second processing module, an acquisition module and a third processing module. The construction module is used to construct a concentration loss prediction model for primary pollutants according to meteorological changes; the acquisition module is used to collect measured meteorological data and measured concentration changes of primary pollutants in the target area; the first processing module is used to input the measured meteorological data into the concentration loss prediction model to obtain the concentration loss of primary pollutants; the second processing module is used to process the concentration loss of primary pollutants and the measured concentration changes to obtain theoretical concentration changes of primary pollutants; the acquisition module is used to obtain the emission factors of primary pollutants generated by each pollution source in the target area; the third processing module is used to input the emission factors and theoretical concentration changes into a chemical mass balance model to obtain the contribution of each pollution source in the target area.

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

[0035] The impact of meteorological conditions on the loss of primary pollutant concentrations is taken into account in order to more accurately identify and quantify the contributions of different pollution sources. The impact of current meteorological conditions on the loss of primary pollutant conversion concentrations can be predicted in real time in order to further improve the accuracy of assessing the contribution of pollution sources. Directly predicting the loss of primary pollutant concentrations under different meteorological conditions not only takes into account the loss when primary pollutants generate secondary pollutants, but also takes into account the loss of primary pollutants caused by other reasons such as diffusion dilution and wet deposition during rain / snow, in order to effectively ensure the accuracy of the theoretical concentration change of primary pollution, thereby improving the accuracy of analyzing the contribution of pollution sources. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 Schematic diagram of a flow chart of an online source analysis method for atmospheric pollutants in this embodiment;

[0037] Figure 2 Schematic diagram of the electrical block diagram of an online source analysis system for atmospheric pollutants in this embodiment. DETAILED DESCRIPTION

[0038] The present invention will be further described below with reference to the accompanying drawings.

[0039] Example 1: This example provides an online source analysis method for atmospheric pollutants, such as Figure 1 As shown, the following operations are included:

[0040] S100. Constructing a model for predicting the concentration loss of primary pollutants as they change with meteorological conditions;

[0041] In this embodiment, the process of constructing the concentration loss prediction model includes:

[0042] S101. Obtain historical information of the sample area; the historical information includes historical meteorological data, and the historical concentration change of the primary pollutant corresponding to the historical meteorological data and the historical concentration generation of the primary pollutant generated by each pollution source;

[0043] For example, during the implementation process, n areas are randomly selected as sample areas, and the n sample areas are marked as 、 、 … Obtain the historical information of each sample area from the historical data and mark the historical information of each sample area accordingly. 、 、 … That is to say, For the sample area Historical information, For the sample area historical information, For the sample area Historical information... For the sample area historical information.

[0044] Historical information includes historical meteorological data, as well as the historical changes in primary pollutant concentrations corresponding to the historical meteorological data and the historical concentrations of primary pollutants generated by each pollution source. Meteorological data includes temperature, humidity, wind speed, etc.

[0045] For example, the sample area Historical information Contains historical weather data , historical concentration changes and historical concentration generation . Historical concentration change For historical meteorological data Next, the sample area The historical concentration change of primary pollutants in the For historical meteorological data Next, the sample area The historical concentration of primary pollutants generated by each pollution source in the sample area. There are m pollution sources in the system, and the historical concentration generation of each pollution source is recorded as 、 、 … The historical concentration generated = + + …+ .in, For the sample area The historical concentration of the first pollutant produced by the first pollution source, For the sample area The historical concentration of primary pollutants generated by the second pollution source, For the sample area The historical concentration of primary pollutants generated by the third pollution source... For the sample area The historical concentration of primary pollutants generated by the mth pollution source.

[0046] Similarly, the sample area Historical information Contains historical weather data , historical concentration changes and historical concentration generation If the sample area There are x pollution sources in the system, then the historical concentration generated is = + + …+ .

[0047] Sample area Historical information Contains historical weather data , historical concentration changes and historical concentration generation If the sample area There are y pollution sources in the system, so the historical concentration generated is = + + …+ .

[0048] Sample area Historical information Contains historical weather data , historical concentration changes and historical concentration generation If the sample area There are l pollution sources in the system, then the historical concentration generated is = + + …+ .

[0049] S102. Processing the historical concentration change and the historical concentration generation to obtain the historical concentration loss of the initial pollutant;

[0050] For example, during implementation, because primary pollutants are converted into secondary pollutants, the historical concentration generation amount will be greater than the historical concentration change amount. Therefore, ideally, the difference between the historical concentration generation amount and the historical concentration change amount can be considered the historical concentration loss of the primary pollutant. The historical concentration loss of the primary pollutant is the historical concentration loss of the primary pollutant during the process of converting the primary pollutant into a secondary pollutant under historical meteorological data.

[0051] For example, the sample area The historical concentration loss is recorded as ,and = - ; Sample area The historical concentration loss is recorded as ,and = - ; Sample area The historical concentration loss is recorded as ,and = - ...sample area The historical concentration loss is recorded as ,and = - .

[0052] S103. Training the historical meteorological data and historical concentration loss to obtain the concentration loss prediction model.

[0053] For example, during the implementation process, the historical meteorological data and historical concentration loss are cleaned to form a data set ={( , ),( , ),( , )……( , )}. Use linear regression model to analyze the data set The concentration loss prediction model can be obtained by training with the historical meteorological data and historical concentration loss.

[0054] S200. Collecting measured meteorological data and measured concentration changes of primary pollutants in the target area;

[0055] For example, during implementation, the area where the initial pollutant online source analysis is to be conducted is marked as a target area. Current meteorological data in the target area is collected and recorded as measured meteorological data. Simultaneously, the current measured concentration change of the initial pollutant is collected.

[0056] S300. Inputting the measured meteorological data into the concentration loss prediction model to obtain the initial pollutant concentration loss;

[0057] For example, during implementation, measured meteorological data is input into the concentration loss prediction model to obtain the primary pollutant concentration loss. This primary pollutant concentration loss is the primary pollutant concentration loss in the target area under the current meteorological conditions. This primary pollutant concentration loss is the predicted value.

[0058] S400. The concentration loss of the initial pollutant and the measured concentration change are processed to obtain the theoretical concentration change of the initial pollutant;

[0059] For example, during implementation, under ideal conditions, the primary pollutant concentration generated by each pollution source in the target area is equal to the sum of the primary pollutant concentration loss and the measured concentration change. Therefore, the theoretical primary pollutant concentration change can be calculated by adding the primary pollutant concentration loss to the measured primary pollutant concentration change. This theoretical primary pollutant concentration change is the primary pollutant concentration generated by each pollution source in the target area under the current meteorological conditions.

[0060] S500. Obtain the emission factors of primary pollutants generated by each pollution source in the target area;

[0061] For example, during implementation, the emission factors of primary pollutants generated by each pollution source in the target area can be obtained from historical emission inventories, industry standards, and literature.

[0062] S600. Input the emission factor and theoretical concentration change into a chemical mass balance model to obtain the contribution of each pollution source in the target area.

[0063] For example, during the implementation process, the emission factors of primary pollutants generated by each pollution source in the target area and the theoretical concentration change of primary pollutants in the target area are input into a chemical mass balance model to obtain the contribution of each pollution source in the target area. On the one hand, the online source apportionment method for atmospheric pollutants in this embodiment considers the impact of meteorological conditions on the loss of primary pollutant concentration, in order to more accurately identify and quantify the contribution of different pollution sources. On the other hand, the online source apportionment method for atmospheric pollutants in this embodiment can predict in real time the impact of current meteorological conditions on the loss of primary pollutant concentration due to conversion, in order to further improve the accuracy of assessing pollution source contributions. On the other hand, the online source apportionment method for atmospheric pollutants in this embodiment directly predicts the concentration loss of primary pollutants under different meteorological conditions, taking into account not only the loss of primary pollutants when they are converted into secondary pollutants, but also the loss of primary pollutants due to other factors such as diffusion dilution and wet deposition during rain / snow, in order to effectively ensure the accuracy of the theoretical concentration change of primary pollutants and thereby improve the accuracy of the estimated pollution source contributions.

[0064] During the actual processing process, certain primary pollutants may only be produced by specific pollution sources. For example, coal-fired power plants may produce characteristic primary pollutants containing mercury (Hg), and oil refineries may produce characteristic primary pollutants containing hydrogen sulfide (H2S). Therefore, in order to clean specific pollution sources, reduce the data processing workload of the chemical mass balance model, and improve the efficiency of analyzing the contribution of each pollution source, in this embodiment, before the emission factors are input into the chemical mass balance model, S000. Pollution source cleaning processing is also performed. The process of cleaning the pollution sources includes:

[0065] S001. Mark the initial pollutants in the target area as Class I pollutants;

[0066] For example, during implementation, it is assumed that the primary pollutants present in the target area include primary pollutant A, primary pollutant B, primary pollutant C, and primary pollutant D, and primary pollutants A, B, C, and D are labeled as a class of pollutants. That is, a class of pollutants M = (A, B, C, D).

[0067] S002. Label the primary pollutants generated by pollution sources in the target area as Class II pollutants;

[0068] For example, during the implementation process, it is assumed that there are 4 pollution sources in the target area, and the first pollution source produces primary pollutants A, B and C. The primary pollutants produced by the first pollution source are recorded as Class II pollutants. ,and = (A, B, C). The second pollution source produces primary pollutants A, B, and D. The primary pollutants produced by the second pollution source are recorded as Class II pollutants. ,and = (A, B, D). The third pollution source produces primary pollutants A, B, C, and D. The primary pollutants produced by the third pollution source are recorded as Class II pollutants. ,and = (A, B, C, D). The fourth pollution source produces primary pollutants A, B, C, and E. The primary pollutants produced by the fourth pollution source are recorded as Class II pollutants. ,and =(A, B, C, E).

[0069] S003. Determine whether the secondary pollutants of the pollution source have a corresponding type of pollutant; if so, retain the emission factor of the pollution source; if not, delete the emission factor of the pollution source.

[0070] For example, during the implementation process, it is determined whether the secondary pollutants of the pollution source are evenly distributed in the corresponding type 1 pollutants. For example, the second type of pollutants generated by the first pollution source = (A, B, C) where A, B, and C all exist in the first category of pollutants. In M = (A, B, C, D), it is considered that the first pollution source in the target area may be producing primary pollutants, so the emission factor of the first pollution source is retained. The second category of pollutants produced by the second pollution source = (A, B, D) where A, B, and D all exist in the same type of pollutant M = (A, B, C, D). If the second pollution source in the target area is considered to be producing primary pollutants, the emission factor of the second pollution source is retained. = (A, B, C, D) A, B, C, D all exist in the first category of pollutants M = (A, B, C, D), and it is considered that the third pollution source in the target area may be producing primary pollutants, so the emission factor of the third pollution source is retained. The second category of pollutants produced by the fourth pollution source = (A, B, C, E) where E does not exist in a class of pollutants M = (A, B, C, D), it is assumed that the fourth pollution source in the target area does not produce primary pollutants, and the emission factor of the fourth pollution source is deleted. By judging whether the primary pollutants in the target area are produced by certain specific pollution sources, it is determined whether there are specific pollution sources in the target area or whether the characteristic pollution sources in the target area are working, and the emission factors of these characteristic pollution sources are deleted or retained. The goal is to reduce the number of emission factors input into the chemical mass balance model, thereby reducing the data processing volume of the chemical mass balance model and improving the efficiency of analyzing the contribution of each pollution source.

[0071] In this embodiment, the process of obtaining the pollution source contribution includes:

[0072] S601. Input the emission factor and theoretical concentration change into the chemical mass balance model to obtain the contribution share combination of the pollution source;

[0073] For example, during the implementation, for example, the primary pollutants present in the target area include primary pollutant A, primary pollutant B, primary pollutant C, and primary pollutant D, and the theoretical concentration change of primary pollutant A in the target area is recorded as The theoretical concentration change of the primary pollutant B is recorded as The theoretical concentration change of the primary pollutant C is recorded as The theoretical concentration change of the primary pollutant D is recorded as .

[0074] The first pollution source in the target area produces primary pollutants A, B and C, and the emission factor of the first pollution source producing primary pollutant A is recorded as , the emission factor of the primary pollutant B produced by the first pollution source is recorded as , the emission factor of the primary pollutant C produced by the first pollution source is recorded as The second pollution source produces primary pollutants A, B and D, and the emission factor of primary pollutant A produced by the second pollution source is recorded as , the emission factor of the primary pollutant B produced by the second pollution source is recorded as , the emission factor of the primary pollutant D produced by the second pollution source is recorded as The third pollution source produces primary pollutants A, B, C and D, and the emission factor of primary pollutant A produced by the third pollution source is recorded as , the emission factor of the primary pollutant B produced by the third pollution source is recorded as , the emission factor of the primary pollutant C produced by the third pollution source is recorded as , the emission factor of the primary pollutant D produced by the third pollution source is recorded as As shown in Table 1:

[0075] Table 1

[0076]

[0077] The theoretical concentration change is recorded as 、 、 and , and emission factors 、 、 、 、 、 、 、 、 and , input into the chemical mass balance model, we can get the relationship shown in formula (1):

[0078] (1),

[0079] in, The contribution share of the first pollution source in the target area; The contribution share of the second pollution source in the target area; This is the contribution share of the third pollution source in the target area.

[0080] Solve the above relationship 、 and , we can get the contribution share of each pollution source in the target area, and get the combination formed by the contribution share of each pollution source in the target area, that is, the contribution share combination ( , , ).

[0081] S602. Determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, clean the contribution share combinations until the number of the contribution share combinations is 1.

[0082] For example, during implementation, when the number of primary pollutant types present in the target area is greater than or equal to the number of pollution sources in the target area, a single contribution share combination will be obtained. When the number of primary pollutant types present in the target area is less than the number of pollution sources in the target area, multiple contribution share combinations may exist.

[0083] Therefore, in order to obtain a more accurate contribution share combination of each pollutant, after obtaining the contribution share combination of each pollution source in the target area, the number of contribution share combinations is also judged to determine whether the number of contribution share combinations is greater than 1.

[0084] When the number of contribution share combinations is 1, the contribution of the pollution source in the contribution share combination is divided by the total contribution share of each pollution source to obtain the contribution degree of the pollution source. For example, the relationship shown in formula (1) can be used to obtain 1 contribution share combination, which is recorded as =( , , ), then the contribution of the first pollution source in the target area is ,and ; The contribution of the second pollution source in the target area is ,and ; The contribution of the third pollution source in the target area is ,and .

[0085] When the number of contribution share combinations is greater than 1, the contribution share combinations are cleaned until the number of contribution share combinations is 1. For example, there are 4 primary pollutants in the target area, including primary pollutant A, primary pollutant B, primary pollutant C and primary pollutant D. There are 5 pollution sources in the target area, and the first pollution source in the target area produces primary pollutant A, primary pollutant B and primary pollutant C, and the emission factor of the first pollution source producing primary pollutant A is recorded as , the emission factor of the primary pollutant B produced by the first pollution source is recorded as , the emission factor of the primary pollutant C produced by the first pollution source is recorded as The second pollution source produces primary pollutants A, B and D, and the emission factor of primary pollutant A produced by the second pollution source is recorded as , the emission factor of the primary pollutant B produced by the second pollution source is recorded as , the emission factor of the primary pollutant D produced by the second pollution source is recorded as The third pollution source produces primary pollutants A, B, C and D, and the emission factor of primary pollutant A produced by the third pollution source is recorded as , the emission factor of the primary pollutant B produced by the third pollution source is recorded as , the emission factor of the primary pollutant C produced by the third pollution source is recorded as , the emission factor of the primary pollutant D produced by the third pollution source is recorded as The fourth pollution source produces primary pollutants A, C, and D, and the emission factor of primary pollutant A produced by the fourth pollution source is recorded as , the emission factor of the primary pollutant C produced by the fourth pollution source is recorded as , the emission factor of the primary pollutant D produced by the fourth pollution source is recorded as The fifth pollution source produces primary pollutants A and D, and the emission factor of primary pollutant A produced by the fifth pollution source is recorded as , the emission factor of the primary pollutant D produced by the fifth pollution source is recorded as As shown in Table 2:

[0086] Table 2

[0087]

[0088] The theoretical concentration change is recorded as 、 、 and , and emission factors 、 、 、 、 、 、 、 、 、 、 、 、 、 and , input into the chemical mass balance model, we can get the relationship shown in formula (2):

[0089] (2),

[0090] The contribution share of the first pollution source in the target area; The contribution share of the second pollution source in the target area; The contribution share of the third pollution source in the target area; The contribution share of the fourth pollution source in the target area; This is the contribution share of the fifth pollution source in the target area.

[0091] According to the relationship shown in formula (2), we can get contribution share combinations, of which >1, and is an integer. The contribution share combinations are marked as 、 、 … ,in, =( , , , , ); =( , , , , ); =( , , , , )…… =( , , , , ). Contribution share portfolio 、 、 … Cleaning is carried out until the There are only 、 、 … Any one of .

[0092] In this embodiment, the process of cleaning the contribution share combination includes:

[0093] S6021. Determine whether any contribution share in the contribution share combination is greater than or equal to 0. If so, retain the contribution share combination; if not, delete the contribution share combination.

[0094] For example, in the implementation process, under ideal conditions, the contribution share of each pollution source does not have a negative number. That is, when at least one contribution share in a contribution share combination is less than 0, the contribution share combination is considered unreasonable and is deleted. For example, when the contribution share combination Contribution share in , , , , When both are greater than or equal to 0, the contribution share combination is considered Reasonable and retain contribution share combination ; When the contribution share combination Contribution share in , , , , If at least one contribution share is less than 0, the contribution share combination is considered Unreasonable, and delete the contribution share combination .

[0095] Traverse all contribution share combinations 、 、 … , if the contribution share combination 、 … If both are unreasonable, then use the contribution share combination The contribution of the pollution source in the target area is divided by the total contribution of each pollution source to obtain the contribution of the pollution source. For example, the contribution of the first pollution source in the target area is ,and ; The contribution of the second pollution source in the target area is ,and ; The contribution of the third pollution source in the target area is ,and ; The contribution of the fourth pollution source in the target area is ,and ; The contribution of the fifth pollution source in the target area is ,and .

[0096] If the contribution share combination … If all are unreasonable, delete the contribution share combination … , and retain the contribution share 、 、 and At this point, the number of contribution share combinations is still greater than 1. Therefore, in order to further obtain a more accurate contribution share combination of each pollutant, in this embodiment, the process of cleaning the contribution share combination further includes:

[0097] S6022. After deleting the contribution share combinations with contribution shares less than 0, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, determine whether any contribution share in the contribution share combination is an integer, if so, retain the contribution share combination; if not, delete the contribution share combination.

[0098] For example, in the implementation process, under ideal conditions, the contribution share of each pollution source is an integer. That is, when at least one contribution share in a contribution share combination is not an integer, the contribution share combination is considered unreasonable and is deleted. For example, when the contribution share combination Contribution share in , , , , If both are integers, the contribution share combination is considered Reasonable and retain contribution share combination ; When the contribution share combination Contribution share in , , , , If at least one contribution share is not an integer, the contribution share combination is considered Unreasonable, and delete the contribution share combination .

[0099] Traverse all contribution share combinations 、 、 and , if the contribution share combination 、 and If both are unreasonable, then use the contribution share combination The contribution share of the pollution source in the output is divided by the total contribution share of all pollution sources to obtain the contribution degree of the pollution source.

[0100] If the contribution share combination and If all are unreasonable, delete the contribution share combination and , and retain the contribution share and At this point, the number of contribution share combinations is still greater than 1. Therefore, in order to further obtain a more accurate contribution share combination of each pollutant, in this embodiment, the process of cleaning the contribution share combination further includes:

[0101] S6023. After deleting the contribution share combinations with non-integer contribution shares, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, calculate the pollution coefficient of the contribution share combination, and process the contribution share of each pollution source and the total contribution share of each pollution source in the contribution share combination with the largest pollution coefficient to obtain the contribution degree of the pollution source.

[0102] In this embodiment, the calculation process of the pollution coefficient includes:

[0103] Obtain the pollution risk value of each contribution share and each primary pollutant in the contribution share combination;

[0104] For example, during the implementation process, the pollution risk value of each primary pollutant can be obtained from historical emission inventories, industry standards, and literature. The pollution risk value of the primary pollutant can be understood as the potential harm that the primary pollutant may cause to human health or the ecosystem. For example, the pollution risk value of the primary pollutant A obtained is recorded as , the pollution risk value of the initial pollutant B is recorded as , the pollution risk value of the initial pollutant C is recorded as , the pollution risk value of the initial pollutant D is recorded as .

[0105] The contribution shares and pollution risk values ​​are calculated to obtain the pollution coefficient of each contribution share combination.

[0106] For example, during implementation, the contribution share combination Contribution share in , , , and , and the pollution risk value of the initial pollutant A is recorded as The pollution risk value of the primary pollutant B is recorded as The pollution risk value of the primary pollutant C is recorded as The pollution risk value of the initial pollutant D is recorded as , perform calculations and obtain the contribution share combination The pollution coefficient .

[0107] = + + + .

[0108] Contribution share portfolio Contribution share in , , , and , and the pollution risk value of the initial pollutant A is recorded as The pollution risk value of the primary pollutant B is recorded as The pollution risk value of the primary pollutant C is recorded as The pollution risk value of the initial pollutant D is recorded as , perform calculations and obtain the contribution share combination The pollution coefficient .

[0109] = + + + .

[0110] Contribution share portfolio The pollution coefficient , contribution share combination The pollution coefficient Sort by and select the contribution share combination with the largest pollution coefficient for contribution calculation. For example, if the contribution share combination The pollution coefficient Greater than contribution share combination The pollution coefficient , then use the contribution share combination The contribution share of the pollution source is divided by the total contribution share of each pollution source to obtain the contribution degree of the pollution source. The pollution coefficient Less than contribution share combination The pollution coefficient , then use the contribution share combination The contribution share of the pollution source is divided by the total contribution share of each pollution source to obtain the contribution degree of the pollution source. The pollution coefficient Equal contribution share combination The pollution coefficient , then use the contribution share combination ( ) is divided by the total contribution share of all pollution sources to obtain the contribution degree of the pollution source.

[0111] Example 2: This example provides an online source analysis system for atmospheric pollutants, which is used to implement the method described in Example 1. Figure 2As shown, the system includes a construction module, a collection module, a first processing module, a second processing module, an acquisition module and a third processing module.

[0112] Among them, the construction module is used to construct a concentration loss prediction model for primary pollutants as the concentration loss changes with meteorological changes; the acquisition module is used to collect measured meteorological data and measured concentration changes of primary pollutants in the target area; the first processing module is used to input the measured meteorological data into the concentration loss prediction model to obtain the concentration loss of primary pollutants; the second processing module is used to process the concentration loss of primary pollutants and the measured concentration changes to obtain the theoretical concentration changes of primary pollutants; the acquisition module is used to obtain the emission factors of primary pollutants generated by each pollution source in the target area; the third processing module is used to input the emission factors and theoretical concentration changes into the chemical mass balance model to obtain the contribution of each pollution source in the target area.

[0113] Although the present invention has been described herein with reference to a number of illustrative embodiments thereof, it will be understood that numerous other modifications and implementations may be devised by those skilled in the art that fall within the scope and spirit of the principles disclosed herein. More specifically, within the scope of the present disclosure, the drawings, and the claims, numerous variations and modifications may be made to the components and / or layout of the subject combination arrangement. In addition to variations and modifications to the components and / or layout, other uses will also be apparent to those skilled in the art.

Claims

1. A method for online source analysis of atmospheric pollutants, characterized in that: The following operations are included: Construct a model to predict the concentration loss of primary pollutants as the concentration loss changes with meteorological changes; Collect measured meteorological data and measured concentration changes of primary pollutants in the target area; Inputting the measured meteorological data into a concentration loss prediction model to obtain the concentration loss of primary pollutants; The concentration loss and measured concentration change of the primary pollutant are processed to obtain the theoretical concentration change of the primary pollutant; Obtain emission factors of primary pollutants generated by each pollution source in the target area; Inputting the emission factors and theoretical concentration changes into a chemical mass balance model to obtain the contribution of each pollution source in the target area; The process of obtaining the contribution of the pollution source includes: Inputting the emission factors and theoretical concentration changes into a chemical mass balance model to obtain a contribution share combination of pollution sources; Determine whether the number of contribution share combinations is greater than 1; if the number of contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of contribution share combinations is greater than 1, clean the contribution share combinations until the number of contribution share combinations is 1.

2. The method according to claim 1, characterized in that The process of constructing the concentration loss prediction model includes: Obtaining historical information of the sample area; the historical information includes historical meteorological data, and historical concentration changes of primary pollutants corresponding to the historical meteorological data and historical concentration generation of primary pollutants generated by each pollution source; Processing the historical concentration change and the historical concentration generation to obtain the historical concentration loss of the primary pollutant; The historical meteorological data and the historical concentration loss are trained to obtain the concentration loss prediction model.

3. The method according to claim 1, characterized in that Before inputting the emission factors into the chemical mass balance model, a pollution source cleaning process is also performed. The process of cleaning the pollution source includes: Marking the initial pollutants in the target area as a type of pollutant; Label the primary pollutants generated by pollution sources in the target area as Class II pollutants; Determine whether the secondary pollutants of the pollution source have a corresponding type of pollutant; if so, retain the emission factor of the pollution source; if not, delete the emission factor of the pollution source.

4. The method according to claim 1, wherein The process of cleaning the contribution share portfolio includes: Determine whether any contribution share in the contribution share combination is greater than or equal to 0. If so, retain the contribution share combination; if not, delete the contribution share combination.

5. The method according to claim 4, characterized in that The process of cleaning the contribution share portfolio also includes: After deleting the contribution share combinations with contribution shares less than 0, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, determine whether any contribution share in the contribution share combination is an integer, and if so, retain the contribution share combination; if not, delete the contribution share combination.

6. The method according to claim 5, characterized in that The process of cleaning the contribution share portfolio also includes: After deleting the contribution share combinations whose contribution shares are non-integer, determine whether the number of the contribution share combinations is greater than 1; if the number of the contribution share combinations is 1, process the contribution share of each pollution source in the contribution share combination and the total contribution share of each pollution source to obtain the contribution degree of the pollution source; if the number of the contribution share combinations is greater than 1, calculate the pollution coefficient of the contribution share combination, and process the contribution share of each pollution source and the total contribution share of each pollution source in the contribution share combination with the largest pollution coefficient to obtain the contribution degree of the pollution source.

7. The method according to claim 6, characterized in that The calculation process of the pollution coefficient includes: Obtain the pollution risk value of each contribution share and each primary pollutant in the contribution share combination; The contribution shares and pollution risk values ​​are calculated to obtain the pollution coefficient of each contribution share combination.

8. An online source analysis system for atmospheric pollutants, characterized in that: The system is used to implement the method according to any one of claims 1 to 7, and the system includes: A construction module, wherein the construction module is used to construct a concentration loss prediction model of the primary pollutant concentration loss as the concentration loss changes with the meteorological changes; A collection module, the collection module is used to collect measured meteorological data and measured concentration changes of primary pollutants in the target area; A first processing module, the first processing module is used to input the measured meteorological data into a concentration loss prediction model to obtain a concentration loss of a primary pollutant; a second processing module, the second processing module being used to process the concentration loss and the measured concentration change of the primary pollutant to obtain a theoretical concentration change of the primary pollutant; An acquisition module, the acquisition module is used to obtain emission factors of primary pollutants generated by each pollution source in the target area; The third processing module is used to input the emission factor and the theoretical concentration change into a chemical mass balance model to obtain the contribution of each pollution source in the target area.

Images