Page 239 - 《高原气象》2026年第1期
P. 239

1 期                      武泽昊等:2016年秋季海口市臭氧污染来源解析模拟研究                                       235

































                                              图1 WRF-Chem模拟区域(填色, 单位: m)
                                       黑色点表示国家基准气候站, 黑色三角表示国控环境空气质量监测站
                  Fig. 1 WRF-Chem simulation area (color represents elevation, unit: m). The black dots indicate the national benchmark
                        climate stations, and the black triangles indicate the state-controlled ambient air quality monitoring stations
                   表1  WRF-Chem模式的物理和化学参数化方案                     制, 模拟期间累计释放 2183个示踪粒子, 单粒子生
                  Table 1  Physicochemical parameterization scheme   命周期设定为 48 h。该配置有效平衡了时空分辨率
                              of WRF-Chem mode
                                                                 与计算效率, 确保对污染物输送路径的动态捕捉
                参数      物理/化学                                    能力。
                                           参数化方案
                类型      模块名称
                                                                 2. 2 过程分析方法
                物理     微物理方案          Purdue Lin(Lin et al, 1983)    过程分析方法将影响 O 的物理化学过程单独
                                                                                           3
                参数     边界层方案           MYJ(Hong et al, 2006)     输出, 通过对单独变量的分析研究O 浓度变化的成
                化方                                                                                3
                     积云参数化方案         GF(Grell and Dévényi, 2002)  因。模式中的污染物浓度的变化主要由以下几个
                 案
                        长波辐射          Goddard(Iacono et al, 2008)  部分构成:
                        短波辐射          Goddard(Iacono et al, 2008)                                          (1)
                                                                            ∆C = ∆C chem + ∆C phy + ∆C emis
                        陆面过程      Noah(Chen et al, 2001; Ek et al, 2003)  式中: ∆C chem 、 ∆C phy 和 ∆C emis 分别代表化学过程、 物
                化学    气相化学方案         RADM2(Stockwell et al, 1990)  理过程和排放变化引起的污染物浓度的变化。物
                参数     气溶胶方案         SORGAM(Schell et al, 2001)  理过程包括平流、 垂直混合和干湿沉降, 化学过程
                化方
                        光解方案       F-TUV(Madronich and Flocke, 1999)  主要为气相化学, 通过将理化过程影响的O 浓度变
                 案                                                                                      3
                       气溶胶反馈                  开启                 化单独输出, 分析O 浓度变化的原因。本文中理化
                                                                                   3
                                                                 过 程 变 量 包 括 : (1)光 化 学 反 应 过 程 净 贡 献
               模拟, 基于后向轨迹, 采用潜在源贡献因子算法                          (Chem_O), 该过程与温度、 辐射等气象因子有
                                                                          3
              (Potential source analysis function, PSCF)和权重轨     关; (2)垂直混合与干沉降过程净贡献(Vmix_O),
                                                                                                            3
               迹分析法(Concentration weight trajectory, CWT)分       该过程与边界层高度和相对湿度等气象因子有关;
               析 海 口 O 污 染 的 潜 在 源 区(Seibert  et  al,  2004;    (3)平 流 输 送 过 程 净 贡 献(ADVT_O = Advh_O +
                        3
                                                                                                             3
                                                                                                    3
               Wang et al, 2009)。WRF-Chem 模式模拟区域 3 的
                                                                 Advz_O), 该过程受风场的影响。式(1)中的 ∆C emis
                                                                        3
               气象场作为 FLEXPART 模式气象驱动场。FLEX‐                      是 dt 时刻的 O 直接排放量, 而 O 作为大气中典型
                                                                                               3
                                                                              3
               PART模式参数设定如下: 以海口主城区(19. 96°N,                    的二次污染物, 实际并不存在直接排放, 因此这里
               110. 27°E)为粒子释放点, 采用小时级连续释放机                      的∆C emis 值为0。
   234   235   236   237   238   239   240   241   242   243   244