CN113674406A - 一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法 - Google Patents

Abstract

本申请公开了一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，涉及爆破效果评价技术领域，包括爆堆岩块形态与爆堆三维形态的确定、爆堆表面局部区域划分与块度统计、爆堆表面岩块三维模拟、爆堆表面岩块三维模拟结果。本发明通过把爆堆表面进行分区，并对每个区域的岩块进行统计分析，获得该区域爆堆块度级配累积曲线和大块率，构成爆堆表面块度局部分布模型；把该模型和爆堆形态信息结合起来，通过改进离散单元构造算法中的向内排列方案形成爆堆表面岩块生成算法，并采用三维可视化技术实现对爆堆形态及表面岩块三维模型模拟功能。

Description

一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法

技术领域

本申请涉及爆破效果评价技术领域，尤其涉及一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法。

背景技术

露天台阶深孔爆破块度分布规律及其直观表达的研究是评价爆破效果的重要手段，也是反馈和改进露天台阶深孔爆破设计和施工方法的重要依据。由于露天深孔台阶爆破块度的检测方法主要为表面块度检测，内部块度是通过间接方法获得，露天台阶深孔爆破块度分布规律研究多集中于表面块度分布规律研究，主要研究内容和方法有：(1)露天深孔台阶爆破块度的检测方法和快速识别方法研究；(2)露天深孔台阶爆破块度分布规律与爆破参数之间的数学模型研究。如：中南大学聂军等人选用基因表达式编程(GEP)算法，建立基于GEP的爆破块度预测模型；中南大学史秀志等人建立了基于LS-SVR的岩石爆破块度预测模型，预测了小样本条件下的矿山岩体爆破块度；天津大学王仁超等人基于LM算法建立了爆破块度神经网络预测模型。

这些研究方法和研究内容存在的主要问题有：(1)露天深孔台阶爆破块度分布规律的数学模型过于抽象，不利于技术人员的掌握和成果的推广应用；(2)露天深孔台阶爆破块度分布规律的数学模型是统计学上的数学模型，反映的是爆堆表面块度的平均级配模型，不能反映爆堆表面块度的实际分布状况。而露天深孔台阶爆破由于同次爆破的工程岩体可能由多种差异性较大的地质体组成，再考虑炸药在爆堆内分布的不均匀性，以及露天台阶深孔爆破作用机理，爆堆表面块度分布存在区域性差异。

为此我们需要提出一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法。

发明内容

本申请实施例提供一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，包括：

步骤S1、爆堆岩块形态与爆堆三维形态的确定；

步骤S2、爆堆表面局部区域划分与块度统计：

对爆堆表面进行区域划分，并对每个区域建立对应的块度分布曲线；按以下方法和步骤进行区域划分和区域块度统计，形成局部块度分布曲线；

步骤1：首先以梯形爆堆的6个表面进行区域划分，然后根据大块在爆堆顶面和坡面的位置及数量进一步划分局部区域；

步骤2：局部区域划分方法：

A1、对获取的图像进行图像处理，获得岩块的尺寸和岩块中心在图像中的位置信息；

A2、局部区域划分，根据图像处理的结果，搜寻大岩块(大岩块指不合格岩块，即块度尺寸超过设计规定尺寸(L_max)，该尺寸根据实际工程确定)的位置，如当块度大于L_max为大块时，则搜索块度大于L_max的岩块并获取该岩块的位置；计算大岩块之间的距离，预设局部区域范围：以大岩块中心为圆心，半径为r(r≥L_max/2，)如果：

情况1：当大岩块均相距大于2r时，以每个大岩块中心为圆心，半径为r，分别划分单独的区域，其他区域构成一个单独区域；

情况2：当存在大岩块与大岩块距离小于2r时，如大岩块A与大岩块B，分别以岩块A与岩块B为圆心，r为半径画圆，用两圆的两条相切线L₁和L₂及两圆的半圆A₁和B₂为边界构成新的区域，若岩块A或岩块B与其他大块距离小于2r，则按照上述方法处理，并把两个区域结合起来构成一个新的区域，其余按照情况1的方法处理；

步骤3：局部区域块度统计；

分别统计爆堆顶面和坡面上各局部区域的块度，绘制局部块度分布曲线；统计爆堆顶面和坡面上的所有块度形成整体爆堆表面块度分布曲线；除了顶面和坡面，爆堆其他面的块度分布假设为整体爆堆表面块度分布曲线；

步骤S3、爆堆表面岩块三维模拟：

由于爆堆整体被简化成了梯形柱体，柱体的12条边和8个顶点分别与两个面和三个面相关，因此必须进行优先处理；而且由与由爆堆整体简化而成的梯形柱体每个面不同，因此对爆堆的6个面、12个边进行规定；

爆堆表面三维模拟的整体算法步骤如下：

B1、依据爆堆局部块度分布模型，分别随机生成8个球体记录半径，并把它随机赋予12条边，即球心位于边的端点；

B2、对12条边分别按照各自箭头的方向进行处理，使之除去端点球半径之后的剩余部分被球体填充，即球心位于边上；在随机生成球体直径时，要判断该球体位于哪个区域，然后依据爆堆局部块度分布模型随机生成球体直径；

B3、根据图3中的规定，分别提取相关边组成6个面，按照上述区域划分的结果对6个面进行区域划分，统计6个面上的球体直径分布，并把面1作为当前面；注意为了形成的封闭链，此时需依据图中信息对相关边进行方向处理，使之整个面上边缘的球体按照顺时针或者逆时针排列，并在此基础上对面上所有球体按照顺序赋值，即每个球体都有相应的序列号；

B4、随机寻找当前面中的一点作为终点，如面的中心点，并对当前面上的球体直径进行统计形成当前面上已有球体尺寸数据，并对每个区域的球体直径进行统计形成该区域的已有球体尺寸数据；

B5、计算构成初始封闭链的球体中距终点最远的球作为球1，并按照序列号寻找其两个相邻球体，把两个相邻球体中距终点最远的球作为球2；

B6、当球1和球2位于同一区域时，则依据该区域块度分布曲线和该区域已有球体尺寸数据生成新的球体的直径；当球1和球2位于不同区域时，则按照小的球体的区域块度分布曲线和区域已有球体尺寸数据生成新的球体的直径；

B7、根据球1和球2和已有的新球体的半径计算新球体的位置，注意球体球心位于面上；

B8、判断新球体与当前面上已有的球的关系，如果存在重合，则利用原先球体的直径减去某个较小的随机数形成新的直径，并返回步骤B7；直至新球体与当前面上已有的球体不重合，进行步骤B9；

B9、在新球体、球1和球2之间会有较小的空隙，可以在中间空隙处填充一个较小的球体，使之与三球体互相切；

B10、根据新球体与球1、球2相邻的球的关系更新封闭链，并把新球体直径分别加入当前面上已有球体尺寸数据和该区域已有球体尺寸数据；如果在更新封闭链时某个球体被剔除，则统计剔除的球体形成总的球体数据，返回步骤5；

B11、当前面被填满，即步骤8生成的球体半径小于某个较小的随机数，并不成功次数超过允许的尝试次数后，当前面填充流程停止；

B12、判断6个面是否都被填满，如果否，则把下个面设置成当前面，返回步骤4；如果是，则进行步骤B13；

B13、对6个面根据所在相应位置进行旋转构成三维爆堆体并进行绘图；

步骤S4、爆堆表面岩块三维模拟结果。

本申请实施例采用下述技术方案：步骤S1中，爆堆岩块形态各异且不规则，很难用一种形状去描述整体爆堆中的岩块；为了减少数值模拟的计算量，缩减计算时间，使之更适用，因此对爆堆岩块进行三维模拟时用球体代替爆堆岩块。

本申请实施例采用下述技术方案：露天深孔台阶爆破爆堆表面凹凸不平，其中爆堆的坡面近似一个斜面，爆堆顶面近似于一个抛物面，坡面和顶面存在较为明显的界限。

本申请实施例采用下述技术方案：为了模型计算的简化，因此把爆堆顶面的爆破漏斗段和坡面进行简化，把爆堆整体近似于以梯形为爆堆剖面的梯形柱体。

本申请实施例采用下述技术方案：步骤S2中，为了使爆堆表面三维模拟更加准确还原爆堆中的岩块，因此提出对爆堆表面进行区域划分，并对每个区域建立对应的块度分布曲线。

本申请实施例采用下述技术方案：步骤B1中，在相接边端点的球相同。

本申请实施例采用下述技术方案：步骤B5中，在寻找球1和球2时，如果出现多个最远球体，则在多个最远球体中随机抽取一个球作为最远球体。

本申请实施例采用下述技术方案：步骤B9中，爆堆会有许多不能识别的细小颗粒，可以对此作为一个补充。

本申请实施例采用下述技术方案：步骤S4包括步骤：对爆堆表面三维模拟的整体算法进行程序编写，并基于获得的爆堆形态主要参数和爆堆块度分布，实现了对爆堆表面的三维模拟。

本申请实施例采用的上述至少一个技术方案能够达到以下有益效果：

本发明通过把爆堆表面进行分区，并对每个区域的岩块进行统计分析，获得该区域爆堆块度级配累积曲线和大块率，构成爆堆表面块度局部分布模型；把该模型和爆堆形态信息结合起来，通过改进离散单元构造算法中的向内排列方案形成爆堆表面岩块生成算法，并采用三维可视化技术实现对爆堆形态及表面岩块三维模型模拟功能。

附图说明

此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

图1为本申请实施例中分割结果图；

图2为本申请实施例中爆堆表面局部区域划分示意图；

图3为本申请实施例中爆堆规定图；

图4为本申请实施例中爆堆表面三维模拟算法简易流程图；

图5为本申请实施例中爆堆表面三维模拟效果图；

图6为本申请实施例中爆堆表面三维模拟空隙填补效果图。

具体实施方式

为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请具体实施例及相应的附图对本申请技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

以下结合附图，详细说明本申请各实施例提供的技术方案。

实施例

一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，包括：

1.爆堆岩块形态与爆堆三维形态的确定

爆堆岩块形态各异且不规则，很难用一种形状去描述整体爆堆中的岩块。为了减少数值模拟的计算量，缩减计算时间，使之更适用，因此对爆堆岩块进行三维模拟时用球体代替爆堆岩块。露天深孔台阶爆破爆堆表面凹凸不平，其中爆堆的坡面近似一个斜面，爆堆顶面近似于一个抛物面，坡面和顶面存在较为明显的界限。为了模型计算的简化，因此把爆堆顶面的爆破漏斗段和坡面进行简化，把爆堆整体近似于以梯形为爆堆剖面的梯形柱体。

2.爆堆表面局部区域划分与块度统计

为了使爆堆表面三维模拟更加准确还原爆堆中的岩块，因此提出对爆堆表面进行区域划分，并对每个区域建立对应的块度分布曲线。按以下方法和步骤进行区域划分和区域块度统计，形成局部块度分布曲线。

步骤1：首先以梯形爆堆的6个表面进行区域划分，然后根据大块在爆堆顶面和坡面的位置及数量进一步划分局部区域。

步骤2：局部区域划分方法。

(1)对获取的图像进行图像处理，获得岩块的尺寸和岩块中心在图像中的位置信息(图1)。

(2)局部区域划分。根据图像处理的结果，搜寻大岩块(大岩块指不合格岩块，即块度尺寸超过设计规定尺寸(L_max)，该尺寸根据实际工程确定)的位置，如当块度大于L_max为大块时，则搜索块度大于L_max的岩块并获取该岩块的位置。计算大岩块之间的距离，预设局部区域范围：以大岩块中心为圆心，半径为r(r≥L_max/2，)如果

情况1：当大岩块均相距大于2r时，以每个大岩块中心为圆心，半径为r，分别划分单独的区域，该面上其他区域构成一个单独区域，如图2(a)所示。

情况2：当存在大岩块与大岩块距离小于2r时，如大岩块A与大岩块B，分别以岩块A与岩块B为圆心，r为半径画圆，用两圆的两条相切线L₁和L₂及两圆的半圆A₁和B₂为边界构成新的区域，若岩块A或岩块B与其他大块距离小于2r，则按照上述方法处理，并把两个区域结合起来构成一个新的区域。其余按照情况1的方法处理。如图2(b)所示。

步骤3：局部区域块度统计。

分别统计爆堆顶面和坡面上各局部区域的块度，形成局部块度分布曲线。

统计爆堆顶面和坡面上的所有块度形成整体爆堆表面块度分布曲线。

除了顶面和坡面，爆堆其他面的块度分布假设为整体爆堆表面块度分布曲线。

3.爆堆表面岩块三维模拟

由于爆堆整体被简化成了梯形柱体，柱体的12条边和8个顶点分别与两个面和三个面相关，因此必须进行优先处理。而且由与由爆堆整体简化而成的梯形柱体每个面不同，因此对爆堆的6个面、12个边进行规定，规定如图3所示；

如图3所示，前面为面1，左面为面2，右面为面3，后面为面4，底面为面5，顶面为面6。箭头为边的方向。图中标示的顶点为该面在二维坐标系中的(0，0)点。

爆堆表面三维模拟的整体算法步骤如下：

(1)依据爆堆局部块度分布模型，分别随机生成8个球体记录半径，并把它随机赋予12条边，即球心位于边的端点。注意在相接边端点的球相同；

(2)对12条边分别按照各自箭头的方向进行处理，使之除去端点球半径之后的剩余部分被球体填充，即球心位于边上。在随机生成球体直径时，要判断该球体位于哪个区域，然后依据爆堆局部块度分布模型随机生成球体直径；

(3)根据图3中的规定，分别提取相关边组成6个面，按照上述区域划分的结果对6个面进行区域划分，统计6个面上的球体直径分布。并把面1作为当前面。注意为了形成的封闭链，此时需依据图中信息对相关边进行方向处理，使之整个面上边缘的球体按照顺时针或者逆时针排列，并在此基础上对面上所有球体按照顺序赋值，即每个球体都有相应的序列号；

(4)随机寻找当前面中的一点作为终点，如面的中心点，并对当前面上的球体直径进行统计形成当前面上已有球体尺寸数据，并对每个区域的球体直径进行统计形成该区域的已有球体尺寸数据；

(5)计算构成初始封闭链的球体中距终点最远的球作为球1，并按照序列号寻找其两个相邻球体，把两个相邻球体中距终点最远的球作为球2。注意，在寻找球1和球2时，如果出现多个最远球体，则在多个最远球体中随机抽取一个球作为最远球体。

(6)当球1和球2位于同一区域时，则依据该区域块度分布曲线和该区域已有球体尺寸数据生成新的球体的直径；当球1和球2位于不同区域时，则按照小的球体的区域块度分布曲线和区域已有球体尺寸数据生成新的球体的直径。

(7)根据球1和球2和已有的新球体的半径计算新球体的位置。注意球体球心位于面上。

(8)判断新球体与当前面上已有的球的关系，如果存在重合，则利用原先球体的直径减去某个较小的随机数形成新的直径，并返回步骤(7)。直至新球体与当前面上已有的球体不重合，进行步骤(9)。

(9)在新球体、球1和球2之间会有较小的空隙，可以在中间空隙处填充一个较小的球体，使之与三球体互相切。(爆堆会有许多不能识别的细小颗粒，可以对此作为一个补充)。

(10)根据新球体与球1、球2相邻的球的关系更新封闭链，并把新球体直径分别加入当前面上已有球体尺寸数据和该区域已有球体尺寸数据。如果在更新封闭链时某个球体被剔除，则统计剔除的球体形成总的球体数据。返回步骤(5)。

(11)当前面被填满，即步骤(8)生成的球体半径小于某个较小的随机数，并不成功次数超过允许的尝试次数后，当前面填充流程停止。

(12)判断6个面是否都被填满，如果否，则把下个面设置成当前面，返回步骤(4)。如果是，则进行步骤(13)。

(13)对6个面根据所在相应位置进行旋转构成三维爆堆体并进行绘图。

4.爆堆表面岩块三维模拟结果

对上述算法进行程序编写，并基于获得的爆堆形态主要参数和爆堆块度分布，实现了对爆堆表面的三维模拟。图5为不同角度的爆堆表面三维模拟效果图。

图5中(a)、(b)、(c)的视角在爆堆外部，(d)的视角在爆堆内部。当对整个爆堆放大观看细节时，就会出现如图5(d)的情形，这是因为即使利用步骤(9)对部分空隙进行填充，但也会存在大部分空隙。考虑到整个爆堆是规则梯形柱体，因此对整个爆堆表面的空隙采用面的形式进行填补。填补效果如图6所示。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。

以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (9)

1.一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，包括：

步骤S1、爆堆岩块形态与爆堆三维形态的确定；

步骤S2、爆堆表面局部区域划分与块度统计：

对爆堆表面进行区域划分，并对每个区域建立对应的块度分布曲线；按以下方法和步骤进行区域划分和区域块度统计，形成局部块度分布曲线；

步骤1：首先以梯形爆堆的6个表面进行区域划分，然后根据大块在爆堆顶面和坡面的位置及数量进一步划分局部区域；

步骤2：局部区域划分方法：

A1、对获取的图像进行图像处理，获得岩块的尺寸和岩块中心在图像中的位置信息；

A2、局部区域划分，根据图像处理的结果，搜寻大岩块(大岩块指不合格岩块，即块度尺寸超过设计规定尺寸(L_max)，该尺寸根据实际工程确定)的位置，如当块度大于L_max为大块时，则搜索块度大于L_max的岩块并获取该岩块的位置；计算大岩块之间的距离，预设局部区域范围：以大岩块中心为圆心，半径为r(r≥L_max/2，)如果：

情况1：当大岩块均相距大于2r时，以每个大岩块中心为圆心，半径为r，分别划分单独的区域，其他区域构成一个单独区域；

情况2：当存在大岩块与大岩块距离小于2r时，如大岩块A与大岩块B，分别以岩块A与岩块B为圆心，r为半径画圆，用两圆的两条相切线L₁和L₂及两圆的半圆A₁和B₂为边界构成新的区域，若岩块A或岩块B与其他大块距离小于2r，则按照上述方法处理，并把两个区域结合起来构成一个新的区域，其余按照情况1的方法处理；

步骤3：局部区域块度统计；

分别统计爆堆顶面和坡面上各局部区域的块度，绘制局部块度分布曲线；统计爆堆顶面和坡面上的所有块度形成整体爆堆表面块度分布曲线；除了顶面和坡面，爆堆其他面的块度分布假设为整体爆堆表面块度分布曲线；

步骤S3、爆堆表面岩块三维模拟：

由于爆堆整体被简化成了梯形柱体，柱体的12条边和8个顶点分别与两个面和三个面相关，因此必须进行优先处理；而且由与由爆堆整体简化而成的梯形柱体每个面不同，因此对爆堆的6个面、12个边进行规定；

爆堆表面三维模拟的整体算法步骤如下：

B1、依据爆堆局部块度分布模型，分别随机生成8个球体记录半径，并把它随机赋予12条边，即球心位于边的端点；

B2、对12条边分别按照各自箭头的方向进行处理，使之除去端点球半径之后的剩余部分被球体填充，即球心位于边上；在随机生成球体直径时，要判断该球体位于哪个区域，然后依据爆堆局部块度分布模型随机生成球体直径；

B3、根据图3中的规定，分别提取相关边组成6个面，按照上述区域划分的结果对6个面进行区域划分，统计6个面上的球体直径分布，并把面1作为当前面；注意为了形成的封闭链，此时需依据图中信息对相关边进行方向处理，使之整个面上边缘的球体按照顺时针或者逆时针排列，并在此基础上对面上所有球体按照顺序赋值，即每个球体都有相应的序列号；

B4、随机寻找当前面中的一点作为终点，如面的中心点，并对当前面上的球体直径进行统计形成当前面上已有球体尺寸数据，并对每个区域的球体直径进行统计形成该区域的已有球体尺寸数据；

B5、计算构成初始封闭链的球体中距终点最远的球作为球1，并按照序列号寻找其两个相邻球体，把两个相邻球体中距终点最远的球作为球2；

B6、当球1和球2位于同一区域时，则依据该区域块度分布曲线和该区域已有球体尺寸数据生成新的球体的直径；当球1和球2位于不同区域时，则按照小的球体的区域块度分布曲线和区域已有球体尺寸数据生成新的球体的直径；

B7、根据球1和球2和已有的新球体的半径计算新球体的位置，注意球体球心位于面上；

B8、判断新球体与当前面上已有的球的关系，如果存在重合，则利用原先球体的直径减去某个较小的随机数形成新的直径，并返回步骤B7；直至新球体与当前面上已有的球体不重合，进行步骤B9；

B9、在新球体、球1和球2之间会有较小的空隙，可以在中间空隙处填充一个较小的球体，使之与三球体互相切；

B10、根据新球体与球1、球2相邻的球的关系更新封闭链，并把新球体直径分别加入当前面上已有球体尺寸数据和该区域已有球体尺寸数据；如果在更新封闭链时某个球体被剔除，则统计剔除的球体形成总的球体数据，返回步骤5；

B11、当前面被填满，即步骤8生成的球体半径小于某个较小的随机数，并不成功次数超过允许的尝试次数后，当前面填充流程停止；

B12、判断6个面是否都被填满，如果否，则把下个面设置成当前面，返回步骤4；如果是，则进行步骤B13；

B13、对6个面根据所在相应位置进行旋转构成三维爆堆体并进行绘图；

步骤S4、爆堆表面岩块三维模拟结果。

2.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤S1中，爆堆岩块形态各异且不规则，很难用一种形状去描述整体爆堆中的岩块；为了减少数值模拟的计算量，缩减计算时间，使之更适用，因此对爆堆岩块进行三维模拟时用球体代替爆堆岩块。

3.根据权利要求2所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，露天深孔台阶爆破爆堆表面凹凸不平，其中爆堆的坡面近似一个斜面，爆堆顶面近似于一个抛物面，坡面和顶面存在较为明显的界限。

4.根据权利要求3所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，为了模型计算的简化，因此把爆堆顶面的爆破漏斗段和坡面进行简化，把爆堆整体近似于以梯形为爆堆剖面的梯形柱体。

5.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤S2中，为了使爆堆表面三维模拟更加准确还原爆堆中的岩块，因此提出对爆堆表面进行区域划分，并对每个区域建立对应的块度分布曲线。

6.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤B1中，在相接边端点的球相同。

7.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤B5中，在寻找球1和球2时，如果出现多个最远球体，则在多个最远球体中随机抽取一个球作为最远球体。

8.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤B9中，爆堆会有许多不能识别的细小颗粒，可以对此作为一个补充。

9.根据权利要求1所述的一种露天深孔台阶爆破爆堆表面岩块的三维模拟方法，其特征在于，步骤S4包括步骤：对爆堆表面三维模拟的整体算法进行程序编写，并基于获得的爆堆形态主要参数和爆堆块度分布，实现了对爆堆表面的三维模拟。

Images