煤礦區(qū)土地利用綜合效益耦合關(guān)系分析

崔藝馨，李永峰，王曉輝

(中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院，江蘇 徐州 221116)

土地是人類生存和發(fā)展的必要和基礎(chǔ)資源，在人地系統(tǒng)得長(zhǎng)期演化過(guò)程中，土地和人類形成了互為反饋的機(jī)制。[1]一方面土地資源利用影響著人類社會(huì)經(jīng)濟(jì)的發(fā)展，而人類社會(huì)經(jīng)濟(jì)發(fā)展反過(guò)來(lái)也會(huì)作用于土地，尤其是工業(yè)化和城鎮(zhèn)化，土地的生態(tài)條件與環(huán)境質(zhì)量在此過(guò)程中受到影響。顯然，在煤礦區(qū)土地?fù)p毀更加嚴(yán)重，生態(tài)環(huán)境問(wèn)題突出，土地資源利用引發(fā)的社會(huì)經(jīng)濟(jì)發(fā)展與生態(tài)間矛盾愈加尖銳，因此研究土地利用效益可以為制定可持續(xù)的土地利用規(guī)劃服務(wù)，通過(guò)解讀土地利用效益耦合關(guān)系趨勢(shì)可以提高礦區(qū)土地復(fù)墾與生態(tài)修復(fù)的及時(shí)性[2]。

1 煤礦區(qū)土地利用綜合評(píng)價(jià)指標(biāo)體系

1.1 煤礦區(qū)土地利用社會(huì)效益與生態(tài)效益耦合機(jī)理

在不同礦井生命周期階段，由于生產(chǎn)組織的重點(diǎn)不同，煤炭開(kāi)采對(duì)礦區(qū)生態(tài)環(huán)境的影響程度和方式也會(huì)有所不同，又由于煤炭和土地具有同位性，土地資源利用方式與狀態(tài)也會(huì)受到不同影響，遵循礦井的生命周期。[3]土地利用效益是土地利用是否可持續(xù)的檢驗(yàn)標(biāo)準(zhǔn)，受到礦區(qū)特點(diǎn)的影響，土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益的影響因素相互滲透、相互促進(jìn)又相互制約，遵循著礦井生命周期的演變規(guī)律(圖1)。

圖1 煤礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益耦合機(jī)理

在規(guī)劃建設(shè)期，煤炭企業(yè)主要進(jìn)行生產(chǎn)和生活區(qū)建設(shè)，對(duì)地表基本不產(chǎn)生影響，此時(shí)企業(yè)為礦區(qū)居民提供了更多就業(yè)崗位，此時(shí)主要表現(xiàn)為對(duì)土地利用社會(huì)經(jīng)濟(jì)的正效益；投產(chǎn)達(dá)產(chǎn)期，礦井開(kāi)始投入生產(chǎn)，且強(qiáng)度和規(guī)模不斷擴(kuò)大，煤炭產(chǎn)量增加，“三廢”排放增加，開(kāi)始出現(xiàn)土地塌陷與占用等，影響土地利用的生態(tài)效益，但是從整體上，煤炭資源開(kāi)發(fā)引導(dǎo)礦區(qū)經(jīng)濟(jì)逐步繁榮，居民生活水平提高，福利增加，提高了土地利用的社會(huì)經(jīng)濟(jì)效益；到達(dá)穩(wěn)產(chǎn)期后，煤炭資源開(kāi)發(fā)對(duì)土地生態(tài)的影響也逐漸趨于最大化，土地?fù)p毀加劇后果就是土地的完整性遭到了破壞，而耕地的耕作條件也被改變，生產(chǎn)力開(kāi)始下降[4]，道路也受到破壞，可能出現(xiàn)廢棄不可利用的荒地。居民生活受到資源開(kāi)采的負(fù)面影響加劇，企業(yè)與地方矛盾突出，此時(shí)煤炭資源開(kāi)發(fā)對(duì)土地利用社會(huì)經(jīng)濟(jì)效益影響的正值已經(jīng)不足以抵消其對(duì)土地利用生態(tài)效益影響的負(fù)值，生態(tài)環(huán)境開(kāi)始制約社會(huì)經(jīng)濟(jì)的發(fā)展；到礦井衰退期，土地利用的社會(huì)經(jīng)濟(jì)效益下降，土地復(fù)墾與生態(tài)修復(fù)的任務(wù)突出，及時(shí)的具有針對(duì)性的措施會(huì)阻止土地生態(tài)環(huán)境的破壞繼續(xù)擴(kuò)大，同時(shí)非煤替代產(chǎn)業(yè)的發(fā)展，都可以緩解生態(tài)環(huán)境破壞對(duì)社會(huì)經(jīng)濟(jì)效益影響的抑制作用，可以加減緩煤炭資源開(kāi)采對(duì)土地利用生態(tài)效益負(fù)值的影響。

1.2 煤礦區(qū)土地利用綜合效益評(píng)價(jià)指標(biāo)體系

基于土地利用社會(huì)效益與生態(tài)效益耦合機(jī)理分析，考慮影響土地利用效益的影響因素，遵循科學(xué)性、合理性、可操作性等原則，構(gòu)建煤礦區(qū)土地利用效益評(píng)價(jià)指標(biāo)體系(表1)。

1.3 評(píng)價(jià)方法與評(píng)價(jià)標(biāo)準(zhǔn)

1.3.1 評(píng)價(jià)模型與方法

1.3.1.1 數(shù)據(jù)標(biāo)準(zhǔn)化處理

由于選取的指標(biāo)來(lái)源不同，單位和數(shù)量級(jí)也不相同，為方便指標(biāo)間比較，對(duì)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理。

(1)

(2)

1.3.1.2 確定評(píng)價(jià)指標(biāo)權(quán)重

熵權(quán)法、層次分析法、德?tīng)柗品ê妥儺愊禂?shù)法是常用的確定權(quán)重的方法，相對(duì)也較為成熟。層次分析法和德?tīng)柗品ㄊ侵饔^賦權(quán)的方法，具有一定的主觀性，煤礦區(qū)土地利用效益評(píng)價(jià)相關(guān)研究極少，沒(méi)有參考，主觀性過(guò)大會(huì)導(dǎo)致誤差增大；而熵權(quán)法和變異系數(shù)法都屬于客觀賦權(quán)的方法，其中變異系數(shù)法更適宜于評(píng)價(jià)指標(biāo)對(duì)于評(píng)價(jià)目標(biāo)比較模糊的情況，且此種方法對(duì)指標(biāo)的經(jīng)濟(jì)意義重視不夠，顯然不適于用于效益評(píng)價(jià)。因此，最終選擇熵權(quán)法確定指標(biāo)權(quán)重。

各指標(biāo)熵值，表示為式(3)、式(4)。定義權(quán)重值，如式(5)所示。

ei=-∑ρimln(ρim),k=1/ln(n)

(3)

(4)

式中，n代表年份，ρit表示第i項(xiàng)指標(biāo)下第t個(gè)方案指標(biāo)值的比重。

(5)

式中，m代表m個(gè)指標(biāo)。

1.3.1.3 耦合模型

通過(guò)式(1)～(5)可以計(jì)算出各個(gè)指標(biāo)的權(quán)重，在此基礎(chǔ)上通過(guò)式(6)，可以計(jì)算出煤礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益I1和生態(tài)效益I2。

(6)

據(jù)分析，煤礦區(qū)土地利用的社會(huì)經(jīng)濟(jì)系統(tǒng)和生態(tài)環(huán)境系統(tǒng)具有明顯的耦合關(guān)系，根據(jù)系統(tǒng)論，可以將整個(gè)煤礦區(qū)的土地系統(tǒng)看作是由社會(huì)經(jīng)濟(jì)系統(tǒng)和生態(tài)系統(tǒng)兩個(gè)系統(tǒng)復(fù)合而成的系統(tǒng)，土地利用效益受到社會(huì)經(jīng)濟(jì)效益和生態(tài)效益兩個(gè)主導(dǎo)因素影響。通過(guò)曲線擬合可以得到兩種效益的演化方程，見(jiàn)式(7)。

表1 煤礦區(qū)土地利用效益評(píng)價(jià)指標(biāo)體系

(7)

式中，A表示土地利用社會(huì)經(jīng)濟(jì)效益子系統(tǒng)的演化方程，B表示生態(tài)效益子系統(tǒng)的演化方程。A與B兩個(gè)子系統(tǒng)相互影響，而A、B與I1、I2相互作用，整個(gè)系統(tǒng)中的所有要素相互關(guān)聯(lián)，相互滲透，共同形成一個(gè)整體。由此通過(guò)對(duì)演化方程(式(8))求導(dǎo)，得到土地利用社會(huì)經(jīng)濟(jì)效益演化速度VA和生態(tài)效益演化速度VB，整個(gè)系統(tǒng)的演化速度V與子系統(tǒng)演化速度VA和VB相關(guān)，由此通過(guò)判斷兩個(gè)子系統(tǒng)直接的耦合度來(lái)判斷整個(gè)系統(tǒng)的協(xié)調(diào)發(fā)展?fàn)顩r。

(8)

在理論上，社會(huì)經(jīng)濟(jì)的發(fā)展符合S型演變理論，雖然煤礦區(qū)具有特殊性，但是結(jié)合礦井生命周期，可以判斷出土地利用的社會(huì)經(jīng)濟(jì)效益符合S型演變理論，而土地利用的生態(tài)效益與社會(huì)經(jīng)濟(jì)相互耦合，因此也應(yīng)該符合S型理論，但是由于煤炭資源開(kāi)采對(duì)生態(tài)環(huán)境影響的累積效應(yīng)與滯后效應(yīng)，導(dǎo)致土地利用生態(tài)效益的演化速度慢于社會(huì)經(jīng)濟(jì)效益。因此，在以VA作為縱軸，VB作為橫軸，以V的變化軌跡構(gòu)成非線性模型時(shí)，形成的圖形應(yīng)該是橢圓型(圖2)。V與橫軸形成一個(gè)夾角α，即為耦合度，可以表示為式(9)。

(9)

圖2 土地利用效益演化速度

1.3.2 評(píng)價(jià)標(biāo)準(zhǔn)

通過(guò)求取α，可以確定煤礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益的耦合程度，以判斷土地利用的協(xié)調(diào)程度。依據(jù)圖2，將耦合度分為低級(jí)共生、協(xié)調(diào)發(fā)展、極限發(fā)展和再生發(fā)展四種類型(表2)。

表2 土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益耦合度判斷標(biāo)準(zhǔn)

表3 煤礦區(qū)土地利用效益指標(biāo)標(biāo)準(zhǔn)化值及其權(quán)重

2 朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益耦合關(guān)系

2.1 研究區(qū)概況

朔南礦區(qū)地處山西朔州市朔城區(qū)，東經(jīng)112°15′～112°34′，北緯39°07′39°25′,礦區(qū)面積580 km2[5],屬于干旱半干旱地區(qū)，生態(tài)脆弱，四季分明。整個(gè)礦區(qū)包括馬營(yíng)堡、豐予、高莊、梵王寺和麻家梁5個(gè)礦井和1個(gè)張家嘴勘查區(qū)。

2.2 指標(biāo)數(shù)據(jù)處理與權(quán)重確定

由于朔城區(qū)區(qū)域劃分以行政邊界為基礎(chǔ)，而礦區(qū)劃分以礦區(qū)的邊界為基礎(chǔ)，而可直接獲取的統(tǒng)計(jì)年鑒數(shù)據(jù)以行政邊界為依據(jù)。因此，本文朔南礦區(qū)社會(huì)經(jīng)濟(jì)數(shù)據(jù)的獲取是以朔州市朔城區(qū)2005～2014年統(tǒng)計(jì)年鑒為基礎(chǔ)，結(jié)合礦區(qū)實(shí)地實(shí)際調(diào)研調(diào)節(jié)而來(lái)。根據(jù)礦區(qū)開(kāi)采的空間布局，針對(duì)不同生命周期階段開(kāi)采后的地表狀況進(jìn)行實(shí)地調(diào)研，包括地表構(gòu)建筑物的毀損情況、地表植被的破壞情況、地表農(nóng)作物的生長(zhǎng)-收成情況，對(duì)照統(tǒng)計(jì)年鑒數(shù)據(jù)，反映礦區(qū)土地利用受采礦影響的程度；同時(shí)，考慮到朔南礦區(qū)雖然屬于朔城區(qū)范圍內(nèi)，但地處農(nóng)村，結(jié)合了統(tǒng)計(jì)年鑒朔城區(qū)農(nóng)民人均收入數(shù)據(jù)，調(diào)研礦區(qū)內(nèi)農(nóng)村居民的收入狀況，為礦山提供勞務(wù)、服務(wù)的情況，礦山對(duì)破壞土地的補(bǔ)償情況等，了解礦山開(kāi)采對(duì)礦區(qū)居民收入、生活的影響。對(duì)照統(tǒng)計(jì)年鑒，反映煤礦開(kāi)采不同階段對(duì)當(dāng)?shù)鼐用裆畹挠绊憽?/p>

本文朔南礦區(qū)煤炭產(chǎn)量主要是以統(tǒng)計(jì)年鑒為基礎(chǔ)，結(jié)合實(shí)地調(diào)研對(duì)缺失數(shù)據(jù)加以補(bǔ)充。土地利用數(shù)據(jù)獲取主要通過(guò)朔城區(qū)與朔南礦區(qū)土地利用現(xiàn)狀。另外，根據(jù)何書(shū)金，蘇光全測(cè)算，井工開(kāi)采每萬(wàn)噸煤造成的土地塌陷面積為0.27 hm2，土地壓占面積為0.005 hm2。[6]據(jù)此推算出朔南礦區(qū)土地塌陷面積與土地壓占面積。在獲取到基礎(chǔ)數(shù)據(jù)的基礎(chǔ)上，通過(guò)相關(guān)計(jì)算計(jì)算指標(biāo)值，對(duì)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，利用熵權(quán)法取得權(quán)重(表3)。

2.3 耦合結(jié)果分析

2.3.1 煤礦區(qū)土地利用綜合效益計(jì)算

通過(guò)式(6)計(jì)算出朔南礦區(qū)2005～2014年土地利用社會(huì)經(jīng)濟(jì)效益和生態(tài)效益的綜合值，演化曲線見(jiàn)圖3。從圖3中可以看出，煤礦區(qū)的社會(huì)經(jīng)濟(jì)效益與生態(tài)效益呈現(xiàn)不同的演變趨勢(shì)，矛盾突出。2009年以前，由于整個(gè)煤礦區(qū)還處于投產(chǎn)初期，煤炭產(chǎn)量多以小煤窯生產(chǎn)為主，產(chǎn)量較少[7]，對(duì)土地的破壞也基本沒(méi)有顯現(xiàn)，此時(shí)土地利用生態(tài)效益波動(dòng)不明顯，并且高于土地利用的社會(huì)經(jīng)濟(jì)效益；隨著整個(gè)礦區(qū)較大型礦井的投產(chǎn)，與生產(chǎn)規(guī)模的擴(kuò)大，土地?fù)p壞開(kāi)始顯現(xiàn)并日益嚴(yán)重，而礦區(qū)以農(nóng)業(yè)為主的產(chǎn)業(yè)結(jié)構(gòu)打破，社會(huì)經(jīng)濟(jì)受煤炭產(chǎn)量影響迅速發(fā)展，因此在2009年以后煤礦區(qū)土地利用的社會(huì)經(jīng)濟(jì)效益超過(guò)其生態(tài)效益，且差距日益拉大，也直觀顯示出煤礦區(qū)土地利用與生態(tài)環(huán)境的矛盾。

2.3.2 土地利用效益耦合度計(jì)算

借助EXCEL對(duì)圖3的社會(huì)經(jīng)濟(jì)效益與生態(tài)效益兩條曲線進(jìn)行擬合，結(jié)合實(shí)際情況，通過(guò)對(duì)比擬合度的大小，選擇適合的趨勢(shì)曲線類型。社會(huì)經(jīng)濟(jì)效益選擇乘冪趨勢(shì)線，擬合度為0.9364，而生態(tài)效益采用多項(xiàng)擬合，擬合度為0.9404。由此得出關(guān)于時(shí)間t的朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益A和生態(tài)效益B的函數(shù)關(guān)系式，見(jiàn)式(10)和式(11)。

A=0.0921t0.7058

(10)

B=0.0014t3-0.0214t2+0.054t+0.3323

(11)

式中，t=1，2，3，…，10，分別對(duì)應(yīng)2005～2014年的10年時(shí)間。

按照式(8)對(duì)土地利用社會(huì)經(jīng)濟(jì)效益A和生態(tài)效益B的函數(shù)方程進(jìn)行求導(dǎo)，得到朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益演化速度VA與生態(tài)效益的演化速度VB的函數(shù)關(guān)系式，見(jiàn)式(12)、式(13)。

VA=0.065t-0.2942

(12)

VB=0.0042t3-0.0656t2+0.0978

(13)

在計(jì)算出VA和VB的基礎(chǔ)上，通過(guò)式(9)計(jì)算出朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益的耦合度α，并繪制耦合度曲線圖，見(jiàn)表4和圖4。

圖3 朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益綜合值

年份VAVBVA/VBα(°)200500650015442211766720060053-00148-358191056200700471-00366-1285612788200800432-005-0864713915200900405-0055-0736114364201000384-00516-0743614336201100367-00398-0921413734201200353-00196-1798911907201300341000937841752201400330046071783567

圖4 朔南礦區(qū)土地利用社會(huì)經(jīng)濟(jì)效益與生態(tài)效益耦合度

根據(jù)圖4可以看出，2005～2014年，朔南礦區(qū)土地利用綜合效益的耦合度呈現(xiàn)倒“U型”，從76.67°到143.64°，又降回到35.67°，經(jīng)歷了協(xié)調(diào)發(fā)展和極限發(fā)展兩個(gè)階段，并且大部分年份處于極限發(fā)展階段，即由于煤炭資源開(kāi)發(fā)對(duì)土地生態(tài)環(huán)境的破壞，土地利用的生態(tài)效益開(kāi)始制約其社會(huì)經(jīng)濟(jì)效益，兩者矛盾日益突出，但后期受到生態(tài)環(huán)境保護(hù)措施的影響，矛盾又有所緩解。2005～2014年朔南礦區(qū)土地利用綜合效益的耦合度大致可以分為三個(gè)階段。

第一階段：協(xié)調(diào)發(fā)展階段，2005年，VA>VB。此時(shí)，煤礦區(qū)土地利用還未有明顯破壞，生態(tài)效益較大，但是受到煤炭產(chǎn)業(yè)的影響，大力拉動(dòng)了社會(huì)經(jīng)濟(jì)發(fā)展，因此社會(huì)經(jīng)濟(jì)效益的增長(zhǎng)速度迅速。說(shuō)明在礦區(qū)投產(chǎn)達(dá)產(chǎn)初期，煤炭產(chǎn)業(yè)帶動(dòng)區(qū)域土地社會(huì)經(jīng)濟(jì)效益快速增長(zhǎng)的同時(shí)，仍然處于土地的生態(tài)承載力范圍之內(nèi)。同時(shí)，在理論上煤炭資源開(kāi)發(fā)對(duì)土地生態(tài)環(huán)境的破壞在加劇，土地利用的社會(huì)經(jīng)濟(jì)效益與生態(tài)效益矛盾也在加劇，但煤炭資源開(kāi)發(fā)帶來(lái)的社會(huì)經(jīng)濟(jì)效益在某種程度上也增大了礦區(qū)的生態(tài)環(huán)境容量，因此生態(tài)效益也在上升[8]。

第二階段：極限發(fā)展階段，2006～2007年，VA>VB。此時(shí)土地利用的社會(huì)經(jīng)濟(jì)效益仍然處于上升狀態(tài)，但煤炭資源開(kāi)發(fā)對(duì)土地生態(tài)環(huán)境的負(fù)效應(yīng)開(kāi)始顯現(xiàn)，開(kāi)始限制土地利用社會(huì)經(jīng)濟(jì)效益的發(fā)展速度。2008～2011年，VAVB。社會(huì)經(jīng)濟(jì)效益的上升速度大于生態(tài)效益的下降速度，但是上升速度和生態(tài)效益下降速度都放緩。一方面，說(shuō)明受到資源有限性與開(kāi)采規(guī)模的限制，加之土地的塌陷與壓占造成的負(fù)面影響出現(xiàn)，土地利用的社會(huì)經(jīng)濟(jì)效益也開(kāi)始放緩；另一方面，說(shuō)明朔城區(qū)委、區(qū)政府“創(chuàng)優(yōu)社會(huì)，發(fā)展環(huán)境，改變作風(fēng)”的三項(xiàng)舉措在礦區(qū)生態(tài)環(huán)境保護(hù)中取得了一定成效。

第三階段：協(xié)調(diào)發(fā)展階段，2013年，VA>VB。2014年，VA

3 結(jié) 論

