Rock grout injection simulation method and system based on solid-liquid coupling

Abstract

The present invention provides a method and system for simulating rock grout injection based on solid-liquid coupling, which relates to the technical field of rock grout injection simulation. The method includes: obtaining a time-varying relationship of density during the solidification process of grout used in grout injection, establishing a grout density-time evolution model; dividing the solidification process of grout used in grout injection into multiple stages and obtaining physical properties of each stage in the solidification process of grout used in grout injection; and adjusting physical parameters of the grout at corresponding times during the rock grout injection simulation based on the obtained grout density-time evolution model and the physical properties of each stage in the solidification process of grout used in grout injection, thereby performing the rock grout injection simulation. The present invention can improve the accuracy of the calculation process for simulating rock and earthworks such as rock grout injection.

Description

[Technical Field]

[0001] (CROSS-REFERENCE TO RELATED APPLICATIONS) This invention claims priority to a Chinese patent application filed with the State Intellectual Property Office of the People's Republic of China on December 17, 2024, bearing application number 202411854086.4 and entitled "Method and system for simulating rock grout injection based on solid-liquid bonding," the entire contents of which are incorporated by reference into and constitute part of the present invention and for all purposes.

[0002] The present invention relates to the technical field of rock grout injection simulation, and in particular to a rock grout injection simulation method and system based on solid-liquid coupling. [Background technology]

[0003] Numerical simulation studies of rock and earth engineering often encounter complex geological conditions involving liquid and solid doping. Take grouting as an example: the rock mass under construction is not a solid rock, but a complex network of crisscrossing fractures. The existence of these fractures reduces the stability of the rock mass, so grouting must be carried out to fill the fractures and reinforce the rock mass.

[0004] Related technologies typically calculate the two states of grout separately: liquid and solid before and after solidification. This disjointed calculation process makes it difficult to determine where the liquid transforms into a solid material and the extent of the liquid's diffusion. In response to this issue, prior art proposed a solution based on grout viscosity, determining whether the grout is liquid or solid based on the grout viscosity. However, this solution ignored errors caused by frictional forces due to rotor rotation in instruments such as rotational viscometers and rotational rheometers when acquiring data such as grout viscosity and shear rheology models, and also ignored errors that affect the calculation results due to the entire process of material transformation. Summary of the Invention [Problem to be solved by the invention]

[0005] In order to solve the deficiencies of the prior art, the present invention provides a rock grout injection simulation method and system based on solid-liquid coupling to improve the accuracy of the calculation process of rock and earthwork simulation such as rock grout injection. [Means for solving the problem]

[0006] To achieve the above object, according to some embodiments, a first aspect of the present invention provides a rock grout injection simulation method based on solid-liquid coupling, the rock grout injection simulation method comprising: To obtain the change relationship of the density of grout used in grout injection with time during its solidification process, and establish the evolution model of grout density and time; Dividing the solidification process of the grout used for grout injection into a plurality of stages, and acquiring the physical properties of each stage in the solidification process of the grout used for grout injection; and adjusting the physical parameters of the grout at a corresponding time during the rock mass 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 in the grout injection, and performing the rock mass grout injection simulation.

[0007] According to a second aspect of the present invention, there is provided a rock grout injection simulation system based on solid-liquid coupling, the rock grout injection simulation system comprising: a model construction module configured to obtain a change relationship of the density of the grout used for grout injection over time during its solidification process, and establish an evolution model of the density of the grout versus time; a grout property acquisition module configured to divide a solidification process of the grout used for grout injection into a plurality of stages and acquire physical properties of each stage in the solidification process of the grout used for grout injection; and a rock mass grout injection simulation module configured to adjust the physical parameters of the grout at a corresponding time during the rock mass 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 to perform the rock mass grout injection simulation. [Effects of the Invention]

[0008] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0009] The present invention provides a method and system for simulating rock grouting based on solid-liquid coupling. In the rock grouting simulation, grout density is used as an index for determining the grout material form, and experiments and simulation calculations are combined to obtain the time-dependent change relationship of grout density throughout the entire process of grout solidifying from liquid to solid through grout transfer experiments. A grout density vs. time evolution model is established to further guide the change of grout physical parameters at different times during the rock grouting simulation, thereby realizing the coupling calculation when the solid-liquid form changes during the simulation. At the same time, compared with grout viscosity, the data acquisition error of grout density is smaller and it is easier to measure when the grout form changes, thereby improving the accuracy of the rock grouting simulation process.

[0010] Additional aspect and advantages of the present invention will be set forth in part in the description that follows, and in part will be obvious from the description, or may be learned through practice of the invention.

[0011] The accompanying drawings of the specification, which form a part of the present invention, are intended to provide a further understanding of the present invention, and the illustrative embodiments of the present invention and their descriptions are intended to interpret the present invention and are not to be construed as unduly limiting the present invention. [Brief explanation of the drawings]

[0012] [Figure 1] 1 is a flowchart of a method according to a first embodiment of the present invention; [Figure 2] FIG. 1 is a schematic diagram of a microscan of the fracture network inside the fracture rock mass during grout injection work in Example 1 of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0013] The invention will be further explained below in connection with the accompanying drawings and examples.

[0014] Example 1 A first embodiment of the present invention provides a rock grout injection simulation method based on solid-liquid coupling, which includes the following steps, as shown in FIG. 1 and FIG. 2 .

[0015] S1, to obtain the change relationship of density over time during the solidification process of grout used in grout injection, and establish the evolution model of grout density and time. S2: Divide the solidification process of the grout used for grout injection into multiple stages, and obtain the physical properties of each stage in the solidification process of the grout used for grout injection. S3, 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 in grout injection, the physical parameters of the grout at the corresponding time during the rock mass grout injection simulation are adjusted, and the rock mass grout injection simulation is performed.

[0016] For simplicity, the method of this embodiment will be described below using a rock and earthwork grouting operation as an example. As shown in Figure 2, the grout 1 used is designated as grout A, and the blending ratio of each material is assumed to be constant. The adopted rock and earthwork numerical calculation software is software B, which can perform numerical calculations and simulations of the rock grouting process. In some embodiments, grout 1 is a cement grout commonly used in rock and earthwork grouting operations. It can be evolved into several other modified grouts, but their various properties must be measurable in a laboratory; otherwise, they are difficult to apply to numerical calculations. Fissures 4 are complex, interwoven fissure network passages present in the fissured rock mass 3. The grout 1 is injected into the fissures 4 of the fissured rock body 3 through the grout injection hole 2 and diffuses into the fissures 4, gradually filling all the fissures 4 of the fissured rock body 3 until the grout 1 no longer flows and completely solidifies, thereby reinforcing the fissured rock body 3.

[0017] In order to accurately represent the diffusion, blocking, and reinforcing effects of grout during the solidification process, numerical calculation methods are usually used to describe the grout injection process. However, as the grout transforms from liquid to solid, its physical and mechanical parameters change, making the numerical calculation process complicated. In addition, calculations can only be performed in stages based on the states before and after solidification, and the two-stage calculation process is difficult to combine.

[0018] In this example, we propose combining experimental data and numerical calculations using grout density as an index for determining grout material morphology. First, we obtain the change law during the grout solidification process through a grout transfer test in the laboratory, establish a grout density versus time evolution model, and obtain the parameters necessary for numerical calculation. Then, based on the grout density versus time evolution model, we determine the relevant parameters required for numerical calculation at different time points and input them into numerical calculation software B for calculation. By taking into account the entire process of grout A solidifying from liquid to solid, and combining the liquid-solid transformation process of grout A, we can determine the final position of liquid grout A flow, i.e., clarify the final diffusion range of liquid grout A. Measuring density using the mass-volume method allows us to use a solution based on grout viscosity, which avoids ignoring errors caused by frictional forces due to rotor rotation of instruments such as rotational viscometers or rotational rheometers when acquiring data such as grout viscosity and shear rheology models. Furthermore, when measuring the viscosity and shear rheology of grout, instruments such as rotational viscometers and rotational rheometers have a disc or cylindrical rotor that can rotate in the grout. When this rotor rotates in the grout, friction occurs between the rotor and the grout. At the same time, the rotor's rotation can disturb the grout itself, potentially altering the progress of the hydration reaction in the grout. However, determining the grout state based on density eliminates measurement errors, thereby improving the accuracy of the simulation calculation. Furthermore, the calculation process in this embodiment calculates the entire process of grout solidifying from liquid to solid, thereby reducing the error caused by ignoring the impact of the solidification process on the calculation results, as occurs in conventional solutions.

