一种利用固定位置污染物探测器辨识城市空间多污染物源的方法技术

技术编号:15299668 阅读:66 留言:0更新日期:2017-05-12 01:50
本发明专利技术提供了一种利用固定位置污染物探测器辨识城市空间多污染物源的方法,属于大气污染物源辨识技术领域。通过城市三维地图建立相应建筑的三维模型,根据气象站提供的城区上风向的主导风速和风向模拟计算出城区的速度场;选择三个污染物探测器读取器位置和污染物浓度数据,通过使用伴随方法计算出可能的污染物源位置和释放强度;对于多污染源情况按照一定的原则继续选择其他探测器,根据每三个污染物探测器的数据辨识污染物源,重复二至三次即可辨识到一个污染物源。然后依据流场稳态情况下污染物传播方程为线性方程这一理论基础将已辨识到的污染物源对其他污染物源辨识过程中产生的影响去除。如此可以快速、准确地辨识出多个污染物源。

Method for identifying multiple pollutant sources in urban space using fixed position pollutant detector

The invention provides a method for identifying a plurality of pollutant sources in a city space by using a fixed position pollutant detector, belonging to the technical field of air pollutant source identification. The establishment of the corresponding 3D model building through the three-dimensional map of the city, according to the meteorological station provided by the city dominated wind direction and wind direction simulation to calculate the velocity field of the city; three pollutant detector reader position and pollutant concentration data, by using the adjoint method to calculate the possible location of pollutant source and release the source of strength; according to certain principles to choose other detectors, according to the source data to identify each of the three pollutants pollutants detector, repeated two to three times can be identified to a pollutant source. Then, according to the equation of pollutant transmission in the steady state of the flow field, the theoretical basis is used to remove the influence of the identified pollutant source on the identification process of other pollutant sources. Thus, a number of pollutant sources can be identified rapidly and accurately.

【技术实现步骤摘要】

