一种基于叠前CRP道集的多尺度各向异性扩散滤波方法

    公开(公告)号:CN103926616B

    公开(公告)日:2016-07-13

    申请号:CN201410144181.5

    申请日:2014-04-11

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种基于叠前CRP道集的多尺度各向异性扩散滤波方法,其包括步骤:1)对输入的叠前CRP道集进行正则化处理;2)采用二维Mallat算法进行多尺度分解,并初始化分解后的每个子剖面;3)经扩散系数法求得扩散门限值,进入步骤6);进行基于扩散张量法的各向异性扩散求得扩散张量参数后,代入非线性各向异性方程进行各向异性扩散滤波,并进入步骤4);4)对每次迭代后的子剖面计算SNR、MSE和PSNR,优选出最佳的地震子剖面;5)对各向异性扩散滤波处理后的各个子剖面进行步骤4)处理,进入步骤6);6)采用二维Mallat重构算法对迭代后的各个子剖面信息进行无损重构叠前CRP道集,重构之后输出最佳的叠前CRP道集。本发明可以广泛用于各种石油勘探地震资料的处理过程中。

    一种地震数据分形保幅方法

    公开(公告)号:CN103630932B

    公开(公告)日:2016-03-30

    申请号:CN201310434890.2

    申请日:2013-09-23

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种地震数据分形保幅方法,其包括以下步骤:1)针对拓频前的地震数据体X0和拓频后的地震数据体Y,采用带通算子Bi分别相应提取它们的局部频段信息,记为BiX0和BiY,i=1,2,3…M,不同的i代表不同的带通频段;2)基于地震信号局部频段成分与整体的特征存在的相似性,在每一对数据体BiX0和BiY之间分别建立保幅关系,并相应提取局部保幅因子Φi;3)对所述步骤2)提取的多个局部保幅因子Φi进行平滑和组合处理,得到整体保幅因子Φ;4)采用整体保幅因子Φ作用于数据体Y,对全频段信号Y实现误差精度允许范围内的保幅。本发明可以广泛应用于对拓频后的地震数据进行振幅保真处理,进一步为后续地震资料岩性反演和高精度解释提供高保真数据。

    一种基于叠前CRP道集的多尺度各向异性扩散滤波方法

    公开(公告)号:CN103926616A

    公开(公告)日:2014-07-16

    申请号:CN201410144181.5

    申请日:2014-04-11

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种基于叠前CRP道集的多尺度各向异性扩散滤波方法,其包括步骤:1)对输入的叠前CRP道集进行正则化处理;2)采用二维Mallat算法进行多尺度分解,并初始化分解后的每个子剖面;3)经扩散系数法求得扩散门限值,进入步骤6);进行基于扩散张量法的各向异性扩散求得扩散张量参数后,代入非线性各向异性方程进行各向异性扩散滤波,并进入步骤4);4)对每次迭代后的子剖面计算SNR、MSE和PSNR,优选出最佳的地震子剖面;5)对各向异性扩散滤波处理后的各个子剖面进行步骤4)处理,进入步骤6);6)采用二维Mallat重构算法对迭代后的各个子剖面信息进行无损重构叠前CRP道集,重构之后输出最佳的叠前CRP道集。本发明可以广泛用于各种石油勘探地震资料的处理过程中。

    一种地震数据分形保幅方法

    公开(公告)号:CN103630932A

    公开(公告)日:2014-03-12

    申请号:CN201310434890.2

    申请日:2013-09-23

    IPC分类号: G01V1/28

    摘要: 本发明涉及一种地震数据分形保幅方法,其包括以下步骤:1)针对拓频前的地震数据体X0和拓频后的地震数据体Y,采用带通算子Bi分别相应提取它们的局部频段信息,记为BiX0和BiY,i=1,2,3…M,不同的i代表不同的带通频段;2)基于地震信号局部频段成分与整体的特征存在的相似性,在每一对数据体BiX0和BiY之间分别建立保幅关系,并相应提取局部保幅因子Φi;3)对所述步骤2)提取的多个局部保幅因子Φi进行平滑和组合处理,得到整体保幅因子Φ;4)采用整体保幅因子Φ作用于数据体Y,对全频段信号Y实现误差精度允许范围内的保幅。本发明可以广泛应用于对拓频后的地震数据进行振幅保真处理,进一步为后续地震资料岩性反演和高精度解释提供高保真数据。

    基于广义瞬时相位及P范数负模的地震弱信号处理方法

    公开(公告)号:CN103645504A

    公开(公告)日:2014-03-19

    申请号:CN201310692354.2

    申请日:2013-12-17

    IPC分类号: G01V1/28 G01V1/36

    摘要: 本发明涉及一种基于广义瞬时相位及P范数负模的地震弱信号处理方法,其包括以下内容:1)对地震信号进行识别,具体过程为:对输入的地震信号进行P阶广义希尔伯特变换并得到广义瞬时相位;求取广义瞬时相位P范数负模;根据广义瞬时相位P范数负模提取地震弱信号的相关属性参数,通过提取的属性参数对地震有效信号进行识别;2)对地震信号的噪声进行去除,得到信噪比比较高的地震信号,具体过程为:对输入的地震信号进行P阶广义希尔伯特变换并得到广义瞬时相位;对计算得到的广义瞬时相位进行平滑滤波处理,得到滤波后的广义瞬时相位信号;采用逐步迭代的方法求取地震信号的瞬时振幅;根据求取的广义瞬时相位信号和瞬时振幅重构地震信号;并对重构地震信号的合理性进行判断。

    一种基于Morlet谱勘测油气的方法

    公开(公告)号: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相比于围岩的吸收系数低百分之十以上,则表示该位置含有气层。

    地震噪声压制方法及装置
    10.
    发明授权

    公开(公告)号:CN110554434B

    公开(公告)日:2020-06-26

    申请号:CN201910768039.0

    申请日:2019-08-20

    IPC分类号: G01V1/36

    摘要: 本发明公开了一种地震噪声压制方法及装置,该方法包括根据噪声的频带范围对地震记录剖面进行带通滤波,确定第一及第二地震记录剖面,对第一地震记录剖面中每个地震道的时频谱进行频移及同步挤压小波变换;将频移时频谱的能量划分为第一能量区域及第二能量区域,确定与第一地震记录剖面对应的反射波信号,确定地震记录剖面的反射波信号。本发明对于同频同相的反射波信号与噪声,利用频移将同一时刻分散于低频频率间隔内的能量集中至高频频率区间,在高频频率区间内利用未被噪声影响的能量区域对噪声所在的能量区域进行拟合预测,从而将反射波信号与噪声较好的分离开来,取得了较好的信噪分离结果,能够提高噪声压制的效果。