CN116047592A - 一种三维地震沿层属性切片定量化色标的校正方法 - Google Patents

一种三维地震沿层属性切片定量化色标的校正方法 Download PDF

Info

Publication number
CN116047592A
CN116047592A CN202111265039.2A CN202111265039A CN116047592A CN 116047592 A CN116047592 A CN 116047592A CN 202111265039 A CN202111265039 A CN 202111265039A CN 116047592 A CN116047592 A CN 116047592A
Authority
CN
China
Prior art keywords
seismic
attribute
reservoir
limit value
slice
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
CN202111265039.2A
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.)
China Petroleum and Chemical Corp
Sinopec North China Oil and Gas Co
Original Assignee
China Petroleum and Chemical Corp
Sinopec North China Oil and Gas Co
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 Petroleum and Chemical Corp, Sinopec North China Oil and Gas Co filed Critical China Petroleum and Chemical Corp
Priority to CN202111265039.2A priority Critical patent/CN116047592A/zh
Publication of CN116047592A publication Critical patent/CN116047592A/zh
Pending legal-status Critical Current

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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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/301—Analysis for determining seismic cross-sections or geostructures
    • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
    • 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

  • 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

一种三维地震沿层属性切片定量化色标的校正方法
技术领域
本发明涉及一种三维地震沿层属性切片定量化色标的校正方法，属于油气田开发技术领域。
背景技术
地震属性切片是用来进行层位和断层等地质信息解释，在油气开发过程中，通常需要借助地震属性切片信息明确储层空间展布特征及有利目标区分布。目前，在利用地震属性进行储层预测时，通常的做法是根据钻井与三维地震进行合成记录标定，明确储层所对应的地震时间范围(时窗)，形成地震属性平面分布图，再根据已完钻井与地震属性开展线性拟合，明确最优属性类型作为三维地震属性预测有效地质属性，根据已完钻井储层特征与地震属性的符合程度，人为选择属性切片色标。
上述地震数据分析过程中，主要是通过人为判定属性切片色标形式对地震属性进行定性评价，该评价过程主观因素强，使得地震属性切片不能准确反映中该范围内的储层特征，不能对属性起到定量评价的效果，储层预测精度难以达到最大化。
发明内容
本发明的目的是提供一种三维地震沿层属性切片定量化色标的校正方法，以解决现有三维地震属性切片色标不能准确反映和预测储层情况的问题。
本发明提供了一种三维地震沿层属性切片定量化色标的校正方法，该方法包括以下步骤：
1)获取目标区三维地震沿层属性切片数据及目标区内测井数据；
2)对获取的三维地震沿层属性切片数据和测井数据建立井震时深关系，确定测井分层顶底与地震切片的相对位置关系，得到测井分层顶底深度数据；
3)根据得到测井分层顶底深度数据确定该深度范围内的储层参数；
4)设置若干地震属性界限值对地震切片属性进行划分，设置若干储层参数界限值对储层参数进行划分；通过交叉计算每个地震属性界限值划分下各储层参数界限值对应的地震属性与储层参数的符合率，选取最优符合率对应的地震属性界限值和储层参数界限值；
5)根据最优符合率的界限值，调整三维地震沿层属性切片色标界限值，实现三维地震沿层属性切片色标校正。
本发明提供的三维地震沿层属性切片定量化色标的校正方法，通过设定相对应的若干三维地震切片属性界限值和若干测井分层储层参数界限值，采用交叉计算每个地震属性界限值划分下各个储层参数界限值对应的地震属性与储层参数的符合率，根据符合率最优的界限值调整三维地震沿层属性切片色标数值，使得校正后的三维地震沿层属性切片色标更能准确反映目标区储层情况。
进一步地，所述步骤1)中地震属性包括振幅类、频率类、相位类这三种类型。
进一步地，为了得到准确的井震时深关系，所述步骤2)中在建立井震时深关系前，需对测井数据进行标准化处理。
进一步地，所述步骤3)中储层参数为储层厚度和储层有效厚度。
进一步地，所述的储层厚度和储层有效厚度根据测井数据确定。
进一步地，为了提高计算效率和准确性，所述步骤4)中地震属性界限值是在对应地震属性值范围中设定区间内选取的，储层参数界限值是在对应储层参数范围中设定区间内选取的。
进一步地，为了提高计算效率和准确性，用于选取地震属性界限值的设定区间是对应地震属性值范围的1/3到2/3，用于选取储层参数界限值的设定区间是对应储层参数范围的1/3到2/3。
进一步地，为了提高交叉计算效率，所述步骤4)中交叉计算通过编程实现。
附图说明
图1是本发明三维地震沿层属性切片色标校正流程图；
图2是本发明实施例中的原始地震振幅属性切片图；
图3是本发明实施例中计算得到的最高符合率组合分布图；
图4是本发明实施例中色标校正后的地震振幅属性切片图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步地说明。
本发明提供了一种三维地震沿层属性切片定量化色标的校正方法，具体流程如图1所示。首先获取目标区的三维地震沿层属性切片数据和目标区测井数据；其次，建立两者之间的井震时深关系，确定测井分层顶底与地震切片的相对位置关系，得到测井分层顶底深度数据，并确定对应深度储层参数范围；之后，设置若干地震属性界限值对地震切片属性进行划分，设置若干储层参数界限值对储层参数进行划分；再通过交叉计算每个地震属性界限值划分下各储层参数界限值对应的地震属性与储层参数的符合率，选取最优符合率对应的地震属性界限值和储层参数界限值；最后，根据最优符合率的界限值，调整三维地震沿层属性切片色标界限值，实现三维地震沿层属性切片色标校正。
步骤1.获取目标区数据
本发明需要获取的数据包括目标区三维地震沿层属性切片数据和目标区测井数据。
其中，三维地震沿层属性切片指的是地震属性的切片数据，其中，地震属性类型包括振幅类(例如瞬时振幅、平均振幅、最大峰值振幅等)，频率类(例如瞬时频率、响应频率、瞬时频率斜率等)、相位类(例如瞬时相位、响应相位、瞬时相位余弦等)。在本实施例中，所选取的地震切片为T9d上10～20ms振幅属性切片，如图2所示。
其中，测井数据用于获取目标区储层参数。测井数据可以确定测井信息和地质信息之间的关系，其包含了各种参数，例如目标区测井储层物性(孔隙度、渗透率)、电性(声波时差、电阻率、密度、中子)、含气性(含气饱和度、全烃)等。在本实施例中，选取了目标区内96口完钻井的测井数据。
步骤2.建立井震时深关系
对步骤1中得到的地震属性切片数据和测井数据，需要确定切片与地层及分层顶底相对位置关系，保证测井分层与地震切片深度对齐。提取地震属性切片所使用的时窗顶底与对应地层的位置关系，通过时深关系计算出对应的分层位置，结合测井分层数据，得到测井分层顶底深度数据。
为保证得到准确的时深关系，在建立时深关系之前需先对测井数据进行标准化处理。由于环境影响或使用仪器不同会导致存在一定误差，需要通过对测井曲线进行标准化处理实现减小或消除误差，使处理后数据能够正确反映测井信息和地质信息之间的关系。其中可以采用频率直方图或趋势面分析对测井曲线进行标准化处理。在测井曲线完成标准化处理后，根据声波时差、密度这些测井曲线合成地震记录标定地层，确定测井分层顶底与地震切片的相对位置关系。
步骤3.确定储层参数范围
根据步骤2确定了测井分层顶底与地震切片的相对位置关系，得到测井分层顶底深度数据，找到测井对应深度的测井数据，并基于找到的测井数据确定对应的储层参数。其中，储层参数为储层厚度或储层有效厚度。在本实施例中，选取的储层参数为储层厚度，得到的储层厚度范围是0-20。
步骤4.计算地震属性和储层参数的符合率
为了得到准确的地震切片色标数值，需对地震属性和储层参数进行符合率计算。通过分别设定界限值对地震属性和储层参数进行划分，计算每个地震属性界限值划分下各储层参数界限值对应的地震属性与储层参数的符合率，得到最优符合率界限值对地震切片色标进行调整。
结合目标地质认识，对地震属性值和储层参数值设置若干界限值，通过设置若干地震属性界限值对地震切片属性进行划分，通过设置若干储层参数界限值对储层参数进行划分。其中，地震属性界限值在地震属性范围的1/3到2/3之间选取，储层参数界限值在储层参数范围的1/3到2/3之间选取。将设置若干界限值后的地震属性数据集分为好、差两类，其中好、差代表地质的好、差；同样将设置若干界限值后的储层参数数据集分成好、差两类，其中好、差代表物探的好、差；其中，根据门槛值判定好与差，一般将小于门槛值的定义为差，大于或等于门槛值的定义为好。
在本实施例中，地震属性值域为-2.8至-39，将地震属性值域的1/3到2/3(-15至-27)等差分成10份得到11个界限值，例如-15、-22、-27等；将地震属性依据每个界限值分成好、差两类，共能得到10组地震属性分类数据，以界限值-22为例，将大于或等于-22的地震属性定义为地质好，将小于-22的地震属性定义为地质差，可得到一组地震属性分类数据。在本实施例中，储层参数值域为0至20，将储层参数值域的1/3到2/3(7至14)等差分成10份得到11个界限值，例如7、10.5、14等；将储层参数依据每个界限值分成好、差两类，共能得到10组储层参数分类数据，以界限值10.5为例，将大于或等于10.5的储层参数定义为物探好，将小于10.5的储层参数定义为物探差，便得到一组储层参数分类数据。对得到的10组地震分类数据和储层参数分类数据进行交叉计算，分别计算每一组地震分类数据与10组储层参数分类数据之间的符合率，得到符合率最高的地震属性界限值和储层参数界限值；以地震属性界限值-22得到的地质好、差数据和储层参数界限值10.5为例得到的物探好、差数据为例，当地质好、物探好时或地质差、物探差时，判定两组数据符合，当地质差、物探好时或地质好、物探差时，判定两组数据不符合；以此为例，通过编程交叉计算10组地震属性分类数据与10组储层参数分类数据的符合率，可以得到100个符合率计算结果，其中最优符合率为89.8％，因此可以得到最优符合率下的组合分布图(如图3所示)，同时可以得到最优符合率下的地震属性界限值为-18和储层参数界限值为7。
步骤4.地震切片色标校正
根据步骤3中得到的最优符合率下对应的地震属性界限值和储层参数界限值，调整对应的地震切片色标数值，实现地震沿层属性切片校正。在本实施例中，根据得到的符合率最高的地震属性界限值为-18，调整地震属性切片色标数值，得到最优地震属性切片图(如图4所示)。

Claims (8)

1.一种三维地震沿层属性切片定量化色标的校正方法，其特征在于，该方法包括以下步骤：
1)获取目标区三维地震沿层属性切片数据及目标区内测井数据；
2)对获取的三维地震沿层属性切片数据和测井数据建立井震时深关系，确定测井分层顶底与地震切片的相对位置关系，得到测井分层顶底深度数据；
3)根据得到测井分层顶底深度数据确定该深度范围内的储层参数；
4)设置若干地震属性界限值对地震切片属性进行划分，设置若干储层参数界限值对储层参数进行划分；通过交叉计算每个地震属性界限值划分下各储层参数界限值对应的地震属性与储层参数的符合率，选取最优符合率对应的地震属性界限值和储层参数界限值；
5)根据最优符合率的界限值，调整三维地震沿层属性切片色标界限值，实现三维地震沿层属性切片色标校正。
2.根据权利要求1所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述步骤1)中地震属性包括振幅类、频率类、相位类这三种类型。
3.根据权利要求1所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述步骤2)中在建立井震时深关系前，需对测井数据进行标准化处理。
4.根据权利要求1所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述步骤3)中储层参数为储层厚度和储层有效厚度。
5.根据权利要求4所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述的储层厚度和储层有效厚度根据测井数据确定。
6.根据权利要求1所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述步骤4)中地震属性界限值是在对应地震属性值范围中设定区间内选取的，储层参数界限值是在对应储层参数范围中设定区间内选取的。
7.根据权利要求6所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，用于选取地震属性界限值的设定区间是对应地震属性值范围的1/3到2/3，用于选取储层参数界限值的设定区间是对应储层参数范围的1/3到2/3。
8.根据权利要求1所述的三维地震沿层属性切片定量化色标的校正方法，其特征在于，所述步骤4)中交叉计算通过编程实现。
CN202111265039.2A 2021-10-28 2021-10-28 一种三维地震沿层属性切片定量化色标的校正方法 Pending CN116047592A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111265039.2A CN116047592A (zh) 2021-10-28 2021-10-28 一种三维地震沿层属性切片定量化色标的校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111265039.2A CN116047592A (zh) 2021-10-28 2021-10-28 一种三维地震沿层属性切片定量化色标的校正方法

Publications (1)

Publication Number Publication Date
CN116047592A true CN116047592A (zh) 2023-05-02

Family

ID=86131849

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111265039.2A Pending CN116047592A (zh) 2021-10-28 2021-10-28 一种三维地震沿层属性切片定量化色标的校正方法

Country Status (1)

Country Link
CN (1) CN116047592A (zh)

Similar Documents

Publication Publication Date Title
CN104047597B (zh) 含油气泥页岩地层测井曲线标准化方法
US20130325349A1 (en) Methods for Generating Depofacies Classifications for Subsurface Oil or Gas Reservoirs or Fields
CA2213125C (en) Determination of formation properties from seismic attributes
CN111596364B (zh) 基于高精度层序地层格架的地震沉积微相组合分析方法
AU2009303613B2 (en) Pseudo logs to improve stratigraphic correlation between sedimentary basins
US9869792B2 (en) Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells
WO2006135567A2 (en) Method of processing seismic data to extract and portray avo information
CN114609675A (zh) 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法
CN112946782B (zh) 一种致密油气储渗体地震精细刻画方法
CN116047592A (zh) 一种三维地震沿层属性切片定量化色标的校正方法
CN113484907B (zh) 一种预测不同类型储层平面上分布的方法
Gastaldi et al. Using Bayesian simulations to predict reservoir thickness under tuning conditions
CN114895373A (zh) 一种基于高频率范围约束的测井曲线标准化方法及系统
CN112859172A (zh) 纵横波微测井资料处理方法及装置
CN113325474B (zh) 生物礁判别方法
CN112287532B (zh) 一种缘控地震驱动建模方法
CN116338775A (zh) 一种薄储层识别方法及装置
Oladapo et al. Correlation and Calibration of 4D Seismic Vintages for Effective Postproduction Intervention
CN117930380A (zh) 基于地震沉积学的碳酸盐岩储层预测方法和系统
CN116520407A (zh) 利用地震多属性基于两次随机森林法预测薄砂体厚度的方法
AU709540B2 (en) Determination of formation properties from seismic attributes
CN116843034A (zh) 基于井数据预处理的机器学习样本制作方法
CN116027417A (zh) 一种基于差异共振的储层预测方法
Cameron Analysis of dipmeter data for sedimentary orientation
CN116027429A (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