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

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

一种基于组分补偿的二维核磁共振测井数据高精度处理方法

二维核磁共振技术能够对储层中各类含氢流体进行无损、快速、定量的测量和表征，但受限于采集方式和参数，核磁共振设备在对页岩油等致密储层中的有机质、沥青等超快弛豫组分进行检测时，经常出现由于信号采集不完整所导致的二维谱中流体组分缺失或不准的问题.本文提出了基于超快弛豫组分补偿技术的T₂-T₁二维谱高精度反演方法，该方法将一维核磁共振前端信号补偿技术进行推广，通过在二维核磁数据反演前对回波数据进行组分补偿，能够有效解决二维核磁共振测井前端信号漏失的问题.实验及测井数据的应用表明，该方法在页岩油等富含快弛豫组分信号的储层中，可以得到更加精准和完整的储层信息.     

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

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

2D NMR技术在石油测井中的应用

近几年,2D NMR技术得到迅速发展,特别是在核磁共振测井领域. 该文将主要介绍2D NMR技术的脉冲序列、弛豫原理以及2D NMR技术在石油测井中应用. 2D NMR技术是在梯度场的作用下,利用一系列回波时间间隔不同的CPMG脉冲进行测量,利用二维的数学反演得到2D NMR. 2D NMR技术可以直接测量自扩散系数、弛豫时间、原油粘度、含油饱和度、可动水饱和度、孔隙度、渗透率等地层流体性质和岩石物性参数. 从2D NMR谱上,可以直观的区分油、气、水,判断储层润湿性,确定内部磁场梯度等. 2D NMR技术为识别流体类型提供了新方法.     

超快速COSY类2D NMR实验弛豫假峰的抑制

在脉冲梯度场基础上研究了进一步缩短COSY类2D NMR实验数据采集时间的新方法.利用主动屏蔽梯度场线圈的梯度场恢复时间仅为100μs数量级这一特点,提出一类基于脉冲梯度场消除弛豫假峰的脉冲序列.并以DQF-COSY脉冲序列为例加以验证.实验结果表明:采用文中的方法,取序列重复时间远小于T₁和T₂是切实可行的;梯度场技术能很好地抑制弛豫假蜂以达到2D NMK数据的快速采集.     

基于改进非线性拟合的核磁共振T2谱多指数反演

核磁共振T₂谱多指数反演算法是开展复杂体系样品核磁共振(NMR)弛豫研究最重要的数学工具. 常用的T₂谱多指数反演算法一般都是事先给出弛豫时间T₂分布的布点, 然后转化为线性拟合问题进行求解. 在求解的T₂谱较为分散的时候, 反演得到的T₂谱精确度不高, 分辨率较低. 非线性拟合是解决这个问题的有效办法. 本文针对分散T₂谱反演利用非线性拟合时遇到的初值依赖及运算复杂问题, 利用线性回归最小二乘方法, 改进了其中的带非负约束非线性优化模型, 将搜索的反演参数从T₂, f 减少为T₂, 加快了收敛速度, 减少了对初值的依赖, 提高了反演精度, 使算法更加稳健. 通过用改进的Levenberg-Marquardt算法和差分进化算法进行计算机模拟反演及实验数据反演, 验证了改进方法在核磁共振T₂ 谱反演中的有效性.     

优化重聚脉冲提高梯度场核磁共振信号强度

缩短射频脉冲宽度, 有助于解决脉冲电力消耗大、样品吸收率高、信噪比低等极端条件核磁共振探测的关键问题. 本文首先分析射频脉冲角度对核磁共振自旋回波信号强度的影响机理, 基于Bloch方程推导了回波信号幅度与扳转角、重聚角的关系. 在特制核磁共振分析仪上采用变脉冲角度技术, 分别在均匀磁场和梯度磁场条件下实现对扳转角和重聚角与回波信号强度关系的数值模拟和实验测量. 结果表明, 梯度场中, 扳转角为90°、重聚角为140°的射频脉冲组合获得最大首波信号强度, 比180°脉冲对应的回波幅值提高13%, 能耗降低至78%. 选用该重聚角(140°) 优化设计饱和恢复脉冲序列探测流体的纵向弛豫时间T₁特性, 准确获得 T₁分布.该结果对于低电力供应、且对信噪比有较高要求的核磁共振测量, 如随钻核磁共振测井和在线核磁共振快速检测等, 具有重要意义. 关键词： 核磁共振 信号强度 重聚脉冲角度 Bloch方程     

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

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

低场核磁共振扩散-横向弛豫二维反演算法

针对低场核磁共振一维反演中无法分辨一维谱中重叠组分和目前报道的扩散-横向弛豫二维反演算法计算量大、计算耗时长的问题,提出了一种计算量小、计算效率高、耗时短的扩散-横向弛豫二维反演算法. 首先对扩散系数D-和横向弛豫时间T₂进行布点;其次根据信号采集条件计算出两个核心矩阵,并分别进行奇异值分解;然后,由所采集信号计算出两个核心矩阵的奇异值截断值,分别对两个核心矩阵的奇异值矩阵进行截断并求其逆矩阵;最后计算出初始反演结果,并添加非负约束经过多次迭代得到最终反演结果. 实验结果证明,提出的扩散-横向弛豫二维反演算法在不影响反演结果准确性的基础上,能极大提高计算效率.     

基于ADI-FDTD模拟的岩石核磁共振横向弛豫特征分析

将复杂的骨架-孔隙系统抽象成等效双孔介质,根据Bloch方程构建数学模型,用交替隐式时域有限差分(ADI-FDTD)和联合反演迭代法(SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振T₂谱的反演,定量研究扩散系数、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响.结果表明:横向宏观磁化矢量衰减速率与扩散系数和微孔隙分量成正比,与孔隙宽度成反比,与表面弛豫速率基本无关.当扩散系数较大、孔隙宽度较小时,核磁共振T₂谱难以直观反映孔隙组分及孔隙结构.应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响.     

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.     

Diffusion of W on A W(211) plane

The diffusion of W on a (211) plane of a W field emitter has been re-examined by means of the fluctuation autocorrelation method. Diffusion along channels yielded E = 16.8 ± 0.5 kcal, D₀ = (3 ± 1) × 10⁻⁵ cm² s⁻¹. For diffusion across channels E =6.6 kcal, D₀ = 4 × 10⁻⁹cm² s⁻¹ at T < 752 K, and E = 24 kcal, D₀ = 5 × 10⁻⁴ cm² s⁻¹ at T > 752 K. The results for diffusion along channels yield E and D₀ values intermediate between recent results by Wang and Ehrlich [Surf. Sci. 206 (1988) 451] using field ion microscopy (E = 19 kcal, D₀ = 7.7 × 10⁻³ cm² s⁻¹) and Tringides and Gomer [J. Chem. Phys. 84 (1986) 4049], using the same method as the present work but a larger slit (E = 13.3 kcal, D₀ = 7 × 10⁻⁷ cm² s⁻¹). The results for cross channel diffus good agreement with those of Tringides and Gomer below 752 K, where these authors stopped. The new high temperature results suggest that the channel wall exchange mechanism postulated by Tringides and Gomer for cross channel diffusion at low T gives way to diffusion by climbing over the channel walls with higher E but also higher D₀ above 752 K. Possible reasons for the discrepancies between these three sets of results and the absence of cross channel diffusion in the work of Wang and Ehrlich are briefly discussed.     

Investigation on the Cs 6S1/2 to 7D electric quadrupole transition via monochromatic two-photon process at 767 nm

We experimentally demonstrate the cesium electric quadrupole transition from the 6S_1/2 ground state to the 7D_3/2,5/2 excited state through a virtual level by using a single laser at 767 nm. The excited state energy level population is characterized by varying the laser power, the temperature of the vapor, and the polarization combinations of the laser beams. The optimized experimental parameters are obtained for a high resolution transition interval identification. The magnetic dipole coupling constant A and electric quadrupole coupling constant B for the 7D_3/2,5/2 states are precisely determined by using the hyperfine levels intervals. The results, A = 7.39 (0.06) MHz, B = −0.19 (0.18) MHz for the 7D_3/2 state, and A = −1.79 (0.05) MHz, B =1.05 (0.29) MHz for the 7D_5/2 state, are in good agreement with the previous reported results. This work is beneficial for the determination of atomic structure information and parity non-conservation, which paves the way for the field of precision measurements and atomic physics.     

Boron diffusion in bcc-Fe studied by first-principles calculations

The diffusion mechanism of boron in bcc-Fe has been studied by first-principles calculations. The diffusion coefficients of the interstitial mechanism, the B–monovacancy complex mechanism, and the B–divacancy complex mechanism have been calculated. The calculated diffusion coefficient of the interstitial mechanism is D_0= 1.05 ×10~(-7)exp(-0.75 e V/k T) m~2· s~(-1), while the diffusion coefficients of the B–monovacancy and the B–divacancy complex mechanisms are D_1= 1.22 × 10~(-6)f1exp(-2.27 e V/k T) m~2· s~(-1)and D_2≈ 8.36 × 10~(-6)exp(-4.81 e V/k T) m~2· s~(-1), respectively. The results indicate that the dominant diffusion mechanism in bcc-Fe is the interstitial mechanism through an octahedral interstitial site instead of the complex mechanism. The calculated diffusion coefficient is in accordance with the reported experiment results measured in Fe–3%Si–B alloy(bcc structure). Since the non-equilibrium segregation of boron is based on the diffusion of the complexes as suggested by the theory, our calculation reasonably explains why the non-equilibrium segregation of boron is not observed in bcc-Fe in experiments.     

Diffusion of silicon in Fe-based amorphous and crystalline alloys

The segregation rate of silicon was measured in three different Fe-based amorphous and crystallized alloys with different silicon contents: 3.5, 5 and 9 at%. Analysis of the segregation kinetics yielded the diffusion activation energies E, as well as the frequency factors D₀. A linear dependence was found between In D₀ and E. In general, the D₀ and E values were lower for the amorphous specimens than for the crystalline ones, were independent of silicon content and are explained in terms of an oversaturated concentration of structural defects. In the crystalline specimens, diffusion behaviour was influenced by silicon content near the solubility limit of silicon in iron.     

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.     

NMR microimaging studies of water diffusivity in saturated, microporous systems

A novel microimaging method is described for measuring the effective diffusion coefficient (D_∞) of water in saturated, microporous materials. A simple multicompartment computer simulation was used to show how the value of D_∞, when combined with measurements of water proton relaxation times and pulsed field gradient spin-echo amplitudes, allows the pore size distribution and pore connectivity to be characterized.     

Diffusion imaging with paired CPMG sequences

We report a method for mapping apparent diffusion coefficients using two interleaved CPMG sequences. Image slice selection is performed prior to each sequence, allowing the use of non-selective “hard” refocussing pulses. Phase cycling of the slice-selection process for each projection minimizes out-of-slice contributions to echo amplitudes. This permits an accurate evaluation of both T₂ and apparent proton self-diffusion coefficients.