利用地震干涉法合成反射地震记录研究.pdfVIP

超声学 利用地震干涉法合成反射地震记录的研究 李磊，陈浩，王秀明 超声物理与探测实验室 摘要：本文简述了地震干涉法的基本原理，对不同检波器位置的地震记录进行相关或褶 积运算而合成新的地震记录。对各向同性的水平层状介质模型进行了模拟，将地震干涉 法利用虚源合成的与采用真实源直接模拟的地震记录和道集波形进行了对比，二者在反 射波的幅度上虽稍有差异，但相位吻合较好，从而验证了地震干涉法合成反射地震记录 的可行性。 关键词：地震干涉法互相关有限差分地震记录 1 引言 物理学中，干涉通常指两个或两个以上不同信号在空间上重叠时形成新信号 的现象，如光的干涉。地震干涉法 SeismicInterferometry 是近年来新兴的一 种地震勘探方法，是干涉思想在地震资料处理中的拓展和应用，研究的对象是地 震信号或声信号。该方法的提出不仅丰富了勘探地震学和天然地震学的理论与应 用研究，且因为具有增强地震资料信噪比和分辨率、不需要震源信息并能将噪声 信号转变为有效信号等特点，其在主动源和被动源地震勘探中都有很好的应用前 景。例如能将多次波转换为一次波加以利用、从噪声信号中提取反射记录等。 该方法的理论研究最早源自1968年C1aerbout的一篇经典论文，文中首次 针对水平层状模型，对自由地表接收到的透射地震记录进行自相关，合成了相应 的自激自收反射地震记录。通过对不同的检波器记录的透射波地震信号进行互相 关，可得到一个新的包含因果与非因果部分的地震信号，其因果部分等效于以其 中一个检波器为虚震源的新的地震记录。地震干涉法最早的提出是基于互相关理 论，后来还发展了褶积型和反褶积型地震干涉技术。 本文以相关型地震干涉法为例，重点介绍了该方法的基本原理并利用模型实 验对其理论进行了验证分析，结果表明地震干涉法合成反射地震记录的可行性。 2基本原理 地震干涉法的基本原理是对不同检波器位置的地震记录在数学上进行相关 或褶积运算而合成新的地震记录。围绕地震干涉法的研究和讨论大多基于互相关 理论，并且利用格林函数来表征检波器和震源之间的场源关系。地震干涉法也被 看作是一种利用地震记录的相关运算获取格林函数的方法。 已知互相关函数是两组数据之间相似性的量度。其表达式为： f r g t+r矽f 1 ^ f 厂 f pg t I 式中。是互相关运算算子。对上式进行傅里叶变换，转换到频域为： ．461． 中国科学院声学研究所纪念建所50周年暨第五届学术交流会 2 H co 2nF co ’G co 2nIF 09 +帜国 ∥旷州 令 1 式中f t 和g t 分别代表两个不同检波器的地震记录，由傅里叶变换中 时域移位和频域相位移的对应关系的性质可知，互相关运算作用予地震记录的一 个重要性质是可以缩短传播时间或传播路径。即地震记录的时间迁移和互相关运 算是等价的。如下图1所示。 互相 l @ 警二次反射波 一 二 和一次反射波 ! P一次反 --- - -- ▲ -_ ＼ W +’V 图1互相关地震干涉法示意图 Schuster,2009 由脉冲点源产生的场叫做格林函数。两个地震检波器间的格林函数相当于在 一个检波器施加了一个万函数的脉冲激励之后，在另外的一个检波器接收到的脉 冲响应。 利用格林函数及其满足的Helmholtz方程，并应用高斯定理和万函数满足的 积分性质，可推导出频域中格林函数表示的相关型互易性方程如下： G 叩 _G 删 工lG川砂掣-G川彳 掣P2x 3 其中A，B为同一研究区域内两个独立声学状态下任意两个震源点，x为区域内任 意接收点，S为区域的边界。 相关型互易方程是地震干涉法的基础方程。方程中的被积函数部分包含单极 源和偶极源两种形式的震源，若直接应用存在困难，特别是对以背景噪声或天然 微震为源的被动源地震数据。为了能更好的将其进行应用，可以采用高频远场近 似，用单极源格林函数近似表示偶极源格林函数，将被积函数简化为只有单一核 函数的积分。高频和远场近似表示偶极源项： A ／r≈jkO xA 4 h[VG xI彳 一OG xI么 ／Or jkG xI彳 一G x B ‘。高频和远场近似可能会产生明显的幅度误差， 同理有hOTG xI曰 +≈一jkG x 但对相位没有影响。将上述远场近似代入到相关型互易方程中： G Bl彳 一G A G B l召 ‘ 2jimll彳 l ≈2fG xⅢ塑兰』竺．dzz 5 $一’ 饥 A d2x ≈2jk B +G xI LG xl 方程左边为新合成地震记录的因果和非因果格林函数部分，右边转换到时域即为 不同位置接收的地震记录的互相关运算。实际中由于震源信息和介质参数的未 ．．462．． 超声学 知，因此波场传播方向也无法确定，需要借助稳相分