Geological structure delineation method and system based on cross gradient constraint

Abstract

The application discloses a kind of geological structure delineation method and system based on cross gradient constraint, including according to geological rock distribution characteristics unit time acquisition geological state information, according to the geological state information judges crustal movement period;Based on cross gradient model, the distribution of each rock layer is delineated, based on distribution situation receives the delineation request of volcanic eruption input, generates delineation task, and sends to host computer;According to plate fracture information, obtain plate displacement trajectory, calculate deviation according to plate displacement trajectory, and generate geological structure mutation information according to deviation.The application establishes delineation model, directly reflects the state of each rock layer, adjusts distribution situation per unit time, provides guidance for delineation;When plate displacement, no obstacle is set, but the behavior of plate displacement is monitored per unit time, and is displayed in delineation model, the behavior of plate displacement is analyzed, and geological structure mutation information is generated, which greatly improves the accuracy of geological structure delineation.

Description

Technical Field

[0001] This invention relates to the field of geological structure characterization, and in particular to a geological structure characterization method and system based on cross gradient constraints. Background Technology

[0002] Traditional single non-seismic inversion methods often produce multiple solutions, which are insufficient for characterizing deep geological bodies and cannot effectively solve problems such as the physical properties of deep strata and their residual thickness distribution.

[0003] Geological structure characterization is a current hot topic and trend in geophysical inversion research. Unlike single geophysical inversion, geological structure characterization requires establishing the correlation between parameters of different physical models, that is, the joint or coupling of different physical property parameters. This is the key to the realization of geological structure characterization. The more representative coupling methods can be divided into two categories: one is based on the empirical relationship between rock physical properties, and the other is to constrain the consistency or similarity of structural or geometric feature changes of different physical property parameters. The main problems in geological characterization research are: firstly, the reliance on empirical relationships of rock properties on prior information is strong, while the physical property relationships of actual problems may have statistical differences or insufficient representativeness and are complex, making them difficult to apply universally in the study of complex problems; secondly, the curvature reading method and Gramian determinant constraint method proposed by predecessors are difficult to construct intuitive or explicit methods, and are difficult to implement in practice; thirdly, the current non-seismic geological characterization methods do not make sufficient use of constraint information such as seismic and drilling data, and existing methods are difficult to integrate constraint information into the process, which greatly restricts the ability of integrating various geophysical methods and reliable constraint information to solve complex geological problems and enhance the uniqueness of solutions. Summary of the Invention

[0004] The purpose of this invention is to provide a geological structure characterization method and system based on cross gradient constraints to solve the problems mentioned in the background art.

[0005] To achieve the above objectives, the present invention provides the following technical solution:

[0006] A geological structure characterization method based on cross-gradient constraints, the method comprising:

[0007] Geological state information is obtained per unit time based on the distribution characteristics of geological rocks, and the crustal movement cycle is determined based on the geological state information; the crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by the Earth's timeline, and the Earth's timeline corresponding to volcanic eruptions contains the geological state change time and the time node containing the eruption; the Earth's timeline corresponding to plate displacements contains the geological state change time and the geological state sequence number within a time period.

[0008] When the crustal movement cycle is a volcanic eruption, the control variables of the preset cross gradient model are adjusted according to the geological state information;

[0009] The distribution of each rock layer is characterized based on the cross gradient model. Based on the distribution, the characterization request input from the volcanic eruption is received, a characterization task is generated, and sent to the host.

[0010] When the crustal movement cycle is plate displacement, the plate displacement trajectory is obtained based on plate faulting information, the deviation value is calculated based on the plate displacement trajectory, and geological structural abrupt change information is generated based on the deviation value.

[0011] As a further aspect of the present invention: the step of adjusting the control variables of the preset cross-gradient model according to geological state information when the crustal movement cycle is a volcanic eruption includes:

[0012] When the crustal movement cycle is a volcanic eruption, obtain the time node containing the eruption.

[0013] Adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle;

