Tổng quan luận án

Luận án tiến sĩ với tiêu đề "Characterizing the Full In-Situ Stress Tensor and Its Applications for Petroleum Activities" (Đặc trưng hóa toàn bộ tenxơ ứng suất tại chỗ và ứng dụng trong hoạt động dầu khí) được thực hiện bởi nghiên cứu sinh Do Quang Khanh dưới sự hướng dẫn của Giáo sư Hyung-Sik Yang tại Khoa Kỹ thuật Năng lượng và Tài nguyên, Trường Sau đại học, Đại học Quốc gia Chonnam (Hàn Quốc), bảo vệ vào tháng 8 năm 2013. Hội đồng chấm luận án gồm có Tiến sĩ So-Keul Chung (Viện Địa chất và Khoáng sản Hàn Quốc - KIGAM), Giáo sư Tam Tran (CNU), Giáo sư Jeong-Hwan Lee (CNU), Giáo sư Piyush Rai (BHU, Ấn Độ) và Giáo sư Hyung-Sik Yang (CNU).

Tính cấp thiết và khoảng trống nghiên cứu

Hiểu biết về trạng thái ứng suất tại chỗ (in-situ stress) là yếu tố quyết định đối với nhiều bài toán khoa học và kỹ thuật trong địa chất, địa vật lý, khai khoáng và kỹ thuật dầu khí, bao gồm: kiến tạo mảng, dự báo động đất, ổn định công trình ngầm, độ ổn định thành giếng khoan, nứt nẻ thủy lực, mô hình thoát vỉa và dòng chảy chất lưu trong vỉa nứt nẻ tự nhiên. Trong các bể trầm tích phục vụ thăm dò và khai thác dầu khí, nhu cầu xác định tenxơ ứng suất từ các tập dữ liệu sẵn có tăng nhanh. Tuy nhiên, việc xác định đầy đủ các thành phần của tenxơ ứng suất thường gặp trở ngại do sự hạn chế của các phương pháp đo trực tiếp truyền thống (như đo nứt nẻ thủy lực hoàn chỉnh hay phương pháp khoan lõi giải tỏa ứng suất - overcoring ít khi được triển khai rộng rãi ở quy mô công nghiệp do chi phí cao). Khoảng trống lớn nhất nằm ở việc khống chế chính xác độ lớn ứng suất nén ngang cực đại ($S_{H\max}$), việc phân tích ứng suất cho các giếng khoan nghiêng bất kỳ, và việc thiếu một công cụ phần mềm tích hợp để phân tích phá hủy cơ học thành giếng nhằm đưa ra các thông số thiết kế thi công giếng tối ưu.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu của luận án được xác định gồm các nội dung sau:

  1. Xây dựng và hoàn thiện phương pháp tích hợp mới nhằm xác định toàn bộ tenxơ ứng suất tại chỗ (hướng và độ lớn của ba ứng suất chính cùng áp suất lỗ rỗng) dựa trên các tập dữ liệu thăm dò dầu khí thông thường sẵn có.
  2. Thiết kế và phát triển bộ chương trình phần mềm mang tên FAoWB (Failure Analysis of Wellbores) trên nền tảng ngôn ngữ lập trình MATLAB để tự động hóa quá trình tính toán ứng suất, phân tích phá hủy thành giếng và thiết lập các biểu đồ rủi ro (risk diagrams).
  3. Đánh giá tính chính xác và độ tin cậy của phần mềm FAoWB thông qua việc kiểm chứng chéo với các công trình nghiên cứu đã được công bố trên thế giới.
  4. Ứng dụng phương pháp tích hợp và phần mềm FAoWB vào các trường hợp nghiên cứu thực tế tại các tầng chứa móng nứt nẻ của Bồn trũng Cửu Long, Việt Nam (mỏ Bạch Hổ và mỏ X), từ đó đưa ra các khuyến nghị kỹ thuật về hướng khoan tối ưu, cửa sổ tỷ trọng dung dịch khoan và khả năng áp dụng kỹ thuật khoan thiếu áp (underbalanced drilling).

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

  • Đối tượng nghiên cứu: Trạng thái tenxơ ứng suất tại chỗ (phương và độ lớn của ứng suất thẳng đứng $S_v$, ứng suất ngang cực tiểu $S_{h\min}$, ứng suất ngang cực đại $S_{H\max}$), áp suất lỗ rỗng ($P_p$), các dạng phá hủy thành giếng (sụp lở thành giếng - breakouts - BOs, và nứt nẻ kéo do khoan tạo - drilling-induced tensile fractures - DITFs).
  • Phạm vi không gian:
    • Dữ liệu kiểm chứng quốc tế: Giếng khoan KTB (Đức) và mỏ Swan Lake (Nam Úc).
    • Dữ liệu ứng dụng thực tế: Tầng móng nứt nẻ tại mỏ Bạch Hổ (trung tâm Bồn trũng Cửu Long) và mỏ X (phía Bắc Bồn trũng Cửu Long), thềm lục địa Việt Nam.
  • Phạm vi thời gian và tài liệu: Luận án hoàn thành năm 2013, sử dụng dữ liệu địa vật lý giếng khoan, thí nghiệm cơ học đá, thí nghiệm đo áp suất ép vỡ rò rỉ (Leak-Off Test - LOT/XLOT) và ảnh giếng khoan (borehole image logs) tích lũy trong quá trình thăm dò, khai thác tại các khu vực nghiên cứu.

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

Luận án đã hệ thống hóa các nghiên cứu nền tảng và phát triển phương pháp xác định ứng suất tại chỗ qua các mốc công trình:

  • Hiện tượng sụp lở thành giếng (BOs): Được Bell và Gough (1979) mô tả lần đầu như một dạng phá hủy cắt do ứng suất nén xung quanh thành giếng. Hiện tượng này sau đó được Zoback và Zoback (1980), Plumb và Cox (1987) phát triển để xác định phương của ứng suất ngang cực đại trên thế giới.
  • Nứt nẻ kéo do khoan (DITFs): Được nhận diện qua tài liệu carota hình ảnh điện/âm môi giếng khoan bởi Aadnoy (1990b), Brudy và Zoback (1993); phương của DITFs định hướng trực giao với các vệt sụp lở BOs và trùng với phương của $S_{H\max}$.
  • Xác định ứng suất thẳng đứng ($S_v$): McGarr và Gay (1978) thiết lập phương pháp tính toán từ tích phân trọng lượng lớp đất đá phủ dựa trên đường cong carota mật độ (density log) và vận tốc sóng địa chấn kiểm tra (checkshot survey), được hoàn thiện bởi Ludwig và cộng sự.
  • Phương pháp nứt nẻ thủy lực và thí nghiệm ép vỡ rò rỉ: Haimson và Fairhurst (1967) đặt nền tảng xác định $S_{h\min}$ và $S_{H\max}$ từ đường cong áp suất - thời gian trong thí nghiệm nứt nẻ thủy lực. Hubbert và Willis (1957) chỉ ra áp suất đóng vết nứt phản ánh độ lớn của $S_{h\min}$. Breckels và van Eeklen (1982) chứng minh thí nghiệm thử rò rỉ (LOT) cung cấp ước lượng cho $S_{h\min}$. Kunze và Steiger (1991) đề xuất quy trình thử rò rỉ mở rộng (XLOT) với nhiều chu kỳ tăng/giảm áp để xác định $S_{h\min}$ với độ tin cậy cao hơn.
  • Khống chế độ lớn $S_{H\max}$: Zoback và cộng sự (1985), Moos và Zoback (1990), Barton và cộng sự (1998) nghiên cứu mối quan hệ giữa độ rộng góc sụp lở ($\theta_b$) và cường độ bền nén đơn trục ($C_0$) để tìm giới hạn dưới của $S_{H\max}$. Sibson (1974), Jaeger và Cook (1979), Byerlee (1978), Zoback và Healy (1984, 1992) phát triển lý thuyết cân bằng ma sát giới hạn trên các mặt đứt gãy thuận lợi để khống chế tỷ số ứng suất hữu hiệu tối đa, làm cơ sở xây dựng đa giác ứng suất (stress polygon).
  • Phân tích ứng suất giếng nghiêng: Aadnoy (1990), Gjonnes và cộng sự (1998), Qian và Pedersen (1991), Mastin (1988), Peska và Zoback (1995) thiết lập các phép biến đổi tọa độ tenxơ để tính toán ứng suất quanh giếng khoan có quỹ đạo không gian bất kỳ và phân tích sự xoay hướng của vết sụp lở/nứt nẻ theo góc nghiêng và góc phương vị giếng.

Khoảng trống luận án giải quyết: Các phương pháp riêng lẻ trước đây thường bộc lộ sai số lớn khi thiếu dữ liệu đo trực tiếp hoặc khi áp dụng cho các cấu trúc địa chất phức tạp có giếng khoan độ nghiêng cao. Tác giả luận án lựa chọn hướng tiếp cận tích hợp toàn diện: kết hợp dữ liệu log mật độ, thí nghiệm LOT/XLOT/thủy lực, dấu vết BOs và DITFs trên ảnh giếng khoan, lý thuyết ma sát giới hạn Byerlee, cùng việc so sánh ba chuẩn bền phá hủy đá (Mohr-Coulomb, Drucker-Prager, Mogi-Coulomb) được đóng gói trong bộ công cụ phần mềm FAoWB.

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

Cơ sở lý thuyết và khung phân tích

  1. Lý thuyết tenxơ ứng suất: Trạng thái ứng suất tại một điểm trong không gian 3 chiều được mô tả bởi tenxơ hạng hai gồm 9 thành phần ứng suất pháp và ứng suất cắt. Trong các bể trầm tích, giả định ứng suất thẳng đứng $S_v$ là một ứng suất chính, tenxơ ứng suất tại chỗ thu gọn về 4 thông số độc lập: phương của $S_{H\max}$, độ lớn $S_v$, độ lớn $S_{h\min}$ và độ lớn $S_{H\max}$.
  2. Phân loại chế độ ứng suất Anderson (1951):
    • Chế độ đứt gãy thuận (Normal Faulting - NF): $S_v > S_{H\max} > S_{h\min}$.
    • Chế độ đứt gãy trượt bằng (Strike-Slip - SS): $S_{H\max} > S_v > S_{h\min}$.
    • Chế độ đứt gãy nghịch (Reverse/Thrust Faulting - RF): $S_{H\max} > S_{h\min} > S_v$.
  3. Nguyên lý ứng suất hữu hiệu Terzaghi (1943): $S' = S - P_p$, trong đó $S'$ là tenxơ ứng suất hữu hiệu, $S$ là ứng suất tổng và $P_p$ là áp suất chất lưu lỗ rỗng đẳng hướng.
  4. Giới hạn ma sát cân bằng trượt (Byerlee, 1978; Jaeger và Cook, 1979): $$\frac{\sigma_1'}{\sigma_3'} = \frac{S_1 - P_p}{S_3 - P_p} \le \left[ \mu + \sqrt{\mu^2 + 1} \right]^2$$ Hệ số ma sát trong vỏ Trái Đất $\mu$ thường nằm trong khoảng 0.6 đến 1.0, tạo thành biên giới hạn trên đa giác ứng suất chuẩn hóa theo độ sâu.
  5. Phương trình Kirsch (1898) và biến đổi tọa độ Peska & Zoback (1995): Sử dụng ma trận quay Euler $R_s$ và $R_b$ để chuyển đổi tenxơ ứng suất từ hệ tọa độ ứng suất chính sang hệ tọa độ địa lý $(X, Y, Z)$ và hệ tọa độ giếng khoan $(x_b, y_b, z_b)$, từ đó xác định các thành phần ứng suất tiếp tuyến $\sigma_{\theta\theta}$, ứng suất trục $\sigma_{zz}$, ứng suất xuyên tâm $\sigma_{rr}$, và ứng suất cắt $\tau_{\theta z}$ trên thành giếng.
Ký hiệu đại lượng Ý nghĩa vật lý / Tên đại lượng Đơn vị / Phạm vi giá trị
$S_v$ Độ lớn ứng suất thẳng đứng (áp lực lớp phủ) $\text{MPa}$ hoặc $\text{SG}$ (tỷ trọng tương đương)
$S_{H\max}$, $S_{h\min}$ Độ lớn ứng suất nén ngang cực đại và cực tiểu $\text{MPa}$
$P_p$, $P_w$ Áp suất lỗ rỗng và áp suất chất lưu (dung dịch) giếng khoan $\text{MPa}$
$\mu$ Hệ số ma sát trong của đá $0.6 - 1.0$
$C_0$ (UCS) Cường độ bền nén đơn trục của đá $\text{MPa}$
$T$ Cường độ bền kéo của đá $\text{MPa}$ (thường quy ước $T \le 0$)
$\theta, \phi$ Góc phương vị giếng và góc nghiêng giếng Độ ($^\circ$)
$\alpha, \beta, \gamma$ Các góc quay Euler chuyển đổi hệ tọa độ Độ ($^\circ$)

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

  • Phương pháp định lượng và xử lý số liệu địa vật lý: Tích phân hàm mật độ khối theo độ sâu từ các đường carota RHOB, hiệu chỉnh DRHO và đường đo đường kính giếng CALI kết hợp số liệu vận tốc âm thanh để tính $S_v$. Phân tích các biểu đồ áp suất theo thời gian của thí nghiệm LOT, XLOT và nứt nẻ thủy lực để xác định áp suất đóng vết nứt ($P_c$), áp suất mở lại ($P_r$) và áp suất phá vỡ ($P_b$), từ đó tính $S_{h\min}$.
  • Phương pháp phân tích hình ảnh giếng khoan (Borehole Image Analysis): Nhận diện và thống kê hướng của các rãnh sụp lở BOs (trên tài liệu 4-arm dipmeter hoặc ảnh điện trở/âm học FMI/UBI) và các vết nứt nẻ kéo DITFs để xác định phương $S_{H\max}$ bằng biểu đồ hoa hồng (rose diagram) và biểu đồ tần số (histogram).
  • Phương pháp mô phỏng số và phân tích nghịch đảo: Xây dựng thuật toán tính toán thuận (forward modeling) trên MATLAB, đối chiếu điều kiện phá hủy theo ba chuẩn bền (Mohr-Coulomb, Drucker-Prager, Mogi-Coulomb) để tính toán dải giá trị cho phép của $S_{H\max}$ và áp suất dung dịch khoan yêu cầu ($P_w$).
  • Nguồn dữ liệu: Bộ số liệu từ các giếng khoan thăm dò và thẩm lượng tại tầng móng kết tinh granitoit mỏ Bạch Hổ và mỏ X thuộc Bồn trũng Cửu Long; tập số liệu công bố quốc tế tại giếng KTB (Đức) và mỏ Swan Lake (Úc).

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

CHAPTER 1: INTRODUCTION

Chương 1 trình bày cơ sở lý luận, triết lý và mục tiêu của đề tài. Tác giả chỉ ra vai trò của ứng suất tại chỗ trong các lĩnh vực địa chất, địa kỹ thuật và khai thác dầu khí. Nội dung chương tổng kết lịch sử phát triển của các kỹ thuật đo ứng suất, chỉ rõ các ưu nhược điểm của từng phương pháp riêng lẻ (như nứt nẻ thủy lực, overcoring, cơ chế tiêu điểm động đất FMS, phân tích carota hình ảnh giếng khoan). Triết lý nghiên cứu được xác lập là tích hợp và hiệu chỉnh các kỹ thuật phù hợp với tập dữ liệu thăm dò thực tế có sẵn. Chương 1 cũng trình bày bố cục tổng thể của luận án gồm 7 chương.

CHAPTER 2: IN-SITU STRESS TENSOR AND ITS RELATING CONCEPTS

Chương 2 thiết lập nền tảng lý thuyết về tenxơ ứng suất tại chỗ trong không gian ba chiều, định nghĩa các thành phần ứng suất pháp và ứng suất cắt. Tác giả phân tích trạng thái ứng suất theo lý thuyết của Anderson (1951) gắn với ba cơ chế đứt gãy tự nhiên (đứt gãy thuận, đứt gãy trượt bằng, đứt gãy nghịch). Chương này trình bày chi tiết khái niệm trạng thái ứng suất tham chiếu (lithostatic và biến dạng đơn trục - uniaxial strain), định luật ứng suất hữu hiệu Terzaghi ($S' = S - P_p$), lý thuyết giới hạn ma sát Byerlee khống chế tỷ số ứng suất hữu hiệu cực đại/cực tiểu qua hàm $f(\mu)$, và cách biểu diễn trạng thái ứng suất trên đa giác ứng suất chuẩn hóa. Cuối cùng, chương tổng hợp các mô hình vòng tròn Mohr 2D, 3D cùng đường bao phá hủy Mohr-Coulomb cho cả mặt yếu có sẵn và khối đá nguyên vẹn ($C \ne 0$).

CHAPTER 3: STRESS AND FAILURE ANALYSIS FOR WELLBORES

Chương 3 đi sâu vào bài toán tập trung ứng suất và cơ chế phá hủy xung quanh thành giếng khoan.

  • Giếng khoan thẳng đứng: Thiết lập hệ phương trình Kirsch (1898) trong hệ tọa độ trụ $(r, \theta, z)$ để xác định ứng suất xuyên tâm $\sigma_{rr}$, ứng suất vòng $\sigma_{\theta\theta}$ và ứng suất trục $\sigma_{zz}$. Dưới tác dụng của trường ứng suất nén không đẳng hướng, ứng suất vòng cực đại ($\sigma_{\theta\theta\max}$) xuất hiện tại vị trí góc $\theta = \pm 90^\circ$ (theo phương $S_{h\min}$), gây ra phá hủy nén cắt (sụp lở BOs) khi vượt quá cường độ bền đá ($C$). Ngược lại, ứng suất vòng cực tiểu ($\sigma_{\theta\theta\min}$) xuất hiện tại $\theta = 0^\circ$ (theo phương $S_{H\max}$), tạo ra nứt nẻ kéo DITFs khi $\sigma_{\theta\theta\min} \le T$.
  • Giếng khoan nghiêng bất kỳ: Trình bày giải pháp chuyển đổi tenxơ qua ba hệ tọa độ (địa lý, ứng suất và giếng khoan) theo Peska và Zoback (1995). Do xuất hiện thành phần ứng suất cắt $\tau_{\theta z}$, các ứng suất chính tại thành giếng xoay một góc $\omega$ so với trục giếng. Tác giả thiết lập biểu thức xác định các ứng suất tiếp tuyến cực đại $\sigma_{t\max}$ và cực tiểu $\sigma_{t\min}$, điều kiện hình thành BOs và DITFs trong không gian 3 chiều, và ứng dụng phương pháp chiếu bán cầu dưới (lower hemisphere projection) để biểu diễn mức độ ổn định của các giếng khoan theo mọi góc nghiêng $\phi$ và góc phương vị $\theta$.

CHAPTER 4: METHODS FOR DETERMINING IN-SITU STRESS

Chương 4 phân tích các phương pháp xác định từng thành phần của tenxơ ứng suất tại chỗ:

  1. Phương pháp đo trực tiếp: Nứt nẻ thủy lực (thiết bị packer, máy bơm cao áp, cảm biến áp suất), kỹ thuật overcoring, cơ chế tiêu điểm động đất (Earth focal mechanism - FMS).
  2. Quy trình tích hợp mới từ dữ liệu dầu khí:
    • Xác định phương $S_{H\max}$: Phân tích thống kê vệt sụp lở BOs và nứt nẻ DITFs trên tài liệu carota 4 cần (four-arm dipmeter) và carota hình ảnh (FMI/UBI).
    • Xác định độ lớn $S_v$: Tích phân số liệu carota mật độ khối (RHOB, DRHO, CALI) theo độ sâu, bổ sung phần nông bị thiếu bằng số liệu vận tốc âm thanh checkshot survey.
    • Xác định độ lớn $S_{h\min}$: Phân tích các điểm gãy áp suất trên biểu đồ áp suất - thời gian hoặc áp suất - căn bậc hai thời gian của các thí nghiệm nứt nẻ thủy lực, XLOT và LOT.
    • Khống chế độ lớn $S_{H\max}$: Sử dụng kết hợp giới hạn ma sát Byerlee trên đa giác ứng suất và mô hình tính toán thuận trạng thái phá hủy thành giếng (sự xuất hiện hoặc không xuất hiện của BOs/DITFs kết hợp với độ bền nén đơn trục $C_0$ và áp suất dung dịch khoan $P_w$).
    • Xác định áp suất lỗ rỗng $P_p$: Từ các phép đo trực tiếp RFT/MDT hoặc tính toán từ log địa vật lý.

CHAPTER 5: MODEL DEVELOPMENT FOR FAILURE ANALYSIS OF WELLBORE (FAOWB)

Chương 5 giới thiệu việc thiết kế, xây dựng và kiểm chứng phần mềm FAoWB được viết bằng ngôn ngữ MATLAB.

  • Cấu trúc phần mềm: Giao diện đồ họa người dùng trực quan gồm màn hình khởi động, màn hình chính, hệ thống menu (File, Input Data, Failure Criteria, Process, Output) và các tab nhập liệu thông số (Description, Stress, Rock properties, Well). Phần mềm cho phép tích hợp các chuẩn bền: Mohr-Coulomb, Drucker-Prager và Mogi-Coulomb.
  • Kiểm chứng tính chính xác qua hai nghiên cứu điển hình:
    • Trường hợp 1 (Case 1): Đối chiếu với nghiên cứu của Barton (1998) tại giếng khoan sâu KTB (Đức) về phân tích phá hủy nén và độ rộng sụp lở thành giếng. Kết quả tính toán phân bố ứng suất, biểu đồ đa giác ứng suất và tỷ trọng dung dịch yêu cầu từ FAoWB hoàn toàn trùng khớp với kết quả đã công bố của Barton.
    • Trường hợp 2 (Case 2): Đối chiếu với nghiên cứu của Meyer (2002) về độ ổn định giếng khoan tại mỏ Swan Lake (Nam Úc). Kết quả biểu đồ rủi ro phá hủy và phân bố ứng suất của FAoWB cho độ tương thích cao, khẳng định tính đúng đắn và độ tin cậy của thuật toán.

CHAPTER 6: CASE STUDIES AND IMPLICATIONS

Chương 6 là phần ứng dụng thực tế của phương pháp tích hợp và phần mềm FAoWB vào hai mỏ dầu khí thuộc Bồn trũng Cửu Long, Việt Nam.

  • Đặc điểm địa chất Bồn trũng Cửu Long: Cấu trúc địa tầng, hệ thống đứt gãy nứt nẻ chính và đặc điểm tầng móng granitoit phong hóa/nứt nẻ.
  • Mỏ Bạch Hổ (White Tiger Field - trung tâm bồn trũng):
    • Phân tích carota hình ảnh tại các khoảng thân giếng trong móng cho thấy sự xuất hiện rõ nét của cả BOs và DITFs. Biểu đồ hoa hồng xác định phương của $S_{H\max}$ định hướng theo phương Tây Bắc - Đông Nam (NW-SE).
    • Xác định đường ứng suất thẳng đứng $S_v$ từ carota mật độ, xác định $S_{h\min}$ từ số liệu nứt nẻ thủy lực và LOT, xác định $P_p$.
    • Sử dụng đa giác ứng suất và phần mềm FAoWB để khống chế dải giá trị của $S_{H\max}$ tại các độ sâu móng 3900 m, 4100 m, 4300 m và 4500 m.
    • Thiết lập biểu đồ phân bố ứng suất và biểu đồ rủi ro (risk diagrams) theo ba chuẩn bền. Kết quả chỉ ra chế độ ứng suất tại móng mỏ Bạch Hổ chuyển dịch giữa chế độ đứt gãy thuận (NF) và trượt bằng (SS).
  • Mỏ X (X Field - phía Bắc bồn trũng):
    • Phân tích số liệu DITFs tại giếng X1 xác định phương $S_{H\max}$ theo hướng NW-SE.
    • Xác định $S_v$, $S_{h\min}$ từ LOTs/FITs và áp suất lỗ rỗng $P_p$.
    • Tính toán tenxơ ứng suất đầy đủ tại độ sâu móng 2300 m; thiết lập biểu đồ rủi ro và đánh giá ổn định cho hai giếng khoan nghiêng thực tế.
    • Đánh giá khả năng áp dụng kỹ thuật khoan thiếu áp ($P_w = 22\text{ MPa}$) tại mỏ X.
  • Ý nghĩa kỹ thuật: Đề xuất quỹ đạo giếng khoan tối ưu định hướng theo phương Đông Bắc - Tây Nam (NE-SW) để hạn chế tối đa rủi ro mất ổn định thành giếng (sụp lở BOs hoặc mất tuần hoàn dung dịch do DITFs), và xác định khoảng tỷ trọng dung dịch khoan an toàn.

CHAPTER 7: CONCLUSIONS AND RECOMMENDATIONS

Chương 7 tổng kết các kết luận khoa học rút ra từ toàn bộ công trình nghiên cứu và đưa ra các khuyến nghị cho công tác thiết kế, thi công giếng khoan và nghiên cứu tiếp theo trong ngành dầu khí.

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

Đóng góp mới về mặt khoa học và phương pháp luận

  1. Luận án đã xây dựng một phương pháp tích hợp hoàn chỉnh cho phép xác định đầy đủ 4 thành phần của tenxơ ứng suất tại chỗ ($S_v$, $S_{h\min}$, $S_{H\max}$ và phương $S_{H\max}$) cùng áp suất lỗ rỗng ($P_p$) từ các tập dữ liệu thăm dò dầu khí sẵn có, giải quyết bài toán khó khăn nhất là khống chế độ lớn ứng suất nén ngang cực đại ($S_{H\max}$).
  2. Phát triển bộ công cụ phần mềm FAoWB bằng ngôn ngữ MATLAB với cấu trúc module hóa, tích hợp các thuật toán biến đổi tọa độ 3D và ba tiêu chuẩn phá hủy đá (Mohr-Coulomb, Drucker-Prager, Mogi-Coulomb), giúp tự động hóa quá trình phân tích ổn định thành giếng và đánh giá rủi ro địa cơ học.
  3. Cung cấp dữ liệu đối chuẩn thực nghiệm có độ tin cậy cao qua việc tái lập chính xác kết quả nghiên cứu quốc tế tại giếng KTB (Đức) và mỏ Swan Lake (Úc).

Đóng góp mới về mặt thực tiễn ứng dụng

  1. Xác định trạng thái ứng suất tầng móng Bồn trũng Cửu Long:
    • Phương của ứng suất ngang cực đại $S_{H\max}$ tại mỏ Bạch Hổ và mỏ X được xác định nhất quán theo phương Tây Bắc - Đông Nam (NW-SE).
    • Chế độ ứng suất đặc trưng trong tầng móng kết tinh granitoit là chế độ đứt gãy thuận (NF: $S_v > S_{H\max} > S_{h\min}$) hoặc trượt bằng (SS: $S_{H\max} > S_v > S_{h\min}$).
  2. Tối ưu hóa thiết kế thi công giếng khoan:
    • Xác lập cơ sở khoa học khuyến nghị quỹ đạo thân giếng tối ưu nên khoan theo hướng Đông Bắc - Tây Nam (NE-SW) - vuông góc với phương $S_{H\max}$ hoặc theo góc nghiêng phù hợp để giảm thiểu nguy cơ sập lở thành giếng và nứt nẻ mất tuần hoàn.
    • Thiết lập cửa sổ tỷ trọng dung dịch khoan an toàn nhằm ngăn ngừa biến dạng dẻo/sụp đổ thành giếng do thiếu áp hoặc gây nứt vỉa mất dung dịch do thừa áp.
    • Đánh giá định lượng tính khả thi của kỹ thuật khoan thiếu áp (UBD) tại mỏ X dưới áp suất kiểm soát $P_w = 22\text{ MPa}$.
Khu vực nghiên cứu Vị trí địa lý trong bồn trũng Khoảng độ sâu phân tích móng Phương $S_{H\max}$ xác định Chế độ ứng suất chủ đạo
Mỏ Bạch Hổ (White Tiger) Trung tâm Bồn trũng Cửu Long $3900\text{ m} - 4500\text{ m}$ Tây Bắc - Đông Nam (NW-SE) Đứt gãy thuận (NF) đến Trượt bằng (SS)
Mỏ X Phía Bắc Bồn trũng Cửu Long $2300\text{ m}$ Tây Bắc - Đông Nam (NW-SE) Đứt gãy thuận (NF) đến Trượt bằng (SS)
Giếng KTB (Kiểm chứng) Đức Thân giếng sâu Xác định theo dữ liệu Barton NF / SS
Mỏ Swan Lake (Kiểm chứng) Nam Úc Các giếng nghiêng Xác định theo dữ liệu Meyer Phân tích ổn định thành giếng

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

Dựa trên phạm vi và các giả định được tác giả trình bày trong luận án:

  • Giả định đàn hồi tuyến tính và môi trường liên tục: Mô hình giải tích Kirsch và các phép biến đổi tenxơ trong FAoWB giả định đá có tính chất đàn hồi, đồng nhất và đẳng hướng. Thực tế đá móng nứt nẻ có tính dị hướng và phi tuyến cao, cần được tích hợp thêm các mô hình cơ học phá hủy phi liên tục.
  • Giả định độ bền kéo của khối đá: Do độ bền kéo của đá thường nhỏ và khối đá chứa nhiều mặt yếu vi nứt nẻ, nghiên cứu giả định độ bền kéo hữu hiệu $T \approx 0$ khi phân tích nứt nẻ DITFs.
  • Đặc tính thấm của thành giếng: Các phương trình cơ bản giả định thành giếng không thấm nhờ hình thành lớp vỏ sét (mud cake) khi khoan cân bằng áp hoặc quá áp; đối với chế độ khoan thiếu áp (UBD), sự xâm nhập chất lưu đòi hỏi phải giả định áp suất lỗ rỗng cục bộ vùng cận đáy giếng.
  • Hướng nghiên cứu tiếp: Tích hợp các module phân tích cơ địa tầng 3D/4D, mở rộng phần mềm FAoWB để mô phỏng tương tác nhiệt - thủy - cơ (THM) trong quá trình bơm ép nứt nẻ thủy lực và khai thác nhiệt đa tầng nứt nẻ.

Giá trị tham khảo

Luận án là tài liệu tham khảo cho:

  • Nghiên cứu sinh và giảng viên chuyên ngành Địa kỹ thuật, Địa cơ học mỏ, Kỹ thuật Dầu khí: Cung cấp khung lý thuyết hoàn chỉnh về cơ chế tập trung ứng suất quanh giếng khoan nghiêng và phương pháp luận tích hợp số liệu carota với cơ học đá.
  • Kỹ sư thiết kế khoan và địa chất công trình dầu khí: Tham khảo quy trình xác định cửa sổ tỷ trọng dung dịch khoan, lựa chọn quỹ đạo giếng khoan tối ưu tại Bồn trũng Cửu Long và các bể trầm tích có điều kiện địa chất tương đồng nhằm hạn chế sự cố sập thành giếng hoặc mất tuần hoàn dung dịch.
  • Các nhà phát triển phần mềm mô phỏng địa cơ học: Tham khảo kiến trúc thiết kế thuật toán, cấu trúc dữ liệu và giao diện tính toán của gói phần mềm FAoWB trên MATLAB.

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

1. Tại sao việc khống chế độ lớn ứng suất nén ngang cực đại ($S_{H\max}$) là khâu khó khăn nhất?

Trong khi ứng suất thẳng đứng $S_v$ có thể tính toán từ tích phân carota mật độ, và ứng suất ngang cực tiểu $S_{h\min}$ có thể đo trực tiếp qua áp suất đóng vết nứt của thí nghiệm nứt nẻ thủy lực hoặc XLOT, thì không có phương pháp đo trực tiếp đơn giản nào để xác định $S_{H\max}$ ở độ sâu lớn. $S_{H\max}$ bắt buộc phải được khống chế gián tiếp thông qua sự kết hợp giữa lý thuyết cân bằng ma sát Byerlee (giới hạn trên của đa giác ứng suất), sự xuất hiện hay không xuất hiện của các vết phá hủy thành giếng (BOs, DITFs) và cường độ bền nén của đá ($C_0$).

2. Sự khác biệt cơ bản giữa cơ chế hình thành rãnh sụp lở (BOs) và nứt nẻ kéo do khoan (DITFs) là gì?

  • Rãnh sụp lở (BOs) là hiện tượng phá hủy nén cắt xảy ra khi ứng suất vòng tiếp tuyến cực đại ($\sigma_{\theta\theta\max}$) vượt quá cường độ bền nén của đá, tập trung tại vị trí góc phương vị của $S_{h\min}$ ($\theta = \pm 90^\circ$ đối với giếng đứng).
  • Nứt nẻ kéo do khoan (DITFs) là hiện tượng phá hủy kéo xảy ra khi ứng suất vòng tiếp tuyến cực tiểu ($\sigma_{\theta\theta\min}$) giảm xuống thấp hơn cường độ bền kéo của đá ($T \le 0$), tập trung tại vị trí góc phương vị của $S_{H\max}$ ($\theta = 0^\circ$ đối với giếng đứng), có hướng vuông góc với vệt BOs.

3. Phần mềm FAoWB được xây dựng trên ngôn ngữ nào và được kiểm chứng bằng những nguồn dữ liệu nào?

Phần mềm FAoWB (Failure Analysis of Wellbores) được lập trình trên nền tảng MATLAB. Độ tin cậy của phần mềm được kiểm chứng qua hai công trình công bố quốc tế: trường hợp nghiên cứu của Barton (1998) tại giếng khoan sâu KTB (Đức) về phân tích độ rộng sụp lở thành giếng, và trường hợp nghiên cứu của Meyer (2002) tại mỏ Swan Lake (Nam Úc) về độ ổn định giếng khoan.

4. Kết quả nghiên cứu ứng suất tại Bồn trũng Cửu Long đưa ra khuyến nghị gì cho công tác khoan giếng?

Kết quả nghiên cứu chỉ ra hướng $S_{H\max}$ trong tầng móng Bồn trũng Cửu Long (mỏ Bạch Hổ và mỏ X) định hướng theo phương NW-SE, với chế độ ứng suất đứt gãy thuận (NF) đến trượt bằng (SS). Luận án khuyến nghị thiết kế quỹ đạo giếng khoan tối ưu theo hướng Đông Bắc - Tây Nam (NE-SW) để đạt độ ổn định cơ học thành giếng cao nhất, đồng thời xác định khoảng dung dịch khoan an toàn nhằm tránh sự cố sụp lở và mất tuần hoàn dung dịch.

Kết luận

Luận án của tác giả Do Quang Khanh đã hệ thống hóa và phát triển phương pháp tích hợp xác định toàn bộ tenxơ ứng suất tại chỗ từ các tập dữ liệu thăm dò dầu khí thông thường, kết hợp xây dựng gói phần mềm FAoWB trên nền tảng MATLAB. Các thuật toán và mô hình phân tích phá hủy giếng khoan đã được kiểm chứng độ chính xác qua các dữ liệu thực tế tại Đức và Úc, đồng thời áp dụng vào tầng móng nứt nẻ tại mỏ Bạch Hổ và mỏ X thuộc Bồn trũng Cửu Long. Kết quả nghiên cứu cung cấp cơ sở khoa học cho việc định hướng quỹ đạo khoan tối ưu theo phương NE-SW, xác lập cửa sổ áp suất dung dịch khoan an toàn và phục vụ phân tích ổn định địa cơ học cho các hoạt động khai thác dầu khí.