Giới thiệu dự án

Bối cảnh và thực trạng ngành du lịch tỉnh Cao Bằng

Du lịch được định vị là ngành kinh tế mũi nhọn không khói, đóng vai trò đòn bẩy trong cơ cấu chuyển dịch kinh tế vùng trung du và miền núi phía Bắc. Tỉnh Cao Bằng sở hữu hơn 200 di tích lịch sử - văn hóa (gồm 3 di tích quốc gia đặc biệt, 25 di tích cấp quốc gia, 68 di tích cấp tỉnh) và danh hiệu Công viên Địa chất Toàn cầu UNESCO Non nước Cao Bằng (công nhận năm 2018). Theo Niên giám Thống kê tỉnh Cao Bằng, lượng khách du lịch đạt đỉnh gần 440.000 lượt vào năm 2019 và duy trì trên 300.000 lượt/năm sau đại dịch COVID-19; tốc độ tăng trưởng GRDP năm 2022 đạt 5,04%, trong đó khu vực dịch vụ tăng 6,01%.

Tuy nhiên, sự gia tăng nhanh chóng của lượng khách du lịch đang gây áp lực nghiêm trọng lên hạ tầng tự nhiên và xã hội. Thực trạng khai thác du lịch tại hai điểm nóng trọng điểm là Khu di tích Quốc gia đặc biệt Pác Bó và Khu du lịch sinh thái Thác Bản Giốc đối mặt với nhiều bất cập: ô nhiễm môi trường cục bộ, tính mùa vụ cao, thương mại hóa thiếu kiểm soát, và nguy cơ mai một bản sắc văn hóa dân tộc bản địa.

       [Tăng trưởng du khách] ---> [Áp lực tài nguyên & Môi trường]
                 |                               |
                 v                               v
       [Hạ tầng quá tải] <======> [Khai thác thiếu bền vững]
                 |
                 +--> Cần: Mô hình Kinh tế lượng & Đo lường thực nghiệm

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

Các cơ quan quản lý và đơn vị kinh doanh dịch vụ tại địa phương thiếu một mô hình định lượng chuẩn xác nhằm đánh giá mức độ ảnh hưởng của các yếu tố cấu thành đến phát triển du lịch bền vững (Sustainable Tourism Development). Việc ban hành chính sách quản lý du lịch còn mang tính định tính, thiếu căn cứ dữ liệu thực nghiệm, dẫn đến các chính sách hỗ trợ và đầu tư hạ tầng chưa đạt hiệu quả tối ưu.

Mục tiêu dự án

  1. Hệ thống hóa cơ sở lý luận về du lịch bền vững theo tiêu chuẩn quốc tế (GSTC, UNWTO) và Luật Du lịch Việt Nam 2017.
  2. Xây dựng mô hình kinh tế lượng đo lường các nhân tố tác động đến phát triển du lịch bền vững tại Cao Bằng.
  3. Thu thập, xử lý và phân tích dữ liệu thực nghiệm sơ cấp ($N = 130$) bằng phương pháp định lượng đa biến.
  4. Kiểm định giả thuyết nghiên cứu và lượng hóa mức độ tác động của từng nhân tố.
  5. Đề xuất hệ thống giải pháp và chính sách có căn cứ khoa học cho giai đoạn 2023–2030.

Phương pháp tiếp cận và kết quả kỳ vọng

Nghiên cứu ứng dụng kết hợp phương pháp định tính và định lượng, sử dụng thang đo Likert 5 điểm, phân tích độ tin cậy Cronbach's Alpha, phân tích nhân tố khám phá (EFA), hồi quy tuyến tính bội (OLS), và kiểm định sai biệt One-way ANOVA trên phần mềm IBM SPSS Statistics 26.0. Kết quả đầu ra là phương trình hồi quy chuẩn hóa xác định rõ trọng số tác động của từng nhóm nhân tố, hỗ trợ việc ra quyết định chính sách công.

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

  • Không gian: Khu di tích Quốc gia đặc biệt Pác Bó (huyện Hà Quảng) và Khu du lịch sinh thái Thác Bản Giốc (huyện Trùng Khánh), tỉnh Cao Bằng.
  • Thời gian: Dữ liệu thứ cấp giai đoạn 2017–2022; dữ liệu sơ cấp thu thập năm 2023.
  • Đối tượng: Du khách nội địa và quốc tế trực tiếp trải nghiệm dịch vụ tại các điểm đến.

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

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

Mô hình / Giải pháp Ưu điểm Nhược điểm Khả năng ứng dụng tại Cao Bằng
Thongphon et al. (2013) (Phuket, Thái Lan) Đánh giá sâu 7 nhân tố hạ tầng, chất thải, an ninh đô thị. Tập trung vào du lịch biển đảo phát triển cao, quy mô mẫu lớn ($N=400$). Thiếu tính đặc thù về tài nguyên di sản và địa hình miền núi hiểm trở.
Nguyễn Thị Kim Ngân & CS (2022) (An Giang) Khung đánh giá 6 yếu tố kinh tế, xã hội, môi trường, hạ tầng. Chưa tích hợp yếu tố cảm nhận rủi ro và giá cả dịch vụ của du khách. Cần điều chỉnh và bổ sung biến quan sát phù hợp với vùng đồng bào thiểu số.
Mô hình đề xuất (Đoàn Hồng Nhung, 2023) Kết hợp nhận thức lợi ích, rủi ro, giá cả, hạ tầng và tài nguyên đặc thù. Giới hạn mẫu khảo sát thực địa ở mức $N=130$ do rào cản địa lý. Tối ưu, phản ánh chính xác thực trạng du lịch sinh thái - lịch sử Cao Bằng.

Ma trận ưu tiên yêu cầu nghiên cứu (MoSCoW)

  • Must have: Kiểm định độ tin cậy thang đo ($\text{Cronbach's } \alpha \ge 0.70$), kiểm định EFA ($\text{KMO} \ge 0.5$, $\text{Sig.} < 0.05$), kiểm định hồi quy OLS không vi phạm giả định.
  • Should have: Phân tích ANOVA theo nhóm nhân khẩu học (công việc, thu nhập, số lần tham quan).
  • Could have: Tích hợp các chỉ số phát triển bền vững theo bộ tiêu chuẩn GSTC.
  • Won't have: Mô hình hóa chuỗi thời gian dài hạn hoặc phương trình cấu trúc tuyến tính SEM bậc cao.

Thiết kế hệ thống biến và mô hình kinh tế lượng

[Nhận thức lợi ích (NTLI)]  (+) H1
[Nhận thức giá cả (NTGC)]   (+) H2
[Nhận thức rủi ro (NTRR)]   (-) H3 ---> [PHÁT TRIỂN DU LỊCH BỀN VỮNG (PTDLBV)]
[Cơ sở hạ tầng (CSHT)]     (+) H4
[Tài nguyên du lịch (TNDL)] (+) H5

Mô hình hồi quy tuyến tính tổng quát:

$$PTDLBV = \beta_0 + \beta_1 NTLI + \beta_2 NTGC + \beta_3 NTRR + \beta_4 CSHT + \beta_5 TNDL + \epsilon$$

Bảng mã hóa các biến quan sát trong thang đo

Nhóm nhân tố Mã hóa biến Nội dung biến quan sát thực nghiệm
Nhận thức lợi ích NTLI1–NTLI6 Cải thiện thu nhập dân cư, tạo việc làm, mở rộng quy mô kinh doanh, nâng cao năng lực cán bộ.
Nhận thức giá cả NTGC1–NTGC5 Xứng đáng với chi phí bỏ ra, giá trị thực sản phẩm, niêm yết công khai minh bạch.
Nhận thức rủi ro NTRR1–NTRR5 Rủi ro ép giá, tệ nạn xã hội trộm cắp, rủi ro đám đông chen lấn, an toàn giao thông đường đèo.
Cơ sở hạ tầng CSHT1–CSHT5 Trở ngại giao thông đèo dốc, dịch vụ vận tải công cộng, mạng lưới viễn thông, công nghệ quảng bá số.
Tài nguyên du lịch TNDL1–TNDL4 Cảnh quan thiên nhiên đặc sắc, lễ hội truyền thống, di tích lịch sử cách mạng, ẩm thực bản địa.
Du lịch bền vững PTDLBV1–PTDLBV3 Đánh giá tổng thể mức độ phát triển bền vững, tính đa dạng tài nguyên và bảo tồn văn hóa cộng đồng.

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

Nghiên cứu áp dụng quy trình 7 bước chuẩn mực theo Ary et al.:

  1. Xác định vấn đề: Khảo sát thực tế điểm nghẽn du lịch Cao Bằng sau đại dịch COVID-19.
  2. Tổng quan tài liệu: Rà soát 15 công trình quốc tế và trong nước về du lịch sinh thái và kinh tế phát triển.
  3. Thiết kế nghiên cứu: Xây dựng bảng hỏi cấu trúc 28 biến quan sát trên thang đo 5 mức Likert.
  4. Thu thập dữ liệu: Phát 140 phiếu điều tra trực tuyến kết hợp trực tiếp tại Pác Bó và Thác Bản Giốc; thu về 131 phiếu hợp lệ ($93,57%$), làm tròn phân tích 130 mẫu chuẩn hóa.
  5. Phân tích dữ liệu: Xử lý dữ liệu định lượng qua pipeline: Làm sạch dữ liệu $\to$ Thống kê mô tả $\to$ Cronbach's Alpha $\to$ EFA $\to$ Tương quan Pearson $\to$ Hồi quy OLS $\to$ ANOVA.
  6. Tổng hợp kết quả: Kiểm định các giả thiết $H_1, H_2, H_3, H_4, H_5$.
  7. Báo cáo & Khuyến nghị: Xây dựng lộ trình chính sách công cho địa phương.

Implementation và kết quả

Quy trình phát triển và xử lý dữ liệu

Quy trình tính toán và kiểm định kinh tế lượng được thực thi qua script phân tích tự động (tương thích Python statsmodels và cú pháp SPSS Syntax):

import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.stats.outliers_influence import variance_inflation_factor

# 1. Load và tiền xử lý dữ liệu khảo sát (N = 130)
df = pd.read_csv("caobang_tourism_survey.csv")
features = ['NTLI', 'NTGC', 'NTRR', 'CSHT', 'TNDL']
target = 'PTDLBV'

# 2. Xây dựng mô hình OLS đa biến
X = df[features]
X_const = sm.add_constant(X)
y = df[target]

model = sm.OLS(y, X_const).fit()

# 3. Kiểm tra hiện tượng đa cộng tuyến (VIF)
vif_data = pd.DataFrame()
vif_data["Feature"] = features
vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(len(features))]

print(model.summary())
print(vif_data)

Kiểm định và đánh giá mô hình (Testing & Validation)

1. Kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Tất cả các thang đo đều đạt chuẩn độ tin cậy ($\alpha > 0.70$) và hệ số tương quan biến - tổng ($Corrected\ Item\text{-}Total\ Correlation > 0.30$):

  • Thang đo Nhận thức lợi ích (NTLI): $\alpha = 0.842$ (6 biến quan sát)
  • Thang đo Nhận thức giá cả (NTGC): $\alpha = 0.789$ (5 biến quan sát)
  • Thang đo Nhận thức rủi ro (NTRR): $\alpha = 0.815$ (5 biến quan sát)
  • Thang đo Cơ sở hạ tầng (CSHT): $\alpha = 0.764$ (5 biến quan sát)
  • Thang đo Tài nguyên du lịch (TNDL): $\alpha = 0.828$ (4 biến quan sát)
  • Thang đo Phát triển DLBV (PTDLBV): $\alpha = 0.803$ (3 biến quan sát)

2. Phân tích nhân tố khám phá (EFA)

  • Hệ số $KMO = 0.812$ ($0.5 < KMO < 1.0$), kiểm định Bartlett có giá trị thống kê $\chi^2$ với mức ý nghĩa $p\text{-value} = 0.000 < 0.05$, chứng minh các biến quan sát có tương quan chặt chẽ.
  • Tổng phương sai trích đạt $63.45% > 50%$; hệ số tải nhân tố (Factor Loading) của tất cả các biến quan sát đều $> 0.55$, trích thành công 5 nhóm nhân tố độc lập qua phép xoay Varimax.

3. Kiểm định các giả định mô hình hồi quy

  • Kiểm định phân phối chuẩn phần dư: Biểu đồ Histogram có đường cong phân phối chuẩn dạng chuông, đồ thị Normal P-P Plot biểu diễn các điểm quan sát phân bố sát đường chéo $45^\circ$.
  • Kiểm định đa cộng tuyến: Hệ số phóng đại phương sai $VIF < 1.45$ (thấp hơn nhiều so với ngưỡng vi phạm $10.0$), khẳng định không xảy ra hiện tượng đa cộng tuyến giữa các biến độc lập.
  • Kiểm định tự tương quan: Hệ số Durbin-Watson đạt $1.892$ (nằm trong khoảng tối ưu từ $1.5$ đến $2.5$).

Kết quả nghiên cứu đạt được

+-----------------------------------------------------------------------+
|                KẾT QUẢ HỒI QUY CÁC NHÂN TỐ TÁC ĐỘNG                   |
+-------------------+----------------+-----------+-----------+----------+
| Biến độc lập      | Hệ số Beta (β) | t-stat    | Sig. (p)  | Kết luận |
+-------------------+----------------+-----------+-----------+----------+
| (Hằng số Const)   | 0.412          | 1.824     | 0.071     |          |
| NTLI              | 0.345          | 4.218     | 0.000     | Chấp nhận|
| TNDL              | 0.289          | 3.845     | 0.000     | Chấp nhận|
| CSHT              | 0.214          | 2.956     | 0.004     | Chấp nhận|
| NTGC              | 0.168          | 2.312     | 0.022     | Chấp nhận|
| NTRR              | -0.142         | -2.085    | 0.039     | Chấp nhận|
+-------------------+----------------+-----------+-----------+----------+
| R² = 0.582        | R² điều chỉnh = 0.565      | F = 34.52 (p = 0.000)|
+-------------------+----------------+-----------+-----------+----------+

Phương trình hồi quy chuẩn hóa:

$$PTDLBV = 0.345 \cdot NTLI + 0.289 \cdot TNDL + 0.214 \cdot CSHT + 0.168 \cdot NTGC - 0.142 \cdot NTRR$$

  • Nhận thức lợi ích ($\beta = 0.345$): Có tác động tích cực mạnh nhất; khi du khách nhận thấy hoạt động du lịch tạo sinh kế và phát triển cộng đồng, mức độ đánh giá phát triển bền vững tăng mạnh nhất.
  • Tài nguyên du lịch ($\beta = 0.289$): Đóng vai trò hạt nhân cốt lõi trong sức hấp dẫn điểm đến.
  • Cơ sở hạ tầng ($\beta = 0.214$): Giao thông kết nối và công nghệ số là điều kiện thúc đẩy bền vững.
  • Nhận thức rủi ro ($\beta = -0.142$): Tác động tiêu cực có ý nghĩa thống kê; rủi ro ép giá, mất trật tự công cộng và nguy hiểm giao thông đèo dốc làm suy giảm tính bền vững của ngành.
  • Kiểm định ANOVA: Có sự khác biệt có ý nghĩa thống kê ($p < 0.05$) về cảm nhận du lịch bền vững giữa các nhóm du khách theo tính chất công việc và tần suất tham quan.

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

  1. Đổi mới mô hình định lượng: Lần đầu tiên tích hợp cấu trúc đa chiều gồm cả yếu tố cản trở (Nhận thức rủi ro - NTRR) và yếu tố thúc đẩy (Nhận thức giá cả - NTGC, Hạ tầng số) trong một mô hình kinh tế lượng cho tỉnh miền núi biên giới Cao Bằng.
  2. Bằng chứng thực nghiệm chính xác: Cung cấp bộ chỉ số đo lường với $R^2 = 58,2%$, vượt trội về độ khớp thực tế so với các mô hình lý thuyết thuần túy trước đây.
  3. Đóng góp chính sách: Cung cấp luận cứ khoa học để tỉnh thực thi Nghị quyết Đại hội Đảng bộ tỉnh lần thứ XIX và Đề án Phát triển kết cấu hạ tầng giai đoạn 2022–2025 (Chương trình số 10-CTr/TU).
Tiêu chí Nghiên cứu Phuket (2013) Nghiên cứu An Giang (2022) Đề tài Cao Bằng (Nhung, 2023)
Phương pháp Định tính kết hợp khảo sát EFA + Hồi quy đa biến EFA + Hồi quy OLS + ANOVA đa tầng
Yếu tố Rủi ro & Giá Không đo lường riêng biệt Không tích hợp Đo lường định lượng cụ thể ($\beta = -0.142$)
Độ giải thích ($R^2$) Không xác định $51.4%$ $58.2%$ (Tối ưu hóa độ chính xác)

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

Tình huống ứng dụng thực tế

  • Khu du lịch sinh thái Thác Bản Giốc: Triển khai Hiệp định hợp tác bảo vệ và khai thác tài nguyên du lịch Thác Bản Giốc (Việt Nam) – Đức Thiên (Trung Quốc) theo mô hình "Du lịch xanh - Thông minh", kiểm soát chặt chẽ sức chứa tải sinh thái ($Carrying\ Capacity$).
  • Khu di tích Pác Bó: Số hóa hướng dẫn viên qua mã QR, nâng cấp hệ thống xe điện thân thiện môi trường, tái đầu tư gần 80 tỷ đồng vào bảo tồn cảnh quan tự nhiên Suối Lê-nin, Núi Các Mác.

Lộ trình triển khai chính sách (Implementation Roadmap)

[2023 - 2024] -----------------> [2024 - 2025] -----------------> [2025 - 2030]
* Chuẩn hóa giá & dịch vụ        * Đẩy nhanh Cao tốc             * Hoàn thiện hệ sinh thái
* Lắp đặt hệ thống QR Code         Đồng Đăng - Trà Lĩnh            du lịch thông minh
* Tập huấn nghiệp vụ người Tày   * Hoàn thiện hạ tầng viễn thông * Trung tâm du lịch Đông Bắc
  • Giai đoạn 2023–2024 (Kiểm soát rủi ro & Chuẩn hóa giá):
    • Niêm yết công khai giá hàng quán tại Pác Bó và Bản Giốc, dẹp bỏ triệt để hiện tượng ép giá ($NTRR \downarrow$).
    • Triển khai cổng thông tin du lịch thông minh, bản đồ số hóa CVĐC Non nước Cao Bằng.
  • Giai đoạn 2024–2025 (Đột phá hạ tầng):
    • Phối hợp đẩy nhanh tiến độ dự án Đường bộ cao tốc Đồng Đăng (Lạng Sơn) – Trà Lĩnh (Cao Bằng).
    • Nâng cấp 148 km đường tỉnh, cải thiện các cung đèo nguy hiểm (Đèo Khau Cốc Chà, Đèo Mã Phục).
  • Giai đoạn 2025–2030 (Phát triển bền vững toàn diện):
    • Nhân rộng mô hình du lịch cộng đồng homestay tại Làng đá Khuổi Ky (Trùng Khánh), bảo tồn kiến trúc nhà sàn truyền thống của người Tày.
    • Tỷ suất sinh lời xã hội ước tính: Tăng chi tiêu bình quân của du khách từ 1,2 triệu đồng lên 2,5 triệu đồng/chuyến, kéo dài thời gian lưu trú từ 1,3 ngày lên 2,5 ngày.

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

  • Hạn chế kỹ thuật: Quy mô mẫu khảo sát $N=130$ chủ yếu tập trung vào khách du lịch nội địa, phương thức phát phiếu trực tuyến chưa tiếp cận sâu nhóm du khách quốc tế cao cấp.
  • Ràng buộc ngoại cảnh: Địa hình miền núi chia cắt phức tạp gây khó khăn cho việc thu thập dữ liệu định kỳ theo chu kỳ mùa vụ.
  • Hướng phát triển tương lai:
    • Mở rộng cỡ mẫu lên $N \ge 500$, ứng dụng mô hình cấu trúc tuyến tính bình phương bé nhất từng phần (PLS-SEM) để phân tích các biến trung gian điều tiết.
    • Ứng dụng phân tích không gian địa lý (GIS) kết hợp dữ liệu lớn từ viễn thông để dự báo mật độ du khách thời gian thực tại các danh thắng.

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

  • Cơ quan quản lý Nhà nước (Sở VH,TT&DL Cao Bằng): Bộ công cụ định lượng hỗ trợ thẩm định dự án đầu tư và hoạch định chính sách du lịch xanh.
  • Doanh nghiệp & Hộ kinh doanh lữ hành: Căn cứ xác định chiến lược định giá, tiêu chuẩn hóa chất lượng dịch vụ và quản trị rủi ro trải nghiệm khách hàng.
  • Cộng đồng cư dân địa phương: Tăng cơ hội việc làm, nâng cao thu nhập bình quân từ chuỗi giá trị du lịch nông nghiệp - làng nghề truyền thống (làng rèn Phúc Sen, hương Phja Thắp).
  • Sinh viên & Giới nghiên cứu: Tài liệu tham khảo ứng dụng kinh tế lượng thực chứng trong phân tích kinh tế du lịch và phát triển bền vững.

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

  1. Mô hình nghiên cứu sử dụng công cụ thống kê nào để xử lý dữ liệu?
    Dữ liệu được làm sạch trên Microsoft Excel 2021 và phân tích chuyên sâu trên phần mềm IBM SPSS Statistics phiên bản 26.0 (thực hiện Cronbach's Alpha, EFA, OLS Regression, ANOVA).

  2. Tại sao nhân tố Nhận thức rủi ro (NTRR) lại có hệ số âm trong mô hình?
    Hệ số $\beta = -0.142$ ($p = 0.039 < 0.05$) phản ánh quy luật kinh tế: khi du khách cảm nhận rủi ro về an toàn giao thông đường đèo, tệ nạn xã hội hoặc ép giá gia tăng, mức độ đánh giá về tính bền vững của điểm đến sẽ suy giảm tương ứng.

  3. Cỡ mẫu 130 quan sát có đảm bảo tính đại diện thống kê không?
    Theo Hair et al. (2014), tỷ lệ mẫu trên một biến quan sát tối thiểu là 5:1. Mô hình có 25 biến quan sát trong thang đo độc lập, cỡ mẫu 130 đạt tỷ lệ $> 5:1$, thỏa mãn hoàn toàn điều kiện kiểm định EFA và hồi quy đa biến.

  4. Dự án giải quyết nút thắt hạ tầng giao thông Cao Bằng như thế nào?
    Đề tài kiến nghị ưu tiên nguồn lực cho dự án cao tốc Đồng Đăng – Trà Lĩnh, cải tạo 6 tuyến quốc lộ huyết mạch (QL3, QL34, QL4A) và phát triển mạng lưới vận tải xe buýt công cộng liên huyện.

  5. Làm sao để người dân địa phương tham gia trực tiếp vào chuỗi giá trị du lịch bền vững?
    Thông qua việc thành lập các Hợp tác xã du lịch cộng đồng, đào tạo kỹ năng đón tiếp du khách, bảo tồn nghề rèn - dệt thổ cẩm bản địa, và phân phối lại doanh thu từ vé tham quan cho công tác bảo tồn di sản.

Kết luận

Khóa luận tốt nghiệp đã hoàn thành xuất sắc mục tiêu nghiên cứu: xác lập cơ sở khoa học và lượng hóa thành công các nhân tố ảnh hưởng đến phát triển du lịch bền vững tại tỉnh Cao Bằng với độ tin cậy $R^2 = 58,2%$. Kết quả khẳng định phát triển du lịch bền vững tại Cao Bằng phải dựa trên kiềng 3 chân: Phát huy lợi ích kinh tế cộng đồng ($NTLI$), Bảo tồn nguyên vẹn tài nguyên thiên nhiên - lịch sử ($TNDL$), và Đột phá nâng cấp kết cấu hạ tầng giao thông ($CSHT$). Công trình cung cấp hệ thống giải pháp khả thi, đóng góp thiết thực vào mục tiêu đưa Cao Bằng trở thành trung tâm du lịch xanh của vùng Đông Bắc Việt Nam vào năm 2030.