CN1186649C - 用于计算自旋-自旋弛豫时间分布的设备和方法 - Google Patents
用于计算自旋-自旋弛豫时间分布的设备和方法 Download PDFInfo
- Publication number
- CN1186649C CN1186649C CNB99123507XA CN99123507A CN1186649C CN 1186649 C CN1186649 C CN 1186649C CN B99123507X A CNB99123507X A CN B99123507XA CN 99123507 A CN99123507 A CN 99123507A CN 1186649 C CN1186649 C CN 1186649C
- Authority
- CN
- China
- Prior art keywords
- spin
- rightarrow
- relaxation time
- magnetic field
- alpha
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000003068 static effect Effects 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 7
- 239000013598 vector Substances 0.000 abstract description 5
- 238000001228 spectrum Methods 0.000 abstract description 4
- 238000002592 echocardiography Methods 0.000 abstract description 3
- 238000000354 decomposition reaction Methods 0.000 abstract description 2
- 230000010354 integration Effects 0.000 abstract description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 41
- 238000005259 measurement Methods 0.000 description 27
- 238000005553 drilling Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 12
- 239000012530 fluid Substances 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 10
- 230000035699 permeability Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 8
- 210000004940 nucleus Anatomy 0.000 description 7
- 239000004215 Carbon black (E152) Substances 0.000 description 6
- 229930195733 hydrocarbon Natural products 0.000 description 6
- 150000002430 hydrocarbons Chemical class 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 5
- 229910052739 hydrogen Inorganic materials 0.000 description 5
- 239000001257 hydrogen Substances 0.000 description 5
- 230000010287 polarization Effects 0.000 description 5
- 238000009738 saturating Methods 0.000 description 5
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000000295 complement effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 239000011734 sodium Substances 0.000 description 4
- 229910052708 sodium Inorganic materials 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 239000004020 conductor Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000005389 magnetism Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 229910000889 permalloy Inorganic materials 0.000 description 3
- 230000010363 phase shift Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000005669 field effect Effects 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000003019 stabilising effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229910000859 α-Fe Inorganic materials 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 229910001030 Iron–nickel alloy Inorganic materials 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- KEAYESYHFKHZAL-IGMARMGPSA-N sodium-23 atom Chemical compound [23Na] KEAYESYHFKHZAL-IGMARMGPSA-N 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000000264 spin echo pulse sequence Methods 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
Images
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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Geophysics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及用于计算自旋-自旋驰豫时间分布的一种方法。所说自旋回波幅值利用一个时间窗口上的接收器电压的硬件积分获得。利用一个线性算子将一个驰豫时间分布与自旋回波相映射,产生所说线性算子的单值分解(SVD),确定所说SVD的矢量,和利用所说矢量压缩所说自旋回波数据。为了消除遥测瓶颈,在井下计算所说T2型谱,并传输到地面。
Description
技术领域
一般来说,本发明涉及用于测量被一个井孔贯穿的地层构造的核磁共振特性的一种设备和方法,更具体地说,本发明涉及用于计算自旋-自旋弛豫时间分布的设备和方法。
背景技术
人们已经清楚地认识到,具有非零核自旋磁矩的地层构造中的原子粒子,例如质子,具有与施加在该构造上的静磁场沿同向排列的趋势。这种磁场可以自然产生,例如地磁场BE。用于施加垂直于BE的第二磁场的射频脉冲在横向平面(垂直于BE)中产生一个磁化分量,该分量以被称为拉莫尔频率ωL的特征谐振频率围绕BE矢量进动,所说拉莫尔频率的高低取决于静磁场的强度和粒子的旋磁比。围绕例如0.5高斯的磁场BE进动的氢原子核(质子)的特征频率大约为2k赫兹。如果使许多氢原子核同相地进动,这些质子的组合磁场能够产生在接收器线圈中产生可检测到的振荡电压,形成本领域技术人员称之为自由感应衰减或自旋回波的状态。存在于岩石孔隙中的水和碳氢化合物的氢原子核产生与由其它固体物质产生的信号不同的核磁共振(NMR)信号。
授予Taicher等人的美国专利US-4717878和授予Kleinberg等人的美国专利US-5055787都记载了NMR测井设备,这些设备应用永磁体使氢原子核极化和产生静磁场B0,应用射频天线激励和探测核磁共振信号以确定一种构造的孔隙率、自由流体比率、和渗透性。原子核以时间常数T1与所施加的磁场B0同向排列。在经过一段时间的极化之后,通过施加垂直于所说静磁场B0的一个射频场B1,可以改变核磁化方向与所施加磁场方向之间的夹角,所说静磁场的拉莫尔频率fL=γB0/2π,其中γ为质子的旋磁比,B0表示静磁场强度。在射频脉冲结束之后,质子在垂直于B0的平面内进动。一个再聚焦射频脉冲序列产生一个自旋回波序列,这些自旋回波在天线中产生可以检测到的NMR信号。
授予Melvin Miller的美国专利US-5280243记载了用于在钻探过程中进行构造评估的一种核磁共振测井仪。这种测井仪包括由设置在沿钻轴外侧纵向延伸的环形凹槽中的一个永磁体和设置在钻轴外侧一个非导磁性套管上的一个天线构成的一个探测器部分。静磁场量值梯度沿径向方向。天线产生基本同时垂直于测井仪纵轴和静磁场方向的射频磁场。在专利‘243所述装置中,磁体的轴向尺寸与其直径相比必须较长,以使所产生的磁场的特性近似于一个2维偶极子。
授予Taicher等人的美国专利US-5757186记载了一种即钻即测型测井仪,这种测井仪包括用于对地层构造进行核磁共振测量的一个检测装置。所说NMR检测装置安装在形成于钻轴外表面的一个环形凹槽中。在一个实施例中,在凹槽中***一个磁力线封闭罩,在磁力线封闭罩的外径向表面上设置有一个磁体。该磁体是由多个径向段构成的,这些磁体段是从测井仪的纵轴径向向外磁化的。磁力线封闭罩用于使磁场具有适合的取向。
在专利‘243和‘186中公开的测井仪存在同样的问题:两者都需要使用非导电磁体,并将磁体设置在钻轴外侧。在专利‘243的测井仪中,钻轴的外表面必须包含一个凹进区域以容纳非导电磁体。在专利‘186的测井仪中,钻轴的外表面必须包含一个凹进区域以容纳磁力线封闭罩、非导电磁体、和天线。因为钻轴的强度是其半径的函数,为仅仅容纳磁体或磁力线封闭罩和天线而减小钻轴外径会导致钻轴具有不可接受的薄弱部分,这一部分在钻探过程中可能发生弯曲或折断。
授予Kleinberg等人的美国专利US-5557201记载了用于在钻探过程中进行构造评估的一种脉冲核磁测井仪。这种测井仪包括一个钻头、一根钻杆柱、和安装在由非磁性合金制成的钻轴内的一个脉冲核磁共振装置。这种测井仪包括形成在钻杆柱和脉冲NMR装置内部的一条通道,钻探泥浆穿过该通道泵入井孔中。所说脉冲NMR装置包括两个管形磁体,它们按照相同磁极彼此面对的方式安装,并且围绕着所说通道,在钻杆柱外表面上位于两个磁体之间安装有一个天线线圈。这种测井仪用于使被本领域技术人员称之为鞍点的测量区域内的原子核发生共振。
授予Sezginer等人的美国专利US-5705927也记载了用于在钻探过程中进行构造评估的一种脉冲核磁测井仪。这种测井仪包括设置在测井仪内部或外部的调整磁体,该磁体通过增大井孔中静磁场的强度抑制井孔流体的磁共振信号,从而使井孔中的拉莫尔频率大于由设置在测井仪凹进部分中的一个射频天线所产生的振荡场的频率。这个调整磁体还减小了探测区域中静磁场的梯度。
发明内容
借助于本发明用于计算自旋-自旋弛豫时间分布的一种设备和方法可以克服现有技术的上述缺点。在井孔所贯穿的构造中施加一个基本轴对称的静磁场。在所说构造中还施加一个振荡磁场。探测所说构造的核磁共振信号并传输到位于井孔中的信号处理器中。所说信号处理器根据所探测的信号计算自旋-自旋弛豫时间分布。可以将所说自旋-自旋弛豫时间传输到所说井孔表面(上孔)。
所说多个信号的幅值为信号加噪声幅值,Aj,它可以用下述关系式表征:
其中ηj为测量值Aj中的噪声,ai为在T2,j时刻所取的T2分布幅值, 表示矩阵X的矩阵元,
其中tw为等待时间,c为一个常数(T1/T2的比值),T1是纵向弛豫时间,T2是自旋-自旋弛豫时间,是Δt为回波间隔,而j=1,2,...N,其中N为在一个脉冲序列中收集的回波数量。用矩阵符号表示,该方程变为
因为噪声η是未知的,所以通过求解本发明定义的函数 的最小值可以近似求得
可以在该函数中加入一调整项,并使用一种适合的迭代求最小值算法求出函数
的最小值。
从以下对于附图的说明可以清楚地了解本发明的优点。应当理解,这些附图仅仅用于说明,而不是用于限定本发明。
附图说明
图1表示一种即钻即测型设备;
图2表示所说低梯度探测器;
图3a-3d表示对应于四个低梯度磁体配置的梯度|B0|的等高线;
图4表示所说高梯度探测器;
图4b表示对应于高梯度磁体配置的梯度|B0|的等高线;
图5表示简单数据采集模式;
图6表示交叉存取数据采集模式;
图7表示脉冲串数据采集模式;和
图8为表示脉冲程序装置的方框示意图。
具体实施方式
参见图1,图中表示了一种核磁共振(NMR)即钻即测型测井仪10。所说测井仪10包括一个钻头12、钻杆柱14、多个射频天线36、38,和至少一个梯度线圈56。所说测井仪10还包括封装在钻轴22内的电路20。所说电路20包括天线36、38的射频谐振电路、一个微处理器、一个数字信号处理器、和一条低压总线。所说测井仪10还包括多个管形磁体30、32、和34,它们都沿平行于测井仪10纵轴的方向极化,但是彼此相对,即以相同的磁极彼此面对。磁体30、32和34可以由导电材料或是非导电材料制成。磁体30、32、和34以及天线36、38的配置形成至少两个具有基本轴对称的静磁场和射频磁场的NMR探测区域60、62。
用于在所说构造中钻井的装置24包括钻头12和钻轴22。所说钻轴22可以包括用于在钻井过程中使测井仪在井孔内的径向运动稳定的一个稳定装置(未示出),但是,所说稳定装置并不是必需的;所以,测井仪10可以在稳定运动或不稳定运动状态下工作。泥浆流动套管28构成了用于运载钻探流体通过钻杆柱14的一条通道90。一个驱动机构26转动所说钻头12和钻杆柱14。在授予Clark等人的美国专利US-4949045中对于这种驱动机构给予了充分的介绍。但是,也可以在钻杆柱中设置一个井下泥浆马达作为驱动机构26。
本发明的实现方案之一是组合N+1个磁体以便在所说地层构造中获得至少N个探查区域。本发明设想的组合形式包括,但是不限于,低梯度-低梯度、高梯度-高梯度、高梯度-低梯度、低梯度-高梯度组合,或者是高梯度、低梯度、和鞍点区域的组合。所说构造中高梯度和低梯度静磁场区域的组合具有若干优点。例如,高梯度区域可能具有较高的信噪比,但是当测井仪10在井孔中作横向运动时可能会产生信号损失。另一方面,当测井仪10处于运动状态时,低梯度区域对于信号损失问题的敏感度较低。而且,在测井仪运动速度适中的情况下,在低梯度区域中可以获得比在高梯度区域中更长的回波脉冲串,从而提供较多的有关渗透率、约束流体和自由流体、以及碳氢化合物类型的信息。此外,将在两种梯度区域获取的数据结合可以提供有关测井仪10的横向运动量的定量信息,并且可以用于对NMR数据进行运动校正,或者,至少可以控制数据的质量。可以将诸如应变仪、加速度计、或磁强计,或这些装置的任意组合构成的装置的测量结果与NMR信息结合在一起以控制数据的质量或者对自旋回波串进行校正。在将高梯度静磁场区域和低梯度静磁场区域结合的情况下,高梯度区域具有较强的扩散效应,所以在确定碳氢化合物类型技术方面比低梯度区域更为有益。最后,低梯度区域的静磁场具有较低的强度量值,所以这个区域由于具有较低的拉莫尔频率而较少受到构造和井孔中流体导电性的影响。
低梯度探测器
参见图2,在下文中称之为低梯度探测器的测井仪的一个部分中,一个中部磁体30与一个下部磁体32沿轴向分开。这些磁体30、32产生其极化方向沿径向的一个基本轴对称的静磁场,而且在合理长的圆柱体形壳层内,所说静磁场具有十分恒定的量值。本发明的实现方案之一是激励构造中分别以不同的射频频率谐振的多个圆柱体形自旋壳层,和用射频脉冲序列顺序地探询每个壳层。
磁体30、32之间的区域适合安置诸如电子器件、射频天线、和其它类似部件。例如,多个电路盒70可以构成泥浆套管28的一个组成部分。这些电路盒70可以安置射频电路(例如,Q-开关、双工器、和前置放大器),优选的是贴近所说射频天线设置。在本发明的一个优选实施例中,所说电路盒70构成透磁部件16的一个组成部分。在这种情况下,为了保持磁场的轴对称性,在每个电路盒70之上设置一个高透磁性的盖罩72。
所说透磁部件16设置在钻轴22内部位于磁体30、32之间。部件16可由单个部件构成或者由在磁体30、32之间的多个部分组合而成。部件16由适合的透磁材料,例如铁氧体、透磁钢或其它铁镍合金、抗腐蚀透磁钢、或在部件设计中用作结构部件的透磁钢,如15-5Ph不锈钢制成。所说透磁部件16使磁场集中,并且也可以运载钻探泥浆通过所说钻杆柱或者为所说钻轴柱提供结构支撑。此外,部件16改善了由磁体30、32产生的静磁场的形状,和使由于在获取NMR信号过程中测井仪的垂直和横向运动所引起的静磁场的变化最小。在磁体30、32之间的管段28可以包括透磁部件16。在这种情况下,位于磁体30、32之下的各个套管段应当由非磁性部件构成。或者,可以由磁体30、32之间包围着管段28的透磁性支架构成部件16。在这种情况下,所说管段可以由磁性或非磁性材料制成。本发明的实现方案之一是将所说支架和管段结合构成部件16。
磁体30、32沿平行于测井仪10纵轴的方向极化,并以相同的磁极彼此面对。对于每个磁体30、32来说,感生磁力线沿着测井仪10的轴向从磁体30、32的一端出发进入地层构造中,并传播进入磁体30、32的另一端。在中部磁体30和下部磁体32之间的区域中,感生磁力线从中心向外传播到构造中,沿基本垂直于测井仪10轴向的一个方向生成一个静磁场。然后,感生磁力线对称地从中部磁体30上方和下部磁体32下方传播进去,并沿纵向会聚在套管28内部。由于磁体之间存在间隔,在中部磁体30与下部磁体32之间的中部区域中的静磁场强度与鞍点区域相比在空间上是均匀的。
磁体30、32之间的间隔量是根据以下几个因素确定的:(1)选择必需的磁场强度和均匀特性;(2)产生在探测区域具有较小径向变化的场,从而使在脉冲序列期间(即CPMG、CPI或其它序列)接收的回波对于测井仪横向运动的敏感度较低;(3)探查深度;和(4)使谐振电路与低压遥测总线之间的干扰最小,以便改善检测NMR信号的接收天线与构造之间的隔离。随着磁体30、32之间的间隔减小,磁场强度逐渐增大,均匀性变差。相反,随着磁体30、32之间间隔增大,磁场强度逐渐减小,并且变得更加均匀。
图2a-2d表示对应于中部磁体30和下部磁体32的四种实验室模拟配置的
的等高线。使用具有预定直径(在模拟所有配置时使用不变的直径)的测井仪计算出这些模拟结果。对应于图2a的配置包括使中部磁体30与下部磁体32分开25英寸的一个非透磁性部件。对应于图2b的配置包括使中部磁体30与下部磁体分开18英寸的一个非透磁性部件。对应于图2c的配置包括使中部磁体30和下部磁体32分开8英寸的一个非透磁性部件。对应于图2d的低梯度探测器包括使中部磁体30与下部磁体32分开25英寸的一个透磁部件16。上述尺寸是仅仅为了说明距离和/或透磁或非透磁材料对于
的影响而模拟的。图3a-3d分别表示对应于图2a-2d所示配置的梯度|B0|的等高线。
在低梯度探测器中,透磁部件16使大部分的磁通量进入测井仪10的中部部分。举例来说,图2d中所示B0场在径向距离测井仪10纵轴大约7英寸处的量值是图2a中所示B0场的两倍,图2a中的B0场是利用由非透磁材料部件隔开的相同磁体配置产生的。此外,所说低梯度探测器沿径向产生更长和更均匀的静磁场。在测井仪的这个部分测得的NMR信号对于测井仪的垂直运动的敏感度是相当小的。参见图3d,利用所说低梯度探测器,在径向距离测井仪纵轴大约7英寸的距离处测得相对较小,大约3高斯/厘米的梯度。这种低梯度使得所测量的NMR信号对于测井仪10的横向运动的敏感度相当小。当运动速度适中时,可以在这个区域获得较长的回波串,从而提供更多的与渗透率、约束和自由流体、和碳氢化合物类型有关的信息。在低梯度探测器的情况下,当采用其它梯度配置时,测井仪10周围的富质子井孔区域仅仅以比施加到探测区域中的频率较高的频率谐振,亦即,没有质子井孔信号。当以相同的射频频率激励时,在钻探泥浆中存在的其它NMR敏感原子核如钠-23与氢相比需要在高得多的静磁场强度下发生谐振。对于低梯度探测器,在可以探测到这种不想要的信号的测井仪周围的井孔区域或接近天线处不产生这种较高的磁场强度。
高梯度探测器
参见图4,在下文中称之为高梯度探测器的测井仪的另一部分中,一个中部磁体30与一个上部磁体34沿轴向分开。磁体30、34沿平行于测井仪10的纵轴的方向极化,并以相同磁极彼此面对。这些磁体30、34产生一个基本轴对称静磁场,其极化方向是径向的,并且在合理长的圆柱体形壳层上,所说静磁场具有相当恒定的量值。本发明的实现方案之一是激励构造中多个圆柱体形自旋壳层,其中每个壳层以不同的射频频率谐振。
如图2c所示,如果磁体30与磁体34之间的间隔大约为8英寸,则静磁场强度的等高线基本是直线,并且
的强度大于低梯度区域的静磁场强度。但是,如图3c所示,梯度|B0|在径向距离测井仪纵轴大约7英寸距离处变大。|B0|的等高线弯曲的,表示梯度沿轴向发生变化。
通过在磁体30、34之间***一个透磁部件16可以改进所说高梯度探测器。图4a表示对应于其中透磁性部件16将所说上部磁体34与中部磁体30分开8英寸的一种配置的
的等高线。图4a所示的等高线表示一个强度略大的磁场,它比图2c所示等高线具有较好的信噪比和沿轴向较小的弯曲度。此外,如图4b所示,所说透磁部件16在轴向产生更为恒定的梯度|B0|,这使得对于受到扩散影响的NMR测量结果的分析简单了。
在高梯度探测器的情况下,当采用其它梯度配置结构时,测井仪10周围富质子井孔区域仅仅以比在探查区域所施加频率高的频率谐振,即没有质子井孔信号。所说高梯度探测器对于井孔流体中少量的钠是敏感的。对于30%NaCl浓度的井孔流体,或许是最不利的情况,由于钠信号造成的孔隙率测量值的误差大约为0.08pu。在低梯度探测器中,钠信号大大小于在高梯度探测器中。因此,对于两种NMR探测器来说钠信号都可以忽略不计。
天线和梯度线圈
参见图2和图4,利用设置在凹进部分50、52中的天线36、38在探测区域中产生一个射频磁场。所说射频磁场可以利用一个或多个射频天线段产生,这些射频天线在测井设备的不同圆周扇面中发射和接收。参见转让给Schlumberger技术公司的美国专利申请No.08/880343和09/094201(律师档案编号No.24784和24784IP)。可取的是,每个天线36、38都包括圆周状缠绕在凹进区域50、52上的一个线圈18。由这样一种线圈结构产生的射频场是基本轴对称的。本发明的实现方案之一是利用天线36、38探测NMR信号。但是,也可以使用独立的天线或接收器探测信号。在凹进区域50、52中天线36、38之下填充有非导电材料54。可取的是,材料54为一种铁氧体,以提高天线36、38的效率。或者,材料54可以包含塑料、橡胶、或加强的环氧树脂合成材料。天线36、38在射频电路激励下谐振,以在探测区域中生成一个射频磁场。
所说凹进区域52在钻轴中形成一个浅槽,而不减少钻轴的内径,这通常用于增加钻轴上外径凹进用以设置天线的区域的强度。凹进区域50的深度大于凹进区域52。由于机械方面的限制,只能在钻轴内径大大减小的位置具有一个深凹区域。本发明的实现方案之一是凹进区域50、52具有基本相同的深度,或者凹进区域52具有大于区域50的深度。
探查区域中圆柱体形自旋壳层可以利用设置在所说凹进区域50和/或52中的至少一个具有方向敏感性的梯度线圈56沿轴向,或者可取地,沿方位角分段。在本发明的一个优选实施例中,三个梯度线圈呈圆周状设置在凹进区域周围,并且彼此之间的角距离为120°。可以采用其它数量的梯度线圈,少于或多于三个,并且这些线圈可以分开不是120°的角度和/或呈非相等角度分布。每个线圈56都由导线环路构成,它们与材料54的外表面曲率一致。由各个梯度线圈56在朝向所说线圈的构造区域中产生的磁场基本平行于由所说磁体产生的静磁场。
如本领域技术人员所知的,在基本NMR测量中,需要向探查构造施加一个脉冲序列。在授予Abdurrahman Sezginer的美国专利US-5596274中和授予Kleinberg等人的美国专利US-5023551中,一个脉冲序列,诸如Carr-Purcell-Meiboom-Gill(CPMG)序列,首先向所说构造施加一个激励脉冲,即一个90°脉冲,该脉冲使所说自旋旋转到横向平面。在所说自旋转动90°并且开始失相之后,再聚焦脉冲,即180°脉冲的载波按照下述关系式相对于90°脉冲序列的载波发生相移:其中括弧中的表示式相对于n=1,2,...,N重复,其中N为在一个CPMG序列中收集的回波数量,回波间隔为表示使所说自旋围绕±x轴转动90°角的一个射频脉冲,如在磁共振测量(相位交替变换)的转动结构中所定义的一样。施加90°脉冲与180°脉冲之间的时间t0小于tcp,即回波间隔的一半。CPMG序列使得能够获得对称测量结果(即,不使用梯度线圈的测量结果)。精确的定时参数,t0、t1和t2依赖于各种因素(例如所施加脉冲的形状)。
在本发明中,施加到梯度线圈56的电流脉冲产生基本平行于所说静磁场的一个附加磁场。所说电流脉冲在所说第一90°和180°反相脉冲之间施加。这个附加场引起自旋的一个附加相移。由于180°反相脉冲对于附加相移没有补偿,所以受到附加场作用的自旋不形成自旋回波。但是,对于没有受到附加场作用的自旋,在时间2tcp产生一个自旋回波,而在每个反相脉冲之后在时间tcp连续产生低波幅的自旋回波。所说脉冲序列为
其中t0 a为90°脉冲与梯度脉冲持续时间δ之间的时间,t0 b为所说梯度脉冲与180°反相脉冲之间的时间,并且
由于随后的180°y脉冲和非均匀磁场的作用,NMR信号的x分量在几个回波时间内衰减。所以,我们仅仅关注信号的y分量。因此,忽略弛豫时间,第一NMR回波信号可以表示为:
其中i为复数虚部单位;γ为旋磁比;Mx 0和My 0为在没有梯度脉冲情况下在第一回波时间在位置r处的x和y磁化分量;G(r)为梯度磁场在同一位置平行于
的分量;δ为梯度脉冲的持续时间;而dc(r)表示NMR探测器的差动灵敏度。
所说梯度线圈56对于获取方位角测量结果具有许多优点。首先,因为使用轴对称天线探测自旋回波,可以在测井仪在井孔中转动的同时记录长回波序列。其次,由于线圈56无需象射频天线36、38那样进行调谐,所以线圈56使NMR-LMD测井仪的结构简化。第三,可以使用相同的天线36、38进行对称和轴对称测量。第四,可以使用所说线圈56获得具有极好的空间分辨率,特别是垂直分辨率的NMR测量结果。
本发明提供了获得方位角NMR测量数据的不同模式。例如,‘简单分解’模式使用至少一个线圈56在选定的象限内分解自旋,其中象限定义为测井仪10周围的一个角距离区段,但是,可以使用更多的线圈56分解多个象限。在任何一种情况下,都获得两种测量结果:具有固定等待时间的一个对称相位交变脉冲序列(PAPS),其后为具有可变等待时间的一个梯度PAPS,通过起动该象限中的线圈56分解选定象限。在一个优选实施例中,使用了上述梯度脉冲序列。从梯度测量结果中减去对称测量结果生成方位角测量结果。在这种模式下,对于每两个PAPS获得一个对称测量数据,对于每8个PAPS获得一次方位角扫描。因为将两种测量方式结合在一起,所以方位角测量的测量噪声大于对称或梯度测量中的噪声。
通过将不同的单个象限分解测量结果结合可以减少噪声影响。例如,通过分解各个象限可以获得四个梯度PAPS测量结果。将这些测量结果结合生成一个合成方位角和对称测量结果。通过将不使用所说梯度线圈56获得的测量结果与使用一个或多个梯度线圈56获得的测量结果相结合,可以生成地层构造的轴向或方位角分辨“图象”。所获得的渗透率和约束流体数据,特别是方位角图象形式的数据对于提高极其偏斜和水平井孔的岩石物理学分析水平和在钻探过程中作出地质井孔位置的确定是非常需要的。
脉冲长度和工作频率的最佳化
在选定的射频工作频率下,对于90°脉冲的持续时间t90,以及180°脉冲的持续时间t180,存在一个最佳持续时间,这个最佳时间可以确保所需的信噪比。在测井仪的主标定过程中可以寻找所说的最佳脉冲长度,从而可以正确地,或者当所说静磁场以一种没有预料到的方式变化时,例如由于在钻探过程中磁性碎屑的积累引起变化时,预置所有脉冲长度。参见转让给Schlumberger技术公司的美国专利申请No.09/031926(律师档案编号No.24786)。这种技术也可以用于选择满足其它条件,例如保持探查深度不变的条件的适合频率。
通过利用至少两个不同的脉冲持续时间和利用与构造的NMR特性无关的一种预定模式测量一种采样的NMR响应可以确定最佳脉冲长度。或者,利用至少两个不同的脉冲持续时间并结合利用从构造的NMR特性计算出来的一种模式可以确定所说最佳脉冲长度。在第一种情况下,积累数据可以提高信噪比,但是,积累过程可能需要很长的时间以获取构造数据。可取的是,在测井仪10暂停钻探操作时的静止时间窗口期间,例如在往钻杆柱上添加新的钻杆段期间积累测量数据。在第二种情况下,如果构造的T2分布是已知的,可以构成一种最好的数据采集模式,这种模式对于一个采集时间单元具有最大的信噪比,并且形成所获得回波的最佳线性组合。实验室模拟表明最佳采集数据模式的最佳时序是当回波序列的持续时间大约等于T2,max即构造的最大T2时间时,和当等待时间tw大约等于2.5×T2,max时(假定T1/T2为常数1.5)时实现的。最佳采集模式可以在百分之几秒到几秒时间的量级内确定最佳脉冲长度。可以采用类似的技术使NMR信号相对于频率最佳化(例如鞍点结构)。所说T2分布有效地帮助了对于测井仪10脉冲长度的有效调节。
数据采集模式
如上所述,测井仪10具有多个天线36、38。在本发明的一个优选实施例中,这些天线36,38并不同时发送或采集数据。可取的是,在一个天线36采集数据之后,另一个天线38在电源重新充电的同时经历最短的等待时间,以便传送下一个脉冲序列。本发明实现方案之一是同时传送或采集数据。此外,本发明还设想出无需等待时间的数据采集技术。
基于这些设计选择,可以采用多种数据采集模式。举例来说,下面记述了NMR数据采集的***性时序:适合于水湿砂岩地带的快速时序,适合于碳酸盐地带的慢速时序,和适合于碳氢化合物矿脉地带(或者油基泥浆渗入区域)的非常慢速时序。在表1中给出这些时序。
表1
等待时间(秒) | 回波间隔(秒) | 回波数量 | |
快速 | 2.3 | 0.5 | 400 |
慢速 | 4.6 | 0.5 | 800 |
非常慢速 | 9.2 | 1.0 | 800 |
对于每种数据采集时序可以使用若干种不同的模式,包括,但是不限于,下列模式:简单模式、交叉存取模式、和脉冲串模式。利用测井仪10获取T2信息的最简单方法是利用两个天线采用同一时序进行CPMG测量。图5表示运用表I所示的快速衰减、慢速衰减、和非常慢速衰减时序时所使用的简单数据采集模式。每个天线36、38交替获取一个长的脉冲序列,这样从每个天线36、38获得有效的渗透率测量数据。
运用交叉存取模式,所说高梯度天线在两个不同频率下测量至少两个圆柱体形壳层,而所说低梯度天线利用单一频率进行测量。图6表示使用表I所示时序中的快速衰减采样、慢速衰减成分、和非常慢速衰减成分时的交叉存取测量。
脉冲串模式提高了信噪比,特别是对于快速衰减成分。此外,脉冲串模式根据约束流体测量结果可以提供有用的T1值。参见转让给Numar公司的国际专利申请WO 98/29639(记载了用于确定纵向弛豫时间T1的一种方法)。还可参见转让给Schlumberger技术公司的美国专利申请No.09/096320(律师档案编号24785)(其中记载了用于极化构造中约束流体的一种方法)。图7表示利用略作改变的表1中时间对于快速衰减采样、慢速衰减成分、和非常慢速衰减成分的脉冲串测量。
除了简单模式、交叉存取模式、和脉冲串模式以外,利用本发明,通过探测井下状态,该操作在钻探作业中形成一个间歇,确定钻探模式,和利用该模式控制数据采集可以使构造评估测量最佳化。标准的转动钻探作业包含许多的自然停顿,在这个期间测井仪保持静止状态:在钻杆柱上添加新的钻杆段的连接时间,泥浆循环和钻杆管可能转动的循环时间,以及钻杆柱卡住和在恢复钻探之前必须拔出的探索或抖动时间。可以利用这些自然停顿进行NMR测量,所说自然停顿不会中断正常的钻探作业,或者故意暂停作业。所说钻探模式包括,但是不限于,钻孔、滑动、下钻、泵浆、探索、短起下钻(上下)和钻杆管连接。测定钻探模式提高了获取需要较长的时间或者要求安静的环境的NMR测量数据,例如,确定T1、T2、天线调谐和碳氢化合物类型的能力。参见转让给Schlumberger技术公司的美国专利申请No.09/031926(律师档案编号24786)。还可以根据环境的变化(例如冲蚀度、盐度等)和/或构造的NMR特性的变化(例如长T1与短T1的比值)调节采集数据模式。
利用在一个时间窗口上接收器电压的硬件积分获得自旋回波幅值。所说测井仪10利用相位敏感检测测量自旋回波信号加噪声幅值的同相和正交分量。可以使用在授予Robert Freedman的美国专利US-5381092中公开的技术在井下计算窗口总数,并将所说窗口总数传送到地面以进行T2转换处理和分析。此外,可以实施授予Abdurrahman Sezginer的美国专利US-5363041中公开的技术,利用一个线性算子映射弛豫时间相对于自旋回波的分布,生成线性算子的单值分解,确定SVD的矢量,和利用所说矢量压缩自旋回波数据。可取的是,在井下计算T2型谱,并传送到地面。这种技术具有消除在将计算T2型谱所需数据传送到地面过程中产生的遥测瓶颈效应的优点。可以使用一个数字信号处理器转换所说T2数据。所说自旋回波的幅值Aj由下列关系式表征:
其中ηj为测量值Aj中的噪声,ai为在T2,i时刻所取的T2分布幅值, 表示矩阵X的矩阵元,
其中tw为等待时间,c为一个常数(T1/T2的比值),T1是纵向弛豫时间,T2是自旋-自旋弛豫时间,Δt为回波间隔,而j=1,2,...N,其中N为在一个脉冲序列中收集的回波数量。用矩阵符号表示,该方程变为
因为噪声η是未知的,所以利用最小二乘法解,即函数 的最小值可以近似求得
这个方程的解受到数据中存在的噪声的严重影响,并且即使所说T2型谱不具有负的分量,这个解也可能具有负值分量。为了解决这个问题,在该函数中加入一调整项,
并在ai≥0,i=1...M约束条件下使用一种适合的迭代求最小值算法(例如共轭梯度投影算法)求出函数
的最小值。参见Ron S.Dembo和Ulrich Tulowitzski所写的文章“On the Minimization ofQuadratic Functions Subject to Box Constrains”,Yale Department ofComputer Science(1984年9月)(记载了共轭梯度投影方法)。利用数字信号处理器进行T2转换的所需时间是非常合理的。例如,假定在所说T2域中有1800个回波和30个采样,则数字信号处理器的转换只需要不到2秒的时间。
脉冲程序装置
在利用测井仪10进行基本NMR测量时,电子设备向探查构造施加一个脉冲序列。测井仪10包括一个脉冲程序装置80,该装置自适应地选择和控制施加到所说构造的脉冲序列。所说脉冲程序装置80利用在测量控制方框82(参见图8)中的信息和测井仪10的工作条件建立所说脉冲序列。可取的是,所说测量控制方框82保存在井下存储器中。方框82的结构是固定的,以使脉冲程序装置80容易在井下自发地适应和修改所说脉冲序列的时序。将方框82的一部分分为多个表84、86和88是有利的。所说的这些表84、86、88不是用于控制基于所说脉冲程序装置输出的脉冲序列的所有测井仪操作,而只是控制所说的这些操作。这使得所说脉冲程序装置80能够在不造成测井仪配置矛盾的前提下改变所说脉冲序列。所说的多个表84、86和88可以包括,但是不限于,下列各种表:描述堆栈缓存器布局的缓存器表、限定累积在缓存器中的获取信号的采集表、规定信号采集所使用的检测滤波器的滤波器系数表、指定各个缓存器所使用的自旋动态校正的自旋动态校正表、和指定从所获取的缓存数据计算出的核磁共振特征值的数据处理表。
所说脉冲程序装置80包括用于产生脉冲序列的一个脉冲序列模板94,它包括依赖于重复和时序变量的一系列状态。可以利用计算方框92由序列配置参数计算出这些变量。计算方框92可以一种可执行或可解释结构实施。根据将要测量的物理量例如T2的不同,可以将时间变量定义为等待时间tw、回波间隔techo、和所获取回波的数量。配置参数包括,但是不限于,t90、脉冲幅值、和脉冲形状。因为这些参数随着测井仪10工作条件的变化而变化,所以可以在测井仪10标定期间或者在测井仪10工作期间周期地计算这些参数。例如,所说脉冲幅值和形状依赖于天线品质因素,进而依赖于测井仪10周围构造的导电性。
通常,在脉冲程序装置80开始产生一个脉冲序列之后,所说序列按照确定方式运行,直至结束为止。为了利用测井仪10实施某些方位角测量模式,所说脉冲程序装置80具有在所说脉冲序列执行过程中改变脉冲序列的性能。程序装置80可以停止所说脉冲序列的执行,并进入一种HALT(暂停)状态,直到一个外加信号在时间tc结束了该状态或者直到经过一个最大时间长度tmax为止。如在本说明书数据采集模式一节中所述,由于可以与所说数据采集时序配合的不同模式中至少一种(交叉存取模式)需要交叉存取若干测量数据,所说程序装置80补偿了在暂停状态下经过的时间。可取的是,补偿是通过将暂停事件分组实现的。例如,一个分组可以包含一对暂停事件,其中一个暂停事件按照如上所述方式运行,另一个暂停事件是持续时间为tmax-tc的一个正常事件。将事件分组使得所说程序装置80可以将具有可变和确定时序的序列相结合。
此外,如在脉冲序列模板94中所定义的,所说状态序列对于序列的各个部分可以包括若干种不同的实施方式。在实践中,有一种实施方式(分支)是根据测井仪的外部条件(例如测井仪的方位)选择的。
为了说明和介绍本发明,上文中已经介绍了本发明的优选和替换实施例。它们并非是穷举或者将本发明限定为所公开的具体形式。显然,对于本领域技术人员来说可以作出许多改变和变化。选择并介绍这些实施例的目的是为了更好地解释本发明的原理及其实际应用,从而使本领域的其他技术人员能够理解本发明的各种实施例和适合于具体用途的各种改进。本发明的范围由所附的权利要求和它们的等价物所限定。
Claims (22)
1、用于测井的一种方法,该方法包括以下步骤:
a)在被井孔贯穿的构造中产生基本轴对称的静磁场;
b)在所说构造中产生一个振荡磁场;
c)探测从所说构造中产生的核磁共振信号;
d)在所说井孔中设置一个信号处理器;和
e)利用所说处理器,根据所探测的信号计算自旋-自旋弛豫时间分布。
2、如权利要求1所述的方法,其特征在于它还包括将自旋-自旋弛豫时间分布传输到所说井孔表面的步骤。
3、如权利要求1所述的方法,其特征在于它还包括向所说构造施加一个射频磁场脉冲序列,和利用自旋-自旋弛豫时间分布确定施加到所说构造中的每个脉冲的最佳时间长度的步骤。
4、如权利要求1所述的方法,其特征在于它还包括探测多个信号的步骤,所说信号的幅值具有信号加噪声的幅值Aj,其中
其中ηj为测量值Aj中的噪声,ai为在T2,i时刻所取的自旋-自旋弛豫时间(T2)分布幅值,
其中tw为等待时间,c为T1/T2的比值,c是一个常数,T1是纵向弛豫时间,T2是自旋-自旋弛豫时间,Δt为回波间隔,而j=1,2,...N,其中N为在一个脉冲序列中收集的回波数量。
5、如权利要求4所述的方法,其特征在于所说信号加噪声幅值 并且还包括在
的约束条件下近似求得
的步骤。
6、如权利要求5所述的方法,其特征在于它还包括确定函数 的最小值的步骤。
7、如权利要求5所述的方法,其特征在于它还包括选择参数,调整项λ,和确定函数
的最小值的步骤。
8、如权利要求6所述的方法,其特征在于它还包括利用一种共扼梯度投影算法求得所说函数最小值的步骤。
9、如权利要求7所述的方法,其特征在于它还包括利用一种共扼梯度投影算法求得所说函数最小值的步骤。
10、如权利要求1所述的方法,其特征在于它还包括向所说构造施加一个射频磁场脉冲序列,和利用所说自旋-自旋弛豫时间分布选择一个最佳工作频率的步骤。
11、如权利要求1所述的方法,其特征在于它还包括向所说构造施加一个射频磁场脉冲序列,和利用所说自旋-自旋弛豫时间分布保持所说构造中基本恒定深度的探查的步骤。
12、用于确定一个井孔周围地层构造的探查区域中核磁共振特性的一种设备,该设备包括:
a)用于在被一个井孔贯穿的构造中产生一个基本轴对称静磁场的装置;
b)用于在所说构造中产生一个振荡磁场的装置;
c)用于探测所说构造中的核磁共振信号的装置;和
d)设置在所说井孔中、用于根据所探测的信号计算自旋-自旋弛豫时间分布的装置。
13、如权利要求12所述的装置,其特征在于它还包括用于将所说自旋-自旋弛豫时间传输到所说井孔表面的装置。
14、如权利要求12所述的装置,其特征在于它还包括用于向所说构造施加一个磁场脉冲序列的装置,和利用所说自旋-自旋弛豫时间分布确定施加到所说构造中的每个脉冲的最佳时间长度的装置。
15、如权利要求12所述的装置,其特征在于它还包括用于向所说构造施加一个磁场脉冲序列的装置,和利用所说自旋-自旋弛豫时间分布选择一个最加工作频率的装置。
16、如权利要求15所述的装置,其特征在于它还包括用于向所说构造施加一个磁场脉冲序列的装置,和利用所说自旋-自旋弛豫时间分布保持所说构造中探查深度基本恒定的装置。
17、如权利要求12所述的装置,其特征在于它还包括用于探测多个信号的装置,所说信号的幅值具有信号加噪声的幅值Aj,
其中
其中ηj为信号中的噪声,ai为在T2,i时刻所取的T2分布幅值, 表示矩阵X的矩阵元,
其中tw为等待时间,c为T1/T2的比值,c是一个常数,T1是纵向弛豫时间,T2是自旋-自旋弛豫时间,Δt为回波间隔,而j=1,2,...N,其中N为在一个脉冲序列中收集的回波数量。
18、如权利要求17所述的装置,其特征在于信号加噪声幅值 并且还包括在
的约束条件下近似求得
的装置。
19、如权利要求18所述的装置,其特征在于它还包括用于确定函数 的最小值的装置。
20、如权利要求18所述的装置,其特征在于它还包括用于选择参数,调整项λ,和确定函数
的最小值的装置。
21、如权利要求19所述的装置,其特征在于它还包括利用一种共扼梯度投影算法求得所说函数最小值的装置。
22、如权利要求20所述的装置,其特征在于它还包括利用一种共扼梯度投影算法求得所说函数最小值的装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/187,130 US6184681B1 (en) | 1998-03-03 | 1998-11-05 | Apparatus and method for computing a distribution of spin-spin relaxation times |
US09/187130 | 1998-11-05 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1253294A CN1253294A (zh) | 2000-05-17 |
CN1186649C true CN1186649C (zh) | 2005-01-26 |
Family
ID=22687726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB99123507XA Expired - Fee Related CN1186649C (zh) | 1998-11-05 | 1999-11-05 | 用于计算自旋-自旋弛豫时间分布的设备和方法 |
Country Status (11)
Country | Link |
---|---|
US (1) | US6184681B1 (zh) |
CN (1) | CN1186649C (zh) |
AR (1) | AR021589A1 (zh) |
AU (1) | AU735770B2 (zh) |
CA (1) | CA2286240C (zh) |
DE (1) | DE19952787A1 (zh) |
EG (1) | EG22135A (zh) |
GB (1) | GB2343959B (zh) |
ID (1) | ID23759A (zh) |
NL (1) | NL1013482C2 (zh) |
NO (1) | NO327711B1 (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6255817B1 (en) * | 1997-06-23 | 2001-07-03 | Schlumberger Technology Corporation | Nuclear magnetic resonance logging with azimuthal resolution |
US6891369B2 (en) * | 1998-08-13 | 2005-05-10 | Schlumberger Technology Corporation | Nuclear magnetic resonance method and logging apparatus for fluid analysis |
US6445180B1 (en) * | 2000-06-28 | 2002-09-03 | Baker Hughes Incorporated | Nuclear magnetic resonance tool with active RF spoiler antenna |
WO2002008795A2 (en) * | 2000-07-21 | 2002-01-31 | Services Petroliers Schlumberger | Nuclear magnetic resonance measurements and methods of analyzing nuclear magnetic resonance data |
WO2003010963A1 (en) * | 2001-07-11 | 2003-02-06 | Daniel Kilbank | System and method for compressing data |
JP2005524260A (ja) * | 2002-04-23 | 2005-08-11 | キルバンク ダニエル | 通信においてマイクロレッツを使用するためのシステムおよび方法 |
US6838875B2 (en) * | 2002-05-10 | 2005-01-04 | Schlumberger Technology Corporation | Processing NMR data in the presence of coherent ringing |
WO2004025308A2 (en) * | 2002-09-11 | 2004-03-25 | Halliburton Energy Services, Inc. | Nmr tool with helical polarization |
US7408486B2 (en) * | 2003-04-21 | 2008-08-05 | Qbit Corporation | System and method for using a microlet-based modem |
US7463027B2 (en) * | 2003-05-02 | 2008-12-09 | Halliburton Energy Services, Inc. | Systems and methods for deep-looking NMR logging |
US20050083061A1 (en) * | 2003-10-17 | 2005-04-21 | Tabanou Jacques R. | Methods and systems for estimating formation resistivity that are less sensitive to skin effects, shoulder-bed effects and formation dips |
CN100373172C (zh) * | 2004-10-21 | 2008-03-05 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
US20080036457A1 (en) * | 2005-03-18 | 2008-02-14 | Baker Hughes Incorporated | NMR Echo Train Compression |
CN100349013C (zh) * | 2005-05-27 | 2007-11-14 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
US8115481B2 (en) * | 2007-07-23 | 2012-02-14 | Baker Hughes Incorporated | Method of discerning water from hydrocarbon fluids using downhole NMR instruments in petroleum reservoirs with formation brine |
US8431497B2 (en) | 2009-08-25 | 2013-04-30 | Berkshire Corporation | Clean room wipes |
JP5481163B2 (ja) * | 2009-10-30 | 2014-04-23 | 株式会社東芝 | 磁気共鳴映像装置 |
US8836328B2 (en) * | 2010-02-03 | 2014-09-16 | Baker Hughes Incorporated | Acoustic excitation with NMR pulse |
US9201159B2 (en) | 2010-06-16 | 2015-12-01 | Halliburton Energy Services, Inc. | Nuclear magnetic resonance logging tool having an array of antennas |
FI123993B (fi) * | 2010-08-31 | 2014-01-31 | Metso Automation Oy | Menetelmä ja laite pulssivälin muuttamiseksi adaptiivisesti NMR-pohjaisessa vesipitoisuusmittauksessa |
MX2014009202A (es) | 2012-01-30 | 2014-10-14 | Schlumberger Technology Bv | Metodo para llevar a cabo correccion de errores de datos rmn. |
US9213124B2 (en) * | 2013-03-22 | 2015-12-15 | Oliden Technology, Llc | Restorable antennae apparatus and system for well logging |
BR112018003876A2 (pt) * | 2015-10-06 | 2018-09-25 | Halliburton Energy Services Inc | ?ferramenta de perfilagem de ressonância magnética nuclear, e, método? |
US11091997B2 (en) * | 2015-11-11 | 2021-08-17 | Schlumberger Technology Corporation | Estimating nuclear magnetic resonance measurement quality |
CN105298486B (zh) * | 2015-11-25 | 2019-04-12 | 中国石油集团钻井工程技术研究院 | 井下可控的随钻核磁共振测井装置 |
CN108590629B (zh) * | 2018-03-23 | 2020-09-18 | 中国石油大学(北京) | 井下三维扫描核磁共振成像测井仪探头、天线、及仪器 |
CN110850491B (zh) * | 2018-08-20 | 2021-11-02 | 中国石油天然气股份有限公司 | T2谱反演方法、装置和存储介质 |
CN113014256B (zh) * | 2021-02-23 | 2023-10-13 | 北京无线电计量测试研究所 | 一种腔耦合原子***制备自旋压缩态的方法 |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4710713A (en) | 1986-03-11 | 1987-12-01 | Numar Corporation | Nuclear magnetic resonance sensing apparatus and techniques |
US5055788A (en) | 1986-08-27 | 1991-10-08 | Schlumberger Technology Corporation | Borehole measurement of NMR characteristics of earth formations |
US4717878A (en) | 1986-09-26 | 1988-01-05 | Numar Corporation | Nuclear magnetic resonance sensing apparatus and techniques |
US4949045A (en) | 1987-10-30 | 1990-08-14 | Schlumberger Technology Corporation | Well logging apparatus having a cylindrical housing with antennas formed in recesses and covered with a waterproof rubber layer |
GB9022145D0 (en) | 1990-10-11 | 1990-11-21 | Oxford Instr Ltd | Magnetic field generating assembly |
US5212447A (en) | 1990-12-03 | 1993-05-18 | Numar Corporation | Apparatus and technique for nmr diffusion measurement |
DE69123260T2 (de) | 1990-12-05 | 1997-06-12 | Numar Corp., Malvern, Pa. | Anordnung zum nmr-messen eines bohrlochs während dessen bohrens |
US5705927A (en) | 1992-07-30 | 1998-01-06 | Schlumberger Technology Corporation | Pulsed nuclear magnetism tool for formation evaluation while drilling including a shortened or truncated CPMG sequence |
EP0581666B1 (en) | 1992-07-30 | 1997-10-01 | Schlumberger Limited | Pulsed nuclear magnetism tool for formation evaluation while drilling |
US5291137A (en) | 1992-11-02 | 1994-03-01 | Schlumberger Technology Corporation | Processing method and apparatus for processing spin echo in-phase and quadrature amplitudes from a pulsed nuclear magnetism tool and producing new output data to be recorded on an output record |
US5486762A (en) | 1992-11-02 | 1996-01-23 | Schlumberger Technology Corp. | Apparatus including multi-wait time pulsed NMR logging method for determining accurate T2-distributions and accurate T1/T2 ratios and generating a more accurate output record using the updated T2-distributions and T1/T2 ratios |
US5596274A (en) | 1992-12-31 | 1997-01-21 | Schlumberger Technology Corporation | Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences |
US5363041A (en) | 1992-12-31 | 1994-11-08 | Schlumberger Technology Corporation | Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences |
US5497087A (en) * | 1994-10-20 | 1996-03-05 | Shell Oil Company | NMR logging of natural gas reservoirs |
AU711508B2 (en) | 1995-03-23 | 1999-10-14 | Schlumberger Technology B.V. | Nuclear magnetic resonance borehole logging apparatus and method |
US5696448A (en) | 1995-06-26 | 1997-12-09 | Numar Corporation | NMR system and method for formation evaluation using diffusion and relaxation log measurements |
US5757186A (en) | 1996-02-23 | 1998-05-26 | Western Atlas International, Inc. | Nuclear magnetic resonance well logging apparatus and method adapted for measurement-while-drilling |
US5828214A (en) | 1996-02-23 | 1998-10-27 | Western Atlas International, Inc. | Method and apparatus for resistivity determination by nuclear magnetic resonance measurement |
MY122012A (en) * | 1996-03-14 | 2006-03-31 | Shell Int Research | Determining a fluid fraction in an earth formation |
US5710511A (en) | 1996-03-14 | 1998-01-20 | Western Atlas International, Inc. | Method and apparatus for eddy current suppression |
US6051973A (en) | 1996-12-30 | 2000-04-18 | Numar Corporation | Method for formation evaluation while drilling |
US6147489A (en) * | 1997-04-09 | 2000-11-14 | Schlumberger Technology Corporation | Method and apparatus for measuring total nuclear magnetic resonance porosity |
US5977768A (en) * | 1997-06-23 | 1999-11-02 | Schlumberger Technology Corporation | Nuclear magnetic resonance logging with azimuthal resolution |
US6049205A (en) | 1997-10-01 | 2000-04-11 | Western Atlas International, Inc. | Multiple frequency method for nuclear magnetic resonance longitudinal relaxation measurement and pulsing sequence for power use optimization |
-
1998
- 1998-11-05 US US09/187,130 patent/US6184681B1/en not_active Expired - Lifetime
-
1999
- 1999-10-06 GB GB9923667A patent/GB2343959B/en not_active Expired - Fee Related
- 1999-10-11 AU AU53588/99A patent/AU735770B2/en not_active Ceased
- 1999-10-13 CA CA002286240A patent/CA2286240C/en not_active Expired - Fee Related
- 1999-10-25 ID IDP990981A patent/ID23759A/id unknown
- 1999-11-02 EG EG137599A patent/EG22135A/xx active
- 1999-11-03 NL NL1013482A patent/NL1013482C2/nl not_active IP Right Cessation
- 1999-11-03 AR ARP990105561A patent/AR021589A1/es not_active Application Discontinuation
- 1999-11-03 DE DE19952787A patent/DE19952787A1/de not_active Withdrawn
- 1999-11-04 NO NO19995393A patent/NO327711B1/no not_active IP Right Cessation
- 1999-11-05 CN CNB99123507XA patent/CN1186649C/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
EG22135A (en) | 2002-08-30 |
CA2286240A1 (en) | 2000-05-05 |
CN1253294A (zh) | 2000-05-17 |
AU735770B2 (en) | 2001-07-12 |
NL1013482C2 (nl) | 2003-03-06 |
AR021589A1 (es) | 2002-07-31 |
ID23759A (id) | 2000-05-11 |
NO995393L (no) | 2000-05-08 |
NL1013482A1 (nl) | 2000-05-09 |
NO327711B1 (no) | 2009-09-14 |
GB9923667D0 (en) | 1999-12-08 |
DE19952787A1 (de) | 2000-05-11 |
CA2286240C (en) | 2002-01-29 |
GB2343959B (en) | 2001-01-17 |
NO995393D0 (no) | 1999-11-04 |
US6184681B1 (en) | 2001-02-06 |
GB2343959A (en) | 2000-05-24 |
AU5358899A (en) | 2000-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1186649C (zh) | 用于计算自旋-自旋弛豫时间分布的设备和方法 | |
CN1253296A (zh) | 用于在钻探过程中获取核磁共振测量数据的设备和方法 | |
CN1253295A (zh) | 用于产生脉冲序列的设备和方法 | |
US6580273B2 (en) | Side-looking NMR probe for oil well logging | |
US7339374B2 (en) | Correction of NMR artifacts due to axial motion and spin-lattice relaxation | |
US9097818B2 (en) | Kerogen porosity volume and pore size distribution using NMR | |
US6018243A (en) | NMR well logging apparatus and method | |
US20090072825A1 (en) | Systems and methods for deep-looking nmr logging | |
CA2415824C (en) | High-resolution high-speed nmr well logging device | |
US20010043066A1 (en) | Nuclear magnetic resonance measurements in well logging using motion triggered pulsing | |
NO345909B1 (en) | Downhole nuclear magnetic resonance (NMR) tool with transversal-dipole antenna configuration | |
CN1249436A (zh) | 使用梯度线圈进行的具有角向分辨率的核磁共振测井 | |
CA2506752A1 (en) | Azimuthal nmr imaging of formation properties from a wellbore | |
NO20231225A1 (en) | Reducing motion effects on nuclear magnetic resonance relaxation data | |
Taicher et al. | A comprehensive approach to studies of porous media (rocks) using a laboratory spectrometer and logging tool with similar operating characteristics | |
CA2476499A1 (en) | Side-looking nmr probe for oil well logging | |
MXPA99009936A (en) | Apparatus and method for calculating a spin-es relaxation time distribution | |
MXPA99009935A (en) | Apparatus and method for obtaining a measurement of nuclear magnetic resonance during perforation | |
MXPA99009828A (en) | Apparatus and method for generating a pu sequence |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050126 Termination date: 20171105 |
|
CF01 | Termination of patent right due to non-payment of annual fee |