CN113722967B - Sedimentary geological action process simulation method and device - Google Patents
Sedimentary geological action process simulation method and device Download PDFInfo
- Publication number
- CN113722967B CN113722967B CN202110885025.4A CN202110885025A CN113722967B CN 113722967 B CN113722967 B CN 113722967B CN 202110885025 A CN202110885025 A CN 202110885025A CN 113722967 B CN113722967 B CN 113722967B
- Authority
- CN
- China
- Prior art keywords
- model
- particles
- updated
- stratum model
- stratum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 230000009471 action Effects 0.000 title claims abstract description 54
- 230000008569 process Effects 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 22
- 239000002245 particle Substances 0.000 claims abstract description 198
- 230000008021 deposition Effects 0.000 claims abstract description 87
- 239000002344 surface layer Substances 0.000 claims abstract description 31
- 238000004062 sedimentation Methods 0.000 claims abstract description 23
- 238000001125 extrusion Methods 0.000 claims abstract description 14
- 230000009467 reduction Effects 0.000 claims abstract description 13
- 239000000758 substrate Substances 0.000 claims abstract description 9
- 230000000694 effects Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 abstract description 30
- 230000007246 mechanism Effects 0.000 abstract description 4
- 238000005137 deposition process Methods 0.000 description 5
- 239000010410 layer Substances 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a sedimentary geological action process simulation method and a device thereof, wherein the method comprises the following specific operations: obtaining a stratum model with initial dead weight balance, and determining a reference plane of sedimentation in the stratum model; changing the position of surface layer particles in the stratum model through extrusion loading and substrate sedimentation to obtain an updated stratum model; generating a model speed linked list at the current moment according to the updated stratum model, and acquiring a deposition lower boundary; generating deposition particles in the updated stratum model according to the deposition lower boundary and the reference plane; obtaining a dense and flush sedimentary formation with the original reference surface by defining a particle generation range and a sedimentary particle velocity; and carrying out speed reduction on the sedimentary stratum model according to the model speed linked list to finish one-time sedimentary action. The invention can truly simulate the sedimentary geological action process, reveal the formation mechanism of different geological structures, and reflect the occurrence process of the sedimentary action from the aspect of fine view.
Description
Technical Field
The invention relates to a sediment geological action process simulation method and a device thereof, belonging to the technical field of geological structure motion analysis.
Background
In geology, geologic structure movement is a complex process, and common geologic forms often develop through lengthy geologic actions, including various geologic action processes such as deposition, ablation, plate extrusion movement, erosion handling, etc. The formation process of the geological morphology is discussed, and the method has extremely important significance for geotechnical engineering, petroleum engineering and the like. Because the scale of the geological structure is often larger, the difficulty of establishing a physical model is higher, and the numerical simulation method has obvious advantages for simulating a large-scale geological engineering model, in particular to a particle discrete element numerical simulation method which can discuss the formation mechanism of the geological structure movement from a microscopic angle. In general, the geological action processes are not independent, and often are the result of the combined action of a plurality of geological actions, so that the simulation of the geological structure movement also often needs to consider the plurality of geological actions, and at present, how to comprehensively simulate various geological action processes by a numerical simulation method still lacks relevant technical support.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a sedimentary geological action process simulation method and a sedimentary geological action process simulation device, which utilize particle discrete elements to simulate sedimentary geological action processes, fully consider the influence of other geological actions in the simulation process, reflect the real and reliable sedimentary action processes, reveal the formation mechanism of geological structures from a microscopic angle, and provide reliable technical support for geotechnical engineering, geological engineering, petroleum engineering and the like.
In order to solve the technical problems, the invention adopts the following technical means:
in a first aspect, the present invention provides a method for simulating a sedimentary geologic process, comprising the steps of:
generating particles in a preset boundary wall, obtaining a stratum model with initial dead weight balance by adjusting the position of the boundary wall, and determining a reference surface for sedimentation in the stratum model;
changing the position of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time step number;
generating a model speed linked list at the current moment according to the updated stratum model, and acquiring a deposition lower boundary;
judging the vertical relative position according to the lower deposition boundary and the reference plane, generating deposition particles in the updated stratum model when the judgment condition is met, and compacting the deposition particles under the action of dead weight by limiting the speed of the deposition particles to obtain a deposition stratum model;
and carrying out speed reduction on the sedimentary stratum model according to the model speed linked list to finish one-time sedimentary action.
With reference to the first aspect, further, the method for obtaining an updated formation model includes the steps of:
removing an upper boundary wall of the stratum model based on the stratum model with initial dead weight balance, applying horizontal opposite speeds to left and right boundary walls of the stratum model, and simulating the extrusion loading effect of the stratum;
applying vertical downward speed to the bottom boundary wall of the stratum model, simulating the sedimentation effect of the stratum, controlling the speeds of left and right vertexes of the bottom boundary wall, and simulating the uneven sedimentation effect of the stratum;
the particles in the stratum model are driven to freely move through extrusion loading action, sedimentation action and uneven sedimentation action, so that the positions of the surface layer particles in the stratum model are changed;
and when the movement time of the stratum model particles reaches the preset simulated movement time step number, obtaining an updated stratum model.
With reference to the first aspect, further, the model speed linked list includes speeds of all particles in the updated formation model at the current time, speeds of the updated formation model left boundary wall, speeds of the updated formation model right boundary wall, and speeds of the updated formation model bottom boundary wall.
With reference to the first aspect, further, the method for obtaining the lower deposition boundary includes:
taking a horizontal range in the horizontal direction, and comparing the particle positions in the horizontal range to obtain particles with the largest vertical positions in the horizontal range, and marking the particles as surface layer particles of the stratum model updated in the horizontal range, wherein the particle diameters of the particles in the horizontal range are smaller than or equal to 2;
continuously selecting a new horizontal range, and acquiring particles with the largest vertical position in each horizontal range until all horizontal particles in the updated stratum model are traversed, so as to acquire all surface layer particles of the updated stratum model;
and sequentially comparing the vertical position coordinates of each surface layer particle, obtaining the surface layer particle with the minimum vertical position, and taking the vertical position coordinates of the surface layer particle as a deposition lower boundary.
With reference to the first aspect, further, the method for determining the vertical relative position according to the deposition lower boundary and the reference plane includes:
calculating the position difference of the deposition lower boundary in the vertical direction according to the vertical position coordinate of the deposition lower boundary and the reference plane to obtain a vertical relative position;
and when the vertical relative position is larger than 2-4 particle diameters, generating deposition particles in the updated stratum model, otherwise, not generating deposition particles in the updated stratum model.
With reference to the first aspect, further, the method for generating the sediment particles in the updated stratum model includes:
acquiring the particle generation range of the updated stratum model according to the deposition lower boundary and the reference surface: (x) 1 ～x 2 ,y 1 ～y 2 ) Wherein x is 1 Representing the horizontal position coordinate, x, of the left boundary wall of the updated stratum model 2 Representing the horizontal position coordinate, y, of the right boundary wall of the updated stratum model 1 Vertical position coordinates representing the lower boundary of the deposit, y 2 Representing the sum of the vertical position coordinates of the reference surface and 2-4 particle diameters;
the deposited particles are generated by adopting a random filling and non-overlapping method in the particle generation range.
With reference to the first aspect, further, the method for obtaining the sedimentary stratum model comprises the following steps:
obtaining deposition calculation time steps according to the diameter and density of the deposited particles;
enabling the deposited particles to move under the action of dead weight, detecting the speed of the deposited particles, and returning the speed of the deposited particles to zero when the speed of the deposited particles is greater than a preset value;
when the movement time of the deposited particles reaches the step number of the deposition calculation time, the deposited particles are densely accumulated under the action of dead weight, and a deposited stratum model is obtained.
With reference to the first aspect, further, the method for performing velocity reduction on the sedimentary stratum model includes:
grouping new deposited particles and existing particles in the updated stratum model to obtain a new deposited particle group and an old deposited particle group;
and reducing the speed of each particle in the old sedimentary particle group according to the model speed linked list, and reducing the speed of a boundary wall in the sedimentary stratum model.
In a second aspect, the present invention provides a sedimentary geologic process simulation device, comprising:
and the stratum model initializing module is used for generating particles in a preset boundary wall, obtaining a stratum model with initial dead weight balance by adjusting the position of the boundary wall, and determining a deposition reference surface in the stratum model.
The stratum model updating module is used for changing the positions of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time steps;
the linked list module is used for generating a model speed linked list at the current moment according to the updated stratum model;
the particle screening module is used for acquiring a deposition lower boundary according to the updated stratum model;
the deposition module is used for judging the vertical relative position according to the lower deposition boundary and the reference plane, generating deposition particles in the updated stratum model when the judgment condition is met, and compacting the deposition particles under the action of dead weight by limiting the speed of the deposition particles to obtain a deposition stratum model;
and the speed reduction module is used for carrying out speed reduction on the sedimentary stratum model according to the model speed linked list so as to finish one-time sedimentary action.
With reference to the second aspect, further, the model velocity linked list includes velocities of all particles in the updated formation model at the current time, velocities of the left boundary wall of the updated formation model, velocities of the right boundary wall of the updated formation model, and velocities of the bottom boundary wall of the updated formation model.
The following advantages can be obtained by adopting the technical means:
the invention provides a sedimentary geological action process simulation method and a device thereof, wherein the sedimentary geological action process is simulated through a particle discrete element numerical model. In the particle deposition process, the invention limits the generation range of deposition particles, ensures that the surface layer of the newly generated deposition particles after each deposition is basically consistent with a reference surface, and ensures that the deposition particles are extruded and compacted by limiting the particle speed without loosening or bouncing off. After each deposition, the method and the device perform speed reduction on the stratum model, so that other structural movements of the whole stratum model can be ensured not to be influenced by the deposition process. The method and the device provide reliable technical support for geotechnical engineering, geological engineering, petroleum engineering and the like.
Drawings
FIG. 1 is a flow chart of steps of a sedimentary geologic process simulation method of the present invention;
FIG. 2 is a schematic view of a formation model according to an embodiment of the present invention;
FIG. 3 is a schematic view of loading a formation model according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of an updated earth model according to an embodiment of the present invention;
FIG. 5 is a schematic illustration of the generation of deposited particles in an embodiment of the present invention;
FIG. 6 is a schematic illustration of dense deposition of deposited particles in an embodiment of the present invention;
FIG. 7 is a graph showing the effect of deposition in an embodiment of the present invention;
FIG. 8 is a schematic diagram of a sedimentary geologic process simulation device according to the present invention;
in the figure, 1 is a left boundary wall, 2 is a right boundary wall, 1 is an upper boundary wall, 4 is a bottom boundary wall, 5 is a reference surface, 6 is surface layer particles of the updated formation model, 7 is a newly deposited particle group, and 8 is an old deposited particle group.
Detailed Description
The technical scheme of the invention is further described below with reference to the accompanying drawings:
the invention provides a sedimentary geological action process simulation method, which is based on particle flow software for carrying out sedimentary geological action simulation, as shown in figure 1, and specifically comprises the following steps:
and A, generating particles in a preset boundary wall, obtaining a stratum model with initial dead weight balance by adjusting the position of the boundary wall, and determining a reference surface of sedimentation in the stratum model.
Step A01, an initial stratum model is composed of internal particles and external (up, down, left and right) boundary walls, and the stratum model is divided into two different materials: the energy dry layer and slip layer, as shown in fig. 2, the light formation represents slip layer and the dark formation represents energy dry layer.
After particles are generated in the boundary wall, parameters are given to the stratum model, the internal stress of the stratum model is uniformly distributed by continuously moving the position of the boundary wall, then the gravity action is applied, the balance convergence is solved, and the stratum model with initial dead weight balance is obtained.
And A02, traversing all particles in the initial dead weight balanced stratum model, judging whether the particles belong to surface particles according to the vertical position coordinates of the particles, and obtaining the surface particles of the initial dead weight balanced stratum model.
In step a03, in the formation model with initial weight balance, the surface layer particles are generally approximately in a horizontal plane, so that the vertical position coordinates of the surface layer particles in the initial state are set as the reference plane for sedimentation in the formation model.
And B, changing the position of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time step number, wherein the specific operation is as follows as shown in figures 2 and 3:
and B01, removing the upper boundary wall 3 of the stratum model based on the stratum model with initial dead weight balance so as to ensure free deformation of the upper boundary of the stratum model. Applying a horizontal facing velocity V to the left and right boundary walls 1, 2 of the earth model 1 、V 2 The extrusion loading effect of the stratum is simulated by the opposite movement of the left and right boundary walls. V (V) 1 、V 2 Is set manually, in general, V 2 And taking 0, namely fixing the right boundary wall, and loading through the left boundary wall.
Step B02, applying a vertical downward speed to the bottom boundary wall 4 of the stratum model, simulating the deposition of the stratum, and simultaneously controlling the speeds V of the left and right vertexes of the bottom boundary wall 3 、V 4 To simulate the settlement of the substrate, V 3 、V 4 When equal, mean that the substrate is uniformly settled, V 3 、V 4 Non-uniform settlement of the formation may be simulated when not equal. V (V) 3 、V 4 Is set manually.
And B03, under the simulation of the steps B01 and B02, driving particles in the stratum model to freely move through the extrusion loading effect, the sedimentation effect and the uneven sedimentation effect, so as to change the positions of surface layer particles in the stratum model.
And B04, when the movement time of the stratum model particles reaches the preset simulated movement time step number (10000 steps are taken in the embodiment of the invention), obtaining an updated stratum model.
To ensure computational efficiency, deposition does not occur every time step, but is typically detected at regular time steps, with the specific number of time steps being related to the size and loading rate of the formation model.
And C, generating a model speed linked list at the current moment according to the updated stratum model, and acquiring a deposition lower boundary.
In order to ensure that other structural movements of the whole formation model are not affected by the deposition process, each velocity in the formation model at the current time is required to be stored before each generation of new deposition particles, and the model velocity linked list includes the velocities of all particles in the formation model updated at the current time, the velocity of the left boundary wall of the formation model updated, the velocity of the right boundary wall of the formation model updated and the velocity of the bottom boundary wall of the formation model updated.
In step C, the method for obtaining the lower deposition boundary is as follows:
step C01, traversing all particles in the updated stratum model, and screening out surface layer particles 6 of the updated stratum model according to the vertical position coordinates of each particle; specific: taking a horizontal range in the horizontal direction, wherein the horizontal range is smaller than or equal to 2 particle diameters, and determining the particles with the largest vertical positions in the horizontal range by searching and comparing the particle positions in the horizontal range, and marking the particles as surface layer particles of the stratum model updated in the horizontal range; and continuously selecting a new horizontal range, searching the particles with the largest vertical positions in each horizontal range until all horizontal particles in the stratum model are traversed, and obtaining all surface particles of the updated stratum model.
And C02, comparing the vertical position coordinates of each surface layer particle in sequence, acquiring the surface layer particle with the minimum vertical position from all the surface layer particles of the updated stratum model, and taking the vertical position coordinates of the surface layer particle as a deposition lower boundary. The relative distance between the surface particles with the smallest vertical position and the reference surface 5 is the largest.
And D, judging the vertical relative positions according to the lower deposition boundary and the reference plane, generating deposition particles in the updated stratum model when the judgment conditions are met, and compacting the deposition particles under the action of dead weight by limiting the speed of the deposition particles to obtain the deposition stratum model.
And D01, calculating the position difference of the vertical position coordinate of the lower deposition boundary and the reference plane in the vertical direction to obtain the vertical relative position.
And D02, when the vertical relative position is larger than 2-4 particle diameters, generating deposition particles in the updated stratum model, otherwise, not generating the deposition particles in the updated stratum model.
Step D03, obtaining updated data according to the lower deposition boundary and the reference planeParticle generation range of the stratum model: (x) 1 ～x 2 ,y 1 ～y 2 ) Wherein x is 1 Representing the horizontal position coordinate, x, of the left boundary wall of the updated stratum model 2 Representing the horizontal position coordinate, y, of the right boundary wall of the updated stratum model 1 Vertical position coordinates representing the lower boundary of the deposit, y 2 The sum of the vertical position coordinates of the reference surface and 2-4 particle diameters is represented.
And D04, setting the speeds of the left boundary wall, the right boundary wall and the bottom boundary wall to be 0, and generating deposited particles in a random filling and non-overlapping mode in a particle generation range, as shown in fig. 5.
And D05, under the condition that other particle microscopic parameters are consistent, the calculated deposition time step number is mainly related to the particle size and density, so that the deposition calculation time step number of the particle flow software PFC is set according to the diameter and density of the deposited particles, and the speed of the deposited particles is limited. To correspond to real geologic motion, the calculated time steps in discrete elements of a particle are typically linked to the actual geologic time, one time step typically representing decades, and the calculated number of time steps represents how many time steps need to be calculated.
Step D06, enabling the deposited particles to move under the action of dead weight, calculating the total force of the deposited particles, detecting the speed of the deposited particles, and when the speed of the deposited particles is larger than a preset value (a smaller value is generally selected, and the embodiment of the invention takes 5 m/s), enabling the speed of the deposited particles to return to zero, so that the speed limit of the deposited particles is realized.
Step D07, repeating step D06 until the movement time of the deposited particles reaches the set deposition calculation time step number (15000 steps are taken in the embodiment of the invention), and at this time, the deposited particles are densely packed under the action of dead weight, as shown in fig. 6, so as to obtain a deposited stratum model.
The invention always limits the speed of the deposited particles below 5m/s, thereby ensuring that the newly generated deposited particles do not generate great collision speed when contacting with the existing particles in the stratum model, enabling the stratum model to reach an equilibrium state quickly, and redefining the newly generated deposited particles into the old deposited particle group after reaching the equilibrium state. In addition, the invention limits the generation range of the deposited particles, so that the surface layer of the newly generated deposited particles is basically consistent with the reference surface after the deposited particles are densely packed.
And E, performing speed reduction on the sedimentary stratum model according to the model speed linked list to finish one-time sedimentary action.
In the particle deposition process, considering that the initial velocity of newly formed deposited particles is zero, the velocity of existing particles in the stratum model develops with the development of geological motion, and there is a velocity difference between the two, so that newly formed deposited particles and previous particles must be separately grouped, i.e., divided into a new deposited particle group 7 and an old deposited particle group 8.
Restoring the particle velocity (old deposited particle group) stored in the model velocity linked list to each particle, and restoring the velocities of the left, right and bottom boundary walls to the velocity V before deposition 1 、V 2 、V 3 、V 4 One deposition is completed as shown in fig. 7.
And (C) circulating the steps B-E, and continuing the operations such as model loading, time step detection and the like, so that a continuous deposition process can be realized.
In the loading process (step B) of the stratum model, the left and right boundary walls have speeds, so that the left and right boundary walls can generate transverse displacement along with the increment of time steps, and when the transverse displacement generated by the left and right boundary walls reaches a specified size, the deposition motion state of the stratum model is recorded once, so that the stratum morphology and the deposition action process of different stages can be conveniently observed and compared. In addition, the invention can also be used as the basis for recording the movement state of the stratum model in other forms, such as time step number, calculated duration and the like.
The invention also provides a sedimentary geological action process simulation device, which comprises a stratum model initialization module, a stratum model updating module, a linked list module, a particle screening module, a sedimentary module and a speed reduction module as shown in fig. 8.
The stratum model initialization module is mainly used for generating particles in a preset boundary wall, obtaining an initial dead weight balanced stratum model by adjusting the position of the boundary wall, and determining a deposition reference surface in the stratum model. The specific operation of the stratigraphic model initialization module is consistent with step a of the method of the present invention.
The stratum model updating module is mainly used for changing the positions of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time steps. The specific operation of the stratigraphic model updating module is consistent with step B of the method of the present invention.
The linked list module is mainly used for generating a model speed linked list at the current moment according to the updated stratum model, wherein the model speed linked list comprises the speeds of all particles in the stratum model updated at the current moment, the speeds of the left boundary wall of the stratum model updated, the speeds of the right boundary wall of the stratum model updated and the speeds of the bottom boundary wall of the stratum model updated; the particle screening module is mainly used for acquiring a deposition lower boundary according to the updated stratum model. The specific operation of the linked list module and the particle screening module is consistent with step C of the method of the present invention.
The deposition module is mainly used for judging vertical relative positions according to a deposition lower boundary and a reference plane, generating deposition particles in the updated stratum model when the judgment conditions are met, and enabling the deposition particles to be piled up tightly under the action of dead weight by limiting the speed of the deposition particles so as to obtain the deposition stratum model. The specific operation of the deposition module corresponds to step D of the method of the invention.
The speed reduction module is mainly used for carrying out speed reduction on the sedimentary stratum model according to the model speed linked list so as to finish one-time sedimentary action.
The method and the device fully consider the influence of other geological actions, can truly simulate the sedimentary geological action process, reveal the formation mechanism of different geological structures, reflect the occurrence process of the sedimentary action from a fine view, and provide reliable technical support for geotechnical engineering, geological engineering, petroleum engineering and the like.
The foregoing is merely a preferred embodiment of the present invention, and it should be noted that modifications and variations could be made by those skilled in the art without departing from the technical principles of the present invention, and such modifications and variations should also be regarded as being within the scope of the invention.
Claims (8)
1. A method for simulating a sedimentary geologic process, comprising the steps of:
generating particles in a preset boundary wall, obtaining a stratum model with initial dead weight balance by adjusting the position of the boundary wall, and determining a reference surface for sedimentation in the stratum model;
changing the position of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time step number;
generating a model speed linked list at the current moment according to the updated stratum model, and acquiring a deposition lower boundary; the model speed linked list comprises the speeds of all particles in the stratum model updated at the current moment, the speeds of the left boundary wall of the stratum model updated, the speeds of the right boundary wall of the stratum model updated and the speeds of the bottom boundary wall of the stratum model updated;
judging the vertical relative position according to the lower deposition boundary and the reference plane, generating deposition particles in the updated stratum model when the judgment condition is met, and compacting the deposition particles under the action of dead weight by limiting the speed of the deposition particles to obtain a deposition stratum model;
and carrying out speed reduction on the sedimentary stratum model according to the model speed linked list to finish one-time sedimentary action.
2. A method of sedimentary geologic process simulation according to claim 1, wherein the method of obtaining an updated earth model comprises the steps of:
removing an upper boundary wall of the stratum model based on the stratum model with initial dead weight balance, applying horizontal opposite speeds to left and right boundary walls of the stratum model, and simulating the extrusion loading effect of the stratum;
applying vertical downward speed to the bottom boundary wall of the stratum model, simulating the sedimentation effect of the stratum, controlling the speeds of left and right vertexes of the bottom boundary wall, and simulating the uneven sedimentation effect of the stratum;
the particles in the stratum model are driven to freely move through extrusion loading action, sedimentation action and uneven sedimentation action, so that the positions of the surface layer particles in the stratum model are changed;
and when the movement time of the stratum model particles reaches the preset simulated movement time step number, obtaining an updated stratum model.
3. The method for simulating a sedimentary geologic process of claim 1, wherein the method for obtaining the sedimentary lower boundary comprises:
taking a horizontal range in the horizontal direction, and comparing the particle positions in the horizontal range to obtain particles with the largest vertical positions in the horizontal range, and marking the particles as surface layer particles of the stratum model updated in the horizontal range, wherein the particle diameters of the particles in the horizontal range are smaller than or equal to 2;
continuously selecting a new horizontal range, and acquiring particles with the largest vertical position in each horizontal range until all horizontal particles in the updated stratum model are traversed, so as to acquire all surface layer particles of the updated stratum model;
and sequentially comparing the vertical position coordinates of each surface layer particle, obtaining the surface layer particle with the minimum vertical position, and taking the vertical position coordinates of the surface layer particle as a deposition lower boundary.
4. The method for simulating a sedimentary geological action process according to claim 1, wherein the method for judging the vertical relative position according to the sedimentary lower boundary and the reference plane is as follows:
calculating the position difference of the deposition lower boundary in the vertical direction according to the vertical position coordinate of the deposition lower boundary and the reference plane to obtain a vertical relative position;
and when the vertical relative position is larger than 2-4 particle diameters, generating deposition particles in the updated stratum model, otherwise, not generating deposition particles in the updated stratum model.
5. A method of modeling a sedimentary geologic process according to claim 1, wherein the method of generating sedimentary particles in the updated earth model is:
acquiring the particle generation range of the updated stratum model according to the deposition lower boundary and the reference surface: (x) 1 ~x 2 ,y 1 ~y 2 ) Wherein x is 1 Representing the horizontal position coordinate, x, of the left boundary wall of the updated stratum model 2 Representing the horizontal position coordinate, y, of the right boundary wall of the updated stratum model 1 Vertical position coordinates representing the lower boundary of the deposit, y 2 Representing the sum of the vertical position coordinates of the reference surface and 2-4 particle diameters;
the deposited particles are generated by adopting a random filling and non-overlapping method in the particle generation range.
6. A sedimentary geologic process simulation method according to claim 1, wherein the sedimentary earth model is obtained by:
obtaining deposition calculation time steps according to the diameter and density of the deposited particles;
enabling the deposited particles to move under the action of dead weight, detecting the speed of the deposited particles, and returning the speed of the deposited particles to zero when the speed of the deposited particles is greater than a preset value;
when the movement time of the deposited particles reaches the step number of the deposition calculation time, the deposited particles are densely accumulated under the action of dead weight, and a deposited stratum model is obtained.
7. The method for simulating a sedimentary geologic process of claim 1, wherein the method for velocity reducing the sedimentary earth model comprises:
grouping new deposited particles and existing particles in the updated stratum model to obtain a new deposited particle group and an old deposited particle group;
and reducing the speed of each particle in the old sedimentary particle group according to the model speed linked list, and reducing the speed of a boundary wall in the sedimentary stratum model.
8. A sedimentary geologic process simulation device, comprising:
the stratum model initializing module is used for generating particles in a preset boundary wall, obtaining a stratum model with initial dead weight balance by adjusting the position of the boundary wall, and determining a reference surface of sedimentation in the stratum model;
the stratum model updating module is used for changing the positions of surface layer particles in the stratum model through extrusion loading and substrate sedimentation, and obtaining an updated stratum model according to the preset simulated motion time steps;
the linked list module is used for generating a model speed linked list at the current moment according to the updated stratum model; the model speed linked list comprises the speeds of all particles in the stratum model updated at the current moment, the speeds of the left boundary wall of the stratum model updated, the speeds of the right boundary wall of the stratum model updated and the speeds of the bottom boundary wall of the stratum model updated;
the particle screening module is used for acquiring a deposition lower boundary according to the updated stratum model;
the deposition module is used for judging the vertical relative position according to the lower deposition boundary and the reference plane, generating deposition particles in the updated stratum model when the judgment condition is met, and compacting the deposition particles under the action of dead weight by limiting the speed of the deposition particles to obtain a deposition stratum model;
and the speed reduction module is used for carrying out speed reduction on the sedimentary stratum model according to the model speed linked list so as to finish one-time sedimentary action.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110885025.4A CN113722967B (en) | 2021-08-03 | 2021-08-03 | Sedimentary geological action process simulation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110885025.4A CN113722967B (en) | 2021-08-03 | 2021-08-03 | Sedimentary geological action process simulation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113722967A CN113722967A (en) | 2021-11-30 |
| CN113722967B true CN113722967B (en) | 2024-03-15 |
Family
ID=78674719
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110885025.4A Active CN113722967B (en) | 2021-08-03 | 2021-08-03 | Sedimentary geological action process simulation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113722967B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204667766U (en) * | 2015-06-02 | 2015-09-23 | 云南国土资源职业学院 | Model is shown in a kind of deposition teaching |
| CN106202749A (en) * | 2016-07-15 | 2016-12-07 | 长江大学 | Deposition analogy method deposition process controlled based on datum level |
| CN106644564A (en) * | 2016-10-12 | 2017-05-10 | 中国地质大学（武汉） | Deep-sea sediment gravity coring test stand and operation method thereof |
| CN106959261A (en) * | 2017-02-28 | 2017-07-18 | 长江大学 | It is a kind of to predict sedimentary particle distribution and the method for proportioning |
| CN107015290A (en) * | 2017-03-13 | 2017-08-04 | 西北大学 | A kind of method that reworked garden basin primary deposit looks are recovered |
| CN111008472A (en) * | 2019-12-02 | 2020-04-14 | 东华理工大学 | Discrete element-based splitting basin extension process simulation method |
| CN111638552A (en) * | 2019-03-01 | 2020-09-08 | 中国石油化工股份有限公司 | Ancient landform restoration method |
-
2021
- 2021-08-03 CN CN202110885025.4A patent/CN113722967B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204667766U (en) * | 2015-06-02 | 2015-09-23 | 云南国土资源职业学院 | Model is shown in a kind of deposition teaching |
| CN106202749A (en) * | 2016-07-15 | 2016-12-07 | 长江大学 | Deposition analogy method deposition process controlled based on datum level |
| CN106644564A (en) * | 2016-10-12 | 2017-05-10 | 中国地质大学（武汉） | Deep-sea sediment gravity coring test stand and operation method thereof |
| CN106959261A (en) * | 2017-02-28 | 2017-07-18 | 长江大学 | It is a kind of to predict sedimentary particle distribution and the method for proportioning |
| CN107015290A (en) * | 2017-03-13 | 2017-08-04 | 西北大学 | A kind of method that reworked garden basin primary deposit looks are recovered |
| CN111638552A (en) * | 2019-03-01 | 2020-09-08 | 中国石油化工股份有限公司 | Ancient landform restoration method |
| CN111008472A (en) * | 2019-12-02 | 2020-04-14 | 东华理工大学 | Discrete element-based splitting basin extension process simulation method |
Non-Patent Citations (3)
| Title |
|---|
| "基于 LBM-IMB-DEM 分析岩土多孔介质的渗透淤堵特性";金 磊 等;《岩土工程学报》;全文 * |
| "基于距离和局部 Delaunay 三角化控制的 颗粒离散元模型填充方法研究 ";王秀菊 等;《岩土力学》;第36卷(第7期);2081-2087 * |
| 基于颗粒离散元法描述岩石断面形貌的方法研究;丛怡;王在泉;张黎明;丛宇;;防灾减灾工程学报(01);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113722967A (en) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105842735B (en) | Region rock mass microseism seismic source location method with complicated VELOCITY DISTRIBUTION | |
| CA2717889C (en) | Method for building a depositional space corresponding to a geological domain | |
| CN106202749B (en) | The deposition analogy method that deposition process is controlled based on datum level | |
| CN105277978B (en) | A kind of method and device for determining near-surface velocity model | |
| CN105319589B (en) | A kind of fully automatic stereo chromatography conversion method using local lineups slope | |
| EP3189357B1 (en) | Sediment transport simulation with parameterized templates for depth profiling | |
| CN105701319B (en) | A kind of microfacies modeling method under horizontal well rule well pattern | |
| CN109188519B (en) | System and method for inverting longitudinal and transverse wave speeds of elastic waves under polar coordinates | |
| CN107894618B (en) | A kind of full waveform inversion gradient preprocess method based on model smoothing algorithm | |
| CN111859748B (en) | Ocean internal wave simulation method based on vertical mixed coordinates | |
| CN108645994A (en) | A kind of geology stochastic inversion methods and device based on Multiple-Point Geostatistics | |
| CN109325250A (en) | A kind of method for numerical simulation and system of HIGH-SPEED LANDSLIDE-clast stream movement etching effect | |
| CN111006987A (en) | Reservoir fracture pore permeability parameter prediction method | |
| CN113722967B (en) | Sedimentary geological action process simulation method and device | |
| CN109376387A (en) | Evaluation method, device and the equipment of the closure of geological fault | |
| CN113221228B (en) | Hydropower station underground cave group surrounding rock mechanical parameter inversion method | |
| CN108680968B (en) | Evaluation method and device for seismic exploration data acquisition observation system in complex structural area | |
| CN113468466A (en) | Multi-working-condition one-dimensional wave equation solving method based on neural network | |
| CN117669411A (en) | Gravel packing sand prevention simulation method based on CFD-DEM | |
| CN110688728B (en) | Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment | |
| CN107507179A (en) | Rock And Soil quantitative analysis method based on GOCAD | |
| CN103236087B (en) | A kind of construction method of triangular prism geologic model | |
| Haeri et al. | Accurate and real-time simulation of rover wheel traction | |
| CN110032788A (en) | A kind of discrete element simulation method, equipment and the system of plate subduction deformation process | |
| CN112733242B (en) | Method for determining large deformation of side slope based on object point method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |