Compaq Visual Fortran与Matlab混合编程实现小波域地震信号多尺度分析
|
摘要
小波函数在时间域和频率域均具有良好的局部化特征,运用小波函数对地震信号进行变换处理,能够提取体现地质体宏观与微观特征的属性参数。在Compaq Visual Fortran环境下,调用Matlab小波工具箱中的函数与命令进行混合编程,能够快速、准确地实现地震信号的多尺度分析,提取体现地质体频率属性的最佳尺度参数。运用Compaq Visual Fortran卓越的计算功能完成较大数据文件的读取,再调用Matlab中丰富的小波函数及其它数学函数库,实现地震信号的各类计算,可有效地减轻程序编写的难度,提高工作效率,缩短研究周期。在含油气的地质体中,尺度参数越大,频率衰减则更加突出,最佳的尺度参数能为油气预测提供重要的参考信息。该方法已在川西坳陷深层须家河组新场地区含气性检测中推广试用,所获得的异常区段与多口测井含气段吻合良好。
As a wavelet function has good localizing characters in time and frequency domain, its application in the transformation of the seismic signal can obtain the attribute parameters being able to embody the macrocosmic and microcosmic characters of the geologic body. In Compaq Visual Fortran language, we design the program by calling the functions and commands of the Matlab wavelet tools, then analyze the multi-scale characters of the seismic signal and obtain the optimal scale parameters of the frequency attributes. The process of program is easier, more efficiency and less time-consuming by using the super calculation power of Compaq Visual Fortran to read the data files, and calling abundance wavelet functions of Matlab. Optimal scale parameters can offer important reference information to predict the gas and oil, the practice of this method to process the seismic data of the deep Xujiahe formation gas reservoirs in west Sichuan depression proves that the abnormal sections are corresponding well to the gassy sections.
As a wavelet function has good localizing characters in time and frequency domain, its application in the transformation of the seismic signal can obtain the attribute parameters being able to embody the macrocosmic and microcosmic characters of the geologic body. In Compaq Visual Fortran language, we design the program by calling the functions and commands of the Matlab wavelet tools, then analyze the multi-scale characters of the seismic signal and obtain the optimal scale parameters of the frequency attributes. The process of program is easier, more efficiency and less time-consuming by using the super calculation power of Compaq Visual Fortran to read the data files, and calling abundance wavelet functions of Matlab. Optimal scale parameters can offer important reference information to predict the gas and oil, the practice of this method to process the seismic data of the deep Xujiahe formation gas reservoirs in west Sichuan depression proves that the abnormal sections are corresponding well to the gassy sections.
引文
[1]葛哲学,陈仲生.Matlab时频分析技术及其应用[M].北京:人民邮电出版社,2006.
[2]郭仕剑,王宝顺,贺志国,等.MATLAB7.X数字信号处理[M].北京:人民邮电出版社,2006.
[3]AVIJITCHAKRABORTY,DAVID OKAYA.Frequency-time decomposition of seismic data using wavelet-based methods[J].Geophysics,1995,60(6):1906.
[4]WANG YANGHUA.Seismic time-frequency spectral decomposition by matching pursuit[J].Geophysics,2007,72(1):13.
[5]WANG B,PANNK.Kirchhoff migration of seismic data compressed by matching pursuit decomposition:66th Annual International Meeting[J].SEG,Expanded Ab-stracts,1996,66:1642.
[6]SCOTT SINGLETON,TURHAN M TANER,SVEN TREITEL.Q Estimation using Gabor-Morlet Joint Time-Frequency Analysis Techniques[M].SEG Annu-al Meeting,2006.
[7]徐伯勋,白旭滨,于常青.地震勘探信息技术提取、分析和预测[M].北京:地质出版社,2001.
[8]SATISH SINHA,PARTHA ROUTH S,PHILL ANNO D,et al.Castagna.Spectral decomposition of seismic data with continuous-wavelet transform[J].Geophys-ics,2005,70(6):19.
[9]刘羽.Matlab与Fortran的混和编程及其应用[J].桂林工学院学报,2005,25(1):76.
[10]XU TIANJI,GAN QIGANG,LI XIANGUI.Research and application of absorptive attenuation filtering tech-nology based on Prony transformation[M].Science Press USA Inc.,Geophysical Solutions for Environ-ment and Engineering,2006.
[11]李显贵,徐天吉,甘其刚.新场气田须家河气藏含气性地震检测研究[J].石油物探,2006,45(2):186.
[12]李显贵,徐天吉,范菊芬.AG技术在XC气田T3x气藏的应用研究[J].物探化探计算技术,2007,29(1):39.
[13]徐天吉,李显贵,范菊芬,等.运用Hartley变换模拟地震波在正交各向异性介质中的传播[J].物探化探计算技术,2006,28(4):299.
[2]郭仕剑,王宝顺,贺志国,等.MATLAB7.X数字信号处理[M].北京:人民邮电出版社,2006.
[3]AVIJITCHAKRABORTY,DAVID OKAYA.Frequency-time decomposition of seismic data using wavelet-based methods[J].Geophysics,1995,60(6):1906.
[4]WANG YANGHUA.Seismic time-frequency spectral decomposition by matching pursuit[J].Geophysics,2007,72(1):13.
[5]WANG B,PANNK.Kirchhoff migration of seismic data compressed by matching pursuit decomposition:66th Annual International Meeting[J].SEG,Expanded Ab-stracts,1996,66:1642.
[6]SCOTT SINGLETON,TURHAN M TANER,SVEN TREITEL.Q Estimation using Gabor-Morlet Joint Time-Frequency Analysis Techniques[M].SEG Annu-al Meeting,2006.
[7]徐伯勋,白旭滨,于常青.地震勘探信息技术提取、分析和预测[M].北京:地质出版社,2001.
[8]SATISH SINHA,PARTHA ROUTH S,PHILL ANNO D,et al.Castagna.Spectral decomposition of seismic data with continuous-wavelet transform[J].Geophys-ics,2005,70(6):19.
[9]刘羽.Matlab与Fortran的混和编程及其应用[J].桂林工学院学报,2005,25(1):76.
[10]XU TIANJI,GAN QIGANG,LI XIANGUI.Research and application of absorptive attenuation filtering tech-nology based on Prony transformation[M].Science Press USA Inc.,Geophysical Solutions for Environ-ment and Engineering,2006.
[11]李显贵,徐天吉,甘其刚.新场气田须家河气藏含气性地震检测研究[J].石油物探,2006,45(2):186.
[12]李显贵,徐天吉,范菊芬.AG技术在XC气田T3x气藏的应用研究[J].物探化探计算技术,2007,29(1):39.
[13]徐天吉,李显贵,范菊芬,等.运用Hartley变换模拟地震波在正交各向异性介质中的传播[J].物探化探计算技术,2006,28(4):299.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心