NGHIÊN CỨU, ỨNG DỤNG PHƯƠNG PHÁP ẢNH ĐIỆN 2D KHẢO SÁT KHẢ NĂNG TÍCH TỤ VÀ DỊCH CHUYỂN Ô NHIỄM TRONG NỀN ĐỊA CHẤT TỪ ÂU THUYỀN THỌ QUANG RA KHU VỰC XUNG QUANH 152 Lương Văn Thọ, Lê Phước Cường, Hồ Hồng[.]
Lương Văn Thọ, Lê Phước Cường, Hồ Hồng Quyên 152 ỨNG DỤNG PHƯƠNG PHÁP ẢNH ĐIỆN 2D TRONG KHẢO SÁT MÔI TRƯỜNG ĐẤT TẠI KHU CÔNG NGHIỆP THỌ QUANG – TP ĐÀ NẴNG GEOLOGICAL SURVEY AT THO QUANG INDUSTRIAL ZONE - DANANG CITY VIA ELECTRICAL GRAPHICS 2-DIMENSIONAL METHOD Lương Văn Thọ1, Lê Phước Cường2, Hồ Hồng Quyên2 Trường Đại học Sư phạm, Đại học Đà Nẵng; tho.luong@gmail.com Trường Đại học Bách khoa, Đại học Đà Nẵng, lpcuong@dut.udn.vn Tóm tắt - Bài báo trình bày kết nghiên cứu thành phần địa chất hai tuyến đo khu công nghiệp Thọ Quang khu dân cư, thành phố Đà Nẵng Tuyến thứ bao gồm 205 điểm liệu tuyến đo có độ dài 165m, khu vực ranh giới khu công nghiệp Thọ Quang âu thuyền Tuyến thứ hai bao gồm 174 điểm liệu thu thập tuyến đo dài 145m theo hướng Bắc-Nam khu vực nghiên cứu khu ranh giới âu thuyền khu vực dân cư Sau xử lý yếu tố gây nhiễu, số liệu định dạng xử lý phần mềm Res2dinv với vòng lặp (tuyến 1) 11 vịng lặp (tuyến 2) thuật tốn sai phân hữu hạn phương pháp bình phương tối thiểu Kết phân tích cho thấy hai tuyến đo phân thành ba tầng địa chất với đặc điểm ô nhiễm khác tuỳ theo độ sâu nghiên cứu Abstract - This paper presents the study results of geological composition measured at two routes between Tho Quang industrial zone and the residential area, Danang city The first route includes 205 measured data of 165m in length which is the boundary between Tho Quang industrial zone and the boat parking area The second route includes 174 measured data of 145m in length north south of the study area which is the boundary between the boat parking area and the residential area After disposing of the confounders, the data was formatted and processed by Res2dinv software with loops (route 1) and 11 loops (route 2) based on finite difference algorithm and the least square methods The analysis results show that both of two routes are divided into geological layers with various contamination characteristics depending on the researched depths Từ khóa - địa chất; ảnh điện 2D; giải đốn; Thọ Quang; ô nhiễm Key words - geology; electrical graphics 2D; interpretation; Tho Quang; contaminated Đặt vấn đề Trong năm gần đây, vấn đề ô nhiễm môi trường nhà máy, khu công nghiệp gây thành phố Đà nẵng quan tâm quan chức nhà khoa học mức độ ảnh hưởng đến chất lượng sống người dân khu vực lân cận Đặc biệt khu công nghiệp dịch vụ thủy sản âu thuyền Thọ Quang, nơi tập trung khu chế xuất hải sản nhà máy sản xuất sản phẩm thủy hải sản xuất Hàng năm trình hoạt động sản xuất thải trực tiếp mơi trường lượng lớn nước thải có nồng độ COD cao chất hữu tồn lưu; chất gây mùi hôi thối phát tán khơng khí, mà qua q trình tích tụ lâu dài làm ô nhiễm môi trường đất nước khu vực xung quanh khu công nghiệp Đặc biệt khu vực âu thuyền xem “bể” chứa chất thải nhà máy, cảng cá, chợ cá khu công nghiệp thủy sản Thọ Quang thải Các độc chất nhiễm có môi trường âu thuyền,theo nước mưa hệ thống mạch nước ngầm tích tụ dịch chuyển khu vực xung quanh, làm ảnh hưởng đến môi trường sinh thái Do đó, lâu dài, để tạo sở cho vấn đề phát triển bền vững khu vực cần phải có khảo sát, điều tra để đánh giá phạm vi, mức độ tích tụ dịch chuyển nhiễm mơi trường đất khu vực ranh giới âu thuyền khu vực dân cư xung quanh khu công nghiệp dịch vụ thủy sản Thọ Quang Từ đề xuất biện pháp quản lý chất lượng môi trường nước, môi trường đất xung quanh khu vực âu thuyền giải pháp xử lý cụ thể, phù hợp với yêu cầu thực tiễn Trước tình hình đó, chúng tơi nhận thấy phương pháp thăm dò ảnh điện 2D tổ hợp phương pháp thăm dò điện nội dung địa vật lý môi trường Đây phương pháp với ưu điểm khơng xâm thực, thiết bị máy móc đo đạc gọn nhẹ dễ thu thập số liệu thực địa Cùng với hỗ trợ thuật tốn chương trình có khả xử lý khối lượng lớn liệu thu thập phương pháp trở thành công cụ hữu hiệu khảo sát, đánh giá mức độ ô nhiễm, tai biến, rủi ro môi trường địa chất [1, 2] Từ mục đích khảo sát mơi trường địa chất, nhằm phát dị thường mức độ tích tụ dịch chuyển yếu tố ô nhiễm môi trường đất xung quanh khu vực âu thuyền Thọ Quang, chúng tơi thực cơng trình nghiên cứu: “Ứng dụng phương pháp ảnh điện 2D khảo sát môi trường đất khu vực âu thuyền Thọ Quang – TP Đà Nẵng” Đối tượng và phương pháp 2.1 Đối tượng Tham số quan trọng thăm dò trường điện không đổi đất đá giá trị điện trở suất (trong hệ đo lường SI, điện trở suất ký hiệu làρ, thứ nguyên Ohm.m) Tham số điện trở suất đặc trưng tính chất dẫn điện mơi trường, tính chất phụ thuộc vào thành phần khoáng vật, thạch học, cấu trúc, điều kiện lịch sử tạo thành nằm đất đá Để có nhìn khách quan tính chất dẫn điện đất, đá mặt đất, Keller, Frischknecht (1966) Daniels, Alberty (1966) đưa bảng số liệu thực nghiệm điện trở suất trình bày Bảng [3] Bảng Điện trở suất số đất, đá, khống vật hóa chất phổ biến [3] Vật liệu Điện trở suất (Ωm) Độ dẫn điện (1/Ωm) Nham thạch đá biến chất: - Granite (đá granit) 5.103 106 10-6 2.10-4 ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(90).2015 - Basalt (đá bazan) - Slate (đá phiến) 6.102 4.107 2,5.10-8 1,7.10-3 - Marble (đá cẩm thạch) 102 2,5.108 4.10-9 10-2 - Quartzite (thạch anh) 10 2.10 5.10-9 10-2 Trầm tích: - Sandstone (sa thạch) 4.103 2,5.10-4 0,125 - Shale (đá phiến sét) 20 2.103 5.10-4 0,05 - Limestone (đá vôi) 50 4.102 2,5.10-3 0,02 Đất nước: - Clay (đất sét) 100 0,01 - Alluvium (đất phù sa) 10 800 1,25.10-3 0,1 - Goundwater (nước ngầm) 10 100 0,01 0,1 - Sea water (nước biển) 0,2 Hóa chất: - Iron (sắt) 9,074.10-8 1,102.107 0,708 - 0,01 phân tử gam KCl 1,413 0,843 - 0,01 phân tử gam NaCl 1,185 6,13 - 0,01M axit Axetic 0,163 6,998.1016 - Xylene 1,429.10-17 103 106 10-6 10-3 Điện trở suất đá xâm nhập biến chất thường có giá trị cao Giá trị điện trở suất loại đá phụ thuộc nhiều vào độ nứt nẻ mức độ chứa nước đới nứt nẻ Giá trị điện trở suất ứng với loại đất đá thay đổi giới hạn rộng, từ hàng triệuΩ.m đến nhỏ mộtΩ.m.Các đá trầm tích thường có độ xốp độ chứa nước cao nên có giá trị điện trở suất thấp so với đá thâm nhập đá biến chất, giá trị điện trở suất đá thường thay đổi khoảng từ10Ω.m đến 10000Ω.m, hầu hết có giá trị nhỏ 1000Ω.m Giá trị điện trở suất phụ thuộc lớn vào độ xốp, độ chứa nước đá đặc biệt độ khoáng hóa nước chứa lỗ rỗng [4] Các trầm tích bở rời khơng gắn kết thường có giá trị điện trở suất thấp đá trầm tích với giá trị thay đổi từ vài Ω.m đến nhỏ 1000Ω.m Giá trị điện trở suất chúng phụ thuộc vào độ xốp (các trầm tích chứa nước bão hịa) hàm lượng khống vật sét, đất sét thường có giá trị điện trở suất thấp đất cát Giá trị điện trở suất nước đất dao động khoảng từ 10 đến 100 Ω.m, phụ thuộc vào hàm lượng muối hồ tan có đất Chú ý rằng, điện trở suất nước biển thấp (khoảng 0,2Ω.m) hàm lượng muối cao Điều giúp cho phương pháp thăm dò điện trở thành kỹ thuật lý tưởng việc đo vẽ đồ xác định ranh giới nhiễm mặn vùng duyên hải [5] Giá trị điện trở suất số loại vật liệu hóa chất nhiễm cơng nghiệp trình bày Bảng Một số kim loại sắt có giá trị điện trở suất thấp Các hoá chất điện phân mạnh KCl NaCl làm giảm cách đáng kể điện trở suất nước đất đến giá trị nhỏ 1Ω.m hóa chất có hàm lượng tương đối thấp Đó đặc tính giúp ta khảo sát khả tích tụ dịch chuyển nhiễm hóa chất cơng nghiệp địa chất [6, 7] 2.2 Phương pháp nghiên cứu 2.2.1 Nghiên cứu thuyết ảnh điện 2D Trong thăm dò ảnh điện 2D, định luật Ohm chi phối truyền dẫn dòng điện môi trường địa chất Dạng vi 153 phân định luật Ohm môi trường liên tục: J = (−U ) = E (1) Trong phương pháp thăm dò điện, người ta thường dùng = điện dẫn suất (hay độ dẫn điện) môi trường, J mật độ dòng diện dẫn điểm xét môi trường, E cường độ điện trường điểm quan sát môi trường Trong hầu hết phương pháp thăm dò điện, nguồn điện thường có dạng nguồn điểm Trong trường hợp này, xét phần tử tích V bao quanh nguồn dịng điện I vị trí (x I , y I , z I ) , phương trình mơ tả quan hệ cường độ dịng mật độ dịng có dạng: I .J = ( x − x I ) ( y − y I ) ( z − z I ) V (2) (2) công thức mà Dey Morrison sử dụng thăm dị ảnh điện 2D năm 1979, đó, hàm Delta Dirac với tính chấtsau: + , x = • ( x) = ; (− x) = ( x) 0, x + + − − ( x)dx = 1; U ( x) ( x)dx = U (0), U ( x) + • U ( x) ( x − x I )dx = U ( x I ), U ( x) − Từ (1) (2) ta viết lại: − (x, y, z )U ( x, y, z ) = I = (x − x I ) ( y − y I ) (z − z I ) V (3) Phương trình (3) phương trình mô tả phân bố điện môi trường nguồn dịng điểm gây Có nhiều kỹ thuật phát triển để giải phương trình cho tốn thuận tốn ngược thăm dị trường điện không đổi 2.2.2 Nghiên cứu thực nghiệm a Nghiên cứu cấu hình thiết bị WennerSchlumberger đo điện trở suất biểu kiến Sự phân bố điện trở suất đặc điểm cấu trúc địa chất môi trường bên mặt đất tạo dáng điệu hay trường điện riêng bên bề mặt Do phương pháp thăm dị ảnh điện 2D trường điện không đổi, muốn biết thông tin môi trường địa chất bên dưới, ta phải tiến hành phép đo điện trở suất biểu kiến bề mặt Phép đo thực cách phát dịng điện khơng đổi có cường độ I vào môi trường địa chất cần khảo sát thông qua điện cực Đối với cấu hình thiết bị Wenner-Schlumberger bốn cực đối xứng (theo Hình 1), hiệu điện hai điện cực thu P1, P2 tính sau: Lương Văn Thọ, Lê Phước Cường, Hồ Hồng Quyên 154 U = U ( P1 ) − U ( P2 ) = = I 2 1 1 − − + rC P rC P rC P r C2 P2 1 11 Trong đó, rC P , rC P , rC P , rC P 1 2 2 (4) - Quy trình đo thực địa hai vị trí Âu thuyền khoảng cách từ điện cực thu P1, P2 đến điện cực dòng C1, C2 Từ cơng thức (4), ta xác định điện trở suất biểu kiến môi trường địa chất bên dưới: a = k Trong đó: k= U I (5) 2 1 1 r − r − r + r C P CP C P CP 1 1 2 165m); âu thuyền khu dân cư (với chiều dài 145m)được trình bày Hình (là tham số hình học phụ thuộc vào xếp điện cực) Điện trở suất biểu kiến điện trở suất thật môi trường địa chất bên Liên hệ giá trị điện trở suất biểu kiến giá trị điện trở suất thật bên môi trường địa chất mối liên hệ phức tạp Việc xác định điện trở suất thật từ giá trị điện trở suất biểu kiến quan sát vấn đề toán ngược thăm dị ảnh điện 2D, dựa vào thuật tốn sai phân hữu hạn Hình Cấu hình thiết bị Wenner-Schlumberger b Nghiên cứu quy trình đo ngoài thực địa cấu hình thiết bị Wenner-Schlumberger các vị trí khu vực âu thuyền Thọ Quang * Tại ranh giới khu công nghiệp dịch vụ thủy sản Thọ Quang âu thuyền Quy trình đo đạc thực nghiệm trình bày theo sơ đồ Hình Cách bố trí điện cực hệ thiết bị giống hệ thiết bị Wenner-alpha, hai cực thu P1, P2 nằm hai cực phát C1, C2 giữ khoảng cách “a” không đổi suốt trình đo mức đo Số phép đo (số điểm đo) mức đo tính theo biểu thức tổng quát [m-(2n+1)] ứng với bước dịch chuyển điểm liệu theo khoảng cách điện cực đơn vị Trước tiên, tiến hành cắm 34 điện cực dọc theo tuyến đo với chiều dài 165m, khoảng cách điện cực đơn vị a = 5m Từ mức đo sâu thứ đến thứ sáu, khoảng cách điện cực a1 = 5m thừa số n thay đổi từ đến 6: Tại mức đo sâu thứ ứng với thiết bị thứ (n = 1), khoảng cách điện cực (C1P1 = P1P2 = P2C2 = a = 5m) giữ không đổi suốt trình đo Các điện cực phát C1, C2 thu P1, P2 kết nối với nguồn máy đo thông qua cuộn cáp, liệu thu thập từ đầu cuối tuyến đo, tổng số điểm liệu ghi mức đo thứ 31 điểm Ở mức đo sâu thứ hai (thừa số n = 2), khoảng cách điện cực (C1P1 = C2P2 = 2a = 10m, P1P2 = a = 5m) giữ không đổi suốt trình đo Các điểm liệu thu thập từ đầu cuối tuyến với tổng số điểm liệu ghi 29 Tương tự, thực thu thập liệu đến mức liệu thứ sáu (n = 6) Đối với mức đo sâu thứ bảy tám, khoảng cách điện cực thiết bị tăng lên a2 = 10m thừa số độ sâu n thay đổi từ đến Ở mức đo sâu thứ bảy (a2 = 10m, n = 3), khoảng cách điện cực C1P1 = C2P2 = 3a2 = 30m P1P2 = a2 = 10m Ở mức có 20 điểm liệu thu thập từ đầu tuyến cuối tuyến Ở mức đo sâu thứ tám (a2 = 10m, n = 4), khoảng cách điện cực C1P1 = C2P2= 4a2 = 40m P1P2 = a2 = 10m, số điểm liệu thu thập từ đầu tuyến cuối tuyến 16 điểm Đối với mức đo sâu cuối cùng, khoảng cách thiết bị a3 = 15m thừa số n = 3, ứng với khoảng cách (C1P1 = C2P2 = 3a3 = 45m, P1P2 = a3 = 15m) số điểm liệu ghi từ đầu cuối tuyến 13 điểm Thiết bị có độ nhạy cho hai cấu trúc: phân bố ngang ứng với thừa số “n” thấp phân bố thẳng đứng ứng với “n” cao, cường độ tín hiệu yếu thiết bị Wenner, lớn thiết bị lưỡng cực gấp hai lần thiết bị Pole-dipole, mức độ bao phủ ngang rộng so với thiết bị Wenner, hẹp so với thiết bị lưỡng cực * Tại khu vực ranh giới âu thuyền khu dân cư Hình Vị trí tuyến đo hai khu vực ranh giới âu thuyền - Vị trí tuyến đo Theo khảo sát thực địa, để có sở đánh giá tích tụ dịch chuyển nhiễm môi trường địa chất từ âu thuyền Thọ Quang khu vực xung quanh, hai tuyến đo lập hai ranh giới bản: khu công nghiệp dịch vụ thủy sản Thọ Quang âu thuyền (với chiều dài Sơ đồ bố trí đo đạc theo Hình 3, số điện cực cắm m = 29 (điện cực) cách với khoảng cách đơn vị a = 5m chiều dài tuyến đo 145m Từ mức liệu thứ đến thứ sáu (ứng với a1 = 5m n = 1, ,6) mức thứ có 27 điểm liệu ghi, mức thứ hai 25 điểm liệu tiếp tục đến mức thứ sáu 17 điểm liệu ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(90).2015 155 Mức liệu thứ bảy tám ứng với a2 = 10m n = 3, Mức thứ bảy có 16 điểm liệu mức thứ tám có 12 điểm liệu ghi Mức liệu thứ chín (ứng với a3 = 15m n = 3) có điểm liệu ghi, mức thứ mười (ứng với a4 = 25m n = 2) ta ghi điểm liệu Hình Kết ảnh điện 2D ranh giới âu thuyền khu dân cư, xử lý phần mềm Res2dinv 3.1.3 Kết hai tuyến biểu diễn hệ trục Để có nhìn bao qt tương quan kết hai tuyến đo, ta biểu diễn mặt cắt ảnh điện hai tuyến hệ trục tọa độ theo Hình Hình Cách bố trí điện cực quy trình thực phép đo để xây mặt cắt ảnh điện 2D cho hệ thiết bị Wenner-Schlumberger Kết nghiên cứu và thảo luận 3.1 Kết nghiên cứu 3.1.1 Tuyến đo khu vực ranh giới khu công nghiệp dịch vụ thủy sản Thọ Quang âu thuyền Tại khu vực này, có khoảng 205 điểm liệu thu thập tuyến đo với chiều dài khoảng 165m Sau xử lý liệu bị nhiễu, liệu ổn định định dạng xử lý phần mềm Res2Dinv với vòng lặp thuật tốn sai phân hữu hạn phương pháp bình phương tối thiểu Kết biểu diễn dạng ảnh điện hai chiều với sai số khoảng 7,2% (Hình 4) Hình Biểu diễn kết hai tuyến đo hệ trục 3.2 Thảo luận giải đoán kết 3.2.1 Tuyến đo (khu vực ranh giới KCN dịch vụ thủy sản Thọ Quang âu thuyền) Kết ảnh điện hai chiều tuyến (Hình 4) cho thấy cấu trúc phân bố địa chất bên (đến độ sâu nghiên cứu) phân thành ba tầng: Hình Mặt cắt ảnh điện 2D khu vực ranh giới khu công nghiệp dịch vụ thủy sản Thọ Quang âu thuyền 3.1.2 Tuyến đo khu vực ranh giới âu thuyền khu dân cư Có 174 điểm liệu thu thập tuyến đo dài 145m theo hướng bắc-nam độ sâu nghiên cứu khoảng 24m Sau xử lý yếu tố gây nhiễu, số liệu định dạng, xử lý phần mềm Res2Dinv thuật tốn sai phân hữu hạn với 11 vịng lặp theo phương pháp bình phương tối thiểu kết ảnh điện thu với sai số 9,7% Kết biểu diễn dạng ảnh điện hai chiều theo Hình + Tầng địa chất thứ phân bố từ mặt đất đến độ sâu khoảng 8m dọc theo tuyến đo, giá trị điện trở suất thay đổi từ 3,32Ωm đến 94Ωm, thành phần vật chất tầng giải đoán đất cát, sạn vụn đất sét (do bồi đắp từ nơi khác đến) Tầng có mật độ chứa nước lớn, chủ yếu phân bố độ sâu khoảng 5m dọc theo tuyến đo Có thể khoảng độ sâu này, số vị trí dọc theo tuyến đo tồn mạch nước ngầm liên thông đến âu thuyền Kết cho thấy nước có dấu hiệu bị nhiễm, cụ thể vị trí khoảng từ 10m đến 21m, 25m đến 43m, 50m đến 58m, 95m đến 105m, 125m đến cuối tuyến đo (trong khoảng giá trị điện trở suất nhỏ 10Ωm, chí giảm đến 3,32Ωm) + Tầng thứ hai phân bố độ sâu từ 8m đến 15,9m có giá trị điện trở suất vào khoảng 100Ωm đến 873Ωm, thành phần gồm đất đá phiến sét, cát kết, sa thạch vụn trộn lẫn bùn đen Tầng có mật độ chứa nước tầng thứ nhất, hình thức tồn chủ yếu nước bề mặt đất đá, dấu hiệu nước ngầm Độ ẩm tầng nhỏ tầng thứ lớn tầng thứ ba Lương Văn Thọ, Lê Phước Cường, Hồ Hồng Quyên 156 + Tầng thứ ba phân bố độ sâu từ 15,9m đến hết độ sâu nghiên cứu (có thể vượt độ sâu nghiên cứu), giá trị điện trở suất lớn vào khoảng 873,04Ω.m-32654,30Ω.m Theo giải đoán, thành phần vật chất chủ yếu đá phiến, đá cuội tảng lớn Có thể tầng đá gốc thuộc hệ thống bờ kè đá hành lang âu thuyền trước đây, có tác dụng ngăn không cho sạt lở đất hay xâm thực dòng chảy ngầm từ âu thuyền 3.2.2 Tuyến đo (khu vực ranh giới âu thuyền khu dân cư) Theo kết khảo sát ảnh điện 2D khu vực ranh giới âu thuyền khu dân cư cấu trúc phân bố địa chất đến độ sâu nghiên cứu phân thành ba tầng: + Tầng địa chất thứ phân bố khoảng từ mặt đất đến độ sâu 6m dọc theo tuyến đo có giá trị điện trở suất thay đổi khoảng từ 0,869Ωm đến 100Ωm Tầng có mật độ chứa nước lớn phân bố gần bề mặt (trong số trường hợp gây nên hiệu ứng bẫy điện, gây nhiễu q trình đo đạc) Ngồi theo miền giá trị điện trở suất thành phần vật chất tầng giải đoán đất sét trộn lẫn đất cát với sạn vụn (phân bố phía trên) đất bùn-sét ngậm nước (phân bố phía dưới, có độ dày vào khoảng 3m đến 6,5m) Điểm đáng ý từ bề mặt đến độ sâu khoảng 3m, ta nhận thấy có dấu hiệu tích tụ yếu tố nhiễm phân bố dọc theo tuyến đo.Cụ thể, vị trí khoảng từ 7m đến 30m, 36m đến 40m, 55m đến 65m, 80m đến 95m, 117m đến 136m dọc theo tuyến đo; khoảng giá trị điện trở suất nhỏ 6,96Ω.m, chí giảm đến 0,869Ω.m Đây yếu tố nhiễm có nguồn nước âu thuyền, thẩm thấu qua hệ thống đê kè ngấm vào lớp địa chất khu vực Có thể giải thích điều rằng, âu thuyền nơi bị ô nhiễm nghiêm trọng nguồn nước thải từ KCN dịch vụ thủy sản Thọ Quang nước thải từ cảng cá, chợ cá khu vực lân cận + Tầng địa chất thứ hai phân bố độ sâu khoảng từ 6,5m đến 15m dọc theo tuyến đo, giá trị điện trở suất vào khoảng 150Ωm đến 400Ωm Theo miền giá trị điện trở suất, thành phần vật chất lớp giải đoán cát kết, sa thạch, trộn lẫn với đá phiến sét Kết lý giải q trình bồi đắp từ nơi khác đến trình san lấp mở rộng hành lang bờ kè âu thuyền để xây dựng khu dân cư Tầng có mật độ nước so với tầng trên, nhiên độ ẩm cao thay đổi mạnh theo mùa (đặc biệt mùa khô mùa mưa) thể tồn lớp nước mặt đất đá tầng địa chất khơng có dấu hiệu nước ngầm Điều lý giải che chắn hệ thống đê kè gia cố tốt hành lang bên (giữa âu thuyền khu vực địa chất độ sâu tương ứng) Điểm đáng ý tầng khoảng từ 55m đến 85m (dọc theo tuyến đo) bị lõm sâu xuống tầng vượt độ sâu nghiên cứu, dấu hiệu sụp lún tầng địa chất vị trí + Tầng địa chất thứ ba phân bố độ sâu khoảng từ 15,9m đến hết độ sâu nghiên cứu, giá trị điện trở suất tầng từ 447Ωm đến 1264Ωm) Thành phần vật chất giải đoán đất đá bazan, đá phiến, đá cuội tảng lớn Có thể tầng đá gốc thuộc hệ thống bờ kè đá hành lang âu thuyền trước đây, có tác dụng ngăn khơng cho sạt lở đất hay xâm thực dòng chảy ngầm âu thuyền trình bồi đắp để quy hoạch xây dựng khu dân cư bị chôn vùi sâu bên Tầng cứng chắc, có độ ẩm thay đổi mạnh theo mùa, nhiên khơng có dấu hiệu nước ngầm liên thơng với âu thuyền 3.2.3 Thảo luận giải đoán tổng hợp hai tuyến Về cấu trúc phân bố địa chất thành phần vật chất hai khu vực khảo sát âu thuyền có tương đồng Tuy nhiên, độ dày tầng địa chất hai vị trí khảo sát tương đối khác Ở tầng địa chất thứ hai khu vực khảo sát có mật độ chứa nước cao nước có dấu hiệu nhiễm Điều làm tăng nguy tích tụ nhiễm môi trường địa chất tầng thứ hai khu vực Theo liệu ảnh điện cho thấy mức độ ô nhiễm tầng địa chất thứ tuyến đo lớn tuyến đo Miền giá trị điện trở suất tầng địa chất thứ tuyến thứ giảm đến 0,869Ωm, tuyến giảm đến 3,32Ωm Điều lý giải tuyến đo thứ (khu vực ranh giới âu thuyềnvà khu dân cư) thực gần với hành lang đê kè âu thuyền so với tuyến 1, nước từ khu vực âu thuyền thẩm thấu trực tiếp qua hành lang âu thuyền dễ dàng Vị trí tuyến đo thứ hai (tại khu vực ranh giới KCN dịch vụ thủy sản Thọ Quang âu thuyền) cách xa hành lang bờ kè âu thuyền hơn, nên yếu tố ô nhiễm từ nước âu thuyền dịch chuyển đến vị trí chủ yếu theo mạch nước ngầm mức độ nhiễm giảm dần theo khoảng cách dịch chuyển Trong tầng địa chất thứ ba, hai vị trí quan sát dấu tích hệ thống đê kè vững âu thuyền khơng có dấu hiệu nước từ âu thuyền thẩm thấu hay xâm thực qua hệ thống Tuy nhiên, vị trí tuyến khảo sát gần âu thuyền nên hành lang đê kè có độ dày lớn vị trí tuyến 3.3 Kiến nghị Kết khảo sát ảnh điện 2D hai khu vực ranh giới âu thuyền cho thấy có dấu hiệu tích tụ dịch chuyển yếu tố nhiễm môi trường đất tầng địa chất gần bề mặt (phân bố độ sâu khoảng từ 2m đến 5m) xung quanh khu vực âu thuyền Các yếu tố nhiễm lan rộng thấm sâu xuống tầng địa chất bên Về lâu dài, khơng có biện pháp xử lý vấn đề ô nhiễm ảnh hưởng đến môi trường sống Đặc biệt vào mùa nắng hạn, khí từ đất có chứa yếu tố ô nhiễm tác động không tốt đến sức khỏe người dân khu dân cư, chưa kể mùi thối khơng khí bốc lên trực tiếp từ âu thuyền Do đó, để khắc phục vấn đề này, cần có biện pháp cụ thể kiểm sốt có hiệu việc xả thải KCN dịch vụ thủy sản Thọ Quang, cảng cá, chợ cá trình hoạt động sản xuất Cần xử lý kịp thời sở vi phạm việc xả thải theo quy định hành bảo vệ môi trường Thực kiểm tra bảo trì, nâng cấp cơng nghệ công suất cho trạm xử lý nước thải KCN Thọ Quang Đồng thời kiểm tra, gia cố lại hệ thống đê kè xung quanh khu vực âu thuyền để ngăn không cho xâm thực, thẩm thấu nước từ âu thuyền khu vực dân cư xung quanh Kết luận Kết nghiên cứu thành phần địa chất khu cơng ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(90).2015 nghiệp Thọ Quang – TP Đà nẵng khu dân cư cho thấy kết cấu địa chất phân thành ba tầng, với đặc điểm ô nhiễm khác tuỳ theo độ sâu nghiên cứu Hai khu vực nghiên cứu (tuyến đo 2) có cấu trúc phân bố địa chất tương đồng độ dày tầng địa chất khác Kết thể rõ phương pháp ảnh điện 2D công cụ hữu hiệu dùng để khảo sát, đánh giá đưa tranh tổng quát môi trường địa chất khu vực nghiên cứu Thơng qua đó, ta biết thành phần đất đá, nước ngầm, đặc điểm cấu trúc, điều kiện hình thành, kiến tạo môi trường địa chất Bên cạnh đó, giá thành đợt khảo sát ảnh điện 2D thấp so với hướng thăm dò khác tổ hợp phương pháp địa vật lý nên cần triển khai mở rộng ứng dụng nhiều Việt Nam Lời cảm ơn Nghiên cứu tài trợ kinh phí từ Bộ Giáo dục Đào tạo Việt Nam khuôn khổ đề tài KHCN mã số B2014-01-16 157 TÀI LIỆU THAM KHẢO [1] Nguyễn Ngọc Thu (2006), Phương pháp thăm dò điện 2D, Liên đoàn Bản đồ Địa chất Miền Nam [2] Nguyễn Thành Vấn, Lê Ngọc Thanh, Nguyễn Minh Anh (2005), “Áp dụng phương pháp ảnh điện để nghiên cứu tính chất bất đồng mơi trường đất đá”, Tạp chí Phát triển Khoa học Công nghệ, tập 8, ĐHQG Tp.HCM, pp 35-42 [3] Dey, A and Morrison, H.F (1979), “Resistivity modelling for arbitrary shaped two dimensional structures”, Geophysical Prospecting, (No.27), pp1020-1036 [4] Loke M.H and Barker R.D (1995), Improvements to the Zohdy method for the inversion of resistivity sounding and pseudesection data, Computers and Geosciences, (Vol21,No.2), pp 321-322 [5] Olayinka A.I and Yaramanci U (2000), Use of block Inversion in the 2D interpretation of apparent resistivity data and its comparision with smoth inversion, Journal of Apply Geophysics, 45, pp 403-416 [6] McGillvray P.R and Oldenburg D.W (1990), Methods for calculating Frechet derivatives and sensitivities for he non-linear inverse problem, A comparative study, Geophysical Prospecting, 38,pp 499-524 [7] Sasaki Y (1992), Resolution of resistivity tomography inferred from numerical simulation, Geophysical Prospecting, 40, pp 453-464 (BBT nhận bài: 16/03/2015, phản biện xong: 21/04/2015) ... (khu vực ranh giới âu thuyền khu dân cư) Theo kết khảo sát ảnh điện 2D khu vực ranh giới âu thuyền khu dân cư cấu trúc phân bố địa chất đến độ sâu nghiên cứu phân thành ba tầng: + Tầng địa chất. .. Tại khu vực ranh giới âu thuyền khu dân cư Hình Vị trí tuyến đo hai khu vực ranh giới âu thuyền - Vị trí tuyến đo Theo khảo sát thực địa, để có sở đánh giá tích tụ dịch chuyển nhiễm môi trường địa. .. giúp ta khảo sát khả tích tụ dịch chuyển ô nhiễm hóa chất công nghiệp địa chất [6, 7] 2.2 Phương pháp nghiên cứu 2.2.1 Nghiên cứu thuyết ảnh điện 2D Trong thăm dò ảnh điện 2D, định luật Ohm chi