使用可变形状标准的核频谱分析中环境校正的制作方法

使用可变形状标准的核频谱分析中环境校正的制作方法
【专利摘要】地层情况可以通过将包含中子源和伽马射线探测器的核频谱分析测井仪放置到井眼中且执行多个环境测量来进行估测。从所述核频谱分析测井仪中发射中子，从而使得所述中子中的一些从与所述核频谱分析测井仪相邻的地层中产生伽马射线，所述中子中的一些从所述核频谱分析测井仪之内的元素中产生伽马射线，且所述中子中的一些从钻井泥浆中的元素中产生伽马射线。可以采用所述测井仪探测由发射的中子诱导的伽马射线的能谱。探测的伽马射线能谱可以使用标准能谱的组合进行分析，其中可以基于所述环境测量而改变至少一个标准能谱的形状以确定环境对伽马射线能谱的影响。
【专利说明】使用可变形状标准的核频谱分析中环境校正
[0001]相关申请
[0002]本申请要求2011年5月24日提交的系列号为61/489，295，名称为“SYSTEMS ANDMETHODS FOR ENVIRONMENTAL CORRECTIONS IN NUCLEAR SPECTROSCOPY USING VARIABLESHAPE STANDARD”的美国临时专利申请的优先权，将该申请的全部内容通过援引加入的方式纳入本文。
【技术领域】
[0003]本公开宽泛地说涉及核频谱分析(nuclear spectroscopy),更具体而言涉及估测和确定(account for)环境对核频谱分析结果的影响的方法。
【背景技术】
[0004]当为了获得石油、天然气或其它材料而钻井通过岩石地层或其它地层时，确定或估测正被钻井通过的地层的类型是有益的。例如，与地层类型有关的信息可用于更精确地解释其它测井测量，从而估测地层孔隙度、含水饱和度、净含烃量、以及地层渗透性和产量。其也可用于基于估测的地层机械特性而作出的钻井决定中。
[0005]存在各种尝试确定或估测地层的方法，其中的一个是随钻测井或LWD。LWD是在挖掘井孔期间、或此后不久通过使用集成到井底钻具组件中的测井仪测量地层的特性。
[0006]LffD的一个已知技术是中子-伽马频谱分析，其使用中子在原子核中产生激发态，然后所述中子可以通过发射一个或多个伽马射线进行衰减。每一个同位素的伽马射线具有可以用于区分特定元素(例如氢)的浓度的特征能谱。结果是，一组同位素的伽马射线提供一组可用于区分相应元素一例如氢、氯、铁、硅、钙、硫、钛、铝、钠、镁、锰和镍一的浓度的特征能谱。
[0007]使用中子-伽马频谱分析的一个缺点是在频谱分析测井仪(spectroscopy tool)中的一个或多个元素可能与地层中的一个或多个元素相同，因而当尝试确定或估测地层中的元素时，所述测井仪可能实际还探测到由地层以外的元素——例如测井仪或钻井泥浆中的元素——产生的伽马射线。就是说，存在于地层和例如测井仪中的元素的能谱实际上是来自地层的元素的能谱与来自测井仪的元素的能谱的总和。来自地层以外的能谱部分通常称为背景能谱。
[0008]另一个缺点是地层的环境条件可能因井而异且在特定井之内还可随深度而变，这可能影响使用LWD中子-伽马频谱分析测量的能谱。不同的地层环境参数的实例是井眼尺寸、岩石地层密度、井眼流体密度、岩石地层的含氢指数、岩石地层的中子慢化长度、岩石地层的热中子俘获截面西格玛、或井眼流体的热中子俘获截面西格玛。
[0009]当尝试确定或估测被钻井通过的地层的类型时，消除、减少这些缺点中任一个的影响、或补偿这些缺点中的任一个都将是有益的。

【发明内容】
[0010]本公开描述了用于确定环境对核频谱分析测量的影响的几个不同的方法、器械、以及系统。
[0011]根据本公开的一个实施方案，使用核频谱分析测井仪估测地层情况的方法包括将包含中子源和伽马射线探测器的核频谱分析测井仪(nuclear spectroscopy tool)放置在井眼中以及执行多个环境测量。从所述核频谱分析测井仪中发射中子，从而使得所述中子中的一些从与所述核频谱分析测井仪相邻的地层中产生伽马射线，所述中子中的一些从所述核频谱分析测井仪之内的元素中产生伽马射线，且所述中子中的一些从钻井泥浆中的元素中产生伽马射线。可以采用测井仪探测由发射的中子诱导的伽马射线的能谱。探测的伽马射线能谱可以使用标准能谱的组合进行分析，其中可以基于所述环境测量而改变至少一个标准能谱的形状以确定环境对伽马射线能谱的影响。
[0012]根据本公开的一个实施方案，用于估测地层情况的核频谱分析测井仪包含这样一种中子源，其将中子发射到与所述测井仪相邻的地层中，从而使得所述中子中的一些可以从岩石地层中产生伽马射线，所述中子中的一些可以从所述核频谱分析测井仪中的元素中产生伽马射线，且所述中子中的一些可以从钻井泥浆中的元素中产生伽马射线。所述测井仪包含伽马射线探测器，其探测由发射的中子诱导的伽马射线的能谱；以及数据处理电路，其使用标准能谱的组合对探测的伽马射线能谱执行分析，其中基于测井仪和地层的环境测量而改变至少一个标准能谱的形状，以确定环境对伽马射线能谱的影响，并且使用所述分析估测所述地层的情况。
[0013]根据本公开的一个实施方案，用于估测地层情况的系统包含这样一种中子源，其将中子发射到与所述测井仪相邻的地层中，从而使得所述中子中的一些可以从该地层中产生伽马射线，所述中子中的一些可以从所述核频谱分析测井仪中的元素中产生伽马射线，且所述中子中的一些可以从钻井泥浆中的元素中产生伽马射线。所述系统包含伽马射线探测器，其探测由发射的中子诱导的伽马射线的能谱；测量装置，其用于执行所述测井仪的环境的测量；以及数据处理电路，其使用标准能谱的组合对探测的伽马射线能谱执行分析，其中基于所述环境测量而改变至少一个标准能谱的形状，用以确定环境对伽马射线能谱的影响，并且使用所述分析估测所述地层的情况。
[0014]其它实施方案包括实现上文描述的方法学实施方案的结构、装置、器械和系统。另夕卜，尽管公开了具有多个元素或方面的多个实施方案，但是对于本领域技术人员来说，从示出且描述了本发明的示例性实施方案的下文详细的说明书中将明显可以得出本发明的其它实施方案、元素、以及方面。因此，附图和详细的说明书将被视为实质上是说明性的而非限制性的。
【专利附图】

【附图说明】
[0015]图1是根据本公开的实施方案的井场系统的示意图。
[0016]图2是根据本公开的实施方案的随钻测井频谱分析测井仪的示意性图例。
[0017]图3和4是根据本公开的实施方案的示意性框图。
[0018]图5是说明了测井仪背景铁比率与地层西格玛的相关性的图表。
[0019]图6是说明了测井仪背景铁比率与地层慢化长度的相关性的图表。
[0020]图7是说明了井眼尺寸如何缓和测井仪背景铁比率与地层西格玛的相关性的图表。
[0021]图8是说明了地层基质铁重量浓度差异的图表。
[0022]图9是说明了背景铁比率的参数化的图表。
[0023]图1OA和IOB是说明了地层基质铁重量浓度的柱状图。
[0024]图1lA和IlB是说明了相对于地层西格玛的参数化的测井数据中铁比率的图表。
[0025]图12是说明了根据本公开的实施方案的方法的流程图。
[0026]图13是说明了测井仪内部背景能谱和测井仪外部背景能谱的实验取得的实例的图表。
[0027]图14是说明了通过建模取得的内部和外部测井仪背景能谱的图表。
[0028]图15是说明了根据本公开的实施方案的方法的流程图。
[0029]图16是说明了来自起源于地层中的伽马射线和来自起源于井眼或泥浆通道中的伽马射线的建模能谱的图表。
[0030]图17是说明了中子捕获的空间分布的蒙特卡罗(Monte Carlo)建模的图表。
[0031]图18是说明了根据本公开的实施方案的方法的流程图。
【具体实施方式】
[0032]包括方法、器械和系统实施方案在内的本公开的一个或多个具体实施方案将在下文中描述。这些描述的实施方案和它们的各种不同元素仅仅是目前公开的技术的实例。公开的器械的结构方面提供装置用于实现公开的方法的不同方面。应该理解，在任何这样实际实施方式的开发中，例如在任何工程或设计项目中，可能作出大量关于实施方式的具体决定用以实现开发者的特定目标，例如遵从系统相关和商业相关的限制，所述目标可能因实施方式而不同。而且，应该理解，这样的开发努力可能是费时的，但对于受益于本公开的本领域技术人员而言将仍然是设计、制造和生产的常规任务。
[0033]当引入本公开的各种不同实施方案的元素时，表述“一”、“一个”以及“该”旨在表示存在一个或多个所述元素。术语“包含”、“包括”、以及“具有”旨在表示包括且表示可能存在所列元素以外的附加元素。此外，应该理解，提及本公开的“一个实施方案”或“实施方案”不旨在解释成排除也包含所列元素的附加实施方案的存在。
[0034]图1说明了井场测量器械、系统和方法的实施方案。图1的井场系统能够在岸上或在海上用以，例如勘探和生产石油、天然气、以及其它资源，它们能够被使用、炼制、以及其它方式加工用于燃料、原料、以及其它目的。在图1的井场系统中，井眼11能够在地下地层(例如岩石地层)中通过使用任意合适技术的旋转钻井来形成。可将钻柱12悬挂在井眼11之内且能够在其下端具有包含钻头105的井底钻具组合100。图1的井场系统的地面系统能够包含布置在井眼11之上的平台和钻塔组件10，其中所述平台和钻塔组件10包含转盘16、方钻杆17、吊钩18以及转环(rotary swivel) 19。钻柱12能够通过转盘16旋转，所述转盘16通过任意合适的装置提供动力且在钻柱12的上端啮合方钻杆17。钻柱12可从吊钩18处悬挂下来通过方钻杆17和转环19，其中所述吊钩18附连到游动滑车(未示出)且所述转环19使得钻柱12可相对于吊钩18旋转。可以替换性地使用顶部驱动系统，其可以是本领域技术人员熟知的顶部驱动系统。
[0035]在图1的井场系统中，所述地面系统也可以包含在井场处形成的存储在池27中的钻井液或泥浆26。泵29可以通过转环19中的端口将钻井液26输送到钻柱12的内部，从而使钻井液如方向箭头8所示向下流动通过钻柱12。钻井液26可以通过钻头105中的端口离开钻柱12，且如方向箭头9所示向上环行通过钻柱12的外部与井眼11的井壁之间的环状区域。以这种方式，钻井液26润滑钻头105且当钻井液26返回到池27用于再循环时将地层钻屑向上携带到地面。
[0036]作为一个例子，图1的井场系统的井底钻具组合100可以包含随钻测井(LWD)模块120、另一个类型的随钻测量(MWD)模块130、旋转式可转向系统(roto-steerable system)和马达150、以及钻头105中的一个或多个。LWD模块120可以安装在特殊类型的钻井钻铤中。还将理解的是，可以使用一个以上的LWD模块或LWD模块之内的测井工具(其通常由附图标记120A表示)。因此，提及LWD模块120还可以替换性地意味着也适用于在120A的位置处的模块。LWD模块120可以具有用于测量、处理和存储信息以及用于与地面设备通信的能力。
[0037]图2呈现了为能谱测井工具121的或包含能谱测井工具121的LWD模块的示例性实施方案的侧视图的示意性框图。所述测井工具121包含底盘122、钻铤123以及延伸通过所述测井工具121的流动导管124。中子源125位于测井工具121中第一位置处，且探测器126 (例如伽马射线探测器)位于与所述中子源125轴向间隔开的第二位置处。中子屏蔽127 (例如硼屏蔽)相对于中子源125径向布置在第二位置处或接近第二位置处。设想到关于这个实施方案以及使用LWD模块120和测井工具122的通常配置或方面的能谱测井仪的其它实施方案的细节与任意合适的运输装置一起使用，所述运输装置例如电缆、挠性管道、随钻测井(LWD)，等等。进一步地，关于环境的信息(例如地层西格玛、泥浆西格玛、密度、井眼尺寸、以及慢化长度)可以使用如下文进一步讨论的额外的设备来获得。
[0038]图3和4呈现LWD模块120的实施方案的侧视图的示意性框图，其中图3和图4彼此之间旋转90度。LWD模块120可以与数据处理电路200连接。尽管LWD模块120和数据处理电路200作为图3中的独立元素描述，但是可以理解的是，数据处理电路200可以全部在LWD模块120之内、在远离LWD模块120的地面处、或部分在LWD模块120之内且部分在地面处实施。通过举例，LWD模块120可以表不Schlumberger Technology Corporation的EcoScope?工具模型。
[0039]LffD模块120的这个实施方案的底盘204可以包括用于发射且探测辐射用以获得频谱分析测量的多种元件和构造。例如，中子源206可以作为发射至少2MeV中子的中子源，用以通过与地层元素的非弹性散射产生伽马射线。通过举例，中子源206可以是电子中子源，例如 Minitron? 装置(其可商购于 Schlumberger Technology Corporation),其可以通过氣核-氣核(d-D)、氣核-氣核(d-T)、氣核-氣核(t-T )或其它合适的反应产生中子脉冲。因此，中子源206可以发射例如大约2MeV或14MeV的中子，或者发射具有宽范围能量的中子，例如由连续的中子源(如241AmBe或252Cf放射性同位素源)所获得的那些。
[0040]在一些实施方案中，频谱分析测井仪可以包含中子监视器208，其可以监视中子从中子源206的发射。通过举例，中子源208可以是塑料闪烁器和光电倍增管，其主要探测直接从中子源206发射的非散射中子，且因此可以提供与来自中子源206的中子输出速率的中子输出率成比例的计数率信号。合适的中子监视器的示例性但非限制性的实例在美国专利号6，884，994中描述，此专利通过援引加入的方式纳入本文。[0041]中子屏蔽210 (其可以包含例如钨、铅或硼)可以提供装置用于大规模地防止来自中子产生器206的中子在内部通过LWD模块120朝向所述屏蔽210的另一侧上的各种不同的福射-探测元件传递。合适的鹤屏蔽材料可商购于115ConstitutionBoulevard, Franklin MA020038的PLANSEE USA LLC。合适的硼屏蔽可以由多种源获得且可以包含几种不同形式的硼(例如金属性硼(metallic boron)、B4C> BN以及其它形式的硼)。在一些实施方案中，使用由kiB同位素富集的硼且其可商购于Ceradyne，P.0.Box798,Quapaw 0Κ74363。
[0042]如图3和4中所示，LWD模块120可以包含两个接近或并置的中子探测器，例如热中子探测器212以及超热中子探测器214。还可以包含两个远离的热中子探测器216Α和216Β且其可以位于比中子探测器212和214距中子产生器206的间隔更远处。例如，接近的中子探测器212和214可以与中子产生器206大约间隔10-14英寸，且远离的中子探测器216Α和216Β可以与中子产生器206间隔18-28英寸。中子探测器可商购于TwinsburgOH 的 GE Reuter Stokes 以及 Houston, TX 的 Schlumberger Technology Corporation。
[0043]短间隔(SS)伽马射线探测器218 (例如使用Nal、LaBr、或GSO的探测器)可以位于接近的中子探测器212和214与远离的中子探测器216A和216B之间。长间隔(LS)伽马射线探测器220可以位于超过远离的中子探测器216A和216B处，位于比伽马射线探测器218距中子产生器206更远处。例如，SS伽马射线探测器218可以与中子产生器206大约间隔10-22英寸，且LS伽马射线探测器220可以与中子产生器206大约间隔24-38英寸。伽马射线探测器可商购于 17900Great Lakes Parkway, Hiram 0H44234-9681 的 Saint-GobainCrystals。LWD模块120的可选实施方案可包含更多或更少的上述的辐射探测器，但是通常可以包含至少一个伽马射线探测器。中子探测器212、214、216A、和/或216B可以是任意合适的中子探测器，例如3He中子探测器。为了探测主要的超热中子，超热中子探测器214可以被热中子屏蔽围绕，而热中子探测器212、216A、和/或216B不被热中子屏蔽围绕。
[0044]伽马射线探测器218和/或220可以是被中子屏蔽围绕的闪烁探测器。例如，中子屏蔽可以包括6Li (例如碳酸锂(Li2CO3)),其可以大体上屏蔽伽马射线探测器218和/或220免于探测热中子且没有产生热中子捕获伽马射线。当来自中子产生器206的快中子非弹性地散射掉周围地层的特定元素时，伽马射线探测器218和220可以探测产生的非弹性伽马射线。
[0045]来自伽马射线探测器218和220的伽马射线的计数率和能谱以及来自中子探测器212、214、216A、和/或216B的中子的计数率可以由数据处理电路200作为数据222来接收。可以是其它提到的元件或结构或者单独的元件或结构的一部分的数据处理电路200提供装置用于接收数据222且执行特定步骤或处理用以确定或估测周围地层的一个或多个特性，例如地层矿物学和本文公开的其它特性。数据处理电路200可以包含处理器224、算法(用于执行步骤用以实现所公开的目标的步骤或指令；未示出)、记忆装置(mem0ry)226、和/或存储器(storage) 228。处理器224可以可操作地耦合到记忆装置226和/或存储器228且用于执行本文中所公开的步骤或算法。本文公开的技术可以由处理器224和/或其它数据处理电路基于可由处理器224执行的相应的指令来执行。上述的指令可以使用任何合适的制品来存储，所述制品可以包含一个或多个有形的、计算机可读的介质用以至少集中(collectively)存储这些指令。所述制品可以包含例如记忆装置226和/或非易失存储器228。记忆装置226和非易失存储器228可以包含任何合适的制品用于存储数据和可执行指令，例如随机存取记忆装置、只读记忆装置、可重写闪存、硬盘、以及光盘。记忆装置226和存储器228在本文中总体地称为数据存储器。处理器的一个实例是在美国专利5，067，090 中描述的 Digital Equipment Corporation F1DP-1l 计算机之内的处理器。其它的处理器可以是在通用计算机中使用的那些，其使用Windows XP、Windows Vista、以及Windows7尤其用于在地面处处理。适合用在随钻测井仪(drilling tool)中的处理器的具体实例有32位浮点处理器，例如Texas Instruments TMS320VC33DSP。
[0046]LffD模块120可以将数据222传送到数据处理电路200，例如通过在测井仪之内的内部连接、遥感勘测系统通信上行线路、和/或通信电缆。数据处理电路200可以确定或估测周围地层的一个或多个特性。通过举例，上述特性可以包含元素的捕获伽马射线的相对能谱产额或地层的元素浓度。此后，数据处理电路200可以输出表明地层的频谱分析测量的信息或报告230。报告230可以存储在记忆装置或存储器中用于后期由电路200或由其它电路进行进一步处理，或者可以通过一个或多个输出装置(例如电子显示器)提供给操作员。
[0047]在中子-捕获频谱分析中，将高能中子发射到围绕测井仪的体积中，在此处它们损失能量且与环境达到热平衡。在一些时间之后，这些“热”中子被原子核捕获，从而使得它们变得受激且用以发射具有特征能谱的伽马射线。捕获元素可以位于地层岩石中、地层孔隙空间中、井眼流体中、或在测井仪自身中。可以识别与每一个元素有关的特征伽马射线能谱，且分析总测量能谱用以取得每一个元素的相对贡献。经常，但不总是，受关注的元素是位于在地层中的那些。在元素——例如S1、Ca、Fe、S、T1、Gd、Mg、Al、K、以及Na——是与储油层共有(common to)的有限组矿物的一部分用于确定或估测测井仪浸入其中的岩石的岩性学的假设下，将来自所述元素的能谱产额组合在一起且总体地分析。该组假设在现有技术中作为闭合模型(例如WALK闭合)已知，且可以更改标准闭合模型用于特定的位置。中子-捕获频谱分析确定或估测的精确度涉及地层的实际元素组成(就是，组成所述地层的不同元素的存在与每一个元素的贡献百分比)与从频谱分析测井仪的测量取得的地层的元素组成之间的差别。
[0048]在一些实施方案中，背景产额可以由在频谱分析测井仪之内也能找到的这些元素中的一个或多个所产生。因此，地层中这些元素的频谱分析测量的精确度可以包含共同元素的背景产额的确定和估测以及那些背景产额的扣除。
[0049]一个背景源是来自测井仪自身之内的铁产额，其可以从总铁中扣除以获得岩石基质铁。背景铁的扣除可以涉及将总测井仪背景分组成两个标准能谱:铁标准能谱，其包含来自测井仪和岩石基质两者的贡献；以及所称的“测井仪背景”标准能谱，其对应于测井仪中除铁之外的其它元素。假设测井仪中铁与其它元素总和的比率是常数。因为这些其它元素一般不出现在通常的井下环境中，所以后一标准是整体背景水平的主要指标。
[0050]在一些实施方案中，围绕测井仪的环境可以影响中子云的尺寸和形状，且可以改变在测井仪自身、测井仪外部的岩石地层、钻井泥浆以及类似物之内热化的中子的分数。地层的环境条件可能因井而异且在特定的井之内还可以随深度而变化。不同的地层环境参数的实例包括但不限于井眼尺寸、地层密度、井眼流体密度、岩石地层的含氢指数、地层的中子慢化长度、地层的热中子俘获截面西格玛、或井眼流体的热中子俘获截面西格玛。本公开描述用于确定环境对核频谱分析测量的影响的几个不同的方法和器械，以及用于确定在地层和测井仪(和/或钻井泥浆或除了地层之外的其它材料)两者中的一个或多个元素的存在。
[0051]用于背景扣除的变量算法(variable algorithm)
[0052]地层的频谱分析测量或估测的精确度通过改善背景产额的扣除的精确度而部分增加。在测井仪的实施方案中，背景的一个源是来自测井仪自身的铁产额。在一些实施方案中，LWD工具或模块，例如前文提到的那些，包含大量的铁(即，由钢制成)。类似地，一些电缆(WL)工具也包含大量的铁(即，同样由钢制成)。在一些实施方案中，基质铁的精确测量或估测是有用的，因为用于确定或估测粘土体积的一些算法由岩石基质中的铁的量来驱动(driven)。用于确定或估测粘土体积的算法或方法论的实例在美国专利号5，786，595中描述。在一些实施方案中，对地层铁基质的测量或估测包括从能谱中的总铁产额中扣除测井仪背景。
[0053]用于扣除背景铁的一个方法可以包括首先使用频谱分析测井仪测得关于地层的测井核频谱分析数据(例如前文描述的那些)，然后通过使用具有两个标准(特征伽马-射线能谱)的总测井仪背景将测井仪背景扣除。所述两个标准能谱可以包含铁标准能谱，其包含来自测井仪和岩石基质两者的贡献；以及所称的“测井仪背景”标准能谱，其对应于测井仪中除铁之外的其它元素。测井仪之内的这些其它元素可以包括铬、镍、锰、钨、来自探测器晶体之内的元素、以及其它痕量贡献。因为这些其它元素一般不出现在测井仪外部的通常井下环境中，所以可以使用来自铁的伽马-射线产额与来自其它测井仪元素的产额之和之间的恒定比率来表示整体背景水平，其随着环境条件而不同。
[0054]测井仪背景铁产额可以由关系式:
[0055]Fe 测井仪/TB Ek,
[0056]推出，其中TB表示“测井仪背景”标准的产额，且k可以是通过实验确定或估测的常数。可以在通常的测井环境中执行上述的实验，其可以使得能谱测井工具的k值为例如
0.55。但是，如果环境的改变引起测井仪背景产额由期待以外的不同元素混合所组成的话，这个值可能不太精确。在那种情况中，真实的铁比率变成k以外的某值，且推出的或估测的基质铁测量(以及取得的粘土估测)是不太精确的。
[0057]在一些情况中，测井环境影响热中子在测井工具周围和内部的空间分布。在一些情况中，且作为结果，测井仪背景的元素组成能谱可以取决于中子在测井仪之内位于何处而不同地采样。
[0058]由于前文提到的原因，捕获原子核的可变位置还可以影响由发射的伽马射线经历的衰减的程度。结合的结果是，来自背景铁的伽马射线产额与来自其它测井仪背景元素的产额之间的比率可以偏离常数值k，且背景铁的扣除的量不够精确。本公开提供了对于给定环境而言更精确的k值，从而可以比现有技术的当前状态更精确地取得岩石基质中的铁量。
[0059]在一些实施方案中，可变背景铁比率的问题可以通过参数化铁比率作为环境条件的函数来解决，其中所述条件由可以在测井仪管柱中执行的其它测量来确定或估测。在一些实施方案中，所述函数的形式和系数可以由频谱分析测井数据、蒙特卡罗建模或其它建模方法而根据经验取得。在一些实施方案中，这两种方法可以结合在混合方法中，例如，使用蒙特卡罗模拟可以提供对环境特性的敏感性，其中所述特性原本难以仅仅由测井数据推出。适合用在本文描述的方法中的蒙特卡罗建模的示例性实例可以在Ellis和Singer的“Well Logging for Earth Scientists”第 2 版(Singer, 2007)中找到。对这些方法以及用于执行所述方法的相应的设备和系统的描述包括在下文的部分中。
[0060]由蒙特卡罗建模取得的自适应铁比率
[0061]在一些实施方案中，蒙特卡罗建模可以用于模拟中子-捕获频谱分析测量和提取背景铁比率对不同的环境参数的依赖性。所述建模可以采用分析来自MCNP的PTRAC输出的MCNP5代码和定制软件(例如，Matlab)来执行。还可以直接更改MCNP代码用以提取需要的信息。MCNP5 代码在例如来自 Los Alamos National Laboratory (2003)的 A GeneralMonte Carlo N-Particle Transport Code, Version5, LA-CP-03-0245 中描述。
[0062]可自蒙特卡罗模拟获得的信息包含中子捕获的空间位置和时间、在其上捕获它的同位素和材料、在其产生处伽马射线的统计权重(statistical weight)和能量、以及入射到探测器上的伽马射线的统计权重和能量。由这个工作的模型产生的受关注的最终数量是由探测器测量的整组元素产额，其在整个能量通道和计时门(timing gate)的合适范围上积分。
[0063]在一些实施方案中，蒙特卡罗模型可以用于产生测井环境的现实范围的预测的铁比率FewwiVTB的数据库。所述比率的功能参数化通过对这些关于有用的环境参数的数据点的变化进行的分析来提取。在这个实施方案中，可变环境参数包括地层的热中子捕获截面西格玛(Σ f)、地层慢化长度(Ls)(其是用于描述热区域之上的中子相互作用的参数)或者含氢指数(HI)、地层密度(Pb)、井眼直径(Dbh)以及钻井泥浆的西格玛值(Σ m)。在一些实施方案中，完整的数据库可以具有超过1000个点。也可能具有其它的环境依赖性。
[0064]地层慢化长度Ls可以通过多种方法确定或估测。如果使用放射性同位素中子源，可以使用至少一个中子探测器或一个伽马-射线探测器来测量Ls。或者，在不同的轴向空间(接近的探测器和远离的探测器)处的热探测器或超热探测器可以用于从接近/远离计数率的比率来确定或估测Ls (或含氢指数)。在另一个方法中，可以使用接近的伽马-射线探测器与远离的伽马-射线探测器之间的计数率比率，或者中子计数率与伽马射线计数率之间的比率。如果中子源是电子源(连续的或脉冲的)，则除了上文的探测器之外，中子监视器也可以用于确定或估测所述源的中子输出。
[0065]可以认为井眼尺寸是被钻头钻井的具有给定直径的井孔的尺寸，或者更精确地它可以通过井眼尺寸测量获得，其可以通过直接机械测径器测量(有时用于电缆工具的)或像超声测径器、密度测径器等尤其适用于LWD/MWD中的井眼尺寸测量的测量来获得。地层密度可以通过传统的伽马-伽马地层密度测量、中子-伽马密度(NGD)、声密度、以及很多其它的测量来获得。孔隙度和/或含氢指数可以由密度测量(如果已知基质密度和地层流体密度)来取得或由中子孔隙度测量来取得。地层和孔眼宏观中子捕获截面(西格玛)的测量是工业中已知的且可通过测量诱导的(捕获)中子伽马射线的衰减或测量中子总数的衰减来获得。井眼西格玛以同样的方式获得。这可以在电缆应用中采用像Schlumberger ThermalDecay Time、Reservoir Saturation Tool 及 Accelerator Porosity Sonde 的工具以及在随钻测井应用中采用Schlumberger EcoScope?工具来实现。
[0066]Feewtt/TB对每一个环境参数的功能依赖性可以通过对数据的肉眼观察来确定，或者可以使用统计的方法例如响应地面方法论。例如，且如上文所述，在一些情况中，由于地层西格玛，铁比率发生强的变化，其中大的西格玛值引起FewweVTB的减小；通过肉眼观察，这个依赖性可以通过递减指数(decreasing exponential)或线性函数来建模。可以添加额外项用以确定其它环境依赖性。在示例性能谱随钻测井工具中铁比率的可行函数(viable function)具有如下形式:
[0067]k ( Σ f，Ls, P b, Dbh, Σ m) -C1
[0068]+c2.exp (-C4.(Dbh - C1).(l_exp(-c3.( Σ f-C2)))
[0069]+C6.exp (_c5.(Dbh-C3) ).( Σ m_C4)
[0070]+c7.(Ls-C8)2
[0071]+C9.(Pb-C5),
[0072]其中独立变量分别是上文提到的五个环境特性，系数 Cl、C2、Οβ Λ C^_ Λ Λ Cg Λ Λ Cg及C9是自由参数，其可以在最小二乘拟合中同时确定或估测，且系数C1X2X3X4以及C5是示例性能谱测井仪所特有的参数。
[0073]图5示出了对地层西格玛的拟合依赖性的实例，其它环境特性保持不变。图6示出了针对慢化长度的相同参数化，且图7说明了孔眼尺寸如何缓和对西格玛所具有的依赖性的强度。这些图中的每一个均示出了总拟合数据库的子集，且每一条曲线均表示上文给出的函数的一维投影。上文给出的指数函数中的每一个项均可以根据构成测量基础的中子输运物理学来解释。例如，随着地层西格玛增加，地层中热中子的总体更快地被捕获，这减小了它们扩散到能谱测井工具部分中的比率。来自地层的向内扩散中子的减少引起测量的FeweVTB比率的整体下降。这个行为与图5中示出的行为相一致且在实验中观察到。慢化长度的变化使得热中子的空间分布发生近似的改变，除有相反影响之外，因为大的Ls允许更多的中子进入探测器附近的钻铤中。
[0074]得到的函数是对测井仪背景铁比率的环境依赖性参数化，且它可以代替常数k来执行比现在可用的背景扣除更精确的背景扣除。在一些实施方案中，使用集成的能谱测井工具，环境的测量可以用于执行实时的自适应扣除。或者，可以依赖常数比率k来实时得到回答，同时对环境变化的校正采用后续处理来完成。
[0075]蒙特卡罗建模的使用相对于用以校正元素(例如铁)的实验或数据驱动方法而言具有优势。虽然蒙特卡罗模拟是统计技术，但是蒙特卡罗模拟可以具有高的统计精度且无需面对实验噪声。模拟的数据因此具有对环境参数的敏感性，所述环境参数的影响太小以致于其不能仅仅从测井数据中提取，但是其仍然存在于真实测量中。建模的另一个优势是独立地改变每一个受关注的环境特性的能力，无需引起其它参数的相关变化。
[0076]例如，在真实地层中，在没有同时影响地层的密度和西格玛的情况下，改变孔隙度、慢化长度或含氢指数可能是困难的。相比之下，蒙特卡罗模拟提供对截面和材料组成的任意控制，从而可以设计在其中改变一个特性而其它特性保持不变的人造地层。这个方法在上文中使用。这个能力是有价值的，因为它使得可能探测环境的任意方面的单独影响。
[0077]由测井数 据而根据经验取得的自适应铁比率
[0078]取得可变铁比率的另一个方法是经验方法，其中将受关注的测井仪的测井数据与独立的数据相比较。为了取得具有足够共性(generality)的函数，这些平行数据集应该覆盖环境变化的真实范围且相关环境测量应该可用。通常的想法是可以为铁比率提供具有一个或多个未知系数的依赖环境的函数，且在迭代过程中调整那些系数直到来自受关注的测井仪的岩石基质数据最接近地匹配独立的、目标数据集。系数的调整可以通过拟合程序(fitting routine)手工或自动地执行。每一次迭代之后,每一个测量点的铁比率的新值可以用于重新计算岩石基质数据。
[0079]作为具体的实例，经验方法可以用于裸眼井，其中不同类型的能谱测井工具的数据是可用的。具有更简单材料组成和更简单构造的测井仪可以具有更少的背景扣除的困难。尽管两个不同测井仪设计的元素产额不可直接比较，但是可以比较重量浓度的最终偏差(derivation)，且特别是岩石基质铁的重量浓度的最终偏差。将较简单的测井工具(例如电缆工具)的数据进行处理，其中使它的基质产额通过闭合模型(例如WALK2，尽管不是必须的)用以产生表示目标数据的净重元素浓度。与此同时，这些数据是来自随钻测井频谱分析工具的数据集，其包含来自频谱分析和众多环境测量的元素产额。这些比较可用于几个井，覆盖宽范围的相关环境条件且总计超过10000英尺。不使用常数铁比率k，可以用一个或多个环境测量的函数来代替。
[0080]图8是由电缆工具与随钻测井工具产生的具有标称背景扣除的基质铁重量浓度之间的差别的实例，示例说明了它对地层西格玛的依赖性。地层西格玛的线性函数是经验数据的可行选择:
[0081]k( Σ f) =B^a2.( Σ f-C6),
[0082]其中，系数ai表示导致两个工具(测井仪)之间对于具体西格玛值的最佳一致性的铁比率，％是描述对铁-测井仪背景比率如何可以作地层西格玛方面校正的负斜率，并且C6是测井仪特异性参数。新的Fewwtt/TB参数化的未知系数通过自动程序而系统地改变，且在每一次迭代处，采用新的背景扣除方程将每一深度处的随钻测井工具产额重新处理且使其通过闭合模型。将每个深度处的基质铁的新的随钻测井工具重量浓度与来自电缆工具的相关铁浓度进行比较。改变所述系数以便在整个多井数据集上平均的所述两个铁浓度之间的差别尽可能地低。图9示出了这个实例中对方位平均的地层西格玛(Σ f)的线性依赖性。可以添加附加的环境依赖性，·但是对于这个数据集，仅仅使用SIFA已经在两个频谱分析测井仪之间的一致性方面产生了改进。图1OA和IOB以比较的方式示例说明了整体改进，示出了在应用环境自适应校正之前与之后的总体铁比较的柱状图。
[0083]应该注意的是，尽管在一些实施方案中多井解决方案更普遍，但是单井优化也是可以的且更易于生产。尽管岩石基质铁比较对于背景铁比率的优化而言是最简单的且最直接，但是两个数据集之间比较的数量还可以扩展到铁之外以包含其它元素浓度。在一些实施方案中，如果可获得的话，核心分析数据(core data)可以用作代替测井数据的独立数据。
[0084]自建模和测井数据二者获得的自适应铁比率
[0085]在一些实施方案中，取得可变铁比率的另一个方法是一种混合方法，其中Fe_tt/TB的环境参数化通过蒙塔卡罗建模进行开发且然后采用测井数据优化或调谐。例如，上文描述的指数函数可以应用在上文描述的经验程序中，其中改变所述系数中的一些或全部用以匹配目标数据。这个方法的简单易懂的应用是根据经验优化所述系数中的一个或多个，而同时剩余的系数能够保留通过蒙特卡罗建模已经确定或估测的值，利用其更加精细的敏感度。新的混合参数化在如上文所述的相同的多井数据集中优化。图1lA和IlB将铁比率相对于由能谱随钻测井工具测量的地层西格玛的原始的和更改的参数化进行比较。在每一西格玛值处的Fewwtt/TB的小扩展由从纯建模取得的其它环境依赖性带来。
[0086]图12是示例说明使用核频谱分析测井仪(例如上文针对图1-4所描述的)用于确定或估测地层情况的方法的流程图。通常如框1210所示，可以将核频谱分析测井仪放置到井眼中。如框1220所示，多个环境测量可以使用核频谱分析测井仪或其它设备执行。环境测量的示例说明性但是非限制性的实例包括井眼尺寸、岩石地层密度、井眼流体密度、岩石地层的含氢指数、地层的中子慢化长度、地层的热中子俘获截面、或井眼流体的热中子俘获截面中的一个或多个。
[0087]如框1230所示，可以从核频谱分析测井仪中发射中子，从而使得所述中子中的一些从与所述核频谱分析测井仪相邻的地层中产生伽马射线，且所述中子中的一些从核频谱分析测井仪之内的元素中产生伽马射线。在一些实施方案中，中子源可以是电子中子发生器(例如脉冲电子中子发生器)或化学源(例如AmBe)。在一些实施方案中，中子可以非弹性或辐射捕获相互作用的方式与地层或测井仪的元素相互作用。测井仪之内的元素可以构成测井仪的一部分，或可以在流动进入以及通过测井仪的材料之内发现。
[0088]如框1240所示，可以探测到由发射的中子诱导的伽马射线的能谱。在一些情况中，所述能谱可以包含具有多个测量的能谱成分的背景。通常如框1250所示，所述背景的至少一个能谱成分与另一个测量的能谱成分之间的背景比率可以根据所述多个环境测量进行确定或估测。在一些实施方案中，背景比率的确定或估测包括以一个或多个测量的环境参数的函数的形式确定或估测来自核频谱分析测井仪之内的铁的铁产额与总的铁背测井仪(iron back tool)背景产额之间的比率。
[0089]在一些实施方案中，用于估测背景比率的算法是基于蒙特卡罗建模、先前的实验测量或蒙特卡罗建模与先前的实验测量的组合。如框1260所示，探测的伽马射线能谱可以使用标准能谱的组合并根据该背景比率扣除所述背景的至少一个能谱成分进行分析。
[0090]在一些实施方案中，被扣除的背景能谱成分包括来自元素(例如水泥环、套管、或井管道中的钙或铁)的伽马射线。在一些实施方案中，被扣除的背景能谱成分包括发射自核频谱分析测井仪之内的铁的伽马射线。测井仪背景铁能谱可以具有相对于测井仪背景的其它能谱成分的特定估测比率。在一些实施方案中，测井仪背景的其它成分可以包含Cr、N1、Mn、W、Mo、Co以及伽马射线探测器的晶体之内的元素中的一个或多个。在一些实施方案中，被扣除的能谱背景成分具有相对于总地层岩石成分而言的特定的估测背景比率。例如，来自钢套管和水泥环的背景铁和钙可以根据相对于总地层元素的比率从频谱分析测量中扣除。
[0091]在一些实施方案中，分析经探测的伽马射线能谱包括使用所测量的环境参数进行确定或估测产生自核频谱分析测井仪之内的元素的所述能谱的部分，以及扣除能谱的那个部分用以确定或估测由岩石地层之内的元素产生的能谱部分。
[0092]在一些实施方案中，且如框1270所示，地层情况可以基于如框1260所执行的分析进行确定或估测。例如，粘土的分数数量可以使用元素-矿物换算(inversion)由元素分数进行确定或估测。由于大部分粘土矿物包含铁，所以这个换算的精确度高度依赖于元素铁分数的精确度。
[0093]使用不同空间区域的两个或多个标准[0094]在一些实施方案中，测井仪背景可以分成多重(就是说，两个或多个标准能谱或部分)用以确定测井仪的贡献，这取决于中子相互作用的不同位置。这个实施方案的一个实例取决于测井仪捕获的中子是在流动导管内部热化还是在测井仪的外部部分周围热化。这些不同的能谱可以从测量或通过蒙特卡罗建模获得。在一些实施方案中，来自这些能谱的贡献的比率可以基于包括井眼尺寸、地层密度、井眼流体密度、孔隙度、地层的含氢指数、地层的中子慢化长度、地层的热中子捕获截面、或井眼流体的热中子捕获截面的环境依赖性进行限制。基于空间区域将单个标准能谱分割成为多个标准的这个方法可以应用于由于环境(例如测井仪背景的实例)其元素组成改变的组合标准。一个特定的实施方案可以包括使用两个不同的能谱来限制变量的数目或自由程度以及任意相应的误差。
[0095]因为中子可按照取决于环境条件的比例在流动导管(或通道)内部以及在测井仪外部周围两处进行热化，所以来自测井仪部件的伽马射线不完全地或最精确地如习惯所示的由单个能谱形状来确定。热中子总体的空间分布受很多环境因素影响，包括但不限于井眼尺寸、地层密度、井眼流体密度、孔隙度、地层的含氢指数、地层的中子慢化长度、地层的热中子捕获截面、或井眼流体的热中子捕获截面。对于在测井仪的内部产生更大比例的热中子的环境，测井仪背景的最终真实能谱将看起来更加类似于内部测井仪部件的材料。在测井仪外部热化相对大分数的中子的环境将产生对外部测井仪部件更稠密(heavily)采样的背景。
[0096]在一些实施方案中，测井仪可以缺乏大部分在径向方向上的均匀性。在一些实施方案中，测井仪可以缺乏方位角和/或轴向上的均匀性。例如，可以在沿着测井仪的轴远离中子源的组成上存在不同，且中子输运可以在流动通道内部与测井仪外部之间不同。结果，测定在泥浆通道中与在测井仪外部热化的中子的相对贡献需要在测井仪长度的不同区域中采样。这个结果的示例说明性的但非限制性的实例是闪烁器晶体，其可以具有不完全屏蔽的中子截面且仅仅存在于距中子源明确(well-defined)的距离处。
[0097]在一些实施方案中，使用单个常数“测井仪背景”标准对测量的能谱进行分析可能造成困难，因为任意单个标准不总是反映测井仪背景的真实能谱形状，这引起其它元素产额中的偏差，且认为是这个其它背景的比率的背景铁方面的校正可能不精确。为了确定测井仪背景中材料的混合，对于在测井仪材料中发现的每一元素可以包括标准能谱，所述材料可以包含Cr、N1、Mn、Fe、Mo、Co、W、以及探测器晶体自身之中的一个或多个。但是，在分析中包括这么多额外的自由程度可能导致不可接受的统计不确定性。
[0098]在一些实施方案中，可以包含两个标准能谱用以确定测井仪背景，其中一个主要包含来自在测井仪内部热化的中子的伽马射线，其主要在流动导管内部，且另一个来自在测井仪的外部部件周围热化的中子。在一些实施方案中，总测井仪背景可以是这两个能谱的线性组合。为了示例说明，进行测量，其中采用富集的kiB围绕测井仪外部缠绕(wrap)用以抑制外部测井仪背景且因此测量“内部”成分，且采用流动管道(没有外部屏蔽)中装载硼砂的泥浆用以提取测井仪背景的“外部”成分。
[0099]图13将两个所得能谱进行比较，其基于被采样的测井仪部件而展现不同的元素组成。图13中，能谱1310表示内部背景标准而能谱1320表示外部背景标准。作为原理的证明，核实了对于三个不同的孔隙度(0、15、以及IOOp.u.)的总的测井仪背景可以精确地与这两个形状的不同组合相匹配。而且，可用于这个匹配的内部-外部背景标准的比率示出了对孔隙度的系统的依赖性，其表明可能不止用于单个测井仪标准且还表明了该比率可以是由其它测量的参数一例如井眼尺寸、岩石地层密度、井眼流体密度、岩石地层的含氢指数、地层的中子慢化长度、地层的热中子捕获截面、或井眼流体的热中子捕获截面——可预测的。
[0100]也可以提取多个背景标准且通过重复上文描述的实验布置的蒙特卡罗建模进行研究。图14示出了内部测井仪背景能谱1410和外部测井仪背景能谱1420。正如在实验测量中，两个假设的标准彼此不同。内部标准1410强调为流动导管材料所特有的高能峰，而钻铤中锰的6-7MeV峰为外部能谱1420中突出的。这些建模的结果与最初的实验研究完全—致。
[0101]在一些实施方案中，使用两个标准提供了至少部分地确定环境影响的优势，其中所述环境影响决定了如何对测井仪背景进行采样。最终的分析将使用更加代表真实测井仪背景的总的背景形状，这将减少背景分数和其它信号产额的提取中产生偏差的可能性。
[0102]第二个优势将来自对于测井仪背景中铁的更加精确的校正，其对于产生精确的粘土估测而言是有用的步骤。铁可以存在于测井仪和地层两者中，且组合的铁贡献可以采用单个铁标准(其和上文讨论的非铁“测井仪背景”相分离)进行分析。为了实现对地层铁的测量，测井仪铁的大的贡献可以从总的铁中扣除。铁扣除的习惯方法采取测井仪铁相对于其它背景元素所具有的固定比率。如果其它测井仪背景产额有偏差(biased)，则将其直接转换成地层铁产额的偏差。而且，正因为所述两个新提出的标准之内的元素混合物十分不同，所以相对于它们的铁比率也存在差别。如图14所示，与为0.55的整体测井仪背景的实验测量值相比较，蒙特卡罗建模对于由内部和外部测井仪背景标准表示的测井仪材料的混合分别预测了 0.33和0.60的Fe/TB比率。这表明对于铁校正使用两个标准可以更加精确。
[0103]这个实施方案应用于具有LWD工具的热中子捕获频谱分析，但是该方法还可以应用于来自非弹性中子相互作用的伽马射线并且还用于电缆工具。非弹性测量易于具有不太严重的环境依赖性，但是它们仍然存在，例如对于密度和孔隙度而言。电缆频谱分析测井仪在它们的材料组成上比用于频谱分析测井仪的一些随钻测井设计更加均匀，所以它们对测井仪背景的环境影响更少。没有流动导管和泥浆通道，中子主要在电缆工具的外部上热化。然而，可以想象其中基于空间区域的多个测井仪背景标准也可以用于电缆工具的情况。
[0104]图15是示例说明使用核频谱分析测井仪用于确定地层情况的方法的流程图。在一些实施方案中，核频谱分析测井仪可以在其材料组成方面具有空间非均质性。在一些实施方案中，核频谱分析测井仪可以是具有沿着电缆工具延伸的弓形弹簧的电缆工具或适应钻井泥浆通道(accommodate passage of drilling mud)的随钻测井工具。
[0105]通常如框1510所示，可以将核频谱分析测井仪放置到井眼中。如框1520所示，可以从核频谱分析测井仪中发射中子，从而使得所述中子中的一些从与所述核频谱分析测井仪相邻的地层中产生伽马射线，所述中子中的一些从核频谱分析测井仪之内的元素中产生伽马射线，且所述中子中的一些从钻井泥浆中的元素中产生伽马射线。在一些实施方案中，中子源可以是电子中子发生器(例如脉冲电子中子发生器)或化学源(例如AmBe或Cf)。在一些实施方案中，中子可以非弹性或辐射捕获相互作用的方式与地层或测井仪的元素相互作用。测井仪之内的元素可以构成测井仪的一部分，或可以在流动进入以及通过测井仪的材料之内发现。[0106]如框1530所示，可以探测到由发射的中子诱导的伽马射线的能谱。在一些情况中，所述能谱可以包含具有多个测量的能谱成分的背景。
[0107]通常如框1540所示，探测的伽马射线能谱可以使用包含至少两个子-标准的标准能谱的组合进行分析，其中所述子-标准表示共同元素或元素组但是其基于中子热化的位置进行区分。在一些实施方案中，所述至少两个子-标准表示来自由核频谱分析测井仪的空间不同部分中的中子热化所产生的测井仪背景或伽马射线能谱的元素的组合。在一些实施方案中，所述至少两个子-标准中的每一个包含在核频谱分析测井仪之内发现的至少两个元素的不同组合。在一些实施方案中，所述至少两个子-标准中的每一个包括Fe、Cr、N1、Mn、W、Co以及Mo的不同组合。
[0108]在其中核频谱分析测井仪是随钻测井工具的一些实施方案中，所述至少两个子-标准分别表示由在所述核频谱分析测井仪之内热化的中子产生的伽马射线和由在所述核频谱分析测井仪外部热化的中子产生的伽马射线。在一些实施方案中，所述至少两个子-标准包含表示由在流动导管内部的钻井泥浆之内热化的中子产生的伽马射线的第一子-标准和表示由在周围井眼或岩石地层中热化的中子产生的伽马射线的第二-子标准。
[0109]在一些实施方案中，蒙特卡罗建模用于产生至少一个存储的基于空间区域的标准能谱。在一些实施方案中，设计为将起源于不同空间区域中的伽马射线分隔的实验测量用于产生至少一个存储的标准能谱。
[0110]在一些实施方案中，所述至少两个子-标准中的每一个表示单个元素且按岩石地层中、核频谱分析测井仪周围的泥浆中、或核频谱分析测井仪内部的流动通道之内的泥浆中起源的区域进行区分。在一些实施方案中，所述单个元素是氢，所述子-标准分别表示岩石地层中的氢和井眼泥浆中的氢。在一些实施方案中，所述单个元素是氯，所述子-标准分别表示岩石地层中的氯和井眼泥浆中的氯。
[0111]在一些实施方案中，区别的步骤包括将由所述至少两个子-标准中的每一个表示的伽马射线衰减进行比较。在一些实施方案中，标准能谱包括一个平均能谱标准和一个或多个表示平均标准与剩余的多个子-标准之间的差别的“差别”标准。在一些实施方案中，标准能谱包括表示一个起源区域的极端之一的能谱、或极端与平均能谱之间的区域的任意线性组合、以及表示相对于这个标称标准的差别的剩余的多个子-标准。在一些实施方案中，取得标准能谱从而使得其线性组合为单个元素提供整体能谱形状的近似值。
[0112]在一些实施方案中，且如框1550所示，地层情况可以基于如框1540所执行的分析进行估测。在一些实施方案中，估测地层情况包括预测地层的组成和基于预测的岩石地层组成作出钻井决定。在一些实施方案中，在预测岩石地层的组成之后提取原油和天然气中的至少一个。例如，从测量的元素浓度确定粘土含量可以用于估测岩石的渗透性且因此用以估测可能存在的任意油的产量。
[0113]可变形状标准
[0114]在核频谱分析中，元素标准能谱用于表示单个元素或以特定比例复合的元素组的特征伽马射线能谱。单个元素标准的实例包括但不限于Fe、Cr、N1、Mn、W、Mo以及Co。复合标准的普通实例是“测井仪背景”标准，其将以被认为是有意义的比例构成测井仪的元素进行组合。传统的分析具有由单个标准能谱表示的每一个单个元素或组合元素组，无论井下环境如何所述单个标准能谱的形状均是固定的常数。如果测井仪浸入其中的环境引起真实的特征能谱具有不同的形状，则不精确的标准的使用可能引起能谱分析中的偏差。
[0115]在一些实施方案中，维持单个标准能谱但是基于对井下环境的其它测量对单个标准能谱的形状进行修改可能是有用的。能谱和其相关的算法可以被称为“自适应标准”或“可变标准”用以对能谱与常数能谱形状的假设进行区别。预测这种形状变化的条件(understanding)可以通过实验获得，例如通过在不同的孔隙度或地层以及井眼含盐量处取得相同的标准。该条件还可以通过使用蒙特卡罗建模和其它建模方法的研究获得。
[0116]自适应标准可以是复合标准(例如测井仪背景)，在这种情况中可以预测测井仪中由中子采样的元素或材料的相对混合物。还可以调整能谱用以确定在其上伽马射线散射的距离。自适应算法可以包含来自测量的信息，包括但不限于井眼尺寸、地层密度、井眼流体密度、孔隙度、地层的含氢指数、地层的中子慢化长度、地层的热中子捕获截面、或井眼流体的热中子捕获截面。自适应变化可以包括但不限于提供拉长能谱波形(smear)或移动能谱的变形(transform)的操作、调整到特定峰的高度和宽度、或者描述康普顿(Compton )散射(其中伽马射线与电子碰撞且将它能量中的一些传递给电子)的固定形状的相对分数或任意元素子-成分的变化。这些自适应可以在井下实时执行或者它们可以在后续处理中施用。
[0117]在一些实施方案中，可变标准也可以是可以受到环境影响的单个元素标准。特别受关注的实例是康普顿散射以及对来自氢和氯的伽马射线的衰减影响，其中氢和氯元素是常存在于井眼和地层两者中的元素。平均而言，井眼中的原子核位于比地层中相同元素的原子核更接近于探测器处。结果，来自地层元素的到达探测器的伽马射线易于比仅仅来自井眼中的相同元素的伽马射线在更长的距离上衰减和散射以及通过更致密的介质。观察到的每一个元素的整体能谱可以是与每一个空间区域有关的不同特征能谱的混合。虽然这些影响对于氢和氯而言是所关注的，但是该问题可适用于占据多个空间区域的任意元素。在一些实施方案中，一个解决方式是针对受关注的元素对地层与井眼成分使用单独的标准(或任意基于空间区域的其它的多个标准)。另一个解决方式是对每一个元素使用一个单个标准，其中可以改变其形状用以确定环境影响。基于康普顿向下散射(意味着在康普顿事件中伽马射线损失它们能量的一部分)以及特定峰的相对高度和宽度的调整对于这个应用而言是特别受关注的，但是也可能是其它变化。该调整可以基于上文描述的近似的环境测量进行。这个方法适于任何元素的标准能谱，包括但不限于H、Cl、S、Ba、Ca、S1、以及Fe。
[0118]在一些实施方案中，可变标准将提供更代表能谱成分的真实形状的形状，其中想要对所述真实形状在给定的环境中进行描述。这减少了能谱拟合程序中所有元素产生偏差的可能性。该方法还提供了统计精度方面的益处，因为单个标准的使用(即没有引入附加的标准)可以最小化拟合中自由程度的数目。
[0119]可变标准形状的方法可适用于来自任意中子能量状态的伽马射线，包括热中子捕获和非弹性中子碰撞。该方法可适用于任何核频谱分析测井仪，包括电缆和随钻测井仪器。
[0120]如上文所述，将单个标准分割成为多个成分标准的方法具有另一个应用，其中当多个标准表示不同的空间区域的相同的单个元素时涉及到该应用。这个方法再一次由环境影响驱动，但是来自不同的角度。然而在环境在其形状被视为常数的能谱中引起不同的元素混合之前，在此所述环境更加简单地影响单个元素的伽马射线衰减的距离。换句话说，先前的方法聚焦在采用来自所述源的中子研究变化对探测器所敏感的不同位置的增殖(populate)的影响，而后面的方法很大程度上关系到由于原子核相互作用场所相对于探测器位置的不同几何形状对能谱的影响的确定。前者更加受中子输运所影响而后者更加受伽马射线输运所影响。
[0121]一个实例是来自由在井眼中的氢相对于在地层中的氢所进行的中子捕获的伽马射线。平均而言，在井眼中的氢位于比在地层中的氢更接近于探测器处。因此来自地层氢的到达探测器的伽马射线比仅仅来自井眼的伽马射线在更长的距离上衰减和散射且通过更致密的介质。这两个不同空间区域的与氢相关的特征伽马射线能谱可以由实验测量或由蒙特卡罗建模进行提取。
[0122]图16示出来自这两个区域的氢的建模的能谱标准，再一次使用能谱随钻测井工具，且图17画出所述两个能谱的空间起源。图16中，将地层中的氢的建模的能谱标注为能谱1610，而将在井眼和泥浆通道中的氢的建模的能谱标注为能谱1620。图17中，将来自地层中的氢的建模的中子捕获标注为区域1710，而将来自井眼和泥浆通道中的氢的建模的中子捕获分别标注为区域1720和1730。由于伽马射线能量的损失平均而言来自通过更大距离以及通过更致密材料的散射，地层H的标准在更低的通道中具有更大的权重。在这个实例中，“井眼”标准包括在测井仪外部周围的井眼与泥浆通道两者。对于随钻测井工具，通过从外部井眼的贡献中分离泥浆通道贡献，也可能取得第三氢标准。
[0123]在一些实施方案中，代替对氢仅仅使用一个标准，增强的频谱分析将使用至少这两个标准。在更小井眼中或在含盐井眼中的测量将观测到来自井眼氢的更少产额，且在更少孔隙地层中的测量将观测到来自地层氢的更少产额。类似这些的环境影响引起探测的氢能谱的整体形状变化。然而通常的单标准分析尝试将常数形状施加在所有的氢上——具有偏差将增殖到所有能谱产额、而不仅仅是氢的风险——双标准分析通过可使H能谱成分以有意义的方式混合而至少部分地确定这些环境影响。
[0124]这个双区域方法的等同仪器使用描述H (包括地层和井眼两者)的平均总能谱的一个标准，以及描述地层和井眼能谱之间的差别的第二标准。两个标准的变化的总和再一次提供确定环境影响的总能谱。在一些实施方案中，两个所述测井仪标准的其它线性组合可以精确地(mathematicalIy)等同且可以使用。
[0125]这个实例讨论了氢的具体情形，但是该方法也可以改善测量出现在地层和井眼两者中的其它元素(最显著的是氯)的精确度。可以因这种方法观测到改进测量的元素包括但不限于氢、氯、硫、钡、钙、硅、以及铁。此外，受关注的空间区域不仅仅限于地层和井眼。使用实验测量或蒙特卡罗建模或其它建模方法，这个技术通常适用于可以基于任意空间区域被分割成多成分的任意单元素标准。
[0126]图18是示例说明使用核频谱分析测井仪来用于确定地层情况的方法的流程图。通常如框1810所示，可以将核频谱分析测井仪放置到井眼中。在一些实施方案中，核频谱分析测井仪可以在其材料组成方面具有空间非均质性。在一些实施方案中，核频谱分析测井仪可以是具有通过或沿着电缆工具延伸的弓形弹簧的电缆工具或者适应钻井泥浆通道的随钻测井工具。
[0127]如框1820所示，多个环境测量可以使用核频谱分析测井仪或其它设备执行。环境测量的示例说明性但是非限制性的实例包括井眼尺寸、地层密度、井眼流体密度、孔隙度、地层的含氢指数、地层的中子慢化长度、地层的热中子俘获截面、或井眼流体的热中子俘获截面中的一个或多个。
[0128]如框1830所示，可以从核频谱分析测井仪中发射中子，从而使得所述中子中的一些从与所述核频谱分析测井仪相邻的地层中产生伽马射线，所述中子中的一些从所述核频谱分析测井仪之内的元素中产生伽马射线且所述中子中的一些从钻井泥浆中的元素中产生伽马射线。在一些实施方案中，中子源可以是电子中子发生器例如脉冲电子中子发生器，或者化学源例如AmBe或Cf。在一些实施方案中，中子可以非弹性或辐射捕获相互作用的方式与地层或测井仪的元素相互作用。测井仪之内的元素可以构成测井仪的一部分，或者可以在流动进入以及通过测井仪的材料之内发现。通常如框1840所示，使用测井仪探测由发射的中子诱导的伽马射线的能谱。
[0129]如框1850所示，可以使用标准能谱的组合对探测的伽马射线能谱进行分析，其中基于所述环境测量改变至少一个标准能谱的形状用以确定环境对伽马射线能谱的影响。在一些实施方案中，具有可变形状的标准表示单个元素，例如氢或氯。在一些实施方案中，被确定的环境影响是伽马射线衰减，如由从岩石地层发射的伽马射线相对于从井眼流体发射的伽马射线的相对的混合物所驱动的。
[0130]在一些实施方案中，具有可变形状的标准表示测井仪背景中元素的集合，例如Fe、Cr、N1、Mn、W、Mo以及Co中的一个或多个。在一些实施方案中，测井仪背景中元素的集合包括在伽马射线探测器的晶体之内的元素。在一些实施方案中，测井仪背景标准的可变形状确定在形成测井仪背景的元素的相对混合物中的改变，如由在给定环境中产生的中子的空间分布所驱动的。在一些实施方案中，至少一个标准能谱的形状和变化通过蒙特卡罗建模、实验测量或蒙特卡罗和实验测量的组合取得。
[0131 ] 在一些实施方案中，且如框1860所示，地层情况可以基于如框1850所执行的分析进行确定或估测。在一些实施方案中，估测地层情况包括预测地层的组成和基于预测的岩石地层的组成作出钻井决定。在一些实施方案中，可以在预测岩石地层的组成之后提取原油和天然气中的至少一个。
[0132]在不脱离本发明的范围的情况下，能够对所描述的示例性实施方案和它们的各种不同特征作出各种修改、添加和组合。例如，尽管上文描述的实施方案涉及具体特征，但是本发明的范围也包括具有特征与实施方案的不同组合的实施方案。
【权利要求】
1.使用核频谱分析测井仪估测地层情况的方法，包括: 将包含中子源和至少一个用于伽马射线的探测器的核频谱分析测井仪放置到井眼中； 对所述测井仪和所述地层的环境执行测量； 将中子从所述测井仪发射到与所述测井仪相邻的岩石地层中，从而使得所述中子中的一些可以从所述岩石地层中产生伽马射线，所述中子中的一些可以从所述核频谱分析测井仪中的元素中产生伽马射线，且所述中子中的一些可以从钻井泥浆中的元素中产生伽马射线.采用所述测井仪探测由发射的中子诱导的伽马射线的能谱； 使用标准能谱的组合对所述探测的伽马射线能谱进行分析，其中基于所述环境测量而改变至少一个标准能谱的形状以确定环境对伽马射线能谱的影响；以及使用所述分析估测所述地层的情况。
2.根据权利要求1所述的方法，其中随所述井眼之内的深度的变化而改变所述至少一个标准能谱的形状。
3.根据权利要求1所述的方法，其中所述中子源包括电子中子发生器或化学源。
4.根据权利要求3所述的方法，其中所述中子源包括脉冲电子中子发生器。
5.根据权利要求3所述的方法，其中所述中子源包括AmBe或Cf。
6.根据权利要求1所述的方法，其中所述中子以非弹性或辐射捕获相互作用的方式与所述岩石地层、所述核频谱分析测井仪中的元素和/或所述钻井泥浆中的元素进行相互作用。`
7.根据权利要求1所述的方法，其中所述核频谱分析测井仪在其材料组成方面具有空间非均质性。
8.根据权利要求7所述的方法，其中所述核频谱分析测井仪包括电缆工具或随钻测井工具。
9.根据权利要求1所述的方法，其中所述环境测量包括对以下参数中的一个或多个进行的测量:井眼尺寸、岩石地层密度、井眼流体密度、岩石地层的含氢指数、岩石地层的中子慢化长度、岩石地层的热中子俘获截面、或井眼流体的热中子俘获截面。
10.根据权利要求1所述的方法，其中具有可变形状的标准表示单个元素。
11.根据权利要求10所述的方法，其中所述单个元素包括氢或氯。
12.根据权利要求11所述的方法，其中被确定的环境影响是伽马射线衰减，由从所述岩石地层发射的伽马射线相对于从井眼流体发射的伽马射线的相对的混合所驱动。
13.根据权利要求1所述的方法，其中具有可变形状的标准表示测井仪背景中元素的
口 O
14.根据权利要求13所述的方法，其中在所述测井仪背景中元素的集合包括Fe、Cr、N1、Mn、W、Mo以及Co中的一个或多个。
15.根据权利要求13所述的方法，其中在所述测井仪背景中元素的集合包括在伽马射线探测器的晶体之内的元素。
16.根据权利要求13所述的方法，其中测井仪背景标准的可变形状确定了形成测井仪背景的元素的相对混合物的改变，其由在给定环境中产生的中子的空间分布所驱动。
17.根据权利要求1所述的方法，其中至少一个标准能谱的形状和变化从蒙特卡罗建模得到且/或从实验测量得到。
18.根据权利要求1所述的方法，其中所述估测地层情况包括预测所述地层的组成。
19.用于估测地层情况的核频谱分析测井仪，包含: 中子源，其将中子发射到与所述测井仪相邻的地层中，从而使得所述中子中的一些可以从岩石地层中产生伽马射线，所述中子中的一些可以从所述核频谱分析测井仪中的元素中产生伽马射线，且所述中子中的一些可以从钻井泥浆中的元素中产生伽马射线； 伽马射线探测器，其探测由发射的中子诱导的伽马射线的能谱；以及 数据处理电路，其使用标准能谱的组合对所述探测的伽马射线能谱执行分析，其中基于所述测井仪和所述地层的环境的测量而改变至少一个标准能谱的形状以确定环境对伽马射线能谱的影响，并且使用所述分析用于估测所述地层情况。
20.用于估测地层情况的系统，包含: 中子源，其将中子发射到与所述测井仪相邻的地层中，从而使得所述中子中的一些可以从所述地层中产生伽马射线，所述中子中的一些可以从所述核频谱分析测井仪中的元素中产生伽马射线，且所述中子中的一些可以从钻井泥浆中的元素中产生伽马射线； 伽马射线探测器，其探测由发射的中子诱导的伽马射线的能谱； 测量装置，其用于执行所述测井仪的环境测量；以及 数据处理电路，其使用标准能谱的组合对所述探测的伽马射线能谱进行分析，其中基于所述环境测量而改变至少一个标准能谱的形状以确定环境对伽马射线能谱的影响，并且其使用所述分析估测所述地层情况。
【文档编号】G01V5/10GK103797383SQ201280036587
【公开日】2014年5月14日 申请日期:2012年5月24日 优先权日:2011年5月24日
【发明者】J·A·格劳, J·迈尔斯, M·贝雷德 申请人:普拉德研究及开发股份有限公司