Tiểu luận nghiên cứu, đánh giá hiện tượng ngưng tụ lỏng sử dụng mô hình mô phỏng tích hợp vỉa giếng

Accelerat ing t he world's research Chuyen de 3-NMQ Hoang Vu Related papers Download a PDF Pack of t he best relat ed papers  T P ĐỒ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á hi n tư ng ngưng t l ng sử d ng mơ hình mơ ph ng 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 CL 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 qu 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 qu mơ ph ng 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 a công ngh mô ph ng phục vụ cho ngành dầu khí, vi c mơ ph ng dòng ch y v a ch a lòng gi ng ch y u thực hi n cách riêng r Lý c a vi c dòng ch y v a ch a lòng gi ng hai hi n 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 tồ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 hi n tượng v t lý khác công vi c ph c t p đ i h i kh tính tốn lớn c a máy tính, điều mà kh th m chí hi n t i chưa có phần mềm chuyên dụng gi i quy t 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 quy t vấn đề gặp ph i q 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 b n c a vi c xây dựng mơ hình tích hợp mơ hình v a ch a tính tố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 c a mơ hình v a c n đáy gi ng t i 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ơ ph ng liên k t kiểm sốt tồn b i mơ hình gi ng Trong q 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 tốn lưu lượng dòng ch y c a pha t i mặt phân cách Vì mơ hình tích hợp ch y u sử dụng để tính tố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 tồ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 cho đ n v n cịn i h n ch Ti m c n hướng nghiên c u có cơng bố c a 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 i tăng, ngăn c n dịng khí chuyển động từ v a vào gi ng Đáng ý nghiên c u c a Zhang [6-8] mơ hình động h c v a gi ng, t p trung vào phân tích mơ ph ng 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 s 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 s thay đổi có tác động tr l i vùng c n đáy gi ng Tuy nhiên k t qu c a lĩnh vực sơ khai, phần mềm mô ph ng ch y u d ng ch ng minh kh (proof of concept) Mỗi phần mềm ch dùng để gi i thích vài hi n 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 hi n tượng, tương tác thực t Tốc độ mô ph ng c a phần mềm cịn ch m, ch 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ô ph ng k t hợp hi n xây dựng dựa nguyên lý “black oil”, không phần mềm sử dụng để mơ ph ng hi n tượng condensate banking Để mô ph ng condensate banking, vi c sử dụng mơ hình thành phần điều bắt buộc Tuy nhiên không ph i vi c kh thi mà vi c áp dụng mơ hình thành phần cho c v a ch a lòng gi ng yêu cầu l i gi i lúc c a số lượng lớn phương trình vi phân; điều làm tăng m c độ ph c t p c a toán tăng th i gian tính tốn lên nhiều lần Bên c nh vi c xây dựng phần mềm để mô ph ng lúc c hai ch độ dòng ch y, tồn t i 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 c a phương pháp vi c xây đựng mối quan h áp suất - lưu lượng thư ng ch 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 qu 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 hi n tượng condensate banking v a hi n 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ô ph ng dòng ch y lòng gi ng Olga sử dụng để đánh giá kh ngưng tụ l ng nh hư ng c a đ n ho t động khai thác c a gi ng Mô t chi ti t c a phương pháp ti n hành k t qu nghiên c u s trình bày phần ti p theo c a chuyên đề Phương pháp nghiên c u Để nghiên c u nh hư ng c a condensate banking đ n q 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ơ ph ng thành phần, có kh mơ ph ng tương đối xác q trình condensate banking diễn v a ch a Tuy nhiên E300 khơng có kh mơ ph ng dịch chuyển bi n đổi c a chất lưu đa pha lòng gi ng, tất c 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 n u v n tốc dịng khí khơng đ lớn phần l ng s 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 c a E300, nghiên c u sinh k t hợp k t qu ch y mô ph ng dòng ch y đa phase v a input đầu vào phần mềm mô ph ng dòng ch y gi ng Olga Đây phần mềm mô ph ng 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 c a 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ô ph ng 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 t i th i điểm cần xem xét Sau thơng số s chuyển thành li u đầu vào c a mơ hình Olga để tính tốn l i lưu lượng chất lưu kh nâng l ng lên mi ng gi ng, qua dự đố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 qu 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 c a v a ch a,  Thành phần tính chất PVT c a 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ô ph ng 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 c a mơ hình tính tốn để có độ lớn tương đương với vùng nh hư ng c a gi ng HT-3P theo k t qu thử v a Để tăng độ xác cho k t qu 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ể hi n 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 c a v a ch a 8439 psia đo t i độ sâu 11677 ft với áp suất điểm sương 7091 psia Giá trị ban đầu c a độ thấm ngang thấm d c 4.5 mD 0.45 mD Độ thấm bi n có giá trị chưa chắn s thay đổi q trình khớp lịch sử Các thơng số b n tính chất v a ch a thể hi n 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ô ph ng 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 c a k t qu mơ ph ng trình khai thác, thành phần chất lưu thay đổi nhanh chóng, đặc bi t t i 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 qu 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 c a chất lưu lấy từ gi ng HT-3P thể hi n 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 ch số suất gi ng (well productivity index) Tốc độ dòng ch y c a (a) Lưu lượng khí (b) Lưu lượng condensate 10 (c) Áp suất đáy Hình K t qu 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 c a tầng MMF30 m H i Th ch khớp lịch sử, mơ hình ch y ti p cho đ n 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 qu dự báo động thái khai thác thể hi n Hình Theo k t qu 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í s tr nên q 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ể hi n 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 v n 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 c a chất lưu khai thác (chất lưu từ v a vào gi ng) theo th i gian thể hi n Hình Từ hình v nh n thấy th i gian đầu, tỷ l mole c a 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 b i trình condensate banking, thành phần nặng có xu hướng tồn t i d ng l ng l i v a, ch có phần nhẹ vào gi ng Trong th i gian sau, tỷ l mole c a 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 cịn 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 qu 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 tốn từ k t qu 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 c a 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 qu dự báo từ mơ hình E300, tính giá trị c a C t i th i điểm cần xem xét B ng thống kê h số C t i 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ơ ph ng sử dụng Olga, ch có PR C sử dụng làm 15 thông số đầu vào Q Pwf s tính tố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 q 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 c a 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 tốn dùng làm li u đầu vào cho mơ hình Olga Hình thể hi n bi n đổi gi n đồ pha c a 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 c a 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ơ ph ng 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 qu mô ph ng lưu lượng khí mơ hình Olga so sánh với k t qu từ Eclipse Hình 10 Từ hình v nh n thấy gi ng HT-3P theo mơ ph ng dòng ch y lòng gi ng dừng khai thác sớm nhiều so với dự báo c a Eclipse Điều lý gi i vào giai đo n sau c a đ i m , lưu lượng khí nh n cho v n tốc khí xuống thấp, khơng cịn kh nâng l ng lên mi ng gi ng Mơ hình Eclipse ch xây dựng để mô ph ng dòng ch y c a chất lưu v a ch a nên không mô ph ng ch 17 4,000 400 350 3,000 300 2,500 250 2,000 200 1,500 150 1,000 100 500 ar Lưu lượng khí Ms f/D 3,500 Áp suất vỉa 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ý c a hi n 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 c a 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 c a chất lưu lòng gi ng t i th i điểm đầu năm 2021 mơ t Hình 11 Bắt đầu t i 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 t i d ng màng m ng; lúc khí v n 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 kick- off 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 18 phần c a gi ng 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ể hi n v n tốc khí l ng t i th i điểm 1.7 ngày sau m gi ng (năm 2021) t i 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í v n di chuyển lên trên, v n tốc khí nh , ch m c 0.2 m/s V n tốc khí nh v y khơng thể nâng tồn lượng chất l ng từ đáy gi ng lên mi ng gi ng, thể hi n 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 ti p theo vi c chất l ng bị dồn l i đáy gi ng gây dừng gi ng khơng thể 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á b n chất y u tố nh hư ng tới trình ngưng tụ l ng v a gi ng thực hi n chi ti t phần độc l p c a tiểu lu n Để làm sáng t vấn đề trình v a gi ng nh hư ng tới th nào, hi n tượng ngưng tụ l ng h thống hoàn ch nh v a - gi ng mô ph ng thực hi n đánh giá Các k t qu tổng hợp sau:  Qua đánh giá tài li u tham kh o th giới thực t nghiên c u c a mình, nh n định cho đ n th i điểm hi n t i chưa có phần mềm thương m i hồn ch nh thành cơng mơ ph ng tích hợp v a-gi ng Modun ROCKX c a phần mềm OLGA modun thương m i giới thi u phát triển để mơ ph ng tích hợp Tuy nhiên thực hi n nghiên c u chuyên sâu NCS nh n thấy modun chưa đ lực để mô ph ng đ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 qu mơ hình mơ ph ng v a ch a mơ ph ng dịng ch y gi ng đề xuất thực hi n thành công Các k t qu mô ph ng 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 qu 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ô ph ng k t hợp cho dự báo phù hợp với thực t so với mô ph ng 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 ch đánh giá khai thác thông qua mô ph ng v a ch a cách nhà thầu 21 thực hi n s 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 cịn 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, ch đánh giá dòng ch y lòng gi ng s cho k t qu th i điểm gi ng bị dừng ngưng tụ l ng t i đá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ô ph ng k t hợp v a gi ng theo cách ti p c n mới, k t qu mô ph ng điểm quan tr ng c a 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 ... 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 qu 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... condensate banking v a hi n 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ơ ph ng dịng ch y lòng gi ng Olga sử dụng để đánh giá kh ngưng tụ l ng nh hư ng c a đ n ho... 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 b n c a vi c xây dựng mơ hình tích