CN112305597A - 储层预测方法及装置 - Google Patents

储层预测方法及装置 Download PDF

Info

Publication number
CN112305597A
CN112305597A CN201910682190.2A CN201910682190A CN112305597A CN 112305597 A CN112305597 A CN 112305597A CN 201910682190 A CN201910682190 A CN 201910682190A CN 112305597 A CN112305597 A CN 112305597A
Authority
CN
China
Prior art keywords
frequency
data
target layer
determining
low
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.)
Granted
Application number
CN201910682190.2A
Other languages
English (en)
Other versions
CN112305597B (zh
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.)
China National Petroleum Corp
BGP Inc
CNPC International Exploration and Production Co Ltd
Original Assignee
China National Petroleum Corp
BGP Inc
CNPC International Exploration and Production Co Ltd
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 China National Petroleum Corp, BGP Inc, CNPC International Exploration and Production Co Ltd filed Critical China National Petroleum Corp
Priority to CN201910682190.2A priority Critical patent/CN112305597B/zh
Publication of CN112305597A publication Critical patent/CN112305597A/zh
Application granted granted Critical
Publication of CN112305597B publication Critical patent/CN112305597B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

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

本发明公开了一种储层预测方法及装置，该方法包括根据目的层地震数据的振幅频谱确定有效频带及峰值频率，对目的层地震数据进行分解得到多频地震数据，将多频地震数据的频率信息按照入射角设置得到叠前AVO道集；确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据；根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；提取目的层流体属性数据的储层预测参数确定目的层的储层分布。本发明分别确定低频频率梯度数据及高频频率梯度数据，进而确定目的层流体属性数据，以实现对目的层的储层预测。据此得到的目的层流体属性数据包含多个频率，可以极大的提高储层预测的精度。

Description

