阅读说明：本技术 一种基于高斯模式的大气污染溯源方法 (Atmospheric pollution tracing method based on Gaussian model ) 是由 马方方 么新 杨宏伟 孙如飞 于 2021-07-14 设计创作，主要内容包括：本发明公开了一种基于高斯模式的大气污染溯源方法,基本理论是基于高斯模式气扩散理论,核心计算方法使用JAVA语言实现,包括4个核心处理模块：污染源信息处理模块、气象数据预处理模块、污染源实时数据处理模块和二段式溯源计算模块；本发明通过克里金插值法将相关要素进行空间插值处理,形成标准的GRID数据文件,避免大量的实时运算,提高模型计算效率,降低溯源时间成本,能够实现污染源分钟级溯源追踪,为行政管理水平的提升提供有力支撑。(The invention discloses an atmospheric pollution tracing method based on a Gaussian model, which is based on a Gaussian model gas diffusion theory, is realized by a core calculation method by using a JAVA language and comprises 4 core processing modules: the system comprises a pollution source information processing module, a meteorological data preprocessing module, a pollution source real-time data processing module and a two-section type tracing calculation module; according to the method, the relevant elements are subjected to spatial interpolation processing through a kriging interpolation method to form a standard GRID data file, so that a large amount of real-time operation is avoided, the model calculation efficiency is improved, the source tracing time cost is reduced, the pollution source minute-level source tracing can be realized, and powerful support is provided for the improvement of the administrative management level.)

1. An atmospheric pollution tracing method based on a Gaussian model is characterized in that: the method comprises the following processing steps: (1) and (3) processing pollution source information: uniformly identifying each pollution source in the research area through a pollution source information processing module, and assigning a pollution source ID; the module needs to identify the attribute of the pollution source, including the pollution emission type and the emission information of the main pollutant type; according to the actual research result of the pollution source in the research area, a GRID format file with standard resolution is manufactured and imported into a database; (2) preprocessing meteorological data: preprocessing meteorological elements by using a meteorological data preprocessing module to form a GRID file with standard resolution; (3) and (3) processing real-time data of a pollution source: acquiring pollution source emission information in real time through an online monitoring system by using a pollution source real-time data processing module, and forming a GRID file with standard resolution through a standardized processing flow; (4) two-stage source tracing calculation: the method is divided into fuzzy tracing and accurate tracing, and tracing calculation of different pollution scenes can be realized.

2. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: the pollutant discharge types comprise point source, line source, area source and source, and the main pollutant types comprise VOCs and PM₁₀、SO₂Nitrogen oxide, the emission information include pipe diameter, row's mouth height, emission efficiency, gas temperature.

3. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: and (3) performing spatial interpolation processing on the related elements by a kriging interpolation method in the step (2) to form a standard GRID data file.

4. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: the pollution source emission information in the step (3) comprises a pollution source ID, an emission position, a type and an emission concentration.

5. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: for a certain pollutant, fuzzy tracing is to preliminarily judge a pollution source by using information of current pollutant types, wind speeds and wind directions and pollutant type distribution after pretreatment, qualitatively analyze which pollution sources are responsible for a pollution scene, and define the pollution sources as fuzzy sources.

6. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: the fuzzy tracing specifically comprises the following steps: determining input pollutant Type and concentration set { Type }_i,C_i}; secondly, calculating the overproof degree of each type of pollutant, and extracting the overproof pollutant as a new set { Over }_i}; ③ according to the wind speed (u), the wind direction (dir) and the pollution source { PS_iAttribute information (pollutant type and position) and the like to acquire fuzzy source (VS) possibly contributing to the current pollution scene_i}

{VS_i}＝Θ(u,dir,{Over_i},{PS_i})

Theta is a judgment function, and PS is judged_iWhether or not it is a source of ambiguity, wherein

7. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 1, wherein: the accurate tracing is in fuzzy tracing { VS_iOn the basis, the spatial pollution distribution characteristics caused by the independent emission of each fuzzy source are calculated by utilizing an atmospheric diffusion theory Gaussian mode formula, and thenRealizing vector superposition calculation of multiple fuzzy sources on the space, and then reversely analyzing contribution coefficients of different fuzzy sources by utilizing concentration information of a current pollution scene in cooperation with a global optimization algorithm; when the current scene has various overproof pollutants, the contribution coefficients of the fuzzy sources of various pollutants are analyzed independently, then comprehensive calculation is carried out according to the weight coefficients, and the comprehensive contribution coefficients of different pollution sources to the current pollution scene are analyzed.

8. The atmospheric pollution tracing method based on the Gaussian mode as claimed in claim 6, wherein: the accurate tracing steps are as follows: (v) for fuzzy source set (VS)_iLet us assume the ith fuzzy source VS_iHas an emission efficiency of Q_VS,iCalculating the pollutant concentration at the monitoring point position as C according to the emission type and relevant basic information_VS,iThe calculation formula is as follows

C_vs,i＝F(Q_vs,i,VS_i,u)

F is a pollutant concentration simulation function, and the atmospheric diffusion process of the pollutant in the specific environment is calculated.

② according to { VS_i}、{Q_vs,i}、{C_vs,iAnd contaminant concentration C_iCarrying out vector superposition calculation on different combination parameters of each fuzzy source, and analyzing the most reliable fuzzy source emission combination { Q) through an optimization algorithm_vs,i}_sSo that the following formula holds

Abs, absolute value function; ε is the confidence value, typically 1%.

Thirdly, an optimal C for a certain pollutant can be obtained according to the first step and the second step_vs,iCan calculate { PS }_iAtContribution of the first contaminant { R_i}

Wherein the fuzzy source is not present { VS_iR of the contamination Source in (1)_i＝0。

For multiple pollutants { Over }_iScene of } for different Over_iRepeating the step I and the step II to obtain the enterprise contribution rate of the ith pollutant (R)_j,i}，R_j,iRepresenting the contribution value of the jth enterprise pollution source to the ith overproof pollutant.

Fifthly, distributing the calculation weights (omega) of different pollutants according to the importance degree of different pollutants and expert experience_iAnd calculating to obtain the comprehensive contribution coefficient of each enterprise pollution source in the current pollution scene

R_jAnd expressing the comprehensive contribution coefficient of the jth enterprise pollution source.

Technical Field

The invention belongs to the technical field of air pollution prevention and treatment, and particularly relates to an atmospheric pollution tracing method based on a Gaussian model.

Background

With the development of social economy, the productivity of various industries is continuously improved, and the total social emission of atmospheric pollutants is increased. Although the national requirements for air pollutant emission are more and more strict, due to enterprise cost and process problems, the phenomenon that the emission is not up to standard and even the pollutants are stolen still exists. In the process of pursuing nice life, the attention of urban residents to air quality is higher and higher, the requirements of relevant supervision law enforcement departments on urban atmospheric pollution treatment are stricter and stricter, and quantitative, accurate and rapid source tracing technical methods are the most urgent needs of the law enforcement departments.

In business practice, the number of gas-related emission enterprises is large, the environment supervision capability is relatively weak, the environment supervision mainly depends on traditional means, the informatization level is not high, the requirement of real environment management cannot be met, and a quick and effective air pollution tracing technology is lacked. More automatic pollution source identification technology and pollution emission source management method are urgently needed, so that more refined emission control is developed.

At present, a numerical model method based on atmospheric physics and atmospheric chemistry theories is applied to an atmospheric pollution emission source identification method by domestic and foreign scholars, the method is mainly used for distinguishing the area and the probability of a pollution occurrence source after a pollution event occurs based on meteorological and environmental pollutant concentration monitoring data, although the method can reduce the range of a pollution site model arrangement enterprise of a pollution site arrangement personnel to a certain extent, the pollution source can be found only by large output results of the method and much labor judgment and time investment.

Disclosure of Invention

In order to solve the above problems, the present invention provides an atmospheric pollution tracing method based on a gaussian model, wherein a core calculation method is implemented by using JAVA language, and comprises 4 core processing modules: the system comprises a pollution source information processing module, a meteorological data preprocessing module, a pollution source real-time data processing module and a two-section type tracing calculation module.

The technical scheme adopted by the invention is as follows:

an atmospheric pollution tracing method based on a Gaussian model is characterized in that: the method comprises the following processing steps: (1) and (3) processing pollution source information: uniformly identifying each pollution source in the research area through a pollution source information processing module, and assigning a pollution source ID; the module needs to identify the attribute of the pollution source, including the pollution emission type and the emission information of the main pollutant type; according to the actual research result of the pollution source in the research area, a GRID format file with standard resolution is manufactured and imported into a database; (2) preprocessing meteorological data: preprocessing meteorological elements by using a meteorological data preprocessing module to form a GRID file with standard resolution; (3) and (3) processing real-time data of a pollution source: acquiring pollution source emission information in real time through an online monitoring system by using a pollution source real-time data processing module, and forming a GRID file with standard resolution through a standardized processing flow; (4) two-stage source tracing calculation: the method is divided into fuzzy tracing and accurate tracing, and tracing calculation of different pollution scenes can be realized.

The basic theory of the invention is based on a Gaussian mode gas diffusion theory, Taylor firstly uses a statistical theory to research atmospheric turbulence, Gauss is analyzed through experiments on the basis of the Taylor research to summarize an atmospheric diffusion mode based on the turbulence statistics, namely a Gaussian model, the Gaussian diffusion model is based on a statistical function, the theoretical basis is that in a flat and open plain, the gas field is assumed to be uniform and stable, the atmospheric diffusion firstly spreads along the main wind direction and then spreads to the periphery, and the distribution rule accords with normal distribution. A large number of researches show that the diffusion concentration of the point source continuous emission atmospheric pollutants accords with normal distribution, the application range of the diffusion concentration is generally less than 50 kilometers, and the establishment of a Gaussian diffusion model is based on four assumed conditions: the pollution point source is uniformly, stably and continuously discharged; the atmospheric pollutants follow mass conservation in the diffusion process; the wind direction of the diffusion area is uniform and stable; the contaminant concentration is normally distributed in the horizontal direction and in the vertical direction. The basic formula is as follows:

wherein q is the atmospheric pollutant discharge amount in unit time, mg/s; c is the concentration of the pollutant at any point in the space, mg/m³；σ_yIs the lateral diffusion coefficient, m; sigma_zIs the longitudinal diffusion coefficient, m; u is the average wind speed, m/s.

The pollutant discharge types comprise point source, line source, area source and source, and the main pollutant types comprise VOCs and PM₁₀、SO₂Nitrogen oxide, the emission information include pipe diameter, row's mouth height, emission efficiency, gas temperature.

The set of all the pollution sources investigated in the pollution source information processing module is recorded as { PS_i}，PS_iIndicating the ith source of contamination, each with numerous attributes of type, location, primary contaminant, exhaust, and so forth.

And (3) performing spatial interpolation processing on the related elements by a kriging interpolation method in the step (2) to form a standard GRID data file.

According to the method, the wind speed (u), the wind direction (dir), the cloud Cover (CA), the radiation (rad), the temperature (T) and the humidity (hum) in meteorological elements have important influence on tracing of pollution sources, the relevant elements are subjected to spatial interpolation processing through a kriging interpolation method to form a standard GRID data file, a large amount of real-time operation is avoided, the model calculation efficiency is improved, the time cost of tracing is reduced, the minute-level tracing of the pollution sources can be realized, and the support level of administrative management is improved.

Wherein sigma_y，σ_zThe wind speed (u), the wind direction (dir), the cloud Cover (CA), the radiation (rad), the temperature (T) and the humidity (hum) are in a correlation relationship, selection and judgment are needed according to different conditions, and optimization is needed in actual application.

The pollution source emission information in the step (3) comprises a pollution source ID, an emission position, a type and an emission concentration.

Each plant area needs to be provided with an online pollution source emission monitoring system, and the elements needing to be monitored comprise an emission rate (u)_ps) Discharge efficiency (q)_ps) Contaminant concentration ({ C)_ps}), discharge temperature (T)_ps) For calculating the lifting height (Delta H) of different discharge openings and the pollutant discharge amount q per unit time_ps。

ΔH＝φ(u,u_ps,T_ps)

And (4) carrying out pollutant type identification on the current pollution scene as module input (pollutant type, concentration value and the like) by combining with an offline research result. The module calculates all the overproof pollutants according to the input conditions, and the source tracing analysis needs to be carried out on each overproof pollutant. For a certain pollutant, fuzzy tracing is to preliminarily judge a pollution Source by using information such as a current pollutant type, a wind speed and a wind direction, and pollutant type distribution after pretreatment, qualitatively analyze which pollution sources are responsible for a pollution scene, and define the pollution sources as fuzzy sources (VS).

The fuzzy tracing specifically comprises the following steps: determining input pollutant Type and concentration set { Type }_i,C_i}; secondly, calculating the overproof degree of each type of pollutant, and extracting the overproof pollutant as a new set { Over }_i}; ③ according to the wind speed (u), the wind direction (dir) and the pollution source { PS_iAttribute information (pollutant type and position) and the like to acquire fuzzy source (VS) possibly contributing to the current pollution scene_i}

{VS_i}＝Θ(u,dir,{Over_i},{PS_i})

Theta is a judgment function, and PS is judged_iWhether or not it is a source of ambiguity, wherein

The accurate tracing is in fuzzy tracing { VS_iBased on the theory of atmospheric diffusion, the Gaussian model is usedCalculating spatial pollution distribution characteristics caused by independent emission of each fuzzy source according to a formula, then realizing vector superposition calculation of a plurality of fuzzy sources on the space, and then reversely analyzing contribution coefficients of different fuzzy sources by utilizing concentration information of a current pollution scene in cooperation with a global optimization algorithm; when the current scene has various overproof pollutants, the contribution coefficients of the fuzzy sources of various pollutants are analyzed independently, then comprehensive calculation is carried out according to the weight coefficients, and the comprehensive contribution coefficients of different pollution sources to the current pollution scene are analyzed.

The accurate tracing steps are as follows: (v) for fuzzy source set (VS)_iLet us assume the ith fuzzy source VS_iHas an emission efficiency of Q_VS,iCalculating the pollutant concentration at the monitoring point position as C according to the emission type and relevant basic information_VS,iThe calculation formula is as follows

C_vs,i＝F(Q_vs,i,VS_i,u)

F is a pollutant concentration simulation function, and the atmospheric diffusion process of the pollutant in the specific environment is calculated.

② according to { VS_i}、{Q_vs,i}、{C_vs,iAnd contaminant concentration C_iCarrying out vector superposition calculation on different combination parameters of each fuzzy source, and analyzing the most reliable fuzzy source emission combination { Q) through an optimization algorithm_vs,i}_sSo that the following formula holds

Abs, absolute value function; ε is the confidence value, typically 1%.

Thirdly, an optimal C for a certain pollutant can be obtained according to the first step and the second step_vs,iCan calculate { PS }_iContribution value at first contaminant { R }_i}

Wherein the fuzzy source is not present { VS_iR of the contamination Source in (1)_i＝0。

For multiple pollutants { Over }_iScene of } for different Over_iRepeating the step I and the step II to obtain the enterprise contribution rate of the ith pollutant (R)_j,i}，R_j,iRepresenting the contribution value of the jth enterprise pollution source to the ith overproof pollutant.

Fifthly, distributing the calculation weights (omega) of different pollutants according to the importance degree of different pollutants and expert experience_iAnd calculating to obtain the comprehensive contribution coefficient of each enterprise pollution source in the current pollution scene

R_jAnd expressing the comprehensive contribution coefficient of the jth enterprise pollution source.

Compared with the prior art, the method has the following advantages:

1. the invention provides an atmospheric pollution tracing method based on a Gaussian model, which is based on a Gaussian model gas diffusion theory, is realized by a core calculation method by using a JAVA language and comprises 4 core processing modules: the system comprises a pollution source information processing module, a meteorological data preprocessing module, a pollution source real-time data processing module and a two-section type tracing calculation module.

2. According to the method, the relevant elements are subjected to spatial interpolation processing through a kriging interpolation method to form a standard GRID data file, so that a large amount of real-time operation is avoided, the model calculation efficiency is improved, the source tracing time cost is reduced, the pollution source minute-level source tracing can be realized, and powerful support is provided for the improvement of the administrative management level.

3. The invention adopts a tracing method combining fuzzy tracing and accurate tracing, can realize tracing and tracing calculation of different pollution scenes, can qualitatively analyze which pollution sources are responsible for the pollution scenes, realizes vector superposition calculation of a plurality of fuzzy sources on the space, and then utilizes concentration information of the current pollution scene to reversely analyze contribution coefficients of different fuzzy sources in cooperation with a global optimization algorithm.

Drawings

FIG. 1 is a flow chart of the atmospheric pollution tracing method based on the Gaussian model

Detailed Description

The invention is further described with reference to specific examples.

Example 1

An atmospheric pollution tracing method based on a Gaussian model is characterized in that: the method comprises the following processing steps: (1) and (3) processing pollution source information: uniformly identifying each pollution source in the research area through a pollution source information processing module, and assigning a pollution source ID; the module needs to identify the attribute of the pollution source, including the pollution emission type and the emission information of the main pollutant type; according to the actual research result of the pollution source in the research area, a GRID format file with standard resolution is manufactured and imported into a database; (2) preprocessing meteorological data: preprocessing meteorological elements by using a meteorological data preprocessing module to form a GRID file with standard resolution; (3) and (3) processing real-time data of a pollution source: acquiring pollution source emission information in real time through an online monitoring system by using a pollution source real-time data processing module, and forming a GRID file with standard resolution through a standardized processing flow; (4) two-stage source tracing calculation: the method is divided into fuzzy tracing and accurate tracing, and tracing calculation of different pollution scenes can be realized.

The pollutant discharge types comprise point source, line source, area source and source, and the main pollutant types comprise VOCs and PM₁₀、SO₂Nitrogen oxide, the emission information include pipe diameter, row's mouth height, emission efficiency, gas temperature.

The set of all the pollution sources investigated in the pollution source information processing module is recorded as { PS_i}，PS_iIndicating the ith source of contamination, each with numerous attributes of type, location, primary contaminant, exhaust, and so forth.

And (3) performing spatial interpolation processing on the related elements by a kriging interpolation method in the step (2) to form a standard GRID data file.

According to the method, the wind speed (u), the wind direction (dir), the cloud Cover (CA), the radiation (rad), the temperature (T) and the humidity (hum) in meteorological elements have important influence on tracing of pollution sources, the relevant elements are subjected to spatial interpolation processing through a kriging interpolation method to form a standard GRID data file, a large amount of real-time operation is avoided, the model calculation efficiency is improved, the time cost of tracing is reduced, the minute-level tracing of the pollution sources can be realized, and the support level of administrative management is improved.

Wherein sigma_y，σ_zThe wind speed (u), the wind direction (dir), the cloud Cover (CA), the radiation (rad), the temperature (T) and the humidity (hum) are in a correlation relationship, selection and judgment are needed according to different conditions, and optimization is needed in actual application.

The pollution source emission information in the step (3) comprises a pollution source ID, an emission position, a type and an emission concentration.

Each plant area needs to be provided with an online pollution source emission monitoring system, and the elements needing to be monitored comprise an emission rate (u)_ps) Discharge efficiency (q)_ps) Contaminant concentration ({ C)_ps}), discharge temperature (T)_ps) For calculating the lifting height (Delta H) of different discharge openings and the pollutant discharge amount q per unit time_ps。

ΔH＝φ(u,u_ps,T_ps)

And (4) carrying out pollutant type identification on the current pollution scene as module input (pollutant type, concentration value and the like) by combining with an offline research result. The module calculates all the overproof pollutants according to the input conditions, and the source tracing analysis needs to be carried out on each overproof pollutant. For a certain pollutant, fuzzy tracing is to preliminarily judge a pollution Source by using information such as a current pollutant type, a wind speed and a wind direction, and pollutant type distribution after pretreatment, qualitatively analyze which pollution sources are responsible for a pollution scene, and define the pollution sources as fuzzy sources (VS).

The fuzzy tracing specifically comprises the following steps: determining input pollutant Type and concentration set { Type }_i,C_i}; secondly, calculating the overproof degree of each type of pollutant, and extracting the overproof pollutant as a new set { Over }_i}; ③ according to the wind speed (u), the wind direction (dir) and the pollution source { PS_iAttribute information (pollutant type and position) and the like to acquire fuzzy source (VS) possibly contributing to the current pollution scene_i}

{VS_i}＝Θ(u,dir,{Over_i},{PS_i})

Theta is a judgment function, and PS is judged_iWhether or not it is a source of ambiguity, wherein

The accurate tracing is in fuzzy tracing { VS_iOn the basis, calculating spatial pollution distribution characteristics caused by independent emission of each fuzzy source by using an atmospheric diffusion theory Gaussian mode formula, then realizing vector superposition calculation of a plurality of fuzzy sources on the space, and then reversely analyzing contribution coefficients of different fuzzy sources by using concentration information of a current pollution scene in cooperation with a global optimization algorithm; when the current scene has various overproof pollutants, the contribution coefficients of the fuzzy sources of various pollutants are analyzed independently, then comprehensive calculation is carried out according to the weight coefficients, and the comprehensive contribution coefficients of different pollution sources to the current pollution scene are analyzed.

The accurate tracing steps are as follows: (v) for fuzzy source set (VS)_iLet us assume the ith fuzzy source VS_iHas an emission efficiency of Q_VS,iCalculating the pollutant concentration at the monitoring point position as C according to the emission type and relevant basic information_VS,iThe calculation formula is as follows

C_vs,i＝F(Q_vs,i,VS_i,u)

F is a pollutant concentration simulation function, and the atmospheric diffusion process of the pollutant in the specific environment is calculated.

② according to { VS_i}、{Q_vs,i}、{C_vs,iAnd contaminant concentration C_iCarrying out vector superposition calculation on different combination parameters of each fuzzy source, and analyzing the most reliable fuzzy source emission combination { Q) through an optimization algorithm_vs,i}_sSo that the following formula holds

Abs, absolute value function; ε is the confidence value, typically 1%.

Thirdly, an optimal C for a certain pollutant can be obtained according to the first step and the second step_vs,iCan calculate { PS }_iContribution value at first contaminant { R }_i}

Wherein the fuzzy source is not present { VS_iR of the contamination Source in (1)_i＝0。

For multiple pollutants { Over }_iScene of } for different Over_iRepeating the step I and the step II to obtain the enterprise contribution rate of the ith pollutant (R)_j,i}，R_j,iRepresenting the contribution value of the jth enterprise pollution source to the ith overproof pollutant.

Fifthly, distributing the calculation weights (omega) of different pollutants according to the importance degree of different pollutants and expert experience_iAnd calculating to obtain the comprehensive contribution coefficient of each enterprise pollution source in the current pollution scene

R_jAnd expressing the comprehensive contribution coefficient of the jth enterprise pollution source.

Application example 1

There is heavy industry garden in certain middle part industry city, and there are multiple type enterprises such as cement plant, steel plant, power plant in the garden, and the pollutant that the garden main monitoring has suspended particulate matter (TSP), sulfur dioxide (SO2), Nitrogen Oxide (NOX), Volatile Organic Compounds (VOCS) etc..

(1) Campus pollutant information collection

Collect 12 enterprise unit basic data in the garden, every enterprise can all produce the pollutant of different grade type to the left lower corner of garden is the origin of coordinates, and the positive east direction is X axle positive direction, and the positive north direction is Y axle positive direction and establishes the plane coordinate system, and detailed information statistics is as follows:

TABLE 1

Serial number	Enterprise number	Major pollutants	Type of discharge	Relative position information
					1	QY001	TSP、VOCs	Point source/exhaust port	(0,0)
2	QY002	TSP、SO2、VOCs	Point source/exhaust port	(240,20)
					3	QY003	TSP、NOX、VOCs	Point source/exhaust port	(260,280)
4	QY004	TSP、SO2、NOX、VOCs	Point source/exhaust port	(500,50)
					5	QY005	SO2、NOX、VOCs	Point source/exhaust port	(560,420)
6	QY006	TSP、SO2、VOCs	Point source/exhaust port	(520,900)
					7	QY007	TSP、SO2、NOX、VOCs	Point source/exhaust port	(480,1320)
