DOAS 差分吸收光谱是怎么算出来的：10 张图

01

吸收与消光：Beer-Lambert 定律的微元积分链

光强不是线性衰减，而是以指数方式急剧被吞噬。微元光程上的吸收累加，构成了所有光学遥测的最底座基石。

Fig 01I(λ) = I₀(λ) · exp[ - ∫ ε(s, λ) ds ] ➔ 负指数衰减

读图光子沿着视线路径每前进一步，就会被分子吸收或粒子散射，呈严格指数衰减。

关键测量量永远是光强比值 $ln(I_0/I)$，而不是绝对光强瓦特数。

挑战在真实开放大气中，光程 $L$ 与入射光强 $I_0$ 是无法放在比色皿里单独标定的！

02

消光三要素：气体吸收、分子瑞利散射与气溶胶米氏消光

总消光不是单一成分。它是痕量气体“指纹吸收”与空气分子瑞利散射、雾霾颗粒米氏散射的严格叠加。

Fig 02ε(λ) = Σ σ_j·c_j + ε_Rayleigh(λ⁻⁴) + ε_Mie(λ⁻α)

光谱动态解耦看板 (波长: 400 nm)

■ 瑞利散射 (Rayleigh): 58.4% (依赖 λ⁻⁴, 极陡)
■ 米氏消光 (Mie): 31.2% (依赖 λ⁻¹, 雾霾宽带)
■ NO₂ 指纹吸收: 特征峰顶 (微弱快速跃迁)

■ 核心物理洞察：频域正交解耦

1. 散射背景（瑞利+米氏）全都是低频慢变函数（随波长平缓漂移）。
2. 目标污染气体全都是高频快速振荡波纹（分子振动能级指纹）。
3. 两者在频域正交——用一个低阶多项式即可彻底剥离全部散射背景！

■ SCANNING: 波长扫描进行中: 350nm ~ 550nm 连续光谱解耦

Rayleigh: 58% | Mie: 31% | NO2: 峰顶

物理大雾天光强暴跌 90%，是因为米氏气溶胶消光，而不是污染气体暴增。

指纹每种气体在紫外可见区有独一无二的电子跃迁谱带，宛如宇宙身份证条形码。

结论只要把低频平滑背景滤掉，就能在浓雾霾天里精准测出几 ppb 的极微弱气体！

03

频域手术刀：快慢变分离的数学解耦

把吸收截面拆为两部分：$sigma(lambda) = sigma_0(lambda) + sigma'(lambda)$。慢变放进多项式背景，快变留下拟合气体浓度。

