CN117252489B

CN117252489B - Ozone pollution source evaluation method, device, medium and electronic equipment

Info

Publication number: CN117252489B
Application number: CN202311543053.3A
Authority: CN
Inventors: 刘雅洁; 淳晶; 周振兴; 王文浚; 王颖; 陶卉婷; 高欢; 张东梅
Original assignee: Hubei Gimbol Environmental Technology Co ltd
Current assignee: Hubei Gimbol Environmental Technology Co ltd
Priority date: 2023-11-20
Filing date: 2023-11-20
Publication date: 2024-03-15
Anticipated expiration: 2043-11-20
Also published as: CN117252489A

Abstract

The application relates to an evaluation method, an evaluation device, a medium and electronic equipment of ozone pollution sources, and relates to the technical field of data analysis, wherein the method comprises the following steps: acquiring a first sensitivity coefficient of a first precursor and a second sensitivity coefficient of a second precursor of each subarea in a preset area within a preset time, and selecting a target subarea with the sum of the coefficients exceeding a first preset value from each subarea; screening first subareas from all target subareas, determining first common industry types in all first subareas within preset time, and determining first contribution quantification values of all first common industry types to ozone pollution; screening the second subareas from the target subareas, determining the second common industry types in the second subareas in a preset time, and determining a second contribution quantification value of the second common industry types to ozone pollution. The method has the effects of determining the distribution condition of ozone pollution sources in specific industries and providing basis for the refined management and control of ozone pollution.

Description

Ozone pollution source evaluation method, device, medium and electronic equipment

Technical Field

The application relates to the technical field of data analysis, in particular to an assessment method, an assessment device, a medium and electronic equipment for ozone pollution sources.

Background

In recent years, along with the continuous promotion of the urban and industrialized processes, ozone pollution in various cities is also remarkable gradually, and particularly in industrial areas in the cities, the ozone pollution mainly comes from production activities of enterprise factories. Ozone, among them, is a gas having strong oxidizing property, and is one of the main components between the stratosphere (about 10-50 km from the ground) and the troposphere (about 0-10 km from the ground) in the atmosphere.

When nitrogen oxides (NOx) and Volatile Organic Compounds (VOCs) mainly generated by production activities of enterprise factories in an industrial area chemically react with oxygen in the atmosphere, ozone is generated, resulting in an increase in ozone concentration, and high concentration ozone affects human health and environment. Therefore, how to analyze the source of ozone pollution and effectively manage and control the ozone pollution is a problem to be solved.

Currently, the means of resolving the source of ozone pollution are generally: based on the model of the atmospheric pollutant source analysis, inputting ozone monitoring data into the model, and finally analyzing to obtain the source of ozone pollution. However, the method can only give the duty ratio of each precursor of the ozone pollution source, and cannot determine the distribution situation of the ozone pollution source in specific industries, so that scientific basis cannot be provided for the precise management and control of the ozone pollution of related departments.

Disclosure of Invention

In order to determine the distribution situation of ozone pollution sources in specific industries, a scientific basis is provided for the refined management and control of ozone pollution of related departments, and the application provides an assessment method, a device, a medium and electronic equipment of the ozone pollution sources.

In a first aspect of the present application, a method for evaluating a source of ozone pollution is provided, comprising:

acquiring a first sensitivity coefficient of a first precursor to ozone generation and a second sensitivity coefficient of a second precursor to ozone generation of each subarea in a preset area within a preset time, determining the sum of the coefficients of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and selecting a target subarea with the sum exceeding a first preset value from each subarea, wherein the first precursor is a volatile organic matter, the second precursor is nitrogen oxide, and the larger the sensitivity coefficient is, the larger the influence on ozone generation is;

screening a first sub-region of each target sub-region for a difference between the first sensitivity coefficient and the second sensitivity coefficient exceeding a second preset value, and determining a first common industry type in which the discharge rate of the first precursor in each first sub-region exceeds a first preset discharge rate within the preset time;

Determining a first contribution quantification value of each first common industry type to ozone pollution in the preset time according to a first enterprise corresponding to each first common industry type in the preset area;

screening a second sub-region of each target sub-region for a difference between the second sensitivity coefficient and the first sensitivity coefficient exceeding a second preset value, and determining a second common industry type in which the discharge rate of the second precursor in each second sub-region exceeds a second preset discharge rate within the preset time;

and determining a second contribution quantification value of each second common industry type to ozone pollution in the preset time according to a second enterprise corresponding to each second common industry type in the preset area.

By adopting the technical scheme, the target subareas are screened out from all subareas in the preset area, the sum of coefficients of the target subareas exceeds a first preset value, which indicates that the target subareas are easier to generate ozone than other subareas, then the first subareas are screened out from all target subareas, the difference of coefficients of the first sensitivity coefficient and the second sensitivity coefficient in the first subareas exceeds a second preset value, which indicates that the first subareas are easier to generate ozone than other target subareas, and the discharged first precursor is easier to generate ozone. The first common industry type is then screened from each first sub-area, i.e. the industry types are not only distributed in a single first sub-area, and the emission rate of the first precursor exceeds the preset emission rate, which means that the emission of the first precursor of the first common industry type is easier for the whole preset area to generate ozone. And finally, determining a first contribution quantification value of the first common industry type to ozone pollution of a preset area in a preset time, so as to determine the specific industry type and distribution condition of the ozone pollution source based on the dimension of the first precursor. Similarly, a second contribution quantification of a second common industry type may be determined, thereby determining a specific industry type and distribution of ozone pollution sources based on the dimensions of the second precursor.

Optionally, before determining that the emission rate of the first precursor in each of the first sub-areas exceeds the first common industry type of the first preset emission rate within the preset time, the method further includes:

acquiring a first contribution value of a first precursor of the preset area to ozone concentration change and a second contribution value of a second precursor to ozone concentration change within the preset time;

if the first contribution value is larger than the second contribution value, reducing a preset initial discharge rate threshold to obtain a first preset discharge rate, and determining the preset initial discharge rate threshold as a second preset discharge rate;

if the first contribution value is smaller than the second contribution value, determining a preset initial emission rate threshold as a first preset emission rate, and reducing the initial emission rate threshold to obtain a second preset emission rate.

By adopting the technical scheme, if the first contribution value is larger than the second contribution value, the fact that the ozone generated by the volatile organic compounds is higher than the ozone generated by the nitrogen oxides is shown, the preset initial emission rate threshold is reduced, and the first preset emission rate is obtained, so that the screening range of the first common industrial type is conveniently enlarged, and the industrial type with larger influence on ozone generation due to the emission of the first precursor is accurately determined; if the first contribution value is smaller than the second contribution value, the fact that the ozone generated by the volatile organic compounds is lower than the ozone generated by the nitrogen oxides is indicated, the preset initial emission rate threshold is reduced, and the second preset emission rate is obtained, so that the screening range of the second common industrial type is conveniently enlarged, and the industrial type with larger influence on ozone generation due to the emission of the second precursor is accurately determined.

Optionally, the determining, according to the first enterprises corresponding to each of the first common industry types in the preset area, a first contribution quantization value of each of the first common industry types to ozone pollution in the preset time specifically includes:

counting the initial discharge rate of the first precursor of each first enterprise corresponding to each first common industry type in the preset area;

obtaining emission influence values of each first common industry type based on first sensitivity coefficients of subareas corresponding to each first enterprise and corresponding first precursor initial emission rates;

and determining a first contribution quantification value of the corresponding first common industry type to the first contribution value within the preset time according to the ratio of the emission influence value of each first common industry type to the sum of all emission influence values.

By adopting the technical scheme, according to the first sensitivity coefficient and the initial emission rate of the first precursor of the subarea of the first enterprises belonging to the first common industry type, the emission influence values of all the first enterprises corresponding to each first common industry type are determined, so that the influence degree of the emission of the first precursor of each first common industry type on ozone generation is represented. And then determining the distribution condition of the first contribution value in the ozone pollution source in each first common industry type in the preset area according to the ratio of the emission influence value of each first common industry type to the sum of the emission influence values of all the first common industry types.

Optionally, the obtaining the emission impact value of each of the first common industry types based on the first sensitivity coefficient of the sub-region corresponding to each of the first enterprises and the corresponding first precursor initial emission rate specifically includes:

multiplying the first sensitivity coefficient of the subarea corresponding to each first enterprise by the corresponding first precursor initial discharge rate to obtain an initial discharge influence value;

counting the number of the first enterprises in the same subarea, and sorting the subareas with the first enterprises according to the number to obtain a sorted subarea set, wherein the more the number is, the more the sorting is;

selecting a preset number of specific subareas from the ordered subarea set according to the sequence from front to back, and matching first correction coefficients corresponding to the specific subareas, wherein the first correction coefficients are larger than 1, and the more the ordering is, the larger the corresponding first correction coefficients are;

and multiplying the initial emission influence value of the first enterprise in each specific subarea by a corresponding first correction coefficient to obtain corrected influence values, and summing the corrected influence values of the first enterprises belonging to the same first common industry type in the specific subarea and the initial emission influence values of other first enterprises belonging to the same first common industry type outside the specific subarea to obtain emission influence values of each first common industry type.

By adopting the technical scheme, the more the number is, the more first enterprises belonging to the first common industry type in the corresponding subarea are, and further the greater the influence of the subarea on the ozone generation of the whole preset area is. And selecting a preset number of specific subareas from the ordered subarea set, and determining a corresponding first correction coefficient, so that the initial emission influence value of a first enterprise in the specific subarea is calibrated, and the finally determined emission influence value of a first common industry type is more objective and accurate.

Optionally, after determining the first contribution quantification value of each of the first common industrial types to ozone pollution in the preset time, the method further includes:

counting the number of the subareas to which each first common industry type belongs;

and matching the second correction coefficients corresponding to the numbers from a preset coefficient matching table, multiplying the second correction coefficients by the first contribution quantized values of the corresponding first common industry types to obtain final first contribution quantized values of the corresponding first common industry types, wherein the second correction coefficients are larger than 1, and the larger the numbers are, the larger the corresponding second correction coefficients are.

By adopting the technical scheme, the more the number is, the more the subareas in the preset area are provided with the first common industry type, which means that the greater the influence of the emission of the first precursor of the first common industry type on the ozone generation in the preset area is, the greater the influence of the emission of the first precursor of the first common industry type on the ozone generation in the whole preset area is, the greater the corresponding second correction coefficient is, and the first contribution quantization value of the first common industry type is multiplied by the corresponding second correction coefficient, so that the first contribution quantization value of the first common industry type is more accurate.

selecting a first common industry type with a first contribution quantization value exceeding a preset quantization value as a main contribution industry type, and determining a main contribution subarea with the largest enterprise existence number corresponding to each main contribution industry type from all subareas;

when the wind speed is greater than a wind speed threshold value, judging whether a densely populated area exists in a wind direction downstream area of the preset area;

if so, judging the diffusion duration of the first precursor discharged by each main contribution subarea to the densely populated area;

Counting a first contribution quantity of a main contribution subarea of which the diffusion duration does not exceed the preset time, and if the first contribution quantity exceeds a first quantity threshold value, determining the preset area as an ozone key monitoring area;

and if the first contribution quantity does not exceed a first quantity threshold value, adjusting the preset time to be the diffusion duration so as to redetermine whether the preset area is an ozone key monitoring area.

By adopting the technical scheme, main contribution subareas with main contribution of the first precursor to ozone pollution are screened from the preset areas, if the wind downstream area is a population dense area, the main contribution subareas are extremely easy to cause ozone pollution to the population dense area, then the first contribution quantity of the main contribution subareas which can diffuse into the population dense area in preset time is screened from the main contribution subareas, if the first contribution quantity exceeds a first quantity threshold value, the first precursor generated by the main contribution subareas is indicated to diffuse into the population dense area, ozone is generated to cause ozone pollution to the population dense area, and then the preset area is determined as an ozone key monitoring area; if the first quantity threshold value is not exceeded, the preset time is adjusted to be the largest diffusion duration in the diffusion durations corresponding to the main contribution subregions, and whether the preset region is determined to be the ozone key monitoring region or not is accurately determined according to the first contribution quantity of the main contribution subregions which can be diffused to the population dense region.

Optionally, after determining that the enterprises corresponding to each main contribution industry type exist in the main contribution subareas with the largest number from the subareas, the method further includes:

when the wind speed is greater than a wind speed threshold value, determining a correlation subarea combination with a wind direction upstream-downstream relation from each main contribution subarea, wherein the correlation subarea combination comprises an upstream subarea and a downstream subarea;

determining a diffusion duration of the first precursor discharged from the upstream sub-region to the downstream sub-region, if the diffusion duration exceeds the preset time, increasing a first contribution quantification value of the main contribution industry type corresponding to the upstream sub-region, and keeping the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-region unchanged;

if the diffusion duration does not exceed the preset time, when the first sensitivity coefficient of the downstream sub-area is greater than the first sensitivity coefficient of the upstream sub-area, reducing the first contribution quantization value of the main contribution industry type corresponding to the downstream sub-area, and increasing the first contribution quantization value of the main contribution industry type corresponding to the upstream sub-area, wherein the reduction amount corresponding to the downstream sub-area is greater than the increase amount corresponding to the upstream sub-area;

And when the first sensitivity coefficient of the downstream sub-region is smaller than that of the upstream sub-region, reducing the first contribution quantization value of the main contribution industry type corresponding to the downstream sub-region, and increasing the first contribution quantization value of the main contribution industry type corresponding to the upstream sub-region, wherein the reduction amount corresponding to the downstream sub-region is smaller than that corresponding to the upstream sub-region.

By adopting the technical scheme, if the diffusion duration of the first precursor of the upstream sub-area to the downstream sub-area in the associated sub-area combination does not exceed the preset time, the first precursor of the upstream sub-area may be diffused into the downstream sub-area to generate ozone, so that the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-area is higher, and therefore, the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-area is reduced, and the first contribution quantification value of the main contribution industry type corresponding to the upstream sub-area is increased. If the diffusion duration exceeds the preset time, the first precursor in the upstream sub-area is not diffused into the downstream sub-area to generate ozone before the preset time is over, the first contribution quantification value of the main contribution industry type corresponding to the upstream sub-area is increased, and the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-area is kept unchanged.

In a second aspect of the present application, there is provided an assessment device for ozone pollution sources, comprising in particular:

the target area determining module is used for obtaining a first sensitivity coefficient of a first precursor to ozone generation and a second sensitivity coefficient of a second precursor to ozone generation in preset time of each subarea in a preset area, determining the sum of the coefficients of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and selecting a target subarea with the sum of the coefficients exceeding a first preset value from each subarea, wherein the first precursor is a volatile organic matter, the second precursor is a nitrogen oxide, and the larger the sensitivity coefficient is, the larger the influence on ozone generation is;

a first industry determination module, configured to screen a first sub-area in which a difference between coefficients obtained by subtracting a second sensitivity coefficient from a first sensitivity coefficient in each target sub-area exceeds a second preset value, and determine a first common industry type in which an emission rate of the first precursor in each first sub-area exceeds a first preset emission rate within the preset time;

the first pollution evaluation module is used for determining a first contribution quantification value of each first common industry type to ozone pollution in the preset time according to a first enterprise corresponding to each first common industry type in the preset area;

A second industry determination module, configured to screen a second sub-area in which a difference between coefficients obtained by subtracting the first sensitivity coefficient from the second sensitivity coefficient exceeds a second preset value, from the target sub-areas, and determine a second common industry type in which the emission rate of the second precursor in each of the second sub-areas exceeds a second preset emission rate within the preset time;

and the second pollution evaluation module is used for determining a second contribution quantification value of each second common industry type to ozone pollution in the preset time according to a second enterprise corresponding to each second common industry type in the preset area.

By adopting the technical scheme, the target area determining module determines target subareas from all subareas in the preset area, the first industry determining module screens first subareas with the difference between the coefficients of the first sensitivity coefficient and the second sensitivity coefficient exceeding a second preset value from all target subareas, determines first common industry types from the first subareas, and then determines first contribution quantization values of all the first common industry types to ozone pollution according to first enterprises corresponding to all the first common industry types in the preset area by the first pollution evaluation module. In addition, the second common industry type is determined by the second industry determination module, and finally the second pollution evaluation module determines a second contribution quantification value of each second common industry type to ozone pollution. Thereby determining the specific industry of the main source of ozone pollution and providing scientific basis for the refined management and control of ozone pollution in related departments.

In summary, the present application includes at least one of the following beneficial technical effects: screening the discharged first precursor from the target sub-area that is more prone to ozone generation is easier to generate the first sub-area of ozone, and then determining contributions of each first common industry type to ozone pollution based on determining the first common industry type from the first sub-area that is more representative of ozone generation effects. Similarly, a second sub-region in which the second precursor is more prone to ozone generation is screened from the target sub-region, and then the contribution of each second common industry type to ozone pollution is determined based on the determination of a first common industry type from the second sub-region that is more representative of ozone generation effects. Thereby determining the distribution of ozone pollution sources in a specific industry.

Drawings

FIG. 1 is a schematic flow chart of an evaluation method of ozone pollution sources according to an embodiment of the present application;

FIG. 2 is a schematic diagram of an evaluation flow of ozone pollution sources according to an embodiment of the present application;

FIG. 3 is a flow chart of another method for evaluating ozone pollution sources according to an embodiment of the present application;

FIG. 4 is a schematic structural diagram of an apparatus for evaluating ozone pollution sources according to an embodiment of the present application;

Fig. 5 is a schematic structural diagram of another ozone pollution source evaluation device according to an embodiment of the present application.

Reference numerals illustrate: 11. a target area determination module; 12. a first industry determination module; 13. a first pollution assessment module; 14. a second industry determination module; 15. a second pollution evaluation module; 16. an emission threshold adjustment module; 17. a contribution quantization correction module; 18. a key region determining module; 19. and an industry quantization adjustment module.

Detailed Description

In order to make the technical solutions in the present specification better understood by those skilled in the art, the technical solutions in the embodiments of the present specification will be clearly and completely described below with reference to the drawings in the embodiments of the present specification, and it is obvious that the described embodiments are only some embodiments of the present application, but not all embodiments.

In the description of embodiments of the present application, words such as "exemplary," "such as" or "for example" are used to mean serving as an example, instance, or illustration. Any embodiment or design described herein as "illustrative," "such as" or "for example" is not necessarily to be construed as preferred or advantageous over other embodiments or designs. Rather, the use of words such as "illustratively," "such as" or "for example," etc., is intended to present related concepts in a concrete fashion.

Referring to fig. 1, an embodiment of the present application discloses a schematic flow chart of an evaluation method of ozone pollution sources, which can be implemented by a computer program or can be run on an evaluation device of ozone pollution sources based on von neumann system. The computer program can be integrated in an application or can be run as a stand-alone tool class application, and specifically comprises:

s101: acquiring a first sensitivity coefficient of a first precursor to ozone generation and a second sensitivity coefficient of a second precursor to ozone generation of each subarea in a preset area within a preset time, determining the sum of the coefficients of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and selecting a target subarea with the sum of the coefficients exceeding a first preset value from each subarea.

Specifically, in the embodiment of the present application, the preset area may be an urban industrial area where ozone pollution sources need to be evaluated, and the subareas may be individual industrial parks in the industrial area. The preset time can be one week or one month. Wherein, the precursor refers to the substance of the previous stage of the product in a certain reaction process or refers to the substance of the previous stage of a certain metabolic intermediate. In the ozone generation process, the first precursor volatile organic compound and the second precursor nitrogen oxide generate ozone through photochemical reaction. The greater the susceptibility coefficient, the greater the impact on ozone generation, e.g., the greater the susceptibility coefficient of the first precursor in a region, the more readily ozone is generated by volatile organic compounds emitted from that region, resulting in ozone pollution.

One possible way to obtain the first and second sensitivity coefficients for each sub-region is: the first and second sensitivity coefficients are determined by a DDM model, which is a mathematical model that can be used to study chemical dynamics and environmental impact in an ecosystem. The model may calculate its impact on the ecosystem based on the input, output and conversion processes of chemicals in the ecosystem. For example, the DDM model determines a first sensitivity coefficient of a first precursor for a subregion based on the emission rate of the first precursor for the subregion for a preset time, the ozone generation amount for the subregion for a preset time. This is the prior art and will not be described in detail here. In other embodiments, the HDDM model may also be used to determine the first and second sensitivity coefficients. The high-order direct decoupling HDDM model is a method to study chemical dynamics and environmental impact in ecosystems. By using a higher-order taylor expansion to represent the response of the simulated concentration c to changes in the vector λ, a change in the sensitivity of c to λj due to a change in λi can be derived. This is also prior art and will not be described in detail here.

After the first sensitivity coefficient and the second sensitivity coefficient of each sub-area are determined, the first sensitivity coefficient and the second sensitivity coefficient of each sub-area are added to obtain a coefficient sum, then the coefficient sum corresponding to each sub-area is compared with a first preset value, a sub-area with the coefficient sum exceeding the first preset value is determined as a target sub-area, and the first preset value is a critical value for measuring the influence degree of a precursor on ozone generation. The larger the sum of the coefficients, the more likely the precursor of the corresponding subregion will generate ozone.

In other embodiments, step S101 further includes: acquiring a first contribution value of a first precursor to ozone concentration change and a second contribution value of a second precursor to ozone concentration change of a preset region within a preset time;

Specifically, the ozone concentration variation of a preset area in preset time is obtained by accessing an environment monitoring APP or a website, then the emission of a first precursor and a second precursor in the preset time, the weather conditions in the preset time and the preset chemical reaction kinetic parameters are input into a preset SARSM-DM model, and finally a first contribution value of the first precursor to the ozone concentration variation and a second contribution value of the second precursor to the ozone concentration variation are obtained. The SARSM-DM model is an atmospheric air chemical transmission model and can be used for simulating the formation and destruction process of ozone in the atmosphere. The model considers the contribution of various precursors to ozone generation, including Volatile Organic Compounds (VOCs) and nitrogen oxides (NOx). This is the prior art and will not be described in detail here.

Comparing the first contribution value with the second contribution value, if the first contribution value is larger than the second contribution value, the method indicates that the ozone generated by the volatile organic compounds is higher than the ozone generated by the nitrogen oxides, then reducing a preset initial emission rate threshold to obtain a first preset emission rate, and in addition, keeping the preset initial emission rate threshold unchanged and determining the preset initial emission rate as a second preset emission rate, so that the screening range of the first common industrial type is conveniently enlarged, and further, the industrial type with larger influence on ozone generation by the emission of the first precursor is accurately determined. If the first contribution value is smaller than the second contribution value, the fact that the ozone generated by the volatile organic compounds is lower than the ozone generated by the nitrogen oxides is indicated, the preset initial emission rate threshold is reduced, and the second preset emission rate is obtained, so that the screening range of the second common industrial type is conveniently enlarged, and the industrial type with larger influence on ozone generation due to the emission of the second precursor is accurately determined. In addition, the preset initial emission rate threshold is maintained and determined as the first preset emission rate. If the first contribution value is equal to the second contribution value, the initial emission rate threshold is maintained unchanged and is determined as both the first preset emission rate and the second preset emission rate.

The reduction of the initial emission rate threshold is obtained by matching from a preset reduction matching table according to the difference between the first contribution value and the second contribution value.

S102: screening the first subareas in which the difference between the first sensitivity coefficient and the second sensitivity coefficient exceeds a second preset value from the target subareas, and determining a first common industry type in which the discharge rate of the first precursor in each first subarea exceeds a first preset discharge rate in a preset time.

Specifically, after the target subareas are determined, the first sensitivity coefficient and the second sensitivity coefficient of each target subarea are subjected to difference to obtain coefficient difference, the coefficient difference is compared with a second preset value, and the second preset value is a positive integer and is used for measuring the critical value of the difference between the first sensitivity coefficient and the second sensitivity coefficient. If the difference between the coefficients exceeds a second preset value, indicating that the volatile organic compounds emitted by the corresponding target subregion are more prone to ozone generation than nitrogen oxides within the preset time, then the corresponding target subregion is determined to be the first subregion.

Next, a first common industry type is determined in which the discharge rate of the first precursor exceeds a first predetermined discharge rate in each first sub-zone for a predetermined time, one possible manner of determination is: the emission rate of the first precursor within a preset time measured by the GC-FID on-line monitoring system of the volatile organic compounds within each factory enterprise in the first sub-area is received. Extracting main camping keywords from official websites of each factory enterprise through a word2vec model, classifying the factory enterprises of the first subareas by utilizing the trained industry classification model according to the main camping keywords of the factory enterprises, adding the discharge rates of first precursors of the factory enterprises of the same industry type, and finally determining the discharge rate of the first precursors of each industry type in a single first subarea within preset time. The industry classification model may employ an XGBoost model.

Further, the high-emission industrial type in which the emission rate of the first precursor in each first sub-area exceeds the first preset emission rate is selected, that is, the industrial type in which the first precursor in each first sub-area has a large influence on ozone generation is determined, and then the high-emission industrial type in each first sub-area is determined, which indicates that the high-emission industrial type may be the industrial type that is easy to produce ozone, and is taken as the first common industrial type, that is, the industrial type in which the first precursor in the preset area has a large influence on ozone generation is determined.

In other embodiments of the application, the determination of the industry type in which the number of occurrences exceeds the preset number in all the first sub-areas may also be referred to as the first common industry type. The number of times of the preset may be 1, and in other embodiments, the number of the first sub-areas may be determined according to the number of the first sub-areas, where the number of the first sub-areas is greater, the number of times of the preset is greater. In the embodiment of the present application, the discharge rate refers to the amount of the precursor discharged per unit time during a certain period.

S103: and determining a first contribution quantification value of each first common industry type to ozone pollution in a preset time according to the first enterprises corresponding to each first common industry type in the preset area.

In one implementation, the first precursor initial discharge rate of each first enterprise corresponding to each first common industry type in the preset area is counted;

obtaining emission influence values of each first common industry type based on the first sensitivity coefficient of the subarea corresponding to each first enterprise and the corresponding first precursor initial emission rate;

and determining a first contribution quantification value of the corresponding first common industry type to the first contribution value in a preset time according to the ratio of the emission influence value of each first common industry type to the sum of all emission influence values.

Specifically, after the first common industry type is determined, the initial emission rate of the first precursor of each first enterprise, namely, the emission rate of the volatile organic compound, is obtained through the GC-FID of each first enterprise belonging to the first industry type in the preset area. And then, acquiring a first sensitivity coefficient of a subarea to which each first enterprise belongs, multiplying the first sensitivity coefficient by the initial emission rate of the first precursor of the first enterprise to obtain a plurality of product results, and finally accumulating the product results corresponding to the first enterprises belonging to the same first common industry type to obtain an emission influence value of each first common industry type, wherein the larger the emission influence value is, the larger the influence of the corresponding first common industry type on ozone generation is, and the more ozone is easy to generate.

In other embodiments, determining emissions impact values for each first common industry type may also be performed as follows: multiplying the first sensitivity coefficient of the subarea corresponding to each first enterprise by the initial discharge rate of the corresponding first precursor to obtain an initial discharge influence value;

counting the number of first enterprises in the same subarea, and sorting the subareas with the first enterprises according to the number to obtain a sorted subarea set, wherein the more the number is, the more the sorting is;

selecting a preset number of specific subareas from the ordered subarea set according to the sequence from front to back, and matching first correction coefficients corresponding to the specific subareas, wherein the first correction coefficients are larger than 1, and the corresponding first correction coefficients are larger when the ordering is more forward;

the initial emission influence values of the first enterprises in the specific subareas are multiplied by the corresponding first correction coefficients to obtain corrected influence values, and the corrected influence values of the first enterprises belonging to the same first common industry type in the specific subareas and the initial emission influence values of other first enterprises belonging to the same first common industry type outside the specific subareas are summed to obtain emission influence values of each first common industry type.

Specifically, after the initial discharge rate of the first precursor of each first enterprise is determined, multiplying the first sensitivity coefficient of the subarea corresponding to the first enterprise by the corresponding initial discharge rate of the first precursor to obtain an initial discharge influence value, so that the influence degree of the discharge of the first precursor of each first enterprise on ozone generation in a preset time is primarily determined. And then counting the number of the first enterprises in each subarea in the preset area, wherein the more the number is, the more the first enterprises belonging to the first common industry type in the corresponding subarea are, and further, the greater the influence of the subarea on the ozone generation of the whole preset area is. And then sequencing the subareas with the first enterprise according to the number to obtain a sequenced subarea set. The more the number, the more forward the corresponding subregions are, and the greater the ozone generation effect on the whole preset region is. And selecting a preset number of specific subareas from the ordered subarea set according to the sequence from front to back, wherein the preset number can be determined according to the number of subareas in the ordered subarea set, and the larger the number of subareas is, the larger the preset number is.

And finally, matching corresponding first correction coefficients from a preset correction coefficient matching table according to the sorting value of each specific subarea, wherein the correction coefficient matching table comprises sorting values and corresponding first correction coefficients, and the sorting values are smaller and the corresponding first correction coefficients are larger when the sorting is closer. For example, if the specific sub-region a is in the first bit in the sorted sub-region set, its sorting value is 1, and then the first correction coefficient corresponding to the specific sub-region a is the largest.

Further, the initial emission influence value of each first enterprise in all the specific subareas is multiplied by a first correction coefficient corresponding to the specific subareas to obtain corrected influence values, and then the corrected influence values of all the first enterprises belonging to the same first common industry type in the specific subareas are summed to obtain a first summation result.

And similarly, summing initial emission influence values of other first enterprises corresponding to the same first common industry type in the subareas outside the specific subareas to obtain a second summation result, and adding the first summation result and the second summation result to obtain accurate and objective emission influence values of each first common industry type.

Further, after determining the emission influence value of each first common industry type, summing the emission influence values to obtain a sum of all emission influence values, dividing the emission influence value of each first common industry type by the sum of all emission influence values, and multiplying the obtained ratio by the first contribution value to obtain a first contribution quantized value of the first common industry type, namely quantizing the contribution condition of the first common industry type to ozone generation in a preset area. Thereby realizing the evaluation of the industrial source of ozone pollution, and facilitating relevant parts to make targeted industry environment-friendly rectifying response, and particularly referring to fig. 2.

S104: screening the second sub-areas of each target sub-area for a difference between the second sensitivity coefficient and the first sensitivity coefficient exceeding a second preset value, and determining a second common industry type in which the discharge rate of the second precursor in each second sub-area exceeds a second preset discharge rate within a preset time.

Specifically, the difference between the coefficients of the second sensitivity coefficient and the first sensitivity coefficient exceeds a second preset value, which indicates that the influence degree of the nitrogen oxides in the corresponding target subarea on the generated ozone is far higher than that of the volatile organic compounds, and further indicates that the nitrogen oxides discharged by the corresponding target subarea in the preset time are easier to generate ozone than the volatile organic compounds. After the second sub-areas are determined, a process of obtaining the second common industry type in each second sub-area is similar to the step S102, which is not described herein.

S105: and determining a second contribution quantification value of each second common industry type to ozone pollution in a preset time according to a second enterprise corresponding to each second common industry type in the preset area.

Specifically, after the second common industry type is determined, obtaining the initial discharge rate of the second precursor of each second enterprise belonging to the second common industry type in a preset area range, then obtaining the second sensitivity coefficient of the subarea to which each second enterprise belongs, multiplying the second sensitivity coefficient by the initial discharge rate of the second precursor of the second enterprise to obtain a plurality of product results, and finally accumulating the product results corresponding to the second enterprises belonging to the same second common industry type to obtain the discharge influence value of each second common industry type, wherein the larger the discharge influence value is, the larger the influence of the corresponding first common industry type on ozone generation is, and the ozone is easier to generate.

Further, after the emission influence value of each second common industry type is determined, summing the emission influence values to obtain a sum of all emission influence values, dividing the emission influence value of each second common industry type by the sum of all emission influence values, and multiplying the obtained ratio by the second contribution value to obtain a second contribution quantized value of the second common industry type, namely quantizing the contribution condition of the second common industry type to ozone generation in a preset area. Thereby realizing the evaluation of the industrial sources of ozone pollution and determining the distribution situation of the ozone pollution sources in specific industries. The details can be seen in step S103, which is not described herein.

In another embodiment, step S102 further includes: determining a pending industry type in which the discharge rate of the first precursor in each first sub-zone does not exceed a first preset discharge rate within a preset time;

and determining the types of the pending industries appearing in the first subareas exceeding the second preset number as the types of the pending common industries, and adding the types of the pending common industries into an ozone pollution to-be-monitored list corresponding to the preset areas.

Specifically, if the emission rate of the first precursor of the industry type in the first sub-area does not exceed the first preset emission rate, indicating that the emission of the first precursor of the industry type in the first sub-area has less influence on ozone generation, the industry type is determined as a pending industry type. And counting the number of the first subareas to which each of the pending industry types belongs, if the number of the first subareas exceeds the second preset number, the fact that although the pending industry types are distributed in more first subareas and have a high probability of having a great influence on ozone pollution although the emission rate of the first precursor is low in preset time is indicated, determining the pending industry type as the pending common industry type, and adding the determined industry type to a corresponding ozone pollution to-be-monitored list in the preset area, so that relevant departments can monitor the emission focus of the first precursor of the pending common industry type, and respond in advance to avoid the fact that the pending common industry type has a great influence on ozone generation later.

In yet another embodiment, after step S103, the method further includes: selecting a first common industry type with a first contribution quantization value exceeding a preset quantization value as a main contribution industry type, and determining a main contribution subarea with the largest enterprise existence number corresponding to each main contribution industry type from all subareas;

if the first precursor is present, judging the diffusion duration of the first precursor discharged by each main contribution subarea to the densely populated area;

counting the first contribution quantity of the main contribution subareas of which the diffusion duration does not exceed the preset time, and if the first contribution quantity exceeds a first quantity threshold value, determining the preset area as an ozone key monitoring area;

and if the first contribution quantity does not exceed the first quantity threshold value, adjusting the preset time to be the maximum diffusion duration so as to redetermine whether the preset area is an ozone key monitoring area.

Specifically, after determining the first contribution quantization value corresponding to each first common industry type, if the first contribution quantization value exceeds a preset quantization value, determining the corresponding first common industry type as a main contribution industry type, which indicates that the main contribution industry type in the first common industry type for discharging the first precursor has the greatest influence on ozone pollution. Then determining that the enterprises in the subareas have the most main contribution industry types, and determining the corresponding subareas as main contribution subareas of the main contribution industry types. Then when the current wind speed is larger than a wind speed threshold value, the diffusion of the first precursor is obvious, and whether a population dense area exists in a wind direction downstream area of the preset area is judged, wherein one feasible judgment mode is as follows: the population density is obtained by dividing the population quantity by the area of each area in the downstream area of the wind direction through the population census website, and if the population density exceeds the density threshold value, the corresponding area is determined to be a densely populated area.

If the population-dense region exists, determining a corresponding diffusion duration according to the distance from each main contribution subarea to the population-dense region and the wind speed, if the diffusion duration of the main contribution subarea does not exceed a preset time, indicating that the first precursor can diffuse to the population-dense region before the end of the statistical duration of the preset time, generating ozone pollution, counting the first contribution quantity of the main contribution subareas, wherein the first contribution quantity exceeds a first quantity threshold value, indicating that the first precursor generated by the existence of more main contribution subareas in the preset region can diffuse to the population-dense region in the preset time, generating serious ozone pollution, and determining the preset region as an ozone key monitoring region.

If the first contribution amount does not exceed the first amount threshold, it is indicated that the first precursor generated by the less main contribution subareas in the preset area diffuses into the population-dense area in the preset time, but it is not necessarily indicated that the preset area does not generate serious ozone pollution in the population-dense area, further, the preset time is prolonged, that is, the time range of ozone pollution source evaluation of the preset area is adjusted, the maximum diffusion duration in the diffusion duration corresponding to each main contribution subarea is selected, the first contribution quantification value of each first common industrial type is redetermined, whether the first contribution amount exceeds the first amount threshold is determined again, and further whether the preset area is an ozone key monitoring area is accurately determined.

In other embodiments, if the first number of contributions does not exceed the first number threshold, then a main contribution sub-region corresponding to the second precursor is determined in a similar manner, and when the second number of contributions of the main contribution sub-region exceeds the second number threshold, the preset region is determined to be an ozone key monitoring region. The second contributing amount is the amount of the discharged second precursor that would spread to the main contributing sub-area of the densely populated area.

In another embodiment, after determining the main contribution subareas corresponding to each main contribution industry type, the method further includes:

when the wind speed is greater than a wind speed threshold value, determining a relevant subarea combination with a wind direction upstream-downstream relation from all main contribution subareas, wherein the relevant subarea combination comprises an upstream subarea and a downstream subarea;

determining diffusion duration of the first precursor discharged from the upstream sub-area to the downstream sub-area, if the diffusion duration exceeds a preset time, increasing a first contribution quantification value of a main contribution industry type corresponding to the upstream sub-area, and keeping the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-area unchanged;

if the diffusion duration does not exceed the preset time, when the first sensitivity coefficient of the downstream sub-area is larger than the first sensitivity coefficient of the upstream sub-area, reducing the first contribution quantization value of the main contribution industry type corresponding to the downstream sub-area, and increasing the first contribution quantization value of the main contribution industry type corresponding to the upstream sub-area, wherein the reduction amount corresponding to the downstream sub-area is larger than the increase amount corresponding to the upstream sub-area;

And when the first sensitivity coefficient of the downstream sub-area is smaller than that of the upstream sub-area, reducing the first contribution quantization value of the main contribution industry type corresponding to the downstream sub-area, and increasing the first contribution quantization value of the main contribution industry type corresponding to the upstream sub-area, wherein the reduction amount corresponding to the downstream sub-area is smaller than that corresponding to the upstream sub-area.

Specifically, an associated sub-region combination with a wind direction upstream-downstream relation is determined from each main contribution sub-region, when the wind speed is greater than the wind speed threshold, if the diffusion duration of the first precursor of the upstream sub-region to the downstream sub-region in the associated sub-region combination does not exceed the preset time, it is indicated that the first precursor of the upstream sub-region may diffuse into the downstream sub-region to generate ozone, resulting in a higher first contribution quantization value of the main contribution industry type corresponding to the downstream sub-region, so that the first contribution quantization value of the main contribution industry type corresponding to the downstream sub-region is reduced, and the first contribution quantization value of the main contribution industry type corresponding to the upstream sub-region is increased, and the reduction amount of the first contribution quantization value is obtained from the reduction amount matching table according to the first sensitivity coefficient of the downstream sub-region, and the reduction amount is larger, the increase amount of the first contribution quantization value is obtained from the increase amount matching table according to the first sensitivity coefficient of the upstream sub-region, and the increase amount is larger. The first sensitivity coefficient of the upstream sub-region is greater than the first sensitivity coefficient of the downstream sub-region by an amount greater than the amount of decrease; the first sensitivity coefficient of the upstream sub-zone is less than the first sensitivity coefficient of the downstream sub-zone, and the increase is less than the decrease.

If the diffusion duration exceeds the preset time, the first precursor in the upstream sub-area is not diffused into the downstream sub-area to generate ozone before the preset time is over, the first contribution quantification value of the main contribution industry type corresponding to the upstream sub-area is increased, and the first contribution quantification value of the main contribution industry type corresponding to the downstream sub-area is kept unchanged.

Referring to fig. 3, an embodiment of the present application discloses a schematic flow chart of another method for evaluating ozone pollution sources, which can be implemented by a computer program or can be run on an evaluation device for evaluating ozone pollution sources based on von neumann system. The computer program can be integrated in an application or can be run as a stand-alone tool class application, and specifically comprises:

s201: acquiring a first sensitivity coefficient of a first precursor to ozone generation and a second sensitivity coefficient of a second precursor to ozone generation of each subarea in a preset area within a preset time, determining the sum of the coefficients of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and selecting a target subarea with the sum of the coefficients exceeding a first preset value from each subarea.

Specifically, reference may be made to step S101, which is not described herein.

S202: screening the first subareas in which the difference between the first sensitivity coefficient and the second sensitivity coefficient exceeds a second preset value from the target subareas, and determining a first common industry type in which the discharge rate of the first precursor in each first subarea exceeds a first preset discharge rate in a preset time.

Specifically, after the first subareas are screened out from the target subareas, the emission rate of the first precursor of each factory enterprise in each first subarea is determined within a preset time, then the name of the factory enterprise is input based on the enterprise information query platform, the industry type of the factory enterprise is determined, the emission rates of the first precursor of the factory enterprise of the same industry type are added, and finally the emission rate of the first precursor corresponding to each industry type in each first subarea is determined.

For each first sub-zone, comparing the emission rate of the first precursor of each industry type therein with a first preset emission rate, and if the emission rate of the first precursor exceeds the first preset emission rate, indicating that the first precursor emitted by the industry type within the first sub-zone has a greater impact on ozone generation, determining the industry type as the target industry type. Finally, counting the number of the first subareas with the same target industry type, if the number exceeds a first preset number, indicating that the target industry type exists in more first subareas in the preset area, wherein the ozone generation of the whole preset area is greatly influenced by the first precursor discharged by the target industry type, and if the target industry type is an industry type easy to generate ozone, determining the target industry type as a first common industry type.

S203: and determining a first contribution quantification value of each first common industry type to ozone pollution in a preset time according to the first enterprises corresponding to each first common industry type in the preset area.

Specifically, reference may be made to step S103, which is not described herein.

S204: the number of sub-areas to which each first common industry type belongs is counted.

S205: and multiplying the second correction coefficient by the first contribution quantized value of the corresponding first common industry type to obtain the final first contribution quantized value of the corresponding first common industry type.

Specifically, after determining the first contribution quantization value corresponding to each first common industry type, determining the second correction coefficient corresponding to the first contribution quantization value of the first common industry type according to the number of sub-areas where each first common industry type exists, where one possible determination method is as follows: the number of second correction coefficients corresponding to the number of matches is greater than 1 from the coefficient matching table, and the greater the number, the greater the influence of the emission of the first precursor on ozone generation and the corresponding second correction coefficient is, which indicates that the first common industry type is in the preset area. And finally multiplying the second correction coefficient by the first contribution quantized value of the corresponding first common industry type to obtain a final first contribution quantized value, thereby realizing calibration of the main industry source evaluation of ozone pollution.

S206: screening the second sub-areas of each target sub-area for a difference between the second sensitivity coefficient and the first sensitivity coefficient exceeding a second preset value, and determining a second common industry type in which the discharge rate of the second precursor in each second sub-area exceeds a second preset discharge rate within a preset time.

S207: and determining a second contribution quantification value of each second common industry type to ozone pollution in a preset time according to a second enterprise corresponding to each second common industry type in the preset area.

Specifically, refer to steps S104-S105, which are not described herein.

The implementation principle of the method for evaluating the ozone pollution source in the embodiment of the application is as follows: screening the discharged first precursor from the target sub-area that is more prone to ozone generation is easier to generate the first sub-area of ozone, and then determining contributions of each first common industry type to ozone pollution based on determining the first common industry type from the first sub-area that is more representative of ozone generation effects. Similarly, a second sub-region in which the second precursor is more prone to ozone generation is screened from the target sub-region, and then the contribution of each second common industry type to ozone pollution is determined based on the determination of a first common industry type from the second sub-region that is more representative of ozone generation effects. Thereby determining the distribution of ozone pollution sources in a specific industry.

The following are device embodiments of the present application, which may be used to perform method embodiments of the present application. For details not disclosed in the device embodiments of the present application, please refer to the method embodiments of the present application.

Fig. 4 is a schematic structural diagram of an apparatus for evaluating ozone pollution sources according to an embodiment of the present application. The evaluation device applied to the ozone pollution source can be implemented as all or part of the device by software, hardware or a combination of both. The apparatus 1 comprises a target area determination module 11, a first industry determination module 12, a first pollution evaluation module 13, a second industry determination module 14 and a second pollution evaluation module 15.

The target area determining module 11 is configured to obtain a first sensitivity coefficient of a first precursor to generate ozone and a second sensitivity coefficient of a second precursor to generate ozone in each subarea within a preset time, determine a sum of coefficients of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and select a target subarea with a coefficient sum exceeding a first preset value from each subarea, where the first precursor is a volatile organic compound, the second precursor is a nitrogen oxide, and the larger the sensitivity coefficient is, the larger the influence on ozone generation is;

A first industry determination module 12, configured to screen first sub-areas in each target sub-area, where a difference between coefficients obtained by subtracting the second sensitivity coefficient from the first sensitivity coefficient exceeds a second preset value, and determine a first common industry type in which a discharge rate of the first precursor in each first sub-area exceeds a first preset discharge rate within a preset time;

the first pollution evaluation module 13 is configured to determine, according to a first enterprise corresponding to each first common industrial type in the preset area, a first contribution quantization value of each first common industrial type to ozone pollution in a preset time;

a second industry determination module 14, configured to screen a second sub-area in which a difference between coefficients obtained by subtracting the first sensitivity coefficient from the second sensitivity coefficient exceeds a second preset value from each target sub-area, and determine a second common industry type in which a discharge rate of a second precursor in each second sub-area exceeds a second preset discharge rate within a preset time;

the second pollution evaluation module 15 is configured to determine, according to a second enterprise corresponding to each second common industrial type in the preset area, a second contribution quantization value of each second common industrial type to ozone pollution in the preset time.

Optionally, as shown in fig. 5, the apparatus 1 further comprises an emission threshold adjustment module 16, in particular for:

Acquiring a first contribution value of a first precursor to ozone concentration change and a second contribution value of a second precursor to ozone concentration change of a preset region within a preset time;

Optionally, the first pollution evaluation module 13 is specifically configured to:

counting the initial discharge rate of the first precursor of each first enterprise corresponding to each first common industry type in a preset area;

Optionally, the first pollution evaluation module 13 is further specifically configured to:

multiplying the first sensitivity coefficient of the subarea corresponding to each first enterprise by the initial discharge rate of the corresponding first precursor to obtain an initial discharge influence value;

the initial emission influence values of the first enterprises in the specific subareas are multiplied by the corresponding first correction coefficients to obtain corrected influence values, and the corrected influence values of the first enterprises belonging to the same first common industry type in the specific subareas and the initial emission influence values of the other first enterprises belonging to the same first common industry type outside the specific subareas are summed to obtain emission influence values of each first common industry type.

Optionally, the apparatus 1 further comprises a contribution quantization correction module 17, in particular for:

Counting the number of subareas to which each first common industry type belongs;

and multiplying the second correction coefficients by the first contribution quantized values of the corresponding first common industry types from the second correction coefficients corresponding to the matching number in the preset coefficient matching table to obtain final first contribution quantized values of the corresponding first common industry types, wherein the second correction coefficients are larger than 1, and the larger the number is, the larger the corresponding second correction coefficients are.

Optionally, the apparatus 1 further includes a key region determining module 18, specifically configured to:

Optionally, the apparatus 1 further comprises an industry quantization adjustment module 19, specifically configured to:

It should be noted that, when the evaluation device for ozone pollution sources provided in the above embodiment performs the evaluation method for ozone pollution sources, only the division of the above functional modules is used as an example, in practical application, the above functional distribution may be performed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules, so as to complete all or part of the functions described above. In addition, the device for evaluating the source of ozone pollution provided in the above embodiment and the method for evaluating the source of ozone pollution are the same conception, and detailed implementation processes are shown in the method embodiment, and are not described herein.

The embodiment of the application also discloses a computer readable storage medium, and the computer readable storage medium stores a computer program, wherein the computer program is executed by a processor, and the method for evaluating the ozone pollution source of the embodiment is adopted.

The computer program may be stored in a computer readable medium, where the computer program includes computer program code, where the computer program code may be in a source code form, an object code form, an executable file form, or some middleware form, etc., and the computer readable medium includes any entity or device capable of carrying the computer program code, a recording medium, a usb disk, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM), a Random Access Memory (RAM), an electrical carrier signal, a telecommunication signal, a software distribution medium, etc., where the computer readable medium includes, but is not limited to, the above components.

The method for evaluating the ozone pollution source of the embodiment is stored in the computer readable storage medium through the computer readable storage medium, and is loaded and executed on a processor so as to facilitate the storage and application of the method.

The embodiment of the application also discloses electronic equipment, wherein a computer program is stored in a computer readable storage medium, and when the computer program is loaded and executed by a processor, the method for evaluating the ozone pollution source is adopted.

The electronic device may be an electronic device such as a desktop computer, a notebook computer, or a cloud server, and the electronic device includes, but is not limited to, a processor and a memory, for example, the electronic device may further include an input/output device, a network access device, a bus, and the like.

The processor may be a Central Processing Unit (CPU), or of course, according to actual use, other general purpose processors, digital Signal Processors (DSP), application Specific Integrated Circuits (ASIC), ready-made programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc., and the general purpose processor may be a microprocessor or any conventional processor, etc., which is not limited in this application.

The memory may be an internal storage unit of the electronic device, for example, a hard disk or a memory of the electronic device, or may be an external storage device of the electronic device, for example, a plug-in hard disk, a Smart Memory Card (SMC), a secure digital card (SD), or a flash memory card (FC) provided on the electronic device, or the like, and may be a combination of the internal storage unit of the electronic device and the external storage device, where the memory is used to store a computer program and other programs and data required by the electronic device, and the memory may be used to temporarily store data that has been output or is to be output, which is not limited in this application.

The method for evaluating the ozone pollution source in the embodiment is stored in the memory of the electronic device and is loaded and executed on the processor of the electronic device, so that the method is convenient to use.

The foregoing is merely exemplary embodiments of the present disclosure and is not intended to limit the scope of the present disclosure. That is, equivalent changes and modifications are contemplated by the teachings of this disclosure, which fall within the scope of the present disclosure. Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure herein. This application is intended to cover any adaptations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a scope and spirit of the disclosure being indicated by the claims.

Claims

1. A method of evaluating a source of ozone pollution, the method comprising:

if the first contribution value is smaller than the second contribution value, determining a preset initial emission rate threshold as a first preset emission rate, and reducing the initial emission rate threshold to obtain a second preset emission rate;

determining a first contribution quantization value of each first common industry type to ozone pollution in the preset time according to a first enterprise corresponding to each first common industry type in the preset area, wherein the determining the first contribution quantization value of each first common industry type to ozone pollution in the preset time according to the first enterprise corresponding to each first common industry type in the preset area specifically comprises the following steps:

determining a first contribution quantification value of the corresponding first common industry type to the first contribution value within the preset time according to the ratio of the emission influence value of each first common industry type to the sum of all emission influence values;

If the first contribution quantity does not exceed a first quantity threshold value, adjusting the preset time to be the maximum diffusion duration so as to redetermine whether the preset area is an ozone key monitoring area;

2. The method according to claim 1, wherein the obtaining the emission impact value of each of the first common industry types based on the first sensitivity coefficient of the corresponding sub-region of each of the first enterprises and the corresponding initial emission rate of the first precursor specifically comprises:

3. The method of claim 1, wherein after determining the first contribution quantification value of each of the first common industry types to ozone pollution for the predetermined time, further comprising:

4. The method according to claim 1, wherein after determining, from the sub-areas, that the most number of main contribution sub-areas exist for each of the main contribution industry types, the method further comprises:

5. An apparatus for evaluating a source of ozone pollution, comprising:

a target area determining module (11) for obtaining a first sensitivity coefficient of a first precursor to ozone generation and a second sensitivity coefficient of a second precursor to ozone generation in a preset time of each subarea in a preset area, determining a coefficient sum of the first sensitivity coefficient and the second sensitivity coefficient corresponding to each subarea, and selecting a target subarea with the coefficient sum exceeding a first preset value from each subarea, wherein the first precursor is a volatile organic matter, the second precursor is a nitrogen oxide, and the larger the sensitivity coefficient is, the larger the influence on ozone generation is;

An emission threshold adjustment module (16) for obtaining a first contribution value of a first precursor to the ozone concentration variation and a second contribution value of a second precursor to the ozone concentration variation of the preset region within the preset time;

a first industry determination module (12) for screening a first sub-region of each target sub-region for which a difference between coefficients of the first sensitivity coefficient minus the second sensitivity coefficient exceeds a second preset value, and determining a first common industry type in which a discharge rate of the first precursor in each of the first sub-regions exceeds a first preset discharge rate within the preset time;

a first pollution evaluation module (13) configured to determine, according to a first enterprise corresponding to each first common industry type in the preset area, a first contribution quantization value of each first common industry type to ozone pollution in the preset time, where the determining, according to a first enterprise corresponding to each first common industry type in the preset area, a first contribution quantization value of each first common industry type to ozone pollution in the preset time specifically includes:

the key region determining module (18) is used for selecting a first common industry type with a first contribution quantization value exceeding a preset quantization value to determine the first common industry type as a main contribution industry type, and determining main contribution subregions with the largest enterprise existence number corresponding to each main contribution industry type from the subregions;

a second industry determination module (14) for screening a second sub-region in which a difference between coefficients of the second sensitivity coefficient subtracted from the first sensitivity coefficient exceeds a second preset value from each of the target sub-regions, and determining a second common industry type in which a discharge rate of the second precursor in each of the second sub-regions exceeds a second preset discharge rate for the preset time;

and the second pollution evaluation module (15) is used for determining a second contribution quantification value of each second common industry type to ozone pollution in the preset time according to a second enterprise corresponding to each second common industry type in the preset area.

6. A computer readable storage medium having a computer program stored therein, characterized in that the method according to any of claims 1-4 is employed when the computer program is loaded and executed by a processor.

7. An electronic device comprising a memory, a processor and a computer program stored in the memory and capable of running on the processor, characterized in that the method according to any of claims 1-4 is used when the computer program is loaded and executed by the processor.