CN111220549B - 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
- CN111220549B CN111220549B CN201811327592.2A CN201811327592A CN111220549B CN 111220549 B CN111220549 B CN 111220549B CN 201811327592 A CN201811327592 A CN 201811327592A CN 111220549 B CN111220549 B CN 111220549B
- Authority
- CN
- China
- Prior art keywords
- discharge surface
- pollutant
- concentration
- area
- sampling
- 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
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 emission 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 applying the optimized spatial concentration distribution model to perform integral calculation on the discharge surface of the optimized spatial concentration distribution model, and measuring and calculating the concentration of the pollutant discharge surface of the area to be measured.
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 actually measured 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 on the basis of the minimum sum of the squared differences of the actually measured line integral concentration and the theoretical line integral concentration, optimizing the pre-constructed pollutant space concentration distribution model to obtain the pollutant discharge surface optimized 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 detected i Is e i ＝∫G(P 1 ,P 2 ,P 3 …)dl i ，
Wherein, G (P) 1 ,P 2 ,P 3 8230) is the pre-constructed pollutant spatial concentration distribution model;
(P 1 ,P 2 ,P 3 8230) are the model parameters to be optimized;
l i for the path of the ith sampling optical path, i ∈ [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 optimized 1 ,P 2 ,P 3 8230in the value range, randomly setting M groups of model parameters (P) 1j ,P 2j ,P 3j 8230j =1 8230M) and were determined to be such that
Minimum model parameter (P) 1a ,P 2a ,P 3a …)；
Step 2, model parameters (P) to be optimized 1 ,P 2 ,P 3 8230in the value range of (P) 1a ,P 2a ,P 3a 8230in the formula) is an initial value
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle 1J ,P 2J ,P 3J 8230The); and
multiple execution of the two steps prevents the optimal parameter (P) from being obtained 1J ,P 2J ,P 3J 8230), trapping in a locally optimal solution,
wherein the content of the first and second substances,
c i the 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(P 1j ,P 2j ,P 3j …)dl i based 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 detected 1j ,P 2j ,P 3j 8230); is determined by the theoretical line integral concentration.
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 real-time measuring line in the discharge surface is obtained by using the optical remote-measuring 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 according to the principle of area sharing directly.
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 modeling an optimized spatial concentration profile of a pollutant discharge surface, according to 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
O monitoring point of OPQR pollutant discharge surface
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 ,P 2 ,P 3 8230wherein (P) 1 ,P 2 ,P 3 8230) 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 of sampling paths is setQuantity N =4, polygon OPQR is the discharge face and O is the monitoring point of the optical telemetry device. Firstly, the position of the sampling optical path is continuously changed by using optical telemetry equipment, the emission surface OPQR is generally sampled, and the position of the sampling optical path with the maximum measured line integral concentration in the emission surface OPQR is determined through data comparison, for example, OA in fig. 4 is the sampling optical 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 S A1 A2 is a polygonal OPA and has an area S A2 . And distributing the number of the A1 and A2 in-plane sampling optical paths N1 and N2 according to the ratio relation of the A1 and A2 areas. Suppose S A1 /S A2 =2, remove the first OA of sampling optical path, the remaining sampling optical path is 3, set up 2 sampling optical paths in exactly A1, set up 1 sampling optical path in A2, N 1 /N 2 And =2, the OA of the boundary between A1 and A2 is another sampling optical path, and the total number is 4. If S A1 /S A2 When the ratio is not equal to 2, but is 2.3 or 2.4, it is still possible to set A1 to include 2 sampling optical paths and A2 to include 1 sampling optical path according to the rough proportional relationship. Aiming at the design problem of 2 sampling optical paths in an 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 equally dividing the 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 area 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, taking 8 sampling optical paths as an example, the polygon SPQR is the discharge surface, OA is the sampling optical path position where the measured line integral concentration is maximum in the discharge surface, and O is the monitoring point of the optical telemetry device. Suppose the area ratio S of A1 and A2 A1 /S A2 Setting up 5 sampling light paths in A1, setting up 2 sampling light paths in A2, and setting up another sampling light path by OA between A1 and A2, wherein the total number is 8; in the area A1, according to the principle OF sharing the area, the sampling optical paths OB, OC, OD, OE and OF are designed, namely the area phases OF polygons OAQB, OBC, OCD, ODE, OERF and OFSEtc.; in the area A2, according to the principle of equally dividing the area, sampling optical paths OG and OH are designed, namely polygons OPH, OHG and OGA are equal in area.
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 according to the principle of area sharing directly.
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 emission characteristics of pollutants on the emission surface more reasonably and comprehensively, and further establish a more accurate optimal spatial 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 method, the step S,
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 ,P 2 ,P 3 …)dl i ；
Wherein, G (P) 1 ,P 2 ,P 3 8230) is the pre-constructed pollutant spatial concentration distribution model;
(P 1 ,P 2 ,P 3 8230) are the model parameters to be optimized; l i For the path of the ith sampling optical path, i ∈ [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 pollutant discharge area optimized spatial concentration distribution model is constructed through optimization processing according to an actually measured line integral concentration of each sampling optical path obtained through the actual measurement, a theoretical line integral concentration of each sampling optical path obtained through the pollutant discharge area optimized spatial concentration distribution model, and a principle that a sum of squared differences 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 optimized 1 ,P 2 ,P 3 8230in the value range, randomly setting M groups of model parameters (P) 1j ,P 2j ,P 3j 8230j) (j =1 8230M) and is determined so that
Minimum model parameter (P) 1a ,P 2a ,P 3a …)；
Step 2, model parameters (P) to be optimized 1 ,P 2 ,P 3 \8230;) value range with (P) 1a ,P 2a ,P 3a 8230treating the crude oil as an initial value
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle 1J ,P 2J ,P 3J 8230; and
multiple execution of the two steps prevents the optimal parameter (P) from being obtained 1J ,P 2J ,P 3J 8230quickly) trapping in the local optimal solution,
wherein the content of the first and second substances,
c i the 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(P 1j ,P 2j ,P 3j …)dl i based 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 detected 1j ,P 2j ,P 3j 8230); etc.).
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 ,P 2J ,P 3J 8230, (G) obtaining an optimized spatial concentration distribution model by introducing into the spatial concentration distribution model of the contaminantsP 1J ,P 2J ,P 3J 8230and G (P) 1J ,P 2J ,P 3J 8230and) performing integral calculation to obtain the concentration of the pollutant discharge surface.
In exemplary embodiments, the contaminants referred to herein may be VOCs S In 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 various media capable of storing program codes.
In addition, any combination of various different implementation manners of the embodiments of the present invention can be made, and the embodiments of the present invention should also be regarded as the disclosure of 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 (8)
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
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 concentration of the pollutant discharge surface of the area to be measured;
the optimal spatial concentration distribution model of the pollutant discharge surface is related to the positions of a plurality of sampling light paths, and each sampling light path in the plurality of sampling light paths on the pollutant discharge surface of the area to be detected is set 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 maximum position of the pollutant actual measurement line integral concentration 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 the two parts into which the discharge surface is divided by the first sampling optical path are uniformly distributed in the part.
2. The method for measuring and calculating the pollutant discharge surface concentration of area to be measured according to claim 1, wherein,
the optimized spatial concentration distribution model is constructed by the following steps:
actually measuring and obtaining 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 actually measured 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 on the basis of the minimum sum of the squared differences of the actually measured line integral concentration and the theoretical line integral concentration, optimizing the pre-constructed pollutant space concentration distribution model to obtain the pollutant discharge surface optimized 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 detected i Comprises the following steps:
e i ＝∫G(P 1 ,P 2 ,P 3 …)dl i ，
wherein the content of the first and second substances,
G(P 1 ,P 2 ,P 3 8230) is the pre-constructed pollutant spatial concentration distribution model;
(P 1 ,P 2 ,P 3 8230) are the model parameters to be optimized;
l i for the path of the ith sampling optical path, i ∈ [1, N ]]And N is the number of sampling optical paths.
5. The method for measuring and calculating the pollutant discharge surface concentration of area to be measured according to claim 2, wherein,
the optimization process obtains an optimized spatial concentration distribution model by the steps of:
step 1, model parameters (P) to be optimized 1 ,P 2 ,P 3 8230in the value range, randomly setting M groups of model parameters (P) 1j ,P 2j ,P 3j 8230j =1 8230M) and were determined to be such that
Minimum model parameter (P) 1a ,P 2a ,P 3a …)；
Step 2, model parameters (P) to be optimized 1 ,P 2 ,P 3 8230in the value range of (P) 1a ,P 2a ,P 3a 8230in the formula) is an initial value
Using a simulated annealing method to obtain the optimal parameter (P) based on the minimum principle 1J ,P 2J ,P 3J 8230; and
repeatedly executing the above two steps to preventThe obtained optimal parameter (P) 1J ,P 2J ,P 3J 8230quickly) trapping in the local optimal solution,
wherein the content of the first and second substances,
c i the measured line integral concentration of the ith sampling light path in the multiple sampling light paths obtained by the optical remote detection equipment is obtained; and
∫G(P 1j ,P 2j ,P 3j …)dl i based on the jth group of model parameters (P) for the ith sampling optical path in the multiple sampling optical paths on the pollutant discharge surface of the region to be detected 1j ,P 2j ,P 3j 8230); is determined by the theoretical line integral concentration.
6. The method for measuring and calculating the pollutant discharge surface concentration of area to be measured according to claim 1, wherein,
when the maximum position of the integral concentration of the actual measurement line of the pollutant in the pollutant discharge surface is difficult to confirm, defining a discharge surface diagonal line as a first sampling light path;
or all sampling light paths are uniformly arranged in the pollutant discharge surface directly according to the principle of area sharing.
7. The method for measuring and calculating the concentration of pollutant discharge surface in area to be measured according to claim 6,
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 less, the monitoring point of the optical telemetering equipment is arranged at the intersection point of the boundary of the discharge surface.
8. 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 CN111220549A (en) | 2020-06-02 |
| CN111220549B true 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 (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
| CN207503015U (en) * | 2017-09-21 | 2018-06-15 | 成都帝图信息技术有限公司 | Decentralized type sewage treating stations intelligent remote supervisory terminal |
| CN108398529A (en) * | 2018-03-06 | 2018-08-14 | 深圳市卡普瑞环境科技有限公司 | A kind of information acquisition method and information acquisition device |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0766080A1 (en) * | 1995-09-29 | 1997-04-02 | FINMECCANICA S.p.A. AZIENDA ANSALDO | System and method for monitoring combustion and pollutants 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 |
| GB0015015D0 (en) * | 2000-06-19 | 2000-08-09 | Sira Electro Optics Ltd | Apparatus and method for analysing the constituents of a gas or a gas like cloud or plume |
| US8235627B2 (en) * | 2005-12-08 | 2012-08-07 | Ellis Mark T | System and method for detecting and remediating contamination |
| US10628757B2 (en) * | 2017-01-26 | 2020-04-21 | International Business Machines Corporation | Dynamic emission discharge reduction |
| CN107703516B (en) * | 2017-10-13 | 2020-12-11 | 天津大学 | Device for remotely measuring ship pollution emission and application method thereof |
-
2018
- 2018-11-08 CN CN201811327592.2A patent/CN111220549B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
| 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 |
| CN108398529A (en) * | 2018-03-06 | 2018-08-14 | 深圳市卡普瑞环境科技有限公司 | A kind of information acquisition method and information acquisition device |
Non-Patent Citations (2)
| Title |
|---|
| OP－FTIＲ 技术在评估炼化装置VOCs排放量中的应用研究;贾润中;《现代化工》;20171031;第37卷(第10期);第177-180页 * |
| 光学气体成像技术在泄漏检测与维修中的应用研究;朱亮;《泄露检测与维修专栏》;20140430;第14卷(第4期);第14-16页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111220549A (en) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jiang et al. | Measuring urban sprawl in Beijing with geo-spatial indices | |
| Kocabas et al. | Assessing cellular automata model behaviour using a sensitivity analysis approach | |
| CN109446696B (en) | CMAQ model-based rapid atmospheric environment capacity measuring and calculating method, storage medium and terminal | |
| CN104732307B (en) | Project work amount acquisition methods and system | |
| Walker et al. | An evaluation of the predictive accuracy of wake effects models for offshore wind farms | |
| CN112197749B (en) | Cross calibration method and device for effective wave height of wave buoy | |
| CN103902451A (en) | Intelligent electricity meter software quality evaluation method | |
| US20180137218A1 (en) | Systems and methods for similarity-based information augmentation | |
| CN103472672A (en) | Correction method of optical proximity correction model | |
| CN111220549B (en) | Method for measuring and calculating pollutant emission surface concentration of area to be measured | |
| CN111222678A (en) | Road surface technical condition prediction method | |
| Hall et al. | The application of isogeometric analysis to the neutron diffusion equation for a pincell problem with an analytic benchmark | |
| CN108182503B (en) | Quality evaluation method and equipment of industrial centrifugal pump | |
| CN111767212B (en) | Software code quality evaluation method, device, equipment and storage medium | |
| CN113486295A (en) | Fourier series-based total ozone change prediction method | |
| Mirats-Tur et al. | Leak detection and localization using models: field results | |
| Wang et al. | The influence and adjust method of hyperparameters’ prior distributions in Bayesian calibration for building stock energy prediction | |
| CN108062435B (en) | Fatigue life calibration method based on nominal stress method | |
| La Sorte et al. | Range maps and species richness patterns: errors of commission and estimates of uncertainty | |
| JP5433559B2 (en) | Source estimation apparatus and source estimation method for diffusion materials | |
| KR101444794B1 (en) | Method for Constructing the Pseudo Hot Pin Power Distribution Using Detector-Signal-Based Planar Radial Peaking Factor in COLSS | |
| Probst | New techniques in noise prediction | |
| Campanelli et al. | Variance gradients and uncertainty budgets for nonlinear measurement functions with independent inputs | |
| US10365230B1 (en) | Scatterometry overlay based on reflection peak locations | |
| Kopmann et al. | Using reliability analysis in morphodynamic simulation with TELEMAC-2D/SISYPHE |
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 |
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 |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |