CN112302657A - 一种露天矿卡车运输外排方法 - Google Patents

一种露天矿卡车运输外排方法 Download PDF

Info

Publication number
CN112302657A
CN112302657A CN202010986549.8A CN202010986549A CN112302657A CN 112302657 A CN112302657 A CN 112302657A CN 202010986549 A CN202010986549 A CN 202010986549A CN 112302657 A CN112302657 A CN 112302657A
Authority
CN
China
Prior art keywords
block
year
discharge
model
constraint
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
CN202010986549.8A
Other languages
English (en)
Other versions
CN112302657B (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.)
Ccteg Shenyang Engineering Co ltd
Original Assignee
Ccteg Shenyang Engineering Co ltd
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 Ccteg Shenyang Engineering Co ltd filed Critical Ccteg Shenyang Engineering Co ltd
Priority to CN202010986549.8A priority Critical patent/CN112302657B/zh
Publication of CN112302657A publication Critical patent/CN112302657A/zh
Application granted granted Critical
Publication of CN112302657B publication Critical patent/CN112302657B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21C—MINING OR QUARRYING
    • E21C47/00—Machines for obtaining or the removal of materials in open-pit mines
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种露天矿卡车运输外排方法，涉及露天矿排土工程领域。包括获取前期数据的步骤，建立外排土场块体模型的步骤，建立0‑1整数规划模型获取外排土场各排弃年份具体堆置形态的步骤，该方法考虑到包括运输费用和占地费用的排弃费用，可保证整个外排期内剥离物排弃总费用最低。同时，传统的卡车运输外排土场设计需要多次计算才能满足排土容量要求，并存在一定的误差，本发明可一次满足排土容量且不存在误差，其外排方法更加的准确和可靠。

Description

