CN110032788B

CN110032788B - Discrete element simulation method, device and system for plate diving deformation process

Info

Publication number: CN110032788B
Application number: CN201910271732.7A
Authority: CN
Inventors: 王一丹; 于福生; 王于恒; 王春英; 王丹丹; 吕旭阳; 王逸群; 杨金月; 王宏杰; 冯佳梦; 刘逸伦
Original assignee: China University of Petroleum Beijing
Current assignee: China University of Petroleum Beijing
Priority date: 2019-04-04
Filing date: 2019-04-04
Publication date: 2021-03-30
Anticipated expiration: 2039-04-04
Also published as: CN110032788A

Abstract

The invention provides a discrete element simulation method, a discrete element simulation system, computer equipment and a computer readable storage medium for a plate diving deformation process, and relates to the technical field of geological exploration. The method comprises the steps of constructing a theoretical experimental model of a plate diving deformation structure according to geological conditions of a zone to be researched; establishing a discrete element numerical simulation model based on the theoretical experimental model; and analyzing the deformation characteristics of the plate diving deformation based on the discrete element numerical simulation model. The method has the advantages of simple model construction mode, easy realization of simulation mode and good simulation result accuracy, and solves the technical problems that the cause research and analysis process of the geological structure cannot be verified, and an effective correlation means is lacked between theory and actual cause.

Description

Discrete element simulation method, device and system for plate diving deformation process

Technical Field

The invention relates to the technical field of geological exploration, in particular to a data simulation technology in geological exploration, and specifically relates to a discrete element simulation method, a system, computer equipment and a computer readable storage medium for a plate diving deformation process.

Background

A dive zone is formed when two rockring slabs collide with each other and one of the slabs descends below the other through the dive process. In the usual case, a dive panel is an ocean panel consisting of ocean hulls, which are more dense and more prone to sinking relative to the earth's hulls because they are made of a silicon-magnesium material. The diving effect can lead to the deformation of the upper disk, forming a diving wedge.

Due to different geological conditions and changing states, the development pattern of the fault layer in the diving wedge can be influenced. At present, modeling to simulate actual geological conditions and changes thereof is one of the main means adopted in fault formation mechanism research. The discrete element method is used as a numerical simulation method based on molecular dynamics, dynamic relaxation solution is carried out by using a center difference method, calculation is simple, convenient and quick, the method is a tool capable of researching mechanical characteristics and motion characteristics of a medium from a fine angle and a microscopic angle, and the discrete element method has particular advantages in the fields of the problem of rupture of a discontinuous medium and the problem of large deformation based on the characteristic of free motion of discrete particles, and is widely applied to structural geological research.

The simulation research on the plate diving deformation process still mainly reflects in the stage of analyzing the fracture development conditions in the stratums with different mechanical properties. For example, Hardy et al studied the development of bi-directional thrust wedges in the upper plate using the discrete element method; jonathan et al studied the development of proliferative wedges in the dive zone using a finite element method.

However, in the modeling and analyzing process of the plate diving deformation process, research on an experimental model and an evolution process of diving deformation is still incomplete, and further research on deformation characteristics of the plate diving deformation process is needed. Moreover, the construction process of the experimental model is generally complex, the related parameters are more, and the matching degree of the theoretical analysis result and the actual situation is difficult to be accurately judged due to the lack of a proper verification relation in the theoretical analysis.

Therefore, how to provide a new scheme to realize more accurate research on the experimental model and the evolution process of the dive deformation is a technical problem to be solved in the field.

Disclosure of Invention

In view of this, the embodiment of the present invention provides a discrete element simulation method, a system, a computer device, and a computer-readable storage medium for a sheet diving deformation process, where a model construction manner is simple, a simulation manner is easy to implement, a simulation result is good in accuracy, and technical problems that a cause research and analysis process of a geological structure cannot be verified, and an effective association means is lacking between a theoretical cause and an actual cause are solved.

One of the objectives of the present invention is to provide a discrete element simulation method for a sheet diving deformation process, which includes:

constructing a theoretical experimental model of a plate diving deformation structure according to the geological conditions of the zone to be researched;

establishing a discrete element numerical simulation model based on the theoretical experimental model;

and analyzing the deformation characteristics of the plate diving deformation based on the discrete element numerical simulation model.

Preferably, the theoretical experimental model for constructing the plate diving deformation structure according to the geological conditions of the zone to be researched comprises the following steps:

determining basic parameters based on the shells of land and ocean and the occurrence of a diving effect;

and establishing a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions.

Preferably, the basic parameters comprise the thickness of the land-crust stratum, the thickness of the ocean-crust stratum, the mechanical parameters of the land-crust stratum and the ocean-crust stratum, the dive amount, the dive speed and the dive angle; the setting principle comprises that the land crust takes a sandstone stratum with rigidity smaller than a first threshold value as a reference, the ocean crust takes a sandstone stratum with deformation resistance larger than a second threshold value and rigidity larger than a third threshold value as a reference, the ocean crust is provided with a movable base with a plurality of telescopic points and has the characteristic of diving in any direction.

Preferably, the establishing of the discrete element numerical simulation model based on the theoretical experimental model includes:

constructing the bottom and lateral boundaries of the diving plate and the land shell plate according to the theoretical experimental model;

setting rock mechanical parameters and substrate and boundary mechanical parameters of the discrete element numerical simulation model;

constructing a simulated stratum adaptive to geological conditions required by the plate diving deformation based on the bottom and the lateral boundary;

dividing the simulated stratum into small layers marked with different colors in equal thickness;

the ocean shell plates have a diving effect towards the land shells at a fixed speed and a specific angle.

Preferably, analyzing the deformation characteristic of the plate diving deformation based on the discrete element numerical simulation model comprises:

forming different simulation fault patterns by changing simulation process parameters of the discrete element numerical simulation model, and obtaining corresponding simulation data;

observing deformation characteristics of the plate depression deformation process corresponding to different depression angles;

and summarizing the differential deformation characteristics of the sheet diving deformation process under different angles based on the simulated formation and fault deformation development rules.

One of the objectives of the present invention is to provide a discrete element simulation system for a sheet diving deformation process, comprising:

the experimental model building module is used for building a theoretical experimental model of the plate diving deformation structure according to the geological conditions of the zone to be researched;

the simulation model building module is used for building a discrete element numerical simulation model based on the theoretical experimental model;

and the plate depression analysis module is used for analyzing the deformation characteristics of plate depression deformation based on the discrete element numerical simulation model.

Preferably, the experimental model building module comprises:

the basic parameter determining module is used for determining basic parameters based on the shells of land and ocean and the occurrence of the diving effect;

and the experimental model establishing module is used for establishing a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions.

Preferably, the simulation model building module includes:

the boundary construction module is used for constructing the bottom and lateral boundaries of the diving plate and the land shell plate according to the theoretical experimental model;

the mechanical parameter setting module is used for setting the rock mechanical parameters and the substrate and boundary mechanical parameters of the discrete element numerical simulation model;

the simulated stratum construction module is used for constructing a simulated stratum which is adaptive to geological conditions required by the plate diving deformation based on the bottom and the lateral boundary;

the simulated formation dividing module is used for dividing the simulated formation into small layers marked with different colors in equal thickness;

and the shell plate diving module is used for enabling the shell plate to generate a diving effect on the land shell at a fixed speed and a specific angle.

Preferably, the plate dive analysis module includes:

the simulation data determining module is used for forming different simulation fault patterns by changing simulation process parameters of the discrete element numerical simulation model and obtaining corresponding simulation data;

the denaturation characteristic observation module is used for observing deformation characteristics of the plate diving deformation process corresponding to different diving angles;

and the degeneration characteristic summarizing module is used for summarizing the differential deformation characteristics of the plate diving deformation process at different angles based on the simulated stratum and fault deformation development rule.

One of the objects of the present invention is to provide a computer apparatus comprising: the system comprises a processor and a storage device, wherein the storage device is used for storing a plurality of instructions, and the instructions are suitable for being loaded by the processor and executing a discrete element simulation method of a slab diving deformation process.

It is an object of the present invention to provide a computer-readable storage medium storing a computer program for performing a discrete element simulation method of a sheet dive deformation process.

The invention has the beneficial effects that the discrete element simulation method, the discrete element simulation system, the computer equipment and the computer readable storage medium for the plate diving deformation process are provided, the model construction mode is simple, the simulation mode is easy to realize, the simulation result accuracy is good, and the technical problems that the cause research and analysis process of the geological structure cannot be verified, and an effective correlation means is lacked between a theory cause and an actual cause are solved.

In order to make the aforementioned and other objects, features and advantages of the invention comprehensible, preferred embodiments accompanied with figures are described in detail below.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.

Fig. 1 is a schematic structural diagram of a discrete element simulation system of a sheet diving deformation process according to an embodiment of the present invention;

fig. 2 is a schematic structural diagram of an experimental model building module in a discrete element simulation system of a sheet diving deformation process according to an embodiment of the present invention;

fig. 3 is a schematic structural diagram of a simulation model building module in a discrete element simulation system of a sheet diving deformation process according to an embodiment of the present invention;

fig. 4 is a schematic structural diagram of a plate dive analysis module in a discrete element simulation system of a plate dive deformation process according to an embodiment of the present invention;

FIG. 5 is a flowchart of a discrete element simulation method for a sheet diving deformation process according to an embodiment of the present invention;

fig. 6 is a detailed flowchart of step S101 in fig. 5;

fig. 7 is a detailed flowchart of step S102 in fig. 5;

fig. 8 is a detailed flowchart of step S103 in fig. 5;

fig. 9(a1) to 9(a2) are schematic diagrams of simulation devices of a discrete element simulation method of a plate diving deformation process with a diving angle of 30 ° for a ocean hull according to an embodiment of the present invention;

fig. 9(b1) to 9(b2) are schematic diagrams of simulation devices of a discrete element simulation method of a plate diving deformation process with a 15 ° diving angle of the ocean hull according to an embodiment of the present invention;

fig. 9(c1) to 9(c2) are schematic diagrams of simulation devices of a discrete element simulation method of a plate diving deformation process with a diving angle of 5 ° for a ocean hull according to an embodiment of the present invention;

fig. 10(a) to 10(f) are cross-sectional evolutions of simulation results of a discrete element simulation method of a plate diving deformation process with a 30 ° diving angle of the ocean hull according to an embodiment of the present invention;

11(a) to 11(f) are cross-sectional evolutions of simulation results of a discrete element simulation method for a plate diving deformation process with a 15 ° hull diving angle according to an embodiment of the present invention;

fig. 12(a) to 12(f) are cross-sectional evolution diagrams of simulation results of a discrete element simulation method of a plate diving deformation process with a diving angle of the ocean shell of 5 ° according to an embodiment of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

As will be appreciated by one skilled in the art, embodiments of the present invention may be embodied as a system, apparatus, method or computer program product. Accordingly, the present disclosure may be embodied in the form of: entirely hardware, entirely software (including firmware, resident software, micro-code, etc.), or a combination of hardware and software.

The principles and spirit of the present invention are explained in detail below with reference to several representative embodiments of the invention.

The invention provides a discrete element simulation system for a plate diving deformation process, which aims to solve the technical problems that the cause research and analysis process of a geological structure cannot be verified, and an effective association means is lacked between a theory and an actual cause.

Specifically, fig. 1 is a schematic structural diagram of a discrete element simulation system of a plate diving deformation process according to an embodiment of the present invention, please refer to fig. 1, where the discrete element simulation system of the plate diving deformation process includes:

the experimental model building module 100 is configured to build a theoretical experimental model of a plate diving deformation structure according to geological conditions of a zone to be studied.

A simulation model building module 200, configured to build a discrete element numerical simulation model based on the theoretical experimental model;

and the plate diving analysis module 300 is configured to analyze deformation characteristics of the plate diving deformation based on the discrete element numerical simulation model.

Fig. 2 is a schematic structural diagram of an experiment model building module according to an embodiment of the present invention, please refer to fig. 2, the experiment model building module includes:

and a basic parameter determining module 101, configured to determine basic parameters based on the shells of land and ocean and for the occurrence of a dive effect.

In one embodiment of the invention, a theoretical experimental model is constructed based on the actual geological conditions of the zone to be studied. Specifically, basic parameters based on the theoretical experimental model are needed for establishing the theoretical experimental model related to plate diving deformation, so that the obtained theoretical experimental model can well reflect actual geological conditions and is used for researching actual geological changes.

The basic parameters of the theoretical experimental model specifically include: the thickness of the land-crust stratum, the thickness of the ocean-crust stratum, the mechanical parameters of the land-crust stratum and the ocean-crust stratum, the dive amount, the dive speed and the dive angle. Based on the basic parameters, a theoretical experimental model for analyzing actual geological change can be established, the parameters for establishing the theoretical experimental model are few, and the theoretical experimental model and the analysis processing process thereof can be simplified.

And the experimental model establishing module 102 is configured to establish a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions.

Specifically, referring to the experimental model and the actual geological conditions, the left side is a thick-layer land-shell stratum, and the middle inclined part and the right side horizontal zone are ocean-shell stratums with certain thicknesses. Correspondingly, the thickness of the land-shell stratum is greater than that of the ocean-shell stratum, the thickness of the land-shell stratum is twice that of the ocean-shell stratum, and the rigidity of the ocean-shell stratum is greater than that of the land-shell stratum, so that the characteristic that the ocean-shell stratum is not easy to deform is embodied. The diving amount and the diving speed of each group of experiments are the same, and only the diving angle is changed.

In an embodiment of the present invention, referring to fig. 9(a1) to 9(c2), the theoretical experimental model for constructing the plate diving deformation structure based on the crustal and oceanic shells and the basic parameters for the diving action specifically includes: and establishing the rock mechanical parameter setting of the differentiation of the land crust and the ocean crust stratum and the theoretical experimental model of the ocean crust diving downwards at different angles by the foundation parameters based on the land crust, the ocean crust and the diving action.

In the invention, the theoretical experimental model comprises the following setting principles:

the crustal stratum is referenced to a sandstone stratum with relatively low rigidity, the ocean shell is referenced to a sandstone stratum with relatively high rigidity and high deformation resistance, the thickness of the crustal stratum is larger than that of the ocean shell, and the ocean shell is provided with a movable base with a plurality of telescopic points and has the characteristic of diving in a specific direction.

Fig. 3 is a schematic structural diagram of a simulation model building module 200, which obtains mechanical parameters and an active base elongation direction that are theoretically required to be set for simulating a rock formation based on a land shell and an ocean shell and a fundamental parameter for generating a diving effect, so as to build a discrete element numerical simulation model, please refer to fig. 3, where the simulation model building module 200 includes:

and the boundary constructing module 201 is used for constructing the bottom and lateral boundaries of the dive plate and the land shell plate according to the theoretical experimental model.

Specifically, the mechanical property of stratum rock and the dive angle are kept unchanged, the length of the land shell base, the length of the ocean shell base and the height of the lateral boundary are reduced to a proper range according to a scale, and the length of the required simulated land shell base, the length of the ocean shell base, the length of the dive zone and the height of the lateral boundary are obtained. For example, down-sized on a scale of 1: 1000000.

And constructing the boundary of the simulation model based on four parameters of the length of the land-shell base, the length of the ocean-shell base, the height of the lateral boundary and the diving angle. The boundary of the constructed simulation model has a diving plate block and a lateral boundary which can move to a specific direction so as to be the same as the actual terrain structure.

Specifically, in order to adapt to structures of different simulated plate depression deformation, the depression angle and the direction of the depression belt can be freely defined, and the depression plate can extend to any position at a specific angle and can move with the boundary in the same horizontal displacement towards the left.

And the mechanical parameter setting module 202 is used for setting rock mechanical parameters, substrate and boundary mechanical parameters of the discrete element numerical simulation model. The friction coefficient, the elastic modulus, the normal direction and the tangential bonding strength of the simulated continental rise stratum and the ocean shell stratum respectively refer to the corresponding mechanical properties of the loose quartz sand, and small and large discontinuous values are taken in a certain effective interval, so that a simulated fault layer combination is formed on the simulated stratum.

Specifically, a theoretical experimental model established by the mechanical properties of the formation rock can obtain the mechanical parameters required to be set for simulating the rock formation, and accordingly, the mechanical parameter setting of the discrete element particles of the simulated formation should be set by referring to the mechanical parameters of physical experimental materials which are verified by a large number of experiments, namely the mechanical parameter setting of the simulated formation refers to the loose quartz sand of the physical experimental materials. Because the discrete element numerical simulation method applies mechanical properties between microscopic discrete element particles and between the microscopic discrete element particles, the mechanical properties are different from macroscopic mechanical parameters measured by an actual material under experimental conditions to a certain extent, but the difference is within a certain range, interval values are taken for the friction coefficient, the elastic modulus, the forward bonding strength and the tangential bonding strength of a simulated ocean-shell stratum and a simulated land-shell stratum within a certain effective range by taking the macroscopic experimental parameters as the reference until the mechanical parameters conforming to the fracture deformation characteristics of the actual material are obtained.

And the simulated stratum construction module 203 is used for constructing a simulated stratum which is adaptive to the geological condition required by the plate diving deformation based on the bottom and the lateral boundary. The simulated stratum is positioned in the bottom and lateral boundaries and is parallel to the constructed bottom boundary, the left land-shell stratum and the right ocean-shell stratum are horizontally spread from bottom to top, the middle ocean-shell stratum is obliquely spread, and the ground layers and the boundaries have contact relations which accord with the mechanical properties of rocks.

Specifically, the simulated formation comprises: left-hand crustal strata, middle-ocean-shelled diving strata, and right-hand ocean-shelled strata. The simulated formation is used to simulate a real formation in actual geological conditions. The simulation stratum is parallel to the constructed bottom boundary and is tightly filled in the model boundary, the left crustal stratum and the right crustal stratum are horizontally laid on the crustal and the crustal base, and the middle crustal stratum is obliquely laid in the boundary. In order to adapt to simulated stratums with different conditions, the thicknesses of different simulated stratums, the stratum rock mechanical properties in different simulated stratums and the contact relation in different simulated stratums can be freely defined as required.

