文 | 劉澤琴 米夢(mèng)凡 楊秀峰 趙子輝 解春苓

1. 北華航天工業(yè)學(xué)院經(jīng)濟(jì)管理學(xué)院 2. 北華航天工業(yè)學(xué)院遙感信息工程學(xué)院 3. 河北省航天遙感信息處理與應(yīng)用協(xié)同創(chuàng)新中心

一、前言

在黨的組織結(jié)構(gòu)和國(guó)家政權(quán)結(jié)構(gòu)中，縣一級(jí)處在承上啟下的關(guān)鍵環(huán)節(jié)，是發(fā)展經(jīng)濟(jì)、保障民生、維護(hù)穩(wěn)定的重要基礎(chǔ)。縣域經(jīng)濟(jì)帶有城鎮(zhèn)經(jīng)濟(jì)和農(nóng)村經(jīng)濟(jì)相結(jié)合的典型特征，在城鎮(zhèn)化發(fā)展進(jìn)程中的土地資源合理利用是具有重要研究?jī)r(jià)值的關(guān)鍵問(wèn)題之一。如何利用遙感、GIS等現(xiàn)代技術(shù)手段，對(duì)縣域土地資源進(jìn)行定性分析、精準(zhǔn)施策、科學(xué)管理，提升土地利用的效益，是縣、鄉(xiāng)政府等決策部門亟待思考的問(wèn)題。

在區(qū)域經(jīng)濟(jì)研究中引入空間地理信息變量是近年來(lái)獲得迅速發(fā)展的經(jīng)濟(jì)地理學(xué)的一個(gè)熱門研究領(lǐng)域。作為經(jīng)濟(jì)學(xué)與遙感學(xué)、地理學(xué)相互融合的一門新興學(xué)科，經(jīng)濟(jì)地理學(xué)在產(chǎn)業(yè)經(jīng)濟(jì)和區(qū)域經(jīng)濟(jì)研究領(lǐng)域受到了越來(lái)越高的重視，其中土地利用與覆被變化（LUCC）屬于熱點(diǎn)領(lǐng)域。這類包括空間變量的模型主要特點(diǎn)是采用遙感信息作為數(shù)據(jù)來(lái)源，并且利用計(jì)算機(jī)技術(shù)和GIS，分析LUCC的空間顯性分布特征和格局構(gòu)成，同時(shí)重視對(duì)影響空間差異性的自然、人類活動(dòng)等驅(qū)動(dòng)因子的分析。目前廣泛采用的分析模型在原理構(gòu)造上大都源于地理元胞自動(dòng)機(jī)（CA）。CLUE-S是源于CA原理的一種比較成熟的模型形式。區(qū)域經(jīng)濟(jì)研究的另一個(gè)熱點(diǎn)領(lǐng)域是土地利用的效益評(píng)價(jià)，涉及生態(tài)（環(huán)境）效益、經(jīng)濟(jì)效益和社會(huì)效益等方面。

空間地理信息分布及格局演化方面的技術(shù)方法目前已取得較大進(jìn)展，空間地理模型分析早已有之，但是將其研究結(jié)果應(yīng)用于地方政府治理和民間應(yīng)用的案例不多，應(yīng)用于較小區(qū)域開展經(jīng)濟(jì)分析的研究數(shù)量更少，尤其是針對(duì)縣域這類國(guó)民經(jīng)濟(jì)基本單元的實(shí)證研究比較鮮見(jiàn)?？h域經(jīng)濟(jì)所處的地位及其結(jié)構(gòu)具有重要研究?jī)r(jià)值，本文力圖在此方面做出探索。

二、研究區(qū)、數(shù)據(jù)及研究方法

1. 研究區(qū)概況

巨鹿縣位于河北省中南部，邢臺(tái)市中部，地處太行山東麓，華北平原的邊緣，古黃河、漳河沖積平原上，地勢(shì)平坦。地理位置為37°N和115°E，四周與邢臺(tái)市7個(gè)縣交界，縣域面積631km。該縣總?cè)丝?2萬(wàn)，下轄6個(gè)鎮(zhèn)，4個(gè)鄉(xiāng)，縣城駐地巨鹿鎮(zhèn)。巨鹿縣交通便捷，公路總里程逾1400km，邢衡高速貫穿全境，國(guó)家二級(jí)公路（省道）南郝線、邢德線橫穿東西，定魏線縱貫?zāi)媳?。?dāng)?shù)剞r(nóng)業(yè)發(fā)展具有特色，是中國(guó)金銀花道地藥材之鄉(xiāng)、中國(guó)枸杞產(chǎn)業(yè)基地、中國(guó)串枝紅杏之鄉(xiāng)，中草藥材和園林水果種植兩項(xiàng)產(chǎn)值在農(nóng)業(yè)總產(chǎn)值中占比近1/4。面對(duì)土地資源以耕地為主、城鎮(zhèn)化帶來(lái)的建設(shè)用地?cái)U(kuò)張、新興的二、三產(chǎn)業(yè)發(fā)展對(duì)土地需求增加等土地利用的多重矛盾，如何實(shí)現(xiàn)土地節(jié)約、集約利用，提升土地利用的經(jīng)濟(jì)效益，是當(dāng)?shù)卣叫杵平獾碾y題。

2. 數(shù)據(jù)

本文基礎(chǔ)數(shù)據(jù)由衛(wèi)星遙感信息和經(jīng)濟(jì)社會(huì)發(fā)展統(tǒng)計(jì)資料兩部分組成。

遙感數(shù)據(jù)用于提取土地利用現(xiàn)狀和交通數(shù)據(jù)，根據(jù)研究目標(biāo)選擇了2010、2014和2018年三期遙感影像，其中2010年遙感數(shù)據(jù)來(lái)源于谷歌影像，數(shù)據(jù)精度7.36m；2014年和2018年遙感影像來(lái)源于高分一號(hào)衛(wèi)星，數(shù)據(jù)精度2m。影像時(shí)間均選取當(dāng)年7-8月份，數(shù)據(jù)質(zhì)量好。土地利用分類處理采用ArcGIS軟件通過(guò)人工解譯進(jìn)行矢量化，土地利用類型包括七種，分別為耕地、林地、城鎮(zhèn)建設(shè)用地、農(nóng)村居民點(diǎn)、水域、工礦企業(yè)用地和其他用地（主要是交通設(shè)施和未利用地）（表1）。圖1為150m×150m尺度的三期土地利用情況?？h域內(nèi)交通設(shè)施包括省道、縣道、鄉(xiāng)村道路三類道路（線狀地理信息）（2014年10月邢衡高速公路邢臺(tái)段通車，鑒于本文所涉時(shí)間范圍自2010年開始，因此未將該高速公路和巨鹿縣高速站口作為交通因素納入研究范圍），以及縣客運(yùn)總站（點(diǎn)狀地理信息）。利用歐式距離法計(jì)算縣域內(nèi)各位置至各交通設(shè)施的距離。

表1 土地利用類型及代碼

圖1 巨鹿縣2010年、2014年、2018年土地利用狀況

經(jīng)濟(jì)社會(huì)發(fā)展數(shù)據(jù)來(lái)源于2011-2018年《河北經(jīng)濟(jì)年鑒》《邢臺(tái)統(tǒng)計(jì)年鑒》和《巨鹿統(tǒng)計(jì)月報(bào)》，涉及國(guó)民經(jīng)濟(jì)核算、農(nóng)村經(jīng)濟(jì)、工業(yè)等部分，所有指標(biāo)均采用現(xiàn)價(jià)（當(dāng)年價(jià)）口徑。在對(duì)2022年交通因素的變化推導(dǎo)中，參考了《巨鹿縣城鄉(xiāng)總體規(guī)劃（2013-2030）》，《巨鹿縣國(guó)民經(jīng)濟(jì)與社會(huì)發(fā)展第十三個(gè)五年規(guī)劃綱要》，以及巨鹿縣近年《政府工作報(bào)告》。

3. 研究方法

（1）土地利用變化模型CLUE-S

CLUE-S模型是由CLUE模型演化而來(lái)的，更加適于開展較小區(qū)域的研究。對(duì)于CLUE-S模型的擬合效果檢驗(yàn)多采用Kappa系數(shù)，一般認(rèn)為該系數(shù)值≥0.75時(shí)，表示模型擬合效果較好。

（2）經(jīng)濟(jì)效益評(píng)價(jià)模型

土地利用經(jīng)濟(jì)效益評(píng)價(jià)主要利用研究區(qū)經(jīng)濟(jì)社會(huì)發(fā)展統(tǒng)計(jì)數(shù)據(jù)，如地區(qū)生產(chǎn)總值及其構(gòu)成、主要產(chǎn)業(yè)產(chǎn)值等，結(jié)合各種土地利用類型的用地量數(shù)據(jù)進(jìn)行簡(jiǎn)單匯總。

三、結(jié)果與分析

1. 2010-2018年巨鹿縣土地利用模擬

（1）Logistic回歸分析結(jié)果

CLUE-S模型中，柵格大小一般應(yīng)在100～500m之間，結(jié)合研究區(qū)的規(guī)模和經(jīng)濟(jì)特征，本文采用150m×150m的柵格進(jìn)行分析。采用SPSS 22中的Binary Logistic回歸模型研究2010-2014年間交通因素對(duì)土地利用類型的驅(qū)動(dòng)關(guān)系，分析得到回歸系數(shù)。就結(jié)果來(lái)看，僅有農(nóng)村居民點(diǎn)的模擬效果欠佳，其他用地類型的ROC曲線下面積均超過(guò)0.6。因此Logistic回歸效果基本可以滿足當(dāng)前研究需要，可將其作為CLUE-S模型分析的位置布局參數(shù)。

（2）地類間轉(zhuǎn)換規(guī)則設(shè)置

針對(duì)巨鹿縣2022年的土地利用情況模擬預(yù)測(cè)，本文擬作三種不同的假設(shè)情境設(shè)定，分別為：“自然發(fā)展”情境，即假設(shè)當(dāng)?shù)乩^續(xù)保持以往的發(fā)展方式，不對(duì)土地利用施加人為干預(yù)；“耕地保護(hù)”情境，即假設(shè)當(dāng)?shù)氐母乜偭勘３植蛔?，并令已有耕地向其他用地類型轉(zhuǎn)換的難度加大，以保持和充分發(fā)揮巨鹿縣在特色農(nóng)產(chǎn)品種植業(yè)方面的優(yōu)勢(shì)；“水域保護(hù)”情境，即假設(shè)當(dāng)?shù)厮蛎娣e保持不變，并令已有耕地向其他用地類型轉(zhuǎn)換的難度加大，以適應(yīng)巨鹿縣“千湖之縣”發(fā)展思路的需要。

土地利用轉(zhuǎn)移彈性用于表示某種用地類型轉(zhuǎn)換為其他用地類型的難易程度，取值為0～1，其中0表示可轉(zhuǎn)換為任意其他地類，1表示不允許轉(zhuǎn)換為其他地類，該值越大表明轉(zhuǎn)換成其他用地類型的難度越高。針對(duì)這三種假設(shè)情境，結(jié)合各種用地類型自身的特征，分別設(shè)置七種地類的土地利用轉(zhuǎn)移彈性。土地利用轉(zhuǎn)移矩陣用于表示某種用地類型轉(zhuǎn)換為其他用地類型的可能性，取值為0或1，其中0表示不可轉(zhuǎn)換，1表示可以轉(zhuǎn)換。

（3）基于CLUE-S的2018年土地利用空間模擬

為了檢驗(yàn)上述參數(shù)設(shè)置的合理性，本文以2018年巨鹿縣遙感影像解譯得到的150m×150m尺度的土地利用情況為目標(biāo)時(shí)期的土地利用需求，見(jiàn)圖1（c），2014年的相應(yīng)數(shù)據(jù)作為基期，通過(guò)線性插值法確定中間年份各種地類的用地量，采用CLUE-S模擬“自然發(fā)展”情境下2018年的土地利用情況（鑒于巨鹿縣當(dāng)?shù)匚丛O(shè)有自然保護(hù)區(qū)或其他明確限制開發(fā)區(qū)域，因此本文在研究中未設(shè)置空間政策與約束條件），得到模擬正確的柵格數(shù)為25544個(gè)，與2018年150m×150m尺度下總柵格數(shù)28019個(gè)相比，得到Kappa系數(shù)為91.17%，表明CLUE-S模型模擬效果好，以上參數(shù)設(shè)置合理，可以將其代入2018-2022年分析周期進(jìn)行2022年土地利用情況的預(yù)測(cè)。

2. 不同假設(shè)情境下2022年巨鹿縣土地利用預(yù)測(cè)

（1）2022年巨鹿縣土地利用需求預(yù)測(cè)

對(duì)土地利用需求的預(yù)測(cè)可以采用簡(jiǎn)單趨勢(shì)外推法，也可以采用形式更為復(fù)雜的經(jīng)濟(jì)模型進(jìn)行分析。本文重點(diǎn)關(guān)注了在縣域經(jīng)濟(jì)發(fā)展中貢獻(xiàn)顯著、同時(shí)用地?cái)?shù)量變化較大的第二產(chǎn)業(yè)發(fā)展情況，以第二產(chǎn)業(yè)增加值作為其發(fā)展水平的測(cè)算指標(biāo)。接下來(lái)本文采用巨鹿縣第二產(chǎn)業(yè)增加值（X）和工礦企業(yè)用地面積（Y）進(jìn)行趨勢(shì)擬合的方式預(yù)測(cè)2022年該土地利用類型的土地需求量。此外，在“耕地保護(hù)”情境下假設(shè)耕地面積保持與2018年相同，在“水域保護(hù)”情境下假設(shè)水域面積保持與2018年相同。在保持土地總量不變的前提下，其他用地類型的土地需求量繼續(xù)延續(xù)原有趨勢(shì)并作適當(dāng)調(diào)整。

（2）2022年巨鹿縣土地利用預(yù)測(cè)

采用交通因素驅(qū)動(dòng)下土地利用類型的Logistic回歸系數(shù)和土地利用轉(zhuǎn)移矩陣，根據(jù)三種假設(shè)情境，分別采用不同的土地利用轉(zhuǎn)移彈性結(jié)合土地利用需求量，以2018年土地利用狀況為基礎(chǔ)，運(yùn)行CLUE-S模型得到不同情境下的2022年巨鹿縣土地利用預(yù)測(cè)，見(jiàn)圖2。

（3）土地利用預(yù)測(cè)結(jié)果對(duì)比分析

圖2 巨鹿縣2022年土地利用預(yù)測(cè)

根據(jù)不同假設(shè)情境開展的2022年巨鹿縣土地利用預(yù)測(cè)結(jié)果顯示，在以下兩個(gè)方面存在比較顯著的差異。

土地利用結(jié)構(gòu)的數(shù)量差異?？v觀2010年以來(lái)的土地利用變化，可以發(fā)現(xiàn)耕地?cái)?shù)量變化的絕對(duì)值最大，對(duì)林地、工礦企業(yè)用地等用地量的增加貢獻(xiàn)最顯著。而在土地需求預(yù)測(cè)中，本文對(duì)工礦企業(yè)用地基于經(jīng)濟(jì)模型做出了數(shù)量規(guī)定，相當(dāng)于設(shè)置了較具剛性的土地需求條件。由于“耕地保護(hù)”情境假設(shè)耕地?cái)?shù)量維持在2018年現(xiàn)狀水平，直接導(dǎo)致的結(jié)果是林地?zé)o法有效增加，因此圖2(b)顯示的結(jié)果與圖1(c)差異最小。在“水域保護(hù)”情境下的情況類似，不過(guò)由于水域面積在土地總量中占比很小（2018年為5.85%），對(duì)其施加的數(shù)量限制條件影響較小，因此該情境與“自然發(fā)展”情境的土地利用需求相似度較高。

土地利用格局的空間差異?！白匀话l(fā)展”和“水域保護(hù)”兩種情境與2018年土地利用現(xiàn)狀的差異最為顯著，共同的變化特征有兩個(gè)：要求將目前縣城北部和東部周邊尤其是東北方向的一些耕地轉(zhuǎn)為林地，形成縣城以南和以北兩個(gè)林木富集區(qū)域；要求在目前縣城西北部和東北部部分地區(qū)安排工礦企業(yè)用地?cái)U(kuò)大，形成連片規(guī)模發(fā)展。此外，三種假設(shè)情境都要求將西部的西郭城鎮(zhèn)靠近與臨縣交界區(qū)（臨近省道S324邢德線）、蘇家營(yíng)鄉(xiāng)與臨縣交界區(qū)和官亭鎮(zhèn)臨近省道S327南郝線部分區(qū)域的農(nóng)村居民點(diǎn)轉(zhuǎn)為工礦企業(yè)用地，充分發(fā)揮交通設(shè)施助力企業(yè)發(fā)展的作用。

3. 巨鹿縣土地利用經(jīng)濟(jì)效益評(píng)價(jià)

考慮到各種土地利用類型與國(guó)民經(jīng)濟(jì)產(chǎn)業(yè)的對(duì)應(yīng)性，以及縣域統(tǒng)計(jì)資料的可得性，本文選取（狹義）農(nóng)業(yè)產(chǎn)值作為耕地的產(chǎn)出指標(biāo)，林業(yè)產(chǎn)值、第三產(chǎn)業(yè)增加值、農(nóng)林牧漁服務(wù)業(yè)產(chǎn)值、漁業(yè)產(chǎn)值和規(guī)模以上工業(yè)總產(chǎn)值分別作為林地、城鎮(zhèn)建設(shè)用地、農(nóng)村居民點(diǎn)、水域和工礦建設(shè)用地的產(chǎn)出指標(biāo)。所有指標(biāo)均采用現(xiàn)價(jià)（當(dāng)年價(jià)）口徑統(tǒng)計(jì)數(shù)據(jù)。根據(jù)2010-2018年統(tǒng)計(jì)資料和相應(yīng)地類的用地量，計(jì)算得到各地類的單位面積經(jīng)濟(jì)效益值，根據(jù)2010-2018年各地類單位面積經(jīng)濟(jì)效益值，通過(guò)趨勢(shì)外推法得到2022年各指標(biāo)預(yù)測(cè)值。以預(yù)測(cè)的2022年單位面積經(jīng)濟(jì)效益值分別結(jié)合三種假設(shè)情境下的各地類用地量，可以計(jì)算得到不同情境下的經(jīng)濟(jì)效益總量預(yù)測(cè)值（表2）。

表2 土地利用經(jīng)濟(jì)效益總量 單位：萬(wàn)元

由表2顯示的三種不同假設(shè)情境下巨鹿縣土地利用經(jīng)濟(jì)效益總量來(lái)看，“自然發(fā)展”情境的經(jīng)濟(jì)效益顯著高于其他兩種情境，這與該情境下單位面積經(jīng)濟(jì)效益值最高的工礦企業(yè)用地比其他二者多有密切關(guān)系。其他兩種假設(shè)情境之間的經(jīng)濟(jì)效益總量差異不顯著，雖然兩種情境下的各地類用地量存在差別，但是結(jié)合各種地類的單位面積經(jīng)濟(jì)效益值之后綜合來(lái)看，經(jīng)濟(jì)效果上是相當(dāng)?shù)?。如果僅從經(jīng)濟(jì)效益總量的角度出發(fā)考慮，當(dāng)?shù)厥紫葢?yīng)當(dāng)考慮選擇“自然發(fā)展”模式，但是如果補(bǔ)充更多的生態(tài)效益和社會(huì)效益資料，則綜合評(píng)價(jià)的結(jié)果可能會(huì)發(fā)生變化。

四、巨鹿縣土地利用及產(chǎn)業(yè)發(fā)展策略啟示

采用交通設(shè)施作為驅(qū)動(dòng)因素，運(yùn)用CLUE-S模型對(duì)2022年巨鹿縣土地利用情況進(jìn)行模擬預(yù)測(cè)，對(duì)比2018年當(dāng)?shù)赝恋乩矛F(xiàn)狀發(fā)現(xiàn)，工礦企業(yè)用地?cái)U(kuò)大和林地增加是預(yù)測(cè)中出現(xiàn)的主要特征，并且工礦企業(yè)的用地指向與當(dāng)?shù)爻青l(xiāng)總體規(guī)劃中指出的未來(lái)中長(zhǎng)期發(fā)展思路是一致的，但是用地類型的布點(diǎn)更加具體，體現(xiàn)了科學(xué)性和應(yīng)用性的結(jié)合。

研究發(fā)現(xiàn)，在“自然發(fā)展”和“水域保護(hù)”兩種假設(shè)情境下，應(yīng)在目前縣城東北方向成規(guī)模地建設(shè)林地，這與2010-2018年間巨鹿縣林地面積顯著增加有關(guān)，也與模型的驅(qū)動(dòng)因素、限制條件等有關(guān)。這一發(fā)展方向與當(dāng)?shù)亟ㄔO(shè)“華北地區(qū)最綠最美平原縣”的發(fā)展目標(biāo)一致，也提示了林地建設(shè)的布局方向。但是就當(dāng)前林業(yè)產(chǎn)值較低的現(xiàn)狀來(lái)看，圍繞林業(yè)資源開發(fā)休閑旅游等第三產(chǎn)業(yè)是可以考慮的發(fā)展策略。

“耕地保護(hù)”和“水域保護(hù)”兩種假設(shè)情境下，經(jīng)濟(jì)效益總量與“自然發(fā)展”情境存在一定差距，這固然與評(píng)價(jià)模型選擇的經(jīng)濟(jì)指標(biāo)有關(guān)，但是也表明將種植業(yè)、漁業(yè)等囿于以往的第一產(chǎn)業(yè)部門是不夠的，應(yīng)該充分挖掘相應(yīng)自然資源、經(jīng)濟(jì)資源的潛力，將第一產(chǎn)業(yè)與工業(yè)、更重要地是與第三產(chǎn)業(yè)充分融合，將巨鹿縣道地中藥材產(chǎn)地的優(yōu)勢(shì)、河湖坑塘密集的優(yōu)勢(shì)充分發(fā)揮出來(lái)，走出一條特色農(nóng)業(yè)發(fā)展之路。

五、結(jié)論

采用空間地理信息尤其是衛(wèi)星遙感信息與經(jīng)濟(jì)效益分析模型相結(jié)合的方法，對(duì)縣域土地利用及其效益開展模擬和評(píng)價(jià)，這一思路具有可行性。這種研究方法突破了以往大量經(jīng)濟(jì)分析中僅僅依靠經(jīng)濟(jì)社會(huì)統(tǒng)計(jì)信息作為數(shù)據(jù)源的思路，將空間地理信息引入分析框架，可以充分發(fā)揮此類信息實(shí)時(shí)性、翔實(shí)性、動(dòng)態(tài)性等優(yōu)勢(shì)，不僅實(shí)現(xiàn)了基礎(chǔ)數(shù)據(jù)的豐富性，更有利于提升經(jīng)濟(jì)分析的質(zhì)量。

在本文所述新的研究框架下，模型分析所提供的結(jié)果更加豐富。根據(jù)不同的假設(shè)情境，可以模擬未來(lái)土地利用的多種格局，并且以圖形的形式呈現(xiàn)出來(lái)，結(jié)果的直觀性更好，比以往常用的經(jīng)濟(jì)研究方法提供的結(jié)果更易于理解。結(jié)合經(jīng)濟(jì)效益分析之后，各種不同假設(shè)情境下的經(jīng)濟(jì)效益及其所對(duì)應(yīng)的產(chǎn)業(yè)發(fā)展策略更加明晰，政策指導(dǎo)性更強(qiáng)。

根據(jù)本文所采用研究方法所提供的研究結(jié)果，對(duì)于地方政府調(diào)整土地利用格局、制定地方經(jīng)濟(jì)社會(huì)發(fā)展規(guī)劃、調(diào)整產(chǎn)業(yè)發(fā)展政策等均具有良好的參考價(jià)值。此類分析的結(jié)果明確、直觀、指向性強(qiáng)，易于轉(zhuǎn)化為具體的政策措施，且研究所耗費(fèi)的時(shí)間、經(jīng)濟(jì)成本均較低，可以提高政府治理效率。