[发明专利]一种基于CWT的大气一次污染物外来源定量分析方法

权利要求书

1.一种基于CWT的大气一次污染物外来源定量分析方法，其特征在于：包括以下步骤：

步骤S1：采用气流轨迹HYSPLIT模型获取研究区后向轨迹数据；

步骤S2：将后向轨迹数据网格化，导入研究区大气污染物逐小时观测值，计算浓度权重轨迹CWT；

其中：C_ij是网格ij上的平均权重浓度；l是轨迹；M是轨迹总数；C_l是轨迹l经过网格ij时对应的所监测的大气污染物浓度值；T_ijl是轨迹l在网格ij停留的时间。

进一步引入权重因子W_ij，减小由于轨迹数据较小时引起的模型不确定性，得到WCWT，公式如下：

WCWT_ij＝W_ij×C_ij (2)

其中，C_ij是网格ij上的平均权重浓度；W_ij为权重函数；m、n分别为CWT网格的总行列数。

步骤S3：进一步计算浓度权重轨迹占比PCWT，高占比区域即为高潜在源区。公式如下：

2.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：步骤S1中气象资料来源为NCEP/NCAR提供的全球资料同化系统再分析资料。

3.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：步骤S1中气象资料选用分辨率0.5度*0.5度的网格数据，使用HYSPLIT软件进行后向轨迹计算，选择研究区的大气监测点坐标作为后向轨迹起始点进行模拟计算。

4.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：步骤S1中模拟时间设置为72h，每天分四个时段模拟，分别为：0点、6点、12点、18点，模拟高度在100-300米之间，将时间划分为春3～5月、夏6～8月、秋9～11月、冬12月～次年2月。

5.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：步骤S2中，网格分辨率选择小于0.5度。

6.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：为了减小网格上经过的所有轨迹数n_ij较小时的不确定性，引入了权重函数W_ij，公式如下：

其中：n_ij为该网格上经过的所有轨迹数，avg为研究区内每个网格内平均轨迹数。

7.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：当P_ij值较大时说明网格ij在总网格浓度权重轨迹WCWT值中占比大，对研究区域的影响较明显。

8.根据权利要求1所述的基于CWT的大气一次污染物外来源定量分析方法，其特征在于：通过行政边界和相对源区方向，统计出不同源区的浓度权重轨迹PCWT值，计算出对研究区影响较明显的区域。