Tổng quan nghiên cứu

Quá trình đô thị hóa toàn cầu đang diễn ra với tốc độ chưa từng có, khi tỷ lệ dân số sống tại các khu vực đô thị tăng từ 30% năm 1950 lên hơn 55% năm 2018 và dự kiến sẽ đạt 68% vào năm 2050 theo báo cáo Triển vọng Đô thị hóa Thế giới. Tại Việt Nam, tỷ lệ đô thị hóa cả nước đã tăng từ 19,6% năm 2009 lên 38,4% vào năm 2018 với 819 đô thị, trong đó Thành phố Hồ Chí Minh là trung tâm kinh tế trọng điểm đặc biệt. Thành phố chiếm 0,6% diện tích tự nhiên và 8,34% dân số cả nước, nhưng đóng góp tới 20,5% tổng sản phẩm quốc nội (GDP) và 27,9% giá trị sản xuất công nghiệp. Sự bùng nổ kinh tế đã thu hút làn sóng di cư mạnh mẽ, đưa mật độ dân số thành phố lên mức kỷ lục hơn 4.029 người/km² trên tổng diện tích 2095,39 km².

Tuy nhiên, tốc độ mở rộng không gian đô thị nhanh chóng và tự phát đã tạo ra áp lực nặng nề lên hạ tầng kỹ thuật, môi trường sinh thái và quỹ đất nông nghiệp. Theo thống kê của Bộ Lao động, Thương binh và Xã hội, bình quân mỗi năm Việt Nam mất khoảng 73.000 ha đất canh tác do quá trình đô thị hóa. Để giải quyết bài toán quản lý không gian và hỗ trợ quy hoạch bền vững, luận văn tập trung phân tích động thái đô thị hóa của Thành phố Hồ Chí Minh trong giai đoạn 43 năm từ năm 1979 đến năm 2022. Nghiên cứu khai thác chuỗi ảnh vệ tinh Landsat đa thời gian kết hợp số liệu dân số qua 9 mốc giám sát (1979, 1989, 1995, 2000, 2005, 2010, 2015, 2020 và 2022). Kết quả cung cấp hệ thống chỉ số định lượng bao gồm tốc độ mở rộng diện tích đô thị, tốc độ gia tăng dân số, tỷ lệ tiêu thụ đất bình quân và hệ số co giãn không gian - dân số phục vụ công tác quy hoạch phát triển đến năm 2030.

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 trên nền tảng lý thuyết đô thị hóa không gian và mô hình tương quan dân số - diện tích đô thị của Julian D. Marshall (năm 2006). Mô hình toán học này xác định mối quan hệ lũy thừa giữa tỷ lệ tăng dân số và tỷ lệ mở rộng không gian xây dựng thông qua hệ số co giãn không gian. Khi hệ số này lớn hơn 1, diện tích tiêu thụ đất của mỗi cư dân mới đến sẽ cao hơn mức bình quân hiện tại, phản ánh xu hướng phát triển dàn trải.

Bên cạnh đó, nghiên cứu áp dụng khung lý thuyết phân loại lớp phủ bề mặt viễn thám dựa trên đặc tính phản xạ phổ của 4 lớp đối tượng chính: Đất trống, Bề mặt công trình xây dựng, Thực vật và Thủy hệ. Các chỉ số đô thị hóa chuẩn hóa bao gồm: Tốc độ mở rộng đô thị (Urban Expansion - UE), Tốc độ gia tăng dân số đô thị (Population Expansion - PE) và Tỷ lệ tiêu thụ đất trên đầu người (Land Consumption Ratio - LCR). Chỉ số LCR được đối chiếu trực tiếp với Quy chuẩn Kỹ thuật Quốc gia về Quy hoạch Xây dựng QCVN 01:2021/BXD ban hành theo Thông tư 01/2021/TT-BXD của Bộ Xây dựng, quy định tiêu chuẩn đất dân dụng bình quân cho đô thị đặc biệt và loại I dao động trong khoảng 45 đến 60 m²/người.

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

Nguồn dữ liệu đầu vào bao gồm chuỗi ảnh vệ tinh Landsat cấp độ Collection 2 Tier 1 do NASA/USGS cung cấp (bao gồm các cảm biến MSS trên Landsat 1-3, TM trên Landsat 4-5, ETM+ trên Landsat 7 và OLI/TIRS trên Landsat 8-9) với độ phân giải không gian 30m. Dữ liệu dân số cấp quận, huyện được thu thập từ Niên giám Thống kê Thành phố Hồ Chí Minh và các kỳ Tổng điều tra Dân số toàn quốc từ năm 1979 đến năm 2022. Ranh giới hành chính dạng Shapefile chuẩn hệ tọa độ EPSG:4326 được chuẩn hóa qua từng thời kỳ từ trung tâm dữ liệu HCMGIS.

Cỡ mẫu huấn luyện và kiểm định bao gồm hơn 3.000 điểm ảnh phân bố đồng đều cho từng thời điểm khảo sát. Nghiên cứu áp dụng phương pháp chọn mẫu ngẫu nhiên phân tầng (Stratified Random Sampling) kết hợp khóa giải đoán ảnh tổ hợp màu và đối chiếu giá trị phổ phản xạ của các kênh sóng quá khứ để đảm bảo tính đại diện cho từng loại lớp phủ.

Lý do lựa chọn thuật toán Rừng ngẫu nhiên (Random Forest - RF) trên nền tảng điện toán đám mây Google Earth Engine (GEE) là nhờ khả năng xử lý dữ liệu phi tuyến tính xuất sắc, hạn chế hiện tượng quá khớp (overfitting) thông qua kỹ thuật kết hợp cây quyết định Bootstrapping và nguyên lý trí tuệ tập thể. Sau phân loại, kỹ thuật lọc không gian bằng Bộ lọc đa số (Majority Filter) được áp dụng để loại bỏ các điểm ảnh nhiễu độc lập, trích xuất chính xác bề mặt khu đô thị tập trung từ lớp công trình xây dựng. Độ tin cậy của mô hình được kiểm định qua Ma trận sai số (Confusion Matrix), mang lại Độ chính xác tổng thể (Overall Accuracy) đạt trên 90% và Hệ số Cohen's Kappa đạt trên 0,85 cho tất cả các mốc thời gian.

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

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

Thứ nhất, diện tích bề mặt đô thị của Thành phố Hồ Chí Minh đã ghi nhận mức tăng trưởng vượt bậc trong hơn bốn thập kỷ. Năm 1979, bề mặt khu đô thị chỉ chiếm khoảng 5% tổng diện tích tự nhiên của thành phố (tương đương khoảng 105 km²). Đến năm 2022, tỷ lệ này đã tăng vọt lên xấp xỉ 32% (khoảng 670 km²), tăng gấp hơn 6 lần so với thời điểm ban đầu.

Thứ hai, chỉ số UE cho thấy sự phân hóa không gian sâu sắc giữa khu vực nội thành và ngoại thành. Tại các quận trung tâm truyền thống như Quận 1, Quận 3 và Quận 5, chỉ số UE tiệm cận mức 0% do quỹ đất xây dựng đã được lấp đầy từ trước năm 1995. Ngược lại, tốc độ mở rộng đô thị bùng nổ mạnh mẽ tại các quận mới và khu vực ven đô như Quận 7, Quận 9, Quận 12, quận Bình Tân, thành phố Thủ Đức và huyện Bình Chánh với mức tăng trưởng diện tích xây dựng đạt từ 150% đến hơn 350% trong các giai đoạn cao điểm từ năm 1995 đến năm 2015.

Thứ ba, sự dịch chuyển của chỉ số PE và LCR phản ánh rõ nét biến động áp lực dân sinh. Dân số thành phố tăng từ khoảng 3,5 triệu người năm 1979 lên gần 9 triệu người vào giai đoạn 2019-2022. Tỷ lệ tiêu thụ đất LCR tại các quận lõi trung tâm duy trì ở mức thấp, dao động từ 15 đến 30 m²/người, thể hiện mô hình đô thị nén cao độ. Trong khi đó, tại các huyện ngoại thành như Hóc Môn, Củ Chi và Bình Chánh, chỉ số LCR tăng nhanh và vượt ngưỡng 80 đến 140 m²/người, vượt xa mức chuẩn 45 đến 60 m²/người theo quy định của QCVN 01:2021/BXD.

Thứ tư, hệ số co giãn không gian - dân số biến động phức tạp theo từng thời kỳ. Giai đoạn 2000-2015, hệ số này tại các huyện vùng ven duy trì ở mức cao từ 1,8 đến 2,6, chứng minh tốc độ tiêu thụ đất vượt xa tốc độ tăng dân số. Tuy nhiên, trong giai đoạn 2020-2022, chỉ số PE có sự sụt giảm nhẹ cục bộ tại một số khu vực ngoại thành do tác động của đại dịch Covid-19 làm đảo chiều dòng di cư lao động.

Thảo luận kết quả

Dữ liệu nghiên cứu được biểu diễn thông qua hệ thống bản đồ chuỗi thời gian phân loại lớp phủ trên Google Earth Engine và các bảng ma trận chuyển dịch không gian đa thời gian. Các biểu đồ đường hồi quy giữa dân số và diện tích xây dựng phản ánh rõ ba mô hình mở rộng đô thị điển hình: mở rộng lan tỏa mép ngoài (edge growth), mở rộng tự phát (spontaneous growth) dọc theo các trục lộ giao thông huyết mạch, và phát triển lấp đầy (infilling growth) tại các khu dân cư hiện hữu.

Nguyên nhân chính thúc đẩy sự giãn nở bề mặt đô thị là quá trình công nghiệp hóa, việc thành lập các khu chế xuất, khu công nghiệp tập trung tại các vùng ven và chính sách chuyển đổi mục đích sử dụng đất nông nghiệp. So sánh với các nghiên cứu tương tự tại Quảng Châu (Trung Quốc) của Shixiao Yu hay Hà Nội của Nông Hữu Dương, mô hình đô thị hóa của Thành phố Hồ Chí Minh mang tính chất phân tán nhanh hơn ở các vùng phụ cận. Sự gia tăng nhanh chóng của bề mặt không thấm đặt ra thách thức nghiêm trọng về tiêu thoát nước, gia tăng nguy cơ ngập úng đô thị và đòi hỏi phải có sự can thiệp quy hoạch đồng bộ.

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

Thứ nhất, kiểm soát chặt chẽ chỉ số tiêu thụ đất và hạn chế chuyển đổi đất nông nghiệp tự phát. Ủy ban nhân dân Thành phố Hồ Chí Minh cùng Sở Quy hoạch - Kiến trúc cần ban hành quy chế quản lý quy hoạch không gian, áp dụng nghiêm ngặt tiêu chuẩn QCVN 01:2021/BXD với chỉ tiêu đất dân dụng bình quân 45-60 m²/người cho các dự án phát triển đô thị mới, phấn đấu kiểm soát tỷ lệ chuyển đổi đất nông nghiệp ngoại thành dưới mức 2.500 ha/năm trong giai đoạn 2025-2030.

Thứ hai, thiết lập nền tảng WebGIS và ứng dụng viễn thám tự động phục vụ giám sát đô thị hóa. Sở Tài nguyên và Môi trường phối hợp cùng Trung tâm HCMGIS tích hợp thuật toán Random Forest trên nền tảng Google Earth Engine vào cổng thông tin địa lý thành phố. Hệ thống cần tự động cập nhật bản đồ biến động lớp phủ định kỳ 12 tháng một lần, rút ngắn thời gian phát hiện sai phạm xây dựng và lấn chiếm đất công xuống dưới 30 ngày kể từ quý 2 năm 2026.

Thứ ba, thúc đẩy mô hình phát triển đô thị định hướng giao thông công cộng (TOD). Sở Giao thông Vận tải và Ban Quản lý Đường sắt Đô thị cần tập trung quy hoạch mật độ cao, tăng hệ số sử dụng đất thêm 20% đến 35% tại các bán kính 500m quanh các nhà ga tuyến Metro số 1 (Bến Thành - Suối Tiên) và Metro số 2 (Bến Thành - Tham Lương), hoàn thành tích hợp quy hoạch chi tiết trước năm 2027 nhằm giảm thiểu áp lực giãn nở dàn trải ra các vùng ven.

Thứ tư, bảo vệ hành lang sinh thái tự nhiên và tăng cường không gian xanh thoát nước. Sở Nông nghiệp và Phát triển Nông thôn phối hợp chính quyền địa phương khoanh vùng bảo vệ nghiêm ngặt 30.000 ha rừng ngập mặn Cần Giờ và các dải đệm ven sông Sài Gòn - Đồng Nai, đảm bảo duy trì diện tích mặt nước và đất cây xanh tối thiểu 15% trong tổng diện tích các khu đô thị mở rộng đến năm 2030.

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

Nghiên cứu mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng trọng tâm:

  1. Cơ quan quản lý nhà nước và nhà hoạch định chính sách đô thị: Cán bộ thuộc Sở Quy hoạch - Kiến trúc, Sở Xây dựng, Sở Tài nguyên và Môi trường và Ủy ban nhân dân các quận, huyện có thể sử dụng trực tiếp số liệu chuỗi thời gian 43 năm và các bản đồ phân vùng UE, LCR để rà soát, điều chỉnh Quy hoạch chung Thành phố Hồ Chí Minh tầm nhìn đến năm 2040 và lập kế hoạch sử dụng đất cấp huyện.

  2. Nhà nghiên cứu, giảng viên và học viên chuyên ngành Địa tin học: Các nhà khoa học trong lĩnh vực Bản đồ, Viễn thám và GIS có thể tham khảo quy trình phân loại học máy Random Forest trên nền tảng đám mây GEE, kết hợp phương pháp lọc Majority Filter và phương pháp chuẩn hóa dữ liệu dân số lịch sử làm tài liệu giảng dạy và mở rộng nghiên cứu tại các đô thị khác.

  3. Tổ chức môi trường, biến đổi khí hậu và phát triển bền vững: Các chuyên gia môi trường có thể khai thác bản đồ bề mặt không thấm đa thời kỳ để đánh giá nguy cơ ngập lụt, hiệu ứng đảo nhiệt đô thị (Urban Heat Island) và mức độ suy thoái hệ sinh thái lưu vực sông Đồng Nai - Sài Gòn.

  4. Doanh nghiệp đầu tư hạ tầng và phát triển bất động sản: Các nhà đầu tư hạ tầng có thể phân tích xu hướng dịch chuyển dân cư và tốc độ mở rộng không gian thực tế tại các huyện vùng ven để xây dựng chiến lược phát triển dự án nhà ở xã hội và trung tâm logistics phù hợp với định hướng phát triển của thành phố.

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

  1. Tại sao nghiên cứu lại chọn dữ liệu ảnh vệ tinh Landsat thay vì các nguồn ảnh khác?
    Dữ liệu Landsat cung cấp chuỗi quan sát mặt đất liên tục dài nhất trong lịch sử viễn thám từ năm 1972 đến nay với độ đồng nhất cao về độ phân giải không gian 30m. Nguồn dữ liệu này hoàn toàn mở, cho phép phân tích biến động liên tục trong 43 năm (1979-2022) mà các vệ tinh hiện đại như Sentinel-2 hay ảnh độ phân giải siêu cao thương mại không thể đáp ứng được về chiều dài lịch sử.

  2. Thuật toán Random Forest có ưu thế gì vượt trội trong phân loại lớp phủ bề mặt?
    Random Forest là giải thuật học máy tập hợp kết hợp hàng trăm cây quyết định thông qua kỹ thuật Bagging và chọn mẫu ngẫu nhiên. Thuật toán xử lý hiệu quả hiện tượng chồng lấn phổ giữa bề mặt công trình xây dựng và đất trống, giảm thiểu sai số ngoại lai và đem lại độ chính xác tổng thể đạt trên 90%, vượt trội hơn so với các phương pháp truyền thống như Maximum Likelihood.

  3. Chỉ số LCR vượt chuẩn tại các huyện ngoại thành phản ánh vấn đề gì?
    Chỉ số LCR tại các huyện Hóc Môn, Bình Chánh đạt từ 80 đến 140 m²/người, cao hơn nhiều so với khung chuẩn 45-60 m²/người của QCVN 01:2021/BXD. Thực tế này phản ánh tình trạng đô thị hóa dàn trải, chuyển đổi đất nông nghiệp thành đất ở mật độ thấp hoặc các khu công nghiệp phân tán có tỷ lệ lấp đầy chưa cao, gây lãng phí tài nguyên đất đai.

  4. Kỹ thuật lọc không gian Majority Filter đóng vai trò gì trong việc xác định ranh giới đô thị?
    Bộ lọc đa số (Majority Filter) dựa trên ma trận pixel lân cận để loại bỏ các điểm ảnh xây dựng đơn lẻ phân tán ở vùng nông thôn và các khoảng hở nhỏ trong khu dân cư. Kỹ thuật này giúp tái cấu trúc và kết nối các bề mặt công trình xây dựng thành một lớp bề mặt đô thị tập trung, đồng nhất và phản ánh chính xác ranh giới vật lý của không gian đô thị.

  5. Việc tính toán chỉ số đô thị hóa tại Thành phố Hồ Chí Minh giai đoạn 2020-2022 có điểm gì đặc biệt?
    Giai đoạn 2020-2022 chịu tác động trực tiếp từ đại dịch Covid-19, dẫn đến sự xáo trộn tạm thời trong phân bố nhân khẩu học. Trong khi diện tích bề mặt đô thị (UE) vẫn tiếp tục mở rộng nhẹ khoảng 2%, thì tốc độ tăng dân số (PE) tại một số quận huyện bị chững lại hoặc suy giảm do làn sóng người lao động tạm thời rời thành phố về quê.

Kết luận

  • Luận văn đã xây dựng thành công cơ sở dữ liệu số hóa hoàn chỉnh về biến động không gian đô thị hóa Thành phố Hồ Chí Minh trong suốt 43 năm (1979-2022) trên nền tảng điện toán đám mây Google Earth Engine.
  • Kết quả định lượng chứng minh tỷ lệ bề mặt đô thị đã tăng từ 5% lên 32% diện tích tự nhiên toàn thành phố, phản ánh sự gia tăng diện tích xây dựng gấp hơn 6 lần nhằm đáp ứng quy mô dân số tăng từ 3,5 triệu lên gần 9 triệu người.
  • Nghiên cứu chỉ rõ quy luật phát triển lệch pha khi khu vực lõi nội thành bão hòa chỉ số UE tiến về 0%, trong khi các quận mới và huyện ngoại thành bùng nổ mở rộng không gian với chỉ số LCR vượt ngưỡng chuẩn 45-60 m²/người theo QCVN 01:2021/BXD.
  • Hệ thống chỉ số UE, PE, LCR và hệ số co giãn Marshall cung cấp luận cứ khoa học chính xác để đánh giá chất lượng tăng trưởng đô thị và nguy cơ suy giảm đất nông nghiệp.
  • Mô hình phân loại Random Forest kết hợp lọc không gian khẳng định độ tin cậy cao với độ chính xác toàn cục trên 90%, sẵn sàng chuyển giao cho các hệ thống giám sát tài nguyên trực tuyến.

Kế hoạch nghiên cứu tiếp theo sẽ tích hợp mạng nơ-ron nhân tạo (ANN) và mô hình Cellular Automata để mô phỏng, dự báo kịch bản mở rộng đô thị của Thành phố Hồ Chí Minh đến năm 2035 và 2050 trước tháng 12 năm 2026. Các nhà quản lý quy hoạch, chuyên gia nghiên cứu và đơn vị phát triển đô thị quan tâm có thể khai thác trực tiếp bộ số liệu và quy trình thuật toán của luận văn để phục vụ công tác quy hoạch không gian bền vững.