D-T2二维核磁共振脉冲序列及反演方法改进设计

面对日益复杂的勘探对象, D-T₂二维核磁共振技术在实际应用中面临无法兼顾扩散系数测量范围和横向弛豫分辨率的困境. 脉冲序列作为D-T₂二维核磁共振数据采集的核心技术, 其性能优劣直接影响应用效果, 在综合对比PFG, STE-PFG, BP-PFG、改良式CPMG, 扩散编程, 多回波间隔CPMG脉冲序列性能的基础上, 有效融合脉冲梯度场、恒定梯度场D-T₂脉冲序列的优点, 本文提出一种基于脉冲梯度场的双变量、两窗口D-T₂脉冲序列改进设计. 针对两个窗口的D-T₂二维核磁共振数据反演, 为突破现有反演方法无法兼顾反演精度和解谱效率的瓶颈, 本着第二个窗口回波信号为主、第一个窗口回波信号为辅的原则, 本文提出一种同时使用两个窗口数据参与解谱的联合TSVD反演方法. 气水、油水、稠油、油气水模型不同信噪比条件下的数值模拟结果表明, 本文提供的D-T₂改进脉冲序列达到了平衡扩散系数测量范围和横向弛豫分辨率的设计要求, 本文提供的联合TSVD反演方法也有效平衡了反演精度要求和解谱效率. 文中的D-T₂改进脉冲序列及联合TSVD反演方法在复杂油气藏流体识别和产能预测中具有广泛的应用前景, 可为促进国内D-T₂二维核磁共振岩心分析技术的发展提供有利条件.     

基于FPGA的低场核磁共振谱仪脉冲编程器设计与实现

脉冲编程器在核磁共振(NMR) 谱仪中非常重要,具有产生脉冲序列所需的射频脉冲,控制回波信号的同步接收等功能. 基于FPGA的同步性能和时序控制功能,设计了低场NMR谱仪的脉冲编程器,可产生任意脉冲序列,实现脉冲频率、相位、幅度的灵活调节. 利用设计的脉冲编程器在0.540 Tesla均匀磁场下进行了CuSO₄溶液的T₁和T₂测量,得到高信噪比的NMR回波串.     

基于蒙特卡罗算法的低场核磁共振测量效率提高方法

核磁共振（NMR）的纵向弛豫时间（T₁）、横向弛豫时间（T₂）、自扩散系数（D₀）,以及T₂-T₁、T₂-D₀测量目前广泛应用于石油测井行业.在测量D₀的SGSE序列中,通过逐渐增大90°和180°脉冲之间的时间间隔（T_d）,可以对液体扩散行为产生的影响进行调节.然而T_d的"起点"、"步进数"和"终点"等参数必须设置得当才能准确测量T₁和D₀.目前参数的设置依赖多次的人工调整和测量人员的经验,耗时且使用门槛较高.本文用蒙特卡罗方法进行大量随机模拟,根据前面若干点的测量结果筛选出满足要求的随机值,预测下一个测量点的位置.该算法可以实时更新参数设置,实现自动化测量,达到降低测量门槛、缩短测量时间的目的.经验证,该算法可以适用于T₁、D₀的测量.     

激光极化129Xe核自旋弛豫率的测量

由弱磁场下光泵自旋交换、强磁场中核磁共振测量方法,获得激光极化低压同位素¹²⁹Xe气体的核自旋弛豫率1/T₁及其与N₂气压力、环境温度的关系.结果表明:在291K的室温下,对于0~11 333.05Pa的N₂气压力范围,1/T₁随着N₂气压力的增加而增大;N₂气压力为1 466.63Pa时,1/T₁随着样品温度从234K到292K增高而减小     

一种改进的核磁共振测井测量脉冲序列

提出了一种可以提高核磁共振测井测量速度的新型测量脉冲序列,它基于对原有CPMG测量脉冲序列的改进.在原有的测量序列中,在一个CPMG脉冲序列结束之后,直接进行等待(T₁延时),以使磁化强度矢量能够完全恢复到平衡位置,然后再进行下一个CPMG脉冲序列的测量.由于测井仪器的连续运动,在测量含有长T₁成分的地层时,等待时间需要取得较长,使得测井速度很低.但在实际地层中,恰恰又是油气的T₁较长,这无疑给核磁共振测井带来了麻烦.作者发现,如果在一个CPMG脉冲序列结束之后在-x方向先施加一个90°脉冲(或在x方向施加一个270°脉冲),再进行等待,就能更好地使长T₁成分得到极化,从而缩短等待时间、提高测量速度.文中对这种改进后的脉冲序列进行了计算和分析,证明它确实具有一定的优越性.     

用于土壤水分测量的磁共振射频线圈宽频匹配方法

准确检测土壤水分对于改善农业用水管理至关重要.为了实现土壤水分的准确测量,本文设计了一种微型inside-out核磁共振（NMR）传感器（测量区域位于传感器周围）.其中磁体由三个轴向磁化的圆柱磁体组成,射频线圈由四个相同的弧面螺旋线圈串联组成.通过调整圆柱体之间的半径和高度,获得了一个圆周对称恒定梯度场.并利用宽频带匹配方法,增加射频线圈的激励区域,从而提高整体信噪比.将本文设计的传感器用于测量含水量不同的小颗粒沙土和大颗粒石英砂土壤模型样品的CPMG（Carr-Purcell-Meiboom-Gill）信号,通过逆拉普拉斯变换获得了T₂分布谱,并将T₂谱峰面积与称重法获得的含水量进行了相关性分析,探讨了该传感器用于测量土壤水分的可行性.     

回波间隔对核磁共振表观孔隙度的影响及矫正方法

本文对具有特定横向弛豫时间（T₂）的硫酸铜溶液进行了多回波间隔（T_E）的核磁共振（NMR）实验,并利用数值模拟对32组具有不同弛豫分量的模型进行了变T_E模拟实验,定量研究了T_E对致密油气、页岩气等低孔低渗储层NMR孔隙度的影响规律.实验结果表明,随着T_E的增大,各T₂弛豫组分NMR孔隙度先维持在100%左右,然后迅速衰减,当T_E增加到一定数值时,趋近于0;不同T₂弛豫组分NMR孔隙度开始迅速衰减及最后变为0的T_E值存在显著差异.根据不同T₂弛豫组分NMR孔隙度与T_E的关系,将整个NMR测量分为无损测量区、快速衰减区、无效参数区和仪器盲区4个区域.对特定弛豫组分而言,在快速衰减区弛豫组分损失量与T_E呈对数关系,本文还给出了该区域NMR孔隙度的校正公式及方法.     

基于差分进化算法的核磁共振T2谱多指数反演

提出了一种基于差分进化(DE)算法的核磁共振弛豫信号多指数反演新方法. 将核磁共振T₂谱反演问题转化为带非负约束的非线性优化问题,不需要预先给定横向弛豫时间T_{2分布,直接利用差分进化算法进行反演计算. 在测量信号低信噪比的情况下,计算机模拟和实验数据反演都表明了该方法在分析处理NMR弛豫信号中的有效性. 关键词： 核磁共振 多指数反演 差分进化 岩心分析}     

非自治物质畸形波的传播操控

研究了(1+1)维的变系数Gross-Pitaevskii方程, 获得了该方程的精确畸形波解. 基于该精确畸形波解, 深入研究了非自治物质畸形波在随时间指数变化的相互作用下的传播动力学行为, 发现非自治畸形波除具有“来无影、去无踪”的不可预测特性外, 也可实现完全激发、抑制激发以及维持激发等操控. 研究表明, 畸形波操控的关键是对累积时间的最大值T_max 与峰值位置T₀ (或T_I,T_II)值大小关系的调节. 当T_max > T₀ (或T_I,T_II)时畸形波被快速地完全激发, 热原子团中的原子增加到凝聚体中. 当T_max = T₀ (或T_I,T_II) 时畸形波激发到最大振幅, 可以维持相当长的时间而不消失, 热原子团中的原子增加到凝聚体中. 当T_max < T₀ (或T_I,T_II)时畸形波没有充足的时间来激发而被抑制甚至消失, 凝聚体中的原子减少. 这些结果在理论和实际应用上具有启迪意义.     

