CN113343441A - 基于筛分法仿真的岩石爆破块度测量方法 - Google Patents
基于筛分法仿真的岩石爆破块度测量方法 Download PDFInfo
- Publication number
- CN113343441A CN113343441A CN202110558798.1A CN202110558798A CN113343441A CN 113343441 A CN113343441 A CN 113343441A CN 202110558798 A CN202110558798 A CN 202110558798A CN 113343441 A CN113343441 A CN 113343441A
- Authority
- CN
- China
- Prior art keywords
- rock
- blasting
- simulation
- screening
- measuring
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000005422 blasting Methods 0.000 title claims abstract description 44
- 238000012216 screening Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 17
- 238000009826 distribution Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 claims abstract description 4
- 239000011148 porous material Substances 0.000 claims 2
- 238000009412 basement excavation Methods 0.000 abstract description 6
- 239000002360 explosive Substances 0.000 abstract description 6
- 238000005065 mining Methods 0.000 abstract description 4
- 238000007619 statistical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 102100021807 ER degradation-enhancing alpha-mannosidase-like protein 1 Human genes 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000005034 decoration Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000009440 infrastructure construction Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于筛分法仿真的岩石爆破块度测量方法,包括如下步骤:1)将岩石样本平铺在空地上,取样量拟按毎次爆破取3~5组,每组的取样量不小于1500kg;2)对每一组岩石样本进行扫描,获得岩石样本的精细化3D‑CAD数据;3)将精细化3D‑CAD数据导入离散元软件中进行3D重构,并对筛分过程进行仿真模拟,获取岩石样本的块度分布;4)多次仿真模拟试验,利用作图法对筛分统计结果进行处理,得出爆堆岩石的块度分布曲线、不均匀系数以及曲率系数。通过对爆堆岩石样本扫描并采用离散元软件对筛分过程进行模拟仿真,解决了现场筛分面临的工作量大的问题,降低了成本,避免了对现场的正常采矿、开挖作业的影响,同时具有精度高等优点。
Description
技术领域
本发明涉及工程爆破技术领域,特别是涉及基于筛分法仿真的岩石爆破块度测量方法。
背景技术
爆破开挖是水利水电工程、矿山、交通等领域中最基本的施工过程,伴随着岩石爆破开挖的进行,大量爆破碎石形成。矿岩爆破效果的好坏直接影响了矿山生产、大坝开挖建设、基础设施建设过程中爆破碎石铲装、运输、二次破碎等后续工艺的效率,同时也对评价破碎方法、研究破碎机理、决定开采方案以及选矿、放矿等至关重要的问题提供指导,为爆破参数调整、网络优化、装药结构改善提供依据。
通常可以用爆堆形状、岩块块度、有无飞石、爆破震动、爆破成本等指标对矿岩爆破效果进行衡量,其中爆堆的块度分布尤为重要。评价岩石爆破块度的方法基本上可划分为直接统计分析法和间接统计分析法。
直接统计分析法包括应用普遍的筛分法、二次爆破岩块统计法、直接测量法。筛分法适用于筛分量小的情况,虽然现场筛分法具有精度高、直观等优点,但由于现场筛分工作量大,容易导致正常采矿、开挖作业受到干扰,成本增加。
发明内容
本发明的目的在于至少解决现有技术中存在的技术问题之一,为此,本发明提出基于筛分法仿真的岩石爆破块度测量方法,包括如下步骤:
1)爆破作业后,采用反铲在爆堆中心取样,并用自卸车装运至一宽敞的空地,将岩石样本平铺在空地上,每个岩石样本之间留有间隙,并保证岩石样本无重叠现象,取样量拟按毎次爆破取3~5组,每组的取样量不小于1500kg;
2)利用高精度三维扫描测量仪对每一组岩石样本进行扫描,获得每一组岩石样本的精细化3D-CAD数据;
3)将岩石样本的精细化3D-CAD数据导入离散元软件中,对岩石样本进行3D重构,获得岩石样本的表面积、比表面积、体积以及质量参数,并对筛分过程进行仿真模拟,获取岩石样本的块度分布;
4)进行多次仿真模拟试验,利用作图法对筛分统计结果进行处理,得出爆堆岩石的块度分布曲线、不均匀系数以及曲率系数。
本发明的有益效果是,通过对爆堆岩石样本扫描并采用离散元软件对筛分过程进行模拟仿真,解决了现场筛分面临的工作量大的问题,降低了成本,避免了对现场的正常采矿、开挖作业的影响,同时具有精度高等优点。
优选地,步骤2)中,采用的离散元软件为3DEC软件。
优选地,步骤3)中,模拟筛分的过程中,构建出筛网的三维模型,筛网具有多个筛孔。
优选地,步骤3)中,构建多个不同尺寸筛孔的筛网的三维模型。
优选地,步骤3)中,筛孔形状为圆形。
优选地,步骤3)中,筛孔的半径大于等于1mm。
优选地,所述步骤3)具体为:将步骤2)导入重构的岩石样本的三维模型加载动力载荷,通过步骤3)中的多个不同尺寸筛孔的筛网分成若干个粒级,根据现场测得的岩石密度,对仿真过程中的岩石模型进行赋值与标定,得到每个粒级的质量百分数。
附图说明
下面结合附图对本发明作进一步说明,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明方法的流程图;
图2是本发明方法具体实施采用的设备及结果处理的流程示意图;
图3是发明具体实施方式中筛孔形状为正方形的多个筛网层叠的示意图;
图4是图3多个不同筛孔尺寸的筛网的示意图;
图5是发明具体实施方式中筛孔形状为圆形的多个筛网层叠的示意图;
图6是图5多个不同筛孔尺寸的筛网的示意图;
图7是已获取相关具体数据的爆堆岩石样品块度分布曲线示意图。
具体实施方式
为了使本领域技术人员更好地理解本发明的技术方案,下面结合附图对本发明进行详细描述,本部分的描述仅是示范性和解释性,不应对本发明的保护范围有任何的限制作用。
需要理解的是,本文中,若有术语“上”、“下”、“左”、“右”、“前”、“后”、“底”、“顶”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。本文中,若有术语“第一”、“第二”等是用于区别不同对象,而不是用于描述特定顺序。
需要说明的是,在本文中,若有术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。本发明的描述中,除非另有明确的限定,设置、安装、连接等词语应做广义理解,所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。本文中,如果有描述到“若干”、“多个”、“多次”,若干的含义是一个或者多个,多个的含义是两个以上,多次的含义是两次以上,大于、小于、超过等理解为不包括本数,以上、以下、以内等理解为包括本数。
本部分将详细描述本发明的具体实施例,本发明之较佳实施例在附图中示出,附图的作用在于用图形补充说明书文字部分的描述,使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案,但其不能理解为对本发明保护范围的限制。
参照图1~图7,本发明实施例一种基于筛分法仿真的岩石爆破块度测量方法,包括如下步骤:
1)爆破作业后,采用反铲在爆堆100中心取样,并用自卸车装运至一宽敞的空地,将岩石样本200平铺在空地上,每个岩石样本200之间留有间隙,并保证岩石样本200无重叠现象,取样量拟按毎次爆破取3~5组;
2)利用高精度三维扫描测量仪300对每一组岩石样本200进行扫描,获得每一组岩石样本200的精细化3D-CAD数据;
3)将岩石样本200的精细化3D-CAD数据导入计算机400的离散元软件中,对岩石样本200进行3D重构,获得岩石样本200的表面积、比表面积、体积以及质量等参数,并对筛分过程进行仿真模拟,获取岩石样本200的块度分布;
4)进行多次仿真模拟试验,利用作图法对筛分统计结果进行处理,得出爆堆岩石的块度分布曲线、不均匀系数以及曲率系数等信息,制作爆堆岩石块度分布曲线图500以及表面积、比表面积、体积等相关参数信息的汇总表600,例如图7所示的为已获取相关具体数据的爆堆岩石样品块度分布曲线示意图。
采用高精度三维扫描仪300对爆堆岩石进行扫描,获得爆堆岩石样本200的全周真实的轮廓及几何形态,不需要解决采用平面摄影所面临的将岩块形体的二维数据转化为三维数据的问题。
优选地,采用的离散元软件为3DEC软件。在本发明其他一些的实施例中,离散元软件为EDEM软件等其他任意软件类型。
继续参照图3至图6,进一步地,模拟筛分的过程中,构建出筛网700的三维模型,筛网700具有多个筛孔710。在本发明一些实施例中,构建多个不同尺寸的筛孔710的筛网700的三维模型。
继续参照图5和图6,优选地,筛孔710形状为圆形。优选地,圆形的筛孔710的半径大于等于1mm,具体的半径尺寸例如为1mm、1.5mm、2mm等
继续参照图3和图4,在本发明其他一些的实施例中,筛孔710形状为正方形。进一步优选地,筛孔710的边长大于等于1mm,具体的边长尺寸例如为1mm、1.5mm、2mm等。
需要理解的是,本领域技术人员有能力根据现行技术规范的指导,设定选择符合相关规范的筛孔形状和尺寸大小。
优选地,所述步骤3)具体为:将步骤2)导入重构的岩石样本200的三维模型加载动力载荷,通过步骤3)中的多个不同尺寸筛孔710的筛网700分成若干个粒级,根据现场测得的岩石密度,对仿真过程中的岩石模型进行赋值与标定,得到每个粒级的质量百分数。
本方法利用离散元软件对筛分过程进行仿真模拟,不仅得到岩石样本200的块度分布,而且也得到岩石样本200的表面积、比表面积、体积、质量等参数。
本方法既既保留了采用现场筛分法精度高、准确度高、直观等优点,又解决了现场筛分法测量岩石爆破块度所面临的工作量大,成本高等问题。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实例的说明只是用于帮助理解本发明的方法及其核心思想。以上仅是本发明的优选实施方式,应当指出,由于文字表达的有限性,而客观上存在无限的具体结构,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进、润饰或变化,也可以将上述技术特征以适当的方式进行组合;这些改进润饰、变化或组合,或未经改进将发明的构思和技术方案直接应用于其它场合的,均应视为本发明的保护范围。
Claims (7)
1.基于筛分法仿真的岩石爆破块度测量方法,其特征在于,包括如下步骤:
1)爆破作业后,采用反铲在爆堆中心取样,并用自卸车装运至一宽敞的空地,将岩石样本平铺在空地上,每个岩石样本之间留有间隙,并保证岩石样本无重叠现象,取样量拟按毎次爆破取3~5组,每组的取样量不小于1500kg;
2)利用高精度三维扫描测量仪对每一组岩石样本进行扫描,获得每一组岩石样本的精细化3D-CAD数据;
3)将岩石样本的精细化3D-CAD数据导入离散元软件中,对岩石样本进行3D重构,获得岩石样本的表面积、比表面积、体积以及质量参数,并对筛分过程进行仿真模拟,获取岩石样本的块度分布;
4)进行多次仿真模拟试验,利用作图法对筛分统计结果进行处理,得出爆堆岩石的块度分布曲线、不均匀系数以及曲率系数。
2.根据权利要求1所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:步骤2)中,采用的离散元软件为3DEC软件。
3.根据权利要求1或2所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:步骤3)中,模拟筛分的过程中,构建出筛网的三维模型,筛网具有多个筛孔。
4.根据权利要求3所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:步骤3)中,构建多个不同尺寸筛孔的筛网的三维模型。
5.根据权利要求4所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:步骤3)中,筛孔形状为圆形。
6.根据权利要求5所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:步骤3)中,筛孔的半径大于等于1mm。
7.根据权利要求6所述的基于筛分法仿真的岩石爆破块度测量方法,其特征在于:所述步骤3)具体为:将步骤2)导入重构的岩石样本的三维模型加载动力载荷,通过步骤3)中的多个不同尺寸筛孔的筛网分成若干个粒级,根据现场测得的岩石密度,对仿真过程中的岩石模型进行赋值与标定,得到每个粒级的质量百分数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110558798.1A CN113343441B (zh) | 2021-05-21 | 2021-05-21 | 基于筛分法仿真的岩石爆破块度测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110558798.1A CN113343441B (zh) | 2021-05-21 | 2021-05-21 | 基于筛分法仿真的岩石爆破块度测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113343441A true CN113343441A (zh) | 2021-09-03 |
CN113343441B CN113343441B (zh) | 2022-10-04 |
Family
ID=77470634
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110558798.1A Active CN113343441B (zh) | 2021-05-21 | 2021-05-21 | 基于筛分法仿真的岩石爆破块度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113343441B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239637A (zh) * | 2014-09-16 | 2014-12-24 | 武汉大学 | 一种离散元爆堆形态模拟方法 |
CN106227913A (zh) * | 2016-07-06 | 2016-12-14 | 四川大学 | 用于水电工程的堆石料爆破级配预测方法 |
CN107489420A (zh) * | 2017-08-11 | 2017-12-19 | 辽宁科技大学 | 无底柱分段崩落法矿岩块度分布远程实时监测系统及方法 |
CN108596881A (zh) * | 2018-04-09 | 2018-09-28 | 中国科学院武汉岩土力学研究所 | 一种爆破后岩石粒度的智能图像统计方法 |
CN109030305A (zh) * | 2018-06-08 | 2018-12-18 | 武汉大学 | 一种基于三维激光扫描和图像处理技术的级配料比表面积确定方法 |
CN110068573A (zh) * | 2019-04-22 | 2019-07-30 | 武汉大学 | 基于摄影测量和颗粒接触理论的爆破块度测量方法 |
US20200089823A1 (en) * | 2018-03-26 | 2020-03-19 | OrePro Holdings LLC | 3d block modelling of a resource boundary in a post-blast muckpile to optimize destination delineation |
CN111814372A (zh) * | 2020-07-02 | 2020-10-23 | 大昌建设集团有限公司 | 一种联合数值计算和块度筛分的爆破块度控制方法 |
CN113344851A (zh) * | 2021-04-28 | 2021-09-03 | 鞍钢矿业爆破有限公司 | 一种摄影法测量爆堆修正函数数据方法 |
-
2021
- 2021-05-21 CN CN202110558798.1A patent/CN113343441B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239637A (zh) * | 2014-09-16 | 2014-12-24 | 武汉大学 | 一种离散元爆堆形态模拟方法 |
CN106227913A (zh) * | 2016-07-06 | 2016-12-14 | 四川大学 | 用于水电工程的堆石料爆破级配预测方法 |
CN107489420A (zh) * | 2017-08-11 | 2017-12-19 | 辽宁科技大学 | 无底柱分段崩落法矿岩块度分布远程实时监测系统及方法 |
US20200089823A1 (en) * | 2018-03-26 | 2020-03-19 | OrePro Holdings LLC | 3d block modelling of a resource boundary in a post-blast muckpile to optimize destination delineation |
CN108596881A (zh) * | 2018-04-09 | 2018-09-28 | 中国科学院武汉岩土力学研究所 | 一种爆破后岩石粒度的智能图像统计方法 |
CN109030305A (zh) * | 2018-06-08 | 2018-12-18 | 武汉大学 | 一种基于三维激光扫描和图像处理技术的级配料比表面积确定方法 |
CN110068573A (zh) * | 2019-04-22 | 2019-07-30 | 武汉大学 | 基于摄影测量和颗粒接触理论的爆破块度测量方法 |
CN111814372A (zh) * | 2020-07-02 | 2020-10-23 | 大昌建设集团有限公司 | 一种联合数值计算和块度筛分的爆破块度控制方法 |
CN113344851A (zh) * | 2021-04-28 | 2021-09-03 | 鞍钢矿业爆破有限公司 | 一种摄影法测量爆堆修正函数数据方法 |
Non-Patent Citations (7)
Title |
---|
ONEDERRA等: "Measuring blast fragmentation at Esperanza mine using high-resolution 3D laser scanning", 《MINING TECHNOLOGY》 * |
XINPING LI等: "Numerical Simulation of Blast Vibration and Crack Forming Effect of Rock-Anchored Beam Excavation in Deep Underground Caverns", 《HINDAWI》 * |
吴亮等: "矿岩爆破破碎机理、块度分布与测量技术研究动态", 《金属矿山》 * |
朱明秀: "第九届国际岩石爆破破岩学术会议论文集摘要(二)", 《工程爆破》 * |
林大泽: "爆堆块度评价方法研究的进展", 《中国安全科学学报》 * |
王鑫瑀等: "岩石爆破破碎块度分布测量方法综述", 《化工矿物与加工》 * |
荆永滨等: "基于块体形状分类的爆堆块度三维几何参数估算", 《爆破》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113343441B (zh) | 2022-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kong et al. | Can 3-D printed gypsum samples replicate natural rocks? An experimental study | |
Elmo et al. | An integrated numerical modelling–discrete fracture network approach applied to the characterisation of rock mass strength of naturally fractured pillars | |
Sterpi et al. | A physical and numerical investigation on the stability of shallow tunnels in strain softening media | |
US10656068B2 (en) | Method for providing a numerical model of a sample of rock | |
CN102023310B (zh) | 三维观测系统均匀性分析方法 | |
CN105913494B (zh) | 多尺度裂缝精细地质建模及数值模拟方法和装置 | |
Gupta | Effects of particle size and confining pressure on breakage factor of rockfill materials using medium triaxial test | |
Liu et al. | Microscale crack propagation in shale samples using focused ion beam scanning electron microscopy and three-dimensional numerical modeling | |
CN106447776A (zh) | 基于3d打印制作的复杂裂隙岩体物理模型及建模方法 | |
CN103235109B (zh) | 一种低重力环境下模拟月壤力学性能的测定方法 | |
CN103940669B (zh) | 一种固体充填材料力学特性测试方法 | |
Jiang et al. | Influence of geometrical distribution of rock joints on deformational behavior of underground opening | |
CN105781543A (zh) | 一种基于裂隙分维指标的采动覆岩质量评价方法 | |
CN114330144B (zh) | 一种边坡危石预警方法、装置、设备及可读存储介质 | |
SA520420062B1 (ar) | اختبارات تثلُّم نانو متري لتمييز الصدوع الهيدروليكية | |
CN110318743A (zh) | 薄互层页岩油储层压裂模拟试验方法及装置 | |
Ray et al. | Characterizing and modeling natural fracture networks in a tight carbonate reservoir in the Middle East: A methodology | |
CN108090283A (zh) | 一种列车动载下隧道—车辆耦合振动的有限元分析方法 | |
CN113343441B (zh) | 基于筛分法仿真的岩石爆破块度测量方法 | |
CN104866682B (zh) | 一种基于地面形变进行页岩气勘探区构造应力场反演的方法 | |
CN103487304A (zh) | 一种土石混合体重塑试样制作方法 | |
Hao et al. | Mesoscopic damage evolution characteristics of jointed sandstone under different loading conditions | |
Sharif et al. | An Investigation of the Factors Controlling the Mechanical Behaviour of Slender Naturally Fractured Pillars | |
CN116305451B (zh) | 连续-非连续地质模型建立方法及装置 | |
CN116226982A (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 |