Giới thiệu dự án

Du lịch văn hóa đang khẳng định vị thế là một trong những trụ cột chiến lược hàng đầu trong định hướng phát triển kinh tế bền vững của Việt Nam. Theo tinh thần Nghị quyết số 08-NQ/TW của Bộ Chính trị về việc phát triển du lịch trở thành ngành kinh tế mũi nhọn và Chiến lược phát triển du lịch Việt Nam đến năm 2030, việc khai thác hiệu quả các giá trị di sản văn hóa không chỉ tạo ra động lực kinh tế đa ngành mà còn giải quyết việc làm và nâng cao sinh kế cho cộng đồng địa phương. Tỉnh Thái Bình sở hữu hệ thống tài nguyên du lịch văn hóa phong phú với 2.969 di tích lịch sử - văn hóa, bao gồm 2 di tích quốc gia đặc biệt (chùa Keo và Khu lăng mộ, đền thờ các vị vua triều Trần), 114 di tích cấp quốc gia và 550 di tích cấp tỉnh.

Tuy nhiên, thực trạng phát triển du lịch văn hóa tại Thái Bình đang đối mặt với bài toán bất cân xứng giữa tiềm năng và hiệu quả khai thác:

  • Hoạt động du lịch tại các điểm đến phần lớn mang tính tự phát, thiếu quy hoạch liên kết đồng bộ.
  • Số ngày lưu trú trung bình của du khách thấp, doanh thu từ các dịch vụ bổ trợ chưa cao.
  • Thiếu các nghiên cứu định lượng chuyên sâu đo lường chính xác các nhân tố chi phối mức độ thỏa mãn của du khách để làm căn cứ hoạch định chính sách.

Đề tài "Những yếu tố tác động đến mức độ hài lòng của du khách đối với du lịch văn hóa tỉnh Thái Bình" (thực hiện bởi sinh viên Nguyễn Thị Thu Hà, GVHD: TS. Nguyễn Thùy Anh, Trường Đại học Kinh tế - ĐHQGHN) được xây dựng nhằm lượng hóa trải nghiệm du khách thông qua các mô hình kinh tế lượng thực nghiệm.

Mục tiêu dự án:
1. Hệ thống hóa cơ sở lý thuyết về sự hài lòng và chất lượng dịch vụ du lịch văn hóa.
2. Xây dựng khung đo lường gồm 5 nhân tố độc lập (25 biến quan sát) và 1 nhân tố phụ thuộc.
3. Thu thập và xử lý tập dữ liệu thực nghiệm (N = 200) bằng phương pháp kiểm định độ tin cậy và hồi quy đa biến.
4. Đo lường mức độ tác động cụ thể của từng yếu tố và đề xuất giải pháp khả thi cho các bên liên quan.

Phạm vi nghiên cứu tập trung vào du khách nội địa đã trải nghiệm tại 5 điểm đến văn hóa tiêu biểu của tỉnh: Chùa Keo, Đền Trần, Đền Tiên La, Đền Đồng Xâm (kèm làng nghề chạm bạc Đồng Xâm) và Nhà thờ Chính Tòa Thái Bình.


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

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

Nghiên cứu so sánh các phương pháp tiếp cận đánh giá chất lượng dịch vụ và sự hài lòng trong ngành du lịch hiện nay:

Tiêu chí Khảo sát định tính truyền thống Mô hình SERVQUAL nguyên bản Khung nghiên cứu kinh tế lượng tích hợp
Độ chính xác đo lường Thấp, mang tính chủ quan Trung bình, tập trung khoảng cách kỳ vọng Cao, lượng hóa hệ số tác động $\beta$ cụ thể
Tính thích ứng đặc thù Linh hoạt nhưng thiếu chuẩn hóa Khó đo lường giá trị di sản phi vật thể Tối ưu riêng cho tài nguyên du lịch văn hóa
Khả năng kiểm định Không có kiểm định thống kê Kiểm định độ tin cậy cơ bản Kiểm định EFA, Pearson, OLS và ANOVA
Giá trị ứng dụng chính sách Mang tính định hướng chung Nhận diện khoảng cách dịch vụ Cung cấp thứ tự ưu tiên phân bổ ngân sách

Để tối ưu hóa bảng hỏi và chỉ số đo lường, các biến số được phân loại theo ma trận MoSCoW:

  • Must-have (Bắt buộc): Thang đo tài nguyên du lịch văn hóa (TN), Cơ sở hạ tầng (HT), Cơ sở vật chất kỹ thuật (KT), Giá cả dịch vụ (GC), An ninh trật tự (AN).
  • Should-have (Nên có): Thống kê nhân khẩu học đa chiều (độ tuổi, thu nhập, học vấn, nghề nghiệp), nguồn thông tin tiếp cận.
  • Could-have (Có thể mở rộng): Đo lường lòng trung thành và ý định quay trở lại (HL2, HL3).
  • Won't-have (Tạm thời loại trừ): Phân tích chi tiết hành vi du khách quốc tế (do giới hạn mẫu trong giai đoạn nghiên cứu).

Thiết kế hệ thống đo lường

Khung lý thuyết của nghiên cứu được xây dựng trên sự kế thừa và tinh chỉnh từ mô hình chất lượng kỹ thuật/chức năng của Grönroos (1984), Chỉ số hài lòng khách hàng Mỹ (ACSI), Chỉ số Châu Âu (ECSI) và mô hình du lịch nông thôn của Phan Việt Đua và cộng sự (2022).

Tech Stack và công cụ phân tích

  • Phần mềm thống kê chính: IBM SPSS Statistics v22.0.
  • Công cụ thu thập dữ liệu: Google Forms API kết hợp bảng hỏi trực tuyến đa kênh.
  • Môi trường xác thực thuật toán bổ trợ: Python v3.10 (pandas v2.0.3, statsmodels v0.14.0, scikit-learn v1.3.0).
  • Cấu trúc thang đo (Data Schema): 25 biến quan sát độc lập và 3 biến quan sát phụ thuộc được mã hóa theo thang đo khoảng Likert 5 mức độ (1: Hoàn toàn không đồng ý $\rightarrow$ 5: Hoàn toàn đồng ý).
Ma trận biến quan sát (Data Schema):
- TN = {TN1: Giá trị độc đáo; TN2: Công tác bảo tồn; TN3: Hàng lưu niệm đặc trưng}
- HT = {HT1: Đường sá rộng rãi; HT2: Mặt đường tốt; HT3: Bãi đỗ xe; HT4: Nhà vệ sinh}
- KT = {KT1: Ăn uống tiện nghi; KT2: Ăn uống sạch sẽ; KT3: Lưu trú tiện nghi; KT4: Lưu trú sạch sẽ; KT5: Cửa hàng lưu niệm; KT6: Điểm vui chơi}
- GC = {GC1: Giá vé vào; GC2: Giá quà lưu niệm; GC3: Giá lưu trú; GC4: Giá ăn uống}
- AN = {AN1: Không trộm cắp; AN2: Không chèo kéo; AN3: Không ăn xin; AN4: Không thách giá; AN5: Vệ sinh môi trường}
- HL = {HL1: Hài lòng chung; HL2: Dự định quay lại; HL3: Sẵn sàng giới thiệu}

Phương pháp nghiên cứu

Nghiên cứu áp dụng quy trình 6 giai đoạn nghiêm ngặt:

  1. Xác định bài toán & Khảo sát sơ bộ: Điều tra nhận diện sơ bộ ($N_1 = 50$) trên các cộng đồng du lịch.
  2. Quy mô mẫu định lượng chính thức: Áp dụng quy tắc kinh nghiệm của Hair et al. (tỷ lệ $5:1$ trên 23-25 biến quan sát $\rightarrow N \ge 115$) và quy tắc hồi quy bội $N \ge 50 + 8p = 50 + 8 \times 5 = 90$. Cỡ mẫu chính thức đạt $N = 200$, đảm bảo tính đại diện và độ tin cậy thống kê.
  3. Thu thập dữ liệu thực địa: Thực hiện từ 20/03/2023 đến 15/04/2023.
  4. Kiểm định độ tin cậy thang đo (Cronbach's Alpha): Tiêu chuẩn $\alpha \ge 0.70$ (chấp nhận mức $\ge 0.60$ cho nghiên cứu khám phá) và hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $> 0.30$.
  5. Phân tích nhân tố khám phá (EFA): Tiêu chuẩn $0.5 \le \text{KMO} \le 1.0$, kiểm định Bartlett có $\text{Sig.} < 0.05$, hệ số Eigenvalue $> 1.0$, tổng phương sai trích (Total Variance Explained) $> 50%$, hệ số tải nhân tố (Factor Loading) $> 0.50$.
  6. Hồi quy tuyến tính OLS & One-way ANOVA: Đánh giá $R^2$, kiểm định đa cộng tuyến qua hệ số phóng đại phương sai $\text{VIF} < 2.0$, kiểm định sự khác biệt giữa các nhóm nhân khẩu học.

Implementation và kết quả

Quy trình phân tích dữ liệu

Thuật toán phân tích kinh tế lượng được thiết lập theo quy trình toán học sau:

  1. Công thức tính hệ số tin cậy Cronbach's Alpha: $$\alpha = \frac{K}{K - 1} \left(1 - \frac{\sum_{i=1}^{K} \sigma_{Y_i}^2}{\sigma_X^2}\right)$$ Trong đó $K$ là số lượng biến quan sát trong thang đo, $\sigma_{Y_i}^2$ là phương sai của biến quan sát thứ $i$, và $\sigma_X^2$ là phương sai của tổng điểm các biến.

  2. Phương trình hồi quy tuyến tính bội OLS: $$HL = \beta_0 + \beta_1 TN + \beta_2 HT + \beta_3 KT + \beta_4 GC + \beta_5 AN + \epsilon$$

Đoạn mã Python mô phỏng quy trình kiểm định thang đo và hồi quy OLS tương đương trên SPSS v22.0:

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

# 1. Tính toán Cronbach's Alpha cho từng thang đo
def cronbach_alpha(df_items: pd.DataFrame) -> float:
    k = df_items.shape[1]
    item_variances = df_items.var(axis=0, ddof=1).sum()
    total_variance = df_items.sum(axis=1).var(ddof=1)
    return (k / (k - 1)) * (1 - (item_variances / total_variance))

# 2. Xử lý mô hình hồi quy tuyến tính đa biến OLS
def run_regression_pipeline(X: pd.DataFrame, y: pd.Series):
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    
    # Tính toán hệ số phóng đại phương sai (VIF)
    vif_data = pd.DataFrame()
    vif_data["Feature"] = X.columns
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(X.shape[1])]
    
    return model.summary(), vif_data

Kiểm định và xác thực mô hình

1. Kết quả kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Tất cả 6 thang đo thành phần đều vượt qua ngưỡng kiểm định tiêu chuẩn:

  • Tài nguyên du lịch văn hóa (TN): $\alpha = 0.812$ (3 biến quan sát, Item-Total Correlation $> 0.58$).
  • Cơ sở hạ tầng (HT): $\alpha = 0.795$ (4 biến quan sát, Item-Total Correlation $> 0.52$).
  • Cơ sở vật chất kỹ thuật (KT): $\alpha = 0.841$ (6 biến quan sát, Item-Total Correlation $> 0.55$).
  • Giá cả dịch vụ (GC): $\alpha = 0.788$ (4 biến quan sát, Item-Total Correlation $> 0.49$).
  • An ninh trật tự, an toàn (AN): $\alpha = 0.824$ (5 biến quan sát, Item-Total Correlation $> 0.53$).
  • Sự hài lòng của du khách (HL): $\alpha = 0.856$ (3 biến quan sát, Item-Total Correlation $> 0.67$).

2. Kết quả phân tích nhân tố khám phá (EFA)

  • Hệ số $\text{KMO} = 0.832$ (thỏa mãn $0.5 < \text{KMO} < 1.0$), kiểm định Bartlett's Test of Sphericity có mức ý nghĩa $\text{Sig.} = 0.000 < 0.05$, chứng minh ma trận tương quan giữa các biến là hoàn toàn phù hợp.
  • Trích xuất 5 nhân tố độc lập tại giá trị $\text{Eigenvalue} > 1.0$, tổng phương sai trích đạt $62.48% > 50%$, chứng minh 5 nhân tố giải thích được $62.48%$ độ biến thiên của dữ liệu.

3. Kết quả hồi quy tuyến tính và phân tích ANOVA

Mô hình hồi quy đạt mức ý nghĩa thống kê cao với $R^2 = 0.584$ và $R^2 \text{ hiệu chỉnh} = 0.573$ ($p < 0.001$). Hệ số $\text{VIF}$ của tất cả các biến đều nằm trong khoảng $1.15 - 1.48 < 2.0$, chứng minh mô hình không vi phạm hiện tượng đa cộng tuyến.

Hệ số hồi quy chuẩn hóa (Standardized Beta):
HL = 0.342*TN + 0.228*AN + 0.195*GC + 0.164*HT + 0.131*KT
(Tất cả các hệ số đều có giá trị p-value < 0.05)

Phân tích One-way ANOVA cho thấy không có sự khác biệt có ý nghĩa thống kê về mức độ hài lòng theo giới tính ($p > 0.05$), nhưng có sự khác biệt đáng kể theo nhóm độ tuổi ($p = 0.021$) và mức thu nhập ($p = 0.034$), trong đó nhóm du khách trung niên và có thu nhập ổn định đánh giá cao hơn về giá trị tài nguyên tâm linh.


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

Nghiên cứu mang lại những đóng góp rõ rệt về cả phương diện học thuật và thực tiễn:

Khía cạnh Các nghiên cứu trước (Gronroos 1984; Phan Thị Dạng 2015) Đề tài nghiên cứu du lịch văn hóa Thái Bình (2023)
Mô hình hóa dữ liệu Đơn biến hoặc khảo sát du lịch sinh thái chung Thiết kế thang đo đặc thù kết hợp văn hóa vật thể, phi vật thể và làng nghề
Độ nhạy của nhân tố Tập trung nhiều vào hạ tầng và cơ sở lưu trú Chứng minh Tài nguyên văn hóa ($\beta = 0.342$) và An ninh trật tự ($\beta = 0.228$) là 2 yếu tố chi phối mạnh nhất
Tối ưu hóa bảng hỏi Bảng hỏi dài, dễ gây nhiễu và sai lệch đo lường Rút gọn chuẩn hóa 25 biến quan sát đạt tổng phương sai trích $62.48%$
Độ tin cậy kỹ thuật Ít kiểm định giả định hồi quy sâu Thực hiện đầy đủ kiểm định đa cộng tuyến (VIF), phần dư chuẩn và One-way ANOVA

Đóng góp thực tiễn cho ngành du lịch địa phương

  • Số liệu định lượng chính xác: Cung cấp bộ chỉ số đo lường thực nghiệm đầu tiên về mức độ hài lòng du lịch văn hóa tại Thái Bình cho Sở Văn hóa, Thể thao và Du lịch.
  • Thứ tự ưu tiên can thiệp: Chỉ ra rằng việc đầu tư bảo tồn tính nguyên bản của di tích mang lại tỷ suất cải thiện sự hài lòng cao gấp 2.6 lần so với việc đơn thuần mở rộng cơ sở vật chất giải trí thương mại.

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

Dựa trên kết quả định lượng, lộ trình can thiệp và triển khai giải pháp được chia thành các trụ cột chiến lược:

1. Nâng cao sức hấp dẫn tài nguyên văn hóa ($\beta = 0.342$)

  • Triển khai gắn mã QR Code tra cứu thông tin lịch sử và thuyết minh tự động (Audio Guide) tại Chùa Keo và Đền Trần.
  • Xây dựng không gian trải nghiệm thực tế quy trình chế tác thủ công tại làng nghề chạm bạc Đồng Xâm thay vì chỉ bày bán thụ động.

2. Tăng cường an ninh trật tự và vệ sinh môi trường ($\beta = 0.228$)

  • Thiết lập đường dây nóng phản ánh tình trạng chèo kéo, ăn xin, nâng giá dịch vụ tại các mùa lễ hội cao điểm.
  • Bố trí các cụm thùng rác phân loại và khu vệ sinh công cộng đạt chuẩn du lịch tại các điểm di tích trọng điểm.

3. Minh bạch hóa giá cả dịch vụ ($\beta = 0.195$)

  • Bắt buộc niêm yết giá công khai đối với tất cả các dịch vụ ăn uống, gửi xe và quà lưu niệm quanh khuôn viên di tích.

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

Hạn chế kỹ thuật

  1. Quy mô và phương pháp chọn mẫu: Mẫu $N = 200$ thu thập chủ yếu qua hình thức trực tuyến và tập trung vào du khách nội địa, chưa bao quát sâu nhóm du khách quốc tế và các đoàn caravan tự túc.
  2. Phạm vi địa lý: Mới khảo sát tại 5 điểm đến tiêu biểu, chưa mở rộng toàn diện sang các quần thể di tích cấp huyện khác như đền A Sào, chùa La Vân.
  3. Mô hình phân tích: Sử dụng hồi quy OLS tuyến tính tĩnh, chưa phân tích các mối quan hệ đa tầng (cấu trúc tuyến tính SEM/PLS-SEM) giữa chất lượng dịch vụ, hình ảnh điểm đến và lòng trung thành.

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

  • Mở rộng quy mô mẫu lên $N \ge 500$, ứng dụng mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) để phân tích tác động điều tiết (moderating effects) của nhân khẩu học.
  • Tích hợp công nghệ phân tích cảm xúc dữ liệu lớn (Big Data Sentiment Analysis) từ các đánh giá trực tuyến trên Google Maps và mạng xã hội du lịch để cập nhật chỉ số CSI theo thời gian thực.

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

+-------------------------------------------------------------------------------+
|                             ĐỐI TƯỢNG HƯỞNG LỢI                               |
+-------------------+-----------------------------------------------------------+
| Nhóm đối tượng    | Giá trị và Lợi ích định lượng                             |
+-------------------+-----------------------------------------------------------+
| Sinh viên /       | Kế thừa khung lý thuyết, bộ câu hỏi 25 biến quan sát và   |
| Học viên          | quy trình phân tích SPSS chuẩn mực làm tài liệu mẫu.      |
+-------------------+-----------------------------------------------------------+
| Cơ quan quản lý   | Căn cứ dữ liệu thực chứng với R2 = 0.584 để phân bổ ngân |
| nhà nước (Sở VHTTDL)| sách phát triển hạ tầng và siết chặt an ninh trật tự.    |
+-------------------+-----------------------------------------------------------+
| Doanh nghiệp      | Định hình các gói tour văn hóa đặc thù, định giá dịch vụ  |
| lữ hành           | hợp lý và liên kết chuỗi cung ứng sản phẩm lưu niệm.      |
+-------------------+-----------------------------------------------------------+
| Cộng đồng địa     | Nâng cao thu nhập nhờ phát triển dịch vụ phụ trợ, gìn giữ |
| phương            | và phát huy giá trị di sản làng nghề truyền thống.        |
+-------------------+-----------------------------------------------------------+

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

1. Yêu cầu dữ liệu tối thiểu để triển khai mô hình hồi quy nghiên cứu sự hài lòng là gì?

Cần tối thiểu 5 mẫu cho mỗi biến quan sát trong phân tích nhân tố EFA (Hair et al., 2006) và tối thiểu $N \ge 50 + 8p$ đối với hồi quy bội (với $p$ là số biến độc lập). Với mô hình 5 nhân tố (25 biến quan sát), quy mô mẫu từ $N = 125$ đến $N = 200$ là tối ưu để đảm bảo độ tin cậy thống kê.

2. Tại sao thang đo Tài nguyên du lịch văn hóa lại có tác động lớn nhất đến sự hài lòng?

Hệ số hồi quy $\beta = 0.342$ cho thấy đây là động lực cốt lõi thu hút du khách đến với Thái Bình. Khách du lịch văn hóa tìm kiếm trải nghiệm chiều sâu tâm linh và lịch sử; do đó, giá trị kiến trúc nguyên bản và không gian di sản đóng vai trò quyết định đến cảm nhận tổng thể.

3. Làm thế nào để kiểm soát hiện tượng đa cộng tuyến trong mô hình?

Cần theo dõi hệ số phóng đại phương sai (VIF). Nếu $\text{VIF} < 2.0$ (hoặc nghiêm ngặt là $< 5.0$), mô hình hoàn toàn không vi phạm hiện tượng đa cộng tuyến giữa các biến độc lập.

4. Chi phí đầu tư để số hóa các điểm đến di tích văn hóa là bao nhiêu?

Chi phí triển khai mã hóa QR, số hóa 3D và hệ thống thuyết minh tự động cơ bản tại 5 điểm đến trọng điểm ước tính dao động từ 150 - 300 triệu VNĐ, mang lại hiệu quả nâng cao trải nghiệm tức thì cho du khách.

5. Kết quả nghiên cứu này có thể áp dụng cho các tỉnh lân cận không?

Bộ công cụ đo lường và thang đo có thể chuyển giao và ứng dụng tốt cho các địa phương có đặc thù tài nguyên văn hóa vùng châu thổ sông Hồng tương đồng như Nam Định, Hà Nam, Hưng Yên.


Kết luận

Nghiên cứu "Những yếu tố tác động đến mức độ hài lòng của du khách đối với du lịch văn hóa tỉnh Thái Bình" đã giải quyết thành công bài toán lượng hóa trải nghiệm du khách thông qua quy trình phân tích kinh tế lượng chuẩn mực. Với hệ số xác định $R^2 = 0.584$, nghiên cứu khẳng định sự hài lòng của du khách chịu sự chi phối đồng thời của 5 nhân tố với thứ tự ưu tiên: Tài nguyên du lịch văn hóa ($\beta = 0.342$), An ninh trật tự ($\beta = 0.228$), Giá cả dịch vụ ($\beta = 0.195$), Cơ sở hạ tầng ($\beta = 0.164$) và Cơ sở vật chất kỹ thuật ($\beta = 0.131$).

Kết quả này cung cấp cơ sở dữ liệu thực chứng có giá trị cao, giúp các nhà quản lý du lịch và chính quyền địa phương chuyển dịch từ tư duy quản lý trực giác sang hoạch định chính sách dựa trên dữ liệu định lượng, tạo nền tảng vững chắc để đưa du lịch văn hóa Thái Bình bứt phá trong giai đoạn 2023 - 2030.