Tổng quan nghiên cứu

Bể trầm tích Phú Khánh nằm trên thềm lục địa miền Trung Việt Nam, sở hữu cấu trúc địa chất trải dài khoảng 250 km theo phương Bắc - Nam (từ 10°30’ đến 15°30’ vĩ độ Bắc) và rộng 50 - 75 km theo phương Đông - Tây (từ 109°20’ đến 112°30’ kinh độ Đông). Địa hình đáy biển có sự phân hóa sâu sắc từ vùng nước nông dưới 300 m ở phía Tây đến các trũng nước sâu từ 300 m đến hơn 4.000 m ở phía Đông. Mặc dù được đánh giá là bồn trũng nước sâu giàu tiềm năng hydrocarbon hàng đầu tại Đông Nam Á, lịch sử thăm dò tại đây còn thưa thớt và chưa đồng bộ. Đến nay, toàn bể mới chỉ ghi nhận 3 giếng khoan thăm dò tại vùng nước nông gồm TH-1X, CMT-1X (khoan năm 2009) và HT-1X (khoan năm 2011), hoàn toàn chưa có giếng khoan nào được thực hiện ở độ sâu vượt quá 1.000 m nước.

Các mô hình mô phỏng bồn trầm tích trước đây chủ yếu sử dụng giả định dòng nhiệt tĩnh cố định trong khoảng 54 - 65 mW/m² dựa trên nhiệt độ đo hiện tại ở các giếng nước nông, tạo ra độ bất định rất lớn khi dự báo thời điểm sinh và phân bố dầu khí ở vùng nước sâu. Luận văn tập trung giải quyết bài toán cốt lõi: Xây dựng mô hình lịch sử dòng nhiệt theo thời gian địa chất và mô phỏng hệ thống dầu khí 1D, 2D dọc tuyến địa chấn dài 200 km sâu đến 13.000 m. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp bức tranh nhiệt động lực học chuẩn xác nhằm định vị các bẫy chứa triển vọng, giảm thiểu rủi ro kinh tế và tối ưu hóa ngân sách hàng chục triệu USD cho các chiến dịch khoan thăm dò nước sâu trong tương lai.

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 kết hợp hai nền tảng nhiệt động lực học và kiến tạo bồn trũng kinh điển: Định luật dẫn nhiệt Fourier mở rộng kết hợp mô hình trung bình điều hòa của Beardsmore và Cull (2001), cùng Mô hình căng giãn đồng nhất nhiệt thạch quyển của McKenzie (1978) kết hợp nguyên lý cân bằng đẳng tĩnh Airy. Độ dẫn nhiệt của các lớp đất đá được chuẩn hóa theo từng kiểu thạch học: cát kết đạt 3,8 W/m/K, sét kết đạt 1,9 W/m/K, bột kết đạt 2,9 W/m/K và đá móng granit đạt 3,2 W/m/K.

Khung lý thuyết vận hành dựa trên 5 khái niệm chuyên ngành cốt lõi:

  • Dòng nhiệt cơ sở: Thông lượng nhiệt dẫn truyền từ manti qua thạch quyển vào đáy cột trầm tích.
  • Giải nén: Quá trình tính toán phục hồi bề dày và độ rỗng ban đầu của từng tập địa tầng trước khi bị nén ép cơ học do chôn vùi trầm tích.
  • Bóc tách tải trọng trầm tích: Phương pháp loại trừ ảnh hưởng của tải trọng trầm tích và cột nước biển để xác định chính xác độ sụt lún kiến tạo thuần túy do nội lực.
  • Hệ số tách giãn beta: Tỷ lệ thể hiện mức độ kéo mỏng của lớp vỏ lục địa trong các pha tách giãn kiến tạo.
  • Nhiệt độ ranh giới bề mặt trầm tích - nước biển và độ sâu mực nước biển cổ: Các điều kiện biên động điều chỉnh gradient nhiệt độ của bồn qua từng kỷ nguyên địa chất.

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

Cơ sở dữ liệu đầu vào bao gồm 14.500 km tuyến địa chấn 2D GeoStreamer, mặt cắt địa chấn 2D tái minh giải với độ phân giải lưới 1 km × 1 km kéo dài 200 km, báo cáo địa chất - địa vật lý của 3 giếng khoan thăm dò (TH-1X, CMT-1X, HT-1X), 42 điểm đo dòng nhiệt toàn cầu chuẩn hóa trong tọa độ từ 108° đến 115° kinh độ Đông và 10° đến 15° vĩ độ Bắc từ cơ sở dữ liệu của Goutorbe và cộng sự (2011), cùng thông số địa hóa (chỉ số TOC từ 4% đến 20%, chỉ số HI từ 300 đến 700 mg HC/g) từ các vết lộ Đầm Thị Nại và các bể lân cận.

Lý do lựa chọn phương pháp phân tích: Luận văn xây dựng thuật toán tự động hóa bằng mã nguồn VBA Macro trong Excel để giải phương trình giải nén tích phân phi tuyến và kỹ thuật bóc tách trầm tích cho 6 tập địa tầng của giếng CMT-1X. Hệ số tách giãn beta được xác định bằng phương pháp hồi quy tuyến tính giữa sụt lún kiến tạo và hàm suy giảm nhiệt thời gian với hằng số thời gian thạch quyển là 62,7 triệu năm. Quy trình này cho phép mô phỏng sự suy giảm nhiệt độ chính xác theo 3 pha tiến hóa kiến tạo của bể Phú Khánh: trước tách giãn, đồng tách giãn (56 đến 15,97 triệu năm trước) và hậu tách giãn (15,97 triệu năm trước đến nay). Mô hình bồn trầm tích 1D và 2D sau đó được mô phỏng tích hợp trên phần mềm chuyên dụng PetroMod kết hợp công cụ PetroRisk để phân tích độ nhạy bằng kỹ thuật lấy mẫu Monte Carlo.

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

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

Thứ nhất, kết quả phân tích 42 điểm dữ liệu bằng kỹ thuật Kriging cho thấy dòng nhiệt hiện tại ở bể Phú Khánh phân bố không đồng đều, dao động từ 52 đến 90 mW/m² với giá trị trung bình là 73 mW/m². Vùng nước nông phía Tây Bắc có giá trị dòng nhiệt cao từ 78 đến 86 mW/m², trong khi vùng nước sâu phía Đông Nam chỉ đạt từ 58 đến 72 mW/m².

Thứ hai, nghiên cứu định lượng được vai trò của hiệu ứng che phủ cách nhiệt từ tầng trầm tích. Tại vùng nước sâu, bề dày trầm tích tích tụ khổng lồ lên tới 10.000 m, cao gấp 3,3 đến 5 lần so với bề dày trầm tích 2.000 - 3.000 m tại các khối nâng giếng TH-1X và CMT-1X. Tầng trầm tích dày này đóng vai trò như một lớp vật liệu cách nhiệt, làm suy giảm từ 20% đến 30% giá trị dòng nhiệt thoát ra bề mặt đáy biển.

Thứ ba, hệ số tách giãn beta tính toán tại giếng CMT-1X đạt giá trị từ 1,35 đến 1,55, phản ánh sự mỏng hóa đáng kể của vỏ thạch quyển trong pha đồng tách giãn Eocene - Miocene Sớm. Lịch sử dòng nhiệt cổ ghi nhận giá trị tăng vọt cực đại ở cuối pha đồng tách giãn (đạt trên 85 mW/m²) trước khi giảm dần theo quy luật hàm mũ trong pha sụt lún nhiệt hậu tách giãn.

Thứ tư, mô hình hóa 1D và 2D làm sáng tỏ quá trình trưởng thành nhiệt của đá mẹ: Tầng đá sét đầm hồ tuổi Oligocene (kerogen loại I và II) tại các trũng sâu đã chính thức bước vào cửa sổ sinh dầu từ giai đoạn Miocene Giữa (khoảng 15 triệu năm trước) và hiện nay đang nằm trong cửa sổ sinh khí cực đại. Tầng đá mẹ châu thổ tuổi Miocene Dưới ở trung tâm bể hiện đang nằm trong cửa sổ sinh dầu thuận lợi.

Thảo luận kết quả

Sự phân hóa dòng nhiệt giữa vùng nước nông và nước sâu phản ánh rõ nét mối tương quan giữa độ dày lớp trầm tích và tốc độ nguội lạnh của thạch quyển. Do bể Phú Khánh là một bồn trũng trẻ (tuổi Kainozoi khoảng 35 triệu năm), thạch quyển vẫn chưa đạt trạng thái cân bằng nhiệt hoàn toàn (vốn cần khoảng 80 triệu năm theo mô hình McKenzie). Vì vậy, giá trị dòng nhiệt nền của toàn bể vẫn duy trì ở mức tương đối cao.

So với công trình trước đây của Nguyễn và cộng sự (2012) vốn áp dụng mức nhiệt cố định 54 - 65 mW/m², việc luận văn định lượng được đỉnh dòng nhiệt cổ giai đoạn đồng tách giãn đã chứng minh đá mẹ Oligocene đạt độ trưởng thành và sinh hydrocarbon sớm hơn từ 2 đến 3 triệu năm. Kết quả dòng nhiệt trung bình 73 mW/m² của nghiên cứu hoàn toàn tương thích và nhất quán với các nghiên cứu quy mô khu vực Biển Đông của He và cộng sự (trung bình 77 mW/m²) và nghiên cứu tại bể Cửa Sông Châu Giang của Li và cộng sự (trung bình 71 mW/m²).

Dữ liệu mô phỏng được trực quan hóa qua biểu đồ phân bố độ sâu giải nén 6 tập địa tầng, đồ thị hồi quy tuyến tính xác định hệ số tách giãn beta, cùng các mặt cắt 2D thể hiện đường đẳng nhiệt và hướng di cư chất lưu. Hydrocarbon sinh ra từ trũng sâu di cư theo phương thẳng đứng qua các đứt gãy kiến tạo sâu và thấm ngang qua các vỉa cát kết, carbonate ám tiêu san hô tuổi Miocene (độ rỗng từ 5% đến 20%) để nạp vào các bẫy cấu trúc vòm nâng và bẫy phi cấu tạo ở cả vùng nước nông lẫn nước sâu.

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

Thứ nhất, tiến hành khoan thăm dò thẩm lượng tại các cấu tạo bẫy vòm móng nứt nẻ và ám tiêu san hô Miocene ở vùng nước sâu từ 500 m đến 1.500 m. Chủ thể thực hiện là Tập đoàn Dầu khí Việt Nam cùng Tổng công ty Thăm dò Khai thác Dầu khí. Mục tiêu hướng tới nâng tỷ lệ thành công của các giếng khoan nước sâu thêm 15% đến 20% trong giai đoạn 2026 - 2028.

Thứ hai, triển khai đo đạc bổ sung nhiệt độ bề mặt đáy biển thời gian thực và khảo sát địa chấn phản xạ đáy giả lập BSR trên toàn bộ 14.500 km tuyến địa chấn vùng nước sâu. Chủ thể thực hiện là Viện Dầu khí Việt Nam phối hợp cùng các nhà thầu quốc tế. Mục tiêu nhằm giảm sai số tính toán dòng nhiệt hiện tại xuống dưới mức 5%, hoàn thành trước quý IV/2027.

Thứ ba, nâng cấp mô hình bồn trầm tích từ dạng 2D sang mô hình 3D động lực học độ phân giải cao cho khu vực trũng Phú Yên và đới cắt trượt Tuy Hòa. Chủ thể thực hiện là nhóm nghiên cứu chuyên sâu thuộc Trường Đại học Bách Khoa kết hợp các đơn vị vận hành mỏ. Mục tiêu là mô phỏng chi tiết 100% các đường di cư bẫy phi cấu tạo phức tạp, lộ trình thực hiện từ năm 2026 đến 2029.

Thứ tư, hoàn thiện khung pháp lý ưu đãi đầu tư và cơ chế phân chia rủi ro tài chính đối với các lô dầu khí nước sâu trên 1.000 m. Chủ thể thực hiện là Bộ Công Thương và các cơ quan quản lý năng lượng. Mục tiêu nhằm thu hút từ 2 đến 3 hợp đồng chia sản phẩm dầu khí mới với tổng vốn đầu tư nước ngoài ước đạt 300 - 500 triệu USD trước năm 2030.

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

Nhóm chuyên viên minh giải địa chất và kỹ sư mô hình hóa bồn trầm tích tại các công ty dầu khí: Ứng dụng quy trình tính toán dòng nhiệt cổ và bộ thông số thạch quyển thực nghiệm để thiết lập và hiệu chỉnh mô hình hệ thống dầu khí cho các lô tìm kiếm thăm dò thực tế.

Nhóm nghiên cứu sinh, giảng viên và học viên cao học chuyên ngành Kỹ thuật Địa chất, Địa vật lý và Kỹ thuật Dầu khí: Khai thác tài liệu như một cẩm nang học thuật mẫu mực về phương pháp giải nén địa tầng, bóc tách trầm tích cân bằng Airy và kỹ thuật lập trình VBA Macro xử lý dữ liệu địa chất tự động.

