Bối cảnh và vấn đề nghiên cứu

Trong kỹ thuật khai thác dầu khí, việc phân tích động học dòng chảy từ vỉa vào đáy giếng và từ đáy giếng lên hệ thống thu gom bề mặt đóng vai trò quyết định đến sản lượng và tuổi thọ công trình khai thác. Khi năng lượng tự nhiên của vỉa suy giảm hoặc xuất hiện các hiện tượng ngập lỏng, sụt giảm áp suất và cát xâm nhập, việc mô phỏng chính xác hệ thống giếng cùng với việc lựa chọn giải pháp khai thác nhân tạo (Artificial Lift) trở thành nhiệm vụ kỹ thuật bắt buộc.

Báo cáo bài tập cuối khóa môn Công nghệ Khai thác Dầu khí với đề tài "Lựa Chọn Hệ Thống Artificial Lift Tối Ưu Trong Mô Phỏng" được thực hiện bởi nhóm sinh viên lớp DC19DK1 (Khoa Kỹ thuật Địa chất & Dầu khí, Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh) dưới sự hướng dẫn của TS. Tạ Quốc Dũng. Đề tài tập trung vào việc ứng dụng phần mềm chuyên dụng PROSPER để thiết lập mô hình, đánh giá hiệu suất khai thác và hiệu chuẩn động học dòng chảy giếng khí trong các điều kiện hoàn thiện và cấu trúc hình học khác nhau.

Nhiệm vụ và mục tiêu nghiên cứu của đồ án bao gồm:

  1. Xây dựng mô hình khai thác giếng khí khô và ướt (Dry and Wet Gas) đơn lẻ, đánh giá đường cong hiệu suất dòng vào (IPR), đường cong hiệu suất nâng thẳng đứng (VLP) và chạy phân tích độ nhạy.
  2. Thiết lập mô hình giếng khoan phân nhánh (Multilateral Well), xác định sản lượng và độ sụt giảm áp suất trên từng nhánh khai thác.
  3. Khảo sát ảnh hưởng của hệ thống đường ống bề mặt (Surface Pipeline) gắn kết với đầu giếng đến áp suất điểm nút và lưu lượng dòng.
  4. Lập mô hình giếng có cấu trúc hoàn thiện phức tạp nhằm kiểm soát cát, bao gồm hoàn thiện nứt vỉa kết hợp chèn sỏi (Frac and Packed well) và hoàn thiện bằng ống lọc chèn sẵn sỏi (Pre-packed screen).
  5. Thực hiện quy trình kiểm tra chất lượng (QC), hiệu chuẩn và khớp mô hình giếng khí (Well Test Matching) với số liệu đo thử giếng thực tế.
  6. Đánh giá cơ sở lựa chọn phương pháp khai thác nhân tạo (Artificial Lift) thích hợp cho giếng khai thác.

Đối tượng nghiên cứu là các mô hình giếng khí với các điều kiện biên thủy lực, đặc tính chất lưu PVT, cấu trúc ống tubing và phương pháp hoàn thiện giếng khác nhau. Phạm vi nghiên cứu giới hạn trong môi trường tính toán mô phỏng số của phần mềm PROSPER kết hợp với các bộ thông số kỹ thuật tiêu chuẩn trong công nghệ khai thác mỏ.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và mô hình kỹ thuật

Đồ án ứng dụng phương pháp phân tích điểm nút (Nodal Analysis) – công cụ nền tảng trong kỹ thuật khai thác dầu khí nhằm xác định điểm làm việc tối ưu giữa nguồn cung cấp chất lưu từ vỉa và khả năng vận chuyển của hệ thống ống ngầm/bề mặt:

  • Mối quan hệ hiệu suất dòng vào (Inflow Performance Relationship - IPR): Biểu diễn mối quan hệ giữa áp suất đáy giếng khi có dòng chảy ($P_{wf}$) và lưu lượng khai thác ($Q$). Các mô hình IPR phản ánh khả năng cấp dòng của vỉa, có xét đến hệ số tổn thất áp suất (Skin factor) do nhiễm bẩn thành hệ, nứt nẻ hoặc cấu trúc kiểm soát cát.
  • Hiệu suất nâng thẳng đứng (Vertical Lift Performance - VLP): Biểu diễn gradient áp suất dọc theo cột ống khai thác (Tubing). Mô hình VLP tính toán tổn thất áp suất do ma sát và thủy tĩnh của dòng chảy đa pha hoặc khí thực, sử dụng các tương quan đa pha chuyên dụng (chẳng hạn như tương quan Petroleum Experts 2).
  • Vận tốc giới hạn: Xác định hai ngưỡng thủy lực quan trọng:
    • Vận tốc tối thiểu Turner ($T$ - Minimum Turner Velocity): Vận tốc dòng khí nhỏ nhất cần thiết để nâng các hạt lỏng (nước hoặc condensate) lên bề mặt, tránh hiện tượng ngập lỏng đáy giếng.
    • Vận tốc xói mòn tối đa ($E$ - Erosional Velocity): Vận tốc dòng chảy cực đại cho phép trong đường ống nhằm tránh gây mài mòn, hư hỏng thiết bị hoặc sụt lở thành hệ.
  • Lưu lượng tiềm năng tuyệt đối (Absolute Open Flow - AOF): Lưu lượng khai thác lý thuyết lớn nhất khi áp suất đáy giếng bằng 0 (không có tổn thất áp suất do thành giếng hay thiết bị).

Phương pháp nghiên cứu và nguồn dữ liệu

Phương pháp nghiên cứu chính là mô phỏng số kỹ thuật khai thác trên phần mềm PROSPER (Petroleum Experts). Quy trình nghiên cứu bao gồm:

  • Thiết lập thuộc tính nhiệt động học chất lưu (PVT): Khai báo các mô hình chất lưu (Dry and Wet Gas), tỷ lệ khí - nước (WGR), tỷ lệ khí - condensate (CGR) và các thông số trạng thái.
  • Xây dựng cấu trúc hình học và vật lý của giếng: Khai báo khảo sát độ lệch giếng (Deviation Survey), thiết bị ngầm (Downhole Equipment), gradient địa nhiệt (Geothermal Gradient) và nhiệt dung trung bình (Average Heat Capacities).
  • Mô phỏng phân tích hệ thống và độ nhạy (Sensitivity Analysis): Đánh giá sự biến thiên của lưu lượng dòng khi thay đổi áp suất đầu giếng, tỷ lệ ngập nước (Water Cut), tỷ lệ khí dầu/khí ngưng tụ (GOR/CGR) và áp suất vỉa.
  • Khớp dữ liệu thử giếng (Well Test Matching): Sử dụng dữ liệu đo đạc thực tế tại giếng (áp suất đầu giếng, nhiệt độ đầu ống, áp suất đo tại độ sâu nhất định, lưu lượng khí, áp suất tĩnh vỉa) để hiệu chỉnh các hệ số tương quan VLP và thông số IPR.

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

Chương 1: Mô phỏng trên phần mềm PROSPER

Chương 1 là phần thực nghiệm kỹ thuật trọng tâm của đồ án, được chia thành các bài thực hành chuyên đề (Hướng dẫn/Tutorial) nhằm giải quyết từng bài toán mô phỏng cụ thể:

1. Hướng dẫn 1: Mô hình sản xuất khí khô và ướt

Nhóm tác giả thiết lập mô hình giếng khí đơn lẻ với tùy chọn hệ thống Dry and Wet Gas. Quá trình mô phỏng bao gồm:

  • Thiết lập đầy đủ dữ liệu PVT, quỹ đạo giếng, thiết bị ngầm, gradient địa nhiệt và nhiệt dung.
  • Xây dựng đường cong IPR có xét đến hệ số Skin.
  • Tính toán lưu lượng dòng chảy tại áp suất đầu giếng quy định 500 psi. Kết quả mô phỏng ban đầu cho thấy lưu lượng khí đạt 19.921 MMscf/d với áp suất điểm nút (Solution Node Pressure) là 2080 psi.
  • Đánh giá đồ thị VLP - IPR xác định điểm làm việc tối ưu là 19.434 Mscf/d. Tại mức lưu lượng khí vượt quá giá trị 33, vận tốc dòng chảy bắt đầu lớn hơn vận tốc xói mòn ($E$), cảnh báo nguy cơ sạt lở thành hệ và đòi hỏi phải hiệu chỉnh lại đường kính ống hoặc thông số VLP.
  • Xây dựng đường cong hiệu suất giếng thông qua đánh giá độ nhạy giữa áp suất đầu giếng và áp suất đáy giếng khi có dòng chảy ($P_{wf}$). Đường cong nâng (Lift curve - VLP) có dạng chữ J đặc trưng, thể hiện rõ vùng giới hạn vận tốc tối thiểu Turner ($T$) và vùng vận tốc xói mòn ($E$).

2. Hướng dẫn 5: Mô hình sản xuất khí khô phân nhánh

Bài toán mở rộng mô hình cho giếng khoan đa nhánh (Multilateral Well) với chất lưu Dry and Wet Gas:

  • Cấu trúc giếng được phân chia qua các điểm nút hình học: Điểm bắt đầu (Tie point), điểm phân nhánh (Joint) và các nhánh riêng biệt (Lateral x, Lateral y).
  • Kết quả IPR cho thấy lưu lượng tiềm năng đạt mức AOF cao. Tuy nhiên, tại giá trị áp suất Tie point hạ xuống 146.5 psi, đường cong IPR xuất hiện hiện tượng quay ngược (backward slope) do vận tốc dòng khí tăng quá nhanh làm gia tăng đột biến tổn thất ma sát trong giếng.
  • Giao điểm VLP - IPR xác định tổng lưu lượng khí khai thác của toàn hệ thống là 56.647 MMscf/d tại áp suất đầu giếng 371 psi.
  • Phân tích đóng góp chi tiết của từng nhánh:
    • Nhánh 1: Đạt lưu lượng 12.863 Mscf/day tại độ sâu đo 1940 ft (tính từ điểm phân nhánh) và áp suất nhánh 379 psi.
    • Nhánh 2: Đạt lưu lượng 43.286 Mscf/day tại độ sâu đo 2483 ft (tính từ điểm phân nhánh) và áp suất nhánh 467 psi.

3. Hướng dẫn 8: Mô hình giếng khí với đường ống

Mô hình mở rộng từ Hướng dẫn 1 bằng cách bổ sung đoạn đường ống bề mặt (Surface Pipeline) dài 5000 ft với đường kính trong (ID) là 6 inch nối từ đầu giếng đến gom góp (Manifold). Khi thiết lập áp suất tại manifold là 350 psi, kết quả tính toán điểm nút cho thấy lưu lượng khai thác tối ưu đạt 19.434 Mscf/day tại mức áp suất manifold tương ứng là 2093 psi.

4. Hướng dẫn 13: Lập mô hình frac và packed well

Nghiên cứu mô phỏng giếng hoàn thiện nứt vỉa thủy lực kết hợp chèn sỏi (Gravel pack) để kiểm soát cát:

  • Kích hoạt tùy chọn nứt vỉa thủy lực trong mô hình vỉa (Hydraulic Fractured Well) và tùy chọn kiểm soát cát (Sand Control) dạng Gravel Pack trong phần IPR.
  • Thiết lập tương quan nâng thẳng đứng Petroleum Experts 2, chọn nút giải pháp tại đáy giếng (Bottom node).
  • Đồ thị hệ thống thể hiện rõ độ sụt áp bổ sung gây ra bởi lớp sỏi chèn ($\Delta P$ do kiểm soát cát) và hệ số Skin tương ứng.

5. Hướng dẫn 14: Lập mô hình giếng với Pre-packed screen

Khảo sát hoàn thiện giếng dạng lỗ mở (Open Hole) hoặc giếng có chống ống (Cased Hole) sử dụng lưới quấn dây chèn sẵn sỏi (Pre-packed screen):

  • Cấu hình kiểm soát cát được cập nhật trực tiếp vào mô hình IPR để xác định độ suy giảm dòng chảy vào.
  • Tính toán giao điểm IPR - VLP với tương quan Petroleum Experts 2 giúp phân tích định lượng giá trị tổng Skin, tổn thất áp suất $\Delta P$ qua màn chắn cát và xác định chính xác lưu lượng làm việc của giếng.

6. Hướng dẫn 19: Thử nghiệm giếng khí khô phù hợp

Quy trình hiệu chuẩn mô hình giếng khí PROSPER dựa trên bộ dữ liệu thử giếng (Well Test Data) được chuẩn hóa theo 4 bước tuần tự:

  1. Kiểm tra dữ liệu thử giếng: Đánh giá tính nhất quán và độ tin cậy của các thông số đo đạc.
  2. Kiểm tra chất lượng PVT: Đánh giá đặc tính chất lưu (với giếng khí khô và condensate hồi lưu, bước này thường không yêu cầu so match phức tạp).
  3. Kiểm tra và hiệu chỉnh tương quan VLP: Lựa chọn tương quan dòng đa pha phù hợp và khớp với phản ứng thực tế của ống khai thác.
  4. Hiệu chỉnh mô hình IPR: Điều chỉnh các thông số vỉa và hệ số Skin để tái tạo chính xác kết quả thử vỉa.
Thông số thử giếng (Well Test) Giá trị thực tế ghi nhận
Ngày thử nghiệm 01/01/2010
Chế độ kiểm tra Kiểm tra luồng chính (Main flow test)
Áp suất đầu giếng ($P_{wh}$) 1000 psig
Nhiệt độ đầu ống chảy ($T_{wh}$) 153 °F
Tỷ lệ khí - nước (WGR) 5 stb/MMscf
Tỷ lệ khí - condensate (CGR) 5 stb/MMscf
Tốc độ dòng khí ($Q_g$) 15 MMscf/d
Độ sâu đặt máy đo 4500 ft
Áp suất đo được tại máy 1920 psig
Áp suất tĩnh vỉa (đỉnh lỗ bắn) 2300 psi

Mô hình sau khi được khớp hoàn chỉnh được sử dụng để ước tính hệ số truyền nhiệt tổng thể ($U$) và chạy các kịch bản độ nhạy theo sự suy giảm áp suất vỉa.

Chương 2: Lựa chọn Artificial Lift phù hợp

Chương 2 của đồ án (trang 60–66 trong tài liệu gốc) trình bày cơ sở lý thuyết và phương pháp luận kỹ thuật để lựa chọn hệ thống khai thác nhân tạo tối ưu khi các giếng tự phun bị suy giảm áp suất hoặc ngập lỏng:

  • Đánh giá khả năng ứng dụng của các phương pháp khai thác nhân tạo phổ biến trong công nghiệp dầu khí (như Gas Lift, bơm điện chìm ESP, bơm trục vít PCP, bơm pít-tông cần hút Sucker Rod Pump).
  • Phân tích các tiêu chí kỹ thuật: độ sâu giếng, lưu lượng kỳ vọng, tỷ lệ khí-dầu (GOR), độ nhớt chất lưu, sự hiện diện của cát/chất rắn và áp suất đầu giếng.
  • Tích hợp đường cong VLP sau khi lắp đặt thiết bị Artificial Lift vào phần mềm PROSPER để dự báo khả năng nâng chất lưu và khôi phục lưu lượng kinh tế cho giếng.

Kết quả và đóng góp

Kết quả kỹ thuật chính

Đồ án đã lượng hóa chi tiết các thông số khai thác thông qua mô phỏng trên PROSPER:

Kịch bản mô hình Thông số đầu vào chính Kết quả lưu lượng Áp suất làm việc / Điểm nút
Giếng khí đơn lẻ (HD 1) $P_{wh} = 500\text{ psi}$, Dry & Wet Gas $19.921\text{ MMscf/d}$ (tính toán)
$19.434\text{ Mscf/d}$ (tối ưu)		 $P_{node} = 2080\text{ psi}$ (Vượt ngưỡng xói mòn tại $Q > 33$)
Giếng đa nhánh (HD 5) $P_{wh} = 200\text{ psi}$, Multilateral $56.647\text{ MMscf/d}$ (tổng giếng)
Nhánh 1: $12.863\text{ Mscf/d}$
Nhánh 2: $43.286\text{ Mscf/d}$		 $P_{wh} = 371\text{ psi}$
Nhánh 1: $379\text{ psi}$ (1940 ft)
Nhánh 2: $467\text{ psi}$ (2483 ft)
Giếng gắn đường ống (HD 8) Tuyến ống $5000\text{ ft}$, $\text{ID} = 6\text{ inch}$, $P_{manifold} = 350\text{ psi}$ $19.434\text{ Mscf/d}$ $P_{manifold} = 2093\text{ psi}$
Khớp thử giếng (HD 19) $P_{stat} = 2300\text{ psi}$, $P_{wh} = 1000\text{ psig}$, độ sâu $4500\text{ ft}$ $15.000\text{ MMscf/d}$ $P_{gauge} = 1920\text{ psig}$ ($T_{wh} = 153\text{ }^\circ\text{F}$)

Đóng góp của đề tài

  1. Chuẩn hóa quy trình mô phỏng Nodal Analysis: Hệ thống hóa các bước thiết lập từ thuộc tính PVT, khảo sát quỹ đạo giếng, gradient nhiệt độ đến tương quan nâng thủy lực cho các cấu trúc giếng khí đa dạng.
  2. Làm rõ các giới hạn thủy lực giếng khí: Xác định trực quan mối liên hệ giữa vận tốc xói mòn ($E$) và vận tốc nâng tối thiểu Turner ($T$), cung cấp cơ sở kỹ thuật để điều chỉnh kích thước ống tubing và kiểm soát áp suất đầu giếng nhằm bảo vệ thành giếng.
  3. Phân tích hiện tượng thủy lực phi tuyến trong giếng đa nhánh: Phát hiện và giải thích cơ chế đường cong IPR bị quay ngược khi áp suất tại điểm kết nối (Tie point) giảm xuống mức quá thấp (146.5 psi).
  4. Quy trình khớp số liệu thử giếng có hệ thống: Xây dựng quy trình 4 bước giúp hiệu chỉnh mô hình số khớp với dữ liệu đo thực địa, nâng cao độ chính xác khi dự báo suy giảm áp suất vỉa.

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

  • Hạn chế:
    • Đề tài được thực hiện trong thời gian ngắn của bài tập cuối khóa; các tác giả tự nhận định kiến thức chuyên môn còn hạn chế nên chưa phân tích sâu toàn diện tất cả các trường hợp biến động vỉa phức tạp.
    • Các bài toán chủ yếu dựa trên các bộ dữ liệu mẫu định hình sẵn trong các module hướng dẫn của PROSPER, chưa phản ánh đầy đủ các dị thường địa chất hoặc cơ chế ngưng tụ phức tạp của các mỏ khí - condensate thực tế tại thềm lục địa Việt Nam.
  • Hướng phát triển tiếp theo:
    • Mở rộng mô phỏng tích hợp toàn mỏ thông qua việc kết nối mô hình giếng PROSPER với mô hình cân bằng vật chất (MBAL) và mô hình mạng lưới thu gom bề mặt (GAP).
    • Phân tích chi tiết tính toán kinh tế - kỹ thuật cho từng hệ thống Artificial Lift cụ thể khi chuyển đổi từ giai đoạn tự phun sang khai thác cơ học trong suốt vòng đời mỏ.

Giá trị tham khảo

Báo cáo bài tập cuối khóa là tài liệu tham khảo thiết thực cho:

  • Sinh viên ngành Kỹ thuật Địa chất, Kỹ thuật Dầu khí và Công nghệ Khai thác: Sử dụng làm tài liệu hướng dẫn từng bước (step-by-step tutorial) để làm chủ phần mềm mô phỏng PROSPER.
  • Kỹ sư khai thác và công nghệ mỏ mới vào nghề: Tham khảo quy trình phân tích điểm nút (Nodal Analysis), phương pháp lập mô hình kiểm soát cát (Gravel Pack, Pre-packed screen), và các bước thực hiện QC/Matching dữ liệu thử giếng.
  • Các đề tài môn học và đồ án tốt nghiệp: Cung cấp mẫu đối sánh số liệu, biểu đồ IPR/VLP và phương pháp luận chạy độ nhạy hệ thống khai thác.

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

1. Ý nghĩa của hai đại lượng $T$ và $E$ trên đường cong VLP trong Hướng dẫn 1 là gì?

$T$ đại diện cho vận tốc tối thiểu Turner (minimum Turner velocity) – tức vận tốc dòng khí nhỏ nhất để cuốn các hạt chất lỏng lên bề mặt, tránh hiện tượng ngập lỏng đáy giếng. $E$ là vận tốc xói mòn tối đa (erosional velocity) – nếu lưu lượng vượt qua ngưỡng này, chất lưu sẽ gây mài mòn, phá hủy đường ống và sạt lở thành hệ giếng.

2. Tại sao đường cong IPR của giếng phân nhánh trong Hướng dẫn 5 lại có hiện tượng quay ngược tại áp suất Tie point = 146.5 psi?

Hiện tượng quay ngược xảy ra do khi áp suất tại Tie point hạ xuống quá thấp (146.5 psi), thể tích riêng của khí nở ra làm vận tốc khí tăng vọt. Vận tốc quá lớn dẫn tới tổn thất áp suất do ma sát tăng đột biến dọc thân giếng, khiến cho lưu lượng thực tế đi vào giếng bị suy giảm thay vì tăng lên.

3. Trong Hướng dẫn 13 và 14, việc hoàn thiện kiểm soát cát ảnh hưởng đến những đại lượng nào trong PROSPER?

Việc hoàn thiện giếng bằng nứt vỉa chèn sỏi (Frac-pack) hoặc ống lọc chèn sẵn sỏi (Pre-packed screen) làm thay đổi đường cong IPR, bổ sung thêm hệ số tổn thất áp suất (Skin kiểm soát cát) và tạo ra độ giảm áp suất phụ trợ ($\Delta P$) qua lớp sỏi chèn tại vùng đáy giếng.

4. Quy trình 4 bước khớp dữ liệu thử giếng (Well Test Matching) trong Hướng dẫn 19 gồm những gì?

Quy trình gồm 4 bước tuần tự:

  1. Xem xét và đánh giá độ tin cậy của số liệu kiểm tra giếng.
  2. Kiểm tra chất lượng và hiệu chuẩn mô hình PVT.
  3. Kiểm tra chất lượng, lựa chọn, hiệu chỉnh và xác nhận tương quan nâng thẳng đứng VLP.
  4. Điều chỉnh các thông số mô hình IPR để tái tạo kết quả thử giếng thực tế.

5. Khi thực hiện phân tích độ nhạy trên màn hình tính toán hệ thống của PROSPER, các giá trị Water Cut và GOR được xử lý như thế nào?

Khi người dùng kích hoạt biến độ nhạy cho Water Cut hoặc GOR trên màn hình tính toán hệ thống, PROSPER sẽ sử dụng trực tiếp các giá trị độ nhạy mới này cho toàn bộ quá trình tính toán và bỏ qua các giá trị cố định đã nhập trước đó trên màn hình IPR hoặc PVT.

Kết luận

Báo cáo bài tập cuối khóa môn Công nghệ Khai thác Dầu khí của nhóm sinh viên Trường Đại học Bách Khoa – ĐHQG-HCM đã xây dựng thành công quy trình mô phỏng toàn diện các mô hình giếng khí trên phần mềm PROSPER. Nghiên cứu đã giải quyết trọn vẹn từ các bài toán giếng khí đơn lẻ, giếng phân nhánh, hệ thống giếng nối đường ống bề mặt đến các cấu trúc hoàn thiện kiểm soát cát phức tạp và hiệu chuẩn mô hình theo số liệu thử giếng thực tế. Kết quả đạt được cung cấp cơ sở dữ liệu kỹ thuật tin cậy phục vụ phân tích thủy lực giếng khoan và định hướng lựa chọn giải pháp khai thác nhân tạo tối ưu trong công nghệ khai thác mỏ.