CN114818399B - Mining parameter optimization method, mining method, storage medium and computer device - Google Patents
Mining parameter optimization method, mining method, storage medium and computer device Download PDFInfo
- Publication number
- CN114818399B CN114818399B CN202110064628.8A CN202110064628A CN114818399B CN 114818399 B CN114818399 B CN 114818399B CN 202110064628 A CN202110064628 A CN 202110064628A CN 114818399 B CN114818399 B CN 114818399B
- Authority
- CN
- China
- Prior art keywords
- overburden
- mining
- parameter
- damage ratio
- stratum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
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
With the progress of technology and equipment, the scale of the coal mine for well engineering in China is developed from the megaton level of ninety years in the last century (100-300 ten thousand tons/year) to the megaton level of the first era (1000-3000 ten thousand tons/year). The intensive level of exploitation is continuously improved, and the method becomes an important way for improving the safety guarantee degree, the resource recovery rate and the economic benefit of the coal mine gradually. According to the intensive development layout, the yield of 14 large coal bases in China accounts for over 95 percent of the whole country.
High-strength mining is an important mode of intensive mining, but simultaneously brings problems of rapid surface subsidence and large deformation, and further causes serious ecological damage. The occurrence of coal in China and the capacity distribution of ecological environment are unbalanced. The yield of Jinshan Mongolian sweet coal in the west represented by the mine area is 3/4 of that of the mine area in China, but the mine area is ecologically fragile and water resource deficient, and desertification of the mine area and the periphery is aggravated by high-strength exploitation. The eastern mining area represented by the two-Huai mining area guarantees the energy supply in the east China area, but the land collapse is caused by the high-strength mining, so that the land resources of mining cities are more tense. The ecological damage of the ground surface of the mining area is caused by overlying strata deformation and fracture caused by coal exploitation, mining damage of lower strata caused by stratum movement and transmission damage of upper strata and ground surface layers caused by the mining damage. Therefore, the root of ecological restoration in mining areas is to furthest reduce the damage to the surface ecology from the exploitation source.
In conclusion, how to realize coordination of efficient restoration of the surface ecology and high-strength exploitation is a great technical problem in the coal industry. However, due to the reasons of rapid high-strength exploitation, large exploitation working surface size, large overlying strata damage range, rapid earth surface sedimentation speed, large deformation, high degree of damage to earth surface and ecology, outstanding contradiction between resource development and ecological mining area construction and the like, the earth surface damage degree is difficult to control from the source. For this reason, optimization of the mining parameters of the working face from the mining source is needed, so that high-intensity coal mining and surface ecological protection are coordinated.
Disclosure of Invention
The invention mainly aims to provide a mining parameter optimization method, a mining method, a storage medium and computer equipment, so as to optimize high-strength mining parameters and realize coordination of high-strength mining and surface ecological protection.
The application provides a mining parameter optimization method which comprises the steps of obtaining the volume of a overburden stratum, obtaining the volume of a rock mass damaged in the overburden stratum caused by mining a mineral deposit under the overburden stratum according to different parameter value combinations of a plurality of mining parameters, taking the ratio of the damaged rock mass volume to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum, establishing a relation model about the relation between the overburden damage ratio of the overburden stratum and the mining parameters according to the corresponding relation between the overburden damage ratio and the different parameter values of the mining parameters for each mining parameter, determining the maximum value of the change speed of the overburden damage ratio according to the relation model, 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 further determining the endpoint value of the parameter value range after the mining parameter is optimized for each mining parameter of the plurality of mining parameters.
In one embodiment, obtaining the volume of rock mass damaged in the overburden formation caused by mining the mineral reserve under the overburden formation according to different parameter value combinations of a plurality of mining parameters comprises dividing the overburden formation into a plurality of unit bodies, obtaining actual shear stress born by each unit body, determining the strain tensor of an elastic matrix in the unit body according to strain analysis aiming at each unit body, determining the theoretical shear stress born by the unit body according to the strain tensor, comparing the actual shear stress born by the unit body with the theoretical shear stress born by the unit body, determining whether the unit body is damaged according to a comparison result, obtaining the volume of each unit body damaged, summing the volumes of each unit body damaged, and taking the sum value as the volume of the rock mass damaged in the overburden formation.
In one embodiment, the theoretical shear stress which can be borne by the unit body is determined according to the strain tensor, which comprises the steps of determining the theoretical normal stress which can be borne by the unit body according to the strain tensor of the elastic matrix in the unit body by utilizing a preset damage constitutive model, and determining the theoretical shear stress which can be borne by the unit body according to the theoretical normal stress which can be borne by each unit body by utilizing a molar coulomb criterion.
In one embodiment, the damage constitutive model is:
Wherein σ ij is the theoretical normal stress that the unit body can bear, G is the shear modulus of the overburden stratum, G 0 is the shear modulus of the elastic matrix in the unit body, E ijkl is the elastic constant tensor of the elastic matrix in the unit body, ε kl is the strain tensor of the elastic matrix in the unit body, E mmkl is the elastic modulus of the elastic matrix in the unit body, and δ ij is Kronecker symbol;
The molar coulomb criterion is:
τ=c+σntanφ
wherein c is the cohesive force between the unit bodies, phi is the internal friction angle between the unit bodies, and sigma n and tau are the theoretical normal stress and the theoretical shear stress which can be born by the unit bodies respectively.
In one embodiment, establishing a relationship model of the relationship between the overburden damage ratio and the production parameter for the overburden formation includes establishing a relationship model of the relationship between the overburden damage ratio and the production parameter for the overburden formation using a Sigmoid function fitting method.
In one embodiment, the plurality of mining parameters comprise a working face mining height, a working face length and a working face advancing speed, a maximum value of a change speed of a cover rock damage ratio is determined according to the relation model, the parameter value of the mining parameter corresponding to the maximum value of the change speed is used as an end point value of a parameter value range after optimization of the mining parameter, the method comprises the steps of determining an increasing speed maximum value of the cover rock damage ratio according to the relation model when the relation model is a relation model of a relation between the cover rock damage ratio and the working face mining height of a cover rock stratum, taking the parameter value of the working face mining height corresponding to the increasing speed maximum value as an upper limit value of the parameter value range after optimization of the working face length, determining an increasing speed maximum value of the cover rock damage ratio according to the relation model when the relation model is a relation model of the relation between the cover rock damage ratio and the working face length of the cover rock damage ratio, taking the parameter value of the parameter value corresponding to the increasing speed maximum value as an upper limit value of the parameter value after optimization of the working face length of the cover rock damage ratio when the relation model is a relation model of the cover rock damage ratio between the cover rock damage ratio and the working face length of the working face is reduced, and the working face advancing speed of the parameter range is determined according to the relation model.
In one embodiment, the relationship model for the relationship between overburden damage ratio and face take height for a overburden formation is:
Wherein D f represents the overburden damage ratio of the overburden stratum, M represents the working face mining height, M 1 represents the parameter value of the working face mining height corresponding to the maximum increase speed of the overburden damage ratio of the overburden stratum, and a, b and c are fitting constants;
the relation model about the relation between the overburden damage ratio and the working surface length of the overburden stratum is:
Wherein D f represents the overburden damage ratio of the overburden formation, L represents the working face length, L 1 represents the parameter value of the working face length corresponding to the maximum value of the rate of increase of the overburden damage ratio of the overburden formation, and a 1、b1 and c 1 are both fitting constants;
the relation model about the relation between the overburden damage ratio and the working surface advancing speed of the overburden stratum is as follows:
Wherein D f represents the overburden damage ratio of the overburden formation, v represents the face advancing speed, v 1 represents the parameter value of the face advancing speed corresponding to the maximum value of the decreasing speed of the overburden damage ratio of the overburden formation, and a 2、b2 and c 2 are fitting constants.
In one embodiment, prior to obtaining the volume of the overburden formation and obtaining the volume of rock mass damaged in the overburden formation resulting from mining the mineral reserve in the overburden formation in accordance with the different parameter value combinations of the plurality of mining parameters, the method further includes the step of constructing a virtual model of the overburden formation after mining the mineral reserve in the overburden formation in accordance with the different parameter value combinations of the plurality of mining parameters using a finite element analysis method.
In a second aspect, the present application provides a method of mining comprising mining a mineral deposit in a overburden formation in accordance with mining parameters optimized by a mining parameter optimization method as described above.
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.
In a fourth aspect, the present application provides a computer device comprising a processor and a storage medium storing program code which, when executed by the processor, performs the steps of the production parameter optimization method as described above.
According to the mining parameter optimization method provided by the application, the mining parameters in the mining process are optimized, the quantitative determination of the mining parameters can be realized, the mining benefit can be ensured, the damage-reducing mining can be realized to the greatest extent, the mining damage can be controlled from the source, and the damage degree of shallow-buried high-strength mining overlying strata can be greatly reduced.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a undue limitation on the application, wherein:
FIG. 1 is a flow chart of a method of optimizing production parameters according to an exemplary embodiment of the application;
FIG. 2 is a flow chart of a method of optimizing production parameters according to an embodiment of the present application;
FIG. 3 is a schematic diagram of matrix and lesion distribution in a cell according to an embodiment of the application;
FIG. 4 is a model of the relationship between overburden damage ratio and face elevation for a overburden formation in accordance with one embodiment of the present application;
FIG. 5 is a model of the relationship between overburden damage ratio and length of a working surface for a overburden formation in accordance with one embodiment of the present application;
FIG. 6 is a model of the relationship between overburden damage ratio and face advance rate for a overburden formation in accordance with one embodiment of the present application.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
Example 1
The present embodiment provides a method for optimizing mining parameters, and fig. 1 is a flowchart of a method for optimizing mining parameters according to an exemplary embodiment of the present application. As shown in fig. 1, the method comprises the steps of:
And S100, acquiring the volume of the overburden stratum, acquiring the damaged rock volume in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of a plurality of mining parameters, and taking the ratio of the damaged rock volume to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum.
The plurality of mining parameters include face mining height, face length and face propulsion speed.
Before acquiring the volume of the overburden formation and acquiring the volume of the damaged rock mass in the overburden formation caused by mining the mineral deposit under the overburden formation in accordance with the different parameter value combinations of the plurality of mining parameters, the method may further include the step of constructing a virtual model of the overburden formation after mining the mineral deposit under the overburden formation in accordance with the different parameter value combinations of the plurality of mining parameters using a finite element analysis method.
The method for obtaining the damaged rock mass volume in the overburden stratum caused by mining the mineral deposit under the overburden stratum according to different parameter value combinations of a plurality of mining parameters can comprise the following steps:
And firstly, dividing the overlying strata into a plurality of unit bodies, and obtaining the actual shear stress born by each unit body.
And secondly, determining the strain tensor of the elastic matrix in each unit body through strain analysis, and determining the theoretical shear stress which can be born by the unit body according to the strain tensor.
In the second step, the theoretical shear stress which can be born by the unit body is determined according to the strain tensor, which can comprise the steps of determining the theoretical normal stress which can be born by the unit body according to the strain tensor of the elastic matrix in the unit body by utilizing a preset damage constitutive model and determining the theoretical shear stress which can be born by the unit body according to the theoretical normal stress which can be born by each unit body by utilizing a molar coulomb criterion.
Wherein, the damage constitutive model can be:
wherein σ ij is the theoretical normal stress that the unit body can bear, G is the shear modulus of the overburden stratum, G 0 is the shear modulus of the elastic matrix in the unit body, E ijkl is the elastic constant tensor of the elastic matrix in the unit body, ε kl is the strain tensor of the elastic matrix in the unit body, E mmkl is the elastic modulus of the elastic matrix in the unit body, and δ ij is the Kronecker symbol.
The molar coulomb criterion may be:
τ=c+σntanφ
wherein c is the cohesive force between the unit bodies, phi is the internal friction angle between the unit bodies, and sigma n and tau are the theoretical normal stress and the theoretical shear stress which can be born by the unit bodies respectively.
And thirdly, comparing the actual shear stress born by the unit body with the theoretical shear stress born by the unit body, and determining whether the unit body is damaged or not according to a comparison result.
And step four, obtaining the volume of each damaged unit body, summing the volumes of each damaged unit body, and taking the sum value as the volume of the damaged rock body in the overburden stratum.
And S200, for each mining parameter, establishing a relation model about the relation between the overburden damage ratio of the overburden stratum and the mining parameter according to the corresponding relation between the overburden damage ratio and different parameter values of the mining parameter.
Wherein establishing a relationship model regarding a relationship between a overburden damage ratio of the overburden formation and the production parameter may include:
And establishing a relation model about the relation between the overburden damage ratio of the overburden stratum and the exploitation parameters by using a Sigmoid function fitting method.
And S300, determining the maximum value of the change speed of the overburden damage ratio according to the relation model, taking the parameter value of the mining parameter corresponding to the maximum value of the change speed as the end point value of the parameter value range after optimizing the mining parameter, and further determining the end point value of the parameter value range after optimizing each mining parameter in the plurality of mining parameters.
Determining a maximum value of the change speed of the overburden damage ratio according to the relation model, taking a parameter value of the mining parameter corresponding to the maximum value of the change speed as an endpoint value of a parameter value range after optimizing the mining parameter, and comprising the following steps:
When the relation model is a relation model of relation between the overburden damage ratio of the overburden stratum and the working face mining height, determining the maximum value of the increasing speed of the overburden damage ratio according to the relation model, and taking the parameter value of the working face mining height corresponding to the maximum value of the increasing speed as the upper limit value of the parameter value range after optimizing the working face mining height.
Wherein, the relation model about the relation between the overburden damage ratio and the working face mining height of the overburden stratum can be:
wherein D f represents the overburden damage ratio of the overburden formation, M represents the working face elevation, M 1 represents the parameter value of the working face elevation corresponding to the maximum rate of increase of the overburden damage ratio of the overburden formation, and a, b and c are fitting constants.
When the relation model is a relation model of relation between the overburden damage ratio and the working face length of the overburden stratum, determining the maximum increase speed of the overburden damage ratio according to the relation model, and taking the parameter value of the working face length corresponding to the maximum increase speed as the upper limit value of the parameter value range after optimizing the working face length.
Wherein, the relation model of the relation between the overburden damage ratio and the working surface length of the overburden stratum can be:
Wherein D f represents the overburden damage ratio of the overburden formation, L represents the working face length, L 1 represents the parameter value of the working face length corresponding to the maximum rate of increase of the overburden damage ratio of the overburden formation, and a 1、b1 and c 1 are both fitting constants.
When the relation model is a relation model of relation between the overburden damage ratio and the working face advancing speed of the overburden stratum, determining a maximum value of the decreasing speed of the overburden damage ratio according to the relation model, and taking a parameter value of the working face advancing speed corresponding to the maximum value of the decreasing speed as a lower limit value of a parameter value range after optimizing the working face advancing speed.
Wherein, the relation model about the relation between the overburden damage ratio and the working surface advancing speed of the overburden stratum can be:
Wherein D f represents the overburden damage ratio of the overburden formation, v represents the face advancing speed, v 1 represents the parameter value of the face advancing speed corresponding to the maximum value of the decreasing speed of the overburden damage ratio of the overburden formation, and a 2、b2 and c 2 are fitting constants.
According to the mining parameter optimization method provided by the application, the mining parameters in the mining process are optimized, the quantitative determination of the mining parameters can be realized, the mining benefit can be ensured, the damage-reducing mining can be realized to the greatest extent, the mining damage can be controlled from the source, and the damage degree of shallow-buried high-strength mining overlying strata can be greatly reduced.
Example two
It is an object of this embodiment to provide a method for optimizing production parameters, and fig. 2 is a flowchart of a method for optimizing production parameters according to an embodiment of the present application.
As shown in fig. 2, the specific steps of this embodiment are as follows:
(1) Defining the damage variable of the unit body inside the rock mass as G is the shear modulus of the rock mass and G 0 is the shear modulus of the elastomeric matrix in the cell.
(2) Construction of damage constitutive model based on elastoplastic mechanics and damage mechanicsWherein σ ij is the theoretical normal stress that the unit body in the rock mass can bear, E ijkl is the elastic constant tensor of the elastic matrix in the unit body, ε kl is the strain tensor of the elastic matrix in the unit body, and E mmkl is the elastic modulus of the elastic matrix in the unit body.
The damage constitutive model can be specifically constructed by the following method:
rock mass engineering practices and related experimental studies show that deformation and strength characteristics of rock mass under most engineering conditions belong to the category of brittle failure. Based on this, it is assumed that (a) the rock mass is composed of two parts of a matrix (crack-free part) and a damaged body (microcrack part), (b) the matrix is an isotropic elastic medium and elastic deformation does not cause damage to the rock mass, (c) the damaged body is a rigid body without yield strength, (d) hydrostatic pressure does not cause damage to the rock mass, and (e) the matrix and the damaged body satisfy deformation coordination, that is, strain equality.
FIG. 3 is a schematic diagram of matrix and lesion distribution in a cell according to an embodiment of the application. As shown in fig. 3, first, the damage variables of the rock mass internal unit body can be defined as:
In expression (1), D is a lesion variable, dV is a unit volume, dV 0 is a unit in-vivo matrix volume, and dV D is a unit in-vivo lesion volume.
The damaged portion in the rock mass is a stress relief zone with zero bias stress. Under the hydrostatic pressure state, the cracks in the rock body tend to be closed, and no obvious damage occurs. It is thus believed that damage to the rock mass is mainly caused by the bias stress, which lies in the same bias stress plane as the matrix stress:
In expressions (2) to (4), J 2 is the second invariant of the partial stress tensor of the rock mass and the matrix, S ij is the partial stress tensor of the rock mass, σ mm is the hydrostatic pressure, and δ ij is the Kronecker symbol.
In the expressions (5) to (6),Representing a second invariant of the deflection stress tensor of the matrix in the unit cell,Is the deflection stress tensor of the matrix in the unit body,Indicating the theoretical normal stress that the matrix in the cell body can withstand.
Substituting the expressions (3) - (6) into the expression (2) to obtain the theoretical normal stress sigma ij which can be born by the unit body, wherein the theoretical normal stress sigma ij is as follows:
according to the assumption that the hydrostatic pressure does not cause damage to the rock mass, so that the hydrostatic pressure is equal to the pressure borne by the matrix, i Further, expression (7) is rewritten as:
based on the assumption that the matrix is an isotropic elastic medium, there are:
In expression (9), E ijkl is the elastic constant tensor of the matrix. Epsilon kl is the strain tensor of the matrix. The assumption that the matrix and the damaged body meet deformation coordination is obtained by:
substituting expressions (9), (10) into expression (8) yields:
the expression (11) is the constitutive model of rock mass damage.
However, the volume of the damaged body in the rock mass is virtual, and cannot be measured, so that the method is not suitable for practical application. The combination of expressions (2), (10) can be obtained:
G=(1-D)G0 (12)
In expressions (12) and (13), G is the shear modulus of the rock mass, which varies with the degree of damage, and G 0 is the shear modulus of the matrix, which is a constant.
In elastoplastic mechanics, the equivalent stress and the equivalent strain have the following relationship:
according to expression (14), a variation curve of the shear modulus G of the rock mass can be obtained by a mechanical test. Substituting expression (13) into expression (11) yields:
Expression (15) is the constitutive model of the final rock mass damage. Expression (15) is more convenient for practical use than expression (11).
(3) And performing secondary development on the Mohr-Coulomb constitutive model of the three-dimensional numerical simulation software based on the damage constitutive model of the damage variable.
(4) Definition of overburden injury ratio suitable for simulation analysisV i represents the volume of one plastic unit and V G the volume of the observed space.
(5) Simulating different working surface mining heights M, different working surface lengths L and different working surface propelling speeds v, and calculating the overburden damage ratio according to simulation results.
The simulation process comprises the following steps of 1) carrying out model construction according to actual working face mining geological conditions, 2) sampling each stratum in a overlying strata to actually measure physical mechanical parameters of each stratum, 3) converting rock sample parameters measured in a laboratory into mechanical parameters of numerical simulation, and 4) carrying out corresponding working face mining simulation.
(6) In the process of calculating the damage ratio of the overlying strata, the total volume of the overlying strata is the volume of all corresponding unit bodies right above the mining working surface, and the damaged volume is the volume of the unit bodies which generate plastic yielding in the total volume.
(7) The states (elastic/inelastic) of the individual units within the observation range were determined using FISH language, and the total lesion volume and the total observation area volume were calculated, respectively.
(8) And (3) determining the theoretical normal stress which can be born by the unit body according to the damage constitutive model in the (2), determining the theoretical shear stress which can be born by the unit body according to the theoretical normal stress which can be born by the unit body by utilizing a molar coulomb criterion, comparing the actual shear stress born by the unit body with the theoretical shear stress, and judging that the unit body is damaged when the actual shear stress born by the unit body is larger than the theoretical shear stress, namely, generating plastic deformation.
Wherein the expression of the molar coulomb criterion is τ=c+σ n tan phi, wherein c is the cohesive force, phi is the internal friction angle, and σ n and τ are the theoretical normal stress and the theoretical shear stress that the unit can withstand, respectively.
And (3) obtaining the overburden damage ratio of different working surface mining heights, different working surface lengths and different working surface propelling speeds by using the FISH language according to the overburden damage ratio calculation formula in the step (4).
(9) The Sigmoid functions between the working face mining height, the working face length, the working face propelling speed and the overburden damage ratio are respectively obtained according to the simulation result:
And
Wherein a, b, c, a 1、b1、c1、a2、b2 and c 2 are both fitting constants.
By taking a typical working surface of the Shendong mining area as an engineering background, corresponding fitting results can be obtained as follows:
Df＝0.0582+0.4095/(1+e-(M-7.8445)/1.9979)
Df＝0.0241+0.9432/(1+e-(L-309.8928)/66.7393)
Df＝0.1719+2.0603/(1+e-(v-11.9551)/2.6909)
Accordingly, a graph of a relation model concerning the relation between the overburden damage ratio and the working face production height of the overburden formation is shown in fig. 4, a graph of a relation model concerning the relation between the overburden damage ratio and the working face length of the overburden formation is shown in fig. 5, and a graph of a relation model concerning the relation between the overburden damage ratio and the working face advancing speed of the overburden formation is shown in fig. 6.
Specifically, the length of the working surface and the advancing speed of the working surface can be set to be common fixed values, so that a Sigmoid function between the mining height of the working surface and the damage ratio of the overburden can be determined. After the optimized working surface mining height is obtained, the working surface mining height can be set to be a certain fixed parameter value of the optimized working surface mining height, the working surface propelling speed is set to be a common fixed value, and then a Sigmoid function between the length of the working surface and the damage ratio of the overlying strata is determined. After the optimized working surface mining height and working surface length are obtained, the working surface mining height and the working surface length can be set to the optimized fixed parameter values, and then a Sigmoid function between the working surface propulsion speed and the overburden damage ratio is determined.
(10) And respectively calculating the working face mining height M 1 and the working face length L 1 corresponding to the maximum value of the damage ratio increasing speed of the overburden and the working face pushing speed v 1 corresponding to the maximum value of the damage ratio decreasing speed by utilizing MATABLE.
Still taking the fitting result obtained under the typical working background of the Shendong mining area as an example, using MATABLE to calculate and obtain the working face mining height corresponding to the maximum value of the damage ratio increasing speed as 7.8445 and the working face length as 309.8928, and simultaneously calculating and obtaining the working face pushing speed corresponding to the maximum value of the cover rock damage ratio decreasing speed as 11.9551.
(11) Working face propulsion speed maximum value is obtained according to the propulsion speed of the coal mining machineWherein t is the daily coal cutting time of the coal cutter, v' is the coal cutting speed of the coal cutter, and d is the cutting depth of the coal cutter.
Still taking the Shendong mining area as an example, the maximum value of the propelling speed of the coal mining machine is 20.58 through calculation.
(12) And finally, determining that the value range of the optimized working surface elevation is 0~M 1, the value range of the optimized working surface length is 0~L 1, and the reasonable range of the optimized working surface propelling speed is v 1～vmax.
Still taking the Shendong mining area as an example, according to the result, the value range of the working face height is finally determined to be 0-7.8445 m, the value range of the working face length is determined to be 0-309.8928 m, and the reasonable range of the working face propelling speed is determined to be 11.9551-20.58 m/d.
For high-strength intensive coal exploitation, the increase of working face exploitation height, length and propelling speed is beneficial to large-scale exploitation, reduces cost and improves efficiency. As shown in fig. 4-6, the overburden damage ratio is positively correlated with the working face elevation and length and negatively correlated with the thrust rate. Therefore, it is preferable to control the damage ratio increase rate to be lower than the maximum value or to control the decrease rate to be higher than the maximum value on the premise of ensuring the production efficiency, so that the working face height and the working face length corresponding to the maximum value of the damage ratio increase rate of the overburden are set as the corresponding upper limit values, and the working face advancing rate corresponding to the maximum value of the damage ratio decrease rate of the overburden is set as the corresponding lower limit values.
The mining parameter optimization method mainly comprises the steps of firstly constructing a damage constitutive model, carrying out secondary development on logarithmic simulation software, secondly providing a overburden damage ratio model for quantitatively analyzing overburden damage degree, and thirdly calculating overburden damage degrees under different mining heights, different working surface lengths and different working surface propelling speeds through a numerical simulation method.
The method of the embodiment adopts Sigmoid function fitting to respectively obtain fitting formulas between the overburden damage ratio and different working surface heights, different working surface lengths and different working surface propelling speeds. The optimized value range of the working surface mining height and the working surface length is smaller than or equal to the corresponding parameter value when the cover rock damage ratio increases at the highest speed. The optimized value range of the working face advancing speed is the corresponding parameter value corresponding to the fastest reduction of the overburden damage ratio and the limit of the coal cutting speed of the coal cutter.
Aiming at severe damage of overburden in shallow high-strength mining of a western mining area, the invention provides a high-strength damage-reducing mining parameter optimization method based on the damage ratio of the overburden to the earth surface, which can ensure the mining benefit, greatly reduce the damage degree of the overburden caused by shallow high-strength mining from the source, reduce the damage of the earth surface from the source and is beneficial to realizing green mining in the ecologically fragile area of the western mining area.
Example III
The embodiment provides a mining method comprising the step of mining the mineral deposits in the overburden formation according to the mining parameters optimized by the mining parameter optimization method as described above.
Example IV
The present embodiment 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.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows.
Storage media, including both permanent and non-permanent, removable and non-removable media, may be implemented in any method or technology for storage of information. The information may be computer readable instructions, data structures, modules of a program, 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 disc read only memory (CD-ROM), digital Versatile Discs (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.
Example five
The present embodiment provides a computer device comprising a processor and a storage medium storing program code which, when executed by the processor, implements the steps of the mining parameter optimization method as described above.
In one embodiment, a computer device includes one or more processors (CPUs), an input/output interface, a network interface, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, etc., such as Read Only Memory (ROM) or flash memory (FLASH FLASH RAM). Memory is an example of computer-readable media.
It is noted that the terms used herein are used merely to describe particular embodiments and are not intended to limit exemplary embodiments in accordance with the present application, when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
It should be understood that the exemplary embodiments in this specification may be embodied in many different forms and should not be construed as limited to only the embodiments set forth herein. These embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of these exemplary embodiments to those skilled in the art, and should not be construed as limiting the application.
Claims (10)
1. A method for optimizing mining parameters, comprising the steps of:
acquiring the volume of a overburden stratum, acquiring the volume of a rock mass damaged in the overburden stratum caused by mining a mineral deposit under the overburden stratum according to different parameter value combinations of a plurality of mining parameters, and taking the ratio of the volume of the rock mass damaged to the volume of the overburden stratum as the overburden damage ratio of the overburden stratum;
for each mining parameter, establishing a relation model of the relation between the overburden damage ratio of the overburden stratum and the mining parameter according to the corresponding relation between the overburden damage ratio and different parameter values of the mining parameter, wherein the plurality of mining parameters comprise working face mining height, working face length and working face propelling speed;
the relation model about the relation between the overburden damage ratio and the working face mining height of the overburden stratum is as follows:
Wherein, Representing the overburden damage ratio of the overburden stratum, M represents the working face mining height, M1 represents the parameter value of the working face mining height corresponding to the maximum value of the increasing speed of the overburden damage ratio of the overburden stratum, and a, b and c are fitting constants;
the relation model about the relation between the overburden damage ratio and the working surface length of the overburden stratum is:
Wherein, The method comprises the steps of representing the overburden damage ratio of a overburden stratum, wherein L represents the length of a working surface, L1 represents the parameter value of the length of the working surface corresponding to the maximum value of the increasing speed of the overburden damage ratio of the overburden stratum, and a1, b1 and c1 are fitting constants;
the relation model about the relation between the overburden damage ratio and the working surface advancing speed of the overburden stratum is as follows:
Wherein, The method comprises the steps of representing the overburden damage ratio of a overburden stratum, v representing the advancing speed of a working surface, v1 representing the parameter value of the advancing speed of the working surface corresponding to the maximum value of the reducing speed of the overburden damage ratio of the overburden stratum, and a2, b2 and c2 being fitting constants;
And determining the maximum value of the change speed of the overburden damage ratio according to the relation model, taking the parameter value of the mining parameter corresponding to the maximum value of the change speed as the end point value of the parameter value range after optimizing the mining parameter, and further determining the end point value of the parameter value range after optimizing each of the plurality of mining parameters.
2. The production parameter optimization method of claim 1, wherein obtaining a damaged rock volume in the overburden formation resulting from production of the mineral reserve under the overburden formation in accordance with different parameter value combinations of the plurality of production parameters comprises:
dividing a overlying strata into a plurality of unit bodies, and acquiring actual shear stress born 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 which can be born by the unit body according to the strain tensor;
Comparing the actual shear stress born by the unit body with the theoretical shear stress capable of being born by the unit body, and determining whether the unit body is damaged or not according to a comparison result;
And obtaining the volume of each damaged unit body, summing the volumes of each damaged unit body, and taking the sum value as the volume of the damaged rock body in the overburden stratum.
3. The mining parameter optimization method of claim 2, wherein determining theoretical shear stress that the unit body can withstand based on the strain tensor comprises:
determining theoretical normal stress which can be born by the unit body according to the strain tensor of the elastic matrix in the unit body by using a preset damage constitutive model;
the theoretical shear stress that each unit can withstand is determined from the theoretical normal stress that the unit can withstand using the molar coulomb criterion.
4. The mining parameter optimization method of claim 3, wherein the impairment constitutive model is:
Wherein, G is the shear modulus of the overburden formation, which is the theoretical normal stress that the unit is capable of withstanding,For the shear modulus of the elastomeric matrix in the cell, eijkl is the elastic constant tensor of the elastomeric matrix in the cell, εkl is the strain tensor of the elastomeric matrix in the cell, emmkl is the elastic modulus of the elastomeric matrix in the cell,Is a Kronecker symbol;
The molar coulomb criterion is:
wherein c is the cohesive force between the unit bodies, phi is the internal friction angle between the unit bodies, and sigma n and tau are the theoretical normal stress and the theoretical shear stress which can be born by the unit bodies respectively.
5. The production parameter optimization method of claim 1, wherein establishing a relational model of a relationship between the overburden damage ratio and the production parameter for the overburden formation comprises:
And establishing a relation model about the relation between the overburden damage ratio of the overburden stratum and the exploitation parameters by using a Sigmoid function fitting method.
6. The method of optimizing production parameters of claim 1,
Determining a maximum value of the change speed of the overburden damage ratio according to the relation model, taking a parameter value of the mining parameter corresponding to the maximum value of the change speed as an endpoint value of a parameter value range after optimizing the mining parameter, and comprising the following steps:
when the relation model is a relation model related to the relation between the overburden damage ratio of the overburden stratum and the working face mining height, determining the maximum value of the increasing speed of the overburden damage ratio according to the relation model, and taking the parameter value of the working face mining height corresponding to the maximum value of the increasing speed as the upper limit value of the parameter value range after optimizing the working face mining height;
When the relation model is a relation model of relation between the overburden damage ratio of the overburden stratum and the length of the working surface, determining the maximum value of the increasing speed of the overburden damage ratio according to the relation model, and taking the parameter value of the length of the working surface corresponding to the maximum value of the increasing speed as the upper limit value of the parameter value range after optimizing the length of the working surface;
When the relation model is a relation model of relation between the overburden damage ratio and the working face advancing speed of the overburden stratum, determining a maximum value of the decreasing speed of the overburden damage ratio according to the relation model, and taking a parameter value of the working face advancing speed corresponding to the maximum value of the decreasing speed as a lower limit value of a parameter value range after optimizing the working face advancing speed.
7. The production parameter optimization method of claim 1, wherein prior to obtaining the volume of the overburden formation and obtaining the volume of rock mass damaged in the overburden formation as a result of producing the mineral reserve under the overburden formation in different parameter value combinations of the plurality of production parameters, the method further comprises the steps of:
And constructing a virtual model of the overburden stratum after mining the mineral deposit under the overburden stratum according to different parameter value combinations of a plurality of mining parameters by using a finite element analysis method.
8. A mining method comprising the steps of:
Mining of a mineral deposit in a overburden formation according to the mining parameters optimized by the mining parameter optimization method as claimed in any one of claims 1 to 7.
9. A storage medium storing a computer program which, when executed by a processor, implements the steps of the production parameter optimization method of any one of claims 1 to 7.
10. A computer device comprising a processor and a storage medium storing program code which, when executed by the processor, implements the steps of the production parameter optimization method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110064628.8A CN114818399B (en) | 2021-01-18 | 2021-01-18 | Mining parameter optimization method, mining method, storage medium and computer device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110064628.8A CN114818399B (en) | 2021-01-18 | 2021-01-18 | Mining parameter optimization method, mining method, storage medium and computer device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114818399A CN114818399A (en) | 2022-07-29 |
| CN114818399B true CN114818399B (en) | 2025-04-29 |
Family
ID=82524784
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110064628.8A Active CN114818399B (en) | 2021-01-18 | 2021-01-18 | Mining parameter optimization method, mining method, storage medium and computer device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114818399B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108104833A (en) * | 2017-12-21 | 2018-06-01 | 山东科技大学 | A kind of underworkings Bolting Parameters design method based on Central Point Method |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8768672B2 (en) * | 2007-08-24 | 2014-07-01 | ExxonMobil. Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
| CN104361211B (en) * | 2014-10-24 | 2018-07-27 | 中冶长天国际工程有限责任公司 | The structure and methods for using them of Statistical Damage Constitutive Model For Rock |
| US10546072B2 (en) * | 2016-03-28 | 2020-01-28 | Baker Huges, A Ge Company, Llc | Obtaining micro- and macro-rock properties with a calibrated rock deformation simulation |
| CN105927217B (en) * | 2016-04-29 | 2019-06-25 | 中国矿业大学 | A kind of nearly total rock up-protective layer mining Design method in seam mining |
| CN107942380A (en) * | 2017-10-23 | 2018-04-20 | 太原理工大学 | A kind of numerical simulation method for considering caving zone goaf |
| CN109577982B (en) * | 2018-12-12 | 2020-09-04 | 中国矿业大学 | Wall-type continuous mining, continuous charging and water-preserving coal mining, water resource migration monitoring, and water disaster warning methods |
| CN109681272B (en) * | 2018-12-24 | 2020-04-28 | 江西理工大学 | Discrimination method of overlying rock instability mutation in goaf of metal mine under the support of cemented ore pillar |
| CN110096787A (en) * | 2019-04-26 | 2019-08-06 | 辽宁工程技术大学 | A kind of optimization method of Waste Dump inside Open Pit Mine casting parameter |
| CN111577280B (en) * | 2020-05-19 | 2021-04-27 | 中国矿业大学 | Cooperative protective mining method of multi-level water resources in mining area |
-
2021
- 2021-01-18 CN CN202110064628.8A patent/CN114818399B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108104833A (en) * | 2017-12-21 | 2018-06-01 | 山东科技大学 | A kind of underworkings Bolting Parameters design method based on Central Point Method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114818399A (en) | 2022-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7137889B2 (en) | Prediction method of the lower limit of oil pool endowment depth of carbonatite in the deep-ultra-deep ocean facies | |
| Yang et al. | Improved PLS and PSO methods-based back analysis for elastic modulus of dam | |
| CA2375412C (en) | An improved simulation method and apparatus | |
| CN103091727B (en) | Based on the stress field DDM inversion method of measured result and linear programming technique | |
| CN108842821B (en) | Calculation method for reasonable buried depth of submarine tunnel constructed by drilling and blasting method | |
| CN106570287A (en) | Method for predicting water inflow of tunnel based on three-dimensional discrete fracture network | |
| CN114818493B (en) | A method for quantitative evaluation of tunnel rock mass integrity | |
| CN117473851A (en) | A method for predicting the floating process of shield tunnel based on random forest | |
| CN113223116A (en) | Landslide geological model based on AutoCAD, construction method and evaluation method | |
| CN111475604A (en) | Data processing method and device | |
| Ma et al. | Prediction of water inflow from fault by particle swarm optimization-based modified grey models | |
| CN117973191B (en) | Surrounding rock comprehensive degradation degree real-time evaluation method based on while-drilling test | |
| CN116307730A (en) | A construction risk assessment method, recording medium and system for a deep buried tunnel | |
| AU2021103225A4 (en) | Embankment dam analysis to reduce flood impact using sensor modelling and gis techniques | |
| CN112883661B (en) | A fracturing simulation method for soft and low-permeability oil and gas reservoirs | |
| CN114818399B (en) | Mining parameter optimization method, mining method, storage medium and computer device | |
| CN117709069A (en) | Three-dimensional caving numerical simulation method for metal mine | |
| CN109902360A (en) | Method, device and machine equipment for optimizing engineering parameters in drilling site operations | |
| CN119712219B (en) | Tunnel surrounding rock stability evaluation method and system based on optical fiber monitoring | |
| Qi et al. | A fast local nonlinear solution technique based on the partitioned finite element and interface element method | |
| CN111767323B (en) | Method for predicting length of stratum fracture in real time based on kernel SMO algorithm | |
| CN118520544A (en) | Deformation prediction method and equipment for adjacent existing tunnels caused by foundation pit construction based on BO-NGBoost-SHAP | |
| CN115788444A (en) | Tunnel intelligent construction system and method based on drilling and blasting method and storage medium | |
| Ariyan et al. | Seismic response of earth dams considering dynamic properties of unsaturated zone | |
| Shi et al. | Back analysis of mechanical parameters based on GPSO-BP neural network and its application |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |