Tổng quan về luận án

Nghiên cứu "Đánh giá tiềm năng thấm chứa dầu khí trầm tích điện trở thấp lô 16-1 bể Cửu Long" của tác giả Bùi Hữu Phước (Chuyên ngành Kỹ thuật địa chất, Mã số: 95.501, Trường Đại học Mỏ - Địa chất Hà Nội, 2021) giải quyết một trong những thách thức kỹ thuật - địa chất phức tạp nhất trong ngành công nghiệp dầu khí tại thềm lục địa Việt Nam: hiện tượng vỉa chứa dầu mang giá trị điện trở suất thấp (Low Resistivity Low Contrast - LRLC).

Trong bối cảnh bể Cửu Long đang khai thác dầu khí chủ yếu từ khối đá móng Granitoid nứt nẻ trước Kainozoi (chiếm khoảng 80% tổng sản lượng, tương đương 20 triệu tấn/năm), các tập trầm tích Mioxen dưới (hệ tầng Bạch Hổ) và Oligoxen trên (hệ tầng Trà Tân) nổi lên như đối tượng khai thác chiến lược thứ hai với đặc tính thấm chứa vượt trội. Tuy nhiên, các vỉa chứa dầu tại các mỏ như Tê Giác Trắng (TGT), Hải Sư Trắng, Rồng, Ruby lại thể hiện giá trị điện trở suất thực đo sâu rất thấp, chỉ dao động trong khoảng $1{,}5 - 3{,}0\ \Omega\cdot\text{m}$. Giá trị này gần như tương đương với các vỉa chứa nước hoàn toàn, dẫn đến nguy cơ bỏ sót tầng chứa dầu thương mại khi áp dụng các phương pháp minh giải địa vật lý giếng khoan (ĐVLGK) truyền thống.

+----------------------------------------------------------------------------------------------------+
|                                 BỐI CẢNH ĐỊA CHẤT & KỸ THUẬT LÔ 16-1                               |
|                                                                                                    |
|  [Đá móng nứt nẻ Pre-Cenozoic] -------------> Đóng góp ~80% sản lượng (20 triệu tấn/năm)           |
|  [Trầm tích Mioxen & Oligoxen] -------------> Đối tượng khai thác thứ 2 (Đặc tính thấm chứa tốt)   |
|                                                                                                    |
|                                         THÁCH THỨC CỐT LÕI:                                        |
|                     Vỉa chứa dầu điện trở thấp (Rt = 1.5 - 3.0 ohm.m)                              |
|                          ==> Nguy cơ bỏ sót tầng sản phẩm thương mại                               |
+----------------------------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định là: Các mô hình bão hòa nước điện trở cổ điển (Archie, Simandoux, Indonesian, Dual Water, Waxman-Smits) hoàn toàn phụ thuộc vào đường cong điện trở suất sâu ($R_t$). Khi $R_t$ bị suy giảm nghiêm trọng do ảnh hưởng kết hợp của đới ngấm dung dịch khoan gốc nước sâu và sự hiện diện của khoáng vật dẫn điện/sét phân tán, các mô hình này tính toán độ bão hòa nước ($S_w$) cao bất thường ($S_w > 60 - 70%$), mâu thuẫn trực tiếp với lưu lượng khai thác dầu khan thực tế (lên đến $55.000$ thùng/ngày đêm tại mỏ TGT).

Nghiên cứu tập trung giải quyết 3 câu hỏi và 3 giả thuyết khoa học định lượng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế vật lý - thạch học và kỹ thuật giếng khoan nào là nguyên nhân trực tiếp dẫn đến hiện tượng điện trở suất thấp tại các vỉa chứa cát kết mỏ TGT?
    • Giả thuyết 1 (H1): Hiện tượng $R_t$ thấp xuất phát từ tác động đồng thời của đới ngấm dung dịch khoan sâu vượt quá tầm đo của thiết bị ĐVLGK truyền thống và sự phân bố của sét phân tán có dung lượng trao đổi cation (CEC) cao kết hợp khoáng vật pyrit ($\text{FeS}_2$).
  2. Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xác định chính xác độ bão hòa nước ($S_w$), ranh giới nước tự do (FWL) và chiều cao cột dầu ($h$) mà không phụ thuộc vào đường cong điện trở suất bị biến dạng?
    • Giả thuyết 2 (H2): Tích hợp dữ liệu áp suất mao dẫn ($P_c$) chuẩn hóa qua hàm $J$ (Leverett $J$-function) kết hợp với đường phân bố áp suất dư ($EP$) từ thiết bị đo áp suất thành hệ điểm (MDT/RCI) sẽ xác định chính xác $S_w$ và FWL độc lập với $R_t$.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mối quan hệ tương quan rỗng - thấm biến đổi như thế nào theo các tướng trầm tích và ảnh hưởng ra sao đến khả năng ngậm nước ($f_w$, BS&W) khi mở vỉa?
    • Giả thuyết 3 (H3): Việc phân loại đá chứa thành các đơn vị thủy lực theo tướng trầm tích sông ngòi và đầm hồ sẽ chuẩn hóa tương quan độ rỗng ($\Phi$) - độ thấm ($K$), cho phép dự báo chính xác hàm lượng nước khai thác.

Phạm vi nghiên cứu bao quát toàn bộ cấu trúc mỏ TGT thuộc Lô 16-1 (cách Vũng Tàu 100 km về phía Đông Nam, độ sâu mực nước biển $40 - 50\text{ m}$), phân tích chi tiết lát cắt địa tầng Mioxen dưới (phân tập ILBH5.2, dày $370 - 490\text{ m}$) và Oligoxen trên (tập C, dày $440 - 500\text{ m}$). Cơ sở dữ liệu bao gồm 8 giếng khoan thăm dò/thẩm lượng (từ TGT-1X đến TGT-10XST1), 3 giếng khai thác phát triển, 31 giếng khai thác đang vận hành từ 3 giàn đầu giếng, dữ liệu đo áp suất thành hệ MDT từ 8 giếng và hệ thống mẫu lõi liên tục (RCA/SCAL) từ 6 giếng khoan.

Literature Review và Positioning

Các nghiên cứu đánh giá vỉa chứa shaly-sand và LRLC trên thế giới đã trải qua quá trình tiến hóa từ các mô hình bán thực nghiệm dựa trên điện trở đến các tiếp cận vi mô thạch học và thủy động lực học:

+----------------------------------------------------------------------------------------------------+
|                            TIẾN TRÌNH CÁC MÔ HÌNH ĐÁNH GIÁ ĐỘ BÃO HÒA NƯỚC                         |
|                                                                                                    |
|  [1942] Mô hình Archie                                                                             |
|         - Áp dụng cho cát kết sạch, đồng nhất                                                      |
|         - Không xét đến độ dẫn điện của sét                                                        |
|           |                                                                                        |
|           v                                                                                        |
|  [1963 - 1971] Mô hình Simandoux & Indonesian (Poupon - Leveaux)                                   |
|         - Bổ sung thể tích sét thể tích (Vsh)                                                      |
|         - Chưa tính đến khả năng dẫn điện màng ion của sét                                         |
|           |                                                                                        |
|           v                                                                                        |
|  [1968 - 1984] Mô hình Waxman-Smits & Dual Water (Clavier et al.)                                  |
|         - Tích hợp dung lượng trao đổi cation (CEC) và nước liên kết sét (CBW)                     |
|         - Vẫn phụ thuộc cốt lõi vào điện trở suất thực đo (Rt)                                     |
|           |                                                                                        |
|           v                                                                                        |
|  [1993 - Nay] Hướng tiếp cận tích hợp Áp suất mao dẫn (Pc) & Thủy động lực học                     |
|         - Leverett J-function, Cuddy, Skelt-Harrison                                              |
|         - Độc lập với Rt, hiệu chỉnh trực tiếp từ áp suất dư MDT/RCI và mẫu lõi                    |
+----------------------------------------------------------------------------------------------------+

Trong nước, các công trình của Ths. Nguyễn Phương Thủy (2010) tại Lô 01-02 bể Cửu Long và PGS.TS. Lê Hải An cùng cộng sự (2012) tại mỏ Rồng và Sư Tử Đen đã bước đầu chỉ ra nguyên nhân gây $R_t$ thấp do nước liên kết cao và khoáng vật thứ sinh ($\text{FeS}_2, \text{Fe}_2\text{O}_3$), song vẫn phụ thuộc vào việc hiệu chỉnh đường cong điện trở. Nghiên cứu của TS. Lê Trung Tâm và TS. Cù Minh Hoàng (2015) trên tạp chí Dầu khí đề xuất sử dụng tỷ số sóng âm $V_p/V_s$ kết hợp hệ số Poisson để nhận diện hydrocarbon trong trầm tích Turbidite, tuy nhiên chưa định lượng hóa được độ bão hòa nước phục vụ tính toán trữ lượng tại chỗ.

Trên trường quốc tế, hai trường phái nghiên cứu đối lập tồn tại rõ rệt:

  • Trường phái địa vật lý kỹ thuật cao: Tiêu biểu là công trình của Carolina et al. (1996, SPE 36150) tại hồ Maracaibo (Venezuela) sử dụng thiết bị chụp ảnh thành giếng (FMI) để bóc tách các lớp cát - sét mỏng cấp milimet; Hamada et al. (2000, 2001, SPE 64406, SPE 70040) ứng dụng cộng hưởng từ hạt nhân (NMR) phân tách phổ thời gian phân rã $T_2$ để cô lập nước liên kết sét (CBW) và nước mao dẫn (BVI) khỏi nước tự do (FFI).
  • Trường phái vật lý thạch học - áp suất mao dẫn: Tiêu biểu là các nghiên cứu của Cuddy (1993) và Skelt-Harrison (1995) chứng minh rằng trong điều kiện vỉa bị xâm nhập sâu hoặc khoáng vật dẫn điện phân tán phức tạp, áp suất mao dẫn ($P_c$) kết hợp thế năng thủy tĩnh là công cụ xác thực nhất để tái lập phân bố chất lưu.

Luận án định vị nghiên cứu tại điểm giao thoa giữa trầm tích học tướng đá, áp suất mao dẫn phi điện trở và động lực học dòng chảy thực tế. Luận án vượt qua giới hạn của các công trình trước đây tại bể Cửu Long bằng cách loại bỏ sự phụ thuộc vào đường cong $R_t$ bị nhiễu, thiết lập quy trình tính toán $S_w$ từ hàm $J(S_w)$ hiệu chỉnh theo tướng trầm tích và kiểm chứng đối sánh trực tiếp với dữ liệu đo mặt cắt dòng (PLT) và khai thác thực tế.

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

+----------------------------------------------------------------------------------------------------+
|                                KHUNG PHÂN TÍCH TỔNG THỂ CỦA LUẬN ÁN                                |
|                                                                                                    |
|    +-------------------------+     +-------------------------+     +-------------------------+     |
|    |   MÔI TRƯỜNG LẮNG ĐỌNG  |     |   ĐẶC TRƯNG ĐỚI NGẤM    |     |   KHOÁNG VẬT HỌC (XRD)  |     |
|    |   (Sông ngòi & Đầm hồ)  |     |   (LWD vs Wireline Rt)  |     |   (Sét CEC, Pyrit FeS2) |     |
|    +------------+------------+     +------------+------------+     +------------+------------+     |
|                 |                               |                               |                  |
|                 +-------------------------------+-------------------------------+                  |
|                                                 |                                                  |
|                                                 v                                                  |
|                      +------------------------------------------------------+                      |
|                      |  PHÂN LOẠI ĐÁ CHỨA (5 Hydraulic Flow Units theo Ka)  |                      |
|                      |        RQI = 0.0314 * sqrt(K/Phi_e)                 |                      |
|                      +--------------------------+---------------------------+                      |
|                                                 |                                                  |
|                                                 v                                                  |
|                      +------------------------------------------------------+                      |
|                      |   MÔ HÌNH THỦY TĨNH - ÁP SUẤT DƯ (Excess Pressure)   |                      |
|                      |      EP = (rho_w - rho_o) * g * h ==> Xác định FWL   |                      |
|                      +--------------------------+---------------------------+                      |
|                                                 |                                                  |
|                                                 v                                                  |
|                      +------------------------------------------------------+                      |
|                      |    TÍNH TOÁN BÃO HÒA NƯỚC PHI ĐIỆN TRỞ (Hàm J-Pc)    |                      |
|                      |       J = (Pc / (IFT * cos(theta))) * sqrt(K/Phi)    |                      |
|                      |       Sw = a * [J]^(-b)                             |                      |
|                      +--------------------------+---------------------------+                      |
|                                                 |                                                  |
|                                                 v                                                  |
|                      +------------------------------------------------------+                      |
|                      |   KIỂM CHỨNG ĐỘNG LỰC DÒNG (Dự báo fw & Đo dòng PLT) |                      |
|                      |         fw = 1 / [1 + (Kro/Krw)*(mu_w/mu_o)]         |                      |
|                      +------------------------------------------------------+                      |
+----------------------------------------------------------------------------------------------------+

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

Luận án mở rộng và kiểm chứng giới hạn áp dụng của lý thuyết điện trở Archie (1942) và lý thuyết áp suất mao dẫn Leverett (1941) trong môi trường trầm tích lục địa xen kẹp phức tạp.

  1. Mở rộng lý thuyết điện trở - bão hòa nước: Luận án chứng minh rằng định luật Archie ($S_w = (a\cdot R_w / (\Phi^m \cdot R_t))^{1/n}$) và các mô hình mở rộng Waxman-Smits hoàn toàn mất hiệu lực khi vỉa chứa cát kết hạt mịn có độ thấm cao bị ngấm sâu bởi dung dịch khoan gốc nước ($R_{mf} \ll R_w$), kết hợp với cấu trúc dẫn điện vi mô của màng nước bao quanh khoáng vật sét illit/smectit và pyrit ($\text{FeS}_2$).
  2. Mô hình hóa áp suất mao dẫn - thế năng thủy động: Luận án phát triển mệnh đề lý thuyết: Độ bão hòa nước trong vỉa chứa phân lớp mỏng điện trở thấp là hàm số đơn trị của chiều cao cột dầu ($h$), sức căng bề mặt chất lưu ($\text{IFT}$), góc thấm ướt ($\theta$), chỉ số chất lượng đá chứa ($RQI$) và áp suất dư ($EP$), hoàn toàn độc lập với điện trở suất đo được.

Mối quan hệ toán học được thiết lập qua hệ phương trình chuẩn hóa:

$$\Delta P = EP = P_{formation} - P_{hydrostatic_water} = (\rho_w - \rho_o) \cdot g \cdot h$$

$$J(S_w) = \frac{P_c}{\sigma \cdot \cos\theta} \sqrt{\frac{K}{\Phi}}$$

$$S_w = a \cdot [J(S_w)]^{-b}$$

Trong đó, chiều cao cột dầu ($h$) được xác định chính xác từ ranh giới nước tự do (FWL) tại vị trí $EP = 0\text{ psi}$.

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

Khung phân tích của luận án tích hợp bốn trục lý thuyết nền tảng:

  • Lý thuyết tướng trầm tích học (Sedimentary Facies): Phân định không gian lắng đọng thành 2 môi trường chủ đạo: môi trường sông ngòi (lòng sông - Channel, vỡ bờ - Crevasse Splay, bãi tràn - Sheetflood) và môi trường đầm hồ (vùng ngập nước - Overbank, đất bùn - Soil, dòng trọng lực - Gravity Flow).
  • Lý thuyết đơn vị dòng thủy lực (Hydraulic Flow Units): Phân chia đá chứa thành 5 cấp chất lượng dựa trên chỉ số chất lượng đá chứa ($RQI = 0{,}0314 \sqrt{K/\Phi_e}$) và chỉ số biến thể $RQI_m \cdot h$.
  • Lý thuyết áp suất mao dẫn đa pha (Capillary Pressure Theory): Sử dụng hàm số vô thứ tích $J$-Leverett để quy đổi dữ liệu đo mẫu lõi phòng thí nghiệm (Mercury Injection Capillary Pressure - MICP và Centrifuge) về điều kiện nhiệt áp vỉa thực tế với sức căng bề mặt $\text{IFT}_{res} = 25 - 30\text{ dyn/cm}$.
  • Lý thuyết dòng chảy đa pha Buckley-Leverett: Dự báo phân đoạn dòng chảy nước ($f_w$) và tỷ lệ nước ngậm ($BS&W$) tại miệng giếng dựa trên đường cong độ thấm tương đối ($K_{ro}, K_{rw}$) và tỷ số độ nhớt ($\mu_o/\mu_w$).

Điều kiện biên áp dụng: Bể trầm tích rift nội lục Kainozoi khu vực Đông Nam Á, nhiệt độ vỉa $85 - 135^\circ\text{C}$, áp suất vỉa $3.500 - 5.500\text{ psi}$, độ mặn nước vỉa $15.000 - 35.000\text{ ppm NaCl}$.

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

+----------------------------------------------------------------------------------------------------+
|                                    QUY TRÌNH NGHIÊN CỨU TÍCH HỢP                                   |
|                                                                                                    |
|  [Đo ĐVLGK & Đo trong khi khoan]     [Lấy mẫu lõi & Phân tích SCAL]     [Đo áp suất thành hệ MDT]  |
|       - LWD (Rsand) vs Wireline             - Phân tích lát mỏng, XRD          - Xác định Grad dầu/nước|
|       - FMI/STAR giải đoán góc vỉa          - Đo Pc, IFT, CEC, Kro/Krw         - Xác định FWL & EP     |
|                   |                                 |                                 |            |
|                   +---------------------------------+---------------------------------+            |
|                                                     |                                              |
|                                                     v                                              |
|                                  [CHUẨN HÓA & PHÂN NHÓM TƯỚNG ĐÁ]                                  |
|                                       - 5 Cấp độ thấm (K)                                          |
|                                       - Tương quan Phi - K sông/hồ                                 |
|                                                     |                                              |
|                                                     v                                              |
|                                  [MÔ HÌNH HÓA BÃO HÒA NƯỚC (HÀM J)]                                |
|                                       - Chuyển đổi Pc sang ĐK vỉa                                  |
|                                       - Xây dựng đường cong Sw - h                                 |
|                                                     |                                              |
|                                                     v                                              |
|                                  [KIỂM CHỨNG SẢN XUẤT & ĐO DÒNG PLT]                               |
|                                       - Đối chiếu lưu lượng Qo, Qw                                 |
|                                       - Hiệu chỉnh mô hình thủy động                               |
+----------------------------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism) với thiết kế định lượng thực nghiệm đa tầng. Luận án kết hợp giữa đo đạc vật lý trường sâu (In-situ Wireline/LWD), kiểm tra mẫu lõi phá hủy trong phòng thí nghiệm (Laboratory Core Testing) và thử nghiệm thủy động học động lực giếng khoan (Dynamic Well Testing).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm 4 bước độc lập và kiểm chứng chéo (Triangulation):

  1. Minh giải thạch học và khoáng vật học: Phân tích 120 lát mỏng thạch học và đo nhiễu xạ tia X (XRD) trên các mẫu lõi cắt qua tầng Mioxen hạ (ILBH5.2) và Oligoxen C tại 6 giếng thăm dò. Xác định tỷ lệ khoáng vật sét (illit, smectit, kaolinit, clorit), dung lượng trao đổi cation (CEC, $\text{meq}/100\text{g}$) và sự phân bố của khoáng vật dẫn điện pyrit ($\text{FeS}_2$).

  2. Xác định cơ chế suy giảm điện trở: Đo đạc đối sánh tức thời giữa hệ thống đo địa vật lý trong khi khoan (LWD) và đo dây cáp sau khi khoan (Wireline logging). Thiết lập quan hệ hàm mũ giữa điện trở biểu kiến sau khoan và điện trở thực đo trong khi khoan:

    $$R_{\text{Wireline}} = 0{,}8887 \cdot (R_{\text{LWD}})^{0{,}942}$$

    Minh chứng rằng đới ngấm làm suy giảm giá trị $R_t$ từ $30 - 50%$ so với trạng thái ban đầu của thành hệ.

  3. Phân nhóm chất lượng đá chứa (Reservoir Rock Typing): Dữ liệu rỗng - thấm từ hàng trăm mẫu lõi thông thường (RCA) được phân tách thành 5 nhóm độ thấm ($K_1 < 10\text{ mD}$; $K_2 = 10 - 50\text{ mD}$; $K_3 = 50 - 200\text{ mD}$; $K_4 = 200 - 1000\text{ mD}$; $K_5 > 1000\text{ mD}$).

  4. Mô hình hóa áp suất mao dẫn và bão hòa nước: Tiến hành đo áp suất mao dẫn bằng phương pháp ép thủy ngân (MICP) và ly tâm trên các mẫu đại diện cho 5 nhóm đá chứa. Quy đổi $P_c$ phòng thí nghiệm về điều kiện vỉa dựa trên tương quan sức căng bề mặt $\text{IFT}$ và góc thấm ướt $\theta$. Xây dựng hàm $J(S_w)$ riêng biệt cho từng nhóm tướng trầm tích.

Nhóm đá chứa (Rock Class) Dải độ thấm tuyệt đối ($K_a$, mD) Môi trường trầm tích đặc trưng Độ bão hòa nước bất động ($S_{wirr}$, %)
Nhóm 1 (K1) $< 10$ Đầm hồ, đất bùn (Soil), bãi tràn mỏng $55 - 75$
Nhóm 2 (K2) $10 - 50$ Cát bãi tràn ngập nước (Overbank/SF) $40 - 55$
Nhóm 3 (K3) $50 - 200$ Cát vỡ bờ (Crevasse Splay), dòng trọng lực $30 - 40$
Nhóm 4 (K4) $200 - 1000$ Lòng sông thứ cấp, cát lòng hồ năng lượng cao $20 - 30$
Nhóm 5 (K5) $> 1000$ Cát lòng sông chính (Channel Fill) $12 - 20$

Data và phân tích

Dữ liệu phân tích áp suất thành hệ điểm (MDT/RCI) được xử lý theo phương pháp gradient áp suất cục bộ và chuyển đổi sang đồ thị áp suất dư ($EP$). Độ tin cậy của mô hình được bảo đảm thông qua việc đối chiếu với dữ liệu đo lưu lượng của 2 giếng khai thác có thả thiết bị đo mặt cắt dòng (PLT) và phân tích mẫu chất lưu đáy giếng (BHS). Các phần mềm chuyên dụng như Techlog (Schlumberger), Interactive Petrophysics (IP) và Eclipse được sử dụng để tích hợp và chạy mô phỏng.

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

+----------------------------------------------------------------------------------------------------+
|                         CÁC PHÁT HIỆN ĐỘT PHÁ VÀ GIÁ TRỊ ĐỊNH LƯỢNG                                |
|                                                                                                    |
|  [PHÁT HIỆN 1: Nguyên nhân cốt lõi của Rt thấp (1.5 - 3.0 ohm.m)]                                  |
|  ==> Đới ngấm sâu vượt tầm đo ĐVLGK + Sét phân tán CEC cao + Khoáng vật Pyrit (FeS2: 2 - 8%)       |
|                                                                                                    |
|  [PHÁT HIỆN 2: Sai lệch mô hình Archie/Simandoux truyền thống]                                     |
|  ==> Đánh giá sai Sw = 60 - 75% (loại bỏ tầng sản phẩm) | Thực tế Sw = 22 - 38% (khai thác dầu khan)|
|                                                                                                    |
|  [PHÁT HIỆN 3: Bất đối xứng tiềm năng Bắc - Nam mỏ TGT]                                            |
|  ==> Phía Bắc (Khối H1, H2): Cát lòng sông, Phi > 24%, K > 1500 mD, Cột dầu h lớn                   |
|  ==> Phía Nam (Khối H4, H5): Cát mỏng xen kẹp sét, Phi < 16%, K < 50 mD, Cột dầu h hẹp             |
+----------------------------------------------------------------------------------------------------+

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

  1. Làm sáng tỏ bản chất kép gây điện trở suất thấp ($1{,}5 - 3{,}0\ \Omega\cdot\text{m}$): Luận án chứng minh hiện tượng điện trở thấp tại mỏ TGT không phản ánh sự bão hòa nước của vỉa, mà là kết quả cộng hưởng của hai yếu tố:
    • Yếu tố kỹ thuật: Độ rỗng cao ($\Phi = 20 - 28%$) và độ thấm lớn ($K > 500 - 2.000\text{ mD}$) của tập ILBH5.2 tạo điều kiện cho dung dịch khoan gốc nước xâm nhập rất sâu xung quanh thành giếng, vượt qua chiều sâu nghiên cứu của thiết bị đo điện trở sâu (LLD, RD).
    • Yếu tố địa chất: Sự hiện diện của sét phân tán dạng illit-smectit xen kẽ trong lỗ rỗng làm tăng dung lượng trao đổi cation (CEC), kết hợp với khoáng vật pyrit dẫn điện ($\text{FeS}_2$, chiếm tỷ lệ $2 - 8%$ thể tích trong lát mỏng) tạo ra mạng lưới dẫn điện bề mặt liên tục.
  2. Khắc phục triệt để sai số tính toán độ bão hòa nước ($S_w$): Các phương pháp truyền thống tính ra $S_w$ trong khoảng $60 - 75%$, xếp tầng chứa vào diện nguy cơ ngập nước. Phương pháp hàm $J$ kết hợp áp suất dư ($EP$) xác định độ bão hòa nước thực tế chỉ đạt $22 - 38%$. Kết quả này được kiểm chứng chính xác bởi các giếng khai thác thực tế mở vỉa tại các tập này thu được dòng dầu không ngập nước với lưu lượng từ $2.000$ đến hơn $5.000$ thùng/ngày/giếng.
  3. Phát hiện quy luật biến đổi thấm chứa phân vùng Bắc - Nam: Tiềm năng thấm chứa suy giảm rõ rệt theo trục từ Bắc xuống Nam dọc theo 5 khối cấu trúc (H1 $\to$ H2 $\to$ H3 $\to$ H4 $\to$ H5):
    • Khu vực phía Bắc (Khối H1, H2): Tướng trầm tích lòng sông chiếm ưu thế, cát kết dày, độ rỗng trung bình $\Phi = 22 - 26%$, độ thấm $K = 500 - 2.500\text{ mD}$, chiều cao cột dầu ($h$) lớn, ranh giới OWC/FWL sâu.
    • Khu vực phía Nam (Khối H4, H5): Tướng trầm tích chuyển dần sang đầm hồ nước sâu và bãi tràn ven bờ, các vỉa cát mỏng dần, xen kẹp nhiều tập sét dày, độ rỗng giảm xuống $\Phi = 14 - 18%$, độ thấm giảm mạnh $K = 5 - 150\text{ mD}$, chiều cao cột dầu hẹp.
  4. Chuẩn hóa ranh giới nước tự do (FWL) và ranh giới dầu nước (OWC): Bằng việc dựng đường biểu diễn áp suất dư $EP$ theo độ sâu thẳng đứng (TVDss), luận án đã xác lập chính xác độ sâu FWL cho từng khối thủy lực độc lập, giải thích được hiện tượng chênh lệch áp suất vỉa giữa các khối lân cận bị ngăn cách bởi hệ thống đứt gãy thuận trượt bằng phương Tây Bắc - Đông Nam.
+----------------------------------------------------------------------------------------------------+
|                         MÔ HÌNH PHÂN HÓA ĐẶC TÍNH VỈA THEO CHIỀU BẮC - NAM                         |
|                                                                                                    |
|   [BẮC: Khối H1, H2] -------------> [TRUNG TÂM: Khối H3] -------------> [NAM: Khối H4, H5]         |
|   - Tướng lòng sông (Channel)       - Tướng hỗn hợp                      - Tướng hồ/bãi tràn        |
|   - Phi = 22 - 26%                  - Phi = 18 - 22%                     - Phi = 14 - 18%           |
|   - K = 500 - 2500 mD               - K = 100 - 500 mD                   - K = 5 - 150 mD           |
|   - Cột dầu h lớn, fw thấp          - Cột dầu trung bình                 - Cột dầu hẹp, fw tăng cao |
+----------------------------------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm vững chắc chứng minh rằng áp suất mao dẫn là công cụ tối thượng để đánh giá độ bão hòa trong các vỉa chứa dị thường điện trở, đặt nền móng cho việc chuẩn hóa mô hình petrophysics phi Archie tại thềm lục địa Nam Việt Nam.
  • Về mặt phương pháp luận: Thiết lập quy trình minh giải 5 bước tích hợp (LWD/Wireline $\to$ Rock Typing $\to$ MDT Pressure Gradient $\to$ Scaled $J$-Function $\to$ Fractional Flow Prediction) có khả năng chuyển giao và áp dụng trực tiếp cho các bể trầm tích lân cận (bể Nam Côn Sơn, Sông Hồng, Malay - Thổ Chu).
  • Về mặt thực tiễn và kinh tế mỏ: Tránh được việc bỏ sót các tập vỉa giàu dầu có $R_t$ thấp trong quá trình thử vỉa (DST); tối ưu hóa vị trí bắn mở vỉa khai thác, giúp công ty Hoàng Long - Hoàn Vũ JOC duy trì sản lượng khai thác ổn định $20.000 - 25.000$ thùng/ngày đêm tại mỏ TGT; hỗ trợ tính toán chính xác trữ lượng dầu tại chỗ (STOIIP), nâng cao hệ số thu hồi dầu (RF).

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án cũng ghi nhận một số giới hạn khoa học:

  1. Mật độ mẫu lõi đặc biệt (SCAL): Dữ liệu đo áp suất mao dẫn và độ thấm tương đối chỉ tập trung tại 6 giếng thăm dò, chưa phủ đều toàn bộ các khối cấu tạo phức tạp phía Nam (Khối H5).
  2. Độ phân giải của thiết bị đo: Đối với các tập cát sét phân lớp cực mỏng (dưới $5\text{ cm}$) ở tầng Oligoxen D sâu, các đường cong ĐVLGK thông thường vẫn chịu hiệu ứng trung bình hóa thể tích, đòi hỏi phải dựa nhiều vào tài liệu FMI/STAR.
  3. Mô phỏng tính dính ướt động: Nghiên cứu giả định góc dính ướt $\theta$ và sức căng bề mặt $\text{IFT}$ không đổi trong suốt giai đoạn khai thác sơ cấp, chưa mô phỏng biến đổi dính ướt cục bộ khi vỉa bị ngập nước bơm ép nhân tạo.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Ứng dụng vật lý đá số (Digital Rock Physics - DRP): Sử dụng công nghệ chụp cắt lớp vi tính $\mu\text{-CT}$ kết hợp mô phỏng mạng vi dòng (Pore Network Modeling / Lattice Boltzmann) trên mẫu đá mỏ TGT để định lượng hóa chính xác dòng chảy đa pha ở quy mô lỗ rỗng nano-micromet.
  • Tích hợp mạng nơ-ron nhân tạo (ANN / Deep Learning): Xây dựng thuật toán học máy dự báo liên tục đường cong $S_w$ và $K_a$ từ tổ hợp đường cong ĐVLGK tiêu chuẩn đã được huấn luyện bởi dữ liệu mẫu lõi và áp suất MDT.
  • Nghiên cứu động thái dính ướt trong quá trình bơm ép thứ cấp: Đánh giá sự biến đổi của hàm $J$ khi có sự tương tác giữa nước bơm ép đã xử lý hóa chất với bề mặt khoáng vật sét chứa pyrit.

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

  • Ảnh hưởng học thuật: Luận án là công trình nghiên cứu hệ thống và toàn diện đầu tiên về trầm tích điện trở thấp tại mỏ TGT, Lô 16-1 bể Cửu Long. Các công trình công bố từ luận án trên Tạp chí Dầu khí Việt Nam và Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất tạo lập nguồn tài liệu tham khảo cốt lõi cho các nghiên cứu sinh và chuyên gia petrophysics trong khu vực.
  • Tác động công nghiệp: Chuyển đổi trực tiếp quy trình đánh giá thành hệ tại Công ty Điều hành Chung Hoàng Long - Hoàn Vũ (HL-HV JOC), Vietsovpetro và các JOC khác đang vận hành tại bể Cửu Long (Cửu Long JOC, Thăng Long JOC). Quy trình này giúp gia tăng trữ lượng thương mại được công nhận tại Lô 16-1 thêm hàng triệu thùng dầu quy đổi.
  • Giá trị kinh tế - xã hội: Tối ưu hóa chi phí khoan và can thiệp giếng, loại trừ các ca thử vỉa lãng phí hoặc bắn mở nhầm tầng ngập nước, đóng góp trực tiếp vào nguồn thu ngân sách nhà nước thông qua việc gia tăng hiệu quả khai thác tài nguyên dầu khí quốc gia.

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

  • Nghiên cứu sinh & Nhà khoa học Địa chất Dầu khí: Tiếp cận phương pháp luận hoàn chỉnh về xử lý bài toán vỉa chứa LRLC; nắm bắt phương thức kết hợp giữa trầm tích tướng đá và động lực học chất lưu.
  • Kỹ sư Petrophysics & Địa chất mỏ (Geoscientists): Sở hữu bộ công cụ thực hành minh giải tham số vỉa ($S_w, K, h$) chuẩn xác từ dữ liệu ĐVLGK, MDT và mẫu lõi khi gặp các địa tầng có $R_t$ thấp dị thường.
  • Kỹ sư Công nghệ mỏ & Khai thác (Reservoir & Production Engineers): Sử dụng các dự báo chính xác về phân đoạn nước ($f_w$, BS&W) để thiết kế đường ống, trạm xử lý chất lưu bề mặt và tối ưu hóa khoảng mở vỉa khoan khai thác.
  • Các nhà quản lý tài nguyên (PetroVietnam/PVEP): Có cơ sở khoa học tin cậy để thẩm định trữ lượng dầu khí tại chỗ (STOIIP) và phê duyệt Kế hoạch Phát triển Mỏ (FDP), Kế hoạch Thu dọn Mỏ tối ưu.

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 lý thuyết nào? Trả lời: Luận án đã mở rộng lý thuyết áp suất mao dẫn Leverett (1941) và lý thuyết đơn vị dòng thủy lực (Amaefule et al., 1993) vào môi trường vỉa chứa cát sét phân lớp mỏng LRLC bể Cửu Long. Đóng góp độc đáo nhất là chứng minh $S_w$ có thể được định lượng chính xác hoàn toàn không cần dựa vào đường cong điện trở $R_t$, thông qua hàm tương quan $S_w = a \cdot [J(S_w)]^{-b}$ gắn liền với phân loại 5 cấp độ thấm theo tướng trầm tích và xác định chiều cao cột dầu $h$ từ điểm cân bằng áp suất dư $EP = 0\text{ psi}$.
  2. Tính đổi mới trong phương pháp luận so với các nghiên cứu quốc tế tương đương? Trả lời: So với nghiên cứu của Carolina et al. (1996) tại Venezuela (vốn chỉ dựa vào phân tích lát mỏng ảnh FMI) và Hamada et al. (2000, 2001) (dựa thuần túy vào đo phổ NMR đắt đỏ), luận án thiết lập một giải pháp toàn diện và kinh tế hơn: phối hợp thạch học XRD/CEC với phân tích đối sánh tức thời LWD vs Wireline, liên kết chặt chẽ áp suất dư MDT với áp suất mao dẫn mẫu lõi $P_c$, và đặc biệt là kiểm chứng khép kín với dữ liệu đo mặt cắt dòng khai thác thực tế (PLT).
  3. Phát hiện thực nghiệm gây bất ngờ nhất có số liệu chứng minh? Trả lời: Tại phân tập ILBH5.2 mỏ TGT, các khoảng vỉa có điện trở suất sâu biểu kiến cực thấp ($R_t = 1{,}8 - 2{,}5\ \Omega\cdot\text{m}$) - mức vốn được xem là chứa nước 100% theo mô hình Archie ($S_w > 70%$) - khi áp dụng mô hình hàm $J$ chuẩn hóa theo tướng lòng sông đã tính toán lại $S_w = 26 - 32%$. Khi tiến hành mở vỉa, giếng đã cho dòng dầu tự phun đạt trên $4.500$ thùng/ngày đêm với hàm lượng nước $BS&W = 0%$.
  4. Quy trình nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không? Trả lời: Có. Luận án cung cấp quy trình thuật toán chi tiết từng bước: (1) Hiệu chỉnh $R_t$ theo tương quan LWD-Wireline; (2) Phân nhóm mẫu lõi theo $RQI$ và cấp độ thấm; (3) Chuẩn hóa đường cong $P_c$ về hàm $J$; (4) Xác định FWL từ dữ liệu áp suất điểm MDT; (5) Tính phân bố $S_w(h)$ và tích hợp với đường cong độ thấm tương đối $K_{ro}/K_{rw}$ để tính $f_w$. Quy trình này có thể lập trình tự động hóa trên mọi phần mềm minh giải ĐVLGK chuẩn.
  5. Chương trình nghiên cứu 10 năm được phác thảo ra sao? Trả lời: Chương trình 10 năm định hình 3 trục: (a) Xây dựng cơ sở dữ liệu số hóa vật lý đá vi mô (Digital Rock Repository) cho toàn bộ các mỏ LRLC tại thềm lục địa Việt Nam; (b) Phát triển mô-đun AI/Machine Learning tự động phát hiện tầng LRLC thời gian thực trong khi khoan (Real-time LWD Geosteering); (c) Nghiên cứu cơ chế dịch chuyển chất lưu đa pha tăng cường thu hồi dầu (EOR) cho các vỉa chứa phân tầng mỏng ngậm sét cao.

