断裂精细切片计算方法、装置、计算机设备和存储介质与流程

本发明涉及地球物理勘探，特别涉及一种断裂精细切片计算方法、装置、计算机设备和存储介质。

背景技术：

1、地震数据中包含大量的地震波传播形态、运动特征的信息，如何利用其中信息揭示地震数据体中的沉积、岩性和储层的特征，是地球物理研究中的重点与难点。国内外学者在如何运用地震属性方法对这些特征进行描述，以开展断裂和裂缝预测研究的方面，已经取得了大量的成果。目前，常用的地震属性有：相干、曲率、蚂蚁追踪、最大似然、afe、边缘检测、张量等。其中相干体技术有利于识别岩层横向非均一性，对断裂和裂缝的预测十分有效。

2、目前运用相干体技术来进行断裂的刻画主要分为c1、c2和c3三代相干算法：(1)第一代相干(c1)算法，由bahorich和farmer(1995)提出，计算量较小且易于实现，但是局限于信噪比的高低；(2)第二代相干(c2)算法，marfurt和kirlin(1998)改善了信噪比问题，但是更易识别大断层，小断层及裂缝识别效果不好；(3)第三代相干(c3)算法，gersztenkom和marfurt(1999)提出，是基于本征结构的相干算法，抗噪能力强且对横向振幅变化更加敏感，但是受计算参数影响大。此后国内外学者多在以上三种算法的基础上进行改进和完善，进一步提高对断裂等构造解释的准确性和合理性。然而，这三种算法对于中小尺度断裂存在成像效果不佳甚至无法成像的问题，常常无法检测和识别中小尺度断裂。

技术实现思路

1、基于此，有必要针对上述技术问题，提供一种断裂精细切片计算方法、装置、计算机设备和存储介质。

2、一种断裂精细切片计算方法，包括：

3、获取原始地震数据；

4、基于所述原始地震数据和互相关公式，获取第一倾角时差和第二倾角时差，其中，所述第一倾角时差为沿主测线方向的相邻道的倾角时差，所述第二倾角时差为联络测线方向的相邻道的倾角时差；

5、基于所述第一倾角时差和第二倾角时差，获取相干数据体；

6、基于所述相干数据体和层位信息，沿层提取均方根振幅信息；

7、利用八方向各向异性滤波器对所述均方根振幅信息进行滤波，获取断裂数据体；

8、对所述断裂数据体进行计算，获取断裂精细切片。

9、在其中一个实施例中，所述互相关公式为：

10、

11、

12、其中，si(t)为沿主测线方向的地震信号，sj(t)为联络测线方向的地震信号，si-1(t)为沿主测线方向的相邻地震信号，sj-1(t)为联络测线方向的相邻地震信号，τ为互相关扫描时差，k为互相关时窗的大小；

13、所述第一倾角时差为：

14、δti＝max(ci)

15、所述第二倾角时差为：

16、δtj＝max(cj)。

17、在其中一个实施例中，所述相干数据体为：

18、

19、其中，为t时刻的相干值，s1(t)为计算的地震道，s2(t)为相邻的地震道，为相邻地震道的倾角时差，l为计算相干的时窗大小。

20、在其中一个实施例中，所述均方根振幅信息为：

21、

22、其中，m为沿层位提取均方根振幅信息所设定的纵向时窗大小。

23、在其中一个实施例中，所述中小尺度断裂数据体为：

24、iitr+1＝iitr+λ(cn▽n(iitr)+cs▽s(iitr)+ce▽e(iitr)+cw▽w(iitr)+cne▽ne(iitr)+cnw▽nw(iitr)+cse▽se(iitr)+csw▽sw(iitr))

25、其中，itr为迭代次数，λ为平滑系数，

26、▽(iitr)为八个方向上的散度，具体展开式为：

27、▽n(iitr)＝iitr(x,y)-iitr(x,y-1)

28、▽s(iitr)＝iitr(x,y)-iitr(x,y+1)

