-
公开(公告)号: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
摘要: 本发明涉及一种基于多次聚焦的噪音压制与成像方法,其特征在于:包括如下步骤: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)。本发明可快速完成速度模型更新,广泛应用于地震波全波形反演方法中。
-
公开(公告)号:CN104898167B
公开(公告)日:2017-03-01
申请号:CN201510340938.2
申请日:2015-06-18
IPC分类号: G01V1/50
摘要: 方法,可以广泛应用于岩石物理弹性参数的提取本发明涉及一种稀疏正则化的岩石物理弹 方法中。性参数提取方法,包括以下步骤:1)叠前地震道集分小、中、大三个角度范围进行叠加,获得小、中两个角度对应的弹性阻抗;2)将小、中两个角度及对应的弹性阻抗代入弹性阻抗方程得到弹性阻抗方程组;3)令则弹性阻抗方程组为欠定方程组:Ax=b;4)对欠定方程组进行运算获得待求解线性方程组: 5)对线性方程组求解得到(xk,zk,sk)T;6)根据计算获得弹性参数纵波速度Vp、横波速度Vs和密(56)对比文件张世鑫等.泌阳凹陷安店区块薄储层精细预测《.石油物探》.2010,第49卷(第1期),David N. Whitcombe et al..Extendedelastic impedance for fluid and lithologyprediction《.GEOPHYSICS》.2002,第67卷(第1期),
-
公开(公告)号:CN102636810A
公开(公告)日:2012-08-15
申请号:CN201210097757.8
申请日:2012-04-05
摘要: 本发明涉及一种基于多次聚焦的噪音压制与成像方法,其特征在于:包括如下步骤: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)根据得到的瞬时频率,进行最佳分辨率瞬时频率分析。本发明具有多分辨率特性,适用于低信噪比资料,对宽频带地震资料可得到高精度瞬时频率。
-
公开(公告)号:CN113156497A
公开(公告)日:2021-07-23
申请号:CN202011143818.0
申请日:2020-10-23
申请人: 中国海洋石油集团有限公司 , 中海油研究总院有限责任公司
摘要: 本发明提供一种在相位空间利用相位解缠绕的地下介质参数反演方法,包括用解缠绕后的相位替换目标函数中采样前的波场,同时将波场方程替换为频率域波动方程的相位形式,由此将由模型空间到波场空间的映射转化为模型空间到相位空间的映射的步骤。本发明其具有以下优点:1.反演的非线性弱;2.对初始速度依赖性较小;3.可以直接反演地下的大尺度构造,分辨率高。
-
公开(公告)号:CN111208562A
公开(公告)日:2020-05-29
申请号:CN202010051543.1
申请日:2020-01-17
申请人: 中国海洋石油集团有限公司 , 中海油研究总院有限责任公司
IPC分类号: G01V1/30
摘要: 本发明公开了一种碳酸盐岩储层中火成岩侵入范围的识别方法,包括初步分析、建过完备字典和判断火成岩的侵入范围等步骤。本发明公开的一种碳酸盐岩储层中火成岩侵入范围的识别方法,根据稀疏分解思想,利用多道匹配追踪从过完备字典中获取的最优匹配子波,对侵入火成岩地层强反射轴进行追踪识别、分离,提高了计算的效率,更能保证原始地震信号的完整性,而且抗噪性增强,更加精确地确定火成岩的分布范围,改善地震资料的可靠性。
-
-
-
-
-
-
-
-
-