用独立分量分析方法实现地震转换波与多次反射波分离.pdfVIP

维普资讯 第 23卷第 1期 防 灾 减 灾 工 程 学 报 Vo1．23No．1 2003年 3月 JournalofDisasterPreventionandM itigationEngineering M ar． 2003 用独立分量分析方法实现地震转换波 与多次反射波分离 刘喜武h，刘 洪 ，郑天愉 (1．中国科学院地质与地球物理研究所 ，中国 北京 100101；2．中国科学院研究生院，中国 北京 100039) 摘要：独立分量分析 (ICA)是新兴的一种统计学方法。其 目的是寻求对非高斯分布数据进行有效表示，使得各个基 分量在统计学意义上独立 ，或者尽最大可能独立。这种表示意在获取数据的基本结构，可有效实现特征提取和信号 分离 。本文概述 ICA 的基本理论和快速算法 ，并在分析地震记录特点的基础上 ，阐明采用 ICA方法可 以实现对沉 积地区远震记录中转换波和多次反射波的分离。研究结果表明，ICA方法可有效地分离地震转换波和多次反射波 ， 并 由此获得较为合理的地下间断面初步解释结果 ，从而有利于扩展接收函数技术的应用 。 关键词 ：独立分量分析；地震转换波 ；多次反射波；信号分离方法 中图分类号：P315．3 文献标识码 ：A 文章编号：1672—2132(2003)01—0011—09 ICA)基于高阶统计特性 ，使分解 出的各个分量之间 0 引言 相互独立[1]。ICA是盲源分离技术近年来的一项重 要进展 ，也是最近在 国际上掀起热潮 的研究领域 。在 地壳 、上地幔 内部具有强速度间断面或过渡带 ， 声音信号分离、生物医学信号处理、图像特征提取、 如莫霍面、410km 间断面和 670km 间断面 。经过这 稀疏编码、财务时间序列分析 以及雷达通讯等方面 些 间断面的远震转换波 (P—S)具有较强信号 。利用 独立分量分析方法都 已得到很好 的应用 。 密集布设 的地震 台站记录大量地震波 ，可 以得到多 ICA 由Jutten和 Herault提 出[3]，目的是寻求 次经过 同一转换 “点”(面元)的记录 ，再利用远震体 统计上独立的非高斯分布数据 的线性变换 。在 BSS 波 的径 向分量反褶积垂 向分量 ，扣除震源和仪器因 问题 中，观测信号是几个相互独立源信号 的线性组 素 ，可以得到接收函数 。接收函数保 留了地下介质信 合 ，寻找一个线性变换 ，使得变换后得到估计信号在 息，通过动校正 ，将经过 同一转换面元 的接收函数叠 统计上尽可能的独立 ，从而仅从线性混合 的观测信 加 ，增强有效信息 ，抑制干扰 。通过追踪这些叠加 的 号 中就可 以分离 出原始的源信号 。ICA实质上是从 转换波震相 ，可 以获得相应间断面 的图像 。然而 ，当 已知 的观测数据 中捕捉数据 的基本结构 。从线性变 地震 台站位于较厚沉积层地 区时 ，由于沉积层 内形 换和线性空间角度 ，原始源信号可 以看作线性空I司 成 的多次地震反射波振幅较强 ，且又与 间断面的转 的基信号 ，且为相互独立 的非高斯分布量 ，混合信号 换波在 同一时段出现，就会将转换波信息淹没掉 。目 则是基信号的线性组合，ICA就是在源信号和线性 前 ，现有 的技术手段 尚无法将混合 的转换波和多次 变换均不可知情况下 ，从观测 的混合 信号 中估计 出 反射波进行有效 的分离 ，从而导致在沉积 区无法从 数据空 间的基本 结构 。 在沉积 区远震转换波和多次反射分离 问题 中， 接收函数 中提取间断面信息 。本文试 图通过 I