Gas pollution leakage tracing method, device and equipment under complex diffusion conditions

Abstract

The present application relates to the field of gas leakage monitoring, and embodiments thereof provide a gas pollution leakage tracing method and device under complex diffusion conditions and electronic equipment. The gas pollution leakage tracing method under complex diffusion conditions comprises: establishing a fingerprint spectrum of surrounding pollution sources, and setting typical factors of each pollution source; the typical factors are used to represent typical pollution components of the pollution sources; obtaining receptor pollution data of a pollution area; calculating an analytical error between a tracing result based on a chemical mass balance model and the receptor pollution data, determining a pollution component as an error factor according to a contribution degree of the analytical error; and determining a variable weight chemical mass balance model according to the typical factors and the error factor, the variable weight chemical mass balance model being used to determine a pollution source. The embodiments of the present application realize accurate tracing of pollutants under complex diffusion conditions which are generally present in a field environment.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of gas leakage monitoring, in particular to a gas pollution leakage tracing method under complex diffusion conditions, a gas pollution leakage tracing device under complex diffusion conditions, an electronic device and a computer readable storage medium. BACKGROUND

[0002] With the rapid development of industrial construction and urbanization, resources and energy are rapidly consumed, and the atmospheric environment will face unprecedented pressure. Understanding the emission characteristics of atmospheric pollutants in key areas, qualitatively and quantitatively elucidating the sources and contributions of characteristic pollutants, is an important condition for fully achieving the "Ten Air Pollution" goals and scientifically and reasonably establishing a long-term mechanism for air pollution prevention and control.

[0003] At present, the chemical mass balance model analysis method is one of the important source analysis techniques. Under the premise of not considering the detection error of the analysis instrument, the tracing error mainly comes from the following three aspects. First, unknown source interference; unknown source interference will cause the receptor data to be unmatched with the established pollution source fingerprint spectrum, thereby causing analysis error. For this problem, the patent pollutant source analysis method provides a solution; second, the emission characteristics of pollution sources change. Not matching with the pollution source fingerprint spectrum established in the early stage; third, the complex conditions on site affect the pollution source diffusion transmission path. The pollution source fingerprint spectrum constructed in the early stage has obvious differences with the pollution characteristics actually diffused to the receptor area. For the above-mentioned second and third error sources, there is currently a lack of effective solutions.

[0004] The chemical mass balance model (CMB) is one of several receptor models applied to air pollution source management. The model was first proposed in 1972 and officially named chemical mass balance method in 1980. The CMB model established by this method is the most studied and most widely used receptor model in the actual work of air particulate matter source analysis. SUMMARY

[0005] The purpose of the embodiments of the present application is to provide a gas pollution leakage tracing method and device under complex diffusion conditions and an electronic device to solve some problems in the prior art.

[0006] In order to achieve the above object, the first aspect of the present application provides a gas pollution leakage source tracing method under complex diffusion conditions, comprising: establishing a fingerprint of a surrounding pollution source, and setting a typical factor of each pollution source; the typical factor is used to represent a typical pollution component of the pollution source; obtaining receptor pollution data of a pollution area; calculating an analytical error between a source tracing result based on a chemical mass balance model and the receptor pollution data, determining a pollution component as an error factor according to a contribution degree of the analytical error; determining a variable weight chemical mass balance model according to the typical factor and the error factor, and the variable weight chemical mass balance model is used to determine a pollution source.

[0007] Preferably, the fingerprint of the surrounding pollution source is established, comprising: taking a pollution emission source that has a pollution impact on a receptor area as the surrounding pollution source; constructing the fingerprint with main pollution components and volume concentration proportions of the pollution source diffused to the receptor area, or constructing the fingerprint with main pollution components and volume concentration proportions of the pollution source emitted.

[0008] Preferably, the typical factor of each pollution source is set, comprising: if a pollution component in the fingerprint of the pollution source meets the following conditions, the pollution component is the typical factor of the pollution source: the volume concentration percentage of the pollution component in the fingerprint of the pollution source exceeds a preset threshold value; exceeds a preset multiple of an average volume concentration proportion of the pollution component; and a pollution volume concentration percentage contribution rate of the pollution source to the receptor point under a specific wind direction exceeds a preset contribution rate threshold value.

[0009] Preferably, the receptor pollution data of the pollution area is obtained, comprising: obtaining an environmental air sample of the receptor area; and analyzing the environmental air sample by using a pollution component qualitative and quantitative analysis device to obtain pollution components and volume concentration proportions.

[0010] Preferably, the analytical error between the source tracing result based on the chemical mass balance model and the receptor pollution data is calculated, comprising: obtaining the analytical error SSE by the following steps:

[0011]

[0012] Wherein, N is the number of surrounding pollution sources, M is the number of pollution components; p ij is a contribution volume concentration of the i th pollution source in the analytical result to the j th pollution component of the receptor area; c j is a volume concentration of the j th pollution component in the receptor pollution data.

[0013] Preferably, the p ij is obtained by the following steps: p ij = F ij · S i ; wherein, F ija volume concentration ratio of the jth pollution component in the fingerprint of the ith pollution source; S i a total volume concentration of the ith pollution source to the pollution component of the receptor region, and the result of the chemical mass balance model-based source tracing.

[0014] Preferably, determining the pollution component as an error factor according to the contribution degree of the analytical error comprises: calculating an analytical sub-error of each pollution component when the analytical error is greater than K times of the baseline error; and setting the pollution component as the error factor if the analytical sub-error of any pollution component exceeds a preset ratio threshold in the analytical error and exceeds a preset multiple of the average of the analytical sub-errors of the pollution components.

[0015] Preferably, calculating the analytical sub-error of each pollution component comprises:

[0016] wherein, SSE j an analytical sub-error corresponding to the jth pollution component; N is the number of the surrounding pollution sources; p ij a contribution volume concentration of the ith pollution source to the jth pollution component in the receptor pollution data in the analytical result; c j a volume concentration of the jth component in the receptor pollution data.

