母傳偉　楊維菁

(中冶沈勘秦皇島工程技術(shù)有限公司礦山設(shè)計院)

?

基于線性規(guī)劃單純形法優(yōu)化礦巖調(diào)運

母傳偉楊維菁

(中冶沈勘秦皇島工程技術(shù)有限公司礦山設(shè)計院)

摘要礦巖調(diào)運是礦山總圖運輸系統(tǒng)的重要組成部分，是構(gòu)成礦山系統(tǒng)工程的關(guān)鍵環(huán)節(jié)，將運籌學(xué)理論和礦巖調(diào)運實際問題相結(jié)合，提出了礦山總圖運輸優(yōu)化的新方法。以某鐵礦為例，將礦山生產(chǎn)運輸實際問題轉(zhuǎn)化為標(biāo)準(zhǔn)化數(shù)學(xué)模型，通過回歸分析求得汽車、鐵路、膠帶機運距和運費的回歸方程，結(jié)合相關(guān)模型參數(shù)，編制運輸平衡表和運價表，按算法程序要求格式構(gòu)造初始單純形表，通過編寫的VB語言計算程序，最終求得最佳礦巖調(diào)運方案。對于礦巖和運輸設(shè)備進(jìn)行統(tǒng)籌安排，為礦山設(shè)計和生產(chǎn)管理找到了成本最低的礦巖運輸方案。

關(guān)鍵詞線性規(guī)劃單純形法礦巖調(diào)運運籌學(xué)礦山系統(tǒng)工程

礦巖調(diào)運是構(gòu)成礦山系統(tǒng)工程的關(guān)鍵環(huán)節(jié)，可根據(jù)露天礦生產(chǎn)規(guī)模，為采礦和剝巖進(jìn)行合理分配運輸，規(guī)劃經(jīng)濟最優(yōu)的流向，使總運輸成本最低。露天礦原礦與巖石運輸規(guī)模巨大，運輸成本占采礦年經(jīng)營成本的30%～50%，降低運費是縮減采礦生產(chǎn)成本的有效途徑。通過科學(xué)方法求得最佳調(diào)運方案，可以最大限度地節(jié)省運輸費用。國內(nèi)大部分露天礦山礦巖生產(chǎn)調(diào)運一般是靠經(jīng)驗，對礦巖的運量和去向分配很少有統(tǒng)一規(guī)劃，隨機性和盲目性大，導(dǎo)致在生產(chǎn)過程中巖石不能按照最優(yōu)線路排到相應(yīng)排土場，巖石運輸功無端增大，造成排巖成本增加，極大地影響企業(yè)經(jīng)濟效益。優(yōu)化礦巖調(diào)運方案，讓礦山生產(chǎn)運輸設(shè)備發(fā)揮最佳效果，對礦山企業(yè)經(jīng)濟意義重大[1-2]。

1　線性規(guī)劃單純形法數(shù)學(xué)模型

單純形法是求解線性規(guī)劃問題的迭代算法，其求解基本思路如圖1。對于一個可行域非空且有最優(yōu)解的線性規(guī)劃問題，通過在可行域的邊界上尋求和變換基本可行解而最終求得最優(yōu)解[3]。

圖1　單純形法求解的基本思路

采用線性規(guī)劃單純形法求解運輸問題時，首先把運輸問題轉(zhuǎn)化成標(biāo)準(zhǔn)矩陣形式的數(shù)學(xué)模型，構(gòu)造(m+1)×(n+1)階矩陣初始單純形表。如果標(biāo)準(zhǔn)化后的目標(biāo)函數(shù)方程右端含有常數(shù)項，在列表時將其從計算一開始就參與矩陣變換。依據(jù)標(biāo)準(zhǔn)化后數(shù)學(xué)模型中的各數(shù)據(jù)所形成的初始單純形表，按表中系數(shù)矩陣中是否包含有一個m階單位子矩陣，再分兩個階段分別計算。鑒別和尋求新的基本可行解是單純形法的核心。當(dāng)鑒別所求的基本可行解為非最優(yōu)解后，就要在此基礎(chǔ)上進(jìn)一步尋求使目標(biāo)函數(shù)值減小的新的基本可行解，通過矩陣的旋轉(zhuǎn)變換，使系數(shù)矩陣中保留一個m階單位子矩陣I0，再令新單純形表中的I0為基，就得一新基本可行解，判斷是否非唯一最優(yōu)解后再求其全部最優(yōu)解。如果已得一個不滿足最優(yōu)性鑒別標(biāo)準(zhǔn)的基本可行解，而又無法再得到一個使目標(biāo)函數(shù)值下降的新基本可行解，即判定該問題有可行解而無最優(yōu)解[4-6]。

