Mineral matrix deformation simulation method, terminal and readable storage medium

Abstract

The invention relates to the technical field of rock mass mechanics, and discloses a mineral matrix deformation simulation method, a terminal and a readable storage medium, the method comprises the following steps: obtaining mineral shape and size data and mineral material attribute data, and establishing a mineral particle accumulation body according to the mineral shape and size data; different stress conditions are set for the mineral particle accumulation body, and the change process of each mineral particle in the mineral particle accumulation body under the acting force is simulated according to the mineral material attribute data; and obtaining a simulation result according to the change process of the mineral particles under the acting force. According to the method, the mineral particle accumulation body consistent with the real mineral in shape and size is constructed, the physical and chemical effects among the mineral particles are considered according to the material attribute data of the mineral, and a more accurate simulation calculation result is obtained.

Description

Technical Field

[0001] The invention relates to the technical field of rock mass mechanics, and in particular to a simulation method, a terminal and a readable storage medium for mineral matrix deformation. Background Art

[0002] The deformation of mineral particles has a wide range of applications in many fields. For example, in the field of oil and gas extraction, it is necessary to simulate the behavior of rock fracture by simulating the deformation of mineral particles. In the study of geological disasters, it is necessary to simulate the possible consequences of rock fracture.

[0003] Traditional simulation methods are usually based on the contact model of spherical mineral particles. Although this method is computationally simple, it has significant limitations in practical applications. The matrix minerals in various backgrounds are diverse and have different properties. The traditional spherical mineral particle model cannot accurately describe the true polyhedral characteristics of mineral particles.

[0004] Secondly, traditional models also fail to fully consider the physical and chemical interactions between mineral particles, such as van der Waals forces and double-layer forces. These interactions have an important influence on the crushing, sliding, and reorganization behaviors between mineral particles at the microscopic scale, which in turn affects the evolution of the pore structure, porosity, and permeability of the matrix. Therefore, traditional simulation methods cannot accurately characterize the actual stress response of mineral particles, especially when multiple minerals coexist and under complex stress conditions, and lack sufficient accuracy and applicability. Summary of the invention

[0005] The purpose of the present invention is to provide a simulation method, terminal and readable storage medium for mineral matrix deformation, aiming to solve the problem that the existing spherical mineral particles cannot accurately characterize the actual stress response of the mineral particles and lack precision and applicability.

[0006] The technical solution adopted by the present invention to solve the technical problem is as follows: The present invention provides a method for simulating the deformation of a mineral matrix, and the method for simulating the deformation of a mineral matrix comprises: Acquire mineral shape and size data and mineral material property data, and establish a mineral particle accumulation body according to the mineral shape and size data; Setting different stress conditions for the mineral particle accumulation body, and simulating the change process of each mineral particle in the mineral particle accumulation body under the action of force according to the mineral material property data; According to the change process of the mineral particles under the force, the position results of the mineral particles and the parameter results of the mineral are obtained.

[0007] Furthermore, simulating the change process of each mineral particle in the mineral particle accumulation body under the action of the force according to the mineral material property data specifically includes: Obtaining collision parameters of each mineral particle according to the mineral particle accumulation body; For each collision of each mineral particle, calculating the force exerted on each mineral particle according to the collision parameters and the mineral material property data; For each mineral particle, the change process of the mineral particle under the force is simulated according to the force acting on the mineral particle.

[0008] Further, the mineral material property data includes Young's modulus and Poisson's ratio; The collision parameters include normal overlap distance and relative velocity at the contact point; The force includes a normal force; The calculating the force exerted on each mineral particle according to the collision parameter and the mineral material property data specifically includes: Calculate the equivalent radius of the collision. If two mineral particles are in contact, the equivalent radius is calculated as: ; in, represents the equivalent radius, Indicates the size of the first of the two mineral particles in contact, Indicates the size of the second of two mineral particles in contact; If the mineral grains are in contact with the boundary, the equivalent radius is calculated as: ; in, Indicates the size of mineral particles; Compute the effective Young's modulus for the collision: ; in, represents the effective Young's modulus. If it is the contact between two mineral particles, then represents the Young's modulus of the first mineral particle of the two contacting mineral particles, It represents the Young's modulus of the second mineral particle in the two contacting mineral particles. represents the Poisson's ratio of the first of the two contacting mineral particles, It represents the Poisson's ratio of the second mineral particle in the two contacting mineral particles. If the mineral particles are in contact with the boundary, then represents the Young's modulus of mineral particles, represents the Young's modulus of the boundary, represents the Poisson's ratio of mineral particles, represents the Poisson's ratio of the boundary; The stiffness coefficient of the mineral is calculated based on the effective Young's modulus and the equivalent radius: ; in, represents the stiffness coefficient; Calculate the effective volume of the collision. If two mineral particles are in contact, the effective mass is calculated as: ; in, represents the effective mass, represents the mass of the first of the two mineral particles in contact, It represents the mass of the second mineral particle among the two mineral particles in contact; If the mineral particles are in contact with the boundary, the effective mass is calculated as: ; in, Indicates the mass of mineral particles; The normal damping coefficient is calculated based on the stiffness coefficient and the effective mass: ; in, represents the normal damping coefficient, represents the damping ratio of the Hertz spring damper; The normal force on each mineral particle is calculated according to the normal overlap distance, the relative velocity at the contact point, the stiffness coefficient and the normal damping coefficient: ; in, represents the normal force, represents the normal overlap distance of the contact, Represents the normal component of the relative velocity at the contact point.

[0009] Furthermore, the force also includes a tangential force; The method further comprises: calculating the normal force on each mineral particle according to the normal overlap distance, the relative speed at the contact point, the stiffness coefficient and the normal damping coefficient; and then: The maximum relative tangential displacement is calculated according to the Poisson's ratio: ; in, represents the maximum relative tangential displacement, represents the friction coefficient; The linear Coulomb spring model parameters are calculated based on the maximum relative tangential displacement: ; in, are the linear Coulomb spring model parameters, represents the tangential relative displacement at the contact point; The tangential force is calculated according to the linear Coulomb spring model parameters, the maximum relative tangential displacement and the normal force: ; in, represents the tangential force, represents the tangential component of the relative velocity at the contact point, represents the tangential damping ratio.

[0010] Further, the force also includes adhesion force; The mineral material property data also includes surface energy; The method of calculating the force exerted on each mineral particle according to the collision parameter and the mineral material property data further includes: The adhesion force is calculated based on the surface energy: ; in, Represents adhesion, represents the surface energy, represents the contact radius; Modify the normal force to: .

[0011] Furthermore, the simulation of the change process of the mineral particles under the force according to the force applied to the mineral particles specifically includes: At each time step, the position and motion state of the mineral particles are updated according to the forces acting on the mineral particles; Calculate the instantaneous value of the specific impact energy received by each mineral particle according to the force received by the mineral particle; Whether the mineral particles are broken is determined based on the instantaneous value of the specific impact energy received by the mineral particles.

[0012] Furthermore, judging whether the mineral particles are broken according to the instantaneous value of the specific impact energy received specifically includes: According to the minimum specific energy value set Calculation of the minimum crushing energy of mineral particles : ; in, Indicates the reference size of the setting. Indicates the size of mineral particles; If the instantaneous value of the impact energy is Greater than minimum crushing energy And the instantaneous value of the impact energy is Greater than the maximum specific impact energy that the mineral particles have received in the previous event step , then the cumulative specific impact energy value of the mineral particles is updated : ; in, It represents the instantaneous value of specific impact energy in the previous time step; According to the cumulative specific impact energy value of mineral particles Calculating the impact strength of mineral particles : ; in, Indicates the material constants set; like If it is greater than the strength of the mineral particles, it is judged that the mineral particles are broken.

[0013] Further, the parameter results include porosity results and effective particle size results; The step of obtaining the position results of the mineral particles and the parameter results of the minerals according to the change process of the mineral particles under the force specifically includes: According to the change process of the mineral particles under the force, the position result of each mineral particle at each time step in the process is obtained; According to the change process of the mineral particles under the force, the porosity result and the effective particle size result of each time step are obtained.

[0014] In addition, to achieve the above-mentioned purpose, the present invention also provides a terminal, which includes: a memory, a processor, and a simulation program for mineral matrix deformation stored in the memory and executable on the processor, wherein the simulation program for mineral matrix deformation, when executed by the processor, controls the terminal to implement the steps of the simulation method for mineral matrix deformation as described above.

[0015] In addition, to achieve the above objectives, the present invention also provides a readable storage medium, which stores a simulation program for mineral matrix deformation. When the simulation program for mineral matrix deformation is executed by a processor, the steps of the simulation method for mineral matrix deformation as described above are implemented.

[0016] The present invention adopts the above technical solution to achieve the following effects: The present invention obtains the mineral shape and size data of the mineral and constructs a mineral particle accumulation body that is consistent with the actual mineral shape and size to replace the original spherical mineral model, thereby overcoming the problem that previous microscopic numerical simulations mainly focused on the interstitial flow of spherical matrix mineral particles. The present invention also obtains the material property data of the mineral and considers the physical and chemical interactions between the mineral particles based on the material property data of the mineral to obtain more accurate simulation calculation results. According to the changes in the simulated mineral particles, the deformation results of the mineral matrix, that is, the mineral accumulation body, are simulated. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 is a flow chart of steps of a method for simulating mineral matrix deformation in a preferred embodiment of the present invention; Figure 2 It is a framework flow chart of a method for simulating mineral matrix deformation in a preferred embodiment of the present invention; Figure 3 is a schematic diagram of an initial accumulation body in a preferred embodiment of the present invention; Figure 4 is a schematic diagram of a stable accumulation body in a preferred embodiment of the present invention; Figure 5 Schematic diagram of the change of effective particle sizes of three mineral particles in a preferred embodiment of the present invention; Figure 6 It is a schematic diagram of the change results of the porosity of three kinds of mineral particles in a preferred embodiment of the present invention; Figure 7 A schematic diagram of an operating environment of a preferred embodiment of a terminal of the present invention. DETAILED DESCRIPTION

[0018] In order to make the purpose, technical solution and advantages of the present invention clearer and more specific, the present invention is further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0019] Embodiment 1 Please refer to Figure 1 and Figure 2 Embodiment 1 of the present invention is a method for simulating deformation of a mineral matrix, comprising the steps of: S1. Obtain mineral shape and size data and mineral material property data, and establish a mineral particle accumulation body according to the mineral shape and size data.

[0020] Specifically, in this embodiment, the material properties of the mineral particles are adjusted according to the collected mineral material property data, mainly including the size, density, surface energy, Poisson's ratio and Young's modulus of the mineral particles.

[0021] In theory, shale minerals are diverse and have different properties. Due to different chemical compositions, the Hamaker constants between shale mineral particles are different, so there are also certain differences in the van der Waals force and double electron layer force between mineral particles, which is an important factor in analyzing the changes in mineral particles. In addition, mineral types have different crystal structures, mechanical properties and other characteristics, and their deformation and crushing mechanisms are also important research topics in the fields of mechanics and engineering.

[0022] In this embodiment, based on the investigation of typical mineral types of shale matrix, three mineral types with higher component contents in shales, namely kaolinite, quartz and calcite, were selected. In this embodiment, various mechanical parameters of these kaolinite, quartz and calcite, such as density, Poisson's ratio and particle size, were collected. The collected data are shown in Table 1.

[0023] Table 1: Mineral physical properties and mechanical parameters

[0024] After adjusting the properties of the mineral particles, import their geometric shape and size with reference to the crystal structure and electron microscope scan image, establish a representative mineral particle model for each mineral in the modeling software, reproduce the mineral particles in three dimensions, and then set the initial conditions, including setting the initial position and speed of the mineral particles, simulate the initial accumulation state in the area, and form an accumulation model.

[0025] In this embodiment, a hexahedral range is given as the boundary, and periodic boundary conditions are used on the four sides to generate mineral particles through continuous injection to form a stable mineral particle accumulation body, such as Figure 3 As shown, Figure 3 In FIG. 1 , the initial accumulation bodies of the three minerals of kaolinite, quartz and calcite of this embodiment are given.

[0026] S2. Setting different stress conditions for the mineral particle accumulation body, and simulating the change process of the mineral particles under the action of force according to the mineral material property data.

[0027] In this embodiment, on the basis of initialization, different stresses are applied to the accumulation body to simulate the change of overlying rock pressure, and for each time step, the changing state of the mineral particles is obtained, and the movement of the mineral particles is simulated by iteratively calculating the force and motion equations between each mineral particle. According to the mineral particles and different overlying rock pressure conditions, as the accumulation body reaches equilibrium again, the stress response and pore structure changes of the shale reservoir are observed, and the changes in the particle size distribution, effective particle size and porosity of the mineral particles after compressive deformation under different stress conditions are statistically analyzed.

[0028] Specifically, in this embodiment, a plate layer is added above the accumulation body, and different stresses of 20 MPa, 30 MPa, 40 MPa, 50 MPa and 60 MPa are applied thereto to simulate the change of overlying stratum pressure.

[0029] Afterwards, during the simulation, at each time step, the interaction forces between all mineral particles are calculated, and the motion state and position of each mineral particle are updated according to the forces.

[0030] Specifically, the interactions in this embodiment include normal force, tangential force and adhesion force, wherein the normal force is calculated using a hysteresis linear spring model, the tangential force is calculated using a linear Coulomb spring model, and the adhesion force is calculated using a Johnson-Kendall-Roberts adhesion model, and the interaction parameters are the model parameters of each model.

[0031] Specifically, in this embodiment, the normal force is calculated according to the following formula: ; in, represents the normal force in Newton. represents the stiffness coefficient, represents the normal damping coefficient, Represents the normal overlap distance of the contact in meters, Represents the normal component of the relative velocity at the contact point in meters per second.

[0032] Among them, the stiffness coefficient Specifically: ; in, represents the equivalent radius, represents the effective Young's modulus.

[0033] If two mineral particles are in contact, the equivalent radius is Calculated according to the following formula: ; in, Indicates the size of the first of the two mineral particles in contact, Represents the size of the second of two mineral particles in contact.

[0034] If the mineral grains are in contact with the boundary, the equivalent radius Calculated according to the following formula: ; in, Indicates the size of mineral particles, , and The unit is meter.

[0035] Effective Young's modulus Calculated by the following formula: ; Among them, if there is contact between two mineral particles, then represents the Young's modulus of the first mineral particle of the two contacting mineral particles, It represents the Young's modulus of the second mineral particle in the two contacting mineral particles. represents the Poisson's ratio of the first of the two contacting mineral particles, It represents the Poisson's ratio of the second mineral particle in the two contacting mineral particles. If the mineral particles are in contact with the boundary, then represents the Young's modulus of mineral particles, represents the Young's modulus of the boundary, represents the Poisson's ratio of mineral particles, represents the Poisson's ratio of the boundary.

[0036] Normal damping coefficient Specifically: ; in, represents the effective mass, Represents the damping ratio of the Hertz spring damper.

[0037] Among them, if two mineral particles are in contact, the effective mass Calculated according to the following formula: ; in, represents the mass of the first of the two mineral particles in contact, Represents the mass of the second mineral particle of two contacting mineral particles.

[0038] If the mineral particles are in contact with the boundary, the effective mass Calculated according to the following formula: ; in, Indicates the mass of mineral particles.

[0039] Damping ratio of Hertz spring damper Specifically: ; in, Represents the damping ratio of the linear spring damper.

[0040] The tangential force is calculated according to the following formula: ; in, represents the tangential force in Newton. represents the friction coefficient, dimensionless, represents the tangential relative displacement at the contact point in meters. represents the tangential component of the relative velocity at the contact point, in meters per second, represents the tangential damping ratio, It represents the maximum relative tangential displacement when the mineral particles begin to slide. are the linear Coulomb spring model parameters, specifically: ; Among them, the maximum relative tangential displacement when the particle starts to slide is Specifically: ; Tangential damping ratio Specifically: ; in, represents the restitution coefficient of the respective material interaction, Represents pi.

[0041] The adhesion force is calculated according to the following formula: ; in, Represents adhesion, represents the surface energy, Represents the contact radius between particles or between particles and the boundary, which is calculated by the normal overlap distance: ; The adhesion force has a corrective effect on the normal force. In this embodiment, the elastic part of the normal force Corrected by the sticky model: .

[0042] in, is the normal elastic value, that is, the normal force is corrected to: .

[0043] The impact received by each mineral particle is then calculated based on the force and movement, and the model is judged whether it has broken or fractured based on the impact. Based on the results of movement, fracture and fracture, the position and movement properties of each mineral particle in the accumulation model are updated.

[0044] Specifically, the impact force of the model is first obtained according to the movement of the mineral particles to determine whether the model is broken or ruptured. Specifically, this embodiment sets a reference minimum specific energy value. For each mineral particle, the minimum specific energy value can be Calculate the minimum crushing energy of the mineral particles : ; in, Indicates the reference dimension of the setting.

[0045] exist At this moment, the instantaneous value of the specific impact energy of the mineral particles is , then in the loading phase of this time step if and only if and When the cumulative impact energy value is updated ,in, is the maximum specific impact energy of the mineral particles in the previous event step, and the cumulative specific impact energy value Updated to: ; in, represents the time step length, express At that moment, it is the instantaneous value of the specific impact energy at the previous time step.

[0046] After completing the update of the cumulative specific impact energy value, the impact strength of the mineral particles is calculated according to the updated cumulative specific impact energy value. : ; in, Indicates the material constant set, when If the strength of the mineral particles is greater than that of the mineral particles, it is judged that the mineral particles are broken and there is a part smaller than Small broken mineral particles, such small broken mineral particles account for for: ; in, The maximum percentage of small broken mineral particles represents breakable material.

[0047] S3. According to the change process of the mineral particles under the force, the position results of the mineral particles and the parameter results of the mineral are obtained.

[0048] According to the mineral particles and different overlying rock pressure conditions, as the accumulation body reaches equilibrium again, the target part is selected to analyze the rotation speed and transformation speed of the mineral particles. When these two tend to 0, it can be judged that the accumulation body has reached stability.

[0049] During the process, simulation results are obtained and output for use in specific application scenarios. For example, in mining or disaster environments, there is a clearer understanding of the rock change process in order to set up corresponding plans. The simulation results are the data available in the simulation, including the three-dimensional simulation model of the mineral particles in the target area at each time step from the beginning to the stable state, which is used to intuitively display the spatial distribution of mineral particles; the key microscopic parameters of each mineral particle, including position coordinates, velocity vector, stress state, deformation degree, etc.; and the changing trend of the overall structural porosity, etc.

[0050] During the process, the overall deformation process can be observed macroscopically, the stress response and pore structure changes of the shale reservoir can be observed, and the changes in particle size distribution, effective particle size and porosity of mineral particles after compressive deformation under different stress conditions can be statistically analyzed.

[0051] According to the above simulation method, the microscopic crushing and deformation law of shale mineral particles can be simulated. In this embodiment, please refer to Figure 4 , Figure 4 This is a schematic diagram of the stable state of three types of rock mineral particles after compression deformation. At the initial moment, all rock mineral particles are large mineral particles. As the simulation process progresses, the pressure head presses down, the gaps between mineral particles gradually decrease, and the porosity of the mineral particle group decreases, causing the large mineral particle model to begin to break, thereby obtaining small mineral particles.

[0052] For microscopic parameter changes, please refer to Figure 5 and Figure 6 ,in, Figure 5 This is a schematic diagram of the change in the effective particle size of three mineral particles. The results show that with the increase in pressure, the effective particle size of the mineral particles shows a decreasing trend. This may be due to the increase in the contact area between the mineral particles, which leads to the enhanced interaction between the mineral particles, thereby causing the mineral particles to compress and deform. Kaolinite has a flaky structure, and its mineral particles may slip or break along the lamellae. Due to its structural characteristics, the reduction in effective particle size is relatively small. Quartz has the highest hardness among the three, has a stronger compressive resistance, a smaller degree of breakage, and a relatively small change in effective particle size. Calcite has a larger initial particle size than quartz, and has a lower Young's modulus, is more susceptible to breakage, and the change in effective particle size is more significant.

[0053] Figure 6This is a schematic diagram of the change in porosity of three mineral particles. The results show that with the increase of stress, the porosity of shale shows a downward trend. This discovery verifies the influence of stress sensitivity on shale porosity. The main reason is that the increase in pressure causes displacement and deformation between rock mineral particles, thereby reducing the pore space. This change is affected by mineral components. For example, quartz is a rigid mineral with strong resistance to compaction, while calcite, as a carbonate mineral, is formed during the evolution of shale after deposition and has weak resistance to compaction. Therefore, under the same pressure, the porosity of quartz mineral particles is higher than that of calcite mineral particles. Although kaolinite undergoes compression and rearrangement between mineral particles, due to its lamellar structure, the degree of porosity decrease is relatively small.

[0054] Embodiment 2 See also Figure 7 Based on the above method, the present invention further provides a terminal, which includes a processor 10, a memory 20 and a display 30. However, it should be understood that it is not required to implement all the components shown, and more or fewer components may be implemented instead.

[0055] In some embodiments, the memory 20 may be an internal storage unit of the terminal, such as a hard disk or memory of the terminal. In other embodiments, the memory 20 may also be an external storage device of the terminal, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the terminal. Further, the memory 20 may also include both an internal storage unit of the terminal and an external storage device. The memory 20 is used to store application software and various types of data installed in the terminal, such as the program code of the installation terminal. The memory 20 may also be used to temporarily store data that has been output or is to be output. In one embodiment, a simulation program 40 for mineral matrix deformation is stored on the memory 20, and the simulation program 40 for mineral matrix deformation can be executed by the processor 10, thereby realizing the terminal in the present application.

[0056] In some embodiments, the processor 10 may be a central processing unit (CPU), a microprocessor or other data processing chip, used to run the program code or process data stored in the memory 20, such as executing the relevant program of the simulation method of mineral matrix deformation.

[0057] In some embodiments, the display 30 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light Emitting Diode) touch device, etc. The display 30 is used to display information on the terminal and to display a visual user interface.

[0058] In one embodiment, when the processor 10 executes the simulation program 40 for mineral matrix deformation in the memory 20 , the steps of the method for simulating mineral matrix deformation as described above are implemented.

[0059] Embodiment 3 This embodiment provides a storage medium, wherein the readable storage medium stores a simulation program for mineral matrix deformation. When the simulation program for mineral matrix deformation is executed by a processor, the steps of the simulation method for mineral matrix deformation as described above are implemented.

[0060] In summary, the present invention obtains the mineral shape and size data of the mineral, constructs a mineral particle accumulation body consistent with the actual mineral shape and size to replace the original spherical mineral model, overcomes the problem that previous microscopic numerical simulations mainly focused on the interstitial flow of spherical matrix mineral particles, and obtains the material property data of the mineral, so as to consider the physical and chemical interactions between the mineral particles based on the material property data of the mineral, and obtain more accurate simulation calculation results.

[0061] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the presence of other identical elements in the process, method, article or terminal including the element.

[0062] Of course, those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing related hardware (such as a processor, a controller, etc.) through a computer program, and the program can be stored in a computer-readable storage medium, and the program can include the processes of the above-mentioned method embodiments when executed. The storage medium can be a memory, a disk, an optical disk, etc.

[0063] It should be understood that the application of the present invention is not limited to the above examples. For ordinary technicians in this field, improvements or changes can be made based on the above description. All these improvements and changes should fall within the scope of protection of the claims attached to the present invention.

Claims

1. A method for simulating deformation of a mineral matrix, characterized in that: The simulation method of the mineral matrix deformation comprises: Acquire mineral shape and size data and mineral material property data, and establish a mineral particle accumulation body according to the mineral shape and size data; Setting different stress conditions for the mineral particle accumulation body, and simulating the change process of each mineral particle in the mineral particle accumulation body under the action of force according to the mineral material property data; According to the change process of the mineral particles under the force, the position results of the mineral particles and the parameter results of the mineral are obtained.

2. A method for simulating deformation of a mineral matrix according to claim 1, characterized in that: The method further comprises: simulating the change process of each mineral particle in the mineral particle accumulation body under the action of the force according to the mineral material property data, specifically comprising: Obtaining collision parameters of each mineral particle according to the mineral particle accumulation body; For each collision of each mineral particle, calculating the force exerted on each mineral particle according to the collision parameters and the mineral material property data; For each mineral particle, the change process of the mineral particle under the force is simulated according to the force acting on the mineral particle.

3. A method for simulating deformation of a mineral matrix according to claim 2, characterized in that: The mineral material property data include Young's modulus and Poisson's ratio; The collision parameters include normal overlap distance and relative velocity at the contact point; The force includes a normal force; The calculating the force exerted on each mineral particle according to the collision parameter and the mineral material property data specifically includes: Calculate the equivalent radius of the collision. If two mineral particles are in contact, the equivalent radius is calculated as: ； in, represents the equivalent radius, Indicates the size of the first of the two mineral particles in contact, Indicates the size of the second of two mineral particles in contact; If the mineral grains are in contact with the boundary, the equivalent radius is calculated as: ； in, Indicates the size of mineral particles; Compute the effective Young's modulus for the collision: ； in, represents the effective Young's modulus. If it is the contact between two mineral particles, then represents the Young's modulus of the first mineral particle of the two contacting mineral particles, It represents the Young's modulus of the second mineral particle in the two contacting mineral particles. represents the Poisson's ratio of the first of the two contacting mineral particles, It represents the Poisson's ratio of the second mineral particle in the two contacting mineral particles. If the mineral particles are in contact with the boundary, then represents the Young's modulus of mineral particles, represents the Young's modulus of the boundary, represents the Poisson's ratio of mineral particles, represents the Poisson's ratio of the boundary; The stiffness coefficient of the mineral is calculated based on the effective Young's modulus and the equivalent radius: ； in, represents the stiffness coefficient; Calculate the effective volume of the collision. If two mineral particles are in contact, the effective mass is calculated as: ； in, represents the effective mass, represents the mass of the first of the two mineral particles in contact, It represents the mass of the second mineral particle among the two mineral particles in contact; If the mineral particles are in contact with the boundary, the effective mass is calculated as: ； in, Indicates the mass of mineral particles; The normal damping coefficient is calculated based on the stiffness coefficient and the effective mass: ； in, represents the normal damping coefficient, represents the damping ratio of the Hertz spring damper; The normal force on each mineral particle is calculated according to the normal overlap distance, the relative velocity at the contact point, the stiffness coefficient and the normal damping coefficient: ； in, represents the normal force, represents the normal overlap distance of the contact, Represents the normal component of the relative velocity at the contact point.

4. A method for simulating deformation of a mineral matrix according to claim 3, characterized in that: The force also includes a tangential force; The method further comprises: calculating the normal force on each mineral particle according to the normal overlap distance, the relative speed at the contact point, the stiffness coefficient and the normal damping coefficient; and then: The maximum relative tangential displacement is calculated according to the Poisson's ratio: ； in, represents the maximum relative tangential displacement, represents the friction coefficient; The linear Coulomb spring model parameters are calculated based on the maximum relative tangential displacement: ； in, are the linear Coulomb spring model parameters, represents the tangential relative displacement at the contact point; The tangential force is calculated according to the linear Coulomb spring model parameters, the maximum relative tangential displacement and the normal force: ； in, represents the tangential force, represents the tangential component of the relative velocity at the contact point, represents the tangential damping ratio.

5. A method for simulating deformation of a mineral matrix according to claim 3, characterized in that: The force also includes adhesion force; The mineral material property data also includes surface energy; The method of calculating the force exerted on each mineral particle according to the collision parameter and the mineral material property data further includes: The adhesion force is calculated based on the surface energy: ； in, Represents adhesion, represents the surface energy, represents the contact radius; Modify the normal force to: 。 6. A method for simulating mineral matrix deformation according to claim 2, characterized in that: The simulation of the change process of the mineral particles under the force according to the force applied to the mineral particles specifically includes: At each time step, the position and motion state of the mineral particles are updated according to the forces acting on the mineral particles; Calculate the instantaneous value of the specific impact energy received by each mineral particle according to the force received by the mineral particle; Whether the mineral particles are broken is determined based on the instantaneous value of the specific impact energy received by the mineral particles.

7. A method for simulating deformation of a mineral matrix according to claim 6, characterized in that: The method of judging whether the mineral particles are broken according to the instantaneous value of the specific impact energy received specifically includes: According to the minimum specific energy value set Calculation of the minimum crushing energy of mineral particles : ； in, Indicates the reference size of the setting. Indicates the size of mineral particles; If the instantaneous value of the impact energy is Greater than minimum crushing energy And the instantaneous value of the impact energy is Greater than the maximum specific impact energy that the mineral particles have received in the previous event step , then the cumulative specific impact energy value of the mineral particles is updated : ； in, It represents the instantaneous value of specific impact energy in the previous time step; According to the cumulative specific impact energy value of mineral particles Calculating the impact strength of mineral particles : ； in, Indicates the material constants set; like If it is greater than the strength of the mineral particles, it is judged that the mineral particles are broken.

8. A method for simulating deformation of a mineral matrix according to claim 1, characterized in that: The parameter results include porosity results and effective particle size results; The step of obtaining the position results of the mineral particles and the parameter results of the minerals according to the change process of the mineral particles under the force specifically includes: According to the change process of the mineral particles under the force, the position result of each mineral particle at each time step in the process is obtained; According to the change process of the mineral particles under the force, the porosity result and the effective particle size result of each time step are obtained.

9. A terminal, characterized in that: The terminal includes: a memory, a processor, and a simulation program for mineral matrix deformation stored in the memory and executable on the processor. When the simulation program for mineral matrix deformation is executed by the processor, the terminal is controlled to implement the steps of the method for simulating mineral matrix deformation as described in any one of claims 1 to 8.

10. A readable storage medium, characterized in that: The readable storage medium stores a simulation program for mineral matrix deformation, and when the simulation program for mineral matrix deformation is executed by a processor, the steps of the method for simulating mineral matrix deformation according to any one of claims 1 to 8 are implemented.

Images