PDF(4622 KB)
基于EMD的声波全波列几种时频分析方法对比
宁琴琴,王祝文,于洋,徐方慧,武焕平
PDF(4622 KB)
PDF(4622 KB)
基于EMD的声波全波列几种时频分析方法对比
Comparison of time-frequency analysis methods for acoustic full waveform based on empirical mode decomposition
时频分析方法是处理阵列声波测井数据的有效手段,但不同时频分析方法特性不同,处理效果也存在差异。为此,利用经验模态分解(EMD)将单点的实测声波全波列分解为有限个固有模态函数(IMF)分量,再采用短时傅里叶变换、Choi-Williams分布、平滑伪Wigner-Ville分布及重排谱图处理原始的声波全波列和前3个IMF分量,并从时间分辨率、频率分辨率和交叉项干扰3方面对比这4种方法的优缺点。结果表明:前3个IMF分量表征了声波全波列中不同的模式波,而4种方法均无交叉项干扰影响识别模式波的问题,其中重排谱图的时间分辨率、频率分辨率最高。引入时频分布的边缘特性,利用基于EMD的4种方法处理干层、水层、油水同层、油层井段的声波全波列,结果表明重排谱图使能量聚焦明显,能较清晰地显示能量很低的纵波,波峰时间、主频的定位更精准,对不同流体性质储集层的区分效果最好,进一步验证了重排谱图的优势。
Time-frequency analysis is an effective way to process non-linear and non-stationary signals such as the array acoustic logging data.The array acoustic logging signal is decomposed into finite intrinsic mode functions (IMFs) by empirical mode decomposition (EMD).Four time-frequency analysis methods,namely the short-time Fourier transform,Choi-Williams distribution,smoothed pseudo Wigner-Ville distribution and reassigned spectrum,were used to process the original array acoustic logging data,as well as IMF1,IMF2,and IMF3.These four methods were compared in terms of time resolution,frequency resolution,and cross term interference.The results showed that IMF1,IMF2,and IMF3 represent different components of the original array acoustic logging data.In fact,the four methods did not exhibit significant cross-term interference.Among them,the time and frequency resolution obtained by the reassigned spectrum method were the highest.By introducing the edge characteristics of time-frequency distribution,the four aforementioned methods based on EMD were also used to process the array acoustic logging data of dry layers,water layers,oil-water layers,and oil layers.The results showed that the reassigned spectrum method produced a clearer energy focus and the compressional waves with low energy were also displayed more clearly.Furthermore,the location of the peak time and main frequency was more accurate,and the performance in distinguishing the dry layer,water layer,oil-water layer,and oil layer was the best,thereby verifying the superior performance of the reassigned spectrum method.
阵列声波测井 / EMD / 短时傅里叶变换 / Choi-Williams分布 / 平滑伪Wigner-Ville分布 / 重排谱图 {{custom_keyword}} /
array acoustic logging; / empirical mode decomposition (EMD); / short-time Fourier transform; / Choi-Williams distribution; / smoothing pseudo-Wigner-Ville distribution; / reassigned spectrogram {{custom_keyword}} /
［1］章成广,王冠贵,刘银斌,等.伪瑞利波和斯通利波的特性及横波首至的研究［J］.地球物理测井,1990,14(6):385-391
ZHANG C G,WANG G G,LIU Y B,et al.Investigation on pseudo-Rayleigh wave and stoneley wave and the first arrival of shear wave［J］.Well Logging Technology,1990,14(6):385-391
［2］章成广,江万哲,潘和平.声波测井原理与应用［M］.北京:石油工业出版社,2009:109-126
ZHANG C G,JIANG W Z,PAN H P.Principle and application of sonic logging［M］.Beijing:Petroleum Industry Press,2009:109-126
［3］陈博涛.EMD分解联合时频分析在阵列声波信号中的应用［D］.长春:吉林大学,2010
CHEN B T.Application of the EMD joint time-frequency analysis in array acoustic signal［D］.Changchun:Jilin University,2010
［4］王文文.阵列声波测井信号分析及其流体识别方法研究［D］.成都:西南石油大学,2010
WANG W W.Research on array acoustic logging signal analysis and fluid recognition method［D］.Chengdu:Southwest Petroleum University,2010
［5］钟剑.阵列声波测井属性信息提取方法研究［D］.成都:西南石油大学,2011
ZHONG J.Research on attribute information extraction method of array acoustic logging［D］.Chengdu:Southwest Petroleum University,2011
［6］ZHANG B X,WANG K X.Theoretical study of perturbation method for acoustic multipole logging in anisotropic formation［J］.Journal of the Acoustical Society of America,1996,99(11):2674-2685
［7］葛哲学,陈仲生.Matlab时频分析技术及其应用［M］.北京:人民邮电出版社,2006:1-354
GE Z X,CHEN Z S.MATLAB time-frequency analysis technology and its application［M］.Beijing:Peoples Posts and Telecom Press,2006:1-354
［8］田鑫,章成广,江万哲.短时傅里叶变换在阵列声波测井信息提取中的应用［J］.国外测井技术,2005,20(3):53-55
TIAN X,ZHANG C G,JIANG W Z.Application of short-time Fourier transform in array acoustic logging information extraction［J］.World Well Logging Technology,2005,20(3):53-55
［9］陶钧.基于时频信号分析的阵列声波测井信号处理研究［D］.天津:天津大学,2009
TAO J.The research of array acoustic logging signal process based on the time-frequency signal analysis［D］.Tianjin:Tianjin University,2009
［10］刘娇,朱强,李其斌,等.基于时频分析的声波时差检测方法研究［J］.微计算机信息,2010,26(10):67-69
LIU J,ZHU Q,LI Q B,et al.Research on the acoustic slowness detection method based on the time-frequency analysis［J］.Microcomputer Information,2010,26(10):67-69
［11］王祝文,刘菁华,聂春燕.时频分析的重排方法及其在声波测井信号处理中的应用［J］.吉林大学学报(地球科学版),2007,37(5):1042-1046
WANG Z W,LIU J H,NIE C Y.The reassignment method based on tyime- frequency analysis and its application in the array acoustic logging data processing［J］.Journal of Jilin University (Earth Science Edition),2007,37 (5):1042-1046
［12］王祝文,刘菁华,聂春燕.基于Choi-Williams时频分布的阵列声波测井信号时频分析［J］.地球物理学进展,2007,22(5):1481-1486
WANG Z W,LIU J H,NIE C Y.Time-frequency analysis of array acoustic logging signal based on Choi-Williams energy distribution［J］.Progress in Geophysics,2007,22(5):1481-1486
［13］张学涛,王祝文,原镜海.利用时频分析方法在阵列声波测井中区分油水层［J］.岩性油气藏,2008,21(1):101-104
ZHANG X T,WANG Z W,YUAN J H.Using time-frequency analysis to identify oil and water layers in array acoustic logging［J］.Lithologic Reservoirs,2008,21(1):101-104
［14］高云娇.阵列声波测井信号的Page时频分布特征［D］.长春:吉林大学,2011
GAO Y J.Page time-frequency distribution characteristics of array acoustic logging signals［D］.Changchun:Jilin University,2011
［15］杨闯,王祝文,向旻,等.基于Choi-Williams时频分布的裂缝性地层时频特征［J］.世界地质,2015,34(3):825-829
YANG C,WANG Z W,XIANG M,et al.Time-frequency characteristics of fractured formation based on Choi-Williams energy distribution［J］.Global Geology,2015,34(3):825-829
［16］王祝文,王晓丽,刘菁华,等.裂缝性地层声波测井的联合时频特征［J］.吉林大学学报(地球科学版),2012,42(4):914-920
WANG Z W,WANG X L,LIU J H,et al.Joint time-frequency characteristics of array acoustic logging signal on fractured formation［J］.Journal of Jilin University (Earth Science Edition),2012,42(4):914-920
［17］WANG Z W,WANG X L,XIANG M,et al.Reservoir information extraction using a fractional Fourier transform and a smooth pseudo Wigner-Ville distribution［J］.Applied Geophysics,2012,9(4):391-400
［18］王晓丽.基于分数阶傅里叶变换的地层裂缝识别方法研究——以中国大陆科学钻探(CCSD)为例［D］.长春:吉林大学,2016
WANG X L.Stratum fracture recognition method based on fractional Fourier transform:A case from Chinese Continental Scientific Drilling(CCSD)［D］.Changchun:Jilin University,2016
［19］XIANG M,WANG Z W,LIU J H.Extracting array acoustic logging signal information by combining fractional Fourier transform and Choi-Williams distribution［J］.Applied Acoustics,2015,90(4):111-115
［20］向旻,帕尔哈提,张峰玮.基于两种时频分析的裂缝性地层阵列声波测井信号时频特征［J］.石油地球物理勘探,2018,53(4):849-857
XIANG M,PARHAT,ZHANG F W.Array acoustic loggings characteristics in fractured formations based on the time-frequency domain analysis［J］.Oil Geophysical Prospecting,2018,53(4):849-857
［21］徐方慧,王祝文,刘菁华,等.基于EMD的声波测井信息提取与火成岩裂缝地层特征分析［J］.石油物探,2018,57(6):936-943
XU F H,WANG Z W,LIU J H,et al.Acoustic logging information extraction and Fractural volcanic formation characteristics based on empirical mode decomposition［J］.Geophysical Prospecting for Petroleum,2018,57 (6):936-943
［22］HUANG N E,ZHENG S,STEVEN R,et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear non-stationary time Series Analysis［J］.Proceeding of the Royal Society of London,Series A,1998,45(4):903-995
［23］FLANDRIN P,RILLING G,GONCALVES P.Empirical mode decomposition as a filter bank［J］.IEEE Signal Processing Letters,2004,11(2):112-114
［24］徐晓刚,徐冠雷,王孝通,等.经验模式分解(EMD)及其应用［J］.电子学报,2009,37(3):581-585
XU X G,XU G L,WANG X T,et al.Empirical Mode Decomposition and its application［J］.Acta Electronica Sinica,2009,37(3):581-585
［25］陈娟.Hilbert-Huang变换及其在信号处理中的应用［D］.大连:大连理工大学,2006
CHEN J.Hilbert-Huang transform and its application in signal processing［D］.Dalian:Dalian University of Technology,2006
［26］张贤达,保铮.非平稳信号分析与处理［M］.北京:国防工业出版社,1998:1-437
ZHANG X D,BAO Z.Nonstationary signal analysis and processing［M］.Beijing:National Defense Industry Press,1998:1-437
［27］张东,吴晓琳.短时傅里叶变换在振动信号处理中的应用［J］.计算机与数字工程,2011,39(8):154-155
ZHANG D,WU X L.Vibration signal processing based on short time Fourier transform［J］.Computer and Digital Engineering,2011,39(8):154-155
［28］严海滔,周怀来,牛聪,等.同步挤压改进短时傅里叶变换分频相干技术在断裂识别中的应用［J］.石油地球物理勘探,2019,54(4):860-866
YAN H T,ZHOU H L,NIU C,et al.Fault identification with short-time Fourier transform frequency-decomposed coherence improved by synchronous extrusion［J］.Oil Geophysical Prospecting,2019,54(4):860-866
［29］毛青春,徐分亮.窗函数及其应用［J］.中国水运(学术版),2007(2):230-232
MAO Q C,XU F L.The windows function and its application［J］.China Shipping (Academic Edition),2007 (2):230-232
［30］赵秋芳,党鹏飞,云美厚,等.不同时频变换谱比法品质因子Q估算效果对比分析［J］.地球物理学进展,2018,33(5):2097-2101
ZHAO Q F,DANG P F,YUN M H,et al.Comparative analysis of quality factor Q estimated by spectral ratio method based on different time-frequency transform［J］.Progress in geophysics,2018,33(5):2097-2101
［31］COHEN L.Generalized phase-space distribution functions［J］.Journal of Mathematical Physics,1966,7(5):781-786
［32］CHOI H I,WILLIAMS W J.Improved time-frequency representation of multi-component signals using exponential kernels［J］.IEEE Transactions on Acoustics,Speech,and Signal Processing,1989,37(6):862-871
［33］熊良才,史铁林,杨叔子.Choi-Williams分布参数优化及其应用［J］.华中科技大学学报(自然科学版),2003,31(1):103-104
XIONG L C,SHI T L,YANG S Z.Parameter optimizing of Choi-Williams distribution and its application［J］.Journal of Huazhong University of Science and Technology (Natural Science Edition),2003,31(1):103-104
［34］LIN P Z.An introduction to time-frequency signal analysis［J］.Sensor Review,1997,17(1):46-53
［35］程柏林,马晓岩,陈蓓.基于时频重排多分量辐射源信号分析研究［J］.系统工程与电子技术,2006,28(5):684-686
CHENG B L,MA X Y,CHEN B.Analysis of multi-component radiator signal based on time-frequency reassignment method［J］.Systems Engineering and Electronics,2006,28(5):684-686
［36］郝志华,关榆君,马孝江.基于经验模式分解的时频重排方法及其应用［J］.测试技术学报,2007,21(1):17-22
HAO Z H,GUAN Y J,MA X J.Time-frequency reassignment method based on empirical mode decomposition and the application to fault diagnosis［J］.Journal of Test and Measurement Technology,2007,21(1):17-22
［37］荣锋,陈宁,郭翠娟,等.一种EEMD-HHT与时频重排结合的转速曲线估计新方法［J］.振动与冲击,2018,37(22):81-87
RONG F,CHEN N,GUO C J,et al.New method for speed curve estimation by using the EEMD-HHT combined with time-frequency reassignment［J］.Journal of Vibration and Shock,2018,37(22):81-87
［38］MALLAT S.信号处理的小波导引［M］.杨力华,戴道清,译.北京:机械工业出版社,2002:80-90
MALLAT S.A wavelet tour of signal processing［M］.YANG L H,DAI D Q,translator.Beijing:China Machine Press,2002:80-90
［39］吴小羊,刘天佑.基于时频重排的地震信号Wigner-Ville分布时频分析［J］.石油地球物理勘探,2009,44(2):201-205
WU X Y,LIU T Y.Time-frequency analysis of Wigner-Ville distribution of seismic signal based on time-frequency rearrangement［J］.Oil Geophysical Prospecting,2009,44(2):201-205
国家自然科学基金项目(41874135)资助。
苏公网安备 32011502010797号/
| 〈 |
|
〉 |
