Tổng quan nghiên cứu

Bồn trũng Cửu Long giữ vị trí chiến lược hàng đầu đối với an ninh năng lượng và sự phát triển kinh tế biển của Việt Nam. Với tổng diện tích khoảng 36.000 km2, bồn trũng này sở hữu tiềm năng tài nguyên ước tính trên 820 triệu m3 dầu quy đổi, trong đó hơn 500 triệu m3 dầu khí đã được phát hiện và khai thác thành công từ các bẫy cấu tạo và tầng móng nứt nẻ. Trong bối cảnh các mỏ dầu truyền thống đang bước vào giai đoạn suy giảm sản lượng tự nhiên, việc đẩy mạnh công tác tìm kiếm, thẩm lượng và đánh giá chính xác trữ lượng tại các tầng chứa mới đóng vai trò sống còn đối với sự phát triển bền vững của ngành công nghiệp dầu khí quốc gia.

Vấn đề cốt lõi đặt ra trong công tác địa chất dầu khí là tính bất định và độ rủi ro rất cao khi xác định thể tích dầu tại chỗ ban đầu tại các tầng chứa trầm tích cát kết vụn có cấu trúc phân lớp phức tạp. Đề tài tập trung giải quyết bài toán định lượng trữ lượng dầu khí cho tập E70 thuộc trầm tích Oligocen sớm tại mỏ X, Lô 09-2/09. Mục tiêu cụ thể của công trình là áp dụng phương pháp thể tích tĩnh kết hợp thuật toán mô phỏng xác suất thống kê nhằm xác lập giá trị trữ lượng dầu tại chỗ ban đầu ở các mức độ tin cậy khác nhau, từ đó cung cấp căn cứ khoa học vững chắc để các nhà điều hành đưa ra quyết định tiếp tục đầu tư khoan mở rộng hoặc xây dựng kế hoạch phát triển mỏ thương mại.

Phạm vi nghiên cứu được thực hiện trên cấu tạo mỏ X, nằm ở phía Bắc Lô 09-2/09 thềm lục địa phía Nam, cách bờ biển Vũng Tàu khoảng 135 km về phía Đông Nam và cách mỏ Rạng Đông khoảng 10 km, tại vùng biển có độ sâu nước dao động từ 50 m đến 70 m. Dữ liệu nghiên cứu bao quát chuỗi khảo sát địa chấn và giếng khoan trong giai đoạn từ năm 2010 đến năm 2018. Kết quả đánh giá cung cấp các chỉ số trữ lượng ở các mức xác suất P90, P50 và P10, giúp tối ưu hóa hệ số thu hồi dầu dự kiến đạt 16% cho tập vỉa E70, đồng thời giảm thiểu đáng kể rủi ro tài chính trị giá hàng chục triệu USD trong các chiến dịch khoan khai thác tiếp theo.

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 hệ thống lý thuyết phân cấp tài nguyên và trữ lượng dầu khí theo tiêu chuẩn quốc tế và quy định hiện hành của ngành dầu khí Việt Nam. Hệ thống này phân định rõ ràng giữa tài nguyên chưa phát hiện và trữ lượng đã phát hiện, trong đó trữ lượng còn lại được phân hạng theo ba cấp độ tin cậy: trữ lượng xác minh 1P với xác suất 90%, trữ lượng có khả năng 2P với xác suất 50% và trữ lượng có thể 3P với xác suất 10%. Ranh giới phân cấp cho các thân dầu được xác định chặt chẽ dựa trên nguyên tắc chia đôi khoảng cách giữa các ranh giới chất lưu đã ghi nhận như điểm dầu cao nhất, điểm khí thấp nhất và điểm tiếp xúc dầu nước.

Cơ sở tính toán cốt lõi là lý thuyết phương pháp thể tích, cho phép xác định tổng lượng dầu tại chỗ ban đầu dựa trên mối tương quan giữa các tham số hình học vỉa và đặc tính thạch học. Các khái niệm và thông số cơ bản bao gồm: thể tích khối đá chứa sản phẩm, tỷ lệ cát hiệu dụng trong vỉa cát sét kết, độ rỗng hiệu dụng của đá chứa, độ bão hòa dầu ban đầu và hệ số thể tích thành hệ của dầu ở điều kiện áp suất vỉa. Bên cạnh đó, luận văn áp dụng lý thuyết xác suất thống kê hiện đại để chuyển đổi các biến số địa chất từ dạng tất định đơn lẻ sang dạng phân bố xác suất liên tục, giúp mô hình hóa toàn diện các yếu tố bất định của môi trường lòng đất.

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

Nguồn dữ liệu đầu vào của đề tài bao gồm khối dữ liệu địa chấn 3D phủ diện tích 500 km2 thu nổ năm 2010 kết hợp với 410 km2 tài liệu địa chấn xử lý năm 2001, cùng tập hợp tài liệu địa vật lý giếng khoan và mẫu lõi từ 3 giếng khoan thăm dò X-1X, X-2X và X-3X. Cỡ mẫu nghiên cứu bao gồm toàn bộ chuỗi số liệu đo carota liên tục qua các tầng địa tầng và 2 mẫu lõi khoan nguyên dạng với tổng chiều dài 37 m thu thập tại giếng X-2X và X-3X. Phương pháp chọn mẫu là chọn mẫu tầng thứ có chủ đích, tập trung toàn diện vào tập trầm tích cát kết E70 ở độ sâu từ 3240 mss đến 3650 mss nhằm phản ánh chính xác nhất đặc tính cơ lý và chất lưu của vỉa.

Quy trình phân tích được thực hiện qua các bước tích hợp công nghệ cao. Trước hết, mô hình cấu trúc 3D được xây dựng trên phần mềm Petrel để minh giải chính xác bề mặt nóc vỉa, đáy vỉa và tính toán thể tích đá chứa. Tiếp theo, phần mềm Bestfit được sử dụng để phân tích kiểm định thống kê và lựa chọn quy luật phân bố xác suất tối ưu cho từng thông số vỉa. Cuối cùng, thuật toán Monte Carlo được triển khai trên phần mềm Crystal Ball với 10.000 vòng lặp ngẫu nhiên để xuất ra đường cong kỳ vọng trữ lượng và phân tích độ nhạy. Lý do lựa chọn tổ hợp phương pháp này là vì tại giai đoạn đầu của mỏ khi chưa có số liệu khai thác dài hạn để áp dụng phương pháp cân bằng vật chất hay đường cong suy giảm sản lượng, phương pháp thể tích kết hợp mô phỏng ngẫu nhiên là công cụ duy nhất đảm bảo tính chính xác và phản ánh trung thực mọi rủi ro địa chất. Toàn bộ chuỗi nghiên cứu và xử lý dữ liệu được tiến hành liên tục từ khi thu thập mẫu đến khi hoàn thiện vào tháng 7 năm 2018.

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

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

Quá trình xử lý và mô hình hóa dữ liệu địa chất đã mang lại 4 phát hiện quan trọng về tập chứa E70 mỏ X:

Thứ nhất, về mặt cấu trúc địa chất, tập E70 được xác định là một nếp lồi dạng vòm khép kín 3 chiều, có dạng hình thoi kéo dài theo phương Đông Bắc - Tây Nam, bị giới hạn bởi đứt gãy chính cùng phương. Diện tích cấu tạo của vỉa đạt 10,9 km2 với biên độ khép kín thẳng đứng lên tới 410 m, trải dài từ đỉnh vỉa ở độ sâu 3240 mss đến điểm tràn cấu tạo ở 3650 mss.

Thứ hai, đặc tính thạch học và thấm chứa của tập E70 thể hiện tiềm năng công nghiệp rất tích cực. Kết quả phân tích mẫu lõi và carota ghi nhận độ rỗng hiệu dụng trung bình dao động từ 10% đến 14%, độ thấm đạt từ hàng chục đến hàng trăm miliDarcy. Thử vỉa DST số 2 tại giếng X-3X cho lưu lượng dòng dầu tự nhiên rất mạnh, đạt từ 1400 đến 2300 thùng/ngày với tỷ trọng dầu nhẹ đạt 35 độ API, cao hơn đáng kể so với dòng dầu 536 thùng/ngày ghi nhận tại tập E40 của giếng X-1X.

Thứ ba, kết quả mô phỏng Monte Carlo xác định trữ lượng dầu tại chỗ ban đầu cho tập E70 tại các mức phân vị xác suất cụ thể: mức P90 tương ứng cấp 1P đạt 28,45 triệu thùng, mức P50 tương ứng cấp 2P đạt 45,82 triệu thùng (tăng trưởng 61,05% so với mức P90), và mức P10 tương ứng cấp 3P đạt 68,16 triệu thùng. Trữ lượng này đóng góp tỷ trọng then chốt vào tổng thể tích dầu tại chỗ ban đầu của toàn bộ mỏ X với mức 2P đạt 158,3 triệu thùng.

Thứ tư, phân tích độ nhạy biến thiên cho thấy thể tích khối đá chứa và độ bão hòa nước ban đầu là hai yếu tố chi phối mạnh mẽ nhất đến kết quả trữ lượng, đóng góp lần lượt khoảng 42% và 26% vào tổng phương sai biến thiên của mô hình.

Thảo luận kết quả

Nguyên nhân tạo nên đặc tính chứa tốt nhưng có sự phân dị lớn của tập E70 bắt nguồn từ môi trường lắng đọng trầm tích đồng bằng châu thổ sông hồ vào thời kỳ Oligocen sớm. Sự biến thiên thạch học giữa các phân lớp cát kết hạt mịn - trung bình xen kẹp sét kết giàu vật liệu hữu cơ đã tạo ra tính bất đồng nhất cao trong không gian vỉa chứa. So sánh với các nghiên cứu trước đây tại các mỏ lân cận trong bồn trũng Cửu Long như mỏ Rạng Đông hay mỏ Cá Ngừ Vàng, các tham số vỉa của tập E70 có độ tương đồng cao về cơ chế bẫy cấu tạo kết hợp nứt nẻ, nhưng cho thấy mức độ bảo tồn độ rỗng nguyên sinh vượt trội hơn khoảng 15% so với các tầng cát kết sâu cùng tuổi.

Về mặt biểu diễn trực quan, toàn bộ kết quả phân tích thống kê và mô phỏng được trình bày rất rõ ràng qua hệ thống bảng biểu đối chiếu thông số đầu vào và các biểu đồ chuyên dụng. Phân bố xác suất tích lũy của trữ lượng dầu được minh họa sinh động bằng biểu đồ hàm phân bố tích lũy hình chữ S, giúp phân định ranh giới giữa các mức P90, P50 và P10 một cách trực quan. Đồng thời, biểu đồ lốc xoáy phân tích độ nhạy thể hiện rõ thứ tự tác động của từng biến số địa chất từ cao xuống thấp, cho phép các kỹ sư nhận diện chính xác các nguồn gây bất định để tập trung xử lý trong các bước thẩm lượng tương lai. Kết quả này chứng minh phương pháp thể tích kết hợp mô phỏng xác suất có độ tin cậy vượt trội so với cách tính tất định truyền thống khi sai số hội tụ được kiểm soát dưới mức 6%.

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

Dựa trên các kết quả nghiên cứu và phân tích độ nhạy của vỉa chứa E70 mỏ X, bốn nhóm giải pháp kỹ thuật và quản lý được đề xuất cụ thể như sau:

Một là, tiến hành khoan thẩm lượng bổ sung 01 đến 02 giếng khoan định hướng tại khu vực trung tâm và sườn Đông Bắc của cấu tạo mỏ X trong giai đoạn 2019 đến 2021. Mục tiêu hành động nhằm xác định chính xác ranh giới tiếp xúc dầu nước và sự phát triển theo phương ngang của thân cát kết, hướng tới chuyển đổi tối thiểu 40% trữ lượng từ cấp P2 sang cấp P1 xác minh. Chủ thể chịu trách nhiệm thực thi là Tổng công ty Thăm dò Khai thác Dầu khí cùng Công ty Điều hành chung.

Hai là, áp dụng giải pháp nứt vỉa thủy lực có kiểm soát và xử lý axít hóa vùng cận đáy giếng cho các phân lớp cát kết có độ thấm dưới 20 miliDarcy. Giải pháp này dự kiến nâng hệ số năng suất vỉa lên từ 25% đến 35%, gia tăng lưu lượng khai thác đơn giếng đạt mức trên 2000 thùng/ngày trong giai đoạn đầu vận hành. Bộ phận công nghệ khai thác và nhà thầu dịch vụ giếng khoan sẽ phối hợp triển khai ngay khi hoàn thiện giếng.

Ba là, xây dựng mô hình mô phỏng khai thác động 3D tích hợp toàn mỏ trong vòng 9 tháng sau khi hoàn tất chiến dịch khoan thẩm lượng. Mô hình này sẽ tối ưu hóa số lượng và vị trí đặt các giếng khai thác ngầm, kết nối hiệu quả với hệ thống thu gom xử lý nổi nhằm giảm 15% chi phí đầu tư xây dựng cơ bản. Nhóm kỹ sư công nghệ mỏ và địa chất chịu trách nhiệm xây dựng và cập nhật mô hình.

Bốn là, chuẩn bị phương án bơm ép nước duy trì áp suất vỉa ngay từ năm thứ hai của quá trình khai thác thương mại. Biện pháp này nhằm bảo tồn năng lượng tầng chứa, ngăn chặn hiện tượng sụt giảm áp suất sớm, qua đó nâng hệ số thu hồi dầu của tập E70 từ mức cơ sở 16% lên mục tiêu 22% đến 25% trong suốt vòng đời mỏ kéo dài 15 đến 20 năm. Đơn vị điều hành khai thác mỏ là cơ quan trực tiếp lập và giám sát phương án vận hành này.

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

Luận văn mang lại giá trị học thuật và thực tiễn sâu sắc cho 4 nhóm đối tượng chính sau:

Nhóm thứ nhất là các kỹ sư địa chất, kỹ sư địa vật lý và chuyên gia công nghệ mỏ đang công tác tại các công ty dầu khí như Tổng công ty Thăm dò Khai thác Dầu khí, Liên doanh Việt - Nga Vietsovpetro hay các công ty điều hành liên doanh. Luận văn cung cấp quy trình mẫu về việc minh giải số liệu địa chấn 3D, chuẩn hóa log giếng khoan và kỹ thuật phân bố xác suất tham số vỉa để ứng dụng trực tiếp vào các đề án thẩm lượng thực tế.

Nhóm thứ hai là các nhà quản lý dự án, chuyên viên kinh tế năng lượng và chuyên viên thẩm định đầu tư mỏ dầu khí. Nghiên cứu mang đến khung đánh giá định lượng rủi ro thông qua các kịch bản trữ lượng P10, P50, P90 và biểu đồ lốc xoáy độ nhạy, giúp các nhà quản trị đưa ra quyết định phê duyệt kế hoạch đại cương phát triển mỏ hoặc phân bổ vốn đầu tư hàng trăm triệu USD một cách an toàn và hiệu quả.

Nhóm thứ ba là giảng viên, nghiên cứu sinh và học viên cao học thuộc các chuyên ngành Kỹ thuật Dầu khí, Địa chất Dầu khí và Quản trị Tài nguyên Năng lượng tại các trường đại học kỹ thuật. Đây là tài liệu tham khảo chuẩn mực về phương pháp luận kết hợp liên hoàn giữa các công cụ phần mềm chuyên dụng hàng đầu thế giới như Petrel, Bestfit và Crystal Ball.

Nhóm thứ tư là các cơ quan quản lý nhà nước về tài nguyên khoáng sản như Tập đoàn Dầu khí Việt Nam và Bộ Công Thương. Luận văn là cơ sở tham chiếu tin cậy để đối soát, kiểm tra và phê duyệt các báo cáo tài nguyên trữ lượng dầu khí theo đúng quy chuẩn phân cấp kỹ thuật quốc gia.

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

Phương pháp thể tích kết hợp mô phỏng xác suất Monte Carlo mang lại ưu thế gì vượt trội cho tập E70 mỏ X? Trong giai đoạn đầu tìm kiếm thăm dò, mỏ X chưa đi vào khai thác nên hoàn toàn thiếu dữ liệu lịch sử sản lượng để áp dụng phương pháp cân bằng vật chất hay phân tích đường cong suy giảm. Việc kết hợp phương pháp thể tích với thuật toán Monte Carlo giúp khắc phục triệt để hạn chế này bằng cách xử lý đồng thời tính bất định của toàn bộ các biến số địa chất như diện tích 10,9 km2, độ rỗng 10% đến 14% qua hàng chục nghìn phép thử ngẫu nhiên, tạo ra dải kết quả có độ tin cậy khoa học cao.

Ý nghĩa thực tiễn của các con số trữ lượng 1P, 2P và 3P đối với việc phát triển mỏ là gì? Các cấp trữ lượng 1P với xác suất 90%, 2P với xác suất 50% và 3P với xác suất 10% đóng vai trò quyết định trong việc hoạch định kinh tế mỏ. Trong đó, mức trữ lượng 2P với 158,3 triệu thùng dầu tại chỗ toàn mỏ và 45,82 triệu thùng tại tập E70 là con số chuẩn mực được các ngân hàng và tổ chức tài chính sử dụng để thẩm định tính khả thi thương mại, từ đó phê duyệt kế hoạch tài chính và thiết kế kích thước giàn công nghệ khai thác.

Bộ dữ liệu mẫu lõi khoan 37 m từ giếng X-2X và X-3X có vai trò như thế nào trong nghiên cứu? Hai mẫu lõi khoan với tổng chiều dài 37 m cung cấp các thông số thực nghiệm đo đạc chính xác về độ rỗng, độ thấm không khí, áp suất mao dẫn và độ bão hòa nước liên kết. Bộ dữ liệu vật lý trực tiếp này đóng vai trò chuẩn hóa và hiệu chỉnh các đường cong địa vật lý giếng khoan, giúp nâng cao độ chính xác khi xác định ranh giới phân lớp và tỷ lệ cát hiệu dụng trong mô hình địa chất 3D.

Thông số địa chất nào ảnh hưởng lớn nhất đến sự biến thiên của trữ lượng dầu tại chỗ tập E70? Kết quả phân tích độ nhạy trên biểu đồ lốc xoáy chỉ ra rằng thể tích khối đá chứa sản phẩm và độ bão hòa nước ban đầu là hai yếu tố chi phối mạnh mẽ nhất, đóng góp hơn 65% vào tổng mức độ biến động của kết quả trữ lượng. Điều này phản ánh rõ nét tính chất phức tạp của cấu trúc vát mỏng vỉa chứa và sự hiện diện của đới chuyển tiếp nước dầu trong trầm tích Oligocen.

Làm thế nào để hiện thực hóa và gia tăng hệ số thu hồi dầu 16% của tập E70 trong thực tế? Hệ số thu hồi 16% được xác lập dựa trên đối sánh đặc tính vỉa tương đồng với các mỏ lân cận như Bạch Hổ và Rạng Đông. Để đạt và nâng cao chỉ số này lên mức 20% đến 25%, đơn vị khai thác cần đồng bộ triển khai kỹ thuật nứt vỉa thủy lực tăng cường lưu lượng đáy giếng và thực hiện bơm ép nước duy trì áp suất vỉa ngay từ giai đoạn đầu khi áp suất thành hệ bắt đầu suy giảm.

Kết luận

Nghiên cứu đánh giá trữ lượng dầu tại chỗ ban đầu tập E70 mỏ X Lô 09-2/09 bằng phương pháp thể tích đã đạt được những kết luận then chốt sau:

  • Xác lập hoàn chỉnh mô hình cấu trúc địa chất 3D của tập vỉa E70 với dạng nếp lồi vòm khép kín diện tích 10,9 km2 và chiều cao khép kín 410 m.
  • Định lượng chính xác các thông số vật lý thạch học vỉa chứa với độ rỗng hiệu dụng từ 10% đến 14%, độ thấm tốt và tỷ trọng dầu nhẹ đạt 35 độ API.
  • Xác định phổ trữ lượng dầu tại chỗ ban đầu bằng mô phỏng ngẫu nhiên Monte Carlo với mức 1P đạt 28,45 triệu thùng, mức 2P đạt 45,82 triệu thùng và mức 3P đạt 68,16 triệu thùng.
  • Nhận diện thành công thể tích khối đá chứa và độ bão hòa nước là hai tham số nhạy cảm trọng yếu chi phối độ bất định của toàn bộ mô hình.
  • Đề xuất hệ thống giải pháp kỹ thuật đồng bộ bao gồm khoan thẩm lượng mở rộng và bơm ép duy trì áp suất vỉa nhằm nâng cao hệ số thu hồi dầu từ 16% lên 25%.

Đóng góp lớn nhất của công trình là xây dựng thành công quy trình đánh giá trữ lượng khoa học, tích hợp chặt chẽ giữa địa chấn 3D, mô hình Petrel và phân tích bất định Crystal Ball, mang lại lời giải tin cậy cho bài toán thăm dò vỉa chứa phức tạp tại bồn trũng Cửu Long. Trong lộ trình 12 đến 24 tháng tới, việc triển khai các giếng khoan thẩm lượng tiếp theo sẽ là bước đi quyết định để chính xác hóa ranh giới mỏ. Độc giả, các kỹ sư và nhà nghiên cứu quan tâm được khuyến khích tham khảo chi tiết toàn bộ luận văn để ứng dụng hiệu quả vào các dự án đánh giá tài nguyên dầu khí tương tự.