1、基金項目：國家自然科學基金面上項目（51774228, 52074205）；陜西省自然科學基金杰青項目(2020JC-44)作者簡介：顧清華(1981- )，男，山東諸城人，教授，博士，主要從事露天礦智能開采、多目標進化算法等研究。E-mail：?；谶M化算法MOEA/D-AU的露天礦多金屬多目標智能配礦研究顧清華1,2，劉思魯1,2，張金龍3（1 西安建筑科技大學 資源工程學院，陜西 西安 710055；2 西安市智慧工業(yè)感知計算與決策重點實驗室，陜西 西安 710055 3 洛陽欒川鉬業(yè)集團股份有限公司，河南 洛陽 471500）摘要：針對露天礦的多金屬多目標配礦問題，提出了基于多目標進化

2、算法的配礦優(yōu)化方法。根據(jù)礦山實際情況，以生產成本、礦石品位偏差和礦石巖性配比偏差最小為優(yōu)化目標建立了露天礦配礦優(yōu)化模型；在基于分解的多目標進化算法（MOEA/D）的基礎上，對算法的更新過程進行了改進，利用種群與權重向量之間的空間位置關系提出了基于角度的更新策略，使算法在求解多目標問題時更好地平衡種群的多樣性與收斂性；由于對選礦因素考慮不充分，無法有效提高礦石的綜合回收率，本文建立了融合氧化率及有害物質參數(shù)的綜合回收率隨機森林預測模型，通過預測模型對算法得到的多組配礦結果進行篩選，獲得一組更加貼合礦山實際情況的配礦計劃。最后以國內某大型鉬鎢銅礦為例進行仿真實驗，實驗結果表明：該配礦計劃在解決多金

3、屬多目標配礦問題時能夠有效提高礦石綜合利用率和企業(yè)經濟效益。關鍵詞：露天礦；配礦優(yōu)化模型；多目標進化算法；綜合回收率；隨機森林中圖分類號：TD854 文獻標志碼：AResearch on Multi-metal and Multi-objective Intelligent Ore Blending in Open-pit Mine Based on Evolutionary Algorithm MOEA/D-AUGU Qinghua1,2, LIU Silu1,2, Zhang Jinlong3(1. School of Resources Engineering, Xian Univers

4、ity of Architecture and Technology, Xian 710055, China; 2. Xian Key Laboratory of Smart Industry Perception Computing and Decision making, Xian 710055, China；3 Luoyang Luanchuan Molybdenum Industry Group Co., Ltd., Luoyang 471500, China)Abstract：Aiming at the problem of multi-metal and multi-objecti

5、ve ore blending in open pit, an optimization method of ore blending based on multi-objective evolutionary algorithm is proposed. According to the actual situation of the mine, the optimization model of open-pit ore blending is established with the minimum production cost, ore grade deviation and ore

6、 lithology ratio deviation as the optimization objectives; Based on the decomposition based multi-objective evolutionary algorithm (MOEA/D), the updating process of the algorithm is improved, and an angle based updating strategy is proposed by using the spatial position relationship between the popu

7、lation and the weight vector, so that the algorithm can better balance the diversity and convergence of the population when solving multi-objective problems; Due to insufficient consideration of mineral processing factors, the comprehensive recovery rate of ore can not be effectively improved. In th

8、is paper, a random forest prediction model of comprehensive recovery rate integrating oxidation rate and harmful material parameters is established. Through the prediction model, multiple groups of ore blending results obtained by the algorithm are screened to obtain a group of ore blending plan mor

9、e suitable for the actual situation of the mine. Finally, taking a large molybdenum-tungsten copper mine in China as an example. The experimental results show that the ore blending plans can effectively improve the comprehensive utilization rate of ore and the economic benefits of enterprises when s

10、olving the problem of multi-metal and multi-objective ore blending.Key words：Open-pit mine, Optimization model of ore blending, Multi-objective evolutionary algorithm, Comprehensive recovery rate, Random forest0 引 言礦產資源作為我國國民經濟的支柱產業(yè)，是社會和經濟發(fā)展的物質基礎。隨著礦產資源利用率降低，品位不斷下降等問題的出現(xiàn)，堅持貧富兼采、合理配礦是提高礦產資源利用率，保證礦石質量

11、的關鍵。由于露天礦山生產條件的日趨復雜，多目標配礦優(yōu)化是大型露天礦山研究的趨勢，編制科學合理的多目標配礦計劃可以有效提高配礦效率，保證礦山企業(yè)可持續(xù)發(fā)展。因此根據(jù)礦山實際生產需求，深入研究多目標配礦問題十分有必要。目前，配礦問題的研究主要集中在計劃模型的構建和優(yōu)化求解算法兩個方面。針對計劃模型的構建而言，當前研究主要側重于配礦指標的選取，根據(jù)不同礦山的生產需求，將礦石品位偏差、生產成本、總產量作為首要考慮指標的研究居多，其中，王李管等1以品位偏差最小為目標函數(shù)建立露天礦配礦優(yōu)化模型，有效解決了多元素多卸礦點配礦而造成的配礦結果粗略等問題；胡乃聯(lián)等2從生產成本的角度出發(fā)，以采掘和運輸成本最小為目

12、標函數(shù)編制露天礦配礦計劃，為礦山實際生產提供依據(jù)；黃啟富等3以總利潤最大化為目標，動態(tài)優(yōu)化解決企業(yè)配礦問題。此外，關于配礦指標的選取，還有對任務量、礦石巖性、氧化率及有害物質方面的研究，其中，柯麗華等4根據(jù)礦石質量要求，以SiO2等有害物質為約束條件構建配礦數(shù)學模型，有效提高了資源利用率；顧清華等5考慮影響選礦回收率的氧化率、有害物質等指標，構建基于采選流程的多金屬配礦優(yōu)化模型。 針對優(yōu)化求解算法而言，主要側重于數(shù)學規(guī)劃、計算機技術與數(shù)學規(guī)劃相結合和智能優(yōu)化算法三個方面，其中，侯煜6針對某露天礦的實際生產，采用線性規(guī)劃原理編制配礦的總體策略與簡化策略，有效解決傳統(tǒng)人工配礦盲目性強且效果較差的問

13、題；劉文博等7從供應鏈角度出發(fā)，分析原礦配礦的生產特點，構建露天礦山供應鏈問題的0-1混合整數(shù)規(guī)劃模型；Eivazy H等8構建了基于混合整數(shù)規(guī)劃的露天礦短期生產計劃模型，提高了生產決策效率；井石滾等9應用先進的地理信息（GIS）技術、全球衛(wèi)星定位（GPS）技術及通用無線分組傳輸（GPRS）技術，設計與實現(xiàn)了露天礦山配礦生產動態(tài)管理系統(tǒng)；吳麗春等10應用0-1整數(shù)規(guī)劃與三維礦業(yè)軟件DIMINE相結合對礦山的配礦進行優(yōu)化，上述研究均為單目標規(guī)劃，但在實際配礦生產作業(yè)中，往往含有多個復雜沖突目標，此時相比于傳統(tǒng)求解方法，智能優(yōu)化算法在求解此類問題上具有明顯優(yōu)勢，其中，李寧等11針對低品位礦產資源的

14、利用率，構建了多目標配礦優(yōu)化模型，并提出求解該模型的混合粒子群算法；顧清華等12針對露天礦多金屬多目標短期配礦問題，建立了基于自適應粒子群算法的配礦優(yōu)化模型。綜上所述，目前配礦研究已經從單目標向多目標發(fā)展，但是對于多目標配礦問題求解而言，大多數(shù)研究還是采用理想點等數(shù)學方法將多目標轉化為單目標，這樣并不能在本質上很好地解決多目標優(yōu)化問題。此外以往研究對選礦因素考慮不充分，僅是將氧化率和有害物質等指標簡單的添加到目標和約束中，但實際生產中上述指標對選礦回收率的影響程度不同，需根據(jù)礦山實際情況具體分析，才能夠顯著改善礦石的綜合回收率。因此本文針對上述問題，在基于分解的多目標進化算法（Multi-ob

15、jective Evolutionary Algorithm Based on Decomposition, MOEA/D）的基礎上，提出了一種改進的多目標進化算法（MOEA/D-AU），在求解多目標配礦優(yōu)化問題上表現(xiàn)優(yōu)異；根據(jù)實際礦山數(shù)據(jù)，利用隨機森林算法建立了綜合回收率預測模型，對進化算法求解的多組配礦方案進行篩選，從而得到一個更加符合實際情況的方案。1 配礦問題描述及模型構建1.1 問題描述露天礦山企業(yè)的實際生產中，配礦作為企業(yè)采選工作的中間環(huán)節(jié)，將不同品位的礦石進行搭配混勻，增加滿足質量標準的礦石產出率。實際生產中，各個出礦點的品位都不相同，即使同一出礦點品位波動也比較大，因此根據(jù)各出

16、礦點的品位變化進行合理配礦至關重要。已知各出礦點的實際品位、各受礦點的目標品位、各出礦點到受礦點的運輸距離以及各出礦點和受礦點的任務量和最大生產能力。多目標優(yōu)化配礦問題可以描述為在實現(xiàn)礦石品位偏差最小、巖性配比偏差最小和生產成本最低的情況下，提高配礦效率，此模型的生產成本包括開采成本、運輸成本及破碎成本。露天礦多目標配礦模型的部分參數(shù)及變量定義如表1所示。表1 基本符號說明Table 1 Description of basic symbols基本符號說明k (k = 1,2,w)礦石中金屬k的種類m (i = 1,2,m)出礦點的數(shù)量n (j = 1,2,n)受礦點的數(shù)量xij表示從出礦點i

17、到受礦點j的運礦量gki表示第i個出礦點礦石k的供礦品位Gkj表示第j個受礦點礦石k的目標品位P1i表示第i個出礦點每噸礦石的開采成本P2ij表示出礦點i到受礦點j的單位運距成本P3j表示第j個受礦點每噸礦石的破碎成本Lij表示從出礦點i到受礦點j的運輸距離Ai表示第i個出礦點的最小任務量Bi表示第i個出礦點的最大生產能力aj表示第j個受礦點的最小任務量bj表示第j個受礦點的最大生產能力C表示出礦點到受礦點的最小運輸量D表示出礦點到受礦點的最大運輸量1.2 模型構建多目標優(yōu)化配礦是一個高復雜性、多約束性的問題，模型根據(jù)礦山實際配礦情況，綜合考慮了礦石品位、生產成本、巖性配比、作業(yè)生產能力以及作

18、業(yè)任務量等多種約束條件，具體模型構建如下： (1) (2) (3)式中，目標函數(shù)（1）表示礦石品位偏差最小化；目標函數(shù)（2）表示采掘、運輸及破碎總成本最小化；目標函數(shù)（3）表示礦石巖性配比偏差最小化，其中,表示礦石巖性，,為期望巖性比例，此外模型將屬于同一巖性的出礦點歸為一類，按不同巖性類別分類排序，即出礦點1到出礦點全部屬于類巖石，出礦點+1到出礦點全部屬于類巖石，以此類推。式中，約束條件（4）表示出礦點最小任務量；約束條件（5）表示出礦點最大生產能力；約束條件（6）表示受礦點最小任務量；約束條件（7）表示受礦點最大生產能力；約束條件（8）出礦點到受礦點的最小、最大運輸量要求。2 多目標進化

19、算法2.1 基于分解的多目標進化算法通過模型構建，可以發(fā)現(xiàn)礦山實際配礦生產作業(yè)為多目標優(yōu)化問題（Multi-objective Optimization Problems, MOPs），該問題各目標之間通常存在沖突，無法同時得到所有目標的最優(yōu)解，而是獲得一組權衡各目標后的折衷解，即Pareto最優(yōu)解集（Pareto Optimal Set, PS）。在求解多目標優(yōu)化問題方面，多目標進化算法（Multi-objective Evolutionary Algorithms, MOEAs）憑借在一次運行下可以獲得一組Pareto最優(yōu)解集的優(yōu)勢而得到廣泛研究，有效解決了傳統(tǒng)多目標優(yōu)化算法魯棒性差，計算

20、效率低等劣勢13。目前MOEAs大致分為三類，不同于基于帕累托支配關系的MOEAs和基于指標的MOEAs這兩類，基于分解的MOEAs在維持種群收斂性與多樣性方面表現(xiàn)優(yōu)越，從而被廣泛應用于求解MOPs，其中最經典的是張青富等14提出的MOEA/D算法。該算法是一種通過聚合函數(shù)將多目標優(yōu)化問題分解成若干個單目標優(yōu)化子問題的方法，常用的分解策略為切比雪夫聚合法（Tchebycheff Approach, TCH），TCH是一種非線性多目標聚合方法，其聚合函數(shù)定義如下： (9)2.2 MOEA/D-AU算法2.2.1 改進動機理想情況下，當MOEA/D算法利用公式（9）定義的每個子問題都取得最優(yōu)解時，

21、種群最終才能均勻的分布在Pareto前沿（Pareto Optimal Front, PF）。但是實際情況并非如此，由于聚合函數(shù)等高線自身的劣勢，導致部分子代解在迭代過程中偏離自己的權重向量方向，從而無法較好的維持種群的多樣性。如圖1所示，二維目標空間下，在6個權重向量（w1, w2, , w6）輔助下分別獲得解（A-F）,圖中虛線分別為解E和解G所對應權重向量w5的聚合函數(shù)等高線。從圖中可以看出，解集并沒有像權重向量一樣分布均勻，尤其是解B和E，雖然取得了不錯的聚合函數(shù)值，但是偏離了它們相應的權重向量（w2和w5）。這是由于在更新過程中，MOEA/D算法僅僅考慮了聚合函數(shù)值，沒有考慮空間位置

22、關系。由圖中聚合函數(shù)等高線可以發(fā)現(xiàn)，解G的聚合函數(shù)值要劣于解E，因此解G被遠離權重向量w5的解E所替換，從而w5所對應的解變?yōu)榻釫。然而,從綜合考慮種群多樣性與收斂性的角度而言，解G要優(yōu)于解E。此外在迭代初期，種群通常遠離PF，更容易出現(xiàn)選擇誤差，這很有可能將搜索局限到PF的某一部分。圖1 二維目標空間下解分布情況示意圖Fig 1 Illustration of the distribution of solutions in the 2-D objective space針對此問題，本文設計了一種基于角度的更新策略，提出了MOEA/D-AU算法。2.2.2 基于角度的更新策略在上述問題的啟發(fā)

23、下，本文在更新過程中不僅僅考慮解的聚合函數(shù)值，還考慮解與所有權重向量的空間位置關系，為此提出了一種基于角度的更新策略，該策略如下：首先每一個子代解都可以根據(jù)公式（10）計算出與所有權重向量之間的夾角余弦值，式中表示F(x)與w之間的銳角，然后找出V個最大的夾角余弦值，使子代解與V個最大夾角余弦值所對應的解進行聚合函數(shù)值的比較。 （10）如圖2所示，假設V為3，子代解Y通過與所有權重向量之間的夾角余弦值進行比較，找到離它最近的3個權重向量w3,w4,w5，Y在更新過程時只需要和w3,w4,w5所對應的解C,D,E進行聚合函數(shù)值的比較即可，由圖2可以發(fā)現(xiàn)，在空間位置的限制下，使子代解Y無法和偏離權

24、重向量較遠的解進行比較，有效避免了聚合函數(shù)本身的劣勢。采用基于角度的更新策略，相比于MOEA/D算法的更新過程，種群的更新不僅僅基于解之間的聚合函數(shù)值，還考慮了空間位置關系，從而使MOEA/D-AU算法在平衡種群多樣性與收斂性方面取得很大改善。圖 2 二維目標空間下解與權重向量間的夾角示意圖Fig 2 Illustration of the acute angle from the solution to the weight vector in the 2-D objective space2.2.3 MOEA/D-AU算法流程為了便于理解，下面給出了本文算法的流程圖。圖3 MOEA/D-A

25、U算法流程圖Fig 3 Flow chart of MOEA/D-AU algorithm2.2.4 MOEA/D-AU與MOEA/D的對比分析為了測試本文算法MOEA/D-AU的性能，將其與原始算法MOEA/D在標準測試問題上進行仿真實驗。WFG系列問題是MOEA研究領域廣泛使用的代表性測試問題，本實驗在WFG系列問題中隨機選取WFG4、WFG6及WFG8為測試問題。此外，本實驗選取國內外通用的多樣性和收斂性評價標準即純分布性PD（Pure Diversity）和世代距離GD（Generational Distance）來衡量算法求解性能15，其中就PD指標而言，數(shù)值越大表示算法性能越優(yōu)異，

26、相反對于GD指標，數(shù)值越小表示算法性能越優(yōu)異。為了更加形象直觀地表示MOEA/D-AU相比于MOEA/D的優(yōu)勢，圖4可視化了PD指標和GD指標下3個目標WFG4、WFG6、WFG8問題的盒圖。其中，盒圖的上下兩條線為樣本的上下四分位數(shù)，中間的線為中位數(shù)，上下虛線的頂端為最大值和最小值，“+”為異常值。由圖4可以發(fā)現(xiàn)，在PD指標盒圖方面，MOEA/D-AU不僅在3個目標WFG4、WFG6、WFG8問題上取得了最優(yōu)的PD值，多樣性方面優(yōu)于MOEA/D，而且發(fā)現(xiàn)該算法的盒子長度是最短的，即算法所對應的PD指標四分位距離是最小的，即該算法的PD指標的最小值和最大值相差不大，表明MOEA/D-AU算法求

27、出的解集整體質量更高，有更好的穩(wěn)定性；在GD指標盒圖方面，MOEA/D-AU在3個目標WFG4、WFG6、WFG8問題上均取得了最優(yōu)的GD值，收斂性方面明顯優(yōu)于MOEA/D。圖4 3個目標WFG4、WFG6、WFG8測試問題上的PD及GD指標盒圖（1: MOEA/D, 2: MOEA/D-AU）Fig 4 The box plots of PD and GD indicators for the three targets WFG4, WFG6 and WFG8 (1:MOEA/D, 2:MOEA/D-AU)3 綜合選礦回收率預測模型3.2 隨機森林算法隨機森林（Random Forest,

28、RF）是在應用統(tǒng)計學采樣原理的基礎上將若干個決策樹集成在一起的一種算法，算法通過引入決策樹的子模型數(shù)和單顆決策樹的最大特征數(shù)兩項重要參數(shù)有效提高了模型的抗噪能力，較好地克服了過擬合問題。此外隨機森林算法即可以應用在分類問題中，也可以應用于數(shù)據(jù)的回歸分析，同時憑借準確率高，參數(shù)調整簡單等優(yōu)勢，使其相比較于SVM，BP神經網絡等機器學習的方法，在默認參數(shù)下就具有較好的回歸預測效果16。隨機森林算法步驟如下：在包含M個樣本的原始訓練集中有放回的隨機抽取M個樣本，進行n次得到所有的采樣數(shù)據(jù)集。對所有采樣數(shù)據(jù)集分別構建完全生長的決策樹進行訓練。最終結果根據(jù)多個決策樹以投票（分類問題）或求平均值（回歸問題

29、）的方式來計算測試樣本的預測值。3.2 隨機森林預測模型構建通過閱讀有關選礦回收率的大量文獻可以發(fā)現(xiàn)，嵌布粒度、氧化率和有害物質等指標和選礦回收率之間存在密切關系，因此將上述指標作為訓練模型的輸入參數(shù)，通過構建的隨機森林模型進行訓練來得到選礦回收率的預測值，模型結構如圖5所示。圖5 選礦回收率預測模型Fig 5 Prediction model of mineral processing recovery4 工程應用以國內某大型金屬露天礦為例進行分析研究。礦區(qū)面積2.01平方公里，包含13個出礦點和3個受礦點，其中礦體中金屬主要有鉬、鎢、銅三種，有害物質主要為二氧化硫，巖性主要有矽卡巖、透輝石

30、長和長英角巖3種。編制合理的配礦計劃，以每個出礦點到受礦點的出礦量為自變量進行算法優(yōu)化，在滿足各個約束條件的同時使總目標偏差最小。礦山每班（每天3班）有1013臺電鏟工作，有3個破碎站進行破碎。礦山配礦的相關參數(shù)如表2、表3、表4、表5所示。表2 配礦相關參數(shù)Table 2 Ore blending related parameters出礦點原礦品位（%）最小任務量（/t）最大生產能力（/t）MoWO3Cu10.0590.0770.005600200020.1560.1500.004600200030.1180.0790.013600200040.0730.0410.005600200050.

31、0800.0760.013600200060.1390.1020.004600200070.1860.1410.011600200080.1010.1010.004600200090.0660.0620.0116002000100.1000.0800.0116002000110.1050.0810.0126002000120.0660.0810.0306002000130.1070.0810.0126002000表3 配礦相關參數(shù)Table 3 Ore blending related parameters名稱目標品位（%）最小任務量（/t）最大破碎量（/t）MoWO3Cu受礦點10.1060

32、.1010.007550021000受礦點20.1010.0780.009700021000受礦點30.1050.0950.013450021000表4 配礦相關參數(shù)Table 4 Ore blending related parameters巖性各受礦點的期望巖性配比（%）受礦點1受礦點2受礦點3矽卡巖78.235.342.6透輝石長35.418.622.0長英角巖10.441.28.5表5 配礦相關參數(shù)Table 5 Ore blending related parameters出礦點出礦點到受礦點的距離（km）受礦點1受礦點2受礦點310.91.661.221.41.060.931.01

33、.31.441.30.81.351.21.451.361.160.91.471.150.91.581.11.11.591.20.90.9101.180.81.5111.21.21.1121.31.21131.41.31.2將上述表中數(shù)據(jù)輸入配礦模型中，MOEA/D-AU算法對其進行求解，得到30組配礦結果。為了提高礦石的綜合回收率，構建綜合選礦回收率預測模型，對30組配礦結果進行篩選，將不滿足各個受礦點綜合回收率指標的配礦結果舍去，從而獲得最貼合實際礦山情況的配礦結果。預測模型根據(jù)此礦山實際情況和選礦廠的選礦試驗，確定將礦石的氧化率和有害物質二氧化硫的濃度作為訓練模型的輸入參數(shù)，通過構建的隨機

34、森林模型進行訓練來得到綜合選礦回收率。最終從篩選后的配礦結果中隨機選取一個配礦計劃，如表6所示。表6 配礦計劃Table 6 Ore blending plan出礦點配礦量（/t）受礦點1受礦點2受礦點31811035320017053031075240080350808955662183135173991252081048009753091010017500110812957121397001301011750通過計算表6中配礦計劃可以發(fā)現(xiàn)，品位偏差和巖性配比偏差都趨近于0，生產成本也優(yōu)于礦山現(xiàn)狀，由此可知，本文所求配礦計劃能夠迅速找到貼合礦山實際情況的配礦結果。為了更加直觀全面的分析所得配

35、礦計劃的優(yōu)勢，從綜合選礦回收率預測模型篩選過后的多組配礦結果中隨機選擇10個配礦計劃，分析10個配礦計劃中各個受礦點所得品位情況。圖6反映了各個受礦點鉬、鎢、銅金屬的品位情況，因選礦廠接受配礦后的礦石品位處于期望品位上下0.05%范圍內，所以圖中描繪了各受礦點品位上下限、各個受礦點期望的礦石品位以及求解結果。從圖6中可以看出，鉬礦石在1號、2號及3號受礦點的計算品位和目標品位都相差不大；對于鎢礦石，在1號、2號受礦點的計算品位多數(shù)略高于期望的目標品位，3號受礦點的計算品位有些波動，但都處于允許的品位波動范圍；對于銅礦石，在1號、2號及3號受礦點品位波動都相對穩(wěn)定，在期望的目標品位上下，取得理想

36、的結果。圖6鉬、鎢、銅三種金屬計算品位情況Fig 6 Calculated grades of molybdenum, tungsten and copper5 結論（1）本文針對露天礦多金屬多目標配礦問題，以生產成本、品位偏差和礦石巖性配比偏差最小為目標函數(shù)構建了多金屬多目標配礦優(yōu)化模型，有效解決了多金屬多目標配礦中廣泛存在的配礦結果粗略等問題。（2）鑒于露天礦多金屬多目標配礦模型的復雜性和約束性，在MOEA/D的基礎上提出了MOEA/D-AU算法求解配礦計劃，并在標準測試集上與原算法MOEA/D進行了對比驗證，結果表明MOEA/D-AU算法在求解復雜多目標優(yōu)化問題上具有一定優(yōu)勢。（3）鑒于

37、對選礦環(huán)節(jié)的各種影響因素考慮不充分，本文根據(jù)礦山實際數(shù)據(jù)構建了融合氧化率及有害物質參數(shù)的綜合回收率隨機森林預測模型，通過預測模型對多目標進化算法求解的多組配礦結果進行篩選，從而得到一個更加符合礦山實際生產情況的配礦計劃。參考文獻（References）：1 王李管,宋華強,畢林,等. 基于目標規(guī)劃的露天礦多元素配礦優(yōu)化J.東北大學學報(自然科學版),2017,38(07):1031-1036.WANG Liguan, SONG Huaqiang, BI Lin, et al. Optimization of Open Pit Multielement Ore Blending Based on

38、 Goal ProgrammingJ. Journal of Northeastern University( Natural Science), 2017, 38(07):1031-1036.2 胡乃聯(lián),李勇,李國清,等.用粒子群算法優(yōu)化編制露天礦生產作業(yè)計劃J.北京科技大學學報,2013,35(04):537-543. HU Nailian, LI Yong, LI Guoqing, et al. Optimization of open-pit-mining operational planning by using a particle swarm algorithmJ. Journa

39、l of University of Science and Technology Beijing, 2013, 35(04):537-543.3 黃啟富,陳建宏.基于多輪PSO算法的中長期動態(tài)優(yōu)化配礦研究J.計算機應用研究,2010,27(10):3711-3717.HUANG Qifu, CHEN Jianhong. Research on middle and long-term optimization of mine ore blending based on multi round PSO algorithmJ. Application Research of Computers,

40、 2010, 27(10):3711-3717.4 柯麗華,何揚揚,葉義成,等.武鋼烏龍泉礦配礦方案研究J.礦冶工程,2016,36(06):22-25.KE Lihua, HE Yangyang, YE Yicheng, et al. Ore Blending Program for Wulongquan Limestone Open Pit MineJ. Mining and Metallurgical Engineering, 2016, 36(06):22-25.5 顧清華,劉海龍,盧才武,等. 基于采選流程下的露天礦多金屬多目標配礦優(yōu)化模型J.金屬礦山,2020(03):56-63.

41、GU Qinghua, LIU Hailong, LU Caiwu, et al. Multi-metal Multi-objective Ore Optimization Model for Open-pit Mine Based on Mining and Beneficiation Process J. Metal Mine, 2020(03):56-63.6 侯煜.線性規(guī)劃法在露天礦配礦領域的應用J.礦業(yè)工程,2006(01):22-24.HOU Yu. The Principle and method for the ore blending in open-pitJ. Mining

42、 Engineering, 2006(01):22-24.7 劉文博,許特,孫德峰.基于混勻配礦的露天礦山供應鏈建模與求解J.同濟大學學報(自然科學版),2017,45(04):609-616.LIU Wenbo, XU Te, SUN Defeng. Modeling and Solution of Open Pit Supply Chain Considering Ore BlendingJ. Journal of Tongji University(Natural Science), 2017,45(04):609-616.8 Eivazy H, Askarinasab H. A mix

43、ed integer linear programming model for short-term open pit mine production schedulingJ. Mining Technology, 2012, 121(2):97-108.9 井石滾,盧才武,李發(fā)本,等.基于GIS/GPS/GPRS的露天礦配礦動態(tài)管理系統(tǒng)J.金屬礦山,2009(06):140-144.JING Shigun, LU Caiwu, LI Faben, et al. Dynamic Management System of the Ore Blending in Open-pit Mine Based on GIS/GPS/GPRSJ. Metal Mine, 2009(06):140-144.10 吳麗春,王李管,彭平安,等.露天礦配礦優(yōu)化方法研究J.礦冶工程,2012,32(04):8-12.WU Chunli