CN109154189A - 具有减少的运动伪影的t2反演 - Google Patents
具有减少的运动伪影的t2反演 Download PDFInfo
- Publication number
- CN109154189A CN109154189A CN201780029808.4A CN201780029808A CN109154189A CN 109154189 A CN109154189 A CN 109154189A CN 201780029808 A CN201780029808 A CN 201780029808A CN 109154189 A CN109154189 A CN 109154189A
- Authority
- CN
- China
- Prior art keywords
- nmr
- echo
- motion artifacts
- measurement data
- correction
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/448—Relaxometry, i.e. quantification of relaxation times or spin density
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Geophysics (AREA)
- Remote Sensing (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geochemistry & Mineralogy (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种用于处理核磁共振(NMR)测量数据的方法,包括:利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据具有受运动伪影影响的回波串,其中所述运动伪影与由于所述NMR工具的运动而在感兴趣的体积中变化的磁场大小相关。所述方法还包括:利用所述处理器通过使用校正反演方法来减小所述运动伪影对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法具有乘法项,所述乘法项具有包括至少一个局部最大值的项和可选衰减项。
Description
相关申请的交叉参考
本申请要求2016年4月4日提交的美国申请号15/089893的权益,该临时申请全部内容以引用的方式并入本文中。
背景技术
土地地层或简称地层可用于各种目的,诸如油气生产、地热生产以及二氧化碳封存。为了最佳地利用地层,通常使用通过穿透地层的钻孔输送的井下工具来表征。
一种类型的井下工具是核磁共振(NMR)工具,其对地层执行NMR测量以确定各种性质,例如像孔隙度。在称为随钻测井的一种应用中,NMR工具联接到钻柱。在钻柱旋转的同时NMR工具执行NMR测量,导致钻头也联接到钻柱以钻出钻孔。然而,钻井过程可能导致钻柱在钻孔中侧向移动,从而连续地改变从NMR工具到被表征的地层的距离。也可能由于在不钻井的情况下工具旋转而发生NMR工具的侧向运动。不幸的是,连续变化的距离可能在所获得的NMR数据中引起运动伪影,导致数据准确性的降低。
发明内容
公开的是一种用于处理核磁共振(NMR)测量数据的方法。所述方法包括:利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据具有受运动伪影影响的回波串,其中运动伪影与由于NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及利用处理器通过使用校正反演方法来减小运动伪影对NMR测量数据的影响,所述校正反演方法对运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法具有乘法项,所述乘法项具有包括至少一个局部最大值的项和可选衰减项。
还公开的是一种用于对土地地层执行核磁共振(NMR)测量的方法。所述方法包括:通过穿透土地地层的钻孔输送NMR工具;利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据具有受运动伪影影响的回波串,其中运动伪影与由于NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及利用处理器通过使用校正反演方法来减小运动伪影对NMR测量数据的影响,所述校正反演方法对运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法具有乘法项,所述乘法项具有包括至少一个局部最大值的项和可选衰减项。
还公开的是一种用于处理核磁共振(NMR)测量数据的系统。所述系统包括具有计算机可执行指令的存储器以及被配置来执行计算机可执行指令的处理器。所述计算机可执行指令包括:接收从NMR工具获得的NMR测量数据,所述NMR测量数据具有受运动伪影影响的回波串,其中运动伪影与由于NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及通过使用校正反演方法来减小运动伪影对NMR测量数据的影响,所述校正反演方法对运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法具有乘法项,所述乘法项具有包括至少一个局部最大值的项和可选衰减项。
附图说明
以下描述不应在任何方面视为限制。参照附图,相同元件用相同数字编号:
图1示出设置在穿透土地的钻孔中的核磁共振(NMR)工具的示例性实施方案的剖视图;
图2描绘理想化的NMR回波串和理想化的波列的各方面以用于说明NMR运动伪影;
图3A-图3C,统称为图3,描绘用于在存在运动伪影的情况下证明现有技术的联合反演的NMR回波串和波列的各方面;
图4描绘由现有技术的联合反演产生的T2分布的各方面;
图5A-图5C,统称为图5,描绘用于证明具有运动校正的联合反演的NMR回波串和波列的各方面;
图6描绘由具有运动校正的联合反演产生的T2分布的各方面;
图7描绘用于证明现有技术的联合反演的具有60ms的短等待时间的NMR长回波串和波列的各方面;
图8描绘由具有长等待时间的NMR长回波串以及具有60ms的短等待时间的波列的现有技术联合反演产生的T2分布的各方面;
图9描绘用于证明具有运动校正的联合反演的具有长等待时间的NMR长回波串和具有60ms的短等待时间的波列的各方面;
图10描绘由具有长等待时间的NMR长回波串以及具有60ms的短等待时间的波列的具有运动校正的联合反演产生的T2分布的各方面;
图11描绘用于证明具有运动校正的联合反演的具有长等待时间的NMR长回波串和具有不同等待时间(一秒和60ms)的两种波列类型的各方面;
图12描绘由具有长等待时间的NMR长回波串以及具有不同等待时间(一秒和60ms)的两种波列类型的具有运动校正的联合反演产生的T2分布的各方面;
图13预设了数值例子的表格结果;
图14是用于使用NMR测量数据估算土地地层的性质的方法的流程图;
图15描绘具有晚期非指数伪影的回波串的各方面;
图16描绘示出晚T2峰值的反演结果的各方面;
图17描绘使用现有技术算法的具有运动校正的反演结果的各方面;
图18示出具有变化振幅的模拟运动伪影的实例;
图19示出具有指数衰减的余弦函数的实例;
图20示出仅在余弦振幅设置为零的情况下指数衰减的实例;
图21示出具有指数衰减和周期变化的余弦函数;
图22描绘具有使用本文公开的增强公式拟合和校正的运动伪影的回波串的各方面;
图23描绘没有错误的晚T2峰值的校正T2分布的各方面;
图24描绘使用增强公式来拟合和校正非指数回波串的另一实例的各方面;并且
图25是用于使用NMR测量数据估算土地地层性质的另一种方法的流程图。
具体实施方式
参考附图,本文通过举例而非限制的方式呈现所公开的设备和方法的一个或多个实施方案的详细描述。
公开了用于处理由核磁共振(NMR)工具执行的测量的方法和设备,所述NMR工具可能在钻孔中经受运动。替代地或与工具运动组合,NMR工具可以具有非轴对称的磁场,使得当工具旋转时,磁场大小在地层中的固定位置处变化。运动可能导致NMR测量不准确地量化地层的性质。由于运动(或非轴对称磁场运动)而引入NMR测量中的不准确性被称为运动伪影。本文公开的处理技术识别运动伪影并将其从NMR测量数据移除以提供更准确地量化地层性质的校正的NMR测量数据。
下面讨论用于实现本文教导的设备。图1示出设置在穿透土地3(包括土地地层4)的钻孔2中的NMR工具10的示例性实施方案的剖视图。NMR工具10被配置来对地层4执行NMR测量。NMR测量包括从地层的原子核(诸如氢原子核)生成NMR信号回波。可以记录长回波串。回波串的衰减由所谓的T2弛豫引起,也称为横向或自旋-自旋弛豫。NMR测量产生横向弛豫时间T2,其是对应于地层4材料的特征或性质的指数衰减时间常数。横向弛豫涉及在NMR测量期间在处理静态磁场时地层4材料中质子的相位相干性的损失。对于地层流体,没有单一的T2值,但是例如在几毫秒到几秒之间的任何地方都有很宽的值分布。T2值的定量分布是NMR工具10的主要输出。相对于钻孔中的深度绘制的T2分布序列可以称为NMR T2分布对数。NMR工具10还可以输出与使地层中的核极化相关联的纵向弛豫时间常数(T1)。
T2(也称为T2)衰减可以通过指数函数之和(多指数近似)来近似,从而导致T2分布。获得此T2分布的过程通常称为T2反演、回声拟合或映射。根据T2分布,可以确定地层中的总孔隙度、部分孔隙度、孔隙尺寸和流体类型-特别感兴趣的性质。长T2分量通常被称为自由流体(FF)分量或可移动的总体积(BVM);中T2分量通常被称为结合水(BW)分量或不可减少的总体积(BVI);并且短T2分量通常被称为粘结合水(CBW)分量。
用于获得T2回波衰减的基本核磁共振方法是接近(或实现完全)平衡极化的长等待时间(TW长,极化时间),随后跟着生成的几百到几千个NMR回波的序列,例如,脉冲回波序列,诸如众所周知的Carr Purcell Meiboom Gill(CPMG)脉冲回波序列。平衡极化对于从回波串的起始振幅获得总孔隙度是有用的。除了具有长TW的长回波串之外,具有短TW短的回波串(称为波列或突发)用于更准确地确定T2分布中的短T2分量。通常但不是必须的,波列比具有长TW长的长回波串更短(即,它们具有更小数量的回波)。
在一个或多个实施方案中,获取并平均多个NMR回波串(具有长TW)。优选地TW足够长以完全极化所有NMR分量。获取并平均多个NMR波列(具有短TW)。通常,波列仅具有少量回波(以节省时间、内存和功率)并且使用短TW(在现有技术中:为了节省时间并因此增加测量的信噪比(SNR);在本公开中以减少运动伪影)。波列TW应足够长以使T2分量完全极化,所述T2分量稍后在反演之后从波列T2分布中提取。平均的波列的数量大于具有长TW的平均回波串的数量,以便更好地确定短T2分量。可以根据诸如下面讨论的那些之一的反演方法来处理(即,反演)这些获取的数据。
NMR工具10中的部件包括使地层材料磁化的静磁场源13以及精确地传输提供振荡磁场的射频能量的定时爆发的天线14。在这些脉冲之间的时间段内,天线从在传输的RF频率下与由静磁场源产生的静磁场共振的那些氢质子接收衰减回波信号。在围绕NMR工具10的环形体积(被称为感兴趣的体积9)中执行NMR测量。因为在质子共振频率与静磁场的强度之间存在线性关系,可以调节传输的射频能量的频率以匹配感兴趣体积中的静磁场。可以理解,NMR工具10可以包括如NMR领域中已知的各种部件和配置。由于NMR工具在本领域中是已知的,不再进一步详细讨论这些工具的部件和配置的具体细节。
NMR工具10由载体5输送通过钻孔2,所述载体可以是钻管,诸如钻柱6。钻头7设置在钻柱6的远端。钻机8被配置来进行钻井操作,诸如旋转钻柱6并因此旋转钻头7以钻出钻孔2。另外,钻机8被配置来将钻井泥浆(即,钻井液)泵送通过钻柱6,以便润滑钻头7并从钻孔2冲洗钻屑。井下电子设备11被配置来操作NMR工具10,处理在井下获得的测量数据,和/或充当遥测的接口以便在井下部件与设置在土地3的表面处的计算机处理系统12之间传递数据或命令。遥测的非限制性实施方案包括用于实时通信的脉冲泥浆和有线钻杆。系统操作和数据处理操作可以由井下电子设备11、计算机处理系统12或其组合执行。在替代实施方案中,载体5可以是铠装电缆,其还可以提供与表面处理系统12的通信。振动传感器14被配置来感测钻柱6的振动并向计算机处理系统12提供振动信号。转速传感器15被配置来感测转速,并因此感测钻柱6的转速变化,并且向计算机处理系统12提供转速信号。地质转向系统16被配置来根据选定的几何形状或迹线使钻柱6进地质转向。钻机控制器17被配置来控制钻机8的操作,例如像,使用从NMR工具10导出的信息来控制地质转向系统16。
接下来,讨论组合长回波串和波列的T2反演。存在可以用于实现这点的若干已知反演方法。参见例如A METHOD FOR INVERTING NMR DATA SETS WITH DIFFERENT SIGNALTO NOISE RATIOS,K.J.Dunn、D.J.Bergman、G.A.LaTorraca、S.M.Stonard和M.B.Crowe;1998年5月26日至29日的SPWLA第39届年度测井研讨会中的单独反演(SI)(也称为拼接技术)和联合反演(JI)(也称为复合数据处理)的变体。该论文在本公开中称为REF1。
在联合反演(JI)技术中,测量数据的多指数近似方程可以表示为:
其中EETi是具有长等待时间TWET的长回波串的时间ti处的第i个回波振幅,并且ETLj是具有短等待时间TWTL的波列的时间tj处的第j个回波振幅。φk是T2分布的受欢迎的T2分量,即,与选择的固定T2k(或T2k仓-T2分类到其中的选定间隔)相关联的指数函数的振幅,其中K从1到T2仓的选择数。φk在反演过程期间被优化以实现与测量的NMR数据的最佳拟合。i的范围从1到长回波串的回波数,而j的范围从1到波列回波的数。理想地,TWET应足够长(例如,>5*T1),其中T1是地层的最长的T1分量,以便完全极化所有NMR分量,在这种情况下,在EETi的公式中的项是1并且可以忽略。参数R需要在JI例程中进行优化。它是地层的T1/T2的测量(并且此外,如稍后将看到的,可以是运动伪影检测器)。换句话说,R不是直接计算为T1/T2,而是通过JI例程优化来拟合回波串数据。回波串和波列可能具有不同的回波间时间TE。长回波串和平均波列需要根据它们的平均数进行加权,这相当于它们反向测量误差的平方(参见REF1)。
在单独的反演(SI)中,平均回波串被反演(即,多指数拟合),给出回波串T2分布。平均波列被反演(与回波串分开,因此名称为SI),给出波列T2分布。产生最终T2分布的原理是在回波串T2分布中用波列T2分布的短T2分量来替换短T2分量。REF1描述了细节,包括用于提高精度的小修改,例如通过Chen和Georgi 1997的方法。
接下来,讨论运动伪影。在NMR回波的采样或接收期间,NMR工具的移动可能在长回波串的衰减中并且小程度地在波列中引起运动伪影。主运动伪影是T2分布中FF分量的减少。然后,主要在BW分量的增加中发现从FF分量中丢失的内容。单独观察长回波串,不可能确定BW分量(或其中的一部分)是真的BW还是运动伪影。
当使用现有技术的反演方法(诸如REF1)来使与具有短等待时间的波列组合的具有长等待时间的长回波串反演时,拟合的T2衰减不完美地拟合两种类型的回波串,或者在两种回波串类型经历相同的运动的情况下R=T1/T2拟合得不切实际地高。如本文公开的,该不拟合用于检测和校正NMR数据中的运动伪影。所公开的反演方法可以找到运动伪影并且移除或至少减少它们,从而导致更准确的NMR数据和从校正的NMR数据导出的更准确的地层性质值。
接下来,参考图2讨论NMR回波串中的运动伪影的表现。如图2所示的这种运动伪影仅在地层中存在自由流体(FF)时才会发生,此时FF的一些部分孔隙度通过运动转换成其中BW区域中有明显的T2的部分孔隙度。图2示出了这些想法。在图2中,横轴是时间轴[s],并且纵轴是NMR振幅。测量的NMR回波衰减可能看起来像迹线FFBW_ET(除了缺乏噪声)。仅从这个回波串,不能确定该回波串开始时的快速衰减是否由于约束水(BW)含量为20%并且剩余部分自由流体(FF)含量为80%,或者BW是否不是真实的并且实际上是运动伪影。在后一种情况下,实际FF将为100%(跟踪FF_ET)。可以看出,波列的外观允许将运动伪影与真实的结合水区分开。观察图中左下方的两个波列迹线,其等待时间(即,极化时间)较短,使得FF分量仅以小程度被极化。迹线FFBW_TL显示了如果NMR样品由80%的自由流体(FF)和20%的BW组成则波列将看起来怎样。FF内容几乎被短的TW完全抑制,而BW内容仅保留了一点衰减。迹线FFmot TL显示了如果NMR样品由100%FF组成并且在迹线FFBW_ET开始时的快速衰减是运动伪影则波列将看起来怎样。在这种情况下,如图所示,波列迹线将是迹线FFBW_ET的早期时间的缩小版本。因此,通过将波列的形状与长回波串的对应开始部分的形状进行比较,可以确定该部分是由于BW还是由于运动伪影。在一个或多个实施方案中,BW和运动伪影可以同时存在。在这种情况下,通过眼睛识别运动伪影可能是具有挑战性的。相反,可以使用如本说明书中进一步描述的算法。
幸运地,在长TW回波串中存在的运动是否在波列中存在并不那么重要。这是因为无论如何,波列显示出非常少的运动伪影,因为这里考虑的运动伪影只能在存在FF NMR信号的情况下产生;然而,在具有短TW的波列中,FF NMR信号受到很大抑制。
如果长TW回波串中的运动伪影看起来类似于快速衰减的结合水分量,则此概念是适用的。如果伪影由轴对称NMR工具的偏心旋转引起或者当其磁场不是完美轴对称时由中心旋转引起,则通常是这种情况。如果运动在回波串中稍后开始(诸如通过突然的冲击),则运动伪影不是问题,因为它不太可能被误认为是错误的BW分量。
接下来,讨论了与运动伪影校正(JIMC)的联合反演。上面呈现的两个可拟合方程(1)用附加的乘法项修改,拟合了运动效果。可能的项是其中Amot和Tmot分别是运动伪影的严重程度(振幅)和特征瞬态时间。T2反演的完整组的可拟合方程现在是:
其中R=T1/T2需要联合反演例程中的优化,并且长回波串和平均波列需要根据它们的平均数进行加权,这相当于它们反向测量误差的平方(参见REF1)。在第一个方程中,假设等待时间TWET足够长以完全极化NMR核。如果不是这种情况,则需要在方程(1)中添加适当的恢复项。新的参数是Amot和Tmot。未知的是:φk、R=T1/T2、Amot和Tmot。将这些可拟合参数约束到真实值≥0是有用的。例如,这可以通过在拟合之前用方程中的这些参数的平方根的平方来代替,或者通过一些其他方式来完成。在一个或多个实施方案中,R被约束为大于一并且Amot被约束在零与一之间。
在上述拟合方程(2)中,假设在长TW回波串和短TW波列中存在相同的运动。实际上,不一定是这种情况但应该是次要的结果,因为波列显示非常少的FF信号,并因此在波列中没有多少FF可以通过运动转换成BW。
接下来,讨论使用长回波串和两种波列类型的具有运动校正的联合反演。方程组(2)可以扩展用于具有完全极化的一个回波串和具有不同等待时间TWTL1和TWTL2的两个波列。然后T2反演的完整组的可拟合方程变成:
可以理解,该方程组可以扩展到任何数量的回波串和波列。具有中等待时间TW(即,大于短TW且小于长TW)或具有两个不同TW的波列对于拟合参数R的实际确定特别有用,注意到R的实际确定是用于获得运动伪影的真实拟合的前提条件。
应注意,方程(2)和(3)的不同回波串不需要具有相同的回波间时间TE。这些方程可以用于具有不同数量的回波NE和不同的回波间时间TE的回波串一起使用。
接下来,呈现了应用与运动伪影校正的联合反演的实例。这些实例基于模拟的NMR回波衰减与模拟的运动伪影。对于本部分,参数是:FF是大于100ms的T2分量,BW是3.3毫秒(ms)与100ms之间的T2分量,并且CBW是小于3.3ms的T2分量。
呈现了自由流体(FF)实例,其中T2的一个驰豫分量=1秒,并且R=T1/T2=1.5。长回波串完全极化,而波列使用60ms的等待时间并且加权96次,噪声降低倍。将现有技术的联合反演(JI)与本发明的具有运动校正的联合反演(JIMC)进行比较。
图3描绘用于在存在运动伪影的情况下证明现有技术的联合反演的NMR回波串和波列的各方面。图3A示出1000个回波(ET)的噪声模拟回波衰减以及波列(左下角)。横轴是时间轴,以秒为单位。纵轴是回波振幅轴,对于没有伪影和噪声且完全极化的回波串的起点归一化为100%。虚线迹线指示无伪影的单指数衰减,特征衰减时间为1秒-预期的无噪声且无运动伪影的回波衰减。噪声迹线与虚线迹线之间的偏差由噪声和运动伪影引起。现有技术的联合反演获得长回波串(迹线FET)和波列(迹线FTL)中的拟合。图3B示出了图3A,但是时间轴被放大(即,放大)以仅显示0秒至0.05秒的时间。图3C也是图3A的放大视图,但示出了波列TL与其拟合FTL的细节。在图3B中可以看出,在ET开始时,反演的拟合(FET)系统性地太低(这导致总孔隙度太低)。类似地,在图3C中,拟合FTL显示与波列TL的系统偏差(即,错配)(即,错误的斜率)。(运动校正拟合应找到真正的运动无伪影回声衰减ETorg。)通过现有技术的联合反演获得的T2分布在图4中给出,其中水平时间仓轴以秒为单位。图4中的Φk的分布提供了对某些性质的确定:
总孔隙度:95.1%(真实为100%);
FF:88.0%(真实为100%);
BW:7.1%(真实为0%);
CBW:0.0%(真实为0%);并且
R=T1/T2:3.2(真实为1.5)。
不切实际地高拟合的R以及图3B和图3C中的错配都是波列TL与回波串ET关于不包括运动伪影的数据处理的不相容性的指示。虽然错配(与要拟合的数据的系统偏差)显而易见或者可以通过适当的算法确定,太高拟合的R可以通过与特定地层类型的可能R进行比较来识别。(具有运动校正的JI,即JIMC,应找到更接近真实R的R)。可以通过在不同地层类型中测试NMR工具或通过分析来确定地层类型的可能或参考R,使得可以将拟合R与可能R进行比较。因此,通过知道正被钻的地层类型的拟合R是否与可能的R相差超过阈值,则可以提供拟合的R与NMR数据不相容的指示。
图5类似于图3并且示出了使用JIMC的噪声长回波串迹线(ET)和对该迹线的拟合(FET)。图5B和图5C是图5A的放大版本(方式与图3B和图3C是图3A的放大版本的方式相同)。FET迹线正确地再现了长回波串的开始处的总孔隙度(图5B),并且图5C的波列的拟合优于图3C的拟合。使用JIMC得到的T2分布在图6中示出。如果绘制了从图6中的上述T2分布产生的多指数衰减,则获得图5A和图5B中的虚线迹线(ETorg)正下方的迹线(ETcor)。这是运动校正拟合,距离无伪影点迹线ETorg不远。与真实T2分量的差异非常小。几乎完全移除了运动伪影。在T2分布中具有两个峰值的FF分量可能看起来不规则。这是因为方程(2)的JIMC尚未包括适当的正则化。在一个或多个实施方案中,将使用这种正则化。图6中的Φk的分布提供了对某些性质的确定:
总孔隙度:99.9%(真实为100%);
FF:99.7%(真实为100%);
BW:0.2%(真实为0%);
CBW:0.0%(真实为0%);
R=T1/T2:1.9(真实为1.5);
Amot:0.13;并且
Tmot:0.033s。
使用JIMC对这些性质的估计比使用上段中的现有技术JI的估算更准确。
使用JI和JIMC的另一个实例关于表征泥质沙子呈现。这是具有三个驰豫分量的模拟实例:T2的60%FF=1s,T2的30%BW=15ms,T2的10%CBW=1.5ms,并且对于所有分量,R=T1/T2=1.5。长回波串完全极化,而波列使用60ms的等待时间并且加权96次,噪声降低倍。实例进一步使用一个或两个波列比较JIMC的输出。本实例中使用的数字等同于自由流体实例中使用的数字。图7和图8涉及使用具有长回波串和等待时间TW=60ms的一个波列的现有技术JI。图8中的Φk的分布提供了对某些性质的确定:
总孔隙度:102.0%(真实为100%);
FF:51.6%(真实为60%);
BW:40.6%(真实为30%);
CBW:9.9%(真实为10%);并且
R=T1/T2:2.7(真实为1.5)。
图9和图10涉及使用具有长回波串和等待时间TW=60ms的一个波列的JIMC。图10中的Φk的分布提供了对某些性质的确定:
总孔隙度:102.0%(真实为100%);
FF:51.8%(真实为60%);
BW:40.2%(真实为30%);
CBW:9.9%(真实为10%);
R=T1/T2:2.6(真实为1.5);
Amot:0.000004;并且
Tmot:0.131s。
注意,在该实例中,使用JIMC的准确度与使用现有技术JI的准确度相当(Amot拟合为几乎零)。然而,找到良好的拟合其中过度R=2.6。由于良好的拟合,很明显,对于该示例,在长回波串和波列中不包含用于运动校正的足够信息。
因为在先前的实例中,在长回波链和波列中没有足够的信息,添加第二波列,其中等待时间TW=1秒,其与第一波列的等待时间(TW=60ms)不同。图11和图12涉及使用具有长回波串的JIMC,具有等待时间TW=60ms的一个波列,以及具有等待时间TW=1s的一个波列。图12中的Φk的分布提供了对某些性质的确定:
总孔隙度:101.0%(真实为100%);
FF:60.4%(真实为60%);
BW:29.7%(真实为30%);
CBW:10.9%(真实为10%);
R=T1/T2:1.4(真实为1.5);
Amot:0.16;并且
Tmot:0.041s。
在具有两个波列的该实例中,JIMC找到正确的R和运动伪影并且正确地移除伪影。此外,参见图11,可以看出,从图12的T2分布导出并且在图11中绘制的所有拟合(FET、FTL1、FTL2)都是优异的。图11中的FET迹线是拟合方程(3)的结果,并且使回波串与其伪影拟合。图11中的ETcor迹线是校正的回波串,使用拟合的T2分布和方程(3)的第一个方程(EET)中的R,而不使用运动项校正的回波串忠实地再现原始回波串而没有运动伪影或噪声。
图13呈现了以上呈现的实例的表格结果。最后一行显示,具有两个波列的JIMC大致获得了正确的T2分量,即使在具有以下性质的泥质沙子的更困难的情况下:总孔隙度=100%;CBW(0至3.3ms)=10%;BW(3.3至100ms)=30%;并且FF(0.1至4秒)=60%。
在前面的描述中描述的运动伪影校正的实施方案是所公开的具有运动伪影校正的联合反演(JIMC)的变体,其是现有技术的联合反演(JI)的修改。然而,应理解,其他实施方案像具有运动伪影校正的单独反演(SIMC)的变体完全在本专利申请的范围内,它们是现有技术的单独反演(SI)的修改。
现在讨论当SI用于T2数据的反演(SIMC)时校正运动伪影的过程的各方面。在SIMC中,平均回波串被反演,给出回波串T2分布。(假设回波串中存在运动伪影,这导致短T2分量的振幅太高而长T2分量的振幅太低。)现在,通过在回波串T2分布中用波列T2分布的短T2分量来替换短T2分量来产生组合的T2分布(目前为止与SI相同)。这种组合减少了替换的短T2分量中的运动伪影,但是长T2分量的振幅仍然太低,并因此总孔隙度也太低。现在为了减少分布的长T2分量中的运动伪影,对波列的短T2分量求和,并且从回波串的总和短T2分量中减去该总和,从而产生差。随后,该差分布在组合的T2分布的长T2分量上。这导致获得最终T2分布,其具有运动伪影减少的短T2和长T2分量以及更正确的总孔隙度。通过使用REF1中描述的Chen和Georgi 1997的方法可以改善结果。
在现有技术SI中,波列的TW常常设置得非常短,仅仅是为了使CBW而不是BW极化。由于运动效果通常在BW区域中表现出来,如上所述,SIMC在减少运动伪影方面可能不是有效的。相反,除了具有正好足够长以使CBW极化的TW的CBW波列之外,还需要具有足够长以使BW分量极化的TW的BW波列。然后如在现有技术SI中使用CBW波列,而BW波列用于运动伪影校正。然后处理的顺序是:通过在回波串T2分布中用CBW波列T2分布的CBW-T2分量来替换CBW-T2分量并且用BW波列T2分布的BW分量替换BW-T2分量来产生组合的T2分布。如果在回波串的被替换的BW-T2分量中存在运动伪影,则现在在组合的T2分布中减少了该运动伪影,但是长T2分量仍具有太低的振幅,并因此总孔隙度也太低。现在为了减少分布的长T2分量中的运动伪影:对BW波列的BW-T2分量求和,并且从回波串的总和BW-T2分量中减去该总和,从而产生差。随后,该差分布在组合的T2分布的长T2分量上。这导致获得最终T2分布,其具有运动伪影减少的短T2和长T2分量以及更正确的总孔隙度。通过使用REF1中描述的Chen和Georgi 1997的方法可以再次改善结果。可以理解,一些运动伪影可能不会影响BW T2分量,但可能影响CBW T2分量或两者。可以使用与以上描述的用于校正BW区域中的运动伪影的技术类似的技术来校正受CBW区域中的运动伪影影响的T2分量。
图14是用于处理核磁共振(NMR)测量数据的方法140的流程图。框141要求利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据受到运动伪影的影响,并且包括利用长极化时间TWET获得的第一回波串和利用比TWET短的短极化时间TWTL获得的第二回波串。方框142要求以下中的至少一个:(i)利用处理器使用第一回波串和第二回波串减小运动伪影对NMR测量数据的影响,以及(ii)利用处理器使用第一回波串和第二回波串识别运动伪影,其中运动伪影与由于NMR工具的运动而在感兴趣的体积中变化的磁场大小相关。运动可能由于NMR工具的径向移动、轴向振动以及磁场的非轴对称旋转中的至少一个。NMR工具可以由诸如钻管的载体输送通过穿透土地地层的钻孔。
关于方法140中的减小,减小可以包括使用校正反演方法,所述方法对运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布。校正反演可以包括使用以下乘法项:
其中Amot表示运动伪影的振幅,Tmot表示运动伪影的瞬态时间常数,并且t表示时间。减小可以包括当第二回波串是单个回波串时使用方程(2),并且当第二回波串包括两个回波串时使用方程(3)。
关于方法140中的识别,识别可以包括确定使用不对运动伪影进行建模的非校正反演方法获得的第一回波串和第二回波串的T2分布的多指数近似是否提供长回波串与短回波串之间的不相容性的指示。不相容性的指示可以由用户确定,该用户检查多指数近似并且向识别处理器输入存在不相容性的指示。在一个实例中,基于感兴趣的土地地层的细节,当R不合理地高时,可能存在不相容性的指示。方程(1)可以用于多指数近似。识别还可以包括执行提供不相容性指示的算法。在算法的一个实例中,算法可以包括方程(1)。如果两个方程(1)的联合拟合较差或拟合的R过高,则指示运动伪影。在算法的另一个实例中,算法可以包括方程(5)。
对于方程(5),在获取EETi时,极化时间被选择为足够长以基本上极化FF,并且在获取ETLj时,极化时间被选择为足够长以基本上极化BW而不是FF。随后,当取决于BW的那些φEk的总和基本上大于(例如,大于多于10%)取决于BW的那些φTk的总和时,检测运动伪影。取决于BW的φEk和φTk是所有φEk和φTk的子集。
方法140还可以包括通过减小运动伪影的影响来提供校正的T2分布并且随后使用校正的T2分布来估算土地地层的性质。
运动伪影检测以及运动伪影校正可以优选地实时地在井下或井上执行,或者在后处理NMR数据时在井上执行。
注意,在某些情况下,使用上面的项(4)作为校正反演的乘法项可能无法校正NMR回波串的后半部分中的运动伪影。图3A中的ETi,1曲线示出了受常规运动影响的模拟NMR回波串。在0.2秒和0.4秒处,一些非指数伪影(局部最大值)是可见的。具有运动校正的拟合曲线(FET)遵循ETi,1曲线的平均值,但不遵循这些伪影。在某些情况下,这些晚期伪影可以通过拟合曲线进行平均,并且仅对反演结果(即,T2分布)具有可忽略的影响。然而,在该领域中获取的数据表明,这种“平均”可能并不总是发生,并且后期伪影可以基本上改变所得到的T2分布。因此,本文公开了增强的公式以校正晚期伪影。
图15描绘了回波串150的各方面,其中在回波串的末端圈出单个晚期非指数伪影(高原)151。在这种情况下,伪影151延伸到回波串的末端。在其他情况下,它可能不会。在大多数情况下,预期晚期伪影151将出现在回波串的后半部分中。
反演过程试图将指数曲线拟合到伪影的此非指数数据。结果,图15中环绕的晚T2峰值出现在T2分布中,因为只有具有长衰减时间的分量可以在数学上“解释”T2曲线在末端处不下降。因此,T2反演的多指数拟合的质量低,使得反演结果可疑。
图17描绘在将校正项(4)应用为校正反演过程中的乘法项之后产生的反演的各方面。在应用校正项(4)之后,注意到晚期T2峰值(环绕的)仍然存在。校正项(4)是在回波串的开始处校正伪影的乘法指数衰减,因此,不能正确拟合图17中所示的非指数晚期伪影。
使用NMR模拟器软件,分析回波串的形状以确定规则的偏心运动对其的影响。该回波串信号形状由余弦函数近似,具有另外的指数衰减,如下所示:
方程(6)具有四个参数,Acos:余弦振幅;Tcos:余弦周期时间;Aexp:指数衰减振幅;以及Texp:指数衰减时间常数。指数部分考虑了去相位,这通常与运动一起发生。由于工具旋转,余弦项解释了周期运动。
图19示出具有以下参数的回波信号曲线的实例:Acos=0.1,Tcos=200ms,Aexp=0.3并且Texp=57ms。当如图20所示余弦振幅Acos被设定为零时,曲线形状与校正项(4)用于校正时相同。
实际NMR数据可能显示另外的失真,因为工具运动可能不是真正规则的,并且反演的回波串是最后一定数量的获取的滚动平均值。因此,通过计算具有不同周期的许多基函数并对它们求平均来生成增强的运动校正曲线形状。这样,余弦最大值变得更宽,并开始相互抵消。以下方程考虑对方程(6)的选定数量的周期时间求平均值:
方程(7)使用以下参数,Acos:余弦振幅;Tcos:余弦周期时间;Aexp:指数衰减振幅;以及Texp:指数衰减时间常数。
余弦周期时间可以由最小和最大周期或具有变化参数的平均周期定义。一种变化是使用频率而不是周期。可替代地,指数衰减项可以由任何其他衰减曲线形状代替。
为了减少要拟合的参数的数量,可以使用方程(7)的简化版本。例如,某些参数可以设置为零或固定值。替代地或另外地,可以将参数设置为固定关系。可能的固定关系是Texp=Tcos/3.5。
使用方程(6)的增强公式,图15中所示的上述实例数据中的高原被正确拟合。使用方程(6)的校正拟合曲线在图22中示出,并且不包括晚期运动伪影。如图23所示,晚期T2峰值消失。在图24中示出用于校正晚期运动伪影的校正拟合曲线的另一实例。
通常,具有短等待时间(例如,十分之几毫秒)的回波串可能不显示这些晚期运动伪影。因此,可以仅使用具有长极化时间(例如,6-20秒)的回波串,使用方程(6)和(7)的增强公式来校正一些晚期运动伪影。在第二步中,诸如在单独反演中,可以考虑短极化时间回波串。
可以理解,方程式(6)和(7)的增强公式可以更广泛地描述为如下的具有衰减项和周期项的方程。
f(t)=(周期项)+(衰减项)
可选地,这些增强的公式可以包括如下的周期项的周期的变化。
f(t)=(具有可变周期的周期项)+(衰减项)
当在识别晚期运动伪影的形状的最佳参数的过程中,可以考虑一些环境参数。例如,在存在冲失的情况下,振幅可以与的振动烈度相关。在存在冲失的情况下,平均周期可以与钻井转速(RPM)相关。周期变化可以与钻井转速变化相关,诸如由粘滑引起的那些钻井转速变化。因此,在一个或多个实施方案中,振动传感器14和/或转速传感器15可以向计算机处理系统12提供信号,其可以被配置来调整增强公式中的周期项的周期以与由振动传感器检测到的振动周期或由转速传感器检测到的转速变化周期一致。
如上所述,可以采用以不同极化时间获取的若干回波串的联合反演,并且期望改善结果。
图25是用于处理核磁共振(NMR)测量数据的方法250的流程图。框251要求利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据具有受运动伪影影响的回波串,其中运动伪影与由于NMR工具的运动而在感兴趣的体积中变化的磁场大小相关。在一个或多个实施方案中,NMR工具可以由钻柱输送通过穿透土地地层的钻孔。在一个或多个实施方案中,回波串包括具有大于或等于一秒的长极化时间的单个回波串。也就是说,在一个或多个实施方案中,一秒是用于区分短等待时间与长等待时间的阈值。在一个或多个实施方案中,回波串包括利用长极化时间TWET获得的第一回波串以及利用比TWET短的短极化时间TWTL获得的第二回波串。框252要求利用处理器通过使用校正反演方法来减小运动伪影对NMR测量数据的影响,所述校正反演方法对运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法包括乘法项,所述乘法项具有(i)具有至少一个局部最大值的项和(ii)可选衰减项。短语“具有至少一个局部最大值的项”旨在包括诸如余弦函数的周期函数或诸如一个周期或更短的周期函数的一部分。
方法250还可以包括从振动传感器接收振动信号并且调整周期时间以对应于振动的周期时间,振动传感器被配置来感测钻柱的振动。方法250还可以包括从转速传感器接收转速信号并且调整周期时间以对应于转速变化的周期时间,转速传感器被配置来感测钻柱的转速。
方法250还可以包括将NMR工具输送通过穿透土地地层的钻孔并且对土地地层执行NMR测量。方法250还可以包括使用与动作相关的设备以及从NMR测量导出的土地地层的性质对地层执行动作。在一个或多个实施方案中,性质是孔隙度。可以根据位置绘制孔隙度以提供孔隙度图。在一个或多个实施方案中,该动作是钻出具有使用孔隙度图选择的几何形状或迹线的钻孔,以便最大化油气生产。在这种情况下,动作相关设备可以是钻机。
以下阐述前述公开内容的一些实施方案:
实施方案1:一种用于处理核磁共振(NMR)测量数据的方法,所述方法包括:利用处理器接收从NMR工具获得的NMR测量数据,所述NMR测量数据包括受运动伪影影响的回波串,其中所述运动伪影与由于所述NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及利用所述处理器通过使用校正反演方法来减小所述运动伪影对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法包括乘法项,所述乘法项包括具有至少一个局部最大值的项和可选衰减项。
实施方案2:根据权利要求1所述的方法,其中所述运动伪影是晚期运动伪影。
实施方案3:根据权利要求2所述的方法,其中所述晚期运动伪影出现在所述回波串的后半部分。
实施方案4:根据权利要求1所述的方法,其中所述衰减项包括:
其中t是时间,Aexp是指数衰减振幅,并且Texp是指数衰减时间常数。
实施方案5:根据权利要求1所述的方法,其中所述周期项包括:
其中t是时间,Acos是余弦振幅,并且Tcos是余弦周期时间。
实施方案6:根据权利要求1所述的方法,其中所述乘法项包括:
其中t是时间,Acos是余弦振幅,Tcos是余弦周期时间,Aexp是指数衰减振幅,并且Texp是指数衰减时间常数。
实施方案7:根据权利要求1所述的方法,其还包括在不同的周期时间内对所述乘法项进行平均以提供平均乘法项。
实施方案8:根据权利要求7所述的方法,其中所述平均乘法项包括:
实施方案9:根据权利要求1所述的方法,其中所述回波串包括具有大于或等于一秒的长极化时间的单个回波串。
实施方案10:根据权利要求1所述的方法,其中所述回波串包括利用长极化时间TWET获得的第一回波串以及利用比TWET短的短极化时间TWTL获得的第二回波串。
实施方案11:根据权利要求1所述的方法,其中所述周期项的周期时间是恒定的。
实施方案12:根据权利要求1所述的方法,其中所述周期项的周期时间变化。
实施方案13:根据权利要求12所述的方法,其还包括从振动传感器接收振动信号并且调整所述周期时间以对应于所述振动的周期时间,所述振动传感器被配置来感测钻柱的振动。
实施方案14:根据权利要求12所述的方法,其还包括从转速传感器接收转速信号并且调整所述周期时间以对应于转速变化的周期时间,所述转速传感器被配置来感测钻柱的转速。
实施方案15:根据权利要求1所述的方法,其中所述NMR测量数据包括对土地地层的NMR测量,并且所述方法还包括使用动作相关设备和具有所述减小效果的所述NMR测量数据对土地地层执行动作。
实施方案16:根据权利要求1所述的方法,其中具有局部最大值的所述项是周期函数。
实施方案17:一种用于对土地地层执行核磁共振(NMR)测量的方法,所述方法包括:通过穿透所述土地地层的钻孔输送NMR工具;利用处理器接收从所述NMR工具获得的NMR测量数据,所述NMR测量数据包括受运动伪影影响的回波串,其中所述运动伪影与由于所述NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及利用所述处理器通过使用校正反演方法来减小所述运动伪影对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法包括乘法项,所述乘法项包括具有至少一个局部最大值的项和可选衰减项。
实施方案18:根据权利要求17所述的方法,其中所述NMR工具联接到钻柱。
实施方案19:一种用于处理核磁共振(NMR)测量数据的系统,所述系统包括:存储器,所述存储器具有计算机可执行指令;以及处理器,所述处理器被配置来执行所述计算机可执行指令,所述计算机可执行指令包括:接收从NMR工具获得的NMR测量数据,所述NMR测量数据包括受运动伪影影响的回波串,其中所述运动伪影与由于所述NMR工具的运动而在感兴趣的体积中变化的磁场大小相关;以及通过使用校正反演方法来减小所述运动伪影对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法包括乘法项,所述乘法项包括具有至少一个局部最大值的项和可选衰减项。
实施方案20:根据权利要求19所述的系统,其中所述NMR工具设置在钻柱上并且所述NMR测量数据包括对土地地层的NMR测量,所述系统还包括被配置来使用具有所述减小效果的所述NMR数据来对所述土地地层执行动作的动作相关设备。
实施方案21:根据权利要求19所述的系统,其还包括被配置来感测所述钻柱的振动的振动传感器,其中所述计算机可执行指令还包括调整所述周期时间以对应于所述振动的周期时间。
实施方案22:根据权利要求19所述的系统,其还包括被配置来感测所述钻柱的转速的转速传感器,其中所述计算机可执行指令还包括调整所述周期时间以对应于所述钻柱的转速振动的周期时间。
为了支持本文的教导,可以使用各种分析部件,包括数字和/或模拟系统。例如,井下电子设备11、计算机处理系统12或NMR工具10可以包括数字和/或模拟系统。系统可以具有部件,诸如处理器、存储介质、存储器、输入、输出、通信链路(有线、无线、脉冲泥浆、光学或其他)、用户接口、软件程序、信号处理器(数字或模拟)以及其他此类部件(诸如电阻器、电容器、电感器等)以便以本领域公知的若干方式中的任何一种方式提供对本文所公开的设备和方法的操作和分析。认为这些教导可以但不必结合存储在非暂态计算机可读介质上的一组计算机可执行指令来实现,所述非暂态计算机可读介质包括存储器(ROM、RAM)、光学(CD-ROM)或磁盘(磁盘、硬盘驱动器)或在执行时致使计算机实现本发明的方法的任何其他类型。除了本公开中描述的功能之外,这些指令还可以提供设备操作、控制、数据收集和分析以及系统设计者、所有者、用户或其他此类人员认为相关的其他功能。
此外,可以包括各种其他部件并且要求它们用于提供本文教导的各方面。例如,可以包括电源(例如,发电机、远程电源和电池中的至少一个)、冷却部件、加热部件、磁体、电磁体、传感器、电极、发送器、接收器、收发器、天线、控制器、光学单元、电气单元或机电单元以支持本文所讨论的各个方面或支持本公开之外的其他功能。
本文使用的术语“载体”表示任何装置、装置部件、装置组合、媒体和/或成员,其可以用于输送、容纳、支持或以其他方式促进另一装置、装置部件、装置组合、媒体和/或成员的使用。其他示例性非限制性载体包括盘管类型的钻柱、接合管类型的钻柱以及其任何组合或部分。其他载体实例包括套管、电缆、电缆探测器、钢丝探测器、落锤、井底组件、钻柱插入件、模块、内部外壳及其基板部分。
本文描述的流程图仅是实例。可以存在对本文描述的这些图或步骤(或操作)的许多变化而不背离本发明的精神。例如,可以按不同的顺序执行所述步骤,或可以增加、删除或修改步骤。所有这些变化都被认为是要求保护的发明的一部分。
已经用冠词“一”或“一个”引入了实施方案的元件。冠词旨在表示存在一个或多个元件。术语“包括”和“具有”旨在是包括性的,使得除了所列要素之外,可能有另外的要素。当与至少两个术语的列表一起使用时,连词“或”旨在表示任何术语或术语组合。术语“第一”、“第二”等不表示特定顺序,而是用于区分不同元件。
虽然已经示出和描述了一个或多个实施方案,但是可以对其进行修改和替换,而不背离本发明的精神和范围。因此,应理解,已经通过说明而非限制的方式描述了本发明。
将认识到,各种部件或技术可以提供某些必要或有益的功能或特征。因此,支持所附权利要求及其变体所需的这些功能和特征被认为是作为本文教导的一部分和所公开的本发明的一部分而固有地包括在内。
虽然已参考示例性实施方案描述本发明,但将了解,可以做出各种改变并且等效物可替代本发明的各元件,而不背离本发明的范围。另外,将理解许多修改以使特定仪器、情况或材料适应本发明的教导,而不背离本发明的基本范围。因此,希望本发明不限于作为用于实施本发明所考虑的最佳方式所公开的具体实施方案,而是本发明将包括落在所附权利要求书的范围内的所有实施方案。
Claims (15)
1.一种用于处理核磁共振(NMR)测量数据的方法(140、250),所述方法(140、250)的特征在于:
利用处理器接收从NMR工具(10)获得的NMR测量数据,所述NMR测量数据包括受运动伪影(151)影响的回波串(150),其中所述运动伪影(151)与由于所述NMR工具(10)的运动而在感兴趣的体积中变化的磁场大小相关;以及
利用所述处理器通过使用校正反演方法(140、250)来减小所述运动伪影(151)对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影(151)进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法(140、250)包括乘法项,所述乘法项包括具有至少一个局部最大值的项和可选衰减项。
2.根据权利要求1所述的方法(140、250),其中所述运动伪影(151)是晚期运动伪影(151)。
3.根据权利要求1所述的方法(140、250),其中所述衰减项包括:
其中t是时间,Aexp是指数衰减振幅,并且Texp是指数衰减时间常数。
4.根据权利要求1所述的方法(140、250),其中所述周期项包括:
其中t是时间,Acos是余弦振幅,并且Tcos是余弦周期时间。
5.根据权利要求1所述的方法(140、250),其中所述乘法项包括:
其中t是时间,Acos是余弦振幅,Tcos是余弦周期时间,Aexp是指数衰减振幅,并且Texp是指数衰减时间常数。
6.根据权利要求1所述的方法(140、250),其还包括在不同的周期时间内对所述乘法项进行平均以提供平均乘法项。
7.根据权利要求6所述的方法(140、250),其中所述平均乘法项包括。
8.根据权利要求1所述的方法(140、250),其中所述回波串(150)包括具有大于或等于一秒的长极化时间的单个回波串(150)。
9.根据权利要求1所述的方法(140、250),其中所述回波串(150)包括利用长极化时间TWET获得的第一回波串(150)以及利用比TWET短的短极化时间TWTL获得的第二回波串(150)。
10.根据权利要求1所述的方法(140、250),其中所述周期项的周期时间变化。
11.根据权利要求10所述的方法(140、250),其还包括从振动传感器(14)接收振动信号并且调整所述周期时间以对应于所述振动的周期时间,所述振动传感器(14)被配置来感测钻柱(6)的振动。
12.根据权利要求10所述的方法(140、250),其还包括从转速传感器(15)接收转速信号并且调整所述周期时间以对应于转速变化的周期时间,所述转速传感器(15)被配置来感测钻柱(6)的转速。
13.根据权利要求1所述的方法(140、250),其中所述NMR测量数据包括对土地地层(4)的NMR测量,并且所述方法(140、250)还包括使用动作相关设备和具有所述减小效果的所述NMR测量数据对土地地层(4)执行动作。
14.根据权利要求1所述的方法(140、250),其还包括通过穿透所述土地地层(4)的钻孔(2)输送所述NMR工具(10)。
15.一种用于处理核磁共振(NMR)测量数据的系统,所述系统的特征在于:
存储器,所述存储器具有计算机可执行指令;以及
处理器,所述处理器被配置来执行所述计算机可执行指令,所述计算机可执行指令包括:
接收从NMR工具(10)获得的NMR测量数据,所述NMR测量数据包括受运动伪影(151)影响的回波串(150),其中所述运动伪影(151)与由于所述NMR工具(10)的运动而在感兴趣的体积中变化的磁场大小相关;以及
通过使用校正反演方法(140、250)来减小所述运动伪影(151)对所述NMR测量数据的影响,所述校正反演方法对所述运动伪影(151)进行建模以提供校正的横向弛豫时间常数(T2)分布,所述校正反演方法(140、250)包括乘法项,所述乘法项包括具有至少一个局部最大值的项和可选衰减项。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/089893 | 2016-04-04 | ||
US15/089,893 US10429536B2 (en) | 2016-04-04 | 2016-04-04 | T2 inversions with reduced motion artifacts |
PCT/US2017/024998 WO2017176548A1 (en) | 2016-04-04 | 2017-03-30 | T2 inversions with reduced motion artifacts |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109154189A true CN109154189A (zh) | 2019-01-04 |
CN109154189B CN109154189B (zh) | 2022-08-26 |
Family
ID=59961451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780029808.4A Active CN109154189B (zh) | 2016-04-04 | 2017-03-30 | 具有减少的运动伪影的t2反演 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10429536B2 (zh) |
EP (1) | EP3440312B1 (zh) |
CN (1) | CN109154189B (zh) |
BR (1) | BR112018070337B1 (zh) |
SA (1) | SA518400163B1 (zh) |
WO (1) | WO2017176548A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11422283B1 (en) * | 2021-07-14 | 2022-08-23 | Halliburton Energy Services, Inc. | Reducing motion effects on nuclear magnetic resonance relaxation data |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8970217B1 (en) | 2010-04-14 | 2015-03-03 | Hypres, Inc. | System and method for noise reduction in magnetic resonance imaging |
US10739489B2 (en) * | 2016-01-15 | 2020-08-11 | Baker Hughes, A Ge Company, Llc | Low gradient magnetic resonance logging for measurement of light hydrocarbon reservoirs |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101063712A (zh) * | 2005-12-20 | 2007-10-31 | 西门子公司 | 用于在磁共振数据中确定横向松弛时间t2*的方法 |
CN102608664A (zh) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | 深度维核磁共振反演获取横向弛豫时间谱的方法及装置 |
CN103080764A (zh) * | 2010-08-30 | 2013-05-01 | 皇家飞利浦电子股份有限公司 | 利用横向弛豫时间测量数据和质子共振频率偏移数据的储水组织和脂肪组织的mri热成像 |
KR20130051084A (ko) * | 2011-11-09 | 2013-05-20 | 삼성전자주식회사 | 고차 확산자기공명 영상의 인공물 저감 장치 및 방법 |
WO2015031149A1 (en) * | 2013-08-25 | 2015-03-05 | Schlumberger Canada Limited | Correction of motion effect in nuclear magnetic resonance (nmr) logging |
US20150241541A1 (en) * | 2014-02-26 | 2015-08-27 | Baker Hughes Incorporated | T2 inversions with reduced motion artifacts |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5428291A (en) | 1993-07-01 | 1995-06-27 | Exxon Research And Engineering Company | Determination of fluid transport properties in porous media by nuclear magnetic resonance measurements of fluid flow |
US6268727B1 (en) | 1997-06-24 | 2001-07-31 | Southwest Research Institute | Measurement of flow fractions flow velocities and flow rates of a multiphase fluid using ESR sensing |
US6069477A (en) | 1997-09-05 | 2000-05-30 | Western Atlas International, Inc. | Method for improving the accuracy of NMR relaxation distribution analysis with two echo trains |
US6566874B1 (en) | 1998-07-30 | 2003-05-20 | Schlumberger Technology Corporation | Detecting tool motion effects on nuclear magnetic resonance measurements |
US6459263B2 (en) | 2000-02-08 | 2002-10-01 | Baker Hughes Incorporated | Nuclear magnetic resonance measurements in well logging using motion triggered pulsing |
US6686738B2 (en) | 2001-04-17 | 2004-02-03 | Baker Hughes Incorporated | Method for determining decay characteristics of multi-component downhole decay data |
BRPI0414998A (pt) | 2003-10-03 | 2006-11-21 | Halliburton Energy Serv Inc | métodos para identificar gás em uma formação geológica, para analisar formações geológicas, e de rmn para analisar formações geológicas, e sistema |
US7268547B2 (en) | 2003-10-07 | 2007-09-11 | Baker Hughes Incorporated | Correction of motion influences in NMR signals |
US7180287B2 (en) | 2004-05-05 | 2007-02-20 | Baker Hughes Incorporated | Reduction of motion artifacts in NMR |
US7196516B2 (en) | 2004-08-16 | 2007-03-27 | Baker Hughes Incorporated | Correction of NMR artifacts due to constant-velocity axial motion and spin-lattice relaxation |
US20060158184A1 (en) | 2005-01-18 | 2006-07-20 | Baker Hughes Incorporated | Multiple echo train inversion |
US20070022244A1 (en) * | 2005-07-25 | 2007-01-25 | Honeywell International Inc. | Methods and systems for refresh and error scrubbing of dynamic memory devices |
US9239367B2 (en) | 2012-09-14 | 2016-01-19 | General Electric Company | Intra-blade filter for motion corrected magnetic resonance data |
US9983288B2 (en) | 2014-04-24 | 2018-05-29 | Ohio State Innovation Foundation | Free-breathing myocardial T2* mapping using single-heartbeat, multi-shot, gradient echo-echo planar imaging (GRE-EPI) and automatic non-rigid motion correction |
-
2016
- 2016-04-04 US US15/089,893 patent/US10429536B2/en active Active
-
2017
- 2017-03-30 WO PCT/US2017/024998 patent/WO2017176548A1/en active Application Filing
- 2017-03-30 EP EP17779559.8A patent/EP3440312B1/en active Active
- 2017-03-30 BR BR112018070337-0A patent/BR112018070337B1/pt active IP Right Grant
- 2017-03-30 CN CN201780029808.4A patent/CN109154189B/zh active Active
-
2018
- 2018-10-03 SA SA518400163A patent/SA518400163B1/ar unknown
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101063712A (zh) * | 2005-12-20 | 2007-10-31 | 西门子公司 | 用于在磁共振数据中确定横向松弛时间t2*的方法 |
CN103080764A (zh) * | 2010-08-30 | 2013-05-01 | 皇家飞利浦电子股份有限公司 | 利用横向弛豫时间测量数据和质子共振频率偏移数据的储水组织和脂肪组织的mri热成像 |
KR20130051084A (ko) * | 2011-11-09 | 2013-05-20 | 삼성전자주식회사 | 고차 확산자기공명 영상의 인공물 저감 장치 및 방법 |
CN102608664A (zh) * | 2012-02-17 | 2012-07-25 | 中国石油大学(北京) | 深度维核磁共振反演获取横向弛豫时间谱的方法及装置 |
WO2015031149A1 (en) * | 2013-08-25 | 2015-03-05 | Schlumberger Canada Limited | Correction of motion effect in nuclear magnetic resonance (nmr) logging |
US20150241541A1 (en) * | 2014-02-26 | 2015-08-27 | Baker Hughes Incorporated | T2 inversions with reduced motion artifacts |
Non-Patent Citations (1)
Title |
---|
JACK HORKOWITZ ET AL: "APPLICATIONS OF A NEW MAGNETIC RESONANCE LOGGING- WHILE-DRILLING TOOL IN A GULF OF MEXICO DEEPWATER DEVELOPMENT PROJECT", 《SPWLA 43 RA ANNUAL LOGGING SYMPOSIUM》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11422283B1 (en) * | 2021-07-14 | 2022-08-23 | Halliburton Energy Services, Inc. | Reducing motion effects on nuclear magnetic resonance relaxation data |
Also Published As
Publication number | Publication date |
---|---|
WO2017176548A1 (en) | 2017-10-12 |
EP3440312A1 (en) | 2019-02-13 |
BR112018070337B1 (pt) | 2023-02-14 |
US20170285214A1 (en) | 2017-10-05 |
SA518400163B1 (ar) | 2023-01-02 |
EP3440312B1 (en) | 2021-12-22 |
BR112018070337A2 (pt) | 2019-01-29 |
EP3440312A4 (en) | 2019-12-11 |
CN109154189B (zh) | 2022-08-26 |
US10429536B2 (en) | 2019-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9671483B2 (en) | T2 inversions with reduced motion artifacts | |
US6255818B1 (en) | Method and apparatus for performing magnetic resonance measurements | |
US6466013B1 (en) | Nuclear magnetic resonance measurements in well logging using an optimized rephasing pulse sequence | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
US7268547B2 (en) | Correction of motion influences in NMR signals | |
Grunewald et al. | Multiecho scheme advances surface NMR for aquifer characterization | |
NO337868B1 (no) | Fremgangsmåte og anordning for behandling av kjernemagnetiske resonans-signaler fra en grunnformasjon | |
WO2015031149A1 (en) | Correction of motion effect in nuclear magnetic resonance (nmr) logging | |
NO336430B1 (no) | Bestemmelse av T1 relaksasjonstid fra multippel ventetid NMR logger frembragt i den samme eller forskjellige passeringer | |
CN109072693A (zh) | 运动检测以及磁共振数据的校正 | |
CN109154189A (zh) | 具有减少的运动伪影的t2反演 | |
US10551521B2 (en) | Magnetic resonance pulse sequences and processing | |
CA2570147C (en) | Phase-alternated carr-purcell nmr echo sequence | |
US10725195B2 (en) | Methods of estimating borehole and formation properties using an electromagnetic induction logging tool having random tool decenter positions during data acquisition | |
US9223048B2 (en) | System and method to detect a fluid flow without a tipping pulse | |
US10895659B2 (en) | Detecting nuclear magnetic resonance logging tool motion | |
Grombacher et al. | Frequency cycling to alleviate unknown frequency offsets for adiabatic half-passage pulses in surface nuclear magnetic resonance | |
US11899158B2 (en) | Logging tool motion error reduction for nuclear magnetic resonance logging via pulse sub-sequences | |
US20230288600A1 (en) | Logging tool motion error reduction for nuclear magnetic resonance logging via pulse angle selection | |
US20230384472A1 (en) | Data Inversion to Reduce Motion Effects on Nuclear Magnetic Resonance Relaxation Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: Texas in the United States Applicant after: BAKER HUGHES Inc. Address before: Texas, USA Applicant before: BAKER HUGHES Inc. |
|
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |