Tổng quan nghiên cứu

Trong kỹ thuật khai thác dầu khí, việc mô phỏng chính xác tổn thất áp suất dọc theo ống nâng đóng vai trò quyết định đến hiệu quả khai thác và tuổi thọ giếng khoan. Thống kê thực tế cho thấy thành phần tổn thất áp suất do cột áp thủy tĩnh chiếm tỷ trọng áp đảo từ 80% đến 95%, trong khi tổn thất do ma sát thành ống chiếm từ 5% đến 20% tổng tổn hao áp suất trong giếng. Tuy nhiên, khi xuất hiện dòng chảy nhiều pha giữa dầu, khí và nước, các tương quan thủy động lực học trở nên vô cùng phức tạp do hiện tượng trượt pha và sự thay đổi liên tục của các chế độ dòng chảy.

Tại mỏ Báo Đen thuộc bồn trũng Cửu Long, việc áp dụng các phần mềm thương mại tiêu chuẩn như PIPESIM để mô phỏng giếng khai thác của Công ty Liên doanh Điều hành Cửu Long thường ghi nhận sai số đáng kể, dao động từ 8% đến 15% so với áp suất đo thực tế. Nguyên nhân chính là các mô hình tương quan thực nghiệm và cơ học tích hợp sẵn được xây dựng dựa trên điều kiện thí nghiệm chuẩn, chưa phản ánh chính xác tính chất cơ lý của dầu thô nhiệt đới giàu paraffin tại thềm lục địa Việt Nam.

Nghiên cứu được thực hiện vào tháng 7 năm 2012 nhằm giải quyết bài toán cấp thiết này với mục tiêu cốt lõi: xây dựng và hiệu chỉnh mô hình dòng chảy cơ học nhiều pha cho ống nâng giếng khai thác mỏ Báo Đen, kết hợp tinh chỉnh các thông số nhiệt động học chất lưu dựa trên kết quả phân tích mẫu thực nghiệm và dữ liệu đo kiểm tra khai thác thực tế. Nghiên cứu tập trung khảo sát chuyên sâu tại giếng BD-7P và áp dụng kiểm chứng độc lập trên giếng BD-1P. Kết quả của luận văn đã giúp giảm sai số dự báo áp suất đáy giếng xuống dưới mức 3.2%, cung cấp công cụ tính toán thủy lực có độ tin cậy cao nhằm tối ưu hóa đường kính ống khai thác và nâng cao sản lượng khai thác cho toàn mỏ.

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 xây dựng trên nền tảng phương trình bảo toàn năng lượng và phương trình động lượng tổng quát cho dòng chảy chất lưu một chiều trong ống nghiêng. Tổng gradient áp suất dọc chiều dài ống được phân tích thành ba thành phần độc lập: gradient do tổn thất ma sát Fanning, gradient do thế năng cột chất lỏng và gradient do gia tốc biến thiên động năng. Trong đó, hệ số ma sát được xác định thông qua chuẩn số Reynolds nhằm phân biệt chính xác trạng thái chảy tầng, chảy chuyển tiếp và chảy rối.

Đối với dòng chảy nhiều pha, luận văn ứng dụng mô hình cơ học cải tiến của Ansari kết hợp khảo sát so sánh với các mô hình tương quan thực nghiệm kinh điển như Hagedorn và Brown, Beggs và Brill, Duns và Ros, Orkiszewski. Mô hình cơ học của Ansari chia dòng chảy thẳng đứng thành 5 chế độ thủy động lực học đặc thù: dòng chảy dạng bọt khí, dạng bọt khí phân tán, dạng nút khí, dạng khuấy và dạng vành xuyến. Khung lý thuyết này vận dụng các khái niệm then chốt bao gồm: tỷ lệ giữ chất lỏng thực tế, vận tốc trượt pha giữa khí và lỏng, động học bọt khí Taylor, tỷ số khí hòa tan trong dầu, hệ số thể tích thành hệ dầu và hệ số lệch khí. Việc mô tả chính xác sự chuyển tiếp giữa các chế độ dòng chảy giúp loại bỏ các giả định hỗn hợp đồng nhất thiếu chính xác trong các phương pháp thực nghiệm trước đây.

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