本专利技术属于大气污染物源辨识
,具体为一种利用固定位置污染物探测器辨识城市空间多污染物源的方法,用于对同时存在的多个污染物源进行位置和释放强度的辨识,从而对污染物进行源头控制。
技术介绍
现有的污染物源辨识方法主要分为两大类:正向计算方法和逆向计算方法。正向计算方法是一种相对比较简单的试错方法,比如GorelickSM等人在文章Identifyingsourcesofgroundwaterpollution:anoptimizationapproach中通过可移动探测器的不断尝试性移动逐渐趋近于污染物源的位置,这样一种方法往往效率比较低。而逆向计算方法则从探测器探测到污染物这一最终状态出发,通过逆向时间和速度场的模拟计算快速识别出污染物源。正向方法往往只能慢慢趋近于污染物源位置,但逆向方法在得到足够多的信息时,可以通过模拟计算快速确定污染源的位置、释放强度甚至污染物释放规律。Skaggs和Kabala在文章“Recoveringthehistoryofagroundwatercontaminantplume:Methodofquasi-reversibility”中利用准可逆方法,通过求解对流扩散方程实现了地下水污染源的辨识,并确定了其准可逆扩散因子求解方法。Zhang和Chen在文章“IdentificationofcontaminantsourcesinenclosedenvironmentsbyinverseCFDmodeling”中将此方法与计算流体力学结合,成功应用于飞机客舱污染物源和污染物释放强度的辨识。但是此方法的一个严重缺陷是需要提前知道污染物源的部分信息,这在实际情况中往往是难以实现的。Wagner在其文章“Simultaneouslyparameterestimationandcontaminantsourcecharacterizationforcoupledgroundwaterflowandcontaminanttransportmodeling”中提出了一套基于概率论的逆向模拟计算方法并将其应用于地下水污染物源的辨识。此方法将地下水流模拟计算与污染物传播模拟计算相结合,利用非线性最大相似理论对污染物源的位置和释放强度进行了逆向辨识。Neupauer和Wilson在其一系列文章“Backwardprobabilisticmodelofgroundwatercontaminationinnon-uniformandtransientflow”等中则将伴随方法引入地下水污染物源位置和释放时间的辨识。此方法能够准确推测出复杂几何环境下地下水污染物的历史特性,即其释放强度、位置和变化规律。Lin在文章“Inversemodelingmethodsforindoorairbornepollutanttracking:literaturereviewandfundamentals”中则通过引入概率分析进一步提高了其计算精确度并将此方法应用于一个实际地下水污染物源辨识的算例中。Liu和Zhai等人通过发表一系列文章,如“Locationidentificationforindoorinstantaneouspointcontaminantsourcebyprobability-basedinversecomputationalfluiddynamicsmodeling”等,则将此理论和方法应用于空气环境,并实现了室内空气污染物源的辨识。其方法能够根据污染物探测器提供的有限的信息实现室内空气污染物源位置和释放强度的快速辨识。前述大量方法仅仅是能成功实现单个污染物源的辨识,而现实中污染物源往往不止一个。目前在多个污染物源辨识方面取得的进展则寥寥可数。具有代表性的是Cai等人在文章Fastidentificationofmultipleindoorconstantcontaminantsourcesbyidealsensors:atheoreticalmodelandnumericalvalidation中使用正向方法进行多污染物源的辨识,结果表明此方法需要耗费大量的计算资源并且要预先知道可能的污染物源位置等很多实际情况下很难得到的数据。本专利技术将基于概率论的伴随方法应用于室外污染物源的辨识,并发展出一整套辨识室外多污染物源的方法。此方法借助有限个固定位置污染物探测器快速确定多个污染物源的具体位置和释放强度。
技术实现思路
城市污染物源的快速辨识有助于城市空气的治理和改善。传统试错法效率低且难以应对多污染源的情况,一般逆向方法需要预知污染源的部分信息,实际情况下难以实现。基于概率的伴随方法能够在不知道污染源任何信息的前提下仅凭探测器探测到的污染物浓度信息实现室内单个污染物源的辨识。本研究进一步改进了伴随方法用于开放空间的多污染物源辨识,并发展了能够仅仅根据有限的固定位置污染物探测器信息辨识多污染物源的方法。本专利技术的技术方案:一种利用固定位置污染物探测器辨识城市空间多污染物源的方法,步骤如下:(1)通过城市三维地图建立相应建筑的三维模型,根据气象站提供的城区上风向的主导风速和风向模拟计算出城区的速度场;(2)默认待测城区只有一个污染物源,在所有污染物探测器中随机挑选三个探测器,所选的探测器确保各个探测器检测到的污染物浓度值差别在20%以上,且探测器的位置不在一条直线上;然后,读取三个探测器检测到的污染物浓度信息、所处的位置坐标以及此时城市气象站记录的主导风向和风速,将数据分别记录为L1、L2和L3;使用从气象站得到的城市主导风向和风速数据模拟城市空间流场,根据L1、L2和L3的数据,辨识到第一个污染物源的第一个可能位置S1(1)和释放强度C1(1),并在研究区域的地图上标示出S1(1)位置(如附图4所示);(3)在地图上从位置S1(1)出发沿着主导风向的方向移动再寻找到三个位置最接近S1(1)的污染物探测器,记录其位置和污染物浓度信息并重复步骤(2)从而辨识到第一个污染物源的第二个可能位置S1(2)及相应的释放强度C1(2),并将S1(2)标注在地图上(如附图5所示)。寻找三个新的探测器时,必须向着一个方向偏移,这样才能使得到的污染物浓度值尽量只来自于一个污染源的影响。(4)如果离S1(1)最近的污染物探测器与离S1(2)最近的污染物探测器不是同一个污染物探测器,则表明没有辨识到污染物源,重复步骤(3)直至离S1(n)与S1(n+1)直至最近的污染物探测器为同一个探测器(n=1,2,3,…)。当离S1(n)与S1(n+1)最近的污染物探测器为同一个探测器时,我们就可以认为S1(n+1)和C1(n+1)即为第一个污染物源的位置和释放强度。(5)当第一个污染物源的位置和释放强度确定之后,根据S1和C1求解公式1-6,从而模拟计算得到由第一个污染源导致的整个区域的污染物浓度分布(如附图6所示),读取各个污染物探测器对应位置处的污染物浓度值。由于在速度场稳定的情况下,公式1-6相对于质量浓度C来说是一个线性方程,因此各个污染物源造成的污染物浓度场可以直接相加或相减。将原有的所有探测器探测到的污染物浓度数据减去由第一个污染物源导致的浓度值,即得到剩余其它污染物源导致的各个污染物探测器处的污染物本文档来自技高网
...
一种利用固定位置污染物探测器辨识城市空间多污染物源的方法

