Luận Văn Thạc Sĩ Về Kỹ Thuật Ước Lượng Phổ Tín Hiệu GPR Trong Mô Hình Các Vật Thể Gần Nhau

Chuyên đề nghiên cứu Kỹ Thuật Ước Lượng Phổ Tín Hiệu GPR Trong Mô Hình Vật ..., cập nhật xu hướng mới, giá trị tham khảo cao cho chuyên gia công nghệ

• Tác giả: Nguyễn Minh Linh
• Người hướng dẫn: TS. Đỗ Hồng Tuấn
• Trường học: Đại học Quốc gia TP. HCM
• Chuyên ngành: Kỹ thuật điện tử
• Đề tài: Kỹ Thuật Ước Lượng Phổ Tín Hiệu GPR Trong Mô Hình Các Vật Thể Đặt Gần Nhau
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• MỤC LỤC HÌNH ẢNH
• CÁC TỪ VIẾT TẮT
• LỜI MỞ ĐẦU
• MỤC ĐÍCH ĐỀ TÀI
• Các vấn đề quan tâm
• Các hướng nghiên cứu hiện nay
• 1. Chương 1: TỔNG QUAN VỀ RADAR XUYÊN ĐẤT
• 1.1. Lịch sử phát triển
• 1.2. Giới thiệu chung về radar xuyên đất
• 1.3. Nguyên lý hoạt động
• 1.4. Các định nghĩa và các đặc điểm của tín hiệu radar xuyên đất
• 1.5. Tiêu chuẩn thiết kế cho GPR
• 2. Chương 2: CƠ SỞ LÝ THUYẾT CỦA RADAR XUYÊN ĐẤT
• 2.1. Sóng điện từ - Phương trình Maxwell
• 2.2. Tính chất của vật liệu
• 2.3. Sự truyền sóng của trường điện từ
• 2.4. Tính chất của sóng
• 2.5. Phản xạ, tán xạ và truyền dẫn tại bề mặt
• 2.6. Độ phân giải và vùng tác động
• 2.7. Sự suy hao tán xạ
• 2.8. Vận tốc truyền
• 3. Chương 3: CÁC PHƯƠNG PHÁP ƯỚC LƯỢNG PHỔ TÍN HIỆU
• 3.1. Lịch sử phát triển của các phương pháp ước lượng phổ
• 3.2. Mô hình FWCW sử dụng trong kỹ thuật ước lượng phổ tín hiệu
• 3.3. Các phương pháp ước lượng phổ không tham số (Nonparametric Methods)
• 3.3.1. Phương pháp Periodogram
• 3.3.2. Phương pháp cửa sổ (Windowing Method)
• 3.4. Phương pháp ước lượng phổ sử dụng không gian con tín hiệu
• 3.4.1. Pisarenko Harmonic Decomposition (PHD)
• 3.4.3. Phương pháp véc tơ riêng (Eigenvector Method)
• 3.5. Phương pháp W-MUSIC
• 3.5.1. Mô hình hóa thuật toán W-MUSIC
• 3.5.2. Chi tiết thuật toán W-MUSIC
• 4. Chương 4: SỬ DỤNG CÁC PHƯƠNG PHÁP ƯỚC LƯỢNG PHỔ TRONG MÔ HÌNH GPR
• 4.1. Mô hình bài toán
• 4.2. Sử dụng phương pháp ước lượng phổ Periodogram
• 4.2.1. Kết hợp thuật toán FFT
• 4.2.2. Phương pháp Periodogram sử dụng các hàm cửa sổ
• 4.3. Phương pháp ước lượng phổ sử dụng không gian con tín hiệu
• 4.3.1. Phương pháp PHD, Eigen Vector, MUSIC và Minimum-Norm
• 4.3.2. Superimposed FFT và Superimposed MUSIC
• 4.4. Phương pháp W-MUSIC
• 5. Chương 5: ĐÁNH GIÁ KẾT QUẢ ĐỀ TÀI
• 5.1. Kết quả đạt được
• 5.2. Giới hạn của luận văn
• 5.3. Hướng phát triển của đề tài
• TÀI LIỆU THAM KHẢO
• PHỤ LỤC: MATLAB CODE

Tóm tắt

I. Tổng quan về GPR

Chương này giới thiệu tổng quan về Radar xuyên đất GPR, bao gồm lịch sử phát triển và nguyên lý hoạt động. GPR là một phương pháp địa vật lý quan trọng, cho phép phát hiện các vật thể dưới bề mặt đất mà không cần phải đào bới. Lịch sử phát triển của GPR bắt đầu từ những năm đầu thế kỷ 20, với những nghiên cứu đầu tiên về việc sử dụng sóng điện từ để xác định vị trí của các vật thể chôn lấp. Các kỹ thuật như FMCW và SFCW đã được phát triển để cải thiện khả năng phát hiện và phân tích tín hiệu. Chương này cũng trình bày sơ đồ tổng quát của hệ thống GPR, giúp người đọc hiểu rõ hơn về cấu trúc và chức năng của các thành phần trong hệ thống.

1.1 Lịch sử phát triển

Việc sử dụng tín hiệu điện từ để xác định sự hiện diện của vật thể đã được nghiên cứu từ đầu thế kỷ 20. Nghiên cứu của Hulsmeyer vào năm 1904 đã mở ra hướng đi mới cho việc phát hiện vật thể bằng sóng điện từ. Sau đó, các nghiên cứu của Leimbach và Lowy đã phát triển kỹ thuật này thành một phương pháp có thể xác định vị trí của các vật thể chôn lấp. Kỹ thuật xung được phát triển từ những năm 1930 đã cho thấy khả năng thăm dò độ sâu trong nhiều loại môi trường khác nhau. Sự phát triển của GPR đã mở rộng ra nhiều lĩnh vực ứng dụng, từ khảo cổ học đến xây dựng và địa kỹ thuật.

1.2 Giới thiệu chung về radar xuyên đất

Radar xuyên đất (GPR) là một công nghệ tiên tiến cho phép phát hiện và phân tích các vật thể dưới bề mặt đất. GPR sử dụng sóng điện từ để tạo ra hình ảnh của các lớp đất và vật thể chôn lấp. Phương pháp này có nhiều ưu điểm, bao gồm khả năng không phá hủy và độ chính xác cao trong việc xác định vị trí của các vật thể. GPR được ứng dụng rộng rãi trong nhiều lĩnh vực như khảo sát địa chất, khảo cổ học và xây dựng. Việc phát triển các kỹ thuật mới trong GPR giúp cải thiện độ phân giải và khả năng phát hiện các vật thể nhỏ nằm gần nhau.

II. Cơ sở lý thuyết của GPR

Chương này trình bày cơ sở lý thuyết của Radar xuyên đất GPR, bao gồm lý thuyết về sóng điện từ và sự truyền sóng trong môi trường vật chất. Các yếu tố như phản xạ, khúc xạ và tán xạ của sóng điện từ sẽ ảnh hưởng đến tín hiệu thu được. Độ phân giải và vận tốc truyền sóng trong các môi trường khác nhau cũng được phân tích. Những kiến thức này là nền tảng để hiểu rõ hơn về cách thức hoạt động của GPR và các yếu tố ảnh hưởng đến chất lượng tín hiệu. Việc nắm vững lý thuyết này sẽ giúp cải thiện khả năng xử lý và phân tích tín hiệu trong các ứng dụng thực tế.

2.1 Sóng điện từ Phương trình Maxwell

Sóng điện từ là một phần quan trọng trong lý thuyết của GPR. Phương trình Maxwell mô tả cách thức sóng điện từ lan truyền trong không gian và tương tác với các vật liệu khác nhau. Hiểu rõ về phương trình này giúp các nhà nghiên cứu có thể dự đoán được hành vi của sóng điện từ khi chúng đi qua các lớp đất khác nhau. Điều này rất quan trọng trong việc thiết kế và tối ưu hóa các hệ thống GPR để đạt được hiệu suất tốt nhất trong việc phát hiện vật thể.

2.2 Tính chất của vật liệu

Tính chất của vật liệu ảnh hưởng lớn đến sự truyền sóng điện từ trong GPR. Các yếu tố như độ ẩm, mật độ và thành phần hóa học của đất sẽ quyết định đến tốc độ truyền sóng và mức độ suy giảm tín hiệu. Việc nghiên cứu các tính chất này giúp cải thiện khả năng phát hiện và phân tích tín hiệu trong các ứng dụng thực tế. Các mô hình lý thuyết về sự tương tác giữa sóng điện từ và vật liệu sẽ được trình bày để làm rõ hơn về vấn đề này.

III. Các phương pháp ước lượng phổ tín hiệu

Chương này giới thiệu các phương pháp ước lượng phổ tín hiệu trong GPR, bao gồm lịch sử phát triển và ứng dụng của các phương pháp này. Các phương pháp như Periodogram, MUSIC, và W-MUSIC sẽ được phân tích chi tiết. Mỗi phương pháp có những ưu điểm và nhược điểm riêng, và việc lựa chọn phương pháp phù hợp sẽ ảnh hưởng đến độ chính xác và hiệu quả của việc phát hiện vật thể. Chương này cũng sẽ trình bày mô hình FWCM sử dụng trong quá trình mô phỏng ước lượng, giúp người đọc hiểu rõ hơn về cách thức hoạt động của các phương pháp này.

3.1 Lịch sử phát triển của các phương pháp ước lượng phổ

Lịch sử phát triển của các phương pháp ước lượng phổ trong GPR bắt đầu từ những năm đầu thế kỷ 20. Các nghiên cứu ban đầu tập trung vào việc phát triển các thuật toán đơn giản để xử lý tín hiệu. Theo thời gian, các phương pháp phức tạp hơn như MUSIC và W-MUSIC đã được phát triển, cho phép cải thiện độ phân giải và khả năng phát hiện các vật thể nhỏ. Việc nghiên cứu và phát triển các phương pháp này đã mở ra nhiều cơ hội mới trong việc ứng dụng GPR trong các lĩnh vực khác nhau.

3.2 Mô hình FWCW sử dụng trong kỹ thuật ước lượng phổ tín hiệu

Mô hình FWCW (Frequency Modulated Continuous Wave) là một trong những mô hình quan trọng trong kỹ thuật ước lượng phổ tín hiệu. Mô hình này cho phép thu thập dữ liệu tín hiệu một cách liên tục và hiệu quả. Việc sử dụng mô hình này trong GPR giúp cải thiện khả năng phát hiện và phân tích tín hiệu, đặc biệt là trong các trường hợp có nhiều vật thể chôn lấp gần nhau. Chương này sẽ trình bày chi tiết về cách thức hoạt động của mô hình FWCW và ứng dụng của nó trong các phương pháp ước lượng phổ.

IV. Sử dụng các phương pháp ước lượng phổ trong mô hình GPR

Chương này trình bày kết quả mô phỏng của các phương pháp ước lượng phổ trong mô hình GPR. Các mô hình đề nghị sẽ được áp dụng cho vật thể đặt dưới lòng đất, đồng thời xét ảnh hưởng của nhiễu có phân bố Gauss cho các mô hình này. Việc phân tích kết quả mô phỏng sẽ giúp đánh giá hiệu quả của từng phương pháp trong việc phát hiện và xác định vị trí của các vật thể. Chương này cũng sẽ thảo luận về các yếu tố ảnh hưởng đến độ chính xác của các phương pháp ước lượng và đề xuất các giải pháp cải thiện.

4.1 Mô hình bài toán

Mô hình bài toán trong GPR được xây dựng dựa trên các yếu tố như loại vật thể, độ sâu chôn lấp và môi trường xung quanh. Việc xác định các thông số này là rất quan trọng để đảm bảo tính chính xác của kết quả mô phỏng. Mô hình sẽ được áp dụng cho các trường hợp cụ thể, giúp người đọc hiểu rõ hơn về cách thức hoạt động của các phương pháp ước lượng trong thực tế. Các kết quả thu được từ mô hình sẽ được so sánh với các phương pháp khác để đánh giá hiệu quả.

4.2 Sử dụng phương pháp ước lượng phổ Periodogram

Phương pháp ước lượng phổ Periodogram là một trong những phương pháp cơ bản trong GPR. Phương pháp này cho phép phân tích tín hiệu thu được và xác định vị trí của các vật thể chôn lấp. Việc kết hợp Periodogram với các thuật toán khác như FFT giúp cải thiện độ chính xác và độ phân giải của kết quả. Chương này sẽ trình bày chi tiết về cách thức hoạt động của phương pháp Periodogram và ứng dụng của nó trong các mô hình GPR.

V. Đánh giá kết quả đạt được

Chương này tóm tắt kết quả đạt được sau khi tiến hành mô phỏng các phương pháp ước lượng trong GPR. Việc đánh giá sẽ dựa trên các tiêu chí như độ chính xác, độ phân giải và khả năng phát hiện các vật thể nhỏ. Đồng thời, chương này cũng sẽ nêu lên những hạn chế của luận văn và hướng phát triển trong tương lai. Các kết quả thu được sẽ được so sánh với các nghiên cứu trước đó để đưa ra những nhận định và đề xuất cải tiến cho các phương pháp ước lượng trong GPR.

5.1 Kết quả đạt được

Kết quả đạt được từ các mô phỏng cho thấy rằng các phương pháp ước lượng phổ trong GPR có khả năng phát hiện và xác định vị trí của các vật thể chôn lấp với độ chính xác cao. Việc áp dụng các thuật toán như W-MUSIC đã cho thấy hiệu quả vượt trội trong việc xử lý tín hiệu, đặc biệt là trong các trường hợp có nhiều vật thể gần nhau. Các kết quả này mở ra nhiều cơ hội mới cho việc ứng dụng GPR trong các lĩnh vực khác nhau.

5.2 Hướng phát triển của đề tài

Hướng phát triển của đề tài sẽ tập trung vào việc cải thiện các phương pháp ước lượng phổ trong GPR. Việc nghiên cứu và phát triển các thuật toán mới có thể giúp nâng cao độ chính xác và hiệu quả của việc phát hiện vật thể. Đồng thời, việc áp dụng công nghệ mới trong GPR cũng sẽ được xem xét để mở rộng khả năng ứng dụng của phương pháp này trong thực tế.

Nội dung chính

Tổng quan nghiên cứu

Kỹ thuật Radar xuyên đất (Ground Penetrating Radar - GPR) là phương pháp địa vật lý hiện đại hoạt động trong dải tần số từ 1 MHz đến 1000 MHz, cho phép khảo sát cấu trúc địa chất và công trình ngầm ở độ sâu từ 40 m đến 100 m mà không cần đào bới hay phá hủy cấu trúc bề mặt. Tuy nhiên, trong môi trường thực tế, khi các công trình ngầm như đường ống, cáp ngầm hoặc dị vật địa chất đặt sát nhau với khoảng cách cự ly dưới 9 cm, các hệ thống GPR truyền thống bị giới hạn nghiêm trọng bởi độ phân giải của phần cứng (thường bị giới hạn ở tần số dưới 2 GHz) và hiện tượng giao thoa sóng phản xạ.

Vấn đề cốt lõi đặt ra là sự suy giảm mạnh của tín hiệu khi truyền qua môi trường đất có độ ẩm và độ dẫn điện cao, kết hợp với hiện tượng chồng lấn phổ khi hai mục tiêu nằm quá gần nhau, khiến thuật toán biến đổi Fourier nhanh (FFT) cổ điển không thể phân tách chính xác các đỉnh phổ riêng biệt. Mục tiêu cụ thể của công trình là nghiên cứu các mô hình tín hiệu GPR, xây dựng và đánh giá hiệu năng các thuật toán ước lượng phổ tham số và phi tham số, từ đó đề xuất giải thuật lai ghép tối ưu Windowed MUSIC (W-MUSIC) nhằm đạt độ phân giải siêu cao trong điều kiện các vật thể đặt gần nhau.

Nghiên cứu được triển khai thực nghiệm mô phỏng trên nền tảng tín hiệu FMCW với môi trường đất đồng nhất có hằng số điện môi tương đối bằng 9 và hệ số suy hao 5 dB/m, hoàn thành vào tháng 09 năm 2013 tại Trường Đại học Bách Khoa – ĐHQG TP.HCM. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc nâng cao độ chính xác nhận dạng dị vật ở dải tỷ số tín hiệu trên nhiễu (SNR) từ 10 dB đến 30 dB, giảm thiểu 100% rủi ro đập phá kết cấu trong kiểm định công trình và tạo cơ sở lý thuyết vững chắc để phát triển các thiết bị GPR nội địa giá thành thấp.

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 hai trụ cột lý thuyết chính: lý thuyết trường điện từ Maxwell và lý thuyết ước lượng phổ hiện đại. Quá trình truyền sóng GPR tuân theo hệ phương trình vi phân Maxwell, xác định sự biến thiên của điện trường và từ trường trong môi trường vật chất. Vận tốc truyền sóng điện từ trong môi trường điện môi không từ tính được xác định bởi công thức liên hệ trực tiếp với hằng số điện môi tương đối:

$$v = \frac{c}{\sqrt{\varepsilon_r}}$$

Với $c = 3 \times 10^8\text{ m/s}$, trong môi trường có hằng số điện môi bằng 9, vận tốc sóng giảm còn $10^8\text{ m/s}$, làm thay đổi trực tiếp độ dài bước sóng và độ phân giải cự ly.

Mô hình tín hiệu radar sử dụng kỹ thuật sóng liên tục điều tần (FMCW), trong đó tín hiệu phát có tần số quét tuyến tính từ tần số bắt đầu đến tần số kết thúc trong chu kỳ quét thời gian. Tín hiệu phản xạ từ các mục tiêu ngầm được trộn với tín hiệu phát để tạo ra tín hiệu phách trong miền thời gian, trong đó tần số của tín hiệu phách tỷ lệ thuận với khoảng cách đến mục tiêu.

Bên cạnh đó, khung lý thuyết phân rã không gian con tín hiệu (Subspace-based Spectral Estimation) được ứng dụng triệt để. Ma trận tự tương quan dữ liệu bậc $M \times M$ được phân tích thành hai không gian trực giao: không gian con tín hiệu tạo bởi các vector riêng ứng với các giá trị riêng lớn nhất và không gian con nhiễu tạo bởi các vector riêng ứng với các giá trị riêng nhỏ nhất. Các khái niệm nền tảng bao gồm: dải động tín hiệu (Dynamic Range) đạt xấp xỉ 96 dB với bộ chuyển đổi ADC 16-bit, hệ số phản xạ Fresnel, bề rộng búp sóng chính và mức triệt tiêu búp sóng phụ của các hàm cửa sổ.

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

Nghiên cứu sử dụng phương pháp mô phỏng số lượng lớn trên công cụ chuyên dụng Matlab, lấy mẫu tín hiệu theo định dạng quét A-scan một chiều theo thời gian. Nguồn dữ liệu mô phỏng được tạo lập với chuỗi dữ liệu gồm 150 mẫu quan sát, chu kỳ quét tín hiệu 20 ms, tần số lấy mẫu đạt 7.5 kHz thỏa mãn hoàn toàn điều kiện định lý lấy mẫu Nyquist. Nhiễu môi trường được mô hình hóa dưới dạng nhiễu trắng Gauss cộng (AWGN) với các mức SNR thay đổi liên tục từ -20 dB đến +30 dB.

Phương pháp chọn mẫu áp dụng kỹ thuật chia phân đoạn tín hiệu có cửa sổ trượt (Windowing Segmentation) với $K$ phân đoạn, mỗi phân đoạn có chiều dài $L$ mẫu và độ dịch chuyển $D = L/2$, tạo ra độ trùng lặp tín hiệu chính xác 50% nhằm bảo toàn năng lượng và đặc tính thống kê của chuỗi dữ liệu.

Lý do lựa chọn phương pháp phân tích: Luận văn tiến hành phân tích so sánh đối chuẩn giữa ba nhóm phương pháp:

• Nhóm phi tham số: Chu kỳ đồ Periodogram kết hợp các hàm cửa sổ Chữ nhật, Hamming, Blackman với số điểm FFT từ 32, 128, 512 đến 1024 điểm.
• Nhóm không gian con: Pisarenko Harmonic Decomposition (PHD), Eigenvector (EV), Minimum-Norm và Multiple Signal Classification (MUSIC).
• Nhóm thuật toán lai: Thuật toán W-MUSIC cải tiến kết hợp biến đổi Fourier phân đoạn và hàm phổ giả MUSIC chuẩn hóa.

Việc lựa chọn thuật toán lai W-MUSIC xuất phát từ hạn chế cố hữu của FFT trong việc phân biệt các búp sóng quá gần nhau và nhược điểm sinh đỉnh phổ giả của MUSIC ở mức SNR thấp. Toàn bộ quy trình nghiên cứu, xây dựng mã nguồn và kiểm chứng thực nghiệm được thực hiện xuyên suốt trong 9 tháng (từ tháng 01/2013 đến tháng 09/2013).

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à khảo sát các kịch bản mục tiêu đặt gần nhau đem lại các phát hiện khoa học quan trọng:

Thứ nhất, phương pháp Periodogram kết hợp thuật toán FFT truyền thống hoàn toàn mất khả năng phân tách khi khoảng cách giữa 2 mục tiêu có độ lệch tần số dưới 100 Hz (tương ứng khoảng cách cự ly dưới 9 cm trong môi trường có hằng số điện môi bằng 9). Khi tăng số điểm lấy mẫu $N_{FFT}$ từ 32 lên 1024 điểm, phổ FFT chỉ làm mượt đường cong hiển thị nhưng hai đỉnh phổ vẫn hoàn toàn dính liền thành một búp sóng duy nhất, dẫn đến sai số định vị mục tiêu lên tới 100%.

Thứ hai, việc ứng dụng các hàm cửa sổ trong phương pháp Periodogram cho thấy sự đánh đổi rõ rệt giữa độ rộng búp chính và biên độ búp phụ. Cửa sổ Blackman triệt tiêu biên độ búp sóng phụ xuống mức 0.0012 (giảm hơn 98.5% so với mức 0.08 của cửa sổ chữ nhật), giúp loại bỏ hoàn toàn hiện tượng rò rỉ phổ. Tuy nhiên, độ rộng búp chính của cửa sổ Blackman lại giãn rộng gấp 3 lần (đạt mức $3/N$ so với $1/N$ của cửa sổ chữ nhật), làm suy giảm nghiêm trọng độ phân giải cự ly khi các mục tiêu nằm sát nhau.

Thứ ba, các phương pháp không gian con như MUSIC và Superimposed MUSIC với bậc ma trận $p = 10$ và $p = 12$ chứng minh khả năng siêu phân giải vượt bậc ở điều kiện SNR bằng 30 dB. Thuật toán phân tách rõ nét hai đỉnh phổ có khoảng cách tần số cực hẹp với độ nhọn đỉnh phổ vượt trội. Tuy nhiên, khi tỷ số SNR giảm xuống mức 10 dB, thuật toán MUSIC nguyên bản bắt đầu xuất hiện các đỉnh phổ giả có biên độ tiệm cận đỉnh thực, gây nhiễu loạn kết quả nhận dạng.

Thứ tư, thuật toán cải tiến W-MUSIC thể hiện ưu thế vượt trội tuyệt đối trong kịch bản ba mục tiêu (S2), trong đó mục tiêu 1 có biên độ 0 dB (tương ứng giá trị 1.0), mục tiêu 2 suy hao -20 dB (biên độ 0.1) đặt gần mục tiêu 1, và mục tiêu 3 ở xa bị suy hao tới -40 dB (biên độ chỉ còn 0.01). W-MUSIC đã phân tách chính xác 100% cả ba vị trí mục tiêu, đồng thời triệt tiêu hoàn toàn các đỉnh phổ giả nhờ cơ chế nhân chập phổ FFT chuẩn hóa và phổ MUSIC cải tiến.

Thảo luận kết quả

Nguyên nhân tạo nên sự vượt trội của thuật toán W-MUSIC nằm ở cơ chế kết hợp thông minh giữa miền tần số Fourier và không gian con trực giao. Phổ chuẩn hóa của biến đổi Fourier $P_n{F}$ đóng vai trò như một bộ lọc không gian dải thông, giúp giữ nguyên năng lượng thực của tín hiệu phản xạ và ức chế các thành phần ngoài dải. Đồng thời, hàm phổ giả MUSIC cải tiến $P_n{M}$ cung cấp độ dốc cực lớn tại các tần số thực của mục tiêu. Khi tích hợp hai hàm này qua phép lấy căn bậc $K$ của tích các phân đoạn, các đỉnh phổ giả vốn không đồng nhất về pha giữa các phân đoạn sẽ bị triệt tiêu triệt để, trong khi các đỉnh mục tiêu thực được cộng hưởng và nâng cao độ sắc nét.

So sánh với các nghiên cứu sử dụng mô hình tự hồi quy Least-Square, phương pháp Least-Square đòi hỏi tỷ số SNR phải luôn duy trì trên 25 dB mới không bị phân rã ma trận sai số, trong khi W-MUSIC vẫn duy trì độ ổn định tuyệt vời ngay cả khi SNR dao động quanh mức 10 dB.

Các tập dữ liệu và kết quả ước lượng phổ trong nghiên cứu được trực quan hóa toàn diện thông qua các biểu đồ phổ 2D thể hiện biên độ theo tần số chuẩn hóa và hệ thống đồ thị không gian 3D quét qua dải SNR từ -20 dB đến +20 dB. Kết quả mô phỏng dạng bảng cũng đã lượng hóa chi tiết sự biến thiên của sai số định vị theo từng bậc ma trận tương quan ($p = 6, 8, 10, 12$), khẳng định bậc $p = 10$ là điểm tối ưu cân bằng giữa độ chính xác và chi phí tính toán vi xử lý.

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

Nhằm chuyển hóa các kết quả nghiên cứu lý thuyết và mô phỏng thành các giải pháp ứng dụng thực tế trong lĩnh vực đo đạc công trình ngầm, luận văn đưa ra bốn khuyến nghị cụ thể:

Thứ nhất, triển khai tích hợp thuật toán W-MUSIC lên các dòng chip xử lý tín hiệu số chuyên dụng DSP (như dòng TMS320C66x) hoặc FPGA nhằm thương mại hóa thiết bị GPR thời gian thực. Mục tiêu hành động là tối ưu hóa mã nguồn thuật toán để đạt tốc độ xử lý dưới 50 ms cho mỗi chu kỳ quét A-scan, do các nhóm nghiên cứu kỹ thuật điện tử và doanh nghiệp sản xuất thiết bị đo lường thực hiện trong giai đoạn 2024-2025.

Thứ hai, thiết lập quy trình chuẩn hóa tiền xử lý tín hiệu phản xạ tự động bù trừ suy hao môi trường. Áp dụng các bộ lọc số thích nghi để bù đắp mức suy hao 5 dB/m của đất ẩm và cân bằng biên độ giữa các mục tiêu nông và sâu trước khi đưa vào khối phân tích W-MUSIC, đặt mục tiêu mở rộng chiều sâu khảo sát chính xác lên trên 5 m, do các viện nghiên cứu địa kỹ thuật chủ trì hoàn thiện trong vòng 6 tháng.

Thứ ba, xây dựng cơ chế tự động nhận dạng bậc không gian con $p$ dựa trên tiêu chuẩn độ dài mô tả tối thiểu (MDL) hoặc tiêu chuẩn thông tin Akaike (AIC). Việc thay thế phương pháp chọn bậc cố định bằng thuật toán ước lượng bậc thích nghi sẽ giúp giảm 30% khối lượng tính toán ma trận trị riêng, hạn chế tối đa nguy cơ chọn thừa bậc gây sinh đỉnh giả khi SNR biến động dưới 10 dB.

Thứ tư, mở rộng thử nghiệm thực địa trên các mô hình địa tầng đa lớp phi đồng nhất và phát triển thuật toán từ định dạng A-scan lên xử lý ảnh quét B-scan và C-scan 3D. Các cơ quan quản lý hạ tầng giao thông và công trình ngầm đô thị cần phối hợp với các trường đại học thực hiện các dự án thử nghiệm trong thời gian 12 tháng nhằm lập bản đồ số hóa hệ thống cáp và cống ngầm đô thị với độ phân giải dưới 5 cm.

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

Luận văn là tài liệu học thuật và kỹ thuật chuyên sâu mang lại giá trị thiết thực cho bốn nhóm đối tượng:

• Kỹ sư nghiên cứu và phát triển phần cứng Radar: Tiếp cận chi tiết mô hình điều chế sóng liên tục FMCW, cấu trúc mạch trộn tần số và phương pháp trích xuất tín hiệu phách, giúp rút ngắn thời gian thiết kế hệ thống thu phát GPR băng thông rộng dưới 2 GHz.
• Chuyên gia địa kỹ thuật, xây dựng và khảo cổ học: Nắm bắt phương pháp luận phân tích dữ liệu radar không phá hủy, hiểu rõ các yếu tố ảnh hưởng như hằng số điện môi và hệ số suy hao để giải đoán chính xác các dị vật ngầm, đường ống kim loại và vết nứt kết cấu bê tông đặt cách nhau dưới 10 cm.
• Học viên cao học và nghiên cứu sinh ngành Xử lý tín hiệu số: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về toán học ứng dụng trong phân tích không gian con, từ thuật toán Pisarenko, Eigenvector, MUSIC cho đến phương pháp xây dựng thuật toán lai W-MUSIC trên Matlab.
• Các đơn vị quản lý hạ tầng đô thị và doanh nghiệp kiểm định: Cung cấp cơ sở khoa học để xây dựng phương án tự chủ công nghệ sản xuất thiết bị dò tìm công trình ngầm, giúp tiết kiệm từ 40% đến 60% chi phí đầu tư trang thiết bị so với việc nhập khẩu máy móc nguyên chiếc từ nước ngoài.

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

Phương pháp Radar xuyên đất GPR có ưu điểm gì vượt trội so với các phương pháp địa vật lý khác?

GPR hoạt động trên dải tần từ 1 MHz đến 1000 MHz, cung cấp khả năng khảo sát cấu trúc ngầm ở độ sâu tới 40 m đến 100 m mà hoàn toàn không xâm lấn hay phá hủy bề mặt. So với phương pháp địa chấn hay đo điện trở suất, GPR có tốc độ thu thập dữ liệu nhanh vượt trội, khả năng cơ động cao và cung cấp độ phân giải hình ảnh không gian chi tiết hơn nhiều lần.

Tại sao thuật toán FFT thông thường không thể phân biệt được các vật thể đặt gần nhau?

Thuật toán FFT bị giới hạn bởi tiêu chuẩn phân giải Rayleigh, phụ thuộc trực tiếp vào độ rộng băng thông phần cứng. Khi khoảng cách giữa 2 dị vật nhỏ hơn 9 cm, tín hiệu phản xạ bị chồng lấn khiến các búp sóng hòa làm một. Dù tăng số điểm lấy mẫu lên 1024 điểm, FFT chỉ làm mịn đường cong biểu diễn mà không thể tách được hai đỉnh tần số riêng biệt.

Thuật toán W-MUSIC khắc phục nhược điểm của thuật toán MUSIC truyền thống bằng cách nào?

Thuật toán MUSIC thông thường dễ tạo ra các đỉnh phổ giả khi tỷ số SNR giảm xuống mức 10 dB. W-MUSIC khắc phục triệt để bằng cách chia tín hiệu thành $K$ phân đoạn có cửa sổ trượt trùng lặp 50%, sau đó nhân phổ Fourier chuẩn hóa với phổ MUSIC cải tiến. Sự kết hợp này triệt tiêu hoàn toàn các đỉnh giả ngẫu nhiên và bảo toàn năng lượng thực của mục tiêu.

Hằng số điện môi của môi trường ảnh hưởng như thế nào đến kết quả ước lượng vị trí vật thể?

Hằng số điện môi tương đối quyết định trực tiếp vận tốc truyền sóng điện từ theo công thức $v = c/\sqrt{\varepsilon_r}$. Với môi trường đất ẩm có hằng số điện môi bằng 9, vận tốc sóng giảm 3 lần so với trong chân không. Nếu xác định sai hằng số điện môi, việc quy đổi từ độ dịch tần số của tín hiệu phách sang khoảng cách cự ly sẽ dẫn đến sai số vị trí nghiêm trọng.

Luận văn sử dụng công cụ và kịch bản nào để kiểm chứng độ chính xác của các thuật toán?

Tác giả sử dụng môi trường mô phỏng Matlab với chuỗi tín hiệu A-scan gồm 150 mẫu, tần số lấy mẫu 7.5 kHz. Hai kịch bản chính được khảo sát: kịch bản S1 gồm 2 mục tiêu có biên độ chênh lệch 20 dB, và kịch bản S2 gồm 3 mục tiêu với mục tiêu xa nhất suy hao tới -40 dB, được thử nghiệm dưới các mức nhiễu Gauss có SNR từ -20 dB đến +30 dB.

Kết luận

• Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết trường điện từ Maxwell và các kỹ thuật ước lượng phổ tín hiệu trong hệ thống Radar xuyên đất FMCW.
• Đã chỉ rõ giới hạn vật lý của phương pháp biến đổi Fourier và các hàm cửa sổ truyền thống khi xử lý các mục tiêu ngầm có khoảng cách tần số hẹp dưới 100 Hz.
• Đánh giá định lượng hiệu năng của các phương pháp không gian con (PHD, Eigenvector, MUSIC, Minimum-Norm) dưới tác động của nhiễu trắng Gauss ở các mức SNR từ 10 dB đến 30 dB.
• Phát triển và kiểm chứng thành công thuật toán lai W-MUSIC với khả năng siêu phân giải, tách bạch chính xác các dị vật đặt cách nhau dưới 9 cm và triệt tiêu 100% đỉnh phổ giả.
• Xác lập lộ trình chuyển giao công nghệ ứng dụng vi xử lý số DSP/FPGA trong vòng 1-2 năm tới nhằm chế tạo thành công thiết bị GPR nội địa phục vụ công tác xây dựng và quản lý đô thị.

Nghiên cứu mở ra hướng đi đầy triển vọng trong việc làm chủ công nghệ xử lý tín hiệu radar tiên tiến tại Việt Nam. Các kỹ sư, nhà khoa học và đơn vị chuyên môn được khuyến khích tiếp tục khai thác các thuật toán đề xuất để ứng dụng vào thực tế kiểm định công trình và số hóa không gian ngầm.

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

Chương 1. TỔNG QUAN VỀ RADAR XUYÊN ĐẤT Đây là chương giới thiệu tổng quan về Radar xuyên đất, về lịch sử phát triển, cho người đọc có cái nhìn tổng quan về Radar xuyên đất. Ngoài ra, chương này cũng đưa ra sơ đồ tổng quát của một hệ thống Radar, các bộ phận cấu thành hệ thống. Sau đó, tác giả sẽ trình bày sơ lược về kĩ thuật Radar xuyên đất, nguyên tắc hoạt động, ứng dụng, tầm quan trọng, lợi ích của nó trong việc dò tìm các vật chất nằm dưới bề mặt trái đất.1 Lịch sử phát triển Việc sử dụng các tín hiệu điện từ đầu tiên để xác định sự hiện diện của vật bằng kim loại từ xa trên đất liền được Hulsmeyer đưa ra vào năm 1904, nhưng sự mô tả đầu tiên của nó sử dụng cho xác định vị trí của các đồ vật được chôn xuất hiện sáu năm sau đó bởi một sáng chế của Leimbach và Lowy người Đức.

Kỹ thuật của họ bao gồm chôn một anten lưỡng cực trong một dãy các lỗ khoan thẳng đứng và so sánh các biên độ của tín hiệu nhận được khi các cặp antenna liên tiếp được sử dụng để truyền và nhận. Bằng cách này, một hình ảnh thô có thể được hình thành trong khu vực bất kỳ trong mảng đó. Các tác giả này mô tả một kỹ thuật thay thế, mà được sử dụng riêng biệt, bề mặt được gắn ăng ten để phát hiện các phản xạ từ một giao diện dưới bề mặt do nước ngầm hoặc quặng của một quặng sắt. Một phần mở rộng của kỹ thuật này đã dẫn đến một biểu thị độ sâu của một giao diện vật được chôn, thông qua kiểm tra giao thoa giữa các sóng phản xạ và rò rỉ trực tiếp giữa các ăng-ten trên mặt đất.

Tính năng chính của việc này, cụ thể là hoạt động sóng liên tục CW (Continuous Wave), sử dụng che chắn hoặc các hiệu ứng nhiễu xạ do các tính năng dưới đất, và sự phụ thuộc vào các biến thể dẫn xuất tạo ra tán xạ, đã được đưa ra trong một số bằng sáng chế khác, bao gồm một số dành cho các ứng dụng ngập nước hoàn toàn trong hầm mỏ. Nghiên cứu của Hiilsenbeck năm 1926 đưa ra đầu tiên sử dụng kỹ thuật xung để xác định cấu trúc các vật được chôn. Ông lưu ý rằng bất kỳ sự thay đổi chất điện môi, HVTH: Nguyễn Minh Linh 5 ƯỚC LƯỢNG PHỔ TÍN HIỆU GPR GVHD: TS. ĐỖ HỒNG TUẤN không nhất thiết liên quan đến dẫn xuất, cũng sẽ tạo ra sự phản xạ dễ dàng hơn trên các nguồn hướng, có lợi thế hơn các phương pháp địa chấn.

Kỹ thuật xung được phát triển từ những năm 1930 trở đi như một phương tiện thăm dò độ sâu đáng kể trong băng, nước ngọt, quặng muối, cát sa mạc và các vật có hình dạng đá. Thăm dò than đá cũng đã được nghiên cứu bởi Cook, Roe và Ellerbruch, mặc dù độ suy giảm cao hơn trong các tài liệu sau này, có nghĩa là có độ sâu lớn hơn vài mét đã không thực tế. Một phần mở rộng hơn lịch sử của Radar xuyên đất GPR (Ground Penetrating Radar) và sự phát triển của nó cho đến giữa thập niên 1970 được đưa ra bởi Nilsson. Sự quan tâm đổi mới trong các chủ đề được đưa ra trong đầu những năm 1970 khi nghiên cứu mặt trăng và đổ bộ lên mặt trăng đã được tiến hành.

Đối với các ứng dụng này, một trong những lợi thế của GPR hơn kỹ thuật địa chấn đã được khai thác, cụ thể là khả năng sử dụng từ xa, bộ chuyển đổi không tiếp xúc của năng lượng bức xạ, thay vì mặt đất liên hệ với các loại cần thiết cho điều tra địa chấn. Từ những năm 1970 cho đến ngày nay, một loạt các ứng dụng đã được mở rộng đều đặn, bao gồm xây dựng và thử nghiệm cấu trúc không phá hủy, khảo cổ học, đường xá và đánh giá chất lượng đường hầm, vị trí của khoảng trống và các vật chứa, đường ống và phát hiện dây cáp, cũng như viễn thám bằng vệ tinh… Mục đích xây dựng thiết bị cho mỗi các ứng dụng này đang được phát triển và người sử dụng hiện nay có một lựa chọn tốt hơn các thiết bị và kỹ thuật.2 Giới thiệu chung về radar xuyên đất Khả năng phát hiện ra các đối tượng được chôn dưới đất từ lâu đã được quan tâm qua nhiều thế kỷ. Một kỹ thuật mà có thể làm sáng tỏ vật dưới bề mặt đất là một thử thách đáng kể hấp dẫn các nhà khoa học nỗ lực nghiên cứu các kỹ thuật để đề ra phương pháp thăm dò phù hợp. Có nhiều phương pháp khác nhau dùng để quan sát các vật dưới lớp bề mặt đất như phương pháp địa chấn, điện trở suất, độ phân cực, khảo sát trọng lực, khảo sát từ, phương pháp nucleonic, bức xạ, thermographic và phương pháp điện từ…, mỗi phương pháp có ưu điểm riêng.

HVTH: Nguyễn Minh Linh 6 ƯỚC LƯỢNG PHỔ TÍN HIỆU GPR GVHD: TS. ĐỖ HỒNG TUẤN Phương pháp Radar xuyên đất GPR (Ground Penetrating Radar) dùng để thăm dò lớp dưới bề mặt trái đất là một lựa chọn đặc biệt thu hút đối với các kỹ sư và các nhà khoa học ứng dụng ở điểm nó bao trùm một loạt các chuyên ngành như truyền sóng điện từ, điện suy hao trong truyền thông, công nghệ anten ultra wideband và thiết kế các hệ thống radar, xử lý tín hiệu dạng sóng và xử lý hình ảnh biệt thức. Hầu hết các GPR là một ứng dụng cụ thể của công nghệ radar ultra wideband impulse. Thuật ngữ “Ground Penetrating Radar (GPR)”, hay “Groud Probing Radar”, hay “sub - surface radar” hoặc “surface-penetrating radar (SPR)” đều đề cập cho một loạt các kỹ thuật điện từ trường được thiết kế chủ yếu cho xác định vị trí của đối tượng hoặc các vật được chôn bên dưới bề mặt của trái đất hoặc nằm trong một cấu trúc nhìn không rõ.

Các công nghệ của GPR phần lớn được ứng dụng theo định hướng và thiết kế triết lý tổng thể, cũng như các phần cứng, thường phụ thuộc vào loại mục tiêu và nguyên liệu của các mục tiêu và môi trường xung quanh nó. Phạm vi của các ứng dụng cho phương pháp GPR rộng và sự tinh tế của các kỹ thuật phục hồi tín hiệu, thiết kế phần cứng và hoạt động thực tiễn gày càng nâng cao khi công nghệ ngày càng hoàn thiện.3 Nguyên lý hoạt động Radar xuyên đất - GPR là một phương pháp địa lý ứng dụng các nguyên lý của sóng điện từ ở dải tần số rất cao (1-1000MHz) để nghiên cứu cấu trúc và các đặc tính của vật chất bên dưới lòng đất mà không phải đào bới, phá hủy cấu trúc của nó. Thiết bị Radar xuyên đất sử dụng các sóng vô tuyến tần số cao để thu thập thông tin phản hồi từ dưới lòng đất. Năng lượng phát ra từ ăngten phát, lan truyền vào trong lòng đất với vận tốc phụ thuộc vào đặc tính điện môi của môi trường, khi gặp vật dị thường sẽ tạo ra các sóng phản xạ và được ăngten thu ghi lại các tín hiệu phản xạ này một cách liên tục, xử lý và tái tạo thành một hình ảnh.

HVTH: Nguyễn Minh Linh 7 ƯỚC LƯỢNG PHỔ TÍN HIỆU GPR GVHD: TS. ĐỖ HỒNG TUẤN Hình 1-1: Sơ đồ hoạt động tổng quát Radar xuyên đất Do các sóng phản xạ này được tạo ra từ những mặt ranh giới trung gian môi trường địa chất, nên các sóng phản xạ thường liên quan đến những điều kiện tạo thành tự nhiên trong cấu trúc địa chất như: ranh giới đá móng, các lớp vật liệu trầm tích có tính vật lý khác nhau, nồng độ sét, những khuyết tật, các khe nứt nẻ, các khối xâm thực cũng như các vật liệu bị chôn vùi do nhân tạo hoặc các khối bê tông, các vật thể không đồng nhất liên quan tới vị trí hang hốc, hàm ếch, tổ mối,… Một trong những vấn đề lớn nhất của Radar xuyên đất là phản xạ từ mặt đất quá lớn, với hằng số điện môi cao giữa mặt đất và không khí chỉ cho phép một lượng nhỏ năng lượng truyền qua mặt phân cách, phản xạ vào mục tiêu và đi qua mặt phân cách đến Ăngten nhận. Độ sâu thẩm thấu của các phương pháp phụ thuộc vào tần số của Ăngten phát - thu tín hiệu và phụ thuộc vào tính chất của đất đá trong môi trường địa chất. Các loại Ăngten thông thường được sử dụng để khảo sát cấu trúc địa chất có tần số là 12,5: 25; 100; 200;400 MHz và độ sâu khảo sát 40 - 100 m.

HVTH: Nguyễn Minh Linh 8 ƯỚC LƯỢNG PHỔ TÍN HIỆU GPR GVHD: TS. ĐỖ HỒNG TUẤN Sau đây là hình ảnh minh họa cho quá trình truyền, bị phản xạ, khúc xạ của sóng điện từ khi đi vào trong lòng đất: Hình 1- 2 Quá trình truyền và nhận của sóng GPR khi đi vào lòng đất Hình dưới đây là đường cong chỉ thời gian truyền từ bộ phát đến bộ nhận của mỗi dạng sóng ở trên. Hình 1- 3 Đường cong thời gian truyền của mỗi dạng sóng Khi sóng được phát từ nguồn phát đi đến các mục tiêu trong lòng đất thì bộ thu sẽ nhận các sóng đến gồm các loại khác nhau như sau: sóng đến trực tiếp từ trong không khí và trong lòng đất, sóng phản xạ và sóng khúc xạ. Đường đi của sóng khúc xạ là xa nhất, và gần nhất là sóng đến trực tiếp.

HVTH: Nguyễn Minh Linh 9 ƯỚC LƯỢNG PHỔ TÍN HIỆU GPR GVHD: TS.4 Các định nghĩa và các đặc điểm của tín hiệu radar xuyên đất 1.1 Dynamic range Bộ nhận radar có khả năng xử lý tín hiệu lớn từ bề mặt phản xạ và các mục tiêu có khoảng cách ngắn, nó cũng có thể phân tích các tín hiệu nhỏ gần tầng nhiễu. Tỉ số giữa tín hiệu có thể nhận được lớn nhất trên tín hiệu nhỏ nhất có thể phân tích được gọi là khoảng động (dynamic range). Dynamic range được xác định bởi biểu thức sau: 𝑉 𝐷𝑦𝑛𝑎𝑚𝑖𝑐𝑅𝑎𝑛𝑔𝑒 = 20𝑙𝑜𝑔 ( 𝑚𝑎𝑥 ) (1.1) 𝑉𝑚𝑖𝑛 Dynamic Range có đơn vị là decibel (dB). Vmax (V) là tín hiệu lớn nhất có thể nhận được.

Vmin (V) là tín hiệu nhỏ nhất có thể phát hiện ra. Dynamic Range của hệ thống sẽ ảnh hưởng tới khoảng cách lớn nhất (maximum range) có thể dò tìm được. Thông thường radar sẽ có dynamic range của hệ thống lớn hơn dynamic range lấy mẫu. Dynamic range của một bộ ADC bằng20𝑙𝑜𝑔(2𝑁), với N là số bit, hay xấp xỉ N lần 6dB.

Do vậy ADC 16 bit sẽ có dynamic range lý thuyết là 96 dB.2 Băng thông Băng thông hệ thống ký hiệu là B được định nghĩa là nghịch đảo của độ rộng xung τp cho radar xung. Băng thông thường được lấy trung tâm và trong thực tế là xấp xỉ bằng với tần số trung tâm GPR xung fc.