一种模拟生物组织内光传播的三维几何蒙特卡洛方法

提出一种几何蒙特卡洛方法(GMC),利用光子位置与物质界面间的几何关系在整个计算区域而非单个网格内计算光子传输。计算区域无网格离散,光在物质界面处的传播严格依据其实际过程进行,消除了网格蒙特卡洛方法(VMC)中的光传播误差并大幅提高了计算速度。对于包含单根血管的情况,采用10μm网格,GMC的计算速度约为VMC的25倍。利用GMC方法计算了激光在多根同轴血管簇皮肤组织模型中的能量沉积,发现血液体积分数不变时,血液能量吸收特性对血管分布的依赖性随着血管数量的增加而降低。当血管数量为20根时,不同血管分布下血液能量吸收最大变动率不超过4%。这表明可以通过假设的血管分布开展研究而无需考虑真实的复杂血管结构,具有重要的实际应用意义。