Filling mining method stope structure parameter optimization method and system based on numerical simulation

Abstract

The invention relates to a filling mining method stope structure parameter optimization method and system based on numerical simulation, and the method comprises the steps: obtaining mine engineering geological conditions, building a mine entity model, and determining a stope structure parameter scheme; establishing a stope model according to the determined stope structure parameters; combining the mine entity model and the stope model, re-grouping the combined model, and assigning values to parameters in the combined model by setting mechanical parameters extracted from initial conditions and mine engineering geological conditions; and carrying out excavation and filling simulation on the combined model, respectively determining the stability indexes of the stope in the states of the original rock, the goaf and the filling body, and screening out an optimal stope structure parameter scheme. The problem that the resource development efficiency is affected due to the fact that an existing stope structure parameter optimization scheme is poor in analysis efficiency and accuracy in mine stability analysis under the multi-stage, multi-time-sequence and multi-scale conditions is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of mining, and specifically to an optimization method and system for stope structure parameters of backfill mining method based on numerical simulation. Background Technique

[0002] The statements in this part only provide background technical information related to the present invention, and do not necessarily constitute prior art.

[0003] Stope structure parameters are parameters used to describe various sizes and configurations of the stope, including the height, width, length of the stope, the width of the room pillar, the height of the sill pillar and the crown pillar, etc., and are important parameters in mining engineering.

[0004] Currently, for the research on stope structure parameters, most rely on engineering analogy method or static theoretical analysis. However, in the complex geological and topographical conditions of mines with multiple time sequences and multiple scales, the efficiency and accuracy of analyzing the rationality of stope structure parameters are lacking, and it cannot be well combined with the occurrence conditions of actual mines. Therefore, it is difficult to select a stope structure parameter scheme that takes into account both economic benefits and safety and stability, resulting in resource waste caused by conservative mining to ensure safety in some mining areas or potential accident hazards caused by ignoring safety to improve economic benefits. Summary of the Invention

[0005] In order to solve the technical problems existing in the above background technique, the present invention provides an optimization method and system for stope structure parameters of backfill mining method based on numerical simulation, which solves the problems that the existing stope structure parameter optimization scheme lacks efficiency and accuracy in the analysis of mine stability under multiple stages, multiple time sequences and multiple scales, and further affects the resource development efficiency.

[0006] In order to achieve the above object, the present invention adopts the following technical solutions:

[0007] The first aspect of the present invention provides an optimization method for stope structure parameters of backfill mining method based on numerical simulation, including the following steps:

[0008] Obtain the engineering geological conditions of the mine, establish a mine entity model, and determine the stope structure parameter scheme;

[0009] According to the determined stope structure parameters, establish a stope model;

[0010] Merge the mine entity model and the stope model, re-group the merged model, and assign values to the parameters in the merged model by setting the initial conditions and the mechanical parameters extracted from the engineering geological conditions of the mine;

[0011] Carry out excavation and filling simulation on the merged model, determine the stability indexes of the stope under the conditions of the original rock, the goaf and the filling body respectively, and screen out the best stope structure parameter scheme.

[0012] As a further implementation method, a mine entity model is established, and a stope structure parameter scheme is determined. Specifically: By obtaining mine engineering geological condition data and on-site investigation, a mine entity model is established. Combining the Mathews graphical method and the engineering analogy method, multiple groups of stope structure parameter schemes are confirmed.

[0013] As a further implementation method, combining the Mathews graphical method and the engineering analogy method, multiple groups of stope structure parameter schemes are confirmed. Specifically: By obtaining mine engineering geological condition data, the parameters required for calculating the stability coefficient N are selected, the stability coefficient of a specific stope is calculated, and referring to the Mathews diagram, the limit shape coefficient S within the stable area is delineated; Combining with the calculation formula of the shape coefficient S, through the engineering analogy method, referring to the stope structure parameter delineation of mines under similar conditions, the safe value ranges of the stope width and length are calculated, and multiple groups of stope structure parameter schemes are formulated within the value range of the stable area.

