Method and system for determining the range of columnar infiltration grouting considering the tortuosity of porous media

Abstract

The present invention relates to the field of ecological restoration technology, and in particular to a method and system for determining the range of columnar infiltration grouting considering the tortuosity of porous media, the method comprising the following steps: collecting porous media characteristic data and groundwater pressure at a selected grouting point, collecting Bingham cement slurry characteristic data and setting grouting parameters at the same time; considering the tortuosity of the porous media, and then determining the seepage velocity of Bingham cement slurry based on the rheological equation of Bingham fluid and the mechanical analysis equation of microelement in the grouting pipe; constructing a columnar infiltration grouting range solution model based on the seepage velocity and the calculation formula of the grouting amount per unit time under the columnar diffusion boundary condition, and solving the quantitative index of the spatial diffusion range of the columnar infiltration grouting according to the grouting parameters to determine the columnar infiltration grouting range. The present invention can improve the quality and reliability of grouting in ecological restoration projects, meet the construction needs under complex ecological engineering environments, and achieve effective restoration and sustainable development of ecosystems.

Description

Technical Field

[0001] The invention relates to the technical field of ecological restoration, and in particular to a method and a system for determining a columnar infiltration grouting range taking into account the tortuosity of a porous medium. Background Art

[0002] At present, in the process of columnar infiltration grouting, the determination of the grouting range often does not fully consider the tortuosity of porous media. Porous media are widely present in ecosystems, such as soil and rock, and their tortuosity has a significant impact on the flow path and diffusion range of slurry. The porous media structure in actual ecological restoration scenarios is complex, and the flow of slurry in them is not a simple linear motion, but is restricted and guided by tortuous channels. Ignoring tortuosity may lead to inaccurate estimation of the grouting range, which in turn affects the effect and sustainability of ecological restoration projects. At the same time, most of the existing grouting range determination methods do not comprehensively consider the influence of Bingham cement slurry's deadweight and grouting angle. As a common grouting material, Bingham cement slurry has different effects on the diffusion of slurry at different grouting angles in different ecological restoration scenarios. For example, in mountain ecological restoration, it may face inclined terrain. At this time, the deadweight of the slurry will work together with the grouting pressure to change the flow direction and diffusion speed of the slurry. However, the existing methods fail to fully integrate these factors, which makes the determination of grouting range in ecological restoration projects have certain limitations. In addition, the geological conditions faced by ecological restoration projects are complex and diverse, and may contain strata with different porosities and permeabilities, and there may also be groundwater flow and interaction with the surrounding ecological environment. In this context, a more accurate and comprehensive method for determining the range of columnar infiltration grouting is urgently needed to improve the quality and reliability of grouting in ecological restoration projects, meet the construction needs in complex ecological engineering environments, and achieve effective restoration and sustainable development of ecosystems. Summary of the invention

[0003] In view of the defects in the prior art, the present invention provides a method and system for determining the range of columnar infiltration grouting taking into account the tortuosity of porous media.

[0004] In order to achieve the above-mentioned purpose, in a first aspect, the present invention provides a method for determining the range of columnar penetration grouting taking into account the tortuosity of porous media, the method comprising the following steps: collecting porous media characteristic data and groundwater pressure at selected grouting point locations, and collecting Bingham cement slurry characteristic data at the same time; setting grouting parameters, the grouting parameters including grouting pressure, grouting time and grouting pipe laying angle; considering the tortuosity of the porous medium, and then determining the seepage velocity of the Bingham cement slurry based on the rheological equation of the Bingham fluid and the mechanical analysis equation of the microelement in the grouting pipe; constructing a columnar penetration grouting range solution model based on the seepage velocity and the calculation formula of the grouting amount per unit time under columnar diffusion boundary conditions; according to the grouting parameters, using the columnar penetration grouting range solution model to solve the quantitative index of the spatial diffusion range of the columnar penetration grouting, so as to determine the columnar penetration grouting range.

[0005] Optionally, the method of considering the tortuosity of the porous medium and determining the seepage velocity of the Bingham cement slurry based on the rheological equation of the Bingham fluid and the mechanical analysis equation of the microelement in the grouting pipe comprises the following steps:

[0006] Performing mechanical analysis on the microelement in the grouting pipe to obtain the mechanical analysis equation, and then combining the rheological equation to obtain the shear rate equation of Bingham cement slurry;

[0007] Integrating the shear rate equation to obtain a Bingham cement slurry flow rate, and calculating the Bingham cement slurry flow rate in the grouting pipe according to the Bingham cement slurry flow rate;

[0008] Calculating the average flow rate of the Bingham cement slurry in the grouting pipe according to the Bingham cement slurry flow rate, and then obtaining the reference seepage velocity of the Bingham cement slurry according to the Dupuit-Forchheimer formula;

[0009] After grouting, the reference seepage velocity is corrected by taking into account the tortuosity of the porous medium to obtain the seepage velocity.

[0010] Optionally, the shear rate equation satisfies the following relationship:

[0011] ,

[0012] in, is the shear rate, is the flow rate of Bingham cement slurry near a certain point in the Bingham cement slurry As the radius r of the microelement in the grouting pipe changes, is the plastic viscosity of Bingham cement slurry, is the yield shear stress of Bingham cement slurry, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

[0013] Optionally, the Bingham cement slurry flow rate and the Bingham cement slurry flow rate respectively satisfy the following relationship:

[0014] ,

[0015] ,

[0016] in, is the Bingham cement slurry flow rate, r is the radius of the microelement in the grouting pipe, is the plastic viscosity of Bingham cement slurry, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the angle between the Bingham cement slurry flow direction and the grouting pipe axis, is the yield shear stress of Bingham cement slurry, is the core radius of Bingham cement slurry, is the radius of the grouting pipe, and Q is the Bingham cement slurry flow rate.

[0017] Optionally, the average flow rate of the Bingham cement slurry satisfies the following relationship,

[0018] ,

[0019] in, is the average flow velocity of the Bingham cement slurry, Q is the flow rate of the Bingham cement slurry, is the radius of the grouting pipe.

[0020] Optionally, the seepage velocity satisfies the following relationship:

[0021] ,

[0022] Wherein, V is the seepage velocity, is the porous medium characteristic data, is the radius of the grouting pipe, is the plastic viscosity of Bingham cement slurry, is the tortuosity of the porous medium, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the yield shear stress of Bingham cement slurry, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

[0023] Optionally, the step of constructing a columnar permeation grouting range solution model based on the seepage velocity and the grouting volume per unit time calculation formula under the columnar diffusion boundary condition comprises the following steps:

[0024] Establishing a pressure gradient equation according to the seepage velocity and the calculation formula of the grouting volume per unit time;

[0025] The pressure gradient equation is integrated and the columnar diffusion boundary condition is considered to obtain the columnar penetration grouting range solution model.

[0026] Optionally, the pressure gradient equation satisfies the following relationship:

[0027] ,

[0028] in, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the tortuosity of the porous medium, is the plastic viscosity of Bingham cement slurry, q is the grouting volume per unit time, is the actual flow path length of Bingham cement slurry in the direction perpendicular to the horizontal plane, h is the columnar diffusion height, is the radius of the grouting pipe, is the porous medium characteristic data, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the yield shear stress of Bingham cement slurry, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

[0029] Optionally, the columnar penetration grouting range solution model satisfies the following relationship:

[0030] ,

[0031] in, is the grouting pressure, is the groundwater pressure, is the plastic viscosity of Bingham cement slurry, is the radius of the grouting pipe, is the quantitative index of the spatial diffusion range of the columnar penetration grouting, t is the grouting time, is the yield shear stress of Bingham cement slurry, is the tortuosity of the porous medium, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

[0032] In summary, this method can fully consider multiple factors such as the tortuosity of porous media, the deadweight and angle of Bingham cement slurry, and accurately determine the columnar penetration grouting range of Bingham cement slurry, so as to improve the quality and reliability of grouting in ecological restoration projects, meet the construction needs in complex ecological engineering environments, and achieve effective restoration and sustainable development of ecosystems.

[0033] In the second aspect, the present invention provides a system for determining the range of columnar infiltration grouting considering the tortuosity of porous media, the system comprising: a data acquisition device, a data output device, a processor and a storage device, the storage device comprising a computer-readable storage medium, the computer-readable storage medium storing a computer program, the computer program comprising program instructions, the program instructions when executed by the processor enable the processor to implement the method for determining the range of columnar infiltration grouting considering the tortuosity of porous media provided by the present invention. The system can stably run the method provided by the present invention, improve the efficiency of determining the range of columnar infiltration grouting, and improve the practicality of the method provided by the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0035] Figure 1 A schematic flow chart of a method for determining a columnar infiltration grouting range taking into account the tortuosity of a porous medium according to an embodiment of the present invention;

[0036] Figure 2 A schematic diagram of mechanical analysis of a microelement in a grouting pipe according to an embodiment of the present invention;

[0037] Figure 3 It is a schematic diagram of the framework of a system for determining a columnar permeation grouting range taking into account the tortuosity of a porous medium according to an embodiment of the present invention. DETAILED DESCRIPTION

[0038] The specific embodiments of the present invention will be described in detail below. It should be noted that the embodiments described herein are only for illustration and are not intended to limit the present invention. In the following description, a large number of specific details are set forth in order to provide a thorough understanding of the present invention. However, it is obvious to those of ordinary skill in the art that these specific details do not need to be adopted to implement the present invention. In other examples, in order to avoid confusing the present invention, known circuits, software or methods are not specifically described.

[0039] Throughout the specification, references to "one embodiment," "an embodiment," "an example," or "an example" mean that a particular feature, structure, or characteristic described in conjunction with the embodiment or example is included in at least one embodiment of the present invention. Therefore, the phrases "in one embodiment," "in an embodiment," "an example," or "an example" appearing in various places throughout the specification do not necessarily all refer to the same embodiment or example. In addition, particular features, structures, or characteristics may be combined in one or more embodiments or examples in any suitable combination and / or subcombination. In addition, it should be understood by those of ordinary skill in the art that the figures provided herein are for illustrative purposes and that the figures are not necessarily drawn to scale.

[0040] It should be noted in advance that, in an optional embodiment, except for independent explanations, the same symbols or letters appearing in all formulas have the same meanings and values.

[0041] In an alternative embodiment, see Figure 1 The present invention provides a method for determining the range of columnar infiltration grouting considering the tortuosity of porous media, the method comprising the following steps:

[0042] S1. Collect porous media characteristic data and groundwater pressure at the selected grouting point, and collect Bingham cement slurry characteristic data at the same time.

[0043] Specifically, in this embodiment, the porous medium characteristic data specifically refers to the formation porosity at the grouting point, which can be measured by density method, gas adsorption method and other methods. The groundwater pressure is measured by arranging a water pressure gauge sensor at the grouting point. The Bingham cement slurry characteristic data includes yield stress and plastic viscosity, and the Bingham cement slurry characteristic data can be measured using professional instruments such as rheometers.

[0044] S2. Setting grouting parameters, wherein the grouting parameters include grouting pressure, grouting time and grouting pipe laying angle.

[0045] Specifically, in this embodiment, the grouting pipe laying angle refers to the angle between the grouting pipe and the horizontal plane. More specifically, the grouting pressure needs to comprehensively consider factors such as the porous medium strength, groundwater pressure, and expected grouting diffusion range at the grouting point location to avoid excessive grouting pressure causing stratum damage or too low grouting pressure to achieve the expected grouting effect; the grouting time needs to be set in combination with the grouting volume, grouting speed, and project progress requirements, and the solidification time of the Bingham cement slurry in the porous medium must also be considered to ensure that the grouting operation is completed before the Bingham cement slurry is completely solidified; the grouting pipe laying angle needs to be reasonably adjusted based on the geological structure around the grouting point location and the expected flow direction of the Bingham cement slurry, so that the Bingham cement slurry can be evenly diffused in the porous medium.

[0046] Furthermore, it can be seen from the above description that the specific values ​​of the grouting parameters need to be set according to actual conditions, and therefore no limitation is made here.

[0047] S3. Considering the tortuosity of the porous medium, the seepage velocity of the Bingham cement slurry is determined based on the rheological equation of the Bingham fluid and the mechanical analysis equation of the microelement in the grouting pipe.

[0048] Wherein, step S3 specifically includes the following steps:

[0049] S31. Perform mechanical analysis on the microelement in the grouting pipe to obtain the mechanical analysis equation, and then combine the rheological equation to obtain the shear rate equation of Bingham cement slurry.

[0050] Specifically, in this embodiment, see Figure 2 , for the grouting pipe with radius r and length The corresponding mechanical analysis equation of the microelement satisfies the following relationship:

[0051] ,

[0052] Where P is the pressure at the front end of the microelement in the grouting pipe, is the change in pressure of the grouting pipe along the length of the microelement, is the shear stress on the inner wall of the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

[0053] By simplifying the mechanical analysis equation, the shear stress can be obtained to satisfy the following relationship:

[0054] ,

[0055] The rheological equation of Bingham fluid is:

[0056] ,

[0057] in, is the plastic viscosity of Bingham cement slurry, is the yield shear stress of Bingham cement slurry, is the shear rate, is the flow rate of Bingham cement slurry near a certain point in the Bingham cement slurry The rate of change of the radius r of the microelement in the grouting pipe.

[0058] Substituting the calculation formula of shear stress into the rheological equation of Bingham fluid, we get the calculation formula of shear rate, that is, the shear rate equation satisfies the following relationship:

[0059] ,

[0060] S32, integrating the shear rate equation to obtain a Bingham cement slurry flow rate, and calculating the Bingham cement slurry flow rate in the grouting pipe according to the Bingham cement slurry flow rate.

[0061] Specifically, in this embodiment, the shear rate equation obtained in step S31 is integrated to obtain the Bingham cement slurry flow rate under different microelement radii that satisfies the following relationship:

[0062] ,

[0063] in, is the core radius of Bingham cement slurry, is the radius of the grouting pipe.

[0064] After obtaining the Bingham cement slurry flow rate under different microelement radii, the Bingham cement slurry flow rate Q in the grouting pipe can be further calculated. The Bingham cement slurry flow rate Q satisfies the following relationship:

[0065] ,

[0066] Among them, W is an algebraic expression set for ease of writing, and .

[0067] S33. Calculate the average flow rate of the Bingham cement slurry in the grouting pipe according to the Bingham cement slurry flow rate, and then obtain the reference seepage velocity of the Bingham cement slurry according to the Dupuit-Forchheimer formula.

[0068] Specifically, in this embodiment, the average flow rate of Bingham cement slurry satisfies the following relationship:

[0069] ,

[0070] in, is the average flow rate of Bingham cement slurry.

[0071] According to the calculation relationship of the Bingham cement slurry flow rate obtained in step S32, it can be further obtained that the average flow rate of the Bingham cement slurry satisfies the following relationship:

[0072] ,

[0073] Furthermore, according to the Dupuit-Forchheimer formula:

[0074] ,

[0075] in, is the reference seepage velocity of Bingham cement slurry, It is the characteristic data of porous media.

[0076] Substituting the calculation formula of the average flow rate of Bingham cement slurry into the Dupuit-Forchheimer formula, we can get:

[0077] ,

[0078] S34. After grouting, the reference seepage velocity is corrected by taking into account the tortuosity of the porous medium to obtain the seepage velocity.

[0079] Specifically, in this embodiment, according to the calculation relationship of the Bingham cement slurry flow rate obtained in step S32, when the Bingham cement slurry flow rate is 0, there is:

[0080] ,

[0081] Solving this relationship, we can obtain:

[0082] ,

[0083] in, is the starting pressure gradient of Bingham cement slurry.

[0084] When grouting begins, the pressure gradient of Bingham cement slurry will be much greater than the starting pressure gradient of Bingham cement slurry, that is:

[0085] ,

[0086] Therefore, the calculation relationship of the reference seepage velocity of Bingham cement slurry can be simplified to obtain:

[0087] ,

[0088] Right now:

[0089] ,

[0090] In the actual seepage process, the infiltration trajectory of Bingham cement slurry is not a straight line. This is because Bingham cement slurry is affected by the internal structure of the soil, so its infiltration trajectory will be a curved path. In view of this, it is necessary to simplify the calculation relationship of the reference seepage velocity Make corrections.

[0091] Since the tortuosity of porous media satisfies the following relationship:

[0092] ,

[0093] ,

[0094] in, is the tortuosity of the porous medium, is the actual flow path length of Bingham cement slurry in porous media, It is the actual flow path length of Bingham cement slurry in the direction perpendicular to the horizontal plane.

[0095] therefore , using the actual change in flow path length Replace the simplified reference seepage velocity calculation formula with And sorting it out, we can get:

[0096] ,

[0097] Where V is the seepage velocity.

[0098] S4. Under the columnar diffusion boundary condition, a columnar penetration grouting range solution model is constructed based on the seepage velocity and the grouting volume per unit time calculation formula.

[0099] Wherein, step S4 specifically includes the following steps:

[0100] S41, establishing a pressure gradient equation according to the seepage velocity and the calculation formula of the grouting volume per unit time.

[0101] Specifically, in this embodiment, the grouting volume per unit time satisfies the following relationship:

[0102] ,

[0103] Among them, q is the grouting volume per unit time, and h is the columnar diffusion height.

[0104] therefore , combined with the calculation formula of the seepage velocity obtained in step S34, the pressure gradient equation satisfies the following relationship:

[0105] ,

[0106] S42, integrating the pressure gradient equation and taking into account the columnar diffusion boundary condition to obtain the columnar penetration grouting range solution model.

[0107] Specifically, in this embodiment, the columnar diffusion boundary condition is considered, that is, when hour, ,when hour, , is the grouting pressure, is the groundwater pressure, It is a quantitative index of the spatial diffusion range of columnar penetration grouting. That is, the diffusion distance of Bingham cement slurry in the columnar penetration grouting space.

[0108] More specifically, integrating the pressure gradient equation and considering the cylindrical diffusion boundary condition, we can obtain:

[0109] ,

[0110] Also because , so the columnar penetration grouting range solution model satisfies the following relationship:

[0111] ,

[0112] S5. According to the grouting parameters, a columnar penetration grouting range solution model is used to solve a quantitative index of the spatial diffusion range of the columnar penetration grouting to determine the columnar penetration grouting range.

[0113] Specifically, in this embodiment, the grouting parameters set in real time are brought into the columnar penetration grouting range solution model to solve the quantitative index of the columnar penetration grouting space diffusion range, and then determine the diffusion distance of Bingham cement slurry in the columnar penetration grouting space.

[0114] It should be noted that, in some cases, the actions described in the specification can be performed in a different order and still achieve the desired results. In this embodiment, the order of steps given is only to make the embodiment appear clearer and easier to explain, rather than to limit it.

[0115] In an alternative embodiment, see Figure 3The present invention also provides a columnar infiltration grouting range determination system considering the tortuosity of porous media, the columnar infiltration grouting range determination system considering the tortuosity of porous media comprises: a data acquisition device 1, a data output device 2, a processor 3 and a storage 4, the storage 4 comprises a computer-readable storage medium, the computer-readable storage medium stores a computer program, the computer program comprises program instructions, and when the program instructions are executed by the processor 3, the processor 3 implements the contents of steps S1 to S5 of the embodiment.

[0116] Specifically, in this embodiment, the data acquisition device 1 includes instruments such as a water pressure gauge sensor and a rheometer and a digital display screen, which are used to collect or input porous medium characteristic data, groundwater pressure, Bingham cement slurry characteristic data and grouting parameters. The data output device 2 also includes a digital display screen for displaying the solved quantitative index of the spatial diffusion range of the columnar infiltration grouting.

[0117] The system can stably operate the method, improve the efficiency of determining the range of columnar penetration grouting, and improve the practicability of the method.

[0118] In summary, this method can fully consider multiple factors such as the tortuosity of porous media, the deadweight and angle of Bingham cement slurry, and accurately determine the columnar penetration grouting range of Bingham cement slurry, so as to improve the quality and reliability of grouting in ecological restoration projects, meet the construction needs under complex ecological engineering environments, and achieve effective restoration and sustainable development of ecosystems. In addition, this system can stably run this method, improve the efficiency of determining the columnar penetration grouting range, and improve the practicality of this method.

[0119] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein by equivalents. These modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be included in the scope of the claims and specification of the present invention.

Claims

1. A method for determining the range of columnar infiltration grouting considering the tortuosity of porous media, characterized in that: The steps include: Collect porous media property data and groundwater pressure at the selected grouting points, and collect Bingham cement slurry property data; Setting grouting parameters, including grouting pressure, grouting time and grouting pipe laying angle; Taking into account the tortuosity of porous media, the seepage velocity of Bingham cement slurry is determined based on the rheological equation of Bingham fluid and the mechanical analysis equation of micro-element in the grouting pipe. Establishing a pressure gradient equation according to the seepage velocity and the grouting volume per unit time calculation formula; The pressure gradient equation satisfies the following relationship: ， in, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the tortuosity of the porous medium, is the plastic viscosity of Bingham cement slurry, q is the grouting volume per unit time, is the actual flow path length of Bingham cement slurry in the direction perpendicular to the horizontal plane, h is the columnar diffusion height, is the radius of the grouting pipe, is the porous medium characteristic data, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the yield shear stress of Bingham cement slurry, is the angle between the Bingham cement slurry flow direction and the grouting pipe axis; Integrating the pressure gradient equation and considering the columnar diffusion boundary condition, a columnar penetration grouting range solution model is obtained; The columnar penetration grouting range solution model satisfies the following relationship: ， in, is the grouting pressure, is the groundwater pressure, is the plastic viscosity of Bingham cement slurry, is the radius of the grouting pipe, is a quantitative index of the spatial diffusion range of columnar penetration grouting, t is the grouting time, is the yield shear stress of Bingham cement slurry, is the tortuosity of the porous medium, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the angle between the Bingham cement slurry flow direction and the grouting pipe axis; According to the grouting parameters, a columnar penetration grouting range solution model is used to solve the quantitative index of the spatial diffusion range of the columnar penetration grouting, so as to determine the columnar penetration grouting range.

2. The method for determining the range of columnar infiltration grouting considering the tortuosity of porous media according to claim 1, characterized in that: The method of considering the tortuosity of the porous medium and determining the seepage velocity of the Bingham cement slurry based on the rheological equation of the Bingham fluid and the mechanical analysis equation of the microelement in the grouting pipe comprises the following steps: Performing mechanical analysis on the microelement in the grouting pipe to obtain the mechanical analysis equation, and then combining the rheological equation to obtain the shear rate equation of Bingham cement slurry; Integrating the shear rate equation to obtain a Bingham cement slurry flow rate, and calculating the Bingham cement slurry flow rate in the grouting pipe according to the Bingham cement slurry flow rate; Calculating the average flow rate of the Bingham cement slurry in the grouting pipe according to the Bingham cement slurry flow rate, and then obtaining the reference seepage velocity of the Bingham cement slurry according to the Dupuit-Forchheimer formula; After grouting, the reference seepage velocity is corrected by taking into account the tortuosity of the porous medium to obtain the seepage velocity.

3. The method for determining the range of columnar infiltration grouting considering the tortuosity of porous media according to claim 2, characterized in that: The shear rate equation satisfies the following relationship: ， in, is the shear rate, is the flow rate of Bingham cement slurry near a certain point in the Bingham cement slurry As the radius r of the microelement in the grouting pipe changes, is the plastic viscosity of Bingham cement slurry, is the yield shear stress of Bingham cement slurry, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

4. The method for determining the range of columnar infiltration grouting considering the tortuosity of porous media according to claim 2, characterized in that: The Bingham cement slurry flow rate and the Bingham cement slurry flow rate respectively satisfy the following relationship: ， ， in, is the Bingham cement slurry flow rate, r is the radius of the microelement in the grouting pipe, is the plastic viscosity of Bingham cement slurry, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the angle between the Bingham cement slurry flow direction and the grouting pipe axis, is the yield shear stress of Bingham cement slurry, is the core radius of Bingham cement slurry, is the radius of the grouting pipe, and Q is the Bingham cement slurry flow rate.

5. The method for determining the range of columnar infiltration grouting considering the tortuosity of porous media according to claim 2, characterized in that: The average flow rate of the Bingham cement slurry satisfies the following relationship: ， in, is the average flow velocity of the Bingham cement slurry, Q is the flow rate of the Bingham cement slurry, is the radius of the grouting pipe.

6. The method for determining the range of columnar infiltration grouting considering the tortuosity of porous media according to claim 2, characterized in that: The seepage velocity satisfies the following relationship: ， Wherein, V is the seepage velocity, is the porous medium characteristic data, is the radius of the grouting pipe, is the plastic viscosity of Bingham cement slurry, is the tortuosity of the porous medium, is the change in pressure of the grouting pipe along the length of the microelement, is the length of the microelement in the grouting pipe, is the density of Bingham cement slurry, is the acceleration due to gravity, is the laying angle of the grouting pipe, is the yield shear stress of Bingham cement slurry, It is the angle between the flow direction of Bingham cement slurry and the axis of the grouting pipe.

7. A columnar infiltration grouting range determination system considering the tortuosity of porous media, characterized in that: The system for determining the range of columnar penetration grouting considering the tortuosity of porous media comprises: a data acquisition device, a data output device, a processor and a storage device, wherein the storage device comprises a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, wherein the computer program comprises program instructions, and when the program instructions are executed by the processor, the processor implements the method for determining the range of columnar penetration grouting considering the tortuosity of porous media as described in any one of claims 1 to 6.

Images