A coal mine underground water flooding simulation method

Abstract

The application discloses a coal mine underground water flooding simulation method, which comprises the following steps: step 1, collecting roadway coordinates and floor elevations obtained by underground measurement, and representing underground connected roadways as a collection of multiple roadway sections; step 2, setting a water outlet position of water flooding simulation, determining a roadway section where the water outlet position is located, and determining a corresponding roadway section name, a starting point coordinate and a tail point coordinate; step 3, projecting a roadway model to obtain a water flow direction and a waterlogging position of water flooding simulation; and step 4, rotating a coordinate axis of the XZ or YZ plane roadway model in step 3 to obtain a ZX or ZY plane roadway model, calculating waterlogging volume by using the ZX or ZY plane roadway model in the rotated two-dimensional coordinate system, and simulating water flooding height. The application realizes quantitative water flooding simulation, and compared with an existing simulation water flooding model, the application realizes a data-driven model, and is more interpretable and universal.

Description

Technical Field

[0001] This invention belongs to the field of mine geological technology and relates to a method for simulating underground water flooding in coal mines. The method can simulate quantitative water outflow accidents in coal mines and achieve intelligent calculation of the flooding range and water accumulation. Background Technology

[0002] Coal mine flooding is one of the five major disasters in coal mines. Once it occurs, it endangers the lives of underground personnel, damages underground equipment, floods roadways, and causes huge economic losses. Current research focuses primarily on hydrological monitoring and flood hazard prediction and early warning to prevent flood accidents and ensure coal mine production safety. However, due to the complexity of hydrogeology, it cannot completely prevent flood accidents. When a flood occurs, the flooding process is particularly important for the escape of underground personnel. The undulating and complex roadways within mines mean that the unknown flooding process affects personnel escape routes. Existing water flow simulation methods mainly simulate the dynamics of water bodies, neglecting the influence of water volume on water level, making the calculated flooding simulation lack certain reliability.

[0003] In view of the above problems, there is an urgent need to study a new method for simulating well flooding that takes into account the influence of water volume on water level. Summary of the Invention

[0004] The purpose of this invention is to provide a method for simulating underground flooding in coal mines, so as to solve the problems of non-quantitative dynamic water simulation and the discrepancy between water level and actual flooding conditions.

[0005] To achieve the above objectives, the present invention employs the following technical solution:

[0006] A method for simulating underground water flooding in coal mines, specifically including the following steps:

[0007] Step 1: Collect the tunnel coordinates and floor elevation obtained from underground measurements, and represent the underground connecting tunnels as a set of multiple tunnel segments. The tunnel segments are numbered sequentially by natural numbers, and the tunnel segment name is Ni. Record the starting point three-dimensional coordinates (Xi, Yi, Zi) and ending point three-dimensional coordinates (Xi+1, Yi+1, Zi+1) of each tunnel segment according to the geodetic coordinate system, where i represents the tunnel segment number, which is a natural number;

[0008] Step 2: Set the water outlet location in the flood simulation, determine the roadway segment where the water outlet location is located in the underground roadway topology structure of Step 1, and determine the corresponding roadway segment name, start point coordinates, and end point coordinates;

[0009] Step 3: Project the tunnel model to obtain the water flow direction and water accumulation location in the flood simulation, including the following sub-steps:

[0010] Step 31: Perform a two-dimensional projection on the underground roadway topology structure of Step 1 to obtain an XZ plane roadway model or a YZ plane roadway model. In the obtained XZ plane roadway model or YZ plane roadway model, number the positions in the roadway in chronological order.

[0011] Step 32: Construct two monotonically decreasing stacks to determine the elevation Z of the road surface. Utilize the characteristics of the monotonically decreasing stacks (last-in, first-out) and the monotonically decreasing nature of the elements within the stack to determine the direction of water flow and the location of water accumulation.

[0012] Step 4: Rotate the coordinate axes of the XZ or YZ plane tunnel model in Step 3 to obtain the ZX or ZY plane tunnel model. Use the rotated ZX or ZY plane tunnel model in the two-dimensional coordinate system to calculate the water accumulation and simulate the flooding height.

[0013] Furthermore, step 32 specifically includes the following process:

[0014] (1) Take the starting point number and ending point number of the tunnel section where the water outlet location is set in step 2 as the bottom elements of two monotonically decreasing stacks respectively. When there is only one element in the stack, the top element of the stack is the bottom element of the stack. The two stacks are called stack A and stack B.

[0015] (2) Determine the tunnel elevation Z corresponding to the top element of stack A and stack B, and determine the water flow direction of the tunnel by the tunnel elevation Z corresponding to the start point number and the end point number. Use the stack with the smaller tunnel elevation Z as the current execution stack.

[0016] (3) Traverse the node numbers in the tunnel along the direction of water flow, and control the node to enter the current execution stack according to the node elevation. Specifically: if the elevation Z of the currently traversed node is less than the elevation Z of the node corresponding to the current top element of the current execution stack, then the current node number is pushed onto the stack, indicating that the water flow direction is the same; otherwise, it indicates that the tunnel has formed a "concave" state, that is, water accumulation has formed. At this time, the current top element of the stack is popped from the stack and marked as a water accumulation position.

[0017] (4) After the current top element is popped from the stack, if the current execution stack is empty, that is, all elements have been popped from the stack, it means that the entire section of the tunnel has been submerged by the water flow. At this time, the current node is put into the current execution stack as the bottom element of the stack, the current execution stack ends its work, and returns to step (2); if the current execution stack is not empty, then continue to traverse all nodes in the direction of the water flow in the tunnel determined in step (2), the current execution stack ends its work, and returns to step (2), until all nodes corresponding to the position numbers in the tunnel have been traversed.

[0018] Furthermore, step 4 specifically includes the following sub-steps:

[0019] Step 41: For two consecutive roadway segments N1 and N2, the starting and ending positions of roadway segment N1 are numbered 1 and 2, and the starting and ending positions of roadway segment N2 are numbered 2 and 3. When the elevation Z2 of the ending point of roadway segment N1 is less than the elevation Z1 of the starting point of roadway segment N1, the elevation Z2 of the starting point of roadway segment N2 is less than the elevation Z3 of the ending point of roadway segment N2, and the elevation Z3 of the ending point of roadway segment N2 is less than the elevation Z1 of the starting point of roadway segment N1, water accumulation occurs at position 2. It means that during the water flow process, when the water surface on the roadway rises from position 2 in roadway segment N2 to position 3, the water flow plane intersects with roadway segment N1, and water accumulation occurs at this point. The polygon formed by roadway segments N1, N2, and the water flow plane is the two-dimensional representation of the water accumulation volume at this point. Rotate the XZ or YZ plane roadway model in step 3, and linearly represent each roadway segment Ni in step 1 as f in the rotated ZX or ZY model Ni , the elevation of each position i in the roadway is represented as Zi. Then, when the water flow rises from position 2 to position 3, the water accumulation volume q at this point in the two-dimensional plane is

[0020]

[0021] where the integral intervals Z2 and Z3 respectively represent the elevation of position 2 and the elevation of position 3 in roadway segment N2, and f N1 represents the linear representation of roadway segment N1 in the ZX or ZY plane roadway model, and f N2 represents the linear representation of roadway segment N2 in the ZX or ZY plane roadway model;

[0022] Step 42: Calculate the water accumulation volume q at the water accumulation positions marked in step 3 in sequence according to the water flow submergence path obtained in step 3. For each calculated water accumulation volume at a water accumulation position, calculate the cumulative value of the water accumulation volume count = ∑q until the calculation stops when Q < count, and mark the next water accumulation position. Q is the water discharge set according to the need of the water flooding simulation. At this time, the cumulative value of the current water accumulation volume exceeds the set water discharge, and this quantitative water flooding simulation ends. Assume that at this time, the water surface intersects with roadway segments Ni and Nj respectively. The starting and ending points of roadway segment Ni are positions n and l, the starting and ending points of roadway segment Nj are positions h and g, and the roadway elevations satisfy Zl < Zh < Zn < Zg. When the water surface is static, it intersects with roadway segment Ni at position m, and calculate the water surface elevation Zm at the intersection point position m

[0023]

[0024] where the integral intervals Zm and Zn respectively represent the elevation of the intersection point position m of the water surface and roadway segment Ni and the elevation of the starting point of roadway segment Ni, and f Ni represents the linear representation of roadway segment Ni in the ZX or ZY plane roadway model, and f Nj represents the linear representation of roadway segment Nj in the ZX or ZY plane roadway model;

[0025] Step 43: Zm is the water surface elevation of roadway segment Ni when the outflow rate is Q. Based on the elevation difference between roadway segments, the water accumulation height at other water accumulation locations within the roadway is calculated. Specifically, in connected roadway segments, the water accumulation height at any water accumulation location p is Zm - Zp; the water accumulation height at any water accumulation location c in non-connected roadway segments is Zm - Zp. 最高 -Zc.

[0026] This invention first constructs a three-dimensional tunnel model that closely approximates reality based on realistic underground data. Then, it maps the simulated water seepage points in the tunnel model to the actual tunnel data. Through data comparison and calculation, it achieves quantitative flooding simulation. Compared to existing flooding simulation models, this invention implements a data-driven model, which is more interpretable and universally applicable.

[0027] This invention can calculate the quantitative water flow spread path and the water level at various locations along the water flow spread path in the roadway after flooding simulation. It can also provide accurate data support for coal mine flooding simulation to achieve a more realistic and evidence-based simulation effect.

[0028] This invention establishes a quantitative underground water flooding simulation model and studies water flooding simulation based on a static roadway model, providing more data support for underground water flooding visualization and thus providing more reliable data for underground personnel to avoid danger. Attached Figure Description

[0029] To more clearly illustrate the embodiments of the present invention, the accompanying drawings used in the present invention will be briefly described below.

[0030] Figure 1 This is a flowchart of the coal mine underground water flooding simulation method of the present invention;

[0031] Figure 2 This is a topological diagram of the actual tunnel section of the present invention;

[0032] Figure 3 This is a two-dimensional projection schematic diagram of the actual tunnel section of the present invention;

[0033] Figure 4 This is a schematic diagram showing the water level at the location of water accumulation within the tunnel section. Detailed Implementation

[0034] This invention provides a method for simulating underground water flooding in coal mines. By calling up actual roadway data, setting the water output and water output location, and inputting the algorithm model constructed by the method of this invention, the user can simulate and calculate the underground water flooding process based on the actual roadway, analyze and calculate the water flow direction, water accumulation location, water flooding path, and water level line of the water accumulation area, and display the water flooding situation of the actual roadway.

[0035] like Figure 1 As shown, the coal mine underground water flooding simulation method of the present invention includes the following steps:

[0036] Step 1: Collect the tunnel coordinates and floor elevation obtained from underground surveys, and represent the connected underground tunnels as a set of multiple tunnel segments. Tunnel segments are numbered sequentially by natural numbers, and the segment name is Ni. Record the starting point's three-dimensional coordinates (Xi, Yi, Zi) and ending point's three-dimensional coordinates (Xi+1, Yi+1, Zi+1) of each tunnel segment according to the geodetic coordinate system, where i represents the tunnel segment number, which is a natural number. The results are summarized in the table below.

[0037] Then, the recorded tunnel segments are connected by their starting and ending points according to their coordinates to form the tunnel's topology.

[0038]

[0039] Step 2: Determine the location of the water outlet in the underground tunnel.

[0040] Before conducting the flooding simulation, the water outlet location is set. The corresponding underground roadway topology obtained in step 1 is used to find the roadway segment where the water outlet location is located, and the corresponding roadway segment name, starting point coordinates, and ending point coordinates are obtained.

[0041] Step 3: Project the tunnel model to simulate flooding.

[0042] Step 31: Perform a two-dimensional projection on the underground roadway topology structure of Step 1 to obtain an XZ plane roadway model or a YZ plane roadway model. In the obtained XZ plane roadway model or YZ plane roadway model, number each location (starting point or ending point of the roadway segment) in the roadway according to the order of occurrence.

[0043] Step 32: Construct two monotonically decreasing stacks to determine the road surface elevation Z. Utilizing the LIFO (Last-In, First-Out) characteristic of monotonically decreasing stacks and the monotonically decreasing nature of their elements, the direction of water flow and the location of water accumulation are determined, thus simulating flooding. The specific process includes the following steps:

[0044] (1) Take the starting point number and ending point number of the tunnel section where the water outlet location is set in step 2 as the bottom elements of two monotonically decreasing stacks (when there is only one element in the stack, the top element of the stack is the bottom element of the stack). The two stacks are called stack A and stack B (or left stack and right stack).

[0045] (2) Determine the tunnel elevation Z corresponding to the top element of stack A and stack B, and determine the water flow direction of the tunnel by the tunnel elevation Z corresponding to the start point number and the end point number. Use the stack with the smaller tunnel elevation Z as the current execution stack.

[0046] (3) Traverse the node numbers in the tunnel along the direction of water flow, and control the node to enter the current execution stack according to the node elevation. Specifically: if the elevation Z of the currently traversed node is less than the elevation Z of the node corresponding to the current top element of the current execution stack, then the current node number is pushed onto the stack, indicating that the water flow direction is the same; otherwise, it indicates that the tunnel has formed a "concave" state, that is, water accumulation has formed. At this time, the current top element of the stack is popped from the stack and marked as a water accumulation position.

[0047] (4) After the current top element is popped from the stack, if the current execution stack is empty, that is, all elements have been popped from the stack, it means that the entire section of the tunnel has been submerged by the water flow. At this time, the current node is put into the current execution stack as the bottom element of the stack, the current execution stack ends its work, and returns to step (2); if the current execution stack is not empty, then continue to traverse all nodes in the direction of the water flow in the tunnel determined in step (2), the current execution stack ends its work, and returns to step (2), until all nodes corresponding to the position numbers in the tunnel have been traversed.

[0048] In the above technical solution, during the judgment process, the popped element of the monotonically decreasing stack represents the water accumulation position in the flood simulation. In step (3), if the current execution stack changes from stack A to stack B or from stack B to stack A, it indicates that the water flow direction in the tunnel has changed. The set of all paths from the set water outlet position to all water accumulation positions is the entire range of the water flow flooding path.

[0049] Step 4: Rotate the coordinate axes of the XZ or YZ plane tunnel model from Step 3 to obtain the ZX or ZY plane tunnel model. Use the rotated ZX or ZY plane tunnel model in the two-dimensional coordinate system to calculate the water accumulation and simulate the flooding height. The calculation process is as follows:

[0050] Step 41: When the tunnel forms a concave shape, water accumulation occurs. Taking two consecutive tunnel segments N1 and N2 as examples, the starting and ending points of tunnel segment N1 are numbered 1 and 2, and the starting and ending points of tunnel segment N2 are numbered 2 and 3. When the elevation Z2 of the ending point of tunnel segment N1 is lower than the elevation Z1 of the starting point of tunnel segment N1, the elevation Z2 of the starting point of tunnel segment N2 is lower than the elevation Z3 of the ending point of tunnel segment N2, and the elevation Z3 of the ending point of tunnel segment N2 is lower than the elevation Z1 of the starting point of tunnel segment N1, water accumulation occurs at position 2. This indicates that during the water flow process, when the water in the tunnel rises from position 2 to position 3 in tunnel segment N2, the water flow plane intersects with tunnel segment N1, and water accumulation occurs at this point. The polygon formed by tunnel segments N1 and N2 and the water flow plane is the two-dimensional representation of the water accumulation volume. Rotate the XZ or YZ plane tunnel model in step 3, and linearly represent each tunnel segment Ni from step 1 in the ZX or ZY model obtained after rotation as f. Ni Let Zi be the elevation of each position i in the tunnel. Then, when the water flow rises from position 2 to position 3, the water accumulation at that point in the two-dimensional plane is q, which can be calculated by definite integral, as follows:

[0051]

[0052] Among them, the integral intervals Z2 and Z3 respectively represent the elevations of position 2 and position 3 in the roadway section N2, and f N1 represents the linear representation of the roadway section N1 in the roadway model of the ZX or ZY plane, and f N2 represents the linear representation of the roadway section N2 in the roadway model of the ZX or ZY plane.

[0053] Step 42: Calculate the accumulated water volume q of the waterlogging positions marked in Step 3 in sequence according to the water flow submergence path obtained in Step 3. The accumulated water volume q represents the accumulated water volume of the roadway section when a certain position in the roadway to another position is "completely" submerged by water. Each time the accumulated water volume of a waterlogging position is obtained, the accumulated value count = ∑q is calculated. The calculation stops and the next waterlogging position is marked until Q < count. Q is the water output set according to the needs of the waterlogging simulation; at this time, the accumulated value of the current accumulated water volume exceeds the set water output, and this quantitative waterlogging simulation ends.

[0054] Assume that at this time the water surface intersects the roadway sections Ni and Nj respectively. The starting point and the end point of the roadway section Ni are position n and position l, the starting point and the end point of the roadway section Nj are position h and position g, and the roadway elevations are Zl < Zh < Zn < Zg. When the water surface is static, it intersects the roadway section Ni at position m, indicating that the end point position l of the roadway section Ni has been submerged by water, the starting point position n of the roadway section Ni has not been submerged by water, and the roadway section Ni is "not completely" submerged by water. Then when calculating the cumulative value of count at this time, only the known roadway positions l, position h, and position n are calculated, and it is considered that all calculated positions are submerged by the water flow, ignoring the state of the roadway section being "not completely" submerged. Therefore, the following formula is used for correction.

[0055] During the accumulation process, when the water surface is static, it intersects the roadway section Ni at position m, and the starting point position n of the roadway section Ni is not submerged by the water flow. The area not submerged by the water flow is expressed as count - Q and is not included in the actual waterlogging simulation waterlogging position. Then the actual non-waterlogged area can be expressed as:

[0056]

[0057] Among them, the integral intervals Zm and Zn respectively represent the elevation of the intersection position m of the water surface and the roadway section Ni and the starting elevation of the roadway section Ni, and f Ni represents the linear representation of the roadway section Ni in the roadway model of the ZX or ZY plane, and f Nj represents the linear representation of the roadway section Nj in the roadway model of the ZX or ZY plane.

[0058] By setting the water output Q, the cumulative value of water accumulation location count, and the linear representation of the planar tunnel model, the horizontal elevation Zm of the intersection location m is calculated using indefinite integral.

[0059] Step 43: Zm is the water surface elevation of roadway section Ni when the water output is Q. The water accumulation height of other water accumulation locations in the roadway is calculated based on the height difference of the roadway sections.

[0060] The approach to calculating the water level at other locations is as follows: (e.g.) Figure 4 As shown, after the quantitative flooding simulation is completed, based on the flooding path, assuming the recorded water flow starts from position d in the tunnel, after the quantitative flooding simulation, the horizontal plane intersects with tunnel segment Ni at position m, forming multiple water accumulation positions i, j, k, l, and h in the tunnel. This generates water accumulation area 1 and water accumulation area 2. The water level height at each water accumulation position in each water accumulation area within the tunnel can be calculated based on the last water accumulation position h after the flooding simulation. Details are as follows:

[0061] In flooded area 2, the water-filled positions (e.g., l) are connected to the last water-filled position h; therefore, flooded area 2 is called the connected roadway segment. In flooded area 1, the water-filled positions (e.g., i, j, k) are not connected to the last water-filled position h; therefore, flooded area 1 is called the non-connected roadway segment. In the connected roadway segment (flooded area 2), the water height at position l is Zm-Zl, and the water height at position h is Zm-Zh. Therefore, it can be deduced that the water height at any water-filled position p in the connected roadway segment is Zm-Zp. In the non-connected roadway segment (flooded area 1), since flooded area 1 is submerged earlier than flooded area 2, the water level in flooded area 1 is independent of the water-filled position m, but is located precisely at the position (k) with the maximum elevation in the connected water-filled sequence (i, j, k); this position is called Z. 最高 At this point, the water level at location k is 0, the water level at location j is Zk-Zj, and the water level at location i is Zk-Zi. Therefore, it can be deduced that the water level at any location c in a non-connected tunnel section is Zk. 最高 -Zc.

[0062] Using this method, the water level values ​​at each water accumulation location in all roadway sections are calculated sequentially, and then flooding simulation is performed based on all water accumulation locations and their water level heights.

[0063] It is more interpretable and can provide accurate data support for coal mine flooding simulation, so as to achieve more realistic and evidence-based simulation results.

[0064] Example:

[0065] Step 1: Collect the coordinates and floor elevation data of the working face roadway obtained from underground surveying, and organize them as shown in the table below. Connect the starting and ending points of the recorded roadway segments according to their coordinates to form a complete underground roadway topology. The working face roadway topology is as follows: Figure 2 As shown.

