Rock Grout Injection Simulation Method and System Based on Solid-Liquid Cohesion

Abstract

This invention provides a rock grout injection simulation method and system based on solid-liquid coupling, relating to the field of rock grout injection simulation technology. [Solution] The present invention includes obtaining the relationship between the time course of density change in the solidification process of grout used for grout injection and establishing an evolutionary model of grout density over time; dividing the solidification process of grout used for grout injection into multiple stages and obtaining the physical properties of each stage in the solidification process of grout used for grout injection; and adjusting the physical parameters of the grout at the relevant point in time during the rock body grout injection simulation based on the obtained evolutionary model of grout density over time and the physical properties of each stage in the solidification process of grout used for grout injection, and then performing the rock body grout injection simulation. The present invention can improve the accuracy of the simulation calculation process for rock and soil construction, such as rock body grout injection.

Description

[Technical Field]

[0001] (Cross-reference of related applications) This invention claims priority to a Chinese patent application filed with the China National Intellectual Property Administration on December 17, 2024, with application number 202411854086.4, and titled "Simulation Method and System for Rock Grout Injection Based on Solid-Liquid Bonding," the entirety of which is incorporated into this invention by reference, constitutes a part of this invention, and is used for all purposes.

[0002] The present invention relates to the technology of rock grout injection simulation, and more specifically to a rock grout injection simulation method and system based on solid-liquid bonding. [Background technology]

[0003] Numerical simulation studies of rockwork often involve complex liquid-solid doping geological conditions. Taking grout injection as an example, the rock mass during actual construction is not a perfect rock, but rather a rock mass containing a complex network of interwoven fissures. The presence of these fissures reduces the stability of the rock mass, so it is necessary to inject grout into the rock mass to fill the fissures and reinforce the rock mass.

[0004] In related technologies, it is common practice to calculate the two states of grout separately: before solidification and after solidification (liquid and solid). As a result, the calculation process is complicated because the two stages of liquid and solid are not combined. Furthermore, in stepwise calculations, it is difficult to determine where the liquid is converted to a solid substance, and it is also difficult to determine the diffusion range of the liquid. In contrast, conventional technologies propose a solution based on grout viscosity, determining whether the grout is liquid or solid based on its viscosity. However, in such approaches, errors due to frictional forces caused by rotor rotation in instruments such as rotational viscometers and rotational rheometers during the acquisition of data such as grout viscosity and shear rheology models are ignored, as are errors that affect the calculation results due to the entire process of material morphology change. [Overview of the project] [Problems that the invention aims to solve]

[0005] To address the shortcomings of conventional technology, the object of the present invention is to provide a rock grout injection simulation method and system based on solid-liquid bonding for improving the accuracy of simulation calculation processes for rock and soil construction work such as rock grout injection. [Means for solving the problem]

[0006] To achieve the above objective, according to several embodiments, according to a first aspect of the present invention, a rock grout injection simulation method based on solid-liquid bonding is provided, and the rock grout injection simulation method is To obtain the relationship between the change in density over time during the solidification process of grout used for grout injection, and to establish an evolutionary model of grout density over time, The solidification process of the grout used for grout injection is divided into multiple stages, and the physical properties of each stage in the solidification process of the grout used for grout injection are obtained. This includes adjusting the physical parameters of the grout at the relevant point in time during the rock grout injection simulation, based on the obtained grout density and time evolution model, and the physical properties of each stage in the solidification process of the grout used for grout injection, and then performing the rock grout injection simulation.

[0007] According to a second aspect of the present invention, a rock grout injection simulation system based on solid-liquid bonding is provided, and the rock grout injection simulation system is A model building module configured to acquire the relationship between the change in density over time during the solidification process of grout used for grout injection, and to establish an evolutionary model of grout density over time, A grout property acquisition module is configured to divide the grout solidification process used for grout injection into multiple stages and acquire the physical properties of each stage in the grout solidification process used for grout injection, The system includes a rock grout injection simulation module configured to adjust the physical parameters of the grout at a given point in time during the rock grout injection simulation, based on the obtained grout density and time evolution model, and the physical properties of each stage in the solidification process of the grout used for grout injection, and then perform the rock grout injection simulation. [Effects of the Invention]

[0008] Compared to the conventional technology, the beneficial effects of the present invention are as follows:

[0009] This invention provides a rock grout injection simulation method and system based on solid-liquid bonding. In rock grout injection simulation, the density of the grout is used as an indicator for determining the grout material morphology. By linking experiments and simulation calculations, the relationship between the change in density over time during the entire process of grout solidification from liquid to solid is obtained through grout transfer experiments. A model of grout density and time evolution is established, further guiding the changes in the physical parameters of the grout at different points in time during rock grout injection simulation. Simultaneously, the error in acquiring data on grout density is smaller compared to grout viscosity, and it is easier to measure when the grout morphology changes, thereby improving the accuracy of the rock grout injection simulation process.

[0010] The advantages of additional aspects of the present invention are given in part in the following description, become apparent from the following description, or are understood through the practice of the present invention.

[0011] The attached drawings of the specification constituting a part of the present invention are for providing a further understanding of the present invention, and the exemplary embodiments and their descriptions of the present invention are for interpreting the present invention and do not constitute an undue limitation of the present invention.

Brief Description of the Drawings

[0012] [Figure 1] It is a flowchart of the method in the first embodiment of the present invention. [Figure 2] It is a schematic diagram of a microscopic scan of the internal fracture network of fractured rock mass during grouting operation in the first embodiment of the present invention.

Modes for Carrying Out the Invention

[0013] Hereinafter, the present invention will be further described in conjunction with the attached drawings and examples.

[0014] Example 1 The first embodiment of the present invention provides a method for simulating rock mass grouting based on solid-liquid combination. As shown in FIGS. 1 to 2, it includes the following steps.

[0015] S1. Obtain the change relationship of the density with time during the solidification process of the grout used for grouting, and establish an evolution model of the density and time of the grout. S2. Divide the solidification process of the grout used for grouting into multiple stages, and respectively obtain the physical properties of each stage during the solidification process of the grout used for grouting. S3. Based on the obtained evolution model of the density and time of the grout and the physical properties of each stage during the solidification process of the grout used for grouting, adjust the physical parameters of the grout at the corresponding time during the simulation of rock mass grouting, and perform the simulation of rock mass grouting.

[0016] To simplify the description, the method according to this embodiment will be explained below using rock and soil construction grout injection work as an example. As shown in Figure 2, assuming that the grout 1 used is Grout A and that the mixing ratio of each material is constant, the rock and soil construction numerical calculation software adopted is Software B, which can perform numerical calculations and simulations of the rock body grout injection process. In some embodiments, Grout 1 is a cement grout commonly used in rock and soil construction grout injection work and can evolve into several other modified grouts, but their various performance characteristics must have features that can be measured in the laboratory; otherwise, it is difficult to apply them to numerical calculations. The fissure 4 is a complex network of fissure passages that are interwoven vertically and horizontally in the fissure rock body 3. Grout 1 is injected into the cracks 4 of the fissure rock body 3 through the grout injection hole 2. As it diffuses into the cracks 4, grout 1 gradually fills all the cracks 4 of the fissure rock body 3 until it stops flowing and solidifies completely, thereby reinforcing the fissure rock body 3.

[0017] In the grout solidification process, numerical calculations are typically used to accurately represent the diffusion, blockage, and reinforcement effects of the grout. However, the transformation of the grout from liquid to solid changes its physical and mechanical parameters, complicating the calculation process. Furthermore, calculations must be performed stepwise based on the pre-solidification and post-solidification states, making it difficult to integrate these two-step calculation processes.

[0018] In this embodiment, we propose linking experimental data and numerical calculations using grout density as an indicator for determining the grout material form. First, we obtain the change discipline in the grout solidification process through grout transfer tests in the laboratory, establish an evolutionary model of grout density and time, and obtain the parameters necessary for numerical calculations. Then, based on the evolutionary model of grout density and time, we determine the relevant parameters necessary for numerical calculations at different points in time and input them into software B for numerical calculations. By considering the entire process of grout A solidifying from liquid to solid in the calculation process, and combining the conversion process between liquid and solid grout A, we can simplify the calculation process, thereby determining the final position of liquid grout A flow, i.e., clarifying the final diffusion range of liquid grout A. By measuring density using the mass-volume method, we can avoid neglecting errors due to frictional force caused by the rotor rotation of instruments such as rotational viscometers or rotational rheometers during the acquisition of data such as grout viscosity and shear rheology models in grout viscosity-based solutions. Furthermore, instruments such as rotational viscometers and rotational rheometers, when measuring the viscosity and shear rheology of grout, have a rotating disc rotor or cylindrical rotor that can rotate within the grout. When this rotor rotates within the grout, a frictional force is generated between the rotor and the grout, and at the same time, the disturbance of the grout itself due to the rotor rotation may alter the progress of the hydration reaction within the grout. However, by determining the grout state based on density, measurement errors are not introduced, thereby improving the accuracy of the simulation calculations. In addition, in the calculation process of this embodiment, the entire process of the grout solidifying from liquid to solid is calculated, reducing errors that ignore the influence of the liquid-to-solidification process in conventional solutions on the calculation results.

[0019] In Step S1, a density measurement device based on the mass-volume method is added to the grout transfer test equipment. By collecting data on the change in density over time throughout the entire process of solidifying from liquid to solid during the grout transfer test, a model of the evolution of grout density over time is established.

[0020] A density measuring device based on the mass-volume method is installed in addition to the grout transfer passage in the grout transfer test equipment. This device calculates real-time instantaneous density data at the point where the measuring device is installed, based on the ratio of grout quality to pipe volume, with a calculation time interval of 1 second.

[0021] In grout transfer tests, density data is acquired in seconds, ultimately obtaining grout density data per second, i.e., data on the change in grout density over time during the test. Based on the obtained data, an evolutionary model of grout density over time is established. Statistical methods are employed to perform discipline and trend analysis on the obtained data, and a mathematical model that best matches the data is identified, found, and fitted.

[0022] An established model of grout density and its evolution over time is a mathematical model of the change in grout density over time. In this example, based on the results of the analysis of the data on the change in grout density over time, we assume that the data fits the mathematical model of the following equations. ρ = a1 + a2t + a3t 2 +...+a n t n-1 Here, ρ is the density of the grout (dimension: kilograms per cubic meter, kg / m³). 3 t represents time (dimensions are seconds, s), and the values ​​of t are a1, a2, a3, ..., a, with the test start time as the zero point. n This is a constant parameter for fitting.

[0023] Based on the obtained data on the time-series changes in grout density, the model is fitted and all constant parameters are solved. However, this model is not unique, and the optimal solution is one in which the obtained mathematical model can contain the most extensive range of time-series density data. If the amount of grout density change data per second included in the obtained mathematical model is relatively small, other mathematical model forms (e.g., exponential, logarithmic, or Fisher distribution functions) can be used for fitting until the optimal solution is found.

[0024] In step S2, the entire grout solidification process is divided into four time points based on the solidification and condensation characteristics of the grout.

[0025] (1) Hydration reaction initiation time t1: The initial stage of the solidification and condensation process of grout is fluid and plastic. The time point at which this process begins is denoted as the hydration reaction initiation time, and the start time of the laboratory test is used as the reference. (2) Initial solidification time t2: Grout gradually loses its fluidity and plasticity over time, and the initial solidification time obtained in the laboratory is used as the standard. (3) Final solidification time t3: After initial solidification, the cement grout gradually hardens until it evolves into a hard solid material, and its strength gradually increases. The final solidification time obtained in the laboratory is used as the standard. (4) Time t4 when standard curing is completed: That is, the time when the solidified grout body after final hardening is cured under standard conditions is recorded as t4.

[0026] Grout density at four time points in the laboratory: ρt 1、 ρt 2、 ρt 3、 Obtain ρt4 for each.

[0027] Based on the density change as a criterion, the solidification process of the grout used for grout injection can be divided into four stages based on the four points in time described above.

[0028] First stage: When the density of the grout is within the range of ρt1~ρt2, the grout in this stage is liquid and still possesses fluidity and plasticity. Second stage: When the grout density is within the range of ρt2~ρt3, the grout in this stage has already lost its fluidity and has acquired a certain strength. Third stage: When the grout density is within the range of ρt3 to ρt4, the strength of the solidified calculus gradually increases during this stage. Fourth stage: Within the test time after the grout density reached ρt4, the grout in this stage had already completely solidified, and the grout had already generated structural strength.

[0029] For the first stage, laboratory tests and studies are conducted on the grout, with data recording intervals of every second. Using the same method as the grout density and time evolution model obtained in step S1, mathematical models are obtained of the correlation between the fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout in its liquid state (liquid-related parameters obtained in the laboratory study, which can be adjusted according to the needs of different numerical simulation software) and time. For convenience of distinction, the fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout in its liquid state will hereafter be referred to as the first physical parameters.

[0030] For the second stage, laboratory tests are conducted on the grout to obtain the grout's strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters at time t2 (parameters obtained from laboratory studies, which can be adjusted according to the needs of different numerical simulation software). For convenience of distinction, the strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters of the grout in its solid state will hereafter be referred to as the second physical parameters.

[0031] For the third stage, laboratory tests are conducted on the grout to obtain the strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters (parameters obtained from laboratory studies, which can be adjusted according to the needs of different numerical simulation software) of the grout at time t3.

[0032] For the fourth stage, laboratory tests are conducted on the grout to obtain the strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters (parameters obtained from laboratory studies, which can be adjusted according to the requirements of different numerical simulation software) of the grout at time t4.

[0033] In step S3, the B software first inputs the grout injection pressure, the surrounding rock pressure in each direction, and the strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters of the fractured rock body that need to be simulated for grout injection, according to the construction requirements. Then, based on the grout density and time evolution model obtained in step S1 and the physical properties of each stage in the solidification process of the grout used for grout injection obtained in step S2, the physical parameters of the grout set in the grout injection simulation process are adjusted, and a rock body grout injection simulation is performed.

[0034] In this embodiment, the rock grout injection simulation is based on a density-time evolution model, and in step S2, a mathematical model is obtained of the correlation between the first physical parameter of the grout in its liquid (first stage) state and time. However, there is an error between the time in the numerical simulation software and the time during the test, and it is not possible to respond every second, so it is necessary to correct the time relationship between the two. The specific operation is as follows.

[0035] Using the time of change of grout fluidity as an example, we record the time at which the grout fluidity reaches the same value in both the test and the numerical simulation, and then record the relationship between these two times. If the grout fluidity is x (its dimension is determined based on the test instrument selected in the laboratory, but its numerical range does not change), then we denote the time at which the grout fluidity reaches x in the test as m, and the time at which the grout fluidity in the numerical simulation is x Let n be the time when this condition is reached. The m and n values ​​change correspondingly for each change in the grout flow rate x.

[0036] In this embodiment, based on the same method as establishing the grout density and time evolution model in step S1, a mathematical model of the change of grout fluidity x in the test with respect to time m within a range of values ​​is fitted and obtained, denoted as equation M, where it is assumed that the data fits the mathematical model of change of the following equation (i.e., equation M), x1 = b1 + b2m + b3m2 +... + b k m k-1 Here, x1 represents the grout fluidity obtained in the test (the dimension is generally mm, but it should be determined based on the measuring instrument actually adopted, and at the same time, attention should be paid to unit conversion during application). The value range of x1 is 0 ≤ x1 ≤ x. m represents time (the dimension is seconds, s), and the value of m starts from the test start time as the zero point and ends at the time when the fluidity reaches the x value. b1, b2, b3,..., b k are the fitting constant parameters, where k ∈ N + is as follows.

[0037] As can be understood by those skilled in the art, based on the obtained data of the change in grout fluidity over time, the model is fitted and all constant parameters are solved. However, the form of this model is not unique, and the optimal solution is that the obtained mathematical model can include the change data of grout fluidity over time in the largest range. When the change data of grout fluidity per second included in the obtained mathematical model is relatively small, other mathematical model forms (such as mathematical models of exponential functions, logarithmic functions, Fisher distribution functions, etc.) can be adopted for fitting until the optimal solution is found.

[0038] Similarly, based on the same method as establishing the evolution model of the density and time of the grout in step S1, the change mathematical model of the grout fluidity x with time and n within the value range in the numerical simulation is obtained by fitting and denoted as equation N. Here, assuming that the data conforms to the change mathematical model of the following equation (i.e., equation N), x2 = c1 + c2n + c3n x2 = c1 + c2n + c3n 2 +... + c k n k-1 Here, x² represents the grout fluidity during numerical simulation (the dimension is generally mm, but should be determined according to the simulation software actually used, and care should be taken with unit conversions during application), and its range is 0 ≤ x² ≤ x, and n represents time (the dimension is seconds, s), and the values ​​of n are c1, c2, c3, ..., c, with the calculation time of the start of the numerical simulation as the zero point and the endpoint of the values ​​when the fluidity reaches the value of x. k k∈N is a constant parameter of the fitting. + That is the case.

[0039] The x-values ​​reached by the grout flow rate in equations M and N are the same, i.e., x1 and x2 Since the range of values ​​is the same, if we let equation M = equation N, we can see that we can obtain the following equation. b1+b2m+b3m 2 +...+b k m k-1 = c1 + c2n + c3n 2 +...+c k n k-1

[0040] Since this equation has only two unknowns, m and n, we can determine the relationship between m and n, and further obtain a mathematical model of the correlation between fluidity and time after correcting for time errors. Time errors for other parameters (permeability, velocity, cohesive force, bulk modulus of the grout in liquid form; strength, stiffness, modulus, void modulus, Poisson's ratio, cohesive force, internal friction angle, and density of the grout in solid form) can be corrected in a similar manner.

[0041] The mathematical model of the time-dependent change in grout density obtained in step S1 is input into software B as density data for the entire grout injection simulation process after time error correction (the data is input into the software using a Matlab program or Fish programming method, and the same applies hereafter). During the grout injection simulation process, physical parameters necessary for the grout injection simulation are set based on the density change (at different stages).

[0042] (1) When the density value is in the first stage, that is, when the density is in the range ρt1~ρt2, the grout is a fluid liquid substance, and a mathematical model of the correlation between the fluidity, permeability, flow velocity, cohesive force, bulk modulus of the grout in its liquid state (after correcting for time errors) and time is input for the flow parameters of the grout. (2) When the density value is in the second stage, i.e., when the density is in the range ρt2~ρt3, the grout is initially solidified and studied based on the solid material, and the strength, stiffness, modulus of elasticity, void ratio, Poisson's ratio, cohesive force, and internal friction angle parameters of the grout at time t2 after correcting for time errors are input. (3) When the density value is in the third stage, i.e., the density is in the range ρt3~ρt4, the grout is finalized, and the strength, stiffness, modulus of elasticity, void ratio, Poisson's ratio, cohesive force, and internal friction angle parameters of the grout at time t3 after correcting for time errors are entered. (4) When the density value is in the fourth stage, i.e., after the density has reached ρt4, the grout is cured, and the strength, stiffness, modulus of elasticity, void modulus, Poisson's ratio, cohesive force, and internal friction angle parameters of the grout at time t4 after correcting for time errors are entered.

[0043] After parameter settings are complete, a grout injection simulation is performed using software. The simulation simulates the injection of grout into the fissure rock body through grout injection holes, causing the grout to flow through the fissures of the fissure rock body. Numerical calculations of rock and soil construction are then performed, and further calculations are completed to show how the material undergoes a solid-liquid morphology change during the numerical calculation.

[0044] After the entire grout injection simulation process is complete, the reliability of the method is verified by comparing the simulation results with the laboratory test results. For example, if the final deformation result of the grout injection simulation and the laboratory test are compared for the deformation amount of a rock body at a fixed location under study, and a relatively large comparison error is found, the model parameters are examined and adjusted until the error falls within an acceptable range. If a relatively small comparison error is found, it is indicated that this simulation method is accurate and reliable.

[0045] Example 2 This embodiment provides a rock grout injection simulation system based on solid-liquid bonding, and the rock grout injection simulation system is A model building module configured to acquire the relationship between the change in density over time during the solidification process of grout used for grout injection, and to establish an evolutionary model of grout density over time, A grout property acquisition module is configured to divide the grout solidification process used for grout injection into multiple stages and acquire the physical properties of each stage in the grout solidification process used for grout injection, The system includes a rock grout injection simulation module configured to adjust the physical parameters of the grout at a given point in time during the rock grout injection simulation, based on the obtained grout density and time evolution model, and the physical properties of each stage in the solidification process of the grout used for grout injection, and then perform the rock grout injection simulation.

[0046] Each module in this embodiment corresponds one-to-one with the steps of the method in Embodiment 1, and the specific implementation process is the same; therefore, a detailed explanation is omitted here.

[0047] The above description represents only preferred embodiments of the present invention and is not intended to limit it. Those skilled in the art will know that various modifications and changes are possible to the present invention. Within the spirit and principles of the present invention, any modifications, equivalent substitutions, or improvements must also be within the scope of protection. [Explanation of Symbols]

[0048] 1. Grout 2 Grout injection holes 3 Fissure rock body 4 fissure

Claims

1. A computer-based simulation method for rock grout injection based on solid-liquid bonding, The goal is to acquire data on the relationship between the change in density over time during the solidification process of grout used for grout injection (the entire process from liquid to solid) using density measuring equipment during grout transfer tests, and to establish an evolutionary model of grout density and time based on the acquired change relationship data. Based on the start of the hydration reaction of the grout used for grout injection, the initial solidification stage, the final solidification stage, and the completion of standard curing conditions, the solidification process of the grout used for grout injection is divided into four stages, and the physical properties of each stage in the solidification process of the grout used for grout injection are obtained. If the current density of the grout used for grout injection is between the density of the grout corresponding to the start of the hydration reaction and the density of the grout corresponding to the initial solidification point, then the grout used for grout injection is in the first stage. Within the first stage, the grout used for grout injection is in a liquid state, and at this time, the physical properties of the grout used for grout injection in the first stage are obtained based on a mathematical model of the correlation between the first physical parameter of the grout used for grout injection and time. The first physical parameters include the fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout used for grout injection when it is in liquid form, and the mathematical model of the correlation between the first physical parameters and time is obtained by fitting the relationship between the changes over time of the fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout obtained under laboratory test conditions. Based on a mathematical model of the correlation between the first physical parameter and time, and the corresponding correlation between the first physical parameter and time in numerical simulations, the relationship of change in the first physical parameter over test time and simulation time is obtained, respectively, and a time error correction is applied to the mathematical model of the correlation between the first physical parameter and time. Based on the obtained grout density and time evolution model, and the physical properties of each stage in the solidification process of the grout used for grout injection, the rock body grout injection simulation is performed by adjusting the physical parameters of the grout at the relevant point in time during the rock body grout injection simulation. A rock body grout injection simulation method characterized by: determining the current stage of the grout to be used for grout injection based on the current grout density output during the rock body grout injection simulation; setting physical parameters for the grout based on the current stage of the grout to be used for grout injection; and, if the grout to be used for grout injection is currently in the first stage, setting grout injection simulation parameters based on a mathematical model of the correlation between the obtained first physical parameters of the grout and time.

2. The aforementioned grout density and time evolution model was obtained by fitting using the following equation: ρ=a １ +a ２ t+a ３ t ２ +...+a ｎ t ｎ－１ Here, ρ represents the density of the grout, t represents time, a １ a ２ a ３ , ... a ｎ The rock grout injection simulation method based on solid-liquid bonding according to claim 1, characterized in that is a constant parameter of fitting.

3. If the current density of the grout used for grout injection lies between the density of the grout corresponding to the initial solidification point and the density of the grout corresponding to the final solidification point, the grout used for grout injection is in a second stage, where, within the second stage, the grout state is solid, and the physical properties of the grout in the second stage are determined based on the second physical parameters of the grout at the initial solidification point. and / or, If the current density of the grout used for grout injection is between the density of the grout corresponding to the final solidification point and the density of the grout corresponding to the completion of standard curing conditions, the grout used for grout injection is in a third stage, where the state of the grout used for grout injection is solid, and the physical properties of the grout used for grout injection in this third stage are determined based on the second physical parameter of the grout used for grout injection at the final solidification point. and / or, When the current density value of the grout used for grout injection reaches the density value of the grout corresponding to the time when standard curing is completed, the grout used for grout injection is in the fourth stage, where the state of the grout used for grout injection is solid, and the physical properties of the grout used for grout injection in the fourth stage are determined based on the second physical parameters at the time when the grout used for grout injection is completed under standard curing conditions. The rock grout injection simulation method based on solid-liquid bonding according to claim 1, characterized in that the second physical parameter includes the strength, stiffness, modulus of elasticity, void ratio, Poisson's ratio, cohesive force, and internal friction angle parameter of the grout used for grout injection when it is in solid form.

4. A rock grout injection simulation system based on solid-liquid bonding, A model building module configured to acquire the relationship between the change in density over time during the solidification process of grout used for grout injection, and to establish an evolutionary model of grout density over time, A grout property acquisition module is configured to acquire the physical properties of grout used for grout injection by dividing the grout solidification process into four stages based on the start of the hydration reaction of the grout used for grout injection, the initial solidification stage, the final solidification stage, and the completion of standard curing conditions, and based on the start of the hydration reaction of the grout used for grout injection, the initial solidification stage, the final solidification stage, and the completion of standard curing conditions, respectively. If the current density of the grout used for grout injection is between the density of the grout corresponding to the start of the hydration reaction and the density of the grout corresponding to the initial solidification point, then the grout used for grout injection is in the first stage. In the first stage, the grout used for grout injection is a liquid, and based on a mathematical model of the correlation between the first physical parameter of the grout used for grout injection and the test time, the physical properties of the grout used for grout injection in the first stage are obtained, and here, The first physical parameters include the fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout used for grout injection when it is in liquid form, and the mathematical model of the correlation between the first physical parameters and time is obtained by fitting the relationship between the changes in fluidity, permeability, flow velocity, cohesive force, and bulk modulus of the grout over time, obtained under laboratory test conditions. A grout property acquisition module that obtains the relationship of change in the first physical parameter over test time and simulation time, based on a mathematical model of the correlation between the first physical parameter and time, and the corresponding correlation with time in numerical simulations of the first physical parameter, and performs time error correction on the mathematical model of the correlation between the first physical parameter and time, A rock grout injection simulation module is configured to perform a rock grout injection simulation by adjusting the physical parameters of the grout at the relevant point in time during the rock grout injection simulation, based on the obtained grout density and time evolution model, and the physical properties of each stage in the solidification process of the grout used for grout injection. A rock grout injection simulation system characterized by including a rock grout injection simulation module that determines the current stage of the grout to be used for grout injection based on the current grout density output during the rock grout injection simulation, sets physical parameters for the grout based on the current stage of the grout to be used for grout injection, and, if the grout to be used for grout injection is currently in the first stage, sets grout injection simulation parameters based on a mathematical model of the correlation between the obtained first physical parameters of the grout and time.

Images