Tổng quan luận án

Trong giai đoạn đầu phát triển các mỏ dầu khí cận biên (marginal fields) ngoài khơi, việc đưa ra quyết định đầu tư và xây dựng kế hoạch phát triển mỏ (Field Development Plan) luôn đối mặt với mức độ rủi ro cao. Mỏ khí ngưng tụ (gas condensate) ST-X thuộc Lô 15.10, bồn trũng Cửu Long là đối tượng có đặc điểm địa chất phức tạp, cấu trúc bị chia cắt bởi nhiều hệ thống đứt gãy, vỉa có độ thấm trung bình thấp, tính bất đồng nhất cao, tồn tại nhiều biên chắn không thấm và phân bố ít nhất 2 vùng đặc tính lưu chất (PVT) khác biệt. Hơn nữa, chi phí khoan và vận hành ngoài khơi rất tốn kém.

Phương pháp mô phỏng vỉa truyền thống khi ứng dụng cho các đối tượng này bộc lộ hạn chế lớn: việc khảo sát các yếu tố không chắc chắn đòi hỏi số lượng kịch bản chạy mô phỏng thủy động lực học rất lớn, tiêu tốn thời gian tính toán và tài nguyên máy tính. Mặt khác, cách tiếp cận truyền thống thường chỉ tìm kiếm một mô hình khớp lịch sử đơn lẻ ("best match"), không phản ánh được đầy đủ khoảng tin cậy của trữ lượng thu hồi và sản lượng khai thác trong tương lai. Khoảng trống nghiên cứu đặt ra là cần xây dựng một quy trình định lượng tổng thể, có khả năng sàng lọc nhanh các thông số nhạy cảm, giảm thiểu số lần chạy mô phỏng số nhưng vẫn xác định chính xác các khoảng xác suất phân bố sản lượng khai thác.

Luận án tiến sĩ kỹ thuật của tác giả Vũ Việt Hưng (2019) với đề tài "Đánh giá các yếu tố không chắc chắn trong giai đoạn đầu phát triển mỏ khí ngưng tụ ST - X (Lô 15.10, ngoài khơi Việt Nam)", chuyên ngành Kỹ thuật dầu khí (Mã số: 62.04), thực hiện dưới sự hướng dẫn của TS. Mai Cao Lân và TS. Nguyễn Chu Chuyên tại Trường Đại học Bách khoa – Đại học Quốc gia TP. Hồ Chí Minh, được triển khai nhằm giải quyết trực tiếp vấn đề này.

Mục tiêu nghiên cứu:

  1. Xác định các yếu tố không chắc chắn và đánh giá định lượng ảnh hưởng của chúng lên kết quả dự báo khai thác, đồng thời đề xuất phương án phát triển mỏ tối ưu cho mỏ ST-X với mục tiêu cực đại hệ số thu hồi và cực tiểu rủi ro.
  2. Đưa ra quy trình hoàn thiện đánh giá các yếu tố không chắc chắn trong công tác phát triển mỏ dầu khí.

Nhiệm vụ nghiên cứu:

  1. Xác định các thông số không chắc chắn và khoảng biến thiên của chúng từ kết quả phân tích thử vỉa (DST - Drill Stem Test).
  2. Khảo sát và phân tích các phương pháp thiết kế thực nghiệm (Experimental Design) ứng dụng trong mô hình mô phỏng vỉa nhằm sàng lọc, đánh giá mức độ ảnh hưởng của các thông số không chắc chắn đến kết quả dự báo.
  3. Khảo sát và xây dựng phương pháp bề mặt phản hồi (Response Surface Modeling - RSM); kết hợp phương pháp mô phỏng Monte Carlo trên bề mặt phản hồi nhằm giảm thiểu thời gian và số lần chạy mô phỏng cho bài toán dự báo khai thác.
  4. Xây dựng mô hình dự báo và tính toán khoảng tin cậy trữ lượng thu hồi ở các cấp P10 – P50 – P90 cho kết quả sản lượng khai thác, từ đó đề xuất phương án phát triển tối ưu mỏ khí ngưng tụ ST-X.

Luận điểm bảo vệ:

  • Luận điểm 1: Hệ phương pháp thiết kế thực nghiệm và phương pháp bề mặt phản hồi ứng dụng trong mô hình mô phỏng vỉa cho phép sàng lọc, đánh giá sự ảnh hưởng của các thông số không chắc chắn đến kết quả dự báo khai thác mỏ khí ngưng tụ ST-X.
  • Luận điểm 2: Phương pháp mô phỏng Monte Carlo trên bề mặt phản hồi cho phép rút ngắn thời gian chạy mô phỏng, từ đó tính toán khoảng tin cậy trữ lượng thu hồi các cấp P10, P50, P90 cho kết quả dự báo sản lượng khai thác trong tương lai mỏ khí ngưng tụ ST-X.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Mỏ khí ngưng tụ ST-X nằm trong Lô 15.10 thuộc bồn trũng Cửu Long, ngoài khơi thềm lục địa phía Nam Việt Nam do Công ty Liên doanh Điều hành chung Cửu Long (Cửu Long JOC) điều hành.
  • Phạm vi nghiên cứu: Xác định và đánh giá định lượng ảnh hưởng của các yếu tố không chắc chắn từ dữ liệu động (thử vỉa DST, MDT, PVT) và dữ liệu tĩnh (địa chất, địa chấn, mẫu lõi) lên phương án phát triển mỏ, trọng tâm là tầng cát kết Oligoxen (tập E và tập F).

Tổng quan tài liệu và vị trí của luận án

Luận án đã tổng hợp, phân tích các công trình nghiên cứu trong và ngoài nước liên quan đến quản lý rủi ro và đánh giá yếu tố không chắc chắn trong công nghệ mỏ:

Các nghiên cứu quốc tế:

  • P. Foley; Behrenbruch và cộng sự (1989): Nghiên cứu rủi ro trong thăm dò và phát triển mỏ bằng phương pháp đường cong ước lượng (Expectation Curve Method). Đây là phương pháp xác suất đơn giản, chưa tích hợp được mô hình mô phỏng vỉa phức tạp.
  • T. Charles và cộng sự (2001): Thực nghiệm định lượng ảnh hưởng của các yếu tố không chắc chắn địa chất và công nghệ mỏ trên 4 mỏ riêng biệt. Hạn chế là chưa hình thành quy trình tổng quát tích hợp từ số liệu tĩnh đến số liệu động.
  • S. Deutsch (2004): Ứng dụng phương pháp địa thống kê và xác suất Bayesian kết hợp tài liệu địa chấn với kết quả giếng khoan để lựa chọn vị trí giếng, nhưng thiếu quy trình liên tục xuyên suốt từ thăm dò đến phát triển.
  • SPE 93280 (2005): Áp dụng thuật toán Neighbourhood Algorithm (NA) kết hợp lấy mẫu Markov Chain Monte Carlo (MCMC) trên mô hình black oil cho mỏ Rigel để cung cấp tập hợp các mô hình hiệu chỉnh tối ưu.
  • SPE 106039 (2007): Ứng dụng thiết kế thực nghiệm Latin Hypercube và Bayesian tuyến tính để hiệu chỉnh mô hình lịch sử cho vỉa chứa gồm 120 giếng, lịch sử khai thác 58 năm, chọn ra 11 mô hình tốt nhất để dự báo P10, P50, P90.
  • Eric Bickel và Reidar B. Bratvold (2008): Khảo sát 494 chuyên gia ngành dầu khí, chỉ ra rằng việc định lượng không chắc chắn cần phải gắn liền trực tiếp với khung công cụ hỗ trợ ra quyết định phát triển mỏ.
  • Mathieu Denis Feraille và Daniel Busby (2009): Xây dựng quy trình quản lý yếu tố không chắc chắn cho mỏ cận biên bằng mô hình bề mặt đáp ứng phi tham số, thiết kế thích ứng và thuật toán MCMC.
  • Ibiada Harrison Itotoi và cộng sự (SPE-140619, 2010): Quản lý yếu tố không chắc chắn trong phát triển mỏ khí bằng thiết kế thực nghiệm, tuy nhiên nghiên cứu chỉ tập trung vào phương án phát triển bề mặt mà chưa đánh giá trong quá trình hiệu chỉnh mô hình.
  • Michael Lev Litvak và cộng sự (SPE-146512, 2011): Tối ưu hóa phát triển mỏ dưới ảnh hưởng của yếu tố không chắc chắn dưới lòng đất, sử dụng kỹ thuật thống kê chuỗi chia nhỏ để giảm thiểu thời gian CPU.
  • Babafemi Anthony Ogunyomi, Christopher J. Jablonowski, Larry Wayne Lake (2011): Xây dựng mô hình tài sản tích hợp (Integrated Asset Model - IAM) phục vụ phân tích bất định và ra quyết định khai thác.
  • Samson Imoh Essien (2012): Hiệu chỉnh mô hình và dự báo khai thác cho tập E mỏ Norne theo phương pháp truyền thống, bộc lộ điểm yếu tốn nhiều thời gian và chỉ cho ra một mô hình đại diện duy nhất.
  • R. Salinas và cộng sự (SPE-169468, 2014): Nhận diện và đánh giá yếu tố không chắc chắn giai đoạn phát triển sớm thông qua phân tích độ nhạy với biểu đồ Tornado và biểu đồ Pareto.

Các nghiên cứu trong nước:

  • Nguyễn Phúc Khải (2011): Đánh giá định lượng các yếu tố ảnh hưởng đến quá trình hiệu chỉnh mô hình vỉa tầng Mioxen dưới mỏ Sóc Đen bằng phương pháp bề mặt đáp ứng kết hợp thiết kế thực nghiệm thừa số 2 cấp, 3 cấp và mô phỏng Monte Carlo.
  • Nguyễn Ngọc Tuấn Anh (2011): Hiệu chỉnh mô hình thủy động lực có xét đến thông số không chắc chắn cho mỏ Sư Tử Đen bằng Latin Hypercube và mô hình hồi quy tuyến tính đa chiều.
  • Nguyễn Phước Long (2016): Dự báo khai thác có xét đến các thông số không chắc chắn cho tầng B mỏ Sư Tử Đen.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Các nghiên cứu trước đây thường thực hiện đơn lẻ ở từng giai đoạn riêng biệt (hoặc chỉ hiệu chỉnh mô hình, hoặc chỉ tối ưu hóa thiết bị bề mặt/vị trí giếng) hoặc sử dụng các phương pháp thống kê truyền thống chưa liên kết chặt chẽ từ khâu giải đoán thử vỉa đến dự báo xác suất. Luận án giải quyết khoảng trống này bằng cách xây dựng một quy trình hoàn chỉnh khép kín: trích xuất miền biến thiên tham số từ tài liệu thử vỉa (DST) $\rightarrow$ sàng lọc các yếu tố nhạy cảm chính bằng thiết kế thực nghiệm (Plackett-Burman / Latin Hypercube) $\rightarrow$ xây dựng phương trình bề mặt phản hồi (RSM) $\rightarrow$ thực hiện 10.000 lần mô phỏng Monte Carlo trực tiếp trên bề mặt phản hồi để tính toán chính xác dải trữ lượng/sản lượng P10-P50-P90 mà không cần chạy mô phỏng thủy động lực học lặp lại nhiều lần.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và các khái niệm sử dụng

  • Yếu tố rủi ro và không chắc chắn trong dầu khí: Luận án phân định rõ yếu tố rủi ro (sự kiện tiềm ẩn gắn liền với tần suất xuất hiện và mức độ tổn thất) và yếu tố không chắc chắn (sự hiểu biết chưa đầy đủ về các tham số địa chất và công nghệ mỏ do giới hạn dữ liệu thu thập). Các yếu tố tĩnh gồm cấu trúc địa chất, ranh giới chất lưu (GOC, OWC), thể tích vỉa; các yếu tố động gồm độ thấm, tính bất đồng nhất, ranh giới đứt gãy, cơ chế năng lượng vỉa và ứng xử PVT.
  • Đặc tính nhiệt động học vỉa khí ngưng tụ: Quá trình suy giảm áp suất vỉa xuống dưới áp suất điểm sương (dew point pressure) làm chất lỏng ngưng tụ (condensate drop-out / condensate banking) lắng đọng xung quanh đáy giếng, làm suy giảm nghiêm trọng độ thấm tương đối của pha khí và hệ số thu hồi chất lưu.

Phương pháp nghiên cứu thực tế

Luận án áp dụng kết hợp các phương pháp toán học thống kê, phân tích dữ liệu mỏ và mô phỏng số:

Phương pháp Kỹ thuật / Công cụ cụ thể Mục đích áp dụng trong luận án
Phân tích thử vỉa (DST) Đồ thị Log-Log, Bán Log (mô hình dòng chảy tỏa tia, đơn tầng, đa tầng, biên chắn song song) Nhận diện tính chất vỉa, độ thấm ($k$), khoảng cách biên đứt gãy, hiện tượng vây condensat và xác định khoảng biến thiên tham số đầu vào
Thiết kế thực nghiệm (ED) Ma trận Hadamard, thiết kế Plackett-Burman bậc 2 (cho 11 thông số), Lấy mẫu Latin Hypercube (LHS), Full Factorial Design Giảm số lần chạy mô phỏng số, sàng lọc các yếu tố có mức độ ảnh hưởng lớn nhất đến hàm mục tiêu
Phân tích độ nhạy Biểu đồ Tornado, Biểu đồ Pareto, Hệ số tương quan Pearson Định lượng mức độ tác động của từng thông số không chắc chắn lên kết quả mô phỏng
Mô hình bề mặt phản hồi (RSM) Hồi quy đa thức bậc hai (Response Surface Methodology), kiểm định đồ thị tương quan chéo (crossplot) Xây dựng phương trình toán học thay thế (proxy model) cho mô hình mô phỏng thủy động lực
Mô phỏng Monte Carlo Chạy mô phỏng ngẫu nhiên 10.000 lần lặp trên phương trình RSM, xác định hàm PDF và hàm CDF Tính toán nhanh phân bố xác suất và khoảng tin cậy của sản lượng khai thác tại các cấp P10, P50, P90
Mô phỏng số vỉa chứa 3D Phần mềm mô phỏng Eclipse, MEPO (Schlumberger) Xây dựng mô hình cấu trúc, mô hình tham số tĩnh (độ rỗng, độ thấm, bão hòa), mô hình động (áp suất mao dẫn, độ thấm tương đối) và hiệu chỉnh lịch sử khai thác

Nguồn tư liệu và dữ liệu

  • Cơ sở dữ liệu địa chất, đo đạc địa vật lý giếng khoan, tài liệu thu nổ và minh giải địa chấn Lô 15.10.
  • Dữ liệu thử vỉa DST, đo áp suất thành hệ MDT, phân tích mẫu lõi đặc biệt, số liệu phân tích PVT từ các giếng khoan thăm dò và thẩm lượng ST-A, ST-B, ST-C, ST-D.
  • Các báo cáo đánh giá trữ lượng, báo cáo kế hoạch phát triển mỏ ST-X từ Liên doanh Cửu Long JOC và Tổng công ty Thăm dò Khai thác Dầu khí (PVEP).
  • Dữ liệu khai thác thử dài hạn tại mỏ ST-X để đối sánh, kiểm chứng mô hình.

Nội dung chính theo từng chương

Chương 1: Tổng quan về lĩnh vực và đối tượng nghiên cứu

Chương này trình bày tổng quan tình hình nghiên cứu lý thuyết, vị trí địa lý, bối cảnh địa chất và kết quả hoạt động thăm dò thẩm lượng tại mỏ khí ngưng tụ ST-X:

  • Vị trí địa lý: Mỏ ST-X nằm ở phía Bắc – Đông Bắc bồn trũng Cửu Long thuộc Lô 15.10, cách Vũng Tàu khoảng 135 km về phía Đông, độ sâu mực nước biển trung bình 56 m. Khí hậu thuộc vùng á xích đạo đại dương, chịu ảnh hưởng của hai mùa gió Đông Bắc và Tây Nam.
  • Đặc điểm địa tầng và kiến tạo: Bồn trũng Cửu Long là bể trầm tích rift kiểu tách giãn. Lịch sử địa chất trải qua 5 giai đoạn: (1) Bắt đầu đứt gãy / lắng đọng tập F; (2) Đỉnh nhô cao của đứt gãy / lắng đọng tập E; (3) Thời kỳ lấp bể / lắng đọng tập D; (4) Giai đoạn xoắn bể / lắng đọng tập C; (5) Nén lún vì nhiệt / lắng đọng tập A và B. Địa tầng từ cổ đến trẻ gồm: Đá móng granitoid trước Kainozoi nứt nẻ; Hệ tầng Cà Cối (Eoxen); Hệ tầng Trà Cú (Eoxen – Oligoxen sớm, gồm tập F và tập E dưới); Hệ tầng Trà Tân (Oligoxen muộn, gồm tập D và C - tầng đá mẹ chính); Hệ tầng Bạch Hổ (Mioxen sớm, tập BI, chứa tầng sét Rotalia đóng vai trò tầng chắn khu vực dày 35–150 m); Hệ tầng Côn Sơn (Mioxen giữa, tập BII); Hệ tầng Đồng Nai (Mioxen muộn, tập BIII); Hệ tầng Biển Đông (Plioxen – Đệ tứ, tập A).
  • Hệ thống dầu khí tập E và F:
    • Tập E: Môi trường sông ngòi xen kẽ ngập lụt hồ; cát kết có độ rỗng từ 8–15%, độ thấm biến thiên từ 0,01 đến 100 mD.
    • Tập F: Trầm tích cát kết lắng đọng trực tiếp trên móng xói mòn, thành phần chủ yếu là arkose và arkose lithic; độ rỗng thường nhỏ hơn 10%, độ thấm nhỏ hơn 10 mD.
  • Kết quả thử vỉa các giếng thăm dò, thẩm lượng:
Giếng khoan Tầng thử vỉa Khoảng thử nghiệm (DST) Lưu lượng khí cực đại Lưu lượng condensat / dầu
ST-A Tập D (DST#1) Cát kết Oligoxen 1,7 triệu bộ khối khí/ngày 590 thùng dầu/ngày
Tập F (DST#2) Cát kết Oligoxen 32,0 triệu bộ khối khí/ngày 3.604 thùng condensat/ngày
Tập E (DST#3) Cát kết Oligoxen 37,7 triệu bộ khối khí/ngày 4.033 thùng condensat/ngày
ST-B Nóc móng Đá móng nứt nẻ 10,0 triệu bộ khối khí/ngày 1.200 thùng condensat/ngày
ST-C Tập F (DST#1) Cát kết Oligoxen 8,7 triệu bộ khối khí/ngày 1.086 thùng condensat/ngày
Tập E (DST#2) Cát kết Oligoxen 10,8 triệu bộ khối khí/ngày 1.372 thùng condensat/ngày
Tập E (DST#2A) Cát kết Oligoxen 5,2 triệu bộ khối khí/ngày 865 thùng condensat/ngày
ST-D Tập E Cát kết Oligoxen 7,9 triệu bộ khối khí/ngày 2.200 thùng condensat/ngày ($38^\circ\text{API}$)

Chương 2: Cơ sở lý thuyết và phương pháp nghiên cứu

Chương 2 thiết lập cơ sở toán học và phương pháp luận cho toàn bộ quy trình đánh giá:

  • Phương pháp thiết kế thực nghiệm (ED): Trình bày ma trận Hadamard, phương pháp Plackett-Burman bậc 2 cho phép khảo sát 11 thông số biến thiên với số lượng mô phỏng cực tiểu, và kỹ thuật lấy mẫu Latin Hypercube (LHS) phân bố đều các giá trị ngẫu nhiên trong không gian đa chiều của các biến vỉa.
  • Mô hình bề mặt phản hồi (RSM): Thiết lập phương trình hàm hồi quy đa thức mô tả mối quan hệ giữa các thông số vỉa đầu vào không chắc chắn ($X_i$) và đại lượng đầu ra cần dự báo ($Y$ – tổng sản lượng khí khai thác mỏ FGPT, tổng sản lượng condensat FOPT). Đánh giá chất lượng bề mặt thông qua sai số hồi quy và biểu đồ tương quan chéo (crossplot).
  • Phân tích độ nhạy và tối ưu hóa: Sử dụng hệ số tương quan Pearson, biểu đồ Tornado và biểu đồ Pareto để định lượng, xếp hạng mức độ ảnh hưởng của từng thông số, từ đó loại bỏ các yếu tố thứ yếu.
  • Mô phỏng Monte Carlo trên RSM: Thay vì chạy mô hình mô phỏng thủy động lực hàng ngàn lần trên máy tính, thuật toán Monte Carlo được thực hiện 10.000 lần lặp trực tiếp trên phương trình bề mặt phản hồi đã xây dựng. Kết quả xuất ra hàm mật độ xác suất (PDF) và hàm phân bố tích lũy (CDF) để xác định chính xác các giá trị phân vị P10 (kịch bản cao/lạc quan), P50 (kịch bản trung bình/kỳ vọng), P90 (kịch bản thấp/thận trọng).

Chương 3: Nhận diện các yếu tố không chắc chắn mỏ ST-X từ phân tích thử vỉa (DST)

Chương 3 tập trung minh giải dữ liệu thử vỉa động (DST) từ các giếng khoan ST-A, ST-B, ST-C, ST-D để phát hiện bản chất các yếu tố bất định:

  • Minh giải đường cong áp suất: Sử dụng đồ thị Log-Log và bán Log cho giai đoạn suy giảm và phục hồi áp suất để nhận diện các chế độ dòng chảy: mô hình đơn tầng, mô hình đa tầng, dòng chảy tỏa tia và ảnh hưởng của các biên chắn song song.
  • Xác định các yếu tố không chắc chắn chính từ DST:
    1. Độ thấm ($k$): Dao động mạnh giữa các tầng và các khối vỉa, giá trị thực tế thấp hơn so với đánh giá ban đầu từ tài liệu địa vật lý giếng khoan.
    2. Khoảng cách tới các biên chắn/đứt gãy không thấm: Sự xuất hiện của các đoạn dốc trên đường đạo hàm áp suất (đồ thị Log-Log) chỉ ra các khối vỉa bị phân chia khép kín.
    3. Hiện tượng vây lắng đọng condensat (condensate bank): Áp suất đáy giếng trong quá trình thử vỉa giảm xuống dưới áp suất bọt khí/điểm sương làm nghẽn dòng chảy chất lưu quanh vùng cận đáy giếng.
    4. Đặc tính chất lưu (PVT): Sự biến thiên của tỷ số khí – ngưng tụ (GCR), hệ số thể tích thành hệ khí ($B_g$), và tỷ số bốc hơi ($r_v$) giữa các khối vỉa.
  • Luận án đã tổng hợp và xác lập bảng miền giá trị biến thiên (giá trị nhỏ nhất, trung bình, lớn nhất) cùng hàm phân bố xác suất cho từng thông số nhận diện được làm dữ liệu đầu vào cho bước mô phỏng.

Chương 4: Quy trình đánh giá ảnh hưởng của các yếu tố không chắc chắn và định hướng phát triển mỏ tối ưu

Chương này trình bày việc áp dụng quy trình tích hợp vào mô hình số thực tế mỏ ST-X:

  • Xây dựng mô hình mô phỏng 3D: Thiết lập mô hình địa chất 3D thể hiện cấu trúc phân bố tướng đá, mô hình phân bố độ rỗng ($\Phi$), độ thấm ($k$) và độ bão hòa nước ban đầu ($S_{wi}$). Đưa vào các đường cong áp suất mao dẫn ($P_c$) và độ thấm tương đối hai pha/ba pha ($k_{ro}, k_{rw}, k_{rg}$).
  • Sàng lọc thông số và rút gọn khoảng biến thiên:
    • Áp dụng thiết kế thực nghiệm Plackett-Burman và Latin Hypercube để chạy các mô hình thử nghiệm ban đầu.
    • Thông qua biểu đồ Tornado và hệ số Pearson, xác định các thông số có ảnh hưởng mang tính chi phối đến sản lượng khai thác và loại bỏ các yếu tố ít tác động.
  • Hiệu chỉnh mô hình (History Matching - HM):
    • Thực hiện quy trình tối ưu hóa hàm mục tiêu sai số.
    • Hiệu chỉnh thành công mô hình với lịch sử thử vỉa và khai thác thực tế của các giếng ST-A (tập E, tập F), ST-C (tập E, tập F) và ST-D (tập E) đối với các đường suy giảm áp suất đáy giếng (WBHP) và lưu lượng khai thác chất lưu (FGPR, FOPR).
  • Xây dựng bề mặt phản hồi (RSM) và mô phỏng Monte Carlo:
    • Thiết lập phương trình bề mặt phản hồi hồi quy bậc hai cho chỉ số tổng sản lượng khí khai thác mỏ cộng dồn (FGPT). Kiểm tra qua biểu đồ crossplot cho thấy mức độ tương thích cao giữa giá trị dự báo từ RSM và mô hình Eclipse.
    • Thực hiện 10.000 lần lặp Monte Carlo trên bề mặt phản hồi để thiết lập biểu đồ Histogram và đường cong tích lũy xác suất CDF.
    • Xuất ra khoảng tin cậy sản lượng khai thác khí và condensat cộng dồn tại các mức P10, P50, P90 tương ứng với các kịch bản phát triển cao, vừa và thấp.
  • Tối ưu hóa phương án phát triển mỏ:
    • Chạy mô phỏng tối ưu hóa số lượng và vị trí các giếng khoan phát triển (các giếng A, B, C, D) nhằm đạt hàm mục tiêu thu hồi cực đại.
    • Cập nhật và so sánh các thông số đánh giá từ mô hình với số liệu khai thác thử dài hạn thực tế tại mỏ ST-X, khẳng định độ tin cậy của quy trình xây dựng.

Kết quả và những đóng góp mới

Những đóng góp mới của luận án

  1. Đóng góp về mặt phương pháp luận: Đã nghiên cứu, xây dựng và hoàn thiện một quy trình tổng quát tích hợp đa phương pháp: từ nhận diện các yếu tố không chắc chắn qua tài liệu thử vỉa (DST), sàng lọc giảm chiều dữ liệu bằng thiết kế thực nghiệm (ED), xây dựng mô hình bề mặt phản hồi (RSM) đến mô phỏng Monte Carlo để định lượng dải sản lượng khai thác mỏ dầu khí.
  2. Tối ưu hóa tài nguyên tính toán: Việc ứng dụng mô phỏng Monte Carlo 10.000 lần lặp trên phương trình bề mặt phản hồi thay thế cho việc chạy trực tiếp mô hình mô phỏng thủy động lực học lớn giúp rút ngắn đáng kể thời gian và khối lượng tính toán trên máy tính, đồng thời đưa ra công cụ tính toán nhanh khi có bất kỳ thông số đầu vào nào thay đổi.
  3. Đóng góp thực tiễn cho mỏ ST-X: Xác định định lượng mức độ ảnh hưởng của các yếu tố bất định đến hệ số thu hồi mỏ ST-X; tính toán chính xác các cấp sản lượng tin cậy P10 – P50 – P90 cho hai tập cát kết E và F; đề xuất phương án bố trí giếng và kế hoạch phát triển mỏ tối ưu giúp cực đại hóa lượng thu hồi khí ngưng tụ và cực tiểu hóa rủi ro đầu tư.

Kiến nghị và giải pháp đề xuất

  • Áp dụng quy trình đánh giá bất định tích hợp đã xây dựng làm quy trình chuẩn để phân tích kinh tế và lập kế hoạch phát triển cho các mỏ khí ngưng tụ và mỏ cận biên có cấu trúc địa chất phức tạp khác tại thềm lục địa Việt Nam.
  • Trong quá trình triển khai khai thác thương mại mỏ ST-X, cần tiếp tục thu thập dữ liệu áp suất đáy giếng và lưu lượng khai thác thực tế liên tục để cập nhật, hiệu chỉnh lại phương trình bề mặt phản hồi theo thời gian thực.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Nghiên cứu được giới hạn dựa trên tập dữ liệu thu thập trong giai đoạn thăm dò – thẩm lượng ban đầu của mỏ ST-X (số lượng giếng khoan còn ít gồm ST-A, ST-B, ST-C, ST-D và thời gian thử vỉa ngắn hạn). Các yếu tố không chắc chắn về cơ chế liên thông thủy động lực học dài hạn giữa các khối đứt gãy xa chưa có đủ dữ liệu khai thác thực tế để kiểm chứng toàn diện.
  • Hướng nghiên cứu tiếp: Mở rộng tích hợp quy trình với mô hình kinh tế động (Dynamic Economic Model) nhằm tối ưu hóa trực tiếp chỉ số NPV và IRR dưới điều kiện biến động của giá dầu khí thế giới; đồng thời cập nhật mô hình khi mỏ bước vào giai đoạn khai thác toàn diện với hệ thống giếng khoan phát triển bổ sung.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên môn hữu ích cho:

  • Nghiên cứu sinh, học viên cao học và giảng viên thuộc các chuyên ngành Kỹ thuật Dầu khí, Địa chất Dầu khí, Kỹ thuật Vỉa và Quản lý Năng lượng: Tham khảo phương pháp tích hợp thiết kế thực nghiệm, mô hình bề mặt phản hồi và mô phỏng xác suất Monte Carlo vào bài toán mô phỏng khai thác mỏ.
  • Kỹ sư công nghệ mỏ, kỹ sư địa chất tại các công ty dầu khí (PVEP, Vietsovpetro, các JOCs): Ứng dụng quy trình sàng lọc tham số nhạy cảm và kỹ thuật hiệu chỉnh mô hình lịch sử (History Matching) để xử lý dữ liệu thử vỉa DST và giảm thiểu rủi ro khi lập kế hoạch phát triển mỏ cận biên (FDP).
  • Các nhà quản trị và người làm chính sách phát triển mỏ: Sử dụng kết quả đánh giá phân vị P10-P50-P90 làm cơ sở phê duyệt trữ lượng thương mại và quyết định mức vốn đầu tư phát triển mỏ hợp lý.

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

1. Các yếu tố không chắc chắn chính của mỏ khí ngưng tụ ST-X được nhận diện từ nguồn dữ liệu nào và bao gồm những tham số gì? Các yếu tố không chắc chắn được nhận diện chủ yếu từ kết quả phân tích số liệu động thử vỉa (DST), đo áp suất MDT và phân tích PVT tại các giếng thăm dò ST-A, ST-B, ST-C, ST-D. Các tham số chính gồm: độ thấm của vỉa ($k$), khoảng cách tới các biên chắn đứt gãy không thấm, bán kính và mức độ lắng đọng vây ngưng tụ condensate quanh đáy giếng, cùng sự biến thiên của các đặc tính lưu chất (GCR, $B_g, r_v$).

2. Phương pháp thiết kế thực nghiệm (ED) đóng vai trò gì trong việc xử lý mô phỏng vỉa mỏ ST-X? Thiết kế thực nghiệm (áp dụng phương pháp Plackett-Burman và Latin Hypercube) đóng vai trò sàng lọc và định lượng mức độ ảnh hưởng của các thông số vỉa đầu vào thông qua biểu đồ Tornado và hệ số tương quan Pearson. Phương pháp này giúp xác định những thông số chi phối nhiều nhất và loại bỏ các thông số ít tác động, qua đó giảm tối đa số lượng ca chạy mô phỏng Eclipse cần thiết.

3. Vì sao tác giả không chạy mô phỏng Monte Carlo trực tiếp trên mô hình mô phỏng vỉa mà phải thông qua bề mặt phản hồi (RSM)? Việc chạy 10.000 ca mô phỏng trực tiếp trên phần mềm mô phỏng thủy động lực học 3D (như Eclipse) đòi hỏi chi phí tính toán và thời gian rất lớn, không khả thi trong thực tế. Bằng cách thiết lập phương trình toán học bề mặt phản hồi (RSM) mô tả chính xác mối quan hệ giữa các biến đầu vào và sản lượng đầu ra, thuật toán Monte Carlo có thể thực hiện 10.000 phép thử ngẫu nhiên trên phương trình này chỉ trong vài giây mà vẫn đảm bảo độ chính xác cao.

4. Khoảng tin cậy P10, P50, P90 trong luận án có ý nghĩa gì đối với công tác phát triển mỏ ST-X? Các giá trị P10, P50, P90 thể hiện các cấp xác suất của tổng sản lượng khí và ngưng tụ khai thác tích lũy (FGPT, FOPT): P90 đại diện cho kịch bản thận trọng (xác suất đạt được ít nhất 90%), P50 đại diện cho kịch bản kỳ vọng/trung bình (xác suất đạt được 50%), và P10 đại diện cho kịch bản lạc quan/tiềm năng cao (xác suất đạt được 10%). Đây là căn cứ khoa học để xây dựng chiến lược khai thác tối ưu và dự phòng rủi ro kinh tế.

5. Kết quả thử vỉa DST ở hai tập cát kết E và F tại giếng khoan thăm dò đầu tiên ST-A cho thấy tiềm năng dòng chảy như thế nào? Tại giếng ST-A, kết quả thử vỉa DST ghi nhận tiềm năng lớn ở cả hai tập: DST#2 trong tập F đạt lưu lượng khí tối đa 32,0 triệu bộ khối khí/ngày và 3.604 thùng condensat/ngày; DST#3 trong tập E đạt lưu lượng khí cực đại 37,7 triệu bộ khối khí/ngày và 4.033 thùng condensat/ngày.

Kết luận

Luận án của tác giả Vũ Việt Hưng đã xây dựng thành công quy trình đánh giá định lượng các yếu tố không chắc chắn trong giai đoạn đầu phát triển mỏ khí ngưng tụ phức tạp thông qua việc tích hợp phân tích thử vỉa DST, thiết kế thực nghiệm, mô hình bề mặt phản hồi và mô phỏng Monte Carlo. Kết quả nghiên cứu đã giải quyết được bài toán cân bằng giữa việc tiết giảm thời gian tính toán mô phỏng số và việc nâng cao độ tin cậy trong dự báo khoảng sản lượng khai thác các cấp P10 – P50 – P90. Đây là công trình nghiên cứu hoàn chỉnh, đóng góp luận cứ khoa học và giải pháp thực tiễn quan trọng phục vụ cho việc lập phương án phát triển tối ưu mỏ khí ngưng tụ ST-X cũng như các mỏ dầu khí cận biên ngoài khơi Việt Nam.