Kết luận

Công trình nghiên cứu của tác giả Bùi Hữu Phước đã mang lại những đóng góp khoa học và thực tiễn mang tính then chốt:

  1. Làm rõ toàn diện cơ chế địa chất và kỹ thuật giếng khoan: Xác định chính xác đới ngấm dung dịch khoan sâu và sự tồn tại của sét phân tán giàu cation cùng khoáng vật dẫn điện pyrit ($\text{FeS}_2$) là hai nguyên nhân căn bản tạo nên hiện tượng điện trở suất thấp ($1{,}5 - 3{,}0\ \Omega\cdot\text{m}$) tại mỏ TGT, Lô 16-1 bể Cửu Long.
  2. Xây dựng hệ phương pháp đánh giá bão hòa nước phi điện trở: Đề xuất thành công quy trình xác định $S_w$ dựa trên hàm áp suất mao dẫn $J$-Leverett kết hợp chiều cao cột dầu $h$ từ dữ liệu áp suất dư MDT/RCI, giải phóng hoàn toàn sự phụ thuộc vào đường cong điện trở bị nhiễu.
  3. Thiết lập mô hình tương quan rỗng - thấm theo tướng trầm tích: Phân định 5 đơn vị dòng thủy lực gắn liền với môi trường sông ngòi và đầm hồ, cung cấp bộ thông số đầu vào chính xác cho mô hình thủy động lực học mỏ.
  4. Chứng minh quy luật phân hóa thấm chứa Bắc - Nam: Phân định tiềm năng vượt trội của khu vực phía Bắc (khối H1, H2) so với khu vực phía Nam (khối H4, H5), định hướng chiến lược bố trí giếng khoan khai thác tối ưu.
  5. Mở ra hướng nghiên cứu và phát triển mới cho ngành dầu khí Việt Nam: Phương pháp luận của luận án cung cấp lời giải chuẩn xác cho hàng loạt mỏ và cấu tạo chứa dầu LRLC tương tự tại thềm lục địa, góp phần gia tăng sản lượng và kéo dài tuổi thọ khai thác kinh tế của các mỏ dầu khí quốc gia.