Mining parameter optimization method, mining method, storage medium and computer device

Abstract

The present invention provides a mining parameter optimization method, a mining method, a storage medium and a computer device. The mining parameter optimization method includes: taking the ratio of the volume of the damaged rock mass to the volume of the overburden formation as the overburden damage ratio of the overburden formation; for each mining parameter, according to the corresponding relationship between the overburden damage ratio and different parameter values ​​of the mining parameter, establishing a relationship model about the relationship between the overburden damage ratio of the overburden formation and the mining parameter; taking the parameter value of the mining parameter corresponding to the maximum value of the change speed as the endpoint value of the parameter value range after the mining parameter is optimized, and then determining the endpoint value of the parameter value range after each mining parameter in the multiple mining parameters is optimized. It can not only ensure the mining efficiency and realize the maximum loss reduction mining, but also control the mining damage from the source.

Description

Technical Field

[0001] The present invention relates to the field of mineral mining parameter design, and in particular to a mining parameter optimization method, a mining method, a storage medium and a computer device. Background Art

[0002] With the advancement of technology and equipment, the scale of my country's underground coal mines has grown from millions of tons (1-3 million tons / year) in the 1990s to tens of millions of tons (10-30 million tons / year) at the beginning of this century. The level of intensive mining has been continuously improved, gradually becoming an important way to improve the safety level of coal mines, resource recovery rate and economic benefits. According to the intensive development layout, the output of my country's 14 large coal bases accounts for more than 95% of the country's total.

[0003] High-intensity mining is an important way of intensive mining, but it also brings about the problems of rapid surface subsidence and large deformation, which in turn leads to serious ecological damage. The distribution of coal endowment and ecological environment capacity in my country is uneven. Among them, the coal production in the western Shanxi, Shaanxi, Inner Mongolia, Ningxia and Gansu represented by the "Shendong" mining area accounts for 3 / 4 of the country, but the ecology in this region is fragile and water resources are scarce. High-intensity mining has aggravated desertification and desertification in the mining area and its surrounding areas. The eastern mining area represented by the "Lianghuai" mining area guarantees the energy supply in East China, but high-intensity mining causes land subsidence, making land resources in mining cities more tense. The damage to the surface ecology in the mining area is caused by the deformation and breakage of the overburden caused by coal mining, the mining damage of the lower rock layer caused by the movement of the rock layer, and the resulting transmission damage of the upper rock layer and the surface layer. Therefore, the fundamental of ecological restoration in mining areas lies in minimizing the damage to the surface ecology from the source of mining.

[0004] In summary, how to coordinate efficient restoration of the surface ecology with high-intensity mining is a major technical challenge in the coal industry. However, due to the rapid advancement of high-intensity mining, the large size of the mining face, the large range of damage to the overlying rock strata, the rapid surface settlement speed and large deformation, the high degree of damage to the surface and ecology, and the prominent contradiction between resource development and ecological mining area construction, it is difficult to control the degree of surface damage from the source. For this reason, it is urgent to optimize the mining parameters of the working face from the source of mining, so as to achieve coordination between high-intensity coal mining and surface ecological protection. Summary of the invention

[0005] The main purpose of the present invention is to provide a mining parameter optimization method, a mining method, a storage medium and a computer device to optimize high-intensity mining parameters and achieve coordination between high-intensity mining and surface ecological protection.

[0006] In a first aspect, the present application provides a mining parameter optimization method, comprising the following steps: obtaining the volume of an overburden stratum, and obtaining the volume of damaged rock mass in the overburden stratum caused by mining a mineral deposit beneath the overburden stratum according to a combination of different parameter values ​​of multiple mining parameters, and taking the ratio of the volume of the damaged rock mass to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum; for each mining parameter, establishing a relationship model regarding the relationship between the overburden damage ratio of the overburden stratum and the mining parameter based on the correspondence between the overburden damage ratio and different parameter values ​​of the mining parameter; determining the maximum value of the rate of change of the overburden damage ratio based on the relationship model, taking the parameter value of the mining parameter corresponding to the maximum rate of change as the endpoint value of the parameter value range after the mining parameter is optimized, and then determining the endpoint value of the parameter value range after the optimization for each of the multiple mining parameters.

[0007] In one embodiment, the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposits under the overburden stratum according to different parameter value combinations of multiple mining parameters is obtained, including: dividing the overburden stratum into a number of unit bodies, and obtaining the actual shear stress borne by each unit body; for each unit body, determining the strain tensor of the elastic matrix in the unit body through strain analysis, and determining the theoretical shear stress that the unit body can withstand according to the strain tensor; comparing the actual shear stress borne by the unit body with the theoretical shear stress that the unit body can withstand, and determining whether the unit body is damaged according to the comparison result; obtaining the volume of each damaged unit body, and summing the volumes of each damaged unit body, and using the sum as the volume of rock mass damaged in the overburden stratum.

[0008] In one embodiment, the theoretical shear stress that the unit body can withstand is determined based on the strain tensor, including: using a preset damage constitutive model to determine the theoretical normal stress that the unit body can withstand based on the strain tensor of the elastic matrix in the unit body; using the Moore-Coulomb criterion to determine the theoretical shear stress that the unit body can withstand based on the theoretical normal stress that each unit body can withstand.

[0009] In one embodiment, the damage constitutive model is:

[0010]

[0011] Among them, σ ij is the theoretical normal stress that the unit can withstand, G is the shear modulus of the overlying stratum, G 0 is the shear modulus of the elastic matrix in the unit cell, E ijkl is the elastic constant tensor of the elastic matrix in the unit cell, ε kl is the strain tensor of the elastic matrix in the unit cell, Emmkl is the elastic modulus of the elastic matrix in the unit cell, δ ij is the Kronecker symbol;

[0012] The Mohr-Coulomb criterion is:

[0013] τ=c+σ n tanφ

[0014] Wherein, c is the cohesion between each unit body, Φ is the internal friction angle between each unit body, σ n and τ are respectively the theoretical normal stress and theoretical shear stress that the unit cell can withstand.

[0015] In one embodiment, a relationship model between the overburden damage ratio of the overburden stratum and the mining parameter is established, including: using a Sigmoid function fitting method to establish a relationship model between the overburden damage ratio of the overburden stratum and the mining parameter.

[0016] In one embodiment, the multiple mining parameters include: working face mining height, working face length and working face advancement speed; determining the maximum value of the change speed of the overburden damage ratio according to the relationship model, and using the parameter value of the mining parameter corresponding to the maximum value of the change speed as the endpoint value of the parameter value range after the mining parameter is optimized, including: when the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden stratum and the working face mining height, determining the maximum value of the increase speed of the overburden damage ratio according to the relationship model, and using the parameter value of the working face mining height corresponding to the maximum value of the increase speed as the upper limit value of the parameter value range after the working face mining height is optimized; when When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face length, the maximum value of the increase rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face length corresponding to the maximum value of the increase rate is used as the upper limit value of the parameter value range after the working face length is optimized; when the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face advancement speed, the maximum value of the decrease rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face advancement speed corresponding to the maximum value of the decrease rate is used as the lower limit value of the parameter value range after the working face advancement speed is optimized.

[0017] In one embodiment, the relationship model between the overburden damage ratio of the overburden strata and the mining height of the working face is:

[0018]

[0019] Among them, D frepresents the overburden damage ratio of the overburden strata, M represents the mining height of the working face, M1 represents the parameter value of the mining height of the working face corresponding to the maximum value of the increase rate of the overburden damage ratio of the overburden strata, and a, b and c are all fitting constants;

[0020] The relationship model between the overburden damage ratio of the overburden strata and the length of the working face is:

[0021]

[0022] Among them, D f represents the overburden damage ratio of the overburden strata, L represents the length of the working face, L1 represents the parameter value of the working face length corresponding to the maximum value of the increase rate of the overburden damage ratio of the overburden strata, and a1, b1 and c1 are all fitting constants;

[0023] The relationship model between the overburden damage ratio of the overburden strata and the working face advancement speed is:

[0024]

[0025] Among them, D f represents the overburden damage ratio of the overburden strata, v represents the working face advancement speed, v1 represents the parameter value of the working face advancement speed corresponding to the maximum value of the reduction speed of the overburden damage ratio of the overburden strata, and a2, b2 and c2 are all fitting constants.

[0026] In one embodiment, before obtaining the volume of the overburden stratum and obtaining the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposits beneath the overburden stratum according to different parameter value combinations of multiple mining parameters, the method also includes the step of: using a finite element analysis method to construct a virtual model of the overburden stratum after mining the mineral deposits beneath the overburden stratum according to different parameter value combinations of multiple mining parameters.

[0027] In a second aspect, the present application provides a mining method, comprising mining a mineral deposit beneath an overburden stratum according to mining parameters optimized by the mining parameter optimization method as described above.

[0028] In a third aspect, the present application provides a storage medium storing a computer program, which, when executed by a processor, implements the steps of the mining parameter optimization method as described above.

[0029] In a fourth aspect, the present application provides a computer device, comprising a processor and a storage medium storing program code, wherein when the program code is executed by the processor, the steps of the mining parameter optimization method described above are implemented.

[0030] The mining parameter optimization method proposed in this application optimizes the mining parameters during the mining process and can realize the quantitative determination of the mining parameters. It can not only ensure the mining efficiency and achieve the maximum loss reduction mining, but also control the mining damage from the source and significantly reduce the degree of overburden damage in shallow buried high-intensity mining. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] The drawings constituting a part of the present application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0032] Figure 1 A flowchart of a mining parameter optimization method according to an exemplary embodiment of the present application;

[0033] Figure 2 It is a flow chart of a mining parameter optimization method according to a specific embodiment of the present application;

[0034] Figure 3 This is a schematic diagram of the distribution of the matrix and the damaged body in a unit body according to a specific embodiment of the present application;

[0035] Figure 4 A relationship model for the relationship between the overburden damage ratio of the overburden stratum and the mining height of the working face according to a specific embodiment of the present application;

[0036] Figure 5 A relationship model for the relationship between the overburden damage ratio of the overburden formation and the working face length according to a specific embodiment of the present application;

[0037] Figure 6 It is a relationship model of the relationship between the overburden damage ratio of the overburden formation and the working face advancement speed according to a specific embodiment of the present application. DETAILED DESCRIPTION

[0038] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0039] Embodiment 1

[0040] This embodiment provides a mining parameter optimization method. Figure 1 FIG. 1 is a flow chart of a mining parameter optimization method according to an exemplary embodiment of the present application. Figure 1 As shown, the method comprises the following steps:

[0041] S100: Obtain the volume of the overburden stratum, and obtain the volume of damaged rock mass in the overburden stratum caused by mining the mineral deposits under the overburden stratum according to different parameter value combinations of multiple mining parameters, and use the ratio of the volume of the damaged rock mass to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum.

[0042] Among them, the multiple mining parameters include: working face mining height, working face length and working face advancement speed.

[0043] Before obtaining the volume of the overburden stratum and obtaining the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposits beneath the overburden stratum according to different parameter value combinations of multiple mining parameters, the method may further include the step of: using a finite element analysis method to construct a virtual model of the overburden stratum after mining the mineral deposits beneath the overburden stratum according to different parameter value combinations of multiple mining parameters.

[0044] Obtaining the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of multiple mining parameters may include:

[0045] The first step is to divide the overburden stratum into several units and obtain the actual shear stress borne by each unit.

[0046] In the second step, for each unit body, the strain tensor of the elastic matrix in the unit body is determined by strain analysis, and the theoretical shear stress that the unit body can withstand is determined based on the strain tensor.

[0047] In the second step, the theoretical shear stress that the unit body can withstand is determined according to the strain tensor, which can include: using a preset damage constitutive model to determine the theoretical normal stress that the unit body can withstand according to the strain tensor of the elastic matrix in the unit body; using the Moore-Coulomb criterion to determine the theoretical shear stress that the unit body can withstand according to the theoretical normal stress that each unit body can withstand.

[0048] Wherein, the damage constitutive model can be:

[0049]

[0050] Among them, σ ij is the theoretical normal stress that the unit can withstand, G is the shear modulus of the overlying stratum, G 0 is the shear modulus of the elastic matrix in the unit cell, E ijkl is the elastic constant tensor of the elastic matrix in the unit cell, ε kl is the strain tensor of the elastic matrix in the unit cell, E mmkl is the elastic modulus of the elastic matrix in the unit cell, δij is the Kronecker symbol.

[0051] The Mohr-Coulomb criterion may be:

[0052] τ=c+σ n tanφ

[0053] Wherein, c is the cohesion between each unit body, Φ is the internal friction angle between each unit body, σ n and τ are respectively the theoretical normal stress and theoretical shear stress that the unit cell can withstand.

[0054] The third step is to compare the actual shear stress borne by the unit body with the theoretical shear stress that the unit body can withstand, and determine whether the unit body is damaged based on the comparison result.

[0055] The fourth step is to obtain the volume of each damaged unit body, and sum the volumes of each damaged unit body, and use the sum as the volume of the damaged rock mass in the overburden stratum.

[0056] S200: For each mining parameter, a relationship model of the relationship between the overburden damage ratio of the overburden stratum and the mining parameter is established according to the corresponding relationship between the overburden damage ratio and different parameter values ​​of the mining parameter.

[0057] Wherein, establishing a relationship model between the overburden damage ratio of the overburden formation and the mining parameter may include:

[0058] A relationship model between the overburden damage ratio of the overburden stratum and the mining parameters is established by using the Sigmoid function fitting method.

[0059] S300: Determine the maximum value of the change rate of the overburden damage ratio according to the relationship model, and use the parameter value of the mining parameter corresponding to the maximum value of the change rate as the endpoint value of the parameter value range after the mining parameter is optimized, and then determine the endpoint value of the parameter value range after the mining parameter is optimized for each of the multiple mining parameters.

[0060] Determining the maximum value of the change rate of the overburden damage ratio according to the relationship model, and taking the parameter value of the mining parameter corresponding to the maximum value of the change rate as the endpoint value of the parameter value range after the mining parameter is optimized, including:

[0061] When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face mining height, the maximum increase rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face mining height corresponding to the maximum increase rate is used as the upper limit value of the parameter value range after the working face mining height is optimized.

[0062] Among them, the relationship model between the overburden damage ratio of the overburden strata and the mining height of the working face can be:

[0063]

[0064] Among them, D f represents the overburden damage ratio of the overburden strata, M represents the mining height of the working face, M1 represents the parameter value of the mining height of the working face corresponding to the maximum value of the increase rate of the overburden damage ratio of the overburden strata, and a, b and c are all fitting constants.

[0065] When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face length, the maximum increase rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face length corresponding to the maximum increase rate is used as the upper limit value of the parameter value range after the working face length is optimized.

[0066] Among them, the relationship model between the overburden damage ratio of the overburden stratum and the working face length can be:

[0067]

[0068] Among them, D f represents the overburden damage ratio of the overburden strata, L represents the length of the working face, L1 represents the parameter value of the working face length corresponding to the maximum increase rate of the overburden damage ratio of the overburden strata, and a1, b1 and c1 are all fitting constants.

[0069] When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face advancement speed, the maximum value of the reduction speed of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face advancement speed corresponding to the maximum value of the reduction speed is used as the lower limit value of the parameter value range after the working face advancement speed is optimized.

[0070] Among them, the relationship model between the overburden damage ratio of the overburden strata and the working face advancement speed can be:

[0071]

[0072] Among them, D f represents the overburden damage ratio of the overburden strata, v represents the working face advancement speed, v1 represents the parameter value of the working face advancement speed corresponding to the maximum value of the reduction speed of the overburden damage ratio of the overburden strata, and a2, b2 and c2 are all fitting constants.

[0073] The mining parameter optimization method proposed in this application optimizes the mining parameters during the mining process and can realize the quantitative determination of the mining parameters. It can not only ensure the mining efficiency and achieve the maximum loss reduction mining, but also control the mining damage from the source and significantly reduce the degree of overburden damage in shallow buried high-intensity mining.

[0074] Embodiment 2

[0075] The purpose of this embodiment is to provide a mining parameter optimization method. Figure 2 The present invention is a flowchart of a mining parameter optimization method according to a specific embodiment of the present application.

[0076] like Figure 2 As shown, the specific steps of this embodiment are as follows:

[0077] (1) Define the damage variable of the unit cell inside the rock mass as G is the shear modulus of the rock mass, G 0 is the shear modulus of the elastic matrix within the unit cell.

[0078] (2) Construction of damage constitutive model based on elastoplastic mechanics and damage mechanics Among them, σ ij is the theoretical normal stress that the unit body in the rock mass can withstand, E ijkl is the elastic constant tensor of the elastic matrix in the unit cell, ε kl is the strain tensor of the elastic matrix in the unit cell, E mmkl is the elastic modulus of the elastic matrix within the unit cell.

[0079] The damage constitutive model can be constructed by the following method:

[0080] Rock mass engineering practice and related experimental studies have shown that under most engineering conditions, the deformation and strength characteristics of rock mass belong to the category of brittle failure. Based on this, the following assumptions are made: (a) The rock mass is composed of two parts: the matrix (the part without cracks) and the damaged body (the part with microcracks); (b) The matrix is ​​an isotropic elastic medium, and elastic deformation does not cause rock mass damage; (c) The damaged body is a rigid-plastic body with no yield strength; (d) Hydrostatic pressure will not cause rock mass damage; (e) The matrix and the damaged body satisfy deformation coordination, that is, the strains are equal.

[0081] Figure 3 FIG. 1 is a schematic diagram showing the distribution of the matrix and the damaged body in a unit according to a specific embodiment of the present application. Figure 3 As shown in Figure 1, first, the damage variable of the unit cell inside the rock mass can be defined as:

[0082]

[0083] In expression (1), D is the damage variable, dV is the unit volume, dV0 is the volume of the matrix inside the unit, and dV D is the volume of the damaged body within the unit cell.

[0084] The damaged part in the rock mass is the stress release zone, and its deviatoric stress is zero. Under the hydrostatic pressure state, the cracks in the rock mass tend to close, and no obvious damage occurs. It can be considered that the damage to the rock mass is mainly caused by deviatoric stress, and the rock mass stress and matrix stress are located in the same deviatoric stress plane:

[0085]

[0086]

[0087]

[0088] In expressions (2) to (4), J2 is the second invariant of the deviatoric stress tensor of the rock mass and matrix, S ij is the deviatoric stress tensor of the rock mass, σ mm is the hydrostatic pressure, δ ij is the Kronecker symbol.

[0089]

[0090]

[0091] In expressions (5) to (6), represents the second invariant of the deviatoric stress tensor of the matrix in the unit cell, is the deviatoric stress tensor of the matrix in the unit cell, It represents the theoretical normal stress that the matrix in the unit cell can withstand.

[0092] Substituting expressions (3) to (6) into expression (2), we can obtain the theoretical normal stress σ that the unit cell can withstand: ij for:

[0093]

[0094] According to the assumption that hydrostatic pressure will not cause damage to the rock mass, the hydrostatic pressure is equal to the pressure on the matrix, that is, Then we rewrite expression (7) as:

[0095]

[0096] Based on the assumption that the matrix is ​​an isotropic elastic medium, we have:

[0097]

[0098] In expression (9), E ijkl is the elastic constant tensor of the matrix. ε kl is the strain tensor of the matrix. Based on the assumption that the matrix and the damaged body satisfy deformation coordination, we can obtain:

[0099]

[0100] Substituting expressions (9) and (10) into expression (8), we obtain:

[0101]

[0102] Expression (11) is the constitutive model of rock damage.

[0103] However, the volume of the damaged body in the rock mass is virtual and cannot be measured, which makes it unsuitable for practical applications. Combining expressions (2) and (10), we can get:

[0104] G=(1-D)G 0 (12)

[0105]

[0106] In expressions (12) and (13), G is the shear modulus of the rock mass, which varies with the degree of damage and is a variable. 0 is the shear modulus of the matrix, which is a constant.

[0107] In elastic-plastic mechanics, the relationship between equivalent stress and equivalent strain is as follows:

[0108]

[0109] According to expression (14), the variation curve of rock mass shear modulus G can be obtained from mechanical tests. Substituting expression (13) into expression (11), we get:

[0110]

[0111] Expression (15) is the final constitutive model of rock damage. Compared with expression (11), expression (15) is more convenient for practical application.

[0112] (3) The damage constitutive model based on damage variables is used to secondary develop the Mohr-Coulomb constitutive model of the three-dimensional numerical simulation software.

[0113] (4) Definition of overburden damage ratio suitable for simulation analysis V i Represents the volume of a plastic unit, V G is the volume of the observed space.

[0114] (5) Simulate different working face mining heights M, different working face lengths L, and different working face advancement speeds v, and calculate the overburden damage ratio based on the simulation results.

[0115] The simulation process includes the following steps: 1) constructing a model based on the actual geological conditions of the working face mining; 2) sampling each rock layer in the overburden stratum to measure the physical and mechanical parameters of each rock layer; 3) converting the rock sample parameters measured in the laboratory into mechanical parameters for numerical simulation; 4) conducting the corresponding working face mining simulation.

[0116] (6) In the calculation process of overburden damage ratio, the total overburden volume is the volume of all units directly above the mining working face, and the damage volume is the volume of the unit that undergoes plastic yield in the total volume.

[0117] (7) The FISH language is used to determine the state (elastic / inelastic) of each unit body within the observation range, and the total damage volume and the total volume of the observation area are calculated respectively.

[0118] (8) Determine the theoretical normal stress that the unit body can withstand based on the damage constitutive model in (2), determine the theoretical shear stress that the unit body can withstand based on the theoretical normal stress that the unit body can withstand using the Mohr-Coulomb criterion, compare the actual shear stress that the unit body withstands with the theoretical shear stress, and determine that the unit body is damaged, i.e., plastic deformation occurs, when the actual shear stress that the unit body withstands is greater than the theoretical shear stress.

[0119] Among them, the expression of the Mohr-Coulomb criterion is τ=c+σ n tanφ, where c is the cohesion, Φ is the internal friction angle, σ n and τ are the theoretical normal stress and theoretical shear stress that the unit cell can withstand, respectively.

[0120] Then, according to the calculation formula of overburden damage ratio in (4), the FISH language is used to obtain the overburden damage ratio of different working face mining heights, different working face lengths and different working face advancement speeds.

[0121] (9) According to the simulation results, the Sigmoid functions between the working face mining height, working face length, working face advancement speed and overburden damage ratio are obtained as follows:

[0122]

[0123] and

[0124] Among them, a, b, c, a1, b1, c1, a2, b2 and c2 are fitting constants.

[0125] Taking the typical working face of Shendong mining area as the engineering background, the corresponding fitting results can be obtained as follows:

[0126] D f =0.0582+0.4095 / (1+e -(M-7.8445) / 1.9979 )

[0127] D f =0.0241+0.9432 / (1+e -(L-309.8928) / 66.7393 )

[0128] D f =0.1719+2.0603 / (1+e -(v-11.9551) / 2.6909 )

[0129] Correspondingly, the curve diagram of the relationship model between the overburden damage ratio of the overburden strata and the mining height of the working face is as follows: Figure 4 As shown in the figure, the curve diagram of the relationship model between the overburden damage ratio of the overburden strata and the working face length is as follows Figure 5 As shown in the figure, the curve diagram of the relationship model between the overburden damage ratio of the overburden strata and the working face advancement speed is as follows Figure 6 shown.

[0130] Specifically, the working face length and the working face advancement speed can be set to commonly used fixed values ​​first, so as to determine the Sigmoid function between the working face mining height and the overburden damage ratio. After obtaining the optimized working face mining height, the working face mining height can be set to a fixed parameter value of the optimized working face mining height, and the working face advancement speed can be set to a commonly used fixed value, so as to determine the Sigmoid function between the working face length and the overburden damage ratio. After obtaining the optimized working face mining height and working face length, both the working face mining height and the working face length can be set to the optimized fixed parameter values, so as to determine the Sigmoid function between the working face advancement speed and the overburden damage ratio.

[0131] (10) MATABLE is used to calculate the working face mining height M1 and working face length L1 corresponding to the maximum value of the overburden damage ratio increase rate, and the working face advancement speed v1 corresponding to the maximum value of the damage ratio decrease rate is calculated.

[0132] Taking the fitting results obtained under the typical working background of Shenhua Mining Area as an example, MATABLE is used to calculate that the maximum value of the damage ratio increase rate corresponds to a working face mining height of 7.8445 and a working face length of 309.8928. At the same time, the maximum value of the overburden damage ratio decrease rate is calculated to correspond to a working face advancement speed of 11.9551.

[0133] (11) Calculate the maximum value of the working face advancement speed based on the advancement speed of the coal mining machine Among them, t is the coal cutting time of the coal shearer every day, v' is the coal cutting speed of the coal shearer, and d is the cutting depth of the coal shearer.

[0134] Still taking the Shenhua mining area as an example, the maximum propulsion speed of the coal mining machine is calculated to be 20.58.

[0135] (12) The optimized working face mining height range is finally determined to be 0~M1, the optimized working face length range is 0~L1, and the optimized working face advancement speed range is v1~v max .

[0136] Still taking the Shenhua mining area as an example, according to the above results, it is finally determined that the range of working face mining height is 0-7.8445m, and the range of working face length is 0-309.8928m; the reasonable range of working face advancement speed is 11.9551-20.58m / d.

[0137] For high-intensity intensive coal mining, increasing the height, length and advancement speed of the working face is conducive to large-scale mining, reducing costs and improving efficiency. Figures 4 to 6 As shown in the figure, the overburden damage ratio is positively correlated with the working face mining height and length, and negatively correlated with the advancement speed. Therefore, it is advisable to control the damage ratio increase speed below the maximum value or the reduction speed above the maximum value under the premise of ensuring mining efficiency. Therefore, the working face mining height and working face length corresponding to the maximum value of the overburden damage ratio increase speed are taken as the corresponding upper limit values, and the working face advancement speed corresponding to the maximum value of the overburden damage ratio reduction speed is taken as the corresponding lower limit value.

[0138] The main steps of the mining parameter optimization method of this embodiment include: first, constructing a damage constitutive model and performing secondary development of the numerical simulation software; second, proposing an overburden damage ratio model to quantitatively analyze the degree of overburden damage; third, calculating the degree of overburden damage under different mining heights, different working face lengths and different working face advancement speeds through numerical simulation methods.

[0139] The method of this embodiment uses Sigmoid function fitting to obtain fitting formulas between the overburden damage ratio and different working face mining heights, different working face lengths, and different working face advancement speeds. The optimized value ranges of the working face mining height and the working face length are less than or equal to the corresponding parameter values ​​corresponding to the fastest increase rate of the overburden damage ratio. The optimized value range of the working face advancement speed is the corresponding parameter value corresponding to the fastest decrease in the overburden damage ratio and the limit of the coal cutting speed of the coal mining machine.

[0140] In view of the severe overburden damage and serious surface destruction caused by shallow buried high-intensity mining in western mining areas, the present invention proposes a high-intensity loss-reducing mining parameter optimization method based on overburden damage ratio. This method can not only ensure mining efficiency, but also significantly reduce the degree of overburden damage caused by shallow buried high-intensity mining from the source, thereby reducing surface damage from the source, which is conducive to realizing green mining in ecologically fragile areas in western mining areas.

[0141] Embodiment 3

[0142] This embodiment provides a mining method, comprising the following steps: mining the mineral deposits under the overburden stratum according to the mining parameters optimized by the mining parameter optimization method as described above.

[0143] Embodiment 4

[0144] This embodiment provides a storage medium storing a computer program. When the computer program is executed by a processor, the steps of the mining parameter optimization method described above are implemented.

[0145] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 The steps of a specified function in a process or multiple processes.

[0146] Storage media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, modules of programs or other data. Examples of storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information that can be accessed by a computing device.

[0147] Embodiment 5

[0148] This embodiment provides a computer device, including a processor and a storage medium storing program codes. When the program codes are executed by the processor, the steps of the mining parameter optimization method described above are implemented.

[0149] In one embodiment, a computer device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0150] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash memory (FLASH RAM). The memory is an example of a computer-readable medium.

[0151] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. When the terms "include" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or their combinations.

[0152] It should be understood that the exemplary embodiments in this specification can be implemented in a variety of different forms and should not be interpreted as being limited to the embodiments described herein. These embodiments are provided to make the disclosure of this application thorough and complete, and to fully convey the concepts of these exemplary embodiments to those of ordinary skill in the art, and should not be interpreted as limiting the present invention.

Claims

1. A mining parameter optimization method, characterized in that: The following steps are involved: Obtaining the volume of the overburden stratum, and obtaining the volume of the rock mass damaged in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of multiple mining parameters, and taking the ratio of the volume of the damaged rock mass to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum; For each mining parameter, a relationship model of the relationship between the overburden damage ratio of the overburden formation and the mining parameter is established according to the corresponding relationship between the overburden damage ratio and different parameter values ​​of the mining parameter, wherein the multiple mining parameters include: working face mining height, working face length and working face advancement speed; The relationship model between the overburden damage ratio of the overburden strata and the mining height of the working face is: in, represents the overburden damage ratio of the overburden strata, M represents the mining height of the working face, M1 represents the parameter value of the mining height of the working face corresponding to the maximum value of the increase rate of the overburden damage ratio of the overburden strata, and a, b and c are all fitting constants; The relationship model between the overburden damage ratio of the overburden strata and the length of the working face is: in, represents the overburden damage ratio of the overburden strata, L represents the length of the working face, L1 represents the parameter value of the working face length corresponding to the maximum value of the increase rate of the overburden damage ratio of the overburden strata, and a1, b1 and c1 are all fitting constants; The relationship model between the overburden damage ratio of the overburden strata and the working face advancement speed is: in, represents the overburden damage ratio of the overburden strata, v represents the working face advancement speed, v1 represents the parameter value of the working face advancement speed corresponding to the maximum value of the overburden damage ratio reduction speed of the overburden strata, and a2, b2 and c2 are all fitting constants; The maximum value of the change rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the mining parameter corresponding to the maximum value of the change rate is used as the endpoint value of the parameter value range after the mining parameter is optimized, thereby determining the endpoint value of the parameter value range after the optimization of each mining parameter in the multiple mining parameters.

2. The mining parameter optimization method according to claim 1, characterized in that: Obtaining the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of multiple mining parameters, including: Divide the overburden strata into several units and obtain the actual shear stress borne by each unit; For each unit body, the strain tensor of the elastic matrix in the unit body is determined by strain analysis, and the theoretical shear stress that the unit body can withstand is determined according to the strain tensor; Compare the actual shear stress borne by the unit body with the theoretical shear stress that the unit body can withstand, and determine whether the unit body is damaged according to the comparison result; The volumes of the damaged unit bodies are obtained, and the volumes of the damaged unit bodies are summed, and the sum is used as the volume of the damaged rock mass in the overburden stratum.

3. The mining parameter optimization method according to claim 2, characterized in that: Determining the theoretical shear stress that the unit cell can withstand according to the strain tensor includes: Using a preset damage constitutive model, according to the strain tensor of the elastic matrix in the unit body, the theoretical normal stress that the unit body can withstand is determined; Using the Mohr-Coulomb criterion, the theoretical shear stress that each unit cell can withstand is determined based on the theoretical normal stress that the unit cell can withstand.

4. The mining parameter optimization method according to claim 3, characterized in that: The damage constitutive model is: in, is the theoretical normal stress that the unit can withstand, G is the shear modulus of the overlying strata, is the shear modulus of the elastic matrix in the unit cell, Eijkl is the elastic constant tensor of the elastic matrix in the unit cell, εkl is the strain tensor of the elastic matrix in the unit cell, Emmkl is the elastic modulus of the elastic matrix in the unit cell, is the Kronecker symbol; The Mohr-Coulomb criterion is: Wherein, c is the cohesion between each of the unit bodies, Φ is the internal friction angle between each of the unit bodies, σn and τ are respectively the theoretical normal stress and theoretical shear stress that the unit body can withstand.

5. The mining parameter optimization method according to claim 1, characterized in that: A relationship model between the overburden damage ratio of the overburden formation and the mining parameters is established, including: A relationship model between the overburden damage ratio of the overburden stratum and the mining parameters is established by using the Sigmoid function fitting method.

6. The mining parameter optimization method according to claim 1, characterized in that: Determining the maximum value of the change rate of the overburden damage ratio according to the relationship model, and taking the parameter value of the mining parameter corresponding to the maximum value of the change rate as the endpoint value of the parameter value range after the mining parameter is optimized, including: When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden stratum and the working face mining height, the maximum value of the increase rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face mining height corresponding to the maximum value of the increase rate is used as the upper limit value of the parameter value range after the working face mining height is optimized; When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden stratum and the working face length, the maximum value of the increase rate of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face length corresponding to the maximum value of the increase rate is used as the upper limit value of the parameter value range after the working face length is optimized; When the relationship model is a relationship model about the relationship between the overburden damage ratio of the overburden strata and the working face advancement speed, the maximum value of the reduction speed of the overburden damage ratio is determined according to the relationship model, and the parameter value of the working face advancement speed corresponding to the maximum value of the reduction speed is used as the lower limit value of the parameter value range after the working face advancement speed is optimized.

7. The mining parameter optimization method according to claim 1, characterized in that: Before obtaining the volume of the overburden stratum and obtaining the volume of rock mass damaged in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of multiple mining parameters, the method further includes the steps of: By using the finite element analysis method, a virtual model of the overburden stratum after the mineral deposit under the overburden stratum is mined according to different parameter value combinations of multiple mining parameters is constructed.

8. A mining method, characterized in that: The following steps are involved: The mineral deposits under the overburden strata are mined according to the mining parameters optimized by the mining parameter optimization method according to any one of claims 1 to 7.

9. A storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the mining parameter optimization method according to any one of claims 1 to 7 are implemented.

10. A computer device comprising a processor and a storage medium storing program codes, wherein when the program codes are executed by the processor, the steps of the mining parameter optimization method according to any one of claims 1 to 7 are implemented.

Images