CN116892390A - Method and device for determining mining range of downward coal of building in coal mining subsidence area - Google Patents
Method and device for determining mining range of downward coal of building in coal mining subsidence area Download PDFInfo
- Publication number
- CN116892390A CN116892390A CN202310994038.4A CN202310994038A CN116892390A CN 116892390 A CN116892390 A CN 116892390A CN 202310994038 A CN202310994038 A CN 202310994038A CN 116892390 A CN116892390 A CN 116892390A
- Authority
- CN
- China
- Prior art keywords
- mining
- building
- coal
- blocks
- determining
- 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.)
- Pending
Links
- 238000005065 mining Methods 0.000 title claims abstract description 266
- 239000003245 coal Substances 0.000 title claims abstract description 192
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000000994 depressogenic effect Effects 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 8
- 230000000717 retained effect Effects 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 16
- 238000004364 calculation method Methods 0.000 description 15
- 238000013178 mathematical model Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 239000002023 wood Substances 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 239000011150 reinforced concrete Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method and a device for determining a mining range of a building under-pressure coal in a coal mining subsidence area, wherein the method comprises the following steps: obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined; and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value. According to the method, the mining range of the coal mine is reserved maximally in a manner of deleting the blocks by determining the movement deformation data of each block to the ground surface of the building, the coal mine mining scheme is accurately and rapidly given, and the maximum mining benefit is realized on the premise of protecting the building.
Description
Technical Field
The invention relates to the technical field of coal mining, in particular to a method and a device for determining a mining range of a building under-pressure coal in a coal mining subsidence area.
Background
The goaf is formed after coal mining, when the goaf area is large enough, the overlying strata of the goaf moves, and the earth surface is deformed, and the process is also called earth surface movement. Ground surface movement can cause different levels of damage to goaf structures.
In order to minimize the loss caused by subsidence in coal mining, it is necessary to predict the movement deformation of the earth's surface that will occur during the mining process prior to coal mining, and to determine the actual coal mining range based on the prediction results.
In recent years, with the continuous perfection and maturity of mathematical models, the automatic processing of actual data of a coal mining area, the automatic mapping of subsidence predicted data and the like are realized through the mathematical models.
Although the workload of researchers on the prediction of the movement deformation of the earth surface caused by coal mining can be reduced through a mathematical model, the researchers still need to continuously adjust when determining the coal mining scheme so as to determine a relatively balanced scheme between the coal mine resource extraction and the influence of the movement deformation of the earth surface caused by coal mining on the building as much as possible, thereby meeting the requirement of maximizing the production benefit.
Disclosure of Invention
The invention provides a method and a device for determining a mining range of a building under-pressure coal in a coal mining subsidence area, which are used for solving the defect that in the prior art, the mining scheme is difficult to balance between mining coal resources and the influence of ground surface movement deformation caused by mining on the building, and realizing the method for determining the mining range of the coal mine, which is convenient for meeting the requirement of maximizing the production benefit.
The invention provides a method for determining a mining range of a building under-pressure coal in a coal mining subsidence area, which comprises the following steps:
obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined;
and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
According to the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area provided by the invention, the step of deleting the blocks in turn until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value comprises the following steps:
sequentially deleting the blocks according to the sequence from big to small of the surface movement deformation data corresponding to the blocks until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value.
According to the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area, the step of determining the ground surface movement deformation data of the building after each block is mined comprises the following steps:
And inputting the geological mining conditions and mining parameters of the blocks into a ground subsidence prediction model to obtain the ground movement deformation data output by the ground subsidence prediction model.
According to the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area provided by the invention, the step of inputting the geological mining conditions and the mining parameters of the blocks into the ground surface subsidence prediction model to obtain the ground surface movement deformation data output by the ground surface subsidence prediction model comprises the following steps:
converting the coordinates of the working surface and the coordinates of the building into a coordinate system established by taking the trend of the coal seam as the positive direction of the x axis and taking the trend of the coal seam as the positive direction of the y axis;
and inputting the coordinates converted by the working face, the coordinates converted by the building, the geological mining conditions and mining parameters of the blocks into a ground surface subsidence prediction model to obtain the ground surface movement deformation data output by the ground surface subsidence prediction model.
According to the method for determining the mining range of the building pressed coal in the coal mining subsidence area provided by the invention, before the step of deleting the blocks in sequence under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, the method further comprises the following steps:
Determining a coal mining damage level corresponding to a protection level of the building, wherein the protection level is pre-associated with the coal mining damage level;
and obtaining the preset threshold according to preset earth surface movement deformation data corresponding to the coal mining damage level, wherein the coal mining damage level is pre-associated with the preset earth surface movement deformation data.
According to the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area provided by the invention, before the step of acquiring the building in the mining subsidence influence range of the working face, the method further comprises the following steps:
determining a rectangular frame according to the maximum value and the minimum value on each coordinate axis in the angular point coordinates of all the working surfaces;
and determining the mining subsidence influence range of all working surfaces according to the rectangular frame.
The invention also provides a device for determining the mining range of the building downward pressure coal in the coal mining subsidence area, which comprises the following steps:
the determining module is used for obtaining a building in the mining subsidence influence range of the working face, dividing the working face into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined;
and the selection module sequentially deletes the blocks under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
The invention also provides electronic equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the method for determining the mining range of the building under pressure in the coal mining subsidence area when executing the program.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method of determining a range of mining of depressed coal of a building in a coal mining subsidence area as described in any one of the above.
The invention also provides a computer program product comprising a computer program which when executed by a processor implements a method of determining a range of mining of a depressed coal of a building in a coal mining subsidence area as defined in any one of the above.
According to the method and the device for determining the mining range of the downward coal of the building in the coal mining subsidence area, the preset working face is divided into the plurality of blocks, the ground surface movement deformation data of the building in the mining subsidence influence range of each block on the coal mining are calculated respectively, and the blocks are deleted in sequence until the sum of the ground surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value under the condition that the ground surface movement deformation data corresponding to the reserved blocks is larger than the preset threshold value, namely, under the condition that the building in the mining subsidence influence range is protected from being damaged due to the coal mining, the coal mining working face is reserved as much as possible, and the mining benefit is maximized.
Drawings
In order to more clearly illustrate the invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for determining the mining range of a building under-pressure coal in a coal mining subsidence area;
FIG. 2 is a schematic diagram mainly used for showing working surfaces and blocks in the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area;
FIG. 3 is a second flow chart of a method for determining a range of mining of a depressed coal of a building in a coal mining subsidence area according to the present invention;
fig. 4 is a schematic diagram of a method for determining a mining range of a depressed coal of a building in a coal mining subsidence area, which is mainly used for displaying a deleted block of a working face.
FIG. 5 is a schematic structural view of a device for determining the mining range of the depressed coal of a building in a coal mining subsidence area;
fig. 6 is a schematic structural diagram of an electronic device provided by the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The following will be described first:
coal is used as a primary energy source, which plays a great role in the development of national economy, but at the same time, the subsidence of the earth surface caused by coal exploitation is also a problem to be solved. When the goaf area caused by coal mining is large enough, the ground surface is caused to move and deform, various environment and resource problems are generated, and even corresponding geological environment disasters can be caused.
In order to reduce the loss caused by mining subsidence to the maximum extent, a correct prediction must be made on the movement deformation condition of the ground surface which will occur in the mining process before coal mining, and the prediction result will be used as an important judgment basis for whether the ground surface building needs reinforcement maintenance or relocation reconstruction, and different treatment measures correspond to different costs.
In the working face mining design related to the building coal pressing, in order to not only sufficiently extract coal resources, but also protect the surface building as much as possible, the mining scheme needs to be continuously adjusted or a corresponding surface building reinforcing method is proposed to achieve the purpose of meeting the two conditions.
In the past, most of researches on technical staff in this respect are mainly aimed at accurately acquiring surface movement deformation data, and the main difference between calculation methods is that mathematical models used are different, and the application conditions of the mathematical models are respectively emphasized by an equivalent transformation line integration algorithm, a trapezoid segmentation algorithm, a triangle segmentation method, a rectangular segmentation algorithm and the like aiming at an irregular working surface. Along with the continuous perfection and maturity of mathematical models in recent years, the preprocessing and the post-processing of optimized data become main work, the functions of automatic processing of measured data, automatic mapping of subsidence predicted data and the like save considerable workload, and the working efficiency is correspondingly improved. However, the research on the automatic design working face is less at present, and researchers always need to continuously try and miss when seeking a reasonable exploitation scheme, and a relative balance is sought between exploiting resources and meeting the ground surface movement deformation so as to meet the requirement of maximizing the production benefit.
The following describes a method for determining a mining range of building depressed coal in a coal mining subsidence area according to the present invention with reference to fig. 1 to 4, and as shown in fig. 1, the method for determining a mining range of building depressed coal in a coal mining subsidence area includes:
step 101, obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined;
the working face is an operation area during coal mining, and operators mine and transport coal on the working face. The working face is defined in advance according to the division condition of the ore to be mined.
After coal mining, the ground surface in a certain range on the periphery of the mining range is moved and deformed, and the ground surface beyond the range is less influenced.
Thus, the area where the coal mining causes the movement deformation of the ground surface is taken as the mining subsidence influence range, and the buildings located in the mining subsidence range are the buildings influenced by the coal mining.
Optionally, determining the mining subsidence influence range after coal mining based on the working face as a calculation reference is used for determining the building located in the mining subsidence influence range.
Alternatively, the working surface may be one or more.
In the case of one face, the production subsidence influence range is determined with reference to the face.
Under the condition that a plurality of working surfaces are provided, a mining influence range can be determined for each working surface, and the mining influence ranges corresponding to all the working surfaces are summarized to be used as a total mining influence range; the preset area can be determined firstly, all the working surfaces are located in the preset area, and then the total mining influence range corresponding to all the working surfaces is determined by taking the preset area as a reference, so that the missing of buildings influenced by coal mining is avoided as much as possible.
After determining the mining subsidence influence range of the working face, dividing the working face into a plurality of blocks, and determining the earth surface movement deformation influence data of each block on the building in the mining subsidence influence range.
Alternatively, the shape of the segments may be the same as the working surface or may be different from the working surface.
For example, in the case where the working surface is rectangular, the shape of the segments may be rectangular having the same aspect ratio as the working surface, or may be square.
Alternatively, the shapes of the segments may be the same or different.
Preferably, as shown in fig. 2, in order to calculate the influence data of each block on the movement deformation of the ground surface at the building, the working surface in the application is rectangular, and the blocks are square and have the same shape.
It should be noted that buildings within the influence of mining subsidence include one or more.
In the case where the number of buildings is plural, the surface movement deformation influence data of each block for each building is calculated and recorded separately.
Optionally, the surface movement deformation data includes surface subsidence data, surface inclination deformation data, and surface level deformation data, and the influence of coal mining on the surface movement deformation at the building is described by a surface subsidence value, a surface inclination deformation value, and a surface level deformation value at the building.
And 102, deleting the blocks in sequence until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to a preset threshold value under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than the preset threshold value.
And the sum of the earth surface movement deformation data corresponding to all the blocks, namely the earth surface movement deformation data of a certain building in the influence range of mining subsidence after mining the coal mine.
When the sum of the earth surface movement deformation data is larger than a preset threshold value, the building is considered to be damaged due to the earth surface movement deformation caused by coal mining, so that all the coal mines in the working surface cannot be mined.
On the basis, partial blocks in the working surface are selected for exploitation. Specifically, the blocks are sequentially deleted according to a preset sequence until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to a preset threshold value.
That is, under the condition that the building is not damaged due to the movement deformation of the ground surface, the maximum mining area is reserved as many blocks as possible, so that the balance between building protection and economic benefit maximization is achieved.
Optionally, the ground movement deformation data is less than or equal to a preset threshold, representing that the ground subsidence data, the ground inclination deformation data and the ground level deformation data at the building are all less than or equal to the corresponding preset thresholds.
Optionally, the preset sequence may be the size sequence of the surface movement deformation data corresponding to the blocks, or may be the position sequence of the blocks.
In one possible implementation, the blocks are deleted sequentially from large to small according to the surface movement deformation data corresponding to the blocks, and the blocks with the largest number are reserved so as to realize the maximization of exploitation benefits under the condition that the building is not destroyed.
In another possible embodiment, the tiles are deleted sequentially according to the location between the tile and the building. For example, a block closer to the building has a larger influence on a dip value in the ground surface movement deformation data at the building, a block closer to the building has a larger influence on inclination deformation data in the ground surface movement deformation data at the building, and therefore, a group of blocks farthest and closest to the building are sequentially deleted until the ground surface movement deformation data corresponding to the remaining blocks is smaller than a preset threshold. The method can remove the blocks to concentrate the reserved block positions as much as possible, is convenient for subsequent exploitation work, reduces exploitation cost and maximizes exploitation benefits.
Optionally, the preset threshold is determined according to building structure and/or building type.
In one possible embodiment, the preset threshold is determined according to the type of building. For example, when the building is a national precious cultural building or an ultra-high-rise building, the protection level is higher and the preset threshold is lower; when the building is a common wood structure bearing house, the protection level is lower, and the preset threshold value is higher.
In another possible embodiment, the preset threshold is determined according to the structure of the building. For example, when the building is of a reinforced concrete structure, the degree of deformation of the ground surface that the building can bear is high, and the corresponding preset threshold value is high; when the building is of a common wood structure, the deformation degree of the ground surface which can be born by the building is low, and the corresponding preset threshold value is low.
In addition, the supporting structure of the building can influence the deformation degree of the ground surface which can be born by the building, and a corresponding preset threshold value can be determined for each building in the mining subsidence influence range according to the actual category and structure of the building.
Optionally, the number of working surfaces is determined according to the actual distribution condition of the coal mine.
On the basis, the application does not limit the number of the blocks, the more the number of the blocks is, the smaller the movement deformation data of each block corresponds to the surface, and the larger the area of the working surface reserved after the blocks are deleted is; the smaller the number of the blocks is, the smaller the surface movement deformation data is corresponding to each block, and the simpler the surface movement deformation data calculation process of each block is.
According to the invention, the preset working face is divided into a plurality of blocks during coal mining, the ground surface movement deformation data of the building in the mining subsidence influence range of each block on the coal mining are calculated respectively, and the blocks are deleted in sequence until the sum of the ground surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value under the condition that the ground surface movement deformation data corresponding to the reserved blocks is larger than the preset threshold value, namely, under the condition that the building is not damaged due to the coal mining in the mining subsidence influence range, the coal mining working face is reserved as much as possible, and the mining benefit is maximized.
In the method for determining the mining range of the building pressed coal in the coal mining subsidence area, the steps of deleting the blocks in turn until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value comprise the following steps:
sequentially deleting the blocks according to the sequence from big to small of the surface movement deformation data corresponding to the blocks until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value.
Specifically, the blocks are deleted in sequence from large to small according to the earth surface movement deformation data corresponding to the blocks, so that the blocks can be reserved as much as possible, and the reserved working surface range is maximized.
It should be noted that, since the surface movement deformation data includes the surface subsidence data, the surface inclination deformation data, and the surface level deformation data, the surface movement deformation data needs to be preprocessed at the time of sorting.
In a possible implementation manner, preprocessing includes sorting the subsurface subsidence data, the surface inclination deformation data and the surface horizontal deformation data corresponding to each block respectively to generate three arrays, and deleting each array in sequence from large to small until the sum of the subsurface subsidence data, the sum of the surface inclination deformation data and the sum of the surface horizontal deformation data corresponding to the reserved blocks do not exceed corresponding preset thresholds.
In another possible embodiment, the preprocessing includes performing weighted average processing on the surface subsidence data, the surface inclination deformation data and the surface level deformation data of each block, and summing the weighted average processing and the surface subsidence deformation data as comprehensive surface level deformation data of each block, and deleting the blocks according to the sequence from the large to the small of the comprehensive surface level deformation data until the comprehensive surface level deformation data corresponding to the reserved blocks is smaller than a corresponding preset threshold value.
According to the invention, the blocks are deleted in sequence according to the sequence of the ground surface movement deformation data from large to small, so that the range of the working surface corresponding to the reserved blocks is maximized, more coal mine resources can be obtained through exploitation, and the exploitation benefit is maximized.
In the method for determining the mining range of the depressed coal of the building in the coal mining subsidence area, the step of determining the movement deformation data of the ground surface of the building after each block is mined comprises the following steps:
and inputting the geological mining conditions and mining parameters of the blocks into a ground subsidence prediction model to obtain the ground movement deformation data output by the ground subsidence prediction model.
Geological mining conditions comprise mining height, maximum mining depth of a working face, inclination angle of a coal bed and inclination tendency of the coal bed, which correspond to the blocks; the mining parameters comprise a sinking coefficient, a main impact angle tangent, a horizontal movement coefficient, a mining impact propagation angle and an inflection point offset distance corresponding to the working surface.
Optionally, the address mining conditions are obtained in advance by geological exploration.
Alternatively, the mining parameters are determined in accordance with a predetermined coal mining scheme.
After the geological mining conditions and mining parameters of the blocks are input into the ground subsidence prediction model, the model outputs ground movement deformation data generated by each block at the building according to an equivalent transformation line integration algorithm.
Optionally, each of the geological mining conditions and mining parameters is stored in an array for later recall and analysis.
Optionally, each block is numbered, and the surface movement deformation data output by the model is associated with each block through the number.
Further, under the condition that the surface recoverable range corresponding to the reserved block is not consistent or the recoverable range is too small, the recovery mode can be modified, for example, longwall recovery is changed into limited-thickness recovery or strip recovery, recovery parameters are determined again and are input into the surface subsidence prediction model, and the reserved block and the recoverable range corresponding to the block are determined until the recoverable range of the surface corresponding to the reserved block achieves the unification of the economic benefit of coal mining and the cost and difficulty of coal mining.
In the method for determining the mining range of the building depressed coal in the coal mining subsidence area, the steps of inputting the geological mining conditions and mining parameters of the blocks into a ground surface subsidence prediction model to obtain the ground surface movement deformation data output by the ground surface subsidence prediction model comprise the following steps:
converting the coordinates of the working surface and the coordinates of the building into a coordinate system established by taking the trend of the coal seam as the positive direction of the x axis and taking the trend of the coal seam as the positive direction of the y axis;
the trend of the coal bed refers to the extending direction of the coal bed on the ground surface, the trend of the coal bed refers to the included angle between the coal bed and the horizontal plane, and the trend of the coal bed are mutually perpendicular.
When geological mining conditions are obtained through exploration, the trend of the coal bed is taken as the positive direction of the x axis, and the coordinate system established by taking the trend of the coal bed as the positive direction of the y axis is recorded to obtain relevant data.
Therefore, in order to facilitate calculation, the coordinates of the working surface and the coordinates of the building are converted from a planar coordinate system to a coordinate system established by taking the trend of the coal seam as the positive direction of the x axis and taking the trend of the coal seam as the positive direction of the y axis, so that the coordinate system used by various data of the input model is unified, and the complex calculation process brought by the trend value parameters is eliminated.
The coordinates of the working face comprise the coordinates of the corner points of the working face and the coordinates of the corner points of each block.
And inputting the coordinates converted by the working face, the coordinates converted by the building, the geological mining conditions and mining parameters of the blocks into a ground surface subsidence prediction model to obtain the ground surface movement deformation data output by the ground surface subsidence prediction model.
And inputting the coordinates converted by the working face, the coordinates converted by the building, the geological mining conditions of the blocks and the mining parameters into the ground surface subsidence prediction model to obtain ground surface movement deformation data output by the ground surface subsidence prediction model, namely ground surface movement deformation data of each block on the building.
It should be noted that, in the present application, when calculating the surface movement deformation data, for each block, the surface movement deformation data of the coal seam of 1 meter thick per mining on the building is calculated and used as the surface movement deformation data corresponding to the block.
Optionally, as shown in fig. 3, the calculation of the surface movement deformation data of each block and the deletion of the block are implemented based on a GIS system (Geographic Information System ). Specifically, through integration of GIS and subsidence prediction mathematical model, according to the relevant analysis and prediction data, the GIS is utilized to carry out space visual display on the calculation result; and the GIS component is utilized to analyze the problem of the ground surface movement deformation rule.
The preparation module is used for data entry, and specifically comprises the number of working faces, coordinates of all corner points of the working faces, mining parameters and geological mining condition data, and all the data are stored in an array mode.
The surface influence range calculation module is used for determining the mining subsidence influence range corresponding to the working surface, and further determining the building in the mining subsidence and the coordinates corresponding to the building.
The working face segmentation module is used for segmenting the working face into a plurality of segments and determining the angular position of each segment.
The coordinate conversion module is used for converting the working face coordinates and the building coordinates.
And in the deformation calculation module, the input data and the converted coordinates are calculated through the deformation calculation module, the ground surface subsidence prediction model in the deformation calculation module outputs ground surface movement deformation data of each layer which is 1m thick per each block and is caused to the building, and the output data and each block are in corresponding relation.
The sorting module is used for sorting the blocks according to the size of the surface movement deformation data, deleting the blocks according to the sequence from large to small under the condition that the sum of the surface movement deformation data corresponding to the blocks is larger than a preset threshold value until the sum of the surface movement deformation data corresponding to the rest blocks is not larger than the preset threshold value, and reserving all the rest blocks as the actual acquirable range of the coal mine.
It should be noted that, because the present application uses the coordinate system corresponding to the trend and trend of the coal seam in the calculation, after the calculation is completed, the corner coordinates of all the working surfaces, the corner coordinates of the blocks and the corner coordinates of the building need to be converted back into the plane coordinate system.
And the output module is used for generating a surface movement deformation data table at each building and drawing a contour map, and the contour map is shown in fig. 4.
Specifically, the large positive direction in fig. 4 represents two working surfaces, the circle represents a building within the influence range of mining subsidence, each small square represents a block, wherein the working surface does not contain the position of the small square, namely the deleted block, and the rest position is the actual mining range.
And the mining of the coal mine is completed within the actual mining range, so that the corresponding building is not damaged.
Alternatively, as shown in fig. 4, after deleting the tiles, the remaining working surface is mostly an irregular polygon, and the working surface is generally rectangular. Thus, the corresponding working surface at the time of actual mining can be determined again on the basis of the reserved blocks so as to facilitate the development of the coal mining work.
Optionally, the area of the corresponding working surface is smaller than the reserved block during actual mining, so that coal mines at positions outside the reserved block are mismined, and the building is damaged.
In the method for determining the mining range of the building under-pressure coal in the coal mining subsidence area, before the step of deleting the blocks in sequence under the condition that the sum of the earth surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, the method further comprises the following steps:
determining a coal mining damage level corresponding to a protection level of the building, wherein the protection level is pre-associated with the coal mining damage level;
Specifically, the protection level of the building refers to the twenty-first item of "building, water body, railway and main roadway coal pillar set and pressed coal exploitation Specification", and the protection level is specifically shown in the following table 1:
TABLE 1
Further, the mining damage degree is divided into four grades according to the twentieth item of the mining standards of building, water body, railway and main roadway coal pillar reservation and coal pressing, wherein the four grades of inclined deformation are respectively grade I0-3 mm/m, grade II 3-6 mm/m, grade III 6-10 mm/m and grade V more than 10mm/m; the four grades of horizontal deformation are respectively grade I0-2 mm/m, grade II 2-4 mm/m, grade III 4-6 mm/m and grade V greater than 6mm/m.
Alternatively, the grades corresponding to the dip value data are empirically classified into four grades.
Optionally, the protection level is empirically pre-associated with the coal mining damage level.
On this basis, each building of the mining subsidence influence range class has its corresponding coal mining damage level.
And obtaining the preset threshold according to preset earth surface movement deformation data corresponding to the coal mining damage level, wherein the coal mining damage level is pre-associated with the preset earth surface movement deformation data.
Further, according to the preset ground surface movement deformation data corresponding to the coal mining damage level, namely the inclination deformation value, the horizontal deformation value and the sinking value corresponding to each level in the specification, a preset threshold value is obtained.
Optionally, since the surface movement deformation data corresponding to each level is a data range, in a feasible embodiment, an intermediate value in the data range is selected as the preset threshold corresponding to the level; in another possible implementation, the minimum value in the data range is selected as the preset threshold corresponding to the level, so as to ensure the safety of the building.
The method for determining the exploitation range of the depressed coal of the building in the coal mining subsidence area, before the step of acquiring the building in the exploitation subsidence influence range of the working face, further comprises the following steps:
determining a rectangular frame according to the maximum value and the minimum value on each coordinate axis in the angular point coordinates of all the working surfaces;
the working surface includes one or more.
Under the condition that a plurality of working surfaces are provided, obtaining the maximum value and the minimum value on each coordinate axis in the angular point coordinates of all the working surfaces, namely x min 、x max 、y min And y max A rectangular box is determined. The coordinate system at this time is a plane coordinate system.
Wherein the coordinates of the four corner points of the rectangular frame are (x) min ,y min )、(x min ,y max )、(x max ,y min ) And (x) max ,y max )。
The rectangular frame determined in the mode can cover all working surfaces, so that the mining subsidence influence range determined based on the rectangular frame is maximum, and a building which is possibly damaged is not missed.
And determining the mining subsidence influence range of all working surfaces according to the rectangular frame.
Optionally, a range within a preset mining influence radius from the rectangular box boundary is taken as the mining subsidence influence range.
Wherein the preset mining influence radius is empirically determined.
The following describes a device for determining a range of exploitation of building under-pressure coal in a coal mining subsidence area, and the device for determining the range of exploitation of building under-pressure coal in the coal mining subsidence area and the method for determining the range of exploitation of building under-pressure coal in the coal mining subsidence area described above can be correspondingly referred to each other.
As shown in fig. 5, the device for determining the mining range of the depressed coal of the building in the coal mining subsidence area comprises a determining module 501 and a selecting module 502:
a determining module 501, configured to obtain a building within a mining subsidence influence range of a working surface, divide the working surface into a plurality of blocks, and determine movement deformation data of each block on a ground surface of the building after coal mining;
The working face is an operation area during coal mining, and operators mine and transport coal on the working face. The working face is defined in advance according to the division condition of the ore to be mined.
After coal mining, the ground surface in a certain range on the periphery of the mining range is moved and deformed, and the ground surface beyond the range is less influenced.
Thus, the area where the coal mining causes the movement deformation of the ground surface is taken as the mining subsidence influence range, and the buildings located in the mining subsidence range are the buildings influenced by the coal mining.
Optionally, determining the mining subsidence influence range after coal mining based on the working face as a calculation reference is used for determining the building located in the mining subsidence influence range.
Alternatively, the working surface may be one or more.
In the case of one face, the production subsidence influence range is determined with reference to the face.
Under the condition that a plurality of working surfaces are provided, a mining influence range can be determined for each working surface, and the mining influence ranges corresponding to all the working surfaces are summarized to be used as a total mining influence range; the preset area can be determined firstly, all the working surfaces are located in the preset area, and then the total mining influence range corresponding to all the working surfaces is determined by taking the preset area as a reference, so that the missing of buildings influenced by coal mining is avoided as much as possible.
After determining the mining subsidence influence range of the working face, dividing the working face into a plurality of blocks, and determining the earth surface movement deformation influence data of each block on the building in the mining subsidence influence range.
Alternatively, the shape of the segments may be the same as the working surface or may be different from the working surface.
For example, in the case where the working surface is rectangular, the shape of the segments may be rectangular having the same aspect ratio as the working surface, or may be square.
Alternatively, the shapes of the segments may be the same or different.
Preferably, as shown in fig. 2, in order to calculate the influence data of each block on the movement deformation of the ground surface at the building, the working surface in the application is rectangular, and the blocks are square and have the same shape.
It should be noted that buildings within the influence of mining subsidence include one or more.
In the case where the number of buildings is plural, the surface movement deformation influence data of each block for each building is calculated and recorded separately.
Optionally, the surface movement deformation data includes surface subsidence data, surface inclination deformation data, and surface level deformation data, and the influence of coal mining on the surface movement deformation at the building is described by a surface subsidence value, a surface inclination deformation value, and a surface level deformation value at the building.
And the selection module 502 is configured to delete the blocks in sequence when the sum of all the surface movement deformation data corresponding to the blocks is greater than a preset threshold value, until the sum of the surface movement deformation data corresponding to the reserved blocks is less than or equal to the preset threshold value.
And the sum of the earth surface movement deformation data corresponding to all the blocks, namely the earth surface movement deformation data of a certain building in the influence range of mining subsidence after mining the coal mine.
When the sum of the earth surface movement deformation data is larger than a preset threshold value, the building is considered to be damaged due to the earth surface movement deformation caused by coal mining, so that all the coal mines in the working surface cannot be mined.
On the basis, partial blocks in the working surface are selected for exploitation. Specifically, the blocks are sequentially deleted according to a preset sequence until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to a preset threshold value.
That is, under the condition that the building is not damaged due to the movement deformation of the ground surface, the maximum mining area is reserved as many blocks as possible, so that the balance between building protection and economic benefit maximization is achieved.
Optionally, the ground movement deformation data is less than or equal to a preset threshold, representing that the ground subsidence data, the ground inclination deformation data and the ground level deformation data at the building are all less than or equal to the corresponding preset thresholds.
Optionally, the preset sequence may be the size sequence of the surface movement deformation data corresponding to the blocks, or may be the position sequence of the blocks.
In one possible implementation, the blocks are deleted sequentially from large to small according to the surface movement deformation data corresponding to the blocks, and the blocks with the largest number are reserved so as to realize the maximization of exploitation benefits under the condition that the building is not destroyed.
In another possible embodiment, the tiles are deleted sequentially according to the location between the tile and the building. For example, a block closer to the building has a larger influence on a dip value in the ground surface movement deformation data at the building, a block closer to the building has a larger influence on inclination deformation data in the ground surface movement deformation data at the building, and therefore, a group of blocks farthest and closest to the building are sequentially deleted until the ground surface movement deformation data corresponding to the remaining blocks is smaller than a preset threshold. The method can remove the blocks to concentrate the reserved block positions as much as possible, is convenient for subsequent exploitation work, reduces exploitation cost and maximizes exploitation benefits.
Optionally, the preset threshold is determined according to building structure and/or building type.
In one possible embodiment, the preset threshold is determined according to the type of building. For example, when the building is a national precious cultural building or an ultra-high-rise building, the protection level is higher and the preset threshold is lower; when the building is a common wood structure bearing house, the protection level is lower, and the preset threshold value is higher.
In another possible embodiment, the preset threshold is determined according to the structure of the building. For example, when the building is of a reinforced concrete structure, the degree of deformation of the ground surface that the building can bear is high, and the corresponding preset threshold value is high; when the building is of a common wood structure, the deformation degree of the ground surface which can be born by the building is low, and the corresponding preset threshold value is low.
In addition, the supporting structure of the building can influence the deformation degree of the ground surface which can be born by the building, and a corresponding preset threshold value can be determined for each building in the mining subsidence influence range according to the actual category and structure of the building.
Optionally, the number of working surfaces is determined according to the actual distribution condition of the coal mine.
On the basis, the application does not limit the number of the blocks, the more the number of the blocks is, the smaller the movement deformation data of each block corresponds to the surface, and the larger the area of the working surface reserved after the blocks are deleted is; the smaller the number of the blocks is, the smaller the surface movement deformation data is corresponding to each block, and the simpler the surface movement deformation data calculation process of each block is.
According to the invention, the preset working face is divided into a plurality of blocks during coal mining, the ground surface movement deformation data of the building in the mining subsidence influence range of each block on the coal mining are calculated respectively, and the blocks are deleted in sequence until the sum of the ground surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value under the condition that the ground surface movement deformation data corresponding to the reserved blocks is larger than the preset threshold value, namely, under the condition that the building is not damaged due to the coal mining in the mining subsidence influence range, the coal mining working face is reserved as much as possible, and the mining benefit is maximized.
Fig. 6 illustrates a physical schematic diagram of an electronic device, as shown in fig. 6, which may include: processor 610, communication interface (Communications Interface) 620, memory 630, and communication bus 640, wherein processor 610, communication interface 620, and memory 630 communicate with each other via communication bus 640. Processor 610 may invoke logic instructions in memory 630 to perform a method of determining a range of mining of a depressed coal of a building within a coal mining subsidence area, the method comprising: obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined; and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
Further, the logic instructions in the memory 630 may be implemented in the form of software functional units and stored in a computer-readable storage medium when sold or used as a stand-alone product. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program, the computer program being storable on a non-transitory computer readable storage medium, the computer program when executed by a processor being capable of performing the method of determining a range of mining of depressed coal of a building in a coal mining subsidence area provided by the methods described above, the method comprising: a method for determining a range of mining of a depressed coal of a building in a coal mining subsidence area, the method comprising: obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined; and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
In yet another aspect, the present invention provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, is implemented to perform the method for determining a range of mining of depressed coal of a building in a coal mining subsidence area provided by the above methods, the method comprising: a method for determining a range of mining of a depressed coal of a building in a coal mining subsidence area, the method comprising: obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined; and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
The apparatus embodiments described above are merely illustrative, wherein the elements illustrated as separate elements may or may not be physically separate, and the elements shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art will understand and implement the present invention without undue burden.
From the above description of the embodiments, it will be apparent to those skilled in the art that the embodiments may be implemented by means of software plus necessary general hardware platforms, or of course may be implemented by means of hardware. Based on this understanding, the foregoing technical solution may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a computer readable storage medium, such as ROM/RAM, a magnetic disk, an optical disk, etc., including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method described in the respective embodiments or some parts of the embodiments.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for determining a range of mining of a depressed coal of a building in a coal mining subsidence area, comprising:
obtaining a building in a mining subsidence influence range of a working surface, dividing the working surface into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined;
and deleting the blocks in turn under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than a preset threshold value, until the reserved sum of the surface movement deformation data corresponding to the blocks is smaller than or equal to the preset threshold value.
2. The method for determining the mining range of the depressed coal of the building in the coal mining subsidence area according to claim 1, wherein the step of sequentially deleting the segments until the sum of the surface movement deformation data corresponding to the retained segments is less than or equal to the preset threshold value comprises:
sequentially deleting the blocks according to the sequence from big to small of the surface movement deformation data corresponding to the blocks until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value.
3. The method of determining a range of mining of a depressed coal of a building in a coal mining subsidence area according to claim 1, wherein said step of determining movement deformation data for a ground surface at said building after each of said segmented coals comprises:
And inputting the geological mining conditions and mining parameters of the blocks into a ground subsidence prediction model to obtain the ground movement deformation data output by the ground subsidence prediction model.
4. A method for determining a range of mining of building depressed coals in a coal mining subsidence area as set forth in claim 3, wherein said step of inputting said segmented geological mining conditions and mining parameters into a ground surface subsidence prediction model to obtain said ground surface movement deformation data outputted from said ground surface subsidence prediction model comprises:
converting the coordinates of the working surface and the coordinates of the building into a coordinate system established by taking the trend of the coal seam as the positive direction of the x axis and taking the trend of the coal seam as the positive direction of the y axis;
and inputting the coordinates converted by the working face, the coordinates converted by the building, the geological mining conditions and mining parameters of the blocks into a ground surface subsidence prediction model to obtain the ground surface movement deformation data output by the ground surface subsidence prediction model.
5. The method for determining a range of mining under-pressure coal in a building in a coal mining subsidence area according to any one of claims 1 to 4, further comprising, before the step of sequentially deleting the segments, if a sum of the surface movement deformation data corresponding to all the segments is greater than a preset threshold:
Determining a coal mining damage level corresponding to a protection level of the building, wherein the protection level is pre-associated with the coal mining damage level;
and obtaining the preset threshold according to preset earth surface movement deformation data corresponding to the coal mining damage level, wherein the coal mining damage level is pre-associated with the preset earth surface movement deformation data.
6. The method for determining the range of mining for a depressed coal of a building in a coal mining subsidence area according to any one of claims 1 to 4, further comprising, prior to said step of acquiring a building within the range of influence of mining subsidence of the face:
determining a rectangular frame according to the maximum value and the minimum value on each coordinate axis in the angular point coordinates of all the working surfaces;
and determining the mining subsidence influence range of all working surfaces according to the rectangular frame.
7. A device for determining a range of mining of depressed coal of a building in a coal mining subsidence area, comprising:
the determining module is used for obtaining a building in the mining subsidence influence range of the working face, dividing the working face into a plurality of blocks, and determining movement deformation data of the ground surface of the building after each block is mined;
And the selection module is used for deleting the blocks in sequence until the sum of the surface movement deformation data corresponding to the reserved blocks is smaller than or equal to the preset threshold value under the condition that the sum of the surface movement deformation data corresponding to all the blocks is larger than the preset threshold value.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor, when executing the program, implements the method of determining a range of mining under building coal in a coal mining subsidence area as claimed in any one of claims 1 to 6.
9. A non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method of determining a range of mining of building downcoal in a coal mining subsidence area as claimed in any one of claims 1 to 6.
10. A computer program product comprising a computer program which, when executed by a processor, implements a method of determining a range of mining of building downcoal in a coal mining subsidence area as claimed in any one of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310994038.4A CN116892390A (en) | 2023-08-08 | 2023-08-08 | Method and device for determining mining range of downward coal of building in coal mining subsidence area |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310994038.4A CN116892390A (en) | 2023-08-08 | 2023-08-08 | Method and device for determining mining range of downward coal of building in coal mining subsidence area |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116892390A true CN116892390A (en) | 2023-10-17 |
Family
ID=88312154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310994038.4A Pending CN116892390A (en) | 2023-08-08 | 2023-08-08 | Method and device for determining mining range of downward coal of building in coal mining subsidence area |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116892390A (en) |
-
2023
- 2023-08-08 CN CN202310994038.4A patent/CN116892390A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114663627B (en) | Mine digital model building method based on three-dimensional point cloud database | |
CN107358654B (en) | Profile reconstruction three-dimensional surface modeling method and system based on polygonal deformation technology | |
CN109993832B (en) | Method for constructing refined three-dimensional model of fully mechanized coal mining face coal seam | |
CN106372297B (en) | Method for determining safe vertical distance between shield and karst cave in sandy soil karst stratum | |
CN1873674A (en) | Estimating mining subsidence under turnout cluster of dedicated railroad, and governing expert system | |
CN113326756A (en) | Method for identifying potential landslide hazard of reservoir bank based on rock mass degradation characteristics | |
CN104166163A (en) | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube | |
Li et al. | Fuzzy probability measures (FPM) based non-symmetric membership function: Engineering examples of ground subsidence due to underground mining | |
CN112443308A (en) | Sandstone uranium mine ground leaching production well network setting method and device and terminal equipment | |
CN110672072A (en) | Information multivariate interpretation and characterization method for fractured rock mass structure | |
CN110555559A (en) | Ground settlement prediction method under insufficient mining condition of thick unconsolidated formation | |
CN112184902A (en) | Underground mining face inversion method for boundary crossing mining identification | |
CN104881724A (en) | GIS-based dynamic prediction method for coal mine gas emission quantity | |
Strzalkowski et al. | Analytical and numerical method assessing the risk of sinkholes formation in mining areas | |
CN112949206A (en) | Anti-aging reinforcing system for building construction foundation | |
CN116892390A (en) | Method and device for determining mining range of downward coal of building in coal mining subsidence area | |
CN117197377A (en) | Three-dimensional dynamic visual modeling method for mining area | |
CN104899392A (en) | GIS-based intelligent analysis method for coal mine gas emission overrun prediction | |
CN110599595A (en) | Three-dimensional characterization method for geological structure | |
CN108062447B (en) | Method for predicting influence of diving precipitation on building settlement | |
DAI et al. | Water distribution extracted from mining subsidence area using Kriging interpolation algorithm | |
CN114676468A (en) | Metal strip mine final boundary optimization method and system | |
Esaki et al. | A GIS-based prediction method to evaluate subsidence-induced damage from coal mining beneath a reservoir, Kyushu, Japan | |
Shi et al. | Study on numerical models in predicting surface deformation caused by underground coal mining | |
CN110700285B (en) | Three-dimensional design method for relieving slope of surface unconsolidated formation of strip mine |
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 |