And the simulated formation dividing module 204 is used for dividing the simulated formation into small layers marked by different colors in equal thickness.

Specifically, the simulated stratums are simulated by adopting discrete element particles corresponding to mechanical parameters of the loose quartz sand, the rigidity of the discrete element particles of the ocean shell stratums is high, after all the particles are filled, the particles are placed still until all the particles are in stable contact, the whole set of simulated stratums are divided by the same thickness, and different stratums are distinguished by different color marks. When the simulated stratum has structural change, the track of the structural change can be clearly reflected by the color mark.

And a hull plate dive module 205 for diving the hull plates at a fixed speed and at a specific angle.

Specifically, after the simulation model is built, a fixed speed is applied to the hull diving plate, so that the diving plate performs diving motion downwards in a specific direction, a speed which is the same as the horizontal component of the side boundary connected with the side boundary is applied to the side boundary, so that the side boundary performs uniform motion towards one side, and the position of the other side boundary is kept unchanged. Correspondingly, under the action of friction force, the ocean-shell stratum and the stratum within the control range of the side boundary are driven to move together, and the ocean-shell stratum drives the overlying land-shell stratum to move. With the increase of the extension, the stratum is fractured and deformed, and the land crust stratum under different diving angles develops different types of structural deformation.

Referring to fig. 9(a1) to 9(c2), in order to ensure that the depth of the ocean hulls at the beginning of diving is the same and the lengths of the ocean hulls and the land hulls on both sides of the diving zone are the same, the lengths of three sets of experimental models are set differently according to the diving angles of the ocean hulls, the experimental model set 1 is set to be 100m long, the experimental model set 2 is set to be 120m long, the experimental model set 3 is set to be 200m long, but the experimental models set in three sets are all 10m high, and are all set at 15% porosity, 2600kg/m³Random packing particle size of 0.125m ℃0.075m two-dimensional disk particles, the number of particles of different sizes obeying a gaussian distribution. The test model is characterized in that the left-side land-shell stratum is divided into 10 layers, the thickness of each layer is 1m, the middle and right-side ocean-shell strata are divided into 5 layers, and the thickness of each layer is 1m, and the layers are distinguished by different colors, so that the fault development condition can be observed conveniently. The right ocean shell substrate is provided with a movable wall body which dives towards the left side at different angles together with the right wall body, and the diving rate is 25%. Except for the difference of key research attributes, other parameters and boundary conditions of the experimental model are kept consistent so as to eliminate the influence of irrelevant factors. Through repeated experiments, the elastic modulus of the land-crust stratum of the experimental main body is 2 multiplied by 10⁷pa, elastic modulus of the Shell formation of 4X 10⁸pa, the characteristic of large rigidity and difficult deformation of the ocean-hull stratum is obviously shown. The ground layer inter-particle property is set to be 0.6, so that the simulated material forms a series of fault layer combination forms meeting the structural recognition on the basis of complying with the Coulomb shear fracture rule.

Fig. 4 is a schematic structural diagram of a plate dive analysis module 300, which analyzes deformation characteristics of a plate dive deformation process based on a simulation result of the discrete element numerical simulation model, please refer to fig. 4, where the plate dive analysis module 300 includes:

the simulation data determining module 301 is configured to form different simulation fault patterns by changing simulation process parameters of the discrete element numerical simulation model, and obtain corresponding simulation data.

The simulation process parameters are specifically as follows: the simulation process comprises a diving mode, a diving amount, a diving speed, a base and boundary mechanical parameter, and the thickness and the mechanical parameter of a land crust stratum and an ocean crust stratum.

And a deformation characteristic observation module 302, configured to observe deformation characteristics of the plate depression deformation process corresponding to different depression angles.

And the degeneration characteristic summarizing module 303 is used for summarizing the differential deformation characteristics of the plate diving deformation process at different angles based on the simulated formation and fault deformation development rules.

In a specific embodiment of the invention, the deformation characteristics of the plate depression deformation process when the depression angle is 30 degrees are observed, the deformation characteristics of the plate depression deformation process when the depression angle is 15 degrees are observed, the deformation characteristics of the plate depression deformation process when the depression angle is 5 degrees are observed, and the differential deformation characteristics of the plate depression deformation process at different angles are summarized based on the simulated formation and fault deformation development rules.

Taking the plate diving experiment model 1 as an example, referring to fig. 9(a1) and fig. 9(a2), experimental settings show that, in the initial stage of model construction, the left and right ocean-hull stratums are horizontally spread, the middle ocean-hull stratum is obliquely spread, the diving angle is 30 °, the boundaries of the right and middle ocean-hull stratums move uniformly from right to left at a horizontal speed of 0.5m/s, and the ocean-hull stratums are driven to perform a squeezing motion on the above-ground-hull stratums to generate corresponding deformation.

Taking the plate diving experimental model 1 as an example, referring to fig. 10(b), in an early stage of extrusion deformation when the diving amount reaches 5%, as the ocean-hull formation moves to the left, the overlying land-hull formation is driven to move to the left synchronously, the upper land-hull formation is extruded to rapidly generate fracture deformation, recoil faults F1, F2, F3 and F4 inclined to the land-hull are formed successively, and as the diving amount increases, faults F5 and F6 inclined to the ocean-hull are formed gradually. The development of the inclined crust fault is earlier than that of the recoil fault, the development mode is forward-extending, the forward-extending characteristic is not obvious because the crust stratum on the right side limits the continuous development of the fault, the recoil fault is developed later, the development mode is forward-extending, the development is gradually carried out towards the crust, and the small branch fault F7 is developed on the forward-extending type, the fault distance is small, and the function of adjusting the extrusion amount of the crust stratum is mainly played.

Taking the plate diving experimental model 1 as an example, referring to fig. 10(c), in the early and medium stage of extrusion deformation when the diving amount reaches 10%, the number of recoil faults is increased, the faults expand from a forward expansion mode to a land crust direction to form a fault F8, and are in a near-parallel relation with the previously formed recoil faults, and then small branch faults are formed; the number of inclined ocean shells has no obvious change, and the inclined ocean shells downwards dive to have strong activity influence, so that the inclined ocean shells extend to the deep part, and the fault distance is obviously increased; overall it forms a nose-down wedge with increased width and height.

Taking the plate diving experimental model 1 as an example, referring to fig. 10(d), in the middle stage of extrusion deformation when the diving amount reaches 15%, the recoil fault continues to develop, and the fault expands towards the direction of the land crust to form a fault F10 which is in a near-parallel relationship with the previously formed recoil fault; the number of inclined ocean hull faults has no obvious change, and the fault distance is continuously increased; its overall form nose-down wedge width and height continue to increase.

Taking the plate diving experimental model 1 as an example, referring to fig. 10(e), in a middle and later stage extrusion deformation stage when the diving amount reaches 20%, the number of faults is obviously increased, the recoil fault continues to develop, the faults F10, F11 and F12 are formed successively, and the three faults control the lifting of the stratum on the left side of the diving wedge; the number of inclined ocean hull faults has no obvious change, and the fault distance is continuously increased; the total width of the dive wedge is increased, the stratum on the left side continues to be lifted, the form of the stratum is changed, and the left side is changed into a wedge with a slightly higher left side than right side from the previous left low and right high.

Taking the plate diving experimental model 1 as an example, as shown in fig. 10(F), in a late stage extrusion deformation stage when the diving amount reaches 25%, the number of faults has no obvious change, the breaking distances of the recoil faults F10, F11 and F12 are increased, and the three faults control the continuous lifting of the stratum on the left side of the diving wedge; the number of inclined ocean hull faults has no obvious change, and the fault distance is continuously increased; the total width of the dive wedge is increased, the left stratum is continuously lifted, and the left side is transformed into a wedge with a height larger than the right side.

Taking the plate diving experimental model 1 as an example, referring to fig. 10(a) to 10(F) showing the plate diving deformation process, the diving angle of the ocean hull is 30 degrees, the diving amount is gradually increased from 0 percent to 25 percent, overall, at the contact part of the land hull and the ocean hull, the land hull fault is gradually inclined to develop in a forward-expanding mode, the fault distance is continuously increased, then the recoil fault develops in a forward-expanding mode and develops towards the interior of the land hull, the small branch fault develops, when the diving amount reaches 13 percent, the first fault F10 controlling the whole form of the diving wedge is formed, the diving amount is gradually increased, when the diving amount reaches 16 percent, the second fault F11 controlling the whole form of the diving wedge is formed, the diving amount is continuously increased, when the whole diving amount reaches 17 percent, the third fault F12 controlling the whole form of the diving wedge is formed, the diving height is continuously increased due to the diving effect, the width is continuously increased, and the shape of the wedge gradually changes from low left to high right to a wedge shape with the left side higher than the right side.

Taking the plate diving experimental model 2 as an example, referring to fig. 9(b1) to 9(b2), in the initial stage of model construction, the experiment sets the diving angle of the ocean hull to be 15 °, the total length of the land hull to be 80m, the total length of the ocean hull to be 60m, the other conditions and parameter settings are the same as those of the experimental model 1, and the boundaries of the right-side and middle ocean hull strata move uniformly from right to left at a horizontal speed of 0.5 m/s.

Taking the plate diving experimental model 2 as an example, referring to the plate diving deformation process shown in fig. 11(a) to 11(f), the hull diving angle is 15 °, and the diving amount is gradually increased from 0% to 25%, and the deformation process is similar to the deformation process under the 30 ° diving angle, but the deformation characteristics are slightly different. At the contact part of the land crust and the ocean crust, the land crust fault is developed in sequence in a forward-unfolding mode, the fault distance is continuously increased, the recoil fault is developed in a sequential-unfolding mode, the fault gradually develops towards the inside of the land crust, the fault is developed in a small branch fault, when the diving amount reaches 17%, a first fault F14 for controlling the integral form of the diving wedge is formed, the diving amount is gradually increased, when the diving amount reaches 22%, a second fault F15 for controlling the integral form of the diving wedge is formed, the diving amount is continuously increased, when the diving amount reaches 25%, a third fault F16 for controlling the integral form of the diving wedge is formed, the width of the diving wedge formed due to the diving effect is continuously increased, the height is continuously increased, and the form of the diving wedge is gradually changed from a left lower part to a right part to a left part and a right part into a more regular wedge with the height being approximately equal to the left side and the.

Taking the plate diving experimental model 3 as an example, referring to fig. 9(c1) and fig. 9(c2), in the initial stage of model construction, the experiment sets the diving angle of the ocean hull to be 5 °, the total length of the land hull to be 160m, the total length of the ocean hull to be 100m, the other conditions and parameter settings are the same as those of the experimental model 1 and the experimental model 2, and the boundaries of the right-side and middle ocean hull strata move uniformly from right to left at a horizontal speed of 0.5 m/s.

Taking the plate diving experimental model 3 as an example, referring to fig. 12(b), in the early stage of the crush deformation when the diving amount reaches 5%, the F1, F2, F3, F4 and F5 inclined to the crust are formed successively, and as the diving amount increases, the small inclined crusty branch faults F7 and F6 develop on the F1 and F3. Because of the small angle of pitch, the impact on the crust formation is relatively uniform, so when the amount of pitch is small, the F1, F2, F3, F4, and F5 faults tending towards the crust develop almost simultaneously, and then as the pitch continues, branch faults tending towards the crust develop in a small scale.

Taking the plate diving experimental model 3 as an example, referring to fig. 12(c), in the middle and early stage extrusion deformation stage when the diving amount reaches 10%, the ocean shells continue to dive downwards, the number of faults is obviously increased, and parallel faults which are close to the ocean shells are developed under the faults F1, F2, F3, F4 and F5 which tend to the land shells, and the closely-spaced parallel faults have similar properties, so that the stratum fault distances controlled by the closely-spaced parallel faults are obviously increased by using the same codes F1, F2, F3, F4 and F5; recoil faults are increased greatly, and development faults are F8, F9, F10, F11 and F12; overall it forms a nose-down wedge with an increased height.

Taking the plate diving experimental model 3 as an example, referring to fig. 12(d), in the middle stage of extrusion deformation when the diving amount reaches 15%, as the diving action continues, the number of faults increases, and the formation fault distances controlled by F1, F2, F3, F4 and F5 continue to increase; continuing to develop the inclined hull branch fault F14, and continuing to increase the recoil fault distance; the height of the formed diving wedge is obviously increased on the whole.

Taking the plate diving experimental model 3 as an example, referring to fig. 12(e), in the middle and later stage extrusion deformation stage when the diving amount reaches 20%, the number of recoil faults is significantly increased, and faults F15 and F16 are formed successively, and the two faults control the lifting of the stratum on the left side of the diving wedge; the number of inclined crust faults has no obvious change, and the fault distance is continuously increased; the formation on the left side of the nose-down wedge continues to lift as a whole.

Taking the plate diving experimental model 3 as an example, referring to fig. 12(F), in a later stage of extrusion deformation when the diving amount reaches 25%, the number of faults does not change obviously, the breaking distances of the recoil faults F15 and F16 are increased obviously, and the two faults control the continuous lifting of the stratum on the left side of the diving wedge; the number of inclined crust faults has no obvious change, the fault distance is continuously increased, and the fault dip angle is steep; the interval diminishes between the different faults on the whole, arranges inseparabler, and dive wedge left side stratum continues the lifting, and the slope is close to the wedge that equals about evolving.

Taking the plate diving experimental model 3 as an example, referring to the plate diving deformation process shown in fig. 12(a) to 12(F), the diving angle of the ocean shell is 5 °, the diving amount is gradually increased from 0% to 25%, overall, firstly, the recoil fault is developed in the land shell stratum, the fault distance is continuously increased, then, the successive development tends to the land shell fault, the fault is gradually developed towards the interior of the land shell, the small branch fault is developed, when the diving amount reaches 20%, the first fault F15 for controlling the whole form of the diving wedge is formed, the diving amount is gradually increased, when the diving amount reaches 23%, the second fault F16 for controlling the whole form of the diving wedge is formed, the height of the diving wedge formed due to the diving effect is continuously increased, and the form of the diving wedge gradually changes from low left to high right to wedge with the left and right slopes nearly equal.

By comparing the structural diagrams formed by the plate diving states in the final three groups of experiments, it can be seen that after the plate diving amount is 25%, the plate diving deformation fracture system is developed in the experiments 1, 2 and 3. The following rules can be summarized for the comparison of three sets of experiments: with the increase of the dive amount, the number of faults is continuously increased, and the integral fault distance is continuously increased; under different diving angles, the lower the diving angle is, the more faults are formed; the lower the dive angle is, the later the kick fault forming time of the integral form of the dive wedge is controlled; the lower the dive angle, the higher the height and the larger the width of the dive wedge, and the different forms of the dive wedge.

Particularly, the geological phenomenon of the depression deformation process of the land shell plate is well reproduced according to the experimental result, and the depression deformation characteristic of the plate is associated with the depression angle. The simulation mode can be used for analyzing the influence of the plate diving angle on the deformation characteristics, so that the understanding of the terrain is further enhanced, and a theoretical guidance basis is provided.

In one embodiment of the invention, the thickness of the land-shell formation is greater than the thickness of the ocean-shell formation to ensure that the nose-down phenomenon is realistic; the ocean-shell stratum and the ocean-shell stratum mechanical parameters are set with reference to the loose quartz sand of the physical experiment material, the ocean-shell stratum rigidity is set to be larger than the land-shell stratum rigidity, and the parameters are properly adjusted based on the discrete element principle.

In an embodiment of the present invention, analyzing deformation characteristics of deformation of a plate during a dive deformation process based on a simulation result of the discrete element numerical simulation model specifically includes:

filling high-density discrete element particles for simulating a left land crust stratum, a middle ocean crust diving zone stratum and a right ocean crust stratum in sequence in the boundary of the model, wherein the whole set of simulated stratum is divided into small layers marked with different colors in equal thickness;

the friction coefficient, the elastic modulus, the normal direction and the tangential bonding strength of the simulated continental rise stratum and the ocean shell stratum respectively refer to the corresponding mechanical properties of the loose quartz sand, and small and large discontinuity values are taken in a certain effective interval so as to form a simulated fault layer combination on the simulated stratum;

changing the simulated dive zone dive angle, and continuously taking discontinuous values for the mechanical properties of the simulated formation until plate dive deformation is formed;

and acquiring the corresponding relation between the simulation process parameters and the plate diving deformation process based on the deformation development rule of the simulated stratum and the fault.

The thickness ratio of the simulated land crust stratum to the ocean crust stratum is controlled to be 2: 1 or so, and the elastic modulus of the simulated land crust stratum is 2 multiplied by 10⁷pa, coefficient of friction 0.6, normal and tangential bond strengths all 3.6 × 10⁵pa, the modulus of elasticity of the simulated crustacean layer is 4 x 10⁸pa, coefficient of friction 0.6, normal and tangential bond strengths all 3.6 × 10⁵pa。

The discrete element simulation system for the plate diving deformation process, provided by the invention, is characterized in that an experimental model is constructed according to actual geological condition parameters, and a discrete element numerical simulation model is constructed by combining theoretical parameters obtained by experimental calculation of the experimental model; the discrete element numerical simulation model simulates and reproduces an actual geological structure, and the influence of the current actual geological condition on the fault development pattern is analyzed according to the corresponding relation between the parameter change and the simulation result in the simulation process, so that a theoretical basis is provided for geological research.

Furthermore, although in the above detailed description several unit modules of the system are mentioned, this division is not mandatory only. Indeed, the features and functions of two or more of the units described above may be embodied in one unit, according to embodiments of the invention. Also, the features and functions of one unit described above may be further divided into embodiments by a plurality of units. The terms "module" and "unit" used above may be software and/or hardware that realizes a predetermined function. While the modules described in the following embodiments are preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated.

Having described the discrete element simulation system of the plate dive deformation process of an exemplary embodiment of the present invention, a method of an exemplary embodiment of the present invention will now be described with reference to the accompanying drawings. The implementation of the method can be referred to the above overall implementation, and repeated details are not repeated.

The invention provides a discrete element simulation method for a plate diving deformation process, which aims to solve the technical problems that the cause research and analysis process of a geological structure cannot be verified, and an effective association means is lacked between a theory and an actual cause.

Specifically, fig. 5 is a schematic flow chart of a discrete element simulation method for a plate dive deformation process according to an embodiment of the present invention, please refer to fig. 5, where the discrete element simulation method for the plate dive deformation process includes:

s101: and constructing a theoretical experimental model of the plate diving deformation structure according to the geological conditions of the zone to be researched.

S102: establishing a discrete element numerical simulation model based on the theoretical experimental model;

s103: and analyzing the deformation characteristics of the plate diving deformation based on the discrete element numerical simulation model.

Fig. 6 is a detailed flowchart of step S101, please refer to fig. 6, which includes the following steps:

s201: the fundamental parameters based on the shells of land and ocean and the occurrence of a dive effect are determined.

S202: and establishing a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions.

Fig. 7 is a specific flowchart of establishing a discrete element numerical simulation model based on the theoretical experimental model, where the mechanical parameters and the moving base elongation direction that need to be set for theoretically simulating a rock formation are obtained based on the land and ocean shells and the fundamental parameters of the occurrence of the dive action to construct the discrete element numerical simulation model, please refer to fig. 7, and establishing the discrete element numerical simulation model based on the theoretical experimental model includes:

s301: and constructing the bottom and lateral boundaries of the diving plate and the land shell plate according to the theoretical experimental model.

S302: and setting the rock mechanical parameters and the substrate and boundary mechanical parameters of the discrete element numerical simulation model. The friction coefficient, the elastic modulus, the normal direction and the tangential bonding strength of the simulated continental rise stratum and the ocean shell stratum respectively refer to the corresponding mechanical properties of the loose quartz sand, and small and large discontinuous values are taken in a certain effective interval, so that a simulated fault layer combination is formed on the simulated stratum.

S303: and constructing a simulated stratum adaptive to the geological condition required by the plate diving deformation based on the bottom and the lateral boundary. The simulated stratum is positioned in the bottom and lateral boundaries and is parallel to the constructed bottom boundary, the left land-shell stratum and the right ocean-shell stratum are horizontally spread from bottom to top, the middle ocean-shell stratum is obliquely spread, and the ground layers and the boundaries have contact relations which accord with the mechanical properties of rocks.

S304: and dividing the simulated stratum into small layers marked by different colors in equal thickness.

S305: so that the ocean shell plate generates a diving effect towards the land shell at a fixed speed and a specific angle.

Referring to fig. 9(a1) to 9(c2), in order to ensure that the depth of the ocean hulls at the beginning of diving is the same and the lengths of the ocean hulls and the land hulls on both sides of the diving zone are the same, the lengths of three sets of experimental models are set differently according to the diving angles of the ocean hulls, the experimental model set 1 is set to be 100m long, the experimental model set 2 is set to be 120m long, the experimental model set 3 is set to be 200m long, but the experimental models set in three sets are all 10m high, and are all set at 15% porosity, 2600kg/m³The density is randomly filled with two-dimensional disc particles with the particle size of 0.125-0.075 m, and the number of the particles with different particle sizes is subjected to Gaussian distribution. The test model is characterized in that the left-side land-shell stratum is divided into 10 layers, the thickness of each layer is 1m, the middle and right-side ocean-shell strata are divided into 5 layers, and the thickness of each layer is 1m, and the layers are distinguished by different colors, so that the fault development condition can be observed conveniently. Right side ocean shell baseThe bottom of the movable wall body is provided with a movable wall body which dives towards the left side at different angles together with the right side wall body, and the diving rate is 25%. Except for the difference of key research attributes, other parameters and boundary conditions of the experimental model are kept consistent so as to eliminate the influence of irrelevant factors. Through repeated experiments, the elastic modulus of the land-crust stratum of the experimental main body is 2 multiplied by 10⁷pa, elastic modulus of the Shell formation of 4X 10⁸pa, the characteristic of large rigidity and difficult deformation of the ocean-hull stratum is obviously shown. The ground layer inter-particle property is set to be 0.6, so that the simulated material forms a series of fault layer combination forms meeting the structural recognition on the basis of complying with the Coulomb shear fracture rule.

Fig. 8 is a schematic flow chart illustrating a process of analyzing deformation characteristics of the plate dive deformation based on the discrete element numerical simulation model, please refer to fig. 8, where analyzing deformation characteristics of the plate dive deformation based on the discrete element numerical simulation model includes:

s401: and forming different simulation fault patterns by changing simulation process parameters of the discrete element numerical simulation model, and obtaining corresponding simulation data.

S402: and observing deformation characteristics of the plate in the depression deformation process corresponding to different depression angles.

S403: and summarizing the differential deformation characteristics of the sheet diving deformation process under different angles based on the simulated formation and fault deformation development rules.

The present invention also provides a computer device comprising: the system comprises a processor and a storage device, wherein the storage device is used for storing a plurality of instructions, and the instructions are suitable for being loaded by the processor and executing a discrete element simulation method of a slab diving deformation process.

The invention also provides a computer-readable storage medium, in which a computer program is stored, the computer program being configured to perform a discrete element simulation method of a sheet dive deformation process.

The beneficial effects of the invention are mainly as follows:

(1) the experimental model is simple in construction mode, and based on the discrete element numerical simulation model constructed by the experimental model, the simulation result can well reproduce the structure of an actual geological fault so as to verify the corresponding relation between the parameter change and the formed fault structure, and the method can be used for analyzing and predicting the influence of the actual current geological condition on the fault development mode;

(2) parameters related to the construction of the experimental model are few, and the calculation process is simple;

(3) the construction mode of the discrete element numerical simulation model is simple and easy to realize, compared with a physical simulation model, the discrete element numerical simulation model is not limited by experimental materials and model design, and the experimental model can be freely defined according to requirements;

(4) during the simulation process and in the final state, the model can be stored at any time, and analytical data such as speed, displacement and the like can be extracted and sliced according to any position;

(5) in the actual simulation process, the change of the simulation condition can be conveniently controlled according to the actual geological condition, the implementation process is simple and easy to operate, and the reproducibility of the actual fault is good.

