Tổng quan nghiên cứu

Hệ thống radar xuyên đất (Ground Penetrating Radar - GPR) là công nghệ thăm dò địa vật lý không phá hủy, hoạt động dựa trên nguyên lý phát xạ sóng điện từ dải siêu rộng (UWB) từ vài MHz đến 10 GHz vào lòng đất và thu nhận tín hiệu phản xạ từ các mặt tiếp giáp điện môi hoặc dị vật. Tuy nhiên, sóng điện từ khi truyền qua môi trường đất đá không đồng nhất phải chịu sự suy hao năng lượng cực lớn, dao động từ 1 dB/m đến trên 100 dB/m ở tần số 100 MHz và có thể vượt ngưỡng 3000 dB/m ở tần số 1 GHz đối với môi trường đất sét ẩm. Thực tế này dẫn đến hai bài toán kỹ thuật hóc búa: bài toán mục tiêu nằm nông sát mặt đất bị che lấp bởi tín hiệu phản xạ mặt đất (clutter) có biên độ lớn gấp 2 đến 6 lần tín hiệu mục tiêu, và bài toán mục tiêu chôn sâu có tín hiệu phản xạ suy giảm mạnh, chìm sâu trong nền nhiễu với tỉ số tín hiệu trên nhiễu (SNR) xuống mức -10 dB đến -20 dB.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Điện tử (mã số 605270) do học viên Lê Thị Thu Thủy thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia TP.HCM dưới sự hướng dẫn của TS. Đỗ Hồng Tuấn, tập trung giải quyết bài toán xác định khoảng cách đến vật thể có tín hiệu phản xạ yếu trong hệ thống GPR sử dụng mô hình sóng liên tục điều tần (FMCW). Mục tiêu cốt lõi của công trình là nghiên cứu, mô phỏng và đánh giá so sánh chuyên sâu ba phương pháp ước lượng phổ tiêu biểu: phương pháp phi tham số Periodogram, phương pháp không gian con MUSIC (Multiple Signal Classification) và phương pháp kết hợp cải tiến WMUSIC (Windowed-MUSIC). Nghiên cứu hướng đến việc tối ưu hóa độ chính xác định vị mục tiêu với sai số khoảng cách dưới 0,05 m và giảm thiểu tối đa tỷ lệ báo động giả trong điều kiện can nhiễu phức tạp, tạo tiền đề khoa học cho việc chế tạo thiết bị GPR nội địa hóa giá thành thấp 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

Nghiên cứu được xây dựng dựa trên hệ thống lý thuyết trường điện từ kinh điển và lý thuyết xử lý tín hiệu thống kê nâng cao:

  • Hệ phương trình Maxwell và lý thuyết truyền sóng điện từ: Mô tả tương tác của sóng điện từ trong môi trường vật chất đặc trưng bởi ba thông số: độ dẫn điện $\sigma$, độ từ thẩm $\mu$ và hằng số điện môi $\varepsilon$. Định luật Snell và hệ số phản xạ Fresnel xác định sự biến đổi biên độ, góc tới và góc phản xạ tại các mặt ranh giới phân lớp địa chất.
  • Mô hình hệ thống radar FMCW: Tín hiệu phát đi có tần số biến thiên tuyến tính theo thời gian. Tín hiệu phản xạ từ mục tiêu sau thời gian trễ $T_d$ được trộn với tín hiệu phát để tạo ra tần số phách (beat frequency) $f_b = T_d \cdot (\Delta f / t_p)$, từ đó xác định khoảng cách chính xác đến mục tiêu qua công thức $R = \frac{c \cdot f_b}{2 \cdot \Delta f \cdot \sqrt{\varepsilon_r}}$, trong đó $c$ là vận tốc ánh sáng và $\varepsilon_r$ là hằng số điện môi tương đối của đất.
  • Định lý Wiener-Khinchine và phân tích phổ không gian con: Mật độ phổ công suất (PSD) của tín hiệu dừng nghĩa rộng (WSS) là biến đổi Fourier của hàm tự tương quan. Mô hình phân rã ma trận tự tương quan $R_x$ kích thước $M \times M$ thành không gian con tín hiệu (chứa $p$ vector riêng ứng với các trị riêng lớn nhất) và không gian con nhiễu (chứa $M-p$ vector riêng ứng với phương sai nhiễu $\sigma_w^2$).
  • Hệ thống các khái niệm then chốt: Periodogram (ước lượng phổ qua biến đổi Fourier nhanh FFT), MUSIC (ước lượng tần số siêu phân giải dựa trên tính trực giao giữa vector dẫn và không gian con nhiễu), WMUSIC (kỹ thuật kết hợp phân đoạn cửa sổ hóa FFT và hàm phổ MUSIC chuẩn hóa), tỷ số tín hiệu trên can nhiễu và nhiễu (SINR), và diện tích phản xạ radar hiệu dụng (RCS).

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

  • Nguồn dữ liệu và môi trường thực nghiệm: Dữ liệu nghiên cứu được xây dựng thông qua mô phỏng toán học chính xác trên phần mềm MATLAB. Mỗi kịch bản khảo sát được kiểm chứng qua 1000 lần lặp độc lập theo phương pháp Monte Carlo để đảm bảo tính hội tụ và loại bỏ triệt để các sai số ngẫu nhiên của nhiễu Gauss trắng.
  • Cỡ mẫu và phương pháp chọn mẫu: Thuật toán xử lý trên chuỗi tín hiệu phách rời rạc theo thời gian thu được sau bộ lọc thông thấp (LPF). Cỡ mẫu tín hiệu $N$ được khảo sát toàn diện qua các mức 32, 64, 128, 256, 512, 1024 và 1500 mẫu. Độ dài biến đổi Fourier $N_{FFT}$ được phân tích từ 128 đến 4096 điểm. Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu đều với chu kỳ cố định, chia nhỏ chuỗi $N$ mẫu thành $K$ phân đoạn có độ dài $L$, dịch chuyển nhau một khoảng offset $D = L/2$ (tương ứng độ trùng lặp 50%).
  • Lý do lựa chọn phương pháp phân tích: Periodogram được chọn làm thuật toán mốc (baseline) đại diện cho nhóm phi tham số nhờ ưu điểm tính toán nhanh nhưng bị giới hạn độ phân giải bởi hiệu ứng búp sóng phụ. Thuật toán MUSIC được chọn để đạt độ phân giải siêu cao nhằm tách hai tần số phách rất gần nhau. Thuật toán WMUSIC được lựa chọn làm giải pháp trọng tâm nhằm kết hợp ưu thế độ sắc nét của MUSIC và khả năng triệt tiêu đỉnh phổ giả của FFT cửa sổ hóa trong điều kiện SNR thấp từ -20 dB đến 30 dB. Toàn bộ tiến trình nghiên cứu và mô phỏng được hoàn thành trong giai đoạn 2012–2014.

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 thực nghiệm 1000 lần lặp cho thấy những kết quả định lượng rõ nét:

  • Khả năng phân tách tín hiệu mục tiêu gần mặt đất (Bài toán 1): Khi tỉ số biên độ can nhiễu mặt đất trên tín hiệu mục tiêu $A_1/A_2$ tăng từ 2 lên 6 tại mức SNR = 0 dB và 20 dB, phương pháp Periodogram hoàn toàn mất khả năng phân tách do búp phổ phụ của tín hiệu mặt đất che khuất mục tiêu. Ngược lại, thuật toán MUSIC và WMUSIC phân tách xuất sắc hai đỉnh phổ cách biệt, duy trì sai số định vị tần số mục tiêu $f_2$ với độ lệch chuẩn cực thấp.
  • Khắc phục triệt để hiện tượng đỉnh phổ giả khi thay đổi bậc ước lượng $p$: Khi khảo sát bậc mô hình $p$ tăng từ 2 đến 12 ở mức SNR = 10 dB và SNR = -10 dB, thuật toán MUSIC thuần túy xuất hiện từ 2 đến 5 đỉnh phổ giả có biên độ nhọn tương đương tín hiệu thật, gây sai lệch nghiêm trọng trong nhận diện mục tiêu. Trong khi đó, thuật toán WMUSIC triệt tiêu hoàn toàn 100% các đỉnh phổ giả này nhờ phép nhân tích chập chuẩn hóa giữa phổ FFT và phổ MUSIC.
  • Độ chính xác và sai số bình phương trung bình (MSE) ở mức SNR thấp (Bài toán 2): Trong điều kiện mục tiêu chôn sâu với SNR suy giảm từ 5 dB xuống -20 dB, sai số MSE của phương pháp WMUSIC luôn thấp hơn Periodogram và MUSIC từ 10 đến 100 lần. Tại mốc SNR = -10 dB, $N = 512$ và $N_{FFT} = 1024$, giá trị MSE của WMUSIC đạt mức xấp xỉ $10^{-4}$, trong khi MUSIC và Periodogram có xu hướng phân kỳ và dao động mạnh.
  • Hiệu quả của việc gia tăng số mẫu $N$ và số điểm $N_{FFT}$: Khi nâng số mẫu $N$ từ 128 lên 512 mẫu, sai số MSE của cả ba phương pháp giảm hơn 60%, giúp chuyển trạng thái từ không phân tách được sang phân tách rõ ràng. Việc tăng $N_{FFT}$ từ 1024 lên 4096 điểm giúp làm mịn đường cong biểu diễn nhưng không thay đổi bản chất giới hạn phân giải vật lý của thuật toán.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp WMUSIC vượt trội là nhờ cơ chế nhân mật độ phổ công suất đã chuẩn hóa của phân đoạn Fourier $P_n{F}$ với phổ nghịch đảo không gian con $P_{su}{M}$. Các thành phần đỉnh phổ giả sinh ra do nhiễu trong không gian con MUSIC không có sự tương đồng về pha và năng lượng trong biến đổi Fourier của các phân đoạn $K$, do đó biên độ tại các tần số giả bị suy giảm triệt để về 0 khi thực hiện phép cộng trên thang đo decibel (dB).

So sánh với các công bố học thuật quốc tế của Jiang W. và cộng sự, nghiên cứu này đã kiểm chứng thành công tính khả thi của thuật toán trong môi trường mô phỏng GPR đặc thù với tín hiệu phản xạ suy hao phi tuyến. Dữ liệu thực nghiệm của luận văn được minh họa sinh động qua hai dạng biểu diễn trực quan:

  • Hệ thống đồ thị đường biểu diễn: Đồ thị quan hệ giữa sai số MSE và tỉ số SNR trải dài từ -20 dB đến 30 dB thể hiện rõ điểm uốn ngưỡng suy giảm hiệu năng của MUSIC thuần túy dưới mốc -5 dB, đồng thời khẳng định độ ổn định của WMUSIC.
  • Biểu đồ mật độ phổ công suất chuẩn hóa PSD (dB): So sánh trực tiếp giữa 3 thuật toán trên cùng một hệ tọa độ tần số chuẩn hóa, cho thấy búp sóng chính của WMUSIC hẹp và sắc nét tương đương MUSIC nhưng hoàn toàn không tồn tại búp sóng phụ lớn như Periodogram (vốn có biên độ búp phụ chỉ thấp hơn búp chính khoảng 13 dB).

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

Để đưa kết quả nghiên cứu vào ứng dụng thực tiễn trong ngành kỹ thuật radar và địa vật lý tại Việt Nam, bốn giải pháp cụ thể được đề xuất:

  • Tích hợp thuật toán WMUSIC vào chip xử lý tín hiệu số DSP/FPGA: Thiết kế kiến trúc tính toán song song ma trận tự tương quan kích thước $M \times M$ và phân rã trị riêng tối ưu hóa, đảm bảo thời gian xử lý dưới 50 ms cho mỗi vết quét (trace) radar. Mục tiêu đặt ra là hoàn thiện khối xử lý nhúng trong vòng 12 tháng, do các nhóm nghiên cứu phần cứng tại các viện và trường đại học kỹ thuật chủ trì.
  • Áp dụng kỹ thuật lọc thích nghi triệt tiêu can nhiễu phản xạ mặt đất: Xây dựng khối tiền xử lý tín hiệu kết hợp giữa phép trừ trung bình nền (background subtraction) và bộ lọc thích nghi LMS trước khi đưa tín hiệu phách vào khối ước lượng phổ, nhằm nâng cao tỉ số SINR thêm tối thiểu 6 dB đến 10 dB. Lộ trình triển khai dự kiến kéo dài 6 tháng do các kỹ sư phát triển thuật toán đảm nhiệm.
  • Xây dựng bảng tra cứu thực nghiệm hằng số điện môi cho các nhóm thổ nhưỡng Việt Nam: Tổ chức đo đạc thực địa và chuẩn hóa thông số hằng số điện môi $\varepsilon_r$ (từ 3 đối với cát khô đến trên 25 đối với đất phù sa ngập nước) tại 5 vùng địa chất trọng điểm. Dự án do các đơn vị trắc địa và địa vật lý phối hợp thực hiện trong 18 tháng nhằm khống chế sai số ước lượng độ sâu thực tế dưới 2%.
  • Nghiên cứu kết hợp kỹ thuật đa bước tần SFCW và thuật toán WMUSIC: Mở rộng mô hình xử lý cho hệ thống radar sóng liên tục bước nhảy tần (SFCW) đối với các mục tiêu chôn sâu trên 5 mét, khai thác lợi thế băng thông hẹp của từng bước nhảy để duy trì SNR cao. Đề tài cần được phát triển trong khuôn khổ nghiên cứu sau đại học với thời gian 24 thá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à thực tiễn sâu sắc, đặc biệt hữu ích cho 4 nhóm đối tượng:

  • Kỹ sư R&D phần cứng và phần mềm Radar: Nắm vững phương pháp thiết kế khối xử lý tín hiệu trung tần cho radar FMCW/SFCW, trực tiếp ứng dụng thuật toán WMUSIC để nâng cao độ nhạy phát hiện mục tiêu có bề mặt phản xạ RCS nhỏ dưới 0,01 m².
  • Nhà khoa học và nghiên cứu sinh ngành Xử lý tín hiệu - Điện tử viễn thông: Sử dụng tài liệu làm cơ sở tham chiếu chuyên sâu về các thuật toán ước lượng phổ không gian con, kỹ thuật phân rã ma trận tự tương quan và các phương pháp triệt tiêu đỉnh phổ giả.
  • Đơn vị khảo sát địa kỹ thuật, công trình ngầm và rà phá bom mìn: Tiếp cận giải pháp công nghệ nâng cao độ chính xác khi định vị đường ống kim loại, cáp viễn thông chôn sâu dưới 3 mét hoặc vật thể chưa nổ trong môi trường đất có độ dẫn điện cao.
  • Giảng viên và học viên cao học tại các trường đại học kỹ thuật: Bổ sung tài liệu tham khảo giảng dạy chất lượng cao cho các môn học chuyên đề như Xử lý tín hiệu số nâng cao, Lý thuyết trường điện từ ứng dụng và Kỹ thuật Radar hiện đại.

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

Tại sao phương pháp Periodogram lại thất bại khi xác định mục tiêu ở gần mặt đất?
Phương pháp Periodogram sử dụng biến đổi Fourier trực tiếp nên độ phân giải bị giới hạn bởi nghịch đảo độ dài lấy mẫu. Khi mục tiêu ở gần mặt đất, tín hiệu phản xạ từ mặt đất có năng lượng vượt trội gấp 2 đến 6 lần tín hiệu mục tiêu, tạo ra búp phổ phụ lớn che lấp hoàn toàn đỉnh phổ của mục tiêu, khiến sai số định vị khoảng cách tăng trên 50%.

Hạn chế lớn nhất của thuật toán siêu phân giải MUSIC thuần túy là gì?
Hạn chế nghiêm trọng nhất của MUSIC là sự xuất hiện của các đỉnh phổ giả khi ước lượng thừa bậc mô hình $p$ (từ 2 lên 12), dẫn đến báo động giả về sự tồn tại của mục tiêu. Ngoài ra, khi tỉ số SNR giảm xuống dưới 0 dB, khả năng phân tách không gian con tín hiệu và nhiễu bị suy giảm mạnh, làm tăng đột biến sai số MSE.

Thuật toán WMUSIC hoạt động như thế nào để loại bỏ hoàn toàn các đỉnh phổ giả?
WMUSIC chia chuỗi dữ liệu $N$ mẫu thành $K$ phân đoạn có độ dài $L$ với độ trùng lặp $D = L/2$, áp dụng hàm cửa sổ và tính toán mật độ phổ công suất chuẩn hóa của cả FFT và MUSIC. Bằng cách nhân hai hàm phổ này trên thang decibel, các đỉnh phổ giả của MUSIC vốn không có sự tích lũy năng lượng nhất quán trên phổ FFT sẽ bị triệt tiêu hoàn toàn về mức 0.

Độ sâu phát hiện hiệu quả của hệ thống FMCW GPR trong nghiên cứu này đạt bao nhiêu?
Khoảng cách thăm dò tối ưu của hệ thống FMCW GPR mô phỏng đạt từ 0,2 mét đến khoảng 3 mét đối với môi trường đất thông thường. Do đất có độ suy hao cao từ 10 dB/m đến trên 100 dB/m, đối với các mục tiêu chôn sâu hơn 5 mét, hệ thống cần chuyển sang sử dụng mô hình SFCW hoặc hạ thấp tần số sóng mang trung tâm.

Tăng số mẫu $N$ và số điểm $N_{FFT}$ mang lại sự khác biệt như thế nào đối với chất lượng ước lượng phổ?
Tăng số mẫu $N$ từ 128 lên 512 mẫu giúp tăng lượng thông tin thống kê thời gian, cải thiện trực tiếp khả năng phân giải hai mục tiêu sát nhau và giảm sai số MSE hơn 60%. Trong khi đó, tăng $N_{FFT}$ từ 256 lên 1024 điểm chỉ làm tăng số điểm nội suy tần số giúp đồ thị mượt mà hơn, nhưng không giúp tách được hai mục tiêu nếu số mẫu $N$ ban đầu quá nhỏ.

Kết luận

Công trình luận văn thạc sĩ đã giải quyết trọn vẹn và tường minh bài toán định vị khoảng cách mục tiêu có tín hiệu phản xạ yếu trong hệ thống radar xuyên đất FMCW GPR:

  • Hệ thống hóa toàn diện cơ sở lý thuyết truyền sóng điện từ trong môi trường suy hao phi tuyến và mô hình toán học của tín hiệu phách radar FMCW.
  • Phân tích và chỉ rõ các hạn chế cố hữu của phương pháp Periodogram trong bài toán can nhiễu mặt đất mạnh và phương pháp MUSIC trong môi trường có tỉ số SNR thấp.
  • Chứng minh thực nghiệm mô phỏng tính ưu việt vượt trội của thuật toán WMUSIC với khả năng triệt tiêu 100% đỉnh phổ giả và duy trì sai số MSE ở mức thấp xấp xỉ $10^{-4}$ ngay cả khi SNR chạm ngưỡng -10 dB.
  • Xác lập bộ thông số kỹ thuật tối ưu cho thuật toán gồm cỡ mẫu $N = 512$, độ phân giải $N_{FFT} = 1024$ và tỷ lệ trùng lặp phân đoạn 50%, đảm bảo sự cân bằng giữa độ chính xác và chi phí tính toán.
  • Đóng góp luận cứ khoa học vững chắc phục vụ lộ trình nội địa hóa và chế tạo thiết bị radar xuyên đất thương mại giá thành hợp lý tại Việt Nam.

Trong giai đoạn 12 đến 24 tháng tới, định hướng phát triển then chốt là tích hợp thuật toán WMUSIC lên phần cứng nhúng DSP/FPGA và triển khai đo đạc kiểm chứng tại hiện trường. Độc giả và các nhà nghiên cứu quan tâm có thể tiếp cận toàn văn luận văn thạc sĩ để khai thác chi tiết hệ thống công thức giải tích và mã nguồn mô phỏng MATLAB chuyên sâu.