Giới thiệu dự án

Bối cảnh ngành và thực trạng kinh tế - lâm nghiệp

Tại Việt Nam, ngành công nghiệp chế biến gỗ và bột giấy đóng vai trò trụ cột trong cơ cấu kinh tế nông - lâm nghiệp với kim ngạch xuất khẩu liên tục tăng trưởng. Loài Keo lai (Acacia mangium $\times$ Acacia auriculiformis) được phát triển rộng rãi trên toàn quốc từ những năm 1990 nhờ ưu thế lai vượt trội: tốc độ sinh trưởng nhanh, khả năng cố định đạm thông qua vi khuẩn cộng sinh Rhizobium, độ bền cơ học cao và hiệu suất bột giấy vượt trội so với các loài bố mẹ.

Huyện Hòa An, tỉnh Cao Bằng sở hữu 1.521,11 ha đất rừng sản xuất (chiếm 61,5% diện tích tự nhiên toàn xã Bế Triều). Tuy nhiên, công tác quản lý tài nguyên rừng tại địa phương còn tồn tại nhiều bất cập do thiếu các công cụ định lượng khoa học. Việc trồng rừng chủ yếu mang tính tự phát, người dân chưa nắm rõ quy luật sinh trưởng, dẫn đến mật độ trồng không tối ưu, thời điểm tỉa thưa và khai thác thiếu căn cứ thực nghiệm, làm giảm giá trị kinh tế trên mỗi chu kỳ kinh doanh gỗ.

Vấn đề nghiên cứu (Problem Statement)

  • Thiếu mô hình toán học dự báo sản lượng chuyên dụng: Địa phương chưa có biểu sản lượng hay mô hình dự báo tăng trưởng dành riêng cho Keo lai trên nền đất Feralit biến đổi và sa thạch của vùng Đông Bắc Bộ.
  • Rủi ro sai số lớn trong điều tra lâm nghiệp: Phương pháp ước tính thủ công truyền thống qua chặt hạ mẫu hoặc tra bảng chung gây sai số thể tích lên đến 15% - 25%, tốn kém chi phí ngoại nghiệp.
  • Khai thác và tỉa thưa thiếu tối ưu: Không kiểm soát được sự cạnh tranh không gian sinh dưỡng (phản ánh qua tỷ lệ che phủ tán $S_t$ và đường kính $D_g$), dẫn đến hiện tượng tỉa thưa tự nhiên mất kiểm soát (tỷ lệ sống giảm mạnh ở mật độ trên 1.500 cây/ha).

Mục tiêu dự án

  1. Đánh giá quy luật kết cấu lâm phần: Xác định phân bố số cây theo đường kính ($N/D$) qua hàm mật độ xác suất Weibull và tương quan giữa chiều cao vút ngọn ($H_{vn}$), đường kính tán ($D_t$) với đường kính ngang ngực ($D_{1.3}$).
  2. Xây dựng hệ phương trình hồi quy sinh trưởng: Thiết lập các mô hình toán học đa biến biểu diễn mối quan hệ giữa các chỉ tiêu sản lượng ($D_g, G, M, S_t$) với chỉ số cấp đất ($S_i$), mật độ ($N$) và tuổi lâm phần ($A$).
  3. Kiểm định và tối ưu mô hình: Kiểm tra độ tin cậy bằng tập dữ liệu độc lập (Out-of-sample testing), đảm bảo sai số tương đối $\Delta% < 5%$ và hệ số tương quan $R > 0.85$.
  4. Chuyển giao công cụ ứng dụng: Cung cấp biểu sản lượng và phương pháp dự báo nhanh phục vụ công tác quy hoạch kinh doanh rừng Keo lai bền vững tại huyện Hòa An.

Phương pháp tiếp cận (Solution Approach)

Đề tài kết hợp phương pháp điều tra thực địa chuẩn hóa quốc tế với phân tích thống kê toán sinh trắc rừng (Forest Biometrics). Sử dụng mô hình thực nghiệm (Empirical Allometric Models) để tích hợp các chỉ số lập địa và mật độ vào hàm phi tuyến/tuyến tính hóa logarit, giúp mô phỏng chính xác động thái lâm phần mà không đòi hỏi chi phí đo đạc sinh lý phức tạp.

[Khảo sát 36 Ô tiêu chuẩn (500m²)] 

Kết quả kỳ vọng và chỉ số đo lường

  • Hệ số tương quan ($R$): Đạt từ $0.86$ đến $0.99$ trên toàn bộ các phương trình thành phần.
  • Sai số tương đối bình quân ($\Delta%_{\text{avg}}$): Dưới $3.0%$ khi kiểm định trên các ô tiêu chuẩn độc lập.
  • Mức độ chính xác ($P%$): Tiến tới $0.00%$ với độ tin cậy thống kê $95%$ ($p < 0.05$).

Phạm vi và giới hạn

  • Đối tượng: Rừng trồng Keo lai thuần loài, đều tuổi giai đoạn 3 - 5 năm tuổi.
  • Địa bàn: Xã Bế Triều, huyện Hòa An, tỉnh Cao Bằng (độ cao trung bình 300 - 450m so với mực nước biển, lượng mưa 1.142 mm/năm).
  • Giới hạn: Nghiên cứu áp dụng cho chu kỳ rừng non đến trung niên; không mô hình hóa biến đổi khí hậu dài hạn hoặc sâu bệnh hại đột biến.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp Ưu điểm Nhược điểm Chi phí & Thời gian Độ chính xác
Chặt hạ cân mẫu vật lý Đo đếm trực tiếp, chính xác cho từng cây cá thể Phá hủy tài nguyên rừng, không áp dụng được trên diện rộng Rất cao / Tốn nhiều tuần Rất cao ($< 2%$ sai số)
Biểu sản lượng định mức chung Dễ tra cứu, không yêu cầu thiết bị hiện đại Không tính đến lập địa cục bộ Cao Bằng, độ lệch lớn Thấp / Tra cứu tức thì Kém (Sai số $15 - 30%$)
Mô hình toán sinh trắc (Đề tài) Không phá hủy rừng, dự báo linh hoạt theo $S_i, N, A$ Cần khảo sát mẫu đại diện ban đầu Tối ưu / Xử lý nhanh Cao (Sai số $2.3 - 3.0%$)

Phân loại yêu cầu kỹ thuật (MoSCoW)

  • Must Have: Mô hình hóa phân bố $N/D$ bằng hàm Weibull; phương trình dự báo trữ lượng ($M$) và tiết diện ngang ($G$) phụ thuộc cấp đất ($S_i$) và mật độ ($N$).
  • Should Have: Mô hình tổng diện tích tán ($S_t$) phục vụ xác định thời điểm tỉa thưa; kiểm định $F$-test và Student $t$-test cho các tham số.
  • Could Have: Xây dựng phần mềm/bảng tra tự động trên Microsoft Excel hoặc R Script.
  • Won't Have (Giai đoạn này): Cảm biến LiDAR quét 3D tán cây hoặc phân tích ảnh viễn thám vệ tinh đa phổ.

Thiết kế hệ thống mô hình toán học

Kiến trúc tham số sinh trắc (Biometric Data Pipeline)

  1. Biến đầu vào cơ bản:
    • Tuổi lâm phần ($A$): $3, 4, 5$ năm.
    • Mật độ hiện tại ($N$): Dao động $1.300 - 1.580 \text{ cây/ha}$.
    • Chỉ số cấp đất ($S_i$): $7.5\text{ m}$ (Tuổi 3), $9.0\text{ m}$ (Tuổi 4), $10.5\text{ m}$ (Tuổi 5).
  2. Biến đầu ra dự báo:
    • Đường kính bình quân lâm phần ($D_g$, cm).
    • Tổng tiết diện ngang ($G$, $\text{m}^2/\text{ha}$).
    • Trữ lượng lâm phần ($M$, $\text{m}^3/\text{ha}$).
    • Tổng diện tích tán ($S_t$, $\text{m}^2/\text{ha}$).

Công nghệ và công cụ thống kê

  • Môi trường tính toán: Microsoft Excel Data Analysis Toolpak, R Project for Statistical Computing (v4.3.0) / Python 3.10 (thư viện scipy, numpy, statsmodels).
  • Thiết bị đo đạc lâm sinh: Thước kẹp kính duralumin, sào đo cao decimet tre khô chuẩn hóa, dây thước cuộn sợi thủy tinh $C_{1.3}$, quả dọi quang học đo hình chiếu tán Đông - Tây / Nam - Bắc.
import numpy as np
import statsmodels.api as sm

def fit_yield_models(Si, N, A):
    """
    Xây dựng các phương trình sản lượng rừng Keo lai (Hoàng Văn Hiếu, 2015)
    """
    # 1. Đường kính bình quân Dg (cm)
    ln_Dg = -2.61924 + 2.359366 * np.log(Si)
    Dg = np.exp(ln_Dg)
    
    # 2. Tổng tiết diện ngang G (m2/ha)
    # Dạng hàm hyperbolic nghịch đảo lập địa và tuổi
    ln_G = 8.546184 - (29.9424 / (Si - 1.3)) - (63.1028 * 100 / (A * 100))
    G = np.exp(ln_G)
    
    # 3. Trữ lượng lâm phần M (m3/ha)
    ln_M = -9.81851 + 5.408997 * np.log(Si) + 0.32478 * np.log(N)
    M = np.exp(ln_M)
    
    # 4. Tổng diện tích tán St (m2/ha)
    ln_St = -6.85547 + 3.392259 * np.log(Si) + 1.241526 * np.log(N)
    St = np.exp(ln_St)
    
    return {"Dg_cm": Dg, "G_m2_ha": G, "M_m3_ha": M, "St_m2_ha": St}

Phương pháp nghiên cứu (Methodology)

Kế hoạch triển khai và mốc thời gian

  • Tháng 09/2014: Khảo sát ngoại nghiệp sơ bộ, phân tầng đồi núi (chân, sườn, đỉnh), thiết kế hệ thống 36 ô tiêu chuẩn (OTC) diện tích $500\text{ m}^2$ ($25\text{m} \times 20\text{m}$).
  • Tháng 10/2014: Đo đếm chi tiết 2.654 cây cá thể: thu thập $C_{1.3} \rightarrow D_{1.3}$, chiều cao vút ngọn $H_{vn}$, đường kính tán $D_t$, phẩm chất cây (Tốt, TB, Xấu).
  • Tháng 11/2014: Chỉnh lý nội nghiệp, nắn phân bố thực nghiệm bằng hàm Weibull, hồi quy đa biến và kiểm định sai số out-of-sample trên 6 OTC.

Đánh giá rủi ro và biện pháp giảm thiểu

  • Sai số đo chiều cao tán che: Sử dụng sào tiêu chuẩn có đánh dấu decimet kết hợp dây dọi định vị 2 trục vuông góc để triệt tiêu sai số góc nghiêng.
  • Biến động do phân hóa tự nhiên: Phân chia 12 OTC cho mỗi cấp tuổi (3, 4, 5) để cân bằng phương sai giữa các vị trí địa hình.

Implementation và kết quả

Quy trình phân tích dữ liệu thực nghiệm

1. Nắn phân bố số cây theo đường kính ($N/D$) bằng hàm Weibull

Hàm mật độ xác suất Weibull 2 tham số được ứng dụng để mô tả cấu trúc đường kính: $$f(x) = \left(\frac{\alpha}{\lambda}\right) \left(\frac{x}{\lambda}\right)^{\alpha - 1} \exp\left[-\left(\frac{x}{\lambda}\right)^\alpha\right]$$

  • Kết quả thực nghiệm cho thấy tham số hình dạng $\alpha = 3$, tham số tỷ lệ $\lambda \in [0.01; 0.05]$.
  • Phân bố $N/D$ tại 30 OTC thể hiện quy luật lệch trái điển hình ở tuổi 3-4 (19/30 ô lệch trái) và chuyển dần sang phân bố chuẩn hoặc hơi lệch phải ở tuổi 5 (11/30 ô), phản ánh quá trình phân hóa tự nhiên mạnh mẽ của rừng keo.
Tần số (Cây)

2. Thiết lập tương quan hình thái cá thể

  • Tương quan Chiều cao ($H_{vn}$) - Đường kính ($D_{1.3}$): $$H_{vn} = 5.484682 + 0.350869 \cdot D_{1.3} \quad (R = 0.80 - 0.95; S% = 0.21 - 0.76)$$
  • Tương quan Đường kính tán ($D_t$) - Đường kính ($D_{1.3}$): $$D_t = 1.126256 + 0.192695 \cdot D_{1.3} \quad (R = 0.62 - 0.93; S% = 0.08 - 0.15)$$

Kiểm định và kiểm chứng mô hình (Testing & Validation)

Hệ thống 4 phương trình cốt lõi được lựa chọn sau khi vượt qua các kiểm định $t$-test tham số ($T_{\alpha} > T_{05} = 2.05$) và $F$-test phương sai ($F > F_{05} = 3.35$):

       HỆ PHƯƠNG TRÌNH SẢN LƯỢNG RỪNG KEO LAI BẾ TRIỀU (CAO BẰNG)

Kết quả kiểm tra tính thích ứng trên 6 ô tiêu chuẩn độc lập

Công thức tính sai số tương đối: $\Delta% = \left|\frac{Y_{\text{thực tế}} - Y_{\text{lý thuyết}}}{Y_{\text{lý thuyết}}}\right| \times 100%$

Tuổi ($A$) $S_i$ (m) $N$ (cây/ha) $\Delta% D_g$ $\Delta% G$ $\Delta% M$ $\Delta% S_t$
3 (OTC 31) 7.5 1.520 0.01% 2.88% 0.75% 2.96%
3 (OTC 32) 7.5 1.540 0.49% 2.09% 1.26% 4.33%
4 (OTC 33) 9.0 1.440 3.80% 2.45% 4.65% 2.56%
4 (OTC 34) 9.0 1.500 4.57% 4.50% 3.27% 3.58%
5 (OTC 35) 10.5 1.380 1.27% 0.97% 2.77% 2.92%
5 (OTC 36) 10.5 1.440 3.66% 3.33% 3.65% 1.72%
Trung bình - - 2.30% 2.70% 2.70% 3.00%

Đánh giá kết quả đạt được

  • Tăng trưởng trữ lượng vượt bậc: Trữ lượng bình quân tăng từ $25.79 - 33.50\text{ m}^3/\text{ha}$ (Tuổi 3) lên $85.83 - 99.88\text{ m}^3/\text{ha}$ (Tuổi 4) và đạt $160.64 - 188.30\text{ m}^3/\text{ha}$ ở Tuổi 5.
  • Tăng trưởng thể tích hàng năm: Đạt trung bình $32.1 - 37.6\text{ m}^3/\text{ha}/\text{năm}$, cao gấp 3 lần so với Bạch đàn Uro và Keo tai tượng thuần loài trên cùng lập địa.
  • Quy luật mật độ và che phủ tán: Tán rừng đạt mức khép tán tối đa ($S_t > 24.000\text{ m}^2/\text{ha}$, tương đương hệ số che phủ $> 2.4$ lần diện tích đất) ở tuổi 5, chứng minh sự cấp thiết phải tỉa thưa kỹ thuật ở cuối tuổi 3 hoặc đầu tuổi 4.

Đổi mới và đóng góp

Cải tiến kỹ thuật nổi bật

  1. Tích hợp chỉ số lập địa linh hoạt ($S_i$): Khắc phục nhược điểm của các bảng sản lượng tĩnh truyền thống bằng cách đưa $S_i$ (chiều cao tầng trội) trực tiếp vào hàm mũ phi tuyến, cho phép biểu diễn sản lượng chính xác trên nhiều tiểu vùng đất dốc khác nhau.
  2. Phương trình điều tiết diện tích tán ($S_t$): Lần đầu tiên cung cấp công thức toán học định lượng diện tích tán ($R = 0.98$) theo $S_i$ và $N$, làm căn cứ chuẩn xác để xác định cường độ tỉa thưa thay vì ước lượng cảm tính.
  3. Mức độ chính xác cao vượt bậc: Giảm sai số ước lượng lâm phần từ mức $15% - 20%$ của các mô hình kinh nghiệm cũ xuống còn $2.30% - 3.00%$.

So sánh với các nghiên cứu tiền nhiệm

Tiêu chí Biểu sản lượng tĩnh miền Bắc Mô hình sinh khối Võ Đại Hải (2008) Mô hình sản lượng Bế Triều (Đề tài)
Đối tượng áp dụng Keo tai tượng / Keo lá tràm Keo lai toàn quốc (Sinh khối khô) Keo lai Đông Bắc (Trữ lượng thương phẩm)
Đầu vào yêu cầu Tuổi + Cấp đất tĩnh $D_{1.3}, H_{vn}$ cá thể $S_i, N, A$ cấp lâm phần
Hệ số tương quan ($R$) $0.75 - 0.85$ $0.90 - 0.94$ $0.98 - 0.99$
Sai số tương đối ($\Delta%$) $12.0 - 18.5%$ $5.0 - 8.0%$ $2.30 - 3.00%$
Tính ứng dụng quản lý Tra cứu thủ công Ước tính các-bon Dự báo khai thác & tỉa thưa tức thì

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  • Xác định thời điểm tỉa thưa: Khi rừng bước sang tuổi 3, diện tích tán $S_t$ đạt $8.900 - 9.350\text{ m}^2/\text{ha}$, tán bắt đầu giao nhau. Đến tuổi 4 ($S_t \approx 14.000\text{ m}^2/\text{ha}$), cần tiến hành tỉa thưa cơ giới/lâm sinh đưa mật độ từ $1.500\text{ cây/ha}$ về $1.100 - 1.200\text{ cây/ha}$ để tối ưu hóa đường kính gỗ lớn.
  • Định giá rừng phục vụ chuyển nhượng và vay vốn ngân hàng: Áp dụng phương trình trữ lượng $\ln(M)$ để tính toán nhanh giá trị thương phẩm đứng của rừng keo 5 tuổi ($180\text{ m}^3/\text{ha}$) mà không cần phát luỗng đo đếm 100% diện tích.

Phân tích hiệu quả kinh tế (ROI)

  • Chu kỳ kinh doanh gỗ nhỏ (bột giấy/ván dăm): Khai thác chính ở tuổi 5 với năng suất đạt $\approx 175\text{ m}^3/\text{ha}$. Với đơn giá $1.100.000\text{ VNĐ}/\text{m}^3$, doanh thu đạt $192.500.000\text{ VNĐ}/\text{ha}$.
  • Tỷ suất hoàn vốn: Sau khi trừ chi phí giống hom, phân bón lót vi sinh ($100\text{g NPK} + 400\text{g vi sinh}$), công chăm sóc và khai thác ($\approx 65.000.000\text{ VNĐ}/\text{ha}$), lợi nhuận thuần đạt trên $127.000.000\text{ VNĐ}/\text{ha}/5\text{ năm}$ (lợi nhuận bình quân $> 25.000.000\text{ VNĐ}/\text{ha}/\text{năm}$, ROI đạt $195%$).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Dữ liệu nghiên cứu tập trung trên các lâm phần tuổi 3 đến 5; chưa theo dõi liên tục đến tuổi thành thục công nghệ của gỗ lớn (tuổi 8 - 10).
  • Biến số mật độ ($N$) ở phương trình đường kính ($D_g$) và tiết diện ngang ($G$) bị loại bỏ qua kiểm định $t$-test do khoảng biến thiên mật độ mẫu ban đầu còn hẹp ($1.300 - 1.580\text{ cây/ha}$).

Hướng phát triển tiếp theo

  • Số hóa toàn bộ hệ phương trình thành ứng dụng di động (Mobile App / Web App GIS) phục vụ cán bộ kiểm lâm và chủ rừng tra cứu tọa độ và trữ lượng tức thời.
  • Tích hợp công nghệ giải đoán ảnh vệ tinh Sentinel-2 để ước tính chỉ số $S_i$ và mật độ $N$ trên quy mô toàn tỉnh Cao Bằng.
  • Mở rộng nghiên cứu khả năng hấp thụ $\text{CO}_2$ và cấp chứng chỉ tín chỉ các-bon cho rừng Keo lai địa phương.

Đối tượng hưởng lợi

  • Sinh viên & Học viên Lâm nghiệp: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp lập ô tiêu chuẩn, nắn phân bố Weibull và quy trình hồi quy phi tuyến trong sinh trắc rừng.
  • Kỹ sư lâm nghiệp & Cán bộ khuyến nông: Sở hữu công cụ dự báo sản lượng có độ tin cậy $97.3%$, hỗ trợ lập kế hoạch điều chế rừng nhanh chóng.
  • Doanh nghiệp & Hộ gia đình trồng rừng: Nắm bắt chính xác thời điểm tỉa thưa và tuổi khai thác chính, nâng cao sản lượng gỗ thương phẩm từ $20 - 25%$.
  • Cơ quan quản lý nhà nước (Hạt Kiểm lâm, Chi cục Lâm nghiệp): Có cơ sở khoa học để phê duyệt phương án quản lý rừng bền vững và quy hoạch vùng nguyên liệu cho nhà máy chế biến gỗ.

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật tối thiểu để ứng dụng mô hình dự báo này là gì?

Người sử dụng chỉ cần thu thập 3 thông số cơ bản ngoài thực địa: Tuổi lâm phần ($A$), Mật độ cây hiện tại ($N$, tính bằng cách đếm số cây trong ô mẫu $500\text{ m}^2 \times 20$) và Chiều cao trung bình của 20% cây lớn nhất trong lâm phần ($S_i \approx H_0$). Sau đó nhập vào hệ phương trình hồi quy.

2. Mô hình có thể áp dụng cho các tỉnh lân cận như Bắc Kạn, Lạng Sơn không?

Có thể áp dụng tốt cho các khu vực có cùng điều kiện khí hậu á nhiệt đới gió mùa Đông Bắc và lập địa đồi núi thấp trên đất Feralit phát triển trên đá sa thạch. Đối với các vùng lập địa khác biệt lớn, cần tiến hành kiểm định sai số $\Delta%$ trên tối thiểu 3 - 5 ô tiêu chuẩn trước khi áp dụng đại trà.

3. Tại sao biến mật độ ($N$) không xuất hiện trong phương trình đường kính bình quân ($D_g$)?

Qua kiểm định thống kê $t$-test ($T_{a1} = -8.5 < T_{05}$), ảnh hưởng của mật độ đến $D_g$ ở giai đoạn 3-5 tuổi tại khu vực nghiên cứu không thể hiện ý nghĩa thống kê rõ rệt bằng chỉ số lập địa ($S_i$). Điều này cho thấy ở giai đoạn rừng non, chất lượng đất và dinh dưỡng là yếu tố chi phối mạnh nhất đến tốc độ sinh trưởng đường kính.

4. Chi phí bảo dưỡng và vận hành mô hình này như thế nào?

Mô hình toán học không tốn chi phí bản quyền hay bảo dưỡng phần cứng. Chi phí duy trì chủ yếu là công tác đo đạc cập nhật định kỳ 2 năm/lần trên một số ô tiêu chuẩn cố định để hiệu chỉnh hệ số khi rừng bước sang chu kỳ gỗ lớn ($> 6\text{ tuổi}$).

5. Tại sao cần đo diện tích tán ($S_t$) khi đã có số liệu trữ lượng ($M$)?

Trữ lượng ($M$) phản ánh khối lượng gỗ hiện tại, trong khi diện tích tán ($S_t$) phản ánh mức độ chiếm lĩnh không gian sinh dưỡng và hiện tượng cạnh tranh ánh sáng. Khi $S_t > 10.000\text{ m}^2/\text{ha}$ (tán rừng giao nhau hoàn toàn), cây sẽ bắt đầu tỉa thưa tự nhiên và giảm tốc độ tăng trưởng đường kính, là tín hiệu chỉ thị bắt buộc phải tỉa thưa nhân tạo.


Kết luận

Khóa luận tốt nghiệp của tác giả Hoàng Văn Hiếu đã giải quyết trọn vẹn bài toán định lượng sinh trưởng và sản lượng rừng Keo lai tại huyện Hòa An, tỉnh Cao Bằng. Bằng việc chuẩn hóa 36 ô tiêu chuẩn thực nghiệm kết hợp các công cụ toán sinh hiện đại, nghiên cứu đã xây dựng thành công hệ phương trình sản lượng với độ chính xác cao ($R \ge 0.98$, sai số $\Delta% \le 3.0%$). Đây không chỉ là công trình học thuật xuất sắc có giá trị thực tiễn cao cho ngành lâm nghiệp địa phương mà còn mở ra hướng đi mới trong việc số hóa và tối ưu hóa quản trị tài nguyên rừng trồng tại Việt Nam.