黃河流域采礦業(yè)集聚的時空演化特征與驅(qū)動力研究

　　摘 要: 黃河流域礦產(chǎn)資源豐富，合理有序開發(fā)礦產(chǎn)資源對黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展意義重大。綜合運(yùn)用區(qū)位商、空間基尼系數(shù)、探索性空間數(shù)據(jù)分析、標(biāo)準(zhǔn)差橢圓等方法分析了 1999— 2018 年黃河流域采礦業(yè)集聚的時空演化特征，并對其驅(qū)動力進(jìn)行地理探測。研究發(fā)現(xiàn): 黃河流域采礦業(yè)集聚程度總體呈先上升、后波動下降、再快速上升的 “N”型發(fā)展態(tài)勢，1998 年和 2016 年兩輪去產(chǎn)能政策顯著提高了黃河流域采礦業(yè)集聚程度; 黃河中上游地區(qū)采礦業(yè)集聚程度明顯高于下游地區(qū)，采礦業(yè)空間相關(guān)性逐漸增強(qiáng)， “高-高”類型城市集中分布在山西和陜西，空間范圍不斷擴(kuò)大; 黃河流域采礦業(yè)空間布局方向性明顯，基本保持 “東—西”方向，研究期內(nèi)采礦業(yè)產(chǎn)業(yè)重心均落入山西省境內(nèi)且總體向西北方向遷移; 黃河流域采礦業(yè)集聚的核心驅(qū)動因素為礦產(chǎn)資源稟賦、政府干預(yù)、經(jīng)濟(jì)發(fā)展水平和下游產(chǎn)業(yè)需求，上、中、下游的核心驅(qū)動因子存在顯著差異。

　　宋梅; 郝旭光; 李丹丹; 常力月， 煤炭工程 發(fā)表時間：2021-09-16

　　關(guān)鍵詞: 黃河流域; 采礦業(yè); 產(chǎn)業(yè)集聚; 時空演化; 驅(qū)動力; 地理探測器

　　黃河流域煤炭、石油、天然氣和有色金屬等礦產(chǎn)資源豐富，是我國重要的能源、化工、原材料和基礎(chǔ)工業(yè)基地[1]。在我國七大江河流域中，黃河流域能源資源具有種類齊全、儲量大、開采條件優(yōu)越等特點(diǎn)，形成了上游水電、中游煤炭、下游石油的分布格局，黃河流域因此也被稱為 “能源流域”[2]。采礦業(yè)作為黃河流域經(jīng)濟(jì)社會發(fā)展的重要支柱產(chǎn)業(yè)，在全國能源和原材料供應(yīng)方面具有重要的戰(zhàn)略地位。自黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展上升為重大國家戰(zhàn)略以來，學(xué)者圍繞黃河流域礦產(chǎn)資源開發(fā)利用開展了大量研究。研究涉及礦區(qū)生態(tài)修復(fù)[3]、礦產(chǎn)開發(fā)空間布局[4]、綜合能源效率[2]、資 源 型 城 市 轉(zhuǎn)型[5]等。黃河流域礦產(chǎn)資源豐富，但流域內(nèi)資源分布并不均衡，礦業(yè)集聚特征明顯。現(xiàn)有關(guān)于產(chǎn)業(yè)集聚的研究主要包括集聚測度和影響因素分析。產(chǎn)業(yè)集聚測度的方法有區(qū)位商、空間基尼系數(shù)、赫芬達(dá)爾指數(shù)和 EG 指數(shù)等。董雯等[6]采用區(qū)位商等方法研究烏魯木齊資源型產(chǎn)業(yè)的演變特征，毛熙彥等[7]借助空間基尼系數(shù)、頡茂華等[8]采用赫芬達(dá)爾-赫希曼指數(shù)和 EG 指數(shù)分別研究了我國市域和省域資源型產(chǎn)業(yè)集聚的時空特征。在測度產(chǎn)業(yè)集聚程度的基礎(chǔ)上，金煜等[9]將影響工業(yè)集聚的因素分為經(jīng)濟(jì)地理、新經(jīng)濟(jì)地理與經(jīng)濟(jì)政策三類并進(jìn)行面板回歸分析，賀燦飛等[10]采用 Tobit 模型發(fā)現(xiàn)資源稟賦、歷史基礎(chǔ)、交通通達(dá)性、市場規(guī)模等對資源密集型產(chǎn)業(yè)分布影響顯著，朱俏俏等[11]對比研究了我國資源型產(chǎn)業(yè)和制造業(yè)的集聚態(tài)勢和影響因素。

　　綜上所述，現(xiàn)有關(guān)于資源型產(chǎn)業(yè)集聚的研究多從全國尺度開展，黃河流域礦產(chǎn)資源豐富，目前較少有文獻(xiàn)從流域視角定量研究其采礦業(yè)集聚的時空特征。本文基于 1999—2018 年黃河流域 90 個地級行政區(qū)采礦業(yè)從業(yè)人員數(shù)據(jù)，借助 ArcGIS 工具的可視化優(yōu)勢，綜合運(yùn)用區(qū)位商、空間基尼系數(shù)、探索性空間數(shù)據(jù)分析、標(biāo)準(zhǔn)差橢圓分析，對近 20 年黃河流域采礦業(yè)集聚的時空演化特征進(jìn)行分析，并運(yùn)用地理探測器模型從全域和局域兩個視角探測其核心驅(qū)動因素，以期為優(yōu)化黃河流域礦產(chǎn)資源開發(fā)布局，推動能源行業(yè)高質(zhì)量發(fā)展提供決策參考。

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

　　1. 1 研究范圍

　　黃河流經(jīng)青海、四川、甘肅、寧夏、內(nèi)蒙古、山西、陜西、河南、山東九省區(qū)，但四川全省已劃入長江經(jīng)濟(jì)帶，蒙東五盟市( 呼倫貝爾市、興安盟、通遼市、赤峰市和錫林郭勒盟) 已納入 《東北振興規(guī)劃》。因此，本研究將黃河流域的空間范圍界定為青海、甘肅、寧夏、內(nèi)蒙古( 不含蒙東五盟市) 、山西、陜西、河南、山東八省區(qū)，共 90 個地級行政區(qū)。參考趙明亮等[12] 的 研 究，將 青 海、甘 肅、寧夏、內(nèi)蒙古所涉及城市均歸為黃河上游，陜西、山西全部城市歸為黃河中游，河南、山東全部城市歸為黃河下游。本研究中采礦業(yè)的范圍與 《國民經(jīng)濟(jì)行業(yè)分類》( GB /T 4754—2017) 中的劃分保持一致，包括煤炭開采和洗選業(yè)、石油和天然氣開采業(yè)、黑色金屬礦采選業(yè)、有色金屬礦采選業(yè)、非金屬礦采選業(yè)、開采專業(yè)及輔助性活動以及其它采礦業(yè)。

　　1. 2 研究方法

　　1. 2. 1 區(qū)位商

　　區(qū)位商是衡量單一行業(yè)在特定區(qū)域?qū)I(yè)化與集聚水平的一項指數(shù)，該方法操作簡單易行，應(yīng)用廣泛，適用于測度研究區(qū)域內(nèi)特定產(chǎn)業(yè)在各地理單元的集聚程度。區(qū)位商越大，表明產(chǎn)業(yè)集中度越高。其計算公式為: LQij = Lij /Li Lj /L ( 1) 其中，LQij為地理單元 i 產(chǎn)業(yè) j 的區(qū)位商，Lij為地理單元 i 產(chǎn)業(yè) j 城鎮(zhèn)單位從業(yè)人員數(shù)，Li 為地理單元 i 全部產(chǎn)業(yè)城鎮(zhèn)單位從業(yè)人員數(shù)，Lj 為全國產(chǎn)業(yè) j 城鎮(zhèn)單位從業(yè)人員數(shù)，L 為全國所有產(chǎn)業(yè)城鎮(zhèn)單位從業(yè)人員數(shù)。

　　1. 2. 2 空間基尼系數(shù)

　　基尼系數(shù)最初被用來衡量收入不均等程度，后來國外學(xué)者將其應(yīng)用于區(qū)域產(chǎn)業(yè)空間分布問題的研究。Krugman 在研究美國制造業(yè)集聚問題時提出了空間基尼系數(shù)，該系數(shù)的取值范圍為[0，1]，值越大表示該產(chǎn)業(yè)在地理上的集聚程度越高。空間基尼系數(shù)考慮了不同地區(qū)地理單元面積大小的差異對產(chǎn)業(yè)集中度的影響，同時考慮了全部行業(yè)的地理分布，這使得不同行業(yè)的空間基尼系數(shù)的計算結(jié)果具有可比性。因此，本文采用空間基尼系數(shù)從全局角度測度黃河流域采礦業(yè)的集聚程度。參考 Krugman [13]的研究，空間基尼系數(shù)的計算公式為: Gj = ∑ n i = 1 ( Sij - Zi ) 2 ( 2) 其中，Gj 為研究區(qū)域產(chǎn)業(yè) j 的空間基尼系數(shù)， Sij為地理單元 i 產(chǎn)業(yè) j 城鎮(zhèn)單位從業(yè)人員數(shù)占研究區(qū)域產(chǎn)業(yè) j 城鎮(zhèn)單位從業(yè)人員數(shù)的比重，Zi 為地理單元 i 城鎮(zhèn)單位從業(yè)人員數(shù)與研究區(qū)域城鎮(zhèn)單位從業(yè)人員數(shù)的比重。

　　1. 2. 3 探索性空間數(shù)據(jù)分析

　　空間自相關(guān)作為有效檢驗空間單元屬性值關(guān)聯(lián)程度的度量方法，分為全局和局域兩類[14]。采用全局 Moran’s I 分析黃河流域采礦業(yè)集聚程度在空間上的相似性。全局 Moran’s I 的范圍為[-1，1]，越接近 1，說明探索數(shù)據(jù)在空間上表現(xiàn)出明顯的正相關(guān)性; 越接近-1，說明探索數(shù)據(jù)在空間上表現(xiàn)出明顯的負(fù)相關(guān)性。全局 Moran’s I 的計算公式為: I = ∑ n i = 1 ∑ n j = 1 wij ( xi - x) ( xj - x) S2 ∑ n i = 1 ∑ n j = 1 wij ( 3) 其中，S2 = ∑ n i = 1 ( xi-x) 2 ，x = 1 n ∑ n i = 1 xi，n 為研究區(qū)域內(nèi)地理單元數(shù)，wij為研究區(qū)域內(nèi)地理單元 i 與 j 的鄰接關(guān)系，xi 是地理單元 i 的觀測值，x 是市域觀測值的平均值，本研究的觀測值為各地理單元采礦業(yè)區(qū)位商。

　　采用局部 Moran’s I 分析黃河流域某一地理單元與其臨近地理單元采礦業(yè)聚集程度之間的相似性。Anselin 等提出局部空間自相關(guān)方法來衡量某一區(qū)域與周圍地區(qū)之間的空間差異程度及顯著性，并結(jié)合 LISA 聚集圖等研究空間分布規(guī)律。LISA 集聚圖中城市可分為 “高-高”“高-低”“低-高”“低-低”四種類型，“高-高”類型是指城市自身的采礦業(yè)集聚度高同時相鄰城市的采礦業(yè)集聚度也高， “高-低” 類型是指城市自身的采礦業(yè)集聚度高而相鄰城市的采礦業(yè)集聚度較低， “低-高”類型是指城市自身的采礦業(yè)集聚度低而相鄰城市的采礦業(yè)集聚度較高， “低-低”類型是指城市自身的采礦業(yè)集聚度低同時相鄰城市的采礦業(yè)集聚度也較低。局部 Moran’s I 的計算公式為: Ii = xi - x S2 ∑ n j = 1 wij ( xj - x) ( 4) 在特定的置信水平下，局部 Moran’s I 大于 0 時，存在局部空間自相關(guān)性，表示相似值發(fā)生聚集; 局部 Moran’s I 小于 0 時，存在負(fù)的局部空間自相關(guān)性，表示不相似的值發(fā)生聚集。

　　1. 2. 4 標(biāo)準(zhǔn)差橢圓分析

　　標(biāo)準(zhǔn)差橢圓是分析空間分布方向性特征的經(jīng)典方法，采用標(biāo)準(zhǔn)差橢圓分析黃河流域采礦業(yè)產(chǎn)業(yè)重心遷移與展布態(tài)勢。標(biāo)準(zhǔn)差橢圓主要元素的計算公式為:其中，xi、yi 分別為第 i 個地理單元的經(jīng)緯度坐標(biāo); X 和 Y 分別為標(biāo)準(zhǔn)差橢圓中心的橫、縱坐標(biāo)，用來確定產(chǎn)業(yè)重心; wi 表示不同地理單元的權(quán)重; θ 為標(biāo)準(zhǔn)差橢圓的方向角度，是由正北沿順時針方向與橢圓長軸的夾角; σx 和 σy 分別表示標(biāo)準(zhǔn)差橢圓的長、短軸的長度。

　　1. 2. 5 地理探測器

　　地理探測器是探測空間分異性以及揭示其背后驅(qū)動力的一組統(tǒng)計學(xué)方法[16]。與傳統(tǒng)方法相比，該模型具有三個優(yōu)勢: 從異質(zhì)性視角探究影響因子的空間分異特征; 對多自變量共線性免疫，并能較好地避免自變量和因變量互為因果的內(nèi)生性問題; 不僅能探測定量因子，也能對多種類型的定性驅(qū)動因子進(jìn)行識別分析[17]。因子探測模型如下: q = 1 - 1 Nσ2∑ L h = 1 Nhσ2 h ( 9) 其中，L 為自變量的分類; Nh 和 σh 分別為層 h 內(nèi)的單元數(shù)和方差; N 和 σ 分別為研究區(qū)整體的單元數(shù)和方差。q 的取值范圍為[0，1]，值越大說明自變量對因變量的解釋力越強(qiáng)。

　　產(chǎn)業(yè)集聚的影響因素相當(dāng)復(fù)雜，參考金煜等[10]的研究，將影響采礦業(yè)集聚的因素分為傳統(tǒng)的經(jīng)濟(jì)地理因素、新經(jīng)濟(jì)地理因素和經(jīng)濟(jì)政策因素三類，選擇 10 個候選因子對黃河流域采礦業(yè)集聚時空分異的驅(qū)動力進(jìn)行探測，黃河流域采礦業(yè)集聚的探測因子見表 1。

　　1. 3 數(shù)據(jù)來源

　　黃河流域 90 個地級行政區(qū)采礦業(yè)城鎮(zhèn)單位從業(yè)人員數(shù)據(jù)及探測因子數(shù)據(jù)來源于 《中國城市統(tǒng)計年鑒( 2000—2019》、黃河流域各省份歷年統(tǒng)計年鑒及各地市統(tǒng)計公報，部分年份缺失的數(shù)據(jù)通過插值法補(bǔ)齊。

　　2 黃河流域采礦業(yè)集聚的時空演化格局分析

　　2. 1 時序演化特征分析

　　1999—2018 年黃河流域采礦業(yè)與制造業(yè)空間基尼系數(shù)如圖 1 所示，黃河流域采礦業(yè)空間基尼系數(shù)呈 “N”型變化趨勢，集聚程度總體呈先上升、后波動下降、再快速上升的發(fā)展態(tài)勢。受 1998 年去產(chǎn)能政策的影響，1999—2002 年黃河流域采礦業(yè)集聚程度提升明顯; 2003—2012 年黃河流域采礦業(yè)集聚程度呈波動下降態(tài)勢，究其原因，2003 年全球資源型產(chǎn)業(yè)進(jìn)入繁榮期，各地大力發(fā)展采礦業(yè)，產(chǎn)業(yè)集聚程度持續(xù)下降，受 2008 年全球金融危機(jī)影響，資源型產(chǎn)品需求低迷，大量煤礦關(guān)停，礦業(yè)集聚程度有所提高，2010 年以后隨著經(jīng)濟(jì)的復(fù)蘇，各地煤礦生產(chǎn)開始恢復(fù)，礦業(yè)集聚程度又出現(xiàn)了下降; 2013 年資源市場繁榮期結(jié)束，煤炭產(chǎn)能過剩嚴(yán)重，為了化解過剩產(chǎn)能，2016 年我國煤炭行業(yè)開始實施供給側(cè)改革，出臺了一系列去產(chǎn)能政策[19，20]。伴隨著我國西南、東部等地區(qū)落后煤炭產(chǎn)能的退出和黃河流域優(yōu)質(zhì)煤炭產(chǎn)能的釋放，黃河流域礦業(yè)集聚程度快速提高。對比采礦業(yè)和制造業(yè)的基尼系數(shù)發(fā)現(xiàn): 黃河流域采礦業(yè)的空間聚集程度遠(yuǎn)高于制造業(yè)，一定程度上反映出制造業(yè)更易于通過交通、資源來源的替代性擺脫空間約束，而采礦業(yè)由于生產(chǎn)投入要素在空間上存在不可移動性，因此更容易形成空間集聚。

　　2. 2 空間演化特征分析

　　2. 2. 1 空間異質(zhì)性分析

　　為明晰黃河流域市域采礦業(yè)集聚的空間異質(zhì)性，借助 ArcGIS10. 3 對采礦業(yè)區(qū)位商進(jìn)行分級可視化，如圖 2 所示，不同等級區(qū)位商的空間分布穩(wěn)中有變，中上游地區(qū)采礦業(yè)集聚程度明顯高于下游地區(qū)。LQ ≤1 說明該城市的采礦業(yè)集聚水平低于或接近全國平均水平，此類城市數(shù)量基本穩(wěn)定，約占黃河流域城市的 50%。LQ>1 說明該城市采礦業(yè)集聚水平高于全國平均水平，特別是當(dāng) LQ>3 時，表現(xiàn)出較高的集聚水平，此類城市表現(xiàn)為中游帶狀集中、上下游零星點(diǎn)狀分布的空間特征，數(shù)量由 1999 年的 18 個增加到 2018 年的 30 個，此類城市基本全部為 《全國資源型城市可持續(xù)發(fā)展規(guī)劃( 2016—2020) 》確定的資源型城市。當(dāng) LQ>5 時，表現(xiàn)出極高的集聚水平，此類城市逐步向中游地區(qū)集中，且以成熟型資源型城市為主。

　　2. 2. 2 空間關(guān)聯(lián)特征分析

　　借助 GeoDa 軟件，對 1999—2018 年黃河流域采礦業(yè)區(qū)位商進(jìn)行全局空間自相關(guān)分析，黃河流域采礦業(yè)區(qū)位商全局 Moran’s I 指數(shù)見表 2。由表 2 可知，黃河流域采礦業(yè)區(qū)位商的 Moran’I 指數(shù)均大于 0，且總體上通過了 10%的顯著性水平檢驗，說明黃河流域采礦業(yè)集聚存在顯著的空間自相關(guān)性。研究期內(nèi)黃河流域采礦業(yè)區(qū)位商全局 Moran’s I 呈波動增長態(tài)勢，空間關(guān)聯(lián)性逐步增強(qiáng)。

　　全局 Moran’s I 只能總體判斷黃河流域采礦業(yè)的空間關(guān)聯(lián)特征，為了檢驗局部地理單元采礦業(yè)集聚與其周邊地理單元的關(guān)聯(lián)性，繪制黃河流域采礦業(yè) LISA 集聚情況如圖 3 所示。由圖 3 可知，不同關(guān)聯(lián)類型城市的數(shù)量和變化趨勢均有較大差異。具體來看， “高-高”類型城市集中分布在山西省和陜西省，空間范圍不斷擴(kuò)大，由 1999 年的山西臨汾逐步增加為山西臨汾、忻州、呂梁、太原、晉中、長治和陜西的榆林和延安; “低-低”類型的城市空間格局中下游相對穩(wěn)定，中上游逐步減少; “低-高”類型的城市由忻州市和安陽市變?yōu)檫\(yùn)城市; “高-低” 類型的城市數(shù)量最少，僅平頂山和東營在 2018 年表現(xiàn)為此類特征。

　　2. 2. 3 產(chǎn)業(yè)重心遷移與展布態(tài)勢

　　根據(jù)黃河流域 90 個城市的地理坐標(biāo)和采礦業(yè)區(qū)位商，運(yùn)用 ArcGIS10. 3 軟件計算黃河流域采礦業(yè)重心及標(biāo)準(zhǔn)差橢圓，并對其進(jìn)行疊加比較，從而得到不同時期黃河流域采礦業(yè)產(chǎn)業(yè)重心遷移路徑與標(biāo)準(zhǔn)差橢圓趨勢，如圖 4 所示。

　　產(chǎn)業(yè)重心遷移路徑方面，黃河流域采礦業(yè) 5 個觀察年份的產(chǎn)業(yè)重心均落入山西省境內(nèi)。黃河流域采礦業(yè)產(chǎn)業(yè)重心在南北方向變化不明顯，東西方向有明顯的遷移，1999—2018 年黃河流域采礦業(yè)產(chǎn)業(yè)重心總體往西北方向偏移，僅 2003 年出現(xiàn)了向東南方向偏移的趨勢，這很可能與當(dāng)時資源產(chǎn)業(yè)進(jìn)入繁榮期，中下游城市發(fā)展采礦業(yè)有關(guān)。

　　標(biāo)準(zhǔn)差橢圓趨勢方面，1999—2018 年方向角范圍為 89. 24° ～ 94. 73°，基本保持東—西方向。黃河流域采礦業(yè)在主軸方向上呈現(xiàn)出持續(xù)收縮—快速擴(kuò)張—逐步收縮的階段性特征。1999—2008 年主軸由 683km 減至 645km 減小了 38km，2008—2013 年，主軸擴(kuò)張了 54km，2013—2018 年，主軸減少了 10km。黃河流域采礦業(yè)在輔軸方向上呈現(xiàn)出持續(xù)收縮的態(tài)勢，輔軸由 1999 年的 298km 下降至 2018 年的 286km，年均收縮 0. 6km。長期來看主軸和輔軸均呈收縮趨勢，進(jìn)一步說明黃河流域采礦業(yè)在空間上更加集聚，產(chǎn)業(yè)集中度不斷提高。

　　3 黃河流域采礦業(yè)集聚的驅(qū)動力解析

　　3. 1 全域因子探測結(jié)果分析

　　因子探測結(jié)果見表 3，各因子對黃河流域采礦業(yè)集聚均有一定的驅(qū)動作用。全域視角下，黃河流域采礦業(yè)集聚的核心驅(qū)動因素是礦產(chǎn)資源稟賦、政府干預(yù)、經(jīng)濟(jì)發(fā)展水平和下游產(chǎn)業(yè)需求。

　　礦產(chǎn)資源稟賦的 q 值遠(yuǎn)大于其他因子，對黃河流域采礦業(yè)集聚起到?jīng)Q定性作用。一方面，比較優(yōu)勢引導(dǎo)產(chǎn)業(yè)區(qū)位，采礦業(yè)的發(fā)展高度依賴礦產(chǎn)資源，豐富的礦產(chǎn)資源是采礦業(yè)集聚的先決條件; 另一方面，資源稟賦初始差異導(dǎo)致采礦業(yè)后續(xù)的投資與技術(shù)積累將集中在特定區(qū)域。

　　政府干預(yù)對黃河流域采礦業(yè)集聚有顯著的影響。大型資源型企業(yè)的布局是國家早期建設(shè)時期為發(fā)展基礎(chǔ)產(chǎn)業(yè)而做出的戰(zhàn)略安排，我國很多資源型城市是在中央政府的支持下建成的，這樣的政府行為至今仍然主導(dǎo)著采礦業(yè)的地理格局。此外，采礦業(yè)稅收貢獻(xiàn)大，政府對城市產(chǎn)業(yè)發(fā)展干預(yù)較高，經(jīng)濟(jì)發(fā)展圍繞資源型產(chǎn)業(yè)，由于路徑依賴和擠出效應(yīng)強(qiáng)化了以采礦業(yè)為核心的資源型產(chǎn)業(yè)的空間集聚。

　　經(jīng)濟(jì)社會發(fā)展水平對黃河流域采礦業(yè)集聚的作用顯著。根據(jù)配第-克拉克定理，隨著經(jīng)濟(jì)的發(fā)展和國民收入水平的提高，勞動力會逐步由第一產(chǎn)業(yè)向第二產(chǎn)業(yè)和第三產(chǎn)業(yè)轉(zhuǎn)移。對于黃河下游地區(qū)，其經(jīng)濟(jì)社會發(fā)展水平較高，產(chǎn)業(yè)結(jié)構(gòu)整體處于第二產(chǎn)業(yè)主導(dǎo)向第三產(chǎn)業(yè)主導(dǎo)的轉(zhuǎn)換階段，在資源枯竭和環(huán)境規(guī)制等因素的共同作用下，采礦業(yè)向中上游轉(zhuǎn)移。而中上游地區(qū)經(jīng)濟(jì)社會發(fā)展相對落后，產(chǎn)業(yè)結(jié)構(gòu)高級度相對較低，礦產(chǎn)資源豐富，當(dāng)?shù)卣邪l(fā)展采礦業(yè)的主觀意愿和客觀條件，使得黃河流域中上游采礦業(yè)集聚程度較高。下游產(chǎn)業(yè)需求對黃河流域采礦業(yè)集聚有明顯的驅(qū)動作用。企業(yè)前后向關(guān)聯(lián)導(dǎo)致的產(chǎn)業(yè)外部性使企業(yè)集聚在其他企業(yè)周圍。資源型產(chǎn)品大多是面向本地市場，相比于勞動密集型產(chǎn)業(yè)的產(chǎn)品，資源型產(chǎn)品具有更高的內(nèi)向度。市場大的地方更容易形成規(guī)模經(jīng)濟(jì)效應(yīng)，帶來大的市場需求，有利于大型資源密集型企業(yè)的發(fā)展。

　　3. 2 局域因子探測結(jié)果分析

　　黃河流域采礦業(yè)集聚驅(qū)動因子的探測結(jié)果見表 3。由表 3 可知，局域視角下各因子對黃河流域上、中、下游對采礦業(yè)集聚的驅(qū)動作用存在顯著的差異。上游采礦業(yè)集聚的核心驅(qū)動因素為經(jīng)濟(jì)社會發(fā)展水平、政府干預(yù)、勞動力成本、礦產(chǎn)資源稟賦、環(huán)境規(guī)制和勞動力素質(zhì)。上游礦產(chǎn)資源稟賦較好而經(jīng)濟(jì)欠發(fā)達(dá)的地區(qū)，農(nóng)牧業(yè)占比較高，當(dāng)?shù)卣Ｍㄟ^增大采礦業(yè)固定資產(chǎn)投資發(fā)展資源型經(jīng)濟(jì)來帶動區(qū)域經(jīng)濟(jì)發(fā)展。與此同時，較低的勞動力成本和受教育程度，為這些地區(qū)采礦業(yè)的發(fā)展提供了豐富的低成本勞動力，加之較低環(huán)境規(guī)制強(qiáng)度，促進(jìn)了采礦業(yè)在這些地區(qū)的集聚。

　　中游采礦業(yè)集聚的核心驅(qū)動因素為礦產(chǎn)資源稟賦、水資源、環(huán)境規(guī)制、下游產(chǎn)業(yè)需求、勞動力成本、經(jīng)濟(jì)發(fā)展水平和勞動力素質(zhì)。黃河中游地區(qū)礦產(chǎn)資源豐富，但處于干旱區(qū)與半干旱區(qū)，水資源匱乏，生態(tài)本底脆弱。與上游和下游地區(qū)相比，除礦產(chǎn)資源稟賦和勞動力成本等因素外，水資源約束和環(huán)境規(guī)制對黃河中游采礦業(yè)的集聚表現(xiàn)更強(qiáng)的驅(qū)動作用。

　　下游采礦業(yè)集聚的核心驅(qū)動因素為經(jīng)濟(jì)發(fā)展水平、下游產(chǎn)業(yè)需求、礦產(chǎn)資源稟賦和勞動力成本。黃河下游地區(qū)經(jīng)濟(jì)發(fā)展水平較高，資源型地區(qū)產(chǎn)業(yè)轉(zhuǎn)型不斷加快，采礦業(yè)在地區(qū)經(jīng)濟(jì)發(fā)展中的比重日益減少。與此同時，下游地區(qū)靠近資源型產(chǎn)品的消費(fèi)市場，能源消費(fèi)的增加和裝備制造等下游關(guān)聯(lián)產(chǎn)業(yè)的快速發(fā)展，強(qiáng)化了采礦業(yè)的集聚。在考慮以上兩個因素的基礎(chǔ)上，礦產(chǎn)資源稟賦較好，勞動力生產(chǎn)成本較低的下游地區(qū)采礦業(yè)集聚程度更高。

　　4 結(jié)論與政策建議

　　4. 1 結(jié)論

　　1) 黃河流域采礦業(yè)空間基尼系數(shù)呈 “N”型變化趨勢，采礦業(yè)集聚程度總體呈先上升、后波動下降、再快速上升的發(fā)展態(tài)勢。1998 年和 2016 年兩輪去產(chǎn)能政策實施后黃河流域采礦業(yè)集聚水平顯著提升，在 2003—2012 年資源市場的繁榮期黃河流域采礦業(yè)集聚水平出現(xiàn)一定程度的下降，期間 2008 年全球金融危機(jī)引起了短期波動。

　　2) 黃河流域采礦業(yè)不同等級區(qū)位商的空間分布穩(wěn)中有變，中上游地區(qū)的采礦業(yè)集聚水平高于下游地區(qū)，流域內(nèi)采礦業(yè)空間相關(guān)性顯著且關(guān)聯(lián)性逐步增強(qiáng)。“高-高”類型城市集中分布在山西和陜西兩省，空間范圍不斷擴(kuò)大; “低-低”類型的城市空間格局中下游相對穩(wěn)定，中上游逐步減少; “低-高” 類型的城市由忻州市和安陽市變?yōu)檫\(yùn)城市; “高-低”類型的城市數(shù)量最少，僅平頂山和東營在 2018 年表現(xiàn)為此類特征。

　　3) 研究期內(nèi)黃河流域采礦業(yè)產(chǎn)業(yè)重心均落入山西省境內(nèi)，產(chǎn)業(yè)重心總體往西北方向遷移。黃河流域采礦業(yè)空間布局方向性明顯，方 向 角 范 圍 為 89. 24° ～94. 73°，基本保持東—西方向。主軸方向上呈現(xiàn)出持續(xù)收縮—快速擴(kuò)張—逐步收縮的階段性特征，輔軸逐漸減小。長期來看主軸和輔軸均呈收縮趨勢，說明黃河流域采礦業(yè)在空間上更加集聚。

　　4) 各探測因子對黃河流域采礦業(yè)集聚均有一定的驅(qū)動作用，但各區(qū)段探測結(jié)果存在顯著差異。全流域的核心驅(qū)動因素是礦產(chǎn)資源稟賦、政府干預(yù)、經(jīng)濟(jì)發(fā)展水平和下游產(chǎn)業(yè)需求，上游的核心驅(qū)動因素為經(jīng)濟(jì)社會發(fā)展水平、政府干預(yù)、勞動力成本、礦產(chǎn)資源稟賦、勞動力素質(zhì)，中游的核心驅(qū)動因素為礦產(chǎn)資源稟賦、水資源、環(huán)境規(guī)制、下游產(chǎn)業(yè)需求、勞動力成本、經(jīng)濟(jì)發(fā)展水平、勞動力素質(zhì)，下游的核心驅(qū)動因素為經(jīng)濟(jì)發(fā)展水平、下游產(chǎn)業(yè)需求、礦產(chǎn)資源稟賦、勞動力成本。

　　4. 2 政策建議

　　1) 生態(tài)優(yōu)先，有序開發(fā)礦產(chǎn)資源。黃河流域礦產(chǎn)資源豐富，隨著國家大型煤炭基地建設(shè)、“西氣東輸”、“西電東送”等重大工程的實施以及煤炭、鋼鐵、水泥等行業(yè)供給側(cè)結(jié)構(gòu)性改革的深入推進(jìn)，黃河流域特別是中上游地區(qū)采礦業(yè)集聚程度日益提高。但該區(qū)域水資源短缺，生態(tài)本底脆弱，礦產(chǎn)資源開發(fā)與資源環(huán)境特別是水資源的矛盾日益突出。黃河流域未來應(yīng)堅持環(huán)境保護(hù)優(yōu)先，綜合考慮地區(qū)資源環(huán)境承載能力、礦產(chǎn)資源稟賦、開發(fā)利用條件和資源對國家安全的重要程度，有序開發(fā)礦產(chǎn)資源。

　　2) 因地制宜，推動礦業(yè)集聚區(qū)經(jīng)濟(jì)高質(zhì)量轉(zhuǎn)型發(fā)展。采礦業(yè)高度集聚容易形成路徑依賴，經(jīng)濟(jì)發(fā)展對資源型產(chǎn)業(yè)過度依賴可能導(dǎo)致 “資源詛咒”，一旦資源枯竭或資源市場進(jìn)入蕭條期，礦業(yè)集聚區(qū)經(jīng)濟(jì)社會發(fā)展將會出現(xiàn)停滯甚至倒退。黃河流域礦業(yè)集 聚 區(qū)，特 別 是 礦 業(yè) “高 - 高” 集 聚 區(qū)，應(yīng) 以 “云大物智移”技術(shù)的興起為契機(jī)，推動智能礦山建設(shè)，升級傳統(tǒng)主導(dǎo)產(chǎn)業(yè)，培育戰(zhàn)略性新興產(chǎn)業(yè)，提高經(jīng)濟(jì)韌性，規(guī)避 “資源詛咒”。

　　3) 分級分類，加快推進(jìn)退出礦井( 區(qū)) 生態(tài)修復(fù)與資源開發(fā)利用。黃河流域采礦業(yè)集聚時空演化的背后蘊(yùn)含著大量新礦井的建設(shè)和老舊礦井的退出。隨著資源的枯竭和去產(chǎn)能政策的實施，黃河流域有大量礦井( 區(qū)) 關(guān)閉退出。廢棄礦井( 區(qū)) 仍賦存一定程度的可利用資源，如工業(yè)廠房、地下空間、礦井水、旅游資源等。應(yīng)盡快建立退出關(guān)閉礦區(qū)可利用資源潛力評價和利用標(biāo)準(zhǔn)，按照典型示范、分類指導(dǎo)、分級治理、逐步推進(jìn)的原則，兼顧生態(tài)、社會與經(jīng)濟(jì)效益，加快推進(jìn)退出礦區(qū)生態(tài)修復(fù)與資源開發(fā)利用。