太赫兹时域光谱技术已经用于皮肤癌、 皮肤烧伤以及皮肤疤痕治疗效果的诊断和过程监测, 通常采用太赫兹时域和频域光谱参数作为区分皮肤组织不同状态的诊断参量。 目前最常用的皮肤活体反射式太赫兹测量装置需将皮肤放置于介质窗上来提高测量的准确性, 这使得皮肤表层水分含量因为阻塞而发生变化, 从而影响太赫兹诊断皮肤的精确性。 随着太赫兹生物应用从离体检测到活体测量方向发展, 需要分析阻塞情况下皮肤太赫兹光谱参数的变化规律。 采用太赫兹反射系统对手臂皮肤在石英介质窗上的阻塞过程进行了连续测量, 然后对不同时刻时域波形的峰峰值、 半波脉宽等13个特征量进行了分析。 结果表明: 时域信号和传递函数的峰峰值、 最大值与最小值之间的拟合斜率等参量随阻塞时间而呈指数衰减; 传递函数的最大值与最小值时间差、 最大值与最小值的斜率等随阻塞时间呈指数增加; 时域波形的半峰全宽、 对数频谱等保持不变。 然后用双德拜模型来模拟皮肤组织在0.2～1 THz频段的介电常数, 采用遗传算法和Levenberg-Marquardt相结合的优化算法得到了皮肤阻塞情况下不同时刻的双德拜参数, 结果表明: ε∞和εs均随着时间呈指数上升, 5 min内增幅分别达到了27.8%和12.5%; ε2、 弛豫时间常数τ1和τ2基本保持不变。 将皮肤组织视为多层媒质结构, 基于Bruggeman有效介质模型, 将计算得到的太赫兹反射率和测量得到的反射率作为目标函数, 再采用上述混合优化算法得出了皮肤含水量随时间的变化规律。 结果表明: 随着阻塞时间的增加, 皮肤角质层中的水分以指数函数的形式增加, 5 min的阻塞时间即可增加23.8%。 因此, 在进行临床应用和研究时, 必须谨慎考虑皮肤与介质窗接触阻塞而导致的参数变化。 该研究有利于提高太赫兹活体检测的准确性和推动太赫兹时域光谱技术的临床应用。

针对非线性最小二乘法的传统算法在皮肤组织太赫兹波段德拜模型参数提取时, 对参数初始值的要求很高, 且所得模型与原始数据偏差较大的问题。将遗传算法应用于非线性最小二乘法对原始数据进行拟合, 并提取参数。遗传算法随机产生初始解, 不需要设置参数初始值, 能简化参数提取过程; 其搜索过程是从多个可能解开始, 而不是单一初始值, 所以不易陷入局部最优, 基本能实现全局优化。以皮肤组织德拜参数提取为例, 将遗传算法应用于非线性最小二乘法进行参数提取时不仅减小了文献中的模型在0.15~0.6THz频段与原始数据的偏差, 而且在整个频率范围与原始数据都能较好吻合, 验证了遗传算法提取德拜参数的有效性和皮肤组织太赫兹波段德拜模型的正确性。