ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC CÔNG NGHỆ BÙI HUY HỒNG NGHIÊN CỨU MỘT SỐ BÀI TỐN TƢƠNG TÁC NƢỚC MẶT – NƢỚC NGẦM VÀ THỬ NGHIỆM ÁP DỤNG VÀO THỰC TẾ Ở VIỆT NAM LUẬN VĂN THẠC SĨ CƠ HỌC KỸ THUẬT Hà Nội - 2012 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC CÔNG NGHỆ BÙI HUY HỒNG NGHIÊN CỨU MỘT SỐ BÀI TỐN TƢƠNG TÁC NƢỚC MẶT – NƢỚC NGẦM VÀ THỬ NGHIỆM ÁP DỤNG VÀO THỰC TẾ Ở VIỆT NAM Ngành: Cơ học kỹ thuật Chuyện ngành: Cơ học kỹ thuật Mã số: 60 52 02 LUẬN VĂN THẠC SĨ CƠ HỌC KỸ THUẬT HƢỚNG DẪN KHOA HỌC: PGS.TS HÀ NGỌC HIẾN Hà Nội - 2012 LỜI CAM ĐOAN Tôi xin cam đoan kết đạt luận văn nghiên cứu, sưu tầm, tổng hợp xếp lại phù hợp với yêu cầu luận văn Tồn điều trình bày khóa luận cá nhân, tham khảo tổng hợp từ nguồn tài liệu khác Tất tài liệu tham khảo, tổng hợp trích dẫn với nguồn gốc rõ ràng Tơi xin chịu hồn tồn trách nhiệm lời cam đoan Nếu có điều sai trái, tơi xin chịu hình thức kỷ luật theo qui định Hà Nội, ngày tháng năm 2012 Học viên Bùi Huy Hoàng MỤC LỤC DANH MỤC HÌNH DANH MỤC BẢNG MỞ ĐẦU CHƢƠNG TỔNG QUAN VỀ CÁC QUÁ TRÌNH LIÊN QUAN ĐẾN TƢƠNG TÁC NƢỚC MẶT NƢỚC NGẦM VÀ CÁC PHƢƠNG PHÁP TIẾP CẬN MƠ HÌNH HĨA 1.1 Quá trình gây dòng ngầm 1.2 Q trình gây dịng mặt 1.3 Q trình gây dịng mặt dịng ngầm 11 1.3.1 Tương tác nước mặt nước ngầm đất ngập nước 11 1.3.2 Bốc 12 1.4 Các mơ hình mô tƣơng tác nƣớc mặt nƣớc ngầm 12 1.4.1 Mơ hình nước ngầm tổng qt 12 1.4.2 Mơ hình nước mặt có tính đến nước ngầm 17 1.4.3 Mơ hình ghép nối nước mặt – nước ngầm 18 CHƢƠNG MƠ HÌNH MODFLOW 21 2.1 Mơ hình tốn học 21 2.2 Phƣơng pháp giải 22 2.3 Phƣơng trình vi phân 24 2.4 Điều kiện biên mơ hình 26 2.4.1 Biên Bốc 26 2.4.2 Biên sông 27 2.4.3 Biên tổng hợp 29 2.4.4 Biên kênh thoát (Drain) 30 2.4.5 Lỗ khoan hút ép nước 30 CHƢƠNG MỘT SỐ BÀI TOÁN MẪU VỀ TƢƠNG TÁC NƢỚC MẶT NƢỚC NGẦM 32 3.1 Tổng quan vể nghiệm giải tích 32 3.2 Kiểm định mơ hình MODFLOW với số tốn tƣơng tác nƣớc mặt-nƣớc ngầm 36 3.3 Sử dụng mơ hình MODFLOW đánh giá đặc trƣng tƣơng tác nƣớc mặt- nƣớc ngầm 38 3.3.1 Sự phụ thuộc vào đặc trưng thấm đáy sông 39 3.3.2 Sự phụ thuộc vào độ rộng sông 41 3.3.3 Sự phụ thuộc vào kích thước vỉa nước ngầm 41 3.4 Kết luận 42 CHƢƠNG ÁP DỤNG MƠ HÌNH MODFLOW TÍNH TOÁN TƢƠNG TÁC NƢỚC MẶT NƢỚC NGẦM THÀNH PHỐ ĐÀ NẴNG 44 4.1 Thiết lập mơ hình 44 4.1.1 Lưới sai phân mơ hình 44 4.1.2 Phân chia lớp mơ hình 45 4.1.3 Dữ liệu yếu tố địa hình địa chất địa chất thủy văn: 46 4.1.4 Dữ liệu khí tượng thủy văn 48 4.1.5 Biên điều kiện biên mơ hình 51 4.1.6 Hệ thống lỗ khoan quan trắc 53 4.2 Tính tốn hiệu chỉnh kiểm định mơ hình 54 4.3 Nghiên cứu tƣơng tác nƣớc mặt nƣớc ngầm vùng Đà Nẵng 56 4.3.1 Phương án 1: Phương án khai thác với lưu lượng khai thác giếng không đổi 58 4.3.2 Phương án 2: Phương án khai thác theo quy hoạch đến năm 2020 thành phố Đà Nẵng 65 4.3.3 Đánh giá trao đổi nước mặt nước ngầm theo hai phương án 70 KẾT LUẬN VÀ KIẾN NGHỊ 74 DANH MỤC CÁC BÀI BÁO 76 TÀI LIỆU THAM KHẢO 77 PHỤ LỤC 82 DANH MỤC HÌNH Hình 1.1 Những q trình liên quan đến tương tác nước ngầm - nước mặt Hình 2.1 Ơ lưới loại mơ hình 23 Hình 2.2 Ơ lưới nút lưới 24 Hình 2.3 Ơ lưới (i,j,k) bên cạnh 24 Hình 2.4 Điều kiện biên bốc cho mơ hình 27 Hình 2.5 Tính tốn lượng bốc 27 Hình 2.6 a) mặt cắt biểu diễn điều kiện biên sơng b) Mơ mơ hình 28 Hình 2.7 Biểu diễn đáy sơng 28 Hình 2.8 Biên tổng hợp 29 Hình 2.9 Điều kiện biên kênh thoát 30 Hình 2.10 Các lưới sai phân hai chiều xung quanh có lỗ khoan 31 Hình 3.1 Sơ đồ hệ thống mơ hình 34 Hình 3.2 Đặt toán Theis (1941) 35 Hình 3.3 Đặt tốn theo Hantush (1965) 35 Hình 3.4 Đặt toán Hunt (1999) 36 Hình 3.5 Đồ thị so sánh kết mơ hình với lời giải giải tích 38 Hình 3.6 Tốc độ suy giảm dòng chảy với trường hợp B với α=1 39 Hình 3.7 Tốc độ suy giảm dịng chảy ứng với trường hợp B α 40 Hình 3.8 Tốc độ suy giảm d ng chảy theo thời gian ứng với độ rộng sông 41 Hình 3.9 Tốc độ suy giảm d ng chảy theo thời gian ứng với kích thước vỉa nước ngầm 42 Hình 4.1 Bản đồ khu vực tính lưới sai phân mơ hình MODFLOW 44 Hình 4.2 Mặt cắt thẳng đứng theo hướng Bắc Nam mơ mơ hình 45 Hình 4.3 Mặt cắt thẳng đứng theo hướng Đơng Tây mơ mơ hình 45 Hình 4.4 Bản đồ đẳng bề mặt địa hình khu vực lập mơ hình 46 Hình 4.5 Bản đồ phân vùng hệ số thấm lớp 47 Hình 4.6 Bản đồ phân vùng hệ số nhả nước lớp 47 Hình 4.7 Biểu đồ phân bố lượng mưa số trạm 48 Hình 4.8 Biều đồ phân bố lượng bốc trung bình tháng Đà Nẵng 49 Hình 4.9 Sơ đồ phân bố lượng bổ cập 50 Hình 4.10 Sơ đồ phân bố lượng bốc 50 Hình 4.11 Sơ đồ phân bố điều kiện biên mực nước tầng chứa nước (qh,qp1), qp1, (-O1) 52 Hình 4.12 Vị trí lỗ khoan quan trắc vùng tính 53 Hình 4.13 So sánh cốt cao mực nước thực đo tính tốn giếng 54 Hình 4.14 Bản đồ mực nước tầng (qh,qp2) tốn chỉnh lý ổn định 55 Hình 4.15 Bản đồ mực nước tầng (qp1) toán chỉnh lý ổn định 55 Hình 4.16 Bản đồ mực nước tầng (ε-o1) toán chỉnh lý ổn định 56 Hình 4.17 Vị trí giếng khai thác khu vực tính tốn mơ mơ hình 57 Hình 4.18 Phân v ng tính tốn mơ hình 58 Hình 4.19 Lưu lượng trao đổi sơng tầng chứa nước v ng 59 Hình 4.20 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa v ng dự báo đến n m 2016 60 Hình 4.21 Lưu lượng trao đổi sông tầng chứa nước v ng 61 Hình 4.22 Biểu đồ lưu lượng trao đổi sông tầng chứa nước v ng dự báo đến n m 2016 62 Hình 4.23 Lưu lượng trao đổi sông tầng chứa nước v ng 63 Hình 4.24.Biểu đồ lưu lượng trao đổi sông tầng chứa nước v ng dự báo đến n m 2016 65 Hình 4.25 Biểu đồ lưu lượng nước trao đổi sông tầng chứa nước v ng dự báo đến n m 2016 67 Hình 4.26.Biểu đồ lưu lượng nước trao đổi sông tầng chứa nước v ng dự báo đến n m 2016 68 Hình 4.27 Biểu đồ lưu lượng nước trao đổi sông tầng chứa nước v ng dự báo đến n m 2016 69 Hình 4.28 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa theo phương án dự báo n m 2016 v ng 72 Hình 4.29 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa theo phương án dự báo n m 2016 v ng 73 Hình 4.30 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa theo phương án dự báo n m 2016 v ng 73 DANH MỤC BẢNG Bảng 1.1 Liệt kê mơ hình theo tác động đề cập 13 Bảng 3.1 Các tham số không thứ nguyên 38 Bảng 4.1 Lượng mưa trung bình tháng n m (mm) 48 Bảng 4.2 Phân bố lượng bốc trung bình tháng Đà Nẵng 49 Bảng 4.3 Vị trí lỗ khoan quan trắc mơ hình 53 Bảng 4.4 Quy hoạch khai thác nguồn nước TP Đà Nẵng đến n m 2020 57 Bảng 4.5 Lưu lượng nước trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 59 Bảng 4.6 Lưu lượng trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 61 Bảng 4.7 Lưu lượng trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 64 Bảng 4.8 Lưu lượng trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 66 Bảng 4.9 Lưu lượng trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 67 Bảng 4.10 Lưu lượng trung bình trao đổi sơng tầng chứa nước theo qu n m theo dự báo v ng (m3/ngày) 68 Bảng 4.11 Lưu lượng trao đổi sông tầng chứa nước phương án v ng tính n m 2016 (m3/ngày) 70 Bảng 4.12 Lưu lượng trao đổi sông tầng chứa nước phương án v ng tính n m 2016 (m3/ngày) 71 MỞ ĐẦU Tính cấp thiết đề tài Nước đất nguồn cung cấp nước quan trọng phục vụ đời sống sinh hoạt hoạt động sản xuất xã hội loài người Nước đất phận chu trình thủy v n hình thành nước mặt xâm nhập vào hệ thống đất đá từ bề mặt đất từ ao hồ sông suối mặt đất Các nguồn cung cấp cho nước đất gồm có: nước mưa, d ng chảy mặt, hồ ao, sông suối phương pháp bổ sung nước nhân tạo Thông thường việc quản l tài nguyên nước thường coi nước mặt nước ngầm hai thành phần tách biệt Một điều phải nước đất nước bề mặt tách rời cô lập chu trình thủy quyển, chúng chu trình kín Cả hai loại tài ngun ảnh hưởng đến số lượng chất lượng lẫn Sự trao đổi nước mặt nước đất xảy nơi sông, hồ, ao, đầm lầy… Sự ảnh hưởng qua lại nước mặt nước đất vấn đề quan trọng môi trường, nhiễm tầng nước mặt gây nhiễm bẩn tầng nước ngầm Trong nhiều trường hợp, nước mặt thu thêm nước chất hòa tan từ hệ thống nước ngầm trường hợp khác nước mặt nguồn cung cấp cho nước ngầm gây thay đổi chất lượng nước ngầm Kết là, mực nước sơng suối cạn gây cạn kiệt nguồn nước ngầm ngược lại, việc bơm hút nước ngầm gây cạn kiệt nguồn nước ao hồ sông suối vùng ngập nước Ơ nhiễm nước mặt gây suy thoái chất lượng nước đất ngược lại nhiễm nước ngầm làm suy thối chất lượng nước mặt Vì vậy, việc quản lý tài nguyên nước cách hiệu yêu cầu hiểu biết rõ ràng mối liên hệ nước đất nước mặt Mục tiêu nghiên cứu - Mục tiêu tổng quát Nghiên cứu chế tương tác nước mặt nước ngầm, ứng dụng mơ hình MODFLOW vào tính tốn mơ cho trường hợp thực tế, làm sở cho việc quy hoạch sử dụng nguồn nước ngầm cách bền vững - Mục tiêu cụ thể + Tìm hiểu trình tương tác nước mặt nước ngầm mơ hình tốn liên quan + Kiểm chứng mơ hình MODFLOW với toán tương tác nước mặt nước ngầm + Áp dụng mơ hình MODFLOW để tính tốn tương tác nước mặt nước ngầm khu vực Đà Nẵng + Trên sở kết mơ hình MODFLOW, đưa phân tích nhận xét tương tác nước mặt nước ngầm, tác động qua lại nước mặt nước ngầm khu vực Đà Nẵng, làm sở cho quy hoạch sử dụng nước Nội dung nghiên cứu cấu trúc luận văn Luận v n có cấu trúc gồm phần (Mở đầu, nội dung nghiên cứu, kết luận - kiến nghị) Phần nội dung nghiên cứu gồm chương sau: Chương Tổng quan trình liên quan đến tương tác nước mặt nước ngầm phương pháp tiếp cận mơ hình hóa: Phần trình bày trình gây nước mặt nước ngầm, tương tác nước mặt nước ngầm, mơ hình sử dụng nghiên cứu trước để tính tốn tương tác nước mặt nước ngầm Chương Mơ hình MODFLOW: Phần trình bày l thuyết mơ hình MODFLOW – mơ hình sử dụng luận v n để tính tốn đánh giá tương tác nước mặt nước ngầm Chương Một số toán mẫu tương tác nước mặt nước ngầm: Phần trình bày số tốn tương tác nước mặt nước ngầm, nghiệm giải tích mơ hình, so sánh kết tính tốn mơ hình giải tích Chương Áp dụng mơ hình MODFLOW tính tốn tương tác nước mặt nước ngầm cho thành phố Đà Nẵng 73 Hình 4.29 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa theo phương án dự báo năm 2016 v ng Hình 4.30 Biểu đồ so sánh lưu lượng trao đổi sông tầng chứa theo phương án dự báo năm 2016 v ng 74 KẾT LUẬN VÀ KIẾN NGHỊ - Kết luận Các kết Luận v n bao gồm: 1) Đã tìm hiểu trình quan trọng liên quan đến tương tác nước ngầm nước mặt cơng cụ mơ hình hóa chúng Các q trình phân loại sau: (1) Q trình gây dịng ngầm; (2) Q trình gây dịng mặt; (3) Q trình gây dịng ngầm dịng mặt Có nhiều phương pháp mơ hình số tính tốn q trình tương tác nước mặt nước ngầm, tính tốn tốc độ suy giảm d ng chảy sông với ưu nhược điểm khác tương ứng với phương pháp mơ hình MODFLOW mơ hình 2) Đã tìm hiểu mơ hình MODFLOW, áp dụng tính tốn so sánh với toán mẫu Hantush, Hunt, Glover Balmer tương tác nước mặt nước ngầm Kết tính tốn mơ hình với kết tính tốn theo giải tích theo lời giải Hantush, Hunt, Glover Balmer, cho ta kết gần Điều chứng tỏ mơ hình MODFLOW sử dụng tốt để nghiên cứu tốn tương tác nước mặt nước ngầm Mơ hình MODFLOW sau sử dụng để tính tốn với số đặc trưng tương tác, sở kết thu đươc, nhận thấy số tham số ảnh hưởng đến tương tác nước mặt với nước ngầm: hệ số rò rỉ dòng (stream leakage B) (liên quan đến độ thấm lớp đáy sông) hệ số ảnh hưởng nhiều đến tốc độ suy giảm d ng chảy sông hay lưu lượng trao đổi sông tầng chứa nước Việc thu thập khảo sát tính tốn tốn cụ thể cần phải hệ số Ngoài độ rộng sơng, kích thước tầng chứa nước ảnh hưởng nhiều đến kết tính tốn 3) Mơ hình MODFLOW sử dụng để tính tốn tương tác nước mặt - nước ngầm thành phố Đà Nẵng Mơ hình hiệu chỉnh với số liệu đo đạc sử dụng để đánh giá phương án khai thác nước ngầm Thành phố Sau mơ hình hiệu chỉnh với số liệu đo đạc, tiến hành tính tốn với kịch 4) Kết tính tốn phương án rằng: - Tại v ng 2, tương tác sông vả tầng chứa nước mạnh khu vực tập trung nhiều bãi giếng khai thác gần sơng, lịng sông sâu rộng, Mực nước ngầm tầng chứa nước khu vực bị giảm mạnh, lượng nước từ vỉa với nguồn nước vào khác không b đắp kịp cho lượng khai thác Điều có khả n ng t ng nguy nhiễm mặn t ng khai thác, mực nước sông vùng bị nhiễm mặn [3] 75 - Tại v ng 3, tương tác xảy yếu Qu 1, nước ngầm cung cấp cho sơng nhiều, có khả n ng t ng khai thác lên qu không 10000m3/ngày - Tại vùng 4, chủ yếu, nước ngầm cung cấp cho sơng Có thể t ng khai thác, song trữ lượng nước ngầm vùng ít, nên việc t ng khai thác khoảng 4000 m3/ngày - Lưu lượng trao đổi vỉa với sông vùng (khu vực cửa sông Hàn) PA giảm khoảng 5% so với PA1, làm giảm nguy ô nhiễm mặn khu vực khai thác nước ngầm - Kiến nghị Bài toán tương tác nước mặt – nước ngầm có nghĩa quan trọng vấn đề quản lý tổng thể nguồn nước lưu vực sông Để tiếp tục kết luận v n, số vấn đề nghiên cứu sau cần tiến hành: - Mở rộng lời giải giải tích với điều kiện thực tế khác Ví dụ trường hợp mực nước sơng thay đổi - Xây dựng mơ hình tương tác nước mặt nước ngầm có tính đến mơ hình dịng chảy lớp khơng bão hịa Về áp dụng với khu vực cụ thể, cụ thể luận v n thành phố Đà Nẵng, số liệu địa chất thủy v n khu vực khơng có nhiều thiếu tính đồng bộ, kết áp dụng cho khu vực thành phố Đà Nẵng luận v n sử dụng công tác quy hoạch nguồn nước với độ tin cậy chấp nhận Đối với toán thiết kế cần bổ sung số liệu quan trắc đồng chi tiết 76 DANH MỤC CÁC BÀI BÁO Hà Ngọc Hiến, B i Huy Hoàng ”Xây dựng mơ hình mơ động thái chất lượng nước ngầm cho thành phố Đà Nẵng”; Tuyển tập công trình – Hội nghị khoa học - Cơ học thủy khí tồn quốc n m 2010 Quy Nhơn 22-24/7, tr.229-238 77 TÀI LIỆU THAM KHẢO Tài liệu tiếng Việt Đỗ Cảnh Dương, Đỗ V n Bình Dự báo quy hoạch khai thác bền vững nguồn nước ngầm thành phố Đà Nẵng sở điều tra chất lượng, trữ lượng, trạng ô nhiễm khả tự bảo vệ nước đất Báo cáo khoa học, Đại học Mỏ Địa chất, 2003 Nguyễn Trường Đỉu, 1993 Báo cáo tìm kiếm đánh giá nước đất tỷ lệ 1:25000 v ng Đà Nẵng – Hội An, Quảng Nam TS Nguyễn Minh Sơn, Nghiên cứu đề uất giải pháp quản l s d ng tổng hợp tài nguyên nước lưu vực sông u ia-Sông Hàn đáp ứng nhu cầu phát triển bền vững thành phố Đà Nẵng Đề tài nghiên cứu khoa học công nghệ, Độc lập cấp nhà nước Quy hoạch phát triển nguồn nước thành phố Đà Nẵng đến n m 2020, UBND Tp Đà Nẵng, 2001 Tài liệu tiếng Anh Abu-Taleb M.F 2003 “Recharge of groundwater through multi-stage reservoirs in a desert basin” Environmental Geology 44: 379–390 Andersen H.E 2004 “Hydrology and nitrogen balance of a seasonally inundated Danish floodplain wetland” Hydrological Processes 18: 415–434 Bair E.S., Springer A.E and Roadcap G.S 1992 CAPZONE – an analytical flow model for simulating confined, leaky confined, and unconfined flow to wells with superposition of regional water levels International Ground Water Modeling Centers, Colorado School of Mines, Golden, Colorado 193 p Bakker M and Anderson E.I 2003 “Steady flow to a well near a system with a leaky bed.” Ground Water 41: 833–840 Bauer P., Held R.J., Zimmermann S., Linn F and Kinzelbach W 2006 “Coupled flow and salinity transport modelling in semi-arid environments: the Shashe River Valley, Botswana.” Journal of Hydrology 316: 163–183 10 Bochever, F.M., 1966, “Evaluation of well-field yield in alluvial aquifers: The impact of a partially penetrating stream”, Proceedings of VODGEO (Hydrogeology), No 13, 84-115, in Russian 11 Bouwer H 2000 “Integrated water management: emerging issues and challenges.” Agricultural Water Management 45: 217–228 78 12 Centre for Natural Resources 1999 IQQM: Integrated quantity and quality model Report No CNR99.025, 2nd edition NSW Department of Land and Water Conservation, Sydney 16 p 13 Chen X., Chen D and Chen X 2006 “Simulation of baseflow accounting for the effect of bank storage and its implications in baseflow separation.” Journal of Hydrology 327: 539–549 14 Conrad, L.P and M.S Beljin, 1996 “Evaluation of an induced infiltration model as applied to glacial aquifer systems.” Water Resour Bull., 32(6), 1209-1220 15 DHI Water and Environment 2005 MIKE SHE – integrated surface water and ground water 16 Dupuit J 1863 Étude théorique et pratique sur le mouvement des eaux dans les canaux découverts et a` travers les terrains perméables 2nd edition Dunod, Paris 17 Fitts, C.R., 1995 TWODAN manual, Scarborough, ME 18 Forchheimer P 1930 Hydraulik Teubner-Verlag, Leipzig 19 Glover R.E and Balmer G.G 1954 “River depletion resulting from pumping a well near a river.” American Geophysical Union Transactions 35: 368–470 20 Harbaugh A.W., Banta E.R., Hill M.C and McDonald M.G 2000 MODFLOW-2000, the USGeological Survey modular ground-water model – user guide to modularization conceptsand the ground-water flow process United States Geological Survey Open-File Report 00-92, Reston, Virginia 121 p 21 Hantush M 1965 “Wells near streams with semipervious beds.” Journal of Geophysical Research 70: 2829–2838 22 Hantush M 1967 “Depletion of flow in right angle stream bends by steady wells.” Water Resources Research 3: 235–240 23 Hantush M 2005 “Modelling stream aquifer interactions with linear response functions.” Journal of Hydrology 311: 59–79 24 Horton R.E 1933 “The role of infiltration in the hydrologic cycle.” American Geophysical Union Transcripts 14: 446–460 25 Hunt B 1999 “Unsteady stream depletion from groundwater pumping.” Ground Water 37: 98– 102 26 Hunt B 2003 “Unsteady stream depletion when pumping from semiconfined aquifer.” Journal of Hydrologic Engineering 8: 12–19 79 27 Jagelke J and Barthel R 2005 “Conceptualisation and implementation of a regional groundwater model for the Necker catchment in the framework of an integrated regional model.” Advances in Geosciences 5: 105–111 28 Jenkins A., Peters N.E and Rodhe A 1994 Hydrology In: Moldan B and Cerny J (editors) “Biochemistry of small catchments: A tool for environmental research.” John Wiley and Sons, Chichester p 31–54 29 Kondolf G.M., Maloney L.M and Williams J.G 1987 “Effects of bank storage and well pumping on base flow, Carmel river, Monterey County, California.” Journal of Hydrology 91: 351–369 30 Larkin, R.G., and Sharp, J.M., Jr., 1992, “On the relationship between river basin geomorphology, aquifer hydraulics, and ground-water flow direction in alluvial aquifers.” Geol Soc America Bull., v 104, p 1608-1620 31 Lin Y.C and Medina Jr M.A 2003 “Incorporating transient storage in conjunctive streamaquifer modeling.” Advances in Water Resources 26: 1001–1019 32 McDonald M.G and Harbaugh A.W 1988 A modular three-dimensional finite-difference ground-water flow model Book 6, Chapter A1 US Government Printing Office, Washington 33 Monninkhoff B 2002 “Coupling of the groundwater model FEFLOW with the hydrodynamic model MIKE11” German-Chinese Conference „Modern methods and instruments for water management and flood protection IWUTagungsberichte, p 161–173 34 Niswonger R.G et al 2006 “Modelling ground-water/surface-water interaction with GSFLOW, a new USGS model, MODFLOW and More; managing ground water systems.” IGWMC, Colorado School of Mines, p 99–103 35 Oxford Scientific Software 2004 A guide to Aquator Version 2.1 Oxford Scientific Software, Oxford 553 p 36 Prickett T.A., and Lonnquist C.G 1971 Selected Digital computer techniques for groundwater resource evaluation Illinois State Water Survey Bulletin 55 Champaign, Illinois 37 Raisin G., Bartley J and Croome R 1999 “Groundwater influence on the water balance and nutrient budget of a small natural wetland in Northeastern Victoria, Australia.” Ecological Engineering 12: 133–147 38 Rassam D., Pagendam D and Hunter H 2005b The Riparian Nitrogen Model: basic theory and conceptualisation CRC Catchment Hydrology and CRC Coastal Zones Technical Report 05/9 Canberra 80 39 Rassam D., Walker G and Knight J 2005a Applicability of the unit response equation to assess salinity impacts of irrigation development in the Mallee region: supplementary analyses Report for Water for a Healthy Country Flagship CSIRO, Townsville 40 Ruiz H.L and Rodriguez J.B.M 2002 “Explicit modeling of a reservoir and an aquifer.” Ingenieria Hidraulica en Mexico 17: 89–97 41 Shaw E.M 1988 Hydrology in practice 2nd edition Van Nostrand Reinhold, London 539 p.Shuttleworth W.J 1992 Evaporation In: Maidment D.R (editor) Handbook of Hydrology.McGraw Hill, New York p 4.1-4.53 42 Šimůnek J., Sejna M and van Genuchten M.T 1999 The HYDRUS-2D software package forsimulating two-dimensional movement of water, heat, and multiple solutes in variablysaturated media Version 2.0, IGWMC-TPS53 International Ground Water ModellingCentre, Colorado School of Mines, Golden, Colorado 43 Singh S.K 2000 “Rate and volume of stream depletion due to pumping.” Journal of Irrigation and Drainage Engineering-ASCE 126: 336–338 44 Strack O.D.L 1989 Groundwater Mechanics Prentice Hall, Englewood Cliffs, New Jersey 45 Swain E.D and Wexler E.J 1996 A coupled surface-water and groundwater flow model (MODBRANCH) for simulation of stream-aquifer interaction Techniques of Water Resources Investigations of the United States Geological Survey US Government Printing Office, Washington 46 Theis C.V 1941 “The effect of a well on the flow of a nearby stream.” American Geophysical Union Transactions 22: 734–738 47 Todd A.K., Buttle J.M and Taylor C.H 2006 “Hydrologic dynamics and linkages in a wetlanddominated basin.” Journal of Hydrology 319: 15–35 48 Victoria University and DSE 2005 REALM user manual REALM version 5.0, December 2005 262 p 49 Waddington J., Roulet N and Hill A 1993 “Runoff mechanisms in a forested groundwater discharge wetland” Journal of Hydrology 147: 37–60 50 Ward R.C and Robinson M 2000 Principles of hydrology 4th edition McGraw-Hill, London 450p 51 Winter T 1999 Relation of streams “Lakes, wetlands to groundwater flow systems.” Hydrogeology Journal 7: 28–45 81 52 Winter T., Harvey J., Franke O and Alley W 1998 Groundwater and surface water: a single resource USGS Circular 1139 US Geological Survey, Denver, Colorado 53 Van Elburg H., Engelen G.B and Hemker C.J 1993 FLOWNET manual Version 5.4 Free University of Amsterdam 82 PHỤ LỤC Hình Bản đồ đẳng đáy lớp Hình Bản đồ đẳng đáy lớp 83 Hình Bản đồ đẳng đáy lớp Hình Bản đồ đẳng đáy lớp 84 Hình Bản đồ đẳng đáy lớp Hình Bản đồ phân vùng hệ số thấm lớp 85 Hình Bản đồ phân vùng hệ số thấm lớp Hình Bản đồ phân vùng hệ số thấm lớp 86 Hình Bản đồ phân vùng hệ số thấm lớp Hình 10 Bản đồ phân vùng hệ số nhả nước lớp 87 Hình 11 Bản đồ phân vùng hệ số nhả nước lớp ... mơ hình sử dụng luận v n để tính tốn đánh giá tương tác nước mặt nước ngầm Chương Một số toán mẫu tương tác nước mặt nước ngầm: Phần trình bày số toán tương tác nước mặt nước ngầm, nghiệm giải... NGHỆ BÙI HUY HỒNG NGHIÊN CỨU MỘT SỐ BÀI TỐN TƢƠNG TÁC NƢỚC MẶT – NƢỚC NGẦM VÀ THỬ NGHIỆM ÁP DỤNG VÀO THỰC TẾ Ở VIỆT NAM Ngành: Cơ học kỹ thuật Chuyện ngành: Cơ học kỹ thuật Mã số: 60 52 02 LUẬN... dòng mặt dòng ngầm 1.3.1 Tương tác nước mặt nước ngầm đất ngập nước Winter et al (1998) [52] yếu tố tác động thuỷ v n đất ngập nước là: mưa, bốc hơi, tương tác nước mặt nước ngầm Lượng nước đầu vào