中國農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的互動(dòng)關(guān)系研究

董菲

摘要 本文梳理了農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)互動(dòng)發(fā)展的階段特征，并基于產(chǎn)業(yè)共生視角建立分時(shí)段Logistic共生演化模型，利用中國2003—2019年相關(guān)數(shù)據(jù)對(duì)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)發(fā)展過程進(jìn)行模擬演化。研究發(fā)現(xiàn)二者發(fā)展存在明顯共生關(guān)系，目前階段農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的促進(jìn)作用更強(qiáng)，屬于非對(duì)稱互惠共生模式。基于研究結(jié)果提出通過提高農(nóng)業(yè)服務(wù)的戰(zhàn)略地位、提升農(nóng)業(yè)產(chǎn)業(yè)化水平等方式促進(jìn)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)由非對(duì)稱互惠共生向?qū)ΨQ互惠共生演進(jìn)。

關(guān)鍵詞 農(nóng)業(yè)服務(wù)；農(nóng)業(yè)生產(chǎn)；共生演化；Logistic共生演化模型

中圖分類號(hào) F326.6 ? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 1007-7731（2023）11-0179-08

Research on the interaction between agricultural service and agricultural production in China

—— Based on Logistic symbiosis evolution model

DONG Fei

（College of Mathematics and Information Science， Xianyang Normal University， Xianyang 712000， China）

Abstract This paper sorts out the stage characteristics of the interactive development of agricultural service and agricultural production， and establishes the time-based Logistic symbiotic evolution model based on the perspective of industrial symbiosis， and simulates the evolution of agricultural service and agricultural production development process by using the relevant Chinese data from 2003 to 2019. The study found that there is an obvious symbiotic relationship between the two companies in their development. At present， agricultural services have a stronger role in promoting agricultural production， which belongs to the asymmetric mutualism mode. Based on the research results， it is proposed that the evolution of agricultural service and agricultural production should be promoted from asymmetric reciprocity to symmetrical mutualism by improving the strategic position of agricultural service and improving the level of agricultural industrialization.

Keywords agricultural services; agricultural production; symbiosis evolution; Logistic symbiosis evolution model

隨著農(nóng)業(yè)分工體系的進(jìn)一步深化以及生產(chǎn)性服務(wù)業(yè)與其他產(chǎn)業(yè)的不斷融合，傳統(tǒng)農(nóng)業(yè)活動(dòng)中的一些非生產(chǎn)環(huán)節(jié)逐漸從農(nóng)業(yè)生產(chǎn)流程剝離出來，發(fā)展成為以農(nóng)業(yè)生產(chǎn)性服務(wù)為主的新興農(nóng)業(yè)服務(wù)行業(yè)，這些行業(yè)同農(nóng)業(yè)生產(chǎn)之間通過要素交換與商品交換相聯(lián)系，二者在長期發(fā)展中形成了穩(wěn)定的相互依賴關(guān)系。

目前，關(guān)于二者關(guān)系的研究主要集中在農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的單向促進(jìn)作用方面，研究認(rèn)為農(nóng)業(yè)服務(wù)將人才、資本、技術(shù)、信息等現(xiàn)代生產(chǎn)要素匯集到農(nóng)業(yè)生產(chǎn)活動(dòng)中，可以促進(jìn)農(nóng)民專業(yè)化水平的提升，提高農(nóng)業(yè)生產(chǎn)的技術(shù)水平[1]及生產(chǎn)效率[2-3]，發(fā)揮農(nóng)地經(jīng)營的規(guī)模經(jīng)濟(jì)效應(yīng)[4]，在產(chǎn)前、產(chǎn)中、產(chǎn)后各個(gè)環(huán)節(jié)促進(jìn)農(nóng)業(yè)生產(chǎn)現(xiàn)代化發(fā)展。此外，還能破解農(nóng)戶兼業(yè)化難題[5]，進(jìn)而帶動(dòng)農(nóng)民收入增長，增加農(nóng)村居民消費(fèi)[6]，縮小城鄉(xiāng)居民收入差距[7]，不僅能夠穩(wěn)固農(nóng)業(yè)基礎(chǔ)，還對(duì)周邊地區(qū)存在著顯著的空間溢出作用[8]。但以往的研究在一定程度上忽略了農(nóng)業(yè)服務(wù)發(fā)展的生產(chǎn)需求基礎(chǔ)、忽視了其與農(nóng)業(yè)生產(chǎn)之間的互動(dòng)規(guī)律及共生特征。隨著技術(shù)、政策以及市場(chǎng)等發(fā)展環(huán)境因素的優(yōu)化，產(chǎn)業(yè)融合與共生現(xiàn)象逐漸普遍，農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間共生發(fā)展已經(jīng)成為重要趨勢(shì)。本文將梳理農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)互動(dòng)發(fā)展的階段特征，進(jìn)而建立農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的共生演化模型，通過共生演化模型揭示農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)協(xié)同發(fā)展的過程以及共生模式。

“共生”在生態(tài)學(xué)中表示不同物種的有機(jī)體在進(jìn)化過程中相互之間形成緊密相伴或相互受益的生態(tài)關(guān)系[9]，也可以描述為2個(gè)不同的種群經(jīng)過長時(shí)間適應(yīng)和改造環(huán)境，形成彼此合作、相互關(guān)聯(lián)的協(xié)同進(jìn)化關(guān)系[10-11]。近年來，隨著共生理論的發(fā)展，共生研究也被廣泛應(yīng)用于經(jīng)濟(jì)、科技、醫(yī)藥等領(lǐng)域，主要基于Logistic方程描述種群共生關(guān)系。例如：徐學(xué)軍、唐強(qiáng)榮、龐博慧、金晟等[12-15]學(xué)者利用Logistic共生演化模型分析服務(wù)業(yè)與制造業(yè)之間的互動(dòng)關(guān)系及影響因素；在此基礎(chǔ)上，于斌斌和胡漢輝[16]基于服務(wù)業(yè)與制造業(yè)互動(dòng)共生視角，探討了產(chǎn)業(yè)集群與城市化的共生演化機(jī)制；吳勇民等[17]發(fā)現(xiàn)金融產(chǎn)業(yè)與高新技術(shù)產(chǎn)業(yè)之間存在非對(duì)稱性互惠共生模式；以及基于Logistic方程的共生模型在創(chuàng)新生態(tài)系統(tǒng)[18]、區(qū)域物流生態(tài)系統(tǒng)[19]、隱性知識(shí)流轉(zhuǎn)網(wǎng)[20]等領(lǐng)域的創(chuàng)新應(yīng)用。除此以外，還有學(xué)者基于投入產(chǎn)出分析法[21]、產(chǎn)業(yè)耦合協(xié)調(diào)度模型[22-23]、協(xié)同度模型[24]等方法刻畫了產(chǎn)業(yè)之間的共生關(guān)系。

縱觀已有研究成果，尚未有研究將共生演化模型應(yīng)用于農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)互動(dòng)關(guān)系方面，本研究有利于從微觀上把握農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)協(xié)同發(fā)展的動(dòng)態(tài)特征，為提高農(nóng)業(yè)生產(chǎn)現(xiàn)代化、農(nóng)業(yè)服務(wù)化提供新的思路。

1 共生演化模型的構(gòu)建

1.1 Logistic模型

農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)共生演化的過程也就是2個(gè)共生單元“種群”數(shù)量不斷增長的過程，從發(fā)展周期來看，共生系統(tǒng)的演化會(huì)經(jīng)歷萌芽期、成長期、發(fā)展期、成熟期、一體化時(shí)期5個(gè)主要階段，“種群”數(shù)量的增加會(huì)受到資源稟賦、社會(huì)環(huán)境、政策和制度環(huán)境等外部環(huán)境的約束，用Logistic模型能夠比較直觀且準(zhǔn)確地描述農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的共生演化特征。

用[Y1（t）]、[Y2（t）]分別表示農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)在[t]時(shí)期的種群密度，自然增長率分別用[n1]、[n2]表示，當(dāng)給定資源要素稟賦和外部環(huán)境（包括技術(shù)水平、市場(chǎng)規(guī)模等）時(shí)，不同種群根據(jù)其各自的發(fā)展特性，在環(huán)境制約下存在一個(gè)最大種群密度，也稱為最大環(huán)境容量，分別用[M1]、[M2]表示。當(dāng)2個(gè)產(chǎn)業(yè)獨(dú)立發(fā)展、互不影響時(shí)，它們的動(dòng)態(tài)演化方程可以表示如下：

[dY1（t）dt=n11-Y1（t）M1Y1（t），Y1（t0）=Y10dY2（t）dt=n21-Y2（t）M2Y2（t），Y2（t0）=Y20] （1）

式中，[Y10]表示農(nóng)業(yè)服務(wù)初始種群密度，[Y20]表示農(nóng)業(yè)生產(chǎn)初始種群密度。種群規(guī)模的擴(kuò)張過程伴隨著對(duì)外部資源的消耗，因此，當(dāng)種群規(guī)模增長到一定程度時(shí)，受資源稟賦約束，其增長速度會(huì)逐漸放慢，分別用[1-Y1（t）M1]和[1-Y2（t）M2]描述農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)發(fā)展過程中的阻滯因子。

假設(shè)同一系統(tǒng)內(nèi)只存在2個(gè)種群相互作用，每一個(gè)種群的發(fā)展既受到特定環(huán)境下自身發(fā)展規(guī)律的影響，還受到另一個(gè)種群密度的影響。基于Logistic種群增長方程，考慮到農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的相互影響，建立2種群相互影響條件下的共生演化模型。龐博慧[25]、吳勇民[17]等學(xué)者在研究中直接引入了另一種群密度作為影響因子，考慮到阻滯因子的形式以及種群密度與環(huán)境容納量的比值能夠更加準(zhǔn)確地反映另一種群對(duì)資源的競(jìng)爭情況，本文在這里引入另一種群密度占最大環(huán)境容量（飽和密度）比重作為交互影響因素，構(gòu)建模型如下：

[dY1（t）dt=n11-Y1（t）M1+λ12Y2（t）M2Y1（t）dY2（t）dt=n21-Y2（t）M2+λ21Y1（t）M1Y2（t）] （2）

式中，[λ12]為農(nóng)業(yè)生產(chǎn)對(duì)農(nóng)業(yè)服務(wù)的共生作用系數(shù)，[λ21]為農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的共生作用系數(shù)，它們的取值范圍界定了農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的共生模式，見表1。

1.2 模型變換

為了更加有效地研究農(nóng)業(yè)生產(chǎn)與農(nóng)業(yè)服務(wù)之間的共生演化關(guān)系，需對(duì)式（2）進(jìn)行變換，令

[ω1=n11+λ12Y2（t）M2， ω2=n21+λ21Y1（t）M1K1=M11+λ12Y2（t）M2，K2=M21+λ21Y1（t）M1] ?（3）

則公式（2）可變換為：

[dY1（t）dt=ω11-Y1（t）K1Y1（t）dY2（t）dt=ω21-Y2（t）K2Y2（t）] （4）

公式（4）在形式上與公式（1）一致，但二者所表示的內(nèi)容與符號(hào)含義不同。公式（1）中的農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)獨(dú)立發(fā)展，互不影響，可以認(rèn)為種群自然增長率與最大環(huán)境容量不會(huì)發(fā)生變化，為常數(shù)；公式（4）中的農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)存在共生關(guān)系，相互影響，種群自然增長率與最大環(huán)境容量隨著另一種群密度的改變而變化。

公式（4）中[ω1]和[ω2]分別為農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)共生演化發(fā)展時(shí)的各自對(duì)應(yīng)的共生自然增長率，[K1]和[K2]為相互影響機(jī)制下各自的共生最大環(huán)境容量。隨著種群規(guī)模的發(fā)展，農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)各自的共生自然增長率及共生最大環(huán)境容量也會(huì)不斷發(fā)生變化，進(jìn)而表現(xiàn)出種群間共生演化的動(dòng)態(tài)趨勢(shì)。

1.3 分時(shí)段疊加Logistic模型

式（4）對(duì)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的共生演化機(jī)制進(jìn)行了描述，為了進(jìn)一步刻畫共生演化的具體過程，并確定模型中各參數(shù)隨時(shí)間變化的數(shù)值，需要對(duì)原模型進(jìn)行分時(shí)段疊加處理。借鑒唐強(qiáng)榮等[26]學(xué)者做法，用年區(qū)分農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的種群演化時(shí)段，則種群密度演化Logistic曲線可以近似看成由各個(gè)較小時(shí)間區(qū)段內(nèi)[t∈ti，ti+1， i=0，1， 2...T]的Logistic曲線連接而成。因此，式（4）在較小時(shí)間區(qū)段[ti，ti+1]內(nèi)的解為：

[Y1（ti+1）=Kt+111+Kt+11-Y1（ti）Y1（ti）e-ωi+11?ΔtY2（ti+1）=Kt+121+Kt+12-Y2（ti）Y2（ti）e-ωi+12?Δt] （5）

取任意相鄰時(shí)間區(qū)間[ti，ti+1]，區(qū)間長度[Δt=ti+1-ti=1]，設(shè)農(nóng)業(yè)服務(wù)在種群密度演化曲線上時(shí)間區(qū)間[ti，ti+1]上的種群密度增量為[ΔY1（ti+1）=Y1（ti+1）-Y1（ti）]，這一時(shí)間區(qū)間內(nèi)平均種群密度為[Y1（ti+1）=Y1（ti+1）+Y1（ti）2]，連接兩端點(diǎn)的直線斜率為[ΔY1（ti+1）Δt]；同理，對(duì)于農(nóng)業(yè)生產(chǎn)在時(shí)間區(qū)間[ti，ti+1]上的種群密度增量為[ΔY2（ti+1）=Y2（ti+1）-Y2（ti）]，這一時(shí)間區(qū)間內(nèi)平均種群密度為[Y2（ti+1）=Y2（ti+1）+Y2（ti）2]，連接兩端點(diǎn)的直線斜率為[ΔY2（ti+1）Δt]。

假設(shè)在較小的時(shí)間間隔內(nèi)種群密度曲線斜率變化不大，曲線上各點(diǎn)的斜率均可以用對(duì)應(yīng)時(shí)間區(qū)間內(nèi)的直線斜率近似替代，因此式（4）可以轉(zhuǎn)換為分時(shí)段疊加Logistic共生演化模型，公式如下：

[ΔY1（ti+1）Δt=ωi+111-Y1（ti+1）Ki+11Y1（ti+1）ΔY2（ti+1）Δt=ωi+121-Y2（ti+1）Ki+12Y2（ti+1）] （6）

1.4 模型參數(shù)識(shí)別

對(duì)式（6）進(jìn)行變換，得到共生最大環(huán)境容量的迭代公式，如下：

[Ki+11=Y1（ti+1）1-ΔY1（ti+1）Δt?ωi+11?Y1（ti+1）Ki+12=Y2（ti+1）1-ΔY2（ti+1）Δt?ωi+12?Y2（ti+1）] （7）

由于共生最大環(huán)境容量必然大于0，即[Ki+11>0，]

[Ki+12>0]，則有

[ωi+11>ΔY1（ti+1）Δt?Y1（ti+1）ωi+12>ΔY2（ti+1）Δt?Y2（ti+1）] （8）

式中，[Y1（ti+1）]、[Y2（ti+1）]、[ΔY1（ti+1）]和[ΔY2（ti+1）]可以由調(diào)查及搜集到的數(shù)據(jù)計(jì)算得出，且[Δt=ti+1-ti=1]，則可以在限定范圍內(nèi)給定一組[ωi+11]、[ωi+12][（i=0，1，2…T）]的估計(jì)值[ωi+11]、[ωi+12][（i=0，1，2…T）]，代入式（7），進(jìn)而可以計(jì)算出各個(gè)時(shí)間區(qū)間內(nèi)的共生最大環(huán)境容量估計(jì)值[Ki+11]、[Ki+12][（i=0，1，2…T）]。將[ωi+11]、[ωi+12]、[Ki+11]、[Ki+12][（i=0，1，2…T）] 4組估計(jì)值代入式（5），可以得到農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)共生時(shí)的種群密度估計(jì)值：

[Y1（ti+1）=Kt+111+Kt+11-Y1（ti）Y1（ti）e-ωi+11?ΔtY2（ti+1）=Kt+121+Kt+12-Y2（ti）Y2（ti）e-ωi+12?Δt] （9）

進(jìn)而獲得農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)在考察期間種群密度的第m次估計(jì)值與實(shí)際值的殘差平方和：

[Δe21（m）=i=0TY1（ti+1）-Y1（ti+1）2Δe22（m）=i=0TY2（ti+1）-Y2（ti+1）2] （10）

通過二次規(guī)劃求解受約束非線性方程組，得到共生條件下分時(shí)段自然增長率的值[ωi+11]、[ωi+12]，將其代入式（7），可以得到考察期內(nèi)共生最大環(huán)境容量的值[Ki+11]、[Ki+12]。通過以上方法，即可獲得分時(shí)段Logistic共生動(dòng)態(tài)演化模型的各參數(shù)。

2 數(shù)據(jù)來源與實(shí)證分析

2.1 數(shù)據(jù)說明

為了檢驗(yàn)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的互動(dòng)關(guān)系是否能夠用分時(shí)段Logistic共生演化模型準(zhǔn)確描述，本文選取中國2003—2019年農(nóng)林牧漁服務(wù)業(yè)增加值、第一產(chǎn)業(yè)增加值分別反映農(nóng)業(yè)服務(wù)增加值與農(nóng)業(yè)生產(chǎn)增加值，并進(jìn)行實(shí)證分析。根據(jù)《國民經(jīng)濟(jì)行業(yè)分類》（GB/T 4754—2011），自2012年起第一產(chǎn)業(yè)增加值不含農(nóng)林牧漁服務(wù)業(yè)，農(nóng)林牧漁增加值包含相關(guān)服務(wù)業(yè)，因此可用農(nóng)林牧漁增加值數(shù)據(jù)減去第一產(chǎn)業(yè)增加值數(shù)據(jù)估算2012—2019年農(nóng)林牧漁服務(wù)業(yè)增加值。對(duì)于2003—2011年農(nóng)業(yè)服務(wù)增加值，結(jié)合農(nóng)林牧漁增加值與第一產(chǎn)業(yè)增加值差額，以及2012—2019年農(nóng)林牧漁服務(wù)業(yè)增加值占農(nóng)林牧漁增加值比重變動(dòng)情況，采用加權(quán)法估算。農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)增加值的變化過程即為二者發(fā)展演化過程，并且增加值變動(dòng)受到資源稟賦、社會(huì)環(huán)境、技術(shù)效率、市場(chǎng)規(guī)模等因素的制約，不能無限增長，與種群密度變化特征相似。因此，本文選取農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的增加值作為各自種群密度指標(biāo)，數(shù)據(jù)來源于2003—2019年的《中國統(tǒng)計(jì)年鑒》及《中國農(nóng)村統(tǒng)計(jì)年鑒》。

2.2 模型檢驗(yàn)

首先假設(shè)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)獨(dú)立發(fā)展，即[λ12=λ21=0]。利用Nelder-Mead Simplex算法對(duì)公式（1）中的參數(shù)[n]、[M]進(jìn)行估計(jì)，可以得到在獨(dú)立發(fā)展假設(shè)下，農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的自然增長率和最大環(huán)境容量，圖1與圖2分別為獨(dú)立發(fā)展假設(shè)下農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的動(dòng)態(tài)演化Logistic曲線。

由圖1可以看出在獨(dú)立發(fā)展假設(shè)下，農(nóng)業(yè)服務(wù)動(dòng)態(tài)演化曲線對(duì)實(shí)際觀測(cè)值的擬合效果較好，通過模型估計(jì)出的獨(dú)立發(fā)展假設(shè)下農(nóng)業(yè)服務(wù)自然增長率為0.170 1，最大環(huán)境容量為5 577.4億元。2019年農(nóng)業(yè)服務(wù)實(shí)際增加值為3 103.3億元，種群密度與最大環(huán)境容量的比值約為55.64%。根據(jù)Logistic曲線S型特征，這一階段的農(nóng)業(yè)服務(wù)發(fā)展進(jìn)入轉(zhuǎn)折后期，增長規(guī)模將逐步減小，從2015—2019年，農(nóng)業(yè)服務(wù)增加值變動(dòng)分別為232.1、233.8、248.7、253.0、289.8億元，增長規(guī)模仍然不斷增加，因此判斷農(nóng)業(yè)服務(wù)發(fā)展應(yīng)處于加速期或轉(zhuǎn)折前期。根據(jù)圖2可以看出農(nóng)業(yè)生產(chǎn)發(fā)展動(dòng)態(tài)模擬曲線呈現(xiàn)S型，模型擬合出的農(nóng)業(yè)生產(chǎn)自然增長率為0.181 3，最大環(huán)境容量為83 552.5億元。2019年農(nóng)業(yè)生產(chǎn)實(shí)際增加值為70 473.6億元，種群密度與最大環(huán)境容量的比值為84.35%，推測(cè)農(nóng)業(yè)生產(chǎn)發(fā)展處于飽和期，增長規(guī)模與增長幅度均逐漸趨近于0。從2015—2019年，農(nóng)業(yè)生產(chǎn)增加值變動(dòng)分別為2 148.3、2 364.6、1 960.3、2 645.7、5 728.4億元，增長幅度分別為3.78%、4.00%、3.20%、4.17%、8.47%。可以看出農(nóng)業(yè)生產(chǎn)經(jīng)歷階段性飽和后，自2017年起重新進(jìn)入新一輪加速期。

根據(jù)上述分析，可以看出獨(dú)立發(fā)展假設(shè)下的農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)發(fā)展動(dòng)態(tài)模擬曲線都不能準(zhǔn)確描述二者實(shí)際動(dòng)態(tài)變化過程。因此，需要考慮二者之間的互動(dòng)共生因素，建立分時(shí)段共生演化模型。

2.3 實(shí)證結(jié)果分析

利用序列二次規(guī)劃法對(duì)公式（6）進(jìn)行參數(shù)識(shí)別，得到2004—2019年農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)在共生發(fā)展假設(shè)下的分時(shí)段共生最大環(huán)境容量與共生自然增長率，具體結(jié)果見表2和表3。

在農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)相互作用下，動(dòng)態(tài)演化曲線如圖3和圖4所示，結(jié)合表2、表3，可以看出由于互動(dòng)共生因素影響，農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)共生最大環(huán)境容量不斷增加，2019年分別達(dá)到8 662.7與194 309.8 億元，遠(yuǎn)高于各自獨(dú)立發(fā)展假設(shè)下的最大環(huán)境容量5 577.4和83 552.5 億元。從圖3可以看出中國農(nóng)業(yè)服務(wù)的增加值迅速增長，2019年種群密度占最大環(huán)境容量的35.82%，說明農(nóng)業(yè)服務(wù)在共生發(fā)展?fàn)顟B(tài)下正處于成長期；根據(jù)圖4可以看出農(nóng)業(yè)生產(chǎn)發(fā)展呈現(xiàn)出階段性成長特征，自2009—2017年經(jīng)歷了一輪成長、發(fā)展、成熟周期后，2018年起又開始保持高速增長，2019年的種群密度占最大環(huán)境容量的36.27%，說明在共生狀態(tài)下農(nóng)業(yè)生產(chǎn)重新進(jìn)入成長期。

綜合獨(dú)立發(fā)展假設(shè)與共生發(fā)展假設(shè)分析結(jié)果，獨(dú)立發(fā)展演化過程不能準(zhǔn)確描述農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的實(shí)際發(fā)展情形，加入互動(dòng)共生因素后，分時(shí)段共生演化模型能夠更加準(zhǔn)確地反映中國農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的增加值歷年動(dòng)態(tài)變化過程。由于農(nóng)業(yè)服務(wù)和農(nóng)業(yè)生產(chǎn)之間存在互動(dòng)影響，因此各自共生最大環(huán)境容量發(fā)生了變化，根據(jù)公式（3），利用遺傳算法得到農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的共生作用系數(shù)[λ12]與[λ21]的估計(jì)值分別為0.500和1.138。共生作用系數(shù)均大于0，且不相等，說明中國農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)之間的共生模式屬于非對(duì)稱互惠共生。此外，農(nóng)業(yè)生產(chǎn)對(duì)農(nóng)業(yè)服務(wù)的共生作用系數(shù)[λ12]小于農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的共生作用系數(shù)[λ21]，說明在二者共生發(fā)展過程中，農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的促進(jìn)作用明顯大于農(nóng)業(yè)生產(chǎn)對(duì)農(nóng)業(yè)服務(wù)的促進(jìn)作用。

3 結(jié)論與建議

本文通過梳理農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)互動(dòng)發(fā)展的階段特征，建立了分時(shí)段Logistic共生演化模型，以此為基礎(chǔ)，搜集相關(guān)數(shù)據(jù)對(duì)2003—2019年農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的發(fā)展過程進(jìn)行模擬演化，發(fā)現(xiàn)二者發(fā)展存在明顯共生關(guān)系，農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的促進(jìn)作用更強(qiáng)，屬于非對(duì)稱互惠共生。但在共生發(fā)展模式中，對(duì)稱互惠共生才是最穩(wěn)定、最高效的共生系統(tǒng)，在對(duì)稱互惠共生階段，共生系統(tǒng)中的各產(chǎn)業(yè)生產(chǎn)邊界模糊，形成一體化融合趨勢(shì)，這也是產(chǎn)業(yè)互動(dòng)發(fā)展的最終方向。因此，需要在本文研究基礎(chǔ)上進(jìn)一步探討中國農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)由非對(duì)稱互惠共生模式向?qū)ΨQ互惠共生模式的演進(jìn)機(jī)制與途徑。

當(dāng)前，中國農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的促進(jìn)作用大于農(nóng)業(yè)生產(chǎn)對(duì)農(nóng)業(yè)服務(wù)的作用，主要原因之一是過去中國農(nóng)業(yè)屬于粗放型生產(chǎn)模式，技術(shù)效應(yīng)與規(guī)模效應(yīng)雙低。隨著近年來農(nóng)業(yè)服務(wù)的不斷發(fā)展，人才、技術(shù)、信息、資本等逐漸匯集到農(nóng)業(yè)生產(chǎn)過程中，使得農(nóng)業(yè)生產(chǎn)的發(fā)展重新進(jìn)入成長期。但當(dāng)前中國農(nóng)業(yè)生產(chǎn)發(fā)展仍處于轉(zhuǎn)型階段，對(duì)農(nóng)業(yè)服務(wù)的反哺主要來源于轉(zhuǎn)型過程中部分生產(chǎn)環(huán)節(jié)對(duì)農(nóng)業(yè)服務(wù)的需求以及勞動(dòng)力轉(zhuǎn)移帶來的微弱促進(jìn)。并且根據(jù)實(shí)證分析結(jié)果，在共生機(jī)制下農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)種群密度各自占最大環(huán)境容量的35.82%和36.27%，二者都尚處于成長期，這一階段的共生模式主要為偏利共生或非對(duì)稱互惠共生，而計(jì)算得到的共生作用系數(shù)為0.500和1.138，屬于偏利共生與非對(duì)稱互惠共生的過渡后期。另一方面，自2003—2019年，中國農(nóng)業(yè)服務(wù)占農(nóng)業(yè)增加值比重由2.4%上升至4.2%，占比偏低且增長緩慢，這說明農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)融合程度不足，農(nóng)業(yè)內(nèi)部產(chǎn)業(yè)結(jié)構(gòu)層次偏低。未來隨著農(nóng)業(yè)生產(chǎn)發(fā)展水平逐步提高，農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)的融合不斷加深，二者聯(lián)系更加緊密、農(nóng)業(yè)服務(wù)對(duì)農(nóng)業(yè)生產(chǎn)的依賴性增強(qiáng)，農(nóng)業(yè)生產(chǎn)的反哺效應(yīng)也將隨之凸顯。

因此，實(shí)現(xiàn)中國農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)由非對(duì)稱互惠共生向?qū)ΨQ互惠共生模式轉(zhuǎn)變，必須重點(diǎn)提高農(nóng)業(yè)生產(chǎn)發(fā)展水平，促進(jìn)農(nóng)業(yè)服務(wù)與農(nóng)業(yè)生產(chǎn)融合，具體建議如下。

（1）提升農(nóng)業(yè)服務(wù)的戰(zhàn)略地位。通過政策扶持、補(bǔ)貼獎(jiǎng)勵(lì)、政府組織等方式，利用專項(xiàng)激勵(lì)措施支持農(nóng)業(yè)服務(wù)公司及組織、專業(yè)技術(shù)協(xié)會(huì)等深入農(nóng)業(yè)產(chǎn)業(yè)鏈，為農(nóng)業(yè)生產(chǎn)提供便利化、高效化、全面化服務(wù)。同時(shí)，對(duì)于外部性較強(qiáng)、盈利能力較低，但在農(nóng)業(yè)生產(chǎn)中發(fā)揮重要作用的服務(wù)環(huán)節(jié)，需要由政府牽頭，成立公益性組織，積極解決農(nóng)業(yè)生產(chǎn)難題，提高農(nóng)業(yè)服務(wù)水平。

（2）提高農(nóng)業(yè)產(chǎn)業(yè)化水平。依托農(nóng)業(yè)產(chǎn)業(yè)鏈、價(jià)值鏈，大力培育龍頭企業(yè)，推動(dòng)產(chǎn)業(yè)集群發(fā)展。加強(qiáng)農(nóng)業(yè)基礎(chǔ)性設(shè)施建設(shè)，確保農(nóng)業(yè)生產(chǎn)信息流通、物流暢通，為農(nóng)業(yè)生產(chǎn)發(fā)展夯實(shí)基礎(chǔ)。完善土地政策與要素配置機(jī)制，搭建土地流轉(zhuǎn)交易平臺(tái)，促進(jìn)土地經(jīng)營權(quán)高效流轉(zhuǎn)，提高農(nóng)業(yè)生產(chǎn)規(guī)模效應(yīng)。鼓勵(lì)各類資金進(jìn)入農(nóng)業(yè)生產(chǎn)領(lǐng)域，培育多元投資主體，支持農(nóng)業(yè)金融服務(wù)創(chuàng)新，提高農(nóng)業(yè)資本配置效率。培育新型農(nóng)業(yè)經(jīng)營主體，完善人才培養(yǎng)機(jī)制和就業(yè)創(chuàng)業(yè)機(jī)制，鼓勵(lì)和引導(dǎo)優(yōu)秀人才返鄉(xiāng)創(chuàng)業(yè)。完善農(nóng)產(chǎn)品市場(chǎng)體系，利用現(xiàn)代信息技術(shù)與物流技術(shù)科學(xué)布局農(nóng)產(chǎn)品流通渠道及定價(jià)機(jī)制。

參考文獻(xiàn)

[1] 郝愛民.農(nóng)業(yè)生產(chǎn)性服務(wù)對(duì)農(nóng)業(yè)技術(shù)進(jìn)步貢獻(xiàn)的影響[J].華南農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2015，14（1）：8-15.

[2] 張忠軍，易中懿.農(nóng)業(yè)生產(chǎn)性服務(wù)外包對(duì)水稻生產(chǎn)率的影響研究——基于358個(gè)農(nóng)戶的實(shí)證分析[J].農(nóng)業(yè)經(jīng)濟(jì)問題，2015，36（10）：69-76.

[3] 陳明，李文秀.生產(chǎn)服務(wù)業(yè)開放對(duì)中國農(nóng)業(yè)生產(chǎn)率的影響[J].華南農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2018，17（5）：12-23.

[4] 劉強(qiáng)，楊萬江.農(nóng)戶行為視角下農(nóng)業(yè)生產(chǎn)性服務(wù)對(duì)土地規(guī)模經(jīng)營的影響[J].中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2016，21（9）：188-197.

[5] 劉振中.我國服務(wù)業(yè)與工業(yè)、農(nóng)業(yè)互動(dòng)發(fā)展的階段性特征與展望[J].經(jīng)濟(jì)縱橫，2018（8）：59-69.

[6] 薛賀香.城鎮(zhèn)化、農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)與農(nóng)村居民消費(fèi)互動(dòng)的實(shí)證研究[J].廣東商學(xué)院學(xué)報(bào)，2013，28（6）：81-88.

[7] 魯釗陽.農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)發(fā)展對(duì)城鄉(xiāng)收入差距的影響[J].南京社會(huì)科學(xué)，2013（2）：23-29.

[8] 張薦華，高軍.發(fā)展農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)會(huì)縮小城鄉(xiāng)居民收入差距嗎？——基于空間溢出和門檻特征的實(shí)證檢驗(yàn)[J].西部論壇，2019，29（1）：45-54.

[9] DE BARY A. Die erscheinung der symbiose[M]. Strassburg：De Gruyter，1879.

[10] AVIL?S L，ABBOT P，CUTTER A D. Population ecology，nonlinear dynamics，and social evolution. I. Associations among Nonrela－tives[J]. The American Naturalist，2002，159（2）：115-127.

[11] NEUHAUSER C，F(xiàn)ARGIONE J E. A mutualism-parasitism continuum model and its application to plant-mycorrhizae interactions[J]. Ecological Modelling，2004，177（3-4）：337-352.

[12] 徐學(xué)軍，唐強(qiáng)榮，樊奇.中國生產(chǎn)性服務(wù)業(yè)與制造業(yè)種群的共生——基于Logistic生長方程的實(shí)證研究[J].管理評(píng)論，2011，23（9）：152-159.

[13] 唐強(qiáng)榮，徐學(xué)軍.基于共生理論的生產(chǎn)性服務(wù)企業(yè)與制造企業(yè)合作關(guān)系的實(shí)證研究[J].工業(yè)技術(shù)經(jīng)濟(jì)，2008，27（12）：81-83.

[14] 龐博慧.中國生產(chǎn)服務(wù)業(yè)與制造業(yè)共生演化模型實(shí)證研究[J].中國管理科學(xué)，2012，20（2）：176-183.

[15] 金晟.生產(chǎn)性服務(wù)業(yè)與制造業(yè)共生演化動(dòng)力機(jī)理探討[J].統(tǒng)計(jì)與決策，2018，34（9）：59-61.

[16] 于斌斌，胡漢輝.產(chǎn)業(yè)集群與城市化共生演化的機(jī)制與路徑——基于制造業(yè)與服務(wù)業(yè)互動(dòng)關(guān)系的視角[J].科學(xué)學(xué)與科學(xué)技術(shù)管理，2014，35（3）：58-68.

[17] 吳勇民，紀(jì)玉山，呂永剛.金融產(chǎn)業(yè)與高新技術(shù)產(chǎn)業(yè)的共生演化研究——來自中國的經(jīng)驗(yàn)證據(jù)[J].經(jīng)濟(jì)學(xué)家，2014（7）：82-92.

[18] 劉平峰，張旺.創(chuàng)新生態(tài)系統(tǒng)共生演化機(jī)制研究[J]. 中國科技論壇，2020（2）：17-27.

[19] 范欽滿，周凌云，樊俊杰，等.區(qū)域物流生態(tài)系統(tǒng)協(xié)同演化模型及穩(wěn)定性分析[J].統(tǒng)計(jì)與決策，2019，35（9）：47-51.

[20] 張寶生，張慶普.隱性知識(shí)流轉(zhuǎn)網(wǎng)成員合作的共生關(guān)系、演化模型及仿真研究[J].現(xiàn)代情報(bào)，2020，40（11）：34-43，53.

[21] 苗林棟，潘文卿. 中國三大增長極生產(chǎn)性服務(wù)業(yè)與制造業(yè)的共生關(guān)系比較[J]. 技術(shù)經(jīng)濟(jì)，2014，33（11）：36-43.

[22] 劉軍躍，李軍鋒，鐘升.生產(chǎn)性服務(wù)業(yè)與裝備制造業(yè)共生關(guān)系研究——基于全國31省市的耦合協(xié)調(diào)度分析[J].湖南科技大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2013，16（1）：111-116.

[23] 張秀艷，白雯，鄭雪.我國區(qū)域經(jīng)濟(jì)韌性的關(guān)聯(lián)識(shí)別與演化特征分析[J].吉林大學(xué)社會(huì)科學(xué)學(xué)報(bào)，2021，61（1）：90-101，237.

[24] 戎陸慶，魏鋒. 基于產(chǎn)業(yè)共生視角的農(nóng)業(yè)與物流業(yè)協(xié)同發(fā)展研究[J]. 江蘇農(nóng)業(yè)科學(xué)，2018，46（5）：332-337.

[25] 龐博慧，郭振. 生產(chǎn)性服務(wù)業(yè)和制造業(yè)共生演化模型研究[J]. 經(jīng)濟(jì)管理，2010，32（9）：28-35.

[26] 唐強(qiáng)榮，徐學(xué)軍，何自力. 生產(chǎn)性服務(wù)業(yè)與制造業(yè)共生發(fā)展模型及實(shí)證研究[J]. 南開管理評(píng)論，2009，12（3）：20-26.

（責(zé)編：何 艷）