Improvements to a technology can clearly be distinguished between hardware improvements (e.g. improvements to the circuit structure of diodes, transistors, switches, etc.) and software improvements (improvements to the process flow). However, as technology advances, many of today's process flow improvements have been seen as direct improvements in hardware circuit architecture. Designers almost always obtain the corresponding hardware circuit structure by programming an improved method flow into the hardware circuit. Thus, it cannot be said that an improvement in the process flow cannot be realized by hardware physical modules. For example, a Programmable Logic Device (PLD), such as a Field Programmable Gate Array (FPGA), is an integrated circuit whose Logic functions are determined by programming the Device by a user. A digital system is "integrated" on a PLD by the designer's own programming without requiring the chip manufacturer to design and fabricate application-specific integrated circuit chips. Furthermore, nowadays, instead of manually making an Integrated Circuit chip, such Programming is often implemented by "logic compiler" software, which is similar to a software compiler used in program development and writing, but the original code before compiling is also written by a specific Programming Language, which is called Hardware Description Language (HDL), and HDL is not only one but many, such as abel (advanced Boolean Expression Language), ahdl (alternate Language Description Language), traffic, pl (core unified Programming Language), HDCal, JHDL (Java Hardware Description Language), langue, Lola, HDL, laspam, hardbyscript Description Language (vhr Description Language), and the like, which are currently used by Hardware compiler-software (Hardware Description Language-software). It will also be apparent to those skilled in the art that hardware circuitry that implements the logical method flows can be readily obtained by merely slightly programming the method flows into an integrated circuit using the hardware description languages described above.

The controller may be implemented in any suitable manner, for example, the controller may take the form of, for example, a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro) processor, logic gates, switches, an Application Specific Integrated Circuit (ASIC), a programmable logic controller, and an embedded microcontroller, examples of which include, but are not limited to, the following microcontrollers: the ARC625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320, the memory controller may also be implemented as part of the control logic for the memory.

Those skilled in the art will also appreciate that, in addition to implementing the controller as pure computer readable program code, the same functionality can be implemented by logically programming method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Such a controller may thus be considered a hardware component, and the means included therein for performing the various functions may also be considered as a structure within the hardware component. Or even means for performing the functions may be regarded as being both a software module for performing the method and a structure within a hardware component.

The systems, devices, modules or units illustrated in the above embodiments may be implemented by a computer chip or an entity, or by a product with certain functions.

For convenience of description, the above devices are described as being divided into various units by function, and are described separately. Of course, the functionality of the units may be implemented in one or more software and/or hardware when implementing the present application.

From the above description of the embodiments, it is clear to those skilled in the art that the present application can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present application may be essentially or partially implemented in the form of software products, which may be stored in a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and include instructions for causing a computer system (which may be a personal computer, a server, or a network system, etc.) to execute the methods described in the embodiments or some parts of the embodiments of the present application.

The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the system embodiment, since it is substantially similar to the method embodiment, the description is simple, and for the relevant points, reference may be made to the partial description of the method embodiment.

The application is operational with numerous general purpose or special purpose computing system environments or configurations. For example: personal computers, server computers, hand-held or portable systems, tablet-type systems, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics systems, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or systems, and the like.

The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing systems that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage systems.

While the present application has been described with examples, those of ordinary skill in the art will appreciate that there are numerous variations and permutations of the present application without departing from the spirit of the application, and it is intended that the appended claims encompass such variations and permutations without departing from the spirit of the application.

Claims

1. A method for discrete element simulation of a sheet nose-down deformation process, the method comprising:

constructing a theoretical experimental model of a plate diving deformation structure according to the geological conditions of the zone to be researched; wherein, according to the geological conditions of the zone needing to be researched, the theoretical experimental model for constructing the plate diving deformation structure comprises the following steps: determining basic parameters based on the shells of land and ocean and the occurrence of a diving effect; establishing a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions; the setting principle comprises that the land crust takes a sandstone stratum with rigidity smaller than a first threshold value as a reference, the ocean crust takes a sandstone stratum with deformation resistance larger than a second threshold value and rigidity larger than a third threshold value as a reference, the ocean crust is provided with a movable base with a plurality of telescopic points and has the characteristic of diving in any direction;

2. The method of claim 1, wherein the base parameters include a thickness of the crust formation, mechanical parameters of the crust formation and the crust formation, a nose-down amount, a nose-down speed, and a nose-down angle.

3. The method of claim 1, wherein building a discrete element numerical simulation model based on the theoretical experimental model comprises:

4. The method of claim 3, wherein analyzing deformation characteristics of the plate dive deformation based on the discrete-element numerical simulation model comprises:

5. A discrete element simulation system for a sheet dive deformation process, the system comprising:

the plate depression analysis module is used for analyzing deformation characteristics of plate depression deformation based on the discrete element numerical simulation model;

wherein, the experimental model construction module comprises:

the experimental model establishing module is used for establishing a theoretical experimental model of plate diving deformation based on a set principle according to the basic parameters and the geological conditions; the setting principle comprises that the land crust takes a sandstone stratum with rigidity smaller than a first threshold value as a reference, the ocean crust takes a sandstone stratum with deformation resistance larger than a second threshold value and rigidity larger than a third threshold value as a reference, the ocean crust is provided with a movable base with a plurality of telescopic points and has the characteristic of diving in any direction.

6. The system of claim 5, wherein the base parameters include a thickness of the crust formation, mechanical parameters of the crust formation and the crust formation, a nose-down amount, a nose-down speed, and a nose-down angle; the setting principle comprises that the land crust takes a sandstone stratum with rigidity smaller than a first threshold value as a reference, the ocean crust takes a sandstone stratum with deformation resistance larger than a second threshold value and rigidity larger than a third threshold value as a reference, the ocean crust is provided with a movable base with a plurality of telescopic points and has the characteristic of diving in any direction.

7. The system of claim 5, wherein the simulation model building module comprises:

8. The system of claim 7, wherein the slab dive analysis module comprises:

9. A computer device, comprising: a processor adapted to implement instructions and a storage device storing instructions adapted to be loaded by the processor and to perform a discrete element simulation method of a sheet dive deformation process according to any of claims 1 to 4.

10. A computer-readable storage medium, characterized in that a computer program is stored for performing a discrete element simulation method of a sheet dive deformation process according to any one of claims 1 to 4.