Numerical simulation method and device for tunnel mud burst and water burst disasters

Abstract

The invention relates to a numerical simulation method and device for tunnel mud burst and water burst disasters. The method comprises the following steps: acquiring a three-dimensional discrete element model of a target area according to geological information of the target area, simulating a mud burst water gushing disaster process caused by tunnel excavation by utilizing the three-dimensional discrete element model based on tunnel information, obtaining a numerical simulation result of the mud burst water gushing disaster, adjusting an angle corresponding to a tunnel crossing fault area in the tunnel information, and obtaining a numerical simulation result of the mud burst water gushing disaster. And obtaining numerical simulation results of the mud burst and water burst disasters under different tunnel information. By adopting the method, the numerical simulation result of the mud burst water gushing disaster simulated by the three-dimensional discrete element model can be obtained at the angles corresponding to different tunnel crossing fault areas, the occurrence process of the tunnel mud burst water gushing can be sufficiently described subsequently, and the accuracy of the simulation result is improved.

Description

Technical Field

[0001] This application relates to the field of tunnel construction technology, and in particular to a numerical simulation method and apparatus for tunnel mudslides and water inrushes. Background Technology

[0002] Mud inrush and water inrush in tunnels are significant hazards affecting safe tunnel construction. Currently, traditional numerical simulation methods for tunnel mud inrush and water inrush disasters rely heavily on theoretical analysis and physical experiments, resulting in low accuracy of simulation results. Summary of the Invention

[0003] Therefore, it is necessary to provide a numerical simulation method and apparatus for tunnel mudslides and water inrushes that can improve the accuracy of simulation results, addressing the aforementioned technical problems.

[0004] Firstly, this application provides a numerical simulation method for tunnel mudslides and water inrushes, the method comprising:

[0005] Based on the geological information of the target area, a three-dimensional discrete element model of the target area is obtained; the three-dimensional discrete element model includes fault areas.

[0006] Based on tunnel information, a three-dimensional discrete element model was used to simulate the mudslide and water inrush disaster process caused by tunnel excavation, and numerical simulation results of the mudslide and water inrush disaster were obtained.

[0007] Adjust the angle corresponding to the fault region crossed by the tunnel in the tunnel information, and based on the adjusted tunnel information, repeat the above steps of simulating the mudslide and water inrush disaster caused by tunnel excavation using a three-dimensional discrete element model based on the tunnel information to obtain the numerical simulation results of the mudslide and water inrush disaster.

[0008] In one embodiment, the three-dimensional discrete element model also includes a surrounding rock region; the geological information includes physical parameters of the fault region, physical parameters of the rock mass in the surrounding rock region, and ductile plating information and brittle joint information of the surrounding rock region.

[0009] In one embodiment, a three-dimensional discrete element model of the target area is obtained based on the geological information of the target area, including:

[0010] Based on geological information, an initial three-dimensional discrete element model of the target area is constructed.

[0011] The initial three-dimensional discrete element model is subjected to geostress equilibrium calculation to obtain the three-dimensional discrete element model.

[0012] In one embodiment, an initial three-dimensional discrete element model of the target area is constructed based on geological information, including:

[0013] Based on the ductile joint information, the instructions for generating ductile structural surfaces for the initial three-dimensional discrete element model are determined, and based on the brittle joint information, the instructions for generating brittle structural surfaces for the initial three-dimensional discrete element model are determined.

[0014] Based on the instructions for generating ductile and brittle structural surfaces, corresponding ductile and brittle structural surface information is generated. Based on the ductile and brittle structural surface information, an initial three-dimensional discrete element model is constructed.

[0015] In one embodiment, geostress equilibrium calculations are performed on the initial three-dimensional discrete element model to obtain a three-dimensional discrete element model, including:

[0016] In the overlying rock mass region of the initial three-dimensional discrete element model, and based on fluid conditions, boundary displacement conditions, and stress conditions, geostress balance calculations are performed to obtain the three-dimensional discrete element model.

[0017] In one embodiment, based on tunnel information, a three-dimensional discrete element model is used to simulate the mudslide and water inrush disaster process caused by tunnel excavation, including:

[0018] The tunnel excavation was simulated using a three-dimensional discrete element model until the tunnel and fault region in the three-dimensional discrete element model were connected.

[0019] Secondly, this application also provides a numerical simulation device for tunnel mudslides and water inrushes, the device comprising:

[0020] The model acquisition module is used to acquire a three-dimensional discrete element model of the target area based on the geological information of the target area; the three-dimensional discrete element model includes fault areas.

[0021] The simulation module is used to simulate the mudslide and water inrush disaster process caused by tunnel excavation using a three-dimensional discrete element model based on tunnel information, and to obtain the numerical simulation results of the mudslide and water inrush disaster.

[0022] The simulation module is also used to adjust the angle corresponding to the fault area crossed by the tunnel in the tunnel information, and based on the adjusted tunnel information, repeat the above steps of simulating the mudslide and water inrush disaster caused by tunnel excavation using a three-dimensional discrete element model based on the tunnel information, and obtain the numerical simulation results of the mudslide and water inrush disaster.

[0023] Thirdly, this application also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described method.

[0024] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-described method.

[0025] Fifthly, this application also provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the above-described method.

[0026] The aforementioned numerical simulation method and apparatus for tunnel mudslide and water inrush disasters obtain a three-dimensional discrete element model of the target area based on the geological information of the target area, and simulates the mudslide and water inrush disaster process caused by tunnel excavation using the three-dimensional discrete element model based on tunnel information, thereby obtaining numerical simulation results of the mudslide and water inrush disaster. By adjusting the angle corresponding to the fault region crossed by the tunnel in the tunnel information, numerical simulation results of the mudslide and water inrush disaster under different tunnel information are obtained. This application obtains numerical simulation results of the mudslide and water inrush disaster simulated by the three-dimensional discrete element model under different angles corresponding to the fault region crossed by the tunnel, which facilitates the subsequent full characterization of the occurrence process of tunnel mudslide and water inrush and improves the accuracy of the simulation results. Attached Figure Description

[0027] To more clearly illustrate the technical solutions in the embodiments of this application or related technologies, the drawings used in the description of the embodiments of this application or related technologies will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0028] Figure 1 This is a flowchart illustrating a numerical simulation method for a tunnel mudslide and water inrush disaster in one embodiment;

[0029] Figure 2 This is a flowchart illustrating a numerical simulation method for tunnel mudslide and water inrush disasters in another embodiment;

[0030] Figure 3 This is a schematic diagram of the geological model corresponding to a 45-degree fault in one embodiment;

[0031] Figure 4 This is a tunnel displacement cloud map of a geological model in one embodiment;

[0032] Figure 5 A cloud map showing the hydraulic aperture of a geological model in one embodiment;

[0033] Figure 6 This is a pore pressure cloud map of a geological model in one embodiment, showing a joint that has been damaged and the flow plane in which it is located.

[0034] Figure 7 This is a structural block diagram of a numerical simulation device for tunnel mudslide and water inrush disaster in one embodiment;

[0035] Figure 8This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0036] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0037] It should be noted that the terms "first," "second," etc., used in this application can be used to describe various elements, but these elements are not limited by these terms. These terms are only used to distinguish the first element from the second element. The terms "comprising" and "having," and any variations thereof, used in this application, are intended to cover non-exclusive inclusion. The term "multiple" used in this application refers to two or more. The term "and / or" used in this application refers to one of the embodiments, or any combination of multiple embodiments.

[0038] Mudslides and water inrushes in tunnels are significant hazards affecting safe tunnel construction. Exploring the mechanisms of these hazards is crucial for improving tunnel construction safety and evaluating the effectiveness of reinforcement measures. Currently, traditional methods in studying the mechanisms of tunnel mudslides and water inrushes largely rely on theoretical analysis and physical experiments, lacking the ability to accurately simulate the dynamic evolution of the hydraulic coupling between fault fracture zones and the tunnel excavation face. Existing numerical models (such as continuous medium models) struggle to accurately characterize discontinuous behaviors like fracture propagation and rock fracturing, leading to insufficient accuracy in predicting water inrush paths—that is, the simulation results are not highly accurate. Furthermore, the seepage-stress coupling mechanism from fault fracture zones to the tunnel excavation face is unclear, existing models cannot effectively reveal the critical conditions for water inrush channel formation, and the mechanisms of mudslides and water inrushes under the geological conditions of brittle and highly sheared surrounding rocks in high-altitude areas have not been explored.

[0039] In one exemplary embodiment, such as Figure 1 As shown, a numerical simulation method for tunnel mudslides and water inrushes is provided. This embodiment illustrates the application of this method to a terminal; it is understood that this method can also be applied to a server, and to a system including both a terminal and a server, and implemented through interaction between the terminal and the server. The terminal can be, but is not limited to, various personal computers, laptops, smartphones, and tablets; the server can be a standalone server or a server cluster composed of multiple servers. In this embodiment, the method includes:

[0040] S102, Based on the geological information of the target area, obtain a three-dimensional discrete element model of the target area; wherein, the three-dimensional discrete element model includes fault areas.

[0041] The target area is the area where a tunnel needs to be excavated. Its specific location can be set according to the actual situation and is not limited in this embodiment.

[0042] Specifically, by obtaining a three-dimensional discrete element model corresponding to the target area based on the geological information of the target area, it can be understood that the target area includes fault areas, and correspondingly, the three-dimensional discrete element model also includes fault areas.

[0043] It should be noted that geological information can be obtained using existing technologies, which will not be elaborated upon in this embodiment. The method of obtaining a three-dimensional discrete element model based on geological information can be set according to the actual situation. In this embodiment, the three-dimensional discrete element model of the target area is obtained by using the discrete element method based on the geological information of the target area. In traditional methods, the application of the discrete element method in fault fracture zone modeling is still limited to rock fracture analysis and has not been extended to the dynamic simulation of hydraulic path evolution. However, this application utilizes the advantages of the discrete element method to break through the continuous medium assumption of traditional numerical models, which can provide a high-precision numerical prediction method for the prevention and control of tunnel mudslides and water inrush disasters under complex geological conditions, and provide theoretical basis and technical support for the support design and safety assessment of similar projects.

[0044] S104. Based on tunnel information, a three-dimensional discrete element model is used to simulate the mudslide and water inrush disaster process caused by tunnel excavation, and numerical simulation results of the mudslide and water inrush disaster are obtained.

[0045] Tunnel information refers to information related to tunnel excavation, such as the depth of tunnel cyclic excavation, the number and length of anchor bolts.

[0046] Specifically, based on tunnel information, a three-dimensional discrete element model is used to simulate tunnel excavation. After the calculation converges, the crack development status is output. The crack is the crack from the tunnel free face to the water-rich fault. If the crack from the tunnel free face to the water-rich fault has been connected, it indicates that a mudslide and water inrush disaster has occurred. Numerical simulation results of the mudslide and water inrush disaster are obtained.

[0047] It should be noted that the numerical simulation results of mudslide and water inrush disasters can include the distance between the tunnel free face and the fault. It can be understood that the distance between the tunnel free face and the fault obtained at the moment before the mudslide and water inrush disaster occurs is the safety critical value. At the same time, the numerical simulation results of mudslide and water inrush disasters can also be used to reflect the seepage path of mudslide and water inrush in the tunnel.

[0048] S106, Adjust the angle corresponding to the fault area crossed by the tunnel in the tunnel information, and based on the adjusted tunnel information, repeat the above steps of simulating the mudslide and water inrush disaster caused by tunnel excavation using a three-dimensional discrete element model based on the tunnel information, and obtain the numerical simulation results of the mudslide and water inrush disaster.

[0049] Specifically, in order to obtain the influence of the angles corresponding to different tunnel crossing fault regions on the occurrence of mudslides and water inrush disasters, the angles corresponding to the tunnel crossing fault regions in the tunnel information are adjusted, and the above steps for obtaining numerical simulation results of mudslides and water inrush disasters are repeated. This facilitates a full characterization of the occurrence process of mudslides and water inrushes in tunnels and improves the accuracy of the simulation results.

[0050] In the above-mentioned numerical simulation method for tunnel mudslide and water inrush disaster, a three-dimensional discrete element model of the target area is obtained based on the geological information of the target area. Based on the tunnel information, the tunnel excavation is simulated using the three-dimensional discrete element model. The numerical simulation results of mudslide and water inrush disaster simulated by the three-dimensional discrete element model are obtained at different angles corresponding to the tunnel crossing fault areas. This facilitates the subsequent full characterization of the occurrence process of tunnel mudslide and water inrush and improves the accuracy of the simulation results.

[0051] In one embodiment, the three-dimensional discrete element model also includes a surrounding rock region; the geological information includes physical parameters of the fault region, physical parameters of the rock mass in the surrounding rock region, and ductile plating information and brittle joint information of the surrounding rock region.

[0052] Specifically, the target area includes the surrounding rock area, and correspondingly, the three-dimensional discrete element model also includes the surrounding rock area. The surrounding rock in the target area can be the ductile-brittle strong shear surrounding rock of the plateau. The physical parameters of the fault area can include the fault angle, and the physical parameters of the rock mass in the surrounding rock area can include the density, Poisson's ratio, Young's modulus, cohesion, internal friction angle and tensile strength of the surrounding rock.

[0053] It should be noted that faults may include rock masses, and geological information also includes physical parameters of the contact surfaces between the fault and the rock masses within the fault, such as normal stiffness, shear stiffness, and cohesion.

[0054] For example, the three-dimensional discrete element model can be grouped into: tunnel region, surrounding rock region, fault region, and rock region in the fault, so that the three-dimensional discrete element model can be manipulated in the future and different regions of the three-dimensional discrete element model can be assigned corresponding physical properties.

[0055] In one embodiment, a three-dimensional discrete element model of the target area is obtained based on the geological information of the target area, including:

[0056] Based on geological information, an initial three-dimensional discrete element model of the target area is constructed.

[0057] The initial three-dimensional discrete element model is subjected to geostress equilibrium calculation to obtain the three-dimensional discrete element model.

[0058] The initial three-dimensional discrete element model can refer to a three-dimensional discrete element numerical model constructed based on the discrete element method.

[0059] Specifically, based on geological information, the discrete element method is used to construct an initial three-dimensional discrete element model of the target area, and the geostress balance calculation is performed on the initial three-dimensional discrete element model to obtain the three-dimensional discrete element model.

[0060] In one embodiment, an initial three-dimensional discrete element model of the target area is constructed based on geological information, including:

[0061] Based on the ductile joint information, the instructions for generating ductile structural surfaces for the initial three-dimensional discrete element model are determined, and based on the brittle joint information, the instructions for generating brittle structural surfaces for the initial three-dimensional discrete element model are determined.

[0062] Based on the instructions for generating ductile and brittle structural surfaces, corresponding ductile and brittle structural surface information is generated. Based on the ductile and brittle structural surface information, an initial three-dimensional discrete element model is constructed.

[0063] The ductile structure generation command refers to the parameterized control command for generating ductile structure surfaces from the initial 3D discrete element model. In practical scenarios, the ductile structure generation command is the fracture command. The brittle structure generation command refers to the parameterized control command for generating brittle structure surfaces from the initial 3D discrete element model. In practical scenarios, the brittle structure generation command is the joint command.

[0064] Among them, ductile structural surface information is used to represent the digital characterization information of ductile plating generated in the three-dimensional discrete element model; brittle structural surface information is used to represent the digital characterization information of brittle joints generated in the three-dimensional discrete element model.

[0065] Specifically, a zero-thickness geometric contact surface is generated according to the brittle structural surface generation command to simulate brittle joints, and a discrete, finite-size unit cell is generated according to the ductile structural surface generation command to simulate ductile joints, thus completing the construction of the initial three-dimensional discrete element model.

[0066] For example, the initial three-dimensional discrete element model can be a coupled model of fault fracture zone-tunnel excavation face.

[0067] In this embodiment, by extracting ductile plating information and brittle joint information from geological information and converting them into corresponding generation instructions, the initial three-dimensional discrete element model can be accurately constructed, enabling the model to fully reproduce the spatial distribution, occurrence characteristics and mechanical properties of ductile plating and brittle joints in the surrounding rock, effectively solving the problem of insufficient representation of the discontinuity of geological structures by traditional models.

[0068] In one embodiment, geostress equilibrium calculations are performed on the initial three-dimensional discrete element model to obtain a three-dimensional discrete element model, including:

[0069] In the overlying rock mass region of the initial three-dimensional discrete element model, and based on fluid conditions, boundary displacement conditions, and stress conditions, geostress balance calculations are performed to obtain the three-dimensional discrete element model.

[0070] The parameters of the rock mass area can be set according to the actual situation. In this embodiment, the coverage is 400 m and the density is 2600 kg / m³. 3 Let's take the rock mass as an example for illustration.

[0071] It should be noted that the fluid conditions, boundary displacement conditions, and stress conditions can all be set according to actual conditions. In this embodiment, the fluid property parameter is a fluid density of 1000 kg / m³. 3 The acceleration due to gravity is 9.81 m. 2 / s, fluid modulus 2 gigapascals, fluid permeability coefficient 10 -8 For example, the fluid conditions can include fluid flowing only in the fracture and the fluid participating in the geostress balance calculation; the boundary displacement conditions can include the normal velocity of each boundary being fixed at 0; the stress conditions can include applying a pore pressure of 4 MPa through the flow node in the fault region.

[0072] Specifically, by applying fluid conditions, boundary displacement conditions, and stress conditions to the initial three-dimensional discrete element model overlying the rock mass region, and then calculating until convergence to achieve geostress equilibrium, a three-dimensional discrete element model is obtained.

[0073] In this embodiment, by overlaying the rock mass region on the initial three-dimensional discrete element model, and performing geostress balance calculations based on fluid conditions, boundary displacement conditions, and stress conditions, the multi-factor coupling criteria of fault water pressure, internal stress, and rock mass failure mode are considered. This facilitates the subsequent accurate determination of the critical value of the minimum safe thickness for preventing mudslides and water inrushes, and solves the limitation of traditional continuous medium models in simulating the hydraulic fracturing of discontinuous rock masses. This provides an analytical and computational basis for the design and construction of high-risk tunnels with mudslides and water inrushes in plateau areas.

[0074] In one embodiment, based on tunnel information, a three-dimensional discrete element model is used to simulate the mudslide and water inrush disaster process caused by tunnel excavation, including:

[0075] The tunnel excavation was simulated using a three-dimensional discrete element model until the tunnel and fault region in the three-dimensional discrete element model were connected.

[0076] Specifically, the tunnel excavation is simulated using a three-dimensional discrete element model. After the three-dimensional discrete element model converges, the development of the fractures is observed. If the fractures from the tunnel face to the water-rich fault are not connected, the tunnel excavation continues to be simulated until the seepage channel is connected from the fault to the tunnel face. That is, at this time, a mudslide and water inrush disaster occurs.

[0077] In this embodiment, the dynamic process of plastic failure of surrounding rock and crack propagation and seepage channel formation driven by water pressure during tunnel construction is simulated. This facilitates subsequent quantitative assessment of the risk level of various factors in tunnel construction under different fault and tunnel spatial location conditions, thereby improving the accuracy of the simulation results.

[0078] To facilitate understanding by those skilled in the art, the numerical simulation method for tunnel mudslides and water inrushes is explained below with a specific example, such as... Figure 2 As shown, the geological information includes geological exploration data, and the initial three-dimensional discrete element model is a geological model, or simply the model.

[0079] Step S1: Based on the previously obtained geological survey data and specifications, establish a geological model containing faults. Use the joint command to generate a geometric contact surface with zero thickness to simulate brittle joints, and use the fracture command to generate discrete, finite-sized unit cells to simulate ductile plating.

[0080] Step S2: Group the model and assign corresponding physical properties to different parts of the model.

[0081] Step S3 involves applying conditions to the model and calculating until convergence to achieve geostress equilibrium.

[0082] Step S4: Execute the tunnel excavation command and observe the crack development after the model converges.

[0083] Step S5: Repeat step S4 until the seepage channel penetrates from the fault to the free face of the tunnel.

[0084] Step S6: Change the angle of intersection between the fault and the tunnel, and repeat steps S2-S5.

[0085] The preliminary geological survey data in step S1 mainly includes the lithology and dip angle of the surrounding rock in typical cross-sections of the ductile-brittle strong shear surrounding rock in the plateau, as well as the dip angle of the ductile joints and brittle joints, which are used to determine the extent of the common surrounding rock and faults and the occurrence of the ductile-brittle joints.

[0086] Specifically, the geological model can be 80*100*50 m in size, the tunnel diameter is 6*7 m, the fault angles are 45°, and the main attitudes of the ductile laminations are: one group dips at 200° to 240° with a dip angle of 45°, and another group dips at 270° to 300° with a dip angle of 30°, with lengths ranging from 1 m to 10 m; the main attitudes of the brittle joints in the ductile-brittle shear zone are: dip angle of 50°, dip direction of 230°, and spacing of 0.5 m; the geological model is as follows: Figure 3 As shown.

[0087] The grouping of the model in step S2 is to facilitate subsequent operations on the model. This includes grouping the tunnel and surrounding rock, grouping the fault and ordinary surrounding rock blocks, and grouping the fault and the contact surfaces between the rock masses in the fault. The physical parameters of the rock mass include the density, Poisson's ratio, Young's modulus, cohesion, internal friction angle, tensile strength, etc. of the surrounding rock. The physical parameters of the contact surface include normal stiffness, shear stiffness, and cohesion.

[0088] Specifically, in this embodiment, the density of the ordinary surrounding rock area is 2500 kg / m³, Young's modulus is 20 MPa, Poisson's ratio is 0.13, cohesion is 200 MPa, internal friction angle is 45 degrees, and tensile strength is 40 MPa; the density of the fault area is 2500 kg / m³, Young's modulus is 20 MPa, Poisson's ratio is 0.13, cohesion is 1 MPa, internal friction angle is 30 degrees, and tensile strength is 1 MPa. The contact normal stiffness between the blocks in the fault area is 2 GPa, the shear stiffness is 1 GPa, the internal friction angle is 20 degrees, and the cohesion is 0.25 MPa; the contact normal stiffness between the blocks in the ordinary surrounding rock is 20 GPa, the shear stiffness is 10 GPa, the internal friction angle is 25 degrees, and the cohesion is 0.5 MPa. By distinguishing between the tunnel and the surrounding rock, the tunnel excavation section is deleted and corresponding support structures are established.

[0089] The conditions applied in step S3 include boundary displacement conditions, stress conditions, and fluid conditions.

[0090] Specifically, in this embodiment, the normal velocity of each boundary is fixed at 0 (boundary displacement condition); the fluid density is 1000 kg / m³. 3 The acceleration due to gravity is 9.81 m. 2 / s, fluid modulus is 2 kPa, fluid permeability coefficient is 10 -8 (Fluid properties); The model overlying thickness is 400 m, and the density is 2600 kg / m³. 3 The rock mass. Before excavation, the command to allow fluid to flow only in the fractures and the command to allow fluid to participate in the calculation are activated. In the fault area, a pore pressure of 4 MPa (stress condition) is applied through the flow node to perform in-situ stress balance calculation. After balance is achieved, the displacement state is cleared.

[0091] The excavation command in step S4 includes the depth of cyclic excavation, the establishment of anchor bolts and lining.

[0092] Specifically, in this embodiment, the excavation depth is 2 m, the anchor bolt length is 4 m, there are 16 anchor bolts per ring, the longitudinal spacing is 1.5 m, the Young's modulus is 600,000 MPa, the tensile strength is 1,000,000 MPa, and the cross-sectional area is 1.57 cm². 2 The bonding strength is 2 million MPa, and the bonding stiffness is 2 million MPa; the lining thickness is 0.2 m, the coupling positive stiffness is 109 MPa, the coupling shear stiffness is 109 MPa, the coupling shear cohesion is 4*103 MPa, and the calculation convergence accuracy is 10. -5 .

[0093] Step S5 occurs because the fracture opening condition failed to penetrate from the water-rich fault to the tunnel free face in a certain calculation step. At this point, the next excavation step is executed. After the calculation converges, the tunnel displacement is observed. Figure 4 As shown, observe the development of the hydraulic aperture, such as Figure 5 As shown, observe the crack propagation status, such as Figure 6 As shown, if the fracture between the tunnel free face and the water-rich fault has been penetrated, it indicates that a mudslide and water inrush disaster has occurred. The safe thickness of the rock wall has reached the critical value in the previous excavation step. In this embodiment, when the distance between the tunnel free face and the fault is at least 23 m, the fracture has been penetrated and destroyed. The 6 MPa water pressure in the fault is also transmitted to the fractures around the tunnel.

[0094] Therefore, this application aims to leverage the advantages of the discrete element method to break through the continuous medium assumption of traditional numerical models, providing a high-precision numerical prediction method for the prevention and control of tunnel water inrush disasters under complex geological conditions, and providing theoretical basis and technical support for support design and safety assessment of similar projects.

[0095] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages in other steps. It is understood that the steps in different embodiments can be freely combined as needed, and all non-contradictory solutions formed by such combinations are within the scope of protection of this application.

[0096] Based on the same inventive concept, this application also provides a numerical simulation device for tunnel mudslide and water inrush disasters, used to implement the numerical simulation method for tunnel mudslide and water inrush disasters described above. The solution provided by this device is similar to the solution described in the above method. Therefore, the specific limitations of one or more embodiments of the numerical simulation device for tunnel mudslide and water inrush disasters provided below can be found in the limitations of the numerical simulation method for tunnel mudslide and water inrush disasters described above, and will not be repeated here.

[0097] In one exemplary embodiment, such as Figure 7 As shown, a numerical simulation device 700 for tunnel mudslide and water inrush disasters is provided. The device 700 includes:

[0098] The model acquisition module 701 is used to acquire a three-dimensional discrete element model of the target area based on the geological information of the target area; wherein, the three-dimensional discrete element model includes fault areas;

[0099] The simulation module 702 is used to simulate the mudslide and water inrush disaster process caused by tunnel excavation using a three-dimensional discrete element model based on tunnel information, and to obtain the numerical simulation results of the mudslide and water inrush disaster.

[0100] The simulation module 702 is also used to adjust the angle corresponding to the fault area crossed by the tunnel in the tunnel information, and based on the adjusted tunnel information, repeat the above steps of simulating the mudslide and water inrush disaster caused by tunnel excavation using a three-dimensional discrete element model based on the tunnel information, and obtain the numerical simulation results of the mudslide and water inrush disaster.

[0101] In one embodiment, the three-dimensional discrete element model also includes a surrounding rock region; the geological information includes physical parameters of the fault region, physical parameters of the rock mass in the surrounding rock region, and ductile plating information and brittle joint information of the surrounding rock region.

[0102] In one embodiment, the model acquisition module 701 is further configured to construct an initial three-dimensional discrete element model of the target area based on geological information;

[0103] The initial three-dimensional discrete element model is subjected to geostress equilibrium calculation to obtain the three-dimensional discrete element model.

[0104] In one embodiment, the model acquisition module 701 is further configured to determine, based on the ductile joint information, a ductile structural surface generation instruction for the initial three-dimensional discrete element model, and based on the brittle joint information, a brittle structural surface generation instruction for the initial three-dimensional discrete element model.

[0105] Based on the instructions for generating ductile and brittle structural surfaces, corresponding ductile and brittle structural surface information is generated. Based on the ductile and brittle structural surface information, an initial three-dimensional discrete element model is constructed.

[0106] In one embodiment, the model acquisition module 701 is further used to apply the initial three-dimensional discrete element model to the overlying rock mass region, and to perform geostress balance calculations based on fluid conditions, boundary displacement conditions, and stress conditions to obtain a three-dimensional discrete element model.

[0107] In one embodiment, the simulation module 702 is further configured to simulate tunnel excavation using a three-dimensional discrete element model until the tunnel and fault region in the three-dimensional discrete element model are interconnected.

[0108] The modules in the aforementioned numerical simulation device for tunnel mudslides and water inrushes can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the computer device's memory as software, so that the processor can call and execute the corresponding operations of each module.

[0109] In one exemplary embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 8 As shown, the computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, Near Field Communication (NFC), or other technologies. When the computer program is executed by the processor, it implements a numerical simulation method for tunnel mudslide and water inrush disasters. The display unit is used to form a visually visible image and can be a display screen, projection device, or virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0110] Those skilled in the art will understand that Figure 8 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0111] In one exemplary embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the above-described numerical simulation method for tunnel mudslide and water inrush disasters.

[0112] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the above-described numerical simulation method for tunnel mudslide and water inrush disasters.

[0113] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the above-described numerical simulation method for tunnel mudslide and water inrush disasters.

[0114] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0115] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, artificial intelligence (AI) processors, etc., and are not limited to these.

[0116] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0117] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A disaster numerical simulation method for tunnel sudden mud gushing, characterized in that, The method comprises: According to the geological information of the target area, a three-dimensional discrete element model of the target area is obtained; wherein the three-dimensional discrete element model comprises a fault area; Based on the tunnel information, the three-dimensional discrete element model is used to simulate the process of sudden mud gushing disaster caused by tunnel excavation, and the numerical simulation result of the sudden mud gushing disaster is obtained; Adjust the angle corresponding to the tunnel passing through the fault area in the tunnel information, and based on the adjusted tunnel information, repeat the above steps of based on the tunnel information, using the three-dimensional discrete element model to simulate the process of sudden mud gushing disaster caused by tunnel excavation, and obtaining the numerical simulation result of the sudden mud gushing disaster.

2. The method of claim 1, wherein, The three-dimensional discrete element model further comprises a surrounding rock area; the geological information comprises physical parameters of the fault area, rock mass physical parameters of the surrounding rock area, and ductile foliation information and brittle joint information of the surrounding rock area.

3. The method of claim 2, wherein, The three-dimensional discrete element model of the target area is obtained according to the geological information of the target area, comprising: According to the geological information, an initial three-dimensional discrete element model of the target area is constructed; The initial three-dimensional discrete element model is subjected to geostress balance calculation to obtain the three-dimensional discrete element model.

4. The method of claim 3, wherein, The initial three-dimensional discrete element model of the target area is constructed according to the geological information, comprising: According to the ductile foliation information, the ductile structure surface generation instruction for the initial three-dimensional discrete element model is determined, and according to the brittle joint information, the brittle structure surface generation instruction for the initial three-dimensional discrete element model is determined; According to the ductile structure surface generation instruction and the brittle structure surface generation instruction, the corresponding ductile structure surface information and brittle structure surface information are generated, and the initial three-dimensional discrete element model is constructed according to the ductile structure surface information and the brittle structure surface information.

5. The method of claim 3, wherein, The initial three-dimensional discrete element model is subjected to geostress balance calculation to obtain the three-dimensional discrete element model, comprising: The overlying rock mass area is calculated on the initial three-dimensional discrete element model, and the geostress balance calculation is carried out based on the fluid condition, the boundary displacement condition and the stress condition, to obtain the three-dimensional discrete element model.

6. The method of claim 1, wherein, The three-dimensional discrete element model is used to simulate the process of sudden mud gushing disaster caused by tunnel excavation based on the tunnel information, comprising: The three-dimensional discrete element model is used to simulate tunnel excavation until the tunnel in the three-dimensional discrete element model and the fault area are mutually penetrated.

7. A device for numerical simulation of disasters of tunnel sudden mud gushing, characterized in that, The device comprises: A model acquisition module is configured to obtain a three-dimensional discrete element model of a target area according to geological information of the target area; wherein the three-dimensional discrete element model comprises a fault area; A simulation module is configured to use the three-dimensional discrete element model to simulate the process of sudden mud gushing disaster caused by tunnel excavation based on tunnel information, and obtain the numerical simulation result of the sudden mud gushing disaster. The simulation module is further configured to adjust an angle corresponding to the tunnel passing through the fault region in the tunnel information, and repeatedly perform the step of simulating the sudden mud gushing disaster process caused by the tunnel excavation based on the tunnel information by using the three-dimensional discrete element model based on the adjusted tunnel information to obtain the numerical simulation result of the sudden mud gushing disaster.

8. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that, The computer program, when executed by the processor, implements the steps of the method of any one of claims 1 to 6.

9. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program, when executed by the processor, implements the steps of the method of any one of claims 1 to 6.

10. A computer program product comprising a computer program, characterized in that, The computer program, when executed by the processor, implements the steps of the method of any one of claims 1 to 6. The computer program, when executed by the processor, implements the steps of the method of any one of claims 1 to 6.

Images