脑水肿检测的开放式核磁共振系统构建,神经病学论文

脑水肿检测的开放式核磁共振系统构建,神经病学论文开放式核磁共振〔opennuclearmagneticreson-ance,ONMR〕系统具有由永磁磁体构成主磁场的磁体、系统体积小、质量小、共振区域位于磁体一侧等特点[1],具备良好的开放性，可用于对材料的分析和检测[2].赵薇等设计了1种由鞍型线圈组成、能够对Halbach磁体的横向磁场进行补偿的柱面匀场线圈[3].胡格丽等采用目的场方式方法设计了1组超导核磁共振成像系统梯度线圈[4].典型的主磁体按构造分为U型磁体和桶状磁体[5]以及Halbach磁体[6-7]等。含水量是材料的一个重要参数。大型油浸变压器中绝缘油和绝缘纸的含水量对油纸绝缘性能的影响极为显着[8].王有元等进行了水分对绝缘纸微观特性影响的分子模拟研究[9]和水分对油纸绝缘热老化速率的影响[10].许渊等采用介质响应技术检测油纸绝缘的含水量[11],张致等采用频域谱方式方法估计固体绝缘的含水量[12],为绝缘诊断提供参考信息。水中的1H原子核是核磁共振〔NMR〕经常测定的磁性原子核之一。核磁共振能够测定冰冻枯燥人红细胞的含水量[13]、石油的含水量[14]、茶叶的含水量[15]和地下水分布[16].颅内含水体积的变化情况能够直接反映脑水肿病情的严重程度。在脑猝死〔脑梗死和脑出血〕、脑炎、颅内肿瘤、脑损伤等病情的同时，通常都伴随着脑水肿[17],颅内含水体积的变化与病情的变化息息相关。针对计算机断层扫描〔CT〕所用的X射线会对人体造成伤害、对软组织分辨率不高和磁共振成像〔MRI〕仪器价格昂贵、体积大等问题，本文设计1种小型化的、价格低廉的、开放式核磁共振分析仪，通过检测目的区域不同深度物体的核磁共振信号和被测物在不同含水体积情况下的核磁共振信号，来检测头颅中不同位置的含水体积。核磁共振系统的关键技术在于主磁体设计以及线圈的优化选择。本文首先分析主磁体构造，使其磁场的梯度、均匀度以及精度知足要求；接着研究不同类型射频〔RF〕线圈的磁场分布规律，挑选适宜搭建实验系统的螺线型射频线圈，使其磁场的平行度和均匀度知足要求；最后搭建核磁共振系统，以位于磁体目的区域不同深度位置的橡皮和不同含水体积的面粉为样品，用实验验证系统的可靠性和准确性，为开放式核磁共振仪进一步应用于临床提供实验支撑。1核磁共振仪主体技术1.1主磁体为了对颅内疾病进行实时监护，应使磁体构造能够包围人体头颅，并在磁体构造的中心区域构成水平方向平坦分布、纵深方向梯度分布的磁场，便于对颅内不同深度区域的含水体积信息进行检测。本文设计得到了磁场条件优良的主磁体。向永磁体制造厂商〔北京中科三环高技术股份有限公司〕购买所需磁体，即磁块尺寸为30mm30mm100mm、相对磁导率r为1.05、剩磁的磁感应强度Br为1.38T的钕铁硼磁块，并采用不导磁的不锈钢材料制作支架，构建实际的磁体模型，其主磁体构造如此图1所示。磁体磁场中心区域的尺寸为50mm50mm50mm.磁场中心区域yoz平面50mm50mm范围磁场磁感应强度Bz的计算值和测量值如此图2所示。从图2能够看出，在yoz平面上的范围为50mm50mm时，即磁场中心区域在y方向的宽度为50mm时，磁场沿水平方向平坦分布。主磁体在磁场中心区域内磁场磁感应强度幅值的测量值为0.063~0.1647T〔计算值为0.063~0.1635T〕，平均值的测量值为0.1077T〔计算值为0.1071T〕，磁场梯度的测量值为2.03T/m〔计算值为2.01T/m〕，磁场梯度较大，可为核磁共振系统提供良好的主磁场。1.2射频线圈射频线圈是装在主磁体内孔里的常规高频脉冲线圈，射频线圈由晶体振荡器和功率放大器提供射频电流工作，核磁共振信号也由同一个射频线圈接收，接收到的信号经过接收器和放大器进行模拟和数字转换、信号平均，最后计算、存储和显示完成样品的成像[18].主磁体的磁场方向必须跟射频线圈的磁场方向垂直，被测物体处于主磁体磁场和射频线圈磁场的共同磁场当中，因而射频线圈的设计不能脱离主磁体磁场，需要根据主磁体磁场的特点进行设计，并搭建核磁共振实验平台验证其可用性和匹配性[19].1.2.1不同射频线圈磁场分析为了选取合适该系统的射频线圈，对笼型线圈、单层八字线圈、双层八字线圈以及螺线型线圈的磁场大小、方向及分布进行了研究，以产生磁场的磁感应强度、均匀度以及平行度为标准进行挑选，最后选定螺线型线圈搭建实验系统。采用外电路添加鼓励的方式进行建模，线圈宽度为0.5mm,电流鼓励为2A,匝间距为0.5mm,线圈外半径为26mm,两线圈之间距离为70mm.螺线型线圈构造及磁感应强度B的大小和方向如此图3所示。从图3能够看出，螺线型线圈产生的磁场均匀度和平行度都较好，可用来设计核磁共振系统。1.2.2螺线型线圈的优化选择利用Maxwell对线圈进行建模拟真，建模需要消耗损费大量的时间，仿真时间也较长，对计算机的要求较高，而影响螺线型线圈磁场的因素较多，通过仿真建模来逐个研究影响因素是不切实际的。本文考虑到这些因素，通过算法编程，设置线圈内半径、线圈外半径、线宽、匝间距、电流大小等参数建立了螺线型线圈的磁场计算方式方法，研究了不同影响因素对磁场的影响规律。研究得到，随着匝间距的增加，磁感应强度变小，磁场均匀度不变；匝间距在0.1~0.5mm范围内时线圈磁场磁感应强度变化不大，匝间距大于0.5mm时，磁场磁感应强度变小；随着2个线圈之间的距离变大，线圈中间位置处的磁感应强度变小，磁场均匀度变化不大；在其他参数固定的条件下，线圈存在使得磁场磁感应强度到达最大值的最优外半径。根据参数优化，设线圈匝数为7匝，线圈外半径为40mm,2个线圈之间的距离为40mm,匝间距为0.3mm,此时磁场磁感应强度到达1mT,其均匀度为0.02mT,优化选择参数后的磁场磁感应强度及其均匀度都有了明显的改善，此时螺线型线圈在中间截面处只要z方向的值。2核磁共振实验2.1实验系统框架与实验条件本文搭建的核磁共振系统构造如此图4所示。在搭建核磁共振实验平台时，射频线圈通过电容匹配电路〔Cm、CT为匹配电容〕与电阻值50?的同轴电缆相连，使得射频线圈的阻抗与实验仪器阻抗相匹配，再将同轴电缆与实验仪器相连，通过操作Prospa软件进行实验。实现核磁共振要知足的基本条件如下：1〕核有自旋〔磁性核〕；2〕有1个稳恒的外磁场B0,能发生能级分裂。主磁场为均匀场，磁场越均匀、Zeeman能级宽度越窄、共振吸收峰越尖锐、信噪比越高都能使得信号越容易观测；3〕有1个与主磁场和总磁矩所组成平面垂直的旋转磁场B1;4〕若线圈磁场角频率00==B〔0B为矢量B0的模〕，则发生核磁共振。华而不实为核的磁旋比，其计算式为【公式】2.2线圈PCB模型根据前面确定的线圈参数，画出螺线型线圈的印刷电路板〔printedcircuitboard,PCB〕模型，如此图5所示。图5中：2个线圈中的A、A、B、B、C、C、D、D为螺杆连接点，都为直径4.5mm的圆洞；AA、BB、CC、DD分别相连接用来固定板子；E、F、E、F为焊接点。在实验中，将2个线圈用塑料螺钉固定，通过焊接点连接。E连接F,E、F与外电路相连。2.3磁体磁场测量与数据计算3维磁场测量平台由1维磁场探针、1维贝尔高斯计、3维运动平台〔由x、y、z这3个方向的高精度步进电机控制〕、运动控制器以及装有控制测量软件的计算机构成，如此图6所示。图6中：磁场探针与高斯计相连，高斯计通过串口与计算机相连，磁场数据由探针测得，共同显示在高斯计和测量软件中；磁场探针固定在3维运动平台上，运动平台由运动控制器控制，通过串口连接在计算机上，由测量软件控制可在1维直线和2维平面上运动；这样通过2个串口便可在终端软件上控制磁场探针运动，并记录测量数据。分别测出在磁场中心区域？25mm、？15mm、?5mm、5mm、15mm和25mm这6个高度上的磁场磁感应强度大小，然后计算得到各位置处的共振频率，再利用阻抗分析仪测得线圈在每个共振频率下的电阻R值、电感L值。将线圈放置在主磁体当中，进行数据测量，每个共振频率位置处测量4组数据，最后求取平均值。测得线圈的阻抗参数，利用算法计算电容匹配电路中电容Cm、CT的值。在磁场中心区域中，不同坐标位置的层面上，主磁体磁场中心区域不同位置高度h处的各种数据〔磁场磁感应强度平均值B、共振频率f、电阻R、电感L、匹配电容Cm、CT〕如表1所示。2.4电路匹配与频率调节根据计算出的匹配电容，完成电容匹配电路的计算，根据前面图4中型等效电路来完成电容匹配电路的焊接。根据计算的参数来完成电路匹配，完成匹配后利用频率分析仪验证匹配的好坏。利用频率分析仪来测量所匹配的电路所对应的共振频率，若此时对应的共振频率的测量值大于计算值，则讲明选择的电容CT过小，所以需要将CT值变大；反之则将CT变小。如此反复调节，直到共振频率的测量值与计算值一致。频率匹配如此图7所示。图7中：f1、f2为幅值下降到总幅值d长度的0.707倍时所对应的频率；Q=f/f为表征线圈匹配好坏的参数；f为共振频率；f=f2?f1为频率差值；Q值越大，讲明线圈的品质因数越好。在共振频率〔计算理论值〕为4.431MHz时匹配电路，在仪器测量中匹配频率到达了4.4187MHz,精度到达了99.71%,f2=4.30901MHz,f1=4.21386MHz,f=0.09515MHz,因而计算得到Q=4.4187MHz/0.09515MHz=46.44,线圈品质因数高。2.5核磁共振系统的建立图8为利用Magritek核磁共振谱仪、Tomco射频功率放大器、半环Halbach磁体、射频线圈以及计算机搭建的开放式核磁共振系统。使用的脉冲信号为自旋回波序列〔carr-purcell-meiboom-gill,CPMG〕射频脉冲序列，为用于测横向弛豫时间2的常用脉冲序列，由90和180的射频脉冲组成,其基本单元90-delay-〔180-delay〕n〔n为回波个数〕如此图9所示[22].a1表示90脉冲的幅值衰减倍数，a2表示180脉冲的幅值衰减倍数，td1表示90和180脉冲的宽度，te表示回波时间，tr表示2个连续的90射频脉冲间的时间间隔。图9中，在90脉冲作用后，每隔回波时间te采用180脉冲重聚焦，并测量其回波信号，理论上回波信号峰值的连线将以横向驰豫时间2为衰减时间常数作指数衰减，可用下式表示为[13]【公式2】式中：y0为测量信号的初始幅值；a为衰减系数；y为随着时间变化的测量信号的幅值。2.5.1目的区域不同深度物体回波信号测量CPMG脉冲序列中90和180脉冲宽度td1为4s,90和180脉冲幅值衰减倍数a1和a2分别为?18dB和？12dB,回波时间te为110s,回波个数n为60,重复次数为2000次，重复时间tr为300ms.将橡皮分别放置在磁场中心区域？25mm〔7.213MHz〕、？15mm〔6.047MHz〕、？5mm〔5.158MHz〕、5mm〔4.168MHz〕、15mm〔3.515MHz〕和25mm〔2.782MHz〕这6个高度上〔射频线圈随样品一起放置在不同高度上〕，单个高度上测得的核磁共振信号如此图10所示。不同高度上橡皮的核磁共振信号如此图11所示，将每个高度上的回波信号叠加。随着纵深距离增加，叠加构成信号的幅值随高度而逐步衰减，原由于主磁体的磁感应强度在纵深方向呈梯度衰减趋势。本实验证明，所设计的开放式核磁共振系统能够获得目的区域不同深度物体的核磁共振信号，通过测定核磁共振信号的峰值频率能够断定样品的放置位置深度。2.5.2目的区域物体不同含水体积信号的测量脑水肿是脑内水增加导致脑容积增大的病理现象，颅内含水体积的变化能够直接反映出脑水肿病情的变化情况。物体含水体积不同，所测得的核磁共振信号也不同，这种不同能够直接通过2的大小来表示。在实验中，用1个容器内部填充可搅拌物质来模拟脑内情况。被测物体为1个塑料容器，华而不实加面粉，通过参加定量清水来改变其含水体积。将被测物体放置在共振频率计算值为4.431MHz位置处，测定2的大小。测得该位置处4组R、L值，其平均值分别为1.54556?、4.53082H,计算得到此时的匹配电容2Cm、CT值分别为99.5742pF和233.6086pF.电容选择稳定性好、绝缘性强、耐压高的陶瓷电容。通过选择不同大小的电容进行组合、电路匹配和频率测量，选取最接近计算频率的电容值。最后选取2Cm、CT的值分别为100pF、〔220+10+3.3〕pF,用Wobble测得匹配的频率大小为4.453MHz,误差仅为0.49%.搭建好整个系统，在核磁共振实验Prospa软件中进行最后Wobble频率匹配。CT值影响匹配频率值的大小，Cm值影响匹配频率的品质因数。图12表示频率匹配好坏的频率扫描图，频率扫描响应图中越细越尖讲明匹配的效果越好，共振吸收峰越尖锐，信噪比越高，越容易观测。从图12能够看出，匹配的效果非常好，中心频率参数设置为4.453MHz.另1个重要参数是脉冲宽度，通过进行脉冲宽度扫描，选择使得信号最大的脉冲宽度值。图13表示不同脉宽下核磁共振信号幅值的大小。从图13能够看到，在脉冲宽度为7s时得到最大核磁共振信号。其他参数通过不断改变调试来选取，选取CpmgAdd信号〔核磁共振信号〕最好的参数值。实验选取CPMG脉冲序列中90和180脉冲宽度td1为4s,幅值衰减倍数a1和a2分别为？21dB和？15dB,回波时间te为200s,回波个数n为4000,重复次数为256次，重复时间tr为5000ms.将被测物体放入系统中间位置，第1组数据为面粉中参加体积为1mL的水，从第2组开场往后每1组参加体积为0.2mL的水，即每1组的含水体积依次为1mL、1.2mL、1.4mL,依此类推，并搅拌均匀，进行实验。图14为参加水体积为2.4mL时的拟合结果，用t/ms表示后1组相对于前1组数据的增量，以此来表征信号的变化。图14〔a〕表示未发生核磁共振时的曲线图〔2=0.1ms〕；图14〔b〕表示发生核磁共振时最后拟合成的曲线图〔2=96.9ms〕。测得不同含水体积情况下表征核磁共振信号大小的2值，2为横向弛豫时间，该值能够表征被测物体内部含有的H1的数量，H1数量越多，测得的核磁共振信号越大，2也会越大，2是核磁共振中一个很重要的参数，如表2所示。从表2中可出看出，随着每次增加体积为0.2mL的水，2的值变大，t为后1组测得的2与前1组测得的2的差值，t的变化范围都在5.5~7.5ms之间。在实际实验中，这个值有一定误差，由于在实验中，随着含水体积的增加