Method and system for analyzing source of one-time artificially discharged VOCs (Volatile Organic Compounds) based on receptor model

Abstract

The invention relates to the technical field of atmosphere volatile organic compound source analysis, in particular to a source analysis method and system for one-time manually discharged VOCs based on a receptor model. Comprising the following steps: constructing a photochemical age parameterized model by using [OH] exposure, initial isoprene emission concentration, [OVOCs] observation concentration, tracer observation concentration, actual [OVOCs] reaction rate constant and tracer reaction rate constant, calculating [OVOCs] primary artificial emission theoretical concentration according to the photochemical age parameterized model, and observing [NMHCs] artificial emission observation concentration, the method comprises the following steps of: obtaining a full-spectrum [VOCs] primary artificial emission observation concentration, carrying out photochemical loss correction and atmospheric diffusion correction on the full-spectrum [VOCs] primary artificial emission observation concentration to obtain a [VOCs] primary artificial emission initial concentration, and carrying out VOCs source analysis by utilizing the [VOCs] primary artificial emission initial concentration according to a PMF receptor model. According to the method, the analysis precision of the PMF receptor model on the VOCs source in the actual atmosphere can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of atmospheric volatile organic compound source analysis, and in particular to a method and system for analyzing the source of primary anthropogenically emitted VOCs based on a receptor model. Background Art

[0002] The Positive Matrix Factorization (PMF) receptor model is a very effective analytical method for evaluating the sources of VOCs.

[0003] The premise of the PMF receptor model to analyze the sources of atmospheric volatile organic compounds (VOCs) is that no photochemical reaction occurs from the emission source to the receptor monitoring point. However, in the actual atmospheric environment, when volatile organic compounds are transmitted from the emission source through the atmosphere to the receptor monitoring point, they will be lost to varying degrees and affected by atmospheric turbulence diffusion, which leads to deviations in the analysis results of the sources of VOCs in the actual atmosphere using the PMF receptor model. Summary of the invention

[0004] The present invention provides a method and system for analyzing the source of VOCs emitted by anthropogenic sources based on a receptor model, the main purpose of which is to improve the accuracy of the PMF receptor model in analyzing the source of VOCs in the actual atmosphere.

[0005] To achieve the above object, the present invention provides a method for analyzing the source of VOCs emitted by anthropogenic sources based on a receptor model, comprising: Obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the reaction rate constant of the tracer, and calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; A photochemical age parameterization model is constructed using the [OH] exposure, the initial isoprene emission concentration, the [OVOCs] observed concentration, the tracer observed concentration, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; The theoretical concentration of anthropogenic emissions of [OVOCs] is calculated according to the photochemical age parameterized model, and the observed concentration of anthropogenic emissions of [NMHCs] is observed using a pre-built observation instrument to obtain the observed concentration of anthropogenic emissions of the full spectrum [VOCs], wherein the observed concentration of anthropogenic emissions of the full spectrum [VOCs] refers to: the theoretical concentration of anthropogenic emissions of [OVOCs] and the observed concentration of anthropogenic emissions of [NMHCs]; The observed concentration of the full spectrum [VOCs] primary anthropogenic emission is corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of the [VOCs] primary anthropogenic emission; According to the pre-built PMF receptor model, the initial concentration of the [VOCs] anthropogenic emission was used to analyze the VOCs source Optionally, the calculation of [OH] exposure, initial isoprene emission concentration, and actual [OVOCs] reaction rate constant includes: Obtaining an emission concentration ratio and an observed concentration ratio of a benzene series species pair, wherein two benzene series species in the benzene series species pair have the same source and different chemical activities; Calculating [OH] exposure using the emission concentration ratio and the observed concentration ratio according to a pre-constructed [OH] exposure formula; The observed concentration of isoprene was obtained, and the initial emission concentration of isoprene was calculated using the following formula according to the [OH] exposure and the observed concentration of isoprene:

[0006] in, represents the initial emission concentration of isoprene, express The observed isoprene concentration at time is the isoprene reaction rate constant, represents [OH] exposure; The photolysis rate of [OVOCs] was calculated using a pre-built photolysis rate formula; Based on the [OVOCs] photolysis rate, the actual [OVOCs] reaction rate constant was calculated using a pre-constructed [OVOCs] reaction rate formula.

[0007] Optionally, the [OH] exposure formula is as follows:

[0008] in, represents the exposure amount of [OH], represents the reaction rate constant of the benzene species with species C and OH radical, represents the reaction rate constant of species B in the benzene series with OH radicals, represents the natural logarithm symbol, Indicates that species C and species B are The emission concentration ratio at the time, Indicates that species C and species B are The observed concentration ratio at time.

[0009] Optionally, the photolysis rate formula is as follows:

[0010] in, represents the photolysis rate of [OVOCs], represents the photolysis rate of [OVOCs] under clear sky conditions, express The photolysis rate under clear sky conditions, Indicates the measured Photolysis rate.

[0011] Optionally, the [OVOCs] reaction rate formula is as follows:

[0012] in, represents the actual [OVOCs] reaction rate constant, represents the reference [OVOCs] reaction rate constant, represents the photolysis rate of [OVOCs], Represents the OH radical concentration.

[0013] Optionally, the photochemical age parameterization model is as follows:

[0014]

[0015]

[0016] in, represents the observed concentration of [OVOCs], represents the observed concentration of the tracer, represents the emission ratio variable of [OVOCs] relative to the tracer, is the tracer reaction rate constant, represents the observed concentration of the tracer, represents the emission ratio of [OVOCs] precursors to tracers, represents the reaction rate variable of [OVOCs] precursor and OH radical, represents the emission ratio variable of [OVOCs] relative to the initial emission concentration of isoprene, Indicates the background concentration of [OVOCs].

[0017] Optionally, the calculating of the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model comprises: The actual [OVOCs] reaction rate constant, [OVOCs] observed concentration, tracer observed concentration, [OH] exposure and isoprene initial emission concentration are nonlinearly fitted using the photochemical age parameterized model, and the emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursors relative to the tracer, the reaction rate constant of [OVOCs] precursors and OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene and the [OVOCs] background concentration are calculated according to a preset least squares method; The emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursor to the tracer, the reaction rate constant of [OVOCs] precursor with OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene, and the background concentration of [OVOCs] are substituted into the photochemical age parameterized model for calculation to obtain the theoretical concentration of [OVOCs] primary anthropogenic emissions.

[0018] Optionally, the photochemical loss correction and atmospheric diffusion correction are performed on the full-spectrum [VOCs] primary anthropogenic emission observation concentration to obtain the [VOCs] primary anthropogenic emission initial concentration, including: Based on the full spectrum [VOCs] primary anthropogenic emission observation concentration, the following formula is used to calculate [ ]Reaction rate fitting value:

