基于最小二乘梯度更新速度模型的地震波全波形反演方法

    公开(公告)号:CN105005076B

    公开(公告)日:2017-05-03

    申请号:CN201510295941.7

    申请日:2015-06-02

    Abstract: 本发明涉及一种基于最小二乘梯度更新速度模型的地震波全波形反演方法,包括以下步骤:1)获得时间域内地震波场的声波方程;2)构建初始速度模型,设置速度模型更新迭代次数N与允许最小误差值ε;3)构造观测波场数据和计算波场数据的波场误差矢量;4)构造目标函数;5)对目标函数计算获得地震波全波形反演的标准方程;6)引入速度模型更新梯度方向gk和更新量步长α;7)利用最小二乘方法求解速度模型更新梯度方向gk;8)对更新量步长α进行插值计算;9)对速度模型进行更新得:mk=mk‑1+αgk;当|αgk|<ε时或速度模型更新次数达到速度模型更新迭代次数N时,速度模型更新结束;否则,进入步骤3)。本发明可快速完成速度模型更新,广泛应用于地震波全波形反演方法中。

    基于最小二乘梯度更新速度模型的地震波全波形反演方法

    公开(公告)号:CN105005076A

    公开(公告)日:2015-10-28

    申请号:CN201510295941.7

    申请日:2015-06-02

    Abstract: 本发明涉及一种基于最小二乘梯度更新速度模型的地震波全波形反演方法,包括以下步骤:1)获得时间域内地震波场的声波方程;2)构建初始速度模型,设置速度模型更新迭代次数N与允许最小误差值ε;3)构造观测波场数据和计算波场数据的波场误差矢量;4)构造目标函数;5)对目标函数计算获得地震波全波形反演的标准方程;6)引入速度模型更新梯度方向gk和更新量步长α;7)利用最小二乘方法求解速度模型更新梯度方向gk;8)对更新量步长α进行插值计算;9)对速度模型进行更新得:mk=mk-1+αgk;当|αgk|<ε时或速度模型更新次数达到速度模型更新迭代次数N时,速度模型更新结束;否则,进入步骤3)。本发明可快速完成速度模型更新,广泛应用于地震波全波形反演方法中。

    一种稀疏正则化的岩石物理弹性参数提取方法

    公开(公告)号:CN104898167A

    公开(公告)日:2015-09-09

    申请号:CN201510340938.2

    申请日:2015-06-18

    Abstract: 本发明涉及一种稀疏正则化的岩石物理弹性参数提取方法,包括以下步骤:1)叠前地震道集分小、中、大三个角度范围进行叠加,获得小、中两个角度对应的弹性阻抗;2)将小、中两个角度及对应的弹性阻抗代入弹性阻抗方程得到弹性阻抗方程组;3)令则弹性阻抗方程组为欠定方程组:Ax=b;4)对欠定方程组进行运算获得待求解线性方程组:5)对线性方程组求解得到(xk,zk,sk)T;6)根据计算获得弹性参数纵波速度Vp、横波速度Vs和密度ρ。本发明是提供一种只利用小、中两个角度的地震道集数据获得的弹性阻抗数据就可以准确提取弹性参数的稀疏正则化的弹性参数提取方法,可以广泛应用于岩石物理弹性参数的提取方法中。

Patent Agency Ranking