基于生態(tài)適宜性分析的平原煤礦縣生態(tài)空間優(yōu)化

張 鵬，李世峰※，劉瑤瑤

（1.中國農(nóng)業(yè)大學水利與土木工程學院，北京100083；

2.中國農(nóng)業(yè)大學農(nóng)業(yè)農(nóng)村部設(shè)施農(nóng)業(yè)工程重點實驗室，北京100083）

0 引 言

煤炭是中國的基本能源，在中國國民經(jīng)濟中占據(jù)著極其重要的戰(zhàn)略位置[1]，煤礦城市在經(jīng)濟發(fā)展的過程中，不僅面臨快速的城鎮(zhèn)化導致城市周邊的生態(tài)空間不斷被蠶食，更需要面對煤礦開采造成了地形地貌破壞等系列問題[2]。平原礦區(qū)由于地下潛水位高，開采過程中不僅造成地面塌陷，更會致使塌陷區(qū)常年積水[3]，據(jù)估算，我國平原煤炭基地開采最終將造成沉陷面積31 813.33 km2，積水面積超過19 088 km2，煤礦塌陷地的治理受到當?shù)卣蛯W者的關(guān)注[4]。而當前一些煤礦城市在生態(tài)建設(shè)往往只關(guān)注塌陷區(qū)的生態(tài)修復，造成塌陷區(qū)生態(tài)修復良好，非塌陷區(qū)生態(tài)建設(shè)不良的局面[5-6]。因此，以塌陷區(qū)為核心，對區(qū)域整體的生態(tài)空間進行合理的評判，將對平原煤礦城市開展生態(tài)文明建設(shè)具有重要意義。

近年來，隨著生態(tài)文明建設(shè)的推進，生態(tài)空間研究迅速發(fā)展，生態(tài)空間是除建設(shè)用地以外的一切自然或人工的植物群落、山林水體及具有綠色潛能的空間等系列生態(tài)用地。學者對全國[7]、省域[8]、市域[9]、縣域[10-11]、甚至社區(qū)[12]等不同尺度的生態(tài)空間進行探討。隨著GIS 技術(shù)的不斷發(fā)展，生態(tài)適宜性分析、景觀格局分析等方法推動著生態(tài)空間的構(gòu)建。王義方等[13]從選取生態(tài)系統(tǒng)穩(wěn)定性、植被覆蓋率、距鐵路距離等共計11 個因子，對濟寧煤礦塌陷區(qū)進行生態(tài)適宜性評價，將研究區(qū)分為生態(tài)核心區(qū)、生態(tài)緩沖區(qū)、生態(tài)敏感區(qū)和綜合建設(shè)區(qū)等區(qū)域，為塌陷區(qū)的生態(tài)建設(shè)提供指導。蒙吉軍等[14]通過對貴陽市景觀類型進行評判選取生態(tài)源點，設(shè)置景觀阻力面運用最小累計阻力模型，構(gòu)建生態(tài)廊道，從而構(gòu)成系統(tǒng)的、有效的、連通的網(wǎng)狀空間結(jié)構(gòu)，連接生態(tài)要素，破解生境破碎化。

本研究以平原煤礦縣—山東省鄆城縣為研究對象，首先結(jié)合研究區(qū)特點，構(gòu)建適宜于平原煤礦縣的生態(tài)適宜性評價體系；其次在生態(tài)適宜性綜合分析的基礎(chǔ)上，引入最小累計阻力模型，不僅對研究區(qū)進行合理的綜合生態(tài)適宜性等級分區(qū)，確保生態(tài)、社會和經(jīng)濟效益的綜合最大化；同時完成生態(tài)網(wǎng)絡(luò)空間的分析，避免生態(tài)分區(qū)的割裂性。通過對不同重要性的分區(qū)結(jié)果實施不同的建設(shè)策略，并對生態(tài)網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化建議，從而為鄆城縣的生態(tài)修復建設(shè)和生態(tài)保護建設(shè)提供方向，為平原煤礦縣域的生態(tài)空間建設(shè)提供參考。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

鄆城縣位于115°40′～116°08′E，北緯35°19′～35°52′N，地處山東省菏澤市東北部，北鄰梁山縣，南接巨野縣、牡丹區(qū)，東與嘉祥縣接壤，西與鄄城縣毗鄰，南北長44.4 km，東西寬35.7 km，總面積1643 km2。

鄆城縣地處黃河沖積平原，全縣地勢平坦，是魯西煤炭基地的重要組成部分，共有趙樓煤礦、郭屯煤礦、彭莊煤礦、鄆城煤礦4 個煤礦處于開采之中，是典型平原煤礦縣（圖1）。自2007 年煤礦開采以來煤炭經(jīng)濟突飛猛進，推動著縣城鎮(zhèn)化的發(fā)展，2018年鄆城縣煤炭銷售收入達到74.38 億元，與此同時地形地貌破壞、生態(tài)空間萎縮、土地利用結(jié)構(gòu)破碎、土地塌陷積水等引發(fā)一系列的社會、經(jīng)濟、生態(tài)問題。截止2017年底，鄆城縣因采煤造成塌陷面積12 km2，其中，造成耕地絕產(chǎn)面積6.59 km2，減產(chǎn)面積5.41 km2，直接致使搬遷村莊27村、12145戶、39480人，波及104村莊居民生產(chǎn)生活，影響規(guī)模大，程度深。2018年鄆城縣為落實《菏澤市采煤塌陷地綜合治理工作實施方案》、《菏澤市人民政府關(guān)于深入推進國土資源節(jié)約用意見》，提出向資源節(jié)約集約利用型，綠色發(fā)展生態(tài)強縣進軍。構(gòu)建以煤礦塌陷區(qū)為生態(tài)修復核心，縣域系統(tǒng)性的生態(tài)空間成為鄆城縣亟待解決的問題。

圖1 鄆城縣礦區(qū)分布圖Fig.1 Yuncheng County mining area distribution map

1.2 數(shù)據(jù)來源

本研究數(shù)據(jù)如下：1）鄆城縣行政區(qū)劃、道路矢量數(shù)據(jù)來源于資源環(huán)境數(shù)據(jù)云平臺（http://www.resdc.cn/）2015年中國縣域行政邊界數(shù)據(jù)；2）2018 年Landsat TM 影像衛(wèi)星數(shù)據(jù)，數(shù)據(jù)分辨率為30 m，影像來源于地理空間數(shù)據(jù)云網(wǎng)站(http://www.gscloud.cn/) ，選擇云量為0 或接近于0 的遙感影像圖，以獲取植物覆蓋度、土地利用類型等數(shù)據(jù)；3）趙樓煤礦、郭屯煤礦、彭莊煤礦、鄆城煤礦4 大煤礦相關(guān)數(shù)據(jù)，由鄆城縣煤礦局2018年1月份統(tǒng)計提供，提取礦區(qū)拐點坐標、邊界范圍、塌陷深度、塌陷面積、工礦點位置等，以充分分析煤礦開采對鄆城縣生態(tài)適宜度的影響。

2 研究方法

生態(tài)適宜性分析由美國麥克哈格教授（McHarg）提出[16]，是生態(tài)規(guī)劃的核心，其基于生態(tài)規(guī)劃理論[17-18]、土地節(jié)約集約利用理論、空間分析理論、景觀生態(tài)學等理論基礎(chǔ)上，從生態(tài)系統(tǒng)與人文系統(tǒng)角度對給定資源開發(fā)利用方式的適宜程度進行評估，選擇有代表性因子進行量化，進而劃分適宜性等級，指導區(qū)域內(nèi)部自然開發(fā)利用、保護發(fā)展方向和環(huán)境功能區(qū)分化[19]。生態(tài)適宜性分析的基本步驟為：明確評價對象與范圍、選取評價指標、量化評價指標、賦予指標權(quán)重、確定評價單元、GIS數(shù)據(jù)庫建立、劃分單因子和綜合生態(tài)適宜性等級[20]。其中評價指標的選取與指標權(quán)重的賦予，直接關(guān)系到評價結(jié)果的科學性和保障性。

基于生態(tài)適宜性分析的生態(tài)網(wǎng)絡(luò)構(gòu)建是選擇部分最適宜區(qū)位生態(tài)源點，并以生態(tài)適宜性綜合指數(shù)倒數(shù)為阻力成本，使用最小累計阻力模型構(gòu)建生態(tài)廊道，以生態(tài)適宜性各等級分區(qū)作為生態(tài)基質(zhì)構(gòu)建生態(tài)空間格局。

2.1 生態(tài)適宜性分析指標的篩選

嚴格遵循指標選取的主導性和綜合性原則、針對性和對比性原則、應(yīng)用性和可行性原則，緊密結(jié)合鄆城縣的實際狀況和生態(tài)建設(shè)目標，即綜合鄆城縣煤炭局、農(nóng)業(yè)局和生態(tài)規(guī)劃專家的意見、參考相關(guān)研究文獻和《生態(tài)環(huán)境狀況評價技術(shù)規(guī)范》（以下簡稱《規(guī)范》）篩選研究區(qū)生態(tài)適宜性分析指標[21-22]。本文從自然生態(tài)、特征限制和人為干擾3個角度共計10個指標來探究研究區(qū)的生態(tài)適宜性（圖2）。

自然生態(tài)：鄆城縣地勢平坦，坡降1/5 000～1/10 000之間，高差9.0 m，故剔除表現(xiàn)地勢特征的高程、坡度和坡向等指標。選取植被覆蓋度、土地利用類型、生物豐度和距天然水體距離4 個因子指標反映全縣自然生態(tài)狀況。植被是生態(tài)系統(tǒng)的重要組成部分，植被覆蓋度是指單位面積內(nèi)植被的垂直投影面積所占比例，是衡量地表植被覆蓋反映當?shù)刂脖簧L狀況的重要指標，植被覆蓋度越高，表示該區(qū)域生態(tài)穩(wěn)定性越高，生態(tài)適宜性越好[23]。土地利用類型反映土地用途、性質(zhì)及其分布規(guī)律，是生態(tài)建設(shè)現(xiàn)狀重要表現(xiàn)，一般而言各土地類型對生態(tài)適宜性貢獻是：林地＞水域＞草地＞耕地＞建設(shè)用地[24]。生態(tài)本位意識的增強，尤其是2006 年起國家環(huán)保局試行《規(guī)范》以來，生物豐度因子作為權(quán)重最大的指標計算，公式（1）成為區(qū)域生態(tài)環(huán)境評價中的重要一環(huán)。河流、湖泊等天然水體，具有水源涵養(yǎng)、廢物處理和生物多樣性保護等重要的生態(tài)功能，對區(qū)域生態(tài)保護局域重要意義。

式中Abio為歸一化系指數(shù)，參考值為511.264；Sf、Sg、Sw、Sc、Scon、Sn和St分別為林地、草地、水域、耕地、建設(shè)用地、未利用地及區(qū)域土地總面積。

特征限制：是區(qū)域生態(tài)保護和生態(tài)修復最為關(guān)鍵的制約因素，本文選用塌陷深度、沉穩(wěn)性、距采礦點距離、距塌陷水體距離4 個因子表現(xiàn)煤礦開采對當?shù)厣鷳B(tài)適宜性的影響。煤礦開采造成塌陷是制約鄆城縣生態(tài)建設(shè)的重要因素。采礦點周邊堆放大量煤矸石，生產(chǎn)過程中“三廢”排放對周邊的土壤、大氣、水質(zhì)產(chǎn)生嚴重的負面影響。沉穩(wěn)性是地質(zhì)表現(xiàn)特征，在一定程度上影響后期的開發(fā)建設(shè)與生態(tài)建設(shè)，未沉穩(wěn)的地區(qū)是未來塌陷區(qū)的主要來源，應(yīng)著重關(guān)注。

人為干擾：選取距居民點距離和距道路距離2 個因子表現(xiàn)居民對鄆城縣的生態(tài)擾動，區(qū)域距離道路和居民點越近，生態(tài)環(huán)境受人類干擾越大。

圖2 生態(tài)適宜性分析指標層次Fig.2 Ecological suitability analysis index level

2.2 生態(tài)適宜性分析指標的量化

生態(tài)適宜性分析是定量表述方法，由于各指標的具有不同的量綱或是定性描述，不便于計算，因此需要對指標無量化處理消除各指標量綱的影響。參考相關(guān)文獻，各指標及適宜性分值見表1。

其中植被覆蓋度和生物豐度使用線性歸一化方法，將是定性描述的因子定量化計算后將原始數(shù)據(jù)等比例縮放至[0,1]范圍內(nèi)。其計算公式為

式中Ai為歸一化數(shù)據(jù)，Xi為第i 種評價因子，Xmax為Xi的最大值，Xmin為Xi的最小值。

土地利用類型利用ArcGIS10.5 平臺上對2018 年Landsat TM遙感影像數(shù)據(jù)進行幾何校正、裁剪，采用監(jiān)督分類最大似然法，對遙感影像圖進行目視解譯，其中塌陷水體、工礦用地來源于鄆城縣煤炭局CAD 文件，塌陷水域為常年積水大于50 mm 的區(qū)域。距天然水體距離、距采礦點距離、距塌陷水體距離、距居民點距離，距道路距離是使用ArcGIS 中“歐式距離”工具得到等級圖層。塌陷深度和沉穩(wěn)性是將CAD 數(shù)據(jù)導入ArcGIS 中使用“重分類”工具繪制得出。

2.3 生態(tài)適宜性分析指標的權(quán)重

本文采用層次分析法（analytic hierarchy process，AHP）對各指標賦予權(quán)重，是目前最為常用的因子權(quán)重分析方法，通過構(gòu)建有層次化的結(jié)構(gòu)模型，各指標兩兩比較構(gòu)建判別矩陣，得到科學有效的決策。

通過對鄆城縣煤炭局、農(nóng)業(yè)局和生態(tài)規(guī)劃專家的訪談獲取他們對各個指標之間的相對重要性的對比，構(gòu)建生態(tài)適宜性判別矩陣。使用yaahp 軟件群決策工具對判斷矩陣進行計算，對比各判斷矩陣進行一致性檢驗，并計算各層次組合權(quán)重值。結(jié)果如下，判斷矩陣均滿足一致性，各指標權(quán)重如表2所示。

表1 指標適宜性分值Table 1 Score value of suitability of evaluation index

2.4 生態(tài)適宜性分析單元劃分方法

評價單元是生態(tài)適宜性分析最基本空間單元，生態(tài)適宜性分析結(jié)果是承載于評價單元之上，其決定生態(tài)適宜性分析結(jié)果精度以及成果的應(yīng)用性。目前國內(nèi)外主要有：基于土壤類型、地貌類型圖斑等單元劃分，將性質(zhì)相對均一的地塊作為評價單元；基于多因素疊加形成的同質(zhì)單元劃分，將疊加后最小的公共斑塊作為評價單元；基于GIS柵格格網(wǎng)進行劃分。本文采用第三種方法進行劃分，采用分辨率為30 m 的柵格圖像，其優(yōu)點為便于表達空間連續(xù)性表達和便于評價模型的運算。

表2 生態(tài)適宜性分析指標權(quán)重Table 2 Ecological suitability analysis index weight

2.5 生態(tài)適宜性分析模型確定方法

本文采用麥克哈格教授（McHarg）提出生態(tài)因子疊置法，采用加權(quán)求和的方式來計算生態(tài)適宜性綜合指數(shù)，其計算公式為

式中Q 表示區(qū)域生態(tài)適宜度的評價結(jié)果；eji為區(qū)內(nèi)第j個生態(tài)因子在評價第i個指標的適宜度值；wji為因子j對于指標i的權(quán)重；ui為指標i在適宜度總體評價中的權(quán)重。

2.6 生態(tài)適宜性等級劃分

生態(tài)適宜性的等級劃分，是將生態(tài)適宜性指數(shù)按照自然斷點法，參考相關(guān)文獻的基礎(chǔ)上將研究區(qū)分為五個等級分區(qū)，分別為最不適宜區(qū)、不適宜區(qū)、勉強適宜區(qū)、適宜區(qū)、最適宜區(qū)，以指導研究區(qū)的生態(tài)建設(shè)。具體分區(qū)標準見表3。

表3 生態(tài)適宜性分區(qū)劃分標準Table 3 standard for ecological suitability division

2.7 MCR模型分析生態(tài)網(wǎng)絡(luò)空間

最小累計阻力模型（minimal cumulative resistance，MCR），最初常用于動物生境保護等生態(tài)領(lǐng)域[25]，引入我國后廣泛應(yīng)用于景觀格局分析[26]。本文借助MCR模型，以生態(tài)適宜性綜合指數(shù)倒數(shù)為生態(tài)阻力面，模擬生態(tài)最適宜區(qū)之間的最短路徑，構(gòu)建生態(tài)網(wǎng)絡(luò)空間。

式中MCR 為最小累積阻力；f是一未知負函數(shù)，為最小累積阻力與生態(tài)適宜性的負相關(guān)關(guān)系；min 為某單元對不同源取累積阻力最小值；Dij為從源j 到源i 的空間距離；Ri為源i對運動過程的阻力系數(shù)。

3 結(jié)果與分析

3.1 生態(tài)適宜性綜合分區(qū)

根據(jù)各指標的等級分區(qū)與指標權(quán)重，利用Arc-GIS10.5軟件“柵格計算器”進行疊加分析，計算得到鄆城縣生態(tài)適宜性綜合指數(shù)得分，區(qū)域生態(tài)適宜度綜合分值在1.66～4.88之間，平均值為3.19。從圖3中可以看出，生態(tài)適宜性較高的區(qū)域集中于縣域北側(cè)黃河周邊，依次向南部遞減；生態(tài)適宜性較低的區(qū)域，集中于城鎮(zhèn)建設(shè)用地，塌陷水域和塌陷深度較深的地域，并向周邊遞增。整體分布趨勢與縣域現(xiàn)狀相符。

圖3 研究區(qū)生態(tài)適宜性指數(shù)Fig.3 Ecological suitability index of study area

結(jié)合生態(tài)適宜的分區(qū)標準，計算各等級面積（表4）。導入鄉(xiāng)鎮(zhèn)矢量邊界文件，將生態(tài)適宜性分區(qū)結(jié)果細化到鄆城縣各鄉(xiāng)鎮(zhèn)（圖4）。

1）最不適宜區(qū)：面積90.66 km2，占據(jù)全縣5.52%，主要位于鄆城縣中東部和南部。生態(tài)最不適宜區(qū)主要集中在鄆州街19.60%、南趙樓鎮(zhèn)19.05%、張營鎮(zhèn)16.89%、丁里長鎮(zhèn)10.86%、郭屯鎮(zhèn)7.17%和潘渡鎮(zhèn)5.86%，6 鎮(zhèn)生態(tài)最不適宜區(qū)面積共計74.05 km2，占據(jù)最不適宜區(qū)面積的81.68%。結(jié)合土地利用分析可以發(fā)現(xiàn)，鄆州街是鄆城縣城所在地，城鎮(zhèn)化建設(shè)高、人口密度大、道路密集是生態(tài)不適宜度高的主要原因，生態(tài)最不適宜面積占據(jù)自身面積28.29%。其余5鎮(zhèn)生態(tài)最不適宜面積偏高的原因主要是塌陷水域和塌陷面積的存在。南趙樓鎮(zhèn)和張營鎮(zhèn)最不適宜分區(qū)面積均超過了10 km2，僅次于鄆州街，面積分別為17.27 km2和15.31 km2，其中南趙樓鎮(zhèn)最不適宜區(qū)面積占據(jù)鄉(xiāng)鎮(zhèn)面積的29.84%，生態(tài)環(huán)境狀況極為嚴峻。

2）不適宜區(qū)：面積247.99 km2，比重16.73%，多位于最不適宜區(qū)與勉強適宜區(qū)的過渡地帶，呈“T”字狀主要分布在潘渡鎮(zhèn)12.71%、張營鎮(zhèn)12.20%、鄆州街11.39%、南趙樓鎮(zhèn)10.80%、楊莊集9.81%、唐廟鎮(zhèn)8.83%、唐塔街8.29%。這些區(qū)域受城鎮(zhèn)化建設(shè)和煤礦開采的雙重影響，生態(tài)環(huán)境處于不適宜狀態(tài)。

3）勉強適宜區(qū)：面積643.42 km2，比重39.15%，面積比重最大，分布最廣。整體呈現(xiàn)片狀分布在全縣各個鄉(xiāng)鎮(zhèn)。其中比重稍重的鄉(xiāng)鎮(zhèn)為黃安鎮(zhèn)7.90%、楊莊集7.33%、隨官屯7.29%、程屯鎮(zhèn)6.87%、武安鎮(zhèn)6.54%、雙橋鎮(zhèn)6.45%。

4）適宜區(qū)：面積530.64 km2，比重32.29%，僅次于勉強適宜區(qū)。呈團塊狀集中在鄆城縣西北部區(qū)域，李集鄉(xiāng)、侯咽集、玉皇廟、武安鎮(zhèn)、張魯集、黃集鄉(xiāng)生態(tài)適宜區(qū)面積比重均超過鄉(xiāng)鎮(zhèn)面積的50%。此區(qū)域城鎮(zhèn)建設(shè)面積較少，多為耕地。

5）最適宜區(qū)：面積103.66 km2，比重6.31%，呈斑塊狀鑲嵌在適宜區(qū)之中，集中分布在玉皇廟鎮(zhèn)、侯咽集鎮(zhèn)、李集鄉(xiāng)、黃集鄉(xiāng)、張魯集鄉(xiāng)，5 鎮(zhèn)生態(tài)最適宜區(qū)面積共計77.67 km2，占據(jù)最適宜區(qū)面積比重的74.93%。其中玉皇廟鎮(zhèn)生態(tài)適宜區(qū)面積最大17.91 km2，占據(jù)鄉(xiāng)鎮(zhèn)面積19.91%，其北與侯咽集鎮(zhèn)，南與水堡鄉(xiāng)，構(gòu)成較為集中的生態(tài)最適宜區(qū)條帶。

表4 生態(tài)適宜性分區(qū)面積及比重Table 4 Area and proportion of ecological suitability division

3.2 生態(tài)網(wǎng)絡(luò)結(jié)構(gòu)

研究區(qū)的生態(tài)源點選擇原則為各鄉(xiāng)鎮(zhèn)中最適宜區(qū)斑塊面積最大為生態(tài)源點，共選取生態(tài)源點14處。在此基礎(chǔ)之上，進行生態(tài)廊道的識別，通過鄆城縣生態(tài)適宜性評價指數(shù)的倒數(shù)構(gòu)成生態(tài)復合阻力面，生成每個生態(tài)源點與其他生態(tài)源點的最短成本路徑并疊加，共計43 條，與生態(tài)適宜性等級分區(qū)構(gòu)建了鄆城縣生態(tài)空間格局如圖5所示。

鄆城縣生態(tài)網(wǎng)絡(luò)特征“點—線—面”框架初現(xiàn)，基本實現(xiàn)全縣連通，但分布不均勻，主要集中于西北部區(qū)域，南部和中部相對欠缺，南北連通優(yōu)于東西連通。由于生態(tài)源點的呈現(xiàn)不均勻狀態(tài)分布，生態(tài)網(wǎng)絡(luò)的疏密也不均。其中唐塔街道和楊莊集鎮(zhèn)的生態(tài)源點壓力極大，其面積小，但卻是溝通中部東西向，北部南北向的極其重要的廊道，一旦無新的生態(tài)源地補充或自身消失，那么鄆城縣的生態(tài)網(wǎng)絡(luò)將會受到嚴重的損害。當前不均勻的網(wǎng)絡(luò)空間結(jié)構(gòu)，不利于提升縣域的生態(tài)水平。

圖4 研究區(qū)生態(tài)適宜性分區(qū)Fig.4 Ecological suitability division of study area

圖5 研究區(qū)生態(tài)網(wǎng)絡(luò)結(jié)構(gòu)Fig.5 Ecological network structure of study area

3.3 生態(tài)空間優(yōu)化措施

通過生態(tài)適宜性和生態(tài)網(wǎng)絡(luò)分析，可以看出由于礦區(qū)塌陷帶來的系列問題致使鄆城縣東西部生態(tài)適宜性、生態(tài)網(wǎng)絡(luò)空間密集度差距較大，中部生態(tài)斑塊缺失，東西連通不佳。當前鄆城縣應(yīng)以最不適宜區(qū)的生態(tài)修復和最適宜區(qū)的保護和拓展為主導，通過對不同等級分區(qū)針對性的生態(tài)工程建設(shè)和生態(tài)網(wǎng)絡(luò)“點—線—面”的優(yōu)化，從而打造層次豐富、功能多樣的復合型網(wǎng)絡(luò)化生態(tài)支撐體系，改進區(qū)域的生態(tài)環(huán)境，完成向煤礦生態(tài)縣的轉(zhuǎn)型。

3.3.1 生態(tài)適宜性分區(qū)優(yōu)化

1）最不適宜區(qū)—生態(tài)修復核心區(qū)

最不適宜區(qū)應(yīng)作為鄆城縣生態(tài)修復核心區(qū)域，是當前生態(tài)修復工程投資建設(shè)的前沿區(qū)域。除城區(qū)外，主要由塌陷水域、塌陷程度較深的區(qū)域組成。研究區(qū)地處黃河沖擊平原，地下水位高，土壤以潮土為主，含水量接近飽和狀態(tài)，當前大面積的常年性積水區(qū)和從四周向中心逐級加深地勢等條件為人工濕地的構(gòu)建提供了基礎(chǔ)條件。鄆城縣應(yīng)以張營鎮(zhèn)、南趙樓的人工濕地工程為先導，打造以兩鎮(zhèn)為生態(tài)修復中心，以丁里長鎮(zhèn)、郭屯鎮(zhèn)、潘渡鎮(zhèn)塌陷積水區(qū)域治理為后續(xù)，構(gòu)建人工濕地走廊。通土壤、植物、動物、水域、微生物的協(xié)調(diào)作用[27]，對塌陷區(qū)水質(zhì)環(huán)境進行凈化治理[28]，不僅變劣勢為優(yōu)勢，更可以成為重要的生態(tài)源地，使得縣域東部區(qū)域的生態(tài)網(wǎng)絡(luò)空間密集型增大；而且可以發(fā)展如濕地旅游、濕地農(nóng)業(yè)、濕地漁業(yè)等多種形式的生態(tài)產(chǎn)業(yè)，增加農(nóng)民收入。

2）不適宜區(qū)—生態(tài)修復過渡區(qū)

不適宜區(qū)應(yīng)作為鄆城縣生態(tài)修復過渡區(qū)域，是生態(tài)修復工程關(guān)注的重點。主要由距離塌陷區(qū)域較近城鎮(zhèn)化水平較高區(qū)域組成。其一應(yīng)連同城區(qū)，通過公共綠地、道路綠地、防護綠地、生產(chǎn)綠地、居住區(qū)綠地等多種形式造林種草增加城鎮(zhèn)居民游憩空間的同時增加中部城鎮(zhèn)建設(shè)區(qū)域生態(tài)斑塊，以提高生態(tài)網(wǎng)絡(luò)穩(wěn)定性；其二做好防災基礎(chǔ)設(shè)施建設(shè)，避免煤礦開采造成新塌陷區(qū)大面積的形成對生態(tài)系統(tǒng)造成劇烈影響。

3）勉強適宜區(qū)—生態(tài)綜合建設(shè)區(qū)

勉強適宜區(qū)應(yīng)作為鄆城縣生態(tài)綜合建設(shè)區(qū)域，當前該區(qū)域整體城鎮(zhèn)化建設(shè)不高，也是鄆城縣城鎮(zhèn)化建設(shè)過程中的理想拓展區(qū)域。在建設(shè)規(guī)劃過程中應(yīng)以生態(tài)化、集約化為建設(shè)原則，避免城鎮(zhèn)建設(shè)散亂對生態(tài)用地的蠶食。

4）適宜區(qū)—生態(tài)保護緩沖區(qū)

適宜區(qū)應(yīng)作為鄆城縣生態(tài)保護緩沖區(qū)域，其主要是由距離最適宜區(qū)較近的耕地組成。應(yīng)當保持現(xiàn)有的良好的農(nóng)田生態(tài)景觀穩(wěn)定性和耕地數(shù)量，改善土壤養(yǎng)分和結(jié)構(gòu)，維護區(qū)域生態(tài)安全和糧食生產(chǎn)安全。以綠色化、有機化為生產(chǎn)過程技術(shù)為基礎(chǔ)，發(fā)展綠色農(nóng)業(yè)，由農(nóng)業(yè)大縣向生態(tài)農(nóng)業(yè)強縣邁進，尤其是對靠近生態(tài)最適宜區(qū)的區(qū)域應(yīng)著重進行生態(tài)改造以起到對核心區(qū)的生態(tài)緩沖作用。

5）最適宜區(qū)—生態(tài)保護核心區(qū)

最適宜區(qū)應(yīng)作為鄆城縣生態(tài)保護核心區(qū)域，大面積集中于臨近北部黃河區(qū)域，天然水體、草地較為集中，生態(tài)系統(tǒng)較其他區(qū)域復雜，生物豐度較高，受人類活動干擾相對較小。但近年來，由于耕地開墾，黃河沿岸草地退化嚴重，林地面積逐漸減少，這也是鄆城縣景觀格局多樣性無法提高、生態(tài)網(wǎng)絡(luò)復雜程度不強的重要原因29]。最適宜區(qū)應(yīng)重點保護，以涵養(yǎng)水源、生物多樣性保護為主，嚴格禁止產(chǎn)生嚴重污染環(huán)境工程建設(shè)；同時打造生態(tài)河道，尤其是構(gòu)建黃河水域生態(tài)帶，形成水—草—林相互協(xié)調(diào)的更為復雜穩(wěn)定的生態(tài)系統(tǒng)。

3.3.2 生態(tài)空間網(wǎng)絡(luò)優(yōu)化

1）點：重點保護生態(tài)源地、強化生態(tài)小斑塊建設(shè)

生態(tài)源地是研究區(qū)的生態(tài)核心，應(yīng)不斷提高生態(tài)源地的生態(tài)功能。增加對玉皇廟鎮(zhèn)、侯咽集鎮(zhèn)、李集鄉(xiāng)、黃集鄉(xiāng)、張魯集鄉(xiāng)和水堡鄉(xiāng)的生態(tài)源地保護力度，6 鎮(zhèn)在鄆城縣占據(jù)著重要的生態(tài)位置，尤其是李集鄉(xiāng)應(yīng)該發(fā)揮極核作用，成為生態(tài)安全格局的核心地和示范地；拓展唐塔街道和楊莊集鎮(zhèn)生態(tài)源地的面積，緩解生態(tài)源地壓力；通過對張營鎮(zhèn)、南趙樓、丁里長鎮(zhèn)、郭屯鎮(zhèn)、潘渡鎮(zhèn)的生態(tài)修復治理增加生態(tài)源地數(shù)量，提高中東部區(qū)域生態(tài)網(wǎng)絡(luò)的復雜性。

2）線：設(shè)置生態(tài)廊道緩沖區(qū)

當前鄆城東部、南部區(qū)域生態(tài)廊道缺少復雜性和連通性，極易收到損害造成生態(tài)廊道的斷裂。為提高縣域生態(tài)空間的整體連通性，建議對北部區(qū)域的廊道設(shè)置生態(tài)緩沖區(qū)，南部東部區(qū)域廊道進行改造增加起墊腳石作用的生態(tài)小斑塊，從而提升廊道之間的聯(lián)系和周邊區(qū)域的交互作用，增加廊道穩(wěn)定性；中部區(qū)域生態(tài)廊道較少，在整體規(guī)劃的過程中可以先構(gòu)建小生態(tài)斑塊，形成“點”基礎(chǔ)。通過對現(xiàn)有廊道資料分析，對生態(tài)廊道設(shè)置1000米左右的生態(tài)廊道寬度[30-31]，嚴禁城鎮(zhèn)化建設(shè)對生態(tài)廊道的破壞。

3）面：全面優(yōu)化生態(tài)基質(zhì)

鄆城縣生態(tài)提升的關(guān)鍵是生態(tài)最不適宜區(qū)的修復和生態(tài)適宜區(qū)的保護，通過對各個生態(tài)適宜分區(qū)的生態(tài)建設(shè)，如通過人工濕地建設(shè)、綠地工程建設(shè)、生態(tài)河道系統(tǒng)全面的提升生態(tài)要素的質(zhì)量，為生態(tài)空間結(jié)構(gòu)的穩(wěn)定性奠定基礎(chǔ)。

4 結(jié)論與討論

煤礦區(qū)塌陷為難點的生態(tài)問題已成為制約平原煤礦地區(qū)社會、經(jīng)濟、生態(tài)可持續(xù)發(fā)展的重要因素。本文在生態(tài)適宜性分析的基礎(chǔ)上引入最小累計阻力模型，通過生態(tài)分區(qū)與生態(tài)網(wǎng)絡(luò)優(yōu)化縣域生態(tài)空間。主要結(jié)論如下：

1）從自然生態(tài)、特征限制、人為干擾3個角度選取10個指標，構(gòu)建平原煤礦地區(qū)生態(tài)適宜性評價體系，將研究區(qū)分為最不適宜區(qū)、不適宜區(qū)、勉強適宜區(qū)、適宜區(qū)、和最適宜區(qū)5 個分區(qū)，有助于對縣域生態(tài)本地把握。當前各分區(qū)面積呈現(xiàn)正態(tài)分布但分布不均勻，最適宜區(qū)集中于縣域北側(cè)，最不適宜區(qū)集中于中南部的塌陷水域和塌陷深度較深的地區(qū)。

2）利用選取各鄉(xiāng)鎮(zhèn)最大生態(tài)最適宜性斑塊為生態(tài)源點，以生態(tài)適宜性綜合指數(shù)倒數(shù)為生態(tài)成本，最小累計阻力模型計算鄆城縣生態(tài)廊道，構(gòu)建生態(tài)網(wǎng)絡(luò)。目前生態(tài)網(wǎng)絡(luò)分布主要集中于西北部區(qū)域，南部和中部相對欠缺，南北連通優(yōu)于東西連通。

3）針對鄆城縣生態(tài)現(xiàn)狀，文章提出應(yīng)以最不適宜區(qū)為生態(tài)修復核心，最適宜區(qū)為生態(tài)保護核心，通過人工濕地建設(shè)、綠地工程建設(shè)、生態(tài)河道建設(shè)等生態(tài)措施對各等級分區(qū)進行優(yōu)化；通過對生態(tài)源地的保護和小生態(tài)斑塊建設(shè)，生態(tài)廊道緩沖區(qū)建設(shè)，生態(tài)基質(zhì)分區(qū)優(yōu)化等措施提升生態(tài)網(wǎng)絡(luò)結(jié)構(gòu)的穩(wěn)定性。從而為區(qū)域的生態(tài)空間優(yōu)化提供依據(jù)。

生態(tài)空間是各種生態(tài)過程以及人類活動在不同尺度上作用的結(jié)果。本文只對2018 年的縣域的生態(tài)適宜性和生態(tài)網(wǎng)絡(luò)進行判別，是一種靜態(tài)分析。不同尺度、多時期動態(tài)分析將會是下一步工作的重點。