CN112834543A - 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 - Google Patents
基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 Download PDFInfo
- Publication number
- CN112834543A CN112834543A CN202010347963.4A CN202010347963A CN112834543A CN 112834543 A CN112834543 A CN 112834543A CN 202010347963 A CN202010347963 A CN 202010347963A CN 112834543 A CN112834543 A CN 112834543A
- Authority
- CN
- China
- Prior art keywords
- gradient
- layer selection
- spectrum
- frequency
- pulse
- 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.)
- Granted
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 70
- 230000005291 magnetic effect Effects 0.000 claims abstract description 59
- 238000001228 spectrum Methods 0.000 claims abstract description 45
- 238000005070 sampling Methods 0.000 claims abstract description 33
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 29
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 21
- 230000008901 benefit Effects 0.000 claims abstract description 19
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 claims abstract description 9
- 238000001816 cooling Methods 0.000 claims abstract description 8
- 230000005284 excitation Effects 0.000 claims description 16
- 230000004044 response Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 5
- 238000010998 test method Methods 0.000 claims description 4
- 238000005057 refrigeration Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 2
- 239000011148 porous material Substances 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 11
- 239000011435 rock Substances 0.000 description 11
- 239000012530 fluid Substances 0.000 description 9
- 239000000126 substance Substances 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 7
- 230000005294 ferromagnetic effect Effects 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 230000018109 developmental process Effects 0.000 description 6
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 230000005298 paramagnetic effect Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000013517 stratification Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 229910052761 rare earth metal Inorganic materials 0.000 description 2
- 150000002910 rare earth metals Chemical class 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- QJVKUMXDEUEQLH-UHFFFAOYSA-N [B].[Fe].[Nd] Chemical compound [B].[Fe].[Nd] QJVKUMXDEUEQLH-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000010997 low field NMR spectroscopy Methods 0.000 description 1
- 229910001172 neodymium magnet Inorganic materials 0.000 description 1
- 231100000956 nontoxicity Toxicity 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了基于脉冲梯度硬件结构的一维空间选层T2谱测试方法，在梯度线圈硬件上增加循环水冷模块，利用成像类的脉冲梯度硬件结构，实现长延时恒定梯度场输出，满足选层测试期间CPMG序列的采样需要，兼容核磁共振分析仪器的磁体平台。针对恒定梯度场选层测试的工作原理，设计适合脉冲梯度硬件结构的恒定梯度模式一维空间选层T2谱测试专用序列。本发明在现有核磁共振分析仪器平台基础上，充分融合了脉冲梯度编码方法、恒定梯度场编码方法各自的优势，针对性的增加和提高部分硬件模块功能实现了一维空间选层T2谱测试功能，既保持了恒定梯度场方法TE不受限制的优点，又保持了脉冲梯度场方法磁体平台结构简单等优点。
Description
技术领域
本发明涉及多孔介质材料实验测量领域，具体涉及一种基于脉冲梯度硬件结构的一维空间选层T2谱测试方法。
背景技术
近几年随着国内大部分油气田开发进入中后期，提高采收率(EOR)成为油气开发的重要工作。油气开发实验中的剩余油空间分布测量一直是EOR实验的重点和难点，填砂管、玻璃刻蚀等传统实验技术难以真实模拟储层渗流状态，给精细化研究提高采收率方案造成困难。进入二十一世纪以来，随着CT、NMR等成像技术在油气勘探开发领域中的应用和普及，为剩余油空间分布测量创造了有利条件。NMR技术因其快速、无损、无毒、只对含氢孔隙流体敏感等特点，已在石油勘探领域已经得到成熟应用，钻井现场的裸眼井测井、随钻测井、地层测试器、岩屑录井等有相应的NMR仪器；近几年，NMR技术在石油开发领域中的应用也越来越广泛。
和勘探领域测量样品整体NMR信号不同的是，EOR实验更需要NMR技术提供岩心的流体空间分布，核磁共振成像技术原理上能够满足该需求。对物体进行核磁共振成像，需要三维梯度定位系统，最常用的是三维正交梯度线圈，实现样品区域内的三维空间梯度编码。梯度场施加后，相对于坐标原点，空间位置
处的采样信号为：
式中
是采样信号；
是成像序列及梯度硬件参数确定的k空间；
是样品中质子密度空间分布，也即岩心中孔隙流体的空间分布。对核磁成像信号
用已知的k空间做傅里叶逆变换即可得到岩心孔隙流体的三维空间分布。
附图1是3种常用核磁共振成像序列示意图：MSE序列在自旋回波(Spin Echo)序列基础上增加了三维梯度空间编码，实现了三维空间的定位成像。FSE序列用了类似CPMG序列的原理，一次激发采集多个后续回波信号，对每次的回波信号都做三维梯度空间编码，把多个回波信号累加提高信噪比，进而降低扫描次数实现快速扫描。GRE序列采用小角度激发，降低重复采样等待时间，缩短累加采样的总时长。因此，FSE、GRE序列一般比MSE序列成像速度快。
岩心等多孔介质，特别是致密类岩心，核磁属性典型特征是弛豫时间短、信号量低，给常用核磁共振成像序列的应用带来严重挑战。核磁共振成像技术起源于医学，测量对象是人体、小动物等中长弛豫时间、高含水体系，因此公式1中没有考虑梯度编码期间的弛豫衰减。针对中短弛豫时间的岩心样品，考虑梯度编码期间弛豫衰减的成像信号表达式如下：
附图1所示的常用成像序列的射频激发脉冲都是软脉冲，软脉冲的特点是频率域带宽小、选层特性好，缺点就是脉宽大，在ms量级。软脉冲成像序列，射频脉冲宽度，加上梯度编码时间(也是ms量级)，导致产生回波信号时的TE至少在5～10ms，导致岩心中小孔隙信号丢失严重。为减少短弛豫信号损失，必须缩短成像序列的回波间隔TE，因此出现了硬脉冲三维成像序列，如附图2所示。HMSE序列利用三维梯度空间编码实现成像的同时，发挥硬脉冲持续时间短(us量级)的巨大优势，将成像的TE缩短至2～3ms，保证了短弛豫信号的采集。
岩心成像除了对TE有更高要求外，SNR受限的影响更大。岩心成像参数：选层数量10层，每层像素点是100*100，假设孔隙流体在岩心内部均匀分布，则图像上每个像素点的信号是整体信号的十万分之一。岩心核磁共振信号本身就是微弱信号，再做三维成像分割，将导致成像像素点的SNR急剧下降。
综上所述，成像序列梯度编码时间占用回波时间，加上岩心样品中短弛豫时间、弱信号典型特征，导致常用的二维或三维核磁共振成像方法很难在致密岩心测试中达到理想效果。
为了弥补三维核磁共振成像方法在SNR方面的缺点，结合EOR实验更关注驱替方向上的剩余油变化，逐渐发展出只做一维选层，层内不再做二维扫描的一维空间选层测试方法：相对于三维核磁共振成像序列十万分之一的信号离散，几层～数十层的一维选层对信号的离散仍能够保证采样数据的SNR。按照一维空间选层方法的不同，可以分为三大类：频率编码方法，相位编码方法，恒定梯度编码方法。
频率编码一维空间选层T2测试序列如附图3所示：在90°和第一个180°脉冲之间进行频率编码，在后续回波信号采集时施加读频率编码梯度，因此每个回波包含所有频率(对应样品不同空间位置)的信号。因此，频率编码一维空间选层T2测试的响应信号表达式如下：
式中S(N·TE,t)是采样信号；TE是回波间隔，N是回波序号；z是选层方向，G是选层方向上的梯度值，t是梯度编码时间。
对每个回波信号S(N·TE,t)先做傅里叶逆变换，得到不同位置处的回波信号，所有回波信号组合得到不同位置处的回波衰减信号M(N·TE,z)，M(N·TE,z)用BRD等T2反演算法解谱即可得到一维空间选层T2谱f(T2,z)，如公式4所示。
相位编码一维空间选层T2测试序列如附图4所示：在90°和第一个180°脉冲之间施加相位编码梯度，后续回波采集使用最短回波间隔TE进行连续采样，改变相位编码梯度，实现对样品一维空间上的完整信号采集。因此，相位编码一维空间选层T2测试的响应信号表达式如下：
式中式中S(N·TE,m)是采样信号；TE是回波间隔，N是回波序号；z是选层方向，Gmax是选层方向上梯度最大值，tp是相位编码时间，m是相位编码步数对应的等差数列。
和频率编码方法的数据处理方法类似，也是先对采样信号S(N·TE,m)做傅里叶逆变换，得到不同位置处的回波衰减信号M(N·TE,z)，M(N·TE,z)用BRD等T2反演算法解谱即可得到一维空间选层T2谱f(T2,z)，如公式6所示。
恒定梯度编码一维空间选层T2测试序列如附图5所示：采样模式和CPMG序列一样，连续SE回波采集回波衰减信号，采样期间存在一个恒定梯度场实现梯度编码选层功能。
和频率编码、相位编码方法最大的区别：梯度编码在整个采样期间完成，因此回波采集一直使用最短TE采样；只有选中的频率位置会被激发采样，其他位置的样品处于待激发采样状态，因此变化位置采样不需要等待时间。因此，恒定梯度编码一维空间选层T2测试的响应信号表达式如下：
式中式中S(N·TE,z)是采样信号；TE是回波间隔，N是回波序号；z是选层方向。恒定梯度编码方式是选择性激发采样，直接对不用位置处的回波信号S(N·TE,z)反演即可得到一维空间选层T2谱f(T2,z)。
从上述三种方法的基本工作原理可知，频率编码、相位编码本质上仍是核磁共振医学成像类的技术，虽然一维梯度编码相比于三维梯度编码的时间短，但是梯度编码依旧占用回波时间，对短弛豫信号采集不利。恒定梯度编码方法类似核磁共振测井仪的多频测量技术，利用主磁场是梯度型磁场的特点，实现选层(变化激发频率)激发采样，因此梯度编码不占用回波时间，最小TE采样充分保证短弛豫信号的采集。成像类的一维选层T2谱测试方法也继承了核磁成像选层参数灵活、硬件结构简单等优点，而恒定梯度编码方法在这方面却不占优势。因此三种方法各有优缺点，其详细对比如表1所示。
表1三种一维空间选层T2谱测试方法的优缺点及性能对比
综上所述，成熟的一维选层T2谱测试方法主要分为两大类：基于脉冲梯度场的成像类方法(频率编码、相位编码)，基于恒定梯度场的方法。这两种方法各有优缺点，目前还没有一种兼容两种方法优点的新技术、新方法。
发明内容
本发明目的是：提供一种基于脉冲梯度硬件结构的一维空间选层T2谱测试方法，兼容上述两大类方法，既保持恒定梯度场方法TE不受限制的优点，又保持了脉冲梯度场方法磁体平台结构简单等优点。
本发明的技术方案是：
基于脉冲梯度硬件结构的一维空间选层T2谱测试方法，包括以下步骤：
(1)使用分析仪器常用的均匀磁场作为主磁场B0，选用成像梯度线圈硬件结构施加梯度场，在CPMG序列回波采样期间施加幅度不变的梯度实现选层功能；
(2)测量不同位置信号通过切换射频激发频率，和该位置恒定梯度场下的共振频率相同，实现层位选择；
(3)在梯度线圈硬件上增加循环水冷模块，通过模块中的制冷设备将梯度线圈产生的热量带走，使其在CPMG回波采样期间保持恒定的梯度输出，且持续时间至少保持在1～2s，达到保护梯度硬件不被烧坏的目的；
(4)选用高功率功放电源，同时通过累加采样时的等待时间TW来控制梯度功放电源输出占空比，保证给梯度线圈供电的功放电源能够持续输出1～2s，达到保护功放电源不过载的目的。
优选的，步骤(3)中在梯度线圈硬件上增加循环水冷模块，利用成像类的脉冲梯度硬件结构，增加循环水冷模块，实现长延时恒定梯度场输出，满足选层测试期间CPMG序列的采样需要，兼容核磁共振分析仪器的磁体平台。
优选的，针对恒定梯度场选层测试的工作原理，设计适合脉冲梯度硬件结构的恒定梯度模式一维空间选层T2谱测试专用序列，完整继承梯度型主磁场恒定梯度一维空间选层T2谱测试方法T2测试下限值小的优点，为满足中短T2弛豫特征的多孔介质一维空间选层T2谱测试创造条件。
优选的，还设计了基于均匀水膜样品测试提前获取核磁共振信号校正常数的方法，解决实际样品测试时，恒定梯度场模式下射频线圈发射、接收中心频率和射频线圈S11曲线谐振点不严格相等造成的核磁信号不一致问题，保证一维空间选层T2谱测试的数据一致性。
优选的，所述的一维空间选层T2谱测试方法，从数据采集角度看，恒定梯度编码方式的响应信号表达式如下：
式中式中S(N·TE,z)是采样信号；TE是回波间隔，N是回波序号；z是选层方向；恒定梯度编码方式是选择性激发采样，直接对不用位置处的回波信号S(N·TE,z)反演即可得到一维空间选层T2谱f(T2,z)。
优选的，所述的一维空间选层T2谱测试方法，从硬件实现角度看，和传统的成像类编码方法类似，梯度场只是在需要的时候施加；并不像恒定梯度场方法一样，主磁场B0就是一个恒定梯度场。因此本发明提出的一维空间选层T2谱测试方法集成了成像类方法的技术优势，主磁场B0仍然是均匀场，保证了磁体平台的通用性，避免了主磁场B0是梯度型磁场导致的测试样品小、磁体结构复杂等缺陷。和传统成像类梯度硬件结构不同的是，需要让梯度场持续输出1～2s，因此硬件结构上需要增加循环水冷模块等。
本发明的优点是：
本发明在现有核磁共振分析仪器平台基础上，充分融合了脉冲梯度编码方法、恒定梯度场编码方法各自的优势，针对性的增加和提高部分硬件模块功能实现了一维空间选层T2谱测试功能，既保持了恒定梯度场方法TE不受限制的优点，又保持了脉冲梯度场方法磁体平台结构简单等优点，因此和传统的三种方法的详细对比如表2所示。
表2本发明提供方法和传统方法的技术对比
附图说明
下面结合附图及实施例对本发明作进一步描述：
图1三种常用核磁共振成像序列；
图2硬脉冲多层自旋回波(HMSE)成像序列；
图3频率编码一维选层T2测试序列；
图4相位编码一维选层T2测试序列；
图5恒定梯度编码一维选层T2测试序列；
图6基于脉冲梯度硬件结构的恒定梯度一维空间选层T2谱测试技术；
图7不同样品的fid及其对应的频谱曲线；
图8多孔介质中局部内部梯度方向示意图；
图9硬脉冲及其对应的频谱曲线；
图10射频线圈S11曲线；
图11射频激发频率、带宽与样品拉莫尔频率、带宽对应关系对测量的影响；
图12恒定梯度场下选层测试信号校正方法示意图；
图13三种一维空间选层T2测试方法正反演数值模拟研究流程图；
图14f(T2,z)构造模型及内部梯度场模型；
图15三种方法对应的f(T2,z)反演谱；
图16三种方法对应的T2谱相对误差及T2下限值；
图17不同弛豫时间均匀水膜三种方法的实测结果。
具体实施方式
本发明具体的实施方式如下。
实现核磁共振的选层功能需要了解“共振”的定义和对测量的影响等基础概念。核磁共振中的“共振”指的是在频率域内，当射频线圈的激发频率等于样品中的质子拉莫尔频率时，样品中的质子会发生共振并吸收射频能量，在射频场撤销后，吸收了射频能量后的高能态质子会弛豫衰减释放能量回到低能态，质子这种弛豫衰减释能过程的快慢受到质子所处环境的影响，因此把质子弛豫衰减释能过程通过射频线圈接收并记录下来，得到弛豫衰减曲线，弛豫衰减曲线经过数学反演得到反映质子所处环境特征的T2谱。
因此，射频激发频率和拉莫尔频率的匹配性会直接影响射频能量的吸收效率，进而影响核磁数据SNR。质子的拉莫尔频率：
f＝γB0/2π (9)
B0是仪器的磁场强度，虽然主磁场的表述经常是均匀场，但是绝对均匀的磁场(即△B0＝0)是做不到的，只是磁场的不均匀性可以做到很小，因此磁体指标里都有磁场均匀性指标(在磁体中心一定空间范围内，最大磁场强度和最小磁场强度的差与平均磁场强度的比值，这个空间范围就是磁场均匀区)。比如低场核磁共振磁体的磁场均匀性一般在30～100ppm范围内，ppm指百万分之一。由于主磁场的不均匀，样品里的质子拉莫尔频率是有带宽的；同时由于多孔介质样品内部梯度G0的影响，也会叠加到样品的拉莫尔频率上去，进一步增大样品拉莫尔频率带宽。因此，主磁场中，多孔介质的拉莫尔频率如下式：
f＝γ(B0+ΔB0+ΔBG0)/2π (10)
多孔介质的内部梯度G0和固体骨架中的铁磁性物质的含量有关，因此顺磁性物质含量越高的样品，G0越大，样品的拉莫尔频率带宽也越大。核磁共振的自由衰减曲线fid可以反映样品的频率带宽，如附图7所示：从附图7可知，不同样品中的质子虽然都在主磁场B0下，但是共振频率带宽是有差异的。纯水里没有固体骨架，所以没有内部梯度G0，因此纯水的频谱曲线能够反映主磁场B0的均匀性。频率带宽一般用半高宽(震幅峰值的一半对应的频率差)表示，如图中水平红色虚线双箭头所示。附图7的仪器主磁场强度B0＝12MHz，纯水的频率带宽是0.6kHz左右，因此该磁体的均匀度＝1000*0.6/12＝50ppm。多孔介质的频率带宽都是大于主磁场的频率带宽的，而且多孔介质的频率带宽与其铁磁性物质含量(准确的讲应该是样品的综合磁化率)成正比。
从附图7可知，铁磁性物质很少的纯净岩心，频率带宽只比主磁场的频率带宽稍大。虽然铁磁性物质的含量越高，对应的频率带宽越大，但是多孔介质的频率带宽是有上限的。这是因为常见的多孔介质材料(特别是天然的多孔介质材料)，非固态质子所处的孔隙在多孔介质中的空间分布是杂乱的，顺磁性物质在骨架中的分布也是不均匀的，所以多孔介质的内部磁场梯度G0(固体骨架与孔隙流体的磁化率不同)很难在孔隙空间内定向累加到很大数值，除非强铁磁性矿物在局部孔隙周围富集，如附图8所示。只有人工合成的稀土材料永磁体能够在间隙区形成强磁场，虽然稀土材料也是从地球矿产中的岩石中提取出来的，但是自然状态下铁磁性矿物的纯度和结构所能产生的磁场强度和人造永磁体的磁场强度相比，要小好几个数量级。
除了多孔介质样品的共振频率带宽，还有射频的激发、接收带宽也会影响核磁共振测量。决定射频激发带宽的两个主要因素是：1、射频脉冲的类型及时间参数；2、射频线圈收发能量的频率特征(S11回波损耗曲线)。
多孔介质核磁共振测试最常用的射频脉冲是硬脉冲，时间域是一个矩形方波，在频率域也是有频率带宽的，如附图9所示。频率域和时间域是倒数关系，所以脉宽越小的脉冲，频率带宽越大。一般低场核磁共振仪器的硬脉冲脉宽在几us～几十us，因此其带宽在几十kHz～几百kHz之间。射频脉冲对应的频率就是频率发生器的发射主频+脉宽对应的频率带宽(附图9中没有标注发射主频，只显示了脉宽对应的频带)。与样品共振频率的带宽类似，硬脉冲对应的频率带宽常用脉宽倒数近似计算。
从附图7、图9可知，硬脉冲的激发频率带宽一般比多孔介质样品共振频率带宽大1到2个数量级，这个条件只是保证射频源提供的能量的频率范围能够充分覆盖多孔介质样品共振吸收能量的频率范围。射频能量能否高效的传递给样品，样品弛豫衰减过程中释放的能量如何采集并记录，这两个过程涉及到核磁共振的关键部件，射频线圈。
螺线管线圈是低场核磁共振仪器最常采用的射频线圈，其最大特征就是自发自收，在射频脉冲激发期间负责发射传递射频能量，回波采样期间则负责感应接收样品的弛豫衰减信号。自发自收射频线圈能量发射和接收性能的评价参数是S11回波损耗曲线，如附图10所示。
S11的定义如下：
因此，S11代表了射频发射射频能量的效率，S11越小，射频脉冲能量进入样品的能量越多，对应的核磁信号量也就越大。因为是自发自收线圈，核磁信号接收过程中，S11代表了射频线圈感应样品能量的效率，S11越小，线圈接收到样品弛豫衰减的能量越多，核磁信号越大。
S11极小值点对应的频率通常称为线圈的谐振频率或谐振点。从S11的定义可知，射频发射和接收能量的效率在S11极小值点对应的频率上最高，其他频率处的效率都比S11极小值点处低。一般的发射和接收电路中，通常定义-3dB为系统的有效工作频率范围，其定义是反射功率是输入功率一半时对应的频率带宽，如下式。因此附图9所示的射频线圈的-3dB带宽在35kHz左右。
从能量高效激发和接收的角度看，频率域中，核磁共振理想的测试条件是：1、样品拉莫尔频率的中心频率fL(主磁场B0)、射频源发射频率的中心频率fO(一般就是采集软件界面上的中心频率)、射频线圈谐振频率fR(S11曲线极小值点对应的频率)三个频率相等，即fL＝fO＝fR。2、样品拉莫尔频率的带宽△fL，射频脉冲的带宽△fO，射频线圈带宽△fR，三个带宽满足△fL<<min(△fO，△fR)。3、采样带宽SW>△fL，一般核磁共振谱仪的接收机采样带宽是可调的，而且调节范围基本都在MHz数量级，完全能够满足测试要求。
上述三个条件中，最关键的是条件2，因为只要射频的激发带宽能够覆盖样品的拉莫尔频率带宽，fL、fO、fR三者不严格相等仍然能够测到核磁信号，只是会影响数据质量和定量测量结果，如附图11所示。图11c对应的情况，就是低场核磁共振仪器中常见的磁体温度变化导致的频率漂移(常用的钕铁硼永磁体的磁场强度B0和温度有线性关系)，简称“频漂”。频漂后(也即fL变化)，如果fL偏移不大，虽不是最佳测量条件，但是依然满足共振条件能够测量，只是采样数据SNR受到一定影响；但是当fL偏移较大导致△fO不能完整覆盖△fL时，就不满足共振条件了，需要重新校准仪器的中心频率，就是把射频源的发射中心频率fO对准到新的拉莫尔频率fL，这一操作简称调频。
通过调频能够把fO实时对准到fL，保证射频源发射的射频能量能够和样品中的质子满足共振条件；射频源的调频很简单也很快，采样软件上可直接操作，但是调频的同时也就意味着射频线圈的谐振点(fR)不等于fO、fL，fR的偏移虽然不会影响质子的共振条件，但是会影响射频线圈的发射和接收效率，进而影响核磁共振信号的一致性。
为了保证核磁共振信号的一致性，同时避免引入射频线圈调谐电路(结构复杂，响应速度慢，会影响测试速度)带来的问题，本发明设计了一种核磁信号校正方法，该方法的示意图如附图12所示：使用一个均匀水膜(长度方向上每个位置处的含氢物质的质量相同)进行一维空间选层T2测试，以频率中心位置处的信号作为基准值，计算两侧不同位置处的信号校正常数。实际样品测试时，通过均匀水膜得到的信号校正常数对不同位置处的信号校正，得到不受射频线圈-3dB带宽影响的核磁共振信号分布。
为了应对不同物性条件岩心的驱替实验测试，更关注的是不同方法对短弛豫信号的识别能力，特别是油田现场开发面临储层致密化的趋势。为直接对比不同方法对短弛豫信号的测量能力，使用数值模拟技术量化不同一维空间选层T2谱测试方法的T2测试下限，具体流程如附图13所示。
为遍历不同孔隙大小对应的T2数值，构造从小到大的T2单峰模型，层厚为2mm，一共50层，T2峰值最小值0.01ms、最大值800ms，如附图14a、14b所示。
储层岩石骨架中含有顺磁性物质或铁磁性物质时，岩石骨架与孔隙流体的磁化率有明显差异，在外部磁场条件下这种磁化率差异将在孔隙空间中形成附加磁场梯度，称为内部梯度场。内部梯度场会通过扩散弛豫加速孔隙流体的T2弛豫衰减，TE越大，内部梯度场影响越大。频率编码、相位编码因为梯度编码需要时间，导致TE较大，内部梯度场的影响不可忽略。内部磁场梯度G0的近似计算公式：
式中，μ0为真空磁导率，B0为仪器磁场强度，Δχ为孔隙流体与骨架的磁化率之差，R为孔隙空间中内部磁场梯度计算点到固液接触面的距离。利用T2与R的正比例关系，简化G0的计算公式如下：
式中，σ为内部梯度强度因子，在仪器主磁场强度B0固定的前提下，σ与岩心骨架顺磁性物质的含量成正比，图14c为σ＝5T·ms/m时的G0与T2的对应关系。
考虑内部梯度场的影响，分别推导出频率编码、相位编码、本发明提供方法对应的一维空间选层T2测试正演响应数学模型，如下式：
式中，TEf为频率编码序列中频率编码窗口的回波间隔；TE1为相位编码序列中梯度编码窗口的回波间隔，TE2为相位编码序列中回波信号采集窗口的回波间隔；G0是内部梯度场大小；Gs是恒定梯度场的梯度大小；noise是随机噪声。
公式15的响应信号处理过程和前述过程一样，频率编码、相位编码方法的正演响应信号先做傅里叶逆变换再做T2反演，本发明提供方法的正演响应信号直接做T2反演。
利用附图14的构造模型和公式15的正演响应数学模型进行数值模拟，重要参数赋值：TEf＝TE1＝2ms，TE2＝TE＝0.1ms，Gs＝0.1T/m，σ＝5T·ms/m，SNR＝500。频率编码、相位编码、本发明提供方法正反演后得到的f(T2,z)反演谱如附图15所示，f(T2,z)反演谱和图14a的f(T2,z)构造谱的相对误差如图16所示。以相对误差5％作为T2下限值的评价指标，在上述模拟条件下，三种方法对应的T2下限值如图16。该模拟条件下，本发明提供方法的T2下限值0.23ms<<相位编码方法的T2下限值4.948ms<频率编码方法的T2下限值13.48ms，说明本发明提供方法在一维空间选层T2测试中具备极大的优势。
为了验证该方法的实际应用效果，制作了三种T2时间的均匀水膜，T2峰值分别为1ms、10ms、30ms，每个均匀水膜的长度20mm，直径相同，含水体积相等；利用相位编码、频率编码、本发明提供方法这三种方法进行一维空间选层T2谱测试实验，对比三种方法的实测结果发现，和数值模拟仿真的结果一致，本发明提供方法的T2下限值最小。
综上，通过理论及实践测试，均证实了本发明提供的基于脉冲梯度硬件结构的恒定梯度一维空间选层T2谱测试技术，实现了传统成像类方法和恒定梯度类方法的优点融合，有利于该技术的市场推广。
需要说明的是，上述实施例仅仅是本发明一部分实施例，而不是全部实施例。基于本发明的技术方法及实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例都属于本发明的保护范围。
Claims (6)
1.基于脉冲梯度硬件结构的一维空间选层T2谱测试方法，其特征在于，包括以下步骤：
(1)使用分析仪器常用的均匀磁场作为主磁场B0，选用成像梯度线圈硬件结构施加梯度场，在CPMG序列回波采样期间施加幅度不变的梯度实现选层功能；
(2)测量不同位置信号通过切换射频激发频率，和该位置恒定梯度场下的共振频率相同，实现层位选择；
(3)在梯度线圈硬件上增加循环水冷模块，通过模块中的制冷设备将梯度线圈产生的热量带走，使其在CPMG回波采样期间保持恒定的梯度输出，且持续时间至少保持在1～2s，达到保护梯度硬件不被烧坏的目的；
(4)选用高功率功放电源，同时通过累加采样时的等待时间TW来控制梯度功放电源输出占空比，保证给梯度线圈供电的功放电源能够持续输出1～2s，达到保护功放电源不过载的目的。
2.根据权利要求1所述的一维空间选层T2谱测试方法，其特征在于：步骤(3)中在梯度线圈硬件上增加循环水冷模块，利用成像类的脉冲梯度硬件结构，增加循环水冷模块，实现长延时恒定梯度场输出，满足选层测试期间CPMG序列的采样需要，兼容核磁共振分析仪器的磁体平台。
3.根据权利要求2所述的一维空间选层T2谱测试方法，其特征在于：针对恒定梯度场选层测试的工作原理，设计适合脉冲梯度硬件结构的恒定梯度模式一维空间选层T2谱测试专用序列，完整继承梯度型主磁场恒定梯度一维空间选层T2谱测试方法T2测试下限值小的优点，为满足中短T2弛豫特征的多孔介质一维空间选层T2谱测试创造条件。
4.根据权利要求3所述的一维空间选层T2谱测试方法，其特征在于：还设计了基于均匀水膜样品测试提前获取核磁共振信号校正常数的方法，解决实际样品测试时，恒定梯度场模式下射频线圈发射、接收中心频率和射频线圈S 11曲线谐振点不严格相等造成的核磁信号不一致问题，保证一维空间选层T2谱测试的数据一致性。
5.根据权利要求1所述的一维空间选层T2谱测试方法，其特征在于：所述的一维空间选层T2谱测试方法，从数据采集角度看，恒定梯度编码方式的响应信号表达式如下：
式中式中S(N·TE,z)是采样信号；TE是回波间隔，N是回波序号；z是选层方向；恒定梯度编码方式是选择性激发采样，直接对不用位置处的回波信号S(N·TE,z)反演即可得到一维空间选层T2谱f(T2,z)。
6.根据权利要求1所述的一维空间选层T2谱测试方法，其特征在于：所述的一维空间选层T2谱测试方法，从硬件实现角度看，梯度场只是在需要的时候施加，主磁场B0仍然是均匀场，保证了磁体平台的通用性。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010347963.4A CN112834543B (zh) | 2020-04-28 | 2020-04-28 | 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010347963.4A CN112834543B (zh) | 2020-04-28 | 2020-04-28 | 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112834543A true CN112834543A (zh) | 2021-05-25 |
| CN112834543B CN112834543B (zh) | 2024-05-14 |
Family
ID=75923112
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010347963.4A Active CN112834543B (zh) | 2020-04-28 | 2020-04-28 | 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112834543B (zh) |
Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0086972A1 (en) * | 1982-02-03 | 1983-08-31 | General Electric Company | Method of NMR imaging which overcomes T2* effects in an inhomogeneous static magnetic field |
| EP0175184A2 (de) * | 1984-09-18 | 1986-03-26 | Bruker Medizintechnik GmbH | Verfahren zum Messen der magnetischen Kernresonanz |
| US20030030436A1 (en) * | 2001-08-08 | 2003-02-13 | Jurgen Hennig | Method for measuring the magnetic resonance (NMR) by steady state signals (SSFP) |
| WO2003040743A1 (en) * | 2001-11-06 | 2003-05-15 | Shell Internationale Research Maatschappij B.V. | Method and apparatus for subterranean formation flow imaging |
| US20030234648A1 (en) * | 2002-06-19 | 2003-12-25 | Krishnamurthy Ganesan | NMR spectroscopy using a gradient field NMR tool |
| US20040189296A1 (en) * | 2003-03-24 | 2004-09-30 | Chevron U.S.A. Inc. | Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins |
| GB2414807A (en) * | 2004-06-04 | 2005-12-07 | Schlumberger Holdings | NMR method for determining a formation fluid property |
| GB2445247A (en) * | 2006-12-22 | 2008-07-02 | Schlumberger Holdings | NMR Apparatus with Sample Chamber |
| EP1999493A2 (en) * | 2006-03-11 | 2008-12-10 | Xigo Nanotools LLC | Compact and portable low-field pulsed nmr dispersion analyzer |
| WO2010003236A1 (en) * | 2008-07-08 | 2010-01-14 | University Of New Brunswick | Spin echo spi methods for quantitative analysis of fluids in porous media |
| JP2010051804A (ja) * | 2008-08-28 | 2010-03-11 | Siemens Ag | 測定シーケンスを形式的に記述して測定シーケンスのプログラミングを自動化するための磁気共鳴シーケンスモデル |
| US20110234220A1 (en) * | 2010-03-24 | 2011-09-29 | Jonathan Mitchell | System and method for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis |
| US20120153952A1 (en) * | 2010-12-16 | 2012-06-21 | General Electric Company | Method and system for gradient linear correction |
| US20120194188A1 (en) * | 2011-02-01 | 2012-08-02 | Zhe Phillip Sun | System and Method for Controlling Apparent Timing Dependencies for T2-Weighted MRI Imaging |
| GB2505232A (en) * | 2012-08-23 | 2014-02-26 | Schlumberger Holdings | Magnetic resonance examination of porous samples |
| US20140253116A1 (en) * | 2013-03-05 | 2014-09-11 | Schlumberger Technology Corporation | System and method for obtaining nuclear magnetic resonance measurements on reservoir fluids for prediction of fluid properties |
| US20140285196A1 (en) * | 2011-09-14 | 2014-09-25 | Petrochina Company Limited | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field |
| CN104280703A (zh) * | 2014-09-28 | 2015-01-14 | 中国石油大学(华东) | 一种基于脉冲梯度磁场的新型d-t2脉冲施加方法 |
| US9864033B1 (en) * | 2013-08-26 | 2018-01-09 | University Of New Brunswick | Free induction decay based magnetic resonance imaging methods |
| US20180364184A1 (en) * | 2017-06-16 | 2018-12-20 | Arcady Reiderman | Apparatus and method for nuclear magnetic resonance measurements on borehole materials |
| US20190346385A1 (en) * | 2018-05-11 | 2019-11-14 | Arcady Reiderman | Method and apparatus for nuclear magnetic resonance measurements on borehole materials |
-
2020
- 2020-04-28 CN CN202010347963.4A patent/CN112834543B/zh active Active
Patent Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0086972A1 (en) * | 1982-02-03 | 1983-08-31 | General Electric Company | Method of NMR imaging which overcomes T2* effects in an inhomogeneous static magnetic field |
| EP0175184A2 (de) * | 1984-09-18 | 1986-03-26 | Bruker Medizintechnik GmbH | Verfahren zum Messen der magnetischen Kernresonanz |
| US20030030436A1 (en) * | 2001-08-08 | 2003-02-13 | Jurgen Hennig | Method for measuring the magnetic resonance (NMR) by steady state signals (SSFP) |
| WO2003040743A1 (en) * | 2001-11-06 | 2003-05-15 | Shell Internationale Research Maatschappij B.V. | Method and apparatus for subterranean formation flow imaging |
| US20030234648A1 (en) * | 2002-06-19 | 2003-12-25 | Krishnamurthy Ganesan | NMR spectroscopy using a gradient field NMR tool |
| US20040189296A1 (en) * | 2003-03-24 | 2004-09-30 | Chevron U.S.A. Inc. | Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins |
| GB2414807A (en) * | 2004-06-04 | 2005-12-07 | Schlumberger Holdings | NMR method for determining a formation fluid property |
| EP1999493A2 (en) * | 2006-03-11 | 2008-12-10 | Xigo Nanotools LLC | Compact and portable low-field pulsed nmr dispersion analyzer |
| GB2445247A (en) * | 2006-12-22 | 2008-07-02 | Schlumberger Holdings | NMR Apparatus with Sample Chamber |
| WO2010003236A1 (en) * | 2008-07-08 | 2010-01-14 | University Of New Brunswick | Spin echo spi methods for quantitative analysis of fluids in porous media |
| JP2010051804A (ja) * | 2008-08-28 | 2010-03-11 | Siemens Ag | 測定シーケンスを形式的に記述して測定シーケンスのプログラミングを自動化するための磁気共鳴シーケンスモデル |
| US20110234220A1 (en) * | 2010-03-24 | 2011-09-29 | Jonathan Mitchell | System and method for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis |
| US20120153952A1 (en) * | 2010-12-16 | 2012-06-21 | General Electric Company | Method and system for gradient linear correction |
| US20120194188A1 (en) * | 2011-02-01 | 2012-08-02 | Zhe Phillip Sun | System and Method for Controlling Apparent Timing Dependencies for T2-Weighted MRI Imaging |
| US20140285196A1 (en) * | 2011-09-14 | 2014-09-25 | Petrochina Company Limited | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field |
| GB2505232A (en) * | 2012-08-23 | 2014-02-26 | Schlumberger Holdings | Magnetic resonance examination of porous samples |
| US20140253116A1 (en) * | 2013-03-05 | 2014-09-11 | Schlumberger Technology Corporation | System and method for obtaining nuclear magnetic resonance measurements on reservoir fluids for prediction of fluid properties |
| US9864033B1 (en) * | 2013-08-26 | 2018-01-09 | University Of New Brunswick | Free induction decay based magnetic resonance imaging methods |
| CN104280703A (zh) * | 2014-09-28 | 2015-01-14 | 中国石油大学(华东) | 一种基于脉冲梯度磁场的新型d-t2脉冲施加方法 |
| US20180364184A1 (en) * | 2017-06-16 | 2018-12-20 | Arcady Reiderman | Apparatus and method for nuclear magnetic resonance measurements on borehole materials |
| US20190346385A1 (en) * | 2018-05-11 | 2019-11-14 | Arcady Reiderman | Method and apparatus for nuclear magnetic resonance measurements on borehole materials |
Non-Patent Citations (3)
| Title |
|---|
| STEFAN A. HERTEL ET AL.: "Fast spatially-resolved T2 measurements with constant-gradient CPMG", 《MAGNETIC RESONANCE IMAGING》, vol. 56, 31 December 2019 (2019-12-31), pages 105 - 107 * |
| 章新友: "《医学图形图像处理 （新世纪第3版）》", 中国中医药出版社, pages: 105 - 107 * |
| 范宜仁;吴飞;李虎;霍宁宁;王要森;邓少贵;杨培强;: "D-T_2二维核磁共振脉冲序列及反演方法改进设计", 物理学报, no. 09, 8 May 2015 (2015-05-08), pages 524 - 539 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112834543B (zh) | 2024-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0098426B1 (en) | Method of eliminating effects of spurious free induction decay nmr signal caused by imperfect 180 degrees rf pulses | |
| FI78988B (fi) | Selektivt foerfarande och anordning foer utfoerande av lokaliserad nmr-spektroskopi. | |
| EP0086972B2 (en) | Method of NMR imaging which overcomes T2* effects in an inhomogeneous static magnetic field | |
| US20060097722A1 (en) | Downhole NMR flow and formation characterization while sampling fluids | |
| JPS5946546A (ja) | 核磁気共鳴による検査方法及び検査装置 | |
| GB2399642A (en) | NMR used to analyse formation fluids | |
| US10746827B2 (en) | MRI imaging system using permanent magnet array | |
| NO342087B1 (no) | Fremgangsmåte og apparat for NMR logging med godt signal/støyforhold | |
| US9389193B1 (en) | Spatially resolved magnetic resonance spin-spin relaxation distribution measurement methods | |
| CN107748126A (zh) | 一种获取岩石孔隙尺寸和孔隙表面弛豫率的核磁共振方法 | |
| De Deene | Review of quantitative MRI principles for gel dosimetry | |
| CN110029990B (zh) | 一种核磁共振测井方法和装置 | |
| Ino et al. | A compact SEOP 3He neutron spin filter with AFP NMR | |
| JPH01113033A (ja) | 核磁気分布決定方法 | |
| CN116930836B (zh) | 多核素同步一体化成像最佳脉冲功率测量方法和系统 | |
| CN106932426A (zh) | 一种测量分子中所有氢‑氢耦合常数的核磁共振方法 | |
| Liao et al. | NMR fluid analyzer applying to petroleum industry | |
| CN103645453A (zh) | 一种获取消除标量耦合调制的单体素一维定域谱的方法 | |
| CN111965719B (zh) | 一种弛豫时间测量方法和装置 | |
| CN112834543B (zh) | 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 | |
| US5317262A (en) | Single shot magnetic resonance method to measure diffusion, flow and/or motion | |
| CA2425254C (en) | Nmr spectroscopy using a gradient field nmr tool | |
| CN101545877B (zh) | 改进非均匀磁场中nmr波谱分辨率的方法和设备 | |
| CN108872291B (zh) | 快速定域多维拉普拉斯核磁共振方法与装置 | |
| US5235280A (en) | Method for determining optimized radio-frequency pulse shapes for selective excitation in magnetic resonance spectroscopy and imaging |
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 |