太赫兹辐射在皮肤中的传输仿真与安全性分析

太赫兹辐射在皮肤中的传输仿真与安全性分析 Transmission simulation and safety analysis of terahertz radiation in skin first-author 刘欣缘（1996-），男，江苏南京人，博士研究生，2019年于武汉大学获得学士学位，主要从事太赫兹辐射传输仿真方面的研究。E-mail： lxy1996@whu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850502&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850505&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850504&type= 22.01333618 32.00399780 刘欣缘 （1996-），男，江苏南京人，博士研究生，2019年于武汉大学获得学士学位，主要从事太赫兹辐射传输仿真方面的研究。E-mail： lxy1996@whu.edu.cn lxy1996@whu.edu.cn corresp E-mail： ls@whu.edu.cn E-mail： ls@whu.edu.cn ls@whu.edu.cn 武汉大学 电子信息学院，湖北 武汉 430079 武汉大学 电子信息学院，湖北 武汉 430079 School of Electronic Information， Wuhan University， Wuhan 430079， China School of Electronic Information， Wuhan University， Wuhan 430079， China 太赫兹辐射中的低频段由于具有很好的穿透性，在太赫兹安检成像系统中得到了广泛的应用，但这增加了人体暴露在太赫兹辐射下的风险，也引发了对太赫兹辐射安全性的担忧。本文建立了皮肤的分层结构模型，计算了皮肤中各层组织的复介电常数和各种光学参数，然后利用蒙特卡洛仿真计算了太赫兹光子在皮肤组织中的传输过程，定量地给出了皮肤各层组织对不同频率太赫兹辐射的衰减倍数。针对穿透性最强的0.1 THz辐射，对光子的穿透深度进行了统计，发现99.9%的光子穿透深度小于0.55 mm，结果表明即使是功率高于损伤阈值的太赫兹辐射也很难传输到皮肤的真皮层下部，证明太赫兹辐射具有较高的安全性，不会对体内器官造成影响。 Low-frequency terahertz radiation has been widely employed in terahertz security imaging systems owing to its excellent penetrability. However， this procedure also increases the risk of human exposure to terahertz radiations and raises safety concerns. In this study， the layered structure model of the skin was first established， following which the complex dielectric constant and various optical parameters of skin tissues were calculated. Next， the process of terahertz photon transmission in the skin tissues was determined via Monte Carlo simulation. Finally， the attenuation of the skin tissues to terahertz radiation at different frequencies was determined quantitatively. The penetration depth of 0.1 terahertz photons was analyzed， and it was found that 99.9% of the penetration depth was within 0.55 mm. These results show that it is difficult to transmit terahertz radiation to the lower dermis of the skin， even when the radiation power is higher than the damage threshold； this proves that terahertz radiation is safe and will not affect internal organs. 光电检测 太赫兹辐射 人体安全性 蒙特卡洛仿真 optoelectronic detection terahertz radiation human safety Monte Carlo simulation 1 引　言 1 太赫兹辐射（0.1~10 THz）介于红外与毫米波之间，具有通信带宽宽、分子特征光谱明显和对非极性物质穿透性强等特点，因此在高速无线通信、体外疾病诊断、非接触式人体安检等领域表现出了广泛的应用前景 ［ 1 1 - 6 6 1-6 ］ 。特别是低频太赫兹（ < 1 THz）常用于安检成像系统，目前被动成像安检仪的应用已经比较成熟，但存在成像分辨率有限，受环境温度影响大的问题。随着太赫兹辐射源输出功率和太赫兹探测器灵敏度的不断提升，分辨率更高、信噪比更优、受环境影响小的主动成像式太赫兹安检系统的发展迅速 ［ 7 7 - 9 9 7-9 ］ 。主动式太赫兹成像系统的普及无疑会增加人体暴露在太赫兹辐射之下的风险，也引发了对太赫兹生物安全性的担忧。已有研究表明，高强度的太赫兹辐射能在多种不同尺度上引起生物材料的变化，常见的效应包括结构蛋白损伤、细胞器功能破坏、细胞死亡和组织凝固等。此外，由于生物组织富含的水分子能有效地吸收太赫兹光子的能量并转换为热能，因此太赫兹辐射在生物组织中的热效应也十分显著 ［ 10 10 ］ 。 正是由于生物组织对太赫兹辐射的强吸收效应，一般认为太赫兹辐射在穿透衣物后会被人体表层皮肤细胞吸收，而不会对人身安全造成威胁 ［ 7 7 ］ 。但这只是一种定性的分析，实际上目前对太赫兹辐射在皮肤中的传输过程与吸收特性还缺乏清晰的认识。Vilagosh等人使用时域三维电磁仿真软件对太赫兹信号在皮肤中的传播进行了仿真 ［ 11 11 ］ ，仿真结果表明毛囊能增加太赫兹辐射在皮肤中的传播距离，但频率为0.45 THz的信号在皮肤中的穿透深度仍小于0.25 mm。Betzalel等人使用CST软件对频率为0.4~0.5 THz的信号在皮肤中的电场分布进行仿真 ［ 12 12 ］ ，仿真结果表明皮肤中汗腺区域的电场强度较高，有助于促进太赫兹信号在皮肤中的传导。相较于从电磁场角度进行模拟的有限时域差分算法，蒙特卡洛方法从粒子角度对光子在介质中的传播进行数值模拟，具有简单、准确的优点 ［ 13 13 - 15 15 13-15 ］ 。本文建立了皮肤分层结构模型，计算了皮肤中各层组织的复介电常数和光学参数，并通过蒙特卡洛仿真模拟太赫兹辐射在皮肤中的传输过程，进而定量计算了太赫兹辐射在皮肤中由吸收与散射引起的衰减。 2 皮肤组织模型 1 2.1　皮肤分层结构 2 人体皮肤组织由表皮、真皮、皮下组织及皮肤附属器构成。综合考虑皮肤的组织结构与电磁特性，在仿真中将皮肤分为4层，如 图1 图1 所示。第一层是去角质的扁平死细胞，厚度约为20 μm，称为角质层。角质层下的棘层与基底层在结构与组成成分上存在比较大的差异，因此又分为两层，其中棘层厚度为30 μm，基底层厚度为50 μm，主要包含脱水细胞、柱状细胞等活细胞以及角质透明蛋白颗粒和黑色素颗粒等 ［ 12 12 ］ 。第四层为真皮层，由于皮肤内血管和毛细血管的分布不均匀，真皮层又可细分为乳头状真皮层（厚度约为150 μm），浅表血管层（厚度约为80 μm）、网状真皮层（厚度约为1 500 μm）和深部血管层（厚度约为170 μm），但真皮层整体的折射率、吸收系数等光学参数比较一致，因此在仿真中作为一层来考虑。真皮层下即为皮下脂肪层，这里不作为皮肤的一部分来考虑。 2.2　皮肤参数计算 2 皮肤组织的介电性能由其复介电常数来表示。复介电常数与组织中传输的电磁波频率相关，对于频率在0.1~2 THz的太赫兹辐射，复介电常数的值可以通过双德拜方程来计算，其中复介电常数的实部和虚部分别表示为： <math display="block" id="M1"><msup><mrow><mi>ε</mi></mrow><mrow><mi>'</mi></mrow></msup><mo>=</mo><msub><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub><mo>+</mo><mfrac><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850516&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850514&type= 55.11800003 10.32933331 ， <math display="block" id="M2"><msup><mrow><mi>ε</mi></mrow><mrow><mo>″</mo></mrow></msup><mo>=</mo><mfrac><mrow><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><mfenced separators="|"><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ε</mi></mrow><mrow><mi mathvariant="normal">∞</mi></mrow></msub></mrow></mfenced><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mi>ω</mi><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850520&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850518&type= 58.16600037 11.51466751 ， 其中： ε ∞ 为高频率极限的介电常数， ε 1 为静态介电常数， ε 2 为中间频率极限值， τ 1 和 τ 2 是松弛时间，下文中将以上参数统称为双德拜系数； ω 为角频率。 由于皮肤组织中各层的双德拜系数并不是已知的，本文根据相关论文中给出的太赫兹离散频率的复介电常数和皮肤整体的双德拜系数进行拟合 ［ 11 11 ， 16 16 ］ ，得到角质层、棘层、基底层和真皮层的双德拜系数，如 表1 表1 所示。 Layers ε ∞ ε 1 ε 2 τ 1 /ps τ 2 /ps Stratum corneum 2.71 17.31 3.56 4.89 0.33 Stratum spinosum 2.73 24.59 3.97 2.65 0.29 Stratum basale 2.64 27.23 4.41 1.84 0.06 Dermis 2.78 27.17 3.67 2.05 0.16 根据表中的系数，利用 式（1） 式（1） 和 式（2） 式（2） 可以计算得到皮肤各层在0.1~2 THz内的复介电常数。 图2 图2 中展示了根据拟合的双德拜系数计算的复介电常数与原始数据的对比，从图中可以看出，无论是实部还是虚部，通过离散点拟合得到的双德拜系数所计算得到的复介电常数与原始数据较为吻合，仅在较低频率范围存在误差。将本文中的拟合结果与Truong ［ 17 17 ］ 和E.Pickwell ［ 18 18 ］ 通过实验测量的双德拜系数计算的复介电常数进行对比，如 图3 图3 ~ 图4 图4 所示。可以看到，由于计算中将皮肤分为4层，各层的成分不同，而实验中则是将皮肤整体当作一层均匀介质，且Truong取用的是肿瘤附近的皮肤样本，而E.Pickwell采用的是手腕处的皮肤样本，因此介电常数在数值上略有区别，但基本相当。特别是皮肤角质层的计算结果的实部与虚部都接近E.Pickwell的实验结果，因此使用这种方法得出的皮肤复介电常数的结果是比较合理的。 通过比较皮肤角质层、棘层、基底层和真皮层的复介电常数可知，无论是实部还是虚部，皮肤各层的复介电常数的数值均随着频率的升高而逐渐下降，其中在0.1~0.4 THz频段内下降速度较快；而0.4~2 THz频段内下降速度趋于平缓。在皮肤各层中，角质层的复介电常数最小，而基底层的复介电常数最大，真皮层的复介电常数略大于棘层的复介电常数。 根据复介电常数可以进一步计算皮肤的折射率、吸收系数、散射系数与散射各向异性因子 ［ 19 19 - 21 21 19-21 ］ ，结果如 图5 图5 ~ 图8 图8 所示。其中，皮肤各层的散射各向异性因子均在10 - 5 量级，说明太赫兹辐射在皮肤中的散射可以看作是各向同性的。 3 蒙特卡洛仿真 1 3.1　仿真原理 2 蒙特卡洛法是一种以概率和统计理论方法为基础的计算方法，使用随机数或更常见的伪随机数来解决很多复杂的计算问题，将所求解的问题同一定的概率模型相联系，用计算机实现统计模拟或抽样，以获得问题的近似解。在某些场合中，传统的经验方法不能逼近真实的物理过程，很难得到满意的结果，而蒙特卡洛法则能够真实地模拟实际物理过程，得到比较准确的结果。而太赫兹辐射在生物组织中的传播就是一个适合用蒙特卡洛法模拟的物理过程 ［ 22 22 - 24 24 22-24 ］ 。 由于生物组织内部散射粒子分布的不确定性，光子在组织中的传播是一个随机迁移的过程，并遵循一定的统计规律。因此，本文采用蒙特卡洛仿真模拟太赫兹光子在介质中的传输。该方法描述了光子传播的局部规则，在最简单的情况下，用概率分布来描述光子在介质体积元素间移动的步长，以及在发生散射事件时光子轨道的偏转角度。三维生物组织模型可分解为若干个体积相同的微小体积元素，每个体积元素的光学参数可以分别设置。当光子在生物组织模型中模拟传播时，需要考虑每一次步长的位移是否离开光子当前所在的体积元素，若离开了当前体积元素进入另一个体积元素，则需要比较两个体积元素的光学参数，考虑路径传播过程中的折射、散射和吸收等相互作用。因此，蒙特卡洛仿真主要考虑生物组织的4个光学参数：折射率、吸收系数、散射系数和散射各向异性因子，根据这些光学参数可以定量地描述光子的传播过程。 3.2　仿真流程 2 蒙特卡洛仿真流程如 图9 图9 所示，该流程描述了光子从产生到消亡的过程，主要可以分为5个步骤： （1）初始化光子的权重 w ，根据光源类型和方向确定光子的初始传播方向； （2）根据当前所在体积元素的吸收系数 μ abs 和散射系数 μ sca ，并利用随机数随机产生一个位移步长 s ： <math display="block" id="M3"><mi>s</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">g</mi><mi>N</mi></mrow><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">b</mi><mi mathvariant="normal">s</mi></mrow></msub><mo>+</mo><msub><mrow><mi>μ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">a</mi></mrow></msub></mrow></mfrac><mo>,</mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mi>N</mi><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mo> </mo><mo> </mo><mn mathvariant="normal">1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850570&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850568&type= 43.01066971 8.46666718 ， 式中 N 为在［0，1］区间内均匀分布的随机数。 （3）根据传播方向计算当前位置到体积元素边界的距离 d ，并与位移步长 s 比较；若 d > s ，则表明该位移仅在当前的体积元素内进行，因此光子按照该传播方向与传播距离移动，根据当前体元的吸收系数 μ abs 更新光子的权重 w = w· exp（ -μ abs · Δ s ）；若 d < s ，则表明该位移将穿越当前的体积元素，因此需要先将光子传播至体元边界，根据当前体元的吸收系数 μ abs 更新光子的权重 w ，将步长 s 减去位移，并根据相邻体元的光学参数计算新的传播方向； （4）判断光子的权重 w 是否小于阈值 w th ，若小于阈值则代表光子濒临消亡，若大于阈值则返回步骤（2）； （5）生存盘赌是为了避免计算过长的传播路径并且对计算结果的影响较小的情况。当光子权重 w 下降到阈值以下时，会以1 -c 的概率终止（通常 c < 0.1），或以 c 的概率提升到 w / c ，因此，所有光子总权重得到很好的保存。 3.3　仿真结果与分析 2 仿真中模拟的皮肤组织结构如 图10 图10 所示（彩图见期刊电子版），使用笛卡尔坐标系进行表示， X - Y 平面与皮肤表面重合， Z 轴垂直于皮肤表面向下。皮肤各层为平行分层结构，各层的厚度依据皮肤组成模型，图中不同的层组织使用不同的颜色表示，从上至下依次为角质层、棘层、基底层和真皮层。图中的黄色平行直线表示太赫兹辐射垂直入射进入人体皮肤组织，范围为半径50 μm的圆形区域。皮肤表面的入射辐照度设定为7.16 W/cm 2 ，这是美国空军研究实验室测得的皮肤损伤阈值，皮肤在这种高强度太赫兹辐射下照射2 s会升温31.72 ℃，并产生不可逆损伤 ［ 25 25 ］ 。 利用上述蒙特卡洛法对0.1~2 THz之间的太赫兹辐射在皮肤中的传播过程进行仿真，结果如 图11 图11 所示。图中绘出了光斑范围内的太赫兹辐射在依次通过皮肤的角质层、棘层与基底层后的平均能量衰减，衰减的能量部分被吸收，部分由于折射和散射的作用改变了传播方向，离开了光斑的中心范围，不足以造成损伤。从图中可以看到，0.1 THz的太赫兹辐射的穿透性最强，三层总衰减只有 - 7.3 dB，而2 THz的太赫兹辐射则衰减了 - 41.2 dB。因此对于高频太赫兹来说，表面三层总厚度为100 μm的皮肤组织已经足以吸收大部分的能量，而对于低频段，能量衰减主要是在真皮层中完成的。仿真结果表明，只有0.1 THz的少部分光子能够进入到皮下脂肪中，如 图12 图12 所示，而0.2 THz以上的光子都无法穿过真皮层继续传播。 从 图12 图12 中可以看出，除了强烈的吸收，太赫兹辐射在皮肤中的散射效应也非常明显，大量的光子经过散射改变了传播方向，降低了局部区域的能量密度，也起到了一定的保护作用。在吸收和散射的共同作用下，即使是穿透能力最强的0.1 THz也很难穿过真皮层继续传播，如 图13 图13 所示，皮下脂肪中的光子数比角质层中相差约10个量级，可以忽略不计。实际上，99.9%的光子的穿透深度都小于0.55 mm，所以即使是很强的太赫兹辐射也很难对真皮层底部造成实质性的伤害。 4 结　论 1 本文建立了皮肤组织的分层模型，通过拟合双德拜系数计算了各层组织在0.1~2 THz之间的复介电常数和各项光学参数，然后使用蒙特卡洛法计算了太赫兹辐射在人体皮肤组织中的传输过程。仿真结果表明，皮肤组织在对太赫兹辐射有强烈吸收的同时，对组织中的细胞和蛋白质等粒子的散射效应也比较显著，在两者的共同作用下，太赫兹辐射经过皮肤表层的角质层、棘层和基底层时均有明显的衰减，且衰减倍数随着频率的增加而增加。频率较高的太赫兹辐射都会在厚度达到mm量级的真皮层完全消失，只有穿透性最强的0.1 THz能有少量光子穿过真皮层达到皮下脂肪，其数量也可以忽略不计。因此，皮肤确实能够保护人体内部免受太赫兹辐射的影响，但由于角质层的衰减作用有限，棘层、基底层和真皮层上部的活性细胞暴露在太赫兹辐射下，仍有引发皮肤病的可能。 参考文献 姚建铨 ， 迟楠 ， 杨鹏飞 ， 等 . 太赫兹通信技术的研究与展望 ［J］. 中国激光 ， 2009 ， 36 （ 9 ）： 2213 - 2233 . 太赫兹通信技术的研究与展望 YAO J Q ， CHI N ， YANG P F ， et al . Study and outlook of terahertz communication technology ［J］. Chinese Journal of Lasers ， 2009 ， 36 （ 9 ）： 2213 - 2233 . （in Chinese） Study and outlook of terahertz communication technology 谢莎 ， 李浩然 ， 李玲香 ， 等 . 太赫兹通信技术综述 ［J］. 通信学报 ， 2020 ， 41 （ 5 ）： 168 - 186 . 太赫兹通信技术综述 XIE SH ， LI H R ， LI L X ， et al . Survey of terahertz communication technology ［J］. Journal on Communications ， 2020 ， 41 （ 5 ）： 168 - 186 . （in Chinese） Survey of terahertz communication technology 孙美玉 ， 祁峰 ， 雷静 ， 等 . 太赫兹技术在医疗领域的研究现状与应用前景 ［J］. 中国医疗设备 ， 2018 ， 33 （ 7 ）： 1 - 8，36 . 太赫兹技术在医疗领域的研究现状与应用前景 SUN M Y ， QI F ， LEI J ， et al . Research status and application prospect of terahertz technology in medical field ［J］. China Medical Devices ， 2018 ， 33 （ 7 ）： 1 - 8，36 . （in Chinese） Research status and application prospect of terahertz technology in medical field 李晓媛 ， 武鹏 ， 刘允 ， 等 . 基于人脸视频的心率参数提取 ［J］. 光学 精密工程 ， 2020 ， 28 （ 3 ）： 548 - 557 . 基于人脸视频的心率参数提取 LI X Y ， WU P ， LIU Y ， et al . Extraction of heart rate parameters from video of human face ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 3 ）： 548 - 557 . （in Chinese） Extraction of heart rate parameters from video of human face 叶麾 ， 郄明蓉 ， 曹寒雨 ， 等 . 太赫兹技术在医学科学中的应用及研究进展 ［J］. 光电工程 ， 2018 ， 45 （ 5 ）： 170528 . 太赫兹技术在医学科学中的应用及研究进展 YE H ， XI M R ， CAO H Y ， et al . Applications of terahertz technology in medical science and research progress ［J］. Opto-Electronic Engineering ， 2018 ， 45 （ 5 ）： 170528 . （in Chinese） Applications of terahertz technology in medical science and research progress 蒋林华 ， 王尉苏 ， 童慧鑫 ， 等 . 太赫兹成像技术在人体安检领域的研究进展 ［J］. 上海理工大学学报 ， 2019 ， 41 （ 1 ）： 46 - 51 . 太赫兹成像技术在人体安检领域的研究进展 JIANG L H ， WANG W S ， TONG H X ， et al . Research progress of terahertz imaging in the field of human security ［J］. Journal of University of Shanghai for Science and Technology ， 2019 ， 41 （ 1 ）： 46 - 51 . （in Chinese） Research progress of terahertz imaging in the field of human security 刘杰 . 基于稀疏线阵MIMO-SAR的 0 .14 THz安检成像系统［D］. 北京 ： 中国工程物理研究院 ， 2018. 基于稀疏线阵MIMO-SAR的 LIU J . A sparse MIMO-SAR-based 0 .14 THz Security Inspection Imaging System［D］. Beijing ： China Academy of Engineering Physics ， 2018. （in Chinese） A sparse MIMO-SAR-based 王与烨 ， 陈霖宇 ， 徐德刚 ， 等 . 太赫兹波三维成像技术研究进展 ［J］. 中国光学 ， 2019 ， 12 （ 1 ）： 1 - 18 . 太赫兹波三维成像技术研究进展 WANG Y Y ， CHEN L Y ， XU D G ， et al . Advances in terahertz three-dimensional imaging techniques ［J］. Chinese Optics ， 2019 ， 12 （ 1 ）： 1 - 18 . （in Chinese） Advances in terahertz three-dimensional imaging techniques 周德亮 . 太赫兹主动成像系统研究 ［D］. 南京 ： 南京大学 ， 2018 . 太赫兹主动成像系统研究 ZHOU D L . Research on Active Terahertz Imaging Systems ［D］. Nanjing ： Nanjing University ， 2018 . （in Chinese） Research on Active Terahertz Imaging Systems WILMINK G J ， GRUNDT J E . Invited review article： current state of research on biological effects of terahertz radiation ［J］. Journal of Infrared ， Millimeter， and Terahertz Waves， 2011 ， 32 （ 10 ）： 1074 - 1122 . Invited review article： current state of research on biological effects of terahertz radiation VILAGOSH Z ， LAJEVARDIPOUR A ， WOOD A W . Modelling terahertz radiation absorption and reflection with computational phantoms of skin and associated appendages ［C］. Nanophotonics Australasia 2017. December 10 - 13 ， 2017 . Melbourne， Australia. SPIE ， 2018 ： 104560 M. Modelling terahertz radiation absorption and reflection with computational phantoms of skin and associated appendages BETZALEL N ， FELDMAN Y ， ISHAI P B . The modeling of the absorbance of sub-THz radiation by human skin ［J］. IEEE Transactions on Terahertz Science and Technology ， 2017 ， 7 （ 5 ）： 521 - 528 . The modeling of the absorbance of sub-THz radiation by human skin IBEY B L ， PAYNE J A ， MIXON D G ， et al . Comparing finite difference time domain and Monte Carlo modeling of human skin interaction with terahertz radiation ［C］. SPIE BiOS. Proc SPIE 6854 ， Optical Interactions With Tissue and Cells XIX， San Jose， California， USA . 2008 ， 6854 ： 68540Z . Comparing finite difference time domain and Monte Carlo modeling of human skin interaction with terahertz radiation 李辉 ， 刘巍 ， 张洋 ， 等 . 激光跟踪仪多基站转站精度模型与误差补偿 ［J］. 光学 精密工程 ， 2019 ， 27 （ 4 ）： 771 - 783 . 激光跟踪仪多基站转站精度模型与误差补偿 LI H ， LIU W ， ZHANG Y ， et al . Model establishment and error compensation of laser tracker station-transfer ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 4 ）： 771 - 783 . （in Chinese） Model establishment and error compensation of laser tracker station-transfer 万振华 ， 赵开春 ， 褚金奎 . 基于偏振成像的方位测量误差建模与分析 ［J］. 光学 精密工程 ， 2019 ， 27 （ 8 ）： 1688 - 1696 . 基于偏振成像的方位测量误差建模与分析 WAN ZH H ， ZHAO K CH ， CHU J K . Modeling and analysis of orientation measurement error based on polarization imaging ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 8 ）： 1688 - 1696 . （in Chinese） Modeling and analysis of orientation measurement error based on polarization imaging REID C . Spectroscopic Methods for Medical Diagnosis at Terahertz Wavelengths ［D］. London ： University College London ， 2009 . Spectroscopic Methods for Medical Diagnosis at Terahertz Wavelengths TRUONG B C Q ， TUAN H D ， WALLACE V P ， et al . The potential of the double Debye parameters to discriminate between basal cell carcinoma and normal skin ［J］. IEEE Transactions on Terahertz Science and Technology ， 2015 ， 5 （ 6 ）： 990 - 998 . The potential of the double Debye parameters to discriminate between basal cell carcinoma and normal skin PICKWELL E ， COLE B E ， FITZGERALD A J ， et al . Simulation of terahertz pulse propagation in biological systems ［J］. Applied Physics Letters ， 2004 ， 84 （ 12 ）： 2190 - 2192 . Simulation of terahertz pulse propagation in biological systems BHANDARI A ， HAMRE B ， FRETTE Ø ， et al . Modeling optical properties of human skin using Mie theory for particles with different size distributions and refractive indices ［J］. Optics Express ， 2011 ， 19 （ 15 ）： 14549 - 14567 . Modeling optical properties of human skin using Mie theory for particles with different size distributions and refractive indices WISCOMBE W J . Improved Mie scattering algorithms ［J］. Applied Optics ， 1980 ， 19 （ 9 ）： 1505 - 1509 . Improved Mie scattering algorithms LI N ， ZHU Y Q ， WANG Z H . A discussion on the applicable condition of Rayleigh scattering ［J］. International Journal of Remote Sensing Applications ， 2015 ， 5 ： 62 . A discussion on the applicable condition of Rayleigh scattering WANG L H ， JACQUES S L ， ZHENG L Q . MCML-Monte Carlo modeling of light transport in multi-layered tissues ［J］. Computer Methods and Programs in Biomedicine ， 1995 ， 47 （ 2 ）： 131 - 146 . MCML-Monte Carlo modeling of light transport in multi-layered tissues ZHANG Y ， CHEN B ， LI D ， et al . Efficient and accurate simulation of light propagation in bio-tissues using the three-dimensional geometric Monte Carlo method ［J］. Numerical Heat Transfer ， Part A： Applications， 2015 ， 68 （ 8 ）： 827 - 846 . Efficient and accurate simulation of light propagation in bio-tissues using the three-dimensional geometric Monte Carlo method LEINO A A ， PULKKINEN A ， TARVAINEN T . ValoMC： a Monte Carlo software and MATLAB toolbox for simulating light transport in biological tissue ［J］. OSA Continuum ， 2019 ， 2 （ 3 ）： 957 - 972 . ValoMC： a Monte Carlo software and MATLAB toolbox for simulating light transport in biological tissue DALZELL D R ， MCQUADE J ， VINCELETTE R ， et al . Damage thresholds for terahertz radiation ［C］. Conference on optical interactions with tissues and cells XXI . Damage thresholds for terahertz radiation Human Effectiveness Directorate ， Air Force Research Laboratory ， Brooks City-Base ， TX；rnHuman Effectiveness Directorate ， Air Force Research Laboratory ， Brooks City-Base ， TX；rnHuman Effectiveness Directorate ， Air Force Research Laboratory ， Brooks City-Base ， 2010 ： 75620 M. 10.37188/OPE.20212905.0999 TB133 A 1004-924X（2021）05-0999-09 刘欣缘,曾昊旻,田昕等.太赫兹辐射在皮肤中的传输仿真与安全性分析[J].光学精密工程,2021,29(05):999-1007. 刘欣缘,曾昊旻,田昕等.太赫兹辐射在皮肤中的传输仿真与安全性分析[J].光学精密工程,2021,29(05):999-1007. LIU Xin-yuan,ZENG Hao-min,TIAN Xin,et al.Transmission simulation and safety analysis of terahertz radiation in skin[J].Optics and Precision Engineering,2021,29(05):999-1007. LIU Xin-yuan,ZENG Hao-min,TIAN Xin,et al.Transmission simulation and safety analysis of terahertz radiation in skin[J].Optics and Precision Engineering,2021,29(05):999-1007. 国防科技创新特区项目 2020-09-18 2020-11-17 2021-05-15 2021-06-24 光学精密工程 5 29