基于EIT测定不同体姿下通气分布实验方法研究

方海东，骆栩潇，吴阳，周童，高贵锋，刘凯，姚佳烽

(1. 南京航空航天大学 机电学院，江苏 南京 210016；2. 苏州大学附属第二医院 呼吸与危重症医学科，江苏 苏州 215004；3. 深圳市安保医疗科技股份有限公司，广东 深圳 518100)

0 引言

电阻抗成像(electrical impedance tomography, EIT)是一种功能性成像技术，通过向外部电极传感阵列施加一定频率的安全电流，并测量边界电压变化，以此重构生物内部阻抗分布图像，实现诊断或监测患者的功能[1]。EIT低成本、实时动态功能成像和对组织功能变化高度敏感等优点，使其非常适用于床边连续监测肺功能[2-3]。在医学领域中，EIT也被积极研究用于乳腺癌检测[4-5]、体外血栓检测[6]和脑检测[7]等。

体姿变化会导致重力依赖区的通气改变，从而导致肺内气体的重新分布以及造成重力依赖区肺排空能力的改变[8]。EIT可以实时地监测肺局部通气量的变化，因此，获得了国内外研究小组的关注。FRERICHS I等[9]通过EIT确定了不规则自主呼吸模式和体姿对无呼吸疾病新生儿通气空间分布的影响。FRERICHS I小组[10]通过对青年和老人在不同体姿下的EIT研究，确定了年龄和体位对局部通气的影响。张超等[11]利用EIT监测和比较不同体姿下局部肺通气量和一秒率的改变情况，反映了体姿改变对肺通气的影响。王冠等[12]统计了不同体姿下健康志愿者EIT的基础范围值，可作为EIT系统研究和临床使用的参考依据。虽然这些研究小组利用EIT开展了不同体姿下的相关临床研究，但仍处于初步探索阶段，仅仅证明了EIT可以反映体姿变化所引起的肺通气变化，缺少量化指标的计算和具体的评价方法。EIT图像所包含的多种生理信息，需要根据预期的临床应用进行进一步的数据处理。

通气中心(CoV)、不均匀性(GI)、局部通气延迟(RVD)3种通气均匀性指标[13]和RVD分布图[14]是在EIT研究中衍生出的重要参数。CoV能够反映肺通气在腹背方向上的均匀性；GI计算了局部和整体的差异，能够更为精确地反映通气均匀性；RVD结合RVD分布图则能反映出肺通气区域内时间尺度上的不均匀性。利用3种通气均匀性指标和RVD分布图像可以从EIT图像中发掘出更深层次的信息，用于比较不同体姿下的肺通气均匀程度。

目前，临床上一般采用俯卧位通气以促进肺泡复张，改善氧合情况[15-16]。本文提出一种无线远程监控电阻抗成像(wireless telemonitoring electrical impedance tomography, WTEIT)实验方法，对不同体姿下的通气分布进行测定，采用3种指标和RVD分布图对通气均匀程度进行评价，用EIT方法证明了俯卧位通气在临床上应用的有效性，也为深入研究提供了动态监测的方法。

1 实验设备与方法

1.1 实验设备

本实验使用的WTEIT系统是自主研发的JLPulmo01-V1.0.0。该系统包括基于Red Pitaya开发板的数据采集模块、可穿戴式弹性缚带多电极传感器和上位机，如图1所示。

图1 电阻抗成像实验设备原理

Red Pitaya发出幅值可调的电压信号，控制压控电流源产生电流信号，电流信号通过模拟多路复用器和多电极传感器对受试者进行激励，采集到的电压信号，经差分放大回到Red Pitaya。在Red Pitaya内完成快速傅里叶变换，并计算出电压幅值，把电压幅值通过无线模块传到上位机，经数据预处理和成像算法进行图像重建、波形显示和数据存储等。电极数量为16电极，采用相邻激励、相邻测量的驱动测量方式，激励电流幅值在1～5 mA范围内可调，频率为122 kHz，成像速度20帧/s。该系统能够满足肺通气功能的床旁连续监测需求。肺功能仪采用日本某公司的AS-507。

1.2 实验方法

选取50名健康志愿者进行试验，其中男性33名，女性17名。选取标准如下：无心肺病史、无吸烟史、无饮酒史及无慢性疾病史；体内无移植性电子设备；无胸、脊柱畸形史；同意进行EIT测量者[12]。年龄、身高、体重和身体质量指数如表1所示。

表1 受检人群基线特征

具体实验步骤如下：

1)根据受试者的胸廓大小选择合适长度的电极带，用酒精擦拭受试者胸部测量区域，在电极上均匀涂抹少量导电膏，将电极带水平固定在第4-第5肋间处，同时保证电极间距均匀、与皮肤接触可靠，且8、9号电极间的基准线正对脊柱，连接EIT数据采集系统，开始监测；

2)受试者在坐姿下进行测试，平静呼吸一段时间后，待EIT采集数据稳定即可开始用力肺活量测试，并同步记录EIT数据，重复3次；

3)受试者依次变换为仰卧和俯卧，重复步骤2)的测试流程。

进行坐姿EIT监测实验如图2所示。

图2 坐姿EIT监测实验

采集到的EIT数据采用如下方式处理：

①上位机自动完成图像重建、波形显示功能；

②根据EIT系统绘制的曲线计算得到一秒率，并与由肺功能仪测得的结果进行比较，验证自主研发设备采集数据的有效性；

③将得到的EIT图像划分为腹侧和背侧两个感兴趣区域，分别统计腹侧和背侧的通气占比。

4)分别计算不同体姿下的3种通气均匀性指标，绘制RVD图，具体方法在下一小节中给出。

1.3 通气均匀性指标

采用通气中心(CoV)、不均匀性(GI)、局部通气延迟(RVD)3种通气均匀性指标和RVD分布图，可以对不同体姿下的肺部通气均匀性进行评估。与通常采用的根据像素计算方式不同，本文是根据节点来计算。

为了计算出这些指标，首先需要表示出功能EIT的潮汐变化(TV)，根据吸气末图像与呼气末图像来计算。一般取1 min为分析周期，并计算平均值以提高信噪比。计算公式如下：

(1)

式中：TVi表示EIT图像中i号节点的TV值；N表示分析周期内包含的呼吸次数；ΔZi,in,n和ΔZi,ex,n分别表示吸气末和呼气末对应的阻抗值。由于这些指标的计算需经过归一化处理，阻抗之差可用吸气末与呼气末的计算电导率差代替。

CoV表示肺背侧通气量与整体通气量之比，比值越大说明肺背侧部分通气越多。计算公式如下：

CoV=[∑(yi×TVi)/∑TVi]×100%

(2)

式中yi表示i号节点在腹背方向上对应的高度，该值被归一化到0～1之间，腹部对应0，背部对应1。

GI是反映各节点TV值差异的通气均匀性指标，是EIT图像中通气区域TV值的中位数与各节点之间差异总和，与通气区域所有节点TV值总和之比。该比值越小，说明肺通气越均匀。计算公式如下：

GI=∑l∈lung|TVl-median(TVl)|/∑l∈lungTVl

(3)

式中归一化后的TV值<0.4部分对应的节点表示在肺通气区域内，其TV值用TVl表示。

RVD是通过计算各节点通气时在时间上的差异得到的参数，其标准差记为RVDSD，RVDSD越小，说明肺通气越均匀。根据EIT肺通气区域不同节点通气情况分别绘制出局部阻抗变化曲线，并由曲线进行计算RVD值。计算公式如下：

(4)

式中：tl,40%表示从吸气开始达到最大阻抗40%所需的时间；tin-ex表示根据全局阻抗曲线计算的吸气总时间；L表示通气区域包含节点的总个数。

另外，可以根据肺通气区内各节点的RVD生成RVD分布图。RVD代表通气延迟程度，RVD越大，代表通气越迟。使通气过程中时间尺度上的不均匀性得到直观反映，具体方法如图3所示，图中以A、B、C3点为例。这3点的RVD指数依次减小，说明在时间尺度上A点通气最迟，实现了区域肺泡通气延迟的可视化。

图3 RVD分布图

2 实验结果

2.1 不同体姿下用力肺活量测试

肺功能检查是对受试者的整体肺通气功能进行诊断和评估的金标准，健康人的一秒率>70%[17-18]。一秒率是其中的重要指标，常用于判断阻塞性气道疾病的严重程度[19]。将EIT系统与肺功能仪测得的一秒率进行比较，验证自主研发设备采集数据的有效性。

某一受试者坐姿下用力肺活量期间EIT测试结果如图4所示。图4(a)中反映了呼吸过程中208个通道的边界电压之和随时间变化的关系。吸气过程中，随着空气的吸入，肺部阻抗增大，边界电压和随之增大；呼气时，边界电压和不断减小。在两次平静呼吸之后，开始进行用力肺活量测试，用力吸气至t=7.4 s时，达到吸气最大值，随后尽力尽快呼出肺内所有气体。这部分呼气量记为用力肺活量(FVC)，7.4 s～8.4 s这1 s内的呼气量记为一秒钟用力呼气容积(FEV1)。FEV1与FVC之比即为一秒率(FEV1%)，此受试者的一秒率为82.8%。图4(b)反映了用力肺活量测试呼气前1 s内的EIT图像。此过程中，肺部截面积不断减小，说明肺部气体不断减少，与图4(a)中的曲线相对应。

图4 受试者坐姿下用力肺活量期间EIT测试结果

在坐姿、仰卧和俯卧3种不同体姿下分别进行用力肺活量测试，其余条件保持相同。用EIT系统和肺功能仪测得的一秒率结果分别用平均值±标准差形式给出，如图5所示。结果表明志愿者不同体姿下的一秒率平均值都在80%以上。EIT系统测得的一秒率结果与用肺功能仪测得的结果具有良好的一致性，说明该系统采集数据的有效性。

图5 不同体姿下EIT系统与肺功能仪一秒率测试结果

2.2 肺部通气均匀性评价

对坐姿、仰卧和俯卧3种不同体姿下腹侧和背侧的通气占比进行统计，结果用平均值±标准差形式给出，如图6所示。3种体姿下，腹背侧肺通气占比分别为(61.25±9.61)%∶(38.75±9.61)%、(66.13±9.27)%∶(33.87±9.27)%和(57.54±7.89)%∶(42.46±7.89)%。其中，俯卧时两者之比最接近1∶1，说明俯卧时的通气分布最均匀。结果表明：EIT能够反映体姿变化导致的肺通气分布和肺排空能力的改变。

图6 不同体姿下腹背侧通气比例结果

根据计算公式，得到不同体姿下3种通气均匀性指标的计算结果如表2所示，用平均值±标准差形式给出。

表2 3种体姿下的通气均匀性指标结果

图7以雷达图的形式给出结果，直观地反映出指标大小关系，并根据肺通气区域内各点的RVD生成RVD分布图。结果表明：CoV在俯卧时最大，代表俯卧时肺背侧通气最多，与通气占比结果吻合；GI在俯卧时最小，代表俯卧时肺通气最均匀；RVDSD在俯卧时最小，也说明了俯卧时通气最均匀。图中的灰度深浅分布清晰地反映出区域通气延迟的程度。俯卧时，肺通气区域范围内的颜色分布变化较为平缓，而对应的RVDSD在俯卧时最小，进一步说明了俯卧时通气最均匀。

图7 不同体姿下通气均匀性指标和RVD分布图

3 讨论

本文采用3种通气均匀性指标对不同体姿下的通气情况进行评价。CoV不需要经过复杂的计算，就能反映重力对通气分布的影响。KARSTEN J等[20]用实验证明：基于EIT的CoV监测可以评估麻醉、气腹和不同PEEP水平的效果。GI与CoV相比，更好地反映了局部与整体的差异。ZHAO Z Q等[21]利用GI来选择与通气均匀性相关的PEEP水平。RVD则体现了肺在通气过程中时间尺度上的局部差异。MUDERS T等[14]利用RVD估计潮气补充量，可能有助于个性化的通气设置。这些指标在临床上已有简单的应用，但仍需要进一步深入研究。

3种体姿的EIT实验证明了俯卧位通气在临床应用上的合理性。不同区域的肺通气量会随体姿改变而变化，且与重力影响有关。俯卧时，腹侧区域通气减少，而背侧区域通气增加，相对另外两种体姿，肺通气分布更均匀。从生理角度分析，主要有以下几点原因[15-16, 22]：1)俯卧时腹侧区域血流增加而背侧区域血流减少，通气血流比改善，分流减少；2)俯卧能减少心脏、腹部对肺的压迫，促进重力依赖区肺复张；3)俯卧时重力依赖区和非重力依赖区变化导致腹背侧顺应性变化，而背侧通气区域顺应性增加相较于腹侧减少明显，使整体顺应性增加。因此，临床上采用俯卧位通气方式。

4 结语

本文提出一种无线远程监控电阻抗成像实验方法，对健康人开展肺通气数据采集，对不同体姿下的通气分布进行测定，采用CoV、GI、RVD3种指标和RVD分布图对通气均匀程度进行评价。得到以下结论。

1)EIT能够反映体姿变化导致的肺通气分布和肺排空能力的改变，俯卧时腹背侧通气占比为(57.54±7.89)%∶(42.46±7.89)%，最接近1∶1。

2)3种指标和RVD分布图在不同体姿下能反映肺通气分布的均匀性。俯卧时CoV最大，为(55.26±1.10)%，说明背侧通气最多；GI最小，为0.25±0.03，说明通气最均匀；RVDSD最小，为5.01±1.43，也说明通气最均匀。

3)俯卧相比坐姿和仰卧，能够使肺部通气更加均匀。因此，临床上采用俯卧位通气以促进肺泡复张，改善氧合情况。