賈祥英，邢萬里，苗媛媛

(1.中國地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083； 2.中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所，北京 100037； 3.中國地質(zhì)科學(xué)院礦產(chǎn)資源研究所自然資源部成礦作用與資源評價重點實驗室，北京 100037)

銅礦是戰(zhàn)略性礦產(chǎn)資源，其產(chǎn)業(yè)涵蓋電力、家電、機械、冶金、交通、輕工、建筑以及國防等諸多領(lǐng)域，對國家經(jīng)濟安全、國防安全以及戰(zhàn)略性新興產(chǎn)業(yè)的發(fā)展至關(guān)重要，是我國僅次于鐵和鋁的重要金屬礦產(chǎn)資源。根據(jù)UN Comtrade、USGS以及自然資源部數(shù)據(jù)統(tǒng)計，近年來中國銅礦進口量居于世界首位，對外依存度高達60%以上；中國銅礦經(jīng)濟可采儲量位居全球前十位，但是我國銅礦存在“貧礦多，富礦少”等特點，致使我國銅礦雖儲量基數(shù)較大，但經(jīng)濟可采儲量卻較少，且2012～2017年期間，中國是唯一銅礦資源經(jīng)濟可采儲量出現(xiàn)負增長的國家；另外，我國精銅產(chǎn)量與消費量均居于世界首位。2017年以來，為響應(yīng)“環(huán)保污染治理”以及“綠水青山”的政策要求，我國減少對于礦產(chǎn)勘查的投入，固廢進口政策收緊等原因，致使我國銅礦進口要求越加嚴格。因此，為解決上述我國銅礦資源保障程度較低，對外依存度較大等問題，需要探尋我國銅礦重要進口來源地及貿(mào)易過程中的關(guān)鍵性影響因素，對制定銅礦進出口政策以及緩解我國銅礦進口壓力具有十分重要的意義。

國內(nèi)外學(xué)者在基于模型的礦產(chǎn)資源貿(mào)易研究方面取得諸多成果。王安建等[1]以人均礦產(chǎn)資源消費與人均GDP S形規(guī)律、礦產(chǎn)資源消費強度變化倒U形規(guī)律和礦產(chǎn)資源需求波次遞進規(guī)律等為基礎(chǔ)，建立能源與礦產(chǎn)資源需求綜合預(yù)測模型，系統(tǒng)預(yù)測了2010～2030年全球及中國一次能源、粗鋼、銅、鋁需求趨勢。柳群義等[2-3]基于S模型對2040年不同類型國家以及全球銅的需求量進行預(yù)測。與需求預(yù)測的其他模型方法(如線性外推法、類比分析法、部門預(yù)測法、數(shù)學(xué)模型法)相比較，S形模型法具有一定的可靠性。賈立文等[4]建立適用于中國和俄羅斯等人均鐵礦石需求增長型國家的面板模型，與灰色模型、協(xié)整模型、ARIMA模型的預(yù)測結(jié)果進行對比，認為四類模型分析鐵礦石需求問題的能力排序為：面板模型>協(xié)整模型>灰色模型>ARIMA模型。上述多為基于模型的對礦產(chǎn)資源供需方面的研究，近年來，許多學(xué)者對基于模型的礦產(chǎn)資源國際貿(mào)易的動態(tài)演變進行了研究[5-9]。

基于引力模型的研究中，國內(nèi)外已在國際商品貿(mào)易的許多方面取得研究成果。其中，朱永光等[10]對傳統(tǒng)引力模型進行改進，建立空間交互模型，對鐵礦石的國際貿(mào)易影響因素進行分析，通過貿(mào)易空間交互過程的分析，提出中國鐵礦石進口策略。黃宇等[11]借鑒庫侖引力的思想，引入地緣經(jīng)濟合作強度的概念和測度方法，基于中國、美國、日本三國2007～2016年各類貿(mào)易數(shù)據(jù)，探討國家間地緣經(jīng)濟競爭與合作的規(guī)律，在考慮距離因素的基礎(chǔ)上，分產(chǎn)品、進口和出口，準確地刻畫了國家間經(jīng)濟合作或競爭的強弱程度。

前人的研究成果為基于模型的礦產(chǎn)資源貿(mào)易研究奠定了良好的基礎(chǔ)，但是關(guān)于中國銅礦進口的關(guān)鍵性影響因素的分析較少。本文在前人研究的基礎(chǔ)上，基于2012～2017年相關(guān)數(shù)據(jù)，通過構(gòu)建引力模型，對影響我國銅礦進口貿(mào)易的因素進行分析，得出的結(jié)論和對策建議可為決策者制定銅礦進出口政策提供參考。

1 研究方法與數(shù)據(jù)來源

1.1 數(shù)據(jù)選取

中國銅礦主要進口來源國為智利、秘魯、澳大利亞、墨西哥、蒙古國，本文將以上六國作為新型貿(mào)易引力模型的研究對象，選取2012～2017年中國、智利、秘魯、澳大利亞、墨西哥、蒙古國的銅礦貿(mào)易數(shù)據(jù)(表1)進行處理分析。王安建等[1]、賈立文等[4]、朱永光等[10]在礦產(chǎn)資源貿(mào)易研究中，認為經(jīng)濟規(guī)模以及人口結(jié)構(gòu)對礦產(chǎn)資源消費產(chǎn)生顯著影響。本文引入市場結(jié)構(gòu)(pop)、經(jīng)濟距離(gdp)等指標構(gòu)建改進貿(mào)易引力模型。將中國從其他五國的銅礦進口量作為解釋變量，市場結(jié)構(gòu)(pop)、經(jīng)濟距離(gdp)、地理距離、出口國銅礦儲量、產(chǎn)量、人均消費量等作為虛擬變量。另外，選取“是否與中國簽署自由貿(mào)易協(xié)議”為政策因素。簽署自由貿(mào)易協(xié)議(以下簡稱“FTA”)的貿(mào)易雙方，為進出口商品關(guān)稅等多方面提供便利。表2為本文所選取基于引力模型的變量的解釋以及數(shù)據(jù)來源。

1.2 研究方法

1.2.1 傳統(tǒng)引力模型

根據(jù)Tinbergen對于引力模型的研究，考慮貿(mào)易雙方經(jīng)濟規(guī)模和地理距離對貿(mào)易的影響，傳統(tǒng)引力模型為式(1)。

表1 2012～2017年中國主要進口來源國銅礦貿(mào)易相關(guān)數(shù)據(jù)Table 1 Data on copper trade of Chinese major importing countries from 2012 to 2017

表2 引力模型變量解釋Table 2 Explanation of gravity model variables

(1)

相應(yīng)的回歸方程為式(2)。

lnYij=β0+β1ln(Mi)+β2ln (Mj)-

β3ln(Dij)+εij

(2)

式(1)和式(2)中：Yij為i從j國進口的貿(mào)易量；G為常量；Mi、Mj為參與貿(mào)易的國家的GDP；D為貿(mào)易兩國之間的距離；?ij為誤差項；β0、β1、β2、β3為相關(guān)系數(shù)。

1.1.2 改進引力模型的構(gòu)建

ANDERSON等[12]指出Tinbergen建立的引力模型只考慮了貿(mào)易雙方的經(jīng)濟規(guī)模和地理距離對貿(mào)易流量的影響，還應(yīng)考慮市場效應(yīng)和地方保護等因素;朱永光等[10]認為引力模型的作用機制包括推力作用、拉力作用以及阻力因素，引入資源儲量、經(jīng)濟距離、產(chǎn)業(yè)結(jié)構(gòu)等指標，構(gòu)建引力模型;石超等[13]設(shè)立進口和出口貿(mào)易額兩個被解釋變量，引入GDP、距離、“是否沿海”“是否臨界”等解釋變量，建立引力模型。在此次計量過程中，筆者引入市場結(jié)構(gòu)、經(jīng)濟距離、地理距離、出口國經(jīng)濟可采儲量、產(chǎn)量、消費量、自由貿(mào)易協(xié)議政策等變量，旨在拓展引力模型的理論范圍，探究影響我國銅礦進口的關(guān)鍵性因素。

引入所列變量后引力模型為式(3)。

lny=β0+β1ln(pop)+β2ln(gdp)+β3ln(r)+

β4ln(c)+β5ln(d)+β6ln(x)+FTA+εij

(3)

2 改進引力模型回歸方程分析

2.1 回歸結(jié)果分析

將上述數(shù)據(jù)構(gòu)建變量數(shù)據(jù)庫，利用Stata工具，得出回歸方程。Stata回歸結(jié)果中(表3)，F(xiàn)檢驗的P值為0，模型總體是統(tǒng)計顯著的，模型較好；R2大于80%，說明模型的擬合度很高，模型較好，說明本文建立的引力模型適用于中國銅礦進口量的研究。

得到回歸方程見式(4)。

lny=34.78404+1.028357ln(pop)-

0.0428637ln(gdp)-1.087128ln(r)+

1.636419ln(c)-0.6334249ln(d)-

0.2579386ln(x)+0.713941FTA

(4)

由式(4)可知，當市場結(jié)構(gòu)(pop)增加時，中國進口量增加人口市場比例增量的102.84%；當經(jīng)濟距離(gdp)增加時，中國銅礦石進口量減少經(jīng)濟規(guī)模增量的4.29%；當出口國銅礦儲量(r)增加時，銅礦進口量減少儲量增值的108.71%；當出口國銅礦產(chǎn)量(c)增加時，中國銅礦進口量增加產(chǎn)量增值的163.64%；當?shù)乩砭嚯x(d)增加時，中國銅礦石進口量減少距離增量的63.34%；當出口國人均精銅消費量(x)增加時，中國銅礦進口量減少消費量增值的25.79%；與中國簽署自由貿(mào)易協(xié)議(FTA)政策為1時，進口量增加FTA從0至1的增量的71.39%。

2.2 實證結(jié)果分析

根據(jù)回歸分析結(jié)果可知，解釋變量pop、c、d、x和FTA在動態(tài)過程中的符號與預(yù)期的符號一致，而解釋變量gdp和r在動態(tài)過程中的符號與預(yù)期的不一致。

1) 解釋變量pop符號為正，與預(yù)期符號一致，說明由人口之比反映的市場差距與中國進口量成正比。但是p值為0.076，大于0.05，表示人口市場差距與中國銅礦進口的相關(guān)性不大。中國為世界人口第一大國，人口之比越大，中國與貿(mào)易另一方國家的人口市場差距越小，消費層次差別越小，貿(mào)易阻力減小，進口量增加。

表3 Stata線性回歸結(jié)果Table 3 Stata linear regression results

2) 解釋變量c符號為正，與預(yù)期符號一致，說明出口國人均銅礦產(chǎn)量與中國進口量成正相關(guān)，出口國人均銅礦產(chǎn)量越高，出口國富余銅礦儲量越多，滿足自身需要的同時，可出口銅礦越多，符合引力模型對出口國人均銅礦產(chǎn)量的理論解釋。同時，回歸分析結(jié)果中，p值為0.000，表明出口國人均銅礦產(chǎn)量對中國進口量的影響十分顯著。

3) 解釋變量d的符號為負，說明兩國之間的地理距離與雙方貿(mào)易量成負相關(guān)。朱永光等[10]指出地理距離與礦產(chǎn)品貿(mào)易存在正U型關(guān)系，起始地理距離越大，貿(mào)易雙方文化差距越大，貿(mào)易運輸成本增加，對雙方貿(mào)易存在阻礙作用；到達一定距離后，距離對礦產(chǎn)品貿(mào)易的影響不再顯著為負?；貧w分析結(jié)果中，地理距離的p值為0.192，大于0.05，且相關(guān)性符號為負，表示地理距離對中國銅礦進口的影響并不顯著，出現(xiàn)由負相關(guān)轉(zhuǎn)為正相關(guān)的可能性，地理距離對中國銅礦石進口的影響不再顯著。

4) 解釋變量x符號為負，與預(yù)期符號一致，說明出口國人均精銅消費量與中國進口量成負相關(guān)，出口國人均精銅消費量越高，本國的銅礦需求越多，可出口銅礦量越少。同時，回歸分析結(jié)果中p值為0.026，小于0.05，表明出口國人均精銅消費量對中國進口量的影響十分顯著。

5) 解釋變量FTA符號為正，與預(yù)期符號一致，說明政策提供的便利促進中國銅礦的進口量。但是p值為0.101，大于0.05，表示政策與中國銅礦進口的相關(guān)性不顯著。筆者認為，多項優(yōu)惠政策的疊加，會促進兩國貿(mào)易，單一的政策促進作用不明顯。如智利與中國同為WTO成員，兩國2005簽署自由貿(mào)易協(xié)議，同時，2015年中國從智利進口金屬礦產(chǎn)品實施零關(guān)稅政策，多項優(yōu)惠政策的疊加，使得中國從智利的銅礦進口量始終排在中國進口量的首位。

6) 解釋變量gdp符號為負，與預(yù)期符號不一致，且p值為0.506，大于0.05，說明市場差距與中國進口量相關(guān)性不顯著。王慶喜等[14]引入經(jīng)濟距離指標時，實證發(fā)現(xiàn)經(jīng)濟差距對貿(mào)易流量的影響存在倒U型關(guān)系，經(jīng)濟差距在一定的范圍內(nèi)促進兩國貿(mào)易，反之，超出范圍后，會產(chǎn)生阻礙作用[11]。說明本文所選國家不屬于經(jīng)濟差距可促進兩國貿(mào)易的范圍內(nèi)。

7) 解釋變量r符號為負，與預(yù)期符號不一致，但p值為0.01，與中國銅礦貿(mào)易量相關(guān)性十分顯著。近年來，中國對于礦產(chǎn)勘查的投入減少，且國內(nèi)環(huán)保污染治理力度不斷增大，為滿足自身發(fā)展和資源儲備的需要，銅礦進口量不斷增加；根據(jù)MinExConsulting資料，全球銅礦資源發(fā)開周期變長，資源轉(zhuǎn)化率走低；此外，本文僅選取2012～2017年數(shù)據(jù)，雖模型擬合度較好，但筆者認為僅能反映研究期內(nèi)出口國的資源儲量與中國進口量的關(guān)系。

3 結(jié)論與建議

本文通過建立改進引力模型，對影響中國銅礦進口的因素進行了研究分析。將影響因素分為出口國影響因素和貿(mào)易空間上的市場結(jié)構(gòu)、經(jīng)濟距離、地理距離等影響因素。通過線性回歸結(jié)果的分析發(fā)現(xiàn)，觀察期內(nèi)，出口國影響因素(出口國銅礦儲量、出口國銅礦產(chǎn)量、出口國銅礦人均消費量)對中國進口影響十分顯著，貿(mào)易空間上的間接影響因素與中國銅礦進口存在相關(guān)性，但影響效果不顯著。

1) 地理距離與經(jīng)濟距離對銅礦石進口存在正U型關(guān)系。在一定的距離內(nèi)，距離對中國銅礦是進口產(chǎn)生阻礙作用；反之，超出一定距離范圍(U型拐點)，距離對中國銅礦進口的影響因子減小，出現(xiàn)由負轉(zhuǎn)正現(xiàn)象，阻礙作用減弱。

2) 在回歸結(jié)果分析中，出口國的銅礦經(jīng)濟可采儲量、產(chǎn)量、人均銅礦消費量對中國銅礦進口的影響十分顯著。根據(jù)中國自然資源部數(shù)據(jù)統(tǒng)計，2015年中國銅礦貿(mào)易對外依存度達到67%，那么針對出口國影響因素制定進口政策，有利于降低我國銅礦進口風(fēng)險。一方面，減少低品位銅礦的進口；另一方面，在符合生態(tài)環(huán)境要求的條件下，應(yīng)增加探礦權(quán)，探明儲量而暫不開發(fā)，降低銅礦進口價格，利于增加銅礦進口及儲備。

3) 貿(mào)易國家間的優(yōu)惠政策會促進銅礦進出口貿(mào)易。我國成為WTO與亞太經(jīng)濟合作組織成員國，與其他國家簽訂自由貿(mào)易協(xié)議，施行“一帶一路”倡議等，在一定程度上促進了我國與其他國家的經(jīng)濟合作，對我國礦產(chǎn)品貿(mào)易起到促進作用。