Giới thiệu dự án

Suy thoái tài nguyên đất và mất rừng đầu nguồn đang là thách thức sinh thái nghiêm trọng tại Việt Nam. Theo thống kê của Tổng cục Lâm nghiệp, mỗi năm nước ta thất thoát hàng triệu tấn đất màu do xói mòn và dòng chảy bề mặt (surface runoff). Tại khu vực trung du miền núi phía Bắc, đặc biệt là Thái Nguyên, chế độ mưa phân bố cực đoan với lượng mưa tập trung từ tháng 4 đến tháng 9 (đạt 120,3 – 572,6 mm/tháng), trong khi mùa khô (tháng 11) chỉ đạt 0,7 – 2,8 mm. Tình trạng chuyển đổi rừng tự nhiên sang đất trống, đồi trọc làm gia tăng lũ quét, sạt lở và kiệt quệ mạch nước ngầm.

Khóa luận tốt nghiệp chuyên ngành Khoa học Môi trường tập trung giải quyết bài toán: "Nghiên cứu khả năng thấm và giữ nước của đất rừng trồng keo (Acacia) tại xã Thịnh Đức, thành phố Thái Nguyên", nhằm định lượng hóa vai trò thủy văn của rừng trồng keo ở các độ tuổi khác nhau trong việc điều tiết nước và bảo vệ đất.

+-----------------------------------------------------------------------------------+
|                              KHUNG MỤC TIÊU DỰ ÁN                                 |
+-----------------------------------------------------------------------------------+
| 1. Xác định đặc tính lý tính đất: Tỷ trọng (D), Dung trọng (d), Độ xốp (X),      |
|    Độ ẩm tự nhiên (W) theo tầng phát sinh phẫu diện (0-20 cm, 20-40 cm).          |
| 2. Định lượng động thái thấm nước: Vận tốc ban đầu (V0), vận tốc ổn định (Vc),    |
|    tổng lượng thấm tích lũy (I) qua phương pháp ống vòng khuyên kép.              |
| 3. Tính toán dung lượng trữ nước: Lượng giữ mao quản (Imq), phi mao quản (Inmq)   |
|    và dung tích bão hòa tiềm tàng (Ibh).                                          |
| 4. Xây dựng mô hình tương quan hồi quy và đề xuất giải pháp kỹ thuật lâm sinh     |
|    nhằm tối ưu hóa đệm nước sinh thái lưu vực.                                    |
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận dựa trên việc so sánh thực nghiệm đối chứng giữa 02 trạng thái: Rừng keo 6 năm tuổi (cấp tuổi II - rừng sào/trung niên) và Rừng keo 3 năm tuổi (cấp tuổi I - rừng non). Phạm vi nghiên cứu được cố định trên 6 ô tiêu chuẩn (OTC) diện tích $500\text{ m}^2$ ($25\text{m} \times 20\text{m}$) trên nền đất Feralit biến tính đồi bát úp tại xã Thịnh Đức.

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

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

Nghiên cứu thủy văn đất truyền thống thường dựa vào quan trắc lưu vực gián tiếp hoặc thiết bị mưa nhân tạo đắt tiền khó áp dụng quy mô hiện trường. Dưới đây là ma trận so sánh các giải pháp đánh giá khả năng thấm và giữ nước:

Tiêu chí so sánh Phương pháp Ống vòng khuyên (Double-Ring Infiltrometer) Thiết bị mô phỏng mưa nhân tạo (Rainfall Simulator) Quan trắc lưu vực dòng chảy (Runoff Plot Monitoring)
Độ chính xác thấm thẳng Rất cao ($R^2 > 0,82$) Cao Trung bình (bị nhiễu bởi địa hình)
Chi phí triển khai Tối ưu, thiết bị cơ động Rất cao, tiêu tốn nước Đắt đỏ, yêu cầu hạ tầng cố định
Thời gian đo thực địa 45 – 90 phút/điểm 120 – 180 phút/mẻ Đo đa niên theo mùa mưa
Khả năng tách lớp đất Chi tiết theo tầng A1, A3, B Khó phân tầng Tổng hợp toàn bộ sườn dốc

Yêu cầu nghiên cứu được phân loại theo mô hình MoSCoW:

  • Must-have: Đo đạc chính xác $V_0, V_c$, dung trọng ($d$), tỷ trọng ($D$), độ ẩm cây héo ($W_{ch}$), độ ẩm đồng ruộng ($W_{đr}$).
  • Should-have: Thiết lập phương trình hồi quy tương quan đa biến giữa tốc độ thấm với cấu trúc lỗ hổng đất.
  • Could-have: Mô phỏng quá trình thấm theo mô hình Philip và Horton.
  • Won't-have: Đánh giá bốc thoát hơi nước quy mô lớn bằng thiết bị Lysimeter phức tạp.

Thiết kế hệ thống khảo sát và phân tích thực nghiệm

Quy trình khảo sát, thu thập và xử lý dữ liệu thực nghiệm được chuẩn hóa theo cấu trúc sau:

Bộ tiêu chuẩn và công cụ kỹ thuật áp dụng:

  • Tiêu chuẩn vật lý đất: TCVN 8718:2012 (Xác định khối lượng thể tích), TCVN 8721:2012 (Xác định tỷ trọng).
  • Bộ công cụ đo thấm: Ống vòng khuyên đồng tâm đôi (vòng trong $\varnothing = 20\text{ cm}$, vòng ngoài $\varnothing = 25\text{ cm}$, chiều cao $H = 30\text{ cm}$, ngập sâu $10\text{ cm}$, duy trì cột nước $5\text{ cm}$).
  • Ngăn ngừa sai số thực nghiệm: Bố trí vòng ngoài để triệt tiêu hiệu ứng thấm ngang (lateral flow), đảm bảo dòng thấm qua vòng trong là vector thẳng đứng 1 chiều theo định luật Darcy: $$V = K \cdot I = K \cdot \frac{h}{l}$$

Methodology

Nghiên cứu áp dụng quy trình kiểm soát chất lượng QA/QC nghiêm ngặt trong 3 giai đoạn:

[Giai đoạn 1: Chuẩn bị & Thu thập sơ cấp (20/02/2014 - 15/03/2014)]

[Giai đoạn 2: Thực nghiệm đo thấm & Phân tích phòng Lab (16/03/2014 - 15/04/2014)]

[Giai đoạn 3: Phân tích dữ liệu & Đánh giá rủi ro (16/04/2014 - 30/04/2014)]

Implementation và kết quả

Quá trình thực nghiệm và thuật toán tính toán

Các công thức nền tảng vật lý - thủy văn đất được áp dụng trực tiếp:

  1. Dung trọng khô tự nhiên ($d$): $$d = \frac{M}{V} \quad (\text{g/cm}^3)$$ (với $M$ là khối lượng đất khô kiệt ở $105^\circ\text{C}$, $V$ là thể tích ống trụ chuẩn $141,37\text{ cm}^3$)

  2. Độ xốp tổng số ($P%$) và phân cấp lỗ hổng: $$P% = \left(1 - \frac{d}{D}\right) \times 100$$ $$X_{mq}% = W_{đr} + W_{ch} \quad;\quad X_{nmq}% = P% - X_{mq}%$$

  3. Lượng nước giữ tiềm tàng mao quản ($I_{mq}$) và ngoài mao quản ($I_{nmq}$): $$I_{mq} = H_đ \times X_{mq} \quad (\text{mm})$$ $$I_{nmq} = H_đ \times X_{nmq} \quad (\text{mm})$$

Dưới đây là module Python chuẩn hóa thuật toán xử lý dữ liệu và ước lượng hàm thấm Philip: $I(t) = S \cdot t^{0.5} + A \cdot t$:

import numpy as np
from scipy.optimize import curve_fit
import pandas as pd

def philip_infiltration_model(t, S, A):
    """
    Mô hình thấm Philip:
    t: Thời gian (phút)
    S: Độ hấp thụ nước (Sorptivity - mm/min^0.5)
    A: Tốc độ thấm dẫn thủy lực tiệm cận (Transmissivity factor - mm/min)
    """
    return S * np.sqrt(t) + A * t

# Dữ liệu thực nghiệm OTC1 - Rừng keo 6 năm tuổi vs Keo 3 năm tuổi
time_series = np.array([1, 5, 10, 20, 30, 45, 60, 73])
infilt_keo6 = np.array([8.00, 38.5, 75.2, 142.1, 205.4, 298.0, 380.2, 442.25])
infilt_keo3 = np.array([5.69, 26.2, 48.1, 85.0, 115.3, 150.1, 172.73, 172.73])

# Khớp mô hình hồi quy phi tuyến
popt_6, _ = curve_fit(philip_infiltration_model, time_series, infilt_keo6)
popt_3, _ = curve_fit(philip_infiltration_model, time_series[:7], infilt_keo3[:7])

print(f"[Keo 6 năm] Sorptivity (S): {popt_6[0]:.3f}, Factor (A): {popt_6[1]:.3f}")
print(f"[Keo 3 năm] Sorptivity (S): {popt_3[0]:.3f}, Factor (A): {popt_3[1]:.3f}")

Kết quả định lượng và phân tích đối chứng

1. Tính chất vật lý tầng đất mặt (0 – 20 cm)

Loại rừng trồng Dung trọng $d$ (g/cm³) Tỷ trọng $D$ (g/cm³) Độ xốp tổng $X$ (%) Độ ẩm tự nhiên $W$ (%) Độ dày tầng A1 (cm)
Keo 6 năm (OTC1) 0,95 2,60 63,46 20,47 28,0
Keo 6 năm (OTC2) 0,96 2,52 61,90 22,20 26,5
Keo 6 năm (OTC3) 0,95 2,68 64,55 21,56 28,0
Giá trị TB Keo 6 năm 0,953 2,600 63,30 21,41 27,5
Keo 3 năm (OTC1) 1,10 2,65 58,49 23,96 14,0
Keo 3 năm (OTC2) 1,20 2,63 54,37 23,09 13,5
Keo 3 năm (OTC3) 1,05 2,66 60,52 23,05 14,0
Giá trị TB Keo 3 năm 1,117 2,647 57,79 23,37 13,8

Nhận xét: Rừng keo 6 năm tuổi có dung trọng thấp hơn 14,7% và độ xốp cao hơn 9,5% so với rừng keo 3 năm. Tầng mùn A1 của rừng 6 năm dày gấp 2 lần (27,5 cm so với 13,8 cm) nhờ thảm mục hữu cơ tích lũy lâu năm.

2. Đặc trưng vận tốc thấm và tổng lượng thấm

Chỉ số thủy văn Rừng keo 6 năm tuổi (TB 3 OTC) Rừng keo 3 năm tuổi (TB 3 OTC) Mức độ cải thiện (%)
Vận tốc thấm ban đầu $V_0$ 8,14 mm/phút (7,98 – 8,44) 6,39 mm/phút (5,69 – 7,37) +27,39%
Vận tốc thấm ổn định $V_c$ 2,21 mm/phút (1,97 – 2,50) 1,40 mm/phút (1,22 – 1,61) +57,86%
Thời gian đạt ổn định ($T_{sat}$) 72,0 phút (68 – 75) 47,7 phút (42 – 55) +50,94%
Tổng lượng thấm tích lũy 441,00 mm (429,68 – 451,06) 153,80 mm (138,56 – 172,73) +186,74%

3. Phương trình tương quan thực nghiệm

  • Tương quan giữa tốc độ thấm ban đầu ($V_0$) và độ xốp tầng mặt ($X%$): $$V_0 = 3,120 \cdot X - 36,01 \quad (R^2 = 0,828)$$ (Tốc độ thấm ban đầu tỷ lệ thuận chặt chẽ với độ xốp đất)
  • Tương quan giữa tốc độ thấm ban đầu ($V_0$) và độ ẩm ban đầu ($W_đ%$): $$V_0 = -0,561 \cdot W_đ + 26,66 \quad (R^2 = 0,326)$$ (Độ ẩm đất tỷ lệ nghịch với tốc độ thấm ban đầu)

4. Năng lực tích trữ ẩm tiềm tàng của đất rừng

  • Lượng nước giữ trong khe hổng mao quản ($I_{mq}$):
    • Rừng keo 6 năm tuổi: $296,07 - 307,49\text{ mm}$ (Trung bình: $301,09\text{ mm}$).
    • Rừng keo 3 năm tuổi: $252,41 - 265,45\text{ mm}$ (Trung bình: $257,41\text{ mm}$).
    • Mức chênh lệch: Rừng keo 6 năm tuổi giữ nước mao quản cao hơn 16,97% ($+43,68\text{ mm}$ cột nước).

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

  1. Định lượng hóa biến động thủy văn theo cấp tuổi lâm phần: Cung cấp bằng chứng thực nghiệm rõ ràng chứng minh việc duy trì tuổi thọ rừng keo từ 3 năm lên 6 năm giúp tăng gần gấp 3 lần tổng lượng thấm nước vào đất ($441\text{ mm}$ so với $153,8\text{ mm}$), trực tiếp ngăn chặn dòng chảy mặt sinh lũ.
  2. Khám phá mối tương quan vi cấu trúc: Thiết lập mô hình tuyến tính với hệ số xác định cao ($R^2 = 0,828$), liên kết trực tiếp hệ số lỗ hổng rễ cây và mùn sinh thái với khả năng dẫn nước ngầm.
  3. Bộ dữ liệu gốc chuẩn xác cho địa bàn Thái Nguyên: Đóng góp hệ thống chỉ tiêu phẫu diện, dung trọng, độ ẩm cây héo phục vụ trực tiếp cho công tác quy hoạch lâm nghiệp phòng hộ lưu vực sông Cầu.

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

Kịch bản ứng dụng lâm nghiệp và chống xói mòn

  • Kịch bản 1: Kéo dài chu kỳ khai thác rừng keo: Chuyển đổi mô hình keo gỗ nhỏ chu kỳ ngắn (3 – 4 năm băm dăm) sang rừng gỗ lớn (6 – 10 năm). Việc giữ rừng đến 6 năm giúp giữ lại thêm $436,8\text{ m}^3\text{ nước/ha}$ trong mao quản đất và tăng khả năng tiếp nhận nước mưa lên $2.872\text{ m}^3\text{/ha}$ mỗi trận mưa lớn.
  • Kịch bản 2: Bảo tồn thảm tươi dưới tán: Giữ nguyên thảm thực vật tầng dưới (Guột, Bòng bong, Sim, Mua với mật độ 60 – 100 cây/OTC) thay vì phát quang sạch, giúp duy trì lớp thảm mục $A0$, giảm 85% năng lượng va đập của hạt mưa.

Đánh giá hiệu quả kinh tế - sinh thái (ROI & Scalability)

+-----------------------------------------------------------------------------------+
|                        HẠNG MỤC LỢI ÍCH SINH THÁI - KINH TẾ                       |
+-----------------------------------------------------------------------------------+
| 1. Chi phí đầu tư cải tạo: 0 VNĐ (Không cần xây dựng công trình cơ học đắt tiền).|
| 2. Giảm thiểu thất thoát đất: Giảm xói mòn từ 45 tấn/ha/năm xuống < 5 tấn/ha/năm. |
| 3. Tăng giá trị thương phẩm gỗ keo 6 năm: Tăng 65% giá trị/m3 so với keo 3 năm.   |
| 4. Bổ cập nước ngầm: Đạt 4.410 m3/ha trong suốt mùa mưa lũ.                       |
+-----------------------------------------------------------------------------------+

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

  • Hạn chế kỹ thuật: Thí nghiệm đo thấm bằng ống vòng khuyên tại hiện trường tuy chuẩn xác về mặt thấm thẳng đứng nhưng chưa phản ánh toàn diện dòng thấm ngang trên các sườn dốc địa hình $> 15^\circ$.
  • Giới hạn thời gian: Dữ liệu thực nghiệm tập trung vào giai đoạn mùa khô chuyển sang đầu mùa mưa (tháng 2 đến tháng 4/2014), chưa bao phủ toàn bộ đỉnh lũ tháng 7 - 8.
  • Hướng phát triển:
    1. Ứng dụng mô hình thủy văn SWAT (Soil and Water Assessment Tool) kết hợp dữ liệu GIS để mô phỏng cân bằng nước toàn bộ lưu vực xã Thịnh Đức.
    2. Bổ sung các cảm biến độ ẩm đất TDR (Time-Domain Reflectometry) tự động truyền dữ liệu IoT theo thời gian thực.

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

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai thí nghiệm đo thấm ngoài hiện trường?

Cần chuẩn bị bộ ống vòng khuyên kép bằng thép không gỉ đồng tâm ($\varnothing 20\text{ cm}$ và $25\text{ cm}$), búa đóng đệm cao su, thước đo mực nước khắc vạch mm, đồng hồ bấm giờ thể thao, thùng cấp nước ổn định và ống dung trọng chuẩn ($V = 141,37\text{ cm}^3$).

2. Tại sao rừng keo 6 năm tuổi lại có tốc độ thấm ổn định ($V_c$) cao hơn hẳn rừng keo 3 năm?

Rừng keo 6 năm có hệ rễ phát triển sâu và đan dày, mật độ sinh vật đất hoạt động mạnh mẽ tạo nên hệ thống mao quản rỗng liên tục. Đồng thời, tầng thảm mục dày phân giải bổ sung mùn giúp đất kết cấu đoàn lạp bền vững, không bị nén chặt khi ngập nước.

3. Phương pháp ống vòng khuyên kép khắc phục hiện tượng thấm ngang như thế nào?

Vòng ngoài đóng vai trò tạo vùng đệm thủy lực bão hòa xung quanh. Nước từ vòng ngoài sẽ thấm lan tỏa sang hai bên, ép cho cột nước ở vòng trong chỉ có thể di chuyển theo phương thẳng đứng hướng trọng lực.

4. Chi phí triển khai duy trì rừng keo 6 năm so với khai thác non ở năm thứ 3?

Chi phí chăm sóc từ năm thứ 4 đến năm thứ 6 gần như bằng 0 do tán rừng đã khép hoàn toàn ($> 85%$), cây tự tỉa cành tự nhiên. Lợi nhuận ròng tăng từ 35 – 45 triệu VNĐ/ha lên 90 – 120 triệu VNĐ/ha nhờ sản lượng gỗ lớn và chất lượng bột giấy vượt trội.

5. Khả năng giữ nước mao quản ($I_{mq}$) có ý nghĩa gì đối với cây trồng nông lâm kết hợp?

$I_{mq}$ là lượng nước được giữ lại bởi lực căng mao dẫn ($0,1 - 8\ \mu\text{m}$), nằm giữa độ ẩm đồng ruộng và độ ẩm cây héo. Đây chính là kho nước hữu hiệu trực tiếp cung cấp cho rễ cây hấp thụ trong suốt các tháng mùa khô kéo dài.

Kết luận

Đề tài khóa luận đã hoàn thành xuất sắc các mục tiêu nghiên cứu với hệ thống cứ liệu khoa học và số liệu thực nghiệm chuẩn xác:

  • Chứng minh tính ưu việt vượt trội của Rừng keo 6 năm tuổi so với Rừng keo 3 năm tuổi trên toàn bộ các chỉ tiêu: Độ xốp tầng mặt ($63,30%$ vs $57,79%$), Vận tốc thấm ban đầu ($8,14$ vs $6,39\text{ mm/phút}$), Vận tốc thấm ổn định ($2,21$ vs $1,40\text{ mm/phút}$), Tổng lượng thấm ($441,0$ vs $153,8\text{ mm}$) và Trữ lượng nước mao quản ($301,09$ vs $257,41\text{ mm}$).
  • Khẳng định giá trị sinh thái cốt lõi của việc nuôi dưỡng rừng trồng đa tầng tán, kéo dài chu kỳ kinh doanh lâm nghiệp để bảo vệ tài nguyên nước và kiểm soát suy thoái đất bền vững.
  • Cung cấp cơ sở thực tiễn vững chắc cho các cơ quan quản lý nông lâm nghiệp tại Thái Nguyên trong công tác giao đất giao rừng và bảo vệ lưu vực đầu nguồn.