4、光在组织中的传输理论：辐射传输方程（RTE）的推导、假设与适用范围

在生物组织光学特性测量与建模中，辐射传输方程（Radiative Transfer Equation, RTE）是描述光在散射、吸收介质中传播行为的核心理论框架。与麦克斯韦方程组直接求解电磁场相比，RTE 忽略了光的波动性（相干效应），将光视为大量无相互作用的粒子（光子），通过能量守恒关系建立微分方程。这一简化使得 RTE 在生物组织这类强散射介质中具有极高的实用价值——组织对光的散射远强于吸收，光子经历多次散射后相位信息迅速丢失，波动效应不再主导能量输运。本章将从物理图像出发，系统推导 RTE 的数学形式，明确其关键假设，并讨论其在组织光学中的适用范围与边界条件。

RTE 推导逻辑与假设边界 基本辐射量 L(r, ŝ, t) — 辐射亮度 [W·m⁻²·sr⁻¹·nm⁻¹] μₐ — 吸收系数 [mm⁻¹] μₛ — 散射系数 [mm⁻¹] p(ŝ, ŝ′) — 散射相函数 n — 组织折射率 (~1.33-1.40) 辐射传输方程 (RTE) ∂L(r,ŝ,t)/∂t · (1/c) + ŝ·∇L(r,ŝ,t) = -(μₐ+μₛ)·L(r,ŝ,t) + μₛ·∫p(ŝ,ŝ′)L(r,ŝ′,t)dΩ′ + S(r,ŝ,t) 能量守恒 · 粒子流平衡 核心假设 ① 无相干效应（非波动） ② 散射粒子独立、稀疏 ③ 折射率均匀（宏观） ④ 偏振不敏感 ⑤ 时间尺度远大于相干时间 ⑥ 边界无反射（或单独处理） 推导路径：从能量守恒到微分方程 步骤1 控制体积元dV 步骤2 通量平衡分析 步骤3 散射积分项 步骤4 源项与边界条件 适用范围与限制 ✓ 适用 扩散层(μₛ′≫μₐ)、近红外 △ 需修正 各向异性强(g>0.9)、薄层 ✗ 不适用 相干成像、超短脉冲(

4.1 RTE 的物理基础与推导过程

4.1.1 核心原理：从辐射度学到输运理论

RTE 的推导始于辐射度学（Radiometry）的基本量定义。辐射亮度 L(r, ŝ, t) 表示在位置 r 处、沿方向 ŝ、单位立体角、单位投影面积、单位时间、单位波长间隔内通过的光功率。在生物组织中，光与组织的相互作用主要表现为吸收（absorption）和散射（scattering）。吸收将光能转化为热能或其他形式的能量，散射则改变光子的传播方向而不改变其能量（弹性散射）。RTE 的本质是光子数守恒方程：在任意微小体积元内，光子的增加率等于流入率减去流出率，再减去吸收损失，加上散射进入该方向的光子，最后加上内部源项。

推导时，我们考虑一个以 r 为中心、体积为 dV 的微小六面体，沿 ŝ 方向取截面 dA、长度 ds 的圆柱体。在时间 dt 内，从 ŝ 方向进入该体积的光子能量为 L(r, ŝ, t)·dA·dΩ·dt。经过 ds 距离后，由于吸收和散射，光强衰减为 L(r+ŝ·ds, ŝ, t+dt)。衰减量正比于消光系数 μₜ = μₐ + μₛ 和传播距离 ds。同时，来自其他方向 ŝ′ 的光子可能被散射到 ŝ 方向，其贡献为 μₛ·∫p(ŝ, ŝ′)·L(r, ŝ′, t)dΩ′·ds·dA·dΩ·dt，其中 p(ŝ, ŝ′) 是散射相函数，描述光子从 ŝ′ 方向散射到 ŝ 方向的概率密度。将上述各项代入能量守恒关系，取极限 ds→0、dt→0，即可得到 RTE 的微分形式。

RTE 标准形式（时间相关）：

∂L(r,ŝ,t)/(c·∂t) + ŝ·∇L(r,ŝ,t) = -(μₐ+μₛ)·L(r,ŝ,t) 
                                  + μₛ·∫₄ₚ p(ŝ,ŝ′)·L(r,ŝ′,t)dΩ′ 
                                  + S(r,ŝ,t)

其中 c 为组织中光速（c = c₀/n），S 为内部源项（如荧光、磷光）。稳态时 ∂L/∂t = 0，方程简化为不含时间导数的形式。

4.1.2 散射相函数与各向异性因子

散射相函数 p(ŝ, ŝ′) 是 RTE 中描述散射方向分布的关键参数。生物组织中，最常用的相函数是 Henyey-Greenstein (HG) 函数：

p(θ) = (1-g²) / [4π(1+g²-2g·cosθ)^(3/2)]

其中 θ 是散射前后方向的夹角，g = ⟨cosθ⟩ 是各向异性因子，表示散射的平均余弦值。生物组织的 g 通常在 0.8~0.95 之间（如皮肤 dermis 约 0.81，脑白质约 0.90），意味着散射强烈前向。HG 函数虽然形式简单，但在后向散射描述上存在偏差，因此有时使用双 HG 函数或 Modified HG 函数以提高拟合精度。各向异性因子的引入使得我们可以定义约化散射系数 μₛ′ = μₛ·(1-g)，它表征各向同性等效散射能力，是扩散理论中的核心参数。

4.1.3 边界条件与折射率失配

在实际组织测量中，光从空气（n≈1.0）入射到组织（n≈1.33~1.40）表面时，会发生菲涅耳反射和折射。RTE 的边界条件需要处理折射率失配问题。常用的边界条件包括：

  • 匹配边界：假设组织与外界折射率相同，边界处辐射亮度连续，适用于浸入式测量（如组织液匹配液）。
  • 失配边界：考虑菲涅耳反射系数 R(θ)，边界处满足 L_out = R(θ)·L_in，内反射光重新进入组织。对于漫反射测量，通常采用折射率失配因子 m = n_org/n_air 来修正边界反射。
  • 外推边界：在扩散近似中，常引入外推长度 z_b = 2·D·(1+R_eff)/(1-R_eff)，其中 R_eff 是有效反射系数，D 是扩散系数。

4.2 RTE 的假设条件与适用范围

4.2.1 核心假设的物理内涵

RTE 的推导基于若干关键假设，理解这些假设是正确应用 RTE 的前提。下表总结了各假设的物理含义及其在生物组织中的有效性：

假设条件 物理含义 生物组织中的有效性
无相干效应 忽略干涉、衍射，光强直接相加 组织强散射，相干长度极短（<10fs），有效
散射粒子独立 散射事件之间无关联，散射截面线性叠加 细胞器间距远大于波长，近似成立
折射率均匀 宏观折射率恒定，不考虑微观折射率起伏 组织平均折射率稳定，但细胞膜处有局部变化
偏振不敏感 忽略偏振态变化，标量方程 多次散射后偏振信息丢失，但浅层测量需注意
时间尺度 光脉冲宽度远大于相干时间 飞秒脉冲需考虑相干效应，皮秒以上适用
无源项 除边界入射外，组织内部无自发辐射 荧光测量需额外加入源项 S

4.2.2 判断方法：何时使用 RTE 而非其他模型

在实际测量中，选择 RTE 还是其近似形式（如扩散方程、蒙特卡罗模拟）取决于光学厚度和测量几何。以下判断流程可帮助决策：

  1. 光学厚度评估：计算光学厚度 τ = μₜ·L（L 为光程长度）。若 τ < 1，为光学薄介质，光子散射次数少，RTE 可简化；若 τ > 10，为光学厚介质，扩散近似可能足够。
  2. 散射-吸收比：计算 μₛ′/μₐ。若比值 > 10，扩散占主导，可用扩散方程；若比值在 1~10 之间，需使用完整 RTE 或蒙特卡罗；若比值 < 1，吸收主导，需考虑 Beer-Lambert 定律修正。
  3. 测量几何：对于透射测量（如组织切片），若厚度远大于输运平均自由程（1/μₜ′），扩散近似有效；对于反射测量（如皮肤漫反射），源-探测器距离需大于 3~5 个输运平均自由程。
  4. 时间分辨需求：若需要时间分辨信息（如 TOF 测量），需使用时间相关 RTE 或蒙特卡罗；若仅需稳态强度，可简化。

经验法则：在近红外窗口（600~1300 nm），典型皮肤组织的 μₐ ≈ 0.01~0.1 mm⁻¹，μₛ′ ≈ 0.5~2 mm⁻¹。对于源-探测器距离大于 5 mm 的反射测量，扩散方程误差小于 5%；对于小于 2 mm 的距离，必须使用完整 RTE 或蒙特卡罗模拟。

4.2.3 案例拆解：皮肤漫反射测量的模型选择

背景与目标：某研究团队欲开发一款便携式皮肤光谱仪，通过测量 600~1000 nm 波段的漫反射光谱来估算皮肤黑色素和血氧含量。探头采用光纤束结构，源-探测器中心距为 1.5 mm 和 3.0 mm 两个通道。需要选择合适的光传输模型进行反演。

初始参数评估：皮肤表皮层厚度约 0.1 mm，真皮层约 1~2 mm。在 700 nm 处，表皮 μₐ ≈ 0.02 mm⁻¹（黑色素吸收），μₛ′ ≈ 1.2 mm⁻¹；真皮 μₐ ≈ 0.03 mm⁻¹（血液吸收），μₛ′ ≈ 0.8 mm⁻¹。输运平均自由程 l* = 1/(μₐ+μₛ′) ≈ 0.8~1.2 mm。

分析判断：对于 1.5 mm 源-探测器距离，该距离约为 1.5~2 个输运平均自由程，处于扩散近似的临界区域。若直接使用扩散方程，误差可能达到 10%~15%。对于 3.0 mm 距离，约为 3~4 个 l*，扩散近似误差可控制在 5% 以内。因此，团队决定：

  • 对于 3.0 mm 通道，采用扩散方程进行快速反演，用于实时监测。
  • 对于 1.5 mm 通道，采用蒙特卡罗模拟生成查找表（LUT），结合 RTE 的离散坐标法进行精确求解。
  • 在反演算法中，将两个通道的数据融合，利用短距离通道对浅层（表皮）敏感、长距离通道对深层（真皮）敏感的特性，实现分层光学参数提取。

结果与复盘：通过对比 20 例志愿者前臂皮肤测量数据，融合算法提取的黑色素含量与皮肤镜评估结果相关系数 r=0.91，血氧饱和度与血氧计测量值偏差小于 3%。复盘发现，短距离通道对表皮厚度变化敏感，若表皮厚度超过 0.15 mm，需引入双层模型修正。此外，探头压力会导致局部血流变化，需在硬件设计中加入压力传感器进行实时校正。

4.3 RTE 的数值求解方法

4.3.1 离散坐标法（SN 近似）

离散坐标法将方向角 ŝ 离散为有限个方向（通常为 S₂、S₄、S₈ 等），将 RTE 转化为一组耦合的偏微分方程。每个方向对应一个权重 wᵢ，积分项用数值求积近似。该方法在光学厚度适中的介质中精度较高，但存在射线效应（ray effect）——即方向离散不足时，点源附近会出现虚假的射线状条纹。生物组织中，由于散射较强，射线效应通常不明显，但需注意边界附近的伪影。

4.3.2 蒙特卡罗模拟

蒙特卡罗（MC）方法通过随机抽样模拟大量光子的随机行走过程，是生物组织光学中最灵活的求解方法。其优势在于：可处理任意复杂几何、任意相函数、多层结构、边界反射等。标准 MC 流程包括：光子发射、步长抽样（基于 μₜ）、散射方向抽样（基于相函数）、边界处理、吸收权重更新、探测统计。MC 的缺点是计算量大，但通过 GPU 并行和方差缩减技术（如重要抽样、分裂-俄罗斯轮盘赌），可在秒级完成典型皮肤模型的模拟。

MC 模拟关键参数设置（以皮肤为例）：

# 组织分层参数（700 nm）
# 层1: 表皮 厚度0.1mm, μₐ=0.02, μₛ=4.0, g=0.85, n=1.40
# 层2: 真皮 厚度1.5mm, μₐ=0.03, μₛ=6.0, g=0.90, n=1.37
# 层3: 皮下 厚度3.0mm, μₐ=0.01, μₛ=8.0, g=0.95, n=1.35

# 光子数: 10^6 ~ 10^7
# 探测半径: 0.5mm 环形探测器
# 时间门: 0~2000 ps, 分辨率 10 ps

4.3.3 扩散近似及其边界

当散射占主导且探测距离足够远时，RTE 可简化为扩散方程。扩散近似的核心假设是辐射亮度在角度上近似各向同性，可展开为 L(r, ŝ) ≈ (1/4π)·Φ(r) + (3/4π)·J(r)·ŝ，其中 Φ 是光通量率，J 是光子流矢量。代入 RTE 后得到：

-∇·(D∇Φ(r)) + μₐ·Φ(r) = S₀(r)

其中 D = 1/[3(μₐ+μₛ′)] 是扩散系数。扩散方程是椭圆型偏微分方程，可用有限元法（FEM）或有限差分法高效求解。但需注意，扩散近似在以下情况失效：① 源-探测器距离小于 3~5 个 l*；② 吸收系数过高（μₐ/μₛ′ > 0.1）；③ 介质中存在低散射区域（如脑脊液层）；④ 时间分辨测量中早期到达的光子（弹道光和蛇形光）。

4.4 常见误区与工程实践要点

4.4.1 误区一：忽略折射率失配的边界效应

许多初学者在应用 RTE 时，直接采用匹配边界条件，忽略了组织-空气界面的菲涅耳反射。实际上，对于 n=1.37 的组织，垂直入射时的反射率约为 2.4%，但大角度入射时反射率急剧上升。在漫反射测量中，边界反射会导致出射光通量减少 10%~20%，若不修正，反演得到的 μₛ′ 会偏低。工程实践中，应在 RTE 求解器中内置折射率失配边界处理，或使用外推边界条件。

4.4.2 误区二：将 HG 相函数视为万能

HG 相函数虽然广泛使用，但其后向散射概率偏低。对于需要精确描述后向散射的测量（如背向漫反射），建议使用双 HG 函数：p(θ) = α·HG(g₁,θ) + (1-α)·HG(g₂,θ)，其中 g₁ 取 0.9 左右描述前向峰，g₂ 取 -0.2~-0.5 描述后向分量，α 为权重。研究表明，双 HG 函数对组织散射的拟合精度可提高 30% 以上，尤其在 120°~180° 散射角范围。

4.4.3 误区三：混淆 μₛ 与 μₛ′

μₛ 是散射系数，表示单位路径长度上发生散射事件的概率；μₛ′ = μₛ(1-g) 是约化散射系数，表示等效各向同性散射能力。在扩散方程中，使用的是 μₛ′ 而非 μₛ。若将 μₛ 直接代入扩散方程，会导致扩散系数偏小，光通量衰减过快。在文献对比时，务必确认参数是 μₛ 还是 μₛ′，以及对应的 g 值。

工程风险提示：在组织光学仪器的设计验证中，建议使用蒙特卡罗模拟作为“金标准”来验证 RTE 或扩散方程求解器的准确性。具体做法：构建与实验一致的三维几何模型，设置相同的光学参数，分别用 MC 和 RTE 求解器计算漫反射率，对比两者差异。若偏差超过 5%，需检查边界条件、相函数选择或离散精度。此外，所有光学参数应注明波长和测量条件，因为 μₐ 和 μₛ′ 随波长变化显著。

4.4.4 总结复盘：RTE 应用决策清单

在实际项目中，建议按以下清单进行 RTE 应用决策：

  1. 明确测量目标：是提取光学参数（μₐ, μₛ′）还是重建光学图像（DOT）？前者对精度要求高，后者更关注空间分辨率。
  2. 评估光学条件：计算 μₛ′/μₐ 比值和光学厚度，判断是否满足扩散近似条件。
  3. 选择求解方法：若需精确解且计算资源充足，选 MC；若需实时反演，选扩散方程或降阶 RTE；若需处理复杂边界，选 FEM 求解扩散方程。
  4. 验证与标定：使用组织模拟液（Intralipid 溶液 + 印度墨水）进行实验验证，确保模拟与实测一致。
  5. 敏感性分析：分析各光学参数对测量信号的敏感性，确定反演算法的稳定性和唯一性。
← 上一章 📖 返回目录 下一章 →