Nghiên cứu và tính toán công nghệ vận chuyển khí đồng hành từ mỏ bạch hổ về bờ

ĐAI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA _oOo NGHIÊN CỨU & TÍNH TOÁN CÔNG NGHỆ VẬN CHUYỂN KHÍ ĐỒNG HÀNH TỪ MỎ BẠCH HỔ VỀ BỜ CHUYÊN NGÀNH: CÔNG NGHỆ HÓA HỌC MÃ SỐ NGÀNH: 10 00 Tp HỒ CHÍ MINH, THÁNG NĂM 2005 CÔNG TRÌNH ĐƯC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: PGS TS Nguyễn Thanh Nam Cán chấm nhận xét 1: Cán chấm nhận xét 2: Luận văn thạc só bảo vệ tại: HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA Ngày…….tháng… năm 2005 TRƯỜNG DẠI HỌC BÁCH KHOA Phòng Đào tạo SĐH Cộng hòa xã hội chủ nghóa Việt nam Độc lập- Tự do- Hạnh Phúc oOo— Tp HCM, ngày tháng năm 2005 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: Vũ Tú Hoài Phái : Ngày tháng năm sinh : 10.07.1967 Nơi sinh: Nam định Chuyên ngành: Công nghệ Hóa học MSHV nam TÊN ĐỀ TÀI: Nghiên cứu tính toán công nghệ vận chuyển khí đồng hành từ mỏ Bạch hổ bờ I NHIỆM VỤ VÀ NỘI DUNG Tính toán suy giảm áp lực dọc tuyến theo phương pháp khác So sánh kết qủa tính với kết qủa chạy phần mềm Hysys 2000 vàsố liệu tuyến ống thực tế • Rút phương pháp tính phù hợp cho tuyến ống dẫn khí thực tế II NGÀY GIAO NHIỆM VU :Ï……………………………………………………………………… • • III NGÀY HOÀN THÀNH NHIỆM VỤ :…………………………………………………… IV HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS TS Nguyễn Thanh Nam CÁN BỘ HƯỚNG DẪN CN BỘ MÔN QL CHUYÊN NGÀNH PGS.TS Nguyễn Thanh Nam Nội dung đề cương luận văn thạc só Hội đồng chuyên ngành thông qua Ngày TRƯỞNG PHÒNG ĐT- SĐH tháng năm 2005 TRƯỞNG KHOA QUẢN LÝ NGÀNH TÓM TẮT LUẬN VĂN Nội dung luận văn tập chung vào việc tính toán công nghệ vận chuyển cho tuyến ống dẫn khí từ mỏ Bạch hổ bờ Tính toán công nghệ vận chuyển dựa sở phương trình dòng Các phương trình xây dựng từ định luật bảo toàn lượng Trong trình xây dựng phương trình dòng, số giả thuyết đưa vào dẫn tới kết qủa tính bị sai lệch Để khắc phục hạn chế này, nhiều phương trình dòng thực nghiệm đề xuất Hầu hết phương trình dòng thực nghiệm có nguồn gốc từ phương trình dòng sở, khác khau việc xử lý hệ số ma sát Các phương trình dòng thực nghiệm sử dụng để tính toán công nghệ cho tuyến ống dẫn khí thực tế luận văn Từ việc phân tích nguyên nhân dẫn đến kết qủa tính thông số công nghệ bị sai lệch kết hợp với kết qủa chạy phần mềm Hysys so sánh kết qủa với số liệu thực tế Từ đó, phương pháp tính phù hợp tìm tuyến ống dẫn khí từ mỏ Bạch hổ bờ LỜI CÁM ƠN Với giúp đỡ nhiệt tình thầy cô bạn, luận văn hoàn thành Lời cảm ơn xin gửi tới thầy PGS.TS Nguyễn Thanh Nam, người bỏ công sức hướng dẫn nhiệt tình chu đáo Ngoài xin gửi lời cảm ơn tới thầy cô giáo thuộc Khoa Hóa Công nghệ Dầu khí, ĐHBK TP HCM đóng góp ý kiến qúy giá kiến thức chuyên môn luận văn Tôi xin chân thành cảm ơn anh chị đồng nghiệp bạn công tác phòng ban chức góp ý chân thành đóng góp kinh nghiệm thực tiễn qúy giá trình thực luận văn: • Phòng Sắc ký khí - Viện CNKH & TK, Vietsovpetro • Phòng Kỹ thuật - XN khí, Vietsovpetro • Ban Lãnh đạo Giàn Nén Lớn, Vietsovpetro • Ban lãnh đạo, Phòng Kỹ thuật – XN Vận chuyển , Công ty Chế biến Kinh doanh Sản phẩm khí • Phòng Kỹ thuật - Nhà máy Xử lý khí Dinh cố • Phòng Đào tạo - Trung tâm Đào tạo Nhân lực Dầu khí Một lần xin chân thành cảm ơn tất cả! -1- MỤC LỤC Trang Lời nói đầu .4 Chương TỔNG QUAN VỂ KHÍ ĐỒNH HÀNH .5 VÀ CÔNG NGHỆ VẬN CHUYỂN KHÍ ĐỒNG HÀNH Giới thiệu Khí đồng hành .6 Phương trình dòng Tham số đầu vào 13 Nhiệm vụ luận văn 14 Chương 15 CƠ SỞ TÍNH TOÁN .15 CÔNG NGHỆ VẬN CHUYỂN KHÍ ĐỒNG HÀNH .15 Phương trình dòng 15 Áp dụng phương trình dòng cho đường ống .19 2.1 2.2 2.3 2.4 Các tham số liên quan phương trình dòng 19 Phương trình dòng đường ống nằm ngang 24 Phương trình dòng thực nghiệm đường ống nằm ngang 28 Hiệu chỉnh phương trình dòng pha cho dòng hai pha 30 Vận tốc cho phép đường ống dẫn khí 39 Tính chất khí 40 4.1 Khí lý tưởng 41 4.2 Khí thực .42 4.2.1 Các phương trình trạng thái khí thực 42 4.2.2 Nhiệt độ áp suất tới hạn 46 HV: Vũ Tú Hoài -2- 4.2.3 4.2.4 Xác định hệ số nén Z 48 Độ nhớt khí 51 Các tham số vật lý liên quan khác 54 5.1 5.2 5.3 5.4 Khối lượng riêng dung dịch lỏng 54 Hệ số hoà tan khí 57 Hệ số quy đổi thể tích 58 Sức căng bề mặt 59 Chương 61 THÔNG SỐ ĐẦU VÀO VÀ ĐIỀU KIỆN BIÊN .61 Thông số tuyến ống dẫn khí 61 Các thông số công nghệ .63 2.1 2.2 Phân tích mẫu phòng thí nghiệm 63 Thông số công nghệ .65 Chương 69 XÂY DỰNG CHƯƠNG TRÌNH 69 TÍNH TOÁN CÔNG NGHỆ VẬN CHUYỂN 69 Giới thiệu .69 Trình tự tính tổn thất áp lực .69 Kết qủa tính áp lực dọc tuyến .72 3.1 3.2 3.3 Tính theo phương trình dòng pha 72 Tính theo phương trình doøng hai pha 74 Tính tổn thất áp lực đầu theo phần mềm Hysys .78 Chương 82 PHÂN TÍCH KẾT QỦA TÍNH VÀ ỨNG DỤNG THỰC TẾ 82 Tổng hợp kết qủa tính 82 Phân tích kết qủa tính 83 Phương pháp tính sửa đổi 87 HV: Vũ Tú Hoài -3- Ứng dụng thực tế 88 KẾT LUẬN 91 TÀI LIỆU THAM KHẢO 93 PHUÏ LỤC A PHỤ LỤC B HV: Vũ Tú Hoài -4- Lời nói đầu Thăm dò khai dầu khí lónh vực đầy tiềm nước ta Để phát triển ngành cách hiệu qủa kinh tế, việc đầu tư cho nghiên cứu công nghệ cần thiết Đối với ngành dầu khí, khâu vận chuyển không phần quan trọng, ví dụ như: đường ống dẫn khí Lan Côn Sơn , PM3-Cà mau, Bạch hổ -Thủ đức, giàn nhiều dự án khác thực tương lai không xa Vận chuyển dầu khí thường trường hợp sau: dẫn dòng pha, dòng hai pha dòng ba pha Hiện dự án kể chủ yếu trường hợp vận chuyển với dòng pha hai pha, dòng ba pha thường đường ống nội mỏ Cho tới giới có nhiều công trình nghiên cứu công bố việc tính toán công nghệ vận chuyển dầu khí với phương pháp rút từ thực nghiệm, áp dụng phương pháp trường hợp cụ thể (kể phần mềm thương mại nay) kết qủa nhận khác khác so với thực tế Vậy nguyên nhân sai khác gì? vấn đề xem xét chi tiết luận văn Trong qúa trình tính toán công nghệ vận chuyển, cần xác định thông số công nghệ (ví dụ áp suất đầu ra), có nhiều tương quan khác đồng thời áp dụng, từ rút tương quan phù hợp Từ mối quan hệ ảnh hưởng qua lại thông số công nghệ, rút nên tác động vào thông số công nghệ (cho trường hợp cụ thể) để có chế độ công nghệ hiệu qủa Phân tích kết qủa tính toán nhận cho khả dự báo tình công nghệ xấu xảy để đề phòng ngăn chặn Kết qủa nghiên cứu tính toán kiểm nghiệm áp dụng cho tuyến ống dẫn khí thực tế từ mỏ Bạch hổ bờ Tuy vấn đề tính toán công nghệ vận chuyển giới mới, ngành dầu khí mẻ Thông qua việc nghiên cứu đề tài này, tác giả muốn đặt tiền đề cho bước vào hướng nghiên cứu công nghệ vận chuyển ngành dầu khí mà nhiều vấn đề thực tế chưa giải HV: Vũ Tú Hoài -5- Chương TỔNG QUAN VỂ KHÍ ĐỒNH HÀNH VÀ CÔNG NGHỆ VẬN CHUYỂN KHÍ ĐỒNG HÀNH Giới thiệu Để công việc vận chuyển dầu khí an toàn hiệu qủa, điều cần phải quan tâm công nghệ vận chuyển Nếu không nắm bắt công nghệ kiểm soát qúa trình vận chuyển Diễn biến qúa trình công nghệ thường phụ thuộc vào thông số công nghệ, cần phải tính toán chi tiết cho thông số công nghệ Việc tính toán thông số công nghệ dựa sở tương quan phương trình dòng Các phương trình dòng thiết lập sở định luật bảo toàn bảo toàn khối lượng, động lượng lượng Số liệu đầu vào phục vụ cho qúa trình tính toán lấy từ tuyến ống dẫn khí từ mỏ Bạch hổ bờ đường kính 406.4 mm dài 110 km Tuyến ống vận chuyển trung bình khoảng 5.5 triệu m3 khí đồng hành/ngày đêm Các số liệu nhiệt độ, áp suất, lưu lượng thành phần khí dòng lưu chất nhận cách đo trực tiếp phân tích mẫu Việc tính toán chênh lệch áp lực dọc tuyến lựa chọn số liệu áp suất hai đầu tuyến ống đo dễ dàng xác Bởi vậy, việc so sánh kết qủa tính thực tế tin cậy Trong thực tế việc tính toán xác thông số công nghệ từ phương trình dòng thường khó thông số liên quan tới nhiều qúa trình hóa lý diễn biến động học phức tạp Bởi vậy, tính toán thông số này, vài giả thiết sử dụng Để tăng mức độ xác tính toán công nghệ, loạt phương trình thực nghiệm có nguồn gốc từ phương trình HV: Vũ Tú Hoài 'String','Giam cao do, m ','BackgroundColor',panelColor, 'horizontalalignment','left','Parent',leftPanel); pipehdv = uicontrol(f,'Style','edit','Units','characters','Position',[43 25 15 1.4], 'String','0','BackgroundColor','white','Parent',leftPanel); opt11 = uicontrol(f,'Style','text','Units','characters','Position',[3 19.5 27 1.4], 'String','Chon TP/TT','foregroundcolor',[0 1],'BackgroundColor',panelColor, 'fontweight','bold','Parent',leftPanel); opt12 = uicontrol(f,'Style','popup','Units','characters','Position',[30 19.5 28 1.4], 'String','Thanh phan|Ty ','BackgroundColor','white','Parent',leftPanel); gas1 = uicontrol(f,'Style','text','Units','characters','Position',[23 17.5 15 1.4], 'fontweight','bold','String','Khi ','BackgroundColor',panelColor,'Parent',leftPanel); liquid1 = uicontrol(f,'Style','text','Units','characters','Position',[43 17.5 15 1.4], 'fontweight','bold','String','Condensate ','BackgroundColor',panelColor,'Parent',leftPanel); gra = uicontrol(f,'Style','text','Units','characters','Position',[3 16 19 1.4], 'String','Ty rieng','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',leftPanel); grag = uicontrol(f,'Style','edit','Units','characters','Position',[23 16 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',leftPanel); gral = uicontrol(f,'Style','edit','Units','characters','Position',[43 16 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',leftPanel); des1 = uicontrol(f,'Style','text','Units','characters','Position',[1 14 13 1.4], 'fontweight','bold','String','Cau tu','BackgroundColor',panelColor,'Parent',leftPanel); frac1 = uicontrol(f,'Style','text','Units','characters','Position',[15 14 13 1.4], 'fontweight','bold','String','Phan mole ','BackgroundColor',panelColor,'Parent',leftPanel); des2 = uicontrol(f,'Style','text','Units','characters','Position',[30 14 13 1.4], 'fontweight','bold','String','Cau tu','BackgroundColor',panelColor,'Parent',leftPanel); frac2 = uicontrol(f,'Style','text','Units','characters','Position',[45 14 13 1.4], 'fontweight','bold','String','Phan mole','BackgroundColor',panelColor,'Parent',leftPanel); ch4 = uicontrol(f,'Style','text','Units','characters','Position',[1 12.5 13 1.4], 'String','CH4','BackgroundColor',panelColor,'Parent',leftPanel); ch4g = uicontrol(f,'Style','edit','Units','characters','Position',[15 12.5 13 1.4], 'String','0.73343','BackgroundColor',panelColor,'Parent',leftPanel); c2h6 = uicontrol(f,'Style','text','Units','characters','Position',[1 11 13 1.4], 'String','C2H6','BackgroundColor',panelColor,'Parent',leftPanel); c2h6g = uicontrol(f,'Style','edit','Units','characters','Position',[15 11 13 1.4], 'String','0.1335','BackgroundColor',panelColor,'Parent',leftPanel); c3h8 = uicontrol(f,'Style','text','Units','characters','Position',[1 9.5 13 1.4], 'String','C3H8','BackgroundColor',panelColor,'Parent',leftPanel); c3h8g = uicontrol(f,'Style','edit','Units','characters','Position',[15 9.5 13 1.4], 'String','0.07258','BackgroundColor',panelColor,'Parent',leftPanel); nc4h10 = uicontrol(f,'Style','text','Units','characters','Position',[1 13 1.4], 'String','n-C4H10','BackgroundColor',panelColor,'Parent',leftPanel); nc4h10g = uicontrol(f,'Style','edit','Units','characters','Position',[15 13 1.4], 'String',' 0.02258','BackgroundColor',panelColor,'Parent',leftPanel); ic4h10 = uicontrol(f,'Style','text','Units','characters','Position',[1 6.5 13 1.4], 'String','i-C4H10','BackgroundColor',panelColor,'Parent',leftPanel); ic4h10g = uicontrol(f,'Style','edit','Units','characters','Position',[15 6.5 13 1.4], 'String','0.01564','BackgroundColor',panelColor,'Parent',leftPanel); nc5h12 = uicontrol(f,'Style','text','Units','characters','Position',[1 13 1.4], 'String','n-C5H12','BackgroundColor',panelColor,'Parent',leftPanel); nc5h12g = uicontrol(f,'Style','edit','Units','characters','Position',[15 13 1.4], 'String','0.00633','BackgroundColor',panelColor,'Parent',leftPanel); ic5h12 = uicontrol(f,'Style','text','Units','characters','Position',[1 3.5 13 1.4], 'String','i-C5H12','BackgroundColor',panelColor,'Parent',leftPanel); ic5h12g = uicontrol(f,'Style','edit','Units','characters','Position',[15 3.5 13 1.4], 'String','0.00565 ','BackgroundColor',panelColor,'Parent',leftPanel); nc6h14 = uicontrol(f,'Style','text','Units','characters','Position',[1 13 1.4], 'String','n-C6H14','BackgroundColor',panelColor,'Parent',leftPanel); nc6h14g = uicontrol(f,'Style','edit','Units','characters','Position',[15 13 1.4], 'String','0.00453','BackgroundColor',panelColor,'Parent',leftPanel); nc7h16 = uicontrol(f,'Style','text','Units','characters','Position',[1 0.5 13 1.4], 'String','n-C7H16','BackgroundColor',panelColor,'Parent',leftPanel); nc7h16g = uicontrol(f,'Style','edit','Units','characters','Position',[15 0.5 13 1.4], 'String','0.00213','BackgroundColor',panelColor,'Parent',leftPanel); nc8h18 = uicontrol(f,'Style','text','Units','characters','Position',[30 12.5 13 1.4], 'String','n-C8H18','BackgroundColor',panelColor,'Parent',leftPanel); nc8h18g = uicontrol(f,'Style','edit','Units','characters','Position',[45 12.5 13 1.4], 'String','0.0095','BackgroundColor',panelColor,'Parent',leftPanel); nc9h20 = uicontrol(f,'Style','text','Units','characters','Position',[30 11 13 1.4], 'String','n-C9H20','BackgroundColor',panelColor,'Parent',leftPanel); nc9h20g = uicontrol(f,'Style','edit','Units','characters','Position',[45 11 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); nc10h22 = uicontrol(f,'Style','text','Units','characters','Position',[30 9.5 13 1.4], 'String','n-C10H22','BackgroundColor',panelColor,'Parent',leftPanel); nc10h22g = uicontrol(f,'Style','edit','Units','characters','Position',[45 9.5 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); n2 = uicontrol(f,'Style','text','Units','characters','Position',[30 13 1.4], 'String','N2','BackgroundColor',panelColor,'Parent',leftPanel); n2g = uicontrol(f,'Style','edit','Units','characters','Position',[45 13 1.4], 'String','0.00181','BackgroundColor',panelColor,'Parent',leftPanel); co2 = uicontrol(f,'Style','text','Units','characters','Position',[30 6.5 13 1.4], 'String','CO2','BackgroundColor',panelColor,'Parent',leftPanel); co2g = uicontrol(f,'Style','edit','Units','characters','Position',[45 6.5 13 1.4], 'String','0.00086','BackgroundColor',panelColor,'Parent',leftPanel); h2s = uicontrol(f,'Style','text','Units','characters','Position',[30 13 1.4], 'String','H2S','BackgroundColor',panelColor,'Parent',leftPanel); h2sg = uicontrol(f,'Style','edit','Units','characters','Position',[45 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); o2 = uicontrol(f,'Style','text','Units','characters','Position',[30 3.5 13 1.4], 'String','O2','BackgroundColor',panelColor,'Parent',leftPanel); o2g = uicontrol(f,'Style','edit','Units','characters','Position',[45 3.5 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); h2 = uicontrol(f,'Style','text','Units','characters','Position',[30 13 1.4], 'String','H2','BackgroundColor',panelColor,'Parent',leftPanel); h2g = uicontrol(f,'Style','edit','Units','characters','Position',[45 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); h2o = uicontrol(f,'Style','text','Units','characters','Position',[30 0.5 13 1.4], 'String','H2O','BackgroundColor',panelColor,'Parent',leftPanel); h2og = uicontrol(f,'Style','edit','Units','characters','Position',[45 0.5 13 1.4], 'String','0','BackgroundColor',panelColor,'Parent',leftPanel); output = uicontrol(f,'Style','text','Units','characters','Position',[20 14 20 1.5], 'fontsize',11,'foregroundcolor',[0 1],'fontweight','bold', 'String','KET QUA','BackgroundColor',panelColor,'Parent',rightupPanel); pres1 = uicontrol(f,'Style','text','Units','characters','Position',[3 12 40 1.4], 'String','Ap suat dau P2, bar','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); pres2 = uicontrol(f,'Style','text','Units','characters','Position',[3 10.5 40 1.4], 'String','Ap suat dau vao P1, bar','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); rate1 = uicontrol(f,'Style','text','Units','characters','Position',[3 40 1.4], 'String','Luu luong, m3/ngd','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); compZ = uicontrol(f,'Style','text','Units','characters','Position',[3 7.5 40 1.4], 'String','He so nen - Z','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); visco = uicontrol(f,'Style','text','Units','characters','Position',[3 40 1.4], 'String','Do nhot , cp','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); fric = uicontrol(f,'Style','text','Units','characters','Position',[3 4.5 40 1.4], 'String','He so ma sat','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightupPanel); inputlable2 = uicontrol(f,'Style','pushbutton','Units','characters','Position',[22.5 4.3 15 2], 'fontsize',13,'fontweight','bold','String','CHAY','foregroundcolor',[1 1], 'BackgroundColor','white','Parent',centerPanel,'callback',@inputlable2callback); center1 = uicontrol(f,'Style','text','Units','characters','Position',[3 2.5 25 1.4], 'fontweight','bold','String','Chon','foregroundcolor',[0 1], 'BackgroundColor',panelColor,'Parent',centerPanel); center2 = uicontrol(f,'Style','text','Units','characters','Position',[32 2.5 25 1.4], 'fontweight','bold','String','Tuong quan','foregroundcolor',[0 1], 'BackgroundColor',panelColor,'Parent',centerPanel); center3 = uicontrol(f,'Style','popup','Units','characters','Position',[3 25 1.4], 'String','Tinh Ap Suat P2|Tinh Ap Suat P1','foregroundcolor',[0 1], 'BackgroundColor','white','Parent',centerPanel); center4 = uicontrol(f,'Style','popup','Units','characters','Position',[32 25 1.4], 'String','Begg & Brin|Flanigane|PanhandleA|PanhandleB|Clinedinst','foregroundcolor',[0 1], 'BackgroundColor','white','Parent',centerPanel); tem11 = uicontrol(f,'Style','text','Units','characters','Position',[3 10 43 1.4], 'String','Nhiet dau vao, oC','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); tem11v = uicontrol(f,'Style','edit','Units','characters','Position',[43 10 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); tem12 = uicontrol(f,'Style','text','Units','characters','Position',[3 8.5 43 1.4], 'String','Nhiet dau ra, oC','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); tem12v = uicontrol(f,'Style','edit','Units','characters','Position',[43 8.5 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); pres11 = uicontrol(f,'Style','text','Units','characters','Position',[3 43 1.4], 'String','Ap suat dau vao P1, bar','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); pres11v = uicontrol(f,'Style','edit','Units','characters','Position',[43 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); pres12 = uicontrol(f,'Style','text','Units','characters','Position',[3 5.5 43 1.4], 'String','Ap suat dau P2, bar','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); pres12v = uicontrol(f,'Style','edit','Units','characters','Position',[43 5.5 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); rate11 = uicontrol(f,'Style','text','Units','characters','Position',[3 43 1.4], 'String','Luu luong khi, m3/ngd','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); rate11v = uicontrol(f,'Style','edit','Units','characters','Position',[43 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); rate12 = uicontrol(f,'Style','text','Units','characters','Position',[3 2.5 43 1.4], 'String','Luu luong long, STB/day','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); rate12v = uicontrol(f,'Style','edit','Units','characters','Position',[43 2.5 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); comZ = uicontrol(f,'Style','text','Units','characters','Position',[3 43 1.4], 'String','Mixing Compresibility - Z','horizontalalignment','left', 'BackgroundColor',panelColor,'Parent',rightdownPanel); comZv= uicontrol(f,'Style','edit','Units','characters','Position',[43 15 1.4], 'String',' ','BackgroundColor',panelColor,'Parent',rightdownPanel); end % inter2 • Tính nhiệt độ áp suất tới hạn function cri_pre_tem % This subrountine finds the critical pressure and critical % temperature of gas from its composition or gravity global PC TC GGRA y OPTION opt= OPTION; if opt==2 GGRA; y(14); y(15); else %Physical constants for typical natural gas constituents %(from Edmister and Lee(1984)) MW=[16.043 30.07 44.097 58.124 58.124 72.151 72.151 86.178 100.205 114.232 128.259 142.286 28.013 44.01 34.076 31.999 2.016 18.015]; PC=[667.8 707.8 616.3 550.7 529.1 488.6 490.4 436.9 396.8 360.6 332.0 304.0 493 1070.9 1306 737.1 188.2 3203.6]; TC = [343.1 549.8 665.7 765.4 734.7 845.4 828.8 913.4 972.5 1023.9 1070.4 1111.8 227.3 547.6 672.4 278.6 59.9 1165.1]; end if (opt==1) GGRA = sum(y.*MW)/28.97; Ppc = sum (y.*PC); Ttc = sum (y.*TC); else %(option==2) Ppc=709.604-58.718*GGRA; Ttc=170.491+307.344*GGRA; end % apply Wichert and Aziz correction for sour gasses A = y(14)+y(15); B = y(15); eps=120*(A^0.9 - A^1.6) + 15*(B^0.5 - B^4); Ttc=Ttc-eps; Ppc=Ppc*Ttc/(Ttc+eps+B*eps*(1-B)); PC=Ppc; TC=Ttc; • Tính hệ số nén Z function compz(P,T) % Calculate gas compressibility Z and reduced gas compressibility CR % in the range 0.2