[发明专利]一种三维分步傅立叶粘滞声波深度偏移方法

说明书

技术领域：

本发明涉及一种三维分步傅立叶粘滞声波深度偏移方法，属于地震勘探数据处理过程中的偏移成像技术，具体说是一种对地震波的吸收和频散进行补偿的深度偏移方法。

背景技术：

现有的商业软件使用的粘弹性补偿是目前工业界仍广泛使用的反Q滤波技术。但地震波的衰减不仅与传播时间有关，与子波频率、速度场分布、射线路径、地下构造展布等都有关系，仅利用反Q滤波不能恢复地震波场的本来特征。随着勘探技术的不断发展，地震成像的目标由完成构造成像任务逐步转移到寻找低幅构造以及岩性异常检测上来。地下介质的粘弹性行为，作为影响地震波传播过程中振幅、频率、相位等变化的主要因素，需要在成像过程中考虑和补偿。

发明内容

本发明的目的是提供一种三维分步傅立叶粘滞声波深度偏移方法，以粘滞性声波波动方程为理论基础，将地震波的吸收和频散考虑到偏移延拓算子中去，在波场延拓同时，补偿地震波的吸收和频散。

本发明的目的是这样实现的：一种三维分步傅立叶粘滞声波深度偏移方法，主要利用三维分步傅立叶粘滞声波深度偏移技术，处理步骤为：

(1)对叠前地震数据进行预处理，建立深度速度模型；

(2)根据实际地质情况，使用频谱比法求取品质因子Q值；

(3)输入叠前或叠后地震数据体、深度速度模型、品质因子Q值模型，输出粘滞声波波动方程叠前或叠后深度偏移结果；

(4)和常规波动方程深度偏移结果进行对比，检验结果正确性；

(5)将对比的两个数据体都转换到时间域，进行频谱分析，检验分析粘滞声波波动方程深度偏移方法提高分辨率的能力。

频率-空间域粘滞声波方程：

∂2P(x,y,z,ω)∂x2+∂2P(x,y,z,ω)∂y2+∂2P(x,y,z,ω)∂z2+ω2M(x,y,z,ω)/ρP(x,y,z,ω)=0]]>

其中ω是角频率，ρ是介质密度，M(x，z，ω)是复粘弹性模量。