生物组织光学性质的测量原理与技术

　第16卷第4期

1997年12月

中　国　生　物　医　学　工　程　学　报

CH I N ESE JOU RNAL O F B I OM ED I CAL EN G I N EER I N G

V o l.16N o.4

D ecem ber1997

生物组织光学性质的测量原理与技术3

谢树森　李　晖

(福建师范大学物理学系,福州350007)

Ch ia T eck Chee

(Schoo l of Science,N anyang T echno logical U niversity,Singapo re1025)本文讨论了组织光学性质参数的测量原理和技术,提出了一种新的测量和计算方法,采用联合测定组织体表面漫反射率和体内光能流率分布,并利用漫射理论和M onte Carlo模型的部分结论,可求出组织的光穿透深度,吸收系数和有效散射系数,以4种猪组织为例,研究了哺乳动物组织的光学性质,这一原理和技术可适用于人体组织光学性质的测量。

关键词:　组织光学;吸收;散射;漫射;M onte Carlo;漫反射率;光能流率

分类号:　R197.39;R318.6

0　前　言

激光医学的进展,尤其是光动力学疗法(PD T)在临床上的深入应用,需要精确了解在一定光照条件下人体组织内的光能分布,以便安排最佳的光治疗方案。其中最关键的问题可归结为如何确定组织体的光学性质基本参数,即吸收系数Λa,散射系数Λs和散射位相函数S(Η)或平均散射余弦g。一旦已知这些光与组织的相互作用参数,在给定的光照方式和边界条件下,光能流率5(r)或其它参量如全反射率R,全透过率T等分布可由有关的数学模型唯一地确定[1,2]。

本文所提出的新方法系采用联合测定组织体表面漫反射率和组织体内部的光能流率分布,并利用漫射理论和M on te Carlo模型的部分结论,可求出组织的光学性质基本参数。

1　组织光学性质参数测量的理论基础

作为电磁波的光在组织中传播行为属于光与组织相互作用问题,在不考虑吸收的情况下,理论上由麦克斯韦方程组及组织体的电磁性质Ε,Λ或折射率,加上边界条件唯一地确定:即在所给定的条件下求解麦克斯韦方程,以得到电矢量在空间中和时间上的分布。其中必然出现一般光学中所有的各种现象,诸如干涉、衍射、反射和偏振等纯粹的物理光学问题。当组织存在光吸收时,应当考虑组织中原子分子的能级结构性质。换言之,此时应采用半经典理论,最严格的处理应使用全量子理论,不难想到,仅由于生物组织折射率的不均匀性,我们就无望获得麦氏方程的数值解,更不用说解析解了。

其实,可以把光在组织体中的传播进而有光能分布的物理实在,用一种粒子的传输过程来

国家自然科学基金和国家教委回国留学人员资助项目

1995年11月27日收稿,1996年4月29日修回

唯象地模拟。粒子的数密度等价为光能。这种假想的粒子无妨称为光子,可以等效于光量子h Τ的集合。同时把组织理解为大量无规则分布的散射粒子和吸收粒子。为此,组织的光学基本参数,可作如下新理解:

Λa :吸收系数,表示组织体中单位长度上一个光子被吸收的几率;

Λs :散射系数,表示组织体中单位长度上一个光子被散射的几率;

S (Η

):散射位相函数,表示散射几率角分布。这里Λa 反映的是组织的原子能级结构性质,而Λs 及S (Η

)则由组织的电磁性质或折射率及其分布决定。这种把光(其本性在经典意义上为电磁波)在组织中的传播看成某种要么被弹性散射要么完全吸收的粒子在组织中传输的方法,称为传输模型理论(T ran spo rt M odel T he 2

o ry )。

该理论中不再出现衍射和偏振等物理光学概念,仅有所谓的可由实验确定的光学性质基本参数。显然,不同波长的光与组织相互作用时,效果不同。这表明,这些参数是光波波长的函数。此外,由于组织体的不均匀性,实验测出的基本参数是一种统计平均值,把组织看成仅有散射和吸收的二种全同粒子,是粗略的考虑。至于更精细的处理,我们准备另文讨论。

有了上述概念后,可以建立Bo ltz m ann 传输方程或M on te Carlo 模型,这是一种结合了实验与数学模型的实际上也是严格的方法。以后,组织光学的所有问题都在此基础上展开研究。111　稳态Bo ltz m ann 传输方程及在一种条件下的近似解

考虑单色的连续光照明下,组织体中的光分布是一种与时间无关的稳定状态。从粒子数守恒出发,可建立辐射亮度的Bo ltz m ann 方程[1,2,4,5]。

s δ L (r ,s δ=-Λt L (r ,s δ)+Λs 4Π∫4ΠP (s δ,s δ’)L (r ,s δ’)d Ξ’(1)

(1)式中:　s 为考察方向的单位矢量;s ’为其它方向的单位矢量;Ξ’为s ’方向的立体角;Λt 为

Λa +Λs ;L (r ,s )为r 处沿s 方向的辐射亮度,单位w c m 2sr ;P (s

δ,s δ’)为沿s δ’入射到r 处向s δ方向散射的几率,该参量满足倒易性,即

P (s δ,s δ’)=P (s δ’,s

δ)且满足归一化条件

14Π∫4ΠP (s δ,s δ’)d Ξ’=1(2) 当P (s

δ,s δ’)沿s 方向为轴对称时,可用相函数分布S (Η)表示S (Η)=∫2Π

0P (s δ,s δ’)sin (Η)d Υ(3)

方程(1)应在一定的边界条件下求解,除了几种简单情形下有解析解外,一般只能有数值

解,具体处理有迭代的方法或其它一些近似的方法[3]。解出L (r ,s

δ),则光能流率分布为5(r )=∫4ΠL (r ,s δ)d Ξ(4)

以下仅考虑宽束平行光垂直入射到半无限大组织体,且散射占主要(Λa νΛs (1-g ))时的解,可把入射光的辐射亮度L 在组织中的传播分成准直(co lli m ated )项L c 和非准直(nonco lli 2m ated )项L d 二部分,即

L (r ,s )=L c (r ,s )+L d (r ,s )

i )准直项L c 满足

823　中国生物医学工程学报16卷