Nghiên cứu mô phỏng và thực nghiệm quá trình khí hóa than ngầm

TÓM TẮT Như biết lượng kỷ XXI vấn đề nóng hổi tồn cầu Khi nguồn nhiên liệu dầu khí đốt dự báo cạn kiệt vòng 50 đến 60 năm tới, dẫn đến giá dầu, khí ngày tăng cao làm cho nhiều ngành sản xuất phụ thuộc nhiều vào nguồn nhiên liệu phải lao đao đặc biệt quốc gia nhập dầu, khí Các nguồn lượng tái tạo : lượng mặt trời, lượng gió, địa nhiệt, lượng biển… năm gần người ta nghiên cứu ứng dụng nhiều, hiệu suất thiết bị thấp, chưa thể đáp ứng nhu cầu sử dụng lượng Trong nguồn nhiên liệu hố thạch than đá với trữ lượng lớn phân bố rộng khắp toàn cầu Để giải vấn đề lượng vài trăm năm tới việc sử dụng than đá giải pháp có ưu Nhưng vấn đề nguồn nhiên liệu sử dụng theo lối truyền thống phát thải lớn điều thời đại ngày Đặc biệt Việt Nam với việc phát mỏ than lòng Đồng Bằng Sông Hồng (ĐBSH) với trữ lượng lớn gần 210 tỉ tấn, vỉa dày - 20m, nằm sâu 120 - 220m, khơng có nước ngầm, nguồn tài nguyên khổng lồ cho ngành lượng ngành cơng nghiệp khác Vì khí hố than đá phương pháp để chuyển than đá thành khí đốt dùng làm ngun liệu tổng hợp hóa chất Phương pháp ứng dụng nhiều năm gần Với ưu điểm lợi ích việc khí hóa than đem lại nên nhóm chọn đề tài "Nghiên cứu mô thực nghiệm q trình khí hố than ngầm" để đóng góp phần nhỏ vào kho tài liệu chuyên nghành khí hóa than góp phần vào cơng đổi nước nhà Em đạt số thành cơng như: - Mơ q trình truyền nhiệt, truyền chất q trình khí hóa than ngầm - Xây dựng thực nghiệm mơ hình khí hóa than thu kết khí hóa mong đợi thu khí cháy, đo nồng độ khí cháy sinh kết khác Tuy cịn gặp nhiều khó khăn nhóm xây dựng phương pháp khí hóa than dù chưa hoàn hảo MỤC LỤC LỜI CÁM ƠN………… ……………………………………………… TÓM TẮT .2 PHỤ LỤC HÌNH ẢNH PHỤ LỤC BẢNG .10 BẢNG THUẬT NGỮ 11 MỞ ĐẦU 14 Tính cấp thiết đề tài 14 Đối tượng phạm vi nghiên cứu 14 Ý nghĩa lí luận thực tiễn đề tài 18 Mục tiêu nghiên cứu 19 Nhiệm vụ nghiên cứu 20 Giới hạn 20 Phương pháp nghiên cứu 20 CHƯƠNG 1: TỔNG QUAN 21 1.1 Tình hình lượng Việt Nam 21 1.2 Tổng quan loại than Việt Nam 22 1.3 Tình hình nghiên cứu UCG nước 24 1.4 Các thí nghiệm liên quan đến UCG 30 1.5 Một số quy trình khí hóa khí than áp dụng giới 32 CHƯƠNG 2: CƠ SỞ LÝ THUYẾT VỀ KHÍ HĨA THAN NGẦM 35 2.1 Đặc điểm nguyên lý khí hóa than 35 2.2 Cơ sở lý thuyết 37 2.3 Mơ hình dịng chảy rối k-ε 47 2.4 Xây dựng phương trình giải pháp điều chỉnh 48 2.5 Tính tốn lý thuyết than sử dụng làm thực nghiệm 52 CHƯƠNG 3: PHƯƠNG PHÁP MÔ PHỎNG QUÁ TRÌNH KHÍ HĨA THAN NGẦM BẰNG PHẦN MỀM COMSOL MULTIPHYSICS 55 3.1.Tổng quan phần mềm Comsol multiphysic 55 3.1.1 Tầm quan trọng phần mềm Comsol multiphysics 3.1.2 Cài đặt phần mềm 3.2 Xây dựng mơ hình mơ 58 3.2.1 Thiết kế mơ hình 3.2.2 Vẽ mơ hình Comsol Multiphysics 3.2.3 Cài đặt điều kiện thông số đầu vào cho mô 3.3 Mô trình truyền nhiệt truyền chất 63 3.3.1 Lựa chọn lời giải phương trình 3.3.2 Thiết lập thơng số điều kiện cho phương trình 3.3.3 Tạo lưới giải mơ hình CHƯƠNG 4: Q TRÌNH THỰC NGHIỆM Q TRÌNH KHÍ HĨA THAN NGẦM 76 4.1 Mục đích ý nghĩa 76 4.2 Mô tả hệ thống 76 4.2.1 Ý nghĩa mơ hình khí hóa thực nghiệm 4.2.2 Ngun liệu than để khí hóa 4.3 Thuyết minh dụng cụ đo 79 4.4 Quá trình thực nghiệm 84 CHƯƠNG KẾT QUẢ MÔ PHỎNG VÀ THỰC NGHIỆM 92 5.1 Kết mô 92 5.1.1 Kết mô khoang nhiệt 92 5.1.2 Kết mơ hình khoang rỗng 94 5.1.3 Kết mô vectơ vận tốc 95 5.1.4 Kết mơ q trình truyền chất 97 5.2 Kết thực nghiệm 103 5.2.1 Kết than sau khí hố 103 5.2.2 Kết đo nhiệt độ 104 5.2.3 Kết đo nhiệt độ so sánh với mô 105 5.2.4 Kết đốt cháy sản phẩm khí 106 5.2.5 Kết phân tích khí cháy 107 CHƯƠNG KẾT LUẬN VÀ KIẾN NGHỊ 115 TÀI LIỆU THAM KHẢO 117 Bài báo 119 Bài báo 124 PHỤ LỤC HÌNH ẢNH Hình Ngun lý khí hóa than ngầm 15 Hình Quá trình khí hóa than ngầm 18 Hình Phương pháp CRIP 19 Hình 1.1 Độ dày so với độ sâu thử nghiệm tiến hành toàn giới (C.Pana, 2009) 26 Hình 1.2 Hình chiếu khoang Hanna II (sửa đổi từ năm 1989) [19] 28 Hình 1.3 Mặt cắt ngang khoang thí nghiệm LBK 30 Hình 2.1 Ngun lí khí hố than ngầm 36 Hình 2.2 Giản đồ mơ hình khí hố than 41 Hình 2.3 Chia lưới miền khí hố 49 Hình 2.4 Điều kiện biên khối lượng 51 Hình 3.1 Hình ảnh thư mục Comsol Multiphysics 55 Hình 3.2 Chọn ngôn ngữ Comsol Multiphysics 56 Hình 3.3 Chọn cài đặt Comsol Multiphysics 56 Hình 3.4 Chọn cài đặt Comsol Multiphysics 57 Hình 3.5 Mơ hình mơ khí hố dạng 2D 59 Hình 3.6 Thiết lập mơ hình mơ (1) 59 Hình 3.7 Thiết lập mơ hình mơ (2) 60 Hình 3.8 Thiết lập mơ hình mơ (3) 60 Hình 3.9 Thiết lập mơ hình mơ (4) 61 Hình 3.10 Mơ hình mơ thiết lập 61 Hình 3.11 Thành phần điều kiện than mơ khí hố 62 Hình 3.12 Chọn Comsol Multiphysics Physics 63 Hình 3.13 Chọn Comsol Multiphysics Study 64 Hình 3.14 Thiết lập điều kiện biên cho Reaction Engineering (1) 64 Hình 3.15 Thiết lập điều kiện biên cho Reaction Engineering (2) 65 Hình 3.16 Thiết lập điều kiện biên cho Reaction Engineering (3) 65 Hình 3.17 Thiết lập điều kiện biên cho Reaction Engineering (4) 66 Hình 3.18 Thiết lập điều kiện biên cho Reaction Engineering (5) 66 Hình 3.19 Thiết lập điều kiện biên cho Reaction Engineering (6) 67 Hình 3.20 Thiết lập điều kiện biên cho Reaction Engineering (7) 68 Hình 3.21 Thiết lập điều kiện biên cho Reaction Engineering (8) 68 Hình 3.22 Thiết lập điều kiện biên cho Reaction Engineering (9) 68 Hình 3.23 Thiết lập điều kiện biên cho Reaction Engineering (10) 69 Hình 3.24 Thiết lập điều kiện biên cho Reaction Engineering (11) 69 Hình 3.25 Thiết lập điều kiện biên cho Reaction Engineering (12) 70 Hình 3.26 Thiết lập điều kiện biên cho Reaction Engineering (13) 70 Hình 3.27 Thiết lập điều kiện biên cho Reaction Engineering (14) 71 Hình 3.28 Thiết lập điều kiện biên cho Materials 71 Hình 3.29 Thiết lập điều kiện biên cho Brinkman Equations (1) 72 Hình 3.30 Thiết lập điều kiện biên cho Brinkman Equations (2) 72 Hình 3.31 Thiết lập điều kiện biên cho Heat Tranfer in Porous Media (1) 73 Hình 3.32 Thiết lập điều kiện biên cho Heat Tranfer in Porous Media (2) 73 Hình 3.33 Thiết lập điều kiện biên cho Transport of Diluted in Porous Media (1) 74 Hình 3.34 Thiết lập điều kiện biên cho Transport of Diluted in Porous Media (2) 74 Hình 3.35 Tạo lưới cho mơ hình 75 Hình 3.36 Thiết lập Study cho giải thuật 75 Hình 4.1: Mơ hình thử nghiệm UCG 76 Hình 4.2 Mơ hình thực nghiệm khí hóa than ngầm 77 Hình 4.3 Hình ảnh than đá dùng để khí hóa 79 Hình 4.4 Điện trở sấy khơng khí 79 Hình 4.5 Điện trở sấy khơng khí 80 Hình 4.6 Hình ảnh thực tế can nhiệt 80 Hình 4.7 Bình Oxy 81 Hình 4.8 Thiết bị phân tích khí Testo 350XL 83 Hình 4.9 Camera nhiệt Fluke Ti9 83 Hình 4.10 Q trình chuẩn bị mơ hình thực nghiệm 85 Hình 4.11 Q trình thử kín mơ hình thực nghiệm sử dụng khí Nitơ 87 Hình 4.12 Kiểm tra thơng số tín hiệu điều khiển hộp điện 87 Hình 4.13 Q trình mồi lửa gia nhiệt khí Oxy 88 Hình 5.1 Hình ảnh mơ khoang nhiệt khí hóa t = 30 92 Hình 5.2 Kết mơ q trình truyền nhiệt lúc tâm phản ứng đạt 5000C 92 Hình 5.3 Kết mơ q trình truyền nhiệt lúc tâm phản ứng đạt 7500C 93 Hình 5.4 Hình dạng khoang rỗng hình thành 94 Hình 5.5 Hình chiếu khoang Hanna II (sửa đổi từ năm 1989) 94 Hình 5.6 Kết mơ trình truyền nhiệt lúc tâm phản ứng vecto vận tốc 95 Hình 5.7: Phân bố nhiệt vecto vận tốc khoang khí hóa dùng phần mềm CFD 96 Hình 5.8 Dạng hình học dự đoán khoang rỗng dùng phần mềm CFD 96 Hình 5.9 Thành phần sản phẩm khí thu cho nồng độ O2 0.10 mol/ m3 97 Hình 5.10 Thành phần sản phẩm khí thu cho nồng độ O2 0.15 mol/ m3 97 Hình 5.11 Thành phần sản phẩm khí thu cho nồng độ O2 0.3 mol/ m3 98 Hình 5.12 Thành phần sản phẩm khí thu cho nồng độ O2 0.45 mol/ m3 99 Hình 5.13 Thành phần sản phẩm khí thu cho nồng độ O2 0.6 mol/ m3 99 Hình 5.14 Thành phần sản phẩm khí thu cho nồng độ O2 0.75 mol/ m3 100 Hình 5.15 Thành phần sản phẩm khí thu cho nồng độ O2 mol/ m3 100 Hình 5.16 Nồng độ CO2 khoang khí hóa t = 1800s 101 Hình 5.17 Nồng độ CO khoang khí hóa t = 1800s 102 Hình 5.18 Nồng độ H2 khoang khí hóa t = 1800s 102 Hình 5.19 Than sau khí hóa, xỉ than 103 Hình 5.20 Hình ảnh so sánh than trước sau khí hố 104 Hình 5.21 Hình ảnh khoang nhiệt khí hố thời điểm t = 15 105 Hình 5.22 Hình ảnh khoang nhiệt khí hố t = 30 105 Hình 5.23 So sánh kết khoang nhiệt mô thực nghiệm 106 Hình 5.24 Hình ảnh que đốm bùm cháy trở lại đưa vào miệng ống thu khí 107 Hình 5.25 So sánh q trình truyền chất mơ thực nghiệm 109 Hình 5.26 Kết từ máy phân tích khí Testo 350 - XL 111 Hình tham khảo:Một số hình ảnh kết phân tích khí máy Testo 350XL 111 PHỤ LỤC BẢNG Bảng 1.1: So sánh số nhiên liệu giá thành sản xuất nhiệt 21 Bảng 1.2 Những thử nghiệm quan trọng lĩnh vực UCG……………………… 27 Bảng 1.3 Đại diện cho nghiên cứu phịng thí nghiệm UCG……………… 31 Bảng 2.1 Thành phần hóa học có than đá sử dụng khí hóa 48 Bảng 3.1 Thành phần chất có than thực nghiệm 62 Bảng 3.2 Thành phần đặc tính than mơ 62 Bảng 3.3 Thành phần oxy mô 62 Bảng 4.1.Thành phần chất có than chưa đốt cháy 79 Bảng 5.1 Kết đo nhiệt độ q trình khí hố than theo thời gian 104 Bảng 5.2 Bảng kết phân tích khí sản phẩm khí hóa theo thời điểm 107 Bảng 5.3 Bảng quy đổi nồng độ chất từ mg/m3 sang mol/m3 108 Bảng 5.4 Bảng tính nhiệt trị theo thể tích sản phẩm khí…………………………… 110 10 BẢNG THUẬT NGỮ Danh pháp Cgi ug D Ks Ý nghĩa Nồng độ mol thành phần i Vận tốc dịng khí dọc trục Hệ số khuếch tán dịng khí Thời gian Tỷ lệ sản phẩm thành phần thứ i phản ứng hóa học Nồng độ mol thành phần I phía đầu kênh khí hố Nhiệt độ pha khí Tổng nồng độ mol pha khí Nhiệt dung riêng chất khí thứ I Tốc độ phản ứng hóa học Nhiệt lượng phản ứng hóa học Hệ số trao đổi nhiệt pha rắn khí Diện tích tiếp xúc pha rắn khí Nhiệt lượng tổn thất pha khí Thành phần phần trăm chất khí thứ i pha khí Nhiệt độ biên pha khí Nhiệt độ ban đầu pha khí Hệ số dẫn nhiệt pha khí Nhiệt độ pha rắn Nhiệt dung riêng pha rắn Nhiệt lượng tổn thất pha rắn Hệ số dẫn nhiệt pha rắn s Khối lượng riêng chất rắn kg/m3 Lưu lượng thể thích thành phần thứ i Nhiệt phản ứng hóa học chất thứ i Tổng số chất Áp suất dịng khí Hệ số dính kết Hệ số hiệu chỉnh Nguồn, hội tụ Hằng số chất khí Tỷ lệ rị rỉ chất khí kg/m3 kJ/mol Si Cgi0 Tg Cg Cpgi Ri Hi A Qg ygi Tggz0, Tgr0 Tg0 Kλg Ts Cs Fs mi ∆Hi N Pg   R k u p  Vận tốc Áp suất Khối lượng riêng Đơn vị mol/m3 m/s m2/s s mol/ m3s mol/m3 K mol/m3 kJ/mol.K mol/ m3.s kJ/mol kJ/( m2K) m2/ m3 kJ K K kW/mK K kJ/kg.K kJ/kg.s.K kW/mK Pa Pa.s J/mol.K m m/s Pa kg/m3 11 hoá, mà xem sản phẩm trình cháy thiếu Oxy Tuy nhiên coi kết bước đầu lĩnh vực nghiên cứu hóa khí than ngầm Việt Nam - Yếu tố thời gian: thời gian thực luận văn có giới hạn, trình thực cần nhiều thời gian để nghiên cứu tìm giải pháp thực mơ cách xác - Độ xác thiết bị đo cịn chưa mang lại tính xác cao - Yếu tố kinh phí ảnh hưởng lớn đến chất lượng luận văn lần thực khí hóa tốn nhiều kinh phí mà em chưa đáp ứng Kiến nghị hướng phát triển đề tài: - Từ bước đầu tìm hiểu tài liệu sử dụng phần mềm Comsol 5.2 phần mang lại kết có tính tương đối, để đạt kết cao xác cần thêm nhiều nỗ lực để nghiên cứu phần mềm mô Comsol - Do điều kiện thực nghiệm khó khăn, kiến thức kỹ thuật hạn chế nên mơ hình thực nghiệm chưa thể đưa vào thực tế Để đưa cơng nghệ khí hóa vào thực tế cần giúp đỡ từ chuyên gia, sử dụng công nghệ tiến tiến cần hỗ trợ đầu tư từ nguồn khác - Bên cạnh tiếp tục phát triển, nâng cấp mơ hình khí hóa than có sẵn để tiếp tục tiến hành thí nghiệm, mang lại số liệu thống kê tốt hơn, nâng cao hiệu suất khí hóa 116 TÀI LIỆU THAM KHẢO [1] Đặng Quốc Phú, Trần Thế Sơn, Trần Văn Phú, Truyền nhiệt , NXB giáo dục, 2004 [2] Patankar S.V, Numerical Heat Transfer and Fluid Flow, McGraw Hill, 1980 [3] Phạm Lê Dzần, Nguyễn Công Hân Công nghệ lò mạng nhiệt Nhà xuất Khoa Học Kỹ Thuật Hà Nội – 2008 [4] PGS.TS Bùi Hải, PGS.TS Trần Thế Sơn, Bài Tập Truyền Nhiệt - Nhiệt Động Và Kỷ Thuật Lạnh, Nhà xuất khoa học kĩ thuật Hà Nội, 2001 [5] Nguyễn Bin, Đỗ Văn Đài, Long Thanh Hùng, Đinh Văn Huỳnh, Nguyễn Trọng Khng, Phạm Văn Thơm, Phạm Xn Tồn, Trần Xoa, Sổ Tay Q Trình Và Thiết Bị Cơng Nghệ Hóa Chất – Tập 1, Nhà xuất khoa học kĩ thuật Hà Nội, 1999 [6] Nguyễn Thanh Quang (ĐHBK - Đà Nẵng), Đặng Thế Hùng (Công ty TNHH Trường Quang II), “Nghiên Cứu Chế Tạo Hệ Thống Hóa Khí Than Tầng Cố Định Ngược Chiều”, Tạp Chí Khoa Học Và Công Nghệ Nhiệt Số 77 [7] Gasification Technology, Technical Issues in the Design of Gasifiers, 1999 [8] Krzysztof Stanczyk, Krzysztof Kapusta, Experimental simulation of hard coal underground gasification for hydrogen production, Fuel, 2012 [9] Krzysztof Stanczyk Hydrogen-oriented underground coal gasification for Europe (HUGE)-Euro commission, 2009 [10] Krzysztof Stanczyk, Krzysztof Kapusta, Pollution of water during underground coal gasification of hard coal and lignite, Fuel, 2011 [11] G X Wang, Semi industrial tests on enhanced underground coal gasification at ZhongLiang-Shan coal mine, Curtin university technology, 2009 [12] M.Wiatowski Kstanczyk, Semi-technical underground coal gasification (UCG) using the shaft method in Experimental Mine “Barbara”, Fuel, 2012 117 [13] Results of the tracer tests during the El-Tremedal underground coal gasification at great depth, Fuel, 2000 [14] Review of Underground Coal Gasificationwith Reference to Alberta's Potential [15] Sateesh Daggupati - Laboratory studies on combustion cavity growth in lignite coal blocks in the context of underground coal gasification, Elsevier, 2009 [16] V Prabu- Simulation of cavity formation in underground coal gasification using bore hole combustion experiments- Elsevier, 2009 [17] Lanhe Yang, Jie Liang, Li Yu, Clean coal technology - Study on the pilot project experiment of underground coal gasification, 2002 [18] Lanhe Yang, Study of the model experiment of blinding - hole UCG, 2002 [19] Ahad Sarraf Shirazi - CFD Simulation of Underground Coal Gasification, 2013 [20] PGS.TS Phạm Lê Dần - Cơng nghệ lị mạng nhiệt.Nhà xuất Khoa Học Kỹ Thuật Hà Nội – 2005 [21] Đồ án anh chị khóa trước [22] COMSOL Mutiphysics Library 118 Review Of Underground Coal Gasification Technologies Nguyen Le Hong Son Nguyen Hoang Anh Hoang Ngoc Dong Da Nang University of Technology Da Nang city, Vietnam Ho Chi Minh City University of Technology and Education Ho Chi Minh city, Vietnam Sonnlh@hcmute.edu.vn , hoanganhskill@gmail.com Abstract—In thewake of increasingchallenges of highprices of oil and gas and uncertaintiesaboutpoliticalstability in many oil and gas producingcountries, coalbecomesmore and more important in the comingyearsforitsvastreserves and widedistributionallover the world The technology of Underground Coal Gasification (UCG), converting in-situ, unmined coal intocombustiblegases, hascontinued to attractworldwideinterestbecause of its ability to exploit coal whichisotherwiseunminablebyconventionalminingtechniquesdue to deepdepositdepths, thin seam thickness or low quality, in an economical, safe and environmentally friendly manner We have recently reviewed the current status of UCG throughout the world and analyzed the criteria for selecting UCG in the Red River Delta (RRD)- Viet Nam This article presents the main results of this work Keywords—underground coal gasification, ucg, syngas, coal seam, unmined coal seams INTRODUCTION not be exploited economically or technically by conventional mining processes The basic principles of UCG are similar to that of conventional surface gasification of coal The UCG process is illustrated in Fig.1 The process of UCG involves the injection of steam and air or oxygen into a coal seam from an injectionwell and recovering combustible fuel gases from the production well These bore wells are drilled into the coal seam from the surface and combustion is initiated at the bottom of one of the bore holes Both bore wells are linked with various available linking techniques to enable gas flow In recent years, people have applied many methods of burning and transfer coal fuel to other fuel types are very effective, it reduce emissions source pollute the environment, as well as transforming coal into liquid fuel, coal washing and especially the coal gasification Coal gasification is a method to transfer coal to gas or used as raw materials for chemical synthesis The quest for energy resources delineates the progress of mankind from the stone age to the modern era For the sustainable development and raising the living standards, man has tried all possible energy resources offered by the nature and used them creatively, for the better quality of life The conversion of available fossil fuel into useful heat energy has been an area of interest for long time This rudimentary practice was encouraged by the abundant availability of the fossil fuels such as coal, natural gas and crude oil The burgeoning energy demand was the main motivation behind the development of various technologies that converted fossil fuel into heat energy The coal gasification is one of the important energy harnessing technologies, which converts coal into useful gaseous fuel The Underground Coal Gasification technology (UCG) is a technology for recovering the energy content of coal reserves by gasifying it in-situ with the application of the skilful utilisation of operating conditions and the geological position of the coal seam This process aims at converting coal into combustible fuel gas by the gasification of the coal seam in the presence of air, oxygen and steam This technology has been proven to be very useful in situations where coal deposits can Figure 1: Principle of underground coal gasification REVIEW OF LITERATURE: Gasification Process- Chemical reactions The gasification process occurs in the coal of the gasifier was shown in Fig Based on the differences in major chemical reactions, the temperature, and the gas compositions, the gasification channel can be divided into three zones: oxidization zone, reduction zone and dry distillation zone 119 Pg    Pg z    Pg  2RTg   r  r r  r  k where, Pg is the fluid pressure, Pa; coefficient, Pa.s;  is the cohesiveness  is a source or convergence item; R is the gas constant, J/mol.K and k is the seepage rate,  m2 Brinkman Equations theory   p    . u   Qbr t u    u    u.     p  t p  Figure 2: The diagram of UCG In the oxidization zone, the multi-phase chemical reactions between oxygen contained in the gasification agent and the carbon in the coal seam occur, producing heat and making the coal seamfull-hot C + O2 → CO2 + 393,8 MJ/kmol (1) 2C + O2 → 2CO + 231,4 MJ/kmol (2) 2CO + O2 → 2CO2 + 571,2 MJ/kmol (3) In the reduction zone, the major reactions are that H2O(g) and CO2 are reduced to H2 and CO under the effect of high temperature, when they meet with the incandescent coal seams C + CO2 → 2CO + 162,4 MJ/kmol (4) C + H2O(g) → CO + H2 + 131,5 MJ/kmol (5) Additionally, under the catalytic action of coal ash and metallic oxides, a certain methanation reaction occurs C + 2H2 → CH4 + 74,9 Heat Transfer and Thermodynamics Species Transport Ts T    Ks s  z  z Ts       K  s r  r     u  F  +  : is the density of the fluid (kg/m3)  : (SI unit) is a dynamic viscosity of the fluid(Pa.s) +  p : is the porosity + k : is the permeability tensor of the porous medium +  F : hệ số Forchheimer + Qbr : is a mass source or mass sink (kg/m3.s)  : is the Laplace When the neglect inertial tern check box is selected, the tern  u.  u p is disabled For incompressible flow, the density stays constant in any fluid particle, which can be expressed Ts    K  s r  FsTs  r  theseequations : + u : is a velocity vector (m/s) + p : is the pressure (Pa) Process- Model of the temperature field of the coal layer Considering the heat conduction within the coal layer, according to the law of conservation of energy, the temperature field equation of the coal layer is obtained as  sCs 1  Q T   p     u   u     .u  I    k 1  br2   p    p  + MJ/kmol (6) In as   p    u.  t   p       u   Qbr reduces to t where and equation .u  Qbr + Cs is the specific heat of the solidphase, KJ/kg.K; + Fs is the heat losscoefficient of the solid phase, KJ/kg.s.K; Heat transfer in a porous media: T   C p u.T  .q  Q  Qvd t q  keff T  C  K  s is the heat conduction coefficient of solid + p eff phase, KW/m.K; + Tsr is the temperature of the external coal layer, K; where: +  : (SI unit: kg/m3) is the fluid density + and rs is the density of the solid phase, kg/m3 Flow equation for the gas phase According to the law of conservation of momentum, for the conditions of UCG, the flow equation of fluid in the gasification channel can be explicitly expressed as + Cp : (SI unit: J/(kg·K)) is the fluid heat capacity at constant pressure + (ρCp)eff (SI unit: J/(m3·K)) is the effective volumetric heat capacity at constant pressure defined by an averaging model to account for both solid matrix and fluid properties 120 + u : (SI unit: m/s) is the fluid velocity field, either an analytic expression or the velocity field from a Fluid Flow interface u should be interpreted as the Darcy velocity, that is, the volume flow rate per unit cross sectional area The average linear velocity (the velocity within the pores) can be calculated as uL = u ∕ θL, where θL is the fluid’s volume fraction, or equivalently the porosity + q : is the conductive heat flux (SI unit: W/m2) + Q : is the heat source (or sink) Add one or several heat sources as separate physics features (W/m3) + k eff Chinese studies have been mainly focused on the production of hydrogen-rich syngas from abandoned shallow coal mines (Depth 300 m ) of lower rank coals and few efforts for using highrank coals (anthartice and bituminous coals) were unsuccessful In the following sub-sections, a brief review of important field trials and their findings are described based on their location Figure summarizes all the UCG field trials performed worldwide in term of their depth and thickness of the coal seam As can be seen, all of the trials except European efforts have been conducted in relatively shallow seams which are not currently targeted because shallow depth of seam limits the application of high operating pressures and increases the possibility of leakage Figure 4: The thickness and the depth of the test is performed all over the world (C.Pana, 2009) Figure 5: Long tunnel configuration in China (Couch (2009)) Yeary and Riggs (1987) calculated the growth for lignite and sub-bituminous coal Recession rate of cavitysidewall has beenshown to have a strongdependenceon the flow rate and temperature of the injected gas Lignite recession rate has been reported to begenerallylowercompared to subbituminous, while having a higher spalling rate (Yeary and Riggs 1987) These observationsareconsistentwithweakmechanicalproperties of lignite and high reactivity of sub-bituminous coal Also, ash has remainedintact in the sidewall for lignite coal, which can be attributed to higher ash content and stronger ash structure in lignite sample A series of experiments have been conducted in European Union, mostly Poland, to investigate the UCG behaviour in lignite coal seams The mainobjective of thisongoingproject is to study the effect of variousparameters on hydrogenconcentration in syngas Experiments wereconducted in a 2.5*0.7*0.7 m coal block with a 10 cmlinkbored in the block Three distinctstages for the UCG process were observable in all of these experiments: ignition, combustion, and gasification with steam As reported by Stańczyk etal (2010), only hr of stablesteamgasificationwasachieved in theirfirst experiments and further gasification was impossible due to low temperature of the coal block and fastdrop in temperature due to high moisturecontent of the experimented lignite (53%) To compensate this, the experiment was continued with high injectionrates of oxygen Replacing oxygen with air was unfeasible due to rapiddecrease in temperature and termination of reactions (Stańczyk et al., 2011) Thus, oxygenrich air was used and optimum oxygen/air ratios for lignite and hard coal UCG were proposedaccordingly Shorter ignition and combustion periods, higher 121 temperatures and lower oxygen/air ratios were found to be attainable for hard coal due to its highercarbon content compared to lignite In order to increase the extent of steam-gasification injection of pure oxygen was continued until high temperature of about 1100 to 1200 oC was reached in the reactionfront (Stańczyket al., 2012) Then, the inlet gas was switched to steam and the heating value of the product syngas increased to around 11.5 MJ/m3 as the injection gas is changed to steam The inlet gas was changed to oxygen whenever temperature dropped to 700-800 °C in the reaction front Instead of cuttingthrough the block, the georadartechnique which is based on the difference in dielectric properties was used in these experiments to identify various zones in the reacted coal Three zones of free-space cavity, partially gasified lignite, and un-reacted dry coal were detected In the experiments performed in IIT Bombay, India, the cavity growth, cavity shape, and product gas composition with respect to various operatingconditions were studied (Daggupati et al., 2010) These experiments were undertaken in a 30*20*25 cm lignite coal block A 3mmdiameter link was drilled to connect the injection and production wells In first series of experiment, cavity growth and finalshape of cavity is of majorconcern Ignition has been initiated using LPG for minutes and gas is switched to pure oxygen to combust the coal Experiments were continued for a maximum of hours with different injection rates and welldistances After termination of experiments, coal block was opened to expose and study the cavity shape Similar to the field trials, the cavity was found to be teardrop shape and symmetricaround the injection point Cavity growth in upstream of injection point is less than width and height of the cavity Common cavity temperature is around 950 to 1000 oC (Daggupati et al., 2010) Increasing flow rate resulted in linear increase in cavity volume in all dimensions This increase is attributable to the effect of increased flow rate on decreasing the mass transfer limitation by removing the ash at coal surface The forward and backward growth and volume of the cavity was found to decrease with increasing the distancebetween the wells, while cavity width and height increase When the cavity reaches the production well, considerable bypass of oxygen terminates the process Stable operation of similar experiments with a mixture of steam and oxygen found to be possible only by the injection of high temperature steam (600-700 oC) in a cyclic manner (Daggupati et al., 2011) Cavity shape was the samebutlarger in case of steam/oxygen injection Coal particles of 1-2 cm in diameter were alsoobservedat the bottom of cavity in this case, while only ash was found in combustion Therefore, it can be concludedthat coal spalling is increased in gasification The best syngas heatingvalue was achieved by a steam/oxygen ratio of 2.5 Prabu and Jaynati (2011) investigated the cavity formation in blocks of wood, coal, and camphor and compared the results In all these cases, teardrop-shape cavities were formed with a more bulbous shape at lower flow rates One of the major challenges for applicability of UCG in steep coal seams is low heating value of produced gas Yang et al (2003) studied the UCG in a 6.8 m*0.25 m*1.1 m coal block with a dipangle of 65° with pure oxygen or mixture of oxygen and steam as injection gases In the case of oxygen/steam mixture, the flow rate of the produced gases decreasedafter 41 hours and injecting the gases close to reaction front was required to sustain the process An auxiliary well between injection and production well was used to supply the gas for the reaction front A steam/oxygen ratio of was found to give the best syngas heating value Using multiple injection points increased the duration of the gasification stage, while shortened the oxidation phase CONCLUSION The, worldwide consumption of coal by UCG is listed below: a Soviet Union: 1500 Ktonnes since 1950 (Cost: US$10 billion) b Europe: less than 10 Ktonnes since 1950 (Cost: US$1 million) c Australia: over 32 Ktonnes between 200-2002 (Cost: US$5 million) UCG trials have been performed over a 50 years period and in different geologicalsettings and various coal ranks The knowledgedeveloped based on these activities can be used to select the appropriate seam and designaneffectiveprocedure Following results can be concluded from the trials: - Depth of the seam dictates the operating pressure which is one of the most important parameters in UCG process Considerable leaks (Up to 20-30 percent of gas) have been reported in shallower seams - Heating value of the produced gas can be increased drastically, by using oxygen-enriched air - Syngases produced from CRIP technology generally have a higher heating value It should be noted that these experiments reviewedhere have only studied shortterm and earlybehaviour of UCG reactor and the results may or may not be representative of a long-term UCG operation In Vietnam, especially in the RRD, unmined coal reserves are very large so coal mining by UCG should be studied in detail References [1] Patankar S.V., Numerical Heat Transfer and Fluid Flow, McGraw Hill, 1980 [2] Smoot D.L., Smith P.J., Coal Combustion and Gasification, Plenum Press, 1985 [3] Gasification Technology- Technical Issues in the Design of Gasifiers1999 [4] Krzysztof Stanczyk, Krzysztof Kapusta-, Experimental simulation of hard coal underground gasification for hydrogen production, Fuel, 2012 [5] Krzysztof Stanczyk Hydrogen-oriented underground coal gasification for Europe (HUGE)-Euro commission 122 [6] Krzysztof Stanczyk, Krzysztof Kapusta, Pollution of water during underground coal gasification of hard coal and lignite, Fuel, 2011 [7] G X Wang, Semi industrial tests on enhanced underground coal gasification at Zhong-Liang-Shan coal mine, Curtin university technology, 2009 [8] M.Wiatowski Kstanczyk, Semi-technical underground coal gasification (UCG) using the shaft method in Experimental Mine “Barbara”, Fuel, 2012 [9] Results of the tracer tests during the El-Tremedal underground coal gasification at great depth, Fuel, 2000 [10] Review of Underground Coal Gasificationwith Reference to Alberta's Potential [11] Prabir Basu- Combusion and gasification- Taylor and Francis [12] V Prabu- Simulation of cavity formation in underground coal gasification using bore hole combustion experiments- Elsevier L H Yang, Study on the Methods of Seepage Combustion in the Process of Underground CoalGasification,pp.72– 75,ChinaUniversityofMiningandTechnologyPress,Xuzhou, 2001 J W Martin, A W Layne, and H J Siriwardane, Thermo-Mechanical Modeling of Ground Movements Associated with Large Underground Coal Gasification Cavities in Thin Coal Seams, Proc 10th UCG Symp., Fallen Leaf Lake, CA, vol 1, pp 295–307, 1984 O.K.Min,FiniteElementModellingofThermoMechanicalResponsesAssociatedwith Underground Coal Conversion, Ph.D thesis, Ohio State University,1982 C B Thorsness and J A Britten, Analysis of Material and Energy Balances for the RockyMountainIUCGFieldTest,CollectionsofTranslatedWorkson Mining,vol.13, pp 7–13,1992 K Guntermann, The Model and Mathematical Models on the Comprehensive Experi- ment of Underground Coal Gasification, Collections of Translated Works on Mining, vol 9, pp 6–10,1988 C W Guo, Review of the Study of Mathematical Simulation on the Underground Coal Gasification, Mining World, vol 15, pp 3– 5,1994 S H Advani, J K Lee, and K S Chen, Geomechanical Modeling Associated with Underground Coal Gasification Processes, Proc 12th UCG Symp., U.S Department of EnergyAssistantSecretaryforFossilEnergy,Washington,DC,vol.1,pp.2 78–286,1986 [7] Gasification Technology, Technical Issues in the Design of Gasifiers, 1999 [8] Krzysztof Stanczyk, Krzysztof Kapusta, Experimental simulation of hard coal underground gasification for hydrogen production, Fuel, 2012 [9] Krzysztof Stanczyk Hydrogen-oriented underground coal gasification for Europe (HUGE)-Euro commission [10] Krzysztof Stanczyk, Krzysztof Kapusta, Pollution of water during underground coal gasification of hard coal and lignite, Fuel, 2011 [11] G X Wang, Semi industrial tests on enhanced underground coal gasification at Zhong-Liang-Shan coal mine, Curtin university technology, 2009 [12] M.Wiatowski Kstanczyk, Semi-technical underground coal gasification UCG) using the shaft method in Experimental Mine “Barbara”, Fuel, 2012 123 NGHIÊN CỨU ỨNG DỤNG PHẦN MỀM MƠ PHỎNG Q TRÌNH TRUYỀN NHIỆT- TRUYỀN CHẤT TRONG KHAI THÁC THAN SỬ DỤNG PHƯƠNG PHÁP KHÍ HÓA THAN NGẦM Applied research numerical simulation on heat and mass transfer in the processof underground coal gasfication Nguyễn Lê Hồng Sơn 1, Nguyễn Hoàng Anh2 Trường Đại Học Sư Phạm Kỹ Thuật Hồ Chí Minh; hongsonnguyenle@gmail.com Trường Cao Đẳng Nghề Đồng Nai; hoanganhskill@gmail.com TĨM TẮT ABSTRACT Khí hóa than ngầm q trình chuyển than đá trực tiếp vỉa than thành khí đốt dùng làm nguyên liệu tổng hợp hoá chất Nghiên cứu thực mơ chi tiết q trình khí hố than ngầm bên vỉa than để sản xuất khí tổng hợp (CO, CO2) Các mô thực cách sử dụng phần mềm Comsol Multiphysics phiên 5.2 để mơ q trình truyền nhiệt truyền chất q trình khí hố than ngầm Các nghiên cứu truyền nhiệt mơ hình dạng khoang rỗng ba chiều: chiều dài, chiều rộng chiều sâu sau q trình khí hố thể vỉa than ngầm Cần lưu ý nhiêt độ bên vỉa than cao nhiều so với nhiệt độ đầu khí sản phẩm, điều tổn thất nhiệt môi trường xung quanh Các khoang rỗng có hình giọt nước, mô dựa chiều rộng, chiều cao, chiều dài phía trước phía sau khoang Điều chỉnh nồng độ O2 cấp vào q trình khí hố ảnh hưởng đến chất lượng khí tổng hợp nghiên cứu Tỷ lệ cấpO2 tốt khoảng 0,15 mol/ m3 Underground Coal Gasification (UCG) is a process in which coal is converted to clean synthetic gas (syngas) insitu This study performs detailed simulations of UCG process inside coal seam to produce syngas (CO, CO2) The simulations are conducted using the Comsol Multiphysics solfware which is simulated the heat and mass transfer during the Underground coal gasification Studies on the heat transfer and the growth of threedimensional cavity geometries in Underground coal gasification are presented in underground coal seam The heat tranfer should be noted that local temperature in the coal seam is higher than the outlet temperature, this might be attributed to gas heat loss to surrounding enviroment The cavity has a tear-drop shape, which could be characterized based on its width, height, forward and backward lenght of the cavity The effects of oxygen injection rate on syngas quality also studied The best oxygen injection rate is about 0,15 mol/ m3 Từ khố: khí hố than ngầm; UCG; truyền nhiệt truyền chất UCG; comsol multiphysics; khí tổng hợp ĐẶT VẤN ĐỀ Q trình khí hóa than ngầm thường diễn thời gian dài, với khối lượng than lớn gấp nhiều lần so với mô hình thực nghiệm Do phải đầu tư khoản chi phí lớn cho q trình xây dựng hệ thống khí hóa phải nghiên cứu kỹ loại than, dự đốn trường hợp xảy Trước u cầu đó, mơ q trình khí hóa than ngầm thiết kế xây dựng gần giống với thực tế nhằm thử nghiệm trước xây dựng hệ thống lớn biết chắn điều Keywords: underground coal gasification, ucg, heat and mass transfer, syngas, oxygen concentration, coal seam, Comsol Multiphysics mang lại lợi nhuận kinh tế cao Với mơ hình này, q trình khí hóa than diễn gần giống lòng đất dựa sở lý thuyết chung cho cơng nghệ khí hóa Với đường ống dẫn hỗn hợp khí + nước đường ống thu hồi khí tổng hợp Qua đó, nghiên cứu xác định thành phần lí hố học, chất lượng khí thương phẩm từ khí sản phẩm thu Trong thời đại khoa học công nghệ ngày nay, công nghệ mô số ngày sử dụng rộng rãi lĩnh vực khoa học kỹ thuật hoạt động 124 Lớp than khoang khí hóa liên tục người Mơ số mang đến cho người lợi ích to lớn tiết kiệm thời gian, kinh phí, nguyên vật liệu, tránh rủi ro điều kiện thực tế, giảm tác động xấu tới môi trường Các mơ hình xây dựng máy tính mơ q trình trước mơ hình thực nghiệm xây dựng Ở đây, nhóm tác giả xin giới thiệu phương pháp mơ q trìnhtruyền nhiệt truyền chất khai thác than sử dụng phương pháp khí hóa than ngầm qua ứng dụng phần mềm mơ Comsol Multiphysics phiên 5.2 Một mơ hình tốn học truyền nhiệt truyền chất q trình khí hóa than xây dựng theo phương pháp bảo tồn lượng khối lượng q trình sản xuất khí hóa than Việc nghiên cứu giá trị mơ kết hợp với q trình thực nghiệm khí hóa than cung cấp sở lý thuyết cần thiết cho nghiên cứu hay dự đốn quy luật diễn q trình khí hóa than việc xác định yếu tố làm ảnh hưởng đến hiệu suất khí hóa Trong đó: CƠ SỞ LÝ THUYẾT Cgi: nồng độ mol thành phần i (mol/m3) Các trình cháy khí hóa lịng đất lớp than phức tạp mặt vật lý phản ứng hóa học với q trình truyền nhiệt truyền chất pha rắn- khí Để dễ dàng q trình tính tốn, ta có điều kiện sau: 2.1 Các phản ứng hóa học khí hóa than: Vg: vận tốc dịng khí dọc trục (m/s) Q trình oxy hóa than: C + O2 = CO2 + 399499 kJ/mol (1) Hình 1: Mơ hình tính tốn học kênh khí hóa Cgi    Cgi    Cgi  Cgi   (Cgi.Vg )  D   D   D  Si (5)  z z  z  r  r  r r D: hệ số khuếch tán dòng khí (trục bán kính khuyếch tán thành phần chất khí) (m2/s) 5: thời gian (s) Si: tỷ lệ sản phẩm thành phần thứ i phản ứng hóa học (mol/ m3 s) Nếu chiều dài kênh khí hóa dài chiều dài bán kính phương trình viết đơn giản sau: Phản ứng Boudouard: C + CO2 = 2CO + 167023 kJ/mol (2) C gi Q trình khí hóa cacbon với nước: C + H2O = CO + H2 –125788 kJ/mol r (3)  C    C gi V g  D gi z  z    S i  Quá trình tạo metan: Điều kiện biên: C + H2 => C H4 (4) 2.2 Phương trình bảo tồn thành phần chất pha khí: Trong q trình đốt khí hố lớp than, nhiều phản ứng vật lí hố học phức tạp xảy lị khí hố với trình truyền nhiệt truyền chất pha rắn- khí Để dễ dàng q trình tính tốn, ta giả định điều kiện sau: Kênh khí hóa làm việc ổn định, thông số vật lý, nhiệt động hệ số dẫn nhiệt, nhiệt dung riêng, hệ số trao đổi nhiệt không thay đổi theo thời gian Bỏ qua trình khuếch tán nhiệt áp suất Bỏ qua trở nhiệt lực liên kết Z=0, Cgi= Cgi0, ≤ ≤ 71 5>0 Cgi0: nồng độ mol thành phần I phía đầu kênh khí hố (mol/m3) Z=L, C gi z 0 ≤ ≤ 71 5>0 L: chiều dài kênh khí hóa, (m) Điều kiện ban đầu trước khí hóa: ≤0, Cgi= Cgi0, ≤ ≤ 71 0≤8≤9 2.3 Phương trình bảo tồn lượng pha khí: 125 Theo định luật bảo tồn chuyển hóa lượng, phương trình bảo tồn lượng cho pha khí:   (V g Tg )  y gi   C gi C pgi r   C pgi C gi z  C g DTg z    Ri H i    A(Tg  Ts )  Q g Tg r = r0, Ts = Tsr0 0≤z≤L ; =0 z = 0, Ts = Tsz0 5>0 5>0 5>0 r1 < r < r0 Trong đó: (6) Tsr0: Nhiệt độ bên ngồi lớp than (tại vị trí r0) số Tại bề mặt pha rắn khí (r= r1), điều kiện biên là: Tg: Nhiệt độ pha khí (K) Cg: Tổng nồng độ mol pha khí (mol/m ) 4•ž Cpgi: Nhiệt dung riêng chất khí thứ I (KJ/mol.K) Ri: Tốc độ phản ứng hóa học (mol/ m s) Hi: Nhiệt lượng phản ứng hóa học (KJ/mol) Ÿ¡ = h(Ts – Tg) - ∑DE E ¢|E (8) Trong đó: mi: Lưu lượng thể thích thành phần thứ i (kg/m3) u: Hệ số trao đổi nhiệt pha rắn khí (KJ/m2) Ÿ ž ¢|E : Nhiệt phản ứng hóa học chất thứ i (KJ/mol) A: Diện tích tiếp xúc pha rắn khí (m / m ) N: tổng số chất Qg: Nhiệt lượng tổn thất pha khí (KJ) Điều kiện ban đầu: = 0, Ts = Ts0 ≤ z ≤ L ≤ r ≤ r1 ygi: Thành phần phần trăm chất khí thứ i pha khí Ts0 : Nhiệt độ ban đầu pha rắn (K) Điều kiện biên: z = 0, Tg = Tgz0; ≤ r ≤ r1 > z = L, ; = 0 ≤ r ≤ r1 > ; 2.5 Phương trình dịng chảy pha khí: Theo định luật bảo tồn động lượng, phương trình dịng chảy kênh khí hóa: 5>0 Pg  = 0,Tg = Tg0 ≤ z ≤ L ≤ r ≤ r1   Pg z    Pg  2RTg  r   (9) r r  r  k Trong đó: Trong đó: Pg  Cg RTg Tggz0, Tgr0: Nhiệt độ biên pha khí (K) Tg0: Nhiệt độ ban đầu pha khí (K) Pg: Áp suất dịng khí (Pa) Kλg: Hệ số dẫn nhiệt pha khí (KW/m K)  : Hệ số dính kết (Pa.s) Ts: Nhiệt độ pha rắn (K)  : Hệ số hiệu chỉnh h: hệ số trao đổi nhiệt đối lưu (KJ/m2 K) R: Hằng số chất khí (J/mol.K) 2.4 Phương trình vùng nhiệt độ lớp than: Theo định luật bảo tồn chuyển hóa lượng, xét q trình dẫn nhiệt lớp than, ta có: T T    T  T    s C s s   K  s s    K s s   K  s s  z  z  r  r  r r  Fs T s (7) Cs: Nhiệt dung riêng pha rắn (KJ/kg.K) Fs: Nhiệt lượng tổn thất pha rắn (KJ/kg.s.K) K  s : Hệ số dẫn nhiệt pha rắn (KW/m K)  s : Khối lượng riêng chất rắn (kg/m3) k: tỷ lệ rị rỉ chất khí (  m 2) MƠ HÌNH MƠ PHỎNG Nhằm thử nghiệm trước đầu tư xây dựng hệ thống lớn với mức chi phí lớn cho q trình xây dựng hệ thống khí hóa phải nghiên cứu xây dựng mơ hình mơ dự đốn trường hợp xảy để chắn nghiên cứu ứng dụng thực tế mang lại lợi nhuận kinh tế cao Với mơ hình này, q trình khí hóa than diễn gần giống lòng đất dựa sở lý thuyết chung cho công nghệ khí hóa 3.1 Mục đích Điều kiện biên: 126 - Mơ q trình truyền nhiệt khí hóa than 3.3 Cài đặt điều kiện thông số đầu vào cho mơ - Mơ hình dạng lỗ rỗng - Xét ảnh hưởng hiệu suất khí hóa phụ thuộc vào nồng độ O2 - Nhiệt độ mơ phỏng: 7500C 3.2 Xây dựng mơ hình - Vận tốc cấp O2 ban đầu từ 0,01 m/s Xét khối than có kích thước 200x100x40(cm) xây dựng phần mềm Comsol phiên 5.2 giả định điều kiện: Bảng 1: Thành phần than - Áp suất mơ từ bar - Kênh khí hóa làm việc ổn định, thông số vật lý, nhiệt động hệ số dẫn nhiệt, nhiệt dung riêng, hệ số trao đổi nhiệt không thay đổi theo thời gian - Bỏ qua trình khuếch tán nhiệt áp suất - Bỏ qua trở nhiệt lực liên kết KẾT QUẢ MÔ PHỎNG - Bỏ qua thay đổi khối lượng dịng khí phản ứng hóa học vùng oxy hóa 4.1 Q trình truyền nhiệt khí hóa than - Lớp than khoang khí hóa liên tục - Trên khối than xây dựng hai đường ống: Đường ống dẫn khí nước vào khối phản ứng có d = 21cm vị trí 50x50x20(cm) khối than Đường ống khí thải sản phẩm khỏi khối khí hóa có d = 34 cm vị trí 150x50x20(cm) khối than Nhiên liệu cấp vào lị mơ khối than tiếp xúc trực tiếp với oxy Lúc đầu, than phản ứng với oxy tạo nhiệt độ cao, trình cháy diễn Sau đó, nhờ có nhiệt độ cao cháy yếm khí mà phản ứng thu nhiệt q trình khí hóa xảy ( Phản ứng chuyển hóa carbon thành CO) Để xác định ảnh hưởng thơng số đến q trình khí hóa tác giả tiến hành mơ theo mơ sau: Hình 3: Kết mơ q trình truyền nhiệt lúc tâm phản ứng đạt 7500C vecto vận tốc Hình Mơ hình mơ khí hóa than dạng 2D 127 Hình 6: Hình dạng khoang rỗng hình giọt nước hình thành Hình4: Phân bố nhiệt độ vecto vận tốc khoang khí hóa dùng phần mềm CFD Kết mơ hình dạng lỗ rỗng hình thành khoang rỗng hình giọt nước Điều phù hợp với kết hình dạng khoang [4]và phù hợp với kết [9] Quá trình mơ giúp ta dễ dàng việc hồn thổ sau khai thác Hình 5: Dạng hình học dự đốn khoang rỗng dùng phần mềm CFD Hình mơ vecto vận tốc q trình truyền nhiệt khoang khí hóa Ta thấy rằng, nhiệt độ tâm q trình khí hóa cao nhất, giảm dần hướng lỗ thu khí Từ đó, ta dự đốn nồng độ Carbon thay đổi nhiệt độ thay đổi vị trí khác vùng khí hóa Comsol Multiphysics thể hình thành dịng chảy phức tạp thơng qua vecto vận tốc khoang khí hóa thể rõ ràng qua hình, điều phù hợp với kết nghiên cứu S Daggupati [15] 4.2 Hình dạng khoang rỗng sau q trình khí hóa than Hình 7: Hình dạng khoang rỗng 4.3 Ảnh hưởng nồng độ O2 đến thành phần sản phẩm khí: Khi ta thay đổi nồng độ Oxy cấp vào nồng độ sản phẩm khí thay đổi theo Theo điều trên, thành phần chất khí CO, CO2 thay đổi theo nồng độ O2 sau: 128 Hình 8: Thành phần sản phẩm khí thu cho nồng độ O2 0,10 mol/ m3 Hình 12: Thành phần sản phẩm khí thu cho nồng độ O2 0,60 mol/ m3 Hình 13: Thành phần sản phẩm khí thu cho nồng độ O2 0,75 mol/ m3 Hình 9: Thành phần sản phẩm khí thu cho nồng độ O2 0,15 mol/ m3 Hình 14: Thành phần sản phẩm khí thu cho nồng độ O2 1,0 mol/ m3 Hình 10: Thành phần sản phẩm khí thu cho nồng độ O2 0,30 mol/ m3 Kết từ hình 8, 9, 10, 11, 12, 13, 14 cho thấy, thực q trình khí hóa, ta tăng nồng độ O2 từ 0,10 - 0,30 mol/s nồng độ sản phẩm khí sinh CO2 CO bắt đầu tăng theo Tại nồng độ O2 0,15 mol/s, hình nồng độ CO cao CO2 thấp Khi nồng độ O2 tăng từ 0,45 - 1,0 mol/s nồng độ CO2 bắt đầu tăng cao CO giảm Điều phù hợp ta tăng lượng oxy lên cao lúc trình cháy diễn lượng CO2 sinh nhiều hơn, lượng CO giảm Quá trình phủ hợp với nghiên cứu Lanhe Yan [14] Hình 11: Thành phần sản phẩm khí thu cho nồng độ O2 0,45 mol/ m3 129 KẾT LUẬN Mơ q trình khí hóa than ngầm dùng phần mềm Comsol Multiphysics ta mơ q trình truyền nhiệt, truyền chất Qua đó, ta hiểu rõ trình truyền nhiệt diễn khoang khí hóa Từ đó, ta dự đốn hình thành khoang rỗng hình giọt nước Quá trình giúp ích cho q trình hồn thổ sau Phương pháp mơ khí hóa than ngầm cho phép cài đặt thông số ban đầu để sản xuất khí than có chứa thành phần khí hóa mong muốn, điều chỉnh thơng số ban đầu để kiểm sốt khí sản phẩm, có lợi dùng vào mục đích khí hóa than ngầm thực tế Thành phần khí than phụ thuộc vào nồng độ O2 cấp vào Sản phẩm có giá trị nhiệt lượng cao nồng độ O2 từ khoảng 0,15- 0,30 mol/s [10] Results of the tracer tests during the El-Tremedal underground coal gasification at great depth, Fuel, 2000 [11] Review of Underground Coal Gasificationwith Reference to Alberta's Potential [12] M.S Blinderman, Forward and reverse combustion linking in underground coal gasification- Combustion and Flame 2007 [13] M.S Blinderman, Forward and reverse combustion linking in underground coal gasification- Combustion and Flame 2007 [14] Lanhe Yang- Study on the model experiment and numerical simulation for underground coal gasification, Elsevier, 2003 [15] Sateesh Daggupati - Laboratory studies on combustion cavity growth in lignite coal blocks in the context of underground coal gasification, Elsevier, 2009 TÀI LIỆU THAM KHẢO [1] Đặng Quốc Phú, Trần Thế Sơn, Trần Văn Phú, Truyền nhiệt – NXB giáo dục – 2004 [2] PGS.TS Bùi Hải – PGS.TS Trần Thế Sơn (2001) Bài Tập Truyền Nhiệt – Nhiệt Động Và Kĩ Thuật Lạnh, Nhà xuất khoa học kĩ thuật Hà Nội [3] Nguyễn Thanh Quang (ĐHBK – Đà Nẵng) – Đặng Thế Hùng (Công ty TNHH Trường Quang II) “Nghiên Cứu Chế Tạo Hệ Thống Hóa Khí Than Tầng Cố Định Ngược Chiều” Tạp Chí Khoa Học Và Cơng Nghệ Nhiệt Số 77 [4] Prabir Basu- Combusion and gasification- Taylor and Francis [5] Gasification Technology- Technical Issues in the Design of Gasifiers-1999 [6] Krzysztof Stanczyk Hydrogen-oriented underground coal gasification for Europe (HUGE)-Euro commission [7] Krzysztof Stanczyk, Krzysztof Kapusta, Pollution of water during underground coal gasification of hard coal and lignite, Fuel, 2011 [8] G X Wang, Semi industrial tests on enhanced underground coal gasification at Zhong-Liang-Shan coal mine, Curtin university technology, 2009 [9] M.Wiatowski Kstanczyk, Semi-technical underground coal gasification (UCG) using the shaft method in Experimental Mine “Barbara”, Fuel, 2012 130 ... phương pháp khí hóa than khác Tùy thuộc kích cỡ than đưa vào lị khí hóa mà áp dụng ba kiểu cơng nghệ khí hóa phổ biến là: khí hóa than tầng cố định; khí hóa than tầng sơi khí hóa than dịng - Than cục... thuyết than sử dụng làm thực nghiệm Mơ q trình khí hố than ngầm phần mềm Comsol Multiphysics Thực nghiệm trình khí hố than ngầm mơ hình thực nghiệm Đánh giá, phân tích kết thu q trình mơ thực nghiệm. .. Đề tài ? ?Nghiên cứu mơ thực nghiệm q trình khí hố than ngầm ’sẽ đưa giải pháp nhằm đánh giá cụ thể việc khí hóa than ngầm ĐBSH sở số liệu nghiên cứu lý thuyết thực nghiệm mơ hình khí hóa cụ thể