王 禎, 吳金華, 李嘉會, 魏金戈
(長安大學 土地工程學院, 陜西 西安 710054)
耕地資源是糧食生產(chǎn)的空間與物質(zhì)基礎(chǔ),是人類賴以生存的重要資源。然而隨著城鎮(zhèn)化進程的加快,耕地數(shù)量減少、耕地質(zhì)量下降、耕地污染等問題已經(jīng)威脅到耕地糧食生產(chǎn)能力[1]。2020年12月16日至18日舉行的中央經(jīng)濟工作會議將耕地問題作為2021年8項重點任務之一,并指出[2]:“保障糧食安全,關(guān)鍵在于落實藏糧于地,藏糧于技戰(zhàn)略。堅決遏制耕地‘非農(nóng)化’,防止‘非糧化’,規(guī)范耕地占補平衡”。這表明耕地與糧食問題越來越成為關(guān)乎國家發(fā)展和人民生活的關(guān)鍵問題[3]。
土地承載力是指一定地區(qū)的土地所能持續(xù)供養(yǎng)的人口數(shù)量,即土地人口承載量。其實質(zhì)是研究人口消費與食物生產(chǎn)、人類需求與資源供給之間的平衡發(fā)展問題[4]。就本質(zhì)而言,糧食由耕地生產(chǎn),因此也有部分學者[1,3,5-12]將土地承載力稱作耕地承載力。耕地承載力研究由來已久,國內(nèi)外學者對耕地承載力的研究主要集中在測算[3,5-12]與預測[3,5-9.13-17]。在測算方面,主要方法包括基于現(xiàn)實糧食產(chǎn)量[3,5-7,9,13-14,18-20]與基于生產(chǎn)潛力的農(nóng)業(yè)生態(tài)區(qū)域法[21-24]。在預測方面,常用的方法包括灰色系統(tǒng)預測模型[5-7,15]、平均增長法[3,13]、一元線性回歸模型[3,14]。目前的觀點認為,中國可以在不給全球可持續(xù)發(fā)展帶來壓力的情況下,在空間容量和糧食生產(chǎn)方面實現(xiàn)基本自給自足[25]。但中國的人糧關(guān)系在空間上存在異質(zhì)性,區(qū)域與區(qū)域的差異明顯。綜合來看,國內(nèi)基于現(xiàn)實糧食產(chǎn)量的耕地承載力研究主要以多時點的測算與預測為主,對區(qū)域內(nèi)耕地承載力的空間特性研究較少;同時現(xiàn)有研究對耕地承載力與耕地承載指數(shù)的影響因素的定量化分析較少。
鑒于此,本文以陜西省為研究區(qū),運用重心遷移模型、地理探測器、灰色系統(tǒng)預測模型等方法,分析陜西省2010—2019年耕地承載指數(shù)的時空變化特征,找出影響耕地資源承載指數(shù)的關(guān)鍵因素,按不同標準對陜西省2020—2025年的耕地承載指數(shù)與人糧關(guān)系進行預測,進而根據(jù)分析結(jié)果,對陜西省11個地級行政區(qū)提出差異化的人糧關(guān)系改善對策。研究結(jié)果可為城市發(fā)展方向、國土空間規(guī)劃、耕地保護政策制定提供參考。
陜西省位于我國西北內(nèi)陸,介于東經(jīng)105°29′—111°15′,北緯31°42′—39°35′,總面積2.06×105km2。地勢南北高,中間低,根據(jù)自然地理環(huán)境與社會文化差異,可將陜西省進一步劃分為陜北、關(guān)中、陜南3個區(qū)域。根據(jù)《陜西省統(tǒng)計年鑒》,2010年以來,全省常住人口增長迅速,其中西安市是人口增長的主要城市,2010—2019年共增長156.91萬人,總增長率為18.60%。2019年末陜西省常用耕地面積3.01×106hm2,與2010年相比增加1.50×105hm2,總增長率為5.24%。陜西省10 a間糧食產(chǎn)量增加6.62×105t,總增長率為5.68%?!吨腥A人民共和國國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和2035年遠景目標綱要》將我國的農(nóng)業(yè)生產(chǎn)劃分為糧食生產(chǎn)功能區(qū)、重要農(nóng)產(chǎn)品保護區(qū)、復種區(qū)3類。陜西省主要是糧食生產(chǎn)功能區(qū),集中分布在關(guān)中平原與榆林市西部,其他地區(qū)為零星分布[26]。
研究使用的陜西省及各地級行政區(qū)的耕地面積、糧食產(chǎn)量、人口、地區(qū)生產(chǎn)總值等數(shù)據(jù)均來源于各年份的《陜西省統(tǒng)計年鑒》。其中,耕地面積使用的是年末常用耕地面積;糧食產(chǎn)量包括谷物、薯類和豆類等,不包括林果、蔬菜與其他經(jīng)濟作物;人口數(shù)據(jù)使用的是常住人口。
2.2.1 耕地資源承載力與耕地資源承載指數(shù) 耕地資源承載力反映的是區(qū)域人口與糧食的關(guān)系,可以用一定人均糧食年消費水平下,區(qū)域糧食生產(chǎn)力所能供養(yǎng)的人口規(guī)模來度量,其計算表達式為:
(1)
式中:LCC表示耕地資源承載力(人);G表示糧食產(chǎn)量(kg);GPC為人均糧食年消費量標準(kg/人)。根據(jù)聯(lián)合國糧農(nóng)組織公布的糧食人均年消費標準可分為溫飽型(400 kg/人)、小康型(450 kg/人)和富裕型(500 kg/人)3級[13],其中400 kg/人可達到人體營養(yǎng)安全要求[18]。因此在計算2011—2019年陜西省耕地資源承載力時采用400 kg/人的標準。
耕地資源承載指數(shù)模型由封志明等[18]構(gòu)建,即現(xiàn)實人口與基于人糧關(guān)系的耕地可承載人口的比值,可以揭示耕地承載力與現(xiàn)實人口的關(guān)系,計算表達式為:
(2)
式中:LCCI為耕地資源承載指數(shù);Pa表示現(xiàn)實人口數(shù)量(人)。
2.2.2 重心遷移模型 運用重心遷移模型來表征陜西省耕地資源承載指數(shù)重心的變化及遷移路徑,以分析耕地資源承載指數(shù)的空間變化規(guī)律。遷移路徑的方向與長度能在一定程度上反映耕地資源承載指數(shù)增加的空間特征與強度,其計算表達式為:
(3)
式中:X,Y分別代表某年某地級行政區(qū)的重心的經(jīng)緯度坐標;n代表地級行政區(qū)數(shù)量; LCCIi代表某年某地級行政區(qū)的耕地資源承載指數(shù);Xi,Yi代表某地級市的幾何中心地理坐標。
2.2.3 分異及因子探測器 運用地理探測器中的因子探測器研究陜西省耕地資源承載指數(shù)時空變化的影響因素。地理探測器是探測空間分異性、揭示其背后驅(qū)動力的一組統(tǒng)計學方法。其核心思想是假設(shè):如果某個自變量對某個因變量有重要影響,那么自變量和因變量的空間分布應該具有相似性[27-28]。分異及因子探測器是其中一種方法,可以探測某因子X多大程度上解釋了屬性Y的空間分異。用q值度量,計算表達式為[29]:
(4)
2.2.4 灰色預測模型 灰色系統(tǒng)由鄧聚龍[30]提出,是指在所研究的對象或數(shù)據(jù)并不完備或確切,具有灰色特性的系統(tǒng)。一般的社會經(jīng)濟系統(tǒng)都是灰色系統(tǒng)。采用灰色系統(tǒng)預測模型GM(1, 1)分別對陜西省各地級行政區(qū)的人口與糧食產(chǎn)量進行預測,計算出各地級行政區(qū)每年的耕地資源承載指數(shù),作為人糧關(guān)系預測的依據(jù)。其基本步驟是先對原始數(shù)據(jù)X0進行累加,使其變?yōu)橛幸?guī)律的生成數(shù)列Xt,用微分方程求解得出預測模型函數(shù)[7]:
(5)
式中:t= 1,2,…,n;u,a為 GM(1, 1)模型的參數(shù)估計值,采用最小二乘法進行估計。構(gòu)建GM(1,1) 模型前需要對原始數(shù)據(jù)進行級比值檢驗,當所有原始數(shù)據(jù)均落于規(guī)定的級比值區(qū)間中,才適合構(gòu)建GM(1,1)模型,否則需要進行平移轉(zhuǎn)換。GM(1, 1)模型要求數(shù)據(jù)序列變化平緩、符合指數(shù)函數(shù)性質(zhì)[31],因此構(gòu)建模型后還需要進行后驗差比檢驗與模型殘差檢驗,一般要求后驗差比小于0.65,相對誤差小于0.20,可達到精度要求。
2.2.5 平均增長法 采用平均增長法對于無法通過后驗差比檢驗與殘差檢驗的GM(1,1)模型進行預測。平均增長法是一種指數(shù)型增長,在人口預測[32]、耕地面積預測[33]中都有應用,其計算公式為:
P=N×(1+R)n
(6)
式中:P為預測值;N為預測初始年份的實際值;R為平均增長率;n為預測年數(shù)。
3.1.1 人糧關(guān)系劃分標準 根據(jù)封志明等[18]的研究,結(jié)合陜西省實際,從糧食供需維度,將陜西省及其各地級行政區(qū)每年的人糧關(guān)系分為供給過剩、供大于需、供需平衡、需大于供、需求過剩5類。各類人糧關(guān)系對應的耕地資源承載指數(shù)詳見表1。
表1 各類人糧關(guān)系對應的耕地資源承載指數(shù)
3.1.2 陜西省耕地承載指數(shù)總體變化 總體來看,10 a間陜西省的糧食產(chǎn)量、人口、耕地承載力都有所上升,耕地承載指數(shù)下降。耕地承載指數(shù)在1.23~1.28之間波動,且波動幅度不大。原因是雖然陜西省的人口持續(xù)增長,但糧食產(chǎn)量一直處于波動的狀態(tài),導致了耕地承載力與耕地資源承載指數(shù)的波動。從人糧關(guān)系的角度,10 a間陜西省人糧關(guān)系始終為需大于供,人糧關(guān)系亟待改善。2010—2019年陜西省耕地承載指數(shù)變化詳見表2。
表2 陜西省2010-2019年耕地承載指數(shù)
3.1.3 陜西省耕地承載指數(shù)時空變化 2010—2019年陜西省各地級行政區(qū)耕地承載指數(shù)變化如圖1所示。從變化方向上來看,根據(jù)各地級行政區(qū)耕地承載指數(shù)的時間變化,可將地級行政區(qū)分為上升型與下降型兩類。下降型城市包括榆林市、銅川市,其余城市均為上升型。從變化趨勢上來看,陜西省所有地級行政區(qū)的耕地承載指數(shù)都具有一定的波動性,變化較為劇烈的時間段為2010—2011年、2014—2015年、2017—2018年。2010年—2011年,各地級行政區(qū)糧食產(chǎn)量大幅下降,導致耕地承載指數(shù)上升。2014—2015年,楊凌示范區(qū)、商洛市、榆林市、延安市的耕地承載指數(shù)發(fā)生了劇烈的變化,其中楊凌示范區(qū)、商洛市為下降,榆林市、延安市為上升。原因是該時間段4個地級行政區(qū)的人口數(shù)量均上升,楊凌示范區(qū)、商洛市糧食產(chǎn)量上升,榆林市、延安市糧食產(chǎn)量下降。2017—2018年,楊凌示范區(qū)、西安市、安康市、商洛市、延安市的耕地承載指數(shù)明顯上升,咸陽市、商洛市、榆林市發(fā)生明顯下降。從分區(qū)角度,陜西省各分區(qū)及其內(nèi)部的耕地承載指數(shù)都具有一定的差異性。2019年陜北、關(guān)中、陜南的總體耕地承載指數(shù)分別為0.68,1.42和1.48,大小關(guān)系表現(xiàn)為:陜南地區(qū)>關(guān)中地區(qū)>陜北地區(qū)。
3.1.4 陜西省人糧關(guān)系時空變化 2010—2019年陜西省各地級行政區(qū)人糧關(guān)系變化如圖2所示。陜西省耕地承載指數(shù)總體上呈現(xiàn)出的南高北低的基本格局,2014年最明顯,從南向北依次為需大于供、需求過剩、供需平衡、供大于需。除咸陽市外,其余地級行政區(qū)的人糧關(guān)系變化都具有一定的方向性,咸陽市的人糧關(guān)系經(jīng)歷了供需平衡→需大于供→供需平衡的變化,體現(xiàn)了耕地承載指數(shù)時間變化上的波動性。從2010年開始,西安市與楊凌示范區(qū)就是陜西省的高值中心,且一直處于需求過剩的狀態(tài);榆林市是低值中心,一直處于供大于需的狀態(tài);到2019年,榆林市的耕地承載指數(shù)為0.51,逼近供給過剩的人糧關(guān)系。
3.1.5 陜西省耕地資源承載指數(shù)重心遷移 陜西省耕地資源承載指數(shù)重心遷移狀況如圖3所示。2010—2019年陜西省的耕地資源承載指數(shù)重心始終位于咸陽市境內(nèi),咸陽市與西安市交界處以北。10 a間重心的總體遷移方向為從東北向西南,指向西安市與楊凌示范區(qū),表明相對于其他方位,西南方向的地級行政區(qū)耕地資源承載指數(shù)增加得更劇烈。重心逐年遷移路徑中,2013—2015年、2017—2018年兩個時間段的路徑存在明顯的突變。2013—2014年,商洛市的耕地承載指數(shù)由1.44上升至1.80,上升比例為24.88%,為該時間段所有地級行政區(qū)最高,導致重心向其所在的東南方向遷移。2014—2015年,商洛市的耕地承載力回歸到之前的水平,而這一時間段西安市、楊凌示范區(qū)減少,延安市、榆林市增加,導致重心向北遷移。2017—2018年,耕地承載指數(shù)大幅上升的楊凌示范區(qū)、西安市、安康市、商洛市位于陜西省南部,而大幅下降的咸陽市、銅川市、榆林市位于陜西省偏北位置,導致該時間段重心長距離南移。
圖1 陜西省2010—2019年耕地承載指數(shù)變化
圖2 陜西省人糧關(guān)系時空變化
根據(jù)孫燕君等[1]的總結(jié),結(jié)合數(shù)據(jù)的可獲取性,選取人均耕地面積、人均地區(qū)生產(chǎn)總值(人均GDP)、糧食復種指數(shù)、地均糧食產(chǎn)量、第一產(chǎn)業(yè)產(chǎn)值占地區(qū)生產(chǎn)總值的比重5個因素作為分異及因子探測器的分析影響因子。將影響因子采用自然斷點法分為3層,根據(jù)公式計算每年各探測因子的q值,結(jié)果詳見表3。
表3 陜西省2010—2019年耕地承載指數(shù)分異及因子探測器計算結(jié)果
圖3 陜西省2010—2019年耕地資源承載指數(shù)重心遷移路徑
5個影響因素中,最大值出現(xiàn)在人均耕地面積,為0.775 6;最小值出現(xiàn)在人均GDP,為0.001 2。各影響因素的平均值大小排序為:人均耕地面積>第一產(chǎn)業(yè)產(chǎn)值占地區(qū)生產(chǎn)總值的比重>糧食復種指數(shù)>地均糧食產(chǎn)量>糧食復種指數(shù)。除人均耕地面積和糧食復種指數(shù)外,其他因素均存在方差較大,即波動較大的問題。q值時間變化情況如圖4所示。人均耕地面積的q值一直穩(wěn)定在較高水平的0.7左右,波動性較小;而地均糧食產(chǎn)量和第一產(chǎn)業(yè)占比分別在2010—2011年、2015—2016年兩個時間段發(fā)生突變;2010—2017年糧食復種指數(shù)的q值穩(wěn)定在0.2左右,但在2017—2018年時間段內(nèi)發(fā)生突變;人均GDP的q值經(jīng)歷過兩次突變,波動性較大。綜上所述,除人均耕地面積外,其他因素或多或少都存在波動問題,且q值較小,認為它們對陜西省耕地資源承載指數(shù)的影響較小或無影響;同時可以認為人均耕地面積是陜西省耕地資源承載指數(shù)的主要影響因素,平均可以解釋69.21%的耕地資源承載指數(shù)時空變化。
圖4 陜西省2010—2019年耕地資源承載指數(shù)q值的變化
運用GM(1, 1)模型分別對各地級行政區(qū)人口與糧食產(chǎn)量進行擬合與預測。經(jīng)檢驗(表4),有7個模型的級比值超過規(guī)定級比值區(qū)間。對7個模型進行平移轉(zhuǎn)換后,銅川市人口、延安市產(chǎn)量2個模型無法通過后驗差比檢驗;楊凌區(qū)產(chǎn)量、榆林市產(chǎn)量2個模型無法通過模型殘差檢驗。對于這4個模型,采用平均增長法進行預測。經(jīng)計算,銅川市人口、延安市產(chǎn)量、楊凌區(qū)產(chǎn)量、榆林市產(chǎn)量的平均增長預測模型分別為:P=N×0.996 1n,P=N×0.987 2n,P=N×1.070 4n,P=N×0.882 0n。
表4 陜西省2010—2019年耕地資源承載指數(shù)GM(1,1)模型擬合精度檢驗
根據(jù)灰色預測模型GM(1,1)的特性,預測2020—2025年的耕地承載指數(shù)及人糧關(guān)系。將各地級行政區(qū)的人口和產(chǎn)量累加得到全省數(shù)據(jù),根據(jù)公式分別按溫飽型(400 kg/人)、小康型(450 kg/人)和富裕型(500 kg/人)3種不同的人均糧食消費標準計算耕地資源承載指數(shù),計算結(jié)果詳見表5。根據(jù)表5,陜西省2020—2025年的耕地資源承載指數(shù)呈下降趨勢。按人均糧食消費水平的最低標準400 kg/人,陜西省2020—2025年的耕地承載指數(shù)將逐漸下降,到2025年時可達1.14,接近人糧平衡的1.125,人糧關(guān)系趨于緩和。但人均糧食消費水平是動態(tài)變化的,會隨著社會經(jīng)濟發(fā)展而改變。當使用450與500 kg/人標準進行預測時,耕地承載指數(shù)仍保持在一個相對較高的值。陜西省人糧關(guān)系時空變化預測結(jié)果如圖5所示。無論以何種人均糧食消費標準,陜西省2020—2025年的人糧關(guān)系都趨于逐步緩和的狀態(tài)。但細化到地級行政區(qū)尺度,仍存在空間差異明顯、人口超載區(qū)眾多等問題。當以400 kg/人作為人均糧食消費標準時,時空變化不顯著;當以450 kg/人作為人均糧食消費標準時,2025年榆林市成為供給過剩區(qū),咸陽市成為需大于供區(qū),安康市、延安市成為需求過剩區(qū)。當以500 kg/人作為人均糧食消費標準時,2020年起將不存在供需平衡區(qū),全省除榆林市外均為預測人口大于理論承載人口的區(qū)域,人口、糧食、耕地資源之間的矛盾十分尖銳。
表5 陜西省耕地資源承載指數(shù)預測結(jié)果
圖5 陜西省人糧關(guān)系時空變化預測結(jié)果
從分地級行政區(qū)層面,2010—2019年陜西省耕地承載指數(shù)的時空變化是人口、糧食產(chǎn)量、人均耕地面積相互作用的結(jié)果。在城市化快速發(fā)展的地區(qū),人口增長和建設(shè)用地擴張是主要特征,耕地保護面臨的主要問題是建設(shè)占用,導致糧食產(chǎn)量降低,如西安市、楊凌示范區(qū)等。而在黃土高原、秦嶺等需要生態(tài)保護的地區(qū),耕地保護與生態(tài)保護會存在一定沖突,若人口增長,人糧關(guān)系就會惡化,如安康市、延安市等。一些發(fā)展農(nóng)業(yè)的地區(qū),耕地保護與糧食生產(chǎn)面臨的主要困境是耕地的非糧化,這種現(xiàn)象在陜西省的水果產(chǎn)區(qū)尤為明顯,如關(guān)中平原、延安市山地、秦嶺北麓與漢江流域等。銅川市是傳統(tǒng)的資源型城市,近年來人口的快速流失是耕地承載指數(shù)降低的主要原因。由于生態(tài)治理成效顯著,加上耕地保護與復墾,使得2010—2019年榆林市耕地面積上升了37.83%,又因為人口增長不快,使得榆林市的人糧關(guān)系快速好轉(zhuǎn)。陜西省各地級行政區(qū)差異化的人糧關(guān)系改善對策詳見表6。
表6 陜西省人糧關(guān)系改善對策
(1) 2010—2019年陜西省總體糧食產(chǎn)量、人口、耕地承載力都有所上升,耕地資源承載指數(shù)小幅度波動下降,人糧關(guān)系稍微緩和;各地級行政區(qū)耕地資源承載指數(shù)波動較大,變化劇烈的時間段為2010—2011年、2014—2015年、2017—2018年,與糧食產(chǎn)量的劇烈變化有關(guān)。
(2) 2010—2019年陜西省耕地承載指數(shù)區(qū)域差異明顯,呈現(xiàn)出南高北低的空間格局,西安市、楊凌示范區(qū)為高值中心,榆林市為低值中心;從分區(qū)角度,耕地資源承載指數(shù)排序為:陜南地區(qū)>關(guān)中地區(qū)>陜北地區(qū);10 a間陜西省的耕地資源承載指數(shù)重心位于咸陽市境內(nèi),總體遷移方向為從東北向西南,指向西安市與楊凌示范區(qū),2013—2015年、2017—2018年兩個時間段的路徑存在明顯的突變,與各地級行政區(qū)耕地承載指數(shù)的消長有關(guān)。
(3) 陜西省耕地資源承載指數(shù)時空變化最主要的影響因素是人均耕地面積,平均解釋率為69.21%。由于糧食作物實際種植面積數(shù)據(jù)的缺失,年末常用耕地面積與基于人糧關(guān)系的耕地承載力指數(shù)并不能完全對應,可能存在耕地建設(shè)占用、從事林果業(yè)等情況。但地理探測器結(jié)果表明,保護耕地面積依舊是保證區(qū)域糧食安全的重要手段。
(4) 運用GM(1,1)、平均增長模型分不同人均糧食消費水平預測陜西省2020—2025年耕地承載指數(shù)與人糧關(guān)系的時空變化。結(jié)果顯示陜西省資源承載指數(shù)會小幅下降,空間格局仍為南高北低,除銅川市—渭南市、榆林市2個低值中心外,其余地級行政區(qū)的人糧關(guān)系較為緊張。
(5) 根據(jù)最新《陜西省統(tǒng)計年鑒》,2020年陜西省常住人口為3 955萬人,糧食產(chǎn)量為1.28×107t。分別以400 kg/人,450 kg/人和500 kg/人為糧食消費標準時,耕地資源承載指數(shù)為1.18,1.33,1.57,與本文的預測值分別相差4.97%,4.79%,5.30%。因此,本研究的預測結(jié)果具有一定參考價值。
(6) 在保障糧食安全的基礎(chǔ)上,《中華人民共和國國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和2035年遠景目標綱要》指出,要優(yōu)化農(nóng)業(yè)生產(chǎn)布局,積極發(fā)展設(shè)施農(nóng)業(yè),因地制宜發(fā)展林果業(yè)[26]。陜西省《全省國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和2035年遠景目標綱要》針對特色優(yōu)勢產(chǎn)業(yè)提出了更加具體的果業(yè)、畜牧業(yè)、設(shè)施農(nóng)業(yè)“3+X”現(xiàn)代農(nóng)業(yè)體系[34]。綜合來看,糧食安全與耕地保護是關(guān)乎國家安全的問題,因此國家與陜西省都將其放在了最重要的位置。目前,耕地保護受到的最主要的威脅就是耕地的“非農(nóng)化”與“非糧化”問題,在陜西省主要表現(xiàn)為占用耕地建房與改變耕地用途發(fā)展林果。要解決耕地保護與特色農(nóng)業(yè)發(fā)展、城鎮(zhèn)化建設(shè)之間的矛盾,就必須明確土地用途,從“防”與“治”兩方面入手?!胺馈钡闹饕侄伟▌澏ㄓ谰没巨r(nóng)田保護紅線,從規(guī)劃與立法層面嚴格限制耕地的用途轉(zhuǎn)變,明確與發(fā)展特色農(nóng)業(yè)發(fā)展和城鎮(zhèn)化建設(shè)的土地相區(qū)分;“治”即通過遙感和人力巡查等方式監(jiān)測耕地用途的改變,并強制整改或處罰。