判斷煤礦區(qū)土地利用綜合效益的耦合度可以將土地利用效益系統(tǒng)分為社會(huì)經(jīng)濟(jì)子系統(tǒng)和生態(tài)子系統(tǒng)，兩個(gè)子系統(tǒng)關(guān)于時(shí)間呈非線性變化，通過(guò)測(cè)算兩者之間的耦合度可以判斷兩種效益的協(xié)調(diào)情況以及煤礦區(qū)土地持續(xù)利用情況[10]。本文以山西省朔南礦區(qū)為例，2005～2014年，朔南礦區(qū)土地利用綜合效益的耦合度經(jīng)歷了協(xié)調(diào)發(fā)展到極限發(fā)展兩個(gè)階段，并發(fā)現(xiàn)土地利用的社會(huì)經(jīng)濟(jì)效益與生態(tài)效益基本符合礦井生命周期規(guī)律，土地利用的社會(huì)經(jīng)濟(jì)效益不斷上升，但隨著煤炭資源開(kāi)采對(duì)土地的破壞，生態(tài)效益下滑，其對(duì)社會(huì)經(jīng)濟(jì)效益的限制愈加明顯，土地利用社會(huì)經(jīng)濟(jì)效益增長(zhǎng)速度變緩。朔南礦區(qū)在干旱半干旱地區(qū)的煤礦區(qū)中具有一定的代表性，說(shuō)明煤礦區(qū)土地利用的社會(huì)經(jīng)濟(jì)效益與生態(tài)效益矛盾較為突出，與非資源型城市土地利用的效益耦合度具有明顯差異[11-12]。因此，應(yīng)結(jié)合礦井生命周期，把握土地利用效益的變化規(guī)律，在促進(jìn)社會(huì)經(jīng)濟(jì)發(fā)展的同時(shí)，及時(shí)發(fā)展非煤產(chǎn)業(yè)或者進(jìn)行產(chǎn)業(yè)升級(jí)，也要合理配置土地資源，把握時(shí)機(jī)進(jìn)行土地整治與生態(tài)修復(fù)。

[1] 談兵.基于社會(huì)經(jīng)濟(jì)與生態(tài)環(huán)境耦合的城市土地利用效益評(píng)價(jià)[D].武漢：武漢科技大學(xué)，2013.

[2] Dewan A M，Yamaguchi Y.Land use and land cover change in Greater Dhaka，Bangladesh：Using remote sensing to promote sustainable urbanization[J].Applied Geography，2009，29(3)：390-401.

[3] 王行風(fēng)，汪云甲，李永峰.基于生命周期理論的煤礦區(qū)土地利用演化模擬[J].地理研究，2009，28(2)：379-390.

[4] 楊燕杰.基于RS/GIS的煤礦區(qū)土地利用變化及其生態(tài)效應(yīng)研究[D].濟(jì)南：山東師范大學(xué)，2008.

[5] 史慧君.山西省朔南礦區(qū)礦山地質(zhì)環(huán)境綜合評(píng)價(jià)[D].烏魯木齊：新疆大學(xué)，2015.

[6] 何書(shū)金，蘇光全.中國(guó)采礦業(yè)的發(fā)展與礦區(qū)土地?fù)p毀預(yù)測(cè)[J].資源科學(xué)，2002，24(2)：17-21.

[7] 宋成舜，崔薛華，柯新利，等.城市土地集約利用綜合效益演化研究——以湖北省咸寧市為例[J].水土保持研究，2013，20(1)：230-234.

[8] 熊征，談兵，宋成舜，等.城市土地利用綜合效益耦合分析——以武漢市為例[J].水土保持研究，2015，22(2)：278-283.

[9] 吳迪.煤礦區(qū)國(guó)土資源利用規(guī)劃協(xié)調(diào)度評(píng)價(jià)研究[D].徐州：中國(guó)礦業(yè)大學(xué)，2014.

[10] 曹堪宏，朱宏偉.基于耦合關(guān)系的土地利用效益評(píng)價(jià)——以廣州和深圳為例[J].中國(guó)農(nóng)村經(jīng)濟(jì)，2010(8)：58-66.

[11] 馬安勝，李誠(chéng)固，劉雙良.城市產(chǎn)業(yè)結(jié)構(gòu)與土地利用耦合關(guān)系研究——以沈陽(yáng)市為例[J].山西財(cái)經(jīng)大學(xué)學(xué)報(bào)，2012，34(4)：29-30.

[12] 趙曉露，高敏華，高軍.烏魯木齊市土地利用效益的耦合關(guān)系分析[J].干旱區(qū)資源與環(huán)境，2011，25(1)：91-95.