Risk-oriented industrial park VOCs comprehensive tracing method and system

Abstract

The invention relates to the technical field of air pollution control, in particular to a risk-oriented industrial park VOCs comprehensive tracing method and system, and the method comprises the steps: setting monitoring distribution points and monitoring elements, and carrying out the monitoring of the monitoring elements through a preset monitoring method; setting risk factors, and evaluating the risk factors according to the monitored monitoring factors; performing risk source quantitative traceability to generate a risk source quantitative traceability result; wherein the risk source quantitative traceability comprises PMF source analysis, receptor and pollution source component spectrum matching, risk source comprehensive contribution evaluation and risk source space matching; performing risk assessment uncertainty and sensitivity analysis on a risk source quantitative traceability result; according to a risk assessment uncertainty and sensitivity analysis result and a risk source quantitative tracing result, adopting IMTs to formulate a comprehensive emission reduction method for VOCs in the industrial park; according to the scheme, VOCs in the industrial park can be precisely traced, the uncertainty is reduced, and the fine control capability of pollutants in the industrial park is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of air pollution control, and in particular to a risk-oriented VOCs comprehensive source tracing method and system for industrial parks. Background Art

[0002] Volatile organic compounds (VOCs) can pose a serious threat to ambient air quality, human health, climate change and ecological environment through direct or indirect effects. In addition, some VOCs substances emitted by industries (such as benzene, organic sulfur, aldehydes, ketones and esters) also produce foul odors, which have become a hot issue for public complaints in recent years. With the development of urban industrialization, industrial production activities have been recognized as the main anthropogenic source of VOCs emissions, accounting for more than 50% of all anthropogenic VOCs emissions. As the main carrier of modern industrial development, industrial parks are characterized by huge and concentrated pollutant emissions and are a very important VOCs emission entity.

[0003] There are many enterprises involved in VOCs emissions in industrial parks, and the sources of VOCs emissions are complicated. The formulation of VOCs control policies can usually only be based on subjective judgments based on information such as corporate environmental protection data, pollution collection and treatment conditions, and the operation of treatment facilities, which is subject to great uncertainty. Obviously, this method can no longer meet the needs of refined VOCs control in industrial parks at this stage. However, there is currently no mature VOCs source tracing technology method that can be applied to the industrial park scale, and the currently commonly used ambient air VOCs source tracing technology can only obtain a rough and large-scale source analysis result, and there is still a certain degree of uncertainty, which makes it difficult for the source analysis results to effectively support the formulation of VOCs control policies in industrial parks.

[0004] Therefore, there is an urgent need for a risk-oriented comprehensive traceability method for VOCs in industrial parks, which can accurately trace the source of VOCs in industrial parks and reduce uncertainty, so as to improve the refined control capabilities of pollutants in industrial parks. Summary of the invention

[0005] One of the purposes of the present invention is to provide a risk-oriented comprehensive traceability method for VOCs in industrial parks, which can accurately trace the source of VOCs in industrial parks, reduce uncertainty, and improve the ability to fine-tune pollutant management and control in industrial parks.

[0006] The basic solution 1 provided by the present invention is a risk-oriented comprehensive VOCs source tracing method for industrial parks, including the following contents:

[0007] S1. Set up monitoring points in the industrial park to be monitored;

[0008] S2. Set monitoring elements of monitoring points, and use preset monitoring methods to monitor the monitoring points and monitor the monitoring elements;

[0009] S3. Set risk factors, and conduct risk factor assessment based on the monitored factors to generate risk factor assessment results;

[0010] S4. Conduct quantitative risk source tracing and generate risk source quantitative tracing results; quantitative risk source tracing includes: PMF source analysis, receptor and pollution source component spectrum matching, risk source comprehensive contribution assessment and risk source spatial matching;

[0011] S5. Conduct risk assessment uncertainty and sensitivity analysis on the risk factor assessment results and risk source quantitative traceability results, and generate risk assessment uncertainty and sensitivity analysis results;

[0012] S6. Based on the risk factor assessment results, the quantitative traceability results of risk sources, and the uncertainty and sensitivity analysis results of risk assessment, IMTs are used to develop a comprehensive VOCs emission reduction method for industrial parks;

[0013] Wherein, the S4 includes: S401, according to the risk factor assessment results, using the PMF model, performing PMF analysis, generating a VOCs source component spectrum, and determining the primary pollution source in the pollution source; the pollution source includes the primary pollution source and the secondary pollution source, the primary pollution source is the industry type; the secondary pollution source is each emission link included in the primary pollution source;

[0014] S403. Use a multivariate linear regression model to perform secondary matching on the VOCs source component spectrum obtained by PMF analysis and the VOCs source component spectrum obtained by measuring the secondary pollution sources included in each primary pollution source, obtain the relationship between the VOCs source component spectrum of each secondary pollution source and the VOCs source component spectrum analyzed by PMF, and obtain the secondary source risk allocation result.

[0015] Beneficial effects of this program: This program is based on the full identification of pollution sources and environmental receptors, and is guided by the mitigation of health risks, environmental risks and odor nuisance risks caused by VOCs emissions in industrial parks. It couples multiple VOCs traceability methods to quantify the source analysis results to the production workshop / production line and other links of the enterprise, and calculates the maximum concentration of VOCs allowed to stay in a specific location (without generating the above three risks) to formulate the emission reduction rate of VOCs species for enterprises in various industries; it provides data support for the precise traceability of VOCs in industrial parks, and at the same time complies with the principles of environmental priority and people-oriented, and promotes the precise control and refined emission reduction of VOCs in industrial parks;

[0016] Specifically, this plan sets up detailed monitoring points and monitoring elements, and sets up corresponding monitoring methods for monitoring. Based on the monitoring results, various subsequent analyses are carried out, including risk factor assessment, quantitative tracing of risk sources, uncertainty and sensitivity analysis, and IMTs are used to formulate a comprehensive VOCs emission reduction method for industrial parks based on the risk factor assessment results, the risk source quantitative tracing results, and the risk assessment uncertainty and sensitivity analysis results; the risk source quantitative tracing also sets up a two-level pollution source analysis. Compared with the traditional analysis that can only determine the industry type (primary pollution source), this plan also uses a multivariate linear regression model to perform a secondary match between the VOCs source component spectrum obtained by PMF analysis and the VOCs source component spectrum obtained by the actual measurement of the secondary pollution sources included in each primary pollution source, to obtain the relationship between the VOCs source component spectrum of each secondary pollution source and the VOCs source component spectrum analyzed by PMF, to obtain the secondary source risk allocation results, and to improve the accuracy and pertinence of risk source analysis in industrial parks;

