06. 痕量气体立体反演与区域通量监测

本章核心目标：聚焦 DOAS 技术在大气环境监测前沿的工程实战应用；掌握典型痕量污染气体（${\text{NO}}_2,{\text{SO}}_2,\text{HCHO},\text{HONO},{\text{O}}_3$）反演窗口精细划分黄金准则；推导差分斜柱浓度（dSCD）向绝对垂直柱浓度（VCD）转换的几何近似与空气质量因子（AMF）辐射传输模型；构建车载移动走航围绕闭合边界估算工业烟羽质量排放通量（Emission Flux）的物理数学模型。

---

1. 典型大气痕量气体反演拟合窗口划分准则

光谱拟合窗口（Fitting Window）的选取是决定 DOAS 反演成败的首要关键。选取拟合窗口时必须严格权衡以下四大黄金工程准则：

1. 吸收截面强度准则：目标气体在该波段具有强烈的离散特征跃迁，差分峰-谷幅值 $\mathrm{\Delta }{\sigma }^{\prime }$ 足够大；
2. 光谱正交性准则：目标气体的指纹结构与窗口内其他共存干扰气体的谱线形状互相关系数极低，确保信息矩阵非奇异；
3. 避免强非线性饱和：避开光强急剧衰减的强饱和吸收区（如深紫外 $<305\text{ nm}$ 的极强臭氧 Hartley 吸收带，或近红外的大水汽饱和吸收带）；
4. 窗口跨度适中：拟合窗口宽度通常限定在 $15\sim 40\text{ nm}$。过窄导致采样自由度不足、易受随机噪声干扰；过宽则低阶多项式无法精确逼近高阶气溶胶消光曲率。

1.1 主流痕量气体拟合窗口标准参数表

根据国际高空与对流层大气成分探测网络（NDACC）与合肥物质科学研究院、德国海德堡大学推荐标准：

目标气体	推荐拟合波长窗口	典型共存拟合组分	科学意义与监测要点
二氧化氮 ${\text{NO}}_2$	$400-450\text{ nm}$ (或 $425-450\text{ nm}$)	${\text{NO}}_2(294\text{ K})$, ${\text{NO}}_2(220\text{ K})$, ${\text{O}}_3$, ${\text{O}}_4$, ${\text{H}}_2\text{O}$, $\text{Ring}$	采用常温/低温双截面解耦对流层与平流层贡献；化石燃烧与尾气核心指标。
二氧化硫 ${\text{SO}}_2$	$305-320\text{ nm}$ (或 $307.5-328\text{ nm}$)	${\text{SO}}_2$, ${\text{O}}_3(223\text{ K})$, ${\text{O}}_3(243\text{ K})$, ${\text{NO}}_2$, $\text{Ring}$	紫外深处信号较弱且受强臭氧遮挡；火电厂、钢铁厂脱硫效率评估与火山喷发预警。
甲醛 $\text{HCHO}$	$336.5-359\text{ nm}$	$\text{HCHO}$, ${\text{O}}_3$, ${\text{NO}}_2$, $\text{BrO}$, ${\text{O}}_4$, $\text{Ring}$	差分吸收极微弱（${\tau }^{\prime }<5\times {10}^{-4}$）；大气 VOCs 氧化活跃度与二次有机气溶胶前体物。
亚硝酸 $\text{HONO}$	$337-370\text{ nm}$ (或 $353-388\text{ nm}$)	$\text{HONO}$, ${\text{NO}}_2$, ${\text{O}}_3$, ${\text{O}}_4$, $\text{HCHO}$, $\text{Ring}$	日落夜间积累、清晨极速光解释放 $\text{OH}$ 自由基，驱动全天光化学烟雾爆发。
臭氧 ${\text{O}}_3$	$315-335\text{ nm}$ (Huggins 带) 或 $550-650\text{ nm}$ (Chappuis 带)	${\text{O}}_3(220\text{ K},273\text{ K})$, ${\text{NO}}_2$, ${\text{O}}_4$, ${\text{H}}_2\text{O}$, $\text{Ring}$	平流层臭氧保护伞与近地面有害臭氧污染双向反演。
氧二聚体 ${\text{O}}_4$	$360\text{ nm}$ / $477\text{ nm}$	${\text{O}}_4$ (${\text{O}}_2-{\text{O}}_2$ 碰撞诱导吸收)	大气的天然光学标尺：在大气中分布仅由气压平方决定，用于反演气溶胶光学厚度与光程修正。

📏 核心物理机理：氧二聚体 ${\text{O}}_4$ 的“天然光学标尺”角色

