CN111220549A - Method for measuring and calculating pollutant emission surface concentration of area to be measured - Google Patents
Method for measuring and calculating pollutant emission surface concentration of area to be measured Download PDFInfo
- Publication number
- CN111220549A CN111220549A CN201811327592.2A CN201811327592A CN111220549A CN 111220549 A CN111220549 A CN 111220549A CN 201811327592 A CN201811327592 A CN 201811327592A CN 111220549 A CN111220549 A CN 111220549A
- Authority
- CN
- China
- Prior art keywords
- discharge surface
- concentration
- pollutant
- sampling
- area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/55—Specular reflectivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The invention provides a method for measuring and calculating the concentration of a pollutant emission surface of a region to be measured, which comprises the following steps: acquiring a discharge surface of the pollutant discharge surface optimized spatial concentration distribution model, which is established based on a pollutant spatial concentration distribution model constructed in advance and actual distribution of pollutant concentration in the pollutant discharge surface; and the emission surface of the optimized space concentration distribution model is subjected to integral calculation by applying the optimized space concentration distribution model, and the pollutant emission surface concentration of the area to be measured is measured and calculated. The invention provides a method for constructing a pollutant spatial concentration distribution model according to the actual distribution characteristics of the pollutant concentration in the pollutant discharge surface of the area to be measured, the constructed model can more comprehensively, truly and accurately reflect the gas concentration distribution characteristics of the pollutant discharge surface of the area to be measured, and the concentration of the pollutant discharge surface of the area to be measured is calculated by using the model, so that the accuracy is higher and the precision is higher.
Description
Technical Field
The invention relates to the field of environmental monitoring and protection, in particular to a method for measuring and calculating the concentration of a pollutant emission surface of a region to be measured.
Background
With the increasingly strict national environmental protection policy, the requirements of environmental protection departments on enterprises involved in pollutant production are also increasingly strict, and 20 national municipalities and municipalities have issued pollutant discharge charging policies currently in China in order to further promote the enterprises to add pollutant treatment facilities and reduce pollutant discharge amount. Under the background, countries and enterprises also put higher requirements on the accuracy of the pollutant discharge amount measuring method.
At present, some pollutant emission amount measuring and calculating methods are designed for many units and enterprises, but when the pollutant emission amount of an area to be measured is measured and calculated, the concentration characteristic of an emission source at a monitoring position is mostly described by adopting the concept of average concentration, the concentration distribution characteristic of an emission surface at the monitoring position and an acquisition method of the surface concentration of the emission surface are not given, the calculating method is simple and rough, and the accuracy of the pollutant emission amount calculation is influenced.
Disclosure of Invention
The invention provides a method for measuring and calculating the concentration of a pollutant emission surface of a region to be measured, which aims to solve the technical problems and at least partially solve the technical problems.
According to an exemplary embodiment of the present disclosure, there is provided a method for measuring and calculating a concentration of a pollutant discharge area of an area to be measured, the method including:
acquiring a discharge surface of the pollutant discharge surface optimized spatial concentration distribution model, which is established based on a pollutant spatial concentration distribution model constructed in advance and actual distribution of pollutant concentration in the pollutant discharge surface; and the emission surface of the optimized space concentration distribution model is subjected to integral calculation by applying the optimized space concentration distribution model, and the pollutant emission surface concentration of the area to be measured is measured and calculated.
Optionally, the optimized spatial concentration distribution model is constructed by the following steps:
and actually measuring and acquiring the actually measured line integral concentration of each sampling light path in the plurality of sampling light paths in the pollutant discharge surface by using optical remote measuring equipment.
And obtaining the theoretical line integral concentration of each of the plurality of sampling light paths on the pollutant discharge surface of the area to be detected by applying a pre-constructed pollutant spatial concentration distribution model according to the path of each of the plurality of sampling light paths.
And according to the actual measurement line integral concentration of each sampling light path obtained by actual measurement, the theoretical line integral concentration of each sampling light path obtained by applying a pre-constructed pollutant space concentration distribution model according to the path of each sampling light path in the plurality of sampling light paths, and by taking the minimum sum of the square differences of the actual measurement line integral concentration and the theoretical line integral concentration as a principle, obtaining the pollutant discharge surface optimized space concentration distribution model by optimizing the pre-constructed pollutant space concentration distribution model.
Optionally, the theoretical line integral concentration of each of the multiple sampling optical paths on the pollutant discharge surface of the region to be measured is obtained by the following steps:
obtaining a path of the sampling light path in the pollutant discharge surface; and
and respectively calculating the theoretical line integral concentration of the sampling light path in the sampling light paths on the pollutant discharge surface of the area to be detected according to the acquired path of the sampling light path and a pre-constructed pollutant space concentration distribution model.
Optionally, the theoretical line integral concentration e of the ith sampling optical path in the multiple sampling optical paths on the pollutant discharge surface of the region to be detectediIs ei＝∫G(P1,P2,P3…)dli，
Wherein, G (P)1,P2,P3…) is the pre-constructed pollutant space concentration distribution model;
(P1,P2,P3…) are the model parameters to be optimized;
lifor the path of the ith sampling optical path, i belongs to [1, N ∈]And N is the number of sampling optical paths.
Optionally, the step of obtaining the optimized spatial concentration distribution model by the optimization processing is as follows:
step 1, model parameters (P) to be optimized1,P2,P3…) in the value range, randomly setting M groups of model parameters (P)1j,P2j,P3j…) (j ═ 1 … M), and determined so that
Minimum model parameter (P)1a,P2a,P3a…)；
Step 2, model parameters (P) to be optimized1,P2,P3…) is in the range of values, in order to (P)1a,P2a,P3a…) is an initial value to
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle1J,P2J,P3J…); and
multiple execution of the two steps prevents the optimal parameter (P) from being obtained1J,P2J,P3J…) falls into a locally optimal solution,
wherein the content of the first and second substances,
cithe actual measurement of the ith sampling optical path in the multiple sampling optical paths obtained by the optical remote-control equipmentLine integral concentration; and
∫G(P1j,P2j,P3j…)dlibased on the jth group model parameter (P) for the ith sampling optical path in the multiple sampling optical paths on the pollutant discharge surface of the region to be detected1j,P2j,P3j…) is calculated.
Optionally, each of the plurality of sampling optical paths on the pollutant discharge surface of the region to be measured is set to be distributed as follows:
and a first sampling light path in the plurality of sampling light paths on the pollutant discharge surface is distributed at the position where the maximum integral concentration of the pollutant actual measurement line in the discharge surface is obtained by using the optical remote-control equipment.
The number of the sampling optical paths on two sides of the first sampling optical path is in direct proportion to the numerical relationship of the area of the two parts of the discharge surface divided by the first sampling optical path.
The sampling optical paths in two portions into which the discharge surface is divided by the first sampling optical path are uniformly distributed in the portions.
Optionally, when the maximum position of the integrated concentration of the actual pollutant measurement line in the pollutant discharge surface is difficult to confirm, defining a discharge surface diagonal line as a first sampling optical path; or all sampling light paths are uniformly arranged in the pollutant discharge surface directly according to the principle of area sharing.
Optionally, the monitoring points of the optical telemetry device are distributed on the boundary of the exhaust surface.
Optionally, when the area of the pollutant discharge surface is larger or the number of sampling optical paths on the pollutant discharge surface is smaller, the monitoring point of the optical telemetry equipment is arranged at the intersection point of the boundary of the pollutant discharge surface.
In another aspect, the present invention provides a machine-readable storage medium having instructions stored thereon, where the instructions are used for causing a machine to execute any one of the methods for measuring and calculating the pollutant emission surface concentration of an area to be measured.
Through the technical scheme, the method for constructing the spatial concentration distribution model of the pollutant discharge surface according to the actual distribution characteristics of the pollutant concentration in the pollutant discharge surface of the area to be measured is provided, the constructed model can reflect the gas concentration distribution characteristics of the pollutant discharge surface of the area to be measured more comprehensively, truly and accurately, and the concentration accuracy and the precision of the pollutant discharge surface of the area to be measured are higher by utilizing the model to measure and calculate. Additional features and advantages of embodiments of the invention will be set forth in the detailed description which follows.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the embodiments of the invention without limiting the embodiments of the invention.
FIG. 1 is a flow chart illustrating a method for measuring a concentration of a pollutant discharge area in a test area according to an exemplary embodiment;
FIG. 2 is a flow chart illustrating a process for constructing an optimized spatial concentration profile of a pollutant discharge surface in accordance with an exemplary embodiment;
FIG. 3 is a flowchart illustrating a method for calculating a theoretical line integral concentration of an ith sampling optical path of a plurality of sampling optical paths on a pollutant discharge surface of a region to be measured according to an exemplary embodiment; FIG. 4 is a schematic diagram of the distribution of 4 sampling optical paths in a discharge plane according to an exemplary embodiment; and
fig. 5 is a schematic diagram showing the distribution of 8 sampling optical paths in a discharge plane according to an exemplary embodiment.
Description of the reference numerals
OPQR pollutant discharge surface O monitoring point
Sampling optical path with maximum actual measurement line integral concentration in OA pollutant discharge surface
Other multiple sampling optical paths in OB, OC, OD, OE, OF, OG and OH pollutant discharge surfaces
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating embodiments of the invention, are given by way of illustration and explanation only, not limitation.
FIG. 1 is a flow chart illustrating a method for measuring a concentration of a pollutant discharge area in a test area according to an exemplary embodiment; as shown in fig. 1, the method includes the following steps.
In step S101, a discharge surface of the pollutant discharge surface optimized spatial concentration distribution model is obtained, the optimized spatial concentration distribution model being established based on a pollutant spatial concentration distribution model constructed in advance and an actual distribution of pollutant concentration within the pollutant discharge surface, the construction of the pollutant discharge surface optimized spatial concentration distribution model being described in detail below with reference to fig. 2.
In step S102, the emission surface of the optimized spatial concentration distribution model is subjected to integral calculation by using the optimized spatial concentration distribution model, and the concentration of the pollutant emission surface of the region to be measured is measured.
Optimized spatial concentration distribution model construction as will be described in detail below in conjunction with fig. 2, in an exemplary embodiment of the present invention, the pollutant discharge surface spatial concentration distribution model may be G (P)1,P2,P3…), wherein (P)1,P2,P3…) are the model parameters to be optimized.
In the exemplary embodiment of the present invention, the actual distribution characteristics of the pollutant concentration in the pollutant discharge surface are described by the optimized spatial concentration distribution model of the pollutant in the discharge surface, and in this embodiment, the optimized spatial concentration distribution model of the pollutant in the discharge surface is mainly related to the positions of the plurality of sampling light paths, so the distribution rationality of the plurality of sampling light paths determines the accuracy of the construction of the optimized spatial concentration distribution model of the pollutant discharge surface.
In an exemplary embodiment of the present invention, each of the plurality of sampling optical paths on the pollutant discharge surface of the region to be measured is configured to be distributed as follows:
a first sampling light path in the plurality of sampling light paths on the pollutant discharge surface is distributed at the position where the maximum integral concentration of the pollutant actual measurement line in the discharge surface is obtained by using optical remote-control equipment; the number of the sampling optical paths on two sides of the first sampling optical path is in direct proportion to the numerical relationship of the area of the two parts of the discharge surface divided by the first sampling optical path; the sampling optical paths in two portions into which the discharge surface is divided by the first sampling optical path are uniformly distributed in the portions.
For example, fig. 4 is a schematic diagram of the distribution of 4 sampling optical paths in the discharge plane according to an exemplary embodiment.
As shown in fig. 4, the number N of sampling optical paths is set to 4, the polygon OPQR is a discharge surface, and O is a monitoring point of the optical telemetry device. Firstly, the optical remote measuring device is used for continuously changing the position of the sampling light path, universal sampling is carried out on the discharge surface OPQR, and the position of the sampling light path with the maximum measured line integral concentration in the discharge surface OPQR is determined through data comparison, for example, OA in figure 4 is the sampling light path with the maximum measured line integral concentration. The sampling optical path OA divides the discharge surface OPQR into two parts a1 and a2, where a1 is a polygon OAQR and has an area SA1A2 is a polygon OPA with an area SA2. The number of the A1 and A2 in-plane sampling optical paths N1 and N2 is distributed according to the ratio relation of the A1 and A2 areas. Suppose SA1/SA2The first sampling optical path OA is removed, the remaining sampling optical paths are 3, exactly 2 sampling optical paths are arranged in a1, and 1 sampling optical path, namely N, is arranged in a21/N2The OA of the boundary between a1 and a2 is another sampling optical path, which is 4 in total. If SA1/SA2When the ratio is not 2 but 2.3 or 2.4, it is still possible to set a1 to include 2 sampling paths and a2 to include 1 sampling path according to the rough ratio. Aiming at the design problem of 2 sampling optical paths in the A1 area, the patent provides the principle that the sampling optical paths in two parts formed by dividing the discharge surface by the first sampling optical path are uniformly distributed in the part, namely the principle of equal area. Designing sampling optical paths OB and OC to divide A1 into 3 parts uniformly, namely polygons OAQB, OBC and OCR have equal areas; similarly, in the region a2, the optical path OD is designed according to the principle of area sharing, so that the areas of the polygons OPD and ODA are equal.
Fig. 5 is a schematic diagram showing the distribution of 8 sampling optical paths in a discharge plane according to an exemplary embodiment. As shown in FIG. 5, for example, 8 sampling optical paths are taken, and the polygon SPQR is the discharge surface and OA is the discharge surfaceAnd the position of the sampling light path with the maximum integral concentration of the internal measurement line is O as a monitoring point of the optical telemetering equipment. Suppose the area ratio S of A1 and A2A1/SA2Setting 5 sampling light paths in A1, setting 2 sampling light paths in A2, and using OA as another sampling light path by the boundary between A1 and A2, wherein the total number of sampling light paths is 8; in the area A1, designing sampling optical paths OB, OC, OD, OE and OF according to the principle OF area sharing, namely, the areas OF polygons OAQB, OBC, OCD, ODE, OERF and OFS are equal; in the area A2, sampling optical paths OG and OH are designed according to the principle of area sharing, namely the areas of polygons OPH, OHG and OGA are equal.
Further, when the position with the maximum actually measured line integral concentration in the pollutant discharge surface is difficult to confirm, a discharge surface diagonal line is defined as a first sampling optical path; or all sampling light paths are uniformly arranged in the pollutant discharge surface directly according to the principle of area sharing.
Furthermore, monitoring points of the optical remote measuring equipment are distributed on the boundary of the discharge surface; when the area of the pollutant discharge surface is larger or the number of sampling light paths on the pollutant discharge surface is smaller, the sampling light paths are arranged at the boundary intersection point of the discharge surface.
The design method of the multiple sampling light paths can monitor the discharge characteristics of pollutants on the discharge surface more reasonably and comprehensively, and further establish a more accurate optimized space concentration distribution model.
Through the technical scheme, the invention provides the method for constructing the optimal spatial concentration distribution model of the pollutant discharge surface according to the actual distribution characteristics of the pollutant concentration in the pollutant discharge surface of the area to be measured, namely according to the reasonably distributed positions of the plurality of sampling light paths and the actually measured line integral concentration.
Fig. 2 is a flow chart illustrating a process for constructing an optimized spatial concentration distribution model of a pollutant discharge surface, as shown in fig. 2, according to an exemplary embodiment, the optimized spatial concentration distribution model is constructed by the following steps:
in step S201, actually measuring and acquiring the actually measured line integral concentration of each of the plurality of sampling optical paths by using an optical telemetry device.
In an exemplary embodiment of the invention, the optical telemetry equipment used is composed of an equipment body and a reflector, a sampling light path is arranged between the equipment body and the reflector, and the optical equipment can directly measure the actual line integral concentration of pollutants on the sampling light path.
In step S202, a theoretical line integral concentration of each of the plurality of sampling optical paths on the pollutant discharge surface of the region to be measured is calculated.
For example, fig. 3 is a flowchart illustrating a method for calculating a theoretical line integral concentration of an ith sampling optical path of a plurality of sampling optical paths on a pollutant discharge surface of a region to be measured according to an exemplary embodiment, where as shown in fig. 3, the theoretical line integral concentration of the ith sampling optical path of the plurality of sampling optical paths on the pollutant discharge surface of the region to be measured includes the following steps:
in a step S301 of the process,
and acquiring the path of the ith sampling optical path in the plurality of sampling optical paths on the pollutant discharge surface.
In step S302, a theoretical line integral concentration of an ith sampling optical path of the multiple sampling optical paths on the pollutant discharge surface of the region to be measured is calculated.
For example, the theoretical line integral concentration of the ith sampling optical path in the multiple sampling optical paths on the pollutant discharge surface of the region to be measured may be ═ G (P)1,P2,P3…)dli；
Wherein, G (P)1,P2,P3…) is the pre-constructed pollutant space concentration distribution model;
(P1,P2,P3…) are the model parameters to be optimized; liFor the path of the ith sampling optical path, i belongs to [1, N ∈]And N is the number of sampling optical paths.
In S203, the optimal spatial concentration distribution model of the pollutant discharge surface is constructed.
In an exemplary embodiment of the present invention, the optimal spatial concentration distribution model of the pollutant discharge surface is constructed through optimization processing according to an actually measured line integral concentration of each sampling light path obtained through the actual measurement, a theoretical line integral concentration of each sampling light path obtained through the optimal spatial concentration distribution model of the pollutant discharge surface, and a principle that a sum of squares of the actually measured line integral concentration and the theoretical line integral concentration is minimum.
For example, the optimization process obtains the optimized spatial concentration distribution model by the following steps:
step 1, model parameters (P) to be optimized1,P2,P3…) in the value range, randomly setting M groups of model parameters (P)1j,P2j,P3j…) (j ═ 1 … M), and determined so that
Minimum model parameter (P)1a,P2a,P3a…)；
Step 2, model parameters (P) to be optimized1,P2,P3…) is in the range of values, in order to (P)1a,P2a,P3a…) as an initial value to
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle1J,P2J,P3J…); and
multiple execution of the two steps prevents the optimal parameter (P) from being obtained1J,P2J,P3J…) falls into a locally optimal solution,
wherein the content of the first and second substances,
cithe measured line integral concentration of the ith sampling optical path in the multiple sampling optical paths obtained by the optical remote-control equipment is obtained; and
∫G(P1j,P2j,P3j…)dlifor the ith of multiple sampling light paths on the pollutant discharge surface of the region to be measuredThe sampling optical path is based on the jth group model parameter (P)1j,P2j,P3j…) is calculated.
Further, in an exemplary embodiment of the present invention, the method for calculating the concentration of the pollutant discharge area of the area to be measured further includes: optimizing the parameter (P)1J,P2J,P3J…) into the pollutant space concentration distribution model to obtain an optimized space concentration distribution model G (P)1J,P2J,P3J…), and for G (P)1J,P2J,P3J…) is calculated, and the pollutant discharge surface concentration is obtained.
In exemplary embodiments, the contaminants referred to herein may be VOCsSIn practical application, the method for measuring and calculating the concentration of the pollutant discharge surface of the area to be measured can be suitable for measuring and calculating the concentration of the pollutant discharge surface of various pollutants.
Through the technical scheme, the method for designing the sampling light path according to the equipartition area principle after the sampling light path with the maximum actually measured line integral concentration in the discharge surface is obtained is provided, the concentration distribution characteristics of pollutants in the discharge surface can be more accurately obtained, and the optimal spatial concentration distribution model of the pollutant discharge surface is more reasonable based on the concentration distribution characteristics of the pollutants in the discharge surface. The model is applied to the emission surface of the optimized space concentration distribution model to measure and calculate the concentration of the pollutant emission surface of the area to be measured, the measuring and calculating method is more reasonable, and the measuring and calculating reliability, precision and accuracy can be effectively improved.
Although the embodiments of the present invention have been described in detail with reference to the accompanying drawings, the embodiments of the present invention are not limited to the details of the above embodiments, and various simple modifications can be made to the technical solutions of the embodiments of the present invention within the technical idea of the embodiments of the present invention, and the simple modifications all belong to the protection scope of the embodiments of the present invention.
It should be noted that the various features described in the above embodiments may be combined in any suitable manner without departing from the scope of the invention. In order to avoid unnecessary repetition, the embodiments of the present invention do not describe every possible combination.
Those skilled in the art will understand that all or part of the steps in the method according to the above embodiments may be implemented by a program, which is stored in a storage medium and includes several instructions to enable a single chip, a chip, or a processor (processor) to execute all or part of the steps in the method according to the embodiments of the present application. 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 addition, any combination of various different implementation manners of the embodiments of the present invention is also possible, and the embodiments of the present invention should be considered as disclosed in the embodiments of the present invention as long as the combination does not depart from the spirit of the embodiments of the present invention.
Claims (9)
1. A method for measuring and calculating the concentration of pollutant emission surface of an area to be measured is characterized by comprising the following steps:
acquiring a discharge surface of the pollutant discharge surface optimized spatial concentration distribution model, which is established based on a pollutant spatial concentration distribution model constructed in advance and actual distribution of pollutant concentration in the pollutant discharge surface; and
and (4) applying the optimized spatial concentration distribution model to the discharge surface of the optimized spatial concentration distribution model to perform integral calculation, and measuring and calculating the pollutant discharge surface concentration of the region to be measured.
2. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 1,
the optimized spatial concentration distribution model is constructed by the following steps:
actually measuring and acquiring the actually measured line integral concentration of each sampling light path in a plurality of sampling light paths arranged in the pollutant discharge surface by using optical remote measuring equipment;
according to the path of each of the plurality of sampling optical paths, obtaining the theoretical line integral concentration of each of the plurality of sampling optical paths on the pollutant discharge surface of the area to be tested by applying a pre-constructed pollutant spatial concentration distribution model; and
and according to the actual measurement line integral concentration of each sampling light path obtained by actual measurement, the theoretical line integral concentration of each sampling light path obtained by applying a pre-constructed pollutant space concentration distribution model according to the path of each sampling light path in the plurality of sampling light paths, and by taking the minimum sum of the square differences of the actual measurement line integral concentration and the theoretical line integral concentration as a principle, obtaining the pollutant discharge surface optimized space concentration distribution model by optimizing the pre-constructed pollutant space concentration distribution model.
3. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 1,
the theoretical line integral concentration of each of a plurality of sampling light paths on the pollutant discharge surface of the area to be measured is obtained by the following steps:
obtaining a path of the sampling light path in the pollutant discharge surface; and
and respectively calculating the theoretical line integral concentration of the sampling light path in the sampling light paths on the pollutant discharge surface of the area to be detected according to the acquired path of the sampling light path and a pre-constructed pollutant space concentration distribution model.
4. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 3,
and the theoretical line integral concentration e of the ith sampling light path in the multiple sampling light paths on the pollutant discharge surface of the area to be detectediComprises the following steps:
ei＝∫G(P1,P2,P3…)dli，
wherein the content of the first and second substances,
G(P1,P2,P3…) is the pre-constructed pollutant space concentration distribution model;
(P1,P2,P3…) are the model parameters to be optimized;
lifor the path of the ith sampling optical path, i belongs to [1, N ∈]And N is the number of sampling optical paths.
5. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 2,
the optimization process obtains an optimized spatial concentration distribution model by the steps of:
step 1, model parameters (P) to be optimized1,P2,P3…) in the value range, randomly setting M groups of model parameters (P)1j,P2j,P3j…) (j ═ 1 … M), and determined so that
Minimum model parameter (P)1a,P2a,P3a…)；
Step 2, model parameters (P) to be optimized1,P2,P3…) is in the range of values, in order to (P)1a,P2a,P3a…) is an initial value to
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle1J,P2J,P3J…); and
multiple execution of the two steps prevents the optimal parameter (P) from being obtained1J,P2J,P3J…) falls into a locally optimal solution,
wherein the content of the first and second substances,
cithe measured line integral concentration of the ith sampling optical path in the multiple sampling optical paths obtained by the optical remote-control equipment is obtained; and
∫G(P1j,P2j,P3j…)dlibased on the jth group model parameter (P) for the ith sampling optical path in the multiple sampling optical paths on the pollutant discharge surface of the region to be detected1j,P2j,P3j…) is calculated.
6. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 1,
each sampling light path in a plurality of sampling light paths on the pollutant discharge surface of the area to be measured is set to be distributed as follows:
a first sampling optical path in the plurality of sampling optical paths on the pollutant discharge surface is distributed at the maximum position of the actually measured line integral concentration of the pollutant in the discharge surface obtained by the optical remote-control equipment;
the number of the sampling optical paths on two sides of the first sampling optical path is in direct proportion to the numerical relationship of the area of the two parts of the discharge surface divided by the first sampling optical path;
the sampling optical paths in two portions into which the discharge surface is divided by the first sampling optical path are uniformly distributed in the portions.
7. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 6,
when the maximum position of the integral concentration of the actual pollutant measuring line in the pollutant discharge surface is difficult to confirm, defining a discharge surface diagonal line as a first sampling optical path;
or all sampling light paths are uniformly arranged in the pollutant discharge surface directly according to the principle of area sharing.
8. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 7,
monitoring points of the optical telemetering equipment are distributed on the boundary of the discharge surface;
when the area of the pollutant discharge surface is larger or the number of sampling optical paths on the pollutant discharge surface is smaller, the monitoring point of the optical telemetering equipment is arranged at the intersection point of the boundary of the pollutant discharge surface.
9. A machine-readable storage medium having instructions stored thereon for causing a machine to perform any one of the methods for measuring and calculating the pollutant emission surface concentration of an area under test described above.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811327592.2A CN111220549B (en) | 2018-11-08 | 2018-11-08 | Method for measuring and calculating pollutant emission surface concentration of area to be measured |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811327592.2A CN111220549B (en) | 2018-11-08 | 2018-11-08 | Method for measuring and calculating pollutant emission surface concentration of area to be measured |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111220549A true CN111220549A (en) | 2020-06-02 |
| CN111220549B CN111220549B (en) | 2023-01-31 |
Family
ID=70828566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811327592.2A Active CN111220549B (en) | 2018-11-08 | 2018-11-08 | Method for measuring and calculating pollutant emission surface concentration of area to be measured |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111220549B (en) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5748318A (en) * | 1996-01-23 | 1998-05-05 | Brown University Research Foundation | Optical stress generator and detector |
| US5813767A (en) * | 1995-09-29 | 1998-09-29 | Finmeccanica S.P.A. Azienda Ansaldo | System and a method for monitoring combustion and pollutants in a thermal plant by means of laser diodes |
| US6251679B1 (en) * | 1998-10-26 | 2001-06-26 | Aerodyne Research, Inc. | Method for measuring pollutant formation |
| EP1292821A1 (en) * | 2000-06-19 | 2003-03-19 | Sira Limited | Apparatus and method for analysing the constituents of a gas or a gas like cloud or plume |
| US6542242B1 (en) * | 1999-05-10 | 2003-04-01 | University Of Washington | Mapping air contaminants using path-integrated optical remote sensing with a non-overlapping variable path length beam geometry |
| US20070131864A1 (en) * | 2005-12-08 | 2007-06-14 | Ellis Mark T | System and method for detecting and remediating contamination |
| CN104881546A (en) * | 2015-06-01 | 2015-09-02 | 中国科学院上海高等研究院 | Method for improving prediction efficiency of atmospheric pollution model |
| CN106383205A (en) * | 2016-10-13 | 2017-02-08 | 北京伟瑞迪科技有限公司 | VOCs (Volatile Organic Compounds) region on-line monitoring and early warning system |
| CN106841087A (en) * | 2017-04-17 | 2017-06-13 | 中国石油化工股份有限公司 | Petrochemical enterprise uncontrollable discharge source VOCS discharges method for measuring and calculating |
| CN106909801A (en) * | 2017-04-17 | 2017-06-30 | 中国神华能源股份有限公司 | Calculate the method and system of the effective small hourly value of gaseous contaminant |
| CN107314983A (en) * | 2016-06-12 | 2017-11-03 | 中国石油化工股份有限公司 | Method based on the infrared measuring and calculating refinery device VOCS discharge capacitys of open type |
| CN107436343A (en) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | It is a kind of to simulate the method for calculating sensitizing range pollutant concentration |
| CN107703516A (en) * | 2017-10-13 | 2018-02-16 | 天津大学 | A kind of long-range measure Pollution From Ships tapping equipment and its application process |
| CN207503015U (en) * | 2017-09-21 | 2018-06-15 | 成都帝图信息技术有限公司 | Decentralized type sewage treating stations intelligent remote supervisory terminal |
| US20180210469A1 (en) * | 2017-01-26 | 2018-07-26 | International Business Machines Corporation | Dynamic emission discharge reduction |
| CN108398529A (en) * | 2018-03-06 | 2018-08-14 | 深圳市卡普瑞环境科技有限公司 | A kind of information acquisition method and information acquisition device |
-
2018
- 2018-11-08 CN CN201811327592.2A patent/CN111220549B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5813767A (en) * | 1995-09-29 | 1998-09-29 | Finmeccanica S.P.A. Azienda Ansaldo | System and a method for monitoring combustion and pollutants in a thermal plant by means of laser diodes |
| US5748318A (en) * | 1996-01-23 | 1998-05-05 | Brown University Research Foundation | Optical stress generator and detector |
| US6251679B1 (en) * | 1998-10-26 | 2001-06-26 | Aerodyne Research, Inc. | Method for measuring pollutant formation |
| US6542242B1 (en) * | 1999-05-10 | 2003-04-01 | University Of Washington | Mapping air contaminants using path-integrated optical remote sensing with a non-overlapping variable path length beam geometry |
| EP1292821A1 (en) * | 2000-06-19 | 2003-03-19 | Sira Limited | Apparatus and method for analysing the constituents of a gas or a gas like cloud or plume |
| US20070131864A1 (en) * | 2005-12-08 | 2007-06-14 | Ellis Mark T | System and method for detecting and remediating contamination |
| CN104881546A (en) * | 2015-06-01 | 2015-09-02 | 中国科学院上海高等研究院 | Method for improving prediction efficiency of atmospheric pollution model |
| CN107314983A (en) * | 2016-06-12 | 2017-11-03 | 中国石油化工股份有限公司 | Method based on the infrared measuring and calculating refinery device VOCS discharge capacitys of open type |
| CN106383205A (en) * | 2016-10-13 | 2017-02-08 | 北京伟瑞迪科技有限公司 | VOCs (Volatile Organic Compounds) region on-line monitoring and early warning system |
| US20180210469A1 (en) * | 2017-01-26 | 2018-07-26 | International Business Machines Corporation | Dynamic emission discharge reduction |
| CN106841087A (en) * | 2017-04-17 | 2017-06-13 | 中国石油化工股份有限公司 | Petrochemical enterprise uncontrollable discharge source VOCS discharges method for measuring and calculating |
| CN106909801A (en) * | 2017-04-17 | 2017-06-30 | 中国神华能源股份有限公司 | Calculate the method and system of the effective small hourly value of gaseous contaminant |
| CN107436343A (en) * | 2017-07-31 | 2017-12-05 | 南京南瑞集团公司 | It is a kind of to simulate the method for calculating sensitizing range pollutant concentration |
| CN207503015U (en) * | 2017-09-21 | 2018-06-15 | 成都帝图信息技术有限公司 | Decentralized type sewage treating stations intelligent remote supervisory terminal |
| CN107703516A (en) * | 2017-10-13 | 2018-02-16 | 天津大学 | A kind of long-range measure Pollution From Ships tapping equipment and its application process |
| CN108398529A (en) * | 2018-03-06 | 2018-08-14 | 深圳市卡普瑞环境科技有限公司 | A kind of information acquisition method and information acquisition device |
Non-Patent Citations (6)
| Title |
|---|
| QU HONGQUAN: "A Dynamic Method to Estimate Source Emission Rate and Predict Contaminant Concentrations", 《2009 2ND INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
| 专祥涛: "《最优化方法基础》", 30 April 2018, 武汉大学出版社 * |
| 刘志明等: "基于红外掩日通量法(SOF)污染源排放气体浓度时空分布反演算法研究", 《物理学报》 * |
| 朱亮: "光学气体成像技术在泄漏检测与维修中的应用研究", 《泄露检测与维修专栏》 * |
| 贾润中: "OP－FTIＲ 技术在评估炼化装置VOCs排放量中的应用研究", 《现代化工》 * |
| 金岭等: "利用SOF-FTIR技术监测化工厂区VOCs排放", 《大气与环境光学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111220549B (en) | 2023-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109446696B (en) | CMAQ model-based rapid atmospheric environment capacity measuring and calculating method, storage medium and terminal | |
| Walker et al. | An evaluation of the predictive accuracy of wake effects models for offshore wind farms | |
| Lei et al. | Roundness error evaluation algorithm based on polar coordinate transform | |
| CN103809290B (en) | A kind of mutual compensation optimizing method of optical system face shape error | |
| CN105813020A (en) | RSSI corrected wireless sensor network positioning algorithm of self-adaptive environment | |
| Moroni et al. | Optimal inspection strategy planning for geometric tolerance verification | |
| Barini et al. | Uncertainty analysis of point-by-point sampling complex surfaces using touch probe CMMs: DOE for complex surfaces verification with CMM | |
| CN103472672A (en) | Correction method of optical proximity correction model | |
| CN116933099B (en) | Automobile part aging test method, system and medium | |
| CN112197749B (en) | Cross calibration method and device for effective wave height of wave buoy | |
| CN109960241B (en) | Automatic test method and system for automobile ECU internal calibration meter | |
| CN109931903A (en) | A kind of cylindricity assessment method based on improvement whale optimization algorithm | |
| CN114241266B (en) | Optical cable intelligent early warning processing method and system | |
| CN109900789B (en) | Imaging method for internal defects of longitudinal section of tree | |
| CN111220549B (en) | Method for measuring and calculating pollutant emission surface concentration of area to be measured | |
| CN111767212B (en) | Software code quality evaluation method, device, equipment and storage medium | |
| CN103439869B (en) | The method of measurement pattern density | |
| CN110146182B (en) | Method, device and system for measuring zero-order vortex light beam topological load | |
| CN106682252B (en) | Coaxiality rapid evaluation method based on minimum entity state | |
| CN109212751B (en) | Method for analyzing tolerance of free-form surface | |
| CN112784392A (en) | Method and system for simulating dynamic distribution of pollutants in area to be tested | |
| CN101013022B (en) | Method for random error elimination in optical element interference sampling data | |
| CN111157038B (en) | Measuring and calculating method for measuring pollutant discharge amount in to-be-measured area | |
| CN103793786B (en) | One remanufactures bent axle circular runout processing quality control method | |
| KR101444794B1 (en) | Method for Constructing the Pseudo Hot Pin Power Distribution Using Detector-Signal-Based Planar Radial Peaking Factor in COLSS |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230105 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Applicant after: CHINA PETROLEUM & CHEMICAL Corp. Applicant after: Sinopec Safety Engineering Research Institute Co.,Ltd. Address before: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Applicant before: CHINA PETROLEUM & CHEMICAL Corp. Applicant before: SINOPEC Research Institute OF SAFETY ENGINEERING |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |