基于扩散近似的参与性介质光学特性反演研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093503 基于扩散近似的参与性介质光学特性反演 王圣刚，阮立明，齐宏，王希影 （哈尔滨工业大学航空航天热物理研究所，黑龙江，哈尔滨150001） (Tel: 0451E-mail: wangsghit@) 摘要：本文基于扩散近似方程，模拟了短脉冲激光在二维非均匀介质中的传输。短脉冲激光在边界处 射入参与性介质，利用多个位置测量的反射和透射信号，结合共轭梯度法重构了介质内的物性及其分 布。计算表明该扩散近似模型能够模拟出光通量的空间分布，而且通过反演计算，较为准确模拟了非 均匀介质内光学参数的变化。 关键词：扩散近似方程；短脉冲激光；共轭梯度法 0 前言 短激光脉冲已经广泛用于生物领域，如激光组织消融，激光手术，光动力学治疗， 光学成像等。由于生物组织的一些生理特性和病变均与组织的参数密切相关，并且活体 组织的光学参数如吸收系数、散射系数等难于直接测量，因此如何获得活体组织的光学 特性参数是一个十分重要的课题。研究表明：组织表面漫反射光（指经过组织内部散射， 再由组织表面透射的光）与组织光学特性有关，并携带有丰富的组织结构信息[1] 。通过 测量组织表面漫反射光的空间分布来反推组织的光学特性参数，是一种可行的方法。 早期基于短脉冲激光的反问题主要集中在对介质内光学特性（如吸收系数、散射系 数等）的反演及成像研究上，均采用扩散近似模型作为其正向计算来模拟介质内光子的 传输。Saquib等[2]讨论了由时间分辨测量扩散光信号重建未知的吸收和散射系数问题。 使用时域有限差分方法对用于测量的扩散模型正问题进行求解。逆问题以贝叶斯公式为 基础。将最大后验估计（MAP-maximum a posteriori ）用于重建计算。Davies等[3]基于反 差分方法，用有限元法对人体组织光子扩散反演问题进行了研究。Arridge等[4]给出一种 使用伴随源项格式直接得到目标函数梯度的简单方法，利用时域信号的Laplace变换数据 和共轭梯度法（CG ）对圆形区域内的吸收系数和散射系数进行重构，并比较CG法、最 速下降法（SD）和代数重构法（ART ）的重构效果，指出CG法效果最佳。高峰等[5]介绍 了基于有限元方法求解作为光子传输模型的扩散方程的图像重建算法。王加俊等[6, 7]在文 献[8]中给出的联合差分算法基础上，研究了针对内部点光学参量的求导方法，给出了一 种基于树形结构的对内部点光学参量求导的策略，但其正模型使用的是稳态辐射传输方 程，丢失了测量信号随时间变化的信息。 近几年来，人们基于完整的辐射传输方程，发展了很多新的方法，如反差分和伴随 矩阵模型等。与广泛使用的扩散近似相比，这些新方法能够改进的重建图像质量，而且 能求解包含低散射区域的介质重建问题。但他们多数使用的是稳态辐射传输方程，丢失 了测量信号随时间变化的信息，而且计算需要耗费大量机时，因此在实际临床应用方面 存在困难，仍需要加以改进。由于大多数生物组织对近红外光表现出高散射、低吸收的 特性，因此光在传播到探测器前经历了多次散射。当吸收系数与约化散射系数相比很小 资助项目：国家自然科学基金 (No. 哈尔滨工业大学优秀青年教师培养计划资助项目 (HITQNJS.2008.021) 或散射长度远小于试验样品的尺寸时，光在样品中的传播及光通量空间分布可以用扩散 方程描述。与Monte Carlo模拟方法相比，虽然扩散近似在某些情况下计算误差较大，但 由于扩散方程计算速度快，物理意义明确，仍然在组织光学模拟及其反演中广泛应用。 本文基于扩散近似方程，模拟了短脉冲激光在二维非均匀介质中的传输。短脉冲激 光在边界处射入参与性介质，利用多个位置测量的反射和透射信号，结合共轭梯度法重 构了介质内的物性及其分布。计算表明该扩散近似模型能够模拟出光通量的空间分布， 而且通过反演计算，较为准确模拟了非均匀介质内光学参数的变化，这就为扩散近似用 于医学成像研究奠定了理论基础。 1 扩散近似模型 扩散（漫射）方程模型（diffusion equation ）是对辐射传输方程的一种简化近似，由 于它最终表示为相