本发明涉及地震波成像,特别地涉及一种成像方法、装置、存储介质和电子设备。
背景技术:
1、小尺度缝洞精细成像是西北超深层高精度成像的核心问题。最小二乘逆时偏移方法通过估计hessian逆算子可以直接对常规偏移结果进行处理,解决图像照明不足、分辨率低等问题。然而,常规的最小二乘成像技术很难突出小尺度缝洞等小尺度异常体。
2、由此可见,如何对小尺度缝洞等小尺度异常体进行精细成像是本领域亟待解决的技术问题。
技术实现思路
1、本发明提供一种成像方法、装置、存储介质和电子设备,解决了一些技术方案中难以进行小尺度缝洞等小尺度异常体的精细成像的技术问题。
2、第一方面,本发明提供了一种成像方法,包括:
3、对逆时偏移像和全局空变散射函数进行区域分块,以得到分块逆时偏移像和分块全局空变散射函数;
4、将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果;
5、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果。
6、在一些实施例中,将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果,包括:
7、通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果。
8、在一些实施例中,通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果,表达式为:
9、
10、其中,i为分块的序号,li为序号为i的分块的最小二乘成像结果,pi为对应分块的逆时偏移像,ki为对应分块的全局空变散射函数,α为正则化因子。
11、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果,包括:
12、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果;
13、对各个分块的散射波最小二乘成像结果进行拼接,得到最终的散射波最小二乘成像结果。
14、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果,表达式为:
15、si=li·g
16、其中,i为分块的序号,si为序号为i的分块的波数域散射波最小二乘成像结果,li为对应分块的逆时偏移像,g为散射傅里叶滤波算子。
17、在一些实施例中,散射傅里叶滤波算子的表达式为:
18、
19、其中,g为散射傅里叶滤波算子,θ为空间波数分布对应的散射角,θ0为成像点散射张角,r为衰减平滑函数。
20、在一些实施例中,空间波数分布对应的散射角的表达式为:
21、
22、其中,β是阻尼因子。
23、第二方面,本发明提供了一种成像装置,包括:
24、分块模块,用于对逆时偏移像和全局空变散射函数进行区域分块,以得到分块逆时偏移像和分块全局空变散射函数;
25、反褶积模块,用于将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果;
26、滤波模块,用于基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果。
27、在一些实施例中,将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果,包括:
28、通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果。
29、在一些实施例中,通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果,表达式为:
30、
31、其中,i为分块的序号,li为序号为i的分块的最小二乘成像结果,pi为对应分块的逆时偏移像,ki为对应分块的全局空变散射函数,α为正则化因子。
32、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果,包括:
33、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果;
34、对各个分块的散射波最小二乘成像结果进行拼接,得到最终的散射波最小二乘成像结果。
35、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果,表达式为:
36、si=li·g
37、其中,i为分块的序号,si为序号为i的分块的波数域散射波最小二乘成像结果,li为对应分块的逆时偏移像,g为散射傅里叶滤波算子。
38、在一些实施例中,散射傅里叶滤波算子的表达式为:
39、
40、其中,g为散射傅里叶滤波算子,θ为空间波数分布对应的散射角,θ0为成像点散射张角,r为衰减平滑函数。
41、在一些实施例中,空间波数分布对应的散射角的表达式为:
42、
43、其中,β是阻尼因子。
44、第三方面,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现如下方法:
45、对逆时偏移像和全局空变散射函数进行区域分块,以得到分块逆时偏移像和分块全局空变散射函数;
46、将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果;
47、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果。
48、在一些实施例中,将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果,包括:
49、通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果。
50、在一些实施例中,通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果,表达式为:
51、
52、其中,i为分块的序号,li为序号为i的分块的最小二乘成像结果,pi为对应分块的逆时偏移像,ki为对应分块的全局空变散射函数,α为正则化因子。
53、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果,包括:
54、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果;
55、对各个分块的散射波最小二乘成像结果进行拼接,得到最终的散射波最小二乘成像结果。
56、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果,表达式为:
57、si=li·g
58、其中,i为分块的序号,si为序号为i的分块的波数域散射波最小二乘成像结果,li为对应分块的逆时偏移像,g为散射傅里叶滤波算子。
59、在一些实施例中,散射傅里叶滤波算子的表达式为:
60、
61、其中,g为散射傅里叶滤波算子,θ为空间波数分布对应的散射角,θ0为成像点散射张角,r为衰减平滑函数。
62、在一些实施例中,空间波数分布对应的散射角的表达式为:
63、
64、其中,β是阻尼因子。
65、第四方面,本发明提供了一种电子设备,包括处理器和存储器,存储器上存储有计算机程序,处理器执行计算机程序时实现如下方法:
66、对逆时偏移像和全局空变散射函数进行区域分块,以得到分块逆时偏移像和分块全局空变散射函数;
67、将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果;
68、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果。
69、在一些实施例中,将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果,包括:
70、通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果。
71、在一些实施例中,通过傅里叶变换将每一分块逆时偏移像和对应的分块全局空变散射函数转到频率波数域,进行空间反褶积,以得到对应分块的最小二乘成像结果,表达式为:
72、
73、其中,i为分块的序号,li为序号为i的分块的最小二乘成像结果,pi为对应分块的逆时偏移像,ki为对应分块的全局空变散射函数,α为正则化因子。
74、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果,包括:
75、基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果;
76、对各个分块的散射波最小二乘成像结果进行拼接,得到最终的散射波最小二乘成像结果。
77、在一些实施例中,基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,以得到各个分块的散射波最小二乘成像结果,表达式为:
78、si=li·g
79、其中,i为分块的序号,si为序号为i的分块的波数域散射波最小二乘成像结果,li为对应分块的逆时偏移像,g为散射傅里叶滤波算子。
80、在一些实施例中,散射傅里叶滤波算子的表达式为:
81、
82、其中,g为散射傅里叶滤波算子,θ为空间波数分布对应的散射角,θ0为成像点散射张角,r为衰减平滑函数。
83、在一些实施例中,空间波数分布对应的散射角的表达式为:
84、
85、其中,β是阻尼因子。
86、本发明提供的一种成像方法、装置、存储介质和电子设备,通过对逆时偏移像和全局空变散射函数进行区域分块,以得到分块逆时偏移像和分块全局空变散射函数;将每一分块逆时偏移像和对应的分块全局空变散射函数进行空间反褶积,以得到对应分块的最小二乘成像结果;基于散射傅里叶滤波算子对各个分块的最小二乘成像结果进行滤波,得到最终的散射波最小二乘成像结果;针对最小二乘成像无法针对小尺度缝洞进行精细成像的关键技术问题,本发明提供一种散射波最小二乘成像方法,在最小二乘反演成像框架下,利用傅里叶散射波滤波分离,实现散射波最小二乘成像,获取散射异常体的高精度成像,提高小尺度缝洞储层的成像效果;解决了一些技术方案中难以进行小尺度缝洞等小尺度异常体的精细成像的技术问题。