(內蒙古科技大學 經濟與管理學院,內蒙古 包頭 014010)
我國“十三五”規劃中明確提出全面節約和高效利用資源的戰略部署。貫徹落實好這一戰略部署,必須科學把握我國的資源國情,加大資源利用管理創新力度,牢固樹立節約集約循環利用的資源觀。同時,各省區間應形成上下聯動、左右互動的資源利用格局,確保全面節約和高效利用資源取得成效。我國資源總量大、人均少,分布不均衡,各省區資源儲量和資源依賴度差異很大。
市場機制與政策調解是能源資源得以跨區域配置的重要手段,因此我們有理由相信,我國資源開發在區域間存在著空間溢出和交互作用。從政府政策層面來看,我國高度重視國家資源安全,適時出臺了資源開發與管理的相關法律法規,并實施“西氣東輸”、“西電東送”等能源工程,意在促進資源開發與利用的區域協調與關聯。從市場層面來看,隨著我國社會主義市場經濟體制向更深層次發展和市場化進程的不斷加深,有利于地區間能源資源的自由流動,增強了相互間的資源聯系。因此,針對各地區在資源開發空間關聯網絡中的不同地位和作用,不斷完善政策調控機制,通過創造更多空間溢出“通道”,加強地區間的關聯程度,提升資源產業的區域協同性已經成為我國資源開發的“新常態”。在此背景下,探尋我國區域資源開發的空間溢出關聯網絡及其來源,對促進我國能源資源的區域協同發展以及實現全面節約和高效利用資源的戰略目標具有十分重要的現實意義。
目前已有學者提出我國資源產業依賴有著明顯的地區差異性和空間關聯特征。劉那日蘇、郝戊[1]在我國東、中、西三大區域層面,運用Dagum基尼系數和Kernel密度估計方法實證考察了我國資源產業依賴的地區差異與演變趨勢,指出我國資源產業依賴具有顯著的空間非均衡特征;薛稚偉等[2]基于Krugman空間基尼系數測度了我國資源產業空間集聚程度,認為我國15個省市區存在資源產業空間集聚現象,這些地區大部分為資源富集的中西部地區;鮑超等[3]運用空間計量模型方法檢驗了河南省用水效率的影響因素,結果表明河南省地級市用水效率存在顯著的空間依賴性,地區用水效率明顯受鄰近地區的影響;劉華軍等[4]利用社會網絡分析方法研究了1995—2012年我國能源消費的空間關聯特征,指出我國能源消費空間關聯呈網絡結構形態,在樣本考察期內空間關聯的緊密程度和網絡穩定性逐步增強;劉那日蘇、袁雪晴[5]運用空間Durbin模型方法從空間經濟學視角實證考察了自然資源開發對經濟增長的空間溢出效應,研究發現自然資源開發在我國各省區間存在顯著關聯。
縱觀現有文獻,空間關聯互動是我國區域資源開發中必須加以考慮的重要決策變量。針對資源開發的空間關聯,雖然現有研究已取得富有價值的研究成果,但仍存在著一定的局限。傳統空間計量方法往往將空間關聯限制在地理上相鄰或經濟上相近地區間,而我國區域資源開發利用的戰略中心是資源匱乏的東部沿海省份與資源富集的中西部省份之間的互聯互動,因此依據傳統空間計量技術得出的結論尚需進一步驗證。
本文利用格蘭杰因果關系檢驗和網絡分析方法識別我國區域資源開發的空間溢出和關聯網絡特征,并采用二次指派程序(QAP)方法探尋空間溢出的影響因素。社會網絡分析法是通過對網絡中關系的分析,刻畫其結構特征。目前該方法被廣泛應用于各類文獻中[6-8]。本文選擇網絡分析法主要基于兩點原因:一是網絡分析法能有效彌補傳統空間計量方法的不足,可全面捕捉空間關聯關系;二是網絡分析法能從整體結構、個體中心性、凝聚子群(板塊)、假設檢驗等角度揭示空間溢出的網絡結構特征與來源。
區域資源開發的空間溢出網絡是區域間自然資源開發利用的關系集,其中各區域是網絡節點,區域間的溢出關系是節點連線,因此準確測度區域資源開發的空間溢出關系是構建區域關聯網絡的關鍵。本文參考相關文獻[9,10],在非結構化的向量自回歸(VAR)模型框架下,依據HQ、SC、AIC、FPE、LR指標來選擇模型最優滯后階數,利用格蘭杰因果檢驗來實現這一目的。即如果A、B兩個區域的資源開發活動相互關聯,而且A是B的格蘭杰因果關系時,說明A省對B省有資源開發空間溢出,因此在關聯網絡中用A指向B的箭頭線將A和B連接,依此方法重復進行,檢驗所有區域兩兩之間的空間溢出關系,繪制出網絡中的節點連線,這樣便可構建我國區域資源開發的一個有向空間溢出網絡。
空間關聯網絡結構是指溢出網絡內部各區域之間相對穩定的關系模式,通常從整體網絡結構和個體網絡結構兩個維度來刻畫關聯網絡結構的特征。整體網絡結構特征通常利用網絡密度、網絡效率、網絡關聯度和網絡等級度等指標來刻畫,而個體網絡結構特征則利用度數中心度、中介中心度和接近中心度等中心性指標表示。此外,網絡分析法通常使用塊模型分析,主要解析空間聯動網絡中凝聚子群(板塊)特征與不同聚類板塊之間的空間聯動關系。本文將采用上述方法分析我國區域資源開發空間溢出網絡結構,并采用QAP回歸方法探尋空間溢出來源。
本文以我國年31個省區市(未包括香港特區、澳門特區和臺灣地區,下同)為基本區域單元開展研究。在相關文獻中,自然資源開發度一般用采礦業就業比重、投資比重、產值比重來衡量。為不失一般性,本文選取采礦業固定資產投資占全社會固定資產投資比重作為資源開發的代理變量。此外,我們從資源稟賦差異和地理近鄰效應兩個角度考察資源開發空間溢出的主要來源,選取區域標準煤產量作為資源稟賦的測量指標,依據空間上是否相鄰和省會城市的球面距離來構建地理近鄰效應矩陣。囿于數據的可得性和統計口徑變化,本文考察的時間范圍為2003—2015年,數據來源于相關年份的《中國統計年鑒》和各省統計年鑒。
序列平穩是進行Granger因果檢驗的必要前提。本文利用ADF單位根檢驗方法對各省資源開發度序列進行了平穩性檢驗,結果顯示一階差分序列在10%的顯著性水平上均是平穩的。由此,本文構造了31個省區市兩兩之間的資源開發VAR模型,在此基礎上利用Granger因果檢驗來識別空間溢出關系,構建空間關聯網絡,并利用Netdrau工具對空間關聯網絡進行可視化,結果見圖1。從圖1可見,圖中每個箭頭代表區域資源開發的空間溢出關系,根據可視化網絡關聯線條,可初步判斷我國區域資源開發普遍存在空間溢出。

圖1 我國區域資源開發空間溢出網絡
整體網絡結構特征:我國區域資源開發空間關聯網絡結構指標測度見表1。從網絡密度來看,31個省區市的實際溢出關系為319個,理論上應存在的可能最大溢出關系數為930個,整體網絡密度為0.3430,說明我國資源開發的區域關聯較普遍。但從網絡密度值來看,關聯關系的緊密程度整體相對不高,全面提高和促進更緊密的資源開發協作尚有較大的空間和潛力。從網絡關聯度來看,資源開發的區域溢出網絡關聯度為1,說明省份之間資源開發的空間關聯可達性較強,存在著普遍的空間溢出效應,不存在“孤立”的省份,每個省的資源開發活動均受到空間關聯網絡的影響。從網絡效率來看,測算值為0.4552,表明區域資源開發的空間溢出存在較多冗余連線,整體網絡較穩定。同時,測算得到的網絡等級度為0,說明區域資源開發空間溢出無等級屬性,在不同資源開發水平上都有可能產生溢出效應。
個體網絡結構特征:區域資源開發溢出是有向關聯網絡,因此有必要考察各區域在網絡中的地位和作用。一般利用度數中心度、接近中心度和中介中心度等指標揭示網絡節點的個體網絡結構特征,本文中心度測度結果見表1。需要說明的是,為了便于相互比較,表1是標準化處理后的相對中心度測度結果。從表1可知,31個省區市在資源開發空間關聯網絡中扮演著不同的角色和作用。從度數中心度看,31個省區市平均度數中心度為57.419。其中,15個超過平均值,從高到低依次是黑龍江、吉林、浙江、寧夏、湖北、青海、新疆、河北、安徽、山東、遼寧、甘肅、天津、湖南、江西,說明這些省份在網絡中處于重要地位,在區域資源開發關聯網絡中與其發生的直接關系最多。同時還發現,上述省區市中除浙江、遼寧、河北、天津外,其他省區市均屬于中西部地區,進一步證實了中西部地區對我國資源開發空間網絡具有舉足輕重的影響。從中介中心度來看,黑龍江、吉林、寧夏、湖北、浙江較高,表明這些省份在資源開發空間關聯網絡中控制其他省份之間資源開發交流的能力較強,而且這5個省份中的中部省份有3個,由此說明中部地區在資源開發空間關聯網絡中起著重要的“橋梁”和“中介”作用。從接近中心度而言,全國平均值為70.908,黑龍江、寧夏、浙江、吉林、新疆、青海、湖北、山東、安徽、河北、甘肅、遼寧、湖南、天津、江西的接近中心度超過了全國平均值,表明上述省份在資源開發空間關聯網絡中扮演著中心行動者的角色。

表1 中國區域資源開發的空間溢出網絡結構特征
我們采用CONCOR方法[11,12](Convergent Correlations),選擇最大切分深度為2,集中標準為0.2,將我國31個省區市歸類為4個聚類板塊,結果見表2。從表2可見,聚類板塊內部溢出關系為85個,板塊間溢出關系為234個,板塊間關系占整體網絡關系的73%,說明聚類板塊間的溢出效應較突出。具體來看,聚類板塊Ⅰ由北京、上海、福建、云南、內蒙古、重慶、廣西組成,向外發出58個溢出關系,接收外部板塊關系35個,板塊內部溢出關系7個,期望內部關系比例為20%,實際內部關系比例為11%,可將其劃分為“凈溢出”聚類板塊。聚類板塊Ⅱ由黑龍江、天津、江蘇、青海、安徽、貴州、遼寧、山東、河南、新疆構成,向外溢出關系95個,接收外部板塊關系65個,內部成員關系42個,實際內部關系比例與期望內部關系比例相當,屬于“雙向溢出”聚類板塊,其成員既發送也接受外部關系,聚類板塊內部省份之間的聯系相對較多。聚類板塊Ⅲ由湖南、吉林、浙江、江西、河北、湖北、寧夏、甘肅構成,內部關系25個,接收外部板塊關系105個,向外溢出關系53個,實際內部關系比例大于期望內部比例,屬于“凈收益”聚類板塊。該聚類板塊既發送也接受外部關系,但后者明顯多于前者。聚類板塊Ⅳ由廣東、山西、海南、陜西、四川、西藏構成,聚類板塊內部關系11個,接收其他板塊關系29個,向外溢出關系28個,內部實際關系比例大于期望內部關系比例,屬于“經紀人”板塊,該板塊成員在網絡中既向外發送關系,也接受外部關系,但板塊內部成員間關系較少,在網絡中主要發揮“中介”和“橋梁”作用。

表2 區域資源開發的空間溢出聚類板塊劃分及其溢出效應
按照網絡密度含義,聚類板塊的網絡密度大于整體網絡密度,說明溢出關系多集中于該聚類板塊。我們以聚類板塊網絡密度大于整體網絡密度的賦值為1,否則為0的原則,將聚類板塊網絡密度矩陣轉化為像矩陣,進而更清晰地刻畫各聚類板塊之間的溢出效應,結果見表3。從表3可見,除聚類板塊Ⅰ外,其他聚類板塊的像矩陣對角線元素均為1,說明板塊內部省份之間存在明顯的資源開發聯動關系,即“俱樂部”效應較突出。聚類板塊Ⅰ向聚類板塊Ⅱ產生溢出。因此,聚類板塊Ⅰ和聚類板塊Ⅳ扮演著重要的“發動機”角色,聚類板塊Ⅱ在資源依賴關聯網絡中同樣發揮著重要的中介和橋梁作用。綜上所述,加強聚類板塊Ⅰ自身關聯、聚類板塊Ⅰ與聚類板塊Ⅳ、聚類板塊Ⅳ與聚類板塊Ⅱ的關聯互動,將對我國資源開發戰略的制定和實現區域協同發展具有重要的決策參考意義。

表3 資源開發空間溢出板塊的密度矩陣與像矩陣
學者們通常將地理距離或空間近鄰視為空間溢出的關鍵影響因素[13],即某區域的資源開發活動可能會影響到與之相鄰或近距離的其他區域。此外,一般而言資源豐裕地區更易于產生資源依賴心理,往往圍繞所謂“資源優勢”開展經濟建設,有意無意被鎖定于資源依賴型發展道路,進而長期熱衷于資源開發活動。因此,我們做出如下理論假設:我國區域資源開發的空間關聯關系R主要受兩個關系因素的影響,分別是由地理位置所決定的近鄰關系S和由資源稟賦差異K。據此,建立如下模型:
R=f(S,K)
(1)
式中,R為資源開發空間溢出矩陣;S為資源稟賦差異矩陣;K為空間距離矩陣。選取鄰接空間矩陣和地理距離矩陣來衡量空間距離效應,若兩省在空間上相鄰則為1,否則為0,依此構建空間鄰接矩陣,表示為I;選取0—500km、500—1000km兩個區間來測度地理距離效應,若省會城市距離在該區間內為1,否則為0,依此構建地理距離矩陣。關于資源稟賦指標S,選取各省煤炭、石油、天然氣預測儲量,并按照標準煤折算系數轉化為標準煤儲量作為其代理變量,然后取實證期間各省平均值的絕對差組建差異矩陣。
由于模型(1)中變量均為關系數據,由于關系數據之間本身可能存在的相關性問題,若直接采用傳統回歸方法,可能會得到有偏估計量[14]。為此,本文選用QAP回歸分析方法對區域資源開發的影響因素進行參數估計。QAP作為一種非參數估計方法,無需假設自變量之間相互獨立,且該方法為網絡分析中常用的方法之一,分析結果更具有穩健性和邏輯說服力。
QAP回歸分析的目的主要是考察多個影響因素矩陣與區域資源開發空間溢出矩陣之間的關聯關系,并對模型可決系數R2的顯著性進行評價。區域資源開發空間溢出矩陣與影響因素差異矩陣的QAP相關性分析結果見表4。其中,空間鄰接矩陣、0—500km、500—1000km距離矩陣的相關系數分別為0.013、0.011和0.025,但均未通過10%的顯著性水平,說明地理位置鄰接性與地理距離近鄰性對區域資源開發的空間溢出并沒有十分顯著的影響。資源稟賦差異矩陣與區域資源開發空間溢出的相關系數為-0.113,且通過了5%的顯著性水平檢驗,說明區域資源稟賦差異對區域資源開發空間溢出具有顯著影響。因此,本文進一步就資源稟賦差異對資源開發空間溢出的影響進行QAP回歸分析。

表4 QAP相關性分析結果
本文選擇5000次隨機置換,結果見表5。其中,概率1表示隨機置換得到的回歸系數不小于實際觀察到的回歸系數的概率,概率2表示隨機置換得到的回歸系數不大于實際觀察到的回歸系數的概率。從表5可知,資源稟賦差異矩陣S的回歸系數為負,且通過了5%的顯著性水平檢驗,表明資源稟賦相似是發生區域資源開發空間溢出的主要原因之一。此外,調整后的可決系數R2為0.132,說明資源稟賦差異可解釋資源開發空間溢出的13.1%。

表5 QAP回歸結果
本文得到以下結論:①我國區域資源開發存在顯著空間溢出特征,呈穩定、多線程的溢出網絡結構狀態,任何區域的資源開發活動均受到整體溢出網絡的影響。②各區域在網絡中具有不同的地位和作用,中西部省份對資源開發空間溢出網絡具有舉足輕重的影響。其中,中部地區省份在網絡中發揮著重要的橋梁作用,同時空間溢出沒有等級屬性,在不同資源開發水平上都有可能發生溢出。③北京、上海、福建、云南、內蒙古、重慶、廣西組成的聚類板塊在資源開發空間溢出網絡中扮演“凈溢出”角色,由黑龍江、天津、江蘇、青海、安徽、貴州、遼寧、山東、河南、新疆組成的聚類板塊在網絡中扮演“雙向溢出”角色,由湖南、吉林、浙江、江西、河北、湖北、寧夏、甘肅組成的聚類板塊在溢出網絡中扮演“凈收益”角色,由廣東、山西、海南、陜西、四川、西藏組成的聚類板塊在網絡中扮演“經紀人”角色。④地理臨近效應對資源開發的空間溢出無顯著影響,而區域資源稟賦差異對資源開發的空間溢出具有顯著的負向影響,即資源稟賦的相似性是區域之間發生資源開發空間溢出的重要原因之一。
首先,本文測度了我國區域資源開發空間溢出網絡的整體結構特征,揭示了各區域在溢出網絡中的中心特征,為能源開發聯合聯動機制的構建提供了重要決策依據。其次,根據資源開發空間溢出網絡特征,我國有必要將省際空間關聯作為區域能源協調的重要決策變量,進行定向調控和精準調控,創造更多的空間關聯“通道”,提高省域間的關聯緊密程度,提升能源資源的空間配置效率。同時,要根據不同板塊功能和網絡個體特征,選擇有針對性的資源發展政策,以提升區域資源開發的空間協同性。第三,各省既要重視與地理位置上近鄰省份間的互聯互動,更要重視空間關聯網絡所呈現的復雜多線程聯動關系,加強資源溢出的空間聯動效應,推動資源合作從“局部”轉向“整體”,從“點”轉向“面”,實現能源資源的區域協同發展。