一种基于非线性挤压S时频变换的弱信号检测方法

    公开(公告)号:CN107918146B

    公开(公告)日:2018-10-30

    申请号:CN201710612806.X

    申请日:2017-07-25

    发明人: 高静怀 李振

    IPC分类号: G01V1/28

    摘要: 本发明公开了一种基于非线性挤压S时频变换的弱信号检测方法,包括:1)计算信号的S变换;2)对上面得到的系数进行阈值处理;3)进行同步挤压操作;4)弱信号检测。该方法从常规的S变换出发,将得到小于阈值的系数进行筛选,然后进行同步挤压操作,得到分辨率更高的时频表示。对于能量比较弱的信号的检测,我们利用提出的非线性挤压S时频变换,就可以在时频域比较清晰地显示出所有的有效信号成分。该技术方案易于实现,可操作性强。

    基于高分辨率时频分析和一致性度量的非均质性刻画方法

    公开(公告)号:CN107272063B

    公开(公告)日:2018-10-30

    申请号:CN201710543985.6

    申请日:2017-07-05

    IPC分类号: G01V1/30

    摘要: 本发明公开了一种基于高分辨率时频分析和一致性度量的非均质性刻画方法,包括步骤1:计算三参数小波变换系数;步骤2:计算重排准则;步骤3:计算重排后的时间‑频率域系数;步骤4:阈值操作,得到有效信号对应系数(有效信号分布空间);步骤5:选取合适的频率,提取单频地震数据体;步骤6:计算步骤5中单频数据体的一致性度量相干体,进行地震数据的非均质性刻画。本发明可以为地震资料处理和解释提供有用的信息,在复杂隐蔽油气藏勘探中,精确的时频表示对于指示细小的地层结构和检测含油气储层至关重要。本发明结合一种高分辨率的时频分析技术(TPW同步挤压变换)和基于一致性度量的相干体技术,进行地震数据非均质性刻画。

    基于二阶同步挤压小波变换的地震资料时频分析方法

    公开(公告)号:CN107132576B

    公开(公告)日:2018-10-30

    申请号:CN201710543582.1

    申请日:2017-07-05

    IPC分类号: G01V1/30

    摘要: 本发明公开一种基于二阶挤压小波变换的地震资料时频分析方法,首次在连续小波变换的基础上,引入二阶同步挤压算子,该变换通过对小波变换的系数在时间‑尺度域进行重排,从而压制了能量扩散,该方法可以得到一种具有更高时频分辨率的分析结果。相比于之前的同步挤压变换,改方法更加适合对地震信号这种带有较强频率调制的信号进行分析,进一步将其应用于实际地震数据的时频分析中,可以更加准确的界定地下碳氢储层的边界,更精细的刻画河道与断层等地质信息,从而有利于地震资料的进一步解释和井位确定。

    一种基于VMD-TFPF压制地震勘探随机噪声方法

    公开(公告)号:CN107589454A

    公开(公告)日:2018-01-16

    申请号:CN201710612826.7

    申请日:2017-07-25

    发明人: 高静怀 李振

    IPC分类号: G01V1/36

    摘要: 本发明公开了一种基于VMD-TFPF压制地震勘探随机噪声方法,包括:1)分解含噪信号;2)计算互相关系数,判断要进行滤波的模态分量;3)对筛选的模态进行滤波;4)计算得到去噪后的信号。本发明能够有效衰减地震随机噪声,并且在滤波过程能够很好地保持有效信号幅值;该技术方案的算法内容易于实现,可操作性好;同时避免了时频峰值滤波窗口大小的选择问题,适用于复杂信号去噪的需求。

    一种信息熵滤波器与地震资料随机噪声衰减方法

    公开(公告)号:CN107589453A

    公开(公告)日:2018-01-16

    申请号:CN201710612809.3

    申请日:2017-07-25

    发明人: 高静怀 杨涛 张兵

    IPC分类号: G01V1/36

    摘要: 本发明公开了一种信息熵滤波器与地震资料随机噪声衰减方法,包括:1)计算梯度结构张量;2)计算地层横向不连续性度量;3)阈值参数的确定;4)结构自适应滤波尺度参数的选择;5)信息熵滤波器处理。本发明基于能同时兼顾脉冲噪声和随机高斯白噪声的衰减出发,首先定义了一种新的滤波器,称为信息熵滤波器,结合梯度结构张量,线型置信度量和结构自适应滤波,能够有效衰减地震资料中脉冲噪声和随机噪声,增强地震同相轴的空间一致性,并且能够保护有效信号和断层、裂缝等地层边缘和细节信号结构,技术方案易于实现,可操作性强。

    一种测井资料统计特性约束的波阻抗反演粒子群优化算法

    公开(公告)号:CN107390269A

    公开(公告)日:2017-11-24

    申请号:CN201710543074.3

    申请日:2017-07-05

    发明人: 高静怀 张兵

    IPC分类号: G01V1/30

    摘要: 本发明公开了一种测井资料统计特性约束的波阻抗反演粒子群优化算法,具体包括:1)对地震数据进行相位校正,使其为零相位记录;2)提取子波振幅谱并得到地震子波;3)基于数据体或剖面内已有的测井资料,对速度信息进行分析和统计,建立用于反演反射系数的可行域,即反射系数的搜索空间;4)用本发明提出的改进的粒子群优化算法反演反射系数;5)基于反演的反射系数计算波阻抗。本发明为克服反射系数反演的多解性,提高了反演结果的分辨率,给反射系数加入概率分布的约束信息,保证了反演的反射系数概率分布的不变性。本发明用于反演的概率分布先验信息是完全来自于实际测井数据,从而保证了反演结果的可信度并且降低了多解性。

    S变换域基于变分原理的Q值估计与补偿方法

    公开(公告)号:CN107356964A

    公开(公告)日:2017-11-17

    申请号:CN201710544007.3

    申请日:2017-07-05

    IPC分类号: G01V1/30

    摘要: 本发明公开一种S变换域基于变分原理的Q值估计与补偿方法,将S变换引入非平稳地震记录褶积模型,首先将地震记录进行S变换,相比于时间域以及频率域方法,在时频域地震记录的非平稳特性可得更加清晰的刻画,基于变分原理进行Q值估计,不需要对地震记录进行加窗处理,可以减小窗函数对Q值估计结果的影响,此外,由于是对整个地震记录进行整体处理,可以避免在处理反射地震数据时吸收与散射的耦合效应,抗噪性能更强,反演的参数精度更高。

    一种保地层结构的地震波衰减补偿方法

    公开(公告)号:CN107132579A

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

    申请号:CN201710543571.3

    申请日:2017-07-05

    IPC分类号: G01V1/36

    CPC分类号: G01V1/362 G01V2210/584

    摘要: 一种保地层结构的地震波衰减补偿方法,采用地震包络局部峰值作为地层结构约束来构造自适应于地层结构的分子时窗,然后采用非线性压缩映射子波振幅谱估计方法从地震道分子分解时频谱中估计正比于时变子波振幅谱的分量,接着计算质心频率,并对质心频率做筛选,减少子波干涉和噪声的影响,然后用质心频率偏移法估计得到稳健的地下介质Q场,并计算得到各时窗对应的吸收补偿滤波器,最后用这些滤波器对相应的时频谱进行加权补偿,再反变换回时间域,就得到衰减补偿后的地震记录。本发明能够在有效补偿地震波衰减、提高地震资料纵向分辨率和深层弱信号识别能力的同时,很好地保持地震资料的横向连续性和复杂结构特征。

    基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法

    公开(公告)号:CN104880730B

    公开(公告)日:2017-04-26

    申请号:CN201510140952.8

    申请日:2015-03-27

    发明人: 高静怀 王平

    IPC分类号: G01V1/28

    摘要: 本发明公开了一种基于Synchrosqueezing变换的地震资料时频分析和衰减估计方法,首次将新的时频分析工具Synchrosqueezing变换用于地震资料时频分析,该变换通过对变换域系数的重排,获得一个更加聚集的时频表示,时频分辨率大大提高,将其用于实际地震资料分析和致密砂岩模型含气性检测,能够准确界定储层的位置,指示河道与断层等地质结构,进而有利于进一步的资料解释和井位确定;提出基于Synchrosqueezing变换的地震衰减估计方法,并给出具体实现流程,对某油田致密砂岩储层三维数据体的衰减估计结果和钻井结果有着较好的一致性,该方法可以帮助地质人员指示含气储层,确定钻井位置。

    一种倾斜层状粘弹性介质中波场正演模拟方法

    公开(公告)号:CN103149586B

    公开(公告)日:2016-01-20

    申请号:CN201310044384.2

    申请日:2013-02-04

    IPC分类号: G01V1/28

    摘要: 本发明公开了一种倾斜层状粘弹性介质中波场正演模拟方法,该方法针对地质模型为含倾斜界面的粘弹性层状介质,将点源激发的球面波分解为平面波,利用矢量波动方程不依赖于坐标系这一性质,在固定坐标系和随倾斜界面变化的动坐标系下交替地讨论波的传播。首先由传统的反射和透射系数推导出单个倾斜界面上平面波的反射和透射系数,得到该界面上的反射及透射波;然后递推求得平面波经过多个倾斜界面时的反射和透射系数。最后由传播后的各平面波合并得到点源激发波场,合并时采用了快速的振荡积分计算方法和一条适合粘弹性介质的积分路径。