Tổng quan nghiên cứu

Bồn trũng Cửu Long là bể trầm tích có tiềm năng dầu khí quan trọng bậc nhất thềm lục địa Việt Nam với diện tích khoảng 36.000 km² và trữ lượng dầu khí tại chỗ ban đầu ước tính từ 3,2 đến 3,4 tỷ tấn dầu quy đổi. Mỏ Gấu Đen nằm ở phía Đông Bắc Lô 16-1 thuộc bồn trũng Cửu Long, cách bờ biển Vũng Tàu khoảng 100 km về phía Đông Nam và cách mỏ Bạch Hổ khoảng 20 km. Dòng dầu thương mại đầu tiên tại mỏ được phát hiện vào năm 2005 qua giếng khoan GD-1X. Tuy nhiên, Gấu Đen thuộc dạng mỏ nhỏ cận biên với cấu trúc địa chất phân khối phức tạp, bị chia cắt mạnh bởi hệ thống đứt gãy dạng bậc thang và tầng chứa có tính bất đồng nhất cao theo cả phương thẳng đứng lẫn phương nằm ngang.

Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình địa chất 3 chiều có độ chính xác cao cho tầng chứa cát kết lục nguyên nội phụ Hệ tầng Bạch Hổ dưới (ILBH 5.2) tuổi Mioxen sớm tại khu vực phía Bắc mỏ Gấu Đen. Mô hình này đóng vai trò nền tảng để thô hóa phục vụ công tác mô phỏng khai thác động, thiết kế quỹ giếng và giảm thiểu rủi ro kinh tế khi phát triển mỏ. Phạm vi nghiên cứu tập trung vào các tích tụ dầu khối H1 và H2, khai thác nguồn dữ liệu tái xử lý địa chấn miền độ sâu PSDM năm 2010 với diện tích 589 km² cùng chuỗi tài liệu từ các giếng khoan thăm dò và thẩm lượng. Kết quả nghiên cứu đã lượng hóa thành công trữ lượng dầu tại chỗ đạt 105 triệu thùng, làm cơ sở khoa học vững chắc cho việc lập kế hoạch phát triển mỏ với hệ thống các giếng khai thác từ GD-1P đến GD-8P.

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 kiến tạo bể rift nội lục điển hình của bồn trũng Cửu Long, trải qua các pha biến dạng kéo tách và nén ép nghịch đảo từ Eoxen đến Mioxen. Cơ sở địa chất mỏ dựa trên lý thuyết đánh giá vỉa chứa cát kết trầm tích sông - hồ và đầm lầy ven bờ.

Khung mô hình hóa ứng dụng lý thuyết đơn vị dòng chảy thủy lực (Hydraulic Flow Units - HU) nhằm khắc phục hạn chế của việc phân loại tướng đá nhị phân truyền thống. Lý thuyết HU tích hợp mối quan hệ giữa chỉ số chất lượng vỉa chứa (Reservoir Quality Index - RQI) và chỉ số đới dòng chảy (Flow Zone Indicator - FZI) rút ra từ phương trình Kozeny-Carman.

Các khái niệm chính được ứng dụng trong mô hình gồm:

  • Độ rỗng hiệu dụng: Tỷ lệ phần trăm thể tích lỗ rỗng liên thông chứa chất lưu di động.
  • Độ thấm tuyệt đối: Khả năng dẫn chất lưu qua môi trường rỗng dưới gradient áp suất, đo bằng đơn vị miliDarcy (mD).
  • Độ bão hòa nước liên kết: Tỷ lệ thể tích nước không thể di động bị giữ lại trong các mao quản nhỏ của đá chứa.
  • Hàm biến sai (Variogram): Đại lượng thống kê không gian xác định phạm vi ảnh hưởng và tính dị hướng của các thuộc tính địa chất theo các phương chính.
  • Biến đổi chuẩn hóa (Normal Score Transform): Kỹ thuật chuyển đổi phân bố dữ liệu thực tế về phân bố chuẩn Gaussian chuẩn tắc phục vụ mô phỏng ngẫu nhiên.

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

Nguồn dữ liệu đầu vào bao gồm 589 km² tài liệu địa chấn 3D tái xử lý miền độ sâu PSDM năm 2010, tập hợp đường cong địa vật lý giếng khoan (gamma ray, mật độ, neutron, điện trở suất, siêu âm) từ 5 giếng khoan thăm dò và thẩm lượng trong khu vực mỏ, kết quả đo áp suất vỉa liên tục (MDT, RCI) và tài liệu phân tích mẫu lõi thường cùng mẫu lõi đặc biệt (SCAL).

Cỡ mẫu nghiên cứu bao gồm hàng trăm mẫu lõi được đo đạc chính xác về độ rỗng, độ thấm và áp suất mao dẫn. Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu tầng theo ranh giới thạch học của các tập vỉa chứa cát kết trong khoảng độ sâu từ 2.600 m đến 3.200 m. Lý do lựa chọn cỡ mẫu này là nhằm đảm bảo tính đại diện cho toàn bộ các tập vỉa ILBH 5.2 có chiều dày biến đổi từ 15 m đến 60 m.

Phương pháp phân tích địa thống kê được lựa chọn vì tính ưu việt trong việc xử lý các thuộc tính bất đồng nhất cao:

  • Thuật toán mô phỏng chỉ định tuần tự (Sequential Indicator Simulation - SIS) được áp dụng để mô hình hóa phân bố không gian rời rạc của các đơn vị dòng chảy HU, bảo tồn hình thái phân bố phi tuyến của các thể cát lòng sông và cát tràn bờ.
  • Thuật toán mô phỏng tuần tự ngẫu nhiên theo Gaussian (Sequential Gaussian Simulation - SGS) được sử dụng để xây dựng mô hình thuộc tính liên tục cho độ rỗng và độ thấm. Phương pháp này đảm bảo tái lập chính xác dải biến thiên thực tế của dữ liệu đo tại giếng khoan mà không làm phẳng cục bộ các giá trị cực trị.
  • Mốc thời gian thực hiện nghiên cứu kéo dài từ giai đoạn thẩm lượng năm 2008 đến khi hoàn thiện mô hình chi tiết và thử nghiệm lịch sử khai thác vào giai đoạn 2011 - 2012.

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

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

Thứ nhất, nghiên cứu đã phân chia thành công tầng chứa cát kết ILBH 5.2 thành 4 đơn vị dòng chảy thủy lực (HU) riêng biệt. Kết quả phân tích từ mẫu lõi cho thấy độ thấm biến thiên mạnh từ 1 mD ở các tướng cát bột mịn đầm hồ (HU 1) đến hơn 250 mD ở các tướng cát lòng sông cổ năng lượng cao (HU 4), trong khi độ rỗng hiệu dụng dao động từ 10% đến 24%.

Thứ hai, mô hình cấu trúc 3 chiều xây dựng trên mạng lưới tọa độ không gian mịn đã làm sáng tỏ hệ thống 12 đứt gãy chính dạng bậc thang phân chia khu vực phía Bắc thành 2 khối tích tụ H1 và H2. Ranh giới dầu nước được xác định chính xác từ dữ liệu đo áp suất MDT/RCI kết hợp tài liệu địa vật lý giếng khoan, giúp xác định thể tích khối đá chứa dầu hiệu dụng với độ tin cậy vượt trội.

Thứ ba, tính toán trữ lượng dầu tại chỗ (STOIIP) từ mô hình địa chất 3 chiều đạt giá trị 105 triệu thùng dầu. Kết quả này có độ tương thích cao khi so sánh với phương pháp thể tích truyền thống nhưng chi tiết hơn về mặt phân bố không gian trữ lượng theo từng khối kiến tạo và từng lớp vỉa nhỏ.

Thứ tư, quá trình thô hóa (Upscaling) đã chuyển đổi mạng lưới địa chất hàng triệu ô lưới về mô hình mô phỏng khai thác với số lượng ô lưới giảm hơn 70%, đồng thời bảo toàn trên 95% thể tích rỗng và hướng dị hướng của độ thấm. Kết quả lặp lại lịch sử thử vỉa (History Matching) trên các giếng thử vỉa đạt độ sai lệch áp suất đáy giếng dưới 3%, minh chứng cho độ tin cậy của mô hình tĩnh.

Thảo luận kết quả

Sự biến thiên phức tạp của các thuộc tính vật lý thạch học tại mỏ Gấu Đen bắt nguồn từ môi trường trầm tích chuyển tiếp giữa đầm hồ nước ngọt và đồng bằng bồi tích ven bờ. Việc áp dụng phương pháp phân loại tướng đá theo đơn vị dòng chảy HU giúp khắc phục triệt để nhược điểm của mô hình nhị phân cát/sét thông thường. Mô hình truyền thống thường đánh giá quá cao tính liên thông của vỉa chứa, dẫn đến dự báo sai thời điểm nước đột phá vào giếng khai thác.

Dữ liệu mô hình hóa được trực quan hóa qua các mặt cắt địa chất 3 chiều hướng Bắc - Nam, bản đồ đẳng sâu nóc vỉa ILBH 5.2 và biểu đồ phân bố tương quan giữa độ rỗng và độ thấm. Biểu đồ áp suất mao dẫn theo chiều cao thân dầu đã mô tả chính xác sự chuyển tiếp của độ bão hòa nước từ vùng nước đáy (Sw = 100%) lên đỉnh vòm tích tụ (Sw giảm xuống dưới 25%). Kết quả này giải thích rõ nguyên nhân lưu lượng khai thác thực tế tại giếng GD-1X và GD-2X đạt mức thương mại cao hơn 15% so với các giếng khoan ở khu vực đứt gãy bị phân cắt mạnh.

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

  • Tối ưu hóa quỹ đạo khoan cho các giếng phát triển từ GD-1P đến GD-8P bằng cách thiết kế góc nghiêng thân giếng đi qua tối đa các thể cát thuộc HU 3 và HU 4, nhằm duy trì sản lượng khai thác ban đầu trên 3.000 thùng dầu/ngày mỗi giếng trong giai đoạn 2013 - 2015, do ban Công nghệ mỏ Công ty Hoàng Long - Hoàn Vũ JOC chủ trì thực hiện.
  • Cập nhật mô hình biến sai Variogram và hiệu chỉnh thuộc tính rỗng thấm định kỳ 6 đến 12 tháng một lần bằng cách tích hợp dữ liệu logging trong khi khoan (LWD) và đo áp suất đáy giếng mới, mục tiêu giảm sai số dự báo trữ lượng xuống dưới 5%, do nhóm Địa chất phát triển mỏ thực hiện.
  • Xây dựng phương án bơm ép nước duy trì áp suất vỉa sau 18 tháng khai thác sơ cấp để ngăn chặn hiện tượng sụt giảm áp suất nhanh của khối H1 và H2, hướng tới mục tiêu nâng cao hệ số thu hồi dầu thêm từ 2% đến 4%, do kỹ sư mô phỏng mỏ PVEP POC đảm nhận.
  • Áp dụng quy trình tích hợp mô hình tướng dòng chảy HU và thuật toán địa thống kê SIS/SGS cho các tầng chứa sâu hơn như Hệ tầng Trà Tân trên (tập C) và Trà Tân giữa (tập D) trong chiến dịch khoan thăm dò mở rộng từ năm 2014, do Hội đồng Khoa học Tổng công ty Thăm dò Khai thác Dầu khí phê duyệt.

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

  • Kỹ sư địa chất và chuyên viên mô hình hóa mỏ: Nắm bắt phương pháp luận chi tiết về xây dựng mô hình cấu trúc 3D, ứng dụng biến sai Variogram và thuật toán SIS/SGS trong việc tái tạo tính bất đồng nhất của tầng chứa cát kết phức tạp.
  • Kỹ sư công nghệ mỏ và chuyên gia mô phỏng khai thác: Tham khảo kỹ thuật thô hóa mô hình mạng lưới (Upscaling) và quy trình lặp lại lịch sử thử vỉa (History Matching) để thiết lập mô hình động dự báo sản lượng chính xác.
  • Chuyên viên thẩm định dự án và quản lý tài nguyên dầu khí: Tiếp cận khung phương pháp luận chuẩn mực trong tính toán trữ lượng dầu tại chỗ (STOIIP), đánh giá rủi ro địa chất và xây dựng phương án phát triển kinh tế cho các mỏ dầu nhỏ cận biên ngoài khơi thềm lục địa.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Địa chất dầu khí: Sử dụng luận văn như một tài liệu học thuật và thực tiễn mẫu mực về tích hợp tài liệu địa chấn 3D PSDM, địa vật lý giếng khoan và mẫu lõi đặc biệt SCAL.

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

