Luận Văn: Áp Dụng Phương Pháp Cộng Hưởng Từ Để Tìm Kiếm Nước Ngầm Ở Miền Bắc Việt Nam

Chuyên khảo phân tích Luận văn nghiên cứu áp dụng phương pháp cộng hưởng từ để tìm kiếm nước ngầm một số vùng trọng điểm, đánh giá các khía cạnh quan trọng, đề xuất hướng nghiên

• Tác giả: Ts. Tăng Đình Nam
• Trường học: Viện Khoa học Địa chất và Khoáng sản
• Chuyên ngành: Địa chất
• Đề tài: Nghiên cứu, áp dụng phương pháp cộng hưởng từ để tìm kiếm nước ngầm một số vùng trọng điểm ở miền Bắc Việt Nam
• Loại tài liệu: báo cáo tổng kết đề tài nghiên cứu khoa học công nghệ cấp bộ
• Năm xuất bản: 2009
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN VỀ ĐO SÂU CỘNG HƯỞNG TỪ
• 1.1. Nguyên lý cơ bản và bài toán thuận của đo sâu cộng hưởng từ
• 1.1.1. Nguyên lý cơ bản của đo sâu cộng hưởng từ
• 1.1.2. Bài toán thuận của đo sâu cộng hưởng từ
• 1.1.3. Môi trường 1D
• 1.2. Bài toán ngược cho môi trường điện trở suất cao, 1D, khung dây đơn
• 1.3. Hàm nhân trong nửa không gian dưới đất cách điện và dẫn điện
• 1.4. Đặc điểm tương đương của đường cong ĐSCHT 1D
• 1.5. Ba mô hình 2D
• 1.6. Tình hình áp dụng đo sâu cộng hưởng từ để khảo sát trực tiếp nước ngầm trên thế giới và ở Việt Nam
• 2. CHƯƠNG 2: THU THẬP, XỬ LÝ, NGHỊCH ĐẢO VÀ GIẢI THÍCH SỐ LIỆU ĐO SÂU CỘNG HƯỞNG TỪ
• 2.1. Thu thập số liệu và các biện pháp chống nhiễu khi đo
• 2.1.1. Đặc điểm nhiễu điện từ của đường dây tải điện
• 2.1.2. Khung dây số tám
• 2.1.3. Khung dây bù
• 2.1.4. Chống nhiễu đỉnh nhọn
• 2.2. Thu thập và đánh giá chất lượng số liệu ĐSCHT
• 2.2.1. Chọn khung dây
• 2.2.2. Chọn các tham số trong thu thập số liệu
• 2.2.3. Đánh giá chất lượng số liệu ĐSCHT
• 2.3. Xử lý, nghịch đảo và giải thích số liệu đo sâu cộng hưởng từ
• 2.3.1. Xử lý số liệu đo sâu cộng hưởng từ
• 2.3.2. Nghịch đảo số liệu đo sâu cộng hưởng từ
• 2.3.3. Giải thích số liệu đo sâu cộng hưởng từ
• 3. CHƯƠNG 3: KẾT QUẢ NGHIÊN CỨU THỬ NGHIỆM VÙNG THỊ XÃ LAI CHÂU
• 3.1. Đặc điểm địa chất thủy văn vùng thị xã Lai Châu
• 3.1.1. Cấu trúc, kiến tạo
• 3.1.2. Đặc điểm địa chất thủy văn
• 3.2. Các phương pháp và kỹ thuật thi công
• 3.2.1. Phương pháp và khối lượng đã thực hiện
• 3.2.2. Mạng lưới khảo sát
• 3.2.3. Kỹ thuật thi công các phương pháp địa vật lý
• 3.3. Ảnh hưởng của gradien trường từ trong diện tích nghiên cứu
• 3.4. Sự thay đổi của trường địa từ B0 trong không gian và từ tính của đất đá trong đo sâu cộng hưởng từ
• 3.5. Đo vẽ trường cảm ứng từ B0 tại phường Quyết Thắng, thị xã Lai Châu
• 3.6. Ảnh hưởng biến thiên trường địa từ theo chu kỳ 24 giờ và giông sét
• 3.6.1. Ảnh hưởng biến thiên trường địa từ theo chu kỳ 24 giờ
• 3.6.2. Ảnh hưởng của giông sét
• 3.7. Các biện pháp chống nhiễu điện từ
• 3.7.1. Chống nhiễu điện từ nhân tạo
• 3.7.2. Chống nhiễu điện từ tự nhiên
• 3.8. Kết quả đo sâu phân cực ở phường Quyết Thắng thị xã Lai Châu
• 3.8.1. Kết quả trên tuyến T3a-0
• 3.8.2. Kết quả trên tuyến T3a-3
• 3.9. Kết quả đo sâu cộng hưởng từ ở phường Quyết Thắng thị xã Lai Châu
• 3.9.1. Điểm đo tại cọc -25 tuyến T3a-0
• 3.9.2. Điểm đo tại cọc -25 tuyến T3a-1
• 3.9.3. Điểm đo tại cọc -51 tuyến T3a-3
• 3.10. Kết quả khoan tại khu phường Quyết thắng
• 3.11. Kết quả đo sâu phân cực ở khu vực xã Nậm Loỏng thị xã Lai Châu
• 3.12. Kết quả đo sâu cộng hưởng từ ở khu vực xã Nậm Loỏng
• 3.12.1. Điểm đo tại cọc 0 tuyến T1
• 3.12.2. Điểm đo tại cọc 5 tuyến Tb1
• 3.13. Kết quả khoan tại cọc 5 tuyến T1b khu vực xã Nậm Loỏng thị xã Lai Châu
• 3.14. Đánh giá hiệu quả của đo sâu cộng hưởng từ ở vùng thị xã Lai Châu
• 4. CHƯƠNG 4: KẾT QUẢ NGHIÊN CỨU THỬ NGHIỆM VÙNG BÓ MƯỜI, THUẬN CHÂU, SƠN LA
• 4.1. Đặc điểm địa chất thủy văn ở vùng Bó Mười, Thuận Châu
• 4.1.1. Cấu trúc kiến tạo
• 4.1.2. Đặc điểm địa chất thủy văn
• 4.2. Các phương pháp địa vật lý và kỹ thuật thi công
• 4.2.1. Phương pháp và khối lượng đã thực hiện
• 4.2.2. Mạng lưới khảo sát
• 4.2.3. Kỹ thuật thi công các phương pháp địa vật lý
• 4.3. Ảnh hưởng của gradien trường từ trong diện tích nghiên cứu
• 4.4. Ảnh hưởng biến thiên trường địa từ theo chu kỳ 24 giờ và giông sét
• 4.5. Kết quả đo sâu phân cực ở Bản Bó, Bó Mười, Thuận Châu
• 4.5.1. Kết quả trên tuyến T2
• 4.5.2. Kết quả trên tuyến T0
• 4.6. Kết quả ĐSCHT ở Bản Bó, Bó Mười, Thuận Châu
• 4.6.1. ĐSCHT S9 tại cọc -5 tuyến T2
• 4.6.2. ĐSCHT S10 và so sánh kết quả nghịch đảo của S9 và S10
• 4.6.3. ĐSCHT S17 tại cọc -20 tuyến T0
• 5. CHƯƠNG 5: ĐÁNH GIÁ KẾT QUẢ XỬ LÝ
• 5.1. ĐSCHT S5 tại cọc 5 tuyến T1b - Nậm Loỏng
• 5.2. ĐSCHT S18 tại Mộc Châu, Sơn La
• 5.3. ĐSCHT S17 tại Thuận Châu, Sơn La
• 5.4. Đánh giá kết quả nghịch đảo các điểm ĐSCH
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về đo sâu cộng hưởng từ

