Giới thiệu dự án

Tốc độ đô thị hóa nhanh chóng tại các đô thị ven biển như thành phố Đà Nẵng thúc đẩy sự bùng nổ của các dự án hạ tầng ngầm và nhà cao tầng. Theo thống kê địa kỹ thuật đô thị, hơn 40% sự cố thi công hố móng sâu và sụt lún công trình bắt nguồn từ việc đánh giá sai lệch cấu trúc địa chất phức tạp, đặc biệt là các thấu kính sét yếu, hang karst ngầm hoặc đới ngập nước mặn ven sông. Công tác khảo sát địa chất công trình truyền thống trước đây chủ yếu phụ thuộc vào khoan thăm dò cơ học đơn lẻ—phương pháp có độ tin cậy điểm cao nhưng chi phí đắt đỏ và để lại "khoảng mù" thông tin giữa các hố khoan.

Vấn đề cốt lõi đặt ra là làm thế nào để xây dựng bản đồ địa chất liên tục 2D và 3D với độ phân giải cao, chi phí tối ưu và có cơ chế kiểm chứng thực tế nhằm phục vụ thiết kế móng công trình sâu. Dự án tập trung giải quyết bài toán này thông qua phương pháp Chụp cắt lớp ảnh điện (Electrical Resistivity Tomography - ERT) đa cực kết hợp đối soát với dữ liệu khoan mẫu tại khu vực nút giao đường Bùi Tá Hán và Chương Dương (dọc sông Đô Tỏa), Quận Ngũ Hành Sơn, TP. Đà Nẵng.

Mục tiêu cụ thể của dự án:

  1. Xây dựng cơ sở lý thuyết truyền trường điện trong môi trường bán không gian bất đẳng hướng nhiều lớp dựa trên phương trình Bessel và hàm thế biến đổi Stefanescu.
  2. Tối ưu hóa cấu hình đo Wenner-Alpha dựa trên tính toán độ nhạy Fréchet 2D/3D và độ sâu khảo sát hiệu dụng ($Z_e \approx 0.519a$).
  3. Thu thập và xử lý số liệu điện trở suất biểu kiến ($\rho_a$) tại thực địa bằng phần mềm chuyên dụng Res2Dmod v3.0 và Res2Dinv v3.59.
  4. Tái tạo cấu trúc phân tầng địa chất 2D-3D và kiểm chứng định lượng với kết quả khoan địa chất thực tế nhằm đạt độ chính xác tương thích trên 90%.

Phạm vi nghiên cứu tập trung vào đới trầm tích Đệ Tứ ven sông Đô Tỏa, độ sâu thăm dò mục tiêu từ $0\text{ m}$ đến $30\text{ m}$, giới hạn bởi các yếu tố can thiệp dòng điện công nghiệp và độ dẫn cao do ranh giới mặn - ngọt.

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

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

Khảo sát địa chất hiện tại dựa vào ba nhóm giải pháp chính với các ưu nhược điểm được tổng hợp trong bảng đối sánh kỹ thuật:

Phương pháp Độ phân giải không gian Độ sâu tiếp cận Chi phí & Thời gian Hạn chế cốt lõi
Khoan thăm dò cơ học (Borehole Drilling) Rất cao tại 1 điểm (1D) Rất sâu ($>100\text{ m}$) Rất cao, tiến độ chậm ($2-3\text{ ngày/hố}$) Bỏ sót bất đồng nhất cục bộ giữa các lỗ khoan.
Đo sâu điện trở truyền thống (1D VES) Trung bình theo phương đứng Sâu ($20-50\text{ m}$) Thấp, thiết bị gọn nhẹ Không phản ánh được biến thiên điện trở suất theo phương ngang.
Chụp cắt lớp ảnh điện 2D/3D (ERT - Đề xuất) Cao liên tục cả phương đứng & ngang Trung bình ($10-60\text{ m}$) Tối ưu, đo tự động nhanh ($1-2\text{ giờ/tuyến}$) Hiện tượng tương đương (Equivalence) khi giải bài toán ngược.

Yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Hệ thu nhận điện trở suất 4 cực đa kênh, thuật toán nghịch đảo bình phương tối thiểu có trọng số (Smoothness-constrained Least-squares), sai số bình phương trung bình (RMS error) $< 5%$.
  • Should have: Khả năng tính toán tự động ma trận Jacobian từ phép cầu phương Gauss (Gauss Quadrature), hiển thị mô hình khối 3D.
  • Could have: Tích hợp trực tiếp dữ liệu địa hình bề mặt bằng thuật toán lưới biến dạng FEM.
  • Won't have: Xử lý ảnh hưởng điện cực phân cực cảm ứng (IP) miền thời gian phức tạp trong giai đoạn 1.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu địa vật lý được xây dựng từ phương trình Poisson dẫn xuất từ định luật Ohm vi phân và định luật bảo toàn điện tích:

$$\nabla \cdot \vec{J} = -\nabla \cdot (\sigma \nabla U) = I \delta(x - x_s)\delta(y - y_s)\delta(z - z_s)$$

Trong đó:

  • $\sigma = \frac{1}{\rho}$ là độ dẫn điện của môi trường $(\Omega^{-1}\cdot\text{m}^{-1})$.
  • $U$ là điện thế tại điểm quan sát $(\text{V})$.
  • $I$ là cường độ dòng điện phát $(\text{A})$.
  • $\delta$ là hàm Dirac delta biểu diễn nguồn dòng điểm tại tọa độ $(x_s, y_s, z_s)$.

Công nghệ sử dụng:

  • Thiết bị phát/thu địa vật lý đa cực: Máy đo ảnh điện tự động đa kênh Syscal Pro (Iris Instruments) kết hợp cáp đa cực 24-48 kênh.
  • Engine mô phỏng bài toán thuận: Res2Dmod v3.0 (Finite Difference Method - FDM và Finite Element Method - FEM).
  • Engine giải bài toán ngược: Res2Dinv v3.59 (Inversion code C/Fortran) áp dụng tối ưu hóa Gauss-Newton hiệu chỉnh Marquardt-Levenberg.
  • Phần mềm đồ họa không gian: Surfer v13.0 và Voxler v4.0 phục vụ nội suy và mô hình hóa thể tích 3D.

Methodology

Phương pháp luận nghiên cứu áp dụng quy trình lặp khoa học thực nghiệm:

  1. Khảo sát tham số tiên nghiệm: Đo kiểm tra điện trở suất mẫu đá, độ khoáng hóa nước ngầm ($M$) thông qua công thức thực nghiệm: $$\rho_w \approx \frac{8.4}{M}\quad (\Omega\cdot\text{m})$$
  2. Thiết kế tuyến đo thực địa: Bố trí tuyến đo trực giao và song song bờ sông Đô Tỏa với khoảng cách điện cực ban đầu $a = 1.0\text{ m} - 2.5\text{ m}$, thừa số độ sâu $n = 1 \div 6$.
  3. Đánh giá rủi ro: Khắc phục điện trở tiếp xúc cao do bề mặt cát khô bằng cách tưới dung dịch điện phân loãng ($\text{NaCl}; 1%$) tại chân cọc đồng; xử lý hiện tượng rò rỉ dòng phát bằng cách đo kiểm tra thế phân cực tự nhiên (SP) trước mỗi chu kỳ phát xung đảo chiều.

Implementation và kết quả

Development process

Cốt lõi của quá trình giải bài toán ngược là thuật toán bình phương tối thiểu cưỡng bức độ trơn (Smoothness-Constrained Least-Squares Inversion):

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

Trong đó:

  • $J$ là ma trận Jacobi chứa các đạo hàm riêng Fréchet $J_{ij} = \frac{\partial f_i}{\partial q_j}$.
  • $g = y - f(q_k)$ là vector độ lệch giữa logarit điện trở suất đo được và giá trị lý thuyết.
  • $\lambda$ là tham số giảm chấn Marquardt (damping factor).
  • $F = \alpha_x C_x^T C_x + \alpha_z C_z^T C_z$ là toán tử làm trơn không gian bậc 1 theo trục $x$ và $z$.

Hàm độ nhạy 3D (Fréchet derivative) cho cấu hình 2 điện cực Pole-Pole:

$$F_{3D}(x, y, z) = \frac{1}{4\pi^2} \frac{x(x - a) + y^2 + z^2}{\left[(x^2 + y^2 + z^2)((x - a)^2 + y^2 + z^2)\right]^{3/2}}$$

Đạo hàm riêng 2D của một khối chữ nhật hữu hạn được tích hợp bằng phép cầu phương Gauss (Gauss Quadrature):

$$\frac{\partial U}{\partial \rho} = \int_{z_1}^{z_2} \int_{x_1}^{x_2} F_{2D}(x, z) dx dz \approx \sum_{l=1}^{n_z} \sum_{k=1}^{n_x} W_l W_k F_{2D}(u_k, v_l)$$

import numpy as np

def calculate_frechet_2d(x, z, a, n_gauss=8):
    """
    Tính hàm độ nhạy Fréchet 2D tích phân theo phương y (-inf đến +inf)
    sử dụng phép cầu phương Gauss-Legendre.
    """
    # Lấy các điểm nút và trọng số Gauss-Legendre
    y_nodes, y_weights = np.polynomial.legendre.leggauss(n_gauss)
    # Scale biến y từ (-inf, +inf) qua biến đổi tọa độ y = tan(theta)
    theta = y_nodes * (np.pi / 2)
    y = np.tan(theta)
    dy = (np.pi / 2) * (1 + y**2) * y_weights
    
    f2d_val = 0.0
    for yi, wy in zip(y, dy):
        r1_sq = x**2 + yi**2 + z**2
        r2_sq = (x - a)**2 + yi**2 + z**2
        numerator = x * (x - a) + yi**2 + z**2
        denominator = (r1_sq * r2_sq)**1.5
        f3d = (1.0 / (4.0 * np.pi**2)) * (numerator / denominator)
        f2d_val += f3d * wy
        
    return f2d_val

def wenner_geometric_factor(a):
    """Hệ số hình học k cho mảng cực Wenner-Alpha (C1 - a - P1 - a - P2 - a - C2)"""
    return 2.0 * np.pi * a

Testing và validation

Hiệu chuẩn mô hình số được thực hiện qua các bài kiểm thử kiểm soát sai số:

  • Độ hội tụ nghịch đảo: Sai số RMS giảm nhanh qua 5 bước lặp từ $14.2%$ xuống còn $2.18%$.
  • Độ nhạy độ sâu khảo sát ($Z_e$): Cấu hình Wenner-Alpha đạt đỉnh hàm độ nhạy 1D ở độ sâu $z = 0.35a$, với chiều sâu khảo sát trung bình hiệu dụng $Z_e = 0.519a$.
BẢNG DỮ LIỆU ĐỐI SOÁT ĐỊA CHẤT TẠI HỐ KHOAN TUYẾN 3:
-------------------------------------------------------------------------------
Độ sâu (m)    Thạch học thực tế (Mẫu khoan)     Điện trở suất ERT (Ω.m)   Trạng thái bão hòa
-------------------------------------------------------------------------------
0.0 - 1.8     Cát hạt mịn, xốp, màu vàng nhạt   180 - 450                 Khô, nằm trên mực nước ngầm
1.8 - 6.5     Cát hạt trung - thô ngậm nước     40 - 85                   Bão hòa nước ngọt
6.5 - 14.2    Sét dẻo mềm, lẫn mùn thực vật      8 - 22                    Độ dẫn cao, bão hòa ion
14.2 - 22.0   Cát pha sét, hạt mịn chặt vừa      25 - 55                   Nhiễm mặn nhẹ ven sông
> 22.0        Đá cát kết nứt nẻ nhẹ (bedrock)   550 - 1200                Tầng nền vững chắc
-------------------------------------------------------------------------------

Kết quả đạt được

  1. Phân tích địa tầng chính xác: Tái tạo thành công mặt cắt địa điện với chiều dài tuyến $120\text{ m}$, chiều sâu nghiên cứu đạt $25\text{ m}$.
  2. Khảo sát đới biến động điện trở: Định vị chính xác tầng sét yếu ở độ sâu từ $6.5\text{ m}$ đến $14.2\text{ m}$ với điện trở suất đặc trưng rất thấp ($\rho < 20,\Omega\cdot\text{m}$), giải thích bởi định luật Archie khi độ xốp $\Phi \approx 45%$ và tỷ lệ sét cao: $$\rho = b \cdot \rho_w \cdot \Phi^{-d}$$
  3. Đối chuẩn khoan mẫu: Ranh giới giữa lớp sét yếu và tầng cát ngậm nước ngầm đo bằng ERT trùng khớp với mẫu lõi khoan thăm dò với sai số độ sâu dưới $0.3\text{ m}$ ($< 4.5%$).

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

  • Tối ưu hóa ma trận Jacobian 2D: Sử dụng phương pháp giải tích kết hợp Gauss Quadrature giúp giảm $42%$ thời gian tính ma trận độ nhạy so với phương pháp FDM thuần túy trên cùng kích thước lưới phần tử.
  • Quy trình kết hợp 2D-3D kết hợp khoan chuẩn: Khắc phục triệt để tính chất bất khả định của bài toán ngược địa vật lý bằng cách đưa thông tin thạch học từ hố khoan vào làm ràng buộc tiên nghiệm (A priori constraint) trong ma trận làm trơn $R_m$.
  • Bản đồ hóa ranh giới nước ngầm và nhiễm mặn: Định lượng hóa mức độ xâm nhập mặn từ sông Đô Tỏa thông qua hệ số dẫn điện ion của môi trường lỗ rỗng, phân tách rõ vùng nước ngọt ($\rho > 40,\Omega\cdot\text{m}$) và vùng ngấm mặn ($\rho < 15,\Omega\cdot\text{m}$).

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

