2、光与生物组织相互作用基础:吸收、散射、折射、反射、透射的基本物理过程及光学参数定义

生物组织光学特性测量的核心,在于理解光子在组织中传播时经历的各类物理事件。当一束光入射到生物组织表面时,并不会简单地沿直线穿透,而是会同时发生反射、折射、散射与吸收。这些过程共同决定了光在组织内的能量分布、传播方向与偏振状态变化。对于从事组织光学建模与测量的工程师而言,必须首先厘清这些基本物理过程的微观机制与宏观表现,并掌握描述它们的光学参数体系。本章将系统梳理光与生物组织相互作用的五大基本过程,明确定义吸收系数、散射系数、各向异性因子、折射率等关键参数,并讨论其在漫反射光谱测量、光学相干层析成像(OCT)等典型应用中的意义。

光与生物组织相互作用——五大物理过程与光学参数 入射光 I₀ 组织表面(空气-组织界面) 镜面反射 R 漫反射 生物组织(散射 + 吸收介质) 折射(n₁→n₂) 散射光 散射 吸收 (能量转化为热/荧光) 透射光 T(衰减后出射) 关键光学参数 μa:吸收系数 (mm⁻¹) μs:散射系数 (mm⁻¹) g:各向异性因子 n:折射率 μs' = μs(1-g): 有效散射系数 μt = μa + μs: 总衰减系数 组织深度 z 光子在组织中经历多次散射与吸收,最终部分透射或反射出组织

上图概括了光与生物组织相互作用的核心物理图景。入射光在组织表面发生反射与折射后进入组织内部,随即在细胞核、线粒体、胶原纤维等微观结构上经历多次散射,同时被血红蛋白、水、黑色素、脂质等发色团吸收。最终,部分光子从组织表面逸出形成漫反射光,部分穿透组织形成透射光。理解这一过程,是建立光在组织中传播的辐射传输方程(RTE)或漫射近似模型的前提。

2.1 核心物理过程与微观机制

2.1.1 吸收(Absorption)

吸收是光能量被组织内特定分子(发色团)转化为其他形式能量(主要是热能,或部分转化为荧光/磷光)的过程。从量子力学角度看,当光子能量恰好匹配分子从基态到激发态的能级差时,光子被吸收。生物组织中主要的吸收体包括:

  • 血红蛋白(HbO₂、Hb):在可见光波段(400-600 nm)有强吸收,是血氧饱和度测量的基础。
  • 水(H₂O):在近红外波段(>900 nm)吸收逐渐增强,限制了光在组织中的穿透深度。
  • 黑色素(Melanin):在紫外-可见波段吸收较强,影响皮肤光学特性。
  • 脂质(Lipid):在近红外波段(930 nm、1200 nm附近)有特征吸收峰。

吸收的宏观度量是吸收系数 μa(单位 mm⁻¹),定义为光子传播单位距离被吸收的概率。根据Beer-Lambert定律,一束强度为 I₀ 的光穿过厚度为 d 的均匀吸收介质后,透射强度为:

I = I₀ · exp(-μa · d)

在组织光学中,μa 的典型值范围从深红波段的 0.01 mm⁻¹(低吸收)到蓝紫波段的 10 mm⁻¹(高吸收)不等。

2.1.2 散射(Scattering)

散射是光子在组织中传播方向发生改变的过程,源于组织内折射率的不均匀性。细胞核(折射率约1.40-1.46)、线粒体(约1.38-1.41)、胶原纤维(约1.47)与细胞质(约1.35-1.37)之间存在折射率失配,形成散射中心。散射的物理机制包括:

  • 瑞利散射(Rayleigh Scattering):当散射体尺寸远小于波长(d << λ)时,散射强度与 λ⁻⁴ 成正比,在短波段占主导。
  • 米氏散射(Mie Scattering):当散射体尺寸与波长相当(d ≈ λ)时,散射具有明显的前向偏好性,是组织散射的主要贡献者。

散射系数 μs(单位 mm⁻¹)定义为单位距离内光子发生散射的概率。生物组织的 μs 通常在 5-50 mm⁻¹ 范围(可见-近红外波段),远大于 μa,因此散射在光与组织相互作用中占据主导地位。

2.1.3 折射与反射(Refraction & Reflection)

当光从空气(n≈1.0)入射到组织表面(n≈1.35-1.45)时,由于折射率不匹配,部分光被反射回空气,部分光折射进入组织。根据菲涅耳公式,垂直入射时的反射率 R 为:

R = ((n₁ - n₂)²) / ((n₁ + n₂)²)

对于空气-组织界面(n₁=1.0, n₂=1.4),R ≈ 2.8%。然而,实际测量中,由于组织表面粗糙且存在漫反射,反射率会更高。折射率 n 不仅影响表面反射,还影响光在组织内的传播速度(v = c/n),进而影响飞行时间测量(如时间分辨光谱)。

2.1.4 透射(Transmission)

透射是光穿过组织后从另一侧出射的过程。对于薄组织样本(如皮肤表皮层、血管壁),透射光携带了组织内部吸收与散射的积分信息。透射率 T 与吸收系数、散射系数及组织厚度密切相关。在漫射近似下,透射光强度可表示为:

T ≈ (1-R)² · exp(-μt · d) · 扩散因子

其中 μt = μa + μs 为总衰减系数,扩散因子反映了多次散射对光程的延长效应。

2.2 光学参数定义与测量方法

