CN111971586A - 地震速度衍生的烃指示 - Google Patents

地震速度衍生的烃指示 Download PDF

Info

Publication number
CN111971586A
CN111971586A CN201980023710.7A CN201980023710A CN111971586A CN 111971586 A CN111971586 A CN 111971586A CN 201980023710 A CN201980023710 A CN 201980023710A CN 111971586 A CN111971586 A CN 111971586A
Authority
CN
China
Prior art keywords
seismic
data
velocity
model
region
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
Application number
CN201980023710.7A
Other languages
English (en)
Inventor
傅康
戴安娜·妮
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.)
BP Corp North America Inc
Original Assignee
BP Corp North America Inc
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 BP Corp North America Inc filed Critical BP Corp North America Inc
Publication of CN111971586A publication Critical patent/CN111971586A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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/282Application of seismic models, synthetic seismograms
    • 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/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • 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/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • 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/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/51Migration
    • G01V2210/512Pre-stack
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (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

地震速度衍生的烃指示
背景技术
本公开总体上涉及分析地震数据,以及更具体地涉及生成一种对地震勘测区域内存在的烃进行成像的属性模型。
本节旨在向读者介绍技术的多种方面,所述技术可能与在下文中描述和/或要求保护的本公开的多种方面有关。相信这个讨论对于给读者提供背景信息以帮助更好地理解本公开的多种方面是有益的。因此,应当理解,要从这个角度来阅读这些陈述,而不是作为现有技术的承认。
地震勘测包括通过向地下发送声能并记录从地下区域内的地质层返回的反射声能来生成地球的地下区域的图像或地图。在地震勘测期间,能量源被放置在地球的表面区域之上或上方的多种位置处,所述区域可能包括烃沉积层(例如,地下区域内包含烃的区域)。每一次激活所述源,所述源生成地震(例如,声波)信号,所述信号向下传播通过地球,被反射,并在它返回时用一个或多个位于地球的地下区域之上或上方的接收器记录。通过接收器记录的地震数据可以随后被用于创建相应的地下区域的图像或轮廓。
随着时间推移,当烃从地球的地下区域被提取出来时,地下区域内烃储池的位置、饱和度以及其他特征和(例如,覆盖层)可能变化。因此,确定地下区域的图像或地图如何随着时间变化可能是有用的,因为这样可以修改与烃提取相关的操作以更有效地从地球的地下区域提取烃。
发明内容
下面阐述本文公开的某些实施方式的概述。应当理解,呈现这些方面仅仅是为了向读者提供这些实施方式的简要概述,并且这些方面不旨在限制本公开的范围。事实上,本公开可以包括可能在下文中未阐述的多种方面。
利用源和接收器的地震获取可能在生成例如地震图像中是有用的。地震图像可以用于例如确定烃沉积层(例如,地下之内含烃的区域)和/或地下钻探危害物。地震图像通常根据由所述源产生的地震波形产生,所述波形从地下内的区域反射出来并被接收器所接收。地震图像可以利用通过例如速度模型构建器应用的高分辨率地震速度模型例如全波形反演(FWI)模型、层析成像模型等来构建。地震速度模型包括指示随着地震波形传播通过地下区域时地震波形速度的变化的数据。地震波形速度的变化指示存在的流体类型(盐水、油或气体)。然而,速度的灵敏度随着地下区域内深度的增加和/或根据各自的地层属性(例如,岩石孔隙度等)而减小。因此,用地震速度模型可能难以确定位于地下增加的深度处和/或位于特定部分的地层之下的烃沉积层的位置。因此,本发明的实施方式涉及用于产生一种基于地震速度模型的替代属性模型的技术,所述替代属性模型包括代表区域的物理属性的数据。
应用校准或其他调整来确定地震速度,以补偿在增加的深度处和/或对地层的特定属性(例如,岩石密度)作出响应而减小的速度灵敏度。速度校准可以表示为类似于流体饱和度的属性,因此被称作基于速度的饱和度(SV)。在一个实施方式中,用地震速度和从测井日志获得的速度数据的组合来确定SV。在另一个实施方式中,SV衍生自一种假设,即地震速度是各自权重为1-SV和SV的岩石在盐水中的沉积速度和充满烃的岩石的沉积速度的几何平均值。在另一个实施方式中,SV是沙或页岩的沉积速度的函数。用这种方法,可以从地震速度获得地层中存在烃的指示。在一个优选的实施方式中,所述速度通过全波形反演(FWI)产生。
附图说明
在阅读以下详细说明并参考附图后,可以更好地理解本公开的多种方面,其中:
图1说明了依据本文所呈现的实施方式的多种处理的流程图,所述处理可以根据通过地震勘测系统获取的地震数据的分析来进行;
图2说明了依据本文所呈现的实施方式的在海洋环境中的海洋勘测系统;
图3说明了依据本文所呈现的实施方式的在非海洋环境中的陆地勘测系统。
图4说明了依据本文所呈现的实施方式的计算系统,所述计算系统可以根据通过图2的海洋勘测系统和/或图3的第二海洋勘测系统获取的数据进行本文所描述的操作。
图5描述了依据本文所呈现的实施方式通过图4的计算系统生成的地震图像的实施例。
图6描述了依据本文所呈现的实施方式通过图4的计算系统生成的地震图像的另一个实施例;且
图7描述了依据本文所呈现的实施方式通过图4的计算系统生成地震图像的方法的流程图。
具体实施方式
下文描述了一种或多种具体实施方式。为了提供这些实施方式的简要说明,没有在说明书中描述实际实施的所有特征。应当理解,在任何这样的实际实施的开发中,如同在任何工程改造或设计项目中,必须作出大量具体实施的决定以实现开发者的具体目标,例如遵循系统相关和业务相关的限制,这可能在实施之间彼此不同。此外,应当理解,这样的开发工作可能是复杂和耗时的,但对那些受益于本公开的普通技术人员而言,这将是常规的设计、制造和生产工作。
地震数据可以提供关于例如地球的地下区域内的烃沉积层的位置和/或变化的说明的有价值信息。此外,可以利用不同的技术来完成地震数据的处理,所述技术可以根据例如待分析和/或成像的地层的地层属性和/或不同深度来选择。因此,在一个实施方式中,可以通过选择特定的处理类型来生成地层的特定部分的更高分辨率图像(例如,更高质量的图像)(即,根据地层属性和/或地层的深度修改所选择的地震数据处理技术以实现地层的特定部分的更高分辨率的图像)。
作为介绍,可以用多种地震勘测系统和技术获取地震数据,参考图2和图3讨论了其中的两种。不管所利用的地震数据收集技术如何,在获取地震数据后,计算系统可以分析所获取的地震数据,并可以用所述地震数据分析的结果(例如,地震记录图,地质地层地图等)来进行烃勘探和生产工业内的多种操作。例如,图1说明了方法10的流程图,所述方法详细描述了根据获取的地震数据的分析可以采用的多种处理。尽管以特定的顺序描述了方法10,但应当注意到可以用任何适合的顺序进行方法10。
现在参考图1,在框12处,可以根据所分析的地震数据来确定与各自地震勘测相关联的地球的地下区域内烃沉积层的位置和性质。在一个实施方式中,可以分析通过一种或多种地震获取技术获取的地震数据,以生成说明地下区域内多种地质地层的地图或轮廓。
根据所鉴定的烃沉积层的位置和性质,在框14处,可以勘探地下区域的特定位置或部分。也就是说,烃勘探机构可以利用烃沉积层的位置来确定在所述地下区域的表面处的位置,以钻探进入地球。因此,烃勘探机构可以利用烃沉积层的位置和性质以及相关联的覆盖层来测定钻探进入地球的路径、如何钻探进入地球等。
在勘探设备被放置在地下区域内之后,在框16处,烃沉积层中储存的烃可以通过天然自喷井、人工举升井等产生。在框18处,所产生的烃可以通过运输工具、管道等被运输至炼油厂。在框20处,可以根据多种精炼程序处理所产生的烃,以利用所述烃来开发不同产品。
应当注意,根据方法10所讨论的处理可以包括其他适合的处理,所述其他适合的处理可以基于烃沉积层的位置和性质,如通过一种或多种地震勘测获取的地震数据所指示的。因此,应当理解,上述处理不旨在描述在确定地下区域内烃沉积层的位置和性质后可进行的处理的详尽列表。
考虑到上述情况,图2说明了海洋勘测系统22(例如,与图1的框12一起使用),所述海洋勘测系统可以用于在海洋环境中获取关于地球的地下区域的地震数据(例如波形)。通常,使用海洋勘测系统22的海洋地震勘测可以在海洋24或其他水体中进行,所述其他水体在位于海底28下方的地球的地下区域26上方。
海洋勘测系统22可以包括船舶30、地震源32、拖缆34、接收器36和/或其他设备,所述其他设备可以帮助获取代表在地球的地下区域26内的地质地层的地震图像。船舶30可以牵引地震源32(例如气枪阵列),所述地震源可以产生导向海底28的能量,例如声波(例如地震波形)。船舶30也可以牵引具有接收器36(例如水听器)的拖缆34,所述接收器可以获取地震波形,所述地震波形代表被地震源32输出后随之被地下区域26内的多种地质地层反射出来的能量。尽管接收器36被描述成被船舶30牵引,但在某些实施方式中,接收器36也可以被部署在海底28的表面上。此外,尽管海洋勘测系统22的说明用一种地震源32(在图2中表示为气枪阵列)和一种接收器36(在图2中表示为多个水听器)来描述,但应当注意,海洋勘测系统22可以包括多种地震源32和多种地震接收器36。同样地,尽管海洋勘测系统22的以上说明用一种地震拖缆34来描述,但应当注意,海洋勘测系统22可以包括多种地震拖缆34。此外,其他船舶30可以包括其他源32、拖缆34等以进行勘测系统的操作。
图3说明了陆地勘测系统38(例如,与图1的框12一起使用),所述陆地勘测系统可以用于在非海洋环境中获得关于地球的地下区域26的信息。陆地勘测系统38可以包括基于陆地的地震源40和基于陆地的接收器44。在某些实施方式中,陆地勘测系统38可以包括一种或多种地震源40和一种或多种接收器44和46。事实上,出于讨论的目的,图3包括基于陆地的地震源40以及两个地震接收器44和46。基于陆地的地震源40(例如地震振动器)可以被放置在地球的表面42上,所述表面在感兴趣的地下区域26的上方。基于陆地的地震源40可以产生导向地球的地下区域26的能量(例如声波、地震波形)。在接触到地下区域26内的多种地质地层(例如盐丘、断层、褶皱)后,基于陆地的地震源40输出的能量可以被地质地层反射出来并被一种或多种基于陆地的接收器(例如44和46)获取或记录。
在某些实施方式中,基于陆地的接收器44和46可以被分散于地球的表面42以形成网格图案。因此,对通过地震源40被导向地下区域26的能量作出响应,每个基于陆地的接收器44或46可以接收所反射的地震波形。在某些情况下,由地震源40产生的一个地震波形可以被不同的地质地层反射出来并被不同的接收器所接收。例如,如图3所示,地震源40可以输出可能导向地下区域26的能量作为地震波形48。第一个接收器44可以接收一个地质地层对地震波形48的反射,且第二个接收器46可以接收不同的地质地层对地震波形48的反射。因此,第一个接收器44可以接收反射的地震波形50,且第二个接收器46可以接收反射的地震波形52。
不管地震数据是如何获取的,计算系统(例如,与图1的框12一起使用)可以分析由基于海洋的接收器36或基于陆地的接收器44和46获取的地震波形,以确定关于在地下区域26内地质结构、烃沉积层的位置和性质等的信息。图4说明了这种计算系统60的实施例,所述计算系统可以进行多种数据分析操作,以分析由接收器36、44或46获取的地震数据,从而确定在地下区域26内的地质地层的结构。
现在参考图4,计算系统60可以包括通信部件62、处理器64、记忆器66、存储器68、输入/输出(I/O)端口70、显示器72等。通信部件62可以是无线或有线的通信部件,所述通信部件可以促进接收器36、44、46,一个或多个数据库74、其他计算装置和其他能够通信的装置之间的通信。在一个实施方式中,计算系统60可以接收接收器数据76(例如,地震数据、地震记录图),所述接收器数据可能之前已经由地震接收器通过网络部件、数据库74等获取。计算系统60的处理器64可以分析或处理接收器数据76,以确定关于地球的地下区域26内的地质地层的多种特征。
处理器64可以是任何类型的能够执行计算机可执行代码的计算机处理器或微处理器。处理器64也可以包括可以进行下述操作的多种处理器。记忆器66和存储器68可以是任何适合的制造品,所述制造品可以用作存储处理器可执行代码、数据等的介质。这些制造品可以代表计算机可读介质(例如,任何适合形式的记忆器或存储器),所述介质可以存储处理器可执行代码,所述代码被处理器64使用以进行本公开的技术。通常,处理器64可以执行软件应用,所述软件应用包括处理地震数据的程序,所述地震数据是依据本文所描述的实施方式通过地震勘测的接收器获取的。
记忆器66和存储器68也可以用于存储数据、数据的分析、软件应用等。记忆器66和存储器68可以代表非暂时性计算机可读介质(例如,任何适合形式的记忆器或存储器),所述介质可以存储处理器可执行代码,所述代码被处理器64使用以进行多种本文所描述的技术。应当注意,非暂时性仅表明所述介质是有形的而不是信号。
I/O端口70可以是可以连接到其他外围部件例如输入装置(例如键盘、鼠标)、传感器、输入/输出(I/O)模块等的接口。I/O端口70可以使计算系统60能够与海洋勘测系统22、陆地勘测系统38等中的其他装置通过I/O端口70通信。
显示器72可以描绘与由处理器64处理的软件或可执行代码相关的可视化。在一个实施方式中,显示器72可以是能够从计算系统60的使用者处接收输入的触屏显示器。显示器72也可以用于观看和分析所获取的地震数据的分析结果,以确定地下区域26内的地质地层、地下区域26内的烃沉积层的位置和性质等。显示器72可以是任何适合类型的显示器,例如液晶显示器(LCD)、等离子显示器,或者例如有机发光二极管(OLED)显示器。除了通过显示器72描绘本文所描述的可视化之外,应当注意,计算系统60也可以通过其他有形元件例如纸(例如,通过打印)等来描绘可视化。
考虑到上述情况,本文描述的当前技术也可以用超级计算机进行,所述超级计算机应用了多个计算系统60、云计算系统等以在多个计算系统中分布要进行的处理。在这个情况下,作为超级计算机的一部分操作的每个计算系统60可以不包括被列作计算系统60的一部分的每个部件。例如,可能不是每个计算系统60都包括显示部件72,因为多个显示部件72对于被设计成连续处理地震数据的超级计算机来说可能是没有用的。
在进行多种类型的地震数据处理后,计算系统60可以将分析的结果存储在一个或多个数据库74中。数据库74可以通信连接到网络中,所述网络可以通过通信部件62将数据传送到计算系统60并从计算系统60接收数据。此外,数据库74可以存储关于地下区域26的信息,例如关于地下区域26的之前的地震记录图、地质样品数据、地震图像等。
尽管已经对于计算系统60讨论了上述的部件,但应当注意,类似的部件可以构成计算系统60。此外,计算系统60也可以是海洋勘测系统22或陆地勘测系统38的一部分,并因此可以监测和控制源32或40、接收器36、44、46等的某些操作。此外,应当注意,所列部件是作为实施例部件提供的,且本文所描述的实施方式不限于参考图4所描述的部件。
在某些实施方式中,计算系统60可以根据通过上述接收器接收的地震数据生成地下区域26的二维显示或三维显示。此外,可以将与多个源/接收器组合相关联的地震数据进行组合以创建地下区域26的近乎连续的可延伸一些距离的轮廓。在二维(2-D)地震勘测中,接收器的位置可以沿着单线放置,而在三维(3-D)勘测中,接收器的位置可以以网格图案分布于表面。因此,2-D地震勘测可以提供直接存在于记录位置下方的地球层的横截面图(纵切)。在另一方面,3-D地震勘测可以创建数据“立方体”或可以对应于地下区域26的3-D图的体积。
此外,4-D(或延时)地震勘测可以包括在3-D勘测期间在不同时间获取的地震数据。使用在不同时间获取的不同地震图像,计算系统60可以比较两个图像以鉴别地下区域26内的变化。
在任何情况下,地震勘测可以由非常大量的单独地震记录或踪迹组成。因此,计算系统60可以用于分析所获取的地震数据,以获得代表地下区域26的图像,并确定烃沉积层的位置和性质。为此,多种地震数据处理算法可以用于从所获取的地震数据中除去噪音、迁移预处理的地震数据、鉴别在多个地震图像之间的偏移、排列多个地震图像等。
在计算系统60分析所获取的地震数据之后,地震数据分析的结果(例如地震记录图、地震图像、地质地层地图等)可以被用于进行在烃勘探和生产工业内的多种操作。例如,如上所述,所获取的地震数据可以用于进行图1的方法10,所述方法详细描述了根据获取的地震数据的分析可以采用的多种处理。
如上所讨论的,计算系统60可以用于分析根据反射的地震波形52获取的地震数据,以获得代表地下区域26的图像,并确定烃沉积层的位置和性质。图5示出了通过用地震速度绘制地震图像而产生的一个区域(例如地下区域26)的图像78(例如速度模型),这是通过全波形反演(FWI)生成的。FWI地震数据处理通常通过使获取的地震数据与模拟(即建模)的地震数据之间的差异最小化来操作以创建地下(速度)模型估计。为了便于讨论图像78,图5包括了两个轴80和82,其分别代表了深度和沿着表面的距离(例如,沿着海底28的表面或地球的表面42的距离)。图像78包括变化的阴影,所述阴影是地震波传播通过地下区域的速度的指示:较深的阴影区域79代表相对较慢的地震波形的速度,且较浅的阴影区域81代表较快的地震波形的速度。例如,沿着线84从图像78的顶侧86向图像78的底侧88移动,图像的阴影通常变得更浅。因此,在更深的地下区域(例如,在图像78中沿着轴80从较深的阴影区域79到较浅的阴影区域81更加向下)中地震波形的速度倾向于增加。
地震波的速度通常随着深度增加而增加,这可能与地下区域内较大深度处增加的压力有关。此外,地震波的速度可能根据地下区域的地质地层的组成(例如,岩石、页岩、沙、岩石孔隙度、液体的存在和气体的存在)而变化。因此,不同的地震波速度是地层类型的指示并可以用于确定烃沉积层的位置(例如,地震波的速度减小而不是增加的地方)。
回来参考图5,图像78包括多个可能是烃沉积层位置的区域90a、90b、90c和90d。具体地,地震波的速度在区域90a、90b、90c和90d中减小。某些烃沉积层的区域比其他更明显(例如,地震波速度减小的量级更大)。此外,区域90a、90b、90c和90d可能基于地震波速度的变化和地震衰减产生“阴影”效应(例如,掩蔽区域92),所述效应可以掩蔽在区域90a、90b、90c和90d的下方可能存在的区域90,并阻止了烃沉积层的另一区域90被定位。此外,如本文所讨论的,速度的灵敏度可能随着深度增加而减小(例如,由于地层的组成的变化)。
本公开涉及用增加的灵敏度在地下区域内的较大深度处鉴定烃沉积层的技术。如上所讨论的,根据FWI速度模型生成图像78;然而图像78和更具体地所述分析可能不能有效地揭示例如地下区域内深处的烃沉积层和/或位于特定地层类型下的烃沉积层。在某些实施方式中,可以用第一速度模型(例如,FWI速度模型或简写为FWI速度)和地下区域的周边区域的物理属性(例如,指示在地下区域26内的地质地层的性质)的组合来生成属性模型。例如,FWI速度可以被写成多个物理参数的平均值(例如,算数平均值、调和平均值、几何平均值等),所述物理参数例如不同岩石地层的速度和/或包括不同流体的岩石的速度。可以从任何适合的方法,例如从井内的井下工具,获得所述参数。应当理解,关于如何得到地震速度的技术是非限制性的。也就是说,从FWI速度或任何其他方法得到的地震速度可以应用于所述技术中。
作为非限制性实施例,FWI速度可以表示为充满盐水的岩石的速度V盐水和充满烃的岩石的速度V的几何平均值,且通过基于速度的饱和度项SV对每个速度(例如,物理属性)进行加权:
VFWI=(V盐水)1-SV(V)SV
然而,应当注意,可以使用除了上述V盐水和V的几何平均值以外的其他平均方法,例如算数平均值、调和平均值等,或者可以将沉积速度拟合为V盐水、V和烃饱和度的函数的任何等式。在某些实施方式中,所述V盐水和V可以是来自测井日志的1维值(例如速度作为深度的函数)。在其他实施方式中,可以根据从适合的方法(例如,测井日志数据)观察或确定的趋势的外推或插值来确定每个速度。在其他实施方式中,V盐水和V可以是2D和/或3D的,这可以进一步提高本文所述技术的精确度。
在一个实施例中,可以重新排列等式以求解SV:
Figure BDA0002709150670000121
本领域普通技术人员应当理解,尽管烃饱和度的物理定义可能与SV不同,但SV属性可以与饱和度直接相关,并校准不同流体(例如,水、盐水、烃、气体)对在大部分深度处的地震波速度的影响。因此,本技术可以应用于现有基础设施附近的勘探。此外,SV的使用可以允许速度模型更好地说明在增加的储池深度处的烃,因为通常随着储池深度的增加,烃对于地震速度模型的影响变得微弱。以这种方式,利用SV作为地震速度模型的校准或调整,以便于用地震速度模型和测井日志数据创建SV(与烃饱和度相关的属性),以突出显示速度被烃影响的区域(例如,下文讨论的图6的区域90a、90b、90c、90d和90e)。
尽管上述SV等式是基于V盐水和V,但应当理解,所用的速度取决于所研究的地质地层。此外,根据SV等式生成在地下区域内特定深度的绘制图像可能是有利的。因此,可以确定更加适合于用SV等式来生成绘制图像的深度阈值。此外,在某些实施方式中,可以根据属性模型生成地震衰减模型或剪切波速度模型。
图6示出了根据本公开,基于图5所示的图像78,使用SV项生成的图像94。也就是说,图像94示出了与图像78相同的地下区域26。深阴影指示低SV值,其指示存在很少或不存在烃,而较浅区域指示烃。例如,图像94说明了图5的相同四个区域90a、90b、90c和90d,并且图6中的这些区域90a、90b、90c和90d指示烃沉积层(例如,包含烃的区域)。如前所述,区域90a、90b、90c和90d对应于图像78中示出的区域(例如,如图5所示),然而图6中的区域90a、90b、90c和90d更加明显(即相比于图像78,图像94中包含烃的区域的可视化更容易明显)。此外,图像94中说明了另一区域90e,而此区域90e与掩蔽区域92重叠(例如,如图5所示)。
在某些实施方式中,计算系统60可以提供区域90a、90b、90c、90d和90e的指示。例如,计算系统60可以确定包括具有高于阈值的指示烃存在的值的像素的图像94的子集。因此,计算系统60可以提供地下26中可能具有烃的区域的指示,并提供区域90a、90b、90c、90d和90e的位置(例如,坐标或相对面积),所述区域包含烃沉积层和/或(在某些情况下)地下钻探危害物。在其他实施方式中,根据当前技术生成的图像(例如图像94)可以通过计算系统60的显示器72来显示,由此帮助计算系统60的使用者定位区域90a、90b、90c、90d和90e。
在某些实施方式中,可以结合地震处理方案,例如图7中说明的方法98来生成图像94。如所说明的,方法98包括地震处理顺序,所述顺序包括步骤100中的地震数据收集、步骤102中的地震数据编辑、步骤104中的初始处理和步骤106中的信号处理、调节、构建速度模型(例如定义衰减模型、剪切波速度模型等)和成像(这可以例如包括产生成像的截面或体积),然后是地震数据的任何解释、任何其他与所需勘探目标一致的图像增强、从所处理的地震数据中生成属性、根据需要重新解释地震数据、以及确定和/或生成钻探计划或其他地震勘测应用。作为处理方法98的结果,可以鉴定到在地下区域26内的烃的位置(例如,在区域90a、90b、90c、90d和90e中)。可以根据改良的或以其它方式校准的FWI速度模型来定位烃沉积层中包含的这些烃。如前面所讨论的,本公开的一个实施方式涉及根据地震速度模型和与地下区域26有关的物理属性的组合来产生属性模型(包括SV)。在一个实施方式中,物理属性是地震波形穿过不同地质地层的速度,所述地质地层例如沙、页岩、充满盐水的岩石和/或充满烃的岩石,所述速度由其他适合的方法(例如,从钻井日志)确定。根据存在的不同地质地层,以及存在的地质地层的组成,速度模型可以表示为每个对地质地层测定的速度的平均值。此外,每个对地质地层测定的速度被给予与基于速度的饱和度SV成比例的权重,所述权重可以通过所测定的速度和地震速度模型来确定。因此,所述属性模型可以提高烃的检测,例如,在地下区域内的更大深度处。此外,上述实施方式也可以用于通过在陆地或海洋环境中使用地震波速来生成复杂的覆盖层的压力预测,这可以促进烃和地下危害物的检测。
上述具体实施方式已用实施例的方式示出,且应当理解,这些实施方式可能易于进行多种修改和替代形式。应当进一步理解,权利要求书并不旨在限制于所公开的特定形式,而是涵盖所有落入本公开的精神和范围内的修改、等同物和替代物。
本文所呈现和要求保护的技术被引用并应用于具有实际性质的实物和具体实施例,所述技术明显改进了现有技术领域且因此不是抽象的、无形的或纯理论的。此外,如果本说明书所附的任何权利要求包含被指定为“用于【进行】……【功能】的方法”或“用于【进行】……【功能】的步骤”的一个或多个要素,这样的要素应在35U.S.C.112(f)下理解。然而,对于包含以任何其他方式指定的要素的任何权利要求,这样的要素不应在35U.S.C.112(f)下理解。

Claims (20)

1.一种方法,所述方法包括:
根据地震波形通过速度模型构建器生成速度模型;
接收代表一个或多个测井周围区域的物理属性测量的数据;
根据所述速度模型和所述数据生成属性模型;以及
根据所述属性模型绘制图像,所述图像用于在包括烃储池并包含有助于烃的存在、迁移或积累的结构或地层特征的地下区域上方的地震勘探。
2.根据权利要求1所述的方法,其中所述数据包括在充满盐水的地下区域内地层的沉积速度。
3.根据权利要求2所述的方法,所述方法包括根据所述一个或多个测井的日志数据生成在充满盐水的地下区域内地层的沉积速度。
4.根据权利要求2所述的方法,其中所述数据包括在充满烃的地下区域内地层的沉积速度。
5.根据权利要求4所述的方法,所述方法包括根据所述测井的日志数据生成在充满烃的地下区域内地层的沉积速度。
6.根据权利要求1所述的方法,所述方法包括与所述速度模型构建器一起利用全波形反演。
7.根据权利要求1所述的方法,所述方法包括与所述速度模型构建器一起利用层析成像。
8.根据权利要求1所述的方法,其中所述属性模型与地下区域内地层的流体饱和度直接相关。
9.根据权利要求1所述的方法,其中所述属性模型包括经校准以反映地下区域内地层的流体饱和度对速度模型的值的影响的速度模型。
10.一种装置,所述装置包括:
处理器,所述处理器被配置为:
根据地震波形,通过地震速度模型构建器,利用全波形反演或层析成像,生成速度模型。
接收代表测井周围区域的第一物理属性测量的第一数据和代表测井周围区域的第二物理属性测量的第二数据;
根据所述速度模型、第一数据和第二数据确定校准因子;以及
根据所述校准因子生成图像,所述图像用于在包括烃储池并包含有助于烃的存在、迁移或积累的结构或地层特征的地下区域上方的地震勘探。
11.根据权利要求10所述的装置,其中所述第一数据包括在充满盐水的地下区域内地层的沉积速度。
12.根据权利要求11所述的装置,其中所述第二数据包括在充满烃的地下区域内地层的沉积速度。
13.根据权利要求12所述的装置,其中所述处理器被配置为对第一数据应用第一权重值以生成加权的第一数据并对第二数据应用第二权重值以生成加权的第二数据。
14.根据权利要求13所述的装置,其中所述处理器被配置为使速度模型的速度等于第一加权数据和第二加权数据。
15.根据权利要求14所述的装置,其中所述处理器被配置为确定作为第一权重值中的项和第二权重值中的项的校准因子,所述校准因子是通过使速度模型中的速度等于第一加权数据和第二加权数据而获得的。
16.根据权利要求10所述的装置,其中所述校准因子与地下区域内地层的流体饱和度直接相关。
17.根据权利要求10所述的装置,其中所述校准因子反映了地下区域内地层的流体饱和度对速度模型的值的影响。
18.一种或多种有形的、非暂时性的、机器可读的介质,所述介质包含指令,所述指令被配置为引起处理器:
根据地震波形通过地震模型构建器产生速度模型;
接收代表测井周围区域的物理属性测量的数据;
根据所述速度模型和所述数据生成属性模型;以及
根据所述属性模型绘制代表地球地下区域内的烃的地震图像。
19.根据权利要求18所述的一种或多种机器可读的介质,其中所述指示包括地球的地下区域内烃的位置。
20.根据权利要求18所述的一种或多种机器可读的介质,所述介质包括指令,所述指令被配置为引起处理器传输指示以显示在显示器上。
CN201980023710.7A 2018-03-30 2019-03-28 地震速度衍生的烃指示 Pending CN111971586A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201862650667P 2018-03-30 2018-03-30
US62/650,667 2018-03-30
PCT/US2019/024486 WO2019191368A1 (en) 2018-03-30 2019-03-28 Seismic velocity derived hydrocarbon indication

Publications (1)

Publication Number Publication Date
CN111971586A true CN111971586A (zh) 2020-11-20

Family

ID=66429470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980023710.7A Pending CN111971586A (zh) 2018-03-30 2019-03-28 地震速度衍生的烃指示

Country Status (9)

Country Link
US (1) US11112515B2 (zh)
EP (1) EP3776013B1 (zh)
CN (1) CN111971586A (zh)
AU (1) AU2019243973A1 (zh)
BR (1) BR112020019836B1 (zh)
CA (1) CA3094120C (zh)
EA (1) EA202092277A1 (zh)
MX (1) MX2020010107A (zh)
WO (1) WO2019191368A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11112515B2 (en) * 2018-03-30 2021-09-07 Bp Corporation North America Inc. Seismic velocity derived hydrocarbon indication

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020028044A1 (en) * 2018-07-31 2020-02-06 Exxonmobil Upstream Research Company Detecting fluid types using petrophysical inversion

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7768870B2 (en) * 2007-06-29 2010-08-03 Institut Francais Du Petrole Method for adjusting a seismic wave velocity model according to information recorded in wells
CN102272631A (zh) * 2009-01-09 2011-12-07 埃克森美孚上游研究公司 用被动地震数据的烃探测
CN103703391A (zh) * 2011-03-30 2014-04-02 埃克森美孚上游研究公司 使用频谱整形的全波场反演的收敛速度
CN104570125A (zh) * 2014-09-26 2015-04-29 郭平 一种利用井数据提高成像速度模型精度的方法
CN106415321A (zh) * 2014-06-23 2017-02-15 雪佛龙美国公司 用于储层建模的基于瞬时等时属性的地质体识别
CN106461810A (zh) * 2014-06-05 2017-02-22 雪佛龙美国公司 构建复杂地球模型的系统和方法
CN106526669A (zh) * 2016-09-19 2017-03-22 中国石油化工股份有限公司 一种页岩油气藏的地震储层预测方法
CN106959467A (zh) * 2017-03-20 2017-07-18 中国科学院地质与地球物理研究所 地震波速度反演方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11112515B2 (en) * 2018-03-30 2021-09-07 Bp Corporation North America Inc. Seismic velocity derived hydrocarbon indication

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7768870B2 (en) * 2007-06-29 2010-08-03 Institut Francais Du Petrole Method for adjusting a seismic wave velocity model according to information recorded in wells
CN102272631A (zh) * 2009-01-09 2011-12-07 埃克森美孚上游研究公司 用被动地震数据的烃探测
CN103703391A (zh) * 2011-03-30 2014-04-02 埃克森美孚上游研究公司 使用频谱整形的全波场反演的收敛速度
CN106461810A (zh) * 2014-06-05 2017-02-22 雪佛龙美国公司 构建复杂地球模型的系统和方法
CN106415321A (zh) * 2014-06-23 2017-02-15 雪佛龙美国公司 用于储层建模的基于瞬时等时属性的地质体识别
CN104570125A (zh) * 2014-09-26 2015-04-29 郭平 一种利用井数据提高成像速度模型精度的方法
CN106526669A (zh) * 2016-09-19 2017-03-22 中国石油化工股份有限公司 一种页岩油气藏的地震储层预测方法
CN106959467A (zh) * 2017-03-20 2017-07-18 中国科学院地质与地球物理研究所 地震波速度反演方法和装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11112515B2 (en) * 2018-03-30 2021-09-07 Bp Corporation North America Inc. Seismic velocity derived hydrocarbon indication

Also Published As

Publication number Publication date
US11112515B2 (en) 2021-09-07
BR112020019836B1 (pt) 2023-03-14
CA3094120C (en) 2023-09-19
WO2019191368A1 (en) 2019-10-03
EP3776013B1 (en) 2024-02-21
CA3094120A1 (en) 2019-10-03
EP3776013A1 (en) 2021-02-17
EA202092277A1 (ru) 2020-12-28
AU2019243973A1 (en) 2020-10-08
MX2020010107A (es) 2021-02-18
BR112020019836A2 (pt) 2021-01-05
US20190302288A1 (en) 2019-10-03

Similar Documents

Publication Publication Date Title
CN112703429B (zh) 基于机器学习的地震属性分析
CA2920499C (en) Stratigraphic function
EP2846175B1 (en) Seismic survey analysis
CA2940406C (en) Characterizing a physical structure using a multidimensional noise model to attenuate noise data
EP3710867B1 (en) Noise attenuation of multiple source seismic data
CN113015926A (zh) 无源地震成像
WO2017160273A1 (en) Dip estimation via a modified structure tensor
CA3094120C (en) Seismic velocity derived hydrocarbon indication
CA3131518C (en) Iterative stochastic seismic inversion
EP3963371B1 (en) 4d time shift and amplitude joint inversion for velocity perturbation
CN111936888B (zh) 用于倾斜正交晶介质的波场传播子
US20220236435A1 (en) Low-Frequency Seismic Survey Design
Penna et al. Geostatistical seismic inversion and 3D modelling of metric flow units, porosity and permeability in Brazilian presalt reservoir
US20240111072A1 (en) Method and Apparatus for Petrophysical Classification, Characterization, and Uncertainty Estimation
US11880008B2 (en) Velocity model construction
EA042380B1 (ru) Обнаружение признаков углеводородов по сейсмическим скоростям
EA043508B1 (ru) Итеративная стохастическая сейсмическая инверсия
EA044564B1 (ru) Построение модели скорости
EA043815B1 (ru) Проектирование низкочастотной сейсмической съемки

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