遼寧省“十四五”森林合理年采伐量測算方法

姜黎黎

（遼寧省林業調查規劃監測院，沈陽 110122）

森林采伐限額制度是《森林法》確立的一項重要的法律制度。森林采伐限額編制是以保護發展森林資源、改善生態環境和推進生態文明、建設美麗中國為宗旨，以森林可持續經營理論為依據，以培育健康、穩定、高效的森林生態系統和提供更多更好的優質林產品為目標。編制限額和實施限額的目的是不斷增加森林資源數量和提升森林資源質量，穩步增強森林生態系統的整體功能。森林合理年采伐量測算是森林采伐限額編制的重點內容、核心內容，如何科學地測算采伐限額主伐合理年伐量，將決定著編限單位森林采伐限額編制的成果質量和準確性，意義十分重大[1]。

遼寧省十四五森林采伐限額編制，以森林資源管理“一張圖”年度更新成果為基礎，嚴格按照采伐限額的編制原則，區分不同森林類別、樹種、林分起源和經營目標，以促進森林結構調整和提高林分質量為目的，合理年采伐量測算是遼寧省“十四五”森林采伐限額編制的重中之重。為提高合理年采伐量測算效率和測算準確性，研究探討了根據國家林草局《關于編制“十四五”期間年森林采伐限額編制工作的通知》（林資發[2019]99 號）的要求，結合遼寧省實際，充分利用國家林草局下發的森林合理年伐量測算系統測算合理年采伐量，通過Py-thon 編寫程序實現森林合理年采伐量測算數據提取及測算方法。

1 界面設計

利用Python 的Tkinter 模塊，設計工具界面，主要使用了Entry、Button、Label 等控件，如圖1。

圖1

2 Access 數據庫設計

數據庫是測算數據的存儲和提取的容器，數據庫中一共需要兩張數據表：測算數據小班表，測算數據中間表。測算數據小班表用來存儲導入的小班表（Excel 格式），測算數據中間表用來存儲測算數據小班表標準化處理后的結果，也是數據提取的基礎表。

3 測算數據提取流程

3.1 測算數據小班表（Excel 格式）導入

將各編限單位按照編限要求填寫好編限單位名稱、編限單位代碼、是否參與編限、編限權屬的小班序列表導入到Access 數據庫中的測算數據小班表中。主要是使用Python 的xlrd 模塊讀取Excel 文件，用pypyodbc 模塊寫入Access 數據庫。xlrd 是Py-thon 讀取Excel 的擴展模塊，可以實現指定表單、指定單元格的讀取。pypyodbc 是一個Python ODBC 庫，具有極大的兼容性、可嵌入性和代碼移植性，可在多平臺下運行。

3.2 測算數據預處理

將測算數據小班表標準化處理，生成數據中間表。

3.3 編限類型劃分

對數據中間表劃分編限類型。

3.3.1 編限類型分類

根據國家林草局下發的森林合理年伐量測算系統需要的6 張數據表，對參與編限的森林資源數據進行編限類型劃分，一共劃分為7 個大類型：一般用材林主伐、短輪伐用材林主伐、公益林更新采伐、撫育采伐、低產(效)林改造、其他采伐、不參與編限，具體劃分情況如表1。

表1 森林資源數據編限類型劃分

3.3.2 確定編限類型劃分條件

根據《森林經營技術規程》DB21/T 706—2013結合遼寧省森林經營實際，確定各編限類型的劃分條件，讓每個參與編限測算的小班數據都能準確的劃分編限類型。

3.3.3 使用Python+pypyodbc，利用SQL 的UPDATE語句對數據中間表劃分編限類型。

如公益林生態疏伐的劃分語句為：

UPDATE XIAOBAN SET編限類型=′44′WHERE森林類別代碼in(′1′,′3′,′4′)And 齡組代碼=′2′And 郁閉度>=0.8

3.4 合理年伐量測算數據提取

按照森林合理年伐量測算系統的要求提取需要的測算數據表，以一般用材林主伐為例（即01.txt），如下：

SELECT * FROM (SELECT 編限單位代碼,′233′as 林種代碼,szdm,numb,起源代碼,林木權屬,

ROUND(SUM(IIF(齡組代碼=′1′,面積,0)),2) AS幼齡林面積,SUM(IIF(齡組代碼=′1′,小班蓄積,0))AS 幼齡林蓄積,

ROUND(SUM(IIF(齡組代碼=′2′,面積,0)),2) AS中齡林面積,SUM(IIF(齡組代碼=′2′,小班蓄積,0))AS 中齡林蓄積,

ROUND(SUM(IIF(齡組代碼=′3′,面積,0)),2) AS近熟林面積,SUM(IIF(齡組代碼=′3′,小班蓄積,0))AS 近熟林蓄積,

ROUND(SUM(IIF(齡組代碼=′4′,面積,0)),2) AS成熟林面積,SUM(IIF(齡組代碼=′4′,小班蓄積,0))AS 成熟林蓄積,

ROUND(SUM(IIF(齡組代碼=′5′,面積,0)),2) AS過熟林面積,SUM(IIF(齡組代碼=′5′,小班蓄積,0))AS 過熟林蓄積,TZMB

FROM (SELECT 編限單位代碼,(select code from codetable where name=trim(優勢樹種)&起源代碼)AS szdm,(select num from codetable where name=trim(優勢樹種) &起源代碼) AS numb,(select mubiao from codetable where name=trim(優勢樹種)&起源代碼)AS tzmb,地類代碼,森林類別代碼,起源代碼,林種代碼,郁閉度,齡組代碼,面積,小班蓄積,林木權屬FROM XIAOBAN WHERE 編限類型in (′10′,′41′,′42′) AND 林種代碼in (′233′) AND 起源代碼=′2′) GROUP BY 編限單位代碼,林木權屬,起源代碼,szdm,numb,TZMB)

WHERE 幼齡林蓄積+中齡林蓄積+近熟林蓄積+成熟林蓄積+過熟林蓄積>0

使用SQL 多表嵌套查詢語句，利用IIF、SUM 函數等提取一般用材林主伐測算數據，最后一個WHERE條件很重要，經過測試發現，如果符合一般用材林主伐條件的某樹種，幼齡林蓄積+中齡林蓄積+近熟林蓄積+成熟林蓄積+過熟林蓄積=0的情況下，導入森林合理年伐量測算系統會報錯，導致導入失敗，所以必須通過WHERE 條件篩選掉等于0的情況。

3.5 合理年采伐量測算數據提取

使用Python 自帶的I/O 接口，利用for 循環構建函數逐行讀取提取結果并導出為文本文件(01-06.txt)。代碼如下：

4 森林合理年伐量測算

將導出的文本文件(01-06.txt)，導入到森林合理年伐量測算系統中進行測算，計算出該單位的合理年伐量。

綜上所述，通過該方法，實現森林合理年采伐量測算數據的一鍵提取，結合森林合理年伐量測算系統的批量測算，極大地提高了遼寧省“十四五”森林合理年采伐量測算的工作效率，最大程度降低新冠疫情對遼寧省“十四五”森林采伐限額編制工作的影響。