CN117808189A - Overhead source pollutant concentration calculation method and device and readable storage medium - Google Patents
Overhead source pollutant concentration calculation method and device and readable storage medium Download PDFInfo
- Publication number
- CN117808189A CN117808189A CN202311525230.5A CN202311525230A CN117808189A CN 117808189 A CN117808189 A CN 117808189A CN 202311525230 A CN202311525230 A CN 202311525230A CN 117808189 A CN117808189 A CN 117808189A
- Authority
- CN
- China
- Prior art keywords
- calculated
- height
- grid
- source
- pollutant concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 147
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 147
- 238000003860 storage Methods 0.000 title claims abstract description 12
- 238000004364 calculation method Methods 0.000 title claims description 24
- 238000012937 correction Methods 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000009826 distribution Methods 0.000 claims abstract description 24
- 239000000779 smoke Substances 0.000 claims description 45
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 33
- 239000003546 flue gas Substances 0.000 claims description 33
- 239000000356 contaminant Substances 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000007599 discharging Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 239000007789 gas Substances 0.000 claims description 2
- 239000003517 fume Substances 0.000 claims 1
- 238000005260 corrosion Methods 0.000 abstract description 11
- 230000007797 corrosion Effects 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 15
- 230000005540 biological transmission Effects 0.000 description 14
- 230000009466 transformation Effects 0.000 description 12
- 238000009792 diffusion process Methods 0.000 description 8
- 230000007613 environmental effect Effects 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 241000220317 Rosa Species 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000011835 investigation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000011109 contamination Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to an overhead source pollutant concentration calculating method and device and a readable storage medium, and belongs to the technical field of pollutant emission in industrial areas. Taking a projection point of an overhead source of an industrial area on the ground as a central point, constructing polar coordinate systems at different heights in the vertical direction of the central point, and dispersing each polar coordinate system to obtain a plurality of grids to be calculated; calculating the pollutant concentration at each grid to be calculated based on the emission intensity of the emission source, the equivalent source height of the emission source, the average wind speed at the chimney outlet and the polar coordinates of each grid to be calculated; and correcting the calculated pollutant concentration by using a pollutant concentration correction model, outputting the corrected pollutant concentration at each grid to be calculated, and obtaining the pollutant concentration grade at each height based on the corrected pollutant concentrations at all grids to be calculated at each height. According to the method, different corrosion protection measures can be adopted for the overhead sources with different heights by acquiring the distribution condition of the pollutant concentration in space.
Description
Technical Field
The invention relates to the technical field of pollutant emission in industrial areas, in particular to an overhead source pollutant concentration calculation method, an overhead source pollutant concentration calculation device and a computer readable storage medium.
Background
Because the power transmission and transformation engineering equipment is mostly directly exposed in the atmosphere, the power transmission and transformation equipment can be corroded, so that the mechanical property of the equipment is reduced, and the service life and the reliability of the equipment are seriously influenced. The waste emission of the industrial area is an important cause for corrosion of power transmission and transformation equipment in the area, so that the research on the pollutant emission distribution rule of the industrial area provides reference for corrosion protection classification of the power transmission and transformation equipment in the area, adopts targeted protection to the power transmission and transformation equipment, and has important significance for improving the reliability and service life of the power transmission and transformation equipment.
The existing method for obtaining the concentration of the pollutants in the industrial area is to measure the concentration of the pollutants in the industrial area by using a gas concentration continuous measuring instrument or a sampler at a certain sampling point, wherein the measuring period is usually several months or one year, and finally, the average value of the concentration of the pollutants measured in the measuring period is used as the concentration of the pollutants in the industrial area. However, because the concentration of the pollutants is influenced by environmental factors such as weather, wind speed and the like, the conventional method directly takes the measurement average value in the measurement period as the concentration of the pollutants, and the influence of different environmental factors on the concentration of the pollutants is not considered, so that the direct taking of the average value in the measurement period as the concentration of the pollutants in an industrial area is not accurate enough; in addition, the existing measuring method cannot acquire the distribution condition of the pollutant concentration in space only by sampling at one sampling point, the pollutant concentrations of different areas with different heights in an industrial area are different, and the heights of power transmission and transformation equipment in the industrial area are not completely the same, so that the same corrosion protection measures are adopted for all the power transmission and transformation equipment in the industrial area based on the pollutant concentration level obtained by one sampling point, the corrosion protection measures for the power transmission and transformation equipment in the area with higher pollutant concentration level are insufficient, and the resource waste is caused by the adoption of the corrosion protection measures with too high level for the power transmission and transformation equipment in the area with lower pollutant concentration level, and the power transmission and transformation equipment in the area with different heights in the industrial area cannot be effectively corroded and protected.
In summary, the existing overhead source pollutant concentration obtaining method does not consider the problem that the obtained pollutant concentration result is inaccurate due to the influence of environmental factors, and the distribution situation of the pollutant concentration in space cannot be obtained only by sampling on sampling points, so that corresponding corrosion protection measures cannot be adopted for power transmission equipment with different heights in an industrial area according to the pollutant concentrations with different heights.
Disclosure of Invention
Therefore, the technical problem to be solved by the invention is to solve the problems that the acquired pollutant concentration result is inaccurate due to the influence of environmental factors not considered in the method for acquiring the overhead source pollutant concentration in the prior art, and the distribution condition of the pollutant concentration in space cannot be obtained only by sampling on sampling points, so that corresponding corrosion protection measures cannot be adopted for power transmission equipment with different heights in an industrial area according to the pollutant concentrations with different heights.
In order to solve the technical problems, the invention provides a method for calculating the concentration of an overhead source pollutant, which comprises the following steps:
taking a projection point of an overhead source in an industrial area on the ground as a central point, constructing polar coordinate systems at different heights in the vertical direction of the central point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and obtaining the polar coordinates of each grid to be calculated;
acquiring the emission intensity of an emission source, the height of a chimney for discharging pollutants, the lifting height of smoke and the average wind speed at a chimney outlet, and calculating the equivalent source height of the emission source based on the height of the chimney and the lifting height of smoke;
calculating a pollutant concentration at each grid to be calculated based on an emission intensity of the emission source, an equivalent source height of the emission source, an average wind speed at the stack outlet, and a polar coordinate of each grid to be calculated;
and constructing a pollutant concentration correction model based on wind directions, atmospheric stability and rainfall correction coefficients in the industrial area, inputting the pollutant concentration of each grid to be calculated into the pollutant concentration correction model, outputting the corrected pollutant concentration of each grid to be calculated, and obtaining the pollutant concentration level at each height in the industrial area based on the corrected pollutant concentrations of all grids to be calculated on a polar coordinate system corresponding to each height.
In one embodiment of the present invention, the calculation formula of the concentration of the contaminant at the grid to be calculated is:
wherein ρ (n) is the pollutant concentration at the nth grid to be calculated, and x is the polar coordinate of the grid to be calculatedThe space coordinate obtained after conversion is that y is the space coordinate obtained after coordinate conversion of the polar coordinate of the grid to be calculated, z is the height of the polar coordinate system of the grid to be calculated from the ground, Q is the emission intensity of the emission source,for the average wind speed at the outlet of the chimney, H is the equivalent source height, sigma y To calculate standard deviation, sigma, of the distribution of contaminants in the horizontal direction at the grid z Exp is an exponential function, which is the standard deviation of the distribution of the contaminants in the vertical direction at the grid to be calculated.
In one embodiment of the invention, the average wind speed at the chimney exitThe calculation formula of (2) is as follows:
wherein,for wind speed at 10m from the ground, H S Is the height of the chimney.
In one embodiment of the invention, the formula of the contaminant concentration correction model is:
ρ(n)′＝k·∑ i (f i ·∑ j ρ(n)·w ij )，
wherein ρ (n)' is the concentration of the pollutant corrected at the nth grid to be calculated, k is the rainfall correction coefficient of the industrial area, f i For the frequency of occurrence of wind direction i within a preset time of an industrial area, ρ (n) is the pollutant concentration at the nth grid to be calculated, w ij The frequency of occurrence of the atmospheric stability j in the wind direction i within the preset time of the industrial area is determined.
In one embodiment of the present invention, the flue gas elevation obtaining process includes:
acquiring the inner diameter of a chimney outlet, the average wind speed at the chimney outlet, the smoke emission rate, the smoke outlet flow velocity and the smoke outlet temperature;
and calculating a smoke heat release rate based on the smoke emission rate and the smoke outlet temperature, and calculating a smoke lifting height based on the chimney height, the smoke heat release rate, an average wind speed at the chimney outlet, the smoke outlet flow speed and the chimney outlet inner diameter.
In one embodiment of the present invention, the formula for calculating the smoke heat release rate is:
Q H ＝C p V 0 ΔT，
wherein Q is H C is the heat release rate of the smoke p Specific heat of flue gas at average constant pressure, V 0 For the flue gas emission rate, Δt is the difference between the flue gas outlet temperature and the ambient temperature.
In one embodiment of the invention, calculating the flue gas heat release rate based on the flue gas emission rate and the flue gas outlet temperature, calculating the flue gas elevation based on the stack height, the flue gas heat release rate, an average wind speed at the stack outlet, the flue gas outlet flow rate, and the stack outlet inner diameter comprises:
when Q is H More than or equal to 21000, and when the delta T is more than or equal to 35K,
when 2100 is less than or equal to Q H When the temperature is less than or equal to 21000 and the delta T is more than or equal to 35K,
when Q is H < 2100, and DeltaT < 35K,
wherein DeltaH 1 For smoke elevation in industrial areas located in cities or hills, ΔH 2 For the flue gas lifting height in the industrial area in plain or rural area, H S For the height of the chimney,for average wind speed at chimney exit, V S And d is the inner diameter of the outlet of the chimney.
In one embodiment of the present invention, taking a projection point of an overhead source in an industrial area on the ground as a center point, constructing polar coordinate systems at different heights in a vertical direction of the center point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height, obtaining a plurality of grids to be calculated of each polar coordinate system, and obtaining polar coordinates of each grid to be calculated includes:
taking a projection point of an overhead source in an industrial area on the ground as a central point, taking a vertical upward direction of the central point as a z axis, taking points on the z axis at different heights from the ground as poles, taking a forward direction of a discharge source as a polar axis, and taking a counterclockwise direction as a forward direction to construct a plurality of polar coordinate systems respectively corresponding to different heights;
and discretizing the polar axis and the polar angle of the polar coordinate system corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, obtaining the polar coordinates of each grid to be calculated on the corresponding polar coordinate system, and taking the height of the polar coordinate system of each grid to be calculated from the ground as the z-axis coordinate of each grid to be calculated.
The invention also provides an overhead source pollutant concentration calculating device, which comprises:
the coordinate system construction module is used for constructing polar coordinate systems at different heights in the vertical direction of a central point by taking the projection point of an overhead source in an industrial area on the ground as the central point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and acquiring the polar coordinates of each grid to be calculated;
the equivalent source height calculation module is used for acquiring the emission intensity of the emission source, the height of a chimney for discharging pollutants, the lifting height of the smoke and the average wind speed at the outlet of the chimney, and calculating the equivalent source height of the emission source based on the height of the chimney and the lifting height of the smoke;
a pollutant concentration calculation module for calculating a pollutant concentration at each grid to be calculated based on an emission intensity of the emission source, an equivalent source height of the emission source, an average wind speed at the chimney exit, and polar coordinates of the grids to be calculated;
the pollutant concentration correction module is used for constructing a pollutant concentration correction model based on wind directions, atmospheric stability and rainfall correction coefficients in the industrial area, inputting the pollutant concentration of each grid to be calculated into the pollutant concentration correction model, outputting the corrected pollutant concentration of each grid to be calculated, and obtaining the pollutant concentration level of each height in the industrial area based on the corrected pollutant concentrations of all grids to be calculated on the polar coordinate system corresponding to each height.
The invention also provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and the computer program realizes the steps of the method for calculating the concentration of the overhead source pollutant when being executed by a processor.
According to the method for calculating the concentration of the overhead source pollutants, provided by the invention, the projection point of the overhead source in the industrial area on the ground is taken as a central point, a polar coordinate system is constructed at different heights in the vertical upward direction of the point, polar axes and polar angles of the polar coordinate system corresponding to each height are scattered to obtain a plurality of grids to be calculated, the grids obtained after each polar coordinate is scattered correspond to a plurality of areas at the same height, the concentration of the pollutants in the different areas at different heights can be calculated by carrying out discrete division on the areas in the industrial area, so that the distribution situation of the concentration of the pollutants in space is obtained, and different corrosion protection measures are adopted for the overhead source at different heights in the industrial area based on the concentration levels of the pollutants at different heights; in addition, the pollutant concentration of each grid to be calculated is calculated based on the emission intensity of the emission source, the equivalent source height of the emission source and the average wind speed at the outlet of the chimney, and then the pollutant concentration of each grid to be calculated in the industrial area is corrected based on the wind direction, the atmospheric stability and the rainfall correction coefficient in the industrial area, so that the influence of environmental factors such as weather, wind speed and the like on the pollutant concentration is considered, and the calculated pollutant concentration is more accurate.
Drawings
In order that the invention may be more readily understood, a more particular description of the invention will be rendered by reference to specific embodiments thereof that are illustrated in the appended drawings, in which
FIG. 1 is a schematic flow chart of a method for calculating the concentration of an overhead source pollutant;
FIG. 2 is a diagram of a computed overhead source SO provided by an embodiment of the invention 2 Schematic diagram of concentration;
FIG. 3 shows an overhead source SO provided in an embodiment of the invention 2 Schematic diagram of concentration correction;
FIG. 4 is a schematic diagram of an industrial area overhead source contaminant concentration research area provided by an embodiment of the present invention;
FIG. 5 is a schematic view of the area of interest shown in FIG. 4 after discretizing the area of interest;
FIG. 6 is a rose diagram of the distribution of primary contamination zones within the investigation region shown in FIG. 4; wherein (a) in fig. 6 is a distribution rose of the first-order pollution area in the investigation region at 0m from the ground, (b) in fig. 6 is a distribution rose of the first-order pollution area in the investigation region at 50m from the ground, and (c) in fig. 6 is a distribution rose of the first-order pollution area in the investigation region at 100m from the ground;
fig. 7 is a schematic structural diagram of an overhead source pollutant concentration calculating device provided by the invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings and specific examples, which are not intended to be limiting, so that those skilled in the art will better understand the invention and practice it.
Referring to fig. 1, fig. 1 is a flowchart of an overhead source pollutant concentration calculation method according to an embodiment of the present application, which specifically includes:
s10: constructing polar coordinate systems at different heights in the vertical direction of a central point by taking projection points of overhead sources in an industrial area on the ground as the central point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and acquiring the polar coordinates of each grid to be calculated;
s20: acquiring the emission intensity of an emission source, the height of a chimney for discharging pollutants, the lifting height of smoke and the average wind speed at the outlet of the chimney, and calculating the equivalent source height of the emission source based on the height of the chimney and the lifting height of smoke;
s30: calculating the pollutant concentration at each grid to be calculated based on the emission intensity of the emission source, the equivalent source height of the emission source, the average wind speed at the chimney outlet and the polar coordinates of each grid to be calculated;
s40: and constructing a pollutant concentration correction model based on wind directions, atmospheric stability and rainfall correction coefficients in the industrial area, inputting the pollutant concentration of each grid to be calculated into the pollutant concentration correction model, outputting the corrected pollutant concentration of each grid to be calculated, and obtaining the pollutant concentration level of each height in the industrial area based on the corrected pollutant concentrations of all grids to be calculated on a polar coordinate system corresponding to each height.
According to the method for calculating the concentration of the overhead source pollutants, the projection point of the overhead source in the industrial area on the ground is taken as a central point, a polar coordinate system is constructed at different heights in the vertical upward direction of the point, polar axes and polar angles of the polar coordinate system corresponding to each height are scattered to obtain a plurality of grids to be calculated, the grids obtained after the scattering of each polar coordinate correspond to a plurality of areas at the same height, the concentration of the pollutants in the different areas at different heights can be calculated by carrying out discrete division on the areas in the industrial area, so that the distribution situation of the concentration of the pollutants in space is obtained, and different corrosion protection measures are adopted for power transmission and transformation equipment at different heights in the industrial area based on the concentration levels of the pollutants at different heights; in addition, the pollutant concentration of each grid to be calculated is calculated based on the emission intensity of the emission source, the equivalent source height of the emission source and the average wind speed at the outlet of the chimney, and then the pollutant concentration of each grid to be calculated in the industrial area is corrected based on the wind direction, the atmospheric stability and the rainfall correction coefficient in the industrial area, so that the influence of environmental factors such as weather, wind speed and the like on the pollutant concentration is considered, and the calculated pollutant concentration is more accurate.
Specifically, in some embodiments of the present application, step S10 specifically includes:
s100: taking a projection point of an overhead source in an industrial area on the ground as a central point, taking a vertical upward direction of the central point as a z axis, taking points on the z axis at different heights from the ground as poles, taking a forward direction of a discharge source as a polar axis, and constructing a plurality of polar coordinate systems respectively corresponding to different heights in a forward direction;
illustratively, taking points 50m and 100m from the ground as a first pole and a second pole on the z-axis, and constructing a polar coordinate system at a height of 50m from the ground based on the first pole with the positive east direction of the emission source as the polar axis and the counterclockwise direction as the positive direction; based on the second pole, a polar coordinate system at a height of 100m from the ground is constructed with the positive east direction of the emission source as the polar axis and the counterclockwise direction as the forward direction.
S110: and discretizing the polar axis and the polar angle of the polar coordinate system corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, obtaining the polar coordinates of each grid to be calculated on the corresponding polar coordinate system, and taking the height of the polar coordinate system of each grid to be calculated from the ground as the z-axis coordinate of each grid to be calculated.
By way of example, as a specific example of the present application, the polar axis r of each polar coordinate system is 20km to cover the whole industrial area, the polar axis r of each polar coordinate system is discretized according to 0.1km for each grid, the polar angle θ of each polar coordinate system is discretized according to 0.5 ° for each grid, a plurality of fan rings with a central angle of 0.5 ° are obtained after each polar coordinate system is discretized, each fan ring is used as one grid to be calculated, and the number of grids to be calculated on each polar coordinate system is 144000.
Further, the process of acquiring the flue gas elevation in step S20 includes:
acquiring the inner diameter of a chimney outlet, the average wind speed at the chimney outlet, the smoke emission rate, the smoke outlet flow velocity and the smoke outlet temperature;
the flue gas heat release rate is calculated based on the flue gas discharge rate and the flue gas outlet temperature, and the flue gas elevation height is calculated based on the chimney height, the flue gas heat release rate, the average wind speed at the chimney outlet, the flue gas outlet flow speed and the chimney outlet inner diameter.
Specifically, the calculation formula of the smoke heat release rate is as follows:
Q H ＝C p V 0 ΔT，
wherein Q is H C is the heat release rate of the smoke p The specific heat of the flue gas is equal to the average constant pressure, and the numerical value is 1.38 KJ/(m) 3 ·K)，V 0 Delta T is the difference between the temperature of the flue gas outlet and the ambient temperature;
wherein Δt=t S -T α ，T S T is the temperature of the flue gas outlet α Is the average temperature of the environment;
T S ＝T C -(H S -H C )*(5/100)，H S is the height of chimney, H C T is the height of flue gas from the ground at any point in the chimney C Is H in chimney C The flue gas temperature is set;
further, the calculation formula of the flue gas lifting height is as follows:
when Q is H More than or equal to 21000, and when the delta T is more than or equal to 35K,
when 2100 is less than or equal to Q H When the temperature is less than or equal to 21000 and the delta T is more than or equal to 35K,
when Q is H < 2100, and DeltaT < 35K,
wherein DeltaH 1 For smoke elevation in industrial areas located in cities or hills, ΔH 2 For the flue gas lifting height in the industrial area in plain or rural area, H S For the height of the chimney,for average wind speed at chimney exit, V S The flow rate of the flue gas outlet is d, and the inner diameter of the chimney outlet is d;
further, the calculation formula of the equivalent source height of the emission source is:
H＝H S +ΔH，
wherein H is the equivalent source height of the emission source, and delta H is the flue gas lifting height.
Specifically, calculating the pollutant concentration at each grid to be calculated in step S30 based on the emission intensity of the emission source, the equivalent source height of the emission source, the average wind speed at the stack outlet, and the polar coordinates of each grid to be calculated includes:
obtaining an overhead continuous point source Gaussian diffusion model based on the continuous point source high-low diffusion model in the unbounded space;
specifically, the formula of the continuous point source gaussian diffusion model in the unbounded space is as follows:
where ρ (x, y, z) represents the concentration of the contaminant at the grid of spatial coordinates (x, y, z), Q represents the emission intensity of the emission source,representing average wind speed, sigma, at the chimney exit y Represents the standard deviation, sigma, of the distribution of contaminants in the horizontal direction at the grid z The standard deviation of the distribution of pollutants at the grid in the vertical direction is represented, and exp is an exponential function;
because the difference between the overhead continuous point source and the unbounded space needs to consider the reflection effect of the ground on the pollutant, the application improves the continuous point source Gaussian diffusion model in the unbounded space to obtain the overhead continuous point source Gaussian diffusion model;
specifically, in the embodiment of the present application, the formula of the overhead continuous point source gaussian diffusion model is:
wherein ρ (n) represents the pollutant concentration at the nth grid to be calculated, x, y, z are the spatial coordinates of the grid to be calculated, respectively, Q is the emission intensity of the emission source,for the average wind speed at the outlet of the chimney, H is the equivalent source height, sigma y To calculate standard deviation, sigma, of the distribution of contaminants in the horizontal direction at the grid z For the standard deviation of the distribution of pollutants at the grid to be calculated in the vertical direction, exp is an exponential function;
inputting the emission intensity of the emission source, the equivalent source height of the emission source, the average wind speed at the outlet of the chimney and the polar coordinates of each grid to be calculated into an overhead continuous point source high-low diffusion model to obtain the pollutant concentration at each grid to be calculated;
because the coordinates constructed in the embodiment of the application are polar coordinates, when the coordinates of the grids to be calculated are input into the model to calculate the concentration of the pollutant, coordinate conversion is needed, that is, x and y are space coordinates obtained after coordinate conversion is respectively performed on the polar coordinates of the grids to be calculated, and z is the height of the polar coordinate system of the grids to be calculated from the ground, that is, the z-axis coordinate of each grid to be calculated.
Preferably, the average wind speed at the outlet of the chimney in the embodiment of the applicationThe calculation formula of (2) is as follows:
wherein,for wind speed at 10m from the ground, H S Is the height of the chimney.
Because the wind power grade can be more accurately determined by utilizing the wind speed which is 10m higher than the ground, the influence of wind power on pollutant diffusion can be reflected, and data support is provided for environmental protection, the pollutant concentration can be more accurately calculated by calculating the average wind speed at the outlet of the chimney through the wind speed which is 10m higher than the ground in the embodiment of the application.
Further, since the influence of the meteorological environmental factors on the pollutant concentration is not considered in the calculation of the pollutant concentration in the above embodiment, the pollutant concentration correction model is constructed according to the wind direction, the atmospheric stability and the rainfall correction coefficient in the industrial area in the embodiment of the present application, specifically, the formula of the pollutant concentration correction model in step S40 is as follows:
ρ(n)′＝k·∑ i (f i ·∑ j ρ(n)·w ij )，
wherein ρ (n)' is the concentration of the pollutant corrected at the nth grid to be calculated, k is the rainfall correction coefficient of the industrial area, i.e. the frequency of occurrence of non-rainy and snowy weather in preset time, f i For a preset time in an industrial areaThe frequency of occurrence of wind direction i, ρ (n) is the concentration of contaminants at the nth grid to be calculated, w ij The frequency of occurrence of the atmospheric stability j in the wind direction i within the preset time of the industrial area is determined.
Below to calculate SO of industrial area overhead source 2 Concentration as an example the above examples are further explained:
referring to fig. 2, fig. 2 shows a computing overhead source SO according to the present embodiment 2 Schematic concentration diagram, firstly, calculating SO at each grid to be calculated based on emission condition analysis 2 Concentration and then calculated SO by using meteorological environment factors 2 Correcting the concentration, and finally, correcting SO based on all grids to be calculated 2 The concentration of SO at different heights 2 Concentration profile, thereby obtaining SO at different heights 2 Concentration grade.
Referring to FIG. 3, FIG. 3 shows the embodiment of the present invention, which provides a method for utilizing weather environment factors to SO 2 The concentration correction schematic diagram, namely the construction principle of the pollutant concentration correction model in the application, is based on the calculated SO (SO) on the basis of different atmospheric stabilities under a single wind direction 2 The concentration is corrected once, and then SO after the primary correction is carried out based on the influences of different wind directions 2 The concentration was corrected twice.
Referring to fig. 4, fig. 4 shows three polar coordinate systems with different heights constructed in this embodiment, and pollutant concentration levels of different areas with different heights can be obtained by calculating pollutant concentrations in the ranges of the three polar coordinate systems; the first polar coordinate system is used for calculating the pollutant concentration at the position with the height of 0m from the ground, the second polar coordinate system is used for calculating the pollutant concentration at the position with the height of 50m from the ground, and the third polar coordinate system is used for calculating the pollutant concentration at the position with the height of 100m from the ground.
Referring to fig. 5, fig. 5 is a schematic view of a region obtained by discretizing the polar coordinate system shown in fig. 4, in which the whole polar coordinate system is divided into 8 large regions according to directions.
SO of different areas in the polar coordinate system corresponding to different heights based on the calculation 2 The concentration can be variedDistribution results of primary contamination zone at height:
referring to fig. 6, fig. 6 shows a first-level pollution-area distribution rose diagram of different areas in three polar coordinate systems calculated by the method provided by the embodiment of the present application; wherein (a) in the figure represents a first-stage pollution area distribution rose at a height of 0m from the ground, (b) in the figure represents a first-stage pollution area distribution rose at a height of 50m from the ground, and (c) in the figure represents a first-stage pollutant distribution rose at a height of 100m from the ground; the gray scale part in the graph is a first-level pollution area, the pollutant concentration levels of different areas with different heights can be obtained based on the distribution rose, and corresponding corrosion protection measures are adopted for the overhead sources of the areas according to the pollutant concentration levels.
Based on the method for calculating the concentration of the overhead source pollutant provided in the foregoing embodiment, the embodiment of the present application further provides an apparatus for calculating the concentration of the overhead source pollutant, as shown in fig. 7, which specifically includes:
the coordinate system construction module 10 is configured to construct a polar coordinate system at different heights in a vertical direction of a center point by taking a projection point of an overhead source in an industrial area on the ground as the center point, and discrete polar axes and polar angles of the polar coordinate system corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and acquire polar coordinates of each grid to be calculated;
the equivalent source height calculation module 20 is configured to obtain an emission intensity of the emission source, a chimney height for discharging pollutants, a flue gas elevation height, and an average wind speed at an outlet of the chimney, and calculate an equivalent source height of the emission source based on the chimney height and the flue gas elevation height;
a pollutant concentration calculation module 30 for calculating the pollutant concentration at each grid to be calculated based on the emission intensity of the emission source, the equivalent source height of the emission source, the average wind speed at the chimney exit and the polar coordinates of the grid to be calculated;
the pollutant concentration correction module 40 is configured to construct a pollutant concentration correction model based on wind direction, atmospheric stability and rainfall correction coefficient in the industrial area, input the pollutant concentration at each grid to be calculated into the pollutant concentration correction model, output the corrected pollutant concentration at each grid to be calculated, and obtain the pollutant concentration level at each height in the industrial area based on the corrected pollutant concentrations at all grids to be calculated on the polar coordinate system corresponding to each height.
Embodiments of the present application also provide a computer readable storage medium having a computer program stored thereon, which when executed by a processor, implements the steps of the above-described elevated source contaminant concentration calculation method.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may 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 instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may 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 which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It is apparent that the above examples are given by way of illustration only and are not limiting of the embodiments. Other variations and modifications of the present invention will be apparent to those of ordinary skill in the art in light of the foregoing description. It is not necessary here nor is it exhaustive of all embodiments. While still being apparent from variations or modifications that may be made by those skilled in the art are within the scope of the invention.
Claims (10)
1. A method of calculating an elevated source contaminant concentration, comprising:
taking a projection point of an overhead source in an industrial area on the ground as a central point, constructing polar coordinate systems at different heights in the vertical direction of the central point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and obtaining the polar coordinates of each grid to be calculated;
acquiring the emission intensity of an emission source, the height of a chimney for discharging pollutants, the lifting height of smoke and the average wind speed at a chimney outlet, and calculating the equivalent source height of the emission source based on the height of the chimney and the lifting height of smoke;
calculating a pollutant concentration at each grid to be calculated based on an emission intensity of the emission source, an equivalent source height of the emission source, an average wind speed at the stack outlet, and a polar coordinate of each grid to be calculated;
and constructing a pollutant concentration correction model based on wind directions, atmospheric stability and rainfall correction coefficients in the industrial area, inputting the pollutant concentration of each grid to be calculated into the pollutant concentration correction model, outputting the corrected pollutant concentration of each grid to be calculated, and obtaining the pollutant concentration level at each height in the industrial area based on the corrected pollutant concentrations of all grids to be calculated on a polar coordinate system corresponding to each height.
2. The method of claim 1, wherein the formula for calculating the concentration of the contaminant at the grid to be calculated is:
wherein ρ (n) is the pollutant concentration at the nth grid to be calculated, x is the space coordinate obtained after coordinate conversion of the polar coordinate of the grid to be calculated, y is the space coordinate obtained after coordinate conversion of the polar coordinate of the grid to be calculated, z is the height of the polar coordinate system of the grid to be calculated from the ground, Q is the emission intensity of the emission source,for the average wind speed at the outlet of the chimney, H is the equivalent source height, sigma y To calculate standard deviation, sigma, of the distribution of contaminants in the horizontal direction at the grid z Exp is an exponential function, which is the standard deviation of the distribution of the contaminants in the vertical direction at the grid to be calculated.
3. The method of claim 2, wherein the average wind speed at the stack outlet is the average wind speed at the stack outletThe calculation formula of (2) is as follows:
wherein,for wind speed at 10m from the ground, H S Is the height of the chimney.
4. The method of claim 1, wherein the formula of the pollutant concentration correction model is:
ρ(n)′＝k·∑ i (f i ·∑ j ρ(n)·w ij )，
wherein ρ (n)' is the concentration of the pollutant corrected at the nth grid to be calculated, k is the rainfall correction coefficient of the industrial area, f i For the frequency of occurrence of wind direction i within a preset time of an industrial area, ρ (n) is the pollutant concentration at the nth grid to be calculated, w ij The frequency of occurrence of the atmospheric stability j in the wind direction i within the preset time of the industrial area is determined.
5. The method of claim 1, wherein the step of obtaining the elevation of the flue gas comprises:
acquiring the inner diameter of a chimney outlet, the average wind speed at the chimney outlet, the smoke emission rate, the smoke outlet flow velocity and the smoke outlet temperature;
and calculating a smoke heat release rate based on the smoke emission rate and the smoke outlet temperature, and calculating a smoke lifting height based on the chimney height, the smoke heat release rate, an average wind speed at the chimney outlet, the smoke outlet flow speed and the chimney outlet inner diameter.
6. The method for calculating the concentration of the elevated source pollutant according to claim 5, wherein the formula for calculating the heat release rate of the flue gas is:
Q H ＝C p V 0 ΔT，
wherein Q is H C is the heat release rate of the smoke p Specific heat of flue gas at average constant pressure, V 0 Is the fume emission rate, delta T is fumeThe difference between the gas outlet temperature and the ambient temperature.
7. The overhead source pollutant concentration calculation method of claim 6, wherein calculating a smoke heat release rate based on the smoke emission rate and the smoke outlet temperature, and calculating a smoke elevation based on the chimney height, the smoke heat release rate, an average wind speed at the chimney outlet, the smoke outlet flow rate, and the chimney outlet inner diameter comprises:
when Q is H More than or equal to 21000, and when the delta T is more than or equal to 35K,
when 2100 is less than or equal to Q H When the temperature is less than or equal to 21000 and the delta T is more than or equal to 35K,
when Q is H < 2100, and DeltaT < 35K,
wherein DeltaH 1 For smoke elevation in industrial areas located in cities or hills, ΔH 2 For the flue gas lifting height in the industrial area in plain or rural area, H S For the height of the chimney,for average wind speed at chimney exit, V S And d is the inner diameter of the outlet of the chimney.
8. The method of claim 1, wherein constructing polar coordinate systems at different heights in a vertical direction of a center point of a projection point of an overhead source on the ground in an industrial area, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height, obtaining a plurality of grids to be calculated of each polar coordinate system, and obtaining polar coordinates of each grid to be calculated comprises:
taking a projection point of an overhead source in an industrial area on the ground as a central point, taking a vertical upward direction of the central point as a z axis, taking points on the z axis at different heights from the ground as poles, taking a forward direction of a discharge source as a polar axis, and taking a counterclockwise direction as a forward direction to construct a plurality of polar coordinate systems respectively corresponding to different heights;
and discretizing the polar axis and the polar angle of the polar coordinate system corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, obtaining the polar coordinates of each grid to be calculated on the corresponding polar coordinate system, and taking the height of the polar coordinate system of each grid to be calculated from the ground as the z-axis coordinate of each grid to be calculated.
9. An overhead source contaminant concentration calculation apparatus, comprising:
the coordinate system construction module is used for constructing polar coordinate systems at different heights in the vertical direction of a central point by taking the projection point of an overhead source in an industrial area on the ground as the central point, dispersing polar axes and polar angles of the polar coordinate systems corresponding to each height to obtain a plurality of grids to be calculated of each polar coordinate system, and acquiring the polar coordinates of each grid to be calculated;
the equivalent source height calculation module is used for acquiring the emission intensity of the emission source, the height of a chimney for discharging pollutants, the lifting height of the smoke and the average wind speed at the outlet of the chimney, and calculating the equivalent source height of the emission source based on the height of the chimney and the lifting height of the smoke;
a pollutant concentration calculation module for calculating a pollutant concentration at each grid to be calculated based on an emission intensity of the emission source, an equivalent source height of the emission source, an average wind speed at the chimney exit, and polar coordinates of the grids to be calculated;
the pollutant concentration correction module is used for constructing a pollutant concentration correction model based on wind directions, atmospheric stability and rainfall correction coefficients in the industrial area, inputting the pollutant concentration of each grid to be calculated into the pollutant concentration correction model, outputting the corrected pollutant concentration of each grid to be calculated, and obtaining the pollutant concentration level of each height in the industrial area based on the corrected pollutant concentrations of all grids to be calculated on the polar coordinate system corresponding to each height.
10. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the elevated source pollutant concentration calculation method of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311525230.5A CN117808189A (en) | 2023-11-15 | 2023-11-15 | Overhead source pollutant concentration calculation method and device and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311525230.5A CN117808189A (en) | 2023-11-15 | 2023-11-15 | Overhead source pollutant concentration calculation method and device and readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117808189A true CN117808189A (en) | 2024-04-02 |
Family
ID=90427644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311525230.5A Pending CN117808189A (en) | 2023-11-15 | 2023-11-15 | Overhead source pollutant concentration calculation method and device and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117808189A (en) |
-
2023
- 2023-11-15 CN CN202311525230.5A patent/CN117808189A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110909483B (en) | Point source atmospheric pollutant emission list verification method based on gridding data | |
| CN108648127B (en) | Method for locking urban air pollution hot spot area | |
| CN109117549B (en) | Pollutant diffusion prediction method and system | |
| CN108489875B (en) | Pollutant tracing system and method based on time period statistical analysis | |
| KR20140024541A (en) | System for providing atmospheric modeling system service | |
| CN116227749B (en) | Method and device for determining pollutant emission amount, storage medium and electronic equipment | |
| CN114329245A (en) | Air pollution tracing method and device, server and storage medium | |
| CN114662344B (en) | Atmospheric pollution source tracing prediction method and system based on continuous online observation data | |
| CN111157688A (en) | Method and device for evaluating influence of pollution source on air quality monitoring station | |
| KR101607878B1 (en) | System and method for providing air quality information by integrating virtual sensor and monitoring sensor | |
| CN113125635A (en) | Atmospheric pollution early warning method and device and readable storage medium | |
| CN114819704A (en) | Tracing method and system for mobile source pollutant emission | |
| CN112132450A (en) | Method for positioning gaseous pollutants | |
| CN115146484A (en) | Environment-friendly monitoring system and monitoring method for detecting environmental parameters | |
| CN113514612B (en) | Tracing method for pollution in area | |
| CN114002381A (en) | Atmospheric pollution source-tracing diffusion analysis method and device | |
| CN116258101B (en) | Method for rapidly calculating emission source intensity based on gas concentration monitoring data | |
| CN117455110A (en) | Carbon emission monitoring and evaluating method considering carbon emission point source diffusion | |
| CN117808189A (en) | Overhead source pollutant concentration calculation method and device and readable storage medium | |
| CN116109323A (en) | Ozone pollution tracing method, device, equipment and storage medium | |
| CN115705510A (en) | Factory gas pollution distribution prediction method and system, electronic equipment and storage medium | |
| KR20010097331A (en) | Information System Of Air Pollution Using Laser Radar And Dispersion model | |
| JP5131400B1 (en) | Search method of unsteady dust source position of falling dust | |
| Diez et al. | Uncertainty propagation of meteorological and emission data in modeling pollutant dispersion in the atmosphere | |
| CN112540052A (en) | Method and system for measuring and calculating pollutant discharge amount of storage tank to be measured |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |