CN107064838B - 一种能够形成变梯度静磁场的磁体系统结构及测量方法 - Google Patents
一种能够形成变梯度静磁场的磁体系统结构及测量方法 Download PDFInfo
- Publication number
- CN107064838B CN107064838B CN201710254930.3A CN201710254930A CN107064838B CN 107064838 B CN107064838 B CN 107064838B CN 201710254930 A CN201710254930 A CN 201710254930A CN 107064838 B CN107064838 B CN 107064838B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- static magnetic
- gradient
- magnet
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims description 7
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 24
- 238000003384 imaging method Methods 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 6
- 229910052704 radon Inorganic materials 0.000 claims abstract description 6
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 6
- 238000009792 diffusion process Methods 0.000 claims description 35
- 230000005415 magnetization Effects 0.000 claims description 32
- 239000013598 vector Substances 0.000 claims description 26
- 239000002184 metal Substances 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 6
- 238000000691 measurement method Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 230000001105 regulatory effect Effects 0.000 claims description 5
- 238000003491 array Methods 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 4
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000036962 time dependent Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 32
- 238000005259 measurement Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 7
- 230000009977 dual effect Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000264 spin echo pulse sequence Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供一种能够形成变梯度静磁场的磁体系统结构,适用于核磁共振检测装置系统。该系统由磁体部分,固定磁体的外壳部分和磁场调节部分组成。本发明提供的形成变梯度静磁场的磁体系统所产生的变梯度磁场的梯度方向和大小在均能调整,在该环境下对样品进行测试后对数据进行拉东变换可作为一种磁共振的成像方式。
Description
技术领域
本发明阐述一种能够在腔体内部产生自带恒定静磁场梯度的磁体系统。通过旋转该磁体系统内环或外环磁体子系统,可调整静磁场梯度在空间中的方向。该磁体系统可用于获取核磁共振中对被测样品的分子自扩散系数测量以及被测样品的二维成像结果。
背景技术
核磁共振检测技术是利用核磁共振原理探测氢原子的技术。通过探测被测物内氢原子的含量和赋存状态,获得被测物内各种成分的信息。其技术手段主要通过核磁共振探测装置的探头内的磁体形成静磁场,通过射频天线向被测物发射射频磁场脉冲、采集共振信号,进而根据采集到的样品不同位置的信号,通过直接测量被测物体中氢核的密度、流体分子弛豫特性、扩散特性等对贮存在被测物内的各成分进行分析。
实现核磁共振测量的必要条件是将被测样品放置在磁场中。磁场的产生可通过永磁体实现。J.A.Jackson等提出将两个轴向磁化的圆柱形磁体北极与北极相对放置,在其外侧建立起相对均匀的圆环状分布静态区域,在其外侧建立起相对均匀的圆环状分布静态区域,相应的射频磁场由位于两磁体中间的螺旋线圈产生,虽然该磁体结构简单,但灵敏区域太小且磁场强度较低。随后,Z.Taicher等使用径向磁化的单个圆柱形磁体建立静态磁场,而射频磁场由缠绕在磁体上的螺旋线圈产生,但其静态磁场分布极为不均,该传感器的灵敏度区域很小,随后,O.Sucre等为了进一步提高灵敏区域内静态磁场的均匀度,将磁体结构改为六个径向磁化的圆柱形磁体,平面矩形线圈被置于传感器中部的磁体表面上。
由此可见,磁体结构的紧凑性设计一直是研究学者们不断追求的目标,都是向传感器的静态磁场的均匀度与灵敏度之间寻找一个优化方案。
然而现有的核磁共振传感器的设计存在螺旋线圈体积过大,而且在被测区域内无法产生静磁场梯度。静磁场梯度对于核磁共振测量分子自扩散系数或实现磁共振成像的空间编码至关重要。因此,缺失静磁场梯度的磁体系统无法对大量存在的各向异性的样品进行测量和研究。
发明内容
针对现有技术的上述缺陷,本发明提供一种能够形成变梯度静磁场的磁体系统结构,具有结构紧密,形成的磁场稳定,磁场梯度的强度和方向可调节,精度高等特点。本发明所要解决的实际问题在于,首先可在磁体系统内部形成静磁场梯度,其次此磁体系统所产生的梯度在空间方向可调整。鉴于上述两点优势,基于此磁体系统组建完成的核磁共振分析仪器可实现对被测样品的自扩散系数测量分析,得到内部流体种类、含量及其各自组分等重要信息。同时,通过改变梯的空间方向,测量得到一系列被测样品的一维成像剖面,进一步对得到的数据进行拉东变换后可实现对被测样品的二维平面核磁共振成像。该实现方法无需额外搭建梯度成像线圈,因此可大大的降低多维核磁共振成像技术在样品分析中的检测成本。
为了实现上述目的,本发明采用以下技术方案:
一种能够形成变梯度静磁场的磁体系统结构,适用于核磁共振检测装置系统,无需额外搭建梯度成像线圈;
进一步地,所述能够形成变梯度静磁场的磁体系统结构由磁体部分,固定磁体的外壳部分和磁场调节部分组成,其中磁体部分由两层环形的磁体阵列组成,分别为内环磁体阵列和外环磁体阵列;所述内环磁体阵列和外环磁体阵列均由沿着柱体轴向方向的n层磁环排列组成,所述n层磁环由2N个磁棒构成(n和N为大于2的自然数),阵列组成存在2N-1对耦合级磁棒,每一对耦合级磁棒按照沿柱体轴对称方式排布,每一层磁环的阵列组成均可相同或不同:所述内环磁体阵列的耦合级磁棒的磁化方向相同,所述外环磁体阵列的耦合级磁棒的磁化方向相反。
进一步地,每一层磁环的阵列组成沿着圆周方向排列,每层磁环中磁棒的厚度和层与层之间的距离,以及每一层磁环的阵列组成的圆周的直径均由迭代算法优化设计得出,所述内环磁体阵列的磁棒的圆周排列方式形成的腔内产生第一静磁场,所述第一静磁场具有特定磁场强度,所述第一静磁场方向与邻近的磁棒磁化方向相同。
进一步地,所述外环磁体阵列的磁棒的圆周排列方式形成的腔内产生第二静磁场,所述第二静磁场具有恒定梯度,所述第一静磁场和所述第二静磁场形成变梯度静磁场的所述磁体系统结构,所述变梯度静磁场通过改变内环磁体阵列与所述外环磁体阵列的相对角度而改变静磁场梯度的空间取向;所述内环磁体阵列与所述外环磁体阵列由各自的内侧金属骨架和外侧金属骨架安装固定。
进一步地,固定磁体的外壳部分由所述金属骨架和温控部分组成,所述温控部分安装于所述内侧金属骨架和外侧金属骨架之间,用于保证所述内环磁体阵列与所述外环磁体阵列的相对温度不变,由高精度温度传感器和温度调节系统组成,能够有效防止磁体的温度漂移导致的磁场变化;调节静磁场梯度的空间取向通过旋转每环磁体阵列的外壳实现自由转动,所述外壳下部有固定孔,用于连接磁场调节部分。
进一步地,所述磁场调节部分由底盘和无磁电机组成,底盘为两层环状结构,分别对应内环磁体阵列与所述外环磁体阵列的外壳部分,无磁电机带动底盘转动,所述磁场调节部分由控制器操作,控制器根据指令动作控制无磁电机旋转,从而能够使所述磁体部分之间相互转动。
本发明还提供了一种能够形成变梯度静磁场的磁体系统结构的测量方法,所述测量方法包括:取N=4,n=4,初始化系统和恒温温度,调整无磁电机至初始位置,此时标记静磁场梯度G的空间取向φ=0,将被测样品放置于所述磁体系统结构的测试腔内,在该梯度方向下进行一维投影,首先对被测样品施加90°射频脉冲,将磁化矢量M0扳转至垂直于静磁场方向的横向平面上,M0由静磁场强度B0、温度等参数决定;由于分子的扩散及静磁场的空间非均匀性等原因,磁化矢量M0发生散相;经历一定时间τ之后,施加180°脉冲;散相后的磁化矢量会在同等时间τ之后实现重聚,形成一个回波信号,然后进行以下步骤:(1)在梯度磁场下通过改变梯度幅值或梯度持续时间,记录回波幅值的变化,得到流体分子的自扩散系数;(2)通过施加成对的频率编码或相位编码梯度,解析被测样品的空间自旋密度信息,实现核磁共振成像。
进一步地,实现对在静磁场中形成的所述磁化矢量M0的扳转操作由射频磁场B1完成,扳转角度为:θ=γB1tp,其中B1为射频磁场强度,tp为射频脉冲的持续时间,γ为质子的旋磁比;通过控制射频脉冲的幅值或持续时间达到改变扳转角的目的。
进一步地,所述步骤(1)中由于存在静磁场非均匀性及分子的扩散时,对横向磁化矢量造成的额外衰减,进一步考虑分子的扩散造成的磁化矢量衰减,在宏观角度上来说,该响应是一个与时间相关的函数,因此在考虑分子扩散的影响后,磁化矢量衰减满足以下公式:其中为静磁场B0的非均匀性,或梯度,
γ/2π=42.58MHz/T,
δ为半回波间隔。
进一步地,所述步骤(2)中通过一种无损的方式快速反映被测样品内部流体的空间分布情况,由于梯度的存在,首先可以获得被测样品某一方向上的一维投影剖面,一维投影剖面p函数是通过对得到的自旋回波信号进行傅立叶变换得到的其中,G为磁场梯度,n为采集自旋回波信号中的点数,td为采集回波中点与点的时间间隔,为每个时间点td下的磁化矢量大小,即核磁共振信号幅度值,Z代表成像位置轴,因此,通过对得到的单个回波整体信号进行傅立叶变换,即可得到一维核磁共振的成像剖面p函数,通过所述磁场调节部分,改变梯度静磁场与均匀静磁场之间的夹角φ,得到一系列的投影剖面p函数,通过反拉东变换,既可以重构得到被测样品的二维平面图像
I(x,y)=∫p(xcosφ+ysinφ,φ)dφ
其中,I(x,y)为二维图像结果,φ为静磁场梯度G与静磁场磁感应强度B0方向的夹角,p函数为在φ夹角下的一维磁共振成像剖面。
固定磁体的外壳部分由金属骨架和温控装置共同组成,金属骨架部分用于固定各环磁体;温控部分用于控制本发明中使用的永磁体的相对温度不变,因此保证核磁共振测量过程中的拉莫尔频率准确。该部分由高精度温度传感器和温度调节系统组成,能够有效监测磁体的温度漂移导致的磁场变化。固定每环磁体的外壳之间均相互独立,可将两个磁体阵列组装在一起形成最终的磁体系统。这样,在内环磁体阵列中产生的均匀静磁场和在外环磁体阵列中产生的自带梯度静磁场叠加后,最终在整体磁体系统内形成具有校正后的梯度静磁场分布。在目前磁体组装状态下,由于在测量区域内部存在梯度静磁场,可对被测样品施加自旋回波脉冲序列,通过改变回波间隔,完成对被测量样品的分子自扩散系数测量。关于自扩散系数的测量本专利将在随后详细介绍,此处不再赘述。
除了完成上述组转后,该系统还可以对内环或外环进行磁场旋转操作,实现梯度方向在空间的变化。磁场调节部分由底盘和无磁电机组成,底盘为双层环状结构,对应每环磁体的外壳部分,相互连接后由无磁电机带动底盘转动,从而能够使各环磁体之间相互转动,调节静磁场梯度方向。通过施加一维成像脉冲序列实现对被测样品在不同投影取向上的一维剖面成像。对得到上述一系列一维结果进行反拉东变化,进而可获得被测样品的二维核磁共振成像分析。
附图说明
图1为本发明提供的能够形成变梯度静磁场的磁体系统结构的整体一双环结构实施例示意图;
图2为本发明提供的双环结构实施例的磁体结构示意图;
图3为本发明提供的双环结构实施例的内环磁体结构和充磁置示意图;
图4为本发明提供的双环结构实施例的内环磁体固定外壳部分示意图;
图5为本发明提供的双环结构实施例的外环磁体结构和充磁置示意图;
图6为本发明提供的双环结构实施例的外环磁体固定外壳部分示意图;
图7为本发明提供的双环结构实施例的两环磁体外壳部分组合的示意图;
图8为本发明提供的双环结构实施例的磁场调节部分结构示意图;
图9为本发明提供的双环结构实施例的磁场调节旋转方式示意图;
图10为本发明提供的双环结构实施例的磁场调节旋转具体情况示意图;
图11为本发明磁场环境下核磁共振测量方法脉冲序列;
实施方式
结合说明书附图说明本发明的具体实施方式。在此需要说明的是,对于这些实施例方式的说明用于帮助理解本发明,但并不构成对本发明的限定。
图1显示说明了本专利一种能够形成变梯度静磁场的双环磁体系统结构,具有结构紧密,形成的磁场稳定,磁场梯度的方向可调节,精度高等特点。该实施例中整个系统由磁体部分1,固定磁体的外壳部分2和磁场调节部分3组成。
本实施例中采用磁体结构为双层环结构,如图2所示,分别为外层磁体阵列4和内层磁体阵列5。考虑到安装便易程度和磁场强度,每个环的层面上按照圆环带有限剖分思路均安置16块磁体,每个永磁体均为八棱柱,两层磁体阵列的组成方式相同。
本实施例中内层磁体阵列5的结构和充磁方式如图3所示。内层磁体阵列5由4个磁体层6,7,8,9组成,其永磁体长度分别为l1,l2,l3和l4。其相互之间的空间间隔为h1,h2,h3。将每个磁棒考虑为类圆横截面情况下的直径为d。通过采用特定的迭代优化算法,对上述参数l1,l2,l3,l4,h1,h2,h3以及d进行优化计算,保证能够在该磁体阵列的磁体腔内产生具有一定均匀度的均匀静磁场B0 h。该均匀静磁场将成为实现本专利所述梯度静磁场的基础和先决条件。最终组装完成后采用具有恒温效果的外壳10对该内磁体阵列5进行封装,如图4所示。
本实施例中外层磁体阵列4的结构和充磁方式如图5所示。外层磁体阵列4同样由4个磁体层l1,l2,l3,l4组成,其永磁体长度分别为L1,L2,L3和L4。其相互之间的空间间隔为H1,H2,H3。将每个磁棒考虑为类圆横截面情况下的直径为D。每个磁体层中各个永磁体块的充磁方向为辐射式充磁。通过采用特定的迭代优化算法,对上述参数L1,L2,L3,L4,H1,H2,H3以及D进行优化计算,保证能够在该磁体阵列的磁体腔内产生具有恒定梯度G的静磁场。该磁体阵列将成为实现本专利所述梯度静磁场的核心。最终组装完成后采用具有恒温效果的外壳15对该外磁体阵列4进行封装,如图6所示。
本实施例中的两环磁体外壳组装的示意图如图7所示。外壳10和15由金属骨架和温控装置共同组成,金属骨架部分用于固定各环磁体,温控部分用于保证磁体的相对温度不变,由高精度温度传感器和温度调节系统组成,位于两环磁体的外壳内部,能够有效放置磁体的温度漂移导致的磁场变化。将两环带有各自磁体系统的磁体外壳套接在一起,即可在最终系统的腔体内部形成具有特定磁场强度B0 h和恒定梯度G的静磁场B0。该磁场可完成对被测样品中分子自扩散系数的测量,得到被测样品内部饱含流体的种类、含量和流体精细组份等重要信息。
由图7可知,固定每环磁体的外壳10和15之间均相互独立,可组装在一起后自由转动,其下部有固定孔,用于连接磁场调节部分,从而能够自由调节静磁场的磁场梯度G的空间方向。
本实施例中磁场调节部分如图8所示,由底盘和无磁电机组成。底盘为双层环状结构,分别包含16和17部分,对应各环磁体的外壳,相互连接后由无磁电机带动底盘转动,从而能够使各环磁体之间相互转动,达到调节静磁场的磁场梯度的效果。本实施例中采用的方式为两环磁体的相对转动,如图9所示。通过旋转底盘16或17,带动磁体阵列4或5的相对空间位置转动,进而改变静磁场梯度G的空间取向。图10显示了固定外层磁体阵列4,转动内层磁体阵列5后磁场梯度G的方向φ的三种情况。每一个φ角度下通过核磁共振成像可以得到被测样品在该梯度方向下的一维投影图像。通过多次改变角度φ并对得到的多组一维投影图像进行数据重构,即可得到被测样品在横向平面上的二维核磁共振图像。
静磁场B0在本专利中由永磁体提供,该值大小决定核磁共振信号的信噪比。被测样品置于静磁场中,自旋系统内发生能级分裂,沿着静磁场方向会产生一个宏观磁化矢量M0。M0由静磁场强度B0,温度等参数决定。因此为了保证核磁共振测量稳定性,需要对永磁体系统设定恒温装置,保证外界环境不会对永磁体的特性产生影响。
射频磁场B1产生的射频脉冲为电磁信号,通常由线圈产生。射频脉冲产生的磁场为射频磁场。射频磁场的方向与静磁场方向垂直,实现对在静磁场中形成的磁化矢量的扳转操作,扳转角度为:θ=γB1tp。其中B1为射频磁场强度,tp为射频脉冲的持续时间。因此可通过控制射频脉冲的幅值或持续时间达到改变扳转角的目的。核磁共振脉冲序列由不同数量和频率属性的射频脉冲按照设定时序组成。通过调整脉冲间时间间隔,脉冲角度及脉冲的频率选择性,实现对被测样品的弛豫、扩散系数和成像等测量。
磁场梯度分为静磁场梯度和脉冲磁场梯度两种。本发明专利只涉及静磁场梯度。通过记录一定时间内分子沿着梯度方向平均扩散位移即可计算出分子的自扩散系数。此方法作为一种有效的自扩散系数测量,应用于流体种类识别及样品标定等领域。同时,梯度静磁场的存在也可以完成对被测样品中流体信号的空间编码,进而实现核磁共振成像测量。
自旋回波为核磁共振测量最常见的一种信号。脉冲序列如图11所示。首先对被测样品施加90°射频脉冲,将磁化矢量M0扳转至垂直于静磁场方向的横向平面上。由于分子的扩散及静磁场的空间非均匀性等原因,磁化矢量M0发生散相。经历一定时间τ之后,施加180°脉冲。散相后的磁化矢量会在同等时间τ之后实现重聚,形成一个回波信号。该回波信号称之为自旋回波信号。结合本专利实现的带有梯度磁场分布的磁体系统,自旋回波技术可实现以下两方面的重要功能:(1)在梯度磁场下通过改变梯度幅值或梯度持续时间,记录自旋回波幅值的变化,可以得到流体分子的自扩散系数;(2)通过施加成对的频率编码或相位编码梯度,解析被测样品的空间自旋密度信息,实现核磁共振成像。
自扩散系数Ds反映分子的扩散快慢程度。由于分子的扩散过程为随机运动,在一定时间之后的扩散传播函数或扩散概率密度符合高斯分布。当分子在梯度磁场中扩散,其一定时间内信号的改变与分子的平均扩散位移有关,可通过这一规律计算得到分子的自扩散系数。由于存在静磁场非均匀性及分子的扩散时,对横向磁化矢量造成的额外衰减需要进一步考虑。由磁场非均匀性造成的磁化矢量衰减可以通过重聚脉冲(180°脉冲)予以消除,但分子的扩散造成的磁化矢量衰减无法消除。在宏观角度上来说,该响应是一个与时间相关的函数。因此在考虑分子扩散的影响后,磁化矢量衰减会同时受到横向弛豫时间T2,以及分子自扩散运动在梯度磁场中的响应量。通常采用脉冲磁场梯度或静磁场梯度实现扩散系数的测量。以自旋回波脉冲序列为例,静磁场梯度情况下回波幅值衰减随实验参数的变化公式分别为(由于选取的半回波间隔
δ小于T2,由T2造成的影响可忽略不计):其中,
γ为质子的旋磁比,
γ/2π=42.58MHz/T,G为静磁场梯度值,δ为半回波间隔,Ds为分子的自扩散系数。
二维核磁共振图像可以通过一种无损的方式快速反映被测样品内部流体的空间分布情况。在本发明中,由于梯度的存在,首先可以获得被测样品某一方向上的一维投影剖面。一维投影剖面p函数是通过对得到的自旋回波信号进行傅立叶变换得到的:其中,G为磁场梯度,n为采集自旋回波信号中的点数,td为采集回波中点与点的时间间隔,为每个时间点td下的磁化矢量大小,即核磁共振信号幅度值,Z代表成像位置轴,因此,通过对得到的单个回波整体信号进行傅立叶变换,即可得到一维核磁共振的成像剖面p函数,通过本专利中阐述的调节装置,改变梯度静磁场与均匀静磁场之间的夹角φ,可得到一系列的投影剖面p函数。通过反拉东变换,既可以重构得到被测样品的二维平面图像:
I(x,y)=∫p(xcosφ+ysinφ,φ)dφ
其中,I(x,y)为二维图像结果,φ为静磁场梯度G与静磁场磁感应强度B0方向的夹角,p函数为在φ夹角下的一维磁共振成像剖面。
需要强调的是以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (2)
1.一种能够形成变梯度静磁场的磁体系统结构的测量方法,其特征在于,所述测量方法基于一种能够形成变梯度静磁场的磁体系统结构;
所述能够形成变梯度静磁场的磁体系统结构适用于核磁共振检测装置系统,无需额外搭建梯度成像线圈;
所述能够形成变梯度静磁场的磁体系统结构由磁体部分,固定磁体的外壳部分和磁场调节部分组成,其中磁体部分由两层环形的磁体阵列组成,分别为内环磁体阵列和外环磁体阵列;所述内环磁体阵列和外环磁体阵列均由沿着柱体轴向方向的n层磁环排列组成,所述n层磁环由2N个磁棒构成,n和N为大于2的自然数,阵列组成存在2N-1对耦合级磁棒,每一对耦合级磁棒按照沿柱体轴对称方式排布,每一层磁环的阵列组成均可相同或不同:所述内环磁体阵列的耦合级磁棒的磁化方向相同,所述外环磁体阵列的耦合级磁棒的磁化方向相反;
每一层磁环的阵列组成沿着圆周方向排列,每层磁环中磁棒的厚度和层与层之间的距离,以及每一层磁环的阵列组成的圆周的直径均由迭代算法优化设计得出,所述内环磁体阵列的磁棒的圆周排列方式形成的腔内产生第一静磁场,所述第一静磁场具有特定磁场强度,所述第一静磁场方向与邻近的磁棒磁化方向相同;
所述外环磁体阵列的磁棒的圆周排列方式形成的腔内产生第二静磁场,所述第二静磁场具有恒定梯度,所述第一静磁场和所述第二静磁场形成变梯度静磁场的所述磁体系统结构,所述变梯度静磁场通过改变内环磁体阵列与所述外环磁体阵列的相对角度而改变静磁场梯度的空间取向;所述内环磁体阵列与所述外环磁体阵列由各自的内侧金属骨架和外侧金属骨架安装固定;
固定磁体的外壳部分由所述金属骨架和温控部分组成,所述温控部分安装于所述内侧金属骨架和外侧金属骨架之间,用于保证所述内环磁体阵列与所述外环磁体阵列的相对温度不变,由高精度温度传感器和温度调节系统组成,能够有效防止磁体的温度漂移导致的磁场变化;
调节静磁场梯度的空间取向通过旋转每环磁体阵列的外壳实现自由转动,所述外壳下部有固定孔,用于连接磁场调节部分;
所述磁场调节部分由底盘和无磁电机组成,底盘为两层环状结构,分别对应内环磁体阵列与所述外环磁体阵列的外壳部分,无磁电机带动底盘转动,所述磁场调节部分由控制器操作,控制器根据指令动作控制无磁电机旋转,从而能够使所述磁体部分之间相互转动;
所述测量方法包括:取N=4,n=4,初始化系统和恒温温度,调整无磁电机至初始位置,此时标记静磁场梯度G的空间取向φ=0,将被测样品放置于所述磁体系统结构的测试腔内,在梯度方向下进行一维投影,首先对被测样品施加90°射频脉冲,将磁化矢量M0扳转至垂直于静磁场方向的横向平面上,M0由静磁场强度B0、温度参数决定;由于分子的扩散及静磁场的空间非均匀性原因,磁化矢量M0发生散相;经历时间τ之后,施加180°脉冲;散相后的磁化矢量会在同等时间τ之后实现重聚,形成一个回波信号,然后进行以下步骤:(1)在梯度磁场下通过改变梯度幅值或梯度持续时间,记录回波幅值的变化,得到流体分子的自扩散系数;(2)通过施加成对的频率编码或相位编码梯度,解析被测样品的空间自旋密度信息,实现核磁共振成像;
实现对在静磁场中形成的所述磁化矢量M0的扳转操作由射频磁场B1完成,扳转角度为:θ=γB1tp,其中B1为射频磁场强度,tp为射频脉冲的持续时间,γ为质子的旋磁比;通过控制射频脉冲的幅值或持续时间达到改变扳转角的目的;
步骤(1)中由于存在静磁场非均匀性及分子的扩散时,对横向磁化矢量造成的额外衰减,考虑分子的扩散造成的磁化矢量衰减,磁化矢量衰减是一个与时间相关的函数,因此在考虑分子扩散的影响后,磁化矢量衰减满足以下公式:其中为静磁场B0的梯度,γ/2π=42.58MHz/T,δ为半回波间隔。
2.根据权利要求1所述的一种能够形成变梯度静磁场的磁体系统结构的测量方法,其特征在于,步骤(2)中通过一种无损的方式快速反映被测样品内部流体的空间分布情况,由于梯度的存在,首先获得被测样品某一方向上的一维投影剖面,一维投影剖面p函数是通过对得到的自旋回波信号进行傅立叶变换得到的其中,为静磁场B0的梯度,n为采集自旋回波信号中的点数,td为采集回波中点与点的时间间隔,为每个时间点td下的磁化矢量大小,即核磁共振信号幅度值,Z代表成像位置轴,因此,通过对得到的单个回波整体信号进行傅立叶变换,即可得到一维核磁共振的成像剖面p函数,通过所述磁场调节部分,改变梯度静磁场与均匀静磁场之间的夹角φ,得到一系列的投影剖面p函数,通过反拉东变换,既可以重构得到被测样品的二维平面图像
I(x,y)=∫p(xcosφ+ysinφ,φ)dφ
其中,I(x,y)为二维图像结果,φ为静磁场梯度G与静磁场磁感应强度B0方向的夹角,p函数为在φ夹角下的一维磁共振成像剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254930.3A CN107064838B (zh) | 2017-04-25 | 2017-04-25 | 一种能够形成变梯度静磁场的磁体系统结构及测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254930.3A CN107064838B (zh) | 2017-04-25 | 2017-04-25 | 一种能够形成变梯度静磁场的磁体系统结构及测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064838A CN107064838A (zh) | 2017-08-18 |
CN107064838B true CN107064838B (zh) | 2023-04-28 |
Family
ID=59599915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710254930.3A Active CN107064838B (zh) | 2017-04-25 | 2017-04-25 | 一种能够形成变梯度静磁场的磁体系统结构及测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064838B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728088B (zh) * | 2017-09-23 | 2019-09-17 | 北京青檬艾柯科技有限公司 | 一种快速核磁共振t1成像方法 |
CN107703174B (zh) * | 2017-10-30 | 2020-05-08 | 中国石油天然气股份有限公司 | 一种核磁共振流体分析仪及其制备方法 |
CN108287173A (zh) * | 2018-01-12 | 2018-07-17 | 中国石油大学(北京) | 便携式核磁共振探测装置 |
CN108761361B (zh) * | 2018-03-28 | 2020-01-31 | 中国石油大学(北京) | 可调式高均匀度核磁共振传感器磁体结构及测量装置 |
CN112330674B (zh) * | 2020-05-07 | 2023-06-30 | 南京信息工程大学 | 一种基于脑部mri三维图像置信度的自适应变尺度卷积核方法 |
CN114983387A (zh) * | 2022-05-31 | 2022-09-02 | 重庆邮电大学 | 一种低成本移动式超低场核磁共振成像装置 |
CN116520206B (zh) * | 2023-06-15 | 2023-09-22 | 潍坊新力超导磁电科技有限公司 | 一种梯度磁场自动测量装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2320597A1 (en) * | 2000-01-06 | 2001-07-06 | Blacklight Power, Inc. | Ion cyclotron power converter and radio and microwave generator |
CN1407351A (zh) * | 2001-09-10 | 2003-04-02 | 施卢默格海外有限公司 | 应用核磁共振在井眼中测量流速的方法和装置及其应用 |
CN1423743A (zh) * | 1999-11-16 | 2003-06-11 | 沃林企业股份有限公司 | 磁性共振分析流量计及流量测量方法 |
CN102330548A (zh) * | 2011-08-08 | 2012-01-25 | 中国石油大学(北京) | 消除振铃的核磁共振回波串获取方法 |
CN102481114A (zh) * | 2009-04-03 | 2012-05-30 | 特西奥普技术有限公司 | 用于磁共振成像的信号获取与处理的方法和设备 |
CN102519999A (zh) * | 2011-11-11 | 2012-06-27 | 中国石油大学(北京) | 核磁共振分析仪和核磁共振测量方法 |
CN103744043A (zh) * | 2014-01-20 | 2014-04-23 | 厦门大学 | 一种在不均匀磁场下获取一维高分辨核磁共振谱图的方法 |
CN105137374A (zh) * | 2014-06-03 | 2015-12-09 | 中国科学院上海微系统与信息技术研究所 | 一种超高分辨率的磁共振成像方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102152774B1 (ko) * | 2013-03-07 | 2020-09-07 | 삼성전자 주식회사 | 반도체 소자 시뮬레이션 시스템 및 이를 이용한 시뮬레이션 방법 |
-
2017
- 2017-04-25 CN CN201710254930.3A patent/CN107064838B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1423743A (zh) * | 1999-11-16 | 2003-06-11 | 沃林企业股份有限公司 | 磁性共振分析流量计及流量测量方法 |
CA2320597A1 (en) * | 2000-01-06 | 2001-07-06 | Blacklight Power, Inc. | Ion cyclotron power converter and radio and microwave generator |
CN1407351A (zh) * | 2001-09-10 | 2003-04-02 | 施卢默格海外有限公司 | 应用核磁共振在井眼中测量流速的方法和装置及其应用 |
CN101126816A (zh) * | 2001-09-10 | 2008-02-20 | 施卢默格海外有限公司 | 应用核磁共振在井眼中测量流速的方法和装置及其应用 |
CN102481114A (zh) * | 2009-04-03 | 2012-05-30 | 特西奥普技术有限公司 | 用于磁共振成像的信号获取与处理的方法和设备 |
CN102330548A (zh) * | 2011-08-08 | 2012-01-25 | 中国石油大学(北京) | 消除振铃的核磁共振回波串获取方法 |
CN102519999A (zh) * | 2011-11-11 | 2012-06-27 | 中国石油大学(北京) | 核磁共振分析仪和核磁共振测量方法 |
CN103744043A (zh) * | 2014-01-20 | 2014-04-23 | 厦门大学 | 一种在不均匀磁场下获取一维高分辨核磁共振谱图的方法 |
CN105137374A (zh) * | 2014-06-03 | 2015-12-09 | 中国科学院上海微系统与信息技术研究所 | 一种超高分辨率的磁共振成像方法及装置 |
Non-Patent Citations (1)
Title |
---|
李新.优化重聚脉冲提高梯度场核磁共振信号强度.《物理学报》.2013,第62卷(第14期),147602-1至147602-7. * |
Also Published As
Publication number | Publication date |
---|---|
CN107064838A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064838B (zh) | 一种能够形成变梯度静磁场的磁体系统结构及测量方法 | |
CA2683411C (en) | Magnetic resonance imaging apparatus and method | |
JP3402916B2 (ja) | 核スピントモグラフ装置のマグネット装置のシム調整方法及び方法を実施する装置 | |
US6674282B2 (en) | Method and apparatus for high resolution ex-situ NMR spectroscopy | |
US10261148B2 (en) | Magnetic resonance imaging apparatus and manufacturing method thereof | |
JP3682627B2 (ja) | 磁気共鳴撮像装置 | |
US8890527B1 (en) | Methods of radio frequency magnetic field mapping | |
Perlo et al. | Single-sided sensor for high-resolution NMR spectroscopy | |
US5229722A (en) | Nqr-imaging | |
Vogel et al. | 3D-Spatial encoding with permanent magnets for ultra-low field magnetic resonance imaging | |
Yang et al. | A low-cost, miniature Halbach magnet designed for portable time domain NMR | |
US20140218025A1 (en) | Transverse volume coils and related magnetic resonance systems and methods | |
Sørland et al. | Artefacts and pitfalls in diffusion measurements by NMR | |
Paulsen et al. | Volume-selective magnetic resonance imaging using an adjustable, single-sided, portable sensor | |
Han et al. | Two-dimensional representation of position, velocity and acceleration by PFG-NMR | |
JP2000262490A (ja) | 磁気共鳴イメージング方法および磁気共鳴イメージング装置 | |
JP5212972B2 (ja) | 計測装置および計測方法 | |
US5317262A (en) | Single shot magnetic resonance method to measure diffusion, flow and/or motion | |
USH1218H (en) | NMR imaging with varying spatial coupling | |
JPH05500169A (ja) | 核磁気共鳴検査に於ける領域選択 | |
Bartusek et al. | Mapping of magnetic field around small coils using the magnetic resonance method | |
CN206906567U (zh) | 一种能够形成变梯度静磁场的磁体系统结构 | |
Rodriguez et al. | Using Proton Nuclear Magnetic Resonance (NMR) as a calibrating reference for magnetic field measurement instruments: Sensitive volume and magnetic field homogeneity | |
US20230280427A1 (en) | Spatial encoding arrangement | |
Lebedev et al. | Localized Gradients |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |