Tổng quan nghiên cứu

Trong ngành công nghiệp khai thác dầu khí toàn cầu, sự cố mất ổn định thành giếng khoan chiếm khoảng 10% đến 15% tổng chi phí khoan, gây thiệt hại hàng trăm triệu USD mỗi năm do các vấn đề kẹt cần khoan, sập thành giếng và mất dung dịch. Vấn đề này càng trở nên nghiêm trọng khi hoạt động thăm dò, khai thác mở rộng sang các vùng nước sâu, các vỉa có nhiệt độ và áp suất cao, hoặc khi thi công các giếng khoan định hướng, giếng khoan ngang có góc nghiêng lớn. Bên cạnh đó, hiện tượng vỡ vụn đất đá thành hệ dẫn đến cát xâm nhập vào dòng chất lưu khai thác làm xói mòn thiết bị, phá hủy ống chống và suy giảm lưu lượng giếng nghiêm trọng.

Nghiên cứu tập trung giải quyết bài toán cốt lõi: mô hình hóa trạng thái ứng suất xung quanh thành giếng khoan nghiêng trong không gian ba chiều và ứng dụng các tiêu chuẩn phá hủy đất đá tiên tiến nhằm tối ưu hóa quỹ đạo khoan cũng như dự báo rủi ro xuất hiện cát. Mục tiêu cụ thể là xây dựng công cụ tính toán giải tích tự động, phân tích sự phân bố ứng suất và đánh giá độ chính xác của các tiêu chuẩn bền truyền thống so với các mô hình phi tuyến hiện đại.

Phạm vi nghiên cứu được thực hiện trên cơ sở dữ liệu địa cơ học tại các độ sâu từ 4.800 ft đến 8.530 ft, bao gồm dữ liệu thực địa từ mỏ Cyrus (thềm lục địa Vương quốc Anh) và mỏ Pagerungan (Indonesia), đồng thời định hướng ứng dụng cho các đối tượng mỏ tại thềm lục địa Việt Nam. Ý nghĩa thực tiễn của công trình thể hiện qua việc giảm thiểu sai số tính toán tỷ trọng dung dịch khoan từ 0,5 đến 0,9 lb/gal, giúp giảm thiểu rủi ro sự cố địa cơ học thành giếng tới 25% đến 30% và kéo dài tuổi thọ khai thác an toàn cho các giếng dầu khí.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết đàn hồi tuyến tính và cơ học đất đá dầu khí, mô tả sự tập trung ứng suất quanh lỗ rỗng hình trụ thông qua lời giải giải tích Kirsch mở rộng trong hệ tọa độ trụ ba chiều (r, $\theta$, z). Trạng thái ứng suất tại chỗ ban đầu bao gồm ba thành phần ứng suất chính trực giao: ứng suất đứng do trọng lượng lớp phủ ($\sigma_v$), ứng suất ngang cực đại ($\sigma_H$) và ứng suất ngang cực tiểu ($\sigma_h$).

Hệ thống lý thuyết tích hợp phân tích ba chế độ ứng suất kiến tạo chính: chế độ đứt gãy thuận ($\sigma_v > \sigma_H > \sigma_h$), chế độ đứt gãy trượt bằng ($\sigma_H > \sigma_v > \sigma_h$) và chế độ đứt gãy nghịch ($\sigma_H > \sigma_h > \sigma_v$). Ba khái niệm cốt lõi được sử dụng để phân tích ranh giới phá hủy dẻo và nứt nẻ của đất đá bao gồm: ứng suất pháp bát diện ($\sigma_{oct}$), ứng suất tiếp bát diện ($\tau_{oct}$) và ứng suất trung bình hiệu dụng phẳng ($\sigma_{m,2}$).

Mô hình nghiên cứu tiến hành đối chiếu và kiểm chứng ba tiêu chuẩn phá hủy đất đá tiêu biểu:

  1. Tiêu chuẩn Mohr-Coulomb: Tiêu chuẩn tuyến tính kinh điển dựa trên lực dính kết ($c$) và góc ma sát trong ($\phi$), giả định ứng suất chính trung gian ($\sigma_2$) không ảnh hưởng đến sự phá hủy.
  2. Tiêu chuẩn Mogi-Coulomb: Tiêu chuẩn ba trục mở rộng xem xét đầy đủ ảnh hưởng gia cường của ứng suất chính trung gian thông qua mối quan hệ tuyến tính giữa $\tau_{oct}$ và $\sigma_{m,2}$.
  3. Tiêu chuẩn Lade sửa đổi (Modified Lade): Tiêu chuẩn phi tuyến tính ba chiều sử dụng các bất biến ứng suất bậc một và bậc ba ($I_1, I_3$) để dự báo áp suất giới hạn của dung dịch khoan.

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

Nghiên cứu sử dụng phương pháp mô hình hóa giải tích kết hợp lập trình thuật toán số trên ngôn ngữ FORTRAN để giải quyết hệ phương phân bố ứng suất và quét không gian góc nghiêng ($0^\circ - 90^\circ$) cùng góc phương vị ($0^\circ - 360^\circ$). Nguồn dữ liệu đầu vào được thu thập từ các thí nghiệm nén ba trục xác định độ bền đá cát với lực cố kết $c = 740\text{ psi}$ đến $1.800\text{ psi}$, góc ma sát trong $\phi = 33^\circ$ đến $43,8^\circ$, hệ số Poisson $\nu = 0,22$ đến $0,25$, kết hợp cùng kết quả thử nghiệm thủy lực đo áp suất nứt vỡ vỉa (Leak-off test).

Cỡ mẫu phân tích bao gồm 3 kịch bản chế độ ứng suất kiến tạo chuẩn hóa tại độ sâu 8.000 ft và 2 tập dữ liệu kiểm chứng mỏ thực tế ngoài khơi với độ sâu vận hành từ 4.800 ft đến 8.530 ft. Phương pháp chọn mẫu tập trung vào các giếng khoan có đầy đủ nhật ký địa vật lý giếng, lịch sử dung dịch khoan thực tế và báo cáo ổn định thành giếng để đảm bảo tính xác thực khi kiểm định. Lý do lựa chọn phương pháp phân tích giải tích kết hợp chương trình máy tính tự lập là khả năng tính toán nhanh, chính xác cao, không đòi hỏi tài nguyên tính toán phức tạp như các phương pháp phần tử hữu hạn nhưng vẫn đảm bảo quét toàn diện hàng nghìn tổ hợp góc quỹ đạo trong thời gian thực hiện nghiên cứu kéo dài 12 tháng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán và kiểm định trên mô hình đã làm sáng tỏ 3 phát hiện khoa học quan trọng:

Thứ nhất, quỹ đạo giếng tối ưu để duy trì ổn định cơ học phụ thuộc chặt chẽ vào chế độ ứng suất kiến tạo tại chỗ thay vì cố định ở một dạng hình học:

  • Trong chế độ ứng suất thuận (NF): Giếng đứng ổn định hơn giếng ngang, tuy nhiên quỹ đạo tối ưu nhất là giếng nghiêng $40^\circ$ theo phương song song với ứng suất ngang cực tiểu $\sigma_h$.
  • Trong chế độ trượt bằng (SS): Giếng khoan ngang có độ ổn định vượt trội hơn giếng đứng và giếng nghiêng trên toàn bộ các phương vị; hướng khoan tối ưu đạt góc phương vị $30^\circ$.
  • Trong chế độ ứng suất nghịch (RF): Giếng khoan theo phương ứng suất ngang cực đại $\sigma_H$ (góc phương vị $0^\circ$) với góc nghiêng $50^\circ$ mang lại độ ổn định cao nhất, giảm thiểu tối đa áp suất trên cân bằng cần thiết.

Thứ hai, có sự chênh lệch rõ rệt giữa các tiêu chuẩn phá hủy đất đá khi tính toán áp suất giếng khoan tối thiểu. Tiêu chuẩn Mohr-Coulomb luôn cho kết quả thận trọng nhất, đưa ra giá trị áp suất dung dịch yêu cầu cao hơn từ 8,5% đến 15% so với thực tế. Tiêu chuẩn Lade sửa đổi đưa ra dự báo lạc quan nhất, trong khi tiêu chuẩn Mogi-Coulomb thể hiện độ tương thích chính xác nhất với dữ liệu khoan thực tế.

Thứ ba, kết quả kiểm chứng tại giếng PGA-2 thuộc mỏ Pagerungan ở độ sâu 6.000 ft cho thấy: Tiêu chuẩn Mohr-Coulomb dự báo giếng sẽ bị sụp đổ nếu tỷ trọng dung dịch nhỏ hơn 11,45 lb/gal (cao hơn 0,9 lb/gal so với thực tế). Tuy nhiên, trên thực tế giếng đã được thi công an toàn tuyệt đối với tỷ trọng dung dịch chỉ 10,55 lb/gal. Mô hình Mogi-Coulomb đã tái hiện chính xác ngưỡng ổn định này ở mức 10,48 lb/gal. Tương tự, tại mỏ Cyrus ở độ sâu 8.530 ft, mô hình Mogi-Coulomb xác nhận giếng ngang hoàn toàn ổn định với dung dịch khoan 9,6 lb/gal, trong khi Mohr-Coulomb dự báo sai lệch rằng giếng sẽ mất ổn định hoàn toàn.

Thảo luận kết quả

Sự sai khác căn bản giữa các tiêu chuẩn bắt nguồn từ việc xem xét hay bỏ qua ứng suất chính trung gian ($\sigma_2$). Tiêu chuẩn Mohr-Coulomb cổ điển giả định $\sigma_2 = \sigma_3$, vô tình loại bỏ hiệu ứng nén giữ của đất đá trong không gian ba chiều, dẫn đến việc đánh giá thấp độ bền thực tế của thành hệ và tính toán thừa trọng lượng dung dịch khoan. Khi tỷ trọng dung dịch bị đẩy lên quá cao một cách không cần thiết, rủi ro nứt vỡ vỉa và mất dung dịch vào tầng chứa sẽ gia tăng đáng kể.

Ngược lại, tiêu chuẩn Mogi-Coulomb và Lade sửa đổi đã tích hợp đại lượng $\sigma_{m,2}$ và $\tau_{oct}$, phản ánh đúng bản chất biến dạng dẻo đa trục của đất đá dưới lòng đất. Trong báo cáo kỹ thuật, dữ liệu phân tích này có thể được biểu diễn trực quan qua đồ thị phân bố ứng suất cực (Polar Plots) thể hiện áp suất trên cân bằng tối thiểu theo góc nghiêng ($0^\circ - 90^\circ$) và phương vị ($0^\circ - 360^\circ$). Bảng so sánh cửa sổ tỷ trọng dung dịch giữa các mô hình sẽ giúp các kỹ sư công nghệ nhận diện ngay biên an toàn cơ học, ngăn ngừa hiện tượng phá hủy cắt gây sụp thành lỗ khoan hoặc vỡ vụn sinh cát khi hạ áp suất đáy giếng trong khai thác.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm tối ưu hóa công tác khoan và khai thác dầu khí:

  1. Chuyển đổi tiêu chuẩn thiết kế tỷ trọng dung dịch khoan: Các bộ phận công nghệ mỏ và khoan cần thay thế tiêu chuẩn Mohr-Coulomb truyền thống bằng tiêu chuẩn Mogi-Coulomb trong tính toán cửa sổ áp suất giếng khoan. Mục tiêu là hiệu chỉnh giảm tỷ trọng dung dịch khoan từ 0,4 đến 0,8 lb/gal, giúp tiết kiệm khoảng 15% chi phí hóa chất dung dịch và hạn chế tối đa nguy cơ ô nhiễm vỉa chứa. Thời gian triển khai áp dụng trong vòng 6 tháng kể từ giai đoạn lập báo cáo thiết kế kỹ thuật giếng.
  2. Thiết kế quỹ đạo giếng bám sát trường ứng suất kiến tạo: Phòng công nghệ khoan cần khảo sát chính xác tensor ứng suất tại chỗ trước khi định hướng giếng. Đối với vùng có cơ chế ứng suất thuận, ưu tiên thiết kế góc nghiêng $40^\circ$ theo hướng $\sigma_h$; đối với vỉa trượt bằng, ưu tiên giếng ngang góc phương vị $30^\circ$. Giải pháp này giúp giảm thiểu 20% đến 30% nguy cơ kẹt cần và sạt lở giếng trong suốt chu kỳ khoan kéo dài từ 30 đến 45 ngày mỗi giếng.
  3. Kiểm soát áp suất dòng đáy giếng để chống xuất hiện cát: Đơn vị khai thác cần thiết lập giới hạn áp suất chênh lệch (drawdown pressure) không vượt quá giới hạn bền cắt của đá thành hệ (duy trì biên độ an toàn áp suất đáy giếng từ 200 đến 300 psi trên ngưỡng tới hạn). Biện pháp này giúp ngăn chặn 90% nguy cơ phá hủy vỉa sinh cát, kéo dài tuổi thọ hệ thống thiết bị khai thác thêm 3 đến 5 năm mà không cần lắp đặt sớm các cụm thiết bị lọc cát tốn kém.
  4. Chuẩn hóa và phát triển phần mềm địa cơ học chuyên dụng: Các viện nghiên cứu và doanh nghiệp dầu khí cần đầu tư nâng cấp chương trình FORTRAN trong luận văn thành phần mềm có giao diện đồ họa thân thiện trong lộ trình 18 tháng, tích hợp trực tiếp vào quy trình thiết kế giếng khoan tại các bể trầm tích Cửu Long và Nam Côn Sơn.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

  1. Kỹ sư khoan và hoàn thiện giếng: Ứng dụng mô hình để xác định chính xác cửa sổ tỷ trọng dung dịch khoan an toàn, lựa chọn chiều sâu đặt ống chống hợp lý và lập quỹ đạo giếng tối ưu nhằm triệt tiêu rủi ro kẹt bộ khoan cụ.
  2. Kỹ sư địa cơ học và quản lý mỏ dầu khí: Tham khảo quy trình phân tích tensor ứng suất 3D, từ đó dự báo chính xác ngưỡng giảm áp suất đáy giếng cho phép nhằm ngăn chặn hiện tượng cát xâm nhập trong quá trình khai thác thứ cấp.
  3. Các nhà nghiên cứu và chuyên gia cơ học kỹ thuật: Khai thác phương pháp giải tích mở rộng, thuật toán lập trình chuyển đổi hệ tọa độ trụ và cơ sở lý thuyết so sánh các tiêu chuẩn phá hủy đa trục (Mohr-Coulomb, Mogi-Coulomb, Modified Lade) để phục vụ các đề tài nghiên cứu chuyên sâu.
  4. Giảng viên, học viên cao học và sinh viên ngành Kỹ thuật Dầu khí / Cơ học: Sử dụng toàn bộ tài liệu, sơ đồ khối thuật toán và mã nguồn chương trình tính toán như một giáo trình tham khảo chuẩn mực về địa cơ học ứng dụng trong công nghiệp năng lượng.

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

Tại sao tiêu chuẩn Mohr-Coulomb truyền thống lại dự báo sai lệch tỷ trọng dung dịch khoan?

Tiêu chuẩn Mohr-Coulomb giả thiết ứng suất chính trung gian không ảnh hưởng đến độ bền đất đá, dẫn đến việc đánh giá thấp khả năng chịu lực thực tế. Tại mỏ Pagerungan, mô hình này dự báo cần tỷ trọng 11,45 lb/gal, trong khi thực tế giếng khoan ổn định hoàn toàn ở mức 10,55 lb/gal, gây ra độ lệch lên tới 0,9 lb/gal.

Việc lựa chọn quỹ đạo giếng khoan tối ưu phụ thuộc vào yếu tố nào?

Quỹ đạo tối ưu phụ thuộc trực tiếp vào trạng thái ứng suất kiến tạo tại chỗ gồm ba chế độ: thuận, trượt bằng và nghịch. Ví dụ, trong chế độ ứng suất thuận ở độ sâu 8.000 ft, giếng nghiêng $40^\circ$ theo hướng ứng suất ngang nhỏ nhất mang lại độ ổn định cao hơn đáng kể so với giếng thẳng đứng hoặc giếng ngang.

Làm thế nào để mô hình hỗ trợ kiểm soát hiện tượng xuất hiện cát khi khai thác?

Mô hình tính toán sự tập trung ứng suất quanh lỗ bắn mở vỉa khi áp suất đáy giếng suy giảm. Bằng cách so sánh ứng suất tiếp xúc với độ bền cắt của vỉa đá cát (lực dính kết khoảng 740 psi), kỹ sư xác định được áp suất chênh lệch tới hạn, từ đó duy trì lưu lượng khai thác an toàn mà không làm sụp đổ cấu trúc đá.

Tiêu chuẩn Mogi-Coulomb có ưu thế gì nổi bật so với tiêu chuẩn Lade sửa đổi?

Tiêu chuẩn Mogi-Coulomb vừa tính đến ảnh hưởng gia cường của ứng suất chính trung gian, vừa giữ nguyên các thông số độ bền vật lý quen thuộc là lực cố kết $c$ và góc ma sát $\phi$. Điều này giúp kết quả tính toán có tính ổn định cao, tránh các dự báo quá lạc quan như tiêu chuẩn Lade sửa đổi.

Thuật toán trong luận văn có thể áp dụng trực tiếp cho các mỏ tại Việt Nam không?

Hoàn toàn có thể áp dụng trực tiếp. Chương trình cho phép nhập các thông số đầu vào đặc thù của thềm lục địa Việt Nam như độ sâu từ 4.000 đến 12.000 ft, hệ số Poisson từ 0,20 đến 0,30 và áp lực lỗ rỗng thực tế để xuất ra đồ thị cửa sổ dung dịch và quỹ đạo tối ưu chỉ trong vài giây.

Kết luận

  • Luận văn đã xây dựng thành công mô hình giải tích phân bố ứng suất 3D quanh thành giếng khoan nghiêng và phát triển chương trình máy tính chuyên dụng.
  • Đánh giá toàn diện và định lượng hóa sự sai khác giữa các tiêu chuẩn phá hủy, chứng minh tiêu chuẩn Mogi-Coulomb phản ánh chính xác nhất điều kiện làm việc thực tế của đất đá dưới lòng giếng khoan.
  • Xác lập quy luật lựa chọn quỹ đạo giếng tối ưu tương ứng với 3 chế độ ứng suất kiến tạo, chỉ ra góc nghiêng tối ưu $40^\circ$ trong chế độ thuận và $50^\circ$ trong chế độ nghịch.
  • Cung cấp phương pháp luận khoa học và công cụ tính toán định lượng phục vụ dự báo nguy cơ xuất hiện cát trong khai thác dầu khí.
  • Hướng phát triển tiếp theo trong 12 đến 24 tháng tới là mở rộng mô hình kết hợp hiệu ứng nhiệt - thủy lực đàn hồi (thermo-poroelasticity) và phi tuyến dẻo cho đất đá yếu.

Các kỹ sư công nghệ và nhà quản lý mỏ nên chủ động áp dụng mô hình phân tích địa cơ học đa trục Mogi-Coulomb vào quy trình thiết kế giếng hiện hành để tối ưu hóa chi phí và đảm bảo an toàn tuyệt đối cho các chiến dịch khoan thăm dò, khai thác dầu khí.