Tổng quan về luận án

Nghiên cứu cấu trúc kiến tạo vỏ Trái Đất tại các vùng chuyển tiếp lục địa - đại dương đóng vai trò nền tảng trong khoa học Địa vật lý và Địa chất biển. Luận án tiến sĩ "Nghiên cứu cấu trúc kiến tạo khu vực thềm lục địa Bắc Miền Trung Việt Nam trên cơ sở xử lý, phân tích minh giải tổng hợp tài liệu địa vật lý" của tác giả Trần Văn Khá (2018), do TS. Hoàng Văn Vượng và PGS. Đỗ Đức Thanh hướng dẫn tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, đại diện cho một công trình tiên phong về mặt phương pháp luận và ứng dụng thực tiễn trong việc làm sáng tỏ bình đồ kiến trúc sâu vùng thềm lục địa Bắc Miền Trung và quần đảo Hoàng Sa.

Về mặt khoa học, khoảng trống nghiên cứu (research gap) lớn nhất được xác định là: "vẫn còn tồn tại một vài vấn đề chưa được giải quyết thỏa đáng trong các nghiên cứu trước đây về việc xác định các hệ thống đứt gãy theo tài liệu địa vật lý cũng như cấu trúc Moho khu vực thềm lục địa Việt Nam... điều này do không có điểm tựa độ sâu bề mặt Moho trên khu vực nghiên cứu". Đa phần các tài liệu khoan dầu khí hiện hữu chỉ giới hạn ở tầng trầm tích Kainozoi nông (vài kilômét), khiến việc giải bài toán ngược mật độ trọng lực bị suy thoái nghiệm nghiêm trọng. Đồng thời, thói quen kỹ thuật chọn hệ số điều hòa $N = 2M/3$ trong phương pháp Gradient chuẩn hóa toàn phần (NFG) truyền thống gây sai số lớn khi xác định độ sâu nguồn dị thường.

Luận án thiết lập 3 câu hỏi và giả thuyết nghiên cứu tường minh:

  1. Câu hỏi 1 (Q1): Quy luật lựa chọn hệ số hài $N$ tối ưu trong chuỗi Fourier của phương pháp NFG 2D để xác định chính xác độ sâu ranh giới mật độ là gì?
    • Giả thuyết 1 (H1): Chuỗi điều hòa bậc $N$ xác lập theo cực trị phổ năng lượng (Aydın, 2007) sẽ loại trừ hiện tượng phân kỳ giả và triệt tiêu hiệu ứng Gibbs tốt hơn công thức kinh nghiệm $2M/3$.
  2. Câu hỏi 2 (Q2): Làm thế nào để khôi phục các đường biên đứt gãy kiến tạo bị đứt đoạn bởi thuật toán Gradient ngang cực đại (HGM) chuẩn của Blakely và Simpson (1986)?
    • Giả thuyết 2 (H2): Mở rộng hệ điều kiện tìm cực đại đa thức bậc 2 trong lân cận lưới $[-d, d]$ sẽ bổ sung đầy đủ các vector mật độ biên bị bỏ sót.
  3. Câu hỏi 3 (Q3): Cấu trúc bề mặt Moho và hình thái vỏ Trái Đất khu vực Bắc Miền Trung – Hoàng Sa biến thiên định lượng ra sao khi được ràng buộc bằng các điểm tựa địa chấn sâu thực tế?
    • Giả thuyết 3 (H3): Khi loại trừ dị thường trọng lực tầng trầm tích biến đổi phi tuyến và tích hợp 2 điểm tựa Moho địa chấn, ranh giới Moho sẽ phản ánh chính xác cơ chế tách giãn vỏ lục địa chuyển tiếp sang vỏ đại dương.

Phạm vi nghiên cứu bao quát không gian từ bể Sông Hồng, bể Nam Hải Nam, bể Phú Khánh tới vùng biển sâu Hoàng Sa, tích hợp tập dữ liệu khảo sát thực địa đồ sộ: hơn 6.000 km tuyến đo địa chấn - trọng lực đồng bộ của PVEP (2009), hàng ngàn kilômét tuyến đo tàu Gagarinsky (1990, 1992), tài liệu địa chấn sâu dự án CSL07-08, cùng lưới trọng lực vệ tinh GEBCO-08 (độ phân giải 30 arc-second) và Sandwell (độ phân giải 1 arc-minute).

Literature Review và Positioning

Lịch sử nghiên cứu kiến tạo và địa vật lý Biển Đông trải qua nhiều giai đoạn với các trường phái tiếp cận khác nhau:

                            TIẾN TRÌNH LỊCH SỬ NGHIÊN CỨU BIỂN ĐÔNG
   Giai đoạn 1950 - 1980       Giai đoạn 1980 - 2000              Giai đoạn 2000 - Nay
   • Khảo sát đơn lẻ           • Ký kết hợp tác Việt - Nga        • Đề tài KC-09-02 (Bùi Công Quế)
   • Chương trình NAGA         • Kulinich (1989): Tách giãn rift  • Đề tài KC.09.20 (Nguyễn Thế Tiệp)
     (Scripps/Mỹ, 1959-1974)   • Tapponnier (1993): Đùn ép trượt  • Mô hình hóa 2D/3D & Vệ tinh
   • Thiếu chiều sâu cấu trúc  • Hayes & Taylor (1983, 1986)      • Luận án Trần Văn Khá (2018)

Nhánh nghiên cứu quốc tế chứng kiến sự đối thoại học thuật sâu sắc giữa hai mô hình địa động lực: mô hình tách giãn tạo rift nhiều pha (Kulinich, 1989; Taylor & Hayes, 1983; Hall, 1996, 2002) và mô hình trượt bằng đùn ép do va chạm Ấn - Á (Brias, Patriat & Tapponnier, 1993). Tại Việt Nam, các công trình của Bùi Công Quế (1990-2005), Hoàng Văn Vượng & Đỗ Chiến Thắng (2005), Phùng Văn Phách (2001, 2014) đã đặt nền móng xử lý trường thế và phân vùng kiến tạo. Tuy nhiên, các công trình này chủ yếu xử lý dị thường khu vực quy mô nhỏ hoặc mang tính định tính trên bình đồ rộng, thiếu mô hình toán học giải bài toán ngược mật độ có kiểm chứng thực nghiệm sâu.

So sánh với các nghiên cứu quốc tế kinh điển:

  • Nghiên cứu Biển Bắc (North Sea Rift System - Murthy & Rao, 1979; Hospers, 1985): Nơi các hàm mật độ mũ được áp dụng thành công nhờ mạng lưới giếng khoan dày đặc.
  • Nghiên cứu rìa thụ động Vịnh Mexico (Blakely & Simpson, 1986; Sandwell et al., 1995): Tiêu chuẩn hóa việc nhận diện ranh giới đứt gãy bằng đạo hàm ngang nhưng bộc lộ khuyết tật đứt đoạn biên tuyến tính trong môi trường trường thế có tỷ số tín hiệu/nhiễu (SNR) thấp.

Luận án định vị chính xác ở điểm giao cắt giữa việc cải tiến giải thuật trường thế toán lý và giải quyết bài toán địa chất thực tế của Biển Đông: chuyển đổi từ phương pháp minh giải định tính sang mô hình định lượng 3 chiều có ràng buộc địa chấn sâu.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết địa vật lý trường thế

Luận án mở rộng và hoàn thiện 3 trụ cột lý thuyết trường thế:

  1. Lý thuyết Gradient chuẩn hóa toàn phần của Berezkin (1967): Xác lập lại toán tử khai triển giải tích xuống nửa không gian dưới: $$G_H(x,z) = \frac{G(x,z)}{G_{CP}(z)} = \frac{\sqrt{V_x^2(x,z) + V_z^2(x,z)}}{\frac{1}{M}\sum_{i=0}^{M-1}\sqrt{V_x^2(x_i,z) + V_z^2(x_i,z)}}$$ Loại bỏ quy tắc kinh nghiệm $N = 2M/3$, luận án chứng minh về mặt toán học rằng số hài $N$ tối ưu phải phụ thuộc vào hàm phân rã phổ năng lượng của dị thường cục bộ, kết hợp hàm làm trơn $Q_n = \left[\frac{\sin(\pi n/N)}{\pi n/N}\right]^\mu$ với chỉ số $\mu = 2$ để triệt tiêu hiệu ứng Gibbs.
  2. Lý thuyết biên Gradient ngang cực đại của Blakely và Simpson (1986): Hệ điều kiện chuẩn đòi hỏi giá trị gradient ngang $H[\Delta g]_{i,j}$ phải thỏa mãn đồng thời hai bất đẳng thức trên 4 hướng lân cận (ngang, dọc, 2 đường chéo). Luận án chứng minh điều kiện này quá khắt khe, dẫn đến việc bỏ sót các điểm cực đại thực tế của đa thức nội suy bậc hai $ax^2 + bx + c$ trên miền đối xứng $[-d, d]$. Việc nới lỏng và bổ sung hệ điều kiện toán học đã giúp khôi phục mạng lưới ranh giới mật độ liên tục.
  3. Lý thuyết mô hình hóa trọng lực 2D/3D (Talwani & Ewing, 1960; Murthy & Rao, 1979): Luận án mở rộng phương trình tích phân đường dọc chu vi đa giác với mật độ biến đổi phi tuyến dạng hàm mũ $\rho(z) = \rho(0)e^{-\beta z}$ và dạng tuyến tính $\rho(z) = \rho(0) + az$.

Khung phân tích tích hợp

     LÝ THUYẾT TRƯỜNG THẾ               ĐỊA CHẤN SÂU & GIẾNG KHOAN             ĐỊA ĐỘNG LỰC HỌC BIỂN
 (Berezkin, Blakely, Talwani)            (Tuyến CSL07-08, 115A-1X)             (Mô hình Rift đa pha)

Khung phân tích thiết lập các điều kiện biên chặt chẽ: ranh giới móng âm học trước Kainozoi có bước nhảy mật độ $\Delta\rho_1 = 0,25 - 0,40\text{ g/cm}^3$; ranh giới vỏ - manti (mặt Moho) có bước nhảy mật độ $\Delta\rho_2 = 0,45\text{ g/cm}^3$.

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

Thiết kế nghiên cứu và quy trình thực nghiệm

Nghiên cứu vận hành trên hệ hình Thực chứng luận thực nghiệm (Empirical Positivism) kết hợp giải bài toán ngược địa vật lý (Geophysical Inverse Theory). Quy trình kết hợp đa tầng dữ liệu được chuẩn hóa nghiêm ngặt:

                          SƠ ĐỒ QUY TRÌNH XỬ LÝ VÀ PHÂN TÍCH DỮ LIỆU

Các thông số thống kê và kiểm định độ tin cậy

Tuyến đo / Đối tượng Chiều dài / Quy mô Phương pháp áp dụng Tham số kỹ thuật then chốt Mục tiêu minh giải
Tuyến AA' 120 km NFG thử nghiệm Lưới mẫu $\Delta x = 0,5\text{ km}$; $N = 20$ Kiểm chứng mô hình lý thuyết lăng trụ 2D
Tuyến PKGE09-09 280 km NFG & Địa chấn sâu Nâng trường $z = 20\text{ km}$; $N = 16$ Xác định ranh giới đáy trầm tích và Moho
Lưới HGM toàn vùng $100.000\text{ km}^2$ HGM cải tiến Nâng trường $0 - 30\text{ km}$; Bước nhảy 5 km Xác định đứt gãy sâu phương KT, TB-ĐB
Điểm chuẩn Moho 1 Giao tuyến N37/09-09 Địa chấn phản xạ sâu Thời gian truyền sóng $t_0$, Vận tốc $V_p$ Cố định độ sâu Moho tuyệt đối điểm 1
Điểm chuẩn Moho 2 Giao tuyến N29/09-09 Địa chấn phản xạ sâu Chuyển đổi Time-to-Depth Cố định độ sâu Moho tuyệt đối điểm 2

Các công cụ phần mềm chuyên dụng được vận hành gồm: Geosoft Oasis Montaj, thuật toán Fortran/MATLAB giải phương trình tích phân Talwani, hệ thống số hóa C-MAP Na Uy, cùng cơ sở dữ liệu GEODAS của NGDC/NOAA.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Quy luật phân bố mật độ phi tuyến của tầng trầm tích Kainozoi: Xây dựng thành công hàm phân bố mật độ trung bình theo chiều dày trầm tích cho toàn vùng Bắc Miền Trung – Hoàng Sa từ số liệu địa chấn và giếng khoan (115A-1X, 121CM-1X): $$\rho(z) = 2,15 + 0,18 \cdot \ln(z + 1)$$ Phát hiện này bác bỏ giả định mật độ đồng nhất ($\rho = 2,30\text{ g/cm}^3$) vốn gây sai lệch tới $15 - 20%$ trong tính toán độ sâu móng trước Kainozoi của các nghiên cứu trước năm 2010.
  2. Minh định hệ thống đứt gãy sâu bằng HGM cải tiến:
    • Ở mức nâng trường $0 - 10\text{ km}$: Nhận diện rõ nét hệ thống đứt gãy đới nông phương Tây Bắc - Đông Nam (song song với đới đứt gãy Sông Hồng) và phương kinh tuyến.
    • Ở mức nâng trường $15 - 30\text{ km}$: Làm lộ diện đới đứt gãy trượt bằng - tách giãn sâu cắt qua vỏ Trái Đất phương á kinh tuyến (gần $109^\circ - 110^\circ\text{E}$), đóng vai trò ranh giới chuyển tiếp kiến tạo giữa khối Indochina và trũng sâu Biển Đông. Thuật toán HGM cải tiến đã bổ sung hơn $30%$ chiều dài các đoạn đứt gãy bị đứt quãng so với thuật toán Blakely & Simpson kinh điển.
         ĐỒ THỊ THỂ HIỆN SỰ THAY ĐỔI ĐỘ SÂU BỀ MẶT MOHO THEO MẶT CẮT TÂY - ĐÔNG
  Độ sâu (km)
     0 km (Bờ biển)                      150 km                   350 km (Biển sâu)
  1. Hình thái bề mặt Moho chính xác cao:
    • Độ sâu Moho tại đới sát bờ và thềm nông Bắc Trung Bộ đạt giá trị lớn nhất: $28 - 31\text{ km}$ (đặc trưng cho vỏ lục địa điển hình).
    • Chiều sâu Moho giảm dần ra phía Đông: tại bồn trũng Hoàng Sa đạt $18 - 22\text{ km}$ (vỏ lục địa bị kéo mỏng / vỏ chuyển tiếp).
    • Tại rìa phía Đông Hoàng Sa giáp trũng sâu Biển Đông, Moho nâng cao mạnh mẽ, chỉ còn $12 - 14\text{ km}$, khẳng định sự xuất hiện của vỏ đại dương Đệ tam.
  2. Tính liên tục tự nhiên của cấu trúc địa chất: Bề dày trầm tích Kainozoi và cấu trúc móng trước Kainozoi trải dài liên tục từ đồng bằng duyên hải Bắc Miền Trung ra tới bể trầm tích Hoàng Sa và Phú Khánh mà không bị chia cắt bởi bất kỳ rãnh hút chìm (subduction trench) nào, chứng minh cấu trúc này là sự kéo dài tự nhiên không thể tách rời của lục địa Việt Nam.

Implications đa chiều

                   CÁC Ý NGHĨA VÀ GIÁ TRỊ ỨNG DỤNG CỦA CÔNG TRÌNH

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn kỹ thuật nội tại:

  1. Mật độ số liệu địa chấn sâu: Mặc dù đã bổ sung được 2 điểm khống chế Moho quan trọng từ tuyến PKGE09-09, mật độ điểm tựa địa chấn sâu vẫn còn thưa thớt trên toàn bộ diện tích hàng trăm nghìn kilômét vuông, buộc mô hình 3D phải sử dụng thuật toán nội suy trọng lực dư thừa mantle.
  2. Độ phân giải tầng móng trước Kainozoi: Sóng địa chấn bị suy giảm năng lượng mạnh qua các tập đá vôi và phun trào basalt Kainozoi, dẫn đến việc xác định ranh giới âm học móng ở một số vùng nước sâu còn độ bất định nhất định.
  3. Giả định 2D trong không gian thực tế 3D: Thuật toán NFG cải tiến trong luận án được xây dựng trên mô hình mặt cắt 2D, chưa tích hợp hoàn chỉnh toán tử NFG 3D tensor toàn phần.

Chương trình nghiên cứu 10 năm tiếp theo (2020 - 2030) cần tập trung:

  • Mở rộng thuật toán NFG và HGM cải tiến lên dạng Ten-xơ trọng lực gradient 3D (Full Tensor Gravity Gradiometry - FTG).
  • Khảo sát các tuyến địa chấn khúc xạ sâu dải rộng (OBS - Ocean Bottom Seismometer) cắt qua quần đảo Hoàng Sa.
  • Mô hình hóa trường nhiệt và đẳng tĩnh (isostasy) kết hợp đối lưu manti dưới lòng Biển Đông.

Tác động và ảnh hưởng

  • Học thuật: Thiết lập chuẩn mực mới trong quy trình minh giải phối hợp trọng lực - địa chấn tại Việt Nam. Các bài báo công bố từ luận án trên Tạp chí Khoa học ĐHQGHN và Tạp chí Các Khoa học về Trái Đất là tài liệu tham khảo cho các nghiên cứu sinh chuyên ngành Vật lý Địa cầu và Địa chất biển.
  • Ngành Dầu khí: Tập đoàn Dầu khí Việt Nam (Petrovietnam), PVEP và Vietsovpetro kế thừa trực tiếp bản đồ ranh giới móng và hệ thống đứt gãy sâu để đánh giá tiềm năng bẫy chứa phi truyền thống (dầu trong móng nứt nẻ granitoid Mesozoi) tại bể Phú Khánh và Sông Hồng.
  • An ninh chủ quyền và Chính sách biển: Đóng góp các tài liệu khoa học cho Ủy ban Biên giới Quốc gia (Bộ Ngoại giao) trong việc củng cố hồ sơ phân định ranh giới thềm lục địa vượt quá 200 hải lý trình Liên Hợp Quốc theo Điều 76 UNCLOS 1982.

Đối tượng hưởng lợi

                                  ĐỐI TƯỢNG HƯỞNG LỢI

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp độc đáo nhất là việc tái thiết lập toán tử tối ưu hóa số hài Fourier $N$ trong phương pháp Gradient chuẩn hóa toàn phần (NFG 2D). Thay vì áp dụng quy ước kinh nghiệm $N = 2M/3$ vốn tồn tại sai số phân kỳ lớn, tác giả đã tích hợp tiêu chuẩn Aydın (2007) để xác lập mối quan hệ hàm số giữa năng lượng phổ điều hòa cực đại với độ sâu tâm nguồn dị thường, giải quyết triệt để hiện tượng nhiễu tần số cao trong khai triển giải tích trường thế.

2. Sự cách tân phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với thuật toán HGM của Blakely & Simpson (1986) (chỉ tìm cực đại địa phương khi thỏa mãn đồng thời 2 điều kiện khắt khe trên 8 điểm nút lân cận), phương pháp HGM cải tiến của luận án đã mở rộng không gian tìm kiếm cực trị của đa thức nội suy bậc 2 trên toàn miền lân cận $[-d, d]$. Cải tiến này giúp khôi phục các đường phương đứt gãy liên tục tại các vùng chuyển tiếp có độ tương phản mật độ yếu mà thuật toán kinh điển bỏ sót.

3. Phát hiện bất ngờ nhất từ dữ liệu thực tế là gì?

Sự tồn tại của ranh giới Moho nâng cao bất đối xứng ($12 - 14\text{ km}$) tại rìa Đông bể Phú Khánh - Hoàng Sa tiếp giáp trũng sâu Biển Đông, đi kèm hệ thống đứt gãy sâu cắm về phía đại dương. Điều này chứng minh quá trình tách giãn vỏ Trái Đất tại khu vực diễn ra theo cơ chế tách giãn bất đối xứng (asymmetric rifting) với sự tham gia của các đứt gãy bóc tách vỏ sâu (detachment faults).

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp không?

Luận án cung cấp tường minh các tham số mô hình hóa: tọa độ các tuyến đo, bước lưới phân giải ($\Delta x, \Delta y$), các bảng thông số hình học lăng trụ 2D/3D (Bảng 2.1 - Bảng 2.8), bảng giá trị vận tốc - mật độ tại các giếng khoan (Bảng 3.1 - Bảng 3.2), cùng hệ phương trình toán học giải tích, cho phép các nhà khoa học độc lập lập trình và tái tạo toàn bộ kết quả tính toán.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo?

Tập trung vào phát triển phần mềm xử lý nghịch đảo tự động dữ liệu trọng trường 3D tích hợp cấu trúc ten-xơ gradient trọng lực, liên kết song song dữ liệu sóng địa chấn phân giải cao và số liệu từ trường biển để xây dựng mô hình địa động lực 4D (biến thiên theo thời gian) của toàn thềm lục địa Việt Nam.

Kết luận

Công trình của Tiến sĩ Trần Văn Khá mang lại 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Cải tiến thành công phương pháp NFG: Tối ưu hóa việc chọn số hài Fourier $N$, loại trừ sai số kinh nghiệm và xác định chuẩn xác độ sâu các ranh giới cấu trúc sâu.
  2. Hoàn thiện thuật toán HGM: Bổ sung hệ điều kiện toán học tìm cực trị gradient ngang, nâng cao khả năng vạch mạng lưới đứt gãy kiến tạo ẩn sâu thêm $30%$ độ liên tục.
  3. Xác lập hàm mật độ trầm tích Kainozoi chuẩn hóa: Cung cấp phương trình thực nghiệm $\rho(z)$ tin cậy cho toàn vùng thềm lục địa Bắc Miền Trung – Hoàng Sa.
  4. Chính xác hóa bản đồ địa hình bề mặt Moho: Xác định độ sâu Moho biến thiên từ $31\text{ km}$ (lục địa) xuống $12\text{ km}$ (biển sâu) dựa trên sự khống chế của 2 điểm tựa địa chấn sâu thực tế.
  5. Cung cấp bằng chứng khoa học vững chắc về sự liên tục tự nhiên của vỏ Trái Đất: Chứng minh sự gắn kết địa chất liền mạch từ đất liền ra vùng biển Hoàng Sa, phục vụ trực tiếp cho công tác quản lý biển đảo và xác lập ranh giới ngoài thềm lục địa theo luật pháp quốc tế.