Slurry diffusion simulation method and device based on disturbance flow distribution

Abstract

The invention provides a slurry diffusion simulation method and device based on disturbance flow distribution, and belongs to the technical field of rock mass engineering. According to the method, for a target fractured rock mass with slurry distribution to be simulated, fissure statistical parameters of the target fractured rock mass are obtained firstly, then a random fissure network model of the target fractured rock mass is established according to the fissure statistical parameters, mesh generation is conducted on the random fissure network model, and a mesh generation result is obtained; then, according to a unified pipe network method and a mesh generation result, obtaining a fractured rock mass permeability result before grouting, and finally, according to a slurry probability diffusion algorithm based on disturbance flow distribution and the fractured rock mass permeability result before grouting, performing fracture network grouting condition simulation on the target fractured rock mass. And obtaining a slurry diffusion simulation result of the target fractured rock mass. According to the application, the interference factors are applied in the flow distribution process, so that the situation of slurry distribution difference possibly caused by the interference of uncertain factors can be simulated.

Description

Technical Field

[0001] The present application relates to the field of rock mass engineering technology, and in particular to a slurry diffusion simulation method and device based on disturbance flow distribution. Background Art

[0002] In recent years, the number of rock mass projects used for water storage, power generation, flood control, transportation, mining and other functions has gradually increased worldwide. However, geological disasters such as seepage, leakage and water inrush pose a potential threat to the safety of the above projects during construction and operation, especially under adverse geological conditions such as karst fissures and faults, the destructive capacity and probability of geological disasters are significantly increased.

[0003] At present, rock mass projects such as pumped storage, tunnel projects, large dams and underground buildings, by grouting into the rock mass, form a waterproof curtain or water-stop curtain around the building facilities and the rock mass as a protective barrier for the project body. Therefore, the effect of grouting directly determines the anti-seepage performance of the curtain. However, due to the hidden nature of the project and many uncertain factors, such as the complex engineering geological conditions (such as grouting splitting, confining pressure and special geological conditions) and changes in construction conditions (such as construction technology, grouting pressure, water-cement ratio, construction time and human economic factors) during the actual grouting process, the flow of slurry inside the rock mass is more likely to be disturbed by the aforementioned uncertain factors, making the quality of the curtain formation also extremely uncertain. At the same time, existing simulation methods are usually constructed based on deterministic factors. The specific distribution of the slurry is obtained by accurately modeling and determining the flow form of the slurry (such as assuming the flow form of the slurry, the flow equation, and the coupling effect of rock, soil and water). When the flow form of the slurry inside the rock and soil and the flow equation it follows change under different geological conditions and construction conditions, the simulation method based on deterministic factors lacks flexibility and the simulation effect is poor.

[0004] Therefore, there is an urgent need for a method that can simulate the distribution of slurry inside the rock mass when the flow state is unknown. Summary of the invention

[0005] The embodiments of the present application provide a method and device for simulating slurry diffusion based on disturbed flow distribution to solve the problem of lack of a method capable of simulating the distribution of slurry inside a rock mass when the flow state is unknown.

[0006] In a first aspect, an embodiment of the present application provides a slurry diffusion simulation method based on disturbance flow distribution, comprising:

[0007] Obtaining fracture statistical parameters of the target fractured rock mass; wherein the fracture statistical parameters include fracture length, fracture inclination and distribution function parameters obeyed by the fracture center point;

[0008] According to the fracture statistical parameters, a random fracture network model of the target fractured rock mass is established; wherein the fracture information of the random fracture network model includes fracture coordinates and fracture numbers;

[0009] Meshing the random fracture network model to obtain meshing results; wherein the meshing results include mesh unit numbers, mesh point coordinates, mesh point numbers, mesh points and edges corresponding to fractures and their numbers;

[0010] According to the unified pipe network method and grid division results, the permeability results of the fractured rock mass before grouting are obtained; among them, the permeability results of the fractured rock mass before grouting include the permeability coefficient, the pressure cloud map of the fractured rock mass, and the flow values ​​at various locations in the fractured rock mass;

[0011] According to the slurry probability diffusion algorithm based on disturbance flow distribution and the permeability results of fractured rock mass before grouting, the fracture network grouting situation of the target fractured rock mass is simulated to obtain the slurry diffusion simulation results of the target fractured rock mass.

[0012] In a second aspect, an embodiment of the present application provides a slurry diffusion simulation device based on disturbance flow distribution, comprising:

[0013] A fracture statistical parameter acquisition module is used to acquire fracture statistical parameters of a target fractured rock mass; wherein the fracture statistical parameters include fracture length, fracture inclination and distribution function parameters obeyed by the fracture center point;

[0014] A model building module is used to build a random fracture network model of the target fractured rock mass according to fracture statistical parameters; wherein the fracture information of the random fracture network model includes fracture coordinates and fracture numbers;

[0015] A mesh generation module is used to perform mesh generation on the random fracture network model to obtain mesh generation results; wherein the mesh generation results include mesh unit numbers, grid point coordinates, grid point numbers, grid points and edges corresponding to fractures and their numbers;

[0016] A generation module is used to obtain the permeability results of the fractured rock mass before grouting according to the unified pipe network method and the grid division results; wherein the permeability results of the fractured rock mass before grouting include the permeability coefficient, the pressure cloud map of the fractured rock mass, and the flow value at each location in the fractured rock mass;

[0017] The slurry diffusion simulation module is used to simulate the fracture network grouting situation of the target fractured rock mass according to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability results of the fractured rock mass before grouting, and obtain the slurry diffusion simulation results of the target fractured rock mass.

[0018] The embodiment of the present application provides a method and device for simulating slurry diffusion based on disturbed flow distribution. For the target fractured rock mass to be simulated for slurry distribution, first obtain the fracture statistical parameters of the target fractured rock mass, then establish a random fracture network model of the target fractured rock mass based on the fracture statistical parameters, then mesh the random fracture network model to obtain the meshing result, then obtain the permeability result of the fractured rock mass before grouting based on the unified pipe network method and the meshing result, finally simulate the fracture network grouting of the target fractured rock mass based on the slurry probability diffusion algorithm based on disturbed flow distribution and the permeability result of the fractured rock mass before grouting, and obtain the slurry diffusion simulation result of the target fractured rock mass. In this way, a method is provided that can simulate the distribution of slurry inside the rock mass when the flow state is unknown. Since the interference factor is applied in the flow distribution process, the slurry distribution difference that may be caused by the interference of the uncertainty factor can be simulated. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0020] Figure 1 It is a schematic diagram of the principle of a distortion model provided in an embodiment of the present application;

[0021] Figure 2 It is a schematic diagram of a fracture generation domain, a interception frame and a fracture rock mass permeability calculation unit provided in an embodiment of the present application;

[0022] Figure 3 It is a schematic diagram of interception and correction of a 50° direction crack model provided in an embodiment of the present application;

[0023] Figure 4 It is a schematic diagram of a 50° direction crack network partitioning result provided in an embodiment of the present application;

[0024] Figure 5 is a schematic diagram of pore pressure provided in an embodiment of the present application;

[0025] Figure 6 is a schematic diagram of a method for naming points in a crack provided in an embodiment of the present application;

[0026] Figure 7 This is a schematic diagram of an equal volume method provided in an embodiment of the present application;

[0027] Figure 8is a schematic diagram of a crack network before and after pretreatment provided in an embodiment of the present application;

[0028] Fig. 9 It is a schematic diagram of flow network division provided in an embodiment of the present application;

[0029] Fig.10 is a schematic diagram of a slurry diffusion simulation at different stages provided in an embodiment of the present application;

[0030] Fig.11 This is a method provided by the embodiment of the present application. 8-15 Schematic diagram of the volume ratio of each part;

[0031] Fig.12 It is a schematic diagram of the effect of disturbance flow distribution on slurry diffusion provided in an embodiment of the present application;

[0032] Fig.13 It is a schematic diagram of a slurry probability diffusion algorithm provided in an embodiment of the present application. DETAILED DESCRIPTION

[0033] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0034] In order to make the purpose, technical solutions and advantages of the present application clearer, specific embodiments will be described below in conjunction with the accompanying drawings.

[0035] As described in the related art, when the flow pattern of slurry in the rock mass and the flow equation it follows change under different geological conditions and construction conditions, the simulation method based on deterministic factors lacks flexibility and has poor simulation effect. Therefore, there is an urgent need for a method that can simulate the distribution of slurry inside the rock mass when the flow state is unknown.

[0036] In order to solve the problems of the prior art, the embodiment of the present application provides a slurry diffusion simulation method and device based on disturbed flow distribution. The slurry diffusion simulation method based on disturbed flow distribution provided by the embodiment of the present application is first introduced below.

[0037] The flow distribution introduced in this application refers to determining the amount of grouting liquid by calculating the proportion of the flow rate of each fracture based on the flow value and the pressure value at both ends of each fracture, while ensuring that the slurry complies with the flow law, and accurately simulating the final distribution of the slurry in the fracture. At the same time, an interference factor is applied during the flow distribution process to simulate the difference in slurry distribution that may be caused by the interference of uncertain factors.

[0038] See also Figure 1 , which shows a flow chart of the implementation of the slurry diffusion simulation method based on disturbance flow distribution provided in an embodiment of the present application, which is described in detail as follows:

[0039] Step 110: Obtain fracture statistical parameters of the target fractured rock mass.

[0040] Specifically, the fracture statistical parameters include the fracture length, fracture inclination and the distribution function parameters obeyed by the fracture center point.

[0041] Step 120: Establish a random fracture network model of the target fractured rock mass according to fracture statistical parameters.

[0042] Specifically, the fracture information of the random fracture network model includes fracture coordinates and fracture numbers.

[0043] Step 130: mesh the random fracture network model to obtain a meshing result.

[0044] Specifically, the mesh generation result includes mesh unit numbers, mesh point coordinates, mesh point numbers, mesh points and edges corresponding to cracks and their numbers.

[0045] Step 140: Obtain the permeability results of the fractured rock mass before grouting based on the unified pipe network method and the grid division results.

[0046] Specifically, the permeability results of the fractured rock mass before grouting include the permeability coefficient, the pressure cloud diagram of the fractured rock mass, and the flow values ​​at various locations in the fractured rock mass;

[0047] Step 150: According to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability result of the fractured rock mass before grouting, the fracture network grouting condition of the target fractured rock mass is simulated to obtain the slurry diffusion simulation result of the target fractured rock mass.

[0048] The following specifically introduces how to simulate the fracture network grouting situation of the target fractured rock mass based on the slurry probabilistic diffusion algorithm of disturbed flow distribution and the permeability results of the fractured rock mass before grouting.

[0049] 1) Preprocess the fracture network of the random fracture network model.

[0050] Specifically, the preprocessing includes breaking up all the fractures in the fracture network into small fractures whose endpoints are connected to each other based on intersections, bifurcation points and blind fracture endpoints.

[0051] 2) Obtain the pressure value of each node in the preprocessed fracture network and the flow value of each fracture.

[0052] 3) According to the pressure value of each node, the flow direction of the slurry in each fracture is obtained, and the blind-end fractures and isolated fractures without slurry flow in the fracture network are deleted.

[0053] 4) According to the flow direction of the slurry in each fracture, the deleted fracture network is divided into flow networks;

[0054] 5) According to the flow distribution principle, the pressure value of each node and the flow value of each fracture, the divided flow network is grouted. At the same time, a disturbance factor is applied during the grouting process to simulate the flow distribution of the slurry under the interference of uncertain factors.

[0055] Specifically, the traffic distribution principle can be:

[0056]

[0057]

[0058]

[0059] Among them, α ij , is the proportionality coefficient, Q ij is the flow value, n is the number of fractures of the same level, ij is the node number at both ends of the fracture, μ is the mean, σ is the standard deviation, V r V is the remaining slurry volume after the upper fracture is filled. ij To distribute the resulting slurry, α ij is the proportionality coefficient.

[0060] Method for controlling disturbance distribution: 1. First, according to formula (1.1), calculate the distribution ratio α of each crack in the same priority level ij ; 2. Establish the slurry disturbance distribution function. Taking the slurry disturbance distribution obeying the normal distribution function as an example, establish the distribution function shown in formula (1.2), where μ is the distribution ratio α ij , σ is the standard deviation; 3. After the distribution function is established, random sampling is used as the new distribution ratio coefficient 4. Substitute into formula (1.3) and recalculate the grout value V obtained in the crack ij In short, it is to control the value of the standard deviation σ in formula (1.2) so that the distribution ratio changes from α to ij The disturbance is Ultimately, the disturbance of slurry flow distribution is achieved.

[0061] In some embodiments, after obtaining the slurry diffusion simulation results of the target fractured rock mass, the permeability results of the fractured rock mass after grouting can also be obtained according to the unified pipe network method and the slurry diffusion simulation results of the target fractured rock mass. Furthermore, the permeability comparison results of the target fractured rock mass before and after grouting can be obtained according to the permeability results of the fractured rock mass after grouting and the permeability results of the fractured rock mass before grouting.

[0062] In order to facilitate understanding of the specific solution of the present application, the implementation process of each step is specifically introduced below.

[0063] 1.1 Collecting fracture statistics

[0064] For the target fractured rock mass for slurry flow simulation, the fracture statistical parameters can be collected by plotting tables: fracture length, fracture inclination and the distribution function parameters obeyed by the fracture center point, as shown in Table 1.

[0065] Table 1 Statistics of crack parameters

[0066]

[0067] 1.2 Establishing a random fracture network model

[0068] A random fracture network can be generated by using a random fracture generation program. Different fracture network models can be randomly generated by controlling the mean and standard deviation of the fracture group spacing distribution function, the mean and variance of the fracture length distribution function, and the mean and variance of the fracture inclination distribution function.

[0069] The specific operation process is as follows:

[0070] 1) Input the parameters in Table 1 into the random crack generation program and arbitrarily set the crack generation domain to 200×200;

[0071] 2) Establish a fractured rock mass permeability calculation model such as Figure 2 (right), output fracture network partition parameter table;

[0072] Among them, the permeability calculation unit of the fractured rock mass is selected, and the process is as follows: a 100×100 unit interception box is established in the center of the fracture generation domain ( Figure 2 The square area in the intercepted box is used as the calculation unit model of fractured rock mass permeability. The fracture permeability model in each direction is mainly to intercept the box and rotate it counterclockwise at a fixed angle, intercept the model, and normalize it, such as Figure 2 As shown in the right picture. Taking the study of the permeability of fractured rock mass at 50° as an example, the corresponding fracture network model is established, as shown in Figure 3 As shown in the picture on the right.

[0073] 1.3 Meshing

[0074] Input the fracture network model in Section 1.2 into the DFN meshing program and output the fracture meshing results as follows: Figure 4 shown.

[0075] 1.4 Calculation of permeability coefficient of fractured rock mass before grouting

[0076] Input the segmentation results, fractured rock mass parameters and slurry fluidity-related parameters in Section 1.3: fractured rock mass permeability (m2), porosity, fracture width (mm), fluid viscosity (Pa·s) and density (kg / m3), and boundary setting conditions (water head difference) into the unified pipe network method calculation program;

[0077] Output the calculation results of crack permeability before grouting - permeability coefficient K, flow value Q of each crack i and pore pressure value P i , and draw the pore pressure cloud map as shown in Figure 5 .

[0078] 1.5 Slurry diffusion simulation based on disturbed flow distribution

[0079] Based on the random fracture network model established in Section 1.2 and the calculation results in Section 1.4, the fracture network grouting situation is simulated.

[0080] The preprocessing is to break up all the cracks into small cracks with only the endpoints connected to each other, such as Figure 6 The crack formed by endpoints 1 and 14 is divided into 6 cracks, namely, cracks 1-6, cracks 6-12 and cracks 12-14. Similarly, the cracks formed by endpoints 3 and 10, endpoints 2 and 15, endpoints 5 and 9, and endpoints 11 and 13 are all broken up into small cracks.

[0081] When simulating grouting, it is assumed that the grout volume is incompressible and the grouting is performed using the equal volume method. At the same time, the movement of the grout follows the flow law. The equal volume method means that for each grouting fracture, the internal volume ratio is as follows: Figure 7 , that is, the ratio of the slurry volume to the volume of the fracture filled is 1:1.

[0082] Slurry flow law refers to:

[0083] 1) The slurry always flows from the area with high pore pressure to the area with low pore pressure. The pressure values ​​at each point are calculated in Section 1.4;

[0084] 2) Only when the fracture is filled will the slurry continue to flow to other fractures intersecting with the fracture.

[0085] by Figure 6For example, fractures 6-7 and 6-12 will be allocated grout only when fractures 1-5 are filled; 3) The grout is allocated step by step within the fracture network. Step by step allocation means that the fractures are divided into different priorities according to the order of contact with the grouting holes or intersecting fractures, and the grout is uniformly allocated according to different priority levels. Figure 7 If the slurry is distributed step by step in the cracks shown, the sequence is shown in Table 2.

[0086] Table 2 Allocation order table

[0087] Rift Level Crack Number First level 1-6、2-4 Level 2 6-5、6-7、6-12、4-7、4-8、4-3 Level 3 7-8、7-10、8-9、8-15、12-11、12-13、12-14 Level 4 8-9、8-15

[0088] The flow distribution principle means that when the upper fracture is filled, the proportion of the remaining slurry flowing to the lower connected fractures must be equal to the proportion of the flow between the fractures, that is, the following formula:

[0089]

[0090]

[0091]

[0092] Among them, α ij , is the proportionality coefficient, Q ij is the flow value, n is the number of fractures of the same level, ij is the node number at both ends of the fracture, μ is the mean, σ is the standard deviation, V r V is the remaining slurry volume after the upper fracture is filled. ij To distribute the resulting slurry, α ij is the proportionality coefficient.

[0093] by Figure 6 Taking the fracture network shown in the figure as an example, the simulated grouting process is as follows:

[0094] 1) For ease of understanding, the crack is named after the node numbers at both ends. ij , where ij means that the pressure at node i is greater than that at node j; the volume of each fracture is named F ij ; The volume of grout injected into each crack is V ij ; The flow rate of each crack is Q ij .like Figure 8 Middle fissure 1-6 is l 1-6 , fracture volume F 1-6 , the volume of slurry poured in is V 1-6 , the flow rate is Q 1-2 ;

[0095] 2) Fracture network preprocessing;

[0096] 3) Get the pressure value P of each node iand the flow rate of each fracture Q ij ;

[0097] 4) According to the pressure value P of each node i , the flow direction of the slurry in each fracture is obtained, and the blind-end fractures and isolated fractures without slurry flow are deleted (taking the blind-end fracture composed of endpoint 3 and node 4 as an example, at this time P3>P4, the slurry does not flow in the fracture, so it is deleted), so that Figure 8 (b);

[0098] 5) Depending on the flow direction, Figure 8 (b) The fracture network is divided into Fig. 9 The No. 1 and No. 2 flow nets are shown. Next, the No. 1 and No. 2 flow nets are grouted separately.

[0099] Take the No. 1 flow net as an example:

[0100] According to the flow distribution principle, that is, the fracture l 1-6 The total slurry volume V is distributed to flow net No. 1 and flow net No. 2 based on the flow ratio of fissures l2-4.

[0101] Assume that the volume of slurry distributed by No. 1 flow net is V 1-6 , crack l 1-6 F 1-6 (Assuming V 1-6 >F 1-6 ), then it flows into the crack l 6-7 , l 6-12 The total volume of the slurry is V 1-6 –F 6-12 = Vr, at this time the fracture perfusion situation is as follows Fig.10 (a); then with crack l 6-7 , l 6-12 Substitute the flow value into formula (1.1) and when no disturbance occurs, Substituting into formula (1.3), we can allocate V r At this time, the perfusion situation is as follows Fig.10 (b); then with crack l 7-8 , l 12-13 , l 12-14 The corresponding flow value is allocated, and the perfusion situation is as follows Fig.10 (c); Finally, the crack is filled 8-15 At this time, the perfusion situation is as follows Fig.10 (d);

[0102] When simulating the injection of the No. 2 flow network, in order to avoid cracks such as 8-15 There are two conflicting allocation orders, and each fracture is given a new parameter - the remaining fillable volume V of the fracture. z , when the fracture network is not grouted, V zThe initial value of is equal to the volume value V of the corresponding small crack ij , as the amount of grouting in the fracture gradually increases, its value gradually tends to 0, and when the fracture is subsequently grouted, it is regarded as a fully grouting fracture;

[0103] Then introduce the remaining fillable volume V z After that, the grouting process of No. 2 flow net becomes:

[0104] 1) Fill the gap 2-4 : Correction of crack l 2-4 The remaining pourable volume is 0;

[0105] 2) Fill the gaps 4-7 and l 4-8 : Correction of crack l 4-7 and l 4-8 The remaining pourable volume is 0;

[0106] 3) Filling the cracks 7-8 and l 8-15 :At this time in l 7-8 Before filling, the crack l 8-15 The slurry comes from the fissure 4-8 , assuming that the crack l 8-15 The slurry volume is F 8-15 , then the crack l 8-15 The remaining pourable volume V z V 8-15 -F 8-15 >0;

[0107] 4) Fill the gaps 7-8 :Fissure l 8-15 The slurry comes from the fissures 7-8 and l 4-8 , and finally 8-15 All slurry sources and their volume proportions are shown as follows Fig.11 shown.

[0108] It should be noted that the realization method of disturbance flow distribution is to adjust α in equation (1.1) during the slurry distribution process. ij Apply a disturbance to become a new Achieve disturbed distribution of slurry.

[0109] The specific method is as follows: Assume that the slurry flow obeys a certain distribution function, set the standard deviation σ to 0.5, and use α ij The value is taken as the mean μ of the distribution function, the distribution function of formula (1.2) is established, and random sampling is performed. The sampling results are Substitute into formula (1.3) and calculate the amount of slurry after redistribution.

[0110] by Fig.12For example, assuming that the slurry flow obeys the normal distribution function, the process of applying disturbance is as follows:

[0111] 1) In Fig.12 In the above example, we assume that the crack l 23 , l 24 The flow ratio is 1:3, and the fracture volumes are 3 and 2 respectively; 23 , l 24 The total amount of slurry is V r =4;

[0112] 2) Without applying the disturbance factor, the crack l 23 , l 24 The amount of slurry obtained is calculated according to the flow distribution by formula 1 and 3;

[0113] 3) With crack l 23 , l 24 The flow rate accounts for 0.25, 0.75 is the mean μ, let σ=0.5 and substitute it into formula (1.2) to establish the sampling function;

[0114] 4) Perform random sampling, take the sampling value as the slurry volume after disturbance, and re-simulate the slurry distribution at this time as follows Fig.12 (b) as shown.

[0115] like Fig.13 As shown, the slurry diffusion simulation process based on disturbed flow distribution is shown.

[0116] 1.6 Calculation of permeability coefficient

[0117] This model uses the unified pipe network method to calculate the permeability coefficient. The pipe network system equates the fracture to a one-dimensional fracture pipe and the rock matrix to a one-dimensional matrix pipe. In the process of establishing the grouting model, only two fluids (water and slurry) are considered to flow in the rock mass, and the hydraulic parameters of the two types of pipes for water or slurry are different. Considering the permeability of the rock mass to water and slurry in actual engineering, it is assumed that water uses fracture pipes and rock matrix pipes as flow channels, while slurry only uses fractures as the only flow channel. Therefore, after grouting, only the permeability of the fracture pipe to water is changed inside the rock mass, while the matrix pipe is not affected.

[0118] After grouting, the permeability coefficient of the fracture tube where the grout is located is modified as follows: the amount of grout allocated to each fracture is assumed to be V j The amount of grouting fluid in the rock mass is V 总 =ΣV jAfter grouting, the cracks can be divided into fully grouted cracks, semi-grouted cracks, and ungrouted cracks. Assume that the permeability coefficient of water through the cracks is K1, and the permeability coefficient of water through the grout stone is K2. Then the permeability coefficient of water through the fully grouted cracks changes from K1 to K2; the ungrouted cracks still maintain the original permeability coefficient K1; while the permeability coefficient of water through the semi-grouted cracks, K3, is determined according to the volume ratio of the grout in the cracks. The calculation of K3 is as follows: Assume that the volume of a certain semi-grouted crack is V, and the volume of the grout in this crack is V 浆 , where (V 浆 <V), then the permeability coefficient of water through this crack is calculated according to the following formula (3):

[0119]

[0120] The calculation process of the permeability coefficient of the unit body is as follows:

[0121] 1) The calculation unit selects a two-dimensional square random crack with a generation domain of 200×200 and a calculation domain of 100×100, as Figure 8 .

[0122] 2) The left boundary of the unit body is the inflow boundary of water and grout; the right boundary is the corresponding outflow boundary; the upper and lower boundaries are impermeable boundaries during the calculation of the permeability coefficient.

[0123] 3) The cracks in the unit body are smooth and straight, and do not expand or generate new cracks under pressure.

[0124] 4) The calculation of the permeability coefficient follows Darcy's law, which is formula (4):

[0125]

[0126] Among them, K is the permeability coefficient, Q is the flow rate, L is the length, ΔP is the pressure difference between the left and right boundaries of the unit body, and A is the area passed by the unit fluid. Note: For a two-dimensional unit body, the unit thickness is taken as 1, so A = L·1.

[0127] In the embodiment of the present application, for the target fractured rock mass to be simulated for slurry distribution, firstly, the fracture statistical parameters of the target fractured rock mass are obtained, then, according to the fracture statistical parameters, a random fracture network model of the target fractured rock mass is established, then the random fracture network model is meshed to obtain the meshing result, then, according to the unified pipe network method and the meshing result, the permeability result of the fractured rock mass before grouting is obtained, and finally, according to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability result of the fractured rock mass before grouting, the fracture network grouting situation of the target fractured rock mass is simulated to obtain the slurry diffusion simulation result of the target fractured rock mass. In this way, a method is provided that can simulate the distribution of slurry inside the rock mass when the flow state is unknown, and because the interference factor is applied in the flow distribution process, the slurry distribution difference that may be caused by the interference of the uncertainty factor can be simulated.

[0128] It should be understood that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0129] The following is an embodiment of the device of the present application. For details not described in detail, please refer to the corresponding method embodiment described above.

[0130] The slurry diffusion simulation device based on disturbance flow distribution in the embodiment of the present application includes:

[0131] A fracture statistical parameter acquisition module is used to acquire fracture statistical parameters of a target fractured rock mass; wherein the fracture statistical parameters include fracture length, fracture inclination and distribution function parameters obeyed by the fracture center point;

[0132] A model building module is used to build a random fracture network model of the target fractured rock mass according to fracture statistical parameters; wherein the fracture information of the random fracture network model includes fracture coordinates and fracture numbers;

[0133] A mesh generation module is used to perform mesh generation on the random fracture network model to obtain mesh generation results; wherein the mesh generation results include mesh unit numbers, grid point coordinates, grid point numbers, grid points and edges corresponding to fractures and their numbers;

[0134] A generation module is used to obtain the permeability results of the fractured rock mass before grouting according to the unified pipe network method and the grid division results; wherein the permeability results of the fractured rock mass before grouting include the permeability coefficient, the pressure cloud map of the fractured rock mass, and the flow value at each location in the fractured rock mass;

[0135] The slurry diffusion simulation module is used to simulate the fracture network grouting situation of the target fractured rock mass according to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability results of the fractured rock mass before grouting, and obtain the slurry diffusion simulation results of the target fractured rock mass.

[0136] In some embodiments, the slurry diffusion simulation module is further configured to:

[0137] Preprocessing the fracture network of the random fracture network model; wherein the preprocessing includes breaking up all fractures of the fracture network into small fractures whose only endpoints are interconnected based on intersections, bifurcations and blind fracture endpoints;

[0138] Obtain the pressure value of each node in the preprocessed fracture network and the flow value of each fracture;

[0139] According to the pressure value of each node, the flow direction of the slurry in each fracture is obtained, and the blind-end fractures and isolated fractures without slurry flow in the fracture network are deleted;

[0140] According to the flow direction of slurry in each fracture, the deleted fracture network is divided into flow networks;

[0141] According to the flow distribution principle, the pressure value of each node and the flow value of each fracture, the divided flow network is grouting, and the disturbance factor is applied during the grouting process to simulate the flow distribution of the slurry under the interference of uncertain factors.

[0142] In some embodiments, traffic distribution principles include:

[0143] When the upper fractures are filled, the ratio of the remaining slurry flowing to the lower connected fractures is equal to the ratio of the flow rates between the lower connected fractures.

[0144] In some embodiments, the slurry diffusion simulation module is specifically used to:

[0145] For the volume of grouting liquid to be filled in each fracture in the flow network, random sampling is performed using the distribution function obeyed by the grout flow, and the obtained sampling value is used as the disturbed grout volume of the corresponding fracture;

[0146] According to the volume of slurry after disturbance in each fracture in the flow network, the flow distribution of slurry under the interference of uncertain factors is simulated.

[0147] In some embodiments, traffic distribution principles include:

[0148] (When no disturbance occurs )

[0149] Among them, α ij , is the proportionality coefficient, Qij is the flow value, n is the number of fractures of the same level, ij is the node number at both ends of the fracture, V r V is the remaining slurry volume after the upper fracture is filled. ij The resulting slurry was partitioned.

[0150] The mean μ of the distribution function obeyed by the slurry flow is the proportional coefficient α ij , the standard deviation σ is a preset value.

[0151] In some embodiments, the distribution function to which the slurry flow is subject includes:

[0152]

[0153] Among them, σ is the standard deviation, and the mean μ is the proportional coefficient α ij .

[0154] In some embodiments, the slurry diffusion simulation device based on disturbance flow distribution further includes a comparison module for:

[0155] According to the unified pipe network method and the slurry diffusion simulation results of the target fractured rock mass, the permeability results of the fractured rock mass after grouting are obtained;

[0156] According to the permeability results of the fractured rock mass after grouting and the permeability results of the fractured rock mass before grouting, the permeability comparison results of the target fractured rock mass before and after grouting are obtained.

[0157] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0158] In addition, the features of the embodiments shown in the drawings of the present application or the various embodiments mentioned in this specification are not necessarily understood as independent embodiments. Instead, each feature described in one example of an embodiment can be combined with one or more other desired features from other embodiments to produce other embodiments not described in words or with reference to the drawings.

[0159] It should also be noted that the exemplary embodiments mentioned in this application describe some methods or systems based on a series of steps or devices. However, this application is not limited to the order of the above steps, that is, the steps can be performed in the order mentioned in the embodiment, or in a different order from the embodiment, or several steps can be performed simultaneously.

[0160] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application 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 make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A slurry diffusion simulation method based on disturbance flow distribution, characterized in that: include: Obtaining fracture statistical parameters of the target fractured rock mass; wherein the fracture statistical parameters include fracture length, fracture inclination and distribution function parameters obeyed by the fracture center point; According to the fracture statistical parameters, a random fracture network model of the target fractured rock mass is established; wherein the fracture information of the random fracture network model includes fracture coordinates and fracture numbers; Meshing the random fracture network model to obtain a meshing result; wherein the meshing result includes a mesh unit number, a mesh point coordinate, a mesh point number, a mesh point and an edge corresponding to a fracture and their numbers; According to the unified pipe network method and the grid division result, the permeability result of the fractured rock mass before grouting is obtained; wherein the permeability result of the fractured rock mass before grouting includes the permeability coefficient, the pressure cloud map of the fractured rock mass, and the flow value at each location in the fractured rock mass; According to the slurry probability diffusion algorithm based on disturbance flow distribution and the permeability result of the fractured rock mass before grouting, the fracture network grouting condition of the target fractured rock mass is simulated to obtain the slurry diffusion simulation result of the target fractured rock mass.

2. The slurry diffusion simulation method based on disturbance flow distribution according to claim 1 is characterized in that: The method of simulating the fracture network grouting of the target fractured rock mass according to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability result of the fractured rock mass before grouting includes: Preprocessing the fracture network of the random fracture network model; wherein the preprocessing includes breaking up all fractures of the fracture network into small fractures whose only endpoints are interconnected based on intersections, bifurcation points and blind fracture endpoints; Obtain the pressure value of each node in the preprocessed fracture network and the flow value of each fracture; According to the pressure value of each node, the flow direction of the slurry in each fracture is obtained, and the blind-end fractures and isolated fractures without slurry flow in the fracture network are deleted; According to the flow direction of slurry in each fracture, the deleted fracture network is divided into flow networks; According to the flow distribution principle, the pressure value of each node and the flow value of each fracture, the divided flow network is grouting, and the disturbance factor is applied during the grouting process to simulate the flow distribution of the slurry under the interference of uncertain factors.

3. The slurry diffusion simulation method based on disturbance flow distribution according to claim 2 is characterized in that: The traffic distribution principles include: When the upper fractures are filled, the ratio of the remaining slurry flowing to the lower connected fractures is equal to the ratio of the flow rates between the lower connected fractures.

4. The slurry diffusion simulation method based on disturbance flow distribution according to claim 3 is characterized in that: The disturbance factor is applied during the grouting process to simulate the flow distribution of the slurry under the interference of uncertain factors, including: For the volume of grouting liquid to be filled in each fracture in the flow network, random sampling is performed using the distribution function obeyed by the grout flow, and the obtained sampling value is used as the disturbed grout volume of the corresponding fracture; According to the volume of slurry after disturbance in each fracture in the flow network, the flow distribution of slurry under the interference of uncertain factors is simulated.

5. The slurry diffusion simulation method based on disturbance flow distribution according to claim 4 is characterized in that: The disturbance flow distribution principles include: Among them, V r is the volume of slurry remaining after the upper fracture is filled, n is the number of fractures at the same level, V ij is the volume of the slurry obtained, ij is the node number at both ends of the crack, α ij is the proportionality coefficient, Q ij is the flow value in the fracture ij, f(x) is the distribution function obeyed by the slurry flow, and μ is α ij , σ is the standard deviation.

6. The slurry diffusion simulation method based on disturbance flow distribution according to claim 1 is characterized in that: The method further comprises: According to the unified pipe network method and the slurry diffusion simulation results of the target fractured rock mass, the permeability results of the fractured rock mass after grouting are obtained; According to the permeability result of the fractured rock mass after grouting and the permeability result of the fractured rock mass before grouting, a permeability comparison result of the target fractured rock mass before and after grouting is obtained.

7. A slurry diffusion simulation device based on disturbance flow distribution, characterized in that: include: A fracture statistical parameter acquisition module, used to acquire fracture statistical parameters of a target fractured rock mass; wherein the fracture statistical parameters include fracture length, fracture inclination and distribution function parameters obeyed by the fracture center point; A model building module, used to build a random fracture network model of the target fractured rock mass according to the fracture statistical parameters; wherein the fracture information of the random fracture network model includes fracture coordinates and fracture numbers; A mesh generation module, used for meshing the random fracture network model to obtain a mesh generation result; wherein the mesh generation result includes a mesh unit number, a mesh point coordinate, a mesh point number, a mesh point and an edge corresponding to a fracture and their numbers; A generation module is used to obtain the permeability results of the fractured rock mass before grouting according to the unified pipe network method and the grid division results; wherein the permeability results of the fractured rock mass before grouting include permeability coefficient, fractured rock mass pressure cloud map, and flow values ​​at various locations in the fractured rock mass; The slurry diffusion simulation module is used to simulate the fracture network grouting situation of the target fractured rock mass according to the slurry probability diffusion algorithm based on the disturbance flow distribution and the permeability result of the fractured rock mass before grouting, so as to obtain the slurry diffusion simulation result of the target fractured rock mass.

8. The slurry diffusion simulation device based on disturbance flow distribution according to claim 7 is characterized in that: The slurry diffusion simulation module is also used to: Preprocessing the fracture network of the random fracture network model; wherein the preprocessing includes breaking up all fractures of the fracture network into small fractures whose only endpoints are interconnected based on intersections, bifurcation points and blind fracture endpoints; Obtain the pressure value of each node in the preprocessed fracture network and the flow value of each fracture; According to the pressure value of each node, the flow direction of the slurry in each fracture is obtained, and the blind-end fractures and isolated fractures without slurry flow in the fracture network are deleted; According to the flow direction of slurry in each fracture, the deleted fracture network is divided into flow networks; According to the flow distribution principle, the pressure value of each node and the flow value of each fracture, the divided flow network is grouting, and the disturbance factor is applied during the grouting process to simulate the flow distribution of the slurry under the interference of uncertain factors.

9. The slurry diffusion simulation device based on disturbance flow distribution according to claim 8, characterized in that: The traffic distribution principles include: When the upper fractures are filled, the ratio of the remaining slurry flowing to the lower connected fractures is equal to the ratio of the flow rates between the lower connected fractures.

10. The slurry diffusion simulation device based on disturbance flow distribution according to claim 7, characterized in that: The device also includes a comparison module, which is used to: According to the unified pipe network method and the slurry diffusion simulation results of the target fractured rock mass, the permeability results of the fractured rock mass after grouting are obtained; According to the permeability result of the fractured rock mass after grouting and the permeability result of the fractured rock mass before grouting, a permeability comparison result of the target fractured rock mass before and after grouting is obtained.

Images