Kịch bản ứng dụng kỹ thuật

  1. Thiết kế hố móng sâu (Deep Excavation): Cung cấp mặt cắt liên tục giúp kỹ sư định vị chính xác vị trí cần hạ mực nước ngầm hoặc cắm cừ larsen chặn tầng cát chảy/bùn sét yếu.
  2. Khảo sát móng cọc nhà cao tầng: Xác định độ sâu tầng đá gốc (bedrock, $\rho > 500,\Omega\cdot\text{m}$) biến thiên phức tạp dọc công trình, tránh rủi ro chiều dài cọc không đồng đều.

Phân tích chi phí và hiệu quả đầu tư (Cost-Benefit Analysis)

So với phương án chỉ dùng khoan thăm dò truyền thống (mật độ $20\text{ m/hố}$), mô hình kết hợp ERT + Khoan kiểm chứng định điểm mang lại lợi ích định lượng:

  • Tiết kiệm chi phí khảo sát: Giảm $45 - 60%$ tổng kinh phí khảo sát nhờ giảm số lượng lỗ khoan không cần thiết từ 10 hố xuống còn 3 hố kiểm chứng.
  • Rút ngắn thời gian thực địa: Giảm $50%$ thời gian thu thập dữ liệu hiện trường (từ 15 ngày xuống 7 ngày cho diện tích $5.000\text{ m}^2$).
  • Hạn chế rủi ro thi công: Giảm thiểu nguy cơ gặp túi bùn ngầm gây nghiêng lệch hố móng trong giai đoạn đào hầm.

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

  • Hạn chế kỹ thuật: Hiện tượng suy giảm độ nhạy theo chiều sâu khiến độ phân giải biên của các tầng địa chất dưới $20\text{ m}$ giảm dần; nhiễu điện từ do mạng lưới dây điện và cọc tiếp địa ngầm trong đô thị gây biến dạng cục bộ đường cong đẳng thế.
  • Ràng buộc tài nguyên: Cần duy trì tiếp xúc điện cực tốt trên bề mặt bê tông/nhựa đường—đòi hỏi kỹ thuật khoan mồi hoặc dùng điện cực tấm phẳng điện dung.
  • Hướng phát triển: Tích hợp kỹ thuật đo điện trở suất 3D thời gian thực (Time-lapse ERT 4D) để quan trắc động lực học thấm nước ngầm xung quanh tường vây hố móng trong suốt quá trình đào đất.

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

  • Sinh viên & Học viên ngành Địa kỹ thuật/Vật lý Địa cầu: Cung cấp tài liệu mẫu mực về chuyển đổi bài toán trường điện Maxwell sang thuật toán số FDM/FEM và thực thi nghịch đảo dữ liệu thực tế.
  • Kỹ sư Địa chất công trình: Nắm vững giải pháp thay thế khoan đơn lẻ, sở hữu công cụ lập mô hình tầng sét/cát liên tục phục vụ tính toán lún móng.
  • Doanh nghiệp Xây dựng & Chủ đầu tư: Tiết kiệm hàng trăm triệu đồng chi phí khoan khảo sát, tối ưu hóa thiết kế cọc móng và ngăn ngừa sự cố sập thành vách tầng hầm.
  • Nhà nghiên cứu địa môi trường: Phương pháp luận chuẩn để ứng dụng quan trắc ranh giới mặn - ngọt tại các vùng ven biển chịu biến đổi khí hậu.

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

1. Cần chuẩn bị trang thiết bị và điều kiện mặt bằng như thế nào để triển khai đo ERT?

Mặt bằng cần có tuyến rải cáp thông thoáng với chiều dài tối thiểu gấp 4 đến 5 lần độ sâu cần khảo sát (do với cấu hình Wenner, độ sâu khảo sát $Z_{max} \approx 0.17L$ với $L$ là chiều dài tuyến). Thiết bị gồm máy đo điện trở chính, cáp đa cực, cọc điện cực thép không gỉ/đồng, ắc quy nguồn $12\text{V}-24\text{V}$ và nước muối để giảm điện trở tiếp xúc chân cọc ($R_{contact} < 2\text{ k}\Omega$).

2. Giới hạn độ sâu phân giải của phương pháp ảnh điện 2D là bao nhiêu?

Độ sâu phân giải phụ thuộc vào chiều dài tuyến đo ($L$), khoảng cách cực nhỏ nhất ($a$) và công suất nguồn phát. Với cấu hình Wenner-Alpha sử dụng cáp 48 điện cực khoảng cách $a = 2.5\text{ m}$ ($L = 117.5\text{ m}$), độ sâu hiệu dụng đạt khoảng $20 - 25\text{ m}$. Muốn sâu hơn ($>50\text{ m}$), cần chuyển đổi cấu hình Wenner-Schlumberger hoặc Pole-Dipole với khoảng cách cực lớn hơn.

3. Phương pháp xử lý ra sao khi gặp công trình có lớp bê tông mặt đường?

Khi đo trên mặt đường bê tông hoặc nhựa, cọc đồng thông thường không thể đóng trực tiếp. Cần sử dụng các giải pháp thay thế: (1) Khoan lỗ mồi nhỏ đường kính $8\text{ mm}$ qua lớp bê tông xuống đất nền; (2) Sử dụng điện cực tấm phẳng (flat-base electrodes) kết hợp gel dẫn điện chuyên dụng hoặc bentonite ướt; (3) Sử dụng phương pháp điện dung không tiếp xúc (Capacitively Coupled Resistivity - CCR).

4. Tại sao vẫn cần khoan thăm dò khi đã có ảnh điện 2D-3D chi tiết?

Phương pháp ảnh điện đo gián tiếp tham số điện trở suất ($\Omega\cdot\text{m}$), vốn chịu ảnh hưởng đồng thời bởi độ ẩm, độ khoáng hóa, thành phần khoáng vật và độ rỗng. Do nguyên lý tương đương (nhiều tổ hợp thạch học khác nhau có thể cho cùng một giá trị điện trở suất), khoan thăm dò đóng vai trò là "chân lý thực địa" (ground-truth) để gán chính xác ranh giới địa tầng và kiểm chuẩn mô hình nghịch đảo.

5. Chi phí đầu tư thiết bị và thời gian hoàn vốn (ROI) cho một đội khảo sát địa kỹ thuật?

Một hệ thống đo ERT đa cực tự động tiêu chuẩn có giá dao động từ $25.000 - $60.000 USD tùy số kênh. Đối với các đơn vị khảo sát địa kỹ thuật thực hiện từ 15 đến 20 dự án công trình cao tầng/năm, khả năng hoàn vốn đạt được trong vòng 12 - 18 tháng nhờ cắt giảm chi phí thuê máy khoan và nâng cao năng lực trúng thầu các dự án hạ tầng lớn.

Kết luận

Đồ án đã chứng minh sự vượt trội của phương pháp chụp cắt lớp ảnh điện 2D-3D trong việc làm sáng tỏ cấu trúc địa tầng ngầm tại khu vực ven sông Ngũ Hành Sơn, Đà Nẵng. Sự kết hợp nhuần nhuyễn giữa nền tảng vật lý lý thuyết (phương trình vi phân Bessel, biến đổi Stefanescu, tích phân đạo hàm Fréchet) và kỹ thuật tính toán hiện đại (Gauss-Newton damping, FDM/FEM) đã tạo nên một quy trình khảo sát địa vật lý chặt chẽ, đạt độ tin cậy cao khi đối chuẩn trực tiếp với mẫu khoan địa chất. Đây là giải pháp công nghệ then chốt, mang lại giá trị kinh tế và độ an toàn kỹ thuật cao, cần được nhân rộng trong công tác quy hoạch, thiết kế nền móng công trình đô thị ven biển hiện đại.