摘要以西藏地區耕地面積為研究對象，結合主要農作物播種面積、產量等時空變化，從經濟、人口、農業發展水平等13個驅動因子，分析了西藏地區耕地面積變化。采用1985~2016年的統計資料，采用統計方法，分析了西藏地區30多年來耕地、糧食生產的動態變化，并從經濟、人口、農業發展水平等方面選取影響耕地面積變化的13個驅動因子進行分析。西藏地區耕地面積總體呈增加態勢，但增長趨勢緩慢;糧食作物單產、總產量在進一步增加，人均耕地面積在逐年減少;糧食主產區糧食總產不穩定，且所占比例有逐年減少趨勢;小麥、油菜播種面積正在逐年增加，青稞播種面積變化最小，但有逐年減少趨勢;全區土地利用程度相對較高;經濟發展、人口增長、農業發展水平是影響西藏耕地面積變化的主要驅動因子。西藏地區土地利用程度相對較高，基礎建設用地占補平衡，農牧民名生活結構正在逐漸發生改變。

　　關鍵詞耕地面積;動態變化;糧食總產量;驅動力;西藏

　　西藏位于中國的西南邊陲，面積120多萬km2，但絕大部分是高山荒漠、雪嶺冰川。西藏是青藏高原的主體部位，素有“世界屋脊”和“第三極”著稱于世。在全球氣候系統普遍變暖的大背景下，青藏高原氣溫也逐年升高，且降雨量在逐漸減少[1]，這對于本來高寒、缺氧、干旱等氣候環境惡劣的青藏高原來講，無疑對農業生產帶來了更大的挑戰。

　　農業方向論文投稿刊物：《農業工程學報》是由中國科學技術協會主管，中國農業工程學會主辦的全國性學術期刊。國際刊號ISSN：1002-6819;國內刊號CN：11-2047/S，郵發代號：18-57。讀者對象為農業工程學科及相關領域的科研、教學及生產科技人員、技術管理及推廣人員和大院校師生。

　　李萍[2]認為西藏糧食作物單產、鄉村勞動力人數和糧食作物播種面積是影響西藏糧食總產的最重要因素，但楊春艷等[2]研究發現，西藏糧食產量低而不穩定，且對氣候變化敏感，王建林等[4]研究發現，1981~2001年期間，西藏人均耕地面積減少了0.0301hm2，且年均遞減率為1.13%，郭志儀等[5]對西藏的人口進行預測，認為西藏人口將保持持續增長的態勢，且人口年齡結構從成年型向老年型轉變，劉合滿等[6]研究認為，西藏人口增長與糧食作物播種面積呈極顯著負相關關系，人均總農作物和糧食作物播種面積低于全國水平，且呈現不斷下降趨勢。人口的快速增加和糧食作物播種面積的不斷下降，將會對西藏糧食安全產生威脅，普布次仁[7]研究認為，以現有播種面積和自給生產為基礎，要達溫飽、小康和富裕生活，西藏平均大約還需提高單產100kg/hm2、200kg/hm2和300kg/hm2。

　　自改革開放以來，國家持續加大對西藏基礎設施建設的投入力度，特別黨的十八大以來，西藏經濟社會快速發展，人口的不斷增加，對土地的需求也將不斷增加，且西藏的后備耕地資源有限，耕地利用和保護面臨著嚴峻的挑戰。治國必治邊、治邊先穩藏。農業是政治家、消費者和環境保護者共同關注的一個焦點。本文以《西藏統計年鑒》為基礎數據，分析了西藏30多年來耕地面積、糧食作物單產、主要糧食作物播種面積、人均耕地面積等變化動態及驅動力分析，旨在為西藏地區農業生產發展提供一定的基礎資料。

　　1研究方法

　　本文以《西藏統計年鑒》為基礎資料，用Excel2013進行數據整理，用SPSS22.0進行統計分析。

　　2西藏地區糧食生產情況

　　可以看出，近30多年來，西藏播種面積一直處于增長趨勢，呈現“快-慢-快”趨勢，表明每年都會有新的土地資源被開發利用;其糧食總產量基本都處于增長趨勢，且1997~2001年處于直線增長趨勢，但播種面積處于緩慢增長動態，說明單產有很大提升。在2002~2010年期間，播種面積處于緩慢增長階段，糧食總產量在逐年的降低，結合周天財等[1]的研究發現，之所以導致總產量降低，是因為在這期間，最高氣溫、年平均溫度在逐年上升，最低氣溫在逐年降低、年平均降雨量在逐年減少所致。在2011年之后，總產量又逐漸提高，分析其原因，除播種面積逐年增加外，年平均氣溫逐年在降低、年降雨量在逐年增多。

　　2.1近30年西藏地區耕地面積及人均耕地面積動態變化實有耕地面積及人均耕地面積動態變化。可以看出，年末實有耕地面積變化基本分為兩個階段，增長階段和減少階段。1987~1999和2004~2016兩個階段，年末實有耕地面積基本處于增長階段，1985~1987和1999~2014階段，處于減少階段，這可能與政府調整結構有關，基礎設施建設等非農業用地在進一步擴大，生態退耕還草還林面積增大所致。由于西藏地區人口還在進一步增加，但耕地面積增長緩慢，人均耕地一直處于下降趨勢。

　　2.2各市區在全區糧食生產中的地位各市區農作物產量及播種面積。可以看出，日喀則、山南、昌都、拉薩是全區的糧食集中產區，阿里、那曲和林芝所占比例非常小。日喀則糧食總產量基本處于“快﹣慢﹣快”增長狀態。日喀則地區糧食作物播種面積一直處于穩步增長狀態，但在2005~2011年之間日喀則地區出現糧食減產現象，主要原因可能是由于在這期間氣溫升高，年降雨量減少有致;拉薩市、昌都、山南糧食總產在近10多年總體變化趨勢不大;那曲、阿里和林芝雖然播種面積、糧食總產所占比重小，其播種面積一直處于增長趨勢，其中那曲地區糧食總產量增幅較大，但阿里糧食總產量較2010年相比，一直處于減少狀態。在近10多年來，全區糧食主產區拉薩市、日喀則、山南糧食總產在全區糧食總產中所占比例不穩定，且基本呈逐年在減少狀態，其中以山南地區變化幅度相對較大，日喀則相對較小，昌都、那曲地區所占比重在逐漸增大。

　　2.3近30多年西藏主要農作物單產、種植面積動態變化主要農作物單位面積產量及播種面積。可以看出，近30多年來，主要農作物單產一直處于穩步增長態勢，其中小麥、冬小麥、青稞2016年單產較1985年增產均超過100%，其中以冬小麥增產幅度最大，達111.15%，豆類作物的增產幅度最小，僅為58.13%，同時表現出每年的單產不穩定。

　　就主要農作物播種面積來講，小麥和油菜的播種面積基本處于增長狀態，但同時表現出每年的播種面積不穩定，變幅較大。青稞每年的播種面積基本趨于穩定，變幅最小，豆類作物的播種面積最為不穩定。總之，小麥和油菜的播種面積有逐年增加的趨勢，青稞的播種面積最為穩定，但呈逐年減少趨勢，分析其原因，可能是人們的生活結構正在逐漸的改變。

　　2.4耕地動態變化分析

　　土地資源數量變化可以用土地利用動態度來表示，而其能夠定量地描述土地利用的變化速率，對于預測未來土地利用變化趨勢有積極作用[8]。參照王秀蘭等[8]方法計算土地利用動態：-1K=100%baaUUUT式中K為研究時段內某種土地利用類型的年變化率，Ua和Ub分別為某種土地利用類型研究初期與研究末期的面積，T為研究時段長。從《西藏統計年鑒》獲得數據如下，全區1985～2000年耕地面積減少了7.30khm2，2000～2016年耕地面積減少了7.28khm2，32年間耕地面積年變化率為0.2%。

　　1985～2000年，昌都、日喀則、林芝三地耕地面積呈減少態勢，其變化率分別為-2.68khm2、-6.61khm2、-2.41khm2，拉薩、山南、那曲、阿里耕地面積呈增加態勢;2000～2016年，山南、日喀則、林芝耕地面積呈減少態勢，其變化率分別為-1.69khm2、-6.61khm2、-0.64khm2，拉薩、昌都、阿里、那曲耕地面積呈增加態勢。1985～2016年期間，全區糧食主產區昌都、山南、日喀則土地利用狀態相對較高，拉薩土地利用動態度呈下降趨勢，這主要可能是近年來拉薩地區的基礎設施建設有關。

　　2.5耕地資源驅動因子分析

　　2.5.1耕地面積變化的影響因子

　　影響耕地面積變化的因素主要有自然和社會經濟因素兩方面，且社會經濟因素更重要，本文重點針對社會經濟因素進行分析[9]。將實有耕地面積計Y，依據統計資料選取總人口(X1)、城鎮人口(X2)、鄉村人口(X3)地區生產總(X4)第一產業總產值(X5)第二產業總產值(X6)第三產業總產值(X7)全社會固定資產投資(X8)農林牧漁業總產值(X9)糧食總產量(X10)城鎮本年完成投資(X11)人均耕地面積(X12)、單位面積產量(X13)等13個因子，并進因子分析，采用因子分析方法中的的最大方差旋轉法得出相關矩陣(表4)，特征值及主成分貢獻率(表5)、主成分載荷矩陣。

　　12個驅動因子間，存在不同程度相關性，總的來看，X12與各驅動因子間成負相關，X10與X6、X7、X8、X9，X10與X11相關性較小外，其他因子間相關性較大，信息重疊，主成分分析結果表明，第一主成分、第二主成分的特征值分別是6.724和6.118，其貢獻率分別為51.721%和47.062%，前兩個主成分累計貢獻率達98.783%，說明已經包含原有13因子的絕大部分信息，符合分析需求。

　　人均耕地面積與耕地面積呈負相關，其余各因子均與耕地面積呈正相關，第一主成分中，城鎮本年完成投資(X11)、全社會固定資產投資(X8)、第二產業總產值(X6)、地區生產總值(X4)和第三產業總產值(X7)相對載荷較高，可以解釋為經濟發展水平;第二主成分中糧食總產量(X10)、單位面積產量(X13)、鄉村人口(X3)、人均耕地面積(X12)、總人口(X1)和城鎮人口(X2)相對載荷較高，可以解釋為農業發展水平和人口。

　　第1和第2主成分已經基本保留了所有指標的原有信息，可以用2個變量Z1和Z2代替原來13個指標(總人口(X1)、城鎮人口(X2)、鄉村人口(X3)地區生產總(X4)第一產業總產值(X5)第二產業總產值(X6)第三產業總產值(X7)全社會固定資產投資(X8)農林牧漁業總產值(X9)糧食總產量(X10)城鎮本年完成投資(X11)人均耕地面積(X12)、單位面積產量(X13))，得出線性組合，其中X1～X13均為標準化后的變量，分別為：Z1=0.224X1+0.261X2+0.183X3+0.342X4+0.285X5+0.354X6+0.341X7+0.363X8+0.294X9+0.098X10+0.364X11﹣0.175X12+0.130X13Z2=0.327X1+0.292X2+0.350X3+0.185X4+0.271X5+0.159X6+0.188X7+0.133X8+0.260X9+0.380X10+0.122X11﹣0.357X12+0.380X13。

　　(1)經濟發展。

　　城鎮本年完成投資由1985年的58408萬元增加到16555043萬元，增加了283倍，全社會固定資產投資由1985年的74940萬元增加到16555043萬元，增加了220倍，第二產業總產值較1985年相比增加了140倍，第三產業總產值增加了105倍，地區生產總值增加了64倍。隨著社會經濟的快速發展，基礎設施的規模逐步擴大，人民生活水平日益提高，不可避免的引起非農業建設對土地的需求，雖然基礎設施投資在不斷擴大，但是西藏地域廣闊，占用農業用地相對較少，同時，西藏按照“占多少補多少”的原則，建設占用耕地時爭取做到“占一補二”或“占一補三”的效果，實現了土地利用“動態平衡”，起到了積極作用，實現了占補平衡，且補大于占。同時，經濟實力的不斷壯大，致使大量的荒地被開墾利用。

　　(2)人口增加。

　　人口的增加是一種持續的外界壓力，對耕地的數量變化起著雙重作用，一方面人口增長，會對生產與生活設施需求增多，大量的土地就會被占用，使人類進一步開墾利用荒地，提高土地利用程度提高，其次，人口增長人均土地將進一步減少[9]。鄉村人口從1985年的178.95萬人，增加到232.83萬人，增加了30.11%，總人口增加65.70%，城鎮人口增加了375.94%，由于人口的快速增長，需要進一步擴大耕地面積來滿足人們的日常生活需求，較1985年相比，年末實有耕地面積增加了6.52%。人均耕地面積由1985年0.107hm2/人，下降到0.072hm2/人，較1985年下降了35.7%，雖然耕地面積有所增加，但增長緩慢，人口增長較快，導致人均耕地面積不斷下降。

　　(3)農業發展水平。

　　建設用地、生態退耕、災毀耕地是西藏地區耕地減少的直接原因。發展科技，提高單產是確保西藏地區糧食安全的重要保障。西藏地區2016年糧食總產量為10.23648萬t，糧食單產5596kg/hm2，分別較1985年增長了48.16%和104.37%，而耕地面積2016年較1985年增加了6.52%，所以，糧食總產量的增加，絕大部分來自于單產的提高。由于西藏自然地理、生態環境等因素制約，西藏耕地面積少、耕地質量較差，隨著農牧民教育程度的逐漸提高，科技投入不斷加大，西藏農業發展正在由增產向提質轉變。

　　3結論

　　本文初步分析了近30來西藏耕地面積及主要農作物時空變化特征，并利用主成分分析法分析了西藏耕地面積變化的驅動因素，主要結論如下：

　　(1)近30多年來，西藏地區耕地面積基本處于增長狀態，但增長幅度較小，近30年來，共增加14.58khm2;人均耕地面積逐年降低，較1985年相比，人均耕地面積降低35.71%;主要農作物單產較1985年有很大的提升，小麥、青稞單產均增產100%以上，豆類增產幅度較小，僅為58.13%;青稞播種面積最為穩定，且有逐年減少趨勢，小麥和油菜有逐年增加趨勢;近10年來，糧食主產區拉薩、日喀則、山南糧食總產在全區糧食總產中所占比例不穩定，且有逐年減少態勢;全區土地利用動態度變化不大，糧食主產區拉薩、日喀則變化相對較大。

　　(2)利用主成分分析將影響西藏耕地面積變化的13個驅動因子進行歸納總結，結果表明經濟發展、人口增加、農業發展水平是影響西藏耕地面積變化的主要驅動力。經濟發展和人口增長，導致城鎮化建設不斷擴大，基礎設施建設采用“占一補一”或“占一補二”的原則，占用耕地相對較少，且許多荒地得以開發利用，實現占補平衡，且補大于占;主要農作物單產水平的提高是糧食總產量增加的主要原因。

　　綜上所述，雖然耕地面積有所增加，但增長幅度較小，隨著社會經濟的發展，基礎設施建設的規模逐漸擴大，人民生活水平的提高，人口不斷地增長，不可避免的引起非農業建設對土地的需求，同時，為了保障生態安全，還林還草的面積有可能進一步增大。西藏后備耕地資源有限，糧食安全問題必須引起足夠重視。

　　參考文獻：

　　[1]周天財,孫建,于歡.1980-2013年青藏高原氣候時空變化特征研究[J].成都信息工程大學學報,2017,32(05):524-531.

　　[2]李萍.西藏糧食生產影響因素的灰色關聯動態分析[J].安徽農業科學,2009,37(27):13249-13251.

　　[3]楊春艷,沈渭壽,李海東.1985-2010年西藏糧食產量對氣候和耕地變化的響應[J].農業工程學報,2015,31(17):261-269.

　　[4]王建林,熊偉,胡單.近20年來西藏耕地與人口變化相關分析[J].資源科學,2004(01):45-51.

　　[5]郭志儀,曹建云.2006-2050年西藏人口發展趨勢預測[J].西藏大學學報(漢文版),2006(04):6-13+27

　　作者：張毅馬躍峰贠民政趙春會王建林