Tổng quan về luận án

Trong bối cảnh tiêu thụ năng lượng toàn cầu dự báo đạt 886,3 triệu tỷ Btu vào năm 2050 theo báo cáo của Cơ quan Thông tin Năng lượng Hoa Kỳ (U.S. EIA, 2021), các giải pháp thu hồi dầu tăng cường (Enhanced Oil Recovery - EOR) tích hợp thu hồi, sử dụng và lưu trữ carbon (CCUS) giữ vai trò chiến lược. Luận án tiến sĩ kỹ thuật địa chất dầu khí với tiêu đề "Study on Enhanced Oil Recovery by CO2 Microbubbles Injection" do nghiên cứu sinh Lê Nguyễn Hải Nam thực hiện tại Khoa Kỹ thuật Tài nguyên Trái đất, Trường Sau đại học Kỹ thuật, Đại học Kyushu, Nhật Bản (Kyushu University, tháng 9/2022) dưới sự hướng dẫn của Giáo sư Yuichi Sugai, đã giải quyết rào cản kỹ thuật cốt lõi trong các dự án CO2-EOR bậc ba.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án nhắm đến là hiện tượng chênh lệch tỷ trọng (gravity segregation), phân luồng khí sớm (viscous fingering/gas channeling) và hiệu suất quét (sweep efficiency) kém khi bơm ép CO2 dạng khí tự do vào các vỉa chứa có tính dị hướng và bất đồng nhất cao (heterogeneous porous media). Dù bọt quy ước (conventional foam) từng được xem là giải pháp điều khiển linh độ (mobility control), sự kém bền nhiệt động học trong môi trường vỉa khắc nghiệt khiến hiệu quả suy giảm nghiêm trọng. Luận án đặt nền tảng tiên phong bằng việc nghiên cứu ứng dụng vi bọt khí CO2 dưới dạng các thể hạt khí keo (Colloidal Gas Aphrons - CGAs) nhằm tối ưu hóa tính chọn lọc nút kín và nâng cao hiệu suất thu hồi dầu.

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

  1. Câu hỏi nghiên cứu 1 (RQ1): Thành phần hóa học (chất hoạt động bề mặt, biopolymer, độ mặn) chi phối như thế nào đến độ bền pha và phân bố kích thước hạt (Bubble Size Distribution - BSD) của hệ vi bọt CO2?
    • Giả thuyết 1 (H1): Sự gia tăng nồng độ biopolymer Xanthan Gum và chất hoạt động bề mặt Sodium Dodecyl Sulfate (SDS) tạo cấu trúc màng đa lớp bảo vệ, giúp thu hẹp kích thước vi bọt và giảm thiểu tốc độ thoát nước theo quy luật hàm Weibull.
  2. Câu hỏi nghiên cứu 2 (RQ2): Quy luật lưu biến học và cơ chế nút kín chọn lọc của vi bọt CO2 trong môi trường xốp đồng nhất và bất đồng nhất diễn ra như thế nào?
    • Giả thuyết 2 (H2): Vi bọt CO2 tuân theo mô hình lưu biến phi Newton Power-law, kích hoạt hiệu ứng Jamin (Jamin effect) để chặn chọn lọc các kênh dẫn có độ thấm cao (high-permeability zones) và chuyển hướng dòng quét sang đới độ thấm thấp.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mức độ gia tăng hệ số thu hồi dầu (Recovery Factor - RF) khi bơm ép vi bọt CO2 trong mô hình vỉa vật lý bất đồng nhất đạt bao nhiêu so với kỹ thuật ngập nước và bơm dung dịch cơ sở?
    • Giả thuyết 3 (H3): Bơm vi bọt CO2 sau giai đoạn ngập nước thứ cấp sẽ gia tăng lượng dầu thu hồi trên 20% lượng dầu nguyên khai tại chỗ (OOIP) và đạt tổng hệ số thu hồi trên 85% trong vỉa bất đồng nhất.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết cấu trúc Aphron đa lớp (Sebba, 1987), cơ chế cản mao dẫn Jamin (Wright, 1933), động học phân rã bọt (Yan et al., 2005) và phương trình khuếch tán Young-Laplace. Phạm vi thực nghiệm được kiểm soát nghiêm ngặt trên các mô hình mô phỏng vỉa cát nén (sandpack models) đơn lõi và song song với độ thấm biến thiên từ 0,08 Darcy đến 4,0 Darcy, tại điều kiện nhiệt độ vỉa 45°C và nồng độ muối khoáng lên đến 20.000 mg/L NaCl.

Literature Review và Positioning

Tổng quan y văn chỉ ra ba dòng nghiên cứu chính trong kỹ thuật kiểm soát độ linh động CO2:

  1. Dòng nghiên cứu bọt CO2 truyền thống (Conventional CO2 Foam): Các công trình của Lake (1989), Du et al. (2007) và Yang et al. (2017) khẳng định bọt khí giúp tăng độ nhớt biểu kiến của pha khí, giảm tỷ số linh độ (mobility ratio). Tuy nhiên, Razavi et al. (2017) chứng minh màng mỏng lamellae của bọt quy ước rất dễ vỡ dưới tác động của lực trượt và dầu béo, dẫn đến mất tính ổn định nhanh chóng.
  2. Dòng nghiên cứu thể khí keo CGAs ứng dụng trong dung dịch khoan và xử lý môi trường: Sebba (1987) đặt nền tảng định nghĩa CGAs bao gồm lõi khí hình cầu đường kính 10–100 μm được bao bọc bởi màng mỏng surfactant ba lớp. Các tác giả Hashim et al. (1998, 2012), Waters et al. (2008), Shivhare et al. (2014) và Pasdar et al. (2019) đã mở rộng ứng dụng CGAs để ngăn xâm nhập dung dịch khoan và kiểm soát nứt nẻ. Bjorndalen et al. (2008, 2014) khảo sát độ bền CGAs dưới áp suất lên đến 13,8 MPa nhưng ghi nhận hiện tượng bất ổn định khi nhiệt độ vượt ngưỡng 50–75°C.
  3. Dòng nghiên cứu vi bọt trong nâng cao thu hồi dầu (Microbubble EOR): Telmadarreie et al. (2016) và Natawijaya et al. (2020) chứng minh tiềm năng cải thiện hiệu suất quét dầu nặng của vi bọt. Shi et al. (2016, 2017) sử dụng mô hình vi mô (micromodel) để chỉ ra khả năng nút kín tạm thời các đới thấm cao.

Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm đối nghịch về ảnh hưởng của chất hoạt động bề mặt lên phân bố kích thước vi bọt. Nhóm nghiên cứu của Xu et al. (2009a) và Bjorndalen et al. (2008) cho rằng việc tăng nồng độ surfactant dưới điểm nồng độ micelle tới hạn (Critical Micelle Concentration - CMC) làm giảm sức căng bề mặt, tạo điều kiện chia cắt bọt mịn hơn và làm giảm đường kính vi bọt. Ngược lại, Arabloo et al. (2014) và Pasdar et al. (2018c) lập luận rằng nồng độ surfactant tăng có thể thúc đẩy sự gia tăng kích thước hạt CGAs do tương tác tạo màng dày hơn.

Luận án của Lê Nguyễn Hải Nam định vị chính xác khoảng trống này bằng việc thực hiện các phân tích thống kê định lượng nghiêm ngặt và kiểm định độ phù hợp phân bố xác suất (goodness of fit). So sánh trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Zhu et al. (2020) chỉ sử dụng phân phối Gaussian (Normal distribution) đơn giản để khớp dữ liệu bọt khoan mà không kiểm định tính chuẩn.
  • Nghiên cứu của Alam et al. (2017) áp dụng phân phối Log-normal cho tế bào tuyển nổi điện hóa. Luận án đã chứng minh vượt trội rằng phân bố đường kính vi bọt CO2 tuân theo hàm phân phối Weibull với mức ý nghĩa thống kê vượt trội, giải quyết triệt để các bất định thực nghiệm trong mô tả kích thước vi bọt phục vụ EOR.

Đó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 thuyết cấu trúc Aphron của Sebba (1987) bằng cách làm sáng tỏ cơ chế gia cố màng bảo vệ nhờ sự kết hợp giữa biopolymer Xanthan Gum (XG) và anionic surfactant Sodium Dodecyl Sulfate (SDS). Nền tảng lý thuyết màng ba lớp (inner monolayer, intermediate aqueous shell, outer bilayer) được bổ sung cơ chế tăng độ nhớt cấu trúc cục bộ và ức chế sự khuếch tán khí CO2 từ nhân bọt ra pha lỏng liên tục:

$$\text{"Microbubbles – Colloidal Gas Aphrons (CGAs) have been reported as unique bubbles with micro-size (10-100 }\mu\text{m) consisting of a multilayer shell of surfactant molecules and a spherical gaseous core."}$$

Nghiên cứu thách thức các giả định trước đây về sự phân rã tự do của bọt khí bằng việc tích hợp phương trình chênh lệch áp suất mao dẫn Young-Laplace trong hiện tượng Ostwald ripening (Disproportionation):

$$\Delta P = P_1 - P_2 = 2\gamma \left(\frac{1}{r_1} - \frac{1}{r_2}\right)$$

Nhờ sự hiện diện của polymer XG ở nồng độ tối ưu, quá trình khuếch tán khí từ các bọt nhỏ ($r_1$) sang các bọt lớn ($r_2$) bị cản trở đáng kể bởi lực cản nhớt, làm giảm tốc độ suy biến kích thước theo thời gian.

Mô hình động học thoát nước tĩnh của Yan et al. (2005) được chứng thực và phát triển thông qua phương trình:

$$V_t = V_F \frac{t^n}{T_{1/2}^n + t^n}$$

Trong đó, hằng số tốc độ thoát nước vi phân $K$ phản ánh trực tiếp độ bền cấu trúc của màng vi bọt. Dữ liệu thực nghiệm của luận án cung cấp bằng chứng định lượng mạnh mẽ: khi nồng độ XG tăng từ 0 lên 5 g/L (tại 3 g/L SDS), hằng số $K$ giảm đột ngột từ $5,2 \times 10^{-5}\text{ mL}^{-1}\text{min}^{-1}$ (giảm 78 lần), thiết lập một bước chuyển dịch mô hình (paradigm shift) từ việc xem vi bọt như một hệ phân tán kém bền sang một hệ lưu chất phức hợp có kiểm soát động học phân rã.

+-------------------------------------------------------------------------------+
|                       KHUNG KHÁI NIỆM HỆ VI BỌT CO2 CGA                       |
+-------------------------------------------------------------------------------+
|                                                                               |
|   +-----------------------------------------------------------------------+   |
|   |                        1. HỆ HÓA CHẤT ĐẦU VÀO                         |   |
|   |   - Biopolymer: Xanthan Gum (XG: 1 - 5 g/L)                           |   |
|   |   - Surfactant: Sodium Dodecyl Sulfate (SDS: 1 - 3 g/L)               |   |
|   |   - Chất điện ly: NaCl (0 - 20 g/L)                                   |   |
|   |   - Pha khí: CO2 tinh khiết 99.9%                                     |   |
|   +-----------------------------------+-----------------------------------+   |
|                                       |                                       |
|                                       v (Cơ chế đồng hóa 8000 rpm, 4 min)     |
|   +-----------------------------------------------------------------------+   |
|   |                 2. CẤU TRÚC VI BỌT VÀ ĐẶC TÍNH ĐỘNG HỌC               |   |
|   |   - Cấu trúc: Nhân khí CO2 + Màng ba lớp polymer/surfactant           |   |
|   |   - Kích thước: Phân bố Weibull (Davg = 47.37 - 63.37 micromet)       |   |
|   |   - Ức chế Ostwald Ripening & Ức chế thoát nước (K giảm 78 lần)       |   |
|   +-----------------------------------+-----------------------------------+   |
|                                       |                                       |
|                                       v (Bơm ép vào môi trường xốp)           |
|   +-----------------------------------------------------------------------+   |
|   |                  3. LƯU BIẾN VÀ ĐỘNG LỰC HỌC DÒNG CHẢY                |   |
|   |   - Dòng chảy Power-law (Ức chế độ nhớt khi cắt / Shear-thinning)     |   |
|   |   - Hiệu ứng Jamin: Nút kín chọn lọc đới thấm cao (High-k zones)       |   |
|   |   - Cân bằng dòng phân đoạn (Fractional Flow 50%:50%)                 |   |
|   +-----------------------------------+-----------------------------------+   |
|                                       |                                       |
|                                       v (Hiệu quả thu hồi dầu)                |
|   +-----------------------------------------------------------------------+   |
|   |                4. TỐI ƯU HÓA EOR VÀ ĐỊA LƯU TRỮ CARBON                |   |
|   |   - Gia tăng hệ số thu hồi dầu đơn lõi: +23.6% sau ngập nước          |   |
|   |   - Thu hồi dầu vỉa bất đồng nhất (1:4): 86.9% OOIP                   |   |
|   |   - Thu hồi dầu vỉa bất đồng nhất (1:2): 93.0% OOIP                   |   |
|   |   - Tăng cường cô lập và lưu trữ địa chất CO2 an toàn                 |   |
|   +-----------------------------------------------------------------------+   |
|                                                                               |
+-------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên ngành ba hệ lý thuyết:

  1. Lý thuyết Hóa lý bề mặt và Hệ phân tán keo: Giải thích tương tác tĩnh điện màng đôi (electric double layer) và lực đẩy tĩnh điện giữa các ion sulfate của SDS.
  2. Lý thuyết Thủy lực học môi trường xốp và Cơ học lưu chất phi Newton: Mô tả ứng xử dòng chảy của hệ vi bọt tuân theo mô hình hàm số mũ (Power-law fluid model): $$\tau = K_p \dot{\gamma}^{n_p}$$ với các chỉ số nhớt $K_p$ và chỉ số dòng chảy $n_p < 1$ (tính chất giả dẻo - shear-thinning).
  3. Lý thuyết Điều khiển độ linh động trong Kỹ thuật mỏ dầu khí (Mobility Control Theory): Mô tả khả năng tái định tuyến dòng quét (flow diversion) dựa trên biến thiên áp suất mao dẫn qua cổ lỗ xốp (pore throats).

