【技术实现步骤摘要】
本申请属于地震波成像，具体涉及一种基于波数调谐的最小二乘偏移成像方法及装置。
技术介绍
1、经典最小二乘偏移成像通过匹配地震观测数据迭代估计地下反射率模型，然而该方法对速度模型和地震子波精度的高度依赖，同时采用了基于波形残差的数据匹配机制，从而导致其在实际地震资料应用中仍然存在较大挑战。成像域最小二乘偏移方法通过计算点扩散函数，近似hessian的逆算子，无需迭代一步获取最小二乘偏移结果，缓解了算法对高精度速度和子波的要求，在实际地震资料处理中已经展示出良好的应用效果。
2、利用波动方程计算全局分布的点扩散函数时，所有散射点同时激发散射波场，因此在波场传播过程中，无法避免相邻点扩散函数之间的干涉作用，此外，由于成像过程中受背向散射噪音的影响，成像结果中的低波数成分存在能量泄漏。如何改善偏移成像的波数分布特征，进而获取高质量的最小二乘成像结果，是急需解决的技术问题。
技术实现思路
1、基于以上技术问题，本申请提出一种基于波数调谐的最小二乘偏移成像方法及装置。
2、第一方面，本申请提出一种基于波数调谐的最小二乘偏移成像方法，包括：
3、获取地震数据的逆时偏移像p和点扩散函数体k；
4、将逆时偏移像p进行均匀分块，得到多个第一块数据，对点扩散函数体k进行均匀分块，得到多个第二块数据；
5、利用第一块数据中的逆时偏移像以及第二块数据中的点扩散函数体计算得到最小二乘成像结果li；
6、采用预构建的波数域调谐算子g，对最小
7、根据波数域散射波成像结果si，得到空间域波数调谐最小二乘成像结果p′i；
8、对所有块数据的空间域波数调谐最小二乘成像结果p′i进行拼接，得到完整的波数调谐的最小二乘偏移成像结果p′。
9、所述利用第一块数据中的逆时偏移像以及第二块数据中的点扩散函数体计算得到最小二乘成像结果li，包括：
10、利用傅里叶变换将逆时偏移像pi和点扩散函数ki转到频率波数域，对频率波数域的逆时偏移像pi和点扩散函数ki进行空间反褶积，得到最小二乘成像结果li，其中，pi为第i个第一块数据的逆时偏移像，ki为第i个第二块数据的点扩散函数体；
11、根据波数域散射波成像结果si得到空间域波数调谐最小二乘成像结果p′i，包括：根据波数域散射波成像结果si进行空间反傅里叶变换，得到空间域波数调谐最小二乘成像结果p′i。
12、所述预构建的波数域调谐算子g，计算式如下：
13、
14、其中，g为预构建的波数域调谐算子，r为分量因子，r0为设定的衰减波数起点，δ为设定的正则化因子。
15、所述分量因子的计算式如下：
16、
17、其中，kx为空间直角坐标系中x轴方向的波数分量，ky为空间直角坐标系中y轴方向的波数分量，kz为空间直角坐标系中z轴方向的波数分量，所述波数分量由地震数据计算得到。
18、所述波数域散射波成像结果si，计算式如下：
19、si＝li·g
20、其中，g为预构建的波数域调谐算子，li为最小二乘成像结果。
21、第二方面，本申请提出一种基于波数调谐的最小二乘偏移成像装置，包括：
22、数据获取模块，用于获取地震数据的逆时偏移像p和点扩散函数体k；
23、数据分块模块，用于将逆时偏移像p进行均匀分块，得到多个第一块数据，对点扩散函数体k进行均匀分块，得到多个第二块数据；
24、第一结果成像模块，用于利用第一块数据中的逆时偏移像以及第二块数据中的点扩散函数体计算得到最小二乘成像结果li；
25、数据滤波模块，用于采用预构建的波数域调谐算子g，对最小二乘成像结果li进行滤波，得到波数域散射波成像结果si；
26、第二结果成像模块，用于根据波数域散射波成像结果si，得到空间域波数调谐最小二乘成像结果pi′；
27、结果拼接模块，用于对所有块数据的空间域波数调谐最小二乘成像结果pi′进行拼接，得到完整的波数调谐的最小二乘偏移成像结果p′。
28、所述第一结果成像模块，包括：频率转换单元、结果成像单元；
29、所述频率转换单元，用于利用傅里叶变换将逆时偏移像pi和点扩散函数ki转到频率波数域；
30、所述结果成像单元，用于对频率波数域的逆时偏移像pi和点扩散函数ki进行空间反褶积，得到最小二乘成像结果li，其中，pi为第一块数据中第i个数据的逆时偏移像，ki为第二块数据中第i个数据的点扩散函数体；
31、所述第二结果成像模块，用于根据波数域散射波成像结果si进行空间反傅里叶变换，得到空间域波数调谐最小二乘成像结果pi′。
32、在数据滤波模块中，所述预构建的波数域调谐算子g，计算式如下：
33、
34、其中，g为预构建的波数域调谐算子，r为分量因子，r0为设定的衰减波数起点，δ为设定的正则化因子。
35、第三方面，本申请提出一种电子设备，包括：存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行所述的一种基于波数调谐的最小二乘偏移成像方法。
36、第四方面，本申请提出一种计算机可读存储介质，其存储有可执行指令，所述指令当被执行时使得处理器执行所述的一种基于波数调谐的最小二乘偏移成像方法。
37、有益效果：
38、本申请提出一种基于波数调谐的最小二乘偏移成像方法及装置，在最小二乘反演框架下，在分块数据内利用波数域谱调谐算子对偏移成像结果进行优化，改善了最小二乘成像低波数成分，压制背向散射导致的低频泄漏问题，提高了最小二乘成像效果。
本文档来自技高网...【技术保护点】
1.一种基于波数调谐的最小二乘偏移成像方法，其特征在于，包括：
2.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述利用第一块数据中的逆时偏移像以及第二块数据中的点扩散函数体计算得到最小二乘成像结果li，包括：
3.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述预构建的波数域调谐算子g，计算式如下：
4.根据权利要求3所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述分量因子的计算式如下：
5.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述波数域散射波成像结果si，计算式如下：
6.一种基于波数调谐的最小二乘偏移成像装置，其特征在于，包括：
7.根据权利要求6所述的基于波数调谐的最小二乘偏移成像装置，其特征在于，所述第一结果成像模块，包括：频率转换单元、结果成像单元；
8.根据权利要求6所述的基于波数调谐的最小二乘偏移成像装置，其特征在于，在数据滤波模块中，所述预构建的波数域调谐算子g，计算式如下：
9
10.一种计算机可读存储介质，其特征在于，其存储有可执行指令，所述指令当被执行时使得处理器执行权利要求1-5任一项所述的基于波数调谐的最小二乘偏移成像方法。
...【技术特征摘要】
1.一种基于波数调谐的最小二乘偏移成像方法，其特征在于，包括：
2.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述利用第一块数据中的逆时偏移像以及第二块数据中的点扩散函数体计算得到最小二乘成像结果li，包括：
3.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述预构建的波数域调谐算子g，计算式如下：
4.根据权利要求3所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述分量因子的计算式如下：
5.根据权利要求1所述的基于波数调谐的最小二乘偏移成像方法，其特征在于，所述波数域散射波成像结果si，计算式如下：
6.一种基于波数调谐的最小二乘偏移成像装置，...
【专利技术属性】
技术研发人员：孙敏傲，向晨，向坤，白英哲，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。