一种侵入岩体的流动方向确定方法、存储介质以及设备与流程

1.本发明涉及构造解析领域，尤其涉及一种侵入岩体的流动方向确定方法、存储介质以及设备。


背景技术：

2.岩浆侵入是内动力地质作用的重要研究内容，是地壳构造变形的重要力学机制之一。国内外研究表明，岩体侵入遵循一定的应力规律，具有定向特征，且与侵入体的几何学样式密切相关。
3.以塔里木盆地为例，塔里木盆地的台盆区(如顺北地区)自古生代以来岩浆活动频繁，地震和钻井资料均揭示了该区深层(古生界)广泛分布走滑断裂体系和大量火成岩侵入体。火成岩侵入体的存在造成走滑断裂成像和储层预测工作十分困难，进而影响区带评价和井位部署。因此，能够建立走滑背景下火成岩侵入体的发育模式并判识其侧向流动方向对于海相盆地深层的选区评价和井轨迹优化都具有重要的意义。
4.目前，关于侵入岩体流动方向的判定技术可总结为以下两种：
5.1)利用岩石样品中的稀有气体同位素(he、ne和ar同位素)变化，判定岩体侵入方向。
6.2)根据不同尺度的剩余重力异常中心的位置与出露岩体中心的位置，判断出岩体的侵入方向。
7.然而，现有的侵入岩体流动方向判定技术适用范围局限，仅针对岩石样品丰富或有岩体出露的地区，而无法适用于岩石样品稀缺且处于盆地深层的侵入岩体流动方向判定。另外，对于广泛发育有高陡走滑断裂体系和大量火成岩的盆区深层区，侵入岩体的存在会造成高陡走滑断裂体系成像和储层预测工作十分困难，进而对侵入岩体流动方向的识别将影响到后续的区带评价以及井位部署等。


技术实现要素：

8.本发明要解决的技术问题是：如何有效确定深层侵入岩体的流动方向，从而能够预测研究的三维区周边非地震资料覆盖区火成岩体展布特征，还可以为火成岩体勘探领域提供依据，实现井位部署和井轨迹优化。
9.为解决上述技术问题，本发明提供了一种侵入岩体的流动方向确定方法、存储介质以及设备。
10.本发明的第一个方面，提供了一种侵入岩体的流动方向确定方法，其包括：
11.识别三维区内的侵入岩体以及识别出各地震剖面内所述侵入岩体的岩床模型；
12.获取不同地震剖面内所述岩床模型的地震反射特征，基于所述地震反射特征确定出目标岩床模型并对所述目标岩床模型进行地震层位追踪，以获取所述目标岩床模型的三维空间形态；
13.基于三维地震资料，确定出所述侵入岩体的邻近层中断裂体系的展布特征，所述
断裂体系包括位于所述侵入岩体上层的逆断层和高陡裂隙；
14.基于所述目标岩床模型的三维空间形态和所述断裂体系的展布特征，确定所述侵入岩体的流动方向。
15.在一些实施例中，所述地震剖面包括时间地震剖面，所述识别三维区内的侵入岩体包括：
16.获取所述三维区内的所述三维地震资料并从所述三维地震资料中提取所述时间地震剖面；
17.基于所述时间地震剖面，根据所述侵入岩体的典型地震反射特征，识别所述三维区内的所述侵入岩体。
18.在一些实施例中，所述识别出各地震剖面内所述侵入岩体的岩床模型，包括：
19.基于发育要素，识别出各时间地震剖面内所述侵入岩体的岩床模型。
20.在一些实施例中，所述岩床模型包括碟形岩床，所述基于发育要素，识别出各时间地震剖面内所述侵入岩体的岩床模型，包括：
21.基于内部岩床、倾斜岩脉和外部岩床，识别出各时间地震剖面内所述侵入岩体的所述岩床模型。
22.在一些实施例中，所述获取不同地震剖面内所述岩床模型的地震反射特征，基于所述地震反射特征确定出目标岩床模型，包括：
23.获取不同时间地震剖面内所述岩床模型的所述地震反射特征，确定出所述地震反射特征最清晰的所述岩床模型，作为所述目标岩床模型。
24.在一些实施例中，对所述目标岩床模型进行地震层位追踪之后，所述方法还包括：对所述目标岩床模型的解释层位进行插值加密处理。
25.在一些实施例中，所述基于三维地震资料，确定出所述侵入岩体的邻近层中断裂体系的展布特征，包括：
26.基于所述三维地震资料提取所述侵入岩体邻近层的相干属性图；
27.基于所述相干属性图确定出所述侵入岩体的邻近层中断裂体系的展布特征。
28.在一些实施例中，基于所述目标岩床模型的三维空间形态和所述断裂体系的展布特征，确定所述侵入岩体的流动方向，包括：
29.基于所述目标岩床模型的三维空间形态确定所述侵入岩体的第一流动方向；
30.根据所述逆断层和所述侵入岩体的流动方向垂直的特征，确定出所述侵入岩体的第二流动方向；
31.基于所述第一流动方向和所述第二流动方向确定出所述侵入岩体的流动方向。
32.本发明的第二个方面，提供了一种存储介质，所述存储介质中存储有计算机程序，所述计算机程序被处理器执行时能够实现如上任意一项所述的侵入岩体的流动方向确定方法。
33.本发明的第三个方面，提供了一种设备，所述设备包括处理器和存储器，所述存储器中存储有计算机程序，所述处理器执行所述存储器中存储的计算机程序时能够实现如上任意一项所述的侵入岩体的流动方向确定方法。
34.与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
35.应用本发明提供的侵入岩体的流动方向确定方法，通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的地震反射特征确定出目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态；另外，基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，分别基于三维空间形态以及断裂体系展布特征和侵入岩体流动方向的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
附图说明
36.通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
37.图1示出了本发明实施例一提供的一种侵入岩体的流动方向确定方法流程示意图；
38.图2示出了碟形岩床“指状”的空间展布示意图；
39.图3示出了火成岩相关构造示意图；
40.图4示出了本发明实施例二提供的一种侵入岩体的流动方向确定方法流程示意图；
41.图5示出了本发明实施例三提供的一种侵入岩体的流动方向确定方法流程示意图；
42.图6a示出了本发明实施例提供的顺北地区三维区典型时间地震剖面示意图；图6b示出了本发明实施例提供的顺北地区三维区精细识别获得的侵入岩体时间地震剖面示意图；
43.图7a示出了碟形岩床完整发育模式的示意图；图7b示出了本发明实施例提供的顺北地区三维区碟形岩床时间地震剖面示意图；
44.图8示出了顺北地区三维区内北西-东南向侵入岩体发育特征示意图；
45.图9a示出了本发明实施例提供的顺北地区三维区碟形岩床立体显示图；图9b示出了沿图9a中a-a’得到的“指状”展布时间地震剖面示意图；
46.图10a示出了本发明实施例提供的顺北地区三维区t
60
界面相干属性图；图10b示出了本发明实施例提供的顺北地区三维区t
60
界面火成岩相关构造精细解释示意图；
47.图11a示出了沿图10b中a-a’得到的顺北地区三维区火成岩相关逆断层典型时间地震剖面示意图；图11b示出了沿图10b中b-b’得到的顺北地区三维区火成岩相关高陡裂隙典型时间地震剖面示意图；
48.图12示出了本发明实施例提供的顺北地区三维区侵入岩体流动方向示意图；
49.图13a示出了三维区玛纳侵入岩体分布位置示意图；图13b示出了a井采集到的花岗斑岩图；
50.图14示出了本发明实施例提供的一种设备示意图。
具体实施方式
51.为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
52.岩浆侵入是内动力地质作用的重要研究内容，是地壳构造变形的重要力学机制之一。国内外研究表明，岩体侵入遵循一定的应力规律，具有定向特征，且与侵入体的几何学样式密切相关。
53.以塔里木盆地为例，塔里木盆地的台盆区(如顺北地区)自古生代以来岩浆活动频繁，地震和钻井资料均揭示了该区深层(古生界)广泛分布走滑断裂体系和大量火成岩侵入体。火成岩侵入体的存在造成走滑断裂成像和储层预测工作十分困难，进而影响区带评价和井位部署。因此，能够建立走滑背景下火成岩侵入体的发育模式并判识其侧向流动方向对于海相盆地深层的选区评价和井轨迹优化都具有重要的意义。
54.目前，关于侵入岩体流动方向的判定技术可总结为以下两种：
55.1)利用岩石样品中的稀有气体同位素(he、ne和ar同位素)变化，判定岩体侵入方向。
56.2)根据不同尺度的剩余重力异常中心的位置与出露岩体中心的位置，判断出岩体的侵入方向。
57.然而，现有的侵入岩体流动方向判定技术适用范围局限，仅针对岩石样品丰富或有岩体出露的地区，而无法适用于岩石样品稀缺且处于盆地深层的侵入岩体流动方向判定。另外，对于广泛发育有高陡走滑断裂体系和大量火成岩的盆区深层区，侵入岩体的存在会造成高陡走滑断裂体系成像和储层预测工作十分困难，进而对侵入岩体流动方向的识别将影响到后续的区带评价以及井位部署等。
58.有鉴于此，本发明提供了一种侵入岩体的流动方向确定方法，通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的地震反射特征确定出目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态；另外，基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，分别基于三维空间形态以及断裂体系展布特征和侵入岩体流动方向的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
59.需要说明的是，本发明提供的侵入岩体的流动方向确定方法可以适用于盆地深层走滑背景下侵入岩体的流动方向确定，具体请参见下文中的描述。
60.实施例一
61.参见图1所示，图1示出了本发明实施例一提供的一种侵入岩体的流动方向确定方法流程示意图，其包括：
62.步骤s101：识别三维区内的侵入岩体以及识别出各地震剖面内侵入岩体的岩床模型。
63.步骤s102：获取不同地震剖面内岩床模型的地震反射特征，基于地震反射特征确定出目标岩床模型并对目标岩床模型进行地震层位追踪，以获取目标岩床模型的三维空间
形态。
64.步骤s103：基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，断裂体系包括位于侵入岩体上层的逆断层和高陡裂隙。
65.步骤s104：基于目标岩床模型的三维空间形态和断裂体系的展布特征，确定侵入岩体的流动方向。
66.在本发明实施例中，步骤s101可以具体为，获取三维区内的三维地震资料并从三维地震资料中提取地震剖面；基于地震剖面，根据侵入岩体的典型地震反射特征，识别三维区内的侵入岩体。
67.在一些实施例中，地震剖面可以为时间地震剖面；在另一些实施例中，地震剖面还可以为深度地震剖面，在下文描述中，将以时间地震剖面为例进行描述。
68.基于时间地震剖面，根据侵入岩体的典型地震反射特征识别三维区内的侵入岩体可以为，根据地震反射特征出现振幅反射异常如振幅增强且终止突然的特征，从时间地震剖面中识别出侵入岩体。
69.另外，在步骤s101中识别各地震剖面内侵入岩体的岩床模型可以为：基于发育要素，识别出各时间地震剖面内侵入岩体的岩床模型。在一些实施例中，基于发育要素识别出各时间地震剖面内侵入岩体的岩床模型可以为，基于内部岩床、倾斜岩脉和外部岩床，识别出各时间地震剖面内侵入岩体的岩床模型。其中，内部岩床、倾斜岩脉和外部岩床位可以为碟形岩床的主要发育要素，基于碟形岩床的主要发育要素，可以从已识别出的侵入岩体中确定出碟形岩床，即岩床模型可以为碟形岩床，具体的，可以根据时间地震剖面中侵入岩体的剖面几何形态和碟形岩床的主要发育要素进行对比，从而识别出侵入岩体中的碟形岩床。在另一些实施例中，还可以基于碟形岩床完整发育模式中的发育要素识别出各时间地震剖面内侵入岩体的岩床模型，其中，碟形岩床完整发育模式中的发育要素可以包括内部岩床、倾斜岩脉、外部岩床和供给岩脉。
70.在本发明实施例中，步骤s102可以具体为，基于三维地震资料，获取不同时间地震剖面内岩床模型的地震反射特征，根据精确度的需要，可以选择满足预设要求的地震反射特征确定出目标岩床模型。
71.在一些实施例中，可以通过获取不同时间地震剖面内岩床模型的地震反射特征，根据精确度需求，选择地震反射特征最清晰的岩床模型，作为目标岩床模型。
72.另外，针对确定出的目标岩床模型还可以进行地震层位追踪，对目标岩床模型进行立体显示，以获取目标岩床模型的三维空间形态。
73.基于目标岩床模型的三维空间形态可以确定出侵入岩体的大致流动方向，参见图2所示，图2示出了碟形岩床“指状”的空间展布示意图，侵入岩体自深部沿着垂向断裂通道向上涌入，在合适的地层界面侵入并在偏离垂向方向的侧向方向上发生流动。在侧向方向上流动过程中，由于流体密度差异、地层岩性和物性差异，侵入岩体并不能沿同一界面均匀铺开，而是逐渐分散成指状流动。因此，可以基于空间展布情况大致识别出，侵入流体的流动方向大致是指向指状分叉的方向。
74.在另一些实施例中，还可以基于三维地震资料，获取不同时间地震剖面内岩床模型的发育规律参量，其中，发育规律参量可以包括侵入层位、发育数量和上覆地层变形强度中的至少一种。通过获取岩床模型的发育规律参量可以获取三维区内侵入岩体的发育规
律，侵入岩体的发育规律和侵入岩体的流动方向相关，例如，侵入岩体可沿多个界面侵入，形成垂向叠置特征，沿侵入岩体的流动方向，岩床模型的数量增多、最高侵入层位增高并且上覆地层变形强度增加。从而，基于发育规律可以确定出侵入岩体的大致流动方向，进一步可实现对基于目标岩床模型的三维空间形态确定出的大致流动方向进行验证。
75.步骤s103可以具体为：基于三维地震资料提取侵入岩体邻近层的相干属性图；基于相干属性图确定出侵入岩体的邻近层中断裂体系的展布特征。
76.其中，侵入岩体侵入时会在侵入岩体上部形成断裂体系，侵入岩体侵入时引起岩床上覆地层变形形成高陡裂隙，同时在岩床边缘形成逆断层，高陡裂隙和逆断层构成火成岩相关构造，也可称为次级断裂体系；侵入岩体的下部为走滑断裂体系。参见图3所述，图3示出了火成岩相关构造示意图。通常，高陡裂隙和逆断层的走向接近垂直，逆断层和侵入岩体的流动方向垂直。在一些实施例中，基于高陡裂隙和逆断层接近垂直的组合方式，还可以验证断裂体系是由侵入岩体的侵入引起的。在本发明实施例中，断裂体系可以包括位于目标岩床模型上部的火成岩相关构造即次级断裂体系，基于三维地震资料提取侵入岩体邻近层的相干属性图，可以为基于三维地震资料提取侵入岩体上部邻近层的相干属性图，基于相干属性图可以确定出侵入岩体的火成岩相关构造，其中，火成岩相关构造可以包括逆断层和高陡裂隙，逆断层发育在岩床模型的倾斜岩脉之上，逆断层的展布特征受控于倾斜岩脉的展布，高陡裂隙发育在岩床之上，通常对应岩脉错断的地方。基于逆断层的展布方向和侵入岩体的流动方向垂直的规律，可以精确确定出侵入岩体的展布方位。
77.进而，在本发明实施例中，步骤s104可以为：基于目标岩床模型的三维空间形态确定侵入岩体的第一流动方向；根据逆断层和侵入岩体的流动方向垂直的特征，确定侵入岩体的第二流动方向；通过结合第一流动方向和第二流动方向确定出侵入岩体的流动方向。其中，第一流动方向可以为侵入岩体的大致流动方向，第二流动方向可以为展布方位。作为一示例，基于目标三维空间形态如指状展布形态可以大致识别出侵入岩体的第一流动方向为自北向南流动，基于目标岩床模型的邻近层中断裂体系的展布特征可以识别出侵入岩体的第二流动方向为与逆断层垂直，综合第一流动方向和第二流动方向即可确定出侵入岩体的流动方向为向南流动过程中与逆断层垂直，从而可以精确确定出侵入岩体的流动方向。
78.以上为本发明实施例一提供的一种侵入岩体的流动方向确定方法，通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的地震反射特征确定出目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态；另外，基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，基于三维空间形态以及断裂体系展布特征和侵入岩体流动的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
79.在本发明实施例中，为了提高确定侵入岩体流动方向的精确度，还可以选择最具有代表性的岩床模型作为目标岩床模型，具体请参见以下实施例二中的描述。
80.实施例二
81.参见图4所示，图4示出了本发明实施例二提供的一种侵入岩体的流动方向确定方法流程示意图，其包括：
82.步骤s201：识别三维区内的侵入岩体以及识别出各时间地震剖面内侵入岩体的岩床模型。
83.步骤s202：获取不同时间地震剖面内岩床模型的地震反射特征，确定出地震反射特征最清晰的岩床模型，作为目标岩床模型并对目标岩床模型进行地震层位追踪，以获取目标岩床模型的三维空间形态。
84.步骤s203：基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，断裂体系包括位于侵入岩体上层的逆断层和高陡裂隙。
85.步骤s204：基于目标岩床模型的三维空间形态和断裂体系的展布特征，确定侵入岩体的流动方向。
86.在本发明实施例中，步骤s201可以具体为，获取三维区内的三维地震资料并从三维地震资料中提取时间地震剖面；基于时间地震剖面，根据侵入岩体的典型地震反射特征，识别三维区内的侵入岩体。其中，为了得到更准确、清晰的时间地震剖面，可以获取三维区内的高精度三维地震资料。
87.基于时间地震剖面，根据侵入岩体的典型地震反射特征识别三维区内的侵入岩体可以为，根据地震反射特征出现振幅反射异常如振幅增强且终止突然的特征，从时间地震剖面中识别出侵入岩体。在另一些实施例中，还可以基于高精度三维地震资料提取深度地震剖面。
88.另外，在步骤s201中识别各时间地震剖面内侵入岩体的岩床模型可以为：基于发育要素，识别出各时间地震剖面内侵入岩体的岩床模型。在一些实施例中，基于发育要素识别出各时间地震剖面内侵入岩体的岩床模型可以为，基于内部岩床、倾斜岩脉和外部岩床，识别出各时间地震剖面内侵入岩体的岩床模型。其中，内部岩床、倾斜岩脉和外部岩床位可以为碟形岩床的主要发育要素，基于碟形岩床的主要发育要素，可以从已识别出的侵入岩体中确定出碟形岩床，即岩床模型可以为碟形岩床，具体的，可以根据时间地震剖面中侵入岩体的剖面几何形态和碟形岩床的主要发育要素进行对比，从而识别出侵入岩体中的碟形岩床。在另一些实施例中，还可以基于碟形岩床完整发育模式中的发育要素识别出各时间地震剖面内侵入岩体的岩床模型，其中，碟形岩床完整发育模式中的发育要素可以包括内部岩床、倾斜岩脉、外部岩床和供给岩脉。
89.在本发明实施例中，步骤s202可以具体为，基于三维地震资料，获取不同时间地震剖面内岩床模型的地震反射特征，通过确定出地震反射特征最清晰的岩床模型，作为目标岩床模型。
90.另外，针对确定出的目标岩床模型还可以进行地震层位追踪，对目标岩床模型进行立体显示，以获取目标岩床模型的三维空间形态。
91.基于目标岩床模型的三维空间形态可以确定出侵入岩体的大致流动方向，参见图2所示，图2示出了碟形岩床“指状”的空间展布示意图，侵入岩体自深部沿着垂向断裂通道向上涌入，在合适的地层界面侵入并偏离垂向方向的侧向方向上发生流动。在侧向方向上流动过程中，由于流体密度差异、地层岩性和物性差异，侵入岩体并不能沿同一界面均匀铺开，而是逐渐分散成指状流动。因此，可以基于空间展布情况大致识别出，侵入流体的流动方向大致是指向指状分叉的方向。
92.在另一些实施例中，还可以基于三维地震资料，获取不同时间地震剖面内岩床模
型的发育规律参量，其中，发育规律参量可以包括侵入层位、发育数量和上覆地层变形强度中的至少一种。通过获取岩床模型的发育规律参量可以获取三维区内侵入岩体的发育规律，侵入岩体的发育规律和侵入岩体的流动方向相关，例如，侵入岩体可沿多个界面侵入，形成垂向叠置特征，沿侵入岩体的流动方向，岩床模型的数量增多、最高侵入层位增高并且上覆地层变形强度增加。从而，基于发育规律可以确定出侵入岩体的大致流动方向，进一步可实现对基于目标岩床模型的三维空间形态确定出的大致流动方向进行验证。
93.步骤s203可以具体为：基于三维地震资料提取侵入岩体邻近层的相干属性图；基于相干属性图确定出侵入岩体的邻近层中断裂体系的展布特征。
94.其中，侵入岩体侵入时会在侵入岩体上部形成断裂体系，侵入岩体侵入时引起岩床上覆地层变形形成高陡裂隙，同时在岩床边缘形成逆断层，高陡裂隙和逆断层构成火成岩相关构造，也可称为次级断裂体系；侵入岩体的下部为走滑断裂体系。参见图3所述，图3示出了火成岩相构造示意图。通常，高陡裂隙和逆断层的走向接近垂直，逆断层和侵入岩体的流动方向垂直。在一些实施例中，基于高陡裂隙和逆断层接近垂直的组合方式，还可以验证断裂体系是由侵入岩体的侵入引起的。在本发明实施例中，断裂体系可以包括位于目标岩床模型上部的火成岩相关构造即次级断裂体系，基于高精度的三维地震资料提取侵入岩体邻近层的相干属性图，可以为基于高精度的三维地震资料提取侵入岩体上部邻近层的相干属性图，基于相干属性图可以确定出侵入岩体的火成岩相关构造，其中，火成岩相构造可以包括逆断层和高陡裂隙，逆断层发育在岩床模型的倾斜岩脉之上，逆断层的展布特征受控于倾斜岩脉的展布，高陡裂隙发育在岩床之上，通常对应岩脉错断的地方。基于逆断层的展布方向和侵入岩体的流动方向垂直的规律，可以精确确定出侵入岩体的展布方位。
95.进而，在本发明实施例中，步骤s204可以为：基于目标岩床模型的三维空间形态确定侵入岩体的第一流动方向；根据逆断层和侵入岩体的流动方向垂直的特征，确定侵入岩体的第二流动方向；通过结合第一流动方向和第二流动方向确定出侵入岩体的流动方向。其中，第一流动方向可以为侵入岩体的大致流动方向，第二流动方向可以为展布方位。作为一示例，基于目标三维空间形态如指状展布形态可以大致识别出侵入岩体的第一流动方向为自北向南流动，基于目标岩床模型的邻近层中断裂体系的展布特征可以识别出侵入岩体的第二流动方向为与逆断层垂直，综合第一流动方向和第二流动方向即可确定出侵入岩体的流动方向为向南流动过程中与逆断层垂直，从而可以精确确定出侵入岩体的流动方向。
96.以上为本发明实施例二提供的一种侵入岩体的流动方向确定方法，通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的时间地震反射特征确定出地震反射特征最清晰的岩床模型，作为目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态，从而通过选择最具有代表性的岩床模型作为目标岩床模型，可以有效提高识别侵入岩体流动方向的准确性。另外，基于高精度三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，基于三维空间形态以及断裂体系展布特征和侵入岩体流动的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
97.在本发明实施例中，为了提高确定侵入岩体流动方向的精确度，还可以精细化目
标岩床模型的三维空间形态，具体请参见以下实施例三中的描述。
98.实施例三
99.参见图5所示，图5示出了本发明实施例三提供的一种侵入岩体的流动方向确定方法流程示意图，其包括：
100.步骤s301：识别三维区内的侵入岩体以及识别出各时间地震剖面内侵入岩体的岩床模型。
101.步骤s302：获取不同时间地震剖面内岩床模型的地震反射特征，确定出地震反射特征最清晰的岩床模型作为目标岩床模型。
102.步骤s303：对目标岩床模型进行地震层位追踪。
103.步骤s304：对目标岩床模型的解释层位进行插值加密处理，以获取插值加密处理后目标岩床模型的三维空间形态。
104.步骤s305：基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，断裂体系包括位于侵入岩体上层的逆断层和高陡裂隙。
105.步骤s306：基于目标岩床模型的三维空间形态和断裂体系的展布特征，确定侵入岩体的流动方向。
106.在本发明实施例中，步骤s301可以具体为，获取三维区内的三维地震资料并从三维地震资料中提取时间地震剖面；基于时间地震剖面，根据侵入岩体的典型地震反射特征，识别三维区内的侵入岩体。其中，为了得到更准确、清晰的时间地震剖面，可以获取三维区内的高精度三维地震资料。
107.基于时间地震剖面，根据侵入岩体的典型地震反射特征识别三维区内的侵入岩体可以为，根据地震反射特征出现振幅反射异常如振幅增强且终止突然的特征，从时间地震剖面中识别出侵入岩体。在另一些实施例中，还可以基于高精度三维地震资料提取深度地震剖面。
108.另外，在步骤s301中识别各时间地震剖面内侵入岩体的岩床模型可以为：基于发育要素，识别出各时间地震剖面内侵入岩体的岩床模型。在一些实施例中，基于发育要素识别出各时间地震剖面内侵入岩体的岩床模型可以为，基于内部岩床、倾斜岩脉和外部岩床，识别出各时间地震剖面内侵入岩体的岩床模型。其中，内部岩床、倾斜岩脉和外部岩床位可以为碟形岩床的主要发育要素，基于碟形岩床的主要发育要素，可以从已识别出的侵入岩体中确定出碟形岩床，即岩床模型可以为碟形岩床，具体的，可以根据时间地震剖面中侵入岩体的剖面几何形态和碟形岩床的主要发育要素进行对比，从而识别出侵入岩体中的碟形岩床。在另一些实施例中，还可以基于碟形岩床完整发育模式中的发育要素识别出各时间地震剖面内侵入岩体的岩床模型，其中，碟形岩床完整发育模式中的发育要素可以包括内部岩床、倾斜岩脉、外部岩床和供给岩脉。
109.在本发明实施例中，步骤s302可以具体为，基于三维地震资料，统计不同时间地震剖面内岩床模型的地震反射特征，通过确定出地震反射特征最清晰的岩床模型，作为目标岩床模型。
110.在步骤s303中，可以采用本领域常规的方法对目标岩床模型进行地震层位追踪。
111.在步骤s304中可以通过对目标岩床模型的解释层位进行插值加密，提高解释精度，同时也有利于获取更精细的目标岩床模型三维空间形态。
112.基于目标岩床模型的三维空间形态可以确定出侵入岩体的大致流动方向，参见图2所示，图2示出了碟形岩床“指状”的空间展布示意图，侵入岩体自深部沿着垂向断裂通道向上涌入，在合适的地层界面侵入并在偏离垂向方向的侧向方向上发生流动。在侧向方向上流动过程中，由于流体密度差异、地层岩性和物性差异，侵入岩体并不能沿同一界面均匀铺开，而是逐渐分散成指状流动。因此，可以基于空间展布情况大致识别出，侵入流体的流动方向大致是指向指状分叉的方向。
113.在另一些实施例中，还可以基于三维地震资料，获取不同时间地震剖面内岩床模型的发育规律参量，其中，发育规律参量可以包括侵入层位、发育数量和上覆地层变形强度中的至少一种。通过获取岩床模型的发育规律参量可以获取三维区内侵入岩体的发育规律，侵入岩体的发育规律和侵入岩体的流动方向相关，例如，侵入岩体可沿多个界面侵入，形成垂向叠置特征，沿侵入岩体的流动方向，岩床模型的数量增多、最高侵入层位增高并且上覆地层变形强度增加。从而，基于发育可以确定出侵入岩体的大致流动方向，进一步可实现对基于目标岩床模型的三维空间形态确定出的大致流动方向进行验证。
114.步骤s305可以具体为：基于三维地震资料提取侵入岩体邻近层的相干属性图；基于相干属性图确定出侵入岩体的邻近层中断裂体系的展布特征。
115.其中，侵入岩体侵入时会在侵入岩体上部形成断裂体系，侵入岩体侵入时引起岩床上覆地层变形形成高陡裂隙，同时在岩床边缘形成逆断层，高陡裂隙和逆断层构成火成岩相关构造，也可称为次级断裂体系；侵入岩体的下部为走滑断裂体系。参见图3所述，图3示出了火成岩相构造示意图。通常，高陡裂隙和逆断层的走向接近垂直，逆断层和侵入岩体的流动方向垂直。在一些实施例中，基于高陡裂隙和逆断层接近垂直的组合方式，还可以验证断裂体系是由侵入岩体的侵入引起的。在本发明实施例中，断裂体系可以包括位于目标岩床模型上部的火成岩相关构造即次级断裂体系，基于高精度的三维地震资料提取侵入岩体邻近层的相干属性图，可以为基于高精度的三维地震资料提取侵入岩体上部邻近层的相干属性图，基于相干属性图可以确定出侵入岩体的火成岩相关构造，其中，火成岩相构造可以包括逆断层和高陡裂隙，逆断层发育在岩床模型的倾斜岩脉之上，逆断层的展布特征受控于倾斜岩脉的展布，高陡裂隙发育在岩床之上，通常对应岩脉错断的地方。基于逆断层的展布方向和侵入岩体的流动方向垂直的规律，可以精确确定出侵入岩体的展布方位。
116.进而，在本发明实施例中，步骤s306可以为：基于目标岩床模型的三维空间形态确定侵入岩体的第一流动方向；根据逆断层和侵入岩体的流动方向垂直的特征，确定侵入岩体的第二流动方向；通过结合第一流动方向和第二流动方向确定出侵入岩体的流动方向。其中，第一流动方向可以为侵入岩体的大致流动方向，第二流动方向可以为展布方位。作为一示例，基于目标三维空间形态如指状展布形态可以大致识别出侵入岩体的第一流动方向为自北向南流动，基于目标岩床模型的邻近层中断裂体系的展布特征可以识别出侵入岩体的第二流动方向为与逆断层垂直，综合第一流动方向和第二流动方向即可确定出侵入岩体的流动方向为向南流动过程中与逆断层垂直，从而可以精确确定出侵入岩体的流动方向。
117.以上为本发明实施例三提供的一种侵入岩体的流动方向确定方法，通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的时间地震反射特征确定出地震反射特征最清晰的岩床模型作为目标岩床模型，通过对目标岩床模型进行地震层位追踪以及对目标岩床模型的解释层位进行插值加密处理，可以提高解释精度，有利于获取更精
细的目标岩床模型三维空间形态。通过选择最具有代表性的岩床模型作为目标岩床模型，可以有效提高识别侵入岩体流动方向的准确性。另外，基于高精度三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，基于三维空间形态以及断裂体系展布特征和侵入岩体流动的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
118.基于本发明实施例提供的侵入岩体的流动方向确定方法，本发明还提供了一具体实施例，请参见实施例四。
119.实施例四
120.以塔里木盆地顺托果勒低隆地区三维覆盖区域为例，对塔里木盆地顺托果勒低隆地区三维覆盖区域进行侵入岩体精细识别、侵入岩体发育特征分析、侵入岩体三维空间形态雕刻、侵入岩体的邻近层中断裂体系解析，综合判识侵入岩体的流动方向。
121.精细识别侵入岩体：基于高精度三维地震资料提取时间地震剖面，根据侵入岩体典型地震反射特征(强振幅异常反射且终止突然)识别三维区内侵入岩体。
122.其中通过精细解释三维区内侵入岩体的剖面几何形态发现侵入岩体可以顺层侵入形成平坦岩床，也可以斜切地层形成倾斜岩脉，具体可参见图6a和图6b所示，图6a示出了本发明实施例提供的顺北地区三维区典型时间地震剖面示意图；图6b示出了本发明实施例提供的顺北地区三维区精细识别获得的侵入岩体时间地震剖面示意图，其中，在侵入岩体可以分布在多个地层如t
74
、t
70
、t
63
和t
60
，即侵入岩体可沿多个界面侵入。
123.识别侵入岩床模型(碟形岩床)：在基于高精度三维地震资料精细识别出侵入岩体特征的基础上，根据碟形岩床发育的主要发育要素，即内部岩床、倾斜岩脉和外部岩床或者根据碟形岩床完整发育的发育要素，即内部岩床、倾斜岩脉、外部岩床和供给岩脉(外部岩床和供给岩脉通常较少发育)，在已识别出的各种几何形态的侵入岩体中识别出碟形岩床。具体可参见图7a和图7b所示，图7a示出了碟形岩床完整发育模式的示意图；图7b示出了本发明实施例提供的顺北地区三维区碟形岩床时间地震剖面示意图。
124.明确三维区侵入岩体变化规律：基于高精度三维地震资料，统计不同时间地震剖面碟形岩床的发育数量、侵入层位和上覆地层变形强度，明确三维区南北向和东西向侵入岩体的发育规律。侵入岩体的发育规律和侵入岩体的流动方向相关，从而，基于发育规律可以确定出侵入岩体的大致流动方向，进一步可实现对基于目标岩床模型的三维空间形态确定出的大致流动方向进行验证。
125.参见图8所示，图8示出了顺北地区三维区内北西-东南向侵入岩体发育特征示意图。其中，图8中的左图是选择顺北地区三维区内的t
74
层平面图，基于t
74
层平面图中的a-a’线得到图8中的右图，即三维区内侵入岩体发育特征示意图。由图8中右图可知，自北向南，碟形岩床的数量增多、最高侵入层位增高且上覆地层变形强度增加。
126.雕刻侵入岩体三维空间形态：基于高精度地震资料获取到反射特征最清晰的碟形岩床作为目标碟形岩床，并对目标碟形岩床进行地震层位追踪。另外，还可以对目标碟形岩床的解释层位做插值加密处理，以提高解释精度，最后对目标碟形岩床进行立体显示，基于目标碟形岩床可以得到侵入岩体的三维空间形态。参见图9a和图9b所示，图9a示出了本发
明实施例提供的顺北地区三维区碟形岩床立体显示图；图9b示出了沿图9a中a-a’得到的“指状”展布时间地震剖面示意图，从图9a中可以看出侵入岩体整体具有北低南高的特征，空间上非均匀铺开，由北向南呈“指状”散开，而“指状”剖面上表现为多个孤立、错断的岩床。
127.侵入岩体的邻近层中断裂体系解释：基于高精度三维地震资料提取侵入岩体邻近层主要界面的相干属性图，精细解析不同层系断裂体系。
128.对于侵入岩体的上部和下部断裂体系具有分层变形特征，侵入岩体的下部为走滑断裂体系，上部为火成岩相关构造即次级断裂，次级断裂可以包括逆断层和高陡裂隙，参见图10a和图10b所示，图10a示出了本发明实施例提供的顺北地区三维区t
60
界面相干属性图；图10b示出了本发明实施例提供的顺北地区三维区t
60
界面火成岩相关构造精细解释示意图，其中，粗线为逆断层，细线为高陡裂隙，在平面上，逆断层和高陡裂隙的走向近乎垂直，形成棋盘格状交错展布特征。
129.参见图11a所示，图11a示出了沿图10b中a-a’得到的顺北地区三维区火成岩相关逆断层典型时间地震剖面示意图，其中，图11a中的左图为沿a-a’得到的时间地震剖面示意图，图11a中的右图为对左图精细解析后得到的逆断层时间地震剖面示意图，可以看出逆断层发育在倾斜岩脉上，其展布特征受控于倾斜岩脉的展布。
130.参见图11b所示，图11b示出了沿图10b中b-b’得到的顺北地区三维区火成岩相关高陡裂隙典型时间地震剖面示意图，其中，图11b中的左图为沿b-b’得到的时间地震剖面示意图，图11b中的右图为对左图精细解析后得到的高陡裂隙时间地震剖面示意图，可以看出高陡裂隙发育在岩床上，通常对应岩脉错断的地方。
131.基于侵入岩体三维空间形态和火成岩相关构造解析，综合判断侵入岩体的流动方向。
132.根据碟形岩床的“指状”展布形态，可以确认出侵入岩体的大致流动方向为自北向南；根据逆断层和高陡裂隙近乎垂直的组合，可以确定出逆断层是由侵入岩体的侵入引起，进而基于逆断层和侵入岩体的流动方向垂直的特征，结合碟形岩床的“指状”展布形态，可以确定出侵入岩体自北向南流动，并在向南流动的过程中与逆断层垂直，从而能够准确的确定出侵入岩体的流动方向，参见图12所示，图12示出了本发明实施例提供的顺北地区三维区侵入岩体流动方向示意图。
133.从而基于本发明实施例提供的侵入岩体的流动方向确定方法，通过精细识别与刻画侵入岩体的三维空间形态，并基于三维空间形态明确侵入岩体的展布规律；以及，基于高精度三维地震资料提取相干属性切片，分析侵入岩体上、下断裂体系的发育特征及其差异；最后根据侵入岩体的展布规律和流动方向的相关性，结合侵入岩体上断裂隙中逆断层和侵入岩体流动方向的相关性，可以有效确定出侵入岩体自深部上涌后的侧向流动方向。
134.最后通过二维测线和a井相关数据，发现在三维区北侧存在巨型侵入体玛纳岩体，该玛纳岩体可向南提供岩浆，进一步验证了本发明提供的侵入岩体流动方向确定方法的准确性。参见图13a和图13b所示，图13a示出了三维区玛纳侵入岩体分布位置示意图；图13b示出了a井采集到的花岗斑岩图。
135.实施例五
136.本发明的另一方面，还提供了一种存储介质，存储介质中存储有计算机程序，计算
机程序被处理器执行时能够实现如上实施例一至实施例四中任意一个实施例所述的侵入岩体的流动方向确定方法。
137.以上描述的处理、功能、方法和/或软件可被记录、存储或固定在一个或多个计算机可读存储介质中，计算机可读存储介质包括程序指令，程序指令将被计算机实现，以使处理器执行所述程序指令。存储介质还可单独包括程序指令、数据文件、数据结构等，或者包括其组合。存储介质或程序指令可被计算机软件领域的技术人员具体设计和理解，或者，存储介质或指令对计算机软件领域的技术人员而言可以是公知和可用的。计算机可读介质的示例包括：磁性介质，例如硬盘、软盘和磁带；光学介质，例如，cdrom盘和dvd；磁光介质，例如，光盘；和硬件装置，具体被配置以存储和执行程序指令，例如，只读存储器(rom)、随机存取存储器(ram)、闪存等。程序指令的示例包括机器代码(例如，由编译器产生的代码)和包含高级代码的文件，可由计算机通过使用解释器来执行所述高级代码。所描述的硬件装置可被配置为用作一个或多个软件模块，以执行以上描述的操作和方法，反之亦然。另外，计算机可读存储介质可分布在联网的计算机系统中，可以分散的方式存储和执行计算机可读代码或程序指令。
138.通过应用本发明实施例提供的存储介质，可以实现与上述实施例一至实施例四中任意一个实施例相同的有益效果，可以通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的地震反射特征确定出目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态；另外，基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，基于三维空间形态以及断裂体系展布特征和侵入岩体流动的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
139.实施例六
140.本发明的另一方面，还提供了一种设备，参见图14所示，图14示出了本发明实施例提供的一种设备示意图。
141.该设备可以包括处理器10和存储器11，存储器11中存储有计算机程序，处理器10执行存储器11中存储的计算机程序时能够实现如上实施例一至实施例四中任意一个实施例所述的侵入岩体的流动方向确定方法。
142.需要说明的是，该设备可以包括一个或多个处理器10以及存储器11，处理器10和存储器11可以通过总线或者其他方式连接。存储器11作为一种非易失性计算机可读存储介质，可用于存储非易失性软件程序、非易失性计算机可执行程序以及模块。处理器10通过运行存储在存储器中的非易失性软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现如上所述的侵入岩体的流动方向确定方法。
143.通过应用本发明实施例提供的设备，可以实现与上述实施例一至实施例四中任意一个实施例相同的有益效果，可以通过识别三维区内的侵入岩体以及侵入岩体中的岩床模型，基于岩床模型的地震反射特征确定出目标岩床模型并对目标岩床模型进行追踪，以获取目标岩床模型的三维空间形态；另外，基于三维地震资料，确定出侵入岩体的邻近层中断裂体系的展布特征，最后通过结合目标岩床模型的三维空间形态和断裂体系的展布特征，
基于三维空间形态以及断裂体系展布特征和侵入岩体流动的相关性，综合确定出侵入岩体的流动方向。该方法可以有效确定深层侵入岩体的流动方向，对预测研究的三维区周边非地震资料覆盖区火成岩体展布特征具有重要意义，还可以为火成岩体勘探领域提供依据，有利于实现井位部署和井轨迹优化。
144.虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。