Optimization method and system for identifying matching degree between backfill and surrounding rock under complex stress disturbance

Abstract

The present invention relates to the field of filling mining engineering, and provides a method and system for optimizing the matching degree of filling body and surrounding rock under complex stress disturbance, the method comprising: monitoring the vibration of the mining site, determining the stress disturbance signal generated by blasting; collecting rock samples, cutting holes in the rock samples to prepare hole models, loading stress disturbance signals into the hole rock samples, and determining the release energy of the hole model under stress disturbance; preparing filling body samples; performing compressive strength tests on the filling body samples to determine the elastic deformation energy of the filling body samples; and determining the matching coefficient between elastic deformation energy and release energy. It is used to solve the defects in the prior art that the filling body is prone to damage, fracture, collapse and other accidents under stress disturbance, and the safety hazards of mining are too great. The matching degree of the filling body and the surrounding rock can be judged under stress disturbance conditions, and the filling ratio can be optimized based on this to improve the safety of mining.

Description

Technical Field

[0001] The invention relates to the technical field of filling mining, and in particular to a method and system for identifying and optimizing the matching degree between a filling body and surrounding rock under complex stress disturbance. Background Art

[0002] Backfill mining technology plays a significant role in controlling ground pressure in mining areas, preventing deep well rock bursts, preventing mining deformation, improving ore recovery rates, and protecting the ecological environment. It is increasingly widely used in metal mines.

[0003] In the process of interaction between the filling body and the surrounding rock, the filling body mainly plays three roles: stress isolation, energy absorption and surrounding rock support. That is, after mining, the surrounding rock unloads energy to release deformation, compresses the filling body in the goaf, and accumulates deformation energy, resulting in energy exchange between the two. However, there is a problem of matching the filling body and the surrounding rock when using filling mining technology. Therefore, it is urgent to study the matching problem in the process of interaction between the surrounding rock and the filling body, so that the filling body can absorb more energy generated by the deformation and destruction of the surrounding rock, so as to achieve the purpose of maintaining the stability of the underground mining field. Scientifically determining the reasonable matching of the filling body and the surrounding rock is a key scientific problem that needs to be solved urgently.

[0004] At present, when studying the matching relationship between backfill proportion and surrounding rock, the usual practice is to specially adjust the backfill material, including ash-sand ratio, tailings grading, concentration, cementitious material ratio, curing time, etc., to ensure the needs of the project. However, it is difficult to meet the actual engineering requirements by only indoor waste rock cementing backfill proportion test. In order to overcome the defects of indoor mechanical tests, some scholars have adopted numerical simulation technology to establish a stope structure model of the combined structure of backfill and surrounding rock, and continuously change the strength of the backfill to determine the optimal ratio parameters of the backfill. Numerical simulation technology plays a huge advantage in dealing with the interaction and matching determination of backfill-surrounding rock at the stope scale. However, numerical models have certain defects in dealing with the heterogeneous properties of engineering rock mass, especially for granular bulk materials such as backfill. Numerical methods cannot yet deal well with the heterogeneity and anisotropy of model materials. At present, whether it is indoor mechanical test methods or numerical calculation methods, when dealing with the matching problem of backfill and surrounding rock, most of the focus is on the static mechanics of the interaction between backfill and surrounding rock. However, the actual mining process is often accompanied by complex stress disturbances such as excavation, mechanical disturbance, blasting, and rock burst.

[0005] Under the action of stress disturbance, the filling body is very likely to suffer from damage, fracture, collapse and other instability phenomena, thus failing to achieve the purpose of stress isolation and supporting the deformation and destruction of the surrounding rock, posing a safety hazard to mining. Summary of the invention

[0006] The present invention provides a method and system for optimizing the matching degree judgment between backfill and surrounding rock under complex stress disturbance, so as to solve the defects in the prior art that the backfill is prone to damage, fracture, collapse and other accidents under stress disturbance, and the potential safety hazards of mining are too great. The matching degree between backfill and surrounding rock can be judged under stress disturbance conditions, and the matching degree can be optimized based on this, thereby improving the safety of mining.

[0007] The present invention provides a method for identifying and optimizing the matching degree between a filling body and surrounding rock under complex stress disturbance, comprising:

[0008] Monitor the blasting stress wave and mining stress in the mining area;

[0009] Collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on blasting stress waves and mining stresses, the form of the stress disturbance loads is the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load. The hole model is used to simulate the mining field;

[0010] The filling body samples are prepared by on-site sampling or indoor material ratio test. The filling body samples are mixed with tailings sand, cementitious materials and water in different proportions;

[0011] Conduct a compressive strength test on the filling body sample, and determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0012] The matching coefficient between the elastic deformation energy and the release energy is determined. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared, and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

[0013] According to the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, the blasting stress wave and mining stress of the stope are monitored, including:

[0014] Set up seismic measuring equipment near the blasting source in the stope to monitor the blasting stress wave;

[0015] Mining stress monitoring equipment is installed in the pillars in the mining area to monitor the size, direction and affected area of ​​the mining stress formed after the surrounding rock is subjected to stress disturbance.

[0016] According to the method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, the blasting stress wave is monitored, including:

[0017] The blasting stress wave is monitored to obtain a waveform curve of the blasting stress wave, which includes peak stress, frequency spectrum, magnitude, particle vibration velocity, vibration displacement and vibration acceleration parameters.

[0018] According to the method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, rock samples are collected, holes are cut in the rock samples to prepare hole models, including:

[0019] According to the stope structure, calculate the volume ratio between holes and rocks in the stope, and determine the volume ratio between surrounding rock and filling body;

[0020] Rock samples are collected from the stope and holes are cut in the rock samples based on the volume ratio of surrounding rock to backfill.

[0021] According to the method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, a stress disturbance load is applied to the hole model based on the blasting stress wave and the mining stress. The stress disturbance load is in the form of the waveform of the blasting stress wave, and the amplitude is the value of the mining stress. The release energy of the hole model under the stress disturbance load is determined, including:

[0022] Load blasting stress wave and mining stress to the hole model until the hole model is destroyed, and record the strain curve of the hole model;

[0023] The release energy of the pore model is determined based on the strain curve.

[0024] According to the method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, a compressive strength test is performed on a backfill sample, and the elastic deformation energy of the backfill sample is determined by the peak stress and peak strength, including:

[0025] Construct damage constitutive model of filling body;

[0026] Static loading is applied to the filling body specimen to determine the peak stress and peak strength of the filling body specimen. The peak stress and peak strength are substituted into the damage constitutive model to determine the elastic deformation energy of the filling body.

[0027] The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention also includes:

[0028] If it is determined that the matching coefficient is greater than or equal to the set threshold, it is determined that the matching degree between the surrounding rock and the filling body is high;

[0029] The working status of the filling body is continuously monitored, and if it is determined that the working status of the filling body is abnormal, an abnormal alarm is issued.

[0030] According to the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention, the matching coefficient between elastic deformation energy and release energy is determined, including:

[0031] Determine the ratio of elastic deformation energy to release energy. If the ratio is greater than or equal to 1, it is determined that the matching degree between the surrounding rock and the filling body is high.

[0032] If the ratio is less than 1, it is determined that the surrounding rock and the filling body do not match.

[0033] The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the present invention also includes:

[0034] The filling ratio is optimized with the matching coefficient, filling cost and filling body strength as the objective function. The factors considered include the filling body ash-sand ratio, concentration and curing time. A third-order polynomial response surface model is established through the three-objective response surface method to optimize the design of the filling body samples.

[0035] The present invention also provides a system for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance, comprising:

[0036] Signal acquisition module, used to monitor blasting stress waves and mining stress in the mining area;

[0037] The first determination module is used to collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on blasting stress waves and mining stresses, the form of the stress disturbance loads is the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load. The hole model is used to simulate underground stopes in metallurgical mines;

[0038] Model building module, used for field sampling or indoor material proportioning test to prepare filling body samples, the filling body samples are made of tailings sand, cementitious materials and water in different proportions;

[0039] The second determination module is used to perform a compressive strength test on the filling body sample, determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0040] The optimization module is used to determine the matching coefficient between elastic deformation energy and release energy. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared and the matching coefficient between elastic deformation energy and release energy is re-determined.

[0041] The present invention also provides an electronic 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, it implements any of the above-mentioned methods for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance.

[0042] The present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the method for optimizing the matching degree between the filling body and the surrounding rock under any of the complex stress disturbances described above is implemented.

[0043] The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements any of the above-mentioned methods for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance.

[0044] In the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided in the present application, the blasting stress wave and mining stress can be determined based on the actual mining site, and then a hole model simulating the surrounding rock and a backfill sample can be constructed respectively. The released energy generated by rock blasting can be calculated by applying the collected blasting stress wave and mining stress to the hole model. The elastic deformation energy of the backfill sample can be determined by performing a compressive strength test on the backfill sample. Subsequently, the matching coefficient between the elastic deformation energy and the release energy can be calculated, which is equivalent to comparing the size relationship between the elastic deformation energy and the release energy. If the matching coefficient is lower than the set threshold, it means that the backfill cannot completely absorb the energy released by rock blasting. In this way, the scheme of the present application can identify the matching degree between the backfill and the surrounding rock under stress disturbance conditions, and optimize the matching degree based on this, thereby improving the safety of mine mining. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0046] Figure 1 It is one of the flow diagrams of the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by an embodiment of the present invention;

[0047] Figure 2 This is the second flow chart of the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by an embodiment of the present invention;

[0048] Figure 3 It is a structural schematic diagram of a system for identifying and optimizing the matching degree between a filling body and surrounding rock under complex stress disturbance provided by an embodiment of the present invention;

[0049] Figure 4 is a schematic diagram of a stress-strain curve of a filling body provided in an embodiment of the present invention;

[0050] Figure 5It is one of the performance evaluation schematic diagrams provided by the embodiment of the present invention;

[0051] Figure 6 This is the second performance evaluation schematic diagram provided by the embodiment of the present invention;

[0052] Figure 7 This is the third performance evaluation schematic diagram provided by the embodiment of the present invention;

[0053] Figure 8 It is a schematic diagram of the physical structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0054] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0055] With the depletion of shallow mineral resources in my country, most mines are facing the technical difficulties of ground pressure, rock burst and large deformation in deep mining. The interaction between filling body and surrounding rock involves dynamic problems. At present, when studying the matching relationship between filling body ratio and surrounding rock, the usual practice is to make special adjustments to the filling materials, including ash-sand ratio, tailings grading, concentration, cementitious material ratio, curing time, etc., to ensure the needs of the project. However, it is difficult to meet the actual engineering requirements only through indoor waste rock cementing filling ratio tests. In order to overcome the defects of indoor mechanical tests, some scholars have adopted numerical simulation technology to establish a stope structure model of the combined structure of filling body and surrounding rock, and continuously change the filling body strength to determine the optimal filling body ratio parameters. Numerical simulation technology plays a huge advantage in dealing with the interaction and matching determination of filling body and surrounding rock at the stope scale. However, numerical models have certain defects in dealing with the heterogeneous properties of engineering rock mass, especially for granular bulk materials such as filling body. Numerical methods cannot yet deal with the heterogeneity and anisotropy of model materials well. At present, whether it is indoor mechanical test methods or numerical calculation methods, when dealing with the matching problem of backfill and surrounding rock, most of the focus is on the static mechanics of the interaction between backfill and surrounding rock. However, the actual mining process is often accompanied by complex stress disturbances such as excavation, mechanical disturbance, blasting, and rock burst.

[0056] The destruction process of the interaction between the filling body and the surrounding rock under stress disturbance is extremely complex, which also brings challenges to the determination and research of the filling body-surrounding rock matching relationship. When the pillar is mined, a shock wave is first generated in the rock body. The shock wave decays quickly in the rock body and gradually decays into a stress wave. At the same time, the action range of the far-field stress wave is large, and the general distance exceeds 120 to 150 times the charge radius. Under the action of stress disturbance, the filling body is very prone to damage, fracture, collapse and other instability phenomena, which can not play the purpose of stress isolation and supporting the deformation and destruction of the surrounding rock, bringing safety hazards to mining. Therefore, it is necessary to study the method of judging and determining the matching degree between the filling body and the surrounding rock under the condition of stress disturbance dynamic load.

[0057] Based on this, the solution of the present application provides a method and system for identifying and optimizing the matching degree between the backfill and the surrounding rock under complex stress disturbance, which can identify the matching degree between the backfill and the surrounding rock under stress disturbance conditions, and optimize the matching degree based on this, thereby improving the safety of mining.

[0058] Figure 1 It is one of the flow charts of the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by an embodiment of the present invention.

[0059] Figure 2 This is the second flow chart of the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by an embodiment of the present invention.

[0060] like Figure 1 and Figure 2 As shown, this embodiment provides a method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance, including:

[0061] Step 101, monitoring the blasting stress wave and mining stress of the stope;

[0062] Step 102, collecting rock samples, cutting holes in the rock samples to prepare a hole model, applying a stress disturbance load to the hole model based on the blasting stress wave and the mining stress, the stress disturbance load is in the form of the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determining the release energy of the hole model under the stress disturbance load, the hole model is used to simulate the underground stope of a metallurgical mine;

[0063] Step 103, sampling on site or preparing a filling body sample by a material ratio test indoors, wherein the filling body sample is made by mixing tailings sand, a cementitious material and water in different proportions;

[0064] Step 104, performing a compressive strength test on the filling body sample, determining the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determining the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0065] Step 105, determining the matching coefficient between the elastic deformation energy and the release energy. If the matching coefficient is lower than a set threshold, re-preparing the filling body sample and re-determining the matching coefficient between the elastic deformation energy and the release energy.

[0066] In practical applications, the stope may specifically include an underground stope in a metallurgical mine.

[0067] In practice, the released energy refers to the total energy released when the hole model ruptures. Since the hole model is used to simulate the surrounding rock after blasting, the total energy calculated when the hole model ruptures can be used to characterize the total energy released when the surrounding rock is destroyed.

[0068] In the implementation, the backfill used can be a tailings cemented backfill. The compressive strength of the tailings cemented backfill is approximately linear before its peak stress, and its stress-strain relationship curve is approximately linear. After exceeding the peak stress, its strength decreases sharply and its bearing capacity is seriously reduced. Therefore, when the compressive strength of the tailings cemented backfill exceeds its peak stress, its internal structure is completely damaged, showing macroscopic damage. Based on this, in this embodiment, the prepared backfill sample can be tested for compressive strength, and the maximum elastic deformation energy of the backfill can be determined. The elastic deformation energy can characterize the maximum stress intensity that the backfill can absorb, that is, the compressive strength of the backfill.

[0069] In practical applications, if it is determined that the matching coefficient between deformation energy and release energy is lower than a set threshold, the filling body can be re-prepared and the matching coefficient can be re-determined until the matching coefficient is higher than or equal to the set threshold, and the re-prepared filling body can be used to replace the original unqualified filling body to improve the safety of the mining process.

[0070] In the method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided in the present embodiment, the blasting stress wave and mining stress can be determined based on the actual mining site, and then a hole model simulating the surrounding rock and a backfill sample can be constructed respectively. The release energy generated by rock blasting can be calculated by applying the collected blasting stress wave and mining stress to the hole model. The elastic deformation energy of the backfill sample can be determined by performing a compressive strength test on the backfill sample. Subsequently, the matching coefficient between the elastic deformation energy and the release energy can be calculated, which is equivalent to comparing the size relationship between the elastic deformation energy and the release energy. If the matching coefficient is lower than the set threshold, it means that the backfill cannot completely absorb the energy released by rock blasting. In this way, the scheme of the present application can identify the matching degree between the backfill and surrounding rock under stress disturbance conditions, and optimize the matching degree based on this, thereby improving the safety of mine mining.

[0071] In an exemplary embodiment, monitoring of blasting stress waves and mining stress in a mining area includes:

[0072] Set up seismic measuring equipment near the blasting source in the stope to monitor the blasting stress wave;

[0073] Mining stress monitoring equipment is installed in the pillars in the mining area to monitor the size, direction and affected area of ​​the mining stress formed after the surrounding rock is subjected to stress disturbance.

[0074] In implementation, the seismic measuring device in this embodiment can be specifically arranged near the blasting source of the mining site, and the seismic measuring device used can be a blasting seismometer, which is used to monitor the blasting stress wave.

[0075] By monitoring the blasting stress wave, the waveform curve of the blasting stress wave can be obtained. The waveform curve includes parameters such as peak stress, frequency spectrum, magnitude, particle vibration velocity, vibration displacement and vibration acceleration.

[0076] The stress monitoring device in this embodiment can be installed at the mining pillar position, and can be used to monitor the activity area, value, direction and frequency of the stress generated during the ore mining process.

[0077] In practical applications, the frequency value of the applied stress disturbance load can be determined by performing wavelet transform on the waveform data of the blasting stress wave, and the loading amplitude of the stress disturbance load can be the maximum value of the mining stress obtained by monitoring.

[0078] The mining pillar generally refers to the ore rock structure supporting the roof of the goaf in the mine. During the ore mining process, the goaf formed during the mining process is mainly maintained by the stability of the surrounding rock itself and the supporting capacity of the pillar. Therefore, the stability of the mining pillar is crucial to the safety of ore mining. In the scheme of this embodiment, by collecting the blasting stress wave of the blasting source and the mining stress at the mining pillar, it can better reflect the blasting mining scenario, provide a reasonable basis for the strength design of the filling body, improve the safety and stability of the filling mining field, and thus improve the safety of rock mining in the entire mine.

[0079] In an exemplary embodiment, a rock sample is collected, and holes are cut into the rock sample to prepare a hole model, including:

[0080] According to the stope structure, calculate the volume ratio between holes and rocks in the stope, and determine the volume ratio between surrounding rock and filling body;

[0081] Rock samples are collected from the stope and holes are cut in the rock samples based on the volume ratio of surrounding rock to backfill.

[0082] In practical applications, the volume ratio of the goaf and pillars in a mine can be calculated, and then the volume ratio of the surrounding rock to the filling body can be determined. In this way, holes can be cut on the collected rock samples to simulate the goaf of the mine. The size of the cut holes can be determined based on the volume ratio of the surrounding rock to the filling body determined in the above process.

[0083] In an exemplary embodiment, collecting a rock sample from a stope includes:

[0084] Rock samples are collected around the explosion source in the mining area, and standard rock samples are obtained through drilling, cutting, grinding and processing.

[0085] In this embodiment, rock samples are collected around the explosion source in the mining area and processed to obtain rock samples, which can reflect the degree of mutual matching between the surrounding rock and the filling body at the surrounding position of the filling mining area, so that the filling ratio parameters can be accurately judged.

[0086] In an exemplary embodiment, a stress disturbance load is applied to the hole model based on the blasting stress wave and the mining stress, the stress disturbance load is in the form of the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and the release energy of the hole model under the stress disturbance load is determined, including:

[0087] Load blasting stress wave and mining stress to the hole model until the hole model is destroyed, and record the strain curve of the hole model;

[0088] The release energy of the pore model is determined based on the strain curve.

[0089] In practical applications, recording the strain curve may specifically include the following steps:

[0090] The compressive stress test of the hole model under the action of complex stress disturbance was carried out, and the rupture events were recorded by acoustic emission technology during the test. The loading path is to input the blasting seismic wave collected on site to load the hole model until it is unstable and damaged; the blasting seismic wave can also be simplified to fatigue load, and the loading frequency is obtained by wavelet transformation of the blasting seismic wave through waveform data processing, and the loading amplitude adopts the maximum value of the mining stress monitored on site. The hole model is loaded until the sample is damaged, and the strain curve of the sample is recorded.

[0091] In implementation, the release energy of the hole model is determined based on the strain curve, which can be specifically determined by integrating the stress-strain curve. The specific formula is as follows:

[0092]

[0093]

[0094]

[0095] In the above formula, U is the total strain energy, U e is the elastic strain energy, n is the number of cycles; the superscript e represents the elastic strain, is the elastic strain at the ith cycle, is the elastic strain at the i+1th cycle; U Rd is the dissipated energy density, that is, the energy released per unit rock mass, where σ i and σ i+1 is the stress at the i-th and i+1-th rings, ε i and ε i+1 is the corresponding deformation at the i-th and i+1-th rings, ε max is the peak strain of the porous rock.

[0096] In an exemplary embodiment, a compressive strength test is performed on a filling body sample, and elastic deformation energy of the filling body sample is determined from peak stress and peak strength, including:

[0097] Construct damage constitutive model of filling body;

[0098] Static loading is applied to the filling body specimen to determine the peak stress and peak strength of the filling body specimen. The peak stress and peak strength are substituted into the damage constitutive model to determine the elastic deformation energy of the filling body.

[0099] In practical applications, the filling body can be regarded as an isotropic continuous medium. According to the strain equivalence principle, it can be obtained that:

[0100] σ=Eε(1-D)(4)

[0101] Where E is the elastic modulus of the filling body, σ is the stress, ε is the strain, and D is the damage factor.

[0102] According to the relationship between the peak stress before the filling failure and the damage factor, the relationship between the two is determined as follows:

[0103] D=A*ε β (5)

[0104] Where A and β are unknown coefficients.

[0105] Combining the above strain equivalence principle and the relationship between stress and damage factor, the constitutive equation of the filling body before the peak value is obtained:

[0106] σ=Eε-EAε β +1(6)

[0107] Then the peak strain ε p By introducing formula (6), we can get the peak stress σ c , at this time the differential corresponding to the stress is 0, thus obtaining the following formula:

[0108] β=σp / (Eε p -σ p )(7)

[0109]

[0110] According to the stress-strain curve of the filling body during compression and the characteristics of compaction and consolidation of the filling body, it is believed that the total energy input from the outside before the peak is completely converted into elastic deformation energy U. Therefore, the calculation formula of the absorbed elastic deformation energy is:

[0111]

[0112] When the peak intensity σ p When , the deformation energy density absorbed by the filling body under pressure is:

[0113]

[0114] In an exemplary embodiment, the stress-strain curve of the filling body can also be obtained by performing indoor tests, and then the deformation energy of the filling body can be obtained by integrating the curve.

[0115] In an exemplary embodiment, it further includes:

[0116] If it is determined that the matching coefficient is greater than or equal to the set threshold, it is determined that the matching degree between the surrounding rock and the filling body is high;

[0117] The working status of the filling body is continuously monitored, and if it is determined that the working status of the filling body is abnormal, an abnormal alarm is issued.

[0118] In practical applications, based on the above formula (3) and formula (10), the matching coefficient calculation formula can be determined as:

[0119]

[0120] Where V CTB is the volume of the filling in the stope, V R is the surrounding rock volume.

[0121] In practical applications, when the optimal ratio is used for mine filling, the long-term service status of the filling body can be inspected through drilling observation, seismic CT and other means to ensure that the filling body can exert its energy absorption, stress isolation and collapse support functions.

[0122] In an exemplary embodiment, determining a matching coefficient between elastic deformation energy and release energy includes:

[0123] Determine the ratio of elastic deformation energy to release energy. If the ratio is greater than or equal to 1, it is determined that the matching degree between the surrounding rock and the filling body is high.

[0124] If the ratio is less than 1, it is determined that the surrounding rock and the filling body do not match.

[0125] In an exemplary embodiment, the method for optimizing the matching degree between the filling body and the surrounding rock under complex stress disturbance further includes:

[0126] The filling ratio is optimized with the matching coefficient K, filling cost and filling body strength as the objective function. The factors considered include the filling body ash-sand ratio, concentration and curing time. The third-order polynomial response surface model is established through the three-objective response surface method to optimize the design of the filling body samples.

[0127] In practical applications, the established third-order polynomial response surface model is as follows:

[0128]

[0129] Among them, y is the target value of each response, x j is the lime-sand ratio, x j is the concentration, x k For maintenance time, and The parameters x are i The lower and upper limits of the range, β0, β i , β ii , β ij , β iij and β ijk is the coefficient to be determined.

[0130] After obtaining the optimal combination parameters of different response values ​​through the response surface method, the same group of parameters with similar different response values ​​are obtained through searching, which is the optimal combination ratio of the filling body.

[0131] The following is a description of the system for identifying and optimizing the matching degree of backfill and surrounding rock under complex stress disturbance provided by the present invention. The system for identifying and optimizing the matching degree of backfill and surrounding rock under complex stress disturbance described below and the method for identifying and optimizing the matching degree of backfill and surrounding rock under complex stress disturbance described above can be referenced to each other.

[0132] Figure 3 It is a structural schematic diagram of a system for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by an embodiment of the present invention.

[0133] like Figure 3 As shown, the matching degree identification and optimization system of filling body and surrounding rock under complex stress disturbance provided by this embodiment includes:

[0134] The signal acquisition module 301 is used to monitor the blasting stress wave and mining stress in the mining area;

[0135] The first determination module 302 is used to collect rock samples, cut holes in the rock samples to prepare a hole model, apply a stress disturbance load to the hole model based on the blasting stress wave and the mining stress, the stress disturbance load is in the form of the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load, and the hole model is used to simulate the underground stope of the metallurgical mine;

[0136] Model construction model 303 is used for field sampling or indoor preparation of filling body samples by material proportioning test. The filling body samples are mixed with tailings sand, cementitious materials and water in different proportions;

[0137] The second determination module 304 is used to perform a compressive strength test on the filling body sample, determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0138] The optimization module 305 is used to determine the matching coefficient between the elastic deformation energy and the release energy. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

[0139] In an exemplary embodiment, the signal acquisition module is further used for:

[0140] Set up seismic measuring equipment near the blasting source in the stope to monitor the blasting stress wave;

[0141] Mining stress monitoring equipment is installed in the pillars in the mining area to monitor the size, direction and affected area of ​​the mining stress formed after the surrounding rock is subjected to stress disturbance.

[0142] In an exemplary embodiment, the first determining module is further configured to:

[0143] According to the stope structure, calculate the volume ratio between holes and rocks in the stope, and determine the volume ratio between surrounding rock and filling body;

[0144] Rock samples are collected from the stope and holes are cut in the rock samples based on the volume ratio of surrounding rock to backfill.

[0145] In an exemplary embodiment, the first determining module is further configured to:

[0146] Load blasting stress wave and mining stress to the hole model until the hole model is destroyed, and record the strain curve of the hole model;

[0147] The release energy of the pore model is determined based on the strain curve.

[0148] In an exemplary embodiment, the second determining module is further configured to:

[0149] Construct damage constitutive model of filling body;

[0150] Static loading is applied to the filling body specimen to determine the peak stress and peak strength of the filling body specimen. The peak stress and peak strength are substituted into the damage constitutive model to determine the elastic deformation energy of the filling body.

[0151] In an exemplary embodiment, a monitoring module is also included, and the monitoring module is specifically used to:

[0152] If it is determined that the matching coefficient is greater than or equal to the set threshold, it is determined that the matching degree between the surrounding rock and the filling body is high;

[0153] The working status of the filling body is continuously monitored, and if it is determined that the working status of the filling body is abnormal, an abnormal alarm is issued.

[0154] In an exemplary embodiment, the monitoring module is specifically used for:

[0155] Determine the ratio of deformation energy to release energy. If the ratio is greater than or equal to 1, it is determined that the matching degree between the surrounding rock and the filling body is high.

[0156] If the ratio is less than 1, it is determined that the surrounding rock and the filling body do not match.

[0157] The specific implementation method of the system for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided in this embodiment can be implemented with reference to the above-mentioned embodiments, which will not be described in detail here.

[0158] The solution of this application also provides a specific application scenario. Taking Huibaoling Iron Mine as an example, filling body and surrounding rock samples were obtained underground. The surrounding rock was taken at five levels with a burial depth of -760m, -690m, -620m, -550m and -480m, and processed into hole models with a hole diameter of 35mm. The filling body samples were obtained by underground drilling. Filling body mechanical tests and rock fatigue tests were carried out indoors.

[0159] The filling body strength is 3.231MPa, the peak strain is 1.036%, the elastic modulus is 645.36MPa, and the stress-strain curve of the filling body is as follows: Figure 4 As shown, Figure 4 is a schematic diagram of a stress-strain curve of a filling body provided in an embodiment of the present invention. Figure 4 It can be determined that the elastic energy of the filling body is 16.123KJ.

[0160] Figure 5 This is one of the performance evaluation schematic diagrams provided in the embodiment of the present invention.

[0161] Figure 6 This is the second performance evaluation schematic diagram provided by the embodiment of the present invention.

[0162] Figure 5The stress-strain curves of the cave rock at four levels of -760m, -690m, -620m and -550m under a confining pressure of 15MPa are shown in Figure 1. Figure 6 This is the calculated release energy curve. The calculated release energies are 14.506 kJ, 11.625 kJ, 7.669 kJ, 5.606 kJ and 3.652 kJ. The calculation shows that the matching coefficient is greater than 1, which indicates that the surrounding rock and the filling body are well matched.

[0163] The solution of the present application also provides another specific embodiment. Specifically, the five levels of -760m, -690m, -620m, -550m and -480m are also used. When the designed filling concentration is 70% and the proportions are 1:8, 1:10, 1:12, 1:15 and 1:20 respectively, the matching coefficient relationship curve between the tailings cemented filling body and the surrounding rock at different burial depths is obtained from the above formula, that is, Figure 7 , Figure 7 This is the third performance evaluation schematic diagram provided by the embodiment of the present invention, such as Figure 7 As shown in the figure, when the matching coefficient is greater than 1, the matching of the surrounding rock and the filling body meets the requirements. Figure 7 The matching coefficients of other burial depths can be calculated to determine the filling ratio.

[0164] Figure 8 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 8 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the method for identifying and optimizing the matching degree between the filling body and the surrounding rock under complex stress disturbance, and the method includes:

[0165] Monitor the blasting stress wave and mining stress in the mining area;

[0166] Collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on blasting stress waves and mining stresses, the form of the stress disturbance loads is the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load. The hole model is used to simulate underground stopes in metallurgical mines;

[0167] The filling body samples are prepared by on-site sampling or indoor material ratio test. The filling body samples are mixed with tailings sand, cementitious materials and water in different proportions;

[0168] Conduct a compressive strength test on the filling body sample, and determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0169] The matching coefficient between the elastic deformation energy and the release energy is determined. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared, and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

[0170] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, RandomAccess Memory), disk or optical disk, etc. Various media that can store program codes.

[0171] On the other hand, the present invention further provides a computer program product, which includes a computer program. The computer program can be stored in a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the method for identifying and optimizing the matching degree between the backfill and the surrounding rock under complex stress disturbance provided by the above methods. The method includes:

[0172] Monitor the blasting stress wave and mining stress in the mining area;

[0173] Collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on blasting stress waves and mining stresses, the form of the stress disturbance loads is the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load. The hole model is used to simulate underground stopes in metallurgical mines;

[0174] The filling body samples are prepared by on-site sampling or indoor material ratio test. The filling body samples are mixed with tailings sand, cementitious materials and water in different proportions;

[0175] Conduct a compressive strength test on the filling body sample, and determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0176] The matching coefficient between the elastic deformation energy and the release energy is determined. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared, and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

[0177] In another aspect, the present invention further provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to execute the method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance provided by the above methods, the method comprising:

[0178] Monitor the blasting stress wave and mining stress in the mining area;

[0179] Collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on blasting stress waves and mining stresses, the form of the stress disturbance loads is the waveform of the blasting stress wave, and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load. The hole model is used to simulate underground stopes in metallurgical mines;

[0180] The filling body samples are prepared by on-site sampling or indoor material ratio test. The filling body samples are mixed with tailings sand, cementitious materials and water in different proportions;

[0181] Conduct a compressive strength test on the filling body sample, and determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve;

[0182] The matching coefficient between the elastic deformation energy and the release energy is determined. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared, and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

[0183] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, i.e., they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those of ordinary skill in the art may understand and implement it without creative effort.

[0184] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods of each embodiment or some parts of the embodiment.

[0185] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for optimizing the matching degree between backfill and surrounding rock under complex stress disturbance, characterized in that: include: Monitor the blasting stress wave and mining stress in the mining area; Collect rock samples, cut holes in the rock samples to prepare hole models, apply stress disturbance loads to the hole models based on the blasting stress wave and the mining stress, wherein the stress disturbance loads are in the form of the waveform of the blasting stress wave and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load, wherein the hole model is used to simulate the stope; The filling body samples are prepared by sampling on site or by material proportion test indoors, wherein the filling body samples are mixed with tailings sand, cementitious materials and water in different proportions; Performing a compressive strength test on the filling body sample, determining the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determining the elastic deformation energy of the filling body sample from the axial compression stress-strain curve; determining a matching coefficient between the elastic deformation energy and the release energy, and if the matching coefficient is lower than a set threshold, re-preparing the filling body sample, and re-determining the matching coefficient between the elastic deformation energy and the release energy; The collecting of rock samples and cutting holes in the rock samples to prepare a hole model comprises: According to the stope structure, the volume ratio between the holes and the rocks in the stope is calculated, and the volume ratio between the surrounding rock and the filling body is determined; collecting a rock sample from the stope, and cutting holes in the rock sample based on a volume ratio of the surrounding rock to the backfill; The step of applying a stress disturbance load to the hole model based on the blasting stress wave and the mining stress, wherein the stress disturbance load is in the form of a waveform of the blasting stress wave and an amplitude is a value of the mining stress, and determining the release energy of the hole model under the stress disturbance load comprises: Loading the blasting stress wave and the mining stress to the hole model until the hole model is destroyed, and recording the strain curve of the hole model; determining a release energy of the pore model based on the strain curve; The release energy of the hole model is determined based on the strain curve by integrating the stress-strain curve. The specific formula is as follows: （1） （2） （3） Among them, in the above formula U is the total strain energy, U e is the elastic strain energy, n is the number of cycles; the superscript e represents the elastic strain, is the elastic strain at the ith cycle, is the elastic strain at the i+1th cycle; is the dissipated energy density, that is, the energy released per unit rock mass, where and is the stress at the i-th and i+1-th rings, and is the corresponding deformation at the i-th and i+1-th rings, is the peak strain of the porous rock.

2. The method for optimizing matching degree between backfill and surrounding rock under complex stress disturbance according to claim 1 is characterized in that: The monitoring of blasting stress waves and mining stress in the stope includes: Setting up seismic measuring equipment near the blasting source of the stope to monitor the blasting stress wave; Mining stress monitoring equipment is arranged in the ore pillars in the stope to monitor the size, direction and affected area of ​​the mining stress formed after the surrounding rock is subjected to stress disturbance.

3. The method for optimizing matching degree between backfill and surrounding rock under complex stress disturbance according to claim 2 is characterized in that: The monitoring of the blasting stress wave comprises: The blasting stress wave is monitored to obtain a waveform curve of the blasting stress wave, wherein the waveform curve includes peak stress, frequency spectrum, magnitude, particle vibration velocity, vibration displacement, and vibration acceleration parameters.

4. The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance according to claim 1 is characterized in that: The method of performing a compressive strength test on the filling body sample and determining the elastic deformation energy of the filling body sample according to the peak stress and the peak strength includes: constructing a damage constitutive model of the filling body; A static loading load is applied to the filling body sample to determine the peak stress and peak strength of the filling body sample, and the peak stress and peak strength are substituted into the damage constitutive model to determine the elastic deformation energy of the filling body.

5. The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance according to any one of claims 1 to 4, characterized in that: Also includes: If it is determined that the matching coefficient is greater than or equal to the set threshold, it is determined that the matching degree between the surrounding rock and the filling body is high; The working state of the filling body is continuously monitored, and if it is determined that the working state of the filling body is abnormal, an abnormal alarm is issued.

6. The method for optimizing matching degree between backfill and surrounding rock under complex stress disturbance according to claim 5 is characterized in that: The determining of the matching coefficient between the elastic deformation energy and the release energy comprises: Determining a ratio of the elastic deformation energy to the release energy, and if the ratio is greater than or equal to 1, determining that the matching degree between the surrounding rock and the filling body is high; If the ratio is less than 1, it is determined that the surrounding rock does not match the filling body.

7. The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance according to any one of claims 1 to 4, characterized in that: Also includes: The filling ratio is optimized with the matching coefficient, filling cost and filling body strength as the objective function. The factors considered include the filling body ash-sand ratio, concentration and curing time. A third-order polynomial response surface model is established through the three-objective response surface method to optimize the design of the filling body samples.

8. The matching degree identification and optimization system between backfill and surrounding rock under complex stress disturbance is characterized by: The method for identifying and optimizing the matching degree between backfill and surrounding rock under complex stress disturbance as described in any one of claims 1 to 7 comprises: Signal acquisition module, used to monitor blasting stress waves and mining stress in the mining area; A first determination module is used to collect rock samples, cut holes in the rock samples to prepare a hole model, apply a stress disturbance load to the hole model based on the blasting stress wave and the mining stress, wherein the stress disturbance load is in the form of a waveform of the blasting stress wave and the amplitude is the value of the mining stress, and determine the release energy of the hole model under the stress disturbance load, wherein the hole model is used to simulate an underground stope of a metallurgical mine; A model building module is used for on-site sampling or indoor material proportioning test to prepare filling body samples, wherein the filling body samples are mixed with tailings sand, cementitious materials and water in different proportions; A second determination module is used to perform a compressive strength test on the filling body sample, determine the elastic deformation energy of the filling body sample from the peak stress and peak strength, or determine the elastic deformation energy of the filling body sample from the axial compression stress-strain curve; The optimization module is used to determine the matching coefficient between the elastic deformation energy and the release energy. If the matching coefficient is lower than a set threshold, the filling body sample is re-prepared and the matching coefficient between the elastic deformation energy and the release energy is re-determined.

Images