CN113917529A - 波阻抗反演方法、装置、存储介质及电子设备 - Google Patents

波阻抗反演方法、装置、存储介质及电子设备 Download PDF

Info

Publication number
CN113917529A
CN113917529A CN202010646937.1A CN202010646937A CN113917529A CN 113917529 A CN113917529 A CN 113917529A CN 202010646937 A CN202010646937 A CN 202010646937A CN 113917529 A CN113917529 A CN 113917529A
Authority
CN
China
Prior art keywords
wave impedance
impedance model
wave
model
wave number
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
CN202010646937.1A
Other languages
English (en)
Other versions
CN113917529B (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 CN202010646937.1A priority Critical patent/CN113917529B/zh
Publication of CN113917529A publication Critical patent/CN113917529A/zh
Application granted granted Critical
Publication of CN113917529B publication Critical patent/CN113917529B/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
    • 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

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)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明涉及油气勘探技术领域，具体涉及一种波阻抗反演方法、装置、存储介质及电子设备，获取宽频带波数，宽频带波数包括低频波数、中频波数以及高频波数；根据低频波数以及初始波阻抗模型，得到第一波阻抗模型；根据中频波数以及第一波阻抗模型，得到第二波阻抗模型；根据高频波数以及第二波阻抗模型，得到第三波阻抗模型，将采集到的地层信息输入第三波阻抗模型进行反演，得到波阻抗反演结果，地层信息表征地层的岩性信息；将地层信息输入到第三波阻抗模型中进行反演，能够反演出分辨率较高的波阻抗反演结果，能够真实的反映储层情况。

Description

波阻抗反演方法、装置、存储介质及电子设备
技术领域
本发明涉及油气勘探技术领域，特别地涉及一种波阻抗反演方法、装置、存储介质及电子设备。
背景技术
储层地质模型即建立储层特征三维分布的数字化模型，其被广泛的应用在地质勘探技术领域，储层建模目的是通过在油气勘探和开发的过程中取得地震、测井、钻井等方面的资料，对储层的各个方面的特性进行描述，达到储层评测的目的。
在现有技术的储层建模在反演的过程中，都涉及多个参数，通过多个参数同时反演，得到储层模型。
在现有技术的反演过程中，由于目标函数对每个参数的敏感性不同，以及每个参数之间的耦合作用等因素，导致多个参数同时反演的难度较大，且不易同时反演出满足要求的波阻抗参数，从而导致现有技术中反演得到的波阻抗参数的准确性得不得保证，不能反映储层的真实情况。
发明内容
针对上述问题，本发明提供一种波阻抗反演方法、装置、存储介质及电子设备，以提高反演得到的波阻抗参数的准确性，能够反映储层的真实情况。
第一方面，本发明提供了一种波阻抗反演方法，所述方法包括：
获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数；
根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型；
根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型；
根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型；
将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
结合第一方面，在第一种可能的实现方式中，所述获取宽频带波数的步骤，包括：
根据地震资料获取所述中频波数；
采集测井内的所述低频波数以及所述高频波数。
结合第一方面，在第二种可能的实现方式中，所述根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型的步骤，包括：
将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果；
计算所述第一输出结果与第一测量结果之间的差值，得到第一残差；
判断所述第一残差是否大于第一预设残差；
若所述第一残差大于所述第一预设残差，则根据所述第一测量结果、所述第一输出结果以及第一对应关系，得到第一模型参数，其中，所述第一对应关系为所述第一测量结果、所述第一输出结果以及所述第一模型参数之间的对应关系；
根据所述第一模型参数调节所述初始波阻抗模型，得到更新后的初始波阻抗模型后，用所述更新后的初始波阻抗模型替代所述初始波阻已修改抗模型，执行所述将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果的步骤，直至所述第一残差不大于所述第一预设残差，得到第一波阻抗模型。
结合第二种可能的实现方式，在第三种可能的实现方式中，所述根据所述第一模型参数调节所述初始波阻抗模型的步骤，包括：
将所述第一模型参数转换为第一三角函数；
将所述第一三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第一调节系数；
根据所述第一调节系数调节所述初始波阻抗模型。
结合第一方面，在第四种可能的实现方式中，所述根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型的步骤，包括：
将所述中频波数输入所述第一波阻抗模型中，得到第二输出结果；
计算所述第二输出结果与第二测量结果之间的差值，得到第二残差；
判断所述第二残差是否大于第二预设残差；
若所述第二残差大于所述第二预设残差，则根据所述第二测量结果、所述第二输出结果以及第二对应关系，得到第二模型参数，其中，所述第二对应关系为所述第二测量结果、所述第二输出结果以及所述第二模型参数之间的对应关系；
根据所述第二模型参数调节所述第一波阻抗模型，得到更新后的第一波阻抗模型后，用所述更新后的第一波阻抗模型替代所述第一波阻抗模型，执行所述将所述低频波数输入所述第一波阻抗模型中，得到第二输出结果的步骤，直至所述第二残差不大于所述第二预设残差，得到第二波阻抗模型。
结合第一方面，在第五种可能的实现方式中，所述根据所述第二模型参数调节所述第一波阻抗模型的步骤，包括：
将所述第二模型参数转换为第二三角函数；
将所述第二三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第二调节系数；
根据所述第二调节系数调节所述第一波阻抗模型。
结合第一方面，在第六种可能的实现方式中，所述根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型的步骤，包括：
将所述高频波数输入所述第二波阻抗模型中，得到第三输出结果；
计算所述第三输出结果与第三测量结果之间的差值，得到第三残差；
判断所述第三残差是否大于第三预设残差；
若所述第三残差大于所述第三预设残差，则根据所述第三测量结果、所述第三输出结果以及第单对应关系，得到第三模型参数，其中，所述第三对应关系为所述第三测量结果、所述第三输出结果以及所述第三模型参数之间的对应关系；
根据所述第三模型参数调节所述第二波阻抗模型，得到更新后的第二波阻抗模型后，用所述更新后的第二波阻抗模型替代所述第二波阻抗模型，执行所述将所述低频波数输入所述第二波阻抗模型中，得到第三输出结果的步骤，直至所述第三残差不大于所述第三预设残差，得到第三波阻抗模型。
结合第六种可能的实现方式，在第七种可能的实现方式中，所述根据所述第三模型参数调节所述第二波阻抗模型的步骤，包括：
将所述第三模型参数转换为第三三角函数；
将所述第三三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第三调节系数；
根据所述第三调节系数调节所述第三波阻抗模型。
第二方面，本发明提供了一种波阻抗反演装置，包括：获取模块、第一处理模块、第二处理模块、第三处理模块和第四处理模块；
所述获取模块，用于获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数；
所述第一处理模块，用于根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型；
所述第二处理模块，用于根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型；
所述第三处理模块，用于根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型；
所述第四处理模块，用于将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
结合第二方面，在第一种可能的实现方式中，获取模块，具体用于根据地震资料获取所述中频波数；
采集测井内的所述低频波数以及所述高频波数。
结合第二方面，在第二种可能的实现方式中，所述第一处理模块，具体用于将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果；
计算所述第一输出结果与第一测量结果之间的差值，得到第一残差；
判断所述第一残差是否大于第一预设残差；
若所述第一残差大于所述第一预设残差，则根据所述测量结果、所述第一输出结果以及第一对应关系，得到第一模型参数，其中，所述第一对应关系为所述测量结果、所述第一输出结果以及所述第一模型参数之间的对应关系；
根据所述第一模型参数调节所述初始波阻抗模型，得到更新后的初始波阻抗模型后，用所述更新后的初始波阻抗模型替代所述初始波阻已修改抗模型，执行所述将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果的步骤，直至所述第一残差不大于所述第一预设残差，得到第一波阻抗模型。
结合第二种可能的实现方式，在第三种可能的实现方式中，所述第一处理模块，还用于将所述第一模型参数转换为第一三角函数；
将所述第一三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第一调节系数；
根据所述第一调节系数调节所述初始波阻抗模型。
结合二方面，在第四种可能的实现方式中，所述第二处理模块，具体用于：
将所述中频波数输入所述第一波阻抗模型中，得到第二输出结果；
计算所述第二输出结果与第二测量结果之间的差值，得到第二残差；
判断所述第二残差是否大于第二预设残差；
若所述第二残差不大于所述第二预设残差，则根据所述测量结果、所述第二输出结果以及第二对应关系，得到第二模型参数，其中，所述第二对应关系为所述测量结果、所述第二输出结果以及所述第二模型参数之间的对应关系；
根据所述第二模型参数调节所述第一波阻抗模型，得到更新后的第一波阻抗模型后，用所述更新后的第一波阻抗模型替代所述第一波阻抗模型，执行所述将所述低频波数输入所述第一波阻抗模型中，得到第二输出结果的步骤，直至所述第二残差不大于所述第二预设残差，得到第二波阻抗模型。
结合第四种可能的实现方式，在第五种可能的实现方式中，所述第二处理模块，还用于将所述第二模型参数转换为第二三角函数；
将所述第二三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第二调节系数；
根据所述第二调节系数调节所述第一波阻抗模型。
结合第二方面，在第六种可能的实现方式中，所述第三处理模块，具体用于将所述高频波数输入所述第二波阻抗模型中，得到第三输出结果；
计算所述第三输出结果与第三测量结果之间的差值，得到第三残差；
判断所述第三残差是否大于第三预设残差；
若所述第三残差不大于所述第三预设残差，则根据所述测量结果、所述第三输出结果以及第单对应关系，得到第三模型参数，其中，所述第三对应关系为所述测量结果、所述第三输出结果以及所述第三模型参数之间的对应关系；
根据所述第三模型参数调节所述第二波阻抗模型，得到更新后的第二波阻抗模型后，用所述更新后的第二波阻抗模型替代所述第二波阻抗模型，执行所述将所述低频波数输入所述第二波阻抗模型中，得到第三输出结果的步骤，直至所述第三残差不大于所述第三预设残差，得到第三波阻抗模型。
结合第六种可能的实现方式，在第七种可能的实现方式中，所述第三处理模块，还用于将所述第三模型参数转换为第三三角函数；
将所述第三三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第三调节系数；
根据所述第三调节系数调节所述第三波阻抗模型。
第三方面，本发明提供一种存储介质，所述存储介质存储有计算机程序，所述存储介质被一个或多个处理器执行时，实现如第一方面所述的波阻抗反演方法。
第四方面，本发明提供一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，执行如第一方面所述的波阻抗反演方法。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
本发明提供的一种波阻抗反演方法、装置、存储介质及电子设备，获取宽频带波数，其中，宽频带波数包括低频波数、中频波数以及高频波数；根据低频波数以及初始波阻抗模型，得到第一波阻抗模型；根据中频波数以及第一波阻抗模型，得到第二波阻抗模型；根据高频波数以及第二波阻抗模型，得到第三波阻抗模型，将采集到的地层信息输入第三波阻抗模型进行反演，得到波阻抗反演结果，其中，地层信息表征地层的岩性信息；依次经过低频波数、中频波数以及高频波数得到的第三波阻抗模型，第三波阻抗模型能够准确的反映储层的真实情况，将地层信息输入到第三波阻抗模型中进行反演，能够反演出分辨率较高的波阻抗反演结果，能够真实的反映储层情况。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图，具体地，在下文中将基于实施例并参考附图来对本发明进行更详细的描述：
图1为本发明实施例提供的一种波阻抗反演方法的流程示意图；
图2为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图3为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图4为本发明实施例提供的第一波阻抗模型的反演结果；
图5为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图6为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图7为本发明实施例提供的第二波阻抗模型的反演结果；
图8为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图9为本发明实施例提供的一种波阻抗反演方法的另一流程示意图；
图10为本发明实施例提供的第三波阻抗模型的反演结果；
图11为本发明实施例提供的一种波阻抗反演装置的连接框图。
附图标记：1-波阻抗反演装置；101-获取模块；102-第一处理模块；103-第二处理模块；104-第三处理模块；105-第四处理模块。
在附图中，相同的部件使用相同的附图标记，附图并未按照实际的比例绘制。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚、明白，下面，将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本申请的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
在本申请的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本申请中的具体含义。
下面结合附图，对本申请的一些实施方式作详细说明。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
实施例一
针对现有技术中存在的上述技术问题，本实施例提供一种波阻抗反演方法，图1为本发明实施例提供的一种波阻抗反演方法的流程示意图，需要说明的是，本发明实施例提供的波阻抗反演方法并不以图1以及以下的具体顺序为限制，应当理解，在其它实施例中，本发明实施例提供的波阻抗反演方法其中部分步骤的顺序可以根据实际需要相互交换，或者其中的部分步骤也可以省略或删除。下面将对图1涉及到的具体流程进行阐述，如图1所示，波阻抗反演方法包括如下步骤：
步骤S1、获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数。
步骤S1具体包括如下子步骤：
子步骤S11、根据地震资料获取所述中频波数。
子步骤S12、采集测井内的所述低频波数以及所述高频波数。
由于地球的吸收和衰减作用，导致获取的地震子波中通常只含有低频波数和中频波数成分，又由于地震检波器自身设计的缘故，地震检波器通常只能检测出8HZ以上的波数，可以理解的是，在检测的过程中过滤了低频波数，因此，在采集的地震资料中通常只提供波阻抗参数的中频波数，因此，需要结合测井内采集的低频波数和高频波数以补全频带，得到宽频带波数，以提高储层预测的分辨率。
步骤S2、根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型。
通过低频波数得到的第一波阻抗模型，能够准确反映低频波数的对应的储层情况。
步骤S3、根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型。
通过中频波数得到的第二波阻抗模型，能够准确反映中频波数的对应的储层情况。
步骤S4、根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型。
通过高频波数得到的第三波阻抗模型，能够准确反映高频波数的对应的储层情况。
在步骤S2和步骤S3中，通过低频波数和中频波数构建储层的大尺度的结构，即整体架构，再通过高频波数完善储层的细节信息，使得得到的第三波阻抗模型能够更真实的反映储层情况。
步骤S5、将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
可选地，地层信息可以包括下述一项或多项的任意组合：地层层位，波阻抗初值，测井速度、测井密度。
本实施例提供的一种波阻抗反演方法，获取宽频带波数，其中，宽频带波数包括低频波数、中频波数以及高频波数；根据低频波数以及初始波阻抗模型，得到第一波阻抗模型；根据中频波数以及第一波阻抗模型，得到第二波阻抗模型；根据高频波数以及第二波阻抗模型，得到第三波阻抗模型，将采集到的地层信息输入第三波阻抗模型进行反演，得到波阻抗反演结果，其中，地层信息表征地层的岩性信息；依次经过低频波数、中频波数以及高频波数得到的第三波阻抗模型，第三波阻抗模型能够准确的反映储层的真实情况，将地层信息输入到第三波阻抗模型中进行反演，能够反演出分辨率较高的波阻抗反演结果，能够真实的反映储层情况。
实施例二
在实施例一的基础上，本实施例提供另一种波阻抗反演方法，具体地，图2为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图2所示，步骤S2包括如下子步骤：
子步骤S21、将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果。
子步骤S22、计算所述第一输出结果与第一测量结果之间的差值，得到第一残差。
测量结果为储层对应的实际采集结果，其包括低频波数对应的第一测量结果、中频波数对应的第二测量结果、以及高频波数对应的第三测量结果。在本实施例中，采用的低频波数得到的第一输出结果应该与第一测量结果进行比较，因此，为了保证得到的残差的准确性，采用滤波的方式将测量结果中的第二测量结果与第三测量结果过滤。可选地，滤波器为低通滤波器。利用低通滤波器将测量结果进行过滤，得到第一测量结果。
子步骤S23、判断所述第一残差是否大于第一预设残差。
如果第一残差不大于第一预设残差(即第一残差小于或等于第一预设残差)，则说明初始波阻抗模型的参数合理，初始波阻抗模型能够真实的反映低频波数对应的储层信息，流程进入步骤S24；如果第一残差大于第一预设残差，则说明初始波阻抗的参数不合理，初始波阻抗模型的分辨率较低，不能真实的反映低频波数对应的储层信息，流程进入步骤S25。
子步骤S24、确定所述初始波阻抗模型为所述第一阻抗模型。
子步骤S25、根据所述第一测量结果、所述第一输出结果以及第一对应关系，得到第一模型参数，其中，所述第一对应关系为所述第一测量结果、所述第一输出结果以及所述第一模型参数之间的对应关系。
其中，第一对应关系可以通过下式表达：
Figure BDA0002573402780000101
Figure BDA0002573402780000102
其中，
Figure BDA0002573402780000103
为中间参数，为了简化公式而设置的，P、p、Q、q、H、h均为累加系数，
Figure BDA0002573402780000104
为第一测量结果，
Figure BDA0002573402780000105
为第一输出结果，Δt为采样时长，m1为第一模型参数，almn为傅里叶系数，Δu1pq为第一残差，T为矩阵的转置。
在本步骤中得到的第一模型参数m1是关于第一测量结果、所述第一输出结果的表达式。
子步骤S26、根据所述第一模型参数调节所述初始波阻抗模型，得到更新后的初始波阻抗模型后，用所述更新后的初始波阻抗模型替代所述初始波阻已修改抗模型，执行所述将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果的步骤，直至所述第一残差不大于所述第一预设残差，得到第一波阻抗模型。
所述根据所述第一模型参数调节所述初始波阻抗模型的步骤，图3为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图3所示，该步骤包括如下子步骤：
子步骤S261、将所述第一模型参数转换为第一三角函数，如cos(narccosf)。
其中，n为系数，f为频率。
子步骤S262、将所述第一三角函数转换为傅里叶级数展开式。
傅里叶展开式如下所示：
Figure BDA0002573402780000111
其中，X0(f)和X1(f)为傅里叶展开式的前两项，Xn-1(f)、Xn(f)和Xn+1(f)为傅里叶展开式的最后三项。
子步骤S263、根据所述傅里叶级数展开式的系数，得到第一调节系数。
傅里叶级数展开式的系数如下所示：
Figure BDA0002573402780000112
Figure BDA0002573402780000113
其中，j、k、N均为参数，j的取值范围为[1,N]之间的整数；c0、cj1均为傅里叶级数展开式的系数，其中，cj1为第一调节系数，f(xk)为傅里叶展开式对应的函数表达式，xk为傅里叶展开式中的一个，π为圆周率。
子步骤S264、根据所述第一调节系数调节所述初始波阻抗模型。
第一调节系数cj1是满足低频波数对应的系数，根据第一调节系数cj1调节初始波阻抗模型，使得调节后的初始波阻抗模型能够更清楚的反映低频波数代表的储层情况。
图4为本发明实施例提供的第一波阻抗模型的反演结果，将地层信息输入第一波阻抗模型，得到的反演结果如图4所示，其中，横轴为D轴，表示距离，单位为千米。纵轴为T轴，表示时间，单位为秒，秒数与深度成正比，秒数越大，则深度越大；秒数越小，则深度越小。T轴对应的数值可以通过时深转换，转换成深度值。
实施例三
在上述实施例的基础上，本实施例提供一种波阻抗反演方法，具体地，图5为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图5所示，步骤S3包括如下子步骤：
子步骤S31、将所述中频波数输入所述第一波阻抗模型中，得到第二输出结果。
子步骤S32、计算所述第二输出结果与第二测量结果之间的差值，得到第二残差。
在本实施例中，采用的中频波数得到的第二输出结果应该与第二测量结果进行比较，因此，为了保证得到的残差的准确性，采用滤波的方式将测量结果中的第一测量结果与第三测量结果过滤。可选地，滤波器为带通滤波器。利用带通滤波器将测量结果进行过滤，得到第二测量结果。
子步骤S33、判断所述第二残差是否大于第二预设残差。
如果第二残差不大于第二预设残差(即第二残差小于或等于第二预设残差)，则说明第一波阻抗模型的参数合理，第一波阻抗模型能够真实的反映低频波数对应的储层信息，流程进入步骤S34；如果第二残差大于第二预设残差，则说明第一波阻抗的参数不合理，第一波阻抗模型的分辨率较低，不能真实的反映中频波数对应的储层信息，流程进入步骤S35。
子步骤S34、确定所述第一波阻抗模型为第二波阻抗模型。
子步骤S35、根据所述第二测量结果、所述第二输出结果以及第二对应关系，得到第二模型参数，其中，所述第二对应关系为所述第二测量结果、所述第二输出结果以及所述第二模型参数之间的对应关系。
其中，第二对应关系可以通过下式表达：
Figure BDA0002573402780000121
Figure BDA0002573402780000122
其中，
Figure BDA0002573402780000123
为中间参数，为了简化公式而设置的，P、p、Q、q、H、h均为累加系数，
Figure BDA0002573402780000124
为第二测量结果，
Figure BDA0002573402780000125
为第二输出结果，Δt为采样时长，m2为第二模型参数，almn为傅里叶系数，Δu2pq为第二残差，T为矩阵的转置。
在本步骤中得到的第二模型参数m2是关于第二测量结果、所述第二输出结果的表达式。
子步骤S36、根据所述第二模型参数调节所述第一波阻抗模型，得到更新后的第一波阻抗模型后，用所述更新后的第一波阻抗模型替代所述第一波阻抗模型，执行所述将所述低频波数输入所述第一波阻抗模型中，得到第二输出结果的步骤，直至所述第二残差不大于所述第二预设残差，得到第二波阻抗模型。
图6为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图6所示，所述根据所述第二模型参数调节所述第一波阻抗模型的步骤，包括：
子步骤S361、将所述第二模型参数转换为第二三角函数。
子步骤S362、将所述第二三角函数转换为傅里叶级数展开式。
子步骤S363、根据所述傅里叶级数展开式的系数，得到第二调节系数。
子步骤S361-S363涉及到的方法与上一实施例中子步骤S261-S263类似，在此不再赘述。
子步骤S364、根据所述第二调节系数调节所述第一波阻抗模型。
第二调节系数是满足中频波数对应的系数，根据第二调节系数调节第一波阻抗模型，使得调节后的第一波阻抗模型能够更清楚的反映中频波数代表的储层情况。
图7为本发明实施例提供的第二波阻抗模型的反演结果，将地层信息输入第二波阻抗模型，得到的反演结果如图7所示，其中，横轴为D轴，表示距离，单位为千米。纵轴为T轴，表示时间，单位为秒，秒数与深度成正比，秒数越大，则深度越大；秒数越小，则深度越小。T轴对应的数值可以通过时深转换，转换成深度值。
实施例四
在上述实施例的基础上，本实施例提供了一种波阻抗反演方法，具体地，图8为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图8所示，步骤S4包括如下子步骤：
子步骤S41、将所述高频波数输入所述第二波阻抗模型中，得到第三输出结果。
子步骤S42、计算所述第三输出结果与第三测量结果之间的差值，得到第三残差。
在本实施例中，采用的高频波数得到的第三输出结果应该与第三测量结果进行比较，因此，为了保证得到的残差的准确性，采用滤波的方式将测量结果中的第一测量结果与第二测量结果过滤。可选地，滤波器为高通滤波器。利用高通滤波器将测量结果进行过滤，得到第三测量结果。
子步骤S43、判断所述第三残差是否大于第三预设残差。
如果第三残差不大于第三预设残差(即第三残差小于或等于第三预设残差)，则说明第二波阻抗模型的参数合理，第二波阻抗模型能够真实的反映高频波数对应的储层信息，流程进入步骤S44；如果第三残差大于第三预设残差，则说明第二波阻抗的参数不合理，第二波阻抗模型的分辨率较低，不能真实的反映高频波数对应的储层信息，流程进入步骤S45。
子步骤S44、确定所述第二波阻抗模型为第三波阻抗模型。
子步骤S45、根据所述第三测量结果、所述第三输出结果以及第单对应关系，得到第三模型参数，其中，所述第三对应关系为所述第三测量结果、所述第三输出结果以及所述第三模型参数之间的对应关。
其中，第三对应关系可以通过下式表达：
Figure BDA0002573402780000141
Figure BDA0002573402780000142
其中，
Figure BDA0002573402780000143
为中间参数，为了简化公式而设置的，P、p、Q、q、H、h均为累加系数，
Figure BDA0002573402780000144
为第三测量结果，
Figure BDA0002573402780000145
为第三输出结果，Δt为采样时长，m3为第三模型参数，almn为傅里叶系数，Δu3pq为第三残差，T为矩阵的转置。
在本步骤中得到的第三模型参数m3是关于第三测量结果、所述第三输出结果的表达式。
子步骤S46、根据所述第三模型参数调节所述第二波阻抗模型，得到更新后的第二波阻抗模型后，用所述更新后的第二波阻抗模型替代所述第二波阻抗模型，执行所述将所述低频波数输入所述第二波阻抗模型中，得到第三输出结果的步骤，直至所述第三残差不大于所述第三预设残差，得到第三波阻抗模型。
9为本发明实施例提供的一种波阻抗反演方法的另一流程示意图，如图9所示，所述根据所述第三模型参数调节所述第二波阻抗模型的步骤，包括：
子步骤S461、将所述第三模型参数转换为第三三角函数。
子步骤S462、将所述第三三角函数转换为傅里叶级数展开式。
子步骤S463、根据所述傅里叶级数展开式的系数，得到第三调节系数。
子步骤S461-S463涉及到的方法与上一实施例中子步骤S261-S263类似，在此不再赘述。
子步骤S464、根据所述第三调节系数调节所述第三波阻抗模型。
第二调节系数是满足中频波数对应的系数，根据第二调节系数调节第一波阻抗模型，使得调节后的第一波阻抗模型能够更清楚的反映中频波数代表的储层情况。
图10为本发明实施例提供的第二波阻抗模型的反演结果，将地层信息输入第二波阻抗模型，得到的反演结果如图10所示，其中，横轴为D轴，表示距离，单位为千米。纵轴为T轴，表示时间，单位为秒，秒数与深度成正比，秒数越大，则深度越大；秒数越小，则深度越小。T轴对应的数值可以通过时深转换，转换成深度值。
从图4中的第一波阻抗模型、图7中的第二波阻抗模型至图10中的第三波阻抗模型，从依次得到的反演结果可以看出，得到的储层剖面包括的地层细节更加清晰，逐渐接近真实的储层剖面。
实施例五
为了执行上述方法类实施例及各个可能的方式中的相应步骤，下面给出一种波阻抗反演装置的实现方式，图11为本发明实施例提供的一种波阻抗反演装置的连接框图，如图11所示。需要说明的是，本实施例所提供的波阻抗反演装置，其基本原理及产生的技术效果和上述实施例相同，为简要描述，本实施例部分未提及之处，可参考上述的实施例中相应内容，波阻抗反演装置1包括：获取模块101、第一处理模块102、第二处理模块103、第三处理模块104以及第四处理模块105。
所述获取模块101，用于获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数；
可以理解地，所述获取模块101用于执行上述步骤S1。
所述第一处理模块102，用于根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型；
可以理解地，所述第一处理模块102用于执行上述步骤S2。
所述第二处理模块103，用于根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型；
可以理解地，所述第二处理模块103用于执行上述步骤S3。
所述第三处理模块104，用于根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型；
可以理解地，所述第三处理模块104用于执行上述步骤S4。
所述第四处理模块105，用于将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
可以理解地，所述第四处理模块105用于执行上述步骤S5。
可选地，获取模块101，具体用于根据地震资料获取所述中频波数；
采集测井内的所述低频波数以及所述高频波数。
可选地，所述第一处理模块102，具体用于将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果；
计算所述第一输出结果与第一测量结果之间的差值，得到第一残差；
判断所述第一残差是否大于第一预设残差；
若所述第一残差大于所述第一预设残差，则根据所述测量结果、所述第一输出结果以及第一对应关系，得到第一模型参数，其中，所述第一对应关系为所述测量结果、所述第一输出结果以及所述第一模型参数之间的对应关系；
根据所述第一模型参数调节所述初始波阻抗模型，得到更新后的初始波阻抗模型后，用所述更新后的初始波阻抗模型替代所述初始波阻已修改抗模型，执行所述将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果的步骤，直至所述第一残差不大于所述第一预设残差，得到第一波阻抗模型。
可选地，所述第一处理模块101，还用于将所述第一模型参数转换为第一三角函数；
将所述第一三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第一调节系数；
根据所述第一调节系数调节所述初始波阻抗模型。
可选地，所述第二处理模块，具体用于：
将所述中频波数输入所述第一波阻抗模型中，得到第二输出结果；
计算所述第二输出结果与第二测量结果之间的差值，得到第二残差；
判断所述第二残差是否大于第二预设残差；
若所述第二残差不大于所述第二预设残差，则根据所述测量结果、所述第二输出结果以及第二对应关系，得到第二模型参数，其中，所述第二对应关系为所述测量结果、所述第二输出结果以及所述第二模型参数之间的对应关系；
根据所述第二模型参数调节所述第一波阻抗模型，得到更新后的第一波阻抗模型后，用所述更新后的第一波阻抗模型替代所述第一波阻抗模型，执行所述将所述低频波数输入所述第一波阻抗模型中，得到第二输出结果的步骤，直至所述第二残差不大于所述第二预设残差，得到第二波阻抗模型。
可选地，所述第二处理模块103，还用于将所述第二模型参数转换为第二三角函数；
将所述第二三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第二调节系数；
根据所述第二调节系数调节所述第一波阻抗模型。
可选地，所述第三处理模块104，具体用于将所述高频波数输入所述第二波阻抗模型中，得到第三输出结果；
计算所述第三输出结果与第三测量结果之间的差值，得到第三残差；
判断所述第三残差是否大于第三预设残差；
若所述第三残差不大于所述第三预设残差，则根据所述测量结果、所述第三输出结果以及第单对应关系，得到第三模型参数，其中，所述第三对应关系为所述测量结果、所述第三输出结果以及所述第三模型参数之间的对应关系；
根据所述第三模型参数调节所述第二波阻抗模型，得到更新后的第二波阻抗模型后，用所述更新后的第二波阻抗模型替代所述第二波阻抗模型，执行所述将所述低频波数输入所述第二波阻抗模型中，得到第三输出结果的步骤，直至所述第三残差不大于所述第三预设残差，得到第三波阻抗模型。
可选地，所述第三处理模块104，还用于将所述第三模型参数转换为第三三角函数；
将所述第三三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第三调节系数；
根据所述第三调节系数调节所述第三波阻抗模型。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，获取模块101、第一处理模块102、第二处理模块103、第三处理模块104以及第四处理模块105的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
实施例六
本实施例提供一种存储介质，所述存储介质存储有计算机程序，所述存储介质被一个或多个处理器执行时，实现如实施例一至实施例四任意一实施例中所述的波阻抗反演方法。
上述存储介质可以是闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等。
实施例七
本实施例提供一种电子设备，所述电子设备包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，执行如实施例一至实施例四任意一实施例中所述的波阻抗反演方法。
处理器可以是专用集成电路(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)，磁存储器，快闪存储器，磁盘或光盘。
综上所述，本发明提供的一种波阻抗反演方法、装置、存储介质及电子设备，获取宽频带波数，其中，宽频带波数包括低频波数、中频波数以及高频波数；根据低频波数以及初始波阻抗模型，得到第一波阻抗模型；根据中频波数以及第一波阻抗模型，得到第二波阻抗模型；根据高频波数以及第二波阻抗模型，得到第三波阻抗模型，将采集到的地层信息输入第三波阻抗模型进行反演，得到波阻抗反演结果，其中，地层信息表征地层的岩性信息；依次经过低频波数、中频波数以及高频波数得到的第三波阻抗模型，第三波阻抗模型能够准确的反映储层的真实情况，将地层信息输入到第三波阻抗模型中进行反演，能够反演出分辨率较高的波阻抗反演结果，能够真实的反映储层情况。
在本发明所提供的实施例中，应该理解到，本发明所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，上述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
需要说明的是，在本发明中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
虽然本发明所揭露的实施方式如上，但上述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (11)

1.一种波阻抗反演方法，其特征在于，所述方法包括：
获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数；
根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型；
根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型；
根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型；
将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
2.根据权利要求1所述的波阻抗反演方法，其特征在于，所述获取宽频带波数的步骤，包括：
根据地震资料获取所述中频波数；
采集测井内的所述低频波数以及所述高频波数。
3.根据权利要求1所述的波阻抗反演方法，其特征在于，所述根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型的步骤，包括：
将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果；
计算所述第一输出结果与第一测量结果之间的差值，得到第一残差；
判断所述第一残差是否大于第一预设残差；
若所述第一残差大于所述第一预设残差，则根据所述第一测量结果、所述第一输出结果以及第一对应关系，得到第一模型参数，其中，所述第一对应关系为所述第一测量结果、所述第一输出结果以及所述第一模型参数之间的对应关系；
根据所述第一模型参数调节所述初始波阻抗模型，得到更新后的初始波阻抗模型后，用所述更新后的初始波阻抗模型替代所述初始波阻已修改抗模型，执行所述将所述低频波数输入所述初始波阻抗模型中，得到第一输出结果的步骤，直至所述第一残差不大于所述第一预设残差，得到第一波阻抗模型。
4.根据权利要求3所述的波阻抗反演方法，其特征在于，所述根据所述第一模型参数调节所述初始波阻抗模型的步骤，包括：
将所述第一模型参数转换为第一三角函数；
将所述第一三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第一调节系数；
根据所述第一调节系数调节所述初始波阻抗模型。
5.根据权利要求1所述的波阻抗反演方法，其特征在于，所述根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型的步骤，包括：
将所述中频波数输入所述第一波阻抗模型中，得到第二输出结果；
计算所述第二输出结果与第二测量结果之间的差值，得到第二残差；
判断所述第二残差是否大于第二预设残差；
若所述第二残差大于所述第二预设残差，则根据所述第二测量结果、所述第二输出结果以及第二对应关系，得到第二模型参数，其中，所述第二对应关系为所述第二测量结果、所述第二输出结果以及所述第二模型参数之间的对应关系；
根据所述第二模型参数调节所述第一波阻抗模型，得到更新后的第一波阻抗模型后，用所述更新后的第一波阻抗模型替代所述第一波阻抗模型，执行所述将所述低频波数输入所述第一波阻抗模型中，得到第二输出结果的步骤，直至所述第二残差不大于所述第二预设残差，得到第二波阻抗模型。
6.根据权利要求5所述的波阻抗反演方法，其特征在于，所述根据所述第二模型参数调节所述第一波阻抗模型的步骤，包括：
将所述第二模型参数转换为第二三角函数；
将所述第二三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第二调节系数；
根据所述第二调节系数调节所述第一波阻抗模型。
7.根据权利要求1所述的波阻抗反演方法，其特征在于，所述根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型的步骤，包括：
将所述高频波数输入所述第二波阻抗模型中，得到第三输出结果；
计算所述第三输出结果与第三测量结果之间的差值，得到第三残差；
判断所述第三残差是否大于第三预设残差；
若所述第三残差大于所述第三预设残差，则根据所述第三测量结果、所述第三输出结果以及第单对应关系，得到第三模型参数，其中，所述第三对应关系为所述第三测量结果、所述第三输出结果以及所述第三模型参数之间的对应关系；
根据所述第三模型参数调节所述第二波阻抗模型，得到更新后的第二波阻抗模型后，用所述更新后的第二波阻抗模型替代所述第二波阻抗模型，执行所述将所述低频波数输入所述第二波阻抗模型中，得到第三输出结果的步骤，直至所述第三残差不大于所述第三预设残差，得到第三波阻抗模型。
8.根据权利要求7所述的波阻抗反演方法，其特征在于，所述根据所述第三模型参数调节所述第二波阻抗模型的步骤，包括：
将所述第三模型参数转换为第三三角函数；
将所述第三三角函数转换为傅里叶级数展开式；
根据所述傅里叶级数展开式的系数，得到第三调节系数；
根据所述第三调节系数调节所述第三波阻抗模型。
9.一种波阻抗反演装置，其特征在于，包括：获取模块、第一处理模块、第二处理模块、第三处理模块和第四处理模块；
所述获取模块，用于获取宽频带波数，其中，所述宽频带波数包括低频波数、中频波数以及高频波数；
所述第一处理模块，用于根据所述低频波数以及初始波阻抗模型，得到第一波阻抗模型；
所述第二处理模块，用于根据所述中频波数以及所述第一波阻抗模型，得到第二波阻抗模型；
所述第三处理模块，用于根据所述高频波数以及所述第二波阻抗模型，得到第三波阻抗模型；
所述第四处理模块，用于将采集到的地层信息输入所述第三波阻抗模型进行反演，得到波阻抗反演结果，其中，所述地层信息表征地层的岩性信息。
10.一种存储介质，其特征在于，所述存储介质存储有计算机程序，所述存储介质被一个或多个处理器执行时，实现如权利要求1-8中任意一项所述的波阻抗反演方法。
11.一种电子设备，其特征在于，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，执行如权利要求1-8中任意一项所述的波阻抗反演方法。
CN202010646937.1A 2020-07-07 2020-07-07 波阻抗反演方法、装置、存储介质及电子设备 Active CN113917529B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010646937.1A CN113917529B (zh) 2020-07-07 2020-07-07 波阻抗反演方法、装置、存储介质及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010646937.1A CN113917529B (zh) 2020-07-07 2020-07-07 波阻抗反演方法、装置、存储介质及电子设备

