中國馬鈴薯區域分布影響因素分析*

楊亞東，羅其友，王道龍，劉 洋，倫閏琪

(中國農業科學院農業資源與農業區劃研究所，北京 100081)

0 引言

馬鈴薯與水稻、玉米和小麥共同作為世界四大糧食作物，每年約被世界上2/3的人口選作主糧[1]。2015年，我國馬鈴薯播種面積和產量分別為552萬hm2和1 897萬t，均占全球1/4左右[2-3]，位居世界首位。2015年1月，農業部開始推進馬鈴薯主糧化戰略[4-5]，我國馬鈴薯產業發展進入整合發展階段。2017年底召開的中央農村工作會議指出“堅持立足國內保障自給的方針，牢牢把握國家糧食安全主動權”[6]。馬鈴薯具有抗旱抗災、穩產高產等特征，作為干旱少雨、土地貧瘠的貧困地區的鐵桿作物，系統研究其區域分布的影響因素，對調節區域糧食安全，促進區域經濟協調發展具有重要意義。

已有文獻關于馬鈴薯種植區域分布影響因素的研究較少，多集中在市場供求的分析和解釋[7-9]、從品種、區域和產業等角度對馬鈴薯發展思路進行研究與比較[10-13]、利用空間分析方法對中國馬鈴薯的空間差異等進行研究[14-16]等。總體看來，相對缺乏馬鈴薯區域分布影響因素的系統關注?；诖?，文章通過定量分析自然、技術、經濟和政策等因素對馬鈴薯種植區域分布的影響及其作用途徑，為馬鈴薯種植的中長期區域布局提供參考。

1 中國馬鈴薯區域分布影響因素辨識及指標體系構建

1.1 影響因素辨識

(1)氣候因素

由于種植業對自然資源的依賴程度較高[17]，降水、光照、溫度和自然災害等情況影響馬鈴薯區域種植潛力[18-19]。1909年以來，中國氣候變暖的速率高于全球平均值，升溫速率為0.9～1.5℃/100年，且具有明顯的區域和季節變化特征[20]。氣候變暖使馬鈴薯最適宜、適宜和不適宜種植區面積縮小，次適宜區和可種植區面積擴大[21]。研究表明，馬鈴薯適宜種植區海拔高度提高了約100～200m，在我國西北地區種植海拔高度上限可達3 000m左右[22]。同樣，氣候變化使得隴中、寧夏中南部、固原中南部和彭陽大部地區可種植面積迅速擴大[23]。該文選取降水量、日照時數、平均氣溫、成災面積作為影響中國馬鈴薯區域分布氣候因素的二級指標，探索氣候因素對馬鈴薯種植區域分布的影響。

(2)技術因素

技術要素的集聚是中國馬鈴薯區域分布的重要推動力。一方面，技術創新提高了馬鈴薯作物的地域適應性，原有的不適宜種植區域馬鈴薯種植面積大幅增加。如廣西等地利用稻草覆蓋免耕栽培技術和稻草包芯高產栽培技術，在冬閑田上開展冬季和早春馬鈴薯生產，播種面積從2009年的幾乎沒有迅速發展到2015年的7.31萬hm2。另一方面，農業科技的研發與推廣提升了耕地綜合使用效益，使得馬鈴薯單產提升進而促進生產效益增加。如西北旱區探索出了馬鈴薯全膜覆蓋壟作集雨溝播栽培技術模式，通過試驗推廣，將有限的天然降水集流到壟溝，使天然降水的利用率達到60%～75%，有效克服了干旱日益嚴重、雨季推遲等不利氣候因素的影響。該文選取科普專職人員數量、科研經費內部支出、機械化水平作為影響中國馬鈴薯區域分布氣候因素的二級指標，探索技術因素對馬鈴薯種植區域分布的影響。

(3)經濟因素

經濟因素是農戶種植行為選擇的主要因素之一。根據《全國農產品成本收益資料匯編(2002-2016)》分析， 2001～2015年，全國馬鈴薯單位面積投入總成本和單位面積生產成本總體上均呈現上升趨勢，但部分地區種植凈利潤則呈現較大波動，如重慶馬鈴薯種植的成本利潤率(成本利潤率=凈利潤/成本費用總額×100%)在2003年和2010年出現了兩次較為明顯的波峰，分別為210.9%和185.32%。此外，與馬鈴薯生長季相似的替代作物，如北方的玉米和南方的油菜，也會影響馬鈴薯農戶的種植意愿。結合分析，該文選取每667m2總成本、每667m2凈利潤、替代作物價格(北方地區玉米、南方地區油菜的每50kg主產品平均出售價格)作為影響中國馬鈴薯區域分布經濟因素的二級指標，探索經濟因素對馬鈴薯種植區域分布的影響。

(4)政策因素

馬鈴薯種植相關鼓勵政策的推出，直接催化了我國馬鈴薯種植格局的明顯變化。在西北一季區和西南一二季作區*北方一季作區細分為東北一季區(黑龍江、吉林、內蒙古東部、遼寧的北部和西部)、華北一季區(內蒙古中西部、河北北部、山西中北部)和西北一季區(甘肅、寧夏、陜西西北部和青海東部)； 中原二季作區主要包括遼寧、河北、河南、山東、江蘇、浙江、安徽、江西和山西南部等地區； 西南一二季混作區主要包括云南、貴州、四川、重慶、西藏、湖南和湖北的西部地區以及陜西的安康市； 南方冬作區主要包括廣西、湖南、廣東、福建、海南、臺灣和江西南部、湖北的東部等省(區)，馬鈴薯產業被作為是旱地地區和貧困山區脫貧致富的重要產業和優勢產業，各級政府先后制定馬鈴薯產業發展的優惠政策。從初級加工、種薯貯藏、原種生產、質量體系和災害補貼等，政策體系基本實現了產前、產中和產后的全面覆蓋，直接激發了農戶的種植欲望。該文從省、市、縣三級政策推行情況，選取政策虛擬變量作為影響中國馬鈴薯區域分布的主要政策因素，探索政策因素對馬鈴薯種植區域分布的影響。

1.2 指標體系構建

該文通過走訪調研和梳理相關文獻，確立了中國馬鈴薯空間格局影響因素評價指標模型。主成分分析變量指標選取及說明如表1，以全國各省市馬鈴薯生產水平作為被因變量，用馬鈴薯種植面積(y)來表示。自變量包含4個子系統和11個指標。

表1 主成分分析變量指標選取及說明

一級指標二級指標 單位數據處理數據來源氣候因素(1～4)1生長季降水量mm研究區域氣象站點相關數據計算得出國家氣象信息中心標準氣象站2生長季日照時數h3生長季平均氣溫℃4馬鈴薯成災面積hm2相關調研數據計算得出國家馬鈴薯產業技術體系相關調研數據；2016年各省(區、市)統計年鑒技術因素(5～7)5科普專職人員數量人《中國科技統計年鑒2016》6科研經費支出萬元7馬鈴薯機械化水平相關調研數據計算得出國家馬鈴薯產業技術體系相關調研數據經濟因素(8～10)8每667m2總成本元《全國農產品成本收益匯編2015》；國家馬鈴薯產業技術體系相關調研數據9每667m2凈利潤元10替代作物價格元北方地區的玉米和南方地區的油菜政策因素(11)11政策虛擬變量(沒有=0，有=1)代表政策作用的虛擬變量

2 中國馬鈴薯區域分布影響因素定量評價

為更深入分析馬鈴薯種植區域分布的影響因素，在充分考慮影響因素的重要性、指標的完整性和數據的可得性等方面的前提下，從氣候、技術、經濟和政策4個方面選取指標，利用全國16省市相關數據進行主成分分析，以客觀反映各指標影響作用。

2.1 研究方法與數據來源

(1)研究方法

主成分回歸分析方法可以通過數據綜合處理，將多個指標減少為1個或幾個能夠反映原來多個指標所代表的絕大部分信息的綜合指標，在互不相關的前提下避免原始指標的重復性，以便于進一步的計算、分析和評價。通過求解觀測變量相關矩陣的特征方程，得到相關數據的特征值和單位特征向量。將特征值按從大到小排序，其分別代表主成分所解釋的觀測變量的方差。依據準則設定提取其中特征值大于1、對解釋原始數據方差信息大于85%貢獻率的主成分予以保留，并測算出綜合得分。

(2)數據來源

該文采用的數據資料主要有4個來源：(1)《全國農產品成本收益資料匯編》中各省(區、市)的馬鈴薯生產成本收益數據； (2)國家馬鈴薯產業技術體系調研數據； (3)國家氣象信息中心的研究區域氣象站點數據； (4)全國各省(區、市)政策性文件。

2.2 結果與分析

在實際應用主成分分析法進行分析問題時，往往會存在原始數據的不同指標的量綱不同，指標之間存在不可比性，這樣會對分析結果造成影響，因此對原始數據進行標準化處理以消除量綱影響。11個影響因素指標的KMO統計量的值為0.577，Bartlett球形檢驗結果顯著，說明11個指標之間并非獨立，即存在相關性，因此數據適合進行主成分分析。

表2 主成分的特征值和貢獻率

成分特征值貢獻率累計貢獻率F14 82934 49334 493F22 94721 04855 541F32 12215 15470 695F41 3419 58180 276

表3 主成分載荷矩陣

指標自然信息外部信息科技信息收益信息生長季降水量0 717-0 054-0 194-0 337生長季日照時數-0 675-0 1520 4590 433生長季平均氣溫0 869-0 143-0 2290 234成災面積-0 5950 520 4990 085科研經費內部支出0 552-0 4060 5420 029科普專職人員數量0 520 0270 5390 265機械化水平0 0040 8530 2050 209每667m2總成本0 5220 2090 3330 651凈利潤0 534-0 0310 036-0 735替代作物價格0 4040 629-0 284-0 505政策-0 587-0 844-0 4530 686

利用SPSS 軟件的主成分分析功能進行分析，計算標準化數據的特征值和方差貢獻率，如表2。根據特征根大于1、累計貢獻率大于80%的主成分選擇原則，選取特征值分別為4.83、2.95、2.12和1.34的4個主成分F1、F2、F3和F4，這4個主成分的累積貢獻率已達80.28%，包含了所有指標的絕大部分信息，而且這4個主成分互不相關，這樣就避免了評價指標信息重復的問題。因此，決定用這4個新變量來代替原來的11個變量。

從主成分載荷矩陣(表3)可以看出，馬鈴薯生長季的降水量、日照時數、平均氣溫和成災面積對第一主成分的載荷較高，這些指標基本反映了馬鈴薯生產的自然信息，這充分說明自然資源條件對馬鈴薯的生產起到了重要作用，良好的自然條件將有利于馬鈴薯種植生產； 機械化水平和政策在第二主成分上有較高載荷，這些指標基本反映了馬鈴薯生產外部信息，這說明馬鈴薯生產的技術條件、環境條件和政策條件影響力較強； 科普專職人員數量、科研經費內部支出在第三主成分上有較高載荷，這些指標基本反映了馬鈴薯生產過程中科技投入信息，表明科技投入能夠直接影響馬鈴薯農戶的生產決策，是影響區域馬鈴薯分布的關鍵因素。每667m2總成本、每667m2凈利潤和替代作物價格在第四主成分上有較高載荷，這些指標基本反映了馬鈴薯生產效益信息，這說明市場信息與馬鈴薯生產緊密相關，比較效益高有助于馬鈴薯生產的擴大，反之則會抑制其發展。

通過計算主成分載荷矩陣中的數據得到主成分表達系數，將系數與標準化后的數據相乘，得到4個主成分的得分。最后以4個主成分所對應的貢獻率與各個主成分得分相乘并求和計算主成分綜合得分模型：

F1=0.326X1-0.307X2+0.395X3-0.271X4+0.251X5+0.237X6+0.002X7+0.238X8+0.243X9+0.184X10-0.267X11

F2=-0.031X1-0.09X2-0.083X3+0.303X4-0.237X5+0.016X6+0.497X7+0.122X8-0.018X9+0.366X10-0.492X11

F3=-0.133X1+0.315X2-0.157X3+0.343X4+0.372X5+0.370X6+0.141X7+0.229X8+0.025X9-0.195X10-0.311X11

F4=-0.291X1+0.374X2+0.202X3+0.073X4+0.025X5+0.229X6+0.180X7+0.562X8-0.635X9-0.436X10+0.592X11

各信息綜合得分及其排名見表4。

表4 各省(區、市)主成分得分

地區自然信息排名外部信息排名科技信息排名收益信息排名總得分排名河北-0 31120 308-0 40110 76120 3511山西-1 3216-0 4016-1 04140 2215-2 5516內蒙古0 7631 0021 5421 1984 502遼寧-0 05100 2610-0 93120 97100 2512黑龍江0 654-0 0814-1 80160 7413-0 4913山東1 4113 1911 2945 07110 951湖北0 575-0 2915-0 19101 7341 829重慶0 4660 9730 0381 7833 246四川0 2380 8850 9651 2373 305貴州-0 18110 9040 7361 1792 648云南0 387-0 02131 5231 2363 127陜西-0 97140 1112-0 1591 4850 4610甘肅0 1990 6561 7410 90113 494寧夏0 8220 4470 4872 1723 913青海-1 22150 1711-1 11150 1716-1 9915新疆-0 81130 299-1 03130 6414-0 9114

由于該文選用的是影響因素指標，因此表4所示排名表示的是各地區馬鈴薯影響信息綜合得分。由表4結果可知，自然信息得分最高的是山東、寧夏、內蒙古、黑龍江，究其原因，這些省(區、市)在馬鈴薯生長季的氣溫、日照和降水相對更加適宜馬鈴薯生產，且自然災害較少； 外部信息得分最高的是山東、內蒙古、重慶、貴州，由于近年來山東馬鈴薯上市期較早，價格較其他作物更具優勢，農戶在收益較高的預期激勵下也增加了機械化方面的投入，部分西北一季區與西南一二級混作區將馬鈴薯作為重要的民生作物，政府更加注重制定優惠政策引導產業發展； 科技信息得分最高的是甘肅、內蒙古、云南，由于西北一季區和西南一二季混作區廣泛推廣集雨溝播、地膜覆蓋、高壟雙行、稻草包芯等栽培技術模式，提高了天然降水的利用率，擴大了可耕種的區域范圍。收益信息得分最高的是山東、寧夏、重慶，以山東為代表的中原二季作區每年鮮薯上市時間最早，起到了引導市場價格年度走向的關鍵作用，寧夏等北方一季作區則是面積最大的種薯、鮮薯產銷區域，價格能夠平衡全年的馬鈴薯走勢，以重慶為代表的西南一二季混作區的馬鈴薯價格售價相對較高。4個信息的綜合加權得分表示各省(區、市)影響信息的整體水平，綜合得分排名靠前的是近年來馬鈴薯生產呈上升趨勢的山東、內蒙古、寧夏、甘肅等地區。值得注意的是，由于較高的種植收益和冬閑田的充分利用，四川、重慶、云南、貴州等西南一二季混作區的馬鈴薯生產呈現良好發展勢頭，加之馬鈴薯增產潛力較大、用途較為廣泛，這些地區在農業結構調整過程中，也將馬鈴薯產業列為發展區域經濟和解決貧困地區農戶脫貧致富的重要措施。

3 結論與建議

3.1 結論

(1)在自然條件方面，中原二季作區和西北一季區具有自然優勢； 在機械化利用與政策促進方面，西北一季區、中原二季作區和西南一二季混作區較為主動； 在科技應用方面，西北一季區和西南一二季混作區效率相對較高； 在種植收益方面，中原二季作區、西北一季區和西南一二季混作區具有明顯優勢。

(2)綜合自然、技術、經濟、政策等4個方面，結合調研情況分析，以山東為代表的中原二季作區和以內蒙古、寧夏、甘肅為代表的北方一季作區作為傳統優勢產區，仍然具有較大的比較優勢。以四川、重慶等省(區、市)為代表的西南一二季混作區后發勢頭明顯。

3.2 建議

(1)以山東為代表的中原二季作區，建議發揮產業比較優勢，研究鮮食專用品種配套栽培技術，選育和推廣抗病專用型和適合外銷型品種，研究建立規?；?、標準化和機械化生產模式，發展以向東出口鮮薯為主導的馬鈴薯生產優勢區。

(2)以甘肅、寧夏為代表的西北一季區，建議壓縮菜用型馬鈴薯種植面積，在條件成熟的地區研究建立淀粉加工和食品加工等加工專用型原料薯生產優勢區。以河北、內蒙古的壩上地區為代表的華北一季區，建議發揮氣候、土壤和區位優勢，增加馬鈴薯生產科研力度，研究建立精品鮮薯和種薯生產優勢區。在黑龍江等東北一季區，雖然馬鈴薯生產并不具有綜合優勢，但這一地區具有土壤肥沃、日照充足和晝夜溫差大的天然優勢，建議研究嘗試建立“北出口、南調運”的有限發展區和種用、鮮食用馬鈴薯戰略儲備區，充實糧食安全總體發展戰略，為保障國家糧食安全提供堅實保障。

(3)以重慶、貴州等為代表的西南一二季混作區地勢復雜、氣候區域差異變化明顯，馬鈴薯種植能夠實現周年生產，建議選育專用品種，推廣新型高效種植模式，推行實用配套增產技術，擴大種植面積，研究建立覆蓋“東盟”、向南出口和輻射秦嶺—淮河以南地區的鮮薯和種薯為主導的馬鈴薯生產優勢區。

[1] 盧肖平. 馬鈴薯主糧化戰略的意義、瓶頸與政策建議.華中農業大學學報(社會科學版), 2015，(3)： 1～7

[2] 陳萌山， 王小虎.中國馬鈴薯主食產業化發展與展望.農業經濟問題， 2015, 36(12): 4～11

[3] 楊炳南， 張小燕，趙鳳敏，等.不同馬鈴薯品種的不同加工產品適宜性評價.農業工程學報， 2015, 31(20): 301～308

[4] 余欣榮. 以科技創新引領馬鈴薯主糧化發展.農業工程技術, 2015，(2)： 15～16

[5] 達存瑩， 岳云，賈秀蘋.甘肅省馬鈴薯產業發展現狀及主糧化發展分析.中國農業資源與區劃， 2016, 37(3): 38～42

[6] 中央農村工作會議報道，http：//www.xinhuanet.com/politics/2017-12/29/c_1122187923.htm， 2017-12-29

[7] 張萌， 羅其友，高明杰，等.馬鈴薯市場研究進展及展望.中國馬鈴薯， 2017, 31(2): 113～118

[8] 米健， 羅其友，高明杰，等.馬鈴薯中長期供求平衡研究.中國農業資源與區劃， 2015, 36(3): 27～34

[9] 高明杰， 張萌，羅其友.我國馬鈴薯價格波動的時空差異性研究.價格理論與實踐, 2017，(10)： 104～107

[10]楊炳南， 張小燕，趙鳳敏，等.不同馬鈴薯品種的不同加工產品適宜性評價.農業工程學報， 2015, 31(20): 301～308

[11]苗百嶺， 侯瓊，梁存柱.基于GIS的陰山旱作區馬鈴薯種植農業氣候區劃.應用生態學報， 2015, 26(1): 278～282

[12]王鶴齡， 王潤元，張強，等.甘肅馬鈴薯種植布局對區域氣候變化的響應.生態學雜志， 2012, 31(5): 1111～1116

[13]謝從華. 馬鈴薯產業的現狀與發展.華中農業大學學報(社會科學版), 2012，(1)： 1～4

[14]呂超， 耿獻輝，張榮敏.中國馬鈴薯種植的時空布局演化研究.統計與決策, 2017，(18)： 136～140

[15]楊亞東， 胡韻菲，栗欣如，等.中國馬鈴薯種植空間格局演變及其驅動因素分析.農業技術經濟, 2017，(8)： 39～47

[16]Peng Yu ，Wang Qinghui，Fan Min.Identification of the key ecological factors influencing vegetation degradation in semi-arid agro-pastoral ecotone considering spatial scales.Acta Oecologica， 2017， 85： 62～68

[17]肖衛東. 中國種植業地理集聚：時空特征、變化趨勢及影響因素.中國農村經濟, 2012，(5)： 19～31

[18]Zhao Junfang，Zhang Yanhong，Qian Yonglan，et al.Coincidence of variation in potato yield and climate in northern China.Original Research Article Science of The Total Environment， 2016， 573： 965～973

[19]N.Ibanez，F.J.Berli，R.W.Masuelli，et al.Influence of altitude and enhanced ultraviolet-B radiation on tuber production，seed viability，leaf pigments and morphology in the wild potato species Solanum kurtzianum Bitter & Wittm collected from an elevational gradient.Original Research Article Plant Science， 2017， 261： 60～68

[20]第三次氣候變化國家評估報告編寫委員會. 第三次氣候變化國家評估報告.北京：科學出版社， 2015.162～165

[21]姚玉璧， 雷俊，牛海洋，等.氣候變暖對半干旱區馬鈴薯產量的影響.生態環境學報， 2016, 25(8): 1264～1270

[22]肖國舉， 仇正躋，張峰舉，等.增溫對西北半干旱區馬鈴薯產量和品質的影響.生態學報， 2015, 35(3): 830～836

[23]Zhang Youliang，Wang Fengxin，Clinton Cleon Shock，et al.Effects of plastic mulch on the radiative and thermal conditions and potato growth under drip irrigation in arid Northwest China.Original Research Article Soil and Tillage Research， 2017， 172： 1～11