一种露天矿卡车运输外排方法
技术领域
本发明涉及露天矿排土工程领域，特别涉及一种露天矿卡车运输外排方法。
背景技术
露天矿排土工程前期无足够内排空间排弃剥离物，需要一定时间的外排期，排弃到规划的外排土场。目前，露天矿外排土场设计通常依据进度计划的内排容量和剥离总量确定外排年份及排弃量，排弃位置设计只是满足了容量要求，例如传统卡车运输外排土场设计确定排弃区域后计算排土体积，小于或大于排弃量，就需要重新确定排弃区域，反复人工调整后才能满足排弃空间要求，该方法费时费力，误差大且排弃费用高。
现有方法无法实现外排土场排弃时间及排弃位置最优，由于外排土场排弃位置只是通过手工设计，没有理论依据，依赖操作者的经验和熟练程度，导致该方法的使用受到限制。
发明内容
针对现有技术存在的不足，本发明的目的是提供一种露天矿卡车运输外排方法。
本发明所采用的技术方案是：一种露天矿卡车运输外排方法，其技术要点是，包括以下步骤：
获取前期数据的步骤，包括根据前期进度计划确定外排年份及考虑松散系数松方量的各年外排量、外排年份采坑出入沟位置、外排土场实体模型、运输单价和占地面积单价；
建立外排土场块体模型的步骤，利用获取的前期数据，建立块体模型，提取块体的质心点数据，由块体三维坐标和块体尺寸计算块体排弃费用为块体编号；
建立0-1整数规划模型，模型编写及求解，获取外排土场各排弃年份具体堆置形态。
上述方案中，所述的块体模型依据外排土场实体模型X、Y、Z三轴范围建立，用外排土场实体模型内部约束建立外排土场块体模型，块体模型平行X、Y轴建立，块体X轴、Y轴方向大小为卡车运输最终平盘宽度，Z轴方向大小为台阶高度。
上述方案中，所述计算块体排弃费用的过程为：把块体模型的各个质心点数据提取出来，根据块体三维坐标区分块体在三轴方向上的位置，定位块体所在台阶位置；同时根据块体三维坐标、块体尺寸及台阶位置计算排弃费用。
上述方案中，所述的计算排弃费用的过程为：根据块体所在台阶位置计算块体质心坐标到本台阶斜坡道水平距离，量取块体本台阶斜坡道到各外排年份采坑出入沟距离；根据运输单价计算块体从采坑出入沟到排弃位置的运输费用，最下块体还要根据块体尺寸和占地面积单价计算占地费用；以外排第1年为基础年，并计算折现后费用；根据运输费用与占地费用求和并折现后的费用确定为排弃费用。
上述方案中，所述的0-1规划模型，其目标函数为：
Figure BDA0002689443940000021
式中，I为台阶数；J为X轴方向块体数量；K为Y轴方向块体数量；T为时段数，a；Fi,j,k,t为第t时段编号(i,j,k)块体排弃费用；Xi,j,k,t为第t时段编号(i,j,k)块体决策变量，取值0或1。
0-1整数规划模型的目标函数为整个外排阶段排弃费用和最小，约束条件包括决策变量0-1约束、时空顺序约束、排弃量约束及排弃约束，决策变量0-1约束表示块体在哪个外排年份排弃，0表示不排弃，1表示排弃；时空顺序约束表示下部台阶周围块体本年份未排弃，上部台阶块体本年份及以前年份不能排弃；排弃量约束表示每年外排块体体积和等于本年外排量；排弃约束表示块体只能排弃一次或不排弃。
本发明的有益效果是：本发明的露天矿卡车运输外排方法，包括获取前期数据的步骤，建立外排土场块体模型的步骤，建立0-1整数规划模型，获取外排土场各排弃年份具体堆置形态的步骤，该方法考虑到包括运输费用和占地费用的排弃费用，可保证整个外排期内剥离物排弃总费用最低。同时，传统的卡车运输外排土场设计需要多次计算才能满足排土容量要求，并存在一定的误差，本发明可一次满足排土容量且不存在误差，其外排方法更加的准确和可靠。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可根据这些附图获得其他的附图。
图1为本发明实施例中一种露天矿卡车运输外排方法流程图。
具体实施方式
使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图1具体实施方式对本发明作进一步详细的说明。
本发明的露天矿卡车运输外排规划方法以某实例为例，具体步骤如下：
步骤1：前期数据准备。包括前期进度计划确定的外排年份及考虑松散系数松方量的各年外排量，用于确定时段数和各时段排弃量，剥离物开挖后比原始状态体积会增大，松散系数为1.1；外排年份采坑出入沟位置，用于确定各时段运距计算起点；外排土场实体模型，用于约束块体模型；运输单价和占地面积单价，用于与块体运距和计算排弃占地面积费用；外排年份及各年外排量如表1所示。
表1外排年份及各年外排量
Figure BDA0002689443940000031
运输单价为1km内基础单价4.73元/m3，每增加1km，增加1.67元/m3。占地面积单价为87.4元/m2。
步骤2：外排土场块体模型建立。在3Dmine软件中依据外排土场实体模型X、Y、Z三轴范围建立块体模型，用外排土场实体模型内部约束建立外排土场块体模型。为方便后期块体编号，块体模型平行X、Y轴建立；为方便后期规划模型条件约束编写，块体X、Y轴方向大小为卡车运输最终平盘宽度，设为40，Z轴方向大小为台阶高度，设为15。
步骤3：块体质心点数据提取。把块体模型的各个质心点数据提取出来，共计6074个点，部分块体质心点三维坐标如表2所示。
表2块体质心点三维坐标
Figure BDA0002689443940000032
根据块体三维坐标区分块体在三轴方向上的位置，定位块体所在台阶位置；同时根据块体三维坐标和块体尺寸进行费用计算。
步骤4：块体排弃费用计算。根据块体所在台阶位置计算块体质心坐标到本台阶斜坡道水平距离，量取块体本台阶斜坡道到各外排年份采坑出入沟距离。根据运输单价计算块体从采坑出入沟到排弃位置的运输费用，最下块体还要根据块体尺寸和占地面积单价计算占地费用。以外排第1年为基础年，考虑折现率，可计算折现后费用，折现率为10％。某块体各外排年份排弃费用见表3。
表3块体各外排年份排弃费用
Figure BDA0002689443940000033
Figure BDA0002689443940000041
步骤5：块体编号。根据块体质心点空间位置进行编号，为方便后期规划模型条件约束编写，虚拟一些块体，与实际排弃块体形成长方体。X、Y、Z轴形成60×66×3块体模型。Z轴方向根据质心点Z值从小到大编号，1代表最下台阶。X轴方向根据质心点X值从大到小编号，1代表XY平面最右侧块体。Y轴方向根据质心点Y值从小到大编号，1代表XY平面最下侧块体。块体质心点三维坐标及编号范围见表4。
表4块体质心点三维坐标及编号范围
Figure BDA0002689443940000042
步骤6：0-1整数规划模型建立。目标函数为整个外排阶段排弃费用和最小。本实施例中的0-1整数规划模型目标函数如下：
Figure BDA0002689443940000043
式中，I为台阶数(Z轴方向块体数量)；J为X轴方向块体数量；K为Y轴方向块体数量；T为时段数，a；Fi,j,k,t为第t时段编号(i,j,k)块体排弃费用；Xi,j,k,t为第t时段编号(i,j,k)块体决策变量，取值0或1。目标函数求解值为整个外排阶段排弃费用和最小。
(1)时空顺序约束，公式如下：
1)此约束限制正下方临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000044
式中：τ表示时段，为实现从第1时段到时段τ的求和，又区别于t，所以加入τ，下式相同。
2)此约束限制下一台阶右侧临近块体排弃完才能排弃本块体。左右侧表示X轴方向。
Figure BDA0002689443940000045
3)此约束限制下一台阶左侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000046
4)此约束限制下一台阶下侧临近块体排弃完才能排弃本块体。上下侧表示Y轴方向。
Figure BDA0002689443940000051
5)此约束限制下一台阶上侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000052
6)此约束限制下一台阶右下侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000053
7此约束限制下一台阶左下侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000054
8)此约束限制下一台阶右上侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000055
9)此约束限制下一台阶左上侧临近块体排弃完才能排弃本块体。
Figure BDA0002689443940000056
(2)排弃量约束，公式如下：
此约束限制各时段排弃量，保证完成各时段排弃计划。
Figure BDA0002689443940000057
式中，Vi,j,k为编号(i,j,k)块体体积；Vt为t时段排弃量。
(3)排弃约束，公式如下：
此约束限制每个块体最多只能排弃一次或不排弃。
Figure BDA0002689443940000058
步骤7：模型求解。根据规划模型在Lingo软件中进行集合段、数据段和目标约束段的编写。集合段编写4个基本集合，其中台阶(Z轴)集合有3个元素，X轴集合有60个元素，Y轴集合有66个元素，时段集合有3个元素。由基本集合组合成派生集合，台阶集合、X轴集合、Y轴集合和时段集合生成3×60×66×3的4维派生集合(表示决策变量和块体排弃费用)；台阶集合、X轴集合和Y轴集合生成3×60×66的3维派生集合(表示块体体积，为减小运算基数，块体体积表示为0或1)。数据段赋值台阶数为3，X轴方向块体数量为60，Y轴方向块体数量为66。目标约束段分时段编写时空顺序约束、排弃量约束(为减小运算基数，排弃量约束转化为块体数量约束)、排弃约束，定义0-1决策变量，编写目标函数，赋值块体排弃费用及块体体积(虚拟块体赋值为0)。运行模型，得到最优解为108118.8万元和决策变量0或1的赋值。各时段排弃块体数量见表5。
表5各时段排弃块体数量
Figure BDA0002689443940000061
步骤8：模拟外排土场时空发展。根据决策变量0或1的赋值，赋值为1的块体即为本年份排弃块体，把质心点数据导入3Dmine软件就可以模拟外排土场各排弃年份具体堆置形态。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (5)

1.一种露天矿卡车运输外排方法，其特征在于，包括以下步骤：
获取前期数据的步骤，包括根据前期进度计划确定外排年份及考虑松散系数松方量的各年外排量、外排年份采坑出入沟位置、外排土场实体模型、运输单价和占地面积单价；
建立外排土场块体模型的步骤，利用获取的前期数据，建立块体模型，提取块体的质心点数据，由块体三维坐标和块体尺寸计算块体排弃费用，为块体编号；
建立0-1整数规划模型，模型编写及求解，获取外排土场各排弃年份具体堆置形态。
2.如权利要求1所述的露天矿卡车运输外排方法，其特征在于，所述的块体模型依据外排土场实体模型X、Y、Z三轴范围建立，用外排土场实体模型内部约束建立外排土场块体模型，块体模型平行X、Y轴建立，块体X轴、Y轴方向大小为卡车运输最终平盘宽度，Z轴方向大小为台阶高度。
3.如权利要求1所述的露天矿卡车运输外排方法，其特征在于,所述计算块体排弃费用的过程为：把块体模型的各个质心点数据提取出来，根据块体三维坐标区分块体在三轴方向上的位置，定位块体所在台阶位置；同时根据块体三维坐标、块体尺寸及台阶位置计算排弃费用。
4.如权利要求3所述的露天矿卡车运输外排方法，其特征在于，所述的计算排弃费用的过程为：根据块体所在台阶位置计算块体质心坐标到本台阶斜坡道水平距离，量取块体本台阶斜坡道到各外排年份采坑出入沟距离；根据运输单价计算块体从采坑出入沟到排弃位置的运输费用，最下块体还要根据块体尺寸和占地面积单价计算占地费用；以外排第1年为基础年，并计算折现后费用；根据运输费用与占地费用求和并折现后的费用确定为排弃费用。
5.如权利要求1所述的露天矿卡车运输外排方法，其特征在于，所述的0-1规划模型，其目标函数为：
Figure FDA0002689443930000011
式中，I为台阶数；J为X轴方向块体数量；K为Y轴方向块体数量；T为时段数，a；Fi,j,k,t为第t时段编号(i,j,k)块体排弃费用；Xi,j,k,t为第t时段编号(i,j,k)块体决策变量，取值0或1。
0-1整数规划模型的目标函数为整个外排阶段排弃费用和最小，约束条件包括决策变量0-1约束、时空顺序约束、排弃量约束及排弃约束，决策变量0-1约束表示块体在哪个外排年份排弃，0表示不排弃，1表示排弃；时空顺序约束表示下部台阶周围块体本年份未排弃，上部台阶块体本年份及以前年份不能排弃；排弃量约束表示每年外排块体体积和等于本年外排量；排弃约束表示块体只能排弃一次或不排弃。
CN202010986549.8A 2020-09-18 2020-09-18 一种露天矿卡车运输外排方法 Active CN112302657B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010986549.8A CN112302657B (zh) 2020-09-18 2020-09-18 一种露天矿卡车运输外排方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010986549.8A CN112302657B (zh) 2020-09-18 2020-09-18 一种露天矿卡车运输外排方法

Publications (2)

Publication Number Publication Date
CN112302657A true CN112302657A (zh) 2021-02-02
CN112302657B CN112302657B (zh) 2024-05-17

Family

ID=74483549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010986549.8A Active CN112302657B (zh) 2020-09-18 2020-09-18 一种露天矿卡车运输外排方法

Country Status (1)

Country Link
CN (1) CN112302657B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114266379A (zh) * 2021-11-22 2022-04-01 中煤科工集团沈阳设计研究院有限公司 一种露天矿卡车运输短期内排方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2294433C1 (ru) * 2005-07-14 2007-02-27 Акционерная компания "АЛРОСА" (Закрытое акционерное общество) Способ отработки двух близкорасположенных рудных тел отдельными карьерами
CN102799744A (zh) * 2012-08-01 2012-11-28 中钢集团马鞍山矿山研究院有限公司 一种露天开采矿山冰碛土的排土工艺技术
CN102966355A (zh) * 2012-12-06 2013-03-13 马钢集团设计研究院有限责任公司 一种创造露天矿内排土条件的方法
CN107687342A (zh) * 2017-08-21 2018-02-13 中国矿业大学 一种露天煤矿拉斗铲工作线分段交替式采区转向方法
CN109284848A (zh) * 2017-07-21 2019-01-29 黄甫 内排露天矿长远排土规划方法及计算系统
CN110905519A (zh) * 2019-11-04 2020-03-24 中煤科工集团沈阳设计研究院有限公司 一种露天矿排土方法
CN111027219A (zh) * 2019-12-12 2020-04-17 神华准格尔能源有限责任公司 相邻同向协同开采露天煤矿公共运输帮的高度确定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2294433C1 (ru) * 2005-07-14 2007-02-27 Акционерная компания "АЛРОСА" (Закрытое акционерное общество) Способ отработки двух близкорасположенных рудных тел отдельными карьерами
CN102799744A (zh) * 2012-08-01 2012-11-28 中钢集团马鞍山矿山研究院有限公司 一种露天开采矿山冰碛土的排土工艺技术
CN102966355A (zh) * 2012-12-06 2013-03-13 马钢集团设计研究院有限责任公司 一种创造露天矿内排土条件的方法
CN109284848A (zh) * 2017-07-21 2019-01-29 黄甫 内排露天矿长远排土规划方法及计算系统
CN107687342A (zh) * 2017-08-21 2018-02-13 中国矿业大学 一种露天煤矿拉斗铲工作线分段交替式采区转向方法
CN110905519A (zh) * 2019-11-04 2020-03-24 中煤科工集团沈阳设计研究院有限公司 一种露天矿排土方法
CN111027219A (zh) * 2019-12-12 2020-04-17 神华准格尔能源有限责任公司 相邻同向协同开采露天煤矿公共运输帮的高度确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘光, 张瑞新: "露天矿排土及土岩流向流量优化模型", 化工矿物与加工, no. 03, 15 June 1991 (1991-06-15), pages 169 - 171 *
孙效玉;邓鹏宏;赵明;: "多时段露天矿生产计划整体优化模型", 东北大学学报(自然科学版), no. 10, 15 October 2016 (2016-10-15), pages 1460 - 1464 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114266379A (zh) * 2021-11-22 2022-04-01 中煤科工集团沈阳设计研究院有限公司 一种露天矿卡车运输短期内排方法