Publications (2)

Publication Number Publication Date
CN113917529A true CN113917529A (zh) 2022-01-11
CN113917529B CN113917529B (zh) 2024-06-18

Family

ID=79231393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010646937.1A Active CN113917529B (zh) 2020-07-07 2020-07-07 波阻抗反演方法、装置、存储介质及电子设备

Country Status (1)

Country Link
CN (1) CN113917529B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4750157A (en) * 1987-05-06 1988-06-07 Standard Oil Production Company Seismic vibrator earth impedance determination and compensation system
CN104155687A (zh) * 2013-05-15 2014-11-19 中国石油化工股份有限公司 一种相控叠后声波阻抗反演方法
CN105182416A (zh) * 2015-09-06 2015-12-23 中国石油天然气股份有限公司 基于分频数据的地震反演方法及其装置
CN105425289A (zh) * 2015-10-29 2016-03-23 中国石油天然气集团公司 确定低频波阻抗的方法和装置
WO2017024702A1 (zh) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 一种射线弹性参数的反演系统
CN108919344A (zh) * 2018-03-30 2018-11-30 北京诺克斯达石油科技有限公司 适用于层状介质的分频构形反演方法
CN110967743A (zh) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 一种分频迭代地震反演方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4750157A (en) * 1987-05-06 1988-06-07 Standard Oil Production Company Seismic vibrator earth impedance determination and compensation system
CN104155687A (zh) * 2013-05-15 2014-11-19 中国石油化工股份有限公司 一种相控叠后声波阻抗反演方法
WO2017024702A1 (zh) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 一种射线弹性参数的反演系统
CN105182416A (zh) * 2015-09-06 2015-12-23 中国石油天然气股份有限公司 基于分频数据的地震反演方法及其装置
CN105425289A (zh) * 2015-10-29 2016-03-23 中国石油天然气集团公司 确定低频波阻抗的方法和装置
CN108919344A (zh) * 2018-03-30 2018-11-30 北京诺克斯达石油科技有限公司 适用于层状介质的分频构形反演方法
CN110967743A (zh) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 一种分频迭代地震反演方法及系统

