Tìm hiểu và nghiên cứu áp dụng phương pháp thăm dò điện đa cực mới vào thực tế

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - NGUYỄN THỊ THOA TÌM HIỂU VÀ NGHIÊN CỨU ÁP DỤNG PHƢƠNG PHÁP THĂM DÒ ĐIỆN ĐA CỰC MỚI VÀO THỰC TẾ LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - NGUYỄN THỊ THOA TÌM HIỂU VÀ NGHIÊN CỨU ÁP DỤNG PHƢƠNG PHÁP THĂM DÒ ĐIỆN ĐA CỰC MỚI VÀO THỰC TẾ Chuyên ngành : Vật lý Địa cầu Mã số : 60440111 LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS.TS Vũ Đức Minh Hà Nội – 2016 Lời cảm ơn Luận văn đƣợc hoàn thành Bộ môn Vật lý Địa cầu, Khoa Vật lý, Trƣờng Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội dƣới hƣớng dẫn khoa học PGS.TS Vũ Đức Minh Học viên xin bày tỏ lòng biết ơn sâu sắc tới PGS.TS Vũ Đức Minh cho học viên đƣợc tham gia đề tài nghiên cứu khoa học PGS, tận tình hƣớng dẫn giúp đỡ học viên suốt q trình thực hồn thành luận văn Học viên xin đƣợc gửi lời cảm ơn tới Ths Đỗ Anh Chung (Viện Sinh thái Bảo vệ cơng trình) giúp đỡ học viên trình thực địa Trong thời gian qua học viên đƣợc Bộ môn Vật lý Địa cầu tạo điều kiện cho đƣợc làm quen sử dụng thiết bị, máy móc Bộ mơn thực tế thu thập số liệu Trong q trình thực hồn thành luận văn, học viên nhận đƣợc quan tâm giúp đỡ, động viên thầy cô giáo Bộ mơn, Khoa Vật lý, Phịng Sau đại học, bạn bè đồng nghiệp Học viên xin đƣợc chân thành cảm ơn Học viên Nguyễn Thị Thoa MỤC LỤC MỞ ĐẦU .1 CHƢƠNG 1: GIỚI THIỆU PHƢƠNG PHÁP THĂM DÒ ĐIỆN ĐA CỰC MỚI .3 1.1 Xây dựng hệ cực đo điện đa cực 1.2 Giới thiệu qui trình khảo sát, thu thập xử lý số liệu phƣơng pháp Thăm dò điện đa cực 1.2.1 File điều khiển .6 1.2.2 Điều kiện tiếp đất 10 1.2.3 Thu thập số liệu thực địa 10 1.3 Phần mềm xử lý số liệu .11 CHƢƠNG 2: CÁCH XÁC ĐỊNH ĐƢỜNG BÃO HÒA THƠNG QUA THUẬT TỐN NỘI SUY VÀ XÁC ĐỊNH ĐIỂM UỐN CỦA ĐƢỜNG CONG ĐIỆN TRỞ SUẤT 12 2.1 Tóm tắt sở lý thuyết số thuật toán nội suy .12 2.1.1 Phƣơng pháp nội suy Polyfit .12 2.1.2 Phƣơng pháp nội suy Lagrange 14 2.1.3 Phƣơng pháp nội suy Spline 15 2.2 Tính tốn thử nghiệm thuật toán 17 2.2.1 Sử dụng phƣơng pháp nội suy Polyfit .17 2.2.2 Sử dụng phƣơng pháp nội suy Lagrange 18 2.2.3 Sử dụng phƣơng pháp nội suy Spline 19 2.3 So sánh phƣơng pháp lựa chọn thuật toán 20 CHƢƠNG 3: KẾT QUẢ ÁP DỤNG THỬ NGHIỆM 21 3.1 Khu vực đối tƣợng nghiên cứu 21 3.2 Phƣơng pháp nội dung nghiên cứu 21 3.3 Kết xác định đƣờng bão hịa đập Đồng Mơ .22 3.3.1 Sơ đồ tuyến khảo sát 22 KẾT LUẬN VÀ KIẾN NGHỊ 35 Danh mục bảng biểu Bảng 1.1: Hệ cực đo đối xứng cải tiến Bảng 2.1: Giá trị đo đạc tuyến theo chiều sâu 17 Bảng 3.1: Kết đo đƣờng bão hòa theo ống Pizomet thời điểm khảo sát 23 Bảng 3.2: Đƣờng bão hịa thân đập Đồng Mơ 32 Hình 3.18: Sơ đồ vị trí đƣờng bão hịa 32 Bảng 3.3: Kết so sánh xác định đƣờng bão hòa pizomet phƣơng pháp Trendline 33 Bảng 3.4: Kết so sánh xác định đƣờng bão hòa pizomet phƣơng pháp nội suy Polyfit 34 Danh mục hình vẽ, đồ thị Hình 2.1: Kết nội suy Polyfit điểm uốn 18 Hình 2.2: Kết nội suy Lagrange 19 Hình 3.1: Sơ đồ tuyến khảo sát pizomet A đập Đồng Mơ .23 Hình 3.2: Kết ảnh điện 2D tuyến DMAT1 24 Hình 3.3: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT1 sử dụng Trendline 24 Hình 3.4: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT1 sử dụng Polyfit 25 Hình 3.5: Kết ảnh điện 2D tuyến DMAT2 25 Hình 3.6: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT2 sử dụng Trendline 26 Hình 3.7: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT2 sử dụng Polyfit 26 Hình 3.8: Kết ảnh điện 2D tuyến DMAT3 27 Hình 3.9: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT3 sử dụng Trendline 27 Hình 3.10: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT3 sử dụng Polyfit 28 Hình 3.11: Kết ảnh điện 2D tuyến DMAT4 .28 Hình 3.12: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT4 sử dụng Trendline 29 Hình 3.13: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT4 sử dụng Polyfit 29 Hình 3.14: Kết ảnh điện 2D tuyến DMAT5 .30 Hình 3.15: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT5 sử dụng Trendline 30 Hình 3.16: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT5 sử dụng Polyfit 31 Hình 3.17: Kết ảnh điện 2D tuyến ngang đập 31 Hình 3.18: Sơ đồ vị trí đƣờng bão hòa 32 MỞ ĐẦU Năm 2001, PGS.TS Vũ Đức Minh đề xuất thành công hệ phƣơng pháp mới, phƣơng pháp Thăm dị điện cải tiến (The Improved Electrical Sounding method - IES), phƣơng pháp đƣợc công bố nhiều báo trƣớc [2, 5, 6, 15, 14] Phƣơng pháp Thăm dò điện đa cực (The Multi-electrode Resistivity Sounding method – MRS) hay gọi phƣơng pháp ảnh điện đa cực (The Multielectrode Resistivity Imaging method - MRI) [12] có ƣu điểm nhƣ: đo liên tục thu đƣợc số liệu tuyến, kết xử lý biểu diễn cho mặt cắt điện trở suất hay độ phân cực PGS.TS Vũ Đức Minh đề xuất phƣơng pháp Thăm dò điện đa cực cải tiến (The Improved Multi-electrode Electrical Sounding method – IMES) [8, 9] vào năm 2010 với việc sử dụng hệ cực đo cải tiến nhƣng thiết lập file điều khiển trình đo nhƣ phƣơng pháp Thăm dò điện đa cực truyền thống áp dụng có hiệu định, có nhiều ƣu việt bật, nhiên hạn chế hoàn thiện khảo sát 1D Để khai thác tối ƣu thiết bị có muốn phƣơng pháp Thăm dò điện đa cực thực đạt hiệu cao việc nghiên cứu môi trƣờng, PGS.TS Vũ Đức Minh nghiên cứu phát triển, tích hợp ƣu việt phƣơng pháp MRS phƣơng pháp IES để tạo phƣơng pháp Thăm dò điện đa cực (The New Advanced Multi-electrode Electrical Sounding method - AMES) đƣợc trình bày tóm tắt luận văn Đa số đập hồ chứa nƣớc ta đập đất hầu hết hệ thống quan trắc đƣờng bão hịa thân đập Mặt khác, đến thời gian sử dụng lâu năm tác động biến đổi khí hậu nên cơng trình bị xuống cấp nghiêm trọng Đƣờng bão hịa thân đập vị trí bề mặt dịng thấm, ln tồn thay đổi theo mực nƣớc hồ, phụ thuộc vào lớp đất đƣợc đắp thân đập… Khi mực nƣớc thƣợng đập dâng cao dẫn đến việc dâng cao đƣờng bão hòa, tăng gradient thấm thân đập, tăng áp lực nƣớc giảm thể tích khối đất khơng bão hịa Việc đồng thời tăng áp lực nƣớc giảm thể tích khối đất khơng bão hịa, tăng gradient thấm thân đập dẫn đến suy giảm cƣờng độ kháng cắt đất Vì vậy, việc xác định đƣợc vị trí đƣờng bão hịa số thời điểm điển hình mực nƣớc hồ so sánh với thiết kế ban đầu để kiểm tra làm sở nhận định, phân tích tình hình thấm thân đập, nhƣ đánh giá ổn định, an toàn đập việc làm cần thiết Với lý nêu trên, học viên mạnh dạn thực đề tài “Tìm hiểu nghiên cứu áp dụng phƣơng pháp Thăm dò điện đa cực vào thực tế” với mục đích làm quen góp phần vào việc áp dụng phƣơng pháp Thăm dị điện đa cực để xác định vị trí đƣờng bão hịa thân đập đất thơng qua việc sử dụng phƣơng pháp nội suy đa thức Cấu trúc luận văn: Mở đầu Chƣơng 1: Giới thiệu phƣơng pháp Thăm dò điện đa cực Chƣơng 2: Cách xác định đƣờng bão hịa thơng qua thuật tốn nội suy xác định điểm uốn đƣờng cong điện trở suất Chƣơng 3: Kết áp dụng thử nghiệm Kết luận kiến nghị CHƢƠNG 1: GIỚI THIỆU PHƢƠNG PHÁP THĂM DÒ ĐIỆN ĐA CỰC MỚI 1.1 Xây dựng hệ cực đo điện đa cực [13] Nhƣ ta biết, với hệ cực đo cải tiến (1D) giữ đƣợc tất ƣu điểm, mà cịn khắc phục đƣợc nhƣợc điểm phƣơng pháp đo sâu điện vi phân trƣớc Các hệ cực đo hoàn toàn tƣơng tự nhƣ hệ cực đo truyền thống (chỉ khác kích thƣớc hệ số hệ cực nên đo gối đầu kích thƣớc, trừ kích thƣớc cuối cùng, đặc biệt đến khoảng cách hệ cực vị trí điện cực phát trùng với vị trí điện cực thu trƣớc đó), đồng thời sử dụng nguyên lý tƣơng hỗ thăm dị điện nên bố trí cặp cực phát bên làm thu ngắn đƣờng dây phát, đơn giản việc thi công thực địa giảm chi phí cho quy trình đo đạc, đồng thời chống rị điện bảo vệ an toàn lao động Mặt khác, hệ cực đo cải tiến có ƣu việt điểm khảo sát cần sử dụng loại hệ cực đo đó, giảm thiểu phép đo mà ta có đƣợc tất thơng tin loại đƣờng cong khác sau xử lý qua phép biến đổi đại số đơn giản, nhƣ lƣợng thông tin thu đƣợc tăng lên gấp bội Đối với hệ cực đo đa cực (2D) điện cực đƣợc cắm lúc tuyến nên dịch chuyển nhiều lần nhƣ hệ cực đo 1D; đo đồng thời với nhiều loại hệ cực đo khác nhau; khoảng cách hệ cực nên số liệu thu đƣợc độ sâu cách làm tăng độ chi tiết môi trƣờng khảo sát sâu so với hệ cực đo truyền thống tăng theo hệ số log Từ nhận xét trên, mục tiêu xây dựng đƣợc hệ cực đo đa cực (gọi tắt hệ cực đo đa cực MC) nhằm phát triển, tích hợp đƣợc ƣu việt hệ cực đo đa cực truyền thống (2D) hệ cực đo cải tiến 1D Trƣớc hết phải tiến hành chuyển đổi từ hệ cực đo cải tiến 1D với khoảng cách hệ cực không thành hệ cực đo có khoảng cách hệ cực 1D để từ ghép thành hệ cực đo đa cực MC (2D) Ta thấy với hệ cực đo cải tiến, với lần phát dòng AB ta thu đƣợc giá trị hiệu từ M1N1 M2N2, với quy luật chung với lần đo khác ln có M1N1/AB M2N2/AB giá trị không đổi (bảng 1.1) Bảng 1.1: Hệ cực đo đối xứng cải tiến AB/2 MN/2 MN/AB 0,5 2,8 5,6 0,5 4,0 0,7 4,0 5,7 0,7 5,6 1,0 5,6 5,6 1,0 8,0 Từ bảng 1.1 ta thấy: M1N1/AB = 5,6 M2N2/AB = Nhƣ vậy, với hệ cực đo cải tiến tăng AB hay MN MN/AB ln ln khơng đổi khoảng cách hệ cực đo tăng theo tỷ lệ log (khơng đều) Cịn hệ cực đo đa cực đo đồng thời với nhiều loại hệ cực khác file điều khiển trình đo mà điện cực thay đổi cách linh hoạt liên tục theo mục tiêu ngƣời sử dụng; đo liên tiếp điểm đo sâu với khoảng cách hệ cực đo tăng theo tỷ lệ số tự nhiên (khoảng cách hệ cực đo đều) Để tích hợp đƣợc ƣu điểm trên, PGS.TS Vũ Đức Minh đề xuất hệ cực đo đa cực MC có khoảng cách hệ cực đo đều, với lần đo AB M1N1/AB = a, M2N2/AB = b (b > a) Lần đo thứ cực AB tăng thay vào vị trí cực M1N1, M1N1 thay vị trí M2N2 trƣớc Với lần đo thứ n AB vị trí cực Mn-1Nn-1, M1N1 vị trí MnNn Trƣớc mắt, lựa chọn hai hệ cực đo thỏa mãn điều kiện nêu nhƣ sau: Hình 3.2: Kết ảnh điện 2D tuyến DMAT1 [7] TT 10 11 12 13 14 15 Chiều sâu (m) 0,00 -0,46 -1,43 -2,50 -3,68 -4,97 -6,39 -7,96 -9,68 -11,57 -13,66 -15,95 -18,47 -21,24 -22,69 ĐTS (Ωm) 346 319 270 498 653 598 469 354 293 264 251 249 257 280 292 ĐTS Trendline Hình 3.3: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT1 sử dụng Trendline[7] Kết xác định có điểm uốn chiều sâu 5,38m với giá trị điện trở suất 563Ωm 24 Hình 3.4: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT1 sử dụng Polyfit Kết xác định có điểm uốn chiều sâu 5,52m với giá trị điện trở suất 547Ωm  Kết khảo sát phƣơng pháp AMES tuyến DMAT2 đƣợc biểu diễn hình 3.5 Hình 3.5: Kết ảnh điện 2D tuyến DMAT2 [7] 25 TT 10 11 12 13 14 Chiều ĐTS sâu (m) (Ωm) 931 - 0,46 820 - 1,43 637 - 2,50 820 - 3,68 992 - 4,97 965 - 6,39 826 - 7,96 683 - 9,68 578 -11,57 512 -13,66 478 -15,95 465 -18,47 469 -21,24 483 ĐTS Trendline Hình 3.6: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT2 sử dụng Trendline[7] Kết xác định có điểm uốn chiều sâu 7,2m với giá trị điện trở suất 761Ωm Hình 3.7: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT2 sử dụng Polyfit Kết xác định có điểm uốn chiều sâu 6,02m với giá trị điện trở suất 858Ωm 26  Kết khảo sát phƣơng pháp AMES tuyến DMAT3 đƣợc biểu diễn hình 3.8  Hình 3.8: Kết ảnh điện 2D tuyến DMAT3 [7] TT 10 11 12 13 14 Chiều ĐTS sâu (m) (Ωm) 0,0 612 -0,5 599 -1,4 573 -2,5 391 -3,7 473 -5,0 597 -6,4 668 -8,0 577 -9,7 431 -11,6 322 -13,7 264 -15,9 245 -18,5 243 -21,2 258 ĐTS Trendline Hình 3.9: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT3 sử dụng Trendline[7] Kết xác định có điểm uốn chiều sâu 8,54m với giá trị điện trở suất 495Ωm 27 Hình 3.10: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT3 sử dụng Polyfit Kết xác định có điểm uốn chiều sâu 8,2m với giá trị điện trở suất 551Ωm  Kết khảo sát phƣơng pháp AMES tuyến DMAT4 đƣợc biểu diễn hình 3.11  Hình 3.11: Kết ảnh điện 2D tuyến DMAT4 [7] 28 TT 10 11 12 13 14 Chiều sâu (m) 0,0 -0,4 -1,4 -2,4 -3,5 -4,8 -6,1 -7,6 -9,3 -11,1 -13,1 -15,2 -17,6 -18,9 ĐTS (Ωm) 711 829 1127 1292 1000 675 478 379 336 325 333 352 376 389 ĐTS Trendline Hình 3.12: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT4 sử dụng Trendline [7] Kết xác định có điểm uốn chiều sâu 6,5m với giá trị điện trở suất 528 Ωm Hình 3.13: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT4 sử dụng Polyfit 29 Kết xác định có điểm uốn chiều sâu 7,9m với giá trị điện trở suất 375Ωm  Kết khảo sát phƣơng pháp AMES tuyến DMAT5 đƣợc biểu diễn hình 3.14 Hình 3.14: Kết ảnh điện 2D tuyến DMAT5 [7] TT 10 11 12 13 14 15 Chiều sâu (m) -0,46 -1,43 -2,50 -3,68 -4,97 -6,39 -7,96 - 9,68 - 11,57 - 13,66 -15,95 -18,47 - 21,24 -22,69 ĐTS (Ωm) 614 657 753 745 692 601 488 387 317 278 263 266 283 308 321 ĐTS Trendline Hình 3.15: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT5 sử dụng Trendline [7] Kết xác định điểm uốn chiều sâu 5,3m giá trị ĐTS 545Ωm 30 Hình 3.16: Đƣờng cong giá trị ĐTS vị trí 54,5m tuyến DMAT5 sử dụng Polyfit Kết xác định điểm uốn chiều sâu 5,38m giá trị ĐTS 567Ωm  Kết khảo sát phƣơng pháp AMES truyến ngang đập đƣợc biểu diễn hình 3.17 Hình 3.17: Kết ảnh điện 2D tuyến ngang đập [7] Qua kết xử lý phân tích nêu ta xác định đƣợc đƣờng bão hòa thân đập Đồng Mơ ba phƣơng pháp nhƣ trình bày bảng 3.2 31 Bảng 3.2: Đƣờng bão hòa thân đập Đồng Mơ Cao trình (m) TT Khoảng cách (m) Thân đập 10 10,5 19 20 26 28 36 38 49,5 51,5 18,95 22,5 25,7 26,2 26,4 25,9 22,9 22,1 18,67 18,6 ĐBH theo pizomet 18,95 ĐBH theo PP hàm dự báo gần 18,95 17,12 18,5 ĐBH theo PP nội suy 18,95 16,98 19,68 18 17,4 17,36 16,4 17,4 15 13,3 13,22 15,4 13,97 Hình 3.18: Sơ đồ vị trí đƣờng bão hịa Để đánh giá sai số đƣờng bão hòa thực (đo pizomet) đƣờng bão hòa đƣợc xác định phƣơng pháp Trendline, phải tiến hành hiệu chỉnh đƣờng bão hòa đƣợc xác định phƣơng pháp Trendline vị trí ống pizomet Muốn vậy, từ kết xác định đƣờng bão hòa phƣơng pháp hàm 32 dự báo gần đo trực tiếp qua ống pizomet (bảng 3.2) ta xây dựng đƣờng dự báo (Trendline - hình 3.18) đƣờng bão hịa thân đập Từ phƣơng trình đƣờng dự báo, cho giá trị vị trí ống pizomet ta xác định đƣợc chiều sâu cao trình đƣờng bão hịa vị trí ống pizomet so với mặt đập (bảng 3.3) Bảng 3.3: Kết so sánh xác định đƣờng bão hịa pizomet phƣơng pháp Trendline Cao trình (m) Chiều sâu (m) Sai lệch (m) Sai số (%) Thân đập ĐBH theo pizomet ĐBH theo PP hàm dự báo gần 26,2 18 18,25 8,2 0,25 26,4 17,4 17,60 9,0 0,2 22,1 15,4 16,00 6,7 0,6 18,67 13,97 13,70 4,7 0,27 TT Từ kết xác định đƣờng bão hòa phƣơng pháp ảnh điện so sánh với kết đo pizomet cho thấy sai số phƣơng pháp ảnh điện 2D đập Đồng Mô lớn khoảng 9% Để đánh giá sai số đƣờng bão hòa thực (đo pizomet) đƣờng bão hòa đƣợc xác định phƣơng pháp nội suy, phải tiến hành hiệu chỉnh đƣờng bão hòa đƣợc xác định phƣơng pháp nội suy vị trí ống pizomet Muốn vậy, từ kết xác định đƣờng bão hòa phƣơng pháp nội suy đo trực tiếp qua ống pizomet (bảng 3.2) ta xây dựng đƣờng nội suy đƣờng bão hòa thân đập Từ phƣơng trình đƣờng nội suy, cho giá trị vị trí ống pizomet ta xác định đƣợc chiều sâu cao trình đƣờng bão hịa vị trí ống pizomet so với mặt đập (bảng 3.4) 33 Bảng 3.4: Kết so sánh xác định đƣờng bão hòa pizomet phƣơng pháp nội suy Polyfit Cao trình (m) TT Chiều Sai sâu lệch (m) (m) Sai số (%) Thân đập ĐBH theo pizomet ĐBH theo PP nội suy 26,2 18 18,06 8,2 0,06 0.7 26,4 17,4 17,31 9,0 0,09 22,1 15,4 15,38 6,7 0,02 0.3 18,67 13,97 13,45 4,7 0,52 11 Từ kết xác định đƣờng bão hòa phƣơng pháp ảnh điện so sánh với kết đo pizomet cho thấy sai số phƣơng pháp nội suy đập Đồng Mô lớn khoảng 1%, cao trình 18,67 thân đập có sai số bất thƣờng (11%) 34 KẾT LUẬN VÀ KIẾN NGHỊ Qua trình thực đề tài luận văn, học viên học hỏi, tiếp thu đƣợc kiến thức bổ sung cho hiểu biết chuyên ngành theo học, đặc biệt kiến thức thực tế, cụ thể nhƣ sau: - Đã tìm hiểu nắm đƣợc nguyên tắc phƣơng pháp Thăm dò điện dòng chiều truyền thống, nhƣ phƣơng pháp Thăm dò điện cải tiến sau - Đƣợc tham gia đề tài nghiên cứu khoa học, làm quen với công việc phải thực từ việc chu n bị, tiến hành đo đạc thu thập số liệu thực địa, xử lý số liệu biểu diễn kết báo cáo Sau hoàn thành luận văn này, học viên thu nhận đƣợc số kết nhƣ sau: - Tìm hiểu thử nghiệm số phƣơng pháp nội suy đa thức nhằm lựa chọn phƣơng pháp phù hợp để áp dụng xử lý số liệu đo sâu điện thu đƣợc phƣơng pháp Thăm dò điện đa cực với mục đích xác định đƣờng bão hịa thân đập đất thơng qua việc sử dụng đạo hàm bậc tìm điểm uốn đƣờng cong điện trở suất - gianh giới thay đổi điện trở suất lớn - Áp dụng phƣơng pháp nội suy đa thức lựa chọn, lập chƣơng trình ngơn ngữ Matlab để tìm đƣờng bão hịa thân đập Đồng Mô - Kết thu đƣợc so sánh với kết hệ thống Pizomet cho thấy hồn tồn phù hợp có độ xác cao trƣớc Từ đó, học viên có số nhận xét kiến nghị: - Phƣơng pháp Thăm dò điện đa cực với việc áp dụng phƣơng pháp nội suy đa thức hồn tồn áp dụng để tìm đƣợc đƣờng bão hịa thân đập đất tốt phƣơng pháp truyền thống trƣớc 35 - Để nâng cao độ xác kết thu đƣợc, cần phải tiếp tục nghiên cứu cải tiến khâu xử lý số liệu mà trƣớc mắt áp dụng xử lý theo phƣơng pháp cải tiến đƣợc đề xuất với thông số mới, đồng thời nên tiến hành đo xử lý 3D mà khuôn khổ luận văn chƣa thể tiến hành đƣợc 36 TÀI LIỆU THAM KHẢO Tiếng Việt Đỗ Anh Chung, Vũ Đức Minh (2013), “Nghiên cứu cải thiện khả tiếp đất điện cực phƣơng pháp điện đa cực cho mơi trƣờng khó tiếp đất”, Tạp chí khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ, 29(2), 57-69 Lê Viết Dƣ Khƣơng, Vũ Đức Minh (2001), “Các phƣơng pháp đo sâu điện trở dùng tổ hợp hệ cực đo hợp lý”, Tạp chí Các Khoa học Trái đất, 23(3), 217 Phan Đăng Cầu – Phan Thị Hà (2002), “Phƣơng pháp số”, Học viện CNBCVT, tr 66-86 Tạ Văn Đ nh (1995), “Phƣơng pháp tính”, Nhà xuất Giáo dục, tr 66-67 Vũ Đức Minh (2002), “Xử lý tài liệu phƣơng pháp đo sâu phân cực kích thích đối xứng cải tiến”, Tạp chí Các Khoa học trái đất, 24(4), 362 Vũ Đức Minh (2005), “Một đề xuất phƣơng pháp đo sâu điện lƣỡng cực cải tiến”, Tuyển tập cơng trình khoa học, Hội nghị khoa học kỹ thuật Địa Vật lý Việt nam lần thứ IV, 449 Vũ Đức Minh, Đỗ Anh Chung (2015), “Một số kết thử nghiệm xác định đƣờng bão hòa thân đập Đồng Mơ phƣơng pháp Thăm dị điện đa cực tiên tiến”, Tạp chí khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ, 31(3), tr 23-37 Vũ Đức Minh (2008), “Nghiên cứu đề xuất thuật tốn xử lý, phân tích tài liệu phƣơng pháp Phân cực kích thích đa cực cải tiến”, Tạp chí khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ, 24(4), 298-304, Hà Nội Vũ Đức Minh (2010), “Phƣơng pháp Thăm dò điện đa cực cải tiến”, Tạp chí khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ, , 26(4), 233-241 37 Tiếng Anh 10 Advanced Geoscienes (2000-2009), “The SuperSting™ with Swift™ automatic resistivity and IP system Instruction Manual”, Advanced Geosciences inc, Austin, Taxas 11 Advanced Geoscienes (2002), “EarthImager 2D resistivity and IP Invesion”, Advanced Geosciences inc, Austin, Taxas 12 Griffiths D.H., Turnbull J (1985), “A multi-electrode array for resistivity surveying”, First Break, 3, 16 13 Vu Duc Minh, Do Anh Chung (2015), “Introduction to the Advanced Multielectrode Electrical Sounding method”, VNU Journal of MathematicsPhysics, 31(3), p 1-14 14 Vu Duc Minh (2001), “Induced-Polarization Sounding methods in a new manner”, Journal of Geology, Series B, No 17-18, 94 15 Vu Duc Minh (2002), “A new approach for document processing in the improved dipole Induced Polarization sounding method”, VNU, Journal of Science, Natural Sciences and Technology, XVIII (3), 40 38