大豆單產(chǎn)水平的時(shí)空演變及收斂特征
——基于東北大豆產(chǎn)區(qū)141個(gè)縣域面板數(shù)據(jù)

辛 戈，張德碩，林德榮

（青島農(nóng)業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院(合作社學(xué)院)，山東 青島 266109）

0 引言

糧食生產(chǎn)能力是糧食生產(chǎn)的核心能力，提升糧食生產(chǎn)能力是實(shí)現(xiàn)國(guó)家糧食安全的關(guān)鍵［1］。大豆作為重要的飼料糧和油料，在滿足居民消費(fèi)升級(jí)需要、保障新時(shí)代國(guó)家糧油安全等方面具有不可替代的作用。近年來，我國(guó)大豆進(jìn)口量持續(xù)攀升，自給率不足20%，對(duì)外依存度高，國(guó)產(chǎn)大豆缺乏市場(chǎng)競(jìng)爭(zhēng)力。導(dǎo)致我國(guó)大豆自給能力不足的重要原因是我國(guó)的大豆單產(chǎn)水平低。數(shù)據(jù)顯示，2020年世界大豆平均單產(chǎn)量為185 kg/667 m2，而我國(guó)大豆單產(chǎn)量?jī)H為132 kg/667 m2，僅為世界平均水平的71%左右。《大豆振興計(jì)劃實(shí)施方案》明確提出推動(dòng)我國(guó)大豆實(shí)現(xiàn)“擴(kuò)面、增產(chǎn)、提效、綠色”的目標(biāo)。2023年中央一號(hào)文件指出深入推進(jìn)大豆和油料產(chǎn)能提升工程。因此，提高大豆單產(chǎn)水平是提高大豆產(chǎn)量和增加自給率的關(guān)鍵路徑。東北地區(qū)是我國(guó)重要的大豆產(chǎn)區(qū)，2021年東北4省（區(qū)）大豆種植面積為513.76萬hm2，約占全國(guó)大豆種植面積的61%；大豆產(chǎn)量達(dá)967.13萬t，約占全國(guó)大豆總產(chǎn)量的59%。由于大豆生產(chǎn)受區(qū)域資源限制性較大，此外東北地區(qū)各省的土地、農(nóng)業(yè)勞動(dòng)力、生產(chǎn)要素等空間分布呈現(xiàn)出明顯的區(qū)域異質(zhì)性，而東北大豆主產(chǎn)區(qū)的大豆單產(chǎn)水平在結(jié)構(gòu)和區(qū)域上的差異特征鮮有報(bào)道。因此，客觀測(cè)度大豆單產(chǎn)水平在東北地區(qū)的現(xiàn)狀和差異特征，探尋提升大豆單產(chǎn)水平的實(shí)現(xiàn)路徑，對(duì)于保障國(guó)家糧食安全和大豆產(chǎn)業(yè)振興具有重要的現(xiàn)實(shí)意義。

1 文獻(xiàn)綜述

近年來，隨著社會(huì)各界對(duì)大豆生產(chǎn)能力的關(guān)注和重視，學(xué)者們針對(duì)大豆生產(chǎn)能力發(fā)展相關(guān)問題開展了卓有成效的研究，主要集中于以下3個(gè)方面。一是關(guān)于大豆產(chǎn)量及影響因素研究：胡國(guó)玉等［2-4］認(rèn)為氣候變化、土壤狀況的不同，會(huì)造成大豆產(chǎn)量差距變大，導(dǎo)致大豆單產(chǎn)水平不一；張曉峰等［5］通過建立GAEZ模型并結(jié)合近50 a中國(guó)氣候變化發(fā)現(xiàn)，中國(guó)大豆生產(chǎn)潛力呈現(xiàn)出由南向北、由西向東增加的趨勢(shì)；朱文博等［6-11］發(fā)現(xiàn)大豆種植面積、每公頃化肥 (折純) 施用量、農(nóng)業(yè)科技水平、灌溉方式等因素都會(huì)影響大豆生產(chǎn)力。二是關(guān)于大豆種植面積的研究：秦婷婷等［12-17］通過研究發(fā)現(xiàn)滯后一期的大豆播種面積、農(nóng)業(yè)比較收益、大豆相關(guān)政策等對(duì)當(dāng)期大豆播種面積都有顯著影響。三是關(guān)于大豆生產(chǎn)效率的研究：李柯逾等［18］通過調(diào)研數(shù)據(jù)發(fā)現(xiàn)，黑龍江省嫩江縣家庭農(nóng)場(chǎng)的大豆生產(chǎn)效率高于合作社、專業(yè)大戶和普通農(nóng)戶；王亞君等［19］通過建立DEA模型測(cè)算了大豆主產(chǎn)區(qū)10個(gè)省份的大豆生產(chǎn)效率，發(fā)現(xiàn)黑龍江、吉林、遼寧和山東等4省的綜合生產(chǎn)效率最高；夏佳佳等［20］測(cè)算了10個(gè)大豆主產(chǎn)區(qū)的全要素生產(chǎn)率，結(jié)果表明東北地區(qū)、西部地區(qū)和黃淮地區(qū)的大豆全要素生產(chǎn)率增長(zhǎng)表現(xiàn)出依次遞減的趨勢(shì)；郭天寶等［21］研究發(fā)現(xiàn)，近年來東北大豆生產(chǎn)的全要素生產(chǎn)率變動(dòng)呈“W”型的波動(dòng)變化趨勢(shì)。

現(xiàn)有研究為本文提供了借鑒和思路，但當(dāng)前研究尺度多為省域與地市視角，在構(gòu)建復(fù)合指標(biāo)體系綜合評(píng)價(jià)大豆生產(chǎn)水平過程中仍存在精準(zhǔn)度不夠的問題。中國(guó)糧食增產(chǎn)依靠的是單產(chǎn)水平的提高［22］。大豆作為我國(guó)重要的糧食作物之一，其單產(chǎn)水平一直備受學(xué)者的關(guān)注。但已有針對(duì)大豆單產(chǎn)的研究，大多集中于外部條件，少有文獻(xiàn)從縣級(jí)層面出發(fā)，探究不同空間尺度下大豆單產(chǎn)水平波動(dòng)的動(dòng)態(tài)時(shí)空特征。基于此，本文運(yùn)用Kernel密度估計(jì)、傳統(tǒng)Markov鏈和空間Markov鏈及收斂模型等方法分析2013—2021年東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平的現(xiàn)狀、時(shí)空演進(jìn)規(guī)律和收斂性特征，以期為保障國(guó)家糧食安全和促進(jìn)大豆產(chǎn)業(yè)可持續(xù)發(fā)展提供參考。

2 研究方法和數(shù)據(jù)來源

2.1 模型構(gòu)建

2.1.1 Markov鏈分析 Markov鏈?zhǔn)且环N時(shí)間和狀態(tài)均為離散的Markov過程。通過構(gòu)建Markov概率轉(zhuǎn)移矩陣，可以分析東北地區(qū)大豆單產(chǎn)水平向更高水平或者向更低水平類型轉(zhuǎn)移的可能性，從而反映出大豆單產(chǎn)水平的變遷規(guī)律。Markov鏈分析方法是一個(gè)隨機(jī)的過程{x(t),t∈T}，令隨機(jī)變量Xt=j，即t年j的狀態(tài)，滿足式（1）：

式（2）中：假設(shè)Pij為某地區(qū)大豆單產(chǎn)水平從t年第i種狀態(tài)轉(zhuǎn)移到t+1年第j種狀態(tài)的概率，ni為樣本期內(nèi)第i種大豆單產(chǎn)水平所出現(xiàn)的總次數(shù)。Pij可以通過極大似然估計(jì)法求得。

2.1.2 空間Markov鏈分析法 空間Markov鏈分析法除了具有與傳統(tǒng)的Markov鏈分析法共性特征外，其所反映出的規(guī)律性特征還可以考察各地區(qū)大豆單產(chǎn)水平動(dòng)態(tài)演進(jìn)過程中空間因素的作用，即相鄰地區(qū)的大豆單產(chǎn)水平對(duì)本地區(qū)大豆單產(chǎn)水平狀態(tài)轉(zhuǎn)移的影響。

2.1.3 收斂模型 收斂模型主要作用是考察落后區(qū)域是否能對(duì)較發(fā)達(dá)區(qū)域?qū)崿F(xiàn)追趕，近年來收斂模型逐漸應(yīng)用于分析地區(qū)差異和演進(jìn)趨勢(shì)的研究中［23］。收斂模型包括σ收斂和β收斂2種傳統(tǒng)收斂類型，β收斂又分為絕對(duì)β收斂和條件β收斂，本文采用上述3種收斂模型探究東北地區(qū)大豆單產(chǎn)水平的斂散特征。

2.1.3.1σ收斂模型 本文采用變異系數(shù)法考察σ收斂特征，建立如下σ收斂檢驗(yàn)公式：

2.1.3.2 絕對(duì)β收斂模型 絕對(duì)β收斂指在其他條件不變的情況下，隨著時(shí)序發(fā)展各縣(市、區(qū)、旗)逐漸收斂至同一水平，其檢驗(yàn)公式為：

2.1.3.3 條件β收斂模型 條件β收斂是指在控制各影響因素后，東北大豆產(chǎn)區(qū)的大豆單產(chǎn)水平將最終收斂至各自的穩(wěn)態(tài)水平。本文將大豆單產(chǎn)的影響因素納入到模型中，借鑒現(xiàn)有研究［9］，選取了如下控制變量：農(nóng)業(yè)地區(qū)產(chǎn)業(yè)結(jié)構(gòu)、單位面積機(jī)械投入量、單位面積化肥使用量、人均GDP。檢驗(yàn)公式為：

式（5）中：C是所加入的各控制變量的合集，γ是其回歸系數(shù)，其他變量含義同上。若β值顯著為負(fù)，則表明東北地區(qū)大豆單產(chǎn)水平存在條件β收斂態(tài)勢(shì)。

2.2 研究區(qū)域界定與數(shù)據(jù)來源

2.2.1 區(qū)域界定 東北大豆產(chǎn)區(qū)的大豆種植主要集中在內(nèi)蒙古自治區(qū)東部和黑龍江省第三積溫帶以北［24］。基于研究完整性和區(qū)域發(fā)展的協(xié)調(diào)性，以大豆生產(chǎn)直接關(guān)聯(lián)性為原則，借鑒已有學(xué)者的研究成果［25］，最終選定黑龍江省(67個(gè)縣)、吉林省(41個(gè)縣)、遼寧省(19個(gè)縣)以及內(nèi)蒙古自治區(qū)(14個(gè)縣)共141個(gè)縣為研究區(qū)域。具體研究區(qū)域如圖1所示。

圖1 東北大豆產(chǎn)區(qū)研究區(qū)域示意圖

2.2.2 數(shù)據(jù)來源 本文以大豆單產(chǎn)水平（t/hm2)作為研究對(duì)象。數(shù)據(jù)來源于2013—2021年各省、各地市的統(tǒng)計(jì)年鑒和縣域國(guó)民經(jīng)濟(jì)統(tǒng)計(jì)公報(bào)。其中黑龍江省和吉林省的縣域數(shù)據(jù)來源于省級(jí)統(tǒng)計(jì)年鑒，遼寧省和內(nèi)蒙古自治區(qū)的省級(jí)統(tǒng)計(jì)年鑒未報(bào)告縣域數(shù)據(jù)，相關(guān)數(shù)據(jù)從地級(jí)市統(tǒng)計(jì)年鑒獲取，個(gè)別缺失數(shù)據(jù)采用線性插值法進(jìn)行補(bǔ)全。

3 實(shí)證分析

3.1 東北地區(qū)大豆單產(chǎn)時(shí)序特征分析

根據(jù)統(tǒng)計(jì)數(shù)據(jù)，測(cè)算了2013—2021年東北大豆產(chǎn)區(qū)141個(gè)縣的大豆單產(chǎn)水平，測(cè)算結(jié)果如圖2所示。結(jié)果表明，東北大豆產(chǎn)區(qū)全樣本黑龍江、吉林、遼寧及內(nèi)蒙古的大豆單產(chǎn)水平在2013—2016年呈現(xiàn)平穩(wěn)波動(dòng)的態(tài)勢(shì)；2016—2017年，東北大豆產(chǎn)區(qū)全樣本、黑龍江、吉林、遼寧、內(nèi)蒙古大豆單產(chǎn)水平呈現(xiàn)出快速上升的趨勢(shì)，2017—2018年，東北大豆產(chǎn)區(qū)全樣本、黑龍江、吉林、遼寧與內(nèi)蒙古大豆單產(chǎn)水平的變化趨勢(shì)開始出現(xiàn)分化，黑龍江、遼寧呈上升趨勢(shì)，而東北大豆產(chǎn)區(qū)全樣本、吉林、內(nèi)蒙古則出現(xiàn)了下降的趨勢(shì)；2019—2020年，東北大豆產(chǎn)區(qū)全樣本、黑龍江、吉林、遼寧的大豆單產(chǎn)水平都呈現(xiàn)上升的趨勢(shì)，內(nèi)蒙古則下降；2020—2021年，所有地區(qū)的大豆單產(chǎn)均下降，甚至遼寧、內(nèi)蒙古大豆單產(chǎn)水平降到近10 a最低值。2017年以前，無論是東北大豆產(chǎn)區(qū)全樣本還是黑龍江、吉林、遼寧、內(nèi)蒙古的大豆單產(chǎn)水平整體都很低，主要原因是大豆的臨時(shí)收儲(chǔ)政策和目標(biāo)價(jià)格補(bǔ)貼并未實(shí)現(xiàn)擴(kuò)面增產(chǎn)的目標(biāo)，甚至大豆播種面積跌破歷史最低水平［14］。2017年實(shí)施大豆補(bǔ)貼政策以后，大豆單產(chǎn)水平有所提高，但是單產(chǎn)水平波動(dòng)很大，可能是因?yàn)閷?shí)施大豆生產(chǎn)者補(bǔ)貼政策后，大豆種植面積有所增加，但大豆面積增長(zhǎng)的連續(xù)性一般［14］，大豆面積的增長(zhǎng)能夠帶動(dòng)單產(chǎn)水平的提高［4］，即大豆單產(chǎn)水平受大豆播種面積的影響較大。其次不同產(chǎn)區(qū)不同品種大豆的品質(zhì)也有一定差異［26-27］，再加上新冠疫情的影響，產(chǎn)生了大豆種植管理不到位等問題。

圖2 2013—2021年東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平

3.2 東北大豆產(chǎn)區(qū)大豆單產(chǎn)時(shí)序演進(jìn)分析

采用Matlab 2021a軟件繪制了東北大豆產(chǎn)區(qū)全樣本及各省的大豆單產(chǎn)水平 Kernel密度估計(jì)圖（圖3）。

圖3 2013—2021年東北大豆產(chǎn)區(qū)全樣本及各省的大豆單產(chǎn)水平Kernel密度估計(jì)圖

