Method, storage medium and device for determining flow direction of intrusive rock mass

Abstract

The present invention discloses a method, storage medium, and device for determining the flow direction of an intrusive rock mass. The method identifies the intrusive rock mass and the bedrock model in the intrusive rock mass within a three-dimensional region, determines the target bedrock model based on the seismic reflection characteristics of the bedrock model, and tracks the target bedrock model to obtain the three-dimensional spatial morphology of the target bedrock model. In addition, based on three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow direction of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses, and can also provide a basis for the exploration of igneous rock masses, facilitating the implementation of well site deployment and well trajectory optimization.

Description

Technical Field

[0001] The present invention relates to the field of structural analysis, and in particular to a method for determining the flow direction of an intrusive rock mass, a storage medium and a device. Background Art

[0002] Magma intrusion is an important research topic in internal geological dynamics and one of the key mechanical mechanisms of crustal deformation. Research both domestically and internationally has shown that rock intrusion follows specific stress patterns, exhibits directional characteristics, and is closely related to the geometric pattern of the intrusive body.

[0003] Taking the Tarim Basin as an example, the platform-basin region of the Tarim Basin (such as the Shunbei area) has experienced frequent magmatic activity since the Paleozoic. Seismic and drilling data reveal widespread strike-slip fault systems and numerous igneous intrusions within the deep (Paleozoic) layers of the region. The presence of igneous intrusions complicates strike-slip fault imaging and reservoir prediction, hindering play evaluation and well placement. Therefore, establishing the developmental patterns of igneous intrusions in a strike-slip setting and identifying their lateral flow directions is crucial for selecting and evaluating deep marine basins and optimizing well trajectories.

[0004] At present, the technologies for determining the flow direction of intrusive rock masses can be summarized into the following two types:

[0005] 1) Use the changes in noble gas isotopes (He, Ne and Ar isotopes) in rock samples to determine the direction of rock intrusion.

[0006] 2) The intrusion direction of the rock mass can be determined based on the position of the center of the residual gravity anomaly at different scales and the position of the center of the exposed rock mass.

[0007] However, existing technologies for determining the flow direction of intrusive rocks are limited in scope, applicable only to areas with abundant rock samples or exposed rock masses. They are not suitable for determining the flow direction of intrusive rocks deep within basins, where rock samples are scarce. Furthermore, in deep basins, where steep strike-slip fault systems and abundant igneous rocks are prevalent, the presence of intrusive rocks complicates imaging of these systems and reservoir prediction. Consequently, identifying the flow direction of intrusive rocks can affect subsequent zone evaluation and well placement. Summary of the Invention

[0008] The technical problem to be solved by the present invention is: how to effectively determine the flow direction of deep intrusive rock masses, so as to predict the distribution characteristics of igneous rock masses in the areas not covered by seismic data around the three-dimensional area under study, and also to provide a basis for the field of igneous rock exploration and realize the optimization of well location and well trajectory.

[0009] In order to solve the above technical problems, the present invention provides a method for determining the flow direction of an intrusive rock mass, a storage medium and a device.

[0010] A first aspect of the present invention provides a method for determining the flow direction of an intrusive rock mass, comprising:

[0011] Identifying intrusive rock masses within a three-dimensional region and identifying bedrock models of the intrusive rock masses within each seismic profile;

[0012] Acquiring seismic reflection characteristics of the bedrock model in different seismic sections, determining a target bedrock model based on the seismic reflection characteristics, and performing seismic layer tracing on the target bedrock model to obtain a three-dimensional spatial shape of the target bedrock model;

[0013] Determining the distribution characteristics of a fault system in a layer adjacent to the intrusive rock mass based on three-dimensional seismic data, wherein the fault system includes reverse faults and high-steep fractures located in an upper layer of the intrusive rock mass;

[0014] The flow direction of the intrusive rock mass is determined based on the three-dimensional spatial form of the target rock bed model and the distribution characteristics of the fracture system.

[0015] In some embodiments, the seismic profile comprises a time seismic profile, and identifying an intrusive rock mass within a three-dimensional region comprises:

[0016] Acquire the three-dimensional seismic data within the three-dimensional region and extract the time seismic profile from the three-dimensional seismic data;

[0017] Based on the time seismic profile, the intrusive rock mass in the three-dimensional region is identified according to typical seismic reflection characteristics of the intrusive rock mass.

[0018] In some embodiments, identifying the bedrock model of the intrusive rock mass in each seismic profile includes:

[0019] Based on the development factors, the bedrock model of the intrusive rock body in each time seismic section is identified.

[0020] In some embodiments, the bedrock model includes a dished bedrock, and identifying the bedrock model of the intrusive rock body in each time seismic profile based on the developmental factors includes:

[0021] The bedrock model of the intrusive rock body in each time seismic section is identified based on the internal bedrock, the inclined dikes and the external bedrock.

[0022] In some embodiments, obtaining seismic reflection characteristics of the bedrock model in different seismic sections and determining a target bedrock model based on the seismic reflection characteristics includes:

[0023] The seismic reflection characteristics of the bedrock model in seismic sections at different times are obtained, and the bedrock model with the clearest seismic reflection characteristics is determined as the target bedrock model.

[0024] In some embodiments, after performing seismic layer tracing on the target bedrock model, the method further includes: performing interpolation encryption processing on the interpreted layers of the target bedrock model.

[0025] In some embodiments, determining the distribution characteristics of the fault system in the adjacent layer of the intrusive rock mass based on the three-dimensional seismic data includes:

[0026] Extracting a coherent attribute map of adjacent layers of the intrusive rock mass based on the three-dimensional seismic data;

[0027] The distribution characteristics of the fracture system in the adjacent layer of the intrusive rock mass are determined based on the coherent attribute map.

[0028] In some embodiments, determining the flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fracture system includes:

[0029] determining a first flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target rock bed model;

[0030] determining a second flow direction of the intrusive rock mass based on a feature that the reverse fault and the flow direction of the intrusive rock mass are perpendicular;

[0031] The flow direction of the intrusive rock mass is determined based on the first flow direction and the second flow direction.

[0032] A second aspect of the present invention provides a storage medium storing a computer program, which, when executed by a processor, can implement any of the above methods for determining the flow direction of an intrusive rock mass.

[0033] A third aspect of the present invention provides a device comprising a processor and a memory, wherein the memory stores a computer program, and when the processor executes the computer program stored in the memory, it can implement the method for determining the flow direction of an intrusive rock mass as described in any one of the above.

[0034] Compared with the prior art, one or more embodiments of the above solutions may have the following advantages or beneficial effects:

[0035] The method for determining the flow direction of an intrusive rock mass provided by the present invention is applied. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass within a three-dimensional area, a target bedrock model is determined based on the seismic reflection characteristics of the bedrock model and the target bedrock model is tracked to obtain the three-dimensional spatial morphology of the target bedrock model; in addition, based on the three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layer of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow direction of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses, which is of great significance for predicting the distribution characteristics of igneous rock masses in areas not covered by seismic data around the three-dimensional area under study. It can also provide a basis for the exploration of igneous rock masses, which is conducive to the realization of well site deployment and well trajectory optimization. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The scope of the present disclosure may be better understood by reading the following detailed description of exemplary embodiments in conjunction with the accompanying drawings, which include:

[0037] Figure 1 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in the first embodiment of the present invention is shown;

[0038] Figure 2 A schematic diagram showing the spatial distribution of the "fingers" of the dish-shaped bedrock is shown;

[0039] Figure 3 A schematic diagram of igneous rock-related structures is shown;

[0040] Figure 4 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in a second embodiment of the present invention is shown;

[0041] Figure 5 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in a third embodiment of the present invention is shown;

[0042] Figure 6a A schematic diagram of a typical time earthquake section of the Shunbei region in three dimensions provided by an embodiment of the present invention is shown; Figure 6b A schematic diagram of a time seismic profile of an intrusive rock mass obtained by fine identification of a three-dimensional area in Shunbei area provided by an embodiment of the present invention is shown;

[0043] Figure 7a A schematic diagram showing the complete development pattern of the dished sill; Figure 7b A schematic diagram of a time seismic cross section of a three-dimensional dished bed in Shunbei area provided by an embodiment of the present invention is shown;

[0044] Figure 8The schematic diagram of the development characteristics of the northwest-southeast intrusive rock mass in the three-dimensional area of ​​Shunbei area is shown;

[0045] Figure 9a It shows a stereoscopic display of a dish-shaped rock bed in the Shunbei area provided by an embodiment of the present invention; Figure 9b Shown along Figure 9a Schematic diagram of the “finger-like” temporal seismic profile obtained from A-A’;

[0046] Figure 10a The three-dimensional area T6 of Shunbei area provided by the embodiment of the present invention is shown. 0 Interface coherence property diagram; Figure 10b The three-dimensional area T6 of Shunbei area provided by the embodiment of the present invention is shown. 0 Schematic diagram of detailed interpretation of interface igneous rock related structures;

[0047] Figure 11a Shown along Figure 10b Schematic diagram of typical time seismic profile of igneous rock-related reverse faults in the Shunbei area obtained by A-A'; Figure 11b Shown along Figure 10b Schematic diagram of typical time seismic profile of high-steep fractures related to igneous rocks in the Shunbei area obtained from B-B';

[0048] Figure 12 A schematic diagram of the flow direction of intrusive rock mass in the three-dimensional area of ​​Shunbei region provided by an embodiment of the present invention is shown;

[0049] Figure 13a A schematic diagram of the distribution of the Mana intrusive rock mass in the three-dimensional area is shown; Figure 13b The granite porphyry map collected in Well A is shown;

[0050] Figure 14 A schematic diagram of a device provided by an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0051] In order to make the objectives, technical solutions and advantages of the present invention clearer, the implementation method of the present invention will be described in detail below with reference to the accompanying drawings and embodiments, so that the implementation process of how the present invention applies technical means to solve technical problems and achieve technical effects can be fully understood and implemented accordingly.

[0052] Magma intrusion is an important research topic in internal geological dynamics and one of the key mechanical mechanisms of crustal deformation. Research both domestically and internationally has shown that rock intrusion follows specific stress patterns, exhibits directional characteristics, and is closely related to the geometric pattern of the intrusive body.

[0053] Taking the Tarim Basin as an example, the platform-basin region of the Tarim Basin (such as the Shunbei area) has experienced frequent magmatic activity since the Paleozoic. Seismic and drilling data reveal widespread strike-slip fault systems and numerous igneous intrusions within the deep (Paleozoic) layers of the region. The presence of igneous intrusions complicates strike-slip fault imaging and reservoir prediction, hindering play evaluation and well placement. Therefore, establishing the developmental patterns of igneous intrusions in a strike-slip setting and identifying their lateral flow directions is crucial for selecting and evaluating deep marine basins and optimizing well trajectories.

[0054] At present, the technologies for determining the flow direction of intrusive rock masses can be summarized into the following two types:

[0055] 1) Use the changes in noble gas isotopes (He, Ne and Ar isotopes) in rock samples to determine the direction of rock intrusion.

[0056] 2) The intrusion direction of the rock mass can be determined based on the position of the center of the residual gravity anomaly at different scales and the position of the center of the exposed rock mass.

[0057] However, existing technologies for determining the flow direction of intrusive rocks are limited in scope, applicable only to areas with abundant rock samples or exposed rock masses. They are not suitable for determining the flow direction of intrusive rocks deep within basins, where rock samples are scarce. Furthermore, in deep basins, where steep strike-slip fault systems and abundant igneous rocks are prevalent, the presence of intrusive rocks complicates imaging of these systems and reservoir prediction. Consequently, identifying the flow direction of intrusive rocks can affect subsequent zone evaluation and well placement.

[0058] In view of this, the present invention provides a method for determining the flow direction of an intrusive rock mass. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass within a three-dimensional area, a target bedrock model is determined based on the seismic reflection characteristics of the bedrock model and the target bedrock model is tracked to obtain the three-dimensional spatial morphology of the target bedrock model. In addition, based on three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow direction of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses, which is of great significance for predicting the distribution characteristics of igneous rock masses in areas not covered by seismic data around the three-dimensional area being studied. It can also provide a basis for the exploration of igneous rock masses and is conducive to the implementation of well site deployment and well trajectory optimization.

[0059] It should be noted that the method for determining the flow direction of an intrusive rock mass provided by the present invention can be applied to determining the flow direction of an intrusive rock mass under the background of deep strike-slip in a basin. Please refer to the following description for details.

[0060] Example 1

[0061] See also Figure 1 As shown, Figure 1 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in the first embodiment of the present invention is shown, which includes:

[0062] Step S101: Identify the intrusive rock mass in the three-dimensional region and identify the bedrock model of the intrusive rock mass in each seismic profile.

[0063] Step S102: obtaining seismic reflection characteristics of the bedrock model in different seismic sections, determining a target bedrock model based on the seismic reflection characteristics, and performing seismic layer tracking on the target bedrock model to obtain a three-dimensional spatial shape of the target bedrock model.

[0064] Step S103: Based on the 3D seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. The fault system includes reverse faults and high-steep fractures located in the upper layers of the intrusive rock mass.

[0065] Step S104: Determine the flow direction of the intrusive rock mass based on the three-dimensional spatial shape of the target bedrock model and the distribution characteristics of the fault system.

[0066] In an embodiment of the present invention, step S101 may specifically include acquiring three-dimensional seismic data within a three-dimensional area and extracting a seismic profile from the three-dimensional seismic data; based on the seismic profile, identifying the intrusive rock mass within the three-dimensional area according to typical seismic reflection characteristics of the intrusive rock mass.

[0067] In some embodiments, the seismic profile may be a time seismic profile; in other embodiments, the seismic profile may also be a depth seismic profile. In the following description, the time seismic profile will be used as an example.

[0068] Based on the time seismic profile, the intrusive rock mass in the three-dimensional area can be identified according to the typical seismic reflection characteristics of the intrusive rock mass. The intrusive rock mass can be identified from the time seismic profile based on the amplitude reflection anomaly such as amplitude enhancement and sudden termination in the seismic reflection characteristics.

[0069] In addition, in step S101, the bedrock model of the intrusive rock mass in each seismic profile may be identified by: identifying the bedrock model of the intrusive rock mass in each time seismic profile based on development factors. In some embodiments, the bedrock model of the intrusive rock mass in each time seismic profile may be identified by: identifying the bedrock model of the intrusive rock mass in each time seismic profile based on internal bedrock, inclined dykes, and external bedrock. The internal bedrock, inclined dykes, and external bedrock may be the main development factors of the dished bedrock. Based on the main development factors of the dished bedrock, the dished bedrock can be determined from the identified intrusive rock mass, that is, the bedrock model may be a dished bedrock. Specifically, the dished bedrock in the intrusive rock mass can be identified by comparing the cross-sectional geometry of the intrusive rock mass in the time seismic profile with the main development factors of the dished bedrock. In other embodiments, the bedrock model of the intrusive rock body in each time seismic profile can also be identified based on the development elements in the complete development model of the dished bedrock, wherein the development elements in the complete development model of the dished bedrock may include internal bedrock, inclined dykes, external bedrock and supply dykes.

[0070] In an embodiment of the present invention, step S102 may specifically be to obtain the seismic reflection characteristics of the bedrock model in the seismic profile at different times based on the three-dimensional seismic data. According to the accuracy requirements, the seismic reflection characteristics that meet the preset requirements may be selected to determine the target bedrock model.

[0071] In some embodiments, seismic reflection characteristics of bedrock models in seismic profiles at different times may be obtained, and based on accuracy requirements, a bedrock model with the clearest seismic reflection characteristics may be selected as the target bedrock model.

[0072] In addition, seismic layer tracking can be performed on the determined target bedrock model, and the target bedrock model can be displayed in stereo to obtain the three-dimensional spatial form of the target bedrock model.

[0073] The approximate flow direction of the intrusive rock mass can be determined based on the three-dimensional spatial morphology of the target bedrock model. Figure 2 As shown, Figure 2 A schematic diagram of the spatial distribution of the "fingers" of the dished bedrock is shown. The intrusive rock mass surges upward from the depths along vertical fault channels, intruding at suitable stratum interfaces and flowing laterally, deviating from the vertical. During this lateral flow, due to differences in fluid density, stratum lithology, and physical properties, the intrusive rock mass does not spread uniformly along the same interface, but instead gradually disperses into finger-like flows. Therefore, based on the spatial distribution, it can be roughly identified that the intrusive fluid is flowing in the direction of the finger-like bifurcation.

[0074] In other embodiments, developmental parameters of the bedrock model within seismic profiles at different times can also be obtained based on three-dimensional seismic data. These developmental parameters can include at least one of the intrusive layer position, the number of layers developed, and the deformation intensity of the overlying strata. By obtaining the developmental parameters of the bedrock model, the developmental patterns of the intrusive rock mass within the three-dimensional region can be determined. The developmental patterns of the intrusive rock mass are related to the flow direction of the intrusive rock mass. For example, an intrusive rock mass may intrude along multiple interfaces, forming a vertically stacked feature. Along the flow direction of the intrusive rock mass, the number of bedrock models increases, the highest intrusive layer increases, and the deformation intensity of the overlying strata increases. Thus, based on the developmental patterns, the approximate flow direction of the intrusive rock mass can be determined, and further verification of the approximate flow direction determined based on the three-dimensional spatial morphology of the target bedrock model can be achieved.

[0075] Step S103 may specifically include: extracting a coherent attribute map of the adjacent layers of the intrusive rock mass based on the 3D seismic data; and determining the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass based on the coherent attribute map.

[0076] Among them, when the intrusive rock mass invades, a fault system will be formed in the upper part of the intrusive rock mass, which will cause the deformation of the overlying strata of the bedrock to form high-steep cracks, and at the same time, reverse faults will be formed at the edge of the bedrock. High-steep cracks and reverse faults constitute igneous rock related structures, which can also be called secondary fault systems; the lower part of the intrusive rock mass is a strike-slip fault system. Figure 3 As stated, Figure 3 A schematic diagram of igneous rock-related structures is shown. Typically, the strike of steep fractures and reverse faults is nearly perpendicular, and the reverse faults are perpendicular to the flow direction of the intrusive rock mass. In some embodiments, based on the combination of steep fractures and reverse faults being nearly perpendicular, it can also be verified that the fracture system is caused by the intrusion of the intrusive rock mass. In an embodiment of the present invention, the fracture system can include igneous rock-related structures located above the target bedrock model, namely a secondary fracture system. A coherent attribute map of the adjacent layers of the intrusive rock mass is extracted based on 3D seismic data. The coherent attribute map can be used to extract a coherent attribute map of the adjacent layers above the intrusive rock mass based on 3D seismic data. Based on the coherent attribute map, the igneous rock-related structures of the intrusive rock mass can be determined. The igneous rock-related structures can include reverse faults and steep fractures. Reverse faults develop above the inclined dykes of the bedrock model. The distribution characteristics of reverse faults are controlled by the distribution of the inclined dykes. Steep fractures develop above the bedrock, usually corresponding to the locations where the dykes are dislocated. Based on the rule that the distribution direction of reverse faults is perpendicular to the flow direction of the intrusive rock mass, the distribution orientation of the intrusive rock mass can be accurately determined.

[0077] Furthermore, in an embodiment of the present invention, step S104 may be: determining the first flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target bedrock model; determining the second flow direction of the intrusive rock mass based on the feature that the reverse fault and the flow direction of the intrusive rock mass are perpendicular; and determining the flow direction of the intrusive rock mass by combining the first flow direction and the second flow direction. The first flow direction may be the approximate flow direction of the intrusive rock mass, and the second flow direction may be the distribution orientation. As an example, based on the target three-dimensional spatial morphology, such as the finger-like distribution morphology, it can be roughly identified that the first flow direction of the intrusive rock mass is flowing from north to south. Based on the distribution characteristics of the fault system in the adjacent layer of the target bedrock model, it can be identified that the second flow direction of the intrusive rock mass is perpendicular to the reverse fault. By combining the first flow direction and the second flow direction, it can be determined that the flow direction of the intrusive rock mass is perpendicular to the reverse fault during the southward flow, thereby accurately determining the flow direction of the intrusive rock mass.

[0078] The above is a method for determining the flow direction of an intrusive rock mass provided in Example 1 of the present invention. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass in a three-dimensional area, the target bedrock model is determined based on the seismic reflection characteristics of the bedrock model and the target bedrock model is tracked to obtain the three-dimensional spatial morphology of the target bedrock model; in addition, based on the three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses, which is of great significance for predicting the distribution characteristics of igneous rock masses in areas not covered by seismic data around the three-dimensional area under study. It can also provide a basis for the field of igneous rock exploration, which is conducive to the realization of well site deployment and well trajectory optimization.

[0079] In the embodiment of the present invention, in order to improve the accuracy of determining the flow direction of the intrusive rock mass, the most representative bedrock model may be selected as the target bedrock model. For details, please refer to the description in the following embodiment 2.

[0080] Example 2

[0081] See also Figure 4 As shown, Figure 4 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in a second embodiment of the present invention is shown, which includes:

[0082] Step S201: Identify the intrusive rock mass in the three-dimensional region and identify the bedrock model of the intrusive rock mass in each time seismic profile.

[0083] Step S202: Obtain the seismic reflection characteristics of the bedrock model in the seismic profile at different times, determine the bedrock model with the clearest seismic reflection characteristics as the target bedrock model, and perform seismic layer tracking on the target bedrock model to obtain the three-dimensional spatial shape of the target bedrock model.

[0084] Step S203: Based on the 3D seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. The fault system includes reverse faults and high-steep fractures located in the upper layers of the intrusive rock mass.

[0085] Step S204: Determine the flow direction of the intrusive rock mass based on the three-dimensional spatial shape of the target bedrock model and the distribution characteristics of the fault system.

[0086] In an embodiment of the present invention, step S201 may specifically include acquiring 3D seismic data within the 3D region and extracting a time seismic profile from the 3D seismic data; based on the time seismic profile, identifying the intrusive rock mass within the 3D region according to the typical seismic reflection characteristics of the intrusive rock mass. To obtain a more accurate and clear time seismic profile, high-precision 3D seismic data within the 3D region may be acquired.

[0087] Based on the time seismic profile, identifying the intrusive rock mass within the three-dimensional region based on the typical seismic reflection characteristics of the intrusive rock mass can be performed by identifying the intrusive rock mass from the time seismic profile based on the presence of amplitude reflection anomalies, such as amplitude enhancement and abrupt termination. In other embodiments, a depth seismic profile can also be extracted based on high-precision three-dimensional seismic data.

[0088] In addition, in step S201, the bedrock model of the intrusive rock mass in each time seismic profile may be identified based on development factors. In some embodiments, the bedrock model of the intrusive rock mass in each time seismic profile may be identified based on internal bedrock, inclined dykes, and external bedrock. The internal bedrock, inclined dykes, and external bedrock may be the main development factors of the dished bedrock. Based on the main development factors of the dished bedrock, the dished bedrock can be determined from the identified intrusive rock mass, that is, the bedrock model may be a dished bedrock. Specifically, the dished bedrock in the intrusive rock mass may be identified by comparing the cross-sectional geometry of the intrusive rock mass in the time seismic profile with the main development factors of the dished bedrock. In other embodiments, the bedrock model of the intrusive rock body in each time seismic profile can also be identified based on the development elements in the complete development model of the dished bedrock, wherein the development elements in the complete development model of the dished bedrock may include internal bedrock, inclined dykes, external bedrock and supply dykes.

[0089] In an embodiment of the present invention, step S202 may specifically be to obtain the seismic reflection characteristics of the bedrock model in the seismic profile at different times based on the three-dimensional seismic data, and determine the bedrock model with the clearest seismic reflection characteristics as the target bedrock model.

[0090] In addition, seismic layer tracking can be performed on the determined target bedrock model, and the target bedrock model can be displayed in stereo to obtain the three-dimensional spatial form of the target bedrock model.

[0091] The approximate flow direction of the intrusive rock mass can be determined based on the three-dimensional spatial morphology of the target bedrock model. Figure 2 As shown, Figure 2 A schematic diagram shows the spatial distribution of the "fingers" of the dished bedrock. The intrusive rock mass surges upward from the depths along vertical fault channels, intruding at suitable stratum interfaces and flowing laterally, deviating from the vertical direction. During this lateral flow, due to differences in fluid density, stratum lithology, and physical properties, the intrusive rock mass does not spread uniformly along the same interface, but instead gradually disperses into finger-like flows. Therefore, based on the spatial distribution, it can be roughly identified that the intrusive fluid is flowing in the direction of the finger-like bifurcation.

[0092] In other embodiments, developmental parameters of the bedrock model within seismic profiles at different times can also be obtained based on three-dimensional seismic data. These developmental parameters can include at least one of the intrusive layer position, the number of layers developed, and the deformation intensity of the overlying strata. By obtaining the developmental parameters of the bedrock model, the developmental patterns of the intrusive rock mass within the three-dimensional region can be determined. The developmental patterns of the intrusive rock mass are related to the flow direction of the intrusive rock mass. For example, an intrusive rock mass may intrude along multiple interfaces, forming a vertically stacked feature. Along the flow direction of the intrusive rock mass, the number of bedrock models increases, the highest intrusive layer increases, and the deformation intensity of the overlying strata increases. Thus, based on the developmental patterns, the approximate flow direction of the intrusive rock mass can be determined, and further verification of the approximate flow direction determined based on the three-dimensional spatial morphology of the target bedrock model can be achieved.

[0093] Step S203 may specifically include: extracting a coherent attribute map of the adjacent layers of the intrusive rock mass based on the 3D seismic data; and determining the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass based on the coherent attribute map.

[0094] Among them, when the intrusive rock mass invades, a fault system will be formed in the upper part of the intrusive rock mass, which will cause the deformation of the overlying strata of the bedrock to form high-steep cracks, and at the same time, reverse faults will be formed at the edge of the bedrock. High-steep cracks and reverse faults constitute igneous rock related structures, which can also be called secondary fault systems; the lower part of the intrusive rock mass is a strike-slip fault system. Figure 3 As stated, Figure 3A schematic diagram of igneous rock facies structures is shown. Typically, high-steep fractures and reverse faults strike nearly perpendicularly, and reverse faults are perpendicular to the flow direction of the intrusive rock mass. In some embodiments, based on the combination of high-steep fractures and reverse faults that are nearly perpendicular, it is also possible to verify that the fracture system is caused by the intrusion of the intrusive rock mass. In an embodiment of the present invention, the fracture system may include igneous rock-related structures located above the target bedrock model, namely a secondary fracture system. A coherence attribute map of the layers adjacent to the intrusive rock mass is extracted based on high-precision 3D seismic data. The coherence attribute map can be used to extract a coherence attribute map of the layers adjacent to the upper intrusive rock mass based on high-precision 3D seismic data. Based on the coherence attribute map, the igneous rock-related structures of the intrusive rock mass can be determined. The igneous rock facies structures may include reverse faults and high-steep fractures. Reverse faults develop above the inclined dykes of the bedrock model, and the distribution characteristics of reverse faults are controlled by the distribution of the inclined dykes. High-steep fractures develop above the bedrock, typically corresponding to locations where the dykes are dislocated. Based on the rule that the distribution direction of reverse faults is perpendicular to the flow direction of the intrusive rock mass, the distribution orientation of the intrusive rock mass can be accurately determined.

[0095] Furthermore, in an embodiment of the present invention, step S204 may be: determining the first flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target bedrock model; determining the second flow direction of the intrusive rock mass based on the feature that the reverse fault and the flow direction of the intrusive rock mass are perpendicular; and determining the flow direction of the intrusive rock mass by combining the first flow direction and the second flow direction. The first flow direction may be the approximate flow direction of the intrusive rock mass, and the second flow direction may be the distribution orientation. As an example, based on the target three-dimensional spatial morphology, such as the finger-like distribution morphology, it can be roughly identified that the first flow direction of the intrusive rock mass is flowing from north to south. Based on the distribution characteristics of the fault system in the adjacent layer of the target bedrock model, it can be identified that the second flow direction of the intrusive rock mass is perpendicular to the reverse fault. By combining the first flow direction and the second flow direction, it can be determined that the flow direction of the intrusive rock mass is perpendicular to the reverse fault during the southward flow, thereby accurately determining the flow direction of the intrusive rock mass.

[0096] The above is a method for determining the flow direction of an intrusive rock mass provided by Example 2 of the present invention. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass within a three-dimensional area, the bedrock model with the clearest seismic reflection characteristics is determined based on the temporal seismic reflection characteristics of the bedrock model. The bedrock model is used as the target bedrock model and the target bedrock model is tracked to obtain the three-dimensional spatial morphology of the target bedrock model. By selecting the most representative bedrock model as the target bedrock model, the accuracy of identifying the flow direction of the intrusive rock mass can be effectively improved. In addition, based on high-precision three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses and is of great significance for predicting the distribution characteristics of igneous rock masses in areas not covered by seismic data in the three-dimensional area of ​​the research. It can also provide a basis for the exploration of igneous rock masses and facilitate the optimization of well location and well trajectory.

[0097] In the embodiment of the present invention, in order to improve the accuracy of determining the flow direction of the intrusive rock mass, the three-dimensional spatial form of the target rock bed model may be further refined. For details, please refer to the description in the following embodiment three.

[0098] Example 3

[0099] See also Figure 5 As shown, Figure 5 A schematic flow chart of a method for determining the flow direction of an intrusive rock mass provided in a third embodiment of the present invention is shown, which includes:

[0100] Step S301: Identify the intrusive rock mass in the three-dimensional region and identify the bedrock model of the intrusive rock mass in each time seismic profile.

[0101] Step S302: obtaining seismic reflection characteristics of bedrock models in seismic profiles at different times, and determining the bedrock model with the clearest seismic reflection characteristics as the target bedrock model.

[0102] Step S303: performing seismic layer tracking on the target bedrock model.

[0103] Step S304: performing interpolation and encryption processing on the interpreted layers of the target bedrock model to obtain the three-dimensional spatial form of the target bedrock model after the interpolation and encryption processing.

[0104] Step S305: Based on the 3D seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. The fault system includes reverse faults and high-steep fractures located in the upper layers of the intrusive rock mass.

[0105] Step S306: Determine the flow direction of the intrusive rock mass based on the three-dimensional spatial shape of the target bedrock model and the distribution characteristics of the fault system.

[0106] In an embodiment of the present invention, step S301 may specifically include acquiring 3D seismic data within the 3D region and extracting a time seismic profile from the 3D seismic data; based on the time seismic profile, identifying the intrusive rock mass within the 3D region according to the typical seismic reflection characteristics of the intrusive rock mass. To obtain a more accurate and clear time seismic profile, high-precision 3D seismic data within the 3D region may be acquired.

[0107] Based on the time seismic profile, identifying the intrusive rock mass within the three-dimensional region based on the typical seismic reflection characteristics of the intrusive rock mass can be performed by identifying the intrusive rock mass from the time seismic profile based on the presence of amplitude reflection anomalies, such as amplitude enhancement and abrupt termination. In other embodiments, a depth seismic profile can also be extracted based on high-precision three-dimensional seismic data.

[0108] In addition, in step S301, the bedrock model of the intrusive rock mass in each time seismic profile may be identified by: identifying the bedrock model of the intrusive rock mass in each time seismic profile based on development factors. In some embodiments, the bedrock model of the intrusive rock mass in each time seismic profile may be identified based on internal bedrock, inclined dykes, and external bedrock. The internal bedrock, inclined dykes, and external bedrock may be the main development factors of the dished bedrock. Based on the main development factors of the dished bedrock, the dished bedrock can be determined from the identified intrusive rock mass, that is, the bedrock model may be a dished bedrock. Specifically, the dished bedrock in the intrusive rock mass may be identified by comparing the cross-sectional geometry of the intrusive rock mass in the time seismic profile with the main development factors of the dished bedrock. In other embodiments, the bedrock model of the intrusive rock body in each time seismic profile can also be identified based on the development elements in the complete development model of the dished bedrock, wherein the development elements in the complete development model of the dished bedrock may include internal bedrock, inclined dykes, external bedrock and supply dykes.

[0109] In an embodiment of the present invention, step S302 may specifically include, based on three-dimensional seismic data, counting the seismic reflection characteristics of the bedrock models in seismic profiles at different times, and determining the bedrock model with the clearest seismic reflection characteristics as the target bedrock model.

[0110] In step S303, conventional methods in the art may be used to perform seismic layer tracking on the target bedrock model.

[0111] In step S304, the interpretation layer of the target bedrock model can be interpolated and encrypted to improve the interpretation accuracy, which is also conducive to obtaining a more refined three-dimensional spatial form of the target bedrock model.

[0112] The approximate flow direction of the intrusive rock mass can be determined based on the three-dimensional spatial morphology of the target bedrock model. Figure 2 As shown, Figure 2 A schematic diagram of the spatial distribution of the "fingers" of the dished bedrock is shown. The intrusive rock mass surges upward from the depths along vertical fault channels, intruding at suitable stratum interfaces and flowing laterally, deviating from the vertical. During this lateral flow, due to differences in fluid density, stratum lithology, and physical properties, the intrusive rock mass does not spread uniformly along the same interface, but instead gradually disperses into finger-like flows. Therefore, based on the spatial distribution, it can be roughly identified that the intrusive fluid is flowing in the direction of the finger-like bifurcation.

[0113] In other embodiments, the development law parameters of the bedrock model within seismic profiles at different times can also be obtained based on three-dimensional seismic data. The development law parameters can include at least one of the intrusive layer position, the number of development layers, and the deformation intensity of the overlying strata. By obtaining the development law parameters of the bedrock model, the development law of the intrusive rock mass within the three-dimensional region can be obtained. The development law of the intrusive rock mass is related to the flow direction of the intrusive rock mass. For example, the intrusive rock mass can intrude along multiple interfaces, forming a vertically stacked feature. Along the flow direction of the intrusive rock mass, the number of bedrock models increases, the highest intrusive layer increases, and the deformation intensity of the overlying strata increases. Therefore, based on the development, the approximate flow direction of the intrusive rock mass can be determined, and the approximate flow direction determined based on the three-dimensional spatial morphology of the target bedrock model can be further verified.

[0114] Step S305 may specifically include: extracting a coherent attribute map of the adjacent layers of the intrusive rock mass based on the 3D seismic data; and determining the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass based on the coherent attribute map.

[0115] Among them, when the intrusive rock mass invades, a fault system will be formed in the upper part of the intrusive rock mass, which will cause the deformation of the overlying strata of the bedrock to form high-steep cracks, and at the same time, reverse faults will be formed at the edge of the bedrock. High-steep cracks and reverse faults constitute igneous rock related structures, which can also be called secondary fault systems; the lower part of the intrusive rock mass is a strike-slip fault system. Figure 3 As stated, Figure 3A schematic diagram of igneous rock facies structures is shown. Typically, high-steep fractures and reverse faults strike nearly perpendicularly, and reverse faults are perpendicular to the flow direction of the intrusive rock mass. In some embodiments, based on the combination of high-steep fractures and reverse faults that are nearly perpendicular, it is also possible to verify that the fracture system is caused by the intrusion of the intrusive rock mass. In an embodiment of the present invention, the fracture system may include igneous rock-related structures located above the target bedrock model, namely a secondary fracture system. A coherence attribute map of the layers adjacent to the intrusive rock mass is extracted based on high-precision 3D seismic data. The coherence attribute map can be used to extract a coherence attribute map of the layers adjacent to the upper intrusive rock mass based on high-precision 3D seismic data. Based on the coherence attribute map, the igneous rock-related structures of the intrusive rock mass can be determined. The igneous rock facies structures may include reverse faults and high-steep fractures. Reverse faults develop above the inclined dykes of the bedrock model, and the distribution characteristics of reverse faults are controlled by the distribution of the inclined dykes. High-steep fractures develop above the bedrock, typically corresponding to locations where the dykes are dislocated. Based on the rule that the distribution direction of reverse faults is perpendicular to the flow direction of the intrusive rock mass, the distribution orientation of the intrusive rock mass can be accurately determined.

[0116] Furthermore, in an embodiment of the present invention, step S306 may be: determining the first flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target bedrock model; determining the second flow direction of the intrusive rock mass based on the feature that the reverse fault and the flow direction of the intrusive rock mass are perpendicular; and determining the flow direction of the intrusive rock mass by combining the first flow direction and the second flow direction. The first flow direction may be the approximate flow direction of the intrusive rock mass, and the second flow direction may be the distribution orientation. As an example, based on the target three-dimensional spatial morphology, such as the finger-like distribution morphology, it can be roughly identified that the first flow direction of the intrusive rock mass is flowing from north to south. Based on the distribution characteristics of the fault system in the adjacent layer of the target bedrock model, it can be identified that the second flow direction of the intrusive rock mass is perpendicular to the reverse fault. By combining the first flow direction and the second flow direction, it can be determined that the flow direction of the intrusive rock mass is perpendicular to the reverse fault during the southward flow, thereby accurately determining the flow direction of the intrusive rock mass.

[0117] The above is a method for determining the flow direction of an intrusive rock mass provided by Example 3 of the present invention. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass in a three-dimensional area, the bedrock model with the clearest seismic reflection characteristics is determined as the target bedrock model based on the temporal seismic reflection characteristics of the bedrock model. By tracking the seismic layer of the target bedrock model and interpolating and encrypting the interpretation layer of the target bedrock model, the interpretation accuracy can be improved, which is conducive to obtaining a more refined three-dimensional spatial form of the target bedrock model. By selecting the most representative bedrock model as the target bedrock model, the accuracy of identifying the flow direction of the intrusive rock mass can be effectively improved. In addition, based on high-precision three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial form of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass is comprehensively determined based on the correlation between the three-dimensional spatial form and the distribution characteristics of the fault system and the flow of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock masses, which is of great significance for predicting the distribution characteristics of igneous rock masses in areas not covered by seismic data around the three-dimensional area under study. It can also provide a basis for the exploration of igneous rock masses and is conducive to the optimization of well location and well trajectory.

[0118] Based on the method for determining the flow direction of an intrusive rock mass provided in an embodiment of the present invention, the present invention also provides a specific embodiment, please refer to Example 4.

[0119] Example 4

[0120] Taking the three-dimensional covered area of ​​Shuntoguole low-rise in the Tarim Basin as an example, the three-dimensional covered area of ​​Shuntoguole low-rise in the Tarim Basin was finely identified intrusive rock masses, the development characteristics of intrusive rock masses were analyzed, the three-dimensional spatial morphology of intrusive rock masses was carved, the fault system in the adjacent layers of the intrusive rock masses was analyzed, and the flow direction of the intrusive rock masses was comprehensively judged.

[0121] Fine identification of intrusive rock masses: Time seismic profiles are extracted based on high-precision 3D seismic data, and intrusive rock masses in the 3D area are identified based on their typical seismic reflection characteristics (strong amplitude abnormal reflection and sudden termination).

[0122] Through the detailed interpretation of the cross-sectional geometry of the intrusive rock mass in the three-dimensional area, it is found that the intrusive rock mass can intrude along the layer to form a flat bed, or it can cut obliquely into the stratum to form an inclined dyke. Figure 6a and Figure 6b As shown, Figure 6a A schematic diagram of a typical time earthquake section of the Shunbei region in three dimensions provided by an embodiment of the present invention is shown; Figure 6b The schematic diagram of the time seismic profile of the intrusive rock mass obtained by fine identification of the three-dimensional area in Shunbei area provided by the embodiment of the present invention is shown, wherein the intrusive rock mass can be distributed in multiple strata such as T7 4、T7 0 、T6 3 and T6 0 , that is, the intrusive rock mass can intrude along multiple interfaces.

[0123] Identification of intrusive bed models (dish-shaped bed): Based on the fine identification of intrusive rock mass characteristics based on high-precision 3D seismic data, dish-shaped bed models are identified in various geometric forms of intrusive rock masses according to the main developmental elements of dish-shaped bed development, namely internal bed, inclined dyke and external bed, or according to the developmental elements of complete dish-shaped bed development, namely internal bed, inclined dyke, external bed and supply dyke (external bed and supply dyke are usually less developed). For details, please refer to Figure 7a and Figure 7b As shown, Figure 7a A schematic diagram showing the complete development pattern of the dished sill; Figure 7b A schematic diagram of a time seismic profile of a three-dimensional dish-shaped bedrock in Shunbei area provided by an embodiment of the present invention is shown.

[0124] Clarifying the evolution of intrusive rock masses in a 3D region: Based on high-precision 3D seismic data, the number of dished bedrocks, the intrusive horizons, and the deformation intensity of overlying strata along seismic profiles at different times are counted to clarify the development patterns of north-south and east-west intrusive rock masses in the 3D region. The development patterns of intrusive rock masses are correlated with their flow direction. Therefore, based on these development patterns, the approximate flow direction of the intrusive rock masses can be determined. This can then be used to verify the approximate flow direction determined based on the 3D spatial morphology of the target bedrock model.

[0125] See also Figure 8 As shown, Figure 8 The following is a schematic diagram showing the development characteristics of northwest-southeast intrusive rock bodies in the three-dimensional area of ​​Shunbei area. Figure 8 The left picture shows T7 in the Shunbei area. 4 Floor plan, based on T7 4 The A-A' line in the floor plan is obtained Figure 8 The right figure in the figure is a schematic diagram of the development characteristics of intrusive rock in the three-dimensional area. Figure 8 As can be seen from the middle right figure, from north to south, the number of dished sills increases, the highest intrusion layer increases, and the deformation intensity of the overlying strata increases.

[0126] Sculpting the 3D spatial morphology of intrusive rock masses: Based on high-precision seismic data, the dished bed with the clearest reflection characteristics is obtained as the target dished bed, and seismic horizon tracking is performed on the target dished bed. In addition, the interpretation horizon of the target dished bed can be interpolated and encrypted to improve the interpretation accuracy. Finally, the target dished bed is displayed in stereo. Based on the target dished bed, the 3D spatial morphology of the intrusive rock mass can be obtained. Figure 9a and Figure 9bAs shown, Figure 9a It shows a stereoscopic display of a dish-shaped rock bed in the Shunbei area provided by an embodiment of the present invention; Figure 9b Shown along Figure 9a The schematic diagram of the “finger-like” time seismic profile obtained from A-A' in Figure 9a It can be seen from the figure that the intrusive rock mass as a whole has the characteristics of being low in the north and high in the south, and is unevenly spread in space, spreading out in a "finger-like" manner from north to south. The "finger-like" section appears as multiple isolated and faulted sills.

[0127] Interpretation of fault systems in adjacent layers of intrusive rock masses: Based on high-precision 3D seismic data, coherent attribute maps of the main interfaces of adjacent layers of intrusive rock masses are extracted to finely analyze the fault systems of different strata.

[0128] The upper and lower fault systems of the intrusive rock mass have layered deformation characteristics. The lower part of the intrusive rock mass is a strike-slip fault system, and the upper part is an igneous rock-related structure, namely a secondary fault. The secondary fault can include reverse faults and high-steep fractures. Figure 10a and Figure 10b As shown, Figure 10a The three-dimensional area T6 of Shunbei area provided by the embodiment of the present invention is shown. 0 Interface coherence property diagram; Figure 10b The three-dimensional area T6 of Shunbei area provided by the embodiment of the present invention is shown. 0 Schematic diagram of detailed interpretation of interface igneous rock-related structures, where thick lines represent reverse faults and thin lines represent steep fractures. In the plane, the strikes of reverse faults and steep fractures are almost perpendicular, forming a checkerboard-like staggered distribution feature.

[0129] See also Figure 11a As shown, Figure 11a Shown along Figure 10b A-A' in the middle shows a typical time seismic profile of the 3D igneous rock-related reverse faults in the Shunbei area, where: Figure 11a The left figure is a schematic diagram of the time seismic section obtained along A-A'. Figure 11a The right picture in the figure is a schematic diagram of the reverse fault time seismic profile obtained after a detailed analysis of the left picture. It can be seen that the reverse fault develops on the inclined dyke, and its distribution characteristics are controlled by the distribution of the inclined dyke.

[0130] See also Figure 11b As shown, Figure 11b Shown along Figure 10b Schematic diagram of typical time seismic profile of high-steep fractures related to igneous rocks in the Shunbei area obtained by B-B', among which: Figure 11b The left figure is a schematic diagram of the time seismic profile obtained along B-B'. Figure 11bThe right picture in the figure is a schematic diagram of the time seismic profile of high-steep fractures obtained after detailed analysis of the left picture. It can be seen that high-steep fractures are developed on the bedrock, usually corresponding to the places where the dykes are broken.

[0131] Based on the three-dimensional spatial morphology of the intrusive rock mass and the analysis of the related structures of the igneous rock, the flow direction of the intrusive rock mass is comprehensively judged.

[0132] Based on the "finger-like" distribution of the dished bedrock, it can be confirmed that the approximate flow direction of the intrusive rock mass is from north to south. Based on the nearly vertical combination of reverse faults and high-steep fractures, it can be determined that the reverse faults are caused by the intrusion of the intrusive rock mass. Furthermore, based on the perpendicular characteristics of the reverse faults and the flow direction of the intrusive rock mass, combined with the "finger-like" distribution of the dished bedrock, it can be determined that the intrusive rock mass flows from north to south and is perpendicular to the reverse faults during the southward flow, thus accurately determining the flow direction of the intrusive rock mass. Figure 12 As shown, Figure 12 A schematic diagram of the flow direction of intrusive rock mass in the three-dimensional area of ​​Shunbei region provided by an embodiment of the present invention is shown.

[0133] Therefore, based on the method for determining the flow direction of an intrusive rock mass provided by an embodiment of the present invention, the three-dimensional spatial morphology of the intrusive rock mass is finely identified and depicted, and the distribution pattern of the intrusive rock mass is clarified based on the three-dimensional spatial morphology; and coherent attribute slices are extracted based on high-precision three-dimensional seismic data to analyze the development characteristics and differences of the upper and lower fault systems of the intrusive rock mass; finally, based on the correlation between the distribution pattern of the intrusive rock mass and the flow direction, combined with the correlation between the reverse faults in the fracture gap above the intrusive rock mass and the flow direction of the intrusive rock mass, the lateral flow direction of the intrusive rock mass after upwelling from the deep can be effectively determined.

[0134] Finally, through the two-dimensional survey line and the relevant data of Well A, it was found that there was a giant intrusive body, the Mana rock body, on the north side of the three-dimensional area. The Mana rock body can provide magma to the south, further verifying the accuracy of the method for determining the flow direction of the intrusive rock body provided by the present invention. Figure 13a and Figure 13b As shown, Figure 13a A schematic diagram of the distribution of the Mana intrusive rock mass in the three-dimensional area is shown; Figure 13b The figure shows the granite porphyry sampled from Well A.

[0135] Example 5

[0136] Another aspect of the present invention provides a storage medium storing a computer program. When the computer program is executed by a processor, the method for determining the flow direction of an intrusive rock mass as described in any one of the first to fourth embodiments above can be implemented.

[0137] The processes, functions, methods, and / or software described above may be recorded, stored, or fixed in one or more computer-readable storage media, which include program instructions that are implemented by a computer so that a processor executes the program instructions. The storage media may also include program instructions, data files, data structures, etc., alone or in combination. The storage media or program instructions may be specifically designed and understood by those skilled in the art of computer software, or the storage media or instructions may be generally known and available to those skilled in the art of computer software. Examples of computer-readable media include: magnetic media, such as hard disks, floppy disks, and magnetic tapes; optical media, such as CDROMs and DVDs; magneto-optical media, such as optical disks; and hardware devices specifically configured to store and execute program instructions, such as read-only memory (ROM), random access memory (RAM), flash memory, etc. Examples of program instructions include machine code (e.g., code generated by a compiler) and files containing higher-level code that can be executed by a computer using an interpreter. The described hardware devices may be configured to function as one or more software modules to perform the operations and methods described above, and vice versa. In addition, the computer-readable storage medium may be distributed among networked computer systems, so that computer-readable codes or program instructions may be stored and executed in a decentralized manner.

[0138] By applying the storage medium provided by the embodiment of the present invention, the same beneficial effects as any one of the embodiments 1 to 4 above can be achieved. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass in a three-dimensional area, the target bedrock model is determined based on the seismic reflection characteristics of the bedrock model and the target bedrock model is tracked to obtain the three-dimensional spatial morphology of the target bedrock model; in addition, based on the three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layer of the intrusive rock mass are determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow of the intrusive rock mass, the flow direction of the intrusive rock mass is comprehensively determined. This method can effectively determine the flow direction of deep intrusive rock mass, which is of great significance for predicting the distribution characteristics of igneous rock mass in the non-seismic data coverage area around the three-dimensional area of ​​the research, and can also provide a basis for the exploration of igneous rock mass, which is conducive to the realization of well location deployment and well trajectory optimization.

[0139] Example 6

[0140] Another aspect of the present invention provides a device, see Figure 14 As shown, Figure 14 A schematic diagram of a device provided by an embodiment of the present invention is shown.

[0141] The device may include a processor 10 and a memory 11, wherein the memory 11 stores a computer program. When the processor 10 executes the computer program stored in the memory 11, the method for determining the flow direction of an intrusive rock mass as described in any one of the first to fourth embodiments above can be implemented.

[0142] It should be noted that the device may include one or more processors 10 and a memory 11, which may be connected via a bus or other means. Memory 11, as a non-volatile computer-readable storage medium, can be used to store non-volatile software programs, non-volatile computer executable programs, and modules. Processor 10 executes the non-volatile software programs, instructions, and modules stored in the memory to perform various functional applications and data processing of the device, thereby implementing the aforementioned method for determining the flow direction of an intrusive rock mass.

[0143] By applying the device provided by the embodiment of the present invention, the same beneficial effects as any one of the embodiments 1 to 4 above can be achieved. By identifying the intrusive rock mass and the bedrock model in the intrusive rock mass in a three-dimensional area, the target bedrock model can be determined based on the seismic reflection characteristics of the bedrock model and the target bedrock model can be tracked to obtain the three-dimensional spatial morphology of the target bedrock model; in addition, based on the three-dimensional seismic data, the distribution characteristics of the fault system in the adjacent layer of the intrusive rock mass can be determined. Finally, by combining the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system, the flow direction of the intrusive rock mass can be comprehensively determined based on the correlation between the three-dimensional spatial morphology and the distribution characteristics of the fault system and the flow of the intrusive rock mass. This method can effectively determine the flow direction of deep intrusive rock mass, which is of great significance for predicting the distribution characteristics of igneous rock mass in the non-seismic data coverage area around the three-dimensional area of ​​the research. It can also provide a basis for the exploration of igneous rock mass, which is conducive to the realization of well site deployment and well trajectory optimization.

[0144] Although the embodiments disclosed herein are as described above, the contents described herein are merely embodiments for facilitating understanding of the present invention and are not intended to limit the present invention. Any person skilled in the art may make any modifications and variations in the form and details of the embodiments without departing from the spirit and scope of the present invention. However, the scope of protection of the present invention shall remain subject to the scope defined by the appended claims.

Claims

1. A method for determining the flow direction of an intrusive rock mass, characterized in that: include: Identifying intrusive rock masses within a three-dimensional region and identifying bedrock models of the intrusive rock masses within each seismic profile; Acquiring seismic reflection characteristics of the bedrock model in different seismic sections, determining a target bedrock model based on the seismic reflection characteristics, and performing seismic layer tracing on the target bedrock model to obtain a three-dimensional spatial shape of the target bedrock model; Determining the distribution characteristics of a fault system in a layer adjacent to the intrusive rock mass based on three-dimensional seismic data, wherein the fault system includes reverse faults and high-steep fractures located in an upper layer of the intrusive rock mass; Determining the flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target bedrock model and the distribution characteristics of the fault system; The determining of the flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target rock bed model and the distribution characteristics of the fracture system includes: Determine the first flow direction of the intrusive rock mass based on the three-dimensional spatial morphology of the target rock bed model; According to the feature that the reverse fault and the flow direction of the intrusive rock mass are perpendicular, the second flow direction of the intrusive rock mass is determined; The flow direction of the intrusive rock mass is determined by combining the first flow direction and the second flow direction.

2. The method according to claim 1, characterized in that The seismic profile includes a time seismic profile, and the identifying of the intrusive rock mass in the three-dimensional region includes: Acquire the three-dimensional seismic data within the three-dimensional region and extract the time seismic profile from the three-dimensional seismic data; Based on the time seismic profile, the intrusive rock mass in the three-dimensional region is identified according to typical seismic reflection characteristics of the intrusive rock mass.

3. The method according to claim 2, characterized in that The identifying of the bedrock model of the intrusive rock mass in each seismic profile includes: Based on the development factors, the bedrock model of the intrusive rock body in each time seismic section is identified.

4. The method according to claim 3, characterized in that The bedrock model includes a dished bedrock, and the bedrock model of the intrusive rock mass in each time seismic profile is identified based on the development factor, including: The bedrock model of the intrusive rock body in each time seismic section is identified based on the internal bedrock, the inclined dikes and the external bedrock.

5. The method according to claim 1 or 4, characterized in that The obtaining of seismic reflection characteristics of the bedrock model in different seismic sections and determining a target bedrock model based on the seismic reflection characteristics includes: The seismic reflection characteristics of the bedrock model in seismic sections at different times are obtained, and the bedrock model with the clearest seismic reflection characteristics is determined as the target bedrock model.

6. The method according to claim 1, characterized in that After performing seismic layer tracing on the target bedrock model, the method further includes: performing interpolation encryption processing on the interpretation layer of the target bedrock model.

7. The method according to claim 1, characterized in that Determining the distribution characteristics of the fault system in the adjacent layers of the intrusive rock mass based on the three-dimensional seismic data includes: Extracting a coherent attribute map of adjacent layers of the intrusive rock mass based on the three-dimensional seismic data; The distribution characteristics of the fracture system in the adjacent layer of the intrusive rock mass are determined based on the coherent attribute map.

8. A storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method for determining the flow direction of an intrusive rock mass as claimed in any one of claims 1 to 7 can be implemented.

9. A device, characterized in that The device includes a processor and a memory, wherein a computer program is stored in the memory. When the processor executes the computer program stored in the memory, the method for determining the flow direction of an intrusive rock mass as described in any one of claims 1 to 7 can be implemented.

Images