资源描述：
《基于谱模拟技术的混合相位地震子波估计方法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第40卷第2期石油物探Vol.40,No.22001年6月GEOPHYSICALPROSPECTINGFORPETROLEUMJune.2001基于谱模拟技术的混合相位地震子波估计方法李鲲鹏李衍达张学工许建华(清华大学自动化系,北京100084)李鲲鹏,李衍达,张学工等,基于谱模拟技术混合相位地震子波估计方法,石油物探,2001,40(2):21~28摘要针对现有的地震子波振幅谱模型的不足提出了新的振幅谱模型及其估计方法。介绍了子波的Z变换在单位园上无零点和子波自相关函数已知的条件下混合相位地
2、震子波估计的方法原理和判别准则,并在此基础上提出了基于地震子波振幅谱模拟技术的混合相位地震子波估计方法。合成和实际数据处理结果均表明,基于地震子波振幅谱模拟技术的混合相位地震子波估计方法具有较好的性能。关键词谱模拟混合相位子波估计最小二乘法最大熵LiKunpengetal.Themixed-phasewaveletestimatemethodbasedonspectralmodeling.Gpp,2001,40(2):21~28ABSTRACTThispaperproposedanewamplitudespectrummodelofseismicwaveletandamet
3、hodforitsparameterestimationafterananalysisofdefectforoldmodel.Themethodologyanddeterminationcriterionformixedphaseseismicwaveletestimationwereintroducedundertheconditionsthattherearenozerosontheunitcircleforz-transformationofwaveletandtheautomationofwaveletisassumedtobeknown.Thenthemethodformixe
4、d-phaseseismicwaveletestimationwasworkoutbasedamplitudespectrummodelingfortheseismicwavelet.Syntheticandrealdataexamplesdemonstratedtheeffective-nessofthemethod.Keywords:spectralmodeling,mixedphase,waveletestimating,leastsquare,maximumentropy引言反射波地震勘探的地震记录与地震子波和反射系数序列之间的关系可用一褶积模型表示x(m)=f(m)*r(m)
5、+n(m),(1)式中,f(m)为地震子波,r(m)为反射系数序列,n(m)为噪音。我们感兴趣的是r(m)。地震子波是由于我们无法产生真正的脉冲激励和无法不失真地接收造成的结果,因此估计和消LiKunpeng,DepartmentofAutomation,TsinghuaUniversity,Beijing100084本文于2000年4月4日收到,修改稿于6月22日收到。22石油物探第40卷除地震子波的影响成为地震数据处理流程中所必需的一个步骤。在地震数据被恰当处理(真振幅和无多次波)以及具有充分的声波和密度测井资料的理想情况下可以求出r(m
6、),继而求得地震子波。如果没有测井资料,要得到地震子波就不得不对地震子波和反射系数序列作出某些假设。通常假设地震子波是最小相位的,而反射系数序列是具有零均值的高斯独立同分布随机过程。实际上,地震子波不是最小相位而是混合相位,反射系数序列的功率谱也不是白[1]色而是具有近似分形性,且其功率谱与频率具有近似关系Pr(f)f。这里之所以说它具有近似分形性,是因为反射系数的功率谱具有强烈的振荡而分形信号的功率谱确实是平滑的。为了估计混合相位地震子波,文献[2]提出了用高阶累量进行子波估计的方法。由于二阶累量实际上就是自相关函数,它不包含任何相位信息,而三阶累量对于象地震反射系数序列这样具有对
7、称分布的随机过程又等于零,因而也不包含相位信息。因此,估计地震子波只能用四阶累量,而四阶累量的估计方差太大使估计的子波不稳定,因而需要补充其他的条件对其进行约束。本文在子波的振幅谱是平滑的、反射系数序列的振幅谱是震荡的假设条件下建立了地震子波的振幅谱模型,提出了用最小二乘法拟合地震子波振幅谱的谱模拟方法。并在此基础上提出以拟合的地震子波振幅谱为约束条件,以反褶积后的地震记录最小熵值最大为判别准则估计地震子波的方法。1地震子波的振幅