氧气（${\text{O}}_2$）在大气中的体积分数极其稳定（$20.95\mathrm{\%}$）。两个氧分子短暂碰撞诱导形成的碰撞复合物 ${\text{O}}_4$ 的数密度与空气压强的平方（$p^2$）成严格正比，其垂直理论剖面完全已知且常年不变。

• 若实测某仰角下 ${\text{O}}_4$ 吸收偏低，说明视线光子遇到了低空浓雾/重霾，光程被气溶胶强烈前向散射提前截断；
• 若 ${\text{O}}_4$ 吸收偏高，说明多重散射使光子在层内反复迂回折返；
• 因此，${\text{O}}_4$ 就像大自然在空气中放置的一把不可篡改的“绝对光程刻度尺”，是反演气溶胶消光廓线的黄金钥匙！

---

2. 从差分斜柱浓度 (dSCD) 到垂直柱浓度 (VCD) 的转换

反演算法直接解出的是沿实际曲折光路积分的差分斜柱浓度（Differential Slant Column Density, dSCD）：

$$\mathrm{\Delta }\text{SCD}(\alpha )=\text{SCD}(\alpha )-{\text{SCD}}_{\text{ref}}$$

然而，斜柱浓度受瞬时太阳位置、观测仰角 $\alpha$、气溶胶浑浊度剧烈影响，无法作为区域环境质量的客观横向评价标准。环境科学真正需要的是垂直柱浓度（Vertical Column Density, VCD）：

$$\text{VCD}\equiv {\int }_0^{\mathrm{\infty }}c(z)dz[\text{molecules}\cdot {\text{cm}}^{-2}]$$

2.1 空气质量因子（Air Mass Factor, AMF）的严密定义

引入无量纲转换系数——空气质量因子（AMF），它表征了视线斜程相较于纯垂直直穿光程的综合光学增益倍数：

核心转换公式：空气质量因子 AMF 与 VCD

$$\text{AMF}\equiv \frac{\text{SCD}}{\text{VCD}}$$

对于差分斜柱浓度 $\mathrm{\Delta }\text{SCD}(\alpha )=\text{SCD}(\alpha )-\text{SCD}({90}^{\circ })$，定义差分空气质量因子（Differential AMF, $\mathrm{\Delta }\text{AMF}$）：

$$\mathrm{\Delta }\text{AMF}(\alpha )\equiv \text{AMF}(\alpha )-\text{AMF}({90}^{\circ })$$

则对流层真实垂直柱浓度 ${\text{VCD}}_{\text{trop}}$ 的严格求解公式为：

$${\text{VCD}}_{\text{trop}}=\frac{\mathrm{\Delta }\text{SCD}(\alpha )}{\mathrm{\Delta }\text{AMF}(\alpha )}=\frac{\text{SCD}(\alpha )-\text{SCD}({90}^{\circ })}{\text{AMF}(\alpha )-\text{AMF}({90}^{\circ })}$$

2.2 几何近似空气质量因子（Geometric AMF）

在能见度较高、晴朗无云、气溶胶光学厚度较小（$\text{AOD}<0.3$）且观测仰角适中（$\alpha \ge {15}^{\circ }$）时，太阳散射光在对流层边界层内的穿行路径可近似视作直线几何：

$${\text{AMF}}_{\text{geo}}(\alpha )\approx \frac{1}{\sin \alpha }$$

天顶方向（${90}^{\circ }$）由于 $\sin ({90}^{\circ })=1$，其 ${\text{AMF}}_{\text{geo}}({90}^{\circ })\approx 1$。由此得到几何近似下的差分空气质量因子：

$$\mathrm{\Delta }{\text{AMF}}_{\text{geo}}(\alpha )=\frac{1}{\sin \alpha }-1$$

代入浓度转换公式：

$${\text{VCD}}_{\text{trop}}\approx \frac{\mathrm{\Delta }\text{SCD}(\alpha )}{\frac{1}{\sin \alpha }-1}$$

🎯 特殊工程几何：α = 30° 的奇妙对齐

当望远镜仰角取 $\alpha ={30}^{\circ }$ 时：

$$\sin ({30}^{\circ })=0.5⟹\mathrm{\Delta }{\text{AMF}}_{\text{geo}}({30}^{\circ })=\frac{1}{0.5}-1=2-1=1.0$$

此时测得的差分斜柱浓度在数值上直接近似等于对流层垂直柱浓度！

$${\text{VCD}}_{\text{trop}}\approx \mathrm{\Delta }\text{SCD}({30}^{\circ })$$

这正是几乎所有 MAX-DOAS 走航车与固定监测站都将 ${30}^{\circ }$ 仰角作为快速直读观测核心角度的深层物理原因。

2.3 辐射传输模型（RTM）下的高精度 AMF 积分

在严重灰霾或近水平低仰角（$1^{\circ }\sim 5^{\circ }$）条件下，强气溶胶多重散射导致直线几何近似失效。此时必须利用严密的**大气辐射传输模型（如 SCIATRAN, VLIDORT, DISORT）**求解微观光子传输过程。

定义微观高度层 $z$ 处的层空气质量因子（Box-AMF） $w(z,\alpha ,\lambda )$：

$$w(z,\alpha ,\lambda )\equiv \frac{\partial \text{SCD}(\alpha )}{\partial \text{VCD}(z)}$$

真实全层空气质量因子由目标气体的归一化垂直分布廓线 $s(z)=c(z)/\text{VCD}$ 与 Box-AMF 积分加权求得：

$$\text{AMF}(\alpha ,\lambda )={\int }_0^{\mathrm{\infty }}w(z,\alpha ,\lambda )\cdot \frac{c(z)}{\text{VCD}}dz$$

Box-AMF 综合考虑了地表反照率（Albedo）、气溶胶光学厚度（AOD）、气溶胶单次散射反照率（SSA）、不对称因子（$g$）以及雷利散射全相位矩阵。

---

3. 基于最优估计方法 (OEM) 的高分辨率垂直廓线反演

仅仅获取垂直柱总量（VCD）往往不足以指导精细化环保施策。污染究竟是滞留在近地面 $100\text{ m}$ 内毒害人体，还是高空 $2\text{ km}$ 处的远距离外来输送气团？

MAX-DOAS 多仰角扫描序列（$1^{\circ },2^{\circ },3^{\circ },5^{\circ },{15}^{\circ },{30}^{\circ }$）通过**最优估计方法（Optimal Estimation Method, OEM）**能够高分辨率反演出气溶胶与痕量气体的垂直剖面。

根据 Clive Rodgers 的现代地球物理反演论，状态向量 $\mathbf{x}$（垂直剖面）满足代价函数极小化：

Rodgers 最优估计反演代价函数

$${\chi }^2(\mathbf{x})={[\mathbf{y}-\mathbf{F}(\mathbf{x})]}^T{\mathbf{S}}_{\epsilon }^{-1}[\mathbf{y}-\mathbf{F}(\mathbf{x})]+{[\mathbf{x}-{\mathbf{x}}_a]}^T{\mathbf{S}}_a^{-1}[\mathbf{x}-{\mathbf{x}}_a]$$

• $\mathbf{y}$ 为各扫描仰角实测的 $\mathrm{\Delta }\text{SCD}$ 向量；
• $\mathbf{F}(\mathbf{x})$ 为由大气辐射传输模型驱动的正向非线性观测算子；
• ${\mathbf{x}}_a$ 与 ${\mathbf{S}}_a$ 分别为先验廓线及其先验误差协方差矩阵；
• ${\mathbf{S}}_{\epsilon }$ 为光谱拟合残差协方差矩阵。

该算法能够在无探空气球或激光雷达的情况下，凭借被动日光散射反演出地面至 $3\text{ km}$ 高度、垂直层厚 $100\text{ m}$ 的污染立体梯队分布。

---

4. 车载移动走航与工业烟羽排放通量 (Emission Flux) 估算

在工业园区偷排核查与突发泄漏溯源中，固定式测站存在监控盲区。将 MAX-DOAS 仪器与高精度 GPS、气象观测仪集成于机动监测车上，构成移动车载 MAX-DOAS（Mobile MAX-DOAS）。

4.1 高斯散度定理与闭合环路质量通量积分

将监测车以恒定车速沿大型工业园区、化工基地或孤立火电厂的外围公路进行围绕闭合环路 $\mathrm{\Gamma }$ 行驶一周。

🎣 生动比喻：大气的无形捕鱼网与高斯通量算子

走航车绕行工厂一周，就像在工厂四周的大气中拉起了一面高耸入云的**“无形光学捕鱼网”**：

• 车顶的望远镜时刻测量垂直穿透该网的光学柱浓度 $\text{VCD}$；
• 气象风速仪测定穿透该网的风矢量大小与角度；
• 根据高等数学中的高斯散度定理（Gauss Divergence Theorem），从网中吹出来的气体总量减去吹进去的自然本底，就是网内源头真实的排放通量！环保执法人员甚至不需要踏进工厂大门半步，就能把排放速率算得清清楚楚。

根据质量守恒与高斯散度定理，当排放源在绕行周期内处于准稳态、且边界内无其他未知汇过程时，从该闭合多边形区域向外净逸散的**气体质量排放通量 $\mathrm{\Phi }$（Emission Flux）**严格等于垂直柱浓度与水平风场沿闭合边界的线积分：

闭合边界工业质量排放通量公式

$$\mathrm{\Phi }={\oint }_{\mathrm{\Gamma }}{\text{VCD}}_{\text{trop}}(s)\cdot ({\mathbf{v}}_{\text{wind}}(z_{\text{plume}})\cdot \mathbf{n}(s))ds$$

• $s$ 为沿监测车行车轨迹的弧长微分；
• ${\text{VCD}}_{\text{trop}}(s)$ 为位置 $s$ 处反演出的对流层垂直柱浓度；
• ${\mathbf{v}}_{\text{wind}}(z_{\text{plume}})$ 为有效烟羽抬升高度 $z_{\text{plume}}$ 处的平均风矢量（包含风速 $u$ 与风向角 ${\theta }_{\text{wind}}$）；
• $\mathbf{n}(s)$ 为行车路径向外的局部二维法向量；
• ${\mathbf{v}}_{\text{wind}}\cdot \mathbf{n}=u\cdot \sin ({\theta }_{\text{drive}}-{\theta }_{\text{wind}})$ 为垂直于行进路线的风速投影分量。

4.2 离散化工程数值算法

在实地采集的一组包含 $N_p$ 个连续 GPS 轨迹点上：

$$\mathrm{\Phi }=\sum _{k=1}^{N_p}({\text{VCD}}_k-{\text{VCD}}_{\text{bkg}})\cdot u_k\cdot \sin ({\theta }_k-{\theta }_{\text{wind}})\cdot \mathrm{\Delta }{L}_k\cdot \frac{M_{\text{mol}}}{N_A}$$

其中：

• ${\text{VCD}}_{\text{bkg}}$：上风向洁净对照段的大气背景基底柱浓度，作差扣除外部区域输入本底；
• $\mathrm{\Delta }{L}_k$：两连续测量点之间的行车几何直线距离（米）；
• $M_{\text{mol}}$：待测气体摩尔质量（例如 ${\text{NO}}_2$ 为 $46.0055\text{ g/mol}$，${\text{SO}}_2$ 为 $64.066\text{ g/mol}$）；
• $N_A=6.02214\times {10}^{23}\text{ molec/mol}$ 为阿伏伽德罗常数。

最终输出标准工业排放速率：$\text{kg/h}$（千克/小时） 或 $\text{t/a}$（吨/年）。

工业实战案例计算

• 某大型煤电厂四周公路周长 $12\text{ km}$，走航车以 $30\text{ km/h}$ 匀速环测用时 $24$ 分钟；

• 烟羽下风向 $2\text{ km}$ 宽度的羽流区内，${\text{NO}}_2$ 垂直柱浓度出现巨型单峰，平均增量 $\mathrm{\Delta }\text{VCD}=3.5\times {10}^{16}{\text{ molec/cm}}^2$；

• 测风激光雷达探测到烟羽中心抬升高度为 $450\text{ m}$，该高度处风速 $u=6.0\text{ m/s}$，垂直于行车路线；

• 计算逸散通量：

  $${\mathrm{\Phi }}_{{\text{NO}}_2}=(3.5\times {10}^{16}{\text{ molec/cm}}^2)\times {10}^4{\text{ cm}}^2/{\text{m}}^2\times (6.0\text{ m/s})\times 2000\text{ m}=4.2\times {10}^{24}\text{ molec/s}$$$${\mathrm{\Phi }}_{\text{mass}}=\frac{4.2\times {10}^{24}\text{ molec/s}}{6.022\times {10}^{23}\text{ molec/mol}}\times 46.0\text{ g/mol}\times {10}^{-3}\text{ kg/g}\approx 0.3208\text{ kg/s}\approx 1155\text{ kg/h}!$$

环保执法机构无需进入工厂大门，仅凭外围公路上行驶一圈，即可在半小时内以无可辩驳的物理光学实测数据，锁定该厂每小时超标排放 $1.15$ 吨 ${\text{NO}}_2$ 的违法事实！

---

5. 专栏全景总结

至此，我们完成了【DOAS 差分吸收光谱技术】从基础微观光学到前沿宏观监测的全链条推导与构建：

DOAS 展现了理论物理、统计应用数学与现代光电工程最精妙的结合：它不回避自然界复杂介质的浑浊与不完美，而是通过深刻把握物理机制的内在频率尺度差异，用数学差分化腐朽为神奇，在纷繁杂乱的光学背景中，为人类守护起清晰洞察大气环境的“火眼金睛”。