Method for predicting surface movement and deformation of mining area based on splitting method

Abstract

In order to solve the problem that there is no means to face the prediction of surface movement and deformation in a mining area with a complex boundary-shaped working face, and it is difficult to perform accurate calculations, a method for predicting surface movement and deformation in a mining area based on a segmentation method that improves the efficiency and accuracy of surface movement and deformation prediction in a mining area is provided. 【Solution means】The method of the present invention first identifies the position information of the mining trace boundary in the mining area, obtains the coordinate information of the mining trace boundary, establishes a functional formula in the horizontal rectangular coordinate system of the mining trace boundary based on the coordinate information, and then divides the mining trace based on whether the boundary of the mining trace can be represented by the same functional formula. Next, each segment after division is calculated using the formula to obtain the surface subsidence value generated by this segment of the mining trace. The subsidence values at the same surface position for all segments of the mining trace are superimposed to obtain the complete surface subsidence value. Finally, based on the relationship between other surface movement and deformation parameters and the surface subsidence value, other parameters of surface movement and deformation in the mining area are obtained.

Description

Technical Field

[0001] The present invention relates to the technical field of surface movement and deformation prediction, and specifically to a method for predicting surface movement and deformation in a mining area based on the segmentation method.

Background Art

[0002] The mining of coal seams causes a certain degree of movement and deformation on the surface, thereby adversely affecting the safety and stability of the surface morphology and related buildings (structures). How to accurately predict the degree of movement and deformation that occurs on the surface after coal mining is extremely important for preventing mining damage and environmental protection. The currently commonly used method for predicting surface movement and deformation in a mining area is the probability integral method. The principle of this method is to regard the strata covered on the mining trace as a plurality of small balls, and take the probability of the small balls falling on the mining trace as the corresponding surface subsidence value. By this method, a relatively accurate surface movement and deformation situation can be obtained.

[0003] Currently, there are relatively many research reports on predicting surface movement and deformation during the mining of a rectangular working face by applying the probability integral method, and there are errors or relatively large errors in predicting surface movement and deformation during the mining of a working face with a complex shape. One of the main reasons is that when predicting surface movement and deformation by the probability integral method in the conventional technology, in the application process from two-dimensional (strata vertical section) to three-dimensional, researchers often ignore the relationship between two orthogonal integral variables and fixedly consider that there is no relationship between the variables in the two integral directions. This does not affect the prediction results when solving the mining of a rectangular working face, but errors occur when predicting the surface deformation during the mining of a non-rectangular working face, which makes it easy to generate prediction errors when solving the prediction of surface movement and deformation in a mining area with a complex boundary shape working face by this method. The second main reason is that predicting the surface movement and deformation in a mining area with a complex boundary shape working face needs to be carried out according to a certain method, but currently, there is still no efficient and accurate prediction method.

[0004] It can be seen that the improvement of the conventional technology is awaited.

Summary of the Invention

Problems to be Solved by the Invention

[0005] In order to solve the problem that there is no means to face the prediction of surface movement and deformation in a mining area with a complex boundary-shaped working face, and accurate calculation is difficult, the present invention provides a method for predicting surface movement and deformation in a mining area based on a segmentation method to improve the efficiency and accuracy of surface movement and deformation prediction in the mining area.

Means for Solving the Problems

[0006] The technical solution of the present invention is as follows.

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

[0008] a. 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 functional formula in the horizontal rectangular coordinate system of the mining trace boundary based on the coordinate information. b. Divide the mining trace based on whether the boundary of the mining trace can be represented by the same functional formula. c. Calculate each segment after division using formula (1) to obtain the surface subsidence value generated by this segment of the mining trace.

[0009] JPEG0007704480000002.jpg25170

[0010] In formula (1), W(x, y) is the subsidence value of the surface of the mining area (x, y), and W cm is the maximum subsidence value of the surface of the mining area, which is a known parameter. H is the burial depth at which the point where the surface of the mining area is located corresponds vertically to 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. t i-1 is the t-direction coordinate of point P on the mining trace boundary. t i is the t-direction coordinate of point P on the mining trace boundary. t it is the coordinate of the point in the t direction, x is the coordinate of the point on the surface of the mining area in the t direction, y is the coordinate of the point on the surface of the mining area in the ρ direction, the ρ direction and the t direction form a plane rectangular coordinate system, π is the pi, 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).

[0011] d. Superimpose the subsidence values at the same surface position for all segment mining traces to obtain the complete surface subsidence value.

[0012] e. Based on the relationship between other surface movement and deformation parameters and the surface subsidence value, obtain other parameter data of the surface movement in the mining area, including inclination, curvature, horizontal movement and horizontal deformation.

[0013] In the above prediction method, the formula (1) includes two sets of rectangular coordinate systems, namely the x,y coordinate system and the t,ρ coordinate system. The x,y coordinate system represents the surface position of the mining area, and the t,ρ coordinate system represents the mining trace position. The projections of the two sets of rectangular coordinate systems onto the horizontal plane overlap.

[0014] In the above prediction method, calculate the formula (1) by a computer.

Advantages of the Invention

[0015] Compared with the prior art, the beneficial technical effects of the present invention are as follows.

[0016] 1. Prediction method for surface movement and deformation of any mining working face based on probability integral method (Patent Publication No. CN106446379B) discloses a prediction method for 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 along the direction into several approximate rectangular mining areas for calculation and integration to obtain the surface movement and deformation value of the entire mining area. It proposes a method that uses Python language programming to realize the prediction of the surface movement and deformation value of any shaped mining area, as well as the drawing of two-dimensional, three-dimensional, and contour maps and the output of calculation data files. Compared with this, the present invention corresponds to proposing a calculation method for 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 surface subsidence prediction formula for any boundary shape working face mining and a relatively simple prediction division method. It does not need to divide the mining trace into a large number of rectangles or triangles, simplifies the calculation flow without reducing the prediction accuracy.

Brief Description of the Drawings

[0018]

Figure 1

Figure 2

Embodiments for Carrying Out the Invention

[0019] Hereinafter, the specific embodiments of the present application will be described in more detail by combining with the drawings and examples. It should be pointed out that all the following detailed descriptions are exemplary and 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 generally 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 horizontal rectangular coordinate system of the mining trace boundary based on the coordinate information.

[0022] Step 2. Segment the mining trace 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] JPEG0007704480000003.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 at which the point where the surface in the mining area is located corresponds vertically to 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 the surface point in the mining area, 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 form a rectangular coordinate system in the plane, π 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 relationships between other surface movement and deformation parameters and the surface subsidence value, obtain other parameters of surface movement and deformation in the mining area, including surface inclination, curvature, horizontal movement, and horizontal deformation.

[0028] Equation (1) includes two sets of rectangular coordinate systems, namely the x,y coordinate system and the t,ρ coordinate system. The origin of the x,y coordinate system is O, and the origin of the t,ρ coordinate system is O’. The x,y coordinate system represents the surface position of the mining area, and the t,ρ coordinate system represents the position of the mining trace. The projections of the two sets of coordinate systems onto the horizontal plane overlap.

[0029] The integration of the above Equation (1) is with respect to the mining trace. Since it is perpendicular to the integration direction, for the boundaries that cannot be expressed by the integration variable, there is no need to calculate. For example, 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 trace has boundaries with complex shapes such as zigzag boundaries and curved boundaries. In Figure 2, this mining trace shape is formed by connecting 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 in sequence. Based on the t direction, establish the functional expressions of the mining trace boundaries, which are 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) respectively. Here, f 1-2 (t) is a line segment perpendicular to the t direction and cannot be expressed by the f(t) function. This boundary is not regarded as a calculation boundary. The connecting line of the boundary points from the third boundary point P3 to the fourth boundary point P4 is originally a continuously differentiable curve. Since it is necessary to divide this curve into two segments with P3 as the boundary according to the taken t direction, this is the reason.

[0031] Based on whether the boundaries of the mining traces can be expressed by the same function formula, the mining traces are divided, and for each region from the smallest to the largest in the t direction, the above formula (1) is used for calculation. Also, when there are multiple integration regions in the same segment (circles 6 and 7 in Figure 2), the above formula (1) is used for calculation one by one.

[0032] Finally, the subsidence values at the same surface positions as the calculation results obtained in all integration regions are superimposed to obtain the complete surface subsidence value. Furthermore, based on the relationships between other surface movement and deformation parameters and the subsidence value, other parameter data of the surface movement and deformation in this mining area, including surface inclination, curvature, horizontal movement, and horizontal deformation, can be obtained.

[0033] What is described above is only a preferred embodiment of the present invention. It should be pointed out that for those skilled in the art of this technology, several more improvements and substitutions can be made without departing from the technical principle of the present invention, and these improvements and substitutions should also be regarded as falling within the protection scope of the present invention.

Explanation of Reference Signs

[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 Origin of the x and y coordinate systems O’ Origin of the t and ρ coordinate systems

Claims

1. A method for predicting surface movement and deformation in a mining area based on the segmentation method, comprising: Step a: identifying the position information of the mining trace boundary in the mining area, obtaining the coordinate information of the mining trace boundary, and establishing a functional equation in the rectangular coordinate system of the horizontal plane of the mining trace boundary based on the coordinate information; Step b: dividing the mining traces based on the criterion of whether the boundaries of the mining traces can be represented by the same functional equation; Step c: calculating for each segment after division using Equation (1) to obtain the surface subsidence value generated by this segment of the mining trace, where In Equation (1), W(x, y) is the subsidence value of the mining area surface (x, y); W ｃｍ - The maximum subsidence value on the surface of the mining area, which is a known parameter, H is the burial depth corresponding vertically to the coal seam or the mining trace at the point where the mining area surface is located; tanβ is the tangent of the main influence angle; t ｉ－１ - Mining trace boundary P ｉ－１ is the t-direction coordinate of the point, t ｉ - Mining trace boundary P ｉ It is the t-direction coordinate of the point, x is the t-direction coordinate of the mining area surface point; y is the ρ-direction coordinate of the mining area surface point, and the ρ-direction and the t-direction constitute a rectangular coordinate system in the plane; π is the ratio of the circumference of a circle to its diameter; e is the base of the natural logarithm; σ is the integration variable; erf(p) is the probability integral function, where p represents a variable; f ２ (t) - functional equation of the boundary of the mining trace two, ρ ２ = f ２ (t), and f １ (t) - A functional equation of the boundary of the mining trace, ρ １ = f １ (t) in step c, and Step d: superimposing the subsidence values of the same surface position for all segment mining traces to obtain the complete surface subsidence value; Step e: obtaining other parameters of the surface movement and deformation in the mining area, including inclination, curvature, horizontal movement, and horizontal deformation, based on the relationship between other surface movement and deformation parameters and the surface subsidence value. The method for predicting surface movement and deformation in a mining area based on the segmentation method is characterized in that Equation (1) includes two sets of rectangular coordinate systems, namely the x, y coordinate system and the t, ρ coordinate system respectively. The x, y coordinate system represents the mining area surface position, the t, ρ coordinate system represents the mining trace position, and the projections of the two sets of rectangular coordinate systems onto the horizontal plane overlap.

2. The method for predicting surface movement and deformation in a mining area based on the segmentation method according to Claim 1, characterized in that Equation (1) is calculated by a computer.

Images