参数 符号 单位 物理意义 典型值(皮肤,650nm)
吸收系数 μa mm⁻¹ 单位距离内光子被吸收的概率 0.01-0.05
散射系数 μs mm⁻¹ 单位距离内光子发生散射的概率 10-40
各向异性因子 g 无量纲 散射角余弦的平均值,g=0各向同性,g→1前向散射 0.8-0.95
有效散射系数 μs' mm⁻¹ μs' = μs(1-g),等效各向同性散射 1-3
折射率 n 无量纲 光在组织中传播速度的减慢因子 1.35-1.45
总衰减系数 μt mm⁻¹ μt = μa + μs,单位距离内光子的总消失概率 10-40

各向异性因子 g 是描述散射方向分布的关键参数,定义为散射角 θ 的余弦平均值:g = <cosθ>。当 g=0 时,散射为各向同性;当 g=0.9 时,光子倾向于以较小角度偏转(前向散射)。生物组织的 g 值通常在 0.7-0.99 之间,这意味着单次散射事件中光子方向改变不大,但经过多次散射后,光子的传播方向仍会趋于随机化。有效散射系数 μs' 将各向异性散射等效为各向同性散射,是漫射近似理论中的核心参数。

2.3 测量方法概述与案例拆解

2.3.1 积分球法测量反射率与透射率

积分球法是最经典的组织光学参数测量方法。其原理是:将组织样本置于积分球入口,利用球内高反射涂层(如硫酸钡,反射率>99%)收集所有方向上的反射光或透射光,通过对比参考标准计算总反射率 R_d 和总透射率 T_t。结合Beer-Lambert定律与逆蒙特卡洛算法,可反推出 μa 和 μs。

案例:猪皮组织光学参数测量

背景与约束:某实验室需要获取猪皮在 632.8 nm(He-Ne激光)下的 μa 和 μs,用于光动力疗法剂量计算。样本厚度为 2.1 mm,表面不平整,需避免水分蒸发影响。

操作流程:

  1. 将猪皮样本夹持在积分球样本架上,确保光斑直径(5 mm)完全落在样本范围内。
  2. 分别测量总反射率 R_d = 0.42,总透射率 T_t = 0.08(已扣除空气-玻璃界面反射)。
  3. 使用逆蒙特卡洛算法(MCML程序),输入 R_d、T_t、样本厚度 d=2.1 mm、折射率 n=1.4,迭代拟合 μa 和 μs。
  4. 拟合结果:μa = 0.023 mm⁻¹,μs = 12.5 mm⁻¹,g 预设为 0.85(文献值)。

结果与验证:将拟合参数代入正向蒙特卡洛模拟,计算得到的 R_d 和 T_t 与实测值偏差小于 2%。后续将该参数用于光剂量计算,预测的坏死深度与动物实验误差在 ±15% 以内。

风险与复盘:样本脱水会导致 μa 和 μs 漂移,需在 30 分钟内完成测量;表面反射对 R_d 影响显著,需用折射率匹配液(如甘油)消除表面镜面反射干扰。

2.3.2 空间分辨漫反射法

对于在体测量,空间分辨漫反射法更为实用。该方法通过光纤探头在组织表面以不同源-探距离(如 0.5-10 mm)采集漫反射光强度,利用漫射近似模型拟合得到 μa 和 μs'。其核心关系式为:

R(ρ) = (1/(4πμt')) · [ (1/ρ) · exp(-μeff·ρ) + (1/μt') · exp(-μeff·ρ) ]

其中 μeff = √(3μa·μt') 为有效衰减系数,ρ 为源-探距离。通过测量多个 ρ 下的 R(ρ),可同时拟合 μa 和 μs'。

2.4 常见误区与工程注意事项

误区一:混淆 μs 与 μs'

在漫射近似模型中,真正起作用的是 μs' 而非 μs。若将 μs 直接代入漫射方程,会导致有效穿透深度计算严重偏差。例如,对于 μs=20 mm⁻¹、g=0.9 的组织,μs'=2 mm⁻¹,两者相差 10 倍。在阅读文献或使用软件时,务必确认参数定义。

误区二:忽略折射率失配的影响

组织-空气界面的折射率失配不仅影响表面反射,还影响内部光子的边界逃逸概率。在蒙特卡洛模拟中,若未正确设置折射率边界条件,会导致漫反射率计算偏差达 10-20%。建议在模拟中显式设置 n_tissue=1.4,n_air=1.0。

经验提醒:

  • 测量 μa 时,尽量选择吸收较弱且散射主导的波段(如 650-800 nm),以降低散射对吸收拟合的干扰。
  • 对于强散射组织(如皮肤真皮层),使用时间分辨或频域方法比稳态方法更可靠,因为可以分离吸收与散射对时间点扩散函数的影响。
  • 在体测量时,探头压力会改变局部血流,导致 μa 变化。建议使用轻接触或非接触探头,并记录施加压力值。

总结复盘

本章系统梳理了光与生物组织相互作用的五大基本物理过程。吸收与散射是决定光在组织内部分布的核心机制,折射与反射控制着光在边界处的耦合效率,透射则提供了组织内部信息的另一种观测窗口。光学参数 μa、μs、g、n 构成了描述这些过程的完整参数集,其中 μs' 和 μt 是工程计算中常用的导出量。在实际测量中,积分球法适用于离体样本的精确标定,空间分辨漫反射法适用于在体快速测量,而时间分辨方法则能提供更丰富的散射与吸收分离信息。理解这些参数的物理含义与测量边界条件,是后续进行光传输建模、光学断层成像和光剂量计算的基础。