CN111199108A - Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model - Google Patents
Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model Download PDFInfo
- Publication number
- CN111199108A CN111199108A CN202010014781.5A CN202010014781A CN111199108A CN 111199108 A CN111199108 A CN 111199108A CN 202010014781 A CN202010014781 A CN 202010014781A CN 111199108 A CN111199108 A CN 111199108A
- Authority
- CN
- China
- Prior art keywords
- blasting
- model
- numerical simulation
- pile
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000005422 blasting Methods 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004088 simulation Methods 0.000 title claims abstract description 21
- 230000011218 segmentation Effects 0.000 title claims abstract description 10
- 238000009826 distribution Methods 0.000 claims abstract description 28
- 239000002360 explosive Substances 0.000 claims description 13
- 238000010586 diagram Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000005520 cutting process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000002699 waste material Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a blasting blockiness prediction method of a gridding segmentation numerical simulation blasting model, which realizes blasting blockiness distribution prediction of the blasting heap by carrying out gridding segmentation on a three-dimensional blasting model and analyzing the blockiness distribution of each segment one by one. The method has the advantages of high precision, thoroughness, no damage, safety, intuition and the like, and has important significance for predicting the blasting block size.
Description
Technical Field
The invention relates to the technical field of engineering blasting, in particular to a blasting blockiness prediction method of a gridding segmentation numerical simulation blasting pile model.
Background
Blasting is one of four major processes of surface mining, the quality of the effect seriously affects the production efficiency and economic benefit of mines, and the blasting block size and the distribution thereof are important indexes for judging the quality of the blasting effect. Too large or too small blasting block degree can restrict mine production, causing cost waste.
Through development for many years, blasting block size prediction methods applying statistical analysis methods, stress wave theories, energy theories or combination modes thereof become mature day by day, but have the defects of large influence by working conditions, time and labor waste and large prediction error. However, the currently developed numerical simulation method is limited to the simulation of a single blast hole due to the limited technical level, but the single blast hole cannot represent the whole blast area due to the complexity of the blasting process, so that the prediction result of the method is limited greatly.
Therefore, how to efficiently, accurately and comprehensively predict the integral blasting lump size and the distribution thereof, the aim of guiding the on-site blasting design and construction is achieved, the problem to be solved urgently in blasting research is solved, and the method has important guiding significance for mine production.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a blasting block size prediction method for a gridding segmentation numerical simulation blasting pile model, which fills the blank of numerical simulation blasting block size prediction of an integral blasting area. The method comprises the steps of carrying out numerical simulation by using GDEM software to obtain a three-dimensional model of the blasting pile, cutting a section of the model by adopting a grid division mode, analyzing the blasting block degree of the section by using image processing and block degree statistics modes, and judging the whole block degree distribution of the blasting pile thoroughly, visually, clearly and accurately.
In order to solve the technical problems, the invention adopts the following technical scheme:
a method for predicting blasting blockiness of a gridding segmentation numerical simulation blasting pile model comprises the following steps:
1) performing numerical simulation by using GDEM software to obtain a three-dimensional explosive stack model;
2) dividing a grid of the three-dimensional explosive stack model to intercept any number of explosive stack transverse and longitudinal section diagrams;
3) performing image processing on the sectional diagram by using MATLAB, wherein the processing procedures comprise gray level image conversion, image enhancement, image segmentation and image edge detection;
4) determining the proportional size based on the three-dimensional model of the blasting pile, quantizing each block degree of the processed image to obtain the size of each block degree, and counting the block degree distribution of the whole profile to obtain the block degree distribution data of the profile;
5) repeating the steps 3) and 4) to obtain the block degree distribution of all the section diagrams, and obtaining the block degree distribution data of the blasting pile by an averaging method;
6) manually analyzing the spatial position distribution rule of each block degree interval;
7) comparing the block size distribution obtained in the steps 5) and 6) with a target, and determining a blasting optimization scheme;
8) and adjusting blasting parameters in GDEM software, and repeating the steps 1) to 7) until the target block size is reached.
Further, the step 1) specifically comprises: and setting model calculation parameters by using GDEM software according to blasting design parameters and geological data, and carrying out numerical simulation on the blasting area to obtain a three-dimensional blasting pile model.
Further, in the step 2), 3 transverse side-looking sections of the three-dimensional explosive stack model are taken from the explosive stack extending direction by dividing the grid pair, and the distances from the positions to the left side edge of the step are respectively 15m, 20m and 35 m; 3 vertical overlooking sections are taken in the vertical direction of the blasting pile, and the distances from the positions to the bottom of the step are respectively 4m, 8m and 10 m.
Compared with the prior art, the invention has the beneficial technical effects that:
1) the method is carried out on the basis of the three-dimensional model of the blasting pile obtained by numerical simulation by using GDEM software, the limitation of the conventional single-hole prediction is broken through, the prediction result is more accurate, and the method is actually fit;
2) the blasting block size is analyzed on the model cutting section by dividing the grids, the size of the grids can be adjusted at will according to requirements, and the internal condition of the blasting pile is comprehensively mastered;
3) compared with the well-known blasting block size prediction method, the method is more thorough, visual and clear, can quickly analyze the spatial position distribution rule of each block size interval, and has a guiding effect on field production.
Drawings
The invention is further illustrated in the following description with reference to the drawings.
FIG. 1 is a three-dimensional model diagram of a blasting pile obtained by the method of the invention;
FIG. 2 is a cross-sectional view of a transverse side elevation section 1 taken by the method of the present invention;
FIG. 3 is a block size diagram of the blasting pile after the section 1 is processed by the method of the invention;
FIG. 4 is the bulk distribution data of the blasting pile of the section 1 of the method of the present invention;
FIG. 5 is a plot of the bulk distribution of the detonation at section 1 of the method of the invention;
FIG. 6 is the overall blasting lump size distribution data of the method of the present invention.
Detailed Description
A method for predicting blasting block size of a gridding segmentation numerical simulation blasting pile model comprises the following steps:
1) setting model calculation parameters by using GDEM software according to blasting design parameters and geological data, and carrying out numerical simulation on a blasting area to obtain a three-dimensional blasting pile model, referring to figure 1;
2) and cutting a section of the three-dimensional explosive stack model by the grid division pair from two directions of a transverse extending direction (a side-looking direction) and a vertical direction (a top-looking direction). In the example, 3 transverse side-looking cross sections are taken in the extending direction of the blasting pile, and the distances from the positions to the left side edge of the step are respectively 15m, 20m and 35 m; taking 3 vertical top-view cross sections in the vertical direction of the explosive pile, wherein the distances from the positions to the bottom of the step are respectively 4m, 8m and 10m, and the transverse side-view cross section 1 refers to fig. 2;
3) importing the section diagram into MATLAB to perform image processing, wherein the processing program comprises image conversion (converting into a gray-scale image), image enhancement, image segmentation and image edge detection, and the processing result of the section 1 refers to FIG. 3;
4) determining a proportional size based on the three-dimensional model of the blasting pile, quantizing each block degree of the processed image to obtain the size of each block degree, counting the block degree distribution of the whole section diagram to obtain the block degree distribution data of the section 1, and referring to fig. 4 and 5;
5) repeating the steps 3) and 4) to obtain the block size distribution of all the section diagrams, and obtaining the block size distribution data of the explosive pile by an averaging value method, referring to the figure 6;
6) manually analyzing the position relation between the part with large block degree ratio in each section and the section in the pile-blasting process, and determining the spatial position distribution rule of each block degree interval;
7. manually comparing the explosive pile block size distribution with a target, analyzing the blasting quality, and determining a blasting optimization scheme;
8. and (4) repeating the steps 1-7 according to the optimized blasting parameters until the target block size is reached.
The method comprises the steps of carrying out numerical simulation by using GDEM software, cutting the model in a grid-dividing mode, carrying out image processing on a profile by using MATLAB, quantizing each block of the processed image, analyzing each section block, and obtaining the overall explosive block distribution by averaging each section block. The method has the advantages of high precision, thoroughness, no damage, safety, intuition and the like, and has important significance for predicting the blasting block size.
The above-described embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solutions of the present invention can be made by those skilled in the art without departing from the spirit of the present invention, and the technical solutions of the present invention are within the scope of the present invention defined by the claims.
Claims (3)
1. A method for predicting blasting blockiness of a gridding segmentation numerical simulation blasting pile model is characterized by comprising the following steps:
1) performing numerical simulation by using GDEM software to obtain a three-dimensional explosive stack model;
2) dividing a grid of the three-dimensional explosive stack model to intercept any number of explosive stack transverse and longitudinal section diagrams;
3) performing image processing on the sectional diagram by using MATLAB, wherein the processing procedures comprise gray level image conversion, image enhancement, image segmentation and image edge detection;
4) determining the proportional size based on the three-dimensional model of the blasting pile, quantizing each block degree of the processed image to obtain the size of each block degree, and counting the block degree distribution of the whole profile to obtain the block degree distribution data of the profile;
5) repeating the steps 3) and 4) to obtain the block degree distribution of all the section diagrams, and obtaining the block degree distribution data of the blasting pile by an averaging method;
6) manually analyzing the spatial position distribution rule of each block degree interval;
7) comparing the block size distribution obtained in the steps 5) and 6) with a target, and determining a blasting optimization scheme;
8) and adjusting blasting parameters in GDEM software, and repeating the steps 1) to 7) until the target block size is reached.
2. The method for predicting the blasting blockiness of the gridding segmentation numerical simulation blasting pile model according to claim 1, wherein the step 1) is specifically as follows: and setting model calculation parameters by using GDEM software according to blasting design parameters and geological data, and carrying out numerical simulation on the blasting area to obtain a three-dimensional blasting pile model.
3. The method for predicting the blasting bulkiness of the grid-cut numerical simulation blasting pile model according to claim 1, wherein in the step 2), 3 transverse side-looking sections are taken from the blasting pile extension direction of the three-dimensional blasting pile model by grid division, and the distances from the positions to the left side edge of the step are respectively 15m, 20m and 35 m; 3 vertical overlooking sections are taken in the vertical direction of the blasting pile, and the distances from the positions to the bottom of the step are respectively 4m, 8m and 10 m.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010014781.5A CN111199108A (en) | 2020-01-07 | 2020-01-07 | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010014781.5A CN111199108A (en) | 2020-01-07 | 2020-01-07 | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111199108A true CN111199108A (en) | 2020-05-26 |
Family
ID=70746538
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010014781.5A Pending CN111199108A (en) | 2020-01-07 | 2020-01-07 | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111199108A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111814372A (en) * | 2020-07-02 | 2020-10-23 | 大昌建设集团有限公司 | Blasting blockiness control method combining numerical calculation and blockiness screening |
| CN112991539A (en) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements |
| CN112991538A (en) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for explosive stack body based on Weibull distribution of discrete element pair stack shape |
| CN113284241A (en) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for surface rock mass of blasting pile with pile shape distributed in Weibull mode |
| CN113674406A (en) * | 2021-03-11 | 2021-11-19 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937397A (en) * | 2012-11-27 | 2013-02-20 | 武汉大学 | Step blasting dynamic design method based on close-range photogrammetry technology |
| CN103699751A (en) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | Sand body reservoir architecture modeling method and system based on space vectors |
| CN105224742A (en) * | 2015-09-29 | 2016-01-06 | 鞍钢集团矿业公司 | A kind of blast action of analyzing is got out of a predicament or an embarrassing situation the method for slope stability |
| CN105891886A (en) * | 2014-12-08 | 2016-08-24 | 北京恒欣泰克科技有限公司 | Geologic body automatic analysis method based on three-dimensional stratum model |
| CN107391837A (en) * | 2017-07-19 | 2017-11-24 | 中南大学 | A kind of aperture charge constitution parameter optimization method of fanhole(s) |
| CN108073774A (en) * | 2017-12-29 | 2018-05-25 | 辽宁科技大学 | A kind of reliable method for verifying quick-fried heap LUMPINESS DISTRIBUTION |
| CN109887073A (en) * | 2018-08-16 | 2019-06-14 | 清能艾科（深圳）能源技术有限公司 | The three-dimensional digital model construction method and device of rock core |
| CN110414341A (en) * | 2019-06-24 | 2019-11-05 | 中国矿业大学 | A kind of Blasting In The Open Mining heap measurement statistical method based on unmanned air vehicle technique |
-
2020
- 2020-01-07 CN CN202010014781.5A patent/CN111199108A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937397A (en) * | 2012-11-27 | 2013-02-20 | 武汉大学 | Step blasting dynamic design method based on close-range photogrammetry technology |
| CN103699751A (en) * | 2013-12-30 | 2014-04-02 | 中国石油大学(北京) | Sand body reservoir architecture modeling method and system based on space vectors |
| CN105891886A (en) * | 2014-12-08 | 2016-08-24 | 北京恒欣泰克科技有限公司 | Geologic body automatic analysis method based on three-dimensional stratum model |
| CN105224742A (en) * | 2015-09-29 | 2016-01-06 | 鞍钢集团矿业公司 | A kind of blast action of analyzing is got out of a predicament or an embarrassing situation the method for slope stability |
| CN107391837A (en) * | 2017-07-19 | 2017-11-24 | 中南大学 | A kind of aperture charge constitution parameter optimization method of fanhole(s) |
| CN108073774A (en) * | 2017-12-29 | 2018-05-25 | 辽宁科技大学 | A kind of reliable method for verifying quick-fried heap LUMPINESS DISTRIBUTION |
| CN109887073A (en) * | 2018-08-16 | 2019-06-14 | 清能艾科（深圳）能源技术有限公司 | The three-dimensional digital model construction method and device of rock core |
| CN110414341A (en) * | 2019-06-24 | 2019-11-05 | 中国矿业大学 | A kind of Blasting In The Open Mining heap measurement statistical method based on unmanned air vehicle technique |
Non-Patent Citations (8)
| Title |
|---|
| 于永江等: "煤体爆堆块度分布的测试", 《煤炭学报》, no. 03 * |
| 吴发名等: "基于原生节理统计和爆破裂纹模拟的堆石料块度分布预测", 《岩石力学与工程学报》, no. 06 * |
| 周志刚等: "某采场大规模爆破施工振动监测分析", 《现代矿业》, no. 09 * |
| 张旭等: "一种新的声速剖面结构参数化方法", 《海洋学报(中文版)》, no. 05 * |
| 荆永滨等: "基于块体形状分类的爆堆块度三维几何参数估算", 《爆破》, no. 01 * |
| 郑炳旭等: "炸药单耗对赤铁矿爆破块度的影响规律数值模拟研究", 《爆破》, no. 03, pages 63 - 66 * |
| 陈柯竹等: "基于CDEM的隧道卸压爆破及岩爆抑制效应模拟", 《工程地质学报》, no. 03 * |
| 高克林等: "台阶炮孔排间毫秒延时爆破爆堆形状的计算机模拟", 《爆破》, no. 04 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111814372A (en) * | 2020-07-02 | 2020-10-23 | 大昌建设集团有限公司 | Blasting blockiness control method combining numerical calculation and blockiness screening |
| CN111814372B (en) * | 2020-07-02 | 2023-12-29 | 大昌建设集团有限公司 | Blasting block control method combining numerical calculation and block screening |
| CN112991539A (en) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for explosive pile and block size distribution based on discrete elements |
| CN112991538A (en) * | 2021-03-11 | 2021-06-18 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for explosive stack body based on Weibull distribution of discrete element pair stack shape |
| CN113284241A (en) * | 2021-03-11 | 2021-08-20 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for surface rock mass of blasting pile with pile shape distributed in Weibull mode |
| CN113674406A (en) * | 2021-03-11 | 2021-11-19 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for surface rock mass of open-air deep hole bench blasting pile |
| CN113284241B (en) * | 2021-03-11 | 2022-03-15 | 广东锡源爆破科技股份有限公司 | Three-dimensional simulation method for surface rock mass of blasting pile with pile shape distributed in Weibull mode |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111199108A (en) | Blasting blockiness prediction method of gridding segmentation numerical simulation blasting pile model | |
| Wang et al. | Analyses and predictions of rock cuttabilities under different confining stresses and rock properties based on rock indentation tests by conical pick | |
| CN110617074A (en) | Incidence relation method for ground settlement and tunneling parameters in shield construction | |
| CN112364560B (en) | Intelligent prediction method for working hours of mine rock drilling equipment | |
| CN107368463B (en) | Roadway nonlinear deformation prediction method based on fiber bragg grating sensor network data | |
| CN103591942B (en) | Method for quickly simplifying and automatically optimizing sound velocity profile based on maximum offset of sound velocity | |
| CN108536979B (en) | Underflow concentration prediction method based on thickener mechanism model | |
| CN113569487B (en) | BP neural network-based method for predicting step blasting throwing effect | |
| CN103761421B (en) | A kind of large-sized mining dredger hoisting mechanism reliability estimation method | |
| CN114386679A (en) | Ball mill energy consumption state monitoring method based on ensemble learning | |
| CN111160490A (en) | Deep learning dangerous rock deformation prediction method and device based on multiple time sequences | |
| CN116049654B (en) | Safety monitoring and early warning method and system for coal preparation equipment | |
| CN110826280A (en) | Process optimization method for improving cylindrical part drawing lug based on finite element simulation | |
| CN116660263A (en) | Grouting effect evaluation method based on digital image processing | |
| CN106599511B (en) | Method for optimizing coal mining method of thin coal seam longwall fully mechanized coal mining face | |
| CN111581771A (en) | Stamping workpiece cracking prediction platform based on artificial intelligence technology | |
| CN113849921B (en) | Method for evaluating pretreatment effect of large data sample of wear of cutter of tunneling machine | |
| CN111144032A (en) | Profile function method for predicting near-horizontal coal seam longwall mining surface subsidence | |
| CN110781615A (en) | CAE simulation accuracy evaluation method | |
| CN116522635A (en) | Stone blasting method, system, terminal equipment and storage medium | |
| CN117647447A (en) | Research method for mechanical properties of cross-fracture sandstone and crack propagation | |
| CN113586157B (en) | Extraction working face salient dangerous area rapid division method based on Kriging interpolation | |
| CN117171920B (en) | Aluminum foil cutlery box stamping forming control method and system | |
| Xing et al. | Equipment Pattern Recognition of Unbalanced Fuel Consumption Data Based on Grouping Multi-BP Neural Network | |
| CN117732886B (en) | Hot rolling quality pre-control method based on cascading intelligent diagnosis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200526 |
|
| RJ01 | Rejection of invention patent application after publication |