[0019] In step S1, a density measurement device based on the mass-volume method is additionally installed in the grout transfer test equipment to collect data on the change in density over time during the entire process of the grout solidifying from liquid to solid during the transfer test, thereby establishing a density and time evolution model of the grout.

[0020] A density measuring device based on the mass-volume method is installed in the grout transfer passage of the grout transfer test device, and this device calculates the instantaneous density data of the location where the measuring device is installed in real time based on the ratio of the grout quality and the piping volume of the location where the measuring device is installed, and the calculation time interval for each calculation is 1 second.

[0021] In the grout transfer test, density data is acquired in units of seconds, and the final data is the grout density data per second, i.e., the time-dependent change data of the grout density during the test. Based on the acquired data, a grout density vs. time evolution model is established. Statistical methods are used to perform discipline and trend analysis on the acquired data, and a mathematical model that is more consistent with the data is found and fitted.

[0022] The established grout density versus time evolution model is a mathematical model of the time-dependent change in grout density. In this embodiment, based on the analysis results of the time-dependent change data of grout density, it is assumed that the data fits the mathematical model of the following equation: ρ=a1+a2t+a3t 2 +...+a n t n-1 where ρ is the density of the grout (dimensions are kilograms per cubic meter, kg / m 3 t represents time (the dimension is seconds, s), and the values ​​of t are a1, a2, a3, ..., a n is a constant parameter of the fitting.

[0023] Based on the obtained data on the change in grout density over time, a model is fitted and all constant parameters are solved. However, this model is not unique, and the optimal solution is determined when the mathematical model can include the widest range of data on the change in density over time. If the data on the change in grout density per second included in the mathematical model is relatively small, other mathematical model forms (such as mathematical models of exponential function, logarithmic function, Fisher distribution function, etc.) can be adopted and fitted until the optimal solution is found.

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

[0025] (1) Hydration reaction start time t1: At the beginning of the solidification and condensation process of grout, it has fluidity and plasticity. The time starting point of this process is recorded as the hydration reaction start time, and is based on the start time of the laboratory test. (2) Initial setting time t2: Grout gradually loses fluidity and plasticity over time, and the initial setting time obtained in the laboratory is used as the basis. (3) Final setting time t3: After initial setting, the cement grout gradually hardens until it evolves into a hard solid material, and its strength gradually increases. This is based on the final setting time obtained in the laboratory. (4) Time t4 when standard condition curing is completed: In other words, the solidified grout stone body after final solidification is cured under standard conditions, and the time when the curing is completed is recorded as t4.

[0026] Density of grout at four points in the laboratory: ρt 1、 ρt 2、 ρt 3、 ρt4 is obtained respectively.

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

[0028] First stage: If the density of the grout is within the range ρt1-ρt2, the grout in this stage is liquid and still has fluidity and plasticity. Second stage: When the density of the grout is within the range of ρt2 to ρt3, the grout in this stage has already lost its fluidity and has a certain strength. Third stage: When the density of grout is within the range of ρt3~ρt4, the strength of the solidified grout stone body in this stage will gradually improve. Fourth stage: During the test time after the density of the grout reaches ρt4, the grout in this stage has already completely solidified and the grout has already developed structural strength.

[0029] In the first stage, laboratory testing is performed 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 for the correlation between time and the grout's liquid fluidity, permeability, flow rate, cohesion, and bulk modulus (liquid-related parameters obtained in the laboratory testing, which can be adjusted according to the needs of different numerical simulation software). For ease of distinction, hereinafter, the liquid fluidity, permeability, flow rate, cohesion, and bulk modulus of the grout are referred to as the first physical parameters.

[0030] In the second stage, laboratory testing is conducted on the grout to obtain the grout's strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters at time t2 (parameters obtained through laboratory testing, which can be adjusted according to the needs of different numerical simulation software). For ease of distinction, hereinafter, the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout when solid are referred to as the second physical parameters.

[0031] In the third stage, laboratory testing is carried out on the grout to obtain the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout at time t3 (the parameters obtained in the laboratory testing can be adjusted according to the needs of different numerical simulation software).

[0032] In the fourth stage, laboratory testing is carried out on the grout to obtain the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout at time t4 (the parameters obtained in the laboratory testing can be adjusted according to the needs of different numerical simulation software).

[0033] In step S3, first, in software B, the grout injection pressure, surrounding rock pressure in each direction, and the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the fractured rock mass for which grout injection simulation is required are input according to the construction needs. 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 then the rock mass grout injection simulation is performed.

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

[0035] Taking the time of change of grout fluidity as an example, the time when the grout fluidity in the test and the numerical simulation reaches the same value is recorded, and then the relationship between the two times is recorded. 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), the time when the grout fluidity in the test reaches x is written as m, and the time when the grout fluidity in the numerical simulation reaches x is written as m. The time when this is reached is denoted as n, and with each change in grout fluidity x, the m and n values ​​change correspondingly.

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

[0037] As will be understood by those skilled in the art, a model is fitted based on the obtained data on the change in grout fluidity over time, and all constant parameters are solved. However, this model format is not unique, and the optimal solution is determined when the obtained mathematical model can include the widest range of data on the change in grout fluidity over time. If the obtained mathematical model contains relatively little data on the change in grout fluidity per second, other mathematical model formats (e.g., mathematical models such as exponential function, logarithmic function, Fisher distribution function, etc.) can be adopted and fitted until an optimal solution is found.

[0038] Similarly, based on the same method as the establishment of the grout density and time evolution model in step S1, the grout fluidity x in the numerical simulation is determined within the range of time and n. is obtained by fitting a mathematical model of change of the following equation (i.e., equation N), where the data is assumed to fit the mathematical model of change of the following equation (i.e., equation N): x2=c1+c2n+c3n 2 +...+c k n k-1 Here, x2 represents the grout fluidity during the numerical simulation (the dimension is generally mm, but should be determined according to the actual simulation software configuration, and attention should be paid to unit conversion during application), and its value range is 0≦x2≦x, and n represents time (the dimension is seconds, s). The value of n is set to zero when the numerical simulation starts, and the end point is when the fluidity reaches the x value, and then c1, c2, c3, ..., c k is a constant parameter of the fitting, k∈N + is.

[0039] The x value at which the grout fluidity reaches in equation M and equation N is the same, i.e., x1 and x2 Since the range of values ​​of M and N is the same, if we set equation M = equation N, 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, i.e., the relationship between m and n can be determined, and a mathematical model of the correlation between fluidity and time after time error correction can be further obtained. The time errors of other parameters (liquid grout permeability, flow velocity, cohesion, bulk modulus, solid grout strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, internal friction angle, and grout density) can be corrected in the same way.

[0041] The mathematical model of the time-dependent change in density of the grout obtained in step S1 is input into software B as the density data for the entire grout injection simulation process after time error correction, and the simulation is performed (the data is input into the software using the Matlab program and the Fish language programming method, and the same applies below). During the grout injection simulation process, the density change (at different stages) is used as the judgment criterion, and the physical parameters required for the grout injection simulation are respectively set.

[0042] (1) When the density value is in the first stage, that is, when the density is in the range ρt1 to ρt2, the grout is a liquid material with fluidity. For the flow parameters of the grout, the obtained time error corrected mathematical model of the correlation between the liquid fluidity, permeability, flow velocity, cohesion, bulk modulus and time of the grout is input. (2) When the density value is in the second stage, that is, the density is in the range ρt2~ρt3, the grout is initially solidified and studied based on the solid material. The strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout obtained at time t2 after time error correction are input. (3) If the density value is in the third stage, that is, if the density is in the range ρt3 to ρt4, the grout is finally solidified, and the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout obtained at time t3 after time error correction are entered. (4) When the density value is in the fourth stage, that is, after the density reaches ρt4, the grout has been cured, and the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout obtained at time t4 after time error correction are entered.

[0043] After the parameter settings are complete, the software performs a grout injection simulation, and the grout is injected into the fractured rock body through the grout injection hole through the simulation, so that the grout flows through the fractures in the fractured rock body, and the rock and earthwork numerical calculation is performed, and the coupling calculation of the solid-liquid transformation of the material during the numerical calculation is completed.

[0044] After the entire grout injection simulation process is completed, the simulation results are compared with the laboratory test results to verify the reliability of the method. For example, for the deformation of a rock mass at a fixed position under study, the final deformation results of the grout injection simulation and the laboratory test are compared, and if it is found that the comparative error is relatively large, the model parameters are checked and adjusted until they reach the error tolerance range, and if it is found that the comparative error is relatively small, it indicates that the simulation method is accurate and reliable.

[0045] Example 2 This embodiment provides a rock grout injection simulation system based on solid-liquid coupling, and the rock grout injection simulation system includes: a model construction module configured to obtain a change relationship of the density of the grout used for grout injection over time during its solidification process, and establish an evolution model of the density of the grout versus time; a grout property acquisition module configured to divide a solidification process of the grout used for grout injection into a plurality of stages and acquire physical properties of each stage in the solidification process of the grout used for grout injection; and a rock mass grout injection simulation module configured to adjust the physical parameters of the grout at a corresponding time during the rock mass 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 to perform the rock mass grout injection simulation.

[0046] It should be noted that each module in this embodiment corresponds one-to-one to the step of the method in the first embodiment, and the specific implementation process is the same, so the description is omitted here.

[0047] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art can make various modifications and changes to the present invention. Any modifications, equivalent replacements, improvements, etc. within the spirit and principle of the present invention shall be included in the protection scope of the present invention. [Explanation of symbols]

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

Claims

1. 1. A computer-implemented method for rock mass grouting simulation based on solid-liquid coupling, comprising: Obtaining change relationship data of the density change over time during the solidification process (the entire process of solidifying from liquid to solid) of the grout used for grout injection by the density measuring device during the grout transfer test, and establishing a grout density and time evolution model based on the obtained change relationship data; Based on the start time of the hydration reaction of the grout used for grout injection, the initial solidification time, the final solidification time, and the time when standard condition curing is completed, 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 value of the grout used for grout injection is between the density value of the grout corresponding to the start of hydration reaction and the density value of the grout corresponding to the initial solidification time, the grout used for grout injection is in the first stage; In the first stage, the state of the grout used for grout injection is liquid, and at this time, the physical properties of the first stage of the grout used for grout injection are obtained based on a mathematical model of the correlation between a first physical parameter of the grout used for grout injection and time, the first physical parameters include the fluidity, permeability, flow rate, cohesion, and bulk modulus of the grout used for the grout injection when it is liquid, and the mathematical model of the correlation between the first physical parameters and time is obtained by fitting the change relationships of the fluidity, permeability, flow rate, cohesion, and bulk modulus of the grout obtained under laboratory test conditions, respectively; Based on the mathematical model of the correlation between the first physical parameter and time and the correlation between the corresponding first physical parameter and time in the numerical simulation, obtain a change relationship between the first physical parameter and the test time and a change relationship between the first physical parameter and the simulation time, respectively, and perform time error correction on the mathematical model of the correlation between the first physical parameter and time; Performing a rock mass grout injection simulation by adjusting each physical parameter of the grout at a corresponding time during the rock mass 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 current stage of the grout used for the grout injection is determined based on the current density of the grout output during the rock mass grout injection simulation; and a physical parameter of the grout is set based on the current stage of the grout used for the grout injection, wherein if the grout used for the grout injection is currently in a first stage, a grout injection simulation parameter is set based on a mathematical model of the correlation between the obtained first physical parameter of the grout and time.

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

3. When the current density value of the grout used for grout injection is between the density value of the grout corresponding to the initial solidification time and the density value of the grout corresponding to the final solidification time, the grout used for grout injection is in a second stage, wherein in the second stage, the state of the grout used for grout injection is solid, and the physical properties of the second stage of the grout used for grout injection are determined based on the second physical parameters of the grout used for grout injection at the initial solidification time; and / or When the current density value of the grout used for grout injection is between the density value of the grout corresponding to the final solidification time and the density of the grout corresponding to the time when standard condition curing is completed, the grout used for grout injection is in a third stage, wherein in the third stage, the state of the grout used for grout injection is solid, and the physical properties of the third stage of the grout used for grout injection are determined based on a second physical parameter of the grout used for grout injection at the final solidification time; and / or When the current density value of the grout used for grout injection reaches a density value of the grout corresponding to the time when standard condition curing is completed, the grout used for grout injection is in a fourth stage, wherein the state of the grout used for grout injection is solid within the fourth stage, and the physical properties of the grout used for grout injection in the fourth stage are determined based on a second physical parameter of the grout used for grout injection at the time when standard condition curing is completed; The method for simulating rock grout injection based on solid-liquid bonding according to claim 1, characterized in that the second physical parameters include the strength, stiffness, elastic modulus, void fraction, Poisson's ratio, cohesion, and internal friction angle parameters of the grout used in the grout injection when solid.

4. A rock mass grouting simulation system based on solid-liquid coupling, comprising: a model construction module configured to obtain a change relationship of the density of the grout used for grout injection over time during its solidification process, and establish an evolution model of the density of the grout versus time; A grout property acquisition module configured to divide the solidification process of the grout used for grout injection into four stages based on the start time of the hydration reaction of the grout used for grout injection, the initial solidification time, the final solidification time, and the time when standard condition curing is completed, and to acquire physical properties of each stage in the solidification process of the grout used for grout injection, If the current density value of the grout used for grout injection is between the density value of the grout corresponding to the start of hydration reaction and the density value of the grout corresponding to the initial solidification time, the grout used for grout injection is in the first stage; In the first stage, the grout used for the grout injection is liquid, and the physical properties of the first stage of the grout used for the grout injection are obtained based on a mathematical model of the correlation between a first physical parameter of the grout used for the grout injection and a test time, wherein: the first physical parameters include the fluidity, permeability, flow rate, cohesion, and bulk modulus of the grout used for the grout injection when it is in a liquid state, and the mathematical model of the correlation between the first physical parameters and time is obtained by fitting the change relationships of the fluidity, permeability, flow rate, cohesion, and bulk modulus of the grout obtained under laboratory test conditions with test time, respectively; a grout property acquisition module that acquires, based on a mathematical model of the correlation between the first physical parameter and time and the correlation between the corresponding first physical parameter and time in the numerical simulation, a change relationship between the first physical parameter and the test time and a change relationship between the first physical parameter and the simulation time, respectively, and performs time error correction on the mathematical model of the correlation between the first physical parameter and time; a rock mass grout injection simulation module configured to perform a rock mass grout injection simulation by adjusting each physical parameter of the grout at a corresponding time during the rock mass 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 mass grout injection simulation module that determines a current stage of the grout to be used for the grout injection based on a current grout density output during the rock mass grout injection simulation, and sets physical parameters of the grout based on the current stage of the grout to be used for the grout injection, and when the grout to be used for the grout injection is currently in a 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