Biến động bề mặt đất bao gồm những dạng chuyển dịch nào?

Biến động bề mặt đất bao gồm hai dạng chuyển dịch hình học chính là chuyển dịch thẳng đứng và chuyển dịch nằm ngang. Chuyển dịch thẳng đứng gồm hiện tượng sụt lún hạ thấp độ cao hoặc nâng cao địa hình, trong khi chuyển dịch nằm ngang biểu hiện qua sự xê dịch vị trí phẳng của các khối đất đá dọc theo các đới đứt gãy hoặc sườn dốc.

Kỹ thuật radar giao thoa (InSAR) quan trắc biến động bề mặt đất như thế nào?

Kỹ thuật InSAR sử dụng hai hoặc nhiều ảnh radar chụp cùng một khu vực ở các thời điểm khác nhau để đo đạc sự sai khác pha của sóng vi ba phản xạ. Thông qua mô hình phân rã pha giao thoa, kỹ thuật này tách biệt thành phần pha biến dạng bề mặt đất chiếu theo hướng ngắm của vệ tinh khỏi các thành phần pha địa hình, sai số quỹ đạo và trễ khí quyển.

Độ chính xác đo đạc biến dạng mặt đất của phương pháp PSI đạt mức nào?

Phương pháp radar tán xạ phản xạ cố định (PSI) có khả năng đo đạc chuyển dịch và biến dạng bề mặt đất với độ chính xác đạt mức milimet qua các chuỗi thời gian quan sát nhiều năm. Kỹ thuật này đạt được độ chính xác cao nhờ tập trung phân tích các mục tiêu phản xạ ổn định trên mặt đất như công trình xây dựng, cầu cống hoặc các mỏm đá lộ đầu.

Nguyên nhân chính gây sụt lún bề mặt đất tại Đồng bằng sông Cửu Long là gì?

Hiện tượng sụt lún bề mặt đất tại Đồng bằng sông Cửu Long phát sinh từ sự kết hợp giữa quá trình nén chặt tự nhiên của các tầng trầm tích trẻ hạt mịn và việc khai thác nước ngầm liên tục phục vụ sinh hoạt, nông nghiệp. Sự suy giảm áp lực nước lỗ rỗng làm gia tăng ứng suất hữu hiệu, dẫn đến cố kết các tầng sét và hạ thấp cao độ bề mặt đất tích lũy qua nhiều năm.

Biến động bề mặt đất (ground surface deformation) là gì?

Biến động bề mặt đất là tập hợp các quá trình thay đổi về hình thái hình học, vị trí không gian và trạng thái động lực của bề mặt vỏ Trái Đất diễn ra theo không gian và thời gian dưới tác động của các hiện tượng tự nhiên hoặc hoạt động nhân sinh. Khái niệm này bao gồm các hiện tượng chuyển dịch thẳng đứng như sụt lún và nâng hạ địa hình, cũng như các chuyển dịch ngang phát sinh từ hoạt động kiến tạo mảng hoặc biến dạng nền đất. Bài viết trình bày hệ thống phân loại, cơ sở mô hình toán học và nguyên lý vật lý của các kỹ thuật quan trắc hiện đại, phân tích dữ liệu thực tiễn và những thách thức kỹ thuật then chốt trong nghiên cứu động lực bề mặt đất.

Bản chất khoa học và phân loại biến động bề mặt đất

Bề mặt Trái Đất không phải là một thực thể tĩnh tại mà liên tục chịu tác động của các trường ứng suất nội sinh và quá trình ngoại sinh. Việc nghiên cứu biến động bề mặt đất đòi hỏi phải phân định rõ ràng giữa các cơ chế phát sinh tự nhiên và những tác động nhân sinh làm gia tốc quá trình biến dạng.

Phân loại theo cơ chế nguồn gốc

Dựa trên nguồn gốc phát sinh năng lượng và lực tác động, biến động bề mặt đất được chia thành ba nhóm cơ bản:

Biến động nội sinh do kiến tạo: Bao gồm chuyển dịch của các khối mảng thạch quyển dọc theo các đới đứt gãy hoạt động, hoạt động nâng hạ đẳng tĩnh sau băng hà, chuyển dịch nền đất trong các trận động đất và biến dạng vòm magma tại các khu vực núi lửa. Các chuyển dịch này thường diễn ra trên quy mô vùng rộng lớn hoặc theo các hành lang kiến tạo xác định.
Biến động ngoại sinh tự nhiên: Phát sinh từ các quá trình địa mạo và địa chất thủy văn như quá trình nén chặt tự nhiên của các tầng trầm tích trẻ hạt mịn châu thổ, hiện tượng hòa tan hình thành hố sụt karst tại các vùng đá vôi, xói mòn bề mặt do dòng chảy mặt, và các chuyển động trượt lở mái dốc vùng đồi núi.
Biến động do tác động nhân sinh: Xuất hiện do các can thiệp cơ học và thủy văn của con người vào vỏ Trái Đất. Các nguyên nhân hàng đầu gồm có việc khai thác quá mức nước dưới đất làm suy giảm áp lực thủy tĩnh trong các tầng chứa nước, khai thác khoáng sản và năng lượng hóa thạch dưới lòng đất, tải trọng tĩnh từ các cụm công trình xây dựng đô thị tập trung, và tải trọng động do hệ thống giao thông vận tải.

Phân loại theo phương thức chuyển dịch hình học

Về mặt động học hình học, chuyển dịch bề mặt đất được phân tích thành các thành phần vector không gian:

Chuyển dịch thẳng đứng: Gồm hiện tượng sụt lún địa hình khi bề mặt hạ thấp độ cao tương đối so với mặt chuẩn trắc địa, và hiện tượng nâng cao khi bề mặt trồi lên do dỡ tải hoặc hoạt động ép trồi địa chất.
Chuyển dịch nằm ngang: Sự dịch chuyển vị trí phẳng của các điểm bề mặt đất, phổ biến dọc theo các vết nứt đứt gãy trượt bằng hoặc hiện tượng dịch chuyển trôi ngang của các khối đất đá trên các bề mặt trượt sườn dốc.

Dạng biến động	Cơ chế vật lý chủ đạo	Tác nhân kích hoạt	Đặc trưng không gian
Sụt lún trầm tích châu thổ	Suy giảm áp lực lỗ rỗng và nén chặt cố kết tầng sét	Khai thác nước dưới đất và nén chặt tự nhiên	Quy mô đồng bằng rộng lớn, tiến triển liên tục
Biến dạng đứt gãy kiến tạo	Tích lũy và giải phóng ứng suất cắt đàn hồi	Chuyển động mảng thạch quyển	Dải tuyến tính dọc đới đứt gãy, gián đoạn chu kỳ
Sụt trũng cục bộ đô thị	Gia tăng ứng suất hữu hiệu kết hợp tải trọng công trình	Công trình xây dựng và khai thác hạ tầng ngầm	Quy mô điểm hoặc cụm đô thị tập trung
Trượt lở mái dốc	Mất cân bằng lực kháng cắt trên mặt trượt sườn đồi	Mưa lớn, ngấm nước và biến động độ dốc	Vùng sườn dốc địa hình đồi núi cục bộ

Cơ sở lý thuyết và mô hình toán học của quan trắc biến động bề mặt

Quan trắc biến động bề mặt đất hiện đại dựa trên nền tảng của kỹ thuật viễn thám radar khẩu độ tổng hợp giao thoa và mạng lưới định vị vệ tinh chính xác cao. Kỹ thuật radar giao thoa vi sai thu nhận hai hoặc nhiều ảnh radar của cùng một khu vực tại các thời điểm khác nhau để đo đạc sự sai khác pha của sóng phản xạ.

Mô hình phân rã pha giao thoa tổng quát

Pha giao thoa tổng hợp thu được giữa hai lần chụp radar vệ tinh là sự kết hợp của nhiều thành phần vật lý khác nhau theo mô hình phân rã tuyến tính:

\Delta \phi = \phi_{\text{topo}} + \phi_{\text{def}} + \phi_{\text{atm}} + \phi_{\text{orb}} + \phi_{\text{noise}}

Trong mô hình toán học này, từng đại lượng được định nghĩa cụ thể:

$\Delta \phi$ là hiệu pha giao thoa tổng hợp đo được từ cặp ảnh radar.
$\phi_{\text{def}}$ là thành phần pha do sự chuyển dịch hoặc biến dạng của bề mặt đất diễn ra giữa hai thời điểm quan sát.
$\phi_{\text{topo}}$ là thành phần pha do độ cao địa hình gây ra, có thể được loại bỏ bằng mô hình số độ cao bề mặt.
$\phi_{\text{atm}}$ là sai số trễ pha do sự thay đổi của hơi nước và chiết suất khí quyển giữa hai lần vệ tinh bay qua.
$\phi_{\text{orb}}$ là sai số pha phát sinh từ độ không chuẩn xác trong việc xác định quỹ đạo chuyển động của vệ tinh.
$\phi_{\text{noise}}$ là thành phần nhiễu pha ngẫu nhiên xuất phát từ hiện tượng mất kết hợp tán xạ hoặc nhiễu điện tử của máy thu.

Mối quan hệ hình học giữa pha và chuyển dịch theo hướng ngắm

Thành phần pha biến dạng bề mặt đất liên hệ trực tiếp với khoảng cách chuyển dịch chiếu theo hướng ngắm của vệ tinh thông qua biểu thức giải tích:

\phi_{\text{def}} = -\frac{4\pi}{\lambda} d_{\text{LOS}}

Trong đó, đại lượng $\lambda$ biểu thị bước sóng điện từ của cảm biến radar trên vệ tinh, và đại lượng $d_{\text{LOS}}$ là độ chuyển dịch của điểm phản xạ trên bề mặt đất chiếu theo hướng ngắm của vệ tinh. Nhân tử thể hiện hành trình truyền sóng vô tuyến khứ hồi hai chiều từ vệ tinh tới mục tiêu bề mặt đất và phản xạ ngược lại.

Theo phân tích tổng quan của Bürgmann và cs. (2000), phương pháp radar giao thoa vi sai chỉ phản ánh thành phần chuyển dịch một chiều chiếu theo hướng ngắm của cảm biến radar. Do đó, để xác định đầy đủ vector chuyển dịch ba chiều trong không gian, cần kết hợp dữ liệu chụp từ các quỹ đạo vệ tinh bay lên và quỹ đạo bay xuống cùng với các trạm định vị mặt đất.

Các phương pháp kỹ thuật quan trắc hiện đại

Việc theo dõi biến động bề mặt đất kết hợp đa dạng các nguồn dữ liệu viễn thám không gian và trắc địa truyền thống nhằm khắc phục nhược điểm của từng phương pháp riêng lẻ.

Kỹ thuật radar tán xạ phản xạ cố định

Kỹ thuật radar giao thoa chuỗi thời gian, tiêu biểu là phương pháp tán xạ phản xạ cố định, tập trung phân tích các mục tiêu trên mặt đất có độ ổn định tán xạ cao theo thời gian như các công trình bê tông, cầu cống, đường sá và các khối đá lộ đầu.

Theo đánh giá chuyên sâu của Crosetto và cs. (2016), kỹ thuật tán xạ phản xạ cố định có khả năng đo đạc chuyển dịch và biến dạng bề mặt đất với độ chính xác đạt mức milimet qua các chuỗi thời gian quan sát dài nhiều năm. Phương pháp này cho phép ước tính và tách biệt chính xác sai số trễ khí quyển khỏi tín hiệu biến dạng thực tế, mang lại bước tiến vượt bậc so với kỹ thuật giao thoa vi sai thông thường.

Phương pháp định vị vệ tinh và thủy chuẩn hình học

Bên cạnh viễn thám radar, các kỹ thuật quan trắc mặt đất đóng vai trò kiểm chuẩn độc lập không thể thay thế:

Mạng lưới định vị vệ tinh liên tục: Cung cấp tọa độ ba chiều liên tục theo thời gian thực với độ chính xác cao tại các điểm mốc cố định, giúp hiệu chỉnh hệ quy chiếu động và kiểm soát chuyển động mảng quy mô lớn.
Phương pháp đo cao hình học: Sử dụng máy thủy bình điện tử độ chính xác cao và mia invar để xác định độ lún đứng tại các tuyến mốc trắc địa quốc gia hoặc công trình chuyên biệt. Phương pháp này có độ tin cậy tuyệt đối về mặt cao độ cục bộ nhưng đòi hỏi kinh phí lớn và hạn chế về mật độ điểm không gian.

Phương pháp quan trắc	Độ chính xác đo đạc	Mật độ không gian	Tần suất chu kỳ	Ưu điểm nổi bật	Hạn chế chính
Radar giao thoa tán xạ cố định	Đạt mức milimet	Rất cao tại vùng đô thị	Chu kỳ theo quỹ đạo vệ tinh	Bao phủ diện rộng, quan sát hồi cứu	Mất kết hợp tại thảm thực vật rậm
Định vị vệ tinh liên tục	Đạt mức milimet đến centimet	Mật độ điểm thưa thớt	Liên tục theo thời gian thực	Đo trực tiếp vector ba chiều	Chi phí xây dựng trạm cao
Đo cao hình học thủy chuẩn	Đạt mức dưới milimet	Điểm đo dọc tuyến mốc	Chu kỳ đo lặp gián đoạn	Độ tin cậy trắc địa tuyệt đối	Tốn công sức, không liên tục diện rộng

Tác động môi trường và thực tiễn nghiên cứu tại Việt Nam

Biến động bề mặt đất, đặc biệt là hiện tượng sụt lún diện rộng, đang là vấn đề môi trường và địa kỹ thuật cấp bách trên toàn cầu cũng như tại Việt Nam.

Đánh giá nguy cơ sụt lún quy mô toàn cầu

