Method for predicting surface displacement and deformation of mining areas based on the division method

Abstract

To address the problem of the lack of means to accurately predict the surface displacement deformation of mining areas in complex boundary-shaped working surfaces, and the difficulty of performing accurate calculations, this invention provides a mining area surface displacement deformation prediction method based on a segmentation method that improves the efficiency and accuracy of mining area surface displacement deformation prediction. [Solution] The method of the present invention first identifies the location information of the mining site boundary in the mining area and obtains the coordinate information of the mining site boundary, establishes a function equation in the horizontal plane rectangular coordinate system of the mining site boundary based on the coordinate information, divides the mining site based on whether the boundary of the mining site can be represented by the same function equation, calculates the ground settlement value occurring in each segment of the mining site using the equation, superimposes the settlement values ​​at the same ground position for all segments of the mining site to obtain the complete ground settlement value, and finally obtains other parameters of the mining area ground movement deformation based on the relationship between other ground movement deformation parameters and the ground settlement value.

Description

[Technical Field]

[0001] The present invention relates to the technology of ground surface displacement deformation prediction, and more specifically to a method for predicting ground surface displacement deformation of mining areas based on a division method. [Background technology]

[0002] Coal seam mining causes a certain degree of displacement and deformation of the ground surface, negatively impacting the safety and stability of the ground surface morphology and related buildings (structures). Accurately predicting the extent of displacement and deformation that occurs on the ground surface after coal mining is crucial for preventing mining damage and protecting the environment. Currently, the most commonly used method for predicting surface displacement and deformation in mining areas is the stochastic integral method. The principle of this method is to consider the geological layers covering the mining site as multiple small spheres, and to use the probability of a sphere falling onto the mining site as the corresponding ground settlement value. This method allows for relatively accurate determination of the ground displacement and deformation.

[0003] Currently, there are a relatively large number of research reports applying stochastic integral methods to predict surface displacement deformation in rectangular work area mining, but errors or relatively large errors exist when predicting surface displacement deformation in complex-shaped work area mining. One of the main reasons is that, in conventional techniques, when predicting surface displacement deformation using stochastic integral methods, researchers often ignore the relationship between two orthogonal integral variables in the process of applying it from two dimensions (longitudinal section of strata) to three dimensions, and it is fixedly assumed that there is no relationship between the two integral direction variables. While this does not affect the prediction results when solving rectangular work area mining, it causes errors when predicting surface deformation in non-rectangular work area mining, making prediction errors likely when solving the prediction of surface displacement deformation in mining areas with complex boundary shapes using this method. The second main reason is that predicting surface displacement deformation in mining areas with complex boundary shapes needs to be done according to a certain method, but currently there is no efficient and accurate prediction method.

[0004] This indicates that improvements to conventional technologies are needed. [Overview of the project] [Problems that the invention aims to solve]

[0005] To address the problem of the lack of means to accurately predict the surface displacement deformation of mining areas in complex boundary-shaped work surfaces and the difficulty of performing accurate calculations, the present invention provides a mining area surface displacement deformation prediction method based on a segmentation method that improves the efficiency and accuracy of mining area surface displacement deformation prediction. [Means for solving the problem]

[0006] The technical proposal of this invention is as follows:

[0007] A method for predicting surface displacement and deformation of a mining area based on a division method, comprising the following steps:

[0008] a. Identify the location information of the mining site boundary in the mining area, obtain the coordinate information of the mining site boundary, and establish a functional equation in the horizontal plane rectangular coordinate system of the mining site boundary based on the said coordinate information. b. The mining area is divided based on whether the boundaries of the mining area can be expressed by the same functional equation. c. For each segment after division, use equation (1) to calculate and obtain the ground settlement value that will occur at the mining site of this segment.

[0009] JPEG2026076913000002.jpg25170

[0010] In equation (1), W(x,y) is the settlement value of the mining area surface (x,y), and W cm ―Maximum settlement value of the mining area surface, a known parameter, H―Reserve depth corresponding perpendicularly to the coal seam (or mining site) at the location of the mining area surface, tanβ―Main influencing angle tangent, t i-1 -Mining site boundary P i-1 The t-coordinate of the point, i -Mining site boundary P iThe t-coordinate of the point is x, the t-coordinate of the mining area surface point is y, the ρ-coordinate of the mining area surface point is ρ and t, the ρ and t directions constitute a plane rectangular coordinate system, π is pi, e is the base of the natural logarithm, σ is the integral variable, erf(p) is the stochastic integral function where p represents a variable, f2(t) is the functional equation of the mining boundary two, ρ2=f2(t), and f1(t) is the functional equation of the mining boundary one, ρ1=f1(t).

[0011] d. The same surface settlement values ​​are superimposed for all segment mining sites to obtain the complete surface settlement value.

[0012] e. Based on the relationship between other surface movement and deformation parameters and surface settlement values, obtain data on mining area surface movement and other parameters, including slope, curvature, horizontal movement, and horizontal deformation.

