中國奶業地理集中演化模式及其成因分析

道日娜，羅燕芳

（1.內蒙古大學經濟管理學院，內蒙古 呼和浩特 010021；2. 長江師范學院財經學院，重慶 涪陵 408100）

中國奶業地理集中演化模式及其成因分析

道日娜1,2，羅燕芳1

（1.內蒙古大學經濟管理學院，內蒙古 呼和浩特 010021；2. 長江師范學院財經學院，重慶 涪陵 408100）

本文利用1990-2015年中國奶業省際面板數據，借助相對熵值法考察了中國奶業空間分布特征，并引入動態面板模型實證了中國奶業空間分布的驅動因素。研究發現：1）1990年以來，中國奶業地理集中演化模式總體呈現“橫向擴散”特征；2）區域層面看，華北和西北奶業呈現“橫向擴散”模式，而大城市、東北和南方奶業凸顯內部擴大模式；3）2011年中國奶業首現衰退惡化跡象，但華北、大城市和南方表現為橫向擴散；4）從區域之間地理集中模式演化看，中國奶業消耗縮減模式存在區域之間梯度轉移現象，其轉移路徑依次為“大城市、華北、東北、西北”；而南方奶業表現出奶牛存欄持續增長的良好發展勢頭；5）計量模型結果表明，地區飼料供給，奶牛存欄一期滯后以及經濟增長促進奶業集聚，而奶牛單產和人口規模遞增驅使奶業擴散。

奶業地理集中模式；奶牛存欄；相對熵值；集聚；分散

Abstract：Based on a state-level panel data set of China’s dairy industry from 1990 to 2015, and applying the dynamic panel model and the relative entropy method, this paper investigated the spatial distribution characteristics of China’s dairy industry and analyzed the driving factors of it. Results show that: 1) since 1990, at the national level,the geographical concentration pattern of China’s dairy industry has been observed as “diffusion”. And at the regional level, both Central North China and Northwest China showed “diffusion” pattern. However, augmentation features were shown in metropolitan areas, Northeast and South China; 2) in 2011, the recession pattern first appeared at the national level, while the pattern of diffusion was maintained in the regions of Central North China, South China and metropolises; 3) in terms of the change of regional geographical concentration patterns, “attrition” pattern shifted orderly from “metropolitan areas to Central North, Northeast, and to Northwest China”. Surprisingly, South China has enjoyed a consistent increase in cow inventories all the time since 1990; and 4) according to the econometric analysis results, regional feed supply, the one year lag of cow inventory and economic growth factors promote dairy cows clustering, while per cow productivity and the total population boost lead to the dispersion of cows.

Key words：dairy geographical concentration pattern; cow inventory; relative entropy; clustering; dispersion

伴隨奶業規模化進程，中國奶業區域格局正在發生深刻變化。一方面，規?；a促進生產專業化和產業鏈縱向合作，降低交易成本，淘汰小規模散養戶，減弱以“小規模散養”為主的奶業產區比重。投入于規?；a的技術進步，突破了氣候、土地資源以及飼料供給等物質要素對奶業生產的先天制約，進而削弱傳統區域的比較優勢。另一方面，政府和乳企對奶業生產的大規模注資，擴充了奶業基礎，直接替代了奶業自我循環積累，助力新興奶業產區迅速崛起。經濟增長的區域差異，導致乳制品消費數量和結構的差異，吸引奶業向大的消費市場或臨近區域轉移。另外，人口密度增加催生政府環境管制，進而驅使奶業向環境管制寬松、人口稀疏區域轉移。綜上，規模經濟、技術進步、經濟增長以及人口密度等因素都會對奶業區域布局產生影響，促進奶業區域重組。

中國奶業區域分布不均衡、“北奶南運”問題早有學者關注[1-2]。然而直到2008年三聚氰胺事件之后，奶業區域分布變遷問題才引起學者們廣泛關注。學者們的研究主要從以下三個方面展開：一是總結中國奶業區域格局演化特征[3]；二是測算中國奶業不同區域的成本比較優勢[3-4]；三是利用生產集中度指標總結中國奶業區域格局演化特征基礎上，基于靜態面板數據模型，實證分析奶業區域格局演化的影響因素[5-8]。這些研究開拓了中國奶業區域分布變遷研究的視野，揭示了這一地理變遷的區域、省際特征及其驅動因素。但是上述研究還存在以下不足：首先，在統計指標選擇方面，基于奶牛存欄或牛奶產量的生產集中度指標，不能揭示長期以來中國奶業在總體層面的時空分布特征，更不能揭示區域內部和區域之間的非均衡演化趨勢。其次，在計量模型選擇方面，基于靜態面板模型的分析，不能克服奶牛存欄與其滯后變量的內生性問題，進而產生有偏差的結果。最后，在理論視角方面，上述研究只籠統地關注奶業區域格局演化的影響因素，并未清晰界定其集聚和分散機制，進而無法比較各個因素的作用方向及其強度對奶業地理分布的影響。

為了彌補上述研究的不足，本文首次將熵值法應用到中國奶業時空分布的研究，以期表示奶業地理分布的非均衡特征及其隨時間變化趨勢。同時，引用動態面板數據模型，從集聚和分散兩個方面，揭示中國奶業地理集中模式演化的影響因素。

1 研究方法和數據來源

1.1 研究方法

1.1.1 Theil的熵值法 在信息論中，熵是對不確定性的一種度量。信息量越大，不確定性就越小，熵也就越小。熵值法是一種客觀賦權法,它根據各項指標觀測值所提供信息量的大小確定指標權重。根據熵的特性，可通過計算熵值判斷一個方案的隨機性、無序性，也可判斷某個指標的離散度。Theil的熵值法由泰爾利用信息理論中的熵概念提出，主要應用于測量收入不平等。Herath等[9]將其作為一種測度行業集中度的指數，研究美國肉豬、肉牛和奶牛的空間和地理分布演化。

熵值有絕對熵值和相對熵值之分。若n代表數據集里中國省份個數，θit代表每一個省份每一年的奶牛存欄比重，則絕對熵值定義如下：

根據定義可知絕對熵值的取值范圍為0≤Н(θ)≤log2n，如果全國的奶牛存欄都集中到一個省，則該省θi=1，Н(θ)=0，表明中國奶業集中度最大或分散度最??；如果全國的奶牛存欄平均地分散到各個省，則每個省的θi值都相等，即Н(θ)=log2n，說明中國奶業集中度最小或分散度最大。絕對熵值能夠有效用來區別具有相同省份個數的區域之間在地理集中度方面隨時間的變化趨勢，但不能夠測度具有不同省份個數的區域之間在地理集中度方面隨時間的變化趨勢。鑒于中國奶業區域劃分，每個區域涵蓋的省份數量不盡相同，所以還需要引入相對熵值分析區域內部和區域之間的變化。相對熵值定義如下：R(θ)=1說明全國奶牛存欄平均地分散到每一個??；R(θ)=0比率趨近于0，說明全國奶牛存欄在地理上逐漸集中；R(θ)比率趨近于1，表明中國奶業區域分布不均衡，差異明顯。

1.1.2 動態面板數據模型方法 本研究旨在探索奶業區域格局演化的影響因素，所用數據包含時間和空間兩個維度，因此，在計量分析時采用動態面板數據模型（Dynamic Panel Data Model）。動態面板數據模型的一個突出優點是通過控制固定效應較好地克服了變量遺漏問題，而且還較好地克服了反向因果性問題[10]。由于解釋變量的內生性問題會導致參數估計的非一致性，故本文采用GMM估計法來估計參數。GMM的基本思想是先求出水平方程的一階差分形式：

式中：i=1, 2, …,N；t=1, 2, …,T，?yit是被解釋變量省際奶牛存欄的變化量，?yit-1是被解釋變量省際奶牛存欄的變化量的滯后期，?xit是解釋變量變化量的矩陣，i，t分別為截面維度和時間維度，a0，…，ak為待估參數矩陣，?εit為隨機擾動項。

對于動態面板數據模型，要增強參數估計結果的可靠性，還需要對模型設定的合理性和工具變量的有效性進行檢驗。一般來說，Sargan檢驗常用來檢驗工具變量的有效性，AR（2）統計量則用來檢驗原模型一階差分后是否存在二階自相關問題。此外，為了評價估計結果和滯后階的穩健性，還需要對估計的模型進行檢驗，常用方法主要是對面板殘差的穩定性進行檢驗。由于IPS檢驗不僅允許殘差項存在序列相關，而且也允許各剖面時間序列具有不同的動態行為和殘差方差，本文選用IPS檢驗來診斷面板殘差的平穩性，進而確認GMM估計是否是偽回歸的結果[10]。

1.2 數據來源

本研究應用的我國省際奶牛存欄和其他奶業相關數據來源于《中國奶業年鑒》（2001-2016），相關解釋變量數據來源于《中國統計年鑒》（2003-2016）、《中國農村統計年鑒》（2003-2016）、《新中國五十五年統計資料匯編》等。數據的時間維度為1990-2015年，空間維度為除重慶、西藏、青海和海南以外的中國大陸27個省市自治區。

2 結果與分析

2.1 基于熵值法的中國奶業地理集中模式分析

2.1.1 地理集中模式分類 地理集中度相對熵值測量方法的一個缺點是不能夠提供各個地理區域之間的奶牛存欄數目的絕對變化。Hubbel和Welsh[11]提出了一種將絕對存欄數量和地理集中度兩者結合起來表征地理集中模式的分類方法。Herath等[9]等將其應用到了奶業領域，具體分類見下表1。由表可得，根據集中度和奶牛存欄數的指標可以把奶業的地理集中模式分為內部擴大、重新分配、消耗縮減、橫向擴散和衰退惡化五種類型。例如，當集中度和奶牛存欄數都增加時對應的地理集中模式屬于內部擴大，其他四種分類以此類推。

表1 地理集中模式分類[9]Table 1 Geographic concentration pattern classification

2.1.2 中國奶業地理集中模式演化 根據中國農業部2010年《中國奶業發展規劃（2009-2013）》，中國大陸奶業版圖（不包括港澳臺）劃分為大城市、華北、東北、西北和南方五大產區。大城市包括北京、上海、天津和重慶四個直轄市；華北包括河北、河南、山西和山東四省；東北包括內蒙古、黑龍江、吉林、遼寧四省；西北包括新疆、陜西、寧夏和甘肅省；南方包括江蘇、浙江、安徽、福建、江西、湖北、湖南、廣東、廣西、海南、四川、貴州、云南和西藏14個省市和自治區?；谌珖膛４鏅诘南鄬刂岛痛鏅诳倲档淖兓闆r，自1990年以來中國奶業發展可以劃分為五個歷史時期，即1990至1999年、1999至2002年、2002至2006年、2006至2011年，2011至2015年。進一步從區域層面考慮，發現只有西北地區在1990至1999年間，其波動較全國奶業波動劇烈，其他區域奶業波動周期均小于全國奶業。因此，這一時段劃分標準，既能表現奶業波動的全局性特征，也能揭示奶業區域異質性特點。具體情況見表2。根據以上區域和時間段劃分，可以得到中國奶業的全國總體以及區域地理集中模式變化（表3和圖1）。

1）1990至2015年中國奶業地理分布變遷在全國層面呈現橫向擴散模式，即奶牛存欄持續增長，而奶業地理集中度在波動中降低（圖1）。可見，經過20多年的發展，中國奶業的地理集聚尚不明顯。區域層面，華北、西北奶業與全國奶業地理變遷模式一致，而大城市、東北和南方奶業則呈現內部擴大模式，即奶牛存欄增長的同時，奶業向區域內部省份集中。具體而言，大城市區域的上海奶牛存欄明顯減少，逐漸向北京和天津集中；東北區域的吉林奶牛存欄大幅減少，向黑龍江、遼寧以及內蒙古集中；南方區域的江蘇、浙江、湖北、廣東、四川、云南等省份的奶牛存欄減少，逐步向安徽、福建、江西、湖南、廣西和貴州等省份集中。

2）1990至1999年是中國奶業內部擴大期，奶牛存欄數持續增長，奶業趨于區域集聚。這一時期，東北和南方奶業都呈現內部擴大模式，華北地區則呈現橫向擴散模式，大城市奶業卻呈現消耗衰減趨勢，可見中國奶業存欄減少現象首先發生在大城市。究其原因，這一時期是中國奶業市場化改革時期，其主要內容是放開牛奶購銷價格、取消牛奶補貼、實行市場競爭、開放乳制品市場[12]。因此，奶業市場的逐步放開，激發了資源富集區域的成本優勢，替代了飼料資源短缺的大城市區域生產。與其他區域不同的是，西北奶業在這一時期呈現兩個階段性特征；首先是1990至1997年橫向擴散模式；其次是1998至1999年的內部擴大模式。即在區域奶牛存欄普遍增長的情況下奶牛先擴散后集中。其具體表現是1990至1997年，西北各省的奶牛存欄普遍增長，而1998年后，陜西、甘肅和寧夏奶牛存欄減少，特別是甘肅，在1997至1998年間，奶牛存欄由21.4萬頭驟降至5萬頭；寧夏也由7.1萬頭降至6.1萬頭，陜西由11.5萬頭降至11.2萬頭，只有新疆奶牛存欄保持增長勢頭，奶業進一步向新疆集中。

3）1999至2002年是中國奶業橫向擴散期，奶牛存欄數和相對熵值雙重遞增，奶業趨于擴散。區域層面，大城市、華北、西北和南方四大區域都呈現橫向擴散模式，與中國奶業的變化模式一致；而東北地區則呈現內部擴大的發展模式。這一時期也是中國奶業高速發展時期，乳品需求增加和乳品企業規模擴張帶動中國乳業進入高速發展階段。全國奶牛存欄數從1999年末的432萬頭增加到2002年末的667.3萬頭，增長了54.46%。

表2 1990-2015年中國奶業總體以及區域相對熵值變化情況Table 2 Changes in the degree of relative entropy of China’s dairy industry at the national and provincial level, 1990 to 2015

表3 1990-2015年全國以及區域奶業地理集中模式變化Table 3 Changes of the national and regional dairy geographic concentration patterns, 1990 to 2015

圖1 1990-2015年全國和區域奶業相對熵值和奶牛存欄趨勢圖Fig.1 Relative entropy and cow inventory trend of both national and regional levels from 1990-2015

4）2002至2006年，中國奶業再度迎來內部擴大期，奶牛存欄持續擴張，而相對熵值持續降低，奶業傾向于集聚。這一時期與1990至1999年中國奶業所呈現的模式是一致的，但是所對應的五大區域內部模式有變化。其中大城市、東北、西北和南方四大區域都呈現內部擴大模式，與中國奶業的變化模式一致；而華北地區則呈現橫向擴散的發展模式。究其原因，自2003年末始，中國奶業養殖環節生產成本遞增，銷售環節乳品企業過度競爭，導致利潤率下降，產業鏈“兩頭不見利”[13]導致奶業內部整合，生產集中度上升。

5）2006至2011年是中國奶業橫向擴散期，奶牛存欄和相對熵值都持續增長，奶業依然趨于擴散。2006年以后，中國奶業進入生產結構調整期，其顯著特征是小規模散養奶戶大量退出，奶業生產向100頭以上規模養殖場集中[14]。由于區域奶業的規模化比重不同，奶業生產結構調整在區域層面表現出差異性。其中，東北和西北兩個區域都呈現橫向擴散模式，而大城市和華北地區則呈現消耗衰減，南方依然延續了2002-2006年的內部擴大模式?？梢?，在奶業生產結構調整中，大城市和華北地區經歷了生產逐步集中和奶牛存欄逐漸減少的發展過程，而南方則迎來了存欄和生產集中度雙重遞增的發展階段。

6）2011年是中國奶業發展的一個重要轉折點，多年來持續增長的奶牛存欄開始下滑，全國奶業進入衰退惡化期。同時，與全國奶業變化趨勢相近的是奶業兩大主產區東北和西北，分別呈現存欄減少、集中度上升的消耗縮減模式和存欄減少、集中度降低的衰退惡化模式。與此相反的是，華北、南方和大城市表現出奶牛存欄增加，地理集中度降低的橫向擴散模式?？梢?，在奶業總體形勢不好的情況下，華北、南方和大城市奶業逆勢而上，表現出強勁的發展勢頭，證實中國奶業區域重組悄然開始。

縱觀中國奶業20多年的地理分布歷史變遷，在全國層面呈現橫向擴散模式。區域層面，華北和西北與全國奶業演化模式一致，而東北、大城市和南方奶業呈現內部擴大模式。2011年中國奶業出現衰退惡化跡象。與此相對應，東北和西北紛紛呈現消耗縮減和衰退惡化模式，而大城市、華北和南方奶業則表現為橫向擴散模式。同時，中國奶業地理分布變遷呈現五個階段性特點；全國總體層面，奶業地理集中模式呈現內部擴大、橫向擴散交替出現，直至2011年開始衰退。區域層面，兩大奶業主產區東北和華北分別以內部擴大和橫向擴散為主導模式，并先后出現了消耗縮減態勢。而中國奶業消耗縮減模式首現于大城市，且經歷過兩次?？偨Y中國奶業區域層面的時空演化，發現中國奶業消耗縮減模式存在區域之間梯度轉移現象，其轉移路徑依次為“大城市、華北、東北、西北”。特別值得關注的是，南方奶業在過去的20多年里，從未出現消耗或衰退跡象，受全國奶業波動影響最小。

2.2 中國奶業地理集中模式演化的決定因素及其作用機理

要素稟賦理論認為，產業的區位選擇由自然資源、地理位置，氣候等外生變量決定，經濟活動的空間分布取決于外生變量在空間上的分散或集中。奶業屬于自然資源轉化產業，其空間分布與區域的飼料資源、地理區位以及氣候條件密切相關，即地區的初始稟賦或第一自然（first nature）對產業分布有決定性影響。許多學者證明區域土地資源稟賦、飼料供給等因素促進奶業集聚[6-8]。在產業空間重組研究文獻中，新經濟地理理論得到了學者們普遍關注，特別是在中國制造業集聚領域得到了最廣泛的應用。突破傳統貿易理論的產業區位選擇外生性假設，新經濟地理學認為產業區位及其變化是一個內生過程。該理論認為根植于規模經濟的生產集中力量和運輸成本引起的生產分散力量的關系決定了產業區位[15]，其強調生產的收益遞增、集聚經濟和累積循環過程對經濟活動空間分布的影響。學者對中國奶業實證研究表明，中國奶業有顯著的路徑依賴特征，產業自我循環積累對產業區域集聚有重要影響，但運輸成本以及規模經濟因素對奶業區域分布影響不顯著[6,8]。

由上述奶業相對熵值分析可以看出，中國奶業仍處于橫向擴散期，產業分散因素對奶業地理分布有重要影響。根據新經濟地理理論，產業擴散是產業集聚的孿生兄弟，產業集聚到一定程度后導致集聚區域內非貿易品價格居高不下，地價上升、環境污染加劇等擁擠成本的產生，從而形成了推動產業向新的地理區位擴散的離心力[16]。奶業空間集聚對生態環境的污染增加，產生人畜矛盾，迫使政府出臺環境規制政策，促使產業轉移。對美國畜牧業研究表明，州政府的環境規制政策促使美國奶牛業、肉牛業以及生豬業等畜牧業由傳統區域向非傳統區域轉移[9]。虞袆等[17]對中國生豬產業的研究，也得到了相同的結論。具體到中國各地奶業，關于環境管制政策信息不完全，我們可以考慮利用人口密度部分替代環境管制政策，分析兩者之間的關系。進一步，對美國奶業的研究表明人口密度、技術進步以及市場力量是決定其空間分布的重要外生力量[18]。同時，道日娜和羅燕芳[8]對中國奶業的研究表明伴隨奶業生產技術進步，奶牛單產水平提升替代奶牛數量增加，尤其在奶業技術發達區域形成奶業驅散力，促使奶牛由高單產區域向低產區域轉移。中國大城市區域出現的消耗縮減模式即是這一地理變遷的證據。因此，基于奶牛單產水平提升的技術進步成為奶業區域發展的離心力。綜上，在新的奶業地理集中模式演化背景下，要素稟賦論和新經濟地理理論相互補充，共同識別奶業集聚與分散的影響因素。具體的模型變量選擇見表4。

