蒙特卡洛建模与模体研究表明全皮肤电子治疗（Total Skin Electron Therapy, TSET）中单位剂量的切伦科夫辐射（Cherenkov emission）发射量是组织光学特性的函数

《MEDICAL PHYSICS》：Monte Carlo modeling and phantom studies show Cherenkov emission per unit dose during total skin electron therapy is a function of tissue optical properties

　　

摘要：背景——切伦科夫成像（Cherenkov imaging）提供了一种无创方法来可视化全皮肤电子治疗（Total Skin Electron Therapy, TSET）中患者体表沉积的辐射剂量。将切伦科夫强度准确转换为辐射剂量对于活体剂量验证（in vivo dosimetry），以评估TSET中的空间剂量分布是必要的。已有研究表明切伦科夫强度与吸收剂量呈线性相关，但组织光学特性对单位剂量切伦科夫发射量的影响尚不清楚。目的——本研究利用蒙特卡洛（Monte Carlo, MC）模拟与实验评估组织光学特性如何影响TSET过程中检测到的单位剂量切伦科夫发射量。方法——采用蒙特卡洛建模模拟TSET期间的切伦科夫产生过程，量化组织光学特性对检测到的切伦科夫发射的影响。研究考察了665 nm波长下临床相关的吸收系数（μa，0.01–1 cm?1）及约化散射系数（μs', 2–40 cm?1）范围，系统评估了组织光学特性、产生深度及切伦科夫发射角分布对组织表面的影响，并将蒙特卡洛结果与一系列固体模体（phantom）测量结果对比，提出了一个拟合组织切伦科夫发射对μa和μs'依赖关系的解析函数。结果——模拟结果表明，切伦科夫发射随组织吸收系数和有效衰减系数增大而减小，随组织散射系数先增后减；80%的表面检出发射来源于皮下0.17–2.0 cm的浅表层；组织中切伦科夫辐射的特定角度产生并未导致表面出现优势出射角，因为大多数切伦科夫光子到达表面前传播方向已被随机化；蒙特卡洛模拟（最大偏差0.23%）与实验相符，标准（最大）百分偏差在2.8%（7.6%）以内。结论——研究人员发现组织光学特性对表面切伦科夫发射有显著影响，单位剂量切伦科夫发射的光学特性依赖性可表达为μa和μs'的二维函数，文中给出了相应的解析表达式；蒙特卡洛模拟与采用一系列组织模拟模体的TSET电子束实验结果吻合。

论文解读：全皮肤电子治疗中切伦科夫单位剂量发射量受组织光学特性调控的蒙特卡洛与模体研究

该研究发表于《Medical Physics》。

一、研究背景与意义

全皮肤电子治疗（Total Skin Electron Therapy, TSET）利用直线加速器产生的电子束治疗全身皮肤疾患，切伦科夫成像（Cherenkov imaging）可通过捕捉高能电子在组织内超光速产生的切伦科夫光子，实时可视化患者体表剂量沉积，有望实现活体剂量验证（in vivo dosimetry）。既往研究发现切伦科夫发射强度与剂量在水等效模体中呈线性关系，但生物组织的吸收系数（absorption coefficient, μ_a）和约化散射系数（reduced scattering coefficient, μ_s' = μ_s(1-g)，其中g为各向异性因子）引起的吸收和散射会改变光子传输，且肤色、色素沉着导致的组织光学特性空间异质性可能使切伦科夫信号差异被误判为剂量不均匀。不同患者及同一患者不同部位的组织光学特性差异大，目前缺乏临床相关范围内将单位剂量切伦科夫发射量表述为μ_a和μ_s'函数的定量关系。因此，研究人员开展本研究，利用蒙特卡洛（Monte Carlo, MC）模拟与固体模体实验，系统阐明组织光学特性对TSET中表面检测切伦科夫发射的影响，并提出解析修正模型。

二、主要关键技术方法

研究人员采用MATLAB编写并经MCmatlab开源代码校验的蒙特卡洛程序，在半无限均匀介质中模拟3.84 MeV电子束（代表临床6 MeV TSET平均能量）产生切伦科夫光子（初始极角44°，方位角均布，波长400–900 nm，初始权重∝λ^?2），追踪光子在给定μ_a（0.01–1 cm^?1@665 nm，基于氧合/脱氧血红蛋白摩尔消光系数计算，StO₂=70%）和μ_s'（2–40 cm^?1@665 nm，μ_s'(λ)=Aλ^?b，b=0.838）下的多重散射与权重衰减，记录逃逸光子随出射角、半径、波长及产生深度z₀的分布，并利用互易定理推算宽束电子入射结果。辅以含100 μm含黑色素表皮层的双层模型评估黑色素影响（通过漫反射光谱Diffuse Reflectance Spectroscopy, DRS拟合获取等效单层μ_a、μ_s'）。实验方面，制作碳-二氧化钛硅胶固体模体（九组μ_a、μ_s'经双导管法标定），在SSD=500 cm、6 MeV HDTSE条件下照射，用门控增强CMOS相机采集切伦科夫图像，线性拟合强度–MU（Monitor Unit）斜率归一化得到实测单位剂量切伦科夫发射，与模拟及解析公式比对。

