CN114496117B - 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
- CN114496117B CN114496117B CN202111476909.0A CN202111476909A CN114496117B CN 114496117 B CN114496117 B CN 114496117B CN 202111476909 A CN202111476909 A CN 202111476909A CN 114496117 B CN114496117 B CN 114496117B
- Authority
- CN
- China
- Prior art keywords
- vertical column
- tropospheric
- spatial average
- hcho
- column concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 46
- WSFSSNUMVMOOMR-UHFFFAOYSA-N Formaldehyde Chemical compound O=C WSFSSNUMVMOOMR-UHFFFAOYSA-N 0.000 claims abstract description 98
- 239000012855 volatile organic compound Substances 0.000 claims abstract description 35
- 239000005436 troposphere Substances 0.000 claims abstract description 32
- 229910002091 carbon monoxide Inorganic materials 0.000 claims abstract description 28
- 238000011160 research Methods 0.000 claims abstract description 5
- 238000012935 Averaging Methods 0.000 claims description 6
- 238000006552 photochemical reaction Methods 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 5
- 230000001105 regulatory effect Effects 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- WSFSSNUMVMOOMR-NJFSPNSNSA-N methanone Chemical compound O=[14CH2] WSFSSNUMVMOOMR-NJFSPNSNSA-N 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 2
- 239000002243 precursor Substances 0.000 abstract description 9
- 239000007789 gas Substances 0.000 description 55
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 24
- 238000010586 diagram Methods 0.000 description 6
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 239000000356 contaminant Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion 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
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000002440 industrial waste Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000003960 organic solvent Substances 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 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
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 and electronic equipment based on satellite hyperspectral remote sensing, wherein the method comprises the following steps: acquiring troposphere vertical column concentration data of each target gas in a research time period of a target area; obtaining the corresponding concentration of the troposphere gas vertical column of each target gas in each theodolite point on the equal longitude and latitude grid; calculating a daily spatial average of tropospheric vertical column concentrations for each target gas for all available picture elements; based on CO, O 3 And daily spatial average data of the tropospheric vertical column concentration of the HCHO, calculating a daily spatial average of the tropospheric vertical column concentration of the HCHO secondary source of all available pels; obtaining VOCs control region threshold and NO x A control region threshold; VOCs control zone threshold and NO based x The control area threshold value divides the troposphere ozone control type of each theodolite point in the target area. Realizes the definition of main ozone precursors in all 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 and electronic equipment based on satellite hyperspectral remote sensing.
Background
Ozone present in the troposphere is a biologically harmful contaminant and is one of the components of photochemical smog. Many human activities involving rapid conversion of chemical energy, such as engine start-up and copier operation, produce ozone, a strong oxidizer that reacts easily with other chemicals to form many toxic oxides that are harmful to human health. Thus, the concentration of tropospheric ozone is closely related to the production activity and physical health of people.
Tropospheric ozone is mainly derived from photochemical reactions-when air mixed with various oxides of nitrogen (NOx), carbon monoxide (CO) and volatile organic compounds (VOCs, such as formaldehyde-HCHO) is exposed to sunlight, ozone is produced. Nitrogen oxides and volatile organics are therefore referred to as "ozone precursors". Automobile exhaust, industrial waste gas and chemical organic solvents are the main artificial emission sources of "ozone precursors". Although these emissions sources are mostly concentrated in cities, some substances (such as nitrogen oxides) can diffuse by wind power to sparse areas of the population, hundreds of kilometers away, where ozone sources are formed.
It is now possible to measure the vertical column (vertical atmospheric column) concentration of various atmospheric contaminants in the troposphere using satellites. Satellite observation can realize observation of a large range of various gas concentrations, and is convenient for wide-area analysis of various gas concentration relations and control of pollutants.
The existing control method for the concentration of the ozone in the troposphere mainly controls the high-concentration troposphere ozone point source, but is influenced by the photochemical reaction of the ozone, and the control for the troposphere ozone point source cannot effectively realize the reduction of the concentration of the ozone in the troposphere 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 realize the large-scale determination of ozone in various places on two precursors VOCs and NO thereof x And to obtain the main influence of what precursor the ozone concentration of the troposphere is subjected to.
In order to achieve the above purpose, the invention provides an ozone control type identification method based on satellite hyperspectral remote sensing, which comprises the following steps:
based on atmospheric spectrum remote sensing data returned by satellites, tropospheric vertical column concentration data of each target gas in a target area in a research time period are obtained, wherein the target gas comprises CO and O 3 HCHO and NO 2 ；
The pixels of each target gas are regulated to equal longitude and latitude grids of the target area, and the corresponding troposphere gas vertical column concentration of each target gas in each theodolite point on the equal longitude and latitude grids is obtained;
acquiring available pixels in the target area, and calculating a daily spatial average value of the tropospheric vertical column concentration of each target gas of all the available pixels;
based on CO, O 3 And the daily spatial average value data of the tropospheric vertical column concentration of HCHO, calculating the daily spatial average value of the tropospheric vertical column concentration of all available pixel HCHO secondary sources, wherein the HCHO secondary sources are the vertical column concentrations generated by photochemical reaction of tropospheric formaldehyde;
for HCHO secondary source, NO 2 Carrying out normalization treatment on the daily space average value of the troposphere vertical column concentration to obtain a corresponding standard value;
based on O 3 Daily spatial average of tropospheric vertical column concentration, daily spatial average of tropospheric vertical column concentration of HCHO secondary source and standard value and NO thereof 2 Daily spatial average value of tropospheric vertical column concentration and standard value thereof, obtaining a threshold value of a VOCs control zone and NO x A control region threshold;
based on the VOCs control region threshold and the NO x And (5) controlling the zone threshold value, and dividing the troposphere ozone control type of each theodolite point in the target zone.
Optionally, the step of regularizing the pixels of each target gas onto the equal longitude and latitude grid of the target area includes:
linearly interpolating longitude and latitude of any two adjacent pixels on the satellite original pixel plane and vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derivative pixel and the original pixel to corresponding theodolite points;
and averaging the vertical column concentration of each target gas pixel in each theodolite point to obtain the vertical column concentration value of each target gas on each theodolite point.
Optionally, the averaging the vertical column concentration of each target gas pixel in each theodolite point is achieved by the following formula:
wherein V is i Represents the vertical column concentration of the tropospheric gas corresponding to the target gas on the ith pixel to be processed on the plane of the original pixel of the satellite after linear interpolation, lon i ，lat i Is the longitude and latitude of the ith pixel, v lin，lat For the concentration of the troposphere vertical column corresponding to the target gas in the corresponding theodolite points on the equal longitude and latitude grid, lon and lat are the longitude and latitude of the ith pixel, and grid is the degree of phase difference between the theodolite points.
Optionally, the acquiring the available pixels in the target area includes:
and removing pixels with cloud quantity larger than 0.5 according to cloud quantity information of each pixel in the target area returned by the satellite, and taking the rest pixels as the available pixels.
Optionally, the calculating of the daily spatial average of the tropospheric vertical column concentration for each target gas for all available picture elements is accomplished by the following formula:
wherein V is x,i For the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O 3 Secondary source of HCHO or NO 2 I represents the ith day in the set time period, { region } is the set of all equal theodolite points in the target area, N is the number of the theodolite points in { region }, lon and lat are the longitude and latitude of a single theodolite point, v x,i,lon,lat The tropospheric gas vertical column concentration on day i in a single theodolite point for the target gas.
Optionally, the calculating of the daily spatial average of tropospheric vertical column concentrations for all available pel HCHO secondary sources is accomplished by the following polynomial fitting operation:
V HCHO ＝(1(V CO ) 1 …(V CO ) n )A+(1(V O3 ) 1 …(V O3 ) n )B+e
wherein V is p,HCHO Daily spatial average of tropospheric vertical column concentration for HCHO secondary source, V HCHO Daily spatial average of total vertical column concentration of tropospheric formaldehyde, V CO 、V O3 CO and O respectively 3 Daily spatial average of tropospheric vertical column concentration,v respectively CO 、V O3 Polynomial coefficients (n)<5) E represents the fit residual.
Optionally, the secondary source of HCHO, NO 2 The daily spatial average value of the tropospheric vertical column concentration is normalized by the following formula:
V x,i,nor ＝V x,i /V x,ref
V x,ref the method is obtained by the following formula:
wherein V is x,i Daily spatial average of tropospheric vertical column concentration for target gas x, x being HCHO secondary source or NO 2 ，V x,i,nor Is V (V) x,i Normalized standard value of V x,ref Daily reference value for tropospheric vertical column concentration of target gas x, d is total number of days of the study time period.
Optionally, obtaining the VOCs control region threshold and the NO x The method for controlling the zone threshold comprises the following steps:
will H on the same dayDaily spatial average of tropospheric vertical column concentration of CHO secondary source with NO 2 The ratio of the daily spatial average of tropospheric vertical column concentrations is used as a marker;
selecting as a first data set all data satisfying said signature less than a first parameter variable reference value from data elements comprising a daily spatial average value of the tropospheric vertical column concentration of O3, a daily spatial average value of the tropospheric vertical column concentration of a HCHO secondary source, a standard value thereof, and NO for each day of said target area during said time period of investigation 2 Daily spatial average of tropospheric vertical column concentration and standard value thereof;
linear fitting of O under the first dataset 3 A first slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration for a secondary source of HCHO, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A second slope of a standard value of a daily spatial average value of tropospheric vertical column concentration, obtaining a first variable reference value corresponding to when the first slope is equal to the second slope, and taking the first variable reference value at that time as the VOCs control zone threshold;
selecting all data meeting the condition that the mark is smaller than a second parameter variable reference value from the data elements as a second data set;
linear fitting of O under the second dataset 3 A third slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration of the HCHO secondary source, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A fourth slope of a standard value of a daily spatial average value of tropospheric vertical column concentration, obtaining a second variable reference value corresponding to when the third slope is equal to the fourth slope, and taking the second variable reference value at that time as the NO x A control region threshold.
Optionally, the control region threshold and the NO based on the VOCs x A control area threshold value for dividing each theodolite point in the target areaA tropospheric ozone control type comprising:
calculating the troposphere vertical column concentration and NO of the HCHO secondary source in each theodolite point in the target area 2 The ratio of tropospheric vertical column concentration;
dividing the ozone control type of the theodolite point with the ratio smaller than the threshold value of the VOCs control area into the VOCs control type;
will have a ratio greater than that of the NO x Ozone control type division of theodolite of control zone threshold into NO x A control type;
will have a ratio greater than the VOCs control region threshold and less than the NO x Ozone control types of theodolite points of the control zone threshold are classified as transition control types.
The invention also proposes an electronic device comprising:
at least one processor; the method comprises the steps of,
a memory communicatively coupled to the at least one processor; wherein,
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 invention firstly obtains each target gas CO, O in the research time period of the target region 3 Tropospheric vertical column concentration data of HCHO and NO2, obtaining corresponding tropospheric vertical column concentrations of each target gas in each theodolite point on the equal longitude and latitude grid, calculating daily spatial average value of tropospheric vertical column concentrations of each target gas of all available pixels, calculating daily spatial average value of tropospheric vertical column concentrations of HCHO secondary sources of all available pixels, and then based on O 3 Daily spatial average of tropospheric vertical column concentration, daily spatial average of tropospheric vertical column concentration of HCHO secondary source and standard value and NO thereof 2 Daily spatial average value of tropospheric vertical column concentration and standard value thereofObtaining the threshold value of the VOCs control region and NO x Control zone threshold, and finally based on VOCs control zone threshold and NO x The threshold value of the control area is used for dividing the troposphere ozone control type of each theodolite point in the target area, thereby realizing a method suitable for satellites and measuring the ozone of each place in a large range for two precursors VOCs and NO thereof x The method for analyzing the sensitivity of the method can clearly study the main influence of the precursor on the ozone concentration of the troposphere at each place in the area, and can specifically control the ozone concentration of different areas 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 foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
Fig. 1 shows a step diagram of an ozone control type identification method based on satellite hyperspectral remote sensing according to the present invention.
Fig. 2 shows a distribution diagram of ozone control types for different phases of the same target area in an ozone control type identification method based on satellite hyperspectral remote sensing according to embodiment 1 of the present 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 illustrated 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 present invention.
As shown in fig. 1, an ozone control type identification method based on satellite hyperspectral remote sensing includes:
step S101: atmospheric spectrum remote sensing data returned based on satellites are used for acquiring tropospheric vertical column concentration data of each target gas in a target area in a research time period, wherein the target gas comprises CO and O 3 HCHO and NO 2 ；
Specifically, various gas troposphere vertical column concentrations are obtained by acquiring various gas secondary products returned by satellites.
Step S102: the pixels of each target gas are regulated to equal longitude and latitude grids of the target area, and the corresponding troposphere gas vertical column concentration of each target gas in each theodolite point on the equal longitude and latitude grids is obtained;
in this step, the pixel of each target gas is regulated to the equal longitude and latitude grid of the target area, including:
linearly interpolating longitude and latitude of any two adjacent pixels on the satellite original pixel plane and vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derivative pixel and the original pixel to corresponding theodolite points;
and averaging the vertical column concentration of each target gas pixel in each theodolite point to obtain the vertical column concentration value of each target gas on each theodolite point.
Specifically, in the data analysis process, multiple gas data are required, and the pixel positions of each gas product returned by the satellite are different. In order to ensure the accuracy of the subsequent process, pixels of various gases need to be regulated to be on equal longitude and latitude grids, and a sampling method is used for processing:
the original pixel plane of the satellite has two dimensions X-Y' corresponding to the CCD direction of the satellite instrument and the satellite travelling direction, and the longitude and latitude and the gas concentration of any two adjacent pixels on the plane are linearly interpolated to obtain the derivative pixel. Dividing the derivative pixels and the original pixels into corresponding theodolite points, and averaging the gas concentration of the pixels in each theodolite point to obtain the gas concentration value of each theodolite point.
In the above process, the vertical column concentration of each target gas pixel in each theodolite point is averaged, and is realized by the following formula:
wherein V is i Represents the vertical column concentration of the tropospheric gas corresponding to the target gas on the ith pixel to be processed on the plane of the original pixel of the satellite after linear interpolation, lon i ，lat i Is the longitude and latitude of the ith pixel, v lon，lat For the concentration of the troposphere vertical column corresponding to the target gas in the corresponding theodolite points on the equal longitude and latitude grid, lon and lat are the longitude and latitude of the ith pixel, and grid is the degree of phase difference between the theodolite points.
Step S103: acquiring available pixels in a target area, and calculating a daily spatial average value of the tropospheric vertical column concentration of each target gas of all the available pixels;
in particular, since cloud cover can seriously affect the accuracy of trace gas inversion, which can affect the sensitivity analysis to ozone, it is desirable to remove unusable pixels affected by cloud cover. In this step, the method for acquiring the available pixels in the target area includes:
and removing pixels with cloud quantity larger than 0.5 according to cloud quantity information of each pixel in a target area returned by the satellite, and taking the rest pixels as available pixels.
After the available pixels are obtained, selecting a target area, and enabling all the available pixels in the target area to be O 3 Secondary source of HCHO, NO 2 The concentration of the troposphere vertical column is averaged and divided into V by day x,i 。
In this step, the daily spatial average of the tropospheric vertical column concentration for each target gas for all available pels is calculated by the following formula:
wherein V is x,i For the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O3, HCHO secondary source or NO 2 I represents the ith day in the set time period, { region } is the set of all equal theodolite points in the target area, N is the number of the theodolite points in { region }, lon and lat are the longitude and latitude of a single theodolite point, v x,i,lon,lat The tropospheric gas vertical column concentration on day i in a single theodolite point for the target gas.
Step S104: based on CO, O 3 And daily spatial average value data of the tropospheric vertical column concentration of HCHO, calculating the daily spatial average value of the tropospheric vertical column concentration of all available pixel HCHO secondary sources, wherein the HCHO secondary sources are the vertical column concentrations generated by photochemical reaction of the tropospheric formaldehyde;
in particular, VOCs and NO in the atmosphere due to photochemical reactions x Will gradually generate O 3 Therefore, the concentration of HCHO and O as the VOC concentration indicator 3 There is a relationship in concentration. Meanwhile, HCHO is also generated in the incomplete combustion process of various biomasses, and CO is generated in the process, so that the HCHO concentration and the CO concentration have a relationship at the same time. The HCHO concentration in the atmosphere can be expressed by the following formula:
V HCHO -V back ＝F(V O3 ,V CO )
wherein V is HCHO For HCHO concentration in the atmosphere, V back Is the background concentration value of atmospheric HCHO, V O3 、V CO O in the atmosphere respectively 3 Concentration and CO concentration.
HCHO, O in a plurality of days in a certain place 3 The concentration of CO is used as data element, and high-order polynomial is used for fitting, and the proportion of various sources of formaldehyde is analyzed through the fitting of the polynomials to establish HCHO and O 3 、CO, and determining the concentration of the troposphere formaldehyde in the vertical column generated by photochemical reaction, namely the concentration of the vertical column of the HCHO secondary source.
The method is realized by the following polynomial fitting operation:
V HCHO ＝(1(V CO ) 1 …(V CO ) n )A+(1(V O3 ) 1 …(V O3 ) n )B+e
wherein V is p,HCHO Daily spatial average of tropospheric vertical column concentration for HCHO secondary source, V HCHO Daily spatial average of total vertical column concentration of tropospheric formaldehyde, V CO 、V O3 CO and O respectively 3 Daily spatial average of tropospheric vertical column concentration,v respectively CO 、V O3 Polynomial coefficients (n)<5) E represents the fit residual.
Step S105: for HCHO secondary source, NO 2 Carrying out normalization treatment on the daily space average value of the troposphere vertical column concentration to obtain a corresponding standard value;
in particular, due to HCHO secondary source, NO 2 The absolute values of the concentrations in the environment are different, and the HCHO secondary source and NO are needed to be processed for the subsequent data processing 2 The daily spatial average value of the tropospheric vertical column concentration is normalized by the following formula:
V x,i,nor ＝V x,i /V x,ref
V x,ref the method is obtained by the following formula:
wherein V is x,i Daily spatial average of tropospheric vertical column concentration for target gas x, x being HCHO secondary source or NO 2 ，V x,i,nor Is V (V) x,i Normalized standard value of V x,ref The daily reference value for the tropospheric vertical column concentration of the target gas x, d is the total number of days of the study time period.
Step S106: based on O 3 Daily spatial average of tropospheric vertical column concentration, daily spatial average of tropospheric vertical column concentration of HCHO secondary source and standard value and NO thereof 2 Daily spatial average value of tropospheric vertical column concentration and standard value thereof, obtaining a threshold value of a VOCs control zone and NO x A control region threshold;
the specific process of the step is as follows:
daily spatial average of tropospheric vertical column concentration of HCHO secondary source on the same day was compared with NO 2 The ratio of the daily spatial average of tropospheric vertical column concentrations is used as a marker;
selecting as a first data set all data satisfying a signature less than a first parameter variable reference value from a data element comprising O for a target area per day during a study time period 3 Daily spatial average of tropospheric vertical column concentration, daily spatial average of tropospheric vertical column concentration of HCHO secondary source and standard value and NO thereof 2 Daily spatial average of tropospheric vertical column concentration and standard value thereof;
linear fitting of O under first dataset 3 A first slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration for a secondary source of HCHO, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A second slope of a standard value of a daily spatial average value of tropospheric vertical column concentration, obtaining a first variable reference value corresponding to the first slope being equal to the second slope, and taking the first variable reference value at that time as a VOCs control zone threshold;
selecting all data meeting the condition that the mark is smaller than the second parameter variable reference value from the data elements as a second data set;
linear fitting of O under second dataset 3 A third slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration of the HCHO secondary source, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A fourth slope of the standard value of the daily spatial average value of the tropospheric vertical column concentration, obtaining a second variable reference value corresponding to the third slope being equal to the fourth slope, and taking the second variable reference value at this time as NO x A control region threshold.
For example, V HCHO，i /V NO2，i As a marker, the marker is smaller than a certain value m 1 All data elements (first variable reference value) are fetched as a first data set, the data elements comprising daily V during the study time period HCHO 、V NO2 、V O3 、V HCHO,nor And V NO2,nor ；
Then, the data set V is linearly fitted O3 For V HCHO,nor Slope (S1) and V O3 For V NO2,nor I.e. for the slope (S2) of the first dataset (V O3 ，V HCHO , nor ) Point set sum (V) O3 ，V NO2,nor ) Respectively performing linear fitting on the point sets to obtain two corresponding primary functions, and adjusting m in the fitting process 1 Can adjust the corresponding slope, when s1=s2, m 1 I.e., the threshold (T1) of the VOCs control zone.
In a similar way, the mark is taken to be larger than a certain value m 2 All data elements (second variable reference value) are taken as a second data set, and V under the data set is linearly fitted O3 For V HCHO,nor Slope (S3) and V O3 For V NO2,nor Is (S4), when s3=s4, m 2 Namely NO x Threshold (T2) of control zone.
Step S107: VOCs control zone threshold and NO based x And controlling the zone threshold value, and dividing the troposphere ozone control type of each theodolite point in the target zone.
The method specifically comprises the following steps:
calculating the concentration and NO of a troposphere vertical column of an HCHO secondary source in each theodolite point in a target area 2 The ratio of tropospheric vertical column concentration;
dividing ozone control types of theodolite points with the ratio smaller than the threshold value of the VOC control area into VOCs control types;
will have a ratio greater than NO x Ozone control type division of theodolite of control zone threshold into NO x A control type;
the ratio is greater than the VOC control zone threshold and less than NO x Ozone control types of theodolite points of the control zone threshold are classified as transition control types.
For example, the tropospheric ozone control types are divided in various places by using the obtained threshold values T1 and T2. V of the lattice point HCHO /v NO2 Less than T1, the grid point is of the VOC control type. V of the lattice point HCHO /v NO2 Greater than T1, the lattice point is NO x Control type. V of the lattice point HCHO /v NO2 Greater than T1 and less than T2, then the trellis point is of the transition control type.
Fig. 2 shows the results of performing different-stage ozone control type identification for a certain region by applying the ozone control type identification method of the present embodiment, wherein the left side is the ozone control type space distribution diagram of the target monitoring region between 14 days and 28 days in 10 months in 2020, the middle is the ozone control type space distribution diagram of the target monitoring region between 29 days and 11 months in 10 months in 2020, and the right side is the ozone control type space distribution diagram of the target monitoring region between 12 days and 30 days in 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 NO 2 The ratio of the troposphere vertical column concentration, the threshold T1 of the VOCs control zone is 1.0, NO x The threshold value T2 of the control area is 2.0, and the ozone control type corresponding to each theodolite point can be accurately obtained from the graph.
According to the method for identifying the ozone control type based on the hyperspectral remote sensing of the satellite, the hyperspectral remote sensing satellite is used for analyzing and classifying the ozone control type of each place in a large range, so that the main ozone precursor of each place can be clarified, and the ozone concentration can be controlled scientifically and effectively.
Example 2
The embodiment provides an electronic device, including:
at least one processor; the method comprises the steps of,
a memory communicatively coupled to the at least one processor; wherein,
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 present 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, which 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) and/or cache memory (cache), and the like. The non-volatile memory may include, for example, read Only Memory (ROM), hard disk, flash memory, and the like.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing and/or instruction execution capabilities, and may control other components in the electronic device to perform the desired functions. In one embodiment of the present disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
It should be understood by those skilled in the art that, in order to solve the technical problem of how to obtain a good user experience effect, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures are also included in the protection scope of the present disclosure.
The detailed description of the present embodiment may refer to the corresponding description in the foregoing embodiments, and will not be repeated herein.
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or 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 various embodiments described.
Claims (10)
1. The method for identifying the ozone control type based on the satellite hyperspectral remote sensing is characterized by comprising the following steps of:
based on atmospheric spectrum remote sensing data returned by satellites, tropospheric vertical column concentration data of each target gas in a target area in a research time period are obtained, wherein the target gas comprises CO and O 3 HCHO and NO 2 ；
The pixels of each target gas are regulated to equal longitude and latitude grids of the target area, and the corresponding troposphere gas vertical column concentration of each target gas in each theodolite point on the equal longitude and latitude grids is obtained;
acquiring available pixels in the target area, and calculating a daily spatial average value of the tropospheric vertical column concentration of each target gas of all the available pixels;
based on CO, O 3 And the daily spatial average value data of the tropospheric vertical column concentration of HCHO, calculating the daily spatial average value of the tropospheric vertical column concentration of all available pixel HCHO secondary sources, wherein the HCHO secondary sources are the vertical column concentrations generated by photochemical reaction of tropospheric formaldehyde;
for HCHO secondary source, NO 2 Carrying out normalization treatment on the daily space average value of the troposphere vertical column concentration to obtain a corresponding standard value;
based on O 3 Daily spatial average of tropospheric vertical column concentration, daily spatial average of tropospheric vertical column concentration of HCHO secondary source and standard value and NO thereof 2 Daily spatial average value of tropospheric vertical column concentration and standard value thereof, obtaining a threshold value of a VOCs control zone and NO x A control region threshold;
based on the VOCs control region threshold and the NO x Control zone threshold, troposphere ozone control dividing each theodolite point in the target zoneTypes.
2. The ozone control type identification method according to claim 1, wherein the step of regularizing the pixels of each target gas on the equal longitude and latitude grid of the target area comprises the steps of:
linearly interpolating longitude and latitude of any two adjacent pixels on the satellite original pixel plane and vertical column concentration of each target gas to obtain a derivative pixel;
dividing the derivative pixel and the original pixel to corresponding theodolite points;
and averaging the vertical column concentration of each target gas pixel in each theodolite point to obtain the vertical column concentration value of each target gas on each theodolite point.
3. The ozone control type identification method according to claim 2, wherein the averaging of the vertical column concentrations of the target gas pixels in each theodolite is achieved by the following formula:
wherein V is i Represents the vertical column concentration of the tropospheric gas corresponding to the target gas on the ith pixel to be processed on the plane of the original pixel of the satellite after linear interpolation, lon i ，lat i Is the longitude and latitude of the ith pixel, v lon，lat For the concentration of the troposphere vertical column corresponding to the target gas in the corresponding theodolite points on the equal longitude and latitude grid, lon and lat are the longitude and latitude of the ith pixel, and grid is the degree of phase difference between the theodolite points.
4. The ozone control type identification method of claim 1, wherein the acquiring available pixels in the target area comprises:
and removing pixels with cloud quantity larger than 0.5 according to cloud quantity 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 according to claim 1, wherein the calculating of the daily spatial average value of the tropospheric vertical column concentration of each target gas for all available picture elements is accomplished by the following formula:
wherein V is x,i For the daily spatial average of the tropospheric vertical column concentration of the target gas, x represents O 3 Secondary source of HCHO or NO 2 I represents the ith day in the set time period, { region } is the set of all equal theodolite points in the target area, N is the number of the theodolite points in { region }, lon and lat are the longitude and latitude of a single theodolite point, v x,i,lon,lat The tropospheric gas vertical column concentration on day i in a single theodolite point for the target gas.
6. The method of claim 1, wherein the calculating of the daily spatial average of the tropospheric vertical column concentrations for all available pel HCHO secondary sources is accomplished by a polynomial fitting operation of:
V HCHO ＝(1(V CO ) 1 …(V CO ) n )A+(1(V O3 ) 1 …(V O3 ) n )B+e
wherein V is p,HCHO Daily spatial average of tropospheric vertical column concentration for HCHO secondary source, V HCHO Daily spatial average of total vertical column concentration of tropospheric formaldehyde, V CO 、V O3 CO and O respectively 3 Each of the tropospheric vertical column concentrations of (2)The daily spatial average value is calculated,v respectively CO 、V O3 Polynomial coefficients (n)<5) E represents the fit residual.
7. The ozone control type identification method according to claim 1, wherein the secondary source of HCHO, NO 2 The daily spatial average value of the tropospheric vertical column concentration is normalized by the following formula:
V x,i,nor ＝V x,i /V x,ref
V x,ref the method is obtained by the following formula:
wherein V is x,i Daily spatial average of tropospheric vertical column concentration for target gas x, x being HCHO secondary source or NO 2 ，V x,i,nor Is V (V) x,i Normalized standard value of V x,ref Daily reference value for tropospheric vertical column concentration of target gas x, d is total number of days of the study time period.
8. The ozone control type identification method according to claim 1, wherein the VOCs control zone threshold and the NO are obtained x The method for controlling the zone threshold comprises the following steps:
daily spatial average of tropospheric vertical column concentration of HCHO secondary source on the same day was compared with NO 2 The ratio of the daily spatial average of tropospheric vertical column concentrations is used as a marker;
selecting as a first data set all data satisfying said signature less than a first parameter variable reference value from a data element comprising O of said target area per day during said investigation time period 3 Daily spatial average of tropospheric vertical column concentration, pairs of HCHO secondary sourcesDaily spatial average value of vertical column concentration of flow layer and standard value and NO thereof 2 Daily spatial average of tropospheric vertical column concentration and standard value thereof;
linear fitting of O under the first dataset 3 A first slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration for a secondary source of HCHO, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A second slope of a standard value of a daily spatial average value of tropospheric vertical column concentration, obtaining a first variable reference value corresponding to when the first slope is equal to the second slope, and taking the first variable reference value at that time as the VOCs control zone threshold;
selecting all data meeting the condition that the mark is smaller than a second parameter variable reference value from the data elements as a second data set;
linear fitting of O under the second dataset 3 A third slope of the daily spatial average of the tropospheric vertical column concentration for a standard value of the daily spatial average of the tropospheric vertical column concentration of the HCHO secondary source, and O 3 Daily spatial average of tropospheric vertical column concentration with NO 2 A fourth slope of a standard value of a daily spatial average value of tropospheric vertical column concentration, obtaining a second variable reference value corresponding to when the third slope is equal to the fourth slope, and taking the second variable reference value at that time as the NO x A control region threshold.
9. The ozone control type identification method according to claim 1, wherein the control zone threshold and the NO are based on the VOCs x A control zone threshold, dividing the tropospheric ozone control type of each theodolite point in the target zone, comprising:
calculating the troposphere vertical column concentration and NO of the HCHO secondary source in each theodolite point in the target area 2 The ratio of tropospheric vertical column concentration;
dividing the ozone control type of the theodolite point with the ratio smaller than the threshold value of the VOCs control area into the VOCs control type;
will have a ratio greater than that of the NO x Ozone control type division of theodolite of control zone threshold into NO x A control type;
will have a ratio greater than the VOCs control region threshold and less than the NO x Ozone control types of theodolite points of the control zone threshold are classified as transition control types.
10. An electronic device, the electronic device comprising:
at least one processor; the method comprises the steps of,
a memory communicatively coupled to the at least one processor; wherein,
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 CN114496117A (en) | 2022-05-13 |
| CN114496117B true 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 |
|---|
| 应静瑶 ; 陈奔 ; 孙兴富 ; .温岭市臭氧浓度的控制区类型研究及其应用.能源环境保护.2018,(05),59-62. * |
| 陈雪萍 ; 巨天珍 ; 张江峪 ; 咸龙 ; 王培玉 ; 张生财 ; .兰州市大气臭氧与其前体物及气象因子相关性研究.环境科学研究.2019,(12),123-131. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114496117A (en) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pendergrass et al. | Surface and top-of-atmosphere radiative feedback kernels for CESM-CAM5 | |
| Jia et al. | Global inorganic nitrogen dry deposition inferred from ground-and space-based measurements | |
| Hu et al. | Evaluation of the new capture vaporizer for aerosol mass spectrometers (AMS) through field studies of inorganic species | |
| Takahama et al. | Quantification of carboxylic and carbonyl functional groups in organic aerosol infrared absorbance spectra | |
| Zhang et al. | Monitoring global tropospheric OH concentrations using satellite observations of atmospheric methane | |
| 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 | |
| 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 | |
| Baresel et al. | Comparison of nitrous oxide (N2O) emissions calculations at a Swedish wastewater treatment plant based on water concentrations versus off-gas concentrations | |
| CN113720972B (en) | Source analysis method and device for atmospheric pollutants | |
| Murphy et al. | Non-methane volatile organic compounds predict odor emitted from five tunnel ventilated broiler sheds | |
| Thorson et al. | Using a low-cost sensor array and machine learning techniques to detect complex pollutant mixtures and identify likely sources | |
| Latombe et al. | Comparison of spatial downscaling methods of general circulation model results to study climate variability during the Last Glacial Maximum | |
| Wei et al. | Development and evaluation of a robust temperature sensitive algorithm for long term NO2 gas sensor network data correction | |
| CN112948761A (en) | River nitrogen pollutant quantitative source analysis system | |
| Vikram et al. | Evaluating and improving the reliability of gas-phase sensor system calibrations across new locations for ambient measurements and personal exposure monitoring | |
| Fang et al. | Observation of atmospheric CO2 and CO at Shangri-La station: results from the only regional station located at southwestern China | |
| Furst et al. | A portable device for methane measurement using a low-cost semiconductor sensor: Development, calibration and environmental applications | |
| Zhang et al. | Impacts of climate change and emissions on atmospheric oxidized nitrogen deposition over East Asia | |
| CN114496117B (en) | Ozone control type identification method based on satellite hyperspectral remote sensing and electronic equipment | |
| Tomicic et al. | Experimental study of H 2 SO 4 aerosol nucleation at high ionization levels | |
| Russell et al. | Enhanced ambient sensing environment—a new method for calibrating low-cost gas sensors | |
| Wang et al. | Chemical Characteristics and Source-Specific Health Risks of the Volatile Organic Compounds in Urban Nanjing, China | |
| Taguem et al. | Quantile regression with a metal oxide sensors array for methane prediction over a municipal solid waste treatment plant | |
| CN115575337B (en) | High-precision atmospheric CO2 concentration observation data calibration method, system and equipment | |
| Ye et al. | HONO chemistry at a suburban site during the EXPLORE-YRD campaign in 2018: formation mechanisms and impacts on O 3 production |
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 |