脉冲梯度场增强的DQF-COSY实验t1噪声研究

脉冲相位的不稳定性是导致二维谱中的t₁噪声的原因之一.本文从理论上和实验上对脉冲相位不稳定性所引起的t₁噪声进行了研究.从角动量算符I_x,I_y和I_z推导中发现,由于磁化矢量在P型和N型DQF-COSY实验的脉冲梯度场中进动方向的内在差异,一部分脉冲相位误差在P型DQF-COSY实验会被重聚掉,在N型DQF-COSY实验中却会被加重.因此P型DQF-COSY实验的t₁噪声比N型小.理论推导还得到t₁噪声在P型和N型中的理论比值为√3:√11.以上理论结果得到我们实验结果的支持.     

Selection of pulse sequences producing maximum tissue contrast in magnetic resonance imaging

The importance of spin density [N(H)] and spin-lattice (T₁) and spin-spin (T₂) relaxation in the characterization of tissue by nuclear magnetic resonance (NMR) is clearly recognized. This work considers which optimized pulse sequences provide the best tissue discrimination between a given pair of tissues. The effects of tissue spin density and machine-imposed minimum rephasing echo times (TE_MIN) for achieving maximum signal tissue contrast are discussed. A long TE_MIN sacrifices T₁-dependent contrast in saturation recovery (SR) and inversion recovery (IR) pulse sequences so that spin-echo (SE) becomes the optimum sequence to provide tissue contrast, due to T₂ relaxation. Pulse sequences providing superior performance may be selected based on spin density and T₁ and T₂ ratios for a given pair of tissues. Selection of the preferred pulse sequence and interpulse delay times to produce maximum tissue contrast is strongly dependent on knowledge of tissue spin densities as well as T₁ and T₂ characteristics. As the spin density ratio increases, IR replaces SR as the preferred sequence and SE replaces IR and SR as the pulse sequence providing superior contrast. To select the optimal pulse sequence and interpulse delay times, an accurate knowledge of tissue spin density, T₁ and T₂ must be known for each tissue.     

在反转恢复测试中磁化矢量的演化特征:辐射阻尼效应

在强辐射阻尼存在下,水样(90% H₂O in D₂O)的反转恢复实验表明:当两脉冲的相位相差180°且反转脉冲角<180°时,或两脉冲相位一致但反转脉冲角>180°时,在检测期观测到的信号强度将不发生从负极大值到正极大值的突变;在同样的条件下,如果存在频率偏置,信号强度存在波动,即beating效应.只有当两脉冲的相位一致而反转脉冲角<180°时,或两脉冲的相位相差180°但反转脉冲角>180°时,在检测期信号强度才发生突变,即jumping效应.这些现象都可通过辐射阻尼理论予以合理地解释.另外,在检测期当磁化矢量运动到-z轴附近(对应于τ=T_rdln{tan[(π±δ)/2]}),信号强度与理论预计的偏差实际上与T₁弛豫效应有关.     

Myocardial proton spin-lattice relaxation times in vitro: Effect of elapsed time after excision

An assumption made in using excised tissue for in vitro nuclear magnetic resonance (NMR) studies is that variables of interest, such as spin-lattice (T₁) relaxation times, remain stable for periods of time after excision sufficient to perform NMR spectroscopy. In this study, we evaluated the changes in T₁ of rat myocardium, measured at two NMR field strengths, at serial time intervals up to 72 hours postmortem. Left ventricular myocardium from six male Sprague-Dawley rats was excised and stored at room temperature in sealed NMR sample tubes. Spin-lattice relaxation times were determined with a modified inversion-recovery pulse sequence immediately postmortem and at intervals up to 72 hours post-excision; NMR studies were performed using 90 MHz and 360 MHz spectrometers. A gradual decrease in T₁ was noted with increasing time post-excision; T₁ was not significantly shorter than baseline until 72 hours postmortem at either field strength. The rate of change of T₁ was similar at the two field strengths. At any given time post-excision, T₁ was significantly higher (p < 0.001) at 360 MHz than at 90 MHz. We conclude that, with proper tissue handling and storage techniques, rat myocardial T₁ is stable postmortem sufficiently long to permit meaningful NMR studies of excised tissue.     

Minimization of diffusive attenuation in T2-weighted NMR images of porous solids using turboSPI

This communication presents the results of T₂-weighted nuclear magnetic resonance imaging (MRI) of a water saturated porous volcanic rock using a fast single point imaging technique (turboSPI). Imaging porous materials with heterogeneously distributed mineral products, and air voids, using conventional imaging methods, which introduce T₂-weighting by increasing the time between the excitation and refocusing pulses, often results in high diffusive signal losses and susceptibility distortion. T₂-weighted images acquired of a water saturated porous rock using turboSPI with an effective echo time of 8.1 ms (actual inter-echo time of 0.9 ms) exhibit significantly decreased diffusive attenuation, compared to analogous images obtained with an inter-echo time of 8.1 ms.     

NMR imaging of white button mushroom (Agaricus bisporis) at various magnetic fields

Nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) have been applied to visualize physiological phenomena in plants and agricultural crops. Imaging sequences that result in contrast of a combination of parameters (e.g., proton density, ) cannot be used for a correct and unique interpretation of the results. In this study multiecho imaging together with monoexponential T₂ decay fitting was applied to determine reliable proton density and T₂ distributions over a mushroom. This was done at three magnetic field strengths (9.4, 4.7, and 0.47 T) because susceptibility inhomogeneities were suspected to influence the T₂ relaxation times negatively, and because the inflences of susceptibility inhomogeneities increase with a rise in magnetic field strength. Electron microscopy was used to understand the different T₂'s for the various tissue types in mushrooms. Large influences of the tissue ultrastructure on the observed T₂ relaxation times were found and explained. Based on the results, it is concluded that imaging mushrooms at low fields (around or below 0.47T) and short echo times has strong advantages over its high-field counterpart, especially with respect to quantitative imaging of the water balance of mushrooms. These conclusions indicate general validity whenever NMR imaging contrast is influenced by susceptibility inhomogeneities.     

Large angle spin-echo imaging

A study was undertaken to assess the use of excitation flip angles greater than 90° for T₁ weighted spin-echo (SE) imaging with a single 180° refocusing pulse and short TR values. Theoretical predictions of signal intensity for SE images with excitation pulse angles of 90–180° were calculated based on the Bloch equations and then measured experimentally from MR images of MnCl₂ phantoms of various concentrations. Liver signal-to-noise ratios (SNR) and liver-spleen contrast-to-noise ratios (CNR) were measured from breathhold MR images of the upper abdomen in 16 patients using 90 and 110° excitation flip angles. The theoretical predictions showed significant improvements in SNR with excitation flip angles >90°, which were more pronounced at small TR values. The phantom studies showed reasonably good agreement with the theoretical predictions in correlating the excitation pulse angle with signal intensity. In the human imaging studies, the 110° excitation pulse angle resulted in a 7.4% (p < .01) increase in liver SNR and an 8.2% (p = .2) increase in liver-spleen CNR compared to the 90° pulse angle at TR = 275 ms. Increased signal intensity resulting from the use of large flip angle excitation pulses with a single echo SE pulse sequence was predicted and confirmed experimentally in phantoms and humans.     

Partial angle inversion recovery (PAIR) MR imaging: spin-echo and snapshot implementation.

