Tổng quan về luận án

Nghiên cứu đặc tính thạch vật lý và động lực học dòng chảy trong các tầng chứa đá cacbonat đóng vai trò then chốt trong ngành công nghiệp dầu khí toàn cầu. Theo các công trình kinh điển của Akbar và cộng sự (2000) cũng như Aljuboori và cộng sự (2019), các bồn trũng chứa cacbonat lưu trữ khoảng 60% đến 70% tổng trữ lượng dầu khí thế giới, tiêu biểu tại các bồn trũng Trung Đông, Caspien và Bắc Mỹ. Tại thềm lục địa Việt Nam, mặc dù các cấu tạo cacbonat có tính phân bố chọn lọc hơn cát kết, các phát hiện quan trọng như mỏ Cá Voi Xanh, Hàm Rồng (bể Sông Hồng), Đại Hùng, Sao Vàng - Đại Nguyệt (bể Nam Côn Sơn) và mỏ CX đã khẳng định tiềm năng kinh tế - kỹ thuật mang tính chiến lược của đối tượng này.

Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Địa vật lý (Mã số: 9520502) của tác giả Nguyễn Trung Dũng, thực hiện tại Trường Đại học Mỏ - Địa chất dưới sự hướng dẫn khoa học của PGS. Phan Thiên Hương và TS. Cù Minh Hoàng, tập trung vào đề tài: "Đặc tính tầng chứa đá cacbonat tuổi Mioxen giữa mỏ CX". Đối tượng nghiên cứu là thể địa chất cacbonat thuộc hệ tầng Tri Tôn phân bố dọc theo đới nâng Tri Tôn – một địa lũy phát triển theo phương Tây Bắc - Đông Nam ở phần phía Nam bể trầm tích Sông Hồng, nằm cách bờ biển Quảng Nam - Quảng Ngãi khoảng 80 km.

Khoảng trống nghiên cứu (Research Gap) lớn nhất xuất phát từ tính dị hướng và bất đồng nhất trầm tích - thạch học cực kỳ phức tạp của đá vôi khối xây rạn san hô (reef build-up) hệ tầng Tri Tôn. Các phương pháp minh giải địa vật lý giếng khoan (ĐVLGK) và thạch vật lý truyền thống áp dụng mô hình điện trở kinh điển Archie (1942) gặp phải sự suy giảm độ chính xác nghiêm trọng. Lý do là hệ số gắn kết xi măng $m$ và hằng số bão hòa $n$ tại mỏ CX biến động rất mạnh trong khoảng rộng ($m$ từ dưới 2,0 đến trên 2,7; $n$ từ 1,9 đến trên 3,0), trong khi quan hệ rỗng - thấm rời rạc không thể biểu diễn bằng một hàm hồi quy đơn biến. Trước luận án này, các công trình tại khu vực đới nâng Tri Tôn (Vũ Ngọc Diệp và nnk, 2011; Nguyễn Xuân Phong và nnk, 2017) chủ yếu tập trung vào tiến hóa sinh tướng, môi trường trầm tích hoặc địa tầng phân tập, hoàn toàn thiếu vắng một chu trình định lượng hóa các tham số thạch vật lý kết hợp kỹ thuật trí tuệ tính toán.

Luận án giải quyết ba câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Làm thế nào để phân tách tối ưu các nhóm đá chứa có cùng đặc tính thủy lực trong đá cacbonat có cấu trúc lỗ rỗng dị thường tại mỏ CX?
  2. RQ2: Mô hình học máy có giám sát nào đạt hiệu năng cao nhất trong việc mô phỏng và liên tục hóa giá trị chỉ số dòng chảy FZI và độ thấm tuyệt đối $K$ dọc thân giếng dựa trên tập hợp đường cong ĐVLGK hạn chế?
  3. RQ3: Làm thế nào để xây dựng mô hình bão hòa nước theo chiều cao (Saturation Height Function - SHF) độc lập với các tham số Archie biến động nhằm nâng cao độ tin cậy khi tính toán thể tích khí tại chỗ (GIIP)?

Các giả thuyết nghiên cứu được kiểm chứng gồm:

  • H1: Tầng chứa cacbonat hệ tầng Tri Tôn mỏ CX cấu thành từ 5 đơn vị dòng chảy thủy lực (HFU) khác biệt có thể bóc tách bằng các thuật toán phân cụm không giám sát kết hợp chỉ số FZI và RQI.
  • H2: Thuật toán hồi quy phi tuyến tính quá trình Gaussian (Gaussian Process Regression) kết hợp phân nhóm HFU mang lại độ chính xác dự báo độ thấm ($R^2 \ge 0,70$) vượt trội so với các mô hình tuyến tính cổ điển.
  • H3: Việc áp dụng hàm bão hòa Skelt-Harrison riêng biệt cho từng HFU giúp loại bỏ sai số tích lũy do độ rỗng thứ sinh và cấu trúc hang hốc, kiểm chứng sự suy giảm độ bão hòa nước $S_w$ xuống dưới 6% tại các vỉa khí dày.

Quy mô thực nghiệm của luận án bao quát toàn bộ tài liệu 4 giếng khoan thăm dò (GK1, GK2, GK3, GK4) với tập hợp đo ĐVLGK đầy đủ (GR, RHOB, NPHI, DTC, DTS, RD, RS, MSFL), hàng trăm mẫu lõi phân tích thông thường (RCAL) và phân tích mẫu đặc biệt (SCAL bao gồm đo áp suất mao dẫn MICP, đo điện trở suất vỉa xác định $m, n$), kết hợp tài liệu thử vỉa hiện đại MDT/RCI và DST. Toàn bộ công trình được cấu trúc thành 127 trang, 71 hình vẽ, 12 bảng biểu số liệu và hệ thống mã nguồn Python xử lý chuyên biệt.

Literature Review và Positioning

Khung lý thuyết phân loại đá chứa cacbonat trong lịch sử địa chất dầu khí phát triển qua ba trường phái tiếp cận chủ đạo:

  • Trường phái Trầm tích - Thạch học: Dẫn đầu bởi Folk (1959), Dunham (1962), và cải tiến của Embry & Klovan (1971), phân loại dựa trên năng lượng môi trường lắng đọng, tỷ lệ bùn/hạt (mudstone, wackestone, packstone, grainstone, boundstone).
  • Trường phái Cấu trúc Lỗ rỗng: Archie (1952), Choquette & Pray (1970) với khái niệm lỗ rỗng chọn lọc cấu trúc (fabric selective) và không chọn lọc cấu trúc (nonfabric selective); Lucia (1983, 1995) phân chia lỗ rỗng giữa hạt (interparticle) và hang hốc (vuggy: separate vugs và touching vugs) gắn với chỉ số kết cấu đá (Rock Fabric Number - RFN 1, 2, 3).
  • Trường phái Thủy lực Đơn vị dòng chảy (Hydraulic Flow Units): Khởi nguồn từ Ebanks (1984), được toán học hóa bởi Amaefule và cộng sự (1993) qua phương trình biến đổi Kozeny-Carman để xác định Chỉ số chất lượng đá (Rock Quality Index - RQI) và Chỉ báo đới dòng chảy (Flow Zone Indicator - FZI).

Trong y văn tồn tại mâu thuẫn sâu sắc giữa phương pháp phân loại thạch học thuần túy và phân loại thủy động lực học. Rebelle & Lalanne (2014) chỉ rõ rằng ranh giới các tướng trầm tích học không nhất thiết trùng khớp với ranh giới phân bố động lực học dòng chảy trong đá cacbonat bị biến đổi thứ sinh mạnh mẽ (dolomit hóa, hòa tan nứt nẻ). Phương pháp Lucia (1983, 1995) bộc lộ nhược điểm nghiêm trọng khi áp dụng cho vỉa chứa có vi lỗ rỗng (microporosity) và hang hốc cô lập (separate vugs) vì không thể đánh giá định lượng chính xác mạng lưới kết nối qua các lát mỏng 2D. Tương tự, phương pháp Winland R35 và Pittman (1992) dựa trên bán kính họng lỗ rỗng đo bằng MICP tại điểm bão hòa thủy ngân 35% bị Spearing và cộng sự (2001) phê phán là mang tính kinh nghiệm giả định, thiếu cơ chế động học tổng quát.

Tại Việt Nam, các nghiên cứu của Lê Hải An và cộng sự (2016) tại bể Phú Khánh, Hà Quang Mẫn và cộng sự (2021, 2023) tại bể Nam Côn Sơn, Lê Trung Tâm và cộng sự (2015) tại bể Sông Hồng đã bắt đầu đưa khái niệm HFU vào đánh giá vỉa. Tuy nhiên, phần lớn các công trình chỉ dừng lại ở các thuật toán cơ bản như hồi quy đa biến hoặc mạng nơ-ron truyền thống (ANN) (Nguyễn Thu Huyền và cộng sự, 2019), chưa xây dựng được chu trình khép kín từ phân cụm không giám sát đa thuật toán đến dự báo độ thấm phi tuyến và mô hình hóa bão hòa nước theo chiều cao.

Luận án định vị tiên phong bằng việc tích hợp cơ học dòng chảy vi mô với trí tuệ nhân tạo hiện đại. Khi so sánh với các công trình quốc tế quy mô lớn tại bồn trũng Caspien (mỏ Tengiz, Kashagan - vỉa cacbonat cổ biến dạng nứt nẻ sâu) và bồn trũng Ả Rập (mỏ Ghawar - vỉa cacbonat hạt thống trị), nghiên cứu của luận án giải quyết thành công bài toán dữ liệu giếng khoan hạn chế (sparse well data) bằng việc tối ưu hóa tổ hợp thuật toán học máy phân cấp và hàm SHF phi Archie, tạo ra một tiêu chuẩn xử lý mới cho các bồn trũng ven rìa Tây Thái Bình Dươ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 phong phú lý thuyết biến đổi Kozeny-Carman và khung phân tích Amaefule (1993) trong môi trường cacbonat sinh quyển nhiệt đới phức tạp. Dựa trên mối quan hệ giữa độ rỗng hữu dụng $\Phi_e$ và độ thấm tuyệt đối $K$, phương trình Kozeny-Carman được chuyển đổi thành mối quan hệ tuyến tính hóa giữa Chỉ số chất lượng đá RQI và Tỷ số thể tích lỗ rỗng $\Phi_z$:

$$\Phi_z = \frac{\Phi_e}{1 - \Phi_e}$$

$$RQI = 0.0314 \sqrt{\frac{K}{\Phi_e}}$$

$$FZI = \frac{RQI}{\Phi_z} = \frac{1}{\sqrt{F_s \tau^2 S_{gv}}}$$

Trong đó, $F_s$ là hệ số hình dạng lỗ rỗng, $\tau$ là độ cong mao dẫn, $S_{gv}$ là diện tích bề mặt hạt trên một đơn vị thể tích hạt. Luận án đã chứng minh trên dữ liệu thực nghiệm rằng tham số $FZI$ là một đại lượng đặc trưng thủy lực bất biến phản ánh bản chất vi cấu trúc hình học khoảng rỗng, độc lập với độ sâu chôn vùi.

log(RQI)

Luận án thách thức trực tiếp mô hình điện trở Archie:

$$S_w = \sqrt[n]{\frac{a \cdot R_w}{\Phi^m \cdot R_t}}$$

Bằng chứng thực nghiệm mỏ CX cho thấy hệ số xi măng hóa $m$ phân tán cực mạnh (trọng tâm 2,2 - 2,5, cực trị 2,7) và $n$ dao động (1,9 - 3,0), phản ánh sự đa dạng của hệ thống hang hốc tách rời (separate vugs) và hang hốc kết nối (touching vugs). Luận án xác lập mệnh đề lý thuyết: Trong đá cacbonat có hệ số Archie bất định, việc mô hình hóa phân bố chất lưu phải được chuyển đổi từ mô hình điện trở thuần túy sang mô hình thủy tĩnh mao dẫn động lực học SHF phân tầng theo từng đơn vị HFU.

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

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

  1. Lý thuyết Thủy lực Vi mô: Phân tích FZI-RQI của Amaefule (1993) nhằm cấu trúc hóa không gian lỗ rỗng thành các tập hợp đẳng thủy lực.
  2. Lý thuyết Cân bằng Áp suất Mao dẫn: Sử dụng mô hình độ cao bão hòa Skelt-Harrison (1995) để mô tả động học phân bố chất lưu:

$$\log(S_w) = - a \cdot \exp\left(-b \cdot h\right) + c \cdot h^{-d}$$

Trong đó $h$ là chiều cao phía trên mực nước tự do (Height Above Free Water Level - HAFWL). Mô hình này vượt trội hoàn toàn so với các hàm Brooks-Corey, Lambda, Leverett J-Function và Thomeer khi khớp với dữ liệu áp suất mao dẫn mỏ CX. 3. Lý thuyết Trí tuệ Tính toán Đa tầng: Tích hợp chuỗi phân cụm không giám sát (Unsupervised Clustering) làm cơ sở gán nhãn, tiếp nối bằng mô hình hồi quy phi tuyến có giám sát (Supervised Regression) xử lý mối quan hệ phi tuyến cao giữa tập tín hiệu ĐVLGK và đặc tính rỗng - thấm.

Điều kiện biên lý thuyết (Boundary Conditions): Khung phân tích này đạt hiệu lực tối ưu trên các đá cacbonat tuổi Kainozoi sinh thành trong môi trường biển nông nhiệt đới, có cấu trúc rạn san hô - tảo, chịu biến đổi thứ sinh hòa tan và nứt nẻ kiến tạo cục bộ.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu tuân thủ chặt chẽ lập trường nhận thức luận Thực chứng (Positivism) kết hợp phương pháp luận quy nạp dựa trên dữ liệu (Data-driven inductive approach) và cơ chế tất định thạch vật lý (Petrophysical determinism). Quy trình nghiên cứu đa mức (multi-level design) được thiết kế từ vĩ mô (mặt cắt địa chấn, phân vùng cấu trúc địa lũy Tri Tôn) đến trung mô (tài liệu đo carota giếng khoan ĐVLGK) và vi mô (phân tích lát mỏng thạch học quang học, chụp ảnh hiển vi điện tử quét SEM, đo rỗng thấm áp suất mao dẫn trên mẫu lõi 1 inch và mẫu lõi nguyên đường kính whole-core).

Quy trình nghiên cứu rigorous

Quy trình kỹ thuật được tiến hành nghiêm ngặt qua các bước tiêu chuẩn công nghiệp:

  1. Kiểm chuẩn và Chuẩn hóa Dữ liệu: Kiểm tra chất lượng (QC) tài liệu giếng khoan tại 4 giếng (GK1, GK2, GK3, GK4), loại bỏ các khoảng sạt lở thành giếng dựa trên đường kính Caliper; thực hiện hiệu chỉnh môi trường (environmental correction) và chuẩn hóa phân bố thống kê (normalization) các đường đo mật độ (RHOB), nơ-tron (NPHI), âm lượng sóng nén/sóng ngang (DTC/DTS), điện trở suất sâu/nông/vi phân (RD/RS/MSFL).
  2. Đồng bộ hóa Không - Thời gian giữa Lõi và Log: Chuyển dịch độ sâu mẫu lõi về độ sâu tuyệt đối theo giếng khoan (Core-to-log depth matching) thông qua việc tương quan đường đo Gamma Ray mẫu lõi (Core Gamma) và đường Gamma Ray đo giếng khoan (GR), loại bỏ hoàn toàn các điểm dị biệt (outliers).
  3. Tam giác đạc Dữ liệu (Data Triangulation): Kết hợp chặt chẽ giữa 3 nguồn dữ liệu độc lập: kết quả thử vỉa động học (MDT/RFT/DST), số liệu phân tích mẫu phòng thí nghiệm (RCAL/SCAL) và dữ liệu địa vật lý trường liên tục. Độ tin cậy của phân tích thạch học được xác nhận qua việc kiểm chứng chéo lát mỏng hiển vi cho từng đới HFU.

Data và phân tích

Tập dữ liệu phân tích mẫu lõi bao quát dải giá trị thạch vật lý rộng: độ rỗng mẫu lõi ($\Phi_{core}$) biến thiên từ 2,0% đến trên 30,0%; độ thấm tuyệt đối ($K_{core}$) trải dài qua 4 bậc độ lớn, từ dưới 1,0 mD đến hơn 2.000 mD. Độ sâu nghiên cứu nằm trong khoảng từ 1.000 m đến trên 4.000 m TVDss.

Để xử lý bài toán phân loại và dự báo, luận án đã triển khai hệ thống công cụ trên nền tảng Python (Scikit-learn, SciPy, Pandas, NumPy, GPflow):

  • Phân nhóm không giám sát: Áp dụng 4 thuật toán gồm K-means clustering, Ward’s hierarchical clustering, Self-Organizing Map (SOM) và Fuzzy C-means (FCM). Phương pháp điểm khuỷu tay (Elbow method) và tiêu chuẩn tổng bình phương khoảng cách nội cụm (Within-Cluster Sum of Squares - WCSS) được dùng để xác lập số cụm tối ưu.
  • Hồi quy có giám sát: So sánh 6 thuật toán tiên tiến: Hồi quy quá trình Gaussian (Gaussian Process Regression - GPR với hàm nhân Exponential), Cây quyết định (Decision Trees), Rừng ngẫu nhiên (Random Forest), Mạng nơ-ron nhân tạo (ANN), Máy vector hỗ trợ (SVM) và Mô hình tổ hợp tối ưu (Ensemble Optimizable).
  • Kiểm định Robustness: Áp dụng phương pháp kiểm định chéo $k$-fold cross-validation ($k=5$) và phân chia tập dữ liệu huấn luyện/kiểm tra độc lập (Train/Test split 80/20) để ngăn ngừa hiện tượng quá khớp (overfitting). Đánh giá qua các chỉ số thống kê chặt chẽ: Hệ số xác định ($R^2$), Sai số tuyệt đối trung bình (MAE) và Căn bậc hai sai số bình phương trung bình (RMSE).

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

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

Luận án mang lại 4 phát hiện khoa học mang tính đột phá trên cơ sở bằng chứng thực nghiệm rõ ràng:

  1. Xác lập Cấu trúc 5 Đơn vị Dòng chảy Thủy lực Tối ưu: Bằng kỹ thuật K-means và Ward's clustering trên tập chỉ số FZI, luận án đã phân chia thành công tầng chứa Mioxen giữa mỏ CX thành 5 nhóm HFU với ranh giới định lượng sắc nét:

    • HFU 1: $FZI < 0,70$ $\mu$m (Đá chứa chất lượng rất xấu, vi lỗ rỗng, bùn vôi chiếm ưu thế, Mudstone/Wackestone nén chặt).
    • HFU 2: $0,70 \le FZI < 1,85$ $\mu$m (Chất lượng xấu - trung bình, Wackestone/Packstone).
    • HFU 3: $1,85 \le FZI < 4,00$ $\mu$m (Chất lượng trung bình - khá, Packstone/Grainstone).
    • HFU 4: $4,00 \le FZI < 8,50$ $\mu$m (Chất lượng tốt, Grainstone giàu sinh vật tiêu biểu).
    • HFU 5: $FZI \ge 8,50$ $\mu$m (Chất lượng rất tốt - tuyệt hảo, đá vôi rạn san hô cấu trúc khung boundstone, hang hốc lớn và nứt nẻ mở rộng).

  2. Hiệu năng Vượt trội của Mô hình Gaussian Process Regression (Exponential): Trong việc liên tục hóa dự báo $FZI$ và độ thấm $K$ từ các đường cong ĐVLGK, mô hình GPR với nhân Exponential đạt hệ số tương quan thực nghiệm $R^2 = 0,70$ trên tập kiểm tra độc lập, vượt xa mô hình hồi quy truyền thống ($R^2 < 0,45$). Thuật toán Ensemble Optimizable cũng thể hiện khả năng phân loại chính xác nhãn HFU với độ tin cậy đạt trên 82%.

  3. Giải mã Hiện tượng Bất thường Thạch vật lý do Nứt nẻ Kiến tạo: Tại các giếng GK2 và GK4, luận án phát hiện các khoảng dị thường nghịch biến: độ rỗng ma trận đo được rất thấp ($\Phi < 5%$) nhưng độ thấm thủy lực thực tế vọt lên trên $2.000$ mD. Phân tích lát mỏng và ảnh quét thành giếng xác nhận sự hiện diện của mạng lưới vi nứt nẻ kiến tạo (tectonic micro-fractures) kết nối các hang hốc thứ sinh, giải thích hiện tượng cho dòng khí lưu lượng cao tại các tầng ngỡ như là đá đặc xít.

  4. Sự Vượt trội của Hàm Bão hòa Skelt-Harrison trong Dự báo Phân bố Nước: Khi kiểm chứng khớp đường cong áp suất mao dẫn qua 5 mô hình (Brooks-Corey, J-function, Lambda, Thomeer, Skelt-Harrison), mô hình Skelt-Harrison đạt sai số bình phương cực tiểu trên toàn bộ 5 HFU. Áp dụng mô hình SHF này chứng minh độ bão hòa nước $S_w$ trong đới chuyển tiếp giảm nhanh và duy trì ở mức siêu thấp ($S_w < 6%$) trên chiều cao cột khí tại các khối xây cacbonat dày ở mỏ CX.

Implications đa chiều

  • Về mặt Lý thuyết: Khẳng định tính tất yếu của việc tích hợp lý thuyết dòng chảy vi mô HFU vào các tầng chứa dị hướng cao, cung cấp cơ sở để tái đánh giá các vỉa cacbonat cổ trên thềm lục địa thềm Đông Á.
  • Về mặt Phương pháp luận: Xác lập chu trình phân tích kết hợp Machine Learning có thể chuyển giao (reproducible workflow), cho phép tính toán độ thấm liên tục dọc thân giếng trong điều kiện thiếu vắng lõi trầm tích.
  • Về mặt Thực tiễn và Khai thác: Giúp các nhà vận hành tối ưu hóa khoảng bắn mở vỉa (perforation intervals), tập trung chính xác vào các tầng HFU 4 và HFU 5; đồng thời tránh bắn nhầm vào các đới chuyển tiếp nước - khí có hiện tượng thấm mao dẫn bất thường.
  • Về mặt Quản lý Tài nguyên: Nâng cao độ chính xác của việc ước tính thể tích khí tại chỗ (GIIP), giảm thiểu độ bất định (uncertainty range) trong báo cáo đánh giá trữ lượng phát triển mỏ trình Hội đồng Đánh giá Trữ lượng Quốc gia.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật nội tại:

  1. Số lượng giếng khoan hạn chế: Toàn bộ công trình xây dựng trên dữ liệu của 4 giếng khoan (GK1 đến GK4) tại mỏ CX. Dù độ bao phủ địa tầng tốt, mật độ giếng thưa tạo ra độ bất định nhất định khi ngoại suy sang các cánh xa của đới nâng Tri Tôn.
  2. Đại diện của mẫu lõi SCAL: Số lượng mẫu đo áp suất mao dẫn đặc biệt và số liệu đo điện trở vỉa xác định $m, n$ còn khiêm tốn do chi phí lấy mẫu cao và điều kiện cơ học vỉa vụn vỡ tại một số đới nứt nẻ.
  3. Hiệu ứng thu phóng quy mô (Scale Effect): Dữ liệu ĐVLGK có độ phân giải hàng chục cm trong khi mẫu lõi plug có kích thước 1 inch. Sự chênh lệch quy mô này đôi khi làm mờ các hiệu ứng dị thường nứt nẻ siêu nhỏ.
  4. Mô hình hóa không gian 3D: Luận án mới tập trung giải quyết bài toán mô phỏng 1D dọc thân giếng khoan, chưa tích hợp trực tiếp trường thuộc tính vận tốc địa chấn đảo ngược (seismic acoustic impedance inversion) vào mô hình rỗng - thấm 3D toàn mỏ.

Chương trình nghiên cứu 5-10 năm tiếp theo được định hình qua 4 hướng đột phá:

  • Tích hợp Deep Learning và Học sâu Xử lý Ảnh: Ứng dụng mạng nơ-ron tích chập (Convolutional Neural Networks - CNN) để phân loại tự động ảnh quét thành giếng độ phân giải cao (FMI/CBIL) và ảnh hiển vi điện tử lát mỏng độ phân giải nano.
  • Mô phỏng Mạng Lưới Lỗ rỗng Số (Digital Rock Physics - DRP): Ứng dụng kỹ thuật chụp cắt lớp vi tính tia X vi mô (Micro-CT scan) 3D để mô phỏng trực tiếp dòng chảy hai pha (khí - nước) qua các hang hốc thực tế.
  • Đồng hóa Địa chấn 4D và Vận tốc Sóng: Mở rộng chu trình học máy bằng cách kết hợp trường thuộc tính địa chấn đàn hồi nhiều thành phần để ngoại suy phân bố HFU trong không gian liên giếng 3D.
  • Đánh giá Động thái Thu hồi Nâng cao (EGR/CCUS): Nghiên cứu tương tác thạch hóa học khi bơm ép $CO_2$ vào tầng chứa cacbonat Mioxen phục vụ lưu trữ carbon và gia tăng thu hồi khí.

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

  • Ảnh hưởng Học thuật: Luận án cung cấp hệ phương pháp luận mẫu mực cho ngành Kỹ thuật Địa vật lý và Vật lý Thạch học Dầu khí tại Việt Nam; tạo tiền đề cho hàng loạt nghiên cứu mở rộng tại Trường Đại học Mỏ - Địa chất, Viện Dầu khí Việt Nam (VPI) và các trường đại học quốc tế trong khu vực Đông Nam Á.
  • Chuyển đổi Công nghiệp Dầu khí: Kết quả nghiên cứu được áp dụng trực tiếp cho các chương trình phát triển mỏ của Tập đoàn Dầu khí Việt Nam (Petrovietnam - PVN), Tổng Công ty Thăm dò Khai thác Dầu khí (PVEP) và các nhà điều hành quốc tế như Vietgazprom, Eastsea Star, Schlumberger, giúp tiết kiệm hàng triệu USD chi phí đo kiểm thử vỉa không cần thiết.
  • Chính sách Năng lượng Quốc gia: Nâng cao độ tin cậy trong việc thẩm lượng quy mô trữ lượng khí tự nhiên tại thềm lục địa miền Trung, phục vụ trực tiếp cho quy hoạch chuỗi dự án khí - điện miền Trung và chiến lược an ninh năng lượng quốc gia đến năm 2030, tầm nhìn 2050.
  • Ý nghĩa Quốc tế: Đóng góp dữ liệu thực nghiệm quý giá về hệ thống bồn trũng nhiệt đới Kainozoi khu vực Biển Đông vào kho tàng cơ sở dữ liệu trầm tích cacbonat toàn cầu.

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

  • Nghiên cứu sinh và Giới Học thuật Chuyên sâu: Tiếp cận một chu trình khoa học hoàn chỉnh kết hợp giữa vật lý địa chất kinh điển và khoa học dữ liệu hiện đại, mở ra các hướng đề tài mới về xử lý tín hiệu ĐVLGK.
  • Chuyên gia Thạch vật lý (Petrophysicists) và Kỹ sư Mô hình Mỏ (Reservoir Engineers): Sở hữu công cụ và mã nguồn tính toán thực tế để tự động hóa quy trình phân tách HFU, dự báo độ thấm liên tục và thiết lập hàm bão hòa nước chính xác.
  • Lãnh đạo Doanh nghiệp và Bộ phận R&D Dầu khí: Có được luận cứ khoa học tin cậy để ra quyết định đầu tư khoan thẩm lượng, tối ưu hóa thiết kế giếng khoan và giảm thiểu rủi ro cạn kiệt sớm của vỉa khí.
  • Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Sử dụng các kết quả đánh giá trữ lượng chuẩn xác làm căn cứ phê duyệt Kế hoạch Phát triển Mỏ (FDP) và Kế hoạch Thu dọn Mỏ.

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 và chuẩn hóa lý thuyết Đơn vị dòng chảy thủy lực (HFU) của Amaefule và cộng sự (1993) trong môi trường khối xây cacbonat phức tạp chịu biến đổi nứt nẻ - hòa tan thứ sinh. Luận án đã tích hợp thành công lý thuyết này với mô hình bão hòa nước phi tuyến Skelt-Harrison, giải quyết triệt để sự thất bại của phương trình điện trở truyền thống Archie vốn bị vô hiệu hóa do hệ số xi măng hóa $m$ và $n$ biến động quá lớn.

2. Đột phá phương pháp luận của luận án khi so sánh với các công trình trước đây là gì?

So với các nghiên cứu trước đây tại Việt Nam (Lê Hải An, 2016; Nguyễn Thu Huyền, 2019; Hà Quang Mẫn, 2021) vốn chỉ áp dụng các mô hình hồi quy tuyến tính đơn giản hoặc mạng ANN cục bộ, đột phá của luận án nằm ở chu trình khép kín đa thuật toán: sử dụng phân cụm không giám sát (K-means, Ward, SOM, FCM) để bóc tách 5 HFU khách quan, tiếp nối bằng mô hình Hồi quy quá trình Gaussian (GPR nhân Exponential) đạt độ tương quan vượt trội $R^2 = 0,70$, và kết thúc bằng việc liên kết vi cấu trúc dòng chảy với hàm bão hòa nước theo chiều cao SHF.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và có bằng chứng dữ liệu ra sao?

Phát hiện bất ngờ nhất là hiện tượng nghịch biến rỗng - thấm tại các giếng GK2 và GK4: các khoảng địa tầng có độ rỗng ma trận cực thấp ($\Phi < 5%$) nhưng độ thấm thực tế vượt mức $2.000$ mD. Dữ liệu lát mỏng hiển vi và thử vỉa MDT/DST đã chứng minh hiện tượng này được kiểm soát bởi các vi nứt nẻ kiến tạo (tectonic fractures) có tính kết nối cao xuyên cắt qua các tập đá vôi hạt mịn, tạo ra các kênh dẫn lưu lượng cao bất chấp ma trận đá bị nén chặt.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết toàn bộ chu trình xử lý dữ liệu từ bước tiền xử lý, chuẩn hóa đường cong ĐVLGK, công thức chuyển đổi thạch vật lý đến mã nguồn Python thực thi (được lưu trữ tại Phụ lục của luận án), cho phép các nhà nghiên cứu độc lập có thể tái lập 100% các kết quả phân cụm và huấn luyện mô hình học máy trên các tập dữ liệu giếng khoan tương tự.

5. Kế hoạch nghiên cứu 10 năm tới được phác thảo như thế nào?

Kế hoạch 10 năm tới định hình việc chuyển dịch từ mô hình thạch vật lý 1D dọc giếng sang mô hình số hóa vỉa chứa 3D/4D: kết hợp Deep Learning xử lý ảnh FMI/lát mỏng, mô phỏng mạng lưới lỗ rỗng số vi mô (Digital Rock Physics) bằng kỹ thuật Micro-CT scan 3D, và tích hợp trường địa chấn đàn hồi nhiều thành phần để đánh giá tiềm năng lưu trữ carbon ($CO_2$ sequestration) và nâng cao hệ số thu hồi khí.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Trung Dũng là một công trình nghiên cứu công phu, nghiêm túc, giải quyết trọn vẹn một trong những bài toán phức tạp nhất của địa vật lý giếng khoan và vật lý thạch học vỉa chứa dầu khí.

Năm đóng góp cốt lõi của công trình bao gồm:

  1. Xây dựng thành công chu trình khoa học khép kín tích hợp các kỹ thuật học máy tiên tiến với vật lý thạch học hiện đại để mô tả định lượng vỉa chứa đá cacbonat có tính bất đồng nhất cao.
  2. Xác lập và bóc tách chính xác 5 Đơn vị Dòng chảy Thủy lực (HFU) cho tầng chứa cacbonat tuổi Mioxen giữa hệ tầng Tri Tôn mỏ CX bằng các thuật toán phân cụm K-means và Ward's clustering.
  3. Phát triển mô hình dự báo độ thấm phi tuyến liên tục dọc thân giếng dựa trên thuật toán Hồi quy quá trình Gaussian (GPR Exponential) đạt độ tin cậy vượt trội ($R^2 = 0,70$) so với các phương pháp truyền thống.
  4. Chuẩn hóa mô hình bão hòa nước theo chiều cao (SHF) sử dụng hàm Skelt-Harrison riêng biệt cho từng HFU, loại bỏ hoàn toàn các sai số bất định do sự biến động của các tham số Archie ($m, n$).
  5. Cung cấp cơ sở dữ liệu và công cụ tính toán mở (Python code) làm tiền đề vững chắc cho công tác lập mô hình địa chất 3D, mô phỏng động lực học dòng chảy và tính toán chính xác thể tích khí tại chỗ cho mỏ CX cũng như các cấu tạo tương đồng trên thềm lục địa Việt Nam.