本发明涉及采矿
技术领域：
，特别地，涉及一种堆积型铝土矿多矿区的开采规划方法。
背景技术：
：堆积型铝土矿分布于峰丛注地、峰林谷地及斜坡地带。矿体平面形态复杂、矿体体积小、分布点多。矿山适用露天开采，公路运输开拓，开采工艺有装载机液压反铲、推土机液压反铲和推土机装载机开采工艺。矿山开采计划是指导矿山生产的依据，计划编制是矿山生产管理中不可缺少的重要部分，如何确定各生产时期的开采对象及其开采量，以及如何实现矿山采场保有资源的更新，将关系到整个矿山的开采效果。目前国内对运用计算机编制采掘生产计划的研究归纳起来主要是两大类：一是运用优化理论，主要是单目标线性规划理论，在露天矿有极少数运用了动态规划理论对采掘生产计划予以优化，并借助于计算机实现优化计算；二是直接运用计算机模拟人工编制的方法，或者逼近的方法进行计划的编制。此两者方法均有目标函数难确定、变量模糊以及求解迭代计算过程复杂等问题。因此，设计一种能够解决上述问题的堆积型铝土矿多矿区的开采规划方法具有重要意义。技术实现要素：本发明的目的在于提供一种堆积型铝土矿多矿区的开采规划方法，本发明按生产时所需考虑的约束条件和指定采场采完优先顺序，解决目标函数难、变量模糊的问题；为解决规划无解的问题，引入惩罚函数，并就目标函数进行变量的分类，降低了整个规划求解迭代计算的复杂程度。具体技术方案如下：一种堆积型铝土矿多矿区的开采规划方法，包括以下步骤：第一步：根据计划目标值，以卸矿平台为目标点，以总运距最小作为目标按运距的远近进行采场选取，最小总运距采用表达式1)计算：s为规划期内总运距，mins为规划期内最小总运距，n为采场个数，sj为第j采场的运距，单位为千米，xj为第j采场的本期实际开采量，单位为吨；第二步：判断合格矿量是否符合要求，具体是：若表达式2)成立，则合格矿量符合要求，否则返回第一步；其中：n为目标总矿量；j为采场编号，取值为1到n的自然数；xj为第j采场目标开采矿量，单位为吨；kj为第j采场可开采矿量；第三步：判断三氧化二铝平均品位是否满足要求，具体是：若表达式3)成立，则满足要求，否则返回第一步；其中：αj为第j采场中三氧化二铝的含量；a1为第一设定值；第四步：判断铝硅比是否满足要求，具体是：若表达式4)成立，则满足要求，否则返回第一步；其中：βj为第j采场中的铝硅比；a2为第二设定值，a3为第三设定值；第五步：判断指定完成采场约束是否符合要求，具体是：若表达式5)成立，则满足要求，否则返回第一步；其中：为本期末第i采场的可采矿量；第六步：通过改进的堆积型铝土矿动态开采计划优化模型得到各采场本期实际开采量，制定开采计划表；执行采矿；改进的堆积型铝土矿动态开采计划模型详见表达式6)：其中：ki为第i采场本期必须参与计划的可采矿量；i为采场编号，取值为1到n的自然数；si为第i采场的运距；σ为罚因子；α(xj)为惩罚函数；f(xj)为加入惩罚函数后的堆积型铝土矿动态开采计划模型；第七步：返回第一步，进入下一轮采矿计划表的制定。以上技术方案中优选的，改进的堆积型铝土矿动态开采计划模型求解过程如下：目标函数求解转换为用于解决约束条件下的最优化问题，通过惩罚函数可以将有约束的目标函数转化为无约束的目标函数；指定完成采场约束转换为目标函数的变量减少约束条件；作为目标函数分解出的变量，可以通过生产经验计算出结果值，将结果值直接参与规划计算，实现了减少目标函数的变量数；根据变量引入σα(xj)解决规划无解的问题，具体是：最优解成立时，xj的取值为规划结果值用于制作开采计划表，其中：minf(x，σ)为最优解；惩罚函数引入求解过程包括矿石量约束和矿山质量约束，矿石量约束详见表达式8)，矿石质量约束详见表达式9)：以上技术方案中优选的，计划目标值包括矿量、品位和含矿率；采场的选取具体是：在一条主干公路上由近及远的顺序按照年平均运距逐年增加的原则进行选取；采场包括系统自动选取的采场和/或生产管理人员指定完成的采场；采场选取重新进行时，能用一个矿量大的采场替换一个矿量小的采场和/或增加新的采场。以上技术方案中优选的，所述开采计划表分期制定分期实施，具体为一周至两周更新一次。应用本发明的方法，具体效果是：1、本发明按生产时所需考虑的约束条件和指定采场采完优先顺序，迭代前分解指定采场采完的计算变量，且分解指定采场采完的约束计算结果参与每一次迭代计算，确保最优解的情况与实际相符，即解决目标函数难、变量模糊的问题，在实际应用过程中方便采用计算机模拟方法调用模型，自动进行规划计算。2、本发明以运距为目标函数，以供矿要求、采场可采矿量、指定开采采场与指定完成采场等为约束条件，采用改进的开采计划模型，同时以生产量统计动态更新开采对象的储量信息进行计划的修正与反馈，实现生产与计划闭环系统。3、为解决规划无解的问题，基于带有约束条件的极值问题求解，改进开采计划编排数学模型(具体是引入惩罚函数)，并就目标函数进行变量的分类，降低了整个规划求解迭代计算的复杂程度。除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照附图，对本发明作进一步详细的说明。附图说明构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：图1是实施例1中堆积型铝土矿多矿区的开采规划方法的逻辑图。具体实施方式以下结合附图对本发明的实施例进行详细说明，但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。实施例1：一种堆积型铝土矿多矿区的开采规划方法，矿山开采计划因素包括：运距、采场可采矿量、供矿要求以及其它特定生产要求等，具体逻辑图详见图1。本实施例的方法具体包括以下步骤：第一步：根据计划目标值(包括矿量、品位和含矿率)，以卸矿平台为目标点，以总运距最小作为目标按运距的远近进行采场选取(在一条主干公路上由近至远的顺序按照年平均运距逐年增加的原则进行选取；采场包括系统自动选取的采场和/或生产管理人员指定完成的采场；采场选取重新进行时，能用一个矿量大的采场替换一个矿量小的采场和/或增加新的采场)，最小总运距采用表达式1)计算：s为规划期内总运距，mins为规划期内最小总运距，n为采场个数，sj为第j采场的运距，单位为千米，xj为第j采场的本期实际开采量，单位为吨；第二步：判断合格矿量是否符合要求，具体是：若表达式2)成立，则合格矿量符合要求，否则返回第一步；其中：n为目标总矿量；j为采场编号，取值为1到n的自然数；xj为第j采场目标开采矿量，单位为吨；kj为第j采场可开采矿量；第三步：判断三氧化二铝(al2o3)平均品位是否满足要求，具体是：若表达式3)成立，则满足要求，否则返回第一步；其中：αj为第j采场中三氧化二铝的含量；a1为第一设定值，可根据企业生产水平确定，如取a1＝53％；第四步：判断铝硅比是否满足要求，具体是：若表达式4)成立，则满足要求，否则返回第一步；其中：βj为第j采场中的铝硅比(al/si)；a2为第二设定值，a3为第三设定值，a2和a3的值均可根据企业生产水平确定，如取a2＝8.5，a3＝9.5；第五步：判断指定完成采场约束是否符合要求，具体是：若表达式5)成立，则满足要求，否则返回第一步；其中：为本期末第i采场的可采矿量；在此步骤之后，一般还有含矿率的判断，含矿率根据实际开采需求确定，若符合则进入下一步，若不符合则返回第一步；第六步：通过改进的堆积型铝土矿动态开采计划优化模型得到各采场本期实际开采量，制定开采计划表，开采计划表分期制定分期实施；执行采矿；改进的堆积型铝土矿动态开采计划模型详见表达式6)：其中：ki为第i采场本期必须参与计划的可采矿量；i为采场编号，取值为1到n的自然数；si为第i采场的运距；σ为罚因子；α(xj)为惩罚函数；f(xj)(构造函数)为加入惩罚函数后的堆积型铝土矿动态开采计划模型；改进的堆积型铝土矿动态开采计划模型求解过程如下：目标函数求解转换为用于解决约束条件下的最优化问题，通过惩罚函数可以将有约束的目标函数转化为无约束的目标函数；指定完成采场约束转换为目标函数的变量减少约束条件；作为目标函数分解出的变量，可以通过生产经验计算出结果值，将结果值直接参与规划计算，实现了减少目标函数的变量数；根据变量引入σα(xj)解决规划无解的问题，具体是：该策略为：对于在无约束的求解过程中企图违反约束的迭代点给予很大的目标函数值，迫使无约束问题的极小点或无限地向可行域d靠近，或保持在可行域d内，直到收敛到原来约束最优化问题的极小点。其中：α(xj)连续；取一个递增且趋于+∞的罚因子序列{σk}：0＜σ1＜σ2＜…＜σk＜σk+1＜…,σk→+∞惩罚函数α(xj)的构造：其中：当α(xj)→0，停止计算，得到近似最小值，表达式7)为最优解；最优解成立时，xj的取值为规划结果值用于制作开采计划表；惩罚函数引入求解过程包括矿石量约束和矿山质量约束，矿石量约束详见表达式8)，矿山质量约束详见表达式9)：第七步：返回第一步，进入下一轮采矿计划表的制定。迭代次数与参与的规划的采场数量有关，以50个采场参与规划为例，指定完成采场约束转换为目标函数的变量减少n个采场，可以减少(50-n)n迭代次数。此处优选给定停止迭代条件如下：第一步迭代停止条件：第二步迭代停止条件：第三步迭代停止条件：应用本实施例的技术方案，效果是：1、本发明按生产时所需考虑的约束条件和指定采场采完优先顺序，迭代前分解指定采场采完的计算变量，且分解指定采场采完的约束计算结果参与每一次迭代计算，确保最优解的情况与实际相符，即解决目标函数难、变量模糊的问题，在实际应用过程中方便采用计算机模拟方法调用模型，自动进行规划计算。2、本发明以运距为目标函数，以供矿要求、采场可采矿量、指定开采采场与指定完成采场等为约束条件，采用改进的开采计划模型，同时以生产量统计动态更新开采对象的储量信息进行计划的修正与反馈，实现生产与计划闭环系统。3、为解决规划无解的问题，基于带有约束条件的极值问题求解，改进开采计划编排数学模型(具体是引入惩罚函数)，并就目标函数进行变量的分类，降低了整个规划求解迭代计算的复杂程度。将停止迭代条件与约束条件进行合理规划，实现求解变量减少的目标，简化求解的过程。实施例2：通过实施例1的方法，进行数据测试应用，应用的结果如下：计算前可供计划的采场数据，如表1：表1采场数据统计表采场编号合格矿量al2o3(％)sio2(％)al/si26-2-8014098.9263.974.8313.2426-2-81129962.656.595.1311.0326-2-8288189.0958.144.9611.7226-2-8310806961.464.2614.4326-2-84206880.960.223.7116.2326-2-8612104.4764.374.5014.3026-2-8760946.1851.617.846.5826-2-882624.9656.488.146.94设定目标值：月目标矿量500000吨，约束条件：(al2o3(％))a1＝59％，(al/si)a2＝12.5、a3＝13.5，其中编号26-2-84为指定开采完成的采场，通过本发明方法计算的结果如表2：表2计算结果计划开采采场编号计划开采合格矿量/吨al2o3(％)sio2(％)al/si运距(km)26-2-8112000056.595.1311.036.4526-2-8310806961.464.2614.436.8526-2-84206880.960.223.7116.237.0026-2-8760946.1851.617.846.587.4526-2-882624.9656.488.146.947.55最终得到实际开采合格矿量为498521.04吨(五个采场的实际开采合格矿量进行加和)，差额为1478.96吨，一般做法是将下一月目标矿量做相应提高来进行补偿，也可以采用别的方式进行补偿，具体根据实际情况确定。以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12