ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM XÂY DỰNG BẢN ĐỒ CHẤT LƯỢNG KHÔNG KHÍ HUYỆN TRẢNG BOM, TỈNH ĐỒNG NAI

Tóm tắt nghiên cứu

  • Câu hỏi nghiên cứu chính: Liệu dữ liệu ảnh vệ tinh quang học thế hệ mới (Landsat 9) kết hợp với hệ thống thông tin địa lý (GIS) có thể mô hình hóa chính xác sự phân bố không gian của chất lượng không khí và đánh giá vai trò của thảm thực vật tại một địa bàn công nghiệp hóa nhanh như huyện Trảng Bom, tỉnh Đồng Nai hay không?
  • Tổng quan phương pháp (Methodology Snapshot): Đề tài tích hợp dữ liệu viễn thám đa phổ từ vệ tinh Landsat 9 (độ phân giải 30m) với 50 điểm quan trắc thực địa đối chứng (đo đạc các chỉ tiêu TSP, $\text{PM}_{10}$, $\text{SO}_2$, $\text{NO}_2$, $\text{CO}$ theo QCVN 05:2013/BTNMT). Nhóm nghiên cứu tiến hành chiết xuất các chỉ số thực vật ($\text{NDVI}$, $\text{VI}$, $\text{TVI}$), tính toán chỉ số ô nhiễm không khí ($\text{API}$), loại trừ sai số môi trường và kiểm định độ tin cậy bằng ma trận sai số kết hợp hệ số Kappa.
  • Phát hiện cốt lõi (Key Findings): Bụi lơ lửng ($\text{TSP}$) là tác nhân gây ô nhiễm chủ đạo với $20/50$ điểm quan trắc vượt quy chuẩn quốc gia, tập trung cao nhất tại các trục giao thông huyết mạch ($\text{TSP}$ trung bình đạt $339,2\ \mu\text{g/m}^3$). Bản đồ $\text{API}$ viễn thám chỉ ra rằng $65,59%$ diện tích toàn huyện duy trì chất lượng không khí trong lành, $31,45%$ ô nhiễm nhẹ, trong khi ô nhiễm vừa và nặng chỉ chiếm dưới $3%$. Nghiên cứu xác lập mối tương quan nghịch biến chặt chẽ giữa chỉ số thực vật $\text{NDVI}$ và chỉ số $\text{API}$ ($R = -0,9174$), khẳng định độ chính xác phân loại của mô hình đạt $73%$ theo kiểm định Kappa.
  • Hàm ý thực tiễn (Implications): Kết quả chứng minh tính khả thi vượt trội của công nghệ viễn thám trong giám sát môi trường diện rộng với chi phí tối ưu, cung cấp cơ sở khoa học then chốt để quy hoạch đô thị sinh thái, thiết lập các vành đai xanh cách ly công nghiệp và định vị mạng lưới quan trắc tự động tại tỉnh Đồng Nai.

Bối cảnh và tầm quan trọng

Thực trạng tri thức và giới hạn phương pháp truyền thống

Giám sát chất lượng môi trường không khí đóng vai trò sống còn trong việc bảo vệ sức khỏe cộng đồng và duy trì cân bằng sinh thái. Phương pháp quan trắc truyền thống thông qua các trạm mặt đất cố định hoặc lấy mẫu định kỳ tuy có độ chính xác cục bộ cao nhưng tồn tại nhược điểm lớn: chi phí đầu tư đắt đỏ, mật độ điểm đo thưa thớt và không thể phản ánh bức tranh phân bố không gian liên tục trên quy mô toàn vùng. Trong khi đó, việc ứng dụng công nghệ viễn thám và GIS đã được nhiều quốc gia tiên tiến thử nghiệm thành công để giải quyết bài toán nội suy không gian và mô phỏng ô nhiễm khí quyển.

Khoảng trống nghiên cứu (Research Gap)

Tại Việt Nam, các nghiên cứu ứng dụng ảnh vệ tinh Landsat trong đánh giá chất lượng không khí chủ yếu tập trung tại vùng miền núi phía Bắc (như Hòa Bình, Sơn La) hoặc các đô thị lớn độc canh dịch vụ. Khu vực kinh tế trọng điểm phía Nam – đặc biệt là các huyện cửa ngõ công nghiệp đang chuyển dịch cơ cấu thần tốc như Trảng Bom (Đồng Nai) – vẫn thiếu vắng một công trình thực nghiệm hoàn chỉnh ứng dụng thế hệ cảm biến Landsat 9 mới phóng để lượng hóa mối tương quan giữa độ che phủ thực vật và mức độ ô nhiễm không khí.

Tính cấp bách và thời điểm nghiên cứu

Huyện Trảng Bom sở hữu vị trí chiến lược dọc tuyến Quốc lộ 1A, giáp ranh thành phố Biên Hòa và nằm trong vùng ảnh hưởng trực tiếp của dự án Cảng hàng không Quốc tế Long Thành. Với 4 Khu công nghiệp quy mô lớn (Hố Nai, Sông Mây, Bàu Xéo, Giang Điền) cùng 8 cụm công nghiệp địa phương thu hút hàng trăm dự án FDI, Trảng Bom đối mặt với áp lực kép: bùng nổ phương tiện cơ giới hạng nặng và tốc độ đô thị hóa nhanh hướng tới mục tiêu nâng cấp lên thị xã vào năm 2030. Ô nhiễm bụi và khí thải công nghiệp đang trở thành điểm nóng xã hội, đòi hỏi giải pháp giám sát diện rộng kịp thời.

Tác động tiềm năng

Công trình nghiên cứu không chỉ lấp đầy khoảng trống dữ liệu môi trường cơ sở tại Đồng Nai mà còn cung cấp công cụ trực quan hóa hiện đại, hỗ trợ các nhà quản lý kiểm soát ô nhiễm liên vùng, tối ưu hóa quy hoạch giao thông và kiến tạo không gian sống bền vững cho hơn 300.000 cư dân địa phương.

Methodology và approach

Thiết kế nghiên cứu và quy trình thu thập dữ liệu

Đề tài áp dụng quy trình nghiên cứu tích hợp đa nguồn dữ liệu, kết hợp chặt chẽ giữa công nghệ viễn thám không gian, hệ thống thông tin địa lý (GIS), kỹ thuật quan trắc hóa lý chuẩn hóa và điều tra xã hội học:

  1. Dữ liệu viễn thám: Sử dụng cảnh ảnh vệ tinh Landsat 9 số hiệu LC09_L1TP_124052_20220112_02_T1 thu nhận ngày 12/01/2022, cung cấp bởi Cục Khảo sát Địa chất Hoa Kỳ (USGS). Ảnh có độ phân giải không gian $30\text{m}$ (các kênh phổ phản xạ) và $15\text{m}$ (kênh toàn sắc Panchromatic), độ phủ mây dưới $1%$, bảo đảm độ tinh khiết quang phổ.
  2. Quan trắc thực địa mặt đất: Bố trí 50 điểm quan trắc hệ thống đại diện cho 4 phân vùng chức năng: Dân cư - chợ ($15$ điểm), Trục đường giao thông ($13$ điểm), Khu công nghiệp ($7$ điểm), và Khu vực Nông - lâm nghiệp ($15$ điểm). Quy trình thu mẫu và phân tích tuân thủ nghiêm ngặt QCVN 05:2013/BTNMT:
    • Bụi tổng $\text{TSP}$: Thu mẫu bằng giấy lọc qua máy hút lưu lượng $250\text{ lít/phút}$ trong $60\text{ phút}$.
    • Khí $\text{CO}$: Thu gom mẫu bằng túi khí chuyên dụng, lưu lượng $0,1\text{ lít/phút}$.
    • Khí $\text{SO}_2$ & $\text{NO}_2$: Hấp thụ hóa học qua dung dịch $\text{K}_2\text{HgCl}_4$ và hỗn hợp Acid Sulfanilic - N-(1-Naphthyl)ethylenediamine.
    • Bụi mịn $\text{PM}_{10}$: Đo đạc trực tiếp bằng thiết bị đo quang học hiện trường.
  3. Điều tra xã hội học: Phỏng vấn bán cấu trúc 50 người dân sinh sống lân cận các điểm quan trắc nhằm ghi nhận tác động của chất lượng không khí đối với sức khỏe và sinh hoạt thường nhật.

Kỹ thuật xử lý ảnh và thuật toán mô hình hóa

  • Tiền xử lý ảnh quang học: Chuyển đổi giá trị số điểm ảnh ($\text{DN}$) sang bức xạ phổ thực tế tại cảm biến ($L_\lambda$), sau đó chuẩn hóa thành phản xạ phổ đỉnh khí quyển ($P_\lambda$) kết hợp nắn chỉnh hình học trên nền hệ quy chiếu WGS84 - UTM Zone 48N.
  • Tính toán chỉ số thực vật: $$\text{NDVI} = \frac{\text{NIR} - \text{RED}}{\text{NIR} + \text{RED}} = \frac{\text{Band 5} - \text{Band 4}}{\text{Band 5} + \text{Band 4}}$$ $$\text{VI} = \rho_{\text{NIR}} - \rho_{\text{RED}}$$ $$\text{TVI} = \sqrt{\text{NDVI} + 0,5}$$
  • Mô hình hóa chỉ số ô nhiễm không khí ($\text{API}$): Sử dụng công cụ Raster Calculator trong ArcGIS để tích hợp các băng phổ nhiệt và cận hồng ngoại. Thuật toán tiến hành loại trừ sai số nền do bốc hơi nước mặt ($\text{API} < 0$) và nhiễu do sinh khối rừng quá đậm đặc ($\text{API} > 325$), quy chuẩn về 5 cấp độ ô nhiễm.

Đánh giá độ tin cậy và kiểm định thống kê

Độ chính xác của bản đồ $\text{API}$ viễn thám được thẩm định thông qua ma trận sai số phân loại ($Error\ Matrix$) đối sánh trực tiếp với 50 điểm thực địa. Hệ số tương đồng Kappa ($K$) được tính toán để xác lập độ tin cậy thống kê: $$K = \frac{T - E}{1 - E}$$ (Trong đó $T$ là độ chính xác toàn cục, $E$ là kỳ vọng phân loại ngẫu nhiên).

Phát hiện chính

1. Bụi lơ lửng ($\text{TSP}$) là tác nhân gây ô nhiễm vượt ngưỡng nghiêm trọng nhất

Kết quả phân tích 50 điểm đo thực nghiệm cho thấy bụi tổng $\text{TSP}$ là yếu tố ô nhiễm đơn lẻ vượt quy chuẩn kỹ thuật quốc gia tại nhiều vị trí ($20/50$ trạm đo, chiếm $40%$). Ô nhiễm tập trung chủ yếu ở các trục hành lang vận tải cơ giới:

  • Khu vực đường giao thông: Giá trị $\text{TSP}$ trung bình đạt $339,2\ \mu\text{g/m}^3$ (vượt ngưỡng cho phép $300\ \mu\text{g/m}^3$ của QCVN 05:2013/BTNMT). Điểm cực đại ghi nhận tại nút giao Quốc lộ 1A - Ngã ba Trị An và đường tránh Võ Nguyên Giáp lên tới $400\ \mu\text{g/m}^3$.
  • Khu dân cư - chợ: Nồng độ $\text{TSP}$ trung bình đạt $234,0\ \mu\text{g/m}^3$ (cực đại $330\ \mu\text{g/m}^3$ tại khu vực chợ trung tâm).
  • Khu công nghiệp: Đạt trung bình $213,6\ \mu\text{g/m}^3$ (cực đại $275\ \mu\text{g/m}^3$).
  • Khu nông - lâm nghiệp: Môi trường duy trì mức độ an toàn cao với $\text{TSP}$ trung bình chỉ $77,0\ \mu\text{g/m}^3$.
Thông số quan trắc Đơn vị Dân cư - Chợ (TB) Giao thông (TB) Khu công nghiệp (TB) Nông - Lâm nghiệp (TB) QCVN 05:2013
Bụi tổng ($\text{TSP}$) $\mu\text{g/m}^3$ $234,0$ $339,2$ (vượt) $213,6$ $77,0$ $300$
Bụi mịn ($\text{PM}_{10}$) $\mu\text{g/m}^3$ $26,5$ $30,2$ $25,7$ $20,7$ $150$
Lưu huỳnh đioxit ($\text{SO}_2$) $\mu\text{g/m}^3$ $18,1$ $21,3$ $16,3$ $14,5$ $350$
Nitơ đioxit ($\text{NO}_2$) $\mu\text{g/m}^3$ $15,3$ $34,9$ $22,1$ $12,9$ $200$
Cacbon monoxit ($\text{CO}$) $\mu\text{g/m}^3$ $450,2$ $465,4$ $447,1$ $403,3$ $30.000$

2. Nồng độ các khí độc hại ($\text{SO}_2$, $\text{NO}2$, $\text{CO}$) và bụi mịn ($\text{PM}{10}$) nằm trong ngưỡng an toàn

Trái ngược với cảnh báo về phát thải công nghiệp nặng, các chỉ tiêu hóa lý dạng khí đều nằm dưới ngưỡng quy chuẩn nhiều lần:

  • $\text{SO}_2$ dao động trong khoảng $14,5 - 21,3\ \mu\text{g/m}^3$ (thấp hơn giới hạn cho phép $16$ lần).
  • $\text{NO}_2$ đạt cao nhất tại dải giao thông ($34,9\ \mu\text{g/m}^3$), thấp hơn quy chuẩn ($200\ \mu\text{g/m}^3$) khoảng $6$ lần.
  • $\text{CO}$ dao động từ $403,3 - 465,4\ \mu\text{g/m}^3$, thấp hơn $60$ lần so với trần quy chuẩn ($30.000\ \mu\text{g/m}^3$).
  • $\text{PM}_{10}$ duy trì ổn định dưới mức $38\ \mu\text{g/m}^3$ (giới hạn an toàn $150\ \mu\text{g/m}^3$).

3. Phân bố không gian từ bản đồ viễn thám: $97%$ diện tích huyện đạt mức an toàn

Bản đồ chất lượng không khí chiết xuất từ ảnh Landsat 9 thể hiện cấu trúc không gian môi trường phân hóa rõ rệt:

  • Cấp 1 - Không khí trong lành: Chiếm diện tích áp đảo với $21.230,6\text{ ha}$ ($65,59%$ diện tích tự nhiên), bao phủ các xã thuần nông và đất lâm nghiệp phía Bắc như Bàu Hàm, Thanh Bình, Cây Gáo, Sông Thao.
  • Cấp 2 - Ô nhiễm nhẹ: Đạt $10.180,1\text{ ha}$ ($31,45%$), phân bố tại vùng đệm nông thôn - đô thị hóa phụ cận thị trấn.
  • Cấp 3 - Ô nhiễm vừa phải: Chiếm $496,1\text{ ha}$ ($1,53%$).
  • Cấp 4 - Ô nhiễm nặng: Chiếm $463,9\text{ ha}$ ($1,43%$), hình thành các "vệt ô nhiễm" tập trung cục bộ dọc tuyến Quốc lộ 1A, Tỉnh lộ DT767, tuyến tránh Võ Nguyên Giáp và các cụm công nghiệp Hố Nai, Sông Mây, Bắc Sơn.

4. Tương quan nghịch cực đại giữa lớp phủ thực vật và chỉ số ô nhiễm ($R = -0,9174$)

Mô hình phân tích hồi quy phi không gian giữa giá trị thực phủ $\text{NDVI}$ và chỉ số ô nhiễm không khí $\text{API}$ ghi nhận hệ số tương quan tuyến tính $R = -0,9174$. Kết quả chứng minh: tại các khu vực có mật độ thảm thực vật che phủ cao, hiện tượng lắng đọng bụi tự nhiên và hấp thu vi khí hậu diễn ra mạnh mẽ, giúp duy trì chỉ số $\text{API}$ ở mức lý tưởng ($0 - 50$). Ngược lại, khu vực bề mặt bê tông hóa cao, thiếu cây xanh luôn đi kèm chỉ số $\text{API}$ vượt ngưỡng $200$.

  Chỉ số API (Ô nhiễm)
    0.0    0.2    0.4    0.6    0.8    1.0
             [ Hệ số tương quan: R = -0.9174 ]

5. Kiểm định độ tin cậy kỹ thuật và phản hồi cộng đồng

  • Mô hình giải đoán $\text{API}$ từ ảnh Landsat 9 đạt hệ số Kappa $73%$, chứng minh độ tương thích cao giữa dữ liệu viễn thám và quan trắc mặt đất.
  • Khảo sát xã hội học cho thấy $75%$ người dân cảm nhận rõ rệt tác động tiêu cực của bụi và tiếng ồn sinh hoạt/giao thông đối với hệ hô hấp, trong khi $25%$ cư dân vùng sinh thái nông nghiệp hoàn toàn hài lòng với chất lượng không khí.

Đóng góp khoa học

Đóng góp về mặt lý thuyết

Nghiên cứu làm phong phú thêm cơ sở lý luận về mối quan hệ tương hỗ giữa đặc tính phản xạ phổ bề mặt (Surface Reflectance), sinh khối thực vật ($\text{NDVI}$) và quá trình khuếch tán hạt bụi lơ lửng trong tầng đối lưu khí quyển tại khu vực nhiệt đới gió mùa cận xích đạo. Công trình xác lập hệ số tương quan định lượng ($R = -0,9174$) có giá trị tham chiếu cao cho các mô hình sinh thái cảnh quan đô thị.

Đổi mới về phương pháp luận

  • Chuẩn hóa quy trình chiết xuất $\text{API}$ đa phổ ứng dụng trực tiếp cảm biến Landsat 9 OLI-2/TIRS-2, tối ưu hóa các bước lọc nhiễu nền phản xạ hơi nước và tán xạ thực vật dày đặc.
  • Thiết lập khung phương pháp tiếp cận đa tầng: kết hợp viễn thám phân giải trung bình ($30\text{m}$) với quan trắc hóa lý chuẩn mực và điều tra xã hội học, tạo ra phương pháp luận hoàn chỉnh, dễ chuyển giao với kinh phí thấp cho các cơ quan quản lý cấp huyện/thị xã.

Ứng dụng thực tiễn và gợi ý chính sách

  1. Quy hoạch đai xanh cách ly: Cung cấp luận cứ để chính quyền huyện Trảng Bom ưu tiên phát triển dải cây xanh cách ly rộng tối thiểu $50 - 100\text{m}$ dọc theo Quốc lộ 1A và xung quanh 4 KCN tập trung, sử dụng các loài cây có chỉ số diện tích lá ($\text{LAI}$) lớn để giữ bụi.
  2. Hiện đại hóa hạ tầng quan trắc: Định vị chính xác các tọa độ điểm nóng ô nhiễm (như ngã ba Trị An, xã Hố Nai 3, xã Bắc Sơn) để ưu tiên lắp đặt trạm quan trắc không khí tự động liên tục ($\text{CEMS}$).
  3. Phân luồng giao thông và kiểm soát khí thải: Ban hành quy định rửa xe tải chuyên dụng xuất phát từ các mỏ vật liệu xây dựng và cụm công nghiệp, đẩy mạnh phân luồng giảm tải cho Quốc lộ 1A vào giờ cao điểm.

Đối tượng quan tâm

  • Các nhà nghiên cứu học thuật (GIS, Viễn thám, Khoa học Môi trường): Kế thừa quy trình tính toán chỉ số $\text{API}$ từ ảnh Landsat 9; ứng dụng tập dữ liệu đối soát thực địa 50 điểm để phát triển các mô hình học máy (Machine Learning / Deep Learning) trong dự báo chất lượng không khí.
  • Cơ quan quản lý nhà nước (Sở TN&MT tỉnh Đồng Nai, UBND huyện Trảng Bom): Sử dụng trực tiếp tập bản đồ số hóa tỷ lệ 1:50.000 phục vụ công tác thanh kiểm tra môi trường, xây dựng báo cáo hiện trạng môi trường 5 năm và thẩm định quy hoạch đô thị Trảng Bom lên thị xã.
  • Doanh nghiệp hạ tầng và Ban Quản lý KCN: Tham chiếu dữ liệu phát tán ô nhiễm cục bộ để tối ưu hóa thiết kế hệ sinh thái công nghiệp sinh thái, nâng cấp hệ thống xử lý khí thải và triển khai các cam kết ESG (Môi trường - Xã hội - Quản trị).
  • Cộng đồng dân cư và tổ chức xã hội: Nâng cao nhận thức về các điểm nóng ô nhiễm bụi, từ đó chủ động các biện pháp bảo vệ sức khỏe hệ hô hấp và tham gia giám sát cộng đồng đối với các nguồn phát thải.

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

1. Phát hiện quan trọng nhất của công trình nghiên cứu này là gì?

Phát hiện cốt lõi là bụi tổng lơ lửng ($\text{TSP}$) là tác nhân gây ô nhiễm duy nhất vượt ngưỡng quy chuẩn, tập trung chủ yếu tại các trục giao thông lớn (trung bình $339,2\ \mu\text{g/m}^3$), trong khi các khí độc hại ($\text{SO}_2$, $\text{NO}_2$, $\text{CO}$) đều ở mức an toàn. Đồng thời, nghiên cứu chứng minh mối tương quan nghịch biến rất cao ($R = -0,9174$) giữa độ che phủ thực vật ($\text{NDVI}$) và chỉ số ô nhiễm không khí ($\text{API}$).

2. Phương pháp nghiên cứu có điểm gì đặc biệt so với các đề tài trước đây?

Nghiên cứu là một trong những công trình tiên phong tại vùng Đông Nam Bộ ứng dụng thành công thế hệ ảnh vệ tinh Landsat 9 (phóng cuối năm 2021) có dải nhạy quang phổ và độ phân giải vô tuyến 14-bit vượt trội, kết hợp thuật toán hiệu chỉnh lọc sai số nền mặt nước và rừng nguyên sinh, đạt độ tin cậy kiểm định Kappa lên tới $73%$.

3. Kết quả nghiên cứu có thể tổng quát hóa (Generalize) cho các địa phương khác không?

Quy trình công nghệ hoàn toàn có thể nhân rộng và áp dụng hiệu quả cho các đô thị công nghiệp có đặc điểm địa hình đồi thoải lượn sóng và khí hậu nhiệt đới cận xích đạo tương đồng tại vùng Đông Nam Bộ như Long Thành, Nhơn Trạch (Đồng Nai), Bến Cát, Tân Uyên (Bình Dương) hay Phú Mỹ (Bà Rịa - Vũng Tàu).

4. Bước đi tiếp theo (Next Steps) để phát triển hướng nghiên cứu này là gì?

Cần mở rộng xây dựng chuỗi dữ liệu đa thời gian (Time-series analysis) nhằm so sánh động thái ô nhiễm giữa mùa mưa và mùa khô; tích hợp thêm tư liệu vệ tinh độ phân giải cao (Sentinel-2, Sentinel-5P TROPOMI chuyên đo nồng độ khí vết) và ứng dụng mạng nơ-ron nhân tạo (ANN) để tự động hóa quá trình hiệu chỉnh sai số khí quyển.

5. Các ứng dụng thực tế có thể triển khai ngay lập tức từ đề tài là gì?

Chính quyền địa phương có thể tích hợp lớp bản đồ $\text{API}$ vào cổng thông tin địa lý WebGIS quản lý môi trường huyện Trảng Bom, sử dụng làm căn cứ pháp lý để yêu cầu các công trình xây dựng hạ tầng giao thông bắt buộc phải có phương án phun sương dập bụi và lập đề án trồng cây xanh cách ly tại các xã có chỉ số ô nhiễm Cấp 3 và Cấp 4.

Kết luận

Báo cáo tổng kết đề tài khoa học công nghệ năm 2022 do Ths. Võ Minh Hoàn chủ nhiệm đã khẳng định tính ưu việt, chuẩn xác và kinh tế của việc ứng dụng công nghệ viễn thám Landsat 9 kết hợp GIS trong việc xây dựng bản đồ chất lượng không khí tại huyện Trảng Bom, tỉnh Đồng Nai. Kết quả nghiên cứu cung cấp bức tranh không gian toàn diện: $97%$ diện tích huyện có chất lượng không khí trong lành và an toàn, nhưng đồng thời cảnh báo tình trạng ô nhiễm bụi lơ lửng cục bộ tại các nút giao thông huyết mạch.

Hướng nghiên cứu tương lai cần tập trung vào việc mô hình hóa chuỗi dữ liệu đa thời gian và tích hợp thuật toán trí tuệ nhân tạo để nâng cao năng lực cảnh báo sớm. Đây là thời điểm then chốt để các nhà quản lý môi trường và quy hoạch đô thị tỉnh Đồng Nai ứng dụng trực tiếp các thành quả khoa học này vào thực tiễn, hướng tới mục tiêu phát triển Trảng Bom thành một thị xã công nghiệp - đô thị sinh thái bền vững.