8	QY008	TSP、SO2、VOCs	Point source/exhaust port	(1100,50)
					9	QY009	TSP、SO2	Point source/exhaust port	(1000,900)
10	QY010	TSP、SO2	Point source/exhaust port	(2500,900)
					11	QY011	TSP、NOX	Point source/exhaust port	(3500,1900)
12	QY012	TSP、SO2、VOCs	Point source/exhaust port	(5650,30)

(2) Meteorological information collection

At 14 days in summer, the weather of the garden is clear and cloudy, and the wind speed is measured to be 2.3m/s at the height of 10 m; the wind direction is southwest; the temperature measured at 1.5 m is 28 deg.C, and the atmospheric stability is relatively stable.

(3) Real-time information acquisition

At 14 days, the online monitoring system displays that values of various indexes monitored by a fixed monitoring point (1000,800) in real time are extracted as follows:

TABLE 2

At 14 days, the exhaust port contaminant types for 12 enterprises were counted as follows:

TABLE 3

(4) Two-stage traceability analysis

TSP and SO according to the measured pollutant index value₂And the overproof phenomenon exists, and the source tracing analysis is carried out aiming at the overproof scenes of two pollutants.

1) Take TSP overproof traceability analysis as an example

Selecting a possible pollution source number set by brushing pollution discharge types of enterprises, wherein the possible pollution source number set comprises QY005

{QY001,QY002,QY003,QY004,QY006,QY007,QY008,QY009,QY0010,QY0011,QY0012}

Secondly, further screening the wind speed, the wind direction and the relative positions of the monitoring points and the pollution sources (enterprise discharge openings), and summarizing the enterprises at the downstream of the wind direction of the monitoring points into wind direction vectors without influencing the current pollution sceneDisplacement vector from some pollutant to monitoring pointHas an inner product of less than 0, i.e.

Taking enterprise QY001 as an example, the wind direction vector isDisplacement vectorComputing the inner product of two vectorsTherefore, it is considered that QY001 has an influence on the excess of TSP in the current scene, that is, QY001 is a fuzzy source.

According to the principle, in the pollutant overproof scene, all fuzzy source sets of TSP overproof are obtained as follows: { VS_i}＝{QY001,QY002,QY003,QY004,QY006,QY008}；

According to the statistics of the monitoring data of TSP in different pollutant source gases, the fuzzy source emission efficiency is extracted as the following table:

TABLE 4

Serial number	Enterprise number	Exhaust port temperature C	TSP discharge efficiency (kg/h)
				1	QY001	120	5300
2	QY002	80	4200
				3	QY003	80	4300
4	QY004	85	2200
				5	QY006	68	3000
6	QY008	110	2800

Thirdly, calculating the concentration C of the pollutant formed by the pollution source QY001 at the monitoring point (1000,800)_QY001Assuming that the height of the monitoring point is Z1.8 m; substituting into formula (4) to obtain concentration value C_QY0011600 in μ g/m³。

Similarly, all fuzzy sources are calculatedWhether the condition is satisfied is judged according to the formula (6). And calculating to obtain:

if the condition is not met, the model calculation parameters need to be adjusted, and the diffusion concentrations of all the fuzzy sources are recalculated; and obtaining the optimized calculated concentration values of all the fuzzy sources through multiple calculations: final acquisition { VS_iThe reasonable calculated concentration of { QY001, QY002, QY003, QY004, QY006, QY008} isCalculating a confidence value:the condition is satisfied.

Fourthly, calculating the contribution values of all pollution sources to the TSP exceeding according to the formula (7)

{R_TSP,j}＝{0.216,0.278,0.33,0.172,0,0.004,0,0,0,0,0,0}

Fifthly, calculating SO in all pollution sources according to the same principle₂The over-standard contribution value:

sixthly, the weight distribution of the two indexes is judged to beω_TSPCalculating to obtain the comprehensive contribution value of all pollution sources to the current pollution scene as 0.4

The above embodiments are only some of the embodiments of the present invention, but the embodiments of the present invention are not limited to the above embodiments, such as various combinations of the embodiments, and any other changes, modifications, substitutions, and combinations which do not depart from the spirit and principle of the present invention should be regarded as equivalent alternatives which are within the scope of the present invention.