TẬP ĐOÀN DẦU KHÍ QUỐC GIA VIỆT NAM VIỆN DẦU KHÍ VIỆT NAM - TIỂU LUẬN Luận án tiến sĩ kỹ thuật "Nghiên cứu, đánh giá tượng ngưng tụ lỏng sử dụng mô hình mô tích hợp vỉa – giếng" Chuyên ngành: Kỹ thuật dầu khí Mã số: 62.52.06.04 Họ tên NCS: Nguyễn Minh Quý Người hướng dẫn 1: TS Phan Ngọc Trung Người hướng dẫn 2: TS Nguyễn Hữu Trung Hà Nội, năm 2017 MỤC LỤC MÔ PHỎNG DÒNG CHẢY TRONG VỈA CHỨA VÀ LÒNG GIẾNG HT-3P2 Phương pháp nghiên cứu Xây dựng mô hình khớp lịch sử 2.1 Cấu trúc ô lưới 2.2 Thông số PVT 2.3 Đường cong thấm pha 2.4 Phục hồi lịch sử 2.5 Dự báo sử dụng phần mềm Eclipse E300 11 Chính xác hóa kết dự báo dùng phần mềm Olga 14 3.1 Phương trình IPR giếng khai thác 14 3.2 Tính chất chất lưu 16 3.3 Kết mô tích hợp vỉa – giếng 17 Tài liệu tham khảo 23 MÔ PHỎNG DÒNG CHẢY TRONG VỈA CHỨA VÀ LÒNG GIẾNG HT3P Trong lịch sử phát triển công nghệ mô phục vụ cho ngành dầu khí, việc mô dòng chảy vỉa chứa lòng giếng chủ yếu thực cách riêng rẽ Lý việc dòng chảy vỉa chứa lòng giếng hai tượng vật lý khác nhau: dòng chảy vỉa chứa trình khuếch tán tuân theo định luật bảo toàn khối lượng định luật Darcy; dòng chảy lòng giếng tuân theo định luật bảo toàn động lượng Việc kết hợp giải lúc phương trình vi phân phần mô tả hai tượng vật lý khác công việc phức tạp đỏi hỏi khả tính toán lớn máy tính, điều mà khứ chí chưa có phần mềm chuyên dụng giải triệt để Việc tìm hiểu tương tác lưu thể vỉa chứa lòng giếng trở nên cần thiết để giải vấn đề gặp phải trình khai thác, ví dụ slugging, lưỡi nước, ngưng tụ lỏng condensate banking Mô hình tích hợp xây dựng cách kết hợp mô hình vỉa chứa khu vực cận đáy giếng mô hình dòng chảy đa pha lòng giếng Nguyên tắc việc xây dựng mô hình tích hợp mô hình vỉa chứa tính toán hệ số độ nhạy cho sản lượng khai thác có tính đến áp suất giếng thiết lập sẵn từ mô hình dòng chảy giếng Tiếp theo, mô hình giếng sử dụng hệ số độ nhạy để giải ngược lại áp suất giếng Hệ số độ nhạy xuất từ ma trận Jacobian mô hình vỉa cận đáy giếng bước lặp cuối Mô hình vỉa xem xét phần tích hợp thêm (plugin) tới mô hình giếng, mô hình mô liên kết kiểm soát toàn mô hình giếng Trong trình mô phỏng, mô hình giếng tạm thời cung cấp ranh giới áp suất tới mô hình mỏ mô hình mỏ tính toán lưu lượng dòng chảy pha mặt phân cách Vì mô hình tích hợp chủ yếu sử dụng để tính toán chuyển tiếp dòng-áp suất tức thời thân giếng khu vực vỉa cận giếng, không cần thiết phải sử dụng toàn mô hình vỉa vùng xa giếng khai thác vỉa chứa có ảnh hưởng đến chuyển tiếp dòng-áp suất tức thời thân giếng khu vực vỉa cận giếng Nghiên cứu mối tương tác dông chảy vỉa giếng ỏi hạn chế Tiệm cận hướng nghiên cứu có công bố Sturn [1] vào năm 2004, Ali [2] Dousi [3] vào năm 2005, Chupin [4] vào năm 2007, Al-Darmaki [5] năm 2008 Zhang [6], [7] & [8] Hầu hết nghiên cứu mối quan hệ vỉa giếng tập trung vào trạng thái dòng chảy tức thời vùng cận đáy giếng giếng Khi áp suất vùng cận đáy giếng giảm xuống giá trị áp suất bão hòa, lượng pha lỏng vùng đáy giếng tăng lên làm áp suất thủy tĩnh tăng, ngăn cản dòng khí chuyển động từ vỉa vào giếng Đáng ý nghiên cứu Zhang [6-8] mô hình động học vỉa giếng, tập trung vào phân tích mô trạng thái dòng chảy áp suất đáy giếng thay đổi Zhang nhận thấy tương tác động học vỉa giếng trình liên tục trình thay đổi xảy vùng cận đáy giếng có ảnh hưởng tức trình vận chuyển lưu chất giếng khai thác, ngược lại chế độ dòng chảy giếng thay đổi mà diễn liquid loading, áp suất đáy giếng thay đổi có tác động trở lại vùng cận đáy giếng Tuy nhiên kết lĩnh vực sơ khai, phần mềm mô chủ yếu dạng chứng minh khả (proof of concept) Mỗi phần mềm dùng để giải thích vài tượng đơn giản, chưa có tác dụng tổng quát áp dụng để giải thích nhiều tượng, tương tác thực tế Tốc độ mô phần mềm chậm, thích hợp với mô hình đơn giản có kích thước nhỏ, khó để áp dụng cho mô hình thực tế với kích thước hàm trăm nghìn hay hàng triệu ô lưới Đặc biệt, phần mềm mô kết hợp xây dựng dựa nguyên lý “black oil”, không phần mềm sử dụng để mô tượng condensate banking Để mô condensate banking, việc sử dụng mô hình thành phần điều bắt buộc Tuy nhiên việc khả thi mà việc áp dụng mô hình thành phần cho vỉa chứa lòng giếng yêu cầu lời giải lúc số lượng lớn phương trình vi phân; điều làm tăng mức độ phức tạp toán tăng thời gian tính toán lên nhiều lần Bên cạnh việc xây dựng phần mềm để mô lúc hai chế độ dòng chảy, tồn phương pháp đơn giản để nghiên cứu tương tác lòng giếng vỉa chứa (Oudeman [9]) Những phương pháp chủ yếu dựa việc kết hợp động thái khai thác đa pha với mối quan hệ áp suất - lưu lượng ống khai thác (vertical flow performance) để dự báo trình ngưng tụ lỏng điều kiện chuyển tiếp Điểm hạn chế phương pháp việc xây đựng mối quan hệ áp suất - lưu lượng thường dựa hàm quan hệ thay mô hình mô phỏng, ví dụ Hargedon & Brown, Beggs & Brill Các hàm quan hệ thường cho kết chấp nhận áp dụng cho điều kiện bình thường cho sai số lớn áp dụng cho điều kiện phức tạp ngưng tụ lỏng Trong chuyên đề này, nghiên cứu sinh nghiên cứu tương tác tượng condensate banking vỉa tượng ngưng tụ lỏng lòng giếng theo hướng khác Thay sử dụng hàm quan hệ, phần mềm mô dòng chảy lòng giếng Olga sử dụng để đánh giá khả ngưng tụ lỏng ảnh hưởng đến hoạt động khai thác giếng Mô tả chi tiết phương pháp tiến hành kết nghiên cứu trình bày phần chuyên đề Phương pháp nghiên cứu Để nghiên cứu ảnh hưởng condensate banking đến trình ngưng tụ lỏng lòng giếng, tác giả sử dụng kết hợp phần mềm Eclipse (E300) phần mềm Olga E300 phần mềm mô thành phần, có khả mô tương đối xác trình condensate banking diễn vỉa chứa Tuy nhiên E300 khả mô dịch chuyển biến đổi chất lưu đa pha lòng giếng, tất chất lưu vào giếng coi có khả lên tới miệng giếng Với cách làm vậy, E300 bỏ qua tương tác pha lỏng khí lòng giếng Trong thực tế vận tốc dòng khí không đủ lớn phần lỏng không nâng lên tới miệng giếng, đọng lại đáy giếng gây dừng hoạt động giếng Để bổ sung cho thiếu hụt E300, nghiên cứu sinh kết hợp kết chạy mô dòng chảy đa phase vỉa input đầu vào phần mềm mô dòng chảy giếng Olga Đây phần mềm mô thành phần tương tác chất lưu lòng giếng Phần mềm có khả tính đến biến đổi tính chất biến đổi pha chất lưu thay đổi áp suất, nhiệt độ dọc theo chiều dài thân giếng (từ khu vực nhiệt độ cao, áp suất cao lên khu vực nhiệt độ thấp, áp suất thấp) Phương thức nghiên cứu sinh sử dụng để kết hợp E300 Olga sau: mô hình E300 sử dụng mô lấy thông số dự báo áp suất vỉa chứa, lưu lượng thành phần chất lưu thời điểm cần xem xét Sau thông số chuyển thành liệu đầu vào mô hình Olga để tính toán lại lưu lượng chất lưu khả nâng lỏng lên miệng giếng, qua dự đoán thời điểm dừng giếng Quy trình cụ thể gồm bước sau: Bước – Xây dựng mô hình E300 khớp lịch sử Bước – Dự báo động thái khai thác sử dụng mô hình E300 khớp lịch sử Bước – Tại thời điểm cần xem xét, sử dụng kết mô hình E300 làm số liệu đầu vào cho mô hình Olga chạy mô hình Các liệu cần phải cung cấp cho mô hình Olga gồm có:  Áp suất vỉa chứa,  Phương trình IPR vỉa chứa,  Thành phần tính chất PVT chất lưu, Áp suất, nhiệt độ dọc theo chiều dài thân giếng Xây dựng mô hình khớp lịch sử 2.1 Cấu trúc ô lưới Để mô dự báo động thái khai thác giếng HT-3P, mô hình 3D cho tầng MMF30 xây dựng sử dụng phần mềm Eclipse E300 Mô hình có kích cỡ 20 x 20 x với kích thước ô lưới 100 ft x 100 ft x 20 ft Thể tích mô hình tính toán để có độ lớn tương đương với vùng ảnh hưởng giếng HT-3P theo kết thử vỉa Để tăng độ xác cho kết mô phỏng, cấu trúc ô lưới chia nhỏ (local grid refinement - LGR) thiết lập vùng xung quanh giếng HT-3P Cấu trúc LGR bao phủ vùng bán kính 250 ft quanh giếng với kích thước ô lưới 10 ft x 10 ft x 20 ft Cấu trúc LGR thứ hai bao phủ vùng bán kính ft quanh giếng với bán kính ô lưới ft x ft x 20 ft Cấu trúc ô lưới thể Hình Hình Cấu trúc ô lưới mô hình mô tầng MMF30 mỏ Hải Thạch Áp suất ban đầu vỉa chứa 8439 psia đo độ sâu 11677 ft với áp suất điểm sương 7091 psia Giá trị ban đầu độ thấm ngang thấm dọc 4.5 mD 0.45 mD Độ thấm biến có giá trị chưa chắn thay đổi trình khớp lịch sử Các thông số tính chất vỉa chứa thể Bảng Bảng Tính chất vỉa chứa mô hình 2.2 Đại lượng Đơn vị Giá trị Độ thấm ngang mD 4.5 Độ thấm dọc mD 0.45 Độ rỗng % 14.795 Áp suất vỉa ban đầu psia 8439 Độ sâu quy chiếu ft 11677 Chiều dày hiệu dụng ft 80 Thông số PVT Khi mô mỏ khí condensate, đặc biệt sử dụng mô hình thành phần, lựa chọn thông số PVT có ý nghĩa tối quan trọng tới độ xác kết mô trình khai thác, thành phần chất lưu thay đổi nhanh chóng, đặc biệt khu vực cận đáy giếng, dẫn tới thay đổi tính chất chất lưu tỉ trọng, độ nhớt, hệ số thể tích dầu tương tác đa pha Bởi lý trên, mô hình PVT sử dụng phương trình trạng thái PengRobinson xây dựng dựa kết nghiên cứu thực nghiệm chất lưu lấy từ giếng HT-3P Thành phần giản đồ pha chất lưu lấy từ giếng HT-3P thể Bảng Hình Bảng Thành phần chất lưu tầng MMF30 mỏ Hải Thạch Hình Giản đồ pha chất lưu tầng MMF30 mỏ Hải Thạch 2.3 Đường cong thấm pha Đối với vỉa chứa đa pha, đặc biệt với mỏ khí condensate có điều kiện phức tạp Hải Thạch, thấm pha có ảnh hưởng lớn đến tính chất dòng chảy số suất giếng (well productivity index) Tốc độ dòng chảy khí lỏng vùng cận đáy giếng phụ thuộc nhiều vào thay đổi thấm pha gây nên trình condensate banking Các đường thấm pha sử dụng mô hình mô xây dựng từ liệu phân tích mẫu đặc biệt giếng HT-3P (Hình 3) Hình Thấm pha chất lưu tầng MMF30 mỏ Hải Thạch 2.4 Phục hồi lịch sử Quá trình khớp lịch sử thực để bảo đảm mô hình mô phản ánh động thái khai thác vỉa Để khớp lịch sử, thông số độ thấm số suất giếng thay đổi Kết khớp lịch sử thể Hình Kết khớp lịch sử tốt làm tăng độ tin cậy khả dự báo mô hình (a) Lưu lượng khí (b) Lưu lượng condensate 10 (c) Áp suất đáy Hình Kết khớp lịch sử giếng HT-3P 2.5 Dự báo sử dụng phần mềm Eclipse E300 Sau mô hình E300 tầng MMF30 mỏ Hải Thạch khớp lịch sử, mô hình chạy tiếp hết năm 2036 để dự báo động thái khai thác giếng HT-3P Các thông số khống chế giai đoạn dự báo gồm có:  Áp suất đáy tối thiểu: 800 psia  Lưu lượng khí tối đa: 3390.21 Mscf/D Với thông số trên, kết dự báo động thái khai thác thể Hình Theo kết dự báo E300 giếng HT-3P khai thác đến cuối năm 2036 áp suất vỉa xuống thấp (~1000 psia) Dự báo có xu hướng lạc quan cần kiểm chứng lại phần mềm Olga nhiều khả vận tốc khí trở nên thấp để nâng phần lỏng lên miệng giếng 11 Hình Dự báo động thái khai thác giếng HT-3P Hồng: áp suất vỉa, đỏ: áp suất đáy, cam: lưu lượng khí, lục: lưu lượng condensate Sự thay đổi áp suất bão hòa condensate vỉa chứa thể Hình Từ hình vẽ nhận thấy với sụt giảm áp suất vỉa chứa theo thời gian, lượng condensate vỉa chứa giảm dần Điều chứng tỏ trình ngưng tụ ngược kết thúc sau thời gian dự báo bắt đầu Tuy nhiên, độ bão hòa condensate vỉa chứa trì mức cao (20%38%), gây khó khăn cho trình sản xuất (a) 1/1/2016 (b) 1/1/2018 12 (c) 1/1/2020 (d) 1/1/2022 (e) 1/1/2024 (f) 1/1/2026 Hình Sự thay đổi áp suất vỉa chứa theo thời gian (a) 1/1/2016 (b) 1/1/2018 (c) 1/1/2020 (d) 1/1/2022 (e) 1/1/2024 (f) 1/1/2026 Hình Sự thay đổi bão hòa condensate theo thời gian 13 Sự biến đổi thành phần chất lưu khai thác (chất lưu từ vỉa vào giếng) theo thời gian thể Hình Từ hình vẽ nhận thấy thời gian đầu, tỷ lệ mole C1 tăng dần theo thời gian, đồng nghĩa với việc chất lưu ngày trở nên nhẹ Điều lý giải trình condensate banking, thành phần nặng có xu hướng tồn dạng lỏng lại vỉa, có phần nhẹ vào giếng Trong thời gian sau, tỷ lệ mole C1 giảm dần cho trình ngưng tụ ngược kết thúc, thành phần nặng không lại vỉa mà vào giếng nhiều Hình Sự biến đổi thành phần chất lưu theo thời gian Đỏ: C1, Lục: C2, C3, CO2 Chính xác hóa kết dự báo dùng phần mềm Olga 3.1 Phương trình IPR giếng khai thác Phương trình IPR sử dụng làm đầu vào cho mô hình Olga tính toán từ kết mô hình Eclipse Trong nghiên cứu tác giả sử dụng phương trình IPR dạng back-pressure, dạng phương trình sử dụng phổ biến cho mỏ khí Công thức chung phương trình IPR có dạng sau: 14 Q = C (PR2 – Pwf2), với  Q: lưu lượng khí khai thác (scf/D)  PR: áp suất vỉa chứa (psia)  Pwf: áp suất đáy (psia)  C: hệ số khả khai thác (scf/(D.psia2)) Ở đây, C tham số chưa biết thay đổi theo thời gian Sử dụng số liệu Q, PR Pwf từ kết dự báo từ mô hình E300, tính giá trị C thời điểm cần xem xét Bảng thống kê hệ số C thời điểm khác thời gian dự báo Bảng Sự biến thiên phương trình IPR theo thời gian Time qg (Mscf/D) PR (WBP9, psia) Pwf (psia) C 1/1/2016 3,390.2 5,241.5 2,657.1 0.166 1/1/2017 3,390.2 4,433.1 2,070.2 0.221 1/1/2018 3,390.2 3,680.2 1,323.9 0.288 1/1/2019 2,990.6 2,983.8 800.0 0.362 1/1/2020 2,125.2 2,476.1 800.0 0.387 1/1/2021 1,544.1 2,120.4 800.0 0.400 1/1/2022 1,151.8 1,864.0 800.0 0.406 1/1/2023 879.7 1,672.2 800.0 0.408 1/1/2024 686.6 1,525.3 800.0 0.407 1/1/2025 546.3 1,410.2 800.0 0.405 1/1/2026 441.7 1,318.2 800.0 0.402 Một điểm quan trọng cần ý mô hình Olga, lưu lượng kiểm soát áp suất miệng giếng thay áp suất đáy giếng mô hình Eclipse Do mô sử dụng Olga, có PR C sử dụng làm 15 thông số đầu vào Q Pwf tính toán lại dựa tương tác khí lỏng xảy lòng giếng 3.2 Tính chất chất lưu Trong trình khai thác, thành phần chất lưu từ vỉa vào giếng thay đổi theo thời gian (Hình 45), tính chất chất lưu thay đổi Sử dụng phần mềm PVT Sim, thay đổi giản đồ pha ứng với thay đổi thành phần chất lưu tính toán dùng làm liệu đầu vào cho mô hình Olga Hình thể biến đổi giản đồ pha chất lưu theo thời gian Từ hình vẽ nhận thấy so với năm 2016, diện tích khu vực hai pha năm sau nhỏ Điều chứng tỏ khí vào giếng chứa nhiều thành phần nhẹ khô hơn, phù hợp với nhận định phần trước (a) 1/1/2016 (b) 1/1/2018 (c) 1/1/2020 (d) 1/1/2022 16 (e) 1/1/2024 (f) 1/1/2026 Hình Sự thay đổi giản đồ pha chất lưu khai thác theo thời gian 3.3 Kết mô tích hợp vỉa – giếng Mô hình mô giếng chạy định kỳ năm lần thời gian dự báo để kiểm tra khả nâng lỏng lòng giếng Phần mềm Olga sử dụng liệu áp suất vỉa phương trình IPR từ mô hình Eclipse để xác hóa lại lưu lượng khí lòng giếng Kết mô lưu lượng khí mô hình Olga so sánh với kết từ Eclipse Hình 10 Từ hình vẽ nhận thấy giếng HT-3P theo mô dòng chảy lòng giếng dừng khai thác sớm nhiều so với dự báo Eclipse Điều lý giải vào giai đoạn sau đời mỏ, lưu lượng khí nhỏ khiến cho vận tốc khí xuống thấp, không khả nâng lỏng lên miệng giếng Mô hình Eclipse xây dựng để mô dòng chảy chất lưu vỉa chứa nên không mô chế 17 4,000 400 Lưu lượng khí (Mscf/D) 3,500 350 3,000 300 2,500 250 2,000 200 1,500 150 1,000 100 500 Áp suất vỉa (bar) qg E300 qg Olga PR 50 0 Hình 10 So sánh sản lượng khí dự báo dùng Olga Eclipse E300 Một điểm đáng ý khác vào thời điểm dừng giếng, áp suất vỉa vào khoảng 150 bar; giá trị thấp giá trị dừng giếng mô hình giếng riêng lẻ (230 bar) Lý tượng mô hình giếng riêng lẻ không tính đến thay đổi tính chất chất lưu theo thời gian chất lưu ngày trở nên nhẹ khô Việc cho thấy tính cần thiết việc kết hợp mô hình vỉa chứa mô hình lòng giếng để có đánh giá xác tình trạng giếng Để hiểu rõ trình dừng giếng, chuyển động chất lưu lòng giếng thời điểm đầu năm 2021 mô tả Hình 11 Bắt đầu thời điểm ngày, giếng đóng giếng không lỏng Sau giếng mở hai pha lỏng, khí bắt đầu vào giếng Tại thời điểm ngày, lượng chất lỏng giếng chưa nhiều tồn dạng màng mỏng; lúc khí chảy tự từ đáy giếng lên miệng giếng Tại thời điểm ngày, lượng chất lỏng giếng bắt đầu nhiều lên tập trung gần khu vực kickoff Lúc chất lỏng chảy theo dạng slug gây ảnh hưởng nhiều đến lưu lượng khí Từ thời điểm ngày trở đi, chất lỏng dồn lại phần giếng 18 khí phải chảy chế độ dòng chảy bọt xuyên qua cột chất lỏng Đến thời điểm ngày, giếng dừng khai thác hoàn toàn (a) ngày (b) ngày (c) ngày (d) ngày (b) ngày (c) ngày Hình 11 Quá trình chết giếng thời điểm đầu năm 2018 Hình 12 thể vận tốc khí lỏng thời điểm 1.7 ngày sau mở giếng (năm 2021) vị trí khau dọc theo chiều dài thân giếng Từ hình vẽ nhận thấy vào lúc này, khí di chuyển lên trên, vận tốc khí nhỏ, mức 0.2 m/s Vận tốc khí nhỏ nâng toàn lượng chất lỏng từ đáy giếng lên miệng giếng, thể việc lượng chất lỏng xuống (vận tốc < 0) nhiều lượng chất lỏng lên (vận tốc > 0) Quá trình tiếp 19 tục thời điểm việc chất lỏng bị dồn lại đáy giếng gây dừng giếng tránh khỏi Hình 12 Vận tốc dầu khí thời điểm 1.7 ngày sau mở giếng (năm 2021) 20 KẾT LUẬN Việc đánh giá chất yếu tố ảnh hưởng tới trình ngưng tụ lỏng vỉa giếng thực chi tiết phần độc lập tiểu luận Để làm sáng tỏ vấn đề trình vỉa giếng ảnh hưởng tới nào, tượng ngưng tụ lỏng hệ thống hoàn chỉnh vỉa - giếng mô thực đánh giá Các kết tổng hợp sau:  Qua đánh giá tài liệu tham khảo giới thực tế nghiên cứu mình, nhận định thời điểm chưa có phần mềm thương mại hoàn chỉnh thành công mô tích hợp vỉa-giếng Modun ROCKX phần mềm OLGA modun thương mại giới thiệu phát triển để mô tích hợp Tuy nhiên thực nghiên cứu chuyên sâu NCS nhận thấy modun chưa đủ lực để mô điều kiện mỏ thực tế Khi liên hệ, nhà sản xuất phần mềm SPT thừa nhận thông báo cố gắng khắc phục, phát triển  Giải pháp tích hợp theo hướng mới: kết hợp kết mô hình mô vỉa chứa mô dòng chảy giếng đề xuất thực thành công Các kết mô vỉa chứa trạng thái pha, thành phần chất lưu động thái lượng vùng cận đáy giếng…được sử dụng làm số liệu đầu vào để xây dựng mô hình dòng chảy giếng  Kết cho thấy với việc tính đến thay đổi động thái pha, thành phần chất lưu động thái lượng vỉa chứa lòng giếng suốt đời mỏ, mô kết hợp cho dự báo phù hợp với thực tế so với mô riêng lẻ vỉa chứa lòng giếng  Kết luận quan trọng cụ thể cho mỏ Hải Thạch rút ra: Khi đánh giá khai thác thông qua mô vỉa chứa cách nhà thầu 21 thực cho tổng sản lượng dự báo khai thác cao thực tế với thời điểm giếng không khả cho dòng mô hình muộn thực tế không tính đến trình ngưng tụ lỏng đáy giếng Ngược lại, đánh giá dòng chảy lòng giếng cho kết thời điểm giếng bị dừng ngưng tụ lỏng đáy giếng sớm so với thực tế không tính đến thay đổi thành phần chất lưu vào giếng ngày khô phần nặng vỉa  Với việc mô kết hợp vỉa giếng theo cách tiếp cận mới, kết mô điểm quan trọng nghiên cứu giúp cho việc dự báo xác động thái khai thác mỏ Hải Thạch nói riêng mỏ khí-condensate nói chung 22 Tài liệu tham khảo Sturn W.L., Belfroid S.P.C., Van Wolfswinkel O., Peters M.C.A.M and Verhelst F.J.P.C.M.G., "Dynamic reservoir well interaction," SPE paper 90108-MS presented at the SPE annual Technical Conference and Exhibition held in Houston TX, USA, 26-29 September 2004 Ali A.M., Falcone G., Hewitt G.F., Bozorgzadeh M., Gringarten A.C., "Experimental investigation of wellbore phase redistribution effects on pressure-transient data," SPE paper 96587-MS presented at by an SPE program Committee following review of information contained in a proposal by the authors, 2005 Dousi N., Veeken C.A.M., and Curri P.K., "Numerical and analytical modeling of the gas well liquid loading process," SPE paper 95282MS presented at the SPE Offshore Europe Aberdeen, 6-9 September 2005 Chupin G., Hu B., Haugset T an Claudel M., "Intergrated wellbore/reservoir models predicts flow transient in liquid loading gas wells," SPE paper 110461-MS presented at the SPE annual Technical Conference and Exhibition held in Anaheim, CA, USA, 2007 Al-Darmaki S., Falcone G., Hale C.P and Hewitt G.F., "Experimental investigation and modeling the effects of rising gas bubbles in a closed pipe," SPE paper 103129-PA presented at the SPE annual Technical Conference and Exhibition held in San Antonio, TX, USA, 24-27 September 2006 He Zhang, Gioia Falcone, Valko P., Catalin Teodoriu, "Numerical modeling of fully transient two phase flow in the near wellbore region during liquid loading in gas well," SPE paper 122785-MS presented at the SPE LACPEC conference held in Cartagena, Colombia, 31 May – June 2009 He Zhang, Gioia Falcone, Catalin Teodoriu, "Modeling fully transient two phase flow in the near wellbore region during liquid loading in gas 23 well," Journal of natural gas science and engineering, vol (2010), pp 122131, 2010 He Zhang, Gioia Falcone, Catalin Teodoriu, "Relative permeability hysteresis effects in near-wellbore region during liquid loading in gas well," SPE paper 139062 presented at the SPE Latin American & Caribean Petroleum engineering conference held in Lima, Peru, 1-3 December 2010 Pieter Oudeman, “Improved Performance”, SPE paper 19103, 1990 24 Prediction of Wet-Gas-Well ... lượng khí, lục: lưu lượng condensate Sự thay đổi áp suất bão hòa condensate vỉa chứa thể Hình Từ hình vẽ nhận thấy với sụt giảm áp suất vỉa chứa theo thời gian, lượng condensate vỉa chứa giảm dần... analytical modeling of the gas well liquid loading process," SPE paper 95282MS presented at the SPE Offshore Europe Aberdeen, 6-9 September 2005 Chupin G., Hu B., Haugset T an Claudel M., "Intergrated... chứa đa pha, đặc biệt với mỏ khí condensate có điều kiện phức tạp Hải Thạch, thấm pha có ảnh hưởng lớn đến tính chất dòng chảy số suất giếng (well productivity index) Tốc độ dòng chảy khí lỏng vùng