[0017] In summary, this solution can accurately trace the source of VOCs in industrial parks, reduce uncertainty, and improve the ability to fine-tune pollutant control in industrial parks.

[0018] The second purpose of the present invention is to provide a risk-oriented comprehensive VOCs traceability system for industrial parks, which can accurately trace the VOCs in industrial parks, reduce uncertainty, and improve the industrial park's ability to fine-tune pollutant control.

[0019] The present invention provides basic solution 2: a risk-oriented comprehensive VOCs traceability system for industrial parks, which adopts the above-mentioned risk-oriented comprehensive VOCs traceability method for industrial parks. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 It is a flow chart of an embodiment of a risk-oriented comprehensive traceability method for VOCs in industrial parks of the present invention;

[0021] Figure 2 A schematic diagram of sampling of VOCs component spectrum at the emission outlet of a pollution source in an embodiment of the risk-oriented VOCs comprehensive source tracing method for an industrial park of the present invention;

[0022] Figure 3 It is a flow chart of a multivariate linear regression model in an embodiment of a risk-oriented comprehensive traceability method for VOCs in an industrial park of the present invention;

[0023] Figure 4 This is an example diagram of the CPF spatial distribution map in the embodiment of the risk-oriented comprehensive traceability method for VOCs in industrial parks of the present invention. DETAILED DESCRIPTION

[0024] The following is further described in detail through specific implementation methods:

[0025] The embodiment is basically as shown in the attached Figure 1 As shown: A risk-oriented comprehensive traceability method for VOCs in industrial parks, including the following:

[0026] S1. Monitoring points are set up in the industrial park to be monitored. The monitoring points include: monitoring points for environmental receptors and monitoring points for pollution sources. The monitoring points for environmental receptors also serve as monitoring points for meteorological data to monitor meteorological data.

[0027] Specifically, the environmental receptor VOCs monitoring in industrial parks adopts high-resolution online monitoring, and the monitoring points are set up in a grid-based manner according to the relevant requirements of the "Technical Specifications for the Layout of Ambient Air Quality Monitoring Points (Trial)" (HJ-664);

[0028] Taking into account the diffusion and transmission laws of atmospheric pollutants and their impact on ambient air quality, the range of monitoring points set up in pollution concentration areas is generally a radius of 100 to 500 meters; when considering the impact of point sources with higher emission altitudes on point concentrations, it can also be expanded to a radius of 500 meters to 4 kilometers.

[0029] The grid monitoring points of environmental receptor VOCs in industrial parks are arranged according to specific needs, mainly considering the following points:

[0030] (1) In areas where pollution is concentrated in industrial parks, pollution monitoring points should be arranged with full consideration of high-concentration areas that are likely to cause harm to human health and areas that may have a significant impact on ambient air quality;

[0031] (2) For grid points with medium and high emission intensity and major pollution projects, they should be set up in the maximum ground concentration area downwind of the dominant wind direction and the second wind direction of the pollution source (generally the dominant risk in the heaviest pollution season);

[0032] (3) Based on the dominant wind direction of the industrial park, background monitoring points are arranged at upwind locations around the industrial park.

[0033] The monitoring of pollution sources selects all the types of industries covered by the industrial park, and selects representative enterprises from each type of industry for monitoring; the sampling methods are divided into organized sampling and unorganized sampling; organized sampling covers all air-related emission outlets of polluting enterprises, and the determination of sampling points and sampling locations refers to the relevant provisions of GB / T 16157 "Methods for Determination of Particle Matter and Sampling of Gaseous Pollutants in Exhaust Gas from Fixed Pollution Sources" and HJ / T 397 "Technical Specifications for Monitoring Exhaust Gas from Fixed Sources"; unorganized emission sampling covers important VOCs unorganized emission points of enterprises, such as workshops, production lines, storage tanks, sewage treatment pools, etc. It is different from the sampling points specified in relevant standards; the purpose of unorganized sampling is to identify the VOCs component spectrum characteristics of important unorganized emission links of enterprises, so unorganized sampling points are required to be set at the direct occurrence of VOCs unorganized emissions, such as workshop ventilation outlets, main unorganized emission points of production lines, storage tank leaks, and the liquid surface of the open sewage treatment pool. In addition, atmospheric VOCs samples can also be collected at the center of the factory boundary to characterize the atmospheric VOCs pollution characteristics of the industry area.

[0034] S2. Set monitoring elements of monitoring points, and use preset monitoring methods to monitor the monitoring points and monitor the monitoring elements;

[0035] Among them, the monitoring elements are based on the environmental risk and health risk characteristics of VOCs components, and refer to the VOCs components recommended in the "2018 Monitoring Plan for Volatile Organic Compounds in Ambient Air in Key Areas" issued by the Ministry of Ecology and Environment for monitoring, including alkanes, alkenes, alkynes, aromatic hydrocarbons, oxygenated volatile organic compounds (OVOCs), halogenated hydrocarbons and other components. Specifically, the detailed list of monitoring elements is as follows Table 1;

[0036] Pre-set monitoring methods, including: environmental receptor monitoring method, pollution source monitoring method;

[0037] Environmental receptor monitoring methods are as follows:

[0038] The VOCs monitoring of environmental receptors in industrial parks uses gas chromatography-mass spectrometry (GCMS-FID) or equivalent equipment to carry out high-resolution continuous and automatic sample collection and analysis in a preset time period; the quality assurance and quality control procedures adopted during the online monitoring process refer to the "Technical Requirements and Test Methods for Gas Chromatography Continuous Monitoring System of Volatile Organic Compounds in Ambient Air" (HJ1010-2018).

[0039] Pollution source monitoring methods include: sample collection methods and sample analysis methods, as follows:

[0040] Sample collection method, specifically:

[0041] Pollution source monitoring uses Summa tanks and 2,4-dinitrophenylhydrazine derivative compound (DNPH) absorption tubes / absorption liquids to collect pollution source discharge outlets and unorganized emission samples;

[0042] When the polluting enterprise is operating under normal conditions, collect emission samples from the sampling port at the rear end of the exhaust gas treatment facility;

[0043] Before sample collection, a portable volatile organic compound detector (PID) was used to detect the initial concentration of TVOCs in the exhaust gas;

[0044] When the TVOCs concentration is higher than the preset concentration index, the sample in the Suma tank needs to be diluted a certain multiple before component analysis can be performed on the GC-MS / FID to avoid excessive TVOCs concentration polluting the analysis system and excessive dilution of VOCs concentration below the method detection limit; the preset concentration index is set according to demand, and the preset concentration index in this embodiment is 10ppmv;

[0045] like Figure 2 As shown, the organized exhaust gas is collected through a sampling gun, which should have the function of removing particles and heating. After the stable organic exhaust gas is filtered and heated, the flow / pressure controller of the Suma tank realizes the timed and quantitative collection of VOCs, and the flow controller of the DNPH absorption liquid realizes the timed and quantitative collection of VOCs;

[0046] For the samples of the organized emission, a preset number range of samples are collected at equal time intervals within a preset time period using a Suma tank and a DNPH absorption liquid; wherein the preset time period and the preset number range are set according to the demand. In this embodiment, the preset time period is 1 hour, and the preset number range is 3 to 4;

[0047] If the organized emission is intermittent, a preset number of samples shall be collected at equal time intervals within the preset time period during the high-intensity emission period; high-intensity emission is considered high-intensity emission when the production load of the enterprise is higher than 80% and the pollution-producing process continues to operate normally;

[0048] The samples of unorganized exhaust gas and VOCs in the central atmosphere of the plant boundary are collected using Suma tanks and DNPH absorption tubes. Under stable production conditions of the enterprise, the Suma tank is equipped with a flow limiting valve and the DNPH absorption tube is equipped with a flow controller for continuous sampling for a preset time. The Suma tank and DNPH absorption tube at each point collect a preset number of samples respectively; the preset time and the preset number are set according to the needs. In this embodiment, the preset time is 1 hour and the preset number is 3;

[0049] Sample analysis method, specifically:

[0050] The samples with too high VOCs concentration collected by the Suma tank need to be diluted a certain multiple in the laboratory before they can be analyzed on the machine. The concentration after dilution is controlled within the preset concentration range. The analysis method and quality assurance and quality control measures refer to "Determination of Volatile Organic Compounds in Ambient Air by Sampling / Gas Chromatography-Mass Spectrometry" HJ 759-2015; the preset concentration range is set according to the needs, which is 0.1-10ppmv in this embodiment;

[0051] VOCs samples collected by DNPH adsorption liquid / adsorption tube need to be pretreated before analysis on the instrument;

[0052] After pretreatment, a preset volume of sample is transferred to a sample injection bottle, and then analyzed using HPLC; the preset volume is set according to demand, and in this embodiment is 2 mL;

[0053] The sample pretreatment method and the working conditions of the HPLC were carried out in accordance with the national standard methods (HJ683-2014 and HJ 1153-2020);

[0054] The pretreated samples were analyzed according to the US Environmental Protection Agency (EPA) TO-11A method and high performance liquid chromatography (HPLC, Shimadzu).

[0055] S3. Set risk factors, and conduct risk factor assessment based on the monitored factors to generate risk factor assessment results;

[0056] The first-level indicators of risk factors are mainly divided into two categories: health risk and environmental effect risk;

[0057] The health risk in this plan refers to the risk of pathogenicity, carcinogenicity, and odor pollution caused by VOCs emitted from industrial parks into the atmosphere, which poses a health threat to the human body;

[0058] Environmental effect risk refers to the direct or indirect effects of VOCs emitted into the atmosphere from industrial parks, which may lead to risks such as deterioration of air quality, climate change, and depletion of the ozone layer.

[0059] Secondary indicators of health risk include: hazard entropy (HI), lifetime cancer risk (LCR), and odor activity value (OAV);

[0060] Secondary indicators of environmental effect risk include: ozone creation potential (OFP) and secondary organic aerosol creation potential (SOAFP);

[0061] The VOCs species and risks involved in each secondary indicator are different. Industrial parks can select secondary indicators as risk assessment factors based on the purpose of VOCs control and the focus of prevention and control. The risk factors involved in the monitored VOCs species are shown in Table 1;

[0062] Table 1: List of VOCs species monitored by environmental receptors and pollution sources and their risk properties

[0063]

[0064]

[0065]

[0066]

[0067]

[0068]

[0069]

[0070]

[0071]

[0072] Among them, risk factor assessment is carried out according to the monitored elements, including:

[0073] OFP represents the maximum contribution of different VOCs species to ozone generation under optimal reaction conditions. The calculation formula is as follows:

[0074] OFP i =C i ×MIR i ;

[0075] In the formula, OFP i represents the ozone generation potential of VOCs species i, dimensionless; C i represents the concentration of VOCs species i during the monitoring period, μg·m-3; MIR i represents the maximum response increment of VOCs species i, MIR i The value of refers to Table 2;

[0076] Table 2: MIR values ​​and SOA yields of VOCs species

[0077]

[0078]

[0079]

[0080]

[0081]

[0082]

[0083] SOAFP represents the formation capacity of secondary organic aerosol (SOA) of different VOCs species, and the calculation formula is as follows:

[0084] SOAFP i =C i ×Y SOA,i ;

[0085] In the formula, SOAFP i represents the secondary organic aerosol generation potential of VOCs species i, μg·m-3; C i represents the concentration of VOCs species i during the monitoring period, μg·m-3; YSOA,i represents the SOA yield of VOCs species i, Y SOA,i Please refer to Table 2 for the value of .

[0086] Hazard entropy (HQ) of specific VOCs species i ) and the cumulative hazard entropy (HI) of all VOCs species are calculated by the following formula:

[0087]

[0088] In the formula, HQ i is the non-cancer inhalation hazard entropy of a specific VOCs species i, usually considered to be HQ i When it is greater than 1, it means that the substance may cause adverse health effects due to inhalation; REL i is the reference exposure level for species i, REL i Refer to Table 3.

[0089] Table 3: Reference Exposure Levels (RELs) for Specific VOCs Species

[0090]

[0091]

[0092] Carcinogenic risk (LCR) of specific VOCs species i ) and the cumulative carcinogenic risk (LCR) of all VOCs species are calculated by the following formula:

[0093]

[0094]

[0095] Where, LCR iis the carcinogenic risk of a specific VOCs species i; DOSEdaily represents the daily inhalation dose of harmful VOCs (unit: mg·kg-1·d-1), which is calculated from the concentration of species i (C i , μg·m-3), normalized respiration rate (BR / BW, L·kg-1·d-1), inhalation absorption factor (A, dimensionless) and exposure frequency (EF, %); CSF i is the inhalation cancer slope factor for species i, dimensionless; ASFa is the sensitivity coefficient, dimensionless; EDa, AT and FAHa are exposure duration (years), lifetime exposure time (years) and time at home (dimensionless), respectively; the parameters involved in the formula refer to Appendix 4 and Appendix 5 for the exposure parameters of each age group a;

[0096] Table 4: Inhalation cancer slope factors for 21 carcinogenic VOCs species

[0097]

[0098] Table 5: Exposure parameters for inhalation carcinogenic risk in residents of different age groups

[0099]

[0100]

[0101] The odor intensity of individual odorants is described by the Odor Activity Value (OAV), which is defined as the ratio between the concentration of a compound and its odor threshold:

[0102]

[0103] In the formula, C i is the concentration of VOCs odor substance i, 1×10–6; OT i is the threshold value of odor substance i, 1×10–6; OAV i OAV is the odor activity value of odor substance i, dimensionless; sum It is the sum of the odor activity values ​​of all odor substances evaluated, dimensionless; the odor threshold of VOCs refers to Table 6:

[0104] Table 6: 54 odorous VOCs and their odor thresholds

[0105]

[0106]

[0107] S4. Conduct quantitative risk source tracing and generate risk source quantitative tracing results; quantitative risk source tracing includes: PMF source analysis, comprehensive risk source contribution assessment, receptor and pollution source component spectrum matching and risk source spatial matching;

[0108] Specifically: According to the risk factor assessment results, the PMF model is used to perform PMF analysis, generate the VOCs source component spectrum, and determine the primary pollution source among the pollution sources; the pollution sources include primary pollution sources and secondary pollution sources, the primary pollution source is the industry type, such as petrochemical, pharmaceutical manufacturing, automobile manufacturing, etc.; the secondary pollution source is the various emission links included in the primary pollution source, that is, the emission environment under the industry type, see Table 7; the primary pollution source among the pollution sources is determined by obtaining the VOCs source component spectrum of each source. The pollutant content in the VOCs source component spectrum of different industries is different, so as to determine the industry to which the VOCs source component spectrum belongs;

[0109] Based on the VOCs source component spectrum of various pollution sources analyzed by the PMF model, as well as the primary pollution sources, the absolute value of the concentration of species in each primary pollution source is used to quantify the contribution of the primary pollution source to the environmental effect risk and health risk;

[0110] A multivariate linear regression model is used to perform a secondary match between the VOCs source component spectrum obtained by PMF analysis and the VOCs source component spectrum obtained by measuring the secondary pollution sources included in each primary pollution source, and the relationship between the VOCs source component spectrum of each secondary pollution source (each emission link of the polluting industry) and the VOCs source component spectrum analyzed by PMF is obtained. The risk source is further refined to obtain the secondary source risk allocation result; the pollution source with the set risk factors is identified as a risk source;

[0111] The conditional bivariate probability function (CBPF) analysis was used to verify the risk source allocation results and identify the impact of risk sources under different wind directions and wind speeds, as well as the spatial distribution of the main sources of different risk factors;

[0112] The risk source allocation results include: the risk allocation results of the primary pollution sources and secondary pollution sources among the determined pollution sources;

[0113] The risk factors are VOCs species with risk factors involved; the risk sources determined by matching receptors with pollution sources are identified based on the coupling relationship between data; CBPF determines the risk sources based on meteorology, the actual spatial distribution of pollution sources, and the spatial distribution of risk factors. The two methods match and corroborate each other to improve the accuracy of risk source determination;

[0114] The specific contents are as follows:

[0115] PMF source analysis:

[0116] The PMF (Orthogonal Matrix Factor) model is a receptor model for multivariate factor analysis. The model operation can decompose the receptor sample matrix (i×j dimension) into a contribution matrix (i×k) and a factor matrix (k×j):

[0117]

[0118] In the formula, x ij is the contribution of the jth species in the ith sample; g ik is the contribution of the kth source in the ith sample; f kj is the concentration of the jth species in the kth source; e ij is the error score;

[0119] The factor distribution graph and contribution are obtained by minimizing the objective function Q using the least squares method:

[0120]

[0121] Where U ij is the uncertainty for each species in each sample.

[0122] To ensure the accuracy and reliability of the model operation results, the VOCs species concentrations input into the PMF model should include the VOCs species covered by the selected risk factors. The indicator parameters of the operation results should meet the Qtrue / Qrobust value between 1 and 1.5, the correlation coefficient (R2) between the monitored and predicted VOCs species concentrations should be greater than 0.6, and the VOCs species residuals should be normally distributed between -3 and 3.

[0123] Receptor and pollution source spectrum matching:

[0124] The PMF model can only obtain preliminary and rough VOCs source allocation results; in order to improve the accuracy and pertinence of risk source analysis in industrial parks, this scheme uses a multivariate linear regression model to perform a secondary match between the VOCs source component spectrum obtained by PMF analysis and the VOCs source component spectrum obtained by actual measurement of pollution sources;

[0125] This scheme quantifies the risk contribution of the secondary pollution sources involved in each primary pollution source (secondary sources include: emissions from workshops, production lines and sewage treatment, as well as all links involving VOCs emissions such as organized emissions and unorganized emissions). At the same time, it can also evaluate the accuracy of the PMF model analysis results (i.e., the PMF model operation results) and the representativeness of the VOCs source component spectrum measurement of the primary pollution source. The subdivision of the secondary pollution source is determined based on the actual emission situation. The more comprehensive the secondary pollution sources covered by the primary pollution source, the higher the match with the source component spectrum analyzed by the receptor.

[0126] The role of the multivariate linear regression model is to evaluate the role and contribution of many independent variables x to the dependent variable y. In the present invention, the VOCs source component spectrum analyzed by the PMF model is used as the dependent variable y, and the VOCs source component spectrum of each secondary pollution source in the industrial park industry is used as the dependent variable x, so as to evaluate the contribution of each secondary pollution source (such as workshops, production lines, sewage treatment, boilers, etc., organized and unorganized emissions, etc.) to the environmental receptors. The relationship between the VOCs source component spectrum of each secondary pollution source covered in the primary pollution source and the source component spectrum analyzed by PMF can be expressed as:

