实战：用WRF-Chem V3.9.1.1模拟一次华北雾霾过程（附完整namelist配置与排放数据处理心得）
实战：用WRF-Chem V3.9.1.1模拟华北雾霾的完整技术指南
华北地区秋冬季雾霾问题一直是环境科学研究的重点。本文将基于WRF-Chem V3.9.1.1版本，详细介绍如何从零搭建一个针对华北雾霾事件的数值模拟系统。不同于基础教程，我们聚焦于实际科研项目中遇到的配置难点和数据处理技巧，特别是化学机理选择、排放源处理以及关键参数优化等核心环节。
1. 模型基础配置与环境准备
1.1 系统依赖与编译选项
WRF-Chem的编译需要特别注意化学模块的兼容性。推荐使用Intel编译器搭配MKL数学库，以下是在Linux系统下的典型环境变量设置：
export NETCDF=/path/to/netcdf export WRF_CHEM=1 # 必须开启化学模块 export WRF_KPP=1 # 如需使用KPP化学机制编译时常见的几个关键选项：
- 32位/64位选择：大型模拟建议使用-dmpar选项进行并行编译
- I/O性能优化：添加-DNETCDF4选项支持NetCDF4格式输出
- 化学机制选择：根据后续要使用的chem_opt参数确定需要编译的化学模块
1.2 计算域设置技巧
华北地区雾霾模拟的domain设计需要兼顾：
- 水平分辨率：外层域建议15-30km，内层域3-5km
- 垂直层数：至少40层，边界层内需加密（近地面10层占50%高度）
- 地形数据：使用USGS 30秒或更高精度数据
典型namelist.wps配置示例：
&geogrid e_we = 120, 151, 211 e_sn = 120, 151, 211 geog_data_res = '10m','2m','30s' dx = 27000, 9000, 3000 dy = 27000, 9000, 3000 /提示：雾霾模拟建议包含整个华北平原及周边山地，以准确捕捉污染物区域输送过程。
2. 化学机制与气溶胶方案选择
2.1 化学机理对比分析
WRF-Chem提供多种化学机制，针对华北雾霾推荐组合方案：
| 化学机制 (chem_opt) | 气溶胶方案 | 适用场景 | 计算开销 |
|---|---|---|---|
| 7 (CBMZ) | MOSAIC 4bin | 常规雾霾 | 中等 |
| 8 (CBMZ) | MOSAIC 8bin | 精细粒径分析 | 较高 |
| 11 (RADM2+MADE) | SORGAM | 硫酸盐主导 | 较低 |
华北雾霾中二次有机气溶胶(SOA)和硝酸盐贡献显著，建议：
chem_opt = 8 # CBMZ+MOSAIC 8bin aer_op_opt = 3 # 精确体积近似光学特性 gas_drydep_opt = 1 # 启用气体干沉降 aer_drydep_opt = 1 # 启用气溶胶干沉降2.2 时间步长协调策略
化学过程的时间步长设置需要与动力过程协调：
chemdt = 5.0 # 化学步长(分钟) photdt = 30.0 # 光解频率 bioemdt = 60.0 # 生物排放更新间隔注意：chemdt通常设置为动力步长的3-6倍，过大会导致数值不稳定，过小则显著增加计算量。
3. 排放源处理实战技巧
3.1 人为排放数据准备
华北地区推荐使用MEIC排放清单，处理流程包括：
- 时间分解：将年排放量按月份、星期、小时模式分解
- 物种映射：匹配排放清单物种与化学机制要求
- 垂直分配：根据行业类型确定排放高度
典型配置参数：
emiss_opt = 4 # CBMZ/MOSAIC排放 emiss_inpt_opt = 102 # RADM2转CBMZ io_style_emissions = 2 # 特定日期排放 auxinput5_interval = 60 # 每小时排放输入3.2 生物源与自然源处理
生物排放配置要点：
- MEGAN模型：适用于植被VOC排放
- GOCART方案：处理沙尘和海盐
bio_emiss_opt = 3 # 在线MEGAN dust_opt = 1 # GOCART沙尘 seas_opt = 1 # 海盐排放 ne_area = 150 # 化学物种数排放预处理工具推荐：
- PREP-CHEM-SRC：处理全球排放数据
- SMOKE：适用于区域高分辨率排放
4. 关键参数优化与调试
4.1 化学-气象耦合设置
气溶胶-辐射反馈是雾霾模拟的关键：
aer_ra_feedback = 1 # 启用辐射反馈 ra_lw_physics = 4 # RRTMG长波 ra_sw_physics = 4 # RRTMG短波4.2 常见问题排查指南
典型错误及解决方法：
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 浓度负值 | 化学步长过大 | 减小chemdt |
| 模式崩溃 | 排放数据缺失 | 检查wrfchemi文件 |
| 异常高值 | 干沉降未启用 | 设置drydep_opt=1 |
调试建议流程：
- 先关闭化学模块运行纯气象模拟
- 逐步开启化学过程（gas→aero→cloud）
- 使用chemdiag=1输出诊断信息
4.3 后处理与可视化
推荐分析工具链：
- NCL/Python：处理NetCDF输出
- VERDI：专业环境数据可视化
- GIS：空间分析与制图
关键诊断变量：
# Python示例：提取PM2.5浓度 import netCDF4 as nc ds = nc.Dataset('wrfout_d01') pm25 = ds['PM2_5_DRY'][:] # 干态PM2.55. 实际案例：2016年华北雾霾事件模拟
5.1 案例配置详情
以2016年12月重污染过程为例：
- 时间范围：2016-12-16至2016-12-21
- 初始场：NCEP FNL 1°×1°数据
- 排放清单：MEIC 2016基准年数据
完整namelist.chem配置片段：
&chem chem_opt = 8 chemdt = 5.0 bioemdt = 60.0 photdt = 30.0 emiss_opt = 4 emiss_inpt_opt = 102 io_style_emissions = 2 /5.2 结果验证方法
建议采用多种观测数据验证：
- 地面站数据：中国环境监测总站小时值
- 卫星反演：MODIS AOD产品
- 激光雷达：气溶胶垂直分布
验证指标计算示例：
# 计算相关系数 from scipy import stats r, p = stats.pearsonr(obs, model) print(f"Correlation: {r:.3f}")在多次实际项目验证中，这套配置对PM2.5模拟的相关系数通常能达到0.6-0.8，能较好再现污染过程的时空演变特征。特别是MOSAIC 8bin方案对二次组分比例的模拟效果明显优于4bin版本。