陳科屹，何友均①，謝和生，李艷娟，趙雨昕 (.中國林業(yè)科學研究院林業(yè)科技信息研究所，北京 0009；.大興安嶺林業(yè)集團公司，黑龍江 大興安嶺 65000)

自1964年8月10日黨中央、國務院正式決定開發(fā)建設大興安嶺以來，大興安嶺國有林區(qū)已經歷58年的持續(xù)開發(fā)歷程。在開發(fā)初期，經濟建設是林區(qū)的核心目標，經濟結構以木材采伐和初級加工占據絕對主導地位，為國家經濟發(fā)展和現代化建設提供了木材支持。由于長期的采育失調，導致可采森林資源銳減，林區(qū)逐漸暴露出資源危機、經濟危困、生態(tài)退化等一系列問題，林區(qū)也由因林而生、因林而興逐步走到了因林而衰。2000年之后，國家先后實施了天然林資源保護工程一、二期等林業(yè)重點工程和全面停止天然林商業(yè)性采伐等政策措施。至此，包括黑龍江省大興安嶺林區(qū)在內的國有林區(qū)開始大力探索林區(qū)轉型發(fā)展路徑，謀求盡快改變獨木支撐的單一經濟發(fā)展格局及推動實現多領域全方位協調發(fā)展成為國有林區(qū)新的出路。

協調發(fā)展是可持續(xù)發(fā)展思想的核心內容，也是新發(fā)展理念的重要組成部分。所謂協調，即正確處理系統組織內外各種關系，系統內部各子系統能夠通過信息、能量和物質交換相互配合，使整個系統綜合功能超過各子系統之和，形成更加有序的系統[1]。國外學者較早開始關注協調發(fā)展問題，涉及經濟學[2]、地理學[3]、生態(tài)學[4]、社會學[5]、醫(yī)學[6]和材料學[7]等多個學科領域。隨著國有林區(qū)轉型發(fā)展進程的持續(xù)推進，近年來已有部分學者關注國有林區(qū)協調發(fā)展問題，主要圍繞生態(tài)建設與經濟發(fā)展[8-10]，以及生態(tài)保護與民生保障[11-12]等兩兩之間的協調關系開展研究，涉及協調發(fā)展的關聯因子、干擾因素和優(yōu)化路徑等方面。王玉芳[13]認為市場、政策失效影響著經濟與生態(tài)之間的協調，技術是解決經濟與生態(tài)協調的主要手段；李朝洪等[14]認為國有林區(qū)是生態(tài)和經濟相互作用、相互依賴、相互制約的復雜系統。也有學者從生態(tài)、經濟和社會3個方面協調發(fā)展狀況對國有林區(qū)發(fā)展情況進行研究。馬文學[15]研究發(fā)現伊春國有林區(qū)生態(tài)、經濟和社會之間的協調關系具有明顯階段性，在不同發(fā)展時期其協調度和協調機制不同；朱曉柯等[16]認為三元系統耦合協調水平的變化與二元耦合協調水平密切相關，在黑龍江省國有林區(qū)，生態(tài)-產業(yè)二元耦合系統和生態(tài)-民生二元耦合系統是制約三元系統耦合協調水平的主要因素。多數學者研究結果表明，我國國有林區(qū)協調發(fā)展狀況總體上呈現持續(xù)向好態(tài)勢，已有研究成果為進一步研究國有林區(qū)耦合協調發(fā)展提供了理論和實踐參考。但是，目前研究缺乏對作為國有林區(qū)森林資源保護發(fā)展主體的森工企業(yè)自身建設的關注，以及林區(qū)基礎設施建設等支撐保障狀況的考慮，而將兩者作為獨立系統層與生態(tài)經濟社會相結合，探究其作用關系和運行規(guī)律的研究還鮮見報道。如果國有林區(qū)缺乏高效的經營管理主體和強大的支撐保障能力，最終也會影響區(qū)域生態(tài)、經濟和社會協調發(fā)展狀況。尤其是2020年之后，黑龍江省大興安嶺國有林區(qū)已順利通過國家改革驗收，組建了新的大興安嶺林業(yè)集團公司，結束了黑龍江省大興安嶺地區(qū)56年的政企合一管理體制。站在新的歷史發(fā)展階段和國家戰(zhàn)略定位，黑龍江省大興安嶺重點國有林區(qū)的發(fā)展面臨新的形勢和要求。

鑒于此，將企業(yè)管理和支撐保障相關內容納入研究范疇，構建基于資源保育、產業(yè)發(fā)展、企業(yè)管理、民生福祉和支撐保障的黑龍江省大興安嶺重點國有林區(qū)復合系統，采用耦合協調模型對2000—2020年國有林區(qū)耦合協調發(fā)展狀況進行定量評價，并采用灰色馬爾科夫模型對2021—2022年復合系統耦合協調發(fā)展趨勢進行預測，以期為推動黑龍江省大興安嶺重點國有林區(qū)持續(xù)深化改革，促進實現全方位和高質量協調發(fā)展提供科學依據和理論支撐。

1 研究區(qū)概況

黑龍江省大興安嶺重點國有林區(qū)主要指大興安嶺林業(yè)集團公司所轄的國有林區(qū)。截至2020年末，大興安嶺林業(yè)集團公司下轄10個林業(yè)局、100個管護區(qū)(林場)、8個國家級自然保護區(qū)，森林覆蓋率為86.14%，林地保有量為798萬hm2，活立木蓄積為6.02億m3；擁有在崗職工35 770人，年末實有離退休人員61 255人，林業(yè)產業(yè)年總產值為534 180 萬元。

2 材料與方法

2.1 構建評價體系

基于協調發(fā)展理論[17]，立足黑龍江省大興安嶺重點國有林區(qū)的自然資源稟賦及其在國家戰(zhàn)略中的角色定位，結合2000—2020年國有林區(qū)轉型發(fā)展成效和大興安嶺林業(yè)集團公司作為新型經營管理主體的主責主業(yè)，參考學者們提出的林業(yè)現代化指標體系[18-19]、高質量發(fā)展指標體系[20]和國有林區(qū)轉型發(fā)展相關研究成果[14]，遵循指標代表性、數據獲得性、結果實用性等原則，構建黑龍江省大興安嶺重點國有林區(qū)協調發(fā)展評價指標體系。

國有林區(qū)的建設發(fā)展主要圍繞林草資源的保護、修復和合理利用問題組織開展相關活動，其建設發(fā)展水平受到區(qū)域資源保育、產業(yè)發(fā)展水平、森工企業(yè)管理、民生福祉和生產生活支撐保障5個方面因素的綜合影響，5大方面共同構成一個復合系統，具有耦合性特征(圖1)。

圖1 黑龍江省大興安嶺重點國有林區(qū)復合系統耦合關系

其中，資源保有量是生態(tài)系統發(fā)揮生態(tài)功能的載體，資源保育狀況是耦合協調發(fā)展的基礎；產業(yè)發(fā)展是持續(xù)推動國有林區(qū)耦合協調發(fā)展的動力源泉，產業(yè)發(fā)展模式和質量直接影響生態(tài)保護成效和民生福祉水平等；企業(yè)管理情況影響著資源保護發(fā)展全過程各個環(huán)節(jié)的效率，是調和各子系統層實現有機互動的關鍵一環(huán)；民生福祉情況關系林區(qū)社會和諧安定，是體現耦合協調發(fā)展成效最直接的反映；支撐保障是保證國有林區(qū)各項事業(yè)穩(wěn)步發(fā)展的后盾。圍繞上述5個子系統層和14個準則層初步選取評價指標因子，經篩選剔除最終得到36個指標因子形成指標層，構建反映黑龍江省大興安嶺重點國有林區(qū)耦合協調發(fā)展狀況的指標體系(表1)。

2.2 數據來源

基于構建的黑龍江省大興安嶺重點國有林區(qū)耦合協調發(fā)展評價指標體系，2000—2020年各指標因子數據通過《中國林業(yè)統計年鑒》《黑龍江統計年鑒》、大興安嶺林業(yè)集團公司及相關單位內部資料、全國森林資源清查資料、中國林業(yè)信息網等來源收集和整理得到。其中，森林覆蓋率、森林蓄積量等源于全國森林資源連續(xù)清查數據的指標，由于每5年清查1次，因此發(fā)布年之外其他年份的數據采用線性插值法進行補充。

2.3 熵權計算

首先，將2001—2020年所有價值量轉化為2000年價格計算的價值量。其中，林業(yè)產業(yè)總產值、林業(yè)產業(yè)增加值采用國內生產總值指數(上年=100)進行平減；在崗職工年均工資、離退休人員年均生活費采用消費者價格指數(上年=100)進行平減；招商引資、林業(yè)固定資產投資、企業(yè)自籌資金采用固定資產投資價格指數(上年=100)進行平減。然后，采用Min-Max標準化方法對所有數據進行標準化處理。最后，計算評價指標信息熵和熵權。上述步驟計算公式[21]為

(1)

(2)

(3)

式(1)～(3)中，x0ij為第i個指標在第j年的原始值或預處理值；xij為第i個指標在第j年的標準化值；ei為第i個指標的熵；wi為第i個指標的熵權；i∈[1,m]，m為指標個數；j∈[1,n]，n為年份。

表1 黑龍江省大興安嶺重點國有林區(qū)耦合協調發(fā)展評價指標體系及權重

自然公園面積=森林公園面積+濕地公園面積；接續(xù)替代產業(yè)產值為林產品種植與采集、花卉及其他觀賞植物種植、陸生野生動物繁育與利用、生物質能源、非木質林產品加工制造、綠色礦業(yè)及林業(yè)旅游與休閑服務6個產業(yè)產值之和；企業(yè)就業(yè)指數=在崗職工/在冊職工人數；從業(yè)人員流動率=(調入人員+調出人員)/從業(yè)人員。

2.4 耦合協調模型構建

采用耦合協調模型度量復合系統內部的良性耦合程度，具體表達式[22]為

(4)

(5)

(6)

式(4)～(6)中，C為耦合度，0≤C≤1；D為耦合協調度，0≤D≤1；T為協調發(fā)展指數；a、b、c、d和e分別為根據式(2)～(3)計算的各子系統層權重(表1)，即a=0.242，b=0.161，c=0.271，d=0.105，e=0.221；Uq為各子系統層(q=1、2、3、4、5)或復合系統(q=0)的發(fā)展指數。

參照王玉芳等[23]和李朝洪等[8]對耦合協調水平的劃分，將其劃分為3個發(fā)展階段和10個發(fā)展類型(表2)。

表2 耦合協調度等級劃分

2.5 趨勢預測

采用灰色馬爾科夫模型對復合系統耦合協調度進行預測，主要包括如下5個步驟[24-25]：

(7)

(8)

式(7)～(8)中，D(0)={d(0)(1),d(0)(2),…,d(0)(n)}為復合系統耦合協調度序列；D(1)={d(1)(1),d(1)(2),…,d(1)(n)}為一次累加序列；t為發(fā)展灰數；u為內生控制灰數；e為自然對數；k=0,1,…,n。

(2)依據灰色預測值相對誤差大小劃分狀態(tài)區(qū)間Zr，計算公式為

(9)

Zr=[fr,gr]。

(10)

式(9)～(10)中，ER為相對誤差；fr和gr分別為狀態(tài)區(qū)間上、下限，r=1,2,…,s，s為區(qū)間劃分數量。

(3)計算k步狀態(tài)轉移概率矩陣P(k)，計算公式為

(11)

(12)

式(11)～(12)中，Nr為狀態(tài)Zr出現頻次；Nrh(k)為狀態(tài)Zr經過k步轉移成Zh出現頻次，h≤s。

(4)對灰色預測值進行修正，得到灰色馬爾科夫預測值(y),計算公式為

(13)

(5)運用殘差檢驗、后驗差檢驗以及關聯分析法檢驗對預測結果進行精度檢驗，計算公式為

(14)

B=S2/S1，

(15)

(16)

3 結果與分析

3.1 系統發(fā)展指數分析

各系統層發(fā)展指數見表3。

表3 各子系統層發(fā)展指數和復合系統綜合發(fā)展指數

由表3可知，總體上各子系統層發(fā)展指數和復合系統綜合發(fā)展指數在一定程度上均有所提升，但具體過程則有所區(qū)別。其中，民生福祉和資源保育2個子系統層發(fā)展指數呈近似線性增長態(tài)勢，年均增幅分別達25%和8%。受苗木產量和育苗面積大幅縮減的影響，資源保育系統層增長趨勢在2004和2007年出現小幅度波動。產業(yè)發(fā)展系統層發(fā)展指數也呈現較好增長趨勢，但2020年受新冠肺炎疫情等因素嚴重影響，發(fā)展指數出現斷崖式下降。企業(yè)管理系統層發(fā)展指數呈近似拋物線趨勢，在2000—2013年期間呈現波動上升趨勢，2013年達到峰值。從2014年開始，黑龍江省大興安嶺重點國有林區(qū)正式全面停止天然林商業(yè)性采伐，而接續(xù)替代產業(yè)尚未發(fā)展壯大，由此引發(fā)企業(yè)收入銳減、招引資金驟降等情況，導致企業(yè)管理系統層發(fā)展指數出現連續(xù)下降。支撐保障系統層發(fā)展指數變化幅度相對最小，與其他系統層低值發(fā)生在早期的現象有所不同，支撐保障系統層低值集中出現在近期，尤其是2014年之后僅能維持在0.288～0.347之間?？梢姡邶埥〈笈d安嶺國有林區(qū)支撐保障薄弱的問題在近年來尤為突出，譬如在冊林業(yè)公路里程數在近21年增幅不足5%，且因管理養(yǎng)護不足，能正常使用的實際里程數不增反減。復合系統綜合發(fā)展指數集中體現了各子系統層的綜合影響效應，呈現增長→下降→再恢復的過程，但由于受2020年新冠肺炎疫情影響，近年來持續(xù)的恢復性增長被中斷。綜合來看，2000—2020年黑龍江省大興安嶺重點國有林區(qū)各子系統發(fā)展方向和發(fā)展勢態(tài)處于動態(tài)變化中，子系統層之間的發(fā)展變化匹配程度和相互作用關系直接影響著復合系統綜合發(fā)展情況。各子系統層均出現過發(fā)展相對滯后的情況，其中，改善效果最明顯的是民生福祉系統層，而相對滯后頻率最高的是企業(yè)管理系統層。從近10年情況來看，相對滯后問題最突出的是支撐保障系統層。

3.2 復合系統耦合協調發(fā)展狀況分析

由表4可知，在耦合度方面，黑龍江省大興安嶺重點國有林區(qū)復合系統耦合度指數為0.656～0.992，始終維持在中高水平，表明5個子系統層間相互依賴、相互作用程度較高，具有較好的耦合性。但是，耦合性高低不能反映子系統層間作用關系的協調取向，需進一步分析子系統層間的耦合協調狀況。在耦合協調度方面，2000—2020年黑龍江省大興安嶺重點國有林區(qū)子系統層間的耦合協調發(fā)展情況大致可劃分為3個時期。第1個時期為2000—2001年，子系統層間處于失調發(fā)展階段，但是隨著時間的推移趨于逐漸脫離輕度失調衰退發(fā)展類型。第2個時期為2002—2009年，子系統層間處于過渡發(fā)展階段，此時已由前一時期的輕度失調衰退發(fā)展類型改善為勉強協調發(fā)展類型。第3個時期為2010年之后，子系統層間的耦合協調發(fā)展階段進入到協調發(fā)展階段，耦合協調度在2013年達到峰值。在2014和2020年，子系統層間的耦合協調度出現2次較大幅度下降，表明此時子系統層間相互掣肘的矛盾再次凸顯。然而即便如此，子系統層間的耦合協調情況也要明顯優(yōu)于2010年及之前的情況。雖然在2010—2020年整個時期，子系統層間的耦合協調狀況均處于協調發(fā)展階段，但耦合協調發(fā)展的水平還不高，除2013年外，僅處于初、中級的低水平發(fā)展類型。

表4 黑龍江省大興安嶺重點國有林區(qū)2000—2020年復合系統耦合協調狀況

3.3 復合系統耦合協調度趨勢預測

經GM(1,1)模型計算，得到黑龍江省大興安嶺重點國有林區(qū)2000—2020年耦合協調度擬合值及其相對誤差(表5)。依據范圍集中原則，將GM(1,1)模型相對誤差劃分為5個狀態(tài)區(qū)間Z1～Z5，各狀態(tài)區(qū)間上、下限依據實際情況進行劃分(表6)。根據GM(1,1)模型相對誤差和狀態(tài)區(qū)間計算狀態(tài)轉移概率矩陣，采用馬爾科夫過程對GM(1,1)模型擬合值進行修正，得到灰色馬爾科夫模型擬合值。通過對比2類模型精度檢驗結果發(fā)現，灰色馬爾科夫模型擬合值在平均相對誤差、后驗差比值和小概率誤差3個方面均優(yōu)于GM(1,1)模型，表明灰色馬爾科夫模型擬合效果更優(yōu)(表7)。通過1～2步狀態(tài)轉移概率矩陣，預測2021—2022年GM(1,1)模型預測值相對誤差的狀態(tài)區(qū)間均為Z3，據此修正2021—2022年GM(1,1)模型耦合協調度預測值得到灰色馬爾科夫模型耦合協調度預測值。

表5 GM(1,1)模型與灰色馬爾科夫模型擬合結果

表6 耦合協調度GM(1,1)模型預測值相對誤差狀態(tài)劃分

由表8可知，經預測，2021—2022年黑龍江省大興安嶺重點國有林區(qū)復合系統耦合協調度總體呈現逐漸優(yōu)化趨勢，到2022年復合系統將恢復至良好協調發(fā)展類型。

表7 耦合協調度預測模型精度對比

表8 2021—2022年黑龍江省大興安嶺重點國有林區(qū)耦合協調度預測情況

4 討論

國家重大政策實施對黑龍江省大興安嶺重點國有林區(qū)耦合協調發(fā)展狀況具有重要影響。自2000年開始，各子系統層發(fā)展指數和復合系統綜合發(fā)展指數在一段時期內均有所增長，但在初期增長幅度有限，耦合協調發(fā)展水平尚未進入協調發(fā)展階段。雖然在天然林資源保護工程一期等一系列舉措的導向下，林區(qū)資源逐步得以恢復，有效推動了國有林區(qū)協調發(fā)展進程。但是，受限于歷史遺留問題疊加、政策紅利未充分釋放等因素制約，耦合協調發(fā)展水平依然處于較低階段。2010年之后，各子系統層發(fā)展指數和復合系統綜合發(fā)展指數均在短期內呈現較大幅度增長，耦合協調發(fā)展水平提升至協調階段。天然林資源保護工程二期、大小興安嶺生態(tài)保護與經濟轉型規(guī)劃(2010—2020年)等政策措施的持續(xù)加持，保障了要素投入水平和政策扶持力度得以逐步提高，進一步夯實了國有林區(qū)發(fā)展基礎，加速了耦合協調發(fā)展進程。但自2014年前后開始，除民生福祉系統層外，其他子系統層發(fā)展指數及復合系統綜合發(fā)展指數均有所回落，尤其是企業(yè)管理和支撐保障系統層下降幅度明顯?？梢姡嫱Ｖ固烊涣稚虡I(yè)性采伐、國有林區(qū)改革指導意見等政策的實施，推動了國有林區(qū)由限采轉為停伐、由政企合一轉為“四分開”，徹底終止了國有林區(qū)傳統主流的生產方式和經營方式，這給國有林區(qū)協調發(fā)展水平的持續(xù)提升帶來暫時的震蕩影響。

各子系統層間的耦合程度與協調程度并非呈線性相關，只有當耦合狀態(tài)達到高度有序時才會形成良性協調，而復合系統耦合協調水平受制于其中最薄弱的子系統層[13]。該研究定量分析結果和筆者實地調研的情況均表明，企業(yè)管理和支撐保障系統層已成為制約黑龍江省大興安嶺國有林區(qū)高質量發(fā)展的短板，尤其是支撐保障方面建設滯后較為嚴重。當前，林區(qū)基礎設施建設不平衡、不充分的問題較為突出，森林經營管護專用道路總量不足、等級不高、線路不暢，部分道路橋涵年久失修、損毀嚴重，已無法通行；通訊信號尚未全覆蓋，監(jiān)測巡護手段較為落后，森林防火等工作依然十分依賴人工瞭望監(jiān)測的傳統方式，高新現代技術、實用技術成果推廣應用不足。此外，國有林區(qū)至今缺乏基礎設施建設方面的專項資金，再加上停伐后大興安嶺林業(yè)集團公司自籌資金能力嚴重不足，導致林區(qū)在支撐保障方面陷入新建不足、有建無管、老舊棄管的局面。如果支撐保障系統層長期得不到改善，將會嚴重影響國有林區(qū)資源的監(jiān)管效率、產業(yè)的轉型發(fā)展、企業(yè)運營管理的革新和民生福祉的持續(xù)改善。另外，值得注意的是，雖然從縱向比較來看，產業(yè)發(fā)展系統層呈現逐年向好勢態(tài)，但通過橫向比較杜鈺瑋等[26]關于其他國有林區(qū)以及王震等[27]關于非國有林區(qū)林業(yè)產業(yè)發(fā)展情況的研究結果來看，黑龍江省大興安嶺國有林區(qū)林業(yè)經濟增長狀況、生產力發(fā)展能力和產業(yè)結構水平皆有待進一步改善。

隨著疫情得到有效控制和國內經濟逐步復蘇，在保持現有各項支持政策穩(wěn)定的前提下，預計到2022年，黑龍江省大興安嶺國有林區(qū)復合系統耦合協調發(fā)展類型將恢復至良好協調發(fā)展狀態(tài)，但距離優(yōu)質耦合協調狀態(tài)尚有較大差距。因此，建議在“十四五”期間繼續(xù)保持天然林保護等各項重點工程和重大政策的持續(xù)投入，確保森林資源和自然生態(tài)系統健康穩(wěn)定；加快產業(yè)結構優(yōu)化步伐，著力破解國有林區(qū)接替產業(yè)同質化發(fā)展存在的隱患[28]，瞄準如中藥材、生態(tài)旅游和森林康養(yǎng)等1～2個潛力較大的接替產業(yè)，集中要素扶持培育，突出產業(yè)和產品特色，主動融入和適應社會主義市場經濟大環(huán)境，增強抵御市場風險的能力，增強產業(yè)發(fā)展后勁。此外，更需要引導大興安嶺林業(yè)集團公司加快構建現代企業(yè)治理體系，完善和創(chuàng)新企業(yè)管理機制和運營機制，逐步提高企業(yè)治理能力和治理效率，持續(xù)釋放國有林區(qū)改革紅利。同時，爭取設立國有林區(qū)支撐保障建設專項資金，并將其分列為新建類專項資金和運營維護類專項資金，強化支撐保障能力。

5 結論

經過2000—2020年21 a的轉型發(fā)展，黑龍江省大興安嶺重點國有林區(qū)復合系統各子系統層發(fā)展指數均有所改善，但是具體發(fā)展變化過程各有差異，其中，民生福祉和資源保育2個子系統層增幅最為明顯；國有林區(qū)復合系統耦合協調發(fā)展階段已由失調階段轉變?yōu)閰f調階段，但是受新冠肺炎疫情等多種因素影響，協調發(fā)展水平不高、狀態(tài)不穩(wěn)；預測到2022年，國有林區(qū)復合系統耦合協調發(fā)展類型將恢復至良好協調發(fā)展狀態(tài)，但距離優(yōu)質耦合協調狀態(tài)尚有較大差距。建議在“十四五”期間，黑龍江省大興安嶺重點國有林區(qū)需要緊緊抓住天然林資源保護修復和接續(xù)替代產業(yè)發(fā)展2大主線，爭取現有扶持政策和要素投入水平保持穩(wěn)定，加快補齊企業(yè)管理和支撐保障系統層的短板，推動國有林區(qū)盡早步入穩(wěn)定的高水平耦合協調發(fā)展狀態(tài)。