CN114757807B - Multi-mode fused online accounting method for actual emission of atmospheric pollutants - Google Patents
Multi-mode fused online accounting method for actual emission of atmospheric pollutants Download PDFInfo
- Publication number
- CN114757807B CN114757807B CN202210660280.3A CN202210660280A CN114757807B CN 114757807 B CN114757807 B CN 114757807B CN 202210660280 A CN202210660280 A CN 202210660280A CN 114757807 B CN114757807 B CN 114757807B
- Authority
- CN
- China
- Prior art keywords
- park
- pollution
- source
- mode
- monitoring
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 33
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 110
- 238000009792 diffusion process Methods 0.000 claims abstract description 57
- 238000005457 optimization Methods 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 230000004927 fusion Effects 0.000 claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 68
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 150000001875 compounds Chemical class 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract 1
- 230000006872 improvement Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 239000008186 active pharmaceutical agent Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an online accounting method for actual emission of atmospheric pollutants by multi-mode fusion, which comprises the steps of constructing a multi-mode atmospheric pollution diffusion model to obtain a visual contribution matrix of pollution source intensity and a monitoring station in a current park; constructing a source parameter inversion algorithm, and carrying out effectiveness evaluation on the constructed multi-mode atmospheric pollution diffusion model; and calculating the total pollutant emission amount of the park based on the pollution source intensity obtained through optimization estimation. The invention integrates the advantages of multiple models through the fusion of the multiple models and the dynamic distribution of the model weight, and improves the accuracy, universality and stability of the set model.
Description
Technical Field
The invention relates to the technical field of environmental monitoring and environmental protection, in particular to an online accounting method for actual emission of atmospheric pollutants by multi-mode fusion.
Background
The industrial park is a centralized place of industrial production, plays a crucial role in the development of local economy and is also a large household for pollutant emission.
In the ecological civilization construction process, the key point for stably improving the local air environment level is to control the emission intensity and the total amount of atmospheric pollutants in an industrial park.
Accurate accounting of pollutant emission in an industrial park is the basis for formulating reasonable emission indexes and scientifically formulating emission reduction plans for the park.
At present, although most organized sewage outlets of major enterprises in industrial parks are provided with online monitoring instruments, the discharge forms such as organized discharge and the like are difficult to monitor through online monitoring means.
In view of the above, it is urgently needed to design and develop a set of scientific method and accounting system for calculating the total discharge amount of the park through the local air quality change monitoring data of the park.
Disclosure of Invention
The invention aims to provide an on-line accounting method for the actual emission of atmospheric pollutants by multi-mode fusion, which is characterized in that the method for directly calculating the total emission of a park through local air quality change monitoring data of the park is realized by acquiring comprehensive meteorological data information, park pollution sources, monitoring sites and topographic building information of the park, calculating the contribution score of the pollution sources to the monitoring sites and calculating the total emission of the park in a certain time, and solving the problems in the prior art. In order to achieve the purpose, the invention provides the following technical scheme:
the online accounting method for the actual emission of the atmospheric pollutants based on multi-mode fusion comprises the following steps:
rolling simulation is carried out to obtain meteorological environment data in a garden area; constructing a simulation grid comprising park information, wherein the park information comprises park monitoring station basic information, park pollution source basic information and park elevation information;
constructing a multi-mode atmospheric pollution diffusion model, simulating and acquiring the concentration distribution relation between the source intensity of the pollution source of the park and the position of a park monitoring station under the diffusion and atmospheric transmission effects according to meteorological environment data and real-time data of the park monitoring station, and acquiring a contribution matrix of the source intensity of the pollution source in the current park to the monitoring station so as to obtain the influence of the emission of the pollution source on the monitoring station;
constructing a source parameter inversion algorithm, and carrying out optimization estimation on the source intensity of the pollution source according to the acquired data of the park monitoring station and the simulation result output by the multi-mode atmospheric pollution diffusion model so as to carry out effectiveness evaluation on the constructed multi-mode atmospheric pollution diffusion model;
and calculating the total pollutant emission amount of the park based on the pollution source intensity obtained through optimization estimation.
As an improvement of the multi-mode fused online accounting method for the actual emission of the atmospheric pollutants, the meteorological environment data in the park area are acquired and comprise meteorological field data which can be directly extracted from a simulation result and meteorological field data which cannot be directly calculated; the simulation grids include a weather simulation grid and a campus simulation grid, wherein,
before acquiring the data of the park meteorological environment, the position index of the park position in the meteorological simulation grid needs to be calculated firstly:
in the formula, i represents a grid index corresponding to a pollution source in a meteorological simulation grid; j represents a grid index corresponding to a monitoring station in the meteorological simulation grid; x and y are respectively expressed as longitude and latitude information of the park or UTM coordinates x and y;
and
respectively representing longitude and latitude or UTM coordinates of a starting point of the meteorological simulation grid;
represents rounding down; d is expressed as the grid resolution.
As an improvement of the multi-mode integrated online accounting method for the actual emission of the atmospheric pollutants, before simulating the concentration distribution relationship between the source intensity of the pollution source of the park and the position of the park monitoring station, the acquired data of the park monitoring station needs to be subjected to unit conversion, and the validity of the acquired data of the park monitoring station is verified, and the method specifically comprises the following steps:
in the formula (I), the compound is shown in the specification,
and
respectively representing the monitoring concentration after unit conversion and before unit conversion, wherein MIN _ NSITE and MIN _ RSITE are the minimum effective monitoring site number and the minimum effective monitoring site proportion set for the configuration file, and Vlag represents whether the unit is effective or not;
is a conversion factor;
the calculation formula of (2) is as follows:
wherein P is a standard atmospheric pressure (Pa); r is an ideal gas constant with the unit of J/mol.K; t is the ambient temperature (K); w m Relative molecular mass (g/mol);
for the final mass unit after conversion to gramCalculating the proportion;
unit proportion before conversion; e.g. from ppm to mg/m 3 When the utility model is used, the water is discharged,
value of 10 3 ，
Value of 10 6 ；
Identifying the data validity of the ith monitoring station; and N is the total number of the monitored stations.
As an improvement of the multi-mode integrated online accounting method for the actual emission of the atmospheric pollutants, the specific method for constructing the park simulation grid is as follows:
in the formula, nx is the number of grids of the park simulation grid in the x direction,
and
maximum and minimum UTM-X coordinates of a pollution source and a monitoring point of the park respectively; offset is the amount of expansion of the grid, and is usually equal to the grid resolution;
representing a UTM area in which the latitude and longitude information is located; lng denotes longitude.
As an improvement of the multi-mode fused online accounting method for the actual emission of the atmospheric pollutants, the specific method for constructing the multi-mode atmospheric pollution diffusion model comprises the following steps:
s1, customizing a linear superposition hypothesis that the data of a single monitoring station of the garden are strong relative to pollution sources in a plurality of gardens to reduce the calculation difficulty:
in the formula (I), the compound is shown in the specification,
is the monitoring result of the monitoring station j, b is the background concentration,
in order to be the emission intensity of the pollution source i,
representing the contribution of the pollution source i to the monitored site j;
and S11, establishing the data of the N monitoring stations of the garden and an M x N x K three-dimensional contribution matrix of the M pollution sources based on a linear superposition hypothesis, wherein K represents the number of the multi-mode atmospheric pollution diffusion models participating in calculation so as to reduce the repeated calling of the multi-mode atmospheric pollution diffusion models in the optimization process, for example, for a Sutton model:
wherein u is the average wind speed; z and h are the height of the receptor site and the height of the source respectively;
and
is a diffusion coefficient related to meteorological conditions, and the brightness variation of the M N K three-dimensional contribution matrix shows the contribution degree phaseFor the difference in size.
As an improvement of the online accounting method for the actual emission of the atmospheric pollutants based on the multi-mode fusion, the method for optimally estimating the source intensity of the pollution source based on the source parameter inversion algorithm comprises the following steps:
s2, estimating the environmental background concentration of the park according to the pollution source intensity and the distribution condition of the monitoring stations and by combining the current wind field data, wherein the specific estimation steps are as follows:
vector pointing to any pollution source i from monitoring site j
Unit vector with real-time wind direction
Performing dot product calculation to obtain vector
In that
Projection in the direction;
when the concentration of the contaminant in the sample is low for any contamination source i, if,
if the monitoring station is a windward station, the monitoring station is considered as an upwind station;
when a plurality of upwind direction sites appear, averaging corresponding monitoring data to be used as the environmental background concentration of the park;
s21, calculating the pseudo inverse of the contribution matrix, obtaining an initial guess of the source intensity of the pollution source, optimizing the distribution of the source intensity of the pollution source to obtain a reliable initial guess, and optimizing the strong distribution of the source intensity of the pollution source, wherein the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is a row vector of 1 x M,
is a row vector of 1 x N,
is that
Pseudo-inverse of the matrix;
is an initial guess of the source strength of the pollution source;
s22, calculating the average relative deviation between the pollution source intensity and the measured value, taking the average relative deviation as the accuracy measurement of the multi-mode atmospheric pollution diffusion model calculation, dynamically generating the weight of the integration of the multi-mode atmospheric pollution diffusion model to avoid the subjective interference caused by artificially setting the integration weight of the multi-mode atmospheric pollution diffusion model, and the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is the average relative deviation of the multi-mode atmospheric pollution diffusion model K,
is a configurable constant.
As an improvement of the online accounting method for the actual emission of the atmospheric pollutants fused in multiple modes, the specific steps for calculating the total emission of the pollution in the park are as follows:
s3, marking the attribute of the pollution source to eliminate the influence of the pollution source outside the garden:
in the formula, Q represents the total amount of the pollutant source discharged outside the park.
As an improvement of the online accounting method for the actual emission of the atmospheric pollutants fused in multiple modes, when the total emission of the pollution in the park is calculated, the maximum total emission Q _ MAX and the maximum relative deviation RME _ MAX are set through a configuration file to verify the validity of the accounting result of the multiple-mode atmospheric pollution diffusion model:
when the fitting result of the multi-mode atmospheric pollution diffusion model is consistent with the park monitoring data, the simulation accounting of the multi-mode atmospheric pollution diffusion model is proved to be effective.
Compared with the prior art, the invention has the following beneficial effects:
1. according to the method, the advantages of multiple models are integrated through fusion of the multiple models and dynamic distribution of model weights, and the accuracy, universality and stability of the set model are improved;
2. by adding the algorithm design of the linear superposition hypothesis and the contribution matrix, the calling of an atmosphere pollution diffusion model with large calculation amount is reduced, the optimization process is simplified, and the timeliness of the algorithm is improved;
3. according to the invention, through setting the virtual source which does not participate in statistics, the influence of factors such as irregular park range, pollution sources outside the park boundary and the like on accounting is reduced.
Drawings
FIG. 1 is a flow diagram of a multi-mode fusion accounting framework proposed in an embodiment of the present invention;
FIG. 2 is a diagram of a campus information configuration interface in an embodiment of the present invention;
FIG. 3 is a plot of the simulated pollutant diffusion contour (left) and contribution matrix (right) for a campus as set forth in an embodiment of the present invention;
fig. 4 is a graph showing the relative magnitude of the emission intensity of the pollution source obtained by the simulation of the park proposed in one embodiment of the present invention (upper) and the comparison of the simulation result with the monitored value (lower).
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments, and all other embodiments obtained by a person of ordinary skill in the art without creative efforts based on the embodiments of the present invention belong to the protection scope of the present invention.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "front", "rear", "left", "right", "top", "bottom", "inner", "outer", and the like, indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are merely for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention.
The present invention will be described in further detail below with reference to the accompanying drawings, but the present invention is not limited thereto.
Referring to fig. 1-4, as an embodiment of the present invention, a multi-mode integrated online accounting method for actual emission of atmospheric pollutants includes a weather simulation module, a park information configuration module, a base station data collection module, a multi-mode atmospheric pollution diffusion model module, a source parameter inversion algorithm module and a total amount accounting module, it should be noted that,
the meteorological simulation module is used for rolling and simulating and calculating the atmospheric physical environment characteristics of the campus area according to the initial field data of the campus; the park information configuration module is used for converting longitude and latitude coordinates of park pollution sources and monitoring sites into local coordinates, processing the park terrain, and dividing to form a park simulation grid, wherein the park simulation grid comprises a meteorological simulation grid and a park simulation grid; a base station data aggregation module: collecting and verifying the effectiveness of the monitoring data of the monitoring station; the multi-mode atmospheric pollution diffusion model module is used for constructing a multi-mode atmospheric pollution diffusion model to simulate the influence of pollution source emission on a monitoring station; the source parameter inversion algorithm module is used for optimizing and estimating to obtain the source intensity information of the pollution source of the industrial park by combining monitoring data of the monitoring station and a simulation result of the multi-mode atmospheric pollution diffusion model module; a total amount accounting module: and calculating the total amount of pollution emission in the garden area in the time period based on the pollution source intensity obtained through optimization estimation.
The method comprises the following steps:
firstly, on the basis of a meteorological simulation module, responding to an API request and returning meteorological environment data of any park place in a simulation range to establish a meteorological simulation grid, wherein the acquired meteorological environment data comprise meteorological field data which can be directly extracted from a simulation result, such as temperature, air pressure, wind direction, wind speed and the like; and data which needs to be further calculated and obtained, such as sunshine intensity, rainfall amount per unit time and the like, wherein the sunshine intensity netRad and the accumulated rainfall raining amount rainCum are calculated in the following way:
in the formula, SW and LW respectively represent short wave and long wave radiation, RAINC and RAINSH respectively represent accumulative precipitation, accumulative grid precipitation and accumulative shallow precipitation, and it needs to be explained that the meteorological simulation module can adopt WRF mode to simulate campus meteorology, especially high altitude meteorology, when in specific implementation, the module can be independently arranged in a high performance Linux server, the calculation efficiency is improved by parallel calculation, and the day-wide meteorological simulation of the same day is completed in a plurality of hours in the morning to support the operation of the total calculation model;
when a user acquires meteorological environment data based on a meteorological simulation module, a form is made according to the specifications, a park information configuration module reads in the form, a local coordinate system is automatically established, longitude and latitude information (the longitude and latitude information of a pollution source and the longitude and latitude information of a monitoring station) is converted into local coordinates, a simulation grid is generated according to grid resolution input by the user, visually displaying on a map, exporting the json file configured in the park by one key during actual use, copying the json file configured in the park to an accounting system to make a directory by a user, and modifying the configuration file of the accounting system, after the park name, the ID and the configuration file name are indicated, the park information configuration can be added into the real-time accounting system, and the park information comprises park monitoring station basic information, park pollution source basic information and park elevation information; the basic information of the park monitoring station comprises the longitude and latitude and the height of the monitoring station, the data of the park monitoring station is the concentration value of the atmospheric pollutants obtained by a monitoring instrument of the monitoring station, and the basic information of the park pollution source comprises the data of the position, the height, the temperature of a discharge port and the like of the pollution source; the campus elevation information refers to altitude information of each location of the campus.
In an embodiment of the present invention, it should be noted that, when the weather simulation module responds to the API request and returns the weather environment data of any park location within the simulation range, the weather simulation module needs to convert the longitude and latitude information (the longitude and latitude information of the pollution source and the longitude and latitude information of the monitoring station) in the API into the UTM coordinate when constructing the model weather simulation grid
Then, calculating the grid index in the meteorological simulation grid corresponding to the park information (park position):
in the formula, i is expressed as a grid index corresponding to a pollution source in the meteorological simulation grid; j represents a grid index corresponding to a monitoring station in the meteorological simulation grid; x and y are respectively expressed as longitude and latitude information of the park or UTM coordinates x and y;
and
respectively representing longitude and latitude or UTM coordinates of a starting point of the meteorological simulation grid;
represents rounding down; d is the grid resolution, and it can be understood that the value "0.5" in the above formula is the offset of the UTM coordinate x, y, and the purpose is to accurately calculate; it can be understood that the method can decouple the meteorological simulation data with large calculation amount and long time consumption from the total accounting data with small calculation amount and high real-time requirement on software and hardware deployment, and the mode of carrying out meteorological simulation calculation by using the high-performance server in the period from sub-night to early morning makes the required total accounting part directly obtain the required meteorological data from the high-performance server through the API, thereby reducing the calculation amount.
Meanwhile, after acquiring the grid index in the meteorological simulation grid corresponding to the park information (park position), the park information configuration module automatically generates the park simulation grid according to the park pollution source position and the monitoring station position and synchronously acquires the geographical elevation information:
in the formula, nx is the number of grids of the park simulation grid in the x direction,
and
maximum and minimum UTM-X coordinates of a pollution source and a monitoring point of the park respectively; offset is the amount of expansion of the grid, usually equal to the grid resolution;
representing the UTM area in which the longitude and latitude information is located; r and c are indexes for searching the geographical elevation file of the park; long and lat respectively represent longitude and latitude, and it can be understood that the value "31" in the above formula is an offset to be added when longitude zone partition is calculated through east longitude in UTM projection; constants such as '36, 60' and the like are offset when the SRTM elevation data file is indexed through latitude and longitude; and 5, dividing the range for the longitude and latitude of the elevation file. It can be understood that the automation level of the park configuration process can be improved through the method, and in practical application, a user only needs to collect conventional data such as park range, pollution sources, monitoring site longitude and latitude and the like, and can automatically convert the data into modules such as calculation grid information and elevation information and the like through the accounting systemProfessional parameters of type.
The second step, for the simulation result deviation that factors such as effectual reduction garden monitoring facilities off-line lead to, the stability of reinforcing accounting system, consequently, still need carry out the unit conversion through garden information configuration module to the garden monitoring station data of acquireing to carry out validity verification to the garden monitoring station data of acquireing, its concrete processing mode is:
in the formula (I), the compound is shown in the specification,
and
respectively representing the monitoring concentration after unit conversion and before unit conversion, wherein MIN _ NSITE and MIN _ RSITE are the minimum effective monitoring site number and the minimum effective monitoring site proportion set for the configuration file, and Vflag represents whether the unit is effective or not;
is a conversion factor; it should be noted that, in the following description,
the calculation formula of (2) is as follows:
wherein P is a standard atmospheric pressure (Pa); r is reasonA desired gas constant; t is the ambient temperature (K); w is a group of m Relative molecular mass (g/mol);
the conversion ratio of the final mass unit and the gram after conversion;
unit proportion before conversion; e.g. from ppm to mg/m 3 When the temperature of the water is higher than the set temperature,
value of 10 3 ，
Value of 10 6 ；
Identifying the data validity of the ith monitoring station; and N is the total number of the monitored stations.
Thirdly, constructing a multi-mode atmospheric pollution diffusion model, wherein topographic and meteorological factors need to be considered when constructing the multi-mode atmospheric pollution diffusion model, the concentration distribution relation between the source intensity of the garden pollution source and the position of the garden monitoring station under the diffusion and atmospheric transmission actions needs to be obtained according to meteorological environment data and real-time data of the garden monitoring station in a simulation manner, for any monitoring station in the simulation range, the contribution degree of the pollution source intensity and the local concentration of the position of the monitoring station can be calculated according to the relationship between the pollution source intensity and the local concentration of the position of the monitoring station set by the multi-mode atmospheric pollution diffusion model, and then forming a contribution matrix by the contribution degree relation of a plurality of pollution sources and a plurality of monitoring stations, (namely obtaining the contribution matrix of the pollution source intensity and the monitoring stations in the current park) in order to obtain the influence of the pollution source emission on the monitoring stations, so that the specific steps of constructing the multi-mode atmospheric pollution diffusion model are as follows:
firstly (S1), in order to solve the problem that the superposition influence of a plurality of sources in a park on a monitoring station is converted into the solution of a multivariate linear equation set, under the premise of definitely optimizing a target and reducing the calculation difficulty, the linear superposition assumption that the data of a single monitoring station in the park is strong relative to the pollution sources in a plurality of parks needs to be customized, so that the calculation difficulty provided in the problem is reduced:
in the formula (I), the compound is shown in the specification,
as a result of monitoring site j, b is background concentration,
in order to be the emission intensity of the pollution source i,
the method is expressed as the contribution of a pollution source i to a monitoring station j, and can be understood that the calling of an atmospheric pollution diffusion model with large calculation amount is reduced through the algorithm design of linear superposition assumption and a contribution matrix, the optimization process is simplified, and the timeliness of the algorithm is improved;
secondly (S11), based on the linear superposition assumption proposed in the step S1, establishing a M x N x K three-dimensional contribution matrix of the data of the N monitoring stations of the park and the M pollution sources, wherein K is expressed as the number of the multi-mode atmospheric pollution diffusion models participating in calculation, so as to reduce the repeated calling of the multi-mode atmospheric pollution diffusion models in the optimization process.
It should be noted that, in an embodiment of the present invention, the present invention preferably takes a Sutton model as an example, that is,
the calculation is as follows:
wherein u is the average wind speed; z and h are the height of the receptor site and the height of the source respectively;
and
the diffusion coefficients are related to meteorological conditions, and the relative sizes of the light and shade change display contribution degrees of the M, N, K three-dimensional contribution matrixes are different, so that the repeated calling of the multi-mode atmospheric pollution diffusion model in the optimization process is reduced, and the algorithm efficiency is improved;
in order to solve the existing problems that the multi-mode atmospheric pollution diffusion model is called repeatedly in the optimization process and the algorithm efficiency is improved, besides the Sutton model, the method can also adopt AERMOD, CALPUFF, Gaussian models and the like, and the reason that the optimized Sutton model is directly calculated through a group of data to obtain the contribution of the pollution source i to the monitoring station j is understood to improve the algorithm efficiency.
Fourthly, constructing a source parameter inversion algorithm, and carrying out optimization estimation on the pollution source intensity according to the acquired data of the park monitoring station and the simulation result output by the multi-mode atmospheric pollution diffusion model so as to evaluate the effectiveness of the constructed multi-mode atmospheric pollution diffusion model,
the method for carrying out optimization estimation on the source intensity of the pollution source based on the source parameter inversion algorithm comprises the following steps:
s2, estimating the environmental background concentration of the park according to the initial guess of the pollution source intensity, the distribution situation of the pollution source intensity and the monitoring stations and the combination of the current wind field data, wherein the specific estimation steps are as follows:
firstly, starting from a monitoring site j, a vector pointing to any pollution source i
Unit vector with real-time wind direction
Performing dot product calculation to obtain vector
In that
Projection in a direction;
second, when for any contamination source i, if,
if the monitoring station is a windward station, the monitoring station is considered as an upwind station;
finally, when a plurality of upwind sites appear, averaging corresponding monitoring data to be used as the environmental background concentration of the park, and it can be understood that by the method, an accounting system can rapidly identify the upwind sites according to the real-time main wind direction through mathematical operation so as to dynamically acquire the background concentration;
s21, calculating a pseudo-inverse of a contribution matrix (M x N x K three-dimensional contribution matrix), after obtaining an initial guess of the source intensity of the pollution source, performing distribution optimization on the source intensity of the pollution source to obtain a reliable initial guess, and then performing optimization on the source intensity distribution of the pollution source, wherein the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is a row vector of 1 x M,
is a row vector of 1 x N,
is that
Pseudo-inverse of the matrix;
the initial guess of the pollution source strength can be understood, and by the method, more reliable initial guess can be obtained, the time required by optimization is reduced, and the situation that the optimization is caused by trapping in a local optimal solution and the optimization effect is poor is reduced;
s22, calculating the average relative deviation between the pollution source intensity and the measured value, taking the average relative deviation as the accuracy measurement of the multi-mode atmospheric pollution diffusion model calculation, dynamically generating the weight of the integration of the multi-mode atmospheric pollution diffusion model, avoiding the subjective interference caused by artificially setting the integration weight of the multi-mode atmospheric pollution diffusion model, and fully utilizing the characteristics of different models in different conditions, wherein the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is the average relative deviation of the multi-mode atmospheric pollution diffusion model K,
the set model is a configurable constant, and can be understood that the advantages of multiple models are integrated through fusion of the multiple models and dynamic distribution of model weights, so that the accuracy, universality and stability of the set model are improved;
fifthly, calculating the garden pollution emission total amount based on the pollution source intensity obtained through optimized estimation, wherein the influence of the pollution source outside the garden needs to be eliminated by marking the pollution source attribute when the garden pollution emission total amount is estimated due to the problems of irregular garden boundary, influence of emission sources outside the garden, such as enterprises outside the garden, living quarters and the like,
s3, the specific estimation step is:
in the formula, Q represents the total amount of the emission of the pollution source outside the garden;
based on the above technical concept, it should be noted that, when the total amount of the pollution emissions of the park is calculated, the validity of the accounting result of the multi-mode atmospheric pollution diffusion model needs to be verified by setting the maximum total amount of the emissions Q _ MAX and the maximum relative deviation RME _ MAX through the configuration file:
it can be understood that, because the pollution source is usually associated with an enterprise parcel, the relative magnitude estimation of the emission intensity in different areas of the park can be indirectly obtained, and the evaluation of the model on the interpretability of the monitoring data can be obtained through the consistency analysis of the simulation result and the monitoring result, that is, when the model fitting result is more consistent with the park monitoring data, the simulation kernel of the model is proved to be effective, and whether the simulation of the model is effective is proved, it can be understood that the model simulation result is not consistent with the monitoring data but is effective if the model simulation result is consistent with the monitoring data, and meanwhile, the model result is not completely consistent with the actual monitoring result.
In the actual calculation, the result of the total discharge amount of the park is the product of the intensity of the pollution source and the discharge time, and the elimination of the interference of the park external source needs to be considered, so that the influence of factors such as irregular park range, the pollution source outside the park boundary and the like on the accounting is reduced by setting the virtual source which does not participate in the statistics.
As an embodiment of the present invention, it should be noted that, in the above technical solution provided by the present invention, except that the park information configuration module in the first step requires a user to perform a survey and check according to the actual situation of the park to obtain accurate information to generate the configuration file, the steps processed by other modules are automatically performed by the accounting system through a timing task, so as to implement hourly accounting and statistics of the actual emission total amount of the atmospheric pollutants in the park.
As an embodiment of the present invention, as shown in fig. 2, a schematic diagram of basic information configuration of an industrial park is shown, in which a dot represents a position of a pollution source of the park, an inverted triangle represents a position of a monitoring site, a white dotted frame represents a boundary of a park simulation grid, and the park grid should cover all the pollution sources and the monitoring sites; at the same time, the user can select the required time,
as shown in fig. 3, a contour diagram (left) of the result of the campus atmospheric pollution diffusion simulation is shown, and it is understood that the horizontal and vertical coordinates in fig. 3 are X and Y coordinates of the UTM coordinate system; a pollution source contributes matrix (right) to the pollution concentration of the monitored site, and the contribution degree is displayed in the form of thermodynamic diagram;
as shown in fig. 4, a plot of the emissions of pollution from a campus is plotted (up), wherein the ordinate is the distance of the pollution source from the geographic center of the campus in the east-west and north-south directions, and the size of the plot indicates the relative magnitude of the emission intensity, and the four largest emissions sources at the current time are plotted; the plot of figure 4 (bottom) is a plot of simulated versus measured concentration of atmospheric pollutants for a campus.
While there have been shown and described the fundamental principles and essential features of the invention and advantages thereof, it will be apparent to those skilled in the art that the invention is not limited to the details of the foregoing exemplary embodiments, but is capable of other specific forms without departing from the spirit or essential characteristics thereof; the present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein, and any reference signs in the claims are not intended to be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.
Claims (7)
1. The online accounting method for the actual emission of the atmospheric pollutants based on multi-mode fusion is characterized by comprising the following steps of:
rolling simulation is carried out to obtain meteorological environment data in a garden area; constructing a simulation grid comprising park information, wherein the park information comprises park monitoring station basic information, park pollution source basic information and park elevation information;
constructing a multi-mode atmospheric pollution diffusion model, simulating and acquiring the concentration distribution relation between the source intensity of the pollution source of the park and the position of a park monitoring station under the diffusion and atmospheric transmission effects according to meteorological environment data and real-time data of the park monitoring station, and acquiring a contribution matrix of the source intensity of the pollution source in the current park to the monitoring station so as to obtain the influence of the emission of the pollution source on the monitoring station;
constructing a source parameter inversion algorithm, and carrying out optimization estimation on the source intensity of the pollution source according to the acquired data of the park monitoring station and the simulation result output by the multi-mode atmospheric pollution diffusion model so as to carry out effectiveness evaluation on the constructed multi-mode atmospheric pollution diffusion model; the method for carrying out optimization estimation on the source intensity of the pollution source based on the source parameter inversion algorithm comprises the following steps:
s2, estimating the environmental background concentration of the park according to the pollution source intensity and the distribution condition of the monitoring stations and by combining the current wind field data, wherein the specific estimation steps are as follows:
vector pointing to any pollution source i from monitoring site j
Unit vector with real-time wind direction
Performing dot product calculation to obtain vector
In that
Projection in a direction;
when for any contamination source i, if
If the monitoring station is a windward station, the monitoring station is considered as an upwind station;
when a plurality of upwind sites appear, averaging corresponding monitoring data to be used as the environmental background concentration of the park;
s21, calculating the pseudo inverse of the contribution matrix, obtaining an initial guess of the source intensity of the pollution source, performing distribution optimization on the source intensity of the pollution source to obtain a reliable initial guess, and then performing optimization on the source intensity distribution of the pollution source, wherein the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is a row vector of 1 x M,
is a row vector of 1 x N,
is that
Pseudo-inverses of the matrices;
is an initial guess of the source strength of the pollution source;
s22, calculating the average relative deviation between the pollution source intensity and the measured value, taking the average relative deviation as the accuracy measurement of the multi-mode atmospheric pollution diffusion model calculation, dynamically generating the weight of the integration of the multi-mode atmospheric pollution diffusion model to avoid the subjective interference caused by artificially setting the integration weight of the multi-mode atmospheric pollution diffusion model, and the specific calculation mode is as follows:
in the formula (I), the compound is shown in the specification,
is the average relative deviation of the multi-mode atmospheric pollution diffusion model K,
is a configurable constant;
and calculating the total pollutant emission amount of the park based on the pollution source intensity obtained through optimization estimation.
2. The multi-mode fused online accounting method for actual emission of atmospheric pollutants as claimed in claim 1, wherein the meteorological environment data comprises directly computable meteorological field data and indirectly computable meteorological field data; the simulation grids comprise a meteorological simulation grid and a campus simulation grid,
wherein, before extracting the data of the park meteorological environment, the position index of the park position in the meteorological simulation grid needs to be calculated firstly:
in the formula, i is expressed as a grid index corresponding to a pollution source in the meteorological simulation grid; j represents a grid index corresponding to a monitoring station in the meteorological simulation grid; x and y are respectively expressed as longitude and latitude information or UTM coordinates of the park;
and
respectively representing longitude and latitude or UTM coordinates of a starting point of the meteorological simulation grid; d is the grid resolution;
indicating a rounding down.
3. The online accounting method for the actual emission of the atmospheric pollutants fused in multiple modes according to claim 1, characterized in that before simulating the concentration distribution relationship between the garden pollution source intensity and the position of the garden monitoring station, the obtained data of the garden monitoring station needs to be converted in units and validated, and the specific processing mode is as follows:
in the formula (I), the compound is shown in the specification,
and
respectively representing the monitoring concentration after unit conversion and before unit conversion, wherein MIN _ NSITE and MIN _ RSITE are the minimum effective monitoring site number and the minimum effective monitoring site proportion set for the configuration file, and Vlag represents whether the unit is effective or not;
is a conversion factor;
the calculation formula of (2) is as follows:
wherein P is a standard atmospheric pressure (Pa); r is an ideal gas constant with the unit of J/mol.K; t is the ambient temperature (K); w m Relative molecular mass (g/mol);
the conversion ratio of the final mass unit and the gram after conversion;
unit proportion before conversion; e.g. from ppm to mg/m 3 When the temperature of the water is higher than the set temperature,
value of 10 3 ，
Value of 10 6 ；
Identifying the data validity of the ith monitoring station; and N is the total number of the monitored stations.
4. The multimode-fused online accounting method for the actual emission of the atmospheric pollutants according to claim 1 or 2, wherein the park simulation grid is constructed by the following specific method:
in the formula, nx is the number of grids of the park simulation grid in the x direction,
and
maximum and minimum UTM-X coordinates of a pollution source and a monitoring point of the park respectively; the offset is the expansion of the grid, and the value is equal to the grid resolution;
representing the UTM area in which the longitude and latitude information is located; lng denotes longitude.
5. The multi-mode fused online accounting method for the actual emission of the atmospheric pollutants as claimed in claim 1, wherein the specific way of constructing the multi-mode atmospheric pollution diffusion model is as follows:
s1, customizing a linear superposition hypothesis that the data of a single monitoring station of the garden are strong relative to pollution sources in a plurality of gardens to reduce the calculation difficulty:
in the formula (I), the compound is shown in the specification,
is the monitoring result of the monitoring station j, b is the background concentration,
in order to be the emission intensity of the pollution source i,
expressed as the contribution of pollution source i to monitored site j;
s11, establishing data of N monitoring stations of the park and an M x N x K three-dimensional contribution matrix of M pollution sources based on a linear superposition hypothesis, wherein K is expressed as the number of the multi-mode atmospheric pollution diffusion models participating in calculation so as to reduce the repeated calling of the multi-mode atmospheric pollution diffusion models in the optimization process:
wherein u is the average wind speed; z and h are the height of the receptor site and the height of the source respectively;
and
is in phase with meteorological conditionsThe relative magnitudes of the contributions are shown by the shading of the three-dimensional contribution matrix of M x N x K, the off-diffusion coefficient.
6. The multimode-fused online accounting method for the actual emission of the atmospheric pollutants according to claim 1, wherein the specific steps of calculating the total emission of the pollution in the park are as follows:
s3, marking the attribute of the pollution source to eliminate the influence of the pollution source outside the garden:
in the formula, Q represents the total amount of the pollutant source emission outside the garden.
7. The multi-mode fused online accounting method for the actual emission of the atmospheric pollutants according to claim 1 or 6, wherein when the total emission of the pollution of the park is calculated, the validity of the accounting result of the multi-mode atmospheric pollution diffusion model is further verified by setting a maximum total emission Q _ MAX and a maximum relative deviation RME _ MAX through a configuration file:
when the fitting result of the multi-mode atmospheric pollution diffusion model is consistent with the park monitoring data, the simulation accounting of the multi-mode atmospheric pollution diffusion model is proved to be effective.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210660280.3A CN114757807B (en) | 2022-06-13 | 2022-06-13 | Multi-mode fused online accounting method for actual emission of atmospheric pollutants |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210660280.3A CN114757807B (en) | 2022-06-13 | 2022-06-13 | Multi-mode fused online accounting method for actual emission of atmospheric pollutants |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114757807A CN114757807A (en) | 2022-07-15 |
| CN114757807B true CN114757807B (en) | 2022-09-13 |
Family
ID=82337231
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210660280.3A Active CN114757807B (en) | 2022-06-13 | 2022-06-13 | Multi-mode fused online accounting method for actual emission of atmospheric pollutants |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114757807B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115144548B (en) * | 2022-08-31 | 2022-11-18 | 天津市环鉴环境检测有限公司 | Harmful gas composition real-time monitoring system and monitoring method thereof |
| CN115600919B (en) * | 2022-09-19 | 2023-08-22 | 江苏蓝创智能科技股份有限公司 | Method for real-time unorganized emission location and total amount of campus emissions calculation |
| CN115545565A (en) * | 2022-11-24 | 2022-12-30 | 江苏省生态环境大数据有限公司 | Method and system for managing and controlling total amount of pollution discharged from park based on atmospheric environment quality |
| CN115792137B (en) * | 2023-01-17 | 2023-06-06 | 河北先河环保科技股份有限公司 | Atmospheric pollution tracing method and device and terminal |
| CN116402408B (en) * | 2023-06-06 | 2023-09-26 | 北京工业大学 | Site concentration difference-based local emission contribution estimation method |
| CN116562712B (en) * | 2023-06-29 | 2023-09-19 | 内江师范学院 | System and method for predicting air quality |
| CN116739224B (en) * | 2023-08-14 | 2023-11-03 | 中科三清科技有限公司 | Method, apparatus and storage medium for determining contribution concentration of heavy point pollution source |
| CN117454729B (en) * | 2023-11-13 | 2024-04-19 | 交通运输部规划研究院 | Atmospheric pollution contribution degree prediction method and device and readable storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111444633A (en) * | 2020-04-20 | 2020-07-24 | 中国科学院大气物理研究所 | Quantitative analysis method and system for atmospheric pollution process |
| CN114155129A (en) * | 2021-11-29 | 2022-03-08 | 中科三清科技有限公司 | Atmospheric environment traceability method and system based on industrial park |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107436343A (en) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | It is a kind of to simulate the method for calculating sensitizing range pollutant concentration |
| CN107871210B (en) * | 2017-11-03 | 2022-03-22 | 南开大学 | Atmospheric environment capacity accounting method |
| CN110472782A (en) * | 2019-08-01 | 2019-11-19 | 软通动力信息技术有限公司 | A kind of data determination method, device, equipment and storage medium |
| CN111537023B (en) * | 2020-05-13 | 2021-05-07 | 浙江大学 | Diffusion simulation and source tracing method for atmospheric pollutants in industrial park |
-
2022
- 2022-06-13 CN CN202210660280.3A patent/CN114757807B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111444633A (en) * | 2020-04-20 | 2020-07-24 | 中国科学院大气物理研究所 | Quantitative analysis method and system for atmospheric pollution process |
| CN114155129A (en) * | 2021-11-29 | 2022-03-08 | 中科三清科技有限公司 | Atmospheric environment traceability method and system based on industrial park |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114757807A (en) | 2022-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114757807B (en) | Multi-mode fused online accounting method for actual emission of atmospheric pollutants | |
| CN109213964A (en) | A kind of satellite AOD product bearing calibration for merging multi-source feature geographic factor | |
| CN111487216B (en) | Carbon dioxide flux inversion method and system | |
| CN114818391B (en) | Pollution source concentration analysis method based on multi-tracing model | |
| CN113436045A (en) | Atmospheric pollution multi-scale fusion traceability method and device and readable storage medium | |
| CN113624921B (en) | Multi-mode comprehensive pollution tracing method | |
| CN113297528A (en) | NO based on multi-source big data2High-resolution space-time distribution calculation method | |
| CN109726867B (en) | High-resolution electric power weather forecasting method based on multi-mode set | |
| CN111027175A (en) | Method for evaluating social and economic influences of flood based on coupling model integrated simulation | |
| CN110766191A (en) | Newly-added PM2.5 fixed monitoring station site selection method based on space-time kriging interpolation | |
| CN112016696B (en) | PM integrating satellite observation and ground observation 1 Concentration inversion method and system | |
| CN114254802B (en) | Prediction method for vegetation coverage space-time change under climate change drive | |
| CN108983324A (en) | A kind of temperature forecast method and system based on Kalman filtering | |
| CN115759488A (en) | Carbon emission monitoring and early warning analysis system and method based on edge calculation | |
| CN105303254A (en) | Method and device for prediction of radiation received by photovoltaic power station | |
| CN114880933A (en) | Atmospheric temperature and humidity profile inversion method and system for non-exploration-site foundation microwave radiometer based on reanalysis data | |
| CN114077798B (en) | Method for estimating grid concentration of atmospheric pollutants in small-scale area | |
| CN113095586A (en) | Short-term multi-meteorological-element forecasting method based on deep neural network | |
| CN104239662A (en) | Method for quantitatively evaluating developing potential of regional wind energy resources | |
| CN114462316A (en) | Air quality optimization simulation method for homeland space planning | |
| Liu et al. | A Bagging-GBDT ensemble learning model for city air pollutant concentration prediction | |
| CN115544706A (en) | Wavelet and XGboost model integrated atmospheric fine particle concentration estimation method | |
| Chen et al. | Variation of gross primary productivity dominated by leaf area index in significantly greening area | |
| CN114878748B (en) | CO (carbon monoxide) 2 Emission monitoring method and emission monitoring system | |
| CN115511192A (en) | Rainfall forecasting method and system based on lightning data assimilation |
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 |