Prediction method and device for residual deformation of surface goaf in mountainous areas

Abstract

The present application relates to the field of residual subsidence monitoring technology, and in particular to a method and device for predicting residual deformation of a goaf area on the surface of a mountainous area, wherein the method comprises: dividing different mountainous terrains and landforms to obtain landform data and actual landform data on the surface of a target mountainous area; obtaining an initial residual deformation value of the surface of a target mountainous area generated by an initial mountainous area surface residual deformation prediction model based on the residual deformation value and landform data caused by flat land mining; correcting the initial residual deformation value to construct a mountainous area surface residual deformation prediction model of the surface of the target mountainous area, and inputting the actual landform data into the pre-constructed mountainous area surface residual deformation prediction model to obtain the residual deformation value of the surface of the target mountainous area. Thus, the problem that the mathematical models of residual subsidence in existing mining areas are mostly established on plain mining areas and cannot meet the requirements for predicting residual deformation of the surface of mountainous areas under complex terrain conditions is solved.

Description

Technical Field

[0001] The present application relates to the technical field of residual subsidence monitoring, and in particular to a method and device for predicting residual deformation of a surface goaf area in a mountainous area. Background Art

[0002] Under the influence of factors such as complex topography and landforms, the surface movement and deformation caused by underground mining of coal resources in mountainous terrain has a more serious impact on surface buildings than in plain areas. At the same time, in terms of the problem of residual deformation of the surface in old mining areas, mountainous areas with relatively complex topography and landforms will also have a certain degree of influence on the residual deformation.

[0003] In the related technology, the equivalent variable mining thickness prediction model can be used to monitor the residual deformation in the plain area, and the traditional methods such as leveling and wire measurement can be used to monitor the residual subsidence of mining in mountainous areas; GNSS (Global Navigation Satellite System) can also be used to realize online monitoring of the residual subsidence of mining in mountainous areas / mining areas.

[0004] However, in the related technologies, the measurement is difficult and the measurement accuracy is low. It is easily disturbed by external factors. In addition, the existing mathematical models of residual subsidence in mining areas are mostly established in plain mining areas. The limited mathematical models hinder the prediction of residual deformation of the surface in mountainous areas and cannot meet the requirements for predicting residual deformation of the surface in mountainous areas under complex terrain conditions. Improvement is urgently needed. Summary of the invention

[0005] The present application provides a method and device for predicting residual deformation of surface goaf in mountainous areas, in order to solve the problems in related technologies, such as high measurement difficulty, low measurement accuracy, and easy interference by external factors. In addition, the existing mathematical models of residual subsidence in mining areas are mostly established in plain mining areas. The limited mathematical models hinder the prediction of residual deformation of the surface in mountainous areas and cannot meet the requirements for predicting residual deformation of the surface in mountainous areas under complex terrain conditions.

[0006] The first aspect of the present application provides a method for predicting residual deformation of a goaf on the surface of a mountainous area, which is applied to a model building stage, wherein the method comprises the following steps: dividing different mountainous terrains to obtain division results, and obtaining at least one landform data in the target mountainous area surface based on the division results; obtaining the initial residual deformation value of the target mountainous area surface by using an initial mountainous area surface residual deformation prediction model through the residual deformation value caused by flat ground mining and the at least one target landform data; correcting the initial residual deformation value, and constructing a mountainous area surface residual deformation prediction model of the target mountainous area surface based on the corrected initial residual deformation value and the initial mountainous area surface residual deformation prediction model, so as to use the mountainous area surface residual deformation prediction model to predict the residual deformation of the goaf on the surface of the target mountainous area, wherein the mountainous area surface residual deformation prediction model includes at least one of a residual subsidence calculation formula, a concave landform residual horizontal movement calculation formula, and a convex landform residual horizontal movement calculation formula.

[0007] Optionally, in one embodiment of the present application, the residual deformation value caused by flat land mining and the at least one target landform data are used to obtain the initial residual deformation value of the target mountain surface using an initial mountain surface residual deformation prediction model, including: obtaining the residual deformation value caused by flat land mining; based on the residual deformation value, inputting the at least one target landform data into the initial mountain surface residual deformation prediction model to obtain the initial residual deformation value of the target mountain surface.

[0008] Optionally, in one embodiment of the present application, the correcting the initial residual deformation value includes: judging whether there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area; if there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area, correcting the initial residual deformation value at the positive and negative slope change points and / or concave-convex landform change points to obtain a corrected initial residual deformation value.

[0009] Optionally, in one embodiment of the present application, the correcting the initial residual deformation value includes: determining whether there is a convex landform on the surface of the target mountain area; if there is a convex landform on the surface of the target mountain area, correcting the initial residual horizontal movement value in the convex landform to obtain a corrected initial residual deformation value.

[0010] Optionally, in one embodiment of the present application, the residual sinking calculation formula may be, but is not limited to:

[0011]

[0012] Among them, W′ 0(x, t) is the residual subsidence value at the surface x in the mountain area at time t, h n is the block height of the left hole area, n 0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, n0′ is the number of small blocks in the right hollow area, h n′ The height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the smoothed surface point, m 1 is the residual activatable thickness in the middle area;

[0013] The expression of the calculation formula for the residual horizontal movement of the concave landform can be, but is not limited to,:

[0014]

[0015] Among them, I b1 ,I z and I b2 are the residual tilts of the left cavity area, the middle cavity area, and the right cavity area, respectively, and U′ 0 (x, t) is the residual horizontal movement value at the surface x in the mountain area at time t;

[0016] The expression of the calculation formula for the residual horizontal movement of the convex landform can be, but is not limited to,:

[0017]

[0018] The second aspect of the present application provides a method for predicting the residual deformation of a mined-out area on the surface of a mountainous area, which is applied to the model usage stage, wherein the method includes the following steps: dividing different mountainous terrains to obtain division results, and obtaining at least one actual landform data on the surface of the target mountainous area based on the division results; inputting the at least one actual landform data into a pre-constructed mountainous surface residual deformation prediction model to obtain the residual deformation value of the target mountainous surface, wherein the pre-constructed mountainous surface residual deformation prediction model is obtained from at least one actual landform data on the surface of the target mountainous area.

[0019] In a third aspect, an embodiment of the present application provides a device for predicting residual deformation of a goaf area on a mountain surface, which is applied to a model building stage, wherein the device comprises: a first acquisition module, which is used to divide different mountain terrains and geomorphologies to obtain a division result, and obtain at least one geomorphological data on the surface of a target mountain area based on the division result to obtain at least one target geomorphological data on the surface of a target mountain area; a first generation module, which is used to obtain an initial residual deformation value of the surface of a target mountain area by using an initial mountain surface residual deformation prediction model through a residual deformation value caused by flat ground mining and the at least one target geomorphological data. The deformation value includes at least one of an initial residual subsidence value and an initial residual horizontal movement value; a construction module is used to correct the initial residual deformation value, and construct a mountain surface residual deformation prediction model of the target mountain surface according to the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model, so as to use the mountain surface residual deformation prediction model to predict the residual deformation of the goaf area of ​​the target mountain surface, wherein the mountain surface residual deformation prediction model includes at least one of a residual subsidence calculation formula, a concave landform residual horizontal movement calculation formula and a convex landform residual horizontal movement calculation formula.

[0020] Optionally, in one embodiment of the present application, the first generation module includes: an acquisition unit for acquiring the residual deformation value caused by flat land mining; a generation unit for inputting the at least one target landform data into the initial mountain surface residual deformation prediction model based on the residual deformation value to obtain the initial residual deformation value of the target mountain surface.

[0021] Optionally, in one embodiment of the present application, the construction module includes: a first judgment unit, used to judge whether there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area; a first correction unit, used to correct the initial residual deformation value at the positive and negative slope change points and / or concave-convex landform change points when there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area, so as to obtain a corrected initial residual deformation value.

[0022] Optionally, in one embodiment of the present application, the construction module includes: a second judgment unit, used to judge whether there is a convex landform on the surface of the target mountain area; a second correction unit, used to correct the initial residual horizontal movement value in the convex landform when there is a convex landform on the surface of the target mountain area, so as to obtain a corrected initial residual deformation value.

[0023] Optionally, in one embodiment of the present application, the residual sinking calculation formula may be, but is not limited to:

[0024]

[0025] Among them, W′ 0 (x, t) is the residual subsidence value at the surface x in the mountain area at time t, h n is the block height of the left hole area, n 0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, n0′ is the number of small blocks in the right hollow area, h n′ The height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the smoothed surface point, m 1 is the residual activatable thickness in the middle area;

[0026] The expression of the calculation formula for the residual horizontal movement of the concave landform can be, but is not limited to,:

[0027]

[0028] Among them, I b1 ,I z and I b2 are the residual tilts of the left cavity area, the middle cavity area, and the right cavity area, respectively, and U′ 0 (x, t) is the residual horizontal movement value at time t at point x on the mountain surface.

[0029] The expression of the calculation formula for the residual horizontal movement of the convex landform can be, but is not limited to,:

[0030]

[0031] The fourth aspect of the present application provides a device for predicting the residual deformation of a mined-out area on the surface of a mountainous area, which is applied to the model use stage, wherein the device includes: a second acquisition module, which is used to divide the terrain and landforms of different mountainous areas, obtain a division result, and obtain at least one actual landform data in the target mountainous area surface based on the division result; a second generation module, which is used to input the at least one actual landform data into a pre-constructed mountainous area surface residual deformation prediction model to obtain the residual deformation value of the target mountainous area surface, wherein the pre-constructed mountainous area surface residual deformation prediction model is obtained from at least one actual landform data of the target mountainous area surface.

[0032] The fifth aspect of the present application provides an electronic device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for predicting residual deformation of surface goaf areas in mountainous areas as described in the above embodiments.

[0033] The sixth aspect of the present application provides a computer-readable storage medium, which stores a computer program. When the program is executed by a processor, it implements the above-mentioned method for predicting residual deformation of surface goaf areas in mountainous areas.

[0034] The seventh aspect of the present application provides a computer program product, including a computer program, which, when executed, implements the above-mentioned method for predicting residual deformation of surface goaf areas in mountainous areas.

[0035] The embodiment of the present application can obtain the actual landform data of the target mountain surface by dividing different mountain landforms, and construct the mountain surface residual deformation prediction model of the target mountain surface through the residual deformation value caused by flat land mining and the initial mountain surface residual deformation prediction model, and then obtain the residual deformation value of the target mountain surface, and solve the problem of residual deformation of the mountain surface by establishing a new method for obtaining the inclination angle of the smoothed terrain trend surface and the mountain surface characteristic coefficient, and provide a calculation model for the residual deformation of the mountain surface goaf. Thus, the problems in the related technology that the measurement is difficult, the measurement accuracy is low, it is easy to be disturbed by external factors, and the mathematical models of the residual subsidence of the existing mining area are mostly established in the plain mining area, and the limited mathematical models hinder the prediction of the residual deformation of the mountain surface, and cannot meet the requirements of the prediction of the residual deformation of the mountain surface under complex terrain conditions.

[0036] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0038] Figure 1 A flowchart of a method for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application;

[0039] Figure 2 A schematic block diagram of a multi-faceted slope provided according to an embodiment of the present application;

[0040] Figure 3 A schematic block diagram of an equivalent variable-thickness sampling and variable-thickness sampling probability integral method model provided according to an embodiment of the present application;

[0041] Figure 4 A schematic block diagram of a cavity area division according to an embodiment of the present application;

[0042] Figure 5A flowchart of the working principle of a method for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application;

[0043] Figure 6 A schematic block diagram of a device for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application;

[0044] Figure 7 A flowchart of another method for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application;

[0045] Figure 8 A schematic diagram of the residual deformation value of the surface of a target mountain area provided according to an embodiment of the present application;

[0046] Fig. 9 A schematic diagram showing a comparison of residual deformation values ​​of a target mountain surface and a flat surface according to an embodiment of the present application;

[0047] Fig.10 A schematic block diagram of another device for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application;

[0048] Fig.11 It is a schematic diagram of the structure of an electronic device provided according to an embodiment of the present application. DETAILED DESCRIPTION

[0049] Embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0050] The following describes the prediction method and device of the residual deformation of the surface goaf of the mountainous area of ​​the embodiment of the present application with reference to the accompanying drawings. In view of the fact that the measurement difficulty and measurement accuracy mentioned in the above background technology are relatively large, and the measurement is very easy to be disturbed by external factors, and the mathematical models of residual subsidence of existing mining areas are mostly established on plain mining areas, the limited mathematical models hinder the prediction of the residual deformation of the mountainous surface, and cannot meet the requirements of the prediction of the residual deformation of the mountainous surface under complex terrain conditions, the present application provides a prediction method of the residual deformation of the surface goaf of the mountainous area, in which the actual landform data of the target mountainous surface can be obtained by dividing different mountainous terrains, and the residual deformation value caused by flat mining and the initial mountainous surface residual deformation prediction model are used to construct the mountainous surface residual deformation prediction model of the target mountainous surface, and then the residual deformation value of the target mountainous surface is obtained, and a new method for obtaining the inclination angle of the smoothed terrain trend surface and the mountainous surface characteristic coefficient is established to solve the problem of residual deformation of the mountainous surface, and a calculation model is provided for the residual deformation of the surface goaf of the mountainous area. As a result, the problems in related technologies such as high measurement difficulty, low measurement accuracy, and susceptibility to interference from external factors are solved, and the existing mathematical models of residual subsidence in mining areas are mostly established in plain mining areas. The limited mathematical models hinder the prediction of residual deformation of the surface in mountainous areas and cannot meet the requirements for predicting residual deformation of the surface in mountainous areas under complex terrain conditions.

[0051] Specifically, Figure 1 The present invention is a flowchart of a method for predicting residual deformation of a surface goaf area in a mountainous area according to an embodiment of the present application.

[0052] like Figure 1 As shown, the prediction method of residual deformation of the surface goaf area in the mountainous area is applied in the model building stage, wherein the method includes the following steps:

[0053] In step S101, different mountainous terrains are divided to obtain a division result, and at least one landform data on the surface of the target mountainous area is obtained based on the division result.

[0054] It can be understood that the embodiments of the present application can divide different mountain terrains and landforms, such as being divided into positive slopes, reverse slopes, and multi-faceted slopes according to different surface inclinations and changes in inclination angles, and the present application does not make specific restrictions. The multi-faceted slope is composed of a certain number of positive and reverse slopes, and the certain number can be set by technicians in this field according to actual conditions, and the present application does not make specific restrictions. It can also be divided into convex landforms and concave landforms according to the morphological characteristics of the terrain and geomorphological principles, and the present application does not make specific restrictions. Convex landforms can be understood as the terrain surface presenting a convex shape relative to the surrounding area, and concave landforms can be understood as the terrain surface presenting a concave shape relative to the surrounding area, and the present application does not make specific restrictions.

[0055] As a possible implementation method, the embodiment of the present application can divide different mountain terrains and landforms, and then obtain at least one landform data in the target mountain surface. Among them, the at least one landform data in the embodiment of the present application can include but is not limited to the number, coordinate area, length, trend surface inclination, etc. of positive slopes and / or reverse slopes, the number, coordinate area, etc. of convex landforms and / or concave landforms in the target image, which can be set by people in this field according to actual conditions, and this application does not make specific restrictions.

[0056] For example, the embodiment of the present application selects Figure 2 A typical multi-faceted slope topography is shown in Figure 2 In the embodiment of the present application, at least one geomorphic data obtained can be, but is not limited to, the left side of the mined coal seam as the coordinate origin, the two sides of the figure are unmined coal pillar areas, the coordinate areas are -300m~0m and 1200m-1500m, and the middle is a limited mining area, the coordinate area is 0m-1200m. Among them, the multi-faceted slope in the embodiment of the present application includes two positive slopes and one reverse slope, a convex landform and a concave landform. Further, the two positive slopes and one reverse slope are 600m long, the positive slope coordinate area is -300m-300m and 900m-1500m, the reverse slope coordinate area is 300m-900m, the convex landform coordinate area is 1400m-1500m, and the concave landform coordinate area is 200m-400m. Furthermore, according to the embodiment of the present application, after measurement and calculation, the inclination angle of the trend surface of the leftmost positive slope is 15°, the inclination angle of the trend surface of the rightmost positive slope is 20°, and the inclination angle of the trend surface of the middle reverse slope is -25°.

[0057] In step S102, the initial residual deformation value of the target mountain surface is obtained by using the residual deformation value caused by flat land mining and at least one target landform data and using the initial mountain surface residual deformation prediction model.

[0058] As a possible implementation method, the embodiment of the present application can obtain the initial residual deformation value of the target mountain surface through the residual deformation value caused by flat land mining and the initial mountain surface residual deformation prediction model based on at least one target landform data.

[0059] In addition, in the embodiments of the present application, the initial residual deformation value may include, but is not limited to, an initial residual subsidence value and an initial residual horizontal movement value, etc., and the present application does not impose any specific limitation.

[0060] Optionally, in one embodiment of the present application, the initial residual deformation value of the target mountain surface is obtained by using the residual deformation value caused by flat land mining and at least one target landform data, using the initial mountain surface residual deformation prediction model, including: obtaining the residual deformation value caused by flat land mining; based on the residual deformation value, inputting at least one target landform data into the initial mountain surface residual deformation prediction model to obtain the initial residual deformation value of the target mountain surface.

[0061] In some embodiments, the present application embodiment can obtain the residual deformation value caused by flat ground mining through the probability integral method prediction model of variable mining thickness, and the specific content is as follows:

[0062] The present application embodiment is based on the analysis of relevant domestic and foreign data, and shows that the residual surface movement deformation value after the working face stops mining is mainly divided into two parts: the hollow area at the boundary of the goaf area and the compaction area in the middle of the goaf area. Figure 3 As shown. Further, the probability integral method is currently recognized as the most commonly used method for calculating surface subsidence. This method considers surface subsidence to be a random event and discretizes the entire mining area into many tiny units. Surface subsidence is the sum of the effects of all mining units on the surface, and the integral method can be used to directly solve it. Therefore, the subsidence value of point x caused by unit mining in the embodiment of the present application can be, but is not limited to:

[0063]

[0064] Among them, r is the main influence radius.

[0065] The related technology divides the goaf into zones and proposes a probability integral method prediction model of equivalent variable mining thickness by using the superposition idea. Usually, the cavity area is simplified to the equivalent mining height m′. Then, according to the calculation formula of any point on the ground by the probability integral method and considering the symmetry of the goaf, the probability integral method model of equivalent variable mining thickness can be obtained, which can be but not limited to:

[0066] W 0 (x) = ∫ AC m′W e (xs)ds+∫ CD m 1 W e (xs)ds+∫ DB m′W e (xs)ds, (2)

[0067]

[0068] U 0 (x) = brI 0 (x), (4)

[0069] Among them, W0 (x) is the residual settlement value of a point on the surface, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, m′ and m 1 Respectively represent the residual activation thickness of the boundary area and the middle area, W e (x) is the unit subsidence caused by unit mining, m is the mining thickness, I 0 (x) is the residual tilt, I b1 ,I z and I b2 They are the residual tilts of the left cavity area, the middle cavity area, and the right cavity area, respectively. 0 (x) is the residual horizontal shift value.

[0070] The reason why the related art simplifies the void area into an equivalent mining height m′ is that the void area thickness in the probability integral method prediction model of equivalent variable mining thickness changes linearly. If the void area thickness is introduced into a linear function with a variable x, it will be impossible to integrate when calculating in the probability integral method model. However, it can still be based on the most basic integration idea. Further, the embodiment of the present application can improve the accuracy by dividing the void area into multiple small areas for calculation, such as Figure 4 As shown, in Figure 4 In the figure, the lower part is the small blocks divided in the boundary hole area, and the width of each block is the same, and the height is the same as the original height corresponding to the center line of the block. Among them, the embodiment of the application can assume that the number of small blocks is n 0 , from which the height of each block can be calculated, where h n is the block height of the left hole area, h n′ is the block height of the right hole area, and its expressions can be:

[0071]

[0072] m 1 =m(1-q), (6)

[0073]

[0074] Among them, n is the block number of the left hole area, from left to right, 1, 2, 3..., m, m 1 are the residual reactivatable thickness of the thick mining area and the middle area, respectively; n′ is the block number of the right hollow area; and n0′ is the number of small blocks in the right hollow area.

[0075] Furthermore, the embodiment itself can refine the cavity area through equations (5) to (7), and use the idea of ​​superposition to construct a probability integral method prediction model of variable mining thickness, and then obtain the residual deformation value caused by flat ground mining. The expression of the probability integral method prediction model of variable mining thickness constructed can be but is not limited to:

[0076]

[0077] In some embodiments, the embodiments of the present application may input at least one target landform data into an initial mountain surface residual deformation prediction model based on the residual deformation value to obtain the initial residual deformation value of the target mountain surface. The specific contents are as follows:

[0078] It can be understood by those skilled in the art that the surface movement of the target mountain area in the embodiment of the present application includes the surface movement caused by coal mining and the surface movement caused by hillside sliding. Therefore, the residual deformation value generated by the surface movement deformation of the target mountain area can be estimated by superimposing the residual deformation value caused by flat land mining and the mountain sliding amount, which can be expressed as but not limited to:

[0079]

[0080]

[0081]

[0082]

[0083] Among them, W(x), U(x), I(x), K(x), and E(x) are the values ​​of flat ground subsidence, horizontal movement, inclination, curvature, and horizontal deformation under the same geological mining conditions, which can be obtained by the flat ground probability integral prediction formula. ΔW(x), ΔU(x), ΔI(x), ΔK(x), and ΔE(x) are the increments of subsidence, horizontal movement, inclination, curvature, and horizontal deformation caused by mining slippage on the slope affected by mining.

[0084] The additional subsidence and additional horizontal movement values ​​caused by mining-induced sliding on the mountain surface can be expressed as, but not limited to:

[0085] ΔW(x)=D x P(x)W(x)tg 2 α x , (13)

[0086] ΔU(x)=D x P(x)W(x)tgα x , (14)

[0087] Among them, ΔW(x) and ΔU(x) are the additional subsidence and additional horizontal movement caused by mining slip, D x is the surface characteristic coefficient of the mountainous area, the value is referenced in Table 1, P(x) is the slip influence function, W(x) is the flat ground subsidence value under the same geological mining conditions, α xis the inclination angle of the terrain trend surface of the smoothed surface point. Wherein, Table 1 is a schematic table of mountain surface characteristic coefficients provided according to an embodiment of the present application.

[0088] Table 1

[0089]

[0090] If there is sufficient actual observation data, the slip influence function P(x) can be expressed as, but not limited to:

[0091]

[0092] or

[0093]

[0094] If there is no sufficient actual observation data, the slip influence function P(x) can be expressed as, but not limited to:

[0095]

[0096] Among them, A, P, and t are slip influence parameters.

[0097] At this time, the calculation formula of the prediction model of surface mining subsidence in mountainous areas can be but is not limited to:

[0098] W′(x)=W(x)+D x P(x)W(x)tan 2 α x , (18)

[0099] U′(x)=U(x)+D x P(x)W(x)tanα x , (19)

[0100] Furthermore, the calculation formula of the mountain surface residual deformation prediction model can be but is not limited to:

[0101] W′ 0 (x) = W 0 (x)+D x P(x)W 0 (x)tan 2 α x , (20)

[0102] U′ 0 (x) = U 0 (x)+D x P(x)W 0 (x)tanα x , (twenty one)

[0103] Furthermore, the embodiment of the present application can substitute equation (8), equation (3) and equation (3) into equation (20) and equation (21), and introduce Knothe time function to characterize the relationship between residual deformation and time, and then the embodiment of the present application can obtain the initial mountain surface residual deformation prediction model to obtain the initial residual deformation value of the target mountain surface. Among them, the calculation formula of the initial mountain surface residual deformation prediction model can be but not limited to:

[0104]

[0105]

[0106] Among them, W′ 0 (x, t) is the residual subsidence value at the surface x in the mountain area at time t, U′ 0 (x, t) is the residual horizontal movement value at time t at point x on the mountain surface.

[0107] In step S103, the initial residual deformation value is corrected, and a mountain surface residual deformation prediction model of the target mountain surface is constructed based on the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model, so as to use the mountain surface residual deformation prediction model to predict the residual deformation of the goaf area of ​​the target mountain surface, wherein the mountain surface residual deformation prediction model includes at least one of the residual subsidence calculation formula, the concave landform residual horizontal movement calculation formula and the convex landform residual horizontal movement calculation formula.

[0108] As a possible implementation method, the embodiment of the present application can correct the initial residual deformation value obtained by the initial mountain surface residual deformation prediction model, and then construct a mountain surface residual deformation prediction model of the target mountain surface based on the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model.

[0109] Among them, in the embodiments of the present application, the mountain surface residual deformation prediction model may include but is not limited to the residual subsidence calculation formula, the concave landform residual horizontal movement calculation formula and the convex landform residual horizontal movement calculation formula, etc., and the present application does not make specific restrictions.

[0110] Optionally, in one embodiment of the present application, the initial residual deformation value is corrected, including: judging whether there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area; if there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area, then correcting the initial residual deformation value at the positive and negative slope change points and / or concave-convex landform change points to obtain a corrected initial residual deformation value.

[0111] It can be understood by those skilled in the art that there are some problems in the present application embodiment after using the initial mountain surface residual deformation prediction model for calculation. When calculating the residual subsidence and residual horizontal movement at the change point of the concave-convex landform and the change point of the trend surface inclination angle, the terrain trend surface inclination angle α x and mountain surface characteristic coefficient D x The change here is large, and the residual subsidence and residual horizontal movement at the change point of the concavo-convex landform and the change point of the trend surface inclination will suddenly change. At this time, it appears as a fault change on the image, which is not in line with reality. Therefore, the embodiment of the present application can correct the initial residual deformation value. The specific process can be:

[0112] The embodiment of the present application can solve the problems existing in the above-mentioned initial mountain surface residual deformation prediction model and smooth the terrain trend surface inclination angle α x and mountain surface characteristic coefficient D x The value of α is created, that is, a positive slope or reverse slope with a slope change of no more than 5°, smoothing a terrain trend surface inclination angle α x , and because of the mutation problem at the positive and negative slope change points and near the concave and convex landform change points, the image is faulted. At this time, the embodiment of the present application can use the nonlinear least squares method or polynomial to fit a smooth arc curve such as a hyperbola, an ellipse or a parabola to solve the above mutation problem; the terrain trend surface inclination angle α can also be calculated at the positive and negative slope change points. x The slope of one trend surface changes evenly to another trend surface. At the point where the concave-convex landform changes, the surface characteristic coefficient D of the mountain area is x The surface characteristic coefficient of one mountain area is uniformly transitioned to the surface characteristic coefficient of another mountain area, thereby solving the above-mentioned mutation problem.

[0113] Optionally, in one embodiment of the present application, the initial residual deformation value is corrected, including: determining whether there is a convex landform on the surface of the target mountain area; if there is a convex landform on the surface of the target mountain area, correcting the initial residual horizontal movement value in the convex landform to obtain a corrected initial residual deformation value.

[0114] It can be understood by those skilled in the art that there are still some problems found after calculation using the initial mountain surface residual deformation prediction model in the embodiment of the present application. The residual horizontal movement in the convex landform also has the problem of not being consistent with the actual law. On the positive slope of the convex landform, the slip residual horizontal movement is calculated to be a negative value, and on the reverse slope of the convex landform, the slip residual horizontal movement is a positive value, which is contrary to the actual situation. At this time, the embodiment of the present application can correct the initial residual horizontal movement value in the convex landform to obtain a corrected initial residual deformation value, and the specific content can be:

[0115] The embodiments of the present application take into account that horizontal movement is mainly related to positive and negative slopes. Therefore, when calculating residual horizontal movement at convex landforms, the calculation formula for residual horizontal movement of the mountain surface can be modified to obtain a law that conforms to actual horizontal movement.

[0116] The modified residual horizontal movement calculation formula can be, but is not limited to, expressed as:

[0117]

[0118] Optionally, in one embodiment of the present application, the residual sinking calculation formula may be, but is not limited to:

[0119]

[0120] Among them, W′ 0 (x, t) is the residual subsidence value at the surface x in the mountain area at time t, h n is the block height of the left hole area, n 0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, n0′ is the number of small blocks in the right hollow area, h n , the height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the smoothed surface point, m 1 It is the residual activatable thickness in the middle area.

[0121] The calculation formula of the residual horizontal movement of concave landform can be expressed as but not limited to:

[0122]

[0123] Among them, I b1 ,I z and I b2 are the residual tilts of the left cavity area, the middle cavity area, and the right cavity area, respectively, and U′ 0 (x, t) is the residual horizontal movement value at time t at point x on the mountain surface.

[0124] The calculation formula of the residual horizontal movement of convex landform can be expressed as but not limited to:

[0125]

[0126] As a possible implementation method, through the above analysis, the embodiment of the present application can obtain a mountain surface residual deformation prediction model of the target mountain surface, wherein the residual subsidence calculation formula can be expressed as but not limited to:

[0127]

[0128] The calculation formula of the residual horizontal movement of concave landform can be expressed as but not limited to:

[0129]

[0130] The calculation formula of the residual horizontal movement of convex landform can be expressed as but not limited to:

[0131]

[0132] The working principle of the method for predicting residual deformation of surface goaf areas in mountainous areas proposed in the embodiment of the present application is described in detail below in conjunction with a specific embodiment.

[0133] in, Figure 5 The present invention is a flow chart of the working principle of a method for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application.

[0134] Step S501: Obtain the residual deformation value caused by flat ground mining through the probability integral method prediction model of variable mining thickness.

[0135] Step S502: obtaining the initial residual deformation value of the target mountain surface through the initial mountain surface residual deformation prediction model.

[0136] Step S503: correcting the initial residual deformation value to construct a mountain surface residual deformation prediction model.

[0137] According to the method for predicting residual deformation of the goaf area on the surface of mountainous area proposed in the embodiment of the present application, the geomorphic data of the surface of the target mountainous area can be obtained by dividing the different mountainous terrains and landforms, and the initial residual deformation value of the surface of the target mountainous area can be obtained by the residual deformation value caused by flat land mining and the initial mountainous surface residual deformation prediction model, and the mountainous surface residual deformation prediction model of the target mountainous area can be constructed, and the problem of residual deformation of the surface of the mountainous area can be solved by establishing a new method for obtaining the inclination angle of the smoothed terrain trend surface and the mountainous surface characteristic coefficient, and a calculation model is provided for the residual deformation of the goaf area on the surface of the mountainous area. Thus, the problems in the related technology that the measurement is difficult, the measurement accuracy is low, and it is easy to be disturbed by external factors, and the existing mathematical models of mining residual sinking in mining areas are mostly established in plain mining areas, and the limited mathematical models hinder the prediction of residual deformation of the surface of the mountainous area, and cannot meet the requirements of the prediction of residual deformation of the surface of the mountainous area under complex terrain conditions.

[0138] Next, a prediction device for residual deformation of a mined-out area on a mountainous surface proposed in accordance with an embodiment of the present application will be described with reference to the accompanying drawings.

[0139] Figure 6 A schematic block diagram of a device for predicting residual deformation of a surface goaf in a mountainous area provided according to an embodiment of the present application.

[0140] like Figure 6 As shown, the prediction device 10 for residual deformation of surface goaf in mountainous area is applied in the model building stage, wherein the device 10 includes: a first acquisition module 100, a first generation module 200 and a construction module 300.

[0141] The first acquisition module 100 is used to divide different mountain terrains and obtain division results, and obtain at least one landform data in the target mountain surface based on the division results to obtain at least one target landform data in the target mountain surface.

[0142] The first generation module 200 is used to obtain the initial residual deformation value of the target mountain surface through the residual deformation value caused by flat land mining and at least one target landform data, using the initial mountain surface residual deformation prediction model, wherein the initial residual deformation value includes at least one of an initial residual subsidence value and an initial residual horizontal movement value.

[0143] Construction module 300 is used to correct the initial residual deformation value, and construct a mountain surface residual deformation prediction model of the target mountain surface based on the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model, so as to use the mountain surface residual deformation prediction model to predict the residual deformation of the goaf area of ​​the target mountain surface, wherein the mountain surface residual deformation prediction model includes at least one of the residual subsidence calculation formula, the concave landform residual horizontal movement calculation formula and the convex landform residual horizontal movement calculation formula.

[0144] Optionally, in one embodiment of the present application, the first generating module 200 includes: an acquiring unit and a generating unit.

[0145] The acquisition unit is used to obtain the residual deformation value caused by flat ground mining.

[0146] The generating unit is used to input at least one target landform data into the initial mountain surface residual deformation prediction model based on the residual deformation value to obtain the initial residual deformation value of the target mountain surface.

[0147] Optionally, in one embodiment of the present application, the construction module 300 includes: a first judgment unit and a first correction unit.

[0148] The first judgment unit is used to judge whether there are positive and negative slope change points and / or concave and convex landform change points on the surface of the target mountain area.

[0149] The first correction unit is used to correct the initial residual deformation value at the positive and negative slope change point and / or the concave-convex landform change point when there are positive and negative slope change points and / or concave-convex landform change points on the surface of the target mountain area to obtain a corrected initial residual deformation value.

[0150] Optionally, in one embodiment of the present application, the construction module 300 includes: a second judgment unit and a second correction unit.

[0151] The second judgment unit is used to judge whether there is a convex landform on the surface of the target mountain area.

[0152] The second correction unit is used to correct the initial residual horizontal movement value in the convex landform when there is a convex landform on the surface of the target mountain area, so as to obtain a corrected initial residual deformation value.

[0153] Optionally, in one embodiment of the present application, the residual sinking calculation formula may be, but is not limited to:

[0154]

[0155] Among them, W′ 0 (x, t) is the residual subsidence value at the surface x in the mountain area at time t, h n is the block height of the left hole area, n 0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, n0′ is the number of small blocks in the right hollow area, h n′ The height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the smoothed surface point, m 1 It is the residual activatable thickness in the middle area.

[0156] The calculation formula of the residual horizontal movement of concave landform can be expressed as but not limited to:

[0157]

[0158] Among them, I b1 ,I z and I b2 are the residual tilts of the left cavity area, the middle cavity area, and the right cavity area, respectively, and U′ 0 (x, t) is the residual horizontal movement value at time t at point x on the mountain surface.

[0159] The calculation formula of the residual horizontal movement of convex landform can be expressed as but not limited to:

[0160]

[0161] The above embodiment describes the model building phase, and the following describes an embodiment of the use phase.

[0162] Figure 7 The present invention is a flowchart of another method for predicting residual deformation of surface goaf areas in mountainous areas provided according to an embodiment of the present application.

[0163] like Figure 7 As shown, the prediction method of residual deformation of the surface goaf area in the mountainous area is applied to the model use stage, wherein the method comprises the following steps:

[0164] In step S701, different mountainous terrains are divided to obtain division results, and at least one actual landform data in the target mountainous area surface is obtained based on the division results.

[0165] It can be understood that the embodiments of the present application can divide different mountain terrains and landforms, such as dividing them into positive slopes, reverse slopes and multi-faceted slopes according to the different surface inclinations and changes in inclination angles, and the present application does not make specific restrictions; it can also be divided into convex landforms and concave landforms according to the morphological characteristics of the terrain and geomorphological principles, and the present application does not make specific restrictions.

[0166] Furthermore, the embodiment of the present application can obtain at least one actual landform data in the target mountain area surface based on the division result.

[0167] For example, the embodiments of the present application may be selected as follows Figure 2 An actual multi-faceted slope topography is shown in Figure 2 In the embodiment of the present application, the left side of the mined coal seam is used as the coordinate origin. The two sides of the figure are unmined coal pillar areas, with coordinate areas of -300m to 0m and 1200m-1500m, and the middle is a limited mining area with a coordinate area of ​​0m-1200m. Among them, the multi-faceted slope in the embodiment of the present application includes two positive slopes and one reverse slope, a convex landform and a concave landform. Further, the two positive slopes and one reverse slope are 600m long, the positive slope coordinate area is -300m-300m and 900m-1500m, the reverse slope coordinate area is 300m-900m, the convex landform coordinate area is 1400m-1500m, and the concave landform coordinate area is 200m-400m. Further, after measurement and calculation, the trend surface inclination angle of the leftmost positive slope in the embodiment of the present application is 15°, the trend surface inclination angle of the rightmost positive slope is 20°, and the trend surface inclination angle of the middle reverse slope is -25°.

[0168] In step S702, at least one actual geomorphic data is input into a pre-constructed mountain surface residual deformation prediction model to obtain a residual deformation value of the target mountain surface, wherein the pre-constructed mountain surface residual deformation prediction model is obtained from at least one actual geomorphic data of the target mountain surface.

[0169] As a possible implementation method, the embodiment of the present application can input at least one actual landform data into a pre-constructed mountain surface residual deformation prediction model to obtain the residual deformation value of the target mountain surface.

[0170] For example, the residual deformation value of the target mountain surface is obtained based on the pre-constructed mountain surface residual deformation prediction model in the embodiment of the present application, and the result diagram can be shown as follows: Figure 8 shown.

[0171] Furthermore, in order to verify the prediction model of residual deformation of mountain surface and the slope angle α of the smoothed terrain trend surface, the embodiment of the present application x and mountain surface characteristic coefficient D x The value method of the target mountain area with concave landform on the left is the surface characteristic coefficient D x The value is -0.3, the surface characteristic coefficient D of the target mountain area with convex landform on the right side x The value is +0.4.

[0172] Furthermore, in order to verify the reliability of the residual deformation prediction model for the mountainous surface, the embodiment of the present application compares and analyzes the residual deformation value of the plain mining with the residual deformation value of the target mountainous surface. Fig. 9 As shown, for residual subsidence, except for the negative subsidence in the horizontal compression part of the valley concave landform, the residual subsidence values ​​of the mountain surface in other places are larger than those of the plain; for residual horizontal movement, under the positive slope condition, the surface horizontal movement caused by the movement of the rock layer and the horizontal component of the slope slip are superimposed, and the residual horizontal movement increases. Under the reverse slope condition, the surface horizontal movement caused by the movement of the rock layer is partially offset by the surface slip, and the residual horizontal movement decreases. It can be concluded that the calculation results of the mountain surface residual deformation prediction model constructed by the embodiment of the present application are consistent with the actual law and can be used in actual mining.

[0173] According to the method for predicting residual deformation of goaf areas on the surface of mountainous areas proposed in the embodiment of the present application, the actual geomorphic data of the target mountainous surface can be obtained by dividing different mountainous terrains and landforms, and the actual geomorphic data can be input into a pre-constructed model for predicting residual deformation of the mountainous surface to obtain the residual deformation value of the target mountainous surface, thereby verifying the stability and reliability of the model, which can well solve the problem of residual deformation of the mountainous surface and provide a practical calculation model for residual deformation of goaf areas on the surface of mountainous areas. Thus, the problems in the related technology such as high measurement difficulty, low measurement accuracy, easy interference by external factors, and the existing mathematical models of residual subsidence in mining areas are mostly established in plain mining areas, and the limited mathematical models hinder the prediction of residual deformation of the mountainous surface and cannot meet the requirements for predicting residual deformation of the mountainous surface under complex terrain conditions are solved.

[0174] Next, a prediction device for residual deformation of a mined-out area on a mountainous surface proposed in accordance with an embodiment of the present application will be described with reference to the accompanying drawings.

[0175] Fig.10 A schematic block diagram of another device for predicting residual deformation of a mined-out area on a mountainous surface provided according to an embodiment of the present application.

[0176] like Fig.10 As shown, the prediction device 20 for residual deformation of surface goaf areas in mountainous areas is applied in the model use stage, wherein the device 20 includes: a second acquisition module 400 and a second generation module 500.

[0177] The second acquisition module 400 is used to divide the landforms of different mountainous areas to obtain a division result, and obtain at least one actual landform data on the surface of the target mountainous area based on the division result.

[0178] The second generation module 500 is used to input at least one actual geomorphic data into a pre-constructed mountain surface residual deformation prediction model to obtain a residual deformation value of the target mountain surface, wherein the pre-constructed mountain surface residual deformation prediction model is obtained from at least one actual geomorphic data of the target mountain surface.

[0179] It should be noted that the above explanation of the embodiment of the method for predicting residual deformation of a goaf area on a mountainous surface is also applicable to the device for predicting residual deformation of a goaf area on a mountainous surface of this embodiment, and will not be repeated here.

[0180] According to the prediction device for residual deformation of goaf areas on the surface of mountainous areas proposed in the embodiment of the present application, the actual geomorphic data of the surface of the target mountainous area can be obtained by dividing the different mountainous terrains and landforms, and the actual geomorphic data can be input into a pre-constructed prediction model for residual deformation of the surface of the mountainous area to obtain the residual deformation value of the surface of the target mountainous area, thereby verifying the stability and reliability of the model, which can well solve the problem of residual deformation of the surface of mountainous areas, and provide a practical calculation model for residual deformation of goaf areas on the surface of mountainous areas. Thus, the problems in the related technology that the measurement is difficult, the measurement accuracy is low, and it is easy to be disturbed by external factors, and the mathematical models of residual subsidence in existing mining areas are mostly established in plain mining areas, and the limited mathematical models hinder the prediction of residual deformation of the surface of mountainous areas, and cannot meet the requirements for the prediction of residual deformation of the surface of mountainous areas under complex terrain conditions.

[0181] Fig.11 This is a schematic diagram of the structure of an electronic device provided according to an embodiment of the present application. The electronic device may include:

[0182] A memory 1101 , a processor 1102 , and a computer program stored in the memory 1101 and executable on the processor 1102 .

[0183] When the processor 1102 executes the program, the method for predicting residual deformation of a mined-out area on a mountainous surface provided in the above-mentioned embodiment is implemented.

[0184] Furthermore, the electronic device further comprises:

[0185] The communication interface 1103 is used for communication between the memory 1101 and the processor 1102 .

[0186] The memory 1101 is used to store computer programs that can be executed on the processor 1102 .

[0187] The memory 1101 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0188] If the memory 1101, the processor 1102 and the communication interface 1103 are implemented independently, the communication interface 1103, the memory 1101 and the processor 1102 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Fig.11 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0189] Optionally, in a specific implementation, if the memory 1101, the processor 1102 and the communication interface 1103 are integrated on a chip, the memory 1101, the processor 1102 and the communication interface 1103 can communicate with each other through an internal interface.

[0190] The processor 1102 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0191] The embodiment of the present application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-mentioned method for predicting residual deformation of a mined-out area on a mountainous surface.

[0192] The embodiment of the present application also provides a computer program product, including a computer program, which, when executed, implements the above-mentioned method for predicting residual deformation of a surface goaf in a mountainous area.

[0193] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0194] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0195] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0196] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or N wirings (electronic devices), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways as necessary and then storing it in a computer memory.

[0197] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented by hardware, as in another embodiment, it can be implemented by any one or a combination of multiple of the following technologies known in the art: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0198] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0199] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0200] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for predicting residual deformation of surface goaf in mountainous areas, characterized in that: Applied to the model building stage, wherein the method comprises the following steps: Dividing different mountainous terrains to obtain a division result, and obtaining at least one landform data on the surface of the target mountainous area based on the division result; Obtaining the initial residual deformation value of the target mountain surface by using the residual deformation value caused by flat land mining and the at least one target landform data and using the initial mountain surface residual deformation prediction model; The initial residual deformation value is corrected, and a mountain surface residual deformation prediction model of the target mountain surface is constructed according to the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model, so as to use the mountain surface residual deformation prediction model to predict the residual deformation of the goaf area of ​​the target mountain surface, wherein the mountain surface residual deformation prediction model includes at least one of a residual subsidence calculation formula, a concave landform residual horizontal movement calculation formula, and a convex landform residual horizontal movement calculation formula; Wherein, the residual subsidence calculation formula is expressed as: Where W0′(x, t) is the residual subsidence value at the surface x in the mountainous area at time t, h n is the block height of the left hole area, n0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, m0′ is the number of small blocks in the right hollow area, h n ′ is the height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the surface point after smoothing, m1 is the residual reactivatable thickness in the middle area; The calculation formula for the residual horizontal movement of the concave landform is expressed as: Among them, I b1 ,I z and I b2 are the residual tilts of the left hole area, the middle hole area, and the right hole area, respectively; U0′(x, t) is the residual horizontal movement value at time t at point x on the mountain surface; The calculation formula for the residual horizontal movement of the convex landform is expressed as:

2. The method for predicting residual deformation of surface goaf in mountainous areas according to claim 1 is characterized in that: The residual deformation value caused by flat land mining and the at least one target landform data are used to obtain the initial residual deformation value of the target mountain surface by using the initial mountain surface residual deformation prediction model, including: Obtain the residual deformation value caused by flat ground mining; Based on the residual deformation value, the at least one target landform data is input into the initial mountain surface residual deformation prediction model to obtain the initial residual deformation value of the target mountain surface.

3. The method for predicting residual deformation of surface goaf in mountainous areas according to claim 1, characterized in that: The correcting the initial residual deformation value comprises: Determine whether there are positive and negative slope change points and / or concave and convex landform change points on the surface of the target mountain area; If there are positive and negative slope change points and / or concavo-convex landform change points on the target mountain surface, the initial residual deformation values ​​at the positive and negative slope change points and / or concavo-convex landform change points are corrected to obtain corrected initial residual deformation values.

4. The method for predicting residual deformation of surface goaf in mountainous areas according to claim 1, characterized in that: The correcting the initial residual deformation value comprises: Determine whether there is a convex landform on the surface of the target mountain area; If there is a convex landform on the surface of the target mountain area, the initial residual horizontal movement value in the convex landform is corrected to obtain a corrected initial residual deformation value.

5. A method for predicting residual deformation of a goaf area on a mountainous surface, using the method for predicting residual deformation of a goaf area on a mountainous surface as described in any one of claims 1 to 4, characterized in that: Applied to the model use phase, wherein the method comprises the following steps: Dividing different mountainous terrains to obtain a division result, and obtaining at least one actual landform data on the surface of the target mountainous area based on the division result; The at least one actual geomorphic data is input into a pre-constructed mountain surface residual deformation prediction model to obtain the residual deformation value of the target mountain surface, wherein the pre-constructed mountain surface residual deformation prediction model is obtained from at least one actual geomorphic data of the target mountain surface.

6. A prediction device for residual deformation of a goaf area on a mountainous surface, using the prediction method for residual deformation of a goaf area on a mountainous surface as described in any one of claims 1 to 4, characterized in that: Applied to the model building stage, wherein the device comprises: An acquisition module, used to divide different mountain terrains and landforms to obtain a division result, and obtain at least one landform data in the target mountain surface based on the division result to obtain at least one target landform data in the target mountain surface; A generating module, configured to obtain an initial residual deformation value of the surface of a target mountain area by using the residual deformation value caused by flat land mining and the at least one target landform data and an initial mountain area surface residual deformation prediction model, wherein the initial residual deformation value includes at least one of an initial residual subsidence value and an initial residual horizontal movement value; A construction module is used to correct the initial residual deformation value, and to construct a mountain surface residual deformation prediction model of the target mountain surface according to the corrected initial residual deformation value and the initial mountain surface residual deformation prediction model, so as to use the mountain surface residual deformation prediction model to predict the residual deformation of the goaf area of ​​the target mountain surface, wherein the mountain surface residual deformation prediction model includes at least one of a residual subsidence calculation formula, a concave landform residual horizontal movement calculation formula, and a convex landform residual horizontal movement calculation formula; Wherein, the residual subsidence calculation formula is expressed as: Where W0′(x, t) is the residual subsidence value at the surface x in the mountainous area at time t, h n is the block height of the left hole area, n0 is the number of small blocks in the left hole area, W e (x) is the unit subsidence caused by unit mining, x is the horizontal coordinate of any point on the surface, s is the horizontal coordinate of the mining unit, n0′ is the number of small blocks in the right hollow area, h n ′ is the height of the block in the right hole area, D x is the surface characteristic coefficient of the mountainous area, P(x) is the slip influence function, α x is the topographic trend surface inclination angle of the surface point after smoothing, m1 is the residual reactivatable thickness in the middle area; The calculation formula for the residual horizontal movement of the concave landform is expressed as: Among them, I b1 ,I z and I b2 are the residual tilts of the left hole area, the middle hole area, and the right hole area, respectively; U0′(x, t) is the residual horizontal movement value at time t at point x on the mountain surface; The calculation formula for the residual horizontal movement of the convex landform is expressed as:

7. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for predicting residual deformation of a surface goaf in a mountainous area as described in any one of claims 1 to 4 or the method for predicting residual deformation of a surface goaf in a mountainous area as described in claim 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement the method for predicting residual deformation of a goaf area on a mountainous surface as described in any one of claims 1 to 4 or the method for predicting residual deformation of a goaf area on a mountainous surface as described in claim 5.

9. A computer program product, characterized in that It comprises a computer program, which, when executed, is used to implement the method for predicting residual deformation of a goaf area on a mountainous surface as claimed in any one of claims 1 to 4 or the method for predicting residual deformation of a goaf area on a mountainous surface as claimed in claim 5.

Images