Tổng quan nghiên cứu

Bồn trũng Nam Côn Sơn có diện tích rộng lớn gần 100.000 km², nằm cách bờ biển Vũng Tàu khoảng 210 km về phía Đông Nam, là một trong những bể trầm tích có tiềm năng dầu khí quan trọng nhất thềm lục địa Việt Nam với tỷ lệ giếng khoan phát hiện đạt khoảng 35%. Trong bối cảnh các mỏ dầu khí lớn đang dần cạn kiệt, việc thăm dò và phát triển các mỏ biên có cấu trúc địa chất phức tạp như mỏ K tại Lô 10 & 11.1 trở thành nhiệm vụ chiến lược cấp thiết. Tầng Miocene giữa của mỏ K chứa tổng cộng 16 vỉa cát kết chứa dầu và khí riêng biệt, nhưng quy mô từng vỉa tương đối nhỏ với trữ lượng dầu tại chỗ tối đa chỉ khoảng 5 triệu thùng cho mỗi vỉa đơn lẻ.

Vấn đề kỹ thuật cốt lõi đặt ra là sự phân tán trữ lượng, tính bất đồng nhất cao của môi trường rỗng nứt nẻ và sự chênh lệch áp suất tự nhiên giữa các tập cát kết, khiến phương án khai thác đơn vỉa truyền thống không đạt hiệu quả kinh tế. Mục tiêu trọng tâm của nghiên cứu là ứng dụng kỹ thuật mô phỏng số thủy động lực học vỉa (Reservoir Simulation) bằng phần mềm chuyên dụng ECLIPSE để xây dựng mô hình mô phỏng 3D dòng chảy đa pha, tích hợp dữ liệu thử vỉa DST từ 2 giếng khoan thăm dò K-1X và K-2X, từ đó dự báo chính xác động thái khai thác và xác lập phương án khai thác hỗn hợp đa vỉa tối ưu.

Nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc quản trị rủi ro địa chất, kiểm soát hiện tượng dòng chảy xuyên giữa các tầng chứa, và đảm bảo mục tiêu thu hồi thương mại cho tổng trữ lượng dầu ban đầu tại chỗ cấp 2P đạt 9,20 triệu thùng (MMSTB) cùng trữ lượng khí tự do và hòa tan cấp 2P đạt 11,40 tỷ bộ khối (BCF), hướng tới hiện thực hóa hệ số thu hồi dầu 25% cho toàn mỏ.

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 cơ học chất lưu trong môi trường lỗ rỗng và lý thuyết mô phỏng vỉa dầu khí hiện đại, bao gồm:

  • Hệ phương trình dòng chảy đa pha: Tích hợp phương trình liên tục bảo toàn khối lượng, định luật Darcy mở rộng cho dòng chảy nhiều pha và phương trình Forchheimer cho dòng chảy vận tốc cao gần đáy giếng. Mối quan hệ tương hỗ giữa đất đá và chất lưu được mô tả qua phương trình nén ép thể tích đá kết hợp với các đường cong độ thấm tương đối hai pha dầu - nước và khí - dầu.
  • Mô hình chất lưu ba pha Black Oil: Mô tả các đặc tính nhiệt động học PVT của dầu, nước và khí thông qua các hàm số áp suất gồm hệ số thể tích thành hệ, độ nhớt động học, tỷ khối và tỷ số hòa tan khí trong dầu.
  • Mô hình giếng và khai thác hỗn hợp đa vỉa: Ứng dụng mô hình giếng Cullender - Smith kết hợp mô hình Leaky-tank để mô phỏng sự tương tác thủy lực giữa lòng giếng và các lớp thấm chứa có gradient áp suất khác nhau, làm sáng tỏ cơ chế xuất hiện hiệu ứng dòng thấm xuyên.
  • Các khái niệm chuyên ngành cốt lõi: Áp suất điểm bọt, độ bão hòa nước liên kết, độ rỗng hiệu dụng, hệ số thu hồi dầu và dòng chảy xuyên nội tại giữa các vỉa.

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

Nghiên cứu sử dụng nguồn dữ liệu địa chất, địa vật lý và công nghệ mỏ thu thập từ 2 giếng khoan thăm dò K-1X và K-2X, bao gồm 246 m tổng chiều dày cát kết hiệu dụng chứa hydrocarbon (trong đó giếng K-1X phát hiện 220 m vỉa từ độ sâu 3091 mMD đến 3680 mMD và giếng K-2X phát hiện 26 m). Phương pháp chọn mẫu là chọn mẫu toàn diện từ các tập mẫu lõi, tài liệu carota giếng khoan, kết quả đo áp suất thành hệ RFT, MDT và các đợt thử vỉa DST chuyên sâu tại các vỉa trọng điểm gồm LTS14, LTS16, LTS17, LTS18 và LTS19. Tiêu chuẩn lọc ngưỡng (cut-off) được xác lập chuẩn xác với thể tích sét dưới 40%, độ rỗng hiệu dụng tối thiểu 10% cho vỉa dầu, 7% cho vỉa khí và độ bão hòa nước dưới 70%.

Phương pháp phân tích chủ đạo là mô phỏng số rời rạc hóa không gian bằng lưới sai phân hữu hạn 3D kết hợp giải thuật lặp ngầm định hoàn toàn trên phần mềm ECLIPSE 100. Lý do lựa chọn phương pháp này là vì tính ưu việt vượt trội trong việc giải hệ phương trình vi phân đạo hàm riêng phi tuyến tính phức tạp, cho phép tái hiện chính xác sự biến thiên áp suất và độ bão hòa chất lưu theo không gian và thời gian trong điều kiện vỉa phân tầng bất đồng nhất cao, điều mà các phương pháp giải tích truyền thống hoàn toàn không thể đáp ứng. Tiến trình nghiên cứu được triển khai mạch lạc qua 5 giai đoạn: xác định mục tiêu, chuẩn hóa dữ liệu đầu vào, xây dựng mô hình số ban đầu, hiệu chỉnh lịch sử với số liệu thử vỉa DST và chạy các kịch bản dự báo khai thác cho giai đoạn phát triển mỏ kéo dài 10 đến 15 năm.

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

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

Quá trình mô phỏng số thủy động lực học tầng Miocene mỏ K đã mang lại các phát hiện khoa học và kỹ thuật nổi bật:

  • Định lượng chuẩn xác trữ lượng dầu khí tại chỗ: Kết quả mô phỏng ngẫu nhiên Monte Carlo kết hợp mô hình thể tích 3D xác định tổng trữ lượng dầu tại chỗ (OIIP) cấp 2P là 9,20 MMSTB (cấp 3P đạt 21,80 MMSTB) và trữ lượng khí ban đầu tại chỗ (GIIP) cấp 2P là 11,40 BCF (cấp 3P đạt 25,60 BCF). Với hệ số thu hồi áp dụng 25%, trữ lượng dầu có thể thu hồi cấp 2P đạt 2,30 MMSTB và cấp 3P đạt 5,45 MMSTB.
  • Xác nhận và lượng hóa hiệu ứng dòng thấm xuyên: Kết quả phân tích đối sánh áp suất vỉa ban đầu chỉ ra sự chênh lệch áp suất đáng kể giữa các tập cát kết độc lập (đặc biệt giữa vỉa LTS8, LTS14 và LTS16). Khi khai thác hỗn hợp qua giếng hoàn thiện kép, hiện tượng dòng thấm xuyên từ vỉa áp suất cao sang vỉa áp suất thấp xuất hiện rõ rệt trong các chu kỳ đóng giếng hoặc thay đổi lưu lượng, làm suy giảm từ 12% đến 18% sản lượng của các vỉa có áp suất đáy thấp hơn nếu không có van kiểm soát lưu lượng chuyên dụng.
  • Đánh giá hiệu quả kịch bản khai thác 2 giếng giả định: Mô hình dự báo sản lượng cho 2 giếng khai thác K-1P và K-2P cho thấy việc khai thác hỗn hợp đồng thời các vỉa trọng tâm cho phép thu hồi trên 85% tổng sản lượng dầu thương mại chỉ trong 5 năm đầu vận hành, đồng thời duy trì độ ngập nước toàn mỏ dưới ngưỡng 30% trong suốt giai đoạn khai thác sơ cấp.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến hiệu ứng dòng thấm xuyên là do sự dị thường về gradient áp suất thành hệ (cao hơn gradient thủy tĩnh chuẩn) kết hợp với độ dẫn thủy lực không đồng đều giữa các tập trầm tích cát kết xen kẹp sét kết tuổi Miocene giữa. Khi nhiều vỉa cùng mở thông vào một ống khai thác, dòng chất lưu có xu hướng tái phân bố theo gradient áp suất cục bộ, dẫn đến việc chất lưu từ vỉa có năng lượng cao thâm nhập ngược vào vỉa có năng lượng thấp hơn.

So sánh với các nghiên cứu tương tự tại các mỏ lân cận trong Bồn trũng Nam Côn Sơn như mỏ Đại Hùng, Rồng Đôi và các cấu tạo tại Lô 12W, hiện tượng dòng thấm xuyên luôn là thách thức lớn đối với các đối tượng vỉa mỏng đa tầng. Giải pháp khai thác hỗn hợp hoàn thiện kép được đề xuất trong luận văn hoàn toàn tương thích và chứng minh được tính ưu việt về mặt kinh tế kỹ thuật so với phương án khoan giếng đơn lẻ tốn kém.

Toàn bộ động thái thủy động lực học trong nghiên cứu được trực quan hóa sinh động thông qua hệ thống biểu đồ suy giảm lưu lượng dầu theo thời gian, đồ thị biến thiên áp suất đáy giếng theo từng phân tầng hoàn thiện và bảng cân bằng vật chất phân lớp chi tiết. Các biểu đồ mặt cắt phân bố độ bão hòa chất lưu 3D còn mô tả rõ nét tiến trình chuyển dịch của ranh giới dầu - nước và hiện tượng dâng nón nước tại đáy vỉa LTS16, cung cấp căn cứ vững chắc cho việc điều chỉnh chế độ khai thác thực tế.

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

Dựa trên kết quả mô phỏng số thủy động lực học, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm tối ưu hóa hoạt động phát triển mỏ K:

  • Thiết kế và triển khai công nghệ hoàn thiện giếng kép tích hợp van kiểm soát dòng thông minh: Ban Công nghệ mỏ và bộ phận thiết kế giếng khoan cần áp dụng cấu trúc hoàn thiện kép (Dual Completion) kết hợp van điều tiết lưu lượng đáy giếng (Inflow Control Valves) cho 2 giếng khai thác K-1P và K-2P ngay trong quý đầu tiên của giai đoạn thi công, nhằm mục tiêu triệt tiêu hoàn toàn hiện tượng dòng thấm xuyên và nâng cao hiệu suất khai thác thêm khoảng 8% đến 10%.
  • Lắp đặt thiết bị đo áp suất và nhiệt độ đáy giếng liên tục: Nhà điều hành mỏ cần trang bị hệ thống cảm biến đáy giếng vĩnh cửu (Permanent Downhole Gauges) với sai số cho phép dưới 0,5 psi trong quá trình hoàn thiện giếng, đảm bảo theo dõi liên tục sự suy giảm áp suất của từng vỉa theo thời gian thực và ngăn chặn rủi ro ngập nước sớm.
  • Thiết lập quy trình cập nhật và hiệu chỉnh mô hình động định kỳ: Nhóm kỹ sư thủy động lực học cần thực hiện cập nhật mô hình số ECLIPSE theo chu kỳ 6 tháng một lần dựa trên dữ liệu khai thác thực tế, hướng tới mục tiêu duy trì độ sai lệch giữa lưu lượng dự báo và thực tế dưới 5%.
  • Xây dựng phương án duy trì áp suất vỉa bằng bơm ép nước thứ cấp: Tổng công ty Thăm dò Khai thác Dầu khí (PVEP) cần lập kế hoạch thử nghiệm bơm ép nước sau 3 đến 5 năm khai thác sơ cấp, với mục tiêu duy trì áp suất vỉa trên mức áp suất điểm bọt và gia tăng hệ số thu hồi dầu toàn mỏ từ mức 25% lên khoảng 28% đến 30%.

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

Công trình nghiên cứu mang giá trị chuyên môn cao, là tài liệu tham khảo đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  • Kỹ sư công nghệ mỏ và kỹ sư mô phỏng thủy động lực học: Tiếp cận toàn diện quy trình thiết lập bộ từ khóa chuẩn mực trên phần mềm ECLIPSE (RUNSPEC, GRID, PROPS, SOLUTION, SUMMARY, SCHEDULE) và phương pháp mô hình hóa toán học cho bài toán khai thác vỉa nhiều lớp có xét đến dòng thấm xuyên.
  • Chuyên viên quản lý dự án và phát triển mỏ dầu khí: Sử dụng phương pháp luận và các chỉ số kinh tế kỹ thuật trong luận văn để xây dựng báo cáo kế hoạch phát triển mỏ (FDP), đánh giá rủi ro thương mại và lựa chọn số lượng giếng khoan tối ưu cho các mỏ biên quy mô nhỏ.
  • Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Địa chất - Dầu khí: Khai thác làm tài liệu giảng dạy, học tập chuyên sâu về phương pháp kết hợp mô phỏng xác suất Monte Carlo, minh giải thử vỉa DST và cơ học vỉa trầm tích Kainozoi thềm lục địa Việt Nam.
  • Chuyên gia phân tích địa chất và địa vật lý mỏ: Tham khảo bộ tiêu chuẩn ngưỡng định lượng (cut-off độ rỗng 10%, thể tích sét 40%), mô hình đường cong thấm tương đối và đặc trưng tướng đá cát kết tầng Miocene tại khu vực Lô 10 & 11.1 Bồn trũng Nam Côn Sơn.

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

Tại sao mỏ K bắt buộc phải áp dụng mô hình khai thác hỗn hợp đa vỉa thay vì khai thác từng vỉa đơn lẻ? Mỏ K bao gồm 16 vỉa cát kết riêng biệt nhưng mỗi vỉa có trữ lượng khá nhỏ, tối đa chỉ khoảng 5 triệu thùng dầu. Nếu đầu tư khoan giếng riêng lẻ cho từng vỉa, chi phí CAPEX và OPEX sẽ vượt quá giá trị kinh tế thu hồi. Khai thác hỗn hợp đồng thời nhiều vỉa qua giếng hoàn thiện kép là giải pháp duy nhất giúp tối ưu hóa chi phí đầu tư, đảm bảo dòng tiền dương và đạt tính khả thi thương mại cho mỏ.

Hiện tượng dòng thấm xuyên (cross-flow) gây ra những rủi ro kỹ thuật nào trong quá trình khai thác? Dòng thấm xuyên xuất hiện do sự chênh lệch áp suất tự nhiên giữa các vỉa chứa. Khi đóng giếng hoặc giảm lưu lượng, chất lưu từ vỉa áp suất cao sẽ chảy ngược vào vỉa áp suất thấp, gây tổn hại thành hệ, làm suy giảm năng lượng vỉa và giảm từ 12% đến 18% sản lượng của các tầng chứa yếu hơn nếu không được kiểm soát bằng van điều tiết lưu lượng chuyên dụng.

Dữ liệu thử vỉa DST đóng vai trò gì khi xây dựng mô hình số cho mỏ chưa đi vào khai thác? Khi mỏ chưa có lịch sử khai thác lâu dài, dữ liệu thử vỉa DST từ các giếng thăm dò K-1X và K-2X là nguồn số liệu thực tế duy nhất về lưu lượng tức thời, áp suất đáy giếng và độ dẫn thủy lực. Dữ liệu này đóng vai trò quyết định trong việc hiệu chỉnh các thông số độ thấm, độ dẫn của lưới mô phỏng ECLIPSE, giúp nâng cao độ tin cậy của mô hình dự báo.

Các thông số cắt mòn (cut-off) nào đã được áp dụng để xác định chiều dày cát chứa hiệu dụng? Nghiên cứu áp dụng các tiêu chuẩn cắt mòn tương đồng với các mỏ lân cận tại Bồn trũng Nam Côn Sơn gồm: hàm lượng sét tối đa 40%, độ rỗng hiệu dụng tối thiểu 10% đối với vỉa dầu và 7% đối với vỉa khí, cùng độ bão hòa nước cực đại 70%. Các ngưỡng này giúp xác định chính xác 246 m chiều dày cát kết chứa hydrocarbon hiệu dụng.

Mô hình mô phỏng trong luận văn có thể mở rộng ứng dụng cho các khu vực nào khác? Phương pháp luận xây dựng mô hình số và giải pháp khai thác hỗn hợp hoàn thiện kép trong nghiên cứu hoàn toàn có thể nhân rộng cho các mỏ dầu khí có đặc điểm địa chất phân tầng mỏng, cấu trúc phức tạp thuộc các lô lân cận trong Bồn trũng Nam Côn Sơn như Lô 12W, Lô 04, Lô 05 cũng như các bồn trũng Cửu Long và Sông Hồng.

Kết luận

  • Xây dựng thành công mô hình số thủy động lực học 3D hoàn chỉnh cho 16 vỉa cát kết tầng Miocene giữa mỏ K bằng phần mềm chuyên dụng ECLIPSE 100.
  • Lượng hóa chính xác trữ lượng dầu ban đầu tại chỗ cấp 2P đạt 9,20 MMSTB và khí đạt 11,40 BCF, đồng thời xác lập mức trữ lượng dầu thu hồi thương mại là 2,30 MMSTB với hệ số thu hồi 25%.
  • Làm sáng tỏ bản chất cơ học chất lưu và định lượng tác động của hiện tượng dòng thấm xuyên trong mô hình khai thác hỗn hợp đa vỉa phân lớp.
  • Khẳng định tính khả thi và hiệu quả kinh tế vượt trội của phương án bố trí 2 giếng khai thác hoàn thiện kép K-1P và K-2P trong việc khai thác tối ưu các vỉa trọng điểm LTS8, LTS14 và LTS16.
  • Cung cấp quy trình chuẩn mực về việc tích hợp dữ liệu thử vỉa DST để hiệu chỉnh mô hình vỉa cho các mỏ biên trước khi phê duyệt kế hoạch phát triển mỏ chính thức.

Đóng góp lớn nhất của luận văn là đã giải quyết trọn vẹn bài toán kỹ thuật phức tạp về khai thác hỗn hợp đa vỉa có hiệu ứng dòng xuyên, cung cấp cơ sở khoa học vững chắc giúp tối ưu hóa hiệu quả đầu tư và giảm thiểu rủi ro vận hành cho mỏ K. Trong kế hoạch tiếp theo kéo dài từ 12 đến 18 tháng tới, các đơn vị vận hành cần tiến hành cập nhật mô hình với số liệu khoan thực địa và thiết kế chi tiết hệ thống hoàn thiện giếng. Các tổ chức, kỹ sư công nghệ mỏ và nhà nghiên cứu dầu khí nên tham khảo và áp dụng ngay khung phương pháp luận này để nâng cao hiệu quả quản lý và phát triển các mỏ dầu khí cận biên tại Việt Nam.