地震反射波成像新方法及其理论基础.pdfVIP

( D 辑) 第30 卷 第6 期 SCIENCE IN CHINA ( Series D ) 2000 年 12 月 地震反射波成像新方法及其理论基础* 黄光远 ( 山东大学数学院信息控制系, 济南250100) 摘要 提出了一个地震反射波的成像公式及其理论基础. 波的传播现象应该由瞬时 谱来刻画, 在数学上表示成三元( 即时间 空间和频率) 复函数. 介质的物理参数, 如反 射系数波和阻抗等均是空间和频率的复函数. 给出了反射波谱和介质反射系数间的关 系, 反演公式中还考虑到多次反射和结构滤波作用. 同时也指出传统的卷积模型和波 动方程反演方法中存在的问题. 各种反演公式均可利用小波变换和三基色方法作出地 下介质参数分布的剖面图, 不同颜色表示不同介质. 关键词 小波 地震勘探 反射系数 剖面图 在文献[1,2] 中已经给出此方法详细阐述和初步分析. 本文将系统地给出此方法的理论基 础和进一步发展. 1 波与介质参数的数学表示 熟知 Fourier 变换不能反映波谱随时间的变化. 但实际信号的谱是时变的. 根据现代小波 [3] 2 , ) 都可以用其二元时-频( 即瞬时)谱重构表示, 即在 理论 , 任何一个波动函数f (t) L ( 每一个瞬时都可得到一条谱线, 它反映了在这一时刻近旁的信号是由哪些简谐波组成. 考虑 到波在空间传播, 描述波现象, 应该用一个时 空和频率的三元(五维)复函数F (t, x , w)来刻画 波现象. 文献[1]中的图 1 就是某点实测地震波的瞬时幅谱. 波在空间的传播过程就是一个随着时 间和空间通道变化的频谱变化过程. 复数F(t, x , w) 的模Fr 和幅角Fa 分别表示在确定位置x 和时 间t 的小邻域内频率为w 的简谐波的波幅和相位. 介质的物理参数也与频率有关. 例如, 当波通过非均匀介质存在反射与透射. 反射系数 和透射系数均与介质性质有关, 但它们也与频率有关. 电磁波段显然如此, 不同的物体反射出 不同的颜色. 在低频段, 如声波也如此, 不同的物体反射出不同音色. 如果仅将反射系数定义 为一个数, 不看成频率的函数, 那么不同的物体反射出仅有强弱差别的同样声音. 这显然不符 合实际. 弹性波也不应例外. 因此, 反射系数应表示成空间和频率有关的二元复函数, 记为 R (x , w), 其模就是入射波与反射波波幅之比, 幅角就是两者相位差. 推广之, 与介质性质有关 的透射系数及波阻抗都应是空间和频率的二元复函数. 事实上, 透射系数就是分层介质作为 一个滤波器的频率特性. [4] 上述两个概念是本文方法的基础 , 也是与现有反射波成像方法区别的基本点. 1998-01-05 收稿, 1999-11-10 收修改稿 * 国家自然科学基金资助项目(批准号: 第6期 黄光远:地震反射波成像新方法及其理论基础 651 2 现有地震波反演方法中存在的问