Luận án nghiên cứu phát triển công nghệ yếm khí cao tải tuần hoàn nội ic (internal circulation)

1 MỞ ĐẦU Tính cấp thiết luận án Hiện nay, khoảng 80% lượng nước thải xử lý công nghệ sinh học [1], bao gồm: nước thải sinh hoạt, nước thải công nghiệp chứa hợp chất hữu sinh hủy được, mức độ định công nghệ sinh học xử lý hóa chất độc Trong xử lý nước thải, cơng nghệ sinh học ln ưu tiên chi phí thấp, sản sinh bùn, sử dụng hóa chất, … đặc biệt thu lượng tái tạo dạng khíbiogas Như vậy, xử lý yếm khí không công cụ để bảo vệ môi trường mà bảo vệ tài nguyên Trong lịch sử phát triển, kỹ thuật vi sinh yếm khí phát triển dần từ hệ đơn giản (bồn kín, bể phốt, làcác cơng trình xử lý sử dụng tượng tự nhiên) đến kỹ thuật mang nhiều yếu tố “bể phản ứng” hiệu Theo bước phát triển công nghệ, có loại kỹ thuật xử lý yếm khí ứng dụng là: (1) bồn phản ứng khuấy trộn – CSTR, (2) kỹ thuật phản ứng tiếp xúc AC, (3) kỹ thuật dòng chảy ngược qua lớp bùn yếm khí – UASB (4) kỹ thuật tầng bùn hạt giãn nở - EGSB vàkỹ thuật tuần hoàn nội - IC Trong khoảng 40 năm trở lại nhiều cơng trình áp dụng kỹ thuật yếm khí cao tải với thiết kế khác áp dụng [2] Các tài liệu tham khảo cho thấy IC với EGSB thiết bị xử lý yếm khí cónăng suất cao tiên tiến nay, kỹ thuật IC có khả chịu tải cónăng suất cao kỹ thuật có [3] Do tốc độ tăng trưởng hệ IC cao Năm 1997 giới có 32 hệ IC, tới năm 2007 giới có tới 2266 hệ, đến năm 2015 ước tính số lượng hệ yếm khícao tải lắp đặt vượt 4000 hệ [4] Một yếu tố quan trọng kỹ thuật xử lý vi sinh khuấy trộn, tăng khuấy trộn làm tăng khả tiếp xúc vi sinh vật với chất ô nhiễm nên hiệu chuyển hóa tăng Các kỹ thuật UASB EGSB sử dụng lượng khí sinh dịng chất lỏng tuần hoàn (bằng bơm tuần hoàn) để tăng khả khuấy trộn Kỹ thuật IC sử dụng khí sinh để tuần hồn nước Đây ưu làm cho IC có hiệu xử lý cao lại tiêu thụ lượng (vìkhơng dùng bơm tuần hồn) Để nghiên cứu chế tạo thành cơng hệ IC, điểm mấu chốt cần nghiên cứu chế tạo cấu tự tuần hồn khơng cần lượng bên ngồi (bơm tuần hồn), sử dụng q trình kéo nước khí tự sinh để tuần hồn hỗn hợp bùn nước thải Hiện nay, tài liệu hướng dẫn thiết kế hệ UASB tương đối phổ biến tài liệu kỹ thuật tuần hồn nước IC cơng bố Ví dụ, sách chuyên khảo công nghệ sinh học xử lý nước thải “Handbook of Biological Wastewater Treatment, 2012” [5] có tới 14 ví dụ tính thiết kế hệ UASB kèm lời giải chi tiết mà khơng có ví dụ hệ tiên tiến kiểu EGSB IC Do vậy, để thúc đẩy ứng dụng hệ IC Việt Nam cần phải có nghiên cứu để xác định thông số phục vụ việc thiết kế chế tạo thiết bị Trong luận án tập trung mơ q trình tuần hồn nội (khí kéo nước) nhằm xác định lượng nước (QN) kéo lên lượng khí (QK) (K = QN/QK) khả khuấy trộn khí sinh nước tuần hồn nhằm phục vụ tính tốn cấu khí kéo nước (airlift) hệ IC Trong luận án trình bày kết thử nghiệm chế tạo vận hành mơ hình hệ IC quy mơ phịng thí nghiệm nhằm xác định lực xử lý hệ IC vận hành hệ thống với nước thải chăn nuôi lợn Mục tiêu nghiên cứu Nghiên cứu ứng dụng hệ xử lý yếm khí tuần hoàn nội (IC) xử lý nước thải giàu hữu cơ, cụ thể là: (i) Xác lập quan hệ thơng số thiết kế thường dùng (ví dụ: tải lượng, tốc độ nước dâng) với thông số thiết kế hệ IC (ví dụ: kích thước ống khí -nước lên, chiều cao vùng phản ứng) (ii) Xác định mối tương quan tải lượng suất xử lý hệ yếm khíIC trường hợp NTCNL (iii) Xác định khả khuấy trộn hệ khí sinh (iv) Xác định giá trị thông số (chiều cao thiết bị, chiều cao vùng phản ứng, đường kính thiết bị, kích thước ống lên) phục vụ thiết kế hệ yếm khíIC Nội dung nghiên cứu Nội dung 1: Thực nghiệm xác định tỷ lệ lượng nước kéo lên lượng khí ( K  QN ) mức ngập nước (H1), chiều cao xả nước (H2) vàtiết diện QK ống dẫn lên khác độ nhớt khối lượng riêng xác định dung dịch Nội dung 2: Thiết lập phương trình biểu diễn mối liên hệ K  QN với QK thông số: mức ngập nước (H1), chiều cao xả nước (H2), tiết diện ống dẫn lên (S) độ nhớt khối lượng riêng xác định dung dịch Nội dung 3: Chế tạo vận hành hệ IC với nước thải chăn nuôi lợn nhằm xác định suất xử lý hệ IC Nội dung 4: Tính tốn xác định khả khuấy trộn khí nội sinh nước tuần hồn để xác định thơng số thiết kế hệ IC Những đóng góp khoa học tính luận án Với mục tiêu nghiên cứu xác định giá trị thông số thiết kế công nghệ hệ thống IC, kết thu từ nghiên cứu gồm: (i) Phương trình lý thuyết tính tốn đại lượng K = QN/QK (tỷ lệ lượng nước kéo lên lượng khí) phụ thuộc vào tiết diện ống lên (S), mức ngập nước (H1), chiều cao xả nước (H2), độ nhớt (µ) vàkhối lượng riêng dung dịch (ρ) (ii) Phương trình kinh nghiệm để tính tốn đại lượng K = QN/QK (tỷ lệ lượng nước kéo lên lượng khí) phụ thuộc vào tiết diện ống lên (S), mức ngập nước (H1), chiều cao xả nước (H2) (iii) Từ việc định lượng khả khuấy trộn hệ xử lý khí nội sinh (ở khí biogas) dịng tuần hồn nội thơng qua cơng suất P (W/m3) vàgradient tốc độ G (s-1) xác định thông số thiết kế hệ IC (chiều cao thiết bị, tải lượng, tốc độ dâng) Giá trị thực tế ứng dụng kết luận án Trong xử lý nước thải xử lý yếm khí ngày đóng vai trị quan trọng khơng chi phí đầu tư vận hành thấp mà cịn lực xử lý ngày cao tiềm thu hồi lượng Hệ yếm khí cao tải tuần hoàn nội (IC) ngày áp dụng rộng rãi Nghiên cứu cung cấp thông tin liên quan đến thiết kế hệ IC khả ứng dụng hệ IC CHƯƠNG I: TỔNG QUAN 1.1 Các q trình vi sinh yếm khí Mặc dù tượng phân hủy yếm khí chất hữu thành khí cháy (ngày gọi biogas) Van Helmont phát từ kỷ 17 tới năm 1890 Omelianski phân lập thành công vi khuẩn tạo hyđrô, axit axetic butyric q trình lên men xellulo Ơng người giả định phản ứng vi sinh tạo metan từ hyđrô carbon đioxit, điều mà tới 1910 Sohngen khẳng định, đồng thời với phản ứng tạo metan từ axit axetic qua đường giải phóng CO2 [6] Vấn đề phân hủy yếm khí bắt đầu quan tâm thực từ năm 1930, nghiên cứu vấn đề dẫn tới phát loại vi khuẩn yếm khí có khả phân hủy chất hữu Tiếp theo nghiên cứu điều kiện mà vi khuẩn yếm khí phát triển sinh trưởng tốt [7] Tới năm 1970, vấn đề địa trị Trung Cận Đông, khủng hoảng lượng kéo theo gia tăng đột biến giá dầu mỏ tiềm năng lượng từ biogas - nguồn lượng tái tạo bắt đầu ý nên chất trình XLYK kỹ thuật XLYK trở thành đề tài nóng giới nghiên cứu Hơn nữa, năm 1970 thời điểm mà Luật môi trường thông qua hàng loạt nước phát triển, động lực để phát triển công nghệ XLYK, với nguồn nước thải giàu hữu Kết nỗ lực hệ thống hóa, lượng hóa chế phân hủy chất hữu vốn phức tạp thành sơ đồ Hình 1.1 [8] Theo Hình 1.1, tác dụng tổ hợp vi sinh vật khác chất hữu phân hủy qua bốn cơng đoạn chuyển hóa: (1) thủy phân, (2) axit hóa, (3) axetat hóa và(4) metan hóa Theo Gujer vàZehnder (1983) [8], 70% metan hình thành từ chất trung gian axetat theo phản ứng (a) 30% metan cịn lại sản phẩm chuyển hóa chất trung gian quan trọng hyđrô theo phản ứng (b) CH3COOH  CH4 + CO2 (a) 4H2 + CO2  CH4 + 2H2O (b) Theo số tác giả khác tỷ lệ có chút khác biệt, cụ thể 72/28 [9] song chất giai đoạn không thay đổi Các phản ứng trung gian khác nhiệt động học q trình yếm khí tham khảo chi tiết [10] Hình 1 Sơ đồ chuyển hóa chất thải hữu điều kiện yếm khí Ngay từ 1982 Mc Carty [11] liệt kê bước tiến hóa cơng nghệ Bắt đầu hệ thống xử lý suất cao “bồn phản ứng khuấy trộn hoàn tồn” bồn phản ứng có kiểm sốt nhiệt độ thường sử dụng cho “nước” thải có nồng độ hữu rắn cao (thực chất bùn thải chất thải hữu cơ), cuối hệ phản ứng có suất cao nhờ mật độ cao vi sinh hoạt tính So với q trình hiếu khí q trình yếm khí tiêu tốn lượng hơn, sinh bùn thải hơn, chi phíxử lý thấp hơn, lại có tiềm thu hồi lượng dạng lượng tái tạo Trong điều kiện yếm khí, chất hữu chuyển đổi thành khíCH4, kg COD loại bỏ sinh khoảng 13,5 MJ CH4 lượng, đem lại 1,5 kWh điện (với hiệu suất chuyển đổi điện đạt 40%) [4] Trong lĩnh vực công nghệ vi sinh xử lý nước thải, động học phản ứng phân hủy yếm khí nói riêng phản ứng vi sinh nói chung tn theo phương trì nh Monod [12]: 𝑟𝑠𝑢 = 𝑘.𝑆 𝐾𝑆 +𝑆 𝑋 (1.1) Trong đó: rsu = tốc độ tiêu thụ chất (chất ô nhiễm), g/m /ngày k = tốc độ tiêu thụ chất riêng tối đa, g chất/g VSS/ngày X = nồng độ sinh khối (vi khuẩn), g VSS/m S = nồng độ chất (quyết định tốc độ), g/m3 Ks = số bán bão hòa = nồng độ chất mà tốc độ tiêu thụ chất riêng 50% tốc độ tiêu thụ chất riêng tối đa, g/m3 Theo phương trình (1.1), tốc độ phản ứng (xử lý) tỷ lệ thuận với tốc độ riêng tối đa k (trong trình XLYK hoạt tính tạo metan riêng - SMA (Specific Methanogenic Activity)) nồng độ vi sinh X hệ phản ứng SMA đặc trưng cho tổ hợp vi khuẩn yếm khí xúc tác cho hệ phản ứng, đại lượng đặc thù phụ thuộc vào nhiều yếu tố [10]: nhiệt độ; khả sinh hủy chất bẩn; có mặt chất độc hay chất ức chế; nồng độ phần VSS nước vào; mức độ axit hóa trước nước thải Như tốc độ phản ứng tỷ lệ thuận với mật độ vi sinh X, nghĩa muốn giảm thể tích thiết bị phải tăng X Tuy nhiên, tăng X cần hệ phản ứng cho phép làm việc X lớn mà khơng ảnh hưởng đến q trình vận hành (bùn bơng lắng có mật độ 6-8 g/L) Vì vậy, kỹ thuật phản ứng muốn có suất xử lý cao phải đồng thời giải u cầu: tăng X khơng thất vi sinh tiếp xúc vi sinh/nước thải tốt (khuấy trộn) Các thiết bị phản ứng “cao tải” (high - rate) hệ (AC, AF, UASB, FB, EGSB, IC, ) trì nh bày sau giải vấn đề 1.2 Lớp đệm bùn vi sinh yếm khí vi sinh dạng hạt Đột phá công nghệ XLYK nước thải giàu hữu phát minh nhóm G Lettinga [13] Theo ý tưởng tác giả, khả lưu giữ bùn hoạt tính bồn phản ứng đảm bảo nhờ hình thành lớp bùn bơng mật độ cao dày vài mét bùn dạng hạt có tốc độ lắng cao kết hợp với việc áp dụng cấu tách ba pha K/L/R có dạng hình nón phần lớp nước Theo G Lettinga [14], ý tưởng UASB kế thừa phát triển nghiên cứu hệ lọc yếm khí với tải lượng 10 kg COD/m3/ngày Trong trình nghiên cứu tác giả ghi nhận xuất bùn hoạt tính cóhoạt tính cao lắng tốt Điều dẫn tới cách tiếp cận kỹ thuật XLYK mà gọi UASB Năm 1971-1972 pilot quy mô6 và30 m3 thử nghiệm, kết thu chứng minh tính khả thi UASB Năm 1976 hệ UASB quy mô 200 m3 đưa vào vận hành nhà máy đường, năm sau hệ UASB quy mô 1000 m3 xây dựng sau nhà máy thứ hai Breda – Hà Lan Thành công UASB kéo theo hàng loạt nghiên cứu vi sinh dạng hạt Sự hình thành vi sinh dạng hạt Trong công nghệ XLYK nước thải, nước thải giàu hữu cơ, tế bào tự tập hợp lại để tạo thành hạt bùn bền, có hoạt tính cao có khả lắng nhanh trở thành đối tượng nghiên cứu hàng đầu [15] Hạt bùn vi sinh yếm khí tập hợp vi sinh có mật độ cao hình thành mối quan hệ tương hỗ với môi trường xung quanh (đặc trưng nước thải, đặc trưng thủy động dịng chảy) Sự phân hủy thành phần nhiễm nước thải thực nhờ lớp vi sinh chức tương ứng có mặt hạt bùn Trên sở quan hệ tương hỗ Fitzpatrick vàcs (1989) [16] cho hình thành bùn dạng hạt vi sinh, (1) tập hợp tế bào (tạo hạt bùn) để hình thành quần thể vi sinh không đồng dạng tập hợp đa tế bào điều kiện sinh lý thuận lợi, điều dẫn tới (2) tương tác tương hỗ vi sinh vật lân cận tạo thuận lợi cho trao đổi yếu tố di truyền, (3) tạo hạt bảo vệ vi sinh khỏi nguyên sinh động vật ăn vi khuẩn yếm khí, (4) thể tích hạt quãng đường khuếch tán sản phẩm lên men trung gian tới vi khuẩn chuyển hóa rút ngắn Điều làm giảm lượng tiêu tốn cho trình vận chuyển chất, tăng hiệu chuyển hóa, (5) điều kiện khơng tốt cho sinh trưởng vi sinh (ví dụ, pH vượt ngưỡng, có mặt chất độc v.v…), mơi trường cục hạt vi sinh bị ảnh hưởng hơn, chuyển hóa bị ảnh hưởng Các lý thuyết/quan điểm hình thành bùn vi sinh yếm khí dạng hạt liên quan đến hoạt động hệ XLYK cao tải, số quan điểm tạo hạt bùn tóm tắt Đối với q trình metan hóa, tự bám dính vi khuẩn hay tượng tự tập hợp vi khuẩn coi bắt đầu q trình tạo hạt bùn Nó hệ tương tác vi khuẩn-vi khuẩn vi khuẩn-bề mặt rắn Về mặt nhiệt động học, vi khuẩn tiến lại gần chúng đồng thời chịu lực tác động ngược chiều nhau, mặt lực đẩy tĩnh điện hạt dấu, mặt khác lực hút phân tử van de Waals, tương tác với phân tử nước (hydrat hóa) hạt hydrat hóa Một số lý thuyết/quan điểm, thuyết hóa lý, xây dựng sở phân tích tượng Lettinga et al (1980) [13] cho có mặt “mầm” hạt vật liệu rắn nhỏ, vi khuẩn bám lên bề mặt “mầm” để tạo thành phôi hạt Tiếp theo, vi khuẩn tiếp tục bám lên phôi hạt tạo hạt bùn thấy mắt Minh chứng cho vai trò mầm rắn chứng minh thực nghiệm Hulshoff Pol cho thêm hạt khoáng chất zeolit hydro-anthraxit có kích cỡ khoảng 100 m vào mơi trường ni vi khuẩn, chúng đẩy nhanh trình tạo hạt bùn yếm khí[17] Hulshoff Pol vàcs (1988) [18] cho tập hợp vi khuẩn thành hạt phản ứng vi khuẩn để chống lại áp lực trình chọn lọc Trong bồn/bể phản ứng kiểu UASB, dịng chất lỏng chảy ngược áp lực q trình chọn lọc, tác dụng dịng chảy bơng tế bào cókhả tạo hạt/bám dính vànhẹ bị rửa trôi theo nước khỏi bồn phản ứng, bồn cịn lại tập hợp vi sinh có khả bám dính có mật độ lớn hơn, nặng Áp lực chọn lọc dòng chảy quan trọng [15] Thực tế cho thấy bùn bơng yếm khí dần chuyển thành tập hợp dạng hạt thời gian lưu thủy lực (HRT) nhỏ Trong trường hợp HRT lớn hay áp lực chọn lọc nhỏ không thấy tượng tạo bùn hạt [17] [18] Như vậy, UASB hệ cao tải khác EGSB IC dịng nước thải ln chảy từ lên, ngược dòng cặn lắng áp lực chọn lọc trình tạo bùn dạng hạt Pereboom Vereijken (1994) [19] cho xuất phát điểm trình tạo hạt q trình va chạm, bám dính hạt chất rắn lơ lửng (SS) có nước thải đầu vào Chuyển động hỗn loạn có định hướng dịng chất lỏng chất khí bồn nguyên nhân va chạm Sự tăng kích thước hạt ban đầu túy khả bám dính vi khuẩn Đã ghi nhận nước thải có hàm lượng SS cao tạo thành nhiều hạt bùn nhỏ với phân bố hẹp, cịn SS thấp khơng có SS hạt bùn sinh có phân bố theo kích thước Ở pH trung tính vi khuẩn hạt tích điện âm Vìvậy, để tăng cường khả bám dính cần cấp thêm ion hóa trị dương cao (cation đa điện tích) cho hạt mầm Nhiều thực nghiệm chứng minh bổ sung ion Ca 2+, Mg2+ Al3+ tăng tốc trình tạo hạt bồn UASB [20] [21] Tuy nhiên, nồng độ Ca2+  500 mg/L có hại trình tạo hạt [21] Các cation đa điện tích tăng cường q trình tạo hạt nhờ tương tác chúng với lớp polyme ngoại bào (Extracellular Polymers-ECPs) vi khuẩn tiết [22] Ca2+ tạo liên kết bắc cầu ECPs với ECPs, vi khuẩn lại bám dính theo ECPs để tạo cấu trúc ba chiều từ tập hợp đa tế bào, vi khuẩn tiếp tục phát triển để tiếp tục tăng kích thước hạt Bản thân ECPs thúc đẩy đồng thời q trình bám dính (adhesion) vi khuẩn bên lực hút nội hạt (cohesion) đem lại độ bền cho hạt bùn ECP làm giảm điện tích âm bề mặt vi khuẩn, lực đẩy tĩnh điện giảm làm chúng hút lực van de Waals tăng cường khả bám dính vi khuẩn lên vật liệu mầm hạt [23] Show & cs (2004a,b) [24] [25] vàWang & cs (2004) [26] nghiên cứu ảnh hưởng chất keo tụ polyme lên trình tạo hạt bùn quátrì nh khởi động UASB cấp độ phịng thí nghiệm Các kết cho thấy với liều polyme mức 20 mg/L thời gian khởi động giảm 46%, đồng thời tải hữu đầu vào tăng thêm 33% Một yếu tố cản trở khả tự tập hợp vi khuẩn lực đẩy tế bào hydrat hóa (bọc phân tử nước) Để giảm điều cần giảm yếu tố hydrat hóa hay dehydrat hóa hay hydrophob hóa (làm kị nước) tế bào Mơ hình kỵ nước hóa tế bào đề xuất J Wilschut vàD Hoekstra [27], theo bề mặt vi khuẩn trở nên kị nước chúng hút để tạo hạt Về mặt nhiệt động, tăng tính kị nước làm giảm lượng tự Gibbs bề mặt tế bào, điều hỗ trợ tương tác tế bào-tế bào tăng cường trình tự tạo hạt vi khuẩn nước [28] Yếu điểm mơ hình/lý thuyết hóa lý coi tế bào 10 hạt vật chất túy Để khắc phục nhược điểm cần xét đến yếu tố sinh học Hướng nghiên cứu thúc đẩy mạnh hai thập kỷ gần Trong hệ xử lý nước thải vi khuẩn thường tập hợp nhiều loài cộng sinh có tác dụng tương hỗ phù hợp với q trình chuyển hóa sinh hóa mà chúng xúc tác để thực Mơ hình tác động tương hỗ cho tác động tương hỗ ban đầu dẫn tới hình thành tập hợp tế bào đầu tiên, vi khuẩn lạc vi hạt bao gồm vi khuẩn nhiều chủng với chức tương hỗ nhau, tăng kích thước ba chiều nhờ sinh trưởng Điều giải thích hoạt tính cao bùn dạng hạt tế bào khơng cần thiết cho chuỗi q trình chuyển hóa tự bị loại bỏ Để trì cấu trúc tối ưu này, tế bào thành phần hạt tự tập hợp theo trật tự hợp lý, nhờ chế thông tin di truyền chúng phát triển để tạo hạt [29] MacLeod & cs (1990) [30] vàGuiot & cs (1988) [31] tác giả đề xuất mô hình đa lớp (multi-layer) Theo đó, hạt bùn có cấu trúc đa lớp, lớp có thành phần tế bào chức chuyển hóa tương ứng Lớp chủ yếu vi khuẩn metan hóa, chúng đồng thời đóng vai trị vi hạt ban đầu để từ phát triển thành hạt bùn Các vi khuẩn tạo H2 tiêu thụ H2 nằm lớp hạt, lớp hỗn hợp vi khuẩn dạng que, dạng cầu dạng sợi Để tạo metan, thứ tự không gian lớp quan trọng Mơ hình cấu trúc đa lớp nhiều tác giả nghiên cứu Dubourgier & cs (1987) [32] kết hợp quan sát hiển vi kết hợp với phép đo hoạt tính đề xuất mơ hình tự tạo hạt Theo bắt đầu vi khuẩn Methanosaeta dạng sợi theo khuẩn lạc vi khuẩn axit hóa dạng cầu que để tạo hạt kích thước 10–50 m Tiếp theo, sợi Methanosaeta nhờ cấu trúc dạng sợi có bề mặt phát triển đóng vai trị cầu nối kéo hạt nhỏ lại với để tạo hạt lớn có kích thước từ 200 m trở lên Tiếp theo chúng dính tiếp vi khuẩn axit hóa vi khuẩn chức phù hợp khác Ý tưởng tác giả sợi Methanosaeta đóng vai trị định việc tạo độ bền cho hạt bùn cách tạo mạng lưới cấu trúc ổn định Các cơng trình sau Morgan &cs (1991a,b) [33] [34] củng cố giả thiết vai trò định Methanosaeta polyme ngoại bào việc tạo hạt bùn 128 QA (lít/phút) 1,5 2,5 3,5 4,5 QA/S (m/h) 291 437 582 728 874 1019 1165 1311 1456 Ktn 6,80 5,60 4,90 4,16 3,60 3,26 2,90 2,62 2,40 K tí nh 6,70 5,31 4,44 3,84 3,39 3,05 2,77 2,54 2,34 -1,53 -5,19 -9,33 -7,68 -5,73 -6,48 -4,60 -3,29 -2,42 2050 3500 4800 5700 6200 7000 7400 7700 8300 Ktn 2,05 2,33 2,40 2,28 2,07 2,00 1,85 1,71 1,66 K tí nh 2,55 2,40 2,25 2,11 1,98 1,86 1,75 1,66 1,57 24,31 2,99 -6,07 -7,35 -4,11 -6,88 -5,16 -3,17 -5,53 2800 4400 5650 6450 7500 8000 8350 8750 9100 Ktn 2,80 2,93 2,83 2,58 2,50 2,29 2,09 1,94 1,82 K tí nh 3,08 2,84 2,62 2,42 2,24 2,08 1,95 1,83 1,72 % sai số 9,99 -3,14 -7,39 -6,39 -10,46 -8,85 -6,73 -6,08 -5,56 5000 5400 6600 7600 8200 8900 9200 9600 9950 Ktn 5,00 3,60 3,30 3,04 2,73 2,54 2,30 2,13 1,99 K tí nh 3,84 3,42 3,07 2,79 2,54 2,34 2,17 2,02 1,89 -23,22 -4,90 -6,85 -8,39 -6,90 -7,90 -5,68 -5,31 -5,04 5000 6800 7950 8400 9200 9800 10400 10500 10800 % sai số 1850 500 % sai số 1850 1850 400 300 % sai số 1850 200 129 QA (lít/phút) 1,5 2,5 3,5 4,5 QA/S (m/h) 291 437 582 728 874 1019 1165 1311 1456 Ktn K tí nh % sai số 1850 Ktn K tí nh % sai số 100 5,0 4,5 4,0 3,4 3,1 2,8 2,6 2,3 2,2 4,94 4,19 3,65 3,23 2,91 2,64 2,42 2,24 2,08 -1,17 -7,52 -8,28 -3,83 -5,25 -5,63 -6,76 -3,98 -3,58 6600 7800 9200 10000 10600 11200 11500 11600 12000 6,6 5,2 4,6 4,0 3,5 3,2 2,9 2,6 2,4 6,49 5,17 4,34 3,76 3,32 2,99 2,71 2,49 2,30 -1,64 -0,52 -5,65 -6,05 -5,92 -6,69 -5,62 -3,47 -4,22 130 PHỤ LỤC MỐI LIÊN HỆ GIỮA ĐỘ NHỚT CỦA DUNG DỊCH ĐƯỜNG SUCROSE VÀ CỦA NƯỚC Ở 20oC 131 Bảng 19 Mối liên hệ độ nhớt dung dịch đường nước 200C [90] Sucrose Sucrose % a 1,000 Sucrose % a % A 1,026 26 2,568 51 17,20 1,053 27 2,703 52 19,30 1,082 28 2,819 53 21,75 1,112 29 3,009 54 24,63 1,144 30 3,181 55 28,02 1,177 31 3,369 56 32,06 1,213 32 3,574 57 36, 88 1,251 33 3,798 58 42,69 1,291 34 4,044 59 49,74 10 1,333 35 4,314 60 58,37 11 1,378 36 4,612 61 69,03 12 1,426 37 4,940 62 82,26 13 1,477 38 5,304 63 98,83 14 1,531 39 5,707 64 119,9 15 1,589 40 6,15 65 146,9 16 1,650 41 6,66 66 181,7 17 1,716 42 7,22 67 227,3 13 1,7S6 43 7,85 68 2S7,9 19 1,861 44 8,56 69 369,4 20 1,941 45 9,36 70 480,6 21 2,027 46 10,28 71 635,0 22 2,120 47 11,31 72 853,2 23 2,220 43 12,49 73 1163, 24 2,326 49 13,84 74 1628, 25 2,442 50 15,40 75 2323, 132 PHỤ LỤC CHƯƠNG TRÌNH XỬ LÝ SỐ LIỆU THỰC NGHIỆM TÌM MƠ HÌNH TỐN BẰNG NGƠN NGỮ LẬP TRÌNH FORTRAN 133 CHƯƠNG TRÌNH XỬ LÝ SỐ LIỆU THỰC NGHIỆM TÌM MƠ HÌNH TỐN BẰNG NGƠN NGỮ LẬP TRÌNH FORTRAN CHARACTER H*20 C C C C C C C C C C C C C C C C CHUONG TRINH haiRegHaiTM.FOR (2018) MO HINH THUC NGHIEM THONG KE BAC CUA YEU TO K=f(Vk,H1,H2) - Vk van toc - H1 muc ngap nuoc - H2 chieu cao xa nuoc Luan an cua Tran Manh Hai - IET VAST Neu IN=0 khong in cac ma tran so lieu vao va ket qua trung gian IN=1 thi in DIMENSION X1(300),X2(300),X3(300),Y(300),YTT(300),A(400),B(20), *TX(20,20),TY(20),DEL(300),SS(300) WRITE(*,*) 'FILE SO LIEU?' READ(*,1) H FORMAT(A) OPEN(1,FILE=H,STATUS='UNKNOWN') WRITE(*,*) 'FILE KET QUA?' READ(*,1) H OPEN(3,FILE=H,STATUS='UNKNOWN') READ (1,*) N,M,IN READ (1,*) (X1(I),X2(I),X3(I),Y(I),I=1,N) WRITE(3,*) 'SO LIEU BAN DAU :' WRITE(3,*) '****************' WRITE (3,2) N,M,IN FORMAT(/5X,'N=',I3,3X,'M=',I2,3X,'IN=',I1) DO 11 I=1,N X1(I)=X1(I)/1000 X2(I)=X2(I)/100 11 X3(I)=X3(I)/100 C C C Xem co in cac ma tran so lieu dau vao khong? IF(IN.EQ.0) GOTO 12 WRITE (3,*) WRITE (3,*)' I,X1(I),X2(I),X3(I) & Y(I) :' WRITE (3,3) (I,X1(I),X2(I),X3(I),Y(I),I=1,N) FORMAT(I5,4F11.4) 12 CONTINUE WRITE(3,*) WRITE(3,*) 'KET QUA TINH TOAN:' WRITE(3,*) '******************' WRITE(3,*) SUM=0 DO 13 I=1,N 13 SUM=SUM+Y(I) YTB=SUM/N SUM=0 DO 14 I=1,N 134 14 SUM=SUM+(Y(I)-YTB)**2 TSS=SUM DO I=1,M DO J=1,M TX(I,J)=0 TY(I)=0 DO 20 I=1,M DO 20 K=1,N CALL TINHC(C,X1,X2,X3,I,K) GOTO (101,102,103,104,105,106,107,108,109,110),I 101 TX(I,1)=TX(I,1)+C*1 102 TX(I,2)=TX(I,2)+C*X1(K) 103 TX(I,3)=TX(I,3)+C*X2(K) 104 TX(I,4)=TX(I,4)+C*X3(K) 105 TX(I,5)=TX(I,5)+C*X1(K)*X2(K) 106 TX(I,6)=TX(I,6)+C*X1(K)*X3(K) 107 TX(I,7)=TX(I,7)+C*X2(K)*X3(K) 108 TX(I,8)=TX(I,8)+C*X1(K)**2 109 TX(I,9)=TX(I,9)+C*X2(K)**2 110 TX(I,10)=TX(I,10)+C*X3(K)**2 20 CONTINUE DO 30 I=1,M DO 30 K=1,N CALL TINHC(C,X1,X2,X3,I,K) 30 TY(I)=TY(I)+C*Y(K) DO 40 I=2,M I1=I-1 DO 40 J=1,I1 40 TX(I,J)=TX(J,I) C C Xem xet co in ket qua trung gian khong? C IF(IN.EQ.0) GOTO 44 WRITE(3,41) ((TX(I,J),J=1,10),TY(I),I=1,10) 41 FORMAT(7E12.6) 44 CONTINUE DO 60 I=1,M DO 50 J=1,M IJ=(I-1)*M+J 50 A(IJ)=TX(I,J) 60 B(I)=TY(I) CALL SIMQ(A,B,M,KS) IF(KS.EQ.0) GOTO 70 WRITE(3,66) KS 66 FORMAT(3X,'KS=',I2,'- VO NGHIEM') GOTO 90 70 WRITE(3,77) (I,B(I),I=1,M) 77 FORMAT(5X,'I=',I2,5X,'B(I)=',E14.6) WRITE(3,*) WRITE(3,*) 'I,X1(I),X2(I),X3(I),Y(I),YTT(I),DEL(I),SS(I):' WRITE(3,*) SUM=0 DO 80 I=1,N YTT(I)=B(1)+B(2)*X1(I)+B(3)*X2(I)+B(4)*X3(I)+B(5)*X1(I)*X2(I)+ *B(6)*X1(I)*X3(I)+B(7)*X2(I)*X3(I)+B(8)*X1(I)**2+B(9)*X2(I)**2+ *B(10)*X3(I)**2 DEL(I)=Y(I)-YTT(I) SS(I)=ABS(DEL(I)/Y(I)*100.) SUM=SUM+DEL(I)**2 80 CONTINUE RSS=SUM 135 R2=1.-RSS/TSS SUM=0 DO 81 I=1,N 81 SUM=SUM+SS(I) SSTB=SUM/N 85 WRITE(3,88) (I,X1(I),X2(I),X3(I),Y(I),YTT(I),DEL(I),SS(I),I=1,N) 88 FORMAT(2X,I3,3X,3F10.4,2X,4F9.4) WRITE(3,89) SSTB 89 FORMAT(/3X,'SSTB=',F6.2) WRITE(3,899) R2 899 FORMAT(/5X,'R2=',F8.5) 90 CONTINUE STOP END SUBROUTINE TINHC(C,X1,X2,X3,I,K) C C TINH MA TRAN XX C DIMENSION X1(300),X2(300),X3(300) GOTO (1,2,3,4,5,6,7,8,9,10),I C=1 RETURN C=X1(K) RETURN C=X2(K) RETURN C=X3(K) RETURN C=X1(K)*X2(K) RETURN C=X1(K)*X3(K) RETURN C=X2(K)*X3(K) RETURN C=X1(K)**2 RETURN C=X2(K)**2 RETURN 10 C=X3(K)**2 RETURN END SUBROUTINE SIMQ(A,B,N,KS) C C GIAI HE PHUONG TRINH CHUAN C DIMENSION A(400),B(20) TOL=0 KS=0 JJ=-N DO 65 J=1,N JY=J+1 JJ=JJ+N+1 BIGA=0 IT=JJ-J DO 30 I=J,N IJ=IT+I IF(ABS(BIGA)-ABS(A(IJ))) 20,30,30 20 BIGA=A(IJ) IMAX=I 30 CONTINUE IF(ABS(BIGA)-TOL) 35,35,40 136 35 KS=1 RETURN 40 II=J+N*(J-2) IT=IMAX-J DO 50 K=J,N II=II+N I2=II+IT SAVE=A(II) A(II)=A(I2) A(I2)=SAVE 50 A(II)=A(II)/BIGA SAVE=B(IMAX) B(IMAX)=B(J) B(J)=SAVE/BIGA IF(J-N) 55,70,55 55 IQS=N*(J-1) DO 65 IX=JY,N IXJ=IQS+IX IT=J-IX DO 60 JX=JY,N IXJX=N*(JX-1)+IX JJX=IXJX+IT 60 A(IXJX)=A(IXJX)-(A(IXJ)*A(JJX)) 65 B(IX)=B(IX)-(B(J)*A(IXJ)) 70 NY=N-1 IT=N*N DO 80 J=1,NY IA=IT-J IB=N-J IC=N DO 80 K=1,J B(IB)=B(IB)-A(IA)*B(IC) IA=IA-N 80 IC=IC-1 RETURN END KẾT QUẢ TÍNH TỐN SO LIEU BAN DAU : **************** N=253 M=10 IN=0 KET QUA TINH TOAN: ****************** I= I= I= I= I= I= I= I= I= I=10 B(I)= 0.788862E+01 B(I)= -0.553385E+01 B(I)= 0.400454E+00 B(I)= -0.135973E+02 B(I)= -0.604282E+00 B(I)= 0.647842E+01 B(I)= 0.710867E+00 B(I)= 0.161286E+01 B(I)= 0.801324E-01 B(I)= 0.378613E+01 137 Bảng 20 Kết tính KMH mức sai số so với Ktn 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 X1(I) (km/h) 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 X2(I) (m) 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 2.85 X3(I) (m) 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 Y(I) 3.70 3.47 3.25 2.80 2.67 2.49 2.25 2.13 4.45 4.27 3.78 3.26 3.00 2.80 2.50 2.27 2.16 5.00 4.47 4.05 3.52 3.20 2.91 2.65 2.40 2.26 6.25 5.20 4.45 3.84 3.37 3.09 2.78 2.60 2.38 7.20 5.87 5.10 4.24 3.73 3.40 YTT(I) 3.81 3.39 3.05 2.77 2.56 2.42 2.35 2.35 4.44 3.93 3.49 3.12 2.81 2.58 2.41 2.31 2.29 5.14 4.54 4.00 3.54 3.14 2.81 2.55 2.36 2.23 5.92 5.22 4.60 4.03 3.54 3.12 2.76 2.48 2.26 6.78 5.98 5.26 4.61 4.02 3.50 DEL(I) -0.11 0.08 0.20 0.03 0.11 0.07 -0.10 -0.22 0.01 0.34 0.29 0.14 0.19 0.22 0.09 -0.04 -0.13 -0.14 -0.07 0.05 -0.02 0.06 0.10 0.10 0.04 0.03 0.33 -0.02 -0.15 -0.19 -0.17 -0.03 0.02 0.12 0.12 0.42 -0.11 -0.16 -0.37 -0.29 -0.10 SS(I) 2.96 2.18 6.15 0.99 4.01 2.68 4.49 10.22 0.29 8.02 7.68 4.39 6.23 7.88 3.49 1.97 5.81 2.82 1.49 1.12 0.49 1.91 3.39 3.77 1.77 1.14 5.28 0.41 3.27 5.04 5.06 0.94 0.61 4.77 5.09 5.90 1.90 3.17 8.62 7.72 3.01 138 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 1.165 1.311 1.456 0.437 0.582 0.728 0.874 1.019 1.165 1.311 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 2.85 2.85 2.85 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 2.65 0.10 0.10 0.10 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 3.05 2.76 2.54 2.93 3.00 2.64 2.53 2.37 2.18 1.98 3.60 3.35 2.92 2.77 2.63 2.33 2.16 2.00 4.70 4.27 3.85 3.40 3.00 2.80 2.50 2.36 2.14 5.90 5.00 4.30 3.68 3.27 2.97 2.70 2.51 2.34 7.20 5.87 4.90 4.20 3.70 3.31 3.05 2.80 2.56 3.05 2.67 2.36 3.21 2.88 2.62 2.43 2.31 2.25 2.27 3.76 3.33 2.98 2.69 2.48 2.33 2.25 2.24 4.97 4.38 3.86 3.41 3.03 2.72 2.48 2.31 2.20 5.76 5.08 4.47 3.92 3.45 3.05 2.71 2.44 2.24 6.63 5.85 5.15 4.51 3.94 3.44 3.01 2.65 2.35 0.00 0.09 0.18 -0.28 0.12 0.02 0.10 0.06 -0.07 -0.29 -0.16 0.02 -0.06 0.08 0.15 0.00 -0.09 -0.24 -0.27 -0.11 -0.01 -0.01 -0.03 0.08 0.02 0.06 -0.06 0.14 -0.08 -0.17 -0.24 -0.18 -0.08 -0.01 0.07 0.10 0.57 0.02 -0.25 -0.31 -0.24 -0.13 0.04 0.15 0.21 0.06 3.25 7.12 9.48 3.96 0.72 3.98 2.66 3.33 14.48 4.31 0.46 2.05 2.76 5.80 0.07 4.08 11.81 5.64 2.54 0.36 0.43 1.13 2.73 0.80 2.33 2.77 2.40 1.55 3.93 6.66 5.49 2.54 0.27 2.88 4.35 7.95 0.31 5.10 7.41 6.52 4.02 1.30 5.49 8.10 139 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 0.437 0.582 0.728 0.874 1.019 1.165 1.311 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.437 0.582 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.45 2.25 2.25 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.50 0.50 2.67 2.80 2.50 2.37 2.20 2.03 1.89 3.65 3.60 3.23 3.04 2.77 2.57 2.33 2.13 2.00 4.60 4.13 3.80 3.40 3.03 2.80 2.48 2.31 2.16 5.60 4.87 4.30 3.68 3.30 3.00 2.68 2.47 2.29 7.20 5.93 4.80 4.04 3.60 3.23 3.00 2.69 2.48 2.73 2.70 3.03 2.72 2.48 2.30 2.20 2.16 2.19 4.06 3.59 3.19 2.85 2.58 2.38 2.25 2.19 2.19 4.80 4.23 3.73 3.30 2.93 2.64 2.42 2.26 2.17 5.60 4.94 4.35 3.82 3.37 2.98 2.66 2.41 2.22 6.49 5.73 5.04 4.42 3.87 3.39 2.98 2.63 2.35 2.85 2.56 -0.36 0.08 0.02 0.07 0.00 -0.13 -0.30 -0.41 0.01 0.04 0.19 0.19 0.19 0.08 -0.06 -0.19 -0.20 -0.10 0.07 0.10 0.10 0.16 0.06 0.05 -0.01 0.00 -0.07 -0.05 -0.14 -0.07 0.02 0.02 0.06 0.07 0.71 0.20 -0.24 -0.38 -0.27 -0.16 0.02 0.06 0.13 -0.12 0.14 13.40 2.91 0.96 2.87 0.12 6.43 15.99 11.34 0.29 1.35 6.28 6.84 7.32 3.41 2.71 9.68 4.26 2.34 1.85 3.00 3.14 5.63 2.57 2.21 0.50 0.06 1.44 1.14 3.88 1.97 0.72 0.83 2.59 2.89 9.90 3.39 5.08 9.48 7.52 4.96 0.83 2.28 5.14 4.54 5.09 140 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 0.728 0.874 1.019 1.165 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.437 0.582 0.728 0.874 1.019 1.165 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.05 2.05 2.05 2.05 2.05 2.05 0.50 0.50 0.50 0.50 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.50 0.50 0.50 0.50 0.50 0.50 2.44 2.23 2.11 1.95 3.27 3.10 2.80 2.53 2.37 2.18 2.00 1.92 4.20 4.00 3.55 3.20 2.87 2.69 2.45 2.27 2.16 5.40 4.67 4.10 3.64 3.23 3.00 2.70 2.49 2.28 7.10 5.57 4.80 4.16 3.73 3.33 3.00 2.71 2.52 2.60 2.60 2.40 2.23 2.14 1.95 2.34 2.18 2.09 2.07 3.43 3.04 2.72 2.47 2.29 2.18 2.13 2.16 4.63 4.08 3.60 3.19 2.84 2.57 2.36 2.22 2.15 5.45 4.81 4.24 3.73 3.29 2.92 2.61 2.38 2.22 6.35 5.61 4.94 4.34 3.81 3.34 2.95 2.62 2.36 2.69 2.41 2.21 2.07 2.00 2.00 0.10 0.05 0.02 -0.12 -0.16 0.06 0.08 0.06 0.08 0.00 -0.13 -0.24 -0.43 -0.08 -0.05 0.01 0.03 0.12 0.09 0.05 0.01 -0.05 -0.14 -0.14 -0.09 -0.06 0.08 0.09 0.11 0.06 0.75 -0.04 -0.14 -0.18 -0.08 -0.01 0.05 0.09 0.16 -0.09 0.19 0.19 0.16 0.14 -0.05 4.20 2.19 0.76 6.40 4.89 1.79 2.69 2.22 3.26 0.05 6.70 12.36 10.31 2.04 1.47 0.38 0.96 4.56 3.71 2.23 0.52 1.02 2.98 3.31 2.38 1.76 2.74 3.16 4.40 2.81 10.53 0.75 3.02 4.36 2.06 0.41 1.79 3.42 6.39 3.34 7.20 8.11 7.32 6.68 2.33 141 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 1.311 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.437 0.582 0.728 0.874 1.019 1.165 0.437 0.582 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 2.05 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 0.50 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.50 0.50 0.50 0.50 0.50 0.50 0.40 0.40 1.80 3.30 3.23 3.00 2.72 2.50 2.40 2.20 2.04 1.88 4.30 4.00 3.75 3.20 2.93 2.71 2.48 2.31 2.12 5.40 4.60 4.10 3.52 3.20 2.97 2.65 2.40 2.24 6.80 5.60 4.90 4.16 3.60 3.26 2.90 2.62 2.40 2.33 2.40 2.28 2.07 2.00 1.85 2.93 2.83 2.06 3.72 3.28 2.91 2.61 2.37 2.21 2.11 2.09 2.13 4.48 3.94 3.48 3.08 2.76 2.50 2.31 2.19 2.13 5.31 4.68 4.13 3.64 3.21 2.86 2.58 2.36 2.21 6.22 5.50 4.85 4.27 3.75 3.30 2.92 2.61 2.37 2.53 2.27 2.08 1.96 1.91 1.92 3.13 2.78 -0.26 -0.42 -0.05 0.09 0.11 0.13 0.19 0.09 -0.05 -0.25 -0.18 0.06 0.27 0.12 0.17 0.21 0.17 0.12 -0.01 0.09 -0.08 -0.03 -0.12 -0.01 0.11 0.07 0.04 0.03 0.58 0.10 0.05 -0.11 -0.15 -0.04 -0.02 0.01 0.03 -0.20 0.13 0.20 0.11 0.09 -0.07 -0.20 0.05 14.59 12.60 1.45 3.04 4.16 5.05 7.91 3.91 2.28 13.16 4.10 1.43 7.18 3.62 5.92 7.79 6.93 5.36 0.62 1.62 1.83 0.69 3.33 0.47 3.59 2.71 1.60 1.14 8.46 1.76 0.98 2.55 4.14 1.34 0.82 0.28 1.18 8.41 5.44 8.79 5.38 4.68 3.92 6.83 1.77 142 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 0.728 0.874 1.019 1.165 1.311 1.456 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 0.291 0.437 0.582 0.728 0.874 1.019 1.165 1.311 1.456 SSTB= 3.95 R2= 0.97722 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 1.85 0.40 0.40 0.40 0.40 0.40 0.40 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.30 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.20 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 2.58 2.50 2.29 2.09 1.94 1.82 3.60 3.30 3.04 2.73 2.54 2.30 2.13 1.99 5.00 4.50 4.00 3.40 3.10 2.80 2.60 2.30 2.20 6.60 5.20 4.60 4.00 3.50 3.20 2.90 2.60 2.40 2.50 2.28 2.13 2.06 2.05 2.10 3.81 3.37 2.99 2.68 2.44 2.26 2.16 2.12 5.18 4.57 4.03 3.55 3.15 2.82 2.55 2.35 2.22 6.10 5.40 4.77 4.20 3.70 3.27 2.91 2.61 2.39 0.08 0.22 0.16 0.03 -0.11 -0.28 -0.21 -0.07 0.05 0.05 0.10 0.04 -0.03 -0.13 -0.18 -0.07 -0.03 -0.15 -0.05 -0.02 0.05 -0.05 -0.02 0.50 -0.20 -0.17 -0.20 -0.20 -0.07 -0.01 -0.01 0.01 3.28 8.81 6.82 1.66 5.44 15.60 5.84 2.00 1.74 1.95 4.06 1.58 1.38 6.72 3.53 1.47 0.69 4.54 1.60 0.56 2.00 2.13 0.88 7.53 3.80 3.60 4.93 5.66 2.19 0.27 0.55 0.38 ... thuật tuần hoàn nội IC sáng chế Vellinga (1986) [54] Thành phần cấu tạo hệ xử lý tuần hoàn nội gồm hệ phân bố nước đầu 19 vào, cấu tuần hoàn nội, vùng phản ứng cao tải phía vùng phản ứng tải lượng...