流域“地-礦”土地水資源利用沖突測度確定及土地整治策略

劉慧芳，畢如田，文 博

（1. 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，太谷 030801；2. 南京農(nóng)業(yè)大學(xué)公共管理學(xué)院，南京 210095）

·土地整理工程·

流域“地-礦”土地水資源利用沖突測度確定及土地整治策略

劉慧芳1,2，畢如田1※，文 博2

（1. 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，太谷 030801；2. 南京農(nóng)業(yè)大學(xué)公共管理學(xué)院，南京 210095）

量化區(qū)域煤炭開采與水土資源相互影響程度，有助于把握土地利用過程中的主要限制因素，提高土地整治的針對性與有效性。該文以山西省澤州縣長河流域?yàn)檠芯繀^(qū)，提出煤炭開采與水土資源利用相互影響測度與土地整治策略的思路和方法。首先提出煤炭開采與水土資源相互影響分析框架，并將土地、水和煤炭作為主要的相互影響對象，將土地利用系統(tǒng)作為影響載體，借助數(shù)學(xué)模型和GIS空間分析技術(shù)，依次評價(jià)土地?fù)p毀累積作用指數(shù)、水資源匱乏指數(shù)和煤炭資源損失指數(shù)；然后確定綜合影響值及其等級空間分布，并據(jù)此將流域劃分為重點(diǎn)生態(tài)修復(fù)區(qū)、重點(diǎn)生態(tài)保護(hù)區(qū)、損毀土地復(fù)墾區(qū)、耕地資源保護(hù)區(qū)、水資源匱乏修復(fù)區(qū)、水資源與耕地資源保護(hù)區(qū)和耕地質(zhì)量提升區(qū)等七大土地利用功能區(qū)，并根據(jù)主要限制因素提出相應(yīng)的土地整治策略。研究結(jié)果表明綜合影響高值區(qū)及中值區(qū)面積占研究區(qū)總面積的分別是20.21%和33.39%，主要分布在流域西部，是煤炭開采活動(dòng)集聚區(qū)，煤炭壓覆及土地?fù)p毀是沖突主要表現(xiàn)特征，劃分為重點(diǎn)生態(tài)修復(fù)區(qū)、重點(diǎn)生態(tài)保護(hù)區(qū)、損毀土地復(fù)墾區(qū)和耕地資源保護(hù)區(qū)；綜合影響低值區(qū)及（基本）無影響區(qū)分別占24.24%和22.16%，主要分布在流域東部，水資源匱乏是沖突主要表現(xiàn)特征，劃分為水資源匱乏修復(fù)區(qū)、水資源與耕地資源保護(hù)區(qū)、耕地質(zhì)量提升區(qū)，以上分區(qū)與土地利用總體規(guī)劃中土地用途分區(qū)基本吻合。綜合影響測度模型有利于定量分析“煤炭-水土”相互利用影響等級與主導(dǎo)限制因素，為制定有效的土地整治策略提供依據(jù)。

土地利用；土地整治；水資源；流域；綜合影響指數(shù)；土地?fù)p毀累積指數(shù)；水資源匱乏指數(shù)；煤炭損失指數(shù)

0 引 言

長期以來，水土資源、煤炭能源分別作為中國社會經(jīng)濟(jì)發(fā)展的基礎(chǔ)物質(zhì)和廉價(jià)動(dòng)力，在國民經(jīng)濟(jì)發(fā)展過程中具有舉足輕重的作用。其中土地資源是社會經(jīng)濟(jì)發(fā)展的重要空間載體，根據(jù)《全國土地利用總體規(guī)劃（2006—2020）》規(guī)劃建設(shè)用地規(guī)模從2005年的3 192.24萬hm2增至2020年3 777.24萬hm2，而2003年以來中國提出土地參與宏觀調(diào)控政策后，更加凸顯土地資源為中國社會經(jīng)濟(jì)發(fā)展提供了空間載體及經(jīng)濟(jì)發(fā)展動(dòng)力。而我國水資源匱乏，人均水資源占有量僅相當(dāng)于世界人均水平的1/4[1]；煤炭是中國重要的生產(chǎn)能源及工業(yè)原料，長期以來一直占中國能源生產(chǎn)和消費(fèi)總量的 65%以上；隨著社會經(jīng)濟(jì)發(fā)展，中國煤炭采量從2000年的10×108t增至2014年的近40×108t，以每年近2億t的速度增加滿足社會經(jīng)濟(jì)增長需求[2]。但隨著煤炭開采及水土資源利用規(guī)模擴(kuò)大，煤炭開采與水土利用相互約束不斷加劇，逐漸成為限制煤炭資源富集區(qū)資源可持續(xù)利用的重要因素。

國外對煤炭開采與水土資源利用相互影響研究多從礦產(chǎn)開采前煤炭開采與土地利用綜合規(guī)劃及管理，確定合理的資源利用方式達(dá)到預(yù)防資源利用沖突的目的。Craynon等[3]基于社會、經(jīng)濟(jì)、環(huán)境3個(gè)方面因素考慮，采用系統(tǒng)工程方法優(yōu)化礦山設(shè)計(jì)，提出運(yùn)用“GIS”技術(shù)來識別存在區(qū)域的煤炭-水-基礎(chǔ)設(shè)施沖突情況，預(yù)先做好礦山規(guī)劃，避免煤炭開采過程中沖突發(fā)生，并以美國Appalachia露天煤礦為例進(jìn)行了驗(yàn)證。Wawryk[4]針對澳大利亞南部保護(hù)區(qū)阿卡魯拉地區(qū)煤炭開采與土地利用相互影響及解決方式進(jìn)行研究，認(rèn)為現(xiàn)存的法律制度具有明顯層次性保護(hù)過程，即具有高價(jià)值的公共土地優(yōu)于私人土地利用，因此，為更好的保護(hù)礦區(qū)范圍內(nèi)的私人土地，提出增加采礦用地壓力，在阿卡魯拉保護(hù)區(qū)建立保護(hù)區(qū)新體系。Hartje等[5]在其研究中指出大部分國家或地區(qū)的采礦者與資源所有者之間簽訂合同時(shí)也規(guī)定了采礦者應(yīng)承擔(dān)的責(zé)任和義務(wù)，并以美國和德國的實(shí)踐為例進(jìn)行了實(shí)證分析，認(rèn)為財(cái)政保證是礦山生態(tài)環(huán)境恢復(fù)治理的重要而有效的策略。國內(nèi)學(xué)者多從資源協(xié)調(diào)利用的角度分析煤炭開采對土地、水資源造成的負(fù)面影響，并著重提出相應(yīng)的調(diào)控措施。婁華君等[6]基于煤、地、水互耗矛盾產(chǎn)生機(jī)理的基礎(chǔ)上，并以南水北調(diào)東線工程為例，提出煤、地、水資源同步利用模型；李鋼[7]認(rèn)為礦地矛盾突出，需通過地礦系統(tǒng)內(nèi)部協(xié)同，更需要通過統(tǒng)籌實(shí)現(xiàn)礦地雙方的共同發(fā)展，主要包括礦地統(tǒng)籌管理、礦地協(xié)調(diào)利用和礦、城、鄉(xiāng)統(tǒng)籌發(fā)展，并提出相應(yīng)技術(shù)體系和統(tǒng)籌方向；陳小芳等[8]從宿州城礦鄉(xiāng)一體化進(jìn)程中發(fā)現(xiàn)城鄉(xiāng)矛盾、城礦鄉(xiāng)之間的沖突即城鎮(zhèn)煤炭壓覆、礦產(chǎn)開采減少城鎮(zhèn)可利用空間等問題，并提出宿州城礦鄉(xiāng)統(tǒng)籌一體化發(fā)展規(guī)劃。目前學(xué)者多從定性角度對煤炭開采與水土資源利用相互影響進(jìn)行研究，但未提出具體的定量評價(jià)方法，缺乏煤炭開采與水土資源利用相互影響程度空間化，以及針對具體類型及程度的有效整治措施。

本文依據(jù)已有研究成果，以山西省澤州縣長河流域?yàn)槔趯?shí)地調(diào)研、資料統(tǒng)計(jì)與GIS技術(shù)，分析土地、水及煤炭資源的分布與利用現(xiàn)狀，對流域煤炭開采與水土資源利用相互影響進(jìn)行定量評價(jià)，并基于此對綜合影響程度分級分區(qū)，實(shí)現(xiàn)流域煤炭開采與水土資源相互影響程度空間差異的可視化；明確分區(qū)中各類資源利用現(xiàn)狀及存在主要問題，更好地識別土地利用中存在的主要限制因素，為進(jìn)一步明確整治對象與生態(tài)修復(fù)提供了空間途徑，為流域生態(tài)恢復(fù)、土地復(fù)墾及土地綜合整治規(guī)劃提供依據(jù)。

1 理論框架

煤炭開采與水土資源相互影響主要指由于土地、水或礦產(chǎn)資源在開發(fā)、利用、保護(hù)和管理過程中，與其他類型資源產(chǎn)生相互不利影響，造成資源互耗現(xiàn)象，具體表現(xiàn)形式主要有土地?fù)p毀、煤炭壓覆、水污染等問題。土地?fù)p毀形成是由于土地和煤層在空間分布上是垂直交錯(cuò)疊置關(guān)系，煤炭開采活動(dòng)在破壞地下巖層空間結(jié)構(gòu)的同時(shí)，對地表造成塌陷、裂縫和壓占的情況，削減土地資源可利用數(shù)量，據(jù)測算，煤炭資源與耕地資源疊置區(qū)域面積占中國耕地總量的40%以上[9]，平均采煤塌陷系數(shù)為0.33 hm2/萬t。另外煤炭開采破壞耕地保水保肥性能，降低耕地質(zhì)量，減弱煤炭采區(qū)的耕地資源空間優(yōu)勢。另外，由于采區(qū)地表不均勻沉降，造成地上建筑物出現(xiàn)裂隙、損毀及倒塌現(xiàn)象，據(jù)山西省抽樣調(diào)查，全省平均每開采1 000萬t煤炭會造成2.05萬m2建筑物損壞[10]；煤炭壓覆問題是農(nóng)村居民點(diǎn)、水庫及城市建設(shè)區(qū)與煤炭開采區(qū)在空間上復(fù)合，導(dǎo)致煤炭資源壓覆，降低煤炭回采率，減弱煤炭資源空間稟賦優(yōu)勢；水污染問題主要包括一方面井工礦或露天礦坑開采過程中抽排大量地下水，導(dǎo)致地下水位下降，無法滿足該區(qū)域生產(chǎn)、生活用水，另一方面煤礦排放的污水、廢水、廢液滲入地下水，污染水質(zhì)。因此，煤炭開采破壞水資源空間分布平衡性及穩(wěn)定性。由于土地、水及礦產(chǎn)資源利用過程中形成互耗，具體體現(xiàn)在資源數(shù)量減少、資源質(zhì)量下降和資源空間配置低效等方面，是煤炭開采與水土資源之間不利影響形成的主要原因（圖 1）。因此在土地綜合整治過程中，應(yīng)該更多地指向資源利用過程中資源數(shù)量與質(zhì)量上的損耗以及資源空間配置中的不協(xié)調(diào)方面。

本文首先以山西省澤州縣長河流域47個(gè)行政村為評價(jià)單元，從資源系統(tǒng)的數(shù)量屬性、質(zhì)量屬性和空間屬性3個(gè)層面構(gòu)建評價(jià)指標(biāo)，采用綜合指數(shù)和法和理想景觀向量模型計(jì)算土地?fù)p毀累積作用指數(shù)、煤炭開采損失程度指數(shù)和水資源匱乏指數(shù)，分析47個(gè)行政村土地、水及煤炭資源的損耗程度；其次，基于資源交互作用的累積效應(yīng)，構(gòu)建加和評價(jià)模型，評價(jià)煤炭開采與水土資源相互作用產(chǎn)生的綜合影響程度，并進(jìn)行影響特征分析和影響程度劃分；最后，基于影響程度及特征構(gòu)建土地利用功能分區(qū)，并提出相應(yīng)的整治措施。

圖1 煤炭開采與水土資源利用相互影響類型及表征Fig.1 Types and representations of interaction between coal mining and water and soil resources utilization

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

2.1 研究區(qū)概況

長河流域位于山西省澤州縣西北部，包括下村鎮(zhèn)、大東溝鎮(zhèn)、川底鄉(xiāng)3個(gè)鄉(xiāng)鎮(zhèn)的47個(gè)行政村（圖2），面積約108 km2，其中耕地面積為6 472.97 hm2，以旱地為主，占研究區(qū)總面積的59.93%，耕作層為30～150 cm；林、草地面積為3 081.17 hm2，占研究區(qū)總面積的59.95%；區(qū)內(nèi)水域面積為62.60 hm2，其中長河是沁河主要支流之一，面積為36.88 hm2，年徑流量為500萬m3左右。研究區(qū)平均海拔860 m；年平均氣溫11 ℃，年降雨量550～600 mm，無霜期170 d。長河流域位于山西沁水煤田界內(nèi)，域內(nèi)井田面積約有58.78 km2，煤炭資源儲量估算約46 702×104t；有海天、華陽、盈盛、岳圣山等11個(gè)年產(chǎn)大于90×104t/a煤礦企業(yè)，主采煤層為山西組 3號、太原組 9號及 15號。隨著社會經(jīng)濟(jì)發(fā)展對煤炭資源需求提高，流域內(nèi)煤炭年開采量總計(jì)達(dá)到550×104t，對當(dāng)?shù)厮临Y源破壞嚴(yán)重；2005—2015年因煤炭開采造成土地破壞面積約為6 313.54 hm2，占研究區(qū)的58.45%；煤炭開采造成地下水位下降，據(jù)農(nóng)戶訪談得知地下水位從400 m降至約600 m；另外，由于鄉(xiāng)鎮(zhèn)、村莊建設(shè)造成煤炭壓覆儲量達(dá)3 058.5×104t，約占資源總儲量的6.5%，影響煤炭企業(yè)生產(chǎn)。流域內(nèi)土地、水與煤炭資源要素在空間上相互疊加交互作用，主要表現(xiàn)為煤炭開采對水、土地的擾動(dòng)，煤炭廢渣、廢液對水資源的點(diǎn)源污染、建設(shè)用地煤炭壓覆等，因此流域煤炭開采與水土資源利用相互影響特征顯著，具有研究的典型性。

圖2 研究區(qū)地理區(qū)位圖Fig.2 Location of study area

2.2 數(shù)據(jù)來源與處理

圖件主要為2015年7月長河流域礦業(yè)廢棄地及損毀土地現(xiàn)狀分布圖、2015年澤州縣土地利用現(xiàn)狀圖及2014年長河流域DEM圖、11個(gè)煤礦的井上下對照圖及煤層資源儲量估算圖、礦區(qū)土地?fù)p毀預(yù)測圖、澤州縣耕地質(zhì)量評價(jià)成果圖（2012）；自然生態(tài)及社會經(jīng)濟(jì)數(shù)據(jù)主要來源于《澤州縣生態(tài)功能區(qū)報(bào)告》、《煤礦區(qū)土地復(fù)墾方案》、《澤州縣土地利用總體規(guī)劃報(bào)告》、《澤州縣第一次全國水利普查成果匯編（2013）》、《澤州縣統(tǒng)計(jì)年鑒（2012—2015）》以及2015年7—8月的農(nóng)戶調(diào)研和村長訪談。數(shù)據(jù)處理首先將11個(gè)煤礦井工上下對照圖及煤層資源儲量估算圖進(jìn)行矢量化、格式轉(zhuǎn)換和坐標(biāo)校準(zhǔn)，并與2015年土地利用現(xiàn)狀圖進(jìn)行地理配準(zhǔn)，獲得村莊、水庫與煤炭儲量空間分布圖；其次根據(jù)各煤礦土地?fù)p毀預(yù)測圖、廢棄工礦地現(xiàn)狀分布圖與研究區(qū)土地利用現(xiàn)狀圖進(jìn)行疊置分析，獲得研究區(qū)損毀土地空間分布圖；最后對社會經(jīng)濟(jì)數(shù)據(jù)進(jìn)行整理分析，通過 Arcgis將社會經(jīng)濟(jì)數(shù)據(jù)與矢量空間數(shù)據(jù)掛接，最終將上述數(shù)據(jù)聚合轉(zhuǎn)換到以行政村為基本單元中，進(jìn)行綜合影響評價(jià)和整治措施制定；其他類型數(shù)據(jù)直接利用已有數(shù)據(jù)成果。

2.3 煤炭開采與水土資源利用相互影響評價(jià)方法

煤炭開采與水土資源利用相互影響程度評價(jià)旨在描述和評估土地、水及礦產(chǎn)資源利用過程中產(chǎn)生的不利作用大小與空間分布，揭示煤炭開采與水土資源脅迫響應(yīng)形成與演變的機(jī)理，為土地綜合整治提供依據(jù)。本研究將煤炭開采與水土資源利用交互作用產(chǎn)生的綜合影響作為一級評價(jià)準(zhǔn)則層，將流域中主要煤、地和水資源作為主要影響對象，以土地利用系統(tǒng)作為載體，將流域中土地?fù)p毀、煤炭壓覆及水量與水質(zhì)下降等作為影響表征，并以此作為定量評價(jià)指標(biāo)選擇的依據(jù)（圖3）。

圖3 煤炭開采與水土資源利用相互影響評價(jià)框架Fig.3 Evaluation framework of interaction between coal mining and water and soil resources utilization

綜合影響指數(shù)（combined influence index，CII）主要通過土地?fù)p毀累積作用指數(shù)（land-destruction cumulative effect index，LCEI）、水資源匱乏指數(shù)(water poverty index，WPI)和煤炭資源損失指數(shù)(coal loss index，CLI)進(jìn)行度量（圖 3）。由于煤炭開采與水土資源利用的相互影響具有累積性和空間擴(kuò)展性[11]，借鑒物質(zhì)-能量流函數(shù)理論[12]，構(gòu)建評價(jià)公式為

式中LCEIi、WPIi、CLIi表示i行政村的土地?fù)p毀累積作用指數(shù)、水資源匱乏指數(shù)和煤炭資源損失指數(shù)；基于數(shù)據(jù)可獲性和完整性，文中i表示以行政村作為評價(jià)單元。

2.3.1 土地?fù)p毀累積作用指數(shù)評價(jià)方法

土地?fù)p毀累積作用指數(shù)（LCEI）用于表征煤炭開采及水資源利用對土地造成損耗程度大小，可通過土地?fù)p毀程度及地形位指數(shù)[13]（terrain niche index，T）進(jìn)行量化。此值越高，說明此行政區(qū)的土地?fù)p毀程度越高，受其他資源的影響程度越高。

式中，wiy表示i行政村土地?fù)p毀第y個(gè)指標(biāo)的權(quán)重，主要采用結(jié)構(gòu)熵權(quán)法[14]計(jì)算。hil表示i行政村土地?fù)p毀評價(jià)的第l個(gè)指標(biāo)的標(biāo)準(zhǔn)化值，此值可采用極差標(biāo)準(zhǔn)化進(jìn)行計(jì)算，評價(jià)指標(biāo)[15-16]見表1。Ti表示i行政村的地形位指數(shù)，Ai和分別為第i行政村的高程值和研究區(qū)域的平均高程值；Bi和B分別為第i行政村的坡度和該研究區(qū)域的平均坡度值，具體在 ArcGIS10.1平臺上從研究區(qū) DEM中提取坡度和高程值，然后在柵格計(jì)算器中根據(jù)公式計(jì)算得出地形位指數(shù)分布。地形位指數(shù)越低，表明該村高程低、坡度小，土地地形條件平坦，有利于生產(chǎn)、生活空間布局，對整個(gè)區(qū)域開發(fā)越重要，若遭到破壞后，造成的土地?fù)p毀越嚴(yán)重，而地形位指數(shù)越高，表明該行政村高程高、坡度大，地形條件較差，土地開發(fā)利用優(yōu)勢較低，相較之土地?fù)p毀程度較低[13]。

表1 土地?fù)p毀指數(shù)評價(jià)指標(biāo)體系Table 1 Indicators for assessing land destruction degree

主要評價(jià)指標(biāo)說明。1）建筑物損毀指數(shù)是建筑物損毀概率與建筑物損毀程度比值，主要評價(jià)煤炭井工開采對建筑物的損毀程度，進(jìn)而反映土地?fù)p毀情況，該值越高，說明煤炭開采對地表建筑物造成的破壞程度越高。建筑物損毀概率是損毀建筑物棟數(shù)/行政村建筑物總棟數(shù)，數(shù)據(jù)來源于澤州縣地籍調(diào)查現(xiàn)狀圖與各個(gè)村實(shí)地調(diào)查；建筑物損毀程度依據(jù)各行政村建筑物的裂縫條數(shù)和裂縫寬度分為高、中、低程度，分別賦值100、70、40。2）耕地資源場勢指數(shù)主要反映耕地資源的集中程度，該值越高，耕地資源集中程度越高，煤炭開采對耕地影響規(guī)模越大，耕地?fù)p毀程度相應(yīng)也越高。3）水土不匹配指數(shù)主要反映各行政村耕地與水資源匹配狀況，該值越高，說明耕地?fù)p毀程度越低。

2.3.2 煤炭開采損失指數(shù)評價(jià)方法

煤炭開采損失指數(shù)表征土地和水資源利用對煤炭可開采量造成損失程度的大小，主要通過煤炭可開采數(shù)量損失、煤炭開采率和煤炭資源空間分布等方面進(jìn)行評價(jià)[2,15,18-21]（表2）。表中各指標(biāo)權(quán)重通過結(jié)構(gòu)熵權(quán)法計(jì)算，指標(biāo)值通過極差標(biāo)準(zhǔn)化后采用理想景觀向量模型[15]（式4）計(jì)算行政村范圍內(nèi)煤炭損失指數(shù)。指數(shù)值越高，表明煤炭可開采量損失程度越高，水土資源利用對煤炭開采影響程度越高。

式中n是指標(biāo)總數(shù)，vic是i行政村第c個(gè)煤炭損失程度評價(jià)指標(biāo)的數(shù)值，MINC是行政村中第c個(gè)評價(jià)指標(biāo)的最低值。

表2 煤炭開采損失程度評價(jià)指標(biāo)體系Table 2 Indicators for assessing the coal loss degree

主要指標(biāo)說明：1）土地、水利用干擾程度主要反映村莊建設(shè)用地及水庫壓覆煤炭情況，該值越高，表明煤炭壓覆情況越嚴(yán)重，煤炭開采損失程度越高。2）煤炭資源場勢指數(shù)也是反映煤炭在空間上的富集情況，該指數(shù)越高，煤炭資源富集程度越高，煤炭壓覆概率越高，造成煤炭可采規(guī)模損失越嚴(yán)重。

2.3.3 水資源匱乏指數(shù)評價(jià)方法

水資源匱乏指數(shù)（water poverty index, WPI）表征煤炭開采及土地資源利用對水資源造成損失程度的大小。主要從水資源數(shù)量損失、水質(zhì)污染和水資源空間分布等方面進(jìn)行評價(jià)，參考相關(guān)文獻(xiàn)確定評價(jià)指標(biāo)[18-20]（表3）。表中各指標(biāo)值權(quán)重和標(biāo)準(zhǔn)化值計(jì)算同上，最后采用理想景觀向量模型式（5）計(jì)算行政村范圍內(nèi)水資源匱乏指數(shù)。該值越高，表明行政村內(nèi)由于其他資源的利用而造成水資源數(shù)量下降、水質(zhì)污染而不能滿足區(qū)域生產(chǎn)、生活用水的缺口也越大。

式中izv是i行政村第z個(gè)評價(jià)指標(biāo)的數(shù)值，MINz是i行政村中第z個(gè)評價(jià)指標(biāo)的最低值。

2.4 土地整治控制策略

受漸變性特征的約束，根據(jù)綜合影響程度分級是土地整治對策與措施制定的基礎(chǔ)。土地整治策略主要是在流域煤炭開采與水土資源交互作用產(chǎn)生的綜合影響評價(jià)結(jié)果的基礎(chǔ)上，針對綜合影響程度及其類型產(chǎn)生的主導(dǎo)限制因素確定土地利用功能分區(qū)（圖 4），并根據(jù)對應(yīng)分區(qū)確定土地整治級別（即優(yōu)先性）及整治策略。

表3 水資源匱乏程度評價(jià)指標(biāo)體系Table 3 Indicators for assessing water poverty degree

圖4 土地整治策略技術(shù)流程Fig.4 Flow chart for land consolidation strategy

根據(jù)綜合影響評價(jià)結(jié)果確定整治的優(yōu)先性和差異性，針對土地利用功能分區(qū)主導(dǎo)影響因素提出修復(fù)、保護(hù)、復(fù)墾和提升 4個(gè)整治措施類型：修復(fù)型主要是針對綜合影響程度高且煤炭、土地和水資源損失較高的區(qū)域，不利于農(nóng)戶生產(chǎn)、生活，對此類區(qū)域在優(yōu)先煤炭開采前提下，實(shí)施適度村莊搬遷或零星居民點(diǎn)整理，釋放壓煤儲量，提高煤炭可開采數(shù)量，并通過村莊搬遷抽稀人口密度以及生態(tài)修復(fù)工程，以緩解資源利用矛盾；保護(hù)型主要是針對綜合影響程度高但水、土資源損毀中等的區(qū)域，適度控制礦產(chǎn)開采規(guī)模，重點(diǎn)進(jìn)行生態(tài)資源保護(hù)；復(fù)墾型主要針對綜合影響程度為中度，壓煤儲量中等，但土地資源損毀重度的地區(qū)，通過合理控制礦產(chǎn)開采，對損毀土地進(jìn)行復(fù)墾使其恢復(fù)到原先可以利用的狀態(tài)；提升型主要是針對綜合影響程度為輕度或（基本）無，且水、土資源空間集聚的區(qū)域，以現(xiàn)代產(chǎn)業(yè)引領(lǐng)型的標(biāo)準(zhǔn)化土地整治為主，提高耕地質(zhì)量等級，為承接現(xiàn)代產(chǎn)業(yè)做準(zhǔn)備。

3 結(jié)果與分析

3.1 土地?fù)p毀累積作用評價(jià)結(jié)果

經(jīng)過實(shí)地調(diào)查，流域采煤塌陷地、居民建筑損毀點(diǎn)、煤矸場以及其他工礦廢棄地，主要分布在區(qū)域西部（圖5a）。通過土地?fù)p毀指數(shù)評價(jià)，并采用GIS中自然裂點(diǎn)法將其分為3類區(qū)域（圖5b）。

圖5 土地?fù)p毀及其累積評價(jià)指數(shù)空間分布圖Fig.5 Distribution for land destruction and its cumulative effect index

1）土地高度損毀區(qū)（>0.137～1.00），涉及15個(gè)行政村，面積為4 775.90 hm2，主要包括：耕地與煤炭采區(qū)空間疊置面積有2 580.85 hm2，占高度損毀區(qū)總面積的54.4%；林、草地有1 247.29 hm2，占高度損毀區(qū)總面積的26.11%。該區(qū)位于長河西部丘陵區(qū)，是研究區(qū)煤礦開采工作面集中區(qū)，因此煤炭開采導(dǎo)致土地大面積裂縫和塌陷以及居民建筑損毀，如常坡村田塊最寬的裂縫達(dá)30 cm，長數(shù)十米，房屋沉陷和裂縫情況嚴(yán)重，以及王虎山村多處房屋沉陷可達(dá)10 cm左右，裂隙高達(dá)約5 m。

2）土地中度損毀區(qū)（>0.099～0.137）涉及 22個(gè)行政村，面積為5 008.005 3 hm2，主要位于：流域中南部，耕地與采區(qū)疊置面積1 975.87 hm2，占中度損毀區(qū)總面積的39.45%；林、草地面積為1 096.99 hm2，占21.90%。該區(qū)長河以西是受海天、華陽等煤礦開采干擾造成塌陷、廢棄壓占土地集中分布，且當(dāng)?shù)胤课萦休p微裂縫，如溝西、中街、原莊等；長河以東是由于鄰近煤礦開采的干擾，耕地出現(xiàn)小規(guī)模裂縫、塌陷等情況，大多可農(nóng)戶自行修復(fù)，但耕地漏肥走水現(xiàn)象普遍，耕地質(zhì)量較低多為4～7等，且植物覆蓋率低，地勢較高處有三級土壤侵蝕，增加土地?fù)p毀程度，主要包括賀坡、庾能、北村等村莊。

3）土地輕度損毀區(qū)（0～0.099）面積為 1 537.795 9 hm2，主要位于：流域中部和南部，包括李山村、劉村、東溝村等，耕地與采區(qū)空間疊置面積4.30 hm2，占0.27%，林、草地面積為417.39 hm2，占27.14%。該區(qū)社會經(jīng)濟(jì)發(fā)展較快，是城鎮(zhèn)發(fā)展和人口聚居中心，且煤炭開采活動(dòng)干擾較低，根據(jù)實(shí)地調(diào)查無土地塌陷損毀及房屋裂縫等情況，地勢較平坦，二級土壤侵蝕，生態(tài)用地比例約為21%，土地?fù)p毀程度較低。

3.2 水資源匱乏評價(jià)結(jié)果

水匱乏指數(shù)（WPI）主要反映由于水資源狀況與行業(yè)需水之間的差異[22]。在本研究中反映由于其他資源利用的干擾下，水資源對生產(chǎn)、生活用水需求的供給能力。

流域地表水資源主要分布在研究區(qū)中西部地區(qū)，主要包括長河、沙溝水庫、寺河水庫及其他 4座水庫，地下水資源的60%左右是埋藏于地面400 m以下的深水層。研究區(qū)缺水率平均為 30%左右，水資源供需比均<1。另外經(jīng)過實(shí)地調(diào)研，煤炭開采對地下水資源擾動(dòng)（平均噸煤開采破壞水量為7 t左右）及農(nóng)田化肥點(diǎn)源污染是研究區(qū)水資源損毀的主要因素。研究區(qū)水匱乏指數(shù)評價(jià)結(jié)果（圖6）顯示，并采用自然裂點(diǎn)法將其分為3類區(qū)域。

水匱乏指數(shù)高值區(qū)（>0.272～1.00）面積為3 077.02 hm2，涉及15個(gè)行政村，此區(qū)無煤炭開采活動(dòng)干擾，無地表水，地下水屬于深層水，主要依靠提水井滿足生產(chǎn)、生活需求。水匱乏指數(shù)中值區(qū)（>0.241～0.272）面積為5 547.64 hm2，涉及22個(gè)行政村，地表水面積為28.67 hm2。該區(qū)緊鄰長河，為沁河的一大支流，常年性河流，年徑流量 500～700萬 m3，除雨季洪水之外一般水量不大，主要滿足于長河附近的農(nóng)田灌溉，輻射范圍較小。水匱乏指數(shù)低值區(qū)（0～0.241）面積為2 697.04 hm2，涉及10個(gè)行政村，地表水面積28.55 hm2。該區(qū)為煤炭集中開采區(qū)域，礦井抽水后地下水位下降導(dǎo)致提水井干枯，但礦井下涌水（平均單井涌水量1 500 m3/d）經(jīng)抽到地面后，經(jīng)過處理，水質(zhì)優(yōu)良，仍能滿足礦區(qū)生產(chǎn)與居民生活用水需要。

圖6 水資源匱乏指數(shù)空間分布圖Fig.6 Distribution of water poverty index (WPI)

3.3 煤炭開采損失評價(jià)結(jié)果

煤炭開采損失主要是由于村莊、水域及水利設(shè)施等壓煤以及煤層充水，限制地下煤炭開采活動(dòng)而導(dǎo)致煤炭開采率下降引起煤炭開采損失。研究區(qū)煤炭開采損失指數(shù)評價(jià)結(jié)果（圖7）顯示，并采用自然裂點(diǎn)法將其分為3類區(qū)域。

煤炭開采損失指數(shù)高值區(qū)（>0.223～1）面積2 286.70 hm2，分布在流域西部，涉及9個(gè)行政村，是煤炭主采區(qū)。根據(jù)《煤炭資源開發(fā)方案》得出該區(qū)煤炭壓覆面積為274.91 hm2，占12.02%，壓煤儲量約2 017.5×104t。主要是因?yàn)榇藚^(qū)含9號（平均厚度1.55 m）和15號（平均厚度 2.16 m）主要開采煤層，由于煤炭壓覆造成煤炭開采損失。煤炭開采損失指數(shù)中值區(qū)（>0.081～0.223）面積為3 340.56 hm2，涉及13個(gè)行政村，煤炭壓覆面積為126.35 hm2，占3.78%。由于該區(qū)雖含有3號、9號及15號煤層但規(guī)模較小，壓煤儲量約932×104t。煤炭開采損失指數(shù)低值區(qū)（0～0.081）面積為5 694.44 hm2，涉及25個(gè)行政村。該區(qū)位于沁水煤田邊緣，僅有有零星壓煤區(qū)域且只為14.90 hm2，占0.26%，壓煤儲量約109×104t。

3.4 綜合影響程度評價(jià)結(jié)果

經(jīng)過式（1）計(jì)算，流域綜合影響指數(shù)值主要在0～3之間，并采用GIS中自然斷裂法將綜合影響程度分為高、中、低和（基本）無等4級（圖8）。

1）綜合影響指數(shù)高值區(qū)（0.599～3），面積為2 286.695 7 hm2（占流域總面積的20.21%），主要分布在流域西北部和西南部，包括①高煤炭開采損失區(qū)域，其內(nèi)煤炭資源豐富，分布9號和15號煤層，同時(shí)壓煤村莊數(shù)量達(dá)9個(gè)左右，導(dǎo)致煤炭壓覆規(guī)模較大；②高土地?fù)p毀區(qū)域占高值區(qū)的81.86%，分布在流域的西北部和西南部，主要是受煤炭開采活動(dòng)集中的影響，出現(xiàn)土地塌陷、裂縫及壓占，以及由此引發(fā)的地質(zhì)災(zāi)害及水土流失等生態(tài)環(huán)境問題，中/低土地?fù)p毀面積較??；③中度水資源匱乏區(qū)域，主要是由于煤炭開采對水資源數(shù)量、質(zhì)量的損耗，但又由于礦井涌水的高效利用，一定程度彌補(bǔ)水資源損耗。

圖7 壓煤村莊（水庫）及其煤炭損失評價(jià)指數(shù)空間分布圖Fig.7 Distribution for coal village (reservoir)and coal loss index (CLI)

2）綜合影響指數(shù)中值區(qū)（>0.439～0.599），面積為3 780.80 hm2（33.39%），主要分布在長河西部， 主要包括①高土地?fù)p毀區(qū)域，主要是西部丘陵區(qū)，區(qū)內(nèi)土地利用類型主要為耕地、村莊和零星的林、草地，土地?fù)p毀主要表現(xiàn)為煤炭開采導(dǎo)致土地塌陷及裂縫、村莊建筑物沉陷及產(chǎn)生的裂縫；②該區(qū)域有沙溝水庫、寺河水庫及其他 4座小型水庫，但地表水蓄量僅能滿足周邊村莊生產(chǎn)、生活用水，供給范圍有限。

圖8 綜合影響空間分布圖Fig.8 Distribution for combined influence index (CII)

3）綜合影響指數(shù)低值區(qū)（>0.387～0.439），面積為2 745.61 hm2（24.24%），主要分布在流域東北部，主要包括①長河?xùn)|側(cè)的耕地，臨近煤炭采掘工作面，受到煤炭開采的擾動(dòng)，但可自行修復(fù)；②東部丘陵區(qū)林地，具有蓄水涵養(yǎng)、水土保持的功能；③該區(qū)無地表水源，且地下水埋深達(dá)400 m，用水成本高。

4）（基本）無綜合影響區(qū)（0～0.387），面積為2 508.60 hm2（22.16%），主要分布在流域東南部，該區(qū)臨近長河且無礦業(yè)生產(chǎn)的擾動(dòng)，長河沿岸有部分水澆地，但以傳統(tǒng)旱作為主。

3.5 基于綜合影響程度的土地整治策略

依據(jù)土地整治技術(shù)流程（圖4），結(jié)合流域綜合影響程度空間分布及土地利用主要問題，將流域土地利用空間分為七大土地利用功能區(qū)（圖9），構(gòu)建煤炭與水土資源綜合利用的整治策略矩陣（表4）。

1）重點(diǎn)生態(tài)修復(fù)區(qū)（Z1）屬于I1級綜合影響控制區(qū)，主要包括高損毀的土地資源及高損失的煤炭資源。這些地區(qū)蘊(yùn)含豐富的煤炭資源，但因村莊及水庫壓煤而造成開采效率下降，同時(shí)這些地區(qū)地勢較高，土地、水資源受到煤炭開采活動(dòng)干擾而易引發(fā)土地塌陷及水土流失。因此，此區(qū)應(yīng)在保證煤炭開采活動(dòng)的前提下，重點(diǎn)采取坡改梯工程、水土保持工程和土壤改良等改良措施，如以蓄水、保土為核心，通過增加表土厚度，平整土地、修筑田坎、修筑蓄水灌排系統(tǒng)，種植護(hù)坡林草等，恢復(fù)地表植被，提高耕地質(zhì)量，維護(hù)生態(tài)環(huán)境。

2）重點(diǎn)生態(tài)保護(hù)區(qū)（Z2）屬于I2級綜合影響控制區(qū)，主要包括受損嚴(yán)重土地。區(qū)內(nèi)土地因地下開采而塌陷、裂縫，水肥流失，地基承載力下降，因煤矸石壓占、廢棄地等威脅土地高效利用，但由于臨近水庫以及礦井涌水彌補(bǔ)，水資源匱乏程度較低，因此該地區(qū)在適度限制煤炭開采活動(dòng)的前提下，通過農(nóng)田平整、植被恢復(fù)及村莊、廢棄地整治等多項(xiàng)措施保護(hù)該區(qū)田、水、林、村綜合生態(tài)系統(tǒng)。

3）損毀土地復(fù)墾區(qū)（Z3）屬于II1等級綜合影響控制區(qū)，主要是土地集中損毀，重點(diǎn)通過地貌重塑、土壤重構(gòu)、林網(wǎng)構(gòu)建及基礎(chǔ)設(shè)施建設(shè)等工程措施進(jìn)行損毀土地的恢復(fù)和重建，可復(fù)墾為農(nóng)田、林草地或休閑農(nóng)業(yè)區(qū)等。

4）耕地資源保護(hù)區(qū)（Z4）屬于II2等級綜合影響控制區(qū)，主要是土地?fù)p毀中度，不影響正常的農(nóng)業(yè)生產(chǎn)和工程建設(shè)，且有一定水源補(bǔ)充，因此在禁止對耕地?fù)p害大的生產(chǎn)建設(shè)活動(dòng)的前提下，本區(qū)重點(diǎn)通過土壤重構(gòu)、土地平整、農(nóng)田水利等措施，改善耕作條件，重點(diǎn)保護(hù)耕地資源。

5）水資源匱乏修復(fù)區(qū)（Z5）屬于III1 等級綜合影響控制區(qū)，主要水資源破壞嚴(yán)重不能滿足生產(chǎn)、生活需求。因此該區(qū)的土地整治以節(jié)水設(shè)施工程建設(shè)為主導(dǎo)，采取渠道防滲、地下管灌、地上膜灌、噴微灌等措施，配套土地平整、田間道路、農(nóng)田防護(hù)與生態(tài)環(huán)境保持等工程，實(shí)現(xiàn)節(jié)水農(nóng)田建設(shè)。

6）水資源與耕地資源保護(hù)區(qū)（Z6）屬于III2 等級綜合影響控制區(qū)，主要位于區(qū)域西南部和東北部，煤炭開采對耕地資源干擾度較小且緊鄰長河，因此該區(qū)主要通過土地整治完善農(nóng)田灌溉排水設(shè)施工程，開展田間道路工程，并改善田塊形狀，開展內(nèi)部零星地類的整治，促進(jìn)農(nóng)田集中連片，以及有效的農(nóng)田節(jié)水工程建設(shè)，盡量滿足水土資源平衡為農(nóng)業(yè)規(guī)模經(jīng)營創(chuàng)造條件。

7）耕地質(zhì)量提升區(qū)（Z7）屬于資源利用協(xié)調(diào)區(qū)，主要位于流域東南部，為低度或基本無損的土地，但水資源較匱乏。此區(qū)煤炭開采干擾較小，耕地資源豐富，因此一方面可通過土地平整工程改善耕作條件、增加農(nóng)田水利設(shè)施、提高耕地質(zhì)量，另一方面可通過土地整治與現(xiàn)代農(nóng)業(yè)發(fā)展相結(jié)合按照產(chǎn)業(yè)化發(fā)展和農(nóng)業(yè)生產(chǎn)總體布局的要求，科學(xué)設(shè)計(jì)農(nóng)田整治方案，為現(xiàn)代農(nóng)業(yè)建設(shè)提供基礎(chǔ)。

圖9 基于綜合影響指數(shù)的土地利用功能分區(qū)Fig.9 Land use performance zoning based on land-mining utilization conflicts

目前，長河流域進(jìn)入煤炭資源產(chǎn)業(yè)轉(zhuǎn)型期，根據(jù)以上結(jié)果，土地整治應(yīng)結(jié)合產(chǎn)業(yè)發(fā)展需求與資源利用情況進(jìn)行統(tǒng)一規(guī)劃，例如，綜合影響較低的區(qū)域以現(xiàn)代產(chǎn)業(yè)用地標(biāo)準(zhǔn)進(jìn)行整治，成為現(xiàn)代產(chǎn)業(yè)承接地；而綜合影響較高的區(qū)域亟需納入生態(tài)恢復(fù)計(jì)劃，滿足當(dāng)?shù)鼐用瘛叭毙枰欢鴮τ诰C合影響中等的區(qū)域，對產(chǎn)生影響的主導(dǎo)因素進(jìn)行有序、有效整治，成為流域產(chǎn)業(yè)轉(zhuǎn)型用地需求的后備力量。

3.6 綜合影響評價(jià)與整治策略的有效性分析

依據(jù)《澤州縣土地利用總體規(guī)劃（2006—2020）》，將劃分的流域土地利用功能分區(qū)圖層與土地用途分區(qū)圖層疊加分析煤炭開采與水土資源利用綜合影響空間分區(qū)的合理性。表 5顯示，以耕地和水資源保護(hù)為主導(dǎo)的區(qū)域（Z5、Z6、Z7）中基本農(nóng)田保護(hù)區(qū)規(guī)劃比例超過40.55%～68.12%，城鎮(zhèn)村建設(shè)用地規(guī)劃比例在 4.95%～12.78%之間；以損毀土地復(fù)墾為主導(dǎo)的區(qū)域（Z3、Z4）中規(guī)劃基本農(nóng)田比例在 25%～37.77%之間，而獨(dú)立工礦用地區(qū)為 7.39%～15.07%之間，高于其他區(qū)同類用地比例；以生態(tài)保護(hù)為主導(dǎo)的區(qū)域（Z1、Z2）中林業(yè)用地區(qū)規(guī)劃比例為26.94%～30.04%，高于其他區(qū)域林地規(guī)劃比例。由此看出基于綜合影響評價(jià)結(jié)果的流域礦區(qū)土地利用功能分區(qū)結(jié)果與土地利用用途分區(qū)規(guī)劃基本吻合，可為長河流域永久性基本農(nóng)田保護(hù)紅線、城鎮(zhèn)建設(shè)邊界及生態(tài)保護(hù)紅線提供一定參考依據(jù)。另據(jù)《澤州縣土地整治規(guī)劃（2011—2015）》中要求長河流域提高基本農(nóng)田、農(nóng)業(yè)現(xiàn)代化及土地復(fù)墾等建設(shè)項(xiàng)目的比例，與面向綜合影響評價(jià)的整治策略基本吻合。

表5 土地利用功能分區(qū)與2020年土地利用規(guī)劃對照表Table 5 Comparison between land use zoning and land use planning result of 2020 %

4 討 論

1）基于流域煤炭開采與水土資源相互影響的類型與表征，建立綜合影響評價(jià)模型，能夠有效地對區(qū)域煤炭開采與水土資源相互影響水平等級評價(jià)，而且能夠?qū)雾?xiàng)資源要素的優(yōu)劣進(jìn)行判斷，從而揭示更多的信息，為土地綜合整治提供科學(xué)依據(jù)。

關(guān)于煤炭開采與水土資源利用關(guān)系定量評價(jià)，劉斌等[23]根據(jù)區(qū)域內(nèi)景觀格局和破碎度，采用生態(tài)指數(shù)、生態(tài)脆弱度指數(shù)和潛在生態(tài)損失度指數(shù)測算景觀生態(tài)風(fēng)險(xiǎn)值。常青等[24-25]從生態(tài)敏感度、土地?fù)p毀累積作用系數(shù)和暴露系數(shù)對礦區(qū)土地?fù)p毀生態(tài)風(fēng)險(xiǎn)進(jìn)行度量。馮宇等[26]則從空間外部壓力、空間暴露性及空間穩(wěn)定性對流域煤炭開采引起的土地利用沖突進(jìn)行研究。部分學(xué)者也采用生態(tài)足跡[27]、生態(tài)梯度[28]、耦合模型[29]等評價(jià)方法對煤炭開采的資源效應(yīng)進(jìn)行評價(jià)。但是多數(shù)學(xué)者側(cè)重煤炭開采對水土資源利用的單向影響，忽視兩者之間的相互影響，另一方面從評價(jià)指標(biāo)體系來看多為宏觀層面指標(biāo)，無法從微觀分析單項(xiàng)影響因素的優(yōu)劣。鑒于此，本研究從立體空間的角度，對地下煤炭資源開采和地上水土資源利用之間相互影響機(jī)理及表征進(jìn)行分析，并以土地、水和煤炭為相互作用對象，從數(shù)量、質(zhì)量和空間等 3方面屬性建立評價(jià)指標(biāo)體系，構(gòu)建綜合影響評價(jià)模型；同時(shí)，為了顯化單項(xiàng)資源的損耗狀態(tài)，構(gòu)建土地?fù)p毀指數(shù)、煤炭損失指數(shù)和水資源匱乏指數(shù)，綜合評價(jià)流域煤炭開采與水土資源利用相互影響程度，豐富了礦區(qū)資源利用效應(yīng)評價(jià)定量研究方法，也為土地整治策略制定提供有效依據(jù)。

2）該研究是在流域土地?fù)p毀評價(jià)基礎(chǔ)上擴(kuò)展成煤炭開采與水土資源相互影響綜合評價(jià)，以數(shù)量、質(zhì)量和空間等屬性作為評價(jià)落腳點(diǎn)，以此制定和探索整治控制策略。然而在資源利用過程中，資源利用相互影響涉及不同時(shí)空尺度，不同時(shí)空尺度所涉及的問題及影響表現(xiàn)形式也有所不同，相應(yīng)的整治規(guī)劃與措施也應(yīng)有所調(diào)整，要實(shí)施有效的資源利用整治措施，必須依據(jù)村級層面的微觀尺度研究結(jié)果[30]。另外，由于煤炭開采與水土資源利用方式形成和變化與研究區(qū)生產(chǎn)、生活與生態(tài)等因素緊密聯(lián)系，即除了土地要素是資源利用相互影響主要受體及整治對象，其他如農(nóng)戶生計(jì)、糧食安全、生態(tài)文明等社會經(jīng)濟(jì)因素也影響資源利用。因此，研究視角放大到“煤炭-水土”資源利用和“社會-經(jīng)濟(jì)”交互復(fù)合系統(tǒng)中，資源利用表現(xiàn)形式和影響范圍將會更加復(fù)雜和廣泛。今后研究中可進(jìn)一步從村級尺度深入探討和分析“煤炭-水土”資源綜合利用在“社會-經(jīng)濟(jì)”復(fù)合系統(tǒng)中的發(fā)展和表現(xiàn)形式及整治策略。

5 結(jié) 論

1）通過實(shí)地調(diào)研及GIS技術(shù)，對煤炭開采與水土資源利用相互影響進(jìn)行綜合評價(jià)，并將其劃分為高、中、低和基本無影響區(qū)，能夠揭示流域單一資源損毀及其資源利用交互作用的區(qū)域差異；綜合土地、水及煤炭資源損毀定量評價(jià)進(jìn)行綜合影響程度分級分區(qū)，能夠更好地揭示不同影響區(qū)內(nèi)主導(dǎo)限制因素，為明確土地整治對象提供依據(jù)。

2）資源利用綜合影響高、中、低值區(qū)及基本（無）影響區(qū)，所占研究區(qū)總面積的比例分別是 20.21%、33.39%、24.24%和22.16%。綜合影響高值區(qū)主要分布在流域西北部和西南部，煤炭利用與水土資源相互影響的特征是煤炭壓覆、高度土地?fù)p毀及中度水資源匱乏；中值區(qū)主要分布在流域西部，高度土地?fù)p毀是主要特征；低值區(qū)主要分布在流域東北部，水資源匱乏是主要特征；無影響區(qū)主要分布在流域東部，無煤炭生產(chǎn)活動(dòng)擾動(dòng)，且緊鄰長河便于農(nóng)田灌溉。

3）通過對以上綜合影響程度異質(zhì)區(qū)進(jìn)行修復(fù)、保護(hù)、復(fù)墾和提升等不同等級的整治策略優(yōu)先排序，可將流域土地功能區(qū)分為重點(diǎn)生態(tài)修復(fù)區(qū)（Z1）、重點(diǎn)生態(tài)保護(hù)區(qū)(Z2)、 損毀土地復(fù)墾區(qū)（Z3）、耕地資源保護(hù)區(qū)（Z4）、水資源匱乏修復(fù)區(qū)（Z5）、水資源與耕地資源保護(hù)區(qū)（Z6）和耕地質(zhì)量提升區(qū)（Z7）等七大功能，其結(jié)果與澤州縣第三輪土地利用總體規(guī)劃中土地用途分區(qū)相比較，對應(yīng)的主要土地用途分區(qū)分別是林業(yè)用地（占Z1的30.04%）、林業(yè)用地（占 Z2的 26.94%）、一般農(nóng)地區(qū)（占 Z3的28.21%）及基本農(nóng)田保護(hù)區(qū)（占Z4、Z5、Z6及Z7區(qū)的37.77%以上），研究結(jié)果與土地規(guī)劃基本吻合，說明此法具有一定的有效性，可為流域土地整治規(guī)劃提供依據(jù)。

該研究是從立體空間角度對地下煤炭資源和地上水土資源之間相互影響評價(jià)并制定相應(yīng)的整治策略的嘗試與探索，該方法有利于土地利用規(guī)劃和土地整治規(guī)劃有機(jī)結(jié)合，提高土地整治的有效性。今后在研究中，可在此基礎(chǔ)上把社會經(jīng)濟(jì)因素考慮進(jìn)來，以進(jìn)一步完善流域資源利用綜合評價(jià)與整治研究。

[1] 謝高地. 自然資源總論[M]. 高等教育出版社，2009，2：220－223.

[2] 蔣冬梅．煤炭城市地-礦沖突分析及其調(diào)控研究[D]．徐州：中國礦業(yè)大學(xué)，2016．Jiang Dongmei. Research on the Land-Mining Conflict and its Control in Coal City[D]. Xuzhou: China University of Mining & Technology, 2016. (in Chinese with English abstract)

[3] Craynon J R, Sarver E A, Ripepi N S, et al. A GIS-based methodology for identifying sustainability conflict areas in mine design-a case study from a surface coal mine in the USA[J]. International Journal of Mining, Reclamation and Environment, 2015: 1－12.

[4] Wawryk A. Conservation and access to land for mining in protected areas: the conflict over mining in south Australia’s Arkaroola wilderness sanctuary[J]. Journal of Environmental Law, 2014(26): 291－317.

[5] Hartje V. US–German experience with funding post mining restoration[C]//2015 Sino-German Symposium Notice“Innovative Planning Methods and Instruments for Post Mining Regions-Key Factors for Sustainable Development”,2015, 10.

[6] 婁華君，莊健鴻．煤礦開采區(qū)水、土地與煤炭資源同步利用模式研究[J]．資源科學(xué)，2007，29(5)：90－96．Lou Huajun, Zhuang Jianhong. The sync-utilization mode of water, land and coal resourcesin coal mining subsidence area[J]. Resources Science, 2007, 29(5): 90－96. (in Chinese with English abstract)

[7] 李鋼．礦地統(tǒng)籌發(fā)展及關(guān)鍵技術(shù)探討[J]．中國土地，2015(6)：55－58．Li Gang. Land-mines development and key technologies discussed[J]. China Land, 2015(6): 55－58. (in Chinese with English abstract)

[8] 陳小芳，鄧?？担剿魉拗菔谐堑V鄉(xiāng)一體化發(fā)展的新模式[J]．赤峰學(xué)院學(xué)報(bào)(自然科學(xué)版)，2014，30(17)：39－41．Chen Xiaofang, Deng Fukang. Explore the city of Suzhou city township development integration of the new mode[J].Journal of Chifeng University (Natural Science Edition),2014, 30(17): 39－41. (in Chinese with English abstract)

[9] 余星滌，趙曉方．礦糧復(fù)合區(qū)資源多用途選擇的探討[J]．中國土地科學(xué)，2006，20(8)：35－38．Yu Xingdi, Zhao Xiaofang. Discussion on multiple choices of resources use in mine-grain mixed zone[J]. China Land Science, 2006, 20(8): 35－38. (in Chinese with English abstract)

[10] 王云．山西省采煤土地塌陷損耗分析與測算[J]．能源環(huán)境保護(hù)，2007，21(6)：10－13．Wang Yun. Coal mining subsidence land loss analysis and measurement in Shanxi[J]. Energy Environmental Protection,2007, 21(6): 10－13. (in Chinese with English abstract)

[11] 王行風(fēng)．煤礦區(qū)生態(tài)環(huán)境累積效應(yīng)研究[D]．徐州：中國礦業(yè)大學(xué)，2010．Wang Xingfeng. Study on eco-environmental cumulative effects in coal mining area[J]. Xuzhou: China University of Mining and Technology, 2010. (in Chinese with English abstract)

[12] 馬修軍，鄔倫，謝昆青．空間動(dòng)態(tài)模型建模方法[J]，北京大學(xué)學(xué)報(bào)(自然科學(xué)版)，2004，40(2)：279－286．Ma Xiujun, Wu Lun, Xie Kunqing. The method for modeling on spatial dynamic process[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2004, 40(2): 279－286. (in Chinese with English abstract)

[13] 饒為民，章家恩，肖紅生，等．基于地形位的城郊景觀分布特征及變化[J]．華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2007，28(3)：67－70.Rao Weimin, Zhang Jiaen, Xiao Hongsheng, et al.Distribution characteristics and dynamic change of suburb landscape based on the terrain niche[J]. Journal of South China Agricultural University, 2007, 28(3): 67－70. (in Chinese with English abstract)

[14] 焦鵬飛，張鳳榮，李燦，等．基于引力模型的縣域中心村空間布局分析：以山西省長治縣為例[J]．資源科學(xué)，2014，36(1)：45－54．Jiao Pengfei, Zhang Fengrong, Li Can, et al. The spatial organization of central villages in Changzhi county based on the gravity model[J]. Resources Science, 2014, 36(1): 45－54.(in Chinese with English abstract)

[15] 常青，劉丹，劉曉文．礦業(yè)城市土地?fù)p毀生態(tài)風(fēng)險(xiǎn)評價(jià)與空間防范策略[J]．農(nóng)業(yè)工程學(xué)報(bào)，2013，29(20)：245－254．Chang Qing, Liu Dan, Liu Xiaowen. Ecological risk assessment and spatial prevention tactic of land destruction in mining city[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013,29(20): 245－254. (in Chinese with English abstract)

[16] 程琳琳，婁尚，劉巒峰，等．礦業(yè)廢棄地再利用空間結(jié)構(gòu)優(yōu)化的技術(shù)體系與方法[J]．農(nóng)業(yè)工程學(xué)報(bào)，2013，29(7)：207－218．Cheng Linlin, Lou Shang, Liu Luanfeng, et al. Technology system and method of spatial structure optimization for mining wasteland reuse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(7): 207－218. (in Chinese with English abstract)

[17] 喬文俊．地質(zhì)塊段法在煤炭資源儲量估算中的應(yīng)用[J]．內(nèi)蒙古煤炭經(jīng)濟(jì)，2013，1：45－47．Qiao Wenjun. Application of geological block method in coal reserve estimate[J]. Inner Mongolia Coal Economy, 2013, 1:45－47. (in Chinese with English abstract)

[18] 周洪，谷樹忠，姚予龍，等．基于資源場勢的我國煤炭空間配置特征研究[J]．自然資源學(xué)報(bào)，2013，28(7): 1094－1105．Zhou Hong, Gu Shuzhong, Yao Yulong, et al. Characteristics of coal resources spatial allocation in China: Based on resource field and potency[J]. Journal of Natural Resources,2013，28(7): 1094－1105. (in Chinese with English abstract)

[19] 楊軍，吳蔚，叢建輝，等．煤炭開采與水資源利用保護(hù)的耦合關(guān)系研究：基于中國 10個(gè)煤炭資源豐裕省份的分析[J]．生態(tài)經(jīng)濟(jì)，2016，32(2)：55－59．Yang Jun, Wu Wei, Cong Jianhui, et al. Coupled research of coal mining and protection and utilization of water resource:Based on the analysis of China 10 provinces in the abundance of coal resources[J]. Ecological Economy, 2016, 32(2): 55－59. (in Chinese with English abstract)

[20] 劉一瑋．高潛水位煤礦區(qū)水土資源協(xié)調(diào)利用研究[D]．徐州：中國礦業(yè)大學(xué)，2014．Liu Yiwei. Study on Coordinated Use of Water and Soil Resources of High Phreatic Water Level in Coal Mine Area[D]. Xuzhou: China University of Mining & Technology,2014. (in Chinese with English abstract)

[21] 邱立新，杜銘華，周田君．西北煤炭開發(fā)對水環(huán)境影響的評價(jià)模型與方法[J]．干旱區(qū)資源與環(huán)境，2008，22(2)：20－25．Qiu Lixin, Du Minhua, Zhou Tianjun. A model and a method for assessment of environmental impact of northwest Chinese coal development on water resources[J]. Journal of Arid Land Resources and Environment, 2008, 22(2): 20－25. (in Chinese with English abstract)

[22] 姜寧，付強(qiáng)，孫穎娜．水匱乏指數(shù)及其在東北地區(qū)水安全評價(jià)中的應(yīng)用[J]．東北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013，44(2)：88－91．Jiang Ning, Fu Qiang, Sun Yingna. Application of water poverty index in water security evaluation in northeast region of China[J]. Journal of Northeast Agricultural University,2013, 44(2): 88－91. (in Chinese with English abstract)

[23] 劉斌，冀巍，丁長春．區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)方法研究[J]．科技創(chuàng)新與應(yīng)用，2013(11)：118－119．

[24] 賈媛，曹玲嫻. 煤炭礦區(qū)生態(tài)風(fēng)險(xiǎn)評價(jià)方法研究[J]．環(huán)境科學(xué)與管理，2011，36(4)：177－182.Jia Yuan, Cao Lingxian. Review on ecological assessment method for coal-mining area[J]. Environmental Science and Management, 2011, 36(4): 177－182. (in Chinese with English abstract)

[25] 常青，邱瑤，謝苗苗，等．基于土地破壞的礦區(qū)生態(tài)風(fēng)險(xiǎn)評價(jià)：理論與方法[J]．生態(tài)學(xué)報(bào)，2012，32(16)：5164－5174.Chang Qing, Qiu Yao, Xie Miaomiao, et al. Theory and method of ecological risk assessment for mining areas based on the land destruction[J]. Acta Ecologica Sinica, 2012,32(16): 5164－5174. (in Chinese with English abstract)

[26] 馮宇，畢如田，王瑾，等．流域礦業(yè)開采引發(fā)的土地利用空間沖突及優(yōu)化配置[J]．中國土地科學(xué)，2016，30(11)：32－40.Feng Yu, Bi Rutian, Wang Jin, et al. Spatial conflict of land use caused by mining exploitation and optimal allocation scheme of land resources in river basin, 2016, 30(11): 32－40.(in Chinese with English abstract)

[27] 劉瑛，楊丹青，楊曉暉．煤礦開采對礦業(yè)城鎮(zhèn)生態(tài)承載力的影響研究[J]. 干旱區(qū)資源與環(huán)境，2009，23（9）：24－28.Liu Ying, Yang Danqing, Yang Xiaohui. Effect of coal mining on town biologically productive capacity[J]. Journal of Arid Land Resources and Environment, 2009, 23(9): 24－28. (in Chinese with English abstract)

[28] 張思鋒，張立，張一恒．基于生態(tài)梯度風(fēng)險(xiǎn)評價(jià)方法的榆林煤炭開采區(qū)生態(tài)風(fēng)險(xiǎn)評[J]. 資源科學(xué)，2011，33(10)：1914－1923.Zhang Sifeng, Zhang Li, Zhang Yiheng. Ecological risk assessment of the Yulin coal mining area based on the PETAR method[J]. Resources Science, 2011, 33(10): 1914－1923. (in Chinese with English abstract)

[29] 楊軍，吳蔚，叢建輝，等. 煤炭開采與水資源利用保護(hù)的耦合關(guān)系研究——基于中國 10個(gè)煤炭資源豐裕省份的分析[J]. 生態(tài)經(jīng)濟(jì)，2016，32(2)：156－159.Yang Jun, Wu Wei, Cong Jianhui, et al. Coupled research of coal mining and protection and utilization of water Resource: based on the analysis of China 10 provinces in the abundance of coal resources[J], Ecological Economy, 2016,32(2): 156－159. (in Chinese with English abstract)

[30] 趙杰，趙士洞．參與性評估法在小尺度區(qū)域土地利用變化研究中的應(yīng)用[J]，資源科學(xué)，2003，25(5)：52－57．Zhao Jie, Zhao Shidong. Application on the participatory rural appraisal method in the research of land use change at local dimension[J]. Resources Science, 2003, 25(5): 52－57.(in Chinese with English abstract)

Determination of conflict between coal mining and water and soil resources utilization and land consolidation strategy of resource destruction in river basin

Liu Huifang1,2, Bi Rutian1※, Wen Bo2
(1. College of Resource & Environment, Shanxi Agricultural University, Taigu030801,China;2.College of Public Administration, Nanjing Agricultural University, Nanjing210095,China)

Quantifying the land-mining utilization conflicts will help to know the major limiting factors and improve the pertinence and effectiveness of land consolidation. Although the land-mining utilization and resource-environment effects and a lot of theoretical frameworks and qualification models have been studied to find ways to solve the resource use conflicts, it is not clear how the spatial heterogeneity of resources utilization conflicts based on resources destruction can be quantitatively visualized, and what kinds of targeted control strategies should be taken in various types and grades of resources utilization conflicts. This study carries out a case study on 47 villages in Changhe Basin of Zezhou County in Shanxi Province. Firstly,the evaluation index system was built from 3 aspects of land-mining resource system: quantity attribute, quality attribute and space attribute. The comprehensive index method and the landscape vector model were adopted to calculate the indices of land damage accumulation function, coal mining loss degree and water poverty for analyzing loss degree of land, water and coal resource. Then based on the accumulating effects of resource utilization conflicts, the sum evaluation model was used to measure the degree of land-mining utilization conflict and classify the conflicts. Thirdly, the consolidation measures were put forward based on different land-mining utilization conflicts degree and regional function division of land use. The results showed that: 1) The field survey and GIS (geographic information system) methods were adopted to overall evaluate the land-mining conflicts and illustrate 4 zones: the high conflict zone, the medium conflict zone, the low conflict zone and the conflict-free zone, which can reveal the regional differences of resources damage and interaction of resource utilization.Combined with land, water and coal resources damage assessment, the major limiting factors in different conflicts areas were revealed, which offered reference to land consolidation. 2) The areas of 4 zones are as follows: the high conflict zone (20.21%of the total area), the medium conflict zone (33.39 % of the total area), the low conflict zone (24.24% of the total area) and the non-conflict zone (22.16% of the total area). The high conflict zone is mainly located in the northwest and southwest parts of basin whose main features are coal being covered, serious land damage and water shortage; the medium conflict zone is distributed in west side of basin with serious land damage; the low conflict zone is mainly located in the northeast of basin,which is lacking of water resources; the non-conflict zone is distributed in east side of basin, which has no disturbance of coal mining and the farmland easy to be irrigated. 3) The zoning matrix of land use for land-mining utilization conflict control strategies was formulated based on conflict zones, which can be used to picture the land use control zone. The zoning matrix of land use included major ecological restoration area (Z1), major ecological restoration area (Z2), damaged land reclamation area (Z3), and cultivated land reserve area (Z4), water deficient restoration area (Z5), cultivated land and water reserve area(Z6), and cultivated land quality improvement area (Z7), which corresponded to forestry land (30.04%), forestry land (26.94%),general agricultural zone (28.21%) and basic farmland protection zone (Z4, Z5, Z6 and Z7 zone were above 37.77%). The results show that this zoning method is effective and can offer reference to land consolidation planning. The study reveals that the quantitative evaluation can not only provide the stereoscopic view of resources utilization conflicts, but also provide an operational approach to integrate land reclamation planning into industrial development planning by making targeted land use conflict prevention decisions, which helps to improve the effectiveness of land consolidation and land management in the mining region.

land use; land consolidation; water resources; river basin; combined influence index; land-destruction cumulative effect index; water poverty index; coal loss index

10.11975/j.issn.1002-6819.2017.14.033

F301.24; TD88

A

1002-6819(2017)-14-0238-12

劉慧芳，畢如田，文 博. 流域“地-礦”土地水資源利用沖突測度確定及土地整治策略[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(14)：238－249.

10.11975/j.issn.1002-6819.2017.14.033 http://www.tcsae.org

Liu Huifang, Bi Rutian, Wen Bo. Determination of conflict between coal mining and water and soil resources utilization and land consolidation strategy of resource destruction in river basin[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2017, 33(14): 238－249. (in Chinese with English abstract)

doi：10.11975/j.issn.1002-6819.2017.14.033 http://www.tcsae.org

2017-02-23

2017-07-05

國土資源部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)課題（201411007）

劉慧芳，女，山西長治人，講師，博士生，主要從事土地可持續(xù)利用方面的研究。太谷 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，030801。

Email：568326027@qq.com

※通信作者：畢如田，男，山西陽泉人，博士生導(dǎo)師。主要研究方向?yàn)橘Y源環(huán)境信息技術(shù)、土地復(fù)墾與規(guī)劃。太谷 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院，030801。Email：brt@sxau.edu.cn