Báo cáo nghiên cứu khoa học: " DÒNG THẤM TRONG ĐẤT KHÔNG BÃO HOÀ"
Chia sẻ: Nguyễn Phương Hà Linh Nguyễn Phương Hà Linh | Ngày: | Loại File: PDF | Số trang:8
lượt xem 8
download
Bài báo đề cập đến các phương trình dòng thấm nước trong toàn bộ nền đất, bao gồm trong đất bão hoà và cả đất không bão hoà, như vậy mới giải quyết toàn bộ bài toán thấm trong nền đất. Trong tính toán, cần thiết lập một phương trình tổng quát cho cả hai trường hợp, như vậy, việc tính toán dòng thấm từ đất không bão hoà đến bão hoà mới diễn ra “êm thuận”. Đối với trường hợp thấm trong đất bão hoà có thể lấy “độ hút dính” bằng 0, và như vậy, hệ số...
Bình luận(0) Đăng nhập để gửi bình luận!
Nội dung Text: Báo cáo nghiên cứu khoa học: " DÒNG THẤM TRONG ĐẤT KHÔNG BÃO HOÀ"
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 DÒNG THẤM TRONG ĐẤT KHÔNG BÃO HOÀ THE PERCOLATION FLOW IN UNSATURATED SOILS Bạch Quốc Tiến Trường Đại học Bách Khoa, Đại học Đà nẵng TÓM TẮT Bài báo đề cập đến các phương trình dòng thấm nước trong toàn bộ nền đất, bao gồm trong đất bão hoà và cả đất không bão hoà, như vậy mới giải quyết toàn bộ bài toán thấm trong nền đất. Trong tính toán, cần thiết lập một phương trình tổng quát cho cả hai trường hợp, như vậy, việc tính toán dòng thấm từ đất không bão hoà đến bão hoà mới diễn ra “êm thuận”. Đối với trường hợp thấm trong đất bão hoà có thể lấy “độ hút dính” bằng 0, và như vậy, hệ số thấm nước của đất không bão hoà kw bằng hệ số thấm nước bão hoà ks. Các phương trình vi phân của dòng thấm có thể giải bằng các phương pháp: sai phân hữu hạn, phần tử hữu hạn . ABSTRACT The paper deals with the equation of the percolation flow in all-soil, including saturated soils and unsaturated soils, so as to completely solve the problem about percolation flow in all types of soils . For computation, it is necessary to make full equations for the above two cases so as to make it easier to compute the percolation flow from unsaturated soils to saturated soils. In case of the percolation flow in the saturated soils, “matric suction” may be equaled zero and so that unsaturated soils permeability coefficient (kw) equals saturated soils permeability coefficient (ks). The equations can be solved by numerical methods such as finite difference and finite element methods. 1. Mở đầu Với các điều kiện về thời tiết và khí hậu trong tự nhiên, việc thấm nước trong các lớp đất không bão hòa trên bề mặt nền đất là một thực tế tất yếu, mặc dầu vậy, trong việc tính toán dòng thấm nước của đất thì thông thường áp dụng cho đất bão hòa mà đã được đề cập đến rất nhiều trong các tài liệu địa chất thuỷ văn, địa kỹ thuật... Công tác tính toán thấm trong đất không bão hòa vẫn còn hạn chế do ít tài liệu đề cập đến vấn đề này. Công tác tính thấm nước trong đất không bão hòa có thể được áp dụng nhiều trong các trường hợp thực tế xây dựng công trình: Bài tóan tính thấm qua đập đất, bài tóan tính lượng nước ngấm trên bề mặt địa hình không ngập nước v.v..., như vậy, để bài toán thấm nước trong toàn bộ nền đất được giải quyết triệt để, phải có phương trình tổng quát đặc trung cho dòng thấm nước áp dụng cho cả trong đất không bão hòa và bão hòa. 2. Định luật Darcy cho đất không bão hòa Dòng thấm nước mô tả bằng định luật Darcy được phát biểu như sau: Vận tốc dòng thấm nước qua khối đất tỷ lệ (quan hệ bậc nhất- tuyến tính) với gradien cột nước thuỷ lực. 23
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 V w = -k w h w y Trong đó: VW: Vận tốc thấm nước. kW: Hệ số thấm đối với pha nước. dhw/dy: Gradien cột nước thuỷ lực theo hướng y (ký hiệu là iw). Hệ số tỷ lệ giữa vận tốc thấm của nước và gradien cột nước thuỷ lực được gọi là hệ số thấm (kw). Đối với đất bão hoà, hệ số thấm có thể xem là hằng số. Dấu âm trong phương trình biểu thị là các dòng thấm nước vận động theo hướng giảm cột nước thuỷ lực. Định luật Darcy cũng có thể áp dụng cho dòng thấm nước qua đất không bão hòa , tuy nhiên, hệ số thấm trong đất không bão hòa thường không được xem là hằng số mà là một biến, biến này chủ yếu là một hàm của độ ẩm hay độ hút dính của đất không bão hòa. 0. 6 Bột đá phiến 0. 5 Độ ẩm thể tích Cát (hạt 0.13mm đến 0.50mm) 0. 4 w ‱ Cát (hạt 0.5mm đến 1.0mm) 0. 3 0. 2 Gradien cột nước thuỷ lực bằng 1. Gradien cột nước thuỷ lực bằng 0.5. 0. 1 0 0. 2 0. 4 0. 6 0. 8 1. 0 1. 2 1. 4 1. 6 1. 8 Hệ số thấm k w ‱. 3‱ / ) (10- m s Hình 1. Xác minh thực nghiệm về định luật thấm Darcy cho dòng thấm nước qua đất không bão hòa (Theo Childs và Collis-Goerge) 3. Dòng thấm ở trạng thái ổn định Dòng thấm không bão hòa cũng được chi phối cùng định luật như đối với đất bão hòa (Định luật Darcy), chỉ khác là đối với đất không bão hòa thì hệ số thấm nước kw là hàm của lực hút dính, độ ẩm thể tích còn đối với đất bão hòa là hằng số, mặt khác, áp lực nước lổ rống trong đất bão hòa thường có giá trị dương còn đối với đất không bão hòa thường có giá trị âm. Mặc dầu có sự khác nhau này, cách thiết lập phương trình vi phân riêng phần cho cả hai trường hợp đều tương tự nhau và cũng có sự chuyển tiếp êm thuận từ trường hợp không bão hòa sang trương hợp bão hòa. 24
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 3.1. Biến đổi hệ số thấm trong đất không bão hòa Trên cơ sở hệ số thấm biến đổi theo các hướng trong không gian, có thể phân ra thành ba trường hợp ứng với các điều kiện khác nhau: + Thấm trạng thái ổn định, không đồng nhất, đẳng hướng (Hình 2a) + Thấm trạng thái ổn định, không đồng nhất, dị hướng (Hình 2b) + Thấm trạng thái ổn định, không đồng nhất, hệ số thấm thay đổi liên tục (Hình 2c) ky kx kx kx A kx B 1 A B kx = k y ky ky ky ky ky ky ky ky ky kx kx kx B B B kx kx kx A A A (a). Không đồng nhất,đẳng hướng (b). Không đồng nhất,dị hướng (c): Hệ số thấm thay đổi l iên tục Hình 2. Biến đổi hệ số thấm nước trong đất không bão hòa. 3.2. Thiết lập phương trình cho dòng thấm một hướng Xét một phân tố đất vô cùng nhỏ trong đất không bão hòa với kích thước dx,dy,dz (Hình 3). Tốc độ thấm vwy giả thiết là dương khi nước thấm lên theo hướng (y) vwy + vwy dy y Phân tố đất có bề dày dz y dy vw y x dx z Hình 3. Dòng thấm nước một hướng qua phân tố đất không bão hòa æ ¶v wy ö ¶v wy dy÷ dxdz - v wy dxdz = 0 hay ç ÷ dxdydz = 0 ç v wy ÷ ÷ ç ¶y ¶y è ø -k w y ( u a -u w ) . h w y Áp dụng định luật Darcy: (1) d yd x d z = 0 y Viết lại theo hướng y như sau (với kwy(ua- uw) viết lại là kwy): k wy h2w + k w y . h w = 0 2 (2) y y y Khi đất bão hòa không đồng nhất thì phương trình thấm có thể viết như sau 25
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 k s h2w + k s . h w = 0 2 (3) y y y Khi đất bão hòa đồng nhất, lúc đó hệ số thấm bão hòa ks là hằng số đối với nền đất. Như vây, phương trình vi phân dòng thấm có thể viết lại là: 2hw = 0 (4) y 2 Phương trình (4) cũng có thể thiết lập cho hướng x và z. Tóm tắt các phương trình vi phân dòng thấm trạng ổn định một hướng cho đất không bão hòa và bão hòa theo bảng sau: Bảng 1. Phương trình thấm ổn định một hướng của đất không bão hòa. Đất không đồng nhất k w y h2w + k w y . h w = 0 2 (2) y y y Bảng 2. Phương trình trạng thái thấm ổn định một hướng của đất bão hòa. Không đồng nhất Đồng nhất k s h2w + k s . h w = 0 ( 3) 2 ( 4) 2h w = 0 y y y y2 3.3. Thiết lập phương trình cho dòng thấm hai hướng: Phân tố đất có bề dày dz y dy vw x vw y x dx z Hình 4. Dòng thấm hai hướng qua phân tố đất không bão hoà Xét một phân tố đất có kích thước dx, dy, dz và giả sử dòng thấm theo hai hướng x và y (Hình 4). Hệ số thấm theo phương x là kwx(ua- uw) và theo phương y là kwy(ua- uw) (Viết tắt là kwx và kwy). Tỷ số kwx/kwy là không đổi tại mọi điểm trong khối đất. Vận tốc thấm vwx và vwy là hướng dương của x và y. Tính liên tục của dòng thấm ổn định hai hướng cho phép biểu thị như sau: 26
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 v w x+ v w x d x -v w x d y d z + v w y+ v w y d y - v w y d x d z = 0 x y Như vậy, lưu lượng thực thấm theo các hướng x,y sẽ là v wx v wy d x d y d z =0 x y Áp dung định luật Darcy: kwx hw kwy hw = 0 (5) + x x y y Phương trình (5) viết lại như sau: k w x h 2w + k w y h2w k wx . h w + k wy . h w (6) 2 2 =0 + x y x x y y Trong trường hợp đẳng hướng (kwx= kwy = kw), khi đó phương trình (6) có thể viết: k w h 2w h 2w + k w . h w + k w . h w (7) 2 2 =0 x y x x y y Khi đất bão hoà, hệ số thấm nước kwx,kwy có thể lấy bằng hệ số thấm nước bão hòa ksx và ksy k s x h 2w k s y h 2w k sx . h w + k sy . h w (8) 2 2 =0 + + x y x x y y Đối với đất bão hòa không đồng nhất, đẳng hướng thì phương trình (8) có thể viết như sau: k s h2w h2w + k s . h w + k s . h w (9) 2 2 =0 x y x x y y Đối với đất bão hòa đồng nhất, dị hướng thì phương trình thấm có thể viết như sau: k s x h 2w k s y h 2w (10) 2 2 =0 + x y Đối với đất bão hòa đồng nhất, đẳng hướng thì phương trình thấm có thể viết như sau: 2h w 2h w (11) =0 + x 2 y2 Tóm tắt các phương trình vi phân dòng thấm trạng ổn định một hướng cho đất không bão hòa và bão hòa theo bảng sau: Bảng 3. Phương trình thấm ổn định hai hướng trong đất không bão hoà Không đồng nhất, dị hướng Không đồng nhất, đẳng hướng h h k h k h 2 hw h k h k h 2 2 2 w w wx w wy w w w w w w k + kwy + =0 k + =0 + . . + . . wx w x y x x y y x y x x y y 2 2 2 2 Bảng 4. Phương trình thấm ổn định hai hướng trong đất bão hoà 27
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 Dị hướng Đẳng hướng Không đồng nhất h h k h k h h h k h k h 2 2 2 2 k k +k + =0 + =0 w w sx w sy w w w s w s w + . . + . . x y x x y y x y x x y y sx sy s 2 2 2 2 Đồng nhất k sx h2w k sy h2w = 0 2hw 2hw = 0 2 2 + + x y x 2 y 2 Như vậy, có thể xem dòng thấm qua đất bão hòa là trường hợp đặc biệt của đất không bão hoà, lúc đó, hệ số thấm nước kw được lấy bằng ks. 3.4. Thiết lập phương trình cho dòng thấm ba hướng Xét một phân tố đất vô cùng nhỏ dx,dy,dz có dòng thấm đi qua. Vận tốc thấm vwx, vwy, vwz theo các hướng x,y,z là dương theo các hướng của trục vwy vwy + dy y vwz vwx y dy v wx v wx + dx dz x vwy x dx z vwz vwz + dz z Hình 5. Dòng thấm nước trạng thái ổn định ba hướng qua phần tử đất không bão hòa. Với dòng thấm ổn định có điều kiện liên tục nên: vwx vwy vwz vwx+ dx-vwx dydz+ vwy+ dy-vwy dxdz+ vwz+ dz-vwz dxdy = 0 x y z v w x + v w y + v w z d x d yd z = 0 Hay x y z Ap dụng định luật Darcy: kwx hw kwy hw kwz hw (12) =0 x x y y z z 2 hw 2 hw 2 hw kwx hw kwy hw kwz hw Hay: kwx 2 kwy 2 kwz 2 (13) . . . =0 x y z x x y y z z Trường hợp đất không đồng nhất, đẳng hướng, khi đó, phương dòng thấm có dạng: 28
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 2 hw 2 hw 2 hw kw hw kw hw kw hw 2 2 kw . . . =0 (14) x y z x x y y z z 2 Trường hợp đất bão hoà, không đồng nhất, dị hướng, phương trình dòng thấm có dạng: 2 hw 2 hw 2 hw ksx hw ksy hw ksz hw ksy 2 ksz 2 (15) ksx . . . =0 x 2 y z x x y y z z Trường hợp đất bão hoà, không đồng nhất, đẳng hướng , phương trình dòng thấm có dạng: 2 hw 2 hw 2 hw ks hw ks hw ks hw 2 2 ks . . . =0 (16) x 2 y z x x y y z z Trường hợp đất bão hoà, đồng nhất, dị hướng, phương trình vi phân dòng thấm có dạng: 2 hw 2 hw 2 hw ksy 2 +ksz 2 = 0 ksx x 2 y z Trường hợp đất bão hoà, đồng nhất, đẳng hướng, phương trình dòng thấm có dạng: 2 hw 2 hw 2 hw =0 + x 2 y 2 z 2 Tóm tắt các phương trình vi phân dòng thấm trạng ổn định ba hướng cho đất không bão hòa và bão hòa theo bảng sau: Bảng 5. Phương trình thấm ổn định ba hướng trong đất không bão hoà Không đồng nhất, dị hướng Không đồng nhất, đẳng hướng 2 hw h h h h k h k h k h h k h k h k h 2 2 2 2 2 w w w wx w wy w wz w w w w w w w w w k k k . . . =0 k . . . =0 wx wy wz w z x x x y z x x y y z z x y y y z z 2 2 2 2 2 2 Bảng 6. Phương trình thấm ổn định ba hướng trong đất bão hoà Dị hướng Đẳng hướng Không đồng nhất 2 hw h h h h k h k h k h h k h k h k h 2 2 2 2 2 w w w sx w sy w sz w w w s w s w s w k k k . . . =0 k . . . =0 sx sy sz s z x x x y z x x y y z z x y y y z z 2 2 2 2 2 2 Đồng nhất 2 hw 2 hw 2 hw 2 hw 2 hw 2 hw ksy 2 +ksz 2 = 0 ksx =0 + x 2 y z x 2 y 2 z 2 29
- TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ, ĐẠI HỌC ĐÀ NẴNG - SỐ 1(30).2009 4. Kết luận Trong việc đánh giá tính thấm của đất không bão hoà, có đưa đến khái niệm độ hút dính (matric suction): ua - uw. Trên cơ sở đó, áp dụng định luật thấm (Darcy) cho đất h . không bão hoà: v = kw(ua - uw) y Có thể xem hiện tượng thấm trong đất bão hòa là một trường hợp của hiện thấm trong đất không bão hoà, hay nói cách khác, thấm trong đất không bão hòa là trường hợp thấm tổng quát cho cả nền đất. TÀI LIỆU THAM KHẢO [1] Trần Văn Cúc, Cơ học chất lỏng, Nhà xuất bản Đại học Quốc gia Hà Nội, Hà Nội, 2004. [2] Phan Ngọc Cừ, Tôn Sĩ Kinh, Động lực nước dưới đất, Nhà xuất bản Đại học và Trung học chuyên nghiệp, Hà Nội, 1981. [3] C.W.Fetter, Địa chất thuỷ văn ứng dụng, T1&T2, Nhà xuất bản Giáo dục, Hà Nội, 2000 [4] D.G.Fredlund, H.Rahardjo, Cơ học đất không bão hòa, T1&T2, Nhà xuất bản Giáo dục, Hà nội, 1998 [5] Doãn Tam Hoè, Phương trình vi phân, Nhà xuất bản Giáo dục, Hà Nội, 2005. 30
CÓ THỂ BẠN MUỐN DOWNLOAD
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU CHẤT LƯỢNG NƯỚC VÀ TÔM TỰ NHIÊN TRONG CÁC MÔ HÌNH TÔM RỪNG Ở CÀ MAU"
12 p | 1363 | 120
-
Báo cáo nghiên cứu khoa học: "Cái tôi trữ tình trong thơ Nguyễn Quang Thiều."
10 p | 614 | 45
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU PHỐI TRỘN CHI TOSAN – GELATI N LÀM MÀNG BAO THỰC PHẨM BAO GÓI BẢO QUẢN PHI LÊ CÁ NGỪ ĐẠI DƯƠNG"
7 p | 518 | 45
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU THỰC NGHIỆM ẢNH HƯỞNG CỦA MƯA AXÍT LÊN TÔM SÚ (PENAEUS MONODON)"
5 p | 454 | 44
-
Báo cáo nghiên cứu khoa học: "ỨNG DỤNG PHƯƠNG PHÁP PCR-GENOTYPI NG (ORF94) TRONG NGHIÊN CỨU VI RÚT GÂY BỆNH ĐỐM TRẮNG TRÊN TÔM SÚ (Penaeus monodon)"
7 p | 378 | 35
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU ĐẶC ĐIỂM SINH HỌC DINH DƯỠNG CÁ ĐỐI (Liza subviridis)"
6 p | 378 | 31
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU ĐẶC ĐIỂM SINH HỌC SINH SẢN CỦA CÁ ĐỐI (Liza subviridis)"
8 p | 331 | 29
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU CẢI TIẾN HỆ THỐNG NUÔI KẾT HỢP LUÂN TRÙNG (Brachionus plicatilis) VỚI BỂ NƯỚC XANH"
11 p | 385 | 29
-
Báo cáo nghiên cứu khoa học: "Quan hệ giữa cấu trúc và ngữ nghĩa câu văn trong tập truyện ngắn “Đêm tái sinh” của tác giả Trần Thuỳ Mai"
10 p | 434 | 24
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU TẠO KHÁNG THỂ ĐƠN DÒNG VI-RÚT GÂY BỆNH HOẠI TỬ CƠ QUAN TẠO MÁU VÀ DƯỚI VỎ (IHHNV) Ở TÔM PENAEID"
6 p | 354 | 23
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU ƯƠNG GIỐNG VÀ NUÔI THƯƠNG PHẨM CÁ THÁT LÁT (Notopterus notopterus Pallas)"
7 p | 306 | 22
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU ĐẶC ĐIỂM SINH HỌC CÁ KẾT (Kryptopterus bleekeri GUNTHER, 1864)"
12 p | 298 | 20
-
Báo cáo nghiên cứu khoa học: "NGHIÊN CỨU DÙNG ARTEMIA ĐỂ HẠN CHẾ SỰ PHÁT TRIỂN CỦA TIÊM MAO TRÙNG (Ciliophora) TRONG HỆ THỐNG NUÔI LUÂN TRÙNG"
10 p | 367 | 18
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU PHÂN VÙNG THỦY VỰC DỰA VÀO QUẦN THỂ ĐỘNG VẬT ĐÁY"
6 p | 347 | 16
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU THIẾT LẬP HỆ THỐNG NUÔI KẾT HỢP LUÂN TRÙNG (Brachionus plicatilis) VỚI BỂ NƯỚC XANH"
10 p | 372 | 16
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU THAY THẾ THỨC ĂN SELCO BẰNG MEN BÁNH MÌ TRONG NUÔI LUÂN TRÙNG (Brachionus plicatilis) THÂM CANH"
10 p | 347 | 15
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU ƯƠNG GIỐNG CÁ KẾT (Micronema bleekeri) BẰNG CÁC LOẠI THỨC ĂN KHÁC NHAU"
9 p | 258 | 9
-
Báo cáo nghiên cứu khoa học: " NGHIÊN CỨU SỰ THÀNH THỤC TRONG AO VÀ KÍCH THÍCH CÁ CÒM (Chitala chitala) SINH SẢN"
8 p | 250 | 7
Chịu trách nhiệm nội dung:
Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA
LIÊN HỆ
Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM
Hotline: 093 303 0098
Email: support@tailieu.vn