CN101568856A - 通过变井源距垂直地震剖面测量进行的盐下反射层析和成像 - Google Patents
通过变井源距垂直地震剖面测量进行的盐下反射层析和成像 Download PDFInfo
- Publication number
- CN101568856A CN101568856A CNA2007800234002A CN200780023400A CN101568856A CN 101568856 A CN101568856 A CN 101568856A CN A2007800234002 A CNA2007800234002 A CN A2007800234002A CN 200780023400 A CN200780023400 A CN 200780023400A CN 101568856 A CN101568856 A CN 101568856A
- Authority
- CN
- China
- Prior art keywords
- area
- travel
- sensor
- stratum
- time
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/42—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators in one well and receivers elsewhere or vice versa
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
使用位于基盐(105)附近的钻孔(101)中的接收器(111a-111d)来进行变井源距VSP测量。来自变井源距VSP的数据(125a-125c)的反射层析反演用于得到地下的速度模型,并可用于使盐下反射成像。
Description
技术领域
本发明涉及提高地震偏移的精确度的地球物理探查方法。特别地,本发明使用变井源距垂直地震剖面测量来确定地下速度并对大地中盐层之下的反射进行成像。
背景技术
在地面地震勘探中,通过地震源传递到大地中的能量从地下地球物理特征反射,并被多个接收器记录。使用可能形成一条线(2维获取)或覆盖一区域(3维获取)的源和接收器配置,将该处理重复多次。使用称为偏移(migration)的步骤对所得到的数据进行处理,以生成反射体的像。
传统的反射波勘探法利用地面源和接收器,来检测由于地下阻抗差而导致的反射。由于源、反射体和接收器之间的长且复杂的传播路径,导致所获得的图像通常在空间精度、分辨率和连贯性方面受到损害。地下的盐层尤其会产生问题。由于盐的压缩波(P波)速度很高(4.48km/s或者14500英尺/秒),因此在盐层的上部和下部存在由于大的速度差而引起的相当大的P波射线弯曲。典型的墨西哥湾(Gulf ofMexico)中的沉积速度可能不大于3km/s。
人们已经采用了许多方法来解决盐下成像的问题。这些方法包括:使用低频率,使用叠前深度偏移(prestack depth migration),使用转换的波,对基盐(base salt)反射进行基准面校正(redatuming),以及地震反演(seismic inversion)。这些方法取得的成功有限。
本发明使用变井源距垂直地震剖面(VSP)测量,通过对反射传播时间进行层析反演,来估计盐下速度。在变井源距VSP测量中,通过响应于与井孔相距多个距离的一或多个地震源的激励使用在钻孔中的多个接收器,来进行测量。然后可以将所估计的速度用于变井源距VSP数据或地面地震数据的偏移。在排水井的钻探(其中已经钻出了能生产或不能生产的初始井)中,这种方法尤其有用。使用本发明的方法,可以对远离初始井的大地的地下成像。
发明内容
本发明的一实施例为一种用于估计地层的特性的方法。所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。该方法包括:在第一区域中布置第一传感器,激活在地面或地面附近的多个源位置处的能量源,以及记录响应于所述源的激活的来自第一传感器的信号。记录信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。在与所述反射相对应的记录信号中拾取传播时间。使用所拾取的传播时间来估计第一区域的下表面的位置、所述界面的位置以及/或者所述第一区域和所述界面之间的速度。第一传感器可以位于界面的下表面附近。第一区域可以包括盐。所拾取的传播时间可以对应于压缩波传播时间。可以在对所拾取的传播时间进行层析反演以给出速度模型的基础上来进行估计。可以进行记录信号的波场分离。所述方法可以还包括识别低速区域。该方法可以还包括使用速度模型偏移记录信号。该方法可以还包括基于所述速度模型进行钻探操作。
本发明的另一实施例为一种用于估计地层的特性的系统,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。所述系统包括被配置为在地面或地面附近的多个源位置处被激活的能量源。第一传感器被配置为生成响应于所述源的激活的信号,其中,所生成的信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。处理器被配置为从与所述反射相对应的记录信号中拾取传播时间,并根据所拾取的传播时间来估计第一区域的下表面的位置、所述界面的位置以及/或者所述第一区域和所述界面之间的速度。第一传感器可以位于第一区域的下表面附近。第一区域可以包括盐。所述处理器可以进一步被配置为拾取与压缩波相对应的传播时间。所述处理器可以被配置为对所拾取的传播时间进行层析反演,以给出速度模型。所述处理器可以进一步被配置为识别地层中的低速区域。所述处理器可以进一步被配置为使用速度模型偏移记录信号。所述系统可以进一步包括被配置为将所述偏移的输出用于进一步的钻探操作的处理器。所述第一传感器可以包括3-C传感器和/或水下地震检波器(hydrophone)。可以使用线缆将第一传感器传送到钻孔中。
本发明的另一实施例为一种与用于估计地层特性的系统一起使用的计算机可读介质,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。所述系统包括被配置为在地面或地面附近的多个源位置处被激活的能量源。所述系统还包括被配置为生成响应于所述源的激活的信号的第一传感器,所生成的信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。所述介质包括使处理器进行以下处理的指令:在与所述反射相对应的记录信号中拾取传播时间,并根据所拾取的传播时间来估计第一区域的下表面的位置、所述界面的位置以及/或者所述第一区域和所述界面之间的速度。所述介质可以包括ROM、EPROM、EAROM、闪存和光盘。
本发明的一实施例为一种用于识别地层中的过压区域的方法。所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。所述方法包括:在第一区域中布置第一传感器,激活在地面或地面附近的多个源位置处的能量源,以及记录响应于所述源的激活、来自第一传感器的信号。记录信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。在与所述反射相对应的记录信号中拾取传播时间。对所拾取的传播时间进行层析反演,以给出速度模型。通过与模型中的低速区域相关联来识别过压区域。第一传感器可以位于界面的下表面附近。第一区域可以包括盐。所拾取的传播时间可以对应于压缩波传播时间。可以在对所拾取的传播时间进行层析反演以给出速度模型的基础上来进行估计。可以进行记录信号的波场分离。可以基于所识别的过压区域来进行钻探操作。
本发明的另一实施例为一种用于识别地层中的过压区域的系统,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。所述系统包括被配置为在地面或地面附近的多个源位置处被激活的能量源。第一传感器被配置为生成响应于所述源的激活的信号,其中,所生成的信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。处理器被配置为从与所述反射相对应的记录信号中拾取传播时间,并对所拾取的传播时间进行层析反演以给出速度模型,以及通过与模型中的低速区域相关联来识别过压区域。第一传感器可以位于第一区域的下表面附近。所述处理器可以进一步被配置为拾取与压缩波相对应的传播时间。所述处理器可以被配置为对所拾取的传播时间进行层析反演,以给出速度模型。所述处理器可以进一步被配置为识别地层中的低速区域。所述处理器可以进一步被配置为使用速度模型偏移记录信号。所述系统可以还包括被配置为将所述偏移的输出用于进一步的钻探操作的处理器。所述第一传感器可以包括3-C传感器和/或水下地震检波器。可以使用线缆将第一传感器传送到钻孔中。
本发明的另一实施例为一种与用于估计地层特性的系统一起使用的计算机可读介质,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差。所述系统包括被配置为在地面或地面附近的多个源位置处被激活的能量源。所述系统还包括被配置为生成响应于所述源的激活的信号的第一传感器,所生成的信号包括来自第一区域的下表面以及第一区域之下的地层中的界面的反射。所述介质包括使处理器进行以下处理的指令:在与所述反射相对应的记录信号中拾取传播时间,并根据所拾取的传播时间来估计第一区域的下表面的位置、所述界面的位置以及/或者所述第一区域和所述界面之间的速度。所述介质可以包括ROM、EPROM、EAROM、闪存和光盘。
附图说明
通过参考附图来最好地理解本发明,其中,相同的附图标记代表相同的部件,以及其中:
图1示出根据本发明的变井源距VSP的数据获取的几何原理;
图2是示出用于变井源距VSP仿真中的速度模型的图;
图3示出图2的模型的示例性射线路径几何图;
图4示出针对图2的模型获得的盐下反射的声照射(insonfication);
图5是针对图2的模型的位于单一深度处的示例性仿真VSP数据;
图6a示出示例性速度模型;
图6b示出用于层析反演的初始模型;
图6c示出对通过使用图6a的速度模型以及图6b所示的初始估计而仿真出的数据进行层析反演的结果;
图7a示出紧接在盐层之下的具有异常低速的示例性速度模型;
图7b示出用于对与图7a相对应的数据进行层析反演的初始模型;
图7c示出对通过使用图7a的速度模型以及图7b所示的初始估计而仿真出的数据进行层析反演的结果;以及
图8示出叠加在速度模型上的偏移后的变井源距VSP数据。
具体实施方式
对于本发明,变井源距VSP的获取几何原理在图1中示出。这里示出地面123上具有钻塔121。它可以是钻探钻塔(drilling rig),或者可以是将线缆传送到钻孔101中的主钻塔(mast rig)。钻孔101穿过层103、105、......。由111a、111b、111c、111d、...表示的地震传感器位于钻孔101中。各传感器可以包括水下地震检波器、单组件地震检波器或者多组件地震检波器。典型地,通过使用诸如位于地面处(或者位于地面处的水体中)的125a的单一地震源来获得单一偏移VSP的数据。出于本发明的目的,将水体的表面考虑为地面。射线127a示出示例性射线路径,该射线路径示出地震能量从源125a到检测器111d的传播,射线127a从位于边界106处的层105的底部被反射,并沿着129a所表示的射线路径到达接收器111d。
在典型的VSP中,将由于单一位置处的源(诸如125a)的运行而产生的数据记录在钻孔中的各接收器111a、111b、111c、111d...中。对反射数据的分析可以提供有关地下的震波速度以及层边界的构造的信息。在变井源距VSP中,针对位于多个源位置处的源(诸如125b、125c...)的运行而重复该处理。在多个检测器处获取来自多个源位置的数据提供对地下区域的冗余采样。这使得能够基于射线在各源与各接收器之间的传播时间来确定地下的速度。将使用传播时间来确定速度称为“层析反演”,并且用于对震波传播时间数据进行该层析反演的多个处理包可在商业上获得。
本发明的一创新点在于将变井源距VSP用于大地中的盐下反射成像的特定问题。其在图2中示出,其中,钻孔203穿过地层。横坐标为距离(该模型具有10000英尺或3.048km的横向范围),纵坐标为深度。该图中的等值线表示地下的P波速度,尤其是,地层201具有超过14000英尺/秒(4.3km/s)的速度。靠近钻孔203的底部配置有传感器阵列205。一或多个传感器位于盐中。
通过使用合适的仿真包,可以生成多个特征(diagnostic)显示。在本发明中,所使用的仿真包为Geo Tomo Inc.的VECONTM包。但这并不意图成为本发明的限制。可以获得的一种显示如图3所示,它示出与所选择的地下边界相对应的射线路径。在所示出的例子中,间隔210为盐层,该岩层的速度远高于紧挨在盐层之上的区域的速度。如本领域技术人员所已知的,这导致在上盐层(top-salt)边界处的大阻抗差。因此,源处所产生的许多地震能量在上盐层处被反射,从而导致向盐之下的层传播的能量减少。高速度的另一后果是在上盐层处的显著的射线弯曲(见211)。该射线弯曲导致在地下的不均匀的覆盖。
图4示出针对相同的模型边界所获得的覆盖(或声照射)。它们为图中的深色线221、231,它们示出:即使源的传播大约为8000英尺(见图3),地下也只有一小部分被声照射。术语“声照射”类似于用于光源的术语“照明”。由于高速度和大阻抗差而导致的问题并不限于盐,对于例如玄武岩的火山岩也可能产生这些问题。
在图5中示出针对图2的模型的位于单一接收器位置处的示例性VSP。纵坐标为时间,横坐标对应于地面处的源的编号。图形化地示出针对每个源位置在所选择的深度处的信号的时间序列。该时间序列是通过弹性波方程的有限差分解仿真得到的。被标识为301的“事件”是从不同源位置的直接到达。事件303是来自盐的(图2中的207)底部(the base of the salt)的反射。来自盐下的边界的反射由305、307和309表示。本发明的目的是绘制这些边界之间的速度的图。虽然图5中的例子针对合成数据,但是在实际的VSP获取领域中可以获得相似的记录。
为了实现该目的,拾取与所选择的传感器深度(图5的显示中所示)的数据中以及其它深度的类似数据中的各“事件”相对应的传播时间。这种所谓的“事件拾取”在本领域是熟知的,并且典型地涉及一震波轨迹与其它震波轨迹或参考信号的互相关。可以将其看作是针对各源位置和各传感器位置而定义传播时间的下列形式的矩阵:
其中,i为发射索引,j是接收器索引。典型地(但并非总是这样),发射位置的数量n大于接收器位置的数量m。然后对传播时间进行层析反演,以定义地下的速度模型。除拾取第一到达时间外,需要拾取(来自于盐下的界面的)盐下反射的时间,以用于反射层析。其可以在波场分离之后进行。如图5所示,可以容易地识别和拾取基盐和盐下反射(303、305、307和309)。
层析源自“截面图”的希腊语。将地下区域划分成单元,则数据被表示为沿着穿过单元的射线路径的线积分。透射层析涉及钻孔到钻孔、地面到钻孔或者地面到地面的观测。反射层析涉及地面到地面的观测(像在传统的反射或折射处理中一样)。在地震层析中,对各单元分配缓慢度(或速度),有时是将衰减系数分配到各单元,并且通过跟踪穿过模型的射线来计算传播时间(以及振幅)。将结果与所观测到的时间(以及振幅)进行比较;然后扰动该模型,并反复重复该处理以使误差最小化。每当假定的速度变化之后,必须重新计算射线路径。
VELMAPTM包包括反射层析反演算法。图6a示出用于生成传播时间矩阵的示例性速度模型。然后使用图6b的初始速度模型对所得到的传播时间矩阵进行反射层析反演。所得到的反射层析反演的输出如图6c所示。初始模型图6b与反演的最终结果之间的最大变化在于351所示的区域,其中,反演后的速度与图6a的实际速度相接近地匹配。由于图3和图4示出该区域中的反射体的最佳声照射,因此这是所希望的。由于射线路径不穿过模型的边缘这一简单原因,因此在图6b和图6c之间模型的边缘处的速度没有变化。
图7a是示例性模型,其中,盐201下面是低速区域203。这样的低速区域通常是由于过压地层而产生的,并且一般位于诸如盐的厚的不可渗透的地层之下。术语“过压”是指地层中的异常高的孔隙流体压力。钻探到过压区域是造成井喷的常见原因。因此,在钻探之前能够预测这样的过压是非常有价值的。
使用图7a的模型来生成传播时间数据。使用图7b的初始速度模型来对传播时间数据进行反射层析反演。其结果如图7c所示。清楚地看到紧挨在井底部之下的低速区域371。如本领域技术人员所知道的,地层中的震波速度是有效应力(覆层应力与地层孔隙压力之差)的函数。在过压地层中,地层孔隙压力很高,从而使得有效应力变小。因此,过压区域中的震波速度很低。识别钻头前面的这样的过压区域对避免井喷非常有用,因为可以采取诸如增大钻探时所使用的泥浆重量的补救措施。
基于所估计的盐下速度分布,对合成的地震数据应用叠前基尔霍夫偏移,以对盐下界面成像。在偏移中,采用精确、稳定的快速行进法(fast marching method)来计算传播时间表(Lou,2006)。图8是对来自变井源距VSP测量的所有接收器集合的合成地震波图的偏移结果。通过将偏移结果叠加到真实的速度模型上,可以看到,该偏移使得三个盐下界面351、353、357正确地、可靠地成像,这证明了本发明的正确性和有效性。特别地,它增强了人们对速度模型以及所解释的过压区域的信心。
在本发明的一实施例中,将钻孔钻到恰好在预期的基盐深度上的深度。将带有多个传感器的线缆传送到钻孔中,并进行变井源距VSP。根据变井源距VSP,如上所述确定速度模型。识别钻孔的底部前面的反射界面,并确定盐下速度。然后可以基于新确定的盐下反射界面通过针对钻孔方向上的预期的盐下孔隙压力状况而适当地选择泥浆重量,来继续钻孔。可选地,在进行钻孔的过程中可以进行变井源距VSP。
可以在通用数字计算机上来实现上述反演和偏移方法。如本领域技术人员所熟知的,计算机指令存在诸如ROM、EPROM、EAROM、闪存和光盘的机器可读存储装置上。这些存储装置可以是计算机的一部分,或者可以通过合适的通信通道链接到计算机上并且甚至可以位于远程位置。类似地,上述类型的地震数据可以存储在计算机上,或者可以通过合适的通信通道链接到计算机上。通信通道可以包括互联网,使得用户可以从一远程位置访问数据并从另一远程位置获得指令,来处理数据。机器可读存储装置上的指令使计算机可以访问多元数据,并根据上述方法对数据进行处理。
尽管前述公开针对本发明的优选实施例,但是对于本领域技术人员来说各种修改是显而易见的。本说明书意在使前述公开包含所附权利要求书范围和精神之内的所有这样的改变。
权利要求书(按照条约第19条的修改)
1.一种用于估计地层的特性的方法,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述方法包括:
(a)在第一区域中布置第一传感器;
(b)激活在地面或地面附近的多个源位置处的能量源;
(c)记录来自第一传感器、响应于所述源的激活的信号,记录信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;以及
(d)根据记录信号来估计至少下列之一:(I)第一区域的下表面的位置,(II)所述界面的位置;以及(III)所述第一区域和所述界面之间的速度;
所述方法的特征在于:
(e)在与所述反射相对应的记录信号中拾取传播时间;以及
(f)使用所拾取的传播时间进行所述估计。
2.根据权利要求1所述的方法,还包括在第一区域的下表面附近布置第一传感器。
3.根据权利要求1所述的方法,其中第一区域包括盐。
4.根据权利要求1所述的方法,其中布置第一传感器进一步包括布置多个传感器。
5.根据权利要求1所述的方法,其中所拾取的传播时间对应于压缩波传播时间。
6.根据权利要求1所述的方法,其中所述估计还包括对所拾取的传播时间进行层析反演以给出速度模型。
7.根据权利要求1所述的方法,还包括在步骤(d)之前进行记录信号的波场分离。
8.根据权利要求6所述的方法,还包括识别地层中的低速区域。
9.根据权利要求6所述的方法,还包括使用所述速度模型来偏移记录信号。
10.根据权利要求9所述的方法,还包括基于所述速度模型进行钻探操作。
11.根据权利要求1所述的方法,所述地层还包括过压区域,所述方法还包括:
(i)记录来自第一传感器、响应于来自过压区域之下的地层中的界面的反射的信号;
(ii)在与所述反射相对应的记录信号中拾取传播时间;
(iii)对所拾取的传播时间进行层析反演,以给出速度模型;以及
(iv)通过与模型中的低速区域相关联来识别过压区域。
12.根据权利要求11所述的方法,还包括在(ii)之前进行记录信号的波场分离。
13.根据权利要求12所述的方法,还包括基于所识别的过压区域来进行钻探操作。
14.一种用于估计地层的特性的系统,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一传感器,其被配置为生成响应于所述源的激活的信号,所生成的信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;
(c)处理器,其被配置为根据记录信号来估计至少下列之一:(I)第一区域的下表面的位置,(II)所述界面的位置,以及(III)所述第一区域和所述界面之间的速度;
所述系统的特征在于,所述处理器被配置为:
(d)从与所述反射相对应的记录信号中拾取传播时间;以及
(e)使用所拾取的传播时间进行所述估计。
15.根据权利要求14所述的系统,其中第一传感器位于第一区域的下表面附近。
16.根据权利要求14所述的系统,其中第一区域包括盐。
17.根据权利要求14所述的系统,其中第一传感器进一步包括多个传感器。
18.根据权利要求14所述的系统,其中所述处理器进一步被配置为拾取与压缩波相对应的传播时间。
19.根据权利要求15所述的系统,其中所述处理器被配置为通过进一步对所拾取的传播时间进行层析反演以给出速度模型,来进行所述估计。
20.根据权利要求19所述的系统,所述处理器还被配置为识别地层中的低速区域。
21.根据权利要求20所述的系统,其中所述处理器还被配置为使用所述速度模型来偏移记录信号。
22.根据权利要求21所述的系统,还包括被配置为使用所述偏移的输出进行进一步的钻探操作的处理器。
23.根据权利要求14所述的系统,其中第一传感器进一步包括下列至少之一:(i)3-C传感器,以及(ii)水下地震检波器。
24.根据权利要求14所述的系统,其中使用线缆将第一传感器传送到钻孔中。
25.根据权利要求14所述的系统,其中所述地层还包括过压区域,所述系统进一步包括:
(i)所述第一传感器还被配置为记录响应于来自过压区域之下的地层中的界面的反射的信号;以及
(ii)所述处理器还被配置为:
(A)在与所述反射相对应的记录信号中拾取传播时间;
(B)对所拾取的传播时间进行层析反演,以给出速度模型;以及
(C)通过与模型中的低速区域相关联来识别过压区域。
26.根据权利要求25所述的系统,其中所述处理器还被配置为在(d)之前进行记录信号的波场分离。
27.根据权利要求25所述的系统,还包括被配置为基于所识别的过压区域来进行钻探操作的处理器。
28.一种与用于估计地层的特性的系统一起使用的计算机可读介质,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一传感器,其被配置为生成响应于所述源的激活的信号,所生成的信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;
所述介质包括使处理器进行以下处理的指令:
(c)根据记录信号来估计至少下列之一:(I)第一区域的下表面的位置,(II)所述界面的位置,以及(III)所述第一区域和所述界面之间的速度;
所述介质的特征在于,所述指令使所述处理器进行以下处理:
(d)从与所述反射相对应的记录信号中拾取传播时间;以及
(e)使用所拾取的传播时间进行所述估计。
29.根据权利要求28所述的介质,进一步包括至少下列之一:(i)ROM,(ii)EPROM,(iii)EAROM,(iv)闪存,以及(v)光盘。
Claims (39)
1.一种用于估计地层的特性的方法,该地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述方法包括:
(a)在第一区域中布置第一传感器;
(b)激活在地面或地面附近的多个源位置处的能量源;
(c)记录来自第一传感器、响应于所述源的激活的信号,记录信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;
(d)在与所述反射相对应的记录信号中拾取传播时间;以及
(e)根据所拾取的传播时间来估计至少下列之一:(i)第一区域的下表面的位置,(ii)所述界面的位置;以及(iii)所述第一区域和所述界面之间的速度。
2.根据权利要求1所述的方法,还包括在第一区域的下表面附近布置第一传感器。
3.根据权利要求1所述的方法,其中第一区域包括盐。
4.根据权利要求1所述的方法,其中布置第一传感器进一步包括布置多个传感器。
5.根据权利要求1所述的方法,其中所拾取的传播时间对应于压缩波传播时间。
6.根据权利要求1所述的方法,其中所述估计进一步包括对所拾取的传播时间进行层析反演以给出速度模型。
7.根据权利要求1所述的方法,还包括在步骤(d)之前进行记录信号的波场分离。
8.根据权利要求6所述的方法,还包括识别地层中的低速区域。
9.根据权利要求6所述的方法,还包括使用所述速度模型来偏移记录信号。
10.根据权利要求9所述的方法,还包括基于所述速度模型进行钻探操作。
11.一种用于估计地层的特性的系统,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一传感器,其被配置为生成响应于所述源的激活的信号,所生成的信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;
(c)处理器,其被配置为:
(A)从与所述反射相对应的记录信号中拾取传播时间;
以及
(B)根据所拾取的传播时间来估计至少下列之一:(I)第一区域的下表面的位置,(II)所述界面的位置,以及(III)所述第一区域和所述界面之间的速度。
12.根据权利要求11所述的系统,其中所述第一传感器位于第一区域的下表面附近。
13.根据权利要求11所述的系统,其中第一区域包括盐。
14.根据权利要求11所述的系统,其中第一传感器进一步包括多个传感器。
15.根据权利要求11所述的系统,其中所述处理器进一步被配置为拾取与压缩波相对应的传播时间。
16.根据权利要求11所述的系统,其中所述处理器被配置为通过进一步对所拾取的传播时间进行层析反演以给出速度模型,来进行所述估计。
17.根据权利要求16所述的系统,所述处理器还被配置为识别地层中的低速区域。
18.根据权利要求17所述的系统,其中所述处理器还被配置为使用所述速度模型来偏移记录信号。
19.根据权利要求18所述的系统,还包括被配置为使用所述偏移的输出进行进一步的钻探操作的处理器。
20.根据权利要求11所述的系统,其中所述第一传感器进一步包括下列的至少之一:(i)3-C传感器,以及(ii)水下地震检波器。
21.根据权利要求11所述的系统,其中使用线缆将第一传感器传送到钻孔中。
22.一种与用于估计地层的特性的系统一起使用的计算机可读介质,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一传感器,其被配置为生成响应于所述源的激活的信号,其中,所生成的信号包括来自下列的反射:(A)第一区域的下表面,以及(B)第一区域之下的地层中的界面;
所述介质包括使处理器进行以下处理的指令:
(c)从与所述反射相对应的记录信号中拾取传播时间;以及
(d)根据所拾取的传播时间来估计至少下列之一:(I)第一区域的下表面的位置,(II)所述界面的位置,以及(III)所述第一区域和所述界面之间的速度。
23.根据权利要求22所述的介质,进一步包括至少下列之一:(i)ROM,(ii)EPROM,(iii)EAROM,(iv)闪存,以及(v)光盘。
24.一种用于识别地层中的过压区域的方法,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述方法包括:
(a)在第一区域中布置第一传感器;
(b)激活在地面或地面附近的多个源位置处的能量源;
(c)记录来自第一传感器、响应于所述源的激活的信号,记录信号包括来自下列的反射:(A)第一区域的下表面,以及(B)过压区域之下的地层中的界面;
(d)在与所述反射相对应的记录信号中拾取传播时间;
(e)对所拾取的传播时间进行层析反演,以给出速度模型;以及
(f)通过与模型中的低速区域相关联来识别过压区域。
25.根据权利要求24所述的方法,还包括在第一区域的下表面附近布置第一传感器。
26.根据权利要求24所述的方法,其中布置第一传感器进一步包括布置多个传感器。
27.根据权利要求24所述的方法,其中所拾取的传播时间对应于压缩波传播时间。
28.根据权利要求24所述的方法,还包括在(d)之前进行记录信号的波场分离。
29.根据权利要求24所述的方法,还包括基于所识别的过压区域来进行钻探操作。
30.一种用于识别地层中的过压区域的系统,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一区域中的第一传感器,其被配置为记录响应于所述源的激活的信号,记录信号包括来自下列的反射:(A)第一区域的下表面,以及(B)过压区域之下的地层中的界面;以及
(c)处理器,其被配置为:
(A)在与所述反射相对应的记录信号中拾取传播时间;
(B)对所拾取的传播时间进行层析反演,以给出速度模型;以及
(C)通过与模型中的低速区域相关联来识别过压区域。
31.根据权利要求30所述的系统,其中第一传感器位于第一区域的下表面附近。
32.根据权利要求30所述的系统,其中第一传感器进一步包括多个传感器。
33.根据权利要求30所述的系统,其中所拾取的传播时间对应于压缩波传播时间。
34.根据权利要求30所述的系统,其中所述处理器还被配置为在(d)之前进行记录信号的波场分离。
35.根据权利要求30所述的系统,还包括被配置为基于所识别的过压区域来进行钻探操作的处理器。
36.根据权利要求30所述的系统,其中第一传感器进一步包括至少下列之一:(i)3-C传感器,以及(ii)水下地震检波器。
37.根据权利要求30所述的系统,其中使用线缆将第一传感器传送到钻孔中。
38.一种与用于识别地层中的过压区域的系统一起使用的计算机可读介质,所述地层包括第一区域,所述第一区域与位于其上的第二区域之间具有大的阻抗差,所述系统包括:
(a)被配置为在地面或地面附近的多个源位置处被激活的能量源;
(b)第一区域中的第一传感器,其被配置为记录响应于所述源的激活的信号,记录信号包括来自下列的反射:(A)第一区域的下表面,以及(B)过压区域之下的地层中的界面;
所述介质包括使处理器进行以下处理的指令:
(c)在与所述反射相对应的记录信号中拾取传播时间;
(d)对所拾取的传播时间进行层析反演,以给出速度模型;以及
(e)通过与模型中的低速区域相关联来识别过压区域。
39.根据权利要求38所述的介质,进一步包括至少下列之一:(i)ROM,(ii)EPROM,(iii)EAROM,(iv)闪存,以及(v)光盘。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US79727706P | 2006-05-03 | 2006-05-03 | |
US60/797,277 | 2006-05-03 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101568856A true CN101568856A (zh) | 2009-10-28 |
Family
ID=38668322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007800234002A Pending CN101568856A (zh) | 2006-05-03 | 2007-05-02 | 通过变井源距垂直地震剖面测量进行的盐下反射层析和成像 |
Country Status (6)
Country | Link |
---|---|
US (2) | US7911878B2 (zh) |
CN (1) | CN101568856A (zh) |
CA (1) | CA2651209C (zh) |
GB (1) | GB2451611A (zh) |
NO (1) | NO20084704L (zh) |
WO (1) | WO2007130550A2 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104781699A (zh) * | 2012-09-17 | 2015-07-15 | 贝克休斯公司 | 床内源垂直地震剖面获取 |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7911878B2 (en) * | 2006-05-03 | 2011-03-22 | Baker Hughes Incorporated | Sub-salt reflection tomography and imaging by walkaway VSP survey |
US7586811B2 (en) * | 2006-06-02 | 2009-09-08 | Westerngeco L.L.C. | Subsalt velocity model building |
US8885440B2 (en) * | 2007-01-05 | 2014-11-11 | Madhumita Sengupta | Constructing velocity models near salt bodies |
US20100118654A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling migration method |
US8750074B2 (en) * | 2008-11-25 | 2014-06-10 | Baker Hughes Incorporated | Anisotropic parameter determination |
US8547794B2 (en) * | 2009-04-16 | 2013-10-01 | Baker Hughes Incorporated | Extending the coverage of VSP/CDP imaging by using first-order downgoing multiples |
US8576663B2 (en) * | 2010-04-30 | 2013-11-05 | Schlumberger Technology Corporation | Multicomponent seismic inversion of VSP data |
CN102053267B (zh) * | 2010-10-22 | 2012-11-07 | 中国石油化工股份有限公司 | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 |
EP2715403A4 (en) * | 2011-05-27 | 2015-01-28 | Conocophillips Co | TARGETED DATA SELECTION VIA A BIDIRECTIONAL WAVE EQUATION FOR IMPROVED PROSPECTS IMAGING COMPLEX GEOLOGICAL STRUCTURES |
US9116255B2 (en) * | 2011-05-27 | 2015-08-25 | Conocophillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
FR2984525B1 (fr) * | 2011-12-15 | 2014-01-17 | Cggveritas Services Sa | Separation de champs d'onde pour des enregistrements sismiques repartis sur des surfaces d'enregistrement non planes |
CN102879805B (zh) * | 2012-10-24 | 2015-06-24 | 北京市市政工程研究院 | 一种基于钻孔与地面相结合的地震波空间探测方法 |
US20140200816A1 (en) * | 2013-01-14 | 2014-07-17 | Westerngeco L.L.C. | Seismic data processing |
AU2014238094A1 (en) | 2013-03-22 | 2015-07-02 | Halliburton Energy Services, Inc. | Migration velocity analysis method for vertical seismic profile data |
CN104570106A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种近地表层析速度分析方法 |
GB2528130A (en) * | 2014-07-11 | 2016-01-13 | Total E&P Uk Ltd | Method of constraining an inversion in the characterisation of the evolution of a subsurface volume |
CN105319598A (zh) * | 2014-07-31 | 2016-02-10 | 中国石油化工股份有限公司 | 一种不需要拾取的叠前地震反射层析方法 |
US9702999B2 (en) | 2014-10-17 | 2017-07-11 | Chevron U.S.A. Inc. | System and method for velocity analysis in the presence of critical reflections |
CN105301639B (zh) * | 2015-10-21 | 2017-09-01 | 中国石油天然气集团公司 | 基于vsp旅行时双加权层析反演速度场的方法及其装置 |
US10444063B2 (en) | 2016-09-23 | 2019-10-15 | Baker Hughes, A Ge Company, Llc | Downhole fiber optic hydrophone |
CN106772609A (zh) * | 2017-03-03 | 2017-05-31 | 中国石油天然气集团公司 | 低速区底界面确定方法和装置 |
US11867855B2 (en) | 2019-12-10 | 2024-01-09 | Baker Hughes Oilfield Operations Llc | Downhole fiber optic hydrophone |
US11573346B2 (en) * | 2021-04-15 | 2023-02-07 | Saudi Arabian Oil Company | Determining a seismic quality factor for subsurface formations for marine vertical seismic profiles |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3371310A (en) * | 1966-03-14 | 1968-02-27 | Pan American Petroleum Corp | Discriminating between primary and multiple seismic reflections |
US4794573A (en) * | 1988-02-11 | 1988-12-27 | Conoco Inc. | Process for separating upgoing and downgoing events on vertical seismic profiles |
US4835745A (en) * | 1988-05-27 | 1989-05-30 | Exxon Production Research Company | Method of radial seismic data collection for improved salt structure imaging |
US4926393A (en) | 1989-01-23 | 1990-05-15 | Conoco Inc. | Multifold vertical seismic profile acquisition method and technique for imaging the flank of a salt dome |
US4918670A (en) * | 1989-02-15 | 1990-04-17 | Conoco Inc. | Method for interval velocity analysis and determination of reflector position from before stack seismic data |
US5170377A (en) * | 1990-03-22 | 1992-12-08 | Western Atlas International, Inc. | 3-D mapping of salt domes |
US5307268A (en) * | 1991-06-07 | 1994-04-26 | Conoco Inc. | Method for using P-S converted waves as a subsalt imaging tool |
US5596548A (en) | 1994-05-12 | 1997-01-21 | Exxon Production Research Company | Seismic imaging using wave equation extrapolation |
US5995446A (en) * | 1998-04-21 | 1999-11-30 | Schlumberger Technology Corporation | Method of conducting drilling operations using vertical seismic profiles |
GB2343951B (en) * | 1998-11-20 | 2003-05-14 | Thomson Marconi Sonar Ltd | Drilling apparatus |
GB0018480D0 (en) * | 2000-07-27 | 2000-09-13 | Geco Prakla Uk Ltd | A method of processing surface seismic data |
FR2821677B1 (fr) * | 2001-03-05 | 2004-04-30 | Geophysique Cie Gle | Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees |
US6751558B2 (en) * | 2001-03-13 | 2004-06-15 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
FR2831961B1 (fr) * | 2001-11-07 | 2004-07-23 | Inst Francais Du Petrole | Methode de traitement de donnees sismiques de puits en amplitude preservee absolue |
US6643590B2 (en) * | 2002-01-04 | 2003-11-04 | Westerngeco, L.L.C. | Method for computing finite-frequency seismic migration traveltimes from monochromatic wavefields |
GB2384053B (en) * | 2002-01-11 | 2005-09-14 | Westerngeco Ltd | A method of and apparatus for processing seismic data |
US6894949B2 (en) * | 2002-10-04 | 2005-05-17 | Baker Hughes Incorporated | Walkaway tomographic monitoring |
CA2577128C (en) * | 2003-08-19 | 2014-10-21 | Tetra Seis Inc. | Method for interpreting seismic data using duplex waves |
US7508733B2 (en) * | 2003-11-14 | 2009-03-24 | Schlumberger Technology Corporation | High-frequency processing of seismic vibrator data |
US7911878B2 (en) * | 2006-05-03 | 2011-03-22 | Baker Hughes Incorporated | Sub-salt reflection tomography and imaging by walkaway VSP survey |
US7586811B2 (en) * | 2006-06-02 | 2009-09-08 | Westerngeco L.L.C. | Subsalt velocity model building |
US7426438B1 (en) * | 2007-03-16 | 2008-09-16 | Westerngeco L.L.C. | Technique to provide seismic data-based products |
-
2007
- 2007-04-05 US US11/697,125 patent/US7911878B2/en not_active Expired - Fee Related
- 2007-04-05 US US11/697,049 patent/US7751279B2/en not_active Expired - Fee Related
- 2007-05-02 CA CA2651209A patent/CA2651209C/en not_active Expired - Fee Related
- 2007-05-02 WO PCT/US2007/010793 patent/WO2007130550A2/en active Application Filing
- 2007-05-02 CN CNA2007800234002A patent/CN101568856A/zh active Pending
-
2008
- 2008-11-10 NO NO20084704A patent/NO20084704L/no not_active Application Discontinuation
- 2008-12-02 GB GB0822039A patent/GB2451611A/en not_active Withdrawn
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104781699A (zh) * | 2012-09-17 | 2015-07-15 | 贝克休斯公司 | 床内源垂直地震剖面获取 |
Also Published As
Publication number | Publication date |
---|---|
GB2451611A (en) | 2009-02-04 |
US20070258324A1 (en) | 2007-11-08 |
NO20084704L (no) | 2009-01-29 |
WO2007130550A2 (en) | 2007-11-15 |
GB2451611A8 (en) | 2009-02-18 |
WO2007130550A3 (en) | 2008-04-24 |
US20070259432A1 (en) | 2007-11-08 |
CA2651209A1 (en) | 2007-11-15 |
US7911878B2 (en) | 2011-03-22 |
GB0822039D0 (en) | 2009-01-07 |
WO2007130550B1 (en) | 2008-10-02 |
CA2651209C (en) | 2013-07-30 |
US7751279B2 (en) | 2010-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101568856A (zh) | 通过变井源距垂直地震剖面测量进行的盐下反射层析和成像 | |
US6928367B2 (en) | Reservoir fracture characterization | |
US8547794B2 (en) | Extending the coverage of VSP/CDP imaging by using first-order downgoing multiples | |
Tang et al. | Processing array acoustic-logging data to image near-borehole geologic structures | |
Bakulin et al. | Smart DAS upholes for simultaneous land near-surface characterization and subsurface imaging | |
US7492664B2 (en) | Method for processing acoustic reflections in array data to image near-borehole geological structure | |
WO2008124759A1 (en) | Location of a microseismic event using the combination of a vsp array and a surface array | |
US6591193B2 (en) | Method and apparatus for acquiring offset checkshot survey data using tube-wave conversion | |
BR112014021827A2 (pt) | extração de dados de cisalhamento sv corte a partir de dados marinhos de onda p | |
US20060136133A1 (en) | Signal thresholding apparatus, systems, and methods | |
CN104781699A (zh) | 床内源垂直地震剖面获取 | |
US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
Pasasa et al. | Prestack Kirchhoff depth migration of shallow seismic data | |
Pugin et al. | First-arrival alignment static corrections applied to shallow seismic reflection data | |
US8857254B2 (en) | Methods and systems for acquiring acceleration waveforms in a borehole | |
US20130094325A1 (en) | Look-ahead seismic while drilling | |
Bergman et al. | High-resolution reflection seismic imaging of the upper crust at Laxemar, southeastern Sweden | |
Urosević et al. | Seismic anisotropy in the upper 500 m of the Southern Sydney Basin | |
US11474272B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
RU2305856C1 (ru) | Способ наземно-скважинной сейсморазведки | |
US20240176038A1 (en) | Method and apparatus for source wavelet estimation | |
Hakim et al. | Subsurface velocity measurement | |
Zhou et al. | Normal and Shear Compliance Estimation for Inclined Fractures Using Full-Waveform Sonic Log Data | |
McLeod | Seismology Near and Far: Local Earthquake Monitoring and Teleseismic Reflection Imaging | |
Adhiansyah et al. | Fit For Purpose Integrated Structural Identification Using Azimuthal Monopole Sonic Waveforms and Vertical Seismic Profiling: The First Case Study From The Vertical Unconventional Well in North Sumatera Basin-Indonesia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20091028 |