CN107391837B - 一种扇形炮孔的孔口装药结构参数优化计算方法 - Google Patents
一种扇形炮孔的孔口装药结构参数优化计算方法 Download PDFInfo
- Publication number
- CN107391837B CN107391837B CN201710590516.XA CN201710590516A CN107391837B CN 107391837 B CN107391837 B CN 107391837B CN 201710590516 A CN201710590516 A CN 201710590516A CN 107391837 B CN107391837 B CN 107391837B
- Authority
- CN
- China
- Prior art keywords
- lumpiness
- value
- scheme
- matrix
- fanhole
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000005457 optimization Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000013467 fragmentation Methods 0.000 claims abstract description 8
- 238000006062 fragmentation reaction Methods 0.000 claims abstract description 8
- 238000005315 distribution function Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims abstract description 3
- 239000000843 powder Substances 0.000 claims description 21
- 239000011435 rock Substances 0.000 claims description 15
- 238000004880 explosion Methods 0.000 claims description 13
- 238000007873 sieving Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 125000006850 spacer group Chemical group 0.000 abstract description 2
- 239000002360 explosive Substances 0.000 description 16
- 238000005422 blasting Methods 0.000 description 8
- 238000009826 distribution Methods 0.000 description 6
- 208000037656 Respiratory Sounds Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 241000237509 Patinopecten sp. Species 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 235000020637 scallop Nutrition 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种扇形炮孔的孔口装药结构参数优化计算方法,包括以下步骤:1)采用穷举算法获得所有常规孔布设方案;从中筛选出可行的方案;计算每种可行方案中非常规孔中的药柱长度;2)对爆区三维模型进行网格划分;3)以原岩块度分布函数作为概率密度函数,连续进行多次随机抽样得到原岩块度矩阵λo;4)求解各个方案下,爆破作用下的块度矩阵λ1;5)将λo中的元素值与λ1中相应的元素一一对比,取其中取较小的值,得到最终的爆破块度矩阵;6)根据块度统计参数对各个方案进行优选,得到相应的扇形炮孔的孔口间隔结构参数,即各个炮孔的装药长度。本发明可自动计算孔口装药结构参数,简单快速,结果精确,有较强的可操作性。
Description
技术领域
本发明涉及矿山及其他工程爆破中使用扇形炮孔爆破施工时,扇形炮孔的孔口装药结构参数优化计算方法。
背景技术
受扇形炮孔几何形态的影响,其炸药在扇面内的面密度势必是不均匀的,如图1所示,随着炮孔深度的增加,孔间距也在增加,使得现场爆破中扇面内的炸药面密度存在不均匀的情况,其炸药面密度依次降低,爆破能量也随之依次降低。一般通过严格控制各炮孔的深度,并使其具有相同的孔底距,来控制采场边界以及远区的岩石破碎块度,但是,近区炮孔的炸药分布就不可能通过炮孔的布设来控制,即使通过增加炮孔填塞长度降低孔口位置的炸药能量密度,或者使用孔底起爆避免爆破过程中孔口位置的应力波叠加,然而,近区位置的炸药面密度偏大的问题依然是存在的。通过改变扇形炮孔装药结构【即炸药药卷在炮孔(眼)内的安置方式】参数可以一定程度上改善孔口的炸药能量密集的问题,目前,大多矿山均采用间隔增加炮孔填塞长度的形式,即选取一些炮孔,装药至孔口,填塞长度等于最小不装药长度(下文简称为常规孔),中间间隔的若干炮孔,增加其填塞长度(下文简称为非常规孔),以减少孔口区域的炸药面密度。但是,目前现场中在计算扇形炮孔装药结构参数,如设计各个炮孔的装药长度,尤其是常规孔的炮孔的数量与位置时,由于缺乏一个合理的确定方法以及精确的计算依据,主要依靠人工经验,主观因素较大,而且很难做到对每一排扇形炮孔都根据炮孔的几何参数与自由面的形态实时调整,以取得更好的爆破效果,因此,需要避免主观经验的干扰,设计一种合理有效的计算方法,能够根据爆破条件的变化,实时、精确的确定最合理的扇形炮孔孔口装药结构参数,获得更好的爆破效果。
发明内容
本发明所解决的技术问题是,针对现有技术的不足,提供一种全新的,精确的,快速的扇形炮孔的孔口装药结构参数优化计算方法,根据每排炮孔的炮孔、自由面状况,分别进行优化计算,尤其适用于多排炮孔的大规模爆破。
本发明所提供的技术方案为:
一种扇形炮孔的孔口装药结构参数优化计算方法,包括以下步骤:
步骤1、采用穷举算法获得所有常规孔布设方案,作为初步方案。首先基于输入的炮孔几何参数(包括各个炮孔的孔口坐标、孔深、孔径和倾角,可通过AutoCAD的图形对象进行直接选取输入)、炸药参数(炸药密度、爆速、绝热指数等常数),根据穷举算法获得所有的常规孔布设方案,每种方案包含常规孔的数量与位置情况;
图2说明了本程序中确定常规孔布置方案的穷举算法孔,如图所示为一组包含10个炮孔的扇形炮孔,d0为孔底距,d1为相邻药柱最小安全距离(药柱影响范围),其中No.1炮孔为已知常规孔,现以此为基准,顺时针方向确定相邻常规孔可行的选取范围,设定常规孔布置方案时的两个基本原则为:一是常规孔的不相邻原则,即任意两个常规孔之间至少应当间隔一个非常规孔,二是常规孔与非常规孔的装药规则,即常规孔的填塞长度应当等于最小允许不装药长度B之外,非常规孔的填塞长度都应当大于最小允许不装药长度B。
首先,由于No.1炮孔为已知常规孔,根据常规孔的不相邻原则,间隔一个炮孔,确定No.3炮孔为距离最近的可选常规孔。
然后,寻找最远的可选边界,如图,由近及远,逐一检验当各个炮孔分别作为常规孔最远的可选边界时,是否满足上述的两个基本原则;当No.7孔作为常规孔时,计算No.1与No.7之间的No.4孔的填塞长度小于等于其最小不装药长度B【根据No.1孔和No.7孔的装药长度及其相邻药柱最小安全距离d1,按照目前普遍采用的邻近药包最小距离原则(即药柱最底端到相邻药柱间的距离均不能小于相邻药柱最小安全距离)[1-3],以d1为半径,作出No.1孔与No.7孔的药柱影响范围,即图2中虚线区域,计算两虚线区域与No.1与No.7之间的各个炮孔的交点,该交点即为各个炮孔的装药段和填塞段的分界点,从而得到No.1与No.7之间的各个炮孔的填塞长度,若得到的某个炮孔的填塞长度小于最小不装药长度B,或者两虚线区域与No.1与No.7之间的某个炮孔无交点,则填塞长度取为B;其中炮孔中的装药长度=炮孔长度-填塞长度,常规孔的填塞长度等于最小不装药长度B】,即此时该炮孔为不符合非常规孔的装药规则,这是不合理的,所以,No.6孔为最远的可选择的小组边界。即No.3~6为以No.1孔为基准的下一个相邻常规孔的选择范围。
同样的,再分别以No.3~6孔中的常规孔作为基准,按照同一方向寻找下一个常规孔的范围,直至扇形炮孔的边界。
步骤2、给定本次爆破的具体工艺约束,包括对边孔的装药长度要求,装药结构有无对称性等要求等;从初步方案中,筛选出满足本次爆破的具体工艺约束的常规孔布设方案;;
步骤3、根据步骤2筛选出的每种方案中常规孔的位置与数量,按照目前普遍采用的邻近药包最小距离原则,计算每种方案中所有非常规孔中的装药长度[1-3];步骤1中计算非常规孔的填塞长度目的在于检验,只要找到不满足设定的两个基本原则的非常规孔,即可停止计算,一般只需计算位于两个假定的常规孔中间位置(危险位置)的一个或两个非常规孔的填塞长度,无需计算所有非常规孔的装药长度;而本步骤则是计算筛选出的每种方案中所有非常规孔中的装药长度;
步骤4、建立爆区三维模型,并对其进行网格划分,设x,y,z三个方向上的网格节点数目分别为nx,ny,nz;导入爆区边界、巷道轮廓线(可通过AutoCAD的图形对象直接选取输入,亦可通过边界函数的形式进行定义)与计算模型厚度信息,建立爆区三维模型,按照给定的x,y,z三个方向上的网格尺寸进行网格划分,x,y,z三个方向上网格尺寸都是各自均匀的;
步骤5、求解爆破前的原岩块度矩阵λo(nx,ny,nz);
式中,fsitu(D)为岩石过筛的比例(体积);D为块度值;D50situ为岩石过筛的比例为50%时的网筛尺寸。
以原岩块度分布函数作为概率密度函数,连续进行nx×ny×nz次随机抽样,得到nx×ny×nz个块度值,分别赋值给爆区三维模型中的各个网格节点,得到原岩块度矩阵λo(nx,ny,nz),其中nx,ny,nz分别表示矩阵λo种x,y,z三个方向上的元素个数,矩阵λo中的元素表示爆破前爆区三维模型中各个网格节点的块度值;
步骤6、求解各个方案下,爆破作用下的块度矩阵λ1(nx,ny,nz);
采用等效叠加原理,使用一系列球状药包单元(高径比小于6的药包)的叠加效果等效柱状药包的爆破作用。
根据各个方案下各炮孔中的药柱长度,将所有药柱分割成一系列直径等于炮孔直径的球状药包单元;基于Harries数学模型,依次计算每种方案下爆区三维模型所有网格节点的块度值,得到爆破作用下的块度矩阵λ1(nx,ny,nz);
步骤7、将λ1(nx,ny,nz)中的元素值λ1(i,j,k)与λo(nx,ny,nz)中相应的元素值λ0(i,j,k)进行对比,取其中取较小的值,得到最终的块度值λfnl(i,j,k),即:
λfnl(i,j,k)=Min(λ0(i,j,k),λ1(i,j,k))
式中,Min表示取较小的值;i=[1,nx],j=[1,ny],k=[1,nz],i,j,k∈N;λfnl(i,j,k)表示最终计算得到的爆区三维模型中各个网格节点的块度值。
如果λ1(i,j,k)小于初始值λ0(i,j,k),则认为该节点主要受到爆破作用,令λfnl(i,j,k)其等于λ1(i,j,k),反之,若λ1(i,j,k)大于初始值λ0(i,j,k),则认为爆破作用不足以影响其原始块度,保持其初始值λ0(i,j,k)不变。经过原岩块度矩阵与爆破作用下的块度矩阵的叠加,得到λfnl(i,j,k),由λfnl(i,j,k)构成最终的爆破块度矩阵λfnl(nx,ny,nz)。
步骤8、得到各方案得到的最终的爆破块度矩阵之后,根据块度统计参数对方案进行优选,得到相应的扇形炮孔的孔口间隔结构参数,即各个炮孔的装药长度。
进一步地,所述步骤5中,输入岩石参数与原岩块度分布参数,按照负指数模型生成原岩块度分布函数:
式中,fsitu(D)为岩石过筛的比例(体积);D为块度值;D50situ为岩石过筛的比例为50%时的网筛尺寸。
进一步地,所述步骤6中,根据各方案下的球状药包单元集合,采用本发明改进的Harries块度模型,计算各网格节点的块度值;
根据原始的Harries数学模型,如图3所示,平面内距离球状药包中心为r的区域可能出现的裂纹条数为:
式中:r为距球状药包中心的距离;T为岩石极限抗拉应变值;b为炮孔半径;ε为孔壁应变值,可由下式计算得出:
其中:ν为岩石泊松比;p为炸药产生的气体压力,单位为MPa;ρ为岩石密度,单位为g/cm3;CP为岩石纵波速度,单位为m/s;K为绝热指数。
则根据裂纹条数,估算距球状药包中心的距离为r的区域的块度值为:
对于柱状药包,根据等效叠加原理,可以将柱状药包的作用等效为一系列球状药包单元共同作用的结果(A.M.Starfield)。如图4所示,柱状药包被分割成m个半径为b的球状药包单元,在点P(μ,θ,h)处叠加得到的可能出现的裂纹条数可由公式(4)得出。
式中,μ为点P到柱状药包中轴线的距离,h为点P到柱状药包孔口的距离;ri为点P至第i个球状药包单元中心的距离(见图4)。
另外,容易得到图中存在如下的几何关系:
ri 2=μ2+(2bi-b-h)2 (6)
结合(4)、(5)和(6),得到距离柱状药包中轴线为μ的网格节点在爆破作用下的块度值为:
进一步地,所述步骤8中,块度统计参数为D10,D50和D90,分别表示能通过10%、50%、90%的爆堆体积的网筛尺寸;其中D50反映了爆堆的平均尺寸,D90-D10表示大部分块度的尺寸波动情况,反映了块度分布的集中程度。在实际生产中,D50涉及各个环节的经济成本,过大则会增加运输、磨矿成本,过小则增加爆破成本,粉矿也可能较多,增加运输过程损失,其取值是一个最优化问题,这里我们认为D50的波动仅仅考虑由于不同的装药结构造成的,其取值落在给定的合理波动区间(D50min,D50max)内即可。而另外一个参数D90-D10,表示爆破块度的块度分布宽度,一定程度上反映了爆破块度的集中程度,其值越小,代表块度更加均匀,粉矿与大块的比例更少。爆区内x,y,z三个方向上网格尺寸都是各自均匀的,所以,对所有网格节点的块度值从小到大进行排序,处在10%,50%以及90%位置的块度值即为D10,D50和D90。
首先,分别将每种方案相应的爆破块度矩阵λfnl(nx,ny,nz)中所有网格节点的块度值从小到大进行排序,处在10%,50%以及90%位置的块度值即为D10,D50和D90;
然后,筛选出所有D50落在合理的区间(D50min,D50max)内的方案,区间(D50min,D50max)的取值根据人为经验给定;
最后,取块度波动范围最小,即D90-D10的值最小的方案,判定为最优方案。
有益效果:
与传统的扇形炮孔参数确定方法相比,本方法有如下几个优点:
1.面向每排炮孔的具体爆破条件,对装药长度进行合理有效的最优化计算;
2.根据输入的炸药、岩石、自由面以及炮孔参数即可完成自动计算,简单方便,有较强的可操作性;
3.在不同的矿山、工程爆破中均可使用,具有广泛地适用性。
附图说明
图1为扇形炮孔几何形态;
图2为一组包含10个炮孔的扇形炮孔;
图3为Harries数学模型原理图
图4为本发明柱状药包爆破模型原理图;
图5为本发明流程图;
图6为本发明实施例中爆区三维模型;
附图标记说明:
1—装药段;2—填塞段;3—巷道;3—巷道;4—爆破区域;5—钻孔边界;6—计算边界。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
大冶有色铜绿山铜铁矿位于湖北省大冶市,年产矿石量约80万吨,采用分段空场嗣后充填法作为主要的采矿方法之一。采场垂直矿体布置,长度为矿体厚度,中段高度60.0m,分4个分层,每个分层高13.0m,采场宽度为8.0m,炮孔直径64mm,扇形炮孔排距1.6m,孔底距1.9m,孔口的最小不装药长度取1.2m。
选用某测试采场顶部分段侧崩区第一排扇形孔为例进行装药结构优化计算。模型的相关参数如图6所示,模型块度为9.0m,高度为13.5m,计算模型边界距炮孔边界0.5m,模型厚度等于排距1.6m。相关的炮孔几何参数、炸药参数以及岩石参数见表1和表2。
表1计算模型炮孔几何参数表
表2计算案例程序输入参数表
输入相关岩石、炸药以及炮孔几何参数,按本发明步骤解算出9种可行的(满足本次爆破的具体工艺约束)装药方案,表3给出了计算得到的每种装药方案的炮孔装药长度、常规孔的分布以及每种装药方案得到的块度矩阵的几个关键的块度分布参数。
表3各方案装药参数及块度参数统计表
标有*的炮孔在相应的方案中作为常规孔
根据关键的块度统计参数D10、D50以及D90,计算得到所有D50落在区间(0.15,,0.30)内的方案中,方案E的D90-D10值最小,则方案E为最优方案。
参考文献:
[1]汪旭光.爆破手册[M].冶金工业出版社.2009
[2]张铁英,王杰.扇形中、深孔间隔装药理论及其应用[J].黄金,1989,(09):15-21.
[3]张铁英,王杰.扇形深孔合理装药长度的计算[J].黄金,1988,(01):12-15.
Claims (4)
1.一种扇形炮孔的孔口装药结构参数优化计算方法,其特征在于,包括以下步骤:
步骤1、采用穷举算法获得所有常规孔布设方案,作为初步方案;每种方案包含常规孔的数量与位置情况;
步骤2、从初步方案中,筛选出满足本次爆破的具体工艺约束的常规孔布设方案;
步骤3、根据步骤2筛选出的每种方案中常规孔的位置与数量,按照邻近药包最小距离原则,即药柱最底端到相邻药柱间的距离均不能小于相邻药柱最小安全距离的原则,计算每种方案中非常规孔中的药柱长度;
步骤4、建立爆区三维模型,并对其进行网格划分,设x,y,z三个方向上的网格节点数目分别为nx,ny,nz;
步骤5、求解爆破前的原岩块度矩阵λo(nx,ny,nz);
以原岩块度分布函数作为概率密度函数,连续进行nx×ny×nz次随机抽样,得到nx×ny×nz个块度值,分别赋值给爆区三维模型中的各个网格节点,得到原岩块度矩阵λo(nx,ny,nz),矩阵λo中的元素表示爆破前爆区三维模型中各个网格节点的块度值;
步骤6、求解步骤2筛选出的各个方案下,爆破作用下的块度矩阵λ1(nx,ny,nz);
根据步骤2筛选出的各个方案下各炮孔中的药柱长度,将所有药柱分割成一系列直径等于炮孔直径的球状药包单元;基于Harries数学模型,依次计算步骤2筛选出的每种方案下爆区三维模型所有网格节点的块度值,得到爆破作用下的块度矩阵λ1(nx,ny,nz);
步骤7、将λ1(nx,ny,nz)中的元素值λ1(i,j,k)与λo(nx,ny,nz)中相应的元素值λ0(i,j,k)进行对比,取其中取较小的值,得到最终的块度值λfnl(i,j,k),即:
λfnl(i,j,k)=Min(λ0(i,j,k),λ1(i,j,k))
式中,Min表示取较小的值;i=[1,nx],j=[1,ny],k=[1,nz],i,j,k∈N;λfnl(i,j,k)表示最终计算得到的爆区三维模型中网格节点(i,j,k)的块度值;
由λfnl(i,j,k)构成最终的爆破块度矩阵λfnl(nx,ny,nz);
步骤8、得到步骤2筛选出的各方案对应的最终的爆破块度矩阵λfnl(nx,ny,nz)之后,根据块度统计参数对步骤2筛选出的各个方案进行优选,得到相应的扇形炮孔的孔口间隔结构参数,即各个炮孔的装药长度。
2.根据权利要求1所述的扇形炮孔的孔口装药结构参数优化计算方法,其特征在于,所述步骤5中,原岩块度分布函数为:
式中,fsitu(D)为岩石过筛的比例;D为块度值;D50situ为岩石过筛的比例为50%时的网筛尺寸。
3.根据权利要求2所述的扇形炮孔的孔口装药结构参数优化计算方法,其特征在于,所述步骤6中,距离柱状药包中轴线为μ的网格节点P在爆破作用下的块度值为:
式中,μ为点P到柱状药包中轴线的距离,h为点P到柱状药包孔口的距离;b为柱状药包分割成的球状药包单元半径,m为柱状药包分割成的球状药包单元的个数;T为岩石极限抗拉应变值;ε为孔壁应变值。
4.根据权利要求1~3中任一项所述的扇形炮孔的孔口装药结构参数优化计算方法,其特征在于,所述步骤8中,块度统计参数为D10,D50和D90,分别表示能通过10%、50%、90%的爆堆体积的网筛尺寸;
首先,分别将每种方案相应的爆破块度矩阵λfnl(nx,ny,nz)中所有网格节点的块度值从小到大进行排序,处在10%,50%以及90%位置的块度值即为D10,D50和D90;
然后,筛选出所有D50落在合理的区间(D50min,D50max)内的方案,区间(D50min,D50max)的取值根据人为经验给定;
最后,取块度波动范围最小,即D90-D10的值最小的方案,判定为最优方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710590516.XA CN107391837B (zh) | 2017-07-19 | 2017-07-19 | 一种扇形炮孔的孔口装药结构参数优化计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710590516.XA CN107391837B (zh) | 2017-07-19 | 2017-07-19 | 一种扇形炮孔的孔口装药结构参数优化计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107391837A CN107391837A (zh) | 2017-11-24 |
CN107391837B true CN107391837B (zh) | 2019-08-16 |
Family
ID=60335697
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710590516.XA Expired - Fee Related CN107391837B (zh) | 2017-07-19 | 2017-07-19 | 一种扇形炮孔的孔口装药结构参数优化计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107391837B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256253B (zh) * | 2018-01-29 | 2022-05-24 | 核工业湖州工程勘察院 | 一种基于Voronoi网格剖分的等面积炮孔布置方法 |
CN108388922A (zh) * | 2018-03-08 | 2018-08-10 | 云南电网有限责任公司电力科学研究院 | 一种用于火电机组主辅设备状态异常的识别方法及系统 |
CN111199108A (zh) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院(有限责任公司) | 一种网格化切分数值模拟爆堆模型的爆破块度预测方法 |
CN111814372B (zh) * | 2020-07-02 | 2023-12-29 | 大昌建设集团有限公司 | 一种联合数值计算和块度筛分的爆破块度控制方法 |
CN114792047B (zh) * | 2022-03-30 | 2024-08-20 | 西安交通大学 | 一种变密度高真空多层绝热结构的优化方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0021843D0 (en) * | 2000-09-06 | 2000-10-18 | Courtney William A | Fluid filled impact absorber |
CN202119900U (zh) * | 2011-06-21 | 2012-01-18 | 北京蓝尊科技有限公司 | 基于爆破震源定位技术的矿山地下开采活动实时监测系统 |
WO2012178182A1 (en) * | 2011-06-23 | 2012-12-27 | Rok Protective Systems, Inc. | Blast resistant pipe protection system and method |
CN102620612B (zh) * | 2012-04-09 | 2016-03-09 | 卓利维(北京)科技有限公司 | 一种井下矿中深孔爆破方法 |
CN103093390B (zh) * | 2013-01-25 | 2016-08-31 | 河南理工大学 | 一种基于联合聚类的煤矿安全评价系统 |
CN106767209B (zh) * | 2017-01-17 | 2018-08-07 | 长沙有色冶金设计研究院有限公司 | 一种炮孔装药及堵塞装置 |
-
2017
- 2017-07-19 CN CN201710590516.XA patent/CN107391837B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107391837A (zh) | 2017-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107391837B (zh) | 一种扇形炮孔的孔口装药结构参数优化计算方法 | |
Monjezi et al. | Optimization of open pit blast parameters using genetic algorithm | |
Wang et al. | Charge design scheme optimization for ring blasting based on the developed Scaled Heelan model | |
CN106650052A (zh) | 一种基于人工神经网络的配料爆破参数智能设计方法 | |
CN113255175A (zh) | 一种层状软岩爆破数值模拟方法 | |
CN106895755B (zh) | 一种露天台阶爆破智能化设计方法 | |
CN108073786A (zh) | 页岩可压性评价方法及装置 | |
CN104929641B (zh) | 用于大倾角煤层顶板灾害的治理方法 | |
CN104390534B (zh) | 一种隧道光面爆破质量控制方法 | |
CN110152771B (zh) | 物料破碎腔结构及多级嵌套式物料破碎腔结构设计方法 | |
CN107169616B (zh) | 矿井未采区构造相对复杂程度的相对熵预测方法 | |
CN110132080A (zh) | 一种露天矿高台阶孔内微差爆破方法 | |
Yari et al. | Developing a mathematical assessment model for blasting patterns management: Sungun copper mine | |
CN102997768B (zh) | 一种爆破岩石的方法 | |
CN108256253A (zh) | 一种基于Voronoi网格剖分的等面积炮孔布置方法 | |
Ding et al. | Distribution characteristics of fragments size and optimization of blasting parameters under blasting impact load in open-pit mine | |
CN106969676A (zh) | 一种基于rtk炮孔定位技术改善爆破效果的布孔微调方法 | |
CN110553561B (zh) | 一种中深孔分段微差塌落的爆破方法 | |
CN111894590A (zh) | 采矿爆破方法 | |
CN116086265A (zh) | 露天台阶的低压损低粉尘爆破方法 | |
CN214666419U (zh) | 一种露天矿山提高延米爆破量的深孔微差挤压爆破结构 | |
CN106599511B (zh) | 一种薄煤层长壁综采工作面采煤方法优选的方法 | |
CN114139438A (zh) | 一种高炉炉料轨迹模型构建方法 | |
CN112668222B (zh) | 一种井下煤仓装煤过程中仓壁受力特征的模拟方法 | |
CN118153456B (zh) | 一种数码延期模块及控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190816 |