CN104350233B - 使用标准化的差分数据的地层体积评价 - Google Patents
使用标准化的差分数据的地层体积评价 Download PDFInfo
- Publication number
- CN104350233B CN104350233B CN201380029177.8A CN201380029177A CN104350233B CN 104350233 B CN104350233 B CN 104350233B CN 201380029177 A CN201380029177 A CN 201380029177A CN 104350233 B CN104350233 B CN 104350233B
- Authority
- CN
- China
- Prior art keywords
- fluid
- point
- data
- geo
- displacement
- 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.)
- Active
Links
- 230000015572 biosynthetic process Effects 0.000 title description 25
- 239000012530 fluid Substances 0.000 claims abstract description 164
- 238000006073 displacement reaction Methods 0.000 claims abstract description 71
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 62
- 239000013598 vector Substances 0.000 claims abstract description 30
- 239000007924 injection Substances 0.000 claims abstract description 26
- 238000002347 injection Methods 0.000 claims abstract description 26
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 47
- 239000004215 Carbon black (E152) Substances 0.000 claims description 21
- 229930195733 hydrocarbon Natural products 0.000 claims description 21
- 150000002430 hydrocarbons Chemical class 0.000 claims description 21
- 230000035699 permeability Effects 0.000 claims description 4
- 230000005251 gamma ray Effects 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 238000005553 drilling Methods 0.000 description 44
- 230000008859 change Effects 0.000 description 38
- 238000005755 formation reaction Methods 0.000 description 24
- 239000000706 filtrate Substances 0.000 description 23
- 238000005516 engineering process Methods 0.000 description 19
- 239000000203 mixture Substances 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 11
- 230000004044 response Effects 0.000 description 11
- 238000007790 scraping Methods 0.000 description 11
- 239000008398 formation water Substances 0.000 description 10
- 239000004615 ingredient Substances 0.000 description 9
- 238000004519 manufacturing process Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 229910052500 inorganic mineral Inorganic materials 0.000 description 6
- 239000011707 mineral Substances 0.000 description 6
- 238000005481 NMR spectroscopy Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000002156 mixing Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000000638 stimulation Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 241001074085 Scophthalmus aquosus Species 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 239000007943 implant Substances 0.000 description 3
- 230000009545 invasion Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000000556 factor analysis Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000002904 solvent Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- NYNKCGWJPNZJMI-UHFFFAOYSA-N Clebopride malate Chemical compound [O-]C(=O)C(O)CC(O)=O.COC1=CC(N)=C(Cl)C=C1C(=O)NC1CC[NH+](CC=2C=CC=CC=2)CC1 NYNKCGWJPNZJMI-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 108010063499 Sigma Factor Proteins 0.000 description 1
- 229920002472 Starch Polymers 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 239000002244 precipitate Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 239000008107 starch Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010257 thawing Methods 0.000 description 1
- 230000000930 thermomechanical effect Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/003—Determining well or borehole volumes
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/164—Injecting CO2 or carbonated water
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
提供了一种基于标准化的差分数据集确定地质地层内的流体的体积数据的方法。该方法包括:基于在相应不同的第一和第二时刻从井眼的测量收集地质地层的第一和第二数据集瞬像,并且使井眼在第一和第二时刻之间经受用于驱替与井眼邻近的地质地层中的流体的流体注入;基于第一和第二数据集瞬像生成差分数据集;标准化差分数据集,以生成标准化的差分数据集;基于标准化的差分数据集确定限定出几何形状且与相应的不同的被驱替流体标记对应的顶点;确定通过代表具有已知的第一特性的第一被驱替流体的第一点且沿着相应的第一矢量指向的第一条线;确定通过代表具有已知的第二特性的第二被驱替流体的第二点且沿着相应的第二矢量指向的第二条线;基于第一条线和第二条线的交点确定与注入流体的特性对应的注入流体点;确定通过注入流体点且沿着与具有至少一种未知特性的另一被驱替流体对应的另一矢量指向的另一条线;基于其它被驱替流体的至少一种已知特性确定沿着其它线的第三点;并且基于差分数据集、第一点、第二点以及第三点确定被驱替流体的体积成分。
Description
背景技术
测井工具可被用在井眼中,例如用于进行地层评价测量以推断围绕井眼的地层的特性以及地层中的流体的特性。常用的测井工具包括:电磁工具、声工具、核工具以及核磁共振(NMR)工具,当然也可以使用其它各种工具类型。
早期的测井工具在井眼被钻好以后设置在电缆上被送入井眼。现代版的这种电缆(WL)工具仍然在广泛使用。然而,对于钻井过程中的实时或准实时信息的需求催生出随钻测量(MWD)工具以及随钻测井(LWD)工具。通过在钻井过程中收集并且处理这些信息,钻井者可以修改或校正井操作的关键步骤使得钻井性能和/或井轨迹最优化。
MWD工具通常提供钻井参数信息,例如钻压、转矩、冲击和振动、温度、压力、每分钟转数(rpm)、泥浆流率、方向及倾角。LWD工具通常提供地层评价测量,例如天然或光谱伽马射线、电阻率、介电质、声速、密度、光电因子、中子孔隙率、西格玛热中子俘获截面(Σ)、各种中子诱发伽马射线谱以及NMR分布。MWD和LWD工具通常具有电缆工具共有的部件(例如发射和接收天线或通常的传感器),但是MWD和LWD工具可以构造成不仅仅能经受而且能在恶劣的钻井环境里运行。术语MWD和LWD通常可以互换使用,并且在本公开中使用的每个术语都可以被理解为包括地层和井眼的信息的收集以及钻井组件的移动和布置的数据的收集。
测井工具可用于确定地层体积特性,也就是量化所研究地层的给定样本中存在的每种和各种组分的体积分数,通常采用百分比表示。地层体积特性分析包括存在组分的鉴别以及在不同的测井测量上为组分指定唯一的标记。当使用相应的地模型时,所有各种组分的正演模型响应被校正,测井测量结果可被转换成组分的体积分数。
发明内容
本发明内容的目的是介绍可供选取的构思,其在下面的详细说明中进一步描述。本发明内容不用于确定要求保护的主题的关键或基本特征,也不用于帮助限制要求保护的主题的范围。
一种用于确定地质地层内的流体的体积数据的方法,所述地质地层中具有井眼,该方法可以包括:基于在相应不同的第一和第二时刻从井眼的测量收集地质地层的第一和第二数据集瞬像,并且使井眼在第一和第二时刻之间经受用于驱替与井眼邻近的地质地层中的流体的流体注入。该方法还可以包括基于第一和第二数据集瞬像生成差分数据集;标准化差分数据集,以生成标准化的差分数据集;基于标准化的差分数据集确定限定出几何形状且与相应的不同的被驱替流体标记对应的顶点。该方法还可以包括确定通过代表具有已知的第一特性的第一被驱替流体的第一点且沿着相应的第一矢量指向的第一条线;确定通过代表具有已知的第二特性的第二被驱替流体的第二点且沿着相应的第二矢量指向的第二条线;基于第一条线和第二条线的交点确定与注入流体的特性对应的注入流体点;确定通过注入流体点且沿着与具有至少一种未知特性的另一被驱替流体对应的另一矢量指向的另一条线。该方法可以附加地包括基于其它被驱替流体的至少一种已知特性确定沿着其它线的第三点;并且基于差分数据集、第一点、第二点以及第三点确定被驱替流体的体积成分。
同样还提供一种相关的测井系统以及非易失性计算机可读介质。
附图说明
图1为可以用于实施一个示例性实施例的井场系统的示意图;
图2和3为描绘根据示例性实施例的地层评价操作的流程图;
图4为与通过流体驱替彼此进行组分替代的单个组分对对应的三维(3D)图。
图5为示出了从时延地质地层瞬像确定差分数据集的示意图。
图6-9为示出了图5的差分数据集的流体驱替标记的3D图。
图10为示出了图9的流体驱替标记被标准化成统一长度的3D图。
图11和12为分别示出了图10的标准化标记点投影在假象球上以及产生的连接各点的测地三角形的示意性3D图。
图13和14为示出了与相同于图4的通过流体驱替彼此替换的单个组分对对应的数据点的3D图,但是这些数据点的相应投影和由此产生的标准化流体标记分别位于水平(X,Y)、垂直正面(Y,Z)以及垂直左面(Z,X)平面上。
图15-17为示出了另一种绘制图12的标记点的方法的二维(2D)图。
图18和19为示出了根据一个示例性实施例的用于确定钻井泥浆滤液以及天然地层烃标记的方法的3D图。
具体实施方式
本说明书参照附图进行说明,附图示出了示例性实施例。然而,也可以使用许多不同的实施例,因此说明书不应该被解释为对这里提出的实施例的限制。与此相反,提供这些实施例以使得本公开充分且完整。在整个说明书中相同的附图标记指代相同的部件。
开始参照图1,首先对可用于实施这里提出的示例性实施例的井场系统进行描述。该井场可以是陆上或海上的。在该示例性的系统中,通过旋转钻井在地下地层106中形成井眼11。本公开的实施例例如还可以使用定向钻井。
钻柱12悬挂在井眼11内且具有底部钻具组合100,该底部钻具组合包括位于其下端的钻头105。地面系统包括设置在井眼11上方的平台和井架组件10,该组件10包括转盘16、方钻杆17、钩18以及转环19。钻柱12被转盘16转动,所述转盘16将方钻杆17接合在钻柱的上端。钻柱12通过方钻杆17和转环19悬挂在被连接到游动滑车(未示出)的钩18上,所述转环19允许钻柱相对于钩旋转。在某些实施例中,还可以使用顶部驱动系统。
在所示的例子中,地面系统进一步示例性地包括储存于形成在井场的池27中的钻井流体或泥浆26。泵29输送钻井流体26使其通过转环19中的端口进入钻柱12的内部,使得钻井流体如方向箭头38所示地向下流动通过钻柱12。钻井流体通过钻头105中的端口离开钻柱12,然后如方向箭头39所示地向上循环通过钻柱的外侧和井眼11的壁之间的环形区域。钻井流体润滑钻头105并且随着其返回池27用于再循环而将地层106的钻屑带到地面。
在多个实施例中,本文公开的系统和方法可以与本领域普通技术人员公知的任何传输方式共同使用。例如,本文公开的系统和方法可与通过电缆、钢丝绳、钻杆传输、连续管钻井和/或随钻传输接口传输的工具或其他电子设备共同使用。图1示出了一种随钻接口,这仅为了示例性目的。然而,本文公开的系统和方法同样可以应用于电缆或任何其他合适的传输平台。所示实施例的底部钻具组合100包括随钻测井(LWD)模块120、随钻测量(MWD)模块130、旋转导向系统和马达以及钻头105。
LWD模块120容纳在钻铤中,且可以包含一种或多种类型的测井工具。还应该理解,可采用一个以上的LWD和/或MWD模块,例如如附图标记120A处所示。(贯穿整个申请中,对120位置处的模块的引述同样可以替换性地表示120A位置处的模块)。LWD模块可具有用于测量、处理和储存信息以及用于与地面仪器、例如所示的测井和控制站160通信的能力。作为举例,LWD模块可以包括以下中的一个或多个:电磁装置、声装置、核磁共振装置、核测量装置(例如,伽马射线、密度、光电因子、希格玛热中子俘获截面、中子孔隙率)等等,当然也可以使用其它测量装置。
MWD模块130也容纳在钻铤中,且可包含一个或多个用于测量钻柱和钻头的特性的装置。MWD工具还可包括用于向井下系统产生电能的装置(未示出)。这典型地可以包括由钻井流体流驱动的泥浆涡轮发电机,应该理解,也可采用其他电能和/或电池系统。MWD模块也可包括以下类型的测量装置中的一种或多种:钻压测量装置、扭矩测量装置、冲击和振动测量装置、温度测量装置、压力测量装置、每分钟转数测量装置、泥浆流率测量装置、方向测量装置以及倾斜测量装置。
根据示例性实施例,上面描述的井眼工具可以用于收集与井眼11相邻的地质地层的测量结果,以确定地质地层106内被驱替的流体的一个或多个特性。处理器170可被提供用于确定这些特性。处理器170可以使用硬件(例如微处理器,等等)以及具有用于执行这里描述的各种操作的可执行计算机指令的非易失性媒介的组合来实施。应该注意到的是,处理器170可以位于井场中,或者可以远程设置。
作为背景,地层评价(FE)的目标之一是地层体积分析,即,研究的地层的给定样本中存在的每种组分的百分比体积分数的量化。地层体积分析的核心是存在的组分的鉴别以及相应的地质模型(有时也被称作“地模型”)。组分在不同的测井测量上被指定一个标记,并且选取的测井测量结果通常被优化以确保每个存在的组分具有唯一的标记。总的来说,实际的考虑,例如技术、操作条件(井几何特征、井眼尺寸、泥浆类型、裸眼井与套管井、温度等等)、HSE方面以及经济学可能限制预期的测井测量结果。此外,基于选取的测量的固有物理性质选取均匀介质“混合定律”,基于执行测量的具体工具类型和设计选取三维几何响应函数。同时考虑地层组分测井测量标记、混合定律以及几何响应函数使得对组分混合物的各种测井测量响应能够正演建模,正演模型的反演然后可将测井测量结果转换回组分体积分数。
特殊地,鉴别和为存在的不同组分指定测井标记(在现场条件下)的操作可能是一种挑战,尤其是在传统的过平衡钻井的情况下存在相对较深的侵入深度而采用具有相对较浅的勘测深度的WL测井工作时,尽管侵入之前获取的LWD测量可能已经太深地进入到地层内和/或欠平衡钻井可用于减缓这些WL特别关注点。然而,尽管通过各种操作可对存在的不同组分的鉴别进行在一定程度修正,但为存在的不同组分指定唯一的标记并不总是具有容易的方法。这可能归因于多种因素。
例如,对在钻井过程和/或泥浆测井操作过程中带回到地面的岩屑的分析通常可以为地质学者和油藏工程师提供有关存在的不同组分的鉴别的重要且提前的线索(这里被称作“地面真相”),其中,具有某些例外(取决于钻井泥浆类型)。可选的钻取岩心操作(其有可能是昂贵且不现实的)为更深一步的操作,用于切割且取出许多英尺的地层整个岩心,以在地面进行进一步的详细分析。同样,井下进行的元素谱测井技术(例如,热中子俘获谱测井、快中子非弹性散射谱测井、元素中子活化谱测井,等等)都能够帮助说明基质组分,并且使地层体积分析挑战简化到仅仅是流体基本体积分数。
此外,可选的地层测试操作(例如,压力梯度、井下流体分析、流体采样,等等)(尽管在沿着井的离散深度点上的这种位置数据的可获取性受到限制)可以被考虑用于测试地层的可生产流体组分。同样,最近引入的先进的多维度NMR测井技术可以帮助鉴别出不同流体组分。
为特殊的组分指定标记的一个先决条件是其通过数学分析与其它组分完全或实质上分离且隔离的体积量(或质量)。对这种样品作出的测量结果接着可被标准化成存在的组分的量,且得到测井标记。应该注意的是,即使当样品在地面被回收时,用于执行与各种井下测井类似的测量的地面仪器也不是轻易能得到或可能的,并且即使这样,在地面执行的测量需要被进一步进行推延到井下压力和温度条件。
这里提供了一种系统方法用于通过单独的测井测量结果来鉴别和校准某些地层组分测井响应。也就是说,不是一次在一个深度寻找存在的各个组分的标记,相反,当替代成对地发生时(即,当一种组分“I”替代另一种组分“J”,而其它所有东西保持一样时)本方法可以寻找由交叉组分(x-组分)替代产生的样式。这有效地成为一种组分与另一种组分对照,并且当其中一种组分测井响应被充分理解时,可以重新构建另一种测井响应。
一种示例性的用于确定地质地层106内的流体的成分数据的实施方式首先总体地参照图2的流程图200进行描述。从方框201开始,在方框202处,该方法示例性地包括基于在相应不同的第一和第二时刻来自井眼11的地质地层106的测量结果收集第一和第二数据集瞬像,并且使井眼经受第一和第二时刻之间的流体注入以驱替与井眼相邻的地质地层中的可移动流体。作为举例,流体注入可以包括各种类型的提高油采收率(EOR)的流体,例如新鲜水、二氧化碳,等等。该方法可进一步包括在方框203基于第一和第二数据集瞬像生成差分数据集,并且在方框204使差分数据集标准化以生成标准化的差分数据集,这将在下面进一步描述。该方法还示例性地包括在方框205基于标准化的差分数据集确定限定一种几何形状并且与相应的不同的被驱替流体标记相对应的顶点,并且在方框206基于来自第二数据集的数据点在几何形状上的位置来确定相对于不同的被驱替流体标记的被驱替流体成分数据,这将在下面进一步详细描述。该方法在方框207示例性地结束。
更特殊地,本方法采用有效的一致性测量。这是通过处理技术、例如侵入校正技术的真正的一致,或者实质上的一致,或者这是因为尽管勘测的实际体积可能不同,仍然在相同类型的地层中读取测量结果。例如,这可能发生在测量结果同时处于它们很少受到侵入影响的情况下,或者处于它们全部被侵入覆盖的情况下。这些测量结果用于探测相同的地层两次或更多次,其中,在不同的探测或瞬像中间期望地层成分发生改变。这允许对已经发生的改变进行表征。应该注意到,对于相同的瞬像,测量结果仅仅需要在彼此之间一致。来自一个瞬像的测量结果与来自另一个瞬像的测量结果不需要一致。
虽然开始的时候可能看起来好像那样的话问题会变得更复杂,但是这也不是必然的情况。例如,对于存在的“Z”组分,就可能会有“Z(Z-1)”组分对交换(远大于Z),但是在自然和实际中,仅有非常小数量的这种对交换会与手头的案例相关。作为举例,当前的储层内的白天天然流体分布(作为在地质时标范围上的流体迁移和替代的结果)以及相对渗透率随着相应流体的饱和度的上升,使得在指定深度下仅有其中一种原始地质天然流体是显著可移动的。也就是说,其它的已经被驱替。此外,扰乱这种初始储层平衡(或者均衡的流体分布)的侵入流体通常被很好地确定是从地面被注入或者产生到地面。
另一方面,通常难以直接隔离开各个流体组分的标记,因为在研究的地层体积中,它们不可能独自出现,或者它们可能没有足够的量可用,不管是否是上面讨论的储层平衡。这通常是过平衡钻井的情形,并且被传统的WL测井加剧。相反,欠平衡钻井应被考虑,或者测井测量应被认为适合于现有的侵入校正技术(例如在Gzara申请的美国专利公开号为No.2009/0177403中描述的方法,其由此通过引用整体被结合到这里),此时情况将不同,并且一种类型的流体组分实际上可覆盖全部的其它组分。然而,即使在这种情况下,缺少存在的流体组分的精确量的信息通常将表现出妨碍该流体标记的导出,但这可以采用这里提出的方法进行克服,这将在下面进一步讨论。
此外,当研究由x-组分替代产生的样式时,其它组分明确地不会起作用,这就减少了否则可能来自于同时尝试求解大量的组分测井测量标记的复杂性。然而,有一种特殊的情况,其中x-组分替代不需要完全地成对发生,而仍然可以实施并且适于采用这里提出的构思。这种特殊的情况就是地下地层具有变化的水盐度,通常是由用于维持储层压力并且维持烃产量而执行的水注入操作导致的。这里,注入的水盐度与原始的地层水(也被称作“原生”水)盐度明显不同,并且这两种盐度在整个储层中以不同比率混合产生不同的水盐度。在这种情况下替代的流体可以被理解为原生地层水、注入水以及未波及的烃的混合物。
这种提出的方法还可以应用于大范围的情形,取决于在不同瞬像之间观察到的地层成分改变的许多可能的源头。事实上,观察到的改变可以是被驱替流体、被驱替碎屑、相变(例如由压力或温度变化引起的)、或者总体上包括有溶解或沉淀的化学反应(例如沥青质沉淀,结垢,盐溶解,酸刺激等)或者总体上压实或压力或应力状态的最终改变的结果。
通常来说,这种改变可能归于各种类别。第一类别为随时间的改变(例如,当在不同时刻探测相同体积的地层时,第一次通常被称作“基础测井”)。关于注入诱发的改变,这些可能包括:小时标,侵入动态(钻道次与刮道次);小时标,储层增产技术(例如与化学反应动态,或溶剂注入关联的侵入);小时标,总体上的测井-注入-测井(LiL)技术(即,多个侵入循环,具有适合用途的侵入流体);以及大时标,储层监测(例如采用注入井)。关于生产诱发的改变,这些可能包括:小时标,欠平衡钻井,或者压力诱发的改变(例如气体膨胀,冷凝水堆积,从溶液出来的气体,气锥,水锥,或漏失带);以及大时标,储层监测(例如采用生产井)。其它的改变为“热机械配置”诱发的改变,其可能包括:小时标,温度诱发的改变(例如冰或水合物的熔化以及融化);大时标,温度诱发的改变(例如当使用热回收技术时,刺激重油特性);以及大时标,应力诱发的改变。
下一个类别包括随径向深度的改变(例如,当仅在一次探测相同地层的越来越深的体积),其对于研究的每个越来越深的体积需要在彼此之间具有不同的一致测量组。关于注入诱发的改变,这些可能包括:小时标,侵入动态(钻过与刮过);小时标,储层增产技术(例如与化学反应动态,或溶剂注入关联的侵入);小时标,总体上的LiL技术(例如,多个具有适合用途的侵入流体的侵入循环)。就生产诱发的改变来说,这些可能包括小时标,欠平衡钻井以及压力诱发的改变(例如冷凝水堆积,或者从溶液出来的气体)。关于全部的“配置”诱发的改变,这些可能包括小时标,温度诱发的改变(例如冰或水合物的熔化以及融化)。
又一个类别包括带间的改变(例如随着深度的改变),其中,具有一种相同的组分并且其参与到全部的预料的x-组分对的替代中。这是一种与直觉稍微有点不同的情况,仅仅当在不同带上存在的相同组分被确定具有相对的置信度时才可以应用。在这种情况下,在指定深度得到的测量结果被与假想的相同的组分占据地层的整个体积的情况对照,这就是该技术如何延展到这种情况。即使当相同组分的性质仅仅是近似已知时,单单我们面对相同组分的事实也足以使该技术起作用。实际中,可以基于仅仅主要地响应于岩石和矿物质的井下测井数据区分相同的岩石矿物学特征,例如(但不限于)先进的元素俘获谱,或者自然伽马射线测井数据。还可以基于地面的观察而进行区分,例如(但不限于)通常的岩心数据,以及泥浆测井数据和特殊的切屑分析。可选地,可以基于仅仅主要地响应于流体的井下测井数据区分相同的流体类型,例如地层测试测井数据。还可以基于地面观察进行区分,例如(但不限于)通常的生产流体分析,以及更特殊的泥浆测井数据和钻井泥浆回流分析。或者还可以简单地确定,因为其可以从地面注入,例如(但不限于)在欠平衡钻井的情况下的钻井泥浆滤液。
当岩石矿物学特征可以明确区分时,流体类型的改变可以被识别,并且当流体类型的改变还伴随有孔隙率的显著变化时,可以现场校准关注的岩石矿物学特征的端点。相反,当流体成分可以明确区分时,岩石矿物学特征的改变可以被识别,并且当岩石矿物学特征的改变还伴随有孔隙率的显著变化时,可以现场校准关注的流体类型的端点。也可以使用前面的各种组合。
应该注意到,与地层评价相比,生产测井或者钻井优化的原理关注于在生产或注入或在钻井过程中井眼自身的容纳物,而不是地层的组分。本领域技术人员可以意识到,这里描述的一些构思例如可以被转用到生产测井或钻井优化的领域(例如钻孔清理以及井涌检测)。
根据第一方面,描述了一种鉴别和分类已经发生的改变的方法。使用与有效一致测量所考虑的m1m2…mαmβ…mn对应的矢量符号并且描述将涉及如的地层的不同瞬像,而不同的地层组分测井标记称为此外,通常表示本身或其任何的线性变换。在某些组分的体积和测井响应已知的为优选的时候,符号还将包括使摆脱这些已知组分的贡献的以产生仅仅单独取决于剩余的未知的“干净的”矢量的这种变换。
在本说明中,这些矢量可替代性地展示为“n”数据点上的曲线,带有值m1m2…mαmβ…mn,在这种情况下,可放弃使用矢量符号并且替代为函数符号以及这就是NMR多分量数据通常如何显示,并且根据相关的曲线已创造出术语“分布”。在本说明书中,测量结果m1m2…mαmβ…mn通过将全部的测量结果相对于每一个固有的噪声量标准化而转换成无单位(或无量纲)的。首先,这有助于将各个测量结果保持在固有噪声水平上方,并且避免噪声与真实信息混杂。其次,当其涉及到显示上面讨论的矢量或函数时在中性或用户独立的刻度上是有用的。应该注意的是,这种测量结果标准化不同于后面介绍的其它标准化,例如标记虚拟标准化,以及标记真实标准化。
在瞬像“i”和“j”之间的的改变接着可以如下表达成全部的矢量的线性组合(假设测量结果具有线性混合定律):
需要注意的是该表达式不是唯一的,因为矢量是互相依存的。在组分“I”和“J”仅仅是成对互换的情况下,更熟悉的表达式如下:
取该差值的好处也可以在现参照图5讨论的例子中呈现出来,图5示出了在过平衡钻井过程中在钻井泥浆滤液侵入的环境中与彼此相减钻道次和刮道次对应的过程。图的上部分“(a)”示出了在钻道次(drill pass)过程中构成基质(-Matrix-)的矿物质(Min-1,Min-2和Min-3)的体积分布,以及填补所考虑的LWD测量所研究的体积内的孔隙空间(-Phi-)的流体(Fld-A,Fld-B和Fld-C)的体积分布。在这种情况下,来自钻道次的LWD测量结果被看作是与存在的这些矿物质和流体中的每一种相对应的相同的测量响应的线性组合,采用它们相应的体积比例进行加权。
图的第二(中间)部分“(b)”示出了在刮道次(wipe pass)过程中构成基质(-Matrix-)的矿物质(Min-1,Min-2和Min-3)的体积分布,以及与填补所考虑的LWD测量所研究的体积内的相同的孔隙空间(-Phi-)的原始流体(Fld-A,Fld-B和Fld-C)近旁的另一种流体(Fld-X)的体积分布。流体Fld-X(例如注入钻井泥浆滤液)表示不是在孔隙空间内原始存在但其现在占据着由流体Fld-A,Fld-B和Fld-C最初占据的孔隙空间的一种新流体。在此,来自刮道次的LWD测量结果被看作是与存在的这些组分中的每一种对应的相同的测量响应的线性组合,采用它们相应的体积比例进行加权。注意到在该例子中,矿物质的体积分布在钻道次和刮道次之间没有改变。
图的最后(下)部分“(c)”示出了与钻道次和刮道次测量结果之差(即差分数据集)对应的体积分布。注意到基质矿物质(以及在钻道次和刮道次间没有移动的任何其它物质)抵消了。再一次,来自钻道次和刮道次的LWD测量结果之差看作是标记的线性组合,其现在不与存在的各个组分对应,而是与彼此交叉替代(Sig-I,Sig-II,Sig-III)的组分对的标记对应。也就是,这是一种比其它组分的标记少的其中一种组分的测井测量标记,采用相应被驱替体积进行加权。
转到参看图6-8,它们与图4类似并且展示出与三种不同的流体替代模式(由点60表示的泥浆滤液替换Fld-A,由点61表示的泥浆滤液替换Fld-B,以及由点62表示的泥浆滤液替换Fld-C)对应的关系,并且在这种情况下表示三个测井测量结果PhiN(表观中子孔隙度),PhiD(表观密度孔隙度),以及PhiΣ(表观Σ孔隙度)。图9示出了在相同的图表上同时展示的全部三个不同的流体替代标记点60-62。
该结果意味着与相同的“I”和“J”对交换对应的数据点将沿着矢量对齐,反之亦然。接着,沿着这些矢量的数据点簇鉴别出哪对地层组分“I”和“J”已经在瞬像"i"和"j"之间彼此进行了替代。为了在实际中有效地区分这些簇,一个方法就是考虑在“n维”空间中每个立体角的数据点柱状图,或者使数据点矢量根据以下公式标准化成为幅度1(即,将它们相对于半径1的n维球体投影):
对于那些位于预设噪声阀值之上的数据点矢量并且其中标准可以采用多种方法限定。这种虚拟标准化明确地揭示出存在的某些x-组分替代模式,其中,这种替代在不同的地层瞬像之间已产生显著的差别。神经网络技术,要素分析,和/或其它统计分析技术此时可以被用于根据掌握的模式自动地划分地层。
在地下地层的水盐度变化的特殊情况下,这通常源于用于维持储层压力和维持烃生产的水注入操作,其中,注入的水盐度明显不同于原始地层的水(即天然水)盐度。以不同比例在储层上将两者混合产生不同的水盐度。一旦已鉴别和/或提取出天然地层水、注入水以及天然地层烃的标记,那么就可以将钻道次和刮道次之间的测井测量结果差值转换成位于被泥浆滤液驱替的地层流体体积内的天然地层水、注入水以及天然地层烃的相应的比例。
以这种方式实现的被驱替流体成分在这里被称作“虚拟成分”。这种虚拟成分单独授予每种流体组分,即,当仅有一种流体被驱替时,虚拟成分将仅仅单独表示这种组分,相反,当一种流体还没有被驱替时,虚拟成分表示不存在该组分。然而,虚拟成分是非线性的并且不会准确地授予多流体间的混合物。虚拟成分自身可以按照多种方式执行,取决于使用的虚拟标准化。一种方式是通过将流体标记定位下面描述的测地三角形内而推导出成分数据,这通过显示出的标记(即,矢量SIG-I,SIG-II以及SIG-III)支持。
虚拟标准化的一种考虑是在相对应的矢量彼此平行的那些情况下,来自不同的x-组分替代模式的数据点簇一旦被标准化后无法彼此区分开。此外,围绕原点“O”聚集并且与具有相似特性的x-组分对(例如被油基泥浆滤液驱替的天然地层油,或者被水基泥浆滤液驱替的天然地层水)对应的数据点簇,可能最后不能与对应于其它x-组分对交换的其它数据点簇区分开,并且当仅仅将那些数据点矢量保持在预设的噪声阀值之上时将无法达标。
参照图10,在此,在图9中所示的三个不同的线和流体替代标记点60-62再次被示出,但是还示出了沿着等于1的距离设置的三条线的相应标准化点70-72(即这些线与半径等于1的球相交的点和/或数据点在半径等于1的球上的投影点)。由于所示的线和相应的与半径为1的球的交叉点之间一一对应,对不同流体替代标记的引述被解释为表示位于半径为1的球上的相应的点70-72。在图11中,只示出了上面讨论的半径为1的球和标准化的点70-72(即,相应的线被移除,其可以被看作位于该点上的冗余信息)。
在图12中,示出了连接不同的标记点或顶点70-72的测地三角形。根据分别被点和两个相对的顶点保持的“立体角”(或区域)与被全部三个顶点70-72保持的立体角之间的比率,在该三角形区域内包含的任何点75实际上都对应于以不同比例替代Fld-A、Fld-B、Fld-C的混合物的泥浆滤液Fld-X的标记。
此外,参照图13-17,示出了将三维(3D)空间中的数据点转换成二维(2D)空间中的相应的表示的过程,在这种情况下,3D空间中的单个点可替代地表现为2D空间中的三角形。参照图13,这示出了与图4所示的单个流体替代标记对应的相同的线和数据点,但是现在具有添加的这些数据点在三个平面XY(水平平面)、YZ(垂直正面平面)以及ZX(垂直左面平面)中的每个上的投影。在图14中,该视图与图13类似但是现在还包括位于半径为1的球上的流体替代标记点70,以及如上所讨论的在三个平面XY、YZ和ZX中的每一个上的相应的投影90-92。
在图15中,通过使平面XY、YZ和ZX上的不同2D投影彼此重叠而将图13和14的3D显示用2D显示取代。在图16中,示出了形成三角形并且连接单个流体替代标记点70的不同投影90-92的线。由此,如图17所示,来自不同数据集的3D数据点可以替代地以2D中的相应的三角形表示。
关于从3D显示向2D显示转换的过程,其中流体替代标记替代地采用三角形以2D代替3D点表示,在某些实施例中,可能更加适合采用这种2D显示。这可能是这种情况:当采用多于三个测井测量结果(即多于三个维度)工作时,在这种情况下,N-维流体替代标记可以可选地被转换成2D标记,该2D标记采用“N x(N-1)/2”多边形表示。
现在另外参照图3的流程图300,在某些实施方式中,可能期望的是通过定位来自由图9-10上的原点“O”和点60-62(如果点60-62也可准确地被鉴别,而不仅仅是上面讨论的虚拟标准化的焦点的点70-72)支撑的四面体内的差分数据集,还同时考虑被驱替流体的真实成分和已经侵入地层的泥浆滤液的体积。这例如在地层水盐度可变的情况下是可能的,因为水是公知的流体。从方框301开始,在方框302中,在相应不同的第一和第二时刻从井眼11收集地质地层的第一和第二数据集瞬像(即,钻探和刮擦瞬像),其中,井眼经历在第一和第二时刻之间用于驱替与井眼相邻的地质地层中的流体的流体注入。与上面讨论的类似,基于第一和第二数据集瞬像生成差分数据集(方框303),差分数据集被标准化以生成标准化的差分数据集(方框304),并且在方框305中,基于标准化的差分数据集确定限定出几何形状并且与相应的不同的被驱替流体标记对应的顶点。
此外,参照图18-19,引入了新的点80-82并且它们分别与点60-62并置,以辨别坐标位于差分数据集参照(以3个轴标签ΔPhiD、ΔPhiN以及ΔPhiΣ示出)中的点60-62和坐标位于第一和第二测量数据集瞬像绝对参照(以3个轴标签PhiD、PhiN以及PhiΣ表示)中的点80-82。这种区分在矢量(以及顶点)的情况下是不需要的,因为矢量将在上述两种参照中保持相同的坐标。还引入了位于差分数据集参照的原点处的点83,并且点80-83坐标分别表示存在的所有流体:天然地层流体Fld-A(例如,地层油)、Fld-B(例如,含盐原生水)、Fld-C(例如,新鲜注入水)以及钻井泥浆滤液Fld-X的特性。
除了图6-17中所用的差分数据集参照(以3个轴标签ΔPhiD、ΔPhiN以及ΔPhiΣ示出),在图18和19中还示出了第一和第二数据集瞬像绝对参照(以3个轴标签PhiD、PhiN以及PhiΣ表示)。采用圆形示出的各个数据点具有不同的坐标系,取决于考虑的差分或绝对参照,而矢量在所述两种参照中都保持相同的坐标。
在示例性的例子中,在方框306中,第一条线101被确定为通过表示具有已知的第一特性(例如,Fld-B)的第一被驱替流体的第一点81,并且沿着相应的第一矢量(例如,Sig-II)指向。此外,在方框307中,第二条线102被确定为通过表示具有已知的第二特性(例如,Fld-C)的第二被驱替流体的第二点82,并且沿着相应的第二矢量(例如,Sig-III)指向。在方框308中,与注入流体(例如Fld-X)的特性对应的注入流体点83基于第一条线101和第二条线102的交点被确定。在方框309中,另一条线100被确定为通过注入流体点83并且沿着另一个与具有未知特性的另一种被驱替流体(例如Fld-A)对应的矢量(例如Sig-I)指向。在方框310中,具有未知特性点80的被驱替流体接着基于被驱替流体的至少一种特性(例如,密度,或者API重力)被确定为沿着线100。这允许在方框311中,被驱替流体的体积成分可以基于差分数据集和点80-83确定。在某些实施例中,在方框312,地层或储层特征(例如,渗透率、相对渗透率、分流量等)还可以基于确定的被驱替流体的体积成分确定,其示例性地结束图3的方法(方框313)。
更特殊地,在得到了原生地层水(例如,Fld-B)和注入水(例如,Fld-C)的盐度后,可以计算相应的测井测量响应81和82。此外,借助于之前描述的通过时延数据采集推导出的与具有泥浆滤液(例如,Sig-II和Sig-III)的x-组分替代的标记对应的两个矢量,我们现在在3D空间中具有两条线101、102。这两条线在泥浆滤液的标记点83处彼此交叉(但两条线在3D空间中不是必然相互交叉,但本领域技术人员可以理解,可以选取误差最小化函数来定位最恰当的点以调用这种相互交叉)。借助于泥浆滤液标记83以及也是通过上述相同的时延数据采集推导出的与具有泥浆滤液的天然地层烃(例如Fld-A)替代对应的矢量(例如Sig-I),我们现在在3D空间中具有一条线100,天然地层烃标记点80位于线100上。因此,如果我们仅知道其中一种准确的天然地层烃特性(例如,密度,因为烃的这个参数通常是公知的),那么相应地也就知道了其它的特性。如上面所提到的,图18示出了如何到达泥浆滤液标记(例如,Fld-X),而图19示出了如何到达天然地层烃标记(例如Fld-A)。也就是说,图18-19示出了在示例的地层水盐度可变的情况下如何到达真实的x-组分替代标记,在那里被驱替流体由三种流体:天然地层烃(Fld-A)、原生地层水(Fld-B)以及注入水(Fld-C)的混合物构成。
一旦以及均被可靠地估计,那么可以按照下面的公式计算被驱替流体真实成分和已经侵入地层的泥浆滤液的体积
现在对水盐度可变的地下地层的一种应用进行讨论,这种应用通常产生于用于维持储层压力并且维持烃生产的水注入操作,其中注入的水盐度明显不同于原始的地层水(即原生水)盐度,并且这两种水以不同比例在储层中进行混合将产生不同的水盐度。使用上面描述的方法,我们现在示出如何鉴别和/或指定与原生地层水,注入水,以及天然地层烃对应的不同的流体x-组分替代标记,以及接下来如何连续地将来自泥浆滤液侵入的测井测量结果差解译(沿着井)为不同比率的原生地层水、注入水以及未清理的烃的混合物。产生的流体比率对照另一种现有的技术进行检测和验证,也就是同时使用电阻率和Σ测量结果来解决在钻孔中存在的水盐度和水体积。
作为对比,本方法聚焦于研究被泥浆滤液驱替的流体混合物的成分(即,流动的是什么),而电阻率和Σ技术聚焦于存在于钻孔内的水(并且不是必须被驱替)。此外,本方法采用线性混合定律使用测量结果,而电阻率和Σ技术使用非线性电阻率混合定律,这额外需要电阻率方程参数,例如所谓的Archie’s“M”和“N”参数的使用和/或调节。此外,本方法不使用任何基质参数,因为当取得钻道次和刮道次的差值时基质对于输入的贡献被抵消掉了,而电阻率和Σ技术需要考虑粘土,等,体积校正并且使用合适的基质Σ。
此外,本方法使用两种道次(例如,钻道次和刮道次),而电阻率和Σ技术基于单个道次。还有,本方法在被驱替流体和泥浆滤液之间具有反差时或者当被驱替流体的特性具有差别时能够分辨,而电阻率和Σ技术在水盐度低时无法分辨。此外,本方法中讨论的x-组分替代标记可以与用于钻井的钻井泥浆协作实现井到井的改变,或者可能不存在或难以鉴别,例如当所有的可移动烃已经被刮走,妨碍了对天然地层油标记的确定。然而,在本方法中,要素分析和/或其它统计分析技术可以直接提取新的标记,无论在钻井泥浆系统中是否发生改变。应该注意的是,使用本方法得到的结果与来自电阻率和Σ技术的结果被确定且对照来自流体样本分析的结果。
一种基于上面描述的方法的示例性解释工作流程如下:
1.获取钻道次;
2.获取刮道次;
3.由钻道次计算地层参数,例如根据一个例子的与流体无关的表观孔隙率;
4.由刮道次计算相同的地层参数;
5.将由钻道次和刮道次计算的相同参数进行对比,用于使刮道次与钻道次深度匹配;
6.在执行完上面的深度匹配操作之后,由钻道次和刮道次重新计算相同的地层参数,以便令人满意地确定钻道次和刮道次处于相对彼此的深度处;
7.计算上面描述的基质校正真实孔隙率(与流体无关的表观孔隙率不同);
8.在输入上执行垂直分辨率匹配,执行电阻率和Σ测井测量的同时反演需要所述输入(输入为电阻率、Σ以及真实孔隙率);
9.执行电阻率和Σ测井测量结果的同时反演(符合要求的初步结果);
10.使用该初步结果鉴别泥浆滤液最可能仅仅驱替原生地层水、仅仅驱替注入水、或者仅仅驱替天然地层烃的带;
11.对有效一致测井测量输入进行平均,在本例子中用于在滑移窗口(例如,10英尺窗口,即在取样率为每英尺2个数据点的情况下在21个数据点上)执行本发明的方法,以求出统计噪声的平均数并且在将钻道次与刮道次彼此相减之前进一步减少钻道次和刮道次之间的任何残余深度失配的影响,而且减弱任何残余测量结果轴向分辨率失配。在本示例实施例中,测井测量输入为表观密度孔隙率、表观中子孔隙率以及表观Σ孔隙率。
12.将钻道次与刮道次彼此相减;
13.根据在步骤10中确定的“带”,划分得到的差分数据集,和/或使用要素分析和/或其它统计分析技术来指定与原生地层水、注入水以及天然地层烃对应的各个流体替代标记。
14.沿着井连续地将测井测量结果差解译为不同比率的原生地层水、注入水、以及未清理烃的混合物;
15.减少在步骤11中提到的10英尺的平均间隔以改善输出结果的垂直分辨率同时监视改进的垂直分辨率和增加的统计噪声之间的此消彼长;
16.根据需要,将该方法的结果与电阻率和Σ测井测量结果的同时反演的结果进行比较,同时记住前者注重于研究被泥浆滤液(即,可移动流体)驱替的流体混合物的成分,而后者注重于研究原始地质水与烃(即占据了整个钻孔空间)。
总体来说,测试结果与电阻率和Σ技术相比是有利的,因为计算的水盐度图达成一致。还可以观察到,被驱替流体成分似乎仅表示出主要的“二元系统”。也就是,被驱替流体成分仅仅是原生水和注入水的混合物,或者仅仅是注入水和天然地层油的混合物,或者仅仅是天然地层油+原生水的混合物。
享有前述描述和相关附图的益处的本领域技术人员能够想到许多修改和其它的实施例。因此,可以理解的是各种修改和实施例都被包含在所附的权利要求的范围内。
Claims (15)
1.一种用于确定地质地层内的流体的体积数据的方法,所述地质地层中具有井眼,该方法包括:
基于在相应不同的第一和第二时刻从井眼的测量收集地质地层的第一和第二数据集瞬像,并且使井眼在第一和第二时刻之间经受用于驱替与井眼邻近的地质地层中的流体的流体注入;
基于第一和第二数据集瞬像生成差分数据集;
标准化差分数据集,以生成标准化的差分数据集;
基于标准化的差分数据集确定限定出几何形状且与相应的不同的被驱替流体标记对应的顶点;
确定通过代表具有已知的第一特性的第一被驱替流体的第一点且沿着相应的第一矢量指向的第一条线;
确定通过代表具有已知的第二特性的第二被驱替流体的第二点且沿着相应的第二矢量指向的第二条线;
基于第一条线和第二条线的交点确定与注入流体的特性对应的注入流体点;
确定通过注入流体点且沿着与具有至少一种未知特性的另一被驱替流体对应的另一矢量指向的另一条线;
基于其它被驱替流体的至少一种已知特性确定沿着其它线的第三点;并且
基于差分数据集、第一点、第二点以及第三点确定被驱替流体的体积成分。
2.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,收集第一和第二数据集瞬像包括:分别在钻道次和刮道次过程中使用随钻测井工具收集第一和第二数据集。
3.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,第一和第二数据集瞬像包括伽马射线测量数据、中子测量数据、密度测量数据以及热中子俘获截面数据中的至少一种。
4.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,标准化包括标准化来自差分数据集的数据点,以与球的表面相符。
5.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,标准化包括标准化来自差分数据集的数据点,以与二维平面的表面相符。
6.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,已知的第一和第二特性中的至少一种包括盐度水平。
7.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,具有所述至少一种未知特性的其它被驱替流体包括烃流体。
8.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,其中,第一被驱替流体包括原生水。
9.如权利要求1所述的用于确定地质地层内的流体的体积数据的方法,进一步包括基于被驱替流体的确定的体积成分确定渗透率、相对流体渗透率以及分流量中的至少一种。
10.一种测井系统,包括:
测井工具,用于在相应不同的第一和第二时刻从地质地层内的井眼收集地质地层的第一和第二数据集瞬像,并且使井眼在第一和第二时刻之间经受用于驱替与井眼邻近的地质地层中的流体的流体注入;以及
处理器,用于:
基于在相应不同的第一和第二时刻从井眼的测量收集地质地层的第一和第二数据集瞬像,其中使井眼在第一和第二时刻之间经受用于驱替与井眼邻近的地质地层中的流体的流体注入;
基于第一和第二数据集瞬像生成差分数据集;
标准化差分数据集,以生成标准化的差分数据集;
基于标准化的差分数据集确定限定出几何形状且与相应的不同的被驱替流体标记对应的顶点;
确定通过代表具有已知的第一特性的第一被驱替流体的第一点且沿着相应的第一矢量指向的第一条线;
确定通过代表具有已知的第二特性的第二被驱替流体的第二点且沿着相应的第二矢量指向的第二条线;
基于第一条线和第二条线的交点确定与注入流体的特性对应的注入流体点;
确定通过注入流体点且沿着与具有至少一种未知特性的另一被驱替流体对应的另一矢量指向的另一条线;
基于其它被驱替流体的至少一种已知特性确定沿着其它线的第三点;并且
基于差分数据集、第一点、第二点以及第三点确定被驱替流体的体积成分。
11.如权利要求10所述的测井系统,其中,所述测井工具包括用于分别在钻道次和刮道次过程中收集第一和第二数据集瞬像的随钻测井工具。
12.如权利要求10所述的测井系统,其中,第一和第二数据集瞬像包括伽马射线测量数据、中子测量数据、密度测量数据以及热中子俘获截面数据中的至少一种。
13.如权利要求10所述的测井系统,其中,所述处理器标准化来自差分数据集的数据点,以与球的表面相符。
14.如权利要求10所述的测井系统,其中,所述处理器标准化来自差分数据集的数据点,以与二维平面的表面相符。
15.如权利要求10所述的测井系统,其中,已知的第一和第二特性中的至少一种包括盐度水平。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261620750P | 2012-04-05 | 2012-04-05 | |
US61/620,750 | 2012-04-05 | ||
US13/837,409 US10385677B2 (en) | 2012-04-05 | 2013-03-15 | Formation volumetric evaluation using normalized differential data |
US13/837,409 | 2013-03-15 | ||
PCT/US2013/035296 WO2013152208A1 (en) | 2012-04-05 | 2013-04-04 | Formation volumetric evaluation using normalized differential data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104350233A CN104350233A (zh) | 2015-02-11 |
CN104350233B true CN104350233B (zh) | 2018-06-08 |
Family
ID=49301052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380029177.8A Active CN104350233B (zh) | 2012-04-05 | 2013-04-04 | 使用标准化的差分数据的地层体积评价 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10385677B2 (zh) |
EP (1) | EP2834682B1 (zh) |
CN (1) | CN104350233B (zh) |
CA (1) | CA2869682A1 (zh) |
MX (1) | MX353195B (zh) |
WO (1) | WO2013152208A1 (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9945181B2 (en) * | 2012-08-31 | 2018-04-17 | Halliburton Energy Services, Inc. | System and method for detecting drilling events using an opto-analytical device |
FR3007165B1 (fr) * | 2013-06-13 | 2016-10-28 | Ifp Energies Now | Methode pour optimiser l’exploitation d’un gisement de fluide par prise en compte d’un terme d’echange geologique et transitoire entre blocs matriciels et fractures |
GB2535343B (en) | 2013-09-16 | 2020-03-25 | Schlumberger Holdings | Formation evaluation using stochastic analysis of log data |
US20160298427A1 (en) | 2015-04-09 | 2016-10-13 | Schlumberger Technology Corporation | Continuum sedimentary basin modeling using particle dynamics simulations |
EP3181806A1 (en) * | 2015-12-17 | 2017-06-21 | Services Pétroliers Schlumberger | Method and system for petrophysical quantity estimation from dielectric and neutron capture cross section measurements |
US10208582B2 (en) * | 2016-08-24 | 2019-02-19 | Saudi Arabian Oil Company | Formation water salinity from borehole measurements |
FR3078739B1 (fr) * | 2018-03-09 | 2020-03-27 | Soletanche Freyssinet | Machine de forage comportant un dispositif de connexion pour un dispositif de mesure de verticalite |
CN111997583A (zh) * | 2020-09-01 | 2020-11-27 | 中国石油集团渤海钻探工程有限公司 | 一种基于模糊神经网络的旋转导向电机调速方法 |
US11719094B2 (en) * | 2020-10-23 | 2023-08-08 | Saudi Arabian Oil Company | Reservoir characterization using rock geochemistry for lithostratigraphic interpretation of a subterranean formation |
CN113685167B (zh) * | 2021-08-04 | 2023-06-02 | 中国石油大学(北京) | 核测井方法、核测井装置和计算机可读存储介质 |
US20230313672A1 (en) * | 2022-03-29 | 2023-10-05 | Halliburton Energy Services, Inc. | Fluid Monitoring In Oil And Gas Wells Using Ultra-Deep Azimuthal Electromagnetic Logging While Drilling Tools |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6216532B1 (en) * | 1996-11-29 | 2001-04-17 | Schlumberger Technology Corporation | Gas flow rate measurement |
CN1896458A (zh) * | 2005-01-11 | 2007-01-17 | 施蓝姆伯格海外股份有限公司 | 导出井下流体的流体性质及其不确定度的系统和方法 |
US7453766B1 (en) * | 2003-11-25 | 2008-11-18 | Michael John Padgett | Method for deriving 3D output volumes using summation along flat spot dip vectors |
WO2011086145A1 (en) * | 2010-01-14 | 2011-07-21 | Services Petroliers Schlumberger | Corrected porosity measurements of underground formations |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3993903A (en) | 1975-03-31 | 1976-11-23 | Chevron Research Company | Low-cost but accurate radioactive logging for determining gas saturation in a reservoir |
US4459480A (en) | 1981-12-04 | 1984-07-10 | Mobil Oil Corporation | Use of pulsed neutron logging to evaluate perforation washing |
DE69223589T2 (de) | 1991-10-22 | 1998-12-10 | Halliburton Energy Services, Inc., Houston, Tex. | Verfahren zum Bohrlochmessen während des Bohrens |
US6585044B2 (en) * | 2000-09-20 | 2003-07-01 | Halliburton Energy Services, Inc. | Method, system and tool for reservoir evaluation and well testing during drilling operations |
US7032661B2 (en) * | 2001-07-20 | 2006-04-25 | Baker Hughes Incorporated | Method and apparatus for combined NMR and formation testing for assessing relative permeability with formation testing and nuclear magnetic resonance testing |
ATE331870T1 (de) | 2002-12-31 | 2006-07-15 | Schlumberger Services Petrol | Verfahren und vorrichtung zur zeitversetzen analyse von ursache und wirkung |
US8023690B2 (en) | 2005-02-04 | 2011-09-20 | Baker Hughes Incorporated | Apparatus and method for imaging fluids downhole |
US7480205B2 (en) | 2005-04-20 | 2009-01-20 | Landmark Graphics Corporation | 3D fast fault restoration |
US8794320B2 (en) * | 2006-03-27 | 2014-08-05 | Shell Oil Company | Water injection systems and methods |
CA2659020C (en) | 2006-06-21 | 2018-08-07 | Terraspark Geosciences, L.P. | Extraction of depositional systems |
US7555390B2 (en) | 2007-03-01 | 2009-06-30 | Schlumberger Technology Corporation | Petrophysical interpretation of multipass array resistivity data obtained while drilling |
US7791017B2 (en) | 2007-07-23 | 2010-09-07 | Schlumberger Technology Corporation | Method to simultaneously determine pore hydrocarbon density and water saturation from pulsed neutron measurements |
US8593140B2 (en) * | 2007-11-02 | 2013-11-26 | Schlumberger Technology Corporation | Formation testing and evaluation using localized injection |
US8005618B2 (en) | 2008-01-09 | 2011-08-23 | Schlumberger Technology Corporation | Logging while drilling system |
US8413726B2 (en) * | 2008-02-04 | 2013-04-09 | Marathon Oil Company | Apparatus, assembly and process for injecting fluid into a subterranean well |
US8214151B2 (en) * | 2008-02-20 | 2012-07-03 | Carbo Ceramics Inc. | Methods of identifying high neutron capture cross section doped proppant in induced subterranean formation fractures |
CA2721376C (en) * | 2008-04-15 | 2016-12-13 | Schlumberger Canada Limited | Formation treatment evaluation |
WO2009142873A1 (en) * | 2008-05-22 | 2009-11-26 | Schlumberger Canada Limited | Downhole measurement of formation characteristics while drilling |
WO2009140738A1 (en) * | 2008-05-23 | 2009-11-26 | The Australian National University | Image data processing |
WO2010009411A2 (en) * | 2008-07-17 | 2010-01-21 | Schlumberger Canada Limited | Hydrocarbon determination in presence of electron and chemical ionization |
EP2261459A1 (en) * | 2009-06-03 | 2010-12-15 | BP Exploration Operating Company Limited | Method and system for configuring crude oil displacement system |
CN102639812A (zh) * | 2009-09-11 | 2012-08-15 | C12能源公司 | 基于流体注射的地下储层分析 |
US8502135B2 (en) | 2010-03-25 | 2013-08-06 | Baker Hughes Incorporated | Method for through-casing 3-phase saturation determination |
CA2734170C (en) * | 2011-03-15 | 2013-09-24 | Exxonmobil Upstream Research Company | Method of injecting solvent into an underground reservoir to aid recovery of hydrocarbons |
-
2013
- 2013-03-15 US US13/837,409 patent/US10385677B2/en active Active
- 2013-04-04 CN CN201380029177.8A patent/CN104350233B/zh active Active
- 2013-04-04 CA CA2869682A patent/CA2869682A1/en not_active Abandoned
- 2013-04-04 EP EP13772513.1A patent/EP2834682B1/en not_active Not-in-force
- 2013-04-04 MX MX2014012042A patent/MX353195B/es active IP Right Grant
- 2013-04-04 WO PCT/US2013/035296 patent/WO2013152208A1/en active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6216532B1 (en) * | 1996-11-29 | 2001-04-17 | Schlumberger Technology Corporation | Gas flow rate measurement |
US7453766B1 (en) * | 2003-11-25 | 2008-11-18 | Michael John Padgett | Method for deriving 3D output volumes using summation along flat spot dip vectors |
CN1896458A (zh) * | 2005-01-11 | 2007-01-17 | 施蓝姆伯格海外股份有限公司 | 导出井下流体的流体性质及其不确定度的系统和方法 |
WO2011086145A1 (en) * | 2010-01-14 | 2011-07-21 | Services Petroliers Schlumberger | Corrected porosity measurements of underground formations |
Non-Patent Citations (1)
Title |
---|
电阻率时间推移测井解释;陈福煊;《天然气工业》;19960131(第1期);第25-28页 * |
Also Published As
Publication number | Publication date |
---|---|
US10385677B2 (en) | 2019-08-20 |
WO2013152208A1 (en) | 2013-10-10 |
CN104350233A (zh) | 2015-02-11 |
MX2014012042A (es) | 2015-01-16 |
EP2834682A4 (en) | 2015-11-04 |
MX353195B (es) | 2018-01-05 |
CA2869682A1 (en) | 2013-10-10 |
US20130338926A1 (en) | 2013-12-19 |
EP2834682A1 (en) | 2015-02-11 |
EP2834682B1 (en) | 2016-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104350233B (zh) | 使用标准化的差分数据的地层体积评价 | |
CN104412130A (zh) | 使用标准化的差分数据的地层成分评价 | |
US8965703B2 (en) | Applications based on fluid properties measured downhole | |
US6844729B2 (en) | Method of using nuclear spectroscopy measurements acquired while drilling | |
US10428642B2 (en) | Transposition of logs onto horizontal wells | |
EP2798376B1 (en) | In-situ characterization of formation constituents | |
WO2013074593A1 (en) | Borehole imaging and formation evaluation while drilling | |
WO2021179288A1 (en) | Surface logging with cuttings-based rock petrophysics analysis | |
US10571600B2 (en) | Determination of formation properties using graphical methods | |
US8441269B2 (en) | Determining formation properties while drilling | |
CA2851779C (en) | A method of and apparatus for calculating ucs and ccs | |
WO2014172563A1 (en) | Petrophysical inversions systems and methods field | |
Massiot et al. | Discussion between a reservoir engineer and a geologist: permeability identification from completion test data and borehole image logs integration | |
Hursan et al. | Utilizing NMR Workflow to Optimize Power Water Injector Placement in the Presence of Tar Barriers | |
Kuptsov et al. | Technology update: Magnetic resonance-while-drilling system improves understanding of complex reservoirs | |
Leveridge | New resistivity-logging tool helps resolve problems of anisotropy, shoulder-bed effects | |
Massiot et al. | Discussion between a reservoir engineer and a geologist | |
Denney | Estimation of Variable Fluid-Mixture Density With 4D-NMR Logging | |
Pereira et al. | Estimation of permeability and permeability anisotropy in horizontal wells through numerical simulation of mud filtrate invasion | |
Maeso et al. | Concentric Saturation Maps for Carbonate Heterogeneity at the Wellbore Scale | |
Walker | ADVANCES IN WIRELINE DATA ACQUISITION IN SYNTHETIC-BASED MUD SYSTEMS AND STRATEGIES FOR IMPROVED FORMATION EVALUATION |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: K*Gesella Inventor after: Jain Vikrant Inventor after: A, P, hope Inventor before: K*Gesella Inventor before: Jain Vikrant Inventor before: HIBLER PATRICK A. |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |