Sequence configuration evolution dynamic simulation method

Abstract

The invention relates to the technical field of evolution dynamic simulation, in particular to a sequence configuration evolution dynamic simulation method. The method comprises the following steps: collecting geological data of a research area, determining a structural evolution sequence of the research area and basic geological data of the research area based on the geological data, and constructing an initial stratum geological model; based on the research area tectonic evolution sequence, performing multi-stage tectonic deformation numerical simulation on the initial stratigraphic geologic model to generate a sequence stratigraphic spatial configuration model; performing seismic forward modeling based on the sequence stratigraphic spatial configuration model to generate a forward modeling seismic data volume; comparing the forward modeling seismic data volume with the acquired actual seismic data volume, and performing iterative correction on simulation parameters of multi-period tectonic deformation numerical simulation according to a comparison result; according to the invention, through sequence configuration evolution dynamic simulation, a sequence configuration dynamic evolution deformation process can be accurately determined.

Description

Technical Field

[0001] This invention relates to the field of evolutionary dynamic simulation technology, and in particular to a method for dynamic simulation of sequence configuration evolution. Background Technology

[0002] The dynamic evolution and deformation process of sequence stratigraphy is a core foundation for understanding basin formation and resource distribution patterns, and has significant guiding significance for the exploration of oil, gas, and other mineral resources. Existing technologies for simulating sequence stratigraphy evolution mostly focus on static simulations of sedimentary processes or single tectonic events, or analyze sequence formation mechanisms through basin subsidence and infill processes. However, existing technologies have significant shortcomings: firstly, they cannot fully reproduce the continuous alteration process of sequence strata by multiple phases of tectonic deformation (such as multi-stage stress, uplift, erosion, and subsidence) across different geological periods, and cannot accurately depict the dynamic evolution trajectory of sequence strata from their early initial state to their current configuration; secondly, there is a lack of effective verification methods to match simulation results with actual geological conditions, making it difficult to guarantee the reliability and accuracy of simulation results, and failing to provide accurate theoretical support for subsequent geological exploration. Therefore, there is an urgent need for a method that can fully simulate multiple phases of tectonic deformation processes and verify simulation results with actual data, thereby accurately clarifying the dynamic evolution and deformation process of sequence stratigraphy. Summary of the Invention

[0003] Therefore, it is necessary to provide a dynamic simulation method for the evolution of hierarchical configurations to solve at least one of the above-mentioned technical problems.

[0004] To achieve the above objectives, a dynamic simulation method for sequence configuration evolution includes the following steps:

[0005] Step S1: Collect geological data of the study area, determine the tectonic evolution sequence and basic geological data of the study area based on the geological data, and construct an initial stratigraphic geological model;

[0006] Step S2: Based on the tectonic evolution sequence of the study area, perform multi-stage tectonic deformation numerical simulations on the initial stratigraphic geological model to generate a sequence stratigraphic spatial configuration model;

[0007] Step S3: Perform forward seismic modeling based on the sequence stratigraphic spatial configuration model to generate forward seismic data volume;

[0008] Step S4: Compare the forward modeling seismic data volume with the acquired actual seismic data volume, and iteratively correct the simulation parameters of the multi-stage tectonic deformation numerical simulation based on the comparison results. When the error between the forward modeling seismic data volume and the actual seismic data volume meets the preset threshold, output the sequence configuration dynamic evolution process data of the study area throughout the geological period.

[0009] The beneficial effects of this invention are as follows: by reproducing key aspects such as stress loading, uplift and erosion, and subsidence in the tectonic evolution process in stages, it can completely simulate the continuous modification process of sequence strata by multiple tectonic deformations, breaking through the limitations of single-stage simulation or static simulation in existing technologies, and realizing the accurate reproduction of the dynamic evolution process of sequence configuration.

[0010] By introducing a step of comparing and verifying seismic forward modeling with actual seismic data, and achieving dynamic optimization of simulation parameters through quantitative comparison, the reliability and accuracy of sequence evolution simulation results are significantly improved, solving the problem of the lack of effective verification methods for existing simulation results.

[0011] It can clarify the correspondence between tectonic deformation parameters (stress direction and magnitude) and sequence configuration changes in different geological periods, providing a precise theoretical basis for a deeper understanding of the tectonic evolution mechanism and sequence development law in the study area, and has important guiding value for the exploration and deployment of oil, gas and other mineral resources. Attached Figure Description

[0012] Figure 1 This is a flowchart illustrating the steps of a dynamic simulation method for sequence configuration evolution.

[0013] Figure 2 for Figure 1 A detailed flowchart illustrating the implementation steps of step S2.

[0014] Figure 3 This is a schematic diagram of geological testing.

[0015] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0016] The technical method of the present invention will now be clearly and completely described with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without inventive effort are within the scope of protection of the present invention.

[0017] Furthermore, the accompanying drawings are merely illustrative of the invention and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and therefore repeated descriptions of them will be omitted. Some block diagrams shown in the drawings are functional entities and do not necessarily correspond to physically or logically independent entities. These functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor methods and / or microcontroller methods.

[0018] It should be understood that although the terms "first," "second," etc., may be used herein to describe various units, these units should not be limited by these terms. These terms are used merely to distinguish one unit from another. For example, without departing from the scope of the exemplary embodiments, a first unit may be referred to as a second unit, and similarly, a second unit may be referred to as a first unit. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.

[0019] To achieve the above objectives, please refer to Figures 1 to 3 A method for dynamic simulation of sequence configuration evolution includes the following steps:

[0020] Step S1: Collect geological data of the study area, determine the tectonic evolution sequence and basic geological data of the study area based on the geological data, and construct an initial stratigraphic geological model;

[0021] Step S2: Based on the tectonic evolution sequence of the study area, perform multi-stage tectonic deformation numerical simulations on the initial stratigraphic geological model to generate a sequence stratigraphic spatial configuration model;

[0022] Step S3: Perform forward seismic modeling based on the sequence stratigraphic spatial configuration model to generate forward seismic data volume;

[0023] Step S4: Compare the forward modeling seismic data volume with the acquired actual seismic data volume, and iteratively correct the simulation parameters of the multi-stage tectonic deformation numerical simulation based on the comparison results. When the error between the forward modeling seismic data volume and the actual seismic data volume meets the preset threshold, output the sequence configuration dynamic evolution process data of the study area throughout the geological period.

[0024] All specific values ​​involved in this embodiment are exemplary parameters used to clearly illustrate the technical operation process and are not the only limitation of the present invention.

[0025] In one embodiment, the study area is a basin, and the collected geological data includes drilling data, well logging data, seismic profile data, core analysis data, paleontological fossil data, and structural geological maps. The collected geological data undergoes data processing and preprocessing, specifically including: standardizing and unifying the data format, converting data from different sources into a standardized format; performing quality checks and outlier removal, identifying and processing outliers and missing values; performing interpolation and spatial reconstruction on the processed geological data, using Kriging interpolation to spatially interpolate discrete data to form a continuous spatial data field; and establishing a data index to generate geological data suitable for geological modeling. Based on the processed geological data, the tectonic evolution sequence of the study area is determined, including identifying tectonic events, determining tectonic phases, and obtaining parameters such as stress field characteristics, erosion rates, and subsidence rates for each tectonic event; simultaneously, basic geological data is determined, including stratigraphic lithology, thickness, distribution range, paleogeographic background, and initial tectonic environment parameters. Based on the tectonic evolution sequence and basic geological data, an initial stratigraphic geological model of the study area was constructed using three-dimensional geological modeling technology. Specifically, this included: establishing a three-dimensional geological grid structure for the study area based on the basic geological data, with the grid size set to 100m × 100m × 10m according to the required research accuracy; determining the initial positions of each sequence boundary according to the tectonic evolution sequence and projecting the sequence boundaries onto the three-dimensional grid; and assigning lithological and physical-mechanical parameters, including density, elastic modulus, and Poisson's ratio, to each sequence unit within the three-dimensional geological grid structure, thus forming the initial stratigraphic geological model. This model can accurately characterize the spatial distribution characteristics and physical-mechanical parameters of the sequence stratigraphy in the early stages of the study area before significant tectonic deformation.

[0026] Based on the tectonic evolution sequence of the study area, the number of multiple tectonic deformation stages and the deformation parameters of each stage were determined. Specifically, this included: extracting the tectonic event sequence from the tectonic evolution sequence of the study area to determine the number of multiple tectonic deformation stages; in this embodiment, the study area experienced three tectonic deformation stages; based on the tectonic event sequence, obtaining the stress field characteristics, erosion rate, and subsidence rate parameters of each tectonic event to form the deformation parameters of each stage in the multiple tectonic deformation stages, including stress direction, stress magnitude, erosion rate, and subsidence rate. Based on the number of multiple tectonic deformation stages and the deformation parameters of each stage, a stage-by-stage tectonic deformation numerical simulation was performed on the initial stratigraphic geological model. Specifically, this included: determining the stage-by-stage tectonic deformation numerical simulation process based on the number of multiple tectonic deformation stages and the deformation parameters of each stage; and performing a stage-by-stage tectonic deformation numerical simulation on the initial stratigraphic geological model based on the stage-by-stage tectonic deformation numerical simulation process. The simulation process employs the finite element method (FEM), specifically including: First-stage stress loading and deformation simulation: Based on the tectonic stress field characteristics of the corresponding geological period, a matching stress direction and magnitude are applied to the initial stratigraphic geological model. The deformation process of the strata under this stress is simulated using the FEM, resulting in the first-stage deformed stratigraphic model. Second-stage uplift and erosion simulation: Based on the tectonic evolution sequence, uplift is applied to the first-stage deformed stratigraphic model. Combined with erosion rate and extent data of the study area, the uplift and erosion process of the strata is simulated, removing the eroded material. The eroded strata are used to obtain a model of the uplifted and eroded strata. Subsidence and secondary stress loading deformation simulations are performed by applying subsidence to the uplifted and eroded strata model to simulate the subsidence process. Then, based on the tectonic stress field characteristics of this geological period, the corresponding stress direction and magnitude are applied again, and the secondary deformation process of the strata is simulated through numerical calculations. This process is repeated in multiple cycles until the tectonic evolution process of all geological periods in the study area is simulated, ultimately obtaining a simulation model representing the current sequence stratigraphic spatial configuration. Stratigraphic deformation data from each stage of the phased tectonic deformation numerical simulation are obtained. Specifically, this includes: obtaining the simulation results for each stage of the phased tectonic deformation numerical simulation process and extracting the stratigraphic deformation data from each stage's simulation results, including stratigraphic deformation, stratigraphic displacement, and stratigraphic strain; integrating the extracted stratigraphic deformation data from each stage's simulation results according to the chronological order of the corresponding deformation stages to obtain the stratigraphic deformation data for each stage's simulation results. Based on the stratigraphic deformation data from each stage of simulation, the results of multiple phases of tectonic deformation numerical simulation are integrated to generate a sequence stratigraphic spatial configuration model. This model characterizes the spatial distribution characteristics and tectonic morphology of the current sequence stratigraphy in the study area.

[0027] Based on the sequence stratigraphic spatial configuration model, the location data of each sequence interface and the physical parameter data of each sequence unit are extracted, and the seismic impedance value of each sequence unit is calculated. Specifically, this includes: extracting the spatial location data of each sequence interface from the sequence stratigraphic spatial configuration model, including the three-dimensional coordinate information of the sequence interface; extracting the physical parameter data of each sequence unit, including parameters such as density and velocity; and calculating the seismic impedance value of each sequence unit based on the density and velocity parameters. The impedance value is the product of density and velocity. Based on the seismic impedance values ​​of each sequence unit, forward modeling of the seismic wavefield is performed to generate a forward model seismic data volume. Specifically, this includes: using the wave equation method for forward seismic modeling, setting the seismic wavelet type to the Ricker wavelet with a dominant frequency of 30Hz; setting observation system parameters, including shot-to-point spacing, receiver-to-point spacing, and record length; constructing a seismic impedance model based on the seismic impedance values ​​of each sequence unit; performing forward modeling of the seismic wavefield using the wave equation numerical simulation method to simulate the propagation process of seismic waves in the medium; and generating and processing the seismic data volume based on the forward modeling results, including stacking and migration, to obtain the forward model seismic data volume. This data volume is consistent with the format of actual seismic acquisition data and includes information such as the reflection characteristics of sequence boundaries, spatial variations in stratigraphic thickness, and structural morphology.

[0028] The generated forward-modeled seismic data volume is compared and analyzed with the actual seismic data volume currently acquired in the study area. The comparison indicators include the reflection characteristics of sequence boundaries, the spatial variation of stratigraphic thickness, and the degree of agreement of structural morphology. The comparison error is calculated using quantitative indicators such as correlation coefficient and root mean square error. The preset threshold is a correlation coefficient greater than 0.8 and a root mean square error less than 0.1. When the comparison error exceeds the preset threshold, the simulation parameters of the multi-stage tectonic deformation numerical simulation are adjusted according to the comparison results, and the process returns to step S2 to re-perform the multi-stage tectonic deformation numerical simulation. Specifically, this includes: analyzing the error sources based on the comparison results, determining the simulation parameters that need to be adjusted, including stress parameters, uplift and erosion parameters, or subsidence parameters; adjusting the simulation parameters for the corresponding geological period, re-performing the multi-stage tectonic deformation numerical simulation, and generating a new sequence stratigraphic spatial configuration model; re-performing the forward seismic simulation based on the new sequence stratigraphic spatial configuration model, and generating a new forward-modeled seismic data volume; comparing the new forward-modeled seismic data volume with the actual seismic data volume to obtain the new comparison error and determine whether it meets the preset threshold. When the comparison error is less than or equal to a preset threshold, model data from each stage of multi-stage tectonic deformation numerical simulation are extracted and integrated to generate dynamic evolution data of sequence configuration throughout the geological period of the study area. Specifically, this includes: extracting stratigraphic model data from each stage of multi-stage tectonic deformation numerical simulation, including stratigraphic deformation, stratigraphic displacement, and stratigraphic strain at each stage; integrating deformation parameters, uplift and erosion, and subsidence information from various geological periods; constructing a complete dynamic evolution trajectory of sequence stratigraphy from its early initial state to its current configuration, clarifying the modification mechanism and influence law of tectonic deformation on sequence configuration in different geological periods; and finally obtaining dynamic evolution deformation process data of the current sequence configuration throughout the geological period of the study area. This data includes stratigraphic models, deformation parameters, evolution trajectories, and other information at each stage, which can be used for the study of sequence configuration evolution laws and oil and gas exploration deployment.

[0029] Please refer to [link / reference needed] for further information. Figure 3 Drilling is a key method for obtaining underground geological information by drilling through multiple geological structures. Drilling involves obtaining basic geological data such as stratigraphic lithology, thickness, and distribution. The collected data is used to determine the tectonic evolution sequence, construct initial stratigraphic geological models, and conduct multi-stage numerical simulations of tectonic deformation.

[0030] Preferably, step S1 includes the following steps:

[0031] Step S11: Collect geological data for the study area;

[0032] Step S12: Organize and preprocess the geological data to obtain geological data for geological modeling;

[0033] Step S13: Determine the tectonic evolution sequence and basic geological data of the study area based on geological data;

[0034] Step S14: Based on the tectonic evolution sequence and basic geological data, construct the initial stratigraphic geological model of the study area using three-dimensional geological modeling technology.

[0035] In this embodiment, the study area is a basin, and the collected geological data includes drilling data, well logging data, seismic profile data, core analysis data, paleontological fossil data, and structural geological maps. Drilling data includes well location coordinates, well depth, stratigraphic layering, and lithological descriptions; well logging data includes natural gamma ray, sonic transit time, density, and resistivity logging curves; seismic profile data includes two-dimensional seismic profiles and three-dimensional seismic data volumes; core analysis data includes core descriptions, thin section identification, and physical property analysis; paleontological fossil data is used to determine stratigraphic age and sedimentary environment; and structural geological maps include structural outline maps, fault distribution maps, and structural profile maps. This geological data provides fundamental data support for subsequent determination of tectonic evolution sequences and geological modeling.

[0036] The collected geological data undergoes data organization and preprocessing, specifically including: standardizing and unifying the data format, converting data from different sources into standard formats, such as converting drilling data into standard tabular data and well logging data into standard curve data; performing quality checks and outlier removal, identifying and processing outliers and missing values ​​in the data, such as removing abnormal well logging values ​​and interpolating missing data; performing interpolation and spatial reconstruction on the processed geological data, using the Kriging interpolation method to spatially interpolate discrete data to form a continuous spatial data field; and establishing a data index to generate geological data that can be used for geological modeling, including stratigraphic interface data, lithological data, and physical property data, providing basic data for 3D geological modeling.

[0037] Based on the processed geological data, the tectonic evolution sequence and basic geological data of the study area were determined. Determining the tectonic evolution sequence included: identifying tectonic events such as fault activity, folding deformation, uplift and erosion, and subsidence and deposition; determining tectonic phases by identifying multiple tectonic deformation stages based on the chronological order and superposition of tectonic events; and acquiring parameters such as stress field characteristics, erosion rate, and subsidence rate for each tectonic event to form the tectonic evolution sequence. Determining the basic geological data included: determining stratigraphic lithology, including sandstone, mudstone, and carbonate rocks; determining stratigraphic thickness, including the thickness distribution of each sequence unit; determining the stratigraphic distribution range, including the spatial distribution of each sequence unit; determining the paleogeographic background, including information on sedimentary environment, paleowater depth, and paleoclimate; and determining initial tectonic environment parameters, including initial stress field and initial temperature field.

[0038] Based on the tectonic evolution sequence and basic geological data, an initial stratigraphic geological model of the study area was constructed using three-dimensional geological modeling technology. Specifically, this included: establishing a three-dimensional geological grid structure for the study area based on the basic geological data; setting the grid size to 100m × 100m × 10m according to the required research accuracy; determining the initial positions of each sequence boundary according to the tectonic evolution sequence; projecting the sequence boundaries onto the three-dimensional grid to form the spatial distribution of the sequence boundaries; and assigning lithological and physical-mechanical parameters, including density, elastic modulus, and Poisson's ratio, to each sequence unit in the three-dimensional geological grid structure to form the initial stratigraphic geological model. This model can accurately characterize the spatial distribution characteristics and physical-mechanical parameters of the sequence stratigraphy in the early stages of the study area before significant tectonic deformation, providing an initial model for subsequent multi-stage tectonic deformation numerical simulations.

[0039] Preferably, step S12 includes the following steps:

[0040] Step S121: Standardize and unify the data format of the geological data, and perform quality checks and outlier removal;

[0041] Step S122: Perform interpolation and spatial reconstruction on the processed geological data, and establish a data index to generate geological data.

[0042] In this embodiment, the collected geological data undergoes data format unification and standardization processing. Specifically, this includes: converting data from different sources into standard formats, such as converting drilling data into standard tabular data, including well location coordinates, well depth, stratigraphic layers, and lithological descriptions; converting logging data into standard curve data, including logging curves for natural gamma, sonic transit time, density, and resistivity; converting seismic profile data into standard seismic data volumes; and converting core analysis data into standard physical property data. Quality checks and outlier removal are then performed, including: identifying and processing outliers and missing values ​​in the data, such as removing abnormal logging values ​​and interpolating missing data; checking data consistency, such as verifying the consistency of stratigraphic layers and the accuracy of lithological descriptions; and checking data completeness, such as checking for missing logging curves and missing seismic profiles. Through data format unification and standardization processing, and quality checks and outlier removal, the quality and consistency of the geological data are ensured, providing reliable foundational data for subsequent data processing.

[0043] The processed geological data undergoes interpolation and spatial reconstruction, specifically including: using Kriging interpolation to spatially interpolate discrete data, forming a continuous spatial data field; for example, interpolating stratigraphic interfaces in drilling data to form continuous stratigraphic interface surfaces; interpolating physical property parameters in well logging data to form a continuous physical property parameter field; and interpolating seismic attributes in seismic data to form a continuous seismic attribute field. Data indexing is then established, including: creating a stratigraphic interface data index, including the spatial location of sequence boundaries and the thickness of sequence units; creating a lithological data index, including the lithological type and distribution of each sequence unit; and creating a physical property data index, including physical and mechanical parameters such as density, elastic modulus, and Poisson's ratio. Geological data, including stratigraphic interface data, lithological data, and physical property data, is generated to provide foundational data for 3D geological modeling. Through interpolation, spatial reconstruction, and data indexing, discrete geological data is converted into a continuous spatial data field, and a data index is established to facilitate subsequent data retrieval and use.

[0044] Preferably, step S14 includes the following steps:

[0045] Step S141: Based on basic geological data, establish a three-dimensional geological grid structure for the study area, and determine the initial positions of each sequence boundary according to the tectonic evolution sequence;

[0046] Step S142: Based on the basic geological data, assign lithological and physical-mechanical parameters to each sequence unit in the three-dimensional geological grid structure to form an initial stratigraphic geological model.

[0047] In this embodiment, a three-dimensional geological grid structure for the study area is established based on fundamental geological data. This includes: determining the three-dimensional spatial extent of the study area, including the X, Y, and Z directions; setting the grid size according to the required research accuracy (in this embodiment, the grid size is set to 100m × 100m × 10m), with the number of grid nodes determined based on the study area's extent; establishing a three-dimensional geological grid structure, with each grid node storing geological attribute information. The initial positions of each sequence boundary are determined based on the tectonic evolution sequence. This includes: determining the depth or elevation of each sequence boundary based on stratigraphic layering information in the tectonic evolution sequence; projecting the sequence boundaries onto the three-dimensional grid to form their spatial distribution; determining the initial morphology of the sequence boundaries, such as horizontal bedding or diagonal bedding, based on the paleogeographic background and sedimentary environment; and determining the initial stress state of the sequence boundaries based on initial tectonic environment parameters. By establishing the three-dimensional geological grid structure and determining the initial positions of each sequence boundary, a spatial framework is provided for assigning subsequent lithological and physical-mechanical parameters.

[0048] Based on fundamental geological data, lithological and physical-mechanical parameters are assigned to each sequence stratigraphic unit within a three-dimensional geological grid structure. Specifically, this includes: assigning lithological types (such as sandstone, mudstone, and carbonate rocks) to each sequence stratigraphic unit based on lithological data; assigning physical-mechanical parameters (including density, elastic modulus, and Poisson's ratio) to each sequence stratigraphic unit based on physical property data; determining the wave impedance value of each sequence stratigraphic unit based on lithological type and physical property parameters (the wave impedance value is the product of density and velocity); determining the initial stress state of each sequence stratigraphic unit based on the tectonic evolution sequence, including the magnitude and direction of initial stress; and determining the initial temperature field of each sequence stratigraphic unit based on the paleogeographic background. This forms an initial stratigraphic geological model, which includes information such as the three-dimensional geological grid structure, the spatial distribution of sequence boundaries, the lithological type, physical-mechanical parameters, wave impedance values, initial stress state, and initial temperature field of each sequence stratigraphic unit. This model can accurately characterize the spatial distribution characteristics and physical-mechanical parameters of sequence stratigraphy in the early stages of the study area before significant tectonic deformation, providing an initial model for subsequent multi-stage tectonic deformation numerical simulations.

[0049] Preferably, step S2 includes the following steps:

[0050] Step S21: Based on the tectonic evolution sequence of the study area, determine the number of multiple tectonic deformation stages and the deformation parameters of each stage;

[0051] Step S22: Based on the number of multiple tectonic deformation stages and the deformation parameters of each stage, perform a stage-by-stage tectonic deformation numerical simulation on the initial stratigraphic geological model;

[0052] Step S23: Obtain the simulation results of each stage in the stage-by-stage structural deformation numerical simulation, and record the formation deformation data of each stage simulation result;

[0053] Step S24: Based on the stratigraphic deformation data from each stage of simulation, integrate the numerical simulation results of tectonic deformation from multiple phases to generate a sequence stratigraphic spatial configuration model.

[0054] In this embodiment, based on the tectonic evolution sequence of the study area, the number of multiple tectonic deformation stages and the deformation parameters of each stage are determined. Specifically, this includes: extracting the tectonic event sequence from the tectonic evolution sequence of the study area to determine the number of multiple tectonic deformation stages. In this embodiment, the study area experienced three tectonic deformation stages: the first, second, and third tectonic deformation stages; based on the tectonic event sequence, obtaining the stress field characteristics, erosion rate, and subsidence rate parameters of each tectonic event to form the deformation parameters of each stage in the multiple tectonic deformation stages, including parameters such as stress direction, stress magnitude, erosion rate, and subsidence rate; determining the duration and deformation rate of each tectonic event according to the tectonic evolution sequence; and determining the temperature field and pressure field parameters of each tectonic event according to the paleogeographic background and sedimentary environment. By determining the number of multiple tectonic deformation stages and the deformation parameters of each stage, parameter input is provided for subsequent stage-by-stage numerical simulation of tectonic deformation.

[0055] Based on the number of tectonic deformation stages and the deformation parameters of each stage, a stage-by-stage numerical simulation of tectonic deformation is performed on the initial stratigraphic geological model. Specifically, this includes: determining the stage-by-stage numerical simulation process based on the number of tectonic deformation stages and the deformation parameters of each stage, including steps such as stress loading, uplift and erosion, subsidence, and secondary stress loading; performing stage-by-stage numerical simulations of tectonic deformation on the initial stratigraphic geological model based on this process, using the finite element method (FEM) to set boundary conditions, initial conditions, and material parameters; and simulating the first stage, stress loading and deformation, by applying matched stress directions and magnitudes to the initial stratigraphic geological model according to the tectonic stress field characteristics of the corresponding geological period, and simulating the deformation process of the strata under this stress using the FEM. The process involves obtaining a stratigraphic model after the first stage of deformation; uplift and erosion simulation: based on the tectonic evolution sequence, uplift is applied to the stratigraphic model after the first stage of deformation. Combining the erosion rate and erosion range data of the study area, the uplift and erosion process of the strata is simulated, and the eroded strata are removed to obtain the uplifted and eroded stratigraphic model; subsidence and secondary stress loading deformation simulation: subsidence is applied to the uplifted and eroded stratigraphic model to simulate the subsidence process of the strata. Subsequently, based on the tectonic stress field characteristics of this geological period, the corresponding stress direction and magnitude are applied again, and the secondary deformation process of the strata is simulated through numerical calculation; multiple cycles are repeated, repeating the above steps of uplift and erosion, subsidence, and stress loading deformation, until the tectonic evolution process of all geological periods in the study area is simulated, and finally, a simulation model representing the current sequence stratigraphic spatial configuration is obtained.

[0056] This process involves acquiring simulation results for each stage of a phased tectonic deformation numerical simulation, and recording the stratigraphic deformation data for each stage. Specifically, this includes: acquiring the simulation results for each stage of the phased tectonic deformation numerical simulation process, including the stratigraphic model after the first stage deformation, the stratigraphic model after uplift and erosion, and the stratigraphic model after subsidence and secondary stress loading deformation; extracting stratigraphic deformation data from the simulation results for each stage, including stratigraphic deformation, stratigraphic displacement, and stratigraphic strain; integrating the extracted stratigraphic deformation data from each stage simulation results according to the chronological order of the corresponding deformation stages to obtain the stratigraphic deformation data for each stage simulation result; and recording the stratigraphic deformation data for each stage simulation result, including storing it as a data file and establishing a data index, to facilitate subsequent data retrieval and use. By acquiring and recording the stratigraphic deformation data from each stage simulation result, data support is provided for the subsequent generation of a sequence stratigraphic spatial configuration model.

[0057] Based on the stratigraphic deformation data from each stage of simulation, a sequence stratigraphic spatial configuration model is generated by integrating numerical simulation results from multiple phases of tectonic deformation. Specifically, this includes: integrating stratigraphic deformation data from each stage of simulation, including data on stratigraphic deformation, displacement, and strain; generating a sequence stratigraphic spatial configuration model based on the integrated data, which characterizes the spatial distribution and tectonic morphology of the current sequence stratigraphy in the study area, including the spatial location of sequence boundaries, the thickness of sequence units, the lithology of sequence units, and their physical and mechanical parameters; post-processing the sequence stratigraphic spatial configuration model, including data smoothing, data interpolation, and data visualization, to improve the model's quality and readability; and outputting the sequence stratigraphic spatial configuration model, including as a data file and as a graphic file, to facilitate subsequent seismic forward modeling and comparative verification. By integrating numerical simulation results from multiple phases of tectonic deformation, a sequence stratigraphic spatial configuration model is generated, providing a foundational model for subsequent seismic forward modeling.

[0058] Preferably, step S21 includes the following steps:

[0059] Step S211: Extract the tectonic event sequence from the tectonic evolution sequence of the study area to determine the number of multiple tectonic deformation stages;

[0060] Step S212: Based on the sequence of tectonic events, obtain the stress field characteristics, erosion rate, and settlement rate parameters of each tectonic event to form the deformation parameters of each stage in the multi-stage tectonic deformation process.

[0061] In this embodiment, the tectonic event sequence in the tectonic evolution sequence of the study area is extracted to determine the number of multiple tectonic deformation stages. Specifically, this includes: identifying tectonic events based on the tectonic evolution sequence of the study area, including faulting, folding deformation, uplift and erosion, and subsidence and deposition; determining the sequence of tectonic events according to their chronological order and superposition relationship, such as the first tectonic event being faulting, the second being uplift and erosion, and the third being subsidence and deposition; determining the number of multiple tectonic deformation stages according to the type and scale of the tectonic events, in this embodiment, the study area experienced three tectonic deformation stages, namely the first, second, and third tectonic deformation stages; determining the duration of each tectonic deformation stage according to the duration of the tectonic events, such as the first tectonic deformation stage lasting 10 Ma, the second lasting 5 Ma, and the third lasting 8 Ma; and determining the chronological order and superposition relationship of each tectonic deformation stage according to the superposition relationship of the tectonic events, such as the first tectonic deformation stage superimposed on the second tectonic deformation stage, and the second tectonic deformation stage superimposed on the third tectonic deformation stage. By extracting the sequence of structural events and determining the number of structural deformation stages in multiple phases, a foundation is provided for obtaining deformation parameters for each subsequent stage.

[0062] Based on the tectonic event sequence, the stress field characteristics, erosion rate, and subsidence rate parameters of each tectonic event are obtained to form deformation parameters for each stage in the multi-stage tectonic deformation process. Specifically, this includes: obtaining the stress field characteristics of each tectonic event based on the tectonic event sequence, including parameters such as stress direction, stress magnitude, and stress type (compression, tension, strike-slip). For example, the stress direction of the first tectonic event is northwest-southeast, the stress magnitude is 50 MPa, and the stress type is compression; the stress direction of the second tectonic event is northeast-southwest, the stress magnitude is 30 MPa, and the stress type is tension; the stress direction of the third tectonic event is near east-west, the stress magnitude is 20 MPa, and the stress type is strike-slip. The erosion rates of each tectonic event were obtained, including parameters such as the magnitude of the erosion rate, the erosion range, and the erosion duration. For example, the erosion rate of the first tectonic event was 100 m / Ma, the erosion range was the entire study area, and the erosion duration was 5 Ma; the erosion rate of the second tectonic event was 50 m / Ma, the erosion range was the eastern part of the study area, and the erosion duration was 3 Ma; the erosion rate of the third tectonic event was 80 m / Ma, the erosion range was the western part of the study area, and the erosion duration was 4 Ma. The subsidence rates of each tectonic event were also obtained, including parameters such as the magnitude of the subsidence rate, the subsidence range, and the subsidence duration. For example, the subsidence rate of the first tectonic event was 50 m / Ma, the subsidence range was the central part of the study area, and the subsidence duration was 8 Ma; the subsidence rate of the second tectonic event was 80 m / Ma, the subsidence range was the southern part of the study area, and the subsidence duration was 6 Ma; the subsidence rate of the third tectonic event was 60 m / Ma, the subsidence range was the northern part of the study area, and the subsidence duration was 7 Ma. Based on the obtained stress field characteristics, erosion rate, and settlement rate parameters, deformation parameters for each stage in the multi-stage tectonic deformation process are formed, including stress parameters, erosion parameters, and settlement parameters, providing parameter inputs for subsequent stage-by-stage numerical simulation of tectonic deformation.

[0063] Preferably, step S22 includes the following steps:

[0064] Step S221: Based on the number of structural deformation stages in multiple phases and the deformation parameters of each stage, determine the numerical simulation process for structural deformation stage by stage;

[0065] Step S222: Based on the step-by-step structural deformation numerical simulation process, perform step-by-step structural deformation numerical simulation on the initial stratigraphic geological model.

[0066] In this embodiment, based on the number of multiple tectonic deformation stages and the deformation parameters of each stage, a stage-by-stage numerical simulation process for tectonic deformation is determined. Specifically, this includes: determining the simulation stages based on the number of multiple tectonic deformation stages (in this embodiment, the study area experienced three tectonic deformation stages, therefore the simulation process includes three stages); determining the simulation steps for each stage based on the deformation parameters of each stage (each stage includes stress loading and deformation simulation, uplift and erosion simulation, and subsidence and secondary stress loading deformation simulation); determining the boundary conditions and initial conditions for the simulation, including fixed boundaries, free boundaries, stress boundaries, and initial stress states and initial temperature fields; determining the material parameters for the simulation, including the density, elastic modulus, and Poisson's ratio of each sequence element; determining the time step and number of iterations for the simulation (the time step is determined based on the duration of the tectonic event, and the number of iterations is determined based on the convergence accuracy); determining the output parameters for the simulation, including formation deformation data such as formation deformation, formation displacement, and formation strain; and determining the termination conditions for the simulation, including stress equilibrium and displacement convergence. By defining the stage-by-stage structural deformation numerical simulation process, we provide process guidance for subsequent stage-by-stage structural deformation numerical simulations.

[0067] Based on a stage-by-stage tectonic deformation numerical simulation process, a stage-by-stage tectonic deformation numerical simulation is performed on the initial stratigraphic geological model. Specifically, this includes: Stage 1: Stress loading and deformation simulation. Based on the tectonic stress field characteristics of the corresponding geological period, a matching stress direction and magnitude are applied to the initial stratigraphic geological model. The deformation process of the strata under this stress is simulated using the finite element method, resulting in the stratigraphic model after the first stage of deformation. Stage 2: Uplift and erosion simulation. According to the tectonic evolution sequence, uplift is applied to the stratigraphic model after the first stage of deformation. Combined with erosion rate and erosion extent data of the study area, the erosion process is simulated. The simulation process involves uplift and erosion of the strata, removing the eroded portion to obtain a model of the uplifted and eroded strata. Subsidence and secondary stress loading deformation simulations are then performed, applying subsidence to the uplifted and eroded strata model to simulate the subsidence process. Subsequently, based on the tectonic stress field characteristics of this geological period, corresponding stress direction and magnitude are applied again, and the secondary deformation process of the strata is simulated through numerical calculations. This process is repeated in multiple cycles until the tectonic evolution process of all geological periods in the study area is simulated, ultimately obtaining a simulation model representing the current spatial configuration of the sequence stratigraphy. The simulation process employs the finite element method (FEM), setting boundary conditions, initial conditions, and material parameters. Iterative calculations are used to solve for strata deformation data such as stress field, displacement field, and strain field, ultimately obtaining the strata deformation data for each stage of the simulation.

[0068] Preferably, step S23 includes the following steps:

[0069] Step S231: Obtain the simulation results of each stage in the numerical simulation process of structural deformation, and extract the formation deformation data from the simulation results of each stage;

[0070] Step S232: Integrate the extracted formation deformation data from each stage of the simulation results according to the chronological order of the corresponding deformation stages, and record the formation deformation data from each stage of the simulation results.

[0071] In this embodiment, the simulation results of each stage in the phased tectonic deformation numerical simulation process are obtained, and the formation deformation data in each stage simulation result is extracted. Specifically, this includes: obtaining the simulation results of each stage in the phased tectonic deformation numerical simulation process, including the formation model after the first stage stress loading and deformation simulation, the formation model after the uplift and erosion simulation, and the formation model after the settlement and secondary stress loading deformation simulation; extracting the formation deformation data in each stage simulation result, including formation deformation, formation displacement, and formation strain, where formation deformation represents the degree of formation deformation, formation displacement represents the magnitude and direction of formation displacement, and formation strain represents the strain state of the formation; extracting relevant data from each stage simulation result, including stress field, displacement field, and strain field; performing quality checks on the extracted formation deformation data, including checking the completeness, consistency, and accuracy of the data to ensure data reliability; and converting the extracted formation deformation data into a standard format data file for easy subsequent data processing and use. By obtaining and extracting the formation deformation data from each stage simulation result, a data foundation is provided for subsequent data integration and processing.

[0072] The extracted formation deformation data from each stage of the simulation results are integrated and processed according to the chronological order of the corresponding deformation stages. This process includes: sorting the extracted formation deformation data according to the chronological order of the corresponding deformation stages, such as first-stage deformation data, second-stage deformation data, third-stage deformation data, etc.; integrating the sorted formation deformation data, including data alignment, data interpolation, and data smoothing, to ensure data consistency and continuity; generating formation deformation data files for each stage of the simulation results based on the integrated formation deformation data, including formation deformation files, formation displacement files, and formation strain files; recording the formation deformation data for each stage of the simulation results, including storing it as a data file, creating a data index, and generating a data report, to facilitate subsequent data retrieval and use; and verifying the recorded formation deformation data, including comparing it with the original simulation results and checking the completeness and accuracy of the data, to ensure data reliability. By integrating and recording the formation deformation data from each stage of the simulation results, data support is provided for the subsequent generation of sequence stratigraphic spatial configuration models.

[0073] Preferably, step S3 includes the following steps:

[0074] Step S31: Based on the sequence stratigraphic spatial configuration model, extract the location data of each sequence interface and the physical parameter data of each sequence element, and calculate the seismic impedance value of each sequence element;

[0075] Step S32: Based on the seismic wave impedance values ​​of each sequence element, perform forward modeling calculations of the seismic wavefield to generate forward modeling seismic data.

[0076] In this embodiment, based on the sequence stratigraphic spatial configuration model, the location data of each sequence interface and the physical parameter data of each sequence unit are extracted, and the seismic impedance value of each sequence unit is calculated. Specifically, this includes: extracting the spatial location data of each sequence interface from the sequence stratigraphic spatial configuration model, including the three-dimensional coordinate information of the sequence interface, the depth or elevation of the sequence interface, and the morphological characteristics of the sequence interface; extracting the physical parameter data of each sequence unit, including physical and mechanical parameters such as density, velocity, elastic modulus, and Poisson's ratio, with density and velocity parameters used to calculate the seismic impedance value; calculating the seismic impedance value of each sequence unit based on the density and velocity parameters, where the impedance value is the product of density and velocity, reflecting the formation's ability to reflect seismic waves; performing a quality check on the calculated seismic impedance value, including checking the data's completeness, consistency, and accuracy to ensure data reliability; and converting the seismic impedance value into a standard format data file to facilitate subsequent seismic forward modeling calculations. By extracting sequence boundary location data, physical parameter data, and calculating seismic wave impedance values, basic data are provided for subsequent forward modeling calculations of the seismic wavefield.

[0077] Based on the seismic impedance values ​​of each sequence unit, forward modeling of the seismic wavefield is performed to generate a forward modeled seismic data volume. Specifically, this includes: using the wave equation method for forward seismic modeling, setting the seismic wavelet type to the Ricker wavelet with a dominant frequency of 30Hz, and determining the wavelet length according to the required accuracy; setting observation system parameters, including shot-to-point spacing, receiver-to-point spacing, recording length, and sampling rate, with shot-to-point and receiver-to-point spacing determined based on the study area and accuracy requirements; constructing a seismic impedance model based on the seismic impedance values ​​of each sequence unit, including the distribution of wave impedance values ​​for each sequence unit; performing forward modeling of the seismic wavefield using the wave equation numerical simulation method to simulate the propagation process of seismic waves in the medium, including excitation, propagation, reflection, and transmission; and generating and processing the seismic data volume based on the forward modeling results, including stacking, migration, and filtering, to obtain the forward modeled seismic data volume. The forward modeled seismic data volume has the same format as the actual seismic acquisition data and includes information such as the reflection characteristics of sequence boundaries, spatial variations in stratigraphic thickness, and structural morphology. Forward modeling of seismic wavefields is used to generate forward modeled seismic data volumes, providing basic data for subsequent comparative verification.

[0078] Preferably, step S4 includes the following steps:

[0079] Step S41: Quantitatively compare the forward modeling seismic data volume with the actual seismic data volume, obtain the comparison error, and determine whether the comparison error is greater than the preset threshold.

[0080] Step S42: When the comparison error is greater than the preset threshold, adjust the simulation parameters of the multi-stage structural deformation numerical simulation according to the comparison results, and repeat the multi-stage structural deformation numerical simulation.

[0081] Step S43: When the comparison error is less than or equal to the preset threshold, extract the model data of each stage of the numerical simulation of tectonic deformation in multiple phases, and integrate them to generate the dynamic evolution process data of sequence configuration in the study area throughout the geological period.

[0082] In this embodiment, the forward model seismic data volume is quantitatively compared with the actual seismic data volume to obtain the comparison error. The comparison error is then determined to be greater than a preset threshold. Specifically, this includes: comparing the forward model seismic data volume generated in step 3 with the actual seismic data volume currently collected in the study area. Comparison indicators include the reflection characteristics of sequence boundaries, spatial variations in stratigraphic thickness, and the degree of agreement of structural morphology. The comparison error is calculated using quantitative indicators such as correlation coefficient and root mean square error (RMSE). The correlation coefficient reflects the similarity between the forward model seismic data volume and the actual seismic data volume, while the RMSE reflects the degree of difference between the two. The preset threshold is a correlation coefficient greater than 0.8 and an RMSE less than 0.1. When the correlation coefficient is greater than 0.8 and the RMSE is less than 0.1, the comparison error is considered to meet the preset threshold. When the correlation coefficient is less than or equal to 0.8 or the RMSE is greater than or equal to 0.1, the comparison error is considered to be greater than the preset threshold. Based on the magnitude of the comparison error, it is determined whether the simulation parameters need to be adjusted and multi-stage numerical simulations of tectonic deformation need to be performed again.

[0083] When the comparison error exceeds a preset threshold, the simulation parameters of the multi-stage tectonic deformation numerical simulation are adjusted based on the comparison results, and the multi-stage tectonic deformation numerical simulation is repeated. Specifically, this includes: analyzing the error sources based on the comparison results, determining the simulation parameters that need adjustment, including stress parameters, uplift and erosion parameters, or subsidence parameters; adjusting the simulation parameters for the corresponding geological period, such as adjusting stress direction, stress magnitude, erosion rate, and subsidence rate; repeating the multi-stage tectonic deformation numerical simulation to generate a new sequence stratigraphic spatial configuration model; based on the new sequence stratigraphic spatial configuration model, performing a new forward seismic modeling simulation to generate a new forward seismic data volume; comparing the new forward seismic data volume with the actual seismic data volume to obtain a new comparison error and determine whether it meets the preset threshold; repeating the above adjustment and simulation process until the comparison error between the forward seismic data volume and the actual seismic data volume meets the preset threshold. By iteratively correcting the simulation parameters, the accuracy and reliability of the simulation results are improved.

[0084] When the comparison error is less than or equal to a preset threshold, model data from each stage of the multi-stage tectonic deformation numerical simulation are extracted and integrated to generate dynamic evolution data of sequence configuration throughout the geological period of the study area. Specifically, this includes: extracting stratigraphic model data from each stage of the multi-stage tectonic deformation numerical simulation process, including stratigraphic deformation, stratigraphic displacement, and stratigraphic strain at each stage; integrating deformation parameters, uplift and erosion, and subsidence information from various geological periods; constructing a complete dynamic evolution trajectory of sequence stratigraphy from its early initial state to its current configuration, clarifying the modification mechanism and influence law of tectonic deformation on sequence configuration in different geological periods; finally obtaining dynamic evolution deformation process data of the current sequence configuration throughout the geological period of the study area, which includes stratigraphic models, deformation parameters, evolution trajectory, and other information at each stage, and can be used for sequence configuration evolution law research and oil and gas exploration deployment; and outputting dynamic evolution process data of sequence configuration throughout the geological period of the study area, including output as data files, output as graphic files, and output as report files, to facilitate subsequent data analysis and application.

[0085] The above description is merely a specific embodiment of the present invention, enabling those skilled in the art to understand or implement the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features of the invention herein.

Claims

1. A method for dynamic simulation of sequence configuration evolution, characterized in that, Includes the following steps: Step S1: Collect geological data of the study area, determine the tectonic evolution sequence and basic geological data of the study area based on the geological data, and construct an initial stratigraphic geological model; Step S2: Based on the tectonic evolution sequence of the study area, perform multi-stage tectonic deformation numerical simulations on the initial stratigraphic geological model to generate a sequence stratigraphic spatial configuration model; Step S3: Perform forward seismic modeling based on the sequence stratigraphic spatial configuration model to generate forward seismic data volume; Step S4: Compare the forward modeling seismic data volume with the acquired actual seismic data volume, and iteratively correct the simulation parameters of the multi-stage tectonic deformation numerical simulation based on the comparison results. When the error between the forward modeling seismic data volume and the actual seismic data volume meets the preset threshold, output the sequence configuration dynamic evolution process data of the study area throughout the geological period.

2. The method for dynamic simulation of sequence configuration evolution according to claim 1, characterized in that, Step S1 includes the following steps: Step S11: Collect geological data for the study area; Step S12: Organize and preprocess the geological data to obtain geological data for geological modeling; Step S13: Determine the tectonic evolution sequence and basic geological data of the study area based on geological data; Step S14: Based on the tectonic evolution sequence and basic geological data, construct the initial stratigraphic geological model of the study area using three-dimensional geological modeling technology.

3. The method for dynamic simulation of sequence configuration evolution according to claim 2, characterized in that, Step S12 includes the following steps: Step S121: Standardize and unify the data format of the geological data, and perform quality checks and outlier removal; Step S122: Perform interpolation and spatial reconstruction on the processed geological data, and establish a data index to generate geological data.

4. The method for dynamic simulation of sequence configuration evolution according to claim 2, characterized in that, Step S14 includes the following steps: Step S141: Based on basic geological data, establish a three-dimensional geological grid structure for the study area, and determine the initial positions of each sequence boundary according to the tectonic evolution sequence; Step S142: Based on the basic geological data, assign lithological and physical-mechanical parameters to each sequence unit in the three-dimensional geological grid structure to form an initial stratigraphic geological model.

5. The method for dynamic simulation of sequence configuration evolution according to claim 3, characterized in that, Step S2 includes the following steps: Step S21: Based on the tectonic evolution sequence of the study area, determine the number of multiple tectonic deformation stages and the deformation parameters of each stage; Step S22: Based on the number of multiple tectonic deformation stages and the deformation parameters of each stage, perform a stage-by-stage tectonic deformation numerical simulation on the initial stratigraphic geological model; Step S23: Obtain the simulation results of each stage in the stage-by-stage structural deformation numerical simulation, and record the formation deformation data of each stage simulation result; Step S24: Based on the stratigraphic deformation data from each stage of simulation, integrate the numerical simulation results of tectonic deformation from multiple phases to generate a sequence stratigraphic spatial configuration model.

6. The method for dynamic simulation of sequence configuration evolution according to claim 4, characterized in that, Step S21 includes the following steps: Step S211: Extract the tectonic event sequence from the tectonic evolution sequence of the study area to determine the number of multiple tectonic deformation stages; Step S212: Based on the sequence of tectonic events, obtain the stress field characteristics, erosion rate, and settlement rate parameters of each tectonic event to form the deformation parameters of each stage in the multi-stage tectonic deformation process.

7. The method for dynamic simulation of sequence configuration evolution according to claim 4, characterized in that, Step S22 includes the following steps: Step S221: Based on the number of structural deformation stages in multiple phases and the deformation parameters of each stage, determine the numerical simulation process for structural deformation stage by stage; Step S222: Based on the step-by-step structural deformation numerical simulation process, perform step-by-step structural deformation numerical simulation on the initial stratigraphic geological model.

8. The method for dynamic simulation of sequence configuration evolution according to claim 4, characterized in that, Step S23 includes the following steps: Step S231: Obtain the simulation results of each stage in the numerical simulation process of structural deformation, and extract the formation deformation data from the simulation results of each stage; Step S232: Integrate the extracted formation deformation data from each stage of the simulation results according to the chronological order of the corresponding deformation stages, and record the formation deformation data from each stage of the simulation results.

9. The method for dynamic simulation of sequence configuration evolution according to claim 1, characterized in that, Step S3 includes the following steps: Step S31: Based on the sequence stratigraphic spatial configuration model, extract the location data of each sequence interface and the physical parameter data of each sequence element, and calculate the seismic impedance value of each sequence element; Step S32: Based on the seismic wave impedance values ​​of each sequence element, perform forward modeling calculations of the seismic wavefield to generate forward modeling seismic data.

10. The method for dynamic simulation of sequence configuration evolution according to claim 1, characterized in that, Step S4 includes the following steps: Step S41: Quantitatively compare the forward modeling seismic data volume with the actual seismic data volume, obtain the comparison error, and determine whether the comparison error is greater than the preset threshold. Step S42: When the comparison error is greater than the preset threshold, adjust the simulation parameters of the multi-stage structural deformation numerical simulation according to the comparison results, and repeat the multi-stage structural deformation numerical simulation. Step S43: When the comparison error is less than or equal to the preset threshold, extract the model data of each stage of the numerical simulation of tectonic deformation in multiple phases, and integrate them to generate the dynamic evolution process data of sequence configuration in the study area throughout the geological period.

Images