一种不依赖子波的叠后波阻抗反演方法及装置

本说明书涉及计算机，尤其是一种不依赖子波的叠后波阻抗反演方法及装置。

背景技术：

1、叠后地震波阻抗反演在储层预测中起着重要作用。作为地震反演算法所需的输入，地震子波将直接影响反演精度，错误的子波可能导致错误的结果。根据是否使用测井数据，常规的地震子波提取方法可以分为两类。在仅从地震数据中提取子波的情况下，通常假定地震子波是最小相位或零相位，然后通过地震记录自相关方法估计子波。在现有技术中，基于地震和测井数据估计子波是地震分析中的常规做法，通常采用最小二乘拟合和频谱匹配方法，并假定地震数据可以用褶积正演模型获取。井震标定过程涉及使用测井信息正演合成地震记录与井旁道观测地震记录进行比较。通常通过合成数据和观测数据之间的相关系数来衡量井震标定的质量，该过程取决于可靠的测井数据和提取的地震子波。在此基础上，往往还需要子波相位旋转来进一步改善匹配。因此，提取子波的准确性对最终合成地震记录的影响很大。目前在生产中通常手动进行井震标定，整个过程涉及多种解释步骤，存在耗时长和容易引入人为误差的局限。

2、因此，现在亟需一种不依赖子波的叠后波阻抗反演方法，消除实际数据地震子波在目标函数中的影响，从而消除子波和噪声干扰的影响，提高反演结果的稳定性。

技术实现思路

1、为解决现有技术中无法实现高精度地震子波提取，因而会直接影响叠后地震反演精度的问题，本说明书实施例提供了一种不依赖子波的叠后波阻抗反演方法及装置，实现了消除实际数据地震子波在目标函数中的影响，从而消除子波和噪声干扰的影响，提高反演结果的稳定性。

2、本说明书实施例提供了一种不依赖子波的叠后波阻抗反演方法，所述方法包括：

3、采集地下结构图像以及地震子波数据，根据所述地震子波数据获取实际数据；

4、根据地下结构图像以及实际数据构建初始褶积模型并进行反演；

5、所述初始褶积模型根据波阻抗以及地震子波数据得到合成数据，根据所述合成数据以及实际数据重构所述初始褶积模型生成褶积正演模型；

6、获取井旁道反射系数以及井旁道地震记录，根据所述井旁道反射系数、井旁道地震记录以及所述实际数据消除所述褶积正演模型中地震子波的影响，生成不依赖子波的褶积模型；

7、根据所述初始褶积模型的反演结果以及所述不依赖子波的褶积模型构建不依赖子波的叠后波阻抗反演目标函数。

8、进一步的，采集地下结构图像以及地震子波数据，根据所述地震子波数据获取实际数据进一步包括，

9、对所述地震子波数据进行去卷积、多次消除以及数据迁移处理，生成所述实际数据。

10、进一步的，根据地下结构图像以及实际数据构建初始褶积模型并进行反演进一步包括，

11、根据所述地下结构图像的每个地震迹线获取地下反射率和波形卷积结果；

12、构建所述初始褶积模型为：

13、wobs*r＝dobs，

14、其中wobs为实际地震数据子波，r为反射系数序列，dobs为观测地震记录；

15、对所述初始褶积模型进行反演，相应的反演结果为：

16、

17、其中，表

18、示向量范数的平方。

19、进一步的，根据所述合成数据以及实际数据重构所述初始褶积模型生成褶积正演模型进一步包括，

20、根据所述合成数据以及实际数据重构所述初始褶积模型生成褶积正演模型为：

21、wsyn*wobs*r＝wsyn*dobs，

22、其中，wsyn为反演所用所述合成数据。

23、进一步的，根据所述井旁道反射系数以及所述实际数据消除所述褶积正演模型中地震子波的影响进一步包括，

24、根据所述井旁道反射系数以及所述实际数据获取井旁道地震记录；

25、根据所述井旁道地震记录消除实际数据的子波影响，表示为，

26、wsyn*wobs*r*rwell＝wsyn*dobs*rwell，

27、其中，rwell为测井数据反射系数，消除所述实际数据以及测井反射系数，得到不依赖子波的褶积模型为，

28、dwell*wsyn*r＝dsyn_well*dobs，

29、其中，dsyn_well为利用测井反射系数与所述合成数据褶积得到的合成井旁道地震记录。

30、进一步的，根据所述褶积模型的反演结果以及所述不依赖子波的褶积模型构建不依赖子波的叠后波阻抗反演目标函数进一步包括，

31、根据所述褶积模型的反演结果以及所述不依赖子波的褶积模型，构建不依赖子波的叠后波阻抗反演目标函数为：

32、

33、另一方面，本说明书实施例还提供了一种不依赖子波的叠后波阻抗反演装置，＝所述装置包括：

34、实际数据获取单元，用于采集地下结构图像以及地震子波数据，根据所述地震子波数据获取实际数据；

35、初始褶积模型构建单元，用于根据地下结构图像以及实际数据构建初始褶积模型并进行反演；

36、褶积正演模型构建单元，用于所述初始褶积模型根据波阻抗以及地震子波数据得到合成数据，根据所述合成数据以及实际数据重构所述初始褶积模型生成褶积正演模型；

37、褶积模型获取单元，用于获取井旁道反射系数以及井旁道地震记录，根据所述井旁道反射系数、井旁道地震记录以及所述实际数据消除所述褶积正演模型中地震子波的影响，生成不依赖子波的褶积模型；

38、目标函数获取单元，用于根据所述初始褶积模型的反演结果以及所述不依赖子波的褶积模型构建不依赖子波的叠后波阻抗反演目标函数。

39、另一方面，本说明书实施例还提供了一种计算机设备，包括存储器、处理器、以及存储在存储器上的计算机程序，处理器执行所述计算机程序时实现上述的方法。

40、另一方面，本说明书实施例还提供了一种计算机存储介质，其上存储有计算机程序，计算机程序被计算机设备的处理器运行时，执行上述的方法。

41、最后，本说明书还提供一种计算机程序产品，所述计算机程序产品被计算机设备的处理器运行时，执行根据上述所述方法的指令。

42、利用本说明书实施例，通过采集地下结构图像以及地震子波数据获取到实际数据，确定了反演用的地震记录；根据地下结构图像以及实际数据构建初始褶积模型并进行反演；根据波阻抗以及地震子波数据得到合成数据，并重构初始褶积模型生成褶积正演模型，获取不依赖子波目标函数中的合成数据，启到带通滤波的作用，从而压制了地震记录中的随机噪声；获取井旁道反射系数以及井旁道地震记录，消除褶积正演模型中地震子波的影响，从而生成不依赖子波的褶积模型，然后构建不依赖子波的叠后波阻抗反演目标函数，根据反演得到的波阻抗结果与测井位置处的井曲线对比检验反演结果的可靠性。实现了消除实际数据地震子波在目标函数中的影响，从而消除子波和噪声干扰的影响，提高反演结果的稳定性。解决了现有技术中无法实现高精度地震子波提取和叠后地震波阻抗反演的技术问题。