CN113554305B - Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters - Google Patents
Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters Download PDFInfo
- Publication number
- CN113554305B CN113554305B CN202110828556.XA CN202110828556A CN113554305B CN 113554305 B CN113554305 B CN 113554305B CN 202110828556 A CN202110828556 A CN 202110828556A CN 113554305 B CN113554305 B CN 113554305B
- Authority
- CN
- China
- Prior art keywords
- straw burning
- concentration
- grid
- concentrated
- influence
- 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.)
- Active
Links
- 239000010902 straw Substances 0.000 title claims abstract description 181
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000010419 fine particle Substances 0.000 claims abstract description 41
- 238000011156 evaluation Methods 0.000 claims abstract description 39
- 238000012544 monitoring process Methods 0.000 claims abstract description 34
- 239000013618 particulate matter Substances 0.000 claims abstract description 22
- 230000008859 change Effects 0.000 claims abstract description 12
- 238000011160 research Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 17
- 210000005104 human peripheral blood lymphocyte Anatomy 0.000 claims description 12
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 230000003203 everyday effect Effects 0.000 description 5
- 230000002354 daily effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 239000000443 aerosol Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Spectroscopy & Molecular Physics (AREA)
Abstract
The invention discloses a method for evaluating the influence of straw burning on the concentration of atmospheric fine particles, which comprises the following steps: s1, counting the number of straw burning fire points of each day in an area according to a satellite remote sensing monitoring result, and selecting and determining a concentrated straw burning time period; s2, according to the satellite remote sensing monitoring result, carrying out inversion to obtain a regional atmosphere fine particulate matter concentration distribution result; s3, carrying out grid division on the research area according to meteorological conditions; s4, calculating a background value of the concentration of the fine particles grid by grid; and S5, quantitatively calculating the contribution percentage of the straw burning in each grid to the air quality influence. According to the method, the atmospheric fine particle concentration gridding evaluation method model is constructed by fusing the atmospheric pollution source and the atmospheric fine particle concentration spatial distribution, so that the evaluation method for the influence of straw burning on the atmospheric environment quality is shown, the straw burning intensity condition of the region can be reflected, and the influence of the straw burning on the change of the atmospheric fine particle concentration can be reflected.
Description
Technical Field
The invention relates to the technical field of atmospheric environment quality evaluation, in particular to a method for evaluating the influence of straw burning on the concentration of atmospheric fine particles.
Background
China is a big agricultural country, a large amount of crop straws are generated every year, and the concentration of atmospheric fine particles is easily and rapidly increased and the atmospheric environment quality is seriously deteriorated due to the fact that a large amount of primary fine particles and secondary aerosol formed by the primary fine particles are discharged from straws burnt in the open air. In recent years, departments such as agriculture and ecological environment in China adopt a series of measures to control the burning prohibition of straws and improve the comprehensive utilization rate, but because the straw burning is generally large in quantity and wide in range, the burning prohibition supervision difficulty is higher, and the accurate evaluation of the influence of the straw burning on the regional atmospheric environment quality is also difficult. Satellite remote sensing is used as a new technology, provides an important means for monitoring atmospheric pollution sources such as straw burning fire points and the like and evaluating the influence of the atmospheric pollution sources on the atmospheric environment quality, has the characteristics of macroscopicity, dynamic property, objectivity, accuracy and the like, has unique advantages in a space range compared with the traditional ground monitoring means in the aspect of information acquisition, and can continuously acquire the straw burning fire points in a large range of areas and the spatial and temporal changes of the concentration distribution of regional atmospheric fine particles, so that the influence degree of the regional straw burning fire points on the concentration of the atmospheric fine particles can be effectively reflected. At present, the mainstream straw burning fire point monitoring at home and abroad can realize the highest 1 time per 10 minutes, the highest spatial resolution can reach 375 meters, the highest resolution of the atmospheric fine particulate matter concentration can reach 1 kilometer, the monitoring frequency can reach once a day, and the requirements of evaluating monthly, quarterly and annual conditions of the influence of the regional straw burning fire point on the atmospheric fine particulate matter concentration can be met.
At present, research and analysis on the influence of straw burning on the concentration of regional atmosphere fine particles are mainly developed based on real-time measurement of ground monitoring sites or atmospheric chemical mode simulation technology, but accurate assessment is difficult to carry out, on one hand, the method is only based on ground limited site data mainly distributed in urban built-up areas, and on the other hand, straw burning mainly occurs in rural areas, and the method is not enough in the coverage degree of the areas; on the other hand, the atmospheric chemical mode is only based on the satellite remote sensing monitoring of the fire point in the part on the evaluation method, and dynamic change of a pollution source list and spatial distribution change of atmospheric fine particulate matter concentration are not considered, so that the traditional evaluation method has the limitation that the contribution of regional atmospheric fine particulate matter concentration cannot be accurately evaluated due to untimely source emission list and insufficient monitoring coverage.
Disclosure of Invention
The invention aims to provide an evaluation method for the influence of straw burning on the concentration of atmospheric fine particles, which constructs a gridding evaluation method model for the concentration of the atmospheric fine particles by fusing an atmospheric pollution source and the spatial distribution of the concentration of the atmospheric fine particles to accurately calculate the influence evaluation method for the straw burning on the quality of the atmospheric environment, can reflect the condition of the straw burning strength of a region and the influence of the straw burning on the change of the concentration of the atmospheric fine particles, and thus overcomes the limitation that the contribution of the concentration of the atmospheric fine particles in the region cannot be accurately evaluated due to untimely source emission list and insufficient monitoring coverage in the traditional evaluation method.
In order to achieve the purpose, the invention provides a method for evaluating the influence of straw burning on the concentration of atmospheric fine particulate matters, which is characterized by comprising the following steps of:
s1, counting the number of straw burning fire points of each day in an area according to a satellite remote sensing monitoring result, and selecting and determining a concentrated straw burning time period;
s2, according to the satellite remote sensing monitoring result, carrying out inversion to obtain a regional atmosphere fine particulate matter concentration distribution result, and carrying out synthetic calculation according to the straw burning time period;
s3, carrying out grid division on the research area according to meteorological conditions;
s4, calculating a background value of the concentration of the fine particles grid by grid;
and S5, quantitatively calculating the contribution percentage of the straw burning in each grid to the air quality influence.
Preferably, the step S1 further includes:
s1.1, the satellite remote sensing monitoring comprises orbit satellite remote sensing detection, multispectral monitoring data of wave bands of 4 micrometers of mid-infrared and 11 micrometers of thermal infrared of polar orbit satellite remote sensing are utilized, and daily thermal anomaly pixel information of an evaluation area is obtained according to a context algorithm;
s1.2 taking each thermal anomaly pixel coordinate as a center to acquire spatial resolution of thermal anomaly multispectral satellite remote sensing dataMaking a thermal anomaly buffer zone with the radius of times;
s1.3, superposing 2m high-resolution satellite images passing within 1 year on the daily thermal abnormal pixel buffer area, and judging the land utilization type within the range of the thermal abnormal pixel buffer area through visual interpretation;
s1.4, counting the number of straw burning fire points in each unit cultivated area of a region every day, and determining a concentrated straw burning period according to the straw burning fire points every day, wherein the concentrated straw burning period comprises a concentrated straw burning early stage, a concentrated straw burning period and a concentrated straw burning later stage;
preferably, the step S2 further includes:
s2.1, the satellite remote sensing monitoring comprises multispectral satellite remote sensing monitoring, and according to multispectral satellite remote sensing monitoring data, a geographical weighted regression method is used for constructing the following satellite remote sensing atmospheric fine particulate concentration remote sensing inversion method model:
lnPM2.5(ui,vi)=β0(ui,vi)+β1(ui,vi)lnAOD+β2(ui,vi)lnHPBL+
β3(ui,vi)ln(1-RH/100)
where β 0(ui, vi) is a regression coefficient of the constant term at the observation point (ui, vi), β 1(ui, vi) is a regression coefficient of the AOD at the observation point (ui, vi), β 2(ui, vi) is a regression coefficient of the HPBL at the observation point (ui, vi), β 3(ui, vi) is a regression coefficient of the RH at the observation point (ui, vi), and PM2.5(ui, vi) is a PM2.5 concentration at the observation point (ui, vi).
S2.2 remote sensing inversion of PM2.5 concentration in an evaluation area:
and obtaining a weight function matrix according to the geographical weighted regression model by adopting a weighted least square method, wherein the weight function matrix is as follows:
β(ui,vi)=[XTW(ui,vi)X]-1XTW(ui,vi)Y
wherein, β is a regression coefficient, W is a weight function matrix, Y is a PM2.5 concentration matrix, X is an input parameter matrix, and the input parameters include: constant terms, AOD, HPBL, and RH;
acquiring the optimal bandwidth of each ground station according to the weight function matrix by adopting a cross verification method, and acquiring the weight function matrix of each input parameter according to the optimal bandwidth;
carrying out spatial matching on the weight function matrix of each input parameter and the geographical position of the ground station corresponding to the input parameter to obtain a spatial function weight matrix;
obtaining regression coefficients of a constant term, AOD, HPBL and RH respectively according to the spatial function weight matrix by a kriging spatial interpolation method;
and according to the AOD, HPBL and RH corresponding to each pixel in the satellite remote sensing data and the regression coefficient, combining the geographical weighted regression model to obtain the PM2.5 concentration of each pixel in the evaluation area.
S2.3 evaluation region PM2.5 concentration synthesis calculation:
according to the PM2.5 concentrations of the regions in the early stage of concentrated straw burning, the period of concentrated straw burning and the later stage of concentrated straw burning, the average PM2.5 concentration of the regions in three periods is synthesized and calculated.
Preferably, the step S3 further includes:
s3.1, the meteorological conditions comprise regional average wind speed, the regional average wind speed is counted, and the average influence radius of straw burning on the air quality is calculated according to the influence range of the wind field in the concentrated straw burning period for 1 hour, wherein the average influence radius of the straw burning on the air quality is as follows:
wherein R is the influence radius (in km), U10For evaluating the wind speed (in m/s) of a 10m wind field latitudinal component on the ground of the region, V10The method comprises the steps that the wind speed of a 10m wind field warp component on the ground of an evaluation area (unit is m/s), M, N is the number of rows and columns of a simulation wind field grid of the evaluation area, and T is the number of simulation days;
and S3.2, dividing the evaluation area into grids of R × R size at equal intervals by taking the influence radius R as the grid size.
4. The method for evaluating the influence of straw burning based on satellite remote sensing on the concentration of the atmospheric fine particulate matters in claim 1, wherein the step S4 further comprises the following steps:
and S4.1, overlapping and evaluating the straw burning fire points of the areas according to the grids defined in the step S3, and counting the number of the straw burning fire points of each area.
S4.2, according to the regional atmosphere fine particle concentration calculated in the step S2 in the early stage of concentrated straw burning, the period of concentrated straw burning and the later stage of concentrated straw burning, superposing the grids defined in the step 3, and respectively counting the average value of the atmosphere fine particle concentration of each grid in different periods;
s4.3, calculating the background value of the concentration of the atmospheric fine particulate matters affecting the grid during the concentrated straw burning, wherein the calculation formula of the background value of the concentration of the atmospheric fine particulate matters is as follows:
wherein the content of the first and second substances,is the atmospheric fine particle concentration background value of the kth grid during the concentrated burning of the straws,andconcentration value, eta, of atmospheric fine particulate matters of kth grid in the early stage and the later stage of concentrated straw burning1And η2The weight factors of the average influence of the early concentrated straw burning stage and the later concentrated straw burning stage on the grid background concentration during the straw burning are determined according to the concentration proportion of the fine particulate matters in the atmosphere of the non-straw burning influence grid, and the calculation formula is as follows:
wherein the content of the first and second substances,in order to influence the concentration of the fine particulate matters in the grid atmosphere during the concentrated straw burning period,andthe concentrations of the atmospheric fine particles influencing the grids are respectively the non-straw burning in the early stage of straw burning and the later stage of straw burning.
Preferably, the step S5 further includes:
calculating the background value of the concentration of the fine particles in the atmosphere of the grid influenced by the straw burning during the concentrated straw burning obtained in the step 4, and calculating the percentage of contribution of the straw burning of each grid to the concentration change of the fine particles in the atmosphere in the area, wherein the calculation formula is as follows:
wherein, CkThe k grid straw burning contributes percentage to the concentration of atmospheric fine particulate matter during the concentrated straw burning,andthe concentration and the background concentration of the fine particulate matters in the grid atmosphere influenced by straw burning during the concentrated straw burning are monitored in a remote sensing mode.
Compared with the prior art, the invention has the beneficial effects that:
a method for evaluating the concentration of atmospheric fine particulate matters by regional straw burning based on satellite remote sensing is constructed, a gridding evaluation method model of the concentration of the atmospheric fine particulate matters is constructed by fusing an atmospheric pollution source and the spatial distribution of the concentration of the atmospheric fine particulate matters to accurately calculate the method for evaluating the influence of the straw burning on the atmospheric environment quality, the method not only can reflect the condition of the regional straw burning strength, but also can reflect the influence of the straw burning on the change of the concentration of the atmospheric fine particulate matters, so that the limitation that the contribution of the regional atmospheric fine particulate matters cannot be accurately evaluated due to untimely source emission list and insufficient monitoring coverage in the traditional evaluation method is overcome, and the method can objectively and quantitatively characterize the influence of the regional straw burning on the atmospheric environment quality.
Drawings
FIG. 1 is a flow chart of a method for evaluating the influence of straw burning on the concentration of atmospheric fine particulate matters.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, a method for evaluating the influence of straw burning on the concentration of atmospheric fine particulate matters is characterized by comprising the following steps:
s1, counting the number of straw burning fire points of each day in an area according to a satellite remote sensing monitoring result, and selecting and determining a concentrated straw burning time period;
s2, according to the satellite remote sensing monitoring result, carrying out inversion to obtain a regional atmosphere fine particulate matter concentration distribution result, and carrying out synthetic calculation according to the straw burning time period;
s3, carrying out grid division on the research area according to meteorological conditions;
s4, calculating a background value of the concentration of the fine particles grid by grid;
and S5, quantitatively calculating the contribution percentage of the straw burning in each grid to the air quality influence.
In this embodiment: the step S1 further includes:
s1.1, acquiring daily thermal abnormal pixel information of an evaluation area by utilizing multispectral monitoring data of mid-infrared (near 4 mu m) and thermal infrared (near 11 mu m) wave bands remotely sensed by polar orbit satellites according to a context algorithm;
s1.2 taking each thermal anomaly pixel coordinate as a center to acquire spatial resolution of thermal anomaly multispectral satellite remote sensing dataMaking a thermal anomaly buffer zone with the radius of times;
s1.3, superposing high-resolution satellite images which are better than 2m and transit within nearly 1 year on a daily thermal abnormal buffer area, judging the main land utilization type within the range of the thermal abnormal pixel buffer area through visual interpretation, and if the farmland or cultivated land area in the thermal abnormal buffer area accounts for more than 80% of the area proportion of the buffer area, judging that the thermal abnormal point is a straw burning fire point; otherwise, the fire is regarded as other fire points;
s1.4, counting the number of straw burning fire points in each unit cultivated land area of a region every day, determining a concentrated straw burning period according to the straw burning fire points every day, and determining that the day is the initial day of the concentrated straw burning period when the average number of the fire points in each hundred hectare cultivated land area per day of the region continuously five days from a certain day exceeds 1; when the number of fire points of the area per hundred hectare of cultivated land per day in the area of five continuous days from a certain day is lower than 0.1, determining that the day is the end day of the straw burning period, setting T days before the initial day of the straw centralized burning period as the early stage of the straw centralized burning according to the time length (T days) of the straw centralized burning period, and setting T days after the end day of the straw centralized burning period as the later stage of the straw centralized burning.
In this embodiment: the step S2 further includes:
s2.1, constructing the following remote sensing inversion method model of the concentration of the satellite remote sensing fine particles by using a geographical weighted regression method according to multispectral satellite remote sensing monitoring data:
lnPM2.5(ui,vi)=β0(ui,vi)+β1(ui,vi)lnAOD+β2(ui,vi)lnHPBL+
β3(ui,vi)ln(1-RH/100)
where β 0(ui, vi) is a regression coefficient of the constant term at the observation point (ui, vi), β 1(ui, vi) is a regression coefficient of the AOD at the observation point (ui, vi), β 2(ui, vi) is a regression coefficient of the HPBL at the observation point (ui, vi), β 3(ui, vi) is a regression coefficient of the RH at the observation point (ui, vi), and PM2.5(ui, vi) is a PM2.5 concentration at the observation point (ui, vi).
S2.2 remote sensing inversion of PM2.5 concentration in an evaluation area:
and obtaining a weight function matrix according to the geographical weighted regression model by adopting a weighted least square method, wherein the weight function matrix is as follows:
β(ui,vi)=[XTW(ui,vi)X]-1XTW(ui,vi)Y
wherein, β is a regression coefficient, W is a weight function matrix, Y is a PM2.5 concentration matrix, X is an input parameter matrix, and the input parameters include: constant terms, AOD, HPBL, and RH;
acquiring the optimal bandwidth of each ground station according to the weight function matrix by adopting a cross verification method, and acquiring the weight function matrix of each input parameter according to the optimal bandwidth;
carrying out spatial matching on the weight function matrix of each input parameter and the geographical position of the ground station corresponding to the input parameter to obtain a spatial function weight matrix;
obtaining regression coefficients of a constant term, AOD, HPBL and RH respectively according to the spatial function weight matrix by a kriging spatial interpolation method;
and according to the AOD, HPBL and RH corresponding to each pixel in the satellite remote sensing data and the regression coefficient, combining the geographical weighted regression model to obtain the PM2.5 concentration of each pixel in the evaluation area.
S2.3 evaluation region PM2.5 concentration synthesis calculation:
according to the PM2.5 concentrations of the regions in the early stage of concentrated straw burning, the period of concentrated straw burning and the later stage of concentrated straw burning, the average PM2.5 concentration of the regions in three periods is synthesized and calculated. The specific synthetic calculation method comprises the following steps: and respectively carrying out statistics on the average value of the effective PM2.5 concentration monitoring results pixel by pixel in each time period according to the remote PM2.5 concentration monitoring results in the region every day in the three time periods.
In this embodiment: the step S3 further includes:
s3.1, counting and evaluating the average wind speed of the area, and calculating the average influence radius of the straw burning on the air quality according to the 1-hour influence range of the wind field in the concentrated straw burning period, wherein the average influence radius is as follows:
wherein R is the influence radius (in km), U10For evaluating the wind speed (in m/s) of a 10m wind field latitudinal component on the ground of the region, V10The method comprises the steps that the wind speed of a 10m wind field warp component on the ground of an evaluation area (unit is m/s), M, N is the number of rows and columns of a simulation wind field grid of the evaluation area, and T is the number of simulation days;
and S3.2, dividing the evaluation area into grids of R × R size at equal intervals by taking the influence radius R as the grid size. In this embodiment: the step S4 further includes:
s4.1, overlapping and evaluating the straw burning fire points of the areas according to the grids defined by the method, counting the number of the straw burning fire points of each area, and if the number of the straw burning fire points in the grids is more than or equal to 3, determining that the atmospheric environment in the grids is greatly influenced by straw burning and determining that the grids are influenced by straw burning; if the number of the straw burning fire points in the grid is less than 3, the influence of straw burning on the atmospheric environment in the grid is considered to be negligible, and the grid is judged to be a non-straw burning influence grid.
S4.2, calculating the concentration of the atmospheric fine particles in the areas at the early stage, the middle stage and the later stage of the concentrated straw burning according to the method, superposing the grids defined by the method, and respectively counting the average value of the concentration of the atmospheric fine particles in each grid at different periods;
s4.3, calculating the background value of the concentration of the atmospheric fine particulate matters affecting the grids during the concentrated straw burning. Calculating the background value of the straw burning influence grid (namely the concentration of the atmospheric fine particulate matters under the condition of assuming that no straw burning fire point occurs) according to the change of the non-straw burning influence grid at different periods of the concentrated straw burning, wherein the calculation method comprises the following steps
Wherein the content of the first and second substances,is the atmospheric fine particle concentration background value of the kth grid during the concentrated burning of the straws,andconcentration value, eta, of atmospheric fine particulate matters of kth grid in the early stage and the later stage of concentrated straw burning1And η2Respectively as the early stage and the later stage of the centralized burning of the straws in the evaluation areaThe grid background concentration average influence weight factor during the straw burning period is determined according to the concentration proportion of the non-straw burning influence grid atmospheric fine particulate matter, and the calculation formula is as follows:
wherein the content of the first and second substances,in order to influence the concentration of the fine particulate matters in the grid atmosphere during the concentrated straw burning period,andthe concentrations of the atmospheric fine particles influencing the grids are respectively the non-straw burning in the early stage of straw burning and the later stage of straw burning.
It can be understood that, the weighted average method adopted by the method is used for calculation, and other weighting coefficients can also be determined according to the situation through a weighting model.
In this embodiment: the step S5 further includes:
calculating the background value of the concentration of the fine particles in the atmosphere of the grid influenced by the straw burning during the concentrated straw burning period obtained by the method, and calculating the percentage contribution of the straw burning of each grid to the concentration change of the fine particles in the atmosphere in the area, wherein the calculation formula is as follows:
wherein, CkThe k grid straw burning contributes percentage to the concentration of atmospheric fine particulate matter during the concentrated straw burning,andthe concentration and the background concentration of the fine particulate matters in the grid atmosphere influenced by straw burning during the concentrated straw burning are monitored in a remote sensing mode.
The method for evaluating the influence of straw burning on the concentration of the atmospheric fine particulate matters based on satellite remote sensing is a comprehensive evaluation method, is a quantitative index for reflecting the space-time distribution characteristics of atmospheric pollution sources such as regional straw burning fire points and the like and the air quality change such as the concentration of the atmospheric fine particulate matters, and reflects the contribution of pollutants discharged by straw burning to the generation of the concentration of the atmospheric fine particulate matters. Therefore, the method overcomes the disadvantages of insufficient traditional monitoring data or incomplete updating of the pollution source list, and makes the contribution evaluation of the straw burning fire point to the concentration of the atmospheric fine particulate matters more comprehensive through satellite remote sensing.
It should be noted that the whole evaluation process can be calculated based on the straw burning fire point of the grid and the concentration value of the atmospheric fine particulate matter thereof, and the evaluation value of each grid is obtained, so that the influence distribution condition of the area is obtained; the average value of all grids can be counted according to the area size, and the comprehensive average influence contribution in the administrative district can be obtained.
In summary, the invention constructs a grid evaluation method model of the concentration of the atmospheric fine particulate matter by fusing the atmospheric pollution source and the spatial distribution of the concentration of the atmospheric fine particulate matter based on the wide-range high-resolution pixel-level straw burning fire point and the spatial and temporal dynamic change of the concentration of the atmospheric fine particulate matter obtained by the satellite monitoring means, so as to accurately calculate the evaluation method of the influence of the straw burning on the atmospheric environment quality, and not only can reflect the condition of the straw burning strength of the region, but also can reflect the contribution of the straw burning on the change of the concentration of the atmospheric fine particulate matter, thereby overcoming the limitation that the contribution of the concentration of the atmospheric fine particulate matter of the region cannot be accurately evaluated due to untimely source emission list and insufficient monitoring coverage in the traditional evaluation method, and objectively quantitatively representing the influence of the straw burning of the region on the atmospheric environment quality. Therefore, the method model for evaluating the influence of straw burning on the concentration of the atmospheric fine particles provides a new effective technical means for comprehensively evaluating the influence level of the regional atmospheric pollution source on the atmospheric environment quality.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
While the invention has been described above with reference to an embodiment, various modifications may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In particular, the various features of the embodiments disclosed herein may be used in any combination, provided that there is no structural conflict, and the combinations are not exhaustively described in this specification merely for the sake of brevity and conservation of resources. Therefore, it is intended that the invention not be limited to the particular embodiments disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.
Claims (5)
1. A method for evaluating the influence of straw burning on the concentration of atmospheric fine particles is characterized by comprising the following steps:
s1, counting the number of straw burning fire points of each day in an area according to a satellite remote sensing monitoring result, and selecting and determining a concentrated straw burning time period;
s2, according to the satellite remote sensing monitoring result, carrying out inversion to obtain a regional atmosphere fine particulate matter concentration distribution result, and carrying out synthetic calculation according to the straw burning time period;
s3, carrying out grid division on the research area according to meteorological conditions;
s4, calculating a background value of the concentration of the fine particles grid by grid;
s5, quantitatively calculating the contribution percentage of straw burning in each grid on the air quality;
the step S4 further includes:
s4.1, overlapping and evaluating the straw burning fire points of the areas according to the grids defined in the step S3, and counting the number of the straw burning fire points of each area;
s4.2, according to the regional atmosphere fine particle concentration calculated in the step S2 in the early stage of concentrated straw burning, the period of concentrated straw burning and the later stage of concentrated straw burning, superposing the grids defined in the step 3, and respectively counting the average value of the atmosphere fine particle concentration of each grid in different periods;
s4.3, calculating the background value of the concentration of the atmospheric fine particulate matters affecting the grid during the concentrated straw burning, wherein the calculation formula of the background value of the concentration of the atmospheric fine particulate matters is as follows:
wherein the content of the first and second substances,is the atmospheric fine particle concentration background value of the kth grid during the concentrated burning of the straws,andconcentration value, eta, of atmospheric fine particulate matters of kth grid in the early stage and the later stage of concentrated straw burning1And η2The weight factors of the average influence of the early concentrated straw burning stage and the later concentrated straw burning stage on the grid background concentration during the straw burning are determined according to the concentration proportion of the fine particulate matters in the atmosphere of the non-straw burning influence grid, and the calculation formula is as follows:
wherein the content of the first and second substances,for non-straw during concentrated burning of strawThe burning affects the concentration of the fine particulate matters in the grid atmosphere,andthe concentrations of the atmospheric fine particles influencing the grids are respectively the non-straw burning in the early stage of straw burning and the later stage of straw burning.
2. The method for evaluating the influence of straw incineration on the concentration of atmospheric fine particulate matter according to claim 1, wherein the step S1 further comprises:
s1.1, the satellite remote sensing monitoring comprises orbit satellite remote sensing detection, multispectral monitoring data of wave bands of 4 micrometers of mid-infrared and 11 micrometers of thermal infrared of polar orbit satellite remote sensing are utilized, and daily thermal anomaly pixel information of an evaluation area is obtained according to a context algorithm;
s1.2 taking each thermal anomaly pixel coordinate as a center to acquire spatial resolution of thermal anomaly multispectral satellite remote sensing dataMaking a thermal anomaly buffer zone with the radius of times;
s1.3, superposing 2m high-resolution satellite images passing within 1 year on the daily thermal abnormal pixel buffer area, and judging the land utilization type within the range of the thermal abnormal pixel buffer area through visual interpretation;
s1.4, counting the number of the straw burning fire points of the unit cultivated area of each day in the area, and determining a concentrated straw burning period according to the straw burning fire points of each day, wherein the concentrated straw burning period comprises a concentrated straw burning early stage, a concentrated straw burning period and a concentrated straw burning later stage.
3. The method for evaluating the influence of straw incineration on the concentration of atmospheric fine particulate matter according to claim 1, wherein the step S2 further comprises:
s2.1, the satellite remote sensing monitoring comprises multispectral satellite remote sensing monitoring, and according to multispectral satellite remote sensing monitoring data, a geographical weighted regression method is used for constructing the following satellite remote sensing atmospheric fine particulate concentration remote sensing inversion method model:
lnPM2.5(ui,vi)=β0(ui,vi)+β1(ui,vi)lnAOD+β2(ui,vi)lnHPBL+β3(ui,vi)ln(1-RH/100)
wherein β 0(ui, vi) is a regression coefficient of the constant term at the observation point (ui, vi), β 1(ui, vi) is a regression coefficient of the AOD at the observation point (ui, vi), β 2(ui, vi) is a regression coefficient of the HPBL at the observation point (ui, vi), β 3(ui, vi) is a regression coefficient of the RH at the observation point (ui, vi), and PM2.5(ui, vi) is a PM2.5 concentration at the observation point (ui, vi);
s2.2 remote sensing inversion of PM2.5 concentration in an evaluation area:
and obtaining a weight function matrix according to the geographical weighted regression model by adopting a weighted least square method, wherein the weight function matrix is as follows:
β(ui,vi)=[XTW(ui,vi)X]-1XTW(ui,vi)Y
wherein, β is a regression coefficient, W is a weight function matrix, Y is a PM2.5 concentration matrix, X is an input parameter matrix, and the input parameters include: constant terms, AOD, HPBL, and RH;
acquiring the optimal bandwidth of each ground station according to the weight function matrix by adopting a cross verification method, and acquiring the weight function matrix of each input parameter according to the optimal bandwidth;
carrying out spatial matching on the weight function matrix of each input parameter and the geographical position of the ground station corresponding to the input parameter to obtain a spatial function weight matrix;
obtaining regression coefficients of a constant term, AOD, HPBL and RH respectively according to the spatial function weight matrix by a kriging spatial interpolation method;
according to the AOD, HPBL and RH corresponding to each pixel in the satellite remote sensing data and the regression coefficient, combining the geographical weighted regression model to obtain the PM2.5 concentration of each pixel in the evaluation area;
s2.3 evaluation region PM2.5 concentration synthesis calculation:
according to the PM2.5 concentrations of the regions in the early stage of concentrated straw burning, the period of concentrated straw burning and the later stage of concentrated straw burning, the average PM2.5 concentration of the regions in three periods is synthesized and calculated.
4. The method for evaluating the influence of straw incineration on the concentration of atmospheric fine particulate matter according to claim 1, wherein the step S3 further comprises:
s3.1, the meteorological conditions comprise regional average wind speed, the regional average wind speed is counted, and the average influence radius of straw burning on the air quality is calculated according to the influence range of the wind field in the concentrated straw burning period for 1 hour, wherein the average influence radius of the straw burning on the air quality is as follows:
wherein R is the influence radius (in km), U10For evaluating the wind speed (in m/s) of a 10m wind field latitudinal component on the ground of the region, V10The method comprises the steps that the wind speed of a 10m wind field warp component on the ground of an evaluation area (unit is m/s), M, N is the number of rows and columns of a simulation wind field grid of the evaluation area, and T is the number of simulation days;
and S3.2, dividing the evaluation area into grids of R × R size at equal intervals by taking the influence radius R as the grid size.
5. The method for evaluating the influence of straw incineration on the concentration of atmospheric fine particulate matter according to claim 1, wherein the step S5 further comprises:
calculating the background value of the concentration of the fine particles in the atmosphere of the grid influenced by the straw burning during the concentrated straw burning obtained in the step 4, and calculating the percentage of contribution of the straw burning of each grid to the concentration change of the fine particles in the atmosphere in the area, wherein the calculation formula is as follows:
wherein, CkThe k grid straw burning contributes percentage to the concentration of atmospheric fine particulate matter during the concentrated straw burning,andthe concentration and the background concentration of the fine particulate matters in the grid atmosphere influenced by straw burning during the concentrated straw burning are monitored in a remote sensing mode.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110828556.XA CN113554305B (en) | 2021-07-22 | 2021-07-22 | Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110828556.XA CN113554305B (en) | 2021-07-22 | 2021-07-22 | Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113554305A CN113554305A (en) | 2021-10-26 |
CN113554305B true CN113554305B (en) | 2022-04-22 |
Family
ID=78103992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110828556.XA Active CN113554305B (en) | 2021-07-22 | 2021-07-22 | Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113554305B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611971A (en) * | 2022-02-09 | 2023-08-18 | 南京浦蓝大气环境研究院有限公司 | County scale biomass combustion source normalization contribution potential quantification system |
CN116363047B (en) * | 2022-08-23 | 2024-01-30 | 生态环境部卫星环境应用中心 | Early warning method for pollution risk of straw incineration atmosphere |
CN116758342B (en) * | 2023-06-01 | 2023-12-15 | 中国地质科学院矿产资源研究所 | Atmospheric pollution grade assessment method and device based on rare earth mineral area |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700095A (en) * | 2015-03-30 | 2015-06-10 | 北京市环境保护监测中心 | Satellite remote sensing monitoring method and processing device for straw burning fire points |
WO2018023406A1 (en) * | 2016-08-02 | 2018-02-08 | 张阳 | Method for collecting data of smog early warning technology and early warning device |
CN110954482A (en) * | 2019-12-02 | 2020-04-03 | 生态环境部卫星环境应用中心 | Atmospheric pollution gridding monitoring method based on static satellite and polar orbit satellite |
CN112378828A (en) * | 2020-12-11 | 2021-02-19 | 中科三清科技有限公司 | Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data |
-
2021
- 2021-07-22 CN CN202110828556.XA patent/CN113554305B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700095A (en) * | 2015-03-30 | 2015-06-10 | 北京市环境保护监测中心 | Satellite remote sensing monitoring method and processing device for straw burning fire points |
WO2018023406A1 (en) * | 2016-08-02 | 2018-02-08 | 张阳 | Method for collecting data of smog early warning technology and early warning device |
CN110954482A (en) * | 2019-12-02 | 2020-04-03 | 生态环境部卫星环境应用中心 | Atmospheric pollution gridding monitoring method based on static satellite and polar orbit satellite |
CN112378828A (en) * | 2020-12-11 | 2021-02-19 | 中科三清科技有限公司 | Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data |
Non-Patent Citations (1)
Title |
---|
基于modis的PM2.5遥感监测及重点关注网格筛选方法研究;陈辉等;《中国环境科学学会科学技术年会论文集》;20180831;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113554305A (en) | 2021-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113554305B (en) | Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters | |
Moazami et al. | A comprehensive evaluation of GPM-IMERG V06 and MRMS with hourly ground-based precipitation observations across Canada | |
Blacutt et al. | Precipitation comparison for the CFSR, MERRA, TRMM3B42 and Combined Scheme datasets in Bolivia | |
Torma et al. | Validation of a high-resolution version of the regional climate model RegCM3 over the Carpathian basin | |
Sen et al. | Projecting climate change, drought conditions and crop productivity in Turkey | |
Walsh et al. | The exceptionally warm winter of 2015/16 in Alaska | |
Manders et al. | Testing the capability of the chemistry transport model LOTOS-EUROS to forecast PM10 levels in the Netherlands | |
CN102435586B (en) | Method and system for generating earth surface albedo product | |
Hong et al. | Inferring vertical variability and diurnal evolution of O3 formation sensitivity based on the vertical distribution of summertime HCHO and NO2 in Guangzhou, China | |
Xin et al. | Surface energy balance closure at ten sites over the Tibetan plateau | |
Rife et al. | Verification of temporal variations in mesoscale numerical wind forecasts | |
CN102779391A (en) | Drought early-warning method and drought early-warning system | |
Žabkar et al. | Evaluation of the high resolution WRF-Chem (v3. 4.1) air quality forecast and its comparison with statistical ozone predictions | |
Bodeker et al. | Aotearoa New Zealand climate change projections guidance: interpreting the latest IPCC WG1 report findings | |
Pierre et al. | Evaluation of catch efficiency transfer functions for unshielded and single-alter-shielded solid precipitation measurements | |
Tuygun et al. | Long-term temporal analysis of the columnar and surface aerosol relationship with planetary boundary layer height at a southern coastal site of Turkey | |
Yarragunta et al. | Validation of lower tropospheric carbon monoxide inferred from MOZART model simulation over India | |
Yassin et al. | Hydrologic-Land Surface Modelling of a Complex System under Precipitation Uncertainty: A Case Study of the Saskatchewan River Basin, Canada | |
Ding et al. | NOx emissions in India derived from OMI satellite observations | |
Zhao et al. | Correcting precipitation measurements made with Geonor T-200B weighing gauges near the August-one ice cap in the Qilian Mountains, Northwest China | |
Alerskans et al. | Optimizing spatial quality control for a dense network of meteorological stations | |
He et al. | Spatiotemporal change patterns and driving factors of land surface temperature in the Yunnan-Kweichow Plateau from 2000 to 2020 | |
Park et al. | Characteristics of atmospheric CO2 fluxes and the estimation of their potential sources around Boseong Standard Weather Observatory (BSWO) | |
Sicardi et al. | Assessment of Kalman filter bias-adjustment technique to improve the simulation of ground-level ozone over Spain | |
Torrez-Rodriguez et al. | Evaluation of temperature and precipitation from CORDEX-CORE South America and Eta-RCM regional climate simulations over the complex terrain of Subtropical Chile |
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 | ||
GR01 | Patent grant |