[0066]

[0067]

[0068]

[0069] Step 2: Set the coordinates of the water outlet location in the roadway as (3923580.526, 36584068.74, 942.469). By comparing with the coordinates of the starting and ending points of the roadway segment, the water outlet location is determined to be on the roadway segment named *** Working Face 6. The starting coordinates of this roadway segment are (3923573.302, 36584029.203, 945.551), and the ending coordinates of this roadway segment are (3923584.678, 36584147.34, 941.717).

[0070] Step 3: Project the tunnel model to simulate flooding.

[0071] The underground roadway topology obtained in Step 1 is projected in two dimensions to obtain a YZ plane roadway model, which is then displayed. The roadway elevation Z at each measurement location (i.e., each node) in the two-dimensional projection of the roadway model is numbered sequentially (starting from 0), as follows: Figure 3 As shown, the water outlet position set in step two is... Figure 3 Between nodes 5 and 6. Next, the elevation Z of the road surface is determined by constructing two monotonically decreasing stacks. Utilizing the LIFO (Last-In, First-Out) property of the monotonically decreasing stacks and the monotonically decreasing nature of its elements, the direction of water flow and the location of water accumulation are determined, thus simulating flooding. The process of using the monotonically decreasing stack method is described below:

[0072] (1) The elevation Z of each measured position on the road surface is used to determine the starting point and ending point (node ​​5 and node 6) of the road section where the water outlet is located. These are respectively used as the bottom of two monotonic stacks, denoted as the left stack and the right stack. The left stack uses node 5 as the bottom and the right stack uses node 6 as the bottom. At this time, there is only one node in each stack. Therefore, the bottom node of the stack is also the top node of the stack.

[0073] (2) The elevation of node 5 is greater than that of node 6, so the water flows towards node 6 and the right stack becomes the current execution stack.

[0074] (3) Compare the elevation of node 7 along the direction of water flow. If it is less than the elevation of node 6 at the top of the right stack, then node 7 is put into the right stack. Subsequently, the elevations of nodes 8, 9, 10, and 11 are all less than the elevation of node 6 at the current top of the right stack. They are pushed onto the stack in sequence. When traversing node 12, the current top of the right stack is 11. The elevation of node 12 is greater than the elevation of node 11. Node 11 is popped from the stack. Node 11 is the first water accumulation position.

[0075] (4) At this time, the top of the right stack is node 10. The elevation of node 12 is less than the elevation of node 10, so push 12 onto the stack. The elevation of node 13 is greater than the elevation of node 12, so pop node 12 from the stack. Then node 12 is the second water accumulation point.

[0076] (5) When traversing to node 17 in this way, the node elevation continuously decreases, so the node number is continuously pushed onto the stack. When traversing to node 24, the elevation of node 24 is greater than that of node 23, so node 23 is popped from the stack as the water accumulation position. At this time, the top of the stack is node 22.

[0077] (6) Traverse all nodes in the manner described above;

[0078] (7) Popping nodes 11, 12, 10, 13, 14, 9, 15, 16, 23, 22, 24, 21, and 25 are the water accumulation locations that appear sequentially in the flood simulation. The water flows from the water outlet location set in step two towards node 25.

[0079] Step 4: Transform the coordinate system and calculate the water volume.

[0080] Specifically: Figure 3 The coordinate axes of the YZ plane tunnel model are rotated 90° to the right to obtain the ZY plane tunnel model. The ZX or ZY plane tunnel model in the rotated two-dimensional coordinate system is then used to calculate the water accumulation along the horizontal line. The calculation process is as follows:

[0081] (1) As can be seen from step three, the first water accumulation position is 11 and the second water accumulation position is 12. When the water flow height rises from position 11 to position 12, the water accumulation volume of this roadway section can be calculated by definite integral simulation, as shown below:

[0082]

[0083] Wherein, the integration intervals Z11 and Z12 represent the elevations of water accumulation positions 11 and 12, respectively, f 11 f represents the linear representation of working face 11 in the ZY plane roadway model. 12 This represents the linear representation of working face 12 of roadway segment *** in the ZY plane roadway model.

[0084] (2) Set the total outflow Q to 50000 m³. 3When water accumulation location 24 is popped from the stack, the cumulative water accumulation count is greater than Q, indicating that the flooding simulation for that volume has ended. This signifies that the water flow is continuously and dynamically moving from the outlet location to location 24, thus ending the flooding and bringing the water level to a standstill. When the water level is at a standstill, it intersects with the *** working face 22 and *** working face 25 of the roadway section. Assuming the water level elevation is x, the un-waterlogged area above the water surface up to the next water accumulation location 21 is represented by count - Q, as shown in the formula below:

[0085]

[0086] Where the integration intervals Zx and Z21 represent the elevation of the water level and the elevation of the next water accumulation location, respectively, f 22 and f 25 This represents the linear representation of two intersecting roadway segments on the ZY plane in the roadway model. Using indefinite integrals, the water surface elevation Zx at the last water accumulation location at the end of the quantitative flooding simulation is calculated to be 922.962.

[0087] (3) The last water accumulation location is the tail point of *** working face 24, with a water surface elevation of 922.962. The water accumulation height at this location is represented as Zx-Z24, which means that the water accumulation height at location 24 of the tunnel road surface after quantitative flooding simulation is 2.15. Similarly, based on the height difference of the tunnel road surface, the water accumulation heights of other connected water accumulation locations are obtained. The locations connected to location 24 include locations 23, 25, 22, and 21, and the calculated water accumulation heights are 4.024, 0.373, 3.965, and 0.913, respectively. Positions 11, 12, 10, 13, 14, 9, 15, and 16 are not connected to the last water accumulation position 24 in the water accumulation position sequence in step three. Therefore, they are calculated from position 17, where the maximum elevation of the connected water accumulation positions is located, and are respectively 16.741, 12.238, 11.286, 5.473, 4.87, 4.482, 4.305, and 2.332.

[0088] By using specific water depth data at various locations on the road surface of mine shafts and tunnels for flood simulation, the results are more interpretable and can provide accurate data support for coal mine flood simulation, thereby achieving more realistic and evidence-based simulation effects.

Claims

1. A method for simulating underground water flooding in coal mines, characterized in that, Specifically, the following steps are included: Step 1: Collect the tunnel coordinates and floor elevation obtained from underground measurements, and represent the underground connecting tunnels as a set of multiple tunnel segments; the tunnel segments are numbered in natural number order, and the tunnel segment name is Ni. Record the starting point three-dimensional coordinates (Xi, Yi, Zi) and ending point three-dimensional coordinates (Xi+1, Yi+1, Zi+1) of each tunnel segment according to the geodetic coordinate system, where i represents the tunnel segment number, which is a natural number; Step 2: Set the water outlet location in the flood simulation, determine the roadway segment where the water outlet location is located in the underground roadway topology structure of Step 1, and determine the corresponding roadway segment name, start point coordinates, and end point coordinates; Step 3: Project the tunnel model to obtain the water flow direction and water accumulation location in the flood simulation, including the following sub-steps: Step 31: Perform a two-dimensional projection on the underground roadway topology structure of Step 1 to obtain an XZ plane roadway model or a YZ plane roadway model. In the obtained XZ plane roadway model or YZ plane roadway model, number the positions in the roadway in chronological order. Step 32: Construct two monotonically decreasing stacks to determine the elevation Z of the road surface. Utilize the LIFO (Last-In, First-Out) property of the monotonically decreasing stacks and the monotonically decreasing nature of the elements within the stack to determine the direction of water flow and the location of water accumulation. The specific process includes the following steps: (1) Take the starting point number and ending point number of the tunnel section where the water outlet location is set in step 2 as the bottom elements of two monotonically decreasing stacks respectively. When there is only one element in the stack, the top element of the stack is the bottom element of the stack. The two stacks are called stack A and stack B. (2) Determine the tunnel elevation Z corresponding to the top element of stack A and stack B, and determine the water flow direction of the tunnel by the tunnel elevation Z corresponding to the starting point number and the ending point number. Take the stack with the smaller tunnel elevation Z as the current execution stack. (3) Traverse the node numbers in the tunnel along the direction of water flow, and control the node to enter the current execution stack according to the node elevation. Specifically: if the elevation Z of the currently traversed node is less than the elevation Z of the node corresponding to the current top element of the current execution stack, then the current node number is pushed onto the stack, indicating that the water flow direction is the same; otherwise, it indicates that the tunnel has formed a "concave" state, that is, water accumulation has formed. At this time, the current top element of the stack is popped from the stack and marked as a water accumulation position. (4) After the current top element is popped from the stack, if the current execution stack is empty, that is, all elements have been popped from the stack, it means that the entire section of the tunnel has been submerged by the water flow. At this time, the current node is pushed into the current execution stack as the bottom element of the stack, the current execution stack ends its work, and returns to step (2); If the current execution stack is not empty, then continue to traverse all nodes in the direction of the water flow in the tunnel determined in step (2), the current execution stack ends its work, and returns to step (2), until all nodes corresponding to the position numbers in the tunnel have been traversed; Step 4: Rotate the coordinate axes of the XZ or YZ plane tunnel model in Step 3 to obtain the ZX or ZY plane tunnel model. Use the rotated ZX or ZY plane tunnel model in the two-dimensional coordinate system to calculate the water accumulation and simulate the flooding height.