Also Published As

Publication number Publication date
CN112302657B (zh) 2024-05-17

Similar Documents

Publication Publication Date Title
CN107704703B (zh) 一种土建bim工程量模型的构建方法
CN109522651A (zh) 一种基于静态场和有偏行走的人群疏散模拟方法
CN107194103A (zh) 一种基于3dsMax的核设施模型辐射场剂量仿真方法
CN104732026B (zh) 一种金矿山设计矿石堆场和排土场的方法
CN112302657A (zh) 一种露天矿卡车运输外排方法
CN110008606A (zh) 一种基于bim风管自动加工方法
CN107563549A (zh) 一种基于bim模型的逃生路径生成方法、装置及设备
CN109446696A (zh) 基于cmaq模型的快速大气环境容量测算方法、存储介质和终端
CN106961671A (zh) 采集室内定位数据的方法和装置
CN112784503B (zh) 一种闸坝调度下水鸟栖息地适宜性计算方法与装置
CN109829195A (zh) 基于bim的土方工程量计算方法
CN107944203A (zh) 一种风速流线可视化的建筑设计方法
CN111666683A (zh) 一种高炉炉内炉料分布状态预测方法及系统
CN109033558A (zh) 基于Revit平台的外脚手架模型建模方法
CN111536843B (zh) 工业电子雷管起爆网路延期时间设计以及爆破方法
CN114596414A (zh) 一种基于无人机测绘三维地貌土方量实时分析方法及系统
CN106651137A (zh) 一种工程数量编码和标准化方法及系统
CN110378023A (zh) 屋面降雨径流模拟方法、装置、计算机设备和存储介质
CN109163708A (zh) 无人机大比例尺地形图房屋线的优化方法、装置与介质
CN110300413B (zh) 一种基于密度算法的通信站点仿真规划方法与系统
CN112407184A (zh) 用于计算船舶重量重心分布的方法及系统
CN112347543A (zh) 一种基于bim的车辆基地场坪土石方平衡设计方法
CN114266379A (zh) 一种露天矿卡车运输短期内排方法
CN111125825B (zh) 古木建筑遗产智能化建模方法及装置
CN114579685A (zh) 一种基于ArcGIS的村镇避难设施的选址方法及其系统

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