Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu về cảm biến thụ động không dây dạng sóng âm bề mặt" của nghiên cứu sinh Nguyễn Thu Hà, thuộc ngành Kỹ thuật điều khiển và tự động hóa (mã số 9520216), được thực hiện tại Đại học Bách khoa Hà Nội dưới sự hướng dẫn khoa học của PGS.TS Hoàng Sĩ Hồng và TS. Cung Thành Long (năm 2023).

Tính cấp thiết và khoảng trống nghiên cứu

Trong xu thế phát triển của mạng cảm biến không dây và Internet kết nối vạn vật (IoT), nhu cầu giám sát thông số vận hành tại các vị trí khắc nghiệt (nhiệt độ cao, áp suất lớn, các bộ phận quay hoặc trong môi trường bức xạ) đặt ra yêu cầu về các loại cảm biến không cần sử dụng pin cấp nguồn và không cần dây dẫn. Cảm biến sóng âm bề mặt (Surface Acoustic Wave - SAW) đáp ứng được các tiêu chuẩn này nhờ khả năng tự thu nhận năng lượng từ sóng vô tuyến thẩm vấn (RF) phát ra từ bộ đọc, kích thước nhỏ gọn do vận tốc sóng âm bề mặt nhỏ hơn sóng điện từ xấp xỉ $10^5$ lần.

Tuy nhiên, việc ứng dụng cảm biến SAW loại đường trễ phản xạ (Reflective Delay Line - RDL) trong các hệ thống đo lường thực tế gặp phải ba khoảng trống kỹ thuật chính:

  1. Mối quan hệ định lượng giữa các thông số hình học của cấu trúc cảm biến (số cặp điện cực IDT, số thanh phản xạ, tỉ số hóa kim, độ dày màng kim loại, vật liệu) đến biên độ sóng phản xạ hồi đáp chưa được khảo sát toàn diện.
  2. Hiện tượng nhập nhằng pha xuất hiện khi đo lường dựa trên độ trễ pha nếu góc trễ pha vượt quá $2\pi$, dẫn đến một giá trị lượng giác tương ứng với nhiều góc pha khác nhau, gây sai lệch thông tin đo.
  3. Khi triển khai hệ đa cảm biến trong cùng một trường thẩm vấn của bộ đọc, hiện tượng chồng lấn tín hiệu đồng thời phản hồi từ các cảm biến làm phát sinh xung đột, trong khi các giải pháp hiện có chỉ giải quyết được ở dải tần số thấp (như 12,5 MHz), chưa đáp ứng được dải tần số cao hàng trăm MHz trong các mạng cảm biến công nghiệp.

Mục tiêu nghiên cứu

Luận án xác định 3 mục tiêu cụ thể:

  1. Nghiên cứu sự ảnh hưởng của các thông số cấu trúc cảm biến SAW đến biên độ sóng phản xạ.
  2. Xây dựng phương pháp tính toán góc pha của tín hiệu phản xạ.
  3. Xây dựng thuật toán đọc thông tin phản hồi đồng thời từ các cảm biến thành phần trong hệ đa cảm biến SAW RDL được mã hóa theo phương pháp tần số trực giao (Orthogonal Frequency Coding - OFC).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Cảm biến thụ động không dây dạng sóng âm bề mặt cấu trúc đường trễ phản xạ (SAW RDL), thuật toán đo pha tín hiệu phản xạ và thuật toán tách tín hiệu chồng lấn trong hệ đa cảm biến mã hóa tần số trực giao (OFC).
  • Phạm vi nghiên cứu:
    • Nghiên cứu tập trung vào cấu trúc cảm biến SAW RDL một cổng gồm một bộ chuyển đổi số (IDT) và các bộ phản xạ đặt cùng một phía trên đế áp điện $YZ\text{-LiNbO}_3$.
    • Khảo sát đáp ứng mô phỏng tại các tần số 125 MHz, 870 MHz và thuật toán giải mã tín hiệu đa cảm biến mở rộng lên tần số 500 MHz.
    • Luận án giới hạn phạm vi mô hình hóa ở cấp độ linh kiện và thuật toán xử lý tín hiệu tại bộ đọc; giả định tín hiệu phản hồi từ IDT được truyền trọn vẹn đến bộ đọc, không nghiên cứu sâu quá trình bức xạ của anten và sự suy hao trên kênh truyền không dây trong môi trường thực tế.

Tổng quan tài liệu và vị trí của luận án

Các công trình nền tảng về sóng âm bề mặt khởi đầu từ phát hiện của Lord Rayleigh (1885) về sự lan truyền sóng đàn hồi trên bề mặt vật rắn, tiếp nối bởi phát hiện về hiệu ứng áp điện thuận và ngược của Pierre và Paul-Jacques Curie (1880–1881). Năm 1965, White và Voltmer chế tạo bộ chuyển đổi số (IDT) trên đế áp điện, đặt nền móng cho việc chuyển đổi tín hiệu điện sang sóng cơ học bề mặt. Đến giai đoạn 1969–1970, Tancrell và cộng sự phát triển các bộ lọc SAW thương mại. Cấu trúc SAW RDL cho các ứng dụng đo lường không dây bắt đầu được nghiên cứu từ năm 1995.

Luận án đã điểm lại các công trình liên quan trực tiếp đến ba hướng nghiên cứu chính:

Hướng nghiên cứu Tác giả, công trình Kết quả đạt được Giới hạn tồn tại mà luận án giải quyết
Ảnh hưởng của thông số cấu trúc đến sóng phản xạ Tác giả tại tài liệu [20] Khảo sát hệ số phản xạ trên đế $\text{LiNbO}_3$ $YX\text{-}128^\circ$, $\text{LiNbO}_3$ $YX\text{-}64^\circ$ và $\text{LiTaO}_3$; tỉ số $h/\lambda$ từ 2% đến 8%, tỉ số hóa kim $a/p$ từ 0,4 đến 0,6 bằng mô hình ghép cặp chế độ (COM). Chưa kiểm soát chính xác độ rộng điện cực IDT; chưa công bố chi tiết kích thước mô hình; chưa đánh giá đầy đủ số lượng điện cực và thời gian kích thích.
Tác giả tại tài liệu [21] Sử dụng IDT một hướng SPUDT để giảm tổn hao, đế $\text{LiNbO}_3$, điện cực nhôm; biên độ phản xạ đầu tiên giảm 50% so với tín hiệu vào. Chỉ khảo sát trên 1 cặp điện cực IDT và 2 bộ phản xạ; phạm vi khảo sát khoảng cách hẹp, chưa bao quát các tham số hình học khác.
Xác định pha và xử lý nhập nhằng pha S. Schuster và cs [23] Ước tính thô nhiệt độ dựa trên độ trễ giữa các cặp phản xạ để xác định góc pha tuyệt đối; dùng hiệu pha $\Delta\varphi_{21}$ và $\Delta\varphi_{31}$. Làm xuất hiện hệ số không đảm bảo, làm giảm tỷ số tín hiệu trên nhiễu (SNR); yêu cầu vị trí bộ phản xạ khắt khe nhưng chưa làm rõ vị trí.
L. Reindl [24] Dùng cận dưới Cramer-Rao (CRLB) và phép biến đổi Fourier rời rạc để ước tính độ trễ từ pha; thêm độ lệch pha tham chiếu. Phải ước tính hệ số trễ nhiệt độ (TCD); giả thiết SNR cố định; bỏ qua thành phần sóng bậc cao, ảnh hưởng đến độ chính xác.
C. Gruber và cs [26] Mô hình hóa quá trình đọc tín hiệu phản hồi từ 3 bộ phản xạ trên phần mềm Matlab, tính được $\tan(\varphi)$. Dữ liệu đầu vào lấy ngẫu nhiên; chưa giải quyết hiện tượng nhập nhằng pha khi góc trễ pha lớn hơn $2\pi$.
Tài liệu [32] Đặt thêm một bộ phản xạ gần IDT nhất có độ trễ pha $< 2\pi$ để xác định số chu kỳ quay pha $n$ ở tần số 2,4 GHz, khoảng cách 1 m. Cấu trúc và vị trí cảm biến phải cố định trước; chỉ thử nghiệm ở cự ly ngắn; chưa xử lý được trường hợp tín hiệu chồng lấn liên tục.
Đọc đồng thời và giải quyết chồng lấn tín hiệu trong hệ đa cảm biến Sorokin và cs [30] Kết hợp OFC và mã hóa vị trí thời gian xung cho hệ 6 cảm biến SAW; xây dựng ma trận thời gian - tần số bằng các bit nhị phân "0" và "1". Tách riêng chức năng đo và nhận dạng làm tăng kích thước cảm biến; chỉ áp dụng cho 6 cảm biến cố định để chủ động tránh xung đột, chưa giải quyết được bài toán chồng lấn thực sự.
Tài liệu [31] Dùng thuật toán Neyman-Pearson và các bộ dò để nhận dạng thẻ SAW dựa trên tổ hợp bit trong khe thời gian. Chỉ ứng dụng cho mục đích nhận dạng thẻ, không trích xuất được thông tin đo lường chi tiết (pha, độ trễ liên tục).
P. Caldero và cs [33] Dùng mảng anten trên bộ đọc, tính hàm tương quan và mật độ xác suất không gian 3D để tách tín hiệu. Giả thiết độ lệch pha tương đối giữa các xung là như nhau; bỏ qua việc tính chính xác độ trễ pha, dễ mất tín hiệu do nhiễu.
Weifeng Liu và cs [34] Xây dựng ma trận nhận dạng, chứng minh định lý hạng ma trận và giải ma trận lượng giác mô tả quan hệ tín hiệu chồng lấn với số bộ phản xạ $N_c$ lên đến 10. Thuật toán chỉ cho kết quả chính xác ở tần số thấp (tần số trung tâm $12,5\text{ MHz}$); khi tần số tăng lên hàng trăm MHz, thuật toán mất tính chính xác.

Vị trí của luận án: Kế thừa khung lý thuyết giải ma trận của Weifeng Liu [34] và nguyên lý mã hóa tần số trực giao của D. Malocha [72, 76], luận án tập trung giải quyết triệt để hai bài toán còn bỏ ngỏ: tính toán chính xác góc trễ pha khi vượt ngưỡng $2\pi$ trên mô hình cấu trúc 3 bộ phản xạ, và xây dựng thuật toán chia khoảng kết hợp lọc nghiệm để giải hệ phương trình ma trận khi nâng tần số hoạt động của hệ đa cảm biến lên 500 MHz.

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

Cơ sở lý thuyết

  1. Hiệu ứng áp điện và lan truyền sóng âm bề mặt: Sự tương tác điện - cơ trong tinh thể áp điện được mô tả qua hệ phương trình cấu thành: $${T} = [c^E]{S} - [e]^t{E}$$ $${D} = [e]{S} + [\varepsilon^S]{E}$$ Trong đó: ${T}$ là vectơ ứng suất cơ, ${S}$ là vectơ biến dạng cơ, ${E}$ là vectơ cường độ điện trường, ${D}$ là vectơ dịch chuyển điện; $[c^E]$, $[\varepsilon^S]$, $[e]$ lần lượt là ma trận độ cứng, ma trận hằng số điện môi và ma trận áp điện. Sóng Rayleigh lan truyền tập trung năng lượng trong độ sâu tương đương một bước sóng bề mặt, với các hạt vật chất chuyển động theo quỹ đạo elip trong mặt phẳng thẳng đứng dọc theo hướng truyền sóng.
  2. Cấu trúc bộ chuyển đổi IDT và điều kiện phản xạ Bragg:
    • Bước sóng âm bề mặt: $\lambda = 4d$ (đối với cấu trúc điện cực đơn hai hướng), trong đó $d$ là độ rộng điện cực. Tần số trung tâm: $f_0 = v/\lambda$, với $v$ là vận tốc lan truyền sóng âm bề mặt.
    • Chiều dài IDT liên quan đến số cặp điện cực $N_p$: $l_{IDT} = N_p \cdot p$, ảnh hưởng trực tiếp đến băng thông của cảm biến.
    • Để đạt hệ số phản xạ cực đại từ các thanh phản xạ, điều kiện phản xạ Bragg được áp dụng: bước sóng cơ học $\lambda = 2p$, với $p$ là chu kỳ của lưới thanh phản xạ.
  3. Mã hóa tần số trực giao (OFC): Tín hiệu phản xạ của cảm biến SAW mã hóa OFC được định nghĩa bởi hàm truyền đạt lý tưởng trên miền thời gian: $$h_{OFC}(t) = \sum_{k=1}^{N_C} a_k \text{rect}\left(\frac{t - \tau_{pk}}{\tau_{chip0}}\right) \cos\left[2\pi f_{chipk}(t - \tau_{pk})\right]$$ Trong đó: $a_k$ là biên độ sóng phản xạ thứ $k$, $\tau_{pk}$ là độ trễ thời gian của sóng phản xạ thứ $k$, $N_C$ là số tần số trực giao thành phần (tương ứng số bộ phản xạ), $\tau_{chip0}$ là độ dài thời gian của một chip, và $f_{chipk}$ là tần số trung tâm của chip thứ $k$.

Lựa chọn vật liệu

  • Vật liệu đế áp điện: Luận án chọn tinh thể Lithium Niobate cắt hướng YZ ($YZ\text{-LiNbO}_3$). Vật liệu này sở hữu hệ số biến đổi điện cơ cao ($K^2 = 4,6%$), vận tốc truyền sóng âm bề mặt lý thuyết là $v = 3480\text{ m/s}$, và hệ số nhiệt độ trễ $\text{TCD} = 94\text{ ppm}/^\circ\text{C}$, phù hợp cho các cấu trúc cảm biến đòi hỏi biên độ tín hiệu phản hồi lớn.
  • Vật liệu chế tạo điện cực IDT và thanh phản xạ: Kim loại Nhôm (Al) được lựa chọn nhờ điện trở suất tương đối thấp ($2,65\ \mu\Omega\cdot\text{cm}$), độ bám dính bề mặt tốt, nhiệt độ sôi thấp ($2792\text{ K}$) giúp quy trình bốc hơi nhiệt và lắng đọng màng mỏng thuận lợi, đồng thời có khối lượng riêng nhỏ giúp giảm tải trọng cơ học lên bề mặt đế.

Phương pháp nghiên cứu và công cụ

Luận án kết hợp nghiên cứu lý thuyết giải tích và mô phỏng số thông qua hai công cụ phần mềm:

  • Phương pháp phần tử hữu hạn (FEM) trên phần mềm Ansys: Xây dựng mô hình 2D thực tế của cảm biến SAW RDL trên đế $YZ\text{-LiNbO}_3$. Phương pháp này giải hệ phương trình vi phân liên kết áp điện, mô tả quá trình lan truyền sóng cơ học từ IDT, phản xạ tại các thanh phản xạ và biến đổi ngược thành điện áp trên IDT.
  • Công cụ tính toán số Matlab/Simulink: Tiếp nhận dữ liệu đáp ứng điện áp xuất ra từ Ansys để xây dựng khối thuật toán đo pha, mô hình hóa bộ lọc giải điều chế đồng pha/vuông pha (I/Q), và lập trình giải hệ phương trình ma trận phi tuyến để tách tín hiệu chồng lấn trong miền tần số và thời gian.

Nội dung chính theo từng chương

Chương 1. Tổng quan về cảm biến thụ động không dây dạng sóng âm bề mặt

Chương 1 hệ thống hóa cơ sở lý thuyết về cảm biến không dây, phân biệt hai nhóm cảm biến tích cực (sử dụng pin) và thụ động (tự cấp nguồn qua sóng RF hoặc năng lượng môi trường). Tác giả so sánh hai cấu trúc cảm biến SAW phổ biến:

Đặc tính so sánh Cấu trúc SAW RDL Cấu trúc SAW cộng hưởng
Khả năng nhúng mã nhận dạng Có khả năng nhúng mã đa dạng Không có
Độ phân giải phép đo Thấp hơn Cao
Cơ chế xử lý tín hiệu hồi đáp Giải mã theo độ trễ thời gian hoặc độ lệch pha của xung phản xạ Giải mã thông qua độ dịch tần số cộng hưởng
Phạm vi ứng dụng Rộng (kết hợp đo lường và nhận dạng ID) Hẹp hơn

Nội dung chương phân tích chi tiết cấu tạo các loại IDT: hai hướng điện cực đơn ($\lambda/4$), hai hướng điện cực kép ($\lambda/8$), và một hướng SPUDT. Đồng thời, chương 1 trình bày ba phương pháp mã hóa trên thẻ SAW RDL: mã hóa vị trí thời gian (Time Position Coding), mã hóa pha (Phase Coding), và mã hóa tần số trực giao (OFC). Với phương pháp OFC, băng thông tín hiệu được trải rộng và các tần số chip thỏa mãn điều kiện trực giao toán học, cho phép nhiều cảm biến cùng hoạt động mà không gây can nhiễu lẫn nhau nếu được phân bổ mã phù hợp.

Chương 2. Xây dựng thuật toán đọc pha của tín hiệu phản hồi

Chương 2 tập trung khảo sát định lượng ảnh hưởng của các tham số hình học cảm biến SAW RDL đến biên độ sóng phản xạ bằng mô phỏng FEM 2D trên Ansys, sau đó xây dựng thuật toán xác định góc pha trên Matlab/Simulink.

Khảo sát ảnh hưởng của các tham số hình học (Mô phỏng FEM trên Ansys)

Mô hình cơ sở gồm 1 IDT và 3 bộ phản xạ hở mạch đặt trên đế $YZ\text{-LiNbO}_3$, được khảo sát chủ yếu tại tần số 870 MHz (bước sóng $\lambda = 4\ \mu\text{m}$, chu kỳ $p = 2\ \mu\text{m}$, độ rộng điện cực $d = 1\ \mu\text{m}$) và đối chứng tại tần số 125 MHz:

  1. Thời gian đặt điện áp kích thích tại IDT: Thời gian cấp xung ảnh hưởng trực tiếp đến số chu kỳ dao động được thiết lập trên IDT. Khảo sát chỉ ra rằng cần duy trì xung kích thích đủ số chu kỳ để dao động đạt trạng thái xác lập, tối ưu hóa năng lượng chuyển đổi sang sóng cơ.
  2. Số cặp điện cực IDT ($N_p$): Khi tăng số điện cực IDT, biên độ sóng phản xạ tăng lên do hiệu suất chuyển đổi điện - cơ tăng và băng thông thu hẹp; tuy nhiên, số lượng điện cực quá lớn sẽ gây hiện tượng tự phản xạ nội tại trong IDT, làm méo dạng sóng.
  3. Số thanh của bộ phản xạ: Khảo sát cho thấy biên độ sóng phản xạ tăng tỉ lệ thuận với số thanh phản xạ trên mỗi bộ cho đến khi đạt giới hạn bão hòa phản xạ Bragg.
  4. Tỉ số hóa kim ($a/p$): Tỉ số giữa khoảng cách điện cực $a$ và chu kỳ $p$ được khảo sát tại các giá trị từ 0,35 đến 0,65. Kết quả ghi nhận tại $a/p = 0,5$ (độ rộng thanh bằng khoảng cách giữa hai thanh, $d = a = \lambda/4$) cho biên độ sóng phản xạ lớn nhất.
  5. Vật liệu thanh phản xạ: So sánh giữa các vật liệu kim loại Nhôm (Al), Vàng (Au), Đồng (Cu), Titan (Ti) và Tungsten. Nhôm cho đáp ứng biên độ ổn định, ít gây suy hao tán xạ do khối lượng lên bề mặt tinh thể.
  6. Độ dày thanh phản xạ ($h_1/\lambda$): Khảo sát tỉ số độ dày màng kim loại trên bước sóng $h_1/\lambda$ từ 1% đến 6% tại 870 MHz. Độ dày tăng làm tăng hệ số phản xạ cơ học nhưng nếu quá dày sẽ gây tán xạ khối và suy giảm vận tốc pha.
Cấu hình cấu trúc tối ưu cho mô hình tính pha:
- Đế áp điện: YZ-LiNbO3 (v = 3480 m/s)
- Cấu trúc: 1 IDT hai hướng điện cực đơn + 3 bộ phản xạ hở mạch
- Vật liệu điện cực và thanh phản xạ: Nhôm (Al)
- Tỉ số hóa kim: a/p = 0,5
- Tần số trung tâm: 870 MHz (λ = 4 µm)

Xây dựng thuật toán đọc pha và khử nhập nhằng pha

Tín hiệu điện áp phản hồi từ 3 bộ phản xạ (ký hiệu 1, 2, 3) tại IDT thu được từ mô phỏng FEM được nhập vào mô hình xử lý tín hiệu trong Matlab/Simulink.

  • Thuật toán trích xuất các thành phần đồng pha $I(t)$ và vuông pha $Q(t)$ thông qua bộ trộn tần và bộ lọc tương tự, từ đó xác định góc pha tức thời: $\varphi_i = \arctan(Q_i / I_i)$.
  • Để giải quyết hiện tượng nhập nhằng pha khi góc trễ pha thực tế lớn hơn $2\pi$, thuật toán sử dụng mối tương quan giữa hai hiệu pha $\Delta\varphi_{21} = \varphi_2 - \varphi_1$ và $\Delta\varphi_{31} = \varphi_3 - \varphi_1$, kết hợp với tỷ số khoảng cách vật lý $L_2 - L_1$ và $L_3 - L_1$ giữa các bộ phản xạ trên đế áp điện.
  • Việc so sánh tỉ lệ chênh lệch thời gian trễ thô với bước nhảy pha $2\pi \cdot n$ cho phép thuật toán khôi phục chính xác số nguyên chu kỳ quay pha $n$, xác định giá trị pha tuyệt đối của tín hiệu phản xạ mà không bị giới hạn trong khoảng $[-\pi, \pi]$ hay $[0, 2\pi]$.

Chương 3. Xây dựng thuật toán đọc đồng thời các tín hiệu phản xạ trong hệ đa cảm biến SAW OFC

Chương 3 giải quyết bài toán cốt lõi khi nhiều cảm biến SAW RDL mã hóa OFC cùng gửi tín hiệu về một bộ đọc, gây ra hiện tượng chồng lấn tín hiệu trong cả miền thời gian và tần số.

Mô hình toán học của tín hiệu chồng lấn

Tín hiệu tổng thu được tại anten của bộ đọc là sự xếp chồng tuyến tính của tín hiệu phản xạ từ $N_t$ cảm biến thành phần trong trường thẩm vấn. Hệ phương trình ma trận mô tả mối quan hệ giữa tín hiệu chồng lấn tại bộ đọc và tín hiệu của từng cảm biến được biểu diễn dưới dạng: $$Y_{KNc} = B_{eNc} \times H_{Nc}$$ Trong đó:

  • $Y_{KNc}$: Ma trận tín hiệu chồng lấn đo được tại bộ đọc sau khi lấy mẫu tại các tần số trực giao.
  • $B_{eNc}$: Ma trận mã hóa biểu diễn sự có mặt của các chip tần số trực giao trên từng cảm biến.
  • $H_{Nc}$: Ma trận chứa ẩn số cần tìm, bao gồm biên độ ($a_{mk}$) và độ trễ pha/thời gian ($\tau_{mk}, \varphi_{mk}$) của từng sóng phản xạ từ cảm biến thứ $m$.

Thuật toán chia khoảng và lọc nghiệm

Khi tần số hoạt động của cảm biến nâng lên mức hàng trăm MHz (chẳng hạn 500 MHz), số nghiệm dao động giả trong phương trình lượng giác tăng mạnh, dẫn đến việc giải hệ phương trình phi tuyến trực tiếp bằng các phương pháp thông thường (như trong công bố của Weifeng Liu ở 12,5 MHz) bị phân kỳ hoặc cho nghiệm sai. Luận án đề xuất:

  1. Thuật toán chia khoảng: Dựa trên bán kính vùng hỏi của bộ đọc ($d_M = 6\text{ m}$) và thời gian trễ lan truyền cực đại $\tau_M$, miền nghiệm thời gian và tần số được chia nhỏ thành các khoảng con xác định.
  2. Lọc nghiệm: Áp dụng điều kiện biên vật lý về khoảng cách giữa các bộ phản xạ trên thẻ cảm biến để loại bỏ các nghiệm ngoại lai không thỏa mãn tính trực giao và tính liên tục thời gian.
  3. Giải tuần tự ma trận $Y_{Nc}$ cho từng khung thời gian tương ứng với từng chip để tách độc lập biên độ và pha của từng cảm biến.

Kết quả kiểm chứng mô phỏng trên Matlab

Thuật toán được kiểm chứng qua các kịch bản mô phỏng số với số lượng cảm biến tăng dần:

  • Hệ 1 cảm biến với 3 bộ phản xạ ($N_c = 3$): Xác nhận tính đúng đắn của thuật toán giải ma trận trong việc tách các đỉnh phản xạ riêng biệt từ tín hiệu tổng.
  • Hệ 3 cảm biến với 3 bộ phản xạ ($N_c = 3$): Phục hồi chính xác biên độ và thời gian trễ của từng cảm biến trong vùng hỏi bán kính 6 m mà không xuất hiện sai số chồng lấn.
  • Hệ 5 cảm biến: Kiểm chứng khả năng phân tách tín hiệu khi mức độ chồng lấn dạng sóng tăng cao.
  • Hệ 10 cảm biến với 4 bộ phản xạ ($N_c = 4$): Mô phỏng toàn diện với 10 thẻ cảm biến hoạt động đồng thời tại tần số trung tâm lên tới 500 MHz. Kết quả sau khi giải hệ phương trình ma trận kết hợp chia khoảng và lọc nghiệm đã nhận dạng chính xác sự có mặt của toàn bộ 10 cảm biến trong vùng hỏi của thiết bị đọc, đồng thời bóc tách thành công biên độ và góc pha của từng sóng phản xạ thành phần.

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học và kỹ thuật

  1. Phương pháp xác định góc trễ pha vượt ngưỡng $2\pi$: Luận án đề xuất phương pháp tính toán góc trễ pha giữa tín hiệu phản xạ và tín hiệu thẩm vấn trên cấu trúc cảm biến SAW RDL 3 bộ phản xạ. Phương pháp kết hợp phân tích tín hiệu trên mô hình Simulink và tỷ số độ trễ hình học, giải quyết triệt để bài toán nhập nhằng pha khi góc lệch pha lớn hơn $2\pi$, đảm bảo tính đơn trị và nâng cao độ chính xác của phép đo.
  2. Thuật toán đọc đồng thời hệ đa cảm biến SAW OFC ở tần số cao: Luận án nghiên cứu và đề xuất thuật toán giải hệ phương trình ma trận liên kết giữa tín hiệu chồng lấn và tín hiệu thành phần bằng phương pháp chia khoảng và lọc nghiệm. Thuật toán cho phép nhận dạng và tách đồng thời tín hiệu phản hồi từ 10 cảm biến SAW mã hóa OFC với số tần số thành phần $N_c = 4$ tại dải tần số hoạt động đạt 500 MHz (vượt qua giới hạn 12,5 MHz của các nghiên cứu tiền nhiệm).
  3. Bộ dữ liệu khảo sát tham số hình học toàn diện: Cung cấp các kết quả định lượng chi tiết về ảnh hưởng của thời gian kích thích IDT, số điện cực IDT, số thanh phản xạ, tỉ số hóa kim ($a/p = 0,5$) và tỉ số độ dày kim loại ($h_1/\lambda$) trên đế $YZ\text{-LiNbO}_3$, thiết lập cơ sở dữ liệu cho việc thiết kế cấu trúc phần cứng SAW RDL tối ưu biên độ sóng phản xạ.

Kiến nghị và giải pháp đề xuất

  • Ứng dụng cấu trúc SAW RDL tối ưu hóa thông số hình học để nâng cao tỷ số tín hiệu trên nhiễu (SNR) cho các hệ thống đo lường không dây trong công nghiệp.
  • Ứng dụng thuật toán ma trận chia khoảng và lọc nghiệm vào phần mềm nhúng của các bộ đọc vô tuyến, phục vụ giám sát nhiệt độ và ứng suất trực tuyến tại các vị trí khó tiếp cận như tiếp điểm tủ điện phân phối, trạm biến áp, và các chi tiết cơ khí chuyển động quay.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Các phân tích và đánh giá ảnh hưởng tham số cấu trúc chủ yếu thực hiện qua mô hình phần tử hữu hạn 2D trên phần mềm Ansys; chưa triển khai mô hình 3D hoàn chỉnh do giới hạn về năng lực tính toán phần cứng, dù mô hình 2D đã phản ánh được các đặc tính lan truyền sóng cơ bản.
  • Luận án giả định tín hiệu phản hồi từ IDT được chuyển giao nguyên vẹn về bộ đọc; chưa mô hình hóa chi tiết ảnh hưởng của anten thu/phát và các hiện tượng fading, suy hao, phản xạ đa đường của kênh truyền sóng vô tuyến trong môi trường thực tế.
  • Phần thực nghiệm chế tạo phần cứng mới dừng lại ở bước khảo sát chế tạo thử nghiệm trên mask, đo đạc thông số tán xạ $S_{11}$ trên máy phân tích mạng véc-tơ (VNA) và thiết kế mạch phối hợp trở kháng (trình bày tại phần Phụ lục), chưa tích hợp thành hệ thống đa cảm biến hoàn chỉnh ngoài hiện trường.

Hướng nghiên cứu tiếp theo

  • Mở rộng mô hình mô phỏng FEM sang không gian 3 chiều để khảo sát sâu hơn hiệu ứng nhiễu xạ bờ và sóng khối sinh ra bên dưới bề mặt đế.
  • Tích hợp mô hình kênh truyền sóng RF thực tế và đặc tính anten vào thuật toán xử lý tín hiệu tại bộ đọc.
  • Chế tạo hoàn thiện mảng cảm biến vật lý và mạch phần cứng bộ đọc hoàn chỉnh để kiểm chứng thực nghiệm thuật toán giải ma trận tách tín hiệu chồng lấn trong môi trường công nghiệp thực tế.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Chuyên ngành Kỹ thuật Điều khiển và Tự động hóa, Kỹ thuật Đo lường và Hệ thống cảm biến, Kỹ thuật Điện tử - Viễn thông. Tài liệu cung cấp phương pháp luận mô hình hóa FEM 2D cảm biến áp điện trên Ansys và kỹ thuật xử lý tín hiệu giải điều chế I/Q, giải ma trận phi tuyến trên Matlab/Simulink.
  • Kỹ sư nghiên cứu và phát triển (R&D): Làm việc trong lĩnh vực vi cơ điện tử (MEMS), thiết kế linh kiện sóng âm bề mặt (bộ lọc, cảm biến SAW), và các hệ thống thẻ nhận dạng RFID thụ động.
  • Các nhà phát triển hệ thống giám sát công nghiệp: Tham khảo cấu trúc mạng cảm biến không dây thụ động để xây dựng giải pháp giám sát nhiệt độ các điểm tiếp xúc phát nhiệt trong tủ điện kín, trạm biến áp cao thế hoặc giám sát trạng thái cơ học của các chi tiết quay trong công nghiệp nặng.

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

1. Tại sao tác giả lựa chọn vật liệu đế là $YZ\text{-LiNbO}_3$ thay vì Thạch anh (Quartz) hay các vật liệu áp điện khác?

Theo số liệu đối sánh tại Bảng 2.2 của luận án, $YZ\text{-LiNbO}_3$ sở hữu hệ số biến đổi điện cơ $K^2$ đạt 4,6%, vượt trội so với Thạch anh ($K^2$ chỉ từ 0% đến 0,16%), Langasite LGS (0,37%) hay $\text{LiTaO}_3$ (0,74%). Hệ số $K^2$ lớn giúp quá trình chuyển đổi năng lượng giữa điện trường và biến dạng cơ học diễn ra hiệu quả, tạo ra sóng âm bề mặt có biên độ lớn, rất quan trọng đối với cảm biến thụ động vốn dĩ không có nguồn pin hỗ trợ. Ngoài ra, vận tốc truyền sóng âm bề mặt lý thuyết của $YZ\text{-LiNbO}_3$ đạt $3480\text{ m/s}$ và hệ số nhiệt độ $\text{TCD} = 94\text{ ppm}/^\circ\text{C}$ phù hợp cho các cấu trúc cảm biến đo lường độ nhạy cao.

2. Thuật toán của luận án giải quyết hiện tượng nhập nhằng pha khi góc trễ pha lớn hơn $2\pi$ bằng nguyên lý nào?

Khi góc trễ pha vượt quá $2\pi$, một giá trị hàm lượng giác ($\arctan$) sẽ tương ứng với vô số góc pha thực tế $\varphi + 2n\pi$. Luận án giải quyết vấn đề này bằng cách thiết kế cấu trúc cảm biến gồm 3 bộ phản xạ với khoảng cách hình học xác định trên bề mặt đế. Bằng cách đo đạc hai hiệu pha $\Delta\varphi_{21} = \varphi_2 - \varphi_1$ và $\Delta\varphi_{31} = \varphi_3 - \varphi_1$, kết hợp với tỷ số trễ thời gian thô giữa các đỉnh phản xạ, thuật toán thiết lập phương trình liên hệ để tính toán chính xác số nguyên chu kỳ quay pha $n$. Từ đó, góc pha tuyệt đối được khôi phục chính xác mà không làm suy giảm tỷ số tín hiệu trên nhiễu (SNR).

3. Thuật toán đọc đồng thời tín hiệu trong hệ đa cảm biến SAW OFC ở Chương 3 đã khắc phục hạn chế gì của nghiên cứu Weifeng Liu [34]?

Nghiên cứu của Weifeng Liu và cộng sự [34] đã thiết lập hệ phương trình ma trận biểu diễn mối quan hệ giữa tín hiệu chồng lấn và tín hiệu thành phần nhưng chỉ giải quyết thành công ở tần số thấp ($12,5\text{ MHz}$). Khi tần số hoạt động nâng lên mức hàng trăm MHz, nghiệm của các hàm lượng giác biến thiên rất nhanh, khiến việc giải phương trình ma trận trực tiếp bị phân kỳ và cho nghiệm sai. Luận án của Nguyễn Thu Hà đã khắc phục triệt để hạn chế này bằng cách đề xuất thuật toán chia khoảng theo miền tần số/thời gian kết hợp điều kiện lọc nghiệm vật lý. Nhờ đó, thuật toán cho phép bóc tách chính xác biên độ và góc pha của 10 cảm biến SAW mã hóa OFC với số tần số $N_c = 4$ ở tần số hoạt động lên tới 500 MHz.

4. Tỉ số hóa kim $a/p$ tối ưu được xác định trong luận án là bao nhiêu và có ý nghĩa gì đối với công nghệ chế tạo?

Qua khảo sát mô phỏng FEM trên Ansys ở cả hai tần số 870 MHz và 125 MHz, tác giả kết luận tỉ số hóa kim tối ưu là $a/p = 0,5$ (trong đó $a$ là khoảng cách giữa hai điện cực/thanh phản xạ, $p$ là chu kỳ lặp lại). Giá trị $a/p = 0,5$ tương ứng với độ rộng ngón tay điện cực bằng khoảng cách giữa các ngón tay ($d = a = \lambda/4$). Kết cấu này giúp hệ thống đạt điều kiện phản xạ Bragg lý tưởng, cho biên độ sóng phản hồi cực đại và đồng thời thuận lợi cho công nghệ quang khắc chế tạo vi cơ điện tử vì kích thước các đường kim loại đồng đều nhau.

Kết luận

Luận án tiến sĩ của Nguyễn Thu Hà đã giải quyết các vấn đề then chốt trong thiết kế cấu trúc và xử lý tín hiệu cho cảm biến sóng âm bề mặt thụ động không dây loại đường trễ phản xạ (SAW RDL). Bằng việc kết hợp phương pháp mô phỏng phần tử hữu hạn FEM trên Ansys và tính toán số trên Matlab/Simulink, công trình đã xác lập các thông số hình học tối ưu để tối đa hóa biên độ sóng phản xạ và đề xuất phương pháp đo pha chính xác khử được nhập nhằng pha khi vượt ngưỡng $2\pi$. Đóng góp trọng tâm của luận án là thuật toán chia khoảng và lọc nghiệm giải hệ phương trình ma trận tín hiệu chồng lấn, cho phép nhận dạng và bóc tách đồng thời tín hiệu từ 10 cảm biến mã hóa tần số trực giao (OFC) ở dải tần số cao 500 MHz, tạo cơ sở kỹ thuật cho việc ứng dụng mạng cảm biến SAW thụ động trong quan trắc công nghiệp.