CN118069958A - Atmospheric environment capacity calculation method and device - Google Patents
Atmospheric environment capacity calculation method and device Download PDFInfo
- Publication number
- CN118069958A CN118069958A CN202410465170.0A CN202410465170A CN118069958A CN 118069958 A CN118069958 A CN 118069958A CN 202410465170 A CN202410465170 A CN 202410465170A CN 118069958 A CN118069958 A CN 118069958A
- Authority
- CN
- China
- Prior art keywords
- target area
- station
- ground
- height
- monitoring data
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 41
- 238000012544 monitoring process Methods 0.000 claims abstract description 57
- 230000007613 environmental effect Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000003746 surface roughness Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 30
- 238000011158 quantitative evaluation Methods 0.000 abstract description 2
- 239000003344 environmental pollutant Substances 0.000 description 15
- 231100000719 pollutant Toxicity 0.000 description 15
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 239000012855 volatile organic compound Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000000356 contaminant Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a device for calculating atmospheric environment capacity, which are used for acquiring key parameters of meteorological monitoring data of a target area; calculating to obtain the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the near-ground station and the high-altitude station; calculating a plurality of groups of wind speed fields at different heights to obtain near-ground rational Charles coefficients; solving the boundary layer height of the target area; and calculating to obtain the atmospheric pollution environment capacity in the target area by combining the horizontal wind speeds of multiple points in the target area and the boundary layer height of the target area. The advantages are that: the method for real-time quantitative rapid calculation of the atmospheric environment capacity of the industrial park and the surrounding area by utilizing the real-time environment weather conventional monitoring data of the industrial park and the surrounding area which are easy to acquire; the method can provide reliable and rapid real-time calculation data of the atmospheric pollution emission amount for the space region where the industrial park is located, so that scientific quantitative evaluation results are provided for the influence of the industrial park on the surrounding environment.
Description
Technical Field
The invention relates to an atmospheric environment capacity calculation method and device, and belongs to the technical field of meteorological data processing.
Background
Atmospheric emissions from industrial parks, such as particulates, hazardous gases, and volatile organic compounds, often have a significant impact on the surrounding environment. And wherein the atmospheric environment capacity is the supportable amount of the atmospheric space and its diffusion and transmission physical characteristics in a certain area capable of holding the atmospheric contaminants while maintaining a certain concentration threshold. The scientific quantitative assessment of the atmospheric environment capacity of the industrial park and the surrounding area is an important scientific theoretical premise and management basis for overall management of the acceptable amount of the atmospheric pollution emission occurring in industrial traffic and other human activities of the industrial park.
Current methods for monitoring or accounting for emissions of specific pollutants from industrial parks have several principle and technical routes:
1. Material balancing method: knowing the amount of feed required for a particular process, knowing at the current state of the art what contaminants and corresponding emissions are produced per unit feed, the pollutant emissions can be estimated; the defects are that: 1) The feeding quantity cannot be monitored in real time by the environment or other management departments under the actual working condition, and only more accurate data in units of years or seasons can be obtained, so that the pollution emission can be monitored only by a value of a year-to-year scale, and the emission change characteristics (such as hour by hour) can not be monitored; 2) The emission intensity of a specific process has a relation number which is estimated according to a normal and complete process flow, but in the actual production process, a large number of illegal operations or other abnormal operations exist in an enterprise main body, so that a larger systematic deviation exists between a theoretical estimated value and an actual situation;
2. Discharge port flow monitoring method: installing monitoring equipment for pollutant concentration and emission flow at the discharge port of the equipment for generating pollutant emission, obtaining pollutant emission data in real time, and performing time integration on the pollutant emission data to obtain total emission; the disadvantage is that only fixed and organized discharge ports can be monitored. Under actual working conditions, most enterprises have a large amount of unorganized emission, such as pipeline running and leaking, or the private part of the enterprises is emitted from an abnormal terminal in order to bypass monitoring, and finally, the monitoring result of the method deviates greatly from the actual monitoring result;
3. Tracing model method: the emission list tracing is carried out by combining a mesoscale atmospheric diffusion mode with a regular matrix method, the spatial distribution of the concentration of the spatial atmospheric pollutants in the region and the time variation characteristics of the spatial atmospheric pollutants in the region are obtained through a numerical mode, and then inversion emission intensity calculation is carried out according to some assumptions. From the general theory, the invention also belongs to this type of method, but the former method has the following drawbacks: 1) At present, the method is based on a mesoscale atmospheric mode, namely resolution is kilometer level, and the simulation area range is hundred kilometers level; the typical space range of the industrial park is only a few kilometers, mainly the street and valley scale, so the precision of the existing mode is too thick; the main physical mechanisms of atmospheric diffusion are different in the mesoscale and the street-valley scale, namely the former related mode can not solve the problem by changing the resolution and the simulation range under the condition of not changing the mode core algorithm and the physical assumption; the smaller the region of interest, the greater the resulting error; 2) The mesoscale principle and method described in 1) requires repeated invocation and solving of the atmospheric diffusion equation in each real-time calculation stage, and the equation is large in calculation amount, so that the method requires a large amount of numerical calculation amount in each real-time calculation stage after superposition of two boxes, and therefore cannot be completed within a specified time (i.e. a 30-minute discharge amount is calculated, the required calculation time is far longer than 30 minutes, generally about 8 hours), and the meaning is lost, and continuous rolling operation is not achieved.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide an atmospheric environment capacity calculation method and device.
In order to solve the technical problems, the invention is realized by adopting the following technical scheme.
In one aspect, the invention provides a method for calculating the capacity of an atmospheric environment, which comprises the following steps:
Acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target area, and meteorological monitoring data of the near-ground station and the high-altitude station;
Performing self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain key parameters of the meteorological monitoring data;
Calculating to obtain the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the meteorological monitoring data;
Calculating wind speed fields on a plurality of groups of different heights according to the rough height of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to solve near-ground rational Charles coefficients;
Solving the boundary layer height of the target area according to the near-ground Consumer coefficient;
and acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating by combining the horizontal wind speeds of the multiple points in the target area and the boundary layer height of the target area to obtain the atmospheric pollution environment capacity in the target area.
Further, the calculation formula of the roughness height of the underlying surface of the target area is as follows:
(1);
Where Z 0 denotes the target area underlying roughness height, Z g denotes the near ground station height, Z h denotes the high altitude station height, B represents a key parameter.
Further, the calculation formula of the near-ground rational Charles coefficient is as follows:
(2);
Wherein R ic represents near-ground Concharrson, R 0 represents surface Raspy, g represents gravitational acceleration, b represents a key parameter, θ v0 represents virtual potential temperature, U z represents wind speed at a high-altitude station site, U g represents wind speed near-ground site, Z g represents near-ground site height, and Z h represents high-altitude station site height;
(3);
Where f 0 denotes the coriolis parameter.
Further, the calculation formula of the boundary layer height of the target area is as follows:
(4);
Wherein: Δu is the difference in wind speed in the vertical direction based on the observation of the park high-altitude station and the ground station; Δz is the difference between the monitored heights of the overhead and ground stations, U z is the wind speed at the top of the atmospheric boundary layer, U g is the wind speed at the ground monitoring station, pbl is the ground height at the top of the atmospheric boundary layer, b 1 is the gradient of the virtual temperature in the vertical direction, and Δθ v is the difference between the virtual temperatures in the vertical direction as observed by the ground and overhead stations.
Further, the calculation formula of the atmospheric pollution environment capacity in the target area is as follows:
(5);
Where E cap represents the atmospheric pollution environmental capacity within the target area, v i,j represents the real-time wind speed at different horizontal coordinates (i, j) within the target area, and pbl is the boundary layer height of the target area.
In a second aspect, the present invention provides an atmospheric environment capacity calculation device, comprising:
The acquisition module is used for acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target area and weather monitoring data of the near-ground station and the high-altitude station;
The fitting module is used for carrying out self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain key parameters of the meteorological monitoring data;
The calculation module is used for calculating and obtaining the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the meteorological monitoring data; calculating wind speed fields on a plurality of groups of different heights according to the rough height of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to solve near-ground rational Charles coefficients; solving the boundary layer height of the target area according to the near-ground Consumer coefficient; and acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating by combining the horizontal wind speeds of the multiple points in the target area and the boundary layer height of the target area to obtain the atmospheric pollution environment capacity in the target area.
The invention has the beneficial effects that:
The invention completes a method for quickly calculating the atmospheric environment capacity of the region in real time quantitatively based on the real-time environmental weather conventional monitoring data of the industrial park and the surrounding region which are easy to acquire. By further combining factors such as meteorological factors, atmospheric pollution concentration, near-ground atmospheric instability, height of a mixed layer and the like, reliable and rapid real-time atmospheric pollution emission calculation data can be provided for a space region where an industrial park is located, so that scientific quantitative evaluation results are provided for the influence of the industrial park on the surrounding environment
Drawings
FIG. 1 is a flow chart of an implementation of the present invention;
FIG. 2 is a schematic illustration of accounting for PM2.5 concentration in an industrial park over time;
FIG. 3 is a schematic illustration of the accounting results of NO2 concentration in an industrial park over time;
FIG. 4 is a schematic illustration of accounting results of SO2 concentration in an industrial park over time;
FIG. 5 is a schematic illustration of accounting results of VOC concentration over time in an industrial park;
FIG. 6 is a schematic illustration of accounting results of PM10 concentration in an industrial park over time;
FIG. 7 is a schematic diagram of the calculation result of the concentration field distribution of pollutants in a certain range at a certain moment in an industrial park, wherein the units are as follows: ug.m -3;
FIG. 8 is a schematic diagram of the results of the analysis and solution of the data from the observation station to obtain the time-series change of the height of the boundary layer of the atmosphere.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and are not intended to limit the scope of the present invention.
In embodiment 1, as shown in fig. 1, the present embodiment introduces a method for calculating the capacity of an atmospheric environment, and the present invention firstly does not leak out the pollutants that are unorganized and are stolen to the atmosphere because the inversion is not based on the flow monitoring result of the discharge port, but according to the concentration of the pollutants in the air of the park; secondly, the method is based on the street-valley scale industrial park, but not based on the mesoscale principle; finally, the technical route avoids a great amount of calculation for solving an atmospheric diffusion equation (because the concentration of each place in the space is not required to be precisely known in practice, only the concentration of the site and the periphery is required to be known and then the concentration space distribution is judged, and the main characteristics are grasped), so that the calculation speed is greatly improved, and the level from real-time monitoring to real-time inversion can be reached.
The specific process of the method of the invention comprises the following steps:
acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target industrial park, and meteorological monitoring data of the near-ground station and the high-altitude station;
Performing self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain a key parameter b of the meteorological monitoring data;
Calculating to obtain the rough height Z 0 of the underlying surface of the target area according to the horizontal coordinate, the height information and the key parameter b of the meteorological monitoring data;
Calculating wind speed fields at different heights according to the rough height Z 0 of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to obtain a near-ground rational Charson coefficient R ic;
And solving the boundary layer height pbl of the target industrial park according to the near-surface Michelson coefficient R ic, and the result is shown in figure 8.
And acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating to obtain the atmospheric pollution environmental capacity E cap in the target industrial park by combining the horizontal wind speeds of the multiple points in the target industrial park and the boundary layer height pbl of the target industrial park.
According to the key parameters b of the meteorological monitoring data of the plurality of monitoring sites in the industrial park, the rough height Z 0 of the regional in the park is solved, and the equivalent relation is as follows:
(1);
Wherein Z 0 represents the rough height of the underlying surface of the target area, Z g represents the height of the near-ground station, Z h represents the height of the high-altitude station, and b is a key parameter, and the rough height is obtained through fitting a large amount of data.
According to the wind field data, the potential temperature deficiency temperature and other data in the industrial park measured by the ground surface stations and the high-altitude stations, a near-ground Consumer coefficient R ic is obtained, and the coefficient is an important physical quantity in an atmospheric physics theory and is used for representing the vertical stability of the atmosphere or the ratio of unstable energy of the atmosphere in a specific vertical layer junction;
(2);
wherein R ic represents near-ground Michelson, R 0 represents surface Raschel number, g represents gravitational acceleration, B represents a key parameter, θ v0 represents a virtual temperature (the virtual temperature is a known quantity, the virtual temperature is a temperature after the latent heat of vapor in an air mass is completely released, obtained by an observation value, and is an input quantity here), U z represents a wind speed of a high-altitude station, U g represents a wind speed of a near-ground station, Z g represents a near-ground station height, and Z h represents a high-altitude station height;
(3);
Where f 0 denotes the coriolis parameter.
(4.1);
(4.2);
(4.3);
Wherein: to obtain the vertical wind speed difference based on the observation of the park high-altitude station and the ground station,/> To generate the vertical height difference of the wind speed difference in the direction, namely the difference between the monitoring heights of the high altitude station and the ground station,/>Wind speed (m/s) at the top of the atmospheric boundary layer,/>The wind speed (m/s) of the ground monitoring station is Z is the height from the ground of the vertical direction at a specific position in space when solving for the atmospheric motion,/>Is the ground height (m) of the ground monitoring station,/>Is the gradient of the virtual potential temperature in the vertical direction,/>Is the difference between the virtual temperature in the vertical direction observed by the ground station and the high-altitude station.
In the formula (4.2) and the formula (4.3)、/>For calculating coefficients, iterative generation using the formula (4.1) equation,/>、/>Is an important intermediate variable for iterative calculation of the vertical change rate of wind speed after the formula is deformed.
Using the obtained Richson numbersSolving the height of the atmosphere boundary layer。
According to the obtained regional wind speed of the monitoring station and the average height of the atmospheric boundary layer, the two-dimensional spatial distribution of the atmospheric environment capacity of the industrial park and the total amount of the two-dimensional spatial distribution in the park range can be approximately obtained;
(5);
In the method, in the process of the invention, Representing regional atmospheric pollution environmental capacity on different horizontal coordinate points in a target industrial park area, wherein the integral of the regional atmospheric environmental capacity in the park on the horizontal space is the total atmospheric environmental capacity in the target park area,/>Is the real-time wind speed at different horizontal coordinates within the target area. The atmospheric pollutant concentration field and the atmospheric environment capacity are combined to calculate the atmospheric pollutant discharge flux of the industrial park, for different types of atmospheric pollutants, the respective discharge intensities can be obtained according to the atmospheric environment capacity and the corresponding monitoring concentration respectively, as shown in fig. 2, the calculation result of the PM2.5 concentration in the corresponding target industrial park changing with time sequence is shown in fig. 3, the calculation result of the NO2 concentration in the corresponding target industrial park changing with time sequence is shown in fig. 4, the calculation result of the SO2 concentration in the corresponding target industrial park changing with time sequence is shown in fig. 5, the calculation result of the VOC concentration in the corresponding target industrial park changing with time sequence is shown in fig. 6, and the calculation result of the PM10 concentration in the corresponding target industrial park changing with time sequence is shown in fig. 6;
(6);
In the method, in the process of the invention, For the atmospheric pollution discharge flux of the industrial park, q represents the atmospheric pollutant concentration field distribution data observed by the industrial park monitoring station, and the real-time near-ground horizontal space distribution result is shown in fig. 7. /(I)For the calculated atmospheric environmental capacity in the industrial park.
Embodiment 2, which is based on the same inventive concept as embodiment 1, presents an atmospheric environment capacity calculation device including:
The acquisition module is used for acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target area and weather monitoring data of the near-ground station and the high-altitude station;
The fitting module is used for carrying out self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain key parameters of the meteorological monitoring data;
The calculation module is used for calculating and obtaining the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the meteorological monitoring data; calculating wind speed fields on a plurality of groups of different heights according to the rough height of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to solve near-ground rational Charles coefficients; solving the boundary layer height of the target area according to the near-ground Consumer coefficient; and acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating by combining the horizontal wind speeds of the multiple points in the target area and the boundary layer height of the target area to obtain the atmospheric pollution environment capacity in the target industrial park.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention 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 invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. 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.
The foregoing is merely a preferred embodiment of the present invention, and it should be noted that modifications and variations could be made by those skilled in the art without departing from the technical principles of the present invention, and such modifications and variations should also be regarded as being within the scope of the invention.
Claims (6)
1. The atmospheric environment capacity calculation method is characterized by comprising the following steps:
Acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target area, and meteorological monitoring data of the near-ground station and the high-altitude station;
Performing self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain key parameters of the meteorological monitoring data;
Calculating to obtain the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the meteorological monitoring data;
Calculating wind speed fields on a plurality of groups of different heights according to the rough height of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to solve near-ground rational Charles coefficients;
Solving the boundary layer height of the target area according to the near-ground Consumer coefficient;
and acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating by combining the horizontal wind speeds of the multiple points in the target area and the boundary layer height of the target area to obtain the atmospheric pollution environment capacity in the target area.
2. The atmospheric environment capacity calculation method according to claim 1, wherein the calculation formula of the underlying surface roughness height of the target area is:
(1);
Where Z 0 denotes the target area underlying roughness height, Z g denotes the near ground station height, Z h denotes the high altitude station height, B represents a key parameter.
3. The atmospheric environment capacity calculation method according to claim 1, wherein the calculation formula of the near-ground richardson coefficient is:
(2);
Wherein R ic represents near-ground Concharrson, R 0 represents surface Raspy, g represents gravitational acceleration, b represents a key parameter, θ v0 represents virtual potential temperature, U z represents wind speed at a high-altitude station site, U g represents wind speed near-ground site, Z g represents near-ground site height, and Z h represents high-altitude station site height;
(3);
Where f 0 denotes the coriolis parameter.
4. The atmospheric environment capacity calculation method according to claim 1, wherein the calculation formula of the boundary layer height of the target area is:
(4);
Wherein: Δu is the difference in wind speed in the vertical direction based on the observation of the park high-altitude station and the ground station; Δz is the difference between the monitored heights of the overhead and ground stations, U z is the wind speed at the top of the atmospheric boundary layer, U g is the wind speed at the ground monitoring station, pbl is the ground height at the top of the atmospheric boundary layer, b 1 is the gradient of the virtual temperature in the vertical direction, and Δθ v is the difference between the virtual temperatures in the vertical direction as observed by the ground and overhead stations.
5. The atmospheric environment capacity calculation method according to claim 1, wherein the calculation formula of the atmospheric pollution environment capacity in the target area is:
(5);
Where E cap represents the atmospheric pollution environmental capacity within the target area, v i,j represents the real-time wind speed at different horizontal coordinates (i, j) within the target area, and pbl is the boundary layer height of the target area.
6. An atmospheric environment capacity calculation device, comprising:
The acquisition module is used for acquiring horizontal coordinates and height information of a near-ground station and a high-altitude station of a target area and weather monitoring data of the near-ground station and the high-altitude station;
The fitting module is used for carrying out self-adaptive fitting according to the meteorological monitoring data of the near-ground station and the high-altitude station to obtain key parameters of the meteorological monitoring data;
The calculation module is used for calculating and obtaining the rough height of the underlying surface of the target area according to the horizontal coordinates, the height information and the key parameters of the meteorological monitoring data; calculating wind speed fields on a plurality of groups of different heights according to the rough height of the underlying surface of the target area and the meteorological monitoring data of the near-ground station and the high-altitude station to solve near-ground rational Charles coefficients; solving the boundary layer height of the target area according to the near-ground Consumer coefficient; and acquiring the horizontal wind speeds of the multiple points in the target area according to the meteorological monitoring data of the near-ground station and the high-altitude station, and calculating by combining the horizontal wind speeds of the multiple points in the target area and the boundary layer height of the target area to obtain the atmospheric pollution environment capacity in the target area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410465170.0A CN118069958B (en) | 2024-04-18 | 2024-04-18 | Atmospheric environment capacity calculation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410465170.0A CN118069958B (en) | 2024-04-18 | 2024-04-18 | Atmospheric environment capacity calculation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN118069958A true CN118069958A (en) | 2024-05-24 |
CN118069958B CN118069958B (en) | 2024-07-19 |
Family
ID=91099467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410465170.0A Active CN118069958B (en) | 2024-04-18 | 2024-04-18 | Atmospheric environment capacity calculation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN118069958B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140324351A1 (en) * | 2008-11-13 | 2014-10-30 | Saint Louis University | Apparatus and method for providing environmental predictive indicators to emergency response managers |
CN105302980A (en) * | 2015-11-10 | 2016-02-03 | 中国科学院遥感与数字地球研究所 | Urban aerodynamic roughness inversion method based on SAR data |
CN113326624A (en) * | 2021-06-07 | 2021-08-31 | 国网河南省电力公司电力科学研究院 | Method and system for predicting height of atmospheric boundary layer in desert area |
CN115146564A (en) * | 2022-08-02 | 2022-10-04 | 上海市气候中心(上海区域气候中心) | Urban ground wind speed refined simulation method based on vertical hierarchical downscaling technology |
CN115545565A (en) * | 2022-11-24 | 2022-12-30 | 江苏省生态环境大数据有限公司 | Method and system for managing and controlling total amount of pollution discharged from park based on atmospheric environment quality |
-
2024
- 2024-04-18 CN CN202410465170.0A patent/CN118069958B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140324351A1 (en) * | 2008-11-13 | 2014-10-30 | Saint Louis University | Apparatus and method for providing environmental predictive indicators to emergency response managers |
CN105302980A (en) * | 2015-11-10 | 2016-02-03 | 中国科学院遥感与数字地球研究所 | Urban aerodynamic roughness inversion method based on SAR data |
CN113326624A (en) * | 2021-06-07 | 2021-08-31 | 国网河南省电力公司电力科学研究院 | Method and system for predicting height of atmospheric boundary layer in desert area |
CN115146564A (en) * | 2022-08-02 | 2022-10-04 | 上海市气候中心(上海区域气候中心) | Urban ground wind speed refined simulation method based on vertical hierarchical downscaling technology |
CN115545565A (en) * | 2022-11-24 | 2022-12-30 | 江苏省生态环境大数据有限公司 | Method and system for managing and controlling total amount of pollution discharged from park based on atmospheric environment quality |
Non-Patent Citations (1)
Title |
---|
吴蓉等: "1961-2010年安徽省大气环境容量系数变化特征分析", 气候变化研究进展, vol. 13, no. 6, 1 November 2017 (2017-11-01), pages 545 * |
Also Published As
Publication number | Publication date |
---|---|
CN118069958B (en) | 2024-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111537023B (en) | Diffusion simulation and source tracing method for atmospheric pollutants in industrial park | |
CN108648127B (en) | Method for locking urban air pollution hot spot area | |
CN110990905A (en) | Auxiliary design system for historical city protection development cooperative control scheme | |
CN107436343A (en) | It is a kind of to simulate the method for calculating sensitizing range pollutant concentration | |
CN113624921B (en) | Multi-mode comprehensive pollution tracing method | |
CN113436045A (en) | Atmospheric pollution multi-scale fusion traceability method and device and readable storage medium | |
CN104834017A (en) | Positioning method for hazardous gas leakage diffusion accident source | |
CN104865353A (en) | Atmospheric pollution data acquisition method for industrial park based on unmanned aerial vehicle | |
CN107145668A (en) | Evaluation method and device of the power plant stack gas pollutant emission to regional atmospheric haze | |
CN107229834A (en) | A kind of complicated landform emergency response air pollution DIFFUSION PREDICTION method | |
CN115420854B (en) | Atmospheric pollutant tracing method based on forward and backward model combination | |
CN104897853A (en) | Thermal power plant pollutant discharging monitoring display method based on tower type diffusion model | |
CN115688636B (en) | Atmospheric environmental pollution remediation simulation method and system based on computational fluid dynamics | |
CN107704641A (en) | Fine simulation of wind method based on outdoor scene vegetation spatial distribution roughness | |
CN114155129A (en) | Atmospheric environment traceability method and system based on industrial park | |
Han et al. | Large eddy simulation for atmospheric boundary layer flow over flat and complex terrains | |
Ruiz et al. | Experimental study of drift deposition from mechanical draft cooling towers in urban environments | |
CN108414692A (en) | A kind of the monitoring point setting method and its relevant device of pollution sources | |
CN112132450A (en) | Method for positioning gaseous pollutants | |
CN118069958B (en) | Atmospheric environment capacity calculation method and device | |
CN111898296B (en) | Multi-scale simulation method and system for nuclear material atmospheric diffusion and sedimentation | |
CN113295589B (en) | Raise dust monitoring method, device and system | |
CN112257354B (en) | Reverse positioning method for air pollution source under dynamic wind field condition | |
CN114722635A (en) | Layout method of atmospheric pollutant monitoring sensors in chemical industrial park | |
CN114441397A (en) | PM2.5 pollution source positioning method based on Gaussian plume and centroid positioning |
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 | ||
GR01 | Patent grant |