一种金属露天矿5d时空动态排产计划模型构建方法
【专利摘要】本发明公开了一种金属露天矿5D时空动态排产计划模型构建方法，主要步骤包括：构建露天矿矿岩体三维块体排产模型；空间块体模型属性赋值；构建矿石价格和采矿及岩体剥离成本数据库；最终构建5D时空动态排产模型；本发明旨在露天矿境界内对矿岩体的开采进行动态模拟开采，利用该5D时空动态排产模型能够精确地模拟出露天矿的中长期开采进度及过程，并动态计算获得的资源净值，使矿山能够在多种中长期开采方案中优选资源净值最大的开采方案，还可进一步考虑贴现率，使得矿山企业优选出净现值NPV最大的中长期进度计划方案，在露天矿开拓优化设计领域有着广泛的应用前景。
【专利说明】
-种金属露天矿5D时空动态排产计划模型构建方法
技术领域
[0001] 本发明属于矿业系统工程及矿山优化技术领域，特别设及一种金属露天矿加时空 动态排产计划模型构建方法。
【背景技术】
[0002] 目前我国大多数金属露天矿山的中长期(通常超过5年W上称为中长期)生产计划 往往存在"一次编制，长期使用"的现象，且只重视某一静态时间点下中长期进度规划，在编 制计划期间主要考虑生产规模、矿石品位配矿、剥采比等采矿参数，但市场环境受多种随机 因素影响存在着诸多不确定性，矿山企业在计划实施过程中往往忽视了矿石市场价格和生 产成本等经济参数的变化，原有的生产计划编制方式往往无法快速反应市场的变化，无法 实现生产计划执行结果的预知、预演。多数金属露天矿山企业没有摆脱W往运种重规划、计 划，不重技术经济动态分析及投资收益动态分析的现象，一方面是由于技术上存在的不足， 另一方面是缺乏此方面的战略规划理念。矿山企业的管理者往往更习惯于W矿山的矿石总 开采量、合理剥采比等为年度生产的考核目标，而不是W投资收益最佳等经济指标进行生 产计划组织与实施，运使得实际的采矿生产背离了矿山企业W市场为导向，追求合理的净 现值为目标的方向，运种静态的中长期生产计划显然缺乏灵活性、实时性和科学性。因此， 为了将采矿生产过程中矿石价格等市场经济参数的变化和时间都反应到金属露天矿中长 期生产计划的动态编制中，本发明提出了一种将时间和矿石产品价格及矿体开采成本集成 到5D时空动态排产模型的构建方法。
【发明内容】
[0003] 为了克服上述现有技术的缺点，本发明的目的在于提供一种金属露天矿加时空动 态排产计划模型构建方法，基于空间块体模型，旨在露天采剥进度计划编制过程中，既为管 理者提供矿石开采量、剥采比、矿石品位等决策信息，又能为管理者提供在不同的开采顺序 下及开采时间安排下矿石价格及成本的变化，形成一种随着进度计划的动态编制，时间、空 间及经济参数动态变化的排产模型，从而为管理者提供一种寻求更好的进度计划或收益最 大的进度计划方案的方法。
[0004] 为了实现上述目的，本发明采用的技术方案是：
[0005] -种金属露天矿5D时空动态排产计划模型构建方法，包括如下步骤：
[0006] (1)在露天矿采剥境界内构建=维块体排产模型
[0007] 在露天开采境界内，将矿体和岩体划分成=维离散的块体，划分后矿岩体形成= 维空间块体模型，每一个块体即采剥时的最小开采单元;利用此块体模型，按照台阶的空间 顺序及时间来实现模拟开采，即形成空间块体排产模型；
[000引（2)块体模型属性赋值
[0009]利用地质勘探资料中的地质钻孔数据对每一个块体的金属品位进行估值，利用直 接赋值法来对每一个块体的比重、矿岩类型进行赋值，结合块体的体积，最终获取每一块体 中所含金属质量；
[0010] (3)构建矿石预测价格和矿体开采及岩体剥离成本数据库
[0011] 根据矿石的历史价格来对矿石在未来中长期的价格进行预测，并对开采成本和剥 离成本进行预测，构建价格和成本数据库；
[0012] (4)构建抓时空动态排产模型
[0013] 在所述=维块体排产模型基础上，加入开采时间为第四维，该开采时间下的预测 价格和成本为第五维，即为空间开采块体进行开采时间和价格、成本赋值，从而形成5D时空 动态排产模型。
[0014] 所述步骤(1)中，=维空间块体模型的块体的大小根据实际需要来进行划分，块体 越大误差越大，块体越小矿体及岩体拟合的越好，误差越小。
[0015] 所述步骤(1)中，在露天矿境界与矿体的交叉处，采用比露天矿境界内的块体尺寸 小的次级块体进行划分。
[0016] 所述步骤(2)中
[0017] 若第i开采块体为矿石，则其质量Hbi=PoVi;
[001引矿石所含的金属量为Mi=m。巧i;
[0019] 若第i开采块体为岩石，则其质量为mwi = PwVi;
[0020] 式中i为开采块体的编号，取值为1，…，n，化为矿石比重，Pw为岩石比重，Vi为第i块 可采块体的体积，gi为直接利用克里格法或距离幕次反比法来对第i块块体的品位的赋值。 [0021 ]所述步骤(3)中，矿石价格数据库建立方法为：
[0022] 收集近十年的矿石价格数据，并对矿石价格的变化趋势进行分析预测，通过分析 前十年的价格变化，分别形成矿石的未来的乐观价、最可能价及悲观价，将乐观价作为最高 价，悲观价作为最低价，预测第j年矿石的期望价格門，并最终计算出每单位金属质量的最 终价格Qj:
[0023]
[0024] Qj = IOOXPjXp
[0025] 其中，T表示露天矿开采总年限，P表示精矿价格的金属含量百分比。
[0026] 所述步骤(3)中，矿体开采及岩体剥离成本数据库建立方法为：
[0027] 采集近十年的采矿及岩体剥离成本数据，计算出矿体开采成本平均增长率ri和岩 体剥离成本平均上涨率n，然后利用上一年的矿体开采成本OCw和岩体剥离成本WCw,计 算下一年度的矿体开采成本和岩体剥离成本，即
[002引 OCj = OCj-I X (1+ri)
[0029] WCj=WCj-I X (l+n)
[0030] 所述步骤(4)中，采用人机交互的方式，按照台阶的开采顺序进行模拟开采，在模 拟开采的过程中加入开采时间和该开采时间下的价格和成本形成5D时空排产模型。
[0031] 根据所述加时空动态排产模型，计算第j年开采第i块开采单元的资源净值Wi,若 开采单元为矿石，贝U
[0032] W〇i = QjXMiXSX(l-f)(l-k)-m〇iXOCj(i = l，...，m)
[0033] 若开采单元为岩石，贝U
[0034]
[0035]
[0036]
[0037]
[00；3 引
[0039] 值为
[0040]
[0041] 上述公式中S为矿石选矿回收率，f为矿石损失率，k为矿石贫化率，m为开采块体为 矿石的总块数，n-m为开采块体为岩石的总块数。与现有技术相比，本发明通过在金属露天 矿=维空间排产模型基础上引入时间（第四维）和矿石价格及成本(第五维）的方法来构建 露天矿中长期进度计划排产模型，能够将原有的不关屯、矿石价格和成本进度计划静态编制 方法改变，形成成为一种随着时间动态变化的露天矿中长期5D时空排产模型，为金属露天 矿的中长期进度计划编制提供了更为精细化的动态编制方法，在露天矿开拓优化设计领域 有着广泛的应用前景。
【附图说明】
[0042] 图1是本发明=维块体模型构建示意图，图中数字为块体划分编号。
[0043] 图2是本发明某矿岩体划分后形成排产模型示意图，图中数字代表放大100倍后得 矿石品位数据。
[0044] 图3是本发明露天矿境界内时空动态排产模型示意图。
【具体实施方式】
[0045] 下面结合附图和实施例详细说明本发明的实施方式。
[0046] 如图1所示，一种金属露天矿5D时空动态排产模型的构建方法，具体的实施步骤如 下：
[0047] 步骤1:在露天开采境界内，将矿体和岩体划分成=维离散的开采块体，块体的大 小根据实际需要来进行划分，通常要依据台阶高度、露天境界大小、模拟开采的精度等要求 来划分，如开采台阶大多数金属露天矿通常为12米，因此境界内块体大小可W采用长10米、 宽10米、高6米的块体来划分，划分后每一个块体即为境界内的最小开采单元;在露天矿境 界与矿体的交叉处，为了提高矿体开采模拟精度，可采用较小的次级块体来划分，如采用长 5米、宽5米、高3米的次级块进行露天矿边界处矿体的划分，运样划分后矿体形成包含n块块 体的=维空间块体排产模型，利用此模型可按照台阶的空间顺序及时间来实现模拟开采， 即形成=维空间排产模型。
[0048] 步骤2:在=维空间块体模型构建之后，利用矿区地质勘探的钻孔品位数据对每一 个块体进行品位估值，可W直接利用克里格法或距离幕次反比法来对第i块块体的品位进 行赋值为gl，根据gl可W将开采块体划分为矿石和岩石两种类型，然后根据=维排产模型中 矿体的空间大小可W求出第i开采块体的体积VI，根据矿岩体的类型，可得到其比重P，如矿 石Pd通常为3.2吨/立方米，岩体Pw通常为2.9吨/立方米，依据如下公式计算：
[0049]若第i开采块体为矿石，则其质量m〇i =化Vi;
[(K)加]矿石所含的金属量为Mi=m。巧i;
[0051] 若第i开采块体为岩石，则其质量为mwi = PwVi。
[0052] 式中i为开采块体的编号，取值为1，…，n。
[0053] 步骤3:收集十年W上的该精矿的市场价格，根据该价格在近十年的变化趋势，分 别形成精矿的未来的乐观价、最可能价及悲观价，利用如下公式计算预测每一年精矿的期 望价格門，如2015年钢精矿平均价格为800元/吨度，根据露天矿初步设计中的开采年限分 别求出下一年度的精矿期望价格，根据精矿价格的金属含量百分比P，可W计算出每单位
(吨）金属席吾的愚给价松化 [0化4]
[0化5]
[0056] 巧9萊4: H义果下牛W王巧失W々的化米W伞W君々君体剥离成本数据，成本数据一 般呈现逐年上涨趋势，根据十年的相关情况算出矿体开采成本平均上涨率ri和岩体剥离成 本平均上涨率K，然后利用上一年的矿体开采成本OCw和岩石岩体剥离成本WCw可W计算 下一年度的矿体开采成本和岩体剥离成本，通常单位为元/吨，即 [0化7] OCj = OCj-I X (1+ri)
[005引 WCj=WCj-I X (l+n)
[0059] 步骤5:依据上述步骤构建的=维空间排产模型和矿石价格及成本数据库，采用人 机交互的方式来按照台阶的开采顺序进行模拟开采，在模拟开采的过程中就依据开采时间 为=维排产模型，加入开采时间（第四维)和该开采时间下的价格和成本(第五维），即进行 时间和价格与成本赋值，形成加时空排产模型，然后可W计算出第j年开采第i块开采单元 的资源净值Wi，若开采的块体单元为矿石，贝U
[0060] Woi = Qj XMiXSX (l-f)(l-k)-m〇iXOCjQ =
[0061] 若开采的块体单元为岩石，即为需要剥离的废石，贝U
[0062] Wwi = -fflwi X wCj (i = I, ,n-m)
[0063] 其中，i为块体数量编号，m为开采单元为矿石的块体总数，n为矿石和岩石块体总 数;j为块体开采年度，T为露天矿开采总年限，j = l，…，T;S为矿石选矿平均回收率，f为矿 石损失率，k为矿石贫化率。
[0064] 根据贴现率可W将其进一步转换为净现值NPV，本发明暂且考虑资源净值Wi。
[0065] 步骤6:根据台阶模拟开采的进行可W依次计算出矿石的开采矿量、品位、体积、岩 石的剥离量、剥采比等，本发明与=维排产模型不同的是在编制的过程中随着块体开采的 进行，根据该块体在排产计划中的开采时间W及在该开采时间约束下的矿石价格和开采成 本，运样就可W随着进度计划编制的过程动态的计算出该进度计划方案下的主要参数，计 算如下：
[0066] 总开采矿i
，111为开采块体为矿石的总块数；
[0067] 总剥离岩石量 n-m为开采块体为岩石的总块数；
[0068] 总平均品位关
[0069] 该方案平均熟
[0070] 该进度计划排产获得的总的资源净值为
[0071]
[0072] 步骤7:利用该时空排产模型，每进行一次进度计划的编制，都会产生该开采计划 编制下的企业所获得的资源净值，若总共有k套不同的进度计划编制方案，则找到资源总净 值最大Wmax,对应的该开采方案即为不考虑贴现率情况下的最优进度计划编制方案。
[007；3] Wmax=Max{Wi，…，Wkl
[0074]本发明将矿体和岩体利用=维空间块体进行划分形成最小开采单元来模拟开采 过程，然后依据开采年限构建了矿石市场价格预测数据和矿石开采成本及岩体剥离成本数 据，在模拟开采的过程中分别将开采时间和该开采时间下的价格和成本数据增加到=维排 场模型，形成5D时空排产模型。上面结合附图对本发明进行了示例性描述，显然本发明具体 实现并不受上述方式的限制，只要采用了本发明的方法构思和技术方案进行的各种非实质 性的改进，或未经改进将本发明的构思和技术方案直接应用于其它场合的，均在本发明的 保护范围之内。
【主权项】
1. 一种金属露天矿时空动态排产计划模型构建方法，其特征在于，包括如下步骤： (1) 在露天矿采剥境界内构建三维块体排产模型 在露天开米境界内，将矿体和岩体划分成三维离散的块体，划分后矿岩体形成三维空 间块体模型，每一个块体即采剥时的最小开采单元;利用此块体模型，按照台阶的空间顺序 及时间来实现模拟开采，即形成空间块体排产模型； (2) 块体模型属性赋值 利用地质勘探资料中的地质钻孔数据对每一个块体的金属品位进行估值，利用直接赋 值法来对每一个块体的比重、矿岩类型进行赋值，结合块体的体积，最终获取每一块体中所 含金属质量； (3) 构建矿石预测价格和矿体开采及岩体剥离成本数据库 根据矿石的历史价格来对矿石在未来中长期的价格进行预测，并对开采成本和剥离成 本进行预测，构建价格和成本数据库； (4) 构建?时空动态排产模型 在所述三维块体排产模型基础上，加入开采时间为第四维，该开采时间下的预测价格 和成本为第五维，即为空间开采块体进行开采时间和价格、成本赋值，从而形成5D时空动态 排产模型。2. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，所 述步骤（1)中，三维空间块体模型的块体的大小根据实际需要来进行划分，块体越大误差越 大，块体越小矿体及岩体拟合的越好，误差越小。3. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，所 述步骤（1)中，在露天矿境界与矿体的交叉处，采用比露天矿境界内的块体尺寸小的次级块 体进行划分。4. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，所 述步骤(2)中 若第i开采块体为矿石，则其质量; 矿石所含的金属量为Mi=nkigi; 若第i开采块体为岩石，则其质量为; 式中i为开采块体的编号，取值为1，···，n，P。为矿石比重，pw为岩石比重，Vi为第i块可采 块体的体积，gl为直接利用克里格法或距离幂次反比法来对第i块块体的品位的赋值。5. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，所 述步骤(3)中，矿石价格数据库建立方法为： 收集近十年的矿石价格数据，并对矿石价格的变化趋势进行分析预测，通过分析前十 年的价格变化，分别形成矿石的未来的乐观价、最可能价及悲观价，将乐观价作为最高价， 悲观价作为最低价，预测第j年矿石的期望价格匕，并最终计算出每单位金属质量的最终价 格Qj: Qj = 100 XPj 八其中，T表;^露天矿开米总年限，p表;^精矿价格的金属含量百分比。6. 根据权利要求1或5所述金属露天矿5D时空动态排产计划模型构建方法，其特征在 于，所述步骤(3)中，矿体开采及岩体剥离成本数据库建立方法为： 采集近十年的采矿及岩体剥离成本数据，计算出矿体开采成本平均增长率rdP岩体剥 离成本平均上涨率r2,然后利用上一年的矿体开采成本OCh和岩体剥离成本WCn，计算下 一年度的矿体开米成本和岩体剥离成本，BP OCj = OCj-iX (1+n) WCj=ffCj-iX (l+r2)7. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，所 述步骤(4)中，采用人机交互的方式，按照台阶的开采顺序进行模拟开采，在模拟开采的过 程中加入开采时间和该开采时间下的价格和成本形成时空排产模型。8. 根据权利要求1所述金属露天矿f5D时空动态排产计划模型构建方法，其特征在于，根 据所述5D时空动态排产模型，计算第j年开采第i块开采单元的资源净值1，若开采单元为 矿石，则 W〇i = QjXMiXSX (l-f)(l-k)-m〇iXOCj(i = l,---,111) 若开采单元为岩石，则 ffwi = ~mwi X wCj (i = 1, · · ·, n-m) 总开采矿量拟0 = Σ[^ι爪; 总剥离岩石量二 总平均品位为《g二 该方案平均剥采比为δ = f ιηο 该进度计划排产获得的总的资源净值为上述公式中S为矿石选矿回收率，f为矿石损失率，k为矿石贫化率，m为开采块体为矿石 的总块数，n-m为开采块体为岩石的总块数。
【文档编号】G06F17/50GK105956928SQ201610300267
【公开日】2016年9月21日
【申请日】2016年5月6日
【发明人】顾清华, 卢才武, 张雯, 聂兴信, 杨震, 郭梨
【申请人】西安建筑科技大学