[0014] As a further implementation method, the mine entity model and the stope model are merged. Specifically: In the numerical simulation software, through the File-grid menu, the Rhino-Griddle plug-in is imported to generate a file in the.f3grid format. After ensuring that the coordinate systems of the two models are consistent and the grid element types are compatible, the stope model is placed into the mine entity model through the zone merge command.

[0015] As a further implementation method, the merged model is regrouped, and the parameters in the merged model are assigned by setting the initial conditions and the mechanical parameters extracted from the mine engineering geological conditions. Specifically: By selecting the elastic model, defining the rock mass material parameters, the model boundary conditions, the self-weight stress and density, and initializing the displacement, the parameters in the merged model are assigned.

[0016] As a further implementation method, excavation and filling simulations are performed on the merged model to determine the stability indicators of the stope under the conditions of the original rock, the goaf, and the filling body respectively. Specifically: According to the zone command combined with the regional definition parameters, the area to be excavated is determined. After excavation, the goaf is filled, a filling area is created, and the filling material is assigned. The defined filling material is assigned to the created filling area.

[0017] As a further implementation method, the best stope structure parameter scheme is screened out. Specifically: Using "maximum principal stress", "minimum principal stress", "plastic zone distribution", and "displacement change amount" as stability indicators, and using "mining ratio", "stope production capacity", "ore loss rate", and "mining cost" as influencing factors of technical and economic indicators, the best stope structure parameter scheme is calculated using the Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS).

[0018] The second aspect of the present invention provides a system required to implement the above method, including:

[0019] A data acquisition module, configured to: obtain the mine engineering geological conditions;

[0020] A model establishment module, configured to: establish a mine entity model, determine a stope structure parameter scheme; establish a stope model according to the determined stope structure parameters;

[0021] A model merging module, configured to: merge the mine entity model and the stope model, regroup the merged model, and assign values to the parameters in the merged model by setting initial conditions and mechanical parameters extracted from the mine engineering geological conditions;

[0022] A parameter optimization module, configured to: perform excavation and filling simulations on the merged model, determine the stability indicators of the stope under the conditions of the original rock, goaf, and filling body respectively, and screen out the best stope structure parameter scheme.

[0023] The third aspect of the present invention provides a computer-readable storage medium.

[0024] A computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, it implements the steps in the above-mentioned method for optimizing stope structure parameters of backfill mining method based on numerical simulation.

[0025] The fourth aspect of the present invention provides a computer device.

[0026] A computer device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, it implements the steps in the above-mentioned method for optimizing stope structure parameters of backfill mining method based on numerical simulation.

[0027] Compared with the prior art, the above one or more technical solutions have the following beneficial effects:

[0028] 1. By combining the Mathews graphical method and the engineering analogy method, it is possible to ensure the stability of the scheme while conforming to the conventional engineering practice, more efficiently and quickly delineate multiple groups of stope structure parameter schemes, and provide a reliable range for scheme optimization.

[0029] 2. By using numerical simulation technology, comprehensively consider the stope structure stability of the stope mined and managed by the backfill mining method in multiple stages, multiple time sequences, and multiple scales such as the original rock, excavated goaf, and filling, which more comprehensively reflects the adaptability of the stope structure parameters to the stope compared with the general static theoretical analysis, ensures the safety of the whole mining process, improves the accuracy and calculation efficiency of the analysis results, and saves the trial-and-error cost.

[0030] 3. By using the approximation ideal method, multiple indicators such as economy, safety, and technology are comprehensively analyzed. Selecting the best from multiple options can not only improve the safety of mine operations but also ensure technical feasibility on the premise of ensuring safety, and improve economic benefits as much as possible; achieving the improvement of the safety and economic benefits of resource recovery and utilization. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] The accompanying drawings forming a part of this invention are used to provide a further understanding of the invention. The schematic embodiments and descriptions thereof of the invention are used to explain the invention and do not constitute an improper limitation to the invention.

[0032] Figure 1 It is a schematic diagram of the framework of the optimization method for stope structure parameters of the filling mining method based on numerical simulation provided by one or more embodiments of the invention;

[0033] Figure 2 It is a schematic diagram of the optimization process of the stope structure parameters of the filling mining method based on numerical simulation provided by one or more embodiments of the invention;

[0034] Figure 3 It is a schematic diagram of the correlation between the Mathews stability coefficient and the shape factor provided by one or more embodiments of the invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0035] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0036] It should be noted that the following detailed descriptions are all exemplary and are intended to provide further explanations of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those of ordinary skill in the technical field to which the present invention belongs.

[0037] Term Explanation:

[0038] The filling mining method refers to the process of backfilling the mined - out area with filling materials (such as tailings, waste rock, cement, etc.) during mining to support the surrounding rock and reduce surface subsidence.

[0039] Rhino8, a 3D modeling software. It is famous for its powerful NURBS (Non - Uniform Rational B - Spline) modeling ability and is widely used in fields such as industrial design, architectural design, jewelry design, ship design, automotive design, 3D printing, animation production, etc.

[0040] FLAC3D, a three - dimensional geotechnical numerical simulation software, mainly used to simulate the mechanical behavior of geotechnical bodies, rocks, and other continuous media under static and dynamic conditions. FLAC3D is widely used in fields such as geotechnical engineering, mining engineering, geological engineering, and underground engineering.

[0041] In this solution, the modeling software involved is not limited to Rhino8, and the numerical simulation software is not limited to FLAC3D. The following embodiments only take Rhino8 and FLAC3D as examples to introduce the specific process.

[0042] In order to solve the technical problems existing in the above-mentioned background technology, the following embodiments provide an optimization method and system for stope structure parameters of filling mining method based on numerical simulation, which solves the problems that the existing optimization schemes for stope structure parameters lack efficiency and accuracy in the analysis of mine stability in multiple stages, multiple time sequences, and multiple scales, thus affecting the resource development efficiency.

[0043] Embodiment 1:

[0044] As Figures 1 - 2 shown, the optimization method for stope structure parameters of filling mining method based on numerical simulation includes the following steps:

[0045] S1: Obtain the mine engineering geological conditions and the basic drawings required for numerical simulation, and determine the modeling research scope;

[0046] S2: Use the mine model range and basic drawings obtained in S1 to establish a mine entity model in Rhino8;

[0047] S3: By analyzing the data related to the occurrence conditions of the ore body obtained in S1, combined with the Mathews graphical method and the engineering analogy method, calculate and draw up multiple sets of stope structure parameter schemes;

[0048] S4: Use the FLAC3D built-in modeling extrusion function to establish stope models with different stope structure parameters set in S3 respectively;

[0049] S5: Use the FLAC3D command stream to merge the mine entity model obtained in S2 with the stope model established in S4;

[0050] S6: Re-group the merged grid model in the model of FLAC3D, set the initial conditions of the entity model through the built-in command stream of the FLAC3D numerical software, and assign values to the model parameters using the mechanical parameters obtained in S1;

[0051] S7: Use the FLAC3D internal command stream to simulate the excavation and filling of the ore body in the stope, and obtain the stability indexes of the stope under the conditions of the original rock, goaf, and filling body respectively;

[0052] S8: Combine the influencing factors such as the stability indexes, technical and economic indexes obtained in S1 and S7, and introduce the Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) to select the best stope structure parameter scheme.

[0053] Among them, S1: Obtain the basic drawings required for the mine engineering geological conditions and numerical simulation, and determine the scope of modeling research; specifically including:

[0054] Through the collection of mine basic data and on-site investigation, obtain engineering geological conditions such as the stress distribution in the mine and the rock mechanics parameters of different strata, and obtain basic model drawings such as borehole data and stratigraphic maps, ore body cross-section maps in the mining area: Simplify the basic drawings appropriately in the auto-cad software for the geometric drawings of the mining area such as the mine topographic map and the ore body cross-section map, and delineate the ore body model range, and save it as a dxf file.

