买专利卖专利找龙图腾,真高效! 查专利查商标用IPTOP,全免费!专利年费监控用IP管家,真方便!
申请/专利权人:哈尔滨工程大学
摘要:本发明涉及一种基于Kaiser时窗积分的OBN资料全波形反演方法,属于地震波形反演领域,所述方法为首先采用Kaiser时窗对震源子波和OBN观测地震记录分别进行积分,然后采用声波方程有限差分正演模拟技术得到正传波场和模拟地震记录,将该模拟地震记录与积分后的观测地震记录的残差记录平方和作为反演的目标函数,并采用声波方程有限差分逆时延拓技术得到反传波场,通过正传波场与反传波场的互相关得到目标函数的梯度,最后对速度模型进行多次迭代更新,从而实现OBN资料全波形反演。本发明方法可有效压制积分处理中的高频旁瓣;且在缺失2Hz以下低频信息的情况下,可为后续的常规全波形反演提供可靠的初始模型,有效避免反演算法陷入局部极值。
主权项:1.一种基于Kaiser时窗积分的OBN资料全波形反演方法,其特征在于,所述方法具体包括以下步骤:1基于Marmousi速度模型,采用海底地震节点的采集方式进行地震数据采集,海底地震节点简称OBN,炮点位于海面,检波点位于海底,然后采用声波方程有限差分正演模拟技术,以震源子波f为扰动,模拟得到OBN观测地震记录Pobs;2采用Kaiser时窗分别对震源子波f和OBN观测地震记录Pobs进行积分,积分次数为n,积分时窗长度为2Tw,得到积分后的震源子波和OBN观测地震记录震源子波的整个积分处理为如下的褶积过程: 式2中,Tw为时窗长度的一半,符号“*”表示褶积,ht为Kaiser时窗序列,其表达式为: 式3中,I0为第一类贝瑟尔函数,β为形状因子,OBN观测地震记录的积分处理过程与震源子波相同;3在全波形反演的每次迭代中,基于初始速度模型,采用声波方程有限差分正演模拟技术,得到以积分后的震源子波为扰动的正传波场和模拟地震记录4建立L2范数目标函数E,其是由模拟地震记录与积分后的OBN观测地震记录的残差记录平方和计算得到;5若目标函数满足给定门限,则输出反演的速度模型,反之,则进行步骤5~9;6采用声波方程有限差分逆时延拓技术得到以残差记录为扰动的反传波场7将正传波场与反传波场进行互相关得到目标函数的梯度g;8计算迭代步长,并采用共轭梯度法对速度模型进行更新;9若达到最大迭代次数,则输出最终反演的速度模型,反之,则重复进行步骤3~9。
全文数据:
权利要求:
百度查询: 哈尔滨工程大学 一种基于Kaiser时窗积分的OBN资料全波形反演方法
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。