Nghiên cứu mô hình hóa toàn cầu của Herrera-García và cs. (2021) chỉ ra có 19% dân số toàn cầu sinh sống tại các khu vực có nguy cơ sụt lún bề mặt đất và dự báo đến năm 2040 đe dọa 635 triệu người sinh sống tại các vùng sụt lún có nguy cơ ngập lụt. Hiện tượng sụt lún làm suy giảm vĩnh viễn cao độ địa hình, làm thay đổi dòng chảy tự nhiên và khuếch đại nghiêm trọng nguy cơ ngập lụt tại các vùng trũng ven biển và các thung lũng nội địa.

Thực tiễn nghiên cứu tại Đồng bằng sông Cửu Long

Đồng bằng sông Cửu Long là một trong những điểm nóng sụt lún bề mặt đất lớn nhất trên thế giới do sự kết hợp giữa trầm tích trẻ nén chặt tự nhiên và việc khai thác nước ngầm gia tăng:

Nghiên cứu của Minderhoud và cs. (2017) chỉ ra quá trình khai thác nước ngầm liên tục trong 25 năm tại Đồng bằng sông Cửu Long dẫn đến độ sụt lún tích lũy trung bình toàn vùng đạt 18 cm. Nghiên cứu cũng xác định tốc độ sụt lún trung bình toàn vùng do khai thác nước ngầm đạt 1,1 cm mỗi năm.
Phân tích dữ liệu giao thoa radar ALOS PALSAR giai đoạn 2006 đến 2010 của Erban và cs. (2014) xác định tốc độ sụt lún bề mặt đất trung bình tại các vị trí giếng quan trắc là 1,6 cm mỗi năm do suy giảm mực nước ngầm và nén chặt tầng trầm tích.

Tốc độ sụt lún bề mặt đất tại đây cao hơn đáng kể so với tốc độ nước biển dâng tương đối, dẫn đến tình trạng ngập lụt gia tăng khi triều cường, sạt lở bờ sông và xâm nhập mặn tiến sâu vào nội đồng.

Ranh giới áp dụng, nguồn sai số và thách thức kỹ thuật

Mặc dù công nghệ quan trắc viễn thám và trắc địa đã đạt nhiều bước tiến, việc xác định biến động bề mặt đất vẫn đối mặt với những giới hạn kỹ thuật nội tại:

Hiện tượng mất kết hợp pha thời gian và không gian: Trên các vùng đất nông nghiệp, rừng cây hoặc mặt nước, tính chất phản xạ sóng vô tuyến thay đổi liên tục theo mùa và thời tiết, khiến tín hiệu giao thoa bị nhiễu hoàn toàn và không thể trích xuất pha biến dạng.
Nhiễu trễ khí quyển: Sự biến thiên ngẫu nhiên của lượng hơi nước trong tầng đối lưu tạo ra các dải pha giả mạo có bước sóng không gian tương đương với biến dạng kiến tạo cục bộ, đòi hỏi phải sử dụng các mô hình dự báo thời tiết hoặc lọc chuỗi thời gian phức tạp để triệt tiêu.
Mơ hồ giải mở pha: Tín hiệu pha thu nhận chỉ nằm trong khoảng tuần hoàn từ âm pi đến dương pi. Tại các khu vực có gradient biến dạng quá dốc như sạt lở đất tức thời hoặc vết nứt đứt gãy lớn, thuật toán giải mở pha có thể mắc sai số bước nhảy bội số của chu kỳ sóng.
Phân tách các thành phần nguyên nhân: Rất khó để phân định rạch ròi bằng dữ liệu viễn thám thuần túy phần sụt lún nào do nén chặt trầm tích tự nhiên và phần nào do khai thác nước ngầm hoặc tải trọng công trình đô thị nếu thiếu các giếng quan trắc áp lực thủy tĩnh và dữ liệu khoan địa kỹ thuật đồng bộ.

Xu hướng phát triển và vấn đề nghiên cứu còn mở

Nghiên cứu biến động bề mặt đất đang chuyển dịch theo các hướng trọng tâm nhằm giải quyết các thách thức trên:

Tích hợp đa nguồn dữ liệu không gian: Đồng hóa dữ liệu radar đa băng sóng với mạng lưới trạm định vị liên tục và dữ liệu đo cao mặt đất nhằm xây dựng trường vector chuyển dịch ba chiều toàn diện và liên tục.
Ứng dụng học máy và điện toán đám mây: Xử lý chuỗi dữ liệu lớn vệ tinh mở trên các nền tảng điện toán đám mây để phát hiện sớm các dị thường biến dạng mái dốc hoặc công trình hạ tầng trọng yếu.
Mô hình hóa liên kết thủy văn địa kỹ thuật: Kết hợp chuỗi thời gian dịch chuyển bề mặt đất từ viễn thám với các mô hình dòng thấm ba chiều để dự báo khả năng phục hồi đàn hồi của tầng chứa nước và hoạch định chiến lược khai thác tài nguyên nước bền vững.