CN107810431B - 校准地震检波器阵列的深度的方法 - Google Patents

校准地震检波器阵列的深度的方法 Download PDF

Info

Publication number
CN107810431B
CN107810431B CN201680036323.3A CN201680036323A CN107810431B CN 107810431 B CN107810431 B CN 107810431B CN 201680036323 A CN201680036323 A CN 201680036323A CN 107810431 B CN107810431 B CN 107810431B
Authority
CN
China
Prior art keywords
depth
geophone
channel
stratum
das
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
Application number
CN201680036323.3A
Other languages
English (en)
Other versions
CN107810431A (zh
Inventor
J·C·霍尔曼
A·A·马特瓦
P·M·兹沃吉斯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shell Internationale Research Maatschappij BV
Original Assignee
Shell Internationale Research Maatschappij BV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shell Internationale Research Maatschappij BV filed Critical Shell Internationale Research Maatschappij BV
Publication of CN107810431A publication Critical patent/CN107810431A/zh
Application granted granted Critical
Publication of CN107810431B publication Critical patent/CN107810431B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/181Geophones
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/168Deployment of receiver elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/162Details
    • G01V1/164Circuits therefore
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/22Transmitting seismic signals to recording or processing apparatus
    • G01V1/226Optoseismic systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种地震检波器阵列具有多个地震检波器通道,每一地震检波器通道耦合到地层中的局部环境。确定每一地震检波器通道的地层物质相关响应,且使所述响应与所述对应的地震检波器通道的假设深度相关联。将作为所述假设深度的函数的所述地层物质相关响应与作为沿着钻孔的深度的函数的所述地层的至少一个岩石物理参数的独立深度测井曲线进行比较。基于所述比较,确定所述假设深度和独立深度测井曲线中的深度之间的一组时滞,所述一组时滞提供地层物质相关响应和地层的至少一个岩石物理参数的独立深度测井曲线之间的最佳相关性。因此,每一地震检波器通道的所述假设深度能够与所述独立深度测井曲线中的对应深度比对。

Description

校准地震检波器阵列的深度的方法
技术领域
本发明涉及一种校准地震检波器阵列中的多个地震检波器通道的深 度的方法,所述地震检波器阵列安置于地层中的钻孔中。
背景技术
已经研发出采用安置于地层中的钻孔中的地震检波器阵列的各种地 震勘探技术。实例包含层析成像技术(例如包含井间地震层析成像)和垂 直地震剖面法(VSP)。
分布式声学传感(DAS)是一种出于地震数据获取目的而在地层中的 钻孔中提供此类地震检波器阵列的有效新颖技术。在Albena Mateeva等人 的文章《利用垂直地震剖面法的油藏监测的分布式声学传感(Distributed acoustic sensing for reservoirmonitoring with vertical seismic profiling)》中 描述了这一技术,所述文章出现在《地球物理勘探(Geophysical Prospecting)》(2014)第62卷第679到692页中。从概念上讲,DAS测 量很简单。DAS询问器单元沿着安置于钻井孔中的光纤发送激光脉冲,并 测量后散射光的信号。可基于沿着光纤的激光脉冲的飞行时间,将光纤细 分到DAS检波器通道中(举例来说,对应于VSP检波器级)。然而,识别 DAS检波器通道相对于地质的精确物理深度并非无足轻重,且其需要相对 于钻孔中的深度对光学深度进行一定校准。
类似地,可能需要对其它类型的地震检波器通道的地震检波器阵列, 例如地音探测仪阵列,进行深度校准。
在US 2013/0279841中描述一种用于确定光纤通道的位置的方法。在 此方法中,一个或多个光纤通道的位置通过以下操作确定:
a)在邻近于通道中的至少一个的已知位置处布置电导体和磁性源;
b)通过电导体传输电流,由此通过磁性源附近的洛伦兹(Lorenz)力 使电导体变形;
c)传送电导体的变形以使邻近通道变形;
d)通过光纤电缆传输光脉冲,并使用通过变形通道反射回来的光脉 冲的变化形式和磁性源的已知位置来确定变形通道的位置。
这一已知的用于确定光纤通道的位置的方法的缺点在于,它需要在钻 井孔中设有额外设备来使在已知深度处的光纤局部变形,从而能够使用通 过变形通道反射回来的光脉冲的变化形式和已知深度来确定变形通道相 对于外部参考系(例如地质)的深度。
发明内容
根据本发明的第一方面,提供一种校准地震检波器阵列中的多个地震 检波器通道的深度的方法,所述地震检波器阵列安置于地层中的钻孔中, 所述方法包括:
-选择安置于地层中的钻孔中的地震检波器阵列,其中所述地震检波 器阵列包括多个地震检波器通道,由此每一地震检波器通道局部耦合到在 邻近于地震检波器通道处存在的地层;
-假设所述多个地震检波器通道中的每一地震检波器通道具有某一 深度;
-确定由传播通过邻近于每一相应的地震检波器通道位置的地层的 地震波诱发的每一地震检波器通道的地层物质相关响应;
-提供作为沿着钻孔的深度的函数的地层的至少一个岩石物理参数 的独立深度测井曲线;
-针对数个地震检波器通道,确定所述假设深度和独立深度测井曲线 中的深度之间的一组时滞,所述一组时滞提供地震检波器通道的地层物质 相关响应和地层的至少一个岩石物理参数的独立深度测井曲线之间的最 佳相关性;以及
-针对所述多个地震检波器通道输出一组经修正深度或修正值,以比 对所述多个地震检波器通道中的每一地震检波器通道的假设深度与独立 深度测井曲线中的对应深度。
附图说明
图1示出了地层中的钻孔的示意性纵截面;
图2示出了DAS缆线的示意图;
图3示出了若干个独立深度测井曲线(3a=γ射线测井曲线;3b=声波 测井曲线;3c=密度测井曲线)和使用DAS缆线测量出的零偏移距VSP (3d);
图4示出了γ射线测井曲线(4a)、声波测井曲线(4b)、分块声波测 井曲线(4c)和DAS通道一致性标量(4d),它们全都是从图3中D1和 D3之间的数据导出;以及
图5示出了γ射线测井曲线(5a)、密度测井曲线(5b)、声波测井曲 线(5c)、分块声波测井曲线(5d)、声学阻抗测井曲线(5e)、分块声学 阻抗测井曲线(5f)和DAS通道一致性标量(5g),它们全都是从图3中 D2和D3之间的数据导出。
这些图不按比例。在不同图中使用的相同附图标号指代类似组件。
具体实施方式
本发明将在下文中仅借助于实例并参考非限制性图式进一步进行说 明。
所属领域的技术人员将容易了解,虽然参考特征和量度的一个或多个 特定组合来说明本发明,但这些特征和量度中有许多在功能上独立于其它 特征和量度,使得其可以同等或类似地独立应用于其它实施例或组合中。
在说明书和权利要求中,术语“地层物质相关响应”是一种受地层的形成物质的弹性特性影响的地震检波器通道中的响应,所述检波器通道耦合到所述地层的形成物质。这是邻近于地震检波器通道的位置的地层的形成物质。适当地,地层物质相关响应可以通道一致性振幅标量形式表达,如可使用本领域中已知的通道一致性标量推导技术来确定。如所属领域的技术人员所理解,通过去除可实现真实响应的其它情况的作用或影响,例如震源特征的差异或震源位置的差异等,通道一致性标量是表征特定通道 的响应的单个数值,如在其耦合到地层的局部环境内所配置。所属领域的 技术人员熟知地震检测中的一致性概念,且所述一致性概念在(举例来说) 具有里程碑意义的M T Taner和F Koehler发表于《地球物理(Geophysics)》 第46卷第1期(1981)第17到22页中的题为《地表一致性修正(Surface consistent corrections)》的论文中得以充分解释。如本说明书和权利要求中所介绍的检波器一致性和/或通道一致性被视为是表面一致性的物理类似 物,其适用于检波器通道被安置在地层中的钻孔中而不是地表上的情况。
此外,在本公开的上下文中,除非另外指定,否则术语“深度”大体 上理解为意指测量深度(MD)。MD是本领域中用于表示沿着钻孔的轨道 或路径的测量长度的术语。此测量结果不同于除了垂直钻孔之外的所有钻 孔的真实垂直深度(TVD)。然而,有时候,将MD转换成TVD很方便, 反之亦然。此类转换在钻孔轨道已知的情况下的可能的。此类转换可便于(举例来说)使地震通道深度与地震速度模型相关联。哪一点被定义为零 深度并不重要。为方便起见,建议将地表处的钻孔的顶端定义为零深度。
在最广泛意义上,术语“数个”是指一个或多个。在特定实施例中, 术语“数个”可意指多个或两个或更多个。
本文中提出一种校准地震检波器阵列中的多个地震检波器通道的深 度的方法,所述地震检波器阵列安置于地层中的钻孔中。地震检波器阵列 包括多个地震检波器通道,其可形成一串互连的地震检波器通道。假设每 一地震检波器通道具有某一深度。每一地震检波器通道局部耦合到在邻近 于地震检波器通道位置处存在的地层。在本公开的上下文中,邻近于地震 检波器通道是指相比于钻孔的纵向方向,地层仅仅在地震检波器通道的位置的横向方向上完全在钻孔外部。确定每一地震检波器通道的地层物质相 关响应。每一地震检波器通道的此类地层物质相关响应适当地由每一地震 检波器通道的单个数值(其可为标量)表示,所述数值优选的是信号振幅 的值。
因此,每一地震检波器通道都将获得一些地层物质相关响应和假设深 度,并且因此,地层物质相关响应可被视为假设深度的函数。作为假设深 度的函数的地层物质相关响应与作为沿着钻孔的深度的函数的地层的至 少一个岩石物理参数的独立深度测井曲线比较。基于所述比较,确定所述 假设深度和独立深度测井曲线中的深度之间的一组时滞,所述一组时滞提 供地层物质相关响应和地层的至少一个岩石物理参数的独立深度测井曲线之间的最佳相关性。
因此,可针对多个地震检波器通道输出一组经修正深度或一组修正 值,以比对所述多个地震检波器通道中的每一地震检波器通道的假设深度 与独立深度测井曲线中的对应深度。通过使地震检波器通道的地层物质相 关响应与独立深度测井曲线相关,基本上整个检波器阵列或至少大部分检 波器阵列可与岩石物理形成特征比对,而不必依赖于例如敲击测试点的若 干个固定点之间的内插和阵列末端比对(end-of-array alignment)。另一优 点在于,这可以在不需要在钻孔中安装额外设备的情况下完成。此外,它 可回顾性地适用于已经使用地震检波器阵列获取的现有地震数据集。
再者,所述方法的优点在于,多个地震检波器通道中的每一地震检波 器通道的假设深度与独立深度测井曲线中的对应深度比对。这可能与真实 深度不同,但是出于多个目的,更为重要的是知道地震检波器通道相对于 地层中的地质特征的相对位置。
所述方法还可适用于时延模式,以最大化地震年份之间的重复性。所 述方法可揭示地震检波器通道随着时间推移的位移。在地震检波器通道是 DAS检波器通道的情况下,所述方法还可用于防范DAS光纤内的潜在改 变,例如折射率随着时间推移的改变,这可在温度改变的影响下出现。此 类温度改变可出现(举例来说)在可(举例来说)涉及蒸汽射流的提高石 油采收率(EOR)操作的过程中。术语“DAS光纤”关于光学耦合到询问 单元的任何类型的光纤,所述询问单元能够将连续的分布式信号集聚到通 道信号中。术语“DAS缆线”关于封装在单个缆线中的一个或多个DAS 光纤,由此DAS光纤由至少一种保护材料包埋和/或围绕。
地层的物质影响地震检波器对入射地震波的响应。如上所陈述,术语 “地层物质相关响应”是一种受地层的物质的局部弹性特性影响的地震检 波器通道中的响应,所述检波器通道耦合到所述地层的物质。利用粒子运 动通过物质传播的地震波可通过运动传感器或应变感应器(例如,地音探 测仪或DAS光纤)接收。响应于地震波的物质变形受物质的弹性张量控 制,且因此与检波器通道耦合的物质的局部弹性特性影响所述响应。所述 物质的弹性特性还影响对通过物质传播的地震波(如压力波)的响应。这 是可以理解的,因为压力脉冲在恒定粒子速度下通过所述物质的声学阻抗 显现。
适当地,独立深度测井曲线中地层的岩石物理参数是取决于地层的弹 性特性的参数。
地层物质相关响应和独立测井曲线的至少一个岩石物理参数优选地 都对地震检波器通道周围的地层的弹性特性敏感。实例包含声波测井曲 线、密度测井曲线和声学阻抗测井曲线。可通过组合声波测井曲线和密度 测井曲线来建构声学阻抗测井曲线。在一些情况下,举例来说,在某些碎 屑环境中,包括砂泥岩地层,γ射线测井曲线可指示声学阻抗,从而同样 适用于此目的。
本文中所描述的方法可用于任何类型的地震检波器阵列,包括地音探 测仪和DAS光纤的阵列。然而,出于简洁性起见,将采用DAS光纤作为 实例来描述说明书的其余部分。DAS光纤的主要优点在于,覆盖范围通常 大于地音探测仪的覆盖范围,和/或地震检波器通道时滞可更低(即,每单 位长度具有更多数目的地震检波器通道)。这两个因素使得更易于使地层 物质相关检波器响应与独立测井曲线相关。所属领域的技术人员能够将相 同原理应用到其它类型的地震检波器阵列。
现参考图1,示出了地层50中的钻孔10的示意性纵截面。钻孔10 穿过数个地质层(由区域51到55示意性地指示),每一地质层具有其自 身的特性,包含弹性特性。DAS光纤20安置于钻孔10中。在多个DAS 检波器通道22中细分(DAS)光纤20。为了看起来清楚,图1中仅描绘 出这些DAS检波器通道22中的几个。实际上,数百个DAS检波器通道 很常见,其中检波器通道间距在2m和10m之间。每一地震检波器通道 22局部耦合到在邻近于所讨论的地震检波器通道22的位置处存在的地层 50。这可通过例如将DAS缆线粘固在井筒后面来完成。为免生疑,在本 公开的上下文中,如果在横向于地震检波器通道和在邻近于所讨论的地震检波器通道的位置处存在的地层之间的钻孔的纵向方向的平面内,地震检 波器通道的位置和所述地层之间存在物理接触路径,无论是完全还是通过 可包含粘固剂和/或井管(例如,套管)的中间材料,那么地震检波器通道 局部耦合到所述地层。
询问器单元30配置在地表35上。询问器单元30连接到DAS光纤20, 且布置成将光学激光脉冲传输到DAS光纤20中并检测受沿着光纤的变形 影响的后散射光信号,所述变形可由冲击地震波产生。由于光纤的任何部 分都可进行变形并被询问地震信息,所以这种测量类型被称为分布式声学 传感。任何变形的位置可根据感测到变形的光学激光脉冲的已知飞行时间 来确定。因此,可基于沿着DAS光纤20的光学激光脉冲的飞行时间,将 DAS光纤20细分到DAS检波器通道22中(举例来说,对应于VSP检波 器级)。
适当地,DAS光纤20封装在DAS缆线40中。图2中示出了实例的 示意图。DAS缆线可包括多个DAS光纤,其在同时使用时可增强信噪比。 作为一实例,图2中示出了两个笔直的纵向光纤20a,但是可提供更多个, 例如五个。或者或另外,可在DAS缆线40中配置一个或多个成螺旋形卷 绕的光纤20b。适当地,这些成螺旋形卷绕的光纤20b可卷绕于核心45 周围。所有事物可包埋于保护材料中,并被一个或多个保护性外层覆盖。 作为实例示出一个保护性外层42。不同于笔直的纵向光纤20a,包括成螺 旋形卷绕的光纤20b的缆线对边射p波(在垂直于缆线的纵向方向延伸的 平面内在缆线与所述平面的相交处具有传播分量的p波)敏感。
可假设多个地震检波器通道中的每一地震检波器通道22具有某一深 度。常见过程是确定DAS光纤20的远端位于钻孔中的位置,并接着导出 DAS检波器通道22在所述远端位置之上的假设深度。DAS检波器通道22 的假设深度还可从光学距离(在关于光纤轨道的假设下)导出,如根据经 反射光学脉冲的飞行时间确定。此光学深度不必与沿着钻孔10的测量深 度一致,因为可能会产生偏差,例如相对于钻孔10的光纤的松弛,或关 于DAS光纤20的折射率的错误假设。
图3的部分3d示出了通过针对部署在钻孔中的五个光纤同时使用 DAS而获取的零偏移距VSP道集的地震道。所述五个光纤可配置在单个 DAS缆线中。使用基于TVD尺度的假设深度,在水平轴上绘制响应时间, 并沿着垂直轴绘制DAS检波器通道。为了进行参考,指定三个深度级(D1、 D2、D3)。可在左侧看到(下行)直达波,且由地质界面的反射产生的上行波也可见。值得注意的是,观察到由在深度D2下方的特征的反射产生 的大量上行波。本发明的实施例允许比对假设深度与更有意义的地质深 度。
为此目的,提供表示与地层50的弹性特性相关的参数且作为沿着钻 孔10的深度的函数的独立深度测井曲线。图3a到3c示出了在与图3d相 同的深度尺度上绘制且在同一钻孔中测量的各种独立测井曲线。水平尺度 是线性尺度。图3a是γ射线测井曲线——水平尺度区间在0.0到150个 API单位([1}由美国石油学会定义的在钻孔中测量的自然γ辐射的标准量 度){2]范围内。图3b是p波的声波速度测井曲线——水平尺度区间在2000 到6000m/s范围内。图3c是密度测井曲线——水平尺度区间在1.45到2.95 g/cm3范围内。仅出于参考目的公开水平尺度;显而易见的是,实行本发 明的实施例不需要绝对值。
为了使图3d的地震道与独立测井曲线相关,通道一致性标量(在当 前情况中:DAS通道一致性标量)是从DAS VSP数据导出,其方式类似 于地表地震数据的地表一致性标量推导(同样参考Taner和Koehler的论 文)。通道一致性标量与每一地震检波器通道中的平均信号振幅成反比。
适当地,相对于每一地震检波器通道的参考信号振幅来确定信号振幅 的值。适当地,信号振幅的值对应于多个相干地震事件的均方根(RMS) 平均振幅,例如对应于见于图3d中的上行波的所有信号。适当地,考虑 多个道集,例如在震源检波器偏移距彼此不同的情况下获得的多个道集, 以导出DAS通道一致性标量。变井源距VSP炮点道集是合适的多个道集 的实例。如此导出的通道一致性标量构成地层物质相关响应的适宜量度。
以任一种此类方式,由通过地层50传播的地震波诱发的每一地震检 波器通道22的地层物质相关响应可因此使用从DAS光纤20后散射的光 的信号的实际测量来确定,且其被视为地震检波器通道的假设深度的函 数,地层物质相关响应通过所述假设深度来确定。
此过程的结果在图4中说明。此图中的图4d示出了基于关于深度的 线性尺度的DAS通道一致性标量的曲线,其从来自利用炸药震源的2D变 井源距VSP勘测的地震道导出。所述深度基于DAS检波器通道的假设深 度。(深度绘制为MD。)仅考虑D1和D3之间的数据。为了易于比较,图 3a的γ射线测井曲线中的一段在图4a中再现,且图3b的声波测井曲线中的一段在图4b中再现。
图4c对应于粗化(upscaling)之后的图4b的声波测井曲线。粗化用 以促进独立测井曲线与DAS检波器通道的检波器一致性标量的比较。这 极其适用于独立深度测井曲线的深度取样比DAS检波器通道取样更精细 (独立深度测井曲线中地层的至少一个岩石物理参数的值——例如p波速 度——可用于每一单位长度的数个深度,所述深度的数目大于每一所述单 位长度的DAS检波器通道的数目)的情况。在所呈现的情况中,已经通 过数据分块执行粗化,其基本上用于将测井值的块表示为单个数据点。适 当地,经粗化深度测井曲线针对每一DAS检波器通道具有地层的至少一 个岩石物理参数的单个值。
如通过比较图4c与图4d可见,图4d的DAS通道一致性标量与图4c 的分块声波测井曲线之间存在显著的相关性。还可基于图4b的未分块数 据观察所述相关性,但是类似性在图4c中可以更为显著地观察到。
一旦相关性建立好,就可针对数个地震检波器通道确定所述假设深度 和独立深度测井曲线中的深度之间的一组时滞。这可基于匹配一些分散的 明显特征,例如通过图4c和4d之间的连接线指示的那些特征,或通过使 用计算技术最大化两个数据集之间的相关性来完成。一种最大化所述相关 性的可能计算技术是在两个数据集之间的滑动窗中挑选互相关的最大值。 最终目标是确定提供地层物质相关响应(例如,地震通道一致性标量)和地层的至少一个岩石物理参数的独立深度测井曲线(例如,声波测井曲线) 之间的最佳相关性的一组时滞。接着可针对多个地震检波器通道建立并输 出一组修正值。接着,所述一组修正值可用于比对多个地震检波器通道中 的每一地震检波器通道的假设深度与独立深度测井曲线中的对应深度。
图4a的γ射线测井曲线和图4d的DAS通道一致性标量之间也存在 相关性,但在本实例中,这是一种弱得多的相关性,因为γ射线测井曲线 仅与地层的弹性特性间接相关。γ射线测井曲线通常与给定的沉积环境的 声学阻抗相关,但是所述相关性并非一直存在。因此,为了深度校准这一 目的,通常优选的是声波测井曲线、密度测井曲线或声学阻抗测井曲线, 它们是地层的各种弹性特性的直接量度。但是,γ射线测井曲线也可用于 一些情形。
作为一实例,图5探索了深度D2和D3之间的复杂区域。如图3c中 可见,此深度范围中的大部分都可获得密度测井曲线,从而允许在声学阻 抗等于介质中的声波速度和介质密度的乘积的条件下,导出此区域的声学 阻抗。图5b中示出了密度测井曲线,并且声波测井曲线中的相关区段在 图5c中再现。为了一般参考,图3a的γ射线测井曲线中的对应区段在图 5a中再现。图5d是分块声波测井曲线(同样,p波速度)。图5e对应于 通过使一深度范围中的图5b和5c相乘而获得的声学阻抗,在所述深度范 围中,可获得声波测井曲线以及密度测井曲线两者的数据。图5f表示使用 分块技术进行粗化之后的图5e的数据。图5g示出了DAS通道一致性标 量。即使在此特定深度范围中的这一复杂地质中,基于多个地震检波器通 道中的每一地震检波器通道的假设深度可进行修正以匹配独立深度测井 曲线中的对应深度,图5g的DAS通道一致性标量和图5d的声波测井曲 线以及图5f的声学阻抗测井曲线之间也有可能存在显著的相关性。
本发明的计算实施方案可涉及使作为深度的函数的地层物质相关响 应与作为深度的函数的土体参数的独立测井曲线相关,并确定最大化相关 性的每一地震检波器通道的时滞。然而,明显地,并不总是需要利用独立 测井曲线计算所有地震检波器通道的时滞。例如,参见图4,它可能足以 用于确定数据中选定数目的清晰的对应特征的时滞,并接着基于数目小得 多的已确定时滞,进行内插以产生修正值的输出然而,在更复杂的如(举 例来说)图5中所示的数据集中,可建立大量的时滞以确保找到最佳相关 性。
到目前为止一直解释的方法有利地采用地层物质相关响应和某些深 度测井曲线,具体地说,声波深度测井曲线和声学阻抗深度测井曲线之间 的所观察到的相关性,通过比对所述深度与独立深度测井曲线来校准地震 检波器通道的深度。
相反,所述所观察到的相关性可用于在相当了解地震检波器通道深度 且所述深度精确的情况下,产生声学阻抗测井曲线或声波测井曲线的代 理。这可有力地补充无法获得任何深度测井曲线或仅可获得较短测井曲线 的情况。因此,可从确实具有深度分辨地震数据的钻孔推断出缺失的岩性 信息。
逆向方法可具有多种有用的应用。举例来说,基于从DAS通道因素 导出的伪声学阻抗测井曲线,低频地表地震数据可与井中的人造物质(well synthetics)相关联。这适用于声波和/或密度测井曲线过短(在垂直意义上) 而不能实现低频地表地震数据的常规连井的那些情况。一个钻孔中的地层 物质相关响应,具体地说,与声学阻抗有关的响应的特征和远距钻孔可获 得的独立测井曲线之间可存在对应关系。这是一种连井的新方式,且可用 于限制通常不可获得测井曲线的地下浅层中的凹陷,或用于检测和/或验证 邻近井之间的地层终止处。
另一实例存在于地质模型的垂直深度校准中。尤其是在DAS光纤获 取的情况下,相比于使用标准VSP数据或比使用基于地震波旅行时的地震 反演的结果,可获得更高的分辨率。
所属领域的技术人员将理解,在不脱离所附权利要求书的范围的情况 下,可以多种不同方式实施本发明。

Claims (13)

1.一种校准地震检波器阵列中的多个地震检波器通道的深度的方法,所述地震检波器阵列安置于地层中的钻孔中,所述方法包括:
选择安置于地层中的钻孔中的地震检波器阵列,其中所述地震检波器阵列包括多个地震检波器通道,由此每一地震检波器通道局部耦合到在邻近于所述地震检波器通道位置处存在的所述地层;
假设所述多个地震检波器通道中的每一地震检波器通道具有某一深度;
确定由传播通过邻近于每一相应的地震检波器通道的所述地层的地震波诱发的每一地震检波器通道的地层物质相关响应,其中将每一地震检波器通道的所述地层物质相关响应确定为信号振幅的值;
提供作为沿着钻孔的深度的函数的地层的至少一个岩石物理参数的独立深度测井曲线;
针对数个地震检波器通道,确定所述假设深度和独立深度测井曲线中的深度之间的一组时滞,所述一组时滞提供地震检波器通道的地层物质相关响应和地层的至少一个岩石物理参数的独立深度测井曲线之间的最佳相关性;以及
针对所述多个地震检波器通道输出一组经修正深度或修正值,以比对所述多个地震检波器通道中的每一地震检波器通道的假设深度与所述独立深度测井曲线中的对应深度。
2.根据权利要求1所述的方法,其中所述独立深度测井曲线包括针对每一单位长度的数个深度的所述地层的所述至少一个岩石物理参数的值,所述深度的数目大于每一所述单位长度的地震检波器通道的数目,所述方法包括使所述独立深度粗化,以获得针对每一地震检波器通道具有所述地层的所述至少一个岩石物理参数的单个值的经粗化深度测井曲线。
3.根据权利要求1所述的方法,其中所述独立深度测井曲线中所述地层的至少一个岩石物理参数取决于所述地层的弹性特性。
4.根据权利要求1所述的方法,其中所述独立深度测井曲线是以下中的一个:声波测井曲线、密度测井曲线和声学阻抗测井曲线。
5.根据权利要求1所述的方法,其中所述信号振幅的值对应于多个地震事件的均方根平均振幅的倒数。
6.根据权利要求1所述的方法,其中使用通道一致性标量技术来确定所述信号振幅的值。
7.根据权利要求1所述的方法,其中所述地震检波器阵列中的所述地震检波器通道形成一串地震检波器通道。
8.根据权利要求1所述的方法,其中通过分布式声学传感(DAS)系统形成所述地震检波器阵列,所述分布式声学传感(DAS)系统在多个DAS检波器通道中细分光纤,由此所述地震检波器通道是所述DAS检波器通道。
9.根据权利要求8所述的方法,其中所述光纤封装在缆线中且用作DAS光纤。
10.根据权利要求9所述的方法,其中多个光纤封装在所述缆线中,且所有光纤都用作DAS光纤。
11.根据权利要求9所述的方法,其中在所述缆线中,所述光纤成螺旋形卷绕于核心周围。
12.根据权利要求9到11中任一项所述的方法,其中使用从所述光纤后散射的光的信号的实际测量来确定每一地震检波器通道的所述地层物质相关响应。
13.根据权利要求12所述的方法,其中所述假设深度对应于通过测量所述后散射光的飞行时间并考虑通过地层的光纤轨道所确定的光学距离。
CN201680036323.3A 2015-06-26 2016-06-24 校准地震检波器阵列的深度的方法 Active CN107810431B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP15174021 2015-06-26
EP15174021.4 2015-06-26
PCT/EP2016/064654 WO2016207341A1 (en) 2015-06-26 2016-06-24 Method of calibrating depths of a seismic receiver array

Publications (2)

Publication Number Publication Date
CN107810431A CN107810431A (zh) 2018-03-16
CN107810431B true CN107810431B (zh) 2019-06-11

Family

ID=53502480

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680036323.3A Active CN107810431B (zh) 2015-06-26 2016-06-24 校准地震检波器阵列的深度的方法

Country Status (8)

Country Link
US (1) US10746887B2 (zh)
EP (1) EP3314308B1 (zh)
CN (1) CN107810431B (zh)
AU (1) AU2016282904B2 (zh)
BR (1) BR112017027708B1 (zh)
CA (1) CA2988953C (zh)
MY (1) MY185535A (zh)
WO (1) WO2016207341A1 (zh)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3020007C (en) 2016-04-07 2023-01-31 Bp Exploration Operating Company Limited Detecting downhole events using acoustic frequency domain features
BR112018070565A2 (pt) 2016-04-07 2019-02-12 Bp Exploration Operating Company Limited detecção de eventos de fundo de poço usando características de domínio da frequência acústicas
AU2017355984B2 (en) 2016-10-06 2020-04-09 Shell Internationale Research Maatschappij B.V. Method of borehole time-lapse monitoring using seismic waves
CA3058256C (en) 2017-03-31 2023-09-12 Bp Exploration Operating Company Limited Well and overburden monitoring using distributed acoustic sensors
EA202090528A1 (ru) 2017-08-23 2020-07-10 Бп Эксплорейшн Оперейтинг Компани Лимитед Обнаружение мест скважинных пескопроявлений
CN111771042A (zh) 2017-10-11 2020-10-13 英国石油勘探运作有限公司 使用声学频域特征来检测事件
CN109164482B (zh) * 2018-07-10 2020-07-28 北京昊锐科技有限公司 基于光纤传感器的井下微地震求解方法、装置及系统
US11073629B2 (en) 2018-10-16 2021-07-27 Halliburton Energy Services, Inc. Method to improve DAS channel location accuracy using global inversion
US20210389486A1 (en) 2018-11-29 2021-12-16 Bp Exploration Operating Company Limited DAS Data Processing to Identify Fluid Inflow Locations and Fluid Type
GB201820331D0 (en) 2018-12-13 2019-01-30 Bp Exploration Operating Co Ltd Distributed acoustic sensing autocalibration
CN109991662A (zh) * 2019-05-15 2019-07-09 中油奥博(成都)科技有限公司 浅部地层二维或三维弹性参数测量和计算的装置及方法
BR112021022662A2 (pt) * 2019-06-11 2021-12-28 Halliburton Energy Services Inc Sistema de furo de poço e método de detecção acústica distribuída
WO2021037586A1 (en) * 2019-08-27 2021-03-04 Bp Exploration Operating Company Limited Depth calibration for distributed acoustic sensors
CA3154435C (en) 2019-10-17 2023-03-28 Lytt Limited Inflow detection using dts features
WO2021073741A1 (en) 2019-10-17 2021-04-22 Lytt Limited Fluid inflow characterization using hybrid das/dts measurements
WO2021093974A1 (en) 2019-11-15 2021-05-20 Lytt Limited Systems and methods for draw down improvements across wellbores
US11041382B2 (en) * 2019-11-25 2021-06-22 Halliburton Energy Services, Inc. Vector strain sensor system for a wellbore
WO2021137846A1 (en) 2019-12-30 2021-07-08 Halliburton Energy Services, Inc. Fiber optic cable depth calibration and downhole applications
CA3180595A1 (en) 2020-06-11 2021-12-16 Lytt Limited Systems and methods for subterranean fluid flow characterization
EP4168647A1 (en) 2020-06-18 2023-04-26 Lytt Limited Event model training using in situ data
US20220307895A1 (en) * 2021-03-24 2022-09-29 Schlumberger Technology Corporation True particle velocity wavefield processing in fiber optics - particle motion sensor hybrid array

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4312040A (en) * 1970-09-09 1982-01-19 Schlumberger Limited Well log depth aligning
CN1114511A (zh) * 1993-09-21 1996-01-03 法国石油研究所 采用统计校准技术导出地质特性的地震记录曲线分析方法
US5966672A (en) * 1997-07-28 1999-10-12 Knupp; Daniel F. Visualization technology method
CN102483347A (zh) * 2009-11-13 2012-05-30 秦内蒂克有限公司 分布式感测的改进
CN103534615A (zh) * 2011-03-16 2014-01-22 光学感应器控股有限公司 使用分布式声传感器的地下监控

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE445855T1 (de) * 2002-12-13 2009-10-15 Schlumberger Technology Bv Verfahren und vorrichtung zur verbesserten tiefenanpassung von bohrlochbildern oder probenbildern
US20060180349A1 (en) * 2005-02-16 2006-08-17 Baker Hughes Incorporated Time and depth correction of MWD and wireline measurements using correlation of surface and downhole measurements
US20060256657A1 (en) * 2005-05-11 2006-11-16 Prism Seismic, Inc. Method for improving the time-depth tie of well log data and seismic data
US9234999B2 (en) 2010-12-21 2016-01-12 Shell Oil Company System and method for making distributed measurements using fiber optic cable
AU2011351365B2 (en) 2010-12-31 2015-05-07 Shell Internationale Research Maatschappij B.V. Method and system for determining the location of a fiber optic channel along the length of a fiber optic cable
US9494461B2 (en) 2011-12-15 2016-11-15 Shell Oil Company Detecting broadside acoustic signals with a fiber optical distrubuted acoustic sensing (DAS) assembly
GB2546937B (en) * 2012-11-02 2017-11-29 Silixa Ltd Combining seismic survey and DAS fluid flow data for improved results

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4312040A (en) * 1970-09-09 1982-01-19 Schlumberger Limited Well log depth aligning
CN1114511A (zh) * 1993-09-21 1996-01-03 法国石油研究所 采用统计校准技术导出地质特性的地震记录曲线分析方法
US5966672A (en) * 1997-07-28 1999-10-12 Knupp; Daniel F. Visualization technology method
CN102483347A (zh) * 2009-11-13 2012-05-30 秦内蒂克有限公司 分布式感测的改进
CN103534615A (zh) * 2011-03-16 2014-01-22 光学感应器控股有限公司 使用分布式声传感器的地下监控

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A practical approach to well-seismic data calibration;VICTOR LINARI;《THE LEADING EDGE》;20040801;第23卷(第8期);第774-775页

