专利名称:与样本间有相对移动的核磁共振装置用的核磁共振脉冲序列的制作方法
技术领域:
本发明一般涉及对在核磁共振(丽R)区域附近的自旋的预处理。
背景技术:
核磁共振(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 Formation Evalution While Drilling Including a Shortened or Truncated CPMG Sequence")中或Miller在美国专利No. 5, 280, 243 ( 1994年1月18 号授权,题为"System For Logging a Well During the Drilling Thereof")中所描述的测井4义之一。
同其它大多数井下地层构造测量方法相比,NMR测量方法具有两 个明显的特征。首先,来自地层的NMR信号是来自一较小的共振空间, 比如通常为一薄的共振空间20a (参见附图2),共振空间20a的径向 厚度与磁场^ (未示)的大小成比例。例如,取决于不同的共振区域 的形状,该空间可以在一个方向小到1毫米(mm)而在另一个方向延 伸长达几英尺。第二, N匿测量可以不是瞬时的。如下文进一步描述, 结合这两点事实使N匿测井仪易于进行移动测量,比如丽R测井仪6 围绕井孔3的周围移动。如下文所述,为了进行腿R测量,NMR测井仪6包括建立称为S。 的静态磁场(未示)的永磁体、发射垂直于磁场S。的时变磁场^的 射频(RF)线圏或天线和接收来自与醒R测量对应的地层构造的自旋 回波的RF线圈或天线。并且这两个线圏可以结合在一单个的发射/接 收天线中。
作为 一个例子,通过发射NMR检测序列以使原子核产生自旋回波, 醒R测井仪6测量地层构造10的氢原子核的自旋-自旋驰豫时间T2。 依次分析自旋回波得到时间T2的分布,从这一分布中可以得到地层构 造的特性。例如,附图4中所描述的Carr-Purcell-Meiboom-Gill (CPMG)序列15就是这样的一种醒R检测序列。应用序列15可以得 到时间T2的分布,并可以应用这一分布来测定并绘制地层构造10的 特性。
应用CPMG序列15来测量时间T2的技术包括下面的步骤。在第 一步中,NMR测井仪6以合适的时间间隔发射磁场&以应用90°激励 脉冲14a来使氢原子核的自旋(最初与磁场S。方向一致)旋转90° 。 虽然没有图示,实际上每个脉沖都是RF载波信号的包络线或脉冲串。 在自旋从磁场S。方向旋转90。后,首先自旋立即开始在与磁场S。垂直 的平面中一致地进动,接着逐渐地失去同步。在第二步中,在NMR脉 冲14a随后的一固定时间T中,丽R测井仪6以脉冲的形式施加磁场^ 并持续一较长时间段(与NMR脉冲14a比较),以应用NMR再聚焦脉 冲14b来旋转进动的自旋一个附加角度180° ,同时使其载波相位移动 ±90° 。 NMR脉冲14b使自旋重新同步,并使自旋在180°再聚焦脉沖 14b后发射一辅助自旋回波(见附图5),该自旋回波的峰值时间大约 等于T。以间隔te (大致为2 ■ T)重复步骤2 "k"次(这里"k"称为 回波数,例如,可以假定是从几百到多达几千中的任何一个值)。对 于第三步骤,在完成自旋回波序列后,在聚集另一序列自旋回波的下 一CPMG序列15开始之前要求经过一等待周期(通常称为等待时间), 以使自旋返回到沿磁场S。的平衡状态。观测每个自旋回波序列的衰减 并用其来导出T2的分布。
时间T2'是在应用90°激励脉冲14a后使自旋具有不再一致地进 动的特点的时间。这里,在90°激励脉冲14a结束时,所有的自旋都 指向与静态磁场夂垂直的一个共同的方向,对于完全均匀的磁场自旋以称为拉莫频率的共振频率进动。拉莫频率可以表述为5。-Y S。,这里 Y为旋磁比, 一个核常数。然而, 一般地磁场S。并不是均匀的,在激 励之后,由于静态磁场S。的非均匀性使自旋具有在T2'内产生的相移。 这种衰减是可逆的,但产生回波的再聚焦脉冲14a可以使这种衰减反 相。此外,还产生不可逆的相移(自旋-自旋驰豫),并可通过T2时 间常数来描述这种相移。这就导致依据T2时间常数在CPMG序列中的 连续回波幅值的衰减。 一般地,应用"内检外测型"NMR测量T2》T2' 的自旋。
如上文所述,可以应用T2时间的分布来测定地层构造的特性。例 如,参考附图6,地层可能包括容纳有约束流体的小孔和容纳有自由的、 可采出的流体的大孔。可以应用分离边界时间T2 (在附图6中被称为 TCUT_。FF)来将T2分布分为两部分 一部分包括比T,—。FF小的时间(表示 约束流体)和一部分包括比T,-附大的时间(表示自由的、可采出的流 体)。
一般地,通过观测由特定的CPMG序列15产生的自旋回波16的 衰减来计算时间T2。不幸的是,钻杆柱5 (参见附图1)可能要经过剧 烈的横向运动。然而,时间T2大致与另一被称为自旋-晶格松驰时间 Tl的时间常数成比例。时间Tl具有使自旋沿磁场S。方向重新返回到 平衡状态的特性,由此,考虑两T1和T2时间,可以认为每个自旋在 Tl恢复过程中以极紧密的螺旋运动形式朝平衡位置运动回来。幸运的 是,时间Tl和T2大致成比例。因此,可以从测量的时间Tl导出时间 T2的分布。事实上,建立约束流体截止点的最初的工作是应用Tl来完 成的。然后根据T2将这些结果表述出来并应用在商业上。参见W. E. Kenyon, J. J. Howard, A.
Sezginer, C. Straley, A. Matteson, K. Horkowtitz, and R. Ehrlich,等 人的 Pore-Size Distribution and 薩R in Microporous Cherty Sandstones,论文LL (在第30届年度记录研讨会(the 30th Annual Logging Symposium )上提交的论文,SWPLA, 1989年6月11-14日)。 基于极化的观、J量可以应用反向恢复序列(inversion recovery sequence)或饱和恢复序歹'J (saturation recovery sequence)。 在 饱和恢复序列情况下,例如应用几个将磁化强度减少为零的90°脉冲 4吏自旋系统饱和。然后在应用监测脉沖或脉冲序列(比如CPMG序列)之前使自旋系统恢复达一可变长度的时间。反向恢复技术要求在原子
核本身的自旋已经与静态磁场成一线之后应用180°脉冲使自旋反向。 随着时间的过去,自旋按照Tl朝平衡方向衰减,但由于180°脉冲在 检测器中并不产生信号,并没有进行测量。然而,在衰减完成之前, 监测脉冲或脉冲序列(比如CPMG序列)中断了衰减,该监测脉冲或脉 冲序列将自旋旋转到测量平面(即,在检测器中产生信号)。有意义 的信息就是在最初的90。"读出"脉冲瞬间之后的信号幅值。这种幅 值完全依赖于在初始的180°脉冲和90°脉冲之间的恢复时间。接下 来测定幅值,允许自旋系统完全驰豫到平衡状态,然后重复脉冲序列。
作为反向恢复序列在井下应用的例子描述在Kleinberg等人的美 国专利No. 5,023, 551 (题为"Nuclear Magnetic Resonance Pulse Sequences For Use With Borehole Logging Tools" , 1991年6月 11曰授权)中。然而描述在美国专利,551中的反向恢复序列并不能应 用绝热脉冲,由此导致仅在一个狭窄的区域中勘测。此外,在存在移 动和"内检外测"的情况下,使一个区域饱和比使其完全反相恢复要 容易。因此,比较可取的做法是使区域饱和。
回过来参考附图2, 一般地,应用基于极化的测量而不应用上文所 描述的基于衰减的测量来测量时间Tl。在基于极化的测量中,每个基 于极化的测量首先包括应用饱和序列来使在共振区域(比如,在附图2 中所描述的圆柱形共振空间20a)中的自旋到达饱和的步骤。接着,经 过一极化时间将共振空间20a极化到静态磁场S。。随后,应用检测序 列(比如CPMG序列)以产生来自地层构造10的自旋回波。然后,分 析开始的少数几个自旋回波的幅值以确定孔隙度分布O (Tl)的极化 加权整数0> ( twait)。因为,只需观测开始的少数几个回波就可以确定 信号的幅值,故可以在比基于衰减测量T2所需的时间更短的持续时间 内测量T1,因而,测量Tl时更不适合于移动NMR测井4义6。以不同的
等待时间连续地重复(在适当的饱和序列之后)检测序列几次以得到 孔隙度分布①(Tl )。
作为一个例子,应用基于极化的测量来测量在位于饱和空间20b 中的共振空间20a(参见附图2)中的氢原子核的时间Tl。在这里,NMR 测井仪6首先使在饱和空间20b中的自旋饱和。然而,极化时间应该 足够长以允许NMR测井仪6能够在井孔内作充分的移动。在这种情况下,NMR测井仪6的移动使共振空间20a移动,以使丽R测井仪6接收 来自移动了的共振空间20a'(参见附图3)中的自旋回波,该共振空间 20a'部分落到最初的饱和空间20b之外。结果,移动了的共振空间20a' 包括一没有饱和自旋的区域(一般称为"移入新自旋"效应)和一带 有饱和自旋的最初的饱和空间20b的部分区域。不幸的是,基于极化 的NMR技术不允许在极化时间中移入"新自旋",因为新自旋可能引 入测量误差。例如,测量可能错误地指示在地层构造中有比实际出现 的更多的约束流体量。
在申请号为PCT/US97/23975的PCT申请(题为"Method For Formation Evaluation While Dri 11 ing,, , 1997年12月29日申请) 中描述了一种使更大区域饱和的一种方法。这一申请指出,在测量开 始时,发射覆盖一个相对较大频率范围和/或特别宽的带宽的一个或多 个射频脉沖或应用具有能够使NMR测井仪周围的圆柱形空间饱和的扫 描频率的一个或多个脉冲。该申请进一步描述了当由于测井仪移动超 出饱和区域引起测量无效时应用加速度峰值来测量。该申请还进一步 描述在测井仪上装备支座以阻止测井仪移动超出饱和区域。
因此,人们一直渴望将在NMR测量装置和所勘测的样本之间的相 对移动引起的误差最小化。
发明内容
本发明公开了一种供在装置和样本之间有相对移动的醒R测量装 置使用的方法。装置、样本或这两者都可能移动。在本发明一个实施 例中,该方法包括发射笫一RF脉冲序列。第一序列具有包络线。在发 射第一序列期间改变包络线,以便样本的第一区域基本饱和。发射第 二 RF脉冲序列以在第一区域内建立共振区域,并测量样本的属性。
在另一实施例中,供在装置和样本之间有相对移动的NMR测量装 置使用的方法包括应用RF载波信号来发射第一 RF脉冲序列。该载波 信号具有一相位。在发射第一序列期间改变包络线,以便样本的第一 区域基本饱和。发射第二 RF脉冲序列以在第一区域内建立共振区域, 并测量样本的属性。
在另一实施例中,在装置和样本之间有相对移动的丽R测量装置 包括至少一个建立静态磁场的磁体、 一第一线图、 一第二线圈和一脉 沖发生器。脉冲发生器与第一线圏和第二线圏连接,并使其适合于应用第一线圈发射第一RF脉冲序列以产生一时变磁场。第一序列包括至 少一再聚焦脉冲以产生至少一个来自样本的共振区域的回波。脉冲发 生器还适合于在发射第一序列期间内应用笫二线圏来至少短暂地改变 静态磁场一次以使一区域饱和度比共振区域更大。
在进一步的实施例中, 一种供在装置和样本之间有相对移动的NMR 测量装置使用的方法包括应用反向恢复序列,该反向恢复序列该包括 至少一个或多个绝热脉冲。
根据本发明的一种用于井孔的丽R测量装置的方法,在该装置和 地层构造的样本之间有相对移动,该方法包括
a) 产生一个静态磁场;
b) 发射第一RF脉冲序列;
c) 在所说的第一序列的发射过程中改变一个参数以使地层构造的 样本的第一区域基本饱和;和
d) 发射第二RF脉冲序列,以便在第一区域内建立共振区域;和
e) 测量样本的特性。
此外,根据本发明的一种用于井孔的丽R测量装置的方法,在该 装置和地层构造的样本之间有相对移动,该方法包括
a) 产生一个静态磁场;
b) 发射第一RF脉冲序列,该第一序列包括至少一个再聚焦脉沖, 该脉冲产生来自地层构造的样本的共振区域的至少一个回波;
c) 在所述序列的发射期间改变所述静态电场至少 一 次,以便使比 共振区域更大的地层构造的区域饱和;和
d) 测量地层构造的样本的特性。
此外,根据本发明的一种用于井孔的NMR测量装置,在该装置和 地层构造的样本之间有相对移动,包括
a) 至少一磁体;
b) 至少一线圈;和
c) 与线圈相连的脉冲发生器,该脉冲发生器适于
i) 通过所说的线圈发射第一RF脉冲序列;
ii) 在所说的第 一序列的发射过程中改变一个参数以使地层构 造的样本的第一区域基本饱和,其中所说的参数是所说的序列的一条 包络线或RF载波信号的相位;iii) 通过所说的线圈发射第二 RF脉冲序列,在所说的笫一区域 内建立共振区域;和
iv) 测量样本的特性。
通过下面的详细描述,本发明的其它实施例将变得更清楚。 附图1为地下井的示意图。
附图2为沿着附图1中的2-2线的井的横截面图。 附图3为在醒R测井仪移动后的井的另一横截面图。 附图4和5为CPMG脉沖序列的波形。 附图6为驰豫时间T2的分布示例。
附图7为依据本发明的一个实施例的基于极化测量流程图。 附图8、 9和10为依据本发明的不同的实施例的NMR测井仪的示 意图。
附图11为沿着附图10中的11-11线的丽R测井仪的横截面图。 附图12为醒R脉冲序列的波形。
附图13、 16、 18和20所示为共振区域饱和的等值线图。
附图14、 15、 17和19和21所示为从NMR测井仪的周围的饱和区 域接收来的对比信号幅值图。
附图22和23所示为应用有和没有交错自由发展时间的不同数目 的脉冲时的共振区域饱和的等值线。
附图24和25所示为应用有和没有交错自由发展时间的不同数目 的脉冲时的共振区域饱和的等值线。
具体实施例方式
参考附图7,易于移动的NMR测量装置(例如NMR测井仪)可以应 用依据本发明实现基于极化的Tl测量的方法实例50。因为所测量的样 本可能会移动,当样本、测量装置或两者都有移动时,可以应用这一 测量方法。方法50包括使样本中的一个区域中的自旋饱和的步骤(块 52),该样本是将要测量其特性的样本;接着,经过预定的时间间隔 (块54),以使至少在该区域中的自旋发生部分极化;随后,方法50 包括应用(块56)检测序列的步骤(例如,基于CPMG的序列)以产生 来自样本的共振区域的自旋回波。如下文进一步描述,当NMR测量装 置移动了时,应用一些技术来使饱和区域的边界和饱和密度最大以使共振区域基本保持在饱和区域内。应用这些的技术可以减少测量误差,
并且(例如)如果应用在较小梯度的几何结构中时,不需要对NMR测 量装置进行稳定的稳定装置。
例如,如下文进一步更详细地描述,测量方法50可以用于绘制地 下地层构造的特性,并且也可以应用到其它的在样本和NMR测量装置 之间产生相对移动的应用中(比如,其它的"内检外测型"的NMR应 用)。在一些实施例中,固R测量装置可能包括电磁场产生元件(比如, 一个线圈、 一个电磁铁和一个永磁体)以产生至少两个磁场 一个称 为S。的磁场(未示)和一与磁场S。基本垂直的称为^的磁场(未示)。 参考附图8,例如,在一些实施例中,NMR测量装置可能是一醒R边钻 井边记录(LWD)的测井仪60,例如,该测井仪60包括建立磁场S。的 环行永磁体32和34和建立时变磁场&的线圈39。在一些实施例中, 磁场^可能(当施加脉冲时)具有称为5。的射频(RF)载波分量。
磁场^的载波频率一般以w。表示。发射磁场^产生一具有一定径 向厚度的共振区域,依据频率,该径向厚度是由在激励区域中的w。和 c^的梯度决定的,这里^为Y ^在磁场^。上的投影。在一些实施例 中,如下文所述,磁场S。还可能(至少部分地)由梯度线圏40和42 产生,以使磁场S。具有一个随着低频变化的分量。丽R测井仪"还可 能包括处理电路,该处理电路包括比如连接到线圏(例如,线圈39、 40和42 )的脉冲发生器65,并适合于以下文将要描述的方式发射磁场 5。和/或磁场51。
大体上,每个基于极化的NMR测量都包括三个标准块52、 54和56 (参见附图7),以及可采用一次或多次测量来获得每个T1值。然而, 可以应用检测序列(即,块52)来实现饱和(即,执行块56的功能), 因此,如果满足了两个条件连续地重复测量(称为"成堆的(stacked )" 试验)和信号检测序列68能够完全消除为进行下一步测量的磁化,就 可以取消块52。如果应用这一技术,则放弃来自第一次测量中的结果, 因为第一次测量是以一不正确的极化时间来执行的。作为一种变型, 仅在应用检测序列之前,通过在共振区域中应用绝热快速变迁 (passage)脉沖来绝热地进4亍激励。
三个基本块52、 54和56还可能有其他变型。作为另外一个例子, 顺序块54——块56——块52可以用来执行每次测量,这种变型从编程角度来看是较优的。当应用第二种变型时,放弃第一次测量。只要实
现了块52、 54和56的功能方法50的其他变型都属可能。
不管是通过显式的饱和序列或还是通过检测序列进行饱和,饱和 的目的都是通过射频(RF)发射使一个较大的区域或空间饱和。如下 文更详细的描述和模拟实验表明,根据特定的实施例,可以选择如下 方式实现饱和应用RF脉冲序列(比如CPMG检测序列),该序列是 适合于应用NMR测井仪60的移动来达到所需的饱和状态;不管NMR测 井仪60有还是没有移动都緩慢改变序列的时间特性;不管固R测井仪 60有还是没有移动都随机地改变序列的特性;或者结合应用这些技术。 一种具有恒定参数的单一CPMG序列在自旋分布区域中产生强烈的 饱和区域,该区域称为"孔穴"(holes)。由于孔穴都彼此完全分开, 虽然实现了深远的孔穴加热,但其仅导致弱的饱和。此外, 一旦破坏 了在孔穴的位置处的磁化,继续施加序列不可能进一步增强饱和。如 下文进一步描述,通过对分布在饱和空间中的这些孔穴进行脉冲"扫 描,,移动NMR测井4义60可以增强饱和密度。
可以修改CPMG检测序列,以将再聚焦脉冲数量增加到测量最初的 回波串幅值所必须的典型的再聚焦脉沖数目(例如10)之上。如果在 极化时间内的麵R测井仪60的移动总是与在检测序列过程中的NMR测 井仪60的移动相联系,则这种方法效果很好。然而,不幸的是,如果 薩R测井仪6 0在检测序列过程中是静止的而在极化时间内是运动的, 则产生的饱和不令人满意。通过模拟(下文讨论)表明,如下文进一 步描述,即使在测井仪没有移动的情况下,也可以通过緩慢改变序列 随时间的特性以扩展饱和区域来避免这个问题。在本文中,例如术语 "序列的特性,, 一般地指序列的包络线或RF载波频率的相位。作为改 变包络线的可能的方法的例子,包络线可能包括脉冲120(参见附图 12),每个脉沖具有一定的持续时间(称为tp),并且脉沖120由称 为te的时间间隔分隔开(从脉冲的中心到脉冲的中心)。这样,如下 文进一步描述,可以改变持续时间tp和/或时间间隔te (作为举例)以 扩展饱和区域。
如下文进一步描述,不仅可以緩慢改变检测序列的特性(即,用 于实现饱和目的的序列),而且还可以以无关联的或随机的方式从一 个脉沖到另一个脉沖改变检测序列的特,性。随机极限值(实现完全随机)是发射不相干噪声。由于序列的相干的、非随机特性占主导作用, 特性中的随机变化与緩慢变化形成对比,在特性的緩慢变化序列中饱 和效果深远。结果,緩慢变化特性可能导致由连续脉冲递增地加热极 度失谐孔穴。在很短的时间间隔内产生的饱和点彼此完全分离。然而, 随机变化特性使序列的相邻脉冲对同一孔穴没有贡献,在4艮短的时间 间隔中产生的饱和能够更均匀地被传播开。因此,随机变化的脉冲通 常形成更一致的饱和密度。如下文所述(和通过模拟说明),可以结 合使用这两种技术以增强序列的性能。如下文进一步描述,如果出现 的移动足够快,以致使能够在几个脉冲周期内扫描在将相邻的孔穴分 开的距离上的孔穴,则序列的相干成分完全被破坏了 ,具有緩慢变化 特性的序列可以与具有随机变化特性的序列类似地执行。
如下文所述,如果结合其它的变化(例如,载波频率相位的变化),
在CPMG序列中的再聚焦脉冲的翻转角度(flip angles)不需要大到 足以产生失谐饱和。因此,通过缩短RF脉冲,可以降低饱和所需的能 量。对于足够短的脉冲,加热孔穴的影响可以忽略。由于这种情况, 可以压缩在脉冲之间的自由发展时间(free evolution period ),这 样就可以在短得多的时间内实现饱和。在很短的脉沖极限内,这种技 术引起发射不相干噪声(可以设计其结构以符合需要)。实际上,有 限的脉冲上升和下落时间设定了脉冲持续时间的下限。如下文所述, 在实现饱和所需的能量和时间和饱和带宽之间有 一种折衷方案。 应用CPMG序列实现饱和
下面,详细讨论应用有和没有緩慢移动引起的5 。的改变的CPMG序 列实现饱和的例子。虽然在本说明书中作为例子特别指CPMG序列,但 上文所述的加热孔穴可以通过各种多相脉冲序列来实现,该多相脉沖 序列具有大量重复脉冲标准块的特点。
在CPMG序列过程中所重复的相干脉沖激励所选择的A 03>>031的
自旋,这里(Oi大致等于共振区域的径向厚度,并且A 0)(频率空间中 的距离)由下式具体定义
△ co = Y 5。-coIf,这里o)rf是用于第一 CPMG序列的磁场S。的RF频率。
随着A W增加激励步幅变得越来越小,但激励从一个脉冲到一个 脉冲都总和起来,在孔穴中激励幅值达到极其大的值。因为横向磁化与T2—致地衰减,所选择的自旋变得"饱和"。通过序列的周期确定 这些孔穴的间隔(称为AC0h)。不可忽视的脉冲持续和失谐效应产生
某些偏差,因此,孔穴的间隔厶on大致由下式描述A^r^,这里
ts是一个再聚焦脉冲的开始到下一个再聚焦脉冲的开始之间的回波间 隔。
结合驰豫,单一CPMG序列技术引起以一定的失谐频率加热孔穴。 在i^加热的孔穴之间是不可能测量的,因为测量区域的宽度A (O延伸 达A ws 2A o),以上,其对于持续时间为tp的180°再聚焦脉冲变为△ (Os 2^ 。因为,te总是比tp大,A(O,厶00h,所以可能有几个孔穴
被加热变成共振区域。为计算信号损失程度,磁场的几何结构、驰豫 时间和检测带宽都必须加以考虑。
为i兌明孔穴的分布,附图13所示为二维等值线(contour plot)
80 (从模拟实验中得出的),它是所计算的孔穴分布的等值线,该孔
穴被加热成在水平轴线上线性变化w。和在垂直轴线上线性变化tp的
Mz=l的纵向磁化。白色区域表示完整地保持磁化,而黑色区域表示对
100%饱和的偏离或反相磁化。在A w-0处应用第一CPMG序列,所示为
在紧接着这种CPMG序列的结尾处对失谐磁化Mz的影响。CPMG脉冲序
列的参数为te-500 m s, tpls。=125 ms, k=1000,这里k为再聚焦脉冲数。
选择驰豫时间较长,但仅是回波串的持续时间的一小部分。在这种模
拟实验中,应用完全矩形脉冲。然而,本发明的实施例可以应用基本
为矩形的脉冲,并且也可以应用基本是非矩形的脉冲。在附图13中,
没有对第一激励脉冲的效果进行仿真(模拟)。
附图14表明对于几个驰豫时间,在移动横坐标A co的频率下第二
次测量可以得到模拟合成的对比信号幅值82 (即Mz/M ),并且当平均
△ co=±0. 75co时,在c^ . t,;r的第一次测量(如前文所述)中饱和减
少该幅值。这就意味着在测量之间载波频率w。已经净皮移动了 △ w。每
个对比信号幅值82对应于不同的时间Tl (例如,等于242)。第二次测量的参数与笫一次测量的参数相同,并且选择脉冲的翻转角度为
180° 。在附图中(以及在模拟实验中),假设^1=0,即,在共振区
、
域附近磁场^的变化可以忽略。对于轴对称的梯度几何结构,横坐标 (△ w/a^)与在第一和第二次测量之间的在(共振区域的)径向上的 差值成比例。当半径的差值比半径小得多时,上述假设co,为常数是一 种有效的近似, 一种证明了在图中选择的恒定的翻转角度是正确的事 实。
正如从附图14所看到的,饱和区域延伸基本不超过2 .紐/ ,即
/欲
两倍共振区域的厚度。因此,如果共振区域的径向移动小于1 ■ ,
则下一测量步骤只是在完全饱和时开始。附图15所示的每个对比信号 振幅84都与在第一序列中大量的再聚焦脉相联系。我们可以看到,在
较小的A Q5处的大部分饱和都发生在开始的IO个回波内了。在这里以
及在下面的例子中,选择11=2 . T2-100亳秒(msec)。
在第一CPMG序列期间测井仪移动可能导致在附近共振区域中增加 损失。例如,附图16所示为在第一序列过程中测井仪的移动速率为-20 cd/s的失谐磁化Mz产生的等值线86。横轴表示失谐频率厶c^与第一 CPMG序列的(脉冲振幅)c^的比率。等值线描述了在第一 CPMG序列 后剩下的相对纵向磁化。假设脉冲幅值恒定。脉冲参数和驰豫时间与 上文的相同。纵轴表示有多少个再聚焦脉冲应用到具有载波o^的第一 CPMG序列上,该数目大致与这种序列的持续时间成比例。再聚焦脉冲 数具有从图顶部的一个再聚焦脉冲(即, 一个块跨越大约500 ms)到 图底部的100再聚焦脉沖(即, 一个块跨越大约50ms)的范围。在这 个例子中,在50ms内,腿R测井仪60移动距离为+lc^ (其大致为半 个外形宽度)。开始时,栽波COkf与△ co-0对应,最后,载波co^与A co=+l . co,对应。如图所示,随着回波数的增加,移动NMR测井仪60 将"扫描"在自旋分布区域上的孔穴,并由此增加了饱和密度。附图17所示为当在整个外形(为模拟的目的采用矩形)宽度为 ±0. 75w,上进行平均时所得的对比信号幅值(即,Mz/M )。从顶部到 底部,幅值88代表15:=1,11,21,31,41,51,61,71,81,91的结果。需注 意的是,如附图14和17所示,随着回波数的增加损失加大,并且对 于超过10个以上的回波,变得比在NMR测井仪60没有移动时的饱和 效果强烈得多。现在饱和区域具有超过5 (Oi的宽度。在与自旋的驰 豫时间可比拟的一段时间内损失增加,并且这种损失对于小的A (O甚 至能够引起负信号。精确的分布取决于自旋集合体的驰豫时间和移动。 对于越小的驰豫时间分布变得越窄。
前文假设在CPMG序列中的脉冲都是完全的矩形脉冲。然而,并不 能达到这种理想的真正的"矩形"脉沖,而是存在有限的上升和下落 时间。这就限制了包含在脉冲中的频谱宽度。在极度失谐状态中,被 加热孔穴的宽度和加热它们的速率与在孔穴的位置处的脉冲的频率分 量的幅值成比例。因此,在一些实施例中,极度失谐加热孔穴比在上
述模拟中的更不显著。
对于在本申请中所讨论的脉沖,具有较宽的频率分布比较有利。
因此,在一些实施例中,具有可能的最短上升和衰减时间常数的矩形 脉冲比较可取。此外,通过变化脉冲包络线的形状以适合于脉冲的频 率内容可以将饱和区域最优化。
一般地,在没有移动时,通过辐射具有变化参数和带宽脉冲的重 复多脉冲序列可以产生深度的饱和。如果在应用该序列的同时緩慢改 变脉冲序列参数,被加热孔穴的位置可以緩慢移过自旋分布区域,并 增强饱和。被改变的脉冲序列参数包括 脉冲间隔t。的变化,
tp的变化,
.比如通过脉冲幅值、磁场方向和栽波频率WKF对&的改变,
^。的变化,和
脉冲相位的变化。
这些参数的综合变化和其它参数的变化都是可能的。可以通过使 磁场力。和A的实际变化(比如,变化磁铁和天线间隔或取向和/或射 频功率)或样本和NMR测井仪60之间的相对移动来改变^。和^。在 这种方式中,样本相对于NMR测井仪60之间的相对移动可以来自于样本的移动(比如流体流动或扩散)或来自于测井仪的移动。
改变^。的另一个方法是借助于电磁铁或梯度线圈来改变静态磁
场。例如,回过来参考附图8,在一些实施例中,NMR测井仪60可能 包括永磁体上部永磁体32和下部永磁体34,该永磁体32和34在NMR 测井仪60的套筒28周围,并产生一径向、轴向对称的磁场&。在平 行于NMR测井仪35的纵轴线的方向将磁体32和34极化以使其相互配 合提供一低的梯度磁场S。。例如,磁体32和34的北极可能彼此相对 以提供一磁场S。,该磁场S。的磁力线从NMR测井仪60的纵轴径向地延 伸出来。在一些实施例中,磁力可穿透的件36围在套筒28的周围, 并且位于上部磁体32和下部磁体34之间。这样排列的结果是磁力可 穿透件36使磁场S。聚集以使磁场S。的梯度最小,因而,在我们所关心 的区域中产生一更均匀的磁场S。。固R测井仪60可能包括或可能不包 括套筒36。在第09, 033, 965号美国专利申请(题为"Nuclear Magnet ic Resonance Apparatus and Method For Generating an Axisymmetric Magnetic Field Having Straight Contour Lines in the Resonance Region" , 1998年3月3日申请)和美国专利US4, 350, 955 (题为 "Magnetic Resonance Apparatus" , 1982年12月21日授权)中都 可以找到对这些装置的更详细的描述,这两份专利(专利申请)在这 里都以引用的方式结合在本申请中。
为改变磁场仏,NMR测井仪35包括梯度线圏,比如线圈40和42, 该线圈也围绕在套筒28的周围。对线圈40和42施加脉冲形式的DC 电流(通过脉冲发生器,比如脉冲发生器65)以在磁场So中产生一附 加磁场分量&。如果在线圈40和42中的电流流动方向相反,则&基 本为径向。可以将线圏40和42设置在磁体32和34之间以便4吏两线 圈40和42都给磁场S。提供一正的磁场分量,该磁场分量与在关注的 区域中的磁场力。的方向可能基本一致或可能不一致,这取决于不同的 实施例。在一些实施例中,线圈40和42可以由一对单压或多匣的电 流回路构成,该回路中的电流大小相等,循环方向相反。例如,该线 圈40和42可以由鞍形线圈构成。
结合径向、轴向对称的磁场》。的设计方案应用线围40和42的其 它实施例都是可能的。例如,参考附图9,在另一 NMR测井仪61中, 用一环行永磁体62替换永磁体32和34,该永磁体62围绕在套筒36的周围,并且位于线圏40和42之间。磁体62产生磁场&,其磁力线 与测井仪61的轴线平行、轴向地延伸。为使磁场&基本上与磁场S。 平行,在线圈40和42中的电流必须以相同的方向流动。作为一个例 子,磁体62的顶部形成磁体62的北极,磁体62的底部形成磁体62 的南极。
除了如上文所述的径向、轴向对称的S。的设计方案外的其它结构 也是可能的。例如,梯度线圈可以应用具有两维(2-D)双极性S。的设 计方案。2-D双极性S。的设计方案的例子在美国专利No.US5, 280,243 (题为 "System For Logging a Well During the Drilling Thereof" 1994年1月18日授予MelvinMiller专利权)中可以找到。按照这种 方式,应用2-D双极性S。设计的NMR测井仪68包括一环行磁铁72, 该磁铁建立一如附图10和11中所示的双极性模式的磁场S。。与其在 NMR测井仪60和61中的对应的部件不同的是,RF线圈73和74不与 醒R测井仪68的纵轴同轴,而是,排列RF线圈73和74以在磁场A 中产生双极性模式,这样就使在共振区域中磁场^的等值线基本垂直 于磁场S。的等值线。测井仪68包括梯度线團76和77,每个线圏都包 括一个或多个矩形环以产生梯度磁场,该梯度磁场与在关注区域中的 由磁铁72所建立的磁场S。的方向一致。
这样,由于前文所描述的结构,自旋以大约^。 + ^。^d^的频率进 动。如果两个矢量平行则产生最大的效应。因而,这种技术的结果就
是改变了 △ CO而不改变0^f。这对于改变COrf有利,因为具有高的品质 因数的天线的带宽限制了 Wrf变化的范围(不需重新调整天线,这至少 在饱和序列过程中如果应用机械开关通过转换电容器来调整是不切实 际的)。在一些实施例中,如果它必须以在整个饱和序列中的不同幅 值来激励,则这种方法的缺点是需要相对较大的能量来驱动电磁铁(与 在成像设备中的应用相比较)。有以下几种方式应用梯度线團(或多 个线圏)
在(磁场^的) 一个脉冲周期内在梯度线圏中产生基本恒定的 电流以有效地改变这个脉冲的共振区域的半径。
-在(磁场^的) 一个脉沖周期内改变在梯度线围中的电流以产 生一 "扫描"(sweep)脉冲而不改变rf脉沖的频率。依据实际的参 数,扫描脉冲可以使特定的区域反相、激励特定的区域或使特定的区域饱和。可以在反相恢复序列中应用这种技术(而不是在饱和序列中)
以使丽R测井仪周围较大的区域反相恢复。
在(磁场^的)脉冲之间激励梯度线圈以破坏可能保留的横向 磁化。如果梯度脉沖持续时间(称为t^d)足够短,那么在整个饱和区 域上的a-6f^"'t^d的变化可以忽略,这类似于随机地改变磁场^的 脉冲的相位。
在与施加磁场^的每个脉冲的同时在梯度线圈中施加脉冲形式 的电流。
可以应用梯度线圈来产生如上文所述的随机或连续变化。 梯度线圈的其它应用都是可能的。 随机变化的CPMG序列
CPMG序列的脉冲串的特性也可以随才;L地改变。例如,可以随机地 改变RF载波脉冲的相位以随机地产生比如0° 、 90° 、 180°和270° 脉冲相位(至少在典型的核磁共振(醒R)波谱仪中可以产生这几种相 位)。参见附图18(所示为不同的回波数的对比信号损失的等值线90) 和附图19 (所示为当在±0.75^的空间厚度上进行平均时不同的回波 数的对比信号损失的等值线92),在所示的实例中的脉沖是随机地产 生的,并且测井仪60和35不移动,除了这种时脉沖相位的随机化外, 所有的自旋和脉冲参数都与上文所述的实施例相同。
我们可以看到,在自旋分布区域中饱和加热出了许多较宽的、完 全分离的条紋区。被饱和区域的宽度比由受移动影响的CPMG序列产生 的区域的宽度小,但饱和分布图比由CPMG序列产生的分布图平坦得多。 这表示在共振区域的范围(应用相干特性)和可靠量的饱和分布(应 用随机特性)之间的一种折衷方案。必须注意的是,对于Ti,2 (这里 100ms) (这里50ms)的自旋(这里L是CPMG序列持续时间), CPMG序列产生的分布具有更光滑的形状。在应用随机相位序列过程中 产生移动能够稍稍增加其性能,但分布仍然保持光滑。
因为不是每个孔穴都是以相同的"速度"被加热,所以产生不完 全饱和的条紋区。如我们所见,取决于不同的位置Aco,甚至某些孔 穴被完全抑制,作为一个例子,在附图20中每第四个孔穴就看不见。 这些未完全饱和的点的位置取决于再聚焦脉冲的持续时间失谐,一持续时间为tp的脉沖使自旋围绕指向3j"方向的"有效转动轴"转 动a(Ao)="〖+ 0、的角度。在a等于2;r的倍数的地方出现未饱和点。 因此,通过改变(d tp,可以^使这些点饱和。
在附图20(所示为对于不同的回波数的对比信号损失的等值线94 ) 和附图21 (所示为当在±0. 7501的径向空间积厚度上进行平均时不同 的回波数的对比信号损失的等值线96)中图解了緩慢增加脉沖长度(在 附图12中以"tp"表示)的实例的这种效果。在这个模拟实验中,脉 沖长度从在第一次再聚焦脉冲中的125^ (180。脉冲)线性地增加到
在第100次再聚焦脉冲中的250卢(360°脉冲),而同时the (脉冲
之间的距离,如在附图12中所描述的)保持不变。所有的其它参数与 在前面的例子中都相同。所得的饱和分布比没有改变脉冲长度所得的 分布更光滑并且稍稍宽些。
再次, 一般地,对于特定移动范围可以通过改变序列的各种参数 以使脉冲序列的饱和效果最优化,比如改变L (L大约与被加热孔穴 的间隔成反比例)、tp和脉冲相位等,以及在相干特性和随机特性之 间的折衷方案。
前面的饱和序列的例子应用极度失谐孔穴加热效应来产生饱和。
如上文所述,持续时间为tp的脉冲总是使失谐的自旋转动比正常的翻
转角度a (0)大的角度a ( △ co )。因而,对于具有a (0) =180°的
再聚焦脉冲(即,180度脉冲),对于失谐它总是保持a (Aco) >180
° 。另一方面,如果a (厶00) = (211+1) '180°则发生最佳的激励,
并且由此产生最佳激励失谐。然后,自旋转离纵轴的有效翻转角度为
6 = 6_,此时对于一给定的△ co最大的有效翻转角度为6B =a ( △ w )
2arctan(鱼)。因而,应用180。脉冲来产生失谐饱和可能浪费能量。 附图22和23说明了饱和分布(在整个共振空间厚度上进行平均) 与应用在序列中的再聚焦脉冲的a (0)的关系。如前文所述随机地改 变相位。在附图22中,所示为自由发展时间tf^(即,在如附图12 中所示的在再聚焦脉冲之间的时间间隔)被设定为375^的对比信号损 失图98,在附图23中,所示为tr^设定为0的对比信号损失图100。在两附图22和23中所示都是a co函数的翻转角度为9。 、 20° 、 30 ° 、 45。 、 90°和180°的脉冲从1到100个的信号损失98和100。 通过改变脉冲持续时间tp来产生不同的翻转角度。如我们所见,对于 不同的时间t^信号损失分布几乎相同,因而,在随机相位变化的情 况下,饱和模式主要由脉冲持续时间决定,而不是由自由发展持续时 间决定。
对于一给定的硬件能够应用的最小脉冲持续时间由脉冲的上升时 间常数(称为U)决定。如果tp〈3tr,则在其截止之前脉冲没有达到 最大值co"当tp进一步减小时,它很快变得几乎无效。对于钻井记录 NMR装置一个较好的估计是tr-5…30戶。
当tp减少时,饱和区域变得更宽。实际上关注的区域主要是符合
<
2tt
(——)-02的区域,即,在两个内部未饱和点内的"(△ wX2;r的
区域。随着a (0的增加最大的翻转角度6^减小。因此,饱和区域越
宽需要在外部区域产生饱和的脉沖越多。如果饱和的时间常数为Ts,
则只有T々L的自旋能够完全饱和。因此,可以实现一种在饱和带宽和
仍然能够实现饱和的、最低的L之间的折衷方案。这还说明,在某些
实施例中,通过将tf^降低到用可行硬件(这里的硬件问题可能包括相位切换时间、脉冲上升和下落时间以及用较长的连续的rf脉冲使RF 电路过载)所能够实现的可能的最低值,以保持序列尽可能地短是有 利的。
附图24和25所示分别为具有(附图24)和没有(附图25) tfr^ 的序列的损失102和104。所示的损失102和104针对不同的驰豫时间。 对于tf^=375^, 100个再聚焦脉冲的序列是40ms长,如果没有自由 发展时间,该序列只有2.4ms长。对于标称的翻转角度a ( 0)-35°时, 两序列都能够使自由流体的驰豫时间(Tl〉50ms)的自旋饱和,但没有 自由发展时间的序列能够使具有比L低20倍的自旋饱和,这对于那些 想解决约束流体内部的自旋分布问题的人是需要的。在这两情况中, 产生饱和所需要的能量是单个180°再聚焦脉冲所需要的能量的
100 il《20倍,这对于井下的核磁共振(NMR)波谱仪来说并不存在
180
严重问题,这种核磁共振波谱仪通常能够在U期间从存储在电容器中 的能量中产生好几百个180°再聚焦脉冲串。在一些实施例中,由具有包括自由发展时间的序列加热的分布比 由连续辐射的序列加热的分布要稍微平滑些。这可能源于在自由发展 时间内产生的附加相位差(在第二种情况中没有),但这不是决定性 的。此外,如果具轴向几何对称磁场的测井仪错位的距离为^,根据 每个自旋在其位置的方位角,每个自旋都要在频率空间中移动不同的 位移A = "0^^ A卩。这就使实际的饱和分布更加显著地光滑。
在模拟实验中应用随机发生器选择四个脉冲相位。因此在不同的 模拟实验中序列性能有稍微的变化。在某些实施例中,可以应用一预 定相位的序列以使饱和性能最优。在某些实施例中,最佳参数的变化
是没有周期性的。
总之,上面描述了对于在NMR共振区域的邻近区域中的自旋的预 处理的示例性技术。即使在NMR装置(比如,NMR测井仪60或35 )相 对样本之间是移动的,这些技术都允许基于极化的Tl测量,并且这些 技术允许至少连同如美国专利申请(申请系列号为NO. 09, 033, 965,如 上所引用)中所描述的具有低的梯度场一起在不稳定地进行钻井的同 时进行基于极化的测量。使测井仪在没有稳定器的情况下能够运行, 这样就使其成为更"与钻机友好",由此极大地增加边钻井边记录(LWD ) 测井仪的使用性。
虽然通过有限数目的实施例描述了本发明,但对于本领域熟练的 技术人员来说,在本申请公开的内容的帮助下能够理解并由此作出各 种改进和变型。附加权利要求所覆盖的所有的这种改进和变型都将落 入本发明的精神和范围内。
权利要求
1.一种用于井孔的NMR测量装置,在该装置和地层构造的样本之间有相对移动,包括a)至少一磁体;b)至少一线圈;和c)与所述线圈相连的脉冲发生器,该脉冲发生器适于i)发射第一RF脉冲序列,以便产生一个时变磁场,所述第一序列包括至少一个再聚焦脉冲,以便产生来自样本的一个共振区域的至少一个回波;ii)建立一个静态磁场,以便建立所述共振区域;以及iii)在所述第一序列的发射期间瞬间改变所述静态电场至少一次,以便使比所述共振区域更大的区域饱和。
2. 如权利要求1所述的装置,其中改变所述静态磁场的步骤进 一步包括发射至少一个其它脉沖。
3. 如权利要求2所述的装置,进一步包括当所述RF脉冲序列未被发射时同步发射至少一个其它脉冲。
4. 如权利要求2所述的装置,进一步包括当至少一个RF脉冲发射时同步发射至少一个其它脉冲。
5. 如权利要求1所述的装置,其中饱和步骤包括发射附加RF 脉沖,以便使比所述共振区域更大的区域饱和。
6. 如权利要求l-5所述的装置,其中由于NMR装置的移动产生所述相对移动。
7. 如权利要求1-5所述的装置,其中由于样本的移动产生所 述相对移动。
8. 如权利要求1 - 5所述的装置,进一步包括产生轴向对称的 静态磁场的步骤,该轴向对称的静态磁场在共振区域中具有等值 线,该等值线在与丽R装置的纵轴基本一致的方向上基本是直线。
9. 如权利要求1 - 5所述的装置,进一步包括在地层构造中钻 井孔的同时测量样本的特性。
全文摘要
本发明公开了一种进行核磁共振测量的方法和装置。在测量过程中NMR测量装置、所测量的样本或这两者都可能移动。依据包络线来调制RF载波信号的包络线以产生一第一RF脉冲序列。在发射第一序列期间可以改变包络线、RF信号的相位和/或静态磁场以使样本的第一区域基本饱和。发射一第二RF脉冲序列以在第一区域内建立一共振区域,并测量样本的特性。
文档编号E21B49/00GK101295012SQ20081000584
公开日2008年10月29日 申请日期1999年12月6日 优先权日1998年12月4日
发明者K·加尼桑, P·斯佩尔 申请人:施卢默格控股有限公司