河南省縣域糧食生產格局變化及其影響因素

薛劍，韓娟，劉玉，張鳳榮

(1.中國農業(yè)大學土地資源管理系，北京100094;2.國土資源部土地整治中心，北京100035;3.地質出版社，北京100083;4.國家農業(yè)信息化工程技術研究中心，北京100097)

糧食安全事關國計民生，一直是全球關注的熱點問題[1]。國內外學者基于全國、流域和省域等不同尺度，從糧食總產、糧食單產、人均糧食占有量、耕地變化等方面圍繞糧食生產與糧食供應安全開展了大量研究[2-3]。宏觀視角的研究表明，我國傳統(tǒng)的“南糧北調”格局被“北糧南調”格局所取代，北方糧食生產全面超越南方[4];糧食生產重心呈現(xiàn)“北進中移”的遷移態(tài)勢，重心移動速率不斷增大[5];案例研究逐步由東部沿海省份和東北地區(qū)向中西部省份轉移，基于縣域單元的省域研究日益增多[6-7]。這些研究為我國糧食生產布局優(yōu)化、糧食發(fā)展政策制定提供了科學依據(jù)和經驗借鑒。然而，現(xiàn)有研究往往忽視空間因素對糧食生產區(qū)域格局形成的影響，研究成果的應用性受到較大限制[8]。隨著GIS技術的發(fā)展，運用探索性空間數(shù)據(jù)分析法(ESDA)探討省域內縣域糧食生產格局變化的研究增多[9-10]，從技術手段上支持了深層次的糧食生產空間格局研究。

河南省是我國13個糧食主產省份之一，也是《全國新增1 000億斤糧食生產能力規(guī)劃(2009—2020年)》中確定的重點區(qū)域之一，對保障我國糧食安全、穩(wěn)定糧食價格發(fā)揮著重要作用。隨著《中原經濟區(qū)規(guī)劃(2012—2020年)》的出臺，中原經濟區(qū)正式上升為國家戰(zhàn)略，河南省將全力打造全國“三化”協(xié)調發(fā)展的示范區(qū)。隨著河南省城市化、工業(yè)化和農業(yè)產業(yè)結構調整進程的加快，耕地變化等對糧食生產的資源性約束也將日益突出。為此，本研究選取1980，1990，2000和2010年為研究斷面，運用探索性空間數(shù)據(jù)分析技術、重心模型等方法，探討河南縣域糧食生產的時空動態(tài)及其驅動機制，以期為糧食生產布局優(yōu)化、糧食生產政策制定提供科學的決策依據(jù)。

1 研究區(qū)概況

河南省地處我國中東部、黃河中下游，位于31°23'N～36°23'N、110°21'E～116°39'E之間，與江蘇、山東、安徽、河北、山西、陜西、湖北相臨，土地面積16.70×104km2;地勢西高東低，地貌類型復雜多樣;屬暖溫帶—亞熱帶、濕潤—半濕潤季風氣候，年平均氣溫12.4℃～16.1℃，年均降雨量432～1 632 mm，氣候溫和，無霜期較長，適宜種植多種農作物。在《農用地分等規(guī)程》(TD/T 1004—2003)中，根據(jù)熱量、水分、地貌以及社會經濟條件，將河南省劃分為燕山太行山山前平原區(qū)、冀魯豫低洼平原區(qū)、黃淮平原區(qū)、鄂豫皖丘陵山地區(qū)、豫西山地丘陵區(qū)5個耕作制度二級類型區(qū)(表1)。全省共設17個設區(qū)市，21個縣級市，88個縣，50個市轄區(qū)。從數(shù)據(jù)的可獲得性和研究需要出發(fā)，以2010年為基準對行政區(qū)劃單元進行修正，將各地級市轄區(qū)進行歸并，如鄭州市中原區(qū)、二七區(qū)、管城區(qū)、金水區(qū)、上街區(qū)、惠濟區(qū)統(tǒng)一歸并為鄭州市區(qū)，最終選取126個單元。1980年的糧食產量等數(shù)據(jù)由中國自然資源數(shù)據(jù)庫提供，1990，2000和2010年數(shù)據(jù)來源于相應年份的《河南統(tǒng)計年鑒》、《河南農村統(tǒng)計年鑒》、《河南改革開放30年》。

表1 河南省耕作制度分區(qū)Tab.1 Division of standard farming system in Henan

2 研究方法

2.1 探索性空間數(shù)據(jù)分析法

探索性空間數(shù)據(jù)分析ESDA是一系列空間分析方法和技術的集合，通過對事物或現(xiàn)象分布格局的描述與可視化探究區(qū)域屬性值的分布模式和空間相對差異?？臻g自相關分析分為全局空間自相關和局部空間自相關分析。本研究采用Global Moran’sI和Local Moran’sI統(tǒng)計量進行空間關聯(lián)分析[11]。

2.2 重心遷移模型

重心遷移模型通過計算特定區(qū)域某一時段糧食產量的分布重心來描述糧食生產的時空分布格局[5]。

式中:Xj，Yj分別表示第j年糧食產量的區(qū)域重心坐標;Xi，Yi分別表示i縣域的重心坐標;Gij表示i縣域第j年的糧食產量;dm為重心遷移距離(km)。

2.3 空間回歸分析法

當統(tǒng)計指標存在空間自相關時，殘差就不再相互獨立，因此，存在空間自相關的數(shù)據(jù)不適合采用普通最小二乘法(OLS)回歸模型。將空間相互關系引入模型，對線性回歸模型通過空間權重矩陣w進行修正的空間計量經濟模型能較好的解決這一問題，其中空間滯后模型(spatial lag model，SLM)和空間誤差模型(spatial error model，SEM)較為常用[11]。

3 結果與分析

3.1 糧食生產總體特征

1980年以來，河南省糧食產量及其在全國糧食總產量中的比例均呈現(xiàn)增長態(tài)勢(圖1)。2010年全省糧食產量為5 437.1萬t，是1980年的2.53倍，年均增長率為3.14%;人均糧食占有量由1980年的295 kg/人增加至2010年的543 kg/人，年均增長率2.05%;糧食產量在全國糧食總產量中的比例也呈現(xiàn)增長趨勢，由1980年的6.88%增加到2000年的8.87%，2001—2005年在9.00%上下徘徊;2006年以后基本維持在10%。

圖1 1980—2010年河南省糧食產量及其占全國比例Fig.1 The grain yield of Henan Province and its proportion in China from 1980 to 2010

從不同耕作區(qū)的糧食產量分布看(圖2)，黃淮平原區(qū)糧食產量所占比例最大，1980年以來占河南糧食總產量的50%以上，31年間比例提高了4.43%;冀魯豫低洼平原區(qū)所占比例由9.48%提高到12.29%;鄂豫皖丘陵山地區(qū)糧食產量比例略有下降;燕山太行山山前平原區(qū)、豫西山地丘陵區(qū)下降明顯，分別降低了2.64%和3.70%，河南省糧食生產向平原區(qū)集中分布的趨勢明顯。

圖2 1980—2010年河南省不同耕作區(qū)糧食產量分布Fig.2 Distribution of grain production in different farming system divisions in Henan Province from 1980 to 2010

1980年糧食產量Gi＞30萬t的縣域為6個，2010年為81個，糧食產量占全省總產量的比例由9.03%增至86.28%。1980—2010年間，Gi＞50萬t的縣域由0個增至60個，糧食產量比例增至72.74%，主要分布在黃淮平原區(qū)、冀魯豫低洼平原區(qū)和鄂豫皖丘陵山地區(qū);40萬～50萬t的縣域由0個增至7個，縣域糧食產量比例增加至5.33%;30萬～40萬t的縣域由6個增至14個，但縣域糧食產量比例由9.03%下降至8.21%，分布較分散;Gi＜30萬t的縣域個數(shù)和糧食產量比例均下降，集中分布在豫西山地丘陵區(qū)以及鄂豫皖丘陵山地區(qū)西部。與1980年相比，2010年118個縣域糧食產量增加，其中增加量＞45萬t的縣域27個，在黃淮平原區(qū)分布密集，糧食生產向產糧大縣集中;8個縣域糧食產量下降(表2)。95個糧食主產縣1980，1990，2000和2010年糧食產量占全省糧食產量的比例分別為81.41%，86.27%，88.01%和89.78%，呈逐年增加的趨勢?？梢姡幽鲜】h域糧食總產量呈現(xiàn)“東高西低”的分布格局，糧食生產向產糧大縣和東部平原區(qū)集中，且隨時間推移，這種分異格局愈加明顯(圖3)。

表2 1980—2010年河南省糧食縣域生產分布Tab.2 Distribution of counties with different grain production in Henan Province from 1980 to 2010

圖3 1980—2010年河南省糧食總產量及其變化Fig.3 Spatial patterns of grain production and changes at county level in Henan Province from 1980 to 2010

3.2 縣域糧食空間差異

采用GeoDa軟件計算了4個年份縣域糧食生產的全局自相關系數(shù)Moran’sI。整個研究期間，全局Moran’sI估計值全部為正，從1980年的0.142 8增加到2010年的0.422 7，總體上呈現(xiàn)擴大趨勢。表明1980年以來河南省各縣域糧食生產存在顯著的正空間自相關特性，即糧食產量水平較高的縣域趨于相鄰，糧食產量水平較低的縣域也趨于相鄰，呈現(xiàn)出空間上高高集聚(H-H)、低低集聚(L-L)的集中分布。同時，隨著時間的推移，這種正的空間自相關增強的趨勢表明同期河南省糧食生產增長差異也存在相似的趨勢。

由于Moran’sI估計值是一種總體統(tǒng)計指標，只能顯示某一區(qū)域與其周邊地區(qū)空間差異的平均程度，不能全面反映區(qū)域糧食的局部空間特征及其空間關系。因此，需要借助Moran散點圖和LISA聚集圖進一步分析縣域糧食生產的局部空間特性。根據(jù)縣域糧食生產數(shù)據(jù)，得到河南省1980，1990，2000和2010年糧食生產的Moran散點圖，統(tǒng)計出在p=0.05顯著性水平下的4種空間聯(lián)系類型的縣域數(shù)。31年間，位于H-H和L-L象限的縣域數(shù)量由75個增至97個，1980，1990，2000和2010年位于H-H和L-L象限的縣域分別占縣域總數(shù)的59.5%，70.6%，73.8%和77.0%，位于H-H象限的縣域糧食產量分別占河南省糧食總產量的54.2%，59.0%，58.7%和61.8%，表明河南省糧食生產的局部空間自相關顯著，且集聚能力不斷增強。

圖4 1980—2010年河南省縣域糧食產量的LISA聚集圖Fig.4 LISA cluster map for grain production at county level from 1980 to 2010 in Henan Province

將Moran散點圖與LISA顯著性水平相結合，繪制p=0.05顯著性水平下的Local Moran’sI圖，進一步分析河南省縣域糧食生產的空間自相關程度及其演化過程(圖4)。①顯著的H-H類型區(qū)在研究期內有一定變化。1980年，H-H類型區(qū)主要分布在周口、南陽和許昌。隨著時間推移，以周口為核心的H-H類型區(qū)呈現(xiàn)“北進南擴”的趨勢，范圍明顯擴大，其他的呈現(xiàn)分散分布且不穩(wěn)定。該區(qū)位于黃淮平原區(qū)，是全國重要的糧食產區(qū)，光熱水土資源匹配較好，農業(yè)基礎設施較為完善，農用地綜合生產能力較高。②顯著的L-L類型區(qū)在洛陽、三門峽和南陽西北部呈現(xiàn)面狀分布，并逐步向平頂山和鄭州擴展。該類型區(qū)處于豫西山地丘陵區(qū)，地勢起伏度較大，水土流失較為嚴重，耕地地力差，農田基礎建設滯后，耕地綜合生產能力低。③顯著的低高(L-H)類型區(qū)和顯著的高低(H-L)類型區(qū)分布的縣域較少，規(guī)律性不甚明顯。

3.3 糧食生產重心變化

利用公式(1)～(3)計算出河南省4個研究年份的糧食生產重心，河南省糧食生產重心在許昌縣和鄢陵縣境內移動。2010年與1980年相比較，向東南方向偏移了26.07 km，經向移動距離遠大于緯向移動距離。這主要是由于河南東部的黃淮平原區(qū)糧食生產占全省的比例較高，且呈現(xiàn)增加趨勢，而豫西山地丘陵區(qū)糧食生產比例持續(xù)下降，共同造成了河南省糧食生產重心向東南方向偏移。重心模型的分析結果與空間自相關分析的結果一致，即河南省糧食生產集中于黃淮平原區(qū)，表現(xiàn)出一定的空間集聚性。

3.4 糧食生產空間格局變化的驅動機制

影響糧食生產空間格局變化的主要因素包括自然因素和經濟社會投入因素?；诂F(xiàn)有研究成果[9，12]，并咨詢相關領域的專家，本研究選擇了農林牧漁業(yè)勞動力、農業(yè)機械總動力等12個主要因素;鑒于河南省縣域糧食生產存在較強的空間自相關性，因此，利用126個縣域糧食產量數(shù)據(jù)與12個主要影響因素進行空間回歸分析。初步計算結果顯示:3個年份SLM和SEM模型的擬合系數(shù)均高于經典回歸模型，且以SEM的擬合系數(shù)最高，因此，以1980年的縣域糧食生產為基礎，以SEM模型的估計結果，分析1990，2000和2010年縣域糧食生產格局演化的影響因素(表3)。

由表3可知，1990年化肥投入量、農村居民人均純收入、農林牧漁勞動力是影響縣域糧食產量格局的主要因素;2000年演變?yōu)檗r林牧漁業(yè)勞動力、復種指數(shù)、種植結構、灌溉面積比率、農用機械總動力、化肥投入量、耕地面積和前期糧食產量等多個因素共同影響的特征;2010年，耕地面積和前期糧食產量的回歸系數(shù)分別提高到0.509和0.486，農林牧漁業(yè)勞動力對糧食產量的影響由正向顯著轉化為負向顯著，表明農林牧漁業(yè)勞動力的過剩阻礙了糧食生產的進一步提高;2010年，產業(yè)結構對糧食生產的影響為負，且通過顯著性檢驗，表明非農產業(yè)的快速發(fā)展對糧食生產的負面影響已經顯現(xiàn)?？傮w而言，隨著糧食播種面積單產提升幅度的降低，糧食總產格局更多地依賴于耕地面積、復種指數(shù)和種植結構，農藥、化肥等的投入因素的影響作用逐漸降低。

表3 空間殘差回歸分析結果Tab.3 Estimated results of SEM in 1990，2000 and 2010

4 結論與討論

1)研究期內河南省糧食總產量及其在全國糧食產量中的比例增加，糧食生產向產糧大縣、黃淮平原區(qū)集中;糧食生產重心東南方向偏移了26.07 km，糧食生產“東高西低”的態(tài)勢更加顯化。

2)1980—2010年河南省縣域糧食生產具有明顯的正相關性，全局空間自相關系數(shù)增加，糧食生產存在顯著的空間聚集分布;1980，1990，2000和2010年位于HH和L-L象限的縣域分別占總單元數(shù)的59.5%，70.6%，73.8%和77.0%，局部集聚格局顯著，集聚能力不斷增強;顯著H-H類型區(qū)主要分布在周口、南陽和許昌，顯著L-L類型區(qū)在洛陽、三門峽和南陽西北部呈現(xiàn)面狀分布。

3)前期糧食產量、耕地面積、復種指數(shù)和種植結構等對河南糧食生產格局的影響日益提高，而農林牧漁勞動力、化肥投入、農藥投入等的影響不顯著或者呈現(xiàn)負影響。在現(xiàn)階段，穩(wěn)定耕地面積并適當提高復種指數(shù)和糧食作物比例，優(yōu)化化肥、農藥等的投入結構，是提高糧食產量的重要途徑。

4)本研究采用空間分析和空間回歸分析等方法，揭示了1980—2010年河南省縣域糧食生產的時空演化特征，并初步分析了糧食生產時空演化的驅動機制。受數(shù)據(jù)資料和文章篇幅的限制，本研究在指標選取上未考慮溫度、降水等氣候因素;只是概略地分析了投入要素的數(shù)量，而對要素投入質量和結構的分析明顯不足。因此，針對不同類型區(qū)選取典型縣域，綜合考慮氣溫、降水等自然因子，深入研究要素投入數(shù)量、質量、結構等對糧食產量的影響，是下一步亟需研究的重要方向。

[1]解宗方，張偉.提高河南省糧食綜合生產能力的戰(zhàn)略途徑[J].農業(yè)現(xiàn)代化研究，2011，32(4):395-399.

[2]解宗方，張偉，滕永忠.基于資源潛力開發(fā)的糧食增產途徑——以河南省為例[J].地域研究與開發(fā)，2012，31(5):126-130.

[3]辛良杰，李秀彬，談明洪.中國區(qū)域糧食生產優(yōu)勢度的演變及分析[J].農業(yè)工程學報，2009，25(2):222-227.

[4]屈寶香，張華，李剛.中國糧食生產布局與結構區(qū)域演變分析[J].中國農業(yè)資源與區(qū)劃，2011，32(1):1-6.

[5]劉彥隨，翟榮新.中國糧食生產時空格局動態(tài)及其優(yōu)化策略探析[J].地域研究與開發(fā)，2009，28(1):1-5.

[6]李裕瑞，劉彥隨，龍花樓.江蘇省糧食生產時空變化的影響機制[J].地理科學進展，2009，28(1):125-131.

[7]石淑芹，陳佑啟，姚艷敏，等.東北地區(qū)耕地變化對糧食生產能力的影響評價[J].地理學報，2008，63(6):574-586.

[8]楊春，陸文聰.基于空間計量經濟模型的縣域糧食生產區(qū)域格局研究[J].農業(yè)技術經濟，2010(5):24-29.

[9]劉玉，郭麗英，劉彥隨.1980—2008年環(huán)渤海地區(qū)縣域糧食的時空動態(tài)變化及分析[J].農業(yè)工程學報，2012，28(12):230-236.

[10]葉長盛，朱傳民.江西省縣域糧食生產地域格局變化研究[J].農業(yè)現(xiàn)代化研究，2011，32(3):315-319.

[11]王法輝.基于GIS的數(shù)量方法與應用[M].姜世國，滕駿華，譯.北京:商務印書館，2009.

[12]賀振.河南省糧食生產動態(tài)變化及驅動機制分析研究[J].干旱區(qū)資源與環(huán)境，2010，24(6):6-9.