2.2.1 被解釋變量 奶業區域格局變化可由區域內新增牛場數量或奶牛存欄比重表示。但是存欄比重變化不僅受本地奶牛存欄量變化的影響，還受到全國總體存欄變化的影響，因此不能直觀表達區域經濟總體行為。Bartik[19]認為利用奶牛存欄水平不僅可以反映區域經濟的總體行為，還可以反映經濟主體不同的經濟決策。現有奶牛存欄規模擴張或縮減、新奶牛場的進入或原有牛場的退出都會影響區域奶牛存欄。因此，本文選擇奶牛存欄數作為被解釋變量進行分析。

2.2.2 解釋變量 1）集聚因素。飼料供給。奶業是自然資源轉化產業，奶牛飼養需要玉米、豆粕以及秸稈等粗飼料和青飼料的投入。事實上，中國奶業區域分布與玉米帶分布天然重合，即表明了這一種養業之間的關系。飼料種植結構決定奶業的物質資料基礎。許多學者的研究表明自然資源是影響奶業區域格局演化的重要因素[6-7]。本文研究選擇省際年玉米產量和年大豆產量作為飼料供給的決定變量，分析飼料種植與奶牛養殖之間的關系。

奶牛存欄滯后變量。由資源優勢決定的奶業主產區通過奶牛的周期性繁殖，形成自我強化和循環積累，推動奶業區域由“第一”優勢（資源）區域向“第二”優勢（資本）區域轉化。有研究表明，上一期的奶牛存欄促進奶業地理集聚[6-7]，也有研究表明奶牛兩年的生產周期是影響牛奶供給的重要因素[20-21]。因此，本文選擇奶牛存欄滯后一期和二期兩個變量，研究奶業自我循環積累機制對地理分布的影響，進而說明中國奶業是否突破了“自我循環”的路徑依賴，轉向新的區域。

表4 模型變量選擇Table 4 Model variable selection

市場潛力。在新經濟地理學文獻中，市場潛力通常用來衡量一個地區市場規模的大小，我們引入這一變量來反映本地市場規模對奶業集聚的影響。有研究表明中國持續的經濟增長促進了奶業發展[22]。因此，我們用省際人均GDP來表示本地市場潛能對地方奶業發展的影響，且這一影響有集聚和分散的兩種可能。中國奶業發展歷史較短，人們的飲奶保健習慣剛剛形成，故奶業從自給自足到商業化發展過程中，經濟增長推動了乳品需求，對地方奶業發展有一定的助推作用。但是隨著經濟快速發展，地價攀升，勞動力工資上漲，其他邊際報酬高的產業對奶業構成擠壓，會驅散奶業在經濟發達地區的集中。經濟增長在不同階段可能起到集聚和分散兩種力量的不同作用，對奶業空間分布產生影響。

2）分散因素。人口規模。人口數量是影響畜牧業產業布局的重要因素。人口規模遞增對奶業布局的影響機制主要表現為兩個方面：一是人口增加推動土地價格上漲，進而增加奶業生產成本，驅散奶業集聚；二是人口規模遞增促使政府加強環境管制，驅使奶業向環境管制寬松地區轉移（污染天堂效應）。美國奶業由東北部人口密集區向西南人口稀疏地區轉移恰好遵循這一邏輯[23]。本文選擇省際人口數量分析地方人口規模變化對奶業區域格局的影響。

奶牛單產。縱觀發達國家奶業現代化歷程，均以奶牛單產水平提升伴隨奶牛存欄減少為特征。因此，奶業發展到一定程度，奶牛單產水平提升替代奶牛存欄數量增長，推動奶業由數量型增長模式向質量型增長模式轉化。同時，奶牛單產水平的空間差異導致奶牛生產率的地區差異，進而導致奶牛由高生產率區域向低生產率區域轉移（大城市向周邊地區），且只要這種空間生產率差異存在，區域之間奶牛遷移就會發生。因此，奶牛單產作為奶業技術進步的重要代理變量，對區域奶牛存欄增長構成一股重要的驅散力量，驅趕奶牛向低產區轉移。

2.3 模型設定與變量描述

根據上述區域奶牛存欄變化的影響因素，模型設定為雙對數線性的計量模型，具體模型如下；

式中：a0是常數項 ，a1、a2、…、a7為待估參數，εit為誤差項。式中，logCowit是被解釋變量，代表各省年奶牛存欄數的對數值；logCowit-1，logCowit-2分別為各省奶牛存欄量的滯后一期和二期對數值。logTopuit為省總人口數量、logCornit為省際玉米年產量、logBeanit為省際大豆年產量、logAgdpit為省際人均GDP、logProdit為省際奶牛單產水平。所有變量都取對數值，以消除可能的異方差。

模型變量的總體性描述分析見表5。

2.4 回歸分析

2.4.1 計量分析 就動態面板數據模型而言，GMM方法又可分為DIF-GMM(差分GMM)和SYS-GMM(系統GMM)兩類。利用STATA12.0軟件，我們分別得到差分GMM和系統GMM兩組模型估計結果，DIF-GMM估計量較易受弱工具變量的影響而產生有限樣本偏誤，而SYS-GMM估計量則結合了差分方程和水平方程，具有更好的有限樣本性質。從估計方法的結果以及相關的檢驗進行對比，系統GMM的Sargan檢驗P值和IPS檢驗P值都優于差分GMM；并且系統GMM的AR（2）檢驗結果說明一階差分后的殘差不存在二階自相關，而在差分GMM中AR（2）檢驗結果拒絕原假設，也即說明擾動項存在自相關。綜合以上的檢驗分析，因此我們以動態面板系統GMM模型的估計結果為依據，分析各因素對奶牛存欄數量的影響（表6）。

表5 變量的總體性描述Table 5 Descriptive statistics of variables

表6 系統GMM模型的回歸結果Table 6 Estimation results of the system GMM model

從表6中可以看出：模型系數聯合顯著性的檢驗值都在1%的水平上顯著；Sargan檢驗結果說明工具變量有效；AR（2）檢驗結果說明一階差分后的殘差是不存在二階自相關的，即本文所設定的模型是合理的；IPS檢驗的p值均小于0.01，驗證面板殘差是平穩的，進而確認本文所選取的系統GMM的估計結果可靠。本文補充多重共線性檢驗，結果表明總體解釋變量方差膨脹因子均值為7.63，小于10是可以接受的。其主要原因是奶牛存欄滯后一期與二期存在共線性，但是由于我們使用的動態面板數據模型，是對變量進行差分分析。同時從回歸結果看二期滯后的存欄數并不顯著，而且引入二期滯后的回歸結果表明并沒有影響到一期滯后的顯著性檢驗，這種共線性是可以接受的。

2.4.2 計量結果 根據表6的結果可以得到：所有的解釋變量均顯著。

1）地區奶牛存欄與大豆年產量顯著正相關，且2011至2015年，大豆在奶業生產中的貢獻不斷增強。這表明區域大豆產量增加促進奶業集聚。該研究結論與李鴻幫等[7]的研究結論一致，也揭示出區域飼料資源稟賦在奶業發展中的基礎作用[4]。

2）地區奶牛存欄與玉米產量負相關，但2011至2015年玉米對奶牛存欄的影響由負轉正。首先，地區奶牛存欄與玉米產量負相關的原因有：玉米產量的大幅度提高受到需求快速增長導致的價格上升的影響，而這部分需求當中用作飼料的占比不高，工業用玉米和其他用途的玉米占有很大的比重。其次，玉米產量受到補貼政策、糧食政策以及農村種植業結構調整因素的影響。然而奶牛存欄數的減少主要受到奶業規?；M程中小規模奶戶退出以及奶粉進口替代的影響。進一步，玉米的供給對散戶養殖的影響要遠比規?；B殖的影響大，而散戶的大規模退出主要是養殖戶與加工企業之間的利益分配矛盾。所以單純從數據的回歸結果上表現出玉米產量與奶牛養殖數量之間存在負相關關系。其次，2011年以后，玉米產量與地區奶牛存欄正相關的原因是，自2011年以后，全國大部分省市奶牛存欄數量與玉米產量集體遞增，可見玉米生產與奶牛存欄變化趨勢一致，玉米成為促進奶牛存欄增加的重要因素。

3）即期奶牛存欄與滯后一期奶牛存欄正相關，而與滯后二期關系不顯著。這一結果表明，中國奶業已經突破原有的路徑依賴模式，由傳統產區向新興區域轉移。西北奶業的衰退以及華北奶業的崛起是這一結果的良好印證。從模型1至模型2，滯后一期的影響在減弱，且滯后二期由負相關變為不相關，表明2011至2015年是中國奶業生產結構調整活躍期，表現為小規模散養戶不斷退出，而大規模養殖場不斷進入。中國奶業統計數據表明，這一時期，1至4頭小規模散養奶戶的奶牛存欄比重由23.99%降至20.3%，而100頭以上規模養殖戶的存欄比重由32.87%增至48.3%，其中特別是1 000頭以上大規模牧場的存欄比重由12.06%增長至23.6%，年均增長率達18.27%。而2008年，1 000頭以上養殖場的存欄比重僅為5.54%。新進入者以資本優勢擴充產業基礎，突破了奶業自我循環積累的歷史發展模式，奶業由傳統區域向新興區域轉移。

4）人均GDP增長促進地區奶牛規模擴張。模型結果表明人均GDP每增加1%，奶牛存欄數增加0.162%，即人均GDP每增長10 000元，奶?？傤^數將增加162頭。經濟增長對奶業的拉動作用主要表現為收入水平的提升帶動乳品消費的增長，進而推動奶業擴張。這與奶業總體長期發展趨勢相吻合，中國奶業迅速崛起與中國經濟增長拉動的消費增加密切相關[20]。

5）總人口對奶業的擠出效應，2011年以后變得顯著且增強。該結果表明，2011年以前，地方奶業環境管制政策不完善，奶業發展尚未實現環境協調發展模式。另外，在奶業規模化發展過程中，大量資本注入沖淡了經濟增長和人口遞增帶來的土地價格上升對奶業的驅散效應。但2011年以后，伴隨奶業規?；嵘龑Νh境要求從嚴，大城市等人口密集區域奶牛存欄大幅減少，表明人口規模對奶業的擠出效應開始顯現。這一變化與美國奶業的區域格局演化趨勢一致[9]。

6）奶牛的單產水平與地區奶牛存欄顯著負相關，且2011至2015年提升奶牛單產對減少奶牛頭數的貢獻在增強。 由模型2的結果可知，省際奶牛單產水平每增加1%，其奶牛存欄數減少0.191%，即奶牛單產水平每增長1 000 kg，省際奶?？傤^數將減少191頭。這一系數表明中國奶業依然處于數量型增長時期，奶牛單產水平的提升將會減少奶牛存欄量。提高奶牛單產是減少地區奶牛存欄，促進奶業區域轉移的重要手段。

3 結論及建議

3.1 結論

過去20多年里，中國奶業地理集中演化模式在全國總體層面呈現橫向擴散特征。區域層面，華北和西北奶業地理集中模式與全國一致，表現為橫向擴散模式，而東北、大城市和南方則呈現為內部擴大模式。2011年是中國奶業發展的重要轉折點，持續增長的奶牛存欄開始下滑，奶業出現衰退惡化跡象。與此相對應，傳統奶業產區東北和西北紛紛呈現消耗縮減和衰退惡化模式，而大城市、華北和南方奶業則表現為橫向擴散模式。總結中國奶業區域層面的時空演化，發現中國奶業消耗縮減模式存在區域之間梯度轉移現象，其轉移路徑依次為“大城市、華北、東北、西北”。令人驚奇的是，南方奶業一路凱歌，在過去的20多年里尚未出現縮減或衰退跡象。

本文計量模型結果表明，區域飼料供給、奶牛存欄滯后一期以及經濟增長促進奶業地理集聚，而奶牛單產和人口數量遞增促使奶業擴散?？梢?，中國奶業地理集中演化模式呈“市場主導，資源為輔”的驅動特征。人口規模遞增對奶業的擠出效應顯著，表明奶業發展過程中的人蓄矛盾或者環境管制因素對奶業發展的約束作用開始顯現。

3.2 建議

在新一輪的奶業區域格局演化過程中，政府應該充分考慮“市場、資源”等集聚因素和“環境、技術”等分散因素對奶業發展的作用，根據五大區域不同特征，引導奶業向資源、環境協調互動方向發展。具體而言，1）對于具有市場和技術優勢，而飼料供給和人口環境壓力大的大城市，應該“以提高單產”為首要目標，同時嚴格環境規制，逐漸減少奶牛存欄；2）東北傳統奶業產區，增大豆減玉米，增加奶牛存欄；3）對于人口稀疏的西北產區，政府引導乳企和大規模牧場布局，并增加飼草料播種面積，增加奶牛存欄，以實現奶業與資源環境的協調；4）鄰近市場，產業基礎好的華北區域，以增加區域大豆種植，并促進飼料外購的方式提升區域奶牛存欄。同時區域經濟增長會有效拉動奶業；5）人口多，飼料短缺但一直表現良好的南方奶業區域，鼓勵發展單產高，符合地方環境氣候條件的奶牛品種，以部分滿足地方市場。同時加強環境規制，以減輕人蓄矛盾。最后，當前“市場主導，資源為輔”的奶業區域分布模式中，政府應該充分考慮奶業發展過程中的環境成本，逐漸引導奶業向人口稀疏、資源環境壓力小的西北區域轉移。

[1]錢貴霞, 解晶. 中國原料奶供求矛盾關系的實證分析[J].農業現代化研究, 2010, 31(3): 272-276.Qian G X, Xie J. Contradictory relationship between raw milk supply and demand in China[J]. Research of Agricultural Modernization, 2010, 31(3): 272-276.

[2]花俊國. 我國奶業生產的波動性及區域特征——基于H-P濾波法的分析[J]. 西北農林科技大學學報(社會科學版),2014, 14(2): 60-67 Hua J G. The volatility and regional characteristics of milk production in China: Based on the analysis of H-P filtering method[J]. Journal of Northwest Agriculture and Forestry University (Social Science Edition), 2014, 14(2): 60-67.

[3]郭利亞, 張娜, 王加啟. 轉型升級期我國奶業生產區域分布特征變化分析[J]. 中國畜牧雜志, 2014, 50(8): 16-20.Guo L Y, Zhang N, Wang J Q. Analysis of varration characteristics of Chinese dairy production during the period of transformation and upgrading[J].Chinese Journal of Animal Science, 2014, 5(8): 16-20.

[4]于海龍, 李秉龍. 中國奶牛養殖的區域優勢分析與對策[J].農業現代化研究, 2012, 33(2): 150-154.Yu H L, Li B L. Analysis and policy proposals on regional advantages of dairy cows production in China[J]. Research of Agricultural Modernization,2012, 33(2): 150-154.

[5]張曼玉, 李彤, 褚林然. 不同養殖規模及主產省間原料奶成本收益分析[J]. 中國畜牧雜志, 2015, 51(24): 54-59.Zhang M Y, Li T, Zhu L R. Analysis on the cost benefit of raw milk production on different scales and main producing provinces[J]. Chinese Journal of Animal Science, 2015, 51(24): 54-59.

[6]衛龍寶, 張菲. 我國奶牛養殖布局變遷及其影響因素研究——基于我國省級面板數據的分析[J]. 中國畜牧雜志,2012, 48(18): 52-56.Wei L B, Zhang F. The distribution changes of dairy cows in China and its influencing factors: Based on the analysis of provincial panel data[J]. Chinese Journal of Animal Science, 2012, 48(18): 52-56.

[7]李幫鴻, 李翠霞, 鄒玉友. 中國原料乳生產區域布局變化及其影響因素研究——基于1995-2011年省際面板數據的分析[J]. 中國畜牧雜志, 2014, 50(8): 51-56.Li B H, Li C X, Zou Y Y. Chinese regional distribution of raw milk production and its influencing factors:Based on 1995-2011 provincial panel data analysis[J].Chinese Journal of Animal Science, 2014, 50(8): 51-56.

[8]道日娜, 羅燕芳. 從資源到資本: 中國奶業區域格局演化成因分析[J]. 農業現代化研究, 2016, 37(2): 205-213.Dao R N, Luo Y F. From resource to capital: Causal analysis of regional structural change in China’s dairy industry[J]. Research of Agricultural Modernization, 2016, 37(2): 205-213.

[9]Herath P, Weersink A J, Carpentier C L. Spatial and temporal changes in the U.S. hog, dairy and fedcattle sectors, 1975-2000[J]. Review of Agricultural Economics, 2004, 27(1): 49-69.

[10]肖干, 徐鯤. 農村金融發展對農業科技進步貢獻率的影響—基于省級動態面板數據模型的實證研究[J]. 農業技術經濟, 2012(8): 87-97.Xiao G, Xu K. Study on of the influnce of rural financial development on the contribution rate of agricultural scientific and technological progress:Based on provincial dynamic panel data model[J].Journal of Agrotechnical Economics, 2012(8): 87-97.

[11]Hubbell B J, Welsh R. An examination of trends in geographic concentration in U.S. hog production 1974-1996[J]. Agricultural and Applied Economics,1998, 30(2): 285-299.

[12]中 國奶 業協 會. 中 國奶 業改 革開放30年[J]. 中 國 奶 牛, 2009(1):Ⅰ-Ⅶ.China Dairy Association. Tracing the past 30 years of China’s dairy industry development after introuducing the reform and opening up policy[J]. China Dairy Cows, 2009(1): Ⅰ-Ⅶ.

[13]李建平, 劉冬梅. 我國奶業發展面臨的新問題與對策建議——對“奶農與乳品加工企業兩頭不見利”現象的思考[J]. 農業經濟問題, 2005(12): 69-72.Li J P, Liu D M. Development of Chinese dairy industry: Problems and solutions[J]. Issues of Agricultural Economy, 2005(12): 69-72.

[14]羅燕芳, 道日娜. 基于馬爾科夫模型的中國牛奶生產結構變遷分析[J]. 前沿, 2015(12): 74-78.Luo Y F, Dao R N. Using Markov Models to analyze the size distribution of China’s dairy farms[J]. Forward Position, 2005(12): 74-78.

[15]Krugman P R. Increasing returns and economic geography[J]. Journal of Political Economy, 1991,99(3): 483-499.

[16]Klmienko M. Comptition, matching, and geograghical clustering at early stages of the industry life cycle[J]. Journal of Economics and Business, 2004, 56(3): 178-195.

[17]虞袆, 張暉, 胡浩. 環境規制對中國生豬生產布局的影響分析[J]. 中國農村經濟, 2011(8): 81-88.Yu W, Zhang H, Hu H. Environmental regulation and pig production distribution in China[J]. Chinese Rural Economy, 2011(8): 81-88.

[18]Abdalla C W, Lanyon L E, Hallberg M C. What we know about historical trends in firm location decisions and regional shifts: Policy issues for an industrializing animal sector[J]. American Journal of Agricultural Economics, 1995, 77(5): 1229-1236.

[19]Bartik T J. The effect of environmental regulation on business location in the United States[J]. Growth and Change, 1988, 19(3): 22-44.

[20]烏云花, 黃季焜, ScottRozelle, 等. 農戶奶牛養殖

Analysis the pattern and causes of geographical concentration change in China’s dairy industry

DAO Ri-na1,2, LUO Yan-fang1
（1. School of Economics and Management, Inner Mongolia University, Hohhot, Inner Mongolia 010021, China;2. School of Finance and Economics, Yangtze Normal University, Fuling, Chongqing 408100, China）

F326.3

A

1000-0275（2017）05-0852-10

國家自然科學基金（71563036）；內蒙古自然科學基金（2014MS0702）;長江師范學院高層次人才引進科研啟動項目（2017KYQD05）。

道日娜（1980-），女（蒙古族），內蒙古赤峰人，管理學博士，講師，主要從事農業經濟方面的研究，E-mail：daorina345@163.com；羅燕芳（1989-），女，內蒙古烏蘭察布人，碩士研究生，主要從事乳業經濟方面的研究，E-mail：xinyufeiyang668@126.com。

2016-12-29，接受日期：2017-08-30

Foundation item:National Natural Science Foundation of China (71563036); Natural Science Foundation of Inner Mongolia (2014MS0702); The high-level introduction of talent research start-up fund of Yangtze Normal University (2017KYQD05).

Corresponding author:DAO Ri-na, E-mail:daorina345@163.com.

Received29 December, 2016;Accepted30 August, 2017

10.13872/j.1000-0275.2017.0088

道日娜, 羅燕芳. 中國奶業地理集中演化模式及其成因分析[J]. 農業現代化研究, 2017, 38(5): 852-861.

Dao R N, Luo Y F. Analysis the pattern and causes of geographical concentration change in China’s dairy industry[J]. Research of Agricultural Modernization, 2017, 38(5): 852-861.