CN114496117A - Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment - Google Patents
Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment Download PDFInfo
- Publication number
- CN114496117A CN114496117A CN202111476909.0A CN202111476909A CN114496117A CN 114496117 A CN114496117 A CN 114496117A CN 202111476909 A CN202111476909 A CN 202111476909A CN 114496117 A CN114496117 A CN 114496117A
- Authority
- CN
- China
- Prior art keywords
- vertical column
- longitude
- tropospheric
- hcho
- spatial average
- 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
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 43
- WSFSSNUMVMOOMR-UHFFFAOYSA-N Formaldehyde Chemical compound O=C WSFSSNUMVMOOMR-UHFFFAOYSA-N 0.000 claims abstract description 72
- 239000005436 troposphere Substances 0.000 claims abstract description 55
- 239000012855 volatile organic compound Substances 0.000 claims abstract description 35
- 238000011160 research Methods 0.000 claims abstract description 6
- 231100001143 noxa Toxicity 0.000 claims abstract description 5
- 239000007789 gas Substances 0.000 claims description 76
- 238000012545 processing Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000006552 photochemical reaction Methods 0.000 claims description 6
- 230000001276 controlling effect Effects 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 4
- 239000011475 Accrington brick Substances 0.000 claims description 3
- WSFSSNUMVMOOMR-NJFSPNSNSA-N methanone Chemical compound O=[14CH2] WSFSSNUMVMOOMR-NJFSPNSNSA-N 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000002243 precursor Substances 0.000 abstract description 8
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 14
- 229910002091 carbon monoxide Inorganic materials 0.000 description 14
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 11
- 238000009826 distribution Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000003344 environmental pollutant Substances 0.000 description 3
- 231100000719 pollutant Toxicity 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003960 organic solvent Substances 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Crystallography & Structural Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses an ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment, wherein the method comprises the following steps: acquiring tropospheric vertical column concentration data of each target gas in a research time phase in a target area; acquiring the corresponding troposphere gas vertical column concentration of each target gas in each longitude and latitude grid point on the equal longitude and latitude grid; calculating the daily spatial average value of the tropospheric vertical column concentration of each target gas of all available picture elements; based on CO and O3And daily spatial average value data of the tropospheric vertical column concentration of the HCHO, and calculating the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source of all the available pixels; obtaining VOCs control region threshold and NOxA control zone threshold; VOCs-based control zone threshold and NOxControl zone thresholding into target zonesTropospheric ozone control types at each longitude and latitude grid point within the domain. The method realizes the definition of main ozone precursors in various places and further scientifically and effectively controls the concentration of ozone.
Description
Technical Field
The invention belongs to the technical field of air quality remote sensing monitoring, and particularly relates to an ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment.
Background
The presence of ozone in the troposphere is a biologically harmful pollutant and is one of the components of photochemical smog. Many human activities involving the rapid conversion of chemical energy, such as the start of internal combustion engines and the operation of copiers, produce ozone, a powerful oxidant that readily reacts with other chemicals to produce many toxic oxides that are harmful to human health. Therefore, tropospheric ozone concentrations are closely linked to production activities and the physical health of people.
Tropospheric ozone is mainly derived from photochemical reactions-ozone is generated when air mixed with various nitrogen oxides (NOx), carbon monoxide (CO) and volatile organic compounds (VOCs, such as formaldehyde-HCHO) is exposed to sunlight. Nitrogen oxides and volatile organics are therefore referred to as "ozone precursors". Automobile exhaust, industrial exhaust gases and chemical organic solvents are the major man-made sources of "ozone precursors". Although these sources are mostly concentrated in cities, some substances (such as nitrogen oxides) can spread by wind to sparsely populated areas hundreds of kilometers away, where they form a source of ozone.
It is now possible to measure vertical column (vertical atmospheric column) concentrations of various atmospheric pollutants in the troposphere using satellites. The satellite observation can realize the observation of a plurality of gas concentrations in a large range, and is convenient for the wide-area analysis of the relationship of the plurality of gas concentrations and the control of pollutants.
The existing control method for the troposphere ozone concentration mainly controls a high-concentration troposphere ozone point source, but is influenced by ozone photochemical reaction, and the control on the troposphere ozone point source cannot effectively reduce the troposphere ozone concentration in a target area, and even can obtain opposite effects.
Disclosure of Invention
The invention aims to provide an ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment, which are used for realizing large-range measurement of ozone in various places on two precursors of VOCs and NOxAnd the sensitivity of the troposphere is analyzed, and the main influence of the troposphere ozone concentration on which preconditions is obtained.
In order to achieve the aim, the invention provides an ozone control type identification method based on satellite hyperspectral remote sensing, which comprises the following steps:
acquiring troposphere vertical column concentration data of target gases in a target region in a research time phase based on atmospheric spectrum remote sensing data returned by a satellite, wherein the target gases comprise CO and O3HCHO and NO2；
Regulating the pixel of each target gas to the equal longitude and latitude grid of the target area to obtain the corresponding troposphere gas vertical column concentration of each target gas in each longitude and latitude grid point on the equal longitude and latitude grid;
obtaining available pixels in the target area, and calculating the daily spatial average value of the troposphere vertical column concentration of each target gas of all the available pixels;
based on CO and O3And daily spatial average value data of the troposphere vertical column concentration of HCHO, and calculating the daily spatial average value of the troposphere vertical column concentration of HCHO secondary sources of all available pixels, wherein the HCHO secondary sources are vertical column concentrations generated by photochemical reaction of troposphere formaldehyde;
for secondary source of HCHO, NO2Carrying out normalization processing on the daily spatial average value of the troposphere vertical column concentration to obtain a corresponding standard value;
based on O3The daily spatial average value of the tropospheric vertical column concentration, the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source, the standard value thereof, and NO2The daily spatial average value and the standard value of the tropospheric vertical column concentration to obtain the threshold value and NO of the VOCs control areaxA control zone threshold;
based on the VOCs control region threshold and the NOxAnd controlling zone threshold values, and dividing troposphere ozone control types of all longitude and latitude grid points in the target zone.
Optionally, the normalizing the pixels of each target gas to the equal longitude and latitude grid of the target area includes:
performing linear interpolation on the longitude and latitude of any two adjacent pixels on the original pixel plane of the satellite and the vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derived pixels and the original pixels into corresponding longitude and latitude grid points;
and averaging the vertical column concentration of each target gas pixel in each longitude and latitude grid point to obtain the vertical column concentration value of each target gas on each longitude and latitude grid point.
Optionally, the vertical column concentration of each target gas pixel in each longitude and latitude grid point is averaged, and the average is implemented by the following formula:
wherein, ViThe vertical column concentration of the troposphere gas, lon, corresponding to the target gas on the ith pixel needing to be processed on the original pixel plane of the satellite after linear interpolationi，latiIs the longitude and latitude, v, of the ith pixellin，latThe vertical column concentration of the troposphere corresponding to the target gas in the longitude and latitude grid points corresponding to the longitude and latitude grid points is obtained, lon and lat are the longitude and latitude of the ith pixel, and grid is the degree of difference between the longitude and latitude grid points.
Optionally, the acquiring the available image elements in the target region includes:
and removing the pixels with the cloud amount larger than 0.5 according to the cloud amount information of each pixel in the target area returned by the satellite, and taking the rest pixels as the available pixels.
Optionally, the calculating a daily spatial average of tropospheric vertical column concentrations for each target gas for all available picture elements is performed by:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O3Secondary sources of HCHO or NO2I represents the ith day in the set period of time, { region } isThe set of all the equal longitude and latitude grid points in the target area is described, N is the number of the longitude and latitude grid points in { region }, lon and lat are the longitude and latitude of a single longitude and latitude grid point, vx,i,lon,latIs the vertical column concentration of tropospheric gas at day i in a single latitude and longitude grid point for the target gas.
Optionally, the calculation of the daily spatial average value of troposphere vertical column concentrations of all available pixel elements HCHO secondary sources is performed by the following multivariate polynomial fitting operation:
VHCHO＝(1(VCO)1…(VCO)n)A+(1(VO3)1…(VO3)n)B+e
wherein Vp,HCHODaily spatial average of tropospheric vertical column concentrations, V, for secondary sources of HCHOHCHOIs the daily spatial average of the total vertical column concentration of tropospheric formaldehyde, VCO、VO3Are respectively CO and O3The daily spatial average of tropospheric vertical column concentrations of,
are each VCO、VO3Polynomial coefficient of (n)<5) And e represents the fitting residual.
Optionally, the pair of secondary sources of HCHO, NO2The normalization process is carried out on the daily spatial average value of the tropospheric vertical column concentration, and is realized by the following formula:
Vx,i,nor＝Vx,i/Vx,ref
Vx,refobtained from the following equation:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of a target gas x, x being HCSecondary source of HO or NO2，Vx,i,norIs a Vx,iNormalized standard value of (V)x,refDaily reference value for tropospheric vertical column concentration of target gas x, d is total days of the study time period.
Optionally, obtaining the VOCs control region threshold and the NOxThe method for controlling the zone threshold value comprises the following steps:
the daily spatial average of tropospheric vertical column concentrations of secondary sources of HCHO on the same day is compared with NO2The ratio of the daily spatial average of tropospheric vertical column concentrations of (a);
selecting as the first data set all data satisfying the flag being less than a first parametric reference value from data elements comprising a daily spatial average of tropospheric vertical pillar concentration of O3, a daily spatial average of tropospheric vertical pillar concentration of secondary HCHO sources and their standard values, and NO for each day of the target region over the study time period2The daily spatial average value and the standard value of the troposphere vertical column concentration;
linearly fitting the first data set to O3A first slope of the daily spatial average of tropospheric vertical column concentrations of (a) to a standard value of the daily spatial average of tropospheric vertical column concentrations of (b) a secondary source of HCHO, and O3Daily spatial average of tropospheric vertical column concentration and NO2Obtaining a first variable reference value corresponding to the first slope being equal to the second slope, and taking the first variable reference value at the moment as the threshold of the VOCs control area;
selecting all data satisfying that the flag is smaller than a second reference variable value from the data elements as a second data set;
linearly fitting the second data set to O3A third slope of the daily spatial average of tropospheric vertical column concentrations to a standard value of the daily spatial average of tropospheric vertical column concentrations of the secondary source of HCHO, and O3Daily spatial average of tropospheric vertical column concentration and NO2Is thick in the vertical column of the troposphereA fourth slope of the standard value of the daily spatial average value of the degree, a second variable reference value corresponding to the third slope being equal to the fourth slope is obtained, and the second variable reference value at that time is taken as the NOxA control zone threshold.
Optionally, the control region threshold based on the VOCs and the NOxA control zone threshold value, which is used for dividing troposphere ozone control types of all longitude and latitude grid points in the target zone and comprises the following steps:
calculating the troposphere vertical column concentration and NO of the HCHO secondary source in each longitude and latitude grid point in the target area2The tropospheric vertical column concentration ratio;
dividing the ozone control type of the longitude and latitude grid point with the ratio smaller than the threshold value of the VOCs control area into VOCs control types;
will be greater than the NOxOzone control type classification into NO for latitude and longitude grid points of control zone thresholdxA control type;
the ratio is larger than the threshold value of the VOCs control area and smaller than the NOxThe ozone control types of the longitude and latitude grid points of the control area threshold are divided into transition control types.
The present invention also proposes an electronic device, comprising:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the ozone control type identification method described above.
The invention has the beneficial effects that:
the method comprises the steps of firstly obtaining target gases CO and O in a target area in a research time stage3The tropospheric vertical column concentration data of HCHO and NO2, then the corresponding tropospheric gas vertical column concentration of each target gas in each longitude and latitude grid point on the equal longitude and latitude grid is obtained, and each target gas of all available pixels is calculatedThe daily spatial average value of the troposphere vertical column concentration is calculated, and then the daily spatial average value of the troposphere vertical column concentration of the HCHO secondary source of all the available pixel elements is calculated and is based on O3The daily spatial average value of the tropospheric vertical column concentration, the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source, the standard value thereof, and NO2The daily spatial average value and the standard value of the tropospheric vertical column concentration to obtain the threshold value and NO of the VOCs control areaxControl zone threshold, final control zone threshold based on VOCs and NOxControlling zone threshold, dividing troposphere ozone control types of all longitude and latitude grid points in the target zone, thereby realizing that the ozone suitable for satellites and capable of measuring the ozone of all regions in a large range can be used for two precursor VOCs and NOxThe sensitivity analysis method can definitely research the main influence of which precursor on the troposphere ozone concentration at each place in the area, and can control the ozone concentration of different areas in a targeted manner after determining the main precursor of the ozone concentration at different places.
The system of the present invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
Fig. 1 shows a step diagram of an ozone control type identification method based on satellite hyperspectral remote sensing according to the invention.
Fig. 2 shows the ozone control type distribution diagram for different stages of the same target area in an ozone control type identification method based on satellite hyperspectral remote sensing according to embodiment 1 of the invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Example 1
Fig. 1 shows a step diagram of an ozone control type identification method based on satellite hyperspectral remote sensing according to the invention.
As shown in fig. 1, a method for identifying ozone control type based on satellite hyperspectral remote sensing comprises the following steps:
step S101: acquiring troposphere vertical column concentration data of target gases in a target region in a research time phase based on atmospheric spectrum remote sensing data returned by a satellite, wherein the target gases comprise CO and O3HCHO and NO2；
Specifically, the concentration of the vertical column of the troposphere of various gases is obtained by obtaining various secondary products of the gases returned by the satellite.
Step S102: the pixel of each target gas is structured on an equal longitude and latitude grid of a target area, and the corresponding troposphere gas vertical column concentration of each target gas in each longitude and latitude grid point on the equal longitude and latitude grid is obtained;
in this step, the process of regulating the pixels of each target gas to the equal longitude and latitude grid of the target area includes:
performing linear interpolation on the longitude and latitude of any two adjacent pixels on the original pixel plane of the satellite and the vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derived pixels and the original pixels into corresponding longitude and latitude grid points;
and averaging the vertical column concentration of each target gas pixel in each longitude and latitude grid point to obtain the vertical column concentration value of each target gas on each longitude and latitude grid point.
Specifically, during data analysis, a plurality of gas data are needed, and the pixel position of each gas product returned by the satellite is different. In order to ensure the accuracy of the subsequent process, the pixels of various gases need to be structured on an equal longitude and latitude grid, and a sampling method is used for processing:
the original pixel plane of the satellite has two dimensions X ^ and Y ^ corresponding to the direction of a satellite instrument CCD and the advancing direction of the satellite, and the longitude and latitude and the gas concentration of any two adjacent pixels on the plane are subjected to linear interpolation to obtain a derivative pixel. And dividing the derived pixels and the original pixels to corresponding longitude and latitude grid points, and averaging the gas concentration of the pixels in each longitude and latitude grid point to obtain the gas concentration value on each longitude and latitude grid point.
In the process, the vertical column concentration of each target gas pixel in each longitude and latitude grid point is averaged, and the average concentration is realized by the following formula:
wherein, ViThe vertical column concentration of the troposphere gas, lon, corresponding to the target gas on the ith pixel needing to be processed on the original pixel plane of the satellite after linear interpolationi，latiIs the longitude and latitude, v, of the ith pixellon，latThe vertical column concentration of the troposphere corresponding to the target gas in the longitude and latitude grid points corresponding to the longitude and latitude grid points is obtained, lon and lat are the longitude and latitude of the ith pixel, and grid is the degree of difference between the longitude and latitude grid points.
Step S103: acquiring available pixels in a target area, and calculating the daily spatial average value of the troposphere vertical column concentration of each target gas of all the available pixels;
in particular, since the accuracy of trace gas inversion is severely affected by cloud layers, which may affect sensitivity analysis for ozone, it is necessary to remove unusable pixels affected by cloud amount. In this step, the method for acquiring the available pixels in the target area includes:
and according to the cloud amount information of each pixel in the target area returned by the satellite, removing the pixels with the cloud amount larger than 0.5, and taking the rest pixels as available pixels.
After obtaining the available pixels, selecting a target area, and adding O of all the available pixels in the target area3Secondary source of HCHO, NO2The tropospheric vertical column concentrations are averaged and divided by day into Vx,i。
In this step, the daily spatial average value of the tropospheric vertical column concentration of each target gas of all available pixels is calculated, which can be realized by the following formula:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O3, HCHO secondary source or NO2I represents the ith day in a set time period, the { region } is a set of all equal longitude and latitude grid points in the target area, N is the number of the longitude and latitude grid points in the { region }, lon and lat are the longitude and latitude of a single longitude and latitude grid point, and v is the longitude and latitude of the single longitude and latitude grid pointx,i,lon,latIs the tropospheric gas vertical column concentration of the target gas at day i in a single longitude and latitude point.
Step S104: based on CO and O3And daily spatial average value data of the troposphere vertical column concentration of HCHO, and calculating the daily spatial average value of the troposphere vertical column concentration of HCHO secondary sources of all available pixels, wherein the HCHO secondary sources are vertical column concentrations generated by photochemical reaction of troposphere formaldehyde;
in particular, atmospheric VOCs and NO due to photochemical reactionsxWill gradually generate O3Thus concentration of HCHO and O as an indicator of VOC concentration3There is a relationship in concentration. Meanwhile, HCHO can be generated in the incomplete combustion process of various biomasses, and CO can be generated in the process, so that the concentration of HCHO is related to the concentration of CO at the same time. The concentration of HCHO in the atmosphere can be represented by the following formula:
VHCHO-Vback＝F(VO3,VCO)
wherein, VHCHOIs the concentration of HCHO in the atmosphere, VbackIs the atmospheric HCHO background concentration value, VO3、VCOAre respectively O in the atmosphere3Concentration and CO concentration.
Can convert HCHO and O on certain places for multiple days3And CO concentration is used as a data element, a high-order polynomial is used for fitting, and various source proportions of formaldehyde are analyzed through multivariate polynomial fitting to establish HCHO and O3And CO, and determining that the troposphere formaldehyde belongs to the vertical column concentration generated by the photochemical reaction, namely the vertical column concentration of the HCHO secondary source.
The method is realized by the following multivariate polynomial fitting operation:
VHCHO＝(1(VCO)1…(VCO)n)A+(1(VO3)1…(VO3)n)B+e
wherein Vp,HCHODaily spatial average of tropospheric vertical column concentrations, V, for secondary sources of HCHOHCHOIs the daily spatial average of the total vertical column concentration of tropospheric formaldehyde, VCO、VO3Are respectively CO and O3The daily spatial average of tropospheric vertical column concentrations of,
are each VCO、VO3Polynomial coefficient of (n)<5) And e represents the fitting residual.
Step S105: for secondary source of HCHO, NO2Carrying out normalization processing on the daily spatial average value of the troposphere vertical column concentration to obtain a corresponding standard value;
in particular, due to secondary sources of HCHO, NO2The absolute values of the concentrations in the environment are different, and in order to facilitate subsequent data processing, HCHO secondary source and NO need to be processed2The normalization processing is carried out on the daily spatial average value of the tropospheric vertical column concentration, and is realized by the following formula:
Vx,i,nor＝Vx,i/Vx,ref
Vx,refobtained from the following equation:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of a target gas x, x being a secondary source of HCHO or NO2，Vx,i,norIs a Vx,iNormalized standard value of (V)x,refDaily reference for tropospheric vertical column concentration of target gas x, d is the total number of days of the study time period.
Step S106: based on O3The daily spatial average value of the tropospheric vertical column concentration, the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source, the standard value thereof, and NO2The daily spatial average value and the standard value of the tropospheric vertical column concentration to obtain the threshold value and NO of the VOCs control areaxA control zone threshold;
the specific process of the step is as follows:
the daily spatial average of tropospheric vertical column concentrations of secondary sources of HCHO on the same day is compared with NO2The ratio of the daily spatial average of tropospheric vertical column concentrations of (a);
all data satisfying the indicia less than the first parametric reference value are selected as a first data set from data elements comprising the O of the target region per day over the study time period3The daily spatial average value of the tropospheric vertical column concentration, the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source, the standard value thereof, and NO2The daily spatial average value and the standard value of the troposphere vertical column concentration;
linear fit first data set lower O3A first slope of the daily spatial average of tropospheric vertical column concentrations of (a) to a standard value of the daily spatial average of tropospheric vertical column concentrations of (b) a secondary source of HCHO, and O3Daily spatial averaging of tropospheric vertical column concentrationsValue and NO2Obtaining a first variable reference value corresponding to the first slope and the second slope being equal, and taking the first variable reference value as a threshold value of the VOCs control area;
selecting all data satisfying the flag smaller than the second reference value from the data elements as a second data set;
linear fit of second data set lower O3A third slope of the daily spatial average of tropospheric vertical column concentrations to a standard value of the daily spatial average of tropospheric vertical column concentrations of the secondary source of HCHO, and O3Daily spatial average of tropospheric vertical column concentration and NO2The fourth slope of the standard value of the daily spatial average value of the tropospheric vertical column concentration of (a), a second variable reference value corresponding to the third slope being equal to the fourth slope is obtained, and the second variable reference value at that time is taken as NOxA control zone threshold.
For example, VHCHO，i/VNO2，iIs used as a flag to mark values less than a certain value m1All data elements of (first variable reference value) are fetched as a first data set, the data elements comprising daily V over the study time periodHCHO、VNO2、VO3、VHCHO,norAnd VNO2,nor；
Then, fitting V under the data set linearlyO3For VHCHO,norThe slope (S1) of (A) and (V)O3For VNO2,norI.e., (V) in the first data set (S2)O3，VHCHO,nor) Set of points and (V)O3，VNO2,nor) Respectively carrying out linear fitting on the point sets to obtain two corresponding linear functions, and adjusting m in the fitting process1Can adjust the corresponding slope, when S1 is equal to S2, m is equal to1I.e., the threshold for the control region of VOCs (T1).
In a similar way, the mark is taken to be larger than a certain value m2(second variable reference value) as a second data set, and linearly fitting V to the data setO3For VHCHO,norThe slope (S3) of (A) and (V)O3For VNO2,norWhen S3 is equal to S4, m is equal to S4 (S4)2Is NOxThreshold of control zone (T2).
Step S107: control zone threshold and NO based on VOCsxAnd controlling zone threshold values, and dividing troposphere ozone control types of all longitude and latitude grid points in the target zone.
The method specifically comprises the following steps:
calculating the troposphere vertical column concentration and NO of the HCHO secondary source in each longitude and latitude grid point in the target area2The tropospheric vertical column concentration ratio;
dividing the ozone control type of the longitude and latitude grid points with the specific value smaller than the threshold value of the VOC control area into VOCs control types;
the ratio is greater than NOxOzone control type classification into NO for latitude and longitude grid points of control zone thresholdxA control type;
the ratio is greater than the VOC control region threshold and less than NOxThe ozone control types of the longitude and latitude grid points of the control area threshold are divided into transition control types.
For example, by using the obtained threshold values T1 and T2, the tropospheric ozone control types are divided in various regions. V if lattice pointHCHO/vNO2Less than T1, the lattice is of the VOC control type. V if lattice pointHCHO/vNO2If greater than T1, the lattice point is NOxThe type of control. V if lattice pointHCHO/vNO2Greater than T1 and less than T2, the grid point is of the transition control type.
Fig. 2 shows the results of performing different stages of ozone control type recognition on a certain area by applying the ozone control type recognition method of the present embodiment, wherein the left side is the spatial distribution map of the ozone control type of the target monitored area between 14 days and 28 days of 10 months in 2020, the middle is the spatial distribution map of the ozone control type of the target monitored area between 29 days and 11 days of 11 months in 2020, and the right side is the spatial distribution map of the ozone control type of the target monitored area between 12 days and 30 days of 11 months in 2020. In the figure, the abscissa is longitude, the ordinate is latitude, and RATIO is HCHO secondary source troposphere vertical column concentration and NO2Tropospheric sagThe ratio of the concentrations of the straight columns, the threshold value T1 of the VOCs control area is 1.0, NOxThe threshold value T2 of the control area is 2.0, and the ozone control type corresponding to each longitude and latitude grid point can be accurately obtained from the graph.
According to the ozone control type identification method based on satellite hyperspectral remote sensing, the ozone control types of all places are analyzed and classified in a large range through the hyperspectral remote sensing satellite, main ozone precursors of all places can be determined, and then the ozone concentration can be controlled scientifically and effectively.
Example 2
This embodiment proposes an electronic device, the electronic device including:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the ozone control type identification method of embodiment 1 above.
An electronic device according to an embodiment of the disclosure includes a memory for storing non-transitory computer readable instructions and a processor. In particular, the memory may include one or more computer program products that may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, Random Access Memory (RAM), cache memory (cache), and/or the like. The non-volatile memory may include, for example, Read Only Memory (ROM), hard disk, flash memory, etc.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing capabilities and/or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures should also be included in the protection scope of the present disclosure.
For the detailed description of the present embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which are not repeated herein.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. An ozone control type identification method based on satellite hyperspectral remote sensing is characterized by comprising the following steps:
acquiring troposphere vertical column concentration data of target gases in a target region in a research time phase based on atmospheric spectrum remote sensing data returned by a satellite, wherein the target gases comprise CO and O3HCHO and NO2；
Regulating the pixel of each target gas to the equal longitude and latitude grid of the target area to obtain the corresponding troposphere gas vertical column concentration of each target gas in each longitude and latitude grid point on the equal longitude and latitude grid;
obtaining available pixels in the target area, and calculating the daily spatial average value of the troposphere vertical column concentration of each target gas of all the available pixels;
based on CO and O3And daily spatial average value data of the troposphere vertical column concentration of HCHO, and calculating the daily spatial average value of the troposphere vertical column concentration of HCHO secondary sources of all available pixels, wherein the HCHO secondary sources are vertical column concentrations generated by photochemical reaction of troposphere formaldehyde;
for secondary source of HCHO, NO2Carrying out normalization processing on the daily spatial average value of the troposphere vertical column concentration to obtain a corresponding standard value;
based on O3Daily spatial average of tropospheric vertical column concentrations, HCHO twoDaily spatial average of tropospheric vertical column concentrations of secondary sources and its standard value and NO2The daily spatial average value and the standard value of the tropospheric vertical column concentration to obtain the threshold value and NO of the VOCs control areaxA control zone threshold;
based on the VOCs control region threshold and the NOxAnd controlling zone threshold values, and dividing troposphere ozone control types of all longitude and latitude grid points in the target zone.
2. The ozone control type identification method of claim 1, wherein the arranging the pixels of each target gas onto the equal latitude and longitude grid of the target area comprises:
performing linear interpolation on the longitude and latitude of any two adjacent pixels on the original pixel plane of the satellite and the vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derived pixels and the original pixels into corresponding longitude and latitude grid points;
and averaging the vertical column concentration of each target gas pixel in each longitude and latitude grid point to obtain the vertical column concentration value of each target gas on each longitude and latitude grid point.
3. The ozone control type identification method of claim 2, wherein the vertical column concentration of each target gas pixel in each longitude and latitude grid point is averaged by the following formula:
wherein, ViThe vertical column concentration of the troposphere gas, lon, corresponding to the target gas on the ith pixel needing to be processed on the original pixel plane of the satellite after linear interpolationi，latiIs the longitude and latitude, v, of the ith pixellon，latThe vertical column concentration of the troposphere corresponding to the target gas in the longitude and latitude grid points corresponding to the longitude and latitude on the equal longitude and latitude grid is represented by lon and lat which are the longitude and latitude of the ith pixel and grid which isThe degrees of difference between the longitude and latitude grid points.
4. The ozone control type identification method of claim 1, wherein the obtaining of the available pixels in the target area comprises:
and removing the pixels with the cloud amount larger than 0.5 according to the cloud amount information of each pixel in the target area returned by the satellite, and taking the rest pixels as the available pixels.
5. The ozone control type identification method of claim 1, wherein the calculating of the daily spatial average of tropospheric vertical column concentrations of each target gas for all available picture elements is performed by the following formula:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O3Secondary sources of HCHO or NO2I represents the ith day in the set time period, the { region } is a set of all the equal longitude and latitude grid points in the target area, N is the number of the longitude and latitude grid points in the { region }, lon and lat are the longitude and latitude of a single longitude and latitude grid point, and v is the longitude and latitude of the single longitude and latitude grid pointx,i,lon,latIs the tropospheric gas vertical column concentration of the target gas at day i in a single longitude and latitude point.
6. The ozone control type identification method of claim 1, wherein the calculation of the daily spatial average of tropospheric vertical column concentrations for all available pixel elements HCHO secondary sources is performed by the following multivariate polynomial fitting operation:
VHCHO＝(1(VCO)1…(VCO)n)A+(1(VO3)1…(VO3)n)B+e
wherein Vp,HCHODaily spatial average of tropospheric vertical column concentrations, V, for secondary sources of HCHOHCHOIs the daily spatial average of the total vertical column concentration of tropospheric formaldehyde, VCO、VO3Are respectively CO and O3The daily spatial average of tropospheric vertical column concentrations of,
are each VCO、VO3Polynomial coefficient of (n)<5) And e represents the fitting residual.
7. The ozone control type identification method of claim 1, wherein the pair of secondary sources of HCHO, NO2The normalization process is carried out on the daily spatial average value of the tropospheric vertical column concentration, and is realized by the following formula:
Vx,i,nor＝Vx,i/Vx,ref
Vx,refobtained from the following equation:
wherein, Vx,iIs the daily spatial average of the tropospheric vertical column concentration of a target gas x, x being a secondary source of HCHO or NO2，Vx,i,norIs a Vx,iNormalized standard value of (V)x,refDaily reference value for tropospheric vertical column concentration of target gas x, d is total days of the study time period.
8. The ozone control type identification method of claim 1, wherein the VOCs control zone threshold and the NO are obtainedxThe method for controlling the zone threshold value comprises the following steps:
daily spatial averaging of tropospheric vertical column concentrations of secondary sources of HCHO on the same dayValue and NO2The ratio of the daily spatial average of tropospheric vertical column concentrations of (a);
selecting as the first data set all data satisfying that the indicia is less than a first parametric reference value from data elements comprising O for each day of the target region over the study time period3The daily spatial average value of the tropospheric vertical column concentration, the daily spatial average value of the tropospheric vertical column concentration of the HCHO secondary source, the standard value thereof, and NO2The daily spatial average value and the standard value of the troposphere vertical column concentration;
linearly fitting the first data set to O3A first slope of the daily spatial average of tropospheric vertical column concentrations of (a) to a standard value of the daily spatial average of tropospheric vertical column concentrations of (b) a secondary source of HCHO, and O3Daily spatial average of tropospheric vertical column concentration and NO2Obtaining a first variable reference value corresponding to the first slope and the second slope being equal, and taking the first variable reference value at the moment as the threshold of the VOCs control area;
selecting all data satisfying that the flag is smaller than a second reference variable value from the data elements as a second data set;
linearly fitting the second data set to O3A third slope of the daily spatial average of tropospheric vertical column concentrations to a standard value of the daily spatial average of tropospheric vertical column concentrations of the secondary source of HCHO, and O3Daily spatial average of tropospheric vertical column concentration and NO2A fourth slope of the standard value of the daily spatial average of tropospheric vertical column concentrations of (a), a second variable reference value corresponding to the third slope being equal to the fourth slope is obtained, and the second variable reference value at that time is taken as the NOxA control zone threshold.
9. The ozone control type identification method of claim 1, wherein the control zone threshold based on the VOCs and the NO arexControl ofAnd the zone threshold value is used for dividing troposphere ozone control types of all longitude and latitude grid points in the target zone, and comprises the following steps:
calculating the troposphere vertical column concentration and NO of the HCHO secondary source in each longitude and latitude grid point in the target area2The tropospheric vertical column concentration ratio;
dividing the ozone control type of the longitude and latitude grid point with the ratio smaller than the threshold value of the VOCs control area into VOCs control types;
will be greater than the NOxOzone control type classification into NO for latitude and longitude grid points of control zone thresholdxA control type;
the ratio is larger than the threshold value of the VOCs control area and smaller than the NOxThe ozone control types of the longitude and latitude grid points of the threshold value of the control area are divided into transition control types.
10. An electronic device, characterized in that the electronic device comprises:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the ozone control type identification method of any one of claims 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111476909.0A CN114496117B (en) | 2021-12-02 | 2021-12-02 | Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111476909.0A CN114496117B (en) | 2021-12-02 | 2021-12-02 | Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114496117A true CN114496117A (en) | 2022-05-13 |
| CN114496117B CN114496117B (en) | 2024-02-23 |
Family
ID=81492405
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111476909.0A Active CN114496117B (en) | 2021-12-02 | 2021-12-02 | Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114496117B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140103229A (en) * | 2013-02-15 | 2014-08-26 | 연세대학교 산학협력단 | Apparatus for analyzing atmosphere |
| KR20140115777A (en) * | 2013-03-22 | 2014-10-01 | 이화여자대학교 산학협력단 | Method for estimating longwave climate feedback using sea surface temperature |
| CN110942049A (en) * | 2019-12-17 | 2020-03-31 | 生态环境部卫星环境应用中心 | Ozone pollution source identification method and system based on satellite remote sensing |
| CN111678880A (en) * | 2020-06-04 | 2020-09-18 | 生态环境部卫星环境应用中心 | Satellite remote sensing identification method and system for stratospheric ozone invading lower layer in troposphere |
| CN112183384A (en) * | 2020-09-29 | 2021-01-05 | 南京信息工程大学 | Aerosol remote sensing classification method based on satellite data |
| CN113176216A (en) * | 2021-01-12 | 2021-07-27 | 生态环境部卫星环境应用中心 | Ozone precursor VOCs high-value area satellite remote sensing identification method |
| CN113552081A (en) * | 2021-07-06 | 2021-10-26 | 中国科学技术大学 | Remote measurement system based on ultra-high spectrum remote sensing non-blind area ozone vertical distribution |
-
2021
- 2021-12-02 CN CN202111476909.0A patent/CN114496117B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140103229A (en) * | 2013-02-15 | 2014-08-26 | 연세대학교 산학협력단 | Apparatus for analyzing atmosphere |
| KR20140115777A (en) * | 2013-03-22 | 2014-10-01 | 이화여자대학교 산학협력단 | Method for estimating longwave climate feedback using sea surface temperature |
| CN110942049A (en) * | 2019-12-17 | 2020-03-31 | 生态环境部卫星环境应用中心 | Ozone pollution source identification method and system based on satellite remote sensing |
| CN111678880A (en) * | 2020-06-04 | 2020-09-18 | 生态环境部卫星环境应用中心 | Satellite remote sensing identification method and system for stratospheric ozone invading lower layer in troposphere |
| CN112183384A (en) * | 2020-09-29 | 2021-01-05 | 南京信息工程大学 | Aerosol remote sensing classification method based on satellite data |
| CN113176216A (en) * | 2021-01-12 | 2021-07-27 | 生态环境部卫星环境应用中心 | Ozone precursor VOCs high-value area satellite remote sensing identification method |
| CN113552081A (en) * | 2021-07-06 | 2021-10-26 | 中国科学技术大学 | Remote measurement system based on ultra-high spectrum remote sensing non-blind area ozone vertical distribution |
Non-Patent Citations (2)
| Title |
|---|
| 应静瑶;陈奔;孙兴富;: "温岭市臭氧浓度的控制区类型研究及其应用", 能源环境保护, no. 05, 15 October 2018 (2018-10-15), pages 59 - 62 * |
| 陈雪萍;巨天珍;张江峪;咸龙;王培玉;张生财;: "兰州市大气臭氧与其前体物及气象因子相关性研究", 环境科学研究, no. 12, 31 December 2019 (2019-12-31), pages 123 - 131 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114496117B (en) | 2024-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pendergrass et al. | Surface and top-of-atmosphere radiative feedback kernels for CESM-CAM5 | |
| Bax et al. | How can odors be measured? An overview of methods and their applications | |
| CN110298560B (en) | Method and device for evaluating atmospheric pollution emission control effect and storage medium | |
| Tadic et al. | Net ozone production and its relationship to nitrogen oxides and volatile organic compounds in the marine boundary layer around the Arabian Peninsula | |
| Capelli et al. | Electronic noses for environmental monitoring applications | |
| Jia et al. | Global inorganic nitrogen dry deposition inferred from ground-and space-based measurements | |
| Ukkola et al. | Reduced streamflow in water-stressed climates consistent with CO2 effects on vegetation | |
| Kong et al. | Evaluation and uncertainty investigation of the NO 2, CO and NH 3 modeling over China under the framework of MICS-Asia III | |
| Mao et al. | High resolution modeling of gaseous methylamines over a polluted region in China: Source-dependent emissions and implications to spatial variations | |
| Thorson et al. | Using a low-cost sensor array and machine learning techniques to detect complex pollutant mixtures and identify likely sources | |
| Huang et al. | Soil nitric oxide emissions from terrestrial ecosystems in China: a synthesis of modeling and measurements | |
| Lama et al. | Quantifying burning efficiency in Megacities using NO 2/CO ratio from the Tropospheric Monitoring Instrument (TROPOMI) | |
| Harris et al. | Tracking nitrous oxide emission processes at a suburban site with semicontinuous, in situ measurements of isotopic composition | |
| CN113155939A (en) | Online volatile organic compound source analysis method, system, equipment and medium | |
| Vikram et al. | Evaluating and improving the reliability of gas-phase sensor system calibrations across new locations for ambient measurements and personal exposure monitoring | |
| Zhou et al. | Kinetics and impacting factors of HO 2 uptake onto submicron atmospheric aerosols during the 2019 Air QUAlity Study (AQUAS) in Yokohama, Japan | |
| Tancev et al. | The relocation problem of field calibrated low-cost sensor systems in air quality monitoring: a sampling bias | |
| Zhang et al. | Impacts of climate change and emissions on atmospheric oxidized nitrogen deposition over East Asia | |
| Tang et al. | Satellite data reveal a common combustion emission pathway for major cities in China | |
| Tomicic et al. | Experimental study of H 2 SO 4 aerosol nucleation at high ionization levels | |
| Zhao et al. | Significant climate impacts of aerosol changes driven by growth in energy use and advances in emission control technology | |
| Kim et al. | Nitrate chemistry in the northeast US–Part 2: Oxygen isotopes reveal differences in particulate and gas-phase formation | |
| CN114496117B (en) | Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment | |
| Denzler et al. | Urban source term estimation for mercury using a boundary-layer budget method | |
| Taguem et al. | Quantile regression with a metal oxide sensors array for methane prediction over a municipal solid waste treatment plant |
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 |