【技术实现步骤摘要】
本专利技术属于地震数据处理，具体涉及一种地震反射波数据外推插值方法、电子设备及存储介质。
技术介绍
1、近年来地震勘探向着更深更复杂的地质条件迈进，为精确定位油气藏，地震算法对地震数据的密度和质量提出了越来越高的要求。但由于实际环境和生产条件的限制，会导致地震数据有所缺失，缺失道的存在通常是由于湖泊、勘探成本有限或建筑物障碍和接收器位置的障碍物造成的。这些数据的缺失，容易产生空间假频，对后续的偏移、层析成像等处理产生不良影响。
2、目前主要有基于变换域的插值方法、基于预测滤波的插值方法、基于稀疏反演的插值方法、基于波动方程的插值方法以及基于深度学习的插值方法。然而上述方法侧重于利用数据在的一次波信息来重建缺失数据，多次波信息通常被当作噪声处理。但多次波通常携带了比一次波更多的信息，可以提供有关地下结构的额外信息，因此越来越多的研究人员开始利用多次波信息，以扩大地下的成像质量和范围。
3、地震干涉技术是近年来一种基于相关型和褶积型互易定理的地震数据处理方法，在初至信号加强、地震波去噪、地震数据插值等方面有着广泛的应用。基于地震干涉的插值方法通过一次波和多次波的相互转换，成功利用多次波实现了缺失数据的重构。传统干涉法是基于相关算法计算不同道之间的时延量，地震子波的负向旁瓣具有很强的相似性，相关计算会导致重构的子波产生两个正向旁瓣。目前主要是基于反褶积法对重构的地震数据进行整形，降低旁瓣的影响。但反褶积法压制旁瓣的效果在很大程度上依赖于子波提取的准确性，在低信噪比数据中，如何准确提取地震子波也是很大的挑战。p>
技术实现思路
1、本专利技术要解决的问题是准确提取地震数据，提出一种地震反射波数据外推插值方法、电子设备及存储介质。
2、为实现上述目的，本专利技术通过以下技术方案实现：
3、一种地震反射波数据外推插值方法，包括如下步骤：
4、s1.采集地震数据，基于炮点s和炮集sn共有的检波点ri选取参考道集；
5、s2.基于步骤s1得到的参考道集通过高阶累积量算法生成虚拟反射道集；
6、s3.将步骤s1得到的参考道集和步骤s2得到的虚拟反射道集进行褶积计算，生成重构的远偏移距反射道集；
7、s4.对步骤s3得到的重构的远偏移距反射道集应用匹配滤波器，生成与原始道集匹配的重构的反射道集。
8、进一步的，步骤s2的具体实现方法为基于步骤s1得到的参考道集重置炮点和检波点的位置，设置检波点ri处为虚拟反射道集的炮点，生成虚拟反射道集，该步骤通过高阶累积量算法实现，表达式为：
9、im[g(rj|ri)v]＝cum{g(rj|sn),g(ri|sn),g(ri|sn)}＝e[g(rj|sn)g(ri|sn)g(ri|sn)]
10、其中，g(rj|ri)v表示重构的虚拟反射道集，g(ri|sn)和g(rj|sn)分别表示由sn激发、ri和rj接收到的反射道集；cum表示累积量计算，e表示均值计算。
11、进一步的，步骤s3的具体实现方法为将一次波转换为多次波，生成重构的远偏移距反射道集g(rj|s)esv，表达式为：
12、g(rj|s)esv＝2ik∫g(rj|ri)vg(ri|s)d2ri
13、其中，i表示虚部。
14、进一步的，步骤s4的具体实现方法为对虚拟道集应用匹配滤波器，基于匹配滤波器，用于压制伪影的影响，生成与原始道集匹配的虚反射道集，表达式为：
15、∈＝|gpred*f-g|+|f|2
16、其中，gpred表示虚拟道集，g表示原始道集，∈表示拟合函数，当∈值最小时，得到最优匹配滤波器f；
17、利用最优匹配滤波器进行匹配，生成与原始道集匹配的重构的反射道集g'，表达式为：
18、gpred*f＝g'。
19、一种电子设备，包括存储器和处理器，存储器存储有计算机程序，所述的处理器执行所述计算机程序时实现所述的一种地震反射波数据外推插值方法的步骤。
20、一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现所述的一种地震反射波数据外推插值方法。
21、本专利技术的有益效果：
22、本专利技术所述的一种地震反射波数据外推插值方法，有效地利用了原始数据数据中的多次波信息，基于多次波和一次波的相互转换重构远偏移距地震数据，丰富了地震数据所包含的信息。
23、本专利技术所述的一种地震反射波数据外推插值方法，将高阶累积量算法应用于干涉插值法中，压制了旁瓣，保证了干涉插值法重构数据的分辨率。
24、本专利技术所述的一种地震反射波数据外推插值方法，对噪声压制有很好的效果，在实现地震数据重构的同时，有效压制了噪声，提高地震数据的信噪比。
本文档来自技高网...【技术保护点】
1.一种地震反射波数据外推插值方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的一种地震反射波数据外推插值方法，其特征在于，步骤S2的具体实现方法为基于步骤S1得到的参考道集重置炮点和检波点的位置，设置检波点Ri处为虚拟反射道集的炮点，生成虚拟反射道集，该步骤通过高阶累积量算法实现，表达式为：
3.根据权利要求2所述的一种地震反射波数据外推插值方法，其特征在于，步骤S3的具体实现方法为将一次波转换为多次波，生成重构的远偏移距反射道集G(Rj|S)esv，表达式为：
4.根据权利要求3所述的一种地震反射波数据外推插值方法，其特征在于，步骤S4的具体实现方法为对虚拟道集应用匹配滤波器，基于匹配滤波器，用于压制伪影的影响，生成与原始道集匹配的虚反射道集，表达式为：
5.一种电子设备，其特征在于，包括存储器和处理器，存储器存储有计算机程序，所述的处理器执行所述计算机程序时实现权利要求1-4任一项所述的一种地震反射波数据外推插值方法的步骤。
6.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理
...【技术特征摘要】
1.一种地震反射波数据外推插值方法，其特征在于，包括如下步骤：
2.根据权利要求1所述的一种地震反射波数据外推插值方法，其特征在于，步骤s2的具体实现方法为基于步骤s1得到的参考道集重置炮点和检波点的位置，设置检波点ri处为虚拟反射道集的炮点，生成虚拟反射道集，该步骤通过高阶累积量算法实现，表达式为：
3.根据权利要求2所述的一种地震反射波数据外推插值方法，其特征在于，步骤s3的具体实现方法为将一次波转换为多次波，生成重构的远偏移距反射道集g(rj|s)esv，表达式为：
4.根据权利要求3...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。