The effects of varying the inversion or excitation RF pulse flip angles on image contrast and imaging time have been investigated in IR imaging theoretically, with phantoms and with normal volunteers. Signal intensity in an IR pulse sequence as a function of excitation, inversion and refocusing pulse flip angles was calculated from the solution to the Bloch equations and was utilized to determine the contrast behavior of a lesion/liver model. Theoretical and experimental results were consistent with each other. With the TI chosen to suppress the fat signal, optimization of the excitation pulse flip angle results in an increase in lesion/liver contrast or allows reduction in imaging time which, in turn, can be traded for an increased number of averages. This, in normal volunteers, improved spleen/liver contrast-to-noise ratio (9.0 vs. 5.7, n = 8, p less than 0.01) and suppressed respiratory ghosts by 33% (p less than 0.01). Reducing or increasing the inversion pulse from 180 degrees results in shorter TI needed to null the signal from the tissue of interest. Although this decreases the contrast-to-noise ratio, it can substantially increase the number of sections which can be imaged per given TR in conventional IR imaging or during breathold in the snapshot IR (turboFLASH) technique. Thus, the optimization of RF pulses is useful in obtaining faster IR images, increasing the contrast and/or increasing the number of imaging planes.     

Novel MR image contrast mechanisms in epilepsy

A range of magnetic resonance (MR) parameters are introduced, which can give rise to image contrast by using suitable pulse sequences, and that can be measured quantitatively. Their relationship to tissue pathology is given as far as possible. Techniques for their measurement, and results from multiple sclerosis, stroke, and epilepsy are given. The parameters are proton density, T₁, T₂, transverse magnetisation decay, which gives estimates of extracellular water and myelin concentrations, magnetisation transfer ratio and T_1sat, and diffusion (including trace and anisotropy measured from the tensor matrix).     

基于T2截止值的幂函数构建毛管压力曲线

岩心核磁共振（NMR）T₂谱和毛管压力曲线都在一定程度上反映了岩石孔隙结构，理论分析表明可利用T₂谱构建毛管压力曲线，由此快速获取储层孔隙结构的信息．本文对18块岩样T₂分布和毛管压力曲线进行了分析，提出将T₂截止值作为幂函数的分段点，采用分段幂函数方法构建岩样的NMR毛管压力曲线，并与采用不分段幂函数方法获得的毛管压力曲线进行了对比．研究结果表明，与不分段幂函数方法相比，分段幂函数方法得到的结果和压汞实验测定的毛管压力曲线吻合度更高，平均拟合度（R²）达到0.943 1，并且论证了T₂截止值作为分段点进行分段幂函数法构建NMR毛管压力曲线的合理性和可靠性．T₂截止值的引入提高了幂函数法构建NMR毛管压力曲线的精度，是利用NMR T₂谱构建岩样毛管压力曲线的有价值探索．     

Nuclear magnetic resonance relaxometry of the normal heart: Relationship between collagen content and relaxation times of the four chambers

The use of nuclear magnetic resonance (NMR) relaxation time measurements for characterization of abnormal cardiac tissue depends upon knowledge of variations of relaxation times of normal myocardium and determinants of these variations. We calculated in vitro NMR T₁ and T₂ relaxation times of canine myocardium from the four cardiac chambers, and determined hydroxyproline concentration (as a measure of collagen) and percent water content of the samples. We found both water content and T₁ relaxation time of the right ventricle to be significantly greater than the left atrium (p < 0.05). T₂ relaxation time of the left ventricle was found to be shorter than each of the other three chambers (p < 0.05). There were significant correlations between the spin-lattice relaxation time and both percent water content (r = 0.58) and hydroxyproline concentration (r = 0.45). A significant correlation was also found between T₂ relaxation time and hydroxyproline concentration (r = 0.49). When T₁ and T₂ were adjusted for water and hydroxyproline content, there was no longer any evidence for significant interchamber differences for either T₁ or T₂. These data suggest that differences in NMR relaxation times exist among the four chambers of the normal canine heart. Furthermore, a major determinant of myocardial spin-lattice relaxation time is tissue water content while both collagen content and percent water content significantly contribute to variability in cardiac chamber T₂ relaxation times.