光在生物组织传输的蒙特卡罗时域模拟.pdf

第 1期 (总第 200期) 机 械 工 程 与 自动 化 No．1 2017年 2月 MECHANICAL ENGINEERING 8L AUTOMAT10N Feb． 文章编号 ：1672—6413(2017)01—0025—02 光在生物组织传输的蒙特卡罗时域模拟 王玉娇 ，王武义 (哈尔滨工业大学 机 电工程学院，黑龙江 哈尔滨 150001) 摘要 ：对现有的蒙特卡罗算法提 出了改进 ，可 以将其应用于生物组织光传输的时域模拟 中；并采用漫射近似 理论模型与仿真结果进行 了对比验证，对其现象和差异进行 了详细的讨论与分析。结果表 明：该蒙特卡罗算 法可 以有效地模拟 出漫反射光谱，且具有较高的精度，对生物组织无损检测具有很高的参考价值 。 关键词：蒙特卡罗模拟；生物组织；无损检测 ；漫射近似理论 中图分类号：TP391．7 文献标识码 ：A 0 引言 布满足柱坐标对称 ，如图2所示 。 在生物组织光子学中，其核心 问题是有效地建立 组织光学模型，即从组织光学特性参数出发 ，结合相应 的数学和物理知识 ，得到 比较准确地反映光与生物组 I发射光子 ’ 织相互作用的解析解或数值解等。然而考虑到组织高 l随机步长 散射随机性 的特点，理论模型往往和实际组织相差较 一 大，不能真实地反映介质 内部的特征信息。本文基于 蒙特卡罗方法，模拟了生物组织的漫反射光谱 ，它可以 厂 I 移动一个步 真实反映出光子在介质 内部的运动情况，这种算法为 生物组织无损检测提供了参考 。 竿圜 1 蒙特卡罗时域算法 l移动到t边界l l 哩! 蒙特卡罗算法名字来源于摩纳哥赌城 Monte Carlo，这是一种概率统计方法 ，其基本思想是在计算 机上实现若干次独立重复实验，统计仿真结果 即可确 定 问题的解L1]。本文将其应用于生物组织的光传输问 Y Y 题 ，其原理是根据现有的传输理论建立一个完备的光 传输概率模型，并将光束等效成若干个光子，通过抽样 ———[二二二 统计每个光子在组织 内的行走轨迹 ，即可得到一个完 整 的光分布模型。 N 蒙特卡罗时域算法流程图如图 1所示 。首先在组 织表面发射一个光子 ，并分配其随机步长 s；光子前进 结束 该步长后，根据概率模型会随机发生散射作用和吸收 图 1 蒙特卡 罗时域算法流程 图 作用 ，散射作用会改变光子的前进方向，吸收作用会降 我们用光通量 (．0，z，)表征组织 内部 的光强分 低光子的权值；光子继续前进一个随机步长 ，并再次发 布情况，其中，P为径 向距离，为组织深度 ，t为光子运 生散射作用和吸收作用 ；根据随机步长长度以及光在 动的任意时刻 ，且坐标原点为光束入射点。此时光通 介质中