[0017] Preferably, determining the variable-weight chemical mass balance model according to the typical factor and the error factor comprises obtaining the variable-weight chemical mass balance model according to the following relationship:

[0018]

[0019] wherein, k ij is the jth pollution component of the ith pollution source, if the jth pollution component is the error factor and the ith pollution source is the typical factor, then k ij is the to-be-optimized weight coefficient, otherwise k ij is a fixed value 1; E j is the analytical concentration error of the jth pollution component.

[0020] Preferably, the method further comprises: correcting the fingerprint of the surrounding pollution source by using the error factor.

[0021] Preferably, determining the pollution source by using the variable-weight chemical mass balance model comprises: taking the minimization of the analytical error SSE of the variable-weight chemical mass balance model as the principle, calculating the to-be-optimized weight coefficient k ij , and further analyzing to obtain a total volume concentration S i of the ith pollution source to the pollution component of the receptor region, and determining the pollution source.

[0022] Preferably, based on the variable weight chemical mass balance model calculation result, the contribution concentration P' of the i th pollution source to the j th pollution component of the receptor region is obtained ij , and the steps are as follows:

[0023] P' ij =k ij ·F ij ·S i ;

[0024] Wherein, F ij is the volume concentration proportion of the j th pollution component in the fingerprint spectrum of the i th pollution source; k ij is the to-be-optimized weight coefficient of the j th pollution component of the i th pollution source; S i is the total volume concentration of the i th pollution source to the pollution component of the receptor region; and the above k ij , S i are the tracing results based on the variable weight chemical mass balance model.

[0025] In the present application, a device for tracing gas pollution leakage under complex diffusion conditions is also provided, which comprises: a typical factor determination module, which establishes a fingerprint spectrum of surrounding pollution sources and sets typical factors of each pollution source; the typical factors are used to represent typical pollution components of the pollution sources; a data acquisition module, which is used to acquire receptor pollution data of a pollution region; an error factor determination module, which is used to calculate analytical errors between tracing results based on a chemical mass balance model and the receptor pollution data, and determine pollution components as error factors according to contribution degrees of the analytical errors; and a model construction module, which is used to determine a variable weight chemical mass balance model according to the typical factors and the error factors, and the variable weight chemical mass balance model is used to determine pollution sources.

[0026] In the present application, an electronic device is also provided, which comprises: at least one processor; and a memory connected with the at least one processor in communication; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to make the at least one processor execute the foregoing gas pollution leakage tracing method under complex diffusion conditions.

[0027] In the present application, a non-transient computer readable storage medium is also provided, which stores computer executable instructions, and the computer executable instructions are used to make the computer execute the foregoing gas pollution leakage tracing method under complex diffusion conditions.

[0028] In a fourth aspect of the present application, a computer readable storage medium is also provided, in which instructions are stored, which when executed on a computer, cause the computer to perform the steps of the aforementioned method for gas pollution leakage source tracing under complex diffusion conditions.

[0029] A computer program product is also provided in the present application, comprising a computer program which, when executed by a processor, implements the aforementioned method for gas pollution leakage source tracing under complex diffusion conditions.

[0030] The aforementioned technical solution has at least the following beneficial effects:

[0031] The embodiment of the present application establishes a variable weight chemical mass balance model based on typical factors and error factors, and provides an overall technical route for the modification of the chemical mass balance model under complex conditions. The problems of low accuracy and large analysis error of existing tracing analysis methods under the condition that the complex conditions affect the diffusion and transmission characteristics of the pollution source and the current pollution source emission characteristics do not match the established fingerprint are solved, and the accuracy and identification efficiency of gas pollution leakage source tracing under complex diffusion conditions commonly existing in the field environment are improved.

[0032] Other features and advantages of the embodiments of the present application will be described in detail in the following specific implementation part. BRIEF DESCRIPTION OF DRAWINGS

[0033] The accompanying drawings are included to provide a further understanding of the embodiments of the present application, and constitute a part of the specification, and are used together with the following specific implementation to explain the embodiments of the present application, but do not constitute a limitation on the embodiments of the present application. In the drawings:

[0034] Figure 1 The steps of the method for gas pollution leakage source tracing under complex diffusion conditions according to the embodiments of the present application are schematically shown;

[0035] Figure 2 The implementation schematic diagram of the method for gas pollution leakage source tracing under complex diffusion conditions according to the embodiments of the present application is schematically shown;

[0036] Figure 3 The relationship between the surrounding pollution source and the receptor area according to the embodiments of the present application is schematically shown;

[0037] Figure 4 The environmental distribution according to the embodiments of the present application is schematically shown;

[0038] Figure 5 The structure schematic diagram of the device for gas pollution leakage source tracing under complex diffusion conditions according to the embodiments of the present application is schematically shown;

[0039] Figure 6A schematic diagram of the first half of the CMB and variable weight CMB backtracking analysis result according to an embodiment of the present application is shown schematically.

[0040] Figure 7 A schematic diagram of the second half of the CMB and variable weight CMB backtracking analysis result according to an embodiment of the present application is shown schematically. DETAILED DESCRIPTION

[0041] The specific embodiments of the embodiments of the present application are described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the embodiments of the present application, and are not used to limit the embodiments of the present application.

[0042] Embodiment 1

[0043] Figure 1 A schematic diagram of the steps of the gas pollution leakage backtracking method under complex diffusion conditions according to an embodiment of the present application is shown schematically, as shown in the figure, a gas pollution leakage backtracking method under complex diffusion conditions, the method comprises: Figure 1

[0044] S01, a fingerprint of a surrounding pollution source is established, and a typical factor of each pollution source is set; the typical factor is used to represent the typical pollution component of the pollution source; specifically, the surrounding pollution source in the embodiment can be a pollution emission source that can have a pollution impact on the receptor region; since the pollution gas generated by the pollution source is easily affected by the complex conditions in the field during diffusion transmission, some substances change the diffusion path, so that the pollution characteristics emitted by the pollution source are different from the actual pollution characteristics affecting the receptor region. In order to accurately analyze the pollution source, it is suggested to construct a fingerprint with the main pollution components and their volume concentration ratios of the pollution source diffused to the receptor region. However, considering the complexity of the field environment, the receptor region is generally simultaneously affected by multiple pollution sources, and it is difficult to determine the main pollution components and their ratios of each single emission source diffused to the receptor region. At this time, the main pollution components and their volume concentration ratios emitted by the pollution source can be used to construct a fingerprint, and the main pollution components and their volume concentration ratios can also be selected at a reasonable position in the diffusion transmission path to establish a fingerprint.

[0045] S02, obtaining receptor pollution data of the pollution area. The receptor pollution data here is the pollution analysis result obtained by analyzing the pollutants collected from the receptor region by an instrument.

[0046] ​S03, calculate the analytical error between the tracing result based on the chemical mass balance model and the receptor pollution data, determine the pollution components as error factors according to the contribution of the analytical error; according to the principle of the chemical mass balance model, if the pollution components and the proportion of each pollution source diffused to the receptor region are the same as the fingerprint spectrum, the qualitative and quantitative analysis equipment does not have detection error, and is not disturbed by other unknown sources, then the analytical result of the receptor pollution data based on the chemical mass balance will not have analytical error. Ignoring the detection error of the qualitative and quantitative analysis equipment and the disturbance of the unknown source, if the analytical error is small, it is considered that the field environment does not affect the diffusion and propagation of the pollution source, and the source analysis result can be directly output at this time. If the analytical error is large, further accurate source analysis is needed.

[0047] S04, determine a variable weight chemical mass balance model according to the typical factor and the error factor, the variable weight chemical mass balance model is used to determine the pollution source. This step sets the weight coefficient to be optimized in the original chemical mass balance model, and solves the weight coefficient to obtain the variable weight chemical mass balance model. And according to the variable weight chemical mass balance model, the pollutant tracing is performed.

[0048] Through the above embodiments, the problems of low accuracy and large analytical error of the existing tracing analysis method under the influence of the diffusion and transmission characteristics of the pollution source under complex conditions and the mismatch between the current pollution source emission characteristics and the established fingerprint spectrum are solved, and the accurate tracing of the pollutants under the complex diffusion conditions commonly existing in the field environment is realized.

[0049] Figure 2 An implementation schematic diagram of a gas pollution leakage tracing method under complex diffusion conditions in an embodiment of the present application is schematically shown. As shown in Figure 2 The main steps include: establishing a surrounding pollution source fingerprint spectrum, and setting a typical factor of each pollution source; obtaining receptor pollution data of a receptor region; calculating the analytical error based on the chemical mass balance model, and determining the error factor; establishing a variable weight chemical mass balance model based on the typical factor and the error factor, and calculating the pollution source.

[0050] Figure 3 A relationship diagram of a surrounding pollution source and a receptor region in an embodiment of the present application is schematically shown. As shown in Figure 3As shown, in the optional embodiment, the fingerprint of the surrounding pollution source is established, including: taking the pollution emission source that has an impact on the receptor area as the surrounding pollution source; and constructing the fingerprint with the main pollution components and their proportions in the receptor area or with the main pollution components and their proportions emitted by the pollution source. The fingerprint of the pollution source in this embodiment is composed of the main pollution components and their proportions emitted by the pollution source, and the fingerprint is constructed with the main pollution components and their proportions in the receptor area if the conditions permit, otherwise the fingerprint is constructed with the main pollution components and their proportions emitted by the pollution source. The number of pollution components contained in the general fingerprint should not be less than 8.

[0051] In some optional embodiments, the typical factor of each pollution source is set, including: if the pollution component in the fingerprint of the pollution source meets the following conditions, the pollution component is the typical factor of the pollution source: the volume concentration percentage of the pollution component in the fingerprint of the pollution source exceeds a preset threshold; exceeds a preset multiple of the average volume concentration percentage of the pollution component; and the pollution source contributes more than a preset contribution rate threshold to the volume concentration percentage of the pollutants at the receptor point in a specific wind direction. The typical factor of the pollution source belongs to the pollution component recorded in the fingerprint of the pollution source, the volume concentration percentage of the pollution component in the fingerprint of the pollution source exceeds a preset threshold, for example, 15%. And exceeds 2 times of 1 / M (i.e. the average volume concentration percentage of the pollution component), M is the number of pollution components, and 2 is a preset multiple. And the pollution source contributes more than a preset contribution rate threshold to the volume concentration percentage of the pollutants at the receptor point in a specific wind direction, for example, 30%; the annual wind frequency of the specific wind direction is not less than 5%. Through the above setting method, the pollution source and its pollution factor that have a significant impact on the pollutant concentration characteristics of the receptor area are determined, and are set as the typical factor of the pollution source. The typical factor is mainly used to solve the problem of the complex conditions affecting the diffusion and transmission characteristics of the pollution source mentioned above. Since any pollution component may change during the diffusion and transmission process, the change of the pollution component with small proportion is not sufficient to affect the tracing calculation result, and therefore it is necessary to screen out the pollution source and its pollution factor that have a significant impact on the pollutant concentration characteristics of the receptor area.

[0052] In some optional embodiments, the method for obtaining receptor pollution data of a pollution area comprises: obtaining an ambient air sample of the receptor area; and analyzing the ambient air sample by using a pollution component qualitative and quantitative analysis device to obtain pollution components and volume concentration proportions. The receptor area is an area to be analyzed to determine the pollution source, and obtaining receptor pollution data is a necessary condition for carrying out source tracing analysis and calculation. The specific method is as follows: obtaining an ambient air sample of the receptor area, and using a pollution component qualitative and quantitative analysis device to determine pollution components and volume concentration proportions; in order to avoid differences in analysis capabilities of different types of devices, it is recommended to obtain receptor pollution data and establish a pollution source fingerprint by using the same type of qualitative and quantitative analysis device; when the receptor area is a point, multiple air samples at different times at the position need to be continuously collected; when the receptor area is a volume, a surface, or a line, multiple air samples at different positions at the same time or at different times in the area need to be continuously collected. Increasing the sample size of the receptor pollution data can improve the accuracy of the source tracing analysis result, and at the same time, accidental situations caused by a single sample can be avoided. Therefore, it is recommended to continuously collect multiple air samples at different times, and for volume, surface, line (such as a factory boundary), and other receptor areas, different spatial positions have different pollution characteristics, and multiple air samples need to be collected at different spatial positions, such as five-point sampling or uniform sampling.

[0053] In some optional embodiments, the method for calculating an analysis error between a source tracing result based on a chemical mass balance model and the receptor pollution data comprises: obtaining the analysis error SSE by the following steps:

[0054] wherein N is the number of surrounding pollution sources, M is the number of pollution components; p ij is the contribution volume concentration of the i th pollution source to the j th pollution component of the receptor area in the analysis result; c j is the volume concentration of the j th pollution component in the receptor pollution data.

[0055] According to the chemical mass balance model, the p ij is obtained by the following steps: p ij = F ij · S i ; wherein F ij is the volume concentration proportion of the j th pollution component in the fingerprint of the i th pollution source; and S i is the total volume concentration of the i th pollution source to the pollution component of the receptor area, which is the source tracing result based on the chemical mass balance model.

[0056] In some optional embodiments, the pollution components are determined as error factors according to the contribution degrees of the analytical errors, comprising: when the analytical error is greater than K times of the reference error, calculating an analytical sub-error of each pollution component; if the analytical sub-error of any pollution component accounts for more than a preset proportion threshold in the analytical error, and more than a preset multiple of the average of the analytical sub-errors of the pollution components, the pollution component is set as an error factor. After the analytical error is calculated, the following method is given to determine whether the analytical error is small enough to represent that the calculation result is accurate enough: a reference error is set, if the SSE is less than K times of the reference error BSSE, it means that the analytical result is good and can be directly output; otherwise, the SSE is greater than K times of the reference error BSSE, that is, the analytical result is poor, and the error factor needs to be further determined, a variable weight CMB model is established, and the analytical result is calculated; the K is recommended to be 5%-20%; the calculation formula of the reference error is as follows:

[0057] In some optional embodiments, the method further comprises: calculating an analytical sub-error of each pollution component, comprising:

[0058]

[0059] Wherein, SSE j is the analytical sub-error of the jth pollution component; N is the number of the surrounding pollution sources; p ij is the contribution volume concentration of the ith pollution source in the analytical result to the jth pollution component in the receptor pollution data; c j is the volume concentration of the jth component in the receptor pollution data.

[0060] In some optional embodiments, the variable weight chemical mass balance model is determined according to the typical factors and the error factors, comprising: the variable weight chemical mass balance model is obtained according to the following relationship:

[0061]

[0062] Wherein, k ij is the to-be-optimized weight coefficient of the jth pollution component of the ith pollution source, if the jth pollution component is an error factor and the ith pollution source is a typical factor, k ij is the to-be-optimized weight coefficient, otherwise k ij is a fixed value 1. E j is the analytical concentration error of the jth pollution component.

[0063] The pollution source is determined by using the variable weight chemical mass balance model, comprising: taking the minimization of the analytical error SSE of the variable weight chemical mass balance model as the principle, calculating the to-be-optimized weight coefficient k ijand further analysis to obtain the total volume concentration S of the i th pollution source to the pollution component of the receptor area i , and further determine the pollution source.

[0064] Based on the calculation results of the variable weight chemical mass balance model, the contribution concentration p of the i th pollution source to the j th pollution component of the receptor area is obtained i ' j , the method is as follows:

[0065] p i ' j =k ij ·F ij ·S i ;

[0066] Wherein, F ij is the volume concentration ratio of the j th pollution component in the fingerprint spectrum of the i th pollution source; k ij is the to-be-optimized weight coefficient of the j th pollution component of the i th pollution source; S i is the total volume concentration of the i th pollution source to the pollution component of the receptor area; the above k ij , S i are the results of the variable weight chemical mass balance model.

[0067] In some optional embodiments, the method further comprises: correcting the fingerprint spectrum of the surrounding pollution source by using the error factor. Specifically, if k ij is less than 1, it means that the j th pollution component actually discharged by the i th pollution source is less than the set ratio in the fingerprint spectrum, or the j th pollution component is seriously attenuated in the pollution gas diffusion process; if k ij is greater than 1, it means that the j th pollution component of the i th pollution source or other pollution sources is more than the set ratio in the fingerprint spectrum; when the above conclusion appears in multiple analysis results, it is recommended to optimize and construct the fingerprint spectrum according to the current latest emission characteristics, and it is recommended to adjust the volume concentration ratio of the j th pollution component (which is determined as the error factor multiple times) in the fingerprint spectrum of the i th pollution source to F ij ×k ij .

[0068] In some optional embodiments, the variable weight chemical mass balance model is used to determine the pollution source, comprising: taking the minimization of the analysis error of the variable weight chemical mass balance model as the principle, calculating the to-be-optimized weight coefficient k ij , and further analyzing to obtain S i (the total volume concentration of the i th pollution source to the pollution component of the receptor area), and determining the pollution source. In this embodiment, a simple variable weight chemical mass balance model is also provided as follows:

[0069]

[0070] wherein k j is the weight coefficient of the jth pollution component; if the jth pollution component is an error factor and is a typical factor of any pollution source, k j is a positive number less than 1, and is generally set to 1 / 6-1 / 3, otherwise k j =1; on this basis, the analytical error SSE minimization principle is used for analytical calculation to obtain s i (the total volume concentration of the i th pollution source to the receptor area pollution component contribution). E j is the analytical concentration error of the jth pollution component.

[0071] Due to the simple variable weight chemical mass balance model, the analytical error weight of the error factor is reduced, which will to some extent cause the deviation of the calculated S i from the actual total contribution concentration of the i th pollution source to the receptor area pollution component, so it is suggested that the contribution concentration of the i th pollution source to the jth pollution component of the receptor area is calculated as follows:

[0072]

[0073] Further, the total contribution concentration of the i th pollution source to the receptor area pollution component is:

[0074] Example 2

[0075] This embodiment mainly provides a determination method of a pollution source typical factor. The pollution source typical factor belongs to a pollution component recorded in a pollution source fingerprint of the pollution source, the volume percentage of the pollution component in the pollution source fingerprint is more than 10%, and more than 2 times of 1 / M, M is the number of pollution components, and the pollution volume percentage contribution rate of the pollution source to the receptor point under a specific wind direction is more than 30%; the annual wind frequency of the specific wind direction is not less than 5%.

[0076] Figure 4 An environmental distribution schematic diagram according to an embodiment of the present application is schematically shown. As Figure 4 shown, a certain refining enterprise includes four pollution emission sources of residual oil, wax oil, reforming and tank area, and a receptor area is arranged in the south reserved land, and the pollution source of the receptor area is analyzed.

[0077] The fingerprint of the above four pollution sources is established, combined with the pollution emission characteristics of each device, and a total of 17 pollution components are included, and the following table shows the example fingerprint of the two pollution sources of residual oil and tank area.

[0078]

[0079]

[0080] Firstly, the typical factor of the pollution source is required to belong to the pollution components recorded in the pollution source fingerprint, the volume percentage of the pollution component in the pollution source fingerprint is more than 10%, and more than 2 times of 1 / 17, that is, the pollution components meeting the conditions of the residue hydrogenation pollution source include isobutane, 1-butene, n-butane, isopentane, and the pollution components meeting the conditions of the tank area pollution source include propylene, n-butane, isopentane, and octane.

[0081] Secondly, the related results of the previous source tracing analysis based on the chemical mass balance model are obtained, it is found that the residue hydrogenation device will have a greater pollution impact on the receptor area under winter conditions (especially under northwest wind conditions), and the maximum pollution contribution rate in the 32 times of source tracing results under this wind direction condition is 37%, and the average contribution rate is 12%, and the annual meteorological condition statistical results show that the northwest wind is the dominant wind direction in winter in this area, and the annual wind frequency is 11%; in summary, the pollution volume percentage contribution rate of the pollution source to the receptor point under the specific wind direction is more than 30%, and the annual wind frequency of the specific wind direction is not less than 5%.

[0082] Therefore, the typical factors of the residue hydrogenation pollution source are set as isobutane, 1-butene, n-butane, and isopentane.

[0083] Example 3

[0084] This embodiment mainly shows the calculation and identification method of the analysis error, the benchmark error, and the error factor based on the chemical mass balance model. A certain refining enterprise includes four pollution emission sources of residue, wax oil, reforming, and tank area, and a receptor area is arranged in the south reserved land. The fingerprint maps of the four pollution emission sources are established, and the typical factors of each emission source are identified, for example, the typical factors of the residue hydrogenation pollution source are isobutane, 1-butene, n-butane, and isopentane.

[0085] The receptor pollution data are obtained, the pollution sources of the receptor area are analyzed, the concentration values of each pollution component contributed by each pollution source are obtained, and the results are as follows.

[0086]

[0087] The analysis error is the error between the fitting results of the analyzed receptor pollution and the actual receptor pollution data, and the calculation formula is as follows:

[0088]

[0089] The analysis error of any jth pollution component is calculated by the following formula:

[0090]

[0091] The benchmark error is calculated as follows:

[0092] Where N is the number of pollution sources, which is 4, and M is the number of pollution components, which is 17.

[0093] Among them, p ij It is the volume concentration of the contribution of the i-th pollution source in the analysis results to the j-th pollutant in the receptor pollution data; for example, the volume concentration of the contribution of the first pollution source, residual oil hydrogenation, to the fifth pollutant in the receptor pollution data, isobutane, is 147 ppb.

[0094] Among them, c j is the volume concentration of the jth component in the receptor pollution data; for example, the volume concentration of the fifth pollutant component, isobutane, is 214 ppb. Calculate the analytical error SSE of each pollutant component j And the overall analytical error SSE, the results are as follows:

[0095]

[0096] After the calculation is completed, the identification process based on the error factor is as follows:

[0097] 1) If the SSE is less than K times the benchmark error BSSE, the analytical result is good and can be output directly; otherwise, if the SSE is greater than K times the benchmark error BSSE, the analytical result is poor and it is necessary to further clarify the error factor, establish a variable weight CMB model, and calculate the analytical result;

[0098] 2) If the analytical error SSE of any contaminant component j If the percentage of the analytical error (SSE) exceeds r and exceeds 2 times 1 / M, where M is the number of contaminants, then the contaminant is set as the error factor. Here, K is 5% and r is 20%.

[0099] According to the calculation results of SSE and BSSE, SSE = 58420, BSSE = 115131, that is, SSE is greater than 5% of the reference error BSSE, that is, the analytical result is poor, and the error factor needs to be further clarified; according to SSE j The calculation results show that the SSE of the sixth pollutant 1-butene j It is 56169, accounting for 96.1% of the analytical error SSE, which is more than 20%, that is, 1-butene is the error factor.

[0100] Example 4

[0101] This embodiment mainly provides a specific implementation method for establishing a variable weight chemical mass balance model based on typical factors and error factors.

[0102] The variable weight chemical mass balance model is known as follows:

[0103]

[0104] wherein k ij is the to-be-optimized weight coefficient of the jth pollution component of the ith pollution source, if the jth pollution component is an error factor and is a typical factor of the ith pollution source, then k ij is set to be the to-be-optimized weight coefficient, otherwise k ij is a fixed value of 1. E j is the analytical concentration error of the jth pollution component. The contribution concentration of the ith pollution source to the jth pollution component of the receptor region is: p ij = k ij · F ij · S i ; wherein k ij is the to-be-optimized weight coefficient of the jth pollution component of the ith pollution source; F ij is the volume concentration proportion of the jth pollution component in the fingerprint spectrum of the ith pollution source; s i is the total volume concentration of the ith pollution source to the pollution component of the receptor region, wherein k ij , S i are the calculation results of the variable weight chemical mass balance model.

