CN112257551B - Method and system for identifying nitrogen oxide pollution source and determining emission - Google Patents

Method and system for identifying nitrogen oxide pollution source and determining emission Download PDF

Info

Publication number
CN112257551B
CN112257551B CN202011119845.4A CN202011119845A CN112257551B CN 112257551 B CN112257551 B CN 112257551B CN 202011119845 A CN202011119845 A CN 202011119845A CN 112257551 B CN112257551 B CN 112257551B
Authority
CN
China
Prior art keywords
emission
area
pollution source
remote sensing
acquiring
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
Application number
CN202011119845.4A
Other languages
Chinese (zh)
Other versions
CN112257551A (en
Inventor
王堃
吴锴
岳涛
童亚莉
高佳佳
马铭晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Beijing Municipal Institute of Labour Protection
Original Assignee
Ocean University of China
Beijing Municipal Institute of Labour Protection
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ocean University of China, Beijing Municipal Institute of Labour Protection filed Critical Ocean University of China
Priority to CN202011119845.4A priority Critical patent/CN112257551B/en
Publication of CN112257551A publication Critical patent/CN112257551A/en
Application granted granted Critical
Publication of CN112257551B publication Critical patent/CN112257551B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/13—Satellite images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
    • 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
    • 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

Abstract

The embodiment of the invention provides a method and a system for identifying a nitrogen oxide pollution source and determining emission, wherein the method comprises the following steps: acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area; acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data; acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area; fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides. The embodiment of the invention can identify NO2The key emission area is used for judging the functional area where the pollution source is positioned and carrying out NO treatmentxAnd estimating the emission, wherein the method is used for identifying the main nitrogen oxide emission source in the target area and calculating the emission.

Description

Method and system for identifying nitrogen oxide pollution source and determining emission
Technical Field
The invention relates to the technical field of environmental monitoring and pollution source identification, in particular to a method and a system for identifying a nitrogen oxide pollution source and determining emission.
Background
Nitrogen Oxides (NO)x) Is an important precursor for generating secondary aerosol and ozone in the atmosphere near the stratum in summer. In order to evaluate the emission of nitrogen oxides in a certain time span and a certain spatial region, the existing method mainly calculates the emission of nitrogen oxides by an emission factor method based on activity level, and estimates the emission of pollutants in different time periods by using a time profile. However, the same kind of sources adopt a uniform time profile, and it is difficult to reflect the time emission characteristic difference existing between the same kind of sources. Meanwhile, uncertainty in the activity level, the emission factor, and the like causes uncertainty in the emission amount.
NO2Column concentration data can be obtained through inversion of satellite sensors such as OMI and TROPOMI, and the spatial distribution characteristics of nitrogen oxides in a target area can be quantitatively reflected, so that ground monitoring sites are supplemented. As the smoke plume diffusion of the point source follows Gaussian distribution in space, the existing method is more for NO2And fitting the linear density of the column concentration on the space and further calculating the nitrogen oxide emission of the pollution source.
However, this approach is mostly aimed at large point sources located remotely or in urban areas that are less affected by other sources of pollution. For areas with complex and concentrated pollution sources, how to identify important pollution sources and carry out emission accounting is a key technical problem of carrying out important pollution source identification and emission accounting by applying remote sensing data.
Disclosure of Invention
The embodiment of the invention provides a method and a system for identifying a nitrogen oxide pollution source and determining emission, which are used for overcoming the defect that the identification of the pollution source and the emission accounting cannot be determined in the prior art and realizing the identification of an important pollution source and the emission accounting.
The embodiment of the invention provides a method for identifying a nitrogen oxide pollution source and determining emission, which comprises the following steps:
acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area;
acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data;
acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
According to the method for determining the emission amount of nitrogen oxides, the obtaining of the preprocessed remote sensing data corresponding to the target area specifically includes:
and preprocessing the remote sensing data corresponding to the target area by a space-time oversampling method to obtain a preprocessed high-spatial-resolution remote sensing image.
According to the method for determining the emission amount of nitrogen oxides, a heavily polluted area in the target area is obtained according to the preprocessed remote sensing data, and the method specifically comprises the following steps:
carrying out image graying on the preprocessed remote sensing data in sequence to obtain grayed image data;
denoising the grayed image data to obtain denoised image data;
and acquiring the heavily polluted area according to the image data subjected to denoising processing and a preset heavily polluted discharge area template.
According to an embodiment of the method for determining the emission amount of nitrogen oxides, the calculation formula of the preset heavy pollution emission area module is as follows:
Figure BDA0002731611740000031
Tx,y＝Wx,y-const，
wherein dstx,yTemplate representing said preset heavy pollutant emission area, Fx,yRepresenting de-noised image data, Tx,yAnd representing a threshold corresponding to each element in the image data after the de-noising processing, const representing an offset value adjustment quantity, x representing a line number, and y representing a column number.
According to the method for determining the emission amount of nitrogen oxides, the method for obtaining the heavy pollution area according to the image data after the denoising processing and the preset heavy pollution emission area template specifically comprises the following steps:
Outx,y＝dstx,y×Cx,y，
and Out represents the heavy pollution area, dst represents the preset heavy pollution discharge area template, and C represents the image data after denoising treatment.
According to the method for determining the emission amount of nitrogen oxides, the NO of each key pollution source and the periphery of each key pollution source in the main wind direction and the vertical wind direction is fitted2Column concentration line concentration, obtaining NO of each key pollution source2The molecular weight, the specific formula is as follows:
Figure BDA0002731611740000032
Figure BDA0002731611740000041
Figure BDA0002731611740000043
wherein A isx,yRepresents said NO2Molecular weight, BxIndicates NO in the X-axis direction2Intercept of line concentration fitting line, ByIndicating NO in the Y-axis direction2Intercept of line concentration fitting line, AxIndicates NO in the X-axis direction2Amplitude of the line concentration fitted line, AyIndicating NO in the Y-axis direction2Amplitude of the line concentration fitted line, fx，y(x, y) is hypothetical NO2Spatial distribution obeyed by column concentration, calculation Ax,yThe required parameters can all be obtained by fitting an exponential modified gaussian distribution.
According to one embodiment of the invention, the method for determining the amount of nitrogen oxides is based on NO2And the ratio of the nitrogen oxide to obtain the emission of the nitrogen oxide specifically comprises the following steps:
Figure BDA0002731611740000042
Figure BDA0002731611740000044
wherein E represents the emission of said nitrogen oxides, Ax,yRepresenting the molecular weight, tau is determined according to the wind speed and the e exponential decay distance of the downwind direction of the target area, alpha is the ratio of J to k, J is NO2First order rate constant of photolysis, k being NO and O3Reaction rate constant, O3Is the ozone concentration.
The embodiment of the invention also provides a system for identifying the nitrogen oxide pollution source and determining the emission, which comprises the following components:
the preprocessing module is used for acquiring preprocessed high-spatial-resolution remote sensing data corresponding to the target area;
the heavily-polluted region identification module is used for acquiring a heavily-polluted region in the target region according to the preprocessed remote sensing data;
the key pollution source identification module is used for acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
the emission calculation module is used for fitting each key pollution source and NO around the key pollution source in the main wind direction and the vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
The embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement any of the steps of the method for determining an emission amount of nitrogen oxide as described above.
Embodiments of the present invention also provide a non-transitory computer readable storage medium, on which a computer program is stored, the computer program, when being executed by a processor, implementing the steps of the method for determining an emission amount of nitrogen oxides as described in any one of the above.
The method and the system for identifying the nitrogen oxide pollution source and determining the emission can identify NO2And in the important emission area, judging the functional area where the pollution source is located and estimating the emission amount of NOx, wherein the method is used for identifying the main nitrogen oxide emission source in the target area and calculating the emission amount of the main nitrogen oxide.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
Fig. 1 is a flowchart of a method for identifying a nitrogen oxide pollution source and determining an emission amount according to an embodiment of the present invention;
FIG. 2 is a schematic view of identifying key pollution areas of nitrogen oxides in Beijing in the embodiment of the present invention;
fig. 3 is a schematic structural diagram of a system for identifying a nitrogen oxide pollution source and determining an emission amount according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Image segmentation techniques are techniques that divide an image into several specific regions with unique properties and propose an object of interest. NO of large point sources or important area sources2The column concentration has a higher value in space than the periphery, and the shape of the smoke plume can be identified and extracted by setting a threshold value. However, due to the different nitrogen oxide levels in different regions, it is difficult to spatially identify the dominant plume in a region if a global single threshold is set. Therefore, how to identify key point sources in different areas according to the background concentration of nitrogen oxides in each area is a key technical problem to be solved at present.
A point of interest (POI) is an important component of map big data, and may be a description of information of a parking lot, a commercial house, a factory enterprise, and the like. The main contents of the method comprise: the name, the category and the longitude and latitude are widely applied to various living fields such as navigation, and have the characteristics of detailed classification, regular updating, high accuracy, easiness in acquisition and the like. The data is widely applied to urban functional area identification research, so that the judgment of the area where the smoke plume is identified based on the interest point provides support for judging the area where the smoke plume is located and the type of a pollution source.
The invention aims to overcome the defects of the prior art and provide a quick and efficient method for identifying and quantifying a nitrogen oxide emission hot spot area. The method can be used for quickly identifying the main pollution source and the functional area of the target area and carrying out the NOx emission accounting of the pollution source.
Fig. 1 is a flowchart of a method for identifying a nitrogen oxide pollution source and determining an emission amount according to an embodiment of the present invention, as shown in fig. 1, the method includes:
s1, acquiring preprocessed high spatial resolution remote sensing data corresponding to the target area;
s2, acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data;
s3, acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
s4, fitting NO of each key pollution source and the periphery of each key pollution source in the main wind direction and the vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
The scheme adopted by the embodiment of the invention comprises remote sensing data preprocessing, hot spot area extraction, hot spot area functional area identification and emission quantitative estimation, and provides a method for quickly identifying an important nitrogen oxide emission source and quantitatively accounting pollutant emission for a target area.
Wherein the remote sensing data preprocessing improves the spatial resolution of the remote sensing data by a space-time oversampling method to reflect NO2A spatial distribution characteristic.
The hot spot area extraction is based on preprocessed remote sensing data, a target area is divided into a plurality of sub-areas, Gaussian window weighted sum of adjacent domain values is selected as a threshold value to generate an extraction template of the emission hot spot area, namely the remote sensing data is higher than the threshold value and is set to be 1, the remote sensing data is lower than the threshold value and is set to be 0, and the hot spot area is a heavily polluted area.
The functional area identification combines the hot spot area with the interest point, the traffic network data and the land utilization data to judge possible pollution sources of the area, the possible pollution sources are key pollution sources in the embodiment of the invention, the key pollution sources of different areas are different, for example, a certain area belongs to a traffic key area, and the possible pollution sources are road moving sources.
Quantitative estimation of emission by fitting the NO of pollution source and its periphery in the main wind direction and its vertical wind direction2Column concentration line concentration to determine NO at background concentration for different contaminant source removal2Molecular weight, and further use of NO2With NOxThe ratio of (2) calculates the amount of nitrogen oxide emissions.
The method for identifying the nitrogen oxide pollution source and determining the emission can identify NO2The key emission area is used for judging the functional area where the pollution source is positioned and carrying out NO treatmentxAnd estimating the emission, wherein the method is used for identifying the main nitrogen oxide emission source in the target area and calculating the emission.
On the basis of the foregoing embodiment, preferably, the acquiring the preprocessed remote sensing data corresponding to the target area specifically includes:
and preprocessing the remote sensing data corresponding to the target area by a space-time oversampling method to obtain a preprocessed high-spatial-resolution remote sensing image.
Using NO over a time span2Column concentration data (such as a month, a quarter and a year) is subjected to space-time oversampling on remote sensing data by utilizing a Gaussian kernel function, and the calculation process is shown as formulas (1) to (5):
Figure BDA0002731611740000081
Figure BDA0002731611740000091
Figure BDA0002731611740000092
Figure BDA0002731611740000093
Figure BDA0002731611740000094
wherein C (j) is the result of oversampling and Ω (i) is NO2Column concentration, S (i, j) is the intersection area of the grid researched by the satellite pixel points, i is the pixel point, j is the grid point, S (x, y) is a two-dimensional Gaussian kernel function, x and y are rectangular coordinates from a satellite track center store, wx、wy、k1、k2And k3Is a kernel function parameter.
On the basis of the foregoing embodiment, preferably, the acquiring a heavily polluted region in the target region according to the preprocessed remote sensing data specifically includes:
carrying out image graying on the preprocessed remote sensing data in sequence to obtain grayed image data;
denoising the grayed image data to obtain denoised image data;
and acquiring the heavily polluted area according to the image data subjected to denoising processing and a preset heavily polluted discharge area template.
Aiming at the preprocessed remote sensing data, extracting a nitrogen oxide hot spot region, namely a heavily polluted region, by adopting an adaptive threshold method based on Gaussian weighted sum:
first, image graying is performed.
Figure BDA0002731611740000095
Wherein, Cx,yIs the space-time oversampling result of formula (1), G is the grayed remote sensing data, x is the line number, y is the column number, CmaxThe value is the maximum value of the remote sensing data after oversampling, and 255 is the maximum value of the gray scale.
Then, image denoising processing is performed.
Fx,y＝filter{Gx-i,y-i},(i,j)∈S， (7)
Wherein, Fx,yThe remote sensing data after noise removal, G, x, y, S and filter are the remote sensing data after graying, x is a line number, y is a column number, S is a template window, and the filter is a filter.
Then, a heavy pollution discharge area template is preset.
Tx,y＝Wx,y-const， (8)
Figure BDA0002731611740000101
Wherein dst is an extraction template, T is a threshold value corresponding to data by data, W is a Gaussian weight sum of a certain neighborhood size, the calculation process is the same as that of the formula (7), and a Gaussian filter is adopted; const is the offset adjustment, x is the row number and y is the column number.
And finally, extracting the hot spot region.
Outx,y＝dstx,y×Cx,y， (10)
Out is the hot spot area after extraction; x is a row number; y is the column number.
The method for extracting the key pollution source of the heavily polluted area in the embodiment of the invention comprises the following steps:
and identifying the hot spot area functional area by using a local emission list and POI data.
And preferably combining the emission list to identify the hot spot area. Such as distinguishing the main pollution sources in the hot spot area extracted in a certain city based on the local emission list.
And judging by adopting POI and a road network for the condition that the non-emission list data basis is adopted or the non-emission list in the hot spot area relates to a point source. Firstly, judging whether industrial factory POI with high nitrogen oxide emission intensity exists in a hot spot area; secondly, judging whether the road network density is a main road motor vehicle emission area or not; again, it is determined whether there are emission sources within the area that are typical non-road movement sources for airports, ports, and the like.
In the embodiment of the invention, the NO of each key pollution source and the periphery of the key pollution source in the main wind direction and the vertical wind direction are fitted2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2The specific calculation process for obtaining the emission of nitrogen oxides based on the ratio of nitrogen oxides to nitrogen oxides is as follows:
the embodiment of the invention provides that an orthogonal coordinate system is set in the main wind direction and the vertical wind direction of the target source, and the source center is set as the origin of coordinates.
NO in two directions perpendicular to coordinate system by referring to Gaussian atmospheric pollutant diffusion model2The distribution of the column concentrations is independent of each other, and the combined distribution reflects NO2Two-dimensional spatial distribution characteristics of column concentration.
Based on two-dimensional NO2The linear concentration (equation (11)) and the amplitude after fitting of the above-mentioned joint distribution reflect the target source NO under the assumption of exponentially modified Gaussian distribution2The sum of the molecular weights is represented by the formula (14).
The embodiment of the invention respectively fits NO2 linear concentration D in the x-axis direction and the y-axis direction by using the exponential correction Gaussian distributionX(x) And DY(y) As shown in equation (12), the function f is obtainedX(x)、fY(y) and its amplitude Ax、AyStandard deviation σx、σyIntercept Bx、By。
Calculating A by taking the joint distribution of the two variables under independent conditions as the product of edge distributionx，yAs shown in formula (15). At the same time, f after fitting at the origin of coordinatesX(0)、fY(0) And fx，y(0, 0) both represent the value of the remote sensing data in the source center, and the embodiment of the invention proposes to adopt fX(0)、fY(0) Average value of (A) and (B)x、ByRespectively estimate fx，y(0，0)、Bx，yTo obtain Ax，yAs shown in formula (16). Finally, the embodiment of the present invention estimates NO by using the Janssen method2Is NOxIs given by equation (17), and further estimates the row using equation (19)And (4) putting.
DX(x)＝1×CX,Y＝0， (11)
Figure BDA0002731611740000121
Figure BDA0002731611740000122
Figure BDA0002731611740000129
Figure BDA0002731611740000124
Figure BDA0002731611740000125
Figure BDA0002731611740000126
Figure BDA0002731611740000127
Figure BDA0002731611740000128
Wherein D isNO2，xNO perpendicular to the main wind direction2Linear density, mole/cm, Ax,yIs molecular weight, E is NOxDischarge amount, ton/hour, t0Is the downwind e exponential decay distance, BxThe method is a baseline value, mostly adopts asymmetric least square regression fitting or a constant value, erf is a Gaussian error function, l is the spatial resolution of the oversampled remote sensing data, and km and 46 are NOx molesR is NO2The atmospheric lifetime is small, typically 5 hours, depending on the local wind speed and t0An estimation is performed.
Figure BDA0002731611740000131
Is NO in steady state for photochemical reaction2Is NOxThe ratio and the relevant parameters can be referred to Janssen method-related research, J is NO2The first order rate constant of photolysis, referenced to, can be 2.2X 10-3s-1K is NO and O3The reaction rate constant, referred to, can be 3.7X 10-4ppb-1s-1(ii) a Alpha is the ratio of J to k, the reference can take 5.95ppb, u is the wind speed, O3Is the ozone concentration, ppb.
The practical applicability of the method is verified in the area of Jingjin Ji and the typical source thereof.
1. And extracting and identifying the hot spot area.
No in Jingjin Ji area after space-time oversampling by the method of the embodiment of the invention2The column concentration may reflect NO2After the hot spot region is extracted by the method of the embodiment of the invention, the main hot spot region in space can be extracted.
Taking Beijing as an example, the method of the embodiment of the invention is used for identifying NO in Beijing2The hot spot areas mainly comprise three areas, fig. 2 is a schematic diagram of identifying key nitrogen oxide polluted areas in Beijing City in the embodiment of the invention, and as shown in fig. 2, the areas in five rings and Chang Ping mainly comprise road motor vehicles and the areas in the first International airport mainly comprise non-road moving sources.
3、NOxAnd (5) checking and calculating the discharge amount.
The embodiment of the invention takes the capital airport and the canadian county industrial park in Beijing city as an example, estimates the emission by adopting summer remote sensing data, and carries out the result comparison of the method and the local emission list. It should be noted that the lag of the emission list in time and the use of only summer remote sensing data in this case are one of the main reasons for the difference of the comparison results.
(1) Beijing capital International airport:
as shown in fig. 2, the hot spot area of the capital airport is calculated by using equations (15) to (18) because the influence of the capital international airport distance on the five-ring intra-urban area is obvious.
In the embodiment of the invention, an estimation orthogonal coordinate system is established in the main wind direction of the extracted capital international airport remote sensing data, and NO in two directions is respectively measured2The column concentrations were fitted and NO was calculated from equations (15) to (18)xThe hourly discharge is about 1.07 t/h, and the annual discharge is estimated to be about 9371 tons.
According to research on water resources and the like, the NO of various mobile source years of international airports in the capital of 2013 is foundxThe discharge amount is about 6287 tons, which is close to the result obtained by the embodiment of the invention.
(2) The Caofen Dian industrial park.
The NOx emission amount of the Caofen Dian industrial park in 2017 is estimated from the existing local atmospheric pollutant emission list data to be about 1.4 ten thousand tons. However, the first steel Jingtang steel is taken as an important pollution source of the industrial park, the ultra-low emission modification acceptance is completed in 2019, and the emission amount of pollutants is greatly reduced.
According to the embodiment of the invention, the emission is calculated by adopting the formulas (15) to (18), an estimation orthogonal coordinate system is established in the main wind direction of the extracted remote sensing data, and NO in two directions is respectively calculated2The column concentrations were fitted and the NOx hourly emissions were calculated according to equations (15) to (18) to be about 0.28 t/hr, estimated to be about 2453 tons per year.
(3) Wuhan steel and its periphery.
In addition, the embodiment of the invention estimates Wuhan steel and peripheral NO based on the formulas (15) to (18)xThe annual emission amount is 13640 tons after 1.56 tons/hour in hour, and the consistency is better than that of the local statistical emission amount of about 1.7 ten thousand tons.
Fig. 3 is a schematic structural diagram of a system for identifying a nitrogen oxide pollution source and determining an emission amount according to an embodiment of the present invention, as shown in fig. 3, the system includes a preprocessing module 301, a heavily polluted area identifying module 302, a key pollution source identifying module 303, and an emission amount calculating module 304, where:
the preprocessing module 301 is configured to obtain preprocessed high spatial resolution remote sensing data corresponding to the target area;
the heavily-polluted region identification module 302 is used for acquiring a heavily-polluted region in the target region according to the preprocessed remote sensing data;
the key pollution source identification module 303 is configured to obtain a key pollution source of the heavily polluted region according to a pollution source emission type in the heavily polluted region;
the emission calculation module 304 is used for fitting each key pollution source and the NO around the key pollution source in the main wind direction and the vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
The present embodiment is a system embodiment corresponding to the above method, and please refer to the above method embodiment for details, which is not described herein again.
Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 4, the electronic device may include: a processor (processor)410, a communication Interface 420, a memory (memory)430 and a communication bus 440, wherein the processor 410, the communication Interface 420 and the memory 430 are communicated with each other via the communication bus 440. The processor 410 may invoke logic instructions in the memory 430 to perform a method for identifying a source of nox pollution and determining an amount of emissions, the method comprising:
acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area;
acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data;
acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
In addition, the logic instructions in the memory 430 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, an embodiment of the present invention further provides a computer program product, where the computer program product includes a computer program stored on a non-transitory computer-readable storage medium, the computer program includes program instructions, and when the program instructions are executed by a computer, the computer can execute a method for identifying a nox pollution source and determining an emission amount, where the method includes:
acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area;
acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data;
acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
In still another aspect, an embodiment of the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored, where the computer program is implemented by a processor to execute the method for identifying and determining the emission amount of the nox pollution source provided in the foregoing embodiments, and the method includes:
acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area;
acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data;
acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (6)

1. A method for identifying a nitrogen oxide emission source and determining an emission amount is characterized by comprising the following steps:
acquiring preprocessed high-spatial-resolution remote sensing data corresponding to a target area;
acquiring a heavily polluted area in the target area according to the preprocessed remote sensing data; the acquiring of the heavily polluted area in the target area according to the preprocessed remote sensing data specifically comprises: carrying out image graying on the preprocessed remote sensing data in sequence to obtain grayed image data; denoising the grayed image data to obtain denoised image data; acquiring the heavily polluted area according to the image data subjected to denoising processing and a preset heavily polluted discharge area template; the calculation formula of the preset heavy pollution emission area module is as follows:
Figure FDA0003192587720000011
Tx,y＝Wx,y-const，
wherein dstx,yTemplate representing said preset heavy pollutant emission area, Fx,yRepresenting de-noised image data, Tx,yRepresenting the threshold corresponding to each element in the image data after de-noising, const representing the adjustment amount of the offset value, and x representing the line numberY represents a column number;
acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
fitting NO of each key pollution source and periphery thereof in main wind direction and vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
2. The method for determining the emission of nitrogen oxides according to claim 1, wherein the obtaining of the preprocessed high spatial resolution remote sensing data corresponding to the target area specifically comprises:
and preprocessing the remote sensing data corresponding to the target area by a space-time oversampling method to obtain a preprocessed remote sensing image.
3. The method for determining the emission amount of nitrogen oxides according to claim 1, wherein the obtaining the heavily polluted region according to the image data after the denoising process and a preset heavily polluted region template specifically comprises:
Outx,y＝dstx,y×Cx,y，
wherein, Outx，yRepresents the hot spot area, dst after extractionx，yA template representing the preset heavy pollution discharge area, Cx，yRepresenting spatio-temporal oversampled processed image data.
4. A system for identifying a source of nitrogen oxide pollution and determining an amount of emission, comprising:
the preprocessing module is used for acquiring preprocessed high-spatial-resolution remote sensing data corresponding to the target area;
the heavily-polluted region identification module is used for acquiring a heavily-polluted region in the target region according to the preprocessed remote sensing data; the acquiring of the heavily polluted area in the target area according to the preprocessed remote sensing data specifically comprises: carrying out image graying on the preprocessed remote sensing data in sequence to obtain grayed image data; denoising the grayed image data to obtain denoised image data; acquiring the heavily polluted area according to the image data subjected to denoising processing and a preset heavily polluted discharge area template; the calculation formula of the preset heavy pollution emission area module is as follows:
Figure FDA0003192587720000021
Tx,y＝Wx,y-const，
wherein dstx,yTemplate representing said preset heavy pollutant emission area, Fx,yRepresenting de-noised image data, Tx,yRepresenting a threshold corresponding to each element in the image data after denoising, const representing an offset value adjustment quantity, x representing a line number, and y representing a column number;
the key pollution source identification module is used for acquiring a key pollution source of the heavily polluted area according to the emission type of the pollution source in the heavily polluted area;
the emission calculation module is used for fitting each key pollution source and NO around the key pollution source in the main wind direction and the vertical wind direction2Column concentration line concentration, obtaining NO of each key pollution source2Molecular weight, and according to NO2And obtaining the emission of the nitrogen oxides according to the proportion of the nitrogen oxides.
5. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method for determining an amount of emission of nitrogen oxides as claimed in any one of claims 1 to 3.
6. A non-transitory computer readable storage medium having a computer program stored thereon, wherein the computer program when executed by a processor implements the steps of the method for determining nox emissions according to any of claims 1 to 3.
CN202011119845.4A 2020-10-19 2020-10-19 Method and system for identifying nitrogen oxide pollution source and determining emission Active CN112257551B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011119845.4A CN112257551B (en) 2020-10-19 2020-10-19 Method and system for identifying nitrogen oxide pollution source and determining emission

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011119845.4A CN112257551B (en) 2020-10-19 2020-10-19 Method and system for identifying nitrogen oxide pollution source and determining emission

Publications (2)

Publication Number Publication Date
CN112257551A CN112257551A (en) 2021-01-22
CN112257551B true CN112257551B (en) 2021-09-10

Family

ID=74245009

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011119845.4A Active CN112257551B (en) 2020-10-19 2020-10-19 Method and system for identifying nitrogen oxide pollution source and determining emission

Country Status (1)

Country Link
CN (1) CN112257551B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113034363B (en) * 2021-03-11 2024-04-09 云南电网有限责任公司电力科学研究院 GEE depth space-time experience Kringing region scale-based nitrogen oxide rapid reduction method
CN113435291B (en) * 2021-06-22 2022-10-28 生态环境部华南环境科学研究所 Sewage outlet investigation workload accounting method, system, device and storage medium
CN114218976B (en) * 2021-09-27 2024-03-08 中国环境科学研究院 Method and system for tracing and positioning refined emission source in high-value area of ozone precursor
CN114460232B (en) * 2022-01-28 2023-02-28 生态环境部卫星环境应用中心 Multipoint source carbon-nitrogen ratio monitoring method and system based on vehicle-mounted measurement system
CN115082546B (en) * 2022-06-22 2023-03-21 中科三清科技有限公司 Method and device for determining pollutant discharge amount, electronic equipment and medium
CN114839327B (en) * 2022-07-01 2023-01-17 北京英视睿达科技股份有限公司 Emission source CO of satellite monitoring target area 2 Hourly discharge method and system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1337608A1 (en) * 2000-09-28 2003-08-27 Clean Diesel Technologies Inc. Low-emissions diesel fuel emulsions
US7071002B1 (en) * 2000-05-09 2006-07-04 Gordon-Darby Systems, Inc. Method and system for vehicle emission testing
CN104200089A (en) * 2014-08-27 2014-12-10 北京市环境保护监测中心 Method for measuring air pollutant emission amount of coal burned at bungalows
CN104865353A (en) * 2015-06-01 2015-08-26 上海交通大学 Atmospheric pollution data acquisition method for industrial park based on unmanned aerial vehicle
CN106153510A (en) * 2015-03-25 2016-11-23 日本电气株式会社 The evaluation method of pollutant and device
CN106845371A (en) * 2016-12-31 2017-06-13 中国科学技术大学 A kind of city road network automotive emission remote sensing monitoring system
CN109523083A (en) * 2018-11-21 2019-03-26 北京市环境保护监测中心 The foundation of air pollution emission Calculating model, measuring method and model foundation device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106650825B (en) * 2016-12-31 2020-05-12 中国科学技术大学 Motor vehicle exhaust emission data fusion system
CN110659838A (en) * 2019-09-27 2020-01-07 佛山科学技术学院 Ecological environment quality assessment method based on remote sensing and GIS

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7071002B1 (en) * 2000-05-09 2006-07-04 Gordon-Darby Systems, Inc. Method and system for vehicle emission testing
EP1337608A1 (en) * 2000-09-28 2003-08-27 Clean Diesel Technologies Inc. Low-emissions diesel fuel emulsions
CN104200089A (en) * 2014-08-27 2014-12-10 北京市环境保护监测中心 Method for measuring air pollutant emission amount of coal burned at bungalows
CN106153510A (en) * 2015-03-25 2016-11-23 日本电气株式会社 The evaluation method of pollutant and device
CN104865353A (en) * 2015-06-01 2015-08-26 上海交通大学 Atmospheric pollution data acquisition method for industrial park based on unmanned aerial vehicle
CN106845371A (en) * 2016-12-31 2017-06-13 中国科学技术大学 A kind of city road network automotive emission remote sensing monitoring system
CN109523083A (en) * 2018-11-21 2019-03-26 北京市环境保护监测中心 The foundation of air pollution emission Calculating model, measuring method and model foundation device

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Remotely sensed canopy nitrogen correlates with nitrous oxide emissions in a lowland tropical rainforest;Soper,FM 等;《ECOLOGY》;20180930;第99卷(第9期);全文 *
刘菲.基于卫星遥感的中国典型人为源氮氧化物排放研究.《中国博士学位论文全文数据库 工程科技Ⅰ辑》.2017,(第3期),第B027-38页. *
基于OMI遥感数据的河北省大气NO2垂直柱浓度时空变化研究;辛名威;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20140131(第1期);全文 *
基于卫星遥感的中国典型人为源氮氧化物排放研究;刘菲;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20170331(第3期);第B027-38页 *
太赫兹成像质量提升方法;张智 等;《红外与激光工程》;20171130;第46卷(第11期);第1-5页 *

Also Published As

Publication number Publication date
CN112257551A (en) 2021-01-22

Similar Documents

Publication Publication Date Title
CN112257551B (en) Method and system for identifying nitrogen oxide pollution source and determining emission
CN108648127B (en) Method for locking urban air pollution hot spot area
CN112182064B (en) Pollutant source analysis method and device, electronic equipment and storage medium
CN109360421B (en) Traffic information prediction method and device based on machine learning and electronic terminal
CN111751789B (en) Method, system, medium and equipment for forecasting radar detection range of artificial satellite
CN112991132B (en) PM2.5Method and device for identifying area cooperatively controlled by ozone
CN110738354B (en) Method and device for predicting particulate matter concentration, storage medium and electronic equipment
CN113538072A (en) Intelligent travel chain identification method and device for freight vehicle and electronic equipment
Kumar et al. An algorithm for automated estimation of road roughness from mobile laser scanning data
Ramezani A note on the normalized definition of Shannon’s diversity index in landscape pattern analysis
CN107909192B (en) Estimation method and device for heavy metal content in soil
CN111612064B (en) PM2.5Method and device for tracing pollution air mass, electronic equipment and storage medium
CN114218976B (en) Method and system for tracing and positioning refined emission source in high-value area of ozone precursor
CN114881381B (en) Urban ponding water level prediction method and system based on improved convolutional neural network
CN112905724B (en) Method and system for searching suspicious digestion site of muck truck
CN111209533B (en) Road surface fine particulate matter inversion method and device, computing equipment and storage medium
CN114002381A (en) Atmospheric pollution source-tracing diffusion analysis method and device
CN115792919B (en) Method for identifying polluted hot spot area through horizontal scanning monitoring of aerosol laser radar
CN115239027B (en) Method and device for forecasting air quality check set
CN111611335A (en) Method, system and storage medium for evaluating applicability of homeland space
CN104200472A (en) Non-local wavelet information based remote sensing image change detection method
CN114414446B (en) Near-ground particulate matter concentration estimation method, device and equipment based on satellite remote sensing
CN110057981A (en) The recognition methods of road network discharge of pollutant sources and system
CN115834838A (en) Method, device and medium for monitoring in tunnel
CN111461163B (en) Urban interior PM2.5 concentration simulation and population exposure evaluation method and device

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