Tổng quan nghiên cứu

Khảo sát công trình ngầm và phát hiện dị thường địa chất là bài toán cấp thiết trong xây dựng hạ tầng đô thị hiện đại. Theo ước tính từ các nghiên cứu địa kỹ thuật, hơn 80% các sự cố sụt lún "hố tử thần" hoặc hư hại đường ống cáp ngầm có thể được phòng tránh nếu phát hiện kịp thời bằng phương pháp không phá hủy. Radar xuyên đất (Ground Penetrating Radar - GPR) là kỹ thuật tiên tiến ứng dụng sóng điện từ tần số cao từ 1 MHz đến 1000 MHz để quét và tái tạo cấu trúc lòng đất ở độ sâu khảo sát từ 0,1 m đến 40 m mà không cần đào bới hay phá vỡ kết cấu hiện hữu.

Tuy nhiên, tín hiệu phản xạ GPR trong thực tế thường xuyên chịu ảnh hưởng nghiêm trọng bởi hiện tượng suy hao năng lượng nhanh, trôi dạt điện áp một chiều (DC offset), hiệu ứng rung giao thoa (ringing) giữa anten phát - thu, và nhiễu nền phức tạp do môi trường phân lớp không đồng nhất gây ra. Mục tiêu cụ thể của luận văn là nghiên cứu các đặc tính truyền sóng, đề xuất và thử nghiệm chuỗi thuật toán xử lý tín hiệu số kết hợp xử lý ảnh chuyên sâu (gồm hiệu chỉnh vị trí zero, lọc Dewow, khuếch đại tự động AGC, giải chập Wiener và kỹ thuật di trú Migration) nhằm nâng cao chất lượng hiển thị trực quan và độ chính xác của ảnh radar.

Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong mốc thời gian từ tháng 01/2013 đến tháng 06/2013. Kết quả mô phỏng trên nền tảng phần mềm chuyên dụng chứng minh hệ thống xử lý tín hiệu đề xuất giúp triệt giảm trên 75% nhiễu nền, nén xung phản xạ và định vị mục tiêu dị thường dưới lòng đất với độ chính xác đạt trên 90%, tạo tiền đề vững chắc cho việc làm chủ công nghệ thiết bị GPR tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Khung lý thuyết của luận văn được xây dựng dựa trên hai nền tảng toán học và vật lý cốt lõi:

Thứ nhất là hệ phương trình Maxwell về điện từ trường trong môi trường đồng nhất và tuyến tính, mô tả quy luật lan truyền và suy giảm của sóng điện từ trong lòng đất. Vận tốc truyền sóng phụ thuộc trực tiếp vào hằng số điện môi tương đối của vật chất, dao động từ 1 (không khí, vận tốc khoảng 300 m/ns) đến 80 (nước ngọt, vận tốc khoảng 33 m/ns; hoặc cát khô từ 3 đến 5 với vận tốc 150 m/ns). Khi sóng gặp mặt phân cách giữa hai môi trường có độ dẫn điện hoặc hằng số điện môi chênh lệch, hiện tượng phản xạ và tán xạ sẽ xảy ra, tạo ra tín hiệu hồi tiếp ghi nhận bởi anten thu.

Thứ hai là lý thuyết xử lý tín hiệu số và mô hình tích chập một chiều/hai chiều. Tín hiệu radar thu được tại bề mặt được mô hình hóa dưới dạng tích chập giữa sóng xung cơ bản của máy phát với hàm đáp ứng xung của các lớp địa chất, cộng với thành phần nhiễu môi trường. Trên cơ sở đó, các bộ lọc tối ưu Wiener, giải chập dự đoán và biến đổi phổ tần số sóng F-K được ứng dụng để tách lọc đáp ứng xung nguyên bản của mục tiêu.

Năm khái niệm then chốt xuyên suốt công trình gồm:

  • Hằng số điện môi tương đối (Relative Permittivity).
  • Độ phân giải phương đứng và phương ngang (Vertical and Horizontal Resolution).
  • Hiện tượng rung anten (Antenna Ringing).
  • Phương pháp điểm sâu chung (Common MidPoint - CMP).
  • Kỹ thuật di trú sóng điện từ (Radar Migration).

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu tổng hợp gồm các chuỗi xung radar phát xạ mô phỏng chuẩn và tập dữ liệu mặt cắt thực tế từ hệ thống thiết bị đo địa vật lý GSSI chuyên dụng. Kích thước mẫu khảo sát bao gồm hàng trăm vết quét (traces) liên tục, mỗi vết quét chứa từ 512 đến 1024 mẫu thời gian rời rạc với tần số lấy mẫu tuân thủ tiêu chuẩn Nyquist (khoảng lấy mẫu thời gian nhỏ hơn hoặc bằng một nửa chu kỳ sóng ngắn nhất).

Phương pháp chọn mẫu dựa trên kỹ thuật quét mặt cắt phản xạ liên tục (Reflection Profiling), di chuyển đồng thời cặp anten phát - thu với khoảng cách cố định dọc theo tuyến khảo sát. Lý do lựa chọn chuỗi phân tích kết hợp giữa miền thời gian và miền tần số là vì dữ liệu GPR thực địa bị suy giảm phi tuyến theo quy luật hàm mũ; việc áp dụng đơn lẻ các bộ lọc tuyến tính không thể phục hồi được các tín hiệu phản xạ ở độ sâu lớn. Timeline nghiên cứu kéo dài liên tục trong 6 tháng (từ tháng 01/2013 đến tháng 06/2013), chia thành các giai đoạn thiết lập mô hình toán, lập trình thuật toán xử lý trên Matlab và đối soát kiểm chứng kết quả trực quan trên các tập dữ liệu thực tế.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và thử nghiệm các thuật toán xử lý tín hiệu đã mang lại bốn phát hiện khoa học quan trọng:

  • Hiệu chỉnh thời gian zero (Time Zero Adjust) và khử DC (DC Removal) đã đồng bộ hóa thành công 100% các vết quét về gốc tọa độ bề mặt thực, loại bỏ hoàn toàn hiện tượng lệch tĩnh biên độ do trôi dạt nhiệt và trễ cáp tín hiệu gây ra.
  • Ứng dụng lọc Dewow kết hợp thuật toán trừ vết quét trung bình (Subtract Mean Trace) giúp triệt tiêu hoàn toàn dải sóng truyền thẳng cực mạnh giữa anten phát và thu, tăng tỷ số tín hiệu trên nhiễu (SNR) của dữ liệu thô lên khoảng 12 dB đến 15 dB mà không làm méo dạng phổ tần số chính.
  • Thuật toán khuếch đại tự động AGC (Automatic Gain Control) theo giá trị biên độ hiệu dụng RMS đã bù đắp xuất sắc mức suy hao năng lượng phân kỳ hình học $1/r^2$. Các tín hiệu phản xạ yếu ở độ sâu từ 2 m đến 5 m được nâng biên độ một cách cân bằng, tăng độ tương phản hình ảnh lên hơn 65% so với ảnh gốc ban đầu.
  • Giải pháp giải chập dự đoán với khoảng dự đoán $n = 20, 50, 80$ kết hợp thuật toán di trú Stolt trong miền F-K đã nén thành công các vệt phản xạ dạng đường cong hyperbol rộng về đúng vị trí không gian thực của vật thể, giúp độ phân giải ngang tăng thêm khoảng 40% đến 50%.

Thảo luận kết quả

Nguyên nhân chính khiến dữ liệu thô GPR bị suy giảm chất lượng là do đất bồi và sét có độ dẫn điện cao (khoảng 1 mS/m đến 1000 mS/m), hoạt động như một bộ lọc thông thấp tự nhiên hấp thụ mạnh các thành phần tần số cao. So sánh với các phương pháp xử lý lọc dải truyền thống trong một nghiên cứu gần đây, chuỗi thuật toán đa tầng trong luận văn chứng minh ưu thế vượt trội khi vừa bảo toàn được thông tin pha của tín hiệu, vừa khôi phục được các phản xạ tầng sâu bị che khuất bởi ringing.

Dữ liệu xử lý được trình bày minh họa trực quan qua biểu đồ A-scan (biến thiên biên độ - thời gian trên từng vết quét đơn lẻ) và ảnh B-scan 2D (mặt cắt quét không gian dọc theo tuyến khảo sát). Bảng so sánh chỉ số năng lượng trước và sau khi giải chập khẳng định năng lượng của xung nguồn đã được nén về dạng xung nhọn cực hẹp tại thời điểm gốc, giảm thiểu trên 70% hiện tượng đa phản xạ giả mạo. Điều này mang ý nghĩa quyết định trong việc giải đoán địa tầng chính xác, giảm thiểu sai số nhận diện vị trí ống dẫn ngầm và vết nứt kết cấu bê tông.

Đề xuất và khuyến nghị

Nhằm tối ưu hóa hiệu quả ứng dụng của công nghệ radar xuyên đất trong thực tế khảo sát tại Việt Nam, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  • Chuẩn hóa quy trình xử lý tín hiệu số đa giai đoạn (tiền xử lý Time Zero - Dewow - AGC - Giải chập Wiener - Di trú F-K) thành quy chuẩn kỹ thuật đồng bộ, hướng đến mục tiêu giảm 80% thời gian xử lý thủ công của kỹ sư trắc địa; lộ trình áp dụng trong vòng 6 tháng tới do các trung tâm đo đạc địa chất và hạ tầng đô thị chủ trì thực hiện.
  • Phát triển module phần mềm tự động hóa giải thuật di trú Stolt và lọc trung vị thích nghi tích hợp trực tiếp trên giao diện Matlab hoặc C++, đặt mục tiêu nâng độ chính xác định vị dị thường lên trên 95% trong vòng 12 tháng tới, giao cho các nhóm nghiên cứu công nghệ thông tin và điện tử đảm trách.
  • Ứng dụng công nghệ radar xuyên đất tần số bước SFCW (Stepped Frequency Continuous Wave) thay thế dần cho hệ thống phát xung truyền thống, nhằm mở rộng dải động lên hơn 20 dB và tăng độ xuyên sâu khảo sát từ 5 m lên 15 m; định hướng thực hiện trong giai đoạn 1 đến 2 năm tới bởi các viện nghiên cứu chuyên ngành điện tử viễn thông.
  • Xây dựng cơ sở dữ liệu số về hằng số điện môi chuẩn cho các loại địa chất đặc thù tại Việt Nam (đất phù sa, đất sét, cát biển, đá granit) với sai số tham số dưới 5%, triển khai trong 18 tháng tới dưới sự phối hợp giữa các trường đại học khối kỹ thuật và các đơn vị kiểm định xây dựng.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chính:

  • Kỹ sư điện tử và chuyên gia xử lý tín hiệu số (DSP): Tiếp cận chi tiết các thuật toán lọc tối ưu Wiener, giải chập dự đoán, khuếch đại AGC và phương pháp biến đổi di trú F-K kèm lưu đồ thuật toán chi tiết để phát triển thuật toán nhúng cho thiết bị radar.
  • Kỹ sư địa chất, trắc địa và khảo sát công trình: Nắm vững các mô hình vật lý lan truyền sóng điện từ, công thức tính toán độ phân giải và phương pháp thu thập dữ liệu mặt cắt để phát hiện chính xác hố sụt ngầm, túi bùn, ranh giới tầng đất đá.
  • Đơn vị quản lý hạ tầng đô thị và nhà thầu xây dựng: Ứng dụng tài liệu để xây dựng quy trình kiểm định không phá hủy chất lượng cầu đường, định vị chính xác mạng lưới đường ống cấp thoát nước, cáp điện ngầm mà không làm ảnh hưởng đến giao thông đô thị.
  • Học viên cao học và sinh viên ngành Kỹ thuật Điện tử - Viễn thông, Địa vật lý: Nguồn tài liệu học tập toàn diện về hệ phương trình Maxwell, tiêu chuẩn lấy mẫu Nyquist trong không gian - thời gian và kỹ thuật lập trình mô phỏng xử lý ảnh GPR 2D trên môi trường Matlab.

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

Hệ thống GPR khảo sát hiệu quả nhất trong dải tần số nào và độ sâu đo đạc là bao nhiêu? Thiết bị GPR thường sử dụng dải tần số vô tuyến từ 1 MHz đến 1000 MHz (phổ biến từ 12,5 MHz đến 400 MHz). Độ sâu khảo sát thực tế dao động từ 0,1 m đến 40 m, phụ thuộc chặt chẽ vào tần số anten và độ dẫn điện của đất; tần số càng cao thì độ phân giải càng sắc nét nhưng độ xuyên sâu giảm đi.

Nguyên nhân gây ra hiện tượng trôi DC và hiệu ứng rung (ringing) trong dữ liệu GPR là gì? Hiện tượng trôi DC xuất hiện do mạch thu khuếch đại làm dịch chuyển điện áp trung bình của tín hiệu so với mức 0. Hiệu ứng rung sinh ra do sóng truyền thẳng cực mạnh giữa anten phát và thu đặt sát nhau, tạo ra các dải phản xạ sọc ngang che lấp tín hiệu địa tầng sâu.

Bộ lọc Dewow có vai trò sống còn như thế nào trong chuỗi xử lý tín hiệu GPR? Bộ lọc Dewow (ứng dụng lọc trung vị median) loại bỏ thành phần tần số rất thấp và hiện tượng bão hòa cảm ứng đầu vào. Việc lọc Dewow giúp đưa dữ liệu về giá trị trung bình bằng 0 chuẩn xác, ngăn ngừa hiện tượng méo dạng phổ trước khi tiến hành các bước giải chập và di trú phức tạp.

Tại sao cần phải áp dụng kỹ thuật di trú (Migration) cho ảnh radar 2D? Trong quá trình quét, một dị vật điểm dưới lòng đất sẽ phản xạ sóng tạo thành một đường cong dạng hyperbol rộng trên ảnh B-scan thô. Kỹ thuật di trú (như Stolt Migration) sẽ di chuyển và gom các năng lượng tán xạ đó về đúng vị trí hình học thực tế, làm sắc nét biên dạng của mục tiêu.

Hệ thống GPR tần số bước SFCW có ưu thế gì vượt trội so với GPR phát xung truyền thống? Hệ thống SFCW tăng dần tần số theo từng bước trên một băng thông cố định, mang lại dải động rộng hơn, khả năng kiểm soát tần số phát tối ưu và tốc độ xử lý nhanh hơn. Công nghệ này cho phép sử dụng bộ lấy mẫu tốc độ vừa phải nhưng vẫn đạt độ phân giải cao và chống nhiễu vượt trội.

Kết luận

  • Luận văn đã hệ thống hóa hoàn chỉnh cơ sở lý thuyết điện từ trường Maxwell, quy luật lan truyền và cơ chế suy hao sóng vô tuyến trong môi trường địa chất phân lớp.
  • Thiết lập và kiểm chứng thành công chuỗi thuật toán xử lý tín hiệu toàn diện từ tiền xử lý (Time Zero, DC Removal, Dewow) đến xử lý nâng cao (AGC, Giải chập Wiener, Di trú Stolt).
  • Nâng cao tỷ số tín hiệu trên nhiễu thêm khoảng 12 dB đến 15 dB và cải thiện độ phân giải hiển thị ảnh cấu trúc ngầm lên hơn 40% so với dữ liệu gốc.
  • Đóng góp giải pháp mô phỏng hoàn chỉnh trên Matlab, mở ra khả năng tự chủ công nghệ phần mềm phân tích dữ liệu địa vật lý tại Việt Nam.
  • Trong thời gian 12 tháng tiếp theo, nhóm nghiên cứu dự kiến sẽ tiếp tục hoàn thiện giải thuật 3D migration và thử nghiệm trên thiết bị phần cứng SFCW băng thông rộng.

Để tìm hiểu sâu hơn về mô hình toán học và các mã nguồn thuật toán xử lý ảnh radar chuyên dụng, bạn đọc hãy tham khảo chi tiết toàn văn luận văn thạc sĩ kỹ thuật của tác giả tại thư viện Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh.