Tổng quan về luận án

Môi trường khí quyển tầng đối lưu đặc trưng bởi tính không đồng nhất và biến động phức tạp theo không gian, thời gian, gây ra các hiện tượng bẻ cong tia sóng, biến đổi pha, tán xạ, hấp thụ và trễ đường truyền vô tuyến. Trong kỹ thuật viễn thông hiện đại, hệ thống thông tin vi ba, radar giám sát, dẫn đường hàng không và định vị vệ tinh toàn cầu (GNSS), việc xác định chính xác trạng thái khúc xạ của tầng đối lưu có ý nghĩa sống còn nhằm tối ưu hóa đường truyền và hiệu chỉnh sai số định vị. Luận án tiến sĩ kỹ thuật viễn thông "Nghiên cứu, đánh giá điều kiện truyền sóng vô tuyến tầng đối lưu khu vực Hà Nội sử dụng các phương pháp cắt lớp vô tuyến và bóng thám không" của nghiên cứu sinh Phạm Chí Công (Trường Đại học Giao thông Vận tải, 2023), dưới sự hướng dẫn khoa học của TS. Trần Hoài Trung và TS. Phạm Xuân Thành, là công trình tiên phong giải quyết bài toán ước lượng cấu trúc chiết suất khí quyển thực nghiệm tại Việt Nam.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án hướng tới là sự thiếu hụt các mô hình cấu trúc chỉ số khúc xạ thực nghiệm theo chiều thẳng đứng tại khu vực nhiệt đới gió mùa Bắc Bộ. Trong nhiều thập kỷ qua, các tính toán thiết kế tuyến vi ba và hiệu chỉnh trễ tầng đối lưu tại Việt Nam chủ yếu dựa vào mô hình khí quyển chuẩn toàn cầu dạng hàm số mũ theo khuyến nghị ITU-R P.453 ($N_0 = 315\text{ N-units}, h_0 = 7,35\text{ km}$), vốn được xây dựng từ dữ liệu của các vùng khí hậu ôn đới. Mô hình này bỏ qua các biến động cực đoan của độ ẩm nhiệt đới và hiện tượng nghịch nhiệt cục bộ, dẫn đến sai số lớn khi dự đoán quỹ đạo tia sóng và xác suất xuất hiện các cơ chế truyền sóng bất thường.

Luận án đặt ra ba câu hỏi nghiên cứu (Research Questions - RQ) và ba giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Quy luật biến đổi theo độ cao và theo mùa của độ khúc xạ vô tuyến $N$, độ dốc khúc xạ $G$ và hệ số bán kính Trái Đất hiệu dụng $k$ tại khu vực Hà Nội diễn ra như thế nào khi trích xuất từ dữ liệu viễn thám cắt lớp vô tuyến (Radio Occultation - RO)?
  • RQ2: Mức độ tương đồng, độ lệch chuẩn và sai số giữa phương pháp đo gián tiếp từ dữ liệu bóng thám không (Radiosonde) tại trạm khí tượng so với phương pháp cắt lớp vô tuyến vệ tinh LEO đạt ngưỡng nào?
  • RQ3: Mô hình chuẩn toàn cầu ITU-R P.453 có phản ánh chính xác cấu trúc tầng đối lưu khu vực Hà Nội hay không, và mức độ sai lệch tuyệt đối $\Delta N$ cũng như sai lệch tương đối $\Delta \varepsilon%$ là bao nhiêu?
  • H1: Chỉ số khúc xạ vô tuyến tầng đối lưu khu vực Hà Nội có sự phân hóa mùa sâu sắc, trong đó độ khúc xạ bề mặt mùa ẩm vượt xa giá trị chuẩn 315 N-units của ITU-R.
  • H2: Kết hợp dữ liệu cắt lớp vô tuyến vệ tinh COSMIC-1 và bóng thám không cho phép tái tạo chính xác cấu trúc thẳng đứng của chiết suất khí quyển với sai số nhỏ hơn các mô hình tham chiếu lý thuyết.
  • H3: Hiện tượng siêu khúc xạ và khúc xạ ống dẫn xuất hiện với tần suất đáng kể trong các lớp nghịch nhiệt sát mặt đất tại Hà Nội, đòi hỏi hiệu chỉnh hệ số $k$-factor thực nghiệm thay vì áp dụng giá trị mặc định $k = 4/3$.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết truyền sóng điện từ trong môi trường không đồng nhất, định luật Snell mở rộng cho môi trường cầu đối xứng, lý thuyết tán xạ Mie – Rayleigh, và các mô hình tính toán chiết suất Bean & Dutton (1968) được chuẩn hóa trong chuỗi khuyến nghị ITU-R P Series (P.453-14, P.834, P.310).

Phạm vi nghiên cứu không gian được giới hạn tại khu vực Hà Nội, xác định xung quanh tọa độ $21,01^\circ\text{N}$ vĩ độ Bắc và $105,80^\circ\text{E}$ kinh độ Đông trong bán kính $2^\circ$ kinh vĩ. Khung thời gian thu thập và phân tích dữ liệu thực nghiệm bao gồm: chuỗi dữ liệu cắt lớp vô tuyến vệ tinh COSMIC-1 giai đoạn 2014–2016 (mặt cắt ẩm wetPrf, Level-2 từ trung tâm CDAAC) và chuỗi dữ liệu thám không tại Trạm Khí tượng Hà Nội giai đoạn 2016–2018 tại thời điểm quan sát chuẩn $00\text{Z}$ ($07:00$ giờ địa phương GMT+7).

Literature Review và Positioning

Nghiên cứu về ảnh hưởng của khí quyển đến truyền sóng vô tuyến khởi đầu từ các công trình kinh điển của Crain (1950), Birnbaum (1950), Bean & Dutton (1968), thiết lập mối quan hệ nhiệt động lực học giữa áp suất toàn phần $P$, nhiệt độ tuyệt đối $T$, và áp suất riêng phần của hơi nước $e$ với chiết suất khí quyển $n$. Trường phái nghiên cứu thực nghiệm quốc tế sau đó phân tách thành hai hướng tiếp cận chính:

  1. Phương pháp đo trực tiếp: Sử dụng máy đo khúc xạ tần số vô tuyến (Microwave Refractometer) dựa trên nguyên lý thay đổi tần số cộng hưởng của hốc cộng hưởng $\Delta f / f \approx -\Delta n$ (Crain, 1950; Vetter, 1968; Hay, 1958). Mặc dù phương pháp này có độ chính xác cao nhất (đạt cấp độ $\pm 0,1\text{ N-unit}$), việc triển khai bị giới hạn nghiêm trọng bởi chi phí thiết bị đắt đỏ, khối lượng lớn và không thể bao phủ mạng lưới diện rộng.
  2. Phương pháp đo gián tiếp và viễn thám vệ tinh: Sử dụng số liệu thám không (Radiosonde) kết hợp mô hình hóa khí động học (Babin et al., 1997; Basha et al., 2013). Bước đột phá xuất hiện từ thập niên 2000 với kỹ thuật cắt lớp vô tuyến GNSS (GNSS Radio Occultation - RO) triển khai trên chòm vệ tinh quỹ đạo thấp LEO như COSMIC (Rocken et al., 2000; Anthes et al., 2008), cho phép trích xuất cấu trúc thẳng đứng của khí quyển toàn cầu với độ phân giải cao và tính khả dụng liên tục trong mọi điều kiện thời tiết.

Trong các tài liệu quốc tế tồn tại cuộc tranh luận kéo dài về độ tin cậy của mô hình tham chiếu toàn cầu ITU-R P.453. Các nghiên cứu tại vùng ôn đới (châu Âu, Bắc Mỹ) cho thấy mô hình hàm số mũ phản ánh tương đối sát giá trị trung bình nhiều năm của độ khúc xạ bề mặt ($N_s \approx 310 - 325\text{ N-units}$). Ngược lại, các công trình khảo sát tại khu vực nhiệt đới gió mùa như Ấn Độ (Sivaramakrishnan, 1989), Nigeria (Willoughby et al., 2002), và vùng Vịnh Ba Tư chứng minh rằng độ khúc xạ bề mặt có thể tăng vọt lên trên $390 - 410\text{ N-units}$ vào mùa mưa ẩm, gây ra hiện tượng ống dẫn sóng bề mặt (surface ducting) với tỷ lệ xuất hiện vượt quá 30% thời gian, điều mà mô hình ITU-R P.453 tiêu chuẩn không thể mô phỏng chính xác.

Luận án định vị nghiên cứu của mình tại giao điểm giữa kỹ thuật viễn thông và vật lý khí quyển nhiệt đới. Bằng việc so sánh trực tiếp với hai chuẩn đối sánh quốc tế lớn—gồm mô hình toàn cầu ITU-R P.453 và cơ sở dữ liệu cắt lớp vô tuyến COSMIC CDAAC của Hoa Kỳ/Đài Loan—luận án đã lấp đầy khoảng trống dữ liệu thực nghiệm tại miền Bắc Việt Nam, đưa ra các đường cong phân bố tham số $N, G, k, M$ đặc thù cho khu vực Hà Nội.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp vào hệ thống lý thuyết lan truyền sóng điện từ trong môi trường điện môi phân lớp thông qua việc mở rộng các mô hình tính toán chiết suất cho khí hậu nhiệt đới đặc thù. Cụ thể, các định luật và công thức nền tảng được hệ thống hóa và kiểm chứng:

  1. Độ khúc xạ vô tuyến $N$: Được định nghĩa từ chỉ số khúc xạ thực $n$ thông qua quan hệ $N = (n - 1) \times 10^6$. Giá trị $N$ được mô hình hóa qua ba thành phần nhiệt động: $$N = 77,6 \frac{P}{T} - 5,6 \frac{e}{T} + 3,75 \times 10^5 \frac{e}{T^2}$$ Trong tính toán kỹ thuật, công thức gần đúng được áp dụng với sai số dưới $0,02%$ trong dải nhiệt độ từ $-50^\circ\text{C}$ đến $+40^\circ\text{C}$: $$N \approx \frac{77,6}{T}\left(P + 4810 \frac{e}{T}\right) = 77,6\frac{P}{T} + 3,732 \times 10^5 \frac{e}{T^2}$$ Trong đó $P$ là áp suất toàn phần (hPa hoặc mbar), $T$ là nhiệt độ tuyệt đối (K), và $e$ là áp suất thành phần của hơi nước (hPa hoặc mbar), được xác định từ độ ẩm tương đối $RH$ và áp suất hơi bão hòa $e_s$: $e = \frac{RH \cdot e_s}{100}$.

  2. Hệ số bán kính Trái Đất hiệu dụng $k$ (k-factor): Được xây dựng dựa trên nguyên lý tương đương hình học giữa tia sóng uốn cong trên Trái Đất thực (bán kính $a \approx 6370\text{ km}$) và tia sóng truyền thẳng trên mô hình Trái Đất tương đương (bán kính $R_e$): $$\frac{1}{R_e} = \frac{1}{a} + \frac{dn}{dh}$$ $$k = \frac{R_e}{a} \approx \frac{1}{1 + \frac{a}{10^6}\frac{dN}{dh}} = \frac{1}{1 + 0,00637 G}$$ với $G = \frac{dN}{dh}$ là độ dốc khúc xạ (refractivity gradient, đơn vị $\text{N-units/km}$).

  3. Mô-đun khúc xạ $M$ (Modified Refractivity): Dùng để phân tích cơ chế khúc xạ ống dẫn sóng trên mô hình Trái Đất phẳng: $$M = N + 10^6 \frac{h}{a} \approx N + 157h\quad (\text{M-units})$$ $$\frac{dM}{dh} = G + 157\quad (\text{M-units/km})$$

Từ hệ thống lý thuyết trên, ba mệnh đề khoa học (Propositions) được xác lập:

  • Proposition 1: Khi $G = -40\text{ N-units/km}$, hệ số $k = 4/3$, tương ứng với điều kiện khúc xạ tiêu chuẩn (Standard Refraction), cự ly truyền sóng tầm nhìn thẳng (LOS) tăng $15%$ so với trường hợp truyền thẳng trong chân không ($k = 1$).
  • Proposition 2: Khi $-157\text{ N-units/km} < G \le -80\text{ N-units/km}$, môi trường chuyển sang trạng thái siêu khúc xạ (Super-refraction), gia tăng nguy cơ suy hao đa đường (multipath fading).
  • Proposition 3: Khi $G < -157\text{ N-units/km}$ (tương ứng $\frac{dM}{dh} < 0$), điều kiện khúc xạ ống dẫn (Ducting) hình thành, bán kính Trái Đất hiệu dụng $R_e < 0$, bẫy sóng điện từ lan truyền vượt đường chân trời (Over-The-Horizon - OTH).

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp ba trụ cột lý thuyết: Lý thuyết truyền sóng tầng đối lưu (Bean & Dutton), Lý thuyết nghịch đảo cắt lớp vô tuyến (Abel Inversion Transform cho tín hiệu GPS dải L1/L2), và Nhiệt động lực học khí quyển. Luận án thiết lập phương pháp đánh giá sai số định lượng thông qua hai chỉ số:

  • Độ lệch tuyệt đối: $\Delta N = N_h - N_s\quad (\text{N-units})$
  • Độ lệch tương đối: $\Delta \varepsilon = 100 \times \frac{\Delta N}{N_s}\quad (%)$

Trong đó $N_h$ là độ khúc xạ thực nghiệm tính từ dữ liệu quan trắc ở độ cao $h$, và $N_s$ là giá trị tham chiếu tính theo mô hình hàm số mũ ITU-R P.453: $$N_s(h) = N_0 \exp\left(-\frac{h}{h_0}\right) = 315 \exp\left(-\frac{h}{7,35}\right)$$

Điều kiện biên của khung phân tích: độ cao khảo sát từ mặt đất ($0\text{ km}$) đến $20\text{ km}$, tần số vô tuyến áp dụng từ VHF, UHF đến dải sóng vi ba ($30\text{ MHz} - 100\text{ GHz}$), độ ẩm tương đối của cảm biến thám không $RH > 15%$.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý thực chứng hậu hiện đại (Post-positivism / Critical Realism), kết hợp đa nguồn dữ liệu thực nghiệm độc lập (Dual-source empirical triangulation). Thiết kế nghiên cứu bao gồm việc thu thập song song hai tập dữ liệu độc lập về không gian và thời gian tại khu vực Hà Nội:

  • Dữ liệu viễn thám vệ tinh không gian: Dữ liệu cắt lớp vô tuyến từ chòm vệ tinh COSMIC-1 (hợp tác NASA/NOAA/Đài Loan).
  • Dữ liệu quan trắc mặt đất tại hiện trường: Dữ liệu bóng thám không vô tuyến (Radiosonde) thả từ Trạm Khí tượng Hà Nội.

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu được thiết kế chặt chẽ qua các giai đoạn chuẩn hóa:

  1. Thu thập dữ liệu cắt lớp vô tuyến COSMIC-1: Tải các tệp dữ liệu profile ẩm (wetPrf, Level-2) từ Trung tâm Lưu trữ và Phân tích Dữ liệu COSMIC (CDAAC). Dữ liệu gốc ở định dạng netCDF (Network Common Data Format) và CDL (Common Data Language). Các trường thông số trích xuất bao gồm: độ cao hình học (Geopotential Height), nhiệt độ khí quyển ($T$), áp suất toàn phần ($P$), áp suất hơi nước ($e$), và chiết suất trích xuất sẵn để đối soát.
  2. Thu thập dữ liệu thám không: Thu thập chuỗi quan trắc của bóng thám không mang máy thăm dò vô tuyến (Radiosonde mang cảm biến nhiệt độ bạch kim, cảm biến độ ẩm điện dung và áp kế màng mỏng) thả định kỳ tại thời điểm $00\text{Z}$ ($07:00$ sáng giờ địa phương). Dữ liệu được giải mã và cấu trúc hóa về định dạng CSV chuẩn.
  3. Phân tích độ nhạy sai số cảm biến khí tượng: Luận án thực hiện phân tích vi phân hàm truyền để đánh giá sự lan truyền sai số cảm biến vào giá trị độ khúc xạ $N$: $$dN = a_i dT + b_i de + c_i dP$$ Với các hệ số độ nhạy theo mô hình khí quyển chuẩn ICAO tại các độ cao khác nhau:
    • Tại mặt đất ($0\text{ km}$): $a_0 = -1,27$, $b_0 = 4,50$, $c_0 = 0,27$.
    • Tại độ cao $1\text{ km}$: $a_1 = -1,09$, $b_1 = 4,72$, $c_1 = 0,28$.
    • Tại độ cao $3\text{ km}$: $a_3 = -0,86$, $b_3 = 5,17$, $c_3 = 0,29$.
    • Tại độ cao $10\text{ km}$: $a_{10} = -0,50$, $b_{10} = 7,52$, $c_{10} = 0,30$.

Phân tích trên chứng minh rằng: ảnh hưởng của áp suất $c_i$ gần như không đổi theo độ cao; ảnh hưởng của nhiệt độ $a_i$ giảm dần theo độ cao; trong khi ảnh hưởng của áp suất hơi nước $b_i$ tăng mạnh theo độ cao (từ $4,50$ lên $7,52$). Do đó, yêu cầu về độ chính xác của cảm biến độ ẩm càng lên cao càng khắt khe. Luận án chỉ ra rằng sai số đo độ khúc xạ $\pm 1,0\text{ N-unit}$ đòi hỏi cảm biến nhiệt độ phải đạt độ chính xác trong ngưỡng $\pm 0,8^\circ\text{C}$, áp suất toàn phần trong ngưỡng $\pm 3,7\text{ mbar}$ và áp suất hơi nước trong ngưỡng $\pm 0,22\text{ mbar}$.

[Dữ liệu COSMIC-1 (wetPrf, netCDF)]       [Dữ liệu Bóng thám không 00Z (CSV)]
    [Trích xuất P, T, e, ZTD]                [Trích xuất P, T, RH, h]
   [Tính N_COSMIC theo Eq. 1.3]              [Tính e & N_Sonde theo Eq. 1.3]
                      [Tính toán Gradient G = dN/dh]
                      [Xác định k-factor & Mô-đun M]
                [Đối sánh với Mô hình chuẩn ITU-R P.453]
             [Tính độ lệch tuyệt đối ΔN & tương đối Δε%]
             [Đánh giá điều kiện truyền sóng: Khúc xạ thường,
              Siêu khúc xạ, Khúc xạ ống dẫn cho khu vực Hà Nội]

Data và phân tích

Luận án xây dựng phần mềm chuyên dụng bằng ngôn ngữ lập trình khoa học để tự động hóa quá trình đọc, lọc, nội suy và xử lý thống kê chuỗi dữ liệu lớn hàng nghìn tệp netCDF và CSV. Các bước kiểm tra tính vững (robustness checks) bao gồm: loại bỏ các điểm dị biệt (outliers) do cảm biến thám không bị trễ thời gian (time-lag constant) khi bay qua các lớp mây có độ dốc độ ẩm cao; đối soát độc lập giữa giá trị $N$ tính gián tiếp qua $(P, T, e)$ với trường chiết suất do thuật toán CDAAC cung cấp.

Phát hiện đột phá và implications

Những phát hiện then chốt

Quá trình phân tích thống kê chuỗi số liệu thực nghiệm giai đoạn 2014–2018 tại khu vực Hà Nội đã mang lại 4 nhóm phát hiện then chốt:

  1. Quy luật biến thiên theo mùa cực đoan của độ khúc xạ bề mặt ($N_s$):

    • Dữ liệu thám không và COSMIC-1 cho thấy độ khúc xạ bề mặt tại Hà Nội dao động từ $310\text{ N-units}$ (vào mùa đông khô lạnh, tháng 12 – tháng 1) lên tới $390 - 415\text{ N-units}$ (vào mùa hè nóng ẩm, tháng 6 – tháng 8). Giá trị trung bình năm của $N_s$ thực nghiệm đạt khoảng $365 - 375\text{ N-units}$, cao hơn đáng kể so với mức chuẩn $315\text{ N-units}$ của mô hình ITU-R P.453.

  2. Mức độ sai lệch lớn của mô hình toàn cầu ITU-R P.453 ở tầng thấp ($0 - 3\text{ km}$):

    • Tại các lớp sát mặt đất ($h < 1\text{ km}$), độ lệch tuyệt đối $\Delta N$ giữa dữ liệu thực nghiệm và mô hình ITU-R dao động từ $+40\text{ N-units}$ đến $+75\text{ N-units}$, tương ứng với độ lệch tương đối $\Delta \varepsilon%$ từ $12%$ đến trên $20%$.
    • Khi lên độ cao trên $5\text{ km}$, độ ẩm khí quyển giảm mạnh, giá trị $N$ thực nghiệm tiệm cận sát đường cong hàm số mũ của ITU-R, với độ lệch $\Delta N < 5\text{ N-units}$ ($\Delta \varepsilon < 2%$). Điều này chứng minh sai số của mô hình ITU-R tại vùng nhiệt đới tập trung hoàn toàn ở lớp biên khí quyển ẩm sát bề mặt.

  3. Phát hiện cấu trúc siêu khúc xạ và khúc xạ ống dẫn trong lớp khí quyển thấp:

    • Thống kê độ dốc khúc xạ $G = \frac{dN}{dh}$ trong lớp $0 - 500\text{ m}$ chỉ ra rằng giá trị trung bình $G$ tại Hà Nội không cố định ở mức chuẩn $-40\text{ N-units/km}$ ($k = 4/3$) mà thường xuyên biến động trong dải từ $-60\text{ N-units/km}$ đến $-110\text{ N-units/km}$ (khúc xạ phụ và siêu khúc xạ).
    • Đặc biệt vào các buổi sáng sớm mùa thu - đông (quan sát lúc $00\text{Z}$ tức $07:00$ GMT+7), do hiện tượng nghịch nhiệt bức xạ ban đêm kết hợp lớp ẩm bề mặt, độ dốc $G$ nhiều lần đạt giá trị âm sâu dưới $-157\text{ N-units/km}$ ($\frac{dM}{dh} < 0$), tạo ra các lớp khúc xạ ống dẫn bề mặt (surface ducts) với độ cao từ $100\text{ m}$ đến $300\text{ m}$.

  4. Độ tương thích cao giữa cắt lớp vô tuyến COSMIC-1 và bóng thám không ở độ cao trên $2\text{ km}$:

    • Phân tích đối sánh giữa dữ liệu vệ tinh COSMIC-1 và bóng thám không cho thấy hệ số tương quan rất cao ($R^2 > 0,95$) ở các tầng đối lưu tự do ($2 - 12\text{ km}$). Sự phân kỳ nhẹ ở tầng dưới $1,5\text{ km}$ được giải thích bởi hạn chế hình học bám tín hiệu của máy thu GPS LEO khi tia sóng cắt qua lớp biên hành tinh dày đặc sol khí và hơi nước.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Kết quả nghiên cứu cung cấp bộ tham số thực nghiệm thực tế về phân bố chiết suất khí quyển nhiệt đới, làm cơ sở để các nhóm công tác của ITU-R (đặc biệt là Study Group 3 - Radiowave Propagation) cập nhật bản đồ số liệu chiết suất toàn cầu cho vùng Đông Nam Á.
  • Về mặt kỹ thuật viễn thông: Cung cấp giá trị hệ số $k$-factor thực nghiệm ($k \approx 1,45 - 1,60$ trong mùa hè thay vì $k = 1,33$) cho các kỹ sư thiết kế tuyến vi ba băng rộng, giúp tính toán chính xác khoảng cự ly che chắn Fresnel, giảm thiểu fading pha sâu và nâng cao độ sẵn sàng của tuyến truyền dẫn.
  • Về định vị vệ tinh và hàng không dân dụng: Kết quả tính toán độ trễ thiên đỉnh tầng đối lưu (Zenith Tropospheric Delay - ZTD) từ dữ liệu thực nghiệm giúp hiệu chỉnh sai số định vị thẳng đứng cho máy bay khi hạ cánh tại Cảng Hàng không Quốc tế Nội Bài, đồng thời hỗ trợ hiệu chỉnh độ trễ đường truyền (Path Delay - PD) với độ chính xác cấp centimet ($1 - 3\text{ cm}$) cho các vệ tinh đo độ cao địa hình và mực nước biển như Sentinel.
  • Về an ninh quốc phòng: Việc xác định tần suất và thời điểm xuất hiện các ống dẫn sóng bề mặt tại Hà Nội giúp các lực lượng tác chiến radar phòng không và hải quân chủ động khai thác hiện tượng truyền sóng vượt đường chân trời (OTH) để mở rộng tầm trinh sát, hoặc điều chỉnh góc quét nhằm khắc phục vùng mù radar (radar holes).

Limitations và Future Research

Luận án thẳng thắn chỉ ra 3 giới hạn nghiên cứu chính:

  1. Độ phân giải tầng thấp của dữ liệu cắt lớp vô tuyến: Do hiện tượng đa đường nghiêm trọng và suy giảm biên độ tín hiệu GPS trong lớp biên khí quyển ẩm nhiệt đới, các phép đo RO COSMIC-1 ở độ cao dưới $1\text{ km}$ đôi khi bị ngắt sớm, dẫn đến mật độ dữ liệu sát mặt đất mỏng hơn so với các tầng trên.
  2. Tần suất mẫu thám không hạn chế: Dữ liệu bóng thám không tại Trạm Khí tượng Hà Nội chỉ thu thập chủ yếu 1 lần/ngày tại thời điểm $00\text{Z}$ ($07:00$ GMT+7), chưa bao quát đầy đủ chu kỳ biến động ngày - đêm (diurnal cycle), đặc biệt là thời điểm đầu giờ chiều khi nhiệt độ bề mặt đạt cực đại.
  3. Phạm vi không gian cục bộ: Nghiên cứu giới hạn trong bán kính $2^\circ$ kinh vĩ quanh Hà Nội, chưa mở rộng ra toàn bộ 7 vùng kinh tế - xã hội của Việt Nam với các điều kiện địa hình phức tạp như ven biển miền Trung hay đồi núi Tây Bắc.

Hướng nghiên cứu tiếp theo được đề xuất:

  • Khai thác chòm vệ tinh thế hệ mới COSMIC-2 (với góc nghiêng quỹ đạo thấp, chuyên sâu cho vùng nhiệt đới) kết hợp hệ thống trạm định vị vệ tinh mặt đất liên tục (VNGEONET) để giải bài toán đồng hóa dữ liệu cắt lớp vô tuyến 3D/4D thời gian thực.
  • Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) dự báo chỉ số $k$-factor và nguy cơ hình thành ống dẫn sóng vô tuyến trước $6 - 24$ giờ dựa trên số liệu dự báo thời tiết số trị (NWR).
  • Mở rộng dải tần nghiên cứu lên băng sóng milimet (EHF, $30 - 300\text{ GHz}$) và bức xạ Terahertz phục vụ mạng di động 6G.

Tác động và ảnh hưởng

