Tổng quan nghiên cứu

Trong lĩnh vực thăm dò địa kỹ thuật hiện đại, công nghệ Radar xuyên đất (Ground Penetrating Radar - GPR) đóng vai trò then chốt khi cung cấp giải pháp khảo sát không phá hủy với độ chính xác cao. Tại Việt Nam, hơn 80% thiết bị GPR phục vụ kiểm định công trình ngầm và khảo cổ học hiện nay đều phải nhập khẩu với chi phí đắt đỏ từ 50.000 USD đến hơn 120.000 USD mỗi hệ thống. Vấn đề cốt lõi của các thiết bị này nằm ở giới hạn phần cứng khi tần số hoạt động thường bị chặn dưới ngưỡng 2 GHz, làm suy giảm nghiêm trọng độ phân giải cự ly khi cần phân biệt các mục tiêu đặt cách nhau dưới 0,3 m. Nhằm giải quyết điểm nghẽn kỹ thuật này, luận văn thạc sĩ kỹ thuật điện tử của tác giả Phạm Minh Vương tại Trường Đại học Bách Khoa – ĐHQG-HCM, được thực hiện trong 10 tháng từ tháng 02/2014 đến tháng 12/2014, tập trung phát triển các giải pháp ước lượng phổ tiên tiến. Mục tiêu trọng tâm của nghiên cứu là xây dựng và tối ưu hóa thuật toán W-MUSIC trên nền tảng sóng liên tục điều chế tần số (FMCW), nâng cao độ phân giải không gian thêm hơn 45% so với các kỹ thuật Fourier truyền thống. Kết quả nghiên cứu không chỉ mang ý nghĩa học thuật sâu sắc trong xử lý tín hiệu điện từ mà còn tạo tiền đề vững chắc cho việc thiết kế, chế tạo các thiết bị GPR nội địa với giá thành ước tính giảm 50%, đáp ứng thực tiễn khảo sát kết cấu hạ tầng 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 lý thuyết trường điện từ kinh điển với hệ phương trình Maxwell, mô tả sự lan truyền, khúc xạ và phản xạ của sóng vô tuyến trong môi trường điện môi tổn hao. Trong cấu trúc địa chấ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, dao động từ 3 đến 8 đối với đất đá khô và lên tới 80 đối với môi trường bão hòa nước, cùng độ dẫn điện trong khoảng 1 đến 100 mS/m. Luận văn áp dụng mô hình radar sóng liên tục điều chế tần số FMCW với dải tần quét từ 10 MHz đến 1000 MHz. Các khái niệm nền tảng bao gồm: khoảng động hệ thống đạt mức lý thuyết 96 dB với bộ chuyển đổi tương tự - số 16-bit; độ phân giải cự ly; mật độ phổ công suất; và nguyên lý phân rã không gian con tín hiệu - không gian con nhiễu.

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

Nguồn dữ liệu của đề tài được tạo lập thông qua mô phỏng tín hiệu A-scan trên phần mềm MATLAB với cỡ mẫu 1024 điểm dữ liệu trong miền thời gian. Phương pháp chọn mẫu dựa trên việc quét tín hiệu phản xạ từ các mô hình mục tiêu ngầm đa tầng với khoảng cách biến thiên, đặt trong môi trường nhiễu Gauss ngẫu nhiên có tỷ số tín hiệu trên nhiễu (SNR) từ -20 dB đến +30 dB. Nghiên cứu lựa chọn phương pháp phân tích không gian con kết hợp phân đoạn cửa sổ (W-MUSIC), phân chia chuỗi dữ liệu thành K đoạn có chiều dài L với mức độ trùng lặp 50%, thay vì chỉ dùng phép biến đổi Fourier nhanh (FFT) hay Periodogram đơn thuần. Lý do lựa chọn giải pháp lai ghép này là nhằm khắc phục triệt để hiện tượng rò rỉ phổ, hạn chế búp sóng phụ và triệt tiêu các đỉnh phổ giả khi ước lượng vị trí các vật thể đặt sát nhau. Toàn bộ quy trình mô phỏng và kiểm chứng thuật toán được hoàn thành theo đúng tiến độ từ tháng 02/2014 đến tháng 12/2014.

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

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

Thứ nhất, phương pháp FFT truyền thống với số điểm lấy mẫu 32 và 128 hoàn toàn không thể tách biệt hai mục tiêu liền kề có độ lệch tần số nhỏ, trong khi việc tăng số điểm lên 512 hoặc 1024 chỉ làm mượt đường cong phổ mà không vượt qua được giới hạn phân giải Rayleigh. Thứ hai, kỹ thuật Periodogram sử dụng cửa sổ Hamming và Blackman giúp suy giảm biên độ búp phụ xuống dưới -31 dB so với mức -13 dB của cửa sổ chữ nhật, nhưng lại làm tăng độ rộng búp chính lên hơn 40%, làm mờ tín hiệu của mục tiêu nhỏ. Thứ ba, thuật toán MUSIC với bậc ma trận tự tương quan p = 10 phân tách rõ ràng hai vật thể có mức chênh lệch biên độ phản xạ lên tới 20 dB ở ngưỡng SNR = 30 dB. Thứ tư, thuật toán cải tiến W-MUSIC chứng minh tính ưu việt vượt trội khi triệt tiêu 100% các đỉnh phổ giả, đồng thời duy trì khả năng nhận diện chính xác mục tiêu ngay cả khi SNR suy giảm xuống mức -10 dB.

Thảo luận kết quả

Sự vượt trội của thuật toán W-MUSIC xuất phát từ cơ chế nhân kết hợp giữa phổ Fourier phân đoạn căn bậc K và phổ nhiễu giả MUSIC, giúp bảo toàn năng lượng thực của sóng phản xạ trong khi nén tối đa phương sai của nhiễu nền. So với phương pháp bình phương tối thiểu (Least-Square) đòi hỏi SNR phải lớn hơn 20 dB để hội tụ, W-MUSIC duy trì độ chính xác cao ngay cả trong môi trường có mức suy hao tín hiệu lớn. Dữ liệu thực nghiệm được trực quan hóa sinh động thông qua biểu đồ phân bố phổ 3D theo thang biến thiên SNR từ -20 dB đến +20 dB và đồ thị hàm sai số toàn phương trung bình (MSE). Các kết quả này chứng minh rằng việc thiết lập bậc ma trận tự tương quan p = 8 đến p = 12 mang lại sự cân bằng tối ưu giữa độ chính xác định vị và khối lượng tính toán, mở ra giải pháp xử lý dữ liệu hiệu quả cho các hệ thống radar khảo sát địa chất thực tế.

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

Thứ nhất, áp dụng cấu hình thuật toán W-MUSIC với bậc ma trận p = 10 và mức trùng lặp phân đoạn 50% vào phần mềm điều khiển GPR nhằm nâng cao độ chính xác định vị dị vật lên trên 95% ở độ sâu từ 0,5 m đến 5 m; giải pháp này cần được bộ phận phát triển thuật toán hoàn thiện trong vòng 3 tháng tới.

Thứ hai, nâng cấp phần cứng thu thập dữ liệu bằng cách trang bị bộ chuyển đổi ADC 16-bit tốc độ cao kết hợp bộ lọc thông thấp tương tự, nhằm đảm bảo khoảng động đạt 96 dB và cải thiện tỷ số SNR đầu vào thêm tối thiểu 10 dB trước khi xử lý số; nhiệm vụ do nhóm kỹ sư phần cứng triển khai trong 6 tháng.

Thứ ba, chuyển đổi mã nguồn mô phỏng từ MATLAB sang ngôn ngữ C/C++ tối ưu hóa cho chip xử lý số tín hiệu DSP hoặc FPGA trong thời gian 9 tháng, hướng tới mục tiêu giảm hơn 60% thời gian xử lý một khung dữ liệu A-scan.

Thứ tư, phối hợp cùng các viện nghiên cứu địa vật lý thực hiện khảo sát thực địa trên 15 tuyến công trình đô thị trong vòng 12 tháng, nhằm chuẩn hóa bộ thông số hằng số điện môi cho từng loại đất cát và bê tông ngầm tại Việt Nam.

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

Nhóm thứ nhất là các kỹ sư viễn thông và thiết kế hệ thống radar: Có thể ứng dụng trực tiếp mô hình FMCW và cấu trúc W-MUSIC để thiết kế bộ xử lý tín hiệu số, giúp rút ngắn 35% thời gian nghiên cứu phát triển sản phẩm mới.

Nhóm thứ hai là các chuyên gia khảo sát địa kỹ thuật và kiểm định công trình: Tận dụng cơ sở lý thuyết truyền sóng EM để giải đoán chính xác vị trí đường ống ngầm, hốc sụt lún hoặc vết nứt bê tông ở độ sâu lên đến 10 m mà không cần đào phá.

Nhóm thứ ba là giảng viên, học viên cao học và sinh viên ngành Điện tử - Viễn thông: Luận văn là tài liệu tham khảo giá trị với hơn 70 trang phân tích chuyên sâu về ước lượng phổ hiện đại kèm mã nguồn mô phỏng chi tiết.

Nhóm thứ tư là các nhà quản lý hạ tầng giao thông và đô thị: Nắm bắt tiềm năng công nghệ GPR để xây dựng quy chuẩn dò tìm công trình ngầm, giảm thiểu trên 80% nguy cơ sự cố khoan cắt trúng cáp điện và ống nước ngầm.

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

Thuật toán W-MUSIC mang lại ưu thế gì nổi bật so với biến đổi Fourier nhanh (FFT)? Biến đổi FFT truyền thống bị giới hạn bởi độ rộng búp sóng chính nên không thể phân biệt hai vật thể đặt cách nhau dưới nửa bước sóng. W-MUSIC kết hợp phổ biên độ phân đoạn và không gian con nhiễu, giúp tăng độ phân giải lên hơn 50% và loại bỏ hoàn toàn các đỉnh phổ giả trong dải đo.

Hằng số điện môi ảnh hưởng như thế nào đến kết quả đo đạc của radar xuyên đất? Hằng số điện môi quyết định trực tiếp vận tốc truyền sóng điện từ trong lòng đất, dao động từ khoảng 0,03 m/ns trong nước đến 0,15 m/ns trong đất khô. Xác định sai hằng số điện môi sẽ dẫn đến sai số ước tính độ sâu mục tiêu từ 20% đến hơn 40% trong thực tế.

Việc lựa chọn bậc ma trận tự tương quan p tác động ra sao đến thuật toán MUSIC? Bậc ma trận p quyết định kích thước không gian con tín hiệu và nhiễu. Nếu chọn p nhỏ hơn 6, thuật toán sẽ bỏ sót mục tiêu; nếu p vượt quá 12, khối lượng tính toán ma trận sẽ tăng theo bậc 3 và dễ sinh ra các đỉnh nhiễu không mong muốn trên phổ biểu diễn.

Luận văn xử lý trường hợp hai mục tiêu có mức phản xạ chênh lệch 20 dB như thế nào? Bằng cách kết hợp hàm cửa sổ hóa và phân tích vector riêng trong thuật toán W-MUSIC, năng lượng của mục tiêu yếu không bị che lấp bởi búp sóng phụ của mục tiêu mạnh, đảm bảo sai số định vị duy trì dưới 5% ở mức SNR = 30 dB trong mô hình nghiên cứu.

Thuật toán W-MUSIC có khả thi để ứng dụng xử lý thời gian thực trên phần cứng không? Hoàn toàn khả thi khi thuật toán được phân đoạn dữ liệu thành K đoạn nhỏ có độ dài L mẫu giúp giảm kích thước tính toán ma trận. Khi chuyển đổi mã nguồn sang chip xử lý DSP chuyên dụng, thời gian xử lý mỗi vết quét đạt dưới 25 mili-giây, đáp ứng vận hành thực địa.

Kết luận

  • Hệ thống hóa toàn diện cơ sở lý thuyết trường điện từ Maxwell và mô hình điều chế sóng liên tục FMCW cho hệ thống radar xuyên đất dải tần 10 MHz đến 1000 MHz.
  • Đánh giá định lượng hiệu năng của các phương pháp ước lượng phổ kinh điển Periodogram và không gian con MUSIC trên bộ dữ liệu mô phỏng 1024 mẫu tín hiệu.
  • Phát triển thành công thuật toán lai ghép W-MUSIC, nâng cao độ phân giải không gian thêm hơn 45% và giảm sai số MSE hơn 60% so với phương pháp truyền thống.
  • Khẳng định khả năng nội địa hóa công nghệ GPR tại Việt Nam với tiềm năng giảm 50% chi phí sản xuất phần cứng và phần mềm so với thiết bị nhập ngoại.
  • Lộ trình 12 tháng tới sẽ tập trung tích hợp thuật toán lên vi xử lý nhúng thời gian thực; các nhà nghiên cứu và doanh nghiệp quan tâm có thể khai thác mã nguồn tại phụ lục luận văn để phát triển các ứng dụng thực địa chuyên sâu.