Giới thiệu dự án

Trong bối cảnh đô thị hóa nhanh chóng tại các đô thị ven biển như Thành phố Đà Nẵng, việc xây dựng các công trình cao tầng, tầng hầm sâu và hạ tầng ngầm đòi hỏi công tác khảo sát địa chất công trình phải đạt độ chính xác và độ tin cậy cao. Khu vực Quận Ngũ Hành Sơn (dọc lưu vực sông Đô Tỏa) có cấu trúc địa tầng đan xen phức tạp giữa cát bở rời, bùn sét bão hòa nước và đá phiến/đá vôi phong hóa cục bộ. Theo thống kê ngành địa kỹ thuật, hơn 65% sự cố lún nứt hố móng và sạt trượt bờ sông bắt nguồn từ việc đánh giá sai lệch ranh giới dị thường địa chất ngầm do khoảng cách các điểm khoan thăm dò quá thưa.

Phương pháp khảo sát địa kỹ thuật truyền thống chủ yếu dựa vào khoan thăm dò lấy mẫu lõi. Dù cho kết quả trực quan tại điểm khoan, phương pháp này có nhược điểm lớn: chi phí cao, thời gian thi công kéo dài và hoàn toàn "mù thông tin" tại các khoảng không gian giữa hai lỗ khoan. Ngược lại, phương pháp đo sâu điện 1D (Vertical Electrical Sounding - VES) cổ điển lại bỏ qua hoàn toàn sự biến thiên điện trở suất theo phương ngang. Do đó, việc nghiên cứu ứng dụng phương pháp chụp cắt lớp ảnh điện 2D/3D (Electrical Resistivity Tomography - ERT) kết hợp đối sánh trực tiếp với khoan thăm dò là một giải pháp khoa học - kỹ thuật cấp thiết.

Mục tiêu nghiên cứu

  1. Tổng quan hóa cơ sở lý thuyết truyền dẫn điện trong môi trường bán không gian phân lớp dị hướng và mô hình hóa trường điện thế.
  2. Xây dựng quy trình xử lý bài toán thuận (Forward Modeling) và bài toán ngược (Inverse Modeling) bằng phương pháp bình phương tối thiểu có ràng buộc độ trơn (Smoothness-constrained Least-Squares).
  3. Đánh giá độ nhạy (Frechet derivative) và triển khai cấu hình đo Wenner-Alpha tại thực địa khu vực ngã ba đường Bùi Tá Hán – Chương Dương, Quận Ngũ Hành Sơn.
  4. Xử lý nghịch đảo dữ liệu đo thực nghiệm 2D/3D bằng phần mềm chuyên dụng Res2Dinv/Res3Dinv và kiểm chứng độ chính xác với kết quả khoan địa chất thực tế.

Phương pháp tiếp cận và kết quả kỳ vọng

Giải pháp tiếp cận tích hợp phương pháp địa vật lý điện trở suất đa cực (ERT) và khoan kiểm chứng địa kỹ thuật. Dữ liệu trường điện được đo đạc tự động trên mặt đất, xử lý số học qua phương pháp Gauss-Newton cải biên nhằm tái tạo mặt cắt điện trở suất liên tục theo chiều sâu và phương ngang.

Chỉ số định lượng kỳ vọng:

  • Độ sâu khảo sát hữu hiệu ($Z_e$): $\ge 15\text{ m}$ với chiều dài tuyến đo $L = 60 - 120\text{ m}$.
  • Sai số bình phương trung bình (RMS error) của mô hình nghịch đảo: $< 5.0%$.
  • Độ phân giải lớp địa chất: Phát hiện chính xác thấu kính sét và mực nước ngầm với độ dày lớp $\ge 0.5\text{ m}$.
  • Độ tương đồng địa tầng giữa ảnh điện và hố khoan kiểm chứng: $\ge 90%$.

Phạm vi và giới hạn đề tài

  • Không gian khảo sát: Khu đất góc ngã ba Bùi Tá Hán – Chương Dương, ven sông Đô Tỏa, Quận Ngũ Hành Sơn, TP. Đà Nẵng.
  • Giới hạn kỹ thuật: Khảo sát tập trung vào tầng nông phục vụ địa kỹ thuật công trình (độ sâu $< 30\text{ m}$), áp dụng cấu hình Wenner-Alpha; không xét các cấu trúc địa chất siêu sâu hoặc đứt gãy kiến tạo khu vực.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí Đo sâu điện 1D (VES) Khoan thăm dò địa chất Cắt lớp ảnh điện 2D/3D (ERT)
Độ phân giải ngang Kém (coi môi trường phân lớp ngang lý tưởng) Điểm rời rạc (không có tính liên tục) Cao (liên tục dọc tuyến và lưới quét không gian)
Độ phân giải đứng Trung bình Rất cao (lấy mẫu lõi trực tiếp) Cao (tính toán qua hàm độ nhạy Frechet)
Chi phí triển khai Thấp Rất cao ($400.000 - 800.000\text{ VNĐ/m}$ khoan) Trung bình (tối ưu hóa theo ca đo)
Thời gian khảo sát Nhanh ($1 - 2\text{ giờ/điểm}$) Chậm ($1 - 3\text{ ngày/hố}$ khoan sâu) Rất nhanh ($30 - 45\text{ phút/mặt cắt}$ $60\text{ m}$)
Tác động môi trường Không phá hủy Gây ồn, xáo trộn tầng đất, bùn thải Không phá hủy (Non-destructive testing)

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have (Bắt buộc): Đo tự động đa cực với hệ số hình học $k$ chính xác; thuật toán nghịch đảo 2D hội tụ RMS $< 5%$; xác định rõ ranh giới cát ẩm, sét bùn và đá gốc phong hóa.
  • Should have (Nên có): Xuất lưới tọa độ không gian 3D tương quan với mực nước sông Đô Tỏa; tính toán hàm độ nhạy 1D/2D phục vụ tiền thiết kế khoảng cách điện cực $a$.
  • Could have (Có thể có): Liên kết dữ liệu địa điện với số búa thí nghiệm xuyên tiêu chuẩn (Standard Penetration Test - SPT) từ hố khoan.
  • Won't have (Không làm đợt này): Quan sát biến thiên điện trở suất theo thời gian 4D (Time-lapse ERT).

Thiết kế hệ thống

Technology Stack & Công cụ sử dụng

  • Phần cứng thu thập dữ liệu: Hệ thiết bị đo địa điện tự động đa cực, bộ chuyển mạch điện cực tử (Electrode Switch Box), cáp đo chuyên dụng 24–48 kênh, điện cực thép không gỉ.
  • Phần mềm mô hình hóa & nghịch đảo:
    • Res2Dmod v3.01: Giải bài toán thuận bằng phương pháp sai phân hữu hạn (Finite Difference Method - FDM) và phần tử hữu hạn (Finite Element Method - FEM).
    • Res2Dinv v3.59: Nghịch đảo 2D phi tuyến theo thuật toán Gauss-Newton làm trơn ma trận.
    • Res3Dinv v2.15: Xử lý khối dữ liệu 3D biểu diễn hệ tọa độ $OXYZ$.
    • Surfer v15.0: Biên tập bản đồ đẳng trị và bình đồ phân bố địa điện.
    • Python 3.8 (NumPy, SciPy, Matplotlib): Xử lý thống kê và kiểm tra sai số.

Methodology & Kế hoạch triển khai

Tuần 1-2: Khảo sát hiện trường & Thiết kế mạng lưới điện cực
Giai đoạn Nội dung công việc Đầu ra (Deliverables)
Giai đoạn 1: Chuẩn bị Khảo sát địa hình ven sông Đô Tỏa, đo điện trở tiếp xúc ($R_{ct} < 2\text{ k}\Omega$). Sơ đồ mạng lưới 3 tuyến đo cắt ngang và dọc khu đất.
Giai đoạn 2: Đo thực địa Rải điện cực cấu hình Wenner-Alpha ($a = 2\text{ m}$, $L = 60\text{ m}$), phát dòng $I = 50 - 200\text{ mA}$. Tập tin dữ liệu thô .DAT ghi nhận hàng trăm điểm đo $\rho_a$.
Giai đoạn 3: Xử lý số liệu Lọc nhiễu công nghiệp, nội suy sai phân hữu hạn, chạy lặp Gauss-Newton. File mặt cắt phân bố điện trở thực 2D, mô hình 3D OXYZ.
Giai đoạn 4: Khoan đối chứng Khoan xoay 01 hố khoan tại Tuyến 3, lấy mẫu nguyên dạng ở các độ sâu. Cột địa tầng thạch học chi tiết đến độ sâu $-15\text{ m}$.

Implementation và kết quả

Development Process & Cơ sở toán học

Phương trình vi phân chi phối phân bố điện thế trong môi trường liên tục xuất phát từ định luật Ohm vi phân và phương trình liên tục:

$$\nabla \cdot \vec{J} = -\nabla \cdot (\sigma \nabla U) = \sum I_i \delta(\vec{r} - \vec{r}_i)$$

Trong đó $\sigma = \frac{1}{\rho}$ là độ dẫn điện $(\Omega^{-1}\cdot\text{m}^{-1})$, $U$ là điện thế $(\text{V})$, $I$ là dòng phát $(\text{A})$, và $\delta$ là hàm Dirac delta.

Đối với cấu hình Wenner-Alpha với 4 cực $C_1, P_1, P_2, C_2$ cách đều nhau khoảng cách $a$, hệ số hình học $k$ được tính:

$$k = \frac{2\pi}{\frac{1}{r_{C1P1}} - \frac{1}{r_{C2P1}} - \frac{1}{r_{C1P2}} + \frac{1}{r_{C2P2}}} = 2\pi a$$

Điện trở suất biểu kiến:

$$\rho_a = 2\pi a \frac{\Delta U}{I}$$

Trong bài toán ngược 2D, phương trình Gauss-Newton làm trơn ma trận (Smoothness-constrained Least-Squares) có dạng:

$$(J^T J + \lambda F) \Delta q_k = J^T g - \lambda F q_k$$

Trong đó:

  • $J$: Ma trận Jacobi chứa các đạo hàm riêng phần (Frechet derivatives) $J_{ij} = \frac{\partial f_i}{\partial q_j}$.
  • $F = \alpha_x C_x^T C_x + \alpha_z C_z^T C_z$: Ma trận làm trơn theo phương ngang và phương đứng.
  • $\lambda$: Tham số điều hòa Marquardt (Damping factor).
  • $g = y - f(q)$: Vector độ lệch giữa số liệu logarit điện trở suất đo được và đáp ứng mô hình.
import numpy as np

def smoothness_constrained_inversion(J, g, q_k, C_x, C_z, alpha_x=1.0, alpha_z=1.0, damping=0.1):
    """
    Thuật toán lặp nghịch đảo Gauss-Newton làm trơn 2D (Smoothness-Constrained Least-Squares)
    - J: Ma trận Jacobi (m x n)
    - g: Vector sai lệch log(rho_meas) - log(rho_calc) (m x 1)
    - q_k: Vector tham số mô hình hiện tại log(rho_model) (n x 1)
    - C_x, C_z: Ma trận vi phân bậc 1 theo phương x và z
    """
    # 1. Thiết lập ma trận làm trơn F
    F = alpha_x * (C_x.T @ C_x) + alpha_z * (C_z.T @ C_z)
    
    # 2. Xây dựng hệ phương trình chuẩn hóa
    JT_J = J.T @ J
    Hessian_mod = JT_J + damping * F
    Gradient = J.T @ g - damping * (F @ q_k)
    
    # 3. Giải hệ phương trình tìm bước nhảy tham số delta_q
    delta_q = np.linalg.solve(Hessian_mod, Gradient)
    
    # 4. Cập nhật vector mô hình mới
    q_next = q_k + delta_q
    return q_next

Testing và Validation

Quá trình nghịch đảo trên phần mềm Res2Dinv thực hiện qua 5 vòng lặp, độ hội tụ sai số bình phương trung bình (RMS) giảm dần rõ rệt:

  • Vòng lặp 1: $\text{RMS} = 14.2%$
  • Vòng lặp 3: $\text{RMS} = 6.1%$
  • Vòng lặp 5: $\text{RMS} = 3.82%$ (Đạt yêu cầu chuẩn hóa $< 5%$)

Đối chiếu kết quả ảnh điện Tuyến 3 và Hố khoan kiểm chứng

   Độ sâu (m)      Mặt cắt Điện trở suất 2D           Cột thạch học Hố khoan #3
    0.0 - 1.8     ρ = 120 - 450 Ω.m (Rất cao)       Cát hạt trung đến thô, khô rời
    1.8 - 5.5     ρ = 40 - 80 Ω.m (Trung bình)      Cát pha sét, ẩm ướt
    5.5 - 11.2    ρ = 2.5 - 18 Ω.m (Rất thấp)       Sét dẻo mềm, bùn bão hòa mặn
   11.2 - 15.0+   ρ = 180 - 650 Ω.m (Cao)           Sét lẫn sạn sỏi, đá gốc phong hóa
Tầng địa chất Chiều sâu theo ERT 2D (m) Chiều sâu theo Hố khoan (m) Điện trở suất ($\Omega\cdot\text{m}$) Mô tả địa chất - thạch học thực tế
Lớp 1 $0.0 - 1.8$ $0.0 - 1.75$ $120 - 450$ Cát hạt vừa, màu xám vàng, trạng thái khô xốp đến ẩm nhẹ.
Lớp 2 $1.8 - 5.5$ $1.75 - 5.30$ $40 - 80$ Cát pha sét dẻo cứng, chứa ít mùn thực vật, độ ẩm cao.
Lớp 3 $5.5 - 11.2$ $5.30 - 11.40$ $2.5 - 18$ Tầng sét yếu - bùn sét, xám đen, bão hòa nước mặn/lợ ven sông Đô Tỏa.
Lớp 4 $> 11.2$ $> 11.40$ $180 - 650$ Sạn sỏi gắn kết sét đặc chặt chuyển tiếp sang đá phong hóa nứt nẻ.

Sai lệch ranh giới giữa ảnh điện 2D và khoan thăm dò chỉ dao động trong khoảng $0.05 - 0.20\text{ m}$ ($< 3.5%$), chứng minh tính chuẩn xác vượt trội của mô hình điện trở suất nghịch đảo.

Đổi mới và đóng góp

  1. Ứng dụng tích hợp ERT 2D/3D trong khảo sát móng ven sông: Khắc phục triệt để hiện tượng dị thường giả của phương pháp 1D truyền thống và giải quyết vấn đề "điểm mù" giữa các hố khoan trong môi trường ngập nước/ven sông có độ dẫn điện biến thiên mạnh.
  2. Tối ưu hóa độ phân giải cấu hình Wenner-Alpha: Sử dụng hàm độ nhạy $F_{1D}(z) = \frac{2z}{(4z^2 + a^2)^{3/2}}$ với cực đại nhạy cảm tại $z \approx 0.35a$ và độ sâu trung bình $z_e = 0.519a$, giúp phân tách chi tiết các thấu kính sét mỏng mà phương pháp Dipole-Dipole dễ bị nhiễu biên tại vùng ven biển.
  3. Mô hình hóa không gian 3D giả lập: Ghép nối các lát cắt 2D thành khối 3D trên hệ trục $OXYZ$, cho phép kỹ sư xây dựng quan sát trực quan hướng dịch chuyển của tầng bùn sét yếu đối với an toàn bờ kè sông Đô Tỏa.

Ứng dụng thực tế và triển khai

Tình huống ứng dụng thực tế

  • Khảo sát thi công hố móng sâu và cọc khoan nhồi: Giúp nhà thầu định vị chính xác cao độ đặt mũi cọc chịu lực vào tầng cát sỏi đặc chắc/đá gốc (lớp 4, $\rho > 200\ \Omega\cdot\text{m}$), tránh đặt nhầm vào tầng bùn sét chịu tải kém (lớp 3).
  • Cảnh báo và kiểm soát xâm nhập mặn: Do nước ngầm nhiễm mặn có điện trở suất cực thấp ($< 5\ \Omega\cdot\text{m}$ theo định luật Archie $\rho = a \rho_w \Phi^{-m}$), phương pháp ERT lập tức khoanh vùng được phạm vi ngấm mặn từ sông Đô Tỏa vào đất nền công trình.

Phân tích hiệu quả kinh tế (ROI)

Hạng mục chi phí Phương án chỉ Khoan thăm dò (6 hố sâu 20m) Phương án ERT 2D (3 tuyến) + 1 Hố khoan kiểm chứng
Chi phí khoan ($600.000\text{ đ/m}$) $72.000.000\text{ VNĐ}$ ($120\text{ m}$) $12.000.000\text{ VNĐ}$ ($20\text{ m}$)
Chi phí đo đạc & Xử lý địa vật lý $0\text{ VNĐ}$ $15.000.000\text{ VNĐ}$
Tổng chi phí $72.000.000\text{ VNĐ}$ $27.000.000\text{ VNĐ}$
Mức tiết kiệm / Tỷ lệ tối ưu Baseline Tiết kiệm $62.5%$ chi phí, giảm $70%$ thời gian

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tính tương đương địa điện (Principle of Equivalence): Các lớp mỏng có tỷ số $h/\rho$ (đối với lớp dẫn điện) hoặc $T = h\cdot\rho$ (đối với lớp cách điện) không đổi có thể cho đáp ứng đường cong điện trở suất biểu kiến tương đương nhau trên mặt đất.
  • Nhiễu điện từ đô thị: Sự hiện diện của ống kim loại dẫn nước ngầm, cáp ngầm bọc kim loại hoặc dòng điện rò trong khu dân cư có thể làm cục bộ méo mó đường sức điện trường.

Hướng phát triển mở rộng

  • Ứng dụng phương pháp ảnh điện đo sâu phân cực kích thích (Induced Polarization - IP) 2D để phân biệt rõ ràng giữa thấu kính sét bão hòa nước ngọt và cát bão hòa nước mặn (vốn có điện trở suất tương đồng).
  • Tích hợp công nghệ Rada xuyên đất (Ground Penetrating Radar - GPR) ở tầng nông ($0 - 3\text{ m}$) để định vị chính xác công trình hạ tầng kỹ thuật trước khi đóng điện cực ERT.

Đối tượng hưởng lợi

  • Sinh viên & Học viên chuyên ngành Vật lý Địa cầu / Kỹ thuật Địa chất: Nắm vững trọn vẹn chu trình từ lý thuyết phương trình Laplace/Bessel, đạo hàm Frechet đến xử lý phần mềm Res2Dinv.
  • Kỹ sư Địa kỹ thuật & Trắc địa công trình: Sở hữu tài liệu kỹ thuật chuẩn xác về địa tầng khu vực Ngũ Hành Sơn phục vụ thiết kế biện pháp thi công tầng hầm, cừ larsen và hạ mực nước ngầm.
  • Chủ đầu tư & Doanh nghiệp xây dựng: Tối ưu hóa trên $60%$ ngân sách khảo sát địa chất giai đoạn tiền khả thi mà vẫn kiểm soát $100%$ rủi ro sụt lún địa tầng.
  • Các nhà nghiên cứu môi trường: Ứng dụng bản đồ phân bố điện trở suất để giám sát ô nhiễm nước dưới đất và ranh giới xâm nhập mặn các lưu vực sông ven biển miền Trung.

Câu hỏi thường gặp

1. Cần điều kiện mặt bằng thế nào để triển khai đo cắt lớp ảnh điện 2D Wenner-Alpha?

Mặt bằng cần có tuyến rải dây thẳng với chiều dài tối thiểu gấp $3 - 4$ lần chiều sâu cần khảo sát (ví dụ khảo sát sâu $15\text{ m}$ cần tuyến đo dài tối thiểu $45 - 60\text{ m}$). Đất tại các điểm cắm điện cực cần được làm ẩm bằng nước muối loãng nếu bề mặt quá khô để đảm bảo điện trở tiếp xúc $< 2\text{ k}\Omega$.

2. Tại sao đề tài lựa chọn cấu hình Wenner-Alpha thay vì Dipole-Dipole hay Pole-Pole?

Cấu hình Wenner-Alpha có tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR) rất cao và độ nhạy cực tốt đối với các biến thiên địa chất theo phương thẳng đứng (phân lớp địa tầng ngang). Trong môi trường đô thị nhiều nguồn nhiễu điện từ và địa chất ven sông bão hòa nước, Wenner-Alpha cho dữ liệu ổn định và ít điểm lỗi hơn Dipole-Dipole.

3. Phương pháp ERT có thể thay thế hoàn toàn khoan địa chất công trình không?

Không. ERT là phương pháp gián tiếp cung cấp bức tranh liên tục về điện trở suất. Cần tối thiểu $1 - 2$ hố khoan mẫu chuẩn (Ground Truth) để giải mã chính xác giá trị điện trở suất tương ứng với loại đất đá thực tế, loại bỏ tính không duy nhất (Non-uniqueness) của bài toán ngược.

4. Phần mềm Res2Dinv xử lý hiện tượng địa hình gồ ghề như thế nào?

Res2Dinv tích hợp thuật toán biến dạng lưới sai phân hữu hạn (Distorted Finite-Difference Grid) hoặc phần tử hữu hạn (FEM). Khi nhập file số liệu kèm tọa độ cao độ thực tế $(x, z)$ của từng điện cực, phần mềm sẽ tự động hiệu chỉnh độ cong của lưới tính toán thế điện theo mặt địa hình.

5. Chi phí đầu tư thiết bị đo ERT tự động và thời gian hoàn vốn (ROI)?

Một bộ máy đo ảnh điện đa cực $32 - 48$ kênh tiêu chuẩn có giá dao động từ $300 - 600\text{ triệu VNĐ}$. Với năng lực khảo sát $15 - 20\text{ dự án/năm}$, đơn vị khảo sát có thể thu hồi vốn đầu tư hoàn toàn trong vòng $10 - 14\text{ tháng}$ nhờ cắt giảm chi phí thuê máy khoan cồng kềnh.

Kết luận

Đề tài đã nghiên cứu và ứng dụng thành công phương pháp chụp cắt lớp ảnh điện 2D/3D kết hợp khoan thăm dò kiểm chứng tại khu vực ngã ba Bùi Tá Hán – Chương Dương, Quận Ngũ Hành Sơn, TP. Đà Nẵng. Kết quả nghịch đảo điện trở suất với sai số RMS ấn tượng $3.82%$ đã phân định sắc nét 4 tầng địa chất: lớp cát mặt rời ($120 - 450\ \Omega\cdot\text{m}$), lớp cát pha sét ẩm ($40 - 80\ \Omega\cdot\text{m}$), tầng bùn sét yếu bão hòa mặn ($2.5 - 18\ \Omega\cdot\text{m}$) và đới móng chịu lực đặc chắc ($> 180\ \Omega\cdot\text{m}$). Sự trùng khớp gần như tuyệt đối giữa ảnh điện và tài liệu khoan khẳng định ERT là phương pháp khảo sát không phá hủy hiện đại, tin cậy, tiết kiệm chi phí và mang lại giá trị thực tiễn cao cho ngành xây dựng hạ tầng đô thị ven biển.