慢度法全波场模拟及其水平慢度积分的实现  

Full-wavefield modeling by slowness method and realization of horizontal slowness integration.

作　　者：张繁昌[1] 印兴耀[1] 赵剑[1] 

机构地区：[1]石油大学(东营)地球资源与信息学院

出　　处：《石油地球物理勘探》2003年第6期597-602,共6页Oil Geophysical Prospecting

摘　　要：叠前地震数据的反演需要以某种正演算法为基础。本文从弹性动力学方程出发,采用慢度法计算层状介质的叠前地震记录。同反射率法一样,慢度法也能够提供地震波场的完全解,包括地震波在地层中传播时形成的一次波、多次波和转换波信息,而且其中间结果容易解释。具体实现步骤为:①通过傅立叶变换和汉克尔交换将本构方程和弹性动力学方程变换到频率-慢度域;②在频率-慢度域得到反射、透射系数矩阵递推公式;③分别对频率和慢度进行积分变换得到时-空域的地震道集。文中还讨论了慢度积分计算中解的稳定性问题。The inversion of pre-stack seismic data necds to take some forward algorithms as a basis. Start-ing from elastic dynamic equations, the paper adopted slovvness method to compute the pre-stack seismic records in layered medium. The slovvness method can provide the complete solu-tion for seismic wavefield like reflectivity method , including the information of primary reflection, multiple and converted wave resulted from propa-gation of seismic wave in beds and inlermediate results can be easily interpreted. The concrete steps of realization is as follows:①the constitu-tive equations and elastic dynamic equations are transformed into frequency-slowness domain by Fourier transform and Hankel transform;②the recursive formulas of reflectivity and transmissivi-ty matrixes are resulted in frequency-slowness domain;③the seismic data gathers can be resulted from integral transformation for frequency and slovvness separately. The stability problem on so-lution in computing slowness is also discussed in the paper.

关 键 词：叠前地震数据 慢度法 水平慢度积分 全波场模拟 地震波 地震勘探 

分 类 号：P631.4[天文地球—地质矿产勘探]