一种双方位拖缆地震资料螺旋角道集获取方法及系统与流程

本发明属于地震数据处理，具体涉及一种双方位拖缆地震资料螺旋角道集获取方法及系统。

背景技术：

1、潜山储层已经成为国内海域油气勘探开发的一个重要领域，裂缝预测对于落实潜山储层有着重要的作用。目前，业界开展叠前裂缝预测主要是针对宽方位地震资料，基于各向异性理论，通过不同方位地震波的传播和衰减理论，建立裂缝型介质的岩石物理模型，利用方位各向异性的椭圆拟合和各向异性反演等方法进行断裂及裂缝储层描述。海上地震勘探一般采用拖缆采集，由于采集方式的局限性，拖缆采集的横纵比通常小于0.05，不具备开展叠前裂缝预测的基础，近几年针对海上探区潜山储层的裂缝预测研究多以叠后为主。

2、随着海上二次三维采集的部署，很多研究区存在两套或以上的不同采集方位的拖缆地震资料，这些地震资料为开展叠前裂缝预测创造了条件。不同时期采集的拖缆地震资料都有其针对性的研究目的，相应的拖缆资料采集参数(包括震源参数、枪缆沉放深度、观测系统等)也存在明显的差异，造成不同采集方位的地震资料存在振幅、频率、相位、时差、面元属性等区块间不一致的问题。基于不同采集方位拖缆资料的处理，首先需要有效消除区块间的各种差异，再采用合理的处理流程，最终获得双方位融合的地震成果。现有的双方位融合处理流程能够挖掘不同方位地震资料的潜力、改善成像效果，但对于研究方位各向异性特征、开展叠前裂缝预测来说并没有起到根本性的改善。

技术实现思路

1、本发明的目的是提供一种双方位拖缆地震资料螺旋角道集获取方法及系统，用以解决现有技术中存在的上述问题。

2、为了实现上述目的，本发明采用以下技术方案：

3、第一方面，提供一种双方位拖缆地震资料螺旋角道集获取方法，包括：

4、a.获取方位一的原始拖缆地震资料和方位二的原始拖缆地震资料，所述方位一对应的有效方位扇区为α1-α2，其中α1和α2分别为方位一的方位角下限和方位角上限，所述方位二对应的有效方位扇区为β1-β2，其中，β1和β2分别为方位二的方位角下限和方位角上限；

5、b.分别对方位一的原始拖缆地震资料和方位二的原始拖缆地震资料进行预处理，得到对应的方位一地震炮集数据a和方位二地震炮集数据b，并分别对方位一地震炮集数据a和方位二地震炮集数据b进行地震数据规则化处理，得到对应的方位一数据体areg和方位二数据体breg；

6、c.对方位一数据体areg和方位二数据体breg开展双方位各向异性叠前深度偏移速度建模，得到双方位融合的各向异性速度v，以及方位一对应的各向异性速度va和方位二对应的各向异性速度vb；

7、d.对方位一地震炮集数据a和方位二地震炮集数据b进行五维数据规则化处理，得到规则化数据集c，所述规则化数据集c对应的有效方位扇区包括α1-α2和β1-β2；

8、e.分别对有效方位扇区α1-α2和β1-β2进行方位扇区划分，得到有效方位扇区α1-α2对应的三个小方位扇区a1、a2和a3，以及有效方位扇区β1-β2对应的三个小方位扇区b1、b2和b3，其中，小方位扇区a1和小方位扇区a3的方位角范围等大且大于小方位扇区a2的方位角范围，小方位扇区b1和小方位扇区b3的方位角范围等大且大于小方位扇区b2的方位角范围；

9、f.根据小方位扇区a1、a2和a3以及小方位扇区b1、b2和b3对规则化数据集c进行划分，得到小方位扇区a1对应的数据体a1，小方位扇区a2对应的数据体a2，小方位扇区a3对应的数据体a3，小方位扇区b1对应的数据体b1，小方位扇区b2对应的数据体b2以及小方位扇区b3对应的数据体b3；

10、g.基于各向异性速度v和各向异性速度va，分别利用数据体a1、a2和a3开展各向异性叠前深度偏移处理，得到对应的叠前深度偏移crp道集pa1、pa2和pa3，基于各向异性速度v和各向异性速度vb，分别利用数据体b1、b2和b3，开展各向异性叠前深度偏移处理，得到对应的叠前深度偏移crp道集pb1、pb2和pb3；

11、h.利用设定的数据一致性判定条件对叠前深度偏移crp道集pa1、pa2和pa3进行数据一致性判定，利用设定的数据一致性判定条件对叠前深度偏移crp道集pb1、pb2和pb3进行数据一致性判定，当叠前深度偏移crp道集pa1、pa2和pa3之间不满足数据一致性判定条件时，更改有效方位扇区α1-α2的方位扇区划分，得到更新后的三个小方位扇区a1、a2和a3，当叠前深度偏移crp道集pb1、pb2和pb3之间不满足数据一致性判定条件时，更改有效方位扇区β1-β2的方位扇区划分，得到更新后的三个小方位扇区b1、b2和b3；

12、i.基于更新后的三个小方位扇区a1、a2和a3，和/或更新后的三个小方位扇区b1、b2和b3，迭代执行步骤f-h，直至叠前深度偏移crp道集pa1、pa2和pa3之间满足数据一致性判定条件，且叠前深度偏移crp道集pb1、pb2和pb3之间满足数据一致性判定条件，得到最终的叠前深度偏移crp道集pa1、pa2和pa3，以及最终的叠前深度偏移crp道集pb1、pb2和pb3；

13、j.从最终的叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3中分别提取井旁地震道集，并获取对应的井正演道集，利用各井正演道集对相应的井旁地震道集进行avo特征评价，确定对应的切除参数，利用对应的切除参数对最终的叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3进行切除处理，并将切除处理后的叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3分别转入射角道集，得到对应的角道集数据体ga1、ga2、ga3、gb1、gb2和gb3；

14、k.将角道集数据体ga1、ga2、ga3、gb1、gb2和gb3合并，按照螺旋角道集方式置道头并进行排序，得到螺旋角道集数据体g。

15、在一个可能的设计中，在得到螺旋角道集数据体g后，所述方法还包括：

16、l.对螺旋角道集数据体g进行保幅处理，得到对应的高信噪比螺旋角道集数据体gden，并输出高信噪比螺旋角道集数据体gden。

17、在一个可能的设计中，所述分别对方位一的原始拖缆地震资料和方位二的原始拖缆地震资料进行预处理，包括：

18、分别对方位一的原始拖缆地震资料和方位二的原始拖缆地震资料进行保真去噪处理、宽频处理和多次波压制处理，以及对方位一的原始拖缆地震资料与方位二的原始拖缆地震资料进行振幅、频率、相位和时差的一致性处理。

19、在一个可能的设计中，所述分别对方位一地震炮集数据a和方位二地震炮集数据b进行地震数据规则化处理，包括：

20、采用基于离散傅里叶变换的数据规则化方法分别对方位一地震炮集数据a和方位二地震炮集数据b进行地震数据规则化处理。

21、在一个可能的设计中，所述对方位一地震炮集数据a和方位二地震炮集数据b进行五维数据规则化处理，包括：

22、采用基于离散傅里叶变换的五维数据规则化处理方法对方位一地震炮集数据a和方位二地震炮集数据b进行五维数据规则化处理。

23、在一个可能的设计中，所述数据一致性判定条件包括：平均覆盖次数差异不大于20％，目的层平均信噪比差异不大于20％，统计时窗内的平均均方根振幅差异不大于10％。

24、在一个可能的设计中，所述利用各井正演道集对相应的井旁地震道集进行avo特征评价，确定对应的切除参数，包括：

25、统计井正演道集的avo变化趋势，并根据井正演道集的avo变化趋势选定有效偏移距范围，利用有效偏移距范围对相应的井旁地震道集进行avo特征评价，确定其各目的层的可用偏移距范围，基于各目的层的可用偏移距范围确定对应的切除参数。

26、第二方面，提供一种双方位拖缆地震资料螺旋角道集获取系统，包括数据获取单元、数据整理单元、速度建模单元、规则处理单元、扇区划分单元、数据划分单元、偏移处理单元、差异判定单元、切除转换单元和数据合并单元，其中：

27、数据获取单元，用于获取方位一的原始拖缆地震资料和方位二的原始拖缆地震资料，所述方位一对应的有效方位扇区为α1-α2，其中α1和α2分别为方位一的方位角下限和方位角上限，所述方位二对应的有效方位扇区为β1-β2，其中，β1和β2分别为方位二的方位角下限和方位角上限；

28、数据整理单元，用于分别对方位一的原始拖缆地震资料和方位二的原始拖缆地震资料进行预处理，得到对应的方位一地震炮集数据a和方位二地震炮集数据b，并分别对方位一地震炮集数据a和方位二地震炮集数据b进行地震数据规则化处理，得到对应的方位一数据体areg和方位二数据体breg；

29、速度建模单元，用于对方位一数据体areg和方位二数据体breg开展双方位各向异性叠前深度偏移速度建模，得到双方位融合的各向异性速度v，以及方位一对应的各向异性速度va和方位二对应的各向异性速度vb；

30、规则处理单元，用于对方位一地震炮集数据a和方位二地震炮集数据b进行五维数据规则化处理，得到规则化数据集c，所述规则化数据集c对应的有效方位扇区包括α1-α2和β1-β2；

31、扇区划分单元，用于分别对有效方位扇区α1-α2和β1-β2进行方位扇区划分，得到有效方位扇区α1-α2对应的三个小方位扇区a1、a2和a3，以及有效方位扇区β1-β2对应的三个小方位扇区b1、b2和b3，其中，小方位扇区a1和小方位扇区a3的方位角范围等大且大于小方位扇区a2的方位角范围，小方位扇区b1和小方位扇区b3的方位角范围等大且大于小方位扇区b2的方位角范围；

32、数据划分单元，用于根据小方位扇区a1、a2和a3以及小方位扇区b1、b2和b3对规则化数据集c进行划分，得到小方位扇区a1对应的数据体a1，小方位扇区a2对应的数据体a2，小方位扇区a3对应的数据体a3，小方位扇区b1对应的数据体b1，小方位扇区b2对应的数据体b2以及小方位扇区b3对应的数据体b3；

33、偏移处理单元，用于基于各向异性速度v和各向异性速度va，分别利用数据体a1、a2和a3开展各向异性叠前深度偏移处理，得到对应的叠前深度偏移crp道集pa1、pa2和pa3，基于各向异性速度v和各向异性速度vb，分别利用数据体b1、b2和b3，开展各向异性叠前深度偏移处理，得到对应的叠前深度偏移crp道集pb1、pb2和pb3；

34、差异判定单元，用于利用设定的数据一致性判定条件对叠前深度偏移crp道集pa1、pa2和pa3进行数据一致性判定，利用设定的数据一致性判定条件对叠前深度偏移crp道集pb1、pb2和pb3进行数据一致性判定，当叠前深度偏移crp道集pa1、pa2和pa3之间不满足数据一致性判定条件时，更改有效方位扇区α1-α2的方位扇区划分，得到更新后的三个小方位扇区a1、a2和a3，当叠前深度偏移crp道集pb1、pb2和pb3之间不满足数据一致性判定条件时，更改有效方位扇区β1-β2的方位扇区划分，得到更新后的三个小方位扇区b1、b2和b3；

35、切除转换单元，用于从叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3中分别提取井旁地震道集，并获取对应的井正演道集，利用各井正演道集对相应的井旁地震道集进行avo特征评价，确定对应的切除参数，利用对应的切除参数对最终的叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3进行切除处理，并将切除处理后的叠前深度偏移crp道集pa1、pa2、pa3、pb1、pb2和pb3分别转入射角道集，得到对应的角道集数据体ga1、ga2、ga3、gb1、gb2和gb3；

36、数据合并单元，用于将角道集数据体ga1、ga2、ga3、gb1、gb2和gb3合并，按照螺旋角道集方式置道头并进行排序，得到螺旋角道集数据体g。

37、在一个可能的设计中，所述系统还包括保幅处理单元，所述保幅处理单元用于对螺旋角道集数据体g进行保幅处理，得到对应的高信噪比螺旋角道集数据体gden，并输出高信噪比螺旋角道集数据体gden。

38、第三方面，提供一种双方位拖缆地震资料螺旋角道集获取系统，包括：

39、存储器，用于存储指令；

40、处理器，用于读取所述存储器中存储的指令，并根据指令执行上述第一方面中任意一种所述的方法。

41、有益效果：

42、1、本发明针对双方位炮集数据开展五维规则化处理，获得有效方位扇区范围内不同面元之间偏移距以及方位角分布更理想的规则化效果。

43、2、本发明通过优选出合理的不等分扇区划分方式，确保了不同方位扇区成果的一致性。

44、3、本发明基于井震一体化分析定量确定最大可用的偏移距范围，确保了方位各向异性属性的可靠性。

45、4、本发明针对螺旋入射角道集开展保幅的提高信噪比处理，确保角道集有较高的信噪比，满足较窄方位数据开展叠前裂缝预测的需求。