29、▽e(iitr)＝iitr(x,y)-iitr(x-1,y)

30、▽w(iitr)＝iitr(x,y)-iitr(x+1,y)

31、▽ne(iitr)＝iitr(x,y)-iitr(x-1,y-1)

32、▽nw(iitr)＝iitr(x,y)-iitr(x+1,y-1)

33、▽se(iitr)＝iitr(x,y)-iitr(x-1,y+1)

34、▽sw(iitr)＝iitr(x,y)-iitr(x+1,y+1)

35、其中，c代表八个方向上的传递系数，为与传递系数相关的扩散系数β，

36、cn＝exp(-||▽n(iitr)||2/β2)

37、cs＝exp(-||▽s(iitr)||2/β2)

38、ce＝exp(-||▽e(iitr)||2/β2)

39、cw＝exp(-||▽w(iitr)||2/β2)。

40、cne＝exp(-||▽ne(iitr)||2/β2)

41、cse＝exp(-||▽se(iitr)||2/β2)

42、csw＝exp(-||▽sw(iitr)||2/β2)

43、在其中一个实施例中，所述断裂精细切片为：

44、

45、其中，tw为半时窗的大小，t0为提取切片所对应的时间深度。

46、在其中一个实施例中，所述均方根振幅信息通过二维图像的方式呈现。

47、一种断裂精细切片计算装置，包括：

48、第一获取模块，用于获取原始地震数据；

49、第二获取模块，用于基于所述原始地震数据和互相关公式，获取第一倾角时差和第二倾角时差，其中，所述第一倾角时差为沿主测线方向的相邻道的倾角时差，所述第二倾角时差为联络测线方向的相邻道的倾角时差；

50、第三获取模块，用于基于所述第一倾角时差和第二倾角时差，获取相干数据体；

51、提取模块，用于基于所述相干数据体和层位信息，沿层提取均方根振幅信息；

52、滤波模块，用于利用八方向各向异性滤波器对所述均方根振幅信息进行滤波，获取断裂数据体；

53、第四获取模块，用于基于所述中小尺度断裂数据体和提取切片公式，获取断裂精细切片。

54、一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：

55、获取原始地震数据；

56、基于所述原始地震数据和互相关公式，获取第一倾角时差和第二倾角时差，其中，所述第一倾角时差为沿主测线方向的相邻道的倾角时差，所述第二倾角时差为联络测线方向的相邻道的倾角时差；

57、基于所述第一倾角时差和第二倾角时差，获取相干数据体；

58、基于所述相干数据体和层位信息，沿层提取均方根振幅信息；

59、利用八方向各向异性滤波器对所述均方根振幅信息进行滤波，获取断裂数据体；

60、对所述断裂数据体进行计算，获取断裂精细切片。

61、一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

62、获取原始地震数据；

63、基于所述原始地震数据和互相关公式，获取第一倾角时差和第二倾角时差，其中，所述第一倾角时差为沿主测线方向的相邻道的倾角时差，所述第二倾角时差为联络测线方向的相邻道的倾角时差；

64、基于所述第一倾角时差和第二倾角时差，获取相干数据体；

65、基于所述相干数据体和层位信息，沿层提取均方根振幅信息；

66、利用八方向各向异性滤波器对所述均方根振幅信息进行滤波，获取断裂数据体；

67、对所述断裂数据体进行计算，获取断裂精细切片。

68、上述断裂精细切片计算方法、装置、计算机设备和存储介质，引入了倾角时差的相干属性，该属性能够降低常规相干引起的横向的背景噪音，并且，还使用了八方向各向异性滤波器，增强了不同方向上的差异性，因此能够对一些中尺度和小尺度的断裂进行增强。通过对整个断裂数据体进行计算后，就能够得到断裂精细切片，所得到的断裂精细切片具有较好的稳定性和抗噪性，相较于传统的断裂切片，有利于中小尺度断裂及裂缝的预测和描述更加精细。