馬彩虹， 安斯文， 文 琦， 滑雨琪， 李聰慧

（寧夏大學(xué)地理科學(xué)與規(guī)劃學(xué)院，寧夏 銀川 750021）

國(guó)土空間系統(tǒng)是土地自然生態(tài)要素和人類(lèi)社會(huì)經(jīng)濟(jì)要素相互作用下形成的動(dòng)態(tài)復(fù)雜巨系統(tǒng)［1］，依據(jù)功能可分為生產(chǎn)、生活和生態(tài)“三生”空間［2］。改革開(kāi)放以來(lái)，隨著人口增長(zhǎng)和快速城市化的推進(jìn)，三生空間沖突日益嚴(yán)重［3-4］。黨的十八大報(bào)告提出了“生產(chǎn)空間集約高效、生活空間宜居適度、生態(tài)空間山清水秀”的目標(biāo)框架，掀起了三生空間的研究熱潮［5-6］。三生空間演變是自然地理因素和人類(lèi)活動(dòng)共同作用的結(jié)果［7］，不同尺度和不同地域類(lèi)型下驅(qū)動(dòng)因子和動(dòng)力機(jī)制差異顯著［8-9］。在流域、省市縣等尺度下，人口增長(zhǎng)、經(jīng)濟(jì)發(fā)展和技術(shù)進(jìn)步是較為普遍的驅(qū)動(dòng)因子，如趙瑞等［10］通過(guò)主成分分析法得出人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展、農(nóng)業(yè)技術(shù)進(jìn)步是推動(dòng)京津冀都市圈三生空間變化的主因；宋永永等［11］運(yùn)用地理探測(cè)器中的因子探測(cè)和交互作用探測(cè)模塊，發(fā)現(xiàn)城鎮(zhèn)化、工業(yè)化及生態(tài)政策是黃河流域國(guó)土空間格局演化的關(guān)鍵動(dòng)力；甘娜等［12］對(duì)川滇藏交界區(qū)域三生空間重構(gòu)的驅(qū)動(dòng)因子進(jìn)行分析，認(rèn)為在生態(tài)脆弱區(qū)或重要生態(tài)功能區(qū)宏觀政策的驅(qū)動(dòng)效應(yīng)顯著；董建紅等［13］運(yùn)用偏最小二乘法揭示三生空間的驅(qū)動(dòng)機(jī)制，發(fā)現(xiàn)生態(tài)政策、人口增長(zhǎng)、對(duì)外貿(mào)易等是甘肅省國(guó)土空間格局調(diào)整的重要驅(qū)動(dòng)力；焦庚英等［14］通過(guò)冗余分析方法對(duì)婺源縣的研究表明，旅游產(chǎn)業(yè)的快速發(fā)展是該縣三生空間格局演變的核心驅(qū)動(dòng)力；王全喜等［15］采用雙變量局部空間自相關(guān)分析法探究了武山縣的國(guó)土空間格局演變，指出區(qū)域一體化發(fā)展在武山縣三生格局演變中起重要作用。村域尺度下三生空間的演變與宏觀和中觀尺度明顯不同，如于正松等［16］運(yùn)用參與式農(nóng)村評(píng)估（PRA）法對(duì)河南省曲溝鎮(zhèn)的研究中發(fā)現(xiàn)，能人的帶動(dòng)效應(yīng)和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整對(duì)于工業(yè)鎮(zhèn)的三生空間演變驅(qū)動(dòng)效應(yīng)顯著；范藝娜等［17］采用半結(jié)構(gòu)式訪談和參與式社區(qū)資源圖法開(kāi)展村域尺度三生空間演變分析，認(rèn)為惠農(nóng)政策和農(nóng)業(yè)技術(shù)創(chuàng)新是影響村域尺度三生空間格局演變的關(guān)鍵驅(qū)動(dòng)力。

綜上所述，學(xué)界采用主成分分析法、地理探測(cè)器、冗余分析法、偏最小二乘法、半結(jié)構(gòu)訪談法等多種方法，開(kāi)展了不同尺度的三生空間演變的驅(qū)動(dòng)力研究，為該領(lǐng)域的研究奠定了良好的基礎(chǔ)［18-21］，但仍存在需要探索的方面：一是這些方法各有優(yōu)勢(shì)，但總體來(lái)說(shuō)屬于推斷性估計(jì)類(lèi)分析法，研究結(jié)果存在不確定性，也給不同時(shí)空尺度的研究結(jié)果比較帶來(lái)一定困難；二是對(duì)驅(qū)動(dòng)因子的識(shí)別較多，動(dòng)力機(jī)制有待加強(qiáng)［22-24］。馬彩虹等［25］提出的土地利用轉(zhuǎn)移流研究法，綜合考慮土地利用變化的矢量屬性，能夠定量追蹤不同時(shí)期土地利用變化的數(shù)量、方向和結(jié)構(gòu)特征［26］，在國(guó)土空間演變驅(qū)動(dòng)力分析方面具有獨(dú)特優(yōu)勢(shì)。本文采用土地利用轉(zhuǎn)移流進(jìn)行寧夏固原市原州區(qū)國(guó)土空間格局演變的驅(qū)動(dòng)力溯源分析，以期為國(guó)土空間格局演變驅(qū)動(dòng)力分析探索新的思路與方法。

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

1.1 研究區(qū)概況

原州區(qū)地處寧夏南部，曾名固原縣（2002 年更名），地理位置位于106°00′～106°30′E、35°45′～36°31′N(xiāo)，轄7 鎮(zhèn)4 鄉(xiāng)3 個(gè)街道辦事處，面積2739.01 km2（圖1），屬半濕潤(rùn)半干旱氣候帶、北方農(nóng)牧交錯(cuò)帶。依據(jù)地形特征可劃分為六盤(pán)山陰濕區(qū)（面積占33.2%）、清水河谷平原區(qū)（20.5%）和黃土丘陵溝壑區(qū)（46.3%）。受陡坡開(kāi)荒、林木采伐、草原過(guò)牧及礦山開(kāi)采等人類(lèi)不合理活動(dòng)的影響，水土流失現(xiàn)象突出，旱災(zāi)、凍害、雹災(zāi)等自然災(zāi)害頻繁。農(nóng)戶生計(jì)脆弱，以苦甲天下著稱(chēng)，人地矛盾突出，被聯(lián)合國(guó)教科文組織認(rèn)為是人類(lèi)不宜生存的地區(qū)，是我國(guó)“三西”（河西、定西、西海固）扶貧的重點(diǎn)區(qū)域，也是國(guó)家精準(zhǔn)扶貧時(shí)期的深度貧困區(qū)縣，于2020年3月脫貧摘帽。改革開(kāi)放以來(lái)，經(jīng)濟(jì)有了較快發(fā)展，城鎮(zhèn)化水平有所提升，但仍存在競(jìng)爭(zhēng)力不高、無(wú)序開(kāi)發(fā)擠占生態(tài)空間、土地利用效率低下等問(wèn)題。

圖1 研究區(qū)概況Fig.1 Overview of the study area

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

2000 年7—8 月的TM 影像數(shù)據(jù)來(lái)源于中國(guó)地理空間數(shù)據(jù)云，使用ArcGIS 10.6軟件對(duì)校正后的數(shù)字影像進(jìn)行拼接、分割和融合處理，結(jié)合野外勘驗(yàn)資料，建立解譯標(biāo)志，人機(jī)交互方式解譯，并對(duì)結(jié)果數(shù)據(jù)進(jìn)行精度檢驗(yàn)，檢驗(yàn)精度達(dá)85%以上；2018 年數(shù)據(jù)為地理國(guó)情矢量數(shù)據(jù)。由于行政區(qū)劃調(diào)整，文中涉及地名和區(qū)劃范圍以原州區(qū)行政區(qū)劃現(xiàn)狀為準(zhǔn)。

1.3 研究方法

1.3.1 三生用地分類(lèi)體系參照主導(dǎo)功能法分類(lèi)體系［2,27-29］，結(jié)合《土地利用現(xiàn)狀分類(lèi)》GB/T 21010-2007 和GB/T 21010-2017、中科院地理所遙感解譯數(shù)據(jù)分類(lèi)系統(tǒng)，將原州區(qū)的土地利用類(lèi)型分為10類(lèi)（表1），地類(lèi)名稱(chēng)以《土地利用現(xiàn)狀分類(lèi)》GB/T 21010-2017作為標(biāo)準(zhǔn)。

表1 原州區(qū)三生用地分類(lèi)體系與土地利用分類(lèi)體系的銜接Tab.1 Connection between classification of production-living-ecological space and classification of land use in Yuanzhou District

1.3.2 定量分析方法基于土地利用轉(zhuǎn)移流的溯源分析法。土地利用轉(zhuǎn)移流具有標(biāo)記流向和數(shù)量的屬性。本文提出了基于土地利用轉(zhuǎn)移流的溯源分析法，追溯土地利用變化的驅(qū)動(dòng)因子及其貢獻(xiàn)率。在土地利用變化中，設(shè)定T期的土地利用代碼為Ct，T+1 期的代碼為Ct+1，將T期和T+1 期的土地利用進(jìn)行疊加分析，通過(guò)代碼運(yùn)算生成土地轉(zhuǎn)移流代碼。如T期某圖斑為旱地，代碼取2，T期中該圖斑中出現(xiàn)2種用地類(lèi)型水澆地和園地，對(duì)應(yīng)代碼編為21和23，采用面積統(tǒng)計(jì)器計(jì)算轉(zhuǎn)移數(shù)量。土地利用轉(zhuǎn)移編碼（Cn）公式如下：

社會(huì)網(wǎng)絡(luò)分析法。采用Ucinet 6.186 進(jìn)行土地利用轉(zhuǎn)移網(wǎng)絡(luò)制圖與分析。對(duì)土地利用類(lèi)型采用代碼方法表示，對(duì)轉(zhuǎn)移方向標(biāo)注箭頭，并依據(jù)土地利用轉(zhuǎn)移流的數(shù)量比例對(duì)線型進(jìn)行量化示意，并標(biāo)記重要轉(zhuǎn)移流的數(shù)量。

核密度分析法。采用ArcGIS密度分析模塊，對(duì)土地利用轉(zhuǎn)移斑塊轉(zhuǎn)換生成點(diǎn)格式數(shù)據(jù)，以斑塊的面積作為value進(jìn)行核密度制圖，生成土地利用轉(zhuǎn)移率（單位為%）柵格圖。

2 結(jié)果與分析

2.1 原州區(qū)三生用地動(dòng)態(tài)變化分析

2.1.1 原州區(qū)三生用地結(jié)構(gòu)變化2000—2018 年，原州區(qū)的生產(chǎn)用地由12.25×104hm2減少為11.27×104hm2，生活用地由0.79×104hm2擴(kuò)張到1.72×104hm2；生態(tài)用地14.33×104hm2增為14.38×104hm2；生產(chǎn)、生活和生態(tài)用地結(jié)構(gòu)由2000 年的44.76%、2.90%、52.34%調(diào)整為2018 年的41.17%、6.29%、52.54%。這種變化與生態(tài)保護(hù)和修復(fù)政策、農(nóng)田水利建設(shè)以及城市化發(fā)展密切相關(guān)。1999 年國(guó)家開(kāi)始了退耕還林政策，原州區(qū)大于15°坡地漸次退耕還林；固海擴(kuò)灌揚(yáng)水工程也于1999年破土動(dòng)工。隨著水利建設(shè)工程的實(shí)施和揚(yáng)黃節(jié)水重大項(xiàng)目的持續(xù)推進(jìn)，清水河谷地的旱地大幅減少，水澆地顯著增加；隨著農(nóng)業(yè)產(chǎn)業(yè)機(jī)構(gòu)的調(diào)整，園地面積略有擴(kuò)張。研究時(shí)段內(nèi)，原州區(qū)的城市化也以較快的速度推進(jìn)，使得生活用地?cái)U(kuò)張迅速。生態(tài)用地總量呈弱增態(tài)勢(shì)，但林地面積增加較為迅速，其他草地大幅減少，生態(tài)用地結(jié)構(gòu)有所優(yōu)化、質(zhì)量有所提升（表2）。

表2 2000—2018年原州區(qū)三生用地面積變化Tab.2 Area change of the production-living-ecological space structure in Yuanzhou District during 2000—2018 /104 hm2

