Tổng quan nghiên cứu

Ngành công nghiệp dầu khí giữ vai trò huyết mạch đối với an ninh năng lượng và sự phát triển kinh tế của Việt Nam. Trước thực trạng suy giảm sản lượng tại các mỏ truyền thống thuộc Miền võng Hà Nội, việc đẩy mạnh công tác tìm kiếm, thăm dò các đối tượng tiềm năng mới tại vùng biển nước nông ngoài khơi trở thành nhiệm vụ chiến lược cấp thiết. Bể trầm tích Sông Hồng có tổng diện tích khoảng 220.000 km2, trong đó phần thuộc chủ quyền thềm lục địa Việt Nam chiếm 126.000 km2, đặc trưng bởi chiều dày trầm tích Kainozoi có nơi đạt tới 17 km. Trong bức tranh cấu trúc rộng lớn này, khu vực lô 103-107 có diện tích ước tính khoảng 13.000 km2, nằm cách bờ biển Hải Phòng khoảng 150 km về phía Đông Nam với độ sâu mực nước biển dao động từ 30 m đến 45 m.

Vấn đề cốt lõi đặt ra là cấu trúc địa chất khu vực Đông Bắc lô 103 vô cùng phức tạp do chịu tác động mạnh mẽ của các pha nghịch đảo kiến tạo thời kỳ Miocen giữa đến Miocen muộn. Số lượng giếng khoan thăm dò tại đây còn thưa thớt và diện tích phủ sóng địa chấn 3D còn hạn chế, dẫn đến những bất định lớn về quy luật phân bố không gian của các thân cát chứa khí. Nghiên cứu tập trung giải quyết mục tiêu làm sáng tỏ đặc điểm cấu trúc kiến tạo, xây dựng mô hình địa chất 3D chi tiết cho các tầng chứa cát kết Miocen (ở khoảng độ sâu 2.200 m đến 3.000 m) tại mỏ khí X, đồng thời lượng hóa trữ lượng khí ban đầu tại chỗ. Việc hoàn thiện mô hình 3D có ý nghĩa thực tiễn to lớn, giúp giảm thiểu rủi ro địa chất, tiết kiệm hàng triệu USD chi phí khoan thẩm lượng và nâng cao hiệu quả thu hồi tài nguyên khí tự nhiên cho đất nước.

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 nền tảng của lý thuyết Địa chất thống kê (Geostatistics) do Matheron và Krige phát triển từ thập niên 1950, kết hợp với trường phái Địa tầng phân tập (Sequence Stratigraphy) hiện đại. Địa thống kê cung cấp công cụ toán học để mô tả các biến ngẫu nhiên không gian có tính liên tục và bất đồng nhất cao như độ rỗng, độ thấm và phân bố tướng đá. Các khái niệm nền tảng bao gồm:

  1. Giả thiết tính dừng (Stationarity): Quy luật phân bố xác suất của các thông số địa chất được xem là bất biến khi tịnh tiến trong một khoảng cách khảo sát nhất định, cho phép ngoại suy dữ liệu từ các vị trí giếng khoan đã biết sang các vùng chưa khoan.
  2. Tính bất đẳng hướng (Anisotropy): Phản ánh tốc độ biến đổi của thuộc tính vỉa chứa theo các phương không gian khác nhau, đặc biệt là sự khác biệt giữa phương song song và phương vuông góc với hướng dòng trầm tích cổ.
  3. Biểu đồ bán phương sai (Variogram): Công cụ định lượng mức độ tương quan không gian giữa các cặp điểm dữ liệu theo khoảng cách lag (h). Mô hình hàm cầu (Spherical model) được ưu tiên ứng dụng vì đảm bảo giá trị bán sai số tiệm cận đạt đúng 100% ngưỡng trần (Sill) tại bán kính ảnh hưởng (Range).
  4. Phương pháp nội suy Kriging và Mô phỏng xác suất Gauss (Gaussian Simulation): Tối ưu hóa trọng số ước tính không thiên lệch và tạo ra các kịch bản ngẫu nhiên có cùng xác suất xuất hiện, phục vụ đánh giá tính bất định của vỉa chứa.

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

Cơ sở dữ liệu của luận văn được tổng hợp từ nguồn tài liệu thực tế thu thập trong giai đoạn từ năm 2010 đến năm 2016 tại Tổng công ty Thăm dò Khai thác Dầu khí (PVEP). Cỡ mẫu nghiên cứu bao gồm toàn bộ tập dữ liệu carota địa vật lý liên tục với độ phân giải đo 0,1524 m tại giếng khoan phát hiện mỏ khí X, kết hợp đối sánh với hơn 10 giếng khoan thăm dò lân cận như 103-T-G-1X, 103-T-H-1X, 107-BAL-1X và 107-PA-1X, cùng hệ thống bản đồ địa chất tỷ lệ 1:1.000. Phương pháp chọn mẫu là chọn mẫu định tầng theo ranh giới địa tầng phân tập và mặt bất chỉnh hợp địa chấn.

Quy trình phân tích sử dụng nền tảng phần mềm chuyên dụng Petrel 2015 của Schlumberger. Lý do lựa chọn tổ hợp phương pháp này là khả năng tích hợp đa biến giữa tài liệu địa vật lý giếng khoan và thuộc tính địa chấn 3D (biên độ cực tiểu, tần số tức thời). Dữ liệu log giếng khoan có độ phân giải siêu cao được thô hóa (upscaling) vào mạng lưới ô mạng 3D có kích thước 50 m x 50 m theo phương ngang và 0,5 m theo phương thẳng đứng. Thuật toán mô phỏng tuần tự chỉ số (Sequential Indicator Simulation - SIS) được dùng để tái tạo phân bố thạch học, kết hợp thuật toán mô phỏng ngẫu nhiên Gauss (GRSF) để phân bố trường độ rỗng và độ thấm, bảo tồn tối đa tính bất đồng nhất tự nhiên của vỉa.

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

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

Quá trình mô hình hóa 3D và minh giải địa chất tại khu vực Đông Bắc lô 103 đã xác lập bốn phát hiện khoa học quan trọng:

  1. Khung cấu trúc mỏ khí X bị chi phối trực tiếp bởi 8 hệ thống đứt gãy kiến tạo chính. Cấu tạo mang hình thái vòm khép kín 3 chiều, trong đó cánh phía Đông bị khống chế bởi một đứt gãy nghịch vòng cung phương Tây Bắc - Đông Nam có biên độ trượt lớn. Một đứt gãy thuận phương Đông - Tây cắt ngang cấu tạo, chia mỏ thành hai khối: khối phía Bắc có đỉnh khép kín 4 chiều ở phần trên trầm tích Miocen giữa, và khối phía Nam khép kín 2 chiều tựa vào mặt đứt gãy.
  2. Tầng chứa chính thuộc thành tạo Miocen (các tập C50, C70 và C90) phân bố ở độ sâu từ 2.200 m đến 3.000 m có chất lượng đá chứa từ trung bình đến rất tốt. Độ rỗng hữu dụng biến thiên phổ biến từ 5% đến 20% (giá trị trung bình toàn mỏ đạt khoảng 13% đến 17%), độ thấm dao động mạnh từ 1 mD đến 100 mD (cục bộ các thể cát châu thổ đạt 120 mD). Tỉ số chiều dày cát chứa hiệu dụng trên tổng chiều dày vỉa (NTG) đạt mức cao từ 60% đến 75%.
  3. Hệ thống đá chắn khu vực và địa phương phát triển hoàn hảo. Các tập sét Miocen giữa và trên có bề dày tích lũy từ 10 m đến hơn 200 m, với hàm lượng khoáng vật sét mịn chiếm từ 65% đến 90%, tạo thành nắp chắn thủy lực kín khít bảo vệ các tích tụ khí.
  4. Kết quả tính toán thể tích trữ lượng khí tự do ban đầu tại chỗ (GIIP) thông qua 200 lần chạy mô phỏng Monte Carlo đã lượng hóa chính xác phổ xác suất trữ lượng ở 3 mức tin cậy: P90 (mức trữ lượng xác minh chắc chắn), P50 (mức trữ lượng có khả năng cao) và P10 (mức trữ lượng kỳ vọng tối đa).

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên sự phân hóa mạnh mẽ về chất lượng vỉa chứa mỏ X bắt nguồn từ môi trường trầm tích châu thổ ven biển kết hợp với các đợt biến đổi thứ sinh do nén ép kiến tạo. Pha trượt bằng phải của hệ đứt gãy Sông Hồng trong Miocen muộn đã tạo ra hiện tượng nâng cao cục bộ, bào mòn cắt xén một phần trầm tích nóc và tái phân bố lại các dòng chảy cát hạt mịn đến trung bình.

Khi so sánh với các tầng chứa Oligocen sâu hơn ở các lô phụ cận (nơi độ rỗng thường bị nén ép xuống dưới 10% và độ thấm dưới 5 mD), cát kết Miocen mỏ X vượt trội hơn từ 30% đến 50% về dung tích chứa và khả năng dẫn dòng. Biểu hiện thử dòng vỉa (DST) cho lưu lượng khí và condensate công nghiệp rất khả quan.

Về mặt hiển thị trực quan, các kết quả nghiên cứu được số hóa hoàn chỉnh trên Petrel: mô hình khung xương cấu trúc được minh họa qua các bản đồ đẳng sâu và đẳng dày tầng C50, C70, C90; mối quan hệ tương quan rỗng - thấm được thể hiện qua các biểu đồ phân tán hàm mũ; và sự phân dị độ bão hòa nước (Sw) theo chiều cao cột chất lưu được biểu diễn sắc nét trên các lát cắt thạch học 3D qua từng thân cát.

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

Dựa trên kết quả mô hình địa chất 3D và đánh giá bất định vỉa chứa, tác giả đưa ra bốn giải pháp mang tính chiến lược:

  1. Tối ưu hóa vị trí giếng khoan thẩm lượng và khai thác: Tổng công ty Thăm dò Khai thác Dầu khí (PVEP) cần ưu tiên thiết kế 2 đến 3 giếng khoan phát triển tập trung vào đỉnh vòm của khối cấu trúc phía Bắc thuộc tầng C70 và C90. Mục tiêu là đón đầu các thể cát có tỷ lệ NTG lớn hơn 70%, nâng cao xác suất thành công thương mại lên trên 85% trong giai đoạn 2017 - 2020.
  2. Tích hợp chuyển giao sang mô hình mô phỏng động lực học khai thác: Bộ phận công nghệ mỏ cần tiếp nhận ngay dữ liệu lưới 3D tĩnh để tiến hành mô phỏng dòng chảy đa pha (Reservoir Simulation). Việc này giúp khớp lịch sử khai thác thử vỉa, dự báo chính xác độ suy giảm áp suất đáy giếng và xác lập sơ đồ bố trí công trình khai thác tối ưu trong vòng đời 10 năm của mỏ.
  3. Mở rộng diện tích thu nổ địa chấn 3D phân giải cao: Đề xuất Tập đoàn Dầu khí Việt Nam triển khai khảo sát địa chấn 3D bổ sung với diện tích khoảng 500 km2 bao trùm phần kéo dài sang các lô lân cận (102, 106, 107). Mục tiêu là nhận diện sớm các bẫy địa tầng dạng thấu kính cát mỏng dưới 10 m vốn chưa thể phân giải trên lưới địa chấn 2D cũ.
  4. Chuẩn hóa quy trình cập nhật mô hình theo thời gian thực: Ban điều hành dự án cần áp dụng quy trình cập nhật địa thống kê thích ứng liên tục (Ensemble Kalman Filter) khi có thêm dữ liệu đo trong khi khoan (LWD/MWD) từ các giếng mới, kiểm soát độ lệch chuẩn tính toán trữ lượng xuống dưới 12% trước khi phê duyệt Kế hoạch phát triển mỏ (FDP).

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn và học thuật cho bốn nhóm đối tượng chính:

  1. Kỹ sư địa chất và địa vật lý dầu khí: Nắm vững quy trình chuẩn từ khâu xử lý tài liệu carota, phân tích thuộc tính địa chấn biên độ cực tiểu đến kỹ thuật thô hóa dữ liệu (upscaling) và mô phỏng thạch học ngẫu nhiên (SIS) trên nền tảng Petrel 2015.
  2. Nhà quản lý dự án và chuyên viên kinh tế mỏ tại PVN, PVEP: Sử dụng bộ thông số trữ lượng P90, P50, P10 và các phân tích rủi ro cấu trúc làm cơ sở vững chắc để xây dựng báo cáo khả thi đầu tư, tối ưu hóa ngân sách khoan thăm dò trị giá hàng chục triệu USD.
  3. Giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Địa chất Dầu khí: Khai thác làm tài liệu tham khảo chuyên sâu về cơ chế nghịch đảo kiến tạo Miocen bồn trũng Sông Hồng và phương pháp ứng dụng hàm bán phương sai Variogram trong mô hình hóa vỉa chứa.
  4. Kỹ sư công nghệ mỏ và khai thác công trình biển: Tham khảo các đặc trưng rỗng - thấm, đường cong áp suất mao dẫn và mô hình phân bố độ bão hòa nước Sw để thiết kế hệ thống giàn khai thác, đường ống dẫn khí và tính toán hệ số thu hồi khí tối đa.

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

  1. Tại sao thành tạo Miocen tại Đông Bắc lô 103 lại được xem là đối tượng tìm kiếm thăm dò khí trọng điểm? Thành tạo Miocen (đặc biệt là Miocen giữa thuộc Hệ tầng Phù Cừ) sở hữu chất lượng đá chứa vượt trội với độ rỗng trung bình từ 13% đến 20% và độ thấm đạt tới 100 - 120 mD. Vỉa chứa nằm trong đới nghịch đảo kiến tạo tạo nên các bẫy cấu trúc vòm rộng lớn, được nạp bởi đá sinh giàu vật chất hữu cơ có chỉ số HI từ 300 đến 694 mg HC/g TOC.

  2. Việc lựa chọn kích thước ô lưới 50 m x 50 m x 0,5 m trong mô hình cấu trúc 3D mang lại lợi thế gì? Kích thước ô lưới này tạo ra sự cân bằng hoàn hảo giữa độ mịn không gian và năng lực xử lý của hệ thống máy tính. Chiều thẳng đứng 0,5 m bảo tồn trọn vẹn tính bất đồng nhất của các lớp kẹp cát sét, trong khi kích thước mặt bằng 50 m x 50 m phản ánh chính xác ranh giới 8 đứt gãy chính mà không làm phát sinh lỗi vặn xoắn ô mạng.

  3. Biểu đồ Variogram đóng vai trò cốt lõi như thế nào trong xây dựng mô hình địa chất? Variogram định lượng hóa mối tương quan không gian giữa các điểm dữ liệu giếng khoan theo khoảng cách và phương vị. Thông qua việc khớp hàm mẫu (mô hình cầu), Variogram xác định bán kính ảnh hưởng và ngưỡng trần phương sai, cung cấp trọng số toán học chính xác cho thuật toán nội suy Kriging và mô phỏng ngẫu nhiên SIS.

  4. Yếu tố nào đóng vai trò là tầng chắn giữ khí tại mỏ X? Đóng vai trò chắn là các tập sét xen tầng dày từ 10 m đến 150 m tuổi Oligocen - Miocen dưới (hàm lượng sét 65% - 90%) làm tầng chắn địa phương, cùng các phức hệ sét Miocen giữa đến Pliocen dày trên 200 m phân bố rộng khắp, tạo nên tầng chắn khu vực hoàn toàn kín khít.

  5. Vì sao cần chạy mô phỏng ngẫu nhiên 200 kịch bản để tính toán trữ lượng khí tại chỗ? Do các thông số thể tích vỉa, độ rỗng và độ bão hòa nước luôn chứa đựng sai số bất định, việc mô phỏng Monte Carlo 200 lần giúp xây dựng đường cong phân bố xác suất lũy tích. Kết quả này cung cấp các kịch bản chuẩn P90, P50 và P10, giúp nhà điều hành định lượng chính xác mức độ rủi ro kinh tế trước khi ra quyết định đầu tư.

Kết luận

  • Luận văn đã làm sáng tỏ cấu trúc kiến tạo phức tạp của mỏ khí X tại Đông Bắc lô 103, chỉ ra vai trò quyết định của 8 đứt gãy chính và pha nghịch đảo Miocen muộn trong việc định hình bẫy chứa dầu khí.
  • Xây dựng thành công mô hình địa chất 3D có độ tin cậy cao trên phần mềm Petrel với kích thước ô lưới 50 m x 50 m x 0,5 m, tích hợp chặt chẽ tài liệu địa vật lý giếng khoan và thuộc tính địa chấn.
  • Lượng hóa các đặc trưng vỉa chứa chính: cát kết Miocen có độ rỗng trung bình 13% - 17%, độ thấm 1 - 100 mD, tỷ lệ cát chứa hiệu dụng NTG đạt 60% - 75%, được che chắn an toàn bởi các tầng sét dày hơn 200 m.
  • Xác lập phổ phân bố trữ lượng khí tại chỗ GIIP theo phương pháp thể tích và phân tích độ nhạy xác suất P90, P50, P10, tạo cơ sở khoa học vững chắc cho việc lập kế hoạch phát triển mỏ.
  • Khuyến nghị PVEP và PVN sớm đưa mô hình tĩnh 3D vào mô phỏng khai thác động lực học và chuẩn bị vị trí khoan 2 đến 3 giếng thẩm lượng trong giai đoạn tiếp theo nhằm hiện thực hóa dòng khí thương mại. Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác mô hình này làm tài liệu chuẩn mực cho các cấu tạo tương tự trong toàn bể trầm tích Sông Hồng.