Nhóm chuyên gia đánh giá trữ lượng và nhà hoạch định chiến lược phát triển mỏ: Sử dụng bản đồ phân bố dòng nhiệt và các đới trưởng thành nhiệt của đá mẹ để khoanh định diện tích triển vọng, xếp hạng mức độ ưu tiên của các cấu trúc bẫy chứa nước sâu.

Nhóm nhà khoa học nghiên cứu địa nhiệt và kiến tạo mảng khu vực Biển Đông: Tham khảo tập dữ liệu 42 điểm đo dòng nhiệt chuẩn hóa và cơ chế tiến hóa manti để đánh giá tiềm năng khai thác năng lượng địa nhiệt sạch tại thềm lục địa Nam Trung Bộ.

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

Tại sao dòng nhiệt ở vùng nước sâu bể Phú Khánh lại thấp hơn vùng nước nông? Hiện tượng này xuất phát từ hiệu ứng che phủ cách nhiệt của lớp phủ trầm tích. Vùng nước sâu có bề dày trầm tích lên tới 10.000 m, hoạt động như một lớp chăn cách nhiệt làm suy giảm từ 20% đến 30% thông lượng nhiệt bề mặt so với trầm tích mỏng 2.000 - 3.000 m ở vùng nước nông.

Mô hình nhiệt McKenzie mang lại ưu thế gì vượt trội so với các cách tiếp cận trước đây? Mô hình McKenzie mô phỏng chính xác sự thay đổi nhiệt độ tức thời trong pha đồng tách giãn và quy luật suy giảm nhiệt tiệm tiến theo hàm mũ suốt 62,7 triệu năm của hằng số thời gian thạch quyển, khắc phục triệt để các sai số do giả định nhiệt độ cố định trước đây.

Đá mẹ nào đóng vai trò phát sinh dầu khí chủ lực tại bể Phú Khánh? Tầng đá sét kết và than đầm hồ tuổi Oligocene là nguồn sinh chính nhờ hàm lượng hữu cơ TOC cao từ 4% đến 20%, chỉ số HI đạt 300 đến 700 mg HC/g, đã hoàn tất pha sinh dầu và đang phát sinh khí mạnh mẽ tại các trũng sâu.

Ý nghĩa của phát hiện dầu tại giếng khoan nước nông CMT-1X đối với nghiên cứu này là gì? Phát hiện dầu trong đá vôi Miocene Dưới tại độ sâu 1.470,7 m của giếng CMT-1X là bằng chứng thực nghiệm khẳng định hệ thống dầu khí đang hoạt động, đồng thời cung cấp số liệu thạch học và nhiệt độ vỉa chuẩn xác để hiệu chỉnh toàn bộ mô hình nhiệt.

Thuật toán giải nén và bóc tách trầm tích trong luận văn được xử lý bằng công cụ nào? Quy trình được tự động hóa thông qua mã nguồn VBA Macro xây dựng trên nền tảng Excel nhằm giải hệ phương trình phi tuyến lặp cho 6 tập địa tầng, trước khi tích hợp vào phần mềm mô phỏng chuyên dụng PetroMod 1D/2D.

Kết luận

  • Xây dựng thành công quy trình khoa học chuẩn hóa để tính toán dòng nhiệt hiện tại và tái lập lịch sử dòng nhiệt cổ cho bồn trũng rift nước sâu.
  • Xác định giá trị dòng nhiệt trung bình toàn bể đạt 73 mW/m², đồng thời làm rõ sự phân dị nhiệt rõ nét giữa vùng nước nông (78 - 86 mW/m²) và vùng nước sâu (58 - 72 mW/m²).
  • Định lượng chính xác hiệu ứng cách nhiệt của cột trầm tích dày 10.000 m làm suy giảm tới 30% dòng nhiệt bề mặt tại vùng nước sâu.
  • Chứng minh đá mẹ Oligocene đã bước vào đỉnh sinh dầu từ 15 triệu năm trước, di cư nạp vào các vỉa chứa cát kết và carbonate san hô có độ rỗng 5% đến 20%.
  • Định hình lộ trình nâng cấp mô hình bồn trầm tích 3D và triển khai các giếng khoan thăm dò trọng điểm trong giai đoạn 2026 - 2030.

Công trình là bước tiến quan trọng trong nghiên cứu địa nhiệt dầu khí nước sâu tại Việt Nam. Các đơn vị thăm dò và nhà nghiên cứu quan tâm có thể ứng dụng trực tiếp quy trình mô hình hóa này để nâng cao hiệu quả đánh giá tiềm năng tài nguyên khoáng sản biển.