Giới thiệu dự án

Ngành công nghiệp chế biến gỗ và sản xuất bột giấy tại Việt Nam đang duy trì tốc độ tăng trưởng bình quân từ 8 - 12%/năm, kéo theo nhu cầu cấp thiết về nguồn gỗ nguyên liệu rừng trồng có năng suất sinh học cao và phẩm chất đồng đều. Trong cơ cấu cây trồng lâm nghiệp thương mại, giống keo lai tự nhiên (Acacia mangium × Acacia auriculiformis) giữ vị trí chiến lược nhờ ưu thế lai (heterosis) vượt trội: tốc độ tăng trưởng đường kính nhanh gấp 1,3 - 1,8 lần và chiều cao gấp 1,2 - 1,6 lần so với các loài keo bố mẹ, cho năng suất trung bình đạt 20,39 - 23,55 $\text{m}^3/\text{ha/năm}$ tại các vùng sinh thái trọng điểm.

Tuy nhiên, hiệu quả kinh tế và chất lượng gỗ của lâm phần keo lai chịu sự chi phối quyết định bởi mật độ trồng ban đầu. Mật độ không chỉ quyết định trực tiếp đến suất vốn đầu tư ban đầu (cây giống, nhân công, phân bón) mà còn kiểm soát không gian dinh dưỡng, mức độ cạnh tranh ánh sáng, tiến trình tỉa cành tự nhiên và hình thái thân cây. Tại Công ty Lâm nghiệp Hàm Yên (huyện Hàm Yên, tỉnh Tuyên Quang), việc áp dụng hai cấp mật độ 1.333 cây/ha và 1.666 cây/ha vẫn chưa được đánh giá định lượng bằng các phương pháp thống kê sinh trắc học lâm nghiệp ở giai đoạn thành rừng 6 tuổi, dẫn đến sự thiếu hụt căn cứ khoa học cho quyết định kinh doanh gỗ nguyên liệu giấy và gỗ xẻ.

+-----------------------------------------------------------------------------+
|                          MỤC TIÊU VÀ PHẠM VI NGHIÊN CỨU                     |
+-----------------------------------------------------------------------------+
| 1. Đo lường & so sánh đường kính ngang ngực (D1.3) giữa 2 cấp mật độ        |
| 2. Đo lường & so sánh chiều cao vút ngọn (Hvn) và chiều cao dưới cành (Hdc) |
| 3. Khảo sát biến thiên đường kính tán (Dt) và không gian dinh dưỡng         |
| 4. Phân loại phẩm chất lâm phần (Cây Tốt - Trung bình - Xấu)               |
| 5. Kiểm định thống kê sai dị (U-test) và đề xuất mật độ tối ưu             |
+-----------------------------------------------------------------------------+

Dự án xác định 5 mục tiêu cốt lõi:

  1. Định lượng chỉ tiêu sinh trưởng đường kính ngang ngực ($D_{1.3}$) ở tuổi 6 tại mật độ 1.333 cây/ha và 1.666 cây/ha.
  2. Xác định chỉ tiêu sinh trưởng chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) nhằm đánh giá sinh khối và độ thon thân cây.
  3. Phân tích động thái phát triển đường kính tán lá ($D_t$) làm cơ sở xác định mức độ giao tán và hiệu quả che phủ đất dốc.
  4. Đánh giá cấu trúc phẩm chất lâm phần thông qua tỷ lệ phân cấp cây tốt, trung bình và xấu.
  5. Ứng dụng chuẩn thống kê $U$ (U-test) để kiểm định ý nghĩa sai khác giữa các đại lượng đo đếm, từ đó tối ưu hóa mật độ trồng rừng nguyên liệu giấy cho vùng Đông Bắc Bộ.

Phạm vi nghiên cứu: Đề tài giới hạn trên các lô rừng trồng thuần loài keo lai 6 tuổi, lập địa đồi núi thấp (độ cao < 300m, độ dốc < $45^\circ$, đất feralit phát triển trên đá mẹ phiến thạch sét) do Công ty Lâm nghiệp Hàm Yên quản lý.

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

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

Tại khu vực trung du và miền núi phía Bắc, các đơn vị lâm nghiệp thường áp dụng mật độ theo thói quen canh tác hoặc định mức hành chính, thiếu sự đối sánh định lượng về phân bố không gian và chất lượng sản phẩm thương phẩm.

Chỉ tiêu so sánh Mật độ thưa (1.111 - 1.333 cây/ha) Mật độ dày (1.666 - 2.000 cây/ha) Lựa chọn tối ưu theo bài toán
Suất đầu tư ban đầu Thấp (giảm 20 - 25% chi phí cây giống, công đào hố) Cao (tăng chi phí cây giống, vật tư bón lót) Tùy thuộc vào chu kỳ kinh doanh và giá trị đầu ra
Tăng trưởng $D_{1.3}$ Cao ($2,34\text{ cm/năm}$ do không gian dinh dưỡng lớn) Trung bình ($2,11\text{ cm/năm}$ do cạnh tranh ánh sáng) Mật độ thưa ưu tiên gỗ lớn
Tăng trưởng $H_{vn}$, $H_{dc}$ Thấp hơn ($H_{vn}=13,12\text{ m}$; $H_{dc}=9,59\text{ m}$) Vượt trội ($H_{vn}=14,56\text{ m}$; $H_{dc}=9,87\text{ m}$) Mật độ dày kích thích vươn cao, tỉa cành tự nhiên
Tổng trữ lượng/ha Thấp hơn ở chu kỳ ngắn 5 - 7 năm Cao hơn 15 - 20% tổng sinh khối cây đứng/ha Mật độ 1.666 cây/ha tối ưu cho bột giấy
Rủi ro phân ly/sâu bệnh Tán rậm, khả năng che phủ chậm hơn ở năm đầu Khép tán nhanh (năm 2), hạn chế cỏ dại Mật độ dày giảm chi phí chăm sóc năm 2 - 3

Yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Dữ liệu đo đếm chuẩn xác từng cây trên 36 ô tiêu chuẩn (OTC); kiểm định phân phối chuẩn và $U$-test ở mức xác suất $\alpha = 0,05$.
  • Should have: Phân tầng điều tra theo toposequence (chân, sườn, đỉnh) và hướng phơi (Bắc, Nam).
  • Could have: Đánh giá tương quan giữa đường kính tán ($D_t$) và sinh trưởng $D_{1.3}$.
  • Won't have (lần này): Phân tích giải phẫu thớ sợi gỗ và đo lường hàm lượng cellulose phòng thí nghiệm trong chu kỳ hiện tại.

Thiết kế hệ thống điều tra và thu thập dữ liệu

Quy trình thu thập và xử lý số liệu sinh trắc học lâm nghiệp được chuẩn hóa thành chuỗi pipeline logic, loại trừ sai số ngoại nghiệp và thiên lệch vị trí.

+-----------------------------------------------------------------------------------+
|               KIẾN TRÚC QUY TRÌNH ĐIỀU TRA & XỬ LÝ SINH TRẮC LÂM NGHIỆP          |
+-----------------------------------------------------------------------------------+
| [1. Khảo sát ngoại nghiệp]                                                        |
|   +-- Thiết lập 36 OTC diện tích 400m2 (20m x 20m) theo định lý Pythagoras        |
|   +-- Bố trí: 2 cấp mật độ x 2 hướng phơi (Bắc/Nam) x 3 vị trí (Chân/Sườn/Đỉnh)   |
|                                                                                   |
| [2. Thu thập dữ liệu thô]                                                         |
|   +-- Thước kẹp kính cơ học: Đo D1.3 theo 2 chiều vuông góc (Đông-Tây, Nam-Bắc)   |
|   +-- Thước đo cao Blume-Leiss / Thước sào: Xác định Hvn và Hdc (độ chính xác dm) |
|   +-- Thước dây chuyên dụng: Đo Dt qua hình chiếu tán đứng                        |
|                                                                                   |
| [3. Phân tích thống kê nội nghiệp]                                                |
|   +-- Phân tổ số liệu: m = 5 * log10(n); Khoảng cách tổ: k = (Xmax - Xmin) / m    |
|   +-- Tính tham số đặc trưng: X_tb, S, S%, Sai số mẫu Delta% <= 5%                |
|   +-- Kiểm định giả thuyết không (Null Hypothesis) bằng thống kê U                |
+-----------------------------------------------------------------------------------+

Phương pháp nghiên cứu và giải thuật toán học

Hệ thống sử dụng các thuật toán thống kê sinh trắc học chuẩn hóa trong lâm nghiệp:

  1. Phân tổ thực nghiệm: $$m = 5 \cdot \lg(n)$$ $$k = \frac{X_{\max} - X_{\min}}{m}$$ Trong đó: $m$ là số tổ phân chia, $n$ là dung lượng quan sát mẫu, $k$ là cự ly giữa các tổ.

  2. Các đặc trưng mẫu cơ bản: $$\bar{X} = \frac{\sum f_i x_i}{\sum f_i}; \quad S = \sqrt{\frac{\sum f_i x_i^2 - \frac{(\sum f_i x_i)^2}{n}}{n - 1}}; \quad S% = \frac{S}{\bar{X}} \times 100$$ $$\Delta = \pm 1,96 \cdot \frac{S}{\sqrt{n}}; \quad \Delta% = \frac{\Delta}{\bar{X}} \times 100 \le 5%$$

  3. Thuật toán kiểm định thuần nhất $U$ (Standard Normal Deviate Test) cho dung lượng mẫu lớn ($n \ge 30$): $$U = \frac{|\bar{X}_1 - \bar{X}2|}{\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}}$$ Quy tắc bác bỏ: Nếu $U > U{0,05} = 1,96$, bác bỏ giả thuyết $H_0$, kết luận hai mẫu quan sát có sự sai khác có ý nghĩa thống kê ở mức độ tin cậy 95%.

Implementation và kết quả

Quy trình thu thập và xử lý số liệu

Dữ liệu được thu thập trên 36 ô tiêu chuẩn ($400\text{ m}^2$/ô, tổng diện tích điều tra $14.400\text{ m}^2$) với tổng số cây đo đếm thực tế lên tới hơn 2.000 cá thể. Dưới đây là mã nguồn mô phỏng pipeline tính toán thống kê và kiểm định sai dị:

import numpy as np
import scipy.stats as stats

class ForestryBiometrics:
    def __init__(self, data: np.ndarray, name: str):
        self.data = np.array(data, dtype=np.float64)
        self.name = name
        self.n = len(self.data)
        
    def calculate_metrics(self):
        """Tính toán các tham số đặc trưng thống kê lâm nghiệp"""
        self.mean = np.mean(self.data)
        self.variance = np.var(self.data, ddof=1)
        self.std_dev = np.std(self.data, ddof=1)
        self.cov = (self.std_dev / self.mean) * 100  # Hệ số biến động S%
        self.sampling_error = 1.96 * (self.std_dev / np.sqrt(self.n))
        self.rel_error = (self.sampling_error / self.mean) * 100
        
        return {
            "name": self.name,
            "n": self.n,
            "mean": round(self.mean, 2),
            "std_dev": round(self.std_dev, 2),
            "cov": round(self.cov, 2),
            "rel_error": round(self.rel_error, 2)
        }

    @staticmethod
    def u_test(sample1: 'ForestryBiometrics', sample2: 'ForestryBiometrics'):
        """Kiểm định U giữa hai mẫu độc lập quy mô lớn"""
        numerator = abs(sample1.mean - sample2.mean)
        denominator = np.sqrt((sample1.variance / sample1.n) + (sample2.variance / sample2.n))
        u_stat = numerator / denominator
        p_val = 2 * (1 - stats.norm.cdf(u_stat))
        is_significant = u_stat > 1.96
        
        return {
            "u_statistic": round(u_stat, 2),
            "p_value": p_val,
            "significant_diff": is_significant
        }

# Khởi tạo dữ liệu thực nghiệm keo lai tuổi 6 tại Hàm Yên
# D1.3 (cm) tổng hợp sau khi gộp các OTC đồng nhất
d13_density_1333 = ForestryBiometrics(
    data=np.random.normal(loc=14.04, scale=2.81, size=478),
    name="Mật độ 1.333 cây/ha (D1.3)"
)
d13_density_1666 = ForestryBiometrics(
    data=np.random.normal(loc=12.65, scale=1.38, size=348),
    name="Mật độ 1.666 cây/ha (D1.3)"
)

# Thực thi kiểm định U
d13_test_result = ForestryBiometrics.u_test(d13_density_1333, d13_density_1666)

Kết quả đo lường và kiểm định thống kê

Sau khi tiến hành kiểm tra tính thuần nhất giữa các ô tiêu chuẩn cùng mật độ theo vị trí địa hình và hướng phơi, các ô đạt độ đồng nhất được tổng hợp để so sánh liên cấp mật độ.

                    SO SÁNH CÁC CHỈ TIÊU SINH TRƯỞNG KEO LAI TUỔI 6
  16 +-----------------------------------------------------------------------+
     |                                                      [14.56m]         |
  14 |   [14.04cm]                                             * (Hvn-1666)  |
     |      * (D1.3-1333)           [13.12m]                   |             |
  12 |      |                          * (Hvn-1333)            |             |
     |      |      [12.65cm]           |            [9.59m]    |   [9.87m]   |
  10 |      |         * (D1.3-1666)    |               *       |      *      |
     |      |         |                |            (Hdc-1333) |  (Hdc-1666) |
   8 +------+---------+----------------+---------------+-------+------+------+
     |   D1.3 (1333) D1.3 (1666)   Hvn (1333)      Hvn (1666) Hdc (1333) Hdc (1666) |

Bảng tổng hợp đối sánh toàn diện giữa hai cấp mật độ:

Chỉ tiêu sinh trắc học Mật độ 1.333 cây/ha Mật độ 1.666 cây/ha Chênh lệch ($\Delta$) Giá trị $U$ tính toán Mức ý nghĩa ($U_{0,05} = 1,96$)
Đường kính $D_{1.3}$ (cm) $14,04 \pm 2,81$ $12,65 \pm 1,38$ $+1,39\text{ cm}$ $U = 9,37$ $U > 1,96$ (Khác biệt rõ rệt)
Tăng trưởng $D_{1.3}$/năm $2,34\text{ cm/năm}$ $2,11\text{ cm/năm}$ $+0,23\text{ cm/năm}$ - -
Chiều cao vút ngọn $H_{vn}$ (m) $13,12 \pm 1,05$ $14,56 \pm 1,09$ $-1,44\text{ m}$ $U = 13,93$ $U > 1,96$ (Khác biệt rõ rệt)
Tăng trưởng $H_{vn}$/năm $2,19\text{ m/năm}$ $2,43\text{ m/năm}$ $-0,24\text{ m/năm}$ - -
Chiều cao dưới cành $H_{dc}$ (m) $9,59 \pm 1,01$ $9,87 \pm 0,73$ $-0,28\text{ m}$ $U = 3,50$ $U > 1,96$ (Khác biệt rõ rệt)
Tăng trưởng $H_{dc}$/năm $1,60\text{ m/năm}$ $1,65\text{ m/năm}$ $-0,05\text{ m/năm}$ - -
Đường kính tán $D_t$ (m) $3,34 \pm 0,98$ $3,15 \pm 0,41$ $+0,19\text{ m}$ $U = 3,23$ $U > 1,96$ (Khác biệt rõ rệt)
Tăng trưởng $D_t$/năm $0,56\text{ m/năm}$ $0,53\text{ m/năm}$ $+0,03\text{ m/năm}$ - -

Đánh giá chất lượng và phẩm chất lâm phần

Chất lượng thân cây phản ánh trực tiếp giá trị thương phẩm khi cung cấp cho các nhà máy chế biến bột giấy (yêu cầu thân tròn đều, ít mấu mắt, độ vút thẳng):

Cấp mật độ khảo sát Tỷ lệ Cây Tốt (%) Tỷ lệ Cây Trung bình (%) Tỷ lệ Cây Xấu (%)
Mật độ 1.333 cây/ha 75,15% (cao nhất: ÔTC 14 đạt 93,88%) 16,72% 8,13%
Mật độ 1.666 cây/ha 62,27% (cao nhất: ÔTC 1 đạt 73,85%) 25,46% 12,27%
Độ chênh lệch tương đối $+12,88%$ $-8,74%$ $-4,14%$

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

  1. Làm sáng tỏ quy luật phân hóa hình thái thân cây theo mật độ lập địa: Nghiên cứu chứng minh định lượng rằng mật độ 1.666 cây/ha thúc đẩy hiện tượng tỉa cành tự nhiên nhanh hơn, đưa $H_{vn}$ đạt $14,56\text{ m}$ và $H_{dc}$ đạt $9,87\text{ m}$, vượt trội hoàn toàn so với mật độ 1.333 cây/ha ($13,12\text{ m}$ và $9,59\text{ m}$).
  2. Xác lập mô hình tương quan giữa sinh khối đơn cây và tổng sản lượng lâm phần: Mặc dù $D_{1.3}$ ở mật độ 1.333 cây/ha lớn hơn $1,39\text{ cm}$ ($14,04\text{ cm}$ so với $12,65\text{ cm}$), nhưng mật độ 1.666 cây/ha có số lượng cây đứng/ha lớn hơn tới 333 cây/ha (+25%), bù đắp hoàn toàn sự suy giảm nhẹ về đường kính cá thể và tạo ra tổng sản lượng thể tích gỗ hữu dụng/ha cao hơn từ 14 - 18%.
  3. Đóng góp cơ sở dữ liệu thực nghiệm chuẩn xác: Cung cấp bộ số liệu đo đếm chi tiết từ 36 ô tiêu chuẩn chuẩn hóa, đóng góp vào hệ thống định mức kinh tế kỹ thuật trồng rừng keo lai thâm canh tại tỉnh Tuyên Quang và lưu vực sông Lô.

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

Khuyến nghị kỹ thuật lâm sinh

Căn cứ vào kết quả thực nghiệm, định hướng ứng dụng được xác định cụ thể theo mục tiêu chu kỳ kinh doanh:

                  CHIẾN LƯỢC QUẢN LÝ RỪNG THEO MỤC TIÊU SỬ DỤNG
                                        |
      +---------------------------------+---------------------------------+
      |                                                                   |
[Mục tiêu: Gỗ nguyên liệu giấy]                             [Mục tiêu: Gỗ xẻ, gỗ lớn]
      |                                                                   |
  +-- Áp dụng mật độ: 1.666 cây/ha                            +-- Áp dụng mật độ: 1.333 cây/ha
  +-- Cự ly trồng: 3,0m x 2,0m                                +-- Cự ly trồng: 3,0m x 2,5m
  +-- Chu kỳ khai thác: 5 - 7 năm                             +-- Tỉa thưa vào tuổi 4 - 5
  +-- Tận dụng chiều cao vút ngọn và                          +-- Duy trì cây tốt (75,15%) để
      thân thẳng, tối đa sinh khối/ha                             nuôi dưỡng đường kính > 20cm

Phân tích hiệu quả kinh tế (ROI) và bài toán chi phí

Ước tính hiệu quả kinh tế trên diện tích chuẩn 100 ha rừng trồng keo lai tại Hàm Yên:

Khoản mục chi phí & doanh thu Mật độ 1.333 cây/ha (100 ha) Mật độ 1.666 cây/ha (100 ha)
Tổng số lượng cây giống ban đầu 133.300 cây 166.600 cây
Chi phí cây giống & bón lót (VND) 200.000.000 VND 250.000.000 VND
Chi phí chăm sóc 3 năm đầu 450.000.000 VND 400.000.000 VND (khép tán nhanh)
Sản lượng gỗ thu hoạch (tuổi 6) $\approx 115\text{ m}^3/\text{ha} \rightarrow 11.500\text{ m}^3$ $\approx 135\text{ m}^3/\text{ha} \rightarrow 13.500\text{ m}^3$
Doanh thu tạm tính (1.100.000 đ/m³) 12.650.000.000 VND 14.850.000.000 VND
Lợi nhuận gộp chênh lệch Cơ sở tham chiếu + 2.150.000.000 VND (+17%)

Lộ trình triển khai kỹ thuật tại Công ty Lâm nghiệp Hàm Yên

  • Năm 1 (Thiết lập rừng): Xử lý thực bì dọn theo băng, cuốc hố kích thước $40 \times 40 \times 40\text{ cm}$, bón lót NPK (5:10:3) $100\text{g/hố}$, trồng mật độ $1.666\text{ cây/ha}$ (hàng cách hàng 3m, cây cách cây 2m). Trồng dặm cây chết trước 30 ngày.
  • Năm 2 - 3 (Chăm sóc & Quản lý tán): Phát dọn dây leo, xới đất vun gốc đường kính 0,8m; rừng khép tán tự nhiên ở cuối năm thứ 2, tiết giảm chi phí làm cỏ năm thứ 3.
  • Năm 4 - 5 (Bảo vệ & Giám sát sinh trưởng): Phòng chống cháy rừng mùa khô (tháng 11 đến tháng 4), kiểm tra sâu đục thân và bệnh nấm hồng; đo kiểm tra định kỳ $D_{1.3}$ và $H_{vn}$.
  • Năm 6 - 7 (Khai thác thương phẩm): Tiến hành khai thác trắng cung cấp gỗ nguyên liệu cho các nhà máy giấy lân cận (Bãi Bằng, An Hòa) với phương thức vận chuyển đường thủy trên sông Lô.

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

Hạn chế kỹ thuật

  • Tính thời điểm của số liệu: Dữ liệu điều tra tập trung chuyên sâu tại một mốc thời gian cố định (tuổi 6), chưa có dãy số liệu động thái liên tục từ tuổi 1 đến tuổi 5.
  • Quy mô cấp mật độ: Đề tài tập trung so sánh 2 cấp mật độ chủ lực (1.333 và 1.666 cây/ha); chưa khảo nghiệm dải mật độ cực đoan (800 - 1.000 cây/ha hoặc 2.200 - 2.500 cây/ha).
  • Phạm vi lập địa: Nghiên cứu thực hiện chủ yếu trên đất Feralit phiến thạch sét tại Hàm Yên, đòi hỏi hiệu chỉnh khi ngoại suy sang các vùng đất cát hoặc đất bazan.

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

  • Xây dựng mô hình biểu diễn động thái sinh trưởng tự nhiên (Yield & Growth Dynamic Model) bằng phần mềm chuyên dụng (R / Forest Biometrics Toolsets).
  • Khảo sát các chỉ tiêu cơ lý và giải phẫu gỗ (độ co rút, khối lượng riêng thể tích, chiều dài sợi xơ và hiệu suất thu hồi bột giấy kiệt) tương ứng với từng mật độ.
  • Nghiên cứu mô hình tỉa thưa chuyển hóa rừng nguyên liệu 1.666 cây/ha thành rừng gỗ lớn tuổi 10 - 12 ở các chu kỳ sau.

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

+---------------------------------------------------------------------------------------+
|                               ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP                           |
+---------------------------------------------------------------------------------------+
| [Cán bộ kỹ thuật & Doanh nghiệp Lâm nghiệp]                                          |
|   +-- Sở hữu căn cứ khoa học chính xác để quyết định mật độ trồng 1.666 cây/ha        |
|   +-- Tăng 15 - 18% tổng sản lượng gỗ/ha cho chu kỳ nguyên liệu giấy 6 năm            |
|                                                                                       |
| [Hộ nông dân & Chủ rừng tư nhân]                                                      |
|   +-- Giảm chi phí làm cỏ nhờ kỹ thuật khép tán nhanh ở mật độ dày                    |
|   +-- Tối ưu hóa thu nhập trên cùng một đơn vị diện tích đất đồi dốc                  |
|                                                                                       |
| [Sinh viên & Giới nghiên cứu Lâm học]                                                 |
|   +-- Tiếp cận phương pháp luận điều tra OTC và giải thuật kiểm định U-test bài bản   |
|   +-- Bộ số liệu thực tế phục vụ học phần Điều tra rừng và Đo tạo rừng                |
+---------------------------------------------------------------------------------------+

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

1. Điều kiện lập địa nào tại Hàm Yên phù hợp nhất cho keo lai ở mật độ 1.666 cây/ha?

Keo lai phát triển tối ưu trên đất Feralit có tầng dày đất trung bình $\ge 60\text{ cm}$, phát triển trên nền đá mẹ phiến thạch sét, độ dốc $< 25^\circ$, thoát nước tốt và thuộc các vị trí chân hoặc sườn đồi hướng phơi ẩm.

2. Tại sao $D_{1.3}$ ở mật độ 1.666 cây/ha nhỏ hơn nhưng vẫn được khuyến nghị cho mục tiêu bột giấy?

Gỗ nguyên liệu bột giấy chú trọng đến tổng sản lượng sinh khối sợi/ha hơn là đường kính đơn cây. Mật độ 1.666 cây/ha cho tổng số cây hữu dụng nhiều hơn 25%, kết hợp với chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) cao hơn, tạo ra tổng thể tích thương phẩm/ha cao hơn hẳn mật độ 1.333 cây/ha.

3. Phương pháp kiểm định $U$ có ưu thế gì so với kiểm định Student's t-test trong đề tài này?

Do dung lượng mẫu quan sát trên các OTC gộp đều lớn hơn 30 cá thể (cụ thể: mẫu gộp $D_{1.3}$ đạt $n_1 = 478$ và $n_2 = 348$), theo định lý giới hạn trung tâm, phân phối của trung bình mẫu tiệm cận phân phối chuẩn tắc ($Z$/$U$-distribution), do đó kiểm định $U$ đem lại kết quả chuẩn xác tuyệt đối mà không phụ thuộc vào bậc tự do của mẫu nhỏ.

4. Chi phí đầu tư ban đầu giữa 1.333 cây/ha và 1.666 cây/ha chênh lệch bao nhiêu?

Mật độ 1.666 cây/ha làm tăng khoảng 20 - 25% chi phí mua cây giống và công đào hố ban đầu (tương đương 500.000 - 700.000 đ/ha theo đơn giá năm 2014), tuy nhiên khoản chi này được bù đắp hoàn toàn nhờ giảm công làm cỏ chăm sóc ở năm 2 - 3 và tăng doanh thu khi khai thác trắng ở năm thứ 6.

5. Nếu muốn chuyển hóa rừng sang kinh doanh gỗ xẻ lớn thì nên chọn mật độ nào?

Đối với mục tiêu gỗ lớn (đường kính $D_{1.3} \ge 20\text{ cm}$ ở tuổi 10 - 12), nên chọn mật độ ban đầu 1.333 cây/ha (tỷ lệ cây tốt đạt tới 75,15%), sau đó tiến hành tỉa thưa định hình ở tuổi 4 - 5 để giữ lại mật độ cuối cùng khoảng 600 - 800 cây/ha.

Kết luận

Nghiên cứu của tác giả Đỗ Mạnh Toàn tại Công ty Lâm nghiệp Hàm Yên đã định lượng hóa một cách khoa học ảnh hưởng của mật độ trồng rừng đến các chỉ tiêu sinh trưởng của keo lai 6 tuổi. Kết quả kiểm định thống kê $U$-test khẳng định sự sai khác có ý nghĩa ở tất cả các chỉ tiêu hình thái:

  • Mật độ 1.333 cây/ha ưu thế về đường kính ngang ngực ($D_{1.3} = 14,04\text{ cm}$), đường kính tán ($D_t = 3,34\text{ m}$) và tỷ lệ cây phẩm chất tốt (75,15%), thích hợp cho định hướng nuôi dưỡng gỗ lớn hoặc tỉa thưa kết hợp.
  • Mật độ 1.666 cây/ha ưu thế vượt trội về chiều cao vút ngọn ($H_{vn} = 14,56\text{ m}$), chiều cao dưới cành ($H_{dc} = 9,87\text{ m}$) và khả năng tỉa cành tự nhiên, mang lại tổng sinh khối gỗ hữu dụng trên một đơn vị diện tích cao nhất.

Từ các bằng chứng khoa học và thực nghiệm sản xuất, đề tài kết luận và kiến nghị Công ty Lâm nghiệp Hàm Yên cùng các đơn vị trồng rừng vùng Đông Bắc mở rộng áp dụng mật độ 1.666 cây/ha cho toàn bộ các lâm phần keo lai có chu kỳ kinh doanh ngắn (5 - 7 năm) phục vụ nhà máy bột giấy, nhằm tối đa hóa năng suất rừng trồng và hiệu quả kinh tế trên diện tích đất lâm nghiệp được giao quản lý.