貴港市土地利用轉(zhuǎn)型及生態(tài)服務(wù)價(jià)值研究

姚富明 韋燕飛 童新華

摘要：隨著我國(guó)城鎮(zhèn)化的快速發(fā)展，目前已經(jīng)趨向于一個(gè)穩(wěn)定的狀態(tài)。怎樣合理使用土地并發(fā)揮生態(tài)效益最大化，是引人深思的問(wèn)題，土地利用轉(zhuǎn)型可以作為重要的突破口。以地處平原丘陵區(qū)的廣西壯族自治區(qū)貴港市為研究區(qū)域，在新型城鎮(zhèn)化背景下，以2010、2014、2018年3期貴港市土地利用現(xiàn)狀數(shù)據(jù)為基礎(chǔ)，運(yùn)用Excel軟件，從土地利用主導(dǎo)功能角度研究2010—2018年貴港市土地利用轉(zhuǎn)型和生態(tài)價(jià)值，結(jié)果表明，其土地利用轉(zhuǎn)型明顯，生態(tài)服務(wù)總價(jià)值逐漸降低，并得出以下結(jié)論：（1）貴港市農(nóng)業(yè)生產(chǎn)用地、林地生態(tài)用地大幅度減少，而城鎮(zhèn)生活用地、工礦生產(chǎn)用地大幅增加，但貴港市仍以農(nóng)業(yè)生產(chǎn)用地、林地生態(tài)用地為主。（2）貴港市的其他生態(tài)用地占比不少，表明貴港市土地資源開(kāi)發(fā)利用程度中等偏上，土地開(kāi)發(fā)后備資源充足。（3）在研究時(shí)間段2010—2018年期間，土地利用綜合動(dòng)態(tài)度提高較大。說(shuō)明貴港市土地利用隨著時(shí)間的延后，其改變速率也會(huì)跟著提升，土地利用轉(zhuǎn)型因人類活動(dòng)干擾而顯著上升。（4）在研究期間，生態(tài)服務(wù)總價(jià)值逐漸降低，貢獻(xiàn)者主要是林地，其次為水體、農(nóng)田。

關(guān)鍵詞：新型城鎮(zhèn)化;土地利用轉(zhuǎn)型;生態(tài)服務(wù)價(jià)值;貴港市;新型城鎮(zhèn)化

中圖分類號(hào)：F301.1;F062.1? 文獻(xiàn)標(biāo)志碼： A

文章編號(hào)：1002-1302（2021）15-0207-08

收稿日期：2020-11-21

基金項(xiàng)目：國(guó)家自然科學(xué)基金（編號(hào)：41861039）。

作者簡(jiǎn)介：姚富明（1998—），男，廣西陸川人，碩士研究生，主要從事土地利用優(yōu)化配置研究。E-mail：yaofmlm123@163.com。

通信作者：韋燕飛，博士，教授，主要從事土地規(guī)劃與評(píng)價(jià)、土地利用優(yōu)化配置研究。E-mail：gxweiyf121@163.com。

在1995年，IHDP（全球變化人文領(lǐng)域計(jì)劃）和IGBP（國(guó)際地圈——生物圈計(jì)劃）共同提出了一個(gè)研究計(jì)劃——LUCC（土地利用和土地覆蓋變化）。不久，這個(gè)計(jì)劃成為地理學(xué)的一個(gè)熱點(diǎn)。近年來(lái)，研究土地利用和土地覆被變化的重要方式就是研究土地利用轉(zhuǎn)型。土地利用轉(zhuǎn)型的表現(xiàn)形式之一就是土地利用變化，即土地利用方式（某個(gè)區(qū)域在特定時(shí)段的主要土地利用類型結(jié)構(gòu)）在時(shí)序上的改變[1]。土地利用轉(zhuǎn)型概念是Grainger第一次提出的;他認(rèn)為，一個(gè)國(guó)家在特定時(shí)期內(nèi)實(shí)際土地覆蓋總體布局就是國(guó)家的土地利用形式[2]。Vitousek等認(rèn)為，影響自然環(huán)境保護(hù)的最重要因素或許是土地利用變化[3]。Izquierdo等也贊成這一看法，并表明建立保護(hù)區(qū)是一項(xiàng)直接的土地利用策略，能避免土地利用變化對(duì)生物多樣性、資源和環(huán)境的潛在不好的侵染[4]。Lambin等是這樣理解土地利用轉(zhuǎn)型概念的，即土地利用不是線性的變化，是伴隨其他社會(huì)與生物系統(tǒng)變化而變化的，也沒(méi)有固定的模式和確定性，其轉(zhuǎn)型來(lái)源主要分為2種：內(nèi)生社會(huì)生態(tài)力、外生社會(huì)經(jīng)濟(jì)因素[5]。國(guó)內(nèi)對(duì)土地利用轉(zhuǎn)型的理解是從土地利用和土地覆被變化研究項(xiàng)目開(kāi)始的。蔡運(yùn)龍認(rèn)為，因?yàn)橥恋乩煤屯恋馗脖蛔兓且粋€(gè)極其龐大而復(fù)雜的項(xiàng)目，為了避免基于片面的觀點(diǎn)，需要找到一種新的綜合方法，包括驅(qū)動(dòng)力、尺度、理論和方法的結(jié)合[6]。龍花樓提出了土地利用轉(zhuǎn)型并首次將這種新方法運(yùn)用于研究長(zhǎng)江沿線樣帶的土地利用轉(zhuǎn)型，指出區(qū)域土地利用轉(zhuǎn)型是指區(qū)域土地利用（一定時(shí)期內(nèi)由主要土地利用類型）在時(shí)間排序上的改變[7]。龍花樓補(bǔ)充了土地利用轉(zhuǎn)型的概念，即在社會(huì)經(jīng)濟(jì)變化和革新的驅(qū)動(dòng)下，某一區(qū)域在一段時(shí)期內(nèi)由一種土地利用形態(tài)（含顯性形態(tài)和隱形形態(tài)）轉(zhuǎn)變?yōu)榱硪环N土地利用形態(tài)的過(guò)程[8]。

從2002年黨的十六大報(bào)告到2017年黨的十九大報(bào)告，提出要建立大中小城市并以城市群為主的協(xié)調(diào)發(fā)展的城市布局，加速推動(dòng)城鎮(zhèn)化歷程的農(nóng)轉(zhuǎn)民與市民化，再到2019年3月全國(guó)人大報(bào)告指出，要提升新型城鎮(zhèn)化品質(zhì)，呈現(xiàn)以人為核心，使城市更適合居住。在新型城鎮(zhèn)化進(jìn)程中，怎樣合理使用土地并發(fā)揮生態(tài)效益最大化，是引人深思的問(wèn)題。新型城鎮(zhèn)化建設(shè)的最佳路徑便是對(duì)土地利用轉(zhuǎn)型的研究。

從國(guó)內(nèi)外研究進(jìn)展上來(lái)看，能很好地理解土地利用轉(zhuǎn)型的起源、實(shí)質(zhì)及理論內(nèi)涵，但對(duì)于應(yīng)用其理論在實(shí)際研究土地利用轉(zhuǎn)型的還甚少，特別是像處于平原丘陵區(qū)的貴港市土地利用轉(zhuǎn)型。本研究將在結(jié)合國(guó)內(nèi)外研究現(xiàn)狀的基礎(chǔ)上，進(jìn)一步對(duì)處于平原丘陵區(qū)的貴港市的土地利用轉(zhuǎn)型進(jìn)行研究，并探討在土地利用轉(zhuǎn)型下，其生態(tài)服務(wù)價(jià)值量的變化關(guān)系，從而揭露二者之間的關(guān)系，并得出相應(yīng)的結(jié)論，提出相對(duì)應(yīng)的建議，既合理利用土地又能發(fā)揮生態(tài)服務(wù)價(jià)值量最大化，從而讓土地更好地為人類發(fā)展服務(wù)。

1 材料與方法

1.1 研究區(qū)概況

1.1.1 區(qū)位特征 廣西壯族自治區(qū)貴港市坐落在廣西東南部，東接梧州市、南鄰玉林市與欽州市、西接南寧市、北鄰來(lái)賓市，擁有華南內(nèi)河第一大港口和廣西最大的沖擊平原潯郁平原，北回歸線穿越中部地理位置較為優(yōu)越。

1.1.2 自然環(huán)境 貴港市行政轄區(qū)10 606 km2，土地資源種類豐富且利用結(jié)構(gòu)復(fù)雜;平原、丘陵是其主要地貌類型;所屬亞熱帶季風(fēng)氣候，溫暖濕潤(rùn)，降水量充足。

1.1.3 社會(huì)經(jīng)濟(jì) 至2018年末，建成區(qū)面積 79 km2，人口561.29萬(wàn)。貴港市2018年全年實(shí)現(xiàn)國(guó)內(nèi)生產(chǎn)總值（GDP）1 169.88億元，人均GDP為26 635元，地區(qū)生產(chǎn)總值農(nóng)業(yè)、工業(yè)、服務(wù)業(yè)比例為 17.0 ∶42.2 ∶40.8，城鎮(zhèn)化率為50.09%。

1.1.4 土地利用現(xiàn)狀 至2018年末，全市土地總面積為1 060 234.27 hm2，其中，農(nóng)用地為 885 461.47 hm2，比例達(dá)83.52%;建設(shè)用地為 90 948.45 hm2，比例為8.58%;其他土地 83 824.35 hm2，比例為7.91%（圖1）。農(nóng)用地中的耕地319 059.82 hm2（圖2）。

1.2 數(shù)據(jù)來(lái)源與處理

1.2.1 數(shù)據(jù)來(lái)源 采用的屬性數(shù)據(jù)主要來(lái)源于《貴港年鑒》（2010—2018年）、《貴港市土地利用年度變更調(diào)查》（2010—2018）、2018年《貴港市國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)》。

1.2.2 數(shù)據(jù)處理 當(dāng)前，對(duì)于“三生”用地的判定還沒(méi)有一致的標(biāo)準(zhǔn)。其劃分的方法主要有2種：（1）與土地利用分類系統(tǒng)二級(jí)地類相結(jié)合的歸并分類法;（2）量化測(cè)算法。本研究采用的是第1種。基于“三生”（生產(chǎn)、生活、生態(tài)）用地理論，分成3種一級(jí)地類，9種二級(jí)地類。從土地利用主導(dǎo)功能角度，對(duì)貴港市二級(jí)地類進(jìn)行再分類（表1）。土地利用主導(dǎo)功能分類是根據(jù)人的主觀用地意圖作為某一類土地的土地利用分類。一種土地利用分類系統(tǒng)的二級(jí)地類可能囊括2種或2種以上的功能，比如同時(shí)具有生產(chǎn)糧食和調(diào)節(jié)氣候的耕地，具有生產(chǎn)和生態(tài)功能。但我們使用耕地是生產(chǎn)糧食，因此被劃分為生產(chǎn)用地。

1.3 研究方法

1.3.1 土地利用結(jié)構(gòu)變化分析 土地利用結(jié)構(gòu)是指某區(qū)域內(nèi)不同利用土地之間的占比關(guān)系，其囊括時(shí)空2個(gè)方面的變化[9]。貴港市土地利用結(jié)構(gòu)變化體現(xiàn)在2010、2014、2018年這3期各個(gè)土地利用類型面積變化情況，可用來(lái)研究該市土地利用變化的趨勢(shì)和規(guī)律。

1.3.2 單一土地利用動(dòng)態(tài)度 單一土地利用動(dòng)態(tài)度指數(shù)是指某區(qū)域在特定時(shí)間限制內(nèi)某種用地類型面積的變化狀況[10]。公式如下：

K=Ub-UaUa×1T×100%。（1）

式中：K為土地利用動(dòng)態(tài)度;Ua、Ub分別為初期和末期某用地類型面積;T為研究時(shí)長(zhǎng)，當(dāng)T設(shè)定為年時(shí)，K值就是某用地類型的年變化率。

1.3.3 土地利用綜合動(dòng)態(tài)度 土地利用綜合動(dòng)態(tài)度是指某區(qū)域在特定時(shí)段內(nèi)整體土地利用變化速率。公式如下[11]：

LC=∑ni=1ΔLUi-j2∑ni=1LUi×1T×100%。（2）

式中：LC為土地利用綜合動(dòng)態(tài)度;LUi為研究開(kāi)始時(shí)i類用地類型面積;ΔLUi-j為研究時(shí)間內(nèi)i類轉(zhuǎn)變?yōu)閖類的用地類型面積;T為研究時(shí)間，當(dāng)T假定為年時(shí)，LC值就是某種用地類型的年變化率。

1.3.4 土地利用程度變化分析 土地利用程度變化分析是依據(jù)土地利用程度的界限來(lái)的。土地利用的頂點(diǎn)是人類不能再對(duì)土地進(jìn)行開(kāi)發(fā)利用，而低點(diǎn)是人類利用土地的始點(diǎn)[11]。莊大方等在中國(guó)資源環(huán)境數(shù)據(jù)庫(kù)中對(duì)4類土地利用的理想狀態(tài)定義為4類土地利用級(jí)，并賦予它本身類別的值（表2）。

土地利用程度可以用土地利用程度綜合指數(shù)體現(xiàn)。公式如下[12-13]：

L=100×∑ni=1Ai×Ci，L∈[100，400]。（3）

式中：L為土地利用程度綜合指數(shù);Ai為第 i級(jí)土地利用程度分級(jí)指數(shù);Ci為第i級(jí)土地利用程度分級(jí)面積比例;n為土地利用類型數(shù)量，本研究將其分為4類，故n取 4。

土地利用程度變化量和變化率表示研究區(qū)土地利用變化走向，公式如下：

ΔLa-b=Lb-La=100×（∑ni=1Ai×Cib-∑ni=1Ai×Cia）;（4）

R=∑ni=1（Ai×Cib）-∑ni=1（Ai×Cia）∑ni=1（Ai×Cib）。（5）

式中：ΔLa-b為土地利用程度變化量;R為土地利用程度變化率;Lb、La為b時(shí)間和a時(shí)間土地利用程度綜合指數(shù);Ai為第i級(jí)土地利用程度分級(jí)指數(shù);Cib、Cia為b時(shí)間和a時(shí)間第i級(jí)土地利用程度面積比例。

1.3.5 土地利用的擴(kuò)展強(qiáng)度指數(shù) 土地利用的擴(kuò)展強(qiáng)度指數(shù)，研究了不同土地類型的擴(kuò)展或收縮特征。擴(kuò)展強(qiáng)度指數(shù)是指特定時(shí)期內(nèi)土地利用擴(kuò)張面積占本身土地總面積的比例，數(shù)字越大，擴(kuò)展越快[14]。其模型如下：

UⅡ=[（ULAi，t+n-ULAi，t）n]×100/TLAi。（6）

式中：UⅡ?yàn)閿U(kuò)展強(qiáng)度指數(shù);ULAi，t+n、ULAi，t為空間單元i在t+n以及t年時(shí)該種土地的利用面積;TLAi為空間單元i的總面積。

1.3.6 生態(tài)服務(wù)價(jià)值計(jì)算模型 生態(tài)服務(wù)是指生態(tài)與生態(tài)系統(tǒng)形成的保持人類賴以生存和發(fā)展的環(huán)境條件和自然效用，即人類從生態(tài)系統(tǒng)中直接或間接獲得的利益，與人類福利密切相關(guān)[15]。評(píng)價(jià)生態(tài)系統(tǒng)服務(wù)功能高低的重要指標(biāo)是生態(tài)服務(wù)價(jià)值，其價(jià)值的變動(dòng)體現(xiàn)人類活動(dòng)對(duì)生態(tài)系的影響。本研究結(jié)合生態(tài)系統(tǒng)服務(wù)價(jià)值的相關(guān)研究，并參考國(guó)內(nèi)學(xué)者的用地分類方法，基于貴港市土地利用的實(shí)際和《土地利用現(xiàn)狀分類》（GB/T 21010—2007）做了歸并（表3）。

Costanza是世界上最早對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行估算的，但Costanza的估算適合西方等發(fā)達(dá)國(guó)家[16]。謝高地等從中國(guó)的實(shí)際情況出發(fā)，并參照Costanza的估算方法，對(duì)中國(guó)200名生態(tài)學(xué)家的問(wèn)卷，制定了中國(guó)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表[17]。本研究在這個(gè)當(dāng)量表的基礎(chǔ)上，根據(jù)貴港市的實(shí)際情況，對(duì)貴港市的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)進(jìn)行了調(diào)整。對(duì)貴港市生態(tài)價(jià)值有不好影響的是建設(shè)用地，對(duì)此其價(jià)值系數(shù)參照段瑞娟等[18]的，水域參照張鳳太等[19]的價(jià)值系數(shù);草地價(jià)值系數(shù)采用張騫[20]的。最后得到結(jié)果見(jiàn)表4。

使用公式（7）對(duì)貴港市的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行估算：

ESV=∑nk=1Ak×VCk。（7）

式中：ESV為生態(tài)系統(tǒng)服務(wù)價(jià)值;Ak為第k類用地類型面積;VCk為第 k 類用地類型生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù);n為土地利用類型數(shù)目。

1.3.7 生態(tài)價(jià)值系數(shù)對(duì)生態(tài)服務(wù)價(jià)值敏感度分析 生態(tài)服務(wù)價(jià)值對(duì)生態(tài)價(jià)值系數(shù)敏感度，即生態(tài)價(jià)值系數(shù)變動(dòng)1%對(duì)生態(tài)服務(wù)價(jià)值的影響狀況。當(dāng)敏感性系數(shù)（Cs）<1，說(shuō)明該生態(tài)服務(wù)類型沒(méi)有彈性，ESV對(duì)VC不敏感，當(dāng)Cs>1，表明該生態(tài)服務(wù)類型是比較有彈性的，ESV對(duì)VC敏感[17]。

根據(jù)公式（8）計(jì)算敏感性度數(shù)（Cs）：

Cs=（ESVj-ESVi）/ESVi（VCjk-VCik）/VCik。（8）

式中：Cs為敏感性度數(shù);ESVi為生態(tài)服務(wù)開(kāi)始的價(jià)值;ESVj為生態(tài)服務(wù)改變后的價(jià)值;VCik為價(jià)值系數(shù)改變前;VCjk為價(jià)值系數(shù)改變后;k為某種類別土地利用。

2 結(jié)果與分析

2.1 土地利用結(jié)構(gòu)變化

土地利用結(jié)構(gòu)是指某區(qū)域內(nèi)不同利用土地之間的占比關(guān)系，其囊括時(shí)空2個(gè)方面的變化[9]。貴港市土地利用結(jié)構(gòu)變化體現(xiàn)了該市在2010、2014、2018年這3期各個(gè)土地利用類型面積變化情況，可用來(lái)研究該市土地利用變化的趨勢(shì)和規(guī)律。運(yùn)用貴港市2010、2014、2018年土地利用現(xiàn)狀數(shù)據(jù)，采用Excel軟件，計(jì)算貴港市的“三生”土地利用類型面積。分析結(jié)果見(jiàn)表5、表6。

從表5可以看出，2010—2014年貴港市“三生”土地利用類型結(jié)構(gòu)情況：（1）農(nóng)業(yè)生產(chǎn)用地和林地生態(tài)用地在2010、2014年的土地利用類型中占據(jù)著很大的比例，說(shuō)明貴港市以這2種用地類型為主。其中，農(nóng)業(yè)生產(chǎn)用地主要分布在平南縣、覃塘區(qū)、桂平市，面積持續(xù)減少，從2010年的377 793.61 hm2減少至2014年的375 856.34 hm2，減少了 1 937.27 hm2。林地生態(tài)用地主要分布在平南縣、桂平市，從2010年的467 811.25 hm2減少至2014年的466 489.76 hm2，減少了1 321.49 hm2，這2種地類面積也是減少最多的。除了這2種地類減少外，還有牧草地生態(tài)用地、其他生態(tài)用地、水域生態(tài)用地，分別減少了479.85、68.09、315.08 hm2。但減少的面積相對(duì)前面2種來(lái)說(shuō)，都不是很大。（2）從

用地類型增加來(lái)看，服務(wù)、工礦這些生產(chǎn)用地和城鎮(zhèn)、農(nóng)村這些生活用地，分別增加了18.77、1 462.00、2 412.10、228.91 hm2。主要以工礦生產(chǎn)用地和城鎮(zhèn)生活用地增加為主。

從表6可以看出，2014—2018年貴港市“三生”土地利用類型結(jié)構(gòu)情況：（1）與表5的情況類似，農(nóng)業(yè)生產(chǎn)用地和林地生態(tài)用地依然是貴港市的主要用地類型，而且二者在2014—2018年減少的面積都很大，也是主要減少的2種用地類型;其他減少的用地類型為牧草地生態(tài)用地、其他生態(tài)用地、水域生態(tài)用地，但減少的面積相對(duì)來(lái)說(shuō)都不是很大。（2）從用地類型增加來(lái)看，也是主要以工礦生產(chǎn)用地和城鎮(zhèn)生活用地為主，另外2種也是和表5一致，分別為服務(wù)以及其他生產(chǎn)用地和農(nóng)村生活用地，面積雖有增加，但增加都不是很大。

綜上來(lái)說(shuō)，2010—2018年貴港市土地利用結(jié)構(gòu)發(fā)生明顯變化，工礦生產(chǎn)用地和城鎮(zhèn)生活用地的面積增加很多，而農(nóng)業(yè)、林地這些生態(tài)用地面積減少很多;服務(wù)以及其他生產(chǎn)用地、農(nóng)村生活用地雖有增加，但增加都不是很明顯;其他生態(tài)用地和水域生態(tài)用地雖有減少，但減少的面積不明顯。對(duì)于其他生態(tài)用地2018年在貴港市占比3.87%，且較為穩(wěn)定，說(shuō)明貴港市土地利用程度中等偏上，土地利用后備資源充足。

2.2 土地利用動(dòng)態(tài)度分析

2.2.1 單一土地利用動(dòng)態(tài)度 根據(jù)公式（1），本研究將以2010—2014年、2014—2018年這2個(gè)時(shí)段，分別研究其單一土地利用動(dòng)態(tài)度（表7）。

從表7可以看出，（1）農(nóng)業(yè)生產(chǎn)用地和林地、牧草地、水域、其他這些生態(tài)用地在這2個(gè)時(shí)段內(nèi)動(dòng)態(tài)度都為負(fù)數(shù)，說(shuō)明這些地類在2010—2018年是處于不斷減少的狀態(tài)。其中，我們可以發(fā)現(xiàn)，牧草地生態(tài)用地動(dòng)態(tài)度最小，說(shuō)明該地類減少的面積程度是最大的;而農(nóng)業(yè)生產(chǎn)用地和水域生態(tài)用地相同，說(shuō)明這2類地類的面積減少的程度相同，其他生態(tài)用地的動(dòng)態(tài)度減少的面積程度稍微少一些。（2）服務(wù)以及其他生產(chǎn)用地、工礦生產(chǎn)用地、城鎮(zhèn)生活用地

和農(nóng)村生活用地這些地類的動(dòng)態(tài)度在這2個(gè)時(shí)段都為正值，說(shuō)明這些地類在2010—2018年都處于增長(zhǎng)狀態(tài)。工礦生產(chǎn)用地和城鎮(zhèn)生活用地的動(dòng)態(tài)度都比較大，說(shuō)明這2類地類面積增加的程度都很大;而服務(wù)以及其他生產(chǎn)用地和農(nóng)村生活用地動(dòng)態(tài)度稍小一些，說(shuō)明這2類土地面積增加的程度要比前面2類少。

2.2.2 土地利用綜合動(dòng)態(tài)度 根據(jù)公式（2），計(jì)算出貴港市2010—2018年的土地利用綜合動(dòng)態(tài)度為0.220，但貴港市在2010—2014年、2014—2018年的綜合動(dòng)態(tài)度為0.097和0.130，土地利用綜合動(dòng)態(tài)度提高明顯。說(shuō)明貴港市土地利用隨著時(shí)間的延后，其改變速率也會(huì)跟著提升，土地利用轉(zhuǎn)型因人類活動(dòng)干擾而顯著上升。

2.2.3 土地利用程度變化 使用貴港市2010、2014、2018年的這3期土地利用現(xiàn)狀數(shù)據(jù)，以區(qū)（縣、市）為基本研究單元，運(yùn)用公式（3）、公式（4）、公式（5），分別計(jì)算出它們?cè)?010—2018年的綜合指數(shù)、變化量和變化率（表8）。

貴港市2010、2014、2018年的土地利用程度綜合指數(shù)分別是247.95、248.58、249.32。貴港市土地利用程度指數(shù)高于中國(guó)科學(xué)院遙感應(yīng)用研究所測(cè)定的中國(guó)總土地利用程度指數(shù)。由此可以得知，貴港市土地利用程度較高，人類活動(dòng)在一定程度上影響了該市的土地利用格局。貴港市2010—2018年指數(shù)提高了1.37，表明其土地利用程度不斷升高。

2.3 貴港市土地利用功能轉(zhuǎn)型分析

2.3.1 基于貴港市3區(qū)1縣1市土地利用功能轉(zhuǎn)型 2010—2018年，貴港市利用類型土地變化情況見(jiàn)圖3。

從圖3可以看出，在研究時(shí)間區(qū)段內(nèi)，貴港市3區(qū)1縣1市土地利用功能轉(zhuǎn)型顯著，他們都表現(xiàn)出共同特點(diǎn)，即生活用地都表現(xiàn)為增加狀態(tài)，其中增加最多的為港北區(qū)，共增加1 685.21 hm2;其次為桂平市增加1 511.82 hm2;最低為港南區(qū)，增加面積為874.44 hm2。由于經(jīng)濟(jì)的快速發(fā)展，人口增加不可避免地占用生產(chǎn)和生態(tài)用地;生態(tài)用地都呈現(xiàn)為下降狀態(tài)，其中桂平市減少的最多，為1 442.11 hm2，其次是平南縣，為1 183.93 hm2，最低港北區(qū)為 625.56 hm2;而生產(chǎn)用地呈現(xiàn)不同的變化狀態(tài)，港北區(qū)、港南區(qū)和覃塘區(qū)、桂平市呈現(xiàn)下降狀態(tài)，其中，港北區(qū)減少的最多，為1 059.65 hm2，其次是覃塘區(qū)為 509.73 hm2。平南縣呈現(xiàn)為上升的狀態(tài)。港北區(qū)、港南區(qū)和覃塘區(qū)都是生產(chǎn)用地和生態(tài)用地轉(zhuǎn)型為生活用地，平南縣則由生態(tài)用地轉(zhuǎn)型為生產(chǎn)用地和生活用地，桂平市則是由生產(chǎn)和生態(tài)用地轉(zhuǎn)型為生活用地。在研究時(shí)間區(qū)段內(nèi)，可以得知，港北區(qū)經(jīng)濟(jì)發(fā)展和人口增加速度是最快的。

2.3.2 貴港市土地利用功能轉(zhuǎn)型總體分析 運(yùn)用公式（6），可計(jì)算出各類用地情況（表9）。

本研究根據(jù)土地利用變化數(shù)據(jù)，計(jì)算出“生產(chǎn)、生活、生態(tài)”用地類型擴(kuò)展強(qiáng)度指數(shù)見(jiàn)表9，其體現(xiàn)的是各類土地的擴(kuò)展情況，農(nóng)業(yè)生產(chǎn)用地和林地、牧草地、其他、水域這些生態(tài)用地在2010—2014年、2014—2018年擴(kuò)展強(qiáng)度指數(shù)為負(fù)值，說(shuō)明這些用地變化特征主要表現(xiàn)為緊縮。其中牧草地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地、水域生態(tài)用地在2010—2014年、2014—2018年擴(kuò)展強(qiáng)度指數(shù)為負(fù)值且較大，說(shuō)明這些地類緊縮很快。服務(wù)以及其他、工礦這些生產(chǎn)用地和城鎮(zhèn)、農(nóng)村這些生活用地的擴(kuò)展指數(shù)都為正值，說(shuō)明這些用地變化特征表現(xiàn)為擴(kuò)展。其中城鎮(zhèn)生活用地和工礦生產(chǎn)用地在2010—2014年、2014—2018年擴(kuò)展強(qiáng)度指數(shù)為正值且較大，說(shuō)明這些用地?cái)U(kuò)展速度快。

2.4 生態(tài)服務(wù)價(jià)值變化分析

2.4.1 生態(tài)服務(wù)總價(jià)值變化 運(yùn)用上面的生態(tài)服務(wù)價(jià)值計(jì)算模型。計(jì)算出2010—2018年貴港市生態(tài)服務(wù)價(jià)值變化（表10）。

從表10可以看出，貴港市總的生態(tài)服務(wù)價(jià)值在2010、2014、2018年分別為14 507.37×106、14 425.81×106、14 362.02×106元，說(shuō)明2010—2018年貴港市生態(tài)服務(wù)價(jià)值呈現(xiàn)下降趨勢(shì)。2010—2018年，貴港市的生態(tài)服務(wù)價(jià)值減少了 145.35×106元，年均減少18.169×106元。主要原因是貴港市經(jīng)濟(jì)發(fā)展迅速，人口的激增，建成區(qū)域面積的變大，林地、農(nóng)田、水體及其他地類向建設(shè)用地轉(zhuǎn)化。

從表10可以看出，2010—2018年6種用地類型的生態(tài)服務(wù)價(jià)值排序情況。林地在貴港市占比是較大的，與農(nóng)田的占比相差不是很大，但林地生態(tài)服務(wù)價(jià)值貢獻(xiàn)卻遠(yuǎn)大于農(nóng)田生態(tài)服務(wù)價(jià)值總量;水體的占比也并沒(méi)有農(nóng)田占比的多，但水體的服務(wù)價(jià)值超過(guò)農(nóng)田的服務(wù)價(jià)值。貴港市生態(tài)價(jià)值主要貢獻(xiàn)者是林地，后為水體、農(nóng)田。然而林地、農(nóng)田在2010—2018年間分別減少了2 474.14、5 486.16 hm2。草地、未利用地對(duì)貴港市的生態(tài)服務(wù)功能價(jià)值貢獻(xiàn)率不大，特別是未利用地。因此，保護(hù)林地、水體、農(nóng)田具有十分重要的意義。

2.5 生態(tài)價(jià)值系數(shù)對(duì)生態(tài)服務(wù)價(jià)值敏感度分析

以公式（8）和表9確定的貴港市生態(tài)價(jià)值系數(shù)為基礎(chǔ)，分別將表4的各用地類生態(tài)價(jià)值系數(shù)增加或減少10%，計(jì)算出調(diào)整后的貴港市2010、2014、2018年生態(tài)服務(wù)價(jià)值量和敏感性度數(shù)（表11）。

從表11可以看出，2010—2018年林地生態(tài)敏感性度數(shù)最大，并出現(xiàn)升高的走向。林地生態(tài)價(jià)值系數(shù)每升高或下降1%，則貴港市生態(tài)服務(wù)總價(jià)值將升高或下降0.63%。從表11可以看出，貴港市從2010—2018年各地類敏感度雖有變化，但其敏感度次序并沒(méi)有發(fā)生變化。特別是建設(shè)用地、草地、未利用地均小于1，而且草地、未利用地接近于0，說(shuō)明所使用的系數(shù)符合貴港市的實(shí)際。貴港市生態(tài)服務(wù)價(jià)值對(duì)生態(tài)系數(shù)變化缺乏彈性，即貴港市生態(tài)服務(wù)價(jià)值對(duì)生態(tài)價(jià)值系數(shù)不敏感。

3 討論與結(jié)論

由于筆者所學(xué)知識(shí)水平有限，不能探討貴港市土地利用轉(zhuǎn)型的驅(qū)動(dòng)機(jī)制，希望未來(lái)在此基礎(chǔ)上，繼續(xù)深入研究貴港市土地利用轉(zhuǎn)型的驅(qū)動(dòng)機(jī)制。

參考的價(jià)值系數(shù)可能還不夠符合實(shí)際情況，計(jì)算出生態(tài)服務(wù)總價(jià)值可能稍微存在偏差。因此，希望在未來(lái)的研究中，確定更加合理符合該市的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)。

利用貴港市2010、2014、2018年這3期土地利用現(xiàn)狀數(shù)據(jù)進(jìn)行分析，按照“三生”分類標(biāo)準(zhǔn)，得出以下結(jié)論：（1）貴港市農(nóng)業(yè)生產(chǎn)用地、林地生態(tài)用地大幅度減少，而城鎮(zhèn)生活用地、工礦生產(chǎn)用地大幅增加，但貴港市仍以農(nóng)業(yè)生產(chǎn)用地、林地生態(tài)用地為主。（2）貴港市的其他生態(tài)用地占比不少，表明貴港市土地資源開(kāi)發(fā)利用程度中等偏上，土地開(kāi)發(fā)后備資源充足。（3）在研究時(shí)間段2010—2018年期間，土地利用綜合動(dòng)態(tài)度提高較大。說(shuō)明貴港市土地利用隨著時(shí)間的延后，其改變速率也會(huì)跟著提升，土地利用轉(zhuǎn)型因人類活動(dòng)干擾而顯著上升。（4）在研究時(shí)間段2010—2018年期間，貴港市生態(tài)服務(wù)總價(jià)值逐漸降低，其貢獻(xiàn)者主要是林地，后為水體、農(nóng)田。

參考文獻(xiàn)：

[1]龍花樓，李秀彬. 長(zhǎng)江沿線樣帶農(nóng)村宅基地轉(zhuǎn)型[J]. 地理學(xué)報(bào)，2005，60（2）：179-188.

[2]Grainger A. The future role of the tropical rain forests in the world forest economy[D]. Oxford：University of Oxford，1986.

[3]Vitousek P M，Harold A M，Lubchenco J，et al. Human domination of earths ecosystems[J]. Science，1997，277（5325）：494-499.

[4]Izquierdo A E，Grau H R. Agriculture adjustment，land-use transition and protected areas in Northwestern Argentina[J]. Journal of Environmental Management，2009，90（2）：858-865.

[5]Lambin E F，Meyfroidt P. Land use transitions：socio-ecological feedback versus socio-economic change[J]. Land Use Policy，2010，27（2）：108-118.

[6]蔡運(yùn)龍. 土地利用/土地覆蓋變化研究：尋求新的綜合途徑[J]. 地理研究，2001，20（6）：645-652.

[7]龍花樓. 土地利用轉(zhuǎn)型——土地利用/覆被變化綜合研究的新途徑[J]. 地理與地理信息科學(xué)，2003，19（1）：87-90.

[8]龍花樓. 論土地利用轉(zhuǎn)型與土地資源管理[J]. 地理研究，2015，34（9）：1607-1618.

[9]李 鑫，歐名豪，劉建生，等. 基于不確定性理論的區(qū)域土地利用結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30（4）：176-184.

[10]謝 菲，舒曉波，廖富強(qiáng)，等. 浮梁縣土地利用變化及驅(qū)動(dòng)力分析[J]. 水土保持研究，2011，18（2）：213-217，221.

[11]莊大方，劉紀(jì)遠(yuǎn). 中國(guó)土地利用程度的區(qū)域分異模型研究[J]. 自然資源學(xué)報(bào)，1997，12（2）：10-16.

[12]胡和兵，劉紅玉，郝敬鋒，等. 南京市九鄉(xiāng)河流域土地利用程度空間異質(zhì)性分析[J]. 地球信息科學(xué)學(xué)報(bào)，2012，14（5）：627-634.

[13]胡和兵，劉紅玉，郝敬鋒，等. 城市化流域生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空分異特征及其對(duì)土地利用程度的響應(yīng)[J]. 生態(tài)學(xué)報(bào)，2013，33（8）：2565-2576.

[14]劉盛和，吳傳鈞，沈洪泉. 基于GIS的北京城市土地利用擴(kuò)展模式[J]. 地理學(xué)報(bào)，2000，55（4）：407-416.

[15]歐陽(yáng)志云，王如松，趙景柱. 生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價(jià)值估算[J]. 應(yīng)用生態(tài)學(xué)報(bào)，1999，10（5）：607-613.

[16]喻建華，高中貴，張 露，等. 昆山市生態(tài)系統(tǒng)服務(wù)價(jià)值變化研究[J]. 長(zhǎng)江流域資源與環(huán)境，2005，14（2）：213-217.

[17]謝高地，甄 霖，魯春霞，等. 一個(gè)基于專家知識(shí)的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法[J]. 自然資源學(xué)報(bào)，2008，23（5）：911-919.