Bối cảnh và vấn đề nghiên cứu

Trong lĩnh vực khảo cổ học và bảo tồn di sản, việc tìm kiếm và nhận dạng các di tích, công trình cổ bị chôn vùi dưới lòng đất đóng vai trò quan trọng trong việc lập quy hoạch, bảo tồn và phát huy giá trị văn hóa lịch sử. Trên thế giới, các phương pháp địa vật lý đã được ứng dụng trong khảo cổ từ những năm 1950 nhằm cung cấp bức tranh khái quát về sự phân bố của các đối tượng ngầm mà không cần tiến hành đào bới hay khai quật diện rộng làm xáo trộn hiện trường. Ngược lại, tại Việt Nam, việc triển khai các kỹ thuật địa vật lý phục vụ khảo cổ học còn rất khiêm tốn, chủ yếu mang tính tự phát hoặc hợp tác cá nhân nhỏ lẻ (như một số khảo sát bằng phương pháp từ và điện trở tại khu di tích Mỹ Sơn). Sự liên kết giữa các nhà nghiên cứu khảo cổ học và các chuyên gia địa vật lý vẫn còn hạn chế.

Trong số các phương pháp địa vật lý, phương pháp thăm dò điện giữ vai trò chủ đạo nhờ khả năng phân giải cấu trúc địa tầng và phân biệt ranh giới dị thường vật chất ở tầng nông. Kỹ thuật chụp cắt lớp ảnh điện hai chiều (2D) và ba chiều (3D) cho phép đo đạc đa cực tự động, xử lý số liệu nhanh chóng qua các phần mềm chuyên dụng và tái tạo mặt cắt phân bố điện trở suất thực. Đề tài khóa luận tốt nghiệp "Tìm kiếm công trình cổ bị chôn vùi bằng phương pháp chụp cắt lớp ảnh điện 2D-3D (có kiểm chứng bằng khoan thăm dò) tại Đà Lạt – Lâm Đồng" do sinh viên Trần Thị Hiếu thực hiện tại Khoa Vật lý, Trường Đại học Sư phạm – Đại học Đà Nẵng (năm 2018) nhằm giải quyết bài toán ứng dụng kỹ thuật ảnh điện vào thực tiễn khảo sát tầng nông.

Mục tiêu nghiên cứu của đề tài gồm 02 nội dung:

  1. Tìm kiếm dấu tích công trình cổ bị chôn vùi tại khu vực Đà Lạt – Lâm Đồng.
  2. Đánh giá độ chính xác và tính hiệu quả của phương pháp ảnh điện 2D trong thăm dò vật chất tầng nông thông qua việc khoan thăm dò kiểm chứng.

Nhiệm vụ nghiên cứu được xác định cụ thể theo 05 bước:

  1. Trình bày tổng quan cơ sở lý thuyết địa chất – vật lý của phương pháp thăm dò điện.
  2. Trình bày cơ sở lý thuyết của phương pháp ảnh điện hai chiều (2D), bao gồm bài toán thuận và bài toán ngược.
  3. Đánh giá độ nhạy của hệ thiết bị Wenner-Alpha để xây dựng bảng thiết bị đo phù hợp với đối tượng khảo sát.
  4. Nghiên cứu quy trình đo đạc thực nghiệm, thu thập, xử lý số liệu và giải đoán kết quả địa điện.
  5. Tiến hành khoan thăm dò kiểm chứng kết quả dị thường trên thực địa.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Sự truyền dẫn dòng điện, các thông số điện trở suất và cấu trúc phân bố địa chất tầng nông chứa dị thường công trình ngầm tại khu vực khảo sát.
  • Phạm vi không gian: Khu vực thực địa tại đèo Prenn, thành phố Đà Lạt, tỉnh Lâm Đồng.
  • Phạm vi thời gian: Quá trình nghiên cứu và khảo sát thực hiện từ tháng 12/2017 đến tháng 02/2018.

Cơ sở lý thuyết và phương pháp

Khóa luận thiết lập nền tảng dựa trên lý thuyết trường điện dừng trong môi trường bán không gian nhiều lớp và lý thuyết nghịch đảo địa vật lý. Dòng điện nhân tạo được phát vào lòng đất qua hệ điện cực nguồn, tạo nên trường điện thế đo được trên mặt đất. Các quy luật vật lý và mô hình toán học chính được sử dụng gồm:

  • Mô hình môi trường bán không gian phân lớp ngang đồng nhất bất đẳng hướng: Mỗi lớp thứ $i$ được đặc trưng bởi bề dày $h_i$, điện trở suất ngang $\rho_{ti}$, điện trở suất dọc $\rho_{ni}$, điện trở suất trung bình nhân $\rho_i = \sqrt{\rho_{ti} \rho_{ni}}$ và hệ số bất đẳng hướng $\lambda_i = \sqrt{\rho_{ni}/\rho_{ti}}$.
  • Phương trình truyền dẫn điện: Xuất phát từ phương trình liên tục $\text{div}\mathbf{J} = 0$, định luật Ohm vi phân $\mathbf{J} = \sigma \mathbf{E} = -\frac{1}{\rho}\text{grad}U$ và chuyển đổi sang hệ tọa độ trụ dẫn đến phương trình vi phân Bessel cấp 0. Hàm thế $U(r)$ trên bề mặt lớp thứ nhất được biểu diễn qua tích phân Fourier-Bessel với hàm nhân $R(\lambda)$ thiết lập theo công thức truy hồi Lipschitz – Lipskaia – Koefoed:

$$U(r) = \frac{I \rho_1}{2\pi} \int_0^\infty R(\lambda) J_0(\lambda r) d\lambda$$

  • Định luật thực nghiệm Archie: Mô tả mối quan hệ giữa điện trở suất của đá xốp $\rho$, điện trở suất của dung dịch chứa trong lỗ rỗng $\rho_w$ và độ rỗng bão hòa $\Phi$ qua biểu thức $\rho = b \cdot \rho_w \cdot \Phi^{-d}$ (với $b \approx 1, d \approx 2$).
  • Bài toán thuận 2D (Forward Modeling): Thiết lập phương trình phân bố thế $\nabla \cdot (\sigma \nabla U) = -I \delta(x-x_s)\delta(y-y_s)\delta(z-z_s)$, giải bằng phương pháp sai phân hữu hạn (Finite Difference Method) hoặc phần tử hữu hạn (Finite Element Method) chia lưới chữ nhật trên phần mềm Res2Dmod.
  • Bài toán ngược 2D (Inversion): Tìm phân bố điện trở suất thực từ các giá trị điện trở suất biểu kiến $\rho_a$ bằng thuật toán bình phương tối thiểu phi tuyến Gauss-Newton cải biên theo phương pháp Marquardt-Levenberg:

$$(\mathbf{J}^T \mathbf{J} + \lambda \mathbf{I})\Delta \mathbf{q} = \mathbf{J}^T \mathbf{g}$$

và phương pháp bình phương tối thiểu cưỡng bức đều đặn (Smoothness-constrained Least-squares) kết hợp tái gia trọng (Robust Inversion) nhằm làm trơn mô hình hoặc phân tách ranh giới khối dị thường sắc nét trên phần mềm Res2Dinv.

Phương pháp nghiên cứu thực tế gồm hai hợp phần:

  1. Phương pháp lý thuyết: Tổng quan cơ sở vật lý địa chất, thiết lập bài toán thế điện trường, phân tích hàm độ nhạy Fréchet 1D, 2D, 3D của các cấu hình điện cực (Wenner-Alpha, Pole-Pole, Dipole-Dipole).
  2. Phương pháp thực nghiệm: Đo đạc thực địa bằng hệ máy thăm dò điện đa cực theo cấu hình Wenner-Alpha dọc theo tuyến đo tại đèo Prenn, giải đoán mặt cắt 2D qua Res2Dinv và tiến hành khoan lấy mẫu lõi địa chất kiểm chứng.
Phương pháp / Kỹ thuật Mục đích sử dụng Công cụ / Phần mềm
Mô phỏng bài toán thuận 2D Tính toán điện trở suất biểu kiến lý thuyết từ mô hình phân lớp lưới chữ nhật Res2Dmod
Nghịch đảo bài toán ngược 2D Tối ưu hóa sai số $\mathbf{g} = \mathbf{y} - \mathbf{f}$, tái tạo mặt cắt điện trở suất thực Res2Dinv
Cấu hình Wenner-Alpha Thu thập số liệu điện thế $\Delta U$ với độ phân giải thẳng đứng cao Hệ máy đo địa điện đa cực
Khoan thăm dò địa kỹ thuật Xác minh cấu trúc địa tầng và đối tượng dị thường tại các vị trí nghi vấn Thiết bị khoan lấy mẫu lõi

Nội dung chính theo từng chương

Chương 1: Cơ sở vật lí - địa chất của phương pháp thăm dò điện

Chương 1 hệ thống hóa bản chất vật lý của sự truyền dẫn dòng điện trong môi trường đất đá và các yếu tố thạch học ảnh hưởng đến tham số điện trở suất.

Vật chất dưới mặt đất được chia thành hai nhóm dẫn điện chính: dẫn điện điện tử (kim loại tự nhiên như Pt, Au, Ag, Cu; khoáng vật sulfua như pirit, galenit, calcopirit; graphit) và dẫn điện điện phân/ion (tất cả nham thạch trầm tích, biến chất, phún xuất và các loại nước tự nhiên). Trong khảo sát tầng nông, cơ chế dẫn điện điện phân chiếm ưu thế tuyệt đối do đất đá tự nhiên luôn chứa dung dịch lỗ rỗng.

Nhóm khoáng vật / Đất đá Dải điện trở suất ($\Omega\cdot\text{m}$)
Vàng, bạch kim, bạc tự nhiên $< 10^{-5}$
Các sunfua (pirit, calcopirit, galenit...), canxiterit, than, grafic $10^{-5} - 1$
Hêmatit, bôcxit, kinôvar, anhydrit, selit $1 - 10^5$
Tràng thạch, thạch anh, calxit, mica, đá granit, đá bazan, dầu $10^5 - 10^{22}$

Điện trở suất của đất đá bị chi phối bởi các yếu tố địa chất và môi trường:

  • Độ ẩm và độ ngậm nước: Đất sét có độ ngậm nước cao nhất, dẫn đến độ dẫn điện lớn ($0{,}01 - 1\ \Omega^{-1}\cdot\text{m}^{-1}$), làm giảm điện trở suất xuống $20 - 2\ \Omega\cdot\text{m}$ so với đá phiến sạch ($6\cdot 10^7\ \Omega\cdot\text{m}$).
  • Độ rỗng và độ nứt nẻ: Mối quan hệ giữa độ rỗng $n_0$ và hệ số rỗng $e_0$ được biểu diễn qua công thức $n_0 = e_0 / (1 + e_0)$. Đất đá nứt nẻ bão hòa nước ngầm có điện trở suất giảm mạnh.
  • Độ khoáng hóa của nước ngầm: Nước ngầm đóng vai trò là chất điện phân mang ion ($\text{Na}^+, \text{K}^+, \text{Cl}^-$). Mối quan hệ giữa độ dẫn điện và nồng độ muối đồng biến ở mức dưới $100.000\ \text{mol/l}$ và chuyển sang nghịch biến khi đạt nồng độ bão hòa cao do các ion bị suy giảm độ linh động. Điện trở suất của nước khoáng quy đổi theo $\text{NaCl}$ có công thức thực nghiệm $\rho_w \approx 8{,}4 / M$ (với $M$ là độ khoáng hóa tính theo $\text{g/l}$).
  • Nhiệt độ: Điện trở suất tỷ lệ nghịch với nhiệt độ; khi nhiệt độ tăng thêm $10^\circ\text{C}$, độ dẫn điện của nước tăng $2 - 3%$, tuân theo quy luật $\rho_t = \rho_{18} / [1 + \alpha(t - 18)]$.

Chương 2: Tổng quan lý thuyết của phương pháp thăm dò điện và ảnh điện hai chiều (2D)

Chương 2 trình bày chi tiết toán học giải bài toán trường điện và quy trình xử lý dữ liệu ảnh điện 2D.

Từ sự phân bố điện thế của cặp cực dòng $C_1, C_2$ và cặp cực thế $P_1, P_2$, hiệu điện thế $\Delta U$ được tính theo biểu thức hình học:

$$\Delta U = \frac{I\rho}{2\pi} \left( \frac{1}{r_{C1P1}} - \frac{1}{r_{C2P1}} - \frac{1}{r_{C1P2}} + \frac{1}{r_{C2P2}} \right)$$

Điện trở suất biểu kiến được xác định qua $\rho_a = k \cdot R = k \cdot (\Delta U / I)$, trong đó $k$ là tham số hình học phụ thuộc vào cấu hình sắp xếp các điện cực. Đối với hệ thiết bị đối xứng Wenner với khoảng cách giữa các điện cực liên tiếp là $a$, công thức rút gọn là $\rho_{app.W} = 2\pi a (\Delta U / I)$.

Khóa luận phân tích sự chuyển tiếp từ đo sâu điện 1D (chỉ phản ánh biến đổi theo chiều sâu) và mặt cắt điện (chỉ phản ánh biến đổi theo phương ngang ở độ sâu cố định) sang phương pháp ảnh điện 2D. Trong bài toán thuận 2D, môi trường được rời rạc hóa thành lưới các khối chữ nhật có kích thước xác định. Bài toán ngược sử dụng ma trận Jacobian $\mathbf{J}$ chứa các đạo hàm riêng phần của đáp ứng mô hình theo từng tham số khối:

$$J_{ij} = \frac{\partial f_i}{\partial q_j}$$

Để giải ma trận Jacobian mà không cần tính lại toàn bộ bài toán thuận ở mỗi bước lặp, tác giả dẫn giải phương pháp tính đạo hàm Fréchet 3D ($F_{3D}$) và 2D ($F_{2D}$) cho cấu hình Pole-Pole theo công thức Loke và Barker (1995), sau đó tích phân số bằng phép cầu phương Gauss (Gauss Quadrature) với số điểm đánh giá hàm $n_x \times n_z$ từ $10 \times 8$ (khi khối ở gần điện cực) giảm dần xuống $4 \times 3$ (khi khối ở xa điện cực). Quá trình tối ưu hóa bình phương tối thiểu có ràng buộc điều kiện làm trơn ma trận vi phân bậc nhất giúp triệt tiêu các dao động kỳ dị của nghiệm.

Chương 3: Độ nhạy và quy trình đo của thiết bị Wenner-Alpha

Chương 3 tập trung vào đánh giá độ nhạy của hệ thiết bị Wenner-Alpha, quy trình triển khai ngoài thực địa tại Đà Lạt và đối sánh với kết quả khoan.

Độ sâu khảo sát đặc trưng được xác định qua hàm độ nhạy 1D $F_{1D}(z)$ theo mô hình Roy và Apparao (1971): hàm đạt cực đại tại độ sâu $z \approx 0{,}35a$ và tiệm cận về 0 khi $z \to \infty$. Chiều sâu khảo sát trung bình ($z_e$) của thiết bị Wenner-Alpha đạt giá trị $z_e = 0{,}519a$ (tương đương $z_e/L \approx 0{,}173$ với $L = 3a$ là chiều dài thiết bị). So với các cấu hình khác như Pole-Pole hay Dipole-Dipole, diện tích dưới đường cong độ nhạy của Wenner-Alpha xung quanh điểm cực đại tập trung hẹp hơn, chứng minh cấu hình này có độ phân giải theo phương thẳng đứng vượt trội và tỷ số tín hiệu trên nhiễu cao.

Quy trình khảo sát thực địa được tiến hành tại khu vực địa chất đèo Prenn (Đà Lạt – Lâm Đồng):

  • Bố trí tuyến đo điện đa cực thẳng hàng với khoảng cách điện cực cơ sở $a$ cố định, sử dụng hệ thiết bị Wenner-Alpha đo tuần tự các mức độ sâu $n = 1, 2, 3...$
  • Thu thập tệp số liệu điện thế và dòng phát, nạp vào phần mềm Res2Dinv để tiến hành lọc nhiễu, nghịch đảo và xuất mặt cắt ảnh điện 2D biểu diễn trường phân bố điện trở suất thực.
  • Dựa trên dị thường điện trở suất trên mặt cắt Tuyến 1, xác định các vị trí ranh giới công trình cổ ngầm và tiến hành khoan kiểm chứng tại hai tọa độ mốc: vị trí khoảng cách $12\ \text{m}$ và vị trí khoảng cách $35\ \text{m}$ trên tuyến đo.

Kết quả và đóng góp

Khóa luận đạt được các kết quả nghiên cứu cụ thể dựa trên số liệu thực nghiệm và tính toán lý thuyết:

  1. Thiết lập cơ sở tính toán độ nhạy cho cấu hình đo tầng nông: Xác lập bảng thông số chiều sâu khảo sát trung bình $z_e = 0{,}519a$ của hệ Wenner-Alpha, cung cấp cơ sở định lượng để lựa chọn khoảng cách điện cực $a$ và chiều dài tuyến đo $L$ phù hợp với độ sâu mục tiêu của các công trình ngầm tầng nông.
  2. Xây dựng quy trình xử lý dữ liệu ảnh điện 2D hoàn chỉnh: Ứng dụng thành công chuỗi phần mềm Res2Dmod và Res2Dinv trong việc mô phỏng và nghịch đảo số liệu địa điện thực tế bằng thuật toán bình phương tối thiểu làm trơn Gauss-Newton, kiểm soát sai số tính toán đạo hàm riêng phần dưới $5%$ so với nghiệm giải tích.
  3. Phát hiện dị thường địa điện và kiểm chứng bằng khoan thực địa:
    • Trên tuyến khảo sát tại khu vực đèo Prenn (Đà Lạt – Lâm Đồng), mặt cắt ảnh điện 2D đã bóc tách rõ sự phân tầng địa chất và khoanh vùng các dị thường điện trở suất cục bộ nghi vấn là tàn tích công trình cổ.
    • Kết quả khoan thăm dò cơ học tại vị trí $12\ \text{m}$ và vị trí $35\ \text{m}$ trên Tuyến 1 đã xác nhận sự hiện diện của các lớp vật chất và ranh giới cấu trúc địa chất ngầm trùng khớp với mô hình nghịch đảo ảnh điện 2D.
  4. Đóng góp về mặt phương pháp: Khóa luận chứng minh tính khả thi, độ tin cậy và sự vượt trội của phương pháp chụp cắt lớp ảnh điện 2D kết hợp kiểm chứng khoan trong điều kiện địa chất đồi núi phức tạp tại Tây Nguyên, tạo tiền đề phương pháp luận phục vụ công tác khảo cổ học không phá hủy tại Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

Dựa trên phạm vi và nội dung trình bày trong văn bản khóa luận, các hạn chế và hướng phát triển được xác định:

  • Hạn chế: Khảo sát thực tế trong khóa luận chủ yếu tập trung vào mặt cắt ảnh điện hai chiều (2D) dọc theo tuyến đơn, chưa triển khai toàn diện mạng lưới đo 3D không gian đầy đủ do hạn chế về thời gian khảo sát thực địa (tháng 12/2017 – 02/2018) và điều kiện thiết bị. Mặt cắt 2D giả định cấu trúc địa chất đồng nhất theo phương vuông góc với tuyến đo, điều này có thể dẫn đến sai số biên đối với các công trình cổ có dạng hình học 3D bất đối xứng.
  • Hướng nghiên cứu tiếp: Mở rộng áp dụng phương pháp ảnh điện 3D hoàn chỉnh trên diện tích khảo sát rộng hơn; kết hợp đồng thời phương pháp ảnh điện với các phương pháp địa vật lý tầng nông khác (như radar xuyên đất - GPR, thăm dò từ trường) nhằm tăng cường độ chính xác giải đoán đối tượng khảo cổ học.

Giá trị tham khảo

Tài liệu khóa luận có giá trị tham khảo chuyên môn đối với các nhóm đối tượng:

  • Sinh viên và học viên ngành Vật lý học, Địa vật lý, Địa kỹ thuật và Địa chất công trình: Nguồn tài liệu tham khảo chi tiết về cơ sở toán - lý của phương pháp điện trở trường dừng, phương pháp tính đạo hàm Fréchet và quy trình lập trình/sử dụng phần mềm mô phỏng Res2Dmod, Res2Dinv.
  • Cán bộ nghiên cứu khảo cổ học và bảo tồn di sản: Cung cấp tài liệu tham khảo thực tế về quy trình khảo sát địa vật lý không phá hủy, cách thức đọc hiểu mặt cắt dị thường điện trở suất 2D và phương pháp định vị điểm khoan kiểm chứng.
  • Phần đáng tham khảo nhất: Mục 2.2 về dẫn giải toán học của bài toán thuận - nghịch 2D và Mục 3.1 về phân tích hàm độ nhạy $F_{1D}, F_{2D}$ của thiết bị Wenner-Alpha.

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

1. Khóa luận sử dụng cấu hình thiết bị đo điện nào trong khảo sát thực địa tại Đà Lạt và lý do lựa chọn là gì?
Tác giả lựa chọn cấu hình thiết bị Wenner-Alpha. Lý do là cấu hình này có đồ thị hàm độ nhạy 1D tập trung hẹp xung quanh vị trí cực đại, mang lại độ phân giải theo phương thẳng đứng cao hơn cấu hình Pole-Pole, tỷ số tín hiệu trên nhiễu tốt và chiều sâu khảo sát trung bình xác định rõ ràng ($z_e = 0{,}519a$).

2. Chiều sâu khảo sát trung bình ($z_e$) của thiết bị Wenner-Alpha được tính theo công thức nào trong văn bản?
Theo kết quả tính toán độ nhạy 1D (dẫn theo Edwards, 1977 và Barker, 1991), chiều sâu khảo sát trung bình của thiết bị Wenner-Alpha bằng $z_e = 0{,}519a$ (với $a$ là khoảng cách giữa hai điện cực liên tiếp), tương ứng với tỉ số $z_e/L \approx 0{,}173$ (với $L = 3a$ là tổng chiều dài thiết bị).

3. Phần mềm nào được sử dụng để giải bài toán thuận và bài toán ngược trong khóa luận?
Khóa luận sử dụng phần mềm Res2Dmod để giải bài toán thuận (tính toán điện trở suất biểu kiến lý thuyết bằng phương pháp sai phân hữu hạn hoặc phần tử hữu hạn trên lưới chữ nhật) và phần mềm Res2Dinv để giải bài toán ngược (tối ưu hóa bình phương tối thiểu Gauss-Newton làm trơn để tái tạo điện trở suất thực 2D).

4. Kết quả khoan thăm dò kiểm chứng được tiến hành tại những vị trí nào trên thực địa?
Khoan kiểm chứng cơ học được thực hiện trực tiếp trên Tuyến 1 tại hai vị trí mốc cụ thể: vị trí khoảng cách $12\ \text{m}$ và vị trí khoảng cách $35\ \text{m}$ tại khu vực đèo Prenn, thành phố Đà Lạt, tỉnh Lâm Đồng.

5. Theo văn bản, độ dẫn điện của đất đá phụ thuộc vào các yếu tố vật lý - địa chất chủ yếu nào?
Độ dẫn điện của đất đá phụ thuộc vào thành phần khoáng vật, kiến trúc bên trong (tính bất đẳng hướng $\lambda$), độ rỗng và độ nứt nẻ ($n_0$), độ ẩm, nhiệt độ/áp suất và đặc biệt là độ ngậm nước cùng độ khoáng hóa của nước ngầm (nồng độ ion muối hòa tan như $\text{NaCl}, \text{KCl}$).

Kết luận

Khóa luận tốt nghiệp của sinh viên Trần Thị Hiếu đã trình bày toàn diện cơ sở lý thuyết toán - lý của phương pháp thăm dò điện và kỹ thuật chụp cắt lớp ảnh điện 2D. Tác giả đã phân tích sâu sắc hàm độ nhạy của thiết bị Wenner-Alpha, làm chủ quy trình mô hình hóa bài toán thuận - nghịch trên phần mềm chuyên dụng và áp dụng thành công vào khảo sát thực địa tại khu vực đèo Prenn (Đà Lạt – Lâm Đồng). Kết quả khoan thăm dò tại vị trí $12\ \text{m}$ và $35\ \text{m}$ trên Tuyến 1 đã khẳng định tính chính xác của mô hình địa điện, chứng minh tính hiệu quả của phương pháp ảnh điện 2D trong việc tìm kiếm các công trình cổ bị chôn vùi ở tầng nông.