Fig 03ln(I₀'/I) = Σ σ'_j · SCD_j + Residual(λ)

数学绝技用一个低阶多项式 $P(lambda)$ 吸收所有不可测的宽带干扰。

比喻就像降噪耳机：把低频轰鸣环境音直接相消，只把歌手清脆的高频唇齿音保留下来。

边界多项式阶数不能太高（通常 3~5 阶）。阶数过高会把气体的吸收指纹当成背景抹杀掉！

04

非线性反演：Levenberg-Marquardt 自适应阻尼收敛

波长漂移（Shift）是非线性的。必须在最速下降法与高斯-牛顿法之间动态平衡，快速求解斜柱浓度 SCD。

Fig 04(J^T·J + μ·I) Δp = J^T·r ➔ 阻尼最小二乘迭代

LM 阻尼自适应流转状态机

拟合残差与收敛指标监控台

4.2 × 10⁻² RMS 拟合残差

当前迭代轮次: 1 / 4
阻尼因子 μ: 10.0 (最速下降探索)
NO₂ 斜柱浓度 (SCD): 5.10 ± 1.20 × 10¹⁶
波长漂移 (Shift): +0.035 nm (动态对齐中)
拟合优度 (Chi-square): 2.840

■ ITERATION: 初始粗估计，阻尼因子较大，以梯度下降法探索

RMS: 4.2e-2 | μ: 10.0

阻尼因子 μ$mu$ 极大时像盲人摸象（稳健但慢），$mu$ 极小时像高斯跳跃（极快但容易发散）。LM 自动在此二者间动态游刃有余。

非线性项主要是由于光谱仪热胀冷缩引起的亚像元波长微小漂移。

判据如果最终拟合残差谱里还能看出一道道明显的波纹，说明漏选了气体吸收截面！

05

Czerny-Turner 光路：色散光栅与微米级狭缝的精密装配

光纤将光线导入 50µm 狭缝，准直镜打平射向衍射光栅，聚焦镜将彩虹光束投射到 2048 像元线阵 CCD。光机电的刚性决定了光谱测量的下限。

Fig 05入射狭缝 ➔ 准直镜 M1 ➔ 全息光栅色散 ➔ 聚焦镜 M2 ➔ 线阵 CCD

色散原理光栅表面每毫米刻有 1200 道精密凹槽，利用多光束干涉将复合波长在空间角度上拉开。

暗电流CCD 必须在深低温（$-50^circ ext{C}$）下运行，把硅晶片热激发的暗电子噪声压制到几乎为零。

热稳定性环境温度变动 1°C，光谱在探测器上就会漂移半个像元！因此密闭机箱必须精密恒温。

06

宇宙条形码：太阳夫琅禾费吸收线与亚像元自校准

太阳光球层的原子蒸气留下了数千道极其锐利的吸收暗线。它不是干扰，而是大自然赐予光谱仪永恒绝对的在线波长校准刻度尺。

Fig 06Fraunhofer Solar Atlas 卷积与亚像元 (< 0.001nm) 锁定

来源太阳内部极热发出的连续白光，穿过太阳自身较冷的外部大气时，被铁、氢、钙等原子吸收留下的刻痕。

自校准无论地面仪器怎么震动温漂，只要对齐夫琅禾费这几根深坑，波长准确度立刻校准回绝对基准。

必须项高精度的 DOAS 拟合绝对不能用出厂标定的波长死数，每次反演都必须动态重对齐！

07

环效应：非弹性旋转拉曼散射的深谷填平

空气中氮氧分子发生非弹性旋转拉曼跃迁，改变了散射光子的波长。连续谱光子跳进夫琅禾费深谷，造成吸收线“被神秘填浅”2%~10%。

Fig 07Rotational Raman Scattering (RRS) 填平伪信号补偿与 Ring 谱

■ DOAS 破解法：专用的 Ring 谱

1. 量子机理：N₂ 与 O₂ 转动能级跃迁（$Delta J = 0, pm 2$），散射光子被赋予微小频移。

2. 致命后果：填平效应伪装成未知气体的“假吸收”，导致 NO₂ 产生严重拟合偏差！

3. 破局解法：理论计算出纯拉曼/瑞利比值谱（Ring 谱），作为**独立拟合通道**加入矩阵相消！

■ RING EFFECT: 非弹性拉曼频移光子向深谷滑落，造成 8% 虚假吸收信号

谷底抬升: +8.2% | Ring 谱补偿状态: 注入矩阵

现象命名1962 年由 Grainger 和 Ring 首次在观测夜空月光与白天蓝天对比时发现，故名 Ring Effect。

量级填平深度达 2%~10%，远大于 NO₂ 本身 0.1% 的弱吸收！不扣除 Ring 谱，反演结果全错。

物理本质它不是仪器问题，而是大自然量子非弹性散射带来的物理光学伴生物。

08

MAX-DOAS 几何放大：低仰角穿透使吸收光程暴增 57 倍

天顶（90°）只能竖直看穿薄薄的污染层；将望远镜倾斜到 1°~3° 贴地平扫，光线在污染边界层内行程被几何倍数拉长。探测灵敏度直接提升近两个数量级。

Fig 08有效光程 L = H / sin(α) ➔ 几何放大因子 AMF_geom

几何公式在无气溶胶纯净视线下，几何空气质量因子简单为 $ ext{AMF}_{geom} = 1 / sin(alpha)$。

差分消背景用同一次循环扫描的 $90^circ$ 天顶谱作为参考谱 $I_{ref}$，直接消除了平流层高空臭氧与二氧化氮的全部背景干扰！

散射限制在重霾天中，由于强烈的米氏多次散射，光线跑不到 57km 就被散射折断，此时需动用辐射传输模型（RTM）求解实际 AMF。

09

氧气二聚体 O₄：大自然的天然光学标尺

两个氧分子碰撞产生的瞬态四聚体 $ ext{O}_4$。它的真实浓度仅取决于大气压力平方，是雷打不动的已知基准。

Fig 09O₂-O₂ 碰撞吸收 ➔ 测定实际光程与气溶胶消光廓线

为什么需要 O₄ 标尺？

1. 雾霾中光子到底在空中被折返了几公里，人类无法直接看到。
2. 但 $ ext{O}_4$ 的浓度剖面是恒定的（严格正比于气压平方 $[O_2]^2$）。
3. 测出的 $ ext{O}_4$ 偏低 ➔ 说明气溶胶多，光线提前弹飞，实际光程变短！

■ 两步反演流水线

Step 1: 反演 $ ext{O}_4$ ➔ 求解实际气溶胶光学厚度 AOD 与消光垂直廓线。
Step 2: 将气溶胶场输入 RTM 模型，精确算出真实空气质量因子 AMF，求得 NO₂ 的真浓度。

■ O4 OPTICAL RULER: 晴空状态：无气溶胶压制，光子深度穿透，O₄ 吸收达 100% 理论基准

O4 吸收强度: 100% | 气溶胶 AOD: 0.12 (晴空)

碰撞吸收单个氧分子在可见光不吸收，但在分子碰撞瞬间，外层电子云重叠跃迁产生吸收，强度正比于 $[O_2]^2$。

不可替代没有 $ ext{O}_4$，人类就无法测量雾霾对光程的缩短效应，反演误差高达 100% 以上。

校正因子国际学术界仍在使用 1.05~1.2 的 O4 截面缩放系数，这是当前大气光学前沿课题之一。

10

闭环实战：从柱浓度 dSCD 到工业园区排放通量（吨/小时）

车载走航围绕工业园区闭合巡游一圈。将光学垂直柱浓度与水平风场做高斯通量线积分，立刻算出污染源每秒钟往天上喷多少气体。

Fig 10Flux = ∮ VCD(s) · v_wind · sin(θ) ds · M ➔ 高斯质量守恒通量

高斯定理只要沿闭合曲线跑一圈，内部所有未知名源的净排放总量就严格等于通过边界流出的通量差。

秒杀原位传统地面固定监测站只能测车屁股边几个 ppb，而 DOAS 车载走航直接把整座城市的立体羽流彻底锁死。

致命因子风速与风向误差直接以线性比例放大通量误差！因此必须配合微波风廓线雷达精准测风。