-
公开(公告)号: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道集。本发明基于量子力学的散射理论,提出深水层间多次波消除技术,填补了地震资料深水层间多次波消除的空白。
-
公开(公告)号: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道集。本发明基于量子力学的散射理论,提出深水层间多次波消除技术,填补了地震资料深水层间多次波消除的空白。
-
公开(公告)号:CN109975868B
公开(公告)日:2020-03-03
申请号:CN201910230946.X
申请日:2019-03-26
申请人: 中国科学院地质与地球物理研究所
摘要: 本发明涉及一种基于泰勒展开的波动方程鬼波压制方法,其特征在于包括以下步骤:(1)将含有鬼波的地震数据表示为一次波和鬼波的叠加,得到由鬼波算子G表示的一次波表达式;(2)将鬼波算子G基于泰勒展开,得到由波场延拓矩阵F表示的一次波表达式;(3)计算得到波场延拓矩阵F,并将其代入步骤(2)中的一次波表达式,得到鬼波压制后的一次波。本发明可以广泛应用于地震数据的鬼波压制领域。
-
公开(公告)号:CN109975868A
公开(公告)日:2019-07-05
申请号:CN201910230946.X
申请日:2019-03-26
申请人: 中国科学院地质与地球物理研究所
摘要: 本发明涉及一种基于泰勒展开的波动方程鬼波压制方法,其特征在于包括以下步骤:(1)将含有鬼波的地震数据表示为一次波和鬼波的叠加,得到由鬼波算子G表示的一次波表达式;(2)将鬼波算子G基于泰勒展开,得到由波场延拓矩阵F表示的一次波表达式;(3)计算得到波场延拓矩阵F,并将其代入步骤(2)中的一次波表达式,得到鬼波压制后的一次波。本发明可以广泛应用于地震数据的鬼波压制领域。
-
公开(公告)号:CN102169188A
公开(公告)日:2011-08-31
申请号:CN201010606151.3
申请日:2010-12-15
申请人: 中国海洋石油总公司 , 中海石油研究中心 , 中国石油大学(北京)
IPC分类号: G01V1/30
摘要: 本发明涉及一种基于Morlet谱勘测油气的方法,其包括以下步骤:1)对地震记录道目的层中的地震信号进行傅里叶分析,得到地震信号的频谱范围;2)利用连续小波变换对地震信号进行分析,得到小波变换系数Wx(b,s);3)定义不同尺度因子s下的瞬时振幅为,对于每一时间点b0,得到不同尺度因子ai所对应的瞬时振幅值{Ai(b0,1/ai)i=0,1,...,N},利用Lagrange样条进行插值,b点处任一频率的振幅值都可由对应的a值代入得到,这些离散点的连线,就构成了点b的频谱,即Morlet谱;4)根据步骤3)中得到的任一时间点处的瞬时振幅,得到每一点的吸收系数;如果该点的吸收系数ξI相比于围岩的吸收系数低百分之十以上,则表示该位置含有气层。
-
公开(公告)号:CN103207409B
公开(公告)日:2016-01-06
申请号:CN201310133893.2
申请日:2013-04-17
IPC分类号: G01V1/28
摘要: 本发明涉及一种频率域全波形反演地震速度建模方法,包括以下步骤:1)获取原始地震炮集记录、震源子波信息和反演所用的初始模型;2)分析步骤1)获得的信息,基于正演算法和优化算法,确定基本的反演参数和从低频到高频的全波形反演框架;3)对不同的频率,计算获得最合适的频域正反演模型网格;4)在低频反演时采用主成分分析方法,压缩参与反演的数据维度;5)判断对应不同频率的投影矩阵维度是否满足阈值转换标准,当满足该转换标准时,进入下一步,否则回到步骤4);6)引入震源编码方法,并利用随机相位编码方法压制串扰噪声;7)判定是否满足迭代截止条件,如果满足迭代截止条件,则进入下一步;否则回到步骤6);8)如果没有完成所有频率的反演,则回到步骤3),直至完成所有频率的反演,获得最后的速度模型,输出速度模型。
-
公开(公告)号:CN103207410A
公开(公告)日:2013-07-17
申请号:CN201310134039.8
申请日:2013-04-17
IPC分类号: G01V1/28
摘要: 本发明涉及一种针对崎岖海底的混合网格模型建立方法,包括以下步骤:1)读入崎岖海底的水深、计算区域的尺度、介质物性参数分布信息;2)根据崎岖海底的水深、计算区域的尺度、介质物性参数分布信息构建出整个模型框架;3)采用正方形网格剖分步骤2)构建出的整个模型框架;4)采用三角网格剖分崎岖海底附近区域;5)根据网格节点邻近网格单元介质的物性参数分布、几何面积信息,采用面积加权平均法计算网格节点上的物性参数;6)将网格节点上的物性参数、计算参数信息存储并输出。本发明能够更准确地刻画崎岖海底的形态,提高计算效率,保证算法的稳定性,同时可以避免边界绕射和角点散射等问题。本发明可以广泛用于崎岖海底条件下的地震正演模拟、层析成像、偏移成像、速度建模等领域中。
-
公开(公告)号:CN104536045A
公开(公告)日:2015-04-22
申请号:CN201510023338.3
申请日:2015-01-16
摘要: 本发明涉及一种基于子波处理的鬼波压制方法,包括以下步骤:1)输入海上地震资料信息;其中,地震资料信息中包含未压制鬼波的地震数据和三种地震子波,三种地震子波包括含电缆鬼波和震源鬼波的地震子波、仅含震源鬼波的地震子波和不含鬼波的地震子波;2)通过地震子波对未压制鬼波的地震数据进行鬼波压制;3)输出同时压制了电缆鬼波和震源鬼波的地震数据。本发明采用了子波匹配的算法压制鬼波,从地震资料处理入手压制鬼波,而不是从采集入手压制鬼波,因此处理过程易于实现,节约成本。本发明能够压制海上地震资料中的震源鬼波和电缆鬼波,可以广泛应用于海洋地震勘探过程中。
-
公开(公告)号:CN102636810B
公开(公告)日:2014-05-14
申请号:CN201210097757.8
申请日:2012-04-05
摘要: 本发明涉及一种基于多次聚焦的噪音压制与成像方法,其特征在于:包括如下步骤:1)利用处理系统提供的分选功能将预处理之后的地震数据分选为共偏移距道集;2)基于共偏移距道集实施多次聚焦属性参数搜索,得到能够描述地震波局部特征的属性参数;3)基于上述属性参数构造多次聚焦局部反射线段;4)基于多次聚焦局部反射线段构造等旅行时面的多次聚焦算子响应;5)共偏移距剖面内所有的多次聚焦算子响应累加构成多次聚焦处理后数据体;6)将多次聚焦处理后数据体进行共中心点叠加,得到多次聚焦零偏移距成像剖面。本发明不但可以提供高质量的零偏移距成像剖面,通过充分利用菲涅耳体内的地震波能量,还可以输出信噪比大幅度提高的叠前道集。
-
公开(公告)号:CN103207409A
公开(公告)日:2013-07-17
申请号:CN201310133893.2
申请日:2013-04-17
IPC分类号: G01V1/28
摘要: 本发明涉及一种频率域全波形反演地震速度建模方法,包括以下步骤:1)获取原始地震炮集记录、震源子波信息和反演所用的初始模型;2)分析步骤1)获得的信息,基于正演算法和优化算法,确定基本的反演参数和从低频到高频的全波形反演框架;3)对不同的频率,计算获得最合适的频域正反演模型网格;4)在低频反演时采用主成分分析方法,压缩参与反演的数据维度;5)判断对应不同频率的投影矩阵维度是否满足阈值转换标准,当满足该转换标准时,进入下一步,否则回到步骤4);6)引入震源编码方法,并利用随机相位编码方法压制串扰噪声;7)判定是否满足迭代截止条件,如果满足迭代截止条件,则进入下一步;否则回到步骤6);8)如果没有完成所有频率的反演,则回到步骤3),直至完成所有频率的反演,获得最后的速度模型,输出速度模型。
-
-
-
-
-
-
-
-
-