范鳳巖,柳群義
(1.中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所,北京 100037;2.中國地質(zhì)科學(xué)院全球礦產(chǎn)資源戰(zhàn)略研究中心,北京 100037;3.中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所自然資源部鹽湖資源與環(huán)境重點實驗室,北京 100037)
錫在日常生活中是一種非常重要的有色金屬原料[1],也是我國發(fā)展戰(zhàn)略性新興產(chǎn)業(yè)的重要礦產(chǎn)資源,應(yīng)用領(lǐng)域幾乎涵蓋我國七大戰(zhàn)略性新興產(chǎn)業(yè)[2]。從全球范圍來看,我國錫礦資源最豐富,但隨著國內(nèi)經(jīng)濟發(fā)展,錫資源需求量也日益增加,導(dǎo)致錫礦的基礎(chǔ)儲量持續(xù)減少,需要借助進口錫礦沙及精煉錫來補充國內(nèi)需求的不足。2017年,中國錫消費量18.34萬t,國內(nèi)礦山錫產(chǎn)量9.3萬t,二次錫回收量2.9萬t,凈進口錫資源6.14萬t(折金屬量),對外依存度達33.5%[3],錫資源的優(yōu)勢正逐漸消失。因此,開展錫資源的供應(yīng)安全動態(tài)評價對于保證錫資源穩(wěn)定供應(yīng),提高錫資源的保障程度具有重要意義。
礦產(chǎn)資源供應(yīng)安全是指一個國家或地區(qū)為了實現(xiàn)保障其經(jīng)濟與社會持續(xù)穩(wěn)定健康發(fā)展的需要,及時、穩(wěn)定、經(jīng)濟和持續(xù)地獲取所需礦產(chǎn)資源,并保證其發(fā)展所依賴的自然基礎(chǔ)和生態(tài)環(huán)境處于可持續(xù)發(fā)展的狀態(tài)[4]。開展礦產(chǎn)資源安全評價作為國家安全戰(zhàn)略的重要一環(huán),受到了學(xué)者們的廣泛關(guān)注。L?SCHEL等[5]將能源安全分為可用性、可承受性、技術(shù)發(fā)展、可持續(xù)性和規(guī)制5個維度20個部分,為分析國家能源安全政策和績效提供了一個綜合的、可行的框架;WANG等[6]采用熵值法計算出安全綜合指數(shù)對中國銅資源的安全態(tài)勢進行評價;ZHOU等[7]采用蒙特卡洛模擬對清潔能源技術(shù)中采用的12種礦產(chǎn)資源的中短期供應(yīng)風(fēng)險進行了定量評估;竇迅等[8]采用模糊綜合評價法建立了電煤供應(yīng)安全預(yù)警機制。朱春華[9]、劉賢信[10]分別采用多元聚合閾值評價法和主成分分析法對我國鐵礦資源供應(yīng)風(fēng)險進行了評價;黃曉兵[11]基于層次分析法,建立中短期安全評價指標(biāo)體系對我國2000~2020年鈷資源安全進行定量分析;林志峰[12]、張佳東等[13]同樣利用層次分析法構(gòu)建對我國鎳資源的安全形勢進行了評價。類似的研究還涉及銦礦[14]、鉻礦[15]、石油[16-17]、天然氣[18]等多個礦產(chǎn)資源。
可以看出,國內(nèi)外學(xué)者的研究對象主要集中在能源礦產(chǎn)以及鐵、銅等部分固體礦產(chǎn)。錫礦一直以來被認(rèn)為是我國的優(yōu)勢礦產(chǎn),其安全問題鮮少受到學(xué)者們的關(guān)注。從研究手段上看,現(xiàn)有的評價方法主要分為兩大類:一類是主觀賦權(quán)法,主要包括層次分析法、模糊綜合評價法等;一類是客觀賦權(quán)法,主要包括熵權(quán)法、主成分分析法等。主觀賦權(quán)法常借助專家經(jīng)驗,很難排除個人因素對于指標(biāo)權(quán)重的影響。其求解過程較為粗糙,對于精度較高問題的解決方案有可能與實際情況并不相符[19]??陀^賦權(quán)法依據(jù)方案的真實數(shù)據(jù)計算,結(jié)果相對客觀,但其權(quán)重的確定完全依賴于客觀數(shù)據(jù)之間的關(guān)系,可能與實際情況相差較大?;谏鲜鲈颍疚膶⒆畛R姷目陀^賦權(quán)法—熵權(quán)法和最常見的主觀賦權(quán)法—層次分析法相結(jié)合,借助郭金維等[20]提出的改進的指標(biāo)權(quán)重綜合求取方法,對我國錫資源的供應(yīng)安全態(tài)勢及演化軌跡進行評價,為完善安全防范機制和政策提供依據(jù)。
錫資源供應(yīng)是否安全受眾多因素的影響,包括國內(nèi)資源的稟賦條件、供需狀況、二次資源的利用與替代情況、抵御國際風(fēng)險的能力等。評價指標(biāo)的選取應(yīng)在遵循科學(xué)性、綜合性、可獲取性、可持續(xù)性、動態(tài)性原則的基礎(chǔ)上,與錫資源的實際情況相吻合。本文針對錫資源開采、開發(fā)、利用特點和面臨的問題,借鑒國內(nèi)外已有的研究成果,確定資源因素、供需因素、風(fēng)險因素、可持續(xù)性因素四類十六個指標(biāo)對錫資源的供應(yīng)安全進行定量評價(圖1)。
1) 品位:指礦石中有用組分的含量,它直接影響選礦效率。作為衡量礦石質(zhì)量好壞的最主要指標(biāo),資源品位在一定程度上反映出一國資源稟賦情況。品位越高,資源稟賦越好,資源供應(yīng)越安全,反之,越不利于資源供應(yīng)。
2) 儲采比:指以當(dāng)前的產(chǎn)量來看,國內(nèi)剩余儲量可供消費的時間長短。儲采比越大,表明國內(nèi)剩余資源越豐富,資源的靜態(tài)保障年限越長,資源供應(yīng)安全程度越高,反之越低;其計算公式為:儲采比=當(dāng)年錫礦基礎(chǔ)儲量/當(dāng)年錫礦產(chǎn)量。
圖1 錫資源供應(yīng)安全評價指標(biāo)體系Fig.1 Index system of tin resources supply safety evaluation
3) 勘查潛力:礦產(chǎn)資源勘查資金的投入能帶來多少新增資源儲量,用當(dāng)年新增查明基礎(chǔ)儲量與勘查資金投入的比值表示,比值越大,勘查資金的投入會帶來更多的新增資源儲量,對資源消耗的補充程度也越高,越有利于資源供應(yīng)安全;其計算公式為:勘查潛力=新增查明資源儲量/錫礦勘查資金投入。
4) 人均資源量:平均每萬人錫礦資源占有量,人均資源量是一個相對指標(biāo),值越大,說明人均資源占有量越多,資源供應(yīng)越安全,反之越不安全;其計算公式為:人均資源量=錫礦基礎(chǔ)儲量/人口。
5) 資源保障度:反映資源生產(chǎn)能在多大程度上滿足消費需求。比值越大說明國內(nèi)資源供應(yīng)能力越強,資源供應(yīng)安全程度越高,反之,資源供應(yīng)安全程度越低;其計算公式為:資源保障度=精煉錫產(chǎn)量/精煉錫消費量。
6) 國際生產(chǎn)集中度:反映全球錫資源生產(chǎn)分布的集中程度。集中度越高,說明錫礦生產(chǎn)集中于少數(shù)國家,容易形成一定程度的市場壟斷,使資源供應(yīng)安全程度降低;其計算公式為:國際生產(chǎn)集中度=前5位錫礦生產(chǎn)國產(chǎn)量/全球錫礦產(chǎn)量總和。
7) 國內(nèi)生產(chǎn)能力:如果一個國家的國內(nèi)生產(chǎn)能力越強,即在全球礦產(chǎn)品的生產(chǎn)中占據(jù)較大的份額,資源供應(yīng)的能力就越強,不易受到國際市場任何一次小波動的影響;其計算公式為:國內(nèi)生產(chǎn)能力=國內(nèi)錫礦產(chǎn)量/全球錫礦產(chǎn)量。
8) 價格波動:價格波動是反映國際市場供需狀況的主要因素之一;價格波動幅度越大,說明國際市場上該資源的供需差距變化越大;價格波動可以用價格波動指數(shù)來表示。為了更真實地反應(yīng)價格波動情況,價格以2000年為基期,采用進出口價格指數(shù)進行折算,其計算公式為:價格波動Pc=∣(Pn-P0)∣/P0,式中:Pn為本期資源價格;P0為上期價格。
9) 進口集中度:主要反映進口來源國集中程度;進口集中度越高,表明進口來源和渠道太集中,不利于風(fēng)險分散,相應(yīng)地,資源供應(yīng)安全程度也越低;其計算公式為:進口集中度=前5位國家進口量的和/全球總進口量。
10) 對外依存度:對外依存度是反映本國礦產(chǎn)資源消費量對國外資源依賴程度的指標(biāo),對外依存度越高,對國際市場依存度越高,抵御國際市場風(fēng)險的能力越低,資源供應(yīng)安全程度越低;反之,資源供應(yīng)安全程度越高。其計算公式為:對外依存度=精煉錫凈進口量/本國年消費量。
11) 進口來源國穩(wěn)定度:該指標(biāo)實際上是用來描述本年度主要進口量中有多少是來源于上一年度的主要進口國。若進口來源國經(jīng)常更換或來自主要進口國的進口量變動較大,則會增加供應(yīng)風(fēng)險。其計算公式為:進口來源國穩(wěn)定度=前5國持續(xù)進口量總和/前5國總進口量。
12) 企業(yè)競爭力:用錫礦企業(yè)的利潤率來表示。企業(yè)的盈利能力越強,表示企業(yè)在國際市場上的競爭能力越強,越能抵御不確定因素帶來的潛在風(fēng)險,越有利于資源的持續(xù)供應(yīng),因而安全程度也越高。其計算公式為:企業(yè)競爭力=營業(yè)利潤/主營業(yè)務(wù)收入。
13) 錫冶煉綜合能耗:由于冶煉環(huán)節(jié)的能源消耗遠高于采選環(huán)節(jié),因此采用冶煉環(huán)節(jié)的綜合能耗反映錫行業(yè)的環(huán)境友好程度。
14) 錫冶煉綜合回收率:反映出資源的開發(fā)利用水平,錫冶煉綜合回收率越高,說明資源的使用越有效率,資源的浪費越少,越有利于資源的安全供應(yīng)。
15) 二次資源供應(yīng)能力:反映了資源重復(fù)利用的程度。資源的重復(fù)利用可以減少不可再生資源的消費。數(shù)值越高,資源的原生供給壓力越小。其計算公式為:二次資源供應(yīng)能力=再生錫產(chǎn)量/精煉錫產(chǎn)量。
16) 資源消費強度:該指標(biāo)反映出錫礦消費對GDP的貢獻度。單位GDP錫礦消耗越高,說明對錫礦的需求也越大,供給壓力越大,越不安全。其計算公式為:資源消費強度=錫礦產(chǎn)量/GDP。
本文所采用的數(shù)據(jù)主要涉及資源、經(jīng)濟、環(huán)境、貿(mào)易等多個方面。品位、礦山錫產(chǎn)量、企業(yè)的利潤、主營業(yè)務(wù)收入(1)《有色金屬工業(yè)統(tǒng)計年鑒》缺少2000年和2001年利潤和收入數(shù)據(jù),用有色金屬行業(yè)利潤率代替。、錫冶煉的綜合能耗及回收率數(shù)據(jù)來源于《中國有色金屬工業(yè)統(tǒng)計年鑒》。資源數(shù)據(jù)中的基礎(chǔ)儲量數(shù)據(jù)來源于《全國礦產(chǎn)資源儲量通報》和中國統(tǒng)計局,新增查明資源儲量數(shù)據(jù)來源于《中國礦產(chǎn)資源年報》和《中國礦產(chǎn)資源報告》,精煉錫產(chǎn)量和消費量數(shù)據(jù)來源于世界金屬統(tǒng)計局(WBMS)。經(jīng)濟數(shù)據(jù)中的人口、GDP、進出口商品價格指數(shù)數(shù)據(jù)來源于中國統(tǒng)計局,錫礦勘查資金投入數(shù)據(jù)來源于《中國國土資源統(tǒng)計年鑒》《中國礦業(yè)年鑒》和《中國礦產(chǎn)資源報告》,價格數(shù)據(jù)來源于SNL。貿(mào)易數(shù)據(jù)中的精煉錫進出口數(shù)據(jù)來源于聯(lián)合國貿(mào)易數(shù)據(jù)庫(UN Comtrade Database)。再生錫產(chǎn)量數(shù)據(jù)來源于國際錫業(yè)協(xié)會(ITRI)。
熵權(quán)法(entropy method)是根據(jù)各指標(biāo)所提供的信息量的大小來確定指標(biāo)權(quán)重的一種方法。若某項指標(biāo)的指標(biāo)值變異程度越大,則信息熵越小,表明其信息的有效價值越大,在綜合評價中對目標(biāo)的影響也就越大,權(quán)重相應(yīng)得也大。反之,某項指標(biāo)值差異程度越小,則信息熵越大,指標(biāo)的權(quán)重也越小。熵權(quán)法計算結(jié)果完全依賴于真實數(shù)據(jù)本身,因其計算簡單,廣泛應(yīng)用于解決多目標(biāo)決策問題。計算過程如下所述。
1) 設(shè)有m個評價對象,n個評價指標(biāo),則形成評價系統(tǒng)的初始數(shù)據(jù)矩陣,見式(1)。
(1)
式中:xij(i=1,2,…,m;j=1,2,…,n)為第i個評價對象在第j項指標(biāo)中的數(shù)值;Xj(j=1,2,…,n)為第j個指標(biāo)的全部評價對象的列向量數(shù)據(jù)。
由于各指標(biāo)的量綱單位均存在差異,為消除因量綱不同對評價結(jié)果造成的影響,需要對各指標(biāo)進行無量綱化處理。常采用極差變換法,公式見式(2)和式(3)。
(適用于正向指標(biāo))
(2)
(適用于負向指標(biāo))
(3)
i=1,2,…,m;j=1,2,…,n
(4)
3) 計算第j項指標(biāo)信息熵的值ej,見式(5)。
(5)
式中,k=1/lnm為非負常數(shù),且0≤ej≤1。
4) 計算第j項指標(biāo)的權(quán)重wj,見式(6)
j=1,2,…,n
(6)
層次分析法(AHP)將復(fù)雜的問題分解為若干個因素,并將這些因素按支配關(guān)系分組,形成一個有序的遞階層次結(jié)構(gòu),然后對同一層次的因素構(gòu)建判斷矩陣,通過兩兩比較的方式確定層次中各因素的相對重要性,最后綜合判斷以確定決策中各因素相對重要性。
其求取指標(biāo)權(quán)重的過程如下所述。
1) 建立問題的層次結(jié)構(gòu)模型。通過對系統(tǒng)的分析,確定系統(tǒng)的總目標(biāo),收集影響目標(biāo)決策的準(zhǔn)則、策略等。確定各個層次的影響因素,上一層次的元素作為準(zhǔn)則對下一層次的元素起支配作用。
2) 構(gòu)建兩兩比較判斷矩陣。通常采用九級標(biāo)度法構(gòu)造判斷矩陣A=(aij)n×n,確定同層子元素相對于上一層元素的相對重要性排序。
3) 由判斷矩陣計算被比較元素的權(quán)重,見式(7)和式(8)。
(7)
(8)
式中:pi為求出的Ai指標(biāo)的權(quán)重;Mi為判斷矩陣同一行各元素的乘積。
4) 檢驗判斷矩陣的一致性,見式(9)~(11)。
CR=CI/RI
(9)
CI=(λmax-n)/(n-1)
(10)
(11)
式中:CR為隨機一致性比率;CI為一致性指標(biāo);RI為平均隨機一致性指標(biāo),其取值見文獻[21];λmax為判斷矩陣的最大特征根,其中(Ap)i為向量Ap的第i個元素,p=(p1,p2,…,pn)T,即由各指標(biāo)權(quán)重組成的列向量。
只有當(dāng)CR<0.10時,認(rèn)為判斷矩陣具有良好。
層次分析法和熵權(quán)法都可能出現(xiàn)指標(biāo)權(quán)重與實際情況偏差較大,部分學(xué)者[22-23]嘗試將兩種方法的權(quán)重相結(jié)合,得到綜合權(quán)重ωi,通常采用式(12)。
(12)
這種方法在求綜合權(quán)重時,只是將兩種方法求得的最底層指標(biāo)權(quán)重的結(jié)果簡單綜合,并沒有將兩種方法求取權(quán)重的中間過程有機地融合起來。改進綜合權(quán)重的求解方法,將兩種方法的中間過程相結(jié)合,得到的指標(biāo)更加客觀準(zhǔn)確,具體步驟如下所述。
1) 設(shè)有m個上層準(zhǔn)則,n個子準(zhǔn)則,通過層次分析法判斷矩陣求得上層準(zhǔn)則權(quán)重B={β1,β2,…,βm},各子準(zhǔn)則的權(quán)重Φ={Φ1,Φ2,…,Φn};熵權(quán)法求得各準(zhǔn)則的權(quán)重A={α1,α2,…,αn}。
最后,利用組合權(quán)重與第二層各因素指標(biāo)中處理后的標(biāo)準(zhǔn)化數(shù)據(jù)相乘再求和,即得到最終各年度的評價值。
本研究依據(jù)實際情況,參考相關(guān)文獻,按線性刻度法來確定錫礦安全評價指標(biāo)分級標(biāo)準(zhǔn)(表1)。
表1 錫資源供應(yīng)安全判斷標(biāo)準(zhǔn)Table 1 Criteria for the supply safety of tin resources
通過上述錫資源安全度評價方法和步驟,可以求出個指標(biāo)權(quán)重(表2)以及資源安全、環(huán)境安全、產(chǎn)業(yè)安全和市場安全度的分值表(表3)和變化趨勢(圖2)。
表2 改進的層次熵權(quán)法各指標(biāo)權(quán)重表Table 2 Index weight of improved hierarchical entropy weight method
表3 各年度評價值表Table 3 Annual index evaluation table
錫資源安全狀況依賴于資源因素、供需因素、風(fēng)險因素和可持續(xù)因素的相互作用。由圖2可以看出,除可持續(xù)因素呈上升趨勢,對錫資源的安全起積極作用外,資源因素、供需因素、風(fēng)險因素都呈下降趨勢,對錫資源安全產(chǎn)生負面影響。在四個一級指標(biāo)的共同影響和作用下,中國錫資源安全狀況整體呈惡化趨勢。具體可分為兩個階段:第一階段(2000~2007年),錫資源安全綜合指數(shù)呈波動下降趨勢,安全等級由安全到危險;第二階段(2007~2017年),系資源安全綜合指數(shù)在波動中趨于平穩(wěn),安全指數(shù)在0.2~0.4之間小幅波動,處于危險狀態(tài)。說明錫資源安全水平?jīng)]有得到改善,處于不安全的狀態(tài)下。
中國錫資源因素評價值在研究期內(nèi)經(jīng)歷了從波動下降到穩(wěn)定下降的趨勢,這與勘查潛力的變動趨勢一致。由于新礦產(chǎn)地的發(fā)現(xiàn)具有偶發(fā)性,勘查資金的投入未必能夠取得相應(yīng)的基礎(chǔ)儲量的增加,因此,受權(quán)重最大的勘查潛力的影響,2006年之前,資源因素評價值波動較大;2006年之后,資源安全呈穩(wěn)定下降的趨勢。在整個研究期,品位、儲采比、人均資源占有量都呈現(xiàn)出不斷下降的趨勢,對資源因素評價值產(chǎn)生了負面影響,這也從一定程度上反映出近十幾年來,我國錫礦資源的補給趕不上資源消耗的速度,所以資源因素評價值不斷下降。
供需因素的評價值在2009年之前震動下降,2009年受金融危機影響,供需因素評價值在該年度下降明顯,2009年之后,隨著中國經(jīng)濟發(fā)展進入新常態(tài),工業(yè)化步入后工業(yè)化發(fā)展階段,供需因素評價值開始震動上升。分析供需因素下的二級指標(biāo)變動趨勢可以發(fā)現(xiàn),國際生產(chǎn)集中度和國內(nèi)生產(chǎn)能力的變動比較平穩(wěn),資源保障程度呈緩慢下降的趨勢,價格波動在2009年之前波動上升,2009年之后波動下降,因此,資源保障程度的趨勢決定了供需因素評價值總體呈下降趨勢。價格波動的拐點決定了供需因素評價值在2009年出現(xiàn)轉(zhuǎn)折,這也說明,供需因素更易受到市場變化的影響(圖3)。
圖2 2000~2017年各指標(biāo)因素評價值和總評價值走勢圖Fig.2 Chart of each indicator evaluation value and total evaluation value from 2000 to 2017
圖3 2000~2017年供需因素各指標(biāo)變動趨勢Fig.3 Trends in various indicators of supply and demand factors from 2000 to 2017
風(fēng)險因素的評價值總體上呈下降趨勢,這主要是受權(quán)重較大的進口集中度的影響。研究期內(nèi),進口集中度基本上呈逐年上升的態(tài)勢,說明資源的進口越來越集中于少數(shù)國家,這不利于風(fēng)險的分散,導(dǎo)致風(fēng)險因素的評價值下降,雖然企業(yè)競爭力對風(fēng)險因素有著積極影響,但也無法彌補進口集中度對風(fēng)險因素評價值產(chǎn)生的負面影響。進口來源國穩(wěn)定度除個別年份變動較大外,整體趨勢比較穩(wěn)定。對外依存度(2)由于對外依存度計算公式有差異,因此造成本文的計算結(jié)果與引言中的引用有所出入。在2007年之后從無到有,也從側(cè)面反映出國內(nèi)錫資源無法滿足本國的需求,我國錫礦的資源優(yōu)勢不在(圖4)。
可持續(xù)因素評價值可以分為三個階段,2005年之前緩慢上升,2005~2010年迅速上升,2010年之后緩慢下降,2005年之后可持續(xù)因素評價值快速上升的主要原因是錫冶煉的綜合能耗和資源消費強度在2005年之后顯著下降。受再生錫產(chǎn)量在2010年達到峰值后下降的影響,可持續(xù)因素評價值也隨之下降。這主要是由于緬甸原生錫礦近幾年低價且容易獲得、電子焊料中錫的使用量下降、制造廠新廢料供應(yīng)減少、冶煉廠來自于越來越嚴(yán)格的環(huán)保壓力導(dǎo)致。錫冶煉回收率由于變化不大,對可持續(xù)因素的影響不大(圖5)。
1) 錫資源供應(yīng)安全狀況在資源因素、供需因素、風(fēng)險因素和可持續(xù)因素的相互作用下整體呈惡化趨勢,安全等級由安全到危險。除可持續(xù)因素呈上升趨勢,對錫資源的安全起積極作用外,資源因素、供需因素、風(fēng)險因素對錫資源安全產(chǎn)生負面影響。說明錫資源供應(yīng)安全狀況堪憂。
圖4 2000~2017年風(fēng)險因素各指標(biāo)變動趨勢Fig.4 Trends in various indicators of risk factors from 2000 to 2017
圖5 2000~2017年可持續(xù)因素各指標(biāo)變動趨勢Fig.5 Trends in various indicators of sustainability factors from 2000 to 2017
2) 錫資源因素評價值在研究期內(nèi)經(jīng)歷了從波動下降到穩(wěn)定下降的趨勢,這與勘查潛力的變動趨勢一致。在整個研究期,品位、儲采比、人均資源占有量都呈現(xiàn)出不斷下降的趨勢,對資源因素評價值產(chǎn)生了負面影響。受市場變化的影響,供需因素的評價值在2009年之前震動下降,2009年受金融危機影響,隨著中國經(jīng)濟發(fā)展進入新常態(tài),供需因素評價值上升。風(fēng)險因素的評價值總體上呈下降趨勢。雖然企業(yè)競爭力對風(fēng)險因素有著積極影響,但也無法彌補權(quán)重較大的進口集中度對風(fēng)險因素評價值產(chǎn)生的負面影響。
3) 可持續(xù)因素經(jīng)歷了由緩慢上升到快速上升再到緩慢下降的變動趨勢。受再生錫產(chǎn)量在2010年達到峰值后下降的影響,可持續(xù)因素評價值在2010年之后開始下降。