Xây dựng chương trình tính toán dòng chảy dầu khí

NGHIÊN CỨU KHOA HỌC Thăm dò, Khai thác Xây dựng chương trình tính tốn dịng chảy dầu - khí lên, tựa dừng bề dày lắng đọng paraffin ống đứng Bài báo trình bày tóm tắt lý thuyết dịng chảy hai pha lỏng - khí ống đứng gồm phương trình bảo tồn khối lượng, động lượng, dòng nhiệt tựa dừng kết tinh, lắng đọng paraffin thành ống Hệ phương trình thu nhận được rời rạc hóa giải số sử dụng phương pháp sai phân Runge-Kutta bậc Chương trình tính tốn cho hệ phương trình rời rạc hóa xây dựng ngơn ngữ lập trình Fortran Các tính tốn sử dụng chương trình thực cho mơ hình giếng khai thác dầu khí Các kết tính tốn tham số dịng chảy bề dày lắng đọng paraffin thu nhận phù hợp tốt mặt định tính với thực tế giếng khai thác dầu khí so sánh với kết tính tốn nhóm nghiên cứu khác công bố quốc tế So sánh cho thấy kết thu nhận sử dụng chương trình phát triển nghiên cứu phù hợp tốt với kết giới công bố Chương trình tính tốn có khả áp dụng tốt vào thực tế có triển vọng tiếp tục hồn thiện, nâng cấp tính tốn cho trường hợp tổng qt với giếng khai thác dầu khí điều kiện không dừng hay với hệ thống đường ống vận chuyển dầu khí nằm ngang Đồng thời cần ý rằng, điều kiện bên cố định, tượng không dừng, xảy giếng khai thác có kiện khởi động hay chuyển đổi chế độ làm việc, tồn khoảng thời gian tương đối ngắn Do chế độ hoạt động thủy động lực học nhiệt giếng khai thác dầu - khí xem chế độ tựa dừng, người ta giả thiết nhiệt độ điểm dọc theo trục thân giếng giống cho pha (ở pha lỏng pha khí) Việc nghiên cứu q trình truyền nhiệt trao đổi khối lượng tựa dừng giếng khai thác cho phép đánh giá ảnh hưởng hiệu ứng nhiệt thủy động lực học khác tới tham số quan trọng nhiệt độ, áp suất điểm phân bố dọc theo trục thân giếng Bằng cách đó, việc phân tích nghiên cứu q trình tựa dừng giếng khai thác khơng có ý nghĩa ứng dụng mà đồng thời cịn có ý nghĩa mặt lý thuyết I Mở đầu Trên giới nhiều mỏ dầu - khí đặc trưng nồng độ cao hydrocacbon nặng (hay gọi paraffin) Chuyển động loại dầu thơ có tính chất từ giếng dầu lên đến hệ thống chứa trở nên phức tạp thường có hình thành lắng đọng paraffin ống dẫn Các kết nghiên cứu lắng đọng paraffin giếng dầu tượng có liên quan trình hóa - lý tự nhiên phức tạp xảy dầu thơ q trình đẩy lên mặt đất Về nhiều khía cạnh, đắn phương pháp phương pháp phát triển nghiên cứu lắng đọng paraffin phụ thuộc vào nghiên cứu chi tiết trình mát paraffin Quá trình liên quan chặt chẽ đến trình khác xảy bên xung quanh giếng dầu Một nghiên cứu có liên đới rộng thực thơng qua mơ hình tốn học liên đới chặt chẽ thể dịng chảy khí - lỏng giếng, trao đổi nhiệt giếng với môi trường đá xung quanh, lắng đọng pha rắn (paraffin) thành bên ống dẫn dầu lên (hay gọi ống hay cần nâng) ([1], [2] - [5]) Tìm hiểu nghiên cứu cơng bố dịng chảy nhiều pha cho thấy chuyển động dòng chảy hỗn hợp hai pha lỏng - khí vấn đề phức tạp tiếp tục nghiên cứu Trong chuyển động dịng chảy hai pha lỏng - khí có xét đến q trình kết tinh, lắng đọng paraffin trao đổi nhiệt với môi trường bên ngồi nói có tính phức tạp đặc biệt Tính đa dạng cấu trúc dịng chảy, biến đổi pha thành phần, tượng phức tạp thủy động lực học… làm cho mơ hình tốn trở nên cồng kềnh với nhiều tham số ràng buộc lẫn làm cho tốn khó giải Trong q trình chuyển động dịng chảy hỗn hợp hai pha lỏng - khí dọc theo ống (trong công nghệ khai thác vận chuyển dầu khí, ống khai thác ống vận chuyển) áp suất nhiệt độ hỗn hợp bị mát, hao tán giảm dần dẫn đến hình thành pha rắn (paraffin) mặt ống pha khí q trình tách khí bay Sự hình thành biến đổi lượng paraffin bám thành ống với trình tách khí làm cho đặc trưng cấu trúc dòng chảy thay đổi theo Nghiên cứu chuyển động dịng chảy hai pha lỏng - khí tốn phổ biến nhiều lĩnh vực có nhiều ứng dụng quan trọng, đặc biệt công nghiệp khai thác vận chuyển dầu khí Bài báo trình bày tóm tắt hệ phương trình dịng chảy hai pha lỏng - khí ống đứng có xét đến kết tinh lắng đọng paraffin thành ống Một số kết mơ phỏng, tính tốn, so sánh phân tích bước đầu thực Trên sở định hướng cho nghiên cứu ứng dụng đề Nội dung báo gồm phần sau: Phần I: Mở đầu; Phần II: Mơ hình tốn học dịng hai pha giếng dầu - khí có xét đến kết tinh lắng đọng paraffin; Phần III: Phương pháp số, xây dựng chương trình số kết tính tốn so sánh; Phần IV: Kết luận Phần cuối danh mục Tài liệu tham khảo Danh mục ký hiệu II Mơ hình tốn học dịng hai pha giếng dầu - khí có xét đến kết tinh lắng đọng paraffin Các trình (cơ học, hóa học, vật lý ) xảy dịng chảy dầu - khí giếng khai thác dầu hoạt động phức tạp, có nhiều dạng khác biến đổi liên tục Vì để nghiên cứu q trình điều người ta phải phát triển phương trình cân phụ thuộc vào nhiều giả thiết đơn giản hóa khác cho trường hợp trình riêng biệt xét đến Các phương trình cuối cùng, sử dụng để xây dựng mơ hình tốn học tương ứng, cần phải xác định hay định nghĩa theo đặc trưng trình nghiên cứu [1], [7] Trong q trình chuyển động dịng lỏng - khí giếng khai thác, kích thước đặc trưng hạt hay phần tử (như bọt khí, hạt rắn chất lưu hay dầu) thường nhỏ nhiều lần kích thước đường kính ống nâng kích thước vĩ mơ khác dịng chảy nghiên cứu Vì thế, thể tích vĩ mơ sở hay cịn gọi control volume lựa chọn đủ lớn để chứa phần tử hai pha Điều cho phép người ta mơ tả hỗn hợp lỏng - khí tổng thể hai mơi trường liên tục lấp đầy thể tích sở nghiên cứu Tại điểm thể tích này, ta sử dụng đại lượng vĩ mơ pha như: vận tốc pha, áp suất pha, mật độ khối lượng pha, tỷ phần thể tích khí, áp dụng phương pháp thơng thường tính tốn vi phân, tích phân cho mơ tả toán học hệ hỗn hợp Đồng thời cần ý rằng, điều kiện bên cố định, tượng không dừng, xảy giếng khai thác có kiện khởi động hay chuyển đổi chế độ làm việc, tồn khoảng thời gian tương đối ngắn Do chế độ hoạt động thủy động lực học nhiệt giếng khai thác dầu - khí xem chế độ tựa dừng, người ta giả thiết nhiệt độ điểm dọc theo trục thân giếng giống cho pha (ở pha lỏng pha khí) Việc nghiên cứu q trình truyền nhiệt trao đổi khối lượng tựa dừng giếng khai thác cho phép đánh giá ảnh hưởng hiệu ứng nhiệt thủy động lực học khác tới tham số quan trọng nhiệt độ, áp suất điểm phân bố dọc theo trục thân giếng Bằng cách đó, việc phân tích nghiên cứu q trình tựa dừng giếng khai thác khơng có ý nghĩa ứng dụng mà đồng thời cịn có ý nghĩa mặt lý thuyết Động học kết tinh lắng đọng paraffin, phương trình tốn học mơ tả q trình lắng đọng paraffin Sự hoạt động vận hành giếng dầu - khí hay giếng khí hóa lỏng thường phức tạp xảy hai loại kết tinh Loại thứ dạng khí ngậm nước (hydrat) Điều kiện cần thiết cho hình thành kết tinh loại có mặt nước (hoặc nước) thành phần hydrocacbon mà chúng hòa tan vào nước đến nhiệt độ áp suất hình thành pha rắn Loại khơng xét đến nghiên cứu ống chứa hai thành phần dầu - khí Loại kết tinh thứ hai, thường quan sát thấy trình vận hành hoạt động giếng dầu - khí, biểu hệ hydrocacbon nặng, mà chúng gọi lắng đọng hay kết tinh paraffin Đối với hình thành kết tinh paraffin, hỗn hợp hai pha phải chứa đựng thành phần nặng (yếu tố trực tiếp tạo kết tinh) Ngoài nhiệt độ giếng (đặc biệt nhiệt độ gần thành giếng) phải khoảng cho phép hình thành kết tinh pha rắn Cụ thể nhiệt độ phải nhỏ hay nhiệt độ kết tinh thành phần nặng Ngoài nghiên cứu thành phần hóa học cho thấy, paraffin kết tinh cơng nghiệp khai thác dầu khí hỗn hợp nhiều thành phần phức tạp mặt hóa học bao gồm chất khống thành phần hữu Về trạng thái gồm pha lỏng pha rắn Ảnh hưởng của pha khí lên q trình lắng đọng paraffin theo kết nghiên cứu phương pháp lấy mẫu paraffin khác Cấu thành chủ yếu lắng đọng paraffin khai thác cơng nghiệp thành phần rắn có chứa nhiều hydrocacbon (paraffin asphanten), ngồi có thành phần lỏng có chứa hydrocacbon (nhựa dầu nặng) Thành phần kết tinh bao gồm paraffin khoảng 10 - 75%, asphanten xấp xỉ - %, nhựa 11 - 30% dầu nặng 60% Thành phần chủ yếu chất khoáng bao gồm: hỗn hợp học chiếm tỷ lệ - 5%, nước chiếm vài phần trăm, lưu huỳnh nhỏ 2% (H2S hợp chất chứa lưu huỳnh) Phụ thuộc vào lượng paraffin hòa tan, loại dầu thơ phân loại hình thức thành nhóm gồm: dầu nhiều paraffin (thành phần theo khối lượng lớn 2%), dầu paraffin (thành phần theo khối lượng từ - %), dầu khơng có paraffin (thành phần theo khối lượng nhỏ 1%) [2] - [5], [7] Dầu thô nhiều mỏ dầu có chứa hàm lượng cao thành phần nặng (các cao phân tử hydrocacbon paraffin) Chuyển động loại dầu thô hệ thống ống từ vỉa đến hệ thống thu gom bị phức tạp hóa q trình hình thành lắng đọng pha rắn cao phân tử thành ống Paraffin lắng bám thành giếng dẫn tới làm giảm lưu lượng giếng, làm tăng chi phí cho thiết bị khai thác, vận chuyển Loại bỏ lắng đọng paraffin thân giếng địi hỏi nhiều cơng sức chi phí [2] - [5] Để đặt sở cho phương pháp sử dụng, hay để đề xuất, đưa phương pháp nhằm dự báo giảm thiểu lắng đọng paraffin, nhiều vấn đề thực tế quan trọng khác nảy sinh công nghiệp khai thác dầu khí địi hỏi sử dụng phương pháp vật lý - tốn để mơ tả q trình xảy vỉa thân giếng Phân tích tài liệu cơng bố hình thành paraffin thân giếng xuất có tối thiểu yếu tố sau: + Trong hỗn hợp dầu thơ chảy thân giếng có chứa cao phân tử hydrocacbon, + Trong thân giếng xuất vùng nhiệt độ gần thành giếng nhỏ nhiệt độ bắt đầu kết tinh paraffin Những nghiên cứu Tronov V.P (trích dẫn [1]) kết tinh chế hình thành paraffin rắn Điều có nghĩa điều kiện dịng chảy lỏng - khí có chứa thành phần nặng, hình thành phát triển lắng đọng paraffin xảy xuất phát triển tinh thể trực tiếp bám dính bề mặt thiết bị đất Sự hình thành paraffin vị trí mà dầu tiếp xúc trực tiếp với thành lạnh thiết bị đất Kết giảm nhiệt độ địa phương dịng lỏng - khí lớp biên gần thành thiết bị giảm khả hòa tan chất lỏng paraffin dẫn đến tách pha rắn lắng đọng bề mặt thiết bị Đối với nồng độ paraffin hòa tan hỗn hợp dầu - khí tồn nhiệt độ Te mà thấp nhiệt độ paraffin rắn hình thành lắng đọng, hay nói cách khác mỏ dầu lắng đọng paraffin xuất tất chế độ làm việc giếng nhiệt độ miệng giếng thấp nhiệt độ bắt đầu kết tinh paraffin Cường độ lắng đọng paraffin xác định qua hoạt động yếu tố khác nhau: làm lạnh thành ống xuống nhiệt độ bắt đầu kết tinh paraffin, nhiệt thủy động lực học hỗn hợp lỏng - khí, yếu tố khí, đặc tính bề mặt thiết bị… Tuy kết tinh trực tiếp bề mặt cứng thiết bị chế hình thành lắng đọng paraffin Điều khẳng định qua quan sát trực tiếp, qua nghiên cứu ảnh hưởng nhiệt độ dịng chảy lỏng - khí thành ống tới cường độ lắng đọng paraffin thiết bị khai thác dầu [1] - [8] Theo mơ hình kết tinh paraffin, trường hợp tổng quát, thu nhỏ thiết diện thân giếng lắng đọng paraffin xuất với hình thành hạt kết tinh pha rắn lòng chất lỏng lớp biên chảy sát thành ống bị làm lạnh xuống nhiệt độ kết tinh thành phần hydrocacbon nặng, hạt rắn lắng đọng kết bám vào thành thiết bị Ngồi ra, việc hình thành phát triển pha rắn thành giếng, trình khuếch tán nhiệt, thành phần hydrocacbon nặng đóng vai trị ảnh hưởng Tuy nhiên mơ hình hóa định tính tượng lắng đọng paraffin lên thành giếng xét đến tất yếu tố việc làm phức tạp, kể đến việc mơ hình hóa chi tiết địi hỏi phải đưa vào mơ hình lượng lớn tham số thực nghiệm Vì ta giả thiết tốc độ phát triển trình lắng đọng paraffin thành cần nâng xác định chủ yếu vào nồng độ thành phần hydrocacbon nặng, vào nhiệt độ chế trao đổi nhiệt - chất dòng bề mặt pha rắn Động học trình kết tinh paraffin mơ tả sử dụng mơ hình khuếch tán, chế độ chảy rối dòng chảy hỗn hợp lỏng - khí, kích cỡ xung động bậc với đường kính ống lớn nhiều so với kích thước khu vực xảy kết tinh lắng đọng paraffin Thêm vào đó, xem xét q trình phát triển hòa tan lắng đọng paraffin ta quy ước bỏ qua hiệu ứng liên kết phân tách [1] Giả thiết hình thành pha rắn xác định chủ yếu giảm tốc độ dòng chảy chất lỏng khu vực gần thành ống chất lỏng bị nhiệt giảm nhiệt độ Giả thiết có nghĩa chất lỏng có chứa thành phần nặng nhiệt độ giảm hình thành pha rắn Sự khuếch tán thành phần khuếch tán nhiệt mang thành phần hydrocacbon nặng đến gần thành giếng tạo điều kiện cho hình thành phát triển pha rắn kết tinh bám dính vào thành giếng Cường độ lắng đọng paraffin phụ thuộc vào vận tốc dịng chảy giếng vào lưu lượng đầu giếng Tốc độ trình lắng đọng paraffin xác định điều kiện cân nhiệt bề mặt lắng đọng pha rắn mối quan hệ tương tác qua lại dòng chảy lỏng - khí paraffin Các tính tốn thực kết so sánh với kết nhóm nghiên cứu Shagapov người khác [1] Các tính tốn sử dụng chương trình phát triển nghiên cứu thực với trường hợp nồng độ pha nặng tương ứng 6,83; 5,0 1,43% Nồng độ pha nặng nghiên cứu Shagapov người khác (và số tham số khác giếng) không đưa giá trị cụ thể so sánh thực cho ba trường hợp kết cho thấy trường hợp nồng độ pha nặng 5% cho kết tính tốn gần với kết nhóm Shagapov Các Hình 3, 4, trình bày kết so sánh cho phân bố nhiệt độ dọc theo trục giếng bề dày lớp paraffin lắng đọng thành ống Kết tính tốn áp suất không đưa nghiên cứu Shagapov [1] nên tính tốn khơng có khơng trình bày Nhận xét: Kết so sánh Hình cho thấy phân bố nhiệt độ dọc theo trục giếng có phù hợp tốt định tính định lượng Cần lưu ý kết tính tốn nhiệt độ nhóm Shagapov cho trường hợp tính tốn thủy động lực học giếng chưa xét đến kết tinh paraffin (kết trình bày tham số giếng cơng bố) Do kết tính tốn khơng thể hồn tồn trùng khớp Đặc biệt kể từ vị trí có kết tinh paraffin kết so sánh nhiệt độ cho thấy khác biệt rõ Khi có kết tinh paraffin ảnh hưởng trường địa nhiệt lên dòng chảy ống thay đổi lớn kết cho thấy phần nhiệt độ giảm chậm nhiều so với khơng có kết tinh paraffin (Shagapov [1]) Kết so sánh bề dày lắng đọng paraffin Hình cho thấy số khác biệt Như cần lưu ý nồng độ pha rắn số tham số khác giếng lớp vật chất khe vành khuyên, hệ số dẫn nhiệt chúng v.v nghiên cứu Shagapov không đưa số liệu cụ thể nên có ảnh hưởng đến kết so sánh Tuy nhiên trường hợp thấy Hình 4, sai số tối đa khoảng 0,5mm (với bề dày tối đa lớp paraffin khoảng 2,0mm), sai số tương đối nhỏ mức chấp nhận Trường hợp tính toán với nồng độ thành phần nặng pha lỏng 6,83%, kết tính tốn áp suất, nhiệt độ kết tinh paraffin dọc giếng theo thời gian trình bày Hình 6, đây: Nhận xét: Tính tốn cho thấy nhiệt độ thành giếng giảm xuống tới nhiệt độ kết tinh paraffin (3030K), giếng xuất kết tinh thành phần nặng pha lỏng (bắt đầu độ cao 2270m) Theo thời gian hoạt động giếng, kết tinh paraffin thành cần nâng xảy ngày nhanh và lấp đầy khoảng 2/3 bán kính miệng giếng sau thời gian hoạt động giếng khoảng 204 (8,5 ngày) Từ vị trí xảy kết tinh paraffin trao đổi nhiệt bề mặt cần nâng với môi trường xảy chậm, không đáng kể (đó nhiệt độ bão hịa paraffin) Để tìm hiểu ảnh hưởng thành phần nặng tới tốc độ kết tinh paraffin dọc thân giếng trình hoạt động giếng, nhóm nghiên cứu tiến hành tính tốn cho mơ hình tốn với nồng độ thành phần nặng thấp (1,43%) Kết tính tốn trình bày Hình Nhận xét: Từ kết tính tốn trình bày Hình 8, cho thấy nồng độ thành phần nặng trọng pha lỏng có ảnh hưởng lớn tới q trình ngưng tụ kết tinh paraffin mặt thành cần cần nâng Với nồng độ thành phần nặng 1,43% sau khoảng 10 ngày giếng hoạt động liên tục bề dày lớp kết tinh paraffin đạt 5mm Trong với nồng độ thành phần nặng 6,83% bề dày kết tinh >20mm Nồng độ thành phần nặng pha lỏng có ảnh hưởng lớn tới phân bố áp suất nhiệt độ giếng Nồng độ thành phần nặng cao nhiệt độ áp suất giảm q trình kết tinh paraffin nhanh thấy Hình 10, 11 cho trường hợp tính tốn so sánh với nồng độ thành phần nặng pha lỏng 1,43%; 5% 6,83% IV Kết luận Nghiên cứu thực bao gồm việc tìm hiểu mơ hình dịng chảy hai pha dầu - khí có lắng đọng paraffin thành ống dẫn Hệ phương trình đóng kín mơ tả tốn, phương pháp số áp dụng cho việc rời rạc hóa tốn giải số hệ phương trình trình bày Trên sở chương trình tính tốn tham số dịng chảy bề dày lắng đọng paraffin xây dựng sử dụng ngơn ngữ lập trình Fortran Một số tính tốn nghiên cứu, so sánh phân tích thực Các kết tính tốn so sánh với kết tính tốn công bố khác giới [1, 13] Các kết so sánh cho thấy có phù hợp tốt mặt định tính định lượng Thực tế cho thấy toán thủy động lực học trao đổi nhiệt - chất dòng chảy nhiều pha tốn khó, đặc biệt tốn động học kết tinh thành phần nặng trình chuyển động pha, thành phần tốn phức tạp lại có ứng dụng rộng rãi thực tế Trong khuôn khổ nghiên cứu này, kết tính tốn nhóm nghiên cứu thực cho trường hợp dòng chảy tựa dừng Các nghiên cứu thực theo hướng: + Cải tiến mơ hình tốn học nghiên cứu chế độ hoạt động giếng trường hợp không dừng tổng quát, + Cải tiến mơ hình tốn học để bao gồm giếng hay hệ thống ống đứng hệ thống nằm ngang, + Mở rộng nghiên cứu cho hệ dầu - nước dầu - khí - nước Tương ứng với mơ hình lý thuyết cải tiến, chương trình tính tốn cần xây dựng phục vụ tính tốn so sánh với số liệu thí nghiệm hay số liệu thực tế Điều đặt nhu cầu xây dựng mơ hình thí nghiệm vật lý với kích cỡ thu nhỏ phịng thí nghiệm để nghiên cứu thực nghiệm q trình kết tinh lắng đọng paraffin hệ thống ống đứng hay nằm ngang, đồng thời cung cấp số liệu thí nghiệm cho mơ tính tốn sử dụng chương trình máy tính Tuy nhiên với mơ hình thí nghiệm vật lý thu nhỏ nghiên cứu tính đồng dạng tham số đặc trưng tốn thí nghiệm thực tế vấn đề cần phải giải Các tính tốn sử dụng chương trình máy tính sau so sánh kiểm nghiệm với kết thí nghiệm áp dụng vào nghiên cứu trường hợp mỏ cụ thể thực tế Trong điều kiện Việt Nam nước có nhiều mỏ khai thác dầu khí, điều có ý nghĩa khoa học ứng dụng to lớn Tài liệu tham khảo Shagapov V.Sh., Musakaev N.G., Khabeev N.S and Bailey S.S., 2004 Mathematical modelling of two-phase flow in a vertical well considering paraffin deposits and external heat exchange International Journal of Heat and Mass Transfer, Vol.47, pp.843 - 851 Ribeiro F S., Mendes P R S., Braga S L., 1997 Obstruction of pipelines due to paraffin deposition during the flow of crude oils Int J Heat Mass Transfer, Vol.40, Issue 18, pp.4319-4328 Hamouda A.A and Davidsen S., 1995 An Approach for Simulation of Paraffin Deposition in Pipelines as a Function of Flow Characteristics With a Reference to Teesside Oil Pipeline SPE International Symposium on Oilfield Chemistry, DOI: 10.2118/28966-MS Brown T.S., Niesen V.G and Erickson D.D., 1993 Measurement and Prediction of the Kinetics of Paraffin Deposition Soc Pet Eng., Paper No.26548, pp.353-368 Noman Shahreyar, 2000 Review of paraffin control and removal in oil wells using Southwestern Petroleum Short Course searchable database A Master Thesis in Petroleum Engineering, Texas Tech University, USA Armand A A., 1950 Resistance in the process of flow of a two - phase system through horizontal tube Vsesoyuznogo Teplotekhn Inst., No Chisholm D., 1983 Two - phase flow in pipelines and heat exchangers Longman Higher Education, London Goodman J., 1958 The heat balance integral and its application to problem involving change of phase Trans Soc Mech Engrs., Vol.80, pp.335 - 442 GS.TSKH Dương Ngọc Hải, 2003, Xây dựng chương trình tính máy tính cá nhân chuyển động không đẳng nhiệt hỗn hợp dầu - khí - nước có xét đến kết tinh paraffin Hà Nội 10 Hewitt G and Butterworth D (Eds.), 1977 Two - phase Flow and Heat Transfer Oxford University Press 11 Nigmatulin R.I (1991), Dynamics of Multiphase Media, Hemi-sphere Publ Corp, NewYork 12 Sloan E.D and Selim M.S., 1989 Heat and mass transfer during the dissociation of hydrates in porous media AIChe Journal, Vol.35, pp.1049 - 1052 13 Shagapov V.Sh., Borodin S.L., Gubaidullin A.A., Dương Ngọc Hải and Musakaev N.G., 2009 Mathematical Modelling of Upward Flow of a Liquid - Gas Mixture in a Vertical Well Proceeding of the international Workshop on Thermal Hydrodynamics Of Multiphase Flows And Applications, pp.161172 14 Encyclopedia Britanica, UK (http://www.britannica.com) 15 Wikipedia Encyclopedia (http://www.wikipedia.org) 16 GS.TSKH Dương Ngọc Hải (chủ nhiệm đề tài),1987 Công nghệ khai thác vận chuyển dầu nhiều paraffin độ nhớt cao mỏ Bạch Hổ Trung tâm Khoa học Tự nhiên Công nghệ Quốc gia, Viện Cơ học, Hà Nội, 88 trang 17 GS.TSKH Dương Ngọc Hải (chủ nhiệm đề tài), 1997 Xây dựng phương pháp lập chương trình tính máy tính cá nhân chuyển động khơng đẳng nhiệt hệ lỏng khí ba pha thành phần ống khai thác thẳng đứng nghiêng Trung tâm Khoa học Tự nhiên Công nghệ Quốc gia, Viện Cơ học, Hà Nội, 138 trang 18 Dương Ngọc Hải, Gubaidullin A., Musakaev N.G Đỗ Văn Chung, 2009 Mơ hình hóa, tính tốn dịng chảy hỗn hợp hai pha lỏng - khí ống đứng Tuyển tập cơng trình khoa học Hội nghị Cơ học toàn quốc - Kỉ niệm 30 năm Viện Cơ học Tạp chí Cơ học, NXB Khoa học Tự nhiên Công nghệ, Hà Nội, trang 421 - 429 19 Nguyễn Duy Trọng (chủ nhiệm đề tài), 2010 Nghiên cứu lý thuyết dịng chảy khơng đẳng nhiệt lỏng - khí ống đứng có kết tinh paraffin Báo cáo đề tài nghiên cứu khoa học 2010 - Viện Cơ học, Viện Khoa học Công nghệ Việt Nam Danh mục ký hiệu: Các số l: Pha lỏng g: Pha khí (l): Thành phần trung bình (s): Thành phần nặng f: Tham số đầu vào đáy (g): Thành phần nhẹ giếng Dương Ngọc Hải Nguyễn Duy Trọng Nguyễn Tất Thắng Viện Cơ học, Viện Khoa học Công nghệ Việt Nam Theo TCDK  ... ngữ lập trình Fortran Các tính tốn sử dụng chương trình thực cho mơ hình giếng khai thác dầu khí Các kết tính tốn tham số dịng chảy bề dày lắng đọng paraffin thu nhận phù hợp tốt mặt định tính với... chương trình tính tốn tham số dịng chảy bề dày lắng đọng paraffin xây dựng sử dụng ngôn ngữ lập trình Fortran Một số tính tốn nghiên cứu, so sánh phân tích thực Các kết tính tốn so sánh với kết tính. .. nằm ngang, + Mở rộng nghiên cứu cho hệ dầu - nước dầu - khí - nước Tương ứng với mơ hình lý thuyết cải tiến, chương trình tính tốn cần xây dựng phục vụ tính tốn so sánh với số liệu thí nghiệm