CN111798052B - Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space - Google Patents
Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space Download PDFInfo
- Publication number
- CN111798052B CN111798052B CN202010637174.4A CN202010637174A CN111798052B CN 111798052 B CN111798052 B CN 111798052B CN 202010637174 A CN202010637174 A CN 202010637174A CN 111798052 B CN111798052 B CN 111798052B
- Authority
- CN
- China
- Prior art keywords
- subsidence
- basin
- volume
- area
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A dynamic prediction method for three-dimensional spatial information of a high-diving-level coal mining subsidence ponding area is suitable for dynamic prediction of a surface mining subsidence ponding area of a high-diving-level mining area. Firstly, acquiring working face information and geological mining conditions of a mining area, dynamically predicting subsidence basins by utilizing a probability integration method according to the mining progress of the working face, and calculating the volumes of the subsidence basins with different contour line sections; establishing a water quantity balance iterative equation according to hydrological water resource data of the coal mining subsidence area; calculating the water accumulation amount of the coal mining subsidence area according to an iterative equation, and establishing a relation model between the volume of the subsidence basin and the volume of water accumulated in the basin; and then predicting three-dimensional dynamic space information of the subsidence basin. The method has simple steps and low cost, can dynamically predict the three-dimensional space information of the coal mining subsidence ponding area along with the advancing of the working face and the lapse of time, is favorable for mastering the dynamic change condition of water resources in the subsidence area, and is suitable for the ecological environment comprehensive treatment of the high-diving-level coal mining subsidence area.
Description
Technical Field
The invention mainly relates to a dynamic prediction method for three-dimensional space information of a high-diving-level coal mining subsidence ponding area, which is suitable for dynamic prediction of water resources and comprehensive treatment of ecological environment of the high-diving-level coal mining subsidence area.
Background
Along with the advance of a working face, the ground surface subsidence basin is influenced by the coupling of factors such as underground water seepage, atmospheric precipitation, the exchange of ground surface rivers and ponding water in the subsidence area, and the ground surface subsidence basin is easy to form large-scale ponding and becomes a typical hydrogeological disaster caused by mining subsidence in the eastern high-diving-level mining area of China. Along with the long-time repeated exploitation, the ponding range of the subsidence area is continuously expanded, the large-scale exploitation subsidence area forms a surface subsidence pond, the ponding in a large range can irreversibly change the original ecological system of the subsidence area, and a new aquatic ecological system and a new land-water composite ecological system are formed. When the original ecological environment is destroyed, the cultivated land in the subsidence area is reduced sharply due to coal mining subsidence, surface buildings and roads are damaged seriously, the burden of the construction land and the human ground shield are aggravated by surface subsidence and water accumulation caused by mining, the time span of the water accumulation in the subsidence area is long, the submerging range is wide, the influence factors are multiple, and a series of problems are brought to the comprehensive treatment and sustainable development of the subsidence area.
At present, the geographical spatial information data of a coal mining subsidence area are mainly obtained by using a method for surveying and mapping on the spot and predicting mining subsidence, however, the surface morphology of the mining subsidence area and the water accumulation condition in a subsidence basin are influenced by a plurality of factors and dynamically change, the observation data obtained by the existing method is lack of timeliness, the three-dimensional spatial information of the water accumulation in the surface subsidence area is difficult to predict in time, and the accuracy of planning and design of the coal mining subsidence area is influenced, so that a series of problems are brought, for example, the ecological management of the water accumulation area is contradicted with the water accumulation range evolution, and a building (structure) is influenced by the water accumulation in the surface subsidence and the subsidence area and even.
Disclosure of Invention
Aiming at the defects of the technology, the dynamic prediction method for the three-dimensional space information of the high-diving-level coal mining subsidence ponding area is simple in process, low in cost and capable of dynamically predicting the three-dimensional space information of the coal mining subsidence ponding area along with the advancing of a working face.
In order to realize the defects of the technology, the method for dynamically predicting the three-dimensional space information of the coal mining subsidence ponding area of the high diving space is provided, and the dynamic prediction is synchronously carried out along with the advancing of a working face until the water resources are balanced, and the method comprises the following steps:
a. acquiring geological mining condition parameters of a subsidence basin area, establishing a three-dimensional coordinate system of a mining subsidence space of the subsidence basin area, and acquiring a dynamic subsidence predicted value of the subsidence basin area by a mining subsidence dynamic prediction method;
b. carrying out interpolation processing on the dynamic subsidence predicted values of the subsidence basins by utilizing a Krigin interpolation method, sequentially generating contour lines of the subsidence basins in the subsidence basin areas, and calculating the volumes of the subsidence basins corresponding to different contour line sections;
c. collecting hydrological water resource data of a coal mining subsidence area, and establishing a water balance iterative equation according to a water balance principle so as to obtain the volume of accumulated water in a subsidence basin of the coal mining subsidence area at different moments;
d. calculating the volume of accumulated water in the subsidence basin at different moments according to a water quantity balance iterative equation, and establishing a relation model between the volume of the subsidence basin and the volume of accumulated water in the basin by combining the volumes of the subsidence basins corresponding to different contour line sections at the moments;
e. according to the solution of the relational model between the volume of the subsidence basin and the accumulated water volume in the basin, the three-dimensional dynamic space information at different moments pushed along with the working surface can be predicted, wherein the three-dimensional dynamic space information comprises the accumulated water depth, the accumulated water area, the accumulated water range, the accumulated water volume and the maximum storage capacity of the subsidence basin.
The method comprises the following specific steps:
s1, obtaining geological mining condition parameters according to the mining area working face information: the method comprises the following steps of (1) mining advancing distance v, main influence angle tangent tan beta, working face average mining thickness m, coal seam inclination angle alpha, sinking coefficient q and main influence radius r; establishing a three-dimensional coordinate system of a mining subsidence space: the projection of the inflection point of the sinking curve of the main section on the ground surface is a coordinate origin o, the x axis points to the direction of the goaf along the ground surface and is parallel to the direction of the coal seam, and the y axis passes through the coordinate origin and is vertical to the x axis and the z axis and downward; adopting a probabilistic integration method based on a Knothe time function to carry out mining subsidence dynamic prediction to obtain a dynamic subsidence prediction value of a subsidence basin;
s2, according to the dynamic subsidence predicted value of the subsidence basin predicted by the probability integration method, carrying out interpolation processing on the dynamic subsidence predicted value of the subsidence basin by using a Krigin interpolation method, sequentially generating a subsidence basin contour line with equal height distance d from the position where the subsidence value is 10mm, wherein the unit mm is obtained by using a formula: d ═ Wm- (10+ nd) calculating the difference d' between the bottom of the sinking basin and the lowest horizontal section in the basin, in mm, n being the number of the horizontal sections except the edge of the sinking basin, WmAt a certain momentCalculating the maximum sinking value of the table by an analytic method, wherein the maximum sinking value of the table is 10, 10+ d, 10+2d1，S2，S3...SnCalculating the volumes of the sinking basins corresponding to different contour line sections by utilizing the section areas of the contour lines of the sinking basins;
s3, analyzing influence factors of the accumulated water range in the subsidence basin, and collecting hydrological water resource data of the coal mining subsidence area, wherein the hydrological water resource data comprises precipitation, evaporation capacity, underground water seepage capacity, artificial water intake and surface river and accumulated water exchange capacity of the subsidence area; establishing a water balance iterative model according to a water balance principle, and calculating the volume of accumulated water in the coal mining subsidence area at different moments T;
s4, according to the accumulated water volume in the subsidence basin at different moments and the subsidence basin volume corresponding to the sections of different isolines at the moments, establishing a relation model between the subsidence basin volume and the accumulated water volume in the basin, solving, calculating to obtain the isoline of the subsidence basin equal to the accumulated water volume in the subsidence basin, and further predicting three-dimensional dynamic space information such as the depth, the area, the range, the volume, the maximum storage capacity and the like of the accumulated water in the subsidence basin at different moments T along with the advancing of the working surface.
According to geological mining condition parameters, using a formula:
calculating the dynamic subsidence predicted value w (x) of subsidence basini,yiT), in which: w (x)i,yiAnd T) is the dynamic subsidence predicted value of any point of the subsidence basin ground surface at the time T, which is unit meter; m is the average mining thickness of the working face and is unit meter; alpha is the coal bed inclination angle; q is a sinking coefficient; r is the major radius of influence in meters, c is the time coefficient, e is a natural constant, about 2.718,
for double integral variables, there is no physical significance and can be interpreted as:
tau is an integral variable in the x direction and the y direction respectively; d is the extent of the sink basin at time T, as does dx, dy in the integral.
Using the formula:
calculating the horizontal cross-sectional area of each contour line of the sink basin, wherein:
is the coordinates of the isopleth point of the boundary of the subsidence basin with the subsidence value of 10+ nd (n is a positive integer and unit mm).
Using the formula:
calculating the corresponding sinking basin volumes of the line sections with different equal values, wherein:
the volume of a sinking basin corresponding to the section of the contour line 10+ (k-1) d; v1Is the volume between contour 10 and the contour 10+ d cross-section; v2Volume between the section of contour 10+ d and contour 10+2d, VnThe volume between the section of contour 10+ (n-1) d and contour 10+ nd; v' is the contour line 10+ nd and the maximum sinking point WmThe volume in between.
The water quantity balance iterative model is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit is the manual extraction amount.
The mathematical model between the volume of the subsidence basin and the volume of the accumulated water in the basin is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
is the corresponding sinking basin volume of the section of the contour line 10+ (k-1) d in unit of cube.
Predicting the three-dimensional space information of the surface water area at different moments T along with the advancing of the working surface, and obtaining the three-dimensional space information by the following formula:
in the formula:
dynamically predicting the depth of the ponding water in the subsidence basin at the time T in units of meters; wmThe maximum sinking value of the sinking basin at the moment T is unit meter;
dynamically predicting the ponding area of the subsidence basin at the time T, wherein the unit square is the unit square;
is the coordinate of the isoline point with the sinking value of 10+ (k-1) d;
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit cube is the manual output;
the maximum accumulated reservoir capacity of the subsidence basin at the moment T is dynamically predicted, and the unit is cubic; v1Is S and S1Volume between, V2Is S1And S2Volume between, VnIs Sn-1And SnV' is SnAnd maximum sinking point WmThe volume in between; the boundary of the water accumulation range is equal to the sinking value of 10+ (k-1) dCoordinate points on the value line
A closed curve is formed by encircling.
Has the advantages that: the method combines the probability integration method dynamic prediction sinking basin volume with the water quantity balance iterative equation, establishes a mathematical model between the sinking basin volume and the water accumulation volume in the basin, dynamically predicts the water accumulation depth, the water accumulation area, the water accumulation volume, the maximum storage capacity and the water accumulation range in the basin along with the advancing of a working surface and the passage of time, has certain timeliness of a prediction result, and can effectively avoid the problem that the comprehensive treatment of the sinking basin cannot follow up or be repeatedly treated in real time due to the lack of timeliness of data.
In the process of dynamically predicting the three-dimensional spatial information of the coal mining subsidence ponding area, the acquired hydraulic information in the basin of the coal mining subsidence ponding area can provide reference for comprehensive treatment of the coal mining subsidence area, can assist in protecting the ecology and the surrounding geographic environment, and the mining subsidence dynamic prediction result can guide reclamation of the mining area; the ponding range and water depth information can provide data for the establishment of lake/wetland ecosystem in the water area of the coal mining subsidence area; the volume of the subsidence basin can provide reference for the establishment of water storage engineering. The method has simple steps and low cost, and can dynamically predict the three-dimensional space information of the coal mining subsidence ponding area along with the advancing of the working face and the lapse of time.
Drawings
FIG. 1 is a schematic diagram of the water circulation of the subsidence ponding area of the high-diving-level coal mining area of the invention.
FIG. 2 is a schematic diagram of a three-dimensional space coordinate system of the high-diving-level coal mining area sinking ponding area probability integration method.
FIG. 3 is a flow chart of dynamic prediction of three-dimensional spatial information of a subsidence ponding area of a high-diving-level coal mining area according to the invention.
Fig. 4 is a line diagram for dynamically predicting three-dimensional spatial information according to an embodiment of the present invention.
FIG. 5 is a schematic diagram of dynamically predicting the variation of the water accumulation range according to the embodiment of the present invention.
In the figure: 1-working surface, 2-mining depth, 3-atmospheric precipitation, 4-average daily evaporation capacity of water surface evaporation, 5-surface runoff yield, 6-water exchange, 7-diving space, 8-space three-dimensional coordinate system, 9-water accumulation depth and 10-inflection offset.
Detailed Description
The invention is further described with reference to the following specific examples:
as shown in FIG. 3, the dynamic prediction method of the three-dimensional spatial information of the coal mining subsidence ponding area of the high-diving space of the invention is synchronously carried out along with the advancing of the working face until the water resources are balanced, and the steps are as follows:
a. acquiring geological mining condition parameters of a subsidence basin area, establishing a three-dimensional coordinate system of a mining subsidence space of the subsidence basin area, and acquiring a dynamic subsidence predicted value of the subsidence basin area by a mining subsidence dynamic prediction method;
b. carrying out interpolation processing on the dynamic subsidence predicted values of the subsidence basins by utilizing a Krigin interpolation method, sequentially generating contour lines of the subsidence basins in the subsidence basin areas, and calculating the volumes of the subsidence basins corresponding to different contour line sections;
c. collecting hydrological water resource data of a coal mining subsidence area, and establishing a water balance iterative equation according to a water balance principle so as to obtain the volume of accumulated water in a subsidence basin of the coal mining subsidence area at different moments;
d. calculating the volume of accumulated water in the subsidence basin at different moments according to a water quantity balance iterative equation, and establishing a relation model between the volume of the subsidence basin and the volume of accumulated water in the basin by combining the volumes of the subsidence basins corresponding to different contour line sections at the moments;
e. according to the solution of the relational model between the volume of the subsidence basin and the accumulated water volume in the basin, the three-dimensional dynamic space information at different moments pushed along with the working surface can be predicted, wherein the three-dimensional dynamic space information comprises the accumulated water depth, the accumulated water area, the accumulated water range, the accumulated water volume and the maximum storage capacity of the subsidence basin.
The method comprises the following specific steps:
s1, obtaining geological mining condition parameters according to the mining area working face information: the method comprises the following steps of (1) mining advancing distance v, main influence angle tangent tan beta, working face average mining thickness m, coal seam inclination angle alpha, sinking coefficient q and main influence radius; establishing a three-dimensional coordinate system of a mining subsidence space: the projection of the inflection point of the sinking curve of the main section on the ground surface is a coordinate origin o, the x axis points to the direction of the goaf along the ground surface and is parallel to the direction of the coal seam, and the y axis passes through the coordinate origin and is vertical to the x axis and the z axis and downward; adopting a probabilistic integration method based on a Knothe time function to carry out mining subsidence dynamic prediction to obtain a dynamic subsidence prediction value of a subsidence basin;
according to geological mining condition parameters, using a formula:
calculating the dynamic subsidence predicted value w (x) of subsidence basini,yiT), in which: w (x)i,yiAnd T) is the dynamic subsidence predicted value of any point of the subsidence basin ground surface at the time T, which is unit meter; m is the average mining thickness of the working face and is unit meter; alpha is the coal bed inclination angle; q is a sinking coefficient; r is the major radius of influence, in meters; c is a time coefficient; e is a natural constant of about 2.718;
d tau is a double integral variable,
tau is an integral variable in the x direction and the y direction respectively; d is the scope of the sink basin at the time T.
S2, according to the dynamic subsidence predicted value of the subsidence basin predicted by the probability integration method, carrying out interpolation processing on the dynamic subsidence predicted value of the subsidence basin by using a Krigin interpolation method, sequentially generating a subsidence basin contour line with equal height distance d from the position where the subsidence value is 10mm, wherein the unit mm is obtained by using a formula: d ═ Wm- (10+ nd) calculating the difference d' between the bottom of the sinking basin and the lowest horizontal section in the basin, in mm, n being the number of the horizontal sections except the edge of the sinking basin, WmFor the maximum subsidence value of the earth surface at a certain time, isolines with subsidence values of 10, 10+ d, 10+2d., 10+ nd, the unit mm and the corresponding section areas S, S of the isolines of different subsidence basins are sequentially calculated by an analytic method1，S2，S3...SnUsing respective sinking basinsCalculating the volumes of sinking basins corresponding to different contour line sections according to the section areas of the contour lines;
using the formula:
calculating the corresponding sinking basin volumes of the line sections with different equal values, wherein: vb kThe volume of a sinking basin corresponding to the section of the contour line 10+ (k-1) d; v1Is the volume between contour 10 and the contour 10+ d cross-section; v2Volume between the section of contour 10+ d and contour 10+2d, VnThe volume between the section of contour 10+ (n-1) d and contour 10+ nd; v' is the contour line 10+ nd and the maximum sinking point WmThe volume in between.
Using the formula:
calculating the horizontal cross-sectional area of each contour line of the sink basin, wherein:
is the coordinates of the isopleth point of the boundary of the subsidence basin with the subsidence value of 10+ nd (n is a positive integer and unit mm).
S3, analyzing influence factors of the accumulated water range in the subsidence basin, and collecting hydrological water resource data of the coal mining subsidence area, wherein the hydrological water resource data comprises precipitation, evaporation capacity, underground water seepage capacity, artificial water intake and surface river and accumulated water exchange capacity of the subsidence area; establishing a water balance iterative model according to a water balance principle, and calculating the volume of accumulated water in the coal mining subsidence area at different moments T;
the water quantity balance iterative model is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient; sb TThe area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit is the manual extraction amount.
S4, according to the accumulated water volume in the subsidence basin at different moments and the subsidence basin volume corresponding to the sections of different isolines at the moments, establishing a relation model between the subsidence basin volume and the accumulated water volume in the basin, solving, calculating to obtain the isoline of the subsidence basin equal to the accumulated water volume in the subsidence basin, and further predicting three-dimensional dynamic space information such as the depth, the area, the range, the volume, the maximum storage capacity and the like of the accumulated water in the subsidence basin at different moments T along with the advancing of the working surface.
The mathematical model between the volume of the subsidence basin and the volume of the accumulated water in the basin is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
is the corresponding sinking basin volume of the section of the contour line 10+ (k-1) d in unit of cube.
Obtaining the contour line of the subsidence basin with the same volume as the water accumulation in the subsidence basin by solving the relation model, and predicting the three-dimensional space information of the surface water accumulation area at different moments T along with the advancing of the working surface to obtain the three-dimensional space information by the following formula:
in the formula:
dynamically predicting the depth of the ponding water in the subsidence basin at the time T in units of meters; wmThe maximum sinking value of the sinking basin at the moment T is unit meter;
dynamically predicting the ponding area of the subsidence basin at the time T, wherein the unit square is the unit square;
is the coordinate of the isoline point with the sinking value of 10+ (k-1) d;
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit cube is the manual output;
the maximum accumulated reservoir capacity of the subsidence basin at the moment T is dynamically predicted, and the unit is cubic; v1Is S and S1Volume between, V2Is S1And S2Volume between, VnIs Sn-1And SnV' is SnAnd maximum sinking point WmThe volume in between; the boundary of the water accumulation range is a coordinate point on an isoline with a sinking value of 10+ (k-1) d
A closed curve is formed by encircling.
The first embodiment,
As shown in fig. 1 and 2, for example, in a high-diving mining area, a sinking basin is formed on the ground surface by a rectangular working surface 1 with a mining burial depth 2, and a water accumulation area with a water accumulation depth 9 is generated in the sinking basin under the combined action of atmospheric precipitation 3, water surface evaporation daily average evaporation 4, ground surface runoff 5 and groundwater seepage flow exchange 6, wherein the water surface is consistent with an underground diving place 7. Establishing a three-dimensional coordinate system 8 of a mining subsidence space, wherein the projection of the inflection point of a subsidence curve of a main section of the strike on the ground surface is a coordinate origin o, an x axis points to the direction of the goaf along the ground surface and is parallel to the strike direction of the coal bed, a y axis passes through the coordinate origin and is vertical to the x axis, and a z axis is vertically downward; and 10 is inflection point offset.
a. Acquiring the working area information and geological mining conditions of the coal mining subsidence area: the working face 1 belongs to a first mining working face of a mining area, a coal seam mined by the working face is a two-fold lower uniform Shanxi group coal seam, the old top lithology of the coal seam is mainly siltstone, and a small amount of siderite exists locally; secondly, a sand-shale interbed; the pseudoroof lithology mainly comprises mudstone and sandy rock, the old bedrock of the coal bed mainly comprises a gray sand-mud interbedded layer, the main rock of siltstone is directly arranged at the bottom, pyrite crystals exist locally, the comprehensive mechanized coal mining is carried out by adopting a strike long-wall caving coal mining method, the length of a working face 1 is 1411m, the cutting length is 178m, the mining depth 2 during mining of the working face 1 is 360m, the mining advancing distance is about 83 m/month, the average mining thickness m of the working face is 3.8m, the coal seam inclination angle α is 8 °, the subsidence coefficient q is 0.98, the main influence angle tan β is 2.16, and the main influence radius r is 367 m. The working face 1 starts to be mined from 10 months in 2013, the harvesting time is 2015 and 3 months, mining subsidence dynamic prediction based on the Knothe time function is carried out on the working face 1 according to the obtained data information, and a subsidence basin dynamic subsidence prediction value w (x) is obtainedi,yi,T)。
b. According to the prediction result of the probability integration method, interpolation and gridding processing are carried out on the subsidence basin prediction data, isolines with the interval d being 10mm are generated, wherein the subsidence value of 10mm is taken as the minimum value of the isolines, namely the subsidence basin boundary, the areas S and S of the sections of the isolines are respectively counted1，S2，S3...SnBy predicting the surface of the mining process of the working surface by the same method and outputting a corresponding contour line, the volume of the subsidence basin corresponding to any contour line section in the basin along with the time can be obtained, as shown in the following formula:
in the formula:
the volume of a sinking basin corresponding to the section of the contour line 10+ (k-1) d is a unit cube; v1Is the volume between contour 10 and the contour 10+ d cross-section; v2Volume between the section of contour 10+ d and contour 10+2d, VnThe volume between the section of contour 10+ (n-1) d and contour 10+ nd; v' is the contour line 10+ nd and the maximum sinking point WmThe volume in between.
c. Analyzing influence factors of a water accumulation range in a subsidence basin, collecting hydrological water resource data of a coal mining subsidence area, establishing a water balance iterative equation according to a water balance principle, and calculating the water accumulation volume of the coal mining subsidence area, wherein the water balance iterative equation is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit cube is the manual output;
d. according to the volume of accumulated water in the subsidence basin at different moments and the volume of the subsidence basin corresponding to different contour line sections at the moments, a relation model between the volume of the subsidence basin and the volume of the accumulated water in the basin is established as follows:
in the formula:
the volume of accumulated water is predicted for sinking the basin,a unit cube;
the volume of a sinking basin corresponding to the contour line 10+ (k-1) d is unit cubic;
obtaining the contour line of the subsidence basin with the same volume as the water accumulation in the subsidence basin by solving the relation model, and further predicting the three-dimensional space information of the surface water accumulation area at different moments T along with the advancing of the working surface, as shown in the following formula:
in the formula:
dynamically predicting the depth of the ponding water in the subsidence basin at the time T in units of meters; wmThe maximum sinking value of the sinking basin at the moment T is unit meter;
dynamically predicting the ponding area of the subsidence basin at the time T, wherein the unit square is the unit square;
is the coordinate of the isoline point with the sinking value of 10+ (k-1) d;
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit cube is the manual output;
the maximum accumulated reservoir capacity of the subsidence basin at the moment T is dynamically predicted, and the unit is cubic; v1Is S and S1Volume between, V2Is S1And S2Volume between, VnIs Sn-1And SnV' is SnAnd maximum sinking point WmThe volume in between.
The boundary of the water accumulation range is a coordinate point on an isoline with a sinking value of 10+ (k-1) d
A closed curve is formed by encircling.
When the first-month surface water accumulation range is calculated, because the surface subsidence is small and no water is accumulated on the surface, the atmospheric precipitation 3 is not considered, only the confluence and leakage generated by the surface runoff 5 of the subsidence basin area are considered, and the calculated first-month surface confluence of the subsidence basin is about 290m3Resulting in about 2901m2The maximum water depth of the ponding area is 0.14m, the first predicted surface ponding area is used as basic data for calculating the ponding range of the next month, the dynamic prediction of the three-dimensional spatial information is carried out on the working surface according to the method, and the three-dimensional spatial information of 38 months after the working surface starts to be pushed is obtained, as shown in fig. 4 and 5.
Claims (4)
1. A dynamic prediction method for three-dimensional spatial information of a high-diving-level coal mining subsidence ponding area is synchronously carried out along with the advancing of a working face until water resources are balanced, and is characterized by comprising the following steps:
a. acquiring geological mining condition parameters of a subsidence basin area, establishing a three-dimensional coordinate system of a mining subsidence space of the subsidence basin area, and acquiring a dynamic subsidence predicted value of the subsidence basin area by a mining subsidence dynamic prediction method;
according to geological mining condition parameters, using a formula:
calculating the dynamic subsidence predicted value w (x) of subsidence basini,yiT), in which: w (x)i,yiAnd T) is the dynamic subsidence predicted value of any point of the subsidence basin ground surface at the time T, which is unit meter; m is the average mining thickness of the working face and is unit meter; alpha is the coal bed inclination angle; q is a sinking coefficient; r is the major radius of influence, in meters, c is the time coefficient, e is a natural constant, about 2.718,
for double integral variables, there is no physical significance and can be interpreted as:
tau is an integral variable in the x direction and the y direction respectively; d is the range of the sink basin at the time T, as is dx and dy in the integral;
b. carrying out interpolation processing on the dynamic subsidence predicted values of the subsidence basins by using a Krigin interpolation method, sequentially generating contour lines of the subsidence basins in the subsidence basin areas, and calculating the volumes of the subsidence basins corresponding to different contour line sections;
using the formula:
calculating the horizontal cross-sectional area of each contour line of the sink basin, wherein:
is the coordinate of the isopleth point of the subsidence basin boundary with the subsidence value of 10+ nd and n is a positive integer and the unit mm;
using the formula:
calculating the corresponding sinking basin volumes of different contour line sections, wherein:
the volume of a sinking basin corresponding to the section of the contour line 10+ (k-1) d; v1Is the volume between contour 10 and the contour 10+ d cross-section; v2Volume between the section of contour 10+ d and contour 10+2d, VnThe volume between the section of contour 10+ (n-1) d and contour 10+ nd; v' is the contour line 10+ nd and the maximum sinking point WmThe volume in between;
c. collecting hydrological water resource data of a coal mining subsidence area, and establishing a water balance iterative equation according to a water balance principle so as to obtain the volume of accumulated water in a subsidence basin of the coal mining subsidence area at different moments;
d. calculating the volume of accumulated water in the subsidence basin at different moments according to a water quantity balance iterative equation, and establishing a relation model between the volume of the subsidence basin and the volume of accumulated water in the basin by combining the volumes of the subsidence basins corresponding to different contour line sections at the moments;
e. according to the solution of the relation model between the volume of the subsidence basin and the accumulated water volume in the basin, the three-dimensional dynamic space information at different moments pushed along with the working surface can be predicted, wherein the three-dimensional dynamic space information comprises the accumulated water depth, the accumulated water area, the accumulated water range, the accumulated water volume and the maximum storage capacity of the subsidence basin;
predicting the three-dimensional space information of the surface water area at different moments T along with the advancing of the working surface, and obtaining the three-dimensional space information by the following formula:
in the formula:
to dynamically predict the depth of water accumulation in the subsidence basin at time T,unit of meter; wmThe maximum sinking value of the sinking basin at the moment T is unit meter;
dynamically predicting the ponding area of the subsidence basin at the time T, wherein the unit square is the unit square;
is the coordinate of the isoline point with the sinking value of 10+ (k-1) d;
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit cube is the manual output;
the maximum accumulated reservoir capacity of the subsidence basin at the moment T is dynamically predicted, and the unit is cubic; v1Is S and S1Volume between, V2Is S1And S2The volume between the two or more of the two,Vnis Sn-1And SnV' is SnAnd maximum sinking point WmThe volume in between; the boundary of the water accumulation range is a coordinate point on an isoline with a sinking value of 10+ (k-1) d
A closed curve is formed by encircling.
2. The dynamic prediction method for the three-dimensional space information of the coal mining subsidence ponding area of the high diving space according to claim 1 is characterized by comprising the following specific steps:
s1, obtaining geological mining condition parameters according to the mining area working face information: the method comprises the following steps of (1) mining advancing distance v, main influence angle tangent tan beta, working face average mining thickness m, coal seam inclination angle alpha, sinking coefficient q and main influence radius r; establishing a three-dimensional coordinate system of a mining subsidence space: the projection of the inflection point of the sinking curve of the main section on the ground surface is a coordinate origin o, the x axis points to the direction of the goaf along the ground surface and is parallel to the direction of the coal seam, and the y axis passes through the coordinate origin and is vertical to the x axis and the z axis and downward; adopting a probabilistic integration method based on a Knothe time function to carry out mining subsidence dynamic prediction to obtain a dynamic subsidence prediction value of a subsidence basin;
s2, according to the dynamic subsidence predicted value of the subsidence basin predicted by the probability integration method, carrying out interpolation processing on the dynamic subsidence predicted value of the subsidence basin by using a Krigin interpolation method, sequentially generating a subsidence basin contour line with equal height distance d from the position where the subsidence value is 10mm, wherein the unit mm is obtained by using a formula: d ═ Wm- (10+ nd) calculating the difference d' between the bottom of the sinking basin and the lowest horizontal section in the basin, in mm, n being the number of the horizontal sections except the edge of the sinking basin, WmFor the maximum subsidence value of the earth surface at a certain time, isolines with subsidence values of 10, 10+ d, 10+2d., 10+ nd, the unit mm and the corresponding section areas S, S of the isolines of different subsidence basins are sequentially calculated by an analytic method1，S2，S3...SnCalculating the volumes of the sinking basins corresponding to different contour line sections by utilizing the section areas of the contour lines of the sinking basins;
s3, analyzing influence factors of the accumulated water range in the subsidence basin, and collecting hydrological water resource data of the coal mining subsidence area, wherein the hydrological water resource data comprises precipitation, evaporation capacity, underground water seepage capacity, artificial water intake and surface river and accumulated water exchange capacity of the subsidence area; establishing a water balance iterative model according to a water balance principle, and calculating the volume of accumulated water in the coal mining subsidence area at different moments T;
s4, according to the accumulated water volume in the subsidence basin at different moments and the subsidence basin volume corresponding to the sections of different isolines at the moments, establishing a relation model between the subsidence basin volume and the accumulated water volume in the basin, solving, calculating to obtain the isoline of the subsidence basin equal to the accumulated water volume in the subsidence basin, and further predicting three-dimensional dynamic space information such as the depth, the area, the range, the volume, the maximum storage capacity and the like of the accumulated water in the subsidence basin at different moments T along with the advancing of the working surface.
3. The method for dynamically predicting the three-dimensional spatial information of the coal mining subsidence ponding area of the high phreatic ground according to claim 2, characterized by comprising the following steps of: the water quantity balance iterative model is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
the unit is the original water volume of the subsidence basin; Δ t is the time interval; delta P is the daily average precipitation in meters; delta E is the daily average evaporation in meters; alpha is the surface runoff coefficient;
the area of the section of the contour line of the boundary of the subsidence basin with the subsidence of 10mm is the unit square;
the unit square is the ponding area of the subsidence basin;
is the groundwater seepage flow, unit cubic; rΔtThe unit cubic is the accumulated water exchange quantity of the surface river and the subsidence area; mΔtThe unit is the manual extraction amount.
4. The method for dynamically predicting the three-dimensional spatial information of the coal mining subsidence ponding area of the high phreatic ground according to claim 2, characterized by comprising the following steps of: the mathematical model between the volume of the subsidence basin and the volume of the accumulated water in the basin is as follows:
in the formula:
predicting the volume of accumulated water for the subsidence basin, unit cubic;
is the corresponding sinking basin volume of the section of the contour line 10+ (k-1) d in unit of cube.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010637174.4A CN111798052B (en) | 2020-07-03 | 2020-07-03 | Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010637174.4A CN111798052B (en) | 2020-07-03 | 2020-07-03 | Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111798052A CN111798052A (en) | 2020-10-20 |
| CN111798052B true CN111798052B (en) | 2021-06-01 |
Family
ID=72810263
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010637174.4A Active CN111798052B (en) | 2020-07-03 | 2020-07-03 | Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111798052B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112632717A (en) * | 2020-12-10 | 2021-04-09 | 扎赉诺尔煤业有限责任公司 | Method for predicting movement deformation of dam in coal mining subsidence area |
| CN113932877B (en) * | 2021-09-30 | 2023-12-22 | 深圳市中金岭南有色金属股份有限公司凡口铅锌矿 | Karst water level prediction method for mining area and terminal equipment |
| CN113836740B (en) * | 2021-09-30 | 2024-03-19 | 安徽大学 | Historical space information calculation method for high-water-level coal mining subsidence water accumulation area |
| CN116703244B (en) * | 2023-08-02 | 2023-10-20 | 中国矿业大学(北京) | Mining subsidence area treatment effect and comprehensive evaluation method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609782B (en) * | 2012-01-10 | 2016-04-13 | 中国矿业大学(北京) | Bian Caibian based on coal-mining subsidence simulation and second ploughing rate answers the method for optimizing on opportunity |
| CN104698507B (en) * | 2015-04-02 | 2017-05-10 | 淮南矿业(集团)有限责任公司 | Quantitative method for water resource effect in coal mining subsidence area |
| CN106021947A (en) * | 2016-05-25 | 2016-10-12 | 安徽农业大学 | Quantitative calculation method for flood storage and detention effect of coal mining subsidence area |
| CN109272250B (en) * | 2018-10-11 | 2020-05-26 | 中国水利水电科学研究院 | Flood storage and waterlogging removal evaluation method for coal mining subsidence area |
-
2020
- 2020-07-03 CN CN202010637174.4A patent/CN111798052B/en active Active
Non-Patent Citations (3)
| Title |
|---|
| A shaft pillar mining subsidence calculation using both probability integral method and numerical simulation;Peixian Li 等;《Computer Modeling in Engineering & Sciences》;20181130;第117卷(第2期);第231-249页 * |
| 高潜水位多煤层开采地表沉陷积水区动态演变预测;高旭光 等;《金属矿山》;20171031;第46卷(第10期);第39-42页 * |
| 高潜水位采煤沉陷区积水范围动态演化规律;陈晓谢 等;《煤田地质与勘探》;20200430;第48卷(第2期);第126-133页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111798052A (en) | 2020-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111798052B (en) | Dynamic prediction method for three-dimensional spatial information of coal mining subsidence ponding area of high diving space | |
| Peng et al. | Characteristics of land subsidence, earth fissures and related disaster chain effects with respect to urban hazards in Xi’an, China | |
| Olesen et al. | Neotectonic deformation in Norway and its implications: a review. | |
| Ketcheson et al. | The hydrological functioning of a constructed fen wetland watershed | |
| Shi et al. | Functional efficiency assessment of the water curtain system in an underground water-sealed oil storage cavern based on time-series monitoring data | |
| Daher et al. | Karst and artificial recharge: Theoretical and practical problems: A preliminary approach to artificial recharge assessment | |
| Xiao et al. | Dynamic subsidence simulation and topsoil removal strategy in high groundwater table and underground coal mining area: a case study in Shandong Province | |
| Li et al. | Innovation for sustainable mining: Integrated planning of underground coal mining and mine reclamation | |
| Zapico et al. | Geomorphic reclamation for reestablishment of landform stability at a watershed scale in mined sites: The Alto Tajo Natural Park, Spain | |
| Ahmed | Using lithologic modeling techniques for aquifer characterization and groundwater flow modeling of the Sohag area, Egypt | |
| Liu et al. | A boundary model of terrain reconstruction in a coal-mining subsidence waterlogged area | |
| CN103235853A (en) | Earthwork balance based dike-pond layout method used for coal mining while refilling in coal-mining subsidence | |
| Guzy et al. | Land subsidence estimation for aquifer drainage induced by underground mining | |
| Zhu et al. | A novel calculation method of subsidence waterlogging spatial information based on remote sensing techniques and surface subsidence prediction | |
| Shi et al. | Spatial variability of the groundwater exploitation potential in an arid alluvial-diluvial plain using GIS-based Dempster-Shafer theory | |
| Mylopoulos et al. | Groundwater modeling and management in a complex lake-aquifer system | |
| Wang et al. | Feasibility prediction analysis of groundwater reservoir construction based on GMS and Monte Carlo analyses: a case study from the Dadougou Coal Mine, Shanxi Province, China | |
| Cooper | Chapter 14 Geological hazards from salt mining, brine extraction and natural salt dissolution in the UK | |
| CN113742995A (en) | Mine water inflow prediction method and system based on coal mine big data | |
| Xue et al. | Determination of hydraulic conductivity and its spatial variability in the Jianghan Plain using a multi-format, multi-method approach | |
| Xu et al. | Impacts of aquitard properties on an overlying unconsolidated aquifer in a mining area of the loess plateau: case study of the Changcun Colliery, Shanxi | |
| Voudouris et al. | Groundwater recharge via deep boreholes in the Patras Industrial Area aquifer system (NW Peloponnesus, Greece) | |
| Marinin et al. | State-of-Art of mine engineering reclamation while developing of steep-dipping ore fields | |
| Donovan et al. | Mine flooding history of a regional below-drainage coalfield dominated by barrier leakage (1970–2014) | |
| Koptieva | Mining Landscapes of the Kryvyi Rih Landscape-Technical System |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |