时间-频率域的极化分析及其应用研究.pdfVIP

中国地球物理2011 ·553· (17)油气田与煤田地球物理勘探 时间一频率域的极化分析及其应用研究 of Domain Polarization ApplicationTime-Frequency Analysis 王焕军’林年添郝继杰孙剑 刘建龙 LinNiantianeta1． WangHuanjun 山东科技大学 山东青岛266510 1．引言 随着多分量地震勘探在油气勘探中的应用越来越广泛，多分量信号的分析研究受到应有的重视。在多分量 地震波场分离的方法中，极化分析是一种有效的方法。目前大部分极化分析的方法只在时间域或者只在频率域 进行，比如协方差矩阵法、频率域的奇异值分解法等。时间域的极化滤波方法并不能处理同一时间内不同频率 成分叠加的情况，同时，频率域的滤波方法也不能处理同一个频带内不同时间的信号。综合利用时间域和频率 域的地震波的特征，可以减少以上问题。本文选取时间—频率极化分析进行研究，并分析其应用效果。 2．基本原理 首先考虑正交分量的傅里叶变换。纵波勘探中，对信号进行傅里叶分解可以得到具有固定频率的时间正弦 曲线的叠加，在三分量信号中，分别对垂直分量、径向分量和切向分量进行傅里叶变换。每个频率含有三 个正交的正弦曲线，在三维空间内，整个信号可以看作是椭圆的叠加。 在固定时窗内(设有N个采样点)，傅里叶变换的公式为： ．v—l F一彳(七)一∑a(j)e(2aJk州 j-o 式中，口∽，彳(国表示傅里叶变换前后地震数据的振幅值。傅里叶变换的振幅和相位频谱是正弦曲线的振幅 和相位。在本文所讲的时间一频率域极化分析方法中，上式中传统的傅里叶变换的振幅谱和相位谱要替换 为椭圆元素的新的傅里叶频谱。通过变换时窗的利用，使得此方法同时依赖于频率域和时间域，避免了完 全在时间域或完全在频率域进行的缺点。 要实现时间一频率域表示(TFRs)极化分析方法，基于时窗的傅里叶变换并不是唯一的方法，也可用 小波变换、S变换(StockweU。1996)来实现。 3．模型数据试算 为了验证本文所述方法的有效性，建立了一个模型进行试算，通过模型生成人工合成地震记录。所建模型 生成的地震记录必须包含多种具有不同极化特征的波型。因为Love波和瑞雷波的相位具有明显不同的极化特 征，在人工合成的地震记录中加入Love波和瑞雷面波。运用本文所述的TFR方法，选取一个时窗，对地震记 录的三个分量(垂向、径向、切向)进行傅里叶分解。处理以后的信号的振幅谱、相位谱等有着不同的特征， 尤其是Love波和瑞雷波呈现出明显的分散效应。结果表明，时间—频率域的极化分析方法具有很好的分离效果。 4．总结 本文采用的时间一频率域极化分析方法，综合利用时间域和频率域内的信息，对信号进行时间一频率 表示，能够解决单一的时间域内滤波和单一的频率域内滤波不能处理的问题。此方法通过模型数据的试算 也得到了验证。由于实现这种滤波方法所用的时窗傅里叶变换、S变换、小波变换方法等，本身也具有局 限性，可以根据目标信号的具体极化特征来设计适应某种信号的极化滤波器，以达到解决实际问题的目的， 更好地促进多波多分量地震在油气勘探中的应用。 本研究由国家自然科学基金项目40839901)资助。 参考文献 C and 【1 Pinnegar polarization fJ】．SEG Abstract。2006，25：1223—1227． J R．Time-frequencyanalysisfilter