光学特性参数对生物组织光分布的影响

肖 郑 颖

(1.莆田学院 机电工程学院，福建 莆田351100；2.福建省激光精密加工工程技术研究中心，福建 莆田351100)



【自然科学基础理论研究】

光学特性参数对生物组织光分布的影响

肖 郑 颖1,2

(1.莆田学院 机电工程学院，福建 莆田351100；2.福建省激光精密加工工程技术研究中心，福建 莆田351100)

利用蒙特卡罗方法模拟准直光束在生物组织中的光分布，分别改变吸收系数、散射系数、衰减系数等光学特性参数，比较分析组织中吸收能量密度、光能流率、漫反射率、透射率、光通量分布的差异。结果表明:增大吸收系数和散射系数均会增大吸收能量密度及其梯度；吸收系数影响光能流率和光通量的大小，散射系数影响二者的衰减梯度；吸收系数和散射系数影响漫反射率和透射率大小，散射系数影响二者的衰减梯度。

组织光学；蒙特卡罗模拟;光学特性参数；光分布

光在生物组织内传输的过程中会发生吸收、散射、反射和透射，但是不同的生物组织结构、形态、成分差异很大，会表现出截然不同的光学特性，这种差异也是临床上实现光学诊断病变组织的依据。[1]光学成像技术是光学方法在生物医学领域应用的重要方向，发展出了光学共焦显微成像、荧光成像、光声成像、光学相干层析成像(Optical Coherence Tomography,OCT)技术等等，与传统的成像技术相比，具有非侵入式、 无损、精度高的优点。

OCT技术通过探测样品的后向散射光实现光学成像，可以通过提取OCT信号中的衰减情况实现对生物组织状态的检测。[2-3]因此，研究生物组织中各参量的衰减情况可以为光学诊断判断生物组织状态提供参考和依据。

蒙特卡罗方法通过光子随机行走的方式模拟光在组织中的传输，已经成为组织光学领域验证其他理论或模型的一种非实验标准。[4]

本文通过蒙特卡罗方法模拟无限窄光束辐照吸收系数、散射系数和衰减系数不同的生物组织，分析生物组织中吸收能量密度、光能流率、漫反射率、透射率、光通量的差异，并分析其产生原因。

1 蒙特卡罗模拟

蒙特卡罗模拟在依据光学特性参数建立生物组织模型的基础上，通过跟踪统计大量光子在介质中的光程，估计出相关物理量。

1.1选取随机变量

1.2各随机变量

图1 移动球坐标系示意图

1.3光子传输模拟

光子垂直组织表面入射后，在组织内部每移动一个步长s，其权重、坐标和方向余弦的变化如表1所示。[5]光子初始权重ω=1，当传输至ω<ωth时，由“轮盘赌”判断是否继续传输。

表1 光子移动一个步长，权重、坐标、方向余弦的变化情况

2 模拟方法

蒙特卡罗方法模拟组织内光分布的流程图如图2所示。模拟初始值采用633 nm下皮肤组织的光学特性参数[6]：μa=2.7 cm-1，μs=187 cm-1，g=0.81。另外,模拟光子数为105个，d=0.1 cm，dz=0.001 cm，dr=0.005 cm,ωth=0.000 01。如表2所示，分别改变μa、μs、μt，获得不同参数下无限窄准直光束入射皮肤组织，吸收能量密度、光能流率、漫反射率、透射率、光通量的分布(图3至图7)，分析光学特性参数对组织中光子传输的影响。图3至图6中的曲线1至曲线5分别对应表2中相应序号的光学特性参数。序号2是将吸收系数改变为初始值的2倍，散射系数不变；序号3是将散射系数改变为初始值的1/2，吸收系数不变；序号4是保持衰减系数不变，将吸收系数改变为初始值的3倍；序号5是保持衰减系数不变，将吸收系数改变为初始值的1/3。上述参数的选择是在兼顾皮肤组织正常的光学特性参数范围和模拟结果差异较明显的前提下选择的。

表2 模拟所采用的光学特性参数

图2 模拟流程图

3 结果与讨论

图3为不同光学特性参数生物组织中吸收能量密度的分布情况。散射系数相同，吸收系数越大(μa2>μa1)，吸收能量密度的衰减梯度越大，在浅表层(r<0.07 cm)，吸收能量密度越大。两曲线在r=0.07 cm附近交叉，在深层组织(r>0.07 cm)，吸收能量密度的分布与浅表层正好相反，这是因为浅表层吸收能量密度较大，更多的光子被吸收，因此组织深层吸收能量密度较小。

吸收系数相同，散射系数越大(μs1>μs3),吸收能量密度的衰减梯度越大，光子的散射能力越强，更多的光子被吸收，因此曲线1在浅表层吸收能量密度更大，二者在r=0.04 cm附近交叉。

衰减系数相同，吸收系数越小(μa5<μa1<μa4)，吸收能量密度衰减梯度越小，在浅表层吸收能量密度越小。

图3 光学特性参数对吸收能量密度A(z)的影响

图4为不同光学特性参数生物组织中光能流率的分布。衰减系数越小(μt3<μt1、μt4、μt5<μt2)，光能流率的衰减梯度越小。衰减系数相同(μt1=μt4=μt5)，光能流率的衰减梯度相同，吸收系数越大(μa4>μa1>μa5),光能流率越小。

图4 光学特性参数对光能流率的影响

散射系数相同，吸收系数越小(μa1<μa2)，光能流率越大。吸收系数相同，散射系数不同(μs3<μs1),两曲线在r=0.04 cm附近交叉，这主要是由于光能流率在数值上等于A(z)与μa的比值[4]，二者吸收系数相同，因此光能流率的分布情况与图3中相应曲线类似。

图5 光学特性参数对径向漫反射率的影响

图5为不同光学特性参数的生物组织中漫反射率的分布，漫反射率在径向r迅速衰减。散射系数相同，吸收系数越大(μa2>μa1)，径向漫反射率越小；衰减系数相同，吸收系数越大(μa4>μa1>μa5)，径向漫反射率越小。这主要是由于更多的光子被吸收，导致由于漫反射从上表面离开组织的光子变少。吸收系数相同，散射系数越小(μs3<μs1),发生散射的光子数越少，径向漫反射率越小，其衰减梯度也明显较小。

图6为不同光学特性参数的生物组织中透射率的分布。散射系数相同，吸收系数越小(μa1<μa2)，透射率越大；衰减系数相同，吸收系数越小(μa5<μa1<μa4)，透射率越大，这主要是由于更多光子未被组织吸收，可以通过组织内部的传输，从下表面离开组织。吸收系数相同，散射系数不同(μs3<μs1),二者曲线差异较大。曲线3由于散射系数较小，更多的光子在中心轴附近通过组织内部，从下表面离开组织，因此透射率近轴明显高于曲线1，且从中心轴开始迅速衰减；而曲线1散射系数较大，因此透射率在近轴并未迅速衰减，甚至出现先增大后减小的现象。

图6 光学特性参数对透射率的影响

图7为不同光学特性参数的生物组织中光通量的二维分布图，图中各序号对应表2中各光学特性参数。散射系数相同，吸收系数越小(μa1<μa2)，同一位置光通量越大，但是光子在径向r上分布得更远，这主要是由于吸收系数小的组织中更多的光子被散射。吸收系数相同，散射系数越小(μs3<μs1),光通量的衰减梯度较小，光子在径向上分布得更远。衰减系数相同，吸收系数越小(μa5<μa1<μa4)，散射系数越大，同一位置光通量越大，光通量的衰减梯度越大，光子在径向上分布得更远。

图7 光学特性参数对光通量的影响

4 结语

同一光源辐照单层生物组织，不同的光学特性参数对组织中光分布的各参量影响各不相同。

(1)增大吸收系数和散射系数均会增大吸收能量密度及其梯度。

(2)光能流率和光通量分布类似，吸收系数影响光能流率和光通量的大小，散射系数影响二者的衰减梯度。

(3)漫反射率与透射率分布类似，吸收系数和散射系数均对漫反射率和透射率大小有影响，散射系数影响二者衰减梯度。

[1] 鱼迎珍.纳米金靶向生物组织的光热治疗时域建模仿真[D].西安：西安电子科技大学硕士学位论文,2014.

[2] 黄义梅,龚玮,谢树森.光子嫩肤热损伤的光学相干层析成像动态监控[J].光电子激光,2009,(5):702-704.

[3] 徐小晖.小鼠皮肤深Ⅱ度烫伤及治疗的光学观察与分析[D].福州：福建师范大学硕士学位论文,2012.

[4] 高迪.光在生物组织中传播的蒙特卡罗研究[D].济南：山东大学硕士学位论文,2010.

[5] 肖郑颖,林金忠,李阳娟.蒙特卡罗方法研究皮肤组织光学特性参数对光子传输的影响[J].中国激光医学杂志,2017,(1):7-11.

[6] Jacques S L.Optical properties of biological tissues:a review[J].Physics in Medicine & Biology,2013,58(11):37-61.

【责任编辑牛怀岗】

InfluenceofOpticalParametersontheLightDistributionofBiologicalTissue

XIAO Zheng-ying1, 2
(1. School of Mechanical and Electrical Engineering,Putian University,Putian 351100, China;
2. Fujian Laser Precision Machining Engineering Technology Research Center,Putian 351100, China)

Monte Carlo Method was employed to simulate the distribution of collimated light beam at the biological issues. With a variety of optical parameters including absorption coefficient, scattering coefficient and attenuation coefficient, the differences distribution of energy absorption intensity, fluence rate, diffuse reflection intensity, transmittance and luminous flux were compared and analyzed. The results show that the increasing of absorption coefficient and scattering coefficient shall increase the energy absorption intensity and its gradient. Absorption coefficient shall influence the fluence rate and luminous flux, and scattering coefficient shall influence the attenuation gradients of these two aforementioned parameters. Absorption coefficient and scattering coefficient shall diffuse reflection intensity and transmittance whose attenuation gradients shall be influenced by the scattering coefficient.

tissue optics; Monte Carlo Simulation; optical parameter; light distribution

Q632

：A

：1009-5128(2017)16-0014-06

2017-04-06

福建省激光精密加工工程技术研究中心开放课题：桌面式3D打印机的研发(20150504)；福建省教育厅A类项目：水库枇杷园区精确灌溉控制系统研究(JA11205);莆田市科技计划项目：鞋模的高速铣削加工技术研究(2016G2006)

肖郑颖(1985—)，女，福建莆田人，莆田学院机电工程学院实验师,理学硕士，主要从事生物医学光学研究。