Tổng quan nghiên cứu

Khảo sát địa vật lý biển đóng vai trò đặc biệt quan trọng trong việc làm sáng tỏ cấu trúc vỏ Trái Đất, phục vụ đánh giá địa chất công trình và tìm kiếm tài nguyên khoáng sản. Tại Việt Nam, phương pháp Địa chấn nông phân giải cao (ĐCNPGC) đã được triển khai từ đầu những năm 1990 với tổng khối lượng đo đạc lũy kế đến nay đạt trên 40.000 km tuyến khảo sát. Riêng trong đề án điều tra khoáng sản rắn vùng biển nông ven bờ giai đoạn 1991–2001, các nhà khoa học đã thực hiện thành công 25.148 km tuyến địa chấn. Tuy nhiên, các tài liệu ghi nhận trong môi trường nước nông ven bờ (độ sâu từ 0 đến 30 m) thường xuyên bị ảnh hưởng bởi nhiều nguồn nhiễu phức tạp như sóng mặt nước, dao động vi địa chấn, nhiễu công nghiệp tần số 50–60 Hz và đặc biệt là nhiễu phản xạ nhiều lần (PXNL) từ bề mặt đáy biển. Những yếu tố này làm suy giảm nghiêm trọng tỷ số tín hiệu trên nhiễu (SNR) và che lấp các ranh giới địa chất có ích ở độ sâu dưới 50 m.

Mục tiêu cốt lõi của nghiên cứu là xây dựng và kiểm chứng tổ hợp quy trình xử lý số liệu địa chấn đơn kênh trên phần mềm chuyên dụng Reflexw, nhằm khôi phục trường sóng phản xạ thực, triệt tiêu phông nhiễu môi trường và nâng cao độ phân giải thẳng đứng xuống dưới mức 1 m. Phạm vi nghiên cứu được thực hiện tại dải bờ biển Cửa Đáy thuộc khu vực giáp ranh giữa tỉnh Ninh Bình và Nam Định, giới hạn bởi tọa độ địa lý từ 106°00' đến 106°20' kinh độ Đông và 19°40' đến 20°00' vĩ độ Bắc. Công trình có ý nghĩa thực tiễn to lớn trong việc cung cấp bức tranh địa tầng Đệ Tứ rõ nét, hỗ trợ trực tiếp cho các dự án nạo vét luồng hàng hải, mở rộng hệ thống cảng biển nước sâu và quy hoạch công trình ven bờ khu vực vịnh Bắc Bộ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết truyền sóng đàn hồi trong môi trường chất lỏng và trầm tích bở rời ngậm nước, bao gồm ba trụ cột chính:

  1. Quy luật hấp thụ và suy giảm năng lượng sóng địa chấn: Khi xung sóng lan truyền qua môi trường trầm tích tầng nông, năng lượng bị tiêu hao mạnh dưới dạng nhiệt lượng theo hàm số mũ $A = A_0 e^{-\alpha z} = A_0 e^{-\rho z / \lambda}$. Với hệ số hấp thụ trung bình $\rho = 0,05$ của trầm tích bở rời, tại độ sâu $z = 170\text{ m}$ (tương đương 100 lần bước sóng ở tần số 1000 Hz và vận tốc 1700 m/s), độ suy giảm biên độ tín hiệu đạt mức $-43,4\text{ dB}$, đòi hỏi phải áp dụng các hàm phục hồi biên độ theo thời gian.
  2. Lý thuyết độ phân giải địa chấn: Độ phân giải thẳng đứng được xác định theo tiêu chuẩn Rayleigh $\Delta h = \lambda / 4 = v / (4f)$. Với tần số chủ đạo $f = 400\text{ Hz}$ và vận tốc truyền sóng $v = 1800\text{ m/s}$, khoảng cách ranh giới nhỏ nhất có thể tách biệt là $1,125\text{ m}$. Độ phân giải ngang phụ thuộc vào bán kính đới Fresnel thứ nhất $r = \sqrt{\frac{v \cdot h}{2f} + \frac{v^2}{16f^2}}$, đòi hỏi tối ưu hóa tần số nguồn phát và khống chế tốc độ di chuyển của tàu để các đới Fresnel gối sát nhau.
  3. Cơ chế hình thành sóng phản xạ nhiều lần (PXNL) và lọc ngược tiên đoán: Do sự chênh lệch âm trở rất lớn giữa môi trường nước biển ($v \approx 1500\text{ m/s}$) và đáy biển ($v \approx 1700 - 1800\text{ m/s}$), mặt nước biển hoạt động như một ranh giới phản xạ hoàn hảo (hệ số phản xạ $K \approx 1$). Sóng PXNL xuất hiện có tính quy luật tại các thời điểm trễ $t + \alpha, t + 2\alpha$, được mô hình hóa và loại trừ qua phương trình tích phân Wiener-Hopf trên miền thời gian - khoảng cách ($T$-$X$).

Các khái niệm then chốt được vận dụng xuyên suốt bao gồm: Địa chấn nông phân giải cao (ĐCNPGC), tỷ số tín hiệu trên nhiễu (SNR), đới Fresnel thứ nhất, nhiễu ngẫu nhiên, nhiễu đồng bộ và phép lọc thích nghi.

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

Tập dữ liệu địa chấn thực địa được thu thập dọc theo các tuyến khảo sát vùng biển Cửa Đáy, tiêu biểu là các tuyến CUADAY-02B và CUADAY-03. Cỡ mẫu khảo sát bao gồm hàng nghìn đường ghi địa chấn liên tục được số hóa chi tiết:

  • Nguồn phát sóng âm: Sử dụng đầu phát Sparker chùm 50 điện cực dạng xương cá, kích nổ với năng lượng 500 J từ bộ tích năng lượng Geontself 97 (công suất cực đại 1000 J), tạo ra dải tần số làm việc hiệu dụng từ 200 Hz đến 1000 Hz.
  • Hệ thống thu nhận: Dải ống thu Hydrophone Benthos gồm 10 phần tử áp điện triệt tiêu gia tốc, chiều dài 12 m, kéo sau đuôi tàu cách nguồn phát từ 4 đến 8 m. Khoảng cách giữa các máy thu được thiết lập chính xác 1 m, đáp ứng định lý lấy mẫu Nyquist cho tần số cao nhất 750 Hz với vận tốc nước biển 1500 m/s.
  • Tham số đo ghi: Chiều dài băng ghi từ 100 ms đến 200 ms, chu kỳ lấy mẫu 1 s, định vị vệ tinh bằng hệ thống Anten DGPS Trimble DSM 1212H cho sai số tọa độ từ 3 đến 4 m ở tỷ lệ bản đồ 1/50.000. Tốc độ hành trình của tàu đo được giữ ổn định ở mức 6–8 hải lý/giờ trong điều kiện sóng biển dưới cấp 4.

Phương pháp phân tích được tiến hành trên phần mềm Reflexw thông qua chuỗi xử lý 5 bước: (1) Tiền xử lý loại bỏ các đường ghi hỏng (trace editing); (2) Phục hồi biên độ bằng hàm khuếch đại mũ Gain Function; (3) Trung bình hóa 5 đường ghi liền kề (running average); (4) Lọc dải thông Bandpass Filter từ 100 Hz đến 1000 Hz; (5) Triệt tiêu nhiễu sóng lặp đáy biển bậc 2 bằng bộ lọc sai số tiên đoán (predictive deconvolution). Lý do lựa chọn chuỗi xử lý này là khả năng tác động chọn lọc trên cả miền tần số lẫn miền $T$-$X$, bảo toàn tối đa tính trung thực của các pha phản xạ sơ cấp.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình áp dụng tổ hợp thuật toán trên mặt cắt địa chấn Cửa Đáy mang lại các kết quả định lượng cụ thể:

  1. Khôi phục năng lượng sóng tầng sâu: Việc áp dụng hàm phục hồi biên độ dạng mũ $g(t) = (1 + at)e^{bt}$ với hệ số lựa chọn $a = 0$ và $b = 0,1$ đã bù đắp hiệu quả phần năng lượng bị tiêu hao theo thời gian truyền sóng hai chiều. Trên mặt cắt CUADAY-03, các phản xạ địa chấn nằm ở khoảng thời gian từ 60 ms trở xuống (độ sâu tương đương trên 50 m) trước đây bị mờ nhạt đã được làm rõ nét, biên độ tín hiệu có ích được tăng cường ước tính khoảng 35% đến 40%.
  2. Tăng cường tính liên tục của ranh giới địa tầng: Thuật toán trung bình hóa 5 đường ghi liền kề (Running Average) đã phát huy tối đa hiệu ứng cộng đồng pha. Các dao động ngẫu nhiên không tương quan bị triệt tiêu lẫn nhau, trong khi các ranh giới phản xạ phẳng và nghiêng tăng tính liên tục theo phương ngang thêm khoảng 30%, khắc phục hiện tượng đứt gãy giả do chập chờn tín hiệu cục bộ.
  3. Thu hẹp phổ tần số và loại bỏ phông nhiễu ngoài dải: Phân tích phổ tần số trước xử lý cho thấy dải tần số bị phân tán rất rộng từ 50 Hz đến 3500 Hz với lượng lớn nhiễu thủy động học (0–50 Hz) và nhiễu điện lưới (50–60 Hz). Sau khi áp dụng bộ lọc dải dạng cos-taper trong khoảng 100–1000 Hz, toàn bộ năng lượng nhiễu tần số thấp và nhiễu cao tần trên 1000 Hz đã bị loại bỏ hoàn toàn, nâng tỷ số SNR tổng thể lên hơn 50%.
  4. Triệt tiêu sóng phản xạ nhiều lần đáy biển: Bằng phương pháp chọn pha (picking) chính xác thời gian xuất hiện sóng phản xạ đáy biển bậc một tại vị trí 24 ms, bộ lọc ngược tiên đoán với tham số trễ thời gian $\text{lag} \le 1,0\text{ ms}$ đã loại bỏ thành công pha sóng PXNL bậc hai xuất hiện tại mốc thời gian 48 ms trên nhóm đường ghi từ 1800 đến 1860 của mặt cắt CUADAY-03.

Thảo luận kết quả

Sự suy giảm chất lượng dữ liệu địa chấn biển nông tại Cửa Đáy bắt nguồn từ đặc điểm địa chất trầm tích Đệ Tứ chứa nhiều bùn sét ngậm nước và cát mịn, kết hợp với địa hình đáy biển thoải gây ra hiện tượng cộng hưởng âm học mạnh trong tầng nước nông từ 10 m đến 30 m. Nếu chỉ áp dụng các bộ lọc tần số 1D đơn lẻ, hệ thống hoàn toàn bất lực trước sóng PXNL đáy biển vì sóng này mang cùng dải phổ tần số (200–1000 Hz) như sóng phản xạ thực từ các tầng sâu.

So sánh với các nghiên cứu địa chấn phản xạ tại vịnh Thái Lan hoặc vùng biển Thổ Chu, việc tích hợp lọc ngược tiên đoán trên miền $T$-$X$ kết hợp lọc dải thông miền tần số là hướng tiếp cận ưu việt nhất. Kết quả minh giải trên các mặt cắt sau xử lý cho phép bóc tách rõ ràng 3 đến 4 phức nếp trầm tích Holocene và Pleistocene, xác định chính xác bề mặt bào mòn cổ và các đới nứt nẻ quy mô nhỏ có biên độ dịch chuyển chỉ từ 1 đến 2 m. Các dữ liệu này có thể được trình bày trực quan thông qua các bảng so sánh chỉ số suy giảm năng lượng theo decibel và các biểu đồ phân bố mật độ phổ biên độ trước - sau xử lý để minh chứng cho độ tin cậy của thuật toán.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu thực nghiệm, bốn nhóm khuyến nghị hành động cụ thể được đề xuất nhằm nâng cao hiệu quả công tác khảo sát địa vật lý biển:

  1. Chuẩn hóa quy trình xử lý 5 bước trên hệ thống phần mềm chuyên dụng: Viện Địa chất và Địa vật lý biển cùng các đơn vị khảo sát địa chất cần ban hành hướng dẫn kỹ thuật chuẩn hóa quy trình 5 bước (Editing - Gain - Running Average - Bandpass - Predictive Deconvolution). Mục tiêu đặt ra là nâng cao tỷ số SNR thêm ít nhất 40% cho toàn bộ các tập dữ liệu địa chấn biển nông đo mới trong giai đoạn 2024–2026.
  2. Nâng cấp cấu hình thiết bị phát thu ngoài thực địa: Các liên đoàn địa chất biển cần chuyển đổi sang sử dụng nguồn phát Sparker đa cực từ 50 đến 100 điện cực, duy trì công suất phát ổn định 500 J kết hợp hệ thống định vị DGPS có sai số dưới 1–2 m. Lộ trình thực hiện cần hoàn thành trước quý IV năm 2026 do Bộ Tài nguyên và Môi trường chỉ đạo.
  3. Tự động hóa khâu chọn pha và tính toán tham số lọc ngược: Các nhóm nghiên cứu tại Trường Đại học Khoa học Tự nhiên cần đẩy mạnh ứng dụng thuật toán học máy (Machine Learning) để tự động hóa khâu bắt pha ranh giới đáy biển (auto-picking) và tối ưu hóa hệ số trễ $\text{lag}$, giúp rút ngắn 50% thời gian xử lý thủ công các mặt cắt địa chấn dài hàng trăm kilomet từ năm 2025.
  4. Xây dựng cơ sở dữ liệu số hóa địa chấn biển nông quốc gia: Tổng cục Biển và Hải đảo chủ trì tập hợp, xử lý lại và chuẩn hóa hơn 40.000 km băng ghi địa chấn lịch sử theo định dạng chuẩn SEG-Y quốc tế, tích hợp bản đồ vận tốc âm học trầm tích đáy biển trên toàn thềm lục địa Việt Nam với mốc hoàn thành vào năm 2027.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu chuyên khảo mang giá trị thực tiễn và học thuật cao cho các nhóm đối tượng:

  1. Kỹ sư và chuyên gia xử lý số liệu địa vật lý: Nắm vững phương pháp thiết lập tham số lọc dải cos-taper và hàm Gain Function trên phần mềm Reflexw để giải quyết bài toán triệt nhiễu trên các băng ghi thực địa phức tạp.
  2. Cơ quan quản lý và quy hoạch công trình cảng - hàng hải: Ứng dụng các lát cắt địa chấn phân giải cao tại khu vực Cửa Đáy để đánh giá hiện trạng bồi lắng luồng tàu, thiết kế nạo vét tuyến giao thông thủy và quy hoạch cụm cảng Ninh Phúc, Tiên Hưng.
  3. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Vật lý Địa cầu, Địa chất biển: Sử dụng làm tài liệu giảng dạy và nghiên cứu chuyên sâu về trường sóng đàn hồi, lý thuyết đới Fresnel và kỹ thuật địa chấn đơn kênh vùng nước nông.
  4. Các doanh nghiệp tư vấn khảo sát địa kỹ thuật công trình ven biển: Tham khảo cấu hình đo đạc nguồn phát Sparker 500 J và Hydrophone 10 phần tử để khảo sát nền móng cho các dự án điện gió ngoài khơi, đê kè chắn sóng và đường ống ngầm vượt biển.

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

Tại sao phương pháp địa chấn nông phân giải cao (ĐCNPGC) lại thích hợp cho vùng biển nông ven bờ? Phương pháp ĐCNPGC sử dụng nguồn phát tần số cao (200–1000 Hz đối với Sparker và lên đến 10.000 Hz đối với Boomer), tạo ra bước sóng ngắn giúp đới Fresnel thứ nhất hẹp lại. Điều này mang lại độ phân giải đứng vượt trội tới vài chục centimet, cho phép phân tách chi tiết các lớp trầm tích mỏng ở độ sâu dưới 50 m mà các hệ thống địa chấn dầu khí sâu không thể thực hiện được.

Nhiễu phản xạ nhiều lần (PXNL) đáy biển phát sinh do đâu và gây tác hại gì? Nhiễu PXNL phát sinh do sự tương phản âm trở rất lớn giữa tầng nước biển ($1500\text{ m/s}$) và lớp trầm tích đáy biển ($1700\text{ m/s}$). Năng lượng sóng bị phản xạ qua lại nhiều lần giữa đáy biển và mặt nước, tạo ra các ranh giới giả tại các mốc thời gian trễ gấp đôi (như 48 ms so với đáy thật ở 24 ms), làm che lấp hoàn toàn các cấu trúc địa chất thực ở tầng sâu.

Vì sao không thể dùng bộ lọc tần số thông thường để triệt tiêu sóng phản xạ nhiều lần? Bản chất của sóng PXNL là sóng phản xạ lặp lại từ chính ranh giới đáy biển nên chúng mang đặc trưng phổ tần số (100–1000 Hz) và biên độ hoàn toàn trùng khớp với sóng phản xạ sơ cấp có ích. Do đó, các bộ lọc dải tần số không thể phân tách được chúng, bắt buộc phải dùng bộ lọc ngược tiên đoán trên miền $T$-$X$ dựa trên tính chu kỳ thời gian.

Hàm phục hồi biên độ $g(t) = e^{0,1t}$ mang lại hiệu quả gì cho mặt cắt CUADAY-03? Khi sóng lan truyền trong trầm tích bở rời, năng lượng bị suy giảm tới $-43,4\text{ dB}$ ở độ sâu 170 m do hiện tượng hấp thụ. Áp dụng hàm khuếch đại mũ với hệ số tối ưu $b = 0,1$ giúp bù trừ chính xác lượng năng lượng tiêu hao theo quy luật hình học và hấp thụ nhiệt, làm nổi bật rõ các ranh giới địa tầng sâu từ 60 ms trở xuống mà không làm bùng phát nhiễu nền.

Công tác thu nổ địa chấn Sparker ngoài biển cần tuân thủ các điều kiện kỹ thuật nào? Khảo sát thực địa cần duy trì tốc độ tàu chạy ổn định từ 6 đến 8 hải lý/giờ, thi công trong điều kiện biển êm dưới cấp 4 để giữ ổn định khoảng cách thu phát từ 4 đến 8 m. Đồng thời, hệ thống cần tích hợp anten DGPS với sai số vị trí dưới 3–4 m và khoảng cách giữa các phần tử Hydrophone chuẩn 1 m nhằm thỏa mãn điều kiện lấy mẫu không gian Nyquist.

Kết luận

  • Xây dựng thành công quy trình xử lý 5 bước tối ưu trên phần mềm Reflexw, nâng cao toàn diện chất lượng tài liệu địa chấn nông phân giải cao tại vùng biển Cửa Đáy.
  • Phục hồi hiệu quả năng lượng sóng phản xạ tầng sâu từ 60 ms trở xuống bằng hàm khuếch đại mũ tối ưu $g(t) = e^{0,1t}$.
  • Tăng cường tính liên tục của các ranh giới địa tầng thêm 30% và triệt tiêu nhiễu ngẫu nhiên nhờ thuật toán trung bình hóa 5 đường ghi liền kề.
  • Loại bỏ hoàn toàn phông nhiễu thủy động học (0–50 Hz) và nhiễu điện vi địa chấn (>1000 Hz) bằng bộ lọc dải thông cos-taper từ 100 Hz đến 1000 Hz.
  • Triệt tiêu dứt điểm sóng phản xạ nhiều lần đáy biển bậc hai ở vị trí 48 ms nhờ ứng dụng chính xác bộ lọc sai số tiên đoán trên miền thời gian - không gian.

Công trình là đóng góp học thuật và công nghệ quan trọng, tạo tiền đề để tái xử lý và nâng cao chất lượng cho hơn 40.000 km tuyến địa chấn thềm lục địa Việt Nam trong các giai đoạn tiếp theo. Quý độc giả, chuyên gia và đơn vị nghiên cứu có thể tham khảo toàn văn luận văn tại Trung tâm Thông tin - Thư viện Đại học Quốc gia Hà Nội để ứng dụng tổ hợp thuật toán này vào các dự án khảo sát địa chất công trình biển.