[0013] In the above prediction method, equation (1) includes two sets of Cartesian coordinate systems, an xy coordinate system and a t,ρ coordinate system, where the xy coordinate system represents the surface location of the mining area and the t,ρ coordinate system represents the location of the mining site, and the projections of the two sets of Cartesian coordinate systems onto the horizontal plane overlap.

[0014] In the above prediction method, equation (1) is calculated by a computer. [Effects of the Invention]

[0015] Compared to conventional technologies, the beneficial technical effects of the present invention are as follows:

[0016] 1. Method for predicting surface movement and deformation of any mining working face based on probability integral method (Patent Publication Number CN106446379B) discloses a method for predicting surface movement and deformation of any mining working face based on probability integral method. This method adopts the vector product method to process the mining area for any shaped mining area, and then divides the processed mining area into several approximately rectangular mining areas along the direction for calculation and integration to obtain the surface movement and deformation value of the whole mining area. A method is proposed, and the prediction of the surface movement and deformation value of any shaped mining area, the drawing of two-dimensional, three-dimensional, and contour maps, and the output of calculation data files are realized by using Python language programming. Compared with this, the present invention corresponds to proposing a method for calculating the surface movement and deformation value of any mining area with different concepts. This method can realize the prediction of the surface movement and deformation of complex boundary shape working face mining by combining formula (1) with a computer. The present invention does not need to simplify the mining trace boundary, always maintains the initial shape, and simplifies the calculation steps.

[0017] 2. The present invention can realize the prediction of surface movement and deformation of complex shaped working face mining by the proposed prediction formula for surface subsidence of any boundary shape working face mining and a relatively simple prediction division method, without the need to subdivide the mining trace into a large number of rectangles or triangles, and simplifies the calculation flow without reducing the prediction accuracy.

Brief Description of Drawings

[0018] [Figure 1] It is a schematic diagram of the usage method of formula (1) of the present invention. [Figure 2] It is a schematic diagram of the prediction of surface movement and deformation in the mining area of the present invention.

Modes for Carrying Out the Invention

[0019] Hereinafter, the specific embodiments of the present application will be described in further detail in conjunction with the drawings and examples. It should be pointed out that all the following detailed descriptions are exemplary and are for the purpose of further explaining the present invention. Unless otherwise specified, the technical terms and scientific terms used in this specification have the same meaning as commonly understood by those skilled in the art.

[0020] A method for predicting surface movement and deformation in a mining area based on the segmentation method, comprising the following steps.

[0021] Step 1. Identify the position information of the mining trace boundary in the mining area, obtain the coordinate information of the mining trace boundary, and establish a function formula in the rectangular coordinate system of the horizontal plane of the mining trace boundary based on the coordinate information.

[0022] Step 2. Segment the mining traces based on the criterion of whether the boundaries of the mining traces can be represented by the same function formula.

[0023] Step 3. Calculate for each segment after segmentation using formula (1) to obtain the surface subsidence value generated by this segment of the mining trace. The calculation of formula (1) can be realized by those skilled in the art using a computer.

[0024] JPEG2026076913000003.jpg24170

[0025] In formula (1), W(x,y) is the subsidence value of the surface (x,y) in the mining area, and W cm is the maximum subsidence value of the surface in the mining area, which is a known parameter, H is the burial depth corresponding vertically to the point where the surface in the mining area is located in the coal seam (or mining trace), tanβ is the tangent of the main influence angle, t i-1 is the t-direction coordinate of point P on the mining trace boundary i-1 is the t-direction coordinate of point P on the mining trace boundary i is the t-direction coordinate of point P on the mining trace boundary i is the t-direction coordinate of point P on the mining trace boundary, x is the t-direction coordinate of the surface point in the mining area, y is the ρ-direction coordinate of the surface point in the mining area, the ρ-direction and the t-direction constitute a rectangular coordinate system, π is the circumference ratio, e is the base of the natural logarithm, σ is the integration variable, erf(p) is the probability integral function, p represents a variable, f2(t) is the function formula of the second mining trace boundary, ρ2 = f2(t), and f1(t) is the function formula of the first mining trace boundary, ρ1 = f1(t).

[0026] Step 4. Superimpose the subsidence values of the same surface position for all segment mining traces to obtain the complete surface subsidence value.

[0027] Step 5. Based on the relationship between other surface displacement deformation parameters and surface settlement values, obtain other parameters for mining area surface displacement deformation, including surface inclination, curvature, horizontal movement, and horizontal deformation.

[0028] Equation (1) includes two sets of Cartesian coordinate systems, an xy coordinate system and a t,ρ coordinate system, where the origin of the x and y coordinate system is O, and the origin of the t and ρ coordinate system is O'. The xy coordinate system represents the surface location of the mining area, and the t,ρ coordinate system represents the location of the mining site, and the projections of the two sets of coordinate systems onto the horizontal plane overlap.

[0029] The integral in equation (1) above is with respect to the mining trace and is perpendicular to the integration direction. Therefore, there is no need to calculate boundaries that cannot be expressed by the integration variable, and the connecting boundary from the second boundary point P2 to the third boundary point P3 in Figure 1 does not need to be calculated.

[0030] Example 1: As shown in Figure 2, this mining site has complex boundary shapes such as zigzag and curved boundaries. In Figure 2, the shape of this mining site is formed by the sequential connection of the first boundary point P1, the second boundary point P2, the third boundary point P3, the fourth boundary point P4, the fifth boundary point P5, the sixth boundary point P6, the seventh boundary point P7, and the eighth boundary point P8. A functional equation of the mining site boundary is established based on the t direction, and each of the f 1-2 (t), f 2-3 (t), f 3-3’ (t), f 3’-4 (t), f 4-5 (t), f 5-6 (t), f 6-7 (t), f 7-8 (t), f 8-1 (t) Here, f 1-2 (t) is a line segment perpendicular to the t-direction and cannot be expressed by the f(t) function. Therefore, this boundary is not used as the computational boundary. The line connecting the third boundary point P3 to the fourth boundary point P4 is originally a continuously differentiable curve, and depending on the chosen t-direction, it is necessary to divide this curve into two segments with P3 as the boundary.

[0031] The mining traces are divided based on whether the boundaries of the mining traces can be expressed by the same functional equation. Then, calculations are performed using equation (1) for each region from smallest to largest in the t-direction. Furthermore, if there are multiple integration regions in the same segment (circles 6 and 7 in Figure 2), calculations are performed one by one using equation (1).

[0032] Finally, by superimposing the settlement value at the same surface location obtained from all integration regions, we can obtain the complete surface settlement value. Furthermore, based on the relationship between the settlement value and other surface movement deformation parameters, we can obtain other parameter data for this mining area's surface movement deformation, including surface inclination, curvature, horizontal movement, and horizontal deformation.

[0033] The foregoing describes only preferred embodiments of the present invention, and it should be noted that those skilled in the art can make several further improvements and substitutions without departing from the technical principles of the present invention, and these improvements and substitutions should also be considered to fall within the scope of protection of the present invention. [Explanation of Symbols]

[0034] P1 First boundary point P2 Second boundary point P3 Third boundary point P4 Fourth boundary point P5 Fifth boundary point P6 Sixth boundary point P7 Seventh boundary point P8 Eighth boundary point O x and the origin of the y coordinate system O't and the origin of the ρ coordinate system

Claims

1. A method for predicting surface displacement and deformation of mining areas based on the division method, Step a involves identifying the location information of the mining site boundary in the mining area, obtaining the coordinate information of the mining site boundary, and establishing a functional equation in the horizontal rectangular coordinate system of the mining site boundary based on the coordinate information. Step b involves dividing the mining area based on whether the boundaries of the mining area can be expressed by the same functional equation, Step c involves calculating the ground settlement value that occurs at the mining site of each segment using formula (1) after division, In equation (1), W(x,y) is the settlement value at the mining area surface (x,y), W ｃｍ —The maximum settlement value of the mining area surface is a known parameter. H - The point on the surface of the mining area corresponds perpendicularly to the deposit depth of the coal seam or mining site. tanβ is the main influencing angle tangent. t ｉ－１ -Mining site boundary P ｉ－１ This is the t-coordinate of the point, t ｉ -Mining site boundary P ｉ This is the t-coordinate of the point, x is the t-coordinate of the mining area surface point, y is the coordinate of the mining area surface point in the ρ direction, and the ρ direction and t direction constitute a plane rectangular coordinate system. π is the ratio of a circle's circumference to its diameter, e is the base of the natural logarithm, σ is the integral variable, erf(p) is a stochastic integral function where p represents a variable. f ２ (t) - Functional equation of the mining site boundary 2, ρ ２ = f ２ (t) and f １ (t) - A functional equation of the boundary of the mining trace, ρ １ = f １ (t) in step c, and Step d involves superimposing the same surface settlement values ​​for all segment mining sites to obtain the complete surface settlement value. A method for predicting mining area surface movement deformation based on a division method, characterized by including step e, which involves obtaining other parameters of mining area surface movement deformation, including slope, curvature, horizontal movement, and horizontal deformation, based on the relationship between other ground surface movement deformation parameters and ground surface settlement values.

2. The method for predicting surface movement deformation of a mining area based on the division method described in 1, characterized in that equation (1) includes two sets of Cartesian coordinate systems, an x, y coordinate system and a t, ρ coordinate system, the x, y coordinate system representing the surface location of the mining area and the t, ρ coordinate system representing the location of the mining site, and the projections of the two sets of Cartesian coordinate systems onto the horizontal plane overlap.

3. A method for predicting surface displacement deformation of a mining area based on the division method according to claim 1, characterized in that equation (1) is calculated by a computer.

Images