季节性剥离露天煤矿工作帮形态动态控制及推进强度优化.pdf

季节性剥离露天煤矿工作帮形态动态控制及推进强度优化刘光伟1，黄云龙1，曹博1，姚勇2（1.辽宁工程技术大学矿业学院,辽宁阜新123000；2.锡林郭勒盟蒙东矿业有限责任公司,内蒙古锡林浩特026000）摘要：季节性剥离露天煤矿由于剥离工程的不连续性使得工作帮超前推进距离巨大，为了减小超前剥离、节省剥离运距，从季节性剥离露天煤矿工作帮形态动态控制、工作帮极限推进强度优化确定进行研究。阐述了季节性剥离露天煤矿工作帮形态及作业程序，分析了工作帮形态及推进强度影响因素。根据季节性剥离露天煤矿工作帮结构形式和采剥作业方式，建立了工作帮形态优化工程模型，论述了工作帮循环推进距离与作业平盘宽度和采场露煤方式的关系；计算了组合台阶循环推进工程量和分层设备生产能力，给出了工作帮循环推进时间以及端帮运输系统可靠性计算方法。以最小化剥离物运输功为优化目标、剥采工程时空接续为约束条件建立了季节性剥离露天煤矿工作帮形态优化数学模型并引入序列二次规划算法（SQP）进行非线性规划求解。提出通过确定剥离台阶组合方式和作业平盘宽度实现季节性剥离露天煤矿工作帮形态动态控制的方法、剥离工作帮极限推进度优化确定方法。以胜利西二号露天煤矿为工程背景，将上述优化控制方法进行应用研究，通过建立剥离工作帮形态优化模型并将西二露天矿实际参数代入进行计算求解。结果表明：当组合台阶中台阶数目为 5、作业平盘宽度为 89m 时，工作帮推进满足剥采工程时空接续要求，且剥离物运输功最小。该形态下工作帮极限推进度为 320m/a，分析确定了台阶组合方式及循环推进距离变化时工作帮极限推进强度变化规律。关键词：露天煤矿；季节性剥离；组合台阶；工作帮；序列二次规划(SQP)；推进度中图分类号：TD824文献标志码：A文章编号：02532336（2023）10004510Dynamic control of working slope shape and optimization of advance speed inseasonally stripped open-pit coal mineLIUGuangwei1,HUANGYunlong1,CAOBo1,YAOYong2（1.Institute of Mining Technology,Liaoning Technical University,Fuxin 123000,China;2.Xilin Gol Mengdong Mining Co.,Ltd.,Xilinhot 026000,China）Abstract:Duetothediscontinuityofthestrippingprojectinseasonalstrippingopen-pitcoalmines,theuntimelyadvancedistanceoftheworkingslopeishuge.Inordertoreducetheadvancestrippingandsavethestrippingdistance,researchwascarriedoutondynamiccon-troloftheshapeoftheworkingslopeandoptimizationofthelimitadvancespeedoftheworkingslopeinseasonalstrippedopen-pitcoalmines.Thispaperexplainstheshapeandoperationprocedureoftheworkingslopeinseasonalstrippedopen-pitcoalminesandanalyzesfactorsthatinfluenceworkingslopeshapeandadvancespeed.Basedonstructuralformofworkingslopeandstrippingoperationmodeofseasonallystrippedopen-pitcoalmines,anengineeringmodelforshapeoptimizationofworkingslopeisestablished.Therelationshipbetweencircularadvancedistanceofworkingslopewithworkingbenchwidthandexposedcoalmodeofstopeisdiscussed.Benchgroupcycleadvanceengineeringquantityandproductioncapacityoflayeredequipment,advancetimeoftheworkingslopecycleandcalcula-tionmethodforreliabilityofsideslopetransportationsystemaregiven.Withoptimizationgoalofminimizingtransportworkofstripping收稿日期：20230427责任编辑：朱恩光DOI：10.13199/ki.cst.2023-0621基金项目：国家自然科学基金资助项目（51974144）；辽宁省揭榜挂帅资助项目（2021JH1/10400011）；辽宁工程技术大学学科创新团队资助项目（LNTU20TD-07）作者简介：刘光伟（1981），男，辽宁沈阳人，教授，博士生导师。E-mail：liu_通讯作者：黄云龙（1994），男，山西朔州人，博士研究生。E-mail：第51卷第10期煤炭科学技术Vol.51No.102023年10月CoalScienceandTechnologyOct.2023刘光伟，黄云龙，曹博，等.季节性剥离露天煤矿工作帮形态动态控制及推进强度优化J.煤炭科学技术，2023，51（10）：4554.LIUGuangwei，HUANGYunlong，CAOBo，et al.Dynamiccontrolofworkingslopeshapeandoptimizationofad-vancespeedinseasonallystrippedopen-pitcoalmineJ.CoalScienceandTechnology，2023，51（10）：4554.45materialandtime-spacecontinuityofstrippingengineeringasconstraintcondition,mathematicalmodelforshapeoptimizationofworkingslopeinseasonallystrippedopen-pitcoalmineisestablished.Sequentialquadraticprogrammingalgorithm(SQP)isintroducedtosolvenonlinearprogramming.Amethodtorealizedynamiccontrolofshapeofworkingslopeofseasonallystrippedopen-pitcoalminebyde-terminingcombinationmodeofstrippingbenchesandworkingbenchwidthandamethodtooptimizelimitadvancespeedofstrippedworkingslopeareproposed.TakingShengliWestNo.2Open-pitCoalMineasanengineeringbackground,above-mentionedoptimizationcontrolmethodisappliedforresearchandcalculationandsolutionarecarriedoutbyestablishingformoptimizationmodelofstrippingworkingslopeandsubstitutingactualparametersofWestNo.2Open-pitMine.Resultsshowthatwhenthenumberofbenchesinbenchgroupis5andworkingbenchwidthis89m,advanceofworkingslopemeetstherequirementsofspace-timecontinuityinstrippingengin-eeringandtransportworkofstrippedobjectsissmallest.Inthisshape,limitadvancespeedoftheworkingslopeis320m/a,andchangeoflimitadvancespeedoftheworkingslopeisanalyzedanddeterminedwhencombinationofbenchesandcyclicadvancedistanceischanged.Key words:open-pitcoalmine；seasonallystripped；benchgroup；workingslope；sequentialquadraticprogramming(SQP)；advancingspeed0引言煤炭是我国的主体能源和重要原料，长期以来为经济社会发展和国家能源安全稳定供应提供了有力保障1-2。近年来，我国大力发展新能源技术，但依然没有改变我国以煤电为主的电力结构3。露天煤炭产量一直保持着良好的增长势头，占全国煤炭总产量的比重稳步升高4。露天煤矿的降本增效开发是实现煤炭行业高质量发展的重要内容。我国幅员辽阔，经度、纬度、地面高程跨度很大。寒区面积占我国陆地面积的 43.5%5-6，主要分布于青藏高原、东北地区、内蒙古东北部及新疆部分地区。高寒地区冬季漫长而寒冷，为有效避开冬季恶劣的剥离施工条件，降低生产成本，高寒地区露天煤矿剥离工程多采用季节性作业。将剥离工程集中在气候条件较好的时期，形成超前剥离7-8，剥离停滞期间采煤作业仍继续进行。我国大部分露天煤矿的剥离工艺为“单斗卡车”间断工艺，广泛采用剥离工程外包生产模式9。为了减小超前剥离，节省剥离物内排运距，季节性剥离露天煤矿工作帮常采用组合台阶循环推进作业。工作帮形式、台阶组合方式、采场露煤方式等直接影响剥离工作帮形态。围绕露天矿工作帮形态及季节性剥离露天煤矿剥采工程优化，周志友10、孙明亮等11分别分析了组合台阶开采工艺特点，阐述了组合台阶开采的重要意义以及组合台阶在露天矿实际生产中的应用；刘宇等12对组合台阶在哈尔乌素露天煤矿应用的具体参数进行了优化确定；白润才等13根据季节性剥离露天煤矿剥离工程的不连续性，提出了在采场和内排土场之间设置临时排土桥的方式来优化露天煤矿内排开拓运输系统；张周爱等7分析了季节性剥离露天矿剥离工程外包的特点，以宝日希勒露天煤矿为工程背景，给出了采场露煤次数及储备煤量的优化方法，提出了合理的采场露煤方案；马忠辉等14利用露天煤矿边坡岩土强度的季节性差异，提高了冬季开采时端帮靠帮开采边坡角，提高了端帮煤炭资源回收量，增加了经济效益。陆翔15研究了寒区露天矿泥岩边坡冻融机理，结合时效边坡理论及露天开采原理，提出了近水平煤层露天矿端帮煤层季节性控制开采方法。笔者基于相关研究成果，围绕季节性剥离露天煤矿工作帮形态动态控制进行研究。分析季节性剥离露天煤矿生产作业特点，将采矿工程问题抽象为相关非线性优化数学模型并引入序列二次规划算法进行求解，提出季节性剥离露天煤矿工作帮形态动态控制方法和工作帮极限推进度优化确定方法。并以胜利西二号露天煤矿为背景进行工程应用研究，确定了最优剥离工作帮形态、工作帮极限推进度。1季节性剥离露天煤矿工作帮形态及影响因素1.1季节性剥离露天煤矿工作帮形态季节性剥离露天煤矿由于其剥离作业集中于气候条件较好的时段，工作帮需要进行超前剥离以揭露足够的可采煤量以及储备煤量。为了推迟剥离、调节生产剥采比，常采用组合台阶作业方式增大工作帮坡角，进而节省运距，降低生产成本。常见的台阶组合方式为：把若干个相邻的开采台阶人为的划分为 1 组，在每一组台阶中同一时间内只设一个工作平盘和只在一个台阶进行采剥作业，这一组台阶就称之为组合台阶，如图 1 所示。图 1中所表示组合台阶为 H1及 H2两组，每组均由 3 个开采台阶组成，每组中同时按 1、2、3 及 1、2、3顺序进行采剥，完成一个循环之后形成图中虚线所示状态，之后在新的位置上按同样的方式进行采剥工2023年第10期煤炭科学技术第51卷46作，一直循环到预定的开采位置为准16。其中组合台阶高度即台阶组合方式以及每次循环推进距离为影响经济效果和技术复杂程度的关键性因素，影响露天矿山生产规模、工作帮年推进度以及工作线长度等。推进方向H1H2123123图1组合台阶构成示意Fig.1Schematicdiagramofcompositionofbenchgroup根据组合台阶作业程序，当各组组合台阶均完成一次作业循环时，工作帮方可向前推进一定距离。因此，组合台阶内台阶数目越多，完成一次循环推进时间越长，剥离作业期内循环作业次数越少，工作帮坡角越大。反之，组合台阶内台阶数目越少，完成一次循环时间越短，剥离作业期内循环作业次数越多，工作帮坡角越小，剥离物端帮运距越大。1.2工作帮形态及推进度影响因素1）采场露煤方式。煤矿露天开采的实质是将煤层上部岩石覆盖层剥离掉，使煤层敞露于地表后进行开采的过程，因此煤层的可采程度决定于剥离工程的进行程度。对于全年剥离作业的露天煤矿，开采方式为“随采随露”，工作帮不存在超前剥离，只需留设少量的储备煤量以满足开拓降深工程及剥离准备工作需要。因此，采煤台阶推进速度略小于剥离工作帮推进速度，也即采煤速度略小于露煤速度；季节性剥离露天煤矿由于剥离作业时间短，需要在剥离作业期间揭露大量的可采煤量来保证生产的接续，所以采煤速度远小于露煤速度。季节性剥离露天矿进行露煤作业时，露煤时长为剥离工作帮推进时长。全年一次露煤时，剥离工作帮沿推进方向完成一次循环推进。需要留设的储备煤量巨大，造成煤层揭露面积过大，不仅增加了煤层自燃的风险，而且大幅提高了剥离物内排运距。全年多次露煤时，可以有效缩短内排运距，降低煤层自燃风险，但剥离工作帮在完成一次循环推进后设备调动频繁，减小了设备实际作业时长，影响了设备能力的发挥。2）作业平盘宽度。组合台阶作业平盘宽度由循环推进距离和非作业平盘宽度构成。平盘宽度过大，达不到提高帮坡角的目的，平盘过窄，不利于设备作业，造成设备效率下降。3）剥离台阶组合方式。剥离台阶组合方式即为一组台阶内台阶数目，台阶数目越多，组合台阶高度越大，工作帮越陡。组合台阶内台阶数目是实行组合台阶推进，提高工作帮坡角的重要依据。4）运输道路可靠性。实行内排作业的露天煤矿在开采过程中工作帮产生剥离物料一般通过端帮道路运输至内排土场，端帮道路的通过能力应与工作帮剥离量相适应。工作帮循环推进距离过大时，通过端帮运输道路的剥离量过大，会导致端帮运输道路行车密度超过规定最大行车密度值。5）年剥离作业时间。季节性剥离露天煤矿须在剥离作业期内完成相应的剥离任务，来保证足够的露煤量，满足产能稳定接续的要求。工作帮形态不不发生变化时，年剥离作业时间越长，工作帮推进能力越大。2组合台阶工作帮推进分析2.1工作帮推进工程模型为了优化季节性剥离露天煤矿剥离工作帮形态，根据组合台阶作业程序及工作帮形态影响因素分析结果。基于以下前提条件：工作帮采用组合台阶推进，各组组合台阶内分层参数一致；采场露煤方式为全年多次露煤；工作帮形式为纵向工作帮。建立季节性剥离露天矿复合煤层工作帮推进工程模型如图 2 所示。假设可采煤层 2 层，分别为 A 煤和 B 煤，煤层近水平赋存。组合台阶作业平盘内划分若干剥离作业区，各作业区布置数组设备，共同完成工作面的剥离任务，各煤层台阶推进速度相同。端帮也为组合台阶，由运输平盘和安全平盘组成，每组台阶内只留设1 个运输平盘，其余均为安全平盘。模型参数定义如下：Ba为工作帮年推进度，m/a；ha为 A 煤厚度，m；hb为 B 煤厚度，m；La为 A 煤底板剥离工作线长度，m；Lb为 B 煤底板剥离工作线长度，m；D 为端帮运输平盘宽度，m；d 为端帮保安平盘宽度，m；c 为最小剥离作业区长度，m；h 为台阶高度，m；B 为组合台阶一次推进距离，m；b 为组合台阶非作业平盘宽度，m；Bp为组合台阶作业平盘宽度，m；n 为组合台阶中台阶数目；为台阶坡面角，（）；M 为台阶组数。以下依据建立的工程模型及定义的模型参数进行分析。2.2循环推进距离及露煤分析2.2.1循环推进距离组合台阶作业的实质为：交替完成将剥离台阶刘光伟等：季节性剥离露天煤矿工作帮形态动态控制及推进强度优化2023年第10期47由工作平盘变换为非工作平盘，非工作平盘变换为工作平盘的作业过程的循环。在此过程当中平盘宽度是动态变化的，组合台阶完成一次循环向前推进距离可表示为式（1），即作业平盘宽度与非作业平盘宽度的差值，循环推进距离的大小取决于组合台阶中作业平盘宽度。B=Bpb（1）2.2.2采场露煤方式与循环推进距离关系采场全年单次露煤时，工作帮循环推进距离为年推进度。采场全年多次露煤时，若每次揭露煤量相同，剥离工作帮每次循环推进距离一致；若每次揭露煤量不同，则剥离工作帮循环推进距离与坑内原煤可采时间有关。设每次揭露煤量为 Mk，k=1，2，s。则揭露原煤量与工作帮循环推进距离的关系为Mk=B(haLpa+hbLpb)（2）式中：Mk为剥离工作帮向前推进距离 B 时采场揭露煤量，t；为原煤容重，t/m；为原煤采出率；Lpa为 A煤平均工作线长度，m；Lpb为 B 煤平均工作线长度，m。露煤可采时间可用式（3）表示，且应满足式（4）中关系。Tk=MkqmNc（3）Tc=sk=1Tk24非等量露煤sTk24等量露煤（4）式中：Tk为露煤可采时间，h；qm为采煤设备实际生产能力，t/（台h）；Nc为采煤设备数量，台；Tc为年采煤作业时间，d。季节性剥离露天矿冬季不进行剥离，坑内可采煤量除正常的接续煤量外，还应留有剥离停滞期储备煤量。因此在实际生产过程中应当在保证采煤接续的情况下优化储备煤量的留设方式，从而尽可能减小每次循环推进工作帮超前剥离量，以实现控制生产剥采比，减小内排运距的目的。23m1112n12n12n12nBpBbhA 煤B 煤hahbbB剥离作业区剥离作业区DdLbLa推进方向剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区剥离作业区BpBbbLGij图2季节性剥离露天煤矿工作帮推进工程模型Fig.2Engineeringmodelforadvancingtheworkingslopeofseasonallystrippedopenpitcoalmines2023年第10期煤炭科学技术第51卷482.3工作帮循环推进时间2.3.1循环推进工程量设工作帮从上至下有多组组合台阶，分别为 i=1，2，m，每组台阶内各分层从上至下编号依次为 j=1，2，n。依据建立的组合台阶工作帮推进工程模型，可计算得出一组组合台阶内单个分层推进一次完成的工程量，如式（5）所示。Qij=BhLGij（5）其中，当前台阶与下部相邻分层剥离工作线长度关系为1）当前分层与端帮运输平盘同水平时：LGijLGij+1=2(D+htan)（6）2）当前分层与端帮保安平盘同水平时：LGijLGij+1=2(d+htan)（7）式中：Qij为第 i 组组合台阶内第 j 个分层向前推进距离 B 完成的工程量，m；LGij为第 i 组组合台阶内第 j个分层的剥离工作线长度，m。2.3.2各分层设备生产能力组合台阶工作帮推进时，每组台阶中仅有一个平盘进行采掘作业，因此生产能力与该平盘设备布置情况有关。已知剥离设备作业区最小长度为 c，则可以计算出平盘上可以布置的最大设备组数如式（8）所示。Ns,max=LGijc（8）式中：Ns为组合台阶作业平盘布置设备组数，组（Ns1，Ns,max）。因此组合台阶作业平盘设备生产能力为Aij=qbNs（9）式中：Aij为第 i 组组合台阶第 j 个作业平盘设备生产能力，m/h；qb为单组剥离设备实际生产能力，m/（组h）。2.3.3循环推进时间组合台阶整体向前推进距离 B 时，完成一次作业循环，相应的采掘设备在组合台阶中各分层之间完成一次走行、调配的循环。设每组台阶完成 1 次循环推进所需时间为 Ti，Ti的计算方法如式（10）所示。Ti=Ti0+Ti1+Ti2+Ti3（10）式中：Ti0为第 i 组组合台阶累计穿爆作业时间，h；Ti1为第 i 组组合台阶各平盘累计采剥作业时间，h；Ti2为第 i 组组合台阶采掘设备在相邻平盘之间走行以及完成一次扩帮作业后回到初始作业水平的时间，h；Ti3为设备停滞时间，h。为了给采掘设备提供有利的工作条件，提高作业效率，需要对坚硬岩石和煤层进行爆破作业。实行爆破作业的露天煤矿，该过程包括 2 个作业环节，即穿孔和爆破。因而穿爆时间可用式（11）表示。Ti0=lzqzNz+t0（11）式中：lz为累计炮孔深度，m；qz为钻机生产能力，m/（台h）；Nz为钻机数量，台；t0为爆破作业时间，包括装药及爆破时间，h。依据开采参数，可以计算出组合台阶进行一次作业循环所完成的工程量，进而得出各平盘设备采剥作业总时间如式（12）所示。Ti1=nj=1QijAij（12）设备在当前台阶完成作业任务后需要移动到下部相邻平盘进行剥离作业，完成一次循环推进后又回到初始作业平盘。因此走行、调配时间可按式（13）计算。Ti2=Ns(t1+t2)（13）t1=h(n1)(v0+v1)1 000Dev0v1（14）t2=2(n1)(B+b+htan)+nj=1LGij4nrk11 000v2（15）式中：t1为设备走行移动坑线所需时间，h；De为移动坑线坡度，8%；v0为采掘设备在移动坑线中向上走行速度，km/h；v1为采掘设备在移动坑线中向下走行速度，km/h；t2为设备平盘走行时间，h；v2为采掘设备在工作平盘上平均走行速度，km/h；k1为路线增长系数；nr为台阶布置坡道数目。为了保证露天煤矿采掘设备能够正常、高效的运转，满足采掘设备的完好率和出动率要求，剥离设备的停机维修和保养至关重要。考虑到设备检修制度和机修设施等因素，设备停滞时间由设备的维修周期和工期决定。Ti3=Ti1+Ti230Txt3+Ti1+Ti230Tyt4（16）式中：Tx为设备小修周期，月；t3为设备小修工期，d；Ty为设备月检周期，月；t4为设备月检工期，d；整理式（16）得 1 组台阶完成 1 次循环推进所需时间为：刘光伟等：季节性剥离露天煤矿工作帮形态动态控制及推进强度优化2023年第10期49Ti=lzqzNz+t0+nj=1BhLGijqbNs+Nsh(n1)(v0+v1)1 000Dev0v1+2(n1)(B+b+htan)+nj=1LGij4nrk11 000v2+Ti1+Ti230Txt3+Ti1+Ti230Tyt4（17）工作帮各组组合台阶均完成 1 次循环推进后，工作帮方可整体向前推进距离 B，以实现采场露煤。因此工作帮完成一次循环推进时间为各组组合台阶完成一次循环推进时间的最大值。即：T=max(T1,T2,Tm)（18）式中：T 为工作帮完成一次循环推进所需时间，h。综上所述，作业平盘布置设备数量越多时，平盘推进时间越短。此外，组合台阶中台阶数目越多时，采掘设备在平盘及移动坑线中走行时间越长，间接的降低了设备的作业效率。因此，需要尽量减小设备非剥离作业时间在一次循环推进时间中的占比。2.4端帮运输道路可靠性煤矿露天开采时，考虑到端帮资源回收、运输系统布置以及边坡稳定性等因素，端帮台阶也常采用组合台阶且与工作帮台阶组合方式不一定相同。因此，工作帮部分平盘产生的剥离物料并不能够通过同水平端帮运输平盘进入内排土场。假设剥离物在工作帮运输过程不存在向下运输，且只经由该平盘上部首个端帮运输平盘运输至内排土场，则通过端帮运输平盘的剥离量由端帮台阶组合方式和工作帮循环推进距离决定。设端帮组合台阶内台阶数目为 nd，端帮道路行车密度计算如式（19）：Nd=BhK3TGK1K2ndLGi*j*2(nd1)(D+htan)2(nd2)(d+htan)（19）式中：Nd为道路行车密度，辆/h；为剥离物比重，t/m；LGi*j*为当前端帮运输平盘同水平剥离工作线长度，m；T为各组组合台阶完成一次循环推进时间的最小值，h；G 为汽车载重，t/辆；K1为时间利用系数，采用三班工作制时为 0.75；K2为汽车载重利用系数，不低于0.9；K3为运输不均衡系数，取 1.11.15。组合台阶循环推进距离增加时，通过端帮运输平盘的剥离量增加，端帮运输道路行车密度加大。为了保证运输道路可靠性，降低汽车运输安全风险，端帮运输道路行车密度不应超过该道路等级规定的最大值。3工作帮形态优化高寒地区季节性剥离露天煤矿冬季不进行剥离作业活动，需要在剥离作业期完成超前剥离任务来保证足够的露煤量以满足采煤接续的要求，与全年均衡剥离作业露天煤矿相比，剥离物端帮运输距离更大。通过优化剥离工作帮形态来实现采剥接续、减小剥离物端帮运距、降低剥离物运输功的目标。3.1优化模型建立3.1.1决策变量选取优化季节性剥离露天煤矿工作帮形态实际上就是优化工作帮台阶组合方式及一次循环作业推进距离，而循环推进距离又取决于作业平盘宽度。确定季节性剥离露天煤矿工作帮形态优化的决策变量为n、Bp。此外工作帮形态还由采掘设备规格、工作面设备布置情况、以及年剥离作业时间等参数共同决定。3.1.2目标函数建立通过优化季节性剥离露天煤矿工作帮形态来达到在保证产能接续、运输系统稳定可靠的基础上剥离物运输功最小的目的。建立目标函数如下：minW=mi=1nj=1ghLGij(Bpb)LGij+LPijK+(Ld+BpbBn)+H1+H2De（20）式中：W 为剥离物运输功，Nm；为加权道路滚动阻力系数；LPij为对应内排平盘工作线长度，m；K 为系数，采用单环内排时，K=2；采用双环内排时，K=4；Ld为工作帮循环推进前，当前台阶剥离物端帮运距，m；Bn为循环推进后内排土场跟进距离，m；H1为当前工作台阶与上部首个端帮运输平盘之间高差，m；H2为端帮运输平盘与内排土场排土平盘之间高差，m。3.1.3约束条件整理1）工作帮循环推进距离：由组合台阶作业方式可知，剥离工作帮循环推进距离应满足B Bp,minb（21）2）剥采工程接续：剥离开工时坑内原煤可采时间应与循环推进时间匹配。Tk1 T（22）3）工作平盘设备布置数量不应超过可布置的最大值1 Ns Ns,max（23）4）端帮运输道路行车密度应小于其道路等级规2023年第10期煤炭科学技术第51卷50定的最大行车密度Nd N（24）式中：Bp,min为最小工作平盘宽度，m；N 为最大行车密度，辆/h；5）其他约束：LGij,Ld,LPij 0（25）分析工作帮形态优化模型目标函数可知，W 为关于 n 和 Bp的非线性函数。通过进行非线性规划求解进而得出设备规格、工作面设备布置情况及年剥离作业时间等参数确定时满足剥采工程时空接续要求且剥离物运输功最小的剥离工作帮台阶组合参数n、Bp。3.2序列二次规划序列二次规划法（SQP）是求解非线性约束优化问题最有效的算法之一17。序列二次规划通过将非线性约束问题转化为简单的二次规划（QP）问题，在每一步迭代中通过求解一个二次规划子问题来确定下一次迭代进行方向，采用线搜索算法确定步长，重复这些步骤直到得到原问题的解18-19。非线性约束数学规划问题的一般形式可以通过式（26）进行表示20：minf(X)s.t.gu(X)0,u=1,2,lhv(X)=0,v=1,2,m（26）f(X)gu(X)hv(X)式中：X 为变量；为目标函数；为不等式约束函数；为等式约束函数；u 为不等式约束函数个数；v 为等式约束函数个数。利用泰勒展开将式（26）中的目标函数转化为二次函数，将约束函数转化为线性函数后得到简单的二次规划问题。二次规划子问题表示如下：minf(X)=12dTkHkdk+f(X)Tdks.t.gu(X)Tdk+gu(X)u=1,2,lhv(X)Tdk+hv(X)v=1,2,m（27）f(X)f(X)gu(X)hv(X)式中：dk为二次规划子问题的最优解，也是迭代进行方向；为的梯度函数；、分别为约束函数的雅各比矩阵；Hk为拉格朗日函数的Hesse 矩阵。原问题的拉格朗日函数为L(X,)=f(X)+gu(X)+hv(X)（28）式中：和 分别为不等式约束与等式约束的乘子。通过一维搜索进行迭代步长确定，步长满足以下条件，Xk+1=Xk+adk（29）式中：Xk为第 k 次迭代值；a 为步长。采用拟牛顿法中 BFGS 公式进行拉格朗日函数Hesse 矩阵计算21，BFGS 公式描述如下。Hk+1=Hk+(1+yTkHkykSTkyk)SkSTkSTkykSkyTkHk+HkykSTkSTkyk（30）Sk=Xk+1Xkyk=L(Xk+1,)L(Xk,)其中：；。序列二次规划法算法求解流程如图 3 所示。线搜索求步长 a更新迭代点 X k+1X k+1 是否满足收敛精度?开始结束Y简化原问题为二次规划子问题N给定 X0,收敛精度；置 k=0,令 H0=I求解二次规划子问题,得到子问题最优解 dk,确定迭代前进方向BFGS 公式修正Hesse 矩阵 Hk+1输出最优解 X*=X k+1 以及目标函数值f*=f(X k+1)图3序列二次规划算法流程Fig.3sequencequadraticprogrammingalgorithmflowchart3.3模型求解实际生产中，为了便于开采设计和工程实现，开采参数常常取整数。因此 n 和 Bp为离散变量。对于离散变量的最优化问题，目前直接的求解算法还不够成熟，通常的处理方法是先将离散变量当作连续变量，用连续变量最优化算法求出连续最优解之后，再作适当的离散化处理，如某种方式的圆整或者取标准值等18。通过给定迭代初值 Bp0=Bpmin、n0=2。根据序列二次规划求解流程对式（20）表述问题进行迭代求解。确定满足约束条件的最优解，将最优解进行取整处刘光伟等：季节性剥离露天煤矿工作帮形态动态控制及推进强度优化2023年第10期51理并比较 W，最终确定最优工作帮形态参数 n、Bp。4工作帮极限推进强度分析露天煤矿正常开采过程中，工作帮各年推进度基本稳定，当采剥工作线缩短、产能增加、煤层赋存条发生变化或者推进至大型地质构造区域时需增大推进度以保证产能的稳定接续。随着当前采掘设备的发展，设备性能大幅提升；采掘设备的小型化，使得采掘设备作业方式更加灵活。在露天煤矿的实际生产中，工作帮可达到的推进度已突破开采设计规范中所建议的推进速度参考值22。特别是对于高寒地区露天煤矿剥离工程非全年进行的情况，确定其在剥离作业期内工作帮可以实现的极限推进度是有必要的。季节性剥离露天煤矿工作帮推进强度受工作帮形态、采掘设备生产能力、工作面设备布置情况以及年剥离作业时间影响。当工作面布置设备数量增加、设备生产能力加大或者年剥离作业时间增加时，工作帮循环推进时间减小，工作帮在剥离作业期内循环推进次数增加，相应的工作帮推进强度加大。工作帮极限推进强度计算如式（31）：Ba,max=BTa,maxTmin（31）式中：Ba,max为当前形态下工作帮极限推进度，m/a；Ta,max为露天煤矿最大年剥离作业时间，h；Tmin为当前工作帮形态下采掘设备生产能力最大时工作帮循环推进时间，h。5工程实例胜利西二号露天煤矿位于内蒙古自治区锡林郭勒盟锡林浩特市西北部宝力根（胜利）苏木境内，设计规模 10Mt/a。采煤工艺为：单斗卡车地面半固定破碎站半连续工艺，剥离为单斗卡车间断工艺。该矿为典型的季节性剥离露天煤矿，剥离作业期在每年的 4 月初至 11 月末。剥离台阶水平划分，剥离模式为外包作业，采用液压反铲直接挖掘的方式进行剥离，沿工作线布置方向水平划分 2 个工程标段，每个工作台阶布置 2 个坡道。剥离物经北端帮运输至内排土场排弃，为单环运输。端帮边坡由运输平盘与安全平盘构成，每隔 1 个运输平盘留设 1 个安全平盘。端帮运输道路为运输支线，属于二级道路。开采境界内 6 号煤层为主采煤层，呈近水平赋存。首采区煤层顶板平均标高+940m，地表平均标高+1096m。采场全年露煤 5 次，每次露煤量 200 万 t，具体指标参数如下：端帮运输平盘宽度D/m40端帮保安平盘宽度d/m5剥离工作帮非作业平盘宽度b/m25剥离台阶高度h/m12剥离台阶坡面角/()65最小剥离作业区宽度c/m250煤层厚度ha/m19.25原煤视密度/(tm3)1.23煤层底板工作线长度La/m1300运岩卡车载重G/t80剥离物比重/(tm3)2.65采掘设备向上走行移动坑线速度v0/(kmh1)2.4采掘设备向下走行移动坑线速度v1/(kmh1)3.5采掘设备平盘移动速度v2/(kmh1)3年剥离作业时间Ta/d240年采煤作业时间Tc/d330剥离设备实际生产能力qb/(万m3台1d1)2.2采煤设备实际生产能力qm/(万m3台1d1)0.76采煤设备数量Nc/台4内排追踪距离lp/m50设备小修周期Tx/d10设备小修工期t3/h1设备月检周期Ty/d20设备月检工期t4/h5工作面剥离设备按最大数量布置时，将上述所列参数代入优化模型中进行求解。结果显示：组合台阶内台阶数目为 5，组合台阶中作业平盘宽度为89m 时，剥离物运输功最小且满足剥采工程时空接续要求。该形态下工作帮循环推进时间为 47.5d，循环推进距离为 64m。端帮运输道路行车密度为141 辆/h，小于厂矿道路设计规范中规定二级道路最大行车密度（170 辆/h）23，计算得工作帮年推进度为 320m/a。工作帮台阶组合方式及循环推进距离变化时，对应极限推进度计算结果分别如图 4、图 5 所示。循环推进距离不变时，组合台阶内台阶数目越多，循环推进时间越长。工作帮在剥离作业期内循环推进次数越少，工作帮年推进度越小。循环推进距离大于 325m 时，工作帮在剥离作业期内无法完成一次循环推进。从图 5 中可以看出，单次循环推进时，工作帮极限推进度最大。多次循环推进时，工作帮极限推进度随循环推进距离的加2023年第10期煤炭科学技术第51卷52大呈线性增加趋势。采场全年多次等量露煤时，工作帮年循环推进次数随循环推进距离的增加呈周期性递减；非等量露煤时，工作帮年推进度为各次循环推进距离的累计值。季节性剥离露天矿实际生产中应根据产能规划和原煤需求以及通过调节工作面设备布置情况等手段综合确定采场露煤方式和循环推进距离，动态控制剥离工作帮形态，优化确定工作帮推进度。6结论1）将露天煤矿工作帮形态控制优化工程问题抽象为非线性规划数学模型，引入序列二次规划算法计算求解决定工作帮形态的关键参数 n、Bp，实现季节性剥离露天煤矿剥离工作帮形态动态控制。2）季节性剥离露天矿采场全年单次露煤时，工作帮极限推进度为循环推进距离。全年多次等量露煤时，工作帮极限推进度随着循环推进距离的变化呈周期性线性变化。3）应用文中提出方法计算确定了胜利西二号露天煤矿满足剥采工程时空接续要求且剥离物运输功最小的组合台阶内台阶数目 n 为 5，作业平盘宽度Bp为 89m，该形态下工作帮极限推进度为 320m/a。参考文献(References)：刘峰，郭林峰，赵路正.双碳背景下煤炭安全区间与绿色低碳技术路径J.煤炭学报，2022，47（1）：115.LIUFeng，GUOLinfeng，ZHAOLuzheng.Researchoncoalsafetyrangeandgreenlow-carbontechnologypathunderthedual-car-bonbackgroundJ.JournalofChinaCoalSociety，20