Negative inversion structure cause analysis method, system, electronic device, and storage medium

Abstract

A method, system, electronic device, and storage medium for analyzing the genesis of negative inversion structures, the method comprising: constructing a geological generation model of negative inversion structures, wherein the geological generation model comprises two strata separated by a pre-existing plate fault; controlling the geological generation model accordingly according to the compression stage and the extension inversion stage, wherein the extension angles of the two strata in the extension inversion stage are adjustable; recording the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage; obtaining the geographical features of the area to be analyzed, matching the geographical features with the geometric and kinematic features, determining whether the area to be analyzed belongs to a negative inversion structure, and if so, obtaining the extension angle of the area to be analyzed in the extension inversion stage. The method can make a genetic analysis of the negative inversion basin in Linhe Depression, analyze its structural sedimentary features in different periods, and provide supplementary suggestions for the geological exploration of the area.

Description

Technical Field

[0001] The present invention relates to the technical field of structural geology, and in particular to a method, system, electronic equipment and storage medium for analyzing the genesis of a negative inversion structure. Background Art

[0002] Discrete element simulation (Distinct / Discrete Element Method for Modeling) is an effective method for studying tectonic dynamics. Discrete element is suitable for studying the rupture problem of non-continuous media, and has been widely used in studying large deformation problems in structural geology.

[0003] Current studies all use single-phase or multi-phase co-directional extension or co-directional compression simulation methods, but there is still a lack of research on the negative inversion structure of early compression superimposed on late extension. Summary of the invention

[0004] In view of the problems existing in the prior art, the present invention provides a method, system, electronic device and storage medium for analyzing the cause of a negative inversion structure.

[0005] The present invention provides a method for analyzing the causes of negative inversion structures, the method comprising:

[0006] Constructing a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault;

[0007] The geological generation model is controlled in sequence according to the compression stage and the extension and reversal stage, wherein the extension angle of the two strata in the extension and reversal stage is adjustable;

[0008] Record the geometric and kinematic characteristics of the geological generation model at different extension angles at each stage;

[0009] The geographical features of the area to be analyzed are obtained, and the geographical features are matched with the geometric and kinematic features to determine whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0010] According to a method for analyzing the causes of negative inversion structures provided by the present invention, the geological generation model includes:

[0011] A stress-bearing box body, wherein particles are filled inside the stress-bearing box body;

[0012] A first rigid base, a deformable base and a second rigid base which are connected in sequence and are arranged at the bottom of the force-bearing box;

[0013] A pre-existing plate fault disposed at the junction of the first rigid base and the deformable base and forming a first angle with the first rigid base;

[0014] Among them, the first rigid base and the pre-existing plate-like faults and the particles therein constitute a first stratum with a fixed position, and the second rigid base, the deformable base and the pre-existing plate-like faults and the particles therein constitute a second stratum that first squeezes toward the first stratum and then stretches in the opposite direction.

[0015] According to a method for analyzing the causes of negative inversion structures provided by the present invention, the property parameters of the filling particles in the load-bearing box are set to a porosity of 0.15 and a density of 2600.0 kg / m 3 , the particle size of spherical particles is 0.03m.

[0016] According to a method for analyzing the causes of negative inversion structures provided by the present invention, the angle between the pre-existing plate fault and the first rigid base is 60°.

[0017] According to a method for analyzing the genesis of a negative inversion structure provided by the present invention, the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage are recorded, including:

[0018] The vertical displacement analysis and horizontal displacement analysis are performed on the constituent particles of the geological generation model to form multiple zones.

[0019] According to a method for analyzing the causes of negative inversion structures provided by the present invention, a cap layer deposition stage is provided between the compression stage and the extension inversion stage.

[0020] The present invention also provides a negative inversion structure cause analysis system, the system comprising:

[0021] a construction module, the construction module constructing a geological generation model of a negative inversion structure, wherein the geological generation model includes two strata separated by a pre-existing plate fault;

[0022] A generation module, wherein the generation module controls the geological generation model in accordance with the compression stage and the extension and reversal stage in sequence, wherein the extension angles of the two strata in the extension and reversal stage are adjustable;

[0023] A recording module, which records the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage;

[0024] A judgment module, wherein the judgment module obtains the geographical features of the area to be analyzed, matches them with the geometric and kinematic features, and determines whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0025] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the negative inversion structure cause analysis method as described in any one of the above items are implemented.

[0026] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the negative inversion structure cause analysis method as described in any one of the above items are implemented.

[0027] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the steps of the negative inversion structure cause analysis method as described in any one of the above items are implemented.

[0028] The negative inversion structure genesis analysis method, system, electronic device and storage medium provided by the present invention can perform a genesis analysis on the negative inversion basin in Linhe Depression, analyze its structural sedimentary features in different periods, and provide supplementary suggestions for geological exploration in the area. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0030] Figure 1 A schematic diagram of the structure of a geological generation model disclosed in the present invention;

[0031] Figure 2 A top view of the initial state of a geological generation model disclosed by the present invention;

[0032] Figure 3 A deformation end top view of a geological generation model disclosed by the present invention;

[0033] Figure 4 This is a diagram of the deformation evolution process of the early compression stage of the negative inversion structure disclosed by the present invention;

[0034] Figure 5 It is a diagram of the total amount of accumulated displacement of particles during the evolution process of the early compression stage of the negative inversion structure disclosed by the present invention;

[0035] Figure 6 A diagram of the vertical displacement of particles during the evolution of the early compression stage of the negative inversion structure disclosed by the present invention;

[0036] Figure 7This is a diagram showing the evolution of the first experimental superposition angle 0° stretching stage disclosed by the present invention;

[0037] Figure 8 A comparison diagram of the simulation results of the anisotropic stretching experiment published by the present invention and the first experiment;

[0038] Fig. 9 A comparison diagram of vertical displacement between the anisotropic stretching experiment and the first experiment disclosed by the present invention;

[0039] Fig.10 A comparison diagram of the top surface structure of the anisotropic stretching experiment and the first experiment disclosed by the present invention;

[0040] Fig.11 A comparison diagram of the slice structure between the anisotropic stretching experiment and the first experiment disclosed by the present invention;

[0041] Fig.12 A schematic diagram of the three-dimensional structure of the Jilantai Sag in the Linhe Depression disclosed in the present invention;

[0042] Fig.13 A schematic diagram of a process for analyzing the cause of negative inversion structure provided by the present invention;

[0043] Fig.14 A schematic diagram of the structure of a negative inversion structure cause analysis system provided by the present invention;

[0044] Fig.15 A schematic diagram of the physical structure of an electronic device provided by the present invention;

[0045] Fig.16 for Figure 6 A local enlarged schematic diagram of the middle A;

[0046] Fig.17 for Figure 6 A partial enlarged schematic diagram of point B in the middle;

[0047] Fig.18 for Figure 6 A partial enlarged schematic diagram of point C in the middle.

[0048] Figure numerals: 1-granite gneiss, 2-schist, 3-conglomerate, 4-sandstone, 5-siltstone, 6-Archean, 7-Proterozoic, 8-lower member of Lisangou Formation, 9-upper member of Lisangou Formation, 10-Guyang Formation, 11-Paleogene, 12-Neogene, 13-pre-existing plate fault, 14-deformable basement, 15-first rigid basement, 16-second rigid basement, 17-thrust fault. DETAILED DESCRIPTION

[0049] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0050] The negative inversion structure cause analysis method provided in the embodiment of the present application is described in detail below through specific embodiments and application scenarios in conjunction with the accompanying drawings.

[0051] Embodiment 1

[0052] This embodiment provides a geological generation model, including a stress-bearing box, a pre-existing plate fault 13, a deformable base 14, a first rigid base 15, and a second rigid base 16. The pre-existing plate fault, the deformable base 14, the first rigid base 15, and the second rigid base 16 are connected to each other, and the box above the base is filled with spherical particles.

[0053] Figure 1 This is a schematic diagram of the structure of a geological generation model disclosed by the present invention, such as Figure 1 As shown in the figure, the initial state of the load-bearing box is set to a cuboid with a length of 8m, a width of 3m, and a height of 1m. The box is filled with 166,586 spherical particles, and the property parameters of the particles are set to porosity 0.15 and density 2600.0kg / m 3 , the particle size of the spherical particles is 0.03 m, and the particle size conforms to the Gaussian distribution. The detailed parameter settings of the geological generation model are shown in Table 1 below.

[0054]

[0055] Table 1

[0056] The pre-existing plate fault is set as a single pre-existing plate fault in the middle of the model, with its strike parallel to the box boundary and its dip angle set to 60°.

[0057] Figure 2 This is a top view of the initial state of a geological generation model disclosed by the present invention, such as Figure 2 As shown, the deformable base is a 1-meter-long deformable base set on the right side of the pre-existing plate fault. Figure 3 This is a top view of the deformation end of a geological generation model disclosed by the present invention, such as Figure 3 As shown, the base can change shape with the movement of the boundary walls on both sides (i.e., the four walls of the load-bearing box), so that the model can be stretched in any direction without breaking.

[0058] The right side base and boundary of the rigid base experimental model are set as rigid baffles, which can move in any direction at any speed.

[0059] Based on the above geological generation model, this embodiment provides a specific method for three-dimensional numerical simulation of negative inversion structure, including the following steps:

[0060] Step 100, set four groups of experiments for extension in different directions, the early compression direction is positive, and the late extension direction is oblique. The experiment is divided into three periods, namely the first period of compression stage, the second period of cover layer deposition stage and the third period of extension reversal stage.

[0061] Step 200: Apply leftward extrusion stress on the right side of the model at an extrusion speed of 0.01 m / s. Stop extrusion when the extrusion volume reaches 0.75 m and the shrinkage rate reaches 15%.

[0062] Step 300, the experiment begins to enter the second period of cover layer deposition. Based on the end of the first period, a 0.7 m thick layer is deposited again.

[0063] Step 400, the experiment begins to enter the third period of stretch reversal stage. This experimental study sets up an experimental model 11 with early forward extrusion and late forward stretching, and the angle between the late stretching direction and the early extrusion direction is 0°; a rightward stretching stress is applied on the right side, and the stretching speed is set to 0.01m / s until the horizontal stretching length of the deformable substrate reaches 0.75m.

[0064] Step 500, the experiment begins to enter the third period of stretch reversal stage. This experimental study sets up an experimental model 12 with early forward extrusion and late oblique stretching. The angle between the late stretching direction and the early extrusion direction is 30°. A rightward stretching stress is applied on the right side, and the stretching speed is set to 0.01 m / s until the horizontal stretching length of the deformable substrate reaches 0.75 m.

[0065] Step 600, the experiment begins to enter the third period of stretch reversal stage. This experimental study sets up an experimental model 13 with early forward extrusion and late oblique stretching. The angle between the late stretching direction and the early extrusion direction is 45°. A rightward tensile stress is applied on the right side, and the stretching speed is set to 0.01 m / s until the horizontal stretching length of the deformable substrate reaches 0.75 m.

[0066] Step 700, the experiment begins to enter the third period of stretching reversal stage. This experimental study sets up an experimental model 14 with early forward extrusion and late oblique stretching. The angle between the late stretching direction and the early extrusion direction is 60°. A rightward tensile stress is applied on the right side, and the stretching speed is set to 0.01 m / s until the stretching length of the deformable substrate reaches 0.75 m.

[0067] The geometric and kinematic characteristics of the deformation model at each step of the experiment are automatically recorded in real time in the software, including the cumulative displacement of the particles and their velocity vectors. Therefore, it is conducive to the qualitative and quantitative study of the deformation process of the negative inversion structure. Furthermore, the three-dimensional experimental model can combine the plane deformation characteristics and the cross-sectional deformation characteristics, which is conducive to the three-dimensional analysis of such negative inversion structures.

[0068] Based on the characteristic data obtained by the above software, this embodiment provides an analysis method, including the following steps:

[0069] Step 800: Analyze the three-dimensional structural characteristics of the extrusion stage.

[0070] Step 800-1: Since the extrusion stages of the first experiment, the second experiment, the third experiment, and the fourth experiment are the same, the extrusion stages of the four groups of experiments are analyzed together in this section. Figure 4 This is a diagram of the deformation evolution process of the early compression stage of the negative inversion structure disclosed by the present invention, such as Figure 4 As shown, the initial model of the experiment has 5 small layers of different colors to represent different strata, all strata are distributed horizontally, and a pre-existing plate fault 13 with a dip angle of 60° is set in the middle of the model.

[0071] Step 800-2, at the initial stage of compression, the model moves to the left due to the compression force, and the stratum undergoes thrust deformation along the pre-existing plate fault plane. When the compression amount reaches 5% (0.25m), it can be clearly observed from the section that the stratum is displaced along the pre-existing plate fault plane, and a corresponding backthrust fault 17 is developed. On the plane, it can be observed that the stratum between the pre-existing plate fault 13 and the backthrust fault 17 begins to rise.

[0072] Step 800-3, when the extrusion amount reaches 10% (0.5m), the thrust structure of the model becomes more obvious. From the section, it can be observed that the reverse fault (reverse fault is a fault in which the upper plate slides upward along the fault plane) caused by the extrusion effect and the fault distance of the back thrust fault are constantly increasing, and on the plane, it can be observed that the lateral deformation area of ​​the stratum is increasing, and at the same time, the uplift height of the stratum formed is also increasing.

[0073] Step 800-4: When the extrusion amount reaches 15% (0.75m), the fault throws of 13 and 17 reach the maximum, and the stratum uplift height also reaches the maximum. During the extrusion period, the overall deformation of the stratum is mainly thrust structure and stratum thickening.

[0074] Step 800-5: Horizontal displacement analysis: Figure 5 This is a diagram of the total cumulative displacement of particles during the evolution of the early compression stage of the negative inversion structure disclosed by the present invention, such as Figure 5As shown in the figure, the model shows different zoning characteristics in the forward compression stage. According to the different cumulative displacements of the particles, the numbers shown in the figure are also different. According to the numbers in the legend, the larger the displacement, the larger the numerical value of this part, and the value of 0 represents the lowest displacement. Therefore, according to the displacement, the model can be divided into three zones from the overall point of view, namely the static zone on the left, the displacement gradient zone in the middle, and the maximum displacement zone on the right. The left area of ​​the pre-existing plate fault in the model remains static because it is not subjected to stress. Because the particles on the far right of the model are closest to the wall of the compressive stress, and the particles also have enough space to move in the horizontal direction, they show a large displacement. The central uplift area presents a displacement gradient zone, and its highest displacement value mainly comes from the topographic uplift formed by the thrust of the stratum along the pre-existing plate fault plane. The stratum on the front side of the pre-existing plate fault has a cumulative displacement reduction due to the adjustment of the impedance during the thrust process. It is shown in the cumulative displacement cloud map as 1 lower than the surrounding digital value. Compared with the structural evolution profile of the compression stage, it is found that it is also the location where the backthrust fault is developed.

[0075] Step 800-6, vertical displacement analysis: Figure 6 This is a graph of the vertical displacement of particles during the evolution of the early compression stage of the negative inversion structure disclosed by the present invention, such as Figure 6 As shown in the vertical displacement diagram, the central uplift in the middle of the model along the pre-existing plate fault is most significantly affected by the compression. As the compression increases, the vertical displacement increases continuously. The vertical displacement is the largest in the middle of the uplift, which is shown as the area with the highest numerical value, while the vertical displacement gradually decreases toward both sides of the uplift.

[0076] Step 900: 3D structural feature analysis during the extension phase

[0077] Step 900-1: Forward superposition experiment result analysis

[0078] Figure 7 This is a diagram showing the evolution of the first experimental superposition angle 0° stretching stage disclosed by the present invention, as shown in FIG. Figure 7 As shown, the images recorded at four stages during the first experiment when the angle between the stretching direction and the early extrusion direction was 0°. Figure 7 Figure a in the middle is an image of the model after the early compression stage and the cover layer deposition stage. The sedimentary cover layer is the thickest stratum at the top of the model. The cover layer is generally horizontally distributed, covering the stratum where thrust structures developed after the compression stage.

[0079] When the extension reaches 0.25m, the strata below the sedimentary cover still retain the characteristics of compressional structure, but it can be observed that the strata begin to develop fault depression areas along the pre-existing plate faults. For easier observation, the entire fault depression area is encircled by a black solid line in the three-dimensional structural evolution diagram. In the model diagram, it can be seen that the top of the strata in the area of ​​​​the black solid line is depressed due to stretching.

[0080] When the extension reaches 0.5m, the fault zone becomes significantly larger. Although the pre-reverse fault throw can still be observed in the original stratum, its reverse fault throw is constantly decreasing. When the extension reaches 0.5m, the width of the fault zone in the horizontal direction reaches the widest, and the degree of fault in the vertical direction reaches the maximum. Because the model is subjected to unidirectional extension force and there is a pre-existing plate fault, the left side of the fault zone, that is, the side of the pre-existing plate fault, basically remains in the same position, and the fault zone shows the characteristic of continuously extending to the right side of the model.

[0081] Step 900-2: Comparative analysis of superposition test results

[0082] Figure 8 This is a comparison chart of the simulation results of the anisotropic stretching experiment published by the present invention and the first experiment, as shown in Figure 8 As shown, the experimental simulation results of the experiments in which the angles between the extension direction and the early compression direction are 0°, 30°, 45°, and 60° respectively, in the extension deformation process, the final extension reaches 0.75m. For the convenience of observation, a total of 5 horizontal marker layers are set on the top surface every 0.5m. The overall deformation characteristics are similar. The reverse faults formed in the compression stage are reactivated along the pre-existing plate fault plane to form negative inversion structures. The fault depression area develops along the pre-existing plate fault, but the degree of fault depression in the deformation area has the most obvious difference. With the increase of the angle between the extension direction and the early compression direction, the deformation degree and range of the stratum will also increase.

[0083] Step 900-3: Comparative analysis of vertical displacement of superimposed test results

[0084] Fig. 9 This is a comparison diagram of the vertical displacement between the anisotropic stretching experiment and the first experiment published by the present invention, as shown in FIG. Fig. 9As shown in the figure, by comparing the vertical displacement graphs of particles in the four experimental models, positive numbers represent positive displacement in the vertical displacement direction, and the larger the numerical value, the larger the positive displacement; negative numbers represent reverse displacement in the vertical displacement direction, and the larger the numerical value, the larger the reverse displacement. Generally speaking, the characteristics reflected by these groups of experiments are similar. It can be observed that the negative values ​​of the strata above the pre-existing plate fault become smaller, proving that the reverse fault distance is reduced and the negative inversion structure is developed. In the figure, the existence of the boundary fault can be clearly divided according to the different numerical values. By comparing the highest vertical displacement of several groups of experiments, it is found that the vertical change of particles in the first experiment is the highest. Secondly, as the superposition angle increases, the vertical displacement of the model as a whole decreases.

[0085] Step 900-4: Comparative analysis of top surface structural features of superimposed test results

[0086] Fig.10 This is a comparison diagram of the top surface structure of the anisotropic stretching experiment and the first experiment published by the present invention, such as Fig.10 As shown in the figure, the final results of the four experimental models, that is, when the extension reaches 0.75m, are further compared and analyzed. By comparing the three experimental models with different extension directions and the first experiment with positive extension force, it can be found that there are obvious differences between them. In the first experiment, because the strata are subjected to positive extension force, the secondary faults formed are parallel to the main faults, and the width of the graben formed in this experiment is the widest among the four experiments. By analyzing the experimental models of anisotropic extension with different angles to the early compression, it is found that their overall structural characteristics are similar. The left main fault develops multiple faults with a certain angle to it. These small faults are formed due to the force of oblique extension. The oblique force has a strike-slip component, and the faults show the nature of right-lateral strike-slip, thus forming a series of small R-type shear faults. In the experimental model with a 45° stacking angle, the number of small faults on the left is more than that of the 30° stacking angle, because its strike-slip component is larger. Compared with the 60° stacking angle experiment, the strike-slip component is larger. Although the number of small faults formed by strike-slip has not increased, their length has increased, which also proves the enhancement of strike-slip at this angle.

[0087] Generally speaking, when there is an overlapping angle between the late extension direction and the early compression direction, the top surface characteristics of the experimental model simultaneously reflect the deformation characteristics of the coexistence of extension and strike-slip. The smaller the angle, the greater the width between the boundary faults on both sides, and the larger the stratum deformation area.

[0088] Step 900-5: Comparative analysis of the structure characteristics of the superimposed experimental results slices

[0089] In order to more comprehensively study and explore the deformation characteristics of negative inversion structures, based on the above experimental understanding, another perspective is used, that is, slicing from 1 / 2 of the model in the positive vertical model direction to analyze the structural characteristics. Fig.11 This is a comparison diagram of the slice structure of the anisotropic stretching experiment and the first experiment disclosed by the present invention, such as Fig.11 As shown in the figure, the profile results of the four groups of experimental models when the extension reaches 0.75m are compared. The profile deformation characteristics of the four groups of experimental models are relatively similar. It can be observed that the deformation area changes with the different superposition angles. When the superposition angle is 0°, the deformation area of ​​the model is the widest. In order to more accurately compare the differences between the groups of models, the fault throws of the left boundary faults of the four groups of models at this extension stage are counted, and the second small layer from the bottom to the top of the model stratum is used as the marker layer. Fig.11 From ad, the stacking angle is from 0° to 60°, and the vertical fault throw of the pre-existing plate fault is 0.218m, 0.252m, 0.271m, and 0.301m, respectively, showing an increasing trend. Since the compression stage is the same, the reverse fault throw of the strata is the same at the beginning of the four groups of experiments. Therefore, after the extension stage, the larger the strata fault throw is, the weaker the degree of negative reversal is. Therefore, the experimental results show that the larger the angle between the late extension direction and the early compression direction, the greater the degree of negative reversal of the fault. After experiencing the same amount of early compression and late extension, the fault throw of the pre-existing plate fault in each model did not reach 0, that is, the strata did not return to the initial horizontal position. This is because the strata developed secondary faults or part of the extension amount was converted into thinning of the strata during the extension stage, and no obvious fracture was seen.

[0090] Specifically, using the above analysis results, this embodiment interprets the origin of the Jilantai Sag in the Linhe Depression, including the following steps:

[0091] Step 1000-1: Fig.12 The schematic diagram of the three-dimensional structure of the Jilantai Sag in the Linhe Depression disclosed in the present invention is as follows: Fig.12 As shown in Figure 2, a 3D structural model of the Jilantai Sag in the Linhe Depression was established based on seismic data interpretation and field outcrop investigation. Fig.12 The solid arrows represent the compression phase, and the dashed arrows represent the extension phase.

[0092] Step 1000-2: Compare the above simulation results with the structural model of Jilantai Sag. The simulation results show that when the angle between the compression stress and the extension stress direction is 0°, the negative inversion faults developed are parallel in direction and no strike-slip faults are developed, which is similar to the structural model of Jilantai Sag. Therefore, it can be considered that Jilantai Sag is a negative inversion oil and gas basin formed under the regional tectonic stress field conditions with the same direction of compression and extension stress in the northwest-southeast direction, and the stress comes from the Helan Mountain tectonic belt on the east.

[0093] The simulation results of the above experiments are compared with the structural model of Jilantai Sag. It can be seen from the three-dimensional structural model of Jilantai Sag that the negative inversion faults developed in the area are parallel to each other, and there are no strike-slip faults. This model is similar to the stratigraphic model formed when the angle between the compression stress and the extension stress direction is 0° in the first experiment. Therefore, it can be considered that Jilantai Sag is a negative inversion oil and gas basin formed under the conditions of regional tectonic stress field with the same direction of compression and extension stress in the northwest-southeast direction, and the stress comes from the Helan Mountain tectonic belt on the east side of the sag. This conclusion can explain the evolutionary mechanism of the genesis of the negative inversion structure in Jilantai Sag in Linhe Depression, analyze its structural sedimentary features in different periods, and put forward additional suggestions for geological exploration in the area.

[0094] Embodiment 2

[0095] Fig.13 A schematic diagram of a negative inversion structure cause analysis method provided by the present invention is shown in FIG. Fig.13 As shown, the present invention provides a method for analyzing the cause of negative inversion structure, the method comprising:

[0096] S1. Constructing a geological generation model of a negative inversion structure, wherein the geological generation model includes two strata separated by a pre-existing plate fault.

[0097] The pre-existing plate fault refers to the plate fault that existed before the compression stage of the negative inversion structure studied in the present invention, and the fault is an important prerequisite for the formation of the negative inversion structure.

[0098] Geological generation model using PFC 3D The discrete element numerical simulation technology is used to quantitatively study the deformation mechanism of negative inversion structures, and more accurately obtain quantitative data on the main controlling factors of the negative inversion structures superimposed in different directions, which provides a reasonable basis for constructing a three-dimensional model of negative inversion structures, and further provides theoretical support for studying the relationship between negative inversion structures and oil and gas accumulation in Linhe Depression in the Hetao Basin. Furthermore, it can be extended to all negative inversion basins similar to Linhe Depression.

[0099] S2. Controlling the geological generation model in accordance with the compression stage and the extension and reversal stage in sequence, wherein the extension angles of the two strata in the extension and reversal stage are adjustable.

[0100] Referring to the four experiments in Example 1, except for the first experiment in which the direction of stratum movement in the extension reversal stage was coaxial with the direction of movement in the compression stage, the second experiment, the third experiment, and the fourth experiment were all non-coaxial.

[0101] Different extension angles will lead to different final geological structure formations.

[0102] The compression stage and the extension and reversal stage are both the movement of one of the two aforementioned strata relative to the other stratum.

[0103] S3. Record the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage.

[0104] Referring to Example 1, the strata in the model are composed of particles, and the geometric and kinematic characteristics include the cumulative displacement of the particles and their velocity vectors. Therefore, it is conducive to the qualitative and quantitative study of the negative inversion structural deformation process. Furthermore, the three-dimensional experimental model can combine the plane deformation characteristics and the profile deformation characteristics.

[0105] S4. Obtain the geographical features of the area to be analyzed, match the geographical features with the geometric and kinematic features, and determine whether the area to be analyzed belongs to a negative inversion structure. If so, obtain the extension angle of the area to be analyzed in the extension inversion stage.

[0106] With reference to Example 1, after experimental simulations of multiple extension angles, it can be learned that after the extension stage, the greater the stratum fault throw at this time, the weaker the degree of negative reversal. The greater the angle between the late extension direction and the early compression direction, the greater the degree of negative reversal of the fault. For the Jilantai Depression, the development of negative reversal faults in the area is parallel to each other, and no strike-slip faults have been developed. This model is similar to the stratigraphic model formed when the angle between the compression stress and the extension stress direction is 0° in the first experiment. Therefore, it can be considered that the Jilantai Depression is a negative reversal oil and gas basin formed under the conditions of a regional tectonic stress field with the same direction of northwest-southeast compression and extension stress, and the stress comes from the Helan Mountain tectonic belt on the east side of the depression.

[0107] This embodiment can analyze the causes of the negative inversion basin in Linhe Depression, analyze its structural sedimentary features in different periods, and provide supplementary suggestions for geological exploration in the area.

[0108] Further, based on the above-mentioned embodiment, this embodiment provides a method for analyzing the causes of negative inversion structures, and the geological generation model includes:

[0109] A stress-bearing box body, wherein particles are filled inside the stress-bearing box body;

[0110] A first rigid base, a deformable base and a second rigid base connected in sequence and arranged at the bottom of the force-bearing box;

[0111] A pre-existing plate fault disposed at the junction of the first rigid base and the deformable base and forming a first angle with the first rigid base;

[0112] Among them, the first rigid base and the pre-existing plate-like faults and the particles therein constitute a first stratum with a fixed position, and the second rigid base, the deformable base and the pre-existing plate-like faults and the particles therein constitute a second stratum that first squeezes toward the first stratum and then stretches in the opposite direction.

[0113] The deformable base is a 1-meter-long deformable base set on the right side of the pre-existing plate fault. The base can change shape with the movement of the boundary walls on both sides, so that the model can be stretched in any direction without breaking.

[0114] The right side base and boundary of the rigid base experimental model are set as rigid baffles, which can move in any direction at any speed.

[0115] The initial state of the force-bearing box is set to a rectangular parallelepiped with a length of 8m, a width of 3m and a height of 1m. 166,586 spherical particles are filled in the box, and the particle size conforms to the Gaussian distribution.

[0116] In this embodiment, filling particles facilitates subsequent analysis of the displacement of the particles, thereby analyzing the zone characteristics, and the simulation of the stretching movement is achieved by setting a deformable substrate.

[0117] Further, based on the above-mentioned embodiment, this embodiment provides a method for analyzing the cause of negative inversion structure, wherein the property parameters of the filling particles in the load-bearing box are set to porosity 0.15, density 2600.0 kg / m 3 , the particle size of spherical particles is 0.03m.

[0118] In this embodiment, the porosity is 0.15 and the density is 2600.0 kg / m 3 , the particle size of spherical particles is 0.03m to simulate the actual composition of the strata in geological movements. The various characteristics generated by the filling particles under this setting after simulation are most consistent with the actual situation, thereby obtaining detailed geological characteristics. The characteristics that cannot be obtained in theoretical derivation are obtained through actual simulation to guide geological exploration.

[0119] Further, based on the above-mentioned embodiment, this embodiment provides a method for analyzing the causes of negative inversion structures, wherein the angle between the pre-existing plate fault and the first rigid base is 60°.

[0120] The pre-existing plate fault is provided with a single pre-existing plate fault in the middle of the model, and its strike is parallel to the boundary of the box.

[0121] This embodiment provides a typical angle between a pre-existing plate fault and a first rigid substrate, which is obtained by means of data analysis and expert experience of the inventor.

[0122] Further, based on the above-mentioned embodiment, this embodiment provides a method for analyzing the genesis of negative inversion structures, which records the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage, including:

[0123] The vertical displacement analysis and horizontal displacement analysis are performed on the constituent particles of the geological generation model to form multiple zones.

[0124] During the extrusion stage, the model shows different zoning characteristics, and multiple zones are divided according to the different cumulative displacements of the particles.

[0125] In the extension stage, more zones are further divided through vertical displacement analysis and horizontal displacement analysis.

[0126] This embodiment performs vertical displacement analysis and horizontal displacement analysis on the constituent particles of the geological generation model, thereby facilitating rapid and accurate analysis of various zones formed in the anti-structure.

[0127] Furthermore, based on the above-mentioned embodiment, this embodiment provides a method for analyzing the causes of negative inversion structures, wherein a cap layer deposition stage is provided between the compression stage and the extension inversion stage.

[0128] During the cap layer deposition stage, based on the end of the compression stage, a 0.7m thick stratum was deposited again.

[0129] The cap layer deposition in this embodiment makes the setting of the entire model closer to the actual formation process of geological structures.

[0130] A negative inversion structure cause analysis system provided by the present invention is described below. The negative inversion structure cause analysis system described below and the negative inversion structure cause analysis method described above can be referenced to each other.

[0131] Fig.14 A schematic diagram of a negative inversion structure cause analysis system provided by the present invention is shown in FIG. Fig.14 As shown, the present invention also provides a negative inversion structure cause analysis system, the system comprising:

[0132] A tectonic module, the tectonic module constructs a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault;

[0133] A generation module, the generation module controls the geological generation model in accordance with the compression stage and the extension and reversal stage, wherein the extension angles of the two strata in the extension and reversal stage are adjustable;

[0134] A recording module, which records the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage;

[0135] The judgment module obtains the geographical features of the area to be analyzed, matches them with the geometric and kinematic features, and determines whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0136] This embodiment can analyze the causes of the negative inversion basin in Linhe Depression, analyze its structural sedimentary features in different periods, and provide supplementary suggestions for geological exploration in the area.

[0137] Fig.15 A schematic diagram of the physical structure of an electronic device provided by the present invention, such as Fig.15 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the negative inversion structure cause analysis method, and the method includes:

[0138] Constructing a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault;

[0139] The geological generation model is controlled in sequence according to the compression stage and the extension and reversal stage, wherein the extension angle of the two strata in the extension and reversal stage is adjustable;

[0140] Record the geometric and kinematic characteristics of the geological generation model at different extension angles at each stage;

[0141] The geographical features of the area to be analyzed are obtained, and the geographical features are matched with the geometric and kinematic features to determine whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0142] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0143] On the other hand, the present invention further provides a computer program product, the computer program product comprising a computer program stored on a non-transitory computer-readable storage medium, the computer program comprising program instructions, when the program instructions are executed by a computer, the computer can execute the negative inversion structure cause analysis method provided by the above methods, the method comprising:

[0144] Constructing a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault;

[0145] The geological generation model is controlled in sequence according to the compression stage and the extension and reversal stage, wherein the extension angle of the two strata in the extension and reversal stage is adjustable;

[0146] Record the geometric and kinematic characteristics of the geological generation model at different extension angles at each stage;

[0147] The geographical features of the area to be analyzed are obtained, and the geographical features are matched with the geometric and kinematic features to determine whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0148] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to perform the above-mentioned negative inversion structure cause analysis method, the method comprising:

[0149] Constructing a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault;

[0150] The geological generation model is controlled in sequence according to the compression stage and the extension and reversal stage, wherein the extension angle of the two strata in the extension and reversal stage is adjustable;

[0151] Record the geometric and kinematic characteristics of the geological generation model at different extension angles at each stage;

[0152] The geographical features of the area to be analyzed are obtained, and the geographical features are matched with the geometric and kinematic features to determine whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

[0153] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0154] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0155] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for analyzing the causes of negative inversion structures, characterized in that: The method comprises: Constructing a geological generative model of a negative inversion structure, wherein the geological generative model includes two strata separated by a pre-existing plate fault; The geological generation model is controlled in sequence according to the compression stage and the extension and reversal stage, wherein the extension angles of the two strata in the extension and reversal stage are adjustable; Recording the geometric and kinematic characteristics of the geological generation model at different extension angles at each stage; Acquire the geographical features of the area to be analyzed, match the geographical features with the geometric and kinematic features, determine whether the area to be analyzed belongs to a negative inversion structure, and if so, obtain the extension angle of the area to be analyzed in the extension inversion stage; The geological generation model includes: A stress-bearing box body, wherein particles are filled inside the stress-bearing box body; A first rigid base, a deformable base and a second rigid base which are connected in sequence and are arranged at the bottom of the force-bearing box; A pre-existing plate fault disposed at the junction of the first rigid base and the deformable base and forming a first angle with the first rigid base; Among them, the first rigid base and the pre-existing plate-like faults and the particles therein constitute a first stratum with a fixed position, and the second rigid base, the deformable base and the pre-existing plate-like faults and the particles therein constitute a second stratum that first squeezes toward the first stratum and then stretches in the opposite direction.

2. The method for analyzing the cause of negative inversion structure according to claim 1, characterized in that: The property parameters of the filling particles in the load-bearing box are set to porosity 0.15, density 2600.0 kg / m 3 , the particle size of spherical particles is 0.03m.

3. The method for analyzing the cause of negative inversion structure according to claim 1, characterized in that: The angle between the pre-existing plate fault and the first rigid substrate is 60°.

4. The method for analyzing the cause of negative inversion structure according to claim 1, characterized in that: The geometric and kinematic characteristics of the geological generation model at different extension angles in each stage are recorded, including: The vertical displacement analysis and horizontal displacement analysis are performed on the constituent particles of the geological generation model to form multiple zones.

5. The method for analyzing the cause of negative inversion structure according to claim 1, characterized in that: A cap layer deposition stage is also provided between the compression stage and the stretching and reversing stage.

6. A negative inversion structure cause analysis system, characterized in that: The system comprises: A construction module, wherein the construction module constructs a geological generation model of a negative inversion structure, wherein the geological generation model includes two strata separated by a pre-existing plate fault; the geological generation model includes a force-bearing box, wherein the interior of the force-bearing box is filled with particles; a first rigid base, a deformable base, and a second rigid base are sequentially connected and arranged at the bottom of the force-bearing box; a pre-existing plate fault is arranged at the connection between the first rigid base and the deformable base and is at a first angle to the first rigid base; wherein the first rigid base and the pre-existing plate fault and the particles therein constitute a first stratum with a fixed position, and the second rigid base, the deformable base, the pre-existing plate fault and the particles therein constitute a second stratum that is first squeezed toward the first stratum and then extends in the opposite direction; A generation module, wherein the generation module controls the geological generation model in accordance with the compression stage and the extension and reversal stage in sequence, wherein the extension angles of the two strata in the extension and reversal stage are adjustable; A recording module, wherein the recording module records the geometric and kinematic characteristics of the geological generation model at different extension angles in each stage; A judgment module, wherein the judgment module obtains the geographical features of the area to be analyzed, matches them with the geometric and kinematic features, and determines whether the area to be analyzed belongs to a negative inversion structure. If so, the extension angle of the area to be analyzed in the extension inversion stage is obtained.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the negative inversion structure cause analysis method according to any one of claims 1 to 5 are implemented.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for analyzing the cause of negative inversion structure according to any one of claims 1 to 5 are implemented.

9. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method for analyzing the cause of negative inversion structure according to any one of claims 1 to 5 are implemented.

Images