Phương pháp đơn vị dòng chảy thủy lực (HU) mang lại lợi ích gì so với mô hình tướng đá truyền thống? Phương pháp HU tích hợp đồng thời độ rỗng và độ thấm thông qua chỉ số FZI và RQI, giúp phân loại vỉa chứa thành các đới có đặc tính thủy lực đồng nhất. Cách tiếp cận này phản ánh chân thực tính bất đồng nhất của vỉa chứa cát kết phức tạp tại mỏ Gấu Đen, hạn chế sai số liên thông thủy lực khi đưa vào mô phỏng khai thác.

Tại sao luận văn sử dụng kết hợp hai thuật toán SIS và SGS trong mô hình hóa? Thuật toán SIS tối ưu cho việc mô phỏng các biến rời rạc như các loại tướng dòng chảy HU, bảo toàn hình thái địa chất của thể cát. Thuật toán SGS phù hợp cho các biến liên tục như độ rỗng và độ thấm, đảm bảo tái hiện chính xác quy luật phân bố xác suất không gian mà không làm mất đi các giá trị cực trị đo được tại giếng khoan.

Ranh giới dầu nước được xác định bằng phương pháp nào trong điều kiện vỉa chứa phân khối phức tạp? Ranh giới dầu nước được xác định chủ yếu thông qua minh giải đường phân bố áp suất chất lưu từ thiết bị đo RCI và MDT kết hợp các đường cong địa vật lý giếng khoan. Tại những khối vỉa tài liệu áp suất chưa rõ ràng, ranh giới được giả định khoa học tại vị trí trung bình giữa điểm gặp dầu sâu nhất và điểm tràn của cấu tạo bẫy chứa.

Quá trình thô hóa mô hình (Upscaling) được kiểm soát chất lượng như thế nào để không làm mất dữ liệu địa chất? Quá trình thô hóa áp dụng phương pháp trung bình số học cho độ rỗng và phương pháp trung bình điều hòa kết hợp hình học cho độ thấm. Chất lượng thô hóa được kiểm soát nghiêm ngặt thông qua việc so sánh biểu đồ phân bố tần suất và kiểm tra độ khớp lịch sử thử vỉa, đảm bảo giữ lại hơn 95% đặc trưng dị hướng của mô hình mịn ban đầu.

Trữ lượng 105 triệu thùng dầu tại chỗ đóng vai trò gì đối với chiến lược phát triển mỏ Gấu Đen? Con số 105 triệu thùng dầu tại chỗ xác nhận tiềm năng thương mại vững chắc của tầng chứa ILBH 5.2 tại khu vực phía Bắc. Đây là căn cứ kỹ thuật then chốt để các nhà đầu tư phê duyệt kế hoạch phát triển mỏ giai đoạn 1, quyết định vị trí hạ đặt giàn khai thác và tối ưu hóa số lượng giếng phát triển từ GD-1P đến GD-8P.

Kết luận

  • Xây dựng thành công mô hình cấu trúc 3 chiều hoàn chỉnh cho tầng chứa ILBH 5.2 mỏ Gấu Đen dựa trên 589 km² địa chấn 3D PSDM và số liệu giếng khoan.
  • Tiên phong ứng dụng phương pháp phân loại tướng dòng chảy thủy lực HU kết hợp chuỗi thuật toán địa thống kê SIS và SGS để mô hình hóa chính xác tính bất đồng nhất vỉa chứa.
  • Lượng hóa chính xác trữ lượng dầu tại chỗ đạt 105 triệu thùng, làm cơ sở khoa học cho công tác thiết kế khai thác mỏ.
  • Hoàn thành xuất sắc quá trình thô hóa mô hình mạng lưới và kiểm chứng thành công qua mô phỏng lặp lại lịch sử thử vỉa với độ tin cậy cao.
  • Cung cấp giải pháp kỹ thuật toàn diện phục vụ kế hoạch đưa mỏ Gấu Đen vào khai thác thương mại an toàn và đạt hiệu quả kinh tế tối ưu trong giai đoạn 2012 - 2015.

Đóng góp lớn nhất của luận văn là đã chuẩn hóa quy trình xây dựng mô hình địa chất 3 chiều tích hợp từ tài liệu địa chấn, địa vật lý giếng khoan và mẫu lõi cho các mỏ nhỏ cận biên tại thềm lục địa Việt Nam. Các kỹ sư và nhà nghiên cứu có thể khai thác phương pháp luận này để áp dụng cho các khối cấu trúc lân cận nhằm nâng cao hiệu quả quản trị và khai thác mỏ.