Giáo trình phát triển và quản lý tài nguyên nước ngầm (dự án tăng cường năng lực đào tạo cho trường đại học thuỷ lợi của chính phủ đan mạch danida

TS NGUYỄN THU HIỀN (Chủ biên) TS HỒ VIỆT HÙNG - TS TRỊNH MINH THỤ PHÁT TRIỂN VÀ QUẢN LÝ TÀI NGUYÊN NƯỚC NGẦM HÀ NỘI - 2007 MỤC LỤC LỜI NÓI ĐẦU CHƯƠNG CƠ SỞ VẬN ĐỘNG VÀ TRUYỀN CHẤT TRONG NƯỚC NGẦM 1.1 Vận động dòng ngầm 1.1.1 1.1.2 1.1.3 1.1.4 nước 1.1.5 1.1.6 Các thành tạo địa chất chứa nước có áp, khơng áp bán áp Cột nước tổn thất cột nước nước ngầm Đặc trưng độ rỗng hệ số nhả nước 10 Định luật Đacxi, hệ số thấm, tính khơng đồng nhất, khơng đẳng hướng độ dẫn 11 Phương trình vi phân liên tục dòng ngầm - Điều kiện ban đầu điều kiện biên 14 Cơ sở vận động dòng ngầm toán đơn giản 16 1.2 Ô nhiễm nước ngầm trình di chuyển vật chất 29 1.2.1 1.2.2 1.2.3 1.2.4 Các nguồn ô nhiễm: tập trung không tập trung 29 Các trình truyền chất: đối lưu, phân tán, hấp thụ phân huỷ 30 Phương trình truyền chất: điều kiện biên điều kiên ban đầu 33 Trình bày tốn học tốn ô nhiễm tập áp dụng 36 CHƯƠNG VẬN ĐỘNG CỦA NƯỚC NGẦM TỚI LỖ KHOAN VÀ CÁC PHƯƠNG PHÁP XÁC ĐỊNH CÁC THÔNG SỐ CỦA TẦNG CHỨA NƯỚC 58 1.3 Dòng thấm hội tụ tới lỗ khoan 58 1.3.1 1.3.2 Vận động ổn định không ổn định tầng chứa nước đồng vô hạn 58 Ảnh hưởng lỗ khoan khơng hồn chỉnh 68 1.4 Dòng thấm tới lỗ khoan vùng gần biên chứa nước 70 1.4.1 1.4.2 1.4.3 Lỗ khoan gần sông 70 Dòng thấm đến lỗ khoan gần biên cách nước 73 Dòng thấm tới lỗ khoan gần biên khác 75 1.5 Sự tác dụng lẫn lỗ khoan 78 1.6 Khái niệm phân loại hút nước thí nghiệm 79 1.6.1 Sơ đồ phương pháp hút nước thí nghiệm xác định thông số ĐCTV tầng chứa nước 83 1.6.2 Sơ đồ chùm thí nghiệm xác định thơng số địa chất thủy văn chuyên môn 89 1.7 Phương pháp tiến hành thí nghiệm thu thập thơng tin 91 1.8 Hiệu suất lỗ khoan 93 1.8.1 1.8.2 Khái niệm 93 Các phương pháp xác định hiệu suất lỗ khoan 94 CHƯƠNG ĐIỀU TRA ĐÁNH GIÁ NƯỚC NGẦM 98 1.9 Tổng quan phương pháp đánh giá nước ngầm 98 1.9.1 1.9.2 1.9.3 Các phương pháp địa chất 98 Phương pháp viễn thám 99 Phương pháp thuỷ văn đồng vị 101 1.10 Phương pháp Địa vật lý 102 1.10.1 1.10.2 Phương pháp Địa vật lý mặt 102 Địa vật lý lỗ khoan 112 1.11 Phương pháp thuỷ văn 121 1.11.1 1.11.2 1.11.3 1.11.4 1.11.5 Thấm lượng bổ cập 121 Phương pháp cân nước 122 Các quan hệ mưa - dòng chảy 123 Các mơ hình dịng chảy ngầm khu vực 125 Nghiên cứu điển hình theo phương pháp cân nước 126 CHƯƠNG MƠ HÌNH TOÁN NƯỚC NGẦM 128 1.12 Tổng quan mơ hình hệ thống nước ngầm 128 1.12.1 1.12.2 1.12.3 1.12.4 1.12.5 1.12.6 1.12.7 1.12.8 Giới thiệu 128 Các bước tiến hành mô hình 129 Mơ hình khái niệm (conceptual model) 132 Thiết kế lưới mơ hình 134 Gán giá trị cho thông số, điều kiện ban đầu điều kiện biên 139 Chạy hiệu chỉnh mơ hình 145 Phân tích độ nhạy 155 Bài toán dự báo 156 1.13 Giới thiệu phần mềm MODFLOW 157 1.13.1 1.13.2 1.13.3 1.13.4 1.13.5 Sai phân hố phương trình cách giải 157 Xác định khoảng cách ô lưới theo phương đứng 158 Các loại mơ lớp mơ hình 159 Tính tốn thành phần hệ số thấm thẳng đứng VCONT 159 Các hệ số dẫn nước nút 160 1.14 Mơ hình lan truyền vật chất chất lượng nước ngầm 161 1.14.1 1.14.2 1.14.3 Đường phần tử vận chuyển đối lưu 161 Ấn định độ phân tán 162 Áp dụng 164 1.15 Thực hành mơ hình Modflow 166 CHƯƠNG QUẢN LÝ NƯỚC NGẦM .188 1.16 Sơ lược điều kiện địa chất thuỷ văn lãnh thổ Việt Nam 188 1.16.1 1.16.2 Phân vùng Địa chất thủy văn lãnh thổ Việt Nam 188 Trữ lượng nước ngầm 191 1.17 Khai thác nước mức vấn đề môi trường liên quan đến khai thác nước ngầm 191 1.17.1 1.17.2 1.17.3 1.17.4 Sụt lún mặt đất khai thác nước mức 191 Xâm nhập nước biển vào nước ngầm 193 Suy thoái chất lượng nước ngầm 202 Các ví dụ đánh giá sơ lún mặt đất xâm nhập mặn 204 1.18 Quản lý lưu vực 207 1.18.1 1.18.2 1.18.3 1.18.4 1.18.5 1.18.6 Khái niệm trữ lượng nước ngầm 207 Một số khái niệm lưu lượng nước ngầm 208 Đánh giá lưu lượng bền vững 211 Quản lý lưu vực việc sử dụng luân phiên nguồn nước 213 Các phương pháp bổ sung nhân tạo cho nước ngầm 218 Quản lý tầng chứa nước quan trắc động thái nước ngầm 222 1.19 Quản lý chất lượng nước ngầm 227 1.19.1 1.19.2 1.19.3 Bảo vệ nước ngầm khỏi bị ô nhiễm 227 Các phương pháp đánh giá ô nhiễm 230 Khắc phục ô nhiễm 234 PHỤ LỤC A 239 PHỤ LỤC B 240 PHỤ LỤC C 243 Phụ lục C .244 TÀI LIỆU THAM KHẢO 245 LỜI NÓI ĐẦU Nước ngầm nguồn nước lớn sẵn có trái đất ổn định có trữ lượng vượt xa so với nguồn nước mặt từ sông, suối, hồ, ao Hiện nay, phát triển ngành kinh tế vấn đề bùng nổ dân số, nhu cầu dùng nước tăng lên không ngừng mâu thuẫn khả cung cấp nước nhu cầu dùng nước ngày gay gắt số lượng chất lượng Vì vậy, nước ngầm lại trở nên gần gũi quan trọng sống người Tuy nhiên, viêc nghiên cứu vấn đề liên quan đến vận động nước ngầm đặc biệt việc phát triển quản lý nguồn tài nguyên nước ngầm quan tâm năm gần Để đáp ứng nhu cầu thiết thực cho việc nghiên cứu lĩnh vực này, Giáo trình “Phát triển quản lý tài nguyên nước ngầm” đề xuất khuôn khổ Tiểu hợp phần 1.3 “Hỗ trợ tăng cường lực cho Trường Đại học Thủy lợi” thuộc dự án hỗ trợ ngành nước (WaterSPS) DANIDA để đưa vào chương trình đào tạo đại học cao học ngành kỹ thuật thủy lợi, xây dựng, giao thông, mỏ địa chất Mục đích giáo trình cung cấp cho sinh viên kiến thức tổng quan qui luật vận động truyền tải vật chất nước dưói đất, thủy lực giếng cách xác định thơng số, ứng dụng mơ hình tốn nước ngầm, kỹ thuật phân tích, đánh giá trữ lượng chất lượng nước ngầm việc áp dụng chúng thực tế với thông tin cập nhật liên quan đến vấn đề phát triển quản lý nước ngầm Nội dung chủ yếu giáo trình gồm chương: Chương 1: Cơ sở dòng chảy truyền chất nước ngầm Chương 2: Vận động nước ngầm tới giếng khoan phương pháp xác định thông số tầng chứa nước Chương 3: Đánh giá trữ lượng nước ngầm Chương 4: Mô hình tốn nước ngầm Chương 5: Quản lý nước ngầm Đề cương giáo trình xây dựng với tư vấn phối hợp chuyên gia tư vấn dự án giảng viên Trường Đại học Thủy lợi biên soạn TS Nguyễn Thu Hiền chủ biên Chương 1, TS Nguyễn Thu Hiền viết, Chương TS Trịnh Minh Thụ viết Chương TS Hồ Việt Hùng viết Trong q trình biên soạn, chúng tơi tham khảo tài liệu ngồi nước chủ yếu tài liệu nước liên quan đến lĩnh vực thủy văn nước ngầm, địa chất thuỷ văn, nhiễm nước ngầm mơ hình tốn nước ngầm cập nhật với phương châm cố gắng giới thiệu nội dung cần thiết mới, tiếp cận vớí quốc tế thích ứng với điều kiện Việt Nam Chúng xin bày tỏ cám ơn đặc biệt tới TS Roger Chenevey - Cố vấn trưởng Tiểu hợp phần 1.3, tới GS.TS Gupta, chuyên gia tư vấn quốc tế xây dựng đề cương giáo trinh tới PGS TS Đoàn Văn Cánh, chuyên gia tư vấn nước xây dựng đề cương giúp đỡ, tư vấn nhiệt tình đặc biệt cung cấp thơng tin q giá thăm dò đánh giá trữ lượng nước ngầm Việt Nam ơng q trình viết giáo trình Chúng xin bày tỏ cám ơn chân thành tới PGS TS Phạm Quí Nhân, chuyên gia phản biện giáo trình với ý kiến đóng góp q báu ơng để đảm bảo chất lượng cho giáo trình Chúng tơi xin cám ơn Ban Giám hiệu Trường Đại học thuỷ lợi Văn phòng Dự án Hỗ trợ ngành nước (WaterSPS) DANIDA giúp đỡ chúng tơi q trình biên soạn Cuốn giáo trình xuất lần đầu thời gian có hạn nên khơng thể tránh khỏi sai sót chưa thực hồn chỉnh Chúng tơi mong nhận ý kiến phê bình đóng góp bạn đọc Mọi ý kiến xin gửi về: Bộ môn Thủy lực, Trường Đại học Thủy lợi , 175, Tây Sơn, Đống Đa, Hà Nội Chúng xin chân thành cám ơn CHƯƠNG 1: CƠ SỞ VẬN ĐỘNG VÀ TRUYỀN CHẤT TRONG NƯỚC NGẦM Sự vận động truyền chất nước ngầm liên quan đến trình vật lý hoá học xảy mặt đất mơi trường địa chất Chương trình bày sở quan trọng vận chuyển q trình truyền chất nước ngầm Đó sở để nghiên cứu vấn để liên quan đến phát triển quản lý tài nguyên nước ngầm 1.1 Vận động dòng ngầm 1.1.1 Các thành tạo địa chất chứa nước có áp, khơng áp bán áp Một thành tạo địa chất chứa lượng nước đáng kể có lỗ rỗng đủ lớn cho khai thác nước từ gọi tầng chứa nước Có nhiều thành tạo địa chất xem tầng chứa nước với khả trữ nước lỗ rỗng khác Lỗ rỗng hình thành đứt gãy, nứt nẻ xếp hạt đất đá Dưới vai trò số loại thành tạo địa chất chứa nước a) Bồi tích (phù sa) Có khoảng 90% tầng chứa nước thuộc loại Chúng bao gồm cuội, sỏi, cát bở rời Những thành tạo chứa nước phân làm bốn loại dựa hình thành nó: lịng sơng suối, thung lũng chơn vùi hay lịng sơng cổ, đồng thung lũng núi Loại thứ bao gồm bồi tích phù sa tạo nên lịng sơng bên cạnh bãi ven sơng Do nằm kề với dịng chảy mặt nên có lượng nước lớn thấm từ sơng ngịi vào đất Loại thứ hai thung lũng chôn vùi hay lịng sơng cổ dịng sơng thay đổi hướng chảy hình thành nên Mặc dù loại gần giống loại thứ nhất, độ thẩm, độ chứa, lượng bổ sung nước ngầm thường Loại thứ ba đồng rộng lớn bồi đắp phù sa Nằm đồng lớp cuội, sỏi cát tạo thành tầng chứa nước quan trọng Loại thứ tư thung lũng núi nơi có nhiều trầm tích chứa nước ngầm lớn Nguồn cung cấp nước chủ yếu nước mưa thấm từ dịng chảy khơng thường xun b) Đá vơi Đá vơi có mật độ, độ rỗng tính thấm nước thay đổi phạm vi lớn, tuỳ thuộc vào cấu tạo phát triển đới nứt nẻ, cáctơ hóa Những lỗ rỗng đá vơi lỗ nhỏ li ti, hang động lớn, hình thành nên dịng sơng ngầm Những mạch nước lớn thường tìm thấy vùng đá vơi Sự hoà tan CaCO3 nước làm cho nước ngầm vùng có độ cứng lớn Cũng hoà tan CaCO3 nước mà hang động, lỗ rỗng đá vôi ngày phát triển Hiện tượng gọi tượng cáctơ (karst) c) Đá hình thành núi lửa ( đá phun trào) Đá hình thành núi lửa tạo thành tầng chứa nước có tính thấm tốt, đặc biệt đá bazan Những lớp cuội sỏi cát vật liệu khác nằm xen kẽ hai lớp dung nham tạo cho đá bazan chứa thấm nước tốt Ngoài ra, khả chứa thấm nước tốt đá bazan cịn tượng phong hố vận động nội sinh gây d) Đá cát kết Đá cát kết đá dăm kết dạng bị xi măng hoá cát cuội sỏi Do vậy, độ rỗng khả sinh nước ngầm chúng bị giảm liên kết xi măng Các tầng chứa nước đá cát kết chứa nước ngầm khe nứt, lỗ hổng song song tồn e) Đá magma biến chất Các dạng đá magma biến chất thường rắn chắc, nứt nẻ nên chúng thường tầng chứa nước Ở nơi loại đá lộ mặt đất, chúng bị phong hoá mạnh phát triển thành tầng chứa nước Lượng nước chứa loại thành tạo tương đối nhỏ đủ dùng cho sinh hoạt số hộ g) Đất sét Đất sét nói chung có độ rỗng tương đối lớn lỗ hổng chúng lại nhỏ đến mức coi chúng vật liệu khơng thấm nước Các tầng đất sét nằm hệ chứa nước tốt hình thành nên thấu kính nước ngầm cục hình thành nên tầng chứa nước bán áp Hình 1.1 Sơ đồ mơ tả loại tầng chứa nước Tầng chứa nước phân chia thành loại sau: Tầng chứa nước khơng áp tầng chứa nước mực nước ngầm mặt tầng bão hoà Mực nước ngầm biến đổi phụ thuộc vào diện tích miền cung cấp nước ngầm, trình khai thác nước ngầm khả thấm tầng chứa nước Tầng chứa nước Hình 1.1 tầng chứa nước không áp Sự dao động mực nước ngầm tương ứng với thay đổi lượng nước tàng trữ tầng chứa nước Để xây dựng đồ mực nước ngầm, ta dựa vào số liệu điều tra giếng vùng Bản đồ đẳng bề mặt nước ngầm gọi đồ thủy đẳng cao Trường hợp đặc biệt tầng chứa nước khơng áp nước thượng tầng (Hình 1.1) Nước thượng tầng thường có vùng trầm tích bở rời, phía lớp cách nước Giếng khoan gặp phải nước thượng tầng thường cung cấp lưu lượng nhỏ mang tính chất tạm thời Tầng chứa nước có áp tầng chứa nước giới hạn tầng không thấm nước áp suất lớn áp suất khí Ở giếng khoan tầng chứa nước có áp, mực nước giếng dâng cao mái cách nước tầng chứa nước Đặc biệt, số trường hợp mực nước cao mặt đất tạo giếng phun Hình 1.1 Miền cung cấp nước cho tầng chứa nước có áp gọi miền cung cấp Sự thay đổi mực nước giếng có áp phụ thuộc chủ yếu vào thay đổi cột nước áp suất Vì thế, coi đường ống dẫn để chuyển nước từ vùng cấp đến vùng thoát Đường thủy áp đường tưởng tượng trùng với đường cột nước thủy tĩnh tầng chứa nước Khi mực nước có áp hạ thấp mái cách nước, trở thành tầng chứa nước khơng áp Tầng chứa nước bán áp tương tự tầng chứa nước có áp, mái có khả thấm xuyên Nước tầng bán áp trao đổi với bên ngồi tùy vào vị trí mực nước ngầm cột nước áp suất tầng chứa nước 1.1.2 Cột nước tổn thất cột nước nước ngầm Dòng chảy ổn định chất lỏng khơng nén tn theo phương trình Becnuli phát biểu tổng cột nước điểm dịng chảy liên tục khơng đổi (1.1) p v2 z+ + = const γ 2g đó: z vị trí điểm nghiên cứu so với mặt chuẩn (m), p áp suất (N/m ), v vận tốc dòng chảy (m/s), γ trọng lượng riêng chất lỏng (N/m3), g gia tốc trọng trường (m/s2) Đối với chất lỏng thực, chất lỏng có tính nhớt sinh số tổn thất cột nước dọc theo dòng chảy Gọi tổn thất cột nước mặt cắt hL (Hình 1.2), phương trình Becnuli hai mặt cắt viết sau: 2 (1.2) p v p V z1 + + = z + + + hL 2g γ 2g γ Tuy nhiên, với trường hợp dịng chảy mơi trường lỗ rỗng, vận tốc dòng chảy thường nhỏ cột nước lưu tốc (v2/2g) bỏ qua Phương trình (1.2) viết thành: p p (1.3) z1 + = z + + hL γ γ cột nước h điểm dòng chảy là: p (1.4) h=Z+ γ Hình 1.2 Phân bố cột nước áp suất tổn thất cột nước qua cột thấm 1.1.3 Đặc trưng độ rỗng hệ số nhả nước Độ rỗng phần thể tích tạo khe hở lỗ rỗng Độ rỗng thường biểu thị theo phần trăm sau: W (1.5) n = w 100% W đó: n độ rỗng, Ww thể tích lỗ rỗng, W tổng thể tích mẫu đất đá Trong vật liệu rời, độ rỗng phụ thuộc vào ba tính chất đất đá: độ nén chặt, hình dạng hạt phân bố kích thước hạt Độ rỗng loại đất đá khác trình bày Bảng 1.1 Hệ số nhả nước trọng lực (Sy) (specific yield) Hệ số nhả nước trọng lực đất đá tỉ số lượng nước (trong đới bão hoà) trọng lực thể tích nó: Wy (1.6) Sy = W Trong Wy thể tích nước Bảng 1.1 Độ rỗng loại đất đá khác (theo Todd Mays 2005) Vật liệu Sỏi thơ Sỏi trung bình Sỏi mịn Cát thơ Cát trung bình Cát mịn Đất phù sa Sét Cát kết hạt mịn Cát kết hạt trung bình Đá vơi Đolomit Cát cồn cát ven biển Độ rỗng (%) 28 32 34 39 39 43 46 42 33 37 30 26 45 Vật liệu Hoàng thổ Than bùn Đá phiến (schist) Bột kết Sét kết Đá phiến sét (shale) Tảng lăn lẫn sét, bột Tảng lăn lẫn cát Tro núi lửa (tuff) Đá Bazan Gabrô bị phong hóa Granit bị phong hố Độ rỗng (%) 49 92 38 35 43 34 31 41 17 43 45 Hệ số giữ nước (Sr) Hệ số giữ nước đất đá tỉ số lượng nước cịn giữ lại sau nước trọng lực thể tích nó: 10 thành chúng thấm qua mơi trường thấm yếu có độ hấp thụ cao chiều dày đủ lớn Việc phục hồi chất lượng nước đất phức tạp tốn từ phải loại trừ tất nguồn gây bẩn CÂU HỎI CUỐI CHƯƠNG Trình bày vấn đề mơi trường liên quan đến khai thác nước ngầm cách mức? Trình bày việc quản lý lưu vực việc sử dụng luân phiên nguồn nước? Trình bày phương pháp bổ sung nhân tạo cho nước ngầm? Trình bày việc quản lý tầng chứa nước quan trắc động thái nước ngầm? Trình bày biện pháp bảo vệ nước ngầm khỏi ô nhiễm? Các phương pháp đánh giá ô nhiễm? Các biện pháp khắc phục ô nhiễm nước ngầm? 238 PHỤ LỤC A Các giá trị hàm sai số Giá trị hàm sai số erf(x) hàm sai số bù erfc(x) cho giá trị dương x x erf(x) erfc(x) x erf(x) erfc(x) 0.00 1.1 0.880205 0.119795 0.05 0.056372 0.943628 1.2 0.910314 0.089686 0.10 0.112463 0.887537 1.3 0.934008 0.065992 0.15 0.167996 0.832004 1.4 0.952285 0.047715 0.20 0.222703 0.777297 1.5 0.966105 0.033895 0.25 0.276326 0.723674 1.6 0.976348 0.023652 0.30 0.328627 0.671373 1.7 0.983790 0.016210 0.35 0.379382 0.620618 1.8 0.989091 0.010909 0.40 0.428392 0.571608 1.9 0.992790 0.007210 0.45 0.475482 0.524518 2.0 0.995322 0.004678 0.50 0.520500 0.479500 2.1 0.997021 0.002979 0.55 0.563323 0.436677 2.2 0.998137 0.001863 0.60 0.603856 0.396144 2.3 0.998857 0.001143 0.65 0.642029 0.357971 2.4 0.999311 0.000689 0.70 0.677801 0.322199 2.5 0.999593 0.000407 0.75 0.711156 0.288844 2.6 0.999764 0.000236 0.80 0.742101 0.257899 2.7 0.999866 0.000134 0.85 0.770668 0.229332 2.8 0.999925 0.000075 0.90 0.796908 0.203092 2.9 0.999959 0.000041 0.95 0.820891 0.179109 3.0 0.999978 0.000022 1.00 0.842701 0.157299 239 PHỤ LỤC B Hàm số Bessel Hàm số biến đổi Bessel loại hai bậc K0(x) x 0.010 0.011 0.012 0.013 0.014 0.015 0.016 0.017 0.018 0.019 0.020 0.021 0.022 0.023 0.024 0.025 0.026 0.027 0.028 0.029 0.030 0.031 0.032 0.033 0.034 0.035 0.036 0.037 0.038 0.039 K0(x) 4.721 4.626 4.539 4.459 4.385 4.316 4.251 4.191 4.134 4.080 4.028 3.980 3.933 3.889 3.846 3.806 3.766 3.729 3.692 3.657 3.623 3.591 3.559 3.528 3.499 3.470 3.442 3.414 3.388 3.362 x 0.040 0.041 0.042 0.043 0.044 0.045 0.046 0.047 0.048 0.049 0.050 0.051 0.052 0.053 0.054 0.055 0.056 0.057 0.058 0.059 0.060 0.061 0.062 0.063 0.064 0.065 0.066 0.067 0.068 0.069 K0(x) 3.336 3.312 3.288 3.264 3.241 3.219 3.197 3.176 3.155 3.134 3.114 3.094 3.075 3.056 3.038 3.019 3.001 2.984 2.967 2.950 2.933 2.916 2.900 2.884 2.869 2.853 2.838 2.823 2.809 2.794 x 0.070 0.071 0.072 0.073 0.074 0.075 0.076 0.077 0.078 0.079 0.080 0.081 0.082 0.083 0.084 0.085 0.086 0.087 0.088 0.089 0.090 0.091 0.092 0.093 0.094 0.095 0.096 0.097 0.098 0.099 K0(x) 2.780 2.766 2.752 2.738 2.725 2.711 2.698 2.685 2.673 2.660 2.647 2.635 2.623 2.611 2.599 2.587 2.576 2.564 2.553 2.542 2.531 2.520 2.509 2.499 2.488 2.478 2.467 2.457 2.447 2.437 240 Phụ lục B: Hàm số biến đổi Bessel loại hai bậc K0(x) (tiếp theo) x 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37 0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48 0.49 0.5 0.51 0.52 0.53 K0(x) 2.427 2.333 2.248 2.169 2.097 2.030 1.967 1.909 1.854 1.802 1.753 1.706 1.662 1.620 1.580 1.541 1.505 1.470 1.436 1.404 1.372 1.342 1.314 1.286 1.259 1.233 1.207 1.183 1.160 1.137 1.114 1.093 1.072 1.052 1.032 1.013 0.994 0.976 0.958 0.941 0.924 0.908 0.892 0.877 x 0.55 0.56 0.57 0.58 0.59 0.60 0.61 0.62 0.63 0.64 0.65 0.66 0.67 0.68 0.69 0.70 0.71 0.72 0.73 0.74 0.75 0.76 0.77 0.78 0.79 0.80 0.81 0.82 0.83 0.84 0.85 0.86 0.87 0.88 0.89 0.90 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 K0(x) 0.847 0.832 0.818 0.804 0.791 0.777 0.765 0.752 0.740 0.728 0.716 0.704 0.693 0.682 0.671 0.660 0.650 0.640 0.630 0.620 0.611 0.601 0.592 0.583 0.574 0.565 0.557 0.548 0.540 0.532 0.524 0.516 0.509 0.501 0.494 0.487 0.480 0.473 0.466 0.459 0.452 0.446 0.440 0.433 x 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5.0 K0(x) 0.421 0.366 0.318 0.278 0.244 0.214 0.188 0.165 0.146 0.129 0.114 0.101 0.0893 0.0791 0.0702 0.0623 0.0554 0.0493 0.0438 0.0390 0.0347 0.0310 0.0276 0.0246 0.0220 0.0196 0.0175 0.0156 0.0140 0.0125 0.0112 0.0100 0.0089 0.0080 0.0071 0.0064 0.0057 0.0051 0.0046 0.0041 0.0037 241 0.54 0.861 0.99 0.427 Phụ lục B: Hàm số biến đổi Bessel loại hai bậc K0(x) (tiếp theo) x 5.0 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 6.0 6.1 6.2 6.3 6.4 6.5 6.6 6.7 6.8 6.9 7.0 7.1 7.2 7.3 74.0 7.5 K0(x) 3.69 x 10-3 3.31 x 10-3 3.00 x 10-3 2.66 x 10-3 2.39 x 10-3 2.14 x 10-3 1.92 x 10-3 1.72 x 10-3 1.54 x 10-3 1.39 x 10-3 1.24 x 10-3 1.12 x 10-3 1.00 x 10-3 9.00 x 10-4 8.08 x 10-4 7.26 x 10-4 6.52 x 10-4 5.86 x 10-4 5.26 x 10-4 4.73 x 10-4 4.25 x 10-4 3.82 x 10-4 343 x 10-4 3.08 x 10-4 2.77 x 10-4 2.49 x 10-4 x 7.60 7.70 7.80 7.90 8.00 8.10 8.20 8.30 8.40 8.50 8.60 8.70 8.80 8.90 9.00 9.10 9.20 9.30 9.40 9.50 9.60 9.70 9.80 9.90 10.00 K0(x) 2.24 x 10-4 2.01x10-4 1.81 x 10-4 1.63 x 10-4 1.46 x 10-4 1.32 x 10-4 1.18 x 10-4 1.07 x 10-4 9.59 x10-5 8.63 x 10-5 7.76 x 10-5 6.98 x10-5 6.28 x10-5 5.65 x 10-5 5.09 x 10-5 4.58 x 10-5 4.12 x 10-5 3.71 x 10-5 3.34 x10-5 3.01 x 10-5 2.71 x 10-5 2.44 x10-5 2.19 x10-5 1.97 x 10-5 1.78 x 10-5 242 PHỤ LỤC C Giá trị hàm W(tD,B) Giá trị hàm W(tD,B) với giá trị khác tD B tD 0.000002 0.000004 0.000006 0.000008 0.00001 0.00002 0.00004 0.00006 0.00008 0.0001 0.0002 0.0004 0.0006 0.0008 0.001 0.002 0.004 0.006 0.008 0.01 0.02 0.04 0.06 0.08 0.1 0.2 0.4 0.6 0.8 0.002 0.004 0.006 0.008 0.01 0.02 0.04 0.06 0.08 12.70 12.10 11.60 11.30 11.00 10.80 10.20 9.52 9.13 8.84 8.62 7.94 7.24 6.84 6.55 6.33 5.64 4.95 4.54 4.26 4.04 3.35 2.68 2.30 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 11.30 11.20 11.10 10.90 10.70 10.60 10.10 9.45 9.08 8.81 8.59 7.92 7.24 6.84 6.55 6.33 5.64 4.95 4.54 4.26 4.04 3.35 2.68 2.30 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 10.50 10.50 10.40 10.40 10.30 10.20 9.84 9.34 9.00 8.75 8.55 7.90 7.22 6.83 6.54 6.32 5.63 4.95 4.54 4.26 4.04 3.35 2.68 2.30 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 9.89 9.89 9.88 9.87 9.84 9.80 9.58 9.19 8.89 8.67 8.48 7.86 7.21 6.82 6.53 6.32 5.63 4.94 4.54 4.26 4.04 3.35 2.68 2.30 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 9.44 9.44 9.44 9.44 9.43 9.42 9.30 9.01 8.77 8.57 8.40 7.82 7.19 6.80 6.52 6.31 5.63 4.94 4.54 4.26 4.04 3.35 2.68 2.30 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 8.06 8.06 8.06 8.06 8.06 8.06 8.06 8.03 7.98 7.91 7.84 7.50 7.01 6.68 6.43 6.23 5.59 4.92 4.53 4.25 4.03 3.35 2.68 2.29 2.03 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 6.67 6.67 6.67 6.67 6.67 6.67 6.67 6.67 6.67 6.67 6.67 6.62 6.45 6.27 6.11 5.97 5.45 4.85 4.48 4.21 4.00 3.34 2.67 2.29 2.02 1.82 1.22 0.702 0.454 0.311 0.219 0.049 0.0038 0.0004 5.87 5.87 5.87 0.87 5.87 5.87 5.87 5.87 5.87 5.87 5.87 5.86 5.83 5.77 5.69 5.61 5.24 4.74 4.41 4.15 3.95 3.31 2.66 2.28 2.02 1.82 1.22 0.702 0.454 0.31 0.219 0.049 0.0038 0.0004 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.29 5.27 5.25 5.21 4.98 4.59 4.30 4.08 3.89 3.28 2.65 2.27 2.01 1.81 1.22 0.701 0.454 0.31 0.219 0.049 0.0038 0.0004 243 Phụ lục C: Giá trị hàm W(tD,B) với giá trị khác tD (tiếp theo) B tD 0.000002 0.000004 0.000006 0.000008 0.00001 0.00002 0.00004 0.00006 0.00008 0.0001 0.0002 0.0004 0.0006 0.0008 0.001 0.002 0.004 0.006 0.008 0.01 0.02 0.04 0.06 0.08 0.1 0.2 0.4 0.6 0.8 0.1 0.2 0.4 0.6 0.8 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.85 4.84 4.83 4.71 4.42 4.18 3.98 3.81 3.24 2.63 2.26 2.00 1.80 1.22 0.70 0.453 0.31 0.219 0.049 0.0038 0.0004 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.51 3.50 3.48 3.43 3.36 3.29 2.95 2.48 2.17 1.94 1.75 1.19 0.693 0.45 0.308 0.218 0.049 0.0038 0.0004 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.23 2.18 2.02 1.85 1.69 1.56 1.11 0.665 0.436 0.301 0.213 0.048 0.0038 0.0004 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.55 1.52 1.46 1.39 1.31 0.996 0.621 0.415 0.289 0.206 0.047 0.0037 0.0004 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.11 1.08 1.05 0.857 0.565 0.387 0.273 0.197 0.046 0.0037 0.0004 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.842 0.228 0.839 0.228 0.832 0.228 0.819 0.228 0.715 0.227 0.502 0.210 0.354 0.177 0.254 0.144 0.185 0.114 0.044 0.034 0.0036 0.0031 0.0004 0.0003 0 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0223 0.0222 0.0218 0.0207 0.011 0.0016 0.0002 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0025 0.0003 0.0021 0.0003 0.0006 0.0002 0.0001 0 244 TÀI LIỆU THAM KHẢO Anderson, M.P 1979 Using models to simulate the movement of contaminants through groundwater flow systems CRC critical reviews in environmental controls vol 9, no 2, 97-156 Anderson, M.P W.W Woesner 2002 Applied groundwater modeling, simulation of flow and avection transport Acedemic Press, California, USA Bear, J.1961 Some experiments on dispersion Journals of Geophysical Research vol 66 no.6, 2455-2467 Bear, J 1972 Dynamics of fluids in porous media New York: American Elsevier Publishing Company, 764 p Blandford, T.N P.S Huỵakorn, 1990 WHPA: a modular semi-analytical model for the delineation of wellhead protection areas, U.S EPA, Office of Ground Water Protection Borden, R.C P.B Bedient 1986 Transport of dissolved hydrocarbons influenced by oxygen-limited biodegredation, Theory development, Water Resources Research 22 no 13, 1983-1990 Boulton, N S 1954 The drawdown of the water-table under non-steady conditions near a pumped well in an unconfined formation, Proc Inst Civil Engrs, v.3, 564579 Boulton, N S 1975 Analysis of data from non-equilibrium pumping tests allowing for delayed yield from storage, Proc Inst Civil Engrs, v.26, 469-482 Boulton, N S., and T.D Streltsova 1975 New equations for determining the formation constants of an aquifer from pumping test data, Water Resources Research, v.11, 148-153 Bouwer, E J P.L McCarty 1984 Modeling of trace organics biotransformation in the subsurface Ground water 22, no 4, 433-440 Carrera, J., 1988 State of the art of the inverse problwm applied to the flow and solute transport equtions, in: Groundwater flow and quality modellling )E Custodio et al eds.) D Reidel Publ Co., 549-583 Carrera, J., G.R Walter, M Kuhn, H Bentley J Fabryka-Martin, 1984 Threedimesional modeling of saline pond leakage calibrated by INVERT-3, a quasithree-dimensional, transient, parameter estimation program, in Pratical applications of ground water models, Nat’l Water Well Assc., Columbus, Ohio, 547-569 Chapelle, F H 1993 Groundwater microbiology and Geochemistry John Wiley and Sons, Inc New York, 424p Cooley, R.L., 1977 A method of estimating parameters and assessing reliability for models of steady state ground water flow Theory and numerical properties, Water resources reseach 13(2), 318-324 245 Cooley, R.L., 1979 A method of esitmating parmeter and assessing reliabitity for modeld of steady state ground water flow, Water resources reseach 15(3), 603617 Cooley, R.L.,1982 Incorporation of prior information on parameters into non-linear regression groundwater flow models Theory, Water resources reseach 18(4), 965-976 Copper, H H., Jr., and C.E Jacob 1946 A generalized graphical method for valuating formation constants and summarizing well-field history, Trans Mer Geophysical Union, v.27, 526-534 Darcy, H 1856 Les fontaines publicque de la ville de Dijon, V Dalmont, Paris, 647p Davies-Smith, A., E L Bolke, and C A Collins, 1988, Geohydrology and digital simulation of the ground-water flow system in the Umatilla Plateau and Horse Heaven Hills area, Oregon and Washington, USGS, Water-Resources Investigations Report 87-4268, 72 p De Josselin De Jong, G 1958 Longitudinal and transverse diffusion in granular deposits Transations, American Geophysical Union39 no 1: 67p Đoàn Văn Cánh Phạm Quý Nhân 2005 Tin học Địa chất thuỷ văn ứng dụng NXB Khoa học Kỹ thuật 220 trang Đồn Văn Cánh Phạm Q Nhân 2003 Tìm kiếm thăm dò đánh giá trữ lượng nước đất NXB Xây dựng Faust, C.R J.W Mercer 1980 Ground water modeling: Recent developments Ground water 18, no 6: 569-577 Faust, C R., P N Sims, C.P Spalding P.F Andersen, 1990 FTWORK: a threedimensional groundwater flow and solute transport code (U), Westinghouse Savannah River Co., WSRC-RP-89-1085, 173p Fenemore, A.D.1988 A three-dimensional model for management of the Waimea Plains Aquifer, Nelson Dept of Science and Industrial Research, Hydrology Centre Publication No 18, 33 Chrischurch, New Zealand, Fetter, C W 1977 Attenuation of waste water elutriated through glacial outwash Ground water 15, no 5, 365-371 Fetter, C.W., 1988 Applied Hydrogeology (2nd edition), Merrill Publishing Co., Columbus, Ohio, 592 p Fetter C W 1999 Contaminant hydrogeology Prentice-Hall, Inc Upper Saddle River, New Jersey, USA, p500 Franz, T N Guiguer, 1990 FLOWPATH, two-dimensional horizontal aquifer simulation model, Waterloo hydrogeologic software, Waterloo, Ontario, 74p Freeze, R A.1969 Theoretical analysis of regional ground water flow Scientific Ser., No Inland Waters Branch, Dept Energy, mine and resources, Ottawa, Canada, 202p Fried, J.J 1975 Groundwater Pollution Amsterdam: Elsevier Scientific Publishing Company 246 Gelhar, L.W., and M A Collins 1971 Gerneral analysis of longitudinal dispersion in uniform flow Water Resources Research no 6,1511-1521 GeoTrans, Incorporated, 1988 SELTRAN: A simple and efficient two-dimensional groundwater flow and transport model, version 2.7, Herndon, VA, 142p Gerhart, S M., and G J Lazorchick, 1988 Evaluation of the ground-water resourcesof the Lower Susquehanna River Basin, Pennsylvania and Maryland, USGS, Water Supply Paper 2284, 128 p Gillham, R W., E A Sudicky, J A Cherry, and E O Frind 1984 An advectiondffusion conceptfor solute transport in heterogeneous unconsolidated geological deposits Water Resources Research 20 no 3, 369-378 Goode, D.J L.F Konikow 1990 Apparent dispersion in transient groundwater flow, Resources Research vol 26, no 10, 2339-2352 Gupta S.K., C.R.Cole, C.T Kincaid A.M Monti, 1987, Couple fluid energy and solute transport (CFEST) model: formulation and user’s manual, BMI/ONWI-660, Office of NuclearWaste Isolation, Battelle Memorial Institute, Columbus, Ohio, 466p Hantush, M S 1956 Analysis of data from pumping tests in leaky aquifers, trans Amer Geophysical Union, v.37, 702-714 Hantush, M S 1964 Drawdown around wells of variable discharge, jour Geophysical research, v.69, 4221-4235 Hantush, M S 1965 Wells near streams with semipervious beds, jour Geophysical research, v.70, 2829-2838 Hantush, M S 1967 Flow to wells in aquifers separated by a semipervious layer Jour Geophysical Research, v.72, 1709-1720 Hantush, M S and Jacob 1955 Non-steady radial flow in an infinite leaky aquifer, trans Amer Geophysical Union, v.36, 95-112 Harvey, R.W., R.L Smith L George 1984 Effect of organic contamination upon microbial distribution and heterotropic uptake in a Cape Cod, Mass., aquifer Applied and environmental microbiology vol 48, 1197-1202 Hill M.C., 1990 MODFLOWP: A computer profgram for estimating parameters of a transient, three-dimensional ground-water flow equations USGS, WaterResources Investgations Report 90-4048, 43p Hollett, K J Geohydrology and Water Resources of the Papago Farms-great Plain Area, Papago Indian Reservation, Arizona and the Upper Rio Sonoyta Area, Sonoram Mexico, US Geological Survey Water Supply Paper 2258, 1985 Isaaks, E H Srivastave, R M 1989 An introduction to applied geostatistics, Oxford University Press, 561p New York, USA Jacob, C E 1946 Radial flow in a leaky artesian aquifer, Trans Amer Geophysical Union, v.27, 198-208 Keys, W.S., 1989 Borehole geophysics applied to ground-water investigations, National Water Well Association, Dublin, OH 247 Keys, W.S & L.M MacCary, 1971 Application of borehole geophysics waterresources investigations, Techniques of Water-resources Investigations Reports, U.S Geological techniques of Water-resources Investigation, Bk 2, Chap E1, http://water.usgs.gov/pubs/twri/twri2-e1, 126p Kipp, K.L., Jr., 1987 HST3D: a computer code for simulatiion of heat and solute transport in three-dimensional ground-water flow systems, USGS, Water Resources Investigations Report 86-4095, 517p Klotz, D., K P Seiler., H Moser, and F Newmaier 1980 Dispersivity and velocity relationshipfrom laboratory and field relationships Journal of Hydrology 45 no 3, 169-84 Konikow, L F J.D Bredehoeft, 1978 Computer model of two dimensional solute transport and dispersion in ground water, Techniques of Water reources Investigations Book 7, Chap C2, USGS, 90p LeBlanc, D R., S P Garabedian, K M Hess, L W Gelhard, R D Quadri, K L Stollenwerk, W W Wood, 1991 Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts, Experimental design and observed tracer movement Water Resources Research 27, no 5, 895-910 Linehan, D S Keith, 1949 Seismic reconnaission for groundwater development, Journal of New England Water Works Assc Vol 63, 76-95 Luckey ,R R., E.D Gutentag, F J Heimes, and J B Weeks, 1986, Digital simulation of ground-water flow in the High Plains aquifer in parts of Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming, USGS, Professional Paper 1400-D, 57p Mann, J.F., Jr 1958 Estimating quantity and quality of groundwater in dry regions using airphotos, Inter Assoc Sci Hydrology Publ 44, 125-134 Miller, C T W J Weber 1984 Modeling organic contamination partitioning in ground-water systems Ground water 22, no 5, 584-92 Mackay , D M., D L Freyberg, P V Roberts, and J A Cherry 1986 A natural gradient experiment on solute transport in a sand aquifer, Approach and overview of plume movement Water Resources Research 22, no 13, 2017-2029 McCarty, P.L., M Reinhard B.E Rittman, 1981 Trace organic in groundwater Environmental science and technology 15, no 1, 40-51 McDonald, M G Harbaugh A W 1988 A modular three dimensional finitedifference groundwater flow model, Technicques of water resources investigations 06-A1, USGS, 576p Mercer, J.W., S.D Thomas, and B Ross, 1982, Parameters and variables appearing in repository siting models, NUREG/CR-3066, U.S Nuclear Regulatory Commission, 244 p Meyer, W.R., 1963 Use of neutron probe to determine the storage coeffiecient of an unconfined aquifer, U.S Geological Servey Professional Paper 450E, E174-E176 Moore, J E 1979 Contributions of groundwater modeling to planning, Journal of hydrology, No 43, 121-128 248 Munns, D N R.L Fox, 1976 The slow reactions that continue after phosphat adsorption: kinetics and equilibrium in tropical soils Soil science Society of America Proceedings 40, 46-51 Neuman, S P 1972 Theory of flow in unconfined aquifers considering delayed response of the water table, Water Resources Research, v.8, 1031-1045 Neuman, S.P., G.R Fogg, and E.A Jacobson, 1980 A statistical approach to the inverse problem of aquifer hydrology Case study, resources reseach Vol 16, No 1, 32-58 Nguyễn Kim Cương 1991 Địa chất thuỷ văn (dùng cho ngành Thuỷ văn) NXB Khoa học Kỹ thuật 268 trang Norris, S.E., 1972 The use of gamma logs in determining the character of unconsolidated sediments and well construction features, Groundwater, v 10, no 6: 14-21 Ophori, D.U J Toth, 1989, Characterization of ground-water flow by field mapping and numerical simulatio, Ross Creek Basin, Alberta, Canada, Ground Water 27(2), 193-201 Picken, John F., and Gerald E Grisak 1981 Scale-dependent dispersion in a stratified Granular aquifer Water Resources Research 17, no 4, 1191-1211 Pinder, George 1973 A Galerkin-finite element simulation of groundwater contamination on Long Island, New York, water Resources Research 9, no 6, 1657-69 Pickens, J F Grisak E G (1981) Scale-depentdent dispersion in a stratified granular aquifer Water Resources Research 17 (4), 1191-1211 Pollock, D.W., 1988 Semianalytical computation of path liens for finite-difference models, Ground Water 16(6), 743-750 Pollock, D.W., 1989 Documentation of computer programs to complete and display pathlines using results from the U.S.G.S Modular three-dimensional finitedifference ground-water model, USGS, Open file report 89-381, 81p Prickett, T.A., T.G Naymik C.G Lonnquist, 1981 A “Random-warlk” solute transport model for selected groundwater quality evaluations Illinois State Water Survey, Bulletin 65, 103p Rice, W A Gorelick S M., Geological inference from “flownet” transmissivity determination: three case studies, Water resources Bulletin, Vol 21, No.6, 919930, 1985 Rorabaugh, M I 1953 Graphical and theoretical analysis of step drawdown test of artesian well, Proc Amer Soc Civil Engrs., v79, 23 p Rumbaugh,J.O., III, and G.M Duffield, 1989 Model Cad TM: Computer-aided design for ground-water modeling, Version 1.0, Geraghty&Miller, Inc., Reston, VA., 52p Shafer, J.M.,1987 Reserve pathline calculation of time related capture zone in nonuniform flow, Ground water 25(3), 283-289 Shafer, J.M., 1990, GWPATH, Version 4.0, J.M Shafer, Champaign, IL 249 Sampler, J., J Carrera, G Galarza, and A Medina, 1990, Application of an automatic calibration technique to modelling an alluvial aquifer, In: Calibration and Reliability in Groundwater Modelling (K Kovar, ed.), IAHS Publ 195, 87-96 Sauty, J P (1978) Identification des parametres du transport hydrodispersif das les aquiferes par interpretation de tracages en ecoulement cyclndrique convergent ou divergent Journal of Hydrology 39, no 3/4,69-103 Sauty, J P 1980 An analysis of hydrodispersive transfer in aquifers Water resources Research 16, no 1, 145-158 Silliman, S.E A.L.Wright, 1988 Stochastic analysis of paths of high hydraulic conductivity in porous media, Water Resources Research 24(11), 1901-1910 Srinivasan, P J.W Mercer 1988 Simulation of biodegradation and sorption processes in ground water Journal of environmental quality 6, 438-443 Sudicky, E A., and J K Cherry 1979 Field observations of tracer dispersion under natural flow conditions in an unconfined sandy aquifer Water Pollution Research, Canada 14, 1-17 Theis, C V 1935 The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage, Trans Amer Geophysical Union, v 16, 519-524 Thiem, G 1906 Hydrologische Methoden, Gehardt, Leipzig, 56p Todd, D K J Bear 1961 Seepage through layered anisotropic porous media Journal of hydraulic division, ASCE, vol 87, no HY3, 31-57 Todd D K Mays L W 2005 Ground water hydrology John Willey and Sones, New York Toth J 1963 A theoretical analysis of groundwater flow in small drainage basins Journal of Geophysical Research, vol 68, 4795-4812 Townley L.R., M.R Davison, 1988, Definition of a capture zone for shalow water table lakes, Journal of Hydrology 104, 53-76 Van Genuchten, M.Th 1981 Analytical solutions for chemical transport with simultaneous adsorption, zero-order production, and first-order decay Journal of Hydrology 49, 213-233 Walton, W C 1996 Aquifer Test Analysis with Windows, CRC Press, Boca Raton, FL Walton, W C 1960 Leaky Artesian Aquifer Conditions in Illinois, Illinois State Water Survey Rept Invest 39, Urbana, 27p Wood, W W., T.P Kramer P.P Hern, Jr 1990 Intergranular diffusion: an important mechanism influencing solute transport in clastic aquifers Science 247 (March 30), 1569-1572 Xu, Moujin and Yoram Eckstein 1995 Use of weighted leastsquares method in evaluation of the relationship between dispersivity and field scale Ground Water 33, no 6, 905-908 Yeh, W., 1986 Review of parameter identification procedure in groundwater hydrology: the inverse problem, Water resources research 22(2), 95-108 250 Zheng, C., 1989 PATH3D, S.S Papadopulos &Assoc., Rockvill, MD Zheng, C., 1990 MT3D, a modular three-dimensional transport model, S.S Papadopulos & Assoc., Rockville, MD Zheng, C., G.D Bennet C.B Andrews, 1991 Analysis of ground water remedial alternatives at a Superfund site Ground Water 29 Zohdy, A.A., G.P Eaton, D.R Mabey, 1974 Application of surface Geophysics to Ground-water investigations, U.S Geological survey Technique of water resources investigations, Chap D1, Bk 2, 116p 251 Thông tin tác giả: Nguyễn Thu Hiền Năm sinh: 1964 Cơ quan công tác: Trường Đại học Thủy lợi Địa liên hệ: Khoa Kỹ thuật tài nguyên nước – Trường Đại học Thủy lợi Hồ Việt Hùng: Năm sinh: 1968 Cơ quan công tác: Trường Đại học Thủy lợi Địa liên hệ: Khoa Kỹ thuật tài nguyên nước – Trường Đại học Thủy lợi Trịnh Minh Thụ Năm sinh: 1964 Cơ quan công tác: Trường Đại học Thủy lợi Địa liên hệ: Khoa Cơng trình – Trường Đại học Thủy lợi Phạm vi đối tượng sử dụng: Cuốn sách sử dụng giáo trình tài liệu tham khảo cho sinh viên trường kỹ thuật thuộc ngành thủy lợi, xây dựng, giao thông, mỏ địa chất… bậc đại học cao học Từ khóa: nước ngầm, dòng chảy, tầng chứa nước, tầng cách nước, hệ số thấm, độ rỗng, truyền chất, giếng, trữ lượng Yêu cầu kiến thức trước học môn này: hồn thành mơn học đại cương thủy lực thủy văn Số lần xuất bản: 01 Nhà xuất bản: Nhà xuất giáo dục 252 ... động nước ngầm đặc biệt việc phát triển quản lý nguồn tài nguyên nước ngầm quan tâm năm gần Để đáp ứng nhu cầu thiết thực cho việc nghiên cứu lĩnh vực này, Giáo trình ? ?Phát triển quản lý tài nguyên. .. nguyên nước ngầm? ?? đề xuất khuôn khổ Tiểu hợp phần 1.3 “Hỗ trợ tăng cường lực cho Trường Đại học Thủy lợi? ?? thuộc dự án hỗ trợ ngành nước (WaterSPS) DANIDA để đưa vào chương trình đào tạo đại học. .. (dispersion) Quá trình phân tán thuỷ động học kết hai trình: trình phân tán học trình khuếch tán phân tử Quá trình phân tán học xảy nước ngầm bị ô nhiễm hoà trộn với nước ngầm không bị ô nhiễm, dẫn