Tổng quan nghiên cứu

Bồn trũng Cửu Long là bể trầm tích dầu khí trọng điểm, đóng góp hơn 75% tổng sản lượng khai thác hydrocarbon tại thềm lục địa Việt Nam. Trong cấu trúc địa chất khu vực, tầng chứa cát kết Oligocene hạ tại mỏ X nằm ở chiều sâu từ 3.200 mét đến hơn 4.000 mét, sở hữu tiềm năng dầu khí to lớn nhưng có cấu trúc địa chất và môi trường lắng đọng vô cùng phức tạp. Trước đây, các nghiên cứu tập trung chủ yếu vào đối tượng đá móng nứt nẻ và cát kết Miocene, khiến tầng chứa Oligocene hạ chưa được xây dựng một mô hình địa chất ba chiều tĩnh hoàn chỉnh, dẫn đến độ bất định cao khi ước lượng trữ lượng tại chỗ ban đầu cũng như thiết kế kế hoạch khai thác.

Mục tiêu cốt lõi của nghiên cứu là thiết lập quy trình chuẩn hóa xây dựng mô hình địa chất cho đá trầm tích vụn, đồng thời khảo sát, so sánh và lựa chọn thuật toán mô hình hóa tướng đá tối ưu nhất cho tầng chứa Oligocene hạ mỏ X. Phạm vi nghiên cứu được thực hiện trên toàn bộ diện tích cấu trúc mỏ X thuộc Bồn trũng Cửu Long, hoàn thành và nghiệm thu vào tháng 12 năm 2011 dựa trên nguồn dữ liệu địa chất, địa vật lý giếng khoan và địa chấn ba chiều hiện có.

Nghiên cứu mang ý nghĩa thực tiễn sâu sắc khi nâng cao độ chính xác trong việc dự báo phân bố không gian vỉa chứa lên trên 85%, giảm thiểu biên độ sai số bất định khi tính toán thể tích dầu ban đầu tại chỗ xuống dưới 15%. Kết quả này đóng vai trò quyết định giúp các kỹ sư công nghệ mỏ tối ưu hóa vị trí các giếng khoan phát triển, nâng cao hệ số thu hồi dầu khí và tiết kiệm hàng triệu USD chi phí rủi ro khoan thẩm lượng.

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 xây dựng trên nền tảng lý thuyết mô hình hóa địa chất tĩnh ba chiều kết hợp trầm tích học vụn và lý thuyết địa tầng phân tập. Khung lý thuyết tích hợp ba trường phái phương pháp luận chính: lý thuyết mô hình xác định dựa trên các thuật toán nội suy hình học, lý thuyết địa thống kê ngẫu nhiên mô phỏng có điều kiện, và lý thuyết cân bằng mao dẫn vỉa chứa thông qua hàm đòn bẩy Leverett J.

Mô hình nghiên cứu vận hành theo chuỗi liên kết tĩnh sang động, chuyển hóa các thuộc tính cấu trúc, tướng thạch học và thông số vật lý vỉa từ không gian ô mạng mịn sang ô mạng thô cho mô phỏng dòng chảy. Năm khái niệm then chốt được phát triển xuyên suốt bao gồm: tướng trầm tích hệ sông, hàm tương quan không gian thực nghiệm Variogram, mô phỏng tuần tự chỉ định, mô phỏng thống kê đa điểm kết hợp tập ảnh đào tạo, và kỹ thuật thô hóa dữ liệu ô lưới. Về mặt số liệu vật lý, các thuộc tính địa chấn trở kháng âm học dưới 9.000 m/s nhân g/cm3 thể hiện xác suất xuất hiện cát kết lòng sông trên 75%, trong khi độ rỗng hiệu dụng của vỉa dao động từ 2% đến 18% tùy thuộc vào từng loại tướng đá cụ thể.

+-------------------------------------------------------------------------+
|                  KHUNG MÔ HÌNH HÓA ĐỊA CHẤT TÍCH HỢP                     |
+-------------------------------------------------------------------------+

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm toàn diện từ 7 giếng khoan thăm dò và khai thác ký hiệu từ HK1 đến HK5 cùng các tuyến liên kết địa chất ngang qua cấu trúc mỏ. Cỡ mẫu phân tích bao gồm hàng nghìn điểm đo địa vật lý giếng khoan thu thập liên tục ở bước phân giải 0,1524 mét đối với các đường cong phổ gamma tia xạ, mật độ thể tích, độ rỗng nơ-tron, sóng âm và điện trở suất sâu. Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu tầng phân lớp có hệ thống qua toàn bộ chiều dày tầng cát kết Oligocene hạ tập A, đối chiếu trực tiếp với các lát cắt mẫu lõi thành giếng và mẫu liền khối nhằm chuẩn hóa các phân loại thạch học.

Phương pháp phân tích số liệu được thực hiện trên nền tảng phần mềm chuyên dụng Irap RMS của hãng Roxar, ứng dụng đa dạng các thuật toán từ nội suy Indicator Kriging đến mô phỏng địa thống kê ngẫu nhiên như Sequential Indicator Simulation, Truncated Gaussian Simulation, Multipoint Statistics và mô phỏng dựa trên đối tượng hình học General Marked Point Process. Lý do lựa chọn tổ hợp phương pháp này là vì môi trường lắng đọng của mỏ X mang đặc trưng châu thổ sông uốn khúc có tính bất đồng nhất cao; các phương pháp nội suy truyền thống luôn làm mờ ranh giới cát kết, trong khi các thuật toán địa thống kê ngẫu nhiên hiện đại cho phép tái hiện chính xác hình học các thân cát lòng sông và các bãi cát vỡ đê. Toàn bộ tiến trình nghiên cứu và xử lý dữ liệu được tiến hành liên tục trong giai đoạn từ năm 2010 đến cuối năm 2011.

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

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

Quá trình phân tích thạch học và môi trường lắng đọng tầng sản phẩm tập A mỏ X xác định hệ thống tướng trầm tích bao gồm ba thành phần chủ đạo: cát kết lòng sông chiếm tỷ phần trung bình 38%, bãi cát vỡ bờ chiếm 18%, và sét bùn đồng bằng ngập lụt chiếm 44% tổng thể tích vỉa.

PHÂN BỐ TỶ PHẦN CÁC LOẠI TƯỚNG TRẦM TÍCH VỈA A

Về mặt so sánh thuật toán, phương pháp Object-based và phương pháp thống kê đa điểm MPS thể hiện sự vượt trội hoàn toàn trong việc tái lập cấu trúc hình học các thân cát lạch sông uốn khúc với bề rộng biến thiên từ 150 mét đến 320 mét và tỷ lệ bề rộng trên bề dày dao động trong khoảng 25 đến 40 lần. Ngược lại, thuật toán Sequential Indicator Simulation chỉ đạt mức độ khớp nối cấu trúc hình học dưới 58%, tạo ra nhiều hiệu ứng nhiễu hạt rời rạc không phù hợp với thực tế địa chất.

Mô hình hóa độ rỗng và độ thấm trên nền tảng tướng đá đã cải thiện vượt bậc hệ số tương quan rỗng - thấm. Khi không phân tách theo tướng, hệ số xác định tương quan chỉ đạt mức 0,52, nhưng khi được khống chế bởi từng loại tướng đá riêng biệt, hệ số tương quan tăng lên đạt mức 0,84 đối với tướng cát lòng sông. Độ rỗng hiệu dụng trung bình trong thân cát đạt 12% đến 16%, trong khi độ thấm biến thiên mạnh mẽ từ 5 mD đến trên 1.200 mD theo phân bố logarit chuẩn.

ĐẶC TRƯNG THỐNG KÊ CỦA CÁC THUẬT TOÁN MÔ HÌNH HÓA TƯỚNG

Kết quả tính toán thể tích rỗng chứa dầu và khí Hydrocarbon Pore Volume cho thấy việc áp dụng mô hình tướng ngẫu nhiên giúp thu hẹp khoảng dao động bất định xuống còn 14%, so với biên độ sai lệch lên đến 35% khi sử dụng các phương pháp tính toán thể tích truyền thống.

Biểu đồ tương quan giữa Độ rỗng và Logarit Độ thấm theo từng tướng:
Log(K) (mD)
  ^
4 |                                 *  *  (Cát lòng sông: R2 = 0.84)
  |                              *  *  *
3 |                           *  *  *
  |                     +  +  +  +        (Cát vỡ bờ: R2 = 0.71)
2 |                  +  +  +
  |            #  #  #                    (Sét ngập lụt: K < 1 mD)
1 |      #  #  #
  +---------------------------------------------> Độ rỗng Phi (%)
  0      4      8      12     16     20

Thảo luận kết quả

Sự khác biệt rõ rệt giữa các phương pháp xuất phát từ khả năng xử lý thông tin không gian phi tuyến tính. Phương pháp xác định Deterministic làm phẳng các dị thường địa chất, dẫn đến đánh giá quá mức tính liên thông tại các đới biên vỉa. Ngược lại, thuật toán MPS tận dụng hiệu quả tập ảnh đào tạo ba chiều và đường cong tỷ phần tướng dọc để ràng buộc không gian, giúp bảo toàn tính dị hướng trầm tích theo phương dòng chảy chính từ Tây sang Đông.

Khi đối chiếu với các công bố khoa học quốc tế về trầm tích sông uốn khúc hệ tầng Mina El Carmen thuộc bồn trũng San Jorge tại Argentina, kết luận của luận văn hoàn toàn tương đồng: phương pháp thống kê đa điểm MPS và mô phỏng đối tượng hình học Object-based là sự lựa chọn duy nhất có khả năng mô tả chính xác các lạch sông có độ uốn khúc cao và bề rộng hẹp dưới 300 mét.

Về mặt hiển thị trực quan, dữ liệu nghiên cứu được biểu diễn tối ưu thông qua hệ thống bản đồ xác suất xuất hiện thân cát hai chiều, biểu đồ tỷ phần tướng dọc thể hiện quy luật trầm tích mịn dần lên trên, và các bảng đối sánh chỉ tiêu vật lý vỉa theo từng kịch bản mô phỏng đồng xác suất. Điều này giúp các nhà địa chất nhận diện rõ ranh giới các đới nước - dầu chuyển tiếp được tính toán thông qua hàm Leverett J với áp suất mao dẫn vỉa chứa biến thiên theo chiều cao cột dầu từ 10 psi đến trên 150 psi.

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

  1. Chuẩn hóa áp dụng thuật toán MPS và Object-based trong xây dựng mô hình vỉa: Tổng công ty Thăm dò Khai thác Dầu khí và các công ty liên doanh điều hành mỏ tại Bồn trũng Cửu Long cần ban hành quy định bắt buộc sử dụng phương pháp thống kê đa điểm MPS hoặc Object-based khi mô hình hóa các tầng cát kết Oligocene hạ trước quý 2 năm 2012, hướng tới mục tiêu nâng cao độ tin cậy mô phỏng vỉa lên trên 90%.
  2. Gia tăng mật độ lấy mẫu lõi và đo ảnh vi điện trở giếng khoan: Các đơn vị khoan và địa vật lý giếng khoan cần thực hiện kế hoạch lấy tối thiểu 30 mét đến 50 mét mẫu lõi nguyên khối kết hợp chụp ảnh FMI tại mỗi giếng khoan thẩm lượng mới trong giai đoạn 2012 đến 2015, nhằm cung cấp dữ liệu kiểm soát góc nghiêng phân lớp và hướng dòng chảy cho tập ảnh đào tạo.
  3. Hoàn thiện quy trình thô hóa ô lưới phục vụ mô phỏng động: Nhóm chuyên trách công nghệ mỏ phối hợp với kỹ sư mô hình địa chất triển khai chuẩn hóa thuật toán trung bình hóa có trọng số để rút giảm số lượng ô lưới từ 3 triệu ô xuống khoảng 250.000 đến 400.000 ô trong vòng 6 tháng trước khi nạp vào phần mềm mô phỏng khai thác, đảm bảo mức độ sai lệch thất thoát thuộc tính vỉa không vượt quá 5%.
  4. Tích hợp khối thuộc tính đảo địa chấn ba chiều vào bản đồ xu hướng tướng: Bộ phận xử lý địa chấn và địa vật lý cần hoàn tất việc chuyển đổi các khối biên độ địa chấn thành khối trở kháng âm học AI và áp dụng bản đồ xác suất thạch học 3D trong năm 2012, nhằm khống chế chặt chẽ tính phân dị không gian của các thân cát nằm giữa các giếng khoan.

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

  • Kỹ sư mô hình địa chất tại các công ty dầu khí: Nhận được quy trình chuẩn hóa chi tiết từng bước xây dựng mô hình tĩnh 3D trên phần mềm Irap RMS hoặc Petrel, giúp giải quyết các bài toán phân bố thân cát phức tạp cho các mỏ dầu khí đang trong giai đoạn phát triển tại Việt Nam.
  • Kỹ sư công nghệ mỏ và mô phỏng vỉa: Nắm vững phương pháp thô hóa ô lưới và cách thiết lập hàm áp suất mao dẫn Leverett J theo từng loại tướng thạch học, phục vụ trực tiếp cho công tác cân bằng lịch sử khai thác và dự báo lưu lượng dòng chảy chất lưu chính xác hơn.
  • Các nhà quản lý dự án và chuyên gia đánh giá trữ lượng: Sử dụng tài liệu làm khung tham chiếu đánh giá rủi ro địa chất, thẩm định độ tin cậy của các báo cáo phân cấp trữ lượng dầu khí tại chỗ theo tiêu chuẩn quốc tế của Hiệp hội Kỹ sư Dầu khí SPE cho các cấp P1, P2 và P3.
  • 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í: Khai thác luận văn như một công trình tổng quan học thuật sâu sắc về lý thuyết địa thống kê ứng dụng, cung cấp các ví dụ thực chứng sinh động phục vụ công tác giảng dạy và nghiên cứu mở rộng trong lĩnh vực khoa học trái đất.

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

Mô hình địa chất ngẫu nhiên Stochastic vượt trội hơn mô hình xác định Deterministic ở những điểm nào?
Mô hình xác định chỉ cho ra một kết quả duy nhất có bề mặt trơn nhẵn, làm triệt tiêu tính bất đồng nhất và phóng đại tính liên thông của vỉa chứa. Ngược lại, mô hình ngẫu nhiên tạo ra hàng trăm kịch bản đồng xác suất có khả năng tích hợp đầy đủ dữ liệu địa chấn và giếng khoan, giúp định lượng chính xác các yếu tố bất định và giảm sai số tính toán thể tích dầu ban đầu tại chỗ xuống dưới 15%.

Tại sao thuật toán thống kê đa điểm MPS được đánh giá là tối ưu nhất cho tầng Oligocene hạ mỏ X?
Tầng cát kết Oligocene hạ mỏ X có môi trường trầm tích châu thổ sông với các lạch uốn khúc hẹp có bề rộng dưới 300 mét. Thuật toán MPS khắc phục hoàn toàn nhược điểm của biến số không gian hai điểm cổ điển bằng cách quét mẫu cấu trúc từ tập ảnh đào tạo ba chiều, giúp tái lập hoàn hảo độ uốn khúc phức tạp và bảo toàn tính liên tục của thân cát.

Hàm Leverett J đóng vai trò gì trong việc xây dựng mô hình phân bố bão hòa nước?
Hàm Leverett J chuyển đổi toàn bộ dữ liệu áp suất mao dẫn từ mẫu lõi thành một đường cong phổ quát đại diện cho từng loại đá chứa. Thông qua việc tích hợp độ rỗng, độ thấm và sức căng bề mặt theo chiều cao cột chất lưu, hàm J cho phép ngoại suy chính xác sự biến đổi của độ bão hòa nước từ 100% tại mực nước tự do về độ bão hòa nước dư trong vỉa dầu.

Quy trình thô hóa ô lưới Upscaling có gây mất mát thông tin địa chất quan trọng hay không?
Quá trình thô hóa chuyển đổi hàng triệu ô mạng tĩnh sang hàng trăm nghìn ô mạng động chắc chắn làm giảm độ chi tiết. Tuy nhiên, nếu áp dụng phương pháp trung bình điều hòa cho phương vuông góc lớp và trung bình số học có trọng số theo tướng cho phương song song lớp, mức độ thất thoát các đặc trưng bất đồng nhất của vỉa sẽ được kiểm soát ở mức an toàn dưới 5%.

Làm thế nào để kiểm soát chất lượng mô hình tướng trước khi mô hình hóa các thông số vật lý vỉa?
Người xây dựng mô hình cần thực hiện kiểm tra chéo thống kê giữa dữ liệu giếng khoan ban đầu và giá trị gán trên ô lưới thông qua biểu đồ phân bố tỷ phần tướng dọc và biểu đồ tần suất. Tỷ phần thể tích của từng loại tướng đá trên toàn bộ mô hình ba chiều bắt buộc phải khớp với số liệu thống kê giếng khoan với sai số cho phép không vượt quá 2% đến 3%.

Kết luận

  • Luận văn đã thiết lập thành công quy trình 9 bước chuẩn hóa trong xây dựng mô hình địa chất tĩnh ba chiều cho đối tượng đá trầm tích vụn tại thềm lục địa Việt Nam.
  • Khẳng định môi trường lắng đọng tầng sản phẩm tập A mỏ X là hệ thống đồng bằng châu thổ sông với tỷ phần cát kết chứa dầu chiếm 38% và cát vỡ bờ chiếm 18%.
  • Chứng minh tính ưu việt vượt trội của thuật toán địa thống kê MPS và Object-based trong việc mô tả hình học lạch sông uốn khúc, vượt trội hơn hẳn so với thuật toán SIS và phương pháp nội suy truyền thống.
  • Xây dựng thành công chuỗi mô hình phân bố thuộc tính vật lý vỉa bao gồm độ rỗng, độ thấm và độ bão hòa nước theo hàm Leverett J, giúp kiểm soát chặt chẽ tính bất đồng nhất của vỉa chứa.
  • Giảm thiểu biên độ bất định trong tính toán trữ lượng dầu khí tại chỗ xuống dưới 14%, cung cấp cơ sở khoa học vững chắc cho công tác lập kế hoạch khai thác tối ưu.

Đóng góp lớn nhất của công trình là cung cấp phương pháp luận thực chứng và giải pháp kỹ thuật tối ưu trong việc lựa chọn thuật toán mô hình hóa tướng đá cho các bể trầm tích phức tạp. Các bước tiếp theo cần triển khai bao gồm việc hoàn thiện tập ảnh đào tạo độ phân giải cao và tích hợp kết quả mô hình tĩnh vào phần mềm mô phỏng dòng chảy động trong giai đoạn 2012 đến 2013. Các đơn vị khai thác dầu khí và các nhóm nghiên cứu địa chất vỉa hãy áp dụng ngay quy trình mô hình hóa tướng ngẫu nhiên này vào các đề án phát triển mỏ thực tế để tối ưu hóa chi phí và nâng cao hiệu quả khai thác tài nguyên năng lượng quốc gia.