电外科血管闭合组织温度的有限元建模与实验测量研究报告

一、研究概况

本研究由密歇根大学机械工程系及生物医学工程系的Roland K. Chen（学生会员，IEEE）、Matthew W. Chastagner、Robert E. Dodde（会员，IEEE）和Albert J. Shih共同完成。其中，Roland K. Chen为通讯作者。该研究成果发表于《IEEE Transactions on Biomedical Engineering》（IEEE生物医学工程学报）2013年2月第60卷第2期，文献编号为10.1109/TBME.2012.2228265。研究获得了美国国家科学基金会（National Science Foundation）资助项目CMMI 0825795的支持。

二、研究背景与目的

本研究属于生物医学工程领域，具体聚焦于电外科（Electrosurgery）血管闭合过程中的生物传热（Bioheat Transfer）与有限元建模（Finite Element Modeling, FEM）。电外科技术通过向组织输送频率范围为300至550 kHz的电流，利用电阻热效应使组织温度升高，导致蛋白质变性（Denature）从而形成血管闭合。尽管电外科设备相比传统血管结扎具有手术时间短、术后住院时间短等优势，但其在靠近神经血管束（Neurovascular Bundle）等精细部位操作时，热扩散（Thermal Spread）可能导致周围神经的热损伤（Thermal Damage），这一问题在外科手术中备受关注。已有研究显示，神经组织在约43°C时即可发生热损伤，而现有电外科设备造成的热损伤范围往往超过电极边缘3毫米。

为预防热损伤、提高患者安全性，深入理解电外科血管闭合过程中的热扩散机制至关重要。然而，血管闭合过程涉及组织蒸发（Evaporation）导致的相变（Phase Change）及潜热（Latent Heat）吸收、水分流失（Water Loss）引起的组织电导率变化，以及血管内壁接触面的接触电阻（Contact Resistance, CR）与组织融合（Tissue Fusion）等复杂现象，这些因素的耦合使得准确预测组织温度分布极具挑战。此前虽有学者建立了电外科烧灼（Cauterizing）的有限元模型，但尚未有能够综合温度、水含量、热损伤、相变和压缩影响的全面模型。因此，本研究旨在开发并实验验证一个先进的有限元模型，以精确预测电外科血管闭合过程中组织温度的时空分布。

三、研究方法与流程

本研究采用活体猪（In-vivo Porcine）模型进行实验测量，并结合有限元建模进行分析。研究流程主要包括实验装置搭建与数据采集、有限元模型构建与求解、以及模型验证与技术效应分析三个主要部分。

（一）实验装置与数据采集 实验选取一头重45千克的三元杂交猪（50%杜洛克、25%约克夏、25%兰德瑞斯）作为模型，实验方案获得密歇根大学医学院动物管理协议批准。动物经麻醉、插管后，通过呼吸机维持，并持续监测血氧、脉搏和呼吸频率。实验对象为一根直径3毫米的股动脉（Femoral Artery），使用一把5毫米双极切割钳（Gyrus ACMI, Olympus #3005PK）进行血管闭合操作。该设备由一个设定为“蒸汽脉冲（Vapor Pulse, VP）1”模式、功率30瓦的发电机（Gyrus ACMI, Olympus Superpulse Generator #744000）供电。设备引线连接至一台以峰值检测模式运行、采样率10 kHz的数字示波器（Agilent #54833A），通过高压探头和电流探头记录电压与电流数据。

组织温度测量是本实验的关键。研究人员设计了一种独特的测温装置：将一系列微型热敏电阻（Microthermistor，Alpha Technics Model #56A1002-C8）安装在一个聚碳酸酯（Polycarbonate）夹具中，并将该夹具用环氧树脂固定在电极顶部，从而精确控制热敏电阻与电极边缘的相对位置。三个尖端暴露深度为0.75毫米的热敏电阻分别被放置在距电极边缘1.0毫米、2.0毫米和3.0毫米处。这些热敏电阻外径仅0.46毫米，时间常数250毫秒，精度在25°C时为±0.1°C，其小巧的尺寸假设不会干扰电流分布。温度数据由采样率40 Hz的数据采集系统记录，并经过3 Hz低通滤波和窗口宽度为15个数据点的移动平均处理以消除电噪声。

（二）有限元模型构建 有限元模型基于热场和电场两个物理场耦合求解。热场遵循Pennes生物传热方程，忽略血液灌注和代谢热源；电场通过求解拉普拉斯方程获得。两者通过电阻热生成项耦合。模型使用COMSOL Multiphysics 4.2商业软件平台求解。

