与样本间有相对移动的核磁共振装置用的核磁共振脉冲序列.pdf

本发明公开了一种进行核磁共振测量的方法和装置。在测量过程中NMR测量装置、所测量的样本或这两者都可能移动。依据包络线来调制RF载波信号的包络线以产生一第一RF脉冲序列。在发射第一序列期间可以改变包络线、RF信号的相位和/或静态磁场以使样本的第一区域基本饱和。发射一第二RF脉冲序列以在第一区域内建立一共振区域，并测量样本的特性。

1.  一种用于井孔的NMR测量装置，在该装置和地层构造的样本之间有相对移动，包括：
a)至少一磁体；
b)至少一线圈；和
c)与所述线圈相连的脉冲发生器，该脉冲发生器适于：
i)发射第一RF脉冲序列，以便产生一个时变磁场，所述第一序列包括至少一个再聚焦脉冲，以便产生来自样本的一个共振区域的至少一个回波；
ii)建立一个静态磁场，以便建立所述共振区域；以及
iii)在所述第一序列的发射期间瞬间改变所述静态电场至少一次，以便使比所述共振区域更大的区域饱和。

2.  如权利要求1所述的装置，其中改变所述静态磁场的步骤进一步包括发射至少一个其它脉冲。

3.  如权利要求2所述的装置，进一步包括当所述RF脉冲序列未被发射时同步发射至少一个其它脉冲。

4.  如权利要求2所述的装置，进一步包括当至少一个RF脉冲发射时同步发射至少一个其它脉冲。

5.  如权利要求1所述的装置，其中饱和步骤包括发射附加RF脉冲，以便使比所述共振区域更大的区域饱和。

6.  如权利要求1-5所述的装置，其中由于NMR装置的移动产生所述相对移动。

7.  如权利要求1-5所述的装置，其中由于样本的移动产生所述相对移动。

8.  如权利要求1-5所述的装置，进一步包括产生轴向对称的静态磁场的步骤，该轴向对称的静态磁场在共振区域中具有等值线，该等值线在与NMR装置的纵轴基本一致的方向上基本是直线。

9.  如权利要求1-5所述的装置，进一步包括在地层构造中钻井孔的同时测量样本的特性。

与样本间有相对移动的核磁共振装置用的核磁共振脉冲序列
本申请是申请号为200410069936.6、申请日为1999年12月6日的、发明名称为“与样本间有相对移动的核磁共振装置用的核磁共振脉冲序列”的发明专利申请的分案中请。
技术领域
本发明一般涉及对在核磁共振(NMR)区域附近的自旋的预处理。
背景技术
核磁共振(NMR)可以用来测定比如形体组织(为医学成像的目的)或地下地层构造(为井下记录的目的)的样本的特性。例如，对于地下地层构造，可以应用NMR来测定并绘制地层构造的孔隙度、地层构造类型、渗透率和油量。
参考附图1，作为一个例子，在边记录边钻井(LWD)操作中应用NMR来绘制地下地层构造10的特性图。这里，成轴对称的NMR测井仪6是钻杆柱5的一部分，应用该钻杆柱5在地层10上钻一井孔3。例如，测井仪6可以是Sezginer等人在美国专利No.5,705,925(1998年1月6号授权，题为“Pulsed Nuclera Magnetism Tool For FormationEvalution While Drilling Including a Shortened or Truncated CPMGSequence”)中或Miller在美国专利No.5,280,243(1994年1月18号授权，题为“System For Logging a Well During the DrillingThereof”)中所描述的测井仪之一。
同其它大多数井下地层构造测量方法相比，NMR测量方法具有两个明显的特征。首先，来自地层的NMR信号是来自一较小的共振空间，比如通常为一薄的共振空间20a(参见附图2)，共振空间20a的径向厚度与磁场(未示)的大小成比例。例如，取决于不同的共振区域的形状，该空间可以在一个方向小到1毫米(mm)而在另一个方向延伸长达几英尺。第二，NMR测量可以不是瞬时的。如下文进一步描述，结合这两点事实使NMR测井仪易于进行移动测量，比如NMR测井仪6围绕井孔3的周围移动。
如下文所述，为了进行NMR测量，NMR测井仪6包括建立称为的静态磁场(未示)的永磁体、发射垂直于磁场的时变磁场的射频(RF)线圈或天线和接收来自与NMR测量对应的地层构造的自旋回波的RF线圈或天线。并且这两个线圈可以结合在一单个的发射/接收天线中。
作为一个例子，通过发射NMR检测序列以使原子核产生自旋回波，NMR测井仪6测量地层构造10的氢原子核的自旋-自旋驰豫时间T2。依次分析自旋回波得到时间T2的分布，从这一分布中可以得到地层构造的特性。例如，附图4中所描述的Carr-Purcell-Meiboom-Gill(CPMG)序列15就是这样的一种NMR检测序列。应用序列15可以得到时间T2的分布，并可以应用这一分布来测定并绘制地层构造10的特性。
应用CPMG序列15来测量时间T2的技术包括下面的步骤。在第一步中，NMR测井仪6以合适的时间间隔发射磁场以应用90°激励脉冲14a来使氢原子核的自旋(最初与磁场方向一致)旋转90°。虽然没有图示，实际上每个脉冲都是RF载波信号的包络线或脉冲串。在自旋从磁场方向旋转90°后，首先自旋立即开始在与磁场垂直的平面中一致地进动，接着逐渐地失去同步。在第二步中，在NMR脉冲14a随后的一固定时间T中，NMR测井仪6以脉冲的形式施加磁场并持续一较长时间段(与NMR脉冲14a比较)，以应用NMR再聚焦脉冲14b来旋转进动的自旋一个附加角度180°，同时使其载波相位移动±90°。NMR脉冲14b使自旋重新同步，并使自旋在180°再聚焦脉冲14b后发射一辅助自旋回波(见附图5)，该自旋回波的峰值时间大约等于T。以间隔t_e(大致为2·T)重复步骤2“k”次(这里“k”称为回波数，例如，可以假定是从几百到多达几千中的任何一个值)。对于第三步骤，在完成自旋回波序列后，在聚集另一序列自旋回波的下一CPMG序列15开始之前要求经过一等待周期(通常称为等待时间)，以使自旋返回到沿磁场的平衡状态。观测每个自旋回波序列的衰减并用其来导出T2的分布。
时间T2^*是在应用90°激励脉冲14a后使自旋具有不再一致地进动的特点的时间。这里，在90°激励脉冲14a结束时，所有的自旋都指向与静态磁场垂直的一个共同的方向，对于完全均匀的磁场自旋以称为拉莫频率的共振频率进动。拉莫频率可以表述为这里γ为旋磁比，一个核常数。然而，一般地磁场并不是均匀的，在激励之后，由于静态磁场的非均匀性使自旋具有在T2^*内产生的相移。这种衰减是可逆的，但产生回波的再聚焦脉冲14a可以使这种衰减反相。此外，还产生不可逆的相移(自旋-自旋驰豫)，并可通过T2时间常数来描述这种相移。这就导致依据T2时间常数在CPMG序列中的连续回波幅值的衰减。一般地，应用“内检外测型”NMR测量T2＞＞T2^*的自旋。
如上文所述，可以应用T2时间的分布来测定地层构造的特性。例如，参考附图6，地层可能包括容纳有约束流体的小孔和容纳有自由的、可采出的流体的大孔。可以应用分离边界时间T2(在附图6中被称为T_CUT-OFF)来将T2分布分为两部分：一部分包括比T_CUT-OFF小的时间(表示约束流体)和一部分包括比T_CUT-OFF大的时间(表示自由的、可采出的流体)。
一般地，通过观测由特定的CPMG序列15产生的自旋回波16的衰减来计算时间T2。不幸的是，钻杆柱5(参见附图1)可能要经过剧烈的横向运动。然而，时间T2大致与另一被称为自旋-晶格松驰时间T1的时间常数成比例。时间T1具有使自旋沿磁场方向重新返回到平衡状态的特性，由此，考虑两T1和T2时间，可以认为每个自旋在T1恢复过程中以极紧密的螺旋运动形式朝平衡位置运动回来。幸运的是，时间T1和T2大致成比例。因此，可以从测量的时间T1导出时间T2的分布。事实上，建立约束流体截止点的最初的工作是应用T1来完成的。然后根据T2将这些结果表述出来并应用在商业上。参见W.E.Kenyon，J.J.                                      Howard，A.Sezginer，C.Straley，A.Matteson，K.Horkowtitz，and R.Ehrlich，等人的Pore-Size Distribution and NMR in Microporous ChertySandstones，论文LL(在第30届年度记录研讨会(the 30th AnnualLogging Sympos ium)上提交的论文，SWPLA，1989年6月11-14日)。
基于极化的测量可以应用反向恢复序列(inversion recoverysequence)或饱和恢复序列(saturation recovery sequence)。在饱和恢复序列情况下，例如应用几个将磁化强度减少为零的90°脉冲使自旋系统饱和。然后在应用监测脉冲或脉冲序列(比如CPMG序列)之前使自旋系统恢复达一可变长度的时间。反向恢复技术要求在原子核本身的自旋已经与静态磁场成一线之后应用180°脉冲使自旋反向。随着时间的过去，自旋按照T1朝平衡方向衰减，但由于180°脉冲在检测器中并不产生信号，并没有进行测量。然而，在衰减完成之前，监测脉冲或脉冲序列(比如CPMG序列)中断了衰减，该监测脉冲或脉冲序列将自旋旋转到测量平面(即，在检测器中产生信号)。有意义的信息就是在最初的90°“读出”脉冲瞬间之后的信号幅值。这种幅值完全依赖于在初始的180°脉冲和90°脉冲之间的恢复时间。接下来测定幅值，允许自旋系统完全驰豫到平衡状态，然后重复脉冲序列。
作为反向恢复序列在井下应用的例子描述在Kleinberg等人的美国专利No.5,023,551(题为“Nuclear Magnetic Resonance PulseSequences For Use With Borehole Logging Tools”，1991年6月11日授权)中。然而描述在美国专利’551中的反向恢复序列并不能应用绝热脉冲，由此导致仅在一个狭窄的区域中勘测。此外，在存在移动和“内检外测”的情况下，使一个区域饱和比使其完全反相恢复要容易。因此，比较可取的做法是使区域饱和。
回过来参考附图2，一般地，应用基于极化的测量而不应用上文所描述的基于衰减的测量来测量时间T1。在基于极化的测量中，每个基于极化的测量首先包括应用饱和序列来使在共振区域(比如，在附图2中所描述的圆柱形共振空间20a)中的自旋到达饱和的步骤。接着，经过一极化时间将共振空间20a极化到静态磁场。随后，应用检测序列(比如CPMG序列)以产生来自地层构造10的自旋回波。然后，分析开始的少数几个自旋回波的幅值以确定孔隙度分布Φ(T1)的极化加权整数Φ(t_wait)。因为，只需观测开始的少数几个回波就可以确定信号的幅值，故可以在比基于衰减测量T2所需的时间更短的持续时间内测量T1，因而，测量T1时更不适合于移动NMR测井仪6。以不同的等待时间连续地重复(在适当的饱和序列之后)检测序列几次以得到孔隙度分布Φ(T1)。
作为一个例子，应用基于极化的测量来测量在位于饱和空间20b中的共振空间20a(参见附图2)中的氢原子核的时间T1。在这里，NMR测井仪6首先使在饱和空间20b中的自旋饱和。然而，极化时间应该足够长以允许NMR测井仪6能够在井孔内作充分的移动。在这种情况下，NMR测井仪6的移动使共振空间20a移动，以使NMR测井仪6接收来自移动了的共振空间20a′(参见附图3)中的自旋回波，该共振空间20a′部分落到最初的饱和空间20b之外。结果，移动了的共振空间20a′包括一没有饱和自旋的区域(一般称为“移入新自旋”效应)和一带有饱和自旋的最初的饱和空间20b的部分区域。不幸的是，基于极化的NMR技术不允许在极化时间中移入“新自旋”，因为新自旋可能引入测量误差。例如，测量可能错误地指示在地层构造中有比实际出现的更多的约束流体量。
在申请号为PCT/US97/23975的PCT申请(题为“Method ForFormation Evaluation While Drilling”，1997年12月29日申请)中描述了一种使更大区域饱和的一种方法。这一申请指出，在测量开始时，发射覆盖一个相对较大频率范围和/或特别宽的带宽的一个或多个射频脉冲或应用具有能够使NMR测井仪周围的圆柱形空间饱和的扫描频率的一个或多个脉冲。该申请进一步描述了当由于测井仪移动超出饱和区域引起测量无效时应用加速度峰值来测量。该申请还进一步描述在测井仪上装备支座以阻止测井仪移动超出饱和区域。
因此，人们一直渴望将在NMR测量装置和所勘测的样本之间的相对移动引起的误差最小化。
发明内容
本发明公开了一种供在装置和样本之间有相对移动的NMR测量装置使用的方法。装置、样本或这两者都可能移动。在本发明一个实施例中，该方法包括发射第一RF脉冲序列。第一序列具有包络线。在发射第一序列期间改变包络线，以便样本的第一区域基本饱和。发射第二RF脉冲序列以在第一区域内建立共振区域，并测量样本的属性。
在另一实施例中，供在装置和样本之间有相对移动的NMR测量装置使用的方法包括应用RF载波信号来发射第一RF脉冲序列。该载波信号具有一相位。在发射第一序列期间改变包络线，以便样本的第一区域基本饱和。发射第二RF脉冲序列以在第一区域内建立共振区域，并测量样本的属性。
在另一实施例中，在装置和样本之间有相对移动的NMR测量装置包括至少一个建立静态磁场的磁体、一第一线圈、一第二线圈和一脉冲发生器。脉冲发生器与第一线圈和第二线圈连接，并使其适合于应用第一线圈发射第一RF脉冲序列以产生一时变磁场。第一序列包括至少一再聚焦脉冲以产生至少一个来自样本的共振区域的回波。脉冲发生器还适合于在发射第一序列期间内应用第二线圈来至少短暂地改变静态磁场一次以使一区域饱和度比共振区域更大。
在进一步的实施例中，一种供在装置和样本之间有相对移动的NMR测量装置使用的方法包括应用反向恢复序列，该反向恢复序列该包括至少一个或多个绝热脉冲。
根据本发明的一种用于井孔的NMR测量装置的方法，在该装置和地层构造的样本之间有相对移动，该方法包括：
a)产生一个静态磁场；
b)发射第一RF脉冲序列；
c)在所说的第一序列的发射过程中改变一个参数以使地层构造的样本的第一区域基本饱和；和
d)发射第二RF脉冲序列，以便在第一区域内建立共振区域；和
e)测量样本的特性。
此外，根据本发明的一种用于井孔的NMR测量装置的方法，在该装置和地层构造的样本之间有相对移动，该方法包括：：
a)产生一个静态磁场；
b)发射第一RF脉冲序列，该第一序列包括至少一个再聚焦脉冲，该脉冲产生来自地层构造的样本的共振区域的至少一个回波；
c)在所述序列的发射期间改变所述静态电场至少一次，以便使比共振区域更大的地层构造的区域饱和；和
d)测量地层构造的样本的特性。
此外，根据本发明的一种用于井孔的NMR测量装置，在该装置和地层构造的样本之间有相对移动，包括：
a)至少一磁体；
b)至少一线圈；和
c)与线圈相连的脉冲发生器，该脉冲发生器适于：
i)通过所说的线圈发射第一RF脉冲序列；
ii)在所说的第一序列的发射过程中改变一个参数以使地层构造的样本的第一区域基本饱和，其中所说的参数是所说的序列的一条包络线或RF载波信号的相位；
iii)通过所说的线圈发射第二RF脉冲序列，在所说的第一区域内建立共振区域；和
iv)测量样本的特性。
附图说明
通过下面的详细描述，本发明的其它实施例将变得更清楚。
附图1为地下井的示意图。
附图2为沿着附图1中的2-2线的井的横截面图。
附图3为在NMR测井仪移动后的井的另一横截面图。
附图4和5为CPMG脉冲序列的波形。
附图6为驰豫时间T2的分布示例。
附图7为依据本发明的一个实施例的基于极化测量流程图。
附图8、9和10为依据本发明的不同的实施例的NMR测井仪的示意图。
附图11为沿着附图10中的11-11线的NMR测井仪的横截面图。
附图12为NMR脉冲序列的波形。
附图13、16、18和20所示为共振区域饱和的等值线图。
附图14、15、17和19和21所示为从NMR测井仪的周围的饱和区域接收来的对比信号幅值图。
附图22和23所示为应用有和没有交错自由发展时间的不同数目的脉冲时的共振区域饱和的等值线。
附图24和25所示为应用有和没有交错自由发展时间的不同数目的脉冲时的共振区域饱和的等值线。
具体实施方式
参考附图7，易于移动的NMR测量装置(例如NMR测井仪)可以应用依据本发明实现基于极化的T1测量的方法实例50。因为所测量的样本可能会移动，当样本、测量装置或两者都有移动时，可以应用这一测量方法。方法50包括使样本中的一个区域中的自旋饱和的步骤(块52)，该样本是将要测量其特性的样本；接着，经过预定的时间间隔(块54)，以使至少在该区域中的自旋发生部分极化；随后，方法50包括应用(块56)检测序列的步骤(例如，基于CPMG的序列)以产生来自样本的共振区域的自旋回波。如下文进一步描述，当NMR测量装置移动了时，应用一些技术来使饱和区域的边界和饱和密度最大以使共振区域基本保持在饱和区域内。应用这些的技术可以减少测量误差，并且(例如)如果应用在较小梯度的几何结构中时，不需要对NMR测量装置进行稳定的稳定装置。
例如，如下文进一步更详细地描述，测量方法50可以用于绘制地下地层构造的特性，并且也可以应用到其它的在样本和NMR测量装置之间产生相对移动的应用中(比如，其它的“内检外测型”的NMR应用)。在一些实施例中，NMR测量装置可能包括电磁场产生元件(比如，一个线圈、一个电磁铁和一个永磁体)以产生至少两个磁场：一个称为的磁场(未示)和一与磁场基本垂直的称为的磁场(未示)。参考附图8，例如，在一些实施例中，NMR测量装置可能是一NMR边钻井边记录(LWD)的测井仪60，例如，该测井仪60包括建立磁场的环行永磁体32和34和建立时变磁场的线圈39。在一些实施例中，磁场可能(当施加脉冲时)具有称为的射频(RF)载波分量。
磁场的载波频率一般以ω₀表示。发射磁场产生一具有一定径向厚度的共振区域，依据频率，该径向厚度是由在激励区域中的ω₀和ω₁的梯度决定的，这里为在磁场上的投影。在一些实施例中，如下文所述，磁场还可能(至少部分地)由梯度线圈40和42产生，以使磁场具有一个随着低频变化的分量。NMR测井仪60还可能包括处理电路，该处理电路包括比如连接到线圈(例如，线圈39、40和42)的脉冲发生器65，并适合于以下文将要描述的方式发射磁场和/或磁场
大体上，每个基于极化的NMR测量都包括三个标准块52、54和56(参见附图7)，以及可采用一次或多次测量来获得每个T1值。然而，可以应用检测序列(即，块52)来实现饱和(即，执行块56的功能)，因此，如果满足了两个条件：连续地重复测量(称为“成堆的(stacked)”试验)和信号检测序列68能够完全消除为进行下一步测量的磁化，就可以取消块52。如果应用这一技术，则放弃来自第一次测量中的结果，因为第一次测量是以一不正确的极化时间来执行的。作为一种变型，仅在应用检测序列之前，通过在共振区域中应用绝热快速变迁(passage)脉冲来绝热地进行激励。
三个基本块52、54和56还可能有其他变型。作为另外一个例子，顺序块54——块56——块52可以用来执行每次测量，这种变型从编程角度来看是较优的。当应用第二种变型时，放弃第一次测量。只要实现了块52、54和56的功能方法50的其他变型都属可能。
不管是通过显式的饱和序列或还是通过检测序列进行饱和，饱和的目的都是通过射频(RF)发射使一个较大的区域或空间饱和。如下文更详细的描述和模拟实验表明，根据特定的实施例，可以选择如下方式实现饱和：应用RF脉冲序列(比如CPMG检测序列)，该序列是适合于应用NMR测井仪60的移动来达到所需的饱和状态；不管NMR测井仪60有还是没有移动都缓慢改变序列的时间特性；不管NMR测井仪60有还是没有移动都随机地改变序列的特性；或者结合应用这些技术。
一种具有恒定参数的单一CPMG序列在自旋分布区域中产生强烈的饱和区域，该区域称为“孔穴”(holes)。由于孔穴都彼此完全分开，虽然实现了深远的孔穴加热，但其仅导致弱的饱和。此外，一旦破坏了在孔穴的位置处的磁化，继续施加序列不可能进一步增强饱和。如下文进一步描述，通过对分布在饱和空间中的这些孔穴进行脉冲“扫描”移动NMR测井仪60可以增强饱和密度。
可以修改CPMG检测序列，以将再聚焦脉冲数量增加到测量最初的回波串幅值所必须的典型的再聚焦脉冲数目(例如10)之上。如果在极化时间内的NMR测井仪60的移动总是与在检测序列过程中的NMR测井仪60的移动相联系，则这种方法效果很好。然而，不幸的是，如果NMR测井仪60在检测序列过程中是静止的而在极化时间内是运动的，则产生的饱和不令人满意。通过模拟(下文讨论)表明，如下文进一步描述，即使在测井仪没有移动的情况下，也可以通过缓慢改变序列随时间的特性以扩展饱和区域来避免这个问题。在本文中，例如术语“序列的特性”一般地指序列的包络线或RF载波频率的相位。作为改变包络线的可能的方法的例子，包络线可能包括脉冲120(参见附图12)，每个脉冲具有一定的持续时间(称为t_p)，并且脉冲120由称为t_e的时间间隔分隔开(从脉冲的中心到脉冲的中心)。这样，如下文进一步描述，可以改变持续时间t_p和/或时间间隔t_e(作为举例)以扩展饱和区域。
如下文进一步描述，不仅可以缓慢改变检测序列的特性(即，用于实现饱和目的的序列)，而且还可以以无关联的或随机的方式从一个脉冲到另一个脉冲改变检测序列的特性。随机极限值(实现完全随机)是发射不相干噪声。由于序列的相干的、非随机特性占主导作用，特性中的随机变化与缓慢变化形成对比，在特性的缓慢变化序列中饱和效果深远。结果，缓慢变化特性可能导致由连续脉冲递增地加热极度失谐孔穴。在很短的时间间隔内产生的饱和点彼此完全分离。然而，随机变化特性使序列的相邻脉冲对同一孔穴没有贡献，在很短的时间间隔中产生的饱和能够更均匀地被传播开。因此，随机变化的脉冲通常形成更一致的饱和密度。如下文所述(和通过模拟说明)，可以结合使用这两种技术以增强序列的性能。如下文进一步描述，如果出现的移动足够快，以致使能够在几个脉冲周期内扫描在将相邻的孔穴分开的距离上的孔穴，则序列的相干成分完全被破坏了，具有缓慢变化特性的序列可以与具有随机变化特性的序列类似地执行。
如下文所述，如果结合其它的变化(例如，载波频率相位的变化)，在CPMG序列中的再聚焦脉冲的翻转角度(flip angles)不需要大到足以产生失谐饱和。因此，通过缩短RF脉冲，可以降低饱和所需的能量。对于足够短的脉冲，加热孔穴的影响可以忽略。由于这种情况，可以压缩在脉冲之间的自由发展时间(free evolution period)，这样就可以在短得多的时间内实现饱和。在很短的脉冲极限内，这种技术引起发射不相干噪声(可以设计其结构以符合需要)。实际上，有限的脉冲上升和下落时间设定了脉冲持续时间的下限。如下文所述，在实现饱和所需的能量和时间和饱和带宽之间有一种折衷方案。
应用CPMG序列实现饱和
下面，详细讨论应用有和没有缓慢移动引起的的改变的CPMG序列实现饱和的例子。虽然在本说明书中作为例子特别指CPMG序列，但上文所述的加热孔穴可以通过各种多相脉冲序列来实现，该多相脉冲序列具有大量重复脉冲标准块的特点。
在CPMG序列过程中所重复的相干脉冲激励所选择的Δω＞＞ω₁的自旋，这里ω₁大致等于共振区域的径向厚度，并且Δω(频率空间中的距离)由下式具体定义：
Δω=γB→0-ωrf,]]>这里ω_rf是用于第一CPMG序列的磁场的RF频率。
随着Δω增加激励步幅变得越来越小，但激励从一个脉冲到一个脉冲都总和起来，在孔穴中激励幅值达到极其大的值。因为横向磁化与T2一致地衰减，所选择的自旋变得“饱和”。通过序列的周期确定这些孔穴的间隔(称为Δω_h)。不可忽视的脉冲持续和失谐效应产生某些偏差，因此，孔穴的间隔Δω_h大致由下式描述：Δωh=2πte,]]>这里t_e是一个再聚焦脉冲的开始到下一个再聚焦脉冲的开始之间的回波间隔。
结合驰豫，单一CPMG序列技术引起以一定的失谐频率加热孔穴。在被加热的孔穴之间是不可能测量的，因为测量区域的宽度Δω延伸达Δω_s≈2Δω₁以上，其对于持续时间为t_p的180°再聚焦脉冲变为Δω_s≈2π/t_p。因为，t_e总是比t_p大，Δω_s＞Δω_h，所以可能有几个孔穴被加热变成共振区域。为计算信号损失程度，磁场的几何结构、驰豫时间和检测带宽都必须加以考虑。
为说明孔穴的分布，附图13所示为二维等值线(contour plot)80(从模拟实验中得出的)，它是所计算的孔穴分布的等值线，该孔穴被加热成在水平轴线上线性变化ω₀和在垂直轴线上线性变化t_p的M_z＝1的纵向磁化。白色区域表示完整地保持磁化，而黑色区域表示对100％饱和的偏离或反相磁化。在Δω＝0处应用第一CPMG序列，所示为在紧接着这种CPMG序列的结尾处对失谐磁化M_z的影响。CPMG脉冲序列的参数为t_e＝500μs，t_p180＝125μs，k＝1000，这里k为再聚焦脉冲数。选择驰豫时间较长，但仅是回波串的持续时间的一小部分。在这种模拟实验中，应用完全矩形脉冲。然而，本发明的实施例可以应用基本为矩形的脉冲，并且也可以应用基本是非矩形的脉冲。在附图13中，没有对第一激励脉冲的效果进行仿真(模拟)。
附图14表明对于几个驰豫时间，在移动横坐标Δω的频率下第二次测量可以得到模拟合成的对比信号幅值82(即M_z/M_∞)，并且当平均Δω＝±0.75ω时，在ω₁·t_p＝π的第一次测量(如前文所述)中饱和减少该幅值。这就意味着在测量之间载波频率ω₀已经被移动了Δω。每个对比信号幅值82对应于不同的时间T1(例如，等于2*T2)。第二次测量的参数与第一次测量的参数相同，并且选择脉冲的翻转角度为180°。在附图中(以及在模拟实验中)，假设dω1dω0=0,]]>即，在共振区域附近磁场的变化可以忽略。对于轴对称的梯度几何结构，横坐标(Δω/ω₁)与在第一和第二次测量之间的在(共振区域的)径向上的差值成比例。当半径的差值比半径小得多时，上述假设ω₁为常数是一种有效的近似，一种证明了在图中选择的恒定的翻转角度是正确的事实。
正如从附图14所看到的，饱和区域延伸基本不超过2·Δω/ω₁，即两倍共振区域的厚度。因此，如果共振区域的径向移动小于1·Δω/ω₁，则下一测量步骤只是在完全饱和时开始。附图15所示的每个对比信号振幅84都与在第一序列中大量的再聚焦脉相联系。我们可以看到，在较小的Δω处的大部分饱和都发生在开始的10个回波内了。在这里以及在下面的例子中，选择T1＝2·T2＝100毫秒(msec)。
在第一CPMG序列期间测井仪移动可能导致在附近共振区域中增加损失。例如，附图16所示为在第一序列过程中测井仪的移动速率为-20ω₁/s的失谐磁化M_z产生的等值线86。横轴表示失谐频率Δω₁与第一CPMG序列的(脉冲振幅)ω₁的比率。等值线描述了在第一CPMG序列后剩下的相对纵向磁化。假设脉冲幅值恒定。脉冲参数和驰豫时间与上文的相同。纵轴表示有多少个再聚焦脉冲应用到具有载波ω_RF的第一CPMG序列上，该数目大致与这种序列的持续时间成比例。再聚焦脉冲数具有从图顶部的一个再聚焦脉冲(即，一个块跨越大约500μs)到图底部的100再聚焦脉冲(即，一个块跨越大约50ms)的范围。在这个例子中，在50ms内，NMR测井仪60移动距离为+1ω₁(其大致为半个外形宽度)。开始时，载波ω_RF与Δω＝0对应，最后，载波ω_RF与Δω＝+1·ω₁对应。如图所示，随着回波数的增加，移动NMR测井仪60将“扫描”在自旋分布区域上的孔穴，并由此增加了饱和密度。
附图17所示为当在整个外形(为模拟的目的采用矩形)宽度为±0.75ω₁上进行平均时所得的对比信号幅值(即，M_z/M_∞)。从顶部到底部，幅值88代表k＝1，11，21，31，41，51，61，71，81，91的结果。需注意的是，如附图14和17所示，随着回波数的增加损失加大，并且对于超过10个以上的回波，变得比在NMR测井仪60没有移动时的饱和效果强烈得多。现在饱和区域具有超过5·ω₁的宽度。在与自旋的驰豫时间可比拟的一段时间内损失增加，并且这种损失对于小的Δω甚至能够引起负信号。精确的分布取决于自旋集合体的驰豫时间和移动。对于越小的驰豫时间分布变得越窄。
前文假设在CPMG序列中的脉冲都是完全的矩形脉冲。然而，并不能达到这种理想的真正的“矩形”脉冲，而是存在有限的上升和下落时间。这就限制了包含在脉冲中的频谱宽度。在极度失谐状态中，被加热孔穴的宽度和加热它们的速率与在孔穴的位置处的脉冲的频率分量的幅值成比例。因此，在一些实施例中，极度失谐加热孔穴比在上述模拟中的更不显著。
对于在本申请中所讨论的脉冲，具有较宽的频率分布比较有利。因此，在一些实施例中，具有可能的最短上升和衰减时间常数的矩形脉冲比较可取。此外，通过变化脉冲包络线的形状以适合于脉冲的频率内容可以将饱和区域最优化。
一般地，在没有移动时，通过辐射具有变化参数和带宽脉冲的重复多脉冲序列可以产生深度的饱和。如果在应用该序列的同时缓慢改变脉冲序列参数，被加热孔穴的位置可以缓慢移过自旋分布区域，并增强饱和。被改变的脉冲序列参数包括：
·脉冲间隔t_e的变化，
·t_p的变化，
·比如通过脉冲幅值、磁场方向和载波频率ω_RF对的改变，·的变化，和
·脉冲相位的变化。
这些参数的综合变化和其它参数的变化都是可能的。可以通过使磁场和的实际变化(比如，变化磁铁和天线间隔或取向和/或射频功率)或样本和NMR测井仪60之间的相对移动来改变和在这种方式中，样本相对于NMR测井仪60之间的相对移动可以来自于样本的移动(比如流体流动或扩散)或来自于测井仪的移动。
改变的另一个方法是借助于电磁铁或梯度线圈来改变静态磁场。例如，回过来参考附图8，在一些实施例中，NMR测井仪60可能包括永磁体上部永磁体32和下部永磁体34，该永磁体32和34在NMR测井仪60的套筒28周围，并产生一径向、轴向对称的磁场在平行于NMR测井仪35的纵轴线的方向将磁体32和34极化以使其相互配合提供一低的梯度磁场例如，磁体32和34的北极可能彼此相对以提供一磁场该磁场的磁力线从NMR测井仪60的纵轴径向地延伸出来。在一些实施例中，磁力可穿透的件36围在套筒28的周围，并且位于上部磁体32和下部磁体34之间。这样排列的结果是磁力可穿透件36使磁场聚集以使磁场的梯度最小，因而，在我们所关心的区域中产生一更均匀的磁场NMR测井仪60可能包括或可能不包括套筒36。在第09,033,965号美国专利申请(题为“Nuclear MagneticResonance Apparatus and Method For Generating an AxisymmetricMagnetic Field Having Straight Contour Lines in the ResonanceRegion”，1998年3月3日申请)和美国专利US4,350,955(题为“Magnetic Resonance Apparatus”，1982年12月21日授权)中都可以找到对这些装置的更详细的描述，这两份专利(专利申请)在这里都以引用的方式结合在本申请中。
为改变磁场NMR测井仪35包括梯度线圈，比如线圈40和42，该线圈也围绕在套筒28的周围。对线圈40和42施加脉冲形式的DC电流(通过脉冲发生器，比如脉冲发生器65)以在磁场中产生一附加磁场分量如果在线圈40和42中的电流流动方向相反，则基本为径向。可以将线圈40和42设置在磁体32和34之间以便使两线圈40和42都给磁场提供一正的磁场分量，该磁场分量与在关注的区域中的磁场的方向可能基本一致或可能不一致，这取决于不同的实施例。在一些实施例中，线圈40和42可以由一对单匝或多匝的电流回路构成，该回路中的电流大小相等，循环方向相反。例如，该线圈40和42可以由鞍形线圈构成。
结合径向、轴向对称的磁场的设计方案应用线圈40和42的其它实施例都是可能的。例如，参考附图9，在另一NMR测井仪61中，用一环行永磁体62替换永磁体32和34，该永磁体62围绕在套筒36的周围，并且位于线圈40和42之间。磁体62产生磁场其磁力线与测井仪61的轴线平行、轴向地延伸。为使磁场基本上与磁场平行，在线圈40和42中的电流必须以相同的方向流动。作为一个例子，磁体62的顶部形成磁体62的北极，磁体62的底部形成磁体62的南极。
除了如上文所述的径向、轴向对称的的设计方案外的其它结构也是可能的。例如，梯度线圈可以应用具有两维(2-D)双极性的设计方案。2-D双极性的设计方案的例子在美国专利No.US 5,280,243(题为“System For Logging a Well During the Drilling Thereof”1994年1月18日授予Melvin Miller专利权)中可以找到。按照这种方式，应用2-D双极性设计的NMR测井仪68包括一环行磁铁72，该磁铁建立一如附图10和11中所示的双极性模式的磁场与其在NMR测井仪60和61中的对应的部件不同的是，RF线圈73和74不与NMR测井仪68的纵轴同轴，而是，排列RF线圈73和74以在磁场中产生双极性模式，这样就使在共振区域中磁场的等值线基本垂直于磁场的等值线。测井仪68包括梯度线圈76和77，每个线圈都包括一个或多个矩形环以产生梯度磁场，该梯度磁场与在关注区域中的由磁铁72所建立的磁场的方向一致。
这样，由于前文所描述的结构，自旋以大约的频率进动。如果两个矢量平行则产生最大的效应。因而，这种技术的结果就是改变了Δω而不改变ω_rf。这对于改变ω_rf有利，因为具有高的品质因数的天线的带宽限制了ω_rf变化的范围(不需重新调整天线，这至少在饱和序列过程中如果应用机械开关通过转换电容器来调整是不切实际的)。在一些实施例中，如果它必须以在整个饱和序列中的不同幅值来激励，则这种方法的缺点是需要相对较大的能量来驱动电磁铁(与在成像设备中的应用相比较)。有以下几种方式应用梯度线圈(或多个线圈)：
·在(磁场的)一个脉冲周期内在梯度线圈中产生基本恒定的电流以有效地改变这个脉冲的共振区域的半径。
·在(磁场的)一个脉冲周期内改变在梯度线圈中的电流以产生一“扫描”(sweep)脉冲而不改变rf脉冲的频率。依据实际的参数，扫描脉冲可以使特定的区域反相、激励特定的区域或使特定的区域饱和。可以在反相恢复序列中应用这种技术(而不是在饱和序列中)以使NMR测井仪周围较大的区域反相恢复。
·在(磁场的)脉冲之间激励梯度线圈以破坏可能保留的横向磁化。如果梯度脉冲持续时间(称为t_grad)足够短，那么在整个饱和区域上的α=ω→0gradienttgrad]]>的变化可以忽略，这类似于随机地改变磁场的脉冲的相位。
·在与施加磁场的每个脉冲的同时在梯度线圈中施加脉冲形式的电流。
·可以应用梯度线圈来产生如上文所述的随机或连续变化。
梯度线圈的其它应用都是可能的。
随机变化的CPMG序列
CPMG序列的脉冲串的特性也可以随机地改变。例如，可以随机地改变RF载波脉冲的相位以随机地产生比如0°、90°、180°和270°脉冲相位(至少在典型的核磁共振(NMR)波谱仪中可以产生这几种相位)。参见附图18(所示为不同的回波数的对比信号损失的等值线90)和附图19(所示为当在±0.75ω₁的空间厚度上进行平均时不同的回波数的对比信号损失的等值线92)，在所示的实例中的脉冲是随机地产生的，并且测井仪60和35不移动。除了这种时脉冲相位的随机化外，所有的自旋和脉冲参数都与上文所述的实施例相同。
我们可以看到，在自旋分布区域中饱和加热出了许多较宽的、完全分离的条纹区。被饱和区域的宽度比由受移动影响的CPMG序列产生的区域的宽度小，但饱和分布图比由CPMG序列产生的分布图平坦得多。这表示在共振区域的范围(应用相干特性)和可靠量的饱和分布(应用随机特性)之间的一种折衷方案。必须注意的是，对于T_1，2(这里100ms)＜＜t_m(这里50ms)的自旋(这里t_m是CPMG序列持续时间)，CPMG序列产生的分布具有更光滑的形状。在应用随机相位序列过程中产生移动能够稍稍增加其性能，但分布仍然保持光滑。
因为不是每个孔穴都是以相同的“速度”被加热，所以产生不完全饱和的条纹区。如我们所见，取决于不同的位置Δω，甚至某些孔穴被完全抑制，作为一个例子，在附图20中每第四个孔穴就看不见。这些未完全饱和的点的位置取决于再聚焦脉冲的持续时间：失谐，一持续时间为tp的脉冲使自旋围绕指向方向的“有效转动轴”转动α(Δω)=ω12+ω2tp]]>的角度。在α等于2π的倍数的地方出现未饱和点。因此，通过改变ω₁·t_p，可以使这些点饱和。
在附图20(所示为对于不同的回波数的对比信号损失的等值线94)和附图21(所示为当在±0.75ω₁的径向空间积厚度上进行平均时不同的回波数的对比信号损失的等值线96)中图解了缓慢增加脉冲长度(在附图12中以“t_p”表示)的实例的这种效果。在这个模拟实验中，脉冲长度从在第一次再聚焦脉冲中的125μs(180°脉冲)线性地增加到在第100次再聚焦脉冲中的250μs(360°脉冲)，而同时t_free(脉冲之间的距离，如在附图12中所描述的)保持不变。所有的其它参数与在前面的例子中都相同。所得的饱和分布比没有改变脉冲长度所得的分布更光滑并且稍稍宽些。
再次，一般地，对于特定移动范围可以通过改变序列的各种参数以使脉冲序列的饱和效果最优化，比如改变t_e(t_e大约与被加热孔穴的间隔成反比例)、t_p和脉冲相位等，以及在相干特性和随机特性之间的折衷方案。
前面的饱和序列的例子应用极度失谐孔穴加热效应来产生饱和。如上文所述，持续时间为t_p的脉冲总是使失谐的自旋转动比正常的翻转角度α(0)大的角度α(Δω)。因而，对于具有α(0)＝180°的再聚焦脉冲(即，180度脉冲)，对于失谐它总是保持α(Δω)＞180°。另一方面，如果α(Δω)＝(2n+1)·180°则发生最佳的激励，并且由此产生最佳激励失谐。然后，自旋转离纵轴的有效翻转角度为θ＝θ_max，此时对于一给定的Δω最大的有效翻转角度为θmax=α(Δω)2arctan(ω1Δω)]]>因而，应用180°脉冲来产生失谐饱和可能浪费能量。
附图22和23说明了饱和分布(在整个共振空间厚度上进行平均)与应用在序列中的再聚焦脉冲的α(0)的关系。如前文所述随机地改变相位。在附图22中，所示为自由发展时间t_free(即，在如附图12中所示的在再聚焦脉冲之间的时间间隔)被设定为375μs的对比信号损失图98，在附图23中，所示为t_free设定为0的对比信号损失图100。在两附图22和23中所示都是Δω函数的翻转角度为9°、20°、30°、45°、90°和180°的脉冲从1到100个的信号损失98和100。通过改变脉冲持续时间t_p来产生不同的翻转角度。如我们所见，对于不同的时间t_free信号损失分布几乎相同，因而，在随机相位变化的情况下，饱和模式主要由脉冲持续时间决定，而不是由自由发展持续时间决定。
对于一给定的硬件能够应用的最小脉冲持续时间由脉冲的上升时间常数(称为t_τ)决定。如果t_p＜3t_τ，则在其截止之前脉冲没有达到最大值ω₁，当t_p进一步减小时，它很快变得几乎无效。对于钻井记录NMR装置一个较好的估计是t_τ＝5...30μs。
当t_p减少时，饱和区域变得更宽。实际上关注的区域主要是符合|Δω|<(2πtp)2-ω12]]>的区域，即，在两个内部未饱和点内的α(Δω)＜2π的区域。随着Δω的增加最大的翻转角度θ_max减小。因此，饱和区域越宽需要在外部区域产生饱和的脉冲越多。如果饱和的时间常数为T_s，则只有T₁＞T_s的自旋能够完全饱和。因此，可以实现一种在饱和带宽和仍然能够实现饱和的、最低的T₁之间的折衷方案。这还说明，在某些实施例中，通过将t_free降低到用可行硬件(这里的硬件问题可能包括相位切换时间、脉冲上升和下落时间以及用较长的连续的rf脉冲使RF电路过载)所能够实现的可能的最低值，以保持序列尽可能地短是有利的。
附图24和25所示分别为具有(附图24)和没有(附图25)t_free的序列的损失102和104。所示的损失102和104针对不同的驰豫时间。对于t_free＝375μs，100个再聚焦脉冲的序列是40ms长，如果没有自由发展时间，该序列只有2.4ms长。对于标称的翻转角度α(0)＝35°时，两序列都能够使自由流体的驰豫时间(T1＞50ms)的自旋饱和，但没有自由发展时间的序列能够使具有比T₁低20倍的自旋饱和，这对于那些想解决约束流体内部的自旋分布问题的人是需要的。在这两情况中，产生饱和所需要的能量是单个180°再聚焦脉冲所需要的能量的100&CenterDot;35180&ap;20]]>倍，这对于井下的核磁共振(NMR)波谱仪来说并不存在严重问题，这种核磁共振波谱仪通常能够在t_w期间从存储在电容器中的能量中产生好几百个180°再聚焦脉冲串。
在一些实施例中，由具有包括自由发展时间的序列加热的分布比由连续辐射的序列加热的分布要稍微平滑些。这可能源于在自由发展时间内产生的附加相位差(在第二种情况中没有)，但这不是决定性的。此外，如果具轴向几何对称磁场的测井仪错位的距离为根据每个自旋在其位置的方位角，每个自旋都要在频率空间中移动不同的位移Δω=dω0/dγ&RightArrow;&CenterDot;Δγ&RightArrow;.]]>这就使实际的饱和分布更加显著地光滑。
在模拟实验中应用随机发生器选择四个脉冲相位。因此在不同的模拟实验中序列性能有稍微的变化。在某些实施例中，可以应用一预定相位的序列以使饱和性能最优。在某些实施例中，最佳参数的变化是没有周期性的。
总之，上面描述了对于在NMR共振区域的邻近区域中的自旋的预处理的示例性技术。即使在NMR装置(比如，NMR测井仪60或35)相对样本之间是移动的，这些技术都允许基于极化的T1测量，并且这些技术允许至少连同如美国专利申请(申请系列号为NO.09,033,965，如上所引用)中所描述的具有低的梯度场一起在不稳定地进行钻井的同时进行基于极化的测量。使测井仪在没有稳定器的情况下能够运行，这样就使其成为更“与钻机友好”，由此极大地增加边钻井边记录(LWD)测井仪的使用性。
虽然通过有限数目的实施例描述了本发明，但对于本领域熟练的技术人员来说，在本申请公开的内容的帮助下能够理解并由此作出各种改进和变型。附加权利要求所覆盖的所有的这种改进和变型都将落入本发明的精神和范围内。