Cấu trúc mật độ của vỏ Trái Đất được xác định bằng phương pháp nào?

Cấu trúc mật độ được xác định chủ yếu qua phương pháp nghịch đảo trường trọng lực Bouguer kết hợp với các mô hình vận tốc sóng đàn hồi thu được từ địa chấn sâu.

Định luật Birch có ý nghĩa gì trong nghiên cứu cấu trúc mật độ?

Định luật Birch thiết lập mối quan hệ tuyến tính giữa vận tốc sóng dọc P và khối lượng riêng của đá dưới áp suất cao, giúp chuyển đổi mô hình địa chấn thành mô hình mật độ ban đầu.

Tại sao không thể chỉ dùng số liệu trọng lực để xác định duy nhất cấu trúc mật độ?

Vì bài toán ngược trường thế trọng lực có tính không duy nhất cố hữu, nghĩa là vô số cấu hình phân bố mật độ ngầm khác nhau có thể tạo ra cùng một trường dị thường trọng lực trên mặt đất.

Mô hình CRUST 5.1 cung cấp thông tin gì về cấu trúc mật độ toàn cầu?

CRUST 5.1 cung cấp dữ liệu số hóa phân tầng về chiều dày, mật độ và vận tốc sóng đàn hồi từ trầm tích, vỏ kết tinh đến manti trên trên mạng lưới ô vuông 5 độ x 5 độ khắp toàn cầu.

Cấu trúc mật độ là gì? Mô hình địa vật lý thạch quyển

Cấu trúc mật độ là sự phân bố không gian ba chiều của khối lượng riêng vật chất bên trong các tầng cấu trúc của Trái Đất hoặc các thiên thể, phản ánh thành phần khoáng vật, trạng thái nhiệt động và các quá trình tiến hóa địa động lực học trong lòng đất. Trong địa vật lý và địa chất học cấu tạo, cấu trúc mật độ của thạch quyển và manti là nguồn gốc trực tiếp tạo nên trường trọng lực dị thường, đồng thời chi phối cơ chế cân bằng đẳng tĩnh, sự đối lưu manti và sự hình thành các bồn trũng trầm tích. Bài viết này trình bày toàn diện về cơ sở lý thuyết trường trọng lực, các quan hệ thực nghiệm chuyển đổi vận tốc sóng - mật độ, phương pháp mô hình hóa nghịch đảo và các nghiên cứu cấu trúc sâu tại khu vực Biển Đông.

Bản chất vật lý và ý nghĩa của cấu trúc mật độ thạch quyển

Vỏ Trái Đất và phần trên cùng của manti tạo thành lớp thạch quyển cứng rắn trôi nổi trên tầng quyển mềm có tính quánh dẻo cao hơn. Khối lượng riêng của đất đá thay đổi liên tục theo chiều thẳng đứng và biến thiên phức tạp theo phương ngang do sự khác biệt về thành phần thạch học, mức độ nén ép do áp suất thạch tĩnh và nhiệt độ ngầm.

Ở tầng nông, lớp trầm tích chưa gắn kết có khối lượng riêng thấp. Khi xuống sâu hơn, quá trình nén chặt và biến chất làm tăng độ đặc khít. Tầng vỏ kết tinh gồm lớp vỏ trên mang tính axit giàu fenspat và thạch anh, lớp vỏ dưới mang tính bazơ giàu pyroxen và plagioclas có khối lượng riêng lớn hơn. Dưới ranh giới Moho, đá siêu bazơ của manti trên như peridotit có khối lượng riêng vượt trội. Sự phân tầng này tạo nên cấu trúc mật độ thẳng đứng cơ bản của thạch quyển.

Cơ sở lý thuyết trường thế và dị thường trọng lực

Để nghiên cứu cấu trúc mật độ bên dưới bề mặt Trái Đất mà không thể khoan tới, các nhà địa vật lý sử dụng phép đo dị thường trọng lực. Theo định luật vạn vật hấp dẫn của Newton, thế hấp dẫn $U$ và véc-tơ gia tốc trọng trường $\mathbf{g}$ tại điểm quan sát ngoài nguồn phân bố mật độ $\rho(\mathbf{r'})$ được biểu diễn qua tích phân thể tích:

\mathbf{g}(\mathbf{r}) = \nabla U(\mathbf{r}) = -G \iiint_{V} \frac{\rho(\mathbf{r'}) (\mathbf{r} - \mathbf{r'})}{|\mathbf{r} - \mathbf{r'}|^3} dV'

Trong phương trình trên, $G$ là hằng số hấp dẫn phổ quát. Khi trừ đi giá trị trọng trường lý thuyết của mô hình Elipxoid quy chiếu chuẩn, ta thu được dị thường trọng lực. Dị thường Bouguer sau khi đã hiệu chỉnh độ cao địa hình và khối nêm địa hình là đại lượng phản ánh trung thực nhất sự phân bố bất đồng nhất của cấu trúc mật độ dưới sâu.

Thuật toán phân tích miền tần số Parker

Một trong những bước đột phá trong phương pháp tính toán trường trọng lực từ các mặt ranh giới mật độ lượn sóng là công trình kinh điển của Parker (1973). Thay vì tính toán tích phân không gian rất chậm và tốn kém tài nguyên máy tính, Parker (1973) đã áp dụng chuỗi biến đổi Fourier để thiết lập mối quan hệ giải tích trực tiếp giữa dị thường trọng trường và độ mấp mô của mặt ranh giới mật độ trong miền tần số không gian:

\mathcal{F}[\Delta g] = 2\pi G \Delta\rho \exp(-|\mathbf{k}| z_0) \sum_{n=1}^{\infty} \frac{|\mathbf{k}|^{n-1}}{n!} \mathcal{F}[h^n(\mathbf{r})]

Trong công thức trên, $\mathcal{F}$ là toán tử biến đổi Fourier hai chiều, $\Delta\rho$ là độ tương phản mật độ qua ranh giới, $\mathbf{k}$ là véc-tơ số sóng không gian, $z_0$ là chiều sâu trung bình của ranh giới và $h(\mathbf{r})$ là độ lệch địa hình bề mặt ranh giới so với mặt chuẩn. Thuật toán này cho phép tái tạo nhanh chóng hình thái ranh giới móng kết tinh và ranh giới Moho từ dữ liệu trọng lực với độ chính xác cao.

Mối quan hệ thực nghiệm giữa vận tốc sóng địa chấn và mật độ

Nghịch đảo trọng lực đơn lẻ gặp phải tính không duy nhất cố hữu của bài toán trường thế: vô số cấu hình phân bố mật độ khác nhau có thể cùng tạo ra một trường dị thường trọng lực duy nhất trên mặt đất. Để khắc phục sự đa nghiệm này, địa vật lý hiện đại liên kết chặt chẽ cấu trúc mật độ với cấu trúc vận tốc truyền sóng đàn hồi đo được từ địa chấn khúc xạ và địa chấn phản xạ sâu.

Quy luật / Mô hình	Dạng phương trình quan hệ	Tác giả và công trình	Đặc trưng ứng dụng thạch học
Định luật Birch	$V_p = a + b \cdot \rho$	Birch (1961), J. Geophys. Res.	Quan hệ tuyến tính áp dụng cho đá có khối lượng nguyên tử trung bình không đổi dưới áp suất nén đến 10 kilobar
Quan hệ Brocher	Hồi quy đa thức bậc cao giữa $V_p$ , $V_s$ và $\rho$	Brocher (2005), Bull. Seismol. Soc. Am.	Tổng hợp chuẩn hóa toàn diện cho các loại đá trầm tích, magma và biến chất trong vỏ Trái Đất
Mô hình CRUST 5.1	Lưới toàn cầu phân tầng mật độ và vận tốc	Mooney và cs. (1998), J. Geophys. Res.	Cung cấp mô hình tham chiếu toàn cầu chi tiết theo các đới kiến tạo và manti trên

Định luật Birch và cơ sở vật lý khoáng vật

Nghiên cứu thực nghiệm đột phá của Birch (1961) đã đo đạc vận tốc sóng dọc $V_p$ trong hàng trăm mẫu đá dưới áp suất nén cao trong phòng thí nghiệm. Birch phát hiện ra rằng khi áp suất thủy tĩnh đạt đến mức đủ để đóng kín các vi nứt nẻ rỗng (khoảng vài kilobar), vận tốc sóng dọc tỷ lệ thuận tuyến tính với khối lượng riêng của mẫu đá, phụ thuộc chặt chẽ vào khối lượng nguyên tử trung bình của các nguyên tố cấu thành khoáng vật. Định luật Birch trở thành nền tảng vật lý cốt lõi để chuyển đổi mô hình sóng địa chấn thành mô hình mật độ ban đầu cho bài toán nghịch đảo trọng lực.

Tổng hợp quan hệ thực nghiệm Brocher

Sau nhiều thập kỷ tích lũy dữ liệu địa chấn và giếng khoan, Brocher (2005) đã công bố công trình tổng kết toàn diện, thiết lập các phương trình hồi quy thực nghiệm đa thức bậc cao giúp chuyển đổi nhất quán giữa các thông số đàn hồi: mật độ khối, vận tốc sóng dọc, vận tốc sóng ngang và hệ số Poisson. Công thức Brocher (2005) cung cấp các hệ thức hồi quy thực nghiệm đa thức giúp chuẩn hóa việc chuyển đổi giữa các thông số đàn hồi cho vỏ lục địa, dù người sử dụng cần lưu ý các tương quan thực nghiệm này không áp dụng trực tiếp cho đá mafic giàu calci hoặc đá biến chất serpentinit.

Mô hình vỏ Trái Đất toàn cầu CRUST 5.1

Nhằm cung cấp một bức tranh toàn cảnh về sự phân bố mật độ và vận tốc sóng trong thạch quyển Trái Đất phục vụ mô phỏng địa chấn toàn cầu và hiệu chỉnh trường trọng lực vệ tinh, Mooney, Laske và Masters (1998) đã xây dựng mô hình số hóa CRUST 5.1. Mô hình này phân chia bề mặt Trái Đất theo mạng lưới ô vuông không gian với bước lưới 5 độ x 5 độ.

Trong mô hình CRUST 5.1 của Mooney và cộng sự (1998), cấu trúc thạch quyển tại mỗi nút lưới được phân chia chi tiết thành các tầng lớp vật lý riêng biệt gồm: tầng nước biển, lớp băng quyển bề mặt, tầng trầm tích mềm tầng nông, tầng trầm tích cứng kết vón, lớp vỏ kết tinh trên, lớp vỏ kết tinh giữa, lớp vỏ kết tinh dưới và tầng đỉnh manti ngay dưới ranh giới Moho. Mỗi lớp được gán giá trị chiều dày, mật độ và vận tốc sóng đàn hồi đại diện, cho phép phân biệt rõ rệt sự tương phản cấu trúc mật độ giữa vỏ lục địa cổ dày và vỏ đại dương trẻ mỏng.

Nghịch đảo liên hợp trọng lực - địa chấn tại khu vực Biển Đông

Biển Đông là một trong những biển rìa lục địa lớn và phức tạp nhất thế giới, được hình thành qua các giai đoạn tách giãn lục địa, mở đáy biển và hút chìm mảng kiến tạo. Việc làm sáng tỏ cấu trúc mật độ sâu của vỏ Trái Đất tại vùng biển này có ý nghĩa khoa học và thực tiễn to lớn đối với việc đánh giá tiềm năng tài nguyên dầu khí và khoáng sản biển sâu.

Trong công trình nghiên cứu hiện đại, Li và cộng sự (2024) công bố trên tạp chí Lithosphere đã ứng dụng phương pháp nghịch đảo liên hợp giữa dữ liệu dị thường trọng lực vệ tinh độ phân giải cao và các tuyến địa chấn sâu để xây dựng các tuyến mặt cắt cấu trúc mật độ thạch quyển hai chiều cắt qua bồn trũng Qiongdongnan thuộc khu vực Tây Bắc Biển Đông. Kết quả nghiên cứu của Li và cs. (2024) đã bóc tách rõ rệt độ tương phản mật độ giữa lớp trầm tích Kainozoi dày đặc nén chặt và móng kết tinh, đồng thời xác định chính xác sự mỏng hóa cực mạnh của vỏ lục địa khi chuyển tiếp từ thềm lục địa nông ra vùng nước sâu, phản ánh cơ chế tách giãn nhiệt và magma nâng trồi trong thời kỳ Kainozoi.

Ứng dụng thực tiễn của nghiên cứu cấu trúc mật độ

Hiểu biết chi tiết về cấu trúc mật độ lòng đất mang lại những giá trị ứng dụng thiết thực:

Thăm dò dầu khí và bể trầm tích: Bản đồ cấu trúc mật độ giúp xác định chiều sâu móng kết tinh, hình thái các khối magma xâm nhập và sự phân bố thể tích các bẫy cấu trúc chứa hydrocacbon trong các bể trầm tích thềm lục địa.
Đánh giá tai biến động đất và địa chấn kiến tạo: Sự tương phản mật độ dọc theo các đứt gãy sinh chấn sâu (như đới đứt gãy Sông Hồng) phản ánh sự tích tụ ứng suất kiến tạo và ranh giới các khối trượt mảng, hỗ trợ phân vùng rủi ro động đất chính xác hơn.
Thăm dò khoáng sản ẩn sâu: Các thân quặng kim loại sulfide đa kim hoặc quặng sắt từ có tỷ trọng rất cao tạo ra dị thường trọng lực cục bộ rõ nét, cho phép phát hiện các mỏ quặng ẩn sâu hàng trăm mét dưới lớp phủ phong hóa.
Khoa học trắc địa và định vị vệ tinh: Mô hình mật độ vỏ Trái Đất là đầu vào bắt buộc để xây dựng mô hình Geoid trọng lực có độ chính xác cao, phục vụ chuyển đổi độ cao thủy chuẩn và hiệu chuẩn quỹ đạo vệ tinh định vị toàn cầu GNSS.