CN113759419A - 一种储层预测方法、装置、存储介质及电子设备 - Google Patents

一种储层预测方法、装置、存储介质及电子设备 Download PDF

Info

Publication number
CN113759419A
CN113759419A CN202010500243.7A CN202010500243A CN113759419A CN 113759419 A CN113759419 A CN 113759419A CN 202010500243 A CN202010500243 A CN 202010500243A CN 113759419 A CN113759419 A CN 113759419A
Authority
CN
China
Prior art keywords
wave impedance
inversion
reservoir
velocity
layer
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
CN202010500243.7A
Other languages
English (en)
Other versions
CN113759419B (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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN202010500243.7A priority Critical patent/CN113759419B/zh
Publication of CN113759419A publication Critical patent/CN113759419A/zh
Application granted granted Critical
Publication of CN113759419B publication Critical patent/CN113759419B/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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • 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/303—Analysis for determining velocity profiles or travel times
    • 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/303—Analysis for determining velocity profiles or travel times
    • G01V1/305—Travel times
    • 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/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • 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/622—Velocity, density or impedance
    • G01V2210/6226—Impedance
    • 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
    • G01V2210/6244—Porosity
    • 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
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

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

一种储层预测方法、装置、存储介质及电子设备
技术领域
本申请涉及石油勘探技术领域，特别地涉及一种储层预测方法、装置、存储介质及电子设备。
背景技术
储层预测是综合应用地震、地质、钻井、测井等各项资料对地下储集层的分布、厚度及岩性和物理性质变化进行追踪和预测的一项先进技术，随着勘探过程越来越复杂，提高储层预测精度也愈发重要。目前，国内外发展了多种储层预测方法，如地震相分析、地震反演、地震模拟分析、地震属性分析等，其中地震反演与模拟是油藏描述中应用最广泛的勘探技术。
目前，在复杂礁滩储层预测中，大多采用通过沉积相变线控制井间插值与外推来进行相控储层的预测，使储层定量预测精度提高，但其沉积相变线通过人工定义存在误差；同时，沉积相变线控制井间插值与外推的相控方法仅为一个平面上控制。但是，由于生物礁气藏是一种特殊地质体，其沉积速度高于周围地层且发育具有随机性，因此，生物礁气藏具有空间形态复杂，横向变化速度快的特征，无法凭借解释层位生成合理的初始模型；且同时生物礁气藏在纵向上存在规模大小的差异，因此用上部沉积相平面图控制下部的内插外推缺乏合理性。除此之外，由于生物礁气藏储层与非储层在阻抗上存在无法完全区分的现象，因此仅凭阻抗反演无法精确刻画出生物礁气藏储层的空间展布。
综上，以一个平面的趋势控制一个体的储层预测会存在误差，且仅凭阻抗反演无法精确刻画生物礁气藏储层的空间展布，导致在利用上述方法对生物礁气藏储层进行精细刻画时存在不准确的问题。
因此，如何准确对生物礁气藏储层进行储层预测是目前亟待解决的问题。
发明内容
本申请提供一种储层预测方法、装置、存储介质及电子设备，解决了现有技术中无法准确对生物礁气藏储层进行储层预测的问题。
第一方面，本申请提供了一种储层预测方法，所述方法包括：
获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的传播速度；
采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；
根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体；
获取电阻率测井曲线；
根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；
根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
根据本申请的实施例，可选的，在上述储层预测方法中，所述获取待测井点的层速度体，包括：
采用迪克斯计算式得到所述地层中每层的层速度，其中，所述迪克斯计算式为：
Figure BDA0002524436050000021
其中，νn为所述地层中第n层的层速度，t1，n为所述地层中第一层到第n层的地震双程垂直反射时间，νR，n为所述地层中第一层到第n层的均方根速度，t1，n-1为所述地层中第一层到第n-1层的地震双程垂直反射时间，νR，n-1为所述地层中第一层到第n-1层的均方根速度。
根据本申请的实施例，可选的，在上述储层预测方法中，所述层速度-阻抗关系式为：
Imp＝a*f(V)+b；
其中，Imp为波阻抗体，f(V)为层速度体，a为第一预设参数，b为第二预设参数。
根据本申请的实施例，可选的，在上述储层预测方法中，根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体，包括：
根据所述波阻抗体和所述测井数据，建立低频模型；
对所述低频模型进行波阻抗反演，以得到波阻抗反演体。
根据本申请的实施例，可选的，在上述储层预测方法中，根据所述电阻率模拟数据体对所述波阻抗反演体进行处理，包括：
根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值；
将所述波阻抗反演体中位于所述储层门槛值以外的阻抗值剔除。
根据本申请的实施例，可选的，在上述储层预测方法中，根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值，包括：
将所述电阻率模拟数据体和所述波阻抗反演体进行交会分析处理，以得到储层门槛值。
根据本申请的实施例，可选的，在上述储层预测方法中，将处理后的波阻抗反演体转换为孔隙度反演体，包括：
建立波阻抗反演体和孔隙度的线性关系；
根据建立的所述线性关系，将处理后的波阻抗反演体转换为孔隙度反演体。
第二方面，本申请提供了一种储层预测装置，所述装置包括：
第一获取模块，用于获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的层速度；
第一转换模块，用于采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；
反演模块，用于根据所述波阻抗体和测井数据，进行波阻抗反演，得到波阻抗反演体；
第二获取模块，用于获取电阻率测井曲线；
波形指示模拟模块，用于根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；
第二转换模块，用于根据所述电阻率模拟数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
第三方面，本申请提供了一种存储介质，该存储介质存储有计算机程序，当该计算机程序被一个或多个处理器执行，实现如上述第一方面中任意一种储层预测方法。
第四方面，本申请提供了一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，实现如上述第一方面中任意一种储层预测方法。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
本申请提供的一种储层预测方法、装置、存储介质及电子设备，所述方法包括：获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的传播速度；采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；根据所述波阻抗体和测井数据，进行波阻抗反演，得到波阻抗反演体；获取电阻率测井曲线；根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。解决了现有技术中无法准确对生物礁气藏储层进行储层预测的问题，避免了现有技术中以一个平面的数据趋势控制一个体的建模存在的误差，从而导致储层刻画不准确的问题，实现了对生物礁气藏储层的精细刻画。
附图说明
在下文中将基于实施例并参考附图来对本申请进行更详细的描述。
图1为本申请实施例一提供的一种储层预测方法的流程示意图。
图2为本申请实施例一提供的一种储层预测方法的另一流程示意图。
图3为本申请实施例二提供的一种层速度连井剖面示意图。
图4为本申请实施例二提供的一种常规波阻抗反演连井剖面示意图。
图5为本申请实施例二提供的一种电阻率模拟数据体连井剖面示意图。
图6为本申请实施例三提供的一种储层预测装置的连接框图。
在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
以下将结合附图及实施例来详细说明本申请的实施方式，借此对本申请如何应用技术手段来解决技术问题，并达到相应技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本申请的保护范围之内。
实施例一
请参阅图1，本实施例提供一种可应用于手机、电脑或平板电脑等电子设备的储层预测方法，该方法主要针对生物礁气藏储层的预测。所述储层预测方法应用于所述电子设备时执行步骤S110至步骤S160。
步骤S110：获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的传播速度
在本实施例中，层速度是指在层状地层中地震波传播的速度，所述层速度体即包括地震波在每层地层的传播速度。所述层速度直接反映地层的岩性，能用来划分地层。
在一些实施例中，所述层速度体还可以理解为一种关于地层层数与层速度两者之间的关系曲线。其中，计算所述层速度具体可以包括以下四种方式：
(1)利用声波测井资料求取层速度。采用这种层速度计算方法求得的层速度准确度和精度都比较高，便于分析速度与岩性之间的对应关系。但是，由于井径和钻井液侵入带对层速度具有不同程度上的影响，而且由于井数量(即待测井点数量)的限制，导致这种声波测井资料十分有限，此外，所述声波测井资料也分布不均，无法进行区域性研究，实用性较低。
(2)利用地震测井资料求取层速度，所述地震测井资料包括VSP(VerticalSeismic Profiling，垂直地震剖面)资料。该种层速度计算方法求取得到的层速度的精度较低，因此，只能只反映一些大段地层的速度差别，但其求得的平均速度精度较高。然而，与所述声波测井资料的缺陷相似，该地震测井资料也十分有限，无法进行区域性研究，实用性较低。
(3)利用叠加速度谱资料求取层速度。速度谱(velocity spectrum)与频谱的概念类似，把地震波的能量相对于波速的变化关系的曲线称为速度谱。虽然叠加速度谱资料精度较低，但由于叠加速度谱资料数量巨大，可控制全盆地，因而，利用这种层速度计算方法可用于研究目的层速度的区域性变化规律，与其他层速度计算方法相比，其实用性得到显著提高。
(4)利用合成速度测井求层速度。这种层速度计算方法求取得到的层速度精度较高，并可以在研究区进行连续速度分析，不仅能用于地震层序层速度计算，还能用于单砂体层速度求取，但是，合成速度测井资料成本的很高，因此，导致研究成本过高。
在本实施例中，以上述第(3)种层速度计算方法得到层速度为例，对本实施例进一步进行说明。具体的：
由于叠加速度与均方根速度存在关联关系，上述第(3)种层速度计算方法也可以理解为根据均方根速度求取层速度。
采用迪克斯计算式得到所述地层中每层的层速度，其中，所述迪克斯计算式为：
Figure BDA0002524436050000061
在上式中，νn为所述地层中第n层的层速度，t1，n为所述地层中第一层到第n层的地震双程垂直反射时间，νR，n为所述地层中第一层到第n层的均方根速度，t1，n-1为所述地层中第一层到第n-1层的地震双程垂直反射时间，νR，n-1为所述地层中第一层到第n-1层的均方根速度，其中，n为大于2的正整数。
示例性的，假设地层中第一层和第一层到第二层的均方根速度，以及地层中第一层和第一层到第二层的地震双程垂直反射时间为已知数，并将这些已知数代入上式中，即可求得第一层的层速度的平方，再通过开方处理，即可求得第一层对应的层速度。
由于叠加速度和均方根速度存在关联关系，因此，可以通过每层地层的叠加速度求取与每层对应的均方根速度，在一些实施例中，叠加速度等于均方根速度，例如，在水平层状介质的情况下(或水平界面覆盖为连续介质)；在一些实施例中，叠加速度不等于均方根速度，例如，在倾斜层状介质的情况下，叠加速度就等于均方根速度乘以倾角的余弦。
在一些实施例中，可以根据共反射点时距曲线求得叠加速度。地震波从震源出发，经过地下介质传播到测线上各观测点的旅行时间，同观测点相对于激发点的距离之间的关系。
需要说明的是，本领域技术人员应当理解，在求取叠加速度时还可以选择本领域其他计算方法求取叠加速度，对此，本实施例不做任何限定。在实际应用中，根据实际需求设定即可。
在本实施例中，先根据共反射点时距曲线求得每层地层对应的叠加速度，进而根据地质的情况(水平层状介质或倾斜层状介质)确定叠加速度与均方根速度的关联关系，再根据叠加速度与均方根速度的关联关系求得每层地层对应的均方根速度，再借助上述迪克斯计算式得到每层地层对应的层速度。
需要说明的是，本领域技术人员应理解，上述层速度计算方法并不限于利用均方根速度求取层速度的方法，其并不造成对本申请的限定。本申请还可以根据实际需求，选取其他层速度计算方法求取层速度。
在本实施例中，所述测井数据一般都是各种不同的物理参数，如电阻率、自然电位、声波速度、岩石体积密度等，可统称为测井信息。而测井数据解释与数字处理的成果，如岩性、泥质含量、含水饱和度、渗透率、地层信息等，可统称为地质信息。
步骤S120：采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体。
在本实施例中，首先，需要建立层速度体与波阻抗体的关联关系，根据该建立的关联关系，即可将层速度体转换为波阻抗体。
可以理解的是，在所述层速度体表征为地层层数与层速度两者之间的关系曲线时，所求取的波阻抗体也可表征为地层层数与波阻抗值的关系曲线。
在本实施例中，所述层速度-阻抗关系式具体为：
Imp＝a*f(V)+b；
在上式中，Imp为波阻抗体，f(V)为层速度体，a为第一预设参数，b为第二预设参数。
需要说明的是，所述第一预设参数和所述第二预设参数根据实际情况设置即可。
在现阶段的地震勘探领域，地震反演已经从幕后走到台前，而且所扮演的角色越来越重要。其以地质为主导，测井为约束，地震为主体，来更精确的了解认识地下岩层结构的宏观空间分布规律，微观岩石物理性质。
由于人类对地球认识，寻找油气等各种生活生产的需要，所以有必要研究地球内部情况，但是以现在的科学和技术水平，只能从地表得到关于地球的信息，只有很小一部分地球内部空间观测数据是来自于钻井。所以只能通过有限的点上的地下观测数据，结合地表地震观测信息，通过反演来获得地球面上或体上的地球深部地质信息。
反演是以各种场反馈的抽象数据为基础，推测地球内部空间特征、定性或定量分析物理特征的过程。地震反演的过程，简言之，是由地表观测的地震反馈数据(例如，层速度)，结合地质情况和测井资料，依据数学物理原理方法，对地球内部结构和物理属性参数进行再认识的过程，也就是说是获得地下岩性分布特征和其物理参数等信息的过程。反演数据比地震振幅数据分辨率更高而且更直观，能够直接从反演剖面上分辨出储层特征，从而地震反演技术能够通过帮助剖析钻井结果的不确定性和评估生产风险来减少油气开发的损失。因此反演技术是目前人类了解认知地下空间，开发利用地球资源的重要手段。
通过地震反演技术，将地震数据和测井数据进行合理的匹配非常重要，搭建地震时间域和测井深度域进行直接对比的桥梁，以此为依据进行地质构造和岩石物性解释，研究储层空间的特征变化规律，为油气的勘探开发提供重要依据。其中，地震波阻抗反演是地震勘探中最常见的方法，地震波阻抗反演是储层预测的有效手段之一，波阻抗与含油气储层有很好的对应性，为油气资源的勘探和开发提供有效的支持。
波阻抗反演实际上是从地震剖面上消除子波影响，留下反射系数，再由反射系数计算出能反映地层物性变化的物理参数波阻抗。
通过将层速度体转换为波阻抗体，并将该波阻抗体运用到反演中，以一个体的数据趋势反映整个生物礁气藏储层的空间展布，有利于精确刻画整个生物礁气藏的储层。
步骤S130：根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体。
在本实施例中，具体利用波阻抗反演方法得到波阻抗反演体。
从广义上讲，所述波阻抗反演方法主要包括叠前波阻抗反演和叠后波阻抗反演。所述叠前波阻抗反演是根据叠前地震资料进行反演；所述叠后波阻抗反演是根据叠后地震资料进行反演。所述叠前波阻抗反演和所述叠后波阻抗反演具有不同的优点。具体的：
所述叠前地震资料中包含更多的频率信息，叠后地震资料相当于经过了滤波处理(信噪比相对较高)，单从频率上看，利用叠前地震资料进行叠前反演得到的反演结果会具有更好的分辨率，但是由于分辨率仅仅只是地震解释的一个方面，而叠前地震资料的信噪比较低，且如果叠前地震资料的信噪比过低，即使获取的反演剖面的分辨率较高，也有可能是虚假的信息，进而导致根据该反演剖面进行地震解释存在一定的误差。因此，由于上述两种反演方法具有不同的优点和缺点，在实际应用中，可根据实际需求，选择对反演结果影响最小的波阻抗反演方法。
在一些实施例中，所述叠前波阻抗反演主要包括基于波动方程的反演和弹性参数波阻抗反演。
在本实施例中，所述叠后波阻抗反演主要包括线性反演和非线性反演。
在本实施例中，所述线性反演主要包括道积分反演、波阻抗递推反演、稀疏脉冲波阻抗反演、宽带约束波阻抗反演和广义线性波阻抗反演。
在本实施例中，所述非线性反演包括模拟退火波阻抗反演、随机模拟反演、随机反演和遗传算法波反演。
需要说明的是，本领域技术人员应当理解，本实施例中，根据所述波阻抗体和所述测井数据进行波阻抗反演的波阻抗反演方法可采用叠后波阻抗反演，也可以采用叠前波阻抗反演。对此，本实施例并不做任何限定，在实际应用中，根据实际需求设定即可。
在本实施例中，可采用弹性参数波阻抗反演方法得到波阻抗反演体。具体的：
弹性参数波阻抗反演是地震数据和测井数据视为互为约束的数据，利用变差函数在测井尺度统计不同沉积微相的空间变化规律，在地震尺度通过协克里金将其与弹性参数进行融合，通过线性加权在进行层序格架的控制下实现在远离井的地方融合地震趋势数据来校正插值数据，其中，这里的地震趋势数据就是本实施例中生成的波阻抗体。
在本实施例中，所述波阻抗体和所述测井数据的约束关系可以为：
Figure BDA0002524436050000091
在上式中，X*(u0)为波阻抗反演体，u0为空间位置点，
Figure BDA0002524436050000092
基于变差分析的加权系数值，
Figure BDA0002524436050000093
为
Figure BDA0002524436050000094
处的测井数据，
Figure BDA0002524436050000095
为测井数据的空间位置点，
Figure BDA0002524436050000096
为基于变差分析的加权系数值，
Figure BDA0002524436050000097
为
Figure BDA0002524436050000098
处的波阻抗值，
Figure BDA0002524436050000099
为波阻抗值的空间位置点，m为正整数，i为正整数，k为正整数。
请参阅图2，所述步骤S130具体可以包括步骤1301和步骤1302：
S1301：根据所述波阻抗体和所述测井数据，建立低频模型。
S1302：对所述低频模型进行波阻抗反演，以得到波阻抗反演体。
在本实施例中，通过对待测井点的波阻抗体进行内插外推，建立波阻抗体的低频模型，再对所述低频模型进行波阻抗反演，以得到波阻抗反演体。
步骤S140：获取电阻率测井曲线。
在地震勘探技术中，测井数据中的自然电位、自然伽马、电阻率以及种子孔隙度等非波阻抗类曲线与地震反射虽没有直接对应关系，但可直接反应地层的岩性和含油性，因此，通过该类曲线和原始波阻抗反演结果进行综合解释，即可对储层进行更为准确的刻画。因此，选取可直接反映生物礁气藏储层的曲线变得尤为重要。
在地震勘探中，测井曲线主要包括电阻率曲线、声波曲线、自然电位曲线、微电极曲线、中子测井曲线、三侧向测井曲线以及微球形聚焦测井曲线。每种测井曲线具有不同的特点。具体的：
(1)电阻率曲线，双侧向是探测不同径向深度电阻率的测井方法。通常情况下，裂缝的存在使双侧向出现差异，模拟实验表明，低角度裂缝的双侧向值呈负差异，而高角度裂缝的双侧向值呈正差异，双侧向幅度差不仅与裂缝的产状有关，而且与裂缝的张开度有关，因此在一些裂缝段也可能无差异。
(2)声波曲线，裂缝在声波曲线上的反映与井筒周围裂缝的产状及发育程度有关。声波曲线对高角度裂缝没有反映，对低角度裂缝或网状裂缝，声波测井值将相应增大；当遇到大的水平裂缝或网状裂缝时，声波能量急剧衰减而产生“周波跳跃”现象。因此利用声波时差可以识别水平裂缝或网状裂缝，但不能用于识别垂直裂缝。声波曲线对裂缝的显示主要取决于裂缝的张开度、发育程度、充填物和流体的性质。
(3)自然电位曲线，其具有的特点包括：a.当地层、泥浆是均匀的，上下围岩岩性相同，自然电位曲线对地层中心对称；b.在地层顶底界面处，自然电位变化最大，当地层较厚(大于四倍井径)时，可用曲线半幅点确定地层界面；c.测量的自然电位幅度为自然电流在井内产生的电压降，它永远小于自然电流回路总的电动势；d.渗透性砂岩的自然电位，对泥岩基线而言，当地层水矿化度大于泥浆滤液矿化度时，自然电位显示为负异常，当地层水矿化度小于泥浆滤液矿化度时，显示为正异常，如果泥浆滤液的矿化度与地层水矿化度大致相等时，自然电位偏转幅度很小，曲线无显示异常。
(4)微电极曲线，其具有的特点是在渗透性地层有幅度差，微电位值大于微梯度值。
(5)感应测井曲线，其具有的特点包括：a.上下围岩相同，单一低电导率地层,当地层厚度大于1.7米时，曲线上可以看到过聚焦产生的局部极值，其厚度小于1.7米时，视电导率曲线呈现一尖峰；b.上下围岩不同，单一低电导率地层,对于厚度大于2米的地层，地层中部的曲线呈倾斜状，地层中心对应于倾斜段的中点，对于厚度小于2米的地层，视电导率曲线偏向与地层电导率差别小的围岩一侧，这是在高低电导率地层，而在中间电导率地层的曲线，对于厚度大于2米的地层，呈比较清楚的台阶状。
(6)中子测井曲线，其具有的特点包括：a.在砂泥岩剖面中，粘土(泥岩)的中子测井计数率最低，致密砂岩的中子测井计数率最高，粉砂岩、泥质砂岩、孔隙中充满液体的砂岩为中等数值；b.气层的中子测井计数率是高值。
(7)三侧向测井曲线，其具有的特点包括：a.高阻层视电阻率曲线对围岩形成高阻异常，异常对称于高阻层中点，异常极大值为视电阻率代表值。如果地层较厚，岩性、电性不均匀，分段取值；b.高阻层界面在三侧向曲线上缺乏明显的特征，但靠近高阻异常的底部；c.深浅三侧向曲线的形态相同，在储集层有幅度差。
(8)微球形聚焦测井曲线，其主要反映冲洗带电阻率，受泥饼和原状地层影响。
在一些实施例中，可以通过对各种测井曲线进行交会分析处理，以得到对所述生物礁气藏的储层最为敏感(最能反映所述生物礁气藏储层的岩性)的测井曲线。
需要说明的是，经过研究，最能反映所述生物礁气藏储层岩性的测井曲线是电阻率测井曲线。
可以理解的是，为了提高储层预测的结果，可对获取到的测井曲线进行预处理。测井曲线预处理是测井解释与数据处理的一项基础工作，它是保证测井解释与数据处理结果精度的重要前提。
在本实施例中，所述预处理包括但不限于深度校正、平滑滤波、环境影响校正等。
在一些实施例中，所述深度校正可采用深度控制曲线进行深度校正，也可以利用相关函数进行深度校正。
在一些实施例中，采用小波变换分析方法进行平滑滤波处理。
步骤S150：根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体。
在地震勘探技术领域中，通过对波阻抗反演数据体和测井数据(例如，电阻率测井曲线)的分析，建立起波阻抗数据与储层参数的相关关系。在此基础上，主要研究了由反演波阻抗数据体，通过一定方式转换，获得所要预测的储层参数数据体，从而为各种油藏描述提高丰富的储层信息，以达到对储层精确预测的目的。
在本实施例中，根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟以得到电阻率模拟数据体具体包括：
建立一个输入数据(波阻抗反演体和电阻率测井曲线)和输出数据(电阻率模拟数据体)的非线性关系，根据该非线性关系，即可获得与所述波阻抗反演体和所述电阻率测井曲线对应的电阻率模拟数据体。
步骤S160：根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
在本实施例中，根据所述电阻率数据体对所述波阻抗反演体进行处理具体包括：
首先，根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值。然后，将所述波阻抗反演体中位于所述储层门槛值以外的阻抗值剔除。
在本实施例中，所述储层门槛值有效地划分了储层和非储层，通过该储层门槛值，可以有效将反演结果中的储层和非储层进行划分，达到精确刻画生物礁气藏储层的空间展布。
在本实施例中，将所述电阻率模拟数据体和所述波阻抗反演体进行交会分析处理，以得到储层门槛值。
需要说明的是，本领域技术人员应当理解，所述储层门槛值还可以通过其他方式获得，对此，本实施例并不做任何限定。在实际应用中，根据实际需求设定即可。
在地震勘探技术领域中，岩石储层物性反演是近年来发展起来的一种反演方法，其目的是更为直接地把地震数据与地质认识结合起来，可以更为直观地为地质人员提供储层解释依据及油气判别依据。
常用的物性反演包括电阻率反演、伽马反演、孔隙度反演、渗透率、饱和度以及其他岩性参数反演。其中，孔隙度反演的主要作用是进行储层量化预测，适用于评价及开发阶段。因此，将处理后的波阻抗反演体转换为孔隙度反演体可更为直观的反映储层信息，便于对储层进行评价。
在本实施例中，将处理后的波阻抗反演体转换为孔隙度反演体具体包括：
首先，建立处理后的波阻抗反演体和孔隙度的线性关系；然后，再根据建立的所述线性关系，将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
在本实施例中，解决了现有技术中无法准确对生物礁气藏储层进行储层预测的问题，整个过程减少了相变线的人工定义，避免了人工定义引起的误差。充分利用了地震资料中的层速度体，进行了数据体层面的储层预测，避免了现有技术中以一个平面的数据趋势控制一个体的建模存在的误差，从而导致储层刻画不准确的问题，实现了对生物礁气藏储层的精细刻画，并提高了生物礁气藏开发经济效益。
实施例二
为了便于理解本申请的方案与效果，本实施例提供一种具体应用实例，以对本方案进行进一步说明。需要说明的是，本领域技术人员应当理解，该示例仅为了便于理解本申请，其任何具体细节并非意在以任何方式限定本申请。具体请参阅图3至图5，具体的：
图3是一种生物礁储层的层速度图剖面示意图，在生物礁发育位置速度有着明显的边界，具有相控的特征。
图4一种常规波阻抗反演连井剖面示意图，待测井点1、待测井点2和待测井点3的储层位置(对应图4中浅色部分)与低阻抗并无对应关系，低阻抗无法反映生物礁储层的发育位置。
图5是一种电阻率模拟数据体连井剖面示意图，该图中低波阻抗值区域与储层能够完全对应，很好的刻画了生物礁储层的空间展布特征。因此，利用电阻率模拟数据体，对波阻抗反演体进行处理，能准确的刻画生物礁储层的空间展布特征。
井名 待测井点1 待测井点2 待测井点3
测试产量(万方/天) 58.97 33.6 23.56
如上表所示，上表为某工区内三口待测井点(包括待测井点1、待测井点2和待测井点3)的测试产能情况。由上表可知，所述待测井点1的测试产能最高，在电阻率数据模拟体上其生物礁储层规模最大，物性最好。而所述待测井点2与所述待测井点3井的测试产能较低，电阻率数据模拟体上其生物礁储层规模较薄较小。
结合通过电阻率数据模拟体和波阻抗反演体对生物礁储层进行预测，并结合各个待测井点的实际测试产能与预测结果进行对比，验证了实施例一中储层预测方法对生物礁气藏储层预测的可靠性。
因此，通过上述应用实例，得到测试产量与刻画出的生物礁储层规模以及储层物性有着较好的对应关系，验证了储层预测方法的可靠性，并提高了生物礁气藏开发效率。
实施例三
本实施例提供一种储层预测装置，请参阅图6，所述装置包括：第一获取模块201、第一转换模块202、反演模块203、第二获取模块204、波形指示模拟模块205以及第二转换模块206。
第一获取模块201，用于获取待测井点的层速度体和测井数据。
在一些实施例中，所述第一获取模块201包括层速度计算模块，所述层速度计算模块用于采用迪克斯计算式得到所述地层中每层的层速度。其中，所述迪克斯计算式为：
Figure BDA0002524436050000141
在上式中，νn为所述地层中第n层的层速度，t1，n为所述地层中第一层到第n层的地震双程垂直反射时间，νR，n为所述地层中第一层到第n层的均方根速度，t1，n-1为所述地层中第一层到n-1层的地震双程垂直反射时间，νR，n-1为所述地层中第一层到第n-1层的均方根速度。
第一转换模块202，用于采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体。
在一些实施例中，所述第一转换模块202具体包括第一建立模块和波阻抗体计算模块。具体的：
在本实施例中，所述第一建立模块具体用于建立层速度体与波阻抗体的层速度-阻抗关系式。所述波阻抗体计算模块具体用于根据第一建立模块建立的层速度-阻抗关系式，将所述层速度体转换为波阻抗体。
所述层速度-阻抗关系式具体为：
Imp＝a*f(V)+b；
在上式中，Imp为波阻抗体，f(V)为层速度体，a为第一预设参数，b为第二预设参数。
需要说明的是，所述第一预设参数和所述第二预设参数根据实际情况设置即可。
反演模块203，用于根据所述波阻抗体，进行波阻抗反演，得到波阻抗反演体。
在一些实施例中，所述反演模块203具体包括第二建立模块和波阻抗反演模块。具体的：
在本实施例中，所述第二建立模块用于根据所述波阻抗体和所述测井数据，建立低频模型。所述波阻抗反演模块用于对第二建立模块所建立的低频模型进行波阻抗反演，以得到波阻抗反演体。
第二获取模块204，用于获取电阻率测井曲线。
在一些实施例中，第二获取模块204具体包括测井曲线获取模块和交会分析模块。具体的：
在本实施例中，所述测井曲线获取模块用于获取多种测井曲线，所述交会分析模块用于对多种测井曲线进行交会分析，以得到反映生物礁气藏储层岩性的测井曲线。
需要说明的是，最能反映生物礁气藏储层岩性的测井曲线是电阻率测井曲线。
可以理解的是，为了提高储层预测的结果，可对获取到的测井曲线进行预处理。测井曲线预处理是测井解释与数据处理的一项基础工作，它是保证测井解释与数据处理结果精度的重要前提。因此，所述第二获取模块具体包括预处理模块，用于对获取多种测井曲线进行预处理。
在本实施例中，所述预处理包括但不限于深度校正、平滑滤波、环境影响校正等。
可以理解的是，在所述预处理包括深度校正和平滑滤波时，所述预处理模块具体包括深度校正模块和平滑滤波模块。具体的：
在本实施例中，所述深度校正模块用于采用深度控制曲线对测井曲线进行对深度校正。
在本实施例中，所述平滑滤波模块用于小波变换分析方法对测井曲线进行对滤波处理。
需要说明的是，本实施例所述预处理包括但不限于深度校正和平滑滤波处理。对此，本实施例并不做任何限定，在实际应用中，可根据实际需求设定其他预处理即可。
波形指示模拟模块205，用于根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体。
在本实施例中，所述波形指示模块205具体包括第三建立模块和模拟模块。
在本实施例中，所述第三建立模块，用于根据波阻抗反演数据和电阻率测井曲线，建立波阻抗-电阻率对应关系。所述模拟模块用于根据所述波阻抗反演体和波阻抗-电阻率对应关系，得到电阻率模拟数据体。
在地震勘探技术领域中，通过对波阻抗反演数据和测井数据(例如，电阻率测井曲线)的分析，建立起波阻抗数据与储层参数的相关关系，在此基础上，主要研究了由反演波阻抗数据体，通过一定方式转换，获得所要预测的储层参数数据体，从而为油藏描述提高丰富的储层信息，达到储层预测的目的。
第二转换模块206，用于根据所述电阻率模拟数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
在本实施例中，所述第二转换模块206具体包括储层门槛值确定模块、剔除模块和孔隙度反演模块。
在本实施例中，所述储层门槛值确定模块，具体用于根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值。在本实施例中，所述储层门槛值确定模块具体用于将所述电阻率模拟数据体和所述波阻抗反演体进行交会分析处理，以得到储层门槛值。
在本实施例中，所述储层门槛值有效地划分了储层和非储层，通过该储层门槛值，可以有效将反演结果中的储层和非储层进行划分，达到精确刻画生物礁气藏储层的空间展布。
需要说明的是，所述储层门槛值还可以通过其他方式获得，对此，本实施例并不做任何限定。
所述剔除模块，具体用于将所述波阻抗反演体中位于所述储层门槛值以外的阻抗值剔除。
所述孔隙度反演模块，具体用于将处理后的波阻抗反演体转换为孔隙度反演体。
在地震勘探技术领域中，岩石储层物性反演是近年来发展起来的一种反演方法，其目的是更为直接地把地震数据与地质认识结合起来，可以更为直观地为地质人员提供储层解释依据及油气判别依据。常用的物性反演包括电阻率反演、伽马反演、孔隙度反演、渗透率、饱和度以及其他岩性参数反演。其中，孔隙度反演的主要作用是进行储层量化预测，适用于评价及开发阶段。因此，将处理后的波阻抗反演体转换为孔隙度反演体可更为直观的反映储层信息，便于对储层进行评价。
在本实施例中，解决了现有技术中无法准确对生物礁气藏储层进行储层预测的问题，整个过程减少了相变线的人工定义，避免了人工定义引起的误差。充分利用了地震资料中的层速度体，进行了数据体层面的储层预测，避免了现有技术中以一个平面的数据趋势控制一个体的建模存在的误差，从而导致储层刻画不准确的问题，实现了对生物礁气藏储层的精细刻画，并提高了生物礁气藏开发经济效益。
实施例四
本实施例还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，所述计算机程序被处理器执行时可以实现如下方法步骤：
步骤S110：获取待测井点的层速度体，其中，所述层速度体包括地震波在每层地层的传播速度。
步骤S120：采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体。
步骤S130：根据所述波阻抗体，进行波阻抗反演，得到波阻抗反演体。
步骤S140：获取电阻率测井曲线。
步骤S150：根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体。
步骤S160：根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
上述方法步骤的具体实施例过程可参见实施例一，本实施例在此不再重复赘述。
实施例五
本实施例提供了一种电子设备，该电子设备可以是手机、电脑或平板电脑等，包括存储器和处理器，所述存储器上存储有计算器程序，该计算机程序被处理器执行时实现如实施例一中所述的储层预测方法。可以理解，电子设备还可以包括，多媒体组件，输入/输出(I/O)接口。
其中，处理器用于执行如实施例一中的储层预测方法中的全部或部分步骤。存储器用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
所述处理器可以是专用集成电路(Application Specific Integrated Circuit，简称ASIC)、数字信号处理器(Digital Signal Processor，简称DSP)、数字信号处理设备(Digital Signal Processing Device，简称DSPD)、可编程逻辑器件(Programmable LogicDevice，简称PLD)、现场可编程门阵列(Field Programmable Gate Array，简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述实施例一中的储层预测方法。
所述存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(Static Random Access Memory，简称SRAM)，电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory，简称EEPROM)，可擦除可编程只读存储器(Erasable Programmable Read-Only Memory，简称EPROM)，可编程只读存储器(Programmable Read-Only Memory，简称PROM)，只读存储器(Read-Only Memory，简称ROM)，磁存储器，快闪存储器，磁盘或光盘。
多媒体组件可以包括屏幕，所述屏幕可以是触摸屏，用于显示各个剖面示意图，例如，层速度体连井剖面示意图。
I/O接口为处理器和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。
综上，本申请提供的一种储层预测方法、装置、存储介质及电子设备，所述方法包括：获取待测井点的层速度体，其中，所述层速度体包括地震波在每层地层的传播速度；采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；根据所述波阻抗体，进行波阻抗反演，得到波阻抗反演体；获取电阻率测井曲线；根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。解决了现有技术中无法准确对生物礁气藏储层进行储层预测的问题，避免了现有技术中以一个平面的数据趋势控制一个体的建模存在的误差，从而导致储层刻画不准确的问题，实现了对生物礁气藏储层的精细刻画。
进一步可以理解的是，为了提高储层预测的结果，可对获取到的测井曲线进行预处理。测井曲线预处理是测井解释与数据处理的一项基础工作，它是保证测井解释与数据处理结果精度的重要前提。
进一步可以理解的是，由于孔隙度反演的主要作用是进行储层量化预测，适用于各种气藏类型的评价及开发阶段。因此，将处理后的波阻抗反演体转换为孔隙度反演体可更为直观的反映生物礁气藏储层的储层信息，便于对储层进行评价。
在本申请实施例所提供的几个实施例中，应该理解到，所揭露的电子设备和方法，也可以通过其它的方式实现。以上所描述的电子设备和方法实施例仅仅是示意性的。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
虽然本申请所揭露的实施方式如上，但所述的内容只是为了便于理解本申请而采用的实施方式，并非用以限定本申请。任何本申请所属技术领域内的技术人员，在不脱离本申请所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本申请的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种储层预测方法，其特征在于，所述方法包括：
获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的传播速度；
采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；
根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体；
获取电阻率测井曲线；
根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；
根据所述电阻率数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
2.如权利要求1所述的方法，其特征在于，所述获取待测井点的层速度体，包括：
采用迪克斯计算式得到所述地层中每层的层速度，其中，所述迪克斯计算式为：
Figure FDA0002524436040000011
其中，νn为所述地层中第n层的层速度，t1，n为所述地层中第一层到第n层的地震双程垂直反射时间，νR，n为所述地层中第一层到第n层的均方根速度，t1，n-1为所述地层中第一层到第n-1层的地震双程垂直反射时间，νR，n-1为所述地层中第一层到第n-1层的均方根速度。
3.如权利要求1所述的方法，其特征在于，所述层速度-阻抗关系式为：
Imp＝a*f(V)+b；
其中，Imp为波阻抗体，f(V)为层速度体，a为第一预设参数，b为第二预设参数。
4.如权利要求1所述的方法，其特征在于，根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体，包括：
根据所述波阻抗体和所述测井数据，建立低频模型；
对所述低频模型进行波阻抗反演，以得到波阻抗反演体。
5.如权利要求1所述的方法，其特征在于，根据所述电阻率模拟数据体对所述波阻抗反演体进行处理，包括：
根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值；
将所述波阻抗反演体中位于所述储层门槛值以外的阻抗值剔除。
6.如权利要求5所述的方法，其特征在于，根据所述电阻率模拟数据体和所述波阻抗反演体，得到储层门槛值，包括：
将所述电阻率模拟数据体和所述波阻抗反演体进行交会分析处理，以得到储层门槛值。
7.如权利要求1所述的方法，其特征在于，将处理后的波阻抗反演体转换为孔隙度反演体，包括：
建立波阻抗反演体和孔隙度的线性关系；
根据建立的所述线性关系，将处理后的波阻抗反演体转换为孔隙度反演体。
8.一种储层预测装置，其特征在于，所述装置包括：
第一获取模块，用于获取待测井点的层速度体和测井数据，其中，所述层速度体包括地震波在每层地层的层速度；
第一转换模块，用于采用层速度-阻抗关系式，将所述层速度体转换为波阻抗体；
反演模块，用于根据所述波阻抗体和所述测井数据，进行波阻抗反演，得到波阻抗反演体；
第二获取模块，用于获取电阻率测井曲线；
波形指示模拟模块，用于根据所述波阻抗反演体和所述电阻率测井曲线，进行波形指示模拟，以得到电阻率模拟数据体；
第二转换模块，用于根据所述电阻率模拟数据体对所述波阻抗反演体进行处理，并将处理后的波阻抗反演体转换为孔隙度反演体，以从所述孔隙度反演体中提取参数用于对储层的相关特性进行预测。
9.一种存储介质，其特征在于，该存储介质存储有计算机程序，当该计算机程序被一个或多个处理器执行，实现如权利要求1-7中任意一项所述的储层预测方法。
10.一种电子设备，其特征在于，包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，实现如权利要求1-7中任意一项所述的储层预测方法。
CN202010500243.7A 2020-06-04 2020-06-04 一种储层预测方法、装置、存储介质及电子设备 Active CN113759419B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010500243.7A CN113759419B (zh) 2020-06-04 2020-06-04 一种储层预测方法、装置、存储介质及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010500243.7A CN113759419B (zh) 2020-06-04 2020-06-04 一种储层预测方法、装置、存储介质及电子设备

Publications (2)

Publication Number Publication Date
CN113759419A true CN113759419A (zh) 2021-12-07
CN113759419B CN113759419B (zh) 2024-06-18

Family

ID=78783704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010500243.7A Active CN113759419B (zh) 2020-06-04 2020-06-04 一种储层预测方法、装置、存储介质及电子设备

Country Status (1)

Country Link
CN (1) CN113759419B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114839671A (zh) * 2022-03-28 2022-08-02 中国矿业大学 地面和地巷瞬变电磁法精细识别煤系地层电性界面方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000000849A2 (en) * 1998-06-18 2000-01-06 Den Norske Stats Oljeselskap A.S Device and method for measurement by guided waves on a metal string in a well
CN101634717A (zh) * 2009-08-26 2010-01-27 中国石油大学(华东) 基于测井和叠前道集地震数据的精细横波阻抗求取技术
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN104516018A (zh) * 2013-09-30 2015-04-15 中国石油化工股份有限公司 一种地球物理勘探中岩性约束下的孔隙度反演方法
CN104570067A (zh) * 2013-10-11 2015-04-29 中国石油化工股份有限公司 一种地球物理勘探中相控地震反演方法
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN104914470A (zh) * 2014-03-14 2015-09-16 恒泰艾普石油天然气技术服务股份有限公司 碳酸盐岩缝洞储层储量校正方法
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
WO2017024702A1 (zh) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 一种射线弹性参数的反演系统
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN108072903A (zh) * 2016-11-09 2018-05-25 中国石油化工股份有限公司 一种测井曲线重构方法
WO2019062655A1 (zh) * 2017-09-27 2019-04-04 中国石油天然气股份有限公司 薄夹层的确定方法和装置
CN111007567A (zh) * 2018-10-08 2020-04-14 中国石油化工股份有限公司 基于地震波形反演的砂泥岩薄互层预测方法及系统

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000000849A2 (en) * 1998-06-18 2000-01-06 Den Norske Stats Oljeselskap A.S Device and method for measurement by guided waves on a metal string in a well
CN101634717A (zh) * 2009-08-26 2010-01-27 中国石油大学(华东) 基于测井和叠前道集地震数据的精细横波阻抗求取技术
CN101980053A (zh) * 2010-10-18 2011-02-23 中国石油化工股份有限公司 一种复杂礁滩储层预测方法
CN104516018A (zh) * 2013-09-30 2015-04-15 中国石油化工股份有限公司 一种地球物理勘探中岩性约束下的孔隙度反演方法
CN104570067A (zh) * 2013-10-11 2015-04-29 中国石油化工股份有限公司 一种地球物理勘探中相控地震反演方法
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN104914470A (zh) * 2014-03-14 2015-09-16 恒泰艾普石油天然气技术服务股份有限公司 碳酸盐岩缝洞储层储量校正方法
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
WO2017024702A1 (zh) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 一种射线弹性参数的反演系统
CN108072903A (zh) * 2016-11-09 2018-05-25 中国石油化工股份有限公司 一种测井曲线重构方法
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
WO2019062655A1 (zh) * 2017-09-27 2019-04-04 中国石油天然气股份有限公司 薄夹层的确定方法和装置
CN111007567A (zh) * 2018-10-08 2020-04-14 中国石油化工股份有限公司 基于地震波形反演的砂泥岩薄互层预测方法及系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114839671A (zh) * 2022-03-28 2022-08-02 中国矿业大学 地面和地巷瞬变电磁法精细识别煤系地层电性界面方法

Also Published As

Publication number Publication date
CN113759419B (zh) 2024-06-18

Similar Documents

Publication Publication Date Title
Corbeanu et al. Detailed internal architecture of a fluvial channel sandstone determined from outcrop, cores, and 3-D ground-penetrating radar: Example from the middle Cretaceous Ferron Sandstone, east-central Utah
CN109425896A (zh) 白云岩油气储层分布预测方法及装置
CN111596364B (zh) 基于高精度层序地层格架的地震沉积微相组合分析方法
Ali et al. Prediction of Cretaceous reservoir zone through petrophysical modeling: Insights from Kadanwari gas field, Middle Indus Basin
Fajana et al. Application of multilayer perceptron neural network and seismic multiattribute transforms in reservoir characterization of Pennay field, Niger Delta
CN112114380A (zh) 一种用于砂砾岩油藏储层的预测方法
Lu et al. Prestack inversion identification of organic reef gas reservoirs of Permian Changxing Formation in Damaoping area, Sichuan Basin, SW China
Hassan et al. An advanced workflow to compress the uncertainties of stochastic distribution of Bahariya reservoir properties using 3D static modeling: an example from Heba Oil Fields, Western Desert, Egypt
MacBeth et al. The 4D seismic signature of oil–water contact movement due to natural production in a stacked turbidite reservoir
Sanda et al. The integrated approach to seismic attributes of lithological characterization of reservoirs: case of the F3 Block, North Sea-Dutch Sector
CN113759419B (zh) 一种储层预测方法、装置、存储介质及电子设备
CN111077578B (zh) 岩层分布预测方法和装置
Okpoli et al. Integration of Well logs and seismic attribute analysis in reservoir identification on PGS field onshore Niger Delta, Nigeria
Paris et al. Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data
Aminzadeh et al. Fundamentals of Petroleum Geophysics
Alshakhs Shale play assessment of the Goldwyer formation in the Canning basin using property modelling
Abu-El Ata et al. Geologic Causatives and Seismic Attributes Heterogeneities for Structural and Stratigraphic Interplays, Sitra Area, North Western Desert, Egypt
Botter et al. Seismic attribute analysis of a fault zone in the Thebe field, Northwest shelf, Australia
Azar et al. Integrated seismic attributes to characterize a widely distributed carbonate clastic deposit system in Khuzestan Province, SW Iran
Alam Estimation of reservoir heterogeneity from the depositional environment in reservoir characterization of a CHOPS Field
Asgari Nezhad et al. Facies Quality Zoning in Shale Gas by Deep Learning Method
Braaksma et al. The origin and nature of seismic reflections of sharp‐based shoreface deposits (upper Jurassic Siliciclastics, northern France)
Ayodele An integrated study of the early cretaceous (Valanginian) reservoir from the Gamtoos Basin, offshore South Africa with special reference to seismic cacies, formation evaluation and static reservoir modeling
Hassan et al. Petroleum Research
Aminzadeh et al. Geophysics for petroleum engineers: Chapter 3. Fundamentals of petroleum geophysics

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