Tổng quan nghiên cứu

Trong giai đoạn 2014 - 2018, tỉnh Bình Dương ghi nhận sự bùng phát dữ dội của bệnh sốt xuất huyết Dengue với tổng số hơn 29.400 ca mắc và 22 trường hợp tử vong. Đáng chú ý, năm 2017 ghi nhận đỉnh dịch với 8.507 ca mắc (tăng 114,6% so với năm 2016), tiếp theo là năm 2018 với 8.244 ca. Sự gia tăng dân số cơ học nhanh chóng tại các khu công nghiệp cùng đặc điểm khí hậu nhiệt đới gió mùa đã khiến dịch bệnh diễn biến phức tạp, khó kiểm soát theo các chu kỳ truyền thống. Vấn đề lớn đối với ngành y tế dự phòng địa phương là hệ thống giám sát chủ yếu dựa trên báo cáo thụ động, thiếu hụt các công cụ phân tích không gian - thời gian trực quan để hỗ trợ cảnh báo sớm và điều phối nguồn lực kịp thời.

Luận văn thạc sĩ chuyên ngành Hệ thống thông tin với đề tài "Biểu diễn và phân tích trực quan dữ liệu dịch bệnh - Áp dụng tập dữ liệu dịch bệnh sốt xuất huyết tỉnh Bình Dương" được thực hiện nhằm giải quyết bài toán cấp thiết trên. Mục tiêu nghiên cứu tập trung vào việc xử lý, đồng bộ hóa và mô hình hóa trực quan tập dữ liệu dịch tễ cùng các thông số khí tượng tại 9 huyện, thị xã, thành phố của tỉnh Bình Dương trong khoảng thời gian 5 năm (từ ngày 01/01/2014 đến ngày 31/12/2018). Nghiên cứu hướng tới việc xây dựng hệ sinh thái biểu đồ đa chiều hỗ trợ chuyên gia dịch tễ rút ngắn khoảng 40% thời gian phân tích số liệu thô, nâng cao khả năng nhận diện ổ dịch trước 2 đến 4 tuần, từ đó tối ưu hóa công tác phân bổ hóa chất dập dịch và giảm thiểu tổn thất cho địa phương có mật độ dân số 769 người/km2.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên ba trụ cột lý thuyết chính trong khoa học dữ liệu và dịch tễ học:

Thứ nhất, lý thuyết trực quan hóa dữ liệu (Data Visualization Theory). Lý thuyết này định nghĩa trực quan hóa là kỹ thuật chuyển đổi dữ liệu thuộc tính thành các đối tượng đồ họa như điểm, đường, thanh và bề mặt biểu diễn. Mục tiêu trọng tâm là khai thác tối đa năng lực xử lý của thị giác và não bộ con người để trích xuất tri thức ẩn, hỗ trợ phương pháp nhìn - hiểu thay vì đọc các bảng số liệu phức tạp.

Thứ hai, mô hình biểu diễn khối không gian - thời gian (Spatio-Temporal Data Cube). Không gian nghiên cứu được ánh xạ trên hệ tọa độ đa chiều: trục mặt phẳng biểu diễn vị trí địa lý thực tế (9 đơn vị hành chính cấp huyện), trục tung biểu diễn chuỗi thời gian (chu kỳ 60 tháng) và các trục giá trị thể hiện số lượng ca mắc cùng các chỉ số môi trường.

Thứ ba, mô hình hồi quy tuyến tính đa biến (Multivariable Linear Regression). Khung toán học thiết lập quan hệ phụ thuộc giữa biến số ca mắc sốt xuất huyết với các biến giải thích độc lập bao gồm nhiệt độ không khí, độ ẩm tương đối và lượng mưa.

Các khái niệm cốt lõi được chuẩn hóa trong nghiên cứu bao gồm: ca mắc sốt xuất huyết (ký hiệu p, đơn vị tính: người), nhiệt độ không khí trung bình (ký hiệu f, đơn vị tính: độ C), độ ẩm không khí tương đối (ký hiệu h, đơn vị tính: phần trăm), lượng mưa tích lũy (ký hiệu e, đơn vị tính: mm) và vector truyền bệnh muỗi Aedes aegypti.

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

Nguồn dữ liệu: Nghiên cứu thu thập và tích hợp hai nguồn dữ liệu thực tế chính thức: số liệu ca bệnh sốt xuất huyết hằng tháng từ Trung tâm Y tế Dự phòng tỉnh Bình Dương và chuỗi dữ liệu khí tượng (nhiệt độ trung bình từ 26°C đến 27°C, độ ẩm từ 66% đến 86%, lượng mưa hàng tháng) từ Trung tâm Quan trắc và Kỹ thuật Tài nguyên Môi trường tỉnh Bình Dương trong giai đoạn từ ngày 01/01/2014 đến ngày 31/12/2018.

Cỡ mẫu và phương pháp chọn mẫu: Cỡ mẫu nghiên cứu gồm toàn bộ 60 quan sát theo chuỗi thời gian 12 tháng liên tục trong 5 năm, bao phủ 100% phạm vi không gian của 9 đơn vị hành chính (thành phố Thủ Dầu Một, thị xã Thuận An, thị xã Dĩ An, thị xã Bến Cát, thị xã Tân Uyên, huyện Bàu Bàng, huyện Bắc Tân Uyên, huyện Dầu Tiếng và huyện Phú Giáo). Phương pháp chọn mẫu toàn bộ theo chuỗi thời gian hồi cứu được áp dụng nhằm loại trừ sai số chọn mẫu cục bộ, đảm bảo tính đại diện và tính liên tục tuyệt đối của dữ liệu dịch tễ.

Phương pháp phân tích: Tác giả phát triển chương trình phân tích trên ngôn ngữ lập trình Python, sử dụng module Pyplot thuộc thư viện Matplotlib kết hợp giao diện tương tác Tkinter để thiết kế bảng điều khiển và dựng biểu đồ 3 chiều đa trục. Đồng thời, nghiên cứu ứng dụng phần mềm SPSS 20 để kiểm định thống kê hồi quy đa biến, phân tích phương sai ANOVA, đánh giá hệ số tự tương quan Durbin-Watson (đạt mức chuẩn 2.294), kiểm tra hiện tượng đa cộng tuyến VIF và đánh giá phân phối phần dư chuẩn. Việc kết hợp lập trình Python và phân tích SPSS được lựa chọn vì giúp kết hợp chặt chẽ giữa kiểm định định lượng toán học chính xác với khả năng biểu diễn trực quan tương tác cao cho người dùng y tế.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, tính quy luật mùa vụ rõ nét và tác động của lượng mưa: Dịch sốt xuất huyết tại Bình Dương bùng phát mạnh mẽ nhất vào các tháng 7, 8, 9 và 10 hằng năm, trùng khớp với thời kỳ cao điểm của mùa mưa. Trong năm 2014, tháng 7 ghi nhận lượng mưa cao nhất trong năm và số ca bệnh tăng vọt, trong khi tháng 2 có lượng mưa thấp nhất tương ứng với số ca nhiễm thấp nhất. Độ ẩm không khí trung bình trong năm dao động từ 76% đến 80%, nhưng tại các tháng đỉnh dịch, độ ẩm đạt ngưỡng cao nhất tới 86% (tháng 9), tạo điều kiện sinh sản tối ưu cho vector truyền bệnh.

Thứ hai, sự phân hóa dịch tễ sâu sắc theo không gian và tốc độ đô thị hóa: Dịch bệnh tập trung cao độ tại các khu vực đô thị phía Nam có mật độ dân cư đông đúc và tập trung nhiều khu công nghiệp. Điển hình trong đợt dịch lớn năm 2017 với 8.507 ca mắc, ba địa bàn thành phố Thủ Dầu Một, thị xã Thuận An và thị xã Dĩ An chiếm trên 65% tổng số ca mắc của toàn tỉnh. Ngược lại, các huyện nông nghiệp phía Bắc như Phú Giáo, Bắc Tân Uyên và Dầu Tiếng duy trì số ca mắc ở mức thấp ổn định.

Thứ ba, sự xuất hiện của hai đợt dịch bùng phát nghiêm trọng: Chuỗi số liệu 5 năm ghi nhận hai đỉnh dịch lớn vào năm 2015 với 5.991 ca mắc (14 ca tử vong, tăng 120,7% so với năm 2014) và năm 2017 với 8.507 ca mắc (3 ca tử vong, tăng 114,6% so với năm 2016).

Thứ tư, kết quả kiểm định mô hình định lượng: Phân tích tương quan Pearson trên SPSS 20 xác nhận biến số ca mắc sốt xuất huyết có mối tương quan tuyến tính có ý nghĩa thống kê cao với nhiệt độ, độ ẩm và lượng mưa với giá trị mức ý nghĩa kiểm định sig đều bằng 0.000. Mô hình hồi quy đạt chuẩn kiểm định chuỗi với giá trị Durbin-Watson bằng 2.294, khẳng định mô hình không bị vi phạm hiện tượng tự tương quan chuỗi bậc nhất.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự gia tăng đột biến ca bệnh vào các tháng 7, 8, 9 và 10 là do lượng mưa dầm tạo ra vô số vũng nước đọng và phế thải chứa nước quanh các khu nhà trọ công nhân. Kết hợp với nền nhiệt độ trung bình thuận lợi từ 26°C đến 27°C, chu kỳ sinh trưởng của ấu trùng muỗi bị rút ngắn đáng kể, làm mật độ vector truyền bệnh tăng nhanh. Tốc độ tăng dân số cơ học nhanh tại Thuận An và Dĩ An tạo điều kiện cho virus lây lan theo cấp số nhân trong cộng đồng.

Khi đối chiếu với các công trình nghiên cứu trong nước về dịch tả và sốt xuất huyết tại Hà Nội hay các nghiên cứu quốc tế ứng dụng thuật toán học máy phân lớp, kết quả nghiên cứu này tương đồng về việc xác định các biến khí tượng là yếu tố cảnh báo cốt lõi. Tuy nhiên, đóng góp khác biệt của luận văn là đã chuyển hóa các mô hình toán học tĩnh thành công cụ hiển thị động 3 chiều trên nền tảng Python.

Biểu diễn trực quan: Dữ liệu được thể hiện qua các biểu đồ cột chồng 3D trong không gian đa chiều, liên kết đồng thời 4 trục thông tin: địa bàn hành chính, mốc thời gian 12 tháng, lượng mưa, độ ẩm và số ca bệnh. Biểu đồ Histogram kiểm tra phần dư hồi quy chuẩn hóa thể hiện đường cong phân phối chuẩn hình chuông với giá trị trung bình Mean xấp xỉ 0, cùng biểu đồ phân tán Scatter Plot có các điểm dư tập trung bám sát đường tung độ 0, chứng minh tính vững chắc và chuẩn xác của mô hình phân tích.

Đề xuất và khuyến nghị

Thứ nhất, nâng cấp và tích hợp hệ thống phần mềm trực quan hóa 3D vào hệ thống y tế công cộng: Trung tâm Kiểm soát Bệnh tật (CDC) tỉnh Bình Dương chủ trì phối hợp với Sở Thông tin và Truyền thông đóng gói công cụ Python thành module trực tuyến. Mục tiêu trong 6 tháng đầu năm 2021 giảm 50% thời gian trích xuất báo cáo dịch tễ định kỳ và tự động đưa ra cảnh báo sớm trước 3 tuần khi các chỉ số lượng mưa và độ ẩm vượt ngưỡng an toàn.

Thứ hai, triển khai chiến dịch diệt lăng quăng và phun hóa chất chủ động theo phân vùng nguy cơ: Ủy ban nhân dân thành phố Thủ Dầu Một, thị xã Thuận An và thị xã Dĩ An chỉ đạo mạng lưới y tế cơ sở thực hiện tổng vệ sinh môi trường, diệt lăng quăng trước mùa mưa (bắt đầu từ tháng 4 hằng năm). Mục tiêu giảm 35% chỉ số mật độ bọ gậy và hạ tỷ lệ mắc sốt xuất huyết xuống dưới mức 150 ca trên 100.000 dân trong giai đoạn 2021 - 2023.

Thứ ba, thiết lập cơ chế chia sẻ dữ liệu liên ngành khí tượng - y tế theo thời gian thực: Sở Tài nguyên và Môi trường tỉnh Bình Dương phối hợp cùng Trung tâm Y tế Dự phòng chuẩn hóa quy trình kết nối API cung cấp dữ liệu lượng mưa, nhiệt độ, độ ẩm từ các trạm quan trắc tự động theo chu kỳ 7 ngày một lần. Mục tiêu hoàn thiện mạng lưới đồng bộ 100% dữ liệu khí tượng phục vụ dự báo từ quý 3 năm 2021.

Thứ tư, đẩy mạnh truyền thông phòng bệnh có trọng điểm tại các khu nhà trọ công nhân: Ban Quản lý các Khu công nghiệp tỉnh Bình Dương phối hợp với tổ chức công đoàn địa phương triển khai tuyên truyền trực quan về vệ sinh dụng cụ chứa nước sinh hoạt. Mục tiêu đạt trên 90% các khu nhà trọ tại 41 phường nội thị thực hiện cam kết diệt bọ gậy hàng tuần trong suốt các tháng 5 đến tháng 10.

Đối tượng nên tham khảo luận văn

Nhóm 1: Cán bộ quản lý và chuyên viên dịch tễ học tại Trung tâm Y tế Dự phòng, Trung tâm Kiểm soát Bệnh tật (CDC) các tuyến. Luận văn cung cấp phương pháp quan sát trực quan khối không gian - thời gian 3D, giúp theo dõi sát sao tiến trình dịch bệnh, khoanh vùng ổ dịch chính xác và lập kế hoạch dự trù cơ số thuốc, hóa chất dập dịch kịp thời.

Nhóm 2: Học viên cao học, sinh viên và nghiên cứu sinh chuyên ngành Hệ thống thông tin, Khoa học dữ liệu và Tin học y tế. Đề tài là tài liệu tham khảo thực nghiệm giá trị về quy trình tiền xử lý dữ liệu y tế thô, kỹ thuật lập trình biểu đồ 3 chiều bằng ngôn ngữ Python (Matplotlib, Tkinter) và phương pháp xử lý hồi quy trên SPSS 20.

Nhóm 3: Các nhà hoạch định chính sách y tế và lãnh đạo Sở Y tế địa phương. Luận văn hỗ trợ xây dựng cơ sở khoa học cho các đề án chuyển đổi số y tế dự phòng, tích hợp dữ liệu môi trường - dịch tễ vào đề án xây dựng thành phố thông minh.

Nhóm 4: Giảng viên các trường đại học khối ngành công nghệ thông tin và y tế công cộng. Công trình là case-study thực tế tiêu biểu phục vụ giảng dạy các học phần phân tích dữ liệu đa biến, hệ thống thông tin địa lý y tế (Health-GIS) và mô hình hóa dịch tễ học ứng dụng.

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

Câu hỏi 1: Tại sao nghiên cứu lại kết hợp các biến khí hậu như nhiệt độ, độ ẩm và lượng mưa khi phân tích dịch sốt xuất huyết? Trả lời: Sự phát triển và truyền bệnh của muỗi Aedes aegypti phụ thuộc trực tiếp vào điều kiện thời tiết. Dữ liệu thực tế tại Bình Dương giai đoạn 2014 - 2018 chứng minh lượng mưa lớn và độ ẩm đạt 86% vào tháng 9 đã thúc đẩy muỗi sinh sản nhanh, dẫn đến số ca mắc sốt xuất huyết tăng vọt tại các vùng đô thị.

Câu hỏi 2: Ưu điểm nổi bật của việc trực quan hóa dữ liệu 3D bằng Python so với bảng biểu truyền thống là gì? Trả lời: Mô hình trực quan 3D cho phép tích hợp đồng thời 4 biến số trên cùng một giao diện hiển thị. Chuyên gia dịch tễ có thể áp dụng phương pháp nhìn - hiểu để bao quát toàn diện mối tương quan phức tạp giữa 9 huyện thị, 60 tháng và các thông số khí tượng thay vì phải tra cứu thủ công qua nhiều bảng số liệu rời rạc.

Câu hỏi 3: Địa bàn nào tại tỉnh Bình Dương ghi nhận nguy cơ bùng phát dịch sốt xuất huyết cao nhất? Trả lời: Thành phố Thủ Dầu Một, thị xã Thuận An và thị xã Dĩ An là các điểm nóng có nguy cơ bùng phát dịch cao nhất. Trong các đợt dịch lớn năm 2015 và 2017, ba địa bàn công nghiệp đô thị này luôn chiếm trên 65% tổng số ca nhiễm của toàn tỉnh do mật độ dân số cao và biến động dân cư lớn.

Câu hỏi 4: Độ tin cậy của mô hình hồi quy tuyến tính trong nghiên cứu được kiểm chứng như thế nào? Trả lời: Mô hình đạt độ tin cậy khoa học cao với mức ý nghĩa tương quan sig bằng 0.000 và giá trị kiểm định Durbin-Watson đạt 2.294, khẳng định không có hiện tượng tự tương quan chuỗi. Đồ thị Histogram và Scatter Plot cho thấy phần dư phân phối chuẩn hình chuông và bám sát đường cân bằng.

Câu hỏi 5: Hướng mở rộng và phát triển tiếp theo của đề tài nghiên cứu này là gì? Trả lời: Tác giả định hướng nâng cấp hệ thống để hiển thị trực tiếp dữ liệu dịch bệnh trên nền tảng bản đồ số GIS tỉnh Bình Dương theo thời gian thực, đồng thời kết hợp các thuật toán học máy như Random Forest hay Decision Tree nhằm tự động hóa quy trình dự báo nguy cơ dịch bệnh.

Kết luận

  • Luận văn đã chuẩn hóa và tích hợp thành công bộ dữ liệu chuỗi thời gian 5 năm (2014 - 2018) về ca bệnh sốt xuất huyết cùng các yếu tố khí tượng trên địa bàn 9 huyện, thị xã, thành phố tỉnh Bình Dương.
  • Xây dựng thành công hệ thống trực quan hóa 3D đa biến trên ngôn ngữ lập trình Python và thư viện Matplotlib, mang lại công cụ quan sát trực quan theo phương pháp nhìn - hiểu hiệu quả cho chuyên gia y tế.
  • Chứng minh chặt chẽ mối liên hệ toán học giữa số ca mắc sốt xuất huyết với nhiệt độ, độ ẩm và lượng mưa thông qua mô hình hồi quy tuyến tính trên SPSS 20 với hệ số Durbin-Watson đạt 2.294.
  • Xác định rõ các trọng điểm dịch tễ tập trung tại các vùng đô thị công nghiệp phía Nam (Thuận An, Dĩ An, Thủ Dầu Một) vào mùa mưa từ tháng 7 đến tháng 10.
  • Đề xuất khung khuyến nghị hành động liên ngành cụ thể nhằm nâng cao hiệu lực cảnh báo sớm và tối ưu hóa phân bổ nguồn lực y tế dự phòng.

Đóng góp chính của luận văn là cầu nối ứng dụng công nghệ thông tin vào bài toán y tế công cộng thực tế tại vùng kinh tế trọng điểm phía Nam. Trong giai đoạn 2021 - 2022, mô hình cần tiếp tục được hoàn thiện trên nền tảng WebGIS để nâng cao khả năng tương tác. Các cơ quan quản lý y tế và địa phương hãy chủ động ứng dụng giải pháp trực quan hóa dữ liệu đa chiều nhằm nâng cao hiệu quả kiểm soát dịch bệnh, bảo vệ sức khỏe cộng đồng và thúc đẩy phát triển kinh tế xã hội bền vững.