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

    公开(公告)号: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)。本发明可快速完成速度模型更新,广泛应用于地震波全波形反演方法中。

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

    公开(公告)号:CN104898167A

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

    申请号:CN201510340938.2

    申请日:2015-06-18

    IPC分类号: G01V1/50

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

    基于多次聚焦的噪音压制与成像方法

    公开(公告)号:CN102636810B

    公开(公告)日:2014-05-14

    申请号:CN201210097757.8

    申请日:2012-04-05

    IPC分类号: G01V1/28 G01V1/36

    摘要: 本发明涉及一种基于多次聚焦的噪音压制与成像方法,其特征在于:包括如下步骤:1)利用处理系统提供的分选功能将预处理之后的地震数据分选为共偏移距道集;2)基于共偏移距道集实施多次聚焦属性参数搜索,得到能够描述地震波局部特征的属性参数;3)基于上述属性参数构造多次聚焦局部反射线段;4)基于多次聚焦局部反射线段构造等旅行时面的多次聚焦算子响应;5)共偏移距剖面内所有的多次聚焦算子响应累加构成多次聚焦处理后数据体;6)将多次聚焦处理后数据体进行共中心点叠加,得到多次聚焦零偏移距成像剖面。本发明不但可以提供高质量的零偏移距成像剖面,通过充分利用菲涅耳体内的地震波能量,还可以输出信噪比大幅度提高的叠前道集。

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

    公开(公告)号: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)。本发明可快速完成速度模型更新,广泛应用于地震波全波形反演方法中。

    基于多次聚焦的噪音压制与成像方法

    公开(公告)号:CN102636810A

    公开(公告)日:2012-08-15

    申请号:CN201210097757.8

    申请日:2012-04-05

    IPC分类号: G01V1/28 G01V1/36

    摘要: 本发明涉及一种基于多次聚焦的噪音压制与成像方法,其特征在于:包括如下步骤:1)利用处理系统提供的分选功能将预处理之后的地震数据分选为共偏移距道集;2)基于共偏移距道集实施多次聚焦属性参数搜索,得到能够描述地震波局部特征的属性参数;3)基于上述属性参数构造多次聚焦局部反射线段;4)基于多次聚焦局部反射线段构造等旅行时面的多次聚焦算子响应;5)共偏移距剖面内所有的多次聚焦算子响应累加构成多次聚焦处理后数据体;6)将多次聚焦处理后数据体进行共中心点叠加,得到多次聚焦零偏移距成像剖面。本发明不但可以提供高质量的零偏移距成像剖面,通过充分利用菲涅耳体内的地震波能量,还可以输出信噪比大幅度提高的叠前道集。

    基于匹配地震子波的物理小波的地震瞬时频率分析方法

    公开(公告)号:CN102353991B

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

    申请号:CN201110154420.1

    申请日:2011-06-09

    IPC分类号: G01V1/30

    摘要: 本发明涉及一种基于匹配地震子波的物理小波的地震瞬时频率分析方法,包括如下步骤:1)获取二维或三维的经偏移或叠加处理后的地震资料;2)根据研究对象对所获取的地震资料进行空间分区,区域里获得测井资料或零偏VSP资料、井旁地震记录、储层的地质构造及其它先验信息;3)通过测井资料或零偏VSP资料,以及井旁地震记录反演地震子波,确定匹配该子波的母物理小波;4)在物理小波域计算地震信号对应的解析信号;5)根据得到的解析信号,基于极平坦滤波器计算瞬时频率;6)根据得到的瞬时频率,进行最佳分辨率瞬时频率分析。本发明具有多分辨率特性,适用于低信噪比资料,对宽频带地震资料可得到高精度瞬时频率。

    基于匹配地震子波的物理小波的地震瞬时频率分析方法

    公开(公告)号:CN102353991A

    公开(公告)日:2012-02-15

    申请号:CN201110154420.1

    申请日:2011-06-09

    IPC分类号: G01V1/30

    摘要: 本发明涉及一种基于匹配地震子波的物理小波的地震瞬时频率分析方法,包括如下步骤:1)获取二维或三维的经偏移或叠加处理后的地震资料;2)根据研究对象对所获取的地震资料进行空间分区,区域里获得测井资料或零偏VSP资料、井旁地震记录、储层的地质构造及其它先验信息;3)通过测井资料或零偏VSP资料,以及井旁地震记录反演地震子波,确定匹配该子波的母物理小波;4)在物理小波域计算地震信号对应的解析信号;5)根据得到的解析信号,基于极平坦滤波器计算瞬时频率;6)根据得到的瞬时频率,进行最佳分辨率瞬时频率分析。本发明具有多分辨率特性,适用于低信噪比资料,对宽频带地震资料可得到高精度瞬时频率。