【技术实现步骤摘要】
本专利技术涉及垂直地震剖面资料处理领域，尤其涉及一种多道反褶积方法、装置、存储介质及电子设备。
技术介绍
1、垂直地震剖面(vertical seismic profile)，即vsp，是一种地震观测方法。垂直地震剖面是相对于地面地震剖面而言的，是在井中观测地震波场，将地震检波器置于井中不同深度来记录地面震源所产生的地震信号。一般情况下，垂直地震剖面的采集方式是震源位于地表，地震检波器位于井中，采集的数据由直达波场和上行波场组成。常规反褶积方法是将直达波场(主要是下行直达波)作为震源子波，在单道基础上求解反褶积算子。然而，常规方法没有考虑井中地震检波器噪音(非白噪)的影响，而且应用反褶积算子后不同道的子波形态差异大。
技术实现思路
1、针对上述问题，本专利技术提供一种多道反褶积方法、装置、存储介质及电子设备，能够有效展宽垂直地震剖面资料的频带并降低环境噪音的影响。
2、第一方面，本专利技术实施例提供一种多道反褶积方法，包括：
3、拾取垂直地震剖面的共炮点道集的直达波旅行时；
4、以共炮点道集中的多个道作为滑动窗口，将当前滑动窗口范围内的各道变换到频率域，基于直达波旅行时进行当前滑动窗口范围的反褶积处理；
5、逐道滑动所述滑动窗口，完成各滑动窗口范围的反褶积处理，直至完成所述共炮点道集中所有道的反褶积处理。
6、在一些实现方式中，本专利技术实施例提供的多道反褶积方法中，所述基于直达波旅行时进行当前滑动窗口范围的反褶积处理，包
7、在所述当前滑动窗口内计算噪音估计和基于直达波旅行时计算信号估计；
8、基于所述噪音估计和所述信号估计计算反褶积算子，并应用到所述当前滑动窗口内的中心道上。
9、在一些实现方式中，本专利技术实施例提供的多道反褶积方法中，在所述滑动窗口内计算噪音估计，包括：
10、在所述滑动窗口内采用第一计算式计算噪音估计，所述第一计算式如下：
11、
12、其中，et(ω)表示噪音估计，所述噪音包括非白噪声；
13、n表示当前滑动窗口内的道数；
14、dm(ω)表示第m道数据。
15、在一些实现方式中，本专利技术实施例提供的多道反褶积方法中，在所述当前滑动窗口内基于直达波旅行时计算信号估计，包括：
16、在所述当前滑动窗口内采用第二计算式，将原始波场拉平以叠加当前滑动窗口内的多道数据，得到信号估计，所述第二计算式如下：
17、
18、其中，表示信号估计；
19、n表示当前滑动窗口内的道数；
20、dm(ω)表示第m道数据；
21、e表示自然常数；
22、复数
23、ω表示时间域的傅里叶变换的圆周频率；
24、tn表示直达波直达旅行时。
25、在一些实现方式中，本专利技术实施例提供的多道反褶积方法中，所述基于所述噪音估计和所述信号估计计算反褶积算子，包括：
26、采用第三计算式计算反褶积算子，所述第三计算式如下：
27、
28、其中，f(ω)表示反褶积算子；
29、表示信号估计；
30、et(ω)表示噪音估计；
31、*表示共轭转置。
32、第二方面，本专利技术实施例提供一种多道反褶积装置，包括：
33、拾取模块，用于拾取垂直地震剖面的共炮点道集的直达波旅行时；
34、反褶积模块，用于以共炮点道集中的多个道作为滑动窗口，将当前滑动窗口范围内的各道变换到频率域，基于直达波旅行时进行当前滑动窗口范围的反褶积处理；逐道滑动所述滑动窗口，完成各滑动窗口范围的反褶积处理，直至完成所述共炮点道集中所有道的反褶积处理。
35、在一些实现方式中，本专利技术实施例提供的多道反褶积装置中，所述反褶积模块进行当前滑动窗口范围的反褶积处理时，用于：
36、在所述当前滑动窗口内计算噪音估计和基于直达波旅行时计算信号估计；
37、基于所述噪音估计和所述信号估计计算反褶积算子，并应用到所述当前滑动窗口内的中心道上。
38、在一些实现方式中，本专利技术实施例提供的多道反褶积装置中，所述反褶积模块采用第三计算式计算反褶积算子，所述第三计算式如下：
39、
40、其中，f(ω)表示反褶积算子；
41、表示信号估计；
42、et(ω)表示噪音估计；
43、*表示共轭转置。
44、第三方面，本专利技术实施例提供一种计算机存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被一个或多个处理器执行时，实现如第一方面所述的方法。
45、第四方面，本专利技术实施例提供一种电子设备，包括存储器和一个或多个处理器，所述存储器上存储有计算机程序，所述计算机程序被所述一个或多个处理器执行时实现如第一方面所述的方法。
46、与现有技术相比，本专利技术的一个或多个实施例至少能够带来如下有益效果：
47、本专利技术提供了一种垂直地震剖面资料处理的技术方案，拾取了垂直地震剖面的共炮点道集的直达波旅行时；以共炮点道集中的多个道作为滑动窗口，将当前滑动窗口范围内的各道变换到频率域，基于直达波旅行时进行当前滑动窗口范围的反褶积处理；逐道滑动该滑动窗口，完成各滑动窗口范围的反褶积处理，直至完成共炮点道集中所有道的反褶积处理。本方案适用于垂直地震剖面多分量或单分量数据处理，能够有效展宽垂直地震剖面资料的频带并降低环境噪音的影响。
本文档来自技高网...【技术保护点】
1.一种多道反褶积方法，其特征在于，包括：
2.根据权利要求1所述的多道反褶积方法，其特征在于，所述基于直达波旅行时进行当前滑动窗口范围的反褶积处理，包括：
3.根据权利要求2所述的多道反褶积方法，其特征在于，在所述滑动窗口内计算噪音估计，包括：
4.根据权利要求2所述的多道反褶积方法，其特征在于，在所述当前滑动窗口内基于直达波旅行时计算信号估计，包括：
5.根据权利要求2所述的多道反褶积方法，其特征在于，所述基于所述噪音估计和所述信号估计计算反褶积算子，包括：
6.一种多道反褶积装置，其特征在于，包括：
7.根据权利要求6所述的多道反褶积装置，其特征在于，所述反褶积模块进行当前滑动窗口范围的反褶积处理时，用于：
8.根据权利要求7所述的多道反褶积装置，其特征在于，所述反褶积模块采用第三计算式计算反褶积算子，所述第三计算式如下：
9.一种计算机存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被一个或多个处理器执行时，实现如权利要求1至5中任一项所述的方法。
10.一种电子设备，其特征在于，包括存储器和一个或多个处理器，所述存储器上存储有计算机程序，所述计算机程序被所述一个或多个处理器执行时实现如权利要求1至5中任一项所述的方法。
...【技术特征摘要】
1.一种多道反褶积方法，其特征在于，包括：
2.根据权利要求1所述的多道反褶积方法，其特征在于，所述基于直达波旅行时进行当前滑动窗口范围的反褶积处理，包括：
3.根据权利要求2所述的多道反褶积方法，其特征在于，在所述滑动窗口内计算噪音估计，包括：
4.根据权利要求2所述的多道反褶积方法，其特征在于，在所述当前滑动窗口内基于直达波旅行时计算信号估计，包括：
5.根据权利要求2所述的多道反褶积方法，其特征在于，所述基于所述噪音估计和所述信号估计计算反褶积算子，包括：
6.一种多道反褶积装置，其特征在于，包括：
7.根据...
【专利技术属性】
技术研发人员：崔树果，吕秋玲，陈占国，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。