基于目标场点法和流函数的磁共振有源匀场线圈设计方法*

黄清明 陈珊珊 张建青 杨洋 郑刚

1) (上海理工大学光电信息与计算机工程学院, 上海 200093)

2) (上海健康医学院医学影像学院, 上海 201318)

3) (上海市分子影像重点实验室, 上海 201318)

4) (上海理工大学医学影像工程研究所, 上海 200093)

磁场均匀性是磁共振系统的重要参数, 提高磁场均匀性有助于磁共振时域信号的检测和磁共振频域信号分辨率的改善.基于有源匀场连续电流密度分布的思想, 采用目标场点法和流函数结合的方法设计匀场线圈, 即由毕奥-萨伐尔定律确定磁场分布与电流密度的关系, 约束线圈半径和设置约束点后, 根据目标场分布逆向求解线圈平面的电流密度分布, 再用流函数将电流密度分布离散化处理, 得到匀场线圈的绕线位置分布.根据电磁仿真计算结果制作包含一阶与二阶匀场线圈应用于磁共振分析仪, 实验验证表明该匀场线圈能有效地改善永磁体核磁共振系统磁场的均匀性.

1 引 言

核磁共振分析仪的永磁磁体由两个平面磁极相对的偶极磁体组成, 其主要作用是产生稳定、均匀的主磁场, 磁场均匀性的提高有助于磁共振时域信号的检测和磁共振频域信号分辨率的改善[1].而影响永磁磁体磁场均匀性的主要因素有磁体本身的非均匀性、周围铁磁性物质的干扰、被检测样品的磁化率突变以及磁体温度波动等, 对磁共振系统的磁体进行匀场可以改善其均匀性[2,3].永磁磁体的匀场包括无源匀场和有源匀场两种方式, 其中有源匀场技术通常是在无源匀场的基础上对磁场均匀性进行优化, 通过放置在磁极中各向通电线圈所产生的磁场来补偿固有磁场的非均匀度, 即在磁体两个平行极面上布置匀场线圈, 匀场电源给不同方向的匀场线圈施加一定大小的电流以产生合适的磁场, 从而对主磁场的非均匀性进行精细修正.

磁场源除了磁场强度要求外, 在一定空间范围内还应满足磁场均匀度要求, 通常将磁场均匀度很高的赫姆霍兹线圈作为标准磁场源.1985年, 南开大学丁守谦[4]发明了一种能产生均匀磁场的马鞍形线圈, 比同半径的赫姆霍兹线圈的均匀度高一个数量级以上.1986年, Turner[5]在设计磁共振成像设备的梯度线圈和屏蔽线圈时提出了目标场法.2000年, Moon和Hatano[6]采用流函数的方法设计凸表面梯度线圈.2002年, Forbes 和 Crozier[7]用采用目标场法设计了超导磁共振系统的圆柱状匀场线圈.2004年, Harvey 等[8]研究了流函数方法对线圈表面电流的优化实现匀场.2005年, 李霞和谢德馨[9]采用流函数方法设计了自屏蔽的单平面梯度线圈.2010年, Liu 等[10]采用目标场法设计了磁共振系统的双平面型匀场线圈.本文采用磁场球谐函数展开级数对主磁场进行数学分析, 选取磁场谐波分量的关键项, 结合目标场点法和流函数设计相应阶次的匀场线圈, 匀场线圈通电产生的磁场抵消主磁场相应的非均匀性分量, 从而实现磁体系统的匀场修正[11].

2 有源匀场理论推导

有源匀场线圈设计的实质是电磁场的计算问题, 电磁场的计算分为正问题和反问题两类, 正问题即已知场源(电流密度等)求解磁场分布, 反问题即已知磁场分布, 求解场源[12].本文的设计方法属于后者, 基于有源匀场连续电流密度分布的思想, 采用目标场点法和流函数结合的方法, 在柱坐标系下, 由毕奥-萨伐尔定律确定磁场分布与电流密度的关系, 约束线圈半径和设置约束点后, 采用目标场分布逆向求解线圈平面的电流密度分布, 再用流函数将密度分布离散化处理, 得到线圈绕线位置分布, 利用数值计算工具编译匀场线圈设计的算法, 从而实现高阶有源匀场线圈的设计.

2.1 基于目标场点法的磁场分析

磁共振系统的磁体空间是没有电荷分布的球形工作区域 (diameter of spherical volume, DSV),根据麦克斯韦方程组可知, 在无电流的情况下, 永磁磁体产生的无源静磁场是一个无旋场, 磁场的分布符合多重级数解析表达式, 与坐标无关的级数项是磁场的均匀分量, 与坐标成线性关系的级数项是一阶非均匀分量, 依此类推[13].

核磁共振分析仪的主磁场B0的方向定义为z轴方向, 射频磁场垂直于z轴方向施加, 设磁场中心为坐标原点, 在球坐标系(r,θ,ϕ)计算得到磁场某一点沿z方向的磁感应强度为

式中µ0为真空磁导率,为连带勒让德函数,为与r,θ,ϕ有关的常系数, 其中α,β分别为Bz的阶次和次序.(1)式说明磁场由谐波分量构成, 在直角坐标系(x,y,z)下表达式为[14−16]

若设计有源匀场线圈使之产生的磁场分别按(2)式中各项形式分布(第一项除外), 且与对应项的磁场大小相等、方向相反, 磁场叠加后相互抵消相应的非均匀性分量, 实现磁共振系统磁场的均匀性调节.例如:设置三个线性匀场线圈X,Y,Z, 消除 (2) 式中的三项, 磁场就可达到线性均匀.通常磁共振成像仪的梯度线圈产生的磁场是X,Y,Z的一阶函数, 常用作一阶有源匀场线圈,即梯度匀场.再增加能够消除等平方分量项的匀场线圈, 磁场就能达到平方均匀, 依此类推.

2.2 有源匀场电流密度函数分析

有源匀场后的磁场均匀度取决于所补偿的项数, 即匀场线圈的数目.匀场线圈的磁场特征是按空间分布, 实际匀场要调整各方向匀场线圈的电流来实现匀场线圈磁场大小的调节.根据(1)式可知,电流分布的影响因素由系数决定.假设线圈的半径为线圈的最大半径,极坐标系 (ρ,φ) 下线圈平面的电流函数的形式为S±(ρ,φ)满足连续流体方程, 则流函数为

表1 Bz 在直角坐标系和球坐标系下的分量表示Table 1.Component representation of Bz in Cartesian and spherical coordinates.

式中Uq为电流密度系数,q表示电流密度展开的项数 (1qQ, 为整数),Sq(ρ) 为流函数展开项.此处l,k分别代表线圈的阶次和次序.依据表1,Bz分量在直角坐标系和球坐标系下的分量表示中, 当k=0时, 电流密度将激发一个纵向(z方向)梯度场或均匀磁场; 当k=1 时, 线圈将激发在x或者y方向的横向线性梯度场; 当k⩾2 时, 即为高阶匀场线圈的情况; 对 于z,z2线 圈 来 说k=0 ; 对 于x,y,xz,yz线圈来说k=1 ; 对于xy线圈k=2.“±”代表线圈的上、下平面, 对于z、xz、yz、z3、z5线圈,上、下平面电流方向相反, 电流密度互为相反数;对于x,y,xyz2,x2+y2,x2−y2,z4线圈, 上、下平面电流方向相同, 电流密度也相同[18,19].

图1 目标场及流函数结合的有源匀场线圈设计算法流程Fig.1.Design algorithm flow of active shimming coil based on target field and flow function.

单一匀场线圈仅产生特定类型(正弦或者余弦) 的谐波场,l和k的值作为目标给出.(i) 为任意非零值, 与对应系数矩阵U成正比.当D为常数时, 求解对应系数矩阵U,U增大仅改变幅度,而不影响其分布.将求解的矩阵U代入(3)式获得流函数S(ρ,φ) , 对连续的流函数进行离散处理, 得到的电流路径即为线圈的布线轨迹.设流函数S的最大值Imax和最小值Imin, 在圆形平面的线圈匝数为N, 则每匝导线上的电流为

(6)式中线圈匝数N可调, 从而实现每匝线圈电流I0的可调, 根据N和I0的不同, 得到流函数的一组等值线:

电流密度函数的等值线图中等值线的位置, 即线圈绕线位置分布[20,21].

图2 各阶次线圈的流函数和绕线分布 (a) X, Y 线圈的流函数和绕线分布; (b) Z 线圈的流函数和绕线分布; (c) XY 线圈的流函数和绕线分布; (d) XZ, YZ 线圈的流函数和绕线分布; (e) Z2 线圈的流函数和绕线分布Fig.2.Flow function distribution and coil winding of each order coils:(a) Flow function distribution and coil winding of X and Y coil; (b) flow function distribution and coil winding of Z coil; (c) flow function distribution and coil winding of XY coil; (d) flow function distribution and coil winding of XZ and YZ coil; (e) flow function distribution and coil winding of Z2 coil.

3 有源匀场线圈仿真设计

3.1 有源匀场线圈设计算法

如图1所示, 目标场及流函数结合的有源匀场线圈设计算法流程.根据毕奥-萨伐尔定理, 采用目标场法和流函数结合的方式求解出线圈的绕线分布, 已知V,D的情况下, 关键在于如何通过方程求解出系数矩阵U, 方程求解的精度决定匀场线圈设计的成败[22].

在MATLAB平台下编译了相应的算法, 算法设计步骤如下:1)确定待设计的匀场线圈的阶次和序次, 确定目标场的阶次和序次; 2)根据核磁共振分析仪磁体尺寸给定双平面型有源匀场线圈圆盘尺寸、线圈位置约束、球形工作区域以及目标场点; 3)在柱坐标下给出有源匀场线圈电流密度、磁感应强度表达式; 4)求解不适定反问题方程得出电流密度表达式系数; 5)改变目标场点, 增加展开式级数, 转 Step5 继续计算; 6)根据计算出的电流密度表达式系数, 得到匀场线圈的电流密度表达式; 7)采用流函数技术进行离散化处理, 输入设定的绕线匝数; 8)终止计算, 输出绕线形状[23].

3.2 有源匀场线圈设计结果

图3 采用有源匀场线圈匀场前后磁共振检测的FID信号和频谱的FWHM (a)有源匀场线圈应用于磁共振系统的实验平台;(b)匀场前FID信号; (c)匀场前频谱的FWHM; (d)二阶匀场线圈匀场后的FID信号; (e)二阶匀场线圈匀场后频谱的FWHMFig.3.FID signal detected by magnetic resonance before and after shimming with active shimming coil:(a) Experimental platform of active shimming coil applied to magnetic resonance system; (b) FID signal before shimming; (c) FWHM of the pre-shimming spectrum; (d) FID signal after the second-order shimming coil shimming; (e) FWHM of after the second-order shimming coil shimming.

表2 有源匀场线圈匀场效果评价技术指标对比Table 2.Comparison of technical indicators for evaluating shimming effect of active shimming coil.

为验证目标场及流函数结合的有源匀场线圈设计算法的有效性, 选取磁体的磁极间距为42 mm的核磁共振分析仪有源匀场线圈作为设计对象, 其工作区直径为27.5 mm球形工作区域, 设计3组有源匀场一阶线圈(X,Y,Z)以及4组二阶线圈(XY,XZ,YZ,Z2), 通过有源匀场线圈产生的磁场对主磁场进行补偿.通过MATLAB实现有源匀场线圈设计仿真, 运行程序后输出7组匀场线圈的绕线分布, 如图2, 各阶次线圈的流函数和绕线分布[24].

4 实验讨论

根据磁共振分析仪主磁场均匀性的调节原理,主磁场越均匀, 磁共振的自由感应衰减信号(free induction decay, FID)的拖尾越长, 信号的积分面积越大, 信号的频谱越窄, 即信号频谱的半高宽(full width half maximum, FWHM)越小, 信号的分辨率越高.如图3所示, 采用有源匀场线圈匀场前后磁共振检测的FID信号和频谱的FWHM, 设计实验如下:采用长的横向弛豫时间(transverse relaxation time, T2)检测蒸馏水样品, 设置射频中心频率设置为17.02 MHz, 将通过算法仿真设计制作出的有源匀场线圈实物安装在核磁共振分析仪磁体两极上, 选择硬脉进行磁共振分析仪的FID信号采集.

有源匀场线圈匀场前后的FID信号积分面积、频谱半高宽、磁场均匀性等技术指标如表2所列, 有源匀场线圈匀场效果技术评价指标对比, 分析可知, 有源匀场线圈匀场后FID信号的积分面积显著增大, 频谱的半高宽明显减小, 说明磁场均匀性得到了很大改善[25,26].

5 结 论

本文根据有源匀场线圈目标场点法和流函数结合的设计理论, 详细推导了磁场的球谐函数展开级数, 对磁共振系统的磁体产生的主磁场进行数学分析, 确定选取一阶、二阶及关键高阶次项进行有源匀场线圈的设计, 利用匀场线圈产生的磁场来补偿主磁场中的非均匀性磁场分量.系统阐述了有源匀场线圈设计算法, 并利用数值计算工具编译实现了软件算法, 通过输入目标参数, 自动输出有源匀场线圈的绕线的形状数据, 根据各阶次线圈的流函数分布及线圈的绕线生产制作了匀场线圈, 应用于核磁共振分析仪的磁体系统, 根据硬脉冲FID信号采集实验结果, 有源匀场线圈匀场后的FID信号积分面积显著增大, 频谱的半高宽明显减小, 说明基于目标场及流函数结合设计制作的有源匀场线圈能够有效补偿主磁场中的不均匀磁场, 使磁共振系统的磁场均匀性得到了很大改善, 能够达到匀场的目的.