Tổng quan luận án

Thông tin thủy âm (Underwater Acoustic - UWA) là kỹ thuật truyền thông đóng vai trò then chốt trong nhiều ứng dụng dân sự và quân sự như quan trắc môi trường biển, thăm dò khai thác dầu khí ngoài khơi, cảnh báo thiên tai và các nhiệm vụ quốc phòng dưới nước. Tuy nhiên, môi trường truyền sóng âm trong môi trường nước nông (độ sâu dưới 200 m, đặc biệt là các vùng nước rất nông dưới 30 m) đặt ra những thách thức vật lý phức tạp:

  • Vận tốc lan truyền sóng âm trong nước rất thấp (xấp xỉ $1500\text{ m/s}$), dẫn đến thời gian trễ lan truyền lớn (hàng chục mili-giây cho cự ly trên 10 mét) và trải trễ đa đường kéo dài làm giảm băng thông kết hợp (coherence bandwidth), biến kênh truyền thủy âm thành kênh băng rộng chọn lọc tần số cao.
  • Suy hao lan truyền lớn và tăng mạnh theo tần số sóng mang do sự kết hợp của suy hao tán xạ hình học, suy hao phân kỳ và suy hao hấp thụ môi trường.
  • Tạp âm môi trường biển không phải là tạp âm trắng Gauss (AWGN) mà là tạp âm màu phụ thuộc tần số với mật độ phổ công suất suy giảm xấp xỉ $18\text{ dB/decade}$.
  • Hiện tượng đa đường biến đổi nhanh theo thời gian do phản xạ mạnh từ mặt biển chuyển động sóng và đáy biển gồ ghề, kết hợp với hiệu ứng dịch tần Doppler rất lớn (tỉ số Doppler tương đối $a = v/c_s$ đáng kể do vận tốc âm bé) gây ra bởi chuyển động của tàu mang đầu phát/thu hoặc sóng mặt nước.

Khoảng trống nghiên cứu lớn trong lĩnh vực này là:

  1. Thiếu một mô hình kênh truyền chuẩn hóa áp dụng chung cho mọi vùng nước do sự khác biệt về cấu trúc địa lý, chu kỳ mùa và thời tiết.
  2. Các mô hình kênh truyền dựa trên hình học (geometry-based) tuy có chi phí tính toán thấp nhưng cho kết quả thống kê sai lệch nhiều so với kênh đo đạc thực tế do giả định phân bố tán xạ không phản ánh đúng môi trường nước nông.
  3. Các mô hình dựa trên đo đạc thực nghiệm (measurement-based) dù phản ánh sát thực tế nhưng lại đòi hỏi quy trình tối ưu hóa phi tuyến phức tạp (như LPNM, ESPRIT, SAGE, INLSA) để ước lượng số lượng lớn tham số kênh.
  4. Việc mô hình hóa phổ công suất Doppler dạng biểu thức giải tích đóng (closed-form) còn thiếu vắng trong các y văn, đặc biệt là việc tách bạch thành phần Doppler do chuyển động của máy phát/máy thu (Tx/Rx) và thành phần do dao động mặt biển.
  5. Các nghiên cứu phân tích hệ thống ghép kênh phân chia theo tần số trực giao (OFDM) dưới nước thường giả định phổ Doppler kinh điển hoặc tạp âm trắng, chưa đánh giá đồng thời tác động của nhiễu xuyên kênh (ICI) trên mô hình kênh thực nghiệm và tạp âm màu môi trường.

Luận án tiến sĩ kỹ thuật viễn thông của tác giả Đỗ Việt Hà (mã số chuyên ngành: 62520208), thực hiện tại Trường Đại học Bách khoa Hà Nội dưới sự hướng dẫn của PGS. Nguyễn Văn Đức và TS. Phạm Văn Tiến, hướng tới ba mục tiêu nghiên cứu cụ thể:

  1. Nghiên cứu, so sánh và đề xuất phương pháp thiết kế bộ mô phỏng kênh truyền thủy âm nước nông (cho cả trường hợp Tx/Rx cố định và Rx di động) có độ chính xác cao, khớp với dữ liệu đo đạc thực tế nhưng rút gọn tối đa độ phức tạp tính toán tham số.
  2. Xây dựng nền tảng lý thuyết và thiết lập biểu thức giải tích dạng đóng cho phổ công suất Doppler của kênh thủy âm nước nông có xét đồng thời chuyển động của bề mặt biển và chuyển động tương đối của trạm thu/phát, được kiểm chứng bằng thực nghiệm.
  3. Sử dụng mô hình kênh truyền xây dựng từ dữ liệu đo đạc thực tế để phân tích giải tích chính xác hiệu năng hệ thống UWA-OFDM (nhiễu ICI, tạp âm môi trường, công suất phát, tỉ số SIR, SINR, dung lượng kênh và hiệu suất phổ), cung cấp các nguyên tắc lựa chọn tham số truyền dẫn tối ưu.

Đối tượng nghiên cứu của luận án là đặc tính truyền sóng, hàm tương quan thời gian - tần số, hàm truyền đạt, phổ công suất Doppler và hiệu năng truyền dẫn OFDM trong môi trường thủy âm nước nông. Phạm vi thực nghiệm của luận án tập trung vào vùng nước nông tại Vịnh Hạ Long, Việt Nam, với các dữ liệu đo đạc thực tế thu thập vào ngày 09/06/2015 ở cự ly truyền dẫn $150\text{ m}$ và tần số sóng mang $12\text{ kHz}$.


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

Trong luận án, tác giả đã tổng hợp và phân loại các hướng tiếp cận mô hình hóa kênh truyền thủy âm trên thế giới thành hai nhóm chính:

  1. Nhóm mô hình dựa trên hình học (Geometry-based modeling):

    • Tiêu biểu là công cụ mô phỏng Bellhop [69] dựa trên lý thuyết tia cho một hình học xác định. Để phản ánh tính ngẫu nhiên, một số nghiên cứu đã kết hợp Bellhop với các điều kiện môi trường như nhiệt độ và độ mặn [90], tốc độ gió [28], hình dạng mặt sóng [38].
    • Hướng tiếp cận hình học kết hợp phương pháp thống kê [14, 18, 27, 55, 73, 75, 105, 106] mô tả kênh truyền thông qua hàm mật độ xác suất (PDF) của góc tới (AOA) và góc xuất phát (AOD). Các giả định phân bố AOA thường dùng gồm phân bố chuẩn [105, 106], phân bố đều [75] hoặc phân bố nửa vòng tròn Rice [14]. Ưu điểm của mô hình hình học là cấu trúc đơn giản, ít tham số và dễ dàng mở rộng sang các điều kiện độ sâu hay cự ly khác nhau, nhưng nhược điểm lớn là các đặc tính thống kê không khớp với các phép đo thực tế do sự biến động không gian - thời gian phức tạp của môi trường nước nông.
  2. Nhóm mô hình dựa trên đo đạc thực nghiệm (Measurement-based modeling):

    • Được nghiên cứu qua các công trình [24, 74, 76, 86, 107]. Dựa trên phân tích dữ liệu đo đạc tại hiện trường, các tác giả xác định phân bố của các tia sóng theo phân bố Rayleigh [24, 86], Rician [76], phân bố K [107], hoặc log-normal [74].
    • Các bộ mô phỏng dạng phát lại (replay-based simulators) [58, 84, 97] tái tạo đáp ứng xung kênh biến đổi theo thời gian (TVCIR) từ dữ liệu đo hoặc sinh TVCIR ngẫu nhiên có cùng đặc tính thống kê. Để trích xuất các tham số (độ lợi tia, dịch Doppler, độ trễ, lệch pha), các thuật toán ước lượng cao cấp thường được sử dụng như ESPRIT [35], SAGE [34], xấp xỉ bình phương tối thiểu phi tuyến lặp INLSA [32], và phương pháp chuẩn $L_p$ (LPNM) [59]. Hạn chế của hướng này là chi phí đo đạc lớn, thuật toán tối ưu hóa phức tạp và mô hình chỉ có giá trị cho chính môi trường đo đạc.
  3. Tình hình nghiên cứu tại Việt Nam:

    • Các nghiên cứu về thông tin thủy âm trong nước còn rất hạn chế [2, 3, 6]. Một số nghiên cứu trước đây bước đầu khảo sát đặc tính lan truyền âm [1, 4, 5, 7, 8]. Trong công trình [6], các tác giả mô phỏng tia truyền sóng thủy âm bằng cách giải phương trình Eikonal cho các điều kiện môi trường định trước, tuy nhiên mô hình này có tính chất bất biến theo thời gian (time-invariant) và chưa đưa ra được mô hình đặc tính thống kê kênh truyền hoàn chỉnh.

Vị trí của luận án: Luận án giải quyết trực tiếp bài toán cân bằng giữa độ chính xác thực tế và độ phức tạp tính toán bằng cách đề xuất phương pháp trích xuất tham số trực tiếp từ biểu đồ trễ công suất (PDP) đo đạc mà không cần chạy các thuật toán tối ưu hóa phi tuyến phức tạp; đồng thời xây dựng mô hình phổ Doppler dạng đóng kết hợp hai thành phần và hoàn thiện khung phân tích ICI cho hệ thống UWA-OFDM.


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

Cơ sở lý thuyết

Luận án dựa trên các nền tảng lý thuyết truyền thông vô tuyến và thủy âm học:

  1. Lý thuyết quá trình ngẫu nhiên và mô hình kênh fading:

    • Quá trình Gauss phức $\mu(t) = \mu_1(t) + j\mu_2(t)$ mô tả kênh fading băng hẹp.
    • Mô hình Rice và Rayleigh: biểu diễn quá trình Gauss thực $\mu_i(t)$ qua tổng vô hạn các hàm điều hòa: $$\mu_i(t) = \lim_{N_i \to \infty} \sum_{n=1}^{N_i} c_{i,n} \cos(2\pi f_{i,n} t + \theta_{i,n})$$ Khi có thành phần trực diện (LOS) $m(t) = \rho e^{j(2\pi f_p t + \theta_p)}$, kênh là quá trình Rice $\xi(t) = |\mu(t) + m(t)|$. Khi không có LOS, kênh trở thành quá trình Rayleigh $\zeta(t) = |\mu_1(t) + j\mu_2(t)|$.
    • Nguyên lý mô hình hóa kênh xác định (deterministic channel modeling): xấp xỉ quá trình ngẫu nhiên bằng số lượng hữu hạn các thành phần điều hòa.
    • Phương pháp Tổng các hàm Sin (Sum-of-Sinusoids - SOS) của Rice [15, 54, 62]: giả định thành phần đồng pha $\mu_1(t)$ và vuông pha $\mu_2(t)$ không tương quan, chỉ mô phỏng được phổ mật độ công suất Doppler (DPSD) đối xứng (tán xạ đẳng hướng).
    • Phương pháp Tổng các hàm số mũ phức (Sum-of-Cisoids - SOC) [79] và mô hình SOC tán xạ không tương quan (SOCUS) [39, 56, 62]: $$\tilde{\mu}(t) = \sum_{n=1}^N c_n e^{j(2\pi f_n t + \theta_n)}$$ Mô hình SOC cho phép mô phỏng cả DPSD đối xứng và bất đối xứng, có cấu trúc tương thích với bản chất lan truyền sóng vật lý, rất phù hợp cho kênh UWA.
  2. Hệ thống hàm đặc tính kênh truyền:

    • Đáp ứng xung kênh biến đổi theo thời gian (TVCIR) $h(\tau, t)$.
    • Hàm truyền đạt kênh biến đổi theo thời gian (TVCTF) $H(f, t) = \int h(\tau, t) e^{-j2\pi f \tau} d\tau$.
    • Hàm tương quan thời gian - tần số (T-FCF): $$R_{HH}(\Delta f, \Delta t) = E{H(f, t) H^*(f + \Delta f, t + \Delta t)}$$
    • Hàm tương quan tần số (FCF) $R_{HH}(\Delta f) = R_{HH}(\Delta f, 0)$ biểu diễn tính chọn lọc tần số.
    • Hàm tương quan thời gian (TCF) $R_{HH}(\Delta t) = R_{HH}(0, \Delta t)$ biểu diễn mức độ biến đổi theo thời gian.
    • Biểu đồ trễ công suất (PDP) $\rho(\tau)$ thu được từ trung bình thống kê $|\hat{h}(\tau, t)|^2$.
  3. Mô hình suy hao lan truyền và hấp thụ trong môi trường nước nông:

    • Suy hao tổng thể: $A(l, f) = \alpha(f)^{l - l_r} (l/l_r)^k$, với $k$ là số mũ suy hao lan truyền ($1 \le k \le 2$).
    • Hệ số suy hao hấp thụ $A_a(D) = 10^{-\frac{D\beta}{20000}}$, trong đó $\beta$ được tính theo công thức thực nghiệm của Francois-Garrison [20] phụ thuộc vào độ mặn $S_a$, tần số sóng mang $f_c$, áp suất thủy tĩnh $P = 1.013(1 + 0.1h)$ và độ sâu $h$.

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

Luận án kết hợp phương pháp giải tích lý thuyết, mô phỏng số và thực nghiệm tại hiện trường:

Hạng mục Phương pháp và công cụ cụ thể
Mô hình hình học Phương pháp Method of Equally Spaced Scatterers (MESS) xác định tọa độ các điểm tán xạ trên mặt và đáy biển.
Mô hình đo đạc truyền thống Tối ưu hóa phi tuyến $L_p$-Norm Method (LPNM) với chuẩn bậc $p = 2$ để khớp hàm FCF/TCF của mô hình mô phỏng với dữ liệu đo.
Phương pháp đề xuất Trích xuất trực tiếp độ lợi đường $c_{n,l} = \sqrt{\hat{\rho}(\tau_l)}$ và trễ truyền dẫn $\tau_l$ từ PDP đo đạc thực nghiệm.
Hệ thống thực nghiệm đo đạc Triển khai tại Vịnh Hạ Long (09/06/2015), cự ly $150\text{ m}$. Tín hiệu dò kênh là chuỗi PN bậc 10 sinh bởi đa thức $G(x) = x^{10} + x^7 + 1$, điều chế BPSK, phát qua transducer. Đầu thu dùng hydrophone, khuếch đại LNA, giải điều chế BPSK và bộ dò tương quan để thu 96 đáp ứng xung kênh $\hat{h}(\tau_l, t_m)$.
Mô hình hóa Doppler Kỹ thuật khớp đường cong (curve fitting) tối ưu hóa các tham số của mô hình kết hợp Spike - Gaussian trên dữ liệu đo 3 kịch bản chuyển động.
Phân tích hệ thống Giải tích tích phân số chính xác công suất nhiễu ICI và công suất tạp âm môi trường màu trên từng sóng mang con của hệ thống UWA-OFDM.

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

Chapter 1: DESIGN OF SHALLOW UWA CHANNEL SIMULATORS (Thiết kế các bộ mô phỏng kênh truyền thủy âm nước nông)

Chương 1 tập trung phân tích, đối sánh hai phương pháp thiết kế bộ mô phỏng kênh UWA (dựa trên hình học và dựa trên đo đạc), từ đó đề xuất phương pháp thiết kế mới cho cả kênh tĩnh và kênh có hiệu ứng Doppler:

  • Mô hình bộ mô phỏng dựa trên hình học:

    • Thiết lập mô hình hình học 2D/3D cho vùng nước nông với độ sâu $h$, khoảng cách $D$, các điểm tán xạ ngẫu nhiên $S_{i,n}$ trên mặt biển ($i=1$) và đáy biển ($i=2$).
    • Đáp ứng xung kênh gồm 3 thành phần: đường trực tiếp $h_0(\tau, t)$, tán xạ mặt $h_1(\tau, t)$, và tán xạ đáy $h_2(\tau, t)$.
    • Tác giả sử dụng phương pháp MESS tính tọa độ tối ưu của các điểm tán xạ: $$x_{i,n}^{\text{opt}} = \frac{D}{N_i}\left(n - \frac{1}{2}\right)$$
    • Áp dụng các tham số hình học mô phỏng: thừa số Rice $c_R = 0$, độ cao/sâu trạm $y_{1}^T = y_{2}^T = y_{1}^R = y_{2}^R = 90\text{ m}$, số tia phản xạ mặt $N_1 = 80$, số tia phản xạ đáy $N_2 = 79$, tần số sóng mang $f_c = 12\text{ kHz}$, cự ly $D = 500\text{ m}$. Kết quả cho thấy FCF của mô hình mô phỏng khớp với mô hình tham chiếu hình học, nhưng khi so sánh với FCF trích xuất từ dữ liệu đo đạc thực tế tại Vịnh Hạ Long, mô hình hình học có sai số rất lớn do các giả định tán xạ đều không phù hợp với thực tế biển đảo.
  • Mô hình bộ mô phỏng dựa trên đo đạc truyền thống:

    • Sử dụng cấu trúc SOCUS với $L$ cụm đường trễ. Áp dụng thuật toán LPNM để tối thiểu hóa hàm chi phí sai số năng lượng FCF/TCF (phương trình chuẩn bậc $p=2$). Kết quả FCF mô phỏng khớp rất chính xác với FCF tham chiếu từ dữ liệu đo đạc, nhưng thuật toán đòi hỏi khối lượng tính toán lớn để tối ưu hóa đồng thời $N_l, c_{n,l}, \tau_l, f_{n,l}$.
  • Phương pháp đề xuất cho kênh UWA tĩnh (Static Channel):

    • Tác giả đề xuất quy trình 3 bước:
      1. Xác định biểu đồ trễ công suất $\hat{\rho}(\tau_l)$ từ 96 mẫu đo CIR thực tế: $\hat{\rho}(\tau_l) = \frac{1}{M}\sum_{m=1}^M |\hat{h}[\tau_l, t_m]|^2$.
      2. Chọn mô hình mô phỏng dạng SOCUS với $f_{n,l} = 0$.
      3. Trích xuất trực tiếp độ trễ $\tau_l$ và độ lợi đường: $$c_{n,l} = \sqrt{\frac{\hat{\rho}(\tau_l)}{N_l}}$$ (với $N_l = 1$, $c_{1,l} = \sqrt{\hat{\rho}(\tau_l)}$).
    • Kết quả so sánh cho thấy hàm FCF của bộ mô phỏng đề xuất trùng khít hoàn toàn với FCF của mô hình đo đạc thực tế mà không cần chạy bất kỳ thuật toán tối ưu hóa phức tạp nào.
  • Phương pháp đề xuất cho trường hợp có hiệu ứng Doppler (Moving Rx):

    • Đối với máy thu chuyển động, tác giả thiết lập hàm tương quan thời gian - tần số T-FCF $R_{HH}(\Delta f, \Delta t)$. Mô hình đề xuất tiếp tục tận dụng các tham số trích xuất trực tiếp từ đo đạc, cho sai số T-FCF rất nhỏ so với mô hình tham chiếu đo đạc và vượt trội so với các mô hình kinh điển.

Chapter 2: MODELING OF DOPPLER POWER SPECTRUM FOR SHALLOW UWA CHANNELS (Mô hình hóa phổ công suất Doppler cho kênh truyền thủy âm nước nông)

Chương 2 giải quyết bài toán mô hình hóa phổ công suất Doppler trong môi trường nước nông có sự kết hợp của hai nguồn gây dịch Doppler: chuyển động cơ học của tàu/đầu dò (Tx/Rx) và chuyển động sóng bề mặt biển:

  • Phân tích lý thuyết thành phần Doppler:

    • Dịch Doppler do chuyển động của trạm thu với vận tốc $V_R$ tạo ra phổ dạng đỉnh nhọn (Spike-shape Doppler spectrum) tại tần số dịch chuyển cực đại $f_D = \frac{V_R f_c}{c_s} \cos(\alpha - \alpha_v^R)$.
    • Dao động ngẫu nhiên của sóng mặt biển tạo ra phổ dạng Gauss (Gaussian-shape Doppler spectrum) do tính chất phân bố thống kê ngẫu nhiên của độ dốc và vận tốc hạt nước trên mặt sóng.
  • Đề xuất mô hình phổ Doppler dạng đóng:

    • Tác giả đề xuất mô hình phổ công suất Doppler tổng thể $S(f)$ dưới dạng tổng của thành phần Spike-shape $S_{\text{spike}}(f)$ và thành phần Gaussian-shape $S_{\text{Gauss}}(f)$: $$S(f) = S_{\text{spike}}(f) + S_{\text{Gauss}}(f)$$
    • Biểu thức giải tích này phản ánh đầy đủ tác động tương hỗ giữa vận tốc di chuyển của thiết bị và điều kiện động lực học mặt nước.
  • Kiểm chứng thực nghiệm tại Vịnh Hạ Long:

    • Thiết lập 3 kịch bản đo đạc Doppler (Scenario 1, Scenario 2, Scenario 3) với các trạng thái chuyển động của tàu thu (trôi tự do theo sóng, chạy với vận tốc ổn định $V_R$, hoặc di chuyển theo quỹ đạo đổi hướng) và các mức độ sóng biển khác nhau.
    • Áp dụng thuật toán tính toán và tối ưu hóa tham số để khớp mô hình đề xuất $S(f)$ với phổ đo đạc tham chiếu $\tilde{S}_n(f)$. Kết quả cho thấy mô hình giải tích đề xuất bám sát hoàn toàn dữ liệu phổ thực nghiệm ở cả 3 kịch bản, xác thực tính đúng đắn của việc kết hợp thành phần Spike và Gaussian.

Chapter 3: UWA-OFDM SYSTEM PERFORMANCE ANALYSIS USING THE MEASUREMENT-BASED UWA CHANNEL MODEL (Phân tích hiệu năng hệ thống UWA-OFDM sử dụng mô hình kênh truyền thủy âm dựa trên đo đạc)

Chương 3 ứng dụng mô hình kênh đo đạc thực tế để đánh giá toàn diện hiệu năng của hệ thống thông tin băng rộng UWA-OFDM dưới tác động đồng thời của nhiễu xuyên kênh (ICI) do Doppler và tạp âm môi trường màu:

  • Phân tích giải tích nhiễu ICI và tạp âm môi trường:

    • Hệ thống OFDM chia băng thông $B$ thành $N$ sóng mang con với khoảng cách sóng mang $\Delta f = B/N$, chu kỳ ký hiệu $T = 1/\Delta f$. Trải trễ lớn được triệt tiêu bằng khoảng bảo vệ (CP), nhưng tính biến đổi theo thời gian và hiệu ứng Doppler làm mất tính trực giao giữa các sóng mang con, gây ra nhiễu ICI nghiêm trọng.
    • Tạp âm môi trường biển được mô hình hóa theo quy luật hàm mật độ phổ công suất suy giảm $18\text{ dB/decade}$ theo tần số, tính toán công suất tạp âm chính xác trên từng băng con của sóng mang.
  • Đánh giá hiệu năng và khảo sát tham số:

    • Luận án thực hiện tính toán chính xác công suất tín hiệu có ích, công suất nhiễu ICI và công suất tạp âm môi trường để xác định các chỉ số:
      • Tỉ số tín hiệu trên can nhiễu (SIR - Signal-to-Interference Ratio).
      • Tỉ số tín hiệu trên can nhiễu cộng tạp âm (SINR - Signal-to-Interference-plus-Noise Ratio).
      • Dung lượng kênh truyền (Capacity) và hiệu suất phổ trung bình (Average Spectral Efficiency).
    • Khảo sát các trường hợp thay đổi số lượng sóng mang con ($N = 256,...$), băng thông tín hiệu ($B = 10\text{ kHz}$), công suất phát $P_T$ và SNR tại đầu thu.
    • Xây dựng hàm phân bố tích lũy (CDF) của SNR, SIR và SINR trên các sóng mang con.
    • Kết quả phân tích chỉ ra sự đánh đổi (trade-off) quan trọng: tăng số sóng mang con $N$ giúp giảm nhiễu xuyên ký hiệu (ISI) nhưng làm tăng mạnh nhiễu ICI do thu hẹp khoảng cách sóng mang $\Delta f$; đồng thời đưa ra các đồ thị hướng dẫn lựa chọn tham số $B, N, P_T$ tối ưu để đạt dung lượng kênh mong muốn với mức công suất phát nhỏ nhất.

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

Luận án đã đạt được các kết quả khoa học và đóng góp mới nổi bật sau:

  1. Phương pháp thiết kế bộ mô phỏng kênh UWA nước nông hiệu quả tính toán:

    • Đã làm rõ các hạn chế căn bản của mô hình hình học (kém chính xác so với thực tế) và mô hình đo đạc truyền thống (tính toán tối ưu hóa phức tạp).
    • Đề xuất phương pháp thiết kế bộ mô phỏng kênh SOCUS trực tiếp từ biểu đồ trễ công suất (PDP) của dữ liệu đo đạc thực nghiệm cho cả kênh tĩnh và kênh di động có Doppler. Phương pháp này loại bỏ hoàn toàn các thuật toán tối ưu hóa phi tuyến phức tạp (như LPNM, ESPRIT), giảm thiểu tối đa khối lượng tính toán nhưng vẫn đảm bảo các hàm đặc tính thống kê (FCF, T-FCF) hoàn toàn trùng khớp với kênh thực tế.
  2. Mô hình giải tích dạng đóng cho phổ công suất Doppler vùng nước nông:

    • Xây dựng thành công cơ sở lý thuyết phân tách hai nguồn gốc vật lý của hiệu ứng Doppler trong môi trường thủy âm nước nông.
    • Thiết lập biểu thức giải tích dạng đóng cho phổ công suất Doppler kết hợp giữa dạng Spike (đặc trưng cho chuyển động máy phát/máy thu) và dạng Gaussian (đặc trưng cho chuyển động mặt nước). Mô hình đã được kiểm chứng và khớp chính xác với dữ liệu đo phổ Doppler thực nghiệm tại Vịnh Hạ Long qua 3 kịch bản chuyển động khác nhau.
  3. Khung phân tích chính xác hiệu năng hệ thống UWA-OFDM:

    • Cung cấp mô hình giải tích định lượng chính xác công suất nhiễu ICI kết hợp với tạp âm môi trường màu (suy giảm theo tần số) trên nền kênh đo đạc thực tế.
    • Đánh giá toàn diện các chỉ số SIR, SINR, dung lượng kênh và hiệu suất phổ theo các tham số hệ thống ($B, N, P_T$).
    • Đưa ra bộ chỉ dẫn thực tiễn (practical guidelines) cho việc thiết kế và tối ưu hóa tham số truyền dẫn của modem UWA-OFDM hoạt động trong môi trường nước nông.

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

Từ các nội dung trình bày trong văn bản luận án, các hạn chế và hướng phát triển tiếp theo được xác định như sau:

  • Hạn chế:

    • Phạm vi thực nghiệm của luận án tập trung tại một vùng nước nông cụ thể (Vịnh Hạ Long, Việt Nam) với cự ly đo đạc lớn nhất của hệ thống thử nghiệm là $150\text{ m}$.
    • Mô hình kênh đo đạc có tính đặc thù cho khu vực khảo sát, chưa bao quát các vùng nước sâu hoặc các vùng biển có cấu trúc địa tầng và điều kiện thủy văn biến đổi theo các mùa khác nhau trong năm.
    • Luận án chưa phân tích chi tiết sự phụ thuộc của suy hao đường truyền và tần số sóng mang theo nhiều khoảng cách truyền dẫn khác nhau do giới hạn của hệ thống đo đạc phần cứng.
  • Hướng nghiên cứu tiếp:

    • Mở rộng chiến dịch đo đạc thu thập dữ liệu tại nhiều vùng biển, nhiều cự ly truyền dẫn và các điều kiện thời tiết, độ sâu khác nhau nhằm xây dựng cơ sở dữ liệu kênh truyền phong phú hơn.
    • Nghiên cứu kết hợp mô hình hình học với dữ liệu thực nghiệm để xây dựng các bộ mô phỏng lai (hybrid simulators) có khả năng thích ứng linh hoạt với việc thay đổi tham số môi trường.
    • Ứng dụng mô hình kênh và kết quả phân tích ICI để thiết kế các thuật toán ước lượng kênh, đồng bộ hóa và cân bằng triệt nhiễu ICI tiên tiến cho hệ thống UWA-OFDM.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu có giá trị cho nhiều đối tượng trong lĩnh vực kỹ thuật vô tuyến và thông tin dưới nước:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Kỹ thuật Viễn thông, Xử lý tín hiệu và Kỹ thuật Thủy âm:

    • Tham khảo phương pháp luận mô hình hóa kênh truyền ngẫu nhiên, phương pháp Sum-of-Cisoids (SOC/SOCUS), quy trình xử lý tín hiệu đo đạc CIR bằng chuỗi PN và kỹ thuật phân tích hàm tương quan thời gian - tần số.
    • Sử dụng biểu thức giải tích phổ Doppler dạng đóng Spike - Gaussian để phục vụ các bài toán mô phỏng kênh truyền UWA.
  • Kỹ sư và chuyên gia phát triển hệ thống thông tin dưới nước (AUV, modem thủy âm, mạng cảm biến ngầm UWSN):

    • Áp dụng trực tiếp phương pháp thiết kế bộ mô phỏng kênh nhanh từ dữ liệu PDP để kiểm thử thuật toán trước khi triển khai phần cứng.
    • Sử dụng các kết quả phân tích và đồ thị đánh đổi tham số trong Chương 3 để lựa chọn băng thông, số sóng mang con OFDM và tính toán mức công suất phát cần thiết nhằm đảm bảo tỉ số SINR và dung lượng kênh mong muốn.

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

1. Tại sao mô hình Sum-of-Cisoids (SOC) lại phù hợp hơn mô hình Sum-of-Sinusoids (SOS) trong mô phỏng kênh thủy âm nước nông?

Mô hình SOS của Rice giả định hai thành phần đồng pha và vuông pha không tương quan, do đó chỉ mô phỏng được các kênh có phổ mật độ công suất Doppler đối xứng dưới điều kiện tán xạ đẳng hướng. Trong khi đó, kênh thủy âm nước nông có tính tán xạ phi đẳng hướng mạnh với phổ Doppler bất đối xứng. Mô hình SOC sử dụng tổng các hàm số mũ phức (cisoids) với cùng một tập tham số cho cả hai thành phần, có cấu trúc tương đồng với sự lan truyền sóng vật lý thực tế, cho phép mô phỏng chính xác cả phổ Doppler đối xứng lẫn bất đối xứng.

2. Cấu hình hệ thống thực nghiệm đo đạc đáp ứng xung kênh tại Vịnh Hạ Long được thiết lập như thế nào?

Thực nghiệm được thực hiện ngày 09/06/2015 ở cự ly $150\text{ m}$. Phía phát đặt trên tàu A gồm máy tính, bộ khuếch đại công suất và transducer phát tín hiệu chuỗi PN bậc 10 (đa thức sinh $G(x) = x^{10} + x^7 + 1$) điều chế BPSK. Phía thu đặt trên thuyền thu gồm hydrophone, bộ khuếch đại tạp âm thấp (LNA), máy tính giải điều chế BPSK và bộ tách tham số tương quan. Dữ liệu thu thập gồm 96 lần đo CIR liên tục để trích xuất $\hat{h}(\tau_l, t_m)$ và tính biểu đồ trễ công suất $\hat{\rho}(\tau_l)$.

3. Điểm cải tiến cốt lõi của phương pháp thiết kế bộ mô phỏng kênh tĩnh đề xuất so với phương pháp LPNM là gì?

Phương pháp LPNM truyền thống đòi hỏi giải bài toán tối ưu hóa phi tuyến phức tạp để tìm các tham số kênh nhằm cực tiểu hóa hàm chi phí sai số FCF. Trong phương pháp đề xuất, tác giả chứng minh được rằng $\hat{\rho}(\tau_l) = \sum_{n=1}^{N_l} c_{n,l}^2$, từ đó trích xuất trực tiếp độ trễ $\tau_l$ từ vị trí đỉnh của PDP đo đạc và gán trực tiếp độ lợi đường $c_{n,l} = \sqrt{\hat{\rho}(\tau_l)/N_l}$ (với $N_l = 1$ thì $c_{1,l} = \sqrt{\hat{\rho}(\tau_l)}$). Phương pháp này cho hàm FCF trùng khớp với thực nghiệm nhưng triệt tiêu hoàn toàn chi phí tính toán tối ưu hóa.

4. Mô hình phổ công suất Doppler đề xuất trong luận án gồm những thành phần nào và mô tả hiện tượng vật lý gì?

Mô hình phổ Doppler đề xuất $S(f)$ là tổng của hai thành phần giải tích:

  • Thành phần dạng Spike ($S_{\text{spike}}(f)$): Mô tả dịch Doppler do chuyển động cơ học tương đối giữa trạm phát và trạm thu (Tx/Rx) với vận tốc $V_R$.
  • Thành phần dạng Gauss ($S_{\text{Gauss}}(f)$): Mô tả độ trải phổ Doppler do dao động ngẫu nhiên của sóng bề mặt biển gây ra khi sóng âm phản xạ qua mặt nước.

5. Khi phân tích hiệu năng hệ thống UWA-OFDM, luận án đã xem xét những yếu tố can nhiễu và tạp âm đặc thù nào?

Luận án xem xét đồng thời hai yếu tố đặc thù của môi trường thủy âm:

  1. Nhiễu xuyên kênh (ICI): Gây ra bởi hiệu ứng Doppler lớn làm mất tính trực giao giữa các sóng mang con khi chu kỳ ký hiệu kéo dài.
  2. Tạp âm môi trường màu: Tạp âm nền của biển (do sóng, gió, chuyển động tàu bè) có mật độ phổ công suất suy giảm $18\text{ dB/decade}$ theo tần số (không phải tạp âm trắng AWGN), được tích phân chính xác trên từng băng tần của sóng mang con.

Kết luận

Luận án tiến sĩ của tác giả Đỗ Việt Hà đã giải quyết toàn diện bài toán mô hình hóa kênh truyền thủy âm vùng nước nông và phân tích hiệu năng hệ thống UWA-OFDM dựa trên số liệu thực nghiệm tại vùng biển Việt Nam. Bằng việc đề xuất phương pháp trích xuất trực tiếp tham số từ biểu đồ trễ công suất và xây dựng biểu thức giải tích phổ Doppler dạng đóng kết hợp Spike - Gaussian, luận án đã dung hòa được yêu cầu về độ chính xác thống kê thực tế và sự đơn giản hóa trong tính toán mô phỏng. Các kết quả phân tích định lượng về tương quan giữa nhiễu ICI, tạp âm môi trường và cấu hình sóng mang cung cấp cơ sở khoa học tin cậy cho việc thiết kế, tối ưu hóa các hệ thống thông tin thủy âm băng rộng hiện đại.