1.一种空气污染物观测站点布局评价方法,其特征在于,包括以下步骤:
步骤1、基于空间代表性计算方法,综合考虑污染物浓度是否达标的情况,计算出污染物观测站点空间代表范围;
步骤2、空气污染物观测站点分布冗余度的综合评估;
步骤3、城市内所有站点空气污染物空间覆盖比例综合评估;
步骤4、城市内所有站点空气污染物人口覆盖比例综合评估;
步骤5、城市内所有站点不同空气污染物空间代表性评估结果一致性;
步骤6、基于空气污染物观测站点分布冗余度、城市内所有站点空气污染物空间覆盖比例、城市内所有站点空气污染物人口覆盖比例及城市内所有站点不同空气污染物空间代表性评估结果一致性,构建空气污染物观测站点布局定量评价指标;
所述步骤1包括以下步骤:
步骤1.1、单个站点单个污染物的空间代表性评估;
给定一个观测站点i,判定另外一个点位(x,y)是否在站点i单个污染物AP的空间代表范围内,采用如下公式:采用如下公式:
单个站点单个污染物的空间代表性评估通过多个指标进行综合评估;首先,采用公式(1)和(2),进行基于点位污染物浓度相似度的单个站点单个污染物空间代表性评估,以站点i为中心,设置50公里×50公里窗口,计算窗口内每个点位的f1i,进而判断待判定点位与观测站点的污染物浓度是否相似;
其次,采用公式(3)和(4),进行基于空气质量日达标程度一致性的单个站点单个污染物空间代表性评估,上式中,i表示的是第i个观测站点,(Xi,Yi)表示的是站点位置,(x,y)表示的是50公里×50公里窗口中一个待判定点位的位置,Nt表示的是观测次数,d表示的是某一空气污染物AP的日达标值,Flag2表示的是tk时刻的该空气污染物AP浓度在上述两个点位是否日达标程度一致,如果两个点位日空气污染物浓度同时达标或者同时不达标,则认为两个点位日达标程度一致,Flag2在tk时刻记为1,f2i(x,y)最终表示的是基于日达标程度一致性的观测站点i空间代表性,该值如果大于或等于0.8,则认为点位(x,y)与站点i基于空气质量日达标程度一致,具有代表性;
采用公式(5),进行基于空气质量年达标程度一致性的单个站点单个污染物空间代表性评估,y表示的是某一空气污染物AP的年达标值,f3i(x,y)表示的是在点位(x,y)与观测站点(Xi,Yi)两处该空气污染物AP平均浓度是否年达标程度一致,如果两个点位年平均空气污染物浓度同时达标或者同时不达标,则认为两个点位年达标程度一致,f3i(x,y)在点位(x,y)处记为1,f3i(x,y)最终表示的是基于年达标程度一致性的单个站点i空间代表性,该值如果为1,则认为点位(x,y)与站点i基于空气质量年达标程度一致,具有代表性;
最后,将以上三个空间代表性范围做交集,得到单个站点单个污染物的综合空间代表范围FAPi(x,y);
步骤1.2、基于步骤1.1,计算城市内单个污染物观测站点的空间代表性;计算城市内其它站点的空间代表性,将这些不同站点的空间代表性做并集,即可计算出城市整体的单个污染物观测站点空间覆盖范围,进而得到城市整体的单个污染物观测站点的空间代表性;
所述步骤2包括以下步骤:
步骤2.1、计算城市内单个污染物观测站点分布的冗余度;计算单个污染物不同站点之间的冗余性指标RAP,旨在量化单个污染物不同站点空间代表范围的重复交叉程度,具体计算公式如下:式(6)中,RAP表示在给定一个城市中单个污染物观测站点的冗余度,其中AP为空气污染物名称,该值范围在0‑1,0表示没有冗余,1表示所有站点的空间代表范围都相同;SAP表示城市单个污染物整体空间代表性所覆盖的面积,SiAP表示由i个站点导致的单个污染物空间代表范围重复交叉面积;
步骤2.2、计算城市内不同空气污染物观测站点分布冗余度的综合指标;空气污染物种类数为Np,综合考虑不同污染物观测站点空间代表范围分布冗余度,得到空气污染物观测站点分布冗余度的综合指标,具体计算公式如下:式(7)中,空气污染物种类数为Np,其中p为大于等于1的正整数,表示第p种空气污染物;
计算一年观测天数Nt中不同污染物为首要大气污染物的天数,记不同空气污染物在观测天数Nt中为首要污染物的天数为tAP,RAP为不同污染物观测站点空间代表范围分布对应的冗余度;
所述步骤3包括以下步骤:
步骤3.1、计算城市内所有站点单个污染物空间覆盖比例;计算城市内所有站点单个污染物空间代表范围的覆盖比例,旨在量化单个污染物所有站点空间代表范围在整个城市的覆盖程度,具体计算公式如下:式(8)中,SAP表示城市内所有站点单个污染物整体空间代表性所覆盖的面积,其中AP为空气污染物名称,S_area表示该城市的总面积,Pro_areaAP表示城市所有站点单个污染物空间覆盖比例;
步骤3.2、计算城市内所有站点不同污染物空间覆盖比例的综合指标;空气污染物种类数为Np,综合考虑不同污染物所有观测站点空间覆盖比例,得到城市内所有站点空气污染物空间覆盖范围比例的综合指标,具体计算公式如下:将不同空气污染物在观测天数Nt中为首要污染物的天数tAP作为权重与各空气污染物的空间覆盖比例Pro_areaAP相乘进行加权求和,并除以权值总和,进而得到城市内所有站点不同空气污染物的空间覆盖比例综合指标Pro_area;
所述步骤4包括以下步骤:
步骤4.1、计算城市内所有站点单个污染物人口覆盖比例;计算了城市内所有站点单个污染物空间代表范围的人口覆盖比例,旨在量化单个污染物所有站点空间代表范围在整个城市的人口覆盖程度,具体计算公式如下:式(10)中,POPAP表示城市内所有站点单个污染物整体空间代表性所覆盖的人口总数,POP表示该城市的总人口数,Pro_popAP表示城市所有站点单个污染物空间代表范围人口覆盖比例;
步骤4.2、空气污染物种类数为Np,综合考虑不同污染物所有观测站点人口覆盖比例,得到城市内所有站点空气污染物人口覆盖比例的综合指标,具体计算公式如下:将不同空气污染物在观测天数Nt中为首要污染物的天数tAP作为权重与各空气污染物的人口覆盖比例Pro_popAP相乘进行加权求和,并除以权值总和,进而得到城市内所有站点不同空气污染物的人口覆盖比例综合指标Pro_pop;
所述步骤5包括以下步骤:
步骤5.1、计算单个站点不同污染物的空间代表性评估结果一致性;空气污染物种类数为Np,给定一个观测站点i,计算单个站点不同污染物的空间代表性评估结果一致性,具体计算公式如下:式(12)中,Ti表示第i个站点不同空气污染物空间代表性评估结果的一致性,该值范围在0‑1,0表示评估结果不一致,1表示单个站点的不同污染物空间代表范围评估结果完全一致;Si表示城市某一站点i不同污染物整体空间代表性所覆盖的面积,Spi表示在某一站点i,p种污染物的空间代表范围重复交叉面积;
步骤5.2、计算所有站点不同污染物的空间代表性评估结果一致性;综合考虑城市内各个站点不同污染物空间代表性评估结果一致性,得到所有站点不同污染物的空间代表性评估结果一致性,具体计算公式如下:式(13)中,T表示不同空气污染物空间代表性评估结果的一致性,Ni表示站点总数,Ti表示第i个站点不同空气污染物空间代表性评估结果的一致性;
所述步骤6包括以下步骤:基于空气污染物观测站点分布冗余度、城市内所有站点空气污染物空间覆盖比例、城市内所有站点空气污染物人口覆盖比例及城市内所有站点不同空气污染物空间代表性评估结果一致性,对城市内的观测站点布局进行评估,由于空间覆盖比例Pro_area、人口覆盖比例Pro_pop和不同空气污染物空间代表性评估结果一致性T均在
0‑1范围,且值越大,分布越合理;而冗余度值R越大,说明冗余程度越大,分布不合理,因此,提出如下综合指标:式(14)中,E为构建的空气污染物观测站点布局定量评价指标,该指标变化范围在0‑1之间,值越大说明站点布局越合理。