倒谱对终止相位识别的方法.docVIP

利用倒谱辅助识别终止相探讨—以9?21集集地震为例 于洪波 （福建厦门 361003 厦门市地震局） 摘要：通过对倒谱理论的分析论证和终止相产生原理的分析，提出应用倒谱对终止相识别的理论依据，并根据厦门地震遥测中心所记录的9?21集集地震进行实例分析，认为倒谱分析是一种可以用来辅助识别终止相的方法，与传统的运用终止相的走时特征和动力学特征识别方法结合使用，进而提高终止相识别的准确度。 关键词：终止相；倒谱；震源；识别；9?21 JiJi集集地震 0引言 作者曾对终止相应用进行过探讨(于洪波，2010)，但在应用过程中，深感终止相识别的困难。终止相的产生过程如下图1所示，在理论上终止相实际上是一个叠加信号，针对此特征给我们运用数学处理方法提供了可能，对于地震远场位移谱解析式中震源的有限性因子使得谱具有振荡特性这一现象，启发我们运用对波形谱再进行某种谱分析来识别终止相。笔者认为，使用倒频谱（cepstrum）的数学处理可以辅助我们较好的识别终止相，近而较为准确的对震源参数进行分析、计算。 1 方法的引进 倒频谱(中国电子仪器仪表学会信号处理学会，1983)分析是近代信号处理科学的一项新技术，它可以处理复杂频谱图上的周期结构，它包括功率倒频谱分析和复倒频谱分析两种主要形式。它对于分析具有同族谐频或异族谐频、多成分边频等复杂信号，找出功率谱上不易发现的问题，非常有效。 某时间序列对数功率谱的富氏变换称为它的倒谱，记为F(r)。即 (1) 随着t→→，有f(t)←→F()[或Φ(w)] ←→f()。自变量与f(t)的自变量具有相同的量纲，称为同态频率或倒频率。倒谱分析是一种非线性信号分析技术，这种技术适用于回波检测、语音分析、地球物理数据处理范畴的种种问题。 2 识别原理 在终止相产生的理论描述中，在波形上我们确认终止相的位置实际上是一个叠加信号。 图1 震源破裂及P波射线 Fig.1 Rupture of earthquake and P waveform travel route 作者简介：于洪波（1972～）,男，工程师，现从事数字化地震遥测台网监测工作 2.1 破裂传播及P波传播 我们以最简单的单侧破裂方式为例，说明产生终止相产生的机理。地震时，假设有一应力积累单元AB，破裂自A点开始并终止于B（图1）。由观测点S边SA和SB，则AS与断层破裂方向之间的夹角为。设破裂传播速度为v，则由始破裂点A发出的P波与终止点B发出的P波的到时差为： （2） 如果把沿AB和BS的信号作为沿AS射线的同源振动的回声看待，那么就是体波远场位移的初动半周期，在波形上的显示就是终止相，它是否适用于倒谱分析呢?事实上，答案同样是肯定的。下面的推导将给予证明。 2.2 倒谱分析与终止相相关性分析 我们可以通过富氏变换进入到频率域来讨论。地震体波远场位移振幅谱可表示为（傅承义，1985）： |U（Q；）|= （3） 式中，G()是震源时间函数的谱，M0是地震矩，K是辐射图型因子，称为有限性因子，其中 X= （4） 这个式子与前面讨论的终止相的成因半周期非常相似（式2）。功率谱可记为：|U（Q；）|2=B?;由于多普勒效应，体波振幅谱出现了一系列频率值为f的极小点。不难发现，第一个极小点频率(n=1)的倒数正好是初动半周期(式2)。事实上，这正是同一效应在时间和频率两个不同领域的表现。 (n=1,2,3,…) (5) 现在进入到倒谱域观察它的另一种表现。对数功率谱具有形式 lg|U（Q；）|2=lgB+2lg (6) 式中B是与频率无关的量。此时可以看出，在对数功率谱中出现了一衰减的振荡项2lg。这使我们很自然地联想到能否对它再进行一次谱分析呢?对该式再进行富氏变换，即 (7) 求解这个积分表达式（魏富胜，2000；魏富胜，2003）。则单侧破裂方式的倒谱可写成： (8) 式中，C表示一常数，是欧拉常数，且 (9) 对比式(2)与(9)