三、研究结果

3 RESULTS——切伦科夫发射谱与总强度对光学特性的依赖

模拟谱显示尽管切伦科夫产生于紫外/蓝光区，但表面检出光因血红蛋白强吸收主要位于"生物透射窗口"600–800 nm（峰值~700 nm）。对400–900 nm积分得总切伦科夫强度I_Che,total，归一化至患者皮肤均值（μ_a,665nm=0.3 cm^?1，μ_s',_665nm=20 cm^?1）。结果显示：I_Che,total随μ_a及有效衰减系数μ_eff=√[3μ_a(μ_a+μ_s')]增大单调下降；随μ_s'先升后降（存在最优散射范围）；与漫反射率R_d及μ_s'呈二维依赖关系。据此拟合得到未归一化切伦科夫强度关于μ_a和μ_s'的二维解析表达式（原文式3），与MC结果标准（最大）偏差为0.6%（1.8%）。

3 RESULTS——切伦科夫发射的角分布

按出射角θ积分得到I_Che,θ，发现各光学特性下出射角分布近似符合Lambertian分布（∝sin2θ=2sinθcosθ），在θ<60°内I_Che,θ/sin2θ波动于±5%，θ>60°明显偏离——表明虽切伦科夫光子初生于特定切伦科夫角（ζ≈44°），但因在典型生物组织中穿越距离远超自由程使其传播方向被完全随机化，不保留原始方向性。

3 RESULTS——模体实验验证

以碳吸收谱的MC模拟及基于phantom光学特性拟合的解析公式（原文式4，标准/最大偏差0.6%/1.0%）与九组硅胶模体实测对比，In-house MC与实验标准（最大）百分偏差为2.8%（7.6%），三者趋势一致，证实光学特性依赖关系的正确性；同时指出若仅μ_a、μ_s'@665 nm相同但吸收谱不同（血红蛋白vs碳），总切伦科夫发射将有显著差异，原式（式3）仅适用于具血红蛋白特征的生物组织。

3 RESULTS——切伦科夫光子的产生深度分布

追踪初生深度z₀（均匀分布0–2 cm）发现表面检出光子数随z₀指数衰减。对各光学条件计算贡献80%表面切伦科夫发射的浅表层厚度，结果为0.17–2.0 cm（均值0.72 cm）：μ_a越大、μ_s'越大，该厚度越薄；低μ_a高μ_s'时可延至模拟上限2.0 cm，说明大部分可用信号源自浅层剂量沉积区。

3 RESULTS——表层黑色素影响的双层模型评估

对fv=0.02（浅色皮肤）、0.15（晒黑皮肤）、0.30（深色皮肤）的100 μm黑色素表皮+底层组织（μ_a=0.1 cm^?1, μ_s'=10 cm^?1），双层MC与通过DRS谱拟合得到的等效单层 bulk光学特性（μ_a略升高）代入单层模型所得总切伦科夫发射偏差分别为?4.3%、?3.9%、+1.3%，证明用单层和患者个体化bulk μ_a、μ_s'可初步校正表层色素影响。

四、讨论与结论翻译

讨论指出：本工作明确了TSET中单位剂量切伦科夫发射是μ_a和μ_s'的二维函数而非仅依赖μ_eff或R_d，给出适用生物组织光谱的解析形式；检出光子多来自皮下<2 cm且经多次散射方向随机化故呈Lambertian出射（θ<60°）；表层黑色素可用等效单层bulk属性近似；强调吸收光谱形状须匹配组织（血红蛋白）方可套用该公式。此框架结合患者皮肤光学特性测量可校正切伦科夫–剂量转换因子空间异质性，提升TSET活体剂量验证精度。

结论翻译：研究人员发现组织光学特性显著影响TSET期间表面切伦科夫发射，单位剂量切伦科夫发射的光学特性依赖性可表达为μ_a,665nm和μ_s',_665nm的二维函数，文中给出了相应解析表达式。切伦科夫发射主要来自浅表组织层——在所研究光学特性范围内，80%的表面切伦科夫发射源于深度0.17–2.0 cm的浅表组织层，平均厚度0.72 cm。当切伦科夫光子出射方向与组织–空气界面法线夹角小于60°时服从Lambertian分布，大于60°时迅速偏离理论分布。组织中切伦科夫辐射的特定角度产生并未导致表面出现优势出射角，因为在到达表面前多数切伦科夫光子的传播方向已被随机化。蒙特卡洛模拟结果与采用一系列组织模拟模体的TSET电子束实验相吻合。