技术领域
本发明属于大气污染监测和源辨识领域,具体涉及一种动态风场条件下的空气污染源逆向定位方法。
背景技术
大气污染已成为现代城市面临的严重问题,污染物的不合理排放不但对环境造成影响,同时威胁着居民的健康。在大气污染事件中,在有限的监测条件下对污染源的位置进行准确识别,对于城市大气污染源的控制管理以及改善城市空气质量意义重大。目前对于大气污染的监测方法不够严谨,监测水平亟待提升,同时缺少应用于实际的有效寻源方式。
现有专利:一种利用固定污染物探测器辨识城市空间多污染物源的方法(授权公告号CN106650017A)采用基于概率论的伴随方法,在稳定流场中通过多个固定位置探测器做到对城市多个污染源的快速辨识。该发明的缺陷在于:需要至少三个固定监测站点提供监测信息,没有考虑到现实情况中由于经济因素或其他条件限制,待测区域不会有众多符合条件的监测站点以供选择,有时一片区域往往只能设置一个固定监测站点,此时上述方法失效。
因此,针对上述情况,本发明进一步丰富了伴随概率方法在城市空间污染物源辨识领域的应用,提出一种利用一个固定污染物监测站在多个气象条件下获得的监测数据来辨识城市空间单个污染源的方法,在保证辨识准确性的条件下,将监测设备进行简化,极大地减少了城市大气污染治理的资金投入,有助于城市空气的治理和改善。
发明内容
本发明的主要目的在于解决城市污染源的快速辨识问题,以及解决专利(授权公告号CN106650017A)没有考虑到待测区域固定监测站的缺陷。提出一种当监测条件有限时,借助单个污染物监测站辨识污染源的方法。
为解决以上问题,本发明提供了一种利用一个固定污染物监测站在多个气象条件下获得的监测数据来辨识城市空间单个污染源的方法,一种动态风场条件下的空气污染源逆向定位方法,步骤如下:
第一步,根据城市卫星地图建立待测区域建筑的三维模型。
第二步,若在一时间段内监测到一固定监测站某种污染物浓度值高于污染物对人体造成危害的限值,将该时段测得的污染物平均浓度值记录为C1,将气象站坐标记录为L;选取气象条件有明显差异的另外两个时段,将这两个时间段该固定监测站监测到的大于空气中本底浓度的该种污染物数据分别记录为C2和C3,要求每个时段内的主导风向和风速较为稳定。
第三步,从气象站处获取待测区域在时段1的主导风向和风速信息,作为速度入口边界条件,使用计算流体力学求解纳维斯托克斯方程,得出待测区域在气象条件1下的流场1。由于每个时段内的主导风向和风速较为稳定,认为每个时段内流场是稳态的。
第四步,在流场1的计算基础上,将监测站坐标和气象条件1的污染物平均浓度信息(L,C1)代入稳态的污染物传播方程的伴随方程:
其中ψ*为位置的伴随概率因子,
通过求解方程得到该气象条件下推得的污染源的潜在位置概率分布,概率最大位置就是污染物源最可能存在的位置。根据一个测点信息在单一气象条件下推得的污染源位置是无限多的。
第五步,重复步骤三至步骤四,以获得另外两个气象条件下的定位结果。
第六步,通过求解如下方程(1-3),将三次辨识结果整合,即可确定唯一一个可能的污染源位置及强度:
其中N为不同气象条件下的辨识次数,
其中
将三个气象条件下的辨识结果代入上述公式(1-3)至公式(1-4),即可确定污染物源唯一的位置P-find和释放强度S-find,并估算出P-find处的浓度C-find。
第七步,前往P-find处确定污染源的真实浓度C-real,若浓度C-real与C-find在一个量级,则寻源成功,定位结束。
以上步骤的排序根据实际情况可做略微调整,调整顺序后的寻源过程也应在本专利保护范围内。
第二步所述的三个时段的气象条件有明显差异,此处的明显差异指主导风向相差至少22.5°。应注意的是,此处提出的风向角度差异是一种优选的方法,若选用其他风向角度差异或者风速大小差异也应在本专利保护范围内。
第二步至第六步所述的记录三个时段的气象参数及污染数据并以此通过伴随概率方法进行源定位。若选用其他数量的时段相应的数据通过本方法进行溯源也应在本专利保护范围内。
所述大气污染物监测站设有仪器可满足对待测地区特定的污染物种类的监测,且保证测试精度。这里检测站内的仪器设置是优选的方案,其他同类型测试仪器的设置方案也在保护范围内。
方法中所依赖的硬件不限于固定监测站点,也可以是其他有相同监测和记录功能的移动装置,受外界条件限制停在某处作为固定测点使用,或是人员携带同等监测和记录功能的便携式移动设备,于某固定位置进行操作。这里固定监测站只是一个优选方案,其他形式的监测设备也在本专利保护范围内。
与现有技术相比,本发明的有益效果是:
本发明计算时不再需要待测区域设有至少三个的固定监测站点,而是仅仅借助一个固定监测站点即可得到同样的辨识效果。在保证辨识准确性的条件下,将设备数量进行简化,极大地降低了监测成本。
附图说明
图1为本发明提供的一种大气污染寻源方法的流程示意图。
图2为本发明实施例提供的监测站L1于气象条件1下的定位结果图。
图3为本发明实施例提供的监测站L1于气象条件2下的定位结果图。
图4为本发明实施例提供的监测站L1于气象条件3下的定位结果图。
图5为本发明实施例提供的监测站L1综合定位结果图。
图6为本发明实施例提供的其他位置监测站的综合定位结果图。
具体实施方法
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
本发明的使用条件是:在气象条件有明显差异的三个时段内固定监测站均监测到大于空气中本底浓度的某种污染物。本发明需要借助计算流体力学对城市空间的流场和污染物场进行逆向模拟计算;再依据污染物监测站提供的浓度数据,使用伴随方法通过求解质量传递方程的伴随方程,求取污染物源位置和释放强度在整个空间区域内的概率分布;最后结合概率理论,综合三个气象条件下的概率分布结果,得到的概率值最大的点所对应的位置和强度为辨识出的污染物源的位置和强度。
以一个三维建筑群模型为实施例,本例中正向的气象数据及污染物浓度数据均由风洞实验获得,实验模型相对于实际建筑的缩小比例为1:500。寻源步骤如下:
第一步,根据风洞模型建立待测区域建筑的三维模型。(省略)
第二步,在时段1,坐标为L1的固定监测站监测到某种污染物浓度值高于污染物对人体造成危害的限值,记录该时段测得的污染物平均浓度值为50ppm,此时主导风向为西北风。在时段2,该监测站检测到污染物浓度位26ppm,此时主导风向为北风。在时段3,该监测站检测到污染物浓度位6ppm,此时主导风向为东北风。三个时段的主导风向相差至少30°。
第三步,从气象站处获取待测区域在时段1的主导风向为北偏西30°,风速3.87m/s。将其作为速度入口边界条件,使用计算流体力学求解纳维斯托克斯方程,得出待测区域在气象条件1下的流场1(参见图2-1)。由于时段1的主导风向和风速是稳定的,认为该时段内流场是稳态的。
第四步,在流场1的计算基础上,将监测站坐标和气象条件1的污染物平均浓度信息(L,C1)代入稳态的污染物传播方程的伴随方程:
其中ψ*为位置的伴随概率因子,
通过求解方程得到该气象条件下推得的污染源的潜在位置概率分布,概率最大位置就是污染物源最可能存在的位置(参见图2-2)。根据一个测点信息在单一气象条件下推得的污染源位置是无限多的。
第五步,从气象站处获取待测区域在时段2的主导风向为正北风,风速3.87m/s;时段3的主导风向为北偏东30°,风速3.87m/s。重复步骤三至步骤四,以获得另外两个气象条件下的流场(参见图3-1和图4-1)和定位结果(参见图3-2和图4-2)。
第六步,通过求解如下方程(1-3),将三次辨识结果整合,即可确定唯一一个可能的污染源位置及强度:
其中N为不同气象条件下的辨识次数,
其中
将三个气象条件下的辨识结果代入上述公式(1-3)至公式(1-4),即可求出污染物源唯一的位置P-find(参见图5)和释放强度275mL/min,并估算出P-find处的浓度C-find=10560ppm。
第七步,前往P-find处确定污染源的真实浓度C-real=7200ppm,浓度C-real与C-find在一个量级,故寻源成功,定位结束。
为证明本方法的适用性,选取其他位置的监测浓度符合条件的固定站点,均可在一定误差范围内定位到实际污染源(参见图6)。
本方法适用于有如下特定情境:
(1)污染物源是释放强度恒定的点源。基于概率的伴随方法只能逆向辨识点源型(或者可以作为点源来考虑)的污染物源,线源和面源不在本研究的讨论范围之内。
(2)污染物是惰性污染物,且气流跟随性较好。为了方便起见,本研究只针对气流跟随性较好的惰性污染物。而如果要进一步考虑可与大气中其它物质反应或气流跟随性较差的颗粒性污染物,只要能够模拟计算准确,此方法同样适用。
(3)在每个时段内流场是稳定的,或每个时段内主导风速风向变化不大,可当作稳态流场考虑。
(4)监测站能够监测到特定种类的污染物并作实时记录。
(5)污染物浓度达到一定的限值才认为有害,才必须搜寻到相应的污染源。
机译: 不利气象条件下空气污染源运行模式的控制方法
机译: 用于操作风场或风场,风场和风场设施的方法。
机译: 用相同方法黑启动风轮机风场并恢复风场和风轮机风场的方法