2.1.2 原州區(qū)三生用地轉(zhuǎn)移網(wǎng)絡(luò)結(jié)構(gòu)特征不同地類(lèi)之間產(chǎn)生了錯(cuò)綜復(fù)雜的轉(zhuǎn)移流，且轉(zhuǎn)移數(shù)量有顯著差異（圖2）。轉(zhuǎn)移流較大的地類(lèi)有旱地、林地、天然牧草地、其他草地，構(gòu)成了原州區(qū)土地利用轉(zhuǎn)移的關(guān)鍵地類(lèi)。這幾種地類(lèi)之間的轉(zhuǎn)移流形成的一級(jí)轉(zhuǎn)移流和二級(jí)轉(zhuǎn)移流是該轉(zhuǎn)移網(wǎng)絡(luò)的基本骨架，決定了原州區(qū)土地利用轉(zhuǎn)移的基本特征。

圖2 2000—2018年原州區(qū)三生用地利用轉(zhuǎn)移網(wǎng)絡(luò)Fig.2 Transfer flow net of production-living-ecological space in Yuanzhou District during 2000—2018

對(duì)占比大于1%的轉(zhuǎn)移流按坡度統(tǒng)計(jì)（圖3），發(fā)現(xiàn)轉(zhuǎn)向高功能生態(tài)用地的地域主要發(fā)生在海拔較高、坡度較大的不宜耕種和不宜居住生活的區(qū)域，體現(xiàn)“生態(tài)空間山清水秀”導(dǎo)向。其中，旱地向林地的轉(zhuǎn)移面積占總面積的7.73%，主要發(fā)生于坡度5°～15°地域；其他草地轉(zhuǎn)林地面積占總面積的8.25%，主要發(fā)生于坡度5°～15°地域；通過(guò)宜林地的植樹(shù)造林和荒山自然封育，天然牧草地向林地的轉(zhuǎn)移面積占總面積的5.40%，其他草地向天然牧草地的轉(zhuǎn)移面積占總面積的4.53%，主要發(fā)生在坡度5°～25°地域。對(duì)坡度適宜區(qū)域進(jìn)行坡改梯田、對(duì)平緩河谷地帶增加灌溉條件，體現(xiàn)“生產(chǎn)空間集約高效”導(dǎo)向。其他草地向旱地的轉(zhuǎn)移面積占總面積的8.06%，主要發(fā)生于坡度5°～15°地域，以梯田建設(shè)為主；旱地改水澆地面積占總面積的4.32%，主要發(fā)生在坡度0～5°地域。生活空間的擴(kuò)張發(fā)生在居住條件更為適宜的地域，體現(xiàn)“生態(tài)空間宜居”導(dǎo)向。鄉(xiāng)村居民點(diǎn)占用旱地面積占總面積的1.09%，城市擴(kuò)張及其他建設(shè)用地占用旱地面積占總面積的1.46%，主要發(fā)生在坡度0～5°地域。

圖3 2000—2018年原州區(qū)不同坡度三生用地主要轉(zhuǎn)移流Fig.3 Transfer flow net of production-living-ecological space in Yuanzhou District during 2000—2018

2.2 原州區(qū)三生空間格局演變特征分析

2.2.1 原州區(qū)三生空間格局變化2000年原州區(qū)的農(nóng)業(yè)生產(chǎn)用地以旱地為主，主要集中分布于中部清水河谷，其他散布于黃土丘陵溝壑區(qū)和六盤(pán)山陰濕區(qū)；生態(tài)用地以其他草地最為廣泛，主要分布于東部的黃土丘陵溝壑區(qū)和西部的六盤(pán)山陰濕區(qū)，天然牧草地主要分布于黃土丘陵溝壑區(qū)北部，林地?cái)?shù)量較少，主要分布于六盤(pán)山陰濕區(qū)。城鎮(zhèn)用地范圍不大，鄉(xiāng)村居民點(diǎn)除了六盤(pán)山陰濕區(qū)較少外，散布于研究區(qū)（圖4）。

圖4 2000年和2018年原州區(qū)三生用地空間格局變化Fig.4 Changes in spatial pattern of production-livingecological space in Yuanzhou District in 2000 and 2018

通過(guò)土地利用結(jié)構(gòu)的調(diào)整和土地整理工程的推進(jìn)，原州區(qū)的三生空間格局明顯優(yōu)化。根據(jù)寧發(fā)改西部［2010］129號(hào)和［2012］233號(hào)文件精神，一是大力實(shí)施鞏固退耕還林“基本口糧田”成果并發(fā)展水澆地項(xiàng)目。隨著固海擴(kuò)灌工程的實(shí)施，清水河谷平原區(qū)部分旱地改水澆地和園地；在水利基礎(chǔ)設(shè)施相對(duì)薄弱原州區(qū)東西部的冬至河、茹河等流域退耕農(nóng)田范圍內(nèi)發(fā)展節(jié)水灌溉農(nóng)業(yè)，以突出水資源高效利用，促進(jìn)農(nóng)業(yè)結(jié)構(gòu)調(diào)整，提高了農(nóng)業(yè)用地的生產(chǎn)集約性，實(shí)現(xiàn)退耕戶農(nóng)業(yè)增效和農(nóng)民增收。與之相應(yīng)，城鄉(xiāng)聚落向中部河谷區(qū)集中態(tài)勢(shì)。二是在中部的緩坡丘陵區(qū)和西南部的山區(qū)緩坡區(qū)大力開(kāi)展土地整理，大量坡耕地和其他草地轉(zhuǎn)為梯田旱地，既保障了居民的口糧田，又提高了水土保持能力。三是隨著退耕還林、荒山造林、自然封育等生態(tài)恢復(fù)措施，形成了環(huán)繞原州區(qū)中心城區(qū)的“U”型森林生態(tài)帶。2018年，城鄉(xiāng)聚落方面以清水河谷平原區(qū)較為密集，六盤(pán)山陰濕區(qū)和黃土丘陵溝壑區(qū)較為稀疏；農(nóng)業(yè)生產(chǎn)用地以清水河谷平原區(qū)為主；林地和草地主要分布于六盤(pán)山陰濕區(qū)和黃土丘陵溝壑區(qū)。

2.2.2 三生空間格局演變的集聚性特征原州區(qū)的三生用地空間格局重構(gòu)，在一定程度上體現(xiàn)了山水林田湖草生命共同體理論，作為“一盤(pán)棋”根據(jù)其適宜性進(jìn)行優(yōu)化配置，三生空間格局轉(zhuǎn)換呈現(xiàn)出明顯空間集聚特征（圖5）。

總體上，土地利用轉(zhuǎn)移流在六盤(pán)山陰濕區(qū)呈現(xiàn)為南北延伸的高值區(qū)，黃土丘陵溝壑區(qū)為散點(diǎn)式高值區(qū)，清水河谷為團(tuán)塊狀高值區(qū)（圖5a）。退耕還林和荒山造林使生態(tài)用地增加且質(zhì)量提升。其中，旱地轉(zhuǎn)林地流的高值區(qū)主要發(fā)生在六盤(pán)山陰濕區(qū)海拔1750～2000 m、坡度＞15°的坡耕地以及黃土丘陵區(qū)（圖5b）。六盤(pán)山陰濕區(qū)南北延伸帶和黃土丘陵溝壑區(qū)存在草地轉(zhuǎn)林地高值區(qū)（圖5c）。黃土丘陵溝壑區(qū)的“坡改梯”行動(dòng)和清水河谷區(qū)的“旱改水”行動(dòng)形成生產(chǎn)用地質(zhì)量提升集聚區(qū)。隨著基本口糧田建設(shè)工程推進(jìn)，對(duì)10°～15°的坡耕地以及其他宜農(nóng)荒地一并整理新建梯田，在炭山鄉(xiāng)、寨科鄉(xiāng)、頭營(yíng)鎮(zhèn)、官?gòu)d鎮(zhèn)、開(kāi)城鎮(zhèn)以及中河鄉(xiāng)、張易鎮(zhèn)等鄉(xiāng)鎮(zhèn)形成高值區(qū)（圖5d）。隨著固海擴(kuò)灌工程的實(shí)施，在清水河谷推行“旱改水”行動(dòng)，在黃鐸堡鎮(zhèn)、三營(yíng)鎮(zhèn)、頭營(yíng)鎮(zhèn)、彭堡鎮(zhèn)等鄉(xiāng)鎮(zhèn)形成高值區(qū)（圖5e）。原州區(qū)在輻射與帶動(dòng)寧夏南部區(qū)縣的發(fā)展中具有重要作用，城鎮(zhèn)人口比例由2000年的16.47%增長(zhǎng)為2018年的51.19%，城區(qū)擴(kuò)張迅速，形成以城區(qū)和重點(diǎn)鄉(xiāng)鎮(zhèn)為核心的高值區(qū)（圖5f）。

圖5 2000—2018年原州區(qū)主要轉(zhuǎn)移流空間集聚格局Fig.5 Spatial agglomeration pattern of main transfer flows in Yuanzhou District during 2000—2018

2.3 原州區(qū)國(guó)土空間格局演變的驅(qū)動(dòng)機(jī)制

原州區(qū)受經(jīng)濟(jì)發(fā)展水平低下和生態(tài)環(huán)境脆弱雙重制約。歷經(jīng)八七扶貧計(jì)劃、東西對(duì)口扶貧協(xié)作、精準(zhǔn)扶貧等多個(gè)階段，向縣外遷移人口7.8×104人，縣內(nèi)由山區(qū)向河谷區(qū)遷移3.67×104人。生態(tài)移民工程推行的主要目的是解決原州區(qū)人口超載問(wèn)題；生態(tài)恢復(fù)和重建則主要是針對(duì)原州區(qū)的生態(tài)脆弱問(wèn)題；扶貧開(kāi)發(fā)主要是建立生態(tài)經(jīng)濟(jì)體系從而為貧困區(qū)縣的發(fā)展賦能。與上述宏觀政策相呼應(yīng)，依據(jù)作用于三生空間的轉(zhuǎn)移流進(jìn)行驅(qū)動(dòng)力歸因，提出了生態(tài)移民與鄉(xiāng)村建設(shè)、城市化與工業(yè)化、水土資源集約利用、生態(tài)建設(shè)等驅(qū)動(dòng)因子。在上述驅(qū)動(dòng)因子的作用下產(chǎn)生了規(guī)模不等的土地利用轉(zhuǎn)移流，不同的轉(zhuǎn)移流作用于相應(yīng)的地域空間對(duì)象層，導(dǎo)致三生空間數(shù)量結(jié)構(gòu)和空間結(jié)構(gòu)變化，基于此，提出了由“動(dòng)力層-響應(yīng)層-結(jié)果層”構(gòu)成國(guó)土空間演變的動(dòng)力機(jī)制（圖6）。

圖6 原州區(qū)國(guó)土空間演變的動(dòng)力機(jī)制及重要轉(zhuǎn)移流占總面積的比例Fig.6 Dinamic mechanism of territorial space evolution and the proportion of important land use transfer flows to total area in Yuanzhou District

采用溯源分析方法，從三生數(shù)量變化和格局變化分析入手，依次梳理導(dǎo)致生活空間、生產(chǎn)空間和生態(tài)空間變化的土地利用轉(zhuǎn)移流，推演驅(qū)動(dòng)因子類(lèi)型，并依據(jù)土地利用轉(zhuǎn)移流所占比例來(lái)定量測(cè)度驅(qū)動(dòng)因子的貢獻(xiàn)率大小。2000—2018 年原州區(qū)的土地利用轉(zhuǎn)移流為15.96×104hm2，占原州區(qū)總面積的58.35%，據(jù)此推演的各類(lèi)驅(qū)動(dòng)因子貢獻(xiàn)率如表3 所示。生態(tài)建設(shè)是最主要的驅(qū)動(dòng)因子，對(duì)原州區(qū)三生空間的驅(qū)動(dòng)貢獻(xiàn)率高達(dá)49.40%。退耕還林、荒山造林、自然封育等生態(tài)建設(shè)行動(dòng)，植被覆蓋度明顯提升，涵養(yǎng)水源能力、生物多樣性保護(hù)功能顯著提升。為解決水土流失嚴(yán)重的問(wèn)題，先后開(kāi)展了茹河一期、馬什好、開(kāi)城、楊達(dá)子溝、冬至河等水土保持綜合治理工程，改善了山區(qū)坡耕地就地?cái)r蓄徑流和泥沙的能力。其次是水土資源集約節(jié)約利用，對(duì)三生空間格局調(diào)整的貢獻(xiàn)率達(dá)到31.01%。固海擴(kuò)灌工程以及近20 a川道旱地轉(zhuǎn)水澆地、坡耕地改寬幅梯田的建設(shè)，有效保持了水土，改善了農(nóng)業(yè)生產(chǎn)及生態(tài)條件，提升了水土資源集約利用水平。生態(tài)移民與鄉(xiāng)村建設(shè)、城市化和工業(yè)化的貢獻(xiàn)率都比較低。生態(tài)移民與鄉(xiāng)村建設(shè)對(duì)三生空間的貢獻(xiàn)率僅為2.94%，城市化與工業(yè)化對(duì)三生空間格局驅(qū)動(dòng)的貢獻(xiàn)率為3.63%。需要注意的是，水土流失仍然是各類(lèi)用地向未利用地轉(zhuǎn)化的重要?jiǎng)恿?，?3.02%，今后仍需加大對(duì)水土流失的治理。

表3 原州區(qū)國(guó)土空間演變驅(qū)動(dòng)力貢獻(xiàn)率Tab.3 Driving force contribution rate of territorial spatial evolution in Yuanzhou District

3 討論

（1）基于土地利用轉(zhuǎn)移流的國(guó)土空間演變溯源分析方法

土地利用轉(zhuǎn)移流是采用土地利用轉(zhuǎn)移矩陣計(jì)算轉(zhuǎn)移數(shù)量，通過(guò)土地利用轉(zhuǎn)移圖斑的源代碼和目標(biāo)代碼，從而實(shí)現(xiàn)對(duì)土地利用轉(zhuǎn)移的來(lái)龍去脈和數(shù)量進(jìn)行定量分析的一種方法。依靠土地利用轉(zhuǎn)移流的流向，采用演繹法可以定性分析驅(qū)動(dòng)該土地利用發(fā)生變化的原因，從而確定推動(dòng)土地利用變化的驅(qū)動(dòng)因子。以土地利用轉(zhuǎn)移流的分析分別向原因端和結(jié)果端推演，建立了“驅(qū)動(dòng)因子-土地利用轉(zhuǎn)移流-三生空間演變”的“動(dòng)力-響應(yīng)-結(jié)果”傳導(dǎo)性驅(qū)動(dòng)機(jī)制。溯源分析發(fā)現(xiàn)，生態(tài)移民與鄉(xiāng)村建設(shè)、城市化與工業(yè)化、水土資源集約利用、生態(tài)建設(shè)、水土流失是推動(dòng)原州區(qū)國(guó)土空間格局演變的驅(qū)動(dòng)力。對(duì)各驅(qū)動(dòng)因子涉及的土地利用轉(zhuǎn)移流的比例進(jìn)行百分制折算，獲得不同驅(qū)動(dòng)因子對(duì)三生空間格局變化的貢獻(xiàn)率。這種基于土地利用轉(zhuǎn)移流溯源分析的驅(qū)動(dòng)機(jī)制分析法，可為國(guó)土空間驅(qū)動(dòng)力研究提供新的視角。

（2）國(guó)土空間格局重構(gòu)的啟示與建議

