Tổng quan về luận án

Nghiên cứu cấu trúc địa chất nông và trung sâu trên đất liền luôn là bài toán then chốt trong điều tra địa chất cơ bản, tìm kiếm khoáng sản ẩn sâu và đánh giá tai biến địa chất. Trong bối cảnh phương pháp địa chấn phản xạ 2D tại Việt Nam trước năm 2010 chủ yếu phục vụ thăm dò dầu khí ở bồn trũng Sông Hồng hoặc thềm lục địa với địa hình tương đối bằng phẳng, việc ứng dụng kỹ thuật này tại các khu vực đồi núi có địa hình phân cắt mạnh và lớp phủ phong hóa phức tạp đối mặt với rào cản kỹ thuật nghiêm trọng. Luận án tiến sĩ địa chất của nghiên cứu sinh Nguyễn Duy Bình (2019) với đề tài "Nghiên cứu một số đặc điểm cấu trúc địa chất trũng Sông Ba và vùng Đông Triều – Quảng Ninh theo tài liệu địa chấn phản xạ" (Mã số chuyên ngành: 9440201, Viện Khoa học Địa chất và Khoáng sản) đã giải quyết triệt để khoảng trống nghiên cứu (research gap) này. Khoảng trống cốt lõi được xác định là: (1) Thiếu vắng quy chuẩn tham số nguồn nổ tối ưu cho trạm ghi số kênh hạn chế trên đất liền; (2) Sự thất bại của các phương pháp hiệu chỉnh tĩnh truyền thống trước đới tốc độ thấp (low-velocity layer - LVL) biến thiên phi tuyến tính; (3) Cấu trúc sâu của đáy bồn trũng Neogen Sông Ba chưa được định lượng chính xác (các tài liệu lỗ khoan nước trước đây chỉ ước đoán chiều dày trung bình 350–400 m, trong khi các dự báo địa chấn nông phân giải cao của Dương Đức Kiêm năm 2006 lại dao động đến 1000 m mà chưa có bằng chứng trường sóng xác thực); (4) Cấu trúc móng và hệ thống đứt gãy khống chế các vỉa than sâu dưới −300 m đến −500 m tại bể than Đông Bắc còn nhiều điểm gián đoạn, giả định gượng ép.

Để giải quyết các khoảng trống trên, luận án đặt ra ba câu hỏi nghiên cứu và ba giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để tối ưu hóa trường sóng kích thích từ nguồn nổ trong hố khoan nhằm gia tăng tỷ số tín hiệu trên nhiễu ($S/N$) và bảo toàn thành phần tần số cao khi sử dụng thiết bị ghi địa chấn đa kênh cấu hình nông?
  • Câu hỏi nghiên cứu 2 (RQ2): Phương pháp hiệu chỉnh tĩnh nào có khả năng triệt tiêu triệt để ảnh hưởng của địa hình mấp mô và sự biến thiên dị thường của đới tốc độ thấp tại vùng đồi núi nhiệt đới?
  • Câu hỏi nghiên cứu 3 (RQ3): Hình thái cấu trúc thực tế của đáy trầm tích Neogen trũng Sông Ba và các hệ thống nếp uốn, đứt gãy sâu tại khối than Mạo Khê – Tràng Bạch (Đông Triều) được biểu hiện như thế nào trên mặt cắt thời gian và mặt cắt độ sâu địa chấn?

Hệ thống giả thuyết khoa học được thiết lập:

  • Giả thuyết 1 (H1): Việc định lượng bán kính vùng biến dạng dẻo $a = k \cdot Q^{1/3}$ kết hợp bố trí chiều sâu điểm nổ $h = 6\text{ m}$ và liều lượng thuốc nổ $Q = 300\text{ g}$ sẽ triệt tiêu sóng mặt (ground roll), sóng âm không khí, đồng thời mở rộng dải phổ tần số cao của sóng phản xạ.
  • Giả thuyết 2 (H2): Phương pháp giao thoa sóng khúc xạ (refraction interferometry statics) vượt trội hơn phương pháp thời gian tương hỗ cổ điển nhờ khả năng trích xuất chính xác mô hình vận tốc lớp mặt mà không bị phụ thuộc vào giả định phẳng của ranh giới khúc xạ.
  • Giả thuyết 3 (H3): Trũng Sông Ba tại khu vực Krông Pa là bồn trũng kiến tạo sụt lún dạng địa hào (graben) kéo tách bất đối xứng có chiều dày tầng Neogen sâu trên 800 m, trong khi địa tầng chứa than Đông Triều bị khống chế bởi hệ thống đứt gãy nghịch và trượt bằng phân khối kéo dài xuống quá mức −500 m.

Khung lý thuyết của nghiên cứu tích hợp ba nền tảng: Thuyết lan truyền sóng đàn hồi trong môi trường dị hướng phân lớp (Elastodynamics in stratified media), Thuyết địa chấn địa tầng học (Seismic Sequence Stratigraphy của Vail và Mitchum) và Thuyết bồn trũng trượt bằng kéo tách (Strike-slip pull-apart basin model). Phạm vi nghiên cứu bao quát tuyến địa chấn Ayun Pa, Krông Pa (Gia Lai) thu nổ bằng máy STRATAVISOR 48 kênh ($\Delta x = 5\text{ m}$, chiều dài chặng cáp 235 m, bội quan sát $N=12$) và tuyến đo Đông Triều (Quảng Ninh) thu nổ bằng trạm Sercel E428XL 410–480 kênh. Kết quả nghiên cứu mang ý nghĩa đột phá khi lần đầu tiên xác lập ranh giới móng Neogen Sông Ba ở độ sâu vượt 800 m, đồng thời chuẩn hóa quy trình công nghệ xử lý địa chấn 2D vùng đồi núi, phục vụ trực tiếp chiến lược điều tra khoáng sản ẩn sâu đến 1000 m của Chính phủ.

Literature Review và Positioning

Lịch sử nghiên cứu bồn trũng Sông Ba và bể than Đông Bắc ghi nhận những bước tiến dài nhưng tồn tại nhiều mâu thuẫn học thuật sâu sắc. Tại trũng Sông Ba, các công trình ban đầu của E. Saurin (1930–1959) và J. Fromaget (1941) chủ yếu khảo sát địa chất bề mặt và xác lập thang cổ sinh ban đầu. Sau năm 1975, Nguyễn Xuân Bao (1976–1981), Nguyễn Đức Dỹ (1980, 1987) và Trịnh Dánh (1982, 2000) đã phân chia hệ tầng Sông Ba ($N_1^3\text{sb}$) và hệ tầng Kon Tum ($N_1\text{kt}$), quy gán tuổi Miocen muộn dựa trên phức hệ bào tử phấn hoa và hóa đá thực vật (Laurus, Dipterocarpus). Ngược lại, Tống Duy Thanh, Vũ Khúc và Nguyễn Địch Dỹ (2005, 2006) trong công trình phân vị địa tầng Việt Nam lại xếp hệ tầng Sông Ba xuống tuổi Oligocen ($E_3$) và Kon Tum lên Miocen ($N_1$). Mâu thuẫn này kéo theo tranh luận về thời điểm khởi động cơ chế tách giãn tạo rift: liệu các đứt gãy trượt bằng bắt đầu từ Paleogen hay Miocen? Về mặt cấu trúc, các lỗ khoan nước nông và lỗ khoan thăm dò năng lượng ENRECA-1 (hợp tác giữa Viện Dầu khí Việt Nam và Sở Địa chất Đan Mạch - Greenland, GEUS, 2001–2004 sâu 500 m tại Krông Pa) mới chỉ bắt gặp 480 m trầm tích Neogen, khiến nhận thức về bề dày bồn trũng bị giới hạn dưới 500 m.

Tại vùng Đông Triều – Quảng Ninh, từ các nghiên cứu sơ khai của Saladin (1882) và Zeiller (1882) đến các đợt thăm dò quy mô lớn của Đovjikov (1965), Đỗ Chí Uy (1970) và Hoàng Văn Cân (1979), hệ tầng chứa than Hòn Gai ($T_3n\text{-}r\text{ hg}$) được phân chia thành 3 phụ hệ tầng ($hg_1, hg_2, hg_3$). Tuy nhiên, cấu trúc móng Paleozoi – Mesozoi sớm và sự phân cắt của hệ thống đứt gãy ngang ($F_A, F_{TB}, F_B, F_{129}, F_2, F_{11}, F_{15}$) dưới mức −300 m đến −500 m vẫn chủ yếu dựa trên phép ngoại suy hình học tuyến khoan thưa, dẫn đến sai số lớn khi tính toán trữ lượng than công nghiệp cấp 334a và 334b.

+-----------------------------------------------------------------------------+
|                          TIẾN TRÌNH NGHIÊN CỨU                              |
+-----------------------------------------------------------------------------+
| Giai đoạn trước 1975: E. Saurin, J. Fromaget -> Khảo sát bề mặt, định tuổi  |
|                                                                             |
| Giai đoạn 1975-2005: Nguyễn Xuân Bao, Trịnh Dánh, Đỗ Chí Uy -> Đo vẽ bản đồ |
|                      và khoan nông (<500 m), tranh cãi địa tầng Oligo-Miocen|
|                                                                             |
| Giai đoạn 2005-2019: Dương Đức Kiêm, ENRECA-1 -> Địa chấn nông, khoan sâu   |
|                      nhưng thiếu lời giải hiệu chỉnh tĩnh địa hình đồi núi  |
|                                                                             |
| LUẬN ÁN (2019):      Nguyễn Duy Bình -> Địa chấn phản xạ 2D + Giao thoa sóng|
|                      khúc xạ -> Xác lập đáy Neogen >800 m & khống chế đứt gãy|
+-----------------------------------------------------------------------------+

Về mặt phương pháp địa vật lý, các công trình quốc tế như nghiên cứu của Steeples và Miller (1990) hay Yilmaz (2001) đã chứng minh vai trò của địa chấn phản xạ nông phân giải cao trong nghiên cứu cấu trúc nội mảng. Tuy nhiên, việc áp dụng phương pháp thời gian tương hỗ (Generalized Reciprocal Method - GRM của Palmer, 1980) tại các khu vực đồi núi nhiệt đới có lớp phong hóa dày thường tạo ra các giả mạo cấu trúc do không giải quyết được biến thiên vận tốc ngang mạnh. Luận án tự định vị là công trình tiên phong tại Việt Nam kết hợp mô hình hóa trường sóng đàn hồi lý thuyết với phương pháp giao thoa sóng khúc xạ (Seismic Interferometry của Schuster, 2009) để hiệu chỉnh tĩnh, vượt qua các giới hạn của dự án ENRECA-1 và mở ra bước đột phá trong minh giải cấu trúc bồn trũng lục địa.

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

Đóng góp cho lý thuyết

Nghiên cứu mang lại những đóng góp mang tính bước ngoặt cho ngành khoa học địa chất và địa chấn học cấu trúc:

  1. Mở rộng thuyết bồn trũng kéo tách trượt bằng (Pull-apart basin theory): Cung cấp bằng chứng địa chấn định lượng chứng minh trũng Neogen Sông Ba không đơn thuần là bồn trũng nông dạng thung lũng sông uốn khúc mà là một bồn rift nội lục sâu, phát triển bất đối xứng dọc đới đứt gãy trượt bằng phương Tây Bắc – Đông Nam với độ sụt lún kiến tạo vượt 800 m tại phân võng Krông Pa.
  2. Xác lập mô hình quan hệ bất chỉnh hợp địa chấn học (Seismic sequence stratigraphy boundaries): Xác định và chuẩn hóa hệ thống 4 ranh giới phản xạ chính trên mặt cắt địa chấn phân giải cao:
    • Ranh giới $R_1$: Bất chỉnh hợp gá đáy (onlap) thể hiện quá trình tích tụ ban đầu phủ trên móng kết tinh tiền Neogen.
    • Ranh giới $R_2$ và $R_4$: Bất chỉnh hợp chống nóc (toplap) phản ánh các pha nâng cao kiến tạo, làm gián đoạn chu kỳ trầm tích.
    • Ranh giới $R_3$: Bất chỉnh hợp bao bọc (downlap) phản ánh tướng trầm tích hồ nước sâu lấp đầy bồn trũng.
Mô hình quan hệ địa tầng địa chấn bồn trũng Neogen:
----------------------------------------------------------------- Mặt đất
       / / / / / / / / / / / / / / / / / / / / / (Trầm tích Đệ Tứ Q)
================================================================= R4 (Toplap)
       ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ (Tập sét - diatomit Kon Tum)
----------------------------------------------------------------- R3 (Downlap)
       . . . . . . . . . . . . . . . . . . . . . (Tập cát-bột kết chứa than N1^3sb)
================================================================= R2 (Toplap)
       * * * * * * * * * * * * * * * * * * * * * (Tập cuội tảng cơ sở N1^3sb)
================================================================= R1 (Onlap)
       v v v v v v v v v v v v v v v v v v v v v (Móng Magma/Biến chất Pre-N)

Khung phân tích độc đáo

Khung phân tích của luận án được xây dựng dựa trên sự liên kết đa ngành chặt chẽ giữa Địa vật lý thực nghiệm – Toán mô hình hóa trường sóng – Địa chất trầm tích và khoáng sản. Phương pháp tiếp cận bao gồm 3 trụ cột tích hợp:

  • Trụ cột 1 - Mô phỏng số truyền sóng đàn hồi (Synthetic wave propagation modeling): Sử dụng các phương trình vi phân đàn hồi giải bằng sai phân hữu hạn để mô phỏng đáp ứng địa chấn trên các mặt cắt cấu trúc phức tạp tuyến Đông Triều (tuyến T-QN), so sánh hiệu quả thu nổ giữa các cấu hình 60 kênh, 120 kênh và 240 kênh nhằm dự báo trước các trường sóng giao thoa và nhiễu sóng phản xạ nhiều lần (multiples).
  • Trụ cột 2 - Công nghệ giao thoa sóng khúc xạ (Refraction wave interferometry): Xử lý hiệu chỉnh tĩnh bằng cách biến đổi sóng khúc xạ đầu tiên (first arrivals) thành các hàm giao thoa ảo, từ đó tính toán chính xác giá trị thời gian trễ (delay time) và mô hình vận tốc lớp vận tốc thấp (RVS, RCS) mà không cần khoan nổ kiểm tra sâu trong lớp phong hóa.
  • Trụ cột 3 - Minh giải địa chất tích hợp tài liệu địa chấn và giếng khoan sâu (Integrated seismic-well calibration): Liên kết trường sóng thời gian thực với cột địa tầng chuẩn của các lỗ khoan sâu (như LK.ENRECA-1, LK.XXI, LK.36) để chuyển đổi từ mặt cắt thời gian ($t_0$) sang mặt cắt độ sâu ($H$), xác định bản chất thạch học của các ranh giới sóng.

Điều kiện biên (Boundary conditions): Quy trình phân tích được thiết kế tối ưu cho địa hình đồi núi có độ dốc sườn $10^\circ \text{–} 50^\circ$, chiều sâu nghiên cứu mục tiêu từ 50 m đến 1000 m, vận tốc truyền sóng lớp phong hóa $V_0 = 300 \text{–} 1200\text{ m/s}$ và vận tốc tầng đá gốc $V_1 > 2200\text{ m/s}$.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng triết lý thực chứng hiện đại (post-positivism) thông qua thiết kế hỗn hợp giữa thực nghiệm tham số dã ngoại định lượng, mô hình hóa toán lý và phân tích thống kê trường sóng số liệu. Đối tượng nghiên cứu được chuẩn hóa chi tiết:

  • Khu vực trũng Sông Ba: Khảo sát trên 2 tuyến địa chấn 2D trọng điểm là tuyến T1 Ayun Pa và tuyến T2 Krông Pa.
  • Khu vực Đông Triều – Quảng Ninh: Khảo sát trên tuyến địa chấn T-QN cắt qua phức nếp lồi Mạo Khê – Tràng Bạch và nếp lõm Trung Lương.
+-----------------------------------------------------------------------------+
|                       QUY TRÌNH THỰC NGHIỆM VÀ XỬ LÝ                        |
+-----------------------------------------------------------------------------+
| THỰC NGHIỆM THU NỔ:  h = 0, 1.4, 3, 6 m | Q = 50, 100, 200, 300 g          |
|                      Chọn tối ưu: h = 6 m, Q = 300 g (STRATAVISOR / E428XL) |
|                                   v                                         |
| XỬ LÝ SỐ LIỆU:       Lọc Bandpass -> F-K Filter -> Deconvolution            |
|                                   v                                         |
| HIỆU CHỈNH TĨNH:     So sánh Phương pháp Thời gian Tương hỗ (GRM)           |
|                      với Phương pháp Giao thoa Sóng khúc xạ (Tối ưu)        |
|                                   v                                         |
| PHÂN TÍCH VẬN TỐC:   Vận tốc RMS (Semblance) -> NMO Correction -> CDP Stack  |
|                                   v                                         |
| MINH GIẢI ĐỊA CHẤT:  Hiệu chỉnh độ sâu -> Liên kết LK.ENRECA-1, LK.36       |
|                      Xác định đứt gãy FA, FB, F129 & Móng Neogen >800 m     |
+-----------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm hiện trường và thu nổ được tiến hành với độ chuẩn xác cao nhằm kiểm soát các biến số ngoại cảnh:

  1. Định lượng tham số nguồn nổ trong hố khoan: Quá trình nổ hình thành sóng đàn hồi phụ thuộc vào bán kính đới biến dạng dẻo $a = k \cdot Q^{1/3}$ (với $k$ là hệ số cơ lý đất đá, đá cứng $k = 1 \div 1.2$, đất dẻo $k = 1.2 \div 1.5$). Các thử nghiệm thực tế trên chặng máy 235 m với cửa sổ quan sát $X_{\min} = 120\text{ m}$ đã so sánh 4 mức độ sâu đặt quả nổ ($h = 0\text{ m}, 1.4\text{ m}, 3\text{ m}, 6\text{ m}$) và 4 mức liều lượng thuốc nổ ($Q = 50\text{ g}, 100\text{ g}, 200\text{ g}, 300\text{ g}$).
    • Kết quả thực nghiệm chiều sâu: Tại $h = 0\text{ m}$ và $1.4\text{ m}$, sóng mặt Rayleigh và sóng âm không khí chiếm ưu thế hoàn toàn, đè bẹp tín hiệu phản xạ. Tại $h = 6\text{ m}$ (đặt dưới đáy lớp phong hóa bở rời), sóng mặt bị triệt tiêu, sóng phản xạ hiện rõ nét với tần số ưu thế cao.
    • Kết quả thực nghiệm lượng nổ: Quả nổ $50\text{ g}$ và $100\text{ g}$ hạn chế nhiễu nhưng biên độ sóng phản xạ ở các kênh xa nguồn ($X_{\max} = 355\text{ m}$) rất yếu, không bảo đảm tỷ số $S/N$. Quả nổ $Q = 300\text{ g}$ tại $h = 6\text{ m}$ cho phổ biên độ rộng, năng lượng đồng đều trên toàn chặng máy 48 kênh.
  2. Kỹ thuật hiệu chỉnh tĩnh bằng giao thoa sóng khúc xạ: Trên các địa hình phân cắt mạnh ở Đông Triều và Sông Ba, lớp vận tốc thấp biến thiên từ vài mét đến hàng chục mét với vận tốc $V_0 = 400 \text{–} 800\text{ m/s}$. Luận án đã so sánh phương pháp thời gian tương hỗ truyền thống và phương pháp giao thoa sóng khúc xạ. Phương pháp giao thoa khúc xạ giải bài toán ngược ma trận thời gian truyền sóng đầu tiên để phân tích thời gian trễ tại từng điểm nổ và máy thu: $$\Delta t_{\text{static}} = -\left(\frac{h_{\text{LVL}}}{V_0} - \frac{h_{\text{LVL}}}{V_{\text{ref}}}\right)$$ Kết quả áp dụng trên tuyến Ayun Pa, Krông Pa và Đông Triều đã loại bỏ hoàn toàn các dạng nhảy pha "hình sin giả" do địa hình, làm thẳng các ranh giới sóng phản xạ sâu.

Data và phân tích

Quy trình xử lý số liệu địa chấn 2D được thực hiện trên hệ thống phần mềm chuyên dụng (ProMAX, Vista) với các bước phân tích nâng cao:

  • Chu trình tiền xử lý: Cắt đầu băng sóng (geometry assignment), lọc nhiễu ngẫu nhiên và nhiễu công nghiệp ($50\text{ Hz}$ notch filter), lọc dải thông (Bandpass filtering $15\text{–}25\text{–}80\text{–}120\text{ Hz}$), triệt tiêu sóng mặt bằng bộ lọc F-K (Frequency-Wavenumber domain).
  • Phân tích vận tốc và hiệu chỉnh động: Xây dựng mô hình phổ vận tốc (velocity semblance analysis) tại từng điểm giữa chung (CMP) với bước nhảy 10–20 CMP. Xác định chính xác trường vận tốc hiệu dụng RMS ($V_{\text{RMS}}$) từ $1800\text{ m/s}$ ở tầng nông lên đến $3800\text{–}4500\text{ m/s}$ ở tầng móng kết tinh.
  • Cộng sóng và dịch chuyển địa chấn: Áp dụng hiệu chỉnh động NMO (Normal Moveout), cộng tín hiệu điểm sâu chung (CDP Stack) với bội quan sát $N = 12$, tiếp nối bởi hiệu chỉnh tĩnh dư (Residual Statics) và dịch chuyển sau cộng (Post-stack Kirchhoff Migration) để đưa các mặt nghiêng và mặt trượt đứt gãy về đúng vị trí không gian thực tế.

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

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

Luận án mang lại 4 phát hiện địa chất – địa vật lý đột phá với độ tin cậy khoa học cao:

+-----------------------------------------------------------------------------+
|                      CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                      |
+-----------------------------------------------------------------------------+
| 1. ĐÁY TRẦM TÍCH NEOGEN SÔNG BA:                                            |
|    - Sâu >800 m tại Krông Pa (vượt xa mức 350-400 m trước đây).             |
|    - Cấu trúc địa hào sụt lún mạnh, chứa 2 phân tập trầm tích phân lớp rõ rệt.|
|                                                                             |
| 2. HỆ THỐNG RANH GIỚI BẤT CHỈNH HỢP ĐỊA CHẤN:                               |
|    - Xác lập R1 (gá đáy), R2 & R4 (chống nóc), R3 (bao bọc).                |
|    - Khống chế các thấu kính than nâu và diatomit không liên tục.          |
|                                                                             |
| 3. KHỐNG CHẾ ĐỨT GÃY SÂU BỂ THAN ĐÔNG TRIỀU:                                |
|    - Làm rõ đứt gãy nghịch FA (góc dốc 70-80 độ, phân chia khối Bắc - Nam). |
|    - Xác lập đứt gãy FTB, FB, F129 với cự ly dịch chuyển trượt 300-350 m.   |
|    - Làm sáng tỏ nếp lõm Trung Lương và nếp lồi Mạo Khê đến mức <-500 m.   |
|                                                                             |
| 4. BƯỚC NHẢY CÔNG NGHỆ HIỆU CHỈNH TĨNH:                                    |
|    - Phương pháp giao thoa sóng khúc xạ giảm sai số tĩnh từ >25 ms về <3 ms.|
+-----------------------------------------------------------------------------+
  1. Xác lập đáy trầm tích Neogen trũng Sông Ba sâu trên 800 m: Đây là phát hiện có ý nghĩa khoa học lớn nhất. Trên mặt cắt địa chấn tuyến Krông Pa, sóng phản xạ từ mặt móng ($R_1$) xuất hiện ở thời gian truyền sóng hai chiều $t_0 \approx 550 \text{–} 600\text{ ms}$, tương ứng độ sâu tuyệt đối vượt 800 m. Cấu trúc lòng bồn trũng có hình thái lồi lõm phức tạp, gồm 2 tập trầm tích: Tập dưới dày trên 450 m gồm cuội kết, sạn kết xen cát bột kết gắn kết yếu; Tập trên dày 350 m gồm cát kết, sét kết diatomit xen kẽ các vỉa than nâu không liên tục đặc trưng bởi các pha sóng đứt đoạn, biên độ thay đổi mạnh.
  2. Khám phá và chuẩn hóa hệ thống ranh giới bất chỉnh hợp địa tầng địa chấn ($R_1 \div R_4$): Liên kết trực tiếp mặt cắt địa chấn Ayun Pa – Krông Pa với tài liệu lỗ khoan LK.ENRECA-1 đã giải mã nguồn gốc các pha sóng: ranh giới $R_4$ là bất chỉnh hợp bào mòn – chống nóc ngăn cách Neogen với Đệ Tứ; $R_3$ là ranh giới bao bọc thể hiện pha ngập lụt lòng hồ cực đại; $R_2$ đánh dấu pha tái cấu trúc kiến tạo giữa Miocen; và $R_1$ là ranh giới móng magma – biến chất Tiền Cambri / Mesozoi.
  3. Chính xác hóa cấu trúc sâu và hệ thống đứt gãy bể than Đông Triều – Quảng Ninh: Tài liệu địa chấn phản xạ 2D trên tuyến Đông Triều đã khống chế rõ ràng cấu trúc phức nếp lồi Mạo Khê – Tràng Bạch và nếp lõm Trung Lương xuống sâu quá −500 m. Đặc biệt, luận án đã làm rõ hình thái hình học của các đứt gãy kiến tạo lớn:
    • Đứt gãy nghịch $F_A$: Cắm Bắc với góc dốc $70^\circ \text{–} 80^\circ$, đới phá hủy rộng 50–100 m, đóng vai trò đứt gãy phân khối chia mỏ Mạo Khê thành khối Bắc và khối Nam.
    • Đứt gãy $F_{TB}$ và $F_{129}$: Là các đứt gãy thuận dịch chuyển ngang với cự ly trượt theo mặt trượt từ 300 m đến 350 m, cự ly dịch chuyển ngang lên tới 250–290 m, làm dịch chuyển và uốn vặn chùm vỉa than 1(36) và 1B(35).
    • Đứt gãy $F_B$: Xác định chính xác ranh giới tiếp xúc giữa phụ hệ tầng Hòn Gai dưới ($hg_1$) và Hòn Gai giữa ($hg_2$), giải quyết dứt điểm các giả định thiếu cơ sở trước đây tại chân đập Yên Dưỡng (khu vực LK.XXI và LK.36).
  4. Chứng minh tính ưu việt tuyệt đối của phương pháp giao thoa sóng khúc xạ: Trên các mặt cắt xử lý so sánh, hiệu chỉnh tĩnh bằng giao thoa sóng khúc xạ đã làm giảm sai số thăng giáng thời gian từ $>25\text{ ms}$ (phương pháp tương hỗ) xuống dưới $3\text{ ms}$, khôi phục hoàn hảo tính liên tục của các phản xạ vỉa than nằm nghiêng $20^\circ \text{–} 50^\circ$.

Implications đa chiều

  • Ý nghĩa lý thuyết: Cung cấp cơ sở định lượng để tái đánh giá lịch sử kiến tạo Kainozoi khu vực Tây Nguyên, khẳng định đới đứt gãy Sông Ba trải qua các pha tách giãn mãnh liệt tạo địa hào sâu trước khi chuyển sang pha nén ép trượt bằng hiện đại.
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình đo đạc và xử lý địa chấn phản xạ 2D cho vùng địa hình đồi núi nhiệt đới, có thể chuyển giao áp dụng cho toàn bộ các bồn than, bồn trũng đứt gãy tại miền Bắc và miền Trung Việt Nam.
  • Ứng dụng thực tiễn và chuyển đổi ngành: Giúp Tập đoàn Than – Khoáng sản Việt Nam (TKV) và Tổng cục Địa chất định vị chính xác đáy tầng chứa than và khoanh vùng các khối cấu trúc an toàn, giảm thiểu rủi ro bục nước ngầm và sập lở lò khi khai thác than hầm lò mức −300 m đến −500 m.
  • Khuyến nghị chính sách: Đóng góp luận cứ khoa học then chốt phục vụ Đề án "Điều tra, đánh giá tổng thể tài nguyên than phần đất liền bể Sông Hồng và bể than Đông Bắc" theo Quyết định của Thủ tướng Chính phủ.

Limitations và Future Research

Luận án thẳng thắn nhìn nhận các giới hạn kỹ thuật không thể tránh khỏi do điều kiện thi công tại thời điểm nghiên cứu:

  1. Mạng lưới tuyến đo 2D đơn lẻ: Do kinh phí và địa hình rừng núi hiểm trở, nghiên cứu mới chỉ thực hiện các tuyến địa chấn 2D riêng biệt (tuyến Krông Pa, Ayun Pa, Đông Triều) mà chưa thể xây dựng mạng lưới không gian 3D, dẫn đến hạn chế trong việc theo dõi liên tục sự biến đổi phương vị của các đứt gãy trượt bằng phân nhánh.
  2. Cấu hình trạm máy đo tại Sông Ba còn hạn chế: Việc sử dụng máy ghi 48 kênh với chiều dài chặng máy 235 m tại trũng Sông Ba khiến bội quan sát chỉ đạt $N = 12$. Mức bội này làm hạn chế khả năng triệt tiêu triệt để sóng phản xạ nhiều lần ở độ sâu lớn hơn 1000 m.
  3. Mật độ giếng khoan khống chế sâu chưa đồng đều: Tại một số đoạn tuyến ở Đông Triều, do thiếu các lỗ khoan sâu dưới −500 m nên việc liên kết đồng danh một số vỉa than mỏng ở cánh nếp uốn vẫn phải dựa vào phép nội suy vận tốc sóng.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển khảo sát địa chấn phản xạ 3D phân giải cao tại các khu vực mỏ trọng điểm để xây dựng mô hình khối cấu trúc không gian 3 chiều.
  • Ứng dụng kỹ thuật Nghịch đảo toàn dạng sóng (Full Waveform Inversion - FWI) và Địa chấn thụ động giao thoa (Passive Ambient Noise Tomography - ANT) để nghiên cứu cấu trúc mỏ sâu mà không cần dùng nguồn nổ hóa học.
  • Khoan kiểm chứng cấu trúc đáy bồn trũng Sông Ba ở độ sâu 800–1000 m để lấy mẫu lõi đáy phục vụ nghiên cứu địa hóa đồng vị và phân tích tướng trầm tích chi tiết.

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

Luận án của NCS Nguyễn Duy Bình đã tạo nên những tác động sâu rộng trên cả phương diện học thuật và kinh tế – xã hội:

  • Tác động học thuật: Mở ra trào lưu ứng dụng địa vật lý địa chấn nông phân giải cao trong nghiên cứu cấu trúc địa chất đất liền tại Việt Nam, trở thành tài liệu tham khảo kinh điển cho các nghiên cứu sinh chuyên ngành Địa chất học và Địa vật lý.
  • Chuyển đổi công nghiệp khai khoáng: Kết quả xác định vị trí các đứt gãy phân khối $F_A, F_{TB}, F_{129}$ đã được các công ty than thuộc TKV (Công ty Than Mạo Khê, Than Tràng Bạch, Than Uông Bí) đưa vào cơ sở dữ liệu địa chất để tối ưu hóa thiết kế mạng lưới đường lò và hạn chế khoan thăm dò lãng phí.
  • Giá trị xã hội và môi trường: Cung cấp thông tin địa chất công trình và đới đứt gãy hoạt động, hỗ trợ đánh giá nguy cơ trượt lở đất và tai biến địa chất phục vụ xây dựng cơ sở hạ tầng bền vững dọc lưu vực sông Ba và vùng kinh tế trọng điểm Quảng Ninh.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận quy trình thực nghiệm tham số nguồn nổ, công thức xác định bán kính vùng phá hủy $a = k \cdot Q^{1/3}$ và kỹ thuật hiệu chỉnh tĩnh bằng giao thoa sóng khúc xạ.
  • Các nhà khoa học địa chất và địa vật lý: Kế thừa tập dữ liệu sóng địa chấn thực tế, mô hình vận tốc RMS và các phát hiện cấu trúc bồn trũng Neogen >800 m để tiếp tục giải quyết bài toán tiến hóa địa động lực Kainozoi Đông Nam Á.
  • Kỹ sư thăm dò mỏ và doanh nghiệp khai khoáng (TKV, PVN): Ứng dụng trực tiếp bản đồ đứt gãy và mặt cắt địa chấn độ sâu để thiết kế khai thác mỏ an toàn và chuẩn xác hóa trữ lượng than mức sâu.
  • Cơ quan quản lý nhà nước (Tổng cục Địa chất và Khoáng sản, Bộ Tài nguyên và Môi trường): Sử dụng kết quả luận án làm cơ sở khoa học phê duyệt các đề án quy hoạch thăm dò tài nguyên khoáng sản quốc gia đến độ sâu 1000 m.

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ì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc chứng minh và xác lập mô hình cấu trúc bồn trũng rift nội lục kéo tách sâu trên 800 m tại Krông Pa, mở rộng Thuyết bồn trũng trượt bằng (Pull-apart basin theory) của Mann và Thuyết địa chấn địa tầng của Vail. Nghiên cứu đã chứng minh sự phát triển đa chu kỳ với 4 ranh giới bất chỉnh hợp địa chấn ($R_1 \div R_4$), giải quyết dứt điểm mâu thuẫn học thuật kéo dài hàng thập kỷ về chiều dày thực tế của hệ tầng Sông Ba.

2. Đổi mới phương pháp luận then chốt so với các nghiên cứu trước đây là gì? Trả lời: So với các nghiên cứu trước đây chỉ dùng phương pháp thời gian tương hỗ (GRM của Palmer) hoặc hiệu chỉnh cao trình đơn giản (Field Statics), luận án đã tiên phong ứng dụng thành công phương pháp Hiệu chỉnh tĩnh bằng giao thoa sóng khúc xạ (Refraction Wave Interferometry). Phương pháp này xử lý chính xác sự bất đồng nhất phi tuyến của đới tốc độ thấp tại vùng đồi núi, khôi phục trung thực trường sóng phản xạ.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng số liệu chứng minh là gì? Trả lời: Phát hiện bất ngờ nhất là chiều sâu móng Neogen tại Krông Pa vượt trên 800 m (trước đây các tài liệu khoan và đo vẽ địa chất chỉ ghi nhận 350–400 m, lỗ khoan ENRECA-1 chỉ sâu 500 m). Bằng chứng thực tế là pha sóng phản xạ $R_1$ xuất hiện rõ nét trên băng địa chấn tại thời gian $t_0 \approx 550 \text{–} 600\text{ ms}$ với vận tốc trung bình lớp phủ $V_{\text{RMS}} \approx 2800\text{ m/s}$.

4. Luận án có cung cấp quy trình nhân bản (Replication protocol) cho các khu vực khác không? Trả lời: Có. Luận án cung cấp bộ quy trình chi tiết gồm: (1) Công thức lựa chọn lượng nổ $Q$ và độ sâu hố nổ $h = 6\text{ m}$; (2) Thiết kế hệ thống quan sát cánh có cửa sổ ($X_{\min}=120\text{ m}, X_{\max}=355\text{ m}, \Delta x=5\text{ m}$); (3) Trình tự 6 bước xử lý tín hiệu và quy trình hiệu chỉnh tĩnh giao thoa sóng khúc xạ chuẩn hóa.

5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Trả lời: Định hướng tập trung vào 3 mũi nhọn: Nâng cấp mạng lưới khảo sát từ 2D lên 3D phân giải cao; Kết hợp địa chấn phản xạ với địa chấn thụ động (Ambient Noise Tomography); Thực hiện các lỗ khoan sâu trên 1000 m tại tâm trũng Krông Pa để nghiên cứu hệ thống dầu khí và khoáng sản ẩn sâu.

Kết luận

Công trình nghiên cứu tiến sĩ của tác giả Nguyễn Duy Bình là một dấu mốc học thuật xuất sắc, kết hợp nhuần nhuyễn giữa địa vật lý ứng dụng và địa chất cấu trúc. Những đóng góp cốt lõi của luận án được đúc kết:

  1. Khẳng định quy mô cấu trúc bồn trũng Neogen Sông Ba: Lần đầu tiên chứng minh bằng tài liệu địa chấn phản xạ chiều sâu đáy bồn trũng tại Krông Pa đạt trên 800 m, làm thay đổi căn bản nhận thức địa chất trầm tích luận Kainozoi Tây Nguyên.
  2. Đột phá công nghệ xử lý hiệu chỉnh tĩnh: Chuẩn hóa phương pháp giao thoa sóng khúc xạ trong xử lý tài liệu địa chấn phản xạ 2D vùng đồi núi, triệt tiêu sai số do lớp vận tốc thấp và địa hình chia cắt phức tạp.
  3. Tối ưu hóa quy trình thu nổ thực địa: Xác lập bộ tham số nổ tối ưu ($h = 6\text{ m}, Q = 300\text{ g}$) cho máy ghi đa kênh trên đất liền, nâng cao vượt bậc tỷ số tín hiệu trên nhiễu ($S/N$).
  4. Làm sáng tỏ cấu trúc sâu bể than Đông Triều – Quảng Ninh: Khống chế chính xác hệ thống nếp uốn và các đứt gãy phân khối lớn ($F_A, F_{TB}, F_B, F_{129}$) đến độ sâu dưới −500 m, phục vụ đắc lực cho công tác thăm dò than an toàn.
  5. Khai phóng các hướng nghiên cứu liên ngành mới: Mở ra các hướng nghiên cứu triển vọng về đánh giá tiềm năng bentonite, than nâu, bẫy cấu trúc chứa nước dưới đất và thăm dò khoáng sản ẩn sâu đến độ sâu 1000 m trên phạm vi toàn quốc.