膈神经电刺激被用作一种人工通气方法，计算模型被用于辅助参数选择。大多数模型假设组织电导率为各向同性，尽管肌肉组织表现出各向异性特征。研究人员旨在研究各向异性肌肉电导率对膈神经激活结果的影响。
为计算电位分布，研究人员使用了解剖学精细的多尺度模型进行颈部非侵入性电刺激，该模型融入了真实的肌纤维走向。膈神经激活阈值通过McIntyre-Richardson-Grill（MRG）神经模型计算。在恒定相应各向同性电导率条件下，分析了从1:1到1:15（横向:纵向电导率）的各向异性比值。额外模拟评估了肌肉体积和电极放置位置的影响，并对颈部其他神经的可能共同激活进行了量化。
增加各向异性比值导致所有轴突直径的膈神经激活阈值持续升高（最高达+90%）。较大的肌肉体积和直接位于肌肉上方的电极位置进一步增强了各向异性效应。考虑各向异性肌肉电导率会增加共同激活神经的数量。
研究建议，在膈神经刺激的研究中采用纳入各向异性电导率的高分辨率模型。
本研究聚焦于非侵入性膈神经电刺激中肌肉组织各向异性电导率对神经激活阈值的影响，为电刺激辅助通气的计算模型优化提供了重要依据。
**一、研究背景与问题**
机械通气是重症监护中挽救生命的治疗手段，但患者吸气肌肉未被主动利用，可能导致呼吸机诱导的膈肌功能障碍。膈神经刺激可维持肌肉活动并促进生理性呼吸运动，其中非侵入性电刺激策略因无需手术植入电极而具有临时干预优势。为实现高效刺激并避免脱靶激活，需仔细选择刺激参数，计算模型和电磁仿真是辅助参数选择的重要工具。然而，现有模型多假设组织电导率为各向同性，忽略了肌肉组织实际具有的各向异性特征——即电流沿肌纤维纵向的传导优于横向传导。文献报道的骨骼肌各向异性比值范围广泛，从接近1（各向同性）到15不等，且测量频率也是影响因素。尽管白质各向异性对脑电图、经颅电刺激等的影响已被证实，肌肉各向异性对膈神经激活阈值影响的系统研究尚缺乏。考虑到颈部解剖结构复杂、肌纤维方向多样，建立纳入真实各向异性肌肉电导率的精细模型对准确预测刺激效果至关重要。
**二、研究方法与关键技术**
研究人员采用Weget等人先前发表并进一步开发的颈部多尺度有限元模型，该模型基于BodyParts3D平台构建，包含13种组织类型，边界为C1至T1椎体水平。研究对象为标准电极配置（位于胸锁乳突肌后缘附近，阴极位于环状软骨水平，阳极位于其上方3.5 cm处，附加接地电极置于颈后部），以及向前后各移动2 cm的两组变位电极。为表征肌肉解剖变异，基于Migotto等人报道的胸锁乳突肌体积变异性，创建了75%、90%、110%和125%四种体积比例的颈部模型。
关键技术创新在于利用COMSOL Multiphysics的曲线坐标接口，通过流场模拟计算肌纤维方向的局部坐标系——以肌肉近端为入口、远端为出口进行流场仿真，生成的流线作为肌纤维走向的近似，从而将各向异性电导率张量与局部纤维方向对齐。纳入各向异性的肌肉包括刺激侧距膈神经和电极20 mm以内的19块肌肉。电导率张量保持体积守恒的椭球约束，各向异性比值设置从1:1（各向同性）递增至1:15。宏观模型的标量电位通过二阶四面体单元计算（约1.27亿自由度），提取的膈神经三束膜外电位导入NEURON环境的MRG有髓轴突微模型，以150 µs单相阴极脉冲为参考，迭代缩放至触发动作电位的最小电流即为激活阈值。选择性评估以激活膈神经所需电流为操作点，对臂丛（3段）、颈袢、耳大神经、胸长神经、迷走神经、锁骨下神经、锁骨上神经（2段）和颈横神经（2段）进行二元激活判定，采用12.8 µm、10 µm和7.3 µm三种代表性纤维直径。
**三、研究结果**
**3.1 膈神经激活与选择性**
标准电极配置和未缩放肌肉体积条件下，各向异性比值增加导致所有纤维直径的激活阈值升高。14 µm纤维在1:1比值时阈值最低（21 mA），7.3 µm纤维在1:15比值时最高（149 mA）。10 µm纤维阈值从37 mA（1:1）升至69 mA（1:15）。相对增量呈比值依赖性，1:2比值时中位增加9.7%，1:15时达83.8%。选择性方面，在12.8 µm膈神经纤维激活阈值下，1:1比值时已招募耳大神经、锁骨上神经、锁骨下神经和颈横神经的全部节段，以及臂丛一节段、颈袢和胸长神经（大直径纤维），共9个脱靶节段；1:15比值时增至11个。较小纤维激活时额外招募更多深部结构，选择性进一步降低。共同激活主要发生在解剖位置表浅、靠近电极的神经。
**3.2 改变电极位置**
电极前移至胸锁乳突肌上方使10 µm纤维阈值从38 mA微增至39 mA（1:1），后移至斜方肌上方则大幅增至90 mA。各向异性效应因位置而异：前位置相对增量12.3%（1:2）至89%（1:15），后位置为12.7%至98.1%，相比标准位置增幅更高。后位置最大四分位距达15.4%。
**3.3 肌肉体积变化**
各向同性条件下体积变化对阈值影响小于1 mA。各向异性条件下，肌肉体积增加使阈值进一步升高：10 µm纤维在1:15比值时，125%体积阈值72 mA，75%体积降至59 mA，相对增量分别为90.3%和58.2%。体积差异效应随各向异性比值增大而增强。
**四、讨论与结论**
研究人员分析表明，肌肉各向异性通过改变电流走向显著影响膈神经激活阈值。纵向电导率增加使电流更趋向沿表浅肌纤维流动，形成"屏蔽效应"，使得到达深部膈神经的有效电流减少。这一机制解释了为何电极位于厚肌群上方、或患者肌肉体积较大时，需要更高刺激强度。临床意义上，简化各向同性模型可能低估实际需求刺激量，而忽视各向异性还可能导致对共同激活风险估计不足——共同激活评估显示，表浅感觉神经首当其冲，但深部运动神经的意外激活亦需警惕；迷走神经因阈值较高而风险相对可控。
研究局限性包括：未纳入对侧肌肉各向异性、基于流场模拟的纤维方向为近似替代（非DTI直接测量）、未建模肌羽状角、以及代表纤维直径的简化处理等。未来研究方向指向个体化模型构建——从基于颈围、皮下脂肪厚度等简易参数的粗略调整，到基于分割MRI和DTI数据的高精度建模，以优化刺激方案设计，例如采用多电极阵列、电流转向、时域干涉刺激等技术提升选择性与安全性。
**研究结论：** 本研究揭示了颈部肌肉各向异性电导率对非侵入性膈神经刺激的重要影响，各向异性比值、肌肉体积和电极位置均显著调控激活阈值和共同激活模式。为最小化建模误差，建议后续模拟中纳入肌肉组织的各向异性电导率。在电刺激辅助通气的研究中，采用考虑各向异性容积传导的高分辨率个体化模型具有重要价值；临床实践中，应避免将电极直接置于肌肉上方，并进一步优化刺激协议以减少共同激活。
