05. 实验标定与数据处理链：校准、平场与环效应

本章核心目标：构建从探测器原始数字电荷量（Raw ADU Counts）到标准科学吸收光谱的严密预处理链；详解电子偏置（Offset）与热激发暗电流（Dark Current）的消除实验；掌握基于太阳夫琅禾费光谱（Fraunhofer Atlas）“天体物理指纹对齐”的亚像元级在线波长对齐算法；深入推导大气非弹性旋转拉曼散射（Rotational Raman Scattering）引起“环效应（Ring Effect）”的量子物理微观机制与校准谱生成流程。

---

1. 探测器原始信号标定与科学预处理链

光谱仪 CCD 探测器输出的原始数据是以模数转换器（ADC）计数值（Analog-to-Digital Units, ADU Counts）存储的无量纲整数。在输入反演方程前，必须对其进行严密的光电信号重构。

1.1 电子偏置（Electronic Offset）扣除

• 物理成因：为了避免探测器模拟放大电路在微弱读出噪声波动时电压跌入负值被 ADC 强制截断为 0 造成统计失真，硬件电路在基准电平上人为施加了一个恒定的正向直流偏置电压（通常对应 $500\sim 1000\text{ ADU}$）；
• 测量规范：

  1. 关闭系统光学快门（或彻底封盖狭缝视窗）；

  2. 将积分曝光时间设置为硬件支持的理论绝对最小值（$t_{\text{exp}}\to 0$，通常为数微秒）；

  3. 连续采集 $M\ge 50$ 帧零曝光数据并进行时间轴平均，以将叠加在偏置上的单次读出随机噪声压低 $\sqrt{M}$ 倍：

     $${\overline{S}}_{\text{offset}}(i)=\frac{1}{M}\sum _{k=1}^M{S}_{\text{offset},k}(i)$$

1.2 芯片热激发暗电流（Dark Current）标定

• 物理成因：即便在无任何光子入射的绝对暗室中，硅半导体价带中的电子由于晶格热振动也会以恒定统计几率被激发跃迁至导带，并在势阱中积聚，形成暗电流；

• 数学规律：在恒温制冷控制系统维持温度 $T$ 稳定的前提下，暗电荷累积量与曝光时间 $t_{\text{exp}}$ 呈现严格的线性正比关系：

  $$S_{\text{dark}}(i,T,t_{\text{exp}})=D_0(i,T)\cdot t_{\text{exp}}$$

  其中 $D_0(i,T)$ 为像元 $i$ 在温度 $T$ 下的固有暗电流率（单位：$\text{ADU}\cdot {\text{s}}^{-1}$）；

• 动态扣除策略： 野外 MAX-DOAS 系统通常随光照强度自适应调整曝光时间（如晴天正午 $0.2\text{ s}$，傍晚阴天 $5\text{ s}$）。系统在日落后或测量间隙自动关闭光阑，针对不同曝光时间阶梯测绘暗电流曲线，或直接在每次科学测量后紧接着采集一幅同等时长同等温度的暗场帧进行逐点相减：

  $$S_{\text{clean}}(i)=S_{\text{raw}}(i)-{\overline{S}}_{\text{offset}}(i)-S_{\text{dark}}(i,T,t_{\text{exp}})$$

---

2. 波长精确校准与仪器线型函数 (Slit Function)

在 DOAS 反演中，波长轴精度的容忍度极其严苛。微小至 $0.005\text{ nm}$ 的波长漂移，在对数差分 $\ln (I/I_{\text{ref}})$ 中即可诱发出高达 ${10}^{-2}$ 量级的强烈伪差分吸收峰，比大多数弱吸收气体信号大出两个数量级以上！

2.1 实验室标准谱线灯离散标定

在实验室受控环境中，使用低压汞（Hg）灯或汞氩（Hg-Ar）复合空心阴极灯照射积分球，获取一系列自然界基准原子发射谱线。

发射元素	理论真空波长	对应空气特征波长
Hg	$313.155\text{ nm}$	$313.064\text{ nm}$
Hg	$365.484\text{ nm}$	$365.381\text{ nm}$
Hg	$404.656\text{ nm}$	$404.542\text{ nm}$
Hg	$435.833\text{ nm}$	$435.711\text{ nm}$
Hg	$546.075\text{ nm}$	$545.922\text{ nm}$

通过高斯多峰拟合精确求得各标准谱线在探测器上的像元峰位中心 $p_k$，然后拟合像元-波长多项式色散方程：

$$\lambda (p)=c_0+c_1\cdot p+c_2\cdot p^2+c_3\cdot p^3$$

同时，原子发射谱线的固有量子自然线宽（$<0.001\text{ nm}$）远小于光谱仪狭缝宽度，因此实测发射线的轮廓即为仪器的狭缝响应函数（Slit Function）$W(\lambda )$。