2. The method for simulating underground water flooding in coal mines as described in claim 1, characterized in that, Step 4 specifically includes the following sub-steps: Step 41: For two consecutive tunnel segments N1 and N2, the starting and ending points of tunnel segment N1 are numbered 1 and 2, and the starting and ending points of tunnel segment N2 are numbered 2 and 3. When the elevation Z2 of the ending point of tunnel segment N1 is less than the elevation Z1 of the starting point of tunnel segment N1, the elevation Z2 of the starting point of tunnel segment N2 is less than the elevation Z3 of the ending point of tunnel segment N2, and the elevation Z3 of the ending point of tunnel segment N2 is less than the elevation Z1 of the starting point of tunnel segment N1, water accumulation occurs at position 2. This indicates that during the water flow process, when the water in the tunnel rises from position 2 to position 3 in tunnel segment N2, the water flow plane intersects with tunnel segment N1, and water accumulation occurs at this point. The polygon formed by tunnel segments N1, N2, and the water flow plane is the two-dimensional representation of the water accumulation at this point. Rotate the XZ or YZ plane tunnel model in step 3, and linearly represent each tunnel segment Ni in step 1 in the ZX or ZY model obtained after rotation as follows: Let Zi be the elevation of each position i in the tunnel. Then, when the water flow rises from position 2 to position 3, the water accumulation q at that point in the two-dimensional plane is: Wherein, the integration intervals Z2 and Z3 represent the elevations of position 2 and position 3 in roadway segment N2, respectively. This represents the linear representation of roadway segment N1 in the ZX or ZY plane roadway model. This represents the linear representation of tunnel segment N2 in the ZX or ZY plane tunnel model. Step 42: Calculate the accumulated water volume q of the waterlogging positions marked in Step 3 in sequence according to the water flow inundation path obtained in Step 3. Each time the accumulated water volume of a waterlogging position is obtained, calculate the accumulated value count of the accumulated water volume , and stop calculating and marking the next waterlogging position until Q < count, where Q is the water discharge set according to the waterlogging simulation requirement; at this time, the accumulated value of the current accumulated water volume exceeds the set water discharge, and the waterlogging simulation ends; assume that the water surface intersects the roadway sections Ni and Nj at this time. The starting point and the end point of the roadway section Ni are positions n and l, the starting point and the end point of the roadway section Nj are positions h and g, and the roadway elevations are Zl < Zh < Zn < Zg. When the water surface is static, it intersects the roadway section Ni at position m, and calculate the water surface elevation Zm of the intersection position m Wherein, the integration intervals Zm and Zn represent the elevation of the intersection point m between the horizontal plane and the roadway segment Ni, and the starting elevation of the roadway segment Ni, respectively. This represents the linear representation of roadway segment Ni in the ZX or ZY plane roadway model. This represents the linear representation of tunnel segment Nj in the ZX or ZY plane tunnel model. Step 43: Zm is the water surface elevation of roadway segment Ni when the outflow rate is Q. Based on the elevation difference between roadway segments, the water accumulation height at other water accumulation locations within the roadway is calculated. Specifically, in connected roadway segments, the water accumulation height at any water accumulation location p is Zm - Zp; the water accumulation height at any water accumulation location c in non-connected roadway segments is Zm - Zp. 最高 -Zc.

Images