【技术保护点】
一种利用固定位置污染物探测器辨识城市空间多污染物源的方法,其特征在于,步骤如下:第一步,根据城市三维地图建立建筑的三维模型,再根据气象站提供的待测城区上风向的主导风速和风向,使用计算流体力学求解纳维斯托克斯方程得出待测城区的速度场;第二步,首先,默认待测城区只有一个污染物源,在所有污染物探测器中随机挑选三个探测器,所选的探测器确保各个探测器检测到的污染物浓度值差别在20%以上,且探测器的位置不在一条直线上;然后,读取三个探测器检测到的污染物浓度信息、所处的位置坐标以及此时城市气象站记录的主导风向和风速,将数据分别记录为L1、L2和L3;对每个探测器探测到的信息,通过求解污染物传播方程的伴随方程:其中,ψ*为伴随概率因子,位置或时间的伴随概率因子;τ为逆向的时间;为探测区域位置矢量;为污染物探测器位置矢量;C表示污染物浓度;Xj为位置坐标轴,其中j=1,2,3分别对应水平、垂直、竖直三个方向;Vj为Xj轴方向上的速度;νC,j表示污染物C在Xj方向上的有效湍流扩散系数;qO为污染物负源的单位体积流量;Γ1,Γ2和Γ3为边界条件;ni为Xj轴方向的单位矢量;为负荷项,其表达式由两个阶跃方程组成:δ(τ)=1,τ=00,τ≠0---(1-3)]]> 得出各个探测器辨识到的污染物源可能存在的位置,上述得出的可能的污染源位置无限多,通过公式(1‑4)求解如下方程将三个探测器的辨识结果整合,即确定唯一一个可能的污染物源:fx(x|C^1,...,C^N;τ0,x→w,1,...,x→w,N,τ1,...,τN)=∫M0Πi=1NP(C^i|M0,x;τ0,x→w,i,τi)fx(x;τ0,x→w,i,τi)dM0∫M0∫xΠi=1NP(C^i|M0,x;τ0,x→w,i,τi)fx(x;τ0,x→w,i,τi)dxdM0---(1-4)]]>其中,N为探测数据的个数;τi、和分别为对应于第i个探测数据的探测位置、探测时间和探测到的污染物浓度;τ0为已知的污染物释放时间;M0为假设的污染物释放强度;为根据第i个探测数据通过公式(1‑1)计算出的污染源概率分布;为根据第i个探测数据求得的相应污染物释放浓度M0和位置x的概率分布;一般将的分布形式定义为正态分布:P(C^i|M0,x;τ0,x→w,i,τi)~N(M0·fx(x;τ0,x→w,i,τi),σϵ2)---(1-5)]]>其中,为对应于第i个探测数据可能的实际污染物浓度;为污染物探测器的测量误差的标准平方差;通过将L1、L2和L3对应的信息代入上述公式(1‑1)至公式(1‑5)求解,辨识到第一个污染物源的第一个可能位置S1(1)和释放强度C1(1)后,在待测城区的地图上标示出S1(1)位置;第三步,在地图上从位置S1(1)出发沿着主导风向的方向移动再寻找三个位置最接近S1(1)的污染物探测器,记录其位置和污染物浓度信息并重复第二步,辨识到第一个污染物源的第二个可能位置S1(2)及相应的释放强度C1(2),并在待测城区的地图上标示出S1(2)位置;寻找三个新的探测器时,必须向同一个方向偏移,使得到的污染物浓度值尽量只来自于一个污染源的影响;第四步,如果离S1(1)最近的污染物探测器与离S1(2)最近的污染物探测器不是同一个污染物探测器,则表明没有辨识到污染物源,重复第三步直至离S1(n)与S1(n+1)最近的污染物探测器为同一个探测器,n=1,2,3,…;当离S1(n)与S1(n+1)最近的污染物探测器为同一个探测器时,S1(n+1)和C1(n+1)即为第一个污染物源的位置和释放强度;第五步,确定第一个污染物源的位置和释放强度后,根据S1(n+1)和C1(n+1)求解污染物状态方程其中,C表示污染物浓度;t为时间;Vj为Xj轴方向上的速度;νC,j表示污染物C在Xj方向上的有效湍流扩散系数;qO为污染物负源的单位体积流量;qI为污染物正源的单位体积流量;CI为相应的污染物体积浓度;SC为所有其它形式的污染物源,则(Sc+qICI‑qOC)为所有外部污染物源的总和;Γ1,Γ2和Γ3为边界条件;C0为各个位置的初始浓度;g1,g2和g3表示已知的污染物边界条件数学表达式;ni为Xj轴方向的单位矢量;得到由第一个污染物源导致的整个区域的污染...

【技术特征摘要】
1.一种利用固定位置污染物探测器辨识城市空间多污染物源的方法,其特征在于,步骤如下:第一步,根据城市三维地图建立建筑的三维模型,再根据气象站提供的待测城区上风向的主导风速和风向,使用计算流体力学求解纳维斯托克斯方程得出待测城区的速度场;第二步,首先,默认待测城区只有一个污染物源,在所有污染物探测器中随机挑选三个探测器,所选的探测器确保各个探测器检测到的污染物浓度值差别在20%以上,且探测器的位置不在一条直线上;然后,读取三个探测器检测到的污染物浓度信息、所处的位置坐标以及此时城市气象站记录的主导风向和风速,将数据分别记录为L1、L2和L3;对每个探测器探测到的信息,通过求解污染物传播方程的伴随方程:其中,ψ*为伴随概率因子,位置或时间的伴随概率因子;τ为逆向的时间;为探测区域位置矢量;为污染物探测器位置矢量;C表示污染物浓度;Xj为位置坐标轴,其中j=1,2,3分别对应水平、垂直、竖直三个方向;Vj为Xj轴方向上的速度;νC,j表示污染物C在Xj方向上的有效湍流扩散系数;qO为污染物负源的单位体积流量;Γ1,Γ2和Γ3为边界条件;ni为Xj轴方向的单位矢量;为负荷项,其表达式由两个阶跃方程组成:δ(τ)=1,τ=00,τ≠0---(1-3)]]>得出各个探测器辨识到的污染物源可能存在的位置,上述得出的可能的污染源位置无限多,通过公式(1-4)求解如下方程将三个探测器的辨识结果整合,即确定唯一一个可能的污染物源:fx(x|C^1,...,C^N;τ0,x→w,1,...,x→w,N,τ1,...,τN)=∫M0Πi=1NP(C^i|M0,x;τ0,x→w,i,τi)fx(x;τ0,x→w,i,τi)dM0∫M0∫xΠi=1NP(C^i|M0,x;τ0,x→w,i,τi)fx(x;τ0,x→w,i,τi)dxdM0---(1-4)]]>其中,N为探测数据的个数;τi、和分别为对应于第i个探测数据的探测位置、探测时间和探测到的污染物浓度;τ0为已知的污染物释放时间;M0为假设的污染物释放强度;为根据第i个探测数据通过公式(1-1)计算出的污染源概率分布;为根据第i个探测数据求得的相应污染物释放浓度M0和位置x的概率分布;一般将的分布形式定...

【专利技术属性】
技术研发人员:翟志强薛雨
申请(专利权)人:大连理工大学
类型:发明
国别省市:辽宁;21

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1