Điều kiện biên (Boundary Conditions): Khung phân tích được xác lập trong giới hạn nhiệt độ vỉa chứa $\le 45^\circ\text{C}$, độ mặn dung dịch không vượt quá 20 g/L NaCl, và tỷ lệ đường kính vi bọt trên kích thước hạt cốt sỏi/cổ lỗ xốp ($D_{\text{bubble}}/D_{\text{pore}}$) nằm trong khoảng tối ưu để kích hoạt lực cản Jamin mà không gây tắc nghẽn vĩnh viễn hệ thống vỉa.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng (positivism) với thiết kế thực nghiệm định lượng đa cấp độ (multi-level quantitative experimental design):

  • Cấp độ vi mô (Micro-scale level): Quan sát hình thái học, cấu trúc bọt và cơ chế va chạm, kết tụ dưới kính hiển vi quang học truyền qua (transmitted-light microscope) kết hợp cảm biến CCD.
  • Cấp độ trung mô (Meso-scale level): Đánh giá động học thoát nước tĩnh (drainage kinetics) và hành vi lưu biến nhớt kế cắt trượt.
  • Cấp độ vĩ mô (Macro-scale coreflood level): Khảo sát động lực học dòng chảy và hệ số thu hồi dầu trên các mô hình mô phỏng vỉa cát nén (sandpack models) đơn lõi và song song (parallel dual-core systems).

Quy trình nghiên cứu rigorous

Quy trình chế tạo và phân tích mẫu được tiêu chuẩn hóa cao độ nhằm loại bỏ sai số hệ thống:

  • Quy trình tạo vi bọt: 200 mL dung dịch cơ sở (gồm XG và SDS hòa tan bằng máy khuấy từ Dlab MS-H280-Pro ở tốc độ 1000 rpm trong 2 giờ) được nạp vào bình chứa dung tích 300 mL. Khí CO2 tinh khiết 99,9% được phân tán qua đĩa xốp khuếch tán đáy (porous stone diffuser) kết hợp cánh khuấy đồng hóa siêu tốc (overhead homogenizer HG-200 Hsiangtai, As One Corporation) ở vận tốc chính xác 8000 rpm trong 4 phút tại nhiệt độ phòng ($25^\circ\text{C}$).
  • Thu nhận và xử lý ảnh số: Mẫu vi bọt được trích xuất ngay tại thời điểm $t = 0\text{ min}$ và $t = 60\text{ min}$, chuyển lên lam kính hiển vi. Ảnh kỹ thuật số độ phân giải cao được chuyển đổi sang định dạng 8-bit và áp dụng thuật toán phân ngưỡng tương phản (contrast thresholding) trên phần mềm ImageJ. Kích thước của tối thiểu 1000 vi bọt trên mỗi mẫu được đo đạc tự động nhằm đảm bảo tính đại diện thống kê tuyệt đối.
  • Kiểm định độ tin cậy và phân phối xác suất: Dữ liệu đường kính vi bọt được xử lý trên nền tảng MATLAB để khớp với ba hàm mật độ xác suất (PDF): Chuẩn (Normal), Log-chuẩn (Log-normal) và Weibull. Độ tin cậy khớp hàm được đánh giá bằng kiểm định Anderson–Darling (AD test) theo Stephens (1974):

$$F_{\text{Weibull}}(x; a, b) = 1 - \exp\left[ -\left(\frac{x}{b}\right)^a \right]$$

với $a$ là thông số hình dạng (shape parameter) và $b$ là thông số tỷ lệ (scale parameter).

+-------------------------------------------------------------------------------+
|                       SƠ ĐỒ QUY TRÌNH THỰC NGHIỆM CHI TIẾT                    |
+-------------------------------------------------------------------------------+
|                                                                               |
|   [Pha chế dung dịch cơ sở]                                                   |
|   - Biopolymer Xanthan Gum (1-5 g/L) + SDS (1-3 g/L) + NaCl (0-20 g/L)        |
|   - Khuấy từ 1000 rpm trong 2 giờ (Máy khuấy Dlab MS-H280-Pro)               |
|                                  |                                            |
|                                  v                                            |
|   [Tạo hệ vi bọt CO2 CGAs]                                                    |
|   - Nạp 200 mL dung dịch vào bình 300 mL                                      |
|   - Sục khí CO2 99.9% qua đĩa phân tán khí xốp                                |
|   - Đồng hóa cắt cao: 8000 rpm trong 4 phút (HG-200 Hsiangtai)                |
|                                  |                                            |
|                                  v                                            |
|   +------------------------------+------------------------------+             |
|   |                                                             |             |
|   v                                                             v             |
|   [Đặc trưng hóa vi mô & độ bền]              [Đánh giá lưu biến học]         |
|   - Đo thoát nước tĩnh (Tính K, T1/2)         - Đo nhớt kế theo tốc độ trượt  |
|   - Kính hiển vi truyền quang + CCD           - Khớp mô hình Power-law        |
|   - ImageJ xử lý ảnh 8-bit (>1000 hạt)        - Khảo sát tỷ lệ Khí : Lỏng     |
|   - Khớp PDF Weibull / Test AD trên MATLAB                                    |
|                                  |                              |             |
|                                  +--------------+---------------+             |
|                                                 |                             |
|                                                 v                             |
|   [Thực nghiệm dòng chảy & Nâng cao thu hồi dầu (Coreflooding)]               |
|   - Mô hình cát nén đơn lõi & Song song hai lõi (Tỷ lệ k = 1:2 và 1:4)        |
|   - Đo áp suất chênh lệch delta-P, kiểm soát dòng phân đoạn (Fractional flow) |
|   - Bơm ngập nước (3 PV) -> Bơm vi bọt CO2 (0.5 PV) -> Bơm rửa đuổi nước      |
|   - Định lượng hệ số thu hồi dầu (% OOIP) tại 45 độ C                         |
|                                                                               |
+-------------------------------------------------------------------------------+

Data và phân tích

Thiết kế mẫu thực nghiệm bao gồm 7 hệ mẫu chuẩn hóa (S1 đến S7) với thành phần biến thiên có hệ thống:

  • S1 (5 g/L XG, 1 g/L SDS), S2 (5 g/L XG, 2 g/L SDS), S3 (5 g/L XG, 3 g/L SDS)
  • S4 (3 g/L XG, 3 g/L SDS), S5 (1 g/L XG, 3 g/L SDS)
  • S6 (5 g/L XG, 3 g/L SDS, 10 g/L NaCl), S7 (5 g/L XG, 3 g/L SDS, 20 g/L NaCl)

Kết quả kiểm định Anderson–Darling (AD) và đồ thị Quantile-Quantile (Q-Q plots) xác nhận:

  • Đối với tất cả 7 mẫu thực nghiệm, hàm phân phối Weibull đều cho giá trị thống kê AD thấp nhất và $P\text{-value} > 0,05$, khẳng định sự tương thích hoàn hảo giữa lý thuyết và thực nghiệm.
  • Ngược lại, hàm phân phối Normal và Log-normal đều cho giá trị AD rất cao và $P\text{-value} < 0,005$, hoàn toàn bác bỏ giả định phân bố chuẩn thông thường.

Hệ thống mô phỏng vỉa cát nén kép (dual-core sandpacks) được thiết lập với các tỷ lệ bất đồng nhất thấm biến thiên (độ thấm cao: $k \approx 4,0\text{ Darcy}$, độ thấm thấp: $k \approx 1,0\text{ Darcy}$ hoặc $2,0\text{ Darcy}$). Quá trình ngập lụt được thực hiện ở lưu lượng kiểm soát nghiêm ngặt thông qua bơm cao áp chuyên dụng và ghi nhận chênh áp liên tục bằng cảm biến áp suất điện tử.

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

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

Thực nghiệm đem lại 5 phát hiện mang tính bước ngoặt với đầy đủ dữ liệu minh chứng:

  1. Quy luật chi phối kích thước hạt của chất hoạt động bề mặt: Tăng nồng độ SDS từ 1 lên 3 g/L (ở mức 5 g/L XG) làm giảm đường kính vi bọt trung bình ($D_{\text{avg}}$) từ 63,37 μm xuống giá trị mịn hơn, đồng thời thu hẹp biên độ phân bố kích thước. Trích dẫn nguyên văn dữ liệu luận án:

    $$\text{"increasing the SDS concentration from 1 to 3 g/L decreased K by 19%."}$$

    Hiện tượng này được lý giải bởi sự suy giảm sức căng bề mặt khí-lỏng, gia tăng mật độ phân tử surfactant tại mặt phân giới, tạo điều kiện chia cắt bọt hiệu quả hơn trong quá trình khuấy cắt cao.

  2. Vai trò cấu trúc sống còn của Biopolymer Xanthan Gum: Tăng nồng độ XG từ 1 lên 5 g/L giúp giảm $D_{\text{avg}}$ từ 55,00 μm xuống 47,37 μm. Quan trọng hơn, XG làm giảm hằng số tốc độ thoát nước $K$ tới 78 lần:

    $$\text{"the K value decreased rapidly from } 5.2 \times 10^{-5}\text{ mL}^{-1}\text{min}^{-1}\text{ (a decrease by a factor of 78) as the XG polymer concentration was increased from 0 to 5 g/L."}$$

    Điều này chứng minh polymer không chỉ đóng vai trò tăng độ nhớt môi trường liên tục mà còn gia cố cơ học cho vỏ bọt Aphron, ngăn ngừa sự hợp giọt (coalescence).

  3. Hiệu ứng phi tuyến của độ mặn khoáng hóa (NaCl): Bổ sung NaCl ở nồng độ thấp (10 g/L) giúp tăng nhẹ độ bền vi bọt do nén lớp điện tích kép, làm giảm lực đẩy tĩnh điện nội tại giữa các nhóm đầu ion âm sulfate. Tuy nhiên, khi tăng độ mặn lên 20 g/L, hiện tượng mất ổn định diễn ra nghiêm trọng: hằng số $K$ tăng hơn 2 lần, đường kính vi bọt trung bình tăng vọt lên 51,80 μm do sự suy giảm độ nhớt của dung dịch polymer và hiệu ứng trọng trường thúc đẩy phân tầng lỏng-khí.

  4. Cơ chế nút kín chọn lọc qua hiệu ứng Jamin: Trong mô hình dòng chảy môi trường xốp, chất lưu vi bọt CO2 thể hiện hành vi Power-law rõ nét. Khả năng cản trở dòng chảy đạt cực đại tại tỷ lệ thể tích Khí : Lỏng bằng 20%. Khi đi vào các kênh dẫn thấm cao, vi bọt bị biến dạng khi đi qua các cổ lỗ xốp hẹp, sinh ra lực cản mao dẫn Jamin ngược chiều chuyển động. Lực cản này tạo ra áp suất chênh lệch đủ lớn để triệt tiêu hiện tượng phân luồng khí sớm, điều chỉnh tỷ lệ dòng phân đoạn trong mô hình vỉa bất đồng nhất hai lõi đạt mức cân bằng lý tưởng 50% : 50%.

  5. Đột phá về hiệu số thu hồi dầu (EOR Recovery Efficiency):

    • Trong mô hình cát nén đơn lõi ở 45°C: Sau khi ngập nước 3 thể tích lỗ rỗng (Pore Volume - PV) chỉ thu hồi được 61,4% OOIP, việc bơm bổ sung 0,5 PV vi bọt CO2 kết hợp ngập nước đuổi tiếp theo đã gia tăng hệ số thu hồi dầu thêm 23,6% OOIP.
    • Trong mô hình vỉa bất đồng nhất song song tỷ lệ thấm 1:4: Bơm ép vi bọt CO2 đạt tổng hệ số thu hồi dầu lên tới 86,9% OOIP, vượt trội hoàn toàn so với mức 75,0% khi chỉ bơm dung dịch polymer cơ sở:

    $$\text{"injection of CO2 microbubbles improved the total oil recovery up to 86.9% compared to the injection of base solution with 75%."}$$

    • Khi tỷ lệ bất đồng nhất giảm xuống 1:2, hệ số thu hồi dầu bằng vi bọt CO2 đạt mức kỷ lục 93,0% OOIP.
+-------------------------------------------------------------------------------+
|             BẢNG SO SÁNH HIỆU QUẢ THU HỒI DẦU TRONG CÁC MÔ HÌNH VỈA           |
+-------------------------------------------------------------------------------+
| Kịch bản thực nghiệm              | Hệ số thu hồi (% OOIP) | Cơ chế chính     |
+-----------------------------------+------------------------+------------------+
| Đơn lõi: Ngập nước cơ sở (3 PV)   |         61.4%          | Đẩy thể tích     |
| Đơn lõi: Ngập nước + Vi bọt CO2   |         85.0% (+23.6%) | Nút kín Jamin    |
| Bất đồng nhất (1:4): Dung dịch base|        75.0%          | Tăng nhớt pha đẩy|
| Bất đồng nhất (1:4): Vi bọt CO2   |         86.9%          | Tái lập dòng quét|
| Bất đồng nhất (1:2): Vi bọt CO2   |         93.0%          | Quét triệt để    |
+-------------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình toán - lý hoàn chỉnh giải thích hành vi động học và phân bố hạt của hệ vi bọt CGA dưới ảnh hưởng đồng thời của chất hoạt động bề mặt, biopolymer và chất điện giải.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình đánh giá vi bọt kết hợp vi ảnh quang học kỹ thuật số, phân tích ảnh tự động ImageJ và kiểm định phân phối thống kê Anderson–Darling.
  • Về mặt thực tiễn vận hành: Thiết lập công thức phối trộn tối ưu ($5\text{ g/L XG} + 3\text{ g/L SDS}$, tỷ lệ khí : lỏng 20%) cho các chiến dịch bơm ép EOR mỏ dầu thực tế, giúp tiết kiệm lượng hóa chất tiêu thụ trong khi vẫn tối đa hóa khả năng kiểm soát độ linh động của khí.
  • Về mặt chính sách và môi trường: Cung cấp lộ trình kỹ thuật khả thi để hiện thực hóa cam kết CCUS kép: vừa tận thu dầu mỏ ở giai đoạn khai thác bậc ba suy kiệt, vừa cố định vĩnh viễn hàng triệu tấn CO2 trong lòng đất dưới dạng bọt bẫy mao dẫn an toàn.

Limitations và Future Research

Luận án thừa nhận khách quan 4 giới hạn nghiên cứu chính:

  1. Giới hạn áp suất và nhiệt độ buồng thí nghiệm: Các thí nghiệm dòng chảy được tiến hành tại 45°C và dải áp suất trung bình. Trong các vỉa dầu sâu thực tế, nhiệt độ có thể vượt quá 80–120°C và áp suất có thể đạt mức siêu tới hạn ($P > 7,38\text{ MPa}$), nơi CO2 chuyển sang trạng thái lỏng/siêu tới hạn (supercritical CO2) với tính chất hòa tan hoàn toàn khác biệt.
  2. Đặc tính môi trường xốp: Mô hình cát nén nhân tạo (sandpacks) có cấu trúc lỗ rỗng đồng đều và tính thấm ướt nước (water-wet) lý tưởng, chưa phản ánh đầy đủ độ phức tạp của đá vỉa carbonate hoặc sa thạch tự nhiên chứa vi nứt nẻ và có tính thấm ướt hỗn hợp (mixed-wet/oil-wet).
  3. Thành phần hóa học lưu chất: Sử dụng muối đơn thành phần NaCl để mô phỏng nước vỉa, chưa xét đến tác động của các ion đa hóa trị ($\text{Ca}^{2+}$, $\text{Mg}^{2+}$, $\text{SO}_4^{2-}$) vốn có khả năng tạo kết tủa hoặc phá vỡ cấu trúc màng polymer sinh học.
  4. Quy mô không gian: Nghiên cứu mới dừng lại ở quy mô mẫu lõi phòng thí nghiệm (lab-scale corefloods), chưa tính đến hiệu ứng suy giảm độ nhớt do ứng suất cắt kéo dài và phân tán khoảng cách xa trong giếng mỏ thực địa.

Chương trình nghiên cứu mở rộng trong 5-10 năm tới cần tập trung vào:

  • Phát triển hệ chất hoạt động bề mặt và polymer chịu nhiệt - chịu mặn cao cấp (như biến tính nano silica hoặc polymer tổng hợp bền nhiệt).
  • Thực hiện thí nghiệm bơm ép vi bọt CO2 ở trạng thái siêu tới hạn trên hệ thiết bị lõi đá thật (High-Pressure High-Temperature - HPHT Coreflooding).
  • Xây dựng mô hình mô phỏng số 3D quy mô vỉa (Reservoir Numerical Simulation) tích hợp động học vi bọt CGAs.
  • Đánh giá khả năng bẫy lưu trữ địa chất lâu dài của CO2 vi bọt sau giai đoạn đóng mỏ.

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

  • Tác động học thuật: Luận án mở ra hướng tiếp cận mới trong cơ học lưu chất mỏ và kỹ thuật phân tán keo. Các bài báo công bố từ luận án trên các tạp chí quốc tế chuyên ngành (Q1/Q2) dự kiến thu hút lượng trích dẫn cao từ cộng đồng nghiên cứu EOR, CCUS và kỹ thuật bề mặt mỏ dầu.
  • Chuyển đổi công nghiệp: Ngành công nghiệp dầu khí tại Đông Nam Á, Nhật Bản và toàn cầu có thể trực tiếp ứng dụng quy trình công nghệ tạo vi bọt để cải tạo các giếng khai thác ngập nước suy kiệt, giảm thiểu chi phí năng lượng bơm nén khí CO2 tự do và hạn chế tối đa rủi ro ăn mòn đường ống do bọt khí phân tầng.
  • Đóng góp môi trường và xã hội: Tối ưu hóa chu trình kinh tế tuần hoàn carbon (Circular Carbon Economy). Phương pháp vi bọt biến CO2 từ một khí thải nhà kính nguy hại thành tác nhân kinh tế hữu ích, hỗ trợ đắc lực cho mục tiêu Net-Zero Carbon 2050 thông qua lưu trữ địa chất kết hợp khai thác dầu bền vững.

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

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Dầu khí/Khoa học Trái đất: Tiếp cận cơ sở dữ liệu thực nghiệm chuẩn xác, phương pháp luận kiểm định phân bố Weibull và khung phân tích tích hợp Jamin - Young-Laplace.
  • Kỹ sư trưởng và Chuyên gia R&D Mỏ Dầu khí: Sở hữu công thức pha chế lưu chất vi bọt tối ưu, các biểu đồ biến thiên áp suất và lưu lượng để thiết kế phương án ngập lụt vi bọt tại hiện trường.
  • Các nhà hoạch định chính sách Năng lượng & Môi trường: Có thêm luận cứ khoa học thực chứng để xây dựng khung pháp lý và chính sách ưu đãi thuế tín chỉ carbon cho các dự án CCUS kết hợp CO2-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 lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Thuyết cấu trúc hạt khí keo Aphron (CGAs) của Sebba (1987) và động học thoát nước của Yan et al. (2005) sang hệ vi bọt CO2-Biopolymer-Surfactant. Luận án đã mô hình hóa toán học thành công cơ chế ức chế hiện tượng Ostwald ripening nhờ sự gia cố màng ba lớp bởi polymer Xanthan Gum, đồng thời chứng minh bằng thống kê rằng phân bố kích thước vi bọt CO2 tuân theo quy luật hàm Weibull thay vì hàm Gaussian hay Log-normal truyền thống.

2. Điểm đổi mới trong phương pháp luận so với ít nhất hai nghiên cứu quốc tế trước đây?

So với nghiên cứu của Zhu et al. (2020) (chỉ áp dụng phân phối Gauss thông thường cho bọt khoan mà không qua kiểm định độ tin cậy) và nghiên cứu của Alam et al. (2017) (sử dụng Log-normal cho bọt tuyển nổi), luận án của Lê Nguyễn Hải Nam đã áp dụng quy trình chụp ảnh truyền quang kỹ thuật số kết hợp xử lý tương phản tự động trên ImageJ với cỡ mẫu lớn ($N \ge 1000$ hạt), sau đó thực hiện kiểm định Anderson–Darling (AD test) trên MATLAB. Kết quả kiểm định với giá trị $P\text{-value} > 0,05$ đã chứng minh tính chuẩn xác vượt trội của hàm phân phối Weibull trong việc mô tả kích thước vi bọt CO2.

3. Phát hiện bất ngờ nhất (counter-intuitive result) thu được từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất nằm ở phản ứng phi tuyến của hệ vi bọt đối với nồng độ muối NaCl. Thông thường, chất điện phân được cho là làm giảm độ bền bọt do phá vỡ lớp điện tích. Tuy nhiên, dữ liệu thực nghiệm chứng minh nồng độ NaCl 10 g/L lại giúp tăng nhẹ độ ổn định và giảm kích thước vi bọt nhờ hiệu ứng nén lớp điện tích kép làm giảm lực đẩy tĩnh điện giữa các đầu ion sulfate. Chỉ khi vượt qua ngưỡng tới hạn lên 20 g/L NaCl, hệ mới bị phân rã nhanh chóng do suy giảm độ nhớt và hiệu ứng trọng trường gây kết tụ bọt.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (replication protocol) chi tiết không?

Có, luận án cung cấp quy trình thực nghiệm chi tiết tuyệt đối: từ thiết bị đồng hóa siêu tốc (HG-200 Hsiangtai ở 8000 rpm trong đúng 4 phút), đĩa khuếch tán khí xốp, thông số khuấy từ (1000 rpm trong 2 giờ), nồng độ hóa chất chính xác từng gam/lít, quy trình phân ngưỡng 8-bit trên ImageJ, đến thiết kế mô hình cát nén đo áp suất và lưu lượng ở 45°C. Mọi phòng thí nghiệm cơ bản về cơ học lưu chất mỏ đều có thể tái lập hoàn toàn quy trình này.

5. Chương trình nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghiên cứu 10 năm tập trung vào 3 trọng tâm:

  1. Nâng cấp hệ hóa chất vi bọt chịu nhiệt độ và độ mặn cao (High-Temperature High-Salinity) bằng cách ghép nối nano-silica hoặc polymer biến tính.
  2. Thử nghiệm bơm ép CO2 vi bọt ở điều kiện áp suất - nhiệt độ siêu tới hạn (Supercritical CO2 HPHT) trên mẫu lõi đá vỉa thực tế (sandstone/carbonate cores).
  3. Phát triển module thuật toán mô phỏng dòng vi bọt phi Newton trong môi trường bất đồng nhất 3D để tích hợp vào các phần mềm thương mại mô phỏng vỉa dầu khí (như Eclipse hoặc CMG).

Kết luận

Luận án tiến sĩ của Lê Nguyễn Hải Nam đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn cốt lõi:

  1. Thiết lập công thức dung dịch vi bọt CO2 tối ưu: Xác định tỷ lệ phối hợp vàng giữa 5 g/L Xanthan Gum và 3 g/L SDS, giúp giảm tốc độ thoát nước tĩnh tới 78 lần và duy trì độ bền cấu trúc vượt trội.
  2. Khám phá và chứng minh phân bố Weibull: Lần đầu tiên chứng minh bằng thống kê thực nghiệm rằng phân bố kích thước vi bọt CO2 tuân theo hàm phân phối xác suất Weibull với độ tin cậy tuyệt đối qua kiểm định Anderson–Darling.
  3. Làm sáng tỏ cơ chế nút kín chọn lọc Jamin: Chứng minh vi bọt CO2 hành xử như một lưu chất Power-law, kích hoạt hiệu ứng cản mao dẫn Jamin tại đới thấm cao và cân bằng dòng phân đoạn về mức 50% : 50% trong mô hình vỉa bất đồng nhất.
  4. Đạt bước đột phá về hệ số thu hồi dầu: Nâng cao hệ số thu hồi dầu thêm 23,6% OOIP sau giai đoạn ngập nước trong mô hình đơn lõi, và đạt tổng hệ số thu hồi dầu từ 86,9% đến 93,0% OOIP trong các mô hình vỉa bất đồng nhất tỷ lệ thấm 1:4 và 1:2.
  5. Mở ra giải pháp kép cho kinh tế carbon: Kết hợp hoàn hảo giữa mục tiêu gia tăng sản lượng khai thác năng lượng thứ cấp/tam cấp và chiến lược lưu trữ CO2 địa chất an toàn, bền vững.

Công trình tạo tiền đề vững chắc mở ra ba nhánh nghiên cứu mới: vi bọt CO2 siêu tới hạn, dung dịch vi bọt nano hóa, và mô phỏng số dòng vi bọt quy mô mỏ, đóng góp thiết thực cho nền khoa học kỹ thuật mỏ dầu khí toàn cầu.