CN108256253A - 一种基于Voronoi网格剖分的等面积炮孔布置方法 - Google Patents

一种基于Voronoi网格剖分的等面积炮孔布置方法 Download PDF

Info

Publication number
CN108256253A
CN108256253A CN201810082888.6A CN201810082888A CN108256253A CN 108256253 A CN108256253 A CN 108256253A CN 201810082888 A CN201810082888 A CN 201810082888A CN 108256253 A CN108256253 A CN 108256253A
Authority
CN
China
Prior art keywords
voronoi
region
target
target demolition
blasthole
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
Application number
CN201810082888.6A
Other languages
English (en)
Other versions
CN108256253B (zh
Inventor
周玉县
郑善喜
黄晓锋
徐清清
李超珍
戚鹏英
邓会芳
张军贤
裴亚兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nuclear Industry Huzhou Survey Planning Design And Research Institute Co ltd
Original Assignee
Nuclear Industry Huzhou Engineering Investigation Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nuclear Industry Huzhou Engineering Investigation Institute filed Critical Nuclear Industry Huzhou Engineering Investigation Institute
Priority to CN201810082888.6A priority Critical patent/CN108256253B/zh
Publication of CN108256253A publication Critical patent/CN108256253A/zh
Application granted granted Critical
Publication of CN108256253B publication Critical patent/CN108256253B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling

Abstract

本发明涉及一种基于Voronoi网格剖分的等面积炮孔布置方法,首先针对目标爆破区域进行三角形网格划分,然后针对目标爆破区域经三角形网格划分后的所有节点,获得目标爆破区域中的各个Voronoi多边形,最后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置;整个设计方法,只需要对非规则爆破区域进行简单的二维三角形网格划分,即可构造出Voronoi多边形区域,能够获得合理的炮孔布置方案,在爆破范围内均匀布置炮孔,保证每个炮孔所控制的爆破区域相等,大大提高了炸药爆炸能量的利用率,改善爆破块度的均匀性,以及爆堆的松散度。

Description

一种基于Voronoi网格剖分的等面积炮孔布置方法
技术领域
本发明涉及一种基于Voronoi网格剖分的等面积炮孔布置方法,属于岩土工程、爆破工程施工技术领域。
背景技术
在岩土工程施工中,爆破是土石方开挖最直接最有效和使用最为广泛的手段,而爆破中炮孔的布置形式在很大程度上决定了爆破效果的好坏,从而影响着后续一系列生产效率。传统的炮孔位置标定方法是利用采集到的爆破区域的参数,在现场进行炮孔位置标定,由于现场实际条件的有限性及复杂性,传统的炮孔布置方法存在很多弊端,例如:炮孔位置标定更多地凭经验完成,炮孔的间距和深度只能粗略估计等。
发明内容
本发明所要解决的技术问题是提供一种基于Voronoi网格剖分的等面积炮孔布置方法,能够获得合理的炮孔布置方案,大大提高了炸药爆炸能量的利用率,改善爆破块度的均匀性,以及爆堆的松散度。
本发明为了解决上述技术问题采用以下技术方案:本发明设计了一种基于Voronoi网格剖分的等面积炮孔布置方法,用于针对目标爆破区域实现炮孔分布的布置,包括如下步骤:
步骤A.利用二维Delaunay化,按边长为a的正三角形,针对目标爆破区域进行三角形网格划分,然后进入步骤B;
步骤B.针对目标爆破区域经三角形网格划分后的所有节点,由最外围的任意一个节点起,按顺时针方向或逆时针方向依次针对各个节点做如下操作:
针对该节点所在的各个正三角网格区域,将该各个正三角网格区域的中心点位置,按顺时针方向或逆时针方向依次串联,构成Voronoi多边形;
完成针对目标爆破区域内所有节点的操作后,获得目标爆破区域中的各个Voronoi多边形,然后进入步骤C;
步骤C.获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
作为本发明的一种优选技术方案:所述步骤A中,利用二维Delaunay化,按预设边长为a的正三角形,针对目标爆破区域进行三角形网格划分。
作为本发明的一种优选技术方案:所述步骤C之后,还包括步骤D如下,执行完步骤C之后,进入步骤D;
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
作为本发明的一种优选技术方案:所述步骤A中,根据目标爆破区域面积S,以及目标爆破区域预设炮孔数目n,按如下公式:
获得针对目标爆破区域进行三角形网格划分中、正三角形的边长a;
所述步骤C之后,还包括步骤D如下,执行完步骤C之后,进入步骤D;
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
作为本发明的一种优选技术方案:在获得目标爆破区域中各个炮孔开挖位置后,根据预设炮孔深度H,按如下公式:
获得各个炮孔所对应的岩石爆破体积V0,进而根据V0,获取各个炮孔的炸药填装量。
本发明所述一种基于Voronoi网格剖分的等面积炮孔布置方法采用以上技术方案与现有技术相比,具有以下技术效果:本发明设计的基于Voronoi网格剖分的等面积炮孔布置方法,只需要对非规则爆破区域进行简单的二维三角形网格划分,即可构造出Voronoi多边形区域,能够获得合理的炮孔布置方案,在爆破范围内均匀布置炮孔,保证每个炮孔所控制的爆破区域相等,大大提高了炸药爆炸能量的利用率,改善爆破块度的均匀性,以及爆堆的松散度。
附图说明
图1是爆破区域结构示意图;
图2是利用二维Delaunay化生成Voronoi多边形原理图;
图3是实例中的对爆破区域进行三角形网格划分;
图4是实例中的利用Voronoi网格生成的炮孔均匀布置图。
具体实施方式
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。
Voronoi网格技术是19世纪50年代,Dirichlet提出的将空间S分解为一系列Voronoi凸多面体的理论。对于空间S内任意给定点集{Pi},i=1,2,…,k,可将空间S分解为一组凸几何体{Vi},i=1,2,…,k,其中,Vi包含Pi,称之为Pi的领域,Pi称为领域Vi的特征点,凸几何体结构称为Dirichlet镶嵌,又称为Voronoi图。Voronoi图满足如下性质,领域Vi靠近Pi最密切,即:
Vi={P:d(P,Pi)<d(P,Pj),i≠j,P∈S}
其中,d(P,Pi)、d(P,Pj)分别表示P与Pi、P与Pj之间的欧几里得距离。
鉴于此,本发明所设计一种基于Voronoi网格剖分的等面积炮孔布置方法,能克服传统的炮孔位置布置方法的弊端,在很大程度上改善爆破效果,提高生产效率。
实际应用中,本发明设计了一种基于Voronoi网格剖分的等面积炮孔布置方法,用于针对目标爆破区域实现炮孔分布的布置,实际应用当中,设定炮孔间、排距大致相等,即要求每个炮孔所控制的平面面积相等;为了保证第一排炮孔爆破时自由面小于间距,将自由面线向已开挖区域法向平移一定距离(炮孔间距的20%),如图1所示。并且就实际应用,具体可以划分为两套方案,其中,第一套方案具体包括如下步骤:
步骤A.利用二维Delaunay化,按预设边长为a的正三角形,针对目标爆破区域进行三角形网格划分,然后进入步骤B。
步骤B.针对目标爆破区域经三角形网格划分后的所有节点,由最外围的任意一个节点起,按顺时针方向或逆时针方向依次针对各个节点做如下操作:
针对该节点所在的各个正三角网格区域,将该各个正三角网格区域的中心点位置,按顺时针方向或逆时针方向依次串联,构成Voronoi多边形;
完成针对目标爆破区域内所有节点的操作后,获得目标爆破区域中的各个Voronoi多边形,然后进入步骤C。
步骤C.获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置,然后进入步骤D。
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
第二套方案具体包括如下步骤:
步骤A.首先根据目标爆破区域面积S,以及目标爆破区域预设炮孔数目n,按如下公式:
获得针对目标爆破区域进行三角形网格划分中、正三角形的边长a;
然后利用二维Delaunay化,按边长为a的正三角形,针对目标爆破区域进行三角形网格划分,然后进入步骤B。
步骤B.针对目标爆破区域经三角形网格划分后的所有节点,由最外围的任意一个节点起,按顺时针方向或逆时针方向依次针对各个节点做如下操作:
针对该节点所在的各个正三角网格区域,将该各个正三角网格区域的中心点位置,按顺时针方向或逆时针方向依次串联,构成Voronoi多边形;
完成针对目标爆破区域内所有节点的操作后,获得目标爆破区域中的各个Voronoi多边形,然后进入步骤C。
步骤C.获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置,然后进入步骤D。
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
基于上述分别所设计的两套技术方案,即可获得目标爆破区域中各个炮孔开挖位置,之后根据预设炮孔深度H,按如下公式:
即可获得各个炮孔所对应的岩石爆破体积V0,进而根据V0,获取各个炮孔的炸药填装量,完成针对目标爆破区域的爆破。
将上述所设计第一套方案应用到实际实施例中,如下:
以某露天矿山为例,岩石为白云岩矿石,硬度大,可爆性差,台阶高度12m,孔深14.5m,超深2.5m,炮孔直径170mm作为实施例,测量整个目标爆破区域沿长度方向尺寸为55.2m,沿宽度方向尺寸为36.8m,整个目标爆破区域面积为2031.36m2,设置炮孔间距参数a为2.0m,并以此参数a作为网格划分的单元尺寸。
步骤A.如图2所示,利用二维Delaunay化,按预设边长为a=2.0m的正三角形,针对目标爆破区域进行三角形网格划分,然后进入步骤B。
步骤B.针对目标爆破区域经三角形网格划分后的所有节点,由最外围的任意一个节点起,按顺时针方向或逆时针方向依次针对各个节点做如下操作:
针对该节点所在的各个正三角网格区域,将该各个正三角网格区域的中心点位置,按顺时针方向或逆时针方向依次串联,构成Voronoi多边形;
完成针对目标爆破区域内所有节点的操作后,获得目标爆破区域中的各个Voronoi多边形,然后进入步骤C。即如图4所示,计算得到Voronoi多边形的边长为1.4m,Voronoi多边形节点数目np为610、单元数目ne为600,其中三角形单元数目为89、四边形数目为0、五边形数目为20、六边形数目为473、七边形数目为18,并记录每个节点的坐标和每个单元形心坐标。
步骤C.获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置,即如图4所示,Voronoi多边形的中心坐标就是要布置炮孔的位置,然后进入步骤D。
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
如图3所示,此时,每个炮孔控制的爆破面积假定炮孔深度为H=14.5m,通过计算可得每个炮孔所爆除的岩石体积以该体积选择单孔炸药装药量。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (5)

1.一种基于Voronoi网格剖分的等面积炮孔布置方法,用于针对目标爆破区域实现炮孔分布的布置,其特征在于,包括如下步骤:
步骤A.利用二维Delaunay化,按边长为a的正三角形,针对目标爆破区域进行三角形网格划分,然后进入步骤B;
步骤B.针对目标爆破区域经三角形网格划分后的所有节点,由最外围的任意一个节点起,按顺时针方向或逆时针方向依次针对各个节点做如下操作:
针对该节点所在的各个正三角网格区域,将该各个正三角网格区域的中心点位置,按顺时针方向或逆时针方向依次串联,构成Voronoi多边形;
完成针对目标爆破区域内所有节点的操作后,获得目标爆破区域中的各个Voronoi多边形,然后进入步骤C;
步骤C.获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
2.根据权利要求1所述一种基于Voronoi网格剖分的等面积炮孔布置方法,其特征在于:所述步骤A中,利用二维Delaunay化,按预设边长为a的正三角形,针对目标爆破区域进行三角形网格划分。
3.根据权利要求2所述一种基于Voronoi网格剖分的等面积炮孔布置方法,其特征在于:所述步骤C之后,还包括步骤D如下,执行完步骤C之后,进入步骤D;
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
4.根据权利要求1所述一种基于Voronoi网格剖分的等面积炮孔布置方法,其特征在于:所述步骤A中,根据目标爆破区域面积S,以及目标爆破区域预设炮孔数目n,按如下公式:
获得针对目标爆破区域进行三角形网格划分中、正三角形的边长a;
所述步骤C之后,还包括步骤D如下,执行完步骤C之后,进入步骤D;
步骤D.针对目标爆破区域中的各个Voronoi多边形,删除其中边数小于4个的Voronoi多边形,更新目标爆破区域中的各个Voronoi多边形,然后获得目标爆破区域中各个Voronoi多边形的中心点位置,即为目标爆破区域中各个炮孔开挖位置。
5.根据权利要求1至4中任意一项所述一种基于Voronoi网格剖分的等面积炮孔布置方法,其特征在于,在获得目标爆破区域中各个炮孔开挖位置后,根据预设炮孔深度H,按如下公式:
获得各个炮孔所对应的岩石爆破体积V0,进而根据V0,获取各个炮孔的炸药填装量。
CN201810082888.6A 2018-01-29 2018-01-29 一种基于Voronoi网格剖分的等面积炮孔布置方法 Active CN108256253B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810082888.6A CN108256253B (zh) 2018-01-29 2018-01-29 一种基于Voronoi网格剖分的等面积炮孔布置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810082888.6A CN108256253B (zh) 2018-01-29 2018-01-29 一种基于Voronoi网格剖分的等面积炮孔布置方法

Publications (2)

Publication Number Publication Date
CN108256253A true CN108256253A (zh) 2018-07-06
CN108256253B CN108256253B (zh) 2022-05-24

Family

ID=62743017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810082888.6A Active CN108256253B (zh) 2018-01-29 2018-01-29 一种基于Voronoi网格剖分的等面积炮孔布置方法

Country Status (1)

Country Link
CN (1) CN108256253B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110047141A (zh) * 2019-03-02 2019-07-23 长沙迪迈数码科技股份有限公司 露天矿山开采现状三维模型更新方法
CN113139319A (zh) * 2021-05-14 2021-07-20 辽宁工程技术大学 一种基于智能岩性识别的炮孔装药量计算方法
CN114508978A (zh) * 2022-01-28 2022-05-17 北京理工大学 炮孔装药量的计算方法、系统及电子设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243680A (zh) * 2011-07-21 2011-11-16 中国科学技术大学 一种网格划分方法及系统
CN102663183A (zh) * 2012-03-31 2012-09-12 浙江大学 数字矿山中的爆破仿真方法
CN103440683A (zh) * 2013-04-28 2013-12-11 大连大学 一种基于三维散乱稠密点云的三角网格重构方法
CN106652036A (zh) * 2016-11-29 2017-05-10 中国科学院自动化研究所 二维三角网格生成方法及系统
US20170132174A1 (en) * 2015-11-06 2017-05-11 Jeffrey W. Holcomb Novel method for the fast derivation of delaunay tesselations
CN107172626A (zh) * 2017-04-18 2017-09-15 中国矿业大学 一种co2地质封存区域的传感网络节点定位优化方法
CN107391837A (zh) * 2017-07-19 2017-11-24 中南大学 一种扇形炮孔的孔口装药结构参数优化计算方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243680A (zh) * 2011-07-21 2011-11-16 中国科学技术大学 一种网格划分方法及系统
CN102663183A (zh) * 2012-03-31 2012-09-12 浙江大学 数字矿山中的爆破仿真方法
CN103440683A (zh) * 2013-04-28 2013-12-11 大连大学 一种基于三维散乱稠密点云的三角网格重构方法
US20170132174A1 (en) * 2015-11-06 2017-05-11 Jeffrey W. Holcomb Novel method for the fast derivation of delaunay tesselations
CN106652036A (zh) * 2016-11-29 2017-05-10 中国科学院自动化研究所 二维三角网格生成方法及系统
CN107172626A (zh) * 2017-04-18 2017-09-15 中国矿业大学 一种co2地质封存区域的传感网络节点定位优化方法
CN107391837A (zh) * 2017-07-19 2017-11-24 中南大学 一种扇形炮孔的孔口装药结构参数优化计算方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
LEONIDAS J. GUIBAS 等: "Randomized incremental construction of Delaunay and Voronoi diagrams", 《ALGORITHMICA》 *
佘红伟: "二维区域网格剖分算法研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 基础科学辑》 *
余东晓 等: "基于voronoi网格技术的逐孔台阶爆破起爆顺序确定方法", 《第一届全国工程安全与防护学术会议》 *
余东晓 等: "逐孔台阶爆破设计与优化软件开发", 《爆破》 *
倪国喜等: "任意区域上的粒子均匀分布方法", 《计算力学学报》 *
卓中文 等: "基于Delaunay三角网生成Voronoi爆破图的研究", 《现代矿业》 *
朱培烨等: "Delaunay非结构网格生成之布点技术", 《航空计算技术》 *
王盛玺等: "基于约束Delaunay三角化的二维非结构网格生成方法", 《计算物理》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110047141A (zh) * 2019-03-02 2019-07-23 长沙迪迈数码科技股份有限公司 露天矿山开采现状三维模型更新方法
CN113139319A (zh) * 2021-05-14 2021-07-20 辽宁工程技术大学 一种基于智能岩性识别的炮孔装药量计算方法
CN114508978A (zh) * 2022-01-28 2022-05-17 北京理工大学 炮孔装药量的计算方法、系统及电子设备

Also Published As

Publication number Publication date
CN108256253B (zh) 2022-05-24

Similar Documents

Publication Publication Date Title
CN108256253A (zh) 一种基于Voronoi网格剖分的等面积炮孔布置方法
CN101344390B (zh) 地形高程的高保真数字建模方法
CN102915560B (zh) 一种阈值无关的机载激光雷达点云滤波方法
CN104156523B (zh) 索网反射面天线的测地线索网生成方法
CN106895755B (zh) 一种露天台阶爆破智能化设计方法
CN105654552B (zh) 一种面向任意分布大规模点云数据的快速Delaunay构网方法
CN111692931B (zh) 一种隧道爆破智能布孔方法
CN107391837B (zh) 一种扇形炮孔的孔口装药结构参数优化计算方法
CN107989612B (zh) 一种高潜水位煤矿区动态预复垦新方法
CN103092933A (zh) 面向海量点云数据的基于矩形拼合的Delaunay三角网并行构网方法
CN112685944B (zh) 一种同类多枚杀爆弹对面目标瞄准点规划方法
CN105474824A (zh) 一种耕作层土壤的剥离及储存方法
CN109555556A (zh) 一种面向耕地保护的固体充填采煤充实率设计方法
CN104322173A (zh) 山地丘陵区缓坡地改造方法
CN106980716A (zh) 一种基于随机块体的地下洞室稳定分析方法
CN106969676B (zh) 一种基于rtk炮孔定位技术改善爆破效果的布孔微调方法
CN102997768B (zh) 一种爆破岩石的方法
CN103823946B (zh) 一种模拟带有不整合面岩体的方法
CN104318594A (zh) 一种平面向量场的四边形网格化方法
Liao et al. Grid partition variable step alpha shapes algorithm
CN102385757A (zh) 一种基于几何空间的语义约束纹理合成方法
CN105977119B (zh) 一种栅控电子枪虚阳极位置的确定方法
CN103444348A (zh) 一种马蹄秧苗及表层泥土去除结构
CN105898779A (zh) 采用可信信息覆盖模型检测无线传感网覆盖空洞的方法
CN206573018U (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
CP03 Change of name, title or address

Address after: No. 581, Fenghuang Road, Huzhou City, Zhejiang Province, 313000

Patentee after: Nuclear industry Huzhou survey, planning, design and Research Institute Co.,Ltd.

Address before: 313000 No. 666 Huanzhu Road, Wuxing District, Huzhou City, Zhejiang Province

Patentee before: NUCLEAR INDUSTRY HUZHOU ENGINEERING INVESTIGATION INSTITUTE

CP03 Change of name, title or address