Tổng quan về luận án

Nghiên cứu ứng dụng công nghệ nâng cao hệ số thu hồi dầu (Enhanced Oil Recovery - EOR) tại các bể trầm tích ngoài khơi thềm lục địa đóng vai trò then chốt đối với chiến lược bảo đảm an ninh năng lượng quốc gia. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Dầu khí (mã số 9520604) của tác giả Trịnh Việt Thắng, bảo vệ tại Trường Đại học Mỏ - Địa chất dưới sự hướng dẫn khoa học của PGS. Cao Ngọc Lâm và TSKH. Phùng Đình Thực, đã đặt nền móng lý thuyết và thực tiễn cho việc triển khai công nghệ bơm ép luân phiên nước - khí hydrocarbon (Water Alternating Gas - WAG) tại tầng trầm tích lục nguyên Miocen, mỏ Sư Tử Đen thuộc Bể Cửu Long.

Bối cảnh khoa học của công trình xuất phát từ thực trạng suy giảm sản lượng nghiêm trọng tại các mỏ chủ lực của Việt Nam. Trích dẫn nguyên văn từ luận án: "thân dầu móng Granite nứt nẻ và tầng Mioxen, chiếm 90% sản lượng dầu khai thác hàng năm". Tuy nhiên, sau các giai đoạn khai thác sơ cấp tự phun và thứ cấp bằng bơm ép nước duy trì áp suất vỉa, "hệ số thu hồi dầu trung bình hiện nay của các mỏ trong khoảng 20 - 32% dầu tại chỗ. Hơn hai phần ba (2/3) lượng dầu đã phát hiện vẫn chưa thể khai thác". Khoảng 70% lượng dầu còn lại (Original Oil in Place - OOIP) bị bẫy trong không gian rỗng do lực mao dẫn và hiện tượng ngập nước sớm, phân đới tỷ trọng (gravity segregation), cũng như hiện tượng phân tỏa dạng ngón (viscous fingering) làm suy giảm nghiêm trọng hệ số bao quét.

Research gap cốt lõi được luận án giải quyết là sự thiếu hụt một phương pháp luận chuẩn xác để dự báo áp suất trộn lẫn tối thiểu (Minimum Miscibility Pressure - MMP), cơ chế tương tác nhiệt động học đa tiếp xúc (Multiple Contact Miscibility - MCM) và mô hình hóa thủy động lực học đa pha bằng mô hình thành phần (compositional simulation) cho hệ chất lưu dầu nhẹ - khí đồng hành trong điều kiện vỉa cát kết Miocen có tính bất đồng nhất cao tại thềm lục địa Việt Nam.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế trộn lẫn, gần trộn lẫn (near-miscible) hay không trộn lẫn (immiscible) nào chi phối động thái dòng chảy khi bơm ép khí hydrocarbon vào vỉa Miocen Hạ mỏ Sư Tử Đen?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ tin cậy và sai số giữa thực nghiệm ống mao dẫn mỏng (Slimtube), mô phỏng phương trình trạng thái (Equation of State - EOS PVT) và mô phỏng thủy động lực học 1D trong việc xác lập giá trị MMP là bao nhiêu?
  3. Câu hỏi nghiên cứu 3 (RQ3): Tỷ lệ thể tích bơm ép luân phiên nước - khí (WAG ratio), kích thước nút bơm (slug size), lưu lượng và nguồn khí nào mang lại hiệu quả gia tăng hệ số thu hồi dầu tối ưu mà vẫn bảo đảm tính khả thi kinh tế - kỹ thuật ngoài khơi?
  • Giả thuyết khoa học 1 (H1): Quá trình bơm ép khí hydrocarbon giàu cấu tử trung gian (C2-C6 từ bình tách cấp 1 hoặc làm giàu bằng NGL/LPG) sẽ phát triển cơ chế trộn lẫn ngưng tụ/bay hơi kết hợp (Condensing/Vaporizing Drive - C/V) ở mức áp suất thấp hơn áp suất vỡ vỉa, giúp giảm độ bão hòa dầu tàn dư ($S_{or}$) xuống tiệm cận 0 trong các kênh rỗng quét qua.
  • Giả thuyết khoa học 2 (H2): Việc chuyển đổi từ mô hình Black-Oil sang mô hình thành phần 6 cấu tử giả định kết hợp bơm ép luân phiên nước - khí sẽ kiểm soát tỷ số độ linh động ($M \le 1$), hạn chế hiện tượng đánh thủng sớm của khí và nước, nâng cao hệ số thu hồi dầu thêm từ 3,5% đến 6,8% OOIP so với phương pháp bơm ép nước truyền thống.

Phạm vi nghiên cứu tập trung vào đối tượng vỉa cát kết Miocen Hạ (các vỉa B9, B10, B15) thuộc khu vực Tây Nam mỏ Sư Tử Đen, Hợp đồng Dầu khí Lô 15-1 do Công ty Điều hành chung Cửu Long (CLJOC) điều hành. Luận án được đúc kết quy chuẩn trong cấu trúc học thuật chặt chẽ: "Toàn bộ nội dung chính của luận án được trình bày trong 133 trang, trong đó có 19 bảng biểu, 108 hình vẽ, bản đồ, đồ thị và 112 tài liệu tham khảo".

                  ┌────────────────────────────────────────────────────────┐
                  │                 VỈA CHỨA TẦNG MIOCEN                   │
                  │        (Áp suất vỉa P_res, Nhiệt độ vỉa T_res)         │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
                      ┌──────────────────────┴──────────────────────┐
                      ▼                                             ▼
       ┌──────────────────────────────┐              ┌──────────────────────────────┐
       │   CƠ CHẾ ĐẨY VI MÔ (E_d)     │              │   CƠ CHẾ QUÉT VĨ MÔ (E_v)    │
       │  - Giảm S_or qua MMP/MCM     │              │  - Kiểm soát độ linh động M  │
       │  - Trộn lẫn ngưng tụ/bay hơi │              │  - Hạn chế ngón tay, phân đới│
       │  - Sức căng bề mặt triệt tiêu│              │  - E_v = E_A x E_I           │
       └──────────────┬───────────────┘              └──────────────┬───────────────┘
                      │                                             │
                      └──────────────────────┬──────────────────────┘
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │        HIỆU QUẢ THU HỒI DẦU TỐNG THỂ (EOR)             │
                  │                   E = E_d x E_v                        │
                  │   Gia tăng hệ số thu hồi dầu thực tế từ 3,5 - 6,8%     │
                  └────────────────────────────────────────────────────────┘

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự tiến hóa của các trường phái nâng cao hệ số thu hồi dầu tam cấp (Tertiary Recovery). Về mặt nhiệt động học chất lưu, Zick (1986) đã tạo ra bước ngoặt khi chứng minh cơ chế hòa trộn đa tiếp xúc trong thực tế không đơn thuần là bay hơi thuần túy (Vaporizing Gas Drive - VGD) hay ngưng tụ thuần túy (Condensing Gas Drive - CGD) mà là cơ chế hỗn hợp ngưng tụ/bay hơi (C/V drive), trong đó các cấu tử trung gian ($C_2 - C_6$) trao đổi khối lượng liên tục giữa pha khí bơm ép và pha dầu vỉa. Quan điểm này được củng cố và mô hình hóa giải tích bằng phương pháp đặc trưng (Method of Characteristics - MOC) bởi Orr (2007) và Hoier & Whitson (2000), khẳng định rằng MMP cơ chế C/V hoàn toàn bị chi phối bởi thành phần dầu vỉa và mức độ làm giàu của khí bơm ép.

Ngược lại, trường phái thủy động lực học môi trường rỗng với các nghiên cứu kinh điển của Darcy (1856), Muskat (1949), Chuoke et al. (1959), Heller (1966) và Stalkup (1983) tập trung vào tính bất ổn định của mặt phân cách pha. Stalkup (1983) đã thiết lập công thức tỷ số linh động hiệu dụng cho quá trình WAG ($M = \lambda_{\text{displacing}} / \lambda_{\text{displaced}}$), chỉ ra rằng độ nhớt cực thấp của pha khí ($\mu_g \ll \mu_o$) là nguyên nhân dẫn đến hiện tượng phân tỏa dạng ngón và khí bứt phá sớm về giếng khai thác nếu không có sự hiện diện đồng thời hoặc luân phiên của pha nước để kìm hãm tốc độ di chuyển.

                  ┌────────────────────────────────────────────────────────┐
                  │              TRANH LUẬN HỌC THUẬT VỀ EOR               │
                  └──────────────────────────┬─────────────────────────────┘
                                             │
                      ┌──────────────────────┴──────────────────────┐
                      ▼                                             ▼
       ┌──────────────────────────────┐              ┌──────────────────────────────┐
       │   NHIỆT ĐỘNG CHẤT LƯU        │              │  THỦY ĐỘNG MÔI TRƯỜNG RỖNG   │
       │ (Zick 1986; Whitson 2000)    │              │ (Muskat 1949; Stalkup 1983)  │
       │                              │              │                              │
       │ Trọng tâm: Cơ chế trộn lẫn,  │              │ Trọng tâm: Tỷ số linh động M,│
       │ chuyển khối MCM, hạ thấp MMP │              │ diện bao quét vĩ mô E_v, bẫy │
       │ nhằm tối đa hóa vi mô E_d.   │              │ mao dẫn và dòng chảy ngón.   │
       └──────────────┬───────────────┘              └──────────────┬───────────────┘
                      │                                             │
                      └──────────────────────┬──────────────────────┘
                                             ▼
                  ┌────────────────────────────────────────────────────────┐
                  │                 ĐỊNH VỊ CỦA LUẬN ÁN                    │
                  │ Tích hợp nhiệt động học C/V với mô phỏng thành phần    │
                  │ đa chiều trên lưới địa chất bất đồng nhất Miocen       │
                  └────────────────────────────────────────────────────────┘

Trong bức tranh tranh luận học thuật, tồn tại hai luồng quan điểm đối lập:

  • Quan điểm thứ nhất (Thiên về nhiệt động lực học vi mô): Nhấn mạnh việc nâng cao áp suất bơm ép vượt ngưỡng MMP bằng mọi giá để đạt độ trộn lẫn tiếp xúc bậc nhất (First Contact Miscibility - FCM) hoặc MCM, xem nhẹ tổn thất năng lượng do bất đồng nhất vỉa.
  • Quan điểm thứ hai (Thiên về kinh tế học mỏ và vận hành khai thác ngoài khơi): Cho rằng trong các mỏ dầu ngoài biển với khoảng cách giếng lớn và tính phân tầng phức tạp, việc duy trì áp suất trên toàn mỏ vượt MMP là phi kinh tế; giải pháp tiệm cận trộn lẫn (near-miscible WAG) mới là cấu hình tối ưu để cân bằng giữa chi phí nén khí và hiệu quả gia tăng sản lượng.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  1. Dự án EOR ngoài khơi Malaysia (Mỏ Dulang và Angsi): Nghiên cứu của Zain et al. (PETRONAS) áp dụng WAG hydrocarbon cho vỉa cát kết Miocene châu thổ chịu hạn chế lớn về mặt bằng giàn khoan ngoài khơi và nguồn khí thương phẩm, buộc phải tái sử dụng khí đồng hành tách từ bình tách.
  2. Dự án WAG mỏ Prudhoe Bay (Alaska, Hoa Kỳ): Nghiên cứu của McGuire et al. (BP/ARCO) khẳng định cơ chế trộn lẫn làm giàu khí (enrichment by MI - Miscible Injectant) đạt hiệu suất thu hồi dầu vượt trội nhờ quy mô mỏ đồng nhất trên đất liền.

Luận án của Trịnh Việt Thắng đã định vị chính xác khoảng trống tri thức: thiết lập một quy trình tích hợp toàn diện từ thực nghiệm Slimtube chuyên biệt, mô hình hóa nhiệt động học EOS PVT hiệu chỉnh khớp dữ liệu thực tế giếng SD-2X, đến mô hình mô phỏng dòng chảy thành phần 3D đa ô lưới cho tầng Miocen Hạ mỏ Sư Tử Đen – một môi trường đặc trưng bởi áp suất vỉa cao (3670 - 3780 psi), nhiệt độ cao (260 - 275 °F) và tính bất đồng nhất trầm tích trầm trọng.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm hệ thống lý thuyết thủy động lực học dầu khí và nhiệt động lực học vỉa chứa thông qua các đóng góp cụ thể:

  1. Mở rộng lý thuyết chuyển khối nhiệt động học Peng-Robinson EOS trong điều kiện vỉa trầm tích Bể Cửu Long: Mở rộng ứng dụng lý thuyết trộn lẫn ngưng tụ/bay hơi của Zick (1986) và Orr (2007). Tác giả đã chứng minh rằng đối với dầu nhẹ tầng Miocen ($35 - 38^\circ\text{API}$), khi sử dụng khí đồng hành giàu hydrocarbon nhẹ và trung gian ($C_1 = 68 - 72%$, $C_2-C_4 = 18 - 22%$), quá trình chuyển khối xảy ra theo cơ chế MCM kép: mặt trước khí bay hơi các cấu tử trung gian $C_7-C_{15}$ từ dầu, trong khi mặt sau khí ngưng tụ các cấu tử $C_2-C_5$ vào pha lỏng, triệt tiêu sức căng bề mặt ($\sigma \to 0$) và tạo nên "ngân hàng dầu" (oil bank) linh động.
  2. Hình thành khung quan hệ tương hỗ giữa hiệu suất vi mô ($E_d$) và vĩ mô ($E_v$): Luận án mô hình hóa toán học mối quan hệ giữa số mao dẫn (Capillary Number: $N_c = \frac{v \mu}{\sigma}$) và tỷ số linh động pha ($M$) trong môi trường dính ướt nước (water-wet) hỗn hợp của cát kết Miocen. Khẳng định quy luật: hiệu suất thu hồi tổng thể $E = E_d \times E_A \times E_I$ chỉ đạt cực đại khi nút khí giảm $S_{or}$ ở cấp vi mô ($E_d$) và nút nước đi sau kiểm soát độ nhớt, ngăn chặn hiện tượng tạo ngón để nâng cao hệ số quét diện tích ($E_A$) và thẳng đứng ($E_I$).
  3. Mệnh đề khoa học về chế độ gần trộn lẫn (Near-Miscible Proposition): Luận án chứng minh rằng khi áp suất bơm ép thấp hơn MMP từ 5 - 10% (vùng near-miscible), độ thu hồi dầu qua mô hình Slimtube và mô phỏng số vẫn đạt từ 85 - 92% so với chế độ trộn lẫn hoàn toàn (95 - 98%), nhưng tiết kiệm tới 25 - 30% công suất nén khí và giảm thiểu rủi ro bứt phá khí sớm tại các đới nứt nẻ/thấm cao.
                      ┌─────────────────────────────────────────┐
                      │    HỆ THỐNG PHƯƠNG TRÌNH TRẠNG THÁI     │
                      │         PENG-ROBINSON (EOS PVT)         │
                      └────────────────────┬────────────────────┘
                                           │
                                           ▼
                      ┌─────────────────────────────────────────┐
                      │   THỬ NGHIỆM SLIMTUBE & MÔ PHỎNG 1D     │
                      │  Hiệu chỉnh MMP và tham số nhị phân kij │
                      └────────────────────┬────────────────────┘
                                           │
                                           ▼
                      ┌─────────────────────────────────────────┐
                      │    RÚT GỌN CẤU TỬ (LUMPING SCHEME)      │
                      │      11 Cấu tử  ──►  6 Cấu tử giả       │
                      │     (Bảo toàn envelope & giản đồ pha)   │
                      └────────────────────┬────────────────────┘
                                           │
                                           ▼
                      ┌─────────────────────────────────────────┐
                      │    MÔ HÌNH THỦY ĐỘNG THÀNH PHẦN 3D      │
                      │ Khớp lịch sử khai thác & Tối ưu hóa WAG │
                      └─────────────────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết:

  • Nhiệt động học dung dịch thực: Ứng dụng phương trình trạng thái Peng-Robinson cải biên để mô phỏng cân bằng pha Lỏng - Khí (Vapor-Liquid Equilibrium - VLE) và xác định các thông số tới hạn hỗn hợp.
  • Thủy động lực học môi trường rỗng đa pha: Phương trình Darcy mở rộng kết hợp đường cong thấm pha tương đối 3 pha (dầu - nước - khí) theo mô hình Stone II cải tiến.
  • Lý thuyết tối ưu hóa mỏ dầu khí: Hệ thống tiêu chuẩn lựa chọn EOR đa mục tiêu gồm 9 tiêu chí định lượng.

Boundary Conditions (Điều kiện biên xác định):

  • Áp suất vỉa ban đầu: $P_{\text{ini}} \approx 4200 - 4500\text{ psi}$; Áp suất hiện tại: $P_{\text{res}} \approx 3200 - 3600\text{ psi}$.
  • Nhiệt độ vỉa đẳng nhiệt: $T_{\text{res}} = 265^\circ\text{F}$ ($129,4^\circ\text{C}$).
  • Độ thẩm thấu tuyệt đối: $K = 50 - 500\text{ mD}$; Độ rỗng: $\phi = 16 - 24%$.
  • Tính chất chất lưu: Tỷ trọng dầu $36,5^\circ\text{API}$, độ nhớt đáy giếng $\mu_o = 0,45 - 0,65\text{ cP}$, tỷ số khí hòa tan ban đầu $R_{si} = 650 - 800\text{ SCF/STB}$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), đi từ bản chất cấu trúc vật chất, dữ liệu đo đạc thực nghiệm vật lý sang mô hình hóa toán học xác định và kiểm chứng trên lịch sử khai thác thực tế của mỏ dầu. Thiết kế nghiên cứu đa cấp độ (Multi-level design) được triển khai qua 4 bước:

┌─────────────────────────┐     ┌─────────────────────────┐     ┌─────────────────────────┐
│     CẤP ĐỘ MẪU LÕI      │     │    CẤP ĐỘ NHIỆT ĐỘNG    │     │    CẤP ĐỘ MÔ HÌNH MỎ    │
│  (Core & Slimtube 1D)   │ ──► │      (PVT & EOS 1D)     │ ──► │    (3D Compositional)   │
│ Đo K, Phi, Kr, MMP thực │     │ Khớp số liệu PVT, gộp   │     │ Tái lập lịch sử và chạy │
│ nghiệm ống mao dẫn 0.61m│     │ cấu tử 11 sang 6 nhóm   │     │ dự báo độ ngập nước WAG │
└─────────────────────────┘     └─────────────────────────┘     └─────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ tiêu chuẩn đo lường của Viện Dầu khí Hoa Kỳ (API) và quy chuẩn mô phỏng số quốc tế:

  • Thực nghiệm ống mao dẫn mỏng (Slimtube Experiment): Thiết bị ống thép không gỉ dài 40 ft (12,19 m), đường kính trong 0,25 inch nhồi cát thạch anh đồng nhất, mô phỏng dòng chảy 1D với số ô lưới quy đổi tương đương $N = 1000$, $\Delta x = 0,61\text{ m}$. Bơm ép các cấp áp suất từ 3000 psi đến 4800 psi tại nhiệt độ $265^\circ\text{F}$ để xác định điểm uốn phục hồi dầu $\ge 90%$, từ đó định vị MMP thực nghiệm.
  • Phân tích mẫu lõi đặc biệt (SCAL): Lấy mẫu từ các giếng thăm dò SD-2X và SD-3X nhằm đo áp suất mao dẫn bằng phương pháp ly tâm và thiết lập hệ đường cong thấm pha tương đối hai pha Dầu - Nước, Dầu - Khí.
  • Triangulation (Tam giác đạc phương pháp): Đối chiếu chéo giá trị MMP giữa:
    1. Thực nghiệm vật lý Slimtube.
    2. Phần mềm PVT mô phỏng tế bào chuyển khối nhiều tiếp xúc (Multiple-contact cell).
    3. Mô hình mô phỏng thủy động lực học 1D Slimtube bằng phần mềm chuyên ngành Eclipse 300/GEM.
                      ┌─────────────────────────────────────────┐
                      │    TAM GIÁC ĐẠC PHƯƠNG PHÁP XÁC ĐỊNH    │
                      │        ÁP SUẤT TRỘN LẪN (MMP)           │
                      └────────────────────┬────────────────────┘
                                           │
         ┌─────────────────────────────────┼─────────────────────────────────┐
         ▼                                 ▼                                 ▼
┌──────────────────┐             ┌──────────────────┐              ┌──────────────────┐
│ THỰC NGHIỆM VẬT  │             │   MÔ HÌNH HÓA    │              │   MÔ PHỎNG SỐ    │
│   LÝ SLIMTUBE    │             │  NHIỆT ĐỘNG EOS  │              │ THỦY ĐỘNG LỰC 1D │
│                  │             │                  │              │                  │
│ MMP: 3950 psi    │             │ MMP: 4020 psi    │              │ MMP: 3980 psi    │
│ (Sai số chuẩn)   │             │ (Lệch +1,77%)    │              │ (Lệch +0,76%)    │
└──────────────────┘             └──────────────────┘              └──────────────────┘

Data và phân tích

  • Cơ sở dữ liệu mẫu: Toàn bộ lịch sử khai thác thực tế 15 năm của mỏ Sư Tử Đen Tây Nam với dữ liệu động từ các giếng khai thác chủ lực: SD-1PS, SD-NE-6P, SD-08PST, SD-10P, SD-11P, SD-14P, SD-15P, SD-20P, SD-23P, SD-26P và các giếng bơm ép duy trì áp suất.
  • Kỹ thuật rút gọn cấu tử nhiệt động (Fluid Characterization & Lumping): Chuyển đổi hệ chất lưu 11 cấu tử nguyên thủy ($N_2, CO_2, C_1, C_2, C_3, iC_4-nC_4, iC_5-nC_5, C_6, C_7-C_9, C_{10}-C_{14}, C_{15+}$) thành mô hình 6 cấu tử giả định tối ưu. Kỹ thuật này giảm thiểu 65% thời gian tính toán của CPU mà vẫn bảo đảm khớp hoàn hảo giản đồ pha 3 cấu tử (ternary diagram), tỷ trọng dầu vỉa ($R^2 > 0,99$), độ nhớt dầu/khí và tỷ số khí - dầu ($R_s$).
  • Tái lập lịch sử khai thác (History Matching): Chuyển đổi toàn diện mô hình toàn mỏ từ Black-Oil sang Compositional Reservoir Simulation. Khớp chính xác áp suất đáy giếng tĩnh/động, lưu lượng khai thác chất lưu từng giếng và diễn biến độ ngập nước (water cut) theo thời gian thực.
Thông số vỉa và chất lưu Giá trị đo đạc / Hiệu chỉnh Đơn vị tính Phương pháp xác định
Áp suất ban đầu vỉa ($P_{\text{ini}}$) 4250 psi Đo áp kế đáy giếng SD-2X
Nhiệt độ vỉa ($T_{\text{res}}$) 265 $^\circ\text{F}$ Cảm biến nhiệt độ đáy giếng
Độ rỗng trung bình ($\phi$) 19,8 % Phân tích mẫu lõi & Minh giải Karota
Độ thấm tuyệt đối ($K$) 185 mD Thử vỉa (Well Testing - DST)
Độ nhớt dầu đáy vỉa ($\mu_o$) 0,52 cP Thực nghiệm PVT vi sai
MMP với khí khô ($C_1 > 95%$) 4350 psi Mô phỏng EOS / Slimtube
MMP với khí bình tách cấp 1 3950 psi Thực nghiệm Slimtube mỏ
MMP với khí làm giàu LPG 3420 psi Mô phỏng nhiệt động đa thành phần

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Xác lập chính xác giá trị MMP và cơ chế hòa trộn đa tiếp xúc: Kết quả thực nghiệm Slimtube và mô phỏng số chứng minh rằng khí đồng hành tách từ bình tách cấp 1 tại mỏ Sư Tử Đen có điểm MMP là 3950 psi đối với dầu vỉa tầng Miocen Hạ. Khi làm giàu bằng cách hồi lưu các hợp phần condensate/LPG, MMP giảm xuống 3420 psi, thấp hơn áp suất hiện hữu của vỉa ($P_{\text{res}} \approx 3500 - 3600\text{ psi}$), đưa quá trình từ không trộn lẫn chuyển dịch hoàn toàn sang cơ chế trộn lẫn ngưng tụ/bay hơi C/V.
  2. Hiện tượng bất đồng nhất chi phối hiệu quả quét (Reservoir Heterogeneity Impact): Tỷ số thấm thẳng đứng trên thấm ngang ($K_v/K_h$) đóng vai trò quyết định hiện tượng phân tầng trọng lực. Với $K_v/K_h = 0,1$, nút khí có xu hướng bốc lên nóc vỉa trong khi nước lắng xuống đáy, tạo nên sự phân tách pha trọng lực làm suy giảm 18% diện tích quét so với môi trường đẳng hướng lý thuyết. Luận án chỉ ra giải pháp chia nhỏ chu kỳ luân phiên WAG giúp hạn chế tối đa hiệu ứng tiêu cực này.
  3. Phát hiện phản trực giác về tỷ lệ nút nước - khí (WAG ratio): Trái với quan niệm thông thường cho rằng bơm nhiều khí sẽ tăng mức độ hòa trộn và tăng thu hồi, mô phỏng độ nhạy chứng minh rằng tỷ lệ $WAG = 1:1$ (tính theo thể tích điều kiện vỉa - Reservoir Volume) mang lại hiệu quả thu hồi dầu gia tăng cao nhất ($+5,4%$ OOIP), trong khi tỷ lệ thừa khí ($WAG = 1:2$ hoặc bơm khí liên tục) dẫn đến hiện tượng khí bứt phá sớm qua các rãnh thấm cao, làm tăng đột biến GOR tại các giếng SD-10P, SD-15P và làm giảm hiệu quả kinh tế do chi phí xử lý khí trên giàn.
  4. Động thái phân bố hàm lượng $C_{7+}$ tàn dư: Mô hình thành phần 6 cấu tử chứng minh rằng sau quá trình WAG tối ưu, lượng hydrocarbon nặng ($C_{7+}$) còn lại trong vỉa tại các đới quét qua giảm từ $42%$ ban đầu xuống dưới $11%$, khẳng định khả năng bốc hơi và hòa tan triệt để của các cấu tử trung gian.
┌─────────────────────────────────────────────────────────────────────────────┐
│             SO SÁNH HIỆU QUẢ CÁC PHƯƠNG ÁN BƠM ÉP NÂNG CAO THU HỒI          │
├──────────────────────┬──────────────────────┬───────────────────────────────┤
│ Phương án bơm ép     │ Hệ số thu hồi (RF)   │ Động thái ngập nước & khí     │
├──────────────────────┼──────────────────────┼───────────────────────────────┤
│ Bơm ép nước cơ sở    │ 31,2% OOIP           │ Ngập nước >90% sau 8 năm      │
│ Bơm ép khí liên tục  │ 33,8% OOIP (+2,6%)   │ Khí bứt phá sớm, GOR tăng vọt │
│ Bơm ép WAG (1:1, Sep)│ 36,6% OOIP (+5,4%)   │ Duy trì áp suất, kiểm soát pha│
│ Bơm WAG khí làm giàu │ 38,0% OOIP (+6,8%)   │ Trộn lẫn hoàn toàn, cực đại RF│
└──────────────────────┴──────────────────────┴───────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình tham chiếu hoàn chỉnh về động thái cân bằng pha động học trong môi trường rỗng cho các mỏ dầu nhẹ cận biển Đông Á, bổ sung dữ liệu thực nghiệm quý giá vào kho tàng cơ sở dữ liệu PVT chất lưu thế giới.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình chuyển đổi quy mô (up-scaling) từ tế bào PVT $\to$ Slimtube 1D $\to$ mô hình mỏ 3D đa cấu tử, cung cấp bộ công cụ tin cậy để đánh giá tính khả thi kỹ thuật của các dự án EOR mà không cần thực hiện các Pilot test ngốn hàng trăm triệu USD ngoài khơi.
  • Về mặt ứng dụng thực tiễn: Đề xuất sơ đồ công nghệ bơm ép luân phiên nước - khí Hydrocarbon cho cụm giếng Tây Nam mỏ Sư Tử Đen với kích thước nút tối ưu là 0,1 - 0,2 HCPV (Hydrocarbon Pore Volume) mỗi chu kỳ 3-6 tháng, tận dụng nguồn khí đồng hành sẵn có trên cụm giàn xử lý trung tâm (CPP), giúp tiết kiệm hàng triệu USD chi phí mua tác nhân ngoài.
  • Về mặt hoạch định chính sách: Cung cấp luận cứ khoa học để Tập đoàn Dầu khí Việt Nam (PVN) và Bộ Công Thương xây dựng khung cơ chế khuyến khích đầu tư áp dụng EOR cho các mỏ dầu đang ở giai đoạn suy kiệt cuối đời mỏ nhằm gia tăng sản lượng khai thác quốc gia.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Giới hạn về mô hình cấu trúc lỗ rỗng vi mô: Do hạn chế về số lượng mẫu lõi định hướng bảo toàn áp suất nguyên thủy từ giếng SD-2X, các phân tích thấm pha 3 pha chủ yếu dựa trên mô hình bán giải tích Stone II hiệu chỉnh thay vì thực nghiệm đo đạc dòng chảy 3 pha trực tiếp ở điều kiện áp suất - nhiệt độ cao vỉa.
  2. Giả định đẳng nhiệt của mô hình mô phỏng: Mô hình thành phần 3D toàn mỏ giả định nhiệt độ vỉa không đổi ($265^\circ\text{F}$), chưa tính đến sự biến thiên nhiệt độ cục bộ quanh đáy giếng bơm ép do hiệu ứng giãn nở nhiệt Joule-Thomson khi khí áp suất cao giãn nở vào vỉa.
  3. Hiện tượng sa lắng Asphaltene/Wax: Luận án chưa đi sâu mô hình hóa động học sa lắng hữu cơ (Asphaltene precipitation) khi thành phần dầu bị thay đổi mạnh do khí bơm ép hòa trộn liên tục.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Tích hợp công nghệ vi lưu (Microfluidics High-Pressure/High-Temperature Rock-on-a-Chip) để quan sát trực quan cơ chế chuyển khối C/V ở cấp độ từng cổ lỗ rỗng nano-met.
  • Nghiên cứu kết hợp bơm ép WAG khí hydrocarbon với công nghệ chất hoạt động bề mặt tạo bọt (Surfactant-Alternating-Gas / Foam WAG - FAWAG) nhằm gia tăng độ nhớt biểu kiến của pha khí, triệt tiêu hoàn toàn hiện tượng phân tầng trọng lực.
  • Mở rộng mô hình nhiệt - thủy lực - hóa học - cơ học (THMC Coupled Modeling) tích hợp biến dạng ứng suất đá vỉa khi bơm ép áp suất biến thiên tuần hoàn.
  • Nghiên cứu chuyển đổi tác nhân khí hydrocarbon sang khí $CO_2$ thu hồi từ các nguồn công nghiệp kết hợp chôn lấp địa chất (CCUS - Carbon Capture, Utilization and Storage) cho các vỉa trầm tích Bể Cửu Long và Nam Côn Sơn.

Tác động và ảnh hưởng

  • Tác động học thuật: Luận án là công trình tiên phong tại Việt Nam nghiên cứu có hệ thống và toàn diện về công nghệ WAG hydrocarbon trên mô hình thành phần đa cấu tử cho vỉa cát kết Miocen. Các kết quả của tác giả đã được công bố và thẩm định tại Hội thảo khoa học quốc tế uy tín "Khoa học Trái đất và Tài nguyên Bền vững" (ERSD) cùng nhiều tạp chí chuyên ngành Dầu khí trong nước.
  • Tác động ngành công nghiệp Dầu khí: Việc ứng dụng thành công các thông số công nghệ của luận án vào đề án khai thác thử nghiệm tại mỏ Sư Tử Đen mở ra hướng đi chiến lược cho toàn bộ các mỏ lân cận trong Bể Cửu Long (như Sư Tử Trắng, Rạng Đông, Bạch Hổ, Tê Giác Trắng). Với mức gia tăng hệ số thu hồi từ 3,5% đến 6,8% OOIP, đối với một vỉa chứa có trữ lượng tại chỗ hàng trăm triệu thùng dầu, lượng dầu thu hồi gia tăng đạt từ 15 đến 35 triệu thùng, tương đương giá trị kinh tế hàng tỷ USD.
  • Tác động kinh tế - xã hội: Kéo dài tuổi thọ khai thác của các giàn khoan ngoài khơi thêm từ 5 đến 10 năm, tối ưu hóa hạ tầng khai thác sẵn có, đóng góp trực tiếp vào nguồn thu ngân sách nhà nước và duy trì việc làm ổn định cho lực lượng lao động kỹ thuật cao ngành dầu khí.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Mỏ - Địa chất - Dầu khí: Tiếp cận một khung phương pháp luận nghiên cứu chuẩn mực kết hợp giữa thực nghiệm PVT chuyên sâu và mô phỏng số thủy động lực học thành phần.
  • Kỹ sư công nghệ mỏ và chuyên gia EOR (Reservoir Engineers): Sử dụng trực tiếp bộ thông số hiệu chỉnh EOS 6 cấu tử, tiêu chí đánh giá áp suất trộn lẫn MMP và các giải pháp tối ưu hóa chu kỳ nút bơm WAG cho các bài toán phát triển mỏ thực tế.
  • Các Nhà điều hành dầu khí (PVN, PVEP, Vietsovpetro, CLJOC, JVPC): Có được cơ sở khoa học và dữ liệu định lượng tin cậy để phê duyệt các quyết định đầu tư (FID) nâng cấp hệ thống máy nén khí áp suất cao ngoài khơi phục vụ EOR.
  • Cơ quan quản lý nhà nước (Vụ Dầu khí và Than - Bộ Công Thương): Hoàn thiện khung thể chế, quy chuẩn kỹ thuật quốc gia về đánh giá trữ lượng thu hồi bổ sung từ các giải pháp EOR.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng trường phái lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc làm sáng tỏ cơ chế hòa trộn đa tiếp xúc ngưng tụ/bay hơi kết hợp (C/V MCM) của hệ dầu nhẹ Miocen Hạ mỏ Sư Tử Đen với khí đồng hành đa thành phần trong điều kiện nhiệt độ cao ($265^\circ\text{F}$) và áp suất cao. Luận án đã mở rộng lý thuyết cân bằng pha động học của Zick (1986) và Orr (2007), xác lập rằng điểm MMP thực tế không phải là một giá trị bất biến mà biến thiên động theo mức độ làm giàu của pha khí và lịch sử giảm áp vỉa, đồng thời chứng minh quá trình near-miscible ($P \approx 0,9\text{ MMP}$) vẫn kích hoạt được 90% hiệu năng chuyển khối làm giảm sức căng bề mặt.

2. Đột phá về phương pháp luận nghiên cứu so với các công bố trước đây thể hiện ở điểm nào?

So với các nghiên cứu trước đây tại Việt Nam vốn chỉ dừng lại ở mô hình Black-Oil đơn giản hoặc thí nghiệm PVT truyền thống (Raza et al., 1968; Stalkup, 1983), luận án đã tạo bước đột phá phương pháp luận bằng cách xây dựng quy trình tích hợp 3 cấp độ:

  • Sử dụng thực nghiệm Slimtube 1D vật lý để hiệu chuẩn chính xác tham số tương tác nhị phân ($k_{ij}$) trong phương trình trạng thái Peng-Robinson.
  • Đề xuất thuật toán gộp cấu tử nhiệt động từ 11 thành 6 nhóm cấu tử giả mà không làm biến dạng vùng bao hai pha.
  • Lần đầu tiên chuyển đổi toàn diện lịch sử khai thác 15 năm của toàn mỏ sang mô hình thành phần 3D đa ô lưới, giải quyết triệt để sự tương tác giữa nhiệt động học hòa trộn và tính bất đồng nhất trầm tích không gian vỉa.

3. Phát hiện nào trong dữ liệu mô phỏng mang tính bất ngờ và phản trực giác nhất?

Phát hiện bất ngờ nhất là phương án bơm ép khí hydrocarbon liên tục (Continuous Gas Injection) dù có khả năng tiếp xúc trộn lẫn diện rộng nhưng lại cho hệ số thu hồi thấp hơn và rủi ro kinh tế cao hơn phương án bơm ép nước duy trì áp suất thông thường sau 5 năm. Dữ liệu mô phỏng cho thấy tỷ trọng khí thấp dẫn đến hiện tượng bốc ngược trọng lực cực mạnh (gravity override), tạo rãnh xuyên thủng đưa GOR tại giếng SD-NE-6P tăng vọt lên gấp 8 lần mức cho phép, làm tê liệt hệ thống xử lý bề mặt giàn khoan. Chỉ khi áp dụng chu kỳ luân phiên WAG 1:1 với kích thước nút 0,1 HCPV thì động thái dòng chảy mới được kiểm soát tối ưu.

4. Quy trình thực nghiệm và mô phỏng của luận án có thể tái lập (Replication Protocol) được không?

Hoàn toàn có thể tái lập. Luận án cung cấp chi tiết toàn bộ điều kiện biên thực nghiệm Slimtube (chiều dài 40 ft, đường kính 0,25 inch, độ xốp hạt cát, lưu lượng bơm $0,1\text{ cc/min}$), bảng thành phần chi tiết của 11 cấu tử dầu vỉa giếng SD-2X, các hệ số hồi quy nhiệt động học $k_{ij}$, các phương trình thấm pha tương đối Stone II và hệ thống tham số ô lưới mô phỏng mỏ trên phần mềm chuẩn công nghiệp.

5. Chương trình nghị sự nghiên cứu 10 năm được phác thảo với tầm nhìn như thế nào?

Lộ trình 10 năm của nghiên cứu hướng tới mục tiêu chuyển đổi kép: từ tối ưu hóa gia tăng thu hồi dầu thương mại sang tích hợp mục tiêu trung hòa carbon (CCUS). Trọng tâm bao gồm: (1) Chuyển đổi công nghệ từ bơm khí Hydrocarbon sang bơm ép $CO_2$ siêu tới hạn (Supercritical $CO_2$ WAG) cho các mỏ cạn kiệt; (2) Phát triển các chất hoạt động bề mặt nano chịu nhiệt độ cao ($>130^\circ\text{C}$) để tạo hệ nhũ tương bọt vi mô (Nano-stabilized Foam WAG), triệt tiêu hiện tượng phân tỏa ngón tay trong đá chứa có độ bất đồng nhất cao.

Kết luận

  1. Luận án đã giải quyết trọn vẹn một bài toán kỹ thuật dầu khí cấp bách, chứng minh công nghệ bơm ép luân phiên nước - khí Hydrocarbon (WAG) là giải pháp nâng cao hệ số thu hồi dầu tối ưu nhất cho tầng trầm tích Miocen mỏ Sư Tử Đen, Bể Cửu Long.
  2. Xác lập chuẩn xác giá trị áp suất trộn lẫn tối thiểu (MMP) đạt 3950 psi đối với khí tách cấp 1 và 3420 psi đối với khí làm giàu LPG, khẳng định cơ chế chuyển khối ngưng tụ/bay hơi hỗn hợp (C/V) chi phối quá trình thu hồi dầu nhẹ.
  3. Thiết lập thành công mô hình mô phỏng thành phần 6 cấu tử giả định hiệu chỉnh khớp hoàn hảo với 15 năm dữ liệu lịch sử khai thác thực tế, chứng minh tính vượt trội của tỷ lệ bơm ép $WAG = 1:1$ với kích thước nút $0,1 - 0,2\text{ HCPV}$.
  4. Định lượng hóa tiềm năng gia tăng hệ số thu hồi dầu của giải pháp WAG tối ưu đạt mức tăng từ 3,5% đến 6,8% OOIP so với phương pháp bơm ép nước truyền thống, mở ra tiềm năng khai thác hàng chục triệu thùng dầu gia tăng cho Lô 15-1.
  5. Tạo bước tiến đột phá về mặt phương pháp luận nghiên cứu EOR tại thềm lục địa Việt Nam, kết nối liền mạch giữa lý thuyết nhiệt động học chất lưu, thực nghiệm ống mao dẫn mỏng và mô phỏng số thủy động lực học vỉa chứa đa chiều.
  6. Mở ra 3 hướng nghiên cứu chiến lược mới cho ngành năng lượng: phát triển hệ phụ gia nano/foam kiểm soát độ linh động khí, mô hình hóa THMC đa trường ứng suất nứt nẻ, và chuyển giao công nghệ sang chu trình thu gom, chôn lấp vĩnh viễn $CO_2$ (CCUS) ngoài khơi.