Tổng quan nghiên cứu

Trong kỹ thuật vỉa và khai thác dầu khí hiện đại, việc xác định chính xác các đặc tính vật lý cũng như hành vi thể tích của lưu chất vỉa theo sự thay đổi áp suất và nhiệt độ là yếu tố then chốt quyết định hiệu quả kinh tế của toàn bộ dự án mỏ. Trong thực tế khai thác ở các vỉa sâu, áp suất chất lưu thường dao động từ 1.500 psia đến hơn 3.500 psia tại nhiệt độ từ 150°F đến 180°F. Tại điều kiện khắc nghiệt này, việc áp dụng phương trình trạng thái khí lý tưởng thông thường dẫn đến sai số nghiêm trọng lên tới 500%, trái ngược hoàn toàn với mức sai số nhỏ chỉ khoảng 2% đến 3% tại điều kiện bề mặt tiêu chuẩn (14,7 psia và 60°F, tương đương 520°R).

Vấn đề nghiên cứu trọng tâm đặt ra là làm thế nào để dự đoán chính xác các thuộc tính nhiệt động lực học và hành vi thể tích (PVT) của hỗn hợp khí tự nhiên khi thiếu hụt dữ liệu đo thực nghiệm trực tiếp từ phòng thí nghiệm. Mục tiêu cụ thể của luận văn là hệ thống hóa cơ sở lý thuyết về khí lý tưởng và khí thực, đánh giá độ tin cậy của hệ số nén khí Z (gas deviation factor), phân tích biểu đồ chuẩn Standing - Katz và kiểm chứng các phương pháp hiệu chỉnh phi hydrocacbon (Wichert - Aziz, Carr - Kobayashi - Burrows) cũng như hiệu chỉnh phân đoạn hydrocacbon nặng C7+ theo mô hình Sutton.

Nghiên cứu được triển khai trong phạm vi chuyên ngành Kỹ thuật Dầu khí tại Trường Đại học Bách khoa - ĐHQG TP.HCM, tập trung vào các hệ khí tự nhiên ngọt, khí chua và khí ngưng tụ (gas-condensate). Kết quả của công trình mang lại giá trị thực tiễn to lớn: giúp kiểm soát sai số tính toán thể tích và khối lượng riêng khí ở mức dưới 1,5%, giảm thiểu 15% chi phí rủi ro khi lập kế hoạch khai thác và nâng cao độ chuẩn xác trong việc tính toán trữ lượng khí ban đầu tại chỗ (GIIP).

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

Khung lý thuyết áp dụng

Khung lý thuyết của công trình được xây dựng trên nền tảng nhiệt động lực học chất lưu dầu khí với hai trụ cột chính:

Thứ nhất, thuyết động học chất khí và phương trình trạng thái khí lý tưởng pV = nRT, trong đó n biểu diễn số mol khí thông qua khối lượng m và trọng lượng phân tử M, với hằng số khí phổ quát R = 10,730 psia.ft³/lb-mole.°R. Từ nền tảng này, các khái niệm cơ sở được chuẩn hóa bao gồm: trọng lượng phân tử biểu kiến (Ma), thể tích tiêu chuẩn (Vsc = 379,4 scf/lb-mol tại 14,7 psia và 60°F), tỷ trọng tương đối của khí (gas specific gravity γg so với không khí có Mair = 28,96 lb/lb-mol), thể tích riêng và khối lượng riêng chất khí.

Thứ hai, lý thuyết khí thực và hệ số nén Z. Để bù đắp sự sai lệch thể tích do lực tương tác giữa các phân tử khí ở áp suất cao, phương trình trạng thái được hiệu chỉnh thành pV = ZnRT. Nghiên cứu áp dụng định luật trạng thái tương ứng thông qua hai đại lượng không thứ nguyên: áp suất giả rút gọn (ppr = p/ppc) và nhiệt độ giả rút gọn (Tpr = T/Tpc), với áp suất giả tới hạn (ppc) và nhiệt độ giả tới hạn (Tpc) được xác định theo quy tắc cộng phần mol của Kay hoặc tính toán từ tỷ trọng khí theo tương quan thực nghiệm Brown và Standing. Bên cạnh đó, các mô hình hiệu chỉnh phi hydrocacbon Wichert - Aziz (sử dụng hệ số hiệu chỉnh nhiệt độ ε dựa trên hàm lượng H2S và CO2), Carr - Kobayashi - Burrows và mô hình trộn phân đoạn nặng Sutton cho thành phần C7+ cấu thành khung lý thuyết hoàn chỉnh cho hệ khí thực.

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

Nghiên cứu sử dụng tập dữ liệu thực nghiệm gồm 45 bộ mẫu chất lưu vỉa đại diện cho các mỏ khí tự nhiên và khí ngưng tụ. Phương pháp chọn mẫu là kỹ thuật lấy mẫu đáy giếng (bottom-hole sampling) kết hợp lấy mẫu tái tổ hợp tại bình tách bề mặt (surface recombination sampling) tuân thủ nghiêm ngặt tiêu chuẩn API RP 44, đảm bảo phản ánh đầy đủ dải biến thiên của tỷ trọng khí từ 0,65 đến 0,85 và dải tạp chất phi hydrocacbon chứa tới 10% H2S cùng 5% CO2.

Phương pháp phân tích kết hợp phân tích sắc ký khí trong phòng thí nghiệm để xác định thành phần phần mol tử số (yi) với phương pháp tính toán giải tích số hóa biểu đồ Standing - Katz. Lý do lựa chọn phương pháp phân tích giải tích thông qua các thuật toán Hall - Yarborough và Dranchuk - Abu - Kassem (với hệ số nén tới hạn zc xấp xỉ 0,27) là nhằm loại bỏ sai số đọc thủ công từ biểu đồ đồ họa, tăng tốc độ xử lý dữ liệu và cho phép tích hợp trực tiếp vào các phần mềm mô phỏng dòng chảy trong vỉa. Toàn bộ quá trình thu thập, xử lý thực nghiệm và mô hình hóa toán học được thực hiện liên tục trong thời gian 12 tháng.

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

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

Nghiên cứu đã làm sáng tỏ bốn phát hiện định lượng cốt lõi:

Thứ nhất, mức độ sai lệch thể tích của khí lý tưởng tăng theo hàm số phi tuyến tính với áp suất vỉa. Tại áp suất khí quyển 14,7 psia, sai số chỉ ở mức 2,1%, nhưng khi áp suất vỉa tăng lên ngưỡng 2.000 psia đến 3.000 psia tại 180°F, sai số thể tích nếu bỏ qua hệ số nén Z lên tới 480% đến 500%. Điều này khẳng định việc áp dụng phương trình khí lý tưởng trong tính toán trữ lượng vỉa là hoàn toàn không khả thi.

Thứ hai, tác động nghiêm trọng của các thành phần phi hydrocacbon lên hệ số nén Z. Đối với khí ngọt chứa dưới 5% CO2 và N2, biểu đồ Standing - Katz cho kết quả đáng tin cậy. Tuy nhiên, đối với hệ khí chua chứa 5% CO2 và 10% H2S ở áp suất 3.500 psia và nhiệt độ 160°F, sự phân cực phân tử làm thay đổi đáng kể hành vi nén, gây sai lệch hệ số Z tới 9,8% nếu không hiệu chỉnh. Khi áp dụng phương pháp Wichert - Aziz để điều chỉnh Tpc và ppc, sai số tính toán khối lượng riêng khí đã giảm từ 9,8% xuống chỉ còn 1,15%.

Thứ ba, hiệu quả vượt trội của quy tắc trộn Sutton đối với các hệ khí có tỷ trọng cao chứa phân đoạn nặng C7+ (với trọng lượng phân tử 161 lb/lb-mol và tỷ trọng 0,80). Việc áp dụng các hệ số hiệu chỉnh FJ, EJ, EK giúp cải thiện độ chính xác tính toán khối lượng riêng khí thêm 18,4% so với việc áp dụng trực tiếp quy tắc cộng thông thường của Kay ở điều kiện 2.000 psi và 150°F.

Thứ tư, độ chính xác của các phương pháp tính giải tích trực tiếp. Thuật toán Dranchuk - Abu - Kassem và Hall - Yarborough tái lập hoàn hảo các đường cong biểu diễn hệ số Z của Standing - Katz với sai số tương đối trung bình chỉ 0,38% trên toàn bộ dải ppr từ 0 đến 15 và Tpr từ 1,05 đến 3,0.

Thảo luận kết quả

Nguyên nhân căn bản của các sai lệch nêu trên bắt nguồn từ lực hút liên phân tử Van der Waals và thể tích chiếm chỗ của các phân tử hydrocacbon khi bị nén ở mật độ cao. Ở áp suất trung bình (ppr từ 1,5 đến 3,0), lực hút chiếm ưu thế làm khí nén nhiều hơn khí lý tưởng (Z < 1,0), trong khi ở áp suất rất cao (ppr > 5,0), lực đẩy chiếm ưu thế khiến khí khó nén hơn (Z > 1,0).

Để trực quan hóa các kết quả này, toàn bộ dữ liệu thực nghiệm và giải tích có thể được biểu diễn trực quan qua biểu đồ tương quan hai chiều: trục hoành thể hiện áp suất giả rút gọn ppr (từ 0 đến 15), trục tung thể hiện hệ số Z (từ 0,2 đến 1,2), với họ các đường cong tham số theo nhiệt độ giả rút gọn Tpr (từ 1,05 đến 3,0). Bảng so sánh đa chỉ tiêu giữa dữ liệu thực nghiệm đo trong tế bào PVT và các mô hình giải tích cho thấy hệ số tương quan tuyến tính R² đạt trên 0,998. Kết quả này chứng minh rằng việc kết hợp hiệu chỉnh Wichert - Aziz hoặc Sutton với thuật toán Dranchuk - Abu - Kassem mang lại độ chính xác tương đương thực nghiệm phòng thí nghiệm với thời gian tính toán giảm hơn 95%.

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

Dựa trên các kết quả đạt được, luận văn đưa ra bốn nhóm giải pháp hành động cụ thể nhằm nâng cao hiệu quả kỹ thuật và kinh tế trong khai thác mỏ:

Thứ nhất, chuẩn hóa quy trình tính toán PVT cho các mỏ khí chua bằng cách áp dụng bắt buộc thuật toán hiệu chỉnh Wichert - Aziz đối với tất cả các dòng khí có tổng nồng độ H2S và CO2 vượt quá 5%. Mục tiêu là giảm thiểu sai số tính toán trữ lượng khí xuống dưới 1,5%. Quy trình này cần được hoàn thiện và ban hành trong vòng 3 tháng bởi Đội ngũ Kỹ sư Công nghệ Mỏ và Khai thác.

Thứ hai, tích hợp thuật toán giải tích Dranchuk - Abu - Kassem và Hall - Yarborough vào hệ thống phần mềm mô phỏng vỉa nội bộ thay thế hoàn toàn cho việc tra cứu đồ thị thủ công. Giải pháp này giúp tăng tốc độ xử lý mô phỏng thủy động lực học lên 300% và loại trừ hoàn toàn sai sót chủ quan của người vận hành, thực hiện trong thời hạn 6 tháng bởi Phòng Công nghệ Thông tin Dầu khí.

Thứ ba, nâng cấp quy trình phân tích mẫu chất lưu tại phòng thí nghiệm thông qua việc kiểm chuẩn định kỳ thiết bị đo tế bào PVT với chu kỳ 6 tháng một lần, đảm bảo độ không đảm bảo đo tỷ trọng khí đạt dưới 0,5% tại điều kiện áp suất cao lên tới 5.000 psia, do Phòng Thí nghiệm Phân tích Lưu chất phụ trách.

Thứ tư, thiết lập hệ thống giám sát tự động nồng độ tạp chất phi hydrocacbon trực tiếp tại đầu giếng với ngưỡng cảnh báo tức thời khi nồng độ CO2 hoặc H2S biến động vượt quá 10% so với giá trị thiết kế ban đầu, thực hiện trong vòng 45 ngày bởi Đội ngũ Kỹ sư Vận hành Giàn khoan.

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

Luận văn là nguồn tài liệu học thuật và kỹ thuật chuyên sâu, đặc biệt giá trị đối với bốn nhóm đối tượng sau:

Thứ nhất, các kỹ sư công nghệ vỉa và kỹ sư khai thác dầu khí tại các công ty năng lượng. Luận văn cung cấp hướng dẫn từng bước về cách tính toán hệ số Z, khối lượng riêng và thể tích riêng của khí ở dải áp suất từ 1.500 đến 3.500 psia, phục vụ trực tiếp công tác lập mô hình mỏ và tối ưu hóa lưu lượng khai thác giếng.

Thứ hai, chuyên viên phân tích dữ liệu chất lưu tại các phòng thí nghiệm PVT. Tài liệu cung cấp cơ sở lý thuyết chuẩn xác để kiểm tra chéo độ tin cậy của các kết quả phân tích mẫu khí tự nhiên và khí ngưng tụ trước khi bàn giao báo cáo cho chủ đầu tư.

Thứ ba, giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Dầu khí và Địa chất Công trình. Đây là tài liệu tham khảo chuẩn mực với hệ thống công thức toán học khép kín từ lý thuyết khí lý tưởng đến các mô hình khí thực phức tạp kèm ví dụ tính toán số liệu cụ thể.

Thứ tư, các nhà quản lý dự án phát triển mỏ và chuyên gia thẩm định trữ lượng. Luận văn giúp đánh giá rủi ro tài chính phát sinh từ sự biến động 5% đến 10% trong tính toán trữ lượng khí do sai lệch thông số PVT.

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

Tại sao không thể sử dụng phương trình khí lý tưởng để tính toán thể tích khí vỉa?

Tại áp suất khí quyển 14,7 psia, phương trình khí lý tưởng chỉ có sai số khoảng 2%. Tuy nhiên, trong điều kiện vỉa sâu với áp suất từ 1.500 psia đến 3.000 psia, lực hút và đẩy giữa các phân tử khí hoạt động rất mạnh, khiến sai số thể tích lên tới 500%, làm sai lệch hoàn toàn tính toán trữ lượng mỏ.

Biểu đồ Standing - Katz có thể áp dụng chính xác cho những loại khí nào?

Biểu đồ Standing - Katz áp dụng chính xác nhất cho các hệ khí ngọt tự nhiên chứa chủ yếu hydrocacbon nhẹ với áp suất giả rút gọn từ 0 đến 15 và nhiệt độ giả rút gọn từ 1,05 đến 3,0. Khi khí chứa trên 5% tạp chất như CO2 hoặc H2S, biểu đồ sẽ xuất hiện sai số tới 10% nếu không được hiệu chỉnh.

Khi nào cần áp dụng phương pháp hiệu chỉnh Wichert - Aziz?

Phương pháp Wichert - Aziz bắt buộc phải áp dụng khi hỗn hợp khí tự nhiên có chứa khí chua với tổng phần mol của CO2 và H2S vượt quá 5%. Phương pháp này sử dụng hệ số điều chỉnh nhiệt độ ε để tính lại nhiệt độ và áp suất giả tới hạn, giúp đưa sai số tính Z-factor về dưới 1,2%.

Phân đoạn nặng C7+ ảnh hưởng như thế nào đến việc xác định tính chất khí?

Phân đoạn heptane-plus (C7+) có trọng lượng phân tử lớn (thường trên 160 lb/lb-mol) và cấu trúc phân tử phức tạp, làm sai lệch đáng kể nhiệt độ và áp suất giả tới hạn của toàn bộ hệ khí. Sử dụng mô hình Sutton với các hệ số FJ, EJ, EK giúp kiểm soát sai số tính tỷ trọng khí dưới 1,1%.

Giữa hai phương pháp Hall - Yarborough và Dranchuk - Abu - Kassem, phương pháp nào vượt trội hơn?

Cả hai phương pháp đều khớp với biểu đồ Standing - Katz với sai số cực nhỏ dưới 0,5%. Tuy nhiên, phương pháp Dranchuk - Abu - Kassem với hệ số nén tới hạn zc xấp xỉ 0,27 có độ ổn định thuật toán cao hơn trên toàn dải áp suất và được tích hợp rộng rãi hơn trong các phần mềm mô phỏng thương mại.

Kết luận

  • Công trình đã làm sáng tỏ giới hạn sai số nghiêm trọng lên tới 500% của phương trình khí lý tưởng tại điều kiện áp suất vỉa cao 3.000 psia, khẳng định tính bắt buộc của hệ số nén Z trong kỹ thuật vỉa.
  • Chuẩn hóa thành công quy trình hiệu chỉnh tính chất giả tới hạn cho các hệ khí chua chứa tạp chất CO2 và H2S bằng phương pháp Wichert - Aziz, đưa sai số tính toán về dưới 1,2%.
  • Chứng minh tính ưu việt của quy tắc trộn Sutton trong việc nâng cao độ chính xác xác định khối lượng riêng khí cho các dòng khí ngưng tụ chứa phân đoạn nặng C7+ thêm 18,4%.
  • Xác nhận độ tin cậy tuyệt đối của giải thuật Dranchuk - Abu - Kassem với độ sai lệch chỉ 0,38% so với biểu đồ chuẩn Standing - Katz phục vụ tự động hóa tính toán.
  • Cung cấp khung lý thuyết hoàn chỉnh và các ví dụ thực tế phục vụ trực tiếp công tác đào tạo học thuật và ứng dụng sản xuất trong ngành dầu khí.

Kế hoạch chuyển giao công nghệ và tích hợp các module giải thuật vào hệ sinh thái phần mềm quản lý khai thác mỏ dự kiến sẽ hoàn tất trong 6 tháng tới. Các kỹ sư và chuyên gia công nghệ mỏ hãy áp dụng ngay các công thức và quy trình hiệu chỉnh giải tích này vào công tác tính toán PVT thực địa để tối ưu hóa hiệu quả kinh tế và độ an toàn cho các dự án khai thác khí.