Also Published As

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

Similar Documents

Publication Publication Date Title
Chen Dip-separated structural filtering using seislet transform and adaptive empirical mode decomposition based dip filter
Kalita et al. Regularized full-waveform inversion with automated salt flooding
Lorentzen et al. History matching the full Norne field model using seismic and production data
CN111123354B (zh) 基于频变反射振幅衰减预测致密气层的方法及设备
AU2019340410B2 (en) Reservoir characterization utilizing inversion of resampled seismic data
CN109001840B (zh) 一种泥页岩成因分析方法及装置
CN101201409B (zh) 一种地震数据变相位校正方法
CN104268412A (zh) 一种角道集射线层析偏移速度分析方法及装置
CN106772586A (zh) 一种基于地震信号奇异性的隐蔽性断裂检测方法
Wang et al. Sample entropy‐based adaptive wavelet de‐noising approach for meteorologic and hydrologic time series
Su et al. Detection of strong precession cycles from the late Pliocene sedimentary records of northeastern Tibetan Plateau
CN112213782B (zh) 分相位地震数据的处理方法、装置和服务器
CN107942374A (zh) 绕射波场提取方法和装置
CN105044773A (zh) 转换波井震标定方法及装置
CN113917529A (zh) 波阻抗反演方法、装置、存储介质及电子设备
CN112649887B (zh) 基于钻井资料定量划分层序地层的方法及装置
CN111830562B (zh) 一种油气储层渗透率预测方法及装置
CN114200522A (zh) 深度域地震子波提取方法和装置、存储介质及电子设备
Heidary et al. A wavelet-based model for determining asphaltene onset pressure
Wang Heterogeneous Seismic Waves Pattern Recognition in Oil Exploration with Spectrum Imaging
CN109782343A (zh) 地层旋回分析方法及装置
Wu et al. Sparse prior-net: A sparse prior-based deep network for seismic data interpolation
CN117310802B (zh) 一种深度域储层地震反演方法
CN212535642U (zh) 一种校正温度漂移的装置和设备
CN117849875B (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