一种地震数据初至波的外推方法技术

本发明专利技术专利提出了一种地震数据初至波的外推方法。首先，在炮集的初至波周围设置时间窗；其次，将参考道和目标道的波场在炮域进行反褶积和叠加，形成虚拟波场；再次，将参考道波场与虚拟波场在检波点域进行褶积和叠加，得到超级虚拟初至波；从次，判断是否采用循环迭代；除此之外，进行初至波外推，目标道位于目标炮集的检波器排列之外；然后，依次在炮域和检波点域进行反褶积叠加及褶积叠加得到外推的超级虚拟初至波；之后，判断是否采用循环迭代。本发明专利技术能有效增强初至波能量，提高初至波信噪比，同时通过初至波外推大大增加可用初至波的数据量，为后续地震数据处理提供更丰富的信息，并能很好地用于复杂近地表条件下采集的地震数据。

【技术实现步骤摘要】
一种地震数据初至波的外推方法
本专利技术涉及地球物理勘探地震数据处理
，特别涉及一种抗噪性强的地震数据初至波外推方法。
技术介绍
地震波速度参数是地震勘探中一个极为重要的物性参数，在地震勘探的各个环节中都发挥着至关重要的作用，不仅为地震数据处理和解释提供主要依据，而且它能够间接的反映地下介质构造和岩性信息。因此，如何精确求取地下介质的速度信息一直是地震勘探领域在不断深入研究的热点问题。初至旅行时层析反演是一种稳定性好且准确高效的速度建模方法，其基于射线理论，利用地震初至波的旅行时信息反演地下介质的速度及界面形态，其反演结果可用于分析不同的岩性分界面，显示地下异常体的位置和规模，识别陡坎、断裂构造等。然而，初至旅行时层析反演的精度与初至拾取质量和射线密度有直接关系，而初至拾取与地震数据的信噪比、各地震道之间的波形变化、近地表条件等有密切关系，在复杂地形条件下，常规初至拾取方法难以获得高精度的初至拾取结果以及数目丰富的高质量初至信息。复杂近地表问题一直是地震勘探所面临的一大难题，如山地、沙漠、戈壁、山前带、黄土塬等地区存在地表本文档来自技高网...

【技术保护点】
1.一种地震数据初至波外推方法，其特征在于，所述方法包括：/n(1)、在原始炮集的初至波周围设置时间窗口，截取包含初至波的波场；/n(2)、对于震源位于X的目标炮，取参考道G(A|X)和目标道G(B|X)，G(A|X)代表震源位于X且在检波器A记录的波场，在傅里叶域，忽略掉振幅项，可表示为

【技术特征摘要】
1.一种地震数据初至波外推方法，其特征在于，所述方法包括：
(1)、在原始炮集的初至波周围设置时间窗口，截取包含初至波的波场；
(2)、对于震源位于X的目标炮，取参考道G(A|X)和目标道G(B|X)，G(A|X)代表震源位于X且在检波器A记录的波场，在傅里叶域，忽略掉振幅项，可表示为其中，ω是角频率，τXA代表初至波旅行时，类似地，G(B|X)代表震源位于X且在检波器B记录的波场，
(3)、在炮域内，将波场G(A|X)和波场G(B|X)进行反褶积，并将每个稳相震源处的反褶积结果叠加，形成检波器A与B之间的虚拟波场G(B|A)vir，在远场近似条件下，计算公式为：



其中，k是波数，Ns是对于检波器对儿A和B而言的稳相震源数目，*代表复共轭，ε为正则化参数；
(4)、在检波点域内，将波场G(A|X)与虚拟波场G(B|A)vir进行褶积，并将每个稳相检波器处的褶积结果叠加，得到震源在X处且在检波器B重构的超级虚拟初至波，计算公式为：



其中，Ng是稳相检波点数目；
(5)、根据超级虚拟初至波的质量决定是否采用循环迭代处理，如果初至波仍难以拾取，使用步骤(4)的输出作为新的输入，并重复步骤(3)和步骤(4)，直到超级初至波的信噪比可被接受；
(6)、对目标炮集所对应的检波器排列中的所有道集都重复步骤(2)到步骤(5)，重构出目标炮集的超级虚拟初至波；
(7)、对震源位于X的目标炮集的初至波进行外推，选取震源位于Y的某炮，取参考道G(A|Y)和目标道G(C|Y)，其中，G(A|Y)的检波器位置A位于目标炮集的检波器排列之内，而G(C|Y)的检波器位置...

【专利技术属性】
技术研发人员：许卓，张凤蛟，韩立国，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22