Độ thấm khí của một vật liệu có lỗ rỗng là một đại lượng đặc trưng cho khả năng cho phép chất khí đi xuyên qua mà không làm thay đổi cấu trúc của chất đó. Tính thấm của vật liệu phụ thuộc vào nhiều thông số như độ rỗng, áp suất khí thẩm thấu, hình dạng, độ nhám và tính kết nối giữa các lỗ rỗng. Đối với trường hợp của đá xi măng, trong cấu trúc lỗ rỗng luôn tồn tại hai pha lỏng và khí.
Tạp chí Khoa học Cơng nghệ Xây dựng NUCE 2018 12 (4): 135–146 ĐÁNH GIÁ THỰC NGHIỆM VÀ MÔ PHỎNG ĐỘ THẤM KHÍ CỦA ĐÁ XI MĂNG CĨ XÉT ĐẾN ẢNH HƯỞNG CỦA ĐỘ BÃO HÒA NƯỚC Từ Sỹ Quâna,∗ a Khoa Kỹ thuật Xây dựng, Trường Đại học Giao thông Vận tải, Số 3, phố Cầu Giấy, phường Láng thượng, quận Đống đa, Hà Nội, Việt Nam Lịch sử viết: Nhận ngày 28/03/2018, Sửa xong 10/5/2018, Chấp nhận đăng 30/5/2018 Tóm tắt Độ thấm khí vật liệu có lỗ rỗng đại lượng đặc trưng cho khả cho phép chất khí xuyên qua mà khơng làm thay đổi cấu trúc chất Tính thấm vật liệu phụ thuộc vào nhiều thơng số độ rỗng, áp suất khí thẩm thấu, hình dạng, độ nhám tính kết nối lỗ rỗng Đối với trường hợp đá xi măng, cấu trúc lỗ rỗng tồn hai pha lỏng khí Sự tồn nước tự hệ thống lỗ rỗng có ảnh hưởng khơng nhỏ đến khả truyền dẫn hấp thụ khí vật liệu Mơ hình thực nghiệm giới thiệu báo phát triển ứng dụng phòng thí nghiệm Xây dựng Cơ học (Laboratoire en Génie Civil et Génie Mécanique) thuộc trường Trung ương thành phố Nantes (Ecole Centrale de Nantes), Cộng hòa Pháp Trên sở phân tích kết đo độ thấm khí thu mẫu hình trụ, dựa lý thuyết mô xây dựng thiết lập, nghiên cứu cho phép thông số cần thiết tương ứng với loại vật liệu xem xét Từ khoá: định luật Darcy; độ thấm nội tại; hệ số thấm thủy lực; độ thấm khí; hiệu ứng Klinkenberg EVALUATION BY EXPERIMENT AND MODELING GAS PERMEABILITY OF CEMENT - BASED MATERIALS CONSIDERING THE EFFECT OF THE WATER SATURATION DEGREE Abstract The gas permeability is a property of the porous media that measures the capacity of a gas to pass through without altering the structure of the substance The permeability of the material is depended on many parameters such as porosity, osmotic pressure, shape, roughness and pore connectivity In the case of hardened cement, two phases, i.e liquid and gas, exist in its pore structure The existence of free water in the pore system has a great influence on the ability of the material to transmit and absorb gas The experimental model introduced in the paper has been developed and applied in the Laboratory of Construction and Mechanics Engineering (Laboratoire en Génie Civil et Génie Mécanique - Ecole Centrale de Nantes), France Based on the analysis of gas permeability measurements obtained on cylindrical samples and many theoretical simulations have been established and set up, the necessary parameters for each type of material considered are indicated Keywords: Darcy’s law; intrinsic permeability; hydraulic conductivity; gas permeability; Klinkenberg effect https://doi.org/10.31814/stce.nuce2018-12(4)-15 © 2018 Trường Đại học Xây dựng (NUCE) Giới thiệu chung Chất lỏng, cụ thể nước từ lâu sử dụng việc đo độ thấm vật liệu có nguồn gốc từ xi măng, vốn có giá trị tương đối nhỏ so với vật liệu rời rạc khác đất hay cát Darcy [1] ∗ Tác giả Địa e-mail: syquantu@gmail.com (Quân, T S.) 135 Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng đề xuất mơ hình thí nghiệm mà theo đánh giá độ thấm dựa khả thẩm thấu lan tỏa nước vào hệ thống rỗng vật liệu Với mẫu có kích thước lớn, thời gian đo kéo dài đến tháng, khiến cơng tác thí nghiệm trở nên tốn thời gian công sức Để khắc phục tình trạng trên, nhiều kỹ thuật đo thấm tiên tiến đề xuất thay thế, chẳng hạn sử dụng sóng siêu âm hay sử dụng khí trơ xun qua mẫu Mơ hình đo thấm trình bày báo dựa nguyên lý thay nước loại khí trơ khơng phản ứng hóa học với bê tơng Về bản, tính chất lý hóa nước khí trơ ln có khác biệt Do vậy, việc sử dụng khí trơ đòi hỏi phải có hiểu biết đặc điểm loại khí cách thức tương tác chúng với vật liệu xem xét 1.1 So sánh nguyên lý đo độ thấm khí độ thấm nước Chất khí thường ưu tiên sử dụng để đo độ thấm nội mẫu trạng thái đóng rắn Điều dựa trên khả thẩm thấu chất khí vào mơi trường rỗng thấp thường tốt chất lỏng Giả sử ta có mẫu hình trụ vật liệu đồng đẳng hướng, có đường kính D, chiều cao H chịu áp suất chất lưu đầu vào P1 , áp suất đầu P2 (Hình 1) Ở thang bậc cấu trúc vi mơ, chất lưu sử dụng thí nghiệm cần thỏa mãn điều kiện sau: ◦ Lực quán tính nhỏ so với lực dòng gây độ nhớt ◦ Dòng chảy thuộc chế độ chảy tầng ◦ Chất lưu qua khơng tạo phản ứng hóa học với mơi trường đo Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 Hình1 Mơhình hìnhhóa hóathí thínghiệm nghiệm đo đo độ độ thấm nước (thấm Hình : Mơ (thấm khí) khí) trongbài bàitốn tốn11chiều chiều Đối với chất lỏng, với giả thiết áp suất tác dụng biến thiên tuyến tính dọc theo chiều Đối với với giả áp mơi suấttrường tác dụng thiên tính thơng dọc theo cao mẫu, độ cao chất mẫu,lỏng, độ thấm nội thiết khibiến tuyến tính qua chiều định luật thấm nộiDarcy mơi trường tính thơng qua định luật Darcy [2] : [2]: µµE E QQ HH rρE EggSS PP -− P2P2 1 KKv v== (1) (1) µE làvới độ nhớt động chất lỏng (Pa.s); ρE khối lượng riêng chất lỏng (kg/m3 ); Q lưu lượng thấm qua mẫu (m3 /s); S diện tích thiết diện mẫu (m2 ); H chiều cao mẫu (m); µ E : Độ nhớt động 2chất lỏng (Pa.s) ; r E : Khối lượng riêng chất lỏng (kg/m3) ; g gia tốc trọng trường (m/s ); P1 áp suất đầu vào mẫu (m H2 O) P2 áp suất đầu : Lưu lượng thấm qua mẫu (m3/s) ; S : Diện tích tiết diện mẫu (m2) ; mẫu (m HQ2 O) nhớt túy Do TrongHthực tế, khí chất lưu; nén dòng khơng phải g : Giachảy tốc trọng trường (m/s ) : Chiều cao trơ củalà mẫu (m) đó, thay chất lỏng chất khí ta khơng thể áp dụng trực tiếp định luật Darcy Ở thang bậc vi P1 : Áp suất đầu vào mẫu (m H2O) ; P2 : Áp suất đầu mẫu (m H2O) 136 Trong thực tế, khí trơ chất lưu nén dòng chảy khơng phải nhớt túy Do đó, thay chất lỏng chất khí ta khơng thể áp dụng trực tiếp định luật Darcy Ở thang bậc vi mô, chế độ dừng, dòng chảy khí trơ kết hợp hệ: dịch chuyển nhớt dịch chuyển trượt Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng mơ, chế độ dừng, dòng chảy khí trơ kết hợp hệ: dịch chuyển nhớt dịch chuyển trượt Dịch chuyển trượt tạo tương tác phân tử khí thành lỗ rỗng Sự trượt trở nên quan trọng kích thước lỗ rỗng mơi trường nhỏ áp suất phun khí thấp Với áp suất khí phun thấp, ứng xử khí trơ xem khí lý tưởng [2]: Ka = Q 2µHP2 S P2 − P2 (2) Ka độ thấm tường chất khí (m2 ); µ độ nhớt động (Pa.s); P1 áp suất đầu vào mẫu (Pa) P2 áp suất đầu mẫu (Pa) Trong thực tế, độ thấm nội phụ thuộc chủ yếu vào cấu trúc lỗ rỗng độ ẩm môi trường Klinkenberg Carman đề xuất cách tính độ thấm nội từ độ thấm khí Tuy nhiên mơ hình Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 Klinkenberg đề xuất thường ưu tiên sử dụng thuận tiện cơng tác thí nghiệm Trong thực tế, độ thấm nội phụ thuộc chủ yếu vào cấu trúc lỗ rỗng độ ẩm môi trường.Klinkenberg Carman đềnội xuấttại cách tính độ thấm nội từ độ thấm 1.2 Hiệu ứng Klinkenberg cách tính độ thấm khí Tuy nhiên mơ hình Klinkenberg đề xuất thường ưu tiên sử dụng thuận tiệnKtrong cơng tác thí nghiệm Độ thấm nội v thực tế đại lượng phụ thuộc vào thân vật liệu Hiệu ứngđiKlinkenberg cách tính độ thấm nội thuộc vào tính chất1.2 khí qua Theo Klinkenberg, độtạithấm nội Kv hàm mà không phụ độ thấm tường Độsuất thấm nội Kbình Ka nghịch đảo áp trung [3]:tế đại lượng phụ thuộc vào thân vật liệu mà v thực Hệ số β phụ khơng phụ thuộc vào tính chất khí qua.Theo Klinkenberg, độ thấm nội K v β hàm độ thấm tường K a đảo áp suất trung bình [3]: Ka = Kv + nghịchvới β hệ số Klinkenberg (Pa) P m ỉ b ÷÷ với b : hệ s Klinkenberg (Pa) (3) K a = K V ỗỗ1 + m ø è Pmôi thuộc vào độ rỗng trường tính chất khí thấm qua Khi (3) áp suất trung β Hệ số b phụ thuộc vào độ rỗng mơi trường tính chất khí thấm qua Khi bình Pm tăng, tỉ số giảm 0, Ka tiến đến b giá trị Kv Độ thấm nội Kv phụ thuộc vào độ ẩm ápPsuất giảm 0, K a tiến đến giá trị K V Độ thấm nội m trung bình Pm tăng, tỉ số vật liệu làm thí nghiệm Giá trị Kv tăngPmkhi hàm lượng nước giảm Với giá trị áp suất trung bình K V phụ thuộc vào độ ẩm vật liệu làm thí nghiệm Giá trị K V tăng hàm lượng định sẵn, β thể đặc tính vật liệu trạng thái độ ẩm cho Lấy xấp xỉ tuyến tính nước giảm.Với giá trị áp suất trung bình định sẵn, b thể đặc tính vật liệu đối giá trị Ka theo áp Lấy suấtxấptrung bình Pmgiá , tatrịcó ướcnghịch lượng với mộtnghịch trạng tháiđảo độ ẩm cho xỉ tuyến tính thể đảođược độ thấm nội K a theo Kv Trong thực tế, áp điểm xáctế,của việc lấy xấp xỉ suất trung bìnhthuộc có thểđộ ướcchảy lượngtầng độ thấmđịnh nội tạiđộ Trong thực Pm , ta chế K Vchính điểm thuộc chế độ chảy tầng định độ xác việc lấy xấp xỉ K V Kv biểu diễn Hình (m2) biểu diễn Hình Kvß/Pm Kv Kv 1/Pm (Pa-1) Hình : Cách xác định gần độ thấm nội theo độ thấm tường Hình Cách xác định gần độ thấm nội theo độ thấm tường 137 Quân, T S / Tạp chí Khoa học Công nghệ Xây dựng Nghiên cứu thực nghiệm đo độ thấm khí in Civil Engineering NUCE 2018 13(5):1-16 Journal of Science and Technology 2.1 Sơ đồ thí nghiệm Nghiên cứu thực nghiệm đo độ thấm khí Thiết bị sử dụng 2.1 Sơ đồ thíphòng nghiệm thí nghiệm bao gồm khoang hình trụ loại Cembureau (tên riêng) sáng chế Kollek [4], mẫu phòng đo độthíthấm khí.baoSơgồm đồ hoạtkhoang độnghình trụ hệ loại thống Thiết bị sử chứa dụng nghiệm Cembureau (tên riêng) sáng chế Kollek [4], chứa mẫu đo độ thấm khí Sơ đồ miêu tả Hình 3: hoạt động hệ thống miêu tả Hình 3: 0-5bars Hình : Sơ đồ thí nghiệm đo độ thấm khí với áp suất phun xác định Hình Sơ đồ thí nghiệm đo độ thấm khí với áp suất phun xác định Theo Cassanas [5], thiết bị đo bao gồm dụng cụ sau: o Bình khí trơ (azote, argon, helium, hiđrogen ) Theo Cassanas [5], bịápđovớibao gồm o thiết Bộ điều áp suất phun cực đạidụng bars.cụ sau (Hình 4): o Khoang chứa mẫu loại Cembureau, cho phép thử với đường kính mẫu thay đổi ◦ Bình khí trơ (azote, argon, hiđrogen ) cm 10helium, cm o Thiết bị điện tử đo lưu lượng khí, cho phép kiểm sốt lưu lượng tối đa 1000 ◦ Bộ điều áp với ápml/min, suất phun bars với ápcực suất đại phun5cực đại 10 bars nhiệt độ 20°C o Bộ phận thu thập liệu đo (lưu lượng, áp suất) o Một tính cố định dùng ghi xửvới lý liệu kính mẫu thay đổi cm 10 cm ◦ Khoang chứa mẫu loạimáy Cembureau, cho để phép thử đường o Một áp kế để đo áp suất khí thời điểm tiến hành thí nghiệm kế đokhí, nhiệtcho độ phòng ◦ Thiết bị điện tửo đoMột lưunhiệt lượng phép kiểm soát lưu lượng tối đa 1000 ml/min, với áp suất ◦ phun cực đại 10 bars nhiệt độ 20 C Bộ phận thu thập liệu đo (lưu lượng, áp suất): ◦ Một máy tính cố định dùng để ghi xử lý liệu ◦ Một áp kế để đo áp suất khí thời điểm tiến hành thí nghiệm ◦ Một nhiệt kế đo nhiệt độ phòng Cụ thể chi tiết khoang chứa mẫu loại Cembureau biểu thị sơ đồ Hình Mẫu thí nghiệm hình trụ, mặt bên ơm sát lớp polime đàn hồi chống thấm, có đường kính đường kính mẫu thử (Hình 4(b)) Bộ phận có tác dụng ngăn cản khí trơ lọt qua thành bên làm thí nghiệm Lớp polime lồng vào buồng khí có dạng hình xăm xe đạp Khi áp suất buồng khí đạt 0,7 MPa, áp suất tác dụng lên mặt bên mẫu đủ lớn để coi thành bên mẫu kín khí Mặt mặt mẫu tiếp xúc với hai đĩa kim loại có rãnh chia Những rãnh chia có tác dụng đảm bảo phân bố đặn áp suất tiết diện mẫu trình phun (Hình 5) 138 Journal of of Science and Technology Civil Engineering 2018 13(5):1-16 Journal Science and in Civil Engineering NUCE 2018 13(5):1-16 Quân, T.Technology S / Tạp chíinKhoa học Cơng nghệ NUCE Xây dựng Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 a) b) (a) (b) Hình : a)thiết Các thiết bị đo áp suất,nhiệt nhiệt độ, lưu lượng đổi tínđổi hiệu; b)hiệu; Khoang Hình : a) Các suất, lưulượng lượng bộquy quyđổi Khoang Hình :4a) Các thiết bịbị đođo ápáp suất, nhiệt độ,độ,lưu vàvàbộquy tíntínhiệu; b)b)Khoang chứa mẫu đo độ thấm khí cho loại mẫu đường kính cm,10 cm Hình a)chứa Các thiết bị đo áp suất, nhiệt độ, lưu lượng quy đổi tín hiệu; b) Khoang chứa mẫu đo chứa mẫu thấmkhí khí chocác cácloại loạimẫu mẫuđường đườngkính kính5cm,10cm 5cm,10cm mẫu đođo độđộ thấm cho Cụ thể độ chithấm tiết Cembureau biểu thị sơ đồ khíkhoang cho cácchứa loại mẫu mẫu loại đường kính cm, 10 cm thể chi khoangchứa chứamẫu mẫuloại loạiCembureau Cembureauđược đượcbiểu biểuthịthịdưới dướisơsơđồđồ CụCụ tiếttiếtcủa khoang ởthể Hình 5: chi sau: sau: Hình Sơ đồ đồ hoạt hoạt động động khoang khoang chứa Hình 55.: Sơ chứa mẫu mẫu thí thí nghiệm nghiệm Mẫu thí nghiệm hình trụ, mặt bên ôm sát lớp polime đàn hồi chống thấm, có đường kính đường kính mẫu thử (Hình 4b).Bộ phận có tác 2.2 Quy trình bị khí mẫu dụng chuẩn ngănHình cản lọt thành bên khoang làm thí nghiệm.Lớp polime lồng Hình : Sơ đồhoạt hoạt động khoang chứa mẫuthíthí nghiệm :5trơ Sơ đồqua động chứa mẫu nghiệm vào buồng khí có dạng hình xăm xe đạp Khi áp suất buồng khí đạt CátMẫu Fontainebleau sauhình sàng ởbên mắt lưới dao từ 300 µm đến mm, trộn có với Mẫu nghiệm hình trụ, mặtbên ơmsát sát lớppơlime pơlime chống thấm, cóxi 0,7MPa, áp suất tác dụng lên mặt bên mẫu đủđộng lớn để coi thành bên mẫu kín thíthí nghiệm trụ, mặt ơm lớp chống thấm, măng pooc lăng với tỉ lệ cấp phối Bảng1 Các mẫu đúc khn hình trụ nhựa khí Mặt mặt mẫu tiếp xúc với hai đĩa kim loại có rãnh chia Những rãnh đường kính đường kínhcủa củamẫu mẫuthử thử(Hình (Hình4b) 4b).Nó Nócócótáctácdụng dụngngăn ngăncản cảnkhí khí đường kính đường kính chia cóđường tác dụng đảm bảo phân bố đều5 đặn ápđó suất trênngâm tiết diện mẫu cứng suốt có kính 10 cm, cao cm sau vào nước suốt trơ lọt qua thành bên làm thí nghiệm Lớp pơlime lồng vào buồng trơ lọt qua thành bên làm thí nghiệm Lớp pơlime lồng vào buồng trình kết phun (Hình 5) q trình đơng dạng hình xămxexe đạp.Khi Khiápápsuất suấttrong trongbuồng buồngkhí khíđạt đạt0,7MPa, 0,7MPa,ápápsuất suấttác tácdụng dụng khíkhí cócó dạng hình xăm đạp Cơng tác chuẩn bị mẫu dùng để đo độ thấm khí tóm lược gồm bước sau [6]: mặt bên mẫu lớn thànhbên bên mẫu kín khí Mặt mặt lênlên mặt bên mẫu sẽsẽ đủđủ lớn đểđể coicoi thành mẫu kín khí Mặt mặt mẫu tiếp xúc với kim loại có rãnh chia Những rãnh chia nàycóCác cótác tác dụng đảmvới tiếp xúc với hai đĩađĩa kim loại cóđường rãnh chia Những rãnh chia ◦mẫu Đo kích thước mẫu 3hai lần (chiều cao, kính) với độ xác 0,1này mm vịdụng trí đođảm tạo bảo phân suấttrên trênthiết thiếtdiện diệnmẫu mẫutrong trongquá quátrình trìnhphun phun(Hình (Hình5).5) ◦đặn bảo phân bốbố ápápsuất góc 120 đặn ◦ Dùng giấy ráp làm nhẵn bề mặt mẫu, cho mặt mặt mẫu song song với vng góc với trục đối xứng mẫu 139 66 Quân, T S / Tạp chí Khoa học Công nghệ Xây dựng ◦ Bọc thành bên mẫu lớp bạc mỏng (mẫu lớn) lớp cao su co nhiệt (mẫu nhỏ) Lớp có tác dụng tăng độ nhẵn bề mặt tiếp xúc mẫu khoang chứa mẫu, tránh tượng khí lọt qua thành bên ◦ Mỗi tổ gồm mẫu đo, ta đánh số mẫu đo theo thứ tự EC%-& chẳng hạn EC1-2, EC3-2, EC4-3 (EC: Echantillon du Ciment; %: tên tổ mẫu; &: thứ tự tổ mẫu) Bảng Tỷ lệ vật liệu theo khối lượng chế tạo mẫu hình trụ có đường kính 10 cm, chiều cao cm Tên tổ mẫu Xi măng/Cát Nước/Xi măng EC1 (4 mẫu) EC2 (4 mẫu) EC3 (4 mẫu) EC4 (4 mẫu) 1,0 1,2 1,4 1,6 0,4 0,4 0,5 0,5 2.3 Quy trình sấy ổn định độ ẩm mẫu nhiệt độ phòng Một nghiên cứu gần đây, Skoczylas [7] loại vữa xi măng có tỉ lệ Nước/Xi măng 0,5, vật liệu tạo thành sau đóng rắn gần đồng sau sấy 60◦ C Sự thay đổi cấu trúc rỗng vật liệu phần lớn gây bay nước tự do: kích thước lỗ rỗng trở nên lớn tính kết nối chúng tăng lên Hai tác nhân kể đến q trình gia nhiệt biến đổi hóa lý kèm theo xuất ứng suất nhiệt: Những bọt khí bị lập bên vật liệu có xu hướng nở ra, làm hình thành nên vi nứt thứ cấp Hệ dẫn đến tăng độ rỗng độ thấm tác dụng nhiệt độ Khi nhiệt độ sấy đạt ngưỡng 105◦ C, vết nứt có xu hướng phát triển số lượng kích thước với bốc nước tự Journal Science and Technology inCivil CivilHiệu Engineering NUCE2018 2018.13(5):1-16 13(5):1-16 phần nước liênofof kết trongand cấu trúc bê in tông ứng Klinkenberg rằng, sấy Journal Science Technology Engineering NUCE với thời gian đủ lớn, nhiệt độ sấy dao động khoảng 60◦ C-105◦ C độ thấm gần khơng đủlớn, lớn,nếu nếunhiệt nhiệtđộ độsấy sấydao daođộng độngtrong trongkhoảng khoảng60°C-105°C 60°C-105°Cthì thìđộ độthấm thấmgần gầnnhư nhưkhông không đủ thay đổi thayđổi đổi thay (a)ở nhiệt độ độ ẩm xác định ; b) Bình thủy (b)tinh ổn định độ ẩm Hình66: :a)a)Tủ Tủsấy sấymẫu mẫu Hình nhiệt độ độ ẩm xác định ; b) Bình thủy tinh ổn định độ ẩm mẫu sau sấy mẫu sau sấy Hình a) Tủ sấy mẫu nhiệt độ độ ẩm xác định; b) Bình thủy tinh ổn định độ ẩm mẫu sau sấy Nhưvậy vậycác cácbước bướcsấy sấyvà vàổn ổnđịnh địnhđộ độẩm ẩmmẫu mẫucó cóthể thểtóm tómlược lượcnhư nhưsau: sau: Như Như cácoobước vàmẫu ổn định độ ẩm mẫu sau: Sấysấy khô mẫu thửnhằm nhằm loại bỏđược nướcthực cótrong mẫu.Đối Đốivới vớimẫu mẫunhỏ, nhỏ, Sấy khơ thử loại bỏ nước tựtựdo có mẫu qtrình trìnhgia gianhiệt nhiệtsẽsẽkhiến khiếnlớp lớpcao caosu suco conhiệt nhiệtôm ômsát sátvào vàothân thânmẫu mẫu(Hình (Hình 140 6a) 6a) Cânđể đểxác xácđịnh địnhkhối khốilượng lượngmẫu mẫuđo đo oo Cân o Ổn định độ ẩm mẫu, đưa nhiệt độ mẫuđo đovề vềnhiệt nhiệtđộ độkhí khíquyển quyểnnhờ nhờbình bìnhthủy thủy o Ổn định độ ẩm mẫu, đưa nhiệt độ mẫu tinh dessicateur chứa hỗn hợp muối hút ẩm, nhằm đảm bảo mẫu có độ ẩm xác tinh dessicateur chứa hỗn hợp muối hút ẩm, nhằm đảm bảo mẫu có độ ẩm xác Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng ◦ Sấy khô mẫu thử nhằm loại bỏ nước tự có mẫu Đối với mẫu nhỏ, q trình gia nhiệt khiến lớp cao su co nhiệt ôm sát vào thân mẫu (Hình 6(a)) ◦ Xác định khối lượng mẫu đo ◦ Ổn định độ ẩm mẫu, đưa nhiệt độ mẫu đo nhiệt độ khí nhờ bình thủy tinh chứa hỗn hợp muối hút ẩm, nhằm đảm bảo mẫu có độ ẩm xác định (Hình 6(b)) Để thiết lập mối quan hệ độ thấm nội độ bão hòa, hàng ngày mẫu cân ổn định nhiệt tiến hành thí nghiệm đo thấm khí, sau lại đưa trở lại buồng sấy Quy trình lặp lặp lại đến khối lượng mẫu không đổi 2.4 Quy trình đo độ thấm khí mẫu Khí azote phun từ lên, qua mẫu đo mặt khoang chứa mẫu có áp suất áp suất khí Áp suất phun dao động từ đến 0,5 MPa, giới hạn cho phép điều áp Giá trị tương ứng với giới hạn lưu lượng cực đại đo thiết bị (500 ml/phút 1000 ml/phút).Trước tiến hành thí nghiệm mẫu, ta cần phải kiểm tra độ kín thiết bị đo thơng qua mẫu tiêu chuẩn kim loại (Hình 4(b)) với áp suất khí phun 0,4 MPa Nếu lưu lượng khí nhỏ 5% lưu lượng cực đại phát thí nghiệm tiến hành Mục đích thí nghiệm để xác định độ thấm nội Kv thơng qua độ thấm khí Ka Để làm điều này, ta cần tác dụng giá trị áp suất khác đầu vào mẫu Với giá trị áp suất phun, ta có giá trị độ thấm tường Ka Quy trình đo độ thấm tóm lược sau: ◦ Tác dụng tối thiểu giá trị áp suất khác nhằm có mối quan hệ tuyến tính giá trị Ka theo nghịch đảo áp suất phun trung bình sơ đồ minh họa Hình Áp suất phun giảm bậc từ cao xuống thấp ngược lại Cách thức chọn giúp đánh giá ảnh hưởng áp suất phun đến kết thu ◦ Với giá trị áp suất tác dụng, giá trị lưu lượng phải lưu khoảng thời gian tối thiểu 15 phút Nếu giá trị lưu lượng thu khác 2%, ta cần phải chờ cho dòng chảy ổn định thu thập lại số liệu từ đầu 2.5 Kết thí nghiệm thu Các mẫu thí nghiệm sấy nhiệt độ 80◦ C đến khối lượng không đổi nhằm đảm bảo q trình gia nhiệt khơng làm thay đổi cấu trúc vật liệu Những thí nghiệm trước rằng, với mẫu hình trụ tương tự có khối lượng từ 700 g đến 900 g, cần thời gian khoảng ngày để đạt khối lượng ổn định không đổi Các mẫu lựa chọn thí nghiệm có đường kính 10 cm, chiều cao cm, khối lượng ban đầu dao động từ 800 g đến 880 g Lượng nước bay chiếm khoảng 8% (từ 60 g đến 75 g) mẫu sau trình xử lý nhiệt Khối lượng mẫu sau sấy dao động từ 710 g đến 790 g (Hình 7) Trong cơng thức (3), độ thấm tường Ka đặc trưng cho khả thẩm thấu khí gas qua hệ thống lỗ rỗng vật liệu Sự bay nước tự làm tăng đường kính lỗ rỗng tính kết nối, kéo theo giảm dần hệ số β tăng lên độ thấm nội Kv Sau bước sấy, mẫu lại đưa đo độ thấm khí Điều cho phép ta thiết lập mối quan hệ hệ số β, độ thấm nội Kv độ bão hòa S w , xác định công thức: Vw Vv 141 Sw = (4) Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng Vv = M sat − Msec ρw (5) Vw thể tích nước trạng thái (cm3 ); Vv thể tích rỗng vật liệu (cm3 ); M sat khối lượng mẫu bão hòa nước (g) Msec khối lượng mẫu sấy khô (g) ρw khối lượng riêng nước (g/cm3 ) Lượng nước/Thời gian Áp dụng công thức (4) cho mẫu đo, ta nhận ΔM (g/ ngày) 40 Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 thấy độ bão hòa ban đầu mẫu dao động từ Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 EC1 trung bình 35 70% - 85% EC2 trung bình VW 30 Nhìn chung, ta đo độ thấm khí EC3 trung bình V SW = W25 SWV= EC4 trung bình độ bão hòa 80% Trên ngưỡng 80%, V (4) VV 20 (4) lượng nhỏ khí trơ xun qua mẫu 15 M M Căn vào biến thiên hệ số Klinkenberg satM -secM VV = sec sat VV = r 10 theo độ bão hòa (Hình 8(a)), mẫu khơ hồn w r (5)(5) w tồn, sai lệch giá trị tương đối thấp Điều với với lý giải giảm tác dụng dịch Thời gian (ngày) chuyển trượt so với dịch chuyển nhớt Trong 3 Thểcó tích nướcđổi trạng thái (cm ) rỗng ; 3) ; VV :VThể tíchtích rỗng củacủa vật vật liệuliệu (cm(cm ) ;3) ; VWđó:Vdo thay mạnh mẽ cấu trúc : Thể tích nước trạng thái (cm : Thể rỗng V W Hình Biến thiên độ nước theo thời vật liệu, quan sát thấy Hình 8(b), độ thấm lượng mẫumẫu bão hòa hòa nướcnước (g) ;(g) ; M secM: Khối lượng mẫumẫu khi4khi sấymẫu khôkhô (g)(g) ; ; gian tổ M sat M: Khối khisựbão : Khối lượng sấy nội sat tại: ởKhối trạnglượng thái khơ có biến thiên lớn sec mẫu đo.riêng nước (g/cm3) ; :r Khối lượng rWgiữa W : Khối lượng riêng nước (g/cm ) ; Hiệu ứng Hiệu Klinkenberg ứng Klinkenberg β (Pa)β (Pa) Kv (m2) Kv (m^2) 1,4E-14 1,4E-14 1200000 1200000 1000000 1000000 800000800000 EC1 trung EC1bình trung bình EC3 trung EC3bình trung bình EC4 trung EC4bình trung bình 1E-14 1E-14 8E-15 8E-15 6E-15 6E-15 400000400000 4E-15 4E-15 200000200000 2E-15 2E-15 0 EC1 trung bình EC1 trung bình EC2 trung bình EC2 trung bình EC3 trung bình EC3 trung bình EC4 trung bình EC4 trung bình 1,2E-14 1,2E-14 EC2 trung EC2bình trung bình 600000600000 Độ thấm nội tại/Độ bão hòa Độ thấm nội tại/Độ bão hòa 00,1 0,10,2 0,20,3 0,30,4 0,40,5 Độ bão (Sw) Độhòa bão hòa (Sw) 0,50,6 0,60,7 0,7 0 0,1 0,1 0,2 a) (a)a) 0,2 0,3 0,3 0,4 0,4 0,5 0,5 0,6 0,6 0,7 0,7 Độ bão (Swhòa ) Độhòa bão (Sw) (b) b) b) Hình Biến thiên bãohòa hòa a)a)Hệ Klinkenberg β;b b);bĐộ thấm nội Knội vtại HìnhHình : Biến thiên theo độtheo bão hòa :của: a) :Hệ số b); Độ thấm nội : Biến thiên theo độđộbão HệsốKlinkenberg số Klinkenberg b) Độ thấm KV KV Áp dụng thức (4) cho đo, tanội nhận độ bão Áplýdụng công thức (4)trong cho các mẫu đo, ta nhận thấy độ bão hòahòa banban đầuđầu củacủa cáccác Các thuyết mơcơng tính tốnmẫu độ thấm tạithấy dao động 85 % mẫu mẫu dao động từ 70từ%70- % 85 -% Chất lỏng lỗ rỗng thường tồn dạng tổ hợp nước lượng nhỏ chất lỏng chung, ta có thể đo độ thấm 80% Trên NhìnNhìn chung, ta đo độ thấm khí khí khi độ độ bãobão hòahòa 80% Trên khơng hòa tan Độ thấm nội trạng thái bão hòa xác định Kv (S w ) phụ thuộc vào độ thấm ngưỡng lượng nhỏ trơthể có thểxuyên xuyên mẫu vào biến ngưỡng 80%,80%, lượng nhỏ khí khí trơ có quaqua mẫu CănCăn vào biến trạng thái nước Kv ảnh hưởng có mặt pha lỏng pha khí Độ thấm trạng củasốhệKlinkenberg số Klinkenberg độ bão (Hình khơ hồn tồn, thiênthiên hệ theotheo độ bão hòa hòa (Hình 8a),8a), khi mẫumẫu khơ hồn tồn, sai sai thái lệch chưagiữa bão hòa vật liệu sở xi măng xác định công thức sau: giá trị tương đối thấp Điều lý giải giảm tác dụng dịch lệch giá trị tương đối thấp Điều lý giải giảm tác dụng dịch chuyển so dịch với dịch chuyển nhớt Trong có dosự cóthay thay mạnh chuyển trượttrượt so với chuyển nhớt Trong đổi đổi mạnh mẽmẽ cấucấu (6) Kv (S e ) = Kđó v Kr (S e ) trúc rỗng vật liệu, quan sát thấy Hình 8b, độ thấm nội trạng thái khơ trúc rỗng vật liệu, quan sát thấy Hình 8b, độ thấm nội trạng thái khơ có có sự lớn mẫu biến biến thiênthiên lớn mẫu đo đo 142 Các lý thuyết mô độ thấm Các lý thuyết mơ tínhtính tốntốn độ thấm nội nội tại lỗ rỗng thường tồndưới dạng tổ hợp nước lượng ChấtChất lỏnglỏng trongtrong lỗ rỗng thường tồn dạng tổ hợp củacủa nước lượng nhỏnhỏ chất lỏng khơng hòa tan Độ thấm nội trạng thái bão hòa xác định chất lỏng khơng hòa tan Độ thấm nội trạng thái bão hòa xác định K K W ) phụ (SV ( S) phụ Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng S e độ bão hòa hiệu dụng; Kv (S e ) độ thấm nội theo độ bão hòa hiệu dụng (m2 ); Kv độ thấm nội trạng thái khô tuyệt đối (m2 ); Kr (S e ) độ thấm tương khí trơ theo độ bão hòa hiệu dụng Độ bão hòa hiệu dụng xác định mối quan hệ sau: Se = Sw − Sr Ss − Sr (7) S s độ bão hòa trạng thái bão hòa (S s ≈ 1); S w độ bão hòa trạng thái bất kỳ; S r độ bão hòa trạng thái khơ (S r ≈ 0) Để đơn giản hóa, q trình tính tốn, ta lấy xấp xỉ sau: S e ≈ S W Vấn đề đặt xác định độ thấm tương đối tương ứng với giá trị độ bão hòa Hai mơ hình ứng dụng nhiều để xác định độ ẩm tương đối mơ hình Burdine [8] Mualem [9] Hai mơ hình khác quan điểm cấu trúc lỗ rỗng Trong mơ hình mình, Burdine mơ hình hóa hệ thống lỗ rỗng nhóm ống mao dẫn song song với nhiều đường kính khác Mơ hình Mualem phức tạp xét đến tương tác kênh lỗ rỗng Đối với vật liệu bê tơng nói chung, mơ hình Mualem thường ưu tiên sử dụng tính tốn độ thấm Tuy nhiên mơ hình có ưu nhược điểm, việc chọn lựa mơ hình tính cần ý cho phù hợp 3.1 Mơ hình Brooks-Corey Corey [10] kết hợp mối quan hệ áp suất độ bão hòa Brooks-Corey thiết lập [11] với mơ hình Burdine để đưa cơng thức sau: 2+λ Kr (S w ) = (1 − S w )2 − S wλ (8) λ số phân bố kích thước lỗ rỗng Brooks-Corey Nghiên cứu số phân bố kích thước lỗ rỗng λ mơ hình áp suất mao dẫn Brooks-Corey đặc trưng cho tính khơng đồng mơi trường rỗng: số phân bố kích thước lỗ rỗng lớn, mơi trường rỗng đồng Trong công thức (8), giá trị λ tiến tới 0, ta có Kr (S w ) = (1 − S w )2 Trái lại, giá trị λ tiến tới vô cùng, ta có Kr (S w ) = (1 − S w )3 Ở trường hợp này, đường cong có xu hướng giảm nhanh so với trường hợp trước Từ ta rút nhận định: Khi giá trị λ lớn, với môi trường rỗng có tính đồng cao, độ thấm nội có xu hướng giảm mạnh tăng lượng nước tự mẫu 3.2 Mơ hình Van Genuchten-Mualem Mualem [9] đưa mơ hình mơ dựa giả thiết lỗ rỗng kết nối với theo chiều dài tỉ lệ tương ứng với bán kính lỗ rỗng Mơ hình giả thiết định luật Poiseuille áp dụng cho lỗ rỗng thành phần Với trạng thái bão hòa xác định, ngoằn ngoèo hệ thống lỗ rỗng, thông số biểu diễn tương tác lỗ rỗng có đường kính khác nhau, diễn giải hàm số mũ độ bão hòa: 1 Kr (S w ) = (1 − S w ) p − S wm 143 2m (9) Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng m thông số thực nghiệm; p thông số hàm mũ Thông số m lấy từ đường cong biểu diễn mối quan hệ áp suất mao dẫn độ bão hòa nước Van Genuchten [12], sở lý thuyết xây dựng Mualem, xác định hệ số này, đồng thời miền ứng dụng vật liệu có nguồn gốc từ xi măng Sau so sánh kết thí nghiệm 45 loại xi măng đất khác nhau, giá trị hàm mũ p dao động khoảng [0,18; 2] Kết mơ thu Trong thí nghiệm đo độ thấm nội tại, giá trị đo khả tín trạng thái đánh giá “khơ tuyệt đối” khơng có tham gia nước Ở trạng thái ẩm phun khí qua mơi trường, nước hấp thụ lượng khí định tạo nên sai số phép đo Sai số lớn mơi trường có độ rỗng cao, tính kết nối lỗ rỗng lớn Do q trình mơ phỏng, mơ hình phải lấy giá trị trạng thái khô làm sở để thiết lập đường cong quan hệ Với giá trị xác định độ bão hòa, thực nghiệm ta thu giá trị độ thấm nội Kv tương ứng So sánh sai số tương đối giá trị giá trị mơ phỏng, ta kết luận mức độ phù hợp phương pháp với loại vật liệu mà ta nghiên cứu Bảng Các tham số hai phương pháp mô tối ưu hóa 16 mẫu đo Tên mẫu EC1-1 EC1-2 EC1-3 EC1-4 EC2-1 EC2-2 EC2-3 EC2-4 EC3-1 EC3-2 EC3-3 EC3-4 EC4-1 EC4-2 EC4-3 EC4-4 Max Min Trung bình Độ lệch chuẩn Trạng thái sấy khơ Brooks-Corey Kv (m2 ) λ p m 3,63E-15 6,08E-15 4,52E-14 2,93E-14 1,24E-14 3,89E-15 2,33E-14 4,36E-14 3,84E-15 1,26E-14 2,87E-14 8,35E-15 1,16E-15 1,22E-14 5,71E-15 9,79E-15 4,52E-14 1,16E-15 1,56E-14 1,37E-14 0,01 0,01 0,01 0,01 0,01 2,92 0,01 3,6 10000 0,01 0,01 10000 10000 10000 1,8 0,01 10000 0,01 2501 4330 0,6 0,4 0,6 1,4 0,8 0,3 0,6 0,5 0,3 0,6 1,4 0,7 0,3 0,2 0,4 1,0 1,40 0,23 0,64 0,35 0,10 0,11 0,13 0,13 0,13 0,15 0,12 0,13 0,15 0,10 0,10 0,97 1,30 0,13 0,14 0,14 1,30 0,10 0,25 0,34 144 Van Genuchten-Mualem Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 Kv(m^2) Độ thấm nội / Độ bão hòa 5,00E-14 4,50E-14 Kết thí nghiệm 4,00E-14 Brooks-Corey Quân, T S / Tạp chí Khoa học Cơng nghệ Xây dựng 4.1 Mơ hình Brooks-Corey 3,50E-14 3,00E-14 Theo Kewen Li [13], giá trị số λ dao động khoảng [0,3; 1,9] Tuy nhiên 2,50E-14 mẫu xem xét, ta nhận thấy số λ nằm khoảng 10−2 , 104 (Bảng 2), giá trị độ thấm 2,00E-14 nội gần không thay đổi rõ rệt Mô hình gần mơ tốt mẫu đo có giá 1,50E-14 trị biến thiên độ thấm nội không lớn 1,00E-14 5,00E-15 4.2 Mô hình Van Genuchten-Mualem 0,00E+00 0,000 0,200 0,400 0,600 0,800 1,000 Trong biểu thức (9), cấu trúc vi mô vật liệu có ảnh hưởng lớnĐộđến giáSwtrị thơng số p bão hòa Điều giải thích khác biệt cấu trúc vi mô bê tơng đất, đặc Hình : So sánh kết thí nghiệm 16 mẫu đo mô theo Brooks-Corey trưng độ rỗng, phân bố hệ thống lỗ rỗng, ngoằn ngoèo, tính kết nối độ nhám 4.2 Mơ hình Van Genuchten-Mualem thành lỗ rỗng Giá trị p đưa J Verdier [14], daobiểu động đốiảnhvới bê tông Trong thứctrong (9), cấukhoảng trúc vi mô[0,18; vật2,0] liệu có hưởng lớnlàđến giá trị thơng Điềutrị nàynày giải thíchsát bởinhất khác mốc so sánh hữu ích Ta số vàop giá để mơ cóbiệt thểtrong cấu giátrúc trị vi mô bê tôngnghiệm đất, đặc trưng độ rỗng, phân bố hệ thống lỗ rỗng, ngoằn ngoèo, độ thấm nội đo thông qua thực kết nối1,4], độthông nhám lỗkhoảng rỗng [0,1; 1,3] (Bảng 2), Với giá trị p khoảngtính [0,23; số thành m Giá trị p đưa Verdierkhớp [14], dao khoảng giá trị mô mơ hình Van Genuchten-Mualem gần nhưJ trùng vớiđộng giá trị đo [0,18 ; 2,0] bê tơng mốc so sánh hữu ích.Ta vào giá trị để mô thực nghiệm Journal of Science and Technology in Civil Engineering NUCE 2018 13(5):1-16 sát giá trị độ thấm nội đo thông qua thực nghiệm Kv(m2) Độ thấm nội / Độ bão hòa Độ thấm nội / Độ bão hòa Kv(m2) 5,00E-14 5,00E-14 4,50E-14 4,00E-14 Kết thí nghiệm 4,50E-14 Kết thí nghiệm Brooks-Corey 4,00E-14 Van Genuchten-Mualem 3,50E-14 3,50E-14 3,00E-14 3,00E-14 2,50E-14 2,50E-14 2,00E-14 2,00E-14 1,50E-14 1,50E-14 1,00E-14 1,00E-14 5,00E-15 5,00E-15 0,00E+00 0,00E+00 0,000 0,200 0,400 0,600 0,800 0,000 1,000 0,200 0,400 0,600 0,800 1,000 Độ bão hòa Sw Độ bão hòa Sw Hình 10 : So sánh kết thí nghiệm 16 mẫu đo mơ theo Van Hình : So Hình sánh kết thí nghiệm trênthí 16 nghiệm mẫu đo vàtrên mơ 16 So sánh kết mẫutheo đoBrooks-Corey Hình 10 So sánh kết thí nghiệm 16 mẫu đo Genuchten-Mualem mô theo Brooks-Corey mô theo Van Genuchten-Mualem 4.2 Mơ hình Van Genuchten-Mualem Trong biểu thức (9), cấu trúc vi mơ vật liệu có ảnh hưởng lớn đến giá trị thông số p Điều giải thích khác biệt cấu trúc vi mô bê Kết luậnbởi chung tông đất, đặc trưng độ rỗng, phân bố hệ thống lỗ rỗng, ngoằn ngoèo, tính kết nối độ nhám thành lỗ rỗng 14 Thông qua việc đo độ thấm khí mẫu nhiều trạng thái độ ẩm khác nhau, ta xác Giá trị p đưa J Verdier [14], dao động khoảng [0,18 ; 2,0] định độ thấm nội tương ứng, đại lượng phụ thuộc vào tính chất vật liệu mà không phụ bê tông mốc so sánh hữu ích.Ta vào giá trị để mơ vào tính chất lưu thấm nó.thơng Đây qua sởnghiệm liệu quan trọng cho phép ta tiến hành bước sát có đặc thể giácủa trị độ thấm nội qua đo thực mô Độ thấm nội / Độ bão hòa Kv(m^2) Đối với mơ hình Brooks-Corey, giá trị λ sử dụng vượt khoảng sở lí thuyết 5,00E-14 nhiều trường hợp, vài giá trị sai số mô lên đến 50% Trái lại trường hợp 4,50E-14 Kết thí nghiệm Van Genuchten-Mualem, sai số tương đối 10%, giới hạn chấp nhận công 4,00E-14 Van Genuchten-Mualem tác thực nghiệm Về mặt tốn học, mơ hình Van Genuchten-Mualem linh hoạt chứa tham 3,50E-14 3,00E-14 2,50E-14 2,00E-14 1,50E-14 1,00E-14 145 Quân, T S / Tạp chí Khoa học Công nghệ Xây dựng số, giá trị p dao động miền giá trị sở lý thuyết cơng trình nghiên cứu trước Đối chiếu giá trị tính tốn theo mơ hình Van Genuchten-Mualem Brooks-Corey với kết thí nghiệm cho phép thiết lập nên miền xác định tham số có liên quan Dựa sở thơng số này, dự đốn (thông qua nội suy ngoại suy) giá trị trung gian độ thấm nội trạng thái khác độ bão hòa mà khó xác định phương pháp thực nghiệm Đối với loại vật liệu nghiên cứu, mơ hình Van Genuchten-Mualem tỏ hiệu thích hợp so với mơ hình Brooks-Corey Điều lý giải giả thiết cấu trúc vi mô vật liệu, đặc biệt vai trò hấp thụ truyền chất hệ thống lỗ rỗng hai mơ hình khác Tài liệu tham khảo [1] Darcy, H (1856) Les fontaines publiques de la ville de dijon Victor Dalmont, Paris [2] Picandet, V (2001) Influence d’un endomagement mécanique sur la perméabilité et sur la diffusivité hydrique des bétons These de Doctorat [3] Klinkenberg, L G (1941) The permeability of porous media to liquids and gases Drill Prod Pract API 200–213 [4] Kollek, J J (1989) The determination of the permeability of concrete to oxygen by the Cembureau method–a recommendation Materials and Structures, 22(3):225–230 [5] Cassanas, M N (2005) Etude expérimentale d’injection de coulis de ciment dans un milieu granulaire Perméabilité au gaz Rapport de stage, Institut de Recherche en Génie Civil et Mécanique (GéM), Ecole centrale de Nantes [6] Bouasker, M (2007-2008) Etude de propagation de coulis dans des sols injectés, mesure de perméabilité Mémoire de DEA Génie Civil (Eau, environnement et génie côtier), Institut de Recherche en Génie Civil et Mécanique (GéM), Ecole centrale de Nantes [7] Lion, M., Skoczylas, F., Lafhaj, Z., Sersar, M (2005) Experimental study on a mortar Temperature effects on porosity and permeability Residual properties or direct measurements under temperature Cement and Concrete Research, 35(10):1937–1942 [8] Burdine, N T (1953) Relative permeability calculations from pore size distribution data Journal of Petroleum Technology, 5(3):71–78 [9] Mualem, Y (1976) A new model for predicting the hydraulic conductivity of unsaturated porous media Water Resources Research, 12(3):513–522 [10] Corey, A T (1994) Mechanics of immiscible fluids in porous media Water Resources Publication [11] Brooks, R H., Corey, A T (1964) Hydraulic properties of porous media Hydrol Pap 3, Colo State Univ., Fort Collins [12] Van Genuchten, M T (1980) A closed-form equation for predicting the hydraulic conductivity of unsaturated soils Soil Science Society of America Journal, 44(5):892–898 [13] Kewen, L (2004) Theoretical development of the Brooks-Corey capillary pressure model from fractal modeling of porous media Fourteenth Symposium on Improved Oil Recovery held in Tulsa, Oklahoma, U.S.A [14] Verdier, J., Carcasses, M., Ollivier, J P (2002) Modelling of a gas flow measurement: Application to nuclear containment vessels Cement and Concrete Research, 32(8):1331–1340 146 ... ta đo độ thấm khí khí khi độ độ bãobão hòahòa 80% Trên khơng hòa tan Độ thấm nội trạng thái bão hòa xác định Kv (S w ) phụ thuộc vào độ thấm ngưỡng lượng nhỏ trơthể có thểxuyên xuyên mẫu vào biến... 0,7 Độ bão (Swhòa ) Đ hòa bão (Sw) (b) b) b) Hình Biến thiên bãohòa hòa a)a)Hệ Klinkenberg β;b b);bĐộ thấm nội Knội vtại HìnhHình : Biến thiên theo độtheo bão hòa :của: a) :Hệ số b); Độ thấm. .. EC2 trung EC2bình trung bình 600000600000 Độ thấm nội tại /Độ bão hòa Độ thấm nội tại /Độ bão hòa 00,1 0,10,2 0,20,3 0,30,4 0,40,5 Độ bão (Sw) Đ hòa bão hòa (Sw) 0,50,6 0,60,7 0,7 0 0,1 0,1 0,2