| 97 | 0 | 298 |
| 下载次数 | 被引频次 | 阅读次数 |
单输入多输出系统提取地震子波的方法存在要求地震道之间的反射系数相同,子波多个且相异(即假设地震记录的差异全部由子波造成,与反射系数无关),以及对噪声敏感等缺陷。为了克服这种假设,本文利用地震数据的频率特征,运用二进正交小波变换Mallat算法的思想,讨论了两种构建伪单输入多输出系统的方法,推导了该系统下噪声子空间提取地震子波的算法,并进行了模型测试和实际资料验证工作。测试和验证结果显示该方法具有较强的抗噪能力。
Abstract:The method of extracting seismic wavelet based on the Single-input Multi-output(SIMO)System,requires reflection coefficients to be the same as the others and wavelets to be different from the others.In other words,it unreasonably assumes that all the difference of seismic record is caused by multiwavelet and has nothing to do with reflection coefficients.In order to avoid these assumptions,Aiming at the character of seismic signal and using the thinking about Mallat algorithm of the Discrete Wavelet Transform,this paper makes a preliminary discussion about two methods of building Pseudo Single-Input Multi-Output System(PSIMO),and then deduces the algorithm of extracting seismic wavelet with noise subspace method in PSIMO system.The results of simulation and application demonstrate that those methods are effective and characterized by strong stability to resist noise.
[1]戴永寿,王俊岭,王伟伟,等.基于高阶累积量ARMA模型线性非线性结合的地震子波提取方法研究[J].地球物理学报,2008,51(6):1851-1859.DAI Yong-Shou,WANG Jun-Ling,WANG Wei-Wei,et al.Seismic wavelet extraction via cumulant-based ARMA model approach with linear and nonlinear combination[J].Chinese J Geophys(in Chinese),2008,51(6):1851-1859.
[2]杨培杰,潘勇,穆星,等.子空间法单输入多输出系统混合相位地震子波提取[J].中国石油大学学报(自然版),2010,34(1):41-45.YANG Pei-Jie,PAN Yong,MU Xing,et al.Mix-phase seismic wavelet extraction of SIMO system by subspace method[J].Journal of China University of Petroleum,2010,34(1):41-45.
[3]李军红,张建中,黄忠来,等.基于倒双谱的地震子波估计方法[J].中国海洋大学学报(自然版),2014,44(11):67-73.LI Jun-Hong,ZHANG Jian-Zhong,HUANG Zhong-Lai,et al.A method for seismic wavelet estimation based on the bicepstrum of seismic data[J].Periodical of Ocean University of China,2014,44(11):67-73.
[4]乐友喜,王才经.任意相位子波的分形反褶积方法[J].石油地球物理勘探,1996,31(6):826-834.Le Y X,Wang C J.Fractal deconvolution of any-phase wavelets[J].Oil Geophysical Prospecting,1996,31(6):826-834.
[5]赵秋亮,李录明,罗省贤.基于分形方法的地震子波提取及应用[J].石油物探,2005,44(1):7-11.ZHANG Qiu-Liang,LI Lu-Ming,LUO Sheng-Xian.Seismic wavelet extraction and application based on fractal technique[J].Geophysical Prospecting for Petroleum,2005,44(1):7-11.
[6]张贤达,保铮.通讯信号处理[M].北京:国防工业出版社,2000.ZHANG Xian-Da,BAO Zheng.Basis for communication precess[M].BeiJing:National Defence of Industry Press,2000.
[7]马建仓.盲信号处理[M].北京:国防工业出版社,2006.MA Jian-Cang,Blind Signal Processing[M].Beijing:National defence of Industry Press,2006.
[8]Zhang X D,Zhang Y S.FIR system identification using high order cumulants alone.IEEE Trans on SP,1994,42(10):2854-2858.
[9]Ulrych T J,Velis D R,Sacchi M D.Wavelet estimation revisited[J].The Leading Edge,1995,14(11):1139-1143.
[10]李国发,彭苏萍,高日胜,等.复赛谱域提取混合相位子波的方法[J].天然气工业,2005,25(1):85-87.LI Guo-Fang,PENG Su-Ping,GAO Ri-Sheng,et al.Method of mixed phase wavelet extraction in cepstrum domain[J].Natural Gas Industry,2005,25(1):85-87.
[11]Petropulu A P et al.Phase reconstruvtion from bispectrum slices[J].IEEE Trans on SP,1998,46(2):527-530.
[12]Pan R L et al.Phase reconstruction in the trispectrum domain[J].IEEETrans on Acoustics,Speech,andSP,1987,36(6):895-897.
[13]李大卫,尹成,熊晓军,等.高阶谱混合方法地震子波估计及处理[J].地球物理学进展,2005,20(1):29-33.LI Da-Wei,YI Cheng,XIONG Xiao-Jun,et al.A hybrid high order spectrum method for seismic wavelet estimation and processing[J].Progress in Geophysics,2005,20(1):29-33.
[14]谢桂生,石玉梅,魏野.双谱地震子波估计[J].西南石油学院学报,2000,22(3):25-28.XIE Gui-Sheng,SHI Yu-Mei,WEI Ye.Wavelet estimation based on bispectrum[J].Journal of Southwest Petroleum Institute,2000,22(3):25-28.
[15]Misra S,Sacchi M.Wavelet estimation by non-linear optimization of all pass operators[J].CSEG RECORDER.2006,38-42.
[16]Porsani M J,Vrsin B.Mixed-phase deconvolution[J].Geophysics,1998,63(2):637-647.
[17]Ursin B,Porsani M J.Estimation of an optimal mixed-phase inverse filter[J].Geophysical Prospecting,2000,48(4):663-676.
[18]Buland Omre H.Bayesian wavelet estimation from seismic and well data[J].Geophysics,2003,68(6):2000-2009.
[19]张广智,刘洪,印兴耀.井旁道地震子波精细提取方法[J].石油地球物理勘探,2005,40(2):158-162.ZHANG Guang-Zhi,LIU Hong,YIN Xing-Yao.Method for fine picking up seismic wavelet at uphole trace[J].Oil Geophysical Prospecting,2005,40(2):158-162.
[20]Lang Tong.Blind identification and equalization based on secondorder statistics:a time domain approach[J].IEEE Transacticins on Information Theory,1994,40(2):340-349.
[21]蒋静.基于子空间的二阶统计量盲信道辨识算法研究[D].郑州:郑州大学,2010.JIANG Jing.Blind Channel Identification Based on Subspace Method of Second-order Statistics[D].Zhengzhou:Zhengzhou University,2010.
[22]杨培杰.地震子波盲提取与非线性反演[D].北京:中国石油大学,2008.YANG PeiJie.Seismic Wavelet Blind Extraction and Non-Liner Inversion[D].Beijing:China University of Petroleum(East China),2008.
[23]邱娜.地震子波分解与重构技术研究[D].青岛:中国海洋大学,2012.QIU Na.Research on Seismic wavelet Decomposition and Reconstruction Technology[D].Qingdao:Ocean University of China,2012.
[24]Eric Moulines,Pierre Duhamel.Subspace Methods for the Blind Identification of Multichannel FIR Filters[J].IEEE Transactions on Signal Processing,1995,43(2):516-525.
[25]何子苏,夏威.现代数字处理及其应用[M].北京:清华大学出版社,2004.HE Zi-Su,XIA Wei.Modern digital processing and its application[M].Beijing:Tsinghua University Press,2004.
[26]Eric Moulines,Pierre Duhamel.Subspace methods for the blind identification of multichannel FIR filters[J].IEEE Transactions on Signal Processing,1995,43(2):516-525.
[27]Xia W,He Z H,Liu B Y,et al.Adaptive multichannel blind identification using manifold optimization[J].Signal Processing,2008(88):1595-1605.
[28]周伟.基于MATLAB的小波分析应用(第2版)[M].西安:西安电子科技大学出版社出版,2010.ZHOU Wei.Wavelet analysis applications based on MATLAB(second edition)[M].XiAn:XiAn Electronic Science&Technology University Press,2010.
[29]葛哲学,沙威.小波分析理论与MATLAB R2007实现[M].北京:电子工业出版社,2007,10:76-85.GE Zhe-Xue,SHA Wei.Wavelet analysis theory and implementation of R2007MATLAB[M].BeiJing:Electronic Industry Press,2007,10:76-85.
基本信息:
DOI:10.16441/j.cnki.hdxb.20150029
中图分类号:P631.4
引用信息:
[1]张汉羽,刘怀山,邢磊,等.基于伪单输入多输出系统的子波提取研究[J],2016,46(06):74-83.DOI:10.16441/j.cnki.hdxb.20150029.
基金信息:
国家自然科学基金项目(41176077、41230318和41304096);; 教育部博士点基金项目(20130132120014);; 国家高技术研究与发展计划项目(2013AA092501)资助~~