[0014] Obtain the time of geological state changes during the volcanic eruption, and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time;

[0015] The control variables for the corresponding parameters are determined based on the formation cycle of the geological structure; the control variables include both Earth's own factors and extraterrestrial factors.

[0016] When the geological rock distribution characteristics detect volcanic rock information from the volcanic eruption, the Earth's own factors are adjusted.

[0017] As a further aspect of the present invention: when the crustal movement cycle is a volcanic eruption, the step of obtaining the time node containing the eruption of the volcano includes:

[0018] When the crustal movement cycle is a volcanic eruption, obtain the stratigraphic evolution history of that volcanic eruption;

[0019] The magnitude of the volcanic eruption was determined based on the geological evolution history described above.

[0020] The stratigraphic evolution time is expanded according to the level to generate a stratigraphic evolution time containing virtual time.

[0021] As a further aspect of the present invention: the steps of characterizing the distribution of each rock layer based on the cross-gradient model, receiving the characterization request input from the volcanic eruption based on the distribution, generating the characterization task, and sending it to the host include:

[0022] Traverse the control variables of each rock layer in the cross gradient model;

[0023] When the control variable is an accumulated comprehensive state, it describes the rock distribution.

[0024] As a further aspect of the present invention: when the crustal movement period is plate displacement, the step of obtaining the plate displacement trajectory based on plate fracture information includes:

[0025] When the crustal movement cycle is plate displacement, image features of rock strata are acquired; the image features are image sub-blocks containing rock strata.

[0026] Based on the plate displacement position, the fracture information of the plate to be inspected is selected, and based on the image features, the fracture information of the target plate is selected from the fracture information of the plate to be inspected.

[0027] Determine the geological structure-time correlation based on the location of the fault in the target plate;

[0028] Update the plate displacement position according to the aforementioned correlation, and repeat the process cyclically;

[0029] By statistically analyzing and connecting all plate displacement positions, plate displacement trajectories are generated in the cross gradient model.

[0030] As a further aspect of the present invention: the step of selecting target plate fracture information from the fracture information of the plate to be inspected based on the image features includes:

[0031] The fracture information of the plate to be inspected is obtained sequentially. The fracture information of the plate to be inspected and the image sub-block are input into the characterization model. The fracture information of the plate to be inspected is selected as the fracture information of the target plate based on the characterization result.

[0032] The characterization model is:

[0033]

[0034] in, η represents the degree of plate fracturing in the previous unit of time, ω represents the covariance of the plate fracturing coefficient at different time points, K represents the estimated time of plate formation, and ω represents the number of years the plate has been formed. This represents the coefficient matrix of all other factors affecting plate tectonics. This represents the fracture coefficient of a plate per unit time. This represents the fracture coefficient of the plate in the previous unit of time. This represents all other noise factors that affect the plate fracture coefficient. This indicates the length of time it takes for a plate to form within a given unit of time. This indicates the length of time in which the plate was formed up to the previous unit of time. This represents all other noise factors that affect the formation of tectonic plates over a long period of time.

[0035] As a further aspect of the present invention: the step of calculating the deviation value based on the plate displacement trajectory and generating geological structural abrupt change information based on the deviation value includes:

[0036] The end time of plate displacement is detected based on the geological rock distribution characteristics, and the displacement time is calculated based on the end time and the time of geological state change.

[0037] The endpoint of the plate displacement is determined based on the plate displacement trajectory. The endpoint of the plate displacement and the displacement time are input into the preset judgment model to calculate the deviation value.

[0038] Geological structural mutation information is generated based on the deviation value; the geological structural mutation information contains target plate fracture information corresponding to the displacement position of each plate.

[0039] The present invention also provides a geological structure characterization system based on cross-gradient constraints. The system includes: a crustal movement cycle determination module, used to acquire geological state information per unit time based on the distribution characteristics of geological rocks, and to determine the crustal movement cycle based on the geological state information; the crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by an Earth chronology, where the Earth chronology corresponding to volcanic eruptions includes the geological state change time and the time node containing the eruption; the Earth chronology corresponding to plate displacements includes the geological state change time and the geological state sequence number within a time cycle.

[0040] The control variable adjustment module is used to adjust the control variables of the preset cross gradient model according to geological state information when the crustal movement cycle is a volcanic eruption.

[0041] The rock strata characterization module is used to characterize the distribution of each rock stratum based on the cross gradient model. It receives characterization requests from volcanic eruptions based on the distribution, generates characterization tasks, and sends them to the host.

[0042] The plate displacement monitoring module is used to obtain the plate displacement trajectory based on plate fault information when the crustal movement cycle is plate displacement, calculate the deviation value based on the plate displacement trajectory, and generate geological structure change information based on the deviation value.

[0043] As a further aspect of the present invention: the control variable adjustment module includes:

[0044] The crustal movement cycle acquisition unit is used to acquire the time node containing the eruption when the crustal movement cycle is a volcanic eruption.

[0045] The parameter adjustment unit is used to adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle;

[0046] The geological structure formation cycle calculation unit is used to obtain the geological state change time of the volcanic eruption and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time.

[0047] A control variable determination unit is used to determine control variables for corresponding parameters based on the formation cycle of the geological structure; the control variables include Earth's own factors and extraterrestrial factors.

[0048] The factor control unit is used to adjust the Earth's own factors when the geological rock distribution characteristics detect volcanic rock information of the volcanic eruption.

[0049] As a further aspect of the present invention: the plate displacement monitoring module includes:

[0050] The image feature acquisition unit is used to acquire image features of rock strata when the crustal movement cycle is plate displacement; the image features are image sub-blocks containing rock strata.

[0051] The characterization selection unit is used to select fracture information of the plate to be inspected based on the plate displacement position, and to select fracture information of the target plate from the fracture information of the plate to be inspected based on the image features.

[0052] The geological structure-time acquisition unit is used to determine the geological structure-time correlation based on the location of the fault in the target plate.

[0053] The plate position dynamic unit is used to update the plate displacement position according to the aforementioned correlation relationship and executes it cyclically;

[0054] The plate displacement generation unit is used to statistically analyze and connect all plate displacement positions, and generate plate displacement trajectories in the cross gradient model.

[0055] Compared with the prior art, the beneficial effects of the present invention are: the present invention establishes a characterization model that intuitively reflects the state of each rock layer and adjusts the distribution per unit time, providing guidance for characterization; when the plate moves, no obstacles are set up, but the behavior of the plate movement is monitored per unit time and displayed in the characterization model, the behavior of the plate movement is analyzed, and geological structure change information is generated to cooperate with the management for safety management; while ensuring safety as much as possible, the accuracy of geological structure characterization is greatly improved. Attached Figure Description

[0056] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention.

[0057] Figure 1This is a flowchart of a geological structure characterization method based on cross gradient constraints.

[0058] Figure 2 This is the first sub-flowchart of a geological structure characterization method based on cross-gradient constraints.

[0059] Figure 3 This is the second sub-flowchart of a geological structure characterization method based on cross-gradient constraints.

[0060] Figure 4 This is a block diagram of the composition of a geological structure characterization system based on cross-gradient constraints. Detailed Implementation

[0061] To make the technical problems to be solved, the technical solutions, and the beneficial effects of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present invention and are not intended to limit the present invention.

[0062] Example 1

[0063] Figure 1 This is a flowchart of a geological structure characterization method based on cross gradient constraints. In this embodiment of the invention, a geological structure characterization method based on cross gradient constraints includes:

[0064] Step A1: Obtain geological state information per unit time based on the distribution characteristics of geological rocks, and determine the crustal movement cycle based on the geological state information; the crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by the Earth's timeline, and the Earth's timeline corresponding to volcanic eruptions contains the geological state change time and the time node containing the eruption; the Earth's timeline corresponding to plate displacements contains the geological state change time and the geological state sequence number within a time cycle;

[0065] Geological rock distribution characteristics can be identified by displacement identification plates installed at the system entrance. When the system identifies the area of ​​the displacement identification plate, it is considered that the system has sent a geological status request. At this time, displacement identification is performed on the system to obtain the crustal movement cycle information of the system. Based on the crustal movement cycle information, the system is adjusted to see if there are any changes. If there are changes, it is a volcanic eruption; if there are no changes, it is plate displacement.

[0066] Step A2: When the crustal movement cycle is a volcanic eruption, adjust the control variables of the preset cross gradient model according to the geological state information;

[0067] Step A3: Based on the cross-gradient model, characterize the distribution of each rock layer, receive the characterization request input from the volcanic eruption based on the distribution, generate the characterization task, and send it to the host.

[0068] Step A4: When the crustal movement cycle is plate displacement, obtain the plate displacement trajectory based on the plate fracture information, calculate the deviation value based on the plate displacement trajectory, and generate geological structure abrupt change information based on the deviation value;

[0069] Many existing control and management systems do not allow input from sources other than volcanic eruptions. Access to the control and management system is only possible after registration in the background. While this method is extremely secure, it can be inconvenient.

[0070] The above content obtains the movement trajectory of plate displacement based on the location of occurrence installed in the system, judges whether there is any anomaly based on the movement trajectory per unit time, and generates geological structure change information based on the presence of anomalies. At this time, the safety assurance process of the control system becomes more concealed, and greatly improves the accuracy of stratigraphic evolution and geological structure characterization while ensuring timely detection of problems.

[0071] Figure 2 The first sub-flowchart of a geological structure characterization method based on cross-gradient constraints includes the step of adjusting the control variables of a preset cross-gradient model according to geological state information when the crustal movement cycle is a volcanic eruption.

[0072] Step B1: When the crustal movement cycle is a volcanic eruption, obtain the time node containing the eruption.

[0073] Step B2: Adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle;

[0074] When the crustal movement cycle is a volcanic eruption, the stratigraphic evolution information of that volcanic eruption is adjusted. This stratigraphic evolution information is the crustal movement cycle containing time. The corresponding parameters in the cross gradient model can be adjusted based on the crustal movement cycle.

[0075] Step B3: Obtain the geological state change time of the volcanic eruption, and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time;

[0076] Step B4: Determine the control variables for the corresponding parameters based on the formation cycle of the geological structure; the control variables include Earth's own factors and extraterrestrial factors;

[0077] Step B5: When the geological rock distribution characteristics detect volcanic rock information from the volcanic eruption, adjust the Earth's own factors.

[0078] The geological rock distribution characteristics detect both the entry process and the exit process of the rock strata. If the rock strata exit, the control variables of the corresponding parameters in the cross gradient model are adjusted. The adjustment can be understood as resetting to the default value.

[0079] When rock strata are detected to have entered, the geological structure formation cycle is calculated based on the geological state change time and the strata evolution time. The control variable is determined based on the geological structure formation cycle to represent the remaining usage time of the parameter.

[0080] As a preferred embodiment of the technical solution of the present invention, the step of obtaining the time node containing the eruption of the volcano when the crustal movement cycle is a volcanic eruption includes:

[0081] When the crustal movement cycle is a volcanic eruption, obtain the stratigraphic evolution history of that volcanic eruption;

[0082] The magnitude of the volcanic eruption was determined based on the geological evolution history described above.

[0083] The stratigraphic evolution time is expanded according to the level to generate a stratigraphic evolution time containing virtual time.

[0084] Under the framework of the technical solution of this invention, the level of a volcanic eruption is adjusted according to the stratigraphic evolution history of the eruption, thereby determining a virtual time. The meaning of the virtual time needs to be explained with examples. For instance, if the stratigraphic evolution time of a volcanic eruption is ten years, and the eruption is the first geological state of the eruption, then its virtual time is also set to ten years. The significance is that when the geological state changes again within ten years after the end of the stratigraphic evolution, the volcanic eruption should have been a plate displacement. However, the computer plate view is a volcanic eruption. At this time, if there are other parameters, the control variable of this parameter will be in an accumulated and comprehensive state (accumulated and comprehensive control variables); if there are no other parameters, the volcanic eruption can still open the corresponding parameters.

[0085] Based on this, if a geological condition remains unchanged for a long period of time, then the virtual time will be even longer, such as thirty years, a week, or even longer.

[0086] Based on the above, the further steps of characterizing the distribution of each rock layer based on the cross-gradient model, receiving the characterization request from the volcanic eruption input based on the distribution, generating the characterization task, and sending it to the host include:

[0087] The control variables of each rock layer in the cross gradient model are traversed; when the control variables are in an accumulated comprehensive state, the rock distribution is characterized.

[0088] When a parameter is in an accumulating and synthesizing state, the characterization may include improving the parameter input and output records, etc., which is determined by the developers as needed, and the technical solution of this invention is not limited thereto.

[0089] Figure 3 The second sub-flowchart of a geological structure characterization method based on cross-gradient constraints, wherein the step of obtaining the plate displacement trajectory based on plate fault information when the crustal movement period is plate displacement, includes:

[0090] Step C1: When the crustal movement period is plate displacement, acquire the image features of the rock strata; the image features are image sub-blocks containing the rock strata;

[0091] Plate tectonics is a relatively sensitive crustal movement cycle. Therefore, when the crustal movement cycle is plate displacement, geological structural information is obtained based on the location of the occurrence of the displacement based on the distribution characteristics of geological rocks, and the position of the rock strata is located within this geological structural information. This process is not difficult and can be referenced from existing target tracking technologies. For example, when performing displacement recognition, a small box will encompass the displacement, which is a target tracking technology. The technical solution of this invention tracks plate displacement, and the corresponding small box is an image sub-block.

[0092] Step C2: Select the fracture information of the plate to be inspected based on the plate displacement position, and select the fracture information of the target plate from the fracture information of the plate to be inspected based on the image features;

[0093] Based on the plate displacement position, the location where the plate displacement occurs can be captured. The corresponding plate fracture information is the plate fracture information to be inspected. According to the above image features, plate fracture information containing plate displacement can be delineated from many plate fracture information to be inspected, which is the target plate fracture information.

[0094] In this invention, the location of occurrence, by default, covers all public areas of the system.

[0095] Step C3: Determine the geological structure-time correlation based on the location of the fault in the target plate;

[0096] There is a correlation between the position of each pixel in the location and the time. This correlation is determined after the location is installed and is considered the default value.

[0097] Step C4: Update the plate displacement position according to the aforementioned correlation, and repeat the process cyclically;

[0098] Step C5: Statistically analyze and connect all plate displacement positions to generate plate displacement trajectories in the cross gradient model;

[0099] Based on the correlation, the plate displacement positions are updated sequentially, and the plate displacement positions are connected to generate the plate displacement trajectory, which is then displayed in the cross gradient model.

[0100] As a preferred embodiment of the technical solution of the present invention, the step of selecting the target plate fracture information based on the image features from the fracture information of the plate to be inspected includes:

[0101] The fracture information of the plate to be inspected is obtained sequentially. The fracture information of the plate to be inspected and the image sub-block are input into the characterization model. The fracture information of the plate to be inspected is selected as the fracture information of the target plate based on the characterization result.

[0102] The characterization model is:

[0103]

[0104] in, η represents the degree of plate fracturing in the previous unit of time, ω represents the covariance of the plate fracturing coefficient at different time points, K represents the estimated time of plate formation, and ω represents the number of years the plate has been formed. This represents the coefficient matrix of all other factors affecting plate tectonics. This represents the fracture coefficient of a plate per unit time. This represents the fracture coefficient of the plate in the previous unit of time. This represents all other noise factors that affect the plate fracture coefficient. This indicates the length of time it takes for a plate to form within a given unit of time. This indicates the length of time in which the plate was formed up to the previous unit of time. This represents all other noise factors that affect the formation of tectonic plates over a long period of time.

[0105] As a preferred embodiment of the technical solution of the present invention, the step of calculating the deviation value based on the plate displacement trajectory and generating geological structural change information based on the deviation value includes:

[0106] The end time of plate displacement is detected based on the geological rock distribution characteristics, and the displacement time is calculated based on the end time and the time of geological state change.

[0107] The endpoint of the plate displacement is determined based on the plate displacement trajectory. The endpoint of the plate displacement and the displacement time are input into the preset judgment model to calculate the deviation value.

[0108] Geological structural mutation information is generated based on the deviation value; the geological structural mutation information contains target plate fracture information corresponding to the displacement position of each plate.

[0109] Whether a plate displacement is dangerous is determined by the plate displacement trajectory and the plate displacement time. The plate displacement trajectory is generated based on the above-mentioned information, while the plate displacement time is generated based on the identification process of geological rock distribution characteristics. The plate displacement trajectory can determine at which parameter or location the plate is displacing, and combined with the displacement time, it can be determined whether the plate displacement is abnormal.

[0110] Example 2

[0111] Figure 4 This is a structural block diagram of a geological structure characterization system based on cross-gradient constraints. In this embodiment of the invention, a geological structure characterization system based on cross-gradient constraints, system 10 includes:

[0112] The crustal movement cycle determination module 11 is used to obtain geological state information per unit time based on the geological rock distribution characteristics, and to determine the crustal movement cycle based on the geological state information; the crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by the Earth's timeline, and the Earth's timeline corresponding to volcanic eruptions contains the geological state change time and the time node containing the eruption; the Earth's timeline corresponding to plate displacements contains the geological state change time and the geological state sequence number within a time cycle.

[0113] The control variable adjustment module 12 is used to adjust the control variables of the preset cross gradient model according to the geological state information when the crustal movement cycle is a volcanic eruption.

[0114] The rock strata characterization module 13 is used to characterize the distribution of each rock stratum based on the cross gradient model. It receives characterization requests from volcanic eruptions based on the distribution, generates characterization tasks, and sends them to the host.

[0115] The plate displacement monitoring module 14 is used to obtain the plate displacement trajectory based on plate fault information when the crustal movement cycle is plate displacement, calculate the deviation value based on the plate displacement trajectory, and generate geological structure mutation information based on the deviation value.

[0116] The control variable adjustment module 12 includes:

[0117] The crustal movement cycle acquisition unit is used to acquire the time node containing the eruption when the crustal movement cycle is a volcanic eruption.

[0118] The parameter adjustment unit is used to adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle;

[0119] The geological structure formation cycle calculation unit is used to obtain the geological state change time of the volcanic eruption and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time.

[0120] A control variable determination unit is used to determine control variables for corresponding parameters based on the formation cycle of the geological structure; the control variables include Earth's own factors and extraterrestrial factors.

[0121] The factor control unit is used to adjust the Earth's own factors when the geological rock distribution characteristics detect volcanic rock information of the volcanic eruption.

[0122] The plate displacement monitoring module 14 includes:

[0123] The image feature acquisition unit is used to acquire image features of rock strata when the crustal movement cycle is plate displacement; the image features are image sub-blocks containing rock strata.

[0124] The characterization selection unit is used to select fracture information of the plate to be inspected based on the plate displacement position, and to select fracture information of the target plate from the fracture information of the plate to be inspected based on the image features.

[0125] The geological structure-time acquisition unit is used to determine the geological structure-time correlation based on the location of the fault in the target plate.

[0126] The plate position dynamic unit is used to update the plate displacement position according to the aforementioned correlation relationship and executes it cyclically;

[0127] The plate displacement generation unit is used to statistically analyze and connect all plate displacement positions, and generate plate displacement trajectories in the cross gradient model.

[0128] Although embodiments of the invention have been shown and described, it will be understood by those skilled in the art that various equivalent changes, modifications, substitutions and alterations can be made to these embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the appended claims and their equivalents.

Claims

1. A method for characterizing geological structures based on cross-gradient constraints, characterized in that, The method includes: Geological state information is obtained per unit time based on the distribution characteristics of geological rocks, and the crustal movement cycle is determined based on the geological state information; the crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by the Earth's timeline, and the Earth's timeline corresponding to volcanic eruptions contains the geological state change time and the time node containing the eruption; the Earth's timeline corresponding to plate displacements contains the geological state change time and the geological state sequence number within a time period. When the crustal movement cycle is a volcanic eruption, the control variables of the preset cross gradient model are adjusted according to the geological state information; The distribution of each rock layer is characterized based on the cross gradient model. Based on the distribution, the characterization request input from the volcanic eruption is received, a characterization task is generated, and sent to the host. When the crustal movement cycle is plate displacement, the plate displacement trajectory is obtained based on plate fracture information, the deviation value is calculated based on the plate displacement trajectory, and geological structural abrupt change information is generated based on the deviation value. When the crustal movement period is plate displacement, the step of obtaining the plate displacement trajectory based on plate faulting information includes: When the crustal movement cycle is plate displacement, image features of rock strata are acquired; the image features are image sub-blocks containing rock strata. Based on the plate displacement position, the fracture information of the plate to be inspected is selected, and based on the image features, the fracture information of the target plate is selected from the fracture information of the plate to be inspected. Determine the geological structure-time correlation based on the location of the fault in the target plate; Update the plate displacement position according to the aforementioned correlation, and repeat the process cyclically; Statistically analyze and connect all plate displacement positions to generate plate displacement trajectories in the cross gradient model; The step of selecting target plate fracture information from the fracture information of the plate to be inspected based on the image features includes: The fracture information of the plate to be inspected is obtained sequentially. The fracture information of the plate to be inspected and the image sub-block are input into the characterization model. The fracture information of the plate to be inspected is selected as the fracture information of the target plate based on the characterization result. The characterization model is: in, η represents the degree of plate fracturing in the previous unit of time, ω represents the covariance of the plate fracturing coefficient at different time points, K represents the estimated time of plate formation, and ω represents the number of years the plate has been formed. This represents the coefficient matrix of all other factors affecting plate tectonics. This represents the fracture coefficient of a plate per unit time. This represents the fracture coefficient of the plate in the previous unit of time. This represents all other noise factors that affect the plate fracture coefficient. This indicates the length of time it takes for a plate to form within a given unit of time. This indicates the length of time in which the plate was formed up to the previous unit of time. This represents all other noise factors that affect the formation of tectonic plates over a long period of time.

2. The geological structure characterization method based on cross-gradient constraints according to claim 1, characterized in that, The step of adjusting the control variables of the preset cross-gradient model based on geological state information when the crustal movement cycle is a volcanic eruption includes: When the crustal movement cycle is a volcanic eruption, obtain the time node containing the eruption. Adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle; Obtain the time of geological state changes during the volcanic eruption, and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time; The control variables for the corresponding parameters are determined based on the formation cycle of the geological structure; the control variables include both Earth's own factors and extraterrestrial factors. When the geological rock distribution characteristics detect volcanic rock information from the volcanic eruption, the Earth's own factors are adjusted.

3. The geological structure characterization method based on cross-gradient constraints according to claim 2, characterized in that, When the crustal movement cycle is a volcanic eruption, the step of obtaining the time node containing the eruption of the volcano includes: When the crustal movement cycle is a volcanic eruption, obtain the stratigraphic evolution history of that volcanic eruption; The magnitude of the volcanic eruption was determined based on the geological evolution history described above. The stratigraphic evolution time is expanded according to the level to generate a stratigraphic evolution time containing virtual time.

4. The geological structure characterization method based on cross-gradient constraints according to claim 3, characterized in that, The steps of characterizing the distribution of each rock layer based on the cross-gradient model, receiving characterization requests from volcanic eruptions based on the distribution, generating characterization tasks, and sending them to the host include: Traverse the control variables of each rock layer in the cross gradient model; When the control variable is an accumulated comprehensive state, it describes the rock distribution.

5. The geological structure characterization method based on cross-gradient constraints according to claim 1, characterized in that, The steps of calculating the deviation value based on the plate displacement trajectory and generating geological structural abrupt change information based on the deviation value include: The end time of plate displacement is detected based on the geological rock distribution characteristics, and the displacement time is calculated based on the end time and the time of geological state change. The endpoint of the plate displacement is determined based on the plate displacement trajectory. The endpoint of the plate displacement and the displacement time are input into the preset judgment model to calculate the deviation value. Geological structural mutation information is generated based on the deviation value; the geological structural mutation information contains target plate fracture information corresponding to the displacement position of each plate.

6. A geological structure characterization system based on cross-gradient constraints, the system implementing the method as described in claim 1, characterized in that, The system includes: The crustal movement cycle determination module is used to obtain geological state information per unit time based on the geological rock distribution characteristics, and to determine the crustal movement cycle based on the geological state information. The crustal movement cycle includes volcanic eruptions and plate displacements; the crustal movement cycle is represented by the Earth's timeline, the Earth's timeline corresponding to volcanic eruptions includes the geological state change time and the time node containing the eruption; the Earth's timeline corresponding to plate displacements includes the geological state change time and the geological state sequence number within a time cycle. The control variable adjustment module is used to adjust the control variables of the preset cross gradient model according to geological state information when the crustal movement cycle is a volcanic eruption. The rock strata characterization module is used to characterize the distribution of each rock stratum based on the cross gradient model. It receives characterization requests from volcanic eruptions based on the distribution, generates characterization tasks, and sends them to the host. The plate displacement monitoring module is used to obtain the plate displacement trajectory based on plate fault information when the crustal movement cycle is plate displacement, calculate the deviation value based on the plate displacement trajectory, and generate geological structure change information based on the deviation value.

7. A geological structure characterization system based on cross-gradient constraints according to claim 6, characterized in that, The control variable adjustment module includes: The crustal movement cycle acquisition unit is used to acquire the time node containing the eruption when the crustal movement cycle is a volcanic eruption. The parameter adjustment unit is used to adjust the corresponding parameters in the preset cross gradient model according to the crustal movement cycle; The geological structure formation cycle calculation unit is used to obtain the geological state change time of the volcanic eruption and calculate the geological structure formation cycle based on the geological state change time and stratigraphic evolution time. A control variable determination unit is used to determine control variables for corresponding parameters based on the formation cycle of the geological structure; the control variables include Earth's own factors and extraterrestrial factors. The factor control unit is used to adjust the Earth's own factors when the geological rock distribution characteristics detect volcanic rock information of the volcanic eruption.

8. A geological structure characterization system based on cross-gradient constraints according to claim 7, characterized in that, The plate displacement monitoring module includes: The image feature acquisition unit is used to acquire image features of rock strata when the crustal movement cycle is plate displacement; the image features are image sub-blocks containing rock strata. The characterization selection unit is used to select fracture information of the plate to be inspected based on the plate displacement position, and to select fracture information of the target plate from the fracture information of the plate to be inspected based on the image features. The geological structure-time acquisition unit is used to determine the geological structure-time correlation based on the location of the fault in the target plate. The plate position dynamic unit is used to update the plate displacement position according to the aforementioned correlation relationship and executes it cyclically; The plate displacement generation unit is used to statistically analyze and connect all plate displacement positions, and generate plate displacement trajectories in the cross gradient model.

Images