[0105] For example, a certain refining enterprise includes four pollution emission sources of residual oil, wax oil, reforming, and tank area, and a receptor region is arranged in the south reserved land, and four pollution emission source fingerprint spectra are established; typical factors of each emission source are identified, such as the typical factors of the residual oil hydrogenation pollution source are isobutane, 1-butene, n-butane, and isopentane; error factors are identified, and since the SSE j of the 6th pollution component 1-butene in the analytical error SSE is 96.1%, which is an error factor.

[0106] Since the 6th pollution component 1-butene is an error factor, and 1-butene is a typical factor of the 1st (residual oil hydrogenation) and 2nd (tank area) pollution sources, then k 1,6 and k 2,6 are to-be-optimized weight coefficients, and the remaining k ij is a fixed value of 1.

[0107] At this time, the variable weight chemical mass balance model is established based on the typical factors and error factors as follows:

[0108]

[0109] The principle of minimizing the analytical error SSE is used to optimize and calculate k 1,6 and k 2,6Value and S i Generally set k ij ∈[1 / 3,3], then according to the variable weight chemical mass balance model, the contribution concentration of the i th pollution source to the j th pollution component of the receptor region is:

[0110] P'ij=k ij ·F ij ·S i .

[0111] Embodiment 5

[0112] Based on the same inventive concept, the application also provides a gas pollution leakage source tracing device under complex diffusion conditions. Figure 5 The structure of the gas pollution leakage source tracing device under complex diffusion conditions according to the embodiment of the application is schematically shown. As shown in the figure, Figure 5 The device comprises a typical factor determination module, which establishes a fingerprint of a surrounding pollution source and sets a typical factor of each pollution source; the typical factor is used to represent a typical pollution component of the pollution source; a data acquisition module, which is used to acquire receptor pollution data of a pollution region; an error factor determination module, which is used to calculate an analytical error between a source tracing result based on a chemical mass balance model and the receptor pollution data, and determine a pollution component as an error factor according to a contribution degree of the analytical error; and a model construction module, which is used to determine a variable weight chemical mass balance model according to the typical factor and the error factor, and the variable weight chemical mass balance model is used to determine a pollution source.

[0113] In some optional embodiments, the establishment of the fingerprint of the surrounding pollution source comprises: taking a pollution emission source that has a pollution impact on a receptor region as the surrounding pollution source; and constructing the fingerprint with main pollution components and their volume concentration proportions of the pollution source diffused to the receptor region, or constructing the fingerprint with main pollution components and their volume concentration proportions emitted by the pollution source.

[0114] In some optional embodiments, the setting of the typical factor of each pollution source comprises: if a pollution component in the fingerprint of a pollution source meets the following conditions, the pollution component is the typical factor of the pollution source: the volume concentration percentage of the pollution component in the fingerprint of the pollution source exceeds a preset threshold value; exceeds a preset multiple of the average volume concentration of the pollution component; and the pollution volume concentration percentage contribution rate of the pollution source to the receptor point under a specific wind direction exceeds a preset contribution rate threshold value.

[0115] In some optional embodiments, the acquisition of the receptor pollution data of the pollution region comprises: acquiring an environmental air sample of the receptor region; and analyzing the environmental air sample by using a pollution component qualitative and quantitative analysis device to obtain pollution components and volume concentration proportions.

[0116] In some optional embodiments, the analytical error between the tracing result based on the chemical mass balance model and the receptor pollution data is calculated, including: obtaining the analytical error SSE by the following steps:

[0117] Wherein, N is the number of surrounding pollution sources, M is the number of pollution components; p ij is the contribution volume concentration of the i th pollution source to the j th pollution component of the receptor area in the analytical result; c j is the volume concentration of the j th pollution component in the receptor pollution data.

[0118] In some optional embodiments, the p ij is obtained by the following steps: p ij = F ij · S i ;

[0119] Wherein, F ij is the volume concentration percentage of the j th pollution component in the i th pollution source according to the fingerprint; S i is the total volume concentration of the i th pollution source to the pollution component of the receptor area, and is the tracing result based on the chemical mass balance model.

[0120] In some optional embodiments, the pollution component is determined as an error factor according to the contribution degree to the analytical error, including: calculating the analytical sub-error of each pollution component when the analytical error is greater than K times of the reference error; if the analytical sub-error of any pollution component accounts for more than a preset proportion threshold in the analytical error, and more than a preset multiple of the volume average, the pollution component is set as an error factor.

[0121] In some optional embodiments, the device further comprises: calculating the analytical sub-error of each pollution component, including:

[0122]

[0123] Wherein, SSE j is the analytical sub-error corresponding to the j th pollution component; N is the number of surrounding pollution sources; p ij is the contribution volume concentration of the i th pollution source to the j th pollution component in the receptor pollution data in the analytical result; c j is the volume concentration of the j th component in the receptor pollution data.

[0124] In some optional embodiments, the variable weight chemical mass balance model is determined according to the typical factor and the error factor, including: obtaining the variable weight chemical mass balance model according to the following relationship:

[0125]

[0126] wherein k ij is a weight coefficient to be optimized for the jthpollutant component of the ithpollution source, and if the jthpollutant component is an error factor and is a typical factor of the ithpollution source, then k ij is a weight coefficient to be optimized, otherwise k ij is a fixed value 1; E j is the analytical concentration error of the jthpollutant component.

[0127] In some optional embodiments, the device further comprises: correcting the fingerprint of the peripheral pollution source by using the error factor.

[0128] In some optional embodiments, determining the pollution source by using the variable weight chemical mass balance model comprises: determining the pollution source based on the principle of minimizing the analytical error SSE of the variable weight chemical mass balance model.

[0129] In some optional embodiments, based on the calculation result of the variable weight chemical mass balance model, the contribution concentration P i ′ j of the ithpollution source to the jthpollutant component of the receptor region is obtained, and the method is as follows:

[0130] P' ij = k ij · F ij · S i ;

[0131] wherein F ij is the volume concentration proportion of the jthpollutant component in the fingerprint of the ithpollution source; k ij is a weight coefficient to be optimized for the jthpollutant component of the ithpollution source; S i is the total volume concentration of the ithpollution source to the pollutant component of the receptor region; and the above k ij , S i are the tracing results based on the variable weight chemical mass balance model.

[0132] The specific definitions of each functional module in the above device for tracing gas pollution leakage under complex diffusion conditions can refer to the definitions of the method for tracing gas pollution leakage under complex diffusion conditions in the above, which will not be repeated here. Each module in the above device can be realized by software, hardware, and a combination thereof, in whole or in part. Each module can be embedded in or independent of the processor in the computer device in hardware form, or can be stored in the memory in the computer device in software form, so as to be called and executed by the processor to perform the operations corresponding to each module.

[0133] Embodiment 6

[0134] To verify the compatibility of CMB model and variable weight chemical mass balance model to the random changes of receptor data, mainly involving the interference error of detection process and the abnormal error in diffusion process, test experiments are specially carried out, and the experimental scheme and process are as follows.

[0135] Based on the fingerprint spectra of 5 types of emission sources of a certain refining and chemical enterprise, i.e. reforming device (source B), wax oil and residual oil device (source C), gas separation device (source D), catalytic device (source E), atmospheric and vacuum distillation device (source F) and external enterprise emission source (source A), the receptor data based on the above emission source fingerprint spectra are simulated and synthesized, and the receptor data are calculated and analyzed by using CMB model and variable weight CMB respectively.

[0136] The specific process is as follows:

[0137] (1) Initial receptor data is generated. The initial receptor data is synthesized by superimposing 6 types of emission source fingerprint spectra according to a certain proportion, that is, the receptor region in the theoretical state is only affected by 6 types of emission sources, without the influence of other emission sources, without the influence of material changes in diffusion process, and the emission characteristics of the emission sources are completely consistent with the fingerprint spectra. In the synthesis process, the contribution rate of external enterprise is set to 70%, and the contribution rates of the other 5 types of emission sources are randomly allocated, and 200 groups of initial receptor data are obtained, that is, the actual value of source A contribution rate of the 200 groups of receptor data is 70%.

[0138] (2) Receptor data is generated. Randomly add (0-20)% error interference to each pollutant of the initial receptor data, and randomly select 3 substances to significantly reduce their concentration values, and the reduction is randomly set to (60-100)%, so as to represent the phenomenon that some pollutants cannot normally diffuse to the tracing region due to complex environmental influences such as weather during the pollution diffusion process, and the influence of the change of emission rules of each emission source on the receptor region data.

[0139] (3) CMB model calculation and analysis. Based on 6 types of emission source fingerprint spectra, the receptor data are calculated and analyzed by using CMB model, and the contribution rate of external enterprise emission source (source A) is obtained.

[0140] (4) Variable weight CMB model calculation and analysis. Based on 6 types of emission source fingerprint spectra, the receptor data are calculated and analyzed by using variable weight CMB model, and the contribution rate of external enterprise emission source (source A) is obtained. Figure 6 The first half part of the schematic diagram of CMB and variable weight CMB tracing analysis results according to the embodiment of the present application is schematically shown; Figure 7 The second half part of the schematic diagram of CMB and variable weight CMB tracing analysis results according to the embodiment of the present application is schematically shown.

[0141] As can be seen from the figure, the known source A true contribution rate is 70%, the statistical CMB model has 23 times of experiment in 200 times of experiments in which the relative error of source A analysis is greater than 10%, that is, the source A analysis contribution rate is less than 63% or greater than 77%, and the number of experiments in which the statistical CMB model has a relative error greater than 10% is only 5 times. At the same time, the maximum error of the CMB model is 47%, and the maximum error of the variable weight CMB model is 59%. That is, the variable weight CMB model has good analysis effect on the interference change of the receptor data caused by the complex environment.

[0142] Embodiment 7

[0143] The embodiment provides a non-volatile computer storage medium, which stores computer executable instructions, and the computer executable instructions can execute the method in any method embodiment and achieve the same technical effect.

[0144] The embodiment also provides a computer program product, which includes a computer program stored on a non-transitory computer readable storage medium, and the computer program includes program instructions, and when the program instructions are executed by a computer, the computer executes the method in the above aspects and achieves the same technical effect.

[0145] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can adopt a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt a computer program product in the form of one or more computer usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer usable program codes.

[0146] The present application is described with reference to flowcharts and / or block diagrams according to the method, device (system), and computer program product of the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to produce a machine, so that the instructions executed by the computer or other programmable data processing devices produce a device that implements the functions specified in the flowcharts and / or block diagrams. Figure 1 The function of one flow or multiple flows and / or blocks Figure 1 The device that implements the function specified in one block or multiple blocks.

[0147] These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the Figure 1 function specified in the flow or flows and / or blocks Figure 1 of the block or blocks.

[0148] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions that are executed on the computer or other programmable apparatus provide steps for implementing the Figure 1 function specified in the flow or flows and / or blocks Figure 1 of the block or blocks.

[0149] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0150] The memory can include non-persistent memory and / or volatile memory, such as a random access memory (RAM) including a cache area for the temporary storage of data. The memory can also include non-volatile memory, such as read only memory (ROM), electrically programmable read only memory (EPROM), electrically erasable programmable read only memory (EEPROM), flash memory, or non-volatile random access memory (NVRAM), for the storage of information, such as data files or program

[0151] Computer readable media includes permanent and non-permanent, removable and non-removable media implemented in any method or technology for the storage of information. The information can be computer readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read only memory (ROM), electrically erasable programmable read only memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital versatile discs (DVDs) or other optical storage, magnetic cassettes, magnetic tapes, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible to a computing device. According to the definition herein, computer readable media does not include transitory media, such as modulated data signals and carrier waves.

[0152] It should also be noted that the terms "comprising", "comprises" or other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but can also include other elements not expressly listed or inherent to such process, method, article, or apparatus.

[0153] The above embodiments are only used to illustrate the present application, but not to limit it. Instead of the above, various modifications and changes can be made to the application by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the application shall fall into the scope of the claims of the application.

Claims

1. A method for tracing the source of gas pollution leakage under complex diffusion conditions, characterized in that: The method comprises: Establish a fingerprint map of surrounding pollution sources and set a typical factor for each pollution source; the typical factor is used to characterize the typical pollutant component of the pollution source; set the typical factor for each pollution source, including: if the pollutant component in the fingerprint map of the pollution source meets the following conditions, then the pollutant component is the typical factor of the pollution source: the volume concentration percentage of the pollutant component in the fingerprint map of the pollution source exceeds a preset threshold; exceeds a preset multiple of the average volume concentration of the pollutant component; and the contribution rate of the pollutant volume concentration percentage of the pollution source to the receptor point under a specific wind direction exceeds a preset contribution rate threshold; Obtain receptor contamination data for polluted areas; Calculating the analytical error between the traceability result based on the chemical mass balance model and the receptor pollution data, and determining the pollution component as an error factor according to the contribution to the analytical error; calculating the analytical error between the traceability result based on the chemical mass balance model and the receptor pollution data, including: obtaining the analytical error SSE by the following steps: Where N is the number of surrounding pollution sources and M is the number of pollution components; is the volume concentration of the jth pollution component contributed by the i-th pollution source to the receptor area in the analytical results; The volume concentration of the jth pollutant component in the receptor pollution data; determining the pollutant component as an error factor based on its contribution to the analytical error, including: when the analytical error is greater than K times the baseline error, calculating the analytical sub-error of each pollutant component; if the analytical sub-error of any pollutant component exceeds a preset contribution threshold and exceeds a preset multiple of the mean of the analytical sub-errors of all pollutants, then setting the pollutant component as the error factor; Determining a variable weight chemical mass balance model based on the typical factors and the error factors, wherein the variable weight chemical mass balance model is used to determine the pollution source; including obtaining the variable weight chemical mass balance model based on the following relationship: in, is the weight coefficient to be optimized for the jth pollution component of the i-th pollution source. If the jth pollution component is an error factor and is a typical factor of the i-th pollution source, then set is the weight coefficient to be optimized, otherwise is a fixed value of 1; is the analytical concentration error of the jth pollutant component; described Obtained through the following steps: ； in, is the volume concentration ratio of the jth pollution component in the fingerprint of the i-th pollution source; is the total volume concentration of the pollution components contributed by the i-th pollution source to the receptor area, which is the source tracing result based on the chemical mass balance model; Calculate the analytical sub-error for each pollution component, including: Among them, SSE j is the analytical sub-error corresponding to the j-th pollution component; N is the number of surrounding pollution sources; is the volume concentration of the contribution of the i-th pollution source to the j-th pollution component in the receptor pollution data in the analytical results; Volume concentration of the jth component in the receptor pollution data.

2. The method according to claim 1, characterized in that Establish a fingerprint of surrounding pollution sources, including: The pollution emission sources that have a pollution impact on the receiving area are referred to as the peripheral pollution sources; A fingerprint map is constructed based on the main pollutant components diffused from the pollution source to the receptor area and their volume concentration ratios, or a fingerprint map is constructed based on the main pollutant components emitted by the pollution source and their volume concentration ratios.

3. The method according to claim 2, characterized in that Obtain receptor contamination data for polluted areas, including: obtaining an ambient air sample in the receptor area; The ambient air sample is analyzed using a pollutant component qualitative and quantitative analysis device to obtain the pollutant components and volume concentration ratios.

4. The method according to claim 1, wherein The method further includes: using the error factor to correct the fingerprint spectrum of the surrounding pollution sources.

5. The method according to claim 1, wherein The variable weight chemical mass balance model is used to determine pollution sources, including: Based on the principle of minimizing the analytical error SSE of the variable weight chemical mass balance model, the weight coefficient to be optimized is calculated. , and further analyze to obtain the total volume concentration of the pollution components contributed by the i-th pollution source to the receptor area S i , and then determine the source of pollution.

6. The method according to claim 5, characterized in that Based on the calculation results of the variable weight chemical mass balance model, the contribution concentration of the i-th pollution source to the j-th pollution component in the receptor area is obtained. , the steps are as follows: ； in, is the volume concentration ratio of the jth pollution component in the fingerprint of the i-th pollution source; is the weight coefficient to be optimized for the jth pollution component of the i-th pollution source; is the total volume concentration of the pollution components contributed by the i-th pollution source to the receptor area; 、 This is the traceability result based on the variable weight chemical mass balance model.

7. A gas pollution leakage tracing device under complex diffusion conditions, used to implement the gas pollution leakage tracing method under complex diffusion conditions as claimed in any one of claims 1 to 6, characterized in that: The device comprises: The typical factor determination module establishes a fingerprint of surrounding pollution sources and sets the typical factor of each pollution source; the typical factor is used to characterize the typical pollution components of the pollution source; Data acquisition module, used to obtain receptor pollution data in the polluted area; an error factor determination module, configured to calculate an analytical error between the traceability result based on the chemical mass balance model and the receptor pollution data, and determine a pollutant component as an error factor based on its contribution to the analytical error; and A model building module is used to determine a variable weight chemical mass balance model based on the typical factors and the error factors, and the variable weight chemical mass balance model is used to determine the pollution source.

8. An electronic device, characterized in that: include: at least one processor; and a memory communicatively coupled to the at least one processor; Wherein, the memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor executes the gas pollution leakage tracing method under complex diffusion conditions as described in any one of claims 1 to 6.

9. A non-transitory computer-readable storage medium, characterized in that The non-transitory computer-readable storage medium stores computer-executable instructions, and the computer-executable instructions are used to enable the computer to execute the gas pollution leakage tracing method under complex diffusion conditions as described in any one of claims 1 to 6.

Images