基于改進(jìn)型引力模型的中國銅礦石進(jìn)口影響因素實(shí)證分析

賈祥英，邢萬里，苗媛媛

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

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

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

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

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

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

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

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

1.2 研究方法

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

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

表1 2012～2017年中國主要進(jìn)口來源國銅礦貿(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國進(jìn)口的貿(mào)易量；G為常量；Mi、Mj為參與貿(mào)易的國家的GDP；D為貿(mào)易兩國之間的距離；?ij為誤差項(xiàng)；β0、β1、β2、β3為相關(guān)系數(shù)。

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

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

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

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

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

(3)

2 改進(jìn)引力模型回歸方程分析

2.1 回歸結(jié)果分析

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

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

2.2 實(shí)證結(jié)果分析

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

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

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

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

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

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

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

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

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

3 結(jié)論與建議

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

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

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

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