CN110032788A

CN110032788A - A kind of discrete element simulation method, equipment and the system of plate subduction deformation process

Info

Publication number: CN110032788A
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: 2019-07-19
Anticipated expiration: 2039-04-04
Also published as: CN110032788B

Abstract

The present invention provides discrete element simulation method, system, computer equipment and the computer readable storage mediums of a kind of plate subduction deformation process, are related to technical field of geological exploration.This method includes the geological conditions in the area studied as needed, constructs the theoretical experimental model of plate subduction deformation structure；Discrete element numerical simulation model is established based on the theoretical experimental model；Deformation behaviour based on the deformation of plate subduction described in the discrete element numerical simulation model analysis.The technical issues of model construction mode of the present invention is simple, and analog form is easily achieved, and analog result accuracy is good, and solving the study on the genesis analytic process of geologic structure cannot verify, lack efficient association means between theory and the practical origin cause of formation.

Description

A kind of discrete element simulation method, equipment and the system of plate subduction deformation process

Technical field

The present invention is about technical field of geological exploration, especially with regard to the analogue technique of data in geological prospecting, specifically Say be a kind of plate subduction deformation process discrete element simulation method, system, computer equipment and computer-readable storage medium Matter.

Background technique

When two lithospheric plates mutually collide, and one of plate drops to another plate by the process dived Under when be formed subduction zone.In general, subducting plate refers to the oceanic plate being made of oceanic crust because oceanic crust by Simatic substance is constituted, and density is larger, is easier to for continental crust.It can lead to the deformation of upper disk by subduction, Form underriding wedge.

Due to different geological conditions and variable condition, the development pattern of underriding wedge interrupting layer will affect.Currently, logical Cross and establish model to simulate practical geological conditions and its variation, be the main means that are used in the formation mechanism study of tomography it One.And Discrete-parcel method carries out dynamic pine using central difference method as a kind of method for numerical simulation based on molecular dynamics It relaxes and solves, calculating is simple and efficient, is that can be based on from the tool of thin microcosmic angle research mechanics feature and motion feature The characteristics of discrete particle free movement, makes it have special advantage in Discontinuous transmission splintering problem, large deformation problem field, It is widely used in tectonic geology research.

Analog study to plate subduction deformation process is still mainly reflected at present in the stratum to different mechanical attributes The stage of the analysis of mature fault situation.For example, Hardy et al. is with Discrete-parcel method to the inverse hair for rushing wedge two-way in upper plate It educates and is studied；Jonathan et al. has studied the development of accretionary wedge in subduction zone with finite element method.

Experiment mould but at present in the modeling and analytic process to plate subduction deformation process, about deformation of diving The research of type, evolutionary process is still not perfect, about plate subduction deformation process deformation behaviour need deeper into research.And And the building process of experimental model is usually more complex, it is more to be related to parameter, and theory analysis is closed due to lacking verifying appropriate System, so that the matching degree of theoretical analysis result and actual conditions is difficult to accurately be judged.

Therefore, a kind of new scheme how is provided, to realize the experimental model, evolutionary process deformed that more accurately dive Research is urgent technical problem to be solved in the field.

Summary of the invention

In view of this, the embodiment of the invention provides a kind of discrete element simulation method of plate subduction deformation process, system, Computer equipment and computer readable storage medium, model construction mode is simple, and analog form is easily achieved, and analog result is quasi- True property is good, and the study on the genesis analytic process for solving geologic structure cannot be verified, lack between theory and the practical origin cause of formation and effectively close Jointly the technical issues of section.

It is an object of the invention to provide a kind of discrete element simulation methods of plate subduction deformation process, comprising:

The geological conditions in the area studied as needed constructs the theoretical experimental model of plate subduction deformation structure；

Discrete element numerical simulation model is established based on the theoretical experimental model；

Deformation behaviour based on the deformation of plate subduction described in the discrete element numerical simulation model analysis.

Preferably, the geological conditions in the area studied as needed, the theoretical experiment mould of building plate subduction deformation structure Type includes:

Determine based on continental crust and oceanic crust and occur the underlying parameter of subduction；

The theory for establishing plate subduction deformation based on setting principle according to the underlying parameter and the geological conditions is real Test model.

Preferably, the underlying parameter includes continental crust formation thickness, oceanic crust formation thickness, continental crust stratum and oceanic crust stratum Mechanics parameter, underriding amount, diving speed and underriding angle；The setting principle include the continental crust with rigidity less than the first threshold The sandstone formation of value is reference, and the oceanic crust is greater than the sandstone that second threshold and rigidity are greater than third threshold value with non-deformability Layer is reference, and the oceanic crust has the moveable substrate of the scalable point in many places, and has the property dived to any direction.

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

According to the theoretical experimental model, bottom and the lateral boundaries of subducting plate and continental crust plate are constructed；

The rock mechanics parameters and substrate and boundary mechanics parameter of discrete element numerical simulation model are set；

Based on the bottom and lateral boundaries, the adaptable simulation of geological conditions needed for building is deformed with the plate subduction Stratum；

The simulated formation is divided into the substratum of different colours label with equal thickness；

To continental crust subduction occurs for oceanic crust plate with fixed speed and special angle.

Preferably, based on plate subduction described in the discrete element numerical simulation model analysis deformation deformation behaviour include:

By changing the simulated process parameter of the discrete element numerical simulation model, to form different simulation tomography samples Formula, and obtain corresponding analogue data；

Observe the deformation behaviour of the corresponding plate subduction deformation process of different dive angles；

Law of development is deformed based on the simulated formation and tomography, is summarised in plate subduction deformation process under different angle Differentiation deformation behaviour.

It is an object of the invention to provide a kind of discrete element simulation systems of plate subduction deformation process, comprising:

Experimental model constructs module, the geological conditions in the area for studying as needed, and building plate subduction deforms structure The theoretical experimental model made；

Simulation model constructs module, for establishing discrete element numerical simulation model based on the theoretical experimental model；

Plate subduction analysis module, for what is deformed based on plate subduction described in the discrete element numerical simulation model analysis Deformation behaviour.

Preferably, the experimental model building module includes:

Underlying parameter determining module, for determining based on continental crust and oceanic crust and occurring the underlying parameter of subduction；

Experimental model establishes module, for being established according to the underlying parameter and the geological conditions based on setting principle The theoretical experimental model of plate subduction deformation.

Preferably, the simulation model building module includes:

Boundary constructs module, for according to the theoretical experimental model, construct subducting plate and continental crust plate bottom and Lateral boundaries；

Mechanics parameter setup module, for rock mechanics parameters and substrate and the boundary of discrete element numerical simulation model to be arranged Mechanics parameter；

Simulated formation constructs module, and for being based on the bottom and lateral boundaries, building deforms institute with the plate subduction The simulated formation for needing geological conditions adaptable；

Simulated formation division module, for the simulated formation to be divided into the substratum of different colours label with equal thickness；

Oceanic crust plate subduction module is made for making oceanic crust plate occur to dive to continental crust with fixed speed and special angle With.

Preferably, the plate subduction analysis module includes:

Analogue data determining module, for passing through the simulated process parameter for changing the discrete element numerical simulation model, with Different simulation fault styles is formed, and obtains corresponding analogue data；

Vertic features Observation Blocks, for observing the deformation behaviour of the corresponding plate subduction deformation process of different dive angles；

Vertic features summarize module, for deforming law of development based on the simulated formation and tomography, are summarised in different angles Spend the differentiation deformation behaviour of lower plate subduction deformation process.

It is an object of the invention to provide a kind of computer equipments, comprising: be adapted for carrying out each instruction processor and Equipment is stored, the storage equipment is stored with a plurality of instruction, and described instruction is bowed suitable for being loaded by processor and being executed a kind of plate Rush the discrete element simulation method of deformation process.

It is an object of the invention to provide a kind of computer readable storage mediums, are stored with computer program, the meter Calculation machine program is used to execute a kind of discrete element simulation method of plate subduction deformation process.

The beneficial effects of the present invention are, provide the discrete element simulation method of plate subduction deformation process a kind of, system, Computer equipment and computer readable storage medium, model construction mode is simple, and analog form is easily achieved, and analog result is quasi- True property is good, and the study on the genesis analytic process for solving geologic structure cannot be verified, lack between theory and the practical origin cause of formation and effectively close Jointly the technical issues of section.

For above and other objects, features and advantages of the invention can be clearer and more comprehensible, preferred embodiment is cited below particularly, And cooperate institute's accompanying drawings, it is described in detail below.

Detailed description of the invention

In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.

Fig. 1 is a kind of structural representation of the discrete element simulation system of plate subduction deformation process provided in an embodiment of the present invention Figure；

Fig. 2 is experimental model in a kind of discrete element simulation system of plate subduction deformation process provided in an embodiment of the present invention Construct the structural schematic diagram of module；

Fig. 3 is simulation model in a kind of discrete element simulation system of plate subduction deformation process provided in an embodiment of the present invention Construct the structural schematic diagram of module；

Fig. 4 is that a kind of discrete element simulation system mesoplate of plate subduction deformation process provided in an embodiment of the present invention dives The structural schematic diagram of analysis module；

Fig. 5 is a kind of flow chart of the discrete element simulation method of plate subduction deformation process provided in an embodiment of the present invention；

Fig. 6 is the specific flow chart of the step S101 in Fig. 5；

Fig. 7 is the specific flow chart of the step S102 in Fig. 5；

Fig. 8 is the specific flow chart of the step S103 in Fig. 5；

Fig. 9 (a1) to Fig. 9 (a2) is that the plate subduction for being 30 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The simulator schematic diagram of the discrete element simulation method of process；

Fig. 9 (b1) to Fig. 9 (b2) is that the plate subduction for being 15 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The simulator schematic diagram of the discrete element simulation method of process；

Fig. 9 (c1) to Fig. 9 (c2) is that the plate subduction for being 5 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The simulator schematic diagram of the discrete element simulation method of process；

Figure 10 (a) to Figure 10 (f) is that the plate subduction for being 30 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The section evolution diagram of the analog result of the discrete element simulation method of process；

Figure 11 (a) to Figure 11 (f) is that the plate subduction for being 15 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The section evolution diagram of the analog result of the discrete element simulation method of process；

Figure 12 (a) to Figure 12 (f) is that the plate subduction for being 5 ° according to subduction of oceanic crust a kind of in embodiment of the present invention angle deforms The section evolution diagram of the analog result of the discrete element simulation method of process.

Specific embodiment

Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.

One skilled in the art will appreciate that embodiments of the present invention can be implemented as a kind of system, device, method or calculating Machine program product.Therefore, disclose can be with specific implementation is as follows by the present invention, it may be assumed that complete hardware, complete software (packet Include firmware, resident software, microcode etc.) or hardware and software combine form.

Below with reference to several representative embodiments of the invention, the principle and spirit of the present invention are explained in detail.

The present invention provides a kind of discrete element simulation system of plate subduction deformation process, is ground with solving the origin cause of formation of geologic structure The technical issues of studying carefully analytic process cannot verify, lack efficient association means between theory and the practical origin cause of formation.

Specifically, Fig. 1 is a kind of discrete element simulation system of plate subduction deformation process provided in an embodiment of the present invention Structural schematic diagram, referring to Figure 1, the discrete element simulation system of the plate subduction deformation process includes:

Experimental model constructs module 100, the geological conditions in the area for studying as needed, building plate subduction deformation The theoretical experimental model of construction.

Simulation model constructs module 200, for establishing discrete element numerical simulation model based on the theoretical experimental model；

Plate subduction analysis module 300, for being become based on plate subduction described in the discrete element numerical simulation model analysis The deformation behaviour of shape.

Fig. 2 is that the experimental model of the embodiment of the present invention constructs the structural schematic diagram of module, referring to Fig. 2, experimental model structure Modeling block includes:

Underlying parameter determining module 101, for determining based on continental crust and oceanic crust and occurring the underlying parameter of subduction.