Công trình nghiên cứu của NCS. Phạm Chí Công mang lại tác động sâu rộng:

  • Tác động học thuật: Đặt nền móng cho hướng nghiên cứu liên ngành giữa Kỹ thuật Vô tuyến và Vật lý Khí quyển tại Việt Nam; cung cấp nguồn dữ liệu đầu vào chuẩn hóa cho các mô hình nghiên cứu tại Viện Vật lý Địa cầu (Viện Hàn lâm Khoa học và Công nghệ Việt Nam).
  • Tác động công nghiệp và hạ tầng số: Hỗ trợ các doanh nghiệp viễn thông lớn (Viettel, VNPT, MobiFone) tối ưu hóa chi phí đầu tư hạ tầng trạm thu phát sóng thông qua việc giảm thiểu dự phòng công suất phát không cần thiết và bố trí độ cao anten vi ba tối ưu.
  • Tác động kinh tế - xã hội: Giảm thiểu chi phí đo đạc khí quyển truyền thống đắt đỏ bằng cách tận dụng nguồn dữ liệu viễn thám mở (COSMIC) và dữ liệu khí tượng sẵn có, tiết kiệm hàng tỷ đồng ngân sách khảo sát.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận xử lý dữ liệu cắt lớp vô tuyến netCDF/CDL và thuật toán ước lượng chiết suất khí quyển chính xác.
  • Kỹ sư thiết kế vô tuyến viễn thông & Radar: Sử dụng các giá trị $k$-factor, độ dốc $G$, và mô-đun $M$ thực nghiệm để cấu hình hệ thống truyền dẫn chính xác tại miền Bắc.
  • Tổng công ty Quản lý Bay Việt Nam (VATM) & Cục Hàng không: Ứng dụng mô hình trễ ZTD chính xác nhằm nâng cao an toàn dẫn đường bay trong điều kiện thời tiết phức tạp.
  • Cơ quan quản lý tần số vô tuyến (Cục Tần số Vô tuyến - Bộ TT&TT): Cơ sở khoa học để hoạch định phân bổ băng tần và kiểm soát can nhiễu đồng kênh giữa các đài phát tầm xa do hiện tượng siêu khúc xạ.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
    Luận án đã chứng minh định lượng rằng mô hình hàm số mũ tham chiếu toàn cầu ITU-R P.453 đánh giá thấp độ khúc xạ tầng thấp tại Hà Nội từ $40 - 75\text{ N-units}$ (sai số tương đối $> 15 - 20%$) trong mùa ẩm, qua đó đề xuất bộ tham số thực nghiệm địa phương hóa cho phương trình biến thiên chiết suất theo chiều thẳng đứng tại vùng nhiệt đới gió mùa Việt Nam.

  2. Tính đổi mới trong phương pháp luận so với các công trình trước đây?
    Khác với các nghiên cứu trước đây tại Việt Nam vốn chỉ dựa vào dữ liệu bề mặt đơn lẻ hoặc mô hình lý thuyết thuần túy, luận án lần đầu tiên tích hợp thành công hai phương pháp thực nghiệm độc lập: phương pháp viễn thám cắt lớp vô tuyến vệ tinh LEO (COSMIC-1) và phương pháp quan trắc bóng thám không vô tuyến (Radiosonde), thiết lập một quy trình xử lý dữ liệu tự động hóa chuẩn mực từ trích xuất định dạng khoa học netCDF/CDL đến ước lượng tham số truyền sóng.

  3. Phát hiện gây bất ngờ nhất về mặt số liệu thực nghiệm là gì?
    Đó là sự xuất hiện với tần suất đáng kể của các lớp khúc xạ ống dẫn sóng ($G < -157\text{ N-units/km}$) và siêu khúc xạ trong lớp khí quyển sát mặt đất ($0 - 300\text{ m}$) vào các thời điểm giao mùa và sáng sớm mùa khô lạnh, trái ngược với giả định truyền sóng tiêu chuẩn ($k = 4/3$) thường được áp dụng mặc định trong các phần mềm quy hoạch mạng viễn thông.

  4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp minh bạch không?
    Toàn bộ cấu trúc tệp dữ liệu, giải thuật trích xuất dữ liệu netCDF của CDAAC, quy trình chuẩn hóa dữ liệu thám không CSV, cùng hệ thống công thức giải tích xác định $N, G, k, M$ và sai số lan truyền đều được trình bày chi tiết trong luận án và các phụ lục kỹ thuật, cho phép các nhà khoa học khác tái lập hoàn toàn quy trình tính toán.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
    Lộ trình 10 năm bao gồm: (1) Mở rộng mạng lưới quan trắc chiết suất cho 7 vùng kinh tế sinh thái của Việt Nam; (2) Tích hợp dữ liệu chòm vệ tinh COSMIC-2 và mạng lưới GNSS đa hệ thống (GPS, GLONASS, Galileo, BeiDou) để xây dựng bản đồ số khúc xạ 4D động theo thời gian thực; (3) Nghiên cứu ảnh hưởng của sol khí và biến đổi khí hậu đến cấu trúc chiết suất tầng đối lưu dải sóng siêu cao tần.

Kết luận

Luận án của TS. Phạm Chí Công là một công trình nghiên cứu công phu, nghiêm túc và có giá trị khoa học cao. Sáu đóng góp cụ thể của luận án bao gồm:

  1. Hệ thống hóa toàn diện cơ sở lý luận về ảnh hưởng của tầng đối lưu nhiệt đới đến lan truyền sóng vô tuyến trong dải tần VHF, UHF và vi ba.
  2. Xây dựng thành công bộ công cụ phần mềm tự động hóa trích xuất và xử lý dữ liệu cắt lớp vô tuyến vệ tinh LEO định dạng netCDF/CDL và dữ liệu thám không.
  3. Xác lập bộ dữ liệu thực nghiệm đầu tiên về phân bố theo độ cao của độ khúc xạ $N$, độ dốc $G$, và mô-đun khúc xạ $M$ cho khu vực Hà Nội giai đoạn nhiều năm.
  4. Chứng minh định lượng sự sai lệch đáng kể của mô hình toàn cầu ITU-R P.453 tại lớp khí quyển thấp nhiệt đới và kiến nghị áp dụng số liệu thực nghiệm.
  5. Phát hiện và cảnh báo quy luật xuất hiện của hiện tượng siêu khúc xạ và khúc xạ ống dẫn sóng tại vùng đồng bằng Bắc Bộ.
  6. Đề xuất quy trình tích hợp dữ liệu viễn thám - khí tượng có hiệu quả kinh tế cao, mở ra hướng ứng dụng thiết thực cho các ngành viễn thông, hàng không, khí tượng và quốc phòng Việt Nam.