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

本发明属于大气污染物源辨识技术领域，具体为一种利用固定位置污染物探测器辨识城市空间多污染物源的方法，用于对同时存在的多个污染物源进行位置和释放强度的辨识，从而对污染物进行源头控制。

背景技术：

现有的污染物源辨识方法主要分为两大类：正向计算方法和逆向计算方法。正向计算方法是一种相对比较简单的试错方法，比如Gorelick S M等人在文章Identifying sources of groundwater pollution:an optimization approach中通过可移动探测器的不断尝试性移动逐渐趋近于污染物源的位置，这样一种方法往往效率比较低。而逆向计算方法则从探测器探测到污染物这一最终状态出发，通过逆向时间和速度场的模拟计算快速识别出污染物源。正向方法往往只能慢慢趋近于污染物源位置，但逆向方法在得到足够多的信息时，可以通过模拟计算快速确定污染源的位置、释放强度甚至污染物释放规律。Skaggs和Kabala在文章“Recovering the history of a groundwater contaminant plume:Method of quasi-reversibility”中利用准可逆方法，通过求解对流扩散方程实现了地下水污染源的辨识，并确定了其准可逆扩散因子求解方法。Zhang和Chen在文章“Identification of contaminant sources in enclosed environments by inverse CFD modeling”中将此方法与计算流体力学结合，成功应用于飞机客舱污染物源和污染物释放强度的辨识。但是此方法的一个严重缺陷是需要提前知道污染物源的部分信息，这在实际情况中往往是难以实现的。Wagner在其文章“Simultaneously parameter estimation and contaminant source characterization for coupled groundwater flow and contaminant transport modeling”中提出了一套基于概率论的逆向模拟计算方法并将其应用于地下水污染物源的辨识。此方法将地下水流模拟计算与污染物传播模拟计算相结合，利用非线性最大相似理论对污染物源的位置和释放强度进行了逆向辨识。Neupauer和Wilson在其一系列文章“Backward probabilistic model of groundwater contamination in non-uniform and transient flow”等中则将伴随方法引入地下水污染物源位置和释放时间的辨识。此方法能够准确推测出复杂几何环境下地下水污染物的历史特性，即其释放强度、位置和变化规律。Lin在文章“Inverse modeling methods for indoor airborne pollutant tracking:literature review and fundamentals”中则通过引入概率分析进一步提高了其计算精确度并将此方法应用于一个实际地下水污染物源辨识的算例中。Liu和Zhai等人通过发表一系列文章，如“Location identification for indoor instantaneous point contaminant source by probability-based inverse computational fluid dynamics modeling”等，则将此理论和方法应用于空气环境，并实现了室内空气污染物源的辨识。其方法能够根据污染物探测器提供的有限的信息实现室内空气污染物源位置和释放强度的快速辨识。前述大量方法仅仅是能成功实现单个污染物源的辨识，而现实中污染物源往往不止一个。目前在多个污染物源辨识方面取得的进展则寥寥可数。具有代表性的是Cai等人在文章Fast identification of multiple indoor constant contaminant sources by ideal sensors:a theoretical model and numerical validation中使用正向方法进行多污染物源的辨识，结果表明此方法需要耗费大量的计算资源并且要预先知道可能的污染物源位置等很多实际情况下很难得到的数据。

本发明将基于概率论的伴随方法应用于室外污染物源的辨识，并发展出一整套辨识室外多污染物源的方法。此方法借助有限个固定位置污染物探测器快速确定多个污染物源的具体位置和释放强度。

技术实现要素：

城市污染物源的快速辨识有助于城市空气的治理和改善。传统试错法效率低且难以应对多污染源的情况，一般逆向方法需要预知污染源的部分信息，实际情况下难以实现。基于概率的伴随方法能够在不知道污染源任何信息的前提下仅凭探测器探测到的污染物浓度信息实现室内单个污染物源的辨识。本研究进一步改进了伴随方法用于开放空间的多污染物源辨识，并发展了能够仅仅根据有限的固定位置污染物探测器信息辨识多污染物源的方法。

本发明的技术方案：

一种利用固定位置污染物探测器辨识城市空间多污染物源的方法，步骤如下：

(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来说是一个线性方程，因此各个污染物源造成的污染物浓度场可以直接相加或相减。将原有的所有探测器探测到的污染物浓度数据减去由第一个污染物源导致的浓度值，即得到剩余其它污染物源导致的各个污染物探测器处的污染物浓度值。

(6)重复第二步至第五步直至所有的污染物探测器的数值都低于污染物对人体造成危害的限值。

为了验证结果的准确性，可以在完成以上步骤辨识到所有污染物源之后，再挑选一个主导风速和风向与原来明显不同的时间段进行重复辨识，以确保所有污染物源都已经被找到。

本发明的有益效果：本发明首次提出有别于传统试错法的多污染物源的快速逆向辨识方法，能够利用有限个固定位置污染物探测器快速辨识出多个污染物源，可直接采用各个城市现有的污染物监测站的数据，初始投资低，速度快，结果准确，能够实现对城市大气环境的实时监测及快速诊断，为城市大气污染的治理和决策提供经济高效的技术方案。为污染物监测站的选址提供指导意见。

附图说明

图1是某城区建筑群三维模型。

图2是城区在主导风速和风向作用下速度场分布。

图3是使用固定位置污染物探测器辨识城市空间多污染物源流程图。

图4是第一次尝试辨识首个污染物源示意图。

图5是第二次尝试辨识首个污染物源示意图。

图6是辨识到的第一个污染物源导致的污染物浓度场分布图。

具体实施方式

以下结合附图和技术方案，进一步说明本发明的具体实施方式。

此发明所使用的方法为逆向方法，它需要将时间逆向考虑，这就决定了此方法需要与模拟计算相结合。使用逆向方法进行城市空间气态污染物源的辨识则需要借助计算流体力学对城市空间的流场和污染物浓度场进行模拟计算。伴随方法通过求解质量传递方程的伴随方程，然后依据测量数据，结合概率理论求取污染物源位置和释放强度在整个空间区域和释放强度坐标上的概率分布，以得到的概率值最大的点所对应的位置和污染物释放强度为辨识出的污染物源的位置和强度。

实施例

一种利用固定位置污染物探测器辨识城市空间多污染物源的方法，步骤如下：

第一步，根据城市三维地图建立建筑的三维模型，再根据气象站提供的待测城区上风向的主导风速和风向，使用计算流体力学求解纳维斯托克斯方程得出待测城区的速度场；

第二步，首先，默认待测城区只有一个污染物源，在所有污染物探测器中随机挑选三个探测器(如附图4所示)，所选的探测器确保各个探测器检测到的污染物浓度值差别在20％以上，且探测器的位置不在一条直线上；然后，读取三个探测器检测到的污染物浓度信息、所处的位置坐标以及此时城市气象站记录的主导风向和风速，将数据分别记录为L1、L2和L3；

对每个探测器探测到的信息，通过求解污染物传播方程的伴随方程：

其中，ψ^*为伴随概率因子(位置或者时间的伴随概率因子)，τ为逆向的时间，为探测区域位置矢量，为污染物探测器位置矢量，C表示污染物浓度，X_j为位置坐标轴，其中j＝1,2,3分别对应水平、垂直、竖直三个方向，V_j为X_j轴方向上的速度，ν_C,j表示污染物C在X_j方向上的有效湍流扩散系数，q_O为污染物负源的单位体积流量，Γ₁,Γ₂和Γ₃为边界条件，n_i为X_j轴方向的单位矢量；为负荷项，其表达式由两个阶跃方程组成：

得出各个探测器辨识到的污染物源可能存在的位置，上述得出的可能的污染源位置无限多，通过公式(1-4)求解如下方程将三个探测器的辨识结果整合，即确定唯一一个可能的污染物源：

其中，N为探测数据的个数，τ_i、和分别为对应于第i个探测数据的探测位置、探测时间(逆向)和探测到的污染物浓度，τ₀为已知的污染物释放时间，M₀为假设的污染物释放强度，为根据第i个探测数据通过公式1-1计算出的污染源分布概率；为根据第i个探测数据求得的相应污染物释放浓度M₀和位置x的概率分布；一般将的分布形式定义为正态分布：

其中，为对应于第i个探测数据可能的实际污染物浓度，为污染物探测器的测量误差的标准平方差；

通过将L1、L2和L3对应的信息代入上述公式(1-1)至公式(1-5)求解，辨识到第一个污染物源的第一个可能位置S1(1)和释放强度C1(1)后，在待测城区的地图上标示出S1(1)位置(如附图4所示)；

第三步，在地图上从位置S1(1)出发沿着主导风向的方向移动再寻找三个位置最接近S1(1)的污染物探测器，记录其位置和污染物浓度信息并重复第二步，辨识到第一个污染物源的第二个可能位置S1(2)及相应的释放强度C1(2)，并在待测城区的地图上标示出S1(2)位置(如附图5所示)；寻找三个新的探测器时，必须向同一个方向偏移，使得到的污染物浓度值尽量只来自于一个污染源的影响；

第四步，如果离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为时间，V_j为X_j轴方向上的速度，ν_C,j表示污染物C在X_j方向上的有效湍流扩散系数，q_O为污染物负源的单位体积流量，q_I为污染物正源的单位体积流量，C_I为相应的污染物体积浓度，S_C为所有其它形式的污染物源，则(Sc+q_IC_I-q_OC)为所有外部污染物源的总和；Γ₁,Γ₂和Γ₃为边界条件，C₀为各个位置的初始浓度，g₁,g₂和g₃表示已知的污染物边界条件数学表达式，n_i为X_j轴方向的单位矢量；

得到由第一个污染物源导致的整个区域的污染物浓度分布(如附图6所示)，读取各个污染物探测器对应位置处的污染物浓度值；由于在速度场稳定的情况下，公式(1-6)相对于质量浓度C来说是一个线性方程，因此各个污染物源造成的污染物浓度场可以直接相加或相减；将原有的所有探测器探测到的污染物浓度数据减去由第一个污染物源导致的浓度值，即得到剩余其它污染物源导致的各个污染物探测器处的污染物浓度值；

第六步，重复第二步至第五步直至所有的污染物探测器的数值都低于污染物对人体造成危害的限值。

本方法适用于如下特定的情境：

(1)流场是稳态的。在进行城市空间污染物源辨识时，本研究将直接使用气象站测量到的气象数据来模拟计算城市空间的流场。气象数据显示城市空间的主导风速和主导风向可能维持数小时不变，我们可以选取这些时段的数据进行污染物源的辨识。

(2)污染物源是释放强度恒定的点源。基于概率的伴随方法只能逆向辨识点源型(或者可以作为点源来考虑)的污染物源，线源和面源不在本研究的讨论范围之内。

(3)污染物是惰性污染物，且气流跟随性较好。为了方便起见，本研究只针对气流跟随性较好的惰性污染物。而如果要进一步考虑可与大气中其它物质反应或气流跟随性较差的颗粒性污染物，只要能够模拟计算准确，此方法同样适用。

(4)污染物探测器能够探测并记录污染物浓度，并且所有探测器的位置都已在地图中标示出来。

(5)有超过三个位置固定的污染物探测器。本方法需要至少三个污染物信息来对污染物源的位置和释放强度进行辨识。之所以需要三个，是因为此方法所要辨识的污染物源的信息包括位置(经度坐标和纬度坐标)以及释放强度三个变量，需要三个数据信息才能确定这三个变量的值。

(6)污染物浓度达到一定的限值才认为有害，才必须搜寻到相应的污染源。