弹性介质波动方程反射波旅行时反演方法

    公开(公告)号:CN108873066A

    公开(公告)日:2018-11-23

    申请号:CN201810666471.4

    申请日:2018-06-26

    IPC分类号: G01V1/28

    摘要: 本发明公开了一种弹性介质波动方程反射波旅行时反演方法。设计新的旅行时残差目标函数;推导新目标函数下的伴随方程和背景参数梯度公式;计算背景参数的梯度;采用拟牛顿法反演算法对梯度进行处理;采用抛物线拟合法求取迭代步长;更新背景参数模型,直到满足收敛条件。本发明的有益效果是通过采用新的旅行时目标函数来提高弹性介质反射波波形反演的精度,获取可靠的纵波速度和横波速度初始模型。

    变密度声波方程时间高阶精度交错网格有限差分方法

    公开(公告)号:CN108279437A

    公开(公告)日:2018-07-13

    申请号:CN201810042714.7

    申请日:2018-01-17

    IPC分类号: G01V1/30 G01V1/36

    摘要: 本发明公开了变密度声波方程时间高阶精度交错网格有限差分方法。首先设计新的差分模板;采用新的差分模板对变密度声波方程进行交错网格有限差分数值离散;推导新差分模板下变密度声波方程的频散关系;采用泰勒级数展开或最优化方法求取新的有限差分系数;将新的差分系数打入相应的离散方程,实现变密度声波方程数值模拟;分析新差分方法的精度和稳定性。本发明的有益效果是具有精确的空间(2M)阶、时间10阶模拟精度。可通过增大时间步长来提高计算效率。在远距离模拟中能更好地保证波形不发生畸变。

    一种基于动校正速度参数的地震资料各向异性建模方法

    公开(公告)号:CN109901221B

    公开(公告)日:2021-03-19

    申请号:CN201910165409.1

    申请日:2019-03-05

    发明人: 任志明 李振春

    IPC分类号: G01V1/28 G01V1/30 G01V1/36

    摘要: 本发明提供了一种基于动校正速度参数的地震资料各向异性建模方法,包括:读取地震观测记录,构建地震观测系统;参数初始模型建立,并转换为新的参数化方式(vn,δ,η,ρ);使用新的参数方式(vn,δ,η,ρ)进行震源波场正向延拓,求取记录残差;使用新的参数方式(vn,δ,η,ρ)震源波场逆向延拓,求取梯度;求取迭代步长,更新vn,δ参数;判断是否满足收敛条件,输出模型。本发明基于敏感性分析采用了新的参数化方式,提高了δ参数的敏感性,有利于恢复地下各向异性参数,改善复杂地质构造的成像和解释结果。

    一种回折波和反射波波动方程旅行时反演方法

    公开(公告)号:CN110187382B

    公开(公告)日:2020-10-13

    申请号:CN201910164911.0

    申请日:2019-03-05

    发明人: 任志明 李振春

    IPC分类号: G01V1/28 G01V1/30

    摘要: 本发明公开了一种回折波和反射波波动方程旅行时反演方法,包括地震资料回折波和反射波分离;回折波和反射波旅行时残差求取;旅行时残差目标函数构建;背景速度与反射系数梯度计算;反射系数模型更新迭代;背景速度模型更新迭代等技术环节。本发明通过采用回折波和反射波进行波动方程旅行时反演来获取可靠的地下速度初始模型,进而提高后续全波形反演和逆时偏移的反演和成像精度。

    基于正则化约束的弹性波最小二乘逆时偏移方法

    公开(公告)号:CN108333628A

    公开(公告)日:2018-07-27

    申请号:CN201810042691.X

    申请日:2018-01-17

    IPC分类号: G01V1/30

    摘要: 本发明公开了基于正则化约束的弹性波最小二乘逆时偏移方法。设计新的目标函数;推导新目标函数下的弹性波反偏移算子和反射系数梯度公式;计算反射系数的梯度;采用共轭梯度法或拟牛顿法反演算法对梯度进行处理;采用抛物线拟合法求取迭代步长;更新反射系数模型,直到满足收敛条件。本发明的有益效果是通过采用新的全变分正则化约束策略来提高成像分辨率及稳定性。

    基于正则化约束的弹性波最小二乘逆时偏移方法

    公开(公告)号:CN108333628B

    公开(公告)日:2019-09-03

    申请号:CN201810042691.X

    申请日:2018-01-17

    IPC分类号: G01V1/30

    摘要: 本发明公开了基于正则化约束的弹性波最小二乘逆时偏移方法。设计新的目标函数;推导新目标函数下的弹性波反偏移算子和反射系数梯度公式;计算反射系数的梯度;采用共轭梯度法或拟牛顿法反演算法对梯度进行处理;采用抛物线拟合法求取迭代步长;更新反射系数模型,直到满足收敛条件。本发明的有益效果是通过采用新的全变分正则化约束策略来提高成像分辨率及稳定性。

    一种VSP地震资料转换波全波形反演方法

    公开(公告)号:CN108845351A

    公开(公告)日:2018-11-20

    申请号:CN201810666468.2

    申请日:2018-06-26

    IPC分类号: G01V1/28 G01V1/30

    CPC分类号: G01V1/282 G01V1/30

    摘要: 本发明公开了一种VSP地震资料转换波全波形反演方法。通过全局优化参数反演法分离上行PS转换波;根据经验公式计算横波速度并对其进行平滑处理得到初始横波速度模型;通过求解纵横波分离的弹性波方程得到每一时间步的正传纵横波波场以及模拟VSP PS波记录;通过求解纵横波分离的伴随弹性波方程得到每一时间步的反传纵横波波场;采用正传纵波波场和反传横波波场计算横波速度梯度;采用三点抛物线拟合法求取迭代步长;通过L-BFGS方法更新横波速度模型;运用维纳滤波从低频到高频进行多尺度反演。本发明的有益效果是通过纵横波分离策略来提高VSP地震资料转换波全波形反演的精度,提高复杂构造的反演效果。

    弹性介质波动方程反射波旅行时反演方法

    公开(公告)号:CN108873066B

    公开(公告)日:2020-03-10

    申请号:CN201810666471.4

    申请日:2018-06-26

    IPC分类号: G01V1/28

    摘要: 本发明公开了一种弹性介质波动方程反射波旅行时反演方法。设计新的旅行时残差目标函数;推导新目标函数下的伴随方程和背景参数梯度公式;计算背景参数的梯度;采用拟牛顿法反演算法对梯度进行处理;采用抛物线拟合法求取迭代步长;更新背景参数模型,直到满足收敛条件。本发明的有益效果是通过采用新的旅行时目标函数来提高弹性介质反射波波形反演的精度,获取可靠的纵波速度和横波速度初始模型。

    一种回折波和反射波波动方程旅行时反演方法

    公开(公告)号:CN110187382A

    公开(公告)日:2019-08-30

    申请号:CN201910164911.0

    申请日:2019-03-05

    发明人: 任志明 李振春

    IPC分类号: G01V1/28 G01V1/30

    摘要: 本发明公开了一种回折波和反射波波动方程旅行时反演方法,包括地震资料回折波和反射波分离;回折波和反射波旅行时残差求取;旅行时残差目标函数构建;背景速度与反射系数梯度计算;反射系数模型更新迭代;背景速度模型更新迭代等技术环节。本发明通过采用回折波和反射波进行波动方程旅行时反演来获取可靠的地下速度初始模型,进而提高后续全波形反演和逆时偏移的反演和成像精度。

    一种基于动校正速度参数的地震资料各向异性建模方法

    公开(公告)号:CN109901221A

    公开(公告)日:2019-06-18

    申请号:CN201910165409.1

    申请日:2019-03-05

    发明人: 任志明 李振春

    IPC分类号: G01V1/28 G01V1/30 G01V1/36

    摘要: 本发明提供了一种基于动校正速度参数的地震资料各向异性建模方法,包括:读取地震观测记录,构建地震观测系统;参数初始模型建立,并转换为新的参数化方式(vn,δ,η,ρ);使用新的参数方式(vn,δ,η,ρ)进行震源波场正向延拓,求取记录残差;使用新的参数方式(vn,δ,η,ρ)震源波场逆向延拓,求取梯度;求取迭代步长,更新vn,δ参数;判断是否满足收敛条件,输出模型。本发明基于敏感性分析采用了新的参数化方式,提高了δ参数的敏感性,有利于恢复地下各向异性参数,改善复杂地质构造的成像和解释结果。