一种基于lowrank有限差分的弹性逆时偏移方法及系统

    公开(公告)号:CN110857998B

    公开(公告)日:2021-11-05

    申请号:CN201810967971.1

    申请日:2018-08-23

    发明人: 韩冬 李博 许璐

    IPC分类号: G01V1/28 G01V1/36

    摘要: 公开了一种基于lowrank有限差分的弹性逆时偏移方法及系统,该方法及系统包括:构建lowrank有限差分系数,基于构建的lowrank有限差分系数的离散化波动方程实现正向及反向弹性波场延拓,利用获得的正向及反向弹性波场,利用标量积成像条件获得PP及PS成像剖面。本发明通过利用平面波理论及lowrank分解方法构建了时空兼顾的有限差分算子,能够满足时间和空间方向高精度弹性波场延拓,并充分利用多分量地震资料,而且在延拓过程中可得矢量的纵波及横波,保持了地下波场的矢量特性,在成像过程中避免了常规弹性逆时偏移中不同波型串扰及转换波极性校正的问题,提高了弹性逆时偏移成像的精度。

    基于低秩分解的粘声逆时偏移方法及系统

    公开(公告)号:CN112649874A

    公开(公告)日:2021-04-13

    申请号:CN201910969338.0

    申请日:2019-10-12

    发明人: 韩冬 许璐 崔月

    IPC分类号: G01V1/36

    摘要: 公开了一种基于低秩分解的粘声逆时偏移方法及系统。该方法可以包括:基于低秩分解的粘声介质正向波场延拓,获得震源波场;基于低秩分解的粘声介质反向波场延拓,获得检波波场;针对检波波场进行滤波处理,获得优化检波波场;根据震源波场与优化检波波场进行归一化互相关成像,获得成像值。本发明基于分数阶粘声波动方程,并借助平面波分解、低秩分解等方法在时间‑空间域中构建兼顾时间和空间方向高精度粘声低秩分解波场延拓算子,并通过自适应滤波算法实现粘声反向波场延拓中高频噪音的有效压制,确保粘声波场延拓精度及反向波场延拓的稳定性,利于偏移成像质量提高。

    基于交错网格低秩有限差分的三维弹性波场数值模拟方法及系统

    公开(公告)号:CN112630823A

    公开(公告)日:2021-04-09

    申请号:CN201910950906.2

    申请日:2019-10-08

    发明人: 崔月 王鹏燕 韩冬

    IPC分类号: G01V1/28

    摘要: 一种基于交错网格低秩有限差分的三维弹性波场数值模拟方法及系统。该方法包括:包括:构建交错网格弹性递归时间积分算子;对交错网格弹性递归时间积分算子进行低秩分解,获得基于低秩分解的弹性递归时间积分算子;基于基于低秩分解的弹性递归时间积分算子,计算基于低秩分解的交错网格有限差分系数;基于所述低秩分解的交错网格有限差分系数获得三维解耦的弹性波方程的数值离散表示,利用离散表示的弹性波方程进行数值模拟,获得纵横波矢量波场。本发明能在高波数处拟合频散关系曲线,从而可实现弹性波场延拓得到三维纵横波矢量波场,其时间和空间方向上的精度高,利于后续的处理解释。

    粘声介质波场延拓和波场分解方法、存储介质及计算设备

    公开(公告)号:CN112444863A

    公开(公告)日:2021-03-05

    申请号:CN201910799470.1

    申请日:2019-08-28

    发明人: 韩冬 李博 崔月 许璐

    IPC分类号: G01V1/30 G01V1/28

    摘要: 本发明涉及粘声介质地震波场延拓和波场分解方法、存储介质及计算设备,方法包括如下步骤:S1:基于时间‑空间域分数阶粘声波动方程,得到时间‑波数域分数阶粘声波动方程;S2:通过对时间‑波数域分数阶粘声波动方程进行求解,得到粘声递归时间方程;S3:通过对粘声递归时间方程中的粘声递归时间积分算子进行低秩分解,得到粘声低秩分解波场延拓算子;S4:将算子代入至粘声递归时间方程并进行反傅里叶变换,得到时间‑空间域的粘声介质地震波场延拓方程;S5:根据粘声介质地震波场延拓方程,构建粘声介质地震波场;S6:利用希尔伯特变换对粘声介质地震波场进行波场分解,以得到上行粘声介质地震波场和下行粘声介质地震波场。

    不规则地震数据全频带重建方法

    公开(公告)号:CN110764135A

    公开(公告)日:2020-02-07

    申请号:CN201810830925.7

    申请日:2018-07-26

    发明人: 许璐 李博 韩冬 张林

    IPC分类号: G01V1/28

    摘要: 本发明公开一种不规则地震数据全频带重建方法,包括步骤:步骤1:采集原始不规则地震数据;步骤2:获得频率域数据体;步骤3:获得频率波数域数据体;步骤4:确定优势频带和阀值;步骤5:将频谱振幅小于阀值的成分置零,在频率波数域进行加权,进行空间方向的傅里叶反变换,得到频率域重建数据体;步骤6:将已采样数据用频率域数据体替换,得到频率域更新数据体;步骤7:将频率域更新数据体重复步骤3到步骤6,直至满足迭代次数或误差终止条件;步骤8:将所有频率成分进行迭代,并输出最终插值后的时间空间域重建道集。本发明方法既保证了计算精度,又减少了计算量,保护了地震数据中的弱反射信号,有效提高小断层和微服构造成像质量。

    吸收衰减介质叠前深度逆时偏移成像方法及系统

    公开(公告)号:CN110658558A

    公开(公告)日:2020-01-07

    申请号:CN201910912183.7

    申请日:2019-09-25

    IPC分类号: G01V1/36

    摘要: 一种吸收衰减介质叠前深度逆时偏移成像方法及系统,方法包括:步骤1:在地震数据的有效频带范围内,给定L个应变拟合参数τσl的初始值;步骤2:基于L个应变拟合参数τσl的初始值计算当前计算频率下的积分函数值;步骤3:根据拟合参数表达式获得当前Q值对应的应力拟合参数τ;步骤4:针对Q值场数据体内的所有Q值点,重复步骤1至步骤3,得到拟合参数数据体;步骤5:获得黏声介质波场正传算子和黏声介质波场反传算子,根据正传算子和反传算子进行波场的正演模拟,获得炮点正传波场与检波点反传波场;步骤6:利用成像条件针对任一时刻的炮点正传波场与检波点反传波场进行成像。本发明在保证成像精度条件下大幅提高叠前深度逆时偏移成像的效率。