[0127] y1=a0+a1x 11 +a2x 12 +…+a p x 1p +ε1;

[0128] y2=a0+a1x 21 +a2x 22 +…+a p x 2p +ε2;

[0129] …

[0130] y n =a0+a1x n1 +a2x n2 +…+a p x np +ε n ;

[0131] In the formula, a0, a1, a2, …, a p represents p+1 parameters to be evaluated, ε i represents the random factor effect on y in the i-th trial i The impact of

[0132] where a0, a1, a2, …, a p By least squares estimation, assuming b0,b1,b2,…,b p They are a0, a1, a2, …, a p Least squares estimation (regression coefficient), then establish the multiple regression equation:

[0133] y1=b0+b1x1+b2x2+…+b p x p

[0134] The regression coefficient in the multiple regression equation is calculated by the following formula:

[0135] B=(X'X) -1 X'Y;

[0136] Where, B=(b0,b1,b2,…,b p ) represents the coefficient matrix of the regression equation, and X' is the transposed matrix of X.

[0137] This scheme uses the stepwise regression method to evaluate the relationship between the VOCs source component spectrum of each secondary pollution source in the primary pollution source and the source component spectrum analyzed by PMF. Specifically, the VOCs source component spectrum of each secondary pollution source is introduced as an emission variable one by one, and when the next emission variable is introduced, the old variables that have been selected for the regression model are tested one by one. The variables with less contribution and significance are deleted from the regression model to ensure that each emission variable in the regression model is meaningful to the source component spectrum analyzed by PMF of the dependent variable. The stepwise regression model is a type of regression model. The steps of the stepwise regression model are as follows Figure 3 As shown:

[0138] The explanation rate parameters of each variable (secondary pollution source) are obtained through the multiple linear regression model to evaluate the contribution of the secondary pollution sources involved in a certain level of pollution source to the VOCs of a certain level of pollution source analyzed from the environmental receptor source. The explanation rate of each dependent variable analyzed by the multiple linear regression equation is shown in Table 7 below:

[0139] Table 7: Example of the degree of interpretation of the VOCs source spectrum of pollution sources on environmental receptors

[0140]

[0141] The regression coefficients are solved by the multivariate linear regression model. The absolute value of each VOCs species in the source component spectrum representing a certain level of pollution source can be calculated by the following formula, and finally the source component spectrum of the reconstructed first-level pollution source can be obtained;

[0142] y i =β0+β1x1+β2x2…+β k x k +ε;

[0143] The correlation analysis between the reconstructed source component spectrum of the primary pollution source and the source component spectrum of the primary pollution source analyzed by environmental receptors can be conducted to evaluate the representativeness of the source component spectrum of the measured industry and the accuracy of its contribution to the source component spectrum analyzed by environmental receptors.

[0144]

[0145] In the formula, r(X sp ,Y sp ) represents the source component spectrum X of a certain reconstructed industry (primary pollution source) sp Source component spectrum Y analyzed by environmental receptors sp The correlation coefficient of Cov(X sp ,Y sp) indicates X sp With Y sp The covariance of Var[X sp ] and Var[Y sp ] represent X sp With Y sp The variance of .

[0146] Comprehensive contribution assessment of risk sources:

[0147] According to the source component spectrum of various pollution sources analyzed by the PMF model, the absolute value of the concentration of species in each pollution source is used to quantify the contribution of pollution sources to environmental risks and health risks:

[0148]

[0149] In the formula, RS k It is the cumulative contribution of risk factor k from a specific pollution source, and the contribution of five risk factors, OFP, SOAFP, HI, LCR and OAV, can be calculated. i is the risk level associated with VOCs species i, V i It is the absolute value of the concentration of VOCs species i in the source component spectrum determined by the PMF model.

[0150] Risk source spatial matching:

[0151] The conditional bivariate probability function (CBPF) analysis is used to verify the risk source allocation results and identify the impact of risk sources under different wind directions and wind speeds, as well as the spatial distribution of the main sources of different risk factors. The calculation is performed using the following formula:

[0152]

[0153] In the formula, CBPF Δθ,Δu represents the probability that the risk source contribution exceeds the threshold x within a specific wind direction interval Δθ and wind speed interval Δu; m Δθ,Δu |C≥x represents the number of instances where the target parameter exceeds the threshold x in the interval Δθ and Δu; m Δθ,Δu Indicates the total number of occurrences in the interval.

[0154] Based on the results calculated by CBPF, the spatial distribution of various risk sources is assigned to the satellite map, and the main sources of various risk factors at the environmental receptor monitoring points can be intuitively found, such as Figure 4 As shown. The method of drawing the spatial distribution map of CPF risk sources: Taking the environmental receptor VOCs monitoring point as the origin, the spatial distribution map of risk sources is drawn covering the entire industrial park; the horizontal and vertical coordinates in the figure represent the four directions of southeast, northwest, and northeast, and the circle (dashed line) represents the wind speed range. The darker the color, the greater the risk contribution of the area to the environmental receptor.

[0155] S5. Conduct risk assessment uncertainty and sensitivity analysis on the risk factor assessment results and risk source quantitative traceability results, and generate risk assessment uncertainty and sensitivity analysis results;

[0156] Risk assessment is inherently uncertain. The uncertainty of risk sources in the calculation process of the theoretical method used in this scheme includes three aspects of uncertainty transfer: (1) the monitoring concentration of volatile organic compound species (Ci); (2) the parameters used in the risk assessment process; (3) the source component spectrum (FPi) analyzed by the PMF model;

[0157] Monte Carlo simulation and Crystal Ball software were used to quantify the uncertainty of the assessment results. The VOCs exposure concentration of environmental receptors was suitable for description by log-normal distribution, and parameters such as OT, REL, EF, AT and ED were modeled by triangular distribution. i The standard deviation comes from the monitoring data of environmental receptors, and the standard deviations of the VOCs source component spectra derived from the risk factors and PMF are set to 20% and 30% of their average values, respectively; each simulation is performed for a preset number of iterations to ensure the convergence and stability of the results; the preset number of iterations is set according to the needs, and is set to 10,000 in this embodiment; 10,000 simulation calculations are performed according to the set parameters to ensure the convergence and stability of the results;

[0158] In addition, sensitivity analysis prioritizes the importance of input assumptions in the risk factor assessment results and identifies sensitive input variables that have the greatest impact on the results, thereby providing guidance for developing the best plan for risk source VOCs control;

[0159] Examples of uncertainty and sensitivity in the health risk assessment of risk sources in industrial parks are shown in Table 8 below:

[0160] Table 8: Examples of uncertainty and sensitivity in health risk assessment of risk sources in industrial parks

[0161]

[0162]

[0163] Among them, A, B, C, D, E and F represent the six types of risk sources analyzed by the comprehensive source analysis method.

[0164] S6. Based on the risk factor assessment results, the quantitative traceability results of risk sources, and the uncertainty and sensitivity analysis results of risk assessment, IMTs are used to develop a comprehensive VOCs emission reduction method for industrial parks;

[0165] Using initial mitigation targets (IMTs), a comprehensive VOCs emission reduction method for industrial parks is formulated, providing quantitative emission reduction targets for VOCs species, pollution sources and important VOCs emission links. IMTs can be regarded as concentrations calculated in reverse based on risk levels, representing the maximum allowable residence concentration of VOCs that pose a risk to environmental receptors without threatening public health.

[0166] IMTs calculation randomly selects risk indicators based on the risk factors that the industrial park itself focuses on;

[0167] The following uses the health risk index and environmental effect risk index as examples to calculate IMTs for the purpose of VOCs emission reduction;

[0168] Taking the health risk index as the purpose of VOCs emission reduction in industrial parks, IMTs is calculated by the following formula:

[0169]

[0170] Where, OAV target , HI target and LCR target Represents the target risk levels of odor pollution, non-carcinogenic risk and carcinogenic risk respectively; OAV target , HI target and LCR target The values ​​of are set to 1, 1, and 1×10 –6 ;

[0171] Taking the environmental effect risk index as the purpose of VOCs emission reduction in industrial parks, IMTs is calculated by the following formula:

[0172]

[0173] In the formula, OFP target and SOAFP target They represent the target risk levels of ozone generation potential and secondary organic aerosol generation potential respectively. Since there are no relevant standards, documents and regulations to define the target values ​​of OFP and SOAPF at this stage, OFP target and SOAFP target The local industrial park environmental management department may set reasonable target values ​​based on actual conditions;

[0174] In calculating the IMT of VOCs species i Afterwards, the observations (Obs) that exceeded these targets during the environmental receptor monitoring period were screened;

[0175] The following formula is used to calculate the emission reduction rate of VOCs species and their sources at the lowest acceptable risk level:

[0176]

[0177] In the formula, R i represents the reduction rate of VOCs species i;

[0178] The share of VOCs emission reduction of each risk source is calculated by the following formula:

[0179] ERS i,j =R i ×P i,j

[0180] Where, ERS i,j represents the emission reduction sharing rate of risk source j to VOCs species i; P i,j It represents the percentage of VOCs species i determined in the PMF model in the source component spectrum of risk source j;

[0181] Examples of VOCs controlled species and their emission reduction sharing rates for major emission sources in industrial parks are shown in Table 9 below:

[0182] Table 9: Examples of VOCs control species and their emission reduction sharing rates for major emission sources in industrial parks

[0183]

[0184]

[0185] Note: A, B, C, D, E and F represent the six types of risk sources analyzed by the comprehensive source analysis method.

[0186] Based on the emission reduction rates of different VOCs species at each risk source obtained by the comprehensive emission reduction method of risk sources, combined with the matching results of the receptor and pollution source component spectra, targeted reduction and control measures are provided for different VOCs emission links of each risk source, that is, reduction and control measures are provided for each secondary pollution source of VOCs emissions from each risk source.

[0187] This embodiment also provides a risk-oriented industrial park VOCs comprehensive traceability system, which adopts the above-mentioned risk-oriented industrial park VOCs comprehensive traceability method.

[0188] The above is only an embodiment of the present invention. The common sense such as the known specific structure and characteristics in the scheme is not described in detail here. The ordinary technicians in the relevant field know all the common technical knowledge in the technical field of the invention before the application date or priority date, can obtain all the existing technologies in the field, and have the ability to apply the conventional experimental means before that date. The ordinary technicians in the relevant field can improve and implement this scheme in combination with their own abilities under the enlightenment given by this application. Some typical known structures or known methods should not become obstacles for ordinary technicians in the relevant field to implement this application. It should be pointed out that for those skilled in the art, without departing from the structure of the present invention, several deformations and improvements can be made, which should also be regarded as the protection scope of the present invention, which will not affect the effect of the implementation of the present invention and the practicality of the patent. The protection scope required by this application shall be based on the content of its claims, and the specific implementation methods and other records in the specification can be used to interpret the content of the claims.

Claims

1. A risk-oriented comprehensive VOCs traceability method for industrial parks, characterized by: It includes the following: S1. Set up monitoring points in the industrial park to be monitored; S2. Set monitoring elements of monitoring points, and use preset monitoring methods to monitor the monitoring points and monitor the monitoring elements; S3. Set risk factors, and conduct risk factor assessment based on the monitored factors to generate risk factor assessment results; S4. Conduct quantitative risk source tracing and generate risk source quantitative tracing results; The quantitative tracing of risk sources includes: PMF source analysis, receptor and pollution source component spectrum matching, risk source comprehensive contribution assessment and risk source spatial matching; S5. Conduct risk assessment uncertainty and sensitivity analysis on the risk factor assessment results and risk source quantitative traceability results, and generate risk assessment uncertainty and sensitivity analysis results; S6. Based on the risk factor assessment results, the quantitative traceability results of risk sources, and the uncertainty and sensitivity analysis results of risk assessment, IMTs are used to develop a comprehensive VOCs emission reduction method for industrial parks; Wherein, the S4 includes: S401, according to the risk factor assessment results, using the PMF model, performing PMF analysis, generating a VOCs source component spectrum, and determining the primary pollution source in the pollution source; the pollution source includes the primary pollution source and the secondary pollution source, the primary pollution source is the industry type; the secondary pollution source is each emission link included in the primary pollution source; S403. Use a multivariate linear regression model to perform secondary matching on the VOCs source component spectrum obtained by PMF analysis and the VOCs source component spectrum obtained by measuring the secondary pollution sources included in each primary pollution source, obtain the relationship between the VOCs source component spectrum of each secondary pollution source and the VOCs source component spectrum analyzed by PMF, and obtain the secondary source risk allocation result.

2. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 1 is characterized in that: The monitoring points include: monitoring points for environmental receptors and monitoring points for pollution sources; Pre-set monitoring methods, including: environmental receptor monitoring methods and pollution source monitoring methods; Environmental receptor monitoring methods, including: Industrial park environmental receptor VOCs monitoring uses gas chromatography-mass spectrometry or equivalent equipment to carry out high-resolution continuous automatic sample collection and analysis for a preset time period; Pollution source monitoring methods, including: sample collection methods and sample analysis methods; Sample collection methods include: monitoring of pollution sources using Suma tanks and DNPH absorption tubes / absorption liquids to collect pollution source discharge outlets and fugitive emission samples; When the polluting enterprise is operating under normal conditions, collect emission samples from the sampling port at the rear end of the exhaust gas treatment facility; Before sample collection, a portable volatile organic compound detector was used to detect the initial concentration of TVOCs in the exhaust gas; When the TVOCs concentration is higher than the preset concentration index, the sample in the Suma tank needs to be diluted a certain number of times before component analysis can be performed on the GC-MS / FID; Organized exhaust gas is collected through a sampling gun. After the organic exhaust gas is filtered and heated, the flow / pressure controller of the Suma tank is used to collect VOCs at a fixed time and quantity. The flow controller of the DNPH absorption liquid is used to collect VOCs at a fixed time and quantity. For samples of organized emissions, a preset number of samples are collected at equal time intervals using a Suma tank and DNPH absorption liquid respectively within a preset time period; If the organized emission is intermittent, a preset number of samples shall be collected at equal time intervals within the preset time period during the high-intensity emission period; The samples of unorganized exhaust gas and VOCs in the central atmosphere of the factory boundary are collected using Suma tanks and DNPH absorption tubes. Under stable production conditions, the Suma tanks are equipped with flow limiting valves and the DNPH absorption tubes are equipped with flow controllers to continuously sample for a preset time. The Suma tanks and DNPH absorption tubes at each point collect a preset number of samples respectively. Sample analysis methods include: samples with too high VOCs concentration collected by the Suma tank need to be diluted a certain multiple in the laboratory before being analyzed on the machine, and the diluted concentration is controlled within the preset concentration range; VOCs samples collected by DNPH adsorption liquid / adsorption tube need to be pretreated before analysis on the instrument; After pretreatment, a preset volume of sample was transferred into a sampling bottle and then analyzed using HPLC; The pretreated samples were analyzed according to the preset method.

3. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 1 is characterized in that: The first-level indicators of the risk factors are divided into two categories: health risk and environmental effect risk; Health risk refers to the risk of pathogenicity, carcinogenicity, and odor pollution from VOCs emitted into the atmosphere from industrial parks posing a health threat to the human body; Environmental effect risk refers to the direct or indirect effects of VOCs emitted into the atmosphere from industrial parks, which may lead to risks such as deterioration of air quality, climate change, and depletion of the ozone layer. The secondary indicators of health risk include: hazard entropy HI, lifetime cancer risk LCR and odor activity value OAV; Secondary indicators of environmental effect risk include: ozone generation potential OFP and secondary organic aerosol generation potential SOAFP; The risk factor assessment is carried out according to the monitored elements, including: OFP represents the maximum contribution of different VOCs species to ozone generation under optimal reaction conditions: OFP i =C i ×MIR i ; In the formula, OFP i represents the ozone generation potential of VOCs species i, dimensionless; C i represents the concentration of VOCs species i during the monitoring period, μg·m-3; MIR i represents the maximum response increment of VOCs species i; SOAFP represents the formation capacity of secondary organic aerosol SOA of different VOCs species: SOAFP i =C i ×Y SOA,i ; In the formula, SOAFP i represents the secondary organic aerosol generation potential of VOCs species i, μg·m-3; C i represents the concentration of VOCs species i during the monitoring period, μg·m-3; YSOA,i represents the SOA yield of VOCs species i; Hazard entropy of specific VOCs species HQ i The cumulative hazard entropy (HI) of all VOCs species is calculated by the following formula: In the formula, HQ i Non-cancer inhalation hazard entropy for specific VOCs species i, HQ i When it is greater than 1, it means that the substance may cause adverse health effects due to inhalation; REL i is the reference exposure level for species i; LCR of carcinogenic risk of specific VOCs species i The cumulative carcinogenic risk LCR of all VOCs species is calculated by the following formula: Where, LCR i is the carcinogenic risk of a specific VOCs species i; DOSEdaily represents the daily inhalation dose of harmful VOCs, in mg·kg-1·d-1, which is calculated by the concentration C of species i. i , unit: μg·m-3, normalized respiration rate BR / BW, L·kg-1·d-1, inhalation absorption factor A, dimensionless, and exposure frequency EF, unit: %; CSF i is the inhalation cancer slope factor for species i, dimensionless; ASFa is the sensitivity coefficient, dimensionless; EDa, AT and FAHa are exposure duration, lifetime exposure time and time at home, respectively; The odor intensity of individual odorants is described by the Odor Activity Value (OAV), which is defined as the ratio between the concentration of a compound and its odor threshold: In the formula, C i is the concentration of VOCs odor substance i, 1×10–6; OT i is the threshold value of odor substance i, 1×10–6; OAV i OAV is the odor activity value of odor substance i, dimensionless; sum It is the sum of the odor activity values ​​of all odor substances evaluated and is dimensionless.

4. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 3 is characterized in that: The PMF model is a receptor model for multivariate factor analysis. The model operation can decompose the receptor sample matrix, i×j dimension, into a contribution matrix i×k and a factor matrix k×j: In the formula, x ij is the contribution of the jth species in the ith sample; g ik is the contribution of the kth source in the ith sample; f kj is the concentration of the jth species in the kth source; e ij is the error score; The factor distribution graph and contribution are obtained by minimizing the objective function Q using the least squares method: Where U ij is the uncertainty for each species in each sample; The VOCs species concentration input into the PMF model includes the VOCs species covered by the selected risk factors; the indicator parameters of the operation results should meet the Qtrue / Qrobust value between 1 and 1.5, the correlation coefficient between the monitored and predicted VOCs species concentrations should be greater than 0.6, and the VOCs species residual normal distribution should be between -3 and 3.

5. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 4 is characterized in that: The multivariate linear regression model is used to evaluate the effects and contributions of multiple independent variables x to the dependent variable y; the VOCs source component spectrum analyzed by the PMF model is used as the dependent variable y, and the VOCs source component spectrum of each secondary pollution source in the primary pollution source of the industrial park is used as the dependent variable x. The multivariate linear regression model evaluates the contribution of each secondary pollution source to the environmental receptor; The relationship between the VOCs source component spectrum of the secondary pollution source of the polluting industry and the source component spectrum analyzed by PMF is expressed as: y1=a0+a1x 11 +a2x 12 +…+a p x 1p +ε1; y2=a0+a1x 21 +a2x 22 +…+a p x 2p +ε2; …… y n =a0+a1x n1 +a2x n2 +…+a p x np +e n ; In the formula, a0, a1, a2, …, a p represents p+1 parameters to be evaluated, ε i represents the random factor effect on y in the i-th trial i The impact of where a0, a1, a2, …, a p By least squares estimation, assuming b0,b1,b2,…,b p They are a0, a1, a2, …, a p Least squares estimation, then establish the multiple regression equation: y1=b0+b1x1+b2x2+…+b p x p ; The regression coefficient in the multiple regression equation is calculated by the following formula: B=(X'X) -1 X'Y; Where, B=(b0,b1,b2,…,b p ) represents the coefficient matrix of the regression equation, and X' is the transposed matrix of X; The regression coefficients are solved by the multivariate linear regression model. The absolute value of each VOCs species in the source component spectrum representing a certain level of pollution source can be calculated by the following formula, and finally the reconstructed industry source component spectrum can be obtained; y i =β0+β1x1+β2x2…+β k x k +e; Conduct correlation analysis between the reconstructed source component spectrum of the primary pollution source and the industry source component spectrum analyzed by environmental receptors to evaluate the representativeness of the measured source component spectrum of the primary pollution source and the accuracy of its contribution to the source component spectrum analyzed by environmental receptors; In the formula, r(X sp ,Y sp ) represents the source component spectrum X of a reconstructed primary pollution source sp Source component spectrum Y analyzed by environmental receptors sp The correlation coefficient of Cov(X sp ,Y sp ) indicates X sp With Y sp The covariance of Var[X sp ] and Var[Y sp ] represent X sp With Y sp The variance of .

6. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 5 is characterized in that: The S4 also includes: S402, according to the VOCs source component spectrum of various pollution sources analyzed by the PMF model, and the primary pollution sources, the absolute value of the concentration of species in each primary pollution source is used to quantify the contribution of the primary pollution source to the environmental effect risk and health risk: In the formula, RS k is the cumulative contribution of risk factor k from a specific pollution source, and the contributions of five risk factors, OFP, SOAFP, HI, LCR and OAV, are calculated; k i is the risk level associated with VOCs species i, V i It is the absolute value of the concentration of VOCs species i in the source component spectrum determined by the PMF model.

7. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 6 is characterized in that: S4 further includes: using CBPF analysis to verify the risk source allocation results, and identifying the impact of risk sources under different wind directions and wind speeds, as well as the spatial distribution of the main sources of different risk factors; wherein the risk source allocation results include: the risk allocation results of the primary pollution sources and secondary sources among the determined pollution sources; wherein the risk factors are VOCs species involving risk factors; Calculate using the following formula: Where, CBPF Δθ,Δu represents the probability that the risk source contribution exceeds the threshold x within a specific wind direction interval Δθ and wind speed interval Δu; m Δθ,Δu |C≥x represents the number of instances where the target parameter exceeds the threshold x in the interval Δθ and Δu; m Δθ,Δu Indicates the total number of occurrences in the interval; Based on the results calculated by CBPF, the spatial distribution of various risk sources is assigned to the satellite map to obtain the source pathways of various risk factors at the environmental receptor monitoring points.

8. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 7 is characterized in that: In S5, Monte Carlo simulation and CrystalBall software were used to quantify the uncertainty of the evaluation results. The VOCs exposure concentration of environmental receptors was described by log-normal distribution, and OT, REL, EF, AT and ED were modeled by triangular distribution. i The standard deviation of comes from the monitoring data of environmental receptors, and the standard deviation of the VOCs source component spectrum derived from the risk assessment parameters and PMF is set to 20% and 30% of its average value, respectively; each simulation is performed for a preset number of iterations; Sensitivity analysis prioritizes the importance of input assumptions in the risk assessment results and identifies sensitive input variables that have the greatest impact on the results.

9. The risk-oriented comprehensive VOCs source tracing method for industrial parks according to claim 8 is characterized in that: The S6 comprises: Taking the health risk index as the purpose of VOCs emission reduction in industrial parks, IMTs is calculated by the following formula: Where, OAV target , HI target and LCR target Represents the target risk levels of odor pollution, non-carcinogenic risk and carcinogenic risk respectively; OAV target , HI target and LCR target The values ​​are set to 1, 1, and 1×10 –6 ; Taking the environmental effect risk index as the purpose of VOCs emission reduction in industrial parks, IMTs is calculated by the following formula: In the formula, OFP target and SOAFP target represent the target risk levels for ozone generation potential and secondary organic aerosol generation potential, respectively; In calculating the IMT of VOCs species i Afterwards, the observations Obs that exceeded these targets during the environmental receptor monitoring period were screened; The following formula is used to calculate the emission reduction rate of VOCs species and their sources at the lowest acceptable risk level: In the formula, R i represents the reduction rate of VOCs species i; The share of VOCs emission reduction of each risk source is calculated by the following formula: ERS i,j =R i ×P i,j Where, ERS i,j represents the emission reduction sharing rate of risk source j to VOCs species i; P i,j It represents the percentage of VOCs species i determined in the PMF model in the source component spectrum of risk source j; Based on the emission reduction rates of different VOCs species at each risk source obtained by the comprehensive emission reduction method of risk sources, combined with the matching results of the receptor and source component spectra of pollution sources, reduction and control measures are provided for each secondary pollution source of VOCs emissions from each risk source.

10. Risk-oriented VOCs comprehensive traceability system for industrial parks, characterized by: A risk-oriented comprehensive VOCs traceability method for industrial parks is adopted as described in any one of claims 1 to 9.

Images