王勇　劉永　龔科

摘要：針對我國當前農(nóng)村勞動力在地域間大規(guī)模轉(zhuǎn)移與產(chǎn)業(yè)集聚水平差距的擴大，本文通過構(gòu)建基于物流園區(qū)空間吸引的引力模型明確引力模型中的“質(zhì)量”與勞動力轉(zhuǎn)移率的關(guān)系，結(jié)合中西部地區(qū)的城鄉(xiāng)“二元結(jié)構(gòu)”形態(tài)構(gòu)建了勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型，從物流園區(qū)空間吸引的視角分析勞動力轉(zhuǎn)移與產(chǎn)業(yè)集聚效應(yīng)，通過模型仿真分析運輸成本、勞動力轉(zhuǎn)移率、勞動力產(chǎn)品收益間接效用函數(shù)冪指數(shù)、不同地區(qū)間制造業(yè)產(chǎn)品產(chǎn)出比關(guān)系，發(fā)現(xiàn)物流園區(qū)空間吸引力與勞動力轉(zhuǎn)移率之間具有正的相關(guān)性，運輸成本與產(chǎn)業(yè)集聚效應(yīng)之間具有負的相關(guān)性，而勞動力轉(zhuǎn)移率對產(chǎn)業(yè)集聚效應(yīng)和勞動力產(chǎn)品收益具有規(guī)模報酬遞增的特征。

關(guān)鍵詞：物流園區(qū)；空間吸引；勞動力轉(zhuǎn)移；產(chǎn)業(yè)集聚；二元結(jié)構(gòu)

中圖分類號：F0615文獻標識碼：A

隨著我國經(jīng)濟的發(fā)展和區(qū)域物流業(yè)發(fā)展的不斷增強，物流園區(qū)空間吸引范圍不斷外延，地區(qū)間收入差距持續(xù)擴大，勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚現(xiàn)象逐漸凸顯。區(qū)域物流的發(fā)展、特別是物流園區(qū)的建設(shè)促進了區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)變，增強了區(qū)域經(jīng)濟競爭力，但也呈現(xiàn)出因地區(qū)經(jīng)濟差異而導致的物流園區(qū)分布構(gòu)建不合理的現(xiàn)象，區(qū)域間缺乏有效的協(xié)調(diào)制約了區(qū)域內(nèi)生產(chǎn)要素的流通，而物流園區(qū)的空間吸引是造成地區(qū)勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚的主要影響因素之一。針對我國中西部地區(qū)城鄉(xiāng)“二元結(jié)構(gòu)”的經(jīng)濟形態(tài)，本文擬從物流園區(qū)空間吸引的角度分析勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚的特征，試圖分析勞動力在地域間的轉(zhuǎn)移問題和產(chǎn)業(yè)集聚效應(yīng)的形成特點，并試圖驗證勞動力轉(zhuǎn)移率對產(chǎn)業(yè)集聚和勞動力產(chǎn)品收益具有正向規(guī)模報酬遞增作用。

一、基于物流園區(qū)空間吸引的勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型

（一） 變量定義

假定存在兩個地區(qū)（地區(qū)1和地區(qū)2）具有對稱的經(jīng)濟結(jié)構(gòu)，也即地區(qū)1和地區(qū)2具有相似的生產(chǎn)函數(shù)類型。本文以上述兩地區(qū)為研究對象，下面介紹涉及到的一些定義和符號。

定義1：地區(qū)1和地區(qū)2均有城市、鄉(xiāng)鎮(zhèn)和農(nóng)村，且地區(qū)1和地區(qū)2具有V、L和H三種生產(chǎn)要素，其中V為可耕種的土地，L為低技能勞動力，H為高技能勞動力，高技能勞動力可以在不同的地區(qū)之間不受限制轉(zhuǎn)移，而低技能勞動力的轉(zhuǎn)移受到條件限制。

定義2 ：D{Duu=1，2，…，m′}為參與評價的專家數(shù)；U{Ui，t|i=1，2；t=1，2，…，r}為地區(qū)i的評價指標屬性；

Xi，u，t，l（i=1，2；u=1，2，…，m′；t=1，2，…，r；l=1，2，…，ν）為專家u對于地區(qū)i的指標因素t下的實物l的綜合評價值，用模糊數(shù)表示為Xi，u，t，l=（ai，u，t，l，bi，u，t，l，ci，u，t，l），W{Wi，u，t|i=1，2；u=1，2，…，m′；t=1，2，…，r}為專家u對地區(qū)i的指標因素t的評價值，并用模糊數(shù)表示為Wi，u，t=（hi，u，t，vi，u，t，si，u，t），也被稱作模糊屬性權(quán)重，Zi，l為多位專家對地區(qū)i的實物l的綜合評價值，用模糊數(shù)表示為：Zi，l=（Ti，l，Oi，l，Hi，l），ωi，l為Zi，l的集成指標函數(shù)值，τ為勞動力轉(zhuǎn)移率的阻尼系數(shù)。

定義3：地區(qū)1和地區(qū)2的三種生產(chǎn)要素歸屬于農(nóng)村部門、鄉(xiāng)鎮(zhèn)企業(yè)部門、物流生產(chǎn)制造業(yè)部門和物流服務(wù)業(yè)部門，其中F為農(nóng)村部門生產(chǎn)的產(chǎn)品，E為鄉(xiāng)鎮(zhèn)企業(yè)生產(chǎn)的產(chǎn)品，Q為物流生產(chǎn)制造業(yè)部門的產(chǎn)品，S為物流服務(wù)業(yè)部門生產(chǎn)的產(chǎn)品，且F型產(chǎn)品和E型產(chǎn)品來自農(nóng)村和鄉(xiāng)鎮(zhèn)，Q型產(chǎn)品和S型產(chǎn)品來自城市。

定義4： d{di|i=1，2}表示地區(qū)i的運輸成本，而1/di表示地區(qū)i的運輸產(chǎn)品的到達率，運輸成本越小越容易到達目的地，運輸產(chǎn)品的價格也就越低，產(chǎn)品的收益也就越高。因此，勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚與運輸成本密切相關(guān)。

定義5 ：F{Fi|i=1，2}表示地區(qū)i的F型產(chǎn)品的產(chǎn)出，LFi表示地區(qū)i利用的低技能勞動力的數(shù)量，Ci表示地區(qū)i的可耕種土地面積。

定義6 ：E{Ei|i=1，2}表示地區(qū)i的E型產(chǎn)品的數(shù)量，Q{Qi|i=1，2}表示地區(qū)i的Q型產(chǎn)品的數(shù)量，S{Si|i=1，2}表示地區(qū)i的S型產(chǎn)品的數(shù)量，LEi表示地區(qū)i的E型產(chǎn)品雇傭的勞動力數(shù)量，LSi表示地區(qū)i的S型產(chǎn)品生產(chǎn)部門雇傭的勞動力數(shù)量，LSi0表示地區(qū)i的城市低技能勞動力初始數(shù)量，ΔLSi表示地區(qū)i的低技能勞動力潛在轉(zhuǎn)移量。由于低技能勞動力轉(zhuǎn)移受到城市相關(guān)條件限制而導致僅有一定比例的勞動力可以從農(nóng)村和鄉(xiāng)鎮(zhèn)轉(zhuǎn)移到城市，設(shè)定λi表示地區(qū)i基于低技能勞動力潛在轉(zhuǎn)移量的轉(zhuǎn)移率（0λi1）。

定義 7： PFi表示地區(qū)i的F型產(chǎn)品價格，PEi表示地區(qū)i的E型產(chǎn)品價格，PQi表示地區(qū)i的Q型產(chǎn)品價格，PSi表示地區(qū)i的S型產(chǎn)品價格，wi表示地區(qū)i的勞動力產(chǎn)品收益，ei表示地區(qū)i的低技能勞動力與勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)。

定義 8 ：Qi，j表示地區(qū)i的第j種Q型產(chǎn)品的生產(chǎn)數(shù)量，Ni表示地區(qū)i的Q型產(chǎn)品的種類，f表示用于Q型產(chǎn)品生產(chǎn)的固定投入，c表示Q型產(chǎn)品生產(chǎn)的邊際投入，gi，j表示地區(qū)i的第j種Q型廠商所雇傭的高技能勞動力數(shù)量。

定義 9 ：γ表示F型產(chǎn)品生產(chǎn)函數(shù)的彈性系數(shù)，ζ表示E型產(chǎn)品的生產(chǎn)相關(guān)系數(shù)，ψi表示地區(qū)i勞動力潛在轉(zhuǎn)移量與勞動力產(chǎn)品收益的相關(guān)系數(shù)，φ表示不同種類的Q型產(chǎn)品生產(chǎn)函數(shù)的替代彈性系數(shù)，σ和β表示S型產(chǎn)品生產(chǎn)函數(shù)的彈性系數(shù)。

（二）基于物流園區(qū)空間吸引的引力模型

物流園區(qū)的空間吸引可產(chǎn)生一定的輻射力，并促進城市在發(fā)展過程中逐漸吸引勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚[1-2]。物流園區(qū)的空間吸引力與其對周圍地區(qū)的吸引強度成正比，與區(qū)域之間的距離成反比。引力模型是用于預測和分析事物空間相互作用形式的數(shù)學模型[3-5]，常用于分析區(qū)域間的信息流、貨物流和人員流動等。本文利用引力模型原理計算城市物流園區(qū)的引力大小，構(gòu)建勞動力轉(zhuǎn)移概率模型，進而構(gòu)建勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型。

假設(shè)在地區(qū)1和地區(qū)2的城市均有物流園區(qū)，而農(nóng)村和城鎮(zhèn)可看做一個整體，簡稱農(nóng)鎮(zhèn)，其沒有物流園區(qū)但擁有可耕種的土地和大量低技能勞動力，物流園區(qū)所在地區(qū)i的城市“質(zhì)量”可表示為mi1，此處的mi1表示為地區(qū)i擁有物流園區(qū)的城市的整體經(jīng)濟實力，而地區(qū)i的農(nóng)鎮(zhèn)的質(zhì)量可表示為mi2，此處的mi2表示為地區(qū)i的農(nóng)村和城鎮(zhèn)的經(jīng)濟基礎(chǔ)，經(jīng)濟實力和經(jīng)濟基礎(chǔ)主要是由交通便利程度、經(jīng)濟社會發(fā)展環(huán)境、地理區(qū)位條件和服務(wù)水平等因素共同組成。若地區(qū)i的城市與農(nóng)鎮(zhèn)的中間地帶有j點（j∈J），其質(zhì)量可表示為mij，其中J表示為城市和農(nóng)鎮(zhèn)之間輻射范圍的邊界線，則地區(qū)i的城市和農(nóng)鎮(zhèn)對j點的吸引力分別為：

Fi1j=kimi1mijli1j2Fi2j=kimi2mijli2j2（1）

式中Fi1j，F(xiàn)i2j表示為地區(qū)i的城市和農(nóng)鎮(zhèn)對j點的吸引力，li1j，li2j表示為地區(qū)i的城市和農(nóng)鎮(zhèn)到j(luò)點的距離，ki為系數(shù)。將地區(qū)i的城市和農(nóng)鎮(zhèn)看作兩點，它們的連線為x軸，兩點之間的直線距離設(shè)為2a，兩點連線的垂直平分線作y軸，交點O取為原點，設(shè)j點的坐標為（x，y），則有：

Fi1j=Fi2jmi1li1j2=mi2li2j2（2）

由上式計算得到：

mi1（（x-a）2+y2）=mi2（（x+a）2+y2）（3）

進一步變形得到：

（x+mi1+mi2mi2-mi1a）2+y2=[（mi1+mi2mi2-mi1）2-1]2 （4）

由上式可知j點軌跡是一個以（-mi1+mi2mi2-mi1a，0）為圓心，以2ami1mi2|mi2-mi1|為半徑的圓。亦即質(zhì)量大的實體點具有較強的空間吸引力，而質(zhì)量小的實體點則具有較弱的空間吸引力。

（三）引力模型中“質(zhì)量”及勞動力轉(zhuǎn)移率的確定

基于引力模型的分析需要確定地區(qū)i的城市質(zhì)量mi1值和農(nóng)鎮(zhèn)質(zhì)量mi2值，本文采用模糊綜合評價方法[6-8]，從交通便利程度（U1）、經(jīng)濟社會發(fā)展環(huán)境（U2）、地理區(qū)位條件（U3）、服務(wù)水平（U4），物流供求狀況（U5）和信息化程度（U6）等指標進行綜合評價?；谏鲜鼍C合評價指標將三角模糊數(shù)和語言變量相結(jié)合（見表1），對不同指標下城市和農(nóng)鎮(zhèn)的滿意度及不同指標選擇的滿意度進行綜合評價。

對上述6個指標進行評價時，語言變量值中的Low被Poor代替，Medium被Fair代替，High被Good代替，表1中的變量縮寫值也作相應(yīng)改變，對應(yīng)的模糊數(shù)值不變。

由定義2可知Xi，u，t，l=（ai，u，t，l，bi，u，t，l，ci，u，t，l），Wi，u，t=（hi，u，t，vi，u，t，si，u，t），則多位專家對于地區(qū)i的實物l在指標因素t上的綜合評價值可表示為：

Zi，l=1m′×r∑m′u=1（（Xi，u，1，lWi，u，1）（Xi，u，2，lWi，u，2）…（Xi，u，r，lWi，u，r））（5）

式中表示向量乘積，表示為向量加和，Zi，t，l表示多位專家對地區(qū)i的實物l在指標因素t上的綜合評價值，m′表示總的專家數(shù)，r表示為總的評價指標數(shù)。

假定R=（b，u，e）表示一個三角模糊數(shù)，這三角模糊數(shù)的集成表示為：P（R）=14（b+2u+e）[6，9]。因此，由定義2可知Zi，l的三角模糊數(shù)表示為Zi，l=（Ti，l，Oi，l，Hi，l），則Zi，l的集成指標函數(shù)值ωi，l可以表示為：

ωi，l=14（Ti，l+2Oi，l+Hi，l）（6）

上式中求得的ωi，l即可對應(yīng)表示地區(qū)i的城市“質(zhì)量”mi1和地區(qū)i的農(nóng)鎮(zhèn)“質(zhì)量”mi2，進而用于引力模型的計算。本文以城市和農(nóng)鎮(zhèn)兩點的連線與軌跡圓的交點j以及阻尼系數(shù)τ來確定勞動力轉(zhuǎn)移率，地區(qū)i的勞動力轉(zhuǎn)移率的計算公式可以表示為：

λi=τi（mi1mi1-mi2-mi1mi2|mi2-mi1|）（7）

（四） 勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型構(gòu)建

基于物流園區(qū)空間吸引的勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型，包含農(nóng)村部門、鄉(xiāng)鎮(zhèn)企業(yè)部門、城市的物流生產(chǎn)制造業(yè)部門和物流服務(wù)業(yè)部門等的生產(chǎn)函數(shù)模型，且生產(chǎn)函數(shù)模型間相互關(guān)聯(lián)和制約[1，3，10]，生產(chǎn)函數(shù)模型依次被介紹如下：

1.農(nóng)村部門的生產(chǎn)函數(shù)模型

農(nóng)村部門利用土地和低技能勞動力生產(chǎn)F型產(chǎn)品，地區(qū)i的F型產(chǎn)品數(shù)量的生產(chǎn)函數(shù)模型為：

Fi=L1-γFiCγi（8）

其中地區(qū)i的F型產(chǎn)品數(shù)量是由地區(qū)i利用的低技能勞動力數(shù)量和地區(qū)i的可耕種土地面積與彈性系數(shù)結(jié)合組成的函數(shù)式。

2.鄉(xiāng)鎮(zhèn)企業(yè)部門的生產(chǎn)函數(shù)模型

鄉(xiāng)鎮(zhèn)企業(yè)生產(chǎn)E型產(chǎn)品，假定E型產(chǎn)品的生產(chǎn)函數(shù)要素只是低技能勞動，則地區(qū)i的E型產(chǎn)品數(shù)量的生產(chǎn)函數(shù)模型為：

Ei=ζLEi（9）

其中地區(qū)i的E型產(chǎn)品數(shù)量是由地區(qū)i的E型產(chǎn)品雇傭的勞動力數(shù)量與生產(chǎn)相關(guān)系數(shù)ζ的乘積組成。

3. 勞動力產(chǎn)品收益的生產(chǎn)函數(shù)模型

假定不同地區(qū)的勞動力產(chǎn)品收益是與F型產(chǎn)品數(shù)量，S型產(chǎn)品數(shù)量，勞動力轉(zhuǎn)移率，運輸成本等相關(guān)的生產(chǎn)函數(shù)，則地區(qū)i的勞動力產(chǎn)品收益生產(chǎn)函數(shù)模型為：

wi=（PQidi）Si2-λiPFiFi2（10）

其中勞動力產(chǎn)品收益的生產(chǎn)函數(shù)模型是F型產(chǎn)品數(shù)量和S型產(chǎn)品數(shù)量的二元二次函數(shù)關(guān)系式。

4.城市低技能勞動力數(shù)量的生產(chǎn)函數(shù)模型

由于受到物流園區(qū)空間吸引的影響，部分低技能勞動力將從農(nóng)村和鄉(xiāng)鎮(zhèn)轉(zhuǎn)移到城市，勞動力轉(zhuǎn)移率是導致城市低技能勞動力數(shù)量變化的主要影響因素，因此地區(qū)i的城市低技能勞動力數(shù)量的生產(chǎn)函數(shù)模型為：

LSi=LSi0+λiΔLSi（11）

其中城市低技能勞動力數(shù)量由初始城市低技能勞動力數(shù)量和低技能勞動力轉(zhuǎn)移量共同組成，城市和農(nóng)村的勞動力產(chǎn)品總收益的不同將導致工資率的差異，進而導致勞動力從農(nóng)村轉(zhuǎn)移到城市。因此，假定低技能勞動力潛在轉(zhuǎn)移量ΔLSi是勞動力產(chǎn)品收益的間接效用函數(shù)，即ΔLSi=ψi|wi|ei。此外，低技能勞動力轉(zhuǎn)移率λi的變化將導致工資率的變化，在勞動力產(chǎn)品總收益一定的情況下，當勞動力轉(zhuǎn)移率λi較小時可導致一個較高的工資率，而當當勞動力轉(zhuǎn)移率λi較大時可導致一個較低的工資率。

5.城市部門相關(guān)產(chǎn)品數(shù)量生產(chǎn)函數(shù)模型

物流生產(chǎn)制造業(yè)部門是一個高技能勞動力密集型部門，高技能勞動力被雇傭用于生產(chǎn)各種類型的Q型產(chǎn)品，這里假定所有的Q型產(chǎn)品具有相同的生產(chǎn)函數(shù)，則地區(qū)i的第j種Q型產(chǎn)品產(chǎn)出數(shù)量的生產(chǎn)函數(shù)模型為：

Qi，j=（gi，j/c）-（f/c） （12）

其中地區(qū)i的第j種Q型產(chǎn)品產(chǎn)出數(shù)量是該部門所雇傭的高技能勞動力數(shù)量，產(chǎn)品生產(chǎn)的固定投入和邊際投入的一階函數(shù)式。根據(jù)Sousa（2012）[11]和Markusen（2013）[12]的研究，地區(qū)i的S型產(chǎn)品產(chǎn)出可以采用以下形式的生產(chǎn)函數(shù)模型：

Si=L1-σ-βeiSi[∑Nij=1Q（φ-1）/φi，j]φ/（φ-1）σEβi（LSi0+λiPFiFi2）β（13）

其中地區(qū)i的S型產(chǎn)品產(chǎn)出是關(guān)于城市低技能勞動力數(shù)量LSi，E型產(chǎn)品數(shù)量，F(xiàn)型產(chǎn)品數(shù)量和各種類型的Q型產(chǎn)品產(chǎn)出的齊次生產(chǎn)函數(shù)。參數(shù)φ表示對不同種類的Q型產(chǎn)品的使用份額，當φ>1時，一定比例的Q型產(chǎn)品使得生產(chǎn)函數(shù)呈現(xiàn)出規(guī)模報酬遞增的特性，而當φ<1時，一定比例的Q型產(chǎn)品使得生產(chǎn)函數(shù)呈現(xiàn)出規(guī)模報酬遞減的特性。

二、仿真計算

（一） 實例相關(guān)數(shù)據(jù)

針對我國中西部地區(qū)典型的城鄉(xiāng)“二元結(jié)構(gòu)”形態(tài)，本文選取西部地區(qū)（如重慶、成都等）包含有對稱經(jīng)濟結(jié)構(gòu)的兩地區(qū)（地區(qū)1和地區(qū)2），進行基于物流園區(qū)空間吸引的勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚特征分析。地區(qū)1和地區(qū)2共同邀請3位專家D={D1，D2，D3}，運用表1中的模糊語言變量值對不同指標下城市和農(nóng)鎮(zhèn)的滿意度及不同指標選擇的滿意度進行綜合評價，見表2至表4。

為了驗證和分析方法求解問題的有效性，本文參考相關(guān)研究文獻[3，10，13-14]對部分參數(shù)進行賦值操作，設(shè)地區(qū)1和地區(qū)2的低技能勞動力數(shù)量為LF1=LF2=1500，地區(qū)1和地區(qū)2的城市初始低技能勞動力數(shù)量分別為LS10=400，LS20=500，高技能勞動力數(shù)量為g1，1=g1，2=g1，3=500，g2，1=g2，2=g2，3=400，地區(qū)1和地區(qū)2可耕種的土地面積為C1=100，C2=200，勞動力轉(zhuǎn)移率的阻尼系數(shù)為τ1=04，τ2=03，地區(qū)1和地區(qū)2的F型產(chǎn)品生產(chǎn)函數(shù)的彈性系數(shù)為γ=015，地區(qū)1和地區(qū)2生產(chǎn)E型產(chǎn)品雇傭的勞動力數(shù)量分別為LE1=400，LE2=500，E型產(chǎn)品的生產(chǎn)相關(guān)系數(shù)為ζ=15，地區(qū)1和地區(qū)2的F型產(chǎn)品和Q型產(chǎn)品的價格分別為PF1=8，PF2=6，PQ1=5，PQ2=4，運輸成本分別為d1=10，d2=8，地區(qū)1和地區(qū)2的勞動力潛在轉(zhuǎn)移量與勞動力產(chǎn)品收益的相關(guān)系數(shù)為ψ1=ψ2=02，地區(qū)1和地區(qū)2的低技能勞動力與勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)為e1=e2=05，地區(qū)1和地區(qū)2的S型產(chǎn)品生產(chǎn)函數(shù)的彈性系數(shù)為σ=β=05，地區(qū)1和地區(qū)2用于Q型產(chǎn)品生產(chǎn)的固定投入和Q型產(chǎn)品生產(chǎn)的邊際投入為f=4，c=2，地區(qū)1和地區(qū)2的Q型產(chǎn)品種類分別為N1=N2=3，地區(qū)1和地區(qū)2的不同種類的Q型產(chǎn)品的使用份額為φ=168。

（二）計算結(jié)果及分析

基于相關(guān)變量的定義，本文應(yīng)用物流園區(qū)空間吸引的引力模型計算得到地區(qū)1的城市質(zhì)量m11和農(nóng)鎮(zhèn)質(zhì)量m12以及地區(qū)2的城市質(zhì)量m21和農(nóng)鎮(zhèn)質(zhì)量m22分別為m11=0648，m12=0498，m21=0574，m22=0467，進而應(yīng)用公式（7）計算得到地區(qū)1和地區(qū)2的勞動力轉(zhuǎn)移率分別為λ1=021，λ2=016。基于計算得到的勞動力轉(zhuǎn)移率和相關(guān)變量的賦值定義，應(yīng)用勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型進行計算，得到如下變量數(shù)值：

地區(qū)1和地區(qū)2的E型產(chǎn)品數(shù)量為：E1=600，E2=750；

地區(qū)1和地區(qū)2的F型產(chǎn)品數(shù)量為：F1=999，F(xiàn)2=1109；

地區(qū)1和地區(qū)2的S型產(chǎn)品產(chǎn)出為：S1=307，S2=76；

地區(qū)1的N1種Q型產(chǎn)品產(chǎn)出數(shù)量為：Q1，1=Q1，2=Q1，3=246；

地區(qū)1的N1種Q型產(chǎn)品產(chǎn)出數(shù)量為：Q2，1=Q2，2=Q2，3=196；

地區(qū)1和地區(qū)2勞動力轉(zhuǎn)移后城市低技能勞動力數(shù)量為：LS1=560，LS2=686。

通過上述計算數(shù)值可發(fā)現(xiàn)地區(qū)1的S型和Q型產(chǎn)品的產(chǎn)出高于地區(qū)2的相應(yīng)產(chǎn)品的產(chǎn)出，而S型和Q型產(chǎn)品的產(chǎn)出規(guī)模是衡量地區(qū)物流發(fā)展水平和產(chǎn)業(yè)集聚程度的關(guān)鍵，進而也是地區(qū)整體經(jīng)濟發(fā)展程度的充分體現(xiàn)。此外，農(nóng)鎮(zhèn)和城市的差別導致農(nóng)鎮(zhèn)勞動力出現(xiàn)轉(zhuǎn)移，且不同地區(qū)的經(jīng)濟發(fā)展差異導致勞動力的轉(zhuǎn)移亦有差別，如地區(qū)1城市低技能勞動力的數(shù)量由初始的400人增加為當前的560人，而地區(qū)2的城市低技能勞動力數(shù)量則由初始的500人增加為當前的686人。

勞動力轉(zhuǎn)移和產(chǎn)業(yè)集聚模型中已經(jīng)充分考慮到受物流園區(qū)空間吸引的影響以及勞動力產(chǎn)品收益情況，農(nóng)鎮(zhèn)將出現(xiàn)向城市轉(zhuǎn)移勞動力，城市將出現(xiàn)產(chǎn)業(yè)集聚的現(xiàn)象，且受到地區(qū)經(jīng)濟發(fā)展水平和條件限制等的影響，勞動力轉(zhuǎn)移率與產(chǎn)業(yè)集聚密切相關(guān)，且勞動轉(zhuǎn)移率的改變對產(chǎn)業(yè)集聚效應(yīng)和地區(qū)差異也將產(chǎn)生一定的影響。

三、敏感度分析

勞動力轉(zhuǎn)移與產(chǎn)業(yè)集聚之間存在著一定的相關(guān)性，在不同運輸成本條件下，勞動力轉(zhuǎn)移率λ的改變將對產(chǎn)業(yè)集聚產(chǎn)生一定的影響。在設(shè)定了有關(guān)參數(shù)的初始值后分析勞動力轉(zhuǎn)移率λ的改變與地區(qū)1和地區(qū)2的S型產(chǎn)品的比率S1/S2的關(guān)系，將有助于探討不同地區(qū)間物流園區(qū)空間吸引的勞動力轉(zhuǎn)移率變化對產(chǎn)業(yè)集聚的影響。如圖1中（a）和（b）所示，x軸表示勞動力轉(zhuǎn)移率λ，y軸表示地區(qū)1和地區(qū)2的S型產(chǎn)品的產(chǎn)出比S1/S2，勞動力轉(zhuǎn)移率從0逐漸增大到1，運輸成本分為四種情況設(shè)定：（1）d1=24，d2=12；（2）d1=27，d2=9；（3）d1=40，d2=8；（4）d1=45，d2=5。

在上述4種運輸成本情況下，圖1（a）中S1的勞動力轉(zhuǎn)移率固定，隨著S2中勞動力轉(zhuǎn)移率的增加，地區(qū)1和地區(qū)2的S型產(chǎn)品的產(chǎn)出比S1/S2呈現(xiàn)緩慢下降的趨勢，這說明了勞動力從農(nóng)鎮(zhèn)轉(zhuǎn)移到城市將造成物流生產(chǎn)制造業(yè)向城市集聚；隨著S2的勞動力轉(zhuǎn)移率增加，集聚效應(yīng)逐漸凸顯，且地區(qū)1和地區(qū)2之間經(jīng)濟發(fā)展水平存在差異（地區(qū)1由于地區(qū)2），這也是導致S1/S2呈現(xiàn)緩慢下降的一個原因。在圖1（b）中，S2中勞動力轉(zhuǎn)移率固定，隨著S1中勞動力轉(zhuǎn)移率的逐漸增加，地區(qū)1和地區(qū)2的S型產(chǎn)品的產(chǎn)出比S1/S2呈現(xiàn)顯著上升趨勢，這同樣說明地區(qū)1的勞動力從農(nóng)鎮(zhèn)轉(zhuǎn)移到城市造成了物流制造業(yè)向城市集聚，而且從曲線的顯著上升程度可以看出地區(qū)1的集聚效應(yīng)總體上優(yōu)于地區(qū)2的集聚效應(yīng)。

此外，隨著勞動力轉(zhuǎn)移率的增加，在圖1的（a）和（b）可以看出地區(qū)1和地區(qū)2的物流運輸成本比率d1/d2越大，產(chǎn)業(yè)集聚效應(yīng)也就越顯著（參見圖1（b）），特別是考慮到地區(qū)1和地區(qū)2之間存在經(jīng)濟差異，產(chǎn)業(yè)集聚效應(yīng)越低，也就越容易導致勞動力的轉(zhuǎn)移，進而造成產(chǎn)業(yè)集聚的萎縮，由此，地區(qū)間勞動力轉(zhuǎn)移的現(xiàn)象也將出現(xiàn)。

勞動力轉(zhuǎn)移與地區(qū)勞動力產(chǎn)品收益也存在一定的相關(guān)性，在設(shè)定不同的地區(qū)勞動力產(chǎn)品收益間接效用函數(shù)冪指數(shù)時，勞動力轉(zhuǎn)移率λ的改變將對勞動力轉(zhuǎn)移的集聚效應(yīng)產(chǎn)生一定的影響。如見圖2所示，x軸表示勞動力轉(zhuǎn)移率λ，y軸表示地區(qū)1和地區(qū)2的城市低技能勞動力數(shù)量比率LS1/LS2，勞動力轉(zhuǎn)移率從0逐漸增大到1，設(shè)定地區(qū)1和地區(qū)2的低技能勞動力與勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)分別為：（1）e1=e2=08；（2）e1=e2=06。

隨著勞動力轉(zhuǎn)移率的增加，由圖2可見地區(qū)1和地區(qū)2的城市低技能勞動力比率將顯著增長，不同地區(qū)的低技能勞動力與勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)也將影響城市低技能勞動力比率的增長幅度，圖2中LS1/LS2曲線具有遞增的正的斜率，且勞動力轉(zhuǎn)移與城市低技能勞動力數(shù)量之間的關(guān)系呈正相關(guān)性，也即隨著勞動力轉(zhuǎn)移率的提高，地區(qū)間的勞動力數(shù)量變化呈現(xiàn)擴大的趨勢。

綜上所述，運輸條件的優(yōu)越性可為物流園區(qū)的空間吸引提供必要的輔助條件，由于物流發(fā)達的地區(qū)具有較好的地理區(qū)位優(yōu)勢，運輸成本則相應(yīng)較低，進而物流制造業(yè)的成本降低，可以說合理的運輸成本是勞動力轉(zhuǎn)移的一個動力和產(chǎn)業(yè)集聚的一個引力。此外，勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)是與產(chǎn)品收益密切相關(guān)的變量，高的冪指數(shù)將導致勞動力轉(zhuǎn)移的增加，進而增加了勞動力供給，這相應(yīng)的削弱了產(chǎn)業(yè)集聚帶來的工資率上升的壓力，從而成為產(chǎn)業(yè)集聚的另一個引力。

四、結(jié)論

本文針對兩個均有農(nóng)鎮(zhèn)和城市且具有對稱經(jīng)濟結(jié)構(gòu)和相似生產(chǎn)函數(shù)類型的地區(qū)展開研究，分析物流園區(qū)空間吸引對勞動力轉(zhuǎn)移率和產(chǎn)業(yè)集聚的影響。結(jié)果表明勞動力轉(zhuǎn)移率與產(chǎn)業(yè)集聚密切相關(guān)，且勞動轉(zhuǎn)移率的改變對產(chǎn)業(yè)集聚效應(yīng)和地區(qū)差異會產(chǎn)生一定的影響；合理的運輸成本是地區(qū)間與同一地區(qū)勞動力轉(zhuǎn)移的動力和產(chǎn)業(yè)集聚的引力，高的勞動力產(chǎn)品收益間接效用函數(shù)的冪指數(shù)也是產(chǎn)業(yè)集聚形成的一個引力。因此，只有加快發(fā)展不同地區(qū)間及同一地區(qū)的城市與鄉(xiāng)鎮(zhèn)產(chǎn)品的物流制造業(yè)和服務(wù)業(yè)，才能提高低技能勞動力的產(chǎn)品收益，進而縮小城鄉(xiāng)二元經(jīng)濟結(jié)構(gòu)造成的地區(qū)差異，防止經(jīng)濟邊緣化現(xiàn)象的出現(xiàn)。

參考文獻：

[1]劉紅梅，李國軍，王克強. 中國農(nóng)業(yè)虛擬水國際貿(mào)易影響因素研究——基于引力模型的分析[J].管理世界，2010（9）：76-87.

[2]楊勝剛，楊建模. 城市化與第三產(chǎn)業(yè)發(fā)展：基于湖南的實證分析[J].商業(yè)研究，2010（3）： 143-147.

[3]李揚. 西部地區(qū)產(chǎn)業(yè)集聚水平測度的實證研究[J].南開經(jīng)濟研究，2009（4）：144-151.

[4]海峰，閻欣，丁燦，邵校. 基于軸輻理論的現(xiàn)代區(qū)域物流網(wǎng)絡(luò)節(jié)點選擇[J].計算機集成制造系統(tǒng)，2012，18（6）：1299-1305.

[5]Nuzzolo，A.，Comi，A. City logistics planning： demand modeling requirements for direct effect forecasting[J].Procedia-Social and Behavioral Science，2014（125）：239-250.

[6]Chou，S. Y.，Chang，Y. H.，Shen，C. Y. A fuzzy simple additive weighting system under group decision making for facility location selection with objective-subjective attributes[J].European Journal of Operational Research，2008，189（1）：132-145.

[7] Wang，Y.，Ma，X.L.，Wang，Y.H.，Mao，H.J.，Zhang，Y. Location optimization of multiple distribution centers under fuzzy environment[J].Journal of Zhejiang University SCIENCE A，2012，13（10）：782-798.

[8]Wang，Y.，Ma，X.L.，Lao，Y.T.，Wang，Y.H. A Fuzzy-based Customer Clustering Approach with Hierarchical Structure for Logistics Network Optimization[J].Expert Systems with Applications，2014，41（2）：521-534.

[9]Kahraman，C.，Ruan，D.，Dogan，I. Fuzzy group decision-making for facility location selection[J].Information Sciences，2003（157）：135-153.

[10]敖榮軍. 制造業(yè)集中、勞動力流動與中部地區(qū)的邊緣化[J].南開經(jīng)濟研究，2005（1）：61-66.

[11]Sousa，J. D.，Mayer，T.， Zignago，S. Market access in global and regional trade[J].Regional Science and Urban Economics，2012，42（6）：1037-1052.

[12]Markusen，J. R. Putting per-capita income back into trade theory[J].Journal of International Economics，2013，90（2）：255-265.

[13]范劍勇，王立軍，沈林潔. 產(chǎn)業(yè)集聚與農(nóng)村勞動力的跨區(qū)域流動[J].管理世界，2004（4）：22-29.

[14]姚林如，李莉. 勞動力轉(zhuǎn)移、產(chǎn)業(yè)集聚與地區(qū)差距[J].財經(jīng)研究，2006，32（8）：135-143.

Abstract： In view of large-scale rural labor interregional migration and the expansion of industrial agglomeration level gap， this paper constructs the gravity model based on logistics park space attraction to determine the relationship between the “quality” in gravity model and labor migration rate， and considers the urban and rural “dual structure” of the Midwestern in China to construct the labor migration and industrial agglomeration models， and studies the labor migration and industrial agglomeration from the perspective of the logistics park space attraction. The relationship among the transportation cost， labor migration rate， labor production earnings indirect utility function power exponent， manufacturing production output ratio between different regions are studied through model simulation， which shows that the relationship between logistics park space attraction and labor migration rate is a positive correlation， and the relationship between transportation cost and industrial agglomeration effect is a negative correlation， and the industrial agglomeration effect and labor production earnings exist the characteristics of increasing returns to scale when labor migration rate increases.

Key words：logistics park； space attraction； labor migration； industrial agglomeration； dual structure

（責任編輯：厲新）