摘要
基于匹配追踪的平滑伪Wigner-Ville分布具有良好的时频分辨率,但由于计算量过大,影响其在实际地震数据处理中的应用。在深入研究基于OpenMP的并行算法的基础上,本文将该并行计算方式应用于匹配追踪法对Morlet子波的寻优过程以及计算多个Morlet子波平滑伪Wigner-Ville分布。实际数据计算结果表明:该并行算法加速性能良好,加速比随线程数增加呈近线性增长,即在确保计算精度的情况下,能显著提高计算效率;而且计算量越大其加速性能体现越充分。
The smoothed-pseudo Wigner-Ville distribution(SPWVD)based on matching pursuit has a high time-frequency resolution,but the calculation is extremely large in the matching pursuit algorithm and SPWVD calculation of all the wavelets.To improve the efficiency of this method,we propose a parallel computing method for the matching pursuit time-frequency analysis based on OpenMP.The searching process in the matching pursuit algorithm and the SPWVD calculation of all the wavelets are achieved in parallel.Field seismic data processing shows that the efficiency of the proposed parallel algorithm is significantly improved,and its computation precision is the same as the serial algorithm.Parallel calculation performance shows that the speedup grows linearly as the number of processors increases,and larger computation amount can fully exploit the advantage of multicore processors.
引文
[1]Mallat S G,Zhang Z.Matching pursuits with timefrequency dictionaries.IEEE Transactions on Signal Processing,1993,41(12):3397-3415.
[2]范兴利,成谷.基于Morlet小波尺度参数寻优的匹配追踪时频分析.中山大学学报自然科学版,2014,53(6):85-92.Fan Xingli,Cheng Gu.Matching pursuit time-frequency analysis based on Morlet wavelet scale parameter optimization.Acta Scientiarum Naturalium Universitatis Sunyatseni,2014,53(6):85-92.
[3]Castagna J P,Sun S,Siegfried R W.Instantaneous spectral analysis:detection of low-frequency shadows associated with hydrocarbons.The Leading Edge,2003,22(2):120-127.
[4]孙万元,张会星,杜艺可.匹配追踪时频分析及其在油气检测中的应用.山东科技大学学报(自然科学版),2011,30(4):51-57.Sun Wanyuan,Zhang Huixing,Du Yike.Matching trace time-frequency analysis and its application in oil gas detection.Journal of Shandong University of Science and Technology(Natural Science),2011,30(4):51-57.
[5]Partyka G,Gridley J,Lopez J.Interpretational applications of spectral decomposition in reservoir characterization.The Leading Edge,1999,18(3):173-184.
[6]邓志文,赵贤正,陈雨红等.自适应波形多道匹配追踪断层识别技术.石油地球物理勘探,2017,52(3):532-537,547.Deng Zhiwen,Zhao Xianzheng,Chen Yuhong et al.Fault identification based on multichannel adaptive waveforms matching pursuit.OGP,2017,52(3):532-537,547.
[7]张在金,张军华,李军等.煤系地层地震强反射剥离方法研究及低频伴影分析.石油地球物理勘探,2016,51(2):376-383.Zhang Zaijin,Zhang Junhua,Li Jun et al.A method for stripping coal seam strong reflection and low-frequency shadow analysis.OGP,2016,51(2):376-383.
[8]段文胜,王鹏,党青宁等.应用匹配追踪傅里叶插值技术实现OVT域连片处理.石油地球物理勘探,2017,52(4):669-677.Duan Wensheng,Wang Peng,Dang Qingning et al.5Ddata regularization based on matching pursuit Fourier interpolation for the OVT domain data merging processing.OGP,2017,52(4):669-677.
[9]Liu J,Wu Y,Han D et al.Time-frequency decomposition based on Ricker wavelet.SEG Technical Program Expanded Abstracts,2004,23:1937-1940.
[10]Liu J,Marfurt K J.Matching pursuit decomposition using Morlet wavelets.SEG Technical Program Expanded Abstracts,2005,24:786-789.
[11]张繁昌,李传辉,印兴耀.基于动态匹配子波库的地震数据快速匹配追踪.石油地球物理勘探,2010,45(5):667-673.Zhang Fanchang,Li Chuanhui,Yin Xingyao.Seismic data fast matching pursuit based on dynamic matching wavelet library.OGP,2010,45(5):667-673.
[12]张繁昌,李传辉.基于正交时频原子的地震信号快速匹配追踪.地球物理学报,2012,55(1):277-283.Zhang Fanchang,Li Chuanhui.Orthogonal time-frequency atom based fast matching pursuit for seismic signal.Chinese Journal of Geophysics,2012,55(1):277-283.
[13]邵君.基于MP的信号稀疏分解算法研究[学位论文].四川成都:西南交通大学,2006.
[14]Wang Y.Seismic time-frequency spectral decomposition by matching pursuit.Geophysics,2007,72(1):V13-V20.
[15]Greg S,Richard B,Yang Xiaoyun.Multicore image processing with OpenMP.IEEE Signal Processing Magazine,2010,27(2):134-138.
[16]邱娜.地震子波分解与重构技术研究[学位论文].山东青岛:中国海洋大学,2012.
[17]Morlet J,Arens G,Fourgeau E et al.Wave propagation and sampling theory partⅡ:Sampling theory and complex waves.Geophysics,1982,47(2):222-236.
[18]Taner M T,Koehler F,Sheriff R E.Complex seismic trace analysis.Geophysics,1979,44(6):1041-1063.
[19]付勋勋,秦启荣,徐峰等.基于Wigne-Ville分布和短时傅立叶变换时频分布计算地震波衰减梯度.新疆石油地质,2012,33(3):353-356.Fu Xunxun,Qin Qirong,Xu Feng et al.Estimation of seismic wave energy attenuation gradient based on WVD and STFT.Xinjiang Petroleum Geology,2012,33(3):353-356.
[20]胡广书.现代信号处理教程(第2版).北京:清华大学出版社,2015.
[21]Auger F,Flandrin P.Improving the readability of timefrequency and time-scale representations by the reassignment method.IEEE Transactions on Signal Processing,1995,43(5):1068-1089.
[22]张显文,韩立国,王宇等.地震信号谱分解匹配追踪快速算法及其应用.石油物探,2010,49(1):1-6.Zhang Xianwen,Han Liguo,Wang Yu et al.Seismic spectral decomposition fast matching pursuit algorithm and its application.GPP,2010,49(1):1-6.
[23]陈永健.OpenMP编译与优化技术研究[学位论文].北京:清华大学,2004.
[24]余涛,朱自强,鲁光银等.基于OpenMP的重力张量并行正演.物探化探计算技术,2013,35(4):446-451.Yu Tao,Zhu Ziqiang,Lu Guangyin et al.Parallel forward modeling of gravity tensor based on OpenMP.Computing Techniques for Geophysical and Geochemical Exploration,2013,35(4):446-451.
[25]Pacheco P S著,邓倩妮译.并行程序设计导论.北京:机械工业出版社,2013.
[26]刘扬,王鹏,杨瑞等.基于OpenMP的遥感影像并行ISODATA聚类研究.计算机工程,2016,42(7):238-243.Liu Yang,Wang Peng,Yang Rui et al.Research on parallel ISODATA clustering for remote sensing image based on OpenMP.Computer Engineering,2016,42(7):238-243.
[27]巫小婷,邓家先.基于OpenMP的压缩感知并行处理算法.计算机应用,2012,32(3):617-619.Wu Xiaoting,Deng Jiaxian.Compressed sensing parallel processing algorithm based on OpenMP.Journal of Computer Applications,2012,32(3):617-619.