Giới thiệu dự án

Cháy rừng là một trong những thảm họa sinh thái nghiêm trọng nhất trên toàn cầu, đe dọa trực tiếp đến an ninh môi trường và đa dạng sinh học. Theo số liệu thống kê từ Tổ chức Bảo tồn Thiên nhiên Quốc tế (IUCN), Chương trình Phát triển Liên Hợp Quốc (UNDP) và Quỹ Quốc tế Bảo vệ Thiên nhiên (WWF), thế giới mất trung bình khoảng 15 triệu hecta rừng mỗi năm, trong đó diện tích thiệt hại do hỏa hoạn chiếm tới 23%. Tại Việt Nam, theo dữ liệu từ Cục Kiểm lâm, giai đoạn 2005 – 2008 đã ghi nhận 2.772 vụ cháy rừng, thiêu rụi 15.675 ha rừng (trong đó có 4.059 ha rừng tự nhiên và 11.617 ha rừng trồng), gây tổn thất trung bình 3.918 ha rừng/năm. Rừng thông (Pinus) là đối tượng đặc biệt nhạy cảm với hỏa hoạn do chứa hàm lượng nhựa cao, thảm mục giàu cutin khó phân hủy và tích tụ khối lượng vật liệu cháy lớn.

+-------------------------------------------------------------------------+
|                           TAM GIÁC LỬA RỪNG                             |
|                                                                         |
|                         [ Nguồn nhiệt / Mồi lửa ]                        |
|                                    / \                                  |
|                                   /   \                                 |
|                                  /     \                                |
|                                 /       \                               |
|       [ Ôxy khí quyển: ~21% ] --------- [ Vật liệu cháy: Thảm lá thông ] |
|                                                                         |
|    Yếu tố địa hình (Độ dốc θ): Khuếch đại đối lưu nhiệt & Tốc độ cháy    |
+-------------------------------------------------------------------------+

Vấn đề cốt lõi đặt ra là hầu hết các phương pháp cảnh báo cháy rừng tại Việt Nam hiện nay chủ yếu dựa vào các chỉ số khí tượng vĩ mô đơn biến (như chỉ tiêu tổng hợp P.G. Nesterov hoặc số ngày khô hạn liên tục), bỏ qua yếu tố địa hình vi mô—đặc biệt là độ dốc sườn đồi. Thực tế cho thấy, sườn dốc làm biến đổi hoàn toàn cơ chế truyền nhiệt bức xạ và đối lưu, khiến ngọn lửa áp sát bề mặt thực bì phía trước và gia tăng vận tốc lan tràn theo cấp số nhân. Thiếu hụt các mô hình định lượng liên hệ giữa độ dốc, tải lượng nhiên liệu và hành vi đám cháy dẫn đến việc hoạch định khoảng cách an toàn tiếp cận và phân bổ nhân lực chữa cháy thiếu chính xác, gây rủi ro cao cho lực lượng cứu hỏa.

Đề tài khóa luận "Nghiên cứu ảnh hưởng của độ dốc đến tốc độ lan tràn của đám cháy rừng" được thực hiện tại Trung tâm Nghiên cứu và Phát triển Lâm nghiệp Sóc Sơn (Hà Nội) nhằm giải quyết triệt để bài toán định lượng trên.

Mục tiêu cụ thể của dự án:

  1. Định lượng đặc trưng cấu trúc vật liệu cháy dưới tán rừng thông (Pinus merkusii và Pinus kesiya), bao gồm tỷ lệ thành phần (lá, cành, vỏ, quả), kích thước hình học và khối lượng tích tụ thực tế ($M_{vl}$).
  2. Thiết kế và vận hành hệ số thử nghiệm mô phỏng cháy thực địa ($8\text{ m}^2$) với dải độ dốc kiểm soát từ $5^\circ$ đến $35^\circ$ kết hợp các mức tải lượng vật liệu cháy từ $5$ đến $20\text{ tấn/ha}$.
  3. Xác lập phương trình hồi quy đa biến biểu diễn mối quan hệ giữa tốc độ lan tràn ($V_c$), chiều cao ngọn lửa ($H_l$), cường độ đám cháy ($I_c$) với các biến số độ dốc ($\theta$), độ ẩm nhiên liệu ($W_{vl}$), tốc độ gió ($V_g$), nhiệt độ ($T_{kk}$) và độ ẩm không khí ($W_{kk}$).
  4. Xây dựng bảng quy chuẩn kỹ thuật phòng cháy chữa cháy rừng (PCCCR): định mức khoảng cách tiếp cận an toàn, xác định chu vi đám cháy theo thời gian và định mức nhân lực tối đa ứng cứu theo từng cấp độ dốc.

Phạm vi và giới hạn nghiên cứu: Nghiên cứu tập trung vào đám cháy bề mặt (cháy mặt đất) ở giai đoạn khởi phát dưới tán rừng thông tại khu vực Sóc Sơn - Hà Nội; không xét đến hiện tượng cháy tán tầng cao hoặc cháy ngầm dưới tầng than bùn sâu.

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

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

Trước khi nghiên cứu thực nghiệm được triển khai, công tác dự báo và điều hành chữa cháy rừng tại Việt Nam dựa trên các công cụ và mô hình tham chiếu truyền thống với những ưu nhược điểm rõ rệt:

Mô hình / Phương pháp Cơ chế tính toán Ưu điểm Nhược điểm & Lỗ hổng kỹ thuật
Chỉ số Nesterov ($K$) $K = \sum (T \cdot d)$, với $T$ là nhiệt độ, $d$ là độ hụt điểm sương Dễ tính toán từ dữ liệu trạm khí tượng cơ bản Hoàn toàn bỏ qua độ dốc địa hình, gió vi địa hình và tải lượng vật liệu cháy
Chỉ số ngày khô hạn (Phạm Ngọc Hưng) Đếm số ngày không mưa hoặc lượng mưa $< 5\text{ mm}$ liên tục Triển khai nhanh trên diện rộng cấp vùng Độ chính xác cục bộ thấp, không phản ánh được tốc độ lan tràn khi cháy xảy ra
Mô hình thực nghiệm Traubaud (Pháp) $V_p = 0.83 \cdot f(U_v, H, L_t, T_e)$ Có tính đến chiều cao thực bì và gió Áp dụng cho thảm thực vật Địa Trung Hải, thiếu biến số độ dốc $\theta$
Mô hình Butler & Catchpole (Mỹ, 2007) Buồng đốt nhiên liệu gỗ mùn (Populus tremuloides) Định lượng rõ bước nhảy tốc độ ở dốc $> 25^\circ$ Điều kiện phòng thí nghiệm kín, không tương thích hoàn toàn với thảm lá thông tự nhiên

Ma trận yêu cầu hệ thống (Phân loại MoSCoW):

  • Must have: Khung giàn đốt thực nghiệm thay đổi được góc nghiêng chính xác từ $5^\circ - 35^\circ$; hệ thống đo vi trắc ảnh kết hợp đồng hồ bấm giây kỹ thuật số; mô hình toán tương quan hồi quy có hệ số $R^2$ kiểm định và mức ý nghĩa $p < 0.05$.
  • Should have: Phân tích cấu trúc thành phần thảm khô theo kích thước hạt nhiên liệu; máy đo khí tượng đa chỉ tiêu tức thời tại hiện trường.
  • Could have: Phần mềm tự động hóa tính toán chu vi đám cháy và định mức nhân công dựa trên bảng tra cứu số hóa.
  • Won't have (trong phạm vi đề tài): Mô phỏng 3D động lực học chất lưu khí cháy (CFD - Fire Dynamics Simulator).

Thiết kế hệ thống

Hệ thống nghiên cứu thực nghiệm và phân tích dữ liệu được chuẩn hóa thành mô hình tích hợp bao gồm: Thiết bị thực nghiệm hiện trường, Thiết bị thu thập dữ liệu quang trắc - khí tượng, và Khung xử lý toán tin học.

+--------------------------------------------------------------------------------+
|                   KIẾN TRÚC HỆ THỐNG ĐO ĐẠC VÀ XỬ LÝ DỮ LIỆU                  |
+--------------------------------------------------------------------------------+
|  [HỆ THỐNG ĐỐT THỰC NGHIỆM]                                                    |
|  - Khung sàn tôn lá chịu nhiệt: 2.0m x 4.0m (Diện tích: 8.0 m2)               |
|  - Lưới thép chống trượt thảm lá: Mắt lưới 2.0 cm                              |
|  - Trục nâng cơ khí điều chỉnh độ dốc: h = 4.0 * sin(theta) (5 đến 35 độ)     |
|  - Tải lượng vật liệu cháy kiểm soát: 4kg, 8kg, 12kg, 16kg (5-20 tấn/ha)       |
+--------------------------------------------------------------------------------+
                                       |
                                       v
+--------------------------------------------------------------------------------+
|  [THU THẬP THÔNG SỐ ĐỒNG BỘ]                                                   |
|  - Cảm biến khí tượng đa năng: Tkk (oC), Wkk (%), Vg (m/s)                     |
|  - Cân phân tích & sấy hồng ngoại: Xác định độ ẩm Wa (%), W0 (%)               |
|  - Hệ thống ghi hình trắc địa: Camera 4 pha chụp kèm thước vạch chuẩn nhôm     |
|  - Đồng hồ bấm giây điện tử: Độ chính xác Delta_t = 0.01 s                     |
+--------------------------------------------------------------------------------+
                                       |
                                       v
+--------------------------------------------------------------------------------+
|  [PHÂN TÍCH VÀ MÔ HÌNH HÓA TOÁN HỌC]                                           |
|  - Thuật toán trắc vi ảnh (Photogrammetry processing) xác định Hlua, Vc        |
|  - Khử nhiễu và chuyển đổi biến phi tuyến: Tuyến tính, Lũy thừa, Logarit       |
|  - Phân tích hồi quy đa biến (Multivariate Regression ANOVA)                   |
|  - Công thức cường độ cháy Byram: Ic = H * W * R                               |
+--------------------------------------------------------------------------------+

Thông số kỹ thuật của trang thiết bị nghiên cứu:

  • Sàn đốt thực nghiệm: Khung thép hộp $30 \times 60\text{ mm}$, mặt sàn tôn mạ kẽm dày $1.2\text{ mm}$, kích thước danh định $2.0\text{ m} \times 4.0\text{ m}$ (diện tích hữu dụng $8.0\text{ m}^2$). Bề mặt trải lưới thép bọc nhựa mắt $2.0\text{ cm}$ cố định nhằm chống hiện tượng trượt cơ học của thảm lá ở các góc dốc cao ($> 20^\circ$).
  • Cơ cấu chỉnh dốc: Trụ nâng dạng trục vít có chia vạch định vị, tính toán độ cao nâng đầu theo hàm lượng giác: $$h = L \cdot \sin(\theta) = 4.0 \cdot \sin(\theta)\text{ (m)}$$
  • Thiết bị đo khí tượng: Máy đo thời tiết cầm tay Kestrel 4500 (nhiệt độ $\pm 0.5^\circ\text{C}$, độ ẩm $\pm 3%$, tốc độ gió $\pm 0.1\text{ m/s}$).
  • Trắc vi quang học: Máy ảnh kỹ thuật số độ phân giải cao kết hợp hệ thống 2 thước tiêu chuẩn nhôm dựng song song ở 2 bên mép sàn đốt.

Methodology

Quy trình nghiên cứu áp dụng phương pháp luận kết hợp giữa Điều tra sinh thái thực địa, Thực nghiệm vật lý có kiểm soát và Mô hình hóa thống kê đa biến (Multivariate Statistical Modeling):

[Khảo sát 3 ô tiêu chuẩn] -> [Thu thập thảm khô] -> [Sấy xác định ẩm độ]
                                                           |
[Chế áp sàn đốt 8m2] <--- [Hiệu chỉnh góc nghiêng 5-35°] <--+
          |
[Đốt thực nghiệm 28 ô] -> [Ghi nhận khí tượng & Vi trắc ảnh]
                                   |
[Hồi quy đơn biến Excel/R] -------> [Hồi quy đa biến ANOVA] -> [Bảng quy chuẩn PCCCR]

Kế hoạch thực hiện và Quản trị rủi ro:

  • Giai đoạn 1 (Tuần 1 - 4): Điều tra cấu trúc thảm cháy tại 3 ô tiêu chuẩn điển hình ($30\text{ m} \times 33.3\text{ m}$) rừng thông 2 lá và 3 lá; lập 4 ô dạng bản ($5\text{ m} \times 5\text{ m}$) và các ô phụ $1\text{ m}^2$.
  • Giai đoạn 2 (Tuần 5 - 8): Gia công thiết bị sàn đốt, căn chỉnh hệ thống quang học vi trắc, thu gom và phơi đồng hóa mẫu lá thông.
  • Giai đoạn 3 (Tuần 9 - 12): Tiến hành 28 lượt đốt thử nghiệm đối chứng (7 mức độ dốc $\times$ 4 mức tải lượng vật liệu), lấy mẫu kiểm tra độ ẩm tức thời.
  • Giai đoạn 4 (Tuần 13 - 16): Xử lý ảnh trắc vi, trích xuất dữ liệu tọa độ lửa, chạy phân tích hồi quy và biên soạn biểu quy chuẩn.
  • Biện pháp giảm thiểu rủi ro: Thiết lập đường băng trắng cản lửa rộng $10\text{ m}$ bao quanh khu vực thử nghiệm tại Sóc Sơn; bố trí bình chữa cháy bột ABC 8kg, máy bơm nước áp lực cao dập lửa cô lập tức thì sau mỗi lần đo.

Implementation và kết quả

Development process

Quá trình đốt thử nghiệm được thực hiện trên 28 ô tiêu chuẩn với thiết kế ma trận biến số trực giao nhằm khảo sát toàn diện các trạng thái kết hợp giữa độ dốc $\theta \in [5^\circ, 35^\circ]$ (bước nhảy $5^\circ$) và khối lượng vật liệu $M_{vl} \in [5, 20]\text{ tấn/ha}$ (bước nhảy $5\text{ tấn/ha}$, tương đương rải $4\text{ kg}, 8\text{ kg}, 12\text{ kg}, 16\text{ kg}$ lá thông trên sàn $8\text{ m}^2$).

Thuật toán tính toán tốc độ lan tràn $V_c$ và chiều cao ngọn lửa $H_l$ dựa trên kỹ thuật phân tích ảnh vi trắc đa pha:

import numpy as np
import pandas as pd

def calculate_fire_behavior(delta_x_meters, delta_t_seconds, flame_pixels, scale_ratio):
    """
    Tính toán các thông số hành vi đám cháy từ dữ liệu vi trắc ảnh
    - delta_x_meters: Quãng đường lửa di chuyển đo được (m)
    - delta_t_seconds: Thời gian giữa 2 pha chụp từ đồng hồ điện tử (s)
    - flame_pixels: Chiều cao ngọn lửa tính bằng pixel trên ảnh
    - scale_ratio: Tỷ lệ quy đổi pixel sang mét thực tế (m/pixel)
    """
    # 1. Vận tốc lan tràn (Rate of Spread - Vc)
    vc = delta_x_meters / delta_t_seconds  # m/s
    
    # 2. Chiều cao ngọn lửa thực tế (Hl)
    hl = flame_pixels * scale_ratio  # m
    
    return vc, hl

def byram_fire_intensity(H_heat_kj_kg, W_fuel_load_ton_ha, R_spread_rate_m_s):
    """
    Tính toán cường độ đám cháy theo mô hình Byram (1959)
    - H_heat_kj_kg: Nhiệt lượng cháy của lá thông (~18,000 kJ/kg)
    - W_fuel_load_ton_ha: Khối lượng vật liệu cháy (tấn/ha) -> đổi sang kg/m2
    - R_spread_rate_m_s: Tốc độ lan tràn (m/s) -> đổi sang m/phút
    """
    w_kg_m2 = W_fuel_load_ton_ha * 0.1
    r_m_min = R_spread_rate_m_s * 60.0
    # Cường độ I (kW/m hoặc kJ/m/s)
    intensity_kw_m = (H_heat_kj_kg * w_kg_m2 * (r_m_min / 60.0))
    return intensity_kw_m

Dữ liệu ẩm độ vật liệu cháy được xác định nghiêm ngặt theo phương pháp cân sấy khối lượng ở $105^\circ\text{C}$:

  • Ẩm độ tương đối: $W_a = \frac{m_1 - m_0}{m_1} \times 100%$
  • Ẩm độ tuyệt đối: $W_0 = \frac{m_1 - m_0}{m_0} \times 100%$ (Trong đó $m_1$ là khối lượng vật liệu tự nhiên ẩm, $m_0$ là khối lượng khô kiệt sau sấy).

Testing và validation

Kết quả khảo sát thành phần thảm khô tại các ô tiêu chuẩn (OTC 01, OTC 02 - Thông 3 lá; OTC 03 - Thông 2 lá) cho thấy thành phần lá kim chiếm ưu thế tuyệt đối:

Thành phần thảm khô dưới tán rừng thông:
+-------------------------------------------------------------+
| [Lá thông: 93.0%]                 [Cành: 4.5%] [Quả: 1.6%] [Vỏ: 0.9%] |
+-------------------------------------------------------------+
(Kích thước lá: Chiều dài trung bình 18.2 - 20.4 cm; Đường kính: 1.0 - 1.2 mm)

Kết quả đo đạc từ 28 lần đốt thử nghiệm trong điều kiện thời tiết thực tế ($T_{kk}: 27.0 - 35.0^\circ\text{C}$, $W_{kk}: 67.0 - 100%$, $V_g: 0.2 - 3.5\text{ m/s}$) cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác:

                    QUAN HỆ GIỮA ĐỘ DỐC VÀ TỐC ĐỘ LAN TRÀN ĐÁM CHÁY
   Vc (m/s)
    0.14 |                                                         * (35°, 20t/ha)
    0.12 |                                                 *
    0.10 |                                         *       *
    0.08 |                                 *       *
    0.06 |                         *       *
    0.04 |         *       *       *
    0.02 | *       *
    0.00 +----------------------------------------------------------------->
           5°      10°     15°     20°     25°     30°     35°     Độ dốc (θ)

Bảng tổng hợp phân tích hồi quy đơn biến:

Cặp nhân tố phân tích Dạng phương trình tối ưu Hệ số $R^2$ Hệ số tương quan $R$ Ý nghĩa vật lý / Nhận xét
$H_l$ theo Khối lượng ($M_{vl}$) $y = 0.214 \cdot x^{0.482}$ $0.624$ $0.790$ Tải lượng càng lớn, ngọn lửa càng cao theo hàm lũy thừa
$H_l$ theo Ẩm độ VL ($W_{vl}$) $y = -0.024x + 1.241$ $0.456$ $0.675$ Ẩm độ kìm hãm mạnh chiều cao ngọn lửa
$H_l$ theo Tốc độ gió ($V_g$) $y = -0.1009x + 0.8281$ $0.1025$ $0.320$ Gió thổi tạt làm nghiêng và giảm chiều cao biểu kiến
$V_c$ theo Tốc độ gió ($V_g$) $y = 0.0046x + 0.0338$ $0.0289$ $0.170$ Gió đẩy ngọn lửa nghiêng về phía trước, gia nhiệt vật liệu
$H_l$ theo Nhiệt độ ($T_{kk}$) $y = 0.0237x + 0.3696$ $0.3699$ $0.608$ Nhiệt độ môi trường cao thúc đẩy bốc thoát hơi nước nhanh
$H_l$ theo Ẩm độ KK ($W_{kk}$) $y = -1.178\ln(x) + 5.8326$ $0.1538$ $0.392$ Tương quan nghịch phi tuyến dạng logarit
$V_c$ theo Độ dốc ($\theta$) $y = 0.0028e^{0.087\theta}$ $0.782$ $0.884$ Tốc độ tăng phi mã khi vượt ngưỡng dốc $20^\circ - 25^\circ$

Kết quả đạt được

Đề tài đã xây dựng thành công các phương trình hồi quy tuyến tính đa biến có ý nghĩa thống kê cao ($p\text{-value} < 0.001$, $F\text{-significance} \ll 0.05$):

1. Phương trình xác định Chiều cao ngọn lửa ($H_l$, đơn vị m): $$H_l = 0.0412 \cdot M_{vl} - 0.0185 \cdot W_{vl} - 0.0824 \cdot V_g + 0.0194 \cdot T_{kk} + 0.0115 \cdot \theta + 0.245$$

2. Phương trình xác định Tốc độ lan tràn đám cháy ($V_c$, đơn vị m/s): $$V_c = 0.0018 \cdot M_{vl} - 0.0009 \cdot W_{vl} + 0.0035 \cdot V_g + 0.0021 \cdot \theta - 0.0124$$

Bảng quy chuẩn an toàn chữa cháy và khoảng cách tiếp cận:

Độ dốc mặt đất ($\theta$) Tải lượng $M_{vl}$ ($5\text{ tấn/ha}$) Tải lượng $M_{vl}$ ($10\text{ tấn/ha}$) Tải lượng $M_{vl}$ ($15\text{ tấn/ha}$) Tải lượng $M_{vl}$ ($20\text{ tấn/ha}$) Phương thức chữa cháy khuyến nghị
$5^\circ - 10^\circ$ $H_l \le 0.6\text{ m}$; $D_{tc} \ge 1.5\text{ m}$ $H_l \le 0.8\text{ m}$; $D_{tc} \ge 2.0\text{ m}$ $H_l \le 1.1\text{ m}$; $D_{tc} \ge 2.5\text{ m}$ $H_l \le 1.3\text{ m}$; $D_{tc} \ge 3.0\text{ m}$ Dập trực tiếp bằng dụng cụ thủ công (bàn cào, vỉ dập)
$15^\circ - 20^\circ$ $H_l \le 0.8\text{ m}$; $D_{tc} \ge 2.0\text{ m}$ $H_l \le 1.1\text{ m}$; $D_{tc} \ge 2.8\text{ m}$ $H_l \le 1.4\text{ m}$; $D_{tc} \ge 3.5\text{ m}$ $H_l \le 1.6\text{ m}$; $D_{tc} \ge 4.2\text{ m}$ Bán cơ giới, máy thổi gió kết hợp tạo đường băng cản
$25^\circ - 35^\circ$ $H_l \le 1.2\text{ m}$; $D_{tc} \ge 3.0\text{ m}$ $H_l \le 1.5\text{ m}$; $D_{tc} \ge 4.0\text{ m}$ $H_l \ge 1.8\text{ m}$; $D_{tc} \ge 5.5\text{ m}$ $H_l \ge 2.2\text{ m}$; Không tiếp cận trực diện Chữa cháy cơ giới từ xa, đốt đón đầu (đốt ngược), mở băng cản

(Ghi chú: $D_{tc}$ là Khoảng cách tiếp cận an toàn tối thiểu của nhân lực).

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

Nghiên cứu mang lại những đóng góp khoa học và giá trị thực tiễn nổi bật cho chuyên ngành Quản lý Tài nguyên Rừng & PCCCR:

+---------------------------------------------------------------------------------+
|                    SO SÁNH ĐÓNG GÓP VỚI CÁC NGHIÊN CỨU TIỀN NHIỆM                |
+---------------------------------------------------------------------------------+
| Chỉ tiêu so sánh      | Mô hình Nesterov / McArthur     | Đề tài (Hoàng Văn Thành, 2017) |
+-----------------------+---------------------------------+--------------------------------+
| Yếu tố vi địa hình    | Bỏ qua độ dốc hoặc định tính    | Định lượng dốc 5° - 35°        |
| Vật liệu cháy cụ thể  | Rừng ôn đới / Bạch đàn Úc       | Rừng thông 2 lá, 3 lá VN       |
| Phương pháp đo lường  | Quan sát mắt thường / Cảm quan  | Vi trắc ảnh kỹ thuật số đa pha |
| Giá trị ứng dụng      | Dự báo cấp cháy chung           | Định mức nhân lực & khoảng cách|
+---------------------------------------------------------------------------------+
  • Làm rõ cơ chế đối lưu góc nghiêng: Chứng minh bằng thực nghiệm rằng khi độ dốc tăng từ $10^\circ$ lên $30^\circ$, vận tốc lan tràn $V_c$ tăng từ $180%$ đến $340%$ tùy thuộc vào tải lượng thảm lá. Lý do là góc nghiêng thu hẹp khoảng cách giữa ngọn lửa và lớp thảm khô chưa cháy phía trước, biến bức xạ nhiệt và dòng khí nóng đối lưu thành nguồn sấy khô vật liệu siêu tốc.
  • Tính toán định mức nhân lực dập lửa chuẩn hóa: Đưa ra công thức tính chu vi đám cháy $P(t) = \pi \cdot \sqrt{2(a^2 + b^2)}$ theo elip phát triển dưới ảnh hưởng của dốc, từ đó tính toán chính xác số người cần thiết cho từng mét chu vi lửa, loại bỏ phương thức điều động cảm tính.

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

Hệ thống kết quả nghiên cứu và bảng tra cứu kỹ thuật có khả năng chuyển giao và áp dụng trực tiếp tại các Ban Quản lý Rừng phòng hộ, Vườn Quốc gia và Chi cục Kiểm lâm:

+-------------------------------------------------------------------------+
|                QUY TRÌNH TRIỂN KHAI VẬN HÀNH TRÊN THỰC ĐỊA              |
|                                                                         |
|  [Trạm quan trắc vi khí tượng] + [Bản đồ độ dốc DEM & Tải lượng rừng]   |
|                                     |                                   |
|                                     v                                   |
|       [Truy xuất phương trình hồi quy Vc & Hl tại tọa độ nghi ngờ]      |
|                                     |                                   |
|                                     v                                   |
|  +-------------------------------------------------------------------+  |
|  | Phân cấp phản ứng:                                                |  |
|  | - Dốc < 15°, Tải < 10 tấn/ha: Đội tuần tra cơ sở (5-7 người)      |  |
|  | - Dốc 15-25°, Tải 10-15 tấn/ha: Đội xung kích + Máy thổi (15-20 ng)|  |
|  | - Dốc > 25°, Tải > 15 tấn/ha: Báo động cấp 2, lập phòng tuyến chặn|  |
|  +-------------------------------------------------------------------+  |
+-------------------------------------------------------------------------+

Kịch bản triển khai tình huống (Use Case Scenario):

  • Địa bàn: Lô rừng thông 3 lá, 30 năm tuổi tại Sóc Sơn, độ dốc bình quân sườn đồi $\theta = 25^\circ$, trữ lượng thảm khô đo đạc $M_{vl} = 15\text{ tấn/ha}$, độ ẩm thảm lá $W_{vl} = 15%$, tốc độ gió $V_g = 2.0\text{ m/s}$.
  • Tính toán tự động:
    • Chiều cao ngọn lửa dự báo: $H_l \approx 1.75\text{ m}$ (Vượt ngưỡng tiếp cận trực diện an toàn).
    • Tốc độ lan tràn dự báo: $V_c \approx 0.078\text{ m/s} = 4.68\text{ m/phút}$.
    • Sau 30 phút phát hiện, đám cháy tiến xa $140.4\text{ m}$, chu vi ước tính $\sim 450\text{ m}$.
  • Quyết định điều hành: Không triển khai dập trực tiếp tại đầu ngọn lửa; điều động 2 tổ cơ động (mỗi tổ 8 người) tiếp cận men theo 2 bên sườn đám cháy (nơi vận tốc lan tràn chỉ bằng $1/3$ tốc độ ngọn), đồng thời triển khai máy cưa xích và máy thổi gió cách tim lửa $300\text{ m}$ về phía đỉnh dốc để mở băng cản lửa trắng rộng $4\text{ m}$.

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

Dù đạt được những kết quả thực nghiệm có độ tin cậy cao, đề tài vẫn tồn tại một số điểm giới hạn kỹ thuật cần được tiếp tục hoàn thiện:

  • Hạn chế kỹ thuật:
    1. Quy mô sàn thực nghiệm $8\text{ m}^2$ ($2\text{ m} \times 4\text{ m}$) tuy phản ánh tốt đám cháy khởi đầu nhưng chưa bao hàm được hiện tượng tạo luồng gió xoáy cục bộ (fire whirl) và hiện tượng tàn lửa bay xa (spotting fire) như trong các đại hỏa hoạn thực tế.
    2. Nghiên cứu mới khảo sát chiều lan tràn ngược dốc (lửa tiến lên dốc); chưa thu thập đủ mẫu thực nghiệm đối với trường hợp lửa lan xuôi dốc (xuống dốc) hoặc lan ngang sườn dốc (flank fire).
  • Hướng phát triển tiếp theo:
    • Tích hợp các hệ số hồi quy thực nghiệm vào hệ thống thông tin địa lý (GIS) và mô hình hóa không gian 3D tự động trên nền tảng WebGIS phục vụ cảnh báo trực tuyến thời gian thực.
    • Mở rộng đối tượng vật liệu cháy cho rừng keo (Acacia), rừng tràm (Melaleuca) và thảm cỏ tranh (Imperata cylindrica).

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

  • Sinh viên & Học viên Lâm nghiệp: Bộ dữ liệu thực nghiệm chuẩn mực về tương quan nhiệt - ẩm - địa hình; phương pháp luận thiết kế thí nghiệm hiện trường mẫu mực có tính ứng dụng cao.
  • Kỹ sư Quản lý Tài nguyên & Kiểm lâm địa phương: Nắm bắt chính xác bảng tra khoảng cách an toàn ($D_{tc}$) và tốc độ lan tràn ($V_c$) để ra quyết định tác chiến bảo đảm tuyệt đối tính mạng cho lực lượng chữa cháy.
  • Chủ rừng & Doanh nghiệp Lâm nghiệp: Cơ sở khoa học để thiết kế mật độ băng cản lửa, đường ranh cản lửa và kế hoạch đốt trước vật liệu cháy có điều khiển dưới tán rừng thông vào đầu mùa khô một cách chủ động, tiết kiệm tới $40%$ ngân sách dọn thực bì hàng năm.
  • Cộng đồng Nghiên cứu Hỏa sinh thái: Bổ sung các hệ số thực nghiệm bản địa của thảm thực vật Việt Nam vào cơ sở dữ liệu mô hình cháy rừng khu vực Đông Nam Á.

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

1. Tại sao độ dốc lại làm tốc độ cháy rừng tăng vọt ngay cả khi không có gió mạnh?

Khi cháy trên mặt phẳng nằm ngang ($0^\circ$), phần lớn nhiệt lượng tỏa ra bay thẳng lên khí quyển theo phương thẳng đứng. Tuy nhiên trên sườn dốc ($\theta > 0^\circ$), khoảng cách hình học giữa ngọn lửa và lớp thảm khô phía trên bị thu hẹp đáng kể. Ngọn lửa nghiêng song song với mặt đất, truyền nhiệt trực tiếp bằng bức xạ và dòng khí đối lưu nóng, sấy khô kiệt lượng ẩm trong lá thông chỉ trong vài giây trước khi ngọn lửa chạm tới, tạo điều kiện cho phản ứng bốc cháy diễn ra gần như đồng thời.

2. Có thể sử dụng bảng tra cứu của đề tài này cho rừng keo hoặc rừng bạch đàn không?

Có thể sử dụng như một tài liệu tham khảo về mặt xu hướng động học (độ dốc làm tăng tốc độ lan tràn). Tuy nhiên, để tính toán chính xác chỉ số $H_l$ và $V_c$, cần hiệu chỉnh lại tham số khối lượng và đặc tính bắt cháy do lá bạch đàn và lá keo có diện tích bề mặt, hàm lượng tinh dầu và độ ẩm phân hủy khác biệt so với lá thông kim.

3. Kỹ thuật vi trắc ảnh trong đề tài có ưu điểm gì so với việc bấm giờ truyền thống?

Phương pháp vi trắc ảnh kết hợp chụp nhiều pha với hệ thống thước chuẩn nhôm và đồng hồ giây kỹ thuật số cho phép loại bỏ hoàn toàn sai số do phản xạ thị giác của con người (thường trễ $0.3 - 0.5\text{ s}$). Kỹ thuật này cho phép đo chính xác biên dạng ngọn lửa ở trạng thái động và bóc tách tọa độ mép lửa sát mặt đất chính xác đến từng milimét.

4. Ẩm độ vật liệu cháy bao nhiêu thì rừng thông bước vào ngưỡng nguy hiểm cực kỳ dễ cháy?

Qua số liệu sấy và đốt thử nghiệm, khi ẩm độ tương đối của lá thông khô ($W_a$) giảm xuống dưới $18%$, thảm vật liệu bước vào trạng thái bén lửa tức thì khi gặp mồi. Khi $W_a < 12%$, tốc độ cháy lan đạt cực đại và ngọn lửa có khả năng tự duy trì mãnh liệt bất chấp độ ẩm không khí xung quanh cao.

5. Tại sao khi khối lượng vật liệu tăng quá cao ở địa hình bằng phẳng, tốc độ lan tràn lại có xu hướng chững lại?

Khi khối lượng thảm khô quá dày ở địa hình bằng, nhu cầu ôxy cung cấp cho quá trình oxy hóa ở lớp đáy sát mặt đất bị thiếu hụt nghiêm trọng do bị lớp thảm phía trên che phủ. Khi đó tốc độ lan tràn chuyển từ trạng thái phụ thuộc vào nhiệt năng sang phụ thuộc vào tốc độ khuếch tán ôxy khí quyển vào vùng phản ứng cháy.

Kết luận

Khóa luận tốt nghiệp "Nghiên cứu ảnh hưởng của độ dốc đến tốc độ lan tràn của đám cháy rừng" của tác giả Hoàng Văn Thành (dưới sự hướng dẫn của GS.TS Vương Văn Quỳnh) đã giải quyết trọn vẹn và chuẩn xác một trong những khoảng trống kỹ thuật quan trọng nhất của công tác phòng chống cháy rừng tại Việt Nam. Bằng việc thực nghiệm hóa 28 kịch bản đốt trên sàn tôn chuyên dụng $8\text{ m}^2$ ở các góc nghiêng $5^\circ - 35^\circ$, công trình đã chứng minh mối tương quan phi tuyến chặt chẽ giữa độ dốc địa hình và hành vi lửa rừng thông, thiết lập hệ thống phương trình hồi quy đa biến và bảng định mức an toàn có độ chính xác cao.

Kết quả của đề tài không chỉ đóng góp giá trị lý luận quan trọng cho ngành Quản lý Tài nguyên Rừng mà còn cung cấp một công cụ tác chiến thực địa thiết thực, giúp các nhà quản lý rừng chủ động dự báo, thiết kế băng cản lửa tối ưu và điều hành lực lượng chữa cháy an toàn, hiệu quả.