改进单程波衰减介质地震波数值模拟方法.pdfVIP

209 SPG/SEG 深圳 2011 国际地球物理会议论文集 改进的单程波衰减介质地震波数值模拟方法 何兵红 吴国忱 （中国石油大学(华东)地球科学与技术学院） 摘要 地球实际并非完全弹性介质，地层吸收致使地震波在地层中传播存在能量衰减及速度频散。本 文采用基于 Futterman 模型的分步傅里叶（SSF）波场延拓方法,实现了单程波动方程非零偏移距衰减 介质地震波传播数值模拟。针对基于检波点波场下延原理进行炮记录计算效率低的问题，我们根据反 射波场物理机制采用反射波正向传播理论，提高了计算的效率。通过数值模拟对比了单程波动方程和 双程波动方程的特点，并探索衰减介质中地震波传播的特点：振幅减小，能量衰减，主频向低频方向 移动，频带变窄，时间分辨率降低。 关键词 单程波 衰减介质 数值模拟 计算效率 的问题，本文对波场记录方式进行改进，采用 1 前言 了反射波场正向传播原理，即：按照波场正向 地震波数值模拟在地震资料处理解释中占 传播方向对反射波场进行波场延拓，一次性接 据重要地位。波动方程数值模拟可以分为双程 收整个检波器排列波场。 波动方程数值模拟和单程波动方程数值模拟。 2 方法原理 双程波动方程会引入多次波，计算速度慢，占 用内存大。单程波数值模拟只含有一次反射波， 2.1 衰减介质的波场延拓算子 计算速度快，有利于衰减介质地震波场特征研 为了提高使用单程波法进行波场延拓对速 究。对于弹性介质地震响应特征研究已经非常 度横向变化的适应性，Stoffa在相移偏移的基础 成熟(Robertsson 1994, Stekl,R.G. Pratt 1998) ，但 上提出了分步Fourier方法[4] 。该方法是基于速 地球介质并非完全弹性介质，诸如孔隙介质、 度场分裂的思想，将整个速度场视为常速背景 各向异性介质、粘弹性介质等各种对地球介质 和变速扰动。对于弹性介质分步傅里叶法分两 近似的介质理论应运而生。本文采用粘弹性理 步进行： 论描述地球介质，利用介质的品质因子Q对吸   P k , z , P k , z , exp iA z 1  x i +1   x i   1 i  收衰减进行定量表示。衰减介质波场延拓算子 (1) P x, z , P x, z , exp iA z          方面，本文采用基于Futterman模型[1] 的衰减介 i1 1 i1 2 i 质分步傅立叶法波场延拓理论。 2 2   2 2