由圖3a可以看出，東北大豆產(chǎn)區(qū)全樣本的大豆單產(chǎn)水平呈現(xiàn)出明顯的多峰分布，且峰值處于上升狀態(tài)，但是波峰寬度縮小，表明其存在顯著的多極化趨勢(shì)和空間非均衡特征；其次，聚集水平下降，大豆單產(chǎn)水平區(qū)域呈現(xiàn)分散狀態(tài)；波峰的寬度在逐漸擴(kuò)大，說明其聚集水平都很低，呈現(xiàn)出內(nèi)部分散的演進(jìn)態(tài)勢(shì)；此外，整個(gè)密度圖靠左分布，單產(chǎn)水平低的地區(qū)數(shù)量多，表明大豆單產(chǎn)還處于較低水平，未來提升潛力較大。從各個(gè)省（區(qū)）看，黑龍江省（圖3b）Kernel密度估計(jì)圖呈現(xiàn)出單峰分布狀態(tài)，表明黑龍江省大豆單產(chǎn)水平較低并呈現(xiàn)出單極化趨勢(shì)，且波峰高度明顯下降，寬度擴(kuò)大，表明黑龍江大豆單產(chǎn)水平逐漸分散，聚集水平降低，內(nèi)部差異隨時(shí)間的變化逐漸擴(kuò)大，整個(gè)密度圖靠左分布，表明黑龍江省大豆單產(chǎn)水平低的地區(qū)數(shù)量多；吉林省（圖3c）Kernel密度估計(jì)圖呈現(xiàn)單峰左移的狀態(tài)，表明吉林省大豆單產(chǎn)水平處于單極化，且單產(chǎn)水平低的地區(qū)數(shù)量多于單產(chǎn)水平高的地區(qū)，同時(shí)波峰高度先上升后下降，寬度擴(kuò)大，表明吉林省大豆單產(chǎn)水平穩(wěn)定性差，有一定聚集效應(yīng)，但是聚集效應(yīng)持續(xù)性弱，大豆單產(chǎn)水平內(nèi)部差異隨時(shí)間的推移逐漸變大；遼寧省Kernel密度估計(jì)圖（圖3d）呈現(xiàn)單峰左移狀態(tài)，這表明遼寧省單產(chǎn)水平低的地區(qū)數(shù)量多，波峰高度下降且寬度擴(kuò)大，說明隨著時(shí)間推移，遼寧省大豆單產(chǎn)水平下降，內(nèi)部差異變大；內(nèi)蒙古自治區(qū)（圖3e）Kernel密度估計(jì)圖呈現(xiàn)多峰分布，表明其大豆單產(chǎn)水平存在顯著的多極分化和空間非均衡特征，波峰高度先上升后下降，說明內(nèi)蒙古自治區(qū)大豆單產(chǎn)水平隨著時(shí)間推移，其穩(wěn)定性變差，聚集水平不斷下降，內(nèi)部差異逐漸擴(kuò)大。

3.3 大豆單產(chǎn)水平空間分布特征分析

為進(jìn)一步反映東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平的空間分布和演進(jìn)特征，以2013年為基準(zhǔn)年使用“四分法”將大豆單產(chǎn)水平劃分成低水平、中低水平、中高水平、高水平等4種不同類型。使用ArcGIS 10.7軟件分別對(duì)2013和2021年的大豆單產(chǎn)水平進(jìn)行空間可視化展示。由圖4可知，東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平的空間差異特征明顯，高水平地區(qū)主要集中在吉林東南部、黑龍江南部和西南部以及遼寧東北部；低水平地區(qū)主要集中在黑龍江東北部、內(nèi)蒙古東北部，在未來一段時(shí)間，低水平區(qū)域仍是大豆單產(chǎn)水平提升的重點(diǎn)區(qū)域。從省區(qū)角度看，內(nèi)蒙古自治區(qū)低水平地區(qū)數(shù)量增多，中低水平地區(qū)數(shù)量減少，負(fù)向溢出效應(yīng)增強(qiáng)。黑龍江省高水平地區(qū)數(shù)量增多，并且由東部地區(qū)轉(zhuǎn)移到西部地區(qū)，中低水平和中高水平地區(qū)數(shù)量減少，低水平地區(qū)數(shù)量增多，內(nèi)部差異化擴(kuò)大。吉林省高水平地區(qū)數(shù)量增加，高水平區(qū)域有一定聚集效應(yīng)，內(nèi)部差異隨時(shí)間推移而逐漸擴(kuò)大。從不同水平上看，東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平的空間分布格局不穩(wěn)定，2013—2021年的空間分布變化較大。高水平地區(qū)數(shù)量增加，由點(diǎn)狀分布向帶狀分布演進(jìn)，隨著時(shí)間推移，聚集效應(yīng)不斷增強(qiáng)。中高水平地區(qū)依然保持點(diǎn)狀分布，但地區(qū)數(shù)量有所下降。中低水平地區(qū)由點(diǎn)狀分布向帶狀分布演進(jìn)。部分中低水平地區(qū)轉(zhuǎn)移到低水平地區(qū)，低水平地區(qū)的數(shù)量增多。

圖4 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平空間分布圖

3.4 大豆單產(chǎn)水平的動(dòng)態(tài)轉(zhuǎn)移特征分析

為探究東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平的流動(dòng)性特征，利用Matlab 2021a軟件的傳統(tǒng)Markov鏈考察東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平分布的動(dòng)態(tài)轉(zhuǎn)移特征（表1）。

表1 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平傳統(tǒng)Markov轉(zhuǎn)移概率矩陣

由表1可知，低水平與高水平所形成的對(duì)角線兩側(cè)的數(shù)并非都小于對(duì)角線上的數(shù)字，這表明東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平受到多種因素的影響，其穩(wěn)定性差。其中，低水平地區(qū)向高水平地區(qū)轉(zhuǎn)移的概率在增大，有跨區(qū)域轉(zhuǎn)移的可能。中高水平地區(qū)向低水平地區(qū)轉(zhuǎn)移的概率增大，即中高水平地區(qū)有可能跨區(qū)域轉(zhuǎn)移到低水平地區(qū)。高水平和中低水平向低于自身水平的區(qū)域轉(zhuǎn)移的概率也在增大，這進(jìn)一步說明東北大豆產(chǎn)區(qū)內(nèi)部的不穩(wěn)定性。

由表2可知，空間Markov概率轉(zhuǎn)移矩陣除了具有與傳統(tǒng)Markov概率轉(zhuǎn)移矩陣的共性外，還能夠反映當(dāng)加入4種不同類型的水平區(qū)域以后，與水平類型高于或者低于自身的單位為鄰時(shí)，向上轉(zhuǎn)移或者向下轉(zhuǎn)移的概率。當(dāng)鄰近類型為中低水平和中高水平時(shí)，自身水平向下轉(zhuǎn)移的概率最大，說明中低水平地區(qū)和中高水平地區(qū)具有負(fù)向溢出作用；當(dāng)鄰近類型為高水平地區(qū)，自身水平跨區(qū)域轉(zhuǎn)移的概率最大，說明高水平地區(qū)具有正向溢出作用。

表2 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平空間Markov轉(zhuǎn)移概率矩陣

3.5 大豆單產(chǎn)水平的收斂性特征分析

運(yùn)用σ收斂、絕對(duì)β收斂、條件β收斂對(duì)東北大豆產(chǎn)區(qū)以及各省份大豆單產(chǎn)水平收斂性進(jìn)行檢驗(yàn)。

3.5.1σ收斂檢驗(yàn) 由表3可知，從整體上看，東北大豆產(chǎn)區(qū)的整體大豆單產(chǎn)水平σ均值為0.2963，最大值0.4372，最小值0.1670，東北大豆產(chǎn)區(qū)以及各省區(qū)的大豆單產(chǎn)水平σ值并沒有隨著時(shí)間的推移而逐漸縮小，表明東北大豆產(chǎn)區(qū)不存在σ收斂，4個(gè)省區(qū)內(nèi)部差異化各有其特殊性；從地區(qū)上看，黑龍江、吉林、遼寧、內(nèi)蒙古的大豆單產(chǎn)水平σ均值分別為0.2945、0.2647、0.2711、0.2413，其中通過黑龍江省的均值最高，說明黑龍江省各縣(市、區(qū))之間大豆單產(chǎn)水平內(nèi)部差異大。內(nèi)蒙古自治區(qū)的σ均值最低，表明內(nèi)蒙古自治區(qū)縣(市、區(qū)、旗)的內(nèi)部差異較小。

表3 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平σ收斂檢驗(yàn)結(jié)果

3.5.2 絕對(duì)β收斂檢驗(yàn) 為準(zhǔn)確地分析東北大豆產(chǎn)區(qū)全樣本及各省大豆單產(chǎn)水平的收斂情況，并檢驗(yàn)各地區(qū)大豆單產(chǎn)水平之間是否存在追趕趨勢(shì)，本文通過Hausman檢驗(yàn)選用固定效應(yīng)模型進(jìn)一步對(duì)東北大豆產(chǎn)區(qū)全樣本及各省份的大豆單產(chǎn)水平進(jìn)行絕對(duì)β收斂檢驗(yàn)（表4）。由表4可知，東北大豆產(chǎn)區(qū)及各省大豆單產(chǎn)水平的收斂系數(shù)β值在1%水平上為負(fù)顯著，存在絕對(duì)β收斂；從4省(區(qū))內(nèi)部來看，黑龍江省的β值為-0.801，是4個(gè)省里面收斂速度最快的省份，內(nèi)蒙古自治區(qū)大豆單產(chǎn)水平的β值為-0.744，表明內(nèi)蒙古的地區(qū)收斂速度最慢，4省收斂速度依次為黑龍江＞吉林＞遼寧＞內(nèi)蒙古。結(jié)果表明東北大豆產(chǎn)區(qū)縣（市、區(qū)、旗）域大豆單產(chǎn)水平的增長(zhǎng)速度與初始水平呈現(xiàn)負(fù)相關(guān)，存在顯著的絕對(duì)β收斂。

表4 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平絕對(duì)β收斂檢驗(yàn)結(jié)果

3.5.3 條件β收斂檢驗(yàn) 通過絕對(duì)β收斂檢驗(yàn)檢驗(yàn)了東北大豆產(chǎn)區(qū)全樣本及其各省間各自的穩(wěn)定狀態(tài)，利用條件β收斂進(jìn)一步檢驗(yàn)東北大豆產(chǎn)區(qū)在各自特定的生產(chǎn)條件和要素稟賦下能否趨向于各自的穩(wěn)定狀態(tài)。由表5可知，在引入4個(gè)控制變量以后，東北大豆產(chǎn)區(qū)全樣本以及各省區(qū)大豆單產(chǎn)水平的收斂系數(shù)β在1%的統(tǒng)計(jì)水平上顯著負(fù)向，說明東北大豆產(chǎn)區(qū)及其各省間存在條件β收斂，即各地區(qū)自身?xiàng)l件的與生產(chǎn)要素稟賦等方面存在明顯的差異，且各自的穩(wěn)定狀態(tài)不同，隨著時(shí)間的推移，各地區(qū)將收斂至穩(wěn)定的水平；加入控制變量以后，各地區(qū)的收斂速度出現(xiàn)了不同的變化，最為明顯的是內(nèi)蒙古的收斂速度，β值為-0.735，4省收斂速度依次為內(nèi)蒙古＞遼寧＞吉林＞黑龍江，這說明如果不考慮要素稟賦與生產(chǎn)條件的絕對(duì)β收斂模型則會(huì)低估了內(nèi)蒙古和遼寧地區(qū)大豆單產(chǎn)水平的收斂性。

表5 東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平條件β收斂檢驗(yàn)結(jié)果

4 小結(jié)與建議

4.1 小結(jié)

本文基于2013—2021年東北大豆產(chǎn)區(qū)141個(gè)縣的面板數(shù)據(jù)，測(cè)算了東北大豆產(chǎn)區(qū)的大豆單水平，并對(duì)其時(shí)空演進(jìn)和收斂性進(jìn)行分析。結(jié)果表明：從時(shí)序角度上看，東北大豆產(chǎn)區(qū)的大豆單產(chǎn)水平呈現(xiàn)出先增后減的時(shí)序特征，東北大豆產(chǎn)區(qū)及各省份的大豆單產(chǎn)水平穩(wěn)定性差，波動(dòng)劇烈，穩(wěn)定性差的地區(qū)未來仍存在較大的提升空間；從空間角度上看，Kernel密度估計(jì)和大豆單產(chǎn)水平空間特征表明東北大豆產(chǎn)區(qū)大豆單產(chǎn)水平呈現(xiàn)出明顯的多峰分布，存在多極化趨勢(shì)和空間非均衡特征，高水平地區(qū)數(shù)量增加，由點(diǎn)狀分布向帶狀分布演進(jìn)，隨著時(shí)間推移，聚集效應(yīng)不斷增強(qiáng)，中低水平地區(qū)和中高水平地區(qū)存在顯著的空間負(fù)相關(guān)，容易導(dǎo)致自身水平以及鄰域水平向下轉(zhuǎn)移，低水平地區(qū)數(shù)量增加，由點(diǎn)狀分布演變成為帶狀分布，進(jìn)一步說明低水平地區(qū)具有負(fù)向溢出效應(yīng)；在收斂性方面，東北大豆產(chǎn)區(qū)全樣本以及各省均不存在σ收斂，但存在顯著的絕對(duì)β收斂和條件β收斂，如果忽略控制變量，則內(nèi)蒙古和遼寧地區(qū)大豆單產(chǎn)水平的收斂性會(huì)被低估，即要素稟賦與生產(chǎn)條件對(duì)于大豆單產(chǎn)水平有顯著的影響。

4.2 建議

基于以上結(jié)論，提出如下建議：

（1）統(tǒng)籌規(guī)劃東北大豆產(chǎn)區(qū)的生產(chǎn)工作。完善東北地區(qū)大豆生產(chǎn)的頂層設(shè)計(jì)，鼓勵(lì)各地區(qū)積極突破行政邊界的外在約束，協(xié)同推進(jìn)大豆生產(chǎn)合理布局和配套基礎(chǔ)設(shè)施建設(shè)。開展中低水平地區(qū)和中高水平地區(qū)的產(chǎn)田改造工作，提高低水平地區(qū)的生產(chǎn)能力。嘗試在不同縣域差異化實(shí)施大豆生產(chǎn)補(bǔ)貼政策，力爭(zhēng)實(shí)現(xiàn)一縣一策，通過提高種植面積穩(wěn)定大豆的單產(chǎn)水平。

（2）加強(qiáng)東北大豆產(chǎn)區(qū)高標(biāo)準(zhǔn)農(nóng)田的建設(shè)工作。建立區(qū)域合作機(jī)制，促進(jìn)高水平地區(qū)與周邊地區(qū)的合作，充分利用高水平地區(qū)的正向溢出效應(yīng)。加大對(duì)高水平地區(qū)的財(cái)政補(bǔ)貼、政策傾斜支持力度，因地制宜地推進(jìn)高水平地區(qū)的良種改造。

（3）加強(qiáng)東北大豆產(chǎn)區(qū)的區(qū)域協(xié)作。各地政府應(yīng)著力破除大豆生產(chǎn)資源要素跨地區(qū)流動(dòng)障礙，促進(jìn)地區(qū)間土地、資金及機(jī)械等要素的高效流動(dòng)和合理集聚。推動(dòng)?xùn)|北產(chǎn)區(qū)整體的大豆生產(chǎn)水平提高，需要加強(qiáng)生產(chǎn)、流通、消費(fèi)環(huán)節(jié)的統(tǒng)籌協(xié)作，充分促動(dòng)鄰近地區(qū)大豆單產(chǎn)水平的協(xié)同提升。