Tổng quan nghiên cứu

Thăm dò địa chấn phản xạ là phương pháp then chốt trong địa vật lý ứng dụng, chiếm tới hơn 80% tổng chi phí khảo sát cấu trúc sâu phục vụ tìm kiếm khoáng sản và dầu khí. Nhiệm vụ cốt lõi của phương pháp là xác định chính xác độ sâu, hình dạng các mặt ranh giới địa tầng và vận tốc truyền sóng đàn hồi trong các lớp đất đá. Mặc dù địa chấn lỗ khoan cung cấp số liệu vận tốc trực tiếp đáng tin cậy, nhưng chi phí thi công giếng khoan sâu rất lớn và không phải khu vực nào cũng có sẵn giếng khoan. Do đó, việc giải bài toán ngược để tính toán vận tốc truyền sóng trực tiếp từ các băng địa chấn và biểu đồ thời khoảng trở thành yêu cầu cấp thiết.

Luận văn tập trung nghiên cứu cơ sở lý thuyết truyền sóng đàn hồi, đánh giá các giải pháp xác định tốc độ truyền sóng hiện hành và đi sâu thử nghiệm phương pháp các điểm tương hỗ kết hợp thuật toán chia đôi. Phạm vi khảo sát được thực hiện trên 5 mô hình địa chất số hai lớp có góc nghiêng ranh giới biến thiên từ âm 5 độ đến dương 5 độ, độ sâu đáy tầng khảo sát dao động từ 30 mét đến 86 mét, với dải tần số xử lý dao động dưới 125 Hz. Ý nghĩa của nghiên cứu thể hiện qua việc giảm thiểu sai số tính toán vận tốc xuống mức dưới 0,5%, rút ngắn thời gian lặp hội tụ nghiệm chỉ từ 12 đến 16 bước, qua đó nâng cao độ chính xác khi minh giải cấu trúc bẫy dầu khí và địa chất công trình.

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 dựa trên cơ sở lý thuyết đàn hồi môi trường liên tục và quang hình học địa chấn. Theo định luật Hooke, môi trường đất đá chịu biến dạng đàn hồi được đặc trưng bởi các tham số như mô đun Young, hệ số Poisson, mật độ đất đá và các hằng số Lame. Quá trình lan truyền sóng dọc (sóng P) và sóng ngang (sóng S) trong môi trường đồng nhất được mô tả qua các biểu thức vận tốc độc lập, trong đó tỷ số vận tốc sóng dọc trên sóng ngang phụ thuộc chặt chẽ vào hệ số Poisson.

Bên cạnh đó, nghiên cứu vận dụng các định luật quang hình học bao gồm định luật Snellius về sự phản xạ - khúc xạ tia sóng, nguyên lý Huygens - Fresnel về nguồn sóng thứ cấp và nguyên lý Fermat về thời gian truyền sóng cực tiểu. Vận tốc biểu kiến được thiết lập dựa trên đạo hàm của biểu đồ thời khoảng theo khoảng cách quan sát. Nghiên cứu cũng hệ thống hóa các khái niệm then chốt như vận tốc hiệu dụng, mặt đẳng thời, trục đồng pha và mối liên hệ hình học của tia sóng theo định luật Bendorff để làm tiền đề cho việc giải bài toán ngược tìm vận tốc lớp.

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

Nguồn dữ liệu của luận văn được tạo lập từ việc giải bài toán thuận trên 5 mô hình địa chất số hai lớp có cấu trúc phẳng và nghiêng. Cỡ mẫu nghiên cứu bao gồm một hệ thống quan sát gồm 13 máy thu tín hiệu đặt cách đều nhau 5 mét dọc tuyến đo, kết hợp 2 điểm nổ mìn tương hỗ đặt tại vị trí 200 mét (điểm nổ trái) và 250 mét (điểm nổ phải) với khoảng nổ cách nhau 50 mét.

Phương pháp chọn mẫu sử dụng kỹ thuật lấy mẫu lưới đều trên tuyến quan sát, đảm bảo ghi nhận liên tục thời gian truyền sóng phản xạ từ các mặt ranh giới. Luận văn lựa chọn phương pháp phân tích các điểm tương hỗ kết hợp thuật toán chia đôi số học để tìm kiếm nghiệm vận tốc. Lý do lựa chọn giải thuật này là vì nó loại bỏ được các giả định đơn giản hóa về môi trường đồng nhất của phương pháp bình phương tọa độ, đồng thời tiết kiệm đáng kể tài nguyên máy tính so với phương pháp quét phổ vận tốc điểm sâu chung vốn đòi hỏi khối lượng xử lý dữ liệu khổng lồ. Kỹ thuật xấp xỉ đạo hàm parabol bậc hai được áp dụng trên các trục đồng pha giúp xác định chính xác góc ló tia sóng tại từng điểm thu.

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

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

Thử nghiệm trên 5 mô hình địa chất số đã mang lại các kết quả định lượng chuẩn xác về khả năng hội tụ và độ chính xác của phương pháp các điểm tương hỗ:

  • Tại Mô hình 1 (cấu trúc hai lớp nằm ngang song song, vận tốc thực lớp 1 là 1000 m/s ở độ sâu 30 mét, lớp 2 là 1500 m/s ở độ sâu 60 mét), thuật toán hội tụ vận tốc lớp 1 đạt 1000,70 m/s sau 16 lần lặp với sai số tương đối chỉ 0,07%. Đối với lớp 2, vận tốc xác định đạt 1500,56 m/s chỉ sau 12 chu kỳ lặp với sai số tương đối cực nhỏ là 0,037%, độ chênh lệch thời gian hội tụ đạt mức dưới 0,000001 giây.
  • Tại Mô hình 2 (lớp 1 nằm ngang, lớp 2 nghiêng 3 độ), vận tốc lớp 1 giữ nguyên độ chính xác 1000,70 m/s, trong khi vận tốc lớp 2 đạt 1499,00 m/s sau 15 lần lặp. Sai số tuyệt đối chỉ là 1,00 m/s, tương ứng với tỷ lệ sai số 0,067%.
  • Tại Mô hình 3 và Mô hình 4 (các cấu trúc bất đối xứng có cả hai lớp cùng nghiêng từ 3 độ đến 5 độ theo chiều thuận và chiều nghịch, độ sâu đáy tầng từ 77 mét đến 86 mét), sai số tuyệt đối của vận tốc lớp thứ hai dao động trong khoảng từ 6,47 m/s đến 6,58 m/s, tương đương sai số tương đối chỉ khoảng 0,43% đến 0,44%.
  • Tốc độ hội tụ của thuật toán chia đôi cho thấy hiệu quả vượt trội, giảm số bước lặp từ hàng trăm bước của phương pháp dò tìm từng bước xuống chỉ còn 12 đến 16 lần lặp trên tất cả các kịch bản mô phỏng.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự gia tăng nhẹ của sai số từ 0,037% ở mô hình phẳng lên 0,44% ở các mô hình nghiêng 5 độ là do sự dịch chuyển hình học của điểm phản xạ thực tế dưới sâu khi ranh giới bị uốn nghiêng. Khi ranh giới có độ nghiêng, góc tới và góc phản xạ của tia sóng thay đổi bất đối xứng, gây ra độ lệch nhỏ trong phép tính đạo hàm bậc hai của biểu đồ thời khoảng tại vị trí các trạm thu. Tuy nhiên, mức sai số dưới 0,5% vẫn vượt trội hơn hẳn so với sai số thông thường từ 3% đến 5% của các phương pháp xác định vận tốc hiệu dụng truyền thống như phương pháp hiệu số cố định hay phương pháp biểu đồ giao nhau.

Các dữ liệu thực nghiệm này có thể được trình bày một cách trực quan thông qua bảng đối sánh sai số vận tốc giữa 5 mô hình địa tầng, kết hợp biểu đồ đường biểu diễn sự suy giảm của độ chênh lệch thời gian theo từng chu kỳ lặp từ 1 đến 16. Đồng thời, đồ thị biểu diễn họ tia sóng đi từ hai điểm nổ tương hỗ cắt nhau tại mặt ranh giới sẽ minh họa rõ ràng độ tin cậy của thuật toán khi xử lý các bề mặt phản xạ có góc dốc phức tạp.

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

Nhằm hoàn thiện và mở rộng khả năng ứng dụng của phương pháp xác định tốc độ truyền sóng trong công tác thăm dò thực tế, luận văn đưa ra 4 khuyến nghị trọng tâm:

  • Tích hợp mô-đun giải thuật chia đôi và phương pháp điểm tương hỗ vào các phần mềm xử lý số liệu địa chấn thương mại trong vòng 6 tháng tới, hướng tới mục tiêu tự động hóa việc tính toán vận tốc lớp cho các trạm ghi số có bước lấy mẫu 2 ms và dải động học từ 100 dB đến 120 dB. Chủ thể thực hiện là các chuyên viên lập trình thuật toán địa vật lý.
  • Chuẩn hóa mạng lưới quan sát thực địa với khoảng cách giữa các đầu thu dao động từ 5 mét đến 10 mét, đảm bảo tối thiểu 12 đến 24 kênh thu sóng rõ nét cho mỗi chặng nổ nhằm tăng tỷ số tín hiệu trên nhiễu lên trên 25% trong quý 3 năm 2026. Chủ thể thực hiện là các đội trưởng đội khảo sát địa chấn thực địa.
  • Triển khai thử nghiệm kết hợp phương pháp điểm tương hỗ với tài liệu carota âm thanh tại ít nhất 3 giếng khoan thăm dò trong vòng 12 tháng, nhằm hiệu chỉnh tham số độ nghiêng ranh giới từ 3 độ đến 15 độ trong môi trường địa chất phức tạp. Chủ thể thực hiện là các kỹ sư minh giải địa vật lý giếng khoan.
  • Mở rộng phạm vi nghiên cứu của mô hình thuật toán từ môi trường 2 lớp sang môi trường đa lớp liên tục có từ 3 ranh giới trở lên và môi trường dị hướng vận tốc trong giai đoạn 2026 - 2027. Chủ thể thực hiện là các nhóm nghiên cứu tại các viện địa chất và trường đại học chuyên ngành.

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

Công trình nghiên cứu mang lại giá trị học thuật và thực tiễn cao cho nhiều nhóm đối tượng chuyên môn:

  • Kỹ sư địa vật lý thăm dò dầu khí: Có thể khai thác các công thức giải bài toán ngược và thuật toán điểm tương hỗ để xác định chính xác mô hình vận tốc lớp, phục vụ chuyển đổi từ mặt cắt địa chấn thời gian sang mặt cắt địa chấn độ sâu, giúp định vị chính xác các bẫy chứa dầu khí ở độ sâu trên 1000 mét.
  • Chuyên gia khảo sát địa chất công trình và khoáng sản: Ứng dụng phương pháp để phân tích các lát cắt địa chấn nông từ 10 mét đến 100 mét, phục vụ đánh giá tính chất cơ lý nền móng công trình xây dựng, đập thủy điện và thăm dò sa khoáng ven biển.
  • Giảng viên và học viên cao học ngành Vật lý địa cầu: Sử dụng luận văn làm tài liệu học tập chuẩn mực về lý thuyết trường sóng đàn hồi, kỹ thuật phân tích biểu đồ thời khoảng và phương pháp mô phỏng số bài toán địa chấn phản xạ.
  • Các nhà phát triển phần mềm khoa học địa chất: Tham khảo sơ đồ khối và thuật toán chia đôi số học lặp 12 đến 16 bước để lập trình các công cụ tự động trích xuất vận tốc và vẽ mặt cắt ranh giới địa tầng.

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

Vận tốc hiệu dụng khác biệt như thế nào so với vận tốc thực của lớp đất đá?

Vận tốc hiệu dụng là đại lượng vận tốc trung bình giả định tính từ bề mặt quan sát đến mặt ranh giới phản xạ dưới sâu, xem môi trường phía trên như một khối đồng nhất. Trong khi đó, vận tốc thực phản ánh chính xác đặc tính cơ lý đàn hồi của riêng từng lớp đất đá cụ thể, giúp xác định đúng bản chất thạch học và ranh giới cấu tạo địa chất.

Tại sao phương pháp các điểm tương hỗ bắt buộc phải sử dụng 2 điểm phát sóng đối xứng?

Phương pháp đòi hỏi 2 điểm phát sóng đặt ở hai đầu tuyến đo cách nhau khoảng 50 mét để tạo ra cặp biểu đồ thời khoảng ngược chiều. Việc ghi nhận sóng phản xạ tương hỗ giúp xác định chính xác thời gian truyền sóng giữa hai điểm cố định, loại bỏ ảnh hưởng bất đồng nhất của đới tốc độ nhỏ gần mặt đất và khử sai số do địa hình gây ra.

Thuật toán chia đôi mang lại ưu thế gì vượt trội so với phương pháp dò tìm từng bước?

Thuật toán chia đôi liên tục thu hẹp khoảng nghiệm chứa vận tốc theo từng nửa khoảng cách, giúp tốc độ hội tụ nghiệm nhanh gấp nhiều lần. Trong các mô hình thử nghiệm, thuật toán chỉ cần từ 12 đến 16 bước lặp để đạt độ chính xác chênh lệch thời gian dưới 0,000001 giây, trong khi phương pháp dò từng bước đòi hỏi hàng trăm phép tính lặp.

Độ nghiêng của mặt ranh giới ảnh hưởng như thế nào đến độ chính xác của vận tốc?

Khi góc nghiêng ranh giới tăng từ 0 độ lên 5 độ, sai số tương đối của phép tính vận tốc lớp thứ hai tăng từ 0,037% lên mức 0,44%. Sự gia tăng này xuất phát từ việc điểm phản xạ thực tế bị trượt dọc theo mặt phẳng nghiêng, làm thay đổi góc ló tia sóng và giá trị đạo hàm của biểu đồ thời khoảng tại điểm thu.

Phương pháp này có áp dụng được cho công tác khảo sát địa chấn nông nền móng công trình không?

Phương pháp hoàn toàn khả thi cho các khảo sát địa chấn nông ở độ sâu từ 10 mét đến 100 mét bằng cách sử dụng nguồn phát không nổ như búa đập trọng lượng 2 đến 3 tấn hoặc nguồn rung. Việc bố trí khoảng cách máy thu dày đặc từ 1 mét đến 5 mét sẽ cung cấp độ phân giải cao để đánh giá tầng phong hóa.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết đàn hồi, quang hình học địa chấn và các giải pháp xác định vận tốc truyền sóng trong địa chấn phản xạ.
  • Phương pháp các điểm tương hỗ kết hợp thuật toán chia đôi được chứng minh là công cụ số học ưu việt, cho phép tính toán trực tiếp vận tốc lớp và tọa độ ranh giới phản xạ.
  • Thử nghiệm trên 5 mô hình địa chất số khẳng định độ chính xác rất cao với sai số vận tốc luôn được khống chế dưới 0,45% và số lần lặp hội tụ tối ưu từ 12 đến 16 chu kỳ.
  • Trong giai đoạn 2026 - 2027, định hướng phát triển trọng tâm là nâng cấp thuật toán cho môi trường đa tầng bất đồng nhất và ứng dụng kiểm chứng trên số liệu địa chấn thực tế ngoài thực địa.
  • Các đơn vị thăm dò địa vật lý và viện nghiên cứu chuyên ngành nên sớm tiếp cận, chuẩn hóa quy trình và tích hợp thuật toán này vào dây chuyền xử lý dữ liệu địa chấn hiện đại.