模型构建的核心是解决血管闭合过程中的三个主要挑战，为此，研究人员开发了三种关键的有限元技术： 1. 技术一：模拟蒸发与比热容变化。 组织约含70%水分，水从液态变为气态所需的潜热是使其从室温升至沸点所需能量的五倍以上，这一巨大吸热效应必须纳入模型。研究人员通过定义一个随温度变化的有效比热容（Effective Specific Heat, c_eff）来实现。c_eff由组织纤维比热容、依赖于温度和水含量的水比热容、以及对应水分流失的潜热贡献三部分加和组成。 2. 技术二：模拟水分流失与电导率变化。 组织温度升高导致水分蒸发，含水量下降，从而使受热损伤的组织导电性显著降低。模型中，当温度低于103°C时，组织电导率随温度以2%/°C的速率增加；一旦组织积累了足够的相变能量，含水量急剧下降，其电导率便被设定为一个极低的恒定值（σ_d = 0.01 S/m），且该过渡过程被平滑处理，以保证数值求解的连续性。 3. 技术三：模拟组织融合与接触电阻。 当双极电极压缩血管时，血管内壁相互接触，在此界面存在不可忽视的接触电阻，它会迫使电流向周围组织扩散，影响温度场。模型中，在血管壁接触处创建了一个厚度仅0.02毫米的薄层作为“接触电阻层”，其初始电导率设为0.005 S/m。当该层温度因电流通过而升高、组织融合过程启动（70-80°C开始）并完成（高于90°C）时，接触电阻层的电导率平滑过渡到与周围受热损伤组织相同的σ_d值，模拟接触电阻的消失。

本模型涉及的物理参数均经过严谨设定。不锈钢电极和血液的材料属性视为常数。血管组织的热导率被设定为温度的函数，并考虑了压缩导致的10%降低。模型几何结构基于实际实验设置进行简化，并利用了两个对称平面（沿血管中心轴的平面A和等分环形电极的中平面B）以减少计算量。计算网格经细化至包含30,022个线性四面体单元，以保证解与网格密度无关。边界条件方面，电学上将实测交流电压转换为均方根值作为准静态直流源施加于活性电极末端，被动电极接地；热学上，对称面设为绝热，远端血管设为36.5°C体温边界，其他外表面设为25°C室温下的自由对流。求解器使用COMSOL的PARDISO直接求解器，采用稳态时间迭代步进方式进行计算，整个计算耗时10至12小时。

四、主要研究结果

有限元模型成功预测了电外科血管闭合过程中组织温度的动态变化趋势，并通过与实验数据的对比得到了验证。根据温度变化特征，闭合过程可分为三个阶段： 1. 第一阶段：融合前期（0-1.5秒）。 此阶段电流主要沿电极边缘的血管壁横向传导，导致该局部区域快速产生热损伤。实验中，1.0毫米处温度缓慢上升，而2.0毫米和3.0毫米处几乎未升温。 2. 第二阶段：融合期（1.5-4.3秒）。 由于接触电阻的存在，电流被迫向更远的区域扩散，导致1.0毫米和2.0毫米处发生显著的电阻加热，温度迅速升高。同时，部分电流穿过接触电阻层，加热血管壁界面，触发并完成主要的组织融合过程。在4.3秒时，1.0毫米处的实验温度达到最大值80.8°C。此阶段结束时间（4.3秒）与文献中形成牢固闭合所需的时间区间（2.06-4.5秒）高度吻合。 3. 第三阶段：融合后期（4.3-15秒）。 电极间组织温度已超过90°C的变性温度，接触电阻层因融合完成而“消失”，电流被约束在两电极片之间的已融合区域内，远离电极的组织不再受到电阻加热。因此，即使发电机持续激活，1.0毫米处温度缓慢下降至74°C，2.0毫米处温度基本稳定，而3.0毫米处温度仅通过热传导缓慢上升。

模型预测与实验测量在三个热电偶位置的平均误差低于7%，最大误差为23.9%。误差主要出现在第一阶段，可能源于热电偶的响应时间延迟以及模型未考虑组织压缩导致的水分挤出效应。尽管如此，模型精确捕捉了各阶段的温升速率和峰值温度时刻，验证了其可靠性。

五、研究结论与价值

本研究提出了一个经过活体实验验证的先进生物传热有限元模型，该模型通过三项关键技术（有效比热容模拟潜热、水分流失对电导率的影响、接触电阻模拟组织融合）的引入，成功地、高精度地预测了电外科血管闭合过程中组织温度的时空分布。对模型技术的单独效应分析进一步证实，若不考虑这些因素，模型将无法预测出真实的温度平台和空间分布。这项研究的科学价值在于，它显著提升了电外科热效应的理论建模水平，揭示了血管闭合过程中“电流路径重分布”与“组织融合”之间的动态耦合机制。其应用价值在于，该模型为理解和定量分析电外科手术中的热扩散风险提供了强大工具，未来可用于辅助开发更安全的电外科设备和控制算法，例如设计带冷却功能的智能电极或优化能量输出时序，从而将热损伤降至最低，提升患者手术安全与术后恢复质量。

六、研究亮点

本研究的亮点主要体现在以下两个方面： 第一，模型方法的独创性与完备性。研究首次系统性地将组织蒸发相变潜热、水分流失导致的电导率剧降、以及接触电阻在组织融合过程中的动态作用这三大关键物理现象集成于一个统一的有限元模型中，显著提升了模型的物理真实性和预测精度。 第二， 实验测量的精巧设计。研究中自制了可精确定位的聚碳酸酯测温夹具，实现了在活体动物模型中对电极周围不同距离组织温度的实时、多点精确测量，为复杂模型的验证提供了可靠且稀缺的体内数据集。