Tổng quan nghiên cứu

Nằm cách bờ biển Vũng Tàu khoảng 320 km về phía Tây Nam thuộc bồn trũng Nam Côn Sơn, cụm mỏ Rồng Đôi - Rồng Đôi Tây giữ vị thế chiến lược trong hạ tầng an ninh năng lượng quốc gia kể từ khi chính thức đón dòng sản phẩm thương mại đầu tiên vào ngày 25 tháng 12 năm 2006. Dự án được thiết kế với công suất vận chuyển khí trung bình 130 triệu bộ khối khí tiêu chuẩn mỗi ngày (MMscfd) và có khả năng đáp ứng lưu lượng đỉnh ngắn hạn đạt 175 MMscfd. Thách thức kỹ thuật lớn nhất của cụm mỏ bắt nguồn từ cấu trúc địa tầng Miocene phân dị phức tạp với 15 vỉa cát mỏng chứa khí và condensate, đòi hỏi việc dự báo lưu lượng khai thác phải đạt độ chính xác cao nhằm tuân thủ nghiêm ngặt các điều khoản trong hợp đồng mua bán khí (GSPA) ký kết với Tập đoàn Dầu khí Việt Nam.

Mục tiêu trọng tâm của nghiên cứu là xây dựng quy trình phân tích và dự báo khả năng khai thác toàn diện cho toàn bộ 6 giếng khai thác hiện hữu (gồm các giếng RT-1P, RT-2P, RT-3P, RT-4P, RT-5P tại mỏ Rồng Đôi và giếng RDT-1P tại mỏ Rồng Đôi Tây) xuyên suốt vòng đời kinh tế 30 năm của dự án. Nghiên cứu khai thác chuỗi dữ liệu vận hành thực tế từ tháng 10 năm 2006 đến hết tháng 12 năm 2010 tại giàn đầu giếng (WHd) và cụm giàn xử lý trung tâm PUQC để lượng hóa các quy luật suy giảm sản lượng.

Về mặt ý nghĩa thực tiễn, việc ứng dụng các mô hình suy giảm chuẩn xác giúp các nhà điều hành duy trì ổn định lưu lượng xuất bán đạt trên 3,77 triệu mét khối khí tiêu chuẩn mỗi ngày (MMSCMD) và hơn 4.160 thùng condensate mỗi ngày (BCPD). Kết quả nghiên cứu cung cấp cơ sở khoa học vững chắc để tối ưu hóa thời điểm can thiệp giếng, lập kế hoạch bắn vỉa bổ sung, khoan nhánh và chuẩn bị lắp đặt hệ thống máy nén khí phụ trợ, qua đó nâng cao hệ số thu hồi hydrocacbon lên trên 65% tổng trữ lượng tại chỗ.

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 lý thuyết phân tích đường cong suy giảm sản lượng (Decline Curve Analysis - DCA) kinh điển của Arps (1945), kết hợp với các mô hình dòng chảy chất lưu và cân bằng vật chất trong vỉa chứa dầu khí do Fetkovich và Horne phát triển. Phương trình vi phân tổng quát của Arps liên kết lưu lượng khai thác tại thời điểm t với lưu lượng ban đầu, tốc độ suy giảm danh nghĩa và hệ số suy giảm b nhằm mô tả chính xác cơ chế động học của vỉa.

Ba dạng đường cong suy giảm cốt lõi được áp dụng trong khung lý thuyết gồm: suy giảm hàm mũ (khi hệ số b bằng 0), suy giảm hyperbol (khi hệ số b nằm trong khoảng từ 0 đến 1, đặc biệt b bằng 0,3 tương ứng với cơ chế khí hòa tan và b bằng 0,5 đối với cơ chế mũ khí hoặc phân dị trọng lực) và suy giảm điều hòa (khi hệ số b bằng 1 tương ứng với cơ chế ngập nước mạnh). Bên cạnh đó, nghiên cứu vận dụng 5 khái niệm chuyên ngành then chốt: tỷ lệ chất lỏng condensate (CLR), độ suy giảm danh nghĩa ban đầu (Di), lưu lượng giới hạn kinh tế cho phép (qa), áp suất vỉa ban đầu dao động từ 3.780 đến 5.426 psia và hệ số thành hệ thể tích chất lưu (Bo). Các mô hình này cho phép ngoại suy sản lượng cộng dồn hữu hạn và xác định chính xác thời điểm giếng chạm ngưỡng giới hạn khai thác.

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

Nguồn dữ liệu đầu vào của đề tài được trích xuất trực tiếp từ hệ thống điều khiển và thu thập dữ liệu giám sát Delta V (DCS) kết hợp cùng báo cáo sản lượng hàng ngày và hàng tháng từ tháng 10 năm 2006 đến tháng 12 năm 2010, tạo nên chuỗi dữ liệu chu kỳ hơn 50 tháng liên tục. Cỡ mẫu nghiên cứu bao trùm 100% số giếng đang hoạt động tại mỏ, gồm 6 giếng khai thác hoàn thiện hỗn hợp (commingled completion) đại diện cho 15 vỉa cát chứa khí phân tầng dày từ 4 đến 15 mét. Phương pháp chọn mẫu phi xác suất toàn diện (purposive census sampling) được áp dụng nhằm đảm bảo phản ánh đầy đủ mọi biến số địa chất và áp suất vận hành của toàn bộ cụm mỏ ngoài khơi.

Quy trình phân tích tích hợp phương pháp tiếp cận kép: xử lý giải tích trên Microsoft Excel và mô phỏng động học chuyên sâu bằng phần mềm MBAL thuộc bộ công cụ tích hợp quản lý mỏ IPM của Petroleum Experts. Lý do lựa chọn giải pháp phân tích này xuất phát từ nhu cầu kiểm chứng chéo độ tin cậy. Bảng tính Excel cung cấp khả năng trực quan hóa các đồ thị bán logarit và log-log giữa CLR, lưu lượng và sản lượng cộng dồn nhằm làm sáng tỏ bản chất vật lý của dòng chảy, trong khi phần mềm MBAL đảm bảo tính toán cân bằng vật chất phi tuyến tính, tích hợp chính xác các thông số nhiệt động học PVT và áp lực vận hành cụm phân dòng từ 80 đến 100 bar. Timeline phân tích được thực hiện qua 4 bước tuần tự từ nhận diện dạng suy giảm, xác định thông số b và Di, đến khớp lịch sử khai thác và ngoại suy tương lai.

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 lịch sử khai thác chỉ ra rằng trong giai đoạn đầu khi áp suất vỉa duy trì ở mức cao từ 4.500 đến 5.426 psia, các giếng như RT-1P và RT-4P suy giảm chủ yếu theo dạng hàm mũ với hệ số b bằng 0. Khi quá trình khai thác kéo dài, áp suất đáy giếng suy giảm và xuất hiện sự phân dị chất lưu, đường cong chuyển dịch rõ rệt sang dạng hyperbol với hệ số b dao động từ 0,3 đến 0,5.

Thứ hai, sản lượng khai thác thực tế của toàn mỏ ghi nhận bước tăng trưởng vượt bậc từ 2,22 MMSCMD khí và 2.685 BCPD condensate (khi vận hành 2 giếng ban đầu vào đầu năm 2007) lên mức ổn định 3,83 MMSCMD khí và 4.163 BCPD condensate vào tháng 4 năm 2008 khi kết nối đủ 6 giếng, tương ứng mức tăng trưởng sản lượng chất lỏng hơn 55%.

Thứ ba, các thông số suy giảm ban đầu Di được lượng hóa chính xác cho từng giếng, dao động trong biên độ từ 15,2% đến 28,4%/năm. Sự suy giảm tỷ lệ condensate (CLR) diễn ra nhanh hơn ở các giếng khai thác tầng cát sâu do hiện tượng ngưng tụ ngược trong vỉa khi áp suất tầng chứa giảm xuống dưới điểm sương.

Thứ tư, sự kết hợp giữa mô phỏng MBAL và phân tích Excel đạt hệ số tương quan khớp lịch sử khai thác vượt trên 92%, chứng minh độ tin cậy cao của bộ thông số dự báo cho phép xác định chính xác thời điểm từng giếng chạm lưu lượng giới hạn kinh tế.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự biến thiên của các quy luật suy giảm tại mỏ Rồng Đôi - Rồng Đôi Tây bắt nguồn từ đặc tính trầm tích điệp Dừa và điệp Thông thuộc thế Miocene sớm đến giữa. Sự xen kẹp dày đặc giữa các lớp cát kết mỏng từ 4 đến 15 mét với các vỉa phiến sét tạo ra sự không đồng nhất lớn về độ thấm, trải rộng từ 1,3 đến 250 md và độ rỗng biến thiên từ 12,5% đến 21,3%. Kết quả này hoàn toàn tương thích với các nghiên cứu của Horne (2003) và Fetkovich về vỉa nứt nẻ và phân tầng, củng cố nhận định rằng các vỉa khí có độ thấm thấp và phân tầng mạnh luôn có xu hướng suy giảm theo dạng hyperbol thay vì hàm mũ thuần túy.

Trong thực tế trình bày báo cáo kỹ thuật, dữ liệu khai thác và dự báo được trực quan hóa tối ưu thông qua hệ thống 6 biểu đồ chuyên dụng: biểu đồ lưu lượng theo thời gian, biểu đồ bán logarit Log(qt) theo thời gian, đồ thị CLR so với sản lượng cộng dồn Np, và đồ thị Log(CLR) so với Np. Trên đồ thị bán logarit, đường xu hướng thể hiện sự thay đổi độ dốc rõ nét tại các mốc giếng bắt đầu xuất hiện nước đồng hành. Bảng đối chiếu tham số giữa MBAL và Excel cho thấy sự nhất quán cao về giá trị Di và b, song mô hình MBAL mang lại ưu thế vượt trội khi cho phép mô phỏng linh hoạt áp suất tại chạc ba đường ống Nam Côn Sơn từ 80 đến 157 bar, giúp dự báo chính xác áp lực đầu giếng theo từng giai đoạn suy giảm năng lượng vỉa.

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

Thứ nhất, ban điều hành mỏ cần thiết lập quy trình kiểm tra lưu lượng định kỳ qua bình tách thử (Test Separator với công suất xử lý 50 triệu bộ khối khí/ngày) với tần suất tối thiểu 14 ngày một lần cho từng giếng. Hoạt động này do Đội ngũ Kỹ sư Công nghệ mỏ chủ trì thực hiện nhằm cập nhật chỉ số suy giảm Di và hệ số b vào cơ sở dữ liệu giám sát trong quý 1 hàng năm, đảm bảo sai số dự báo không vượt quá 5%.

Thứ hai, tiến hành triển khai dự án lắp đặt hệ thống máy nén khí phụ trợ (booster compressor) trên giàn PUQC vào mốc thời gian sau 14 năm khai thác thương mại. Ban Quản lý Kỹ thuật và Dự án cần hoàn tất báo cáo khả thi chi tiết trước năm vận hành thứ 10 để duy trì áp suất xuất khí ổn định ở mức 80 đến 160 barg, đáp ứng yêu cầu áp lực tiếp nhận tại đường ống Nam Côn Sơn khi áp suất đáy giếng suy giảm.

Thứ ba, thực hiện các chiến dịch can thiệp giếng chuyên sâu bằng công nghệ ống cuộn xoắn (Coiled Tubing) và bắn vỉa mở rộng tại các tập cát mỏng chưa khai thác triệt để thuộc các vỉa N/O, P và S1. Phòng Khoan và Hoàn thiện giếng cần chủ động lập phương án kỹ thuật trong giai đoạn 2012 - 2015 nhằm mục tiêu gia tăng hệ số thu hồi thêm từ 8% đến 12% sản lượng condensate tích lũy.

Thứ tư, nâng cấp hệ thống kết nối dữ liệu tự động giữa phần mềm SCADA Delta V và công cụ mô phỏng động lực học MBAL để thực hiện hiệu chỉnh lịch sử khai thác theo thời gian thực 24/7. Nhóm Kỹ sư Khai thác phối hợp bộ phận CNTT vận hành hệ thống này nhằm kiểm soát chặt chẽ tỷ lệ ngập nước (Water-cut) luôn ở mức an toàn dưới 50%, ngăn ngừa nguy cơ ngập úng giếng gây dừng khai thác đột ngột.

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

Nhóm kỹ sư công nghệ mỏ và kỹ sư khai thác dầu khí: Luận văn cung cấp cẩm nang thực hành chi tiết về ứng dụng DCA và phần mềm MBAL để xử lý chuỗi dữ liệu khai thác thực tế, hỗ trợ trực tiếp công tác lập kế hoạch can thiệp giếng và tối ưu hóa chế độ mở côn cho từng giếng khoan ngoài khơi.

Nhóm nhà quản lý dự án và chuyên viên phân tích thương mại năng lượng: Cung cấp các công cụ định lượng chuẩn xác về sản lượng khí và condensate theo từng mốc thời gian, làm cơ sở xây dựng kế hoạch dòng tiền, quản trị rủi ro hợp đồng mua bán khí GSPA và hoạch định chi phí đầu tư giàn khai thác.

Nhóm giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Kỹ thuật Dầu khí: Tài liệu là nguồn tham khảo học thuật giá trị cao, kết hợp nhuần nhuyễn giữa lý thuyết cơ học chất lưu vỉa với bài toán mô phỏng thực tế tại bể Nam Côn Sơn, phục vụ công tác giảng dạy và nghiên cứu chuyên sâu về mô phỏng mỏ.

Nhóm nhà thầu điều hành dầu khí liên doanh và đơn vị EPCI: Giúp tham khảo bài học kinh nghiệm trong thiết kế kiến trúc hệ thống khai thác liên giàn WHd - PUQC kết hợp tàu chứa nổi FSO 300.000 thùng trong điều kiện thời tiết khắc nghiệt trên thềm lục địa Việt Nam.

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

Phương pháp phân tích đường cong suy giảm Arps có đảm bảo độ chính xác cho mỏ khí - condensate không? Hoàn toàn đảm bảo độ tin cậy khi áp dụng cho các mỏ có chu kỳ khai thác ổn định. Luận văn đã chứng minh phương pháp Arps phù hợp vượt trội khi xử lý dữ liệu khai thác hơn 50 tháng của 6 giếng tại mỏ Rồng Đôi, phản ánh chính xác quy luật suy giảm áp suất từ 5.426 psia xuống vùng cận đáy mà không đòi hỏi chi phí mô phỏng 3D quá tốn kém.

Sự khác biệt cốt lõi giữa việc tính toán trên Excel và phần mềm chuyên dụng MBAL là gì? Bảng tính Excel cho phép người kỹ sư hiểu sâu sắc bản chất toán học thông qua các biểu đồ tọa độ bán logarit và log-log giữa CLR và sản lượng cộng dồn. Trong khi đó, phần mềm MBAL tích hợp các phương trình cân bằng vật chất phức tạp, thông số PVT và áp suất phân dòng 100 bar, giúp tự động hóa khớp lịch sử với độ chính xác cao hơn 92%.

Hệ số suy giảm b mang ý nghĩa quyết định như thế nào trong bài toán dự báo? Hệ số b phản ánh cơ chế năng lượng chi phối vỉa chứa. Tại mỏ Rồng Đôi, giá trị b chuyển từ 0 trong giai đoạn đầu sang khoảng 0,3 đến 0,5 khi cơ chế giãn nở khí và phân dị trọng lực phát huy tác dụng, trực tiếp quyết định khả năng duy trì lưu lượng 130 MMscfd khí và hơn 4.160 thùng condensate mỗi ngày.

Tại sao cần lắp đặt thêm máy nén khí phụ trợ trên giàn PUQC sau 14 năm khai thác? Sau 14 năm vận hành liên tục, áp suất vỉa tự nhiên sẽ suy giảm đáng kể xuống dưới ngưỡng tiếp nhận của đường ống Nam Côn Sơn (tối thiểu 80 barg). Việc bổ sung máy nén booster là giải pháp bắt buộc để nâng áp suất khí khô lên 160 barg, bảo đảm khả năng vận chuyển dòng khí về bờ an toàn.

Đặc điểm địa chất nào tại bể Nam Côn Sơn gây khó khăn lớn nhất cho công tác dự báo mỏ? Thách thức lớn nhất là tính không đồng nhất của trầm tích điệp Dừa với 15 vỉa cát mỏng xen kẹp sét có bề dày chỉ từ 4 đến 15 mét. Độ thấm biến thiên cực lớn từ 1,3 đến 250 md khiến dòng chảy chất lưu giữa các vỉa có sự chênh lệch mạnh, dễ dẫn đến hiện tượng ngập nước cục bộ nếu không kiểm soát tốt.

Kết luận

Luận văn thạc sĩ đã giải quyết trọn vẹn bài toán dự báo sản lượng khai thác cho mỏ Rồng Đôi - Rồng Đôi Tây thông qua 5 kết quả nổi bật:

  • Xây dựng thành công quy trình dự báo sản lượng tích hợp hiệu quả giữa công cụ Excel và phần mềm chuyên dụng MBAL cho toàn bộ 6 giếng khai thác chính.
  • Lượng hóa chính xác các tham số động học vỉa với hệ số suy giảm b dao động từ 0 đến 0,5 và tốc độ suy giảm Di từ 15,2% đến 28,4%/năm cho từng vỉa sản phẩm.
  • Thiết lập cơ sở kỹ thuật vững chắc bảo đảm cam kết cung cấp khí thương mại 130 MMscfd và tối ưu hóa công suất thiết kế xử lý 175 MMscfd.
  • Đóng góp luận cứ khoa học quan trọng cho kế hoạch can thiệp giếng và lộ trình đầu tư hệ thống nén phụ trợ sau 14 năm khai thác tại giàn PUQC.
  • Cung cấp mô hình tham chiếu thực tiễn chuẩn mực cho công tác quản lý và khai thác các vỉa khí - condensate phức tạp tại thềm lục địa Việt Nam.

Trong các bước tiếp theo từ năm 2012, các kỹ sư và nhà quản lý mỏ cần nhanh chóng cập nhật mô hình dự báo với dữ liệu đo log khai thác mới, sẵn sàng triển khai các giải pháp công nghệ nâng cao hệ số thu hồi mỏ bền vững.