In an embodiment of the invention, the practical geological conditions in the area studied as needed constructs theoretical real Test model.Specifically, it establishes the theoretical experimental model deformed about plate subduction and needs the basis ginseng based on theoretical experimental model Number makes obtained theoretical experimental model can be good at reflecting practical geological conditions, the research for moving to practically qualitative change.

The underlying parameter of theoretical experimental model specifically includes: continental crust formation thickness, oceanic crust formation thickness, continental crust stratum and ocean The mechanics parameter on shell stratum, underriding amount, diving speed, underriding angle.Based on these underlying parameters, can establish for analyzing The practically dynamic theoretical experimental model of qualitative change, the parameter that the experimental model that theorizes is related to is few, can simplify theoretical experimental model And it analyzes treatment process.

Experimental model establishes module 102, for being based on setting principle according to the underlying parameter and the geological conditions Establish the theoretical experimental model of plate subduction deformation.

Specifically, referring to experimental model and practical geological conditions, left side for thick-layer continental crust stratum, intermediate, inclined position and It is certain thickness oceanic crust stratum in the horizontal stripes of right side.Correspondingly, the continental crust formation thickness is greater than oceanic crust formation thickness, herein Be set as continental crust formation thickness is twice of oceanic crust formation thickness, and oceanic crust stiffness layer should be greater than continental crust stratum, to embody ocean The characteristics of shell stratum is unlikely to deform.The underriding amount of every group of experiment, diving speed are all identical, only change underriding angle.

In one particular embodiment of the present invention, referring to Fig. 9 (a1) to Fig. 9 (c2), based on continental crust and oceanic crust and generation The theoretical experimental model of the underlying parameter of subduction, building plate subduction deformation structure specifically includes: by described based on continental crust And oceanic crust and generation subduction underlying parameter, the rock mechanics parameters for establishing the continental crust and oceanic crust stratum differentiation are set, And the theoretical experimental model to be dived with oceanic crust with different angle downwards.

In the present invention, theoretical experimental model includes following setting principle:

Continental crust stratum is reference with the relatively small sandstone formation of rigidity, and oceanic crust is larger with the biggish rigidity of non-deformability Sandstone formation be reference, and continental crust formation thickness be greater than oceanic crust stratum, oceanic crust setting have the scalable point in many places movable base Bottom has the property dived to specific direction.

Fig. 3 is the structural schematic diagram that simulation model constructs module 200, by based on continental crust and oceanic crust and generation subduction base Plinth parameter obtains the mechanics parameter and operation base prolonging direction set needed for theoretical analog lithostratigraphy, discrete to construct First numerical simulator, refers to Fig. 3, and the simulation model building module 200 includes:

Boundary constructs module 201, for constructing the bottom of subducting plate and continental crust plate according to the theoretical experimental model And lateral boundaries.

Specifically, keep stratum rock mechanics attribute and underriding angle constant, by Continental Crust Basement length, oceanic crust base length And proportionally ruler is contracted to suitable range to lateral boundaries height, simulates Continental Crust Basement length, oceanic crust base needed for obtaining Length, underriding strip length and lateral boundaries height.Such as the scale compression size according to 1:1000000.

Based on four Continental Crust Basement length, oceanic crust base length, lateral boundaries height and underriding angle parameters, mould is constructed The boundary of analog model.Constructed simulation model boundary have can to the movable subducting plate of specific direction and lateral boundaries, with It is identical as actual landform structure.

It specifically, is the structure for adapting to different simulation plate subduction deformations, the underriding angle and direction of subduction zone can be free Definition, subducting plate can extend to any position with special angle, and horizontal displacement same movement to the left is done together with boundary.

Mechanics parameter setup module 202, for be arranged discrete element numerical simulation model rock mechanics parameters and substrate with Boundary mechanics parameter.Coefficient of friction, elasticity modulus, normal direction and the tangential adhesion strength on the simulation continental crust stratum and oceanic crust stratum The corresponding mechanical attribute for referring to loose quartz sand respectively, does from childhood and big interruption value in certain valid interval, so that institute It states and forms simulation fault combination on simulated formation.

Specifically, the theoretical experimental model established by formation rock mechanical attribute, can obtain needed for simulation rock stratum The mechanics parameter of setting, correspondingly, the discrete element granular mechanics parameter setting of simulated formation should be with reference to by lot of experiment validation The mechanics parameter of Physical Experiment selection set, i.e., the described simulated formation mechanics parameter setting is with reference to physics experimental material pine Dissipate quartz sand.Since mechanical attribute is applied between microcosmic discrete element particle and particle by discrete element numerical simulation method, with reality Measured macroscopical mechanical parameters have different to material under experimental conditions, but difference is in a certain range, thus use with On the basis of macroscopical measuring parameter, to coefficient of friction, the bullet on simulation oceanic crust stratum and continental crust stratum in certain effective range Property modulus, forward direction and tangential adhesion strength do interval value, until obtain meet real material failure by rupture feature mechanics ginseng Number.

Simulated formation constructs module 203, and for being based on the bottom and lateral boundaries, building is deformed with the plate subduction The adaptable simulated formation of required geological conditions.The simulated formation is located within the bottom and lateral boundaries, simulated formation Parallel with constructed bottom boundary, left side continental crust stratum and right side oceanic crust stratum are in horizontal spread, middle part oceanic crust from bottom to top The oblique spread in stratum, each stratum and boundary have the contact relation for meeting mechanical properties of rock.

Specifically, the simulated formation includes: left side continental crust stratum, intermediate subduction of oceanic crust stratum and right side oceanic crust stratum. The simulated formation is used to simulate the true stratum in practical geological conditions.Simulated formation is parallel and tight with constructed bottom boundary Close to be filled within model boundary, left side continental crust stratum and the tiling of right side oceanic crust water flooding are set on continental crust and oceanic crust substrate, Middle part oceanic crust stratum is obliquely layed in boundary.For adapt to different condition simulated formation, the thickness of different simulated formations, Formation rock mechanical attribute in different simulated formations, the contact relation in different simulated formations can be defined freely on demand.

Simulated formation division module 204, for the simulated formation to be divided into the small of different colours label with equal thickness Layer.

Specifically, simulated formation is all made of is simulated corresponding to the discrete element particle of loose quartz sand mechanics parameter, and Oceanic crust stratum discrete element particle rigidity is larger, it is all it is particles filled after, particles at rest until all particles to reach contact steady It is fixed, a whole set of simulated formation is divided with same thickness, Different Strata marks be distinguish in different colors.When simulated formation is tied When structure changes, the track of structure change can be clearly reflected by the color mark.

Oceanic crust plate subduction module 205, for making oceanic crust plate dive with fixed speed and special angle to continental crust Effect.

Specifically, after simulation model is established, a fixed speed is applied to subduction of oceanic crust plate, makes subducting plate Do underriding movement downwards with specific direction, the side boundary to connect with it applies a speed identical with its horizontal component, makes It moves with uniform velocity to side, and other side boundary position remains unchanged.Correspondingly, under frictional force effect, oceanic crust stratum and one Lateral boundaries drive the stratum in its control range to move together, and oceanic crust stratum will drive and cover continental crust ground movement thereon.With Failure by rupture occurs for the increase of span, stratum, and the structural deformation of different patterns is developed on the continental crust stratum under different underriding angles.

Referring to Fig. 9 (a1) to Fig. 9 (c2), in order to which the depth for guaranteeing that oceanic crust has just started when diving is identical, and subduction zone two sides Oceanic crust is identical with the length of continental crust, therefore with the difference of subduction of oceanic crust angle, difference is arranged in the length of three groups of experimental models, the 1 group of experimental model is set as long 100m, and the 2nd group of experimental model is set as long 120m, and the 3rd group of experimental model is set as long 200m, But three groups of experimental model height are all 10m, all with 15% porosity, 2600kg/m³Density fill at random partial size 0.125m~ The two-dimensional disc particle of 0.075m, different-grain diameter particle quantity Gaussian distributed.Continental crust stratum is divided on the left of experimental model 10 layers, every thickness degree 1m, middle part and right side oceanic crust stratum are divided into 5 layers, and every thickness degree 1m is distinguish in different colors, side Just mature fault situation is observed.Right side oceanic crust substrate is equipped with movable wall body, is bowed to the left with different angle together with the wall of right side Punching, underriding rate 25%.In addition to key studies attribute difference, remaining parameter and boundary condition are consistent experimental model, with row Except uncorrelated factor influences.It is 2 × 10 by subject continental crust stratum elasticity modulus by testing repeatedly⁷Pa, oceanic crust stratum Elasticity modulus is 4 × 10⁸Pa, with the very much not easily-deformable feature of with significantly showing oceanic crust stiffness layer.Attribute is all provided between formation particles It is set to 0.6 coefficient of friction, makes simulation material on the basis of deferring to coulomb shear fracture criterion, forms a series of constructions that meet and recognize The fault combination form of knowledge.

Fig. 4 is the structural schematic diagram of plate subduction analysis module 300, the simulation based on the discrete element numerical simulation model As a result, the deformation behaviour of analysis plate underriding deformation process, refers to Fig. 4, the plate subduction analysis module 300 includes:

Analogue data determining module 301, for passing through the simulated process parameter for changing the discrete element numerical simulation model, To form different simulation fault styles, and obtain corresponding analogue data.

Simulated process parameter specifically: underriding mode, underriding amount, diving speed, substrate and the boundary of the simulation process Mechanics parameter, continental crust stratum and oceanic crust stratum thickness and mechanics parameter.

Vertic features Observation Blocks 302, the deformation for observing the corresponding plate subduction deformation process of different dive angles are special Sign.

Vertic features summarize module 303, for deforming law of development based on the simulated formation and tomography, are summarised in difference The differentiation deformation behaviour of plate subduction deformation process under angle.

In a specific embodiment of the invention, the deformation of plate subduction deformation process when observation dive angle is 30 ° Feature, the deformation behaviour of plate subduction deformation process when observation dive angle is 15 °, plate subduction deforms when observation dive angle is 5 ° The deformation behaviour of process deforms law of development based on the simulated formation and tomography, is summarised in plate subduction under different angle and becomes The differentiation deformation behaviour of shape process.

By taking plate subduction experimental model 1 as an example, referring to shown in Fig. 9 (a1) and Fig. 9 (a2), experimental setup shows model The initial stage of building, left and right side oceanic crust stratum are in a horizontal state spread, and the inclined spread in middle part oceanic crust stratum is bowed The angle of attack is 30 °, and right side and middle part oceanic crust stratigraphic boundary drive ocean jointly with the horizontal velocity right-to-left uniform motion of 0.5m/s Compressional movement is done on the above continental crust stratum in shell stratum, generates corresponding deformation.

By taking plate subduction experimental model 1 as an example, referring to shown in Figure 10 (b), it is squeezed and deformed early stage underriding amount reaches 5% Stage covers continental crust stratum and synchronizes and move to left as oceanic crust stratum moves to the left in drive, upper continental crust stratum is squeezed, fast Speed produces fracture deformation, and recoil fault F1, F2, F3 and the F4 for being inclined to continental crust are sequentially formed, with the increase of underriding amount, gradually shape At tomography F5, F6 of tendency oceanic crust.Oceanic crust mature fault is inclined to for entirety earlier than to recoil fault, and developing mode is preceding exhibition Formula continues to develop because the oceanic crust stratum on the right side of it limits tomography, so its preceding exhibition formula feature is unobvious, then development recoil Tomography, development mode are also preceding exhibition formula, are gradually developed to continental crust direction, and develop small-sized branch fault F7 thereon, and turn-off is small, main Play the role of adjusting continental crust stratum amount of compression.

By taking plate subduction experimental model 1 as an example, referring to shown in Figure 10 (c), squeezed in the early metaphase that underriding amount reaches 10% Deformation stage, recoil fault quantity increase, Zhan Shixiang continental crust Directional Extension before tomography, form tomography F8, and be previously formed Recoil fault be in nearly parallel relation, subsequently form small-sized branch fault；Oceanic crust tomography number is inclined to without significant change, and with Dive downwards stronger activity of oceanic crust influences, and tomography extends to deep, and its turn-off significantly increases；It forms underriding wedge on the whole Width and height increase.

By taking plate subduction experimental model 1 as an example, referring to shown in Figure 10 (d), squeezes and become in the mid-term that underriding amount reaches 15% Shape stage, recoil fault continue to develop, and tomography forms tomography F10, be in the recoil fault being previously formed to continental crust Directional Extension Nearly parallel relation；Oceanic crust tomography number is inclined to without significant change, turn-off continues to increase；On the whole its formed underriding wedge width and Height continues to increase.

By taking plate subduction experimental model 1 as an example, referring to shown in Figure 10 (e), squeezed in the middle and later periods that underriding amount reaches 20% Deformation stage, tomography quantity increased significantly, and recoil fault continues to develop, and sequentially form tomography F10, F11, F12, and this three disconnected Layer controls the lifting on stratum on the left of underriding wedge；Oceanic crust tomography number is inclined to without significant change, turn-off continues to increase；On the whole Underriding wedge width becomes larger, and left side stratum continues to be lifted, and form changes, and the left low and right high from before develops into left side than right The slightly higher wedge shape in side.

By taking plate subduction experimental model 1 as an example, referring to shown in Figure 10 (f), squeezes and become advanced stage underriding amount reaches 25% In the shape stage, without significant change, recoil fault F10, F11, F12 turn-off increases tomography quantity, this three faults controls underriding wedge Continue to be lifted in left side stratum；Oceanic crust tomography number is inclined to without significant change, turn-off continues to increase；Underriding wedge width on the whole Become larger, left side stratum continues to be lifted, and develops into the left side wedge shape higher than right side.

By taking plate subduction experimental model 1 as an example, referring to the oceanic crust of plate subduction deformation process shown in Figure 10 (a) to Figure 10 (f) Dive angle is 30 °, and underriding amount progressively increases to 25% from 0%, on the whole for, it is first former in continental crust and oceanic crust contact site Exhibition formula successively develops tendency continental crust tomography, and fault throw constantly increases, then before exhibition formula develop recoil fault in succession, tomography by Gradually to continental crust in-house development, small-sized branch fault is developed, when underriding amount reaches 13%, first control underriding wedge configuration Tomography F10 formed, be gradually increased underriding amount, when underriding amount reaches 16%, Article 2 controls the disconnected of underriding wedge configuration Layer F11 is formed, and continues to increase underriding amount, and when underriding amount reaches 17%, Article 3 controls the tomography F12 of underriding wedge configuration It is formed, since the underriding wedge height that subduction is formed constantly increases, width constantly becomes larger, and its form is gradually by left low and right high Develop into the left side wedge shape higher than right side.

By taking plate subduction experimental model 2 as an example, referring to shown in Fig. 9 (b1) to Fig. 9 (b2), the initial stage of model construction, Experimental setup subduction of oceanic crust angle is 15 °, and continental crust total length is 80m, and oceanic crust total length is 60m, other conditions and parameter setting with Experimental model 1 is identical, and right side and middle part oceanic crust stratigraphic boundary are jointly with the horizontal velocity right-to-left uniform motion of 0.5m/s.

By taking plate subduction experimental model 2 as an example, referring to plate subduction deformation process shown in Figure 11 (a) to Figure 11 (f), oceanic crust Dive angle is 15 °, and underriding amount progressively increases to 25% from 0%, and deformation process is similar with deformation process under 30 ° of dive angles, but Some differences of its deformation behaviour.In continental crust and oceanic crust contact site, first exhibition formula successively develops tendency continental crust tomography in the past, and tomography is disconnected Increase away from continuous, then before exhibition formula develop recoil fault in succession, tomography is gradually to continental crust in-house development, and small-sized branch fault Development, when underriding amount reaches 17%, the tomography F14 of first control underriding wedge configuration is formed, and is gradually increased underriding amount, When underriding amount reaches 22%, the tomography F15 of Article 2 control underriding wedge configuration is formed, and continues to increase underriding amount, when bowing When momentum reaches 25%, the tomography F16 of Article 3 control underriding wedge configuration is formed, the underriding wedge formed due to subduction Width constantly increases, height constantly increase, and its form gradually by left low and right high develop into the left and right sides height almost it is equal More regular wedge shape.

By taking plate subduction experimental model 3 as an example, referring to shown in Fig. 9 (c1) and Fig. 9 (c2), the initial rank of model construction Section, experimental setup subduction of oceanic crust angle are 5 °, and continental crust total length is 160m, and oceanic crust total length is 100m, and other conditions and parameter are set Set identical as experimental model 1 and experimental model 2, right side and middle part oceanic crust stratigraphic boundary are jointly with the horizontal velocity of 0.5m/s from right Uniform motion to the left.

By taking plate subduction experimental model 3 as an example, referring to shown in Figure 12 (b), it is squeezed and deformed early stage underriding amount reaches 5% In the stage, F1, F2, F3, F4, F5 tomography for being inclined to continental crust sequentially form, and with the increase of underriding amount, develop and small-sized incline on F1, F3 To oceanic crust branch fault F7, F6.Because its underriding angle is smaller, the influence to continental crust stratum is more uniform, so in underriding amount When smaller, be inclined to F1, F2, F3, F4, F5 tomography of continental crust near and meanwhile develop, secondly with the continuation of subduction, development There is the branch fault of small-sized tendency oceanic crust.

By taking plate subduction experimental model 3 as an example, referring to shown in Figure 12 (c), squeezed in the middle early stage that underriding amount reaches 10% Deformation stage, oceanic crust continue to dive downwards, and tomography quantity increased significantly, and be inclined to hair below F1, F2, F3, F4, F5 tomography of continental crust Educate away from its close parallel fault, the parallel fault of these close proximities has similar quality, so with same code name F1, F2, F3, F4, F5 indicate that the stratigraphic separation of both parallel faults of close proximity co- controlling significantly increases；Recoil fault largely increases It is more, develop tomography F8, F9, F10, F11, F12；It forms underriding wedge height and increases on the whole.

By taking plate subduction experimental model 3 as an example, referring to shown in Figure 12 (d), squeezes and become in the mid-term that underriding amount reaches 15% In the shape stage, with the continuous progress of subduction, tomography quantity increases, and the stratigraphic separation of F1, F2, F3, F4, F5 control continues to increase Greatly；Continue development tendency oceanic crust branch fault F14, recoil fault turn-off continues to increase；It is bright to form underriding wedge height for it on the whole It is aobvious to increase.

By taking plate subduction experimental model 3 as an example, referring to shown in Figure 12 (e), squeezed in the middle and later periods that underriding amount reaches 20% Deformation stage, recoil fault quantity increased significantly, and sequentially form tomography F15, F16, and this two faults controls underriding wedge is left The lifting on side stratum；Continental crust tomography quantity is inclined to without significant change, turn-off continues to increase；On the whole on the left of underriding wedge stratum after Continuous lifting.

By taking plate subduction experimental model 3 as an example, referring to shown in Figure 12 (f), squeezes and become advanced stage underriding amount reaches 25% Shape stage, tomography quantity significantly increases without significant change, recoil fault F15, F16 turn-off, this two faults controls underriding wedge Continue to be lifted in left side stratum；Continental crust tomography quantity is inclined to without significant change, turn-off continues to increase, and co-hade steepening；It is whole Spacing becomes smaller between different tomographies on body, and arrangement is even closer, and stratum continues to be lifted on the left of underriding wedge, and it is close to develop into the left and right gradient In equal wedge shape.

By taking plate subduction experimental model 3 as an example, referring to plate subduction deformation process shown in Figure 12 (a) to Figure 12 (f), oceanic crust Dive angle be 5 °, underriding amount progressively increases to 25% from 0%, on the whole for, first develop recoil fault, tomography on continental crust stratum Turn-off constantly increases, and then continental crust tomography is inclined in development in succession, and tomography is gradually to continental crust in-house development, small-sized branch fault hair It educates, when underriding amount reaches 20%, the tomography F15 of first control underriding wedge configuration is formed, it is gradually increased underriding amount, when When underriding amount reaches 23%, the tomography F16 of Article 2 control underriding wedge configuration is formed, the underriding formed due to subduction Wedge height constantly increases, and its form gradually develops into the left and right gradient by left low and right high and is bordering on equal wedge shape.

The structural map formed by comparing the final plate dive status of three groups of experiments, it can be seen that experienced 25% After underriding amount, plate subduction deformation fracture system has been developed in experiment 1, experiment 2, experiment 3.Three groups of Experimental comparisons can be summarized such as Lower rule: with the increase of underriding amount, tomography quantity is increasing, and whole turn-off is continuously increased；Under different underriding angles, bow Angle of attack degree is lower, and it is more to form tomography quantity；Underriding angle is lower, and the recoil fault of control underriding wedge configuration forms the time It is more late；Underriding angle is lower, and the height for forming underriding wedge is higher, and width is bigger, and the form for forming underriding wedge is different.

Specifically, this experimental result has preferably reappeared continental crust oceanic crust plate subduction deformation process geological phenomenon, makes plate Underriding deformation behaviour is associated with underriding angle.Such analog form can be used in analysis plate underriding angle to deformation behaviour It influences, to further enhance the understanding to this landform, provides theoretical direction foundation.

In an embodiment of the invention, continental crust formation thickness is greater than oceanic crust formation thickness, to guarantee the underriding Phenomenon meets reality；The oceanic crust stratum and oceanic crust stratum dynamics parameter setting refer to the loose quartz sand of physics experimental material, and With setting oceanic crust stiffness layer greater than continental crust stiffness layer, and be based on discrete element principle and appropriate adjustment is done to parameter.

In an embodiment of the invention, the analog result based on the discrete element numerical simulation model, analysis plates The deformation behaviour of block underriding deformation process deformation specifically includes:

Continental crust stratum on the left of simulation, intermediate subduction of oceanic crust band stratum and right side oceanic crust are sequentially filled in the model boundary The high density discrete element particle on stratum, the whole series simulated formation are divided into the substratum of different colours label with equal thickness；

Coefficient of friction, elasticity modulus, normal direction and the tangential adhesion strength on the simulation continental crust stratum and oceanic crust stratum are distinguished With reference to the corresponding mechanical attribute of loose quartz sand, done in certain valid interval from childhood and big interruption value, so that the mould Simulation fault combination is formed on quasi- stratum；

Change the simulation subduction zone underriding angle, and continue to do the mechanical attribute of the simulated formation interruption value, Until forming plate subduction deformation；

Deformation law of development based on the simulated formation and tomography obtains the simulated process parameter and bows with the plate Rush the corresponding relationship of deformation process.

The simulation continental crust stratum and the control of oceanic crust formation thickness ratio are in 2:1 or so, the simulation continental crust stratum elasticity Modulus is 2 × 10⁷Pa, coefficient of friction 0.6, normal direction and tangential adhesion strength are 3.6 × 10⁵Pa, the simulation oceanic crust stratum Elasticity modulus is 4 × 10⁸Pa, coefficient of friction 0.6, normal direction and tangential adhesion strength are 3.6 × 10⁵pa。

It as above is a kind of discrete element simulation system of plate subduction deformation process provided by the invention, according to actual Geological condition parameter constructs experimental model, and combines and construct discrete element by the theoretical parameter that the experimental model experimental calculation obtains Numerical simulator；Practical geologic structure is reappeared by discrete element numerical simulation modeling, is become by each parameter in simulation process Change the corresponding relationship with analog result, analyze influence of the current practical geological conditions to mature fault pattern, to be geology Research provides theoretical foundation.