[0019] Among them, [VOCs i ] t represents the observed concentration of anthropogenic emissions of the i-th [VOCs] at time t, ER i represents the emission ratio of the i-th [VOCs] relative to the preset full-spectrum tracer, k model,i represents the reaction rate fitting value of the i-th [VOCs]; Get [VOCs i ] Theoretical reaction rate constant, according to the [VOCs i ]Reaction rate theoretical constant, [VOCs i The reaction rate fitting value is calculated using the following formula to calculate the empirical coefficient:

[0020] in, Represents the empirical coefficient, K theoretical,i represents the theoretical constant of the reaction rate of the i-th [VOCs]; According to the empirical coefficient, the initial concentration of [VOCs] primary anthropogenic emissions is calculated using the following formula:

[0021] Among them, [VOCs i ]0 represents the initial concentration of the first anthropogenic emission of the i-th [VOCs].

[0022] Optionally, after calculating the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model, the method further comprises: Extract atmospheric volatile organic compounds in a pre-constructed atmospheric volatile organic compound set in turn, and obtain multiple groups of tracer-[OVOCs] training sets and multiple groups of tracer-[OVOCs] verification sets of the atmospheric volatile organic compounds, wherein the tracer training concentration in the tracer-[OVOCs] training set is an independent variable, and the [OVOCs] observed training concentration is a dependent variable, each group of tracer-[OVOCs] training sets and each group of tracer-[OVOCs] verification sets correspond to a tracer category, and the multiple groups of tracer-[OVOCs] verification sets are calculated and obtained according to the photochemical age parameterization model, and the atmospheric volatile organic compounds belong to non-methane hydrocarbons or oxygen-containing volatile organic compounds; Extracting initial machine learning algorithms in sequence from a pre-built machine learning algorithm set, wherein the initial machine learning algorithms include: decision tree algorithm, random forest algorithm, extreme random tree algorithm, gradient boosting decision tree algorithm, Adaboost algorithm, Catboost algorithm, Xgboost algorithm, and LightGBM algorithm; Extracting tracer-[OVOCs] training sets in sequence from the plurality of tracer-[OVOCs] training sets; Training the initial machine learning algorithm using the tracer-[OVOCs] training set to obtain a target machine learning algorithm; Using the tracer-[OVOCs] validation set to validate the target machine learning algorithm, obtain a prediction accuracy set, and calculate the prediction accuracy mean of the prediction accuracy set, wherein the prediction accuracy in the prediction accuracy set represents the voting score value of the target machine learning algorithm for the tracer category; Summarize the prediction accuracy means corresponding to each group of tracer-[OVOCs] training sets to obtain the prediction accuracy mean set; Summarize the prediction accuracy mean sets corresponding to each target machine learning algorithm to obtain multiple groups of prediction accuracy mean sets; The following formula is used to calculate the tracer score value according to the plurality of prediction accuracy mean sets to obtain the tracer score value set:

[0023] in, represents the tracer score value of the p-th tracer, Q represents the number of target machine learning algorithms, represents the mean prediction accuracy of the qth target machine learning algorithm for the pth tracer; Extracting a maximum tracer score from the tracer score value set, and identifying a tracer category corresponding to the maximum tracer score; The tracer category corresponding to the maximum tracer score is used as the target tracer of the atmospheric volatile organic compounds. To achieve the above purpose, the present invention also provides a source analysis system for primary anthropogenic emission of VOCs based on a receptor model, comprising: A photochemical age parameterization model construction module is used to obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the tracer reaction rate constant, calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; and construct a photochemical age parameterization model using the [OH] exposure, the initial emission concentration of isoprene, the observed concentration of [OVOCs], the observed concentration of the tracer, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; [OVOCs] primary anthropogenic emission observed concentration calculation module, used to calculate the [OVOCs] primary anthropogenic emission theoretical concentration according to the photochemical age parameterized model, and observe the [NMHCs] anthropogenic emission observed concentration using a pre-built observation instrument to obtain the full-spectrum [VOCs] primary anthropogenic emission observed concentration, wherein the full-spectrum [VOCs] primary anthropogenic emission observed concentration refers to: [OVOCs] primary anthropogenic emission theoretical concentration and [NMHCs] anthropogenic emission observed concentration; [VOCs] initial concentration calculation module for anthropogenic emission, used to perform photochemical loss correction and atmospheric diffusion correction on the full-spectrum [VOCs] initial concentration of anthropogenic emission to obtain the initial concentration of anthropogenic emission of [VOCs]; The [VOCs] source analysis module is used to perform VOCs source analysis based on the pre-built PMF receptor model using the [VOCs] initial concentration of anthropogenic emissions.

[0024] In order to solve the above problem, the present invention further provides an electronic device, the electronic device comprising: a memory storing at least one instruction; and The processor executes the instructions stored in the memory to implement the above-mentioned source analysis method of primary anthropogenic VOCs emission based on the receptor model.

[0025] In order to solve the above problems, the present invention also provides a computer-readable storage medium, in which at least one instruction is stored. The at least one instruction is executed by a processor in an electronic device to implement the above-mentioned source analysis method of primary anthropogenic VOC emissions based on a receptor model.

[0026] In order to solve the problem described in the background technology, the present invention needs to first construct a photochemical age parameterization model, and then calculate the theoretical concentration of [OVOCs] anthropogenic emissions according to the photochemical age parameterization model. Since the photochemical age parameterization model includes known parameters, multiple unknown variables, multiple groups of measurable parameters and multiple groups of calculable parameters, among which the known parameters are: tracer reaction rate constant; the multiple groups of unknown variables include: [OVOCs] relative to tracer emission ratio variable, [OVOCs] precursor relative to tracer emission ratio variable, [OV The reaction rate variables of [OCs] ​​precursors and OH free radicals, the emission ratio variables of [OVOCs] relative to the initial emission concentration of isoprene, and the background concentration of [OVOCs]; multiple groups of measurable parameters include: [OVOCs] observed concentration, tracer observed concentration; multiple groups of calculable parameters include: [OH] exposure, reference [OVOCs] reaction rate constant, and initial emission concentration of isoprene. Therefore, the observed concentration of [OVOCs], the observed concentration of tracers, and the tracer reaction rate constant can be obtained first, and then the [OH] exposure amount, initial emission concentration of isoprene and actual [OVOCs] reaction rate constant. At this time, a photochemical age parameterized model can be constructed according to the parameters obtained and calculated as above and multiple groups of unknown variables. Then, the theoretical concentration of [OVOCs] one-time anthropogenic emission is calculated according to the photochemical age parameterized model, and the observed concentration of [NMHCs] anthropogenic emission is observed using a pre-constructed observation instrument to obtain the observed concentration of [VOCs] one-time anthropogenic emission of the full spectrum, wherein the observed concentration of [VOCs] one-time anthropogenic emission refers to: the theoretical concentration of [OVOCs] one-time anthropogenic emission and the observed concentration of [NMHCs] anthropogenic emission. Due to the existence of atmospheric photochemical loss and atmospheric diffusion, the observed concentration of the full spectrum [VOCs] one-time anthropogenic emission needs to be corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of [VOCs] one-time anthropogenic emission. Finally, according to the pre-constructed PMF receptor model, the initial concentration of [VOCs] one-time anthropogenic emission is used to perform VOCs source analysis, and the source analysis of VOCs emitted by anthropogenic emission based on the receptor model is completed. Therefore, the present invention can improve the analysis accuracy of the PMF receptor model for the sources of VOCs in the actual atmosphere. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] Figure 1 A schematic flow chart of a method for analyzing the source of VOCs emitted by anthropogenic sources based on a receptor model provided in one embodiment of the present invention; Figure 2 A functional module diagram of a source analysis system for anthropogenically emitted VOCs based on a receptor model provided in one embodiment of the present invention; Figure 3 A schematic diagram of the structure of an electronic device for implementing the method for analyzing the source of primary anthropogenically emitted VOCs based on a receptor model provided in one embodiment of the present invention.

[0028] The realization of the purpose, functional features and advantages of the present invention will be further explained in conjunction with embodiments and with reference to the accompanying drawings. DETAILED DESCRIPTION

[0029] It should be understood that the specific embodiments described herein are only used to explain the present invention, and are not used to limit the present invention.

[0030] The embodiment of the present application provides a method for analyzing the sources of VOCs emitted by anthropogenic primary emission based on a receptor model. The execution subject of the method for analyzing the sources of VOCs emitted by anthropogenic primary emission based on a receptor model includes but is not limited to at least one of the electronic devices such as a server and a terminal that can be configured to execute the method provided by the embodiment of the present application. In other words, the method for analyzing the sources of VOCs emitted by anthropogenic primary emission based on a receptor model can be executed by software or hardware installed on a terminal device or a server device, and the software can be a blockchain platform. The server includes but is not limited to: a single server, a server cluster, a cloud server or a cloud server cluster, etc.

[0031] Reference Figure 1 FIG. 1 is a flow chart of a method for analyzing the source of VOCs emitted by anthropogenic sources based on a receptor model according to an embodiment of the present invention. In this embodiment, the method for analyzing the source of VOCs emitted by anthropogenic sources based on a receptor model includes: S1. Obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the tracer reaction rate constant, and calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant.

[0032] It is understandable that the observed concentration of [OVOCs] refers to the observed concentration of a certain oxygenated volatile organic compound (OVOCs) at a monitoring site. For example, the observed concentration of acetaldehyde at a certain site is The tracer observed concentration refers to the observed concentration of compounds at the monitoring site that can specifically indicate or track the source and distribution of the oxygen-containing volatile organic compounds, generally acetylene, , benzene or toluene. It should be noted that in the embodiment of the present invention, volatile organic compounds refer to non-methane hydrocarbons (NMHCs) and oxygen-containing volatile organic compounds.

[0033] Furthermore, the tracer reaction rate constant refers to the reaction rate constant of the tracer and the OH radical, and the [OH] exposure can be expressed as ,in represents the air mass reaction time. The initial emission concentration of isoprene refers to the concentration of isoprene emitted from the emission source. The actual [OVOCs] reaction rate constant refers to the actual reaction rate constant of oxygen-containing volatile organic compounds with OH free radicals.

[0034] In the embodiment of the present invention, after calculating the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model, the method further comprises: Extract atmospheric volatile organic compounds in a pre-constructed atmospheric volatile organic compound set in turn, and obtain multiple groups of tracer-[OVOCs] training sets and multiple groups of tracer-[OVOCs] verification sets of the atmospheric volatile organic compounds, wherein the tracer training concentration in the tracer-[OVOCs] training set is an independent variable, and the [OVOCs] observed training concentration is a dependent variable, each group of tracer-[OVOCs] training sets and each group of tracer-[OVOCs] verification sets correspond to a tracer category, and the multiple groups of tracer-[OVOCs] verification sets are calculated and obtained according to the photochemical age parameterization model, and the atmospheric volatile organic compounds belong to non-methane hydrocarbons or oxygen-containing volatile organic compounds; Extracting initial machine learning algorithms in sequence from a pre-built machine learning algorithm set, wherein the initial machine learning algorithms include: decision tree algorithm, random forest algorithm, extreme random tree algorithm, gradient boosting decision tree algorithm, Adaboost algorithm, Catboost algorithm, Xgboost algorithm, and LightGBM algorithm; Extracting tracer-[OVOCs] training sets in sequence from the plurality of tracer-[OVOCs] training sets; Training the initial machine learning algorithm using the tracer-[OVOCs] training set to obtain a target machine learning algorithm; Using the tracer-[OVOCs] validation set to validate the target machine learning algorithm, obtain a prediction accuracy set, and calculate the prediction accuracy mean of the prediction accuracy set, wherein the prediction accuracy in the prediction accuracy set represents the voting score value of the target machine learning algorithm for the tracer category; Summarize the prediction accuracy means corresponding to each group of tracer-[OVOCs] training sets to obtain the prediction accuracy mean set; Summarize the prediction accuracy mean sets corresponding to each target machine learning algorithm to obtain multiple groups of prediction accuracy mean sets; The following formula is used to calculate the tracer score value according to the plurality of prediction accuracy mean sets to obtain the tracer score value set:

[0035] in, represents the tracer score value of the p-th tracer, Q represents the number of target machine learning algorithms, represents the mean prediction accuracy of the qth target machine learning algorithm for the pth tracer; Extracting a maximum tracer score from the tracer score value set, and identifying a tracer category corresponding to the maximum tracer score; The tracer category corresponding to the maximum tracer score is used as the target tracer of the atmospheric volatile organic compounds.

[0036] In detail, the atmospheric volatile organic compound set may be: acetaldehyde, acetone, methanol, xylene, styrene and toluene, etc. The multiple groups of tracer-[OVOCs] training sets refer to the sets of multiple groups of tracer training concentrations and [OVOCs] observed training concentrations, the tracer training concentrations refer to the tracer concentrations input into the initial machine learning algorithm for training the initial machine learning algorithm, the [OVOCs] observed training concentrations refer to the [OVOCs] observed concentrations observed at the monitoring site under the tracer training concentrations, the tracer-[OVOCs] verification set refers to the sets of multiple groups of tracer verification concentrations and [OVOCs] theoretical verification concentrations, wherein the [OVOCs] theoretical verification concentration refers to the theoretical concentration of [OVOCs] primary anthropogenic emissions calculated according to the photochemical age parameterization model under the tracer verification concentrations.

[0037] Furthermore, the prediction accuracy set refers to the prediction error of each group of tracer verification concentrations and [OVOCs] theoretical verification concentrations in the tracer-[OVOCs] verification set under the initial machine learning algorithm. When performing verification, the tracer verification concentrations can be first input into the target machine learning algorithm, so as to obtain the [OVOCs] theoretical prediction concentration according to the target machine learning algorithm, and then the prediction accuracy is calculated according to the [OVOCs] theoretical verification concentration and the [OVOCs] theoretical prediction concentration. The prediction accuracy can be the reciprocal of the difference between the [OVOCs] theoretical verification concentration and the [OVOCs] theoretical prediction concentration. The prediction accuracy mean refers to the average prediction accuracy of the prediction accuracy set.

[0038] It can be understood that each group of tracer-[OVOCs] training sets corresponds to a tracer category, and therefore, multiple tracer categories correspond to multiple prediction accuracy means, namely the prediction accuracy mean sets. Each prediction accuracy mean set in the multiple prediction accuracy mean sets corresponds to a target machine learning algorithm.

[0039] For example, when the prediction accuracy mean set corresponding to the decision tree algorithm is 2, 4, 7, and 9 (the corresponding tracers are , benzene, toluene and acetylene), the prediction accuracy mean set corresponding to the random forest algorithm is 1, 2, 4, 6 (the corresponding tracers are acetylene, , benzene or toluene), the prediction accuracy mean set corresponding to the extreme random tree algorithm is 3, 4, 6, and 7 (the corresponding tracers are acetylene, , benzene or toluene), the first tracer has a tracer score of 2+1+3, the second tracer has a tracer score of 4+2+4, the third tracer has a tracer score of 7+4+6, and the fourth tracer has a tracer score of 9+6+7. At this time, the fourth tracer, acetylene, is the target tracer of the atmospheric volatile organic compounds. The target tracer refers to the most suitable tracer for the atmospheric volatile organic compounds.

[0040] In the embodiment of the present invention, the calculation of [OH] exposure, initial isoprene emission concentration and actual [OVOCs] reaction rate constant includes: Obtaining an emission concentration ratio and an observed concentration ratio of a benzene series species pair, wherein two benzene series species in the benzene series species pair have the same source and different chemical activities; Calculating [OH] exposure using the emission concentration ratio and the observed concentration ratio according to a pre-constructed [OH] exposure formula; The observed concentration of isoprene was obtained, and the initial emission concentration of isoprene was calculated using the following formula according to the [OH] exposure and the observed concentration of isoprene:

[0041] in, represents the initial emission concentration of isoprene, express The observed isoprene concentration at time is the isoprene reaction rate constant, represents [OH] exposure; The photolysis rate of [OVOCs] was calculated using a pre-built photolysis rate formula; Based on the [OVOCs] photolysis rate, the actual [OVOCs] reaction rate constant was calculated using a pre-constructed [OVOCs] reaction rate formula.

[0042] It is understandable that the benzene series species pair is composed of two benzene series (benzene, toluene and xylene, etc.), for example: ethylbenzene and m-p-xylene, ethylbenzene and o-xylene, benzene and toluene, etc. The two selected benzene series are generally from the same source and have different atmospheric activities. The emission concentration ratio refers to the initial concentration emission ratio of the benzene series species pair at the start of monitoring, which can be obtained by combining the maximum slope concentration ratio of the benzene series species pair actually observed during the period when no photochemical reaction occurs at night (wherein the maximum slope concentration ratio is the ratio of the maximum concentration slope of the more active benzene series to the less active benzene series in the benzene series species pair), or by referring to relevant documents on emission inventories. The observed concentration ratio refers to the ratio of the observed concentrations of the benzene series species pair over a period of time from the start of monitoring.

[0043] In detail, the isoprene observed concentration refers to the observed concentration of isoprene at a monitoring site at a certain moment. The isoprene initial emission concentration refers to the initial emission concentration of isoprene at the emission source. The isoprene reaction rate constant refers to the reaction rate constant of isoprene and OH radicals, the [OVOCs] photolysis rate refers to the photolysis rate of a certain oxygen-containing volatile organic compound (OVOCs), and the actual [OVOCs] reaction rate constant refers to the corrected reaction rate constant of OVOCs and OH radicals.

[0044] In the embodiment of the present invention, the [OH] exposure amount formula is as follows:

[0045] in, represents the exposure amount of [OH], represents the reaction rate constant of the benzene species with species C and OH radical, represents the reaction rate constant of species B in the benzene series with OH radicals, represents the natural logarithm symbol, Indicates that species C and species B are The emission concentration ratio at the time, Indicates that species C and species B are The observed concentration ratio at time.

[0046] Furthermore, the sampling location of the emission concentration ratio is the emission source, and the observation location of the observed concentration ratio is the monitoring site.

[0047] In the embodiment of the present invention, the photolysis rate formula is as follows:

[0048] in, represents the photolysis rate of [OVOCs], represents the photolysis rate of [OVOCs] under clear sky conditions, express The photolysis rate under clear sky conditions, Indicates the measured Photolysis rate. The photolysis rate formula is based on the assumption that factors such as cloud cover have an impact on OVOCs and The effect of the photolysis rate is constructed under the same premise.

[0049] Understandably, and It can be calculated according to the parameterization method in the chemical reaction mechanism MCN (V3.1.1). The specific formula is as follows:

[0050] in, , , Three parameters representing different photolysis reactions are obtained by N, photochemical flux, absorption cross section and quantum yield under clear sky at an altitude of 0.5 km, obtained by fitting using the least squares method, is the solar zenith angle, express or .

[0051] In the embodiment of the present invention, the [OVOCs] reaction rate formula is as follows:

[0052] in, represents the actual [OVOCs] reaction rate constant, represents the reference [OVOCs] reaction rate constant, represents the photolysis rate of [OVOCs], Represents the OH radical concentration.

[0053] It should be understood that the reference [OVOCs] reaction rate constant refers to the reference reaction rate constant of OVOCs with OH radicals.

[0054] S2. Construct a photochemical age parameterized model using the [OH] exposure, initial isoprene emission concentration, [OVOCs] observed concentration, tracer observed concentration, actual [OVOCs] reaction rate constant, and tracer reaction rate constant.

[0055] It can be understood that the photochemical age parameterization model refers to the method for estimating atmospheric OVOCs determined by de Gouw et al. The photochemical age parameterization model decomposes the atmospheric concentration of OVOCs into four parts: [OVOCs] primary anthropogenic emission theoretical concentration, secondary anthropogenic generation concentration, natural source emission concentration and [OVOCs] background concentration. See the following embodiments for details.

[0056] S3. Calculate the theoretical concentration of [OVOCs] primary anthropogenic emissions according to the photochemical age parameterized model, and use a pre-built observation instrument to observe the observed concentration of [NMHCs] anthropogenic emissions to obtain the full-spectrum [VOCs] primary anthropogenic emission observed concentration, wherein the full-spectrum [VOCs] primary anthropogenic emission observed concentration refers to: the theoretical concentration of [OVOCs] primary anthropogenic emissions and the observed concentration of [NMHCs] anthropogenic emissions.

[0057] It is understandable that the theoretical concentration of anthropogenic emissions of [OVOCs] refers to the theoretical concentration of anthropogenic emissions of OVOCs at the emission source calculated according to the photochemical age parameterization model, and the observed concentration of anthropogenic emissions of [NMHCs] refers to the concentration of NMHCs observed at the monitoring site. The observed concentration of anthropogenic emissions of full-spectrum [VOCs] refers to the observed concentration of anthropogenic non-methane hydrocarbons (NMHCs) and oxygen-containing volatile organic compounds at the monitoring site, where atmospheric volatile organic compounds include non-methane hydrocarbons (NMHCs) and oxygen-containing volatile organic compounds (OVOCs).

[0058] Furthermore, since the NMHCs at the monitoring sites are directly derived from primary emissions, and part of the OVOCs are directly emitted from anthropogenic sources, and part are generated by the conversion of NMHCs through photochemical reactions, the theoretical concentration of [OVOCs] primary anthropogenic emissions does not include the concentration of OVOCs generated from secondary NMHCs.

[0059] In the embodiment of the present invention, the photochemical age parameterization model is as follows:

[0060]

[0061]

[0062] in, represents the observed concentration of [OVOCs], represents the observed concentration of the tracer, represents the emission ratio variable of [OVOCs] relative to the tracer, is the tracer reaction rate constant, represents the observed concentration of the tracer, represents the emission ratio of [OVOCs] precursors to tracers, represents the reaction rate variable of [OVOCs] precursor and OH radical, represents the emission ratio variable of [OVOCs] relative to the initial emission concentration of isoprene, Indicates the background concentration of [OVOCs].

[0063] It can be understood that the emission ratio variable of [OVOCs] relative to the tracer refers to the unknown variable of the ratio of the theoretical concentration of anthropogenic emissions of [OVOCs] to the observed concentration of the tracer. The observed concentration of the tracer refers to the observed concentration of the tracer at the monitoring site. The emission ratio variable of [OVOCs] precursor relative to the tracer refers to the unknown variable of the ratio of the emission concentration of the OVOCs precursor to the observed concentration of the tracer. The emission ratio variable of [OVOCs] relative to the initial emission concentration of isoprene refers to the unknown variable of the ratio of the theoretical concentration of anthropogenic emissions of [OVOCs] to the initial emission concentration of isoprene.

[0064] In the embodiment of the present invention, the calculation of the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model includes: The actual [OVOCs] reaction rate constant, [OVOCs] observed concentration, tracer observed concentration, [OH] exposure and isoprene initial emission concentration are nonlinearly fitted using the photochemical age parameterized model, and the emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursors relative to the tracer, the reaction rate constant of [OVOCs] precursors and OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene and the [OVOCs] background concentration are calculated according to a preset least squares method; The emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursor to the tracer, the reaction rate constant of [OVOCs] precursor with OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene, and the background concentration of [OVOCs] are substituted into the photochemical age parameterized model for calculation to obtain the theoretical concentration of [OVOCs] primary anthropogenic emissions.

[0065] It can be understood that the actual [OVOCs] reaction rate constant, [OVOCs] observed concentration, tracer observed concentration, [OH] exposure and isoprene initial emission concentration are substituted into the photochemical age parameterization model, and the difference between the [OVOCs] observed concentration and the [OVOCs] theoretical concentration is minimized, and the least squares method is used to solve , , as well as , and then , , as well as By substituting the specific values ​​of into the photochemical age parameterization model, the theoretical concentration of primary anthropogenic emissions, secondary anthropogenic generation concentration, natural source emission concentration and [OVOCs] background concentration can be calculated respectively.

[0066] S4. Perform photochemical loss correction and atmospheric diffusion correction on the observed concentration of the full-spectrum [VOCs] primary anthropogenic emission to obtain the initial concentration of the [VOCs] primary anthropogenic emission.

[0067] It can be understood that the photochemical loss correction refers to the correction of atmospheric photochemical loss of the full spectrum [VOCs] observed concentration of anthropogenic emissions, and the atmospheric diffusion correction refers to the correction of atmospheric diffusion of the full spectrum [VOCs] observed concentration of anthropogenic emissions.

[0068] In the embodiment of the present invention, the photochemical loss correction and atmospheric diffusion correction are performed on the full spectrum [VOCs] primary anthropogenic emission observation concentration to obtain the [VOCs] primary anthropogenic emission initial concentration, including: Based on the full spectrum [VOCs] primary anthropogenic emission observation concentration, the following formula is used to calculate [ ]Reaction rate fitting value:

[0069] Among them, [VOCs i ] t represents the observed concentration of anthropogenic emissions of the i-th [VOCs] at time t, ER i represents the emission ratio of the i-th [VOCs] relative to the preset full-spectrum tracer, k model,i represents the reaction rate fitting value of the i-th [VOCs]; Get [VOCs i ] Theoretical reaction rate constant, according to the [VOCs i ]Reaction rate theoretical constant, [VOCs i The reaction rate fitting value is calculated using the following formula to calculate the empirical coefficient:

[0070] in, Represents the empirical coefficient, K theoretical,i represents the theoretical constant of the reaction rate of the i-th [VOCs]; According to the empirical coefficient, the initial concentration of [VOCs] primary anthropogenic emissions is calculated using the following formula:

[0071] Among them, [VOCs i ]0 represents the initial concentration of the first anthropogenic emission of the i-th [VOCs].

[0072] Furthermore, the full spectrum tracer refers to the target tracer of the i-th full spectrum [VOCs] selected according to the target machine learning algorithm. The i-th [VOCs] reaction rate fitting value refers to the fitting value of the reaction rate constant of the i-th [VOCs] with the OH free radical, and the [ ] The theoretical reaction rate constant refers to the reaction rate constant of the i-th [VOCs] and OH free radicals measured in the laboratory. The initial concentration of the i-th [VOCs] anthropogenic emission refers to the initial concentration of the i-th [VOCs] emitted from the emission source.

[0073] S5. According to the pre-constructed PMF receptor model, the initial concentration of the [VOCs] primary anthropogenic emission is used to perform VOCs source analysis, and the source analysis of the primary anthropogenic emission of VOCs based on the receptor model is completed.

[0074] It is understandable that the PMF receptor model refers to the Positive Matrix Factorization (PMF) model. Since it is an existing technology to analyze the source of VOCs based on the initial concentration of anthropogenic emissions of [VOCs] and the PMF receptor model, it will not be repeated here.

[0075] In order to solve the problem described in the background technology, the present invention needs to first construct a photochemical age parameterization model, and then calculate the theoretical concentration of [OVOCs] anthropogenic emissions according to the photochemical age parameterization model. Since the photochemical age parameterization model includes known parameters, multiple unknown variables, multiple groups of measurable parameters and multiple groups of calculable parameters, among which the known parameters are: tracer reaction rate constant; the multiple groups of unknown variables include: [OVOCs] relative to tracer emission ratio variable, [OVOCs] precursor relative to tracer emission ratio variable, [OV The reaction rate variables of [OCs] ​​precursors and OH free radicals, the emission ratio variables of [OVOCs] relative to the initial emission concentration of isoprene, and the background concentration of [OVOCs]; multiple groups of measurable parameters include: [OVOCs] observed concentration, tracer observed concentration; multiple groups of calculable parameters include: [OH] exposure, reference [OVOCs] reaction rate constant, and initial emission concentration of isoprene. Therefore, the observed concentration of [OVOCs], the observed concentration of tracers, and the tracer reaction rate constant can be obtained first, and then the [OH] exposure amount, initial emission concentration of isoprene and actual [OVOCs] reaction rate constant. At this time, a photochemical age parameterized model can be constructed according to the parameters obtained and calculated as above and multiple groups of unknown variables. Then, the theoretical concentration of [OVOCs] one-time anthropogenic emission is calculated according to the photochemical age parameterized model, and the observed concentration of [NMHCs] anthropogenic emission is observed using a pre-constructed observation instrument to obtain the observed concentration of [VOCs] one-time anthropogenic emission of the full spectrum, wherein the observed concentration of [VOCs] one-time anthropogenic emission refers to: the theoretical concentration of [OVOCs] one-time anthropogenic emission and the observed concentration of [NMHCs] anthropogenic emission. Due to the existence of atmospheric photochemical loss and atmospheric diffusion, the observed concentration of the full spectrum [VOCs] one-time anthropogenic emission needs to be corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of [VOCs] one-time anthropogenic emission. Finally, according to the pre-constructed PMF receptor model, the initial concentration of [VOCs] one-time anthropogenic emission is used to perform VOCs source analysis, and the source analysis of VOCs emitted by anthropogenic emission based on the receptor model is completed. Therefore, the present invention can improve the analysis accuracy of the PMF receptor model for the sources of VOCs in the actual atmosphere.

[0076] like Figure 2 , which is a functional module diagram of a source analysis system for primary anthropogenic VOC emissions based on a receptor model provided by an embodiment of the present invention.

[0077] The source analysis system 100 for anthropogenic VOCs emission based on a receptor model of the present invention can be installed in an electronic device. According to the functions implemented, the source analysis system 100 for anthropogenic VOCs emission based on a receptor model can include a photochemical age parameterization model construction module 101, a [OVOCs] anthropogenic emission observation concentration calculation module 102, a [VOCs] anthropogenic emission initial concentration calculation module 103 and a [VOCs] source analysis module 104. The module of the present invention can also be referred to as a unit, which refers to a series of computer program segments that can be executed by an electronic device processor and can complete fixed functions, which are stored in the memory of the electronic device.

[0078] The photochemical age parameterized model construction module 101 is used to obtain the observed concentration of [OVOCs], the observed concentration of the tracer and the tracer reaction rate constant, calculate the [OH] exposure, the initial emission concentration of isoprene and the actual [OVOCs] reaction rate constant; and construct a photochemical age parameterized model using the [OH] exposure, the initial emission concentration of isoprene, the observed concentration of [OVOCs], the observed concentration of the tracer, the actual [OVOCs] reaction rate constant and the tracer reaction rate constant; The [OVOCs] primary anthropogenic emission observed concentration calculation module 102 is used to calculate the [OVOCs] primary anthropogenic emission theoretical concentration according to the photochemical age parameterization model, and use a pre-built observation instrument to observe the [NMHCs] anthropogenic emission observed concentration to obtain the full-spectrum [VOCs] primary anthropogenic emission observed concentration, wherein the full-spectrum [VOCs] primary anthropogenic emission observed concentration refers to: the [OVOCs] primary anthropogenic emission theoretical concentration and the [NMHCs] anthropogenic emission observed concentration; The [VOCs] primary anthropogenic emission initial concentration calculation module 103 is used to perform photochemical loss correction and atmospheric diffusion correction on the full-spectrum [VOCs] primary anthropogenic emission observation concentration to obtain the [VOCs] primary anthropogenic emission initial concentration; The [VOCs] source analysis module 104 is used to perform VOCs source analysis based on the pre-constructed PMF receptor model using the [VOCs] initial concentration of anthropogenic emissions.

[0079] In detail, each module in the source analysis system 100 of primary anthropogenic VOCs emission based on the receptor model in the embodiment of the present invention is used in the same manner as described above. Figure 1 The same technical means as the source analysis method of primary anthropogenic VOC emissions based on the receptor model described in the text can produce the same technical effects, so it will not be repeated here.

[0080] like Figure 3, is a schematic diagram of the structure of an electronic device for implementing a method for analyzing the source of primary anthropogenically emitted VOCs based on a receptor model provided by an embodiment of the present invention.

[0081] The electronic device 1 may include a processor 10, a memory 11 and a bus 12, and may also include a computer program stored in the memory 11 and executable on the processor 10, such as a source analysis method program for primary anthropogenic VOC emissions based on a receptor model.

[0082] Among them, the memory 11 includes at least one type of readable storage medium, and the readable storage medium includes flash memory, mobile hard disk, multimedia card, card-type memory (for example: SD or DX memory, etc.), magnetic memory, disk, optical disk, etc. In some embodiments, the memory 11 can be an internal storage unit of the electronic device 1, such as a mobile hard disk of the electronic device 1. In other embodiments, the memory 11 can also be an external storage device of the electronic device 1, such as a plug-in mobile hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (SecureDigital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device 1. Further, the memory 11 also includes an internal storage unit of the electronic device 1 and an external storage device. The memory 11 can not only be used to store application software and various types of data installed in the electronic device 1, such as the code of the source analysis method program of a primary anthropogenic emission of VOCs based on a receptor model, but can also be used to temporarily store data that has been output or is to be output.

[0083] In some embodiments, the processor 10 may be composed of an integrated circuit, for example, a single packaged integrated circuit, or a plurality of packaged integrated circuits with the same or different functions, including one or more central processing units (CPUs), microprocessors, digital processing chips, graphics processors, and combinations of various control chips. The processor 10 is the control core (Control Unit) of the electronic device, and uses various interfaces and lines to connect various components of the entire electronic device, and executes or executes programs or modules stored in the memory 11 (such as a source analysis method program for a primary anthropogenic emission of VOCs based on a receptor model, etc.), and calls data stored in the memory 11 to execute various functions of the electronic device 1 and process data.

[0084] The bus 12 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The bus 12 may be divided into an address bus, a data bus, a control bus, etc. The bus 12 is configured to realize connection and communication between the memory 11 and at least one processor 10, etc.

[0085] Figure 3 Only an electronic device with components is shown, and those skilled in the art will understand that Figure 3 The structure shown does not constitute a limitation on the electronic device 1, and may include fewer or more components than shown in the figure, or combine certain components, or arrange the components differently.

[0086] For example, although not shown, the electronic device 1 may also include a power source (such as a battery) for supplying power to various components. Preferably, the power source may be logically connected to the at least one processor 10 through a power management device, so that the power management device can realize functions such as charging management, discharging management, and power consumption management. The power source may also include one or more DC or AC power sources, recharging devices, power failure detection circuits, power converters or inverters, power status indicators, and other arbitrary components. The electronic device 1 may also include a variety of sensors, Bluetooth modules, Wi-Fi modules, etc., which will not be repeated here.

[0087] Furthermore, the electronic device 1 may also include a network interface. Optionally, the network interface may include a wired interface and / or a wireless interface (such as a WI-FI interface, a Bluetooth interface, etc.), which is generally used to establish a communication connection between the electronic device 1 and other electronic devices.

[0088] Optionally, the electronic device 1 may further include a user interface, which may be a display, an input unit (such as a keyboard), or a standard wired interface or a wireless interface. Optionally, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, and an OLED (Organic Light-Emitting Diode) touch device. The display may also be appropriately referred to as a display screen or a display unit, which is used to display information processed in the electronic device 1 and to display a visual user interface.

[0089] The source analysis method program of primary anthropogenic VOCs emission based on the receptor model stored in the memory 11 of the electronic device 1 is a combination of multiple instructions. When running in the processor 10, it can achieve: Obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the reaction rate constant of the tracer, and calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; A photochemical age parameterization model is constructed using the [OH] exposure, the initial isoprene emission concentration, the [OVOCs] observed concentration, the tracer observed concentration, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; The theoretical concentration of anthropogenic emissions of [OVOCs] is calculated according to the photochemical age parameterized model, and the observed concentration of anthropogenic emissions of [NMHCs] is observed using a pre-built observation instrument to obtain the observed concentration of anthropogenic emissions of the full spectrum [VOCs], wherein the observed concentration of anthropogenic emissions of the full spectrum [VOCs] refers to: the theoretical concentration of anthropogenic emissions of [OVOCs] and the observed concentration of anthropogenic emissions of [NMHCs]; The observed concentration of the full spectrum [VOCs] primary anthropogenic emission is corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of the [VOCs] primary anthropogenic emission; According to the pre-constructed PMF receptor model, the initial concentration of the [VOCs] primary anthropogenic emission is used to perform VOCs source analysis, completing the source analysis of the primary anthropogenic emission of VOCs based on the receptor model.

[0090] Specifically, the specific implementation method of the processor 10 for the above instructions can refer to Figures 1 to 3 The description of the relevant steps in the corresponding embodiments will not be repeated here.

[0091] Furthermore, if the module / unit integrated in the electronic device 1 is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. The computer-readable storage medium can be volatile or non-volatile. For example, the computer-readable medium can include: any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a mobile hard disk, a magnetic disk, an optical disk, a computer memory, and a read-only memory (ROM).

[0092] The present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor of an electronic device, the computer program can implement: Obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the reaction rate constant of the tracer, and calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; A photochemical age parameterization model is constructed using the [OH] exposure, the initial isoprene emission concentration, the [OVOCs] observed concentration, the tracer observed concentration, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; The theoretical concentration of anthropogenic emissions of [OVOCs] is calculated according to the photochemical age parameterized model, and the observed concentration of anthropogenic emissions of [NMHCs] is observed using a pre-built observation instrument to obtain the observed concentration of anthropogenic emissions of the full spectrum [VOCs], wherein the observed concentration of anthropogenic emissions of the full spectrum [VOCs] refers to: the theoretical concentration of anthropogenic emissions of [OVOCs] and the observed concentration of anthropogenic emissions of [NMHCs]; The observed concentration of the full spectrum [VOCs] primary anthropogenic emission is corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of the [VOCs] primary anthropogenic emission; According to the pre-constructed PMF receptor model, the initial concentration of the [VOCs] primary anthropogenic emission is used to perform VOCs source analysis, completing the source analysis of the primary anthropogenic emission of VOCs based on the receptor model.

[0093] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, systems and methods can be implemented in other ways. For example, the system embodiments described above are only illustrative, and actual implementation may have other division methods.

[0094] The modules described as separate components may or may not be physically separated, and the components shown as modules may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0095] In addition, each functional module in each embodiment of the present invention may be integrated into one processing unit, each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of hardware plus software functional modules.

[0096] It is obvious to those skilled in the art that the present invention is not limited to the details of the above exemplary embodiments, and that the present invention can be implemented in other specific forms without departing from the spirit or essential characteristics of the present invention.

[0097] Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solution of the present invention can be modified or replaced by equivalents without departing from the spirit and scope of the technical solution of the present invention.

Claims

1. A method for analyzing the source of anthropogenic VOCs based on a receptor model, characterized in that: The method comprises: Obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the reaction rate constant of the tracer, and calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; A photochemical age parameterization model is constructed using the [OH] exposure, the initial isoprene emission concentration, the [OVOCs] observed concentration, the tracer observed concentration, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; The theoretical concentration of anthropogenic emissions of [OVOCs] is calculated according to the photochemical age parameterized model, and the observed concentration of anthropogenic emissions of [NMHCs] is observed using a pre-built observation instrument to obtain the observed concentration of anthropogenic emissions of the full spectrum [VOCs], wherein the observed concentration of anthropogenic emissions of the full spectrum [VOCs] refers to: the theoretical concentration of anthropogenic emissions of [OVOCs] and the observed concentration of anthropogenic emissions of [NMHCs]; The observed concentration of the full spectrum [VOCs] primary anthropogenic emission is corrected for photochemical loss and atmospheric diffusion to obtain the initial concentration of the [VOCs] primary anthropogenic emission; According to the pre-constructed PMF receptor model, the initial concentration of the [VOCs] primary anthropogenic emission is used to perform VOCs source analysis, completing the source analysis of the primary anthropogenic emission of VOCs based on the receptor model.

2. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 1, characterized in that: The calculation of [OH] exposure, initial isoprene emission concentration and actual [OVOCs] reaction rate constant includes: Obtaining an emission concentration ratio and an observed concentration ratio of a benzene series species pair, wherein two benzene series species in the benzene series species pair have the same source and different chemical activities; Calculating [OH] exposure using the emission concentration ratio and the observed concentration ratio according to a pre-constructed [OH] exposure formula; The observed concentration of isoprene was obtained, and the initial emission concentration of isoprene was calculated using the following formula according to the [OH] exposure and the observed concentration of isoprene: in, represents the initial emission concentration of isoprene, express The observed isoprene concentration at time is the isoprene reaction rate constant, represents [OH] exposure; The photolysis rate of [OVOCs] was calculated using the pre-built photolysis rate formula; Based on the [OVOCs] photolysis rate, the actual [OVOCs] reaction rate constant was calculated using a pre-constructed [OVOCs] reaction rate formula.

3. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 2, characterized in that: The [OH] exposure formula is as follows: in, represents the exposure amount of [OH], represents the reaction rate constant of the benzene series species C and the OH radical, represents the reaction rate constant of species B in the benzene series with OH radicals, represents the natural logarithm symbol, Indicates that species C and species B are The emission concentration ratio at the time, Indicates that species C and species B are The observed concentration ratio at time.

4. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 3, characterized in that: The photolysis rate formula is as follows: in, represents the photolysis rate of [OVOCs], represents the photolysis rate of [OVOCs] under clear sky conditions, express The photolysis rate under clear sky conditions, Indicates the measured Photolysis rate.

5. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 4, characterized in that: The [OVOCs] reaction rate formula is as follows: in, represents the actual [OVOCs] reaction rate constant, represents the reference [OVOCs] reaction rate constant, represents the photolysis rate of [OVOCs], Represents the OH radical concentration.

6. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 5, characterized in that: The photochemical age parameterization model is as follows: in, represents the observed concentration of [OVOCs], represents the observed concentration of the tracer, represents the emission ratio variable of [OVOCs] relative to the tracer, is the tracer reaction rate constant, represents the observed concentration of the tracer, represents the emission ratio of [OVOCs] precursors to tracers, represents the reaction rate variable of [OVOCs] precursor and OH radical, represents the emission ratio variable of [OVOCs] relative to the initial emission concentration of isoprene, Indicates the background concentration of [OVOCs].

7. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model according to claim 6, characterized in that: The method of calculating the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model comprises: The actual [OVOCs] reaction rate constant, [OVOCs] observed concentration, tracer observed concentration, [OH] exposure and isoprene initial emission concentration are nonlinearly fitted using the photochemical age parameterized model, and the emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursors relative to the tracer, the reaction rate constant of [OVOCs] precursors and OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene and the [OVOCs] background concentration are calculated according to a preset least squares method; The emission ratio of [OVOCs] relative to the tracer, the emission ratio of [OVOCs] precursor to the tracer, the reaction rate constant of [OVOCs] precursor with OH radicals, the emission ratio of [OVOCs] relative to the initial emission concentration of isoprene, and the background concentration of [OVOCs] are substituted into the photochemical age parameterized model for calculation to obtain the theoretical concentration of [OVOCs] primary anthropogenic emissions.

8. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 7, characterized in that: The photochemical loss correction and atmospheric diffusion correction are performed on the observed concentration of the full spectrum [VOCs] primary anthropogenic emission to obtain the initial concentration of the [VOCs] primary anthropogenic emission, including: Based on the full spectrum [VOCs] primary anthropogenic emission observation concentration, [VOCs] was calculated using the following formula i ]Reaction rate fitting value: Among them, [VOCs i ] t represents the observed concentration of anthropogenic emissions of the i-th [VOCs] at time t, ER i represents the emission ratio of the i-th [VOCs] relative to the preset full-spectrum tracer, k model,i represents the reaction rate fitting value of the i-th [VOCs]; Get [VOCs i ] Theoretical reaction rate constant, according to the [VOCs i ]Reaction rate theoretical constant, [VOCs i The reaction rate fitting value is calculated using the following formula to calculate the empirical coefficient: Where D represents the empirical coefficient, K theoretical,i represents the theoretical constant of the reaction rate of the i-th [VOCs]; According to the empirical coefficient, the initial concentration of [VOCs] primary anthropogenic emissions is calculated using the following formula: Among them, [VOCs i ]0 represents the initial concentration of the first anthropogenic emission of the i-th [VOCs].

9. The method for analyzing the source of primary anthropogenic VOCs emission based on a receptor model as claimed in claim 8, characterized in that: After calculating the theoretical concentration of anthropogenic emissions of [OVOCs] according to the photochemical age parameterized model, the method further comprises: Extract atmospheric volatile organic compounds in a pre-constructed atmospheric volatile organic compound set in turn, and obtain multiple groups of tracer-[OVOCs] training sets and multiple groups of tracer-[OVOCs] verification sets of the atmospheric volatile organic compounds, wherein the tracer training concentration in the tracer-[OVOCs] training set is an independent variable, and the [OVOCs] observed training concentration is a dependent variable, each group of tracer-[OVOCs] training sets and each group of tracer-[OVOCs] verification sets correspond to a tracer category, and the multiple groups of tracer-[OVOCs] verification sets are calculated and obtained according to the photochemical age parameterization model, and the atmospheric volatile organic compounds belong to non-methane hydrocarbons or oxygen-containing volatile organic compounds; Extracting initial machine learning algorithms in sequence from a pre-built machine learning algorithm set, wherein the initial machine learning algorithms include: decision tree algorithm, random forest algorithm, extreme random tree algorithm, gradient boosting decision tree algorithm, Adaboost algorithm, Catboost algorithm, Xgboost algorithm, and LightGBM algorithm; Extracting tracer-[OVOCs] training sets in sequence from the plurality of tracer-[OVOCs] training sets; Training the initial machine learning algorithm using the tracer-[OVOCs] training set to obtain a target machine learning algorithm; Using the tracer-[OVOCs] validation set to validate the target machine learning algorithm, obtain a prediction accuracy set, and calculate the prediction accuracy mean of the prediction accuracy set, wherein the prediction accuracy in the prediction accuracy set represents the voting score value of the target machine learning algorithm for the tracer category; Summarize the prediction accuracy means corresponding to each group of tracer-[OVOCs] training sets to obtain the prediction accuracy mean set; Summarize the prediction accuracy mean sets corresponding to each target machine learning algorithm to obtain multiple groups of prediction accuracy mean sets; The following formula is used to calculate the tracer score value according to the plurality of prediction accuracy mean sets to obtain the tracer score value set: in, represents the tracer score value of the p-th tracer, Q represents the number of target machine learning algorithms, represents the mean prediction accuracy of the qth target machine learning algorithm for the pth tracer; Extracting a maximum tracer score from the tracer score value set, and identifying a tracer category corresponding to the maximum tracer score; The tracer category corresponding to the maximum tracer score is used as the target tracer of the atmospheric volatile organic compounds.

10. A source analysis system for primary anthropogenic VOCs emissions based on a receptor model, characterized in that: The system comprises: A photochemical age parameterization model construction module is used to obtain the observed concentration of [OVOCs], the observed concentration of the tracer, and the tracer reaction rate constant, calculate the [OH] exposure, the initial emission concentration of isoprene, and the actual [OVOCs] reaction rate constant; and construct a photochemical age parameterization model using the [OH] exposure, the initial emission concentration of isoprene, the observed concentration of [OVOCs], the observed concentration of the tracer, the actual [OVOCs] reaction rate constant, and the tracer reaction rate constant; [OVOCs] primary anthropogenic emission observed concentration calculation module, used to calculate the [OVOCs] primary anthropogenic emission theoretical concentration according to the photochemical age parameterized model, and observe the [NMHCs] anthropogenic emission observed concentration using a pre-built observation instrument to obtain the full-spectrum [VOCs] primary anthropogenic emission observed concentration, wherein the full-spectrum [VOCs] primary anthropogenic emission observed concentration refers to: [OVOCs] primary anthropogenic emission theoretical concentration and [NMHCs] anthropogenic emission observed concentration; [VOCs] initial concentration calculation module for anthropogenic emission, used to perform photochemical loss correction and atmospheric diffusion correction on the full-spectrum [VOCs] initial concentration of anthropogenic emission to obtain the initial concentration of anthropogenic emission of [VOCs]; The [VOCs] source analysis module is used to perform VOCs source analysis based on the pre-built PMF receptor model using the [VOCs] initial concentration of anthropogenic emissions.

Images