储层预测方法及装置
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及储层预测方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
随着三维地震勘探技术的不断进步和三维地震资料在油田勘探开发领域应用的不断深入，地震属性解释的重要性日益提高。地震属性主要包括时间、振幅、频率、相位、相干和衰减几大类属性。在油田勘探开发领域，地震属性主要应用于油藏储层厚度、孔隙度、渗透率等储层参数的预测，其中地震资料的振幅、频率等单项地震属性的应用比较广泛。
流体活动属性(又称为流度属性)，是一种逐渐广泛应用的地震属性，主要是通过低频段的两个频率计算低频的频谱梯度来实现，但两个频率的选择需要地质专家具有较好的经验，导致现有的储层预测存在不够客观的缺陷。另外，通过地质专家依靠经验选择该两个低频频率，将严重降低储层预测的精度。
因此，现有的储层预测存在预测结果不客观及预测结果精度低的问题。
发明内容
本发明实施例提供一种储层预测方法，用以避免储层预测结果的不客观性，提高储层预测结果的精度，该方法包括：
在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率；
在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据；
将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集；
分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据；
根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度；
在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
本发明实施例还提供一种储层预测装置，用以避免储层预测结果的不客观性，提高储层预测结果的精度，该装置包括：
频带频率确定模块，用于在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率；
分解模块，用于在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据；
道头置换模块，用于将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集；
梯度数据确定模块，用于分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据；
流体属性数据确定模块，用于根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度；
储层预测模块，用于在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述储层预测方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述储层预测方法的计算机程序。
本发明实施例中，将多频地震数据的频率信息按照入射角设置以进行道头置换得到叠前AVO道集，分别确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据，进而根据低频频率梯度数据及高频频率梯度数据确定目的层流体属性数据，以实现对目的层的储层预测。据此得到的目的层流体属性数据包含低频和高频的多个频率，可以极大的提高储层预测的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例提供的储层预测方法的实现流程图；
图2为本发明实施例提供的储层预测方法中步骤101的实现流程图；
图3为本发明实施例提供的储层预测方法中步骤104的实现流程图；
图4为本发明实施例提供的储层预测装置的功能模块图；
图5为本发明实施例提供的储层预测装置中频带频率确定模块401的结构框图；
图6为本发明实施例提供的储层预测装置中梯度数据确定模块404的结构框图；
图7为本发明实施例提供的实际储层地震数据的合成地震记录标定剖面示意图；
图8为本发明实施例提供的实际储层地震数据的振幅频谱示意图；
图9为本发明实施例提供的多频地震数据的多道时频分解剖面示意图；
图10为本发明实施例提供的低频频率梯度数据剖面示意图；
图11为本发明实施例提供的高频频率梯度数据剖面示意图；
图12为本发明实施例提供的目的层流体属性数据剖面示意图；
图13为本发明实施例提供的从目的层流体属性数据中提取的流体属性平面示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1示出了本发明实施例提供的储层预测方法的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
如图1所示，储层预测方法，其包括：
步骤101，在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率；
步骤102，在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据；
步骤103，将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集；
步骤104，分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据；
步骤105，根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度；
步骤106，在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
在本发明实施例中，目的层是钻井的目的层，即所要研究的对象或者区域。可以通过对目的层地震数据进行井震时间深度匹配，确定目的层地震数据在时间域对应的时窗范围。进而在该时窗范围内根据目的层地震数据的振幅频谱，确定目的层地震数据振幅频谱的有效频带；以及根据目的层地震数据的振幅频谱，确定目的层地震数据振幅频谱的峰值频率。在本发明的一实施例中，该目的层地震数据以地震数据体的形式体现。
在根据目的层地震数据的振幅频谱确定其有效频带后，在有效频带内以预先设置的频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据。其中，在对目的层地震数据进行分解时，根据有效频带及预设频段间隔，将目的层地震数据在过完备时频原子库中进行分解，将目的层地震数据自适应的表示为匹配原子的线性组合形式，在三维地震数据的基础上增加频率轴，进而得到四维的多频地震数据。在本发明的一实施例中，多频地震数据以多频地震数据体的形式体现。
在本发明的一实施例中，预设频段间隔为预先设定的频段间隔，本领域技术人员可以理解的是，可以根据实际情况和具体需求预先设置该预设频段间隔。例如，将预设频段间隔预先设置为5赫兹，本领域技术人员可以理解的是，还可以将该预设频段间隔预先设置为除上述5赫兹之外的其他频率间隔，例如，2赫兹、4赫兹、8赫兹或者10赫兹等，本发明实施例对此不做特别的限制。
在得到包含频率信息的多频地震数据后，对多频地震数据进行道头置换。具体的，在对多频地震数据进行道头置换时，此时道头附加频率信息，按照叠前道集格式对多频地震数据进行加载，将多频地震数据的频率信息按照入射角设置，以模拟叠前AVO道集，即得到目的层地震数据对应的叠前AVO道集。
在得到目的层地震数据对应的叠前AVO道集后，以目的层地震数据振幅频谱的峰值频率为门槛值，目的层地震数据对应的叠前AVO道集包括两个部分：频率小于峰值频率的低频叠前AVO道集，以及频率不小于峰值频率的高频叠前AVO道集。此时，基于低频叠前AVO道集及高频叠前AVO道集，分别确定低频叠前AVO道集对应的低频频率梯度数据，以及高频叠前AVO道集对应的高频频率梯度数据。
在得到低频叠前AVO道集对应的低频频率梯度数据，以及高频叠前AVO道集对应的高频频率梯度数据后，进而可以根据低频频率梯度数据和高频频率梯度数据的差异度参数确定目的层流体属性数据。低频频率梯度数据和高频频率梯度数据的差异度参数，反映了低频频率梯度数据和高频频率梯度数据差异程度。在本发明的一实施例中，该目的层流体属性数据以地震数据体的形式体现。
在本发明的一实施例中，低频频率梯度数据和高频频率梯度数据的差异度参数至少包括低频频率梯度数据和高频频率梯度数据的差值。即在确定目的层流体属性数据时，具体的，将低频频率梯度数据和高频频率梯度数据的差值，确定为目的层流体属性数据。另外，还可以将高频频率梯度数据和低频频率梯度数据的差值，确定为目的层流体属性数据。
在得到目的层流体属性数据后，在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，以对目的层进行储层预测。此时得到的目的层流体属性数据包含低频和高频的多个频率，可以极大的提高储层预测的精度。
在本发明的一实施例中，目的层流体属性数据的储层预测参数包括目的层流体属性数据的均方根振幅。具体的，在进行储层预测时，在目的层地震数据对应的时窗范围内提取目的层流体属性数据的均方根振幅，以实现对目的层的储层预测。
在本发明实施例中，将多频地震数据的频率信息按照入射角设置以进行道头置换得到叠前AVO道集，分别确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据，进而根据低频频率梯度数据及高频频率梯度数据确定目的层流体属性数据，以实现对目的层的储层预测。据此得到的目的层流体属性数据包含低频和高频的多个频率，可以极大的提高储层预测的精度。
图2示出了本发明实施例提供的储层预测方法中步骤101的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，如图2所示，步骤101，在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率，包括：
步骤201，利用测井数据及地震子波确定合成地震记录；
步骤202，对合成地震记录进行层位标定，确定目的层地震数据对应的时窗范围；
步骤203，在时窗范围内提取目的层地震数据的振幅频谱；
步骤204，根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率。
测井数据是指在测井是形成的能够反映出不同岩性、层位特征的数据。在本发明实施例中，测井数据主要包括声波时差测井曲线、密度测井曲线等。地震子波是一段具有确定的起始时间、能量有限且具有一定延续长度的信号，它是地震记录中的基本单元。在确定合成地震记录时，主要通过勘探区域已知井的钻井分层、声波时差测井曲线、密度测井曲线及地震子波激发制作合成地震记录。
在激发制作得到合成地震记录后，对合成地震记录进行层位标定，将钻井的目的层的分层在地震数据上标定的位置确定为目的层。利用井震时间深度匹配，用地震剖面中目的层对应的时间范围作为目的层地震数据属性提取的时间窗口，确定目的层地震数据对应的时窗范围。
在确定目的层地震数据对应的时窗范围后，在上述时窗范围内提取目的层地震数据的振幅频谱，进而通过目的层地震数据的振幅频谱，确定目的层地震数据振幅频谱的有效频带(即目的层地震数据振幅频谱的带宽)及峰值频率。
在本发明实施例中，对合成地震记录进行层位标定确定目的层地震数据对应的时窗范围，在时窗范围内提取目的层地震数据的振幅频谱，进而根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率，可以进一步提高储层预测的精度。
图3示出了本发明实施例提供的储层预测方法中步骤104的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，如图3所示，步骤104，分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据，包括：
步骤301，在时窗范围内获取频率小于频率峰值的低频叠前AVO道集的AVO属性，确定低频频率梯度数据；
步骤302，在时窗范围内获取频率不小于频率峰值的高频叠前AVO道集的AVO属性，确定高频频率梯度数据。
在确定叠前AVO道集对应的频率梯度数据时，分别获取低频叠前AVO道集的频率梯度数据及高频叠前AVO道集的频率梯度数据。具体的，以峰值频率为门槛值，在时窗范围内计算、获取低频叠前AVO道集的AVO属性，进而确定低频叠前AVO道集对应的低频频率梯度数据。同样以峰值频率为门槛值，在时窗范围内计算、获取高频叠前AVO道集的AVO属性，进而确定高频叠前AVO道集对应的低频频率梯度数据。至此可以分别确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据。
在本发明实施例中，分别在时窗范围内获取低频叠前AVO道集对应的低频频率梯度数据，及在时窗范围内获取高频叠前AVO道集对应的高频频率梯度数据，可以在提高储层预测精度的同时，提高储层预测的效率。
本发明实施例中还提供了一种储层预测装置，如下面的实施例所述。由于这些装置解决问题的原理与储层预测方法相似，因此这些装置的实施可以参见方法的实施，重复之处不再赘述。
图4示出了本发明实施例提供的储层预测装置的功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
参考图4，所述储层预测装置所包含的各个模块用于执行图1对应实施例中的各个步骤，具体请参阅图1以及图1对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述储层预测装置包括频带频率确定模块401、分解模块402、道头置换模块403、梯度数据确定模块404、流体属性数据确定模块405及储层预测模块406。
频带频率确定模块401，用于在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率。
分解模块402，用于在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据。
道头置换模块403，用于将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集。
梯度数据确定模块404，用于分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据。
流体属性数据确定模块405，用于根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度。
储层预测模块406，用于在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
在本发明的一实施例中，低频频率梯度数据和高频频率梯度数据的差异度参数至少包括低频频率梯度数据和高频频率梯度数据的差值。即流体属性数据确定模块405在确定目的层流体属性数据时，具体的，将低频频率梯度数据和高频频率梯度数据的差值，确定为目的层流体属性数据。另外，还可以将高频频率梯度数据和低频频率梯度数据的差值，确定为目的层流体属性数据。
在本发明的一实施例中，目的层流体属性数据的储层预测参数包括目的层流体属性数据的均方根振幅。具体的，在进行储层预测时，储层预测模块406在目的层地震数据对应的时窗范围内提取目的层流体属性数据的均方根振幅，以实现对目的层的储层预测。
在本发明实施例中，道头置换模块403将多频地震数据的频率信息按照入射角设置以进行道头置换得到叠前AVO道集，梯度数据确定模块404分别确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据，进而流体属性数据确定模块405根据低频频率梯度数据及高频频率梯度数据确定目的层流体属性数据，储层预测模块406以实现对目的层的储层预测。据此得到的目的层流体属性数据包含低频和高频的多个频率，可以极大的提高储层预测的精度。
图5示出了本发明实施例提供的储层预测装置中频带频率确定模块401的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，参考图5，所述频带频率确定模块所包含的各个单元用于执行图2对应实施例中的各个步骤，具体请参阅图2以及图2对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述频带频率确定模块包括合成地震记录确定单元501、层位标定单元502、振幅频谱提取单元503及频带频率确定单元504。
合成地震记录确定单元501，用于利用测井数据及地震子波确定合成地震记录；
层位标定单元502，用于对合成地震记录进行层位标定，确定目的层地震数据对应的时窗范围；
振幅频谱提取单元503，用于在时窗范围内提取目的层地震数据的振幅频谱；
频带频率确定单元504，用于根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率。
在本发明实施例中，层位标定单元502对合成地震记录进行层位标定确定目的层地震数据对应的时窗范围，振幅频谱提取单元503在时窗范围内提取目的层地震数据的振幅频谱，进而频带频率确定单元504根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率，可以进一步提高储层预测的精度。
图6示出了本发明实施例提供的储层预测装置中梯度数据确定模块404的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
在本发明的一实施例中，参考图6，所述梯度数据确定模块404所包含的各个单元用于执行图3对应实施例中的各个步骤，具体请参阅图3以及图3对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述梯度数据确定模块404包括低频频率梯度数据确定单元601和高频频率梯度数据确定单元602。
低频频率梯度数据确定单元601，用于在时窗范围内获取频率小于频率峰值的低频叠前AVO道集的AVO属性，确定低频频率梯度数据。
高频频率梯度数据确定单元602，用于在时窗范围内获取频率不小于频率峰值的高频叠前AVO道集的AVO属性，确定高频频率梯度数据。
在本发明实施例中，低频频率梯度数据确定单元601在时窗范围内获取低频叠前AVO道集对应的低频频率梯度数据，及高频频率梯度数据确定单元602在时窗范围内获取高频叠前AVO道集对应的高频频率梯度数据，可以在提高储层预测精度的同时，提高储层预测的效率。
图7示出了本发明实施例提供的实际储层地震数据的合成地震记录标定剖面示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
根据图7对合成地震记录进行的层位标定，确定的目的层在地震剖面中的时窗范围约为1580ms-1700ms，约为120ms的时间范围。
图8示出了本发明实施例提供的实际储层地震数据的振幅频谱示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
如图8所示，实际储层地震数据的振幅频谱的有效频带为5赫兹-30赫兹，峰值频率约为16赫兹，作为后续的频率门槛值。
图9示出了本发明实施例提供的多频地震数据的多道时频分解剖面示意，图10示出了本发明实施例提供的低频频率梯度数据剖面示意，图11示出了本发明实施例提供的高频频率梯度数据剖面示意，图12示出了本发明实施例提供的目的层流体属性数据剖面示意，图13示出了本发明实施例提供的从目的层流体属性数据中提取的流体属性平面示意。
鉴于据此得到的实际储层地震数据对应的流体属性数据包含了低频和高频的多个频率，因此，能够有效的提高储层预测的精度。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述储层预测方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述储层预测方法的计算机程序。
综上所述，本发明实施例中，将多频地震数据的频率信息按照入射角设置以进行道头置换得到叠前AVO道集，分别确定低频叠前AVO道集对应的低频频率梯度数据及高频叠前AVO道集对应的高频频率梯度数据，进而根据低频频率梯度数据及高频频率梯度数据确定目的层流体属性数据，以实现对目的层的储层预测。据此得到的目的层流体属性数据包含低频和高频的多个频率，可以极大的提高储层预测的精度。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种储层预测方法，其特征在于，包括：
在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率；
在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据；
将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集；
分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据；
根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度；
在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
2.如权利要求1所述的储层预测方法，其特征在于，低频频率梯度数据和高频频率梯度数据的差异度参数至少包括低频频率梯度数据和高频频率梯度数据的差值。
3.如权利要求1所述的储层预测方法，其特征在于，目的层流体属性数据的储层预测参数包括目的层流体属性数据的均方根振幅。
4.如权利要求1所述的储层预测方法，其特征在于，在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率，包括：
利用测井数据及地震子波确定合成地震记录；
对合成地震记录进行层位标定，确定目的层地震数据对应的时窗范围；
在时窗范围内提取目的层地震数据的振幅频谱；
根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率。
5.如权利要求1所述的储层预测方法，其特征在于，分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据，包括：
在时窗范围内获取频率小于频率峰值的低频叠前AVO道集的AVO属性，确定低频频率梯度数据；
在时窗范围内获取频率不小于频率峰值的高频叠前AVO道集的AVO属性，确定高频频率梯度数据。
6.一种储层预测装置，其特征在于，包括：
频带频率确定模块，用于在目的层地震数据对应的时窗范围内，根据目的层地震数据的振幅频谱确定目的层地震数据的有效频带及峰值频率；
分解模块，用于在有效频带内以预设频段间隔对目的层地震数据进行分解，得到包含频率信息的多频地震数据；
道头置换模块，用于将多频地震数据的频率信息按照入射角设置以进行道头置换，得到目的层地震数据对应的叠前AVO道集；
梯度数据确定模块，用于分别确定频率小于峰值频率的低频叠前AVO道集对应的低频频率梯度数据及频率不小于峰值频率的高频叠前AVO道集对应的高频频率梯度数据；
流体属性数据确定模块，用于根据低频频率梯度数据和高频频率梯度数据的差异度参数，确定目的层流体属性数据；差异度参数反映低频频率梯度数据和高频频率梯度数据差异程度；
储层预测模块，用于在目的层地震数据对应的时窗范围内提取目的层流体属性数据的储层预测参数，确定目的层的储层分布。
7.如权利要求6所述的储层预测装置，其特征在于，低频频率梯度数据和高频频率梯度数据的差异度参数至少包括低频频率梯度数据和高频频率梯度数据的差值。
8.如权利要求6所述的储层预测装置，其特征在于，目的层流体属性数据的储层预测参数包括目的层流体属性数据的均方根振幅。
9.如权利要求6所述的储层预测装置，其特征在于，频带频率确定模块包括：
合成地震记录确定单元，用于利用测井数据及地震子波确定合成地震记录；
层位标定单元，用于对合成地震记录进行层位标定，确定目的层地震数据对应的时窗范围；
振幅频谱提取单元，用于在时窗范围内提取目的层地震数据的振幅频谱；
频带频率确定单元，用于根据目的层地震数据的振幅频谱确定目的层地震数据振幅频谱的有效频带及峰值频率。
10.如权利要求6所述的储层预测装置，其特征在于，梯度数据确定模块包括：
低频频率梯度数据确定单元，用于在时窗范围内获取频率小于频率峰值的低频叠前AVO道集的AVO属性，确定低频频率梯度数据；
高频频率梯度数据确定单元，用于在时窗范围内获取频率不小于频率峰值的高频叠前AVO道集的AVO属性，确定高频频率梯度数据。
11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
CN201910682190.2A 2019-07-26 2019-07-26 储层预测方法及装置 Active CN112305597B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910682190.2A CN112305597B (zh) 2019-07-26 2019-07-26 储层预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910682190.2A CN112305597B (zh) 2019-07-26 2019-07-26 储层预测方法及装置

Publications (2)

Publication Number Publication Date
CN112305597A true CN112305597A (zh) 2021-02-02
CN112305597B CN112305597B (zh) 2024-09-27

Family

ID=74328773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910682190.2A Active CN112305597B (zh) 2019-07-26 2019-07-26 储层预测方法及装置

Country Status (1)

Country Link
CN (1) CN112305597B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114402A (zh) * 2021-12-21 2022-03-01 中国石油大学(北京) 储层流体检测方法、装置和电子设备

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6058073A (en) * 1999-03-30 2000-05-02 Atlantic Richfield Company Elastic impedance estimation for inversion of far offset seismic sections
US20020185329A1 (en) * 2001-04-30 2002-12-12 The Regents Of The University Of California Frequency-dependent processing and interpretation (FDPI) of seismic data for identifiying, imaging and monitoring fluid-saturated underground reservoirs
US20040122596A1 (en) * 2002-12-19 2004-06-24 Core Laboratories, Inc. Method for high frequency restoration of seismic data
US20050075790A1 (en) * 2003-09-24 2005-04-07 Taner M. Turhan Method for processing seismic data to identify anomalous absorption zones
US20090132169A1 (en) * 2007-11-19 2009-05-21 Schlumberger Technology Corporation Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements
US20110083844A1 (en) * 2008-06-18 2011-04-14 Oppert Shauna K Seismic Fluid Prediction Via Expanded AVO Anomalies
CN102109611A (zh) * 2009-12-23 2011-06-29 中国石油天然气集团公司 一种利用低频地震属性预测油藏优质储层的方法
CN103728654A (zh) * 2012-10-16 2014-04-16 中国石油化工股份有限公司 一种碳酸盐岩油气储层预测的方法
CN104391324A (zh) * 2014-12-03 2015-03-04 成都理工大学 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术
CN104977612A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 一种利用地震能量差异属性预测储层的方法
CN106405642A (zh) * 2016-07-19 2017-02-15 西安石油大学 基于脱压实声波速度的地震反演储层预测方法
CN106405639A (zh) * 2015-07-30 2017-02-15 中国石油化工股份有限公司 一种叠前地震储层岩性参数的反演方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6058073A (en) * 1999-03-30 2000-05-02 Atlantic Richfield Company Elastic impedance estimation for inversion of far offset seismic sections
US20020185329A1 (en) * 2001-04-30 2002-12-12 The Regents Of The University Of California Frequency-dependent processing and interpretation (FDPI) of seismic data for identifiying, imaging and monitoring fluid-saturated underground reservoirs
US20040122596A1 (en) * 2002-12-19 2004-06-24 Core Laboratories, Inc. Method for high frequency restoration of seismic data
US20050075790A1 (en) * 2003-09-24 2005-04-07 Taner M. Turhan Method for processing seismic data to identify anomalous absorption zones
US20090132169A1 (en) * 2007-11-19 2009-05-21 Schlumberger Technology Corporation Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements
US20110083844A1 (en) * 2008-06-18 2011-04-14 Oppert Shauna K Seismic Fluid Prediction Via Expanded AVO Anomalies
CN102109611A (zh) * 2009-12-23 2011-06-29 中国石油天然气集团公司 一种利用低频地震属性预测油藏优质储层的方法
CN103728654A (zh) * 2012-10-16 2014-04-16 中国石油化工股份有限公司 一种碳酸盐岩油气储层预测的方法
CN104977612A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 一种利用地震能量差异属性预测储层的方法
CN104391324A (zh) * 2014-12-03 2015-03-04 成都理工大学 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术
CN106405639A (zh) * 2015-07-30 2017-02-15 中国石油化工股份有限公司 一种叠前地震储层岩性参数的反演方法
CN106405642A (zh) * 2016-07-19 2017-02-15 西安石油大学 基于脱压实声波速度的地震反演储层预测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
XIAODONG WEI等: "SEG Technical Program Expanded Abstracts", THE GRADIENT OF THE AMPLITUDE SPECTRUM OF SEISMIC DATA AND ITS APPLICATION IN RESERVOIR PREDICTION, 14 October 2009 (2009-10-14), pages 1102 - 1106 *
王大兴;赵玉华;周义军;王永刚;强敏;程思检;童庆佳;: "苏里格气田多波地震处理与储层预测技术研究及应用", 中国石油勘探, vol. 16, no. 1, 30 September 2011 (2011-09-30), pages 95 - 102 *
祝旭双;龚福华;张世荣;夏训文;: "叠前弹性参数反演在鄂尔多斯Su76区块储层预测中的应用", 长江大学学报(自然科学版), vol. 08, no. 12, 31 December 2011 (2011-12-31), pages 43 - 45 *
董银萍等: "基于匹配追踪分解的流体活动因子预测煤层气甜点区", 煤田地质与勘探, vol. 46, no. 05, 31 July 2018 (2018-07-31), pages 90 - 96 *
谭辉煌;徐明华;彭宇;周长友;: "基于S变换的分频解释技术在三维地震资料中的运用", 科技传播, no. 05, 31 March 2011 (2011-03-31), pages 173 - 174 *
陈志刚等: ""两宽一高"地震资料的敏感方位油气检测技术在乍得潜山油藏描述的应用", 地球物理学进展, vol. 32, no. 03, 30 June 2017 (2017-06-30), pages 1114 - 1120 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114402A (zh) * 2021-12-21 2022-03-01 中国石油大学(北京) 储层流体检测方法、装置和电子设备

Also Published As

Publication number Publication date
CN112305597B (zh) 2024-09-27

Similar Documents

Publication Publication Date Title
CN105353407B (zh) 一种叠后地震波阻抗反演方法
US8923093B2 (en) Determining the quality of a seismic inversion
CN109061764B (zh) 一种分频融合波阻抗反演方法
US20150168573A1 (en) Geologic quality factor inversion method
WO2009067330A2 (en) Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements
CN109188520B (zh) 薄储层厚度预测方法及装置
WO2016135505A1 (en) A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation
RU2631407C1 (ru) Способ и устройство для обработки сейсмических сигналов
US10087733B2 (en) Fracture mapping using vertical seismic profiling wave data
EP2075597A2 (en) Spectral conditioning for surface seismic data
CN110095811B (zh) 膏岩层速度模型构建处理方法及装置
CN106707344A (zh) 一种地层层序划分方法及装置
CN102928875B (zh) 基于分数阶傅里叶域的子波提取方法
CN112305597B (zh) 储层预测方法及装置
CN113391353B (zh) 地震数据的处理方法及装置
CN112147677B (zh) 油气储层参数标签数据生成方法及装置
CN112649887A (zh) 基于钻井资料定量划分层序地层的方法及装置
CN114114421B (zh) 基于深度学习的导向自学习地震数据去噪方法及装置
CN112130209B (zh) 岩溶储层预测方法及装置
CN113589365B (zh) 基于时频域信息的储层尖灭线描述方法
CN112649862B (zh) 基于地层结构信息分离的断溶体识别方法和装置
CN114563816A (zh) 油气藏评价阶段建立地震解释速度模型的方法及装置
CN112379429A (zh) 地震数据的振幅补偿方法及装置
CN105989240A (zh) 砂层预测方法及装置
CN112305616A (zh) 光纤井中地震数据剖面的获取方法及装置

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