Giới thiệu dự án

Khu di tích lịch sử Pác Bó (xã Trường Hà, huyện Hà Quảng, tỉnh Cao Bằng) là Di tích Quốc gia Đặc biệt lưu giữ 42 di tích gốc gắn liền với giai đoạn hoạt động cách mạng 1941–1945 của Chủ tịch Hồ Chí Minh. Với tài nguyên tự nhiên đặc sắc gồm hệ thống núi đá vôi Karst, rừng đặc dụng che phủ $91{,}46%$ và dòng thủy văn suối Lê Nin, Pác Bó đón trên $25.000$ lượt khách/năm trong tổng số $650.000$ lượt khách toàn tỉnh Cao Bằng (tăng trưởng $18%$/năm, doanh thu đạt 115 tỷ VNĐ năm 2015). Sự gia tăng nhanh chóng của lượng khách du lịch và các hoạt động nhân sinh đặt ra áp lực lớn lên chất lượng môi trường sinh thái (MTST).

+-------------------------------------------------------------------------+
|                  KHU DI TÍCH QUỐC GIA ĐẶC BIỆT PÁC BÓ                   |
| (Quy hoạch: 500 ha | Sử dụng: 295 ha | 42 di tích gốc | >25.000 khách)  |
+------------------------------------+------------------------------------+
                                     |
                +--------------------+--------------------+
                | Áp lực nhân sinh đa chiều tác động     |
                +--------------------+--------------------+
                                     |
      +------------------------------+------------------------------+
      |                              |                              |
+-----v---------------+    +---------v-----------+    +-------------v-----+
| Hoạt động Du lịch   |    | Sản xuất Nông nghiệp|    | Dân sinh & Thủy sản|
| - Rác thải, nước thải|    | - Dư lượng N, P     |    | - Khai thác gỗ củi|
| - Biến động cảnh quan|    | - Thuốc BVTV        |    | - Đánh bắt trái phép|
| - Áp lực quá tải     |    | - Xói mòn đất dốc   |    | - Xả thải sinh hoạt |
+---------------------+    +---------------------+    +-------------------+

Vấn đề nghiên cứu và Hiện trạng cấp bách

  1. Ô nhiễm cục bộ nguồn nước mặt: Lưu lượng xả thải chưa qua xử lý tập trung từ các điểm dịch vụ và 122 hộ dân bản Pác Bó trực tiếp đe dọa chỉ số Oxy hòa tan (DO) và gia tăng Nhu cầu Oxy hóa học (COD), Nhu cầu Oxy sinh hóa ($BOD_5$) tại suối Lê Nin.
  2. Biến động sử dụng đất và suy thoái thổ nhưỡng: Diện tích đất sản xuất nông nghiệp giảm $1{,}14\text{ ha}$ (giai đoạn 2014–2015), thúc đẩy áp lực thâm canh hóa học và chuyển dịch sinh kế sang khai thác tài nguyên rừng đặc dụng ($2.527{,}96\text{ ha}$).
  3. Xung đột giữa bảo tồn di sản và thương mại hóa du lịch: Hạ tầng thu gom chất thải rắn (CTR) chưa đồng bộ với tỷ lệ thu gom phân loại tại nguồn đạt dưới $60%$, gây nguy cơ suy thoái thạch quyển và hiện tượng phú dưỡng hóa thủy vực.

Mục tiêu đề tài

  1. Đánh giá toàn diện hiện trạng chất lượng MTST (đất, nước mặt suối Lê Nin, chất thải rắn) tại Khu di tích lịch sử Pác Bó.
  2. Định lượng các luồng phát thải và tác động tiêu cực từ hoạt động du lịch, canh tác nông nghiệp, lâm nghiệp và khai thác thủy sản.
  3. Đánh giá thực trạng hiệu quả công tác quản lý môi trường sinh thái của Ban quản lý khu di tích và chính quyền địa phương.
  4. Xây dựng hệ thống giải pháp kỹ thuật, chính sách và sinh kế cộng đồng nhằm giảm thiểu suy thoái môi trường, duy trì sức tải sinh thái tự nhiên.

Phạm vi và Giới hạn nghiên cứu

  • Không gian: Toàn bộ ranh giới quy hoạch $500\text{ ha}$ (diện tích vùng lõi sử dụng $295\text{ ha}$) thuộc xã Trường Hà, huyện Hà Quảng, tỉnh Cao Bằng.
  • Thời gian khảo sát: Khảo sát thực địa từ tháng 09/2015 đến tháng 11/2015; tổng hợp, mô hình hóa số liệu hoàn thành trong tháng 05/2016.
  • Đối tượng: Môi trường vật lý (nước mặt, thổ nhưỡng, CTR), đa dạng sinh học rừng đặc dụng và cộng đồng 122 hộ dân ($513$ nhân khẩu) thôn Pác Bó.

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

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

Tiêu chí phân tích Phương pháp truyền thống (Báo cáo định kỳ BQL) Mô hình nghiên cứu của đề tài (Đa chiều tích hợp)
Phạm vi giám sát Thống kê số lượng khách và thu gom rác đơn thuần Tích hợp giám sát lý - hóa nước mặt, biến động sử dụng đất và xã hội học
Chuẩn đối chiếu Đánh giá cảm quan, thiếu chỉ tiêu kỹ thuật Đối chuẩn $QCVN\ 08:2008/BTNMT$ (Cột A1/B1) và chỉ số $BOD_5/COD$
Xử lý số liệu Thống kê mô tả thủ công qua sổ ghi chép Phân tích hồi quy, chuẩn hóa dữ liệu thực địa trên SPSS v22.0 và Excel 2016
Khả năng dự báo Không có mô hình dự báo sức chứa môi trường Thiết lập công thức ước tính tải lượng ô nhiễm hữu cơ và sức chứa du lịch

Đánh giá yêu cầu theo mô hình MoSCoW

  • Must have (Bắt buộc): Thiết lập trạm quan trắc định kỳ chỉ số DO, COD, $BOD_5$ tại đập tràn suối Lê Nin; quy hoạch hệ thống thùng rác phân loại 03 ngăn đạt chuẩn tại $100%$ các điểm di tích gốc (Hang Cốc Bó, Bàn đá, Cột mốc 108, Lán Khuổi Nặm).
  • Should have (Nên có): Xây dựng hệ thống mương gom và bể tự hoại 3 ngăn hợp quy cho $100%$ hộ kinh doanh dịch vụ ăn uống; tập huấn kỹ thuật giảm dư lượng Nitrate, Phosphate trong canh tác nông nghiệp dốc.
  • Could have (Có thể có): Ứng dụng hệ thống GIS (ArcGIS v10.3) lập bản đồ phân vùng nhạy cảm sinh thái; lắp đặt bảng thông tin điện tử cảnh báo sức chứa du lịch tức thời.
  • Won't have (Chưa thực hiện): Xây dựng nhà máy xử lý nước thải tập trung công suất lớn (không khả thi do địa hình núi đá chia cắt mạnh và nguồn vốn hạn chế).

Thiết kế hệ thống đánh giá và quan trắc

+-------------------------------------------------------------------------+
|                KIẾN TRÚC QUAN TRẮC VÀ ĐÁNH GIÁ TẢI LƯỢNG                |
+-------------------------------------------------------------------------+

Công thức đánh giá sức chứa tải lượng du lịch

Sức chứa thực tế ($RCC$) của khu di tích được xác định thông qua Sức chứa vật lý ($PCC$) kết hợp hệ số hiệu chỉnh môi trường sinh thái ($C_f$):

$$PCC = A \times \frac{V}{a} \times R_f$$

$$RCC = PCC \times (1 - C_{f1}) \times (1 - C_{f2}) \times \dots \times (1 - C_{fn})$$

Trong đó:

  • $A$: Diện tích khu vực tham quan hữu dụng ($m^2$).
  • $V/a$: Định mức không gian cho 1 khách ($1\text{ khách}/5m^2$).
  • $R_f$: Hệ số luân chuyển theo thời gian mở cửa ($1\text{ ngày} = 8\text{ giờ} / 2\text{ giờ/lượt} = 4\text{ lượt/ngày}$).
  • $C_{f1}, C_{f2}$: Các hệ số suy giảm do lượng mưa tập trung ($80%$ vào tháng 7–9), xói mòn sườn dốc Karst và bảo vệ di tích gốc.

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

Giai đoạn 1: Điều tra thứ cấp (09/2015)
Giai đoạn 2: Khảo sát thực địa & Đo đạc (10/2015 - 11/2015)
Giai đoạn 3: Phân tích & Xử lý số liệu (04/2016 - 05/2016)
  • Đảm bảo chất lượng dữ liệu (QA/QC): Tất cả mẫu nước được bảo quản lạnh ở $4^\circ\text{C}$ và phân tích theo phương pháp tiêu chuẩn Standard Methods for the Examination of Water and Wastewater (22nd Edition); phiếu điều tra được làm sạch và kiểm tra chéo (Cross-check) loại bỏ phương sai ngoại lai trước khi chạy thống kê.

Implementation và kết quả

Quy trình phát triển công cụ phân tích và xử lý dữ liệu

Mã nguồn Python sử dụng các thư viện pandas và numpy để chuẩn hóa, tính toán chỉ số tỷ lệ $BOD_5/COD$ và phân loại tự động chất lượng nước mặt suối Lê Nin dựa trên ngưỡng $QCVN\ 08:2008/BTNMT$:

import numpy as np
import pandas as pd

# Dữ liệu quan trắc nước mặt suối Lê Nin (Khu vực Đập Tràn)
water_monitoring_data = {
    "Year": [2013, 2014, 2015],
    "DO_mg_L": [7.82, 7.45, 6.95],
    "BOD5_mg_L": [2.10, 3.40, 4.85],
    "COD_mg_L": [4.20, 6.80, 9.60],
    "TSS_mg_L": [12.5, 18.2, 24.6],
}

df_water = pd.DataFrame(water_monitoring_data)

# Tính toán tỷ số Biodegradability Index (BOD5 / COD)
df_water["Biodegradability_Index"] = df_water["BOD5_mg_L"] / df_water["COD_mg_L"]


# Hàm đánh giá chất lượng nguồn nước mặt theo QCVN 08:2008/BTNMT (Cột A1)
def evaluate_water_quality(row):
    # QCVN 08:2008/BTNMT Cột A1: DO >= 6 mg/L, BOD5 <= 4 mg/L, COD <= 10 mg/L
    if row["DO_mg_L"] >= 6.0 and row["BOD5_mg_L"] <= 4.0 and row["COD_mg_L"] <= 10.0:
        return "Loại A1 (Bảo tồn thủy sinh & Cấp nước sinh hoạt)"
    # QCVN 08:2008/BTNMT Cột A2: DO >= 5 mg/L, BOD5 <= 6 mg/L, COD <= 15 mg/L
    elif (
        row["DO_mg_L"] >= 5.0 and row["BOD5_mg_L"] <= 6.0 and row["COD_mg_L"] <= 15.0
    ):
        return "Loại A2 (Cấp nước sinh hoạt sau xử lý phù hợp)"
    else:
        return "Loại B (Cảnh báo ô nhiễm hữu cơ nhẹ)"


df_water["Quality_Status"] = df_water.apply(evaluate_water_quality, axis=1)
print(
    df_water[
        [
            "Year",
            "DO_mg_L",
            "BOD5_mg_L",
            "COD_mg_L",
            "Biodegradability_Index",
            "Quality_Status",
        ]
    ]
)

Kiểm định và đánh giá thực nghiệm

Kết quả phân tích mẫu nước tại đập tràn suối Lê Nin qua 3 năm (2013–2015) ghi nhận xu hướng gia tăng rõ rệt của hàm lượng chất hữu cơ do áp lực xả thải sinh hoạt và du lịch:

Nồng độ (mg/L)
10 |                                                [COD: 9.60]
 8 |  [DO: 7.82]               [DO: 7.45]
 6 |              [COD: 6.80]               [DO: 6.95]
 4 |  [COD: 4.20]              [BOD5: 3.40]
 2 |  [BOD5: 2.10]                                  [BOD5: 4.85]
         Năm 2013                  Năm 2014             Năm 2015
  • Chỉ số DO (Oxy hòa tan): Giảm từ $7{,}82\text{ mg/L}$ (2013) xuống $6{,}95\text{ mg/L}$ (2015) (giảm $11{,}12%$).
  • Chỉ số $BOD_5$: Tăng từ $2{,}10\text{ mg/L}$ lên $4{,}85\text{ mg/L}$ (tăng $130{,}95%$, vượt ngưỡng Cột A1 của $QCVN\ 08:2008/BTNMT$).
  • Chỉ số COD: Tăng từ $4{,}20\text{ mg/L}$ lên $9{,}60\text{ mg/L}$ (tăng $128{,}57%$).

Khảo sát định lượng ý kiến cộng đồng và du khách

Kết quả phỏng vấn cấu trúc $30$ hộ dân, $10$ cán bộ quản lý và $30$ du khách tại Pác Bó:

Nhóm tác động được khảo sát Tỷ lệ đánh giá "Ảnh hưởng đáng kể" (%) Tác nhân chính được xác định
Ô nhiễm Rác thải sinh hoạt $76{,}67%$ ($23/30$ hộ dân) Du khách vứt rác tùy tiện; thiếu thùng rác dọc đường mòn
Biến đổi chất lượng nước $63{,}33%$ ($19/30$ hộ dân) Rửa đồ dùng sinh hoạt, nước thải từ hàng quán ăn uống
Xâm lấn thảm thực vật $53{,}33%$ ($16/30$ hộ dân) Khai thác củi đốt, dẫm đạp cây tái sinh tại khu rừng đặc dụng
Ô nhiễm tiếng ồn & Giao thông $43{,}33%$ ($13/30$ hộ dân) Phương tiện chở khách tập trung vào các dịp lễ cao điểm

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

  1. Phương pháp tiếp cận liên ngành (Socio-Ecological Framework): Đề tài là công trình đầu tiên tại Pác Bó kết hợp giữa phân tích hóa lý môi trường ($BOD_5$, COD, DO), dữ liệu biến động trắc địa đất đai ($GIS/Cadastral\ data$) và điều tra xã hội học cộng đồng bản địa Tày, Nùng.
  2. Định lượng hóa biến động sử dụng đất gắn với áp lực môi trường: Chứng minh mối quan hệ giữa việc giảm $1{,}14\text{ ha}$ đất sản xuất nông nghiệp với xu hướng gia tăng thâm canh phân bón hóa học ($NPK$) và áp lực chặt phá lấy lâm sản phụ tại $2.189{,}39\text{ ha}$ rừng tự nhiên đặc dụng.
  3. Mô hình quản lý chất thải rắn phân tán thích ứng địa hình Karst: Đề xuất mô hình thu gom phân kỳ kết hợp ủ phân hữu cơ vi sinh tại chỗ ($Composting$), giúp giảm tải $45%$ khối lượng CTR vận chuyển ra khỏi thung lũng Pác Bó.
+--------------------------------------------------------------------------+
|                      MÔ HÌNH SO SÁNH GIẢI PHÁP                           |
+------------------------------------+-------------------------------------+
| Nghiên cứu trước đây               | Đề tài Lê Thị Xuân Hồng (2016)      |
| (Phạm Trung Lương 2010;            |                                     |
|  Đặng Văn Huyến 2003)              |                                     |
+------------------------------------+-------------------------------------+
| - Đánh giá định tính chung         | - Định lượng hóa nồng độ BOD5/COD   |
| - Quy mô cấp vùng/toàn tỉnh        | - Dữ liệu thực địa vi mô cấp thôn   |
| - Giải pháp hành chính đơn thuần   | - Tích hợp giải pháp kinh tế - kỹ thuật|
+------------------------------------+-------------------------------------+

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

Kịch bản triển khai phân vùng bảo vệ môi trường Pác Bó

|               VÙNG LÕI (295 ha): BẢO TỒN NGUYÊN TRẠNG                    |
| - Cấm tuyệt đối xây dựng mới | Phục hồi thảm thực vật tự nhiên          |
| - Giới hạn khách tham quan <500 người/thời điểm                         |
|          VÙNG ĐỆM DỊCH VỤ DU LỊCH & KHU DÂN CƯ (BẢN PÁC BÓ)             |
| - Lắp đặt bể tự hoại 3 ngăn | Điểm thu gom CTR tập trung                |
| - Chuyển đổi mô hình canh tác hữu cơ | Tuyến xe điện trung chuyển        |

Dự toán chi phí và Hiệu quả kinh tế - sinh thái (Cost-Benefit Analysis)

  • Chi phí đầu tư hạ tầng xanh ban đầu: Ước tính khoảng $450\text{ triệu VNĐ}$ (gồm hệ thống thùng rác compost composite, trạm quan trắc nhanh cầm tay, biển chỉ dẫn giáo dục môi trường).
  • Lợi ích sinh thái - xã hội:
    • Giảm $60%$ nguy cơ ô nhiễm nguồn nước cấp sinh hoạt hạ lưu suối Lê Nin.
    • Tăng thời gian lưu trú của khách du lịch sinh thái từ $1{,}5\text{ giờ}$ lên $4\text{ giờ}$, nâng mức chi tiêu bình quân của du khách tại bản địa thêm $25%$.
    • Bảo tồn nguyên vẹn $100\text{ ha}$ đất rừng do các hộ dân trực tiếp quản lý theo Quyết định 126/QĐ-TTg.

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

Hạn chế kỹ thuật và nguồn lực

  • Chuỗi số liệu quan trắc ngắn: Dữ liệu quan trắc hóa lý nước mặt chỉ bao gồm 3 năm liên tiếp với tần suất $2\text{ đợt/năm}$, chưa phản ánh hết sự biến thiên theo chu kỳ thủy văn mùa lũ cực đoan.
  • Thiếu thiết bị đo tự động: Các chỉ tiêu khí thải và tiếng ồn chưa được đo đạc bằng máy đo chuyên dụng mà chủ yếu dựa trên khảo sát cảm quan và lưu lượng phương tiện.

Hướng phát triển đề xuất

  1. Ứng dụng công nghệ IoT: Lắp đặt hệ thống sensor quan trắc nước mặt tự động (DO, pH, độ đục, TDS) truyền dữ liệu thời gian thực về Trung tâm điều hành BQL.
  2. Xây dựng bản đồ số nhạy cảm sinh thái (GIS-based WebApp): Phân tích không gian đa tiêu chí (MCDA) để tối ưu hóa tuyến đường mòn du lịch sinh thái, tránh các vùng đất có độ dốc $>25^\circ$ dễ sạt lở.

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

+-------------------------------------------------------------------------+
|                         CÁC BÊN THỤ HƯỞNG                               |
+--------------------+--------------------+-------------------------------+
| Sinh viên & Học viên| - Bộ dữ liệu thực địa mẫu về MTST khu di tích     |
|                    | - Phương pháp luận đánh giá tác động du lịch      |
+--------------------+--------------------+-------------------------------+
| Ban quản lý Di tích| - Bộ chỉ số giám sát chất lượng nước suối Lê Nin  |
| & Địa phương       | - Khung quy chế quản lý chất thải rắn phân tán    |
+--------------------+--------------------+-------------------------------+
| Doanh nghiệp       | - Cơ sở quy hoạch các tour du lịch sinh thái bền  |
| Lữ hành            |   vững, không gây quá tải môi trường              |
+--------------------+--------------------+-------------------------------+
| Cộng đồng cư dân   | - Nâng cao thu nhập qua dịch vụ homestay sinh thái|
| bản địa (Tày, Nùng)| - Bảo tồn nguồn tài nguyên nước và rừng bền vững  |
+--------------------+--------------------+-------------------------------+

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

1. Nguồn thải nào gây nguy cơ lớn nhất đối với chất lượng nước suối Lê Nin?

Nước thải sinh hoạt chưa qua xử lý từ 122 hộ dân bản Pác Bó và các cơ sở ăn uống dọc tuyến di tích, kết hợp với dư lượng phân bón hóa học Nitrate/Phosphate bị rửa trôi từ $228{,}89\text{ ha}$ đất canh tác nông nghiệp khi có mưa lớn tập trung (chiếm $80%$ lượng mưa cả năm).

2. Biến động diện tích đất nông nghiệp ảnh hưởng thế nào đến môi trường sinh thái?

Việc đất sản xuất nông nghiệp giảm $1{,}14\text{ ha}$ trong giai đoạn 2014–2015 làm tăng tỷ lệ lao động dư thừa trong nông nghiệp, thúc đẩy nguy cơ khai thác lâm sản trái phép tại $2.527{,}96\text{ ha}$ rừng đặc dụng và gia tăng mức độ sử dụng thuốc BVTV trên phần diện tích còn lại để bù đắp năng suất.

3. Tỷ số $BOD_5/COD$ phản ánh điều gì về chất lượng nước suối Lê Nin?

Tỷ số $BOD_5/COD$ đạt mức $\approx 0{,}50$ vào năm 2015 cho thấy nước mặt có tỷ lệ chất hữu cơ dễ phân hủy sinh học cao. Đây là dấu hiệu đặc trưng của ô nhiễm từ nguồn nước thải sinh hoạt và dịch vụ du lịch, đòi hỏi can thiệp bằng các biện pháp xử lý sinh học hiếu khí/kỵ khí tự nhiên.

4. Giải pháp nào có tính khả thi cao nhất cho công tác quản lý chất thải rắn tại Pác Bó?

Mô hình thùng rác phân loại 3 ngăn gắn biển hướng dẫn song ngữ/biểu tượng trực quan tại 42 điểm di tích gốc, kết hợp chuyển giao kỹ thuật ủ rác hữu cơ vi sinh tại các hộ dân bản Pác Bó để giảm tải chi phí vận chuyển ngoại vi.

Kết luận

Đề tài khóa luận tốt nghiệp của tác giả Lê Thị Xuân Hồng đã cung cấp bức tranh toàn cảnh, thực chứng về mối quan hệ biện chứng giữa phát triển du lịch, sản xuất dân sinh và chất lượng môi trường sinh thái tại Di tích Quốc gia Đặc biệt Pác Bó. Các số liệu phân tích hóa lý nước suối Lê Nin ($BOD_5$ tăng $130{,}95%$, DO giảm $11{,}12%$) cùng các phát hiện về biến động đất đai là cơ sở khoa học quan trọng phục vụ công tác quy hoạch du lịch bền vững. Việc triển khai đồng bộ các giải pháp công nghệ sinh thái, phân vùng bảo tồn và phát triển sinh kế cộng đồng sẽ bảo vệ vững chắc cảnh quan "núi Các Mác, suối Lê Nin", giữ trọn vẹn giá trị lịch sử và sinh thái cho thế hệ tương lai.