Nghiên cứu này trình bày phương pháp xác định độ lỗ rỗng hữu hiệu bằng phương pháp thí nghiệm ngoài thực địa ở khu vực gần ranh giới mặn nhạt NDĐ tầng qp1 tại xã Cẩm Hoàng, huyện Cẩm Giàng, Tỉnh Hải Dương, sau đó chỉnh lý tính toán bằng phương pháp mô hình số sử dụng mô hình Modflow và MT3D. Mời các bạn cùng tham khảo!
Bài báo khoa học Xác định đỗ lỗ rỗng hữu hiệu tầng chứa nước Pleistocen khu vực huyện Cẩm Giàng, tỉnh Hải Dương phương pháp thực nghiệm mô hình số Hồng Văn Duy1, Nguyễn Trung Hiếu2, Tống Thanh Tùng3, Đoàn Thu Hà4* Viện Khoa học Tài nguyên nước; hoangduydctv@gmail.com; Trường Đại học Thủy lợi; trunghieu.ma@hotmail.com; Liên đoàn quy hoạch điều tra tài nguyên nước miền Bắc; thanhtungtv51@gmail.com; Trường Đại học Thủy lợi; thuha_ctn@tlu.edu.vn *Tác giả liên hệ: thuha_ctn@tlu.edu.vn; Tel.: +84–948172299 Ban Biên tập nhận bài: 12/5/2021; Ngày phản biện xong: 29/7/2021; Ngày đăng bài: 25/10/2021 Tóm tắt: Độ lỗ hổng hữu hiệu (n0) thông số địa chất thủy văn chun mơn quan trọng, có ý nghĩa xác định mức độ lan truyển nhanh hay chậm vật chất gây ô nhiễm, nhiễm mặn nước đất Khu vực Hải Dương với mức độ mặn nhạt nước đất đan xen phức tạp, lưu lượng khai thác nước đất tăng lên theo thời gian, nguyên nhân gây xâm nhập mặn, lựa chọn khu vực nghiên cứu Sử dụng phương pháp thí nghiệm bơm nước ép dung dịch muối gần ranh giới mặn–nhạt nước đất cách sử dụng mơ hình Modflow MT3D xác định độ lỗ hổng hữu hiệu tầng Pleistocen (qp1) khu vực xã Cẩm Hoàng, huyện Cẩm Giàng, tỉnh Hải Dương 0,23 Từ khố: Nhiễm mặn; Hút nước thí nghiệm; Dung dịch chất thị; Độ lỗ rỗng hữu hiệu; Mô hình số Đặt vấn đề Một thơng số địa chất thuỷ văn quan trọng cần xác định điều tra đánh giá tài nguyên nước đất độ lỗ hổng hữu hiệu (n0), phần lỗ hổng đất đá chứa nước có khả cho nước qua Nghiên cứu xác định độ lỗ hổng hữu hiệu có ý nghĩa quan trọng việc tính tốn dịch chuyển ranh giới mặn nhạt tầng chứa nước, tính tốn khả lan truyền vật chất môi trường nước đất nhằm đánh giá khả năng, tốc độ nhiễm bẩn Khu vực thành thành phố Hải Dương, nơi nước đất có nguy xâm nhập mặn nhiễm bẩn cao, việc nghiên cứu xác định độ lỗ hổng hữu hiệu việc làm thiếu điều tra địa chất thuỷ văn nhằm đánh giá tài nguyên nước đất Xác định độ lỗ hổng hữu hiệu thực thí nghiệm thực địa lý tính tốn văn phịng phương pháp khác Vùng Hải Dương thưc thí nghiệm cụm lỗ khoan thí nghiệm CHD3 gần ranh giới mặn nhạt tầng chứa nước pleistocen (qp) xã Cẩm Hồng, huyện Cẩm Giàng (Hình 1) Chất thị muối ăn (NaCl) ép vào lỗ khoan, tiến hành bơm nước lỗ khoan khác Việc chỉnh lý xác định độ lỗ hổng hữu hiệu thực phương pháp mơ hình số Nước đất (NDĐ) khu vực Hải Dương khai thác phục vụ nhu cầu kinh tế xã hội khu vực từ xa xưa, nhiên mức độ mặn nhạt đan xen phức tạp ảnh hưởng việc khai thác ngày tăng làm tăng nguy xâm nhặp mặn Cùng với thời gian Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 http://tapchikttv.vn/ Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 23 khai thác, diện tích phễu hạ thấp mực nước toàn vùng ngày mở rộng tiến gần đến ranh giới mặn nhạt tầng chứa nước Pleistocen khu vực xã Cẩm Hoàng, huyện Cẩm Giàng, Tỉnh Hải Dương [1–6] Chất thị muối ăn (NaCl) ép vào cụm lỗ khoan thí nghiệm CHD3 nhằm xác định độ lỗ rỗng hữu hiệu n0, thông số ảnh hưởng đến lan truyền vật chất nước đất cần thiết phải xây dựng mơ hình dự báo q trình dịch chuyển ranh giới mặn nhạt nước đất phía trung tâm huyện Cẩm Giàng theo hướng dịng chảy Sơ đồ ranh giới mặn nhạt tầng chứa nước Pleistocen (qp1) thể Hình Nghiên cứu trình bày phương pháp xác định độ lỗ rỗng hữu hiệu phương pháp thí nghiệm ngồi thực địa khu vực gần ranh giới mặn nhạt NDĐ tầng qp1 xã Cẩm Hoàng, huyện Cẩm Giàng, Tỉnh Hải Dương, sau chỉnh lý tính tốn phương pháp mơ hình số sử dụng mơ hình Modflow MT3D Hình Sơ đồ vị trí cụm thí nghiệm ngồi thực địa Phương pháp nghiên cứu số liệu thu thập 2.1 Giới thiệu khu vực nghiên cứu Với đặc điểm địa chất thủy văn vùng Hải Dương, kết hợp với nhu cầu khai thác ngày tăng, diễn biễn chất lượng nước ranh giới mặn nhạt phức tạp Nghiên cứu xác định độ lỗ hổng hữu hiệu (n0) thực sử dụng hệ giếng có sẵn thuộc thơn Ngọc Lâu, xã Cẩm Hồng, huyện Cẩm Giàng, tỉnh Hải Dương Sơ đồ vị trí cụm thí nghiệm ngồi thực địa thể Hình Độ lỗ hổng hữu hiệu (n0) xác định cho tầng chứa nước khu vực gần ranh giới mặn nhạt NDĐ tầng qp1 Cấu trúc địa chất thủy văn địa điểm nghiên cứu thể Hình [1–11] Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 24 Hình Sơ đồ cấu trúc địa chất thủy văn khu vực Hải Dương 2.2 Cơ sở khoa học Phương trình lan truyền vật chất mơ tả chế dịch chuyển (đối lưu) phân tán thủy động lực NDĐ có hướng dịng chảy theo phương x không gian chiều (x, y) viết [7–8]: Dx 2C C C Ux R x x t (1) Trong Dx hệ số phân tán thủy động lực học theo hướng x (L2/T); C nồng độ chất hoà tan nước (M/L3); Ux (U = V/nhh) vận tốc thực dòng nước tương ứng theo hướng x y (M/T); V vận tốc Đắc-xi; nhh độ lỗ rỗng hữu hiệu; R hệ số chậm trễ; t thời gian (T) Hệ số phân tán thuỷ động lực xác định theo công thức sau [7–15]: Dx=D’x +D*d ; D’x=aLU (2) * Trong D’x hệ số phân tán học tương ứng theo hướng x y (L /T); D d hệ số khuếch tán phân tử môi trường rỗng (L2/T); aL độ phân tán dọc Phương trình (1) có lời giải có đầy đủ điều kiện ban đầu điều kiện biên mô tả sau: Điều kiện ban đầu phân bố nồng độ vật chất xem xét vào thời điểm ban đầu tùy ý t=t0 vị trí miền tính tốn: c co (x, y) (3) Các điều kiện biên đồng thời dạng sau: - Biên có nồng độ biết: C = Cc c (4) - Biên có gradient nồng độ pháp tuyến với đường biên biết (biên Neumann): c D xy Vn (cq c) q c qc (5) n - Biên có dịng vật chất khuếch tánlôi pháp tuyến với biên biết (biên Cauchi): Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 25 C V0Cv biên qc (6) n n Trong V0, Cv tương ứng dịng nước nồng độ chất hoà tan nước qua biên Vn C Dn 2.3 Phương pháp nghiên cứu 2.3.1 Phương pháp mơ hình hố xác định thơng số lan truyền Đối với trường hợp hút nước từ LK trung tâm ép dung dịch vào LK gần miền lan truyền có chế lan truyền [11] thể hình Hình Hai miền có chế lan truyền thí nghiệm hút nước ép chất thị Chiều rộng W đới thu nước từ phía thượng lưu dịng chảy vào LK ép nước đới cấp nước từ LK ép nước cho đới hạ lưu dịng chảy [20] có giá trị W 4rI hệ số thấm phần xung quanh LK nhỏ hệ số thấm tầng chứa nước Trong thí nghiệm cho nồng độ chất thị NDĐ LK hút nước, mơ hình chiều theo tuyến qua LK hút nước ép nước lại cho kết nồng độ chất thị rìa LK hút nước Vì sử dụng nồng độ tươngđối có dạng thể hình nồng độ chất thị quan trắc LK hút nước nồng độ chất thị xác định theo mơ hình chiều Như hai đường cong phân bố nồng độ tương đối LK hút nước điểm mép LK hút nước mơ hình chiều mặt lý thuyết trùng Nếu ký hiệu nồng độ chất thị nước LK hút nước CLKvới giá trị cực đại CLKmax cực tiểu CLKmin nồng độ xác định theo mơ hình chiều mép LK hút nước C1D với giá trị cực đại C1Dmax cực tiểu C1Dmin ta có nồng độ tương đối LK hút nước mép LK hút nước mơ hình chiều sau: C C LK C LK ; C LKmax C LK C C1D C1D C1D max C1D (7) CLKmin C1Dmax (a) Thời gian C=1 Nồng độ tương đối Nồng độ tuyệt đối CLKmax (b) C=0 Thời gian Hình Nồng độ tuyệt đối (a) tương đối (b) NDĐ LK hút nước mép LK hút nước mơ hình chiều Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 26 2.3.2 Sơ đồ thí nghiệm Sơ đồ thí nghiệm gồm ba lỗ khoan gồm lỗ khoan trung tâm CHD3 lỗ khoan hút nước hai lỗ khoan quan sát kết hợp ép dung dịch chất thị CHD3–1 CHD3–2 tạo thành tuyến cách 8m (Hình 5) Dung dịch thỉ sử dụng muối NaCl Hình Sơ đồ vị trí lỗ khoan thí nghiệm Lỗ khoan hút nước có đường kính 168 mm, lỗ khoan ép dung dịch chất thị có đường kính 90 mm Thời gian hút nước ép dung dịch thí nghiệm 168h Hút nước ép dung dịch muối bắt đầu lúc Lưu lượng hút nước 705 m3/ngày (8,15 l/s), lưu lượng ép dung dịch chất thị 43,2 m3/ngày (0,5 l/s) Nồng độ chất thị xác định qua công thức biểu diễn mối tương quan nồng độ muối ăn độ dẫn điện qui đổi điều kiện nhiệt độ 25oC Độ tổng khoáng hoá NDĐ tự nhiên 0,7 g/l, nồng độ muối ăn hoà vào NDĐ g/l tạo nên dung dịch có độ tổng khống hố 5,71 g/l dùng để ép vào lỗ khoan Sau hút nước thí nghiệm gần đạt trạng thái ổn định bắt đầu ép muối vào lỗ khoan quan sát CHD3–1 Nước lỗ khoan ép dung dịch xáo trộn ống ép dung dịch thường xuyên nâng lên hạ xuống toàn chiều dày cột nước lỗ khoan Đầu đo độ dẫn điện tự động đặt tầng chứa nước [1] Kết nghiên cứu thảo luận 3.1 Kết chỉnh lý mơ hình dịng chảy nước đất khu vực nghiên cứu 3.1.1 Nhập liệu lưới bước lưới mô hình Dựa vào bán kính ảnh hưởng lỗ khoan hút nước thí nghiệm chùm thí nghiệm CHD3 xác lập diện tích vùng lập mơ hình 16km2 Chia tồn vùng nghiên cứu thành lưới với số lượng 40 hàng 36 cột (Hình 6) Kích thước lưới 1x1m phù hợp để chạy tốn dịch chuyển vật chất Hình Thiết kế lưới sai phân diện tích vùng nghiên cứu Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 27 3.1.2 Nhập liệu phân tầng địa chất thủy văn Dựa mặt cắt cột địa tầng lỗ khoan chùm thí nghiêm xây dựng mơ hình cấu trúc địa chất thủy văn gồm lớp làm tảng để xây dựng mơ hình dịng chảy mơ hình dịch chuyển chất bẩn - Lớp 1: Lớp thấm nước yếu bề mặt (LCNBM) - Lớp 2: Tầng chứa nước lỗ hổng trầm tích Holocene (qh) - Lớp 3: Lớp cách nước thuộc hệ tầng Vĩnh Phúc aQ13vp2 (LCNVP) - Lớp 4: Tầng chứa nước lỗ hổng trầm tích Pleistocene (qp2) - Lớp 5: Lớp cách nước thuộc hệ tầng Hà Nội Q12–3hn2 (LCNHN) - Lớp 6: Tầng chứa nước lỗ hổng trầm tích Pleistocene (qp1) - Lớp 7: Lớp cách nước thuộc hệ tầng Vĩnh Bảo N2vb2 (LCNVB) 3.1.3 Nhập liệu thông số địa chất thủy văn Các thông số địa chất thủy văn tầng chứa nước chủ yếu dựa theo tài liệu hút nước thí nghiệm chùm CHD3 giai đoạn trước tính tốn theo điều kiện vận động khơng ổn định nhờ phần mềm Aquifer Test, có kết Bảng [1] Bảng Thông số địa chất thủy văn tầng chứa nước TT Số hiệu lỗ khoan Hệ số dẫn nước (Km, m /ngày) Hệ số thấm Hệ số nhả nước Ghi (K, m/ngày) đàn hồi () CHD3 1520 152 0,234 CHD3–1 1410 141 0,143 CHD3–2 1450 145 0,0114 3.1.4 Nhập liệu điều kiện biên mơ hình - Biên bổ cập (RECHARGE): Lượng bổ cập cho nước ngầm thời điểm bơm chùm thí nghiệm - Biên bốc (Evapo): Lượng bốc thời điểm bơm chùm thí nghiệm - Biên mực nước (CHD) gán cho tầng chứa nước qh, qp2, qp1 Trị số mực nước biên xác định theo tài liệu lỗ khoan quan sát CHD3–2 3.1.5 Nhập liệu trạng khai thác: Coi lỗ khoan hút nước cơng trình khai thác nước có cơng suất 705 m3/ngày 3.1.6 Nhập liệu mực nước ban đầu: Dữ liệu mực nước ban đầu lựa chọn thời điểm trước trước bơm chùm thí nghiệm 3.2 Chỉnh lý đánh giá sai số mơ hình 3.2.1 Kết chỉnh lý mơ hình dịng chảy nước đất Để chỉnh lý mơ hình theo tốn ổn định tiến hành dựa theo tài liệu quan trắc mực nước thực tế các lỗ khoan quan sát Mức độ tin cậy mơ hình xác định qua sai số mực nước mơ hình lỗ khoan quan sát thực tế Kết lý mơ hình cho thấy sai số lớn 0,155 m lỗ khoan CHD3–1, sai số nhỏ 0,011 m lỗ khoan CHD3–2, sai số trung bình 0,095 m sai số trung bình tuyệt đối 0,095 m (Hình 7) Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 28 Hình Đồ thị chỉnh lý mực nước 3.2.2 Kết chỉnh lý mơ hình dịch chuyển vật chất Dựa kết thí nghiệm đổ muối để tiến hành chỉnh lý thông số độ lỗ hổng hữu hiệu, hệ số phân tán, tỷ số khuếch tán dọc khuếch tán ngang Tại chùm thí nghiệm CHD3 tiến hành ép dung dịch muối ăn lỗ khoan CHD3–2, sở thiết lập thơng số cho mơ sau: - Các điều kiện biên nồng độ nguồn gây ô nhiễm Nồng độ không đổi lỗ khoan CHD3–1 với EC=2500 mg/l thể Hình Hình Gán biên nồng độ không đổi cho nguồn gây ô nhiễm Đối với điều kiện biên nồng độ ban đầu tầng chứa nước lựa chọn nồng độ trung bình EC0 = 1200 mg/l Điều kiện biên bổ cập mưa cho tầng chứa nước EC = mg/l - Các thông sên bổ cập mưa cho tầng chứa nước EC = mg/l Các thông số hệ số phân tán thủy động lực mơ hình lựa trọn sau: + Tỷ số hệ số phân tán ngang/ hệ số phân tán dọc lấy 0,1 (hay nói cách khác hệ số phân tán dọc αL = 10 αT) + Dd hệ số khuyếch tán phân tử cát, cuội, sỏi theo hướng dẫn chuyên môn 0,3 – 0,7x10–5 cm2/s, lựa chọn Dd = 0,0017m2/ngày Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 29 Bảng Các thông số dịch chuyển mơ hình Hệ số phân tán ngang (αL) 0,1 Hệ số phân tán dọc (αT) 0,01 Hệ số khuyếch tán phân tử (m2/ngày) 0,00011232 0,1 0,01 0,00011232 0,1 0,01 0,00011232 0,1 0,01 0,00011232 0,1 0,01 0,00011232 0,1 0,08 0,0017 0,1 0,01 0,00011232 Lớp Ghi (m2/ngày) - Khuyếch tán phân tử h.t, cuội, Chỉnh lý mơ hình dịch chuyển vật chất tiến hành dựa theo tài liệu quan trắc nồng độ EC các lỗ khoan quan sát Mức độ tin cậy mơ hình xác định qua sai số mực nước mơ hình lỗ khoan quan sát thực tế (Hình 9) Hình Đồ thị chỉnh lý nồng độ EC lỗ khoan CHD3 Theo kết nghiên cứu [16–20] độ lỗ rỗng tổng cộng cát cuội sỏi có độ chọn lọc tốt 0,25–0,5 cát lẫn cuội sởi 0,2–0,35 Đối với đất đá bở rời hạt thô độ lỗ rỗng hữu hiệu có giá trị tương đương với độ lỗ rỗng tồn phần khơng có tồn lỗ rỗng khơng thơng Vì giá trị độ lỗ rỗng hữu hiệu 0,23 xác định phương pháp trình bày viết theo kết thí nghiệm phù hợp, khơng có mâu thuẫn, nằm giới hạn giá trị thực cát cuội sỏi tầng chứa nước Pleistocen Kết luận Từ kết phân tích số liệu thí nghiệm trường mô phỏng, hiệu chỉnh phương pháp mơ hình số dịng chảy Modlfow dịch chuyển chất bẩn MT3D xác định thông số lan truyền độ lỗ rỗng hữu hiệu tầng chứa nước qp1 rút số kết luận sau: - Qua phân tích, sử dụng mơ hình số Modflow + MT3D kết thí nghiệm trường cơng cụ tốt để xác định độ lỗ rỗng hữu hiệu tầng chứa nước lỗ hổng phục vụ dự báo mức độ nhiễm mặn nước đất - Với việc thí nghiệm quan trắc nồng độ chất thị muối ăn (NaCl) LK hút nước phù hợp với thực tế khu vực nghiên cứu để đánh giá đỗ lỗ hữu hiệu tầng chứa nước; Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 30 - Phương pháp mơ hình số lan truyền chất hồ tan NDĐ xác định thông số theo phương pháp tổng hiệu bình phương hiệu tốn này; - Kết nghiên cứu khu vực Hải Dương xác định độ lỗ rỗng hữu hiệu tầng chứa nước qp1 0,23 phù hợp, khơng có mâu thuẫn, nằm giới hạn giá trị thực cát cuội sỏi tầng chứa nước Pleistocen Đóng góp cho nghiên cứu: Xây dựng ý tưởng nghiên cứu: H.V.D.; Lựa chọn phương pháp nghiên cứu: H.V.D., T.T.T.; Thu thập, phân tích, xử lý số liệu: H.V.D., N.T.H.; Viết thảo báo: H.V.D., N.T.H.; Chỉnh sửa báo: Đ.T.H Lời cảm ơn: Bài báo hoàn thành khuôn khổ thực Đề tài nghiên cứu theo Nghị định thư Bộ Khoa học Công nghệ Việt Nam Bộ Nghiên cứu giáo dục Cộng hòa liên bang Đức, “Nghiên cứu ứng dụng phát triển công nghệ khai thác nước thấm từ sông Việt Nam phục vụ sinh hoạt sản xuất”, mã số đề tài 60.GER–19 Lời cam đoan: Tập thể tác giả cam đoan báo cơng trình nghiên cứu tập thể tác giả, chưa công bố đâu, không chép từ nghiên cứu trước đây; khơng có tranh chấp lợi ích nhóm tác giả Tài liệu tham khảo 10 11 12 Tùng, T.T cs Báo cáo kết dự án: Bảo vệ nước đất đô thị lớn – Đô thị Hải Dương Lưu trữ Cục Quản lý tài nguyên nước, 2016 Quỳnh, C.V cs Báo cáo điều tra địa chất đô thị thành phố Hải Dương Lưu trữ Liên đoàn Quy hoạch Điều tra tài nguyên nước miền Bắc, Hà Nội, 1999 Quỳnh, C.V cs Báo cáo kết thăm dò tỷ mỷ nước đất vùng Cẩm Giàng, Hải Dương – Mỹ Văn, Hưng Yên, Lưu trữ Cục Quản lý tài nguyên nước, Hà Nội, 1997 Tuyển, C.T cs Báo cáo tìm kiếm nước đất vùng Văn Lâm – Văn Giang tỉ lệ 1:50.000, Lưu trữ Liên đoàn Quy hoạch Điều tra tài nguyên nước miền Bắc, Hà Nội, 1990 Bốn, C.V cs Báo cáo tìm kiếm nước đất vùng Nam Sách – Hải Hưng tỉ lệ 1:25.000, Lưu trữ Liên đoàn Quy hoạch Điều tra tài nguyên nước miền Bắc, Hà Nội, 1982 Quyền, P.B cs Báo cáo Điều tra, đánh giá tài nguyên nước vùng thủ Hà Nội, 2016 Hồng, N.V Mơ hình lan truyền chất ô nhiễm môi trường nước Nhà xuất Viện Hàn lâm Khoa học Công nghệ Việt Nam, Hà Nội, 2018 Hoàng, N.V Nghiên cứu xây dựng phần mềm mơ hình phần tử hữu hạn mơ chuyển động lan truyền chất ô nhiễm nhiễm mặn môi trường nước đất-ứng dụng cho khu vực ven biển miền Trung, 2017 Thanh, T.N cs Phương pháp luận xác định độ lỗ rỗng hữu hiệu độ phân tán tầng chứa nước áp dụng cho thí nghiệm trường phía Nam Hà Nội, 2018 Hạ, N.T Sự hình thành thành phần hóa học tầng chứa nước Pleistocen vùng Đồng Bằng Bắc Bộ, 2001 Hạ, N.T Sự hình thành thành phần hóa học nước đất trầm tích Đệ tứ vùng đồng Bắc ý nghĩa cung cấp nước, 2006 Duy, K Nghiên cứu tác động biến đổi khí hậu nước biển dâng đến xâm nhập mặn tầng chứa nước ven biển trầm tích Đệ Tứ đề xuất giải pháp thích ứng Tạp chí Khí tượng Thủy văn 2021, 730, 22-31; doi:10.36335/VNJHM.2021(730).22-31 13 14 15 16 17 18 19 20 31 vùng Hậu Lộc, tỉnh Thanh Hóa, 2015 Quy, P.N cs Lựa chọn khai thác tối ưu nhằm tránh xâm nhập mặn lỗ khoan khai thác nước đất tầng chứa nước ven biển, lấy vị dụ vùng tỉnh Ninh Thuận, 2020 Tường, P.V Đặc điểm thủy địa hóa thực trạng nhiễm mặn nước đất trầm tích Đệ Tứ vùng ven biển Hà Tĩnh, 2014 Bảng, P.V Nước đất vùng đồng Bắc 1984 Mantoglou, A Pumping management of coastal aquifers using analytical models of saltwater intrusion Water Resour Res 2003, 391(12), 1335 Brouyère, S Modeling tracer injection and well–aquifer interactions: a new mathematical and numerical approach Hydrogeology, Department of Georesources, Geotechnologies and Building Materials, University of Liège, Building B52, 4000 Sart Tilman, Belgium, 2008 Bear, J.; Verruijt, A Modeling groundwater flow and pollution D Reidel Publishing Company, Dordrecht, Holand, 1987 Fetter, C.W Applied Hydrogeology Prentice Hall Inc New Jersey 07458, 2001 Zlotnik, V.A.; Logan, J.D Boundary Conditions for Convergent Radial Tracer Tests and Effect of Well Bore Mixing Volume Water Resour Res 1996, 32(7), 2323– 2328 https://doi.org/10.1029/96wr01103 Ditermining the effective porosity of Pleistocene aquifer at Cam Hoang – Hai Duong by using experiment methods and numerical models Hoang Van Duy1, Nguyen Trung Hieu2, Tong Thanh Tung3 , Doan Thu Ha4* Institute of Water Resources Science; hoangduydctv@gmail.com Thuyloi University; trunghieu.ma@hotmail.com Divison of Northern Water Resources Planning and Investigation; thanhtungtv51@gmail.com Thuyloi University; thuha_ctn@tlu.edu.vn Abstract: The effective porosity (n0) is one of the important specialized hydrogeological parameters, which are significant in determining the fast or slow transmission of polluting and saline substances in groundwater In Hai Duong area, the groundwater salinity level is very complicated, the groundwater quantity extraction increases over time, which is the cause of saline intrusion Effective porosity of 0.23 was determined for saline–pale boundary area at Cam Hoang–Hai Duong Pleistocene aquifer (qp1) by using pumping and salt pressing experiment methods and the Modflow and MT3D models Keywords: Salinity; Laboratory water absorption; Indicator solution; Effective porosity; Numerical model ... tâm huyện Cẩm Giàng theo hướng dòng chảy Sơ đồ ranh giới mặn nhạt tầng chứa nước Pleistocen (qp1) thể Hình Nghiên cứu trình bày phương pháp xác định độ lỗ rỗng hữu hiệu phương pháp thí nghiệm thực. .. thông số lan truyền độ lỗ rỗng hữu hiệu tầng chứa nước qp1 rút số kết luận sau: - Qua phân tích, sử dụng mơ hình số Modflow + MT3D kết thí nghiệm trường công cụ tốt để xác định độ lỗ rỗng hữu hiệu. .. phương pháp tổng hiệu bình phương hiệu toán này; - Kết nghiên cứu khu vực Hải Dương xác định độ lỗ rỗng hữu hiệu tầng chứa nước qp1 0,23 phù hợp, khơng có mâu thuẫn, nằm giới hạn giá trị thực cát