一种适应于海量地震数据时频分析的快速算法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
为了克服MATLAB时频工具箱在计算大容量的地震数据的分频处理时,存在计算速度慢、内存不足以及FFT的计算长度严重影响时频分辨率等问题,使其适应于实际海量地震数据的时频分析。对MATLAB时频工具箱中基于离散短时傅里叶变换的时频算法进行了改进,采用Fortran语言在微机上编程实现了该算法。理论实算表明,改进后的Fortran程序能够高效地计算大容量、短时窗地震数据的分频处理。
As issues of low calculation speed,out of memory,and the calculation length of FFT seriously affecting the capacity of time-frequency resolution,etc.,happened during calculating of frequency decomposition on large-capacity seismic data.In order to overcome those issues,this paper improved MATLAB time-frequency toolbox's time-frequency algorithm,which based on the discrete short time Fourier transform,implemented the algorithm by applying of Fortran language in microcomputer programming.The theory calculation tests show that the improved Fortran program can efficiently calculate the frequency decomposition of a large-capacity,short-term window seismic data.
引文
[1]PANERAS D E,MANI R,NAWAB S H.STFT computation usingpruned FFT algorithms[J].IEEE Signal Processing Letters,1994,1(4):61-63.
    [2]黄德济,赵宪生.时频域中主参数剖面的计算及应用[J].物探与化探计算技术,1994,16(3):197-204.
    [3]STOCKWELL R G,MANSINHA L,LOWE R P.Localization of thecomplex spectrum:the S transform[J].IEEE Trans on Signal Pro-cessing,1996,44(4):998-1001.
    [4]徐天吉,程冰洁.瞬时谱分析技术与气水识别方法研究[J].石油天然气学报,2008,30(4):73-78.
    [5]刘喜武,刘婉莹,刘洪,等.地震信号广义时频分析及其数值实现[J].物探与化探计算技术,2007,29(5):386-390.
    [6]MARFURT K J,KIRLIN R L.Narrow-band spectral analysis andthin-bed tuning[J].Geophysics,2001,66(4):1274-1283.
    [7]田仁飞,曹俊兴,方磊.时频分析技术及其在识别含气层中的应用[J].天然气地球科学,2008,19(2):72-75.
    [8]向富强.地震相干分析和时频分析方法及其在储层描述中的应用[D].成都:成都理工大学,2008.
    [9]何坤,周激流,李健,等.STFT与图像压缩的探讨[J].计算机应用研究,2005,22(3):158-163.
    [10]葛哲学,陈仲生.MATLAB时频分析技术及其应用[M].北京:人民邮电出版社,2007.
    [11]王宏禹,陈喆.非平稳随机信号分析与处理[M].2版.北京:国防工业出版社,2008.
    [12]丁玉美.数字信号处理[M].西安:西安电子科技大学出版社,2005.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心