Đo sâu cộng hưởng từ (ĐSCHT) là phương pháp địa vật lý duy nhất trên thế giới được sử dụng để khảo sát trực tiếp nước ngầm từ mặt đất. Phương pháp này dựa trên hiện tượng cộng hưởng từ hạt nhân, đặc biệt là proton trong phân tử nước. Các phương pháp địa vật lý khác chỉ có thể khảo sát nước ngầm một cách gián tiếp thông qua các tham số vật lý như điện trở suất, hằng số điện môi, hoặc tốc độ truyền sóng địa chấn. ĐSCHT được phát triển từ những năm 1980 bởi các nhà khoa học Nga, với máy Hydroscope là thiết bị đầu tiên được sử dụng. Phương pháp này đã được áp dụng rộng rãi ở nhiều quốc gia như Ấn Độ, Đức, Mỹ, và Trung Quốc. Tại Việt Nam, ĐSCHT được nghiên cứu và áp dụng tại các vùng trọng điểm ở miền Bắc Việt Nam, đặc biệt là trong các khu vực đá vôi karst.

1.1 Nguyên lý cơ bản của đo sâu cộng hưởng từ

Nguyên lý cơ bản của ĐSCHT dựa trên hiện tượng cộng hưởng từ hạt nhân của proton trong phân tử nước. Proton có momen từ và spin, khi đặt trong trường từ tĩnh B0, các momen từ này sẽ định hướng theo phương của trường từ. Khi phát một xung điện từ với tần số Larmor, các momen từ sẽ lệch khỏi phương ban đầu và tạo ra tín hiệu cộng hưởng từ. Tín hiệu này được thu nhận và phân tích để xác định sự hiện diện và phân bố của nước ngầm. Phương pháp này cho phép xác định trực tiếp lượng nước tự do trong lòng đất, loại bỏ các yếu tố nhiễu từ nước liên kết hoặc chất rắn.

1.2 Bài toán thuận của đo sâu cộng hưởng từ

Bài toán thuận của ĐSCHT liên quan đến việc mô hình hóa tín hiệu cộng hưởng từ trong môi trường dẫn điện. Các nghiên cứu của Weichman và cộng sự (1999, 2000) đã đưa ra biểu thức tổng quát cho tín hiệu cộng hưởng từ trong môi trường 3D, tính đến sự phân cực elip của trường từ kích thích. Tín hiệu cộng hưởng từ phụ thuộc vào các tham số như độ từ khuynh, điện trở suất, và vị trí địa lý. Các kết quả nghiên cứu lý thuyết này đã mở đường cho việc áp dụng ĐSCHT trong các điều kiện địa chất phức tạp, đặc biệt là ở các vùng đá vôi karst.

II. Thu thập xử lý và giải thích số liệu đo sâu cộng hưởng từ

Quá trình thu thập và xử lý số liệu ĐSCHT đòi hỏi sự chính xác cao để đảm bảo kết quả khảo sát nước ngầm được tin cậy. Các bước chính bao gồm thu thập số liệu, chống nhiễu điện từ, xử lý và nghịch đảo số liệu. Nhiễu điện từ là một trong những thách thức lớn nhất khi áp dụng ĐSCHT, đặc biệt là trong các khu vực có mật độ dân cư cao hoặc gần các công trình công nghiệp. Các biện pháp chống nhiễu như sử dụng khung dây bù và lọc tín hiệu đã được áp dụng để cải thiện chất lượng số liệu. Sau khi thu thập, số liệu được xử lý và nghịch đảo để xác định phân bố nước ngầm trong lòng đất.

2.1 Thu thập số liệu và chống nhiễu

Thu thập số liệu ĐSCHT đòi hỏi sự chính xác trong việc lựa chọn khung dây và các tham số đo. Nhiễu điện từ, đặc biệt là từ các đường dây tải điện và giông sét, có thể ảnh hưởng nghiêm trọng đến chất lượng số liệu. Các biện pháp chống nhiễu như sử dụng khung dây bù và lọc tín hiệu đỉnh nhọn đã được áp dụng để giảm thiểu nhiễu. Việc đánh giá chất lượng số liệu ngay tại hiện trường là bước quan trọng để đảm bảo kết quả khảo sát chính xác.

2.2 Xử lý và nghịch đảo số liệu

Sau khi thu thập, số liệu ĐSCHT được xử lý để loại bỏ nhiễu và ngoại suy các giá trị ban đầu. Quá trình nghịch đảo số liệu được thực hiện để xác định phân bố nước ngầm trong lòng đất. Các phần mềm như SAMOVAR được sử dụng để giải bài toán nghịch trong môi trường 1D. Tuy nhiên, trong các môi trường phức tạp như đá vôi karst, việc áp dụng các mô hình 2D/3D vẫn đang trong giai đoạn nghiên cứu và phát triển.

III. Kết quả nghiên cứu thử nghiệm tại miền Bắc Việt Nam

Nghiên cứu đã tiến hành thử nghiệm ĐSCHT tại hai khu vực trọng điểm ở miền Bắc Việt Nam: thị xã Lai Châu và Bó Mười, Thuận Châu, Sơn La. Các kết quả cho thấy hiệu quả của phương pháp trong việc xác định nước ngầm tại các vùng đá vôi karst. Tại thị xã Lai Châu, ĐSCHT đã xác định được các tầng chứa nước với độ sâu và phân bố rõ ràng. Tại Bó Mười, phương pháp này cũng cho kết quả khả quan, giúp xác định các vị trí khoan nước ngầm tiềm năng. Các kết quả này đã chứng minh tính ứng dụng cao của ĐSCHT trong việc tìm kiếm và khai thác tài nguyên nước tại các khu vực địa chất phức tạp.

3.1 Kết quả tại thị xã Lai Châu

Tại thị xã Lai Châu, ĐSCHT đã được áp dụng để khảo sát nước ngầm trong điều kiện địa chất phức tạp. Các kết quả đo sâu phân cực và cộng hưởng từ đã xác định được các tầng chứa nước với độ sâu từ 20 đến 50 mét. Các kết quả khoan thực tế đã xác nhận tính chính xác của phương pháp, mở ra tiềm năng khai thác nước ngầm phục vụ nhu cầu kinh tế - xã hội của địa phương.

3.2 Kết quả tại Bó Mười Thuận Châu Sơn La

Tại Bó Mười, ĐSCHT đã được áp dụng để khảo sát nước ngầm trong các tầng đá vôi karst. Các kết quả đo sâu cộng hưởng từ đã xác định được các vị trí có tiềm năng chứa nước ngầm, với độ sâu từ 30 đến 60 mét. Các kết quả này đã được xác nhận thông qua các đợt khoan thực tế, chứng minh hiệu quả của phương pháp trong việc tìm kiếm nước ngầm tại các khu vực địa chất phức tạp.

Nội dung chính

Nghiên Cứu Áp Dụng Phương Pháp Đo Sâu Cộng Hưởng Từ (MRS) Tìm Kiếm Nước Ngầm Trọng Điểm Miền Bắc Việt Nam

---

Tóm tắt nghiên cứu

• Câu hỏi nghiên cứu chính (Research Question): Phương pháp đo sâu cộng hưởng từ (Magnetic Resonance Sounding - MRS) có khả năng phát hiện, định lượng trực tiếp và nâng cao độ chính xác trong công tác thăm dò nguồn nước dưới đất tại các cấu trúc địa chất phức tạp, đặc biệt là địa hình đá vôi karst tại vùng trọng điểm miền Bắc Việt Nam hay không?
• Phương pháp nghiên cứu (Methodology Snapshot): Đề tài ứng dụng hệ thống thiết bị đo sâu cộng hưởng từ hạt nhân NUMIS Plus kết hợp phần mềm nghịch đảo SAMOVAR (dựa trên thuật toán điều chỉnh hóa Tikhonov cho mô hình 1D). Nghiên cứu triển khai các cấu hình khung dây phát - thu chuyên dụng (khung dây số 8, khung dây bù) nhằm triệt tiêu nhiễu điện từ công nghiệp và tự nhiên, đồng thời tích hợp đo sâu phân cực kích thích (IP) và kiểm chứng đối sánh bằng các lỗ khoan thực nghiệm tại thị xã Lai Châu và khu vực Bó Mười (Thuận Châu, Sơn La).
• Kết quả then chốt (Key Findings): MRS chứng minh tính ưu việt vượt trội khi bắt được tín hiệu trực tiếp từ proton ($^1\text{H}$) của phân tử nước tự do, ước lượng chính xác độ sâu đỉnh tầng chứa nước và hàm lượng nước ($w%$) với sai số dưới 15–20% tại các đới karst hóa. Đồng thời, nghiên cứu làm rõ các giới hạn vật lý của mô hình nghịch đảo 1D khi áp dụng vào dị thường 2D/3D (hiện tượng lớp nước ảo, hiệu ứng màn chắn).
• Ý nghĩa thực tiễn (Implications): Cung cấp cơ sở khoa học và quy trình công nghệ chuẩn mực để chuyển giao kỹ thuật thăm dò nước ngầm không xâm lấn, mở ra giải pháp giải quyết triệt để bài toán thiếu nước sinh hoạt và sản xuất tại các vùng núi đá vôi hiểm trở phía Bắc.

---

Bối cảnh và tầm quan trọng

flowchart TD
    A["Thách thức: Thiếu nước vùng núi đá vôi miền Bắc"] --> B["Hạn chế phương pháp truyền thống: Chỉ đo gián tiếp (Điện trở, Tốc độ truyền sóng)"]
    B --> C["Giải pháp: Đo sâu cộng hưởng từ (MRS)"]
    C --> D["Nguyên lý: Tương tác Proton H1 trong nước tự do"]
    D --> E["Kết quả: Xác định trực tiếp vị trí, độ sâu, thể tích tầng chứa nước"]

Hiện trạng tri thức và hạn chế của các phương pháp truyền thống

Trong suốt nhiều thập kỷ, công tác địa vật lý tìm kiếm nước dưới đất tại Việt Nam và trên thế giới phụ thuộc chủ yếu vào các phương pháp đo gián tiếp như: đo sâu điện trở suất (DC), phân cực kích thích (IP), điện từ (EM) hoặc địa chấn khúc xạ. Các phương pháp này đo lường tham số môi trường đất đá (điện trở suất, độ phân cực, tốc độ truyền sóng đàn hồi) thay vì đo trực tiếp lượng nước. Do đó, trong môi trường địa chất phức tạp như vùng đá vôi karst hóa cao, sự nhầm lẫn giữa tầng sét dẫn điện, đới dập vỡ chứa bùn ẩm và hang karst chứa nước thực thụ thường xuyên xảy ra, dẫn đến tỷ lệ khoan thăm dò thất bại cao.

Khoảng trống nghiên cứu được nhận diện

Kể từ khi phương pháp đo sâu cộng hưởng từ bề mặt (MRS/SNMR) ra đời vào cuối thập niên 1980 tại Nga và phát triển mạnh tại Pháp, Đức, Mỹ, công nghệ này đã chứng minh khả năng đo đạc trực tiếp tín hiệu từ hạt nhân hydro trong nước. Tuy nhiên, phần lớn các ứng dụng quốc tế tập trung vào môi trường trầm tích bở rời (chiếm khoảng 39% các đợt khảo sát toàn cầu). Việc áp dụng MRS cho môi trường đá nứt nẻ và karst (chỉ chiếm 9% trên thế giới) đối mặt với rào cản kỹ thuật cực lớn:

• Cấu trúc chứa nước dạng ống/hang ngầm mang tính chất 2D/3D phi đối xứng.
• Tỷ số tín hiệu trên nhiễu ($S/N$) cực thấp do tín hiệu cộng hưởng từ chỉ ở mức nanovolt ($\text{nV}$).
• Địa hình hiểm trở và nhiễu điện từ tự nhiên/nhân tạo tại miền núi phía Bắc Việt Nam.

Tính cấp thiết và tác động tiềm năng

Miền núi Tây Bắc (Lai Châu, Sơn La) là những khu vực thường xuyên chịu hạn hán khốc liệt trong mùa khô, địa hình karst gây thất thoát nước mặt nghiêm trọng. Đề tài cấp Bộ do Viện Khoa học Địa chất và Khoáng sản chủ trì có tính thời sự cấp bách: vừa làm chủ công nghệ đo sâu tiên tiến nhất thế giới, vừa xây dựng quy trình kỹ thuật hoàn chỉnh nhằm định vị chính xác túi nước ngầm, phục vụ an sinh xã hội và phát triển kinh tế bền vững.

---

Phương pháp luận và cách tiếp cận (Methodology & Approach)

graph LR
    subgraph Thu_thap [1. Thu thập dữ liệu thực địa]
        A[Xác định trường từ tĩnh B0] --> B[Phát xung xoay chiều tần số Larmor]
        B --> C[Chống nhiễu: Khung số 8 / Khung bù]
    end
    subgraph Xu_ly [2. Xử lý & Nghịch đảo]
        C --> D[Bộ tách sóng đồng bộ & Lọc thông thấp]
        D --> E[Nghịch đảo Tikhonov 1D trên SAMOVAR]
    end
    subgraph Kiem_chung [3. Đánh giá & Kiểm chứng]
        E --> F[Mô hình hóa số 2D/3D]
        F --> G[Khoan kiểm chứng thực tế]
    end

Cơ sở vật lý và nguyên lý kích thích cộng hưởng từ

Phương pháp MRS dựa trên mômen từ hạt nhân của proton ($^1\text{H}$) trong phân tử nước tự do:

1. Trạng thái cân bằng: Dưới tác dụng của trường địa từ tĩnh ($\vec{B}_0$), các spin proton định hướng song song với $\vec{B}_0$, quay đảo quanh trục từ trường với tần số góc Larmor $\omega_0 = \gamma |\vec{B}_0|$ (trong đó $\gamma$ là tỷ số từ hồi chuyển của proton, $\gamma \approx 2{,}675 \times 10^8 \text{ rad}\cdot\text{s}^{-1}\cdot\text{T}^{-1}$).
2. Kích thích cộng hưởng: Máy phát dòng điện xoay chiều $I(t) = I_0 \cos(\omega_0 t)$ vào khung dây trên mặt đất với thời gian phát $\tau_p$, tạo mômen xung $q = I_0 \tau_p$ sinh ra từ trường kích thích $\vec{B}T$. Thành phần $\vec{B}{T\bot}$ vuông góc với $\vec{B}_0$ làm lệch vectơ từ hóa của proton khỏi trục cân bằng một góc $\theta$.
3. Suy giảm cảm ứng tự do (FID): Khi ngắt xung kích thích, các proton hồi phục về trạng thái ban đầu, phát ra tín hiệu điện từ cảm ứng suy giảm theo thời gian với hằng số $T_2^*$ (suy giảm ngang) và $T_1$ (suy giảm dọc). Biên độ ban đầu $E_0$ tỷ lệ thuận trực tiếp với số lượng hạt nhân hydro (hàm lượng nước $w%$).

Thiết kế thu thập số liệu và kỹ thuật triệt tiêu nhiễu

Tín hiệu cộng hưởng từ rất nhỏ (từ $10\text{ nV}$ đến vài trăm $\text{nV}$), dễ bị vùi lấp bởi nhiễu công nghiệp và giông sét. Đề tài áp dụng hệ thống NUMIS Plus với các giải pháp kỹ thuật tối ưu:

• Khung dây số 8 (Figure-eight loop): Gồm 2 vòng dây đấu ngược chiều. Nhiễu điện từ từ nguồn xa tạo ra 2 suất điện động cảm ứng ngược pha tự triệt tiêu lẫn nhau, trong khi tín hiệu MRS cục bộ được giữ lại và cộng gộp, nâng tỷ số $S/N$ lên gấp nhiều lần.
• Khung dây bù (Compensation loop - Cx): Đặt khung phụ cách khung phát một khoảng $2{,}5L$ theo phương xác định để thu riêng tín hiệu nhiễu và trừ trực tiếp khỏi tín hiệu khung chính.
• Tách sóng đồng bộ (Synchronous Detection) & Lọc thông thấp: Tách tín hiệu tần số Larmor dao động hẹp ($\pm 1\text{ Hz}$) ra khỏi dải phổ nhiễu ngẫu nhiên.
• Số lần cộng dồn ($n$): Tối ưu từ $n = 32$ đến $n = 128$ chu kỳ đo tùy thuộc vào biên độ nhiễu thực tế tại hiện trường.

Kỹ thuật xử lý nghịch đảo và mô hình hóa toán học

• Thuật toán nghịch đảo 1D: Ứng dụng phương pháp điều chỉnh hóa Tikhonov (Tikhonov regularization) tích hợp trong phần mềm SAMOVAR để giải bài toán ngược phi tuyến:

$$\min_{w} \left{ | \mathbf{A} w - e_0 |{L_2}^2 + \eta | \mathbf{W}\eta |_{L_2}^2 \right}$$

Trong đó ma trận hạt nhân $\mathbf{A}$ biểu thị độ nhạy không gian của khung dây, $e_0$ là tập số liệu đo thực nghiệm, và $\eta$ là tham số điều chỉnh độ trơn của phân bố hàm lượng nước $w(z)$.

• Phân tích mô hình số 2D nâng cao: Sử dụng các mô hình cấu trúc 2D (thấu kính nước, lớp nước vát nhọn, ống dẫn karst) để đánh giá độ chệch và nhận diện các hiệu ứng tương đương hình học khi áp dụng thuật toán nghịch đảo 1D.

---

Phát hiện chính của nghiên cứu

stateDiagram-v2
    [*] --> ThuThapTinHieu: Do tai mat dat
    ThuThapTinHieu --> KiemTraChatLuong: Danh gia ty so S/N va EN/IN
    KiemTraChatLuong --> TachSongVaNghichDao: S/N > 2, Loc so bo
    TachSongVaNghichDao --> PhanTichHangSoT2: T2* > 400-500 ms (Dau hieu Karst)
    TachSongVaNghichDao --> XacDinhDoSau: Dinh lop nuoc & Ham luong w%
    PhanTichHangSoT2 --> DoiChieuKhoanKiemChung: Khang dinh vi tri hang ngam
    XacDinhDoSau --> DoiChieuKhoanKiemChung
    DoiChieuKhoanKiemChung --> [*]

1. Khả năng phát hiện trực tiếp và chính xác nguồn nước ngầm karst

Tại phường Quyết Thắng và xã Nậm Loỏng (thị xã Lai Châu), các điểm đo MRS (như cọc 5 tuyến T1b) đã ghi nhận biên độ tín hiệu ban đầu rõ nét với hằng số suy giảm $T_2^* > 400\text{ ms}$ – $500\text{ ms}$. Đây là dấu hiệu đặc trưng khẳng định nước tồn tại trong các khoảng rỗng lớn (nước tự do trong hang karst/khe nứt lớn) chứ không phải nước liên kết trong sét ($T_2^* < 30\text{ ms}$). Kết quả khoan kiểm chứng thực tế tại cọc 5 tuyến T1b đã bắt gặp đới karst chứa nước áp lực ở chiều sâu dự báo, lưu lượng khai thác thực tế đạt giá trị cao.

2. Định lượng ranh giới độ sâu và hàm lượng nước

Nghiên cứu chỉ ra rằng MRS xác định độ sâu nóc (đỉnh) tầng chứa nước với độ tin cậy rất cao (sai số $\varepsilon < 15%$). Tại các điểm đo thử nghiệm ở Bản Bó (Thuận Châu, Sơn La), các đường cong phân bố $w(z)$ nghịch đảo từ SAMOVAR thể hiện rõ nét vị trí phân tầng của các tầng ngậm nước nông (đới phong hóa nứt nẻ $5\text{ m} - 15\text{ m}$) và tầng karst sâu ($30\text{ m} - 60\text{ m}$).

3. Phát hiện hiện tượng "Tầng chứa nước ảo" (Phantom Aquifer) do dịch chuyển tần số

Một đóng góp vật lý quan trọng là việc làm rõ tác động của độ dịch chuyển tần số Larmor $\Delta\omega = \omega_0 - \omega_{rf}$ (do biến thiên nhật triều của từ trường Trái Đất hoặc gradient từ địa phương) kết hợp với các sóng hài bậc cao của xung kích thích. Khi $\Delta\omega \neq 0$, nếu chỉ giải thích đơn thuần đường cong biên độ ở các mômen xung lớn ($q > 5000\text{ A}\cdot\text{ms}$) mà không kiểm tra độ lệch pha âm, người xử lý rất dễ nhận định sai lầm về sự tồn tại của một "lớp nước thứ hai" dưới sâu.

4. Nhận diện các giới hạn của mô hình 1D trong môi trường 2D/3D

Qua mô hình hóa lý thuyết và đối sánh thực địa:

• Đối với dị thường chứa nước dạng hình hộp/ống karst hẹp có chiều rộng nhỏ hơn đường kính khung dây ($W < d$), nghịch đảo 1D cho ra giá trị hàm lượng nước $w%$ nhỏ hơn thực tế khoảng 2 lần.
• Hiệu ứng màn chắn (Shielding effect): Tầng nước bão hòa nằm nông có hàm lượng nước lớn gấp $\ge 5$ lần tầng sâu sẽ triệt tiêu phần lớn đáp ứng của tầng sâu, khiến thuật toán 1D không thể phân giải được tầng chứa nước khe nứt bên dưới.

---

Đóng góp khoa học và giá trị gia tăng

Khía cạnh đóng góp	Chi tiết giá trị khoa học & Thực tiễn
Đóng góp lý thuyết	• Hệ thống hóa cơ sở lý thuyết truyền sóng điện từ phân cực elip trong môi trường dẫn điện 3D. • Phân tích bản chất vectơ của tín hiệu MRS khi độ từ khuynh nghiêng ($I \neq 90^\circ$), chứng minh sự phá vỡ nguyên lý tương hoán (reciprocity principle).
Đổi mới phương pháp luận	• Xây dựng quy trình phối hợp tối ưu giữa đo điện trở suất/IP truyền thống và đo sâu MRS. • Xác lập tiêu chuẩn lựa chọn khung dây (khung số 8 vs khung bù) và thông số cộng dồn $n$ thích ứng với mức ồn địa phương.
Ứng dụng thực tiễn	• Biên soạn dự thảo "Quy trình công nghệ đo sâu cộng hưởng từ tìm kiếm nước ngầm vùng đá vôi karst" (Tiêu chuẩn cơ sở) gồm 10 phụ lục hướng dẫn chi tiết. • Giảm thiểu rủi ro và chi phí khoan thăm dò địa chất tại các tỉnh miền núi phía Bắc.
Hàm ý chính sách & Quản lý	• Cung cấp công cụ khoa học chính xác cho Bộ Tài nguyên và Môi trường trong việc lập bản đồ quy hoạch tài nguyên nước ngầm. • Đề xuất chuẩn hóa phương pháp MRS vào danh mục khảo sát địa vật lý bắt buộc tại các vùng khan hiếm nước đặc thù.

---

Đối tượng quan tâm và giá trị ứng dụng

• Các nhà nghiên cứu Địa vật lý & Địa chất thủy văn: Tiếp cận hệ thống lý thuyết chuẩn xác về nghịch đảo cộng hưởng từ hạt nhân bề mặt, các thuật toán xử lý bài toán thuận/nghịch đa chiều và phương pháp phân tích hằng số suy giảm $T_1$, $T_2^*$.
• Kỹ sư khảo sát & Doanh nghiệp thăm dò khai thác nước ngầm: Nắm vững quy trình vận hành thực địa máy NUMIS Plus, kỹ thuật bố trí khung dây chống nhiễu, cách lựa chọn bước xung $q$ và phương pháp đọc biểu đồ minh giải để định vị chính xác vị trí đặt giếng khoan.
• Cơ quan quản lý Nhà nước (Bộ TN&MT, Sở TN&MT các tỉnh miền núi): Sở hữu căn cứ khoa học vững chắc để phê duyệt các dự án cấp nước sinh hoạt nông thôn, tối ưu hóa nguồn vốn ngân sách đầu tư cho hạ tầng thủy lợi và nước sạch.

---

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

1. Điểm khác biệt cốt lõi giữa đo sâu cộng hưởng từ (MRS) và các phương pháp địa điện truyền thống là gì?

Các phương pháp địa điện (đo điện trở, đo sâu phân cực) chỉ đo gián tiếp các tính chất vật lý của đất đá (vốn dễ bị nhầm lẫn giữa tầng sét ướt và tầng chứa nước). Ngược lại, MRS kích thích trực tiếp hạt nhân nguyên tử hydro ($^1\text{H}$) trong phân tử nước tự do, mang lại tín hiệu đáp ứng chỉ xuất hiện khi có nước, loại bỏ hoàn toàn hiện tượng "dị thường giả" do đất đá dẫn điện.

2. Tại sao khảo sát nước trong đá karst lại khó khăn hơn nhiều so với trầm tích bở rời?

Trầm tích bở rời (cát, cuội sỏi) có cấu trúc phân lớp nằm ngang tương đối đồng nhất (mô hình 1D). Trong khi đó, nước ngầm trong đá karst phân bố dạng hang hốc, khe nứt bất đối xứng (mô hình 2D/3D), lượng nước cục bộ nhỏ và tín hiệu suy giảm rất phức tạp, đòi hỏi phải áp dụng các kỹ thuật chống nhiễu chuyên sâu và mô hình hóa nâng cao.

3. Kết quả nghiên cứu của đề tài có thể khái quát hóa (generalize) cho các vùng địa chất khác không?

Có. Mặc dù các thử nghiệm thực tế tập trung tại Lai Châu và Sơn La, quy trình công nghệ, phương pháp chống nhiễu bằng khung dây số 8 và nguyên lý giải tích dữ liệu trong báo cáo hoàn toàn có thể áp dụng cho mọi vùng địa hình karst, đới phong hóa đá magma (granit, bazan) hoặc bồn trũng trầm tích trên toàn lãnh thổ Việt Nam.

4. Thiết bị NUMIS Plus có phát hiện được nước ở độ sâu bao nhiêu và có phân biệt được nước mặn/nhạt không?

Độ sâu khảo sát tối đa của NUMIS Plus đạt khoảng $100\text{ m} - 150\text{ m}$ (tùy thuộc kích thước khung dây và mức độ dẫn điện của đất đá). Ở giai đoạn công nghệ hiện tại, phương pháp MRS bề mặt chưa phân biệt được độ khoáng hóa (nước mặn hay nước ngọt), cần kết hợp với phương pháp đo điện trở suất để đánh giá chất lượng nước.

5. Hướng phát triển tiếp theo của công nghệ này tại Việt Nam là gì?

Cần nâng cấp nghiên cứu thuật toán nghịch đảo 2D/3D đa kênh (multi-channel MRS) với hệ thống đa khung thu, tích hợp xử lý đa hằng số suy giảm (multi-exponential decay) để đánh giá đồng thời cả thể tích trữ lượng và độ thấm thủy lực của tầng chứa nước.

---

Kết luận

Báo cáo nghiên cứu khoa học của TS. Tăng Đình Nam và nhóm tác giả Viện Khoa học Địa chất và Khoáng sản là một công trình mang tính tiên phong, đặt nền móng vững chắc cho việc ứng dụng công nghệ đo sâu cộng hưởng từ hạt nhân tại Việt Nam. Đề tài không chỉ làm chủ thiết bị hiện đại NUMIS Plus mà còn xây dựng thành công quy trình công nghệ chuẩn mực, giải quyết xuất sắc bài toán tìm kiếm trực tiếp nguồn nước ngầm trong cấu trúc địa chất karst hiểm trở tại Tây Bắc.

Trong bối cảnh biến đổi khí hậu và an ninh nguồn nước đang trở thành thách thức sống còn, việc nhân rộng quy trình công nghệ MRS vào các chương trình điều tra cơ bản địa chất và cấp nước sinh hoạt quốc gia là hướng đi chiến lược, mang lại hiệu quả kinh tế – xã hội to lớn và lâu dài.

Trích đoạn nội dung tài liệu

CHƯƠNG 1 TỔNG QUAN VỀ ĐO SÂU CỘNG HƯỞNG TỪ Đo sâu cộng hưởng từ (viết tắt: ĐSCHT; tiếng Anh: Magnetic Resonance Sounding, viết tắt: MRS) dựa trên hiện tượng cộng hưởng từ hạt nhân và hiện là phương pháp địa vật lý duy nhất trên thế giới được dùng để khảo sát trực tiếp nước ngầm từ mặt đất. Tất cả các phương pháp địa vật lý khác dùng trong tìm kiếm nước ngầm đều là gián tiếp theo các tham số vật lý như điện trở suất, các tham số phân cực kích thích, hằng số điện môi, tốc độ truyền sóng địa chấn,….G và các cộng tác viên ở Nga đã nghiên cứu phương pháp, giải bài toán thuận và nghịch cho ĐSCHT bằng khung dây đơn (một khung dây làm chức năng phát trong quá trình phát sóng điều hoà đơn sắc, sau đó nhờ bộ chuyển mạch để làm chức năng thu), trong môi trường điện trở suất cao, sáng chế máy cộng hưởng từ đầu tiên phục vụ khảo sát nước ngầm mang tên Hydroscope, nhận bằng sáng chế ở Liên Xô năm 1988 [61] và ở Anh năm 1989 [62]. Đến năm 1999 và 2000 Weichman và nnk công bố công trình nghiên cứu lý thuyết tổng quát về ĐSCHT trong môi trường dẫn điện, tính tới hiện tượng phân cực elip cho loại khung dây bất kỳ (khung dây đơn, khung dây tách biệt,. Lần đầu tiên việc đánh giá các tính chất của các khung dây tách biệt và các tính chất của các đường cong ĐSCHT tương ứng được trình bày vào năm 2003 [19] và thể hiện đầy đủ trong luận án tiến sỹ của Hertrich M năm 2005 [17], cho thấy độ phân giải cao trong khảo sát môi trường 2D nông.

Những nghiên cứu lý thuyết về ĐSCHT trên thế giới tập trung chủ yếu ở bốn trung tâm là Nga, Pháp, Đức và Mỹ. Sau đây sẽ trình bày cơ sở vật lý của ĐSCHT , cơ sở toán học - kết quả nghiên cứu lý thuyết bằng mô hình hoá toán học dựa vào những công trình của bốn trung tâm nêu trên cho đến năm 2008. Một số bài báo tiếng Việt đã giới thiệu ĐSCHT [10,45,46] và kết quả thử nghiệm đầu tiên ở Việt Nam [47].1 Nguyên lý cơ bản và bài toán thuận của đo sâu cộng hưởng từ 1. Nguyên lý cơ bản của đo sâu cộng hưởng từ Hydro có 3 đồng vị: 1H, 2H và 3H, trong đó 1H phong phú nhất, chiếm 99.

Hydro có trong phân tử nước. Trong hạt nhân của 1H chỉ có một proton và đó là đối tượng khảo sát trực tiếp nước ngầm của ĐSCHT. Proton là hạt tích điện dương, có khối lượng đặc trưng, có momen từ µ (ký hiệu chữ hạt 4 tích điện dương, có khối lượng đặc trưng, có momen từ µ (ký hiệu chữ đậm là vectơ) và có spin. Spin là đại lượng mang tính lượng tử, đặc trưng cho hạt vi mô, có thể coi đó là momen động lượng quay nội tại của hạt vi mô, thường gọi là momen spin, nó mô tả bản chất trường từ của hạt vi mô.

Momen từ liên hệ cộng tuyến với momen spin I theo biểu thức: µ = γh I (1.1) Trong đó : h I là xung lượng góc của hạt nhân ( h là hằng số Plank rút gọn: h = h/2 π = 1,054592. 10-34 J s, h là hằng số Plank), γ là tỉ số từ hồi chuyển (có giá trị đặc trưng cho mỗi hạt nhân, đối với hạt nhân hydro 8 -1 -1 γ = 2,675. Khi các proton nằm trong trường từ tĩnh bên ngoài, trong ĐSCHT là trường địa từ B0, thì các momen từ định hướng theo phương z của trường từ bên ngoài, các momen từ quay đảo quanh trục z với tần số vòng Larmor ω0 (tần số Larmor f 0 ) và ở trạng thái cân bằng thì tạo ra độ từ hoá M0 được xác định theo các biểu thức sau: (1.3) Trong đó kB là hằng số Boltzmann, T là nhiệt độ tuyệt đối, N/V là số spin trong đơn vị thể tích. Dấu âm ở biểu thức của tần số vòng Larmor ω0 thể hiện sự quay đảo của các momen từ quanh trục z theo chiều kim đồng hồ (Levitt, 1997).

Khi các momen từ đã nằm trong trường từ tĩnh B0 người ta phát dòng xoay chiều với tần số vòng Larmor: I(t) = I0cos( ω0 t) vào khung dây phát trong một thời gian ngắn τ p (xung điện từ) để tạo ra trường từ điều hoà cộng hưởng BT. Dưới tác dụng của BT, các momen từ quay đảo sẽ lệch khỏi phương z của trường từ tĩnh B0 một góc θ về phương của trường BT. Do trường điện từ điều hoà truyền trong môi trường đất đá dẫn điện nên trường từ kích thích BT sẽ phân cực elip và tạo ra lệch pha giữa trường đã truyền đi (transmitted field) so với trường kích thích do sự dẫn điện của môi trường. Các nghiên cứu trước đây của Trushkin và nnk (1995) và Shushakov (1996) đã xem xét điện trở suất của đất đá trong khi tính biên độ và pha của tín hiệu cộng hưởng từ và mô tả tín hiệu cộng hưởng từ có giá trị phức, nhưng lại không khảo sát phân cực elip mà dùng phép xấp xỉ cho trường hợp độ từ khuynh 900 thích hợp cho vùng gần cực địa từ ở Xiberi, Nga.

Weichman và nnk (1999, 2000) đã khảo sát đầy đủ sự phân cực elip của trường kích thích và biểu thức của tín hiệu phức cho loại khung dây và vị trí địa lý bất kỳ. 5 Weichman và nnk (2000) đã chứng tỏ rằng chỉ có thành phần BT⊥ của trường kích thích BT vuông góc với trường địa từ B0 mới có tác dụng tới các spin. Trong quá trình tính toán B T⊥ phân cực elip được phân tích thành hai phân cực tròn với các biên độ khác nhau và quay trái chiều nhau B T⊥ và B T_ , trong đó B T+ quay theo chiều kim đồng hồ và như vậy là quay cùng chiều với momen từ, còn B -T quay ngược chiều kim đồng hồ quanh trường từ tĩnh B0. Trong quá trình ĐSCHT người ta tăng dần độ sâu khảo sát tại mỗi điểm đo bằng cách tăng momen xung q (q = I0 τ p ).

Momen xung q trong ĐSCHT đóng vai trò gần tương tự như khoảng cách giữa 2 điện cực phát AB trong đo sâu điện thẳng đứng. Sau khi ngắt xung, các momen từ hồi phục, suy giảm dần để trở lại trạng thái ban đầu theo phương của trường địa từ B0 và tạo ra tín hiệu cộng hưởng từ (điện áp) trong khung dây thu. Điện áp đo được ngay sau ngắt xung càng lớn khi hạt nhân hydro càng nhiều. Chúng có trong phân tử nước, dầu khí (đối tượng khảo sát của carota cộng hưởng từ) và cả trong phân tử ở thể rắn.

Chúng được phân biệt dựa vào đặc điểm hồi phục khác nhau, theo hằng số thời gian suy giảm T khác nhau do mức độ linh động khác nhau. Có hai loại T tương ứng với hai loại vectơ từ hoá, đó là T1 đối với vectơ từ hoá dọc (dọc theo phương của trường địa từ B0) và T2 đối với vectơ từ hoá ngang (vuông góc với phương của trường địa từ B0). Sự suy giảm đó được gọi là suy giảm cảm ứng tự do (Free Induction Decay, FID) và được dùng trong ĐSCHT. Trong thực tế thường hay gặp loại suy giảm cảm ứng tự do nhanh chóng hơn, gọi là suy giảm T2* do tính bất đồng nhất của trường vật lý cũng như sự biến đổi độ cảm từ trong môi trường.

Các hằng số thời gian suy giảm tuân theo quan hệ sau : T2* ≤ T2 ≤ T1 [12,29]. Hạt nhân trong phân tử chất rắn có mức độ linh động rất hạn chế và suy giảm rất nhanh (T2 rất nhỏ), như bitum (một loại hydrocarbon ở thể rắn) có T2 ≤ 0,1 ms. Trong đất đá lỗ hổng như cát kết, trầm tích vỡ vụn giàu thạch anh, các hạt nhân trong phân tử chất lỏng liên kết với bề mặt chất rắn (nước liên kết) hoặc trong lỗ hổng nhỏ khi dao động luôn va chạm với bề mặt chất rắn, tiêu hao năng lượng nhiều nên năng lượng suy giảm nhanh (T2 nhỏ), hạt nhân ở trong hổng lớn hơn (nước tự do) có T2 lớn hơn [11,58]. Máy đo có thời gian chết (“dead time”, máy chưa đo được) 30 ÷ 40 ms, nên chỉ đo được tín hiệu của nước tự do.

Đó là cơ sở vật lý của việc xác định trực tiếp nước ngầm của ĐSCHT ở mặt đất. Tín hiệu dao động với tần số vòng Larmor, suy giảm theo quy luật hàm mũ với hằng số thời gian T2* , lệch pha ϕ 0 so với trường phát và phụ thuộc vào momen xung q như sau: 6 (1.4) Tín hiệu ban đầu E0 và ϕ 0 (ngay sau khi ngắt xung: t = 0) không đo được do có thời gian chết τ d của máy, nên phải ngoại suy để tính ra theo các biểu thức sau (Legchenko A.6) Trong đó: ∆ω ( ∆ω = ω0 − ωrf ) được gọi là dịch chuyển tần số (frequency shift hoặc frequency offset), là hiệu số giữa tần số Larmor vòng ω0 của trường kích thích và tần số vòng của tín hiệu cộng hưởng từ đo được ω rf. Bốn tham số E( τ d ,q), T2* (q), ϕ ( τ d ,q) và ∆ω (q) được tính theo các thành phần cùng pha và lệch pha trong quá trình xử lý số liệu (Legchenko A. Các tham số đó được dùng để ngoại suy theo các biểu thức (1.5) từ các giá trị đo được sau thời gian chết τ d thành giá trị biên độ (1.5) ban đầu ở thời điểm t = 0.

Đường ghi biên độ suy giảm với thời gian theo hàm mũ, còn pha thay đổi tuyến tính. Các kết quả khảo sát đều dùng biên độ và pha ban đầu. Thế phức = liên hệ tuyến tính với phân bố hàm lượng nước f(r) qua hàm nhân : (1.7) Trong đó hàm nhân phức phụ thuộc vào độ từ khuynh, điện trở suất ρ (V ) trong thể tích khảo sát, phương tuyến đo so với hệ toạ độ từ và hệ quan sát. Bài toán thuận của đo sâu cộng hưởng từ 1.

Tín hiệu cộng hưởng từ tổng quát Bài toán thuận của ĐSCHT tổng quát cho môi trường dẫn điện, độ từ thẩm bằng độ từ thẩm của không khí, ba chiều (3D) với các khung dây bất kỳ được trình bầy đầu tiên trong công trình của Weichman và nnk (1999, 2000). Sau đây là biểu thức thế V của tín hiệu ĐSCHT (Weichman và nnk, 2000) ở thời điểm t = 0: 7 (1.