Nguồn dữ liệu nghiên cứu bao gồm bộ số liệu phân tích nhiệt động học chất lưu PVT từ phòng thí nghiệm và chuỗi dữ liệu áp suất, nhiệt độ theo chiều sâu đo bằng thiết bị đo kiểm tra khai thác PLT tại giếng BD-7P và giếng BD-1P thuộc mỏ Báo Đen. Phương pháp chọn mẫu tập trung vào các giếng khai thác đại diện có điều kiện làm việc ổn định, với chiều sâu ống nâng khảo sát hơn 3000 mét và tập hợp hàng chục điểm đo áp suất thực tế dọc thân giếng.

Quy trình phân tích và xử lý số liệu được tiến hành qua hai giai đoạn chặt chẽ:

  • Giai đoạn 1: Hiệu chỉnh các tương quan nhiệt động học PVT bao gồm tỷ số khí hòa tan, hệ số lệch khí, hệ số thể tích thành hệ và độ nhớt chất lưu bằng phương pháp thiết kế thực nghiệm đa yếu tố 2 mức và 3 mức kết hợp phân tích phương sai ANOVA trên phần mềm Minitab.
  • Giai đoạn 2: Ứng dụng phương pháp bình phương tối thiểu và thuật toán tối ưu hóa trong phần mềm PIPESIM cùng Excel để hiệu chỉnh các hệ số ma sát của nút chất lỏng, hệ số tổn thất áp suất và tỷ lệ chiều dài bọt khí Taylor trong mô hình cơ học. Lý do lựa chọn phương pháp hiệu chỉnh hai bước tách biệt này là nhằm cô lập sai số nhiệt động học của chất lưu trước khi tối ưu hóa phương trình thủy lực, đảm bảo mô hình đạt độ chính xác vật lý cao nhất.

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

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

Quá trình phân tích dữ liệu thực nghiệm và tính toán mô phỏng thủy lực tại mỏ Báo Đen đã mang lại 4 phát hiện quan trọng:

  • Thứ nhất, phân tích cơ cấu tổn thất áp suất khẳng định thành phần cột áp thủy tĩnh chiếm tỷ trọng từ 82% đến 94% tổng tổn thất trong ống nâng giếng BD-7P, trong khi tổn thất do ma sát thành ống chỉ chiếm từ 6% đến 18% và tổn thất do biến thiên động năng chiếm dưới 1%.
  • Thứ hai, các công thức thực nghiệm PVT mặc định của Standing và Beggs - Robinson tạo ra độ lệch từ 15% đến 25% đối với tỷ số khí hòa tan và độ nhớt dầu thô mỏ Báo Đen. Sau khi thực hiện thiết kế thực nghiệm đa yếu tố và phân tích hồi quy, mô hình PVT mới đã giảm sai số tương đối xuống dưới 2.5% so với số liệu đo trong phòng thí nghiệm.
  • Thứ ba, mô hình cơ học Ansari nguyên bản khi đưa vào phần mềm PIPESIM cho thấy sai số tính toán phân bố áp suất dọc thân giếng BD-7P dao động từ 8.5% đến 12.3% so với dữ liệu đo PLT do đánh giá chưa chuẩn xác hình thái dòng chảy dạng nút.
  • Thứ tư, sau khi hoàn tất hiệu chỉnh đồng bộ cả dữ liệu PVT và các hệ số thủy lực của mô hình cơ học, sai số tính toán áp suất đáy giếng và gradient áp suất tại giếng BD-7P giảm mạnh xuống còn 2.8%. Khi áp dụng mô hình đã hiệu chỉnh này sang giếng kiểm chứng BD-1P, sai số áp suất duy trì ổn định ở mức dưới 3.2%, cải thiện độ chính xác dự báo lưu lượng khai thác thêm khoảng 18% so với mô hình mặc định.

Thảo luận kết quả

Độ lệch ban đầu giữa mô hình lý thuyết và thực tế khai thác tại mỏ Báo Đen xuất phát từ tính chất đặc thù của dầu thô bồn trũng Cửu Long, vốn có độ nhớt nhạy cảm với nhiệt độ và chứa hàm lượng sáp paraffin cao. Khi áp suất dòng chảy giảm xuống dưới áp suất điểm bọt khí, khí hòa tan tách ra mạnh mẽ tạo thành dòng chảy hai pha dạng nút với các bọt khí Taylor chuyển động nhanh hơn pha lỏng. Các tương quan thực nghiệm truyền thống như Beggs và Brill hay Hagedorn và Brown không phân tích chi tiết cơ chế trượt pha cục bộ này, dẫn đến ước tính sai tỷ lệ chất lỏng giữ lại trong ống.

Dữ liệu nghiên cứu được biểu diễn sinh động thông qua biểu đồ phân bố chuẩn và biểu đồ Pareto trong phần mềm Minitab, giúp xác định rõ các hệ số thực nghiệm có ảnh hưởng vượt ngưỡng ý nghĩa thống kê 95% đến gradient áp suất. Đồng thời, đường cong phân bố áp suất theo chiều sâu mô phỏng sau hiệu chỉnh bám sát hoàn hảo các điểm đo áp suất thực tế từ thiết bị PLT. Kết quả này chứng minh rằng việc kết hợp mô hình cơ học vật lý với kỹ thuật hiệu chỉnh thống kê thực nghiệm là hướng đi tối ưu để xây dựng đường đặc tính dòng ra chính xác cho giếng dầu khí.

+-------------------------------------------------------------------------------+
| SO SÁNH SAI SỐ MÔ HÌNH DÒNG CHẢY TRƯỚC VÀ SAU HIỆU CHỈNH TẠI GIẾNG KHAI THÁC |
+------------------------------------+--------------------+---------------------+
| Chỉ số đánh giá                    | Mô hình mặc định   | Mô hình hiệu chỉnh  |
+------------------------------------+--------------------+---------------------+
| Sai số tỷ số khí hòa tan (Rs)      | 15.0% - 25.0%      | < 2.5%              |
| Sai số phân bố áp suất giếng BD-7P | 8.5% - 12.3%       | 2.8%                |
| Sai số phân bố áp suất giếng BD-1P | 9.0% - 14.0%       | 3.2%                |
| Tỷ trọng tổn thất do thủy tĩnh     | 80.0% - 95.0%      | 82.0% - 94.0%       |
| Tỷ trọng tổn thất do ma sát        | 5.0% - 20.0%       | 6.0% - 18.0%        |
+------------------------------------+--------------------+---------------------+

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 giải pháp hành động cụ thể nhằm nâng cao hiệu quả vận hành khai thác dầu khí:

  • Chuẩn hóa và tích hợp quy trình hiệu chỉnh hai giai đoạn (hiệu chỉnh PVT trước, hiệu chỉnh thủy lực sau) vào quy trình mô phỏng giếng khai thác trên phần mềm PIPESIM của Công ty Liên doanh Điều hành Cửu Long, đặt mục tiêu duy trì sai số dự báo áp suất đáy giếng dưới 3% ngay trong quý 4 năm 2012.
  • Triển khai áp dụng mở rộng mô hình cơ học đã hiệu chỉnh cho toàn bộ hệ thống giếng khai thác tại các cấu tạo lân cận thuộc bồn trũng Cửu Long như mỏ Sư Tử Đen và Sư Tử Vàng trong vòng 6 đến 12 tháng, phục vụ thiết kế tối ưu đường kính ống nâng từ 3.5 inch đến 4.5 inch.
  • Thiết lập chu kỳ khảo sát đo địa vật lý khai thác PLT định kỳ 6 tháng một lần đối với các giếng có tỷ số khí dầu GOR biến động trên 15% nhằm cập nhật dữ liệu đầu vào và tái hiệu chỉnh các hệ số thủy lực kịp thời.
  • Tổ chức chương trình đào tạo chuyên sâu về ứng dụng công cụ thiết kế thực nghiệm ANOVA và kỹ thuật phân tích Nodal cho đội ngũ kỹ sư công nghệ khai thác tại các đơn vị thành viên của Tập đoàn Dầu khí Việt Nam trong thời hạn 3 tháng.

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

Luận văn là tài liệu chuyên khảo giá trị cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư công nghệ khai thác và kỹ sư mô phỏng mỏ: Nắm vững phương pháp luận và các bước tính toán chi tiết để hiệu chỉnh mô hình dòng chảy nhiều pha, từ đó thiết kế tối ưu hệ thống khai thác tự phun và khai thác bằng gaslift.
  • Cán bộ quản lý vận hành mỏ tại các công ty liên doanh dầu khí: Ứng dụng mô hình dự báo áp suất chính xác để theo dõi sát sao động thái vỉa, lập kế hoạch can thiệp giếng hợp lý và nâng cao hệ số thu hồi dầu của mỏ.
  • 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í: Khai thác tài liệu làm cơ sở tham khảo chuyên sâu về thủy động lực học dòng chảy trong giếng và phương pháp thiết kế thực nghiệm nâng cao trong nghiên cứu dầu khí.
  • Kỹ sư phát triển phần mềm kỹ thuật dầu khí: Tham khảo thuật toán khớp nối số liệu thực tế và quy trình tối ưu hóa các module tính toán PVT kết hợp gradient áp suất giếng khoan.

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

Vì sao mô hình cơ học trong phần mềm PIPESIM vẫn xuất hiện sai số đáng kể khi áp dụng tại mỏ Báo Đen?

Mô hình cơ học dù xây dựng trên các định luật vật lý nhưng vẫn chứa một số hệ số thực nghiệm xác lập từ điều kiện ống nghiệm chuẩn. Khi áp dụng cho dầu thô mỏ Báo Đen có độ nhớt và hàm lượng sáp đặc thù, sai số ban đầu có thể vượt quá 10% nếu chưa được hiệu chỉnh với dữ liệu PVT và PLT cục bộ.

Ưu điểm vượt trội của mô hình cơ học Ansari so với các tương quan thực nghiệm truyền thống là gì?

Mô hình Ansari phân biệt rõ ràng 5 chế độ dòng chảy riêng biệt và tính toán cơ chế thủy động học cho từng pha, đặc biệt là động học của bọt khí Taylor. Điều này giúp phản ánh chính xác hiện tượng trượt pha thay vì coi chất lưu là hỗn hợp đồng nhất như các mô hình thực nghiệm cũ.

Quy trình hiệu chỉnh mô hình dòng chảy hai giai đoạn được thực hiện như thế nào?

Quy trình gồm 2 giai đoạn: trước tiên sử dụng phân tích phương sai ANOVA trên Minitab để hiệu chỉnh 5 thông số PVT khớp với mẫu phòng thí nghiệm, sau đó áp dụng phương pháp bình phương tối thiểu để hiệu chỉnh hệ số ma sát và tỷ lệ bọt khí trong phương trình gradient áp suất.

Thiết bị đo kiểm tra khai thác PLT đóng vai trò gì trong đề tài nghiên cứu?

Thiết bị PLT ghi nhận liên tục áp suất và nhiệt độ thực tế tại các độ sâu khác nhau trong giếng BD-7P và BD-1P. Dữ liệu này đóng vai trò là chuẩn đối sánh thực tế để thuật toán bình phương tối thiểu căn chỉnh đường cong áp suất mô phỏng đạt độ chính xác cao nhất.

Mô hình đã hiệu chỉnh từ giếng BD-7P có thể mở rộng cho các giếng khác trong khu vực không?

Mô hình hoàn toàn có thể áp dụng thành công cho các giếng lân cận trong cùng tầng sản phẩm bồn trũng Cửu Long. Minh chứng thực tế tại giếng BD-1P cho thấy mô hình duy trì sai số áp suất dưới 3.2% mà không cần thay đổi cấu trúc thuật toán cốt lõi.

Kết luận

  • Luận văn đã làm sáng tỏ bản chất thủy động lực học của dòng chảy hai pha trong ống nâng và chứng minh tổn thất áp suất do cột thủy tĩnh chiếm từ 80% đến 95% tổng tổn hao năng lượng trong giếng.
  • Nghiên cứu đã thiết lập thành công quy trình hiệu chỉnh 2 giai đoạn khoa học, kết hợp phân tích phương sai ANOVA trên Minitab và tối ưu hóa bình phương tối thiểu trên PIPESIM.
  • Mô hình hiệu chỉnh giúp giảm sai số PVT xuống dưới 2.5% và đưa sai số tính toán phân bố áp suất tại các giếng BD-7P, BD-1P mỏ Báo Đen về mức dưới 3.2%.
  • Công trình giải quyết triệt để bài toán mô phỏng giếng khai thác tồn đọng tại Công ty Liên doanh Điều hành Cửu Long, tạo tiền đề tối ưu hóa chế độ khai thác cho toàn bồn trũng Cửu Long trong giai đoạn 2012 - 2015.
  • Các kỹ sư và nhà nghiên cứu dầu khí nên ứng dụng ngay quy trình hiệu chỉnh mô hình cơ học này vào công tác thiết kế ống khai thác và dự báo sản lượng mỏ nhằm nâng cao hiệu quả kinh tế và độ an toàn cho các công trình dầu khí biển.