起搏器携带人员在MRI环境中的轨迹特异性运动感应电动势评估

《MEDICAL PHYSICS》：Trajectory-specific motion-induced electromotive force assessment of pacemaker-carrying workers in MRI environments

　　

目的：携带主动植入式医疗器械（AIMD，如起搏器（PM）和植入式心律转复除颤器（ICD））的个体在磁共振成像（MRI）系统内部及周边空间变化的静态磁场中运动时，会在导电导线中感应出电动势（EMF）。该机制在标准MRI安全评估中未被明确涵盖。本研究提出一种基于测量的计算框架，用于在与研究临床实践相关的代表性暴露条件下估算运动感应EMF。方法：在3 T和7 T MRI系统的边缘区域测量的静态磁场分量（Bx，Bz）于规则网格上插值。推导空间梯度并用于计算代表PM/ICD导线配置的等效导电回路中的运动感应EMF。考虑了三种运动类型：沿扫描仪轴（z）的线性运动、垂直于孔径（x）的线性运动以及回路的原位旋转。使用法拉第定律计算EMF，假设回路面积为0.0225 m2、参考线速度为1 m/s、角速度为1 rad/s。获得每种运动场景的空间EMF分布和时间波形。结果：感应EMF随运动方向和空间位置显著变化。对于线性运动，最高值出现在孔径入口附近空间场变化强烈的区域。在3 T时，最大EMF达到160 mV（z向运动）和114 mV（x向运动），而在7 T时相应值为81 mV和104 mV（针对所分析轨迹）。旋转运动产生较低且更局部的EMF值，3 T和7 T的最大值分别为30 mV和35 mV。尽管标称场强更高，7 T系统并未一致产生更高的EMF，这反映了空间场梯度的差异。结论：MRI边缘场中的运动感应EMF取决于局部空间场梯度和运动特性，而非仅取决于标称场强。所提出的基于测量的框架提供了空间分辨的、特定于场景的EMF估计，并在给定适当场图的情况下可应用于不同的MRI系统。该方法通过解决标准评估中未明确涵盖的运动相关暴露条件，补充了现有的安全评估。

研究背景方面，主动植入式医疗器械（AIMD，Active Implantable Medical Devices）包括心脏起搏器（PM，Pacemaker）和植入式心律转复除颤器（ICD，Implantable Cardioverter Defibrillator）在患者和医护人员群体中日益普遍，同时磁共振成像（MRI，Magnetic Resonance Imaging）的临床应用持续扩展，高场及超高场（3 T和7 T）系统不断部署。虽然MR条件设备和安全协议显著降低扫描患者风险，但MRI设施的更广泛电磁环境对携带AIMD个体尤其在孔径外边缘场区域仍是关切点。欧洲指令2013/35/EU规定雇主须对被视为“特别风险工人”的AIMD携带工人进行个体EMF（Electromotive Force，电动势）风险评估，国际标准和指南也定义了MRI环境暴露限值和程序。现有AIMD安全评估多聚焦接受成像的患者，关注射频加热、梯度感应电压和装置故障，大多数PM/ICD系统标记为MR条件但仅基于患者扫描风险评估，不能外推至在MRI环境内工作的场景。医护人员在静态磁场梯度中移动可能产生时变磁通并在导电回路（含植入导线系统）中感应EMF，源于身体或装置运动与静态磁场空间梯度的相互作用，依赖局部场梯度、回路几何、取向和运动动力学（平移或旋转），感应电压在装置输入级可能干扰导致不当感知、起搏抑制或意外刺激，对常规接近孔径、俯身或移动的操作者尤其相关。以往体内测量复杂，研究多用体外模型和数值模拟，更高场MRI系统引入使定量表征边缘场梯度及相关运动感应EMF更重要，标称场强不足以决定风险，因静态场空间分布及梯度依赖磁体设计和屏蔽配置，故需结合现场特定场映射与基于物理的EMF估计来评估高暴露场景。该研究发表于《MEDICAL PHYSICS》，研究人员提出使用测量静态场图估算PM导线回路中运动感应EMF的计算框架，结合空间分辨磁场数据与轨迹依赖运动模型（含线性行走路径和原位躯干旋转），在两种MRI环境（3 T和7 T）应用以量化代表性高暴露场景的感应EMF并识别操作者移动相关的空间风险区。

研究人员开展研究用到的主要关键技术方法包括：获取胸高水平面静态磁场分量B_x(x,z)和B_z(x,z)的高分辨率二维图作为计算框架输入，PM导线建模为符合植入安全建模的等效刚性导电回路，回路面积取ISO 14117附录L推荐的225 cm²平面最坏情况；基于法拉第定律一般形式计算感应EMF，假设回路表面磁场局部均匀，磁通积分简化为回路面积乘垂直回路平面的磁场分量，用中心有限差分在测量网格计算磁场分量空间导数，测量数据经立方插值至0.01 m网格后求导以提高数值稳定性；线性运动恒定速度下EMF简化为沿运动方向垂直回路平面的磁场分量方向导数，旋转运动绕垂直轴建模引入时变回路法向量由转角ωt定义、回路中心固定；采用三维轴向霍尔探头磁力计在MRI室实验测量静态磁场输入，于胸高x–z平面以0.1 m间距规则网格采样边缘场区，测得数据在右手笛卡尔坐标系（原点为扫描仪等中心点，z轴沿孔径和患者床方向，x轴水平左右向，y轴垂直向）下表达并立方插值为0.01 m精细规则网格后求空间梯度；开发MATLAB交互仿真环境集成场图与图形界面供定义运动条件（类型、方向、速度、回路面积等）并自动计算EMF；框架应用于临床3 T和超高场7 T系统，独立对每数据集做相同数值流水线（预处理、插值、梯度计算）及相同几何与运动参数下的线性平移与旋转运动仿真以直接比较不同静场强与边缘场梯度结构的EMF水平与空间分布。

研究结果部分保留原文小标题并说明研究与结论如下：

3.1 Spatial field and gradient characterization（空间场与梯度表征）：研究人员测量两系统胸高静态磁场图显示距孔径增加场衰减，3 T系统在机架区近场分量局域空间变化更陡，7 T系统分布更广衰减更平滑；由插值分量图计算磁场幅度梯度空间图确认两扫描仪近孔边缘区局域最高空间变化，空间范围与峰值不同。梯度图仅支持运动感应耦合的物理解释，EMF计算严格基于回路表面垂直磁场分量空间导数。文献指出峰值梯度可能在测量网格外区域更高取决于系统设计。

3.2 Motion-induced EMF—Linear translations（运动感应EMF—线性平移）：研究人员模拟PM回路以1 m/s在测量边缘场线性运动，得到相同参数下峰EMF因运动方向与系统而异：沿扫描仪轴（z向）3 T为159.1 mV，7 T为80.8 mV；垂直扫描仪轴（x向）3 T为113.5 mV，7 T为103.7 mV。3 T系统z向运动峰EMF高于x向，与观测更大空间梯度max |?B_z/?z|=7.07 T/m vs max |?B_x/?x|=5.04 T/m一致；7 T系统x向高于z向，与max |?B_z/?z|=3.59 T/m vs max |?B_x/?x|=4.61 T/m一致。高暴露EMF位置自动识别于机架近区空间梯度最大处；EMF图反映测量磁场分量及其空间导数局域结构，高EMF区对应相关分量梯度快速变化区，受非均匀边缘场分布与测量插值影响，但整体趋势（孔口及离轴高危区）稳健一致于底层测量梯度。

3.3 Motion-induced EMF—rotational motion（运动感应EMF—旋转运动）：研究人员以1 rad/s绕垂直轴旋转躯干固定回路中心，相同面积下峰EMF为3 T 29.8 mV，7 T 35.5 mV。旋转EMF不直接随边缘梯度幅值缩放，而主要取决于选定位置局域矢量场组成，故7 T略高于3 T；旋转期间EMF时间历程呈正弦行为，与磁场矢量同旋转回路法向点积的理论依赖一致。

3.4 Representative high-exposure scenario spatial mapping（代表性高暴露场景空间映射）：研究人员总结固定参数下代表性高暴露场景感应EMF及空间坐标，EMF值局域于孔口近区高梯度边缘场区，坐标自动提取叠加于场与梯度图以追溯空间场结构与感应电压极大关系；空间EMF图显示对运动方向强各向异性，证实暴露评估需联合建模回路取向与运动方向。等效时间磁场变化率dB/dt=EMF/A在3 T线性运动约5–7 T/s，7 T约3–5 T/s，旋转两系统约1–2 T/s，为代表高估算非绝对上限。

3.5 Cross-scanner comparison under matched motion conditions（匹配运动条件下跨扫描仪比较）：研究人员在相同运动参数与回路几何下直接比较3 T和7 T系统，发现更高标称场强不一定产生更高运动感应EMF，感应电压主要由局域边缘场空间梯度与矢量场组成决定而非仅B₀场强；本测量数据中3 T扫描仪两主方向线性平移峰EMF更高，而7 T旋转峰EMF在测试位置较高；7 T系统可及区边缘场更受限且空间梯度更低，致尽管标称静场更强但运动感应EMF减小；相同旋转运动学下两系统EMF时间进程因各自局域磁场矢量组成不同而幅相特性各异，旋转EMF由场矢量在时变回路法向投影决定。

讨论部分总结：研究人员结果表明PM等效回路在MRI环境内运动感应EMF强烈依赖运动类型与局域磁场空间结构；7 T场景中垂直患者床（x向）线性运动峰EMF最高，其次平行（z向），3 T场景中平行最高其次垂直，两场景旋转运动在同等面积与速度下幅值显著更低；关键发现为7 T系统未系统性产生高于3 T的感应EMF，z向线性运动峰EMF约低一半，旋转幅值相当，证实运动感应EMF主要由局域空间梯度与矢量场组成决定而非仅标称B₀场强，3 T安装中更强边缘场空间变化产生更大有效dB/ds值，与文献记载1.5、3、7 T临床系统梯度上限约19、17、16 T/m一致，说明低场系统因磁体与屏蔽配置可能具可比拟甚至更高运动感应EMF，评估应基于系统特定空间场分布而非仅标称场强，本测量框架可推广至1.5 T等系统；x与z平移方向依赖符合法拉第定律，垂直回路表面磁通驱动分量为法向磁场分量，z向运动耦合?B_z/?z，x向耦合?B_x/?x（躯干转向使回路平面旋转），7 T中x向更高EMF表明测量图中横向梯度陡于纵向；旋转运动峰EMF更低因其源于取向变化而非空间梯度遍历，EMF随角速度与局域面内磁场矢量幅值缩放，对局域梯度热点敏感性降低；需注意成像梯度切换产生时变磁场感应的EMF通常高于静态场中运动感应EMF，但梯度暴露受控且标准明确涵盖，而运动感应EMF取决于B₀场中移动且可能发生于非活跃成像期（如患者搬运），二者互补，运动感应EMF评估提供标准未完全覆盖的暴露条件洞察；本研究用等效单回路模型与有效面积抽象捕捉主导磁通耦合，忽略实际导线走行与返回路径细节，可按植入配置调整；假设回路表面磁场局部均匀使磁通积分简化为面积乘法向分量，对典型PM/ICD回路尺寸与MRI边缘场梯度一般满足，局域高梯度区或有二阶修正；磁场输入为胸高二维图，未显式建模垂直变化与全三维效应，结果应视为基于可用测量的场景特定评估而非广义保守估计；未建模组织电导、导线阻抗与装置电路响应，依ISO 14117、EN 50527-1、EN 50527-2-1标准估算开路EMF而非感应电流；运动建模为刚体恒速平移或旋转，未来可纳入复杂生物力学轨迹；设备安全角度计算峰EMF在数十至数百毫伏（面积225 cm²、1 m/s、1 rad/s），可与高于典型心脏植入电子设备感知阈值比拟，提示边缘场区快速移动可能为不可忽略耦合机制，现代PM导线含滤波保护但仍需运动感知暴露评估；感应EMF频率约1 Hz低于标准免疫测试起始16.6 Hz，难直接与PM/ICD编程灵敏度阈值（单极约2 mV，双极约0.3 mV）及低频免疫电平2 mV（峰峰值至1 kHz恒定）比较，后者仅作初步基准参考；双极感知感应面积较单极小约17倍，但框架预测EMF仍远高于灵敏度设置与免疫电平，确认高空间变化静态场中运动可诱发装置故障（过感知），与文献体外结果一致；框架基于测量场图与闭合法拉第定律，无需全波求解器、解剖模型与专有机密信息，交互实现允许非专家选位与运动获时变EMF，适合作MRI设施筛选与比较风险评估工具以识别高危区与模式；方法基于地图且特定于扫描仪，支持现场依赖评估与更新；未来可精细化用多高度场测量与三维插值重建回路平面磁场分布做表面积分通量评估以估空间非均匀性，需更密体积场映射；可耦合受控体模实验（导电躯干体模与仪器化导线）关联预测EMF与装置行为（感知扰动等）以验证并推向装置级风险表征；本研究聚焦确定性轨迹以可控解释关系，未捕获真实操作者运动变异性，未来可拓展至轨迹集合与概率运动模型以定义暴露包络与统计描述（如EMF百分位数）；当前优先正面区（主要交互区），侧与后区需专门测量，弯腰等耦合平移旋转需全三维场映射；计算运动感应EMF是PM携带者在MRI环境安全评估中间步骤，实际风险还含导线配置、装置易感性与滤波特性，结果应解释为电磁水平暴露定量估计，可作高阶装置响应与临床风险模型输入，未来将耦合本EMF框架与概率轨迹及装置级响应模型以更全面评估职业安全。

结论部分原文翻译：本研究提出一种计算框架以估算MRI边缘场环境中PM/ICD等效导线回路内的运动感应EMF。结果显示感应EMF主要取决于运动类型、轨迹及局域磁场空间梯度，而非仅依赖标称扫描仪场强。所提方法结合测量场图与感应模型及交互式运动定义，能够以最小计算复杂度实现位置分辨与轨迹特定的EMF估计。该方法特定于扫描仪、透明且易实施，可作为MRI设施的实用筛选与比较风险评估工具，支持识别AIMD携带者及暴露操作者的更高风险运动模式与空间区域。框架基于一阶等效回路近似，但易于扩展至更高分辨率场积分与实验体模验证，这些特性使该工具成为理论感应建模与现代高场装置中操作性MRI安全评估间的有用桥梁。