本发明涉及地球物理勘探，更具体地，涉及一种四分量地表一致性反褶积方法、装置、电子设备及介质。
背景技术：
1、随着油气勘探技术的发展和勘探开发程度的不断深入，复杂、微幅构造和深层构造的精确成像越来越受到关注。由于多数勘探地区地震地质条件复杂，造成地震资料的信噪比低，分辨率不高，对小断层和小幅度构造以及复杂构造成像不清楚，尤其是岩性储集层，其体积一般都较小，这时常规的地震勘探方法已不能满足这种要求，必须进行高精度勘探，高精度处理也随着高精度采集的实施相应地发展起来。高精度数据处理技术，总的标是追求高信噪比、高分辨率、高保真度，即“三高”处理，因此应对提高信噪比、分辨率、保真度的各种处理方法进行深入透彻的研究，为下一步的精确解释打下良好的基础。这些技术对寻找复杂断块或簿互层、岩性圈闭等隐蔽型油气藏发挥着非常重要的作用，同时，也使勘探风险降低到最小，提高勘探效益。高分辨率是高精度数据处理的基本要求，反褶积技术是提高地震资料分辨率的主要手段之一。
2、反褶积是通过压缩基本地震子波以提高地震资料的时间分辨率的过程。反褶积方法都基于地震波的传播过程是一个线性系统，符合褶积模型，即地震数据是由震源子波和地层反射系数序列的褶积，加上一些随机噪声组成的。反射系数序列本身具有足够的分辨率，只要去掉子波的影响，就能达到提高数据纵向分辨率的目的。
3、常规使用的单道脉冲反褶积和单道预测反褶积为各道用本身的数据作自相关或互相关，以求得本道的反褶积因子，然后用这个反褶积因子和本道数据褶积形成褶积的结果。单道反褶积方法有2个假设前提：一是反射系数序列是白噪化的随机序列；二是输人子波为最小相位。但实际输入的地震道并不完全满足这两个条件由于陆上地震的每个激发点条件(如震源耦合、炸药量大小、地表岩性等)和接收点条件(地表耦合、低速带变化)均不相同，使得相邻地震道的特征产生差异，计算出的反褶积因子变的不稳定。另一方面，随机噪声的影响加重了这种不稳定性。
4、目前尚需开发一种四分量地表一致性反褶积方法。
5、公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
1、本发明提出了一种四分量地表一致性反褶积方法、装置、电子设备及介质，其能够提升地震资料的纵向分辨率，同时保持地震资料的振幅、频率和相位的一致性。
2、第一方面，本公开实施例提供了一种四分量地表一致性反褶积方法，包括：
3、根据地震数据进行一致性谱分析；
4、进行四分量地表一致性谱分解；
5、计算四分量地表一致性反褶积算子；
6、将所述反褶积算子应用于所述地震数据，得到处理后的数据。
7、作为本公开实施例的一种具体实现方式，根据地震数据进行一致性谱分析，进行四分量地表一致性谱分解包括：
8、建立地震记录褶积模型；
9、根据地表一致性假设，计算地震子波；
10、针对所述地震子波进行傅立叶变换，分解为振幅谱和相位谱；
11、针对所述振幅谱取对数，获得对数振幅谱方程，进而获得对数振幅谱方程组。
12、作为本公开实施例的一种具体实现方式，所述地震记录褶积模型为：
13、x(t)＝w(t)*y(t)+n(t)
14、其中，x(t)为地震记录，w(t)为地震子波，y(t)为地层反射系数，n(t)为随机噪声。
15、作为本公开实施例的一种具体实现方式，所述地震子波为：
16、wij(t)＝si(t)*rj(t)*ck(t)*ol(t)
17、其中，wij(t)为地震子波，i为炮点，j为接收点，k＝(i+j)/2，l＝(i-j)/2，si(t)为带有炮点影响的子波分量，rj(t)为带有检波点影响的子波分量，ck(t)为与共中心点有关的子波分量，ol(t)为与偏移距有关的子波分量。
18、作为本公开实施例的一种具体实现方式，所述振幅谱为：
19、
20、其中，aij为振幅谱，为带有炮点影响的振幅谱分量，为带有检波点影响的振幅谱分量，为与共中心点有关的振幅谱分量，为与偏移距有关的振幅谱分量。
21、作为本公开实施例的一种具体实现方式，计算四分量地表一致性反褶积算子包括：
22、通过共轭梯度法求解所述对数振幅谱方程组，计算反褶积因子；
23、对所述反褶积因子取幂、平方后做反付氏变换得到其自相关函数；
24、根据所述自相关函数构造预测滤波方程，求解所述预测滤波方程获得预测滤波算子，进而得到所述反褶积算子。
25、作为本公开实施例的一种具体实现方式，所述对数振幅谱方程组为：
26、
27、
28、
29、
30、其中，i为炮点，j为接收点，k＝(i+j)/2，l＝(i-j)/2，aij为振幅谱，为带有炮点影响的振幅谱分量，为带有检波点影响的振幅谱分量，为与共中心点有关的振幅谱分量，为与偏移距有关的振幅谱分量。
31、第二方面，本公开实施例还提供了一种四分量地表一致性反褶积装置，包括：
32、分析模块，根据地震数据进行一致性谱分析；
33、分解模块，进行四分量地表一致性谱分解；
34、计算模块，计算四分量地表一致性反褶积算子；
35、应用模块，将所述反褶积算子应用于所述地震数据，得到处理后的数据。
36、作为本公开实施例的一种具体实现方式，根据地震数据进行一致性谱分析，进行四分量地表一致性谱分解包括：
37、建立地震记录褶积模型；
38、根据地表一致性假设，计算地震子波；
39、针对所述地震子波进行傅立叶变换，分解为振幅谱和相位谱；
40、针对所述振幅谱取对数，获得对数振幅谱方程，进而获得对数振幅谱方程组。
41、作为本公开实施例的一种具体实现方式，所述地震记录褶积模型为：
42、x(t)＝w(t)*y(t)+n(t)
43、其中，x(t)为地震记录，w(t)为地震子波，y(t)为地层反射系数，n(t)为随机噪声。
44、作为本公开实施例的一种具体实现方式，所述地震子波为：
45、wij(t)＝si(t)*rj(t)*ck(t)*ol(t)
46、其中，wij(t)为地震子波，i为炮点，j为接收点，k＝(i+j)/2，l＝(i-j)/2，si(t)为带有炮点影响的子波分量，rj(t)为带有检波点影响的子波分量，ck(t)为与共中心点有关的子波分量，ol(t)为与偏移距有关的子波分量。
47、作为本公开实施例的一种具体实现方式，所述振幅谱为：
48、
49、其中，aij为振幅谱，为带有炮点影响的振幅谱分量，为带有检波点影响的振幅谱分量，为与共中心点有关的振幅谱分量，为与偏移距有关的振幅谱分量。
50、作为本公开实施例的一种具体实现方式，计算四分量地表一致性反褶积算子包括：
51、通过共轭梯度法求解所述对数振幅谱方程组，计算反褶积因子；
52、对所述反褶积因子取幂、平方后做反付氏变换得到其自相关函数；
53、根据所述自相关函数构造预测滤波方程，求解所述预测滤波方程获得预测滤波算子，进而得到所述反褶积算子。
54、作为本公开实施例的一种具体实现方式，所述对数振幅谱方程组为：
55、
56、
57、
58、
59、其中，i为炮点，j为接收点，k＝(i+j)/2，l＝(i-j)/2，aij为振幅谱，为带有炮点影响的振幅谱分量，为带有检波点影响的振幅谱分量，为与共中心点有关的振幅谱分量，为与偏移距有关的振幅谱分量。
60、第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
61、存储器，存储有可执行指令；
62、处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的四分量地表一致性反褶积方法。
63、第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的四分量地表一致性反褶积方法。
64、其有益效果在于：
65、本发明基于地震子波可以被分解为共炮点、共接收点、共偏移距、等多种成份的思想，消除了炮点、检波点、共深度点和炮检距几个方向上滤波器的混合效应，求出的反褶积因子比较平稳，褶积效果使得地震记录的振幅、频率、相位的一致性好，能更好地压缩地震子波，更利于提高叠加和偏移效果。
66、本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。