瑞利分馏模型（Rayleigh fractionation model）是稳定同位素地球化学中最经典、最广泛应用的开放系统分馏模型。它描述了一个有限储库（reservoir）在连续移除产物（product）时，反应物（reactant） 同位素组成随反应进展（或剩余分数 f）发生指数式变化的过程。

该模型最初由Lord Rayleigh于1902年为分馏蒸馏推导，后来广泛用于稳定同位素领域，如：

• 水汽从云中逐步凝结降水（雨出过程，导致剩余水汽越来越富集重同位素）

• 蒸发过程（剩余液态水富集1⁸O和2H）

• 有机污染物或烃类在微生物降解中的动力学分馏

• 矿物从溶液中逐步沉淀

• 冰盖形成过程中的海水同位素演化

核心假设（最严格形式）

1. 系统为开放系统：产物一旦生成，就立即从系统中移除，不再与剩余反应物发生同位素交换（无回流）。

2. 分馏因子 α（或 ε）在整个过程中恒定（温度、压力等条件不变）。

3. 反应物储库均匀混合（well-mixed）。

4. 过程为连续的、微分式的（instantaneous removal）。

基本推导与公式

考虑反应物 R（剩余相），产物 P（瞬时移除相）。

瞬时分馏：α = R_P / R_R  （重/轻同位素比值之比，通常 α < 1 表示轻同位素优先进入产物）

在微小时间 dt 内，移走 dN 个分子（N 为剩余分子数），其中重同位素移走 dN_h = (R_P) × dN = α R_R dN

剩余重同位素变化：d(R_R N) = - α R_R dN

展开：R_R dN + N dR_R = - α R_R dN

整理得：dR_R / R_R = (α - 1) dN / N

积分两边（从初始 R₀、N₀ 到当前 R、N）：

ln(R / R₀) = (α - 1) ln(f)

其中 f = N / N₀ = 剩余分数（0 < f ≤ 1）

因此经典瑞利方程：

R = R₀ f^(α - 1)

或写成 δ 值形式（最常用，适用于分馏较小时近似）：

δ = [δ₀ + 1000] f^(α - 1) - 1000 ‰

（精确形式需用 (1000 + δ)/1000 代入 R）

三条典型曲线（同一张图中最常见的表达）

• 剩余反应物（residual reactant）：δ 随 f 减小而指数上升（越来越富集重同位素）

• 瞬时产物（instantaneous product）：δ_P = δ_R + 1000(α - 1) ≈ δ_R + ε （始终比当前剩余反应物轻 ε ‰）

• 累积产物（accumulated / integrated product）：从 f=1 到当前 f 的所有瞬时产物的质量加权平均，通常最轻，随 f 减小而逐渐变重但始终轻于剩余反应物。

实际应用举例

过程	α 示例（25℃）	ε (‰)	典型行为	意义
水 → 水蒸气 (1⁸O)	~1.0094	~+9.4	剩余水 δ1⁸O 上升	解释云中雨出导致降水 δ1⁸O 随高度降低
水 → 水蒸气 (2H)	~1.080	~+80	剩余水 δD 急剧上升	极地冰芯 δD 极负
有机物降解 (13C)	1.005–1.030	+5～+30	剩余底物 δ13C 变重	估算降解程度 B = (δ_t - δ₀)/ε
甲烷生成 (13C)	~1.060–1.080	+60～+80	产物 CH₄ 极轻	区分生物 vs 热成因甲烷

局限性与扩展

• α 必须恒定（实际中温度变化、动力学 vs 平衡会变）

• 无混合、无回流（真实大气中常有再蒸发、混合，需要修正模型）

• 仅适用于单一过程（多步需分段或耦合模型）

• 现代扩展：考虑扩散限制、表面反应动力学、衰变链等（如污染物降解链）

一句话总结：瑞利分馏模型用一个简单的指数方程，完美捕捉了开放系统中“轻同位素先跑、重同位素滞留”的指数放大效应，成为解释自然界从云雨到污染物降解、从冰芯到石油成熟度等无数“偏差放大”故事的基石工具。

转载本文请联系原作者获取授权，同时请注明本文来自李中平科学网博客。
链接地址：https://blog.sciencenet.cn/blog-3549522-1525922.html

上一篇：为什么稳定同位素地球化学本质上是“偏差的科学”，都围绕“偏差如何随过程放大”展开