2.2 基于太阳夫琅禾费谱线 (Fraunhofer Atlas) 的亚像元在线对齐

在长期野外监测中，仪器不可避免地受到环境气温热胀冷缩与机械蠕变影响。实验室预校准的波长曲线会随时间漂移。

🌌 形象比喻：天体物理指纹对齐 / 宇宙条形码锁频

太阳自身就是一个永恒不变且充满丰富吸收特征的天然校准光源！在太阳约 $5778\text{ K}$ 的光球层中，铁（$\text{Fe}$）、钙（$\text{Ca}$）、镁（$\text{Mg}$）、钠（$\text{Na}$）和氢（$\text{H}$）等原子的共振吸收，在连续黑体辐射谱上刻下了成千上万道极其深邃锐利的暗线——夫琅禾费吸收线（Fraunhofer Lines）。 这些谱线宛如大自然刻在阳光里的**“宇宙条形码”**，其物理波长由量子力学第一性原理决定，万年不变。我们只需将实测光谱与超高分辨率太阳物理谱库（如 Kurucz Atlas）进行亚像元卷积对齐，就能随时随地把波长轴锁死在 $0.001\text{ nm}$ 的绝对精度上！

---

3. 非弹性旋转拉曼散射与“环效应 (Ring Effect)”

在被动散射日光 DOAS 观测的早期历史上，科学家经常遭遇一个令人费解的怪象：即便扣除了所有已知污染气体的吸收谱，天顶散射光与直接太阳光的拟合残差中依然残留着幅度高达数个百分点的巨大特征振荡，且振荡峰位与太阳夫琅禾费吸收线的位置严丝合缝地重合。

这一物理现象由 J. F. Grainger 和 J. Ring 于 1962 年首次在文献中报道，历史上被命名为环效应（Ring Effect）。

3.1 环效应的微观量子物理机理

• 误区排除：环效应不是某种未知化学气体的吸收，而是日光光子在大气传输中经历的微观非弹性散射（Inelastic Scattering）；

• 旋转拉曼跃迁（Rotational Raman Scattering, RRS）： 光子与大气主成分对称双原子分子（${\text{N}}_2$ 和 ${\text{O}}_2$）碰撞时，约 $96\mathrm{\%}$ 发生保频的弹性瑞利散射（$\mathrm{\Delta }J=0$）；但约有 $4\mathrm{\%}$ 的光子会诱导分子转动能级跃迁，伴随能量交换：

  $$\mathrm{\Delta }J=0,\pm 2$$
  • 斯托克斯跃迁（Stokes, $\mathrm{\Delta }J=+2$）：入射光子激发分子进入更高转动激发态，光子损失能量，出射波长发生红移（$\lambda >{\lambda }_0$）；
  • 反斯托克斯跃迁（Anti-Stokes, $\mathrm{\Delta }J=-2$）：处于转动激发态的分子将能量转移给光子回到基态，出射波长发生蓝移（$\lambda <{\lambda }_0$）。

在常温（$300\text{ K}$）下，氮氧分子的转动能级间距较窄，引发的光子频移在紫外-可见波段对应于 $\mathrm{\Delta }\lambda \approx \pm 0.1\sim 2.5\text{ nm}$ 的微小波长位移。

🎢 生动比喻：光子跳楼机与深谷填平效应

太阳本身的夫琅禾费线就像一个极其陡峭昏暗的大峡谷：峡谷两侧是光强极其刺眼的万丈高台（强连续谱），而峡谷底部的光子几乎被太阳本身的原核吸收殆尽，漆黑一片。 当这束日光穿过地球厚厚的氮氧大气时，非弹性旋转拉曼散射就像一台微观“光子跳楼机”：高台上极少一部分（4%）的光子在跳楼机上释放或吸收转动动能，波长发生微小的左右偏移，不偏不倚地掉落并灌注进了漆黑狭窄的峡谷谷底！ 结果：从地面测量时，原本极深的太阳吸收谷底被非物理地“填浅”了 2%~10%！

3.2 吸收线填充（Filling-in of Fraunhofer Lines）与差分伪信号

太阳原本的夫琅禾费吸收线中心本极为昏暗甚至接近零光强，而吸收线两侧则是极明亮的连续谱。

当这束日光穿过地球厚厚的大气层时，吸收线两侧数万倍强度的连续谱光子中，有一小部分（$4\mathrm{\%}$）由于非弹性旋转拉曼散射改变了波长，被“抛掷”并“浇注”进了暗弱的吸收线谷底！

• 宏观测量结果：天空中散射日光的光谱线谷底被非物理地填浅了 $2\mathrm{\%}\sim 10\mathrm{\%}$；

• 对 DOAS 差分拟合的灾难：

  $${\tau }_{\text{apparent}}=\ln \left(\frac{I_{\text{meas}}}{I_{\text{ref}}}\right)$$

  若测量谱与参考谱的光程不同或散射几何不同（例如大太阳天顶角时斜程拉曼散射远强于小天顶角），线深度的填充差异在对数除法中会直接演变为强达 $\mathrm{\Delta }\tau \approx 0.01\sim 0.05$ 的高频假象特征！

• 相比之下，对流层微弱痕量气体（如 $\text{HCHO}$ 或 $\text{HONO}$）的差分光学厚度通常小于 $0.001$。环效应引起的伪差分信号比真实气体信号高出 10 到 50 倍！

---

4. Ring 效应修正谱的严密数学生成与拟合

为了彻底剔除这一假象，现代 DOAS 理论将环效应转化为一个“等效吸收成分”，在反演拟合方程中进行正交消除。

4.1 基于量子拉曼截面的正向物理建模法

根据 Chance & Spurr 以及 Bussemer 等人的理论，针对给定的高分辨率太阳参考通量光谱 $I_0(\lambda )$，考虑 ${\text{N}}_2$（$78.08\mathrm{\%}$）与 ${\text{O}}_2$（$20.95\mathrm{\%}$）在环境温度 $T$ 下的量子转动跃迁几率，建立微观非弹性散射截面张量 ${\sigma }_{\text{RRS}}({\lambda }^{\prime }\to \lambda ,T)$。

计算穿行大气后的纯旋转拉曼非弹性光谱 $I_{\text{Raman}}(\lambda )$：

$$I_{\text{Raman}}(\lambda )={\int }_{-\mathrm{\infty }}^{+\mathrm{\infty }}{I}_0({\lambda }^{\prime })\cdot {\sigma }_{\text{RRS}}({\lambda }^{\prime }\to \lambda ,T)d{\lambda }^{\prime }$$

将纯拉曼散射谱与弹性未充填太阳谱分别与仪器狭缝函数 $W(\lambda )$ 卷积，定义标准无量纲 Ring 差分截面 ${\sigma }_{\text{Ring}}(\lambda )$：

核心公式：标准 Ring 差分参考截面生成

$${\sigma }_{\text{Ring}}(\lambda )\equiv \frac{I_{\text{Raman}}(\lambda )\otimes W(\lambda )}{I_0(\lambda )\otimes W(\lambda )}-P_{\text{poly}}(\lambda )$$

• 分母为仪器狭缝平滑后的弹性太阳参考光谱；
• 分子为经过旋转拉曼非弹性卷积的纯散射光谱；
• $P_{\text{poly}}(\lambda )$ 用于扣除生成的比值谱中随波长平滑变化的低频基底，保留纯快变高频填充结构。

4.2 将 Ring 谱纳入 DOAS 联合拟合方程

在反演软件（如 QDOAS）中，将计算生成的 ${\sigma }_{\text{Ring}}(\lambda )$ 与二氧化氮、臭氧等真实气体的吸收截面并列，作为一个**“等效虚拟吸收截面（Pseudo-absorber）”**加入拟合方程矩阵：

包含 Ring 效应校正项的 DOAS 拟合方程

$$\ln \left(\frac{I({\lambda }_i)}{I_{\text{ref}}({\lambda }_i)}\right)=-\sum _{j=1}^M{\sigma }_j^{\prime }({\lambda }_i)\cdot \mathrm{\Delta }{\text{SCD}}_j-{\alpha }_{\text{Ring}}\cdot {\sigma }_{\text{Ring}}({\lambda }_i)-P_m({\lambda }_i)+{\epsilon }_i$$

• ${\alpha }_{\text{Ring}}$：反演求解的环效应填充权重系数；
• 拟合效果：原本高达 ${10}^{-2}$ 的环效应虚假共轭峰被 ${\alpha }_{\text{Ring}}{\sigma }_{\text{Ring}}$ 项彻底吸纳抵消；
• 残差变化：拟合残差 RMS 瞬间从 ${10}^{-2}$ 暴降至 $3\times {10}^{-4}$，隐藏在深处的极微弱气体吸收指纹得以完美现形！

---

5. 本章总结与下一章预告

本章完整建立了从探测器底层电平校正（Dark / Offset）、太阳夫琅禾费谱线宇宙条形码亚像元波长锁定，到利用量子旋转拉曼散射理论生成 Ring 修正谱的完整数据处理中枢。

有了精准的数据处理链，我们便能全面迈向真实的大气监测战场：如何针对不同污染气体（${\text{NO}}_2,{\text{SO}}_2,\text{HCHO},\text{HONO}$）划分最佳反演窗口？如何利用辐射传输模型将斜柱浓度（dSCD）转化为真正具备环境评估意义的垂直柱浓度（VCD）？又如何通过车载移动走航计算工业园区的宏观排放通量？

请进入最终实战篇：06. 痕量气体立体反演与通量监测。