In addition, although being referred to several unit modules of system in the above detailed description, it is this to divide only simultaneously Non-imposed.In fact, embodiment according to the present invention, the feature and function of two or more above-described units can To embody in a unit.Equally, the feature and function of an above-described unit can also be served as reasons with further division Multiple units embody.Terms used above " module " and " unit ", can be realize predetermined function software and/or Hardware.Although module described in following embodiment is preferably realized with software, the group of hardware or software and hardware The realization of conjunction is also that may and be contemplated.

After the discrete element simulation system for the plate subduction deformation process for describing exemplary embodiment of the invention, connect Get off, is introduced with reference to method of the attached drawing to exemplary embodiment of the invention.The implementation of this method may refer to above-mentioned whole The implementation of body, overlaps will not be repeated.

The present invention provides a kind of discrete element simulation method of plate subduction deformation process, is ground with solving the origin cause of formation of geologic structure The technical issues of studying carefully analytic process cannot verify, lack efficient association means between theory and the practical origin cause of formation.

Specifically, Fig. 5 is a kind of discrete element simulation method of plate subduction deformation process provided in an embodiment of the present invention Flow diagram, refers to Fig. 5, and the discrete element simulation method of the plate subduction deformation process includes:

S101: the geological conditions in the area studied as needed constructs the theoretical experimental model of plate subduction deformation structure.

S102: discrete element numerical simulation model is established based on the theoretical experimental model；

S103: the deformation behaviour based on the deformation of plate subduction described in the discrete element numerical simulation model analysis.

Fig. 6 is the specific flow chart of step S101, referring to Fig. 6, the step includes:

S201: determine based on continental crust and oceanic crust and occur the underlying parameter of subduction.

S202: the reason of plate subduction deformation is established based on setting principle according to the underlying parameter and the geological conditions By experimental model.

Fig. 7 is the specific flow chart that discrete element numerical simulation model is established based on the theoretical experimental model, by being based on land Shell and oceanic crust and generation subduction underlying parameter obtain the mechanics parameter set needed for theoretical analog lithostratigraphy and activity Base prolonging direction refers to Fig. 7 constructing discrete element numerical simulation model, based on the theoretical experimental model establish from Dissipating member numerical simulator includes:

S301: according to the theoretical experimental model, bottom and the lateral boundaries of subducting plate and continental crust plate are constructed.

S302: the rock mechanics parameters and substrate and boundary mechanics parameter of setting discrete element numerical simulation model.The mould Coefficient of friction, elasticity modulus, normal direction and the tangential adhesion strength on quasi- continental crust stratum and oceanic crust stratum refer to loose quartz sand respectively Corresponding mechanical attribute, done in certain valid interval from childhood and big interruption value, so as to be formed on the simulated formation Simulate fault combination.

S303: being based on the bottom and lateral boundaries, and geological conditions needed for building is deformed with the plate subduction is adapted Simulated formation.The simulated formation is located within the bottom and lateral boundaries, simulated formation and constructed bottom boundary phase In parallel, left side continental crust stratum and right side oceanic crust stratum are in horizontal spread from bottom to top, the oblique spread in middle part oceanic crust stratum, respectively Stratum and boundary have the contact relation for meeting mechanical properties of rock.

S304: the simulated formation is divided into the substratum of different colours label with equal thickness.

S305: make oceanic crust plate that subduction occur to continental crust with fixed speed and special angle.

Fig. 8 is that the process of the deformation behaviour deformed based on plate subduction described in the discrete element numerical simulation model analysis is shown It is intended to, refers to Fig. 8, the deformation behaviour based on the deformation of plate subduction described in the discrete element numerical simulation model analysis includes:

S401: the simulated process parameter by changing the discrete element numerical simulation model, it is disconnected to form different simulations Layer pattern, and obtain corresponding analogue data.

S402: the deformation behaviour of the corresponding plate subduction deformation process of different dive angles is observed.

S403: deforming law of development based on the simulated formation and tomography, is summarised in plate subduction under different angle and deforms The differentiation deformation behaviour of process.

The present invention also provides a kind of computer equipments, comprising: it is adapted for carrying out the processor and storage equipment of each instruction, The storage equipment is stored with a plurality of instruction, and described instruction is suitable for being loaded by processor and executing a kind of plate subduction deformation process Discrete element simulation method.

The present invention also provides a kind of computer readable storage mediums, are stored with computer program, the computer program For executing a kind of discrete element simulation method of plate subduction deformation process.

Beneficial effects of the present invention are mainly as follows:

(1) building mode of experimental model is simple, and the discrete element numerical simulation model based on experimental model building, mould Quasi- result can be good at reappearing the structure of practical geological fault, pair for being changed and being formed by between fault structure with certificate parameter It should be related to, can be used in the impact analysis and prediction according to actual current geological conditions to mature fault mode；

(2) parameter involved in building experimental model is few, and calculating process is simple；

(3) building mode of discrete element numerical simulation model is simple, it is easy to accomplish, compared to physical analogy model, not by The limitation of experimental material and modelling can freely define experimental model according to demand；

(4) in simulation process and final state model can be saved at any time, extract speed and displacement etc. analysis data, It is sliced by any position；

(5) variation of simulated conditions can easily be controlled according to actual geological conditions during realistic simulation, implemented Process is simple to operation, to the favorable reproducibility of actual tomography.

It is improvement on hardware (for example, to diode, crystal that the improvement of one technology, which can be distinguished clearly, Pipe, switch etc. circuit structures improvement) or software on improvement (improvement for method flow).However, with technology The improvement of development, current many method flows can be considered as directly improving for hardware circuit.Designer is almost All corresponding hardware circuit is obtained by the way that improved method flow to be programmed into hardware circuit.Therefore, it cannot be said that one The improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device (Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate Array, FPGA)) just It is such a integrated circuit, logic function determines device programming by user.It is voluntarily programmed by designer Lai one Dedicated integrated circuit is designed without asking chip maker and made to a digital display circuit " integrated " on a piece of PLD Chip.Moreover, nowadays, substitution manually makes IC chip, and " logic compiler (logic is also used in this programming instead mostly Compiler) " software realizes that when it writes with program development software compiler used is similar, and before compiling Also handy specific programming language is write for source code, this is referred to as hardware description language (Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kinds, such as ABEL (Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、 CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language) etc., VHDL (Very-High-Speed Integrated Circuit is most generally used at present Hardware Description Language) and Verilog2.Those skilled in the art also will be apparent to the skilled artisan that only needs will be square Method process slightly programming in logic and is programmed into integrated circuit with above-mentioned several hardware description languages, so that it may be readily available reality The now hardware circuit of the logical method process.

Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing The computer for the computer readable program code (such as software or firmware) that device and storage can be executed by (micro-) processor can Read medium, logic gate, switch, specific integrated circuit (Application Specific Integrated Circuit, ASIC), the form of programmable logic controller (PLC) and insertion microcontroller, the example of controller includes but is not limited to following microcontroller Device: ARC625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320 are deposited Memory controller is also implemented as a part of the control logic of memory.

It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind Hardware component, and the structure that the device for realizing various functions for including in it can also be considered as in hardware component.Or Even, can will be considered as realizing the device of various functions either the software module of implementation method can be Hardware Subdivision again Structure in part.

System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity, Or it is realized by the product with certain function.

For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this The function of each unit can be realized in the same or multiple software and or hardware when application.

As seen through the above description of the embodiments, those skilled in the art can be understood that the application can It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer system (can be personal computer, server or network system etc.) executes the certain of each embodiment of the application or embodiment Method described in part.

All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.

The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, clothes Business device computer, hand system or portable system, plate system, multicomputer system, microprocessor-based system, set Top box, programmable consumer electronics system, network PC, minicomputer, mainframe computer including any of the above system or system Distributed computing environment etc..

The application can describe in the general context of computer-executable instructions executed by a computer, such as program Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by Task is executed by the connected teleprocessing system of communication network.In a distributed computing environment, program module can be with In the local and remote computer storage media including storage system.

Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's Spirit.

Claims

1. a kind of discrete element simulation method of plate subduction deformation process, which is characterized in that the described method includes:

2. the method according to claim 1, wherein the geological conditions in the area studied as needed, structure The theoretical experimental model for building plate subduction deformation structure includes:

The theoretical experiment mould of plate subduction deformation is established based on setting principle according to the underlying parameter and the geological conditions Type.

3. according to the method described in claim 2, it is characterized in that, the underlying parameter includes continental crust formation thickness, oceanic crust Thickness degree, continental crust stratum and oceanic crust stratum mechanics parameter, underriding amount, diving speed and underriding angle；The setting principle The sandstone formation of first threshold is less than for reference with rigidity including the continental crust, the oceanic crust is greater than the second threshold with non-deformability The sandstone formation that value and rigidity are greater than third threshold value is reference, and the oceanic crust has the moveable substrate of the scalable point in many places, and has The property that oriented any direction dives.

4. according to the method described in claim 2, it is characterized in that, establishing discrete element Numerical-Mode based on the theoretical experimental model Analog model includes:

Based on the bottom and lateral boundaries, the adaptable simulation of geological conditions needed for building is deformed with the plate subduction Layer；

5. according to the method described in claim 4, it is characterized in that, based on plate described in the discrete element numerical simulation model analysis Block dive deformation deformation behaviour include:

By changing the simulated process parameter of the discrete element numerical simulation model, to form different simulation fault styles, and Obtain corresponding analogue data；

Law of development is deformed based on the simulated formation and tomography, is summarised in the difference of plate subduction deformation process under different angle Change deformation behaviour.

6. a kind of discrete element simulation system of plate subduction deformation process, which is characterized in that the system comprises:

Experimental model constructs module, and the geological conditions in the area for studying as needed constructs plate subduction deformation structure Theoretical experimental model；

Plate subduction analysis module, the deformation for being deformed based on plate subduction described in the discrete element numerical simulation model analysis Feature.

7. system according to claim 6, which is characterized in that the experimental model constructs module and includes:

Experimental model establishes module, establishes plate for being based on setting principle according to the underlying parameter and the geological conditions Dive the theoretical experimental model deformed.

8. system according to claim 7, which is characterized in that the underlying parameter includes continental crust formation thickness, oceanic crust Thickness degree, continental crust stratum and oceanic crust stratum mechanics parameter, underriding amount, diving speed and underriding angle；The setting principle The sandstone formation of first threshold is less than for reference with rigidity including the continental crust, the oceanic crust is greater than the second threshold with non-deformability The sandstone formation that value and rigidity are greater than third threshold value is reference, and the oceanic crust has the moveable substrate of the scalable point in many places, and has The property that oriented any direction dives.

9. system according to claim 7, which is characterized in that the simulation model constructs module and includes:

Boundary constructs module, for according to the theoretical experimental model, constructing the bottom of subducting plate and continental crust plate and lateral Boundary；

Mechanics parameter setup module, for the rock mechanics parameters and substrate and boundary mechanics of discrete element numerical simulation model to be arranged Parameter；

Simulated formation construct module, for be based on the bottom and lateral boundaries, building with the plate subduction deformation needed to The adaptable simulated formation of matter condition；

Oceanic crust plate subduction module, for making oceanic crust plate that subduction occur to continental crust with fixed speed and special angle.

10. system according to claim 9, which is characterized in that the plate subduction analysis module includes:

Analogue data determining module, for the simulated process parameter by changing the discrete element numerical simulation model, to be formed Different simulation fault styles, and obtain corresponding analogue data；

Vertic features summarize module, for deforming law of development based on the simulated formation and tomography, are summarised under different angle The differentiation deformation behaviour of plate subduction deformation process.

11. a kind of computer equipment characterized by comprising it is adapted for carrying out the processor and storage equipment of each instruction, it is described Storage equipment is stored with a plurality of instruction, and described instruction is suitable for being loaded by processor and being executed such as claim 1 to 5 any one institute A kind of discrete element simulation method for the plate subduction deformation process stated.

12. a kind of computer readable storage medium, which is characterized in that be stored with computer program, the computer program is used for Execute a kind of discrete element simulation method of plate subduction deformation process as described in claim 1 to 5 any one.