基于广义S变换的时频域谱模拟反褶积

摘要

常规提高地震分辨率的方法,如反褶积通常假设反射系数是白噪声(Yu和Zhao,2002)且子波是不随时间和空间变化的.然而由于薄互层的多周期沉积旋回特征,使得实际的反射系数并不是随机序列.因此,简单地把地震数据的谱当做子波的谱来拓宽地震信号的频带是不准确的(Zhao et al,1996,TangBowen et al,2010),也不满足对薄互层和岩性储层勘探的技术要求.实际上,地震剖面通常被认为是子波和反射系数(Reflection Coefficient,RC)的褶积.为了消除地震信号中子波对反射系数的影响,需要从地震信号中提取子波,进而消除子波的影响.(Yun和Ding,2005;Sun,2000a;Li et al,2007;Lietal,2005)等人提出的子波提取方法都可以用于消除子波的影响，(Rosa和Ulrych,1991)在假设了波的振幅谱是平滑的基础上，通过快速傅立叶变换，提出了利用谱模拟方法从地震信号中提取地震子波的振幅谱。考虑到了波随时间变化的特点，提出一个新的方法来代替傅立叶变换方法，即基于广义S变换在时频域中获得更加准确有效的时空变子波的谱模拟反褶积方法。通过理论模型和实际数据处理验证了方法的可行性和实际资料处理效果。这种方法具有更好的时频分辨率，可以很好的模拟子波随着时间空间的变化特征，有效提高地震信号的分辨率。合成记录模型试验和实际的地震数据处理都证明了这种方法的正确性和实用性，其结果能显示出地震数据的更多细节信息。