CN102967689B - A kind of Pollution Source Monitoring points distributing method - Google Patents

A kind of Pollution Source Monitoring points distributing method Download PDF

Info

Publication number
CN102967689B
CN102967689B CN201210479898.6A CN201210479898A CN102967689B CN 102967689 B CN102967689 B CN 102967689B CN 201210479898 A CN201210479898 A CN 201210479898A CN 102967689 B CN102967689 B CN 102967689B
Authority
CN
China
Prior art keywords
pollution source
region
pollution
source
unknown
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.)
Expired - Fee Related
Application number
CN201210479898.6A
Other languages
Chinese (zh)
Other versions
CN102967689A (en
Inventor
尤学一
王淋淋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201210479898.6A priority Critical patent/CN102967689B/en
Publication of CN102967689A publication Critical patent/CN102967689A/en
Application granted granted Critical
Publication of CN102967689B publication Critical patent/CN102967689B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of pollution source discrimination method based on related coefficient, comprise and build storehouse and identify two processes, build storehouse and comprise Sensor monitoring point is set; Area of space segmentation and the supposition of unknown source position; Simulation unknown source dispose procedure; Obtain the concentration information of all unknown sources release; Set up unknown source position and concentration information database; Identifying adopts concentration information degree of correlation coupling and Nearest Neighbor Method, and when true pollution occurs, based on pollutant levels information, utilize correlation analysis method to find out the maximum location point of the degree of correlation in database, this location point is pollution source point of release; Pollution Source Monitoring points distributing method then adopts different principles of layouting for strong, weak convective region characteristic area, compared with prior art, pollution source discrimination method based on related coefficient improves recognition efficiency and accuracy, quantity of compared with prior art layouting Pollution Source Monitoring points distributing method reduces, and is efficient points distributing method; And lay a good foundation for pollution source discrimination method of the present invention.

Description

A kind of Pollution Source Monitoring points distributing method
Technical field
The present invention relates to indirect problem and know source domain, particularly relate to the identification field of unknown pollution source off-position in the interior space.
Background technology
To the identification of unknown source, more at area researches such as heat transfer, groundwater contamination transmission, atmospheric pollution transmission.Along with the leakage of the epiphytotics outbursts such as SARS and harmful gas, air pollution problem inside rooms is come into one's own day by day, to the identification of indoor pollutant release information, become the emphasis of research, determine pollution source release information fast and accurately, to taking counter-measure effectively, alleviate personnel injury and property loss significant.
More research is had to the identification of the unknown pollution source in indoor, Sohn et al. (list of references [1]) and Sreedharan et al. (list of references [2]) discharge pollutant in the space be made up of multiple room of three layers, adopt the position of Bayesian probability statistical method identification pollution source, intensity and duration.The method can be used for identifying the release information of pollution source, but need adopt and organize sensor more, and will be built with the sampling database of thousands of groups of sampled points, cost is high in a particular application, workload is large simultaneously.
In addition, Liu et al. (list of references [3]) based on the multizone modeling method of followed probability, the pollutant position discharged in the interior space that fast track is formed at multiple compartment.Arvelo et al. (list of references [4]) etc. utilizes genetic algorithm to come the position of the pollution source that identification discharges in the space be made up of nine compartments and porch.Said method is when identifying pollutant release information, and all adopt multizone mathematical model, computational accuracy is more coarse, only can obtain the macroscopic information of pollutant release.Therefore, when pollution source identification, need to build and a kind ofly can meet lower workload and less operating cost, the method polluting accuracy of identification can be improved again.
Fluid Mechanics Computation (CFD) simulation is a kind of simple, quick, economic method, by the simulation of CFD, can obtain more accurate, detailed information of flow and pollute transmission diffuse information.Therefore, it is more and more for pollution sources information identification, such as: Liu et al. (list of references [5]) etc., by CFD simulation means, oppositely solve the off-position of transient state point source by employing united state probabilistic method.Zhang et al. is according to against solving method and the reversible solving method of puppet, being carried out the pollution source (list of references [6,7]) of identification single-point, instantaneous relase by CFD inverse.Cai et al. (list of references [8]) supposes that pollution source are known, sets up pollution source position and concentration data storehouse, introduce the method for similar parameter, identify unknown source of release position and intensity by CFD simulation.The method applied in the present invention is also based on CFD simulated implementation.
In pollution sources information identification, pollution source off-position is vital, also be difficult to identify (list of references [9]), instantaneous or the sustained release problem for single point source, once pollution source position is determined, just can based on CFD method simulating pollution thing dispose procedure and CONCENTRATION DISTRIBUTION situation, by the linear relationship between pollution source intensity and pollutant levels, utilize linear scale relation can determine actual pollutant release strength (list of references [7,8]).Based on above-mentioned analysis, emphasis of the present invention inquires into the identification of pollution source off-position.
List of references:
[1]Sohn M,Reynolds P,Singh N,et al.Rapidly locating and characterizing pollutant releases in buildings[J].Journal of Air and Waste Management Association,2002,52:1422-1432
[2]Sreedharan P,Sohn M,Gadgil A,et al.Systems approach to evaluating sensor characteristics forreal-time monitoring of high-risk indoor contaminant releases[J].Atmospheric Environment,2006,40(19):3490-3502
[3]Liu X,Zhai Z.Prompt tracking of indoor airborne contaminant source location with probablity-basedinverse multi-zone modeling[J].Building and Environment,2009,44(6):1135-1143.
[4]Arvelo J,Brandt A,Roger R P,et al.2002.An enhancement multizone model and its application tooptimum placement of CBW sensors[J].ASHRAE Trans.108(2):818-825
[5]Liu X,Zhai Z.Location identification for indoor instantaneous point contaminant source byprobability-based inverse computational fluid dynamics modeling[J].Indoor Air,2008,18(1):2-11
[6]Zhang T,Chen Q.Identification of contaminant sources in enclosed environments by inverse CFDmodeling[J].Indoor Air,2007,17(3):167-177
[7]Zhang T,Chen Q.Identification of contaminant sources in enclosed spaces by a single sensor[J].Indoor Air,2007,17(6):439-449.
[8]Cai H.,Li X,Long W.,et al.Real-time determination of indoor contaminant source location andstrength,Part I:with one sensor[J].Proceedings:Building Simulation,2007,1272-1278
[9] Zhang Tengfei. the indirect problem modeling [J] in identification room air pollution source. building heat-energy air conditioner, 2008,27 (6): 18-23
Summary of the invention
Based on above-mentioned technical matters, the present invention proposes a kind of pollution source discrimination method based on related coefficient, utilize the simulation of unknown pollution source to carry out the degree of correlation with the actual monitoring result occurred when polluting and mate the related coefficient obtained, realize pollution source accurate recognition, and, carry out according to strong and weak convective region the method that pollution monitoring layouts.
The present invention proposes a kind of pollution source discrimination method based on related coefficient, comprise and build storehouse and identify two processes, it is characterized in that, the method comprises the following steps:
Step one, arranges Sensor monitoring point in pollution detection environment, and simulation unknown source dispose procedure, obtains pollutant levels message sample data; Carry out area of space segmentation and the supposition of unknown pollution source position of pollution detection environment, by unknown pollution source position and pollutant levels information building database;
Step 2, adopts concentration information degree of correlation coupling and Nearest Neighbor Method, when true contamination accident occurs, based on the pollutant levels information of monitoring, utilize correlation analysis method, find out the location point that the degree of correlation in database is maximum, this location point is the pollution source point of release of wanted identification.
Described Sensor monitoring point adopts single-sensor, and its position is set in the center in interior space given row air port; When there being multiple discharge air port, should select the center of air quantity maximum exhaust outlet that sensor is set; When multiple exhaust outlet air quantity is consistent, chosen distance location of air supply inlet exhaust outlet center farthest arranges sensor.
The concrete process of the segmentation of described area of space and the supposition of unknown pollution source position comprises: interior space region is evenly divided into physical dimension, little space cell P that shape is consistent; If interior space out-of-shape, with the volume shared by each space cell for partitioning standards, and ensure that shape is each other more consistent; After marking off junior unit, at the center of each junior unit, unknown pollution source are set, pollution source are designated as successively
S={S 1,S 2,...,S N}。
Described unknown source release and identification process specifically comprise following process:
Adopt the release of numerical simulation software simulating pollution source, the concentration information of each unknown source obtained at monitoring point place, is designated as M={ α successively 1, α 2..., α n, and the concentration information monitored in true dispose procedure is designated as α 0, the concrete monitor value of each concentration information is designated as
α i={c i1,c i2,...,c im},0≤i≤N,1≤j≤m,m=τ/Δt,m>2
τ is the time period of monitoring, and Δ t is the time interval of monitoring;
By reality, the concentration information polluted in the concentration information of monitoring when occurring and unknown pollution source database carries out Controlling UEP, adopts following relatedness computation formula:
R ( α i , α 0 ) = Σ j = 1 m ( c ij - c ij ‾ ) ( c 0 j - c 0 j ‾ ) Σ j = 1 m ( c ij - c ij ‾ ) 2 Σ j = 1 m ( c 0 j - c 0 j ‾ ) 2 , 1 ≤ i ≤ N - - - ( 1 )
with represent the mean value of concentration sequence respectively;
Obtain the maximum pollution source location point of relevance degree by formula (1), this point is most probable point of release.
The invention allows for a kind of Pollution Source Monitoring points distributing method, the method comprises the following steps:
For interior space region, the mesh Reynolds number in flow field is obtained by the simulation of CFD method, according to the height of mesh Reynolds number, flow field, space is divided into strong convection district and two, weak convective region characteristic area, different principles of layouting is adopted for each characteristic area, namely adopt encryption to layout principle to strong convection district, sparse principle of layouting is adopted to weak convective region.
The mesh Reynolds number maximal value in flow field, described space is designated as Re max, then get
R 0 = 1 100 Re max - - - ( 2 )
To Re < R 0region, be defined as weak convective region, to Re > R 0region, be defined as strong convection district.
To strong convection district: according to delimited region, become by this Region dividing physical dimension, shape close to consistent little space cell Pa, wherein Pa < P; To weak convective region: become by this Region dividing physical dimension, shape close to consistent little space cell Pb, Pb > P, wherein interior space region is evenly divided into physical dimension, little space cell that shape is consistent by P.
Compared with prior art, pollution source discrimination method based on related coefficient of the present invention is known in source domain in indirect problem, improve recognition efficiency, increase the accuracy of identification, on the other hand, Pollution Source Monitoring points distributing method of the present invention is compared with uniform stationing technical method of the prior art, and quantity of layouting reduces; Be not only efficient points distributing method.And lay a good foundation for pollution source discrimination method of the present invention, two methods achieve jointly finds pollution source point of release fast and accurately, and then takes measure timely and effectively to carry out harm reduction, reduce the loss.
Accompanying drawing explanation
Fig. 1 meeting room structural representation: 1-1,3-D view; 1-2,2 d plane picture
Strong convection district under Fig. 2 two dimensional model schematic diagram: the 2-1 that layouts represents sample point and point of release schematic diagram; 2-2 represents the legend mark of all point of releases;
Region-wide uniform stationing schematic diagram under Fig. 3 two dimensional model;
The weak convective region of Fig. 4 adopts the schematic diagram after sparse layouting and legend mark;
Fig. 5 is the pollution source discrimination method schematic flow sheet based on related coefficient of the present invention.
Reference numeral: 1-lamp; 2-air outlet; 3-air inlet; 4-cupboard; 5-people; 6-desk; 7-cupboard; S1-S8-pollution source point of release; R-Sensor monitoring point;
Embodiment
Illustrate below in conjunction with accompanying drawing and concrete scene of implementing, the technical scheme based on the pollution source identification of related coefficient and the recognition methods of efficiently layouting of improvement of the present invention be described in detail:
Such as, being identified object is a three-dimensional meeting room.Set up meeting room three-dimensional model as shown in Figure 1, its long (X) × wide (Y) × height (Z) is 9m × 4m × 3.2m, and have two air inlets 3 (0.4m × 0.4m) at the top in room, air output is 0.128m 3/ s; Be air outlet 2 (0.8m × 0.4m) in the centre position at top, room, interior space wall be set to the adiabatic wall of non-slippage; In the interior space, suppose there is eight members participated in a conference, be all sitting in fixing position, be highly 1m during sitting, all having corresponding desk in face of everyone, is highly 0.6m, simulate the oral area that potential viral communication source is positioned at people, it is designated as S1-S8 successively.Assuming that only have a pollution source releasing virus at every turn, release strength is 50units/s, only has a sensor R in the interior space, and its coordinate is (X, Y, Z)=(4.5,2,3.5) m, is positioned at the center position of air outlet 2, and detecting of sensor is limited to 1ppm.The position of potential pollution source is as shown in table 1.
Table 1 potential pollution source off-position
Embodiment is as follows:
Utilize CFD method, by unsteady state circumstance, pulsatile flow field is obtained to Indoor Flow Field;
The Numerical-Mode polluted before occurring plans to build the storehouse stage: eight in this example potential pollution source location points, with the source strength of 100units/s release pollutant, are collected the concentration information of pollutant, set up pollution source position and concentration information Sample Storehouse by sensor.
Simulate the dispose procedure of point to be identified: be divided into and launch simulation 8 times, with the recognition effect of verification method at diverse location place, the source strength of each release is 50units/s, at the sensing station place of simulation, setting detects and is limited to 1ppm, obtain the concentration information of pollutant, set up to be identified some Sample Storehouse.
Calculate the relevance degree R of concentration information in the concentration data and Sample Storehouse of monitoring, select three time periods such as 180s, 240s, 300s to analyze recognition result.
Can find out within selected tri-time periods of 180s, 240s, 300s by the recognition result of table 2, eight diverse location release pollution source all can accurately identify.
The recognition result of unknown pollution source in table 2 different time sections
In order to illustrate that the efficient fabric point methods of improvement is to the beneficial effect identified, show implementation process intuitively simultaneously, be described for two dimensional model below, in Fig. 2,2-1 is a two dimensional model, region in dotted line frame is strong convection district, this part is encrypted and layouts, corresponding recognition result is analyzed by not encrypting to layout and increase in dotted line to layout in contrast dotted line frame, the position of three points that the degree of correlation obtained as can be seen from table 3 is maximum, encryption is adopted to layout (1 ' to strong convection district, 2 ', , , 9 ') after, the relevance degree identified comparatively unencryption increases before layouting, the position identified is also more accurate.
Recognition result after table 3 strong convection district adopts encryption to layout
To two dimensional model as shown in Figure 3, strong convection district is roughly on the upside of straight line L, downside is in weak convective region, and Fig. 3 adopts overall uniform stationing to know method for distinguishing, and the sample point of laying is more, Fig. 4 adopts sparse points distributing method to weak convective region, lay sample point less, table 4 give weak convective region adopt sparse layout after recognition result, can find out, still can effectively identify after the point of release being positioned at weak convective region adopts sparse layouting, and the degree of correlation is still higher.And under this sparse condition of layouting, although the point of release in strong convection district also can identify preferably, but the relevance degree of institute's identification point is obviously lower, and the reliability of identification, than the weak effect of weak convective region, can find out that this points distributing method can carry out effective identification to the pollution source of weak convective region.
The weak convective region of table 4 adopts the recognition result after sparse layouting

Claims (2)

1. a Pollution Source Monitoring points distributing method, is characterized in that, the method comprises the following steps:
For interior space region, the mesh Reynolds number in flow field is obtained by the simulation of CFD method, according to the height of mesh Reynolds number, flow field, space is divided into strong convection district and two, weak convective region characteristic area, different principles of layouting is adopted for each characteristic area, namely adopt encryption to layout to strong convection district, sparse layouting is adopted to weak convective region;
The mesh Reynolds number maximal value in flow field, described space is designated as Re max, then get
R 0 = 1 100 Re max
To Re < R 0region, be defined as weak convective region, to Re > R 0region, be defined as strong convection district.
2. Pollution Source Monitoring points distributing method according to claim 1, is characterized in that, to strong convection district: according to delimited region, is become by this Region dividing physical dimension, shape close to consistent little space cell Pa, wherein Pa < P; To weak convective region: become by this Region dividing physical dimension, shape close to consistent little space cell Pb, Pb > P, wherein interior space region is evenly divided into physical dimension, little space cell that shape is consistent by P.
CN201210479898.6A 2012-11-22 2012-11-22 A kind of Pollution Source Monitoring points distributing method Expired - Fee Related CN102967689B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210479898.6A CN102967689B (en) 2012-11-22 2012-11-22 A kind of Pollution Source Monitoring points distributing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210479898.6A CN102967689B (en) 2012-11-22 2012-11-22 A kind of Pollution Source Monitoring points distributing method

Publications (2)

Publication Number Publication Date
CN102967689A CN102967689A (en) 2013-03-13
CN102967689B true CN102967689B (en) 2015-08-19

Family

ID=47797975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210479898.6A Expired - Fee Related CN102967689B (en) 2012-11-22 2012-11-22 A kind of Pollution Source Monitoring points distributing method

Country Status (1)

Country Link
CN (1) CN102967689B (en)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104697908A (en) * 2015-03-12 2015-06-10 国家海洋局天津海水淡化与综合利用研究所 Method for monitoring drifting salt deposition of seawater cooling tower
CN105528753B (en) * 2015-12-14 2018-04-20 北京邮电大学 A kind of air pollution source tracing method based on Network Virus Propagation model
JP6870442B2 (en) * 2017-04-06 2021-05-12 富士電機株式会社 Source analyzer and source analysis system
CN107315047A (en) * 2017-05-26 2017-11-03 内蒙古科技大学 A kind of method for reviewing pollution source of water body
CN108364103B (en) * 2018-02-23 2023-04-18 山东农业大学 Arrangement method of orchard sensor
CN109885804A (en) * 2019-01-23 2019-06-14 大连理工大学 A kind of air monitoring and source discrimination method based on monitoring car
CN110007054A (en) * 2019-04-30 2019-07-12 华北电力大学 A kind of coke-oven plant's polluter uncontrollable discharge source position determining device and method
CN111912937B (en) * 2020-07-20 2022-02-11 清华大学 Intelligent real-time tracing method and system for indoor sudden pollutants
CN112345698B (en) * 2020-10-30 2022-04-12 大连理工大学 Gridding arrangement method for air pollutant monitoring sites
CN112417553A (en) * 2020-11-16 2021-02-26 东南大学 Limited sensor arrangement strategy and intelligent ventilation system of industrial building
CN112432317A (en) * 2020-11-16 2021-03-02 东南大学 Sensor optimal arrangement method for classroom and ventilation monitoring system thereof
CN112712169A (en) * 2021-01-08 2021-04-27 中国科学院地理科学与资源研究所 Model building method and application of full residual depth network based on graph convolution
CN113627246B (en) * 2021-07-01 2024-04-16 同济大学 Building environment sudden pollution source identification method combining visual detection and concentration sensing
CN114359002B (en) * 2022-03-21 2022-05-20 四川国蓝中天环境科技集团有限公司 Atmospheric pollution small-scale tracing method and system based on mobile monitoring trend mining
CN117129638B (en) * 2023-10-26 2024-01-12 江西怡杉环保股份有限公司 Regional air environment quality monitoring method and system
CN117708617B (en) * 2024-02-05 2024-04-30 四川国蓝中天环境科技集团有限公司 Atmospheric pollution tracing method based on multi-source big data and pollution characteristic space-time matching

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6726406B2 (en) * 1999-10-28 2004-04-27 Battelle Memorial Institute In situ formation of reactive barriers for pollution control
US6591702B2 (en) * 2000-12-04 2003-07-15 Gas Technology Institute Method for identifying sources of rapidly released contaminants at contaminated sites
CN1834656A (en) * 2006-03-20 2006-09-20 天津大学 Monitoring point return selection method for deposit sediment of river channel
CN201368845Y (en) * 2008-12-12 2009-12-23 廖德兵 Vehicle mobile water quality organic pollutant surveillance analysis device
CN101566569B (en) * 2009-04-30 2011-09-14 西南科技大学 System and method for identifying a plurality of fluorescence spectrum mixed materials through characteristic parameter

Also Published As

Publication number Publication date
CN102967689A (en) 2013-03-13

Similar Documents

Publication Publication Date Title
CN102967689B (en) A kind of Pollution Source Monitoring points distributing method
CN103914622B (en) A kind of chemical leakage fast prediction alarm emergency response decision-making method
Jin et al. Simulating natural ventilation in and around buildings by fast fluid dynamics
CN203533216U (en) Real-time positioning analysis system of gas pipe network leak sources
CN105868444B (en) A kind of gas ductwork leakage detection and location method based on characteristic image
Yi et al. Assessing effects of wind speed and wind direction on discharge coefficient of sidewall opening in a dairy building model–A numerical study
CN106023561A (en) &#39;Four-meter-in-one&#39; energy management integrated service system
CN111912937B (en) Intelligent real-time tracing method and system for indoor sudden pollutants
Freire et al. On the improvement of natural ventilation models
CN104091001A (en) Three-dimensional simulation method and emergency evacuation method for gas leakage and diffusion
Papakonstantinou et al. Computational analysis of thermal comfort: the case of the archaeological museum of Athens
Zhang et al. Convective heat transfer prediction in large rectangular cross-sectional area Earth-to-Air Heat Exchangers
Zhai et al. Identifying decaying contaminant source location in building HVAC system using the adjoint probability method
Shao et al. Long-term prediction of dynamic distribution of passive contaminant in complex recirculating ventilation system
Chang et al. Computational fluid dynamics simulation of air exhaust dispersion from negative isolation wards of hospitals
Cui et al. Numerical study of the fire-smoke temperature law in the shaft of a high-rise building under the chimney effect in winter
Basok et al. Smart geothermal ventilation system
Chow et al. Scale modeling on natural smoke filling in an atrium
Jayaraman et al. Tracer gas transport under mixed convection conditions in an experimental atrium: Comparison between experiments and CFD predictions
Dama et al. Wind and buoyancy driven natural ventilation in double skin façades
Oladokun et al. Indoor mould growth prediction using coupled computational fluid dynamics and mould growth model
Benhamza et al. An experimental and a numerical study of horizontal earth-air heat exchanger in a hot climate
Zheng et al. Study on indoor air movement under the crowding conditions by numerical simulations and PIV small-scale experiments
Popescu et al. Particle penetration research through buildings’ cracks
Ahmed et al. Efficient Sensor Location for HVAC Systems Using Lagrangian Coherent Structures

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150819

Termination date: 20211122