周 敏，齊增湘*，呂婧瑋，劉 鑫

(1.南華大學(xué) 建筑學(xué)院，湖南 衡陽 421001；2.湖南省健康城市營造工程技術(shù)研究中心，湖南 衡陽 421001；3.衡陽市科技局 生態(tài)型區(qū)域/城市規(guī)劃與管理重點實驗室，湖南 衡陽 421001)

工礦城鎮(zhèn)不僅在我國經(jīng)濟發(fā)展中做出了巨大的貢獻(xiàn)[1]，而且對加快現(xiàn)代化建設(shè)和推動工業(yè)化進(jìn)程起了有力的支撐作用。但由于缺乏統(tǒng)籌規(guī)劃和對資源粗放型的開發(fā)，致使城市建設(shè)用地與非城市用地發(fā)展不協(xié)調(diào)，從而帶來了生態(tài)環(huán)境破壞、資源過度開發(fā)、產(chǎn)業(yè)衰退等城市問題[2]，其中廢棄工礦區(qū)成為了制約城市化進(jìn)程與可持續(xù)發(fā)展的主導(dǎo)因素。1977年生態(tài)修復(fù)被作為生態(tài)學(xué)分支進(jìn)行研究后，其逐漸被認(rèn)為對促進(jìn)生態(tài)文明建設(shè)、廢棄礦區(qū)土地再利用、城鎮(zhèn)可持續(xù)發(fā)展具有重要支撐作用[3]。尤其在黨的“十八大”全面深化改革和“十九大”報告的重大戰(zhàn)略部署下，“生態(tài)修復(fù)，城市修補”形成了城市整體小干預(yù)、漸進(jìn)性的治理模式[4-5]。因此，歷史遺留的廢棄工礦區(qū)及生態(tài)環(huán)境問題在新的城市治理模式下進(jìn)行生態(tài)修復(fù)以及土地再利用具有迫切的現(xiàn)實意義。

在當(dāng)前生態(tài)文明發(fā)展下，生態(tài)安全格局的構(gòu)建被認(rèn)為是實現(xiàn)區(qū)域生態(tài)安全與土地可持續(xù)發(fā)展的基本保障和必要途徑[6]。其理論通過識別和保護(hù)關(guān)鍵的局部、點、空間關(guān)系來維護(hù)區(qū)域生態(tài)系統(tǒng)的穩(wěn)定發(fā)展[7]。R.T.T.Forman[8]首先系統(tǒng)性總結(jié)了景觀格局的優(yōu)化方法，俞孔堅[7]則是在此基礎(chǔ)上提出了景觀生態(tài)安全格局的概念，此后成為國內(nèi)的研究熱點，引發(fā)了學(xué)者在城市的景觀安全格局[9]、風(fēng)景名勝區(qū)[10]、地質(zhì)災(zāi)害區(qū)[11]、農(nóng)村居民點[12]等方面都做了廣泛研究。目前，除了少量基于“格局—過程”的區(qū)域生態(tài)安全格局構(gòu)建模式[13]外，“源地識別—阻力面構(gòu)建—提取廊道”的方法已成為構(gòu)建生態(tài)安全格局的基本范式[11]。近些年， J.P.knaapen[14]提出的最小累積阻力模型(minimum cumulative resistance model，MCR)成為了主要的模型方法研究之一，其針對于景觀對空間運動過程的阻礙作用展現(xiàn)出良好的模擬效果，相比傳統(tǒng)模型能更好地表達(dá)景觀安全格局與生態(tài)過程的相互關(guān)系。目前景觀安全格局研究中存有不足，首先，在研究范圍方面，雖已形成鎮(zhèn)域、縣域、市域等多個不同尺度和區(qū)域的研究，但針對典型的工礦區(qū)城鎮(zhèn)研究鮮少。其次，在源的構(gòu)建與識別方面，局限于單一源，尚少考慮“多源”綜合下所表現(xiàn)出的內(nèi)外空間分布差異，源的識別方法上以直接識別為主，缺乏客觀性。最后，在構(gòu)建生態(tài)阻力面上僅考慮景觀類型的空間位置卻忽略景觀類型間的相互影響，因此，很難達(dá)到滿意的研究效果，降低了研究方法的實踐價值。

鑒于當(dāng)前研究不足，以水口山鎮(zhèn)為研究區(qū)，綜合考慮城市和生態(tài)兩者的擴張過程，基于遙感和ArcGIS軟件，選取土壤重金屬綜合污染指數(shù)、土地類型、歸一化植被指數(shù)等9個因子構(gòu)建阻力面，模擬城市空間和生態(tài)源地的水平競合過程，量化工礦區(qū)構(gòu)建生態(tài)用地的潛力，提出生態(tài)修復(fù)分區(qū)和治理對策，為城市的空間布局及生態(tài)治理提供技術(shù)支持。

1 研究區(qū)與研究數(shù)據(jù)

1.1 研究區(qū)概況

水口山鎮(zhèn)位于常寧市東北角，北傍湘江、東依舂陵河、南與煙洲鎮(zhèn)相連、西與蓬塘鄉(xiāng)、柏坊鎮(zhèn)接址，總面積86.19 km2，區(qū)內(nèi)以平原、低山丘陵區(qū)地貌為主。氣候?qū)儆趤啛釒О霛駶櫦撅L(fēng)氣候，年平均氣溫17.9℃，年平均風(fēng)速3.5 m/s，年平均降雨量1 394.9 mm，年平均日照時數(shù)122.9 h。鎮(zhèn)內(nèi)蘊藏著豐富的鉛、鋅、金等20多種礦產(chǎn)資源，其鉛鋅儲量居全國第1位，素有“世界鉛都”“有色金屬之鄉(xiāng)“的美譽，在國家經(jīng)濟發(fā)展中占居重要的地位。100多a的開采歷史背景下，私挖亂采的活動頻繁，部分礦產(chǎn)資源開發(fā)利用至今都保持著粗放的開發(fā)模式，未經(jīng)處理的污水排放，對生態(tài)環(huán)境造成嚴(yán)重的破壞，裸露的礦區(qū)及淘汰的工業(yè)基地達(dá)到588.69 hm2，嚴(yán)重制約了經(jīng)濟的可持續(xù)發(fā)展(圖1)。

圖1 工礦區(qū)分布示意圖

1.2 數(shù)據(jù)獲取和處理

采用的數(shù)據(jù)主要包括研究區(qū)2018年3月Landsat 8 OLI遙感影像數(shù)據(jù)、DEM 30 m分辨率的數(shù)字高程數(shù)據(jù)(www.gscloud.cn)、土地利用類型數(shù)據(jù)利用遙感影像，結(jié)合常寧市水口山鎮(zhèn)(工業(yè)新城)總體規(guī)劃(2015-2030)及實地調(diào)查和定位，在ENVI5.4軟件中目視解譯為水田、旱地、有林地、灌木林、疏林地、其他林地、高覆蓋度草地、河渠、水庫坑塘、城鎮(zhèn)用地、農(nóng)村居民點、其他建設(shè)用地共12種用地類型。土壤重金屬綜合污染指數(shù)于2020年6月采用隨機與網(wǎng)格布點(500 m×500 m)結(jié)合的方法，嚴(yán)格按《土壤環(huán)境監(jiān)測技術(shù)規(guī)范》(HJ/T166-2006)完成實地采樣，并用X-MET 8000手持式分析儀器測定土壤中Pb、Cu、Zn、As、Cd和Hg 6種重金屬含量，運用SPSS軟件對實測數(shù)據(jù)進(jìn)行正態(tài)檢驗及處理，最后在ArcGIS中采用地統(tǒng)計學(xué)的Kriging法插值生成。NDVI在ENVI中通過波段公式獲得。數(shù)據(jù)柵格大小全部為30 m×30 m，統(tǒng)一投影坐標(biāo)為WGS 1984 UTM ZONE 49N。

2 研究方法

2.1 歸一化植被指數(shù)

為了消除部分輻射誤差與大氣的影響，用ENVI對多光譜遙感影像進(jìn)行輻射定標(biāo)和大氣校正后提取歸一化植被指數(shù)(normalized difference vegetation index，NDVI)[15]。

(1)

式中，NIR為紅外光譜反射值，R為紅光反射值。

2.2 土壤重金屬綜合污染指數(shù)

根據(jù)《土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風(fēng)險管控標(biāo)準(zhǔn)》(GB 15618-2018)和《土壤環(huán)境監(jiān)測技術(shù)規(guī)范》(HJ/T166-2006)[16-17]，采用內(nèi)梅羅法計算土壤重金屬綜合污染指數(shù)。

(2)

式中，PImean和PImax分別為平均單項污染指數(shù)和最大單項污染指數(shù)。

2.3 水平競合模型構(gòu)建

城市空間演變可看作是城市生態(tài)用地與城市建設(shè)用地兩者相互對弈的過程。為了實現(xiàn)生態(tài)服務(wù)功能與社會經(jīng)濟效益兩者的最大化，城市的生態(tài)用地和建設(shè)用地兩者同樣需要擴張，將工礦區(qū)生態(tài)修復(fù)結(jié)果看成是城市建設(shè)與生態(tài)用地競合的過程。引入MCR模型形成建設(shè)用地和生態(tài)擴張的累積阻力面，并以兩者的差值量化模擬城市空間和生態(tài)用地的競合過程，以此來評價工礦區(qū)的生態(tài)修復(fù)潛力。

2.3.1 最小累積阻力模型 最小累積阻力模型(MCR)是耗費距離模型的衍生應(yīng)用，其實質(zhì)上是通過“源”克服其擴散或遷徙“阻力”的耗費，反映的是一種可達(dá)性[18]。模型考慮源、距離和阻力面特征三個要素，其公式[19]：

(3)

式中，MCR為最小累積模型；f為一個未知的正函數(shù)，反映的是空間中任何一點的最小阻力與其到所有源的距離和景觀基面特征的正相關(guān)關(guān)系。Dij代表從空間中某單元i到源j的實地距離，Ri是某一景觀i的阻力值。

本研究借鑒類似方法，將工礦區(qū)修復(fù)為城市建設(shè)和生態(tài)用地兩種用途，由此建立兩種用途的最小累積阻力差值模型，模擬工礦區(qū)生態(tài)修復(fù)的潛力并劃分優(yōu)先級：

ΔMCR=MCR生態(tài)用地-MCR建設(shè)用地#

(4)

式中，MCR生態(tài)用地為生態(tài)用地綜合最小累積阻力表面，MCR建設(shè)用地為城市擴張用地綜合最小累積阻力表面。ΔMCR為生態(tài)用地與城市用地兩者綜合最小累積阻力表面的差值結(jié)果。

3 結(jié)果與分析

3.1 源的選取

在城市空間擴張源地識別方面，城市集中建設(shè)的區(qū)域是其重點發(fā)展的對象，可作為城市空間擴張過程的源地?；?015年水口山鎮(zhèn)土地利用變更調(diào)查與實地調(diào)研，提取了集中分布的中心區(qū)及大型交通用地，總面積達(dá)到355.95 hm2(圖2)。

在生態(tài)擴張源地識別方面，源地為事物或事物向外起點和基地[20]，因此必須是生態(tài)服務(wù)功能較強、生境質(zhì)量服務(wù)較高的區(qū)域[19]。借鑒相關(guān)研究[21]，采用生態(tài)服務(wù)價值評價與生境質(zhì)量評價等權(quán)疊加的方法，得到水口山鎮(zhèn)生態(tài)用地綜合評價結(jié)果，再用自然斷點法選取等級最高的作為生態(tài)用地擴張的源地，面積為3 880 hm2，結(jié)構(gòu)上看源地分布均勻(圖2)。

3.2 水平競合分析

生態(tài)過程中，源地是通過克服阻力對空間水平方向的競爭性控制和覆蓋過程來建立空間聯(lián)系[20]。該阻力是通過不同因子構(gòu)成的阻力面來實現(xiàn)的，其反映了源在空間運行時受到阻力強度空間的分布與生態(tài)用地空間連通的趨勢[20]。本研究的MCR模型是建立在城市和生態(tài)用地競合過程之上，為了建立可比性，需要建立統(tǒng)一的阻力面體系?；谙嚓P(guān)文獻(xiàn)[22-23]，結(jié)合鎮(zhèn)區(qū)現(xiàn)狀，選取了土壤重金屬綜合污染指數(shù)、土地類型、高程、坡度、歸一化植被指數(shù)、交通、河流、距城鎮(zhèn)距離、距礦區(qū)距離9個指標(biāo)為評價因子來構(gòu)建阻力表面評價體系(表1)，為使評價因子具有可比性，在衡量阻力因子時，通過查閱資料，將各項阻力因素用5分制進(jìn)行等級劃分(圖3)，最后各阻力因子按權(quán)重[20,24-25]疊加構(gòu)成城鎮(zhèn)和生態(tài)擴張的綜合因子阻力面(圖4)。

圖2 城市與生態(tài)擴張源地分布

一般而言高差大，植被指數(shù)高、坡度大的適合生態(tài)用地，反之適合城市建設(shè)用地。鑒于本研究的生態(tài)修復(fù)潛力與優(yōu)先級，是在城市空間和生態(tài)兩者用地擴張競合量化工礦區(qū)生態(tài)修復(fù)范圍的基礎(chǔ)上，因此需要采用構(gòu)建綜合生態(tài)用地MCR阻力表面減去綜合城市建設(shè)用地MCR阻力表面，得到兩過程的綜合ΔMCR差值表面。

3.3 最小累積阻力表面的計算

根據(jù)構(gòu)建阻力可否直接判斷，可將阻力分為顯性與隱形阻力兩類[26]。以傳統(tǒng)構(gòu)建方法利用MCR模型構(gòu)建出生態(tài)用地和城鎮(zhèn)兩者擴展源的顯性阻力表面，其次為減少單一阻力面導(dǎo)致阻力表面的構(gòu)建與實際情況偏差，考慮物質(zhì)與能量之間的相互影響，采用地統(tǒng)計學(xué)中的Kriging以空間自相關(guān)為基礎(chǔ)構(gòu)建出兩者的隱性阻力面，參考已有的研究成果[26]，確定顯性與隱性兩者阻力面權(quán)重值為0.7和0.3，構(gòu)建出兩者綜合MCR阻力表面(圖5)。最后采用式(4)計算出生態(tài)與城市用地擴張ΔMCR差值表面(圖6)，其呈現(xiàn)出大值區(qū)分布在北部的志輝冶煉廠區(qū)，而小值區(qū)則集中在南部的裸露礦區(qū)，并由城市中心區(qū)域向周邊減少。

3.4 生態(tài)修復(fù)用地潛力

根據(jù)最小累積阻力表面的差值結(jié)果，選取ΔMCR<0時宜為生態(tài)用地的值，用自然斷點法將鎮(zhèn)區(qū)劃分Ⅰ～Ⅴ級的潛力區(qū)(圖7)，其阻力相差值越大，則潛力區(qū)級別越大，修復(fù)為城市生態(tài)用地的后備潛力越高而生態(tài)修復(fù)優(yōu)先性越低。通過計算(表2)鎮(zhèn)區(qū)適宜地后備潛力面積為5 603.22 hm2，Ⅰ與Ⅱ級潛力區(qū)占到生態(tài)適宜地后備潛力總面積的28.65%和26.57%，Ⅲ與Ⅳ級潛力區(qū)占到23.97%和15.59%，Ⅴ級潛力區(qū)達(dá)到5.22%。

表1 阻力因子評價指標(biāo)體系

依據(jù)工礦區(qū)范圍內(nèi)不同區(qū)域優(yōu)化對生態(tài)系統(tǒng)功能發(fā)揮的作用有不同影響，區(qū)內(nèi)4種潛力區(qū)的修復(fù)優(yōu)先性分別為Ⅰ區(qū)>Ⅱ區(qū)>Ⅲ區(qū)>Ⅳ區(qū)(圖7)。通過統(tǒng)計(表2)工礦生態(tài)適宜地后備潛力面積為333.09 hm2，Ⅰ級與Ⅱ級潛力區(qū)面積為166.41 hm2和90.9 hm2，分別占到工礦區(qū)生態(tài)適宜地后備潛力總面積的49.96%與 27.29%；Ⅲ級與Ⅳ級潛力區(qū)面積為73.98 hm2和1.8 hm2，分別占到工礦區(qū)生態(tài)適宜地后備潛力總面積的22.21%與0.54%；其中金銅項目冶煉等集中的工業(yè)園區(qū)、康家灣礦區(qū)、水口山金屬有色公司和部分裸露的采礦區(qū)以及尾礦填埋區(qū)為重點潛力區(qū)。在今后規(guī)劃設(shè)計中，可以將這些區(qū)域作為重點生態(tài)項目建設(shè)區(qū)。

4 結(jié)論與討論

基于ΔMCR差值表面分析和修復(fù)優(yōu)先級分區(qū)，工礦區(qū)內(nèi)修復(fù)生態(tài)用地的面積占工礦用地面積的59.62%，可為城市生態(tài)用地建設(shè)提供333.09 hm2用地。就工礦城鎮(zhèn)而言，對廢棄工礦區(qū)生態(tài)修復(fù)不僅可以維護(hù)生態(tài)系統(tǒng)的穩(wěn)定性也一定程度上緩解了城市問題，將城市和生態(tài)兩者擴張因素納入生態(tài)修復(fù)的評價體系中，對生態(tài)修復(fù)依據(jù)更具有科學(xué)性和指導(dǎo)性價值。

根據(jù)不同修復(fù)優(yōu)先級分區(qū)的治理重點，可采用有區(qū)別的措施與規(guī)劃策略進(jìn)行生態(tài)修復(fù)。Ⅰ級潛力區(qū)是最主要的生態(tài)用地修復(fù)區(qū)，占到工礦區(qū)生態(tài)適宜地后備潛力總面積的49.96%，其主要分布在金銅項目冶煉等集中的工業(yè)園區(qū)，在治理方面該區(qū)域是污染物擴散的源頭，重金屬復(fù)合污染程度強，工業(yè)冶煉活動產(chǎn)生的“三廢”導(dǎo)致重金屬的沉降和集聚，降低土壤中的重金屬是生態(tài)修復(fù)的核心?？刹捎玫统杀厩矣行У木G色修復(fù)技術(shù)—微生物修復(fù)或植物修復(fù)，通過微生物的吸附、轉(zhuǎn)化來降低土壤重金屬，或通過超積累植物利用其對污染元素特殊的富集能力，將重金屬從土壤中移除。如狼把草(Bidenstripartita)、龍葵(Solanumnigrum)對Cd和Zn有富集作用，蜈蚣草(Eremochloaciliaris)可修復(fù)萃取Pb-As、Zn-As、Cu-As等復(fù)合污染土壤。在規(guī)劃策略上，由于距城鎮(zhèn)中心區(qū)較近，可納入城市用地的開敞空間來協(xié)調(diào)發(fā)展，在城鎮(zhèn)轉(zhuǎn)型發(fā)展期，依托現(xiàn)存工廠遺址，在保留工業(yè)文化記憶基礎(chǔ)上發(fā)展第3產(chǎn)業(yè)。

Ⅱ級潛力區(qū)占到工礦區(qū)生態(tài)適宜地后備潛力總面積的27.29%，其主要分布在尾礦填埋區(qū)、水口山金屬有限公司及部分裸露礦區(qū)。長期的選礦及尾礦填埋是區(qū)域的主要活動，選礦廢水的排放以及坑塘洼地內(nèi)大面積廢棄尾礦的堆積，對水資源造成了嚴(yán)重污染。

圖3 單因素阻力面

圖4 生態(tài)擴張與城鎮(zhèn)擴張綜合阻力面

水質(zhì)的治理是該區(qū)域生態(tài)修復(fù)的重點，選礦區(qū)應(yīng)主要以治污水、抓節(jié)水為目標(biāo)?？刹扇∥鬯至鞔胧┣壹薪ㄔO(shè)污水處理系統(tǒng)，達(dá)到污水的排放標(biāo)準(zhǔn)后循環(huán)使用，從根本上治理了超標(biāo)污水的排放且降低了治水成本和保護(hù)了水資源，同時加大水系淤泥清理和河堤工程建設(shè)力度。在規(guī)劃策略上，可將部分修復(fù)成小型的綠色場地，如經(jīng)濟型農(nóng)林用地等。

Ⅲ級和Ⅳ級潛力區(qū)面積分別占到工礦區(qū)生態(tài)適宜地后備潛力總面積的22.21%與0.54%，由于長期的采礦活動，區(qū)域植被覆蓋度低，土地結(jié)構(gòu)損壞嚴(yán)重，在治理方面，植被恢復(fù)是區(qū)域生態(tài)修復(fù)的側(cè)重點，針對自然環(huán)境及露天礦區(qū)構(gòu)成特點采用綜合治理的方式進(jìn)行植被恢復(fù)，在極端的環(huán)境下，應(yīng)選擇耐干旱瘠薄、抗逆性好、優(yōu)良土壤能力、根系發(fā)達(dá)、成活率強的植物種類，如紫穗槐(Amorphafruticosa)、小葉榕(F.microcarpa)、車桑子(Dodonaeaviscosa)、類蘆(Neyraudiareynaudiana)、滇合歡(Albiziasimeonis)等，在群落配置中，以草灌植物為主，突出以綠色為基調(diào)的層次感。為彰顯地方特色，適當(dāng)選用一定量的當(dāng)?shù)叵蠕h喬木，如樟樹(Cinnamomumcamphora)等。喬木、灌木和草本的配置面積比例分別是10%～20%、20%～30%、50%～70%。為提高存活率，根據(jù)不同的地類和植被類型進(jìn)行整地，采取覆土措施、自然與人為增肥措施恢復(fù)土壤基質(zhì)和肥力。在規(guī)劃策略上，建議結(jié)合部分裸露的礦區(qū)修復(fù)為城市的郊野公園與半開放的生態(tài)保育區(qū)。

圖5 生態(tài)與城市兩者擴張綜合最小累積阻力表面

圖6 基于ΔMCR模型結(jié)果

本研究探索了在工礦區(qū)小尺度上生態(tài)修復(fù)潛力的評價方法，其阻力因子體系及模型的構(gòu)建可為類似的研究提供借鑒。但仍有幾個方面需要強調(diào)：1)最小阻力模型建立是在指標(biāo)體系構(gòu)建的基礎(chǔ)之上，模型中阻力因子參數(shù)是相對值而非絕對值且評價體系無量綱，其權(quán)重的設(shè)定都是基于相關(guān)研究成果進(jìn)行的，有待深入探討；2)模型構(gòu)建方面，只考慮單元有促進(jìn)或阻礙作用以及社會效益最大化和生態(tài)受到最小破壞的用地競合，而未考慮政策制度、技術(shù)水平、開發(fā)理念等一些非指標(biāo)因素的影響；3)在內(nèi)梅羅綜合污染指數(shù)方面，只在0～20 cm土壤取樣測算，可進(jìn)一步探討土壤縱向深度各金屬的變化，對重金屬綜合污染分布研究進(jìn)行完善。今后可在完善阻力因子及模型體系構(gòu)建基礎(chǔ)上進(jìn)一步對工礦區(qū)的修復(fù)利用做深入的探究，如建立完整的工礦城鎮(zhèn)修復(fù)體系、建立工礦城鎮(zhèn)多元化修復(fù)模式以及多產(chǎn)業(yè)重組模式等。但本研究結(jié)果可為水口山鎮(zhèn)規(guī)劃設(shè)計以及修復(fù)項目在空間尺度上落地和科學(xué)治理提供客觀依據(jù)，以期促進(jìn)工礦城鎮(zhèn)的可持續(xù)發(fā)展，并為同類型礦區(qū)的生態(tài)修復(fù)提供借鑒。

注:1.志輝冶煉；2.松柏化肥廠；3.宏緯度冶煉廠；4.化工四廠；5.金銅項目冶煉集中工業(yè)園區(qū)；6.康家灣礦區(qū)；7.水口山金屬有限公司；8.裸露礦區(qū)；9.尾礦填埋區(qū)；10.康華化肥公司。

Fig.7 Spatial distribution of ecological restoration potential grades in towns and industrial and mining areas