34(2), 97-106 Tạp chí CÁC KHOA HỌC VỀ TRÁI ĐẤT 6-2012 ĐẶC ĐIỂM ĐỊA CHẤT TẦNG NÔNG DẢI VEN SÔNG HỒNG - TÂY HỒ - HÀ NỘI THEO TÀI LIỆU ĐỊA VẬT LÝ NGUYỄN VĂN GIẢNG1, NOBORU HIDA2, MAKSIM BANO3 E-mail: giangnv@igp-vast.vn Viện Vật lý Địa cầu - Viện KH&CN Việt Nam Akita University, Japan EOST - Strasbourg University, France Ngày nhận bài: - - 2012 Mở đầu Nghiên cứu cấu trúc địa chất - địa chất thủy văn dải ven sông Hồng thuộc thành phố Hà Nội tài liệu khảo sát địa vật lý góp phần làm sở khoa học cho việc triển khai đề án thấm lọc quy hoạch quản lý khai thác nguồn nuớc đất có ý nghĩa khoa học thực tiễn Như biết, vùng ven sông, hồ có điều kiện thuận lợi để xây dựng công trình khai thác thấm Nguyên lý làm việc công trình khai thác thấm ven bờ khai thác, mực nước đất hạ thấp xuống mực nước sông hồ, từ nguồn nước mặt thấm xuyên qua lớp đất đá bổ sung cho công trình khai thác [3] Thực chất dạng bổ sung tự nhiên đơn giản, thuận lợi mặt lợi dụng nguồn nước mặt sẵn có mà không cần phải xây dựng bồn chứa nhân tạo cho nước thấm, mặt khác công trình khai thác thường có lưu lượng lớn chất lượng tốt Sông Hồng, đoạn từ Việt Trì đến hết địa phận thành phố Hà Nội có nhiều điều kiện tự nhiên, điều kiện địa chất - địa chất thủy văn triển vọng tốt để xây dựng công trình khai thác thấm [4, 13] Để góp phần minh chứng cho luận điểm đây, tổ hợp phương pháp địa vật lý thực đo sâu điện đối xứng (VES) đo sâu mặt cắt điện (EP) kiểu Schlumberger; đo điện từ tần số thấp (VLF); đo Georadar (GPR), nhằm tìm cấu trúc địa chất - địa chất thủy văn tầng nông có khả đáp ứng cho công tác bổ sung tầng chứa nước đất dải ven sông Hồng, Tây Hồ, Hà Nội mà kết nghiên cứu khảo sát trình bày Vì muốn tiếp cận tối đa, trực tiếp đến cấu trúc dải ven sông dải cát lòng sông, nên phương pháp đo địa vật lý phải tiến hành thời kỳ mực nước sông hạ xuống thấp vào mùa khô Do đó, đợt khảo sát địa vật lý thực từ 2006 đến 2010 Sơ lược điều kiện tự nhiên vùng nghiên cứu Vùng nghiên cứu dải ven sông Hồng từ cầu Thăng Long đến cầu Long Biên thuộc địa phận quận Tây Hồ, thành phố Hà Nội Đây khu vực phẳng thuộc đồng aluvi trẻ tuổi Holocen hình thành từ 4000 năm trước [4, 13], độ cao tuyệt đối 6-8m có đặc điểm khí hậu chung với toàn đồng châu thổ sông Hồng nhiệt đới ẩm có gió mùa với nhiệt độ không khí bình quân năm 23,4°C Mùa nóng đồng thời mùa mưa kéo dài từ tháng đến tháng 10, chiếm khoảng 85% lượng mưa năm (lượng mưa trung bình năm 1532mm) nên thời kỳ có độ dư nước lớn Mùa khô từ tháng 11 đến tháng năm sau mùa mưa, lượng mưa lượng bốc hơi; thời kỳ hụt nước [7, 13] Sông Hồng sông lớn thứ hai Việt Nam bắt nguồn từ Trung Quốc chảy qua Việt Nam dài 510km Chiều rộng sông thay đổi theo thời kỳ năm từ 480 đến 1440m với lưu lượng bình quân thời kỳ từ 1990 đến 2640m3/s 97 Về mùa lũ mực nước sông dâng cao, thường 10m tức cao bề mặt đồng đê Mùa khô nước lại bị khai thác mạnh mẽ để tưới Nước sông Hồng nói chung tốt thành phần hoá học, riêng độ đục lớn đặc biệt mùa lũ Độ cao tuyệt đối đáy sông (chỗ sâu nhất) từ 0,2 đến 2,0m [4, 5] Khái quát đặc điểm địa chất Khu vực nghiên cứu nằm đới Hà Nội có lịch sử phát triển địa chất gồm giai đoạn lớn: Neoproterozoi - Cambri sớm, Mesozoi Neogen Đệ tứ Thời kỳ Neogen, tái hoạt động mạnh mẽ đứt gãy sâu Sông Lô, Sông Chảy, Vĩnh Ninh, vùng trũng Hà Nội hình thành cấu trúc dạng khối tảng Dọc đứt gãy sâu diễn trình tách giãn tạo nên địa hào, địa hình bị lún chìm, nước biển lấn sâu vào lục địa để lắng đọng hệ tầng Vĩnh Bảo (N2vb) dày khoảng 50-350m Đến cuối Pliocen lại nâng cao chịu bóc mòn Đến Đệ tứ trình lún chìm nâng cao diễn có tính chu kỳ dẫn đến biển tiến với trình lắng đọng trầm tích với kiểu nguồn gốc khác để hình thành hệ tầng Lệ Chi, Hà Nội, Vĩnh Phúc, Hải Hưng, Thái Bình [13] Thang địa tầng trầm tích Đệ tứ phân chia bảng Bảng Thang địa tầng Đệ tứ Giới Hệ Thống Holocen Kainozoi Đệ tứ Pleistocen Neogen Toàn diện tích khu vực nghiên cứu phủ trầm tích Đệ tứ phân làm phân vị với tuổi nguồn gốc khác gồm: hệ tầng Lệ Chi tuổi Pleistocen sớm (aQ11lc), hệ tầng Hà Nội tuổi Pleistocen - muộn (a, apQ12-3hn), hệ tầng Vĩnh Phúc tuổi Pleistocen muộn (a,lbQ13vp), hệ tầng Hải Hưng tuổi Holocen sớm - (lb,mQ21-2hh), hệ tầng Thái Bình tuổi Holocen muộn (a1; alb; a2Q23tb) Các trầm tích hệ tầng Hà Nội không lộ mặt mà bắt gặp nhờ lỗ khoan độ sâu khoảng 20m trở xuống, chiều dày đến 34m Mặt cắt hệ tầng phân làm lớp từ lên gồm cuội lẫn sỏi, sỏi lẫn cát hạt thô, bột sét lẫn cát hạt mịn Hệ tầng Thái Bình (Q23tb) trầm tích đại tạo thành sau biển lùi, mực nước biển hạ thấp, vai trò sông Hồng lớn dần trình hình thành đồng sông Hồng mà giai đoạn có trình xâm thực ngang chiếm ưu Trầm tích đại chủ yếu có nguồn gốc sông phân bố dọc hai bờ Hệ thống đê ven sông thiết lập dẫn tới phần trầm tích đê bị ngừng bồi đắp phù sa, đê hàng năm vào mùa lũ, bãi bồi lại phủ lớp mỏng phù sa cát bột sét màu mỡ Như vậy, lịch sử phát triển trầm tích Đệ tứ vùng nghiên cứu trải qua giai đoạn: 98 Phụ thống Ký hiệu Thượng Trung Hạ Thượng Trung Hạ Q2 Q2 Q2 Q1 Q1 Q1 N2 Niên đại tuyệt đối (năm ) 4000 6000 10.000 125.000 700.000 1.600.000 - Giai đoạn đầu: địa hình bị phân dị mạnh, dẫn tới trình xâm thực sâu, sông đào lòng mạnh đưa đến việc bồi lắng vật liệu hạt thô tạo nên trầm tích hệ tầng Lệ Chi, Hà Nội, chứa nước phong phú đối tượng cung cấp nước - Giai đoạn giữa: mức độ phân dị địa hình giảm trước, xâm thực sâu giảm Cuối giai đoạn trình xâm thực cân với trình tích tụ Nét đặc trưng trầm tích hệ tầng Vĩnh Phúc vào cuối giai đoạn bị phong hoá mạnh dẫn tới hình thành bề mặt dạng sét loang lổ - Giai đoạn cuối: nhiều biến cố địa chất xảy việc mở đầu tích tụ hồ đầm lầy ven biển với tích tụ than bùn Quá trình phân dị địa hình yếu, trình bồi tụ tăng cường Môi trường địa chất ngày chịu tác động trình nội sinh, ngoại sinh người Số liệu khảo sát địa vật lý thực địa Các tuyến đo sâu điện chạy dọc theo đường bờ sông cắt ngang qua sông Hồng xác định nguồn tài liệu phục vụ cho nghiên cứu cấu trúc địa chất thủy văn Trong thời gian thực phép đo điện điện từ vào mùa khô từ 2006 đến 2010, mực nước sông Hồng có lúc xuống thấp kỷ lục nên điều kiện thuận lợi để hoàn thành số tuyến đo Georadar cắt ngang qua sông dọc theo bờ phải sông Hồng Tất điểm đo sâu điện theo kiểu đo đối xứng Schlumberger (VES) với khoảng cách cực đại hai cực phát dòng AB=1000m tương ứng với hai cực thu 100m, đồng thời kiểu đo mặt cắt (EP) với khoảng cách hai cực phát cố định 400m sử dụng để thu thập số liệu địa điện thiết bị SAS 4000 hãng ABEM chế tạo [1, 9, 16, 17] Mật độ điểm đo tuyến VES 200m/điểm, EP 50m/điểm trừ đoạn tuyến cắt ngang qua dòng chảy lòng sông Hình trình bày sơ đồ vị trí tuyến đo VES EP dải ven sông Hồng Số lượng điểm đo cụ thể cho tuyến tổng hợp bảng sau đây: Bảng Tổng hợp điểm đo VES EP tuyến đo Tên tuyến đo Số điểm VES Số điểm EP Ghi SH1 16 68 Bỏ 600m khúc lòng sông SH2 25 100 Tuyến dọc bờ phải sông Hồng Hai tuyến đo Georadar thực khu vực nghiên cứu, tuyến thứ GPR1 cắt ngang qua sông dài 750m tuyến thứ hai GPR2 dọc theo bờ phải sông Hồng dài 1000m (hình 1, 2) Công việc đo Georadar tiến hành vào thời điểm nước sông Hồng xuống mức cạt kiệt môi trường mặt đất khô hoàn toàn Mục đích đo radar nhằm tìm cấu trúc địa chất phần nông thiết bị Pulse Ekko 100A với dải tần số anten thu-phát 100 50MHz [2, 6, 8, 12] Cũng có vị trí với tuyến GPR2 2000m tuyến đo điện từ tần số thấp (VLF) với bước đo 5m (hình 2) nhằm tìm đới cấu trúc có khả ngậm nước dập vỡ thiết bị VLF-Wadi [15] Trên sở điểm đo chập VES [12] 200m tuyến đo chập Georadar VLF [15] cho phép xác định sai số phép đo tài liệu VES 8%, tài liệu VLF 5% radar 2% Như vậy, đánh giá tài liệu đo đạt chất lượng tốt đủ điều kiện làm nguồn tài liệu đầu vào cho minh giải cấu trúc địa chất phần Hình Sơ đồ vị trí tuyến đo địa vật lý dải ven sông Hồng - Tây Hồ, Hà Nội 99 Hình Tuyến đo địa vật lý dọc bờ phải sông Hồng Kết minh giải tài liệu địa vật lý Mỗi loại đất đá có độ dẫn điện khác người ta đo giá trị điện trở suất đất đá ρ[Ωm] để xác định độ dẫn điện chúng Điện trở suất phụ thuộc vào thành phần thạch học, độ rỗng, độ bão hòa chất lưu đất đá, nhiệt độ áp suất nên loại đất đá khoáng vật chúng nằm dải giá trị định [14, 16, 17] Để có sở cho việc minh giải địa chất tài liệu địa điện, tác giả thực việc nghiên cứu xác định giá trị điện trở suất số mẫu vật cấu trúc địa chất khu vực nghiên cứu nhận kết trình bày bảng đây: Việc phân tích đường cong đo sâu điện dựa sở xây dựng mô hình lý thuyết để tiếp cận đến đường cong quan sát thực tế nhằm tìm lớp cấu trúc có bề dày giá trị điện trở suất biểu kiến tương ứng đặc trưng cho môi trường nghiên cứu Các tác giả sử dụng phương pháp phân tích lựa chọn [10] để thay đổi tham số mô hình tiên nghiệm cho đường cong lý thuyết mô hình trùng với đường cong đo đạc với sai số chấp nhận Hay nói khác cực tiểu hóa phiếm hàm sau: n {[ ] } G ( p) = ∑ ρ ei (ri , p) − ρ k (ri ) / ρ k (ri ) i =1 Trong đó: Bảng Điện trở suất số loại vật liệu Tây Hồ - Hà Nội Vật liệu Nước mặt Nước giếng khoan Đất trồng khô Đất trồng ướt Đất sét khô Đất sét ướt Cát khô Cát ướt 100 Điện trở suất [Ωm] 75 120 200 18 500 10 1000 50 p = {ρ1 , ρ , , ρ n , h1 , h2 , , hn −1 } tham số điện trở suất bề dày lớp mô hình ρei(ri,p) đường cong lý thuyết mô hình p ρk(ri) đường cong đo đạc thực tế Quá trình cực tiểu hóa phiếm hàm G ( p ) thực máy tính trình lặp, cho đường cong lý thuyết sát với đường cong đo đạc thực tế Do vậy, toán nghiệm Muốn có nghiệm ổn định phải tìm mô hình ban đầu sát với thực tế [10] phiếm hàm G ( p ) thay phiếm hàm điều chỉnh: Mα ( p ) = G ( p ) − αΩ( p ) toán Ω( p ) toán tử ổn định nghiệm Chúng ta thay đổi tham số lát cắt, cho đường cong mô hình lý thuyết trùng với đường cong đo đạc thực trực tiếp hình máy tính [5] Hình mặt cắt địa điện theo tài liệu tuyến SH1 cắt ngang qua sông Hồng với độ dài 4200m, cắt qua hai khúc lòng sông có hướng từ bờ hữu sang bờ tả Nhìn vào hình thái lớp phân chia phần nông mặt cắt cho thấy có đan xen số ranh giới Điều chứng tỏ cấu tạo đáy sông tầng bên bờ phải có liên kết với nước sông có thông thủy với tầng chứa nước ngầm bên bờ phải Hình trình bày mặt cắt địa điện tuyến SH2 dọc theo ven sông phía bờ phải theo hướng đông - tây với độ dài 5200m Các lớp cấu trúc phân chia theo mức độ dẫn điện chúng giá trị điện trở suất tính Ohm.m Hình thái lớp không phẳng, đoạn cuối tuyến (phía Tây), không biểu đường ranh giới liên tục mà bị cắt phần cấu trúc nông (đến 30m) Đặc điểm sở để liên tưởng đến kiểu cửa sổ địa chất thủy văn cấu trúc địa chất tầng nông Hình Mặt cắt địa điện theo tài liệu VES tuyến SH1 Hình Mặt cắt địa điện theo tài liệu VES tuyến SH2 101 Đối với số liệu Georadar mục tiêu nghiên cứu khảo sát xây dựng mặt cắt phản ánh trạng cấu trúc vật chất bên mặt đất với độ phân giải cao Dựa sở lý thuyết sóng radar cho trường hợp sử dụng kiểu đo phản xạ góc rộng cho thấy sóng phản xạ thu giản đồ tạo từ mặt ranh giới trung gian môi trường địa chất nên sóng phản xạ thường liên quan đến điều kiện thành tạo tự nhiên cấu trúc địa chất như: ranh giới lớp vật liệu trầm tích có tính chất vật lý khác nhau, khe nứt nẻ, khối xâm thực [2, 6, 8] Nghĩa là, phát sóng xuống đất dải tần số 50MHz sóng truyền sâu xuống môi trường, truyền mà phía trước có gặp mặt ranh giới (mặt trung gian) hai loại vật chất có tính chất điện khác môi trường biến đổi hệ số phản xạ R hệ số truyền qua T xác định Tính chất điện (độ dẫn, số điện môi) môi trường khảo sát có ảnh hưởng trực tiếp đến độ sâu nghiên cứu Khi độ dẫn môi trường khảo sát lớn 10mS/m (điện trở suất nhỏ 100Ohm.m) GPR bị hạn chế độ xuyên sâu Lúc muốn sâu tần số anten định Trong vật chất môi trường khảo sát nước thành phần quan trọng định tính chất điện môi trường [3, 8] Đối với môi trường địa chất ven sông Hồng, thấy có hai loại vật liệu đặc trưng, môi trường cát có độ dẫn khoảng 0,01mS/m môi trường sét, bùn sét có độ dẫn khoảng >10mS/m Trong xử lý tín hiệu radar, phải chọn tín hiệu rõ với tiêu chí hệ số phản xạ phải đạt giá trị tối thiểu 0,01 Do tuyến đo bố trí qua sông dọc theo ven sông với địa hình có nhiều vật cản nên chế độ thu thập số liệu chọn chế độ rời rạc Nghĩa là, người điều khiển trình thu thập số liệu cách ấn phím cách bàn phím nút bấm Trigger Odometer nối với khối điện tử Bước đo 0,5m cho loại anten có tần số 50MHz Như vậy, tổng số điểm đo tuyến lớn để bảo đảm điều kiện sử dụng chương trình Reflexw minh giải số liệu radar [11] Sau chuyển số liệu từ PulseEKKO sang dạng format Reflexw, ta bắt đầu bước xử lý chức lọc: 1D, 2D Trong đó, lọc 1D dùng để xử lý bước theo tia sóng (trace) riêng biệt; lọc 2D dùng để loại trace bị nhiễu minh giải số liệu [11, 12] Mặt cắt cấu trúc dạng cuối mà Reflexw đưa trình bày theo tuyến phản ánh trạng cấu trúc địa chất đến độ sâu 30m Kết giải đoán số liệu rada theo tuyến GPR1 GPR2 chia mặt cắt cấu trúc lớp khác (hình 5, 6) có loại vật liệu tương ứng như: - Đất bồi, - Cát pha bột sét, - Cát pha sét, - Sét pha cát, - Sét Hình thái lớp cấu trúc đa dạng, không theo dạng phân lớp ngang mà chúng đan cắt nhau, rõ lớp thứ thứ Kết tương đồng với kết mặt cắt cấu trúc theo số liệu VES hình mô tả phần Hình Mặt cắt Georadar tuyến GPR1 dài 750m tờ phải đến bờ trái sông Hồng (N-B) Ký hiệu: = Đất bồi, = Cát pha bột sét, = Cát pha sét, = Sét pha cát, = Sét 102 Hình Mặt cắt Georadar tuyến GPR2 dài 1000m theo hướng đông - tây (Đ-T) thuộc bờ phải sông Hồng Ký hiệu: = Đất bồi, = Cát pha bột sét, = Cát pha sét, = Sét pha cát, = Sét Do độ sâu nghiên cứu số liệu radar môi trường địa chất khu vực bị hạn chế (30m), nên muốn tìm hiểu hình thái cấu trúc sâu hơn, tác giả sử dụng số liệu đo mặt cắt điện theo kiểu EP để xây dựng đường cong điện trở suất biểu kiến ρa cho tuyến SH1 (hình 7) SH2 (hình 8) Hình 100 500 1000 1500 2000 2500 3000 3500 4000 (m) 3 2 9 7 6 5 4 3 100 ← Hình Đường cong ρa theo tài liệu EP tuyến SH1 500 1000 1500 2000 → Hình Đường cong ρa theo tài liệu EP tuyến SH2 ρa(Ωm) ρa(Ωm) dáng đường cong phản ánh cấu trúc bất đẳng hướng theo hai tuyến đo Như vậy, cấu trúc 30m tạo mặt lồi lõm rõ ràng điều kiện thuận lợi cho việc chứa liên thông tầng chứa nước ngầm vùng nghiên cứu [14] 2500 500 1000 3000 1500 3500 2000 4000 2500 3000 3500 4000 4500 5000 3 2 9 8 7 6 100 100 5 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 (m) 103 Đồng thời tác giả sử dụng số liệu VLF để nghiên cứu cấu trúc địa chất nằm gần mặt đất, xác định vị trí đới cấu trúc bị phá hủy, cấu trúc bất đồng nhất, đặc biệt ứng dụng vào việc xác định đánh giá trạng tầng chứa nước ngầm theo tuyến dài 2000m dọc bờ phải sông Hồng, có vị trí với tuyến GPR2 (hình 2) Trong dải tần số thấp, thông thường từ 15 đến 30kHz, sóng VLF truyền từ đài phát gặp bề mặt phân cách (ví dụ không khí/đất chẳng hạn) xuất dòng thứ cấp đất (tương tự dòng cảm ứng) Dòng tạo trường thứ cấp ngược dấu với trường sơ cấp Thiết bị Wadi ghi lại tỷ số (tỷ số phần trăm) cường độ trường theo phương thẳng đứng theo phương nằm ngang bề mặt Trái Đất Vì trường sơ cấp phát từ đài phát theo phương nằm ngang, điều chứng tỏ bình thường tỷ số xấp xỉ = Nó xấp xỉ = không tồn vùng dẫn môi trường địa chất đồng Sự lệch khỏi trạng thái thông thường chứng tỏ có dị thường xuất Khi tiến hành khảo sát thiết bị Wadi ta đồng thời xác định giá trị biên độ giá trị lệch pha trường thứ cấp so với trường sơ cấp ban đầu Với VLF-Wadi, đại lượng xác định thông qua phần thực phần ảo tương ứng trường đo [15] Thiết bị Wadi cho phép ghi nhận giá trị phần thực lẫn phần ảo (hình 9a), song thông tin phản ánh thực trạng cấu trúc môi trường thông qua giá trị phần thực Việc xử lý, minh giải dị thường VLF tiến hành qua bước với trợ giúp chương trình tính toán, xây dựng mô hình cấu trúc theo độ dẫn điện môi trường khảo sát Trên sở giá trị phân bố điện trở suất ta giải đoán trạng cấu trúc môi trường địa chất Mặt cắt theo độ dẫn điện thành lập số liệu VLF cho tuyến ven sông Hồng trình bày hình 9b, ta thấy đới có giá trị trị độ dẫn tương phản xuất xen kẽ mặt cắt, điều lần khẳng định vùng nghiên cứu có cấu trúc tầng nông bất đẳng hướng bất đồng rõ rệt a b Hình Mặt cắt theo độ dẫn điện tuyến VLF tuyến ven sông Hồng a) Số liệu VLF xử lý độ sâu 10m; b) Mặt cắt theo độ dẫn điện 104 Kết luận Số liệu đo sâu điện, đo Georadar đo VLF cho phép ta xây dựng mặt cắt cấu trúc tầng nông vùng nghiên cứu với độ xác độ phân giải chấp nhận, giúp ta liên kết để làm rõ đặc điểm cấu trúc bất đẳng hướng bất đồng dải ven bờ phải sông Hồng đáy sông thuộc Tây Hồ - Hà Nội Đến độ sâu 30m, mặt cắt cấu trúc phân chia lớp với hình thái xen kẹp nhau, lớp 3, tương ứng với độ sâu đáy sông vật liệu có khả thấm nước tốt tiền đề cho bổ cập nước sông cho tầng chứa nước ngầm khai thác Hình thái cấu trúc lớp độ sâu lớn 30m có bề mặt đa dạng điều kiện thuận lợi cho việc chứa liên thông tầng chứa nước ngầm vùng nghiên cứu Việc hình thành đới cấu trúc phân chia theo độ dẫn điện sở cho việc phân chia theo cấu trúc địa chất thủy văn, cấu trúc độ sâu đến gần 100m Tình hình khai thác nước đất tập trung Hà Nội tăng lên mạnh mẽ theo thời gian, khoảng 700.000m3/ngày năm 2011 Như vậy, công tác điều tra nghiên cứu đánh giá tài nguyên nước nói chung nước đất nói riêng phải đẩy mạnh, việc đánh giá nguồn bổ cập làm sở cho giải pháp phát triển bền vững nguồn nước kết nghiên cứu có ý nghĩa thực tiễn Lời cảm ơn: Nghiên cứu tài trợ Quỹ Phát triển Khoa học Công nghệ Quốc gia (NAFOSTED); đề tài mã số 105.04-2011.05 TÀI LIỆU DẪN [1] B Allred, N Fausey, J Daniels, R Ehsani, 2003: Applications of near-surface geophysical methods to midwest USA agriculture Proceedings of the EAGE 65th Annual Conference & Exhibition, Stavanger, Norway, 2-5 June 2003 [2] J.L Davis, A.P Annan, 1989: Ground penetrating radar for high resolution mapping of soil and rock stratigraphy Geophysical Prospecting, 37(5), p.531-551 [3] P De Beukelaar, et al, 2004: Monitoring of water infiltration using GPR data Near surface 2004 - 10th European Meeting of Environmental and Engineering Geophysics, Utrecht, Netherlands A027 [4] Nguyễn Văn Đản, Tống Ngọc Thanh, 2000: Về khả xây dựng công trình khai thác thấm lọc ven sông Hồng cung cấp nước cho thành phố Hà Nội TC Địa chất, A/260; tr.43-49 [5] Nguyễn Văn Giảng, 1998: Kết bước đầu nghiên cứu biến động môi trường nước ngầm Hà Nội phương pháp địa điện TC Các Khoa học Trái Đất, T 20, 1, tr.21-26 [6] Nguyễn Văn Giảng, 2000: Khả áp dụng radar xuyên đất địa kỹ thuật môi trường Việt Nam Tạp chí Địa chất, A/257 (3-4), tr.23-32 [7] Lê Văn Hiển (chủ biên), 2000: Nước đất đồng Bắc Bộ Cục Địa chất Khoáng sản Việt Nam, Hà Nội [8] S.Mc Greary, J.F Daly, 1998: Highresolution imaging of Quaternary coastal stratigraphy using GPR Proceedings of GPR’98 pp.273-277 [9] Michael S.Z., 1994: The geoelectrical methods in Geophysical exploration.Chapter Elsevier publishers [10] Nguyễn Trọng Nga, 1997: Thăm dò điện phân giải cao (Dùng cho học viên cao học Địa vật lý) Đại học mỏ Địa chất - Hà Nội [11] K.J Sandmeier REFLEX version 4.2 Copyright 1998 [12] Sensors & Software, 1996: PulseEKKO100RUN, User’s Guide, Version 1.2 Technical Manual 25 [13] Vũ Nhật Thắng (chủ biên), 2003: Địa chất tài nguyên khoáng sản thành phố Hà Nội Cục Địa chất Khoáng sản Việt Nam, Hà Nội [14] D.V Tuyen, T Canh, A Weller, 2000: The Application of Electrical tomography for solving Hydro-Engineering Geology problems in Vietnam, Advances in Natural Science, vol 1, No.3, pp.101-108 [15] VLF WADI instrument introduction, 1995: ABEM Sweden manual [16] Zdanov M.X., 1986: Thăm dò điện Nhà xuất Nhedra, Matxcova (tiếng Nga) [17] Zonge Engineering & Reseach Organization, Inc., 1999: The application of Surface Electrical Geophysics to Ground Water Problems p.27 105 SUMMARY The characteristics of shallow geological structure for Red River side - Tayho - Hanoi area by geophysical data The study of shallow geological structure by geophysical data is interested for many scientists by advanced geophysical technology nowadays The exploiting groundwater is too much by time for domestic and social-economic development The complex of geophysical methods as well as electrical sounding and profiling and electromagnetic methods such as Georadar and VLF were carried out by longitudinal and cross-lines in the study area for near-surface structure investigation The measured geophysical data were processed and interpreted by modern geophysical tools The results are clearly shown anisotropy and inhomongeneous of shallow structure by morphology for the study area The nearsurface structure consists of layers upto 30 m of depth revealed by GPR data The depths of the second and third layers in different locations are alternately corresponding to the river-bed levels This is a good condition for linked assessment between aquifers and possibility for recharge by the Red River water to groundwater on right side of the river 106