[0055] S2: Use the mine model range and basic drawings obtained in S1 to establish a mine entity model in Rhino8; specifically including:

[0056] Import the cad terrain file (contour data) obtained in S1 into Rhino8 (professional modeling software) for surface modeling, and import the pre-ore body and stratigraphic basic drawings into Rhino8 to generate entities. After merging all the curves and surfaces in the model into a whole, use the GSurf and GVol functions of the external Griddle plugin in Rhino8 to perform mesh division on the entity model. After the model generalization is completed, save the entity model as a.f3grid format.

[0057] S3: By analyzing the relevant data on the occurrence conditions of the ore body obtained in S1, and combining the Mathews graphical method and the engineering analogy method, draw up multiple groups of stope structure parameter schemes; specifically including:

[0058] The stope span is one of the important factors in the stope structure parameters. The Mathews graphical method can simply and quickly judge the stability of the mining area. Analyze the relevant data on the occurrence conditions of the ore body obtained in S1, select the parameters required to calculate the stability coefficient N, calculate the stability coefficient of a specific stope, and circle the limit shape coefficient S in the stable area according to the Mathews diagram. Combining with the calculation formula of the shape coefficient S, through the engineering analogy method, consult and refer to the stope structure parameter delineation of mines with similar conditions, calculate the safe value range of the stope width and length, and draw up multiple groups of stope structure parameter schemes within the stable area value range.

[0059] The main calculation steps of the Mathews graphical method are as follows: Comprehensively consider the rock mass quality, stress state, the orientation of the stope exposed surface, etc. for mine development and tunneling engineering, initially determine the rock mass stability index N, calculate the stope exposed area shape coefficient S value, and according to the N value and S, use the predicted stable area, potentially unstable area, and caving area divided in the stability chart to preliminarily judge the stability of the stope. The specific calculation formula is as follows:

[0060] N = Q′·A·B·C;

[0061] Where: N is the Mathews stability coefficient; Q′ is the modified Q value; A is the rock stress coefficient; B is the joint orientation coefficient; C is the gravity adjustment coefficient.

[0062] In this embodiment, the correlation between the Mathews stability coefficient and the hydraulic radius is as Figure 3 shown. Regarding the Q value, the Mathews stability diagram method uses the modified NGI tunnel quality index Q′. Different from the Q value, the stress reduction factor (SRF) and the joint water inflow reduction factor (J w ) in the Q′ value are both 1.0, making the Mathews stability diagram method have a wider applicability in different hydrogeological and in-situ stress environments. It can be applied to various types of underground projects, such as stope, tunnel, etc., without having to readjust complex parameters for each specific situation. At the same time, the calculation is relatively simplified and more convenient for rapid application in practical engineering. Therefore, the calculation formula for the Q′ value is:

[0063]

[0064] Where: RQD is the rock mass quality index, J n is the number of joint sets, J r is the joint roughness coefficient, J a is the joint alteration coefficient.

[0065] The A value takes into account the influence of high stress on reducing the rock mass stability. A is the ratio of the uniaxial compressive strength σ c value of the rock at the boundary of the evaluated goaf face to the induced compressive stress σ i . The relational expression for the A value is:

[0066]

[0067] The B value is obtained from Table 1 according to the angle α between the dip angle of the main joint set and the dip angle of the exposed face.

[0068] Table 1 Corresponding relationship between B value and angle α

[0069] included angle α 90° 60° 40° 20° 0° value of B 1.0 0.8 0.4 0.2 0.3

[0070] Calculation formula for the stope exposed face orientation correction coefficient C value:

[0071] C = 8 - 7cosα;

[0072] Where: α is the angle between the exposed face and the horizontal plane.

[0073] According to the value range of the shape coefficient S of the stope in a stable state corresponding to the calculated N value, the shape coefficient S can be defined as the ratio of the area to the perimeter, and the relational expression is as follows:

[0074]

[0075] In the formula: X is the stope width, and Y is the stope length. When Y / X > 4, the shape coefficient S basically remains unchanged, that is, it is considered that the stability of the exposed surface is mainly controlled by the exposed surface width at this time. By calculation, the safe value ranges of the stope width and length are obtained. Through the engineering analogy method, the stope structure parameter delineation situations of mines under similar conditions are consulted and referred to, and multiple groups of stope structure parameter schemes are drawn up within the value range calculated by the Mathews graphical method for subsequent optimization.

[0076] S4: Use the FLAC3D built-in modeling extrusion function to establish geometric models with different stope structure parameters selected in S3 respectively. After mesh division, further generate solid models.

[0077] S5: Use the FLAC3D command stream to merge the mine solid model obtained in S2 with the stope model established in S4; specifically including:

[0078] In FLAC3D, import the.f3grid format file generated by the Rhino-Griddle plugin through the File-grid menu. After ensuring that the coordinate systems of the two models are consistent and the mesh element types are compatible, input the zone merge command to place the stope model into the mine solid model.

[0079] S6: Re-group the merged mesh model in the model of FLAC3D, set the initial conditions of the solid model through the built-in command stream of the FLAC3D numerical software, and use the mechanical parameters obtained in S1 to assign parameters to different groups of the model respectively; specifically including:

[0080] Select the elastic model;

[0081] zone cmodel mohr-coulomb;

[0082] Define the rock mass material parameters;

[0083] zone initialize density Assign the density value range group 'group name AA' union group 'group name O-A';

[0084] zone property bulk Assign the bulk modulus value shear Assign the shear modulus value cohesion Assign the corresponding value of the cohesion friction Assign the internal friction angle value range group 'group name BB' union group 'group name O-B';

[0085] Define the model boundary conditions;

[0086] zone gridpoint fix velocity-z range position-z boundary interval;

[0087] zone gridpoint fix velocity-x range position-x boundary interval;

[0088] zone gridpoint fix velocity-x range position-x boundary interval;

[0089] zone gridpoint fix velocity-y range position-y boundary interval;

[0090] zone gridpoint fix velocity-y range position-y boundary interval;

[0091] zone face apply stress-z apply stress assignment range position-z interval assignment position-x interval assignment position-y interval assignment;

[0092] Define self-weight stress and density;

[0093] model gravity 00-10;

[0094] zone property density density assignment;

[0095] Initialize displacement;

[0096] model initialize.

[0097] S7: Use the internal command stream of FLAC3D to simulate the excavation and filling of the ore body in the stope, and obtain the stability indicators of the stope under the conditions of the original rock, goaf, and filling body respectively. Specifically, it includes:

[0098] Run the model after initialization of S6 to obtain a simulation of the excavation process. First, define the excavation area: Use the "zone" command in combination with the area definition parameters to determine the area to be excavated. For example, "zone brick x-range 5 10 y-range 3 8 z-range 2 6" defines a cuboid area with x ranging from 5 to 10, y ranging from 3 to 8, and z ranging from 2 to 6 as the excavation area. Delete the excavation area: Use the "zone delete" command to delete the defined excavation area. For example, "zone delete id 101 - 200", assuming the element IDs of the excavation area are from 101 to 200.

[0099] After excavation is completed, backfill the mined - out area: First, create a backfill material model and set its properties. For example, "model material create elastic id 2 bulk 2.0e9 shear 1.2e9" creates an elastic backfill material and sets the bulk modulus and shear modulus. Subsequently, create the backfill area: Use the "zone" command to create an area with the same shape and position as the excavation area for backfilling. For example, "zone brick x-range 5 10 y-range 3 8 z-range 2 6". Assign the backfill material: Use the "zone property material 2" command to assign the defined backfill material to the created backfill area.

[0100] The stope stability indicators include: maximum principal stress, minimum principal stress, displacement, and plastic zone. After the simulation calculations for different steps are completed, they are respectively called and displayed using command streams:

[0101] Maximum principal stress: plot zone contour stress principal maximum;

[0102] Minimum principal stress: plot zone contour stress principal minimum;

[0103] Displacement contour map call: plot contour displacement

[0104] Plastic zone contour map call: plot zone color by state.

[0105] S8: Comprehensively select the influencing factors such as the safety indicators, technical and economic indicators of the mine obtained in S1 and S7 to achieve a more comprehensive evaluation of the proposed solutions. Evaluate and analyze different solutions in combination with the influencing factors. There are various evaluation methods, such as the analytic hierarchy process, grey relational method, fuzzy decision method, and technique for order preference by similarity to an ideal solution (TOPSIS), etc. TOPSIS can comprehensively consider the above multiple conflicting and restrictive factors. By establishing a multi-objective function, it can more comprehensively evaluate the advantages and disadvantages of different stope structure parameter solutions, which is more objective, easier to calculate, can make full use of data, handle fuzziness more objectively, has strong comprehensiveness and low cost compared with other methods. Therefore, introduce TOPSIS to select the best stope structure parameter solution; specifically include:

[0106] Comprehensively select four technical and economic indicators of stope production capacity, ore recovery rate, mining and cutting ratio, and mining cost, as well as four safety indicators of "maximum principal stress, minimum principal stress, displacement, and plastic zone" obtained from Flac3d numerical simulation. Among them, stope production capacity and ore recovery rate are benefit-type indicators, and the remaining 6 indicators are cost-type indicators. Use the vector normalization method to obtain the normalized decision matrix. The specific formula is as follows:

[0107] Benefit-type attribute

[0108] Cost-type attribute

[0109] In the formula: r ij is the element of the normalized decision matrix; m is the number of solutions; a ij is the j-th indicator of the i-th solution in the decision matrix.

[0110] Determine the positive ideal solution Z + and the negative ideal solution Z - , and calculate the distances from each solution to the ideal solution and the negative ideal solution. The specific formula is as follows:

[0111]

[0112] In the formula: is the distance from the i-th solution to the positive ideal solution; is the distance from the i-th solution to the negative ideal solution; Z + is the positive ideal solution, and the optimal value of each indicator (the maximum value for benefit-type indicators and the minimum value for cost-type indicators); Z - is the negative ideal solution, which is the worst value of each indicator (the minimum value for benefit-type indicators and the maximum value for cost-type indicators); Z ij is the value representing the j-th indicator of the i-th solution.

[0113] Determine the relative closeness C i , and the calculation formula is:

[0114]

[0115] According to C i Sort the advantages and disadvantages of each plan according to the size. The closer to 1, the closer the structural parameter plan is to the positive ideal solution and the better; the closer to 0, the closer the plan is to the negative ideal solution and the worse, and finally obtain the best plan.

[0116] In this embodiment, for the part of determining the model conditions and parameters: collect the basic data of the mine, conduct on-site process geological exploration, obtain the basic maps of the mining area, rock structural plane parameters, and stress distribution, and obtain the rock mechanical parameters through rock mechanics experiments and rock mass quality evaluations on rock samples, and initially establish a geometric model and a mechanical model.

[0117] In this embodiment, for the part of model establishment and numerical analysis: select the mining method in combination with the occurrence conditions of the ore body and the actual project, calculate the exposed area of the stope, and use the Mathews graphical method and engineering analogy method to draw up multiple groups of stope structural parameter plans, and simulate the stope operation through FLAC3D software to solve the stability of each plan;

[0118] In this embodiment, for the part of result analysis and plan optimization: select stability indexes such as "maximum principal stress", "minimum principal stress", "plastic zone distribution", "displacement change amount", and technical and economic index influencing factors such as "development ratio", "stope production capacity", "ore loss rate", "mining cost", etc., and introduce the technique for order preference by similarity to ideal solution to calculate the best stope structural plan.

[0119] In this embodiment, by combining the Mathews graphical method and the engineering analogy method, it can ensure that the proposed parameter range will not damage the stope stability on the basis of the theoretical calculation of stope stability. At the same time, the engineering analogy method further screens out the plans that meet the conventional engineering practice, and more efficiently and quickly delimits multiple groups of stope structural parameter plans, providing a reliable range for plan optimization.

[0120] This embodiment uses numerical simulation technology. When the mine adopts the backfill mining method for mining operations through Flac3d and other software simulation and restoration, the model is easy to adjust and has a high degree of restoration. It only needs to assign values ​​through command streams, apply boundary conditions and other functions to change the mining area's occurrence conditions, and has a wide range of applications. Through data evolution, the stability of mining sites in different regions in multiple stages, multiple time series and multiple scales, such as the original rock state, the empty area state after excavation, and the empty area filling state, can be analyzed. The study of the applicability of the mining site structural parameters to the mining site during the entire mining process is realized, which can more comprehensively identify potential risks than general static theoretical analysis and provide a reliable basis for scheme optimization; compared with the similar material simulation method, it is not affected by factors such as materials, sites, and time, and the trial and error cost is low; under low-cost test conditions, the safety of the entire mining process is ensured and the accuracy and calculation efficiency of the analysis results are improved, so as to achieve truly economical, efficient and accurate pre-mining preview analysis.

[0121] This embodiment combines the qualitative and quantitative approach to the ideal method to comprehensively consider the impact of multiple indicators related to mine production benefits such as economy, safety, technology, etc. on the selection of stope structure parameter schemes, thereby reducing the limitations of a single indicator and achieving the organic unity of economic benefits, social benefits, and environmental benefits; selecting the best from a variety of schemes, and maximizing the economic benefits while ensuring the safety of mine operations and technical feasibility; achieving a more comprehensive improvement in the safety and economic benefits of resource recycling from multiple perspectives, and providing technical support for the realization of resource-saving and environmentally friendly construction of mines.

[0122] Embodiment 2:

[0123] A system for implementing the above method comprises:

[0124] The data acquisition module is configured to: obtain mining engineering geological conditions;

[0125] The model building module is configured to: build a mine entity model and determine a stope structure parameter scheme; and build a stope model according to the determined stope structure parameters;

[0126] The model merging module is configured to: merge the mine entity model and the stope model, regroup the merged model, and assign values ​​to the parameters in the merged model by setting the initial conditions and mechanical parameters extracted from the mine engineering geological conditions;

[0127] The parameter optimization module is configured to: perform excavation and filling simulation on the merged model, determine the stability index of the stope under the original rock, void and filling state respectively, and select the best stope structure parameter scheme.

[0128] By combining the Mathews graphical method with the engineering analogy method, the scheme can be ensured to be stable while complying with conventional engineering practices. It can also more efficiently and quickly identify multiple groups of stope structure parameter schemes, providing a reliable range for scheme optimization.

[0129] Numerical simulation technology is used to comprehensively consider the stability of the mine structure in multiple stages, multiple time series and multiple scales such as original rock, excavation and filling in the mine, which is mined and managed by the backfill mining method. It reflects the adaptability of the mine structure parameters to the mine more comprehensively than the general static theoretical analysis, ensures the safety of the entire mining process, improves the accuracy and calculation efficiency of the analysis results, and saves the cost of trial and error.

[0130] By approximating the ideal method, a comprehensive analysis of multiple indicators such as economy, safety, and technology is conducted. Selecting the best among multiple options can not only improve the safety of mine operations, but also ensure the feasibility of technology while ensuring safety, thereby maximizing economic benefits and improving the safety and economic benefits of resource recycling.

[0131] Embodiment three:

[0132] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When the program is executed by a processor, the steps in the method for optimizing stope structure parameters of the backfill mining method based on numerical simulation as described in the first embodiment above are implemented.

[0133] Embodiment 4:

[0134] This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the method for optimizing the structural parameters of a filling mining method based on numerical simulation as described in the first embodiment above are implemented.

[0135] The steps or modules involved in the above embodiments 2 to 4 correspond to those in embodiment 1. For the specific implementation, please refer to the relevant description of embodiment 1. The term "computer-readable storage medium" should be understood as a single medium or multiple media including one or more instruction sets; it should also be understood to include any medium that can store, encode or carry an instruction set for execution by a processor and enable the processor to execute any method in the present invention.

[0136] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for optimizing stope structure parameters of backfill mining method based on numerical simulation, characterized in that: The following steps are involved: Obtain mining engineering geological conditions, establish mining entity model, and determine stope structure parameter plan; Establish a stope model based on the determined stope structural parameters; Merge the mine entity model and the stope model, regroup the merged model, and assign values ​​to the parameters in the merged model by setting the initial conditions and mechanical parameters extracted from the mine engineering geological conditions; The excavation and filling simulations were carried out on the merged model to determine the stability indicators of the stope under the original rock, void and filling states, and to screen out the optimal stope structure parameter scheme.

2. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 1, characterized in that: Establish a mine entity model and determine the stope structure parameter plan. Specifically, establish a mine entity model by obtaining mine engineering geological conditions data and on-site investigation, and combine the Mat hews graphical method with the engineering analogy method to confirm multiple sets of stope structure parameter plans.

3. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 2, characterized in that: Combining the Mathews graphical method with the engineering analogy method, multiple groups of stope structural parameter schemes are confirmed. Specifically, by obtaining the engineering geological conditions data of the mine, the parameters required to calculate the stability coefficient N are selected, the stability coefficient of the specific stope is calculated, and the limit shape coefficient S in the stable zone is delineated by comparing the Mathews diagram; combined with the calculation formula of the shape coefficient S, the stope structural parameter delineation of mines with similar conditions is consulted and referred to by the engineering analogy method, the safe value range of the stope width and length is calculated, and multiple groups of stope structural parameter schemes are formulated within the value range of the stable zone.

4. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 1, characterized in that: Merge the mine entity model and the stope model. Specifically, in the numerical simulation software, import the Rhino-Griddle plug-in through the File-grid menu to generate a .f3grid format file. After ensuring that the coordinate systems of the two models are consistent and the grid unit types are compatible, use the zone merge command to place the stope model into the mine entity model.

5. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 1, characterized in that: The merged models are regrouped, and the parameters in the merged model are assigned by setting the initial conditions and the mechanical parameters extracted from the mining engineering geological conditions. Specifically, the parameters in the merged model are assigned by selecting the elastic model, defining the rock material parameters, model boundary conditions, self-weight stress and density, and initializing the displacement.

6. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 1, characterized in that: The merged model is simulated for excavation and filling to determine the stability indicators of the mine in the original rock, void and filling states respectively. Specifically, the area to be excavated is determined according to the zone command combined with the area definition parameters. After the excavation is completed, the void is filled, a filling area is created, and the filling material is assigned. The defined filling material is assigned to the created filling area.

7. The method for optimizing stope structure parameters of backfill mining method based on numerical simulation according to claim 1, characterized in that: The optimal stope structure parameter scheme was screened out, specifically: taking "maximum principal stress", "minimum principal stress", "plastic zone distribution" and "displacement change" as stability indicators, and "mining-cutting ratio", "stope production capacity", "ore loss rate" and "mining cost" as technical and economic indicator influencing factors, and using the approximate ideal solution method to calculate the optimal stope structure parameter scheme.

8. The stope structure parameter optimization system of filling mining method based on numerical simulation is characterized by: The data acquisition module is configured to: obtain mining engineering geological conditions; The model building module is configured to: build a mine entity model and determine a stope structure parameter scheme; and build a stope model according to the determined stope structure parameters; The model merging module is configured to: merge the mine entity model and the stope model, regroup the merged model, and assign values ​​to the parameters in the merged model by setting the initial conditions and mechanical parameters extracted from the mine engineering geological conditions; The parameter optimization module is configured to: perform excavation and filling simulation on the merged model, determine the stability index of the stope under the original rock, void and filling state respectively, and select the best stope structure parameter scheme.

9. A computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the method for optimizing stope structure parameters of the backfill mining method based on numerical simulation as described in any one of claims 1 to 7.

10. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the method for optimizing the structural parameters of a filling mining method based on numerical simulation are implemented as described in any one of claims 1 to 7.

Images