基于局部自适应凸化方法的多尺度地震全波形反演方法

    公开(公告)号:CN107422379A

    公开(公告)日:2017-12-01

    申请号:CN201710623459.0

    申请日:2017-07-27

    IPC分类号: G01V1/30

    CPC分类号: G01V1/303 G01V2210/6222

    摘要: 本发明涉及一种基于局部自适应凸化方法的多尺度地震全波形反演方法,其步骤:预处理,以零值序列作为初始值,初始速度模型作为起始值;截取直达波信息;得到正演模拟直达波;建立震源函数反演的目标函数;得到直达波残差、直达波残差反传波场;计算震源函数的更新梯度、更新方向并寻找步长;输出反演得到的高精度震源函数;进行衰减时窗处理;模拟时窗内的地震数据;进行局部凸化处理和分离处理;去除观测记录中的高频成分;建立最小二乘目标函数,得到波场残差;得到模型空间的残差反传波场;得到模型更新梯度;计算模型更新方向,并寻找步长;输出多尺度地震全波形反演的反演结果;输出最终反演结果。本发明在地震勘探技术领域中广泛应用。

    一种频率域全波形反演地震速度建模方法

    公开(公告)号:CN103207409A

    公开(公告)日:2013-07-17

    申请号:CN201310133893.2

    申请日:2013-04-17

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种频率域全波形反演地震速度建模方法,包括以下步骤:1)获取原始地震炮集记录、震源子波信息和反演所用的初始模型;2)分析步骤1)获得的信息,基于正演算法和优化算法,确定基本的反演参数和从低频到高频的全波形反演框架;3)对不同的频率,计算获得最合适的频域正反演模型网格;4)在低频反演时采用主成分分析方法,压缩参与反演的数据维度;5)判断对应不同频率的投影矩阵维度是否满足阈值转换标准,当满足该转换标准时,进入下一步,否则回到步骤4);6)引入震源编码方法,并利用随机相位编码方法压制串扰噪声;7)判定是否满足迭代截止条件,如果满足迭代截止条件,则进入下一步;否则回到步骤6);8)如果没有完成所有频率的反演,则回到步骤3),直至完成所有频率的反演,获得最后的速度模型,输出速度模型。

    基于小生境主从式并行遗传算法的双相介质参数反演方法

    公开(公告)号:CN103235879A

    公开(公告)日:2013-08-07

    申请号:CN201310134062.7

    申请日:2013-04-17

    IPC分类号: G06F19/00 G06N3/12

    摘要: 本发明涉及一种基于小生境主从式并行遗传算法的双相各向异性介质储层参数反演方法,是利用小生境主从式并行遗传算法求解双相各向异性介质储层参数,其核心思想是系统分为一个主处理器和若干从处理器,主处理器监控整个种群,在适应度的计算阶段,由主处理器将适应度的计算分配到各个从处理器上进行,计算完成之后再由主处理器收集结果,然后由主处理器进行小生境淘汰、选择、交叉、变异等操作,并由此产生新一代种群,从而完成一次循环,极大地提高了储层参数反演的计算效率。本发明在求解储层参数进化过程中引入共享度的概念,通过调整各个个体的适应度来限制个别个体的大量增长,造就小生境的进化环境,提高了遗传算法求解多峰值储层参数反演优化问题的能力和求解质量。本发明可以广泛用于油气藏参数反演过程中。

    一种频率域全波形反演地震速度建模方法

    公开(公告)号:CN103207409B

    公开(公告)日:2016-01-06

    申请号:CN201310133893.2

    申请日:2013-04-17

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种频率域全波形反演地震速度建模方法,包括以下步骤:1)获取原始地震炮集记录、震源子波信息和反演所用的初始模型;2)分析步骤1)获得的信息,基于正演算法和优化算法,确定基本的反演参数和从低频到高频的全波形反演框架;3)对不同的频率,计算获得最合适的频域正反演模型网格;4)在低频反演时采用主成分分析方法,压缩参与反演的数据维度;5)判断对应不同频率的投影矩阵维度是否满足阈值转换标准,当满足该转换标准时,进入下一步,否则回到步骤4);6)引入震源编码方法,并利用随机相位编码方法压制串扰噪声;7)判定是否满足迭代截止条件,如果满足迭代截止条件,则进入下一步;否则回到步骤6);8)如果没有完成所有频率的反演,则回到步骤3),直至完成所有频率的反演,获得最后的速度模型,输出速度模型。

    一种基于复地震道分解和重构的地震相位和频率校正方法

    公开(公告)号:CN103852788A

    公开(公告)日:2014-06-11

    申请号:CN201410068809.8

    申请日:2014-02-27

    IPC分类号: G01V1/28 G01V1/36

    摘要: 本发明涉及一种基于复地震道分解和重构的地震相位和频率校正方法,其包括以下步骤:1)根据已有地震记录的有效频带范围,建立遍历有效频带范围的子波库,及非平稳褶积模型;2)将子波库通过Hilbert变换变为复数,建立复子波库;3)利用复子波库将非平稳褶积模型改写成复数形式,进而通过线性运算及L1范数约束的稀疏反演求解复反射系数4)将复反射系数代入复数形式的非平稳褶积模型中,并取其实数部分,实现对地震剖面的重构,进而通过在重构前改变复子波库中地震子波的频率及求得的复反射系数中的相位信息,实现对原地震记录进行频率、相位以频率和相位同时校正的目的。本发明可以广泛在石油勘探过程中应用。

    基于局部自适应凸化方法的多尺度地震全波形反演方法

    公开(公告)号:CN107422379B

    公开(公告)日:2019-01-11

    申请号:CN201710623459.0

    申请日:2017-07-27

    IPC分类号: G01V1/30

    摘要: 本发明涉及一种基于局部自适应凸化方法的多尺度地震全波形反演方法,其步骤:预处理,以零值序列作为初始值,初始速度模型作为起始值;截取直达波信息;得到正演模拟直达波;建立震源函数反演的目标函数;得到直达波残差、直达波残差反传波场;计算震源函数的更新梯度、更新方向并寻找步长;输出反演得到的高精度震源函数;进行衰减时窗处理;模拟时窗内的地震数据;进行局部凸化处理和分离处理;去除观测记录中的高频成分;建立最小二乘目标函数,得到波场残差;得到模型空间的残差反传波场;得到模型更新梯度;计算模型更新方向,并寻找步长;输出多尺度地震全波形反演的反演结果;输出最终反演结果。本发明在地震勘探技术领域中广泛应用。