一种叠前地震贝叶斯反演方法和装置制造方法及图纸

本发明专利技术公开了一种叠前地震贝叶斯反演方法和装置。该方法利用地震数据中的时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程。进一步，以贝叶斯反演模型为基础，利用第一正演方程和第二正演方程，建立时频域联合叠前地震反演目标函数，预先设定的参数约束条件的约束下，确定时频域联合叠前地震反演目标函数的最优解，得到地层弹性参数。本发明专利技术充分利用叠前地震数据中的时间域及频率域数据，使得反演结果的分辨率比现有方法有较大提高。

【技术实现步骤摘要】

本专利技术涉及地震监测领域，更具体的说是涉及一种叠前地震贝叶斯反演方法和装 置。
技术介绍
地震反演是获取地下介质弹性参数的有效途径。根据所采用地震资料的不同，地 震反演可分为叠后反演和叠前反演。 在叠前地震反演中，贝叶斯反演在考虑模型及数据参数不完备性方面更有优势。 贝叶斯反演理论在地球物理反演中应用广泛，该方法分别假设模型参数和似然函数服从某 种先验分布，通过求解后验概率密度函数实现反演估计。 现有的叠前地震贝叶斯反演的实现主要在时间域进行，反演得到地层介弹性参 数，其反演结果的分辨率低。
技术实现思路
有鉴于此，本专利技术提供一种叠前地震贝叶斯反演方法和装置，采用时频域数据联 合反演的方式，得到地层弹性参数，提高了反演结果的分辨率。 为实现上述目的，本专利技术提供如下技术方案： 一种叠前地震贝叶斯反演方法，包括： 获取时间域叠前地震数据和频率域叠前地震数据； 分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒， 得到第一正演方程和第二正演方程； 以贝叶斯反演模型为基础，利用所述第一正演方程和所述第二正演方程，建立时 频域联合叠前地震反演目标函数； 基于预先设定的参数约束条件，确定所述时频域联合叠前地震反演目标函数的最 优解，得到地层弹性参数。 优选的，所述分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行 正演算推倒，得到第一正演方程和第二正演方程，之后还包括： 利用预设算法，对所述第一正演方程和所述第二正演方程中的参数进行去相关处 理。 优选的，所述获取时间域叠前地震数据和频率域叠前地震数据，包括： 获取时间域叠前地震数据； 对所述时间域叠前地震数据进行傅里叶变换，得到所述频率域叠前地震数据。 优选的，所述时间域叠前地震数据为至少三个角叠加地震数据。 一种叠前地震贝叶斯反演装置，包括： 数据采集单元，用于获取时间域叠前地震数据和频率域叠前地震数据； 正演算单元，用于分别利用所述时间域叠前地震数据和所述频率域叠前地震数据 进行正演算推倒，得到第一正演方程和第二正演方程； 反演目标函数建立单元，用于以贝叶斯反演模型为基础，利用所述第一正演方程 和所述第二正演方程，建立时频域联合叠前地震反演目标函数； 地层弹性参数计算单元，用于基于预先设定的参数约束条件，确定所述时频域联 合叠前地震反演目标函数的最优解，得到地层弹性参数。 优选的，还包括：连接与所述正演单元和所述反演目标函数建立单元之间的参数 去相关处理单元，用于利用预设算法，对所述第一正演方程和所述第二正演方程中的参数 进行去相关处理。 优选的，所述数据采集单元包括： 数据采集子单元，用于获取时间域叠前地震数据； 傅里叶变换单元，用于对所述时间域叠前地震数据进行傅里叶变换，得到所述频 率域叠前地震数据。 经由上述的技术方案可知，与现有技术相比，本专利技术公开了一种叠前地震贝叶斯 反演方法和装置。该方法利用地震数据中的时间域叠前地震数据和所述频率域叠前地震数 据进行正演算推倒，得到第一正演方程和第二正演方程。进一步，以贝叶斯反演模型为基 础，利用第一正演方程和第二正演方程，建立时频域联合叠前地震反演目标函数，预先设定 的参数约束条件的约束下，确定时频域联合叠前地震反演目标函数的最优解，得到地层弹 性参数。本专利技术充分利用叠前地震数据中的时间域及频率域数据，使得反演结果的分辨率 比现有方法有$父大提尚。【附图说明】 为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 专利技术的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 提供的附图获得其他的附图。 图1示出了本专利技术一实施例公开的一种地震贝叶斯反演方法的流程示意图； 图2示出了本专利技术另一实施例公开的一种地震贝叶斯反演方法的流程示意图； 图3示出了本专利技术一实施例公开的一种地震贝叶斯反演装置的结构示意图。【具体实施方式】 下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于 本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本专利技术保护的范围。 参见图1示出了本专利技术一实施例公开的一种地震贝叶斯反演方法的流程示意图。 由图1可知，该方法包括： 101 :获取时间域叠前地震数据和频率域叠前地震数据。 可选的，在本实施例中对通过所述时间域叠前地震数据进行傅里叶变换，得到所 述频率域叠前地震数据。 102:分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算 推倒，得到第一正演方程和第二正演方程。 需要说明的是，该正演算推倒为现有技术手段，其具体实现方式如下： 地震反射系数Rpp与地层弹性参数（纵波速度Vp、横波速度Vs及密度P )之间的 关系为：【主权项】1. 一种叠前地震贝叶斯反演方法，其特征在于，包括： 获取时间域叠前地震数据和频率域叠前地震数据； 分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到 第一正演方程和第二正演方程； 以贝叶斯反演模型为基础，利用所述第一正演方程和所述第二正演方程，建立时频域 联合叠前地震反演目标函数； 基于预先设定的参数约束条件，确定所述时频域联合叠前地震反演目标函数的最优 解，得到地层弹性参数。2. 根据权利要求1所述的方法，其特征在于，所述分别利用所述时间域叠前地震数据 和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程，之后还 包括： 利用预设算法，对所述第一正演方程和所述第二正演方程中的参数进行去相关处理。3. 根据权利要求1所述的方法，其特征在于，所述获取时间域叠前地震数据和频率域 叠前地震数据，包括： 获取时间域叠前地震数据； 对所述时间域叠前地震数据进行傅里叶变换，得到所述频率域叠前地震数据。4. 根据权利要求3所述的方法，其特征在于，所述时间域叠前地震数据为至少三个角 叠加地震数据。5. -种叠前地震贝叶斯反演装置，其特征在于，包括： 数据采集单元，用于获取时间域叠前地震数据和频率域叠前地震数据； 正演算单元，用于分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行 正演算推倒，得到第一正演方程和第二正演方程； 反演目标函数建立单元，用于以贝叶斯反演模型为基础，利用所述第一正演方程和所 述第二正演方程，建立时频域联合叠前地震反演目标函数； 地层弹性参数计算单元，用于基于预先设定的参数约束条件，确定所述时频域联合叠 前地震反演目标函数的最优解，得到地层弹性参数。6. 根据权利要求5所述的装置，其特征在于，还包括：连接与所述正演单元和所述反演 目标函数建立单元之间的参数去相关处理单元，用于利用预设算法，对所述第一正演方程 和所述第二正演方程中的参数进行去相关处理。7. 根据权利要求5所述的装置，其特征在于，所述数据采集单元包括： 数据采集子单元，用于获取时间域叠前地震数据； 傅里叶变换单元，用于对所述时间域叠前地震数据进行傅里叶变换，得到所述频率域 叠前地震数据。【专利摘要】本专利技术公开了本文档来自技高网...

【技术保护点】
一种叠前地震贝叶斯反演方法，其特征在于，包括：获取时间域叠前地震数据和频率域叠前地震数据；分别利用所述时间域叠前地震数据和所述频率域叠前地震数据进行正演算推倒，得到第一正演方程和第二正演方程；以贝叶斯反演模型为基础，利用所述第一正演方程和所述第二正演方程，建立时频域联合叠前地震反演目标函数；基于预先设定的参数约束条件，确定所述时频域联合叠前地震反演目标函数的最优解，得到地层弹性参数。

【技术特征摘要】

【专利技术属性】
技术研发人员：宗兆云，印兴耀，吴国忱，李坤，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37