2　礦巖調(diào)運方案及相關(guān)參數(shù)確定

2.1礦山運輸問題解決方案

單純形法所需要的數(shù)學(xué)參數(shù)是剝巖量、運輸量和單位運輸成本，要先求出相關(guān)參數(shù)或給出相關(guān)參數(shù)的求解方程，再編制運輸平衡表和運價表，明確生產(chǎn)水平、剝巖量、運輸量和單位成本，將實際生產(chǎn)問題轉(zhuǎn)化成為數(shù)學(xué)模型，根據(jù)給定的約束條件，建立線性規(guī)劃單純形法約束方程，按照其算法步驟編制計算程序，求出方程最優(yōu)解，即可得到總運費最小的礦巖調(diào)運方案。

礦山運輸所要解決的問題是以最低的運輸成本把開采礦石或剝離巖石運往相應(yīng)地點。礦山各開采水平的巖石產(chǎn)量一定，各個排土場的容量一定，而且剝離巖石必須全部運往相應(yīng)排土場，從各開采水平到各排土場的運費和運量呈線性關(guān)系。假定礦山采場有m個開采水平剝離巖石，將巖石運往n個排土場，剝巖量與排棄量相等，屬于運籌學(xué)中的平衡運輸問題。按圖1求解思路把所求解的平衡運輸問題轉(zhuǎn)化為線性規(guī)劃單純性法數(shù)學(xué)模型，求出最優(yōu)解，找到使總運費最小而且滿足剝離巖石與排棄巖石相等的最佳運輸方案。

假設(shè)采場開采水平i運往排土場j的單位運價為Cij，i=1，2，…，m；j=1，2，…,n，Xij為開采水平i至排土場j的最佳運量，Z為最小總運費，則運用線性規(guī)劃原理構(gòu)造數(shù)學(xué)模型。

目標(biāo)函數(shù)

(1)

約束條件1

(2)

約束條件2

(3)

非負(fù)條件

(4)

2.2有關(guān)參數(shù)的確定

線性規(guī)劃單純形法求解的關(guān)鍵步驟是根據(jù)采場開拓運輸方式和運輸排土方式編制平衡表和運價表，相關(guān)計算參數(shù)的計算是求解的關(guān)鍵，各參數(shù)應(yīng)根據(jù)其影響因素分別求解。

2.2.1采場汽車運距的確定

采場汽車運輸距離分為固定線和移動線兩部分，其表現(xiàn)形式都是具體數(shù)據(jù)陣列。固定線運距可圖上或現(xiàn)場實測，移動線運距則與諸多因素有關(guān)，如采礦方法是橫向推進(jìn)還是縱向推進(jìn)、采用何種運輸方式、各開采年份的采剝量等，需根據(jù)每年及上一年的采場設(shè)計年末圖進(jìn)行計算和量測。采場移動線運距的確定工作復(fù)雜繁重，數(shù)據(jù)量大，對計算過程和最終結(jié)果影響最大。

2.2.2汽車效率回歸方程

根據(jù)運輸距離，結(jié)合現(xiàn)場運輸設(shè)備型號規(guī)格、設(shè)備性能、道路坡度、路況以及現(xiàn)場實際統(tǒng)計的樣本資料，求出汽車運輸效率與運距的回歸方程。采用Excel的回歸分析功能，將收集到的汽車生產(chǎn)月報表相對于運距的現(xiàn)場樣本數(shù)值輸入Excel表格，做出散點圖，根據(jù)數(shù)據(jù)與散點圖進(jìn)行回歸分析，得到汽車效率與運距的回歸方程曲線，見圖2。

圖2　汽車運距-運輸效率回歸方程

2.2.3汽車運費回歸方程

汽車運費指標(biāo)應(yīng)根據(jù)運輸距離，結(jié)合現(xiàn)場運輸設(shè)備型號規(guī)格、設(shè)備性能、耗油指標(biāo)、輪胎消耗、運輸線路情況、當(dāng)前所燃材料和人工費用價格等因素，根據(jù)現(xiàn)場實際統(tǒng)計的樣本資料，求出汽車運費與運輸距離的回歸方程。根據(jù)數(shù)據(jù)與散點圖進(jìn)行回歸分析，得到汽車運費與運距的回歸方程式5，汽車運距-運費回歸方程曲線見圖3。

圖3　汽車運距-運費回歸方程曲線

85 t礦用電動輪自卸汽車的運費回歸方程：

R1=2.42-0.16L+0.036L2-0.003L3,

(5)

式中,L為汽車運距，km；R1為汽車運費，元/(t·km)。

2.2.4鐵路運費回歸方程

同理，結(jié)合現(xiàn)場鐵路機車和車輛設(shè)備型號規(guī)格、性能、線路坡度等各種因素，根據(jù)現(xiàn)場實際統(tǒng)計的樣本資料，求出鐵路運費與運輸距離的回歸方程及曲線圖(圖4)?；貧w方程為

R2=L/(2.185 97L-0.289 95) ,

(6)

式中,L為鐵路運距，km；R2為鐵路運費，元/(t·km)。

圖4　鐵路運距-運費回歸方程曲線

2.2.5膠帶機運費回歸方程

結(jié)合現(xiàn)場膠帶機設(shè)備型號規(guī)格、性能、耗電指標(biāo)、材料消耗、膠帶機坡度等各種因素，根據(jù)現(xiàn)場實際統(tǒng)計的樣本資料，求出膠帶機運費與運輸距離的回歸方程及曲線圖(圖5)。

圖5　膠帶機運距-運費回歸方程曲線

回歸方程

(7)

式中,L為膠帶機運距，km；R3為膠帶機運費，元/(t·km)。

2.2.6運輸平衡表和運價表的編制

當(dāng)采場各水平汽車運距和各種運輸方式的單位巖石量運輸經(jīng)營費確定后，可編制巖石運輸平衡表。平衡表上的各行表示采場各開采水平編號、各水平的汽車運距及其巖石量，各列表示去各排土場的運輸方式或排土場編號，最低行各列為每個排土場的容量。

運距表編制的主要依據(jù)是計算運距和各運輸方式的運輸經(jīng)營費，單位經(jīng)營費Cij應(yīng)根據(jù)所采用的設(shè)備情況及經(jīng)營費指標(biāo),分別計算汽車、鐵路和膠帶機經(jīng)營費指標(biāo)，應(yīng)考慮資金時間價值的影響。按以上方法計算的Cij為靜態(tài)值，沒考慮費用的時間價值，應(yīng)進(jìn)行修正。而各項費用應(yīng)按動態(tài)計算，即計算每一項費用時，考慮到該項費用的支出時間，將該年的費用支出貼現(xiàn)到計算基準(zhǔn)年，即

(8)

(9)

式中,NPV為凈現(xiàn)值;D為礦山總壽命;Fy為第r年的現(xiàn)金流量;Mr為一次償付復(fù)利系數(shù);1/Mr為一次償付現(xiàn)值因數(shù);P8為投資基準(zhǔn)收益率;Er為開采年份。

按通常的設(shè)計程序在巖石運輸系統(tǒng)設(shè)計時，一般不能給出開采進(jìn)度計劃，僅能確定巖石產(chǎn)量及各水平的分層巖石量，所以，要精確計算Er是比較困難的，只能根據(jù)已知的分層巖石量近似計算各水平的開采年份Er。

(10)

式中,Er、Er-1分別為r水平和r-1水平的開采年份；Xr為r水平的分層巖石量，萬t;P3為巖石年產(chǎn)量，萬t/a。

修正后的運價為：Cij=Cij/Mr(元/t)。

3　實例分析

3.1礦山生產(chǎn)運輸情況簡介

某鐵礦擴建工程設(shè)計中，礦石規(guī)模1 800萬t/a，采場剝離巖石3 700萬t/a，設(shè)計4種排土方式，即河?xùn)|鐵路排土(1#排土場)、河?xùn)|膠帶排土機排土(2#排土場)、河西膠帶排土機排土(3#排土場)和汽車排土(4#排土場)。各巖石運輸排土方式：①全汽車運輸，運輸經(jīng)營費包括汽車運費和推土機排土費用；②汽車轉(zhuǎn)鐵路運輸，運輸經(jīng)營費包括汽車運費、鐵路運費、汽車轉(zhuǎn)鐵路轉(zhuǎn)載費和排土場電鏟倒裝費；③汽車轉(zhuǎn)膠帶機運輸，運輸經(jīng)營費包括汽車運費、膠帶機運費、破碎費用和排土機排土費。露天采礦場多個臺階同時生產(chǎn)，將露天采場各開采水平的巖石運到相應(yīng)排土場堆排，尋求使整個露天礦巖石運輸系統(tǒng)費用總和最低，即最優(yōu)的巖石調(diào)配方案。

3.2建立運輸平衡表

各排土場排棄的巖石量等于采場各生產(chǎn)水平剝離的巖石量。在程序中輸入露天采場各開采水平剝離的巖石量、各巖石運輸方式接收的巖石量，建立運輸平衡表，見表1。

表1　運輸平衡表

3.3構(gòu)造約束方程組系數(shù)矩陣

約束方程組的系數(shù)為式(11)，系數(shù)矩陣形式見表2。矩陣中每個元素非0即1，且每一列元素只有兩個1，其余均為0，前m行之和等于后n行之和，矩陣的秩為m+n+1。

(11)

3.4計算單位巖量運輸經(jīng)營費

汽車、鐵路和膠帶機運輸單位巖量經(jīng)營費計算[7]：

(12)

各噸位汽車、鐵路、膠帶機運費及推土機、電鏟倒裝費計算[8-20]:

(13)

式中,Cq、Ct、Cj分別為采用汽車、鐵路、膠帶運輸方式巖石單位運費；Rq、Rt、Rj分別為汽車、鐵路、膠帶運輸噸公里費用；Ft(0.6元/t)、Fd、Fp(0.854元/t)、Fp(0.5元/t)分別為單位礦巖推土機排土費、電鏟倒裝費、破碎費、排土機排土費；L、Lt、Lj分別為汽車、鐵路、膠帶運距；R118、R85、R77、R42分別為118，85，77，42t汽車運費；Rt、Rj分別為鐵路、膠帶機運費；F1、F2分別為電鏟第一次、第二次倒裝費，F(xiàn)1=F2=1.372元/t。

采場內(nèi)汽車固定線運距可以從設(shè)計運輸線路圖中量取，移動線運距采用加權(quán)平均，從采場年末圖中分臺階量取后經(jīng)計算求得。把巖石單位運輸經(jīng)營費和運輸平衡表合為表3。

3.5計算及結(jié)果分析

通過程序運行計算，計算結(jié)果經(jīng)檢驗無誤即得最優(yōu)運輸方案，采場各水平巖石量以最優(yōu)(即總運費最低)方式分別排放到相應(yīng)的排土場。用VB語言編制的計算程序與用Excel軟件無縫結(jié)合，利用Excel軟件的強大數(shù)據(jù)庫處理和分析功能，使原始數(shù)據(jù)輸入和計算結(jié)果輸出通過Excel給出分析圖線，通過圖線分析，可以直觀地看出各年度數(shù)據(jù)。如哪一年計算出現(xiàn)偏差，可直接顯示出來，然后再去調(diào)整程序或者檢查原始數(shù)據(jù)的問題。

表3　巖石單位運輸經(jīng)營費和運輸平衡表

4　結(jié)　論

(1)礦山運輸問題存在多個可行運輸方案，利用線形規(guī)劃單純形法原理建立數(shù)學(xué)模型，編制計算程序求解，是優(yōu)化礦巖調(diào)運的一條實用方法。借助單純形法計算程序找到礦巖調(diào)運最佳方案，最大限度地節(jié)約運輸費用。

(2)計算程序求解速度快，精度高，省時省力，還可以根據(jù)現(xiàn)場實際參數(shù)的變化及時調(diào)整，能夠有效適應(yīng)礦山實際生產(chǎn)條件的變化。

(3)線性規(guī)劃單純形法也是為生產(chǎn)管理者服務(wù)的決策方法，在解決包含多個相關(guān)約束變量問題時給出最優(yōu)決策。在給定資源條件下如何發(fā)揮最大經(jīng)濟效益的問題，也可采用此方法求解。

參考文獻(xiàn)

[1]母傳偉.宏大鐵礦露天采場開拓運輸優(yōu)化設(shè)計與研究[J].金屬礦山,2014(12):26-29.

[2]母傳偉.影響因子量化加權(quán)法優(yōu)選選礦廠廠址和運輸方式[J].金屬礦山,2014(10):27-30.

[3]云慶夏.采礦系統(tǒng)工程[M].西安：陜西科學(xué)技術(shù)出版社，1992.

[4]桂中岳.礦業(yè)運籌學(xué)[M].北京：冶金工業(yè)出版社，1978.

[5]李寶祥.金屬礦床露天開采[M].北京：冶金工業(yè)出版社，1991.

[6]李德,錢頌迪.運籌學(xué)[M].北京：清華大學(xué)出版社,1982.

[7]母傳偉.線性規(guī)劃單純形法解決礦山礦巖調(diào)運問題應(yīng)用實例[J].金屬礦山,2006(8):122-123.

[8]楊玉新.深凹露天礦提高運輸系統(tǒng)及關(guān)鍵技術(shù)探討[J].礦業(yè)工程,2004(10):28-30.

[9]趙娜.單純形法解線性規(guī)劃問題的算法探究[J].吉林廣播電視大學(xué)學(xué)報,2011(3):112-115.

[10]包麗君.基于線性規(guī)劃單純形法計算運輸問題最優(yōu)解的研究[J].寧波廣播電視大學(xué)學(xué)報,2012(2):128-130.

[11]劉桂蓮.如何在Excel中建立并求解線性規(guī)劃模型[J].承德民族職業(yè)技術(shù)學(xué)院學(xué)報,2005(2):76-79.

[12]李事捷.線性規(guī)劃單純形法在采礦計劃中的應(yīng)用[J].金屬礦山,2009(4):20-25.

[13]侯煜.線性規(guī)劃單純形法在露天礦配礦領(lǐng)域的應(yīng)用[J].礦業(yè)工程,2006(4):22-24.

[14]包麗君.線性規(guī)劃單純形法在物資調(diào)運問題中的應(yīng)用[J].寧波廣播電視大學(xué)學(xué)報,2010(6):127-128.

[15]陳冬英.線性規(guī)劃單純形法在運輸問題中的應(yīng)用[J].科技廣場,2008(11):41-44.

[16]竇永平.線性規(guī)劃模型建立的教學(xué)思路[J].蘭州商學(xué)院學(xué)報,1990(2):66-69.

[17]耿德志.線性規(guī)劃求解算法研究[J].軟件導(dǎo)刊,2011(6):48-50.

[18]何季陽.線性規(guī)劃問題矩陣調(diào)配求解法初探[J].吉林工學(xué)院學(xué)報,1983(4):6-15.

[19]胡成江.運用線性規(guī)劃單純形法實現(xiàn)產(chǎn)品的最優(yōu)組合[J].企業(yè)工作研究,2007(5):69-70.

[20]袁祜田.運用線性規(guī)劃單純形法進(jìn)行土方調(diào)運[J].工業(yè)建筑,1982(4):8-22.

(收稿日期2016-05-13)

Optimization of Allocation and Transportation of Ores and Rocks Based on the Simplex Method of Linear Programming

Mu ChuanweiYang Weijing

(Mine Design Institute,Shenkan Qinhuangdao Engineering & Technology Co.,Ltd.,MCC Group)

AbstractAllocation and Transportation of ores and rocks is the important component of the mine general transportation system and it is the key link in the process of mining complex system.A new optimization design method of mining general transportation is proposed based on operational research theory and the actual situation of the allocation and transportation of ores and rocks.Taking a iron mine as the research example,the actual problem of mine production transportation is converted to standardized mathematical model,the regression equation of the transportation distance and transportation costs of automobile,railway and sealing-tape machine is obtained.Combing with the parameters of related models,the transportation balance sheet and transportation costs sheet are complied,the initial simplex tableau is complied based on the format that are required by algorithm program.The final optimal dispatching scheme of ores and rocks is obtained by using the designed VB language computation program.The ores and rocks and transportation devices are given overall consideration to provide the ores and rocks transportation scheme with the lowest costs for mine design and production management.

KeywordsSimplex method of linear programming,Transportation and allocation of ores and rocks,Operations research,Mine system engineering

母傳偉(1970—)，男，副院長，總工程師，教授級高級工程師，066001 河北省秦皇島市。