環京津貧困區耕地資源與農民收入的關系
——基于保定市7個貧困縣的面板數據分析

趙 麗，張蓬濤，朱永明

（河北農業大學 國土資源學院，河北 保定071001）

耕地作為人類社會生存和發展的重要物質資源，其利用和保護已引起政府和學者們的普遍關注。一些學者針對耕地資源變化與經濟因素的相關性進行了研究，如曲福田、吳麗梅［1］提出經濟增長與耕地非農化的庫茲涅茨曲線假說并對中國經濟發達典型區進行了驗證；許廣月［2］從中國省級面板數據角度對耕地資源與經濟增長之間的關系進行了實證分析；李兆富、楊桂山［3］對蘇州市近50a的耕地資源變化與經濟發展關系進行了研究，學者們［4－7］針對不同區域對耕地資源與經濟增長關系展開了相關研究，但針對耕地數量與農民收入的關系研究較少，而涉及環京津貧困區的該項研究就更少。環京津貧困區受經濟、社會、生態、政策等多方面的影響，土地利用變化更為劇烈，同時農民賴以生存的耕地資源面臨著來自城鄉建設、糧食生產、環境保護等多方面的需求壓力，因此研究該區的耕地資源變化與農民收入的關系意義更為重大。鑒于此，本研究以保定市7個貧困縣作為研究對象，對其耕地數量與農民收入之間的關系展開研究，以期為該區域政府制定相關發展政策提供實證參考依據。

1 數據來源與研究方法

1.1 數據來源與研究方法

考慮到數據的可獲取性，選取耕地面積數量和農民人均年純收入進行研究。數據來自《河北經濟年鑒（1997—2011年）》和《河北農村統計年鑒（1997—2011年）》。為了便于分析，用gd表示耕地面積，用sr表示農民人均年純收入，為消除二者數值之間差距，分別對2組數據取對數，取對數后分別記做lngd和lnsr。

考慮到保定市7個貧困縣的自然條件及區位差異，首先采用系統聚類法［8－9］，以2010年7縣的耕地面積數量和農民人均年純收入為研究對象對其進行聚類分析，按其聚類結果并結合區域實際情況將其分為3個區域。在此基礎上，對該3個區域進行耕地數量和農民收入的關系研究。根據對已有數據和現有文獻的研究，采用面板數據模型進行實證研究，主要分3步進行：首先，利用面板數據單位根檢驗方法對其進行單位根檢驗，采用3種面板數據單位根檢驗方法，即 LLC（Levin—Lin—Chu）檢驗、Fisher—ADF檢驗和Fisher—PP檢驗［10］，如果在檢驗中均拒絕存在單位根的原假設，則該序列是平穩的，反之則不平穩；在此基礎上，利用面板數據協整檢驗方法，對其進行長期均衡關系檢驗，即協整檢驗。利用EG兩步法的思想，將時間序列的協整檢驗擴展到面板數據中，分別稱為Pedroni檢驗和 Kao檢驗［11－12］，兩者在實證中均有相當廣泛的應用。采用Panel—ADF，Group—ADF，Panel—PP，Group—PP統計量和Kao檢驗對其進行協整檢驗［13－14］，最后根據面板數據的模型估計理論對3個區域的耕地數量與農民收入關系建立回歸模型。模型估計具體理論可參見文獻［15］。

1.2 研究區概況

保定市位于河北省中部，太行山北部東麓，冀中平原西部，北鄰北京市和張家口市，東接廊坊市和滄州市，南與石家莊市和衡水市相連，西部與山西省接壤。地處京、津、石三角腹地，市中心北距北京140 km，東距天津市145km，西南距河北省會石家莊市125km，直接可達首都機場、正定機場及天津、秦皇島、黃驊等海港。保定地區在地貌類型上是兼有平原、湖泊、濕地、丘陵、山地、亞高山草甸的地區。研究區7個縣分別是淶源縣、順平縣、唐縣、阜平縣、淶水縣、曲陽縣和易縣。其中淶源縣、順平縣、唐縣、阜平縣被定為國家級貧困縣，淶水縣、曲陽縣和易縣同上述4個縣共同被列入國家燕山—太行山連片特殊困難地區。2010年7個縣的耕地總面積為1.77×105hm2，人均耕地面積為599.84m2；國內生產總值共2.94×1010元，人均國內生產總值9 957.30元；財政收入共120 041萬元，人均財政收入為406.92元；農民人均年純收入為2 844.14元，經濟發展水平較低，與河北省的平均水平相比有很大差距，與北京市經濟發展水平相比則相差更遠。7縣的人均國內生產總值、人均財政收入和農民人均年純收入分別是河北省的34.7%，22.0%和47.7%，分別是北京市的13.1%，3.4%和21.4%，分別是天津市的13.6%，4.9%和28.2%。

2 結果與分析

2.1 研究區耕地數量與農民收入變化

2.1.1 耕地數量變化 從研究區耕地面積曲線圖（圖1）可以看出，在1996—2010年保定市7個貧困縣的耕地面積變化大致可分為3個階段：1996—2001年為變化平穩期，7個縣的耕地面積曲線均表現出變化平穩，在這期間各縣的耕地面積變化不大；2001—2008年為變化起伏期，7縣的耕地面積均表現出一定的波動起伏，波動的變化曲線在這期間表現出較高的一致性，即自2002年耕地面積開始減少，降至2003年，之后至2005年面積有一定程度的增加，2005年之后面積出現新一輪下降。其中淶源縣、曲陽縣和易縣3縣的波動幅度比其它縣較大。這期間的面積起伏變化與相應的國家政策和經濟發展有一定關系。經對其耕地來源及去向分析，耕地面積的減少主要源于兩個方面：一是2002年以來國家推行的退耕還林政策的全面實施，使一部分耕地向林地和園地轉移；二是非農建設用地占用。耕地面積的增加則主要來自對農村居民點整理和未利用地的開墾。2008—2010年為變化相對平穩期，2008年之后各縣的耕地面積變化幅度較上一時段平穩。除順平縣外，其余縣耕地面積在這期間均有一定程度的增加，但增長幅度不大。由于2008年各縣的退耕工作基本完成，補償標準相對較低（黃河流域及北方地區每1hm2退耕地每年補助現金1 050元。原每1hm2退耕地每年300元生活補助費，繼續直接補助給退耕農戶，并與管護任務掛鉤），由于利益驅動，有部分的退耕地重新被復耕，使耕地面積有所增加。

圖1 1996－2010年保定市7個貧困縣耕地面積變化

2.1.2 農民收入變化分析 研究區各縣的農民人均純收入均呈上升趨勢，但上升速度均相對較慢（圖2）。尤其淶源縣1996年農民人均純收入為1 119元，至2010年收入為1 998元，15a間農民人均純收入僅增長了879元，年均增長率4.23%。7縣中年均增長率最大的為阜平縣，由1996年的858元上升至2010年的2 460元，年均增長率7.81%。其余縣的農民人均純收入年均增長率基本保持在4%～6%之間，農民人均收入增長緩慢。

圖2 1996－2010年保定市7個貧困縣農民人均純收入變化

2.2 耕地數量與農民收入的關系分析

2.2.1 聚類分析 根據研究重點和考慮到數據的時效性，選取保定市7個貧困縣的耕地面積和農民人均純收入兩個指標，對其2010年指標數據進行類別劃分。采用SPSS 13.0軟件聚類分析法。采用類間距離系統聚類進行分析處理。根據分類結果（圖3）將研究區分為3個類別區，或稱3個區域：順平縣、阜平縣、淶源縣為一類（山區）；唐縣、曲陽縣為一類（環京區域）；淶水縣和易縣為一類（近京區域）。

圖3 保定市7個貧困縣重新調整距離后的聚類結果

2.2.2 面板數據單位根檢驗 利用Eviews 6.0軟件分別對各區域面板數據進行單位根檢驗。選擇了LLC檢驗、Fisher—ADF檢驗和Fisher—PP檢驗，得出變量的單位根檢驗結果詳見表1。從表1的檢驗結果可以看出，各區域lngd和lnsr在經過一階差分后在5%顯著性水平下都拒絕原假設。即耕地面積和農民人均純收入經過一階差分后都是平穩的，可以進行后續檢驗。

2.2.3 協整檢驗 為了直接驗證保定7個貧困縣耕地面積和農民人均純收入兩變量之間是否存在協整關系，采用Pedroni協整檢驗方法和Kao檢驗法。由于所涉及的時間跨度少于20a，故主要依靠Panel—ADF，Group—ADF，Panel—PP和 Group—PP這4個統計量來做基本判斷，并結合Kao檢驗。1996—2010年保定7個貧困縣耕地面積和農民人均純收入兩變量的對數值之間的Pedroni協整檢驗和Kao檢驗結果詳見表2。根據表2中的協整檢驗結果，模型1和模型2的Panel—ADF和Group—ADF的統計量指標和Kao檢驗法統計指標全部在5%的顯著性水平下拒絕了“不存在協整關系”的原假設。據此判斷，1996年以來，保定市7個縣的耕地面積與農民人均年純收入之間存在協整關系。

2.2.4 模型估計 基于上面的分析結果，模型的時間序列變量均通過平穩性檢驗，耕地面積與農民人均純收入之間存在著長期的均衡關系，則可利用面板模型估計各地區二者之間關系。通過F檢驗和Hausman檢驗最終確定用個體時點固定效應模型［16］。回歸采用個體和時期固定影響變截距模型。

表1 各類區域面板單位根檢驗結果

表2 各類區域協整檢驗結果

（1）山區區域耕地面積與農民人均純收入關系分析。從表3的回歸結果可以看到，山區區域常數項在1%的顯著水平下通過顯著性檢驗，可決系數R2較高，表示模型擬合優度較好。回歸模型為：

山區區域：lnsr＝3.958 5＋0.362 6lngd（1）

模型估計結果顯示，山區（淶源縣、順平縣、阜平縣）農民人均純收入與耕地面積呈正相關關系，環京和近京區域人均純收入與耕地面積呈負相關關系。模型（1）表明，山區區域耕地面積的增加會促使農民人均純收入增加，耕地面積每增加1%，農民人均純收入增加0.362 6%。從表4的個體影響中可以看出，阜平縣的截距項最大，其次順平縣，淶源縣的截距項最小，說明在山區區域的這3個縣中，阜平縣農民收入與耕地面積無關的部分最大，其次順平縣，淶源縣最小。從時期影響可以看出（表5），隨著年份的增加，截距項越來越大，說明農民人均年純收入在逐年增長。

分析該區域實際情況，農民主要收入來源為耕地，故對耕地的開發利用較為重視，收入的提高激勵農民開發出更多的耕地，這一區域在2004和2005年期間通過開發未利用地新增耕地2 539.83hm2，同時注意到區域的后備土地資源有限和建設用地占用，僅靠增加耕地來增加收入的渠道并不可行，在后期逐漸顯現，農民人均純收入的增長速度減緩，在2009年表現更為突出，增加農民收入還需要多渠道的收入來源。由于該區域典型地貌特征為山區，尤其淶源縣為全山區縣，因此耕地在數量、質量上均受地理因素限制，從而影響農民收入的增長速度。

表3 山區區域模型回歸分析結果

表4 山區區域模型個體影響：各縣在總體常數項基礎上的變化程度

（2）環京區域與近京區域耕地面積與農民人均純收入關系分析。從表6—7的回歸結果可以看到，環京區域常數項在5%的顯著水平下，近京區域常數項在1%的顯著水平下，均通過了顯著性檢驗。

回歸模型為：

表5 山區區域模型時期影響：各年份在總體常數項基礎上的變化程度

表6 環京區域模型回歸分析結果

表7 近京區域模型回歸分析結果

模型（2）—（3）表明，環京區域（唐縣、曲陽縣）和近京區域（淶水縣、易縣）耕地面積的增加未能帶來農民人均純收入的提高，或者說耕地面積減少反而會增加農民收入。

通過對該兩個區域實際情況分析，一是由于其相對有利的區位條件和相關產業的發展，促進了農民的出外打工；二是由于耕地面積的不斷減少，迫使農村剩余勞動力向城市轉移。農民的收入來源、結構在發生著變化，說明這兩個區域僅僅靠增加耕地面積來提高農民人均純收入已行不通。從另一側面也反映出該區域由于經濟建設對耕地的占用，農民必然要尋找新的收入來源。

從回歸分析的個體影響來看（表8—9），在環京區域，曲陽縣相對唐縣農民收入與耕地面積無關的部分較大；近京區域，則表現為易縣相對淶水縣農民收入與耕地面積無關的部分較大。其農民收入與耕地面積的相關性大小除與區域的農民收入來源和收入結構有一定關系外，在一定程度上也可以體現出該區域的經濟發展狀況。

表8 環京區域模型個體影響：各縣在總體常數項基礎上的變化程度

表9 近京區域模型個體影響：各縣在總體常數項基礎上的變化程度

從回歸分析結果的時期影響來看（表10），環京區域隨著年份的增加，截距項盡管總體趨勢是上升的，但其變化波動較大，說明農民人均年純收入的增長并不穩定。主要原因在于這兩個縣耕地面積減少，農民努力擴展就業渠道來增加收入，但受文化、技術水平等因素影響，呈現出就業乏力，難以使收入維持在一個相對穩定的水平，所從事的工作受季節、外部環境等因素影響較大。

表10 環京區域模型時期影響：各年份在總體常數項基礎上的變化程度

近京區域（表11）隨著年份的增加，截距項直線增加，其增長幅度比另外兩個區域高出很多，說明該區域的農民人均純收入在逐年快速增長，并且這種增長幅度高于山區區域和環京區域的增長幅度。相比而言，近京區域的地理區位優勢更為突出，為農民收入來源的多樣化提供了更為有利的區位條件。但應注意受政策因素、基礎條件等因素影響，農民人均純收入仍維持在一個相對較低的水平上，淶水縣農民人均年純收入最高水平未達4 000元，而易縣則未達3 000元。

表11 近京區域模型時期影響：各年份在總體常數項基礎上的變化程度

3 結果討論

（1）耕地面積變化與農民人均純收入之間均存在著長期的均衡關系。但不同區域在兩者相互影響的方向上存在差異，山區區域結果顯示兩者存在正相關關系，而環京和近京區域結果顯示兩者存在負相關關系。耕地數量變化對農民收入的正負向影響與農民的收入來源及收入結構關系密切。山區區域與另外兩個區域相比，更為貧困，農民長期以來保持著以糧食生產為主的收入結構，因此耕地面積的增加會促進農民收入的增加，兩者呈正向影響，但這也是農業勞動力轉移進程落后、非農收入增長緩慢的結果［17］。根據農業邊際收益論解釋，當農戶有機會選擇收益率更高的生產經營活動時，就可能減少或放棄耕種土地［18］。環京和近京區域區位優勢相對山區區域較為明顯，其非農就業機會相對較多，農業外部的就業機會增加和比較利益的驅動會促進農民向非農領域轉移。此時，因其收入結構的變化，耕地面積的減少并不會造成農民收入的下降。另外，由于經濟作物的比較利益，農業結構內部調整也會造成耕地數量變化與農民收入的負相關關系形成。因此，為減少耕地數量對農民收入的制約，應充分利用環京津的區位優勢，積極發展經濟，拓展農民的就業渠道，使農民能夠通過非農就業來提高收入水平。

（2）從各區域內部耕地數量與農民人均年純收入關系來看，在山區區域的3個縣中，阜平縣農民收入與耕地面積無關的部分最大，其次順平縣，淶源縣最小；在環京區域，曲陽縣相對唐縣農民收入與耕地面積無關的部分較大；近京區域，則表現為易縣相對淶水縣農民收入與耕地面積無關的部分較大。其農民收入與耕地面積的相關性大小除與區域的農民收入來源和收入結構有一定關系外，在一定程度上也可以體現出該區域的經濟發展狀況，相對發達的經濟條件會促進農業投資，提高單位面積產量，降低耕地面積的限制，從而使農民增收［19］。

（3）由于研究區3個不同區域自然條件稟賦和社會經濟條件的差異，使得農民收入的變化幅度差異較大。尤其表現在近京區域，由于其相對優越的自然、區位、經濟等條件，在耕地面積逐年減少的趨勢下，農民收入仍能保持相對穩定的增長，其截距項的變化印證了近京區域的農民收入增長優勢。而經濟條件相對較差的山區區域則表現出農民收入的增長乏力。這一研究對于如何協調耕地、農民收入與經濟發展的關系具有指導意義。

（4）根據對研究區的數據資料分析和實地調查，生態退耕和建設占用 （耕地非農化）是造成研究區耕地面積減少的主要原因。自2002年在該區域廣泛開展退耕還林工程，至2006年累計退耕面積2.44×104hm2。從近期來看，該區域過多承擔了維護和改善京津兩地生態環境的責任，但在生態補償力度上還很欠缺［20］，也是造成該區域經濟相對落后的原因之一。為此，京津兩地政府應當通過財政轉移支付、產業帶動等方式對該區域加以補償，以保證京津與環京津地區的協調發展。同時，加大對環京津貧困區在市場導向、資金融資等方面的政策扶持力度，以促進該區域的經濟發展，為農民提供就業渠道。

4 結論

山區區域、環京、近京區域3個區域的耕地面積與農民人均純收入之間均存在著長期的均衡關系。這意味著研究區域在相當長的一段時間內，農民人均純收入不可避免的要受到耕地面積變化的影響。

耕地面積變化對農民人均純收入的影響方向因區域而異。山區3縣分析模型顯示農民人均純收入與耕地面積呈正相關關系，耕地面積每增加1%，農民人均純收入增加0.362 6%。而環京區域和近京區域分析模型則顯示農民人均純收入與耕地面積呈負相關關系，耕地面積每增加1%，農民人均純收入分別下降0.501 7%和0.132 8%。耕地數量與農民收入的正負向關系主要與農民的收入來源及收入結構有關，同時區域的自然條件限制也是不可忽略的因素。

區域內部耕地數量與農民人均年純收入關系存在差異。耕地數量與農民收入的相關性大小在一定程度上體現出該區域的經濟發展水平。經濟發展條件較好的地區，農民收入受耕地面積限制較小，因此各區域應當積極尋求促進區域經濟增長的措施，提高單位面積耕地的投入產出，同時拓展農民就業渠道，從而保障農民收入的提高。

時期影響分析中，農民人均純收入均隨著時間的推移，截距項越來越大，說明農民人均年純收入在逐年增長。但增長幅度各區域間有明顯差異。山區區域表現出總體上升，個別年份有回落；環京區域表現為波動上升；而近京區域表現為直線上升，增長幅度明顯超過前兩個區域，這主要是由耕地資源稟賦和區位條件差異所致。

［1］曲福田，吳麗梅.經濟增長與耕地非農化的庫茲涅茨曲線假說及驗證［J］.資源科學，2004，26（5）：61－67.

［2］許廣月.耕地資源與經濟的增長關系：基于中國省級面板數據的實證分析［J］.中國農村經濟，2009（10）：21－30.

［3］李兆富，楊桂山.蘇州市近50年耕地資源變化過程與經濟發展關系研究［J］.資源科學，2005，27（4）：50－55.

［4］龍開勝，陳利根.耕地資源數量與經濟發展和城鄉收入關系的計量分析［J］.資源科學，2007，29（4）：139－145.

［5］郭貫成.耕地面積變化與經濟發展水平的相關分析：以江蘇十三個市為例［J］.長江流域資源與環境，2001，10（5）：440－447.

［6］葉浩，濮勵杰.江蘇省耕地面積變化與經濟增長的協整性與因果關系分析［J］.自然資源學報，2007，22（5）：766－774.

［7］范輝，陳超，余國忠.河南省耕地資源與城市化發展的協調性分析［J］.水土保持通報，2010，30（4）：184－189.

［8］鮑艷，胡振琪，柏玉，等.主成分聚類分析在土地利用生態安全評價中的應用［J］.農業工程學報，2006，22（8）：87－90.

［9］周嘉，張洪峰，尚金城，等.模糊綜合評判法在生態旅游戰略環境評價中的應用［J］.東北林業大學學報，2004，32（2）：52－54.

［10］陳華，王穩.中國保險發展對技術創新長短期影響效應的實證研究［J］.中國軟科學，2011（3）：82－87.

［11］Pedron I P.Critical values for cointegration tests in heterogeneous panels with multiple regressors［J］.Oxford Bulletin of Economics and Statistics，1999，61（S1）：653－678.

［12］尹希果.計量經濟學原理與操作［M］.重慶：重慶大學出版社，2009.

［13］郭軍華.中國城市化對城鄉收入差距的影響：基于東、中、西部面板數據的實證研究［J］.經濟問題探索，2009（12）：1－7.

［14］Luciano G.On the power of panel cointegration tests：A monte carlo comparison［J］.Economic Letters，2003，80（1）：105－111.

［15］胡莉莉，牛叔文，馬莉，等.基于面板數據模型的中國農業生產用能與農業經濟增長關系［J］.農業工程學報，2011，27（6）：1－6.

［16］丁洪建.基于面板數據的省會城市居住地價的差異及成因研究［J］.中國土地科學，2009，23（6）：14－19.

［17］陸益龍.農戶的耕地使用行為及其影響：基于2006CGSS的實證分析［J］.江蘇社會科學，2012（2）：110－117.

［18］姜天龍，郭慶.農戶收入結構支撐下的種糧積極性及可持續性分析：以吉林省為例［J］.農業經濟問題，2012（6）：14－20.

［19］張雯，張廣勝.農民增收背景下的現代農業發展途徑探討［J］.農業經濟，2012（1）：11－12.

［20］張蓬濤，張貴軍，崔海寧.基于退耕的環京津貧困地區生態補償標準研究［J］.中國水土保持，2011（6）：9－12.