基于尺度差異的山西省采煤土地損毀影響因素分析

曹銀貴，白中科，周 偉，王金滿，潘 健，胡 飛，景 明

（1.中國地質大學 （北京）土地科學技術學院，北京100083；2.國土資源部土地整治重點實驗室，北京100035；3.湖北省安陸市國土資源局，湖北 孝感432600）

采煤過程中土地損毀類型主要包括挖損、塌陷和壓占土地[1－3]。井工開采對土地的損毀以塌陷和矸石壓占為主，而露天開采則以挖損和排土場壓占為主[4]。山西省是中國的能源重化工基地，煤炭開采導致土地嚴重損毀。山西省1949～2010年每萬噸煤平均塌陷土地面積為0.15hm2[5]，塌陷土地中耕地的塌陷損毀現象嚴重[6]，其中塌陷損毀的土地中40%為耕地[7]。由于采煤影響，大量水澆地變成旱地，并且大量耕地絕產被棄耕[8－9]。根據環境經濟學理論，礦產資源的開發利用引起的生態環境問題，屬于外部不經濟性問題，應通過內部化來解決[10]。從土地損毀的角度來看，即要探尋煤炭開采過程中影響土地損毀的主要因素及其作用效應。

從宏觀經濟學的角度來看，經濟增長促進能源消耗，通過經濟結構調整，提高經濟運行效率，可以促進能源節約[11]。趙夢楠等[11]采用格蘭杰因果關系檢驗法說明經濟增長是拉動煤炭消耗的原因。從技術經濟學的角度來看，煤炭生產的技術效率受到居多因素的影響[12]。劉新梅[12]運用隨機前沿分析模型驗證煤炭技術效率的區域差異較大，并且地質條件對煤炭開采技術效率有負面影響，而運輸條件、人力資源、開采技術等因素對區域煤炭開采技術效率有正面影響，其中人力資源和開采技術影響較大。土地損毀最直接的影響因素是煤炭資源的開采規模，這是土地損毀的直接驅動力。張金鎖等認為煤炭開采規模的影響因素主要有煤炭儲量、累計開采量、開采率、價格、成本、技術水平等[10]。李振林等[13]認為煤炭資源開發利用總量受到國家煤炭開發總量調控、國民經濟發展、開采技術條件、礦井建設程序及建設時間、運輸條件等因素的影響。

針對煤炭資源最優開采規模的研究從20世紀70年代開始，在研究方法和影響因子選擇層面取得一些重要進展，建立了動態規劃模型[14]、一般競爭性模型[15]、礦產資源優化動態模型[16]、礦產資源開發價值動態決策模型[17]、礦產資源跨期優化模型[16]等。引入了煤炭資源開采的破壞效應影響因素[14]，如技術進步對資源稀缺性的測度[15]、煤炭資源開采的生態約束[16]、技術不確定性的影響[18]、煤炭資源儲量的影響[16]、煤炭資源開發價值的影響[17]、煤炭開采成本的影響[16]等。

上述研究成果主要是分析煤炭資源開采規模的影響因素與研究方法，這些影響因素從類型上可分為正相關性因素和負相關性因素，從影響因素內涵上可分為技術進步型因素、生態約束型因素、開發價值型因素、生產成本型因素、煤炭資源型因素。本次研究以煤炭資源開采規模為紐帶，從山西省尺度、煤炭基地尺度和縣（市、區）煤炭產地尺度探尋土地損毀面積變化的影響因素及其作用效應與差異。

1 研究區概況

山西省重點煤炭基地包括晉北煤炭基地、晉中煤炭基地和晉東煤炭基地。晉北煤炭基地地處山西省太原以北地區，區域面積占山西省總面積的1／3以上，區域總面積4.29萬km2。從行政區劃來看包含大同市、朔州市、忻州市、太原市、呂梁市所屬的縣（市、區）。從煤炭資源分布來看包含大同煤田、寧武煤田全部，河東煤田北部地區，劃分為大同、平朔、朔南、軒崗、河保偏和嵐縣6個礦區。晉中煤炭基地位于山西省中部地區，區域總面積3.81萬km2。從行政區劃來看包含太原市、呂梁市、晉中市、臨汾市、長治市所屬的縣（市、區），從煤炭資源分布來看包含西山煤田、河東煤田（中、南部分）、霍西煤田和沁水煤田（西北部分），劃分為西山、東山、離柳、汾西、霍州、鄉寧、霍東、石隰8個礦區。晉東煤炭基地地處山西省東南部地區，區域總面積2.81萬km2。從行政區劃來看包含陽泉市、晉中市、長治市、晉城市和臨汾市所屬的縣（市、區），從煤炭資源分布來看包含沁水煤田東南部，劃分為陽泉、武夏、潞安和晉城4個礦區。從煤炭開采的工藝分布來看，2000～2010年期間，山西省井工開采的煤炭產量占全省煤炭總產量的92%，露天開采的煤炭產量占全省煤炭總產量的8%[20]。

2 數據來源與研究方法

2.1 數據來源

本次研究所涉及的煤炭產量、耕地面積、地區生產總值、煤炭工業投資等數據主要來源于2001～2011年的《山西省統計年鑒》，煤炭基礎儲量數據來源于2005～2011年的《全國礦產資源儲量通報》。

2.2 研究方法

2.2.1 土地損毀系數確定的方法

土地損毀系數是指在一段時間范圍內每采萬噸煤所導致的土地損毀面積[21]。從煤炭開采對土地的損毀來看，土地損毀系數與地貌類型、采礦工藝和損毀土地的類型緊密相關[22]。從目前已有的研究來看，土地損毀系數主要從2個方面來確定。

1）基于地貌類型的土地損毀系數確定。基于山西省井工煤礦土地復墾方案統計分析，山區萬噸煤土地損毀系數的范圍為0.09～0.34hm2，平均萬噸煤土地損毀系數為0.25hm2；丘陵區萬噸煤土地損毀系數的范圍是0.02～0.57hm2，平均萬噸煤土地損毀系數為0.23hm2；平原區萬噸煤土地損毀系數為0.19～0.43hm2，平均萬噸煤土地損毀系數為0.31hm2。經地貌類型加權求和分析，山西省井工煤礦生產萬噸煤土地損毀系數平均為0.23hm2。基于山西省露天煤礦土地復墾方案統計分析，丘陵區萬噸煤土地損毀系數的范圍為0.04～0.55hm2，平均萬噸煤土地損毀系數為0.24hm2；低山侵蝕區萬噸煤土地損毀系數的范圍為0.09～0.39hm2，平均萬噸煤土地損毀系數為0.24hm2；中高山區萬噸煤土地損毀系數的范圍為0.14～0.22hm2，平均萬噸煤土地損毀系數為0.20hm2。經地貌類型加權求和分析，山西省露天煤礦平均萬噸煤土地損毀系數為0.24hm2。

2）基于損毀類型的土地損毀系數確定。由于煤礦區土地損毀的類型主要包括挖損、塌陷和壓占3個類型，從這3個類型來分析土地損毀可以使土地損毀的過程更加清晰。鑒于此，何書金、蘇光全[21]在大量實證研究的基礎上，從土地損毀的3個類型出發，確定了土地損毀系數計算模型。其計算公式分別為式（1）、式（2）、式（3）。

式中：Sw為挖損土地面積，hm2；St為塌陷土地面積，hm2；Sz為壓占土地面積，hm2；O為煤炭產量，萬t；Lr為露天開采率，%；Wr為露天開采萬噸煤產量的挖損土地面積，hm2；Jr為井工開采率，%；Tr為井工開采萬噸煤產量的塌陷土地面積，hm2。Zl為露天開采萬噸煤產量的壓占土地面積，hm2；Zj為井工開采萬噸煤產量的壓占土地面積，hm2。

何書金、蘇光全研究表明，露天開采每萬噸煤挖損土地面積為0.104hm2，井工開采每萬噸煤塌陷土地面積為0.270hm2，露天開采每萬噸煤壓占土地面積為0.107hm2，井工開采每萬噸煤壓占土地面積為0.005hm2[21]。

2.2.2 土地損毀面積與影響因素關系的研究方法

1）散點擬合分析法。從不同尺度出發，利用Excel 2003的散點擬合工具分別構建散點擬合模型來分析土地損毀面積與單一影響因素之間的關系。為提高擬合精度，將土地損毀面積與影響因素取對數（以log10為底的對數）進行擬合分析。為了提高擬合函數的可比性，結合擬合系數R2的大小，統一擬合成二次函數。

2）多元回歸分析法。從不同尺度出發，利用SPSS16.0的多元回歸分析工具分別研究土地損毀面積與多個影響因素之間的關系。根據F檢驗和t檢驗分別評價回歸模型的顯著性和影響因素系數的顯著性。

3 結果與分析

3.1 基于不同損毀系數的土地損毀面積對比

利用基于地貌類型的土地損毀系數和煤炭產量計算得出山西省2000～2010年總土地損毀面積、露天開采土地損毀面積和井工開采土地損毀面積。2000～2010年間，山西省總土地損毀面積逐年遞增，累計增加了1.95倍，其年均增長率為11.43%。但在2009年，由于煤炭產量的下降，總土地損毀面積略微有所下降。在山西省，露天開采的煤炭產量比例較低，相比井工開采而言，其土地損毀面積遠低于井工開采土地損毀面積，但是二者之間土地損毀面積差在縮小。2000年露天開采土地損毀面積211.53hm2，井工開采土地損毀面積5582.24hm2，井工開采土地損毀面積是露天開采土地損毀面積的26.39倍，2010年，二者之間的倍數變為11.03倍。

利用基于損毀類型的土地損毀系數和煤炭產量計算得出山西省2000～2010年總土地損毀面積、露天開采土地損毀面積和井工開采土地損毀面積。2000～2010年間，山西省總土地損毀面積在逐年遞增，累計增加了1.91倍，其年均增長率為11.29%。露天開采土地損毀面積和井工開采土地損毀面積也在逐年遞增。從具體的損毀類型來看，挖損土地面積最少，其次是壓占土地面積，塌陷土地面積最大。2000年挖損土地面積91.66hm2，壓占土地面積是挖損土地面積的2.35倍，塌陷土地面積是挖損土地面積的71.49倍。2010年挖損土地面積615.85hm2，壓占土地面積是挖損土地面積的1.58倍，塌陷土地面積是挖損土地面積的29.89倍。從各類型的年均增長率來看，挖損土地面積的增長率最快，挖損土地面積年均增長率為20.98%，塌陷土地年均增長率為10.88%，壓占土地面積年均增長率為16.28%。

通過對比土地損毀面積差發現，基于地貌類型的土地損毀系數計算所得的總土地損毀面積和井工開采土地損毀面積要小于基于損毀類型的土地損毀系數計算所得的總土地損毀面積和井工開采土地損毀面積，而露天開采土地損毀面積則反之。為了提高土地損毀面積的準確度，將2種損毀系數計算所得的土地損毀面積求和平均，得到平均土地損毀面積、露天開采平均土地損毀面積和井工開采平均土地損毀面積。

利用平均土地損毀面積除以總的煤炭產量、露天開采平均土地損毀面積除以露天開采煤炭產量、井工開采平均土地損毀面積除以井工開采煤炭產量分別得出相應的平均土地損毀系數。2000～2010年期間山西省采煤土地損毀系數范圍為0.2503～0.2515hm2／萬t煤，其平均土地損毀系數為0.2509 hm2／萬t煤。露天開采平均土地損毀系數為0.2255hm2／萬t煤，井工開采平均土地損毀系數為0.2525hm2／萬t煤。利用平均土地損毀系數計算土地損毀面積，并探討土地損毀面積與影響因子之間的關系。

3.2 土地損毀面積影響因素分析

煤炭資源基礎儲量屬于煤炭資源型因素，反映了一定經濟技術條件下能夠開采的煤炭資源總量。地區生產總值屬于開發價值型因素，在山西省，煤炭工業產值是地區生產總值的重要組成部分，其煤炭資源開采的價值量一定程度上會影響地區生產總值。煤炭工業投資屬于生產成本型因素，增加對煤炭工業的投資，改善煤炭開采工藝，提高煤炭生產的技術水平，勢必會減少煤炭開采過程中的土地損毀。耕地面積屬于生態約束型因素，耕地保護是中國的基本國策，必須得嚴守1.5億hm2耕地紅線。山西省煤炭產地近40%屬于煤糧復合區，煤炭開采對耕地的損毀非常嚴重。保護耕地，防止耕地在煤炭開采過程中受損有利于礦區生態環境保護。

在山西省尺度，以2000～2010年為研究期，選擇了煤炭資源基礎儲量、地區生產總值、煤炭工業投資和耕地面積4個因素探尋與土地損毀面積之間的關系。在煤炭基地尺度和縣（市、區）煤炭產地尺度，以2000～2008年為研究期，選擇了地區生產總值和耕地面積2個因素探尋與土地損毀面積之間的關系。

3.2.1 山西省尺度

1）土地損毀面積與煤炭基礎儲量的關系。在2000～2010年內，山西省土地損毀面積與煤炭基礎儲量對數之間的擬合度高達0.8869（圖1（a））。通過對擬合函數的最值分析，當煤炭基礎儲量對數達到2.92億t時，即煤炭基礎儲量達到839.40億t時，土地損毀面積達到最大值1.88萬hm2。由于該擬合函數曲線開口向下（圖1（a）），并且曲線位于擬合函數對稱軸的右側，表明土地損毀面積隨著煤炭基礎儲量對數的增加而減少，即當煤炭基礎儲量大于839.40億t時，土地損毀面積隨煤炭基礎儲量的增加而減少。

2）土地損毀面積與地區生產總值的關系。在2000～2010年內，山西省土地損毀面積與地區生產總值對數之間的擬合度高達0.9747（圖1（b））。通過對擬合函數的最值分析，當地區生產總值對數為8.18萬元時，即地區生產總值為15241.80億元時，土地損毀面積達到最大值1.81萬hm2。由于該擬合函數曲線開口向下（圖1（b）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積隨著地區生產總值對數的增加而增加，即當地區生產總值小于15241.80億元時，土地損毀面積隨地區生產總值的增加而增加。

3）土地損毀面積與煤炭工業投資的關系。在2000～2010年內，山西省土地損毀面積與煤炭工業投資對數之間的擬合度高達0.9729（圖1（c））。通過對擬合函數的最值分析，當煤炭工業投資對數為7.88萬元時，即煤炭工業投資為7593.98億元時，土地損毀面積達到最大值1.98萬hm2。由于該擬合函數曲線開口向下（圖1（c）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積隨著煤炭工業投資對數的增加而增加，即當煤炭工業投資小于7593.98億元時，土地損毀面積隨煤炭工業投資的增加而增加。

4）土地損毀面積與耕地面積的關系。在2000～2010年內，山西省土地損毀面積與耕地面積對數之間的擬合度高達0.8263（圖1（d））。通過對擬合函數的最值分析，當耕地面積對數為2.54萬hm2時，即耕地面積為350.39萬hm2時，土地損毀面積達到最大值1.57萬hm2。由于該擬合函數曲線開口向下（圖1（d）），并且曲線位于擬合函數對稱軸的右側，表明土地損毀面積隨著耕地面積對數的增加而減少，即當耕地面積大于350.30萬hm2時，土地損毀面積隨耕地面積的增加而減少。

在研究期內就單一因素來看，山西省土地損毀面積隨著煤炭基礎儲量、耕地面積的增加而減少，而隨著地區生產總值、煤炭工業投資的增加而增加。但是煤炭開采過程中土地損毀面積的變化是多個因素綜合作用的結果，因此需要從多因素共同作用層面來分析土地損毀面積的變化。

5）土地損毀面積多因素影響分析。采用多元回歸分析法，建立了土地損毀面積與煤炭基礎儲量、地區生產總值、煤炭工業投資及耕地面積之間的多元回歸模型。多元回歸模型為y＝2.97－1.43x1＋0.66x2＋0.26x3－0.94x4，其回歸關系擬合程度高達0.9940，并且回歸關系能通過F檢驗的要求，Sig＜0.001。另外標準化系數均能達到t檢驗的要求，Sig＜0.05。可見，多元回歸模型模擬精度較高。并且在多因素綜合作用的結果下，煤炭基礎儲量、耕地面積是影響土地損毀面積的負相關因素，而地區生產總值、煤炭工業投資是影響土地損毀面積的正相關因素。

圖1 2000～2010年山西省土地損毀面積與各影響因素之間的關系、

3.2.2 煤炭基地尺度

1）土地損毀面積與地區生產總值的關系。晉北煤炭基地在2000～2008年內，土地損毀面積對數與地區生產總值對數擬合度高達0.9449（圖2（a））。通過對擬合函數的最值分析，當地區生產總值對數為6.69萬元時，即地區生產總值為490.23億元時，土地損毀面積對數達最大值3.77hm2，即土地損毀面積達到最大值5886.91hm2。由于該擬合函數曲線開口向下（圖2（a）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積對數隨地區生產總值對數的增加而增加，即當地區生產總值小于490.23億元時，土地損毀面積隨地區生產總值的增加而增加。

晉中煤炭基地在2000～2008年內，土地損毀面積對數與地區生產總值對數擬合度高達0.9573（圖2（b））。通過對擬合函數的最值分析，當地區生產總值對數為7.52萬元時，即地區生產總值為3283.66億元時，土地損毀面積對數達到最大值3.57hm2，即土地損毀面積達到最大值3751.43hm2。由于該擬合函數曲線開口向下（圖2（b）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積對數隨地區生產總值對數的增加而增加，即當地區生產總值小于3283.66億元時，土地損毀面積隨地區生產總值的增加而增加。

晉東煤炭基地在2000～2008年內，土地損毀面積對數與地區生產總值對數擬合度高達0.9771（圖2（c））。通過對擬合函數的最值分析，當地區生產總值對數為7.24萬元時，即地區生產總值為1755.92億元時，土地損毀面積對數達到最大值3.68hm2，即土地損毀面積達到最大值4831.45hm2。由于該擬合函數曲線開口向下（圖2（c）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積對數隨地區生產總值對數的增加而增加，即當地區生產總值小于1755.92億元時，土地損毀面積隨地區生產總值的增加而增加。

2）土地損毀面積與耕地面積的關系。晉北煤炭基地在2000～2008年內，損毀土地面積對數與耕地面積對數擬合度高達0.8945（圖2（d））。通過對擬合函數的最值分析，當耕地面積對數為1.92萬hm2時，即耕地面積為83.64萬hm2時，土地損毀面積對數達到最大值3.74hm2，即土地損毀面積達到最大值5479.63hm2。由于該擬合函數曲線開口向下（圖2（d）），并且曲線位于擬合函數對稱軸的右側，表明土地損毀面積對數隨耕地面積對數的增加而減少，即當耕地面積大于83.64萬hm2時，土地損毀面積隨耕地面積的增加而減少。

晉中煤炭基地在2000～2008年內，土地損毀面積對數與耕地面積對數擬合度高達0.8282（圖2（e））。通過對擬合函數的最值分析，當耕地面積對數為1.99萬hm2時，即耕地面積為99.38萬hm2時，土地損毀面積對數達到最小值3.11hm2，即土地損毀面積達到最小值1293.51hm2。由于該擬合函數曲線開口向上（圖2（e）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積對數隨耕地面積對數的增加而減少，即當耕地面積小于99.38萬hm2時，土地損毀面積隨耕地面積的增加而減少。

晉東煤炭基地在2000～2008年內，土地損毀面積對數與耕地面積對數擬合度高達0.9216（圖2（f））。通過對擬合函數的最值分析，當耕地面積對數為1.76萬hm2時，即耕地面積為57.61萬hm2時，土地損毀面積對數達到最小值3.05hm2，即土地損毀面積達到最小值1109.93hm2。由于該擬合函數曲線開口向上（圖2（f）），并且曲線位于擬合函數對稱軸的左側，表明土地損毀面積對數隨耕地面積對數的增加而減少，即當耕地面積小于57.61萬hm2時，土地損毀面積隨耕地面積的增加而減少。

總的來看，在3大煤炭基地，土地損毀面積均隨地區生產總值的增加而增加，隨耕地面積的增加而減少。在各煤炭基地，地區生產總值與耕地面積對土地損毀面積的共同影響需要進一步分析。

3）土地損毀面積多因素影響分析。在3大煤炭基地，采用多元回歸分析法，分別建立了土地損毀面積與地區生產總值、耕地面積之間的多元回歸模型。晉北基地、晉中基地、晉東基地的多元回歸模型分別是：y＝4.52＋0.18x2－3.70x4，y＝2.75＋0.39x2－2.14x4，y＝3.71＋0.31x2－3.01x4，3大煤炭基地土地損毀面積與影響因素之間回歸關系擬合程度高，其擬合系數分別是0.9930、0.9320、0.9950，并且回歸關系均能通過F檢驗的要求，Sig＜0.001。另外標準化系數均能達到t檢驗的要求，Sig＜0.05。可見，多元回歸模型模擬精度較高。并且在2個因素共同作用的情況下，3大煤炭基地都能體現出地區生產總值與土地損毀面積之間的正相關性，耕地面積與土地損毀面積之間的負相關性。

圖2 2000～2008年3大煤炭基地土地損毀面積與各影響因素之間的關系

3.2.3 縣（市、區）煤炭產地尺度

1）土地損毀面積單因素關系分析。在2000～2008年期間，相比山西尺度和煤炭基地尺度，在縣（市、區）煤炭產地尺度土地損毀面積對數值與地區生產總值對數值、耕地面積對數值之間的擬合程度不高。地區生產總值對數與土地損毀面積對數的擬合曲線發展趨勢相對一致（圖3（a）、圖3（b）、圖3（c）），土地損毀面積對數隨地區生產總值對數的增加而增加。耕地面積對數與土地損毀面積對數的擬合曲線發展趨勢一致性差，耕地面積對數與土地損毀面積對數之間的相關性體現不明顯（圖3（d）、圖3（e）、圖3（f））。

2）土地損毀面積多因素影響分析。在3大煤炭基地內，分別建立3大煤炭基地各縣（市、區）煤炭產地土地損毀面積與地區生產總值、耕地面積之間的多元回歸模型。晉北基地煤炭產地、晉中基地煤炭產地、晉東基地煤炭產地的多元回歸模型分別是：y＝－0.03＋0.53x2＋0.50x4，y＝1.06＋0.11x2－0.18x4，y＝0.35＋0.34x2＋0.31x4。回歸關系雖均能通過F檢驗的要求，Sig＜0.001，但是回歸關系擬合程度不高，擬和系數均低于0.5。可見，多元回歸模型模擬精度低。另外，在2個因素共同作用的情況下，3大煤炭基地各縣（市、區）煤炭產地都能體現出地區生產總值與土地損毀面積之間的正相關性，但耕地面積與土地損毀面積之間的負相關性體現不明顯，其中晉北煤炭基地各縣（市、區）煤炭產地和晉東煤炭基地各縣（市、區）煤炭產地表現出了正相關性，而晉中煤炭基地各縣（市、區）煤炭產地表現出了負相關性。

圖3 2000～2008年山西省3大煤炭基地各縣（市、區）煤炭產地土地損毀面積與各影響因素之間的關系

4 結語

1）在研究期內，山西省在煤炭開采過程中土地損毀具有明顯的階段性和動態性。由于土地損毀具有階段性和動態性，精確的土地損毀系數適合于單個礦區的研究，在省域尺度、煤炭基地尺度及縣（市、區）煤炭產地尺度很難確定精確的土地損毀系數，因此需要引入平均土地損毀系數研究大尺度層面土地損毀情況。本文基于地貌類型土地損毀系數和基于損毀類型土地損毀類型系數計算了山西省2000～2010年內的平均土地損毀系數，其值為0.2509hm2／萬t煤。露天開采平均土地損毀系數為0.2255hm2／萬t煤，井工開采平均土地損毀系數為0.2525hm2／萬t煤。另外，確定平均土地損毀系數的研究期不宜太長。在煤炭開采的初期，每采萬噸煤的土地損毀面積（土地損毀系數）相對比較大，包括工業場地占地、排土場及矸石山的下墊面占地及首采煤層形成的塌陷。隨著開采過程的深入，工業場地、排土場、矸石山的下墊面已經相對穩定，并且采坑形成了內排土場，這樣土地損毀系數會逐漸減小。另外井工礦隨所采煤層的深度變化，一定程度上會加深塌陷的深度，新的土地損毀面積會有所減小，因此也可以降低土地損毀系數。

2）在煤炭開采過程中，土地損毀面積的影響因素比較多，包括新的資源探明、技術革新、煤炭賦存條件、沉積層的松散程度等等。從統計學的角度來看，這些數據之間肯定存在一定的關系，因此有必要在研究過程中進一步明晰彼此之間的真實聯系。就一個尺度情況下，影響因素與土地損毀面積表現出的關系需要進一步分析。因此需要從多個尺度出發，利用單因素擬合分析和多因素回歸分析，在相關檢驗通過的情況下，進一步分析影響因素與土地損毀面積之間表現出的某種聯系的趨同性與差異性。影響因素對土地損毀面積作用的過程是復雜的，綜合的，既有趨同性，又有差異性。在山西省尺度，煤炭基礎儲量、耕地面積是影響土地損毀面積的負相關性因素，而地區生產總值、煤炭工業投資是影響土地損毀面積的正相關性因素。在煤炭基地尺度，體現出地區生產總值對土地損毀面積的正相關性，耕地面積對土地損毀面積的負相關性。在縣（市、區）煤炭產地尺度，基本上能體現出地區生產總值對土地損毀面積的正相關性，但耕地面積對土地損毀面積的負相關性體現不明顯。山西省尺度和煤炭基地尺度是利用地區生產總值和耕地面積控制采煤土地損毀面積的最佳尺度。

3）在山西省尺度，煤炭基礎儲量在研究期內與土地損毀面積之間是負相關的，由于山西省煤炭開采時間長，煤炭基礎儲量上升空間有限。期初煤炭基礎儲量高，開采相對容易，造成的土地損毀面積相對較小，隨著煤炭開采速度加快，在采煤技術水平提升有限的情況下，煤炭基礎儲量下降，采煤的難度逐漸加大，勢必加大土地損毀。地區生產總值促進土地損毀。由于山西省經濟產業結構相對單一，經濟發展對煤炭工業的依賴程度高，在追求經濟發展的過程中，煤炭的產量逐步提升，勢必影響煤炭開采的規模，導致土地損毀面積增加。在煤炭基地尺度和縣（市、區）煤炭產地尺度均能體現地區生產總值對土地損毀面積的正相關性。在山西省尺度，煤炭工業投資促進了土地損毀，可見投資的方向主要在煤炭開采層面。因此，要合理調配煤炭工業投資的方向，一方面控制生產性的投資，另一方面加大對技術革新及生態環境保護方面的投資，減輕煤炭開采對土地的損毀。在山西省尺度和煤炭基地尺度耕地面積對土地損毀面積的負相關性是明顯的，而在縣（市、區）煤炭產地尺度體現的不是很明顯。相對山西省尺度和煤炭基地尺度而言，各縣（市、區）煤炭產地尺度耕地面積及土地損毀面積偏低，負相關性不明顯，在大統計樣本情況下很難建立二者之間的關系。總之，針對采煤土地損毀面積的控制應該從不同的尺度出發，制定針對性的控制措施。

致謝：感謝陳思、孫琦、盧元清、宋雪姣等碩士研究生提供的數據支持！

[1]陜永杰，郝蓉，白中科，等.礦區復合生態系統中土壤演替和植被演替的相互影響[J].煤礦環境保護，2001，15（5）：28－30.

[2]胡振琪.半干旱地區煤殲石山綠化技術研究[J].煤炭學報，1995，20（3）：322－324.

[3]何書金，蘇光全.礦區廢棄土地復墾潛力評價方法與應用實例[J].地理研究，2000，19（2）：163－171.

[4]白中科，郭青霞，史源英，等.安家嶺露天煤礦土地利用結構預測[J].煤炭學報，1999，24（2）：207－211.

[5]李枝榮.立足存量挖潛，集約利用土地，實現耕地總量的動態平衡[J].改革先聲，1997，（2）：21－22.

[6]高慧卿，高敏.山西省耕地資源評價與可持續利用對策[J].山西農業科學，2011，39（6）：595－598.

[7]李素清，張金屯.山西生態環境破壞對可持續發展的影響及對策研究[J].干旱區資源與環境，2005，19（2）：56－61.

[8]李錚，王晉民.山西省土肥水資源利用問題剖析及應對糧食安全的對策[J].山西農業科學，2005，33（4）：3－6.

[9]王國梁.建國以來山西省耕地變化態勢與耕地安全[J].經濟地理，2010，30（9）：1542－1545.

[10]張金鎖，王喜蓮.煤炭開采規模的影響因素及模型研究[J].能源技術與管理，2007，（5）：68－71.

[11]趙夢楠，周德群.中國煤炭消耗與經濟增長的結構變化及因果關系研究[J].管理學報，2008，5（5）：717－721.

[12]劉新梅，許中偉，徐潤芳.基于SFA模型的我國區域煤炭生產技術效率及其影響因素的實證研究[J].軟科學，2008，22（3）：37－40.

[13]李振林，馮耀挺，贠培琪.咸陽市煤炭資源開發利用總量調控影響因素分析[J].中國煤炭，2010，36（4）：30－33.

[14]Schulze W D.The optimal use of non－renewable resources：The theory of extraction[J].Journal of Environmental Economics and Management，1974，1（1）：53－73.

[15]Farzin Y H.Technological change and the dynamics of resource scarcity measures[J].Journal of Environmental Economics and Management，1995，29（1）：105－120.

[16]張金鎖.煤礦礦井技術改造方案專家決策系統研究[J].西安科技學院學報，2000，20（3）：196－200.

[17]魏曉平，王新宇.礦產資源最適耗竭經濟分析[J].中國管理科學，2002，10（5）：79－82.

[18]Anne E，Aude P.Optimally eating a stochastic cake：A recursive utility approach [J].Resource and Energy Economics，2003，25（2）：129－139.

[19]蔡登谷.關于山西礦區復墾的考察報告[J].林業經濟，2008，（4）：36－38.

[20]山西省統計局編.山西統計年鑒[M].北京：中國統計出版社，2011.

[21]何書金，蘇光全.中國采礦業的發展與礦區土地損毀預測[J].資源科學，2002，24（2）：17－21.

[22]周偉，曹銀貴，白中科，等.礦區土地復墾監測指標體系[J].中國土地科學，2012，26（11），68－73.