深水层间多次波消除方法

    公开(公告)号:CN102169189B

    公开(公告)日:2015-01-14

    申请号:CN201110021675.0

    申请日:2011-01-19

    IPC分类号: G01V1/36

    摘要: 本发明涉及深水层间多次波消除方法,应用于油气勘探地震资料处理中,它包括1、消除自由表面多次波;2、对原始数据抽道集;3、将抽道集得到的数据(xs,xg;t)进行二维傅立叶变换得到频率-波数域数据(ks,kg;ω);4、利用逆散射级数和频率-波数域数据D(ks,kg;ω)推导出n阶散射波场Dn;5、根据Dn得到消除与自由表面相关多次波后的波场D0;6、依据散射波场形成层间多次波的两个约束条件,及下行散射点的深度Zg比相邻的两个散射点Zg,和Zg,迭代计算出Vn,得到有关层间多次波的级数项;7、将逆散射级数项进行反傅立叶变换,得到d1(xs,xg;t);8、将d(xs,xg;t)与d1(xs,xg;t)数据匹配相减,得到消除层间多次波以后的CMP道集。本发明基于量子力学的散射理论,提出深水层间多次波消除技术,填补了地震资料深水层间多次波消除的空白。

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

    公开(公告)号:CN105005076B

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

    申请号:CN201510295941.7

    申请日:2015-06-02

    IPC分类号: G01V1/30

    摘要: 本发明涉及一种基于最小二乘梯度更新速度模型的地震波全波形反演方法,包括以下步骤: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

    IPC分类号: G01V1/30

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

    深水层间多次波消除方法

    公开(公告)号:CN102169189A

    公开(公告)日:2011-08-31

    申请号:CN201110021675.0

    申请日:2011-01-19

    IPC分类号: G01V1/36

    摘要: 本发明涉及深水层间多次波消除方法,应用于油气勘探地震资料处理中,它包括1、消除自由表面多次波;2、对原始数据抽道集;3、将抽道集得到的数据(xs,xg;t)进行二维傅立叶变换得到频率-波数域数据(ks,kg;ω);4、利用逆散射级数和频率-波数域数据D(ks,kg;ω)推导出n阶散射波场Dn;5、根据Dn得到消除与自由表面相关多次波后的波场D0;6、依据散射波场形成层间多次波的两个约束条件,及下行散射点的深度Zg比相邻的两个散射点Zg,和Zg,,小,迭代计算出Vn,得到有关层间多次波的级数项;7、将逆散射级数项进行反傅立叶变换,得到d1(xs,xg;t);8、将d(xs,xg;t)与d1(xs,xg;t)数据匹配相减,得到消除层间多次波以后的CMP道集。本发明基于量子力学的散射理论,提出深水层间多次波消除技术,填补了地震资料深水层间多次波消除的空白。