倪紅紅，馬 強，卜元坤，楊小軒，李衛(wèi)忠

（1.西北農(nóng)林科技大學(xué)林學(xué)院，陜西 楊凌 712100；2.安徽科技學(xué)院資源與環(huán)境學(xué)院，安徽 滁州 233100）

林業(yè)在國民經(jīng)濟及社會可持續(xù)發(fā)展戰(zhàn)略中占據(jù)極其重要的地位［1］，具有經(jīng)濟、社會、生態(tài)和文化等價值。建國以來，我國林業(yè)產(chǎn)值快速增長，由1949年的2.39×109元上升至2020年的7.55×1012元［2］。林業(yè)也是規(guī)模最大的綠色經(jīng)濟體，在拉動內(nèi)需、推動經(jīng)濟綠色增長、積累生態(tài)資產(chǎn)方面具有重大作用。作為林業(yè)的主體，森林生態(tài)系統(tǒng)是最大的陸地儲碳庫，具有很強的儲碳能力，為我國應(yīng)對氣候變化提供了戰(zhàn)略支撐［3-4］。因此，林業(yè)建設(shè)是事關(guān)經(jīng)濟社會可持續(xù)發(fā)展的根本性問題，發(fā)展林業(yè)不僅有巨大的經(jīng)濟利益，也有巨大的生態(tài)功能和社會功能［5-6］。

林業(yè)產(chǎn)業(yè)作為林業(yè)經(jīng)濟發(fā)展的主體，其類型多樣，除了直接從事林業(yè)生產(chǎn)、經(jīng)營管理以外，還包括森林康養(yǎng)、生態(tài)旅游、森林培育、木材加工及木制產(chǎn)品制造、林業(yè)專業(yè)技術(shù)服務(wù)和其他相關(guān)活動等。林業(yè)企業(yè)則是以森林資源、環(huán)境及其產(chǎn)品為主要經(jīng)營對象的企業(yè)，其經(jīng)營活動對經(jīng)濟、生態(tài)和社會等方面都將產(chǎn)生積極影響［7］。一方面可提供社會需要的產(chǎn)品、實現(xiàn)企業(yè)健康可持續(xù)發(fā)展，另一方面，在促進經(jīng)濟發(fā)展和生態(tài)環(huán)境建設(shè)方面也發(fā)揮著重要作用［8］。林業(yè)企業(yè)發(fā)展狀況關(guān)系到地區(qū)經(jīng)濟增長和生態(tài)建設(shè)，對林業(yè)企業(yè)空間格局演變研究及其影響因素分析尤為重要。

關(guān)于林業(yè)企業(yè)的研究主要集中在林業(yè)經(jīng)營與管理、林業(yè)企業(yè)社會價值、林業(yè)產(chǎn)業(yè)集聚等方面。在林業(yè)經(jīng)營與管理方面，諸多學(xué)者就林業(yè)企業(yè)資本結(jié)構(gòu)、績效、杠桿率、稅收政策等［9］進行了研究，提出了若干從企業(yè)結(jié)構(gòu)上優(yōu)化發(fā)展過程的建議。在林業(yè)企業(yè)社會價值方面，學(xué)者們通過對林業(yè)企業(yè)碳匯項目［10］、生態(tài)建設(shè)［11］、社會責任［12］等的研究，認為林業(yè)企業(yè)具有區(qū)別于一般企業(yè)的特性，即在綠色可持續(xù)發(fā)展方面承擔著特殊的社會責任。此外，學(xué)者們就林業(yè)產(chǎn)業(yè)集聚狀況和相關(guān)因素進行了分析和探討。例如，Munnich 等［13］運用波特的“鉆石模型”產(chǎn)業(yè)集群理論分析了美國明尼蘇達州東北部地區(qū)林業(yè)產(chǎn)業(yè)集聚狀況以及該地所存在的競爭優(yōu)勢。Todd 等［14］基于勞動力、運輸基礎(chǔ)設(shè)施、原材料等因素對密西西比州林業(yè)企業(yè)進行空間分析，以確定潛在的林業(yè)企業(yè)集群。韓瑩［15］以全國30 個省份的面板數(shù)據(jù)，探討了我國林業(yè)產(chǎn)業(yè)集聚時空分異，并使用空間計量模型對林業(yè)產(chǎn)業(yè)經(jīng)濟效應(yīng)進行了研究，提出了林業(yè)產(chǎn)業(yè)的發(fā)展對策。

目前，企業(yè)空間格局研究是經(jīng)濟學(xué)、地理學(xué)及管理學(xué)中的研究熱點之一［16］，主要集中在科技型企業(yè)［17-18］、服務(wù)型企業(yè)［19］、創(chuàng)新型企業(yè)［20-21］等，大量學(xué)者利用GIS空間分析技術(shù)探究了企業(yè)的空間演變規(guī)律，使用回歸模型探究相關(guān)影響因素的作用，并提出優(yōu)化資源配置的方案［22-24］。對于農(nóng)林類企業(yè)等以第一產(chǎn)業(yè)為主，第二、第三產(chǎn)業(yè)為輔的行業(yè)，企業(yè)空間格局研究比較少。熊友云等［25］基于相關(guān)分析、生態(tài)分布格局理論，分析了中國農(nóng)業(yè)龍頭企業(yè)在宏觀、微觀地域及分行業(yè)空間布局特征。白如山等［26］運用最近鄰點距離指數(shù)、核密度估計等現(xiàn)代地學(xué)方法，從類型-空間二維，探究阜陽市農(nóng)業(yè)龍頭企業(yè)空間集聚格局及影響因素。蔣輝等［27］探究中國農(nóng)業(yè)產(chǎn)業(yè)化龍頭企業(yè)的空間分布特征，并通過多元線性回歸分析影響農(nóng)業(yè)龍頭企業(yè)空間集聚的因素。

上述研究表明，當前企業(yè)布局研究傾向于利用區(qū)位理論探究制造業(yè)、科技企業(yè)等第二、第三產(chǎn)業(yè)的空間格局，忽略了對于農(nóng)業(yè)、林業(yè)等第一產(chǎn)業(yè)布局的研究。而在以林業(yè)企業(yè)為主體的研究中，多關(guān)注林業(yè)企業(yè)的經(jīng)營與管理［28］以及林業(yè)企業(yè)的社會價值［29］。對于林業(yè)企業(yè)空間布局演變及影響因素的研究較少，且主要基于面板數(shù)據(jù)，從宏觀角度使用統(tǒng)計學(xué)方法對產(chǎn)業(yè)集聚進行分析。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

陜西省既是“一帶一路”建設(shè)的重要樞紐，也是新時代西部大開發(fā)、黃河流域生態(tài)保護和高質(zhì)量發(fā)展等國家重大發(fā)展戰(zhàn)略的重要區(qū)域?！笆濉逼陂g，陜西省林業(yè)產(chǎn)業(yè)有序健康發(fā)展，林業(yè)產(chǎn)業(yè)總產(chǎn)值達到1.48×1011元，從業(yè)人數(shù)超過3×106人。依托森林資源稟賦，發(fā)展綠色富民產(chǎn)業(yè)，形成了以綠色產(chǎn)業(yè)促生態(tài)建設(shè)、以生態(tài)建設(shè)帶動脫貧攻堅的生態(tài)經(jīng)濟體系，生態(tài)旅游、森林體驗、森林康養(yǎng)等服務(wù)型產(chǎn)業(yè)蓬勃發(fā)展。同時，陜西省林業(yè)產(chǎn)業(yè)也存在三大產(chǎn)業(yè)融合不夠緊密，產(chǎn)業(yè)鏈尚未完全貫通，企業(yè)單體規(guī)模偏小，市場競爭能力和抗風險能力差等問題［30］。

本研究獲取陜西省行政區(qū)內(nèi)的林業(yè)企業(yè)10773家，企業(yè)成立時間截止2020 年12 月31 日（圖1），并選取2000、2005、2010、2015、2020 年5 個時間節(jié)點研究其空間演化過程。由于企業(yè)地址為文本數(shù)據(jù)，利用高德API 進行批量地址解析，通過“地址”字段獲得相應(yīng)的經(jīng)緯度坐標，建立林業(yè)企業(yè)空間數(shù)據(jù)庫。

圖1 陜西省林業(yè)企業(yè)分布示意圖Fig.1 Distribution diagram of forestry enterprises in Shaanxi Province

1.2 數(shù)據(jù)來源及篩選

本文以陜西省2020 年林業(yè)相關(guān)企業(yè)為研究對象，企業(yè)樣本數(shù)據(jù)主要來源于全國商業(yè)查詢平臺“啟信寶”，通過對企業(yè)名、經(jīng)營范圍進行“林業(yè)產(chǎn)業(yè)”“林業(yè)企業(yè)”“林業(yè)”關(guān)鍵詞的篩選，獲取陜西省內(nèi)林業(yè)企業(yè)數(shù)據(jù)，該數(shù)據(jù)包括“企業(yè)名稱”“登記狀態(tài)”“法定代表人”“注冊時間”“注冊資本”“地址”“企業(yè)類型”“電話”“郵箱”等。然后對林業(yè)企業(yè)原始數(shù)據(jù)進行清洗，剔除企業(yè)狀態(tài)為“撤銷”“注銷”“吊銷”及基本信息不明確的企業(yè)。由于使用關(guān)鍵詞爬取會采集到一部分與林業(yè)企業(yè)無關(guān)，但名稱相近的企業(yè)數(shù)據(jù)，通過人工識別進一步清洗數(shù)據(jù)，基于企業(yè)名稱及類別標簽進行語義分析，判斷企業(yè)的主營業(yè)務(wù)，剔除不相關(guān)的企業(yè)數(shù)據(jù)。

影響因素數(shù)據(jù)中，2020年公路基礎(chǔ)數(shù)據(jù)來源于OpenStreetMap 志愿者地理信息網(wǎng)站（http://download.geofabrik.de/），并通過投影、分割、計算獲得公路密度；2020年常住人口數(shù)量、地區(qū)生產(chǎn)總值、第一產(chǎn)業(yè)增加值占比、社會消費品零售總額數(shù)據(jù)基于《中國縣域統(tǒng)計年鑒（縣市卷）（2021年）》整理獲得；2020 年企業(yè)注冊資本平均值數(shù)據(jù)基于林業(yè)企業(yè)的屬性數(shù)據(jù)計算整理獲得；2020 年林地面積、園地面積數(shù)據(jù)來源于中國國土調(diào)查成果共享應(yīng)用服務(wù)平臺（https://gtdc.mnr.gov.cn/）。

1.3 研究方法

1.3.1 平均最近鄰分析法平均最近鄰分析通過計算企業(yè)之間的平均距離，表現(xiàn)林業(yè)企業(yè)之間的臨近程度，從而刻畫出林業(yè)企業(yè)數(shù)據(jù)點的集聚程度。最近鄰指數(shù)（ANN）的計算公式如下：

1.3.2 標準差橢圓標準差橢圓及重心可以有效地反映林業(yè)企業(yè)分布離散程度。橢圓的大小反映空間格局總體要素的集中程度，長半軸反映分布的主導(dǎo)方向。標準差橢圓的公式如下：

式中：xi和yi為每個要素的空間位置坐標；(,)為要素的算數(shù)平均中心；n為要素總數(shù)；SDEx和SDEy為林業(yè)企業(yè)標準差橢圓重心z的經(jīng)緯度坐標。

1.3.3 核密度分析法核密度分析用于計算要素在其周圍鄰域中的密度，可以直觀的反映林業(yè)企業(yè)的集聚-擴散程度。計算公式如下：

式中：f(x)為核密度估計值；n為林業(yè)企業(yè)個數(shù)；h為帶寬，h＞0；k為核函數(shù)；(x-xi)為估計點x到樣本點xi的歐式距離。f(x)越高說明企業(yè)集聚程度越高，反之就越低。

1.3.4 普通最小二乘法（OLS）模型與地理加權(quán)回歸（GWR）模型多變量線性回歸分析包括全局線性回歸分析與局部線性回歸分析2 種主要類型，全局線性回歸方法中最傳統(tǒng)的模型為建立在最小二乘法基礎(chǔ)上的回歸模型——OLS模型：

式中：y為觀測的因變量；β0為常數(shù)項估計值；xk為第k個自變量的值；βk為第k個自變量的系數(shù)估計值；n為自變量個數(shù)；ε為誤差項。

OLS 模型只是對樣本和變量進行全局的估計，沒有考慮空間相關(guān)性的影響。GWR 模型是對線性回歸模型的一個拓展，可以對具有空間自相關(guān)的樣本進行很好地擬合，反映其空間特征。GWR模型設(shè)定如下：

式中：yi為n×1 維被解釋變量；xik為n×k維解釋變量矩陣；βk(ui,vi)為因素k在回歸點i的回歸系數(shù)；k為自變量個數(shù)；(ui,vi)為第i個觀察點的經(jīng)緯度坐標；εi為分布的隨機誤差項。

2 結(jié)果與分析

2.1 林業(yè)企業(yè)時空格局演變分析

2.1.1 陜西省林業(yè)企業(yè)平均最近鄰分析2000—2020年陜西省林業(yè)企業(yè)最近鄰指數(shù)均小于1（表1），且逐漸減小，表明陜西省林業(yè)企業(yè)空間分布呈現(xiàn)集聚的特征，且空間集聚態(tài)勢不斷加強。2000—2015年中國林業(yè)企業(yè)處于初步發(fā)展階段，由于陜西省最先頒布封山禁牧的政策，同時得益于豐厚的森林自然資源，形成了規(guī)模較大的林業(yè)產(chǎn)業(yè)集聚中心。林業(yè)企業(yè)集聚分布能夠有效整合人力、資本、技術(shù)等條件要素，建設(shè)符合地方特色的產(chǎn)業(yè)集群，打造規(guī)模效應(yīng)，促進產(chǎn)業(yè)效率提升和區(qū)域經(jīng)濟發(fā)展。2015年后，中國林業(yè)產(chǎn)業(yè)進入快速發(fā)展階段，人口集聚、技術(shù)的進步、自然資源的發(fā)展、配套設(shè)施建設(shè)、政策支持等多項利好因素，拉動更多林業(yè)企業(yè)空間分布向集聚趨勢發(fā)展。

表1 2000—2020年陜西省林業(yè)企業(yè)平均最近鄰分析結(jié)果Tab.1 Average nearest neighbor analysis results of forestry enterprises in Shaanxi Province during 2000—2020

2.1.2 陜西省林業(yè)企業(yè)時空分布標準差橢圓通過不同時間節(jié)點標準差橢圓來描繪陜西省林業(yè)企業(yè)空間方向差異及重心變化（圖2）。（1）在橢圓重心方面，重心軌跡變化為先向東北遷移，后往東南遷移，但重心均位于咸陽市東南部，區(qū)間變化不大，重心分布區(qū)域較為穩(wěn)定，表明林業(yè)企業(yè)呈集中分布。（2）在長短軸方面，2000—2020 年橢圓長軸由2.008 增長至2.169，短軸由1.157 縮短至1.096，林業(yè)企業(yè)在東北—西南方向呈擴張態(tài)勢，在西北—東南方向上呈收縮態(tài)勢。（3）在旋轉(zhuǎn)角度方面，方位角度數(shù)由26.299縮短至20.218，逐漸減小，林業(yè)企業(yè)數(shù)量在西北—東南方向增長較快，但空間分布仍沿東北—西南方向。在西安市、咸陽市等關(guān)中地區(qū)，林業(yè)企業(yè)密度的增加更為迅速，該區(qū)域經(jīng)濟發(fā)達，有利于林業(yè)企業(yè)的生存和發(fā)展。2000—2020 年林業(yè)企業(yè)標準差橢圓重心雖有小幅遷移，但一直位于該地區(qū)。

圖2 2000—2020年陜西省林業(yè)企業(yè)標準差橢圓分析結(jié)果Fig.2 Standard deviation ellipse analysis results of forestry enterprises in Shaanxi Province during 2000—2020

2.1.3 陜西省林業(yè)企業(yè)時空分布密度變化通過核密度分析法對陜西省林業(yè)企業(yè)空間分布進行分析，將結(jié)果按照自然斷點法分為5 類（圖3）。在各個時間節(jié)點下陜西省林業(yè)企業(yè)的分布密度均呈現(xiàn)出顯著的空間集聚性。2000—2020年，林業(yè)企業(yè)數(shù)量不斷增加，核密度值不斷提升，最大值由1974 提升至75381，增長了約37 倍。分布空間由2000 年的相對均衡過渡到2020年的絕對集中，呈現(xiàn)出以西安市為中心的塊狀集聚，且范圍逐漸增大。

圖3 2000—2020年陜西省林業(yè)企業(yè)空間分布密度Fig.3 Spatial distribution density of forestry enterprises in Shaanxi Province during 2000—2020

2000 年陜西省林業(yè)企業(yè)在西安市、漢中市、寶雞市、延安市等地呈現(xiàn)出多個分散的塊狀集聚區(qū)，林業(yè)企業(yè)集聚區(qū)與黃土高原水土保持區(qū)、關(guān)中平原生態(tài)協(xié)同發(fā)展區(qū)和秦巴水源涵養(yǎng)區(qū)相對應(yīng)。林業(yè)企業(yè)在空間分布上形成塊狀集聚、多中心發(fā)展的格局，發(fā)展相對均衡（圖3a）。2005年各市熱點區(qū)域的數(shù)量明顯下降，集聚中心以各市市區(qū)為中心，更加凸顯。熱點區(qū)域主要集中在關(guān)中平原生態(tài)協(xié)同發(fā)展區(qū)，這一階段，資源型林業(yè)企業(yè)數(shù)量減少，林業(yè)企業(yè)向市區(qū)集聚程度明顯上升（圖3b）。2010 年陜西省林業(yè)企業(yè)小型集聚中心基本消失，大型集聚中心范圍略有擴大，以西安市為最高等級集聚中心，向西北、東北方向擴散（圖3c）。2015年陜西省林業(yè)企業(yè)在安康市、商洛市、榆林市出現(xiàn)新的小型集聚中心（圖3d）。2020 年西安市和咸陽市的臨近交界區(qū)域成為陜西省林業(yè)企業(yè)集聚中心，西咸一體化進程的推進，吸引更多的新注冊企業(yè)在此發(fā)展，呈現(xiàn)沿渭河的集中連片帶狀分布。榆林市在毛烏素沙地生態(tài)修復(fù)區(qū)新增一個集聚區(qū)?？傮w表明，林業(yè)企業(yè)發(fā)展主要依賴于資本和市場，其次受林業(yè)資源影響較大（圖3e）。

2.2 林業(yè)企業(yè)發(fā)展狀況

2.2.1 陜西省林業(yè)企業(yè)經(jīng)營范圍變化從經(jīng)營范圍詞云圖來看，各個時間點上林業(yè)類企業(yè)的業(yè)務(wù)范圍既有共同點，又具有不同階段的鮮明特征（圖4）。如數(shù)量排名前10的關(guān)鍵詞基本一致，以林業(yè)類相關(guān)的基本產(chǎn)品服務(wù)關(guān)鍵詞為主，具有較大的共性。主要包括經(jīng)營屬性類關(guān)鍵詞：“林場”“林業(yè)局”“有限公司”等；運營方式類關(guān)鍵詞：“銷售”“服務(wù)”“管理”“批發(fā)”“工程”等；林業(yè)相關(guān)類關(guān)鍵詞：“林業(yè)”“農(nóng)業(yè)”“林木”等。其中“林業(yè)”“銷售”“服務(wù)”3 個關(guān)鍵詞的數(shù)量始終是前三名，林業(yè)類企業(yè)的業(yè)務(wù)范圍鮮明，即與林業(yè)相關(guān)的產(chǎn)品銷售或者衍生服務(wù)，如林業(yè)產(chǎn)品銷售、林業(yè)機械服務(wù)、栽培服務(wù)、病蟲害防治服務(wù)等。而后的關(guān)鍵詞則呈現(xiàn)出各個時間節(jié)點上的不同特征，從“苗圃”“種植”到“生物防治”“設(shè)備”再到“科技”“技術(shù)”，體現(xiàn)出林業(yè)企業(yè)的經(jīng)營范圍從初加工、再加工到深加工的深入發(fā)展歷程，林業(yè)企業(yè)在保持自身產(chǎn)業(yè)特點的同時，科技含量不斷增加，并擴增相關(guān)服務(wù)項目。

圖4 2000—2020年陜西省林業(yè)企業(yè)經(jīng)營范圍詞云圖Fig.4 Word cloud charts of forestry enterprise business scope in Shaanxi Province during 2000—2020

分時間段來看（圖5），2000 年林業(yè)類企業(yè)的國有化特征明顯，關(guān)鍵詞“林業(yè)局”“林場”出現(xiàn)次數(shù)高，均高于“有限公司”的數(shù)量。表明林業(yè)類企業(yè)的市場化進程仍在進行之中，國有類傳統(tǒng)林業(yè)企業(yè)占據(jù)主導(dǎo)地位，而林業(yè)類民營企業(yè)較少。2005年，“有限公司”關(guān)鍵詞數(shù)量開始超過“林業(yè)局”，隨著市場經(jīng)濟的興起，國家出臺《關(guān)于加快林業(yè)發(fā)展的決定》，鼓勵非公有林業(yè)的發(fā)展，資本開始進入林業(yè)產(chǎn)業(yè)，民營企業(yè)增多，并出現(xiàn)了“投資”關(guān)鍵詞。2010年開始，林業(yè)產(chǎn)品銷售和衍生服務(wù)蓬勃興起，“生物防治”“生態(tài)”“設(shè)備”等林業(yè)的衍生產(chǎn)品增多，林業(yè)企業(yè)的經(jīng)營范圍逐漸壯大，并集中在和政策、產(chǎn)業(yè)升級等相關(guān)的領(lǐng)域。2015年“技術(shù)咨詢”“項目”“生態(tài)”等關(guān)鍵詞快速增長，但“農(nóng)業(yè)”“種植”關(guān)鍵詞仍保持在前幾名，表明林業(yè)企業(yè)的發(fā)展是基數(shù)的擴大和增量的特色化。2020 年關(guān)鍵詞數(shù)量增長幅度最大，表明2015—2020是近20 a陜西省林業(yè)產(chǎn)業(yè)發(fā)展最快的時間段，期間林業(yè)產(chǎn)業(yè)規(guī)模迅速壯大，產(chǎn)品貿(mào)易十分旺盛，“銷售”“服務(wù)”關(guān)鍵詞激增，處于關(guān)鍵詞前兩位。同時，林業(yè)產(chǎn)業(yè)結(jié)構(gòu)進一步優(yōu)化，產(chǎn)業(yè)類型豐富化，包括“花卉”“技術(shù)”“生物防治”“園林綠化”“項目”等。

圖5 2000—2020年陜西省林業(yè)企業(yè)經(jīng)營范圍關(guān)鍵詞Fig.5 Keywords for business scope of forestry enterprises in Shaanxi Province during 2000—2020

2.2.2 陜西省林業(yè)企業(yè)分布數(shù)量影響因素分析通過對陜西省林業(yè)企業(yè)的空間集聚和分布方向分析，表明其具有空間異質(zhì)性，因此考慮采用地理加權(quán)回歸模型對其影響因素進行實證分析。

根據(jù)已有研究成果［31-32］，同時考慮到數(shù)據(jù)的可獲取性，本研究選取交通便利性、勞動力數(shù)量、企業(yè)發(fā)展狀況、經(jīng)濟發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、生產(chǎn)條件、市場規(guī)模等陜西省林業(yè)企業(yè)空間分布的影響因素進行分析（表2）。

表2 林業(yè)企業(yè)區(qū)位選擇的影響因素及定義Tab.2 Influencing factors and definition of location choice of forestry enterprises

（1）OLS模型

地理加權(quán)回歸模型使用的前提條件是因變量（各縣區(qū)林業(yè)企業(yè)數(shù)量）具有空間自相關(guān)性，如果存在空間自相關(guān)性，則OLS 模型估計有偏，而更適用GWR模型［33］。

首先使用OLS 模型測定各影響因素對于陜西省林業(yè)企業(yè)分布的影響，由于各項指標單位不同，數(shù)值差距大，為盡量消除異方差的影響，對各項變量數(shù)值進行歸一化處理，同時對所選指標進行合理性檢驗（表3）?；貧w結(jié)果顯示，決定系數(shù)（R2）=0.878，調(diào) 整R2=0.869，赤 池 信 息 準 則（AICc）=1187.205，全局擬合效果良好，常住人口數(shù)量、地區(qū)生產(chǎn)總值、第一產(chǎn)業(yè)增加值占比、林地面積、社會消費品零售總額5 項指標對林業(yè)企業(yè)分布起正向作用，而公路密度、企業(yè)注冊資本平均值、園地面積3項指標對林業(yè)企業(yè)分布起負向作用，8 項指標的共線關(guān)系較低，所選指標合理。

表3 普通最小二乘法（OLS）模型估計結(jié)果Tab.3 Estimation results of OLS model

（2）GWR模型

OLS 模型的Koenker (BP)統(tǒng)計量為0.036×10-3，表明模型中未考慮到空間異質(zhì)性等因素，可以利用GWR模型處理空間異質(zhì)性和異方差性問題，深入探究各個因素對于林業(yè)企業(yè)空間分布的影響。

GWR 模型結(jié)果顯示，R2=0.900、調(diào)整R2=0.886、AICc=1178.622，根據(jù)表中對回歸系數(shù)中位數(shù)、最小值、最大值和平均值的統(tǒng)計，系數(shù)值有明顯的浮動變化，表明不同縣區(qū)林業(yè)企業(yè)分布影響因素存在顯著差異，而調(diào)整后的R2更大，AICc 差值大于3，說明GWR模型更優(yōu)，能更好地解釋各變量對林業(yè)企業(yè)分布數(shù)量的影響。

從GWR模型回歸系數(shù)來看（表4），社會消費品零售總額、地區(qū)生產(chǎn)總值、常住人口數(shù)量、第一產(chǎn)業(yè)增加值占比、林地面積5 項因素的回歸系數(shù)在全部縣區(qū)均為正，園地面積、公路密度兩項因素的回歸系數(shù)在全部縣區(qū)均為負，企業(yè)注冊資本平均值的回歸系數(shù)在不同縣區(qū)呈現(xiàn)正、負的差異性。具體表現(xiàn)在：

表4 地理加權(quán)回歸（GWR）模型系數(shù)Tab.4 GWR model coefficients

（1）社會消費品零售總額因素對于林業(yè)企業(yè)分布數(shù)量的影響最為顯著。林業(yè)產(chǎn)業(yè)是涵蓋一、二、三產(chǎn)業(yè)的復(fù)合產(chǎn)業(yè)群體，社會消費品零售總額越高，市場規(guī)模越大、商業(yè)發(fā)達程度越高，越有利于林產(chǎn)品的線上、線下銷售。對于林業(yè)產(chǎn)業(yè)鏈中的第三產(chǎn)業(yè)吸引力更大，使得以林產(chǎn)品銷售、提供林業(yè)服務(wù)的小型林業(yè)企業(yè)匯聚，導(dǎo)致林業(yè)企業(yè)分布數(shù)量更多，因此回歸系數(shù)均為正。由于陜南林業(yè)產(chǎn)業(yè)以第一產(chǎn)業(yè)為主，市場規(guī)模相對較小，而陜北榆林市經(jīng)濟相對發(fā)達，市場也更廣闊，因此回歸系數(shù)呈現(xiàn)由南向北遞增的趨勢（圖6a）。

圖6 陜西省林業(yè)企業(yè)分布數(shù)量與各影響因素的回歸系數(shù)Fig.6 Regression coefficient between distribution quantity of forestry enterprises in Shaanxi Province and various influencing factors

（2）地區(qū)生產(chǎn)總值因素對于林業(yè)企業(yè)分布數(shù)量的影響較為顯著。GDP 是地區(qū)經(jīng)濟實力的代表性因素，一般來說，GDP越高，越能夠保障地區(qū)產(chǎn)業(yè)發(fā)展的可持續(xù)性，并提供良好的企業(yè)發(fā)展外部環(huán)境，有利于吸引企業(yè)入駐，因此回歸系數(shù)均為正，GDP與林業(yè)企業(yè)分布數(shù)量呈現(xiàn)正相關(guān)關(guān)系。陜南地區(qū)安康市和漢中市屬于GDP單核城市，市區(qū)GDP遠大于其他縣區(qū)，林業(yè)企業(yè)也在市區(qū)形成集聚，因此陜南地區(qū)的回歸系數(shù)最高（圖6b）。

（3）常住人口數(shù)量因素的回歸系數(shù)較為顯著，且均為正相關(guān)。陜西省林業(yè)產(chǎn)業(yè)體系發(fā)展態(tài)勢良好，以第一、第二產(chǎn)業(yè)為主，經(jīng)營粗放，大而不強，第三產(chǎn)業(yè)占比偏小，以勞動密集型產(chǎn)業(yè)為主，對于勞動力的需求大。高值區(qū)域主要位于漢中市和寶雞市2個林業(yè)大市，但整體上系數(shù)偏差不大，影響強度較為接近（圖6c）。

（4）第一產(chǎn)業(yè)增加值占比因素的回歸系數(shù)較為顯著，且均為正相關(guān)。陜西省林業(yè)產(chǎn)業(yè)結(jié)構(gòu)體系中經(jīng)濟林、林果業(yè)等第一產(chǎn)業(yè)占有較大比重，第一產(chǎn)業(yè)增加值占比高的地區(qū)對于林業(yè)產(chǎn)業(yè)的依賴更高，如漢中市是陜西省森林產(chǎn)業(yè)重心區(qū)，林業(yè)企業(yè)數(shù)量也就更多。回歸系數(shù)從陜西省西南向東北過渡，與陜西省林地資源分布也較為一致（圖6d）。

（5）林地面積因素對于林業(yè)企業(yè)分布數(shù)量的影響較為顯著，且均為正。林地面積是刻畫林業(yè)資源豐富程度的重要指標，林業(yè)產(chǎn)業(yè)的特殊性，使得林業(yè)產(chǎn)業(yè)中的第一、第二產(chǎn)業(yè)都更加依賴林木資源，因此林地面積較多的地區(qū)，林業(yè)企業(yè)分布數(shù)量也更多。從回歸系數(shù)變化來看，在漢中市、寶雞市等林地資源豐富的地區(qū)影響程度強，而在林地資源相對匱乏的榆林市影響程度相對較弱（圖6e）。

（6）園地面積因素的回歸系數(shù)較為顯著，且均為負相關(guān)。陜西省園地約為林地面積的10%，主要分布在渭南市、榆林市、延安市、咸陽市，園地以果園為主，而林果產(chǎn)業(yè)具有占地面積大，相對分散的特點，以小型、個體企業(yè)為主，因此回歸系數(shù)呈負值，但絕對值較小。陜南園地面積較小，回歸系數(shù)絕對值較高（圖6f）。

（7）公路密度因素的回歸系數(shù)均為負，與預(yù)期期望不一致。推測為公路密度越高，城市化水平越高，林業(yè)資源反而并不豐富，不利于林業(yè)第一產(chǎn)業(yè)的布局和發(fā)展。此外，陜西省地形地貌特殊，除中部關(guān)中平原道路密度較大外，陜北、陜南均為高原或山區(qū)，兩地所有縣區(qū)的道路密度均較低且接近。造成陜北、陜南林業(yè)企業(yè)分布數(shù)量多的地區(qū)，道路密度比較低，公路密度與林業(yè)企業(yè)分布數(shù)量呈現(xiàn)負相關(guān)關(guān)系（圖6g）。

（8）企業(yè)注冊資本平均值因素在陜西省大部分地區(qū)呈負向作用。林業(yè)企業(yè)多數(shù)為資源依賴型企業(yè)，而地方的林地面積和森林資源是有限的，企業(yè)注冊資本平均值越高，企業(yè)的平均規(guī)模越大，市場化競爭更加激烈，不利于新的企業(yè)入駐?；貧w系數(shù)從南向北呈現(xiàn)遞減的態(tài)勢，表明企業(yè)注冊資本平均值對林業(yè)企業(yè)分布數(shù)量的相對壓力逐漸降低?；貧w系數(shù)為正的縣區(qū)位于榆林市，該市地處毛烏素沙地腹地，是三北工程的重要節(jié)點，林業(yè)產(chǎn)業(yè)與其他地區(qū)較為不同，除經(jīng)濟林果產(chǎn)業(yè)外，森林康養(yǎng)、林下經(jīng)濟、生態(tài)旅游、林業(yè)工程也較為發(fā)達，政府引導(dǎo)和支持力度比較大，林業(yè)企業(yè)數(shù)量增加的同時，企業(yè)注冊資本也維持在高值（圖6h）。

3 討論

在以往的研究中，企業(yè)的空間格局和驅(qū)動因素分析已經(jīng)較為成熟［17-23］，但關(guān)于林業(yè)等第一產(chǎn)業(yè)的企業(yè)格局演化，還鮮有探索，研究基礎(chǔ)較為薄弱［25-27］。本文結(jié)合其他產(chǎn)業(yè)研究中的方法，采用GIS 空間分析方法，定量研究了陜西省林業(yè)企業(yè)的時空演變特征。同時，考慮到林業(yè)產(chǎn)業(yè)的特點和不同，從交通便利性、勞動力數(shù)量、企業(yè)發(fā)展狀況、經(jīng)濟發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、生產(chǎn)條件、市場規(guī)模等方面，利用OLS和GWR 模型探究了陜西省林業(yè)企業(yè)空間分布影響因素的作用和強度。陜西省林業(yè)企業(yè)的空間分布格局表明，林業(yè)企業(yè)作為以第一、第二產(chǎn)業(yè)為主的產(chǎn)業(yè)，其遷移速度較為緩慢，遷移過程是連續(xù)的。從林業(yè)經(jīng)營范圍關(guān)鍵詞來看，陜西省林業(yè)產(chǎn)業(yè)的結(jié)構(gòu)調(diào)整比較明顯，林業(yè)服務(wù)類企業(yè)等林業(yè)第三產(chǎn)業(yè)中的企業(yè)數(shù)量明顯提升。而陜西省林業(yè)企業(yè)分布數(shù)量因素的測度結(jié)果表明，社會消費品零售總額、常住人口數(shù)量、地區(qū)生產(chǎn)總值、第一產(chǎn)業(yè)增加值占比、林地面積5項因素具有較為顯著的正相關(guān)作用，與陜西省林業(yè)產(chǎn)業(yè)第一、第二產(chǎn)業(yè)占比較大，資源依賴型和勞動密集型的特點顯著相關(guān)。企業(yè)注冊資本平均值、公路密度、園地面積等因素的回歸系數(shù)主要為負，和林業(yè)企業(yè)的分布數(shù)量呈現(xiàn)負相關(guān)作用。其中，陜西省特殊的地貌特征使得除關(guān)中平原外的縣區(qū)道路密度都比較低，但漢中市、安康市等林業(yè)資源豐富的地區(qū)又存在較多的林業(yè)企業(yè)，從而產(chǎn)生了和以往企業(yè)分布研究不一樣的影響效果。

因此，建議政府部門從以下3 個方面優(yōu)化林業(yè)資源的配置。（1）做精第一產(chǎn)業(yè)：陜西省林果產(chǎn)業(yè)發(fā)達，應(yīng)利用好林業(yè)資源和人力資源，抓住以蘋果、核桃、獼猴桃等為主的特色經(jīng)濟林產(chǎn)業(yè)，促進產(chǎn)品的提煉與深加工。（2）做強第二產(chǎn)業(yè)：發(fā)揮市場配置資源的決定性作用和產(chǎn)業(yè)集聚區(qū)的溢出輻射效應(yīng)，促進林業(yè)產(chǎn)業(yè)向優(yōu)勢區(qū)集聚，強化產(chǎn)業(yè)聚集發(fā)展，形成第一、第二產(chǎn)業(yè)的產(chǎn)業(yè)鏈融合，平衡好資源和市場優(yōu)勢。（3）做大第三產(chǎn)業(yè)：充分挖掘陜西省林業(yè)產(chǎn)業(yè)中的旅游、康養(yǎng)、文化價值，培育新業(yè)態(tài)新產(chǎn)業(yè)，利用好資本優(yōu)勢，打造一批特色突出的服務(wù)品牌。

4 結(jié)論

陜西省林業(yè)資源豐富多樣，其林業(yè)企業(yè)的空間演化特征及分布數(shù)量的影響因素探究具備一定代表性。主要結(jié)論如下：

（1）從時間變化上看，2000—2020年，陜西省林業(yè)企業(yè)數(shù)量增長的速度逐漸加快，林業(yè)企業(yè)發(fā)展勢頭較好。陜西省林業(yè)企業(yè)最近鄰指數(shù)呈減小趨勢，空間分布類型均為聚集型，且集聚程度不斷增強。

（2）從空間分布上看，陜西省林業(yè)企業(yè)分布標準差橢圓大致呈西南—東北方向，空間分布重心均位于咸陽市東部，較為穩(wěn)定。核密度熱點區(qū)始終以西安市區(qū)為中心且無較大變化，但次級熱點區(qū)范圍則逐步擴大，沿連霍（連云港—霍爾果斯）高速向西擴展至楊凌示范區(qū)。

（3）從林業(yè)企業(yè)的經(jīng)營范圍來看，林業(yè)類企業(yè)的業(yè)務(wù)范圍鮮明但經(jīng)營范圍逐漸擴大，主要以林業(yè)相關(guān)的產(chǎn)品銷售或者衍生服務(wù)為主，產(chǎn)業(yè)結(jié)構(gòu)進一步優(yōu)化，產(chǎn)業(yè)類型逐漸豐富。表明陜西省林業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整成效顯著，第一、第二產(chǎn)業(yè)加快轉(zhuǎn)型升級，第三產(chǎn)業(yè)貢獻率明顯提高。

（4）從林業(yè)企業(yè)分布數(shù)量的影響因素來看，社會消費品零售總額、地區(qū)生產(chǎn)總值、常住人口數(shù)量、第一產(chǎn)業(yè)增加值占比、林地面積5 項因素與林業(yè)企業(yè)分布數(shù)量呈正相關(guān)性，其中社會消費品零售總額的影響強度最大；園地面積、公路密度、企業(yè)注冊資本平均值與林業(yè)企業(yè)分布數(shù)量呈現(xiàn)負相關(guān)性。表明林業(yè)企業(yè)數(shù)量分布需要在同類型企業(yè)的競爭壓力與地區(qū)市場容量、經(jīng)濟發(fā)展水平、資源優(yōu)勢之間找到平衡。