Also Published As

Publication number Publication date
WO2016207341A1 (en) 2016-12-29
BR112017027708A2 (pt) 2018-09-11
AU2016282904B2 (en) 2019-07-11
EP3314308A1 (en) 2018-05-02
EP3314308B1 (en) 2019-04-10
US20190004195A1 (en) 2019-01-03
US10746887B2 (en) 2020-08-18
AU2016282904A1 (en) 2017-12-14
CA2988953A1 (en) 2016-12-29
CA2988953C (en) 2023-05-23
MY185535A (en) 2021-05-19
BR112017027708B1 (pt) 2021-03-02
CN107810431A (zh) 2018-03-16

Similar Documents

Publication Publication Date Title
CN107810431B (zh) 校准地震检波器阵列的深度的方法
CA3064870C (en) Angular response compensation for das vsp
US9880047B2 (en) Fiber optic distributed vibration sensing with directional sensitivity
Yu et al. Walkaway VSP using multimode optical fibers in a hybrid wireline
EP2686709B1 (en) Subsurface monitoring using distributed acoustic sensors
US9465126B2 (en) Seismic geophysical surveying
Hull et al. Monitoring horizontal well hydraulic stimulations and geomechanical deformation processes in the unconventional shales of the Midland Basin using fiber-based time-lapse VSPs, microseismic, and strain data
CN111239798A (zh) 光纤声波传感井地地震数据联采系统和井驱数据处理方法
CN209911570U (zh) 井中光纤时频电磁和四分量地震数据采集装置
White et al. Toward improved distributed acoustic sensing sensitivity for surface-based reflection seismics: Configuration tests at the Aquistore CO2 storage site
CN109804273B (zh) 使用地震波进行井眼时移监测的方法
CA2961201C (en) Method and system for acquisition of seismic data
US20220179112A1 (en) Detecting and monitoring formation features with an optical fiber
Hiltunen et al. Continuation of down-hole geophysical testing for rock sockets.
Nasralla Distributed Acoustic Sensing (DAS) As a New Tool For Subsurface Imaging
Jiang Borehole seismic imaging: A full waveform inversion approach

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant