lead transfer function

概念

植入导线在 MRI RF 场中的电极处发热，不能用单个集总回路描述。导线长度与组织中波长同量级（lead-resonant-length），沿线各处拾取的电场对电极的贡献既不等大、相位也不同。传递函数法把这层关系写成一个线性泛函。

分布模型

把导线视为浸在有耗介质中的传输线（方程本身、通解与反射系数见 transmission-line-equations）：每一微段 $\mathrm{d}z$ 有串联阻抗（导体电阻与自感）和对组织的并联导纳（绝缘层电容与组织电导）；入射切向电场沿线逐段注入分布电动势。两端由边界条件封闭——近端是 IPG 外壳对组织的阻抗，远端是电极与组织的界面阻抗。

与单回路集总模型的区别在并联支路：有了逐段对组织的返回路径，导线各段才不必共用同一个回路电流，“某一段拾取的场对电极贡献多少”才成为可问的问题。

传递函数

在该线性模型下，电极处的复响应 $A$ 可写成沿路径的加权积分：

$$A={\int }_0^{L}S(z)E_{\tan }(z)\mathrm{d}z$$

逐符号说明：

• $z$：沿导线路径的弧长坐标，单位 m；$z=0$ 取近端，$z=L$ 取远端电极。
• $E_{\tan }(z)$：入射切向电场沿该路径的复振幅（幅值与相位），单位 V/m；“入射”指未放入器件时该路径上的场。
• $S(z)$：复值传递函数，表示 $z$ 处单位切向电场对电极响应的贡献；量纲随 $A$ 的取法而定。
• $A$：电极处的复响应，可取为流入组织的电极电流，或与电极附近局部 SAR 平方根成正比的量。局部功率沉积与温升正比于 $|A{|}^2$。

为什么用它

$S(z)$ 只由器件本身及其周围介质决定，与 MRI 线圈、患者体型、导线走行无关。测定一次即可对任意临床路径的 $E_{\tan }(z)$ 复用，把”人体里的场”与”器件的响应”这两件难事解耦。ISO/TS 10974 的 Tier 3 正是这样分工：人体—线圈仿真给出各临床路径的复数 $E_{\tan }(z)$，AIMD 的等效电磁模型给出响应，二者结合得到各热点的功率沉积分布（iso10974-rf-heating-2-tiers）。

获取 $S(z)$ 的途径有二：数值上用矩量法计算；实验上沿导线逐点局部激励扫描测量。无论哪条，器件等效模型都须另行验证后才可用于评估。

谐振长度效应藏在哪里

写成积分形式后，“某些长度耦合最强”这件事并没有消失。与激励无关的是 $S(z)$ 里的公共往返因子 $1/\left(1-{\rho }_p{\rho }_t{\mathrm{e}}^{-2\gamma \ell }\right)$——它只由线本身与两端负载决定，其取极大的条件 $2\beta \ell \equiv \arg ({\rho }_p{\rho }_t)(\mathrm{mod}2\pi )$ 就是谐振长度。该因子、终端电流的闭式解与长度扫描算例见 field-excited-transmission-line，本卡不重复。

但峰是否出现、出现在哪，还要看 $E_{\tan }(z)$ 的相位分布与该因子如何配合——这一层才是传递函数方法要处理的，下面用两种激励对照说明。

field-excited-transmission-line 里那条长度扫描曲线是”均匀同相激励”专属的结果，不是导线的固有属性。 把激励换成带相位推进的行波场 $E_{\tan }(z)=E_0{\mathrm{e}}^{-\mathrm{j}{k}_{t}z}$（64 MHz、${\epsilon }_r\approx 70$ 时 $k_t\approx 11.2$ rad/m，而线上 $\beta \approx 5.09$ rad/m，二者失配约 2.2 倍），同一根线上：

长度 (m)	均匀同相 $|I_{\mathrm{tip}}|$ (mA per V/m)	行波激励 $|I_{\mathrm{tip}}|$
0.34（均匀场下的峰）	2.87	1.91
0.70（均匀场下的谷）	0.18	0.86
0.958（均匀场下的峰）	1.37	0.29
1.57（均匀场下的峰）	1.56	1.32

有的峰塌掉、有的谷填起来：相位失配使沿线贡献部分相消，改变的不只是幅值，还包括峰的存在与位置。所以均匀同相假设在 64 MHz、数十厘米量级的植入导线上站不住——沿 0.4–0.6 m 的组织波相位推进已超过一整圈。它的正当用途只是教学示例：用最简的激励把 $S(z)$ 的结构解出来看。真实评估必须用各临床路径的复数 $E_{\tan }(z)$。

上表两列均由同一组方程解出并经 Wolfram 独立数值核验；行波激励一列的闭式解与符号解逐项相符。

适用边界

• 线性、细长、可参数化为一维路径的结构。分叉导线、大金属体、多导体间强耦合需另作处理。
• 相位不能丢。$E_{\tan }(z)$ 是复量，不同段的贡献会相消；按幅值相加会系统性高估。用均匀同相场代替真实场也不保守——它既可能高估也可能低估，最坏相位分布下的耦合比它高出数倍。
• $S(z)$ 不是器件的不变量。改变端接阻抗、加入滤波器、更换绝缘层都会改写 $S(z)$。
• 它给出的是电极处的电磁驱动量，从该量到组织温升还须经热传导与灌注模型。

与滤波器位置的关系

在导线某处串入高阻抗，会压低该点近端一侧的 $S(z)$：那一段拾取的场不再能驱动电极电流，只能经各段自身对组织的电容返回。于是滤波器越靠近电极，$S(z)$ 中仍能驱动电极的那一段越短，电极电流越小。

这条位置依赖在只有单个串联回路的集总模型里根本无法表达——串联阻抗放在同一回路的任何位置，效果完全相同。

同一根 50 cm 导线的三种配置（LTspice 原理图）。每一节含串联 $L^{\prime }\mathrm{d}z$、该节自身的感应电动势 $E_{\tan }\mathrm{d}z$、以及一只回到组织的并联电容 $C^{\prime }\mathrm{d}z$；左端是 IPG 外壳对组织的电容，右端是电极—组织界面电阻。A 无滤波器；B 带阻串在电极端；C 带阻串在 IPG 端。元件值均为假设值，图用于说明并联返回路径如何使”位置”成为一个可问的问题。

参考

• Park S-M, Kamondetdacha R, Nyenhuis JA. “Calculation of MRI-induced heating of an implanted medical lead wire with an electric field transfer function.” J Magn Reson Imaging 2007;26(5):1278–1285. doi:10.1002/jmri.21159 —— 定义电场传递函数，以切向电场与传递函数沿导线长度的乘积积分求电极附近电场；传递函数由矩量法数值求得，并对 64 MHz 下不同长度的裸导线与端部裸露的绝缘导线给出结果。
• Yeung CJ, Susil RC, Atalar E. “RF safety of wires in interventional MRI: using a safety index.” Magn Reson Med 2002. PubMed 11754458 —— 给出可预测导线尖端加热绝对值的理论模型，用矩量法计算绝缘与裸导线上感应电流造成的局部 SAR 放大。
• Agrawal AK, Price HJ, Gurbaxani SH. “Transient response of multiconductor transmission lines excited by a nonuniform electromagnetic field.” IEEE Trans Electromagn Compat 1980;EMC-22(2):119–129. doi:10.1109/TEMC.1980.303824 —— 场对传输线耦合的经典公式化之一，由 Maxwell 方程导出非均匀入射场激励下的传输线方程，其中入射切向场以分布串联源的形式进入。它提供的是本卡所解的那组方程（单导体为其 $n=1$ 情形），不提供上面的闭式解与反射系数分解——后者是本卡对该方程组在均匀场、两端集总负载下自行求解并核验的结果。
• iso10974-rf-heating-2-tiers —— Tier 3 的输入即各临床路径的复数 $E_{\tan }(z)$ 与经验证的 AIMD 等效模型；本卡描述的是该 tier 所依赖的方法族。

关联

• 相关概念：lead-resonant-length、transmission-line-equations
• 相关源卡：US7363090B2_Greatbatch、US9999764B2_Medtronic
• hazard:rf-heating · solution:filtering