一种多分量地震资料偏移成像方法技术

本发明专利技术公开一种多分量地震资料偏移成像方法

【技术实现步骤摘要】
一种多分量地震资料偏移成像方法、系统、设备及介质


[0001]本专利技术涉及地震波传播数值模拟和成像领域，特别是涉及一种多分量地震资料偏移成像方法
、
系统
、
设备及介质
。

技术介绍

[0002]起伏地表条件是地震勘探面临的第一个亟待解决的难题
。
剧烈的地形起伏会直接导致数据采集困难
、
信噪比低
、
散射噪声严重
、
静校正不准确
、
成像精度低等一系列问题
。
为消除起伏地表的影响，目前工业界最常用的做法是在偏移成像之前对地震资料进行静校正等预处理，受相关预处理方法误差的影响，该过程会对地震资料的运动学及动力学特征带来畸变，从而影响后续偏移成像的准确性
。
基于起伏地表直接进行偏移成像是解决复杂近地表地区地震成像问题的根本方法，而实现起伏地表条件下的高效高精度地震波数值模拟是其核心内容
。
尽管当前起伏地表地震波数值模拟已经取得了较为丰富的研究成果，但其在实际应用中依然本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种多分量地震资料偏移成像方法，其特征在于，包括：获取工区地质
‑
地球物理资料和工区地表高程；根据所述工区地质
‑
地球物理资料和所述工区地表高程对工区的规则网格深度域模型进行离散节点剖分，得到所述规则网格深度域模型对应的初始离散节点位置坐标；对所述初始离散节点位置坐标利用离散节点控制方程进行节点位置自适应移动，得到所述规则网格深度域模型对应的最终离散节点位置坐标；根据所述最终离散节点位置坐标利用径向基函数有限差分方法构建弹性波方程离散节点有限差分数值求解格式；利用滤波器根据所述最终离散节点位置坐标对所述弹性波方程离散节点有限差分数值求解格式进行求解，得到每一炮每一时刻的震源端波场和检波器波场；根据所述震源端波场和检波器波场利用离散节点弹性波场纵横波分离径向基函数有限差分数值计算格式确定离散节点的纯净纵波场矢量和纯净横波场矢量；基于时间一致性原理，根据所述离散节点的纯净纵波场矢量和纯净横波场矢量利用偏移成像条件对震源端的纯净波场矢量和检波器端的纯净波场矢量进行成像，得到每一炮的单炮偏移剖面
。2.
根据权利要求1所述的多分量地震资料偏移成像方法，其特征在于，根据所述工区地质
‑
地球物理资料和所述工区地表高程对工区的规则网格深度域模型进行离散节点剖分，得到所述规则网格深度域模型对应的初始离散节点位置坐标，具体包括：利用所述工区地质
‑
地球物理资料和所述工区地表高程基于工区的规则网格深度域模型的物理边界，逐层进行离散节点剖分，得到物理边界上的离散节点初始空间位置坐标
、
离散节点处的纵波速度
、
横波速度和密度；根据所述物理边界上的离散节点初始空间位置坐标
、
离散节点处的纵波速度
、
横波速度和密度基于自适应无网格剖分方法，得到所述规则网格深度域模型对应的初始离散节点位置坐标
。3.
根据权利要求1所述的多分量地震资料偏移成像方法，其特征在于，对所述初始离散节点位置坐标利用离散节点控制方程进行节点位置自适应移动，得到所述规则网格深度域模型对应的最终离散节点位置坐标，具体包括：对所述初始离散节点位置坐标利用离散节点控制方程，通过
Gauss
‑
Seidel
迭代法进行求解，得到所述规则网格深度域模型对应的最终离散节点位置坐标
。4.
根据权利要求1所述的多分量地震资料偏移成像方法，其特征在于，所述离散节点弹性波场纵横波分离径向基函数有限差分数值计算格式的表达式为：性波场纵横波分离径向基函数有限差分数值计算格式的表达式为：其中，为离散节点的纯净纵波场矢量，
v
p
为最终离散节点位置所对应的纵波速度，为水平二阶导数系数，
u
j
为弹性波场的水平分量，为混合二级导数的径向基函数有限差
分系数，
w
j
为弹性波场的垂直分量，为笛卡尔坐标系下
x
方向的单位矢量，为垂直二阶导数系数，为笛卡尔坐标系下
z
方向的单位矢量，为离散节点的纯净横波场矢量，
v
s
为最终离散节点位置所对应的横波速度，
i
和
j
为均为索引号，
n
为参与差分计算的节点总数
。5.
一种多分量地震资料偏移成像系统，其特征在于，包括：获取模块，用于获取工区地质
‑
地球物理资料和工区地表高程；...

【专利技术属性】
技术研发人员：谷丙洛，乔运聪，刘兴农，张鑫怡，黄建平，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：