原州區(qū)地處西北干旱地區(qū)，降水稀少，植被覆蓋率低，生態(tài)系統(tǒng)脆弱敏感，長(zhǎng)期以來(lái)存在“貧困-墾殖-生態(tài)惡化-貧困”的惡性循環(huán)。近20 a 來(lái)，原州區(qū)在國(guó)土空間格局重構(gòu)中，遵循“山水林田湖草”生命共同體理論，基于“人-地-業(yè)”系統(tǒng)“一盤(pán)棋”的調(diào)控思路，成效顯著。針對(duì)人口超載問(wèn)題，“減壓”和“擴(kuò)容”策略兼用，針對(duì)脆弱生態(tài)系統(tǒng)，主要采取了生態(tài)恢復(fù)和重建策略；針對(duì)貧困問(wèn)題，首先在人口問(wèn)題上采用移民搬遷和勞務(wù)輸出主要手段，為區(qū)域減貧作出積極貢獻(xiàn)；針對(duì)低下的生態(tài)經(jīng)濟(jì)體系和民生欠賬問(wèn)題，通過(guò)建立生態(tài)經(jīng)濟(jì)體系為區(qū)縣發(fā)展賦能。總體上，以移民搬遷作為破解人地矛盾的突破口，兼顧糧食生產(chǎn)和生態(tài)保護(hù)雙重需求，促進(jìn)國(guó)土空間向生產(chǎn)空間增效、生活空間擴(kuò)容、生態(tài)空間提質(zhì)方向演化。今后需要進(jìn)一步落實(shí)“兩山論”轉(zhuǎn)化需求，加快生態(tài)產(chǎn)品價(jià)值實(shí)現(xiàn)機(jī)制的創(chuàng)新，從而有效保障和提升三生空間格局的進(jìn)一步優(yōu)化和提質(zhì)增效。主要措施建議如下，一是重點(diǎn)探索如何建立綠色地理標(biāo)志產(chǎn)品的研發(fā)與規(guī)模化發(fā)展；二是重點(diǎn)探索如何提供優(yōu)質(zhì)生態(tài)系統(tǒng)調(diào)節(jié)與服務(wù)產(chǎn)品；三是要著力考慮如何建立有效的生態(tài)產(chǎn)品交易平臺(tái)與高質(zhì)量的物流服務(wù)體系。

4 結(jié)論

采用土地利用轉(zhuǎn)移流、社會(huì)網(wǎng)絡(luò)分析方法以及核密度分析法，探究了2000—2018年原州區(qū)國(guó)土空間格局演變及其動(dòng)力機(jī)制。結(jié)論如下：

（1）原州區(qū)三生空間向生產(chǎn)空間增效、生活空間擴(kuò)容、生態(tài)空間提質(zhì)的方向演化。2000—2018年生產(chǎn)、生活和生態(tài)用地結(jié)構(gòu)由2000 年的44.76%、2.90%、52.34%調(diào)整為2018 年的41.17%、6.29%、52.54%，呈現(xiàn)出生產(chǎn)用地減少、生活用地?cái)U(kuò)張比例較大、生態(tài)用地弱增的特征。

（2）原州區(qū)三生用地轉(zhuǎn)移趨勢(shì)與現(xiàn)行國(guó)土空間優(yōu)化目標(biāo)有較高的契合度。轉(zhuǎn)向生態(tài)用地的地域主要發(fā)生在不宜耕作或不宜居住區(qū)，體現(xiàn)出“生態(tài)空間山清水秀”導(dǎo)向；通過(guò)坡改梯田、提供灌溉條件以提升水土資源集約利用水平，體現(xiàn)出“生產(chǎn)空間集約高效”導(dǎo)向；生活空間的擴(kuò)張主要發(fā)生在居住條件更為適宜的河谷區(qū)，體現(xiàn)出“生態(tài)空間宜居”導(dǎo)向。

（3）三生空間演變的空間集聚特征突出。土地利用轉(zhuǎn)移流在六盤(pán)山陰濕區(qū)呈現(xiàn)為南北延伸帶狀高值區(qū)、黃土丘陵溝壑區(qū)為散點(diǎn)高值區(qū)，清水河谷鄉(xiāng)鎮(zhèn)和縣城為團(tuán)塊狀高值區(qū)。不同關(guān)鍵地類(lèi)的轉(zhuǎn)移流在空間上呈現(xiàn)出不同的空間集聚特征。

（4）基于土地利用轉(zhuǎn)移流溯源法，提出了原州區(qū)國(guó)土空間格局演變的“動(dòng)力-響應(yīng)-結(jié)果”傳導(dǎo)性動(dòng)力機(jī)制。生態(tài)建設(shè)是原州區(qū)國(guó)土空間格局演變的最重要的驅(qū)動(dòng)力，其次是水土資源集約利用，二者對(duì)原州區(qū)三生空間演變貢獻(xiàn)率分別為49.40%和31.01%；水土流失仍然是各類(lèi)用地向未利用地轉(zhuǎn)化的重要?jiǎng)恿?，?3.02%，今后仍需加大對(duì)水土流失的治理。