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ề mặt vật lý, nơi các đại lượng nhiệt độ ($T$), áp suất khí quyển ($P$), và áp suất thành phần hơi nước ($e$) biến thiên liên tục theo không gian và thời gian. Khi sóng điện từ truyền qua môi trường này, sự biến đổi của chiết suất môi trường (chỉ số khúc xạ vô tuyến $n$) gây ra hiện tượng bẻ cong quỹ đạo tia sóng theo định luật Snell, tạo ra trễ pha, tán xạ, hấp thụ năng lượng và fading đa đường (multipath fading). Luận án tiến sĩ kỹ thuật viễn thông của nghiên cứu sinh Phạm Chí Công (2023) với đề tài "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" đã thiết lập một hệ thống nghiên cứu thực nghiệm tiên phong tại Việt Nam, kết hợp đa nguồn số liệu viễn thám không gian và quan trắc thám không tại chỗ để mô hình hóa cấu trúc khúc xạ khí quyển tầng thấp.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc thiết kế các tuyến thông tin vô tuyến (VHF, UHF, sóng vi ba SHF/EHF), hệ thống định vị vệ tinh GNSS (GPS, GLONASS, BeiDou, Galileo), ra-đa giám sát và đo cao vệ tinh (satellite altimetry) tại Việt Nam phần lớn vẫn đang áp dụng mô hình khí quyển chuẩn toàn cầu của Liên minh Viễn thông Quốc tế (ITU-R P.453). Tuy nhiên, các mô hình chuẩn này được xây dựng dựa trên số liệu trung bình thống kê của khí hậu ôn đới, hoàn toàn không phản ánh chính xác các đặc thù khí hậu nhiệt đới gió mùa có 4 mùa rõ rệt, độ ẩm biến động mạnh và tồn tại các lớp nghịch nhiệt cục bộ như khu vực Hà Nội.

Các câu hỏi nghiên cứu và giả thuyết được xác định rõ ràng:

  • RQ1: Cấu trúc phân bố không gian và biến thiên theo mùa của độ khúc xạ vô tuyến $N$ và gradient độ dốc thẳng đứng $G = dN/dh$ tại khu vực Hà Nội diễn ra theo quy luật thực nghiệm nào?
  • RQ2: Mức độ sai lệch giữa mô hình hàm số mũ toàn cầu ITU-R P.453 ($N(h) = N_0 e^{-h/h_0}$) so với số liệu thực nghiệm đo đạc từ vệ tinh và thám không tại chỗ là bao nhiêu?
  • RQ3: Phương pháp cắt lớp vô tuyến (Radio Occultation - RO) và bóng thám không (Radiosonde) có độ tương thích và khả năng phối hợp như thế nào trong việc lập bản đồ điều kiện truyền sóng tầng đối lưu?
  • H1: Độ dốc khúc xạ $G$ và hệ số bán kính Trái Đất hiệu dụng $k$ tại tầng biên khí quyển Hà Nội (0–2 km) có độ lệch đáng kể so với giá trị khúc xạ chuẩn ($G = -40\text{ N-units/km}, k = 4/3$), làm gia tăng xác suất xuất hiện siêu khúc xạ và khúc xạ ống dẫn (ducting) trong các tháng giao mùa nóng ẩm.
  • H2: Tích hợp dữ liệu cắt lớp vô tuyến COSMIC-1 và bóng thám không cho phép hiệu chỉnh chính xác độ trễ thiên đỉnh tầng đối lưu (ZTD) và xác lập bản đồ tham số truyền sóng thực nghiệm với độ tin cậy cao hơn hẳn mô hình tham chiếu toàn cầu.

Khung lý thuyết của luận án được xây dựng vững chắc dựa trên Lý thuyết Khúc xạ Khí quyển Tầng đối lưu (Bean & Dutton, 1968), Khuyến nghị ITU-R P.453-14 / ITU-R P.834, Lý thuyết Tán xạ Đàn hồi Mie – Rayleigh và Nguyên lý Cắt lớp Vô tuyến GNSS-LEO. Về phạm vi không gian, luận án tập trung khảo sát khu vực Hà Nội với tọa độ trung tâm $21.01^\circ\text{N}$, $105.80^\circ\text{E}$ trong bán kính $2^\circ$ kinh vĩ. Dữ liệu thực nghiệm bao gồm chuỗi quan trắc đa năm (2014–2018) từ chòm vệ tinh COSMIC-1 (level-2 wetPrf) và trạm bóng thám không khí tượng Hà Nội (mã trạm 48820), mang lại ý nghĩa đột phá trong việc tối ưu hóa đường truyền viễn thông và bảo đảm độ chính xác định vị hàng không, trinh sát ra-đa quốc phòng.

Literature Review và Positioning

Lý thuyết truyền sóng vô tuyến trong khí quyển đã được hình thành và phát triển qua nhiều thập kỷ, khởi đầu từ các công trình kinh điển của Bean & Dutton (1968), Hall (1979), Freeman (2007) và liên tục được chuẩn hóa thông qua các khuyến nghị của ITU-R (P Series). Nghiên cứu chỉ số khúc xạ vô tuyến $n$ và độ khúc xạ $N$ (với $N = (n - 1) \times 10^6$) được xác định thông qua phương trình bán thực nghiệm ba thành phần:

$$N = 77{,}6 \frac{P}{T} - 5{,}6 \frac{e}{T} + 3{,}75 \times 10^5 \frac{e}{T^2}$$

hoặc dạng rút gọn với sai số dưới $0{,}02%$ trong dải nhiệt độ $-50^\circ\text{C}$ đến $+40^\circ\text{C}$ (Smith & Weintraub, 1953; ITU-R P.453-14):

$$N = 77{,}6 \frac{P}{T} + 3{,}73 \times 10^5 \frac{e}{T^2}$$

Trong y văn thế giới, tồn tại hai luồng quan điểm và tranh luận kỹ thuật chính:

  1. Tranh luận về phương pháp thu thập dữ liệu: Trường phái đo trực tiếp sử dụng khúc xạ kế vi sóng khoang cộng hưởng (Crain, 1950; Birnbaum, 1950; Vetter, 1967) có ưu thế về độ chính xác tuyệt đối do chỉ phụ thuộc vào một đại lượng tần số cộng hưởng ($\Delta f / f \approx -\Delta n$), nhưng thiết bị cồng kềnh, chi phí đắt đỏ và không thể triển khai trên diện rộng. Ngược lại, trường phái đo gián tiếp qua thông số khí tượng ($P, T, RH$) được hỗ trợ bởi mạng lưới bóng thám không (radiosonde) và vệ tinh cắt lớp vô tuyến GNSS-LEO (Ware et al., 1996; Rocken et al., 2000; Kursinski et al., 1997) cho phép bao phủ không gian rộng lớn nhưng đòi hỏi thuật toán chuyển đổi phức tạp và kiểm soát lan truyền sai số nghiêm ngặt.
  2. Tranh luận về tính phổ quát của mô hình ITU-R P.453: Mô hình hàm số mũ tham chiếu chuẩn $N(h) = N_0 e^{-h/h_0}$ (với $N_0 = 315\text{ N-units}$, $h_0 = 7{,}35\text{ km}$) được giả định áp dụng trên toàn cầu. Tuy nhiên, các nghiên cứu độc lập tại các vùng khí hậu khác nhau đã chỉ ra những sai lệch cục bộ nghiêm trọng.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Grabner & Kvicera (2003) tại Cộng hòa Séc (khí hậu ôn đới lục địa) cho thấy biến thiên độ dốc $G$ duy trì độ ổn định cao quanh giá trị $-40\text{ N-units/km}$, hiện tượng khúc xạ ống dẫn gần như chỉ xuất hiện vào các đợt nghịch nhiệt bức xạ ban đêm mùa đông ở độ cao sát mặt đất.
  • Nghiên cứu của Adediji & Ajewole (2008) tại Akure, Nigeria (khí hậu nhiệt đới xích đạo ẩm) ghi nhận độ khúc xạ bề mặt $N_s$ dao động rất lớn giữa mùa khô ($320\text{ N-units}$) và mùa mưa ($385\text{ N-units}$), với sự xuất hiện dày đặc của các lớp siêu khúc xạ ($G < -80\text{ N-units/km}$) ở tầng thấp do hàm lượng hơi nước tăng vọt.

Luận án của Phạm Chí Công định vị nghiên cứu của mình tại giao điểm của viễn thám vô tuyến hiện đại và khí tượng nhiệt đới gió mùa Đông Nam Á. Công trình đã lấp đầy khoảng trống dữ liệu thực nghiệm tại vùng đồng bằng Bắc Bộ Việt Nam—nơi có sự chuyển giao phức tạp giữa gió mùa Đông Bắc khô lạnh và gió mùa Đông Nam nóng ẩm—cung cấp bằng chứng thực nghiệm đối sánh định lượng giữa mô hình toàn cầu ITU-R và dữ liệu viễn thám thực tế.

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

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

Luận án đã mở rộng ứng dụng của lý thuyết quang hình và điện từ học truyền sóng vô tuyến trong môi trường khí quyển không đồng nhất nhiệt đới thông qua các đóng góp lý thuyết cụ thể:

  1. Kiểm chứng và xác lập giới hạn biên (boundary conditions) của mô hình ITU-R P.453: Chứng minh rằng mô hình tham chiếu lũy thừa toàn cầu chỉ có giá trị xấp xỉ ở các tầng khí quyển cao ($h > 3\text{ km}$), trong khi ở tầng biên đối lưu ($0 \le h \le 2\text{ km}$), sự phân tầng của hơi nước làm cho sai số tuyệt đối $\Delta N$ và sai số tương đối $\Delta R$ lệch chuẩn đáng kể so với mô hình lý thuyết.
  2. Phát triển mệnh đề lý thuyết về hệ số bán kính Trái Đất hiệu dụng thực nghiệm:
    • Mệnh đề 1 (P1): Tại khu vực Hà Nội, hệ số $k$-factor không phải là một hằng số $4/3 \approx 1{,}33$ như lý thuyết khúc xạ tiêu chuẩn, mà là một biến ngẫu nhiên phụ thuộc theo mùa, dao động trong dải từ $k = 1{,}15$ (khúc xạ phụ) đến $k > 2{,}0$ hoặc $k < 0$ (siêu khúc xạ và khúc xạ ống dẫn).
    • Mệnh đề 2 (P2): Độ dốc khúc xạ thẳng đứng $G = dN/dh$ chịu sự chi phối áp đảo của gradient áp suất hơi nước $de/dh$ ở độ cao dưới $1\text{ km}$, vượt trội so với ảnh hưởng của gradient nhiệt độ $dT/dh$ và áp suất $dP/dh$.
+-------------------------------------------------------------------------------+
|                       KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN                     |
+-------------------------------------------------------------------------------+
|   [Lý thuyết Khúc xạ Khí quyển]  [Nhiệt động lực học Khí quyển]  [Quang hình học LEO-RO]  |
|          (ITU-R P.453-14)            (Clausius-Clapeyron)        (Biến đổi ngược Abel)   |
+-------------------------------------------------------------------------------+
+-------------------------------------------------------------------------------+
|                           TÍCH HỢP ĐA NGUỒN DỮ LIỆU                           |
+-------------------------------------------------------------------------------+
+-------------------------------------------------------------------------------+
|                          XỬ LÝ VÀ TÍNH TOÁN THAM SỐ                           |
+-------------------------------------------------------------------------------+
+-------------------------------------------------------------------------------+
|                      ĐỐI SOÁT VÀ XÁC LẬP MÔ HÌNH KHU VỰC                      |
|       So sánh đối chuẩn với mô hình ITU-R P.453 (Delta N, Delta R)             |
+-------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên ngành ba hệ lý thuyết cốt lõi:

  • Lý thuyết nhiệt động lực học chất khí ẩm: Mô tả sự phụ thuộc của áp suất hơi bão hòa $e_s$ vào nhiệt độ qua phương trình thực nghiệm nâng cao với hệ số tăng cường $EF$ (Enhancement Factor) cho pha lỏng: $$e_s = EF \cdot a \cdot \exp\left(\frac{b \cdot t}{t + c}\right)$$ với các hằng số khí quyển chuẩn $a = 6{,}1115$, $b = 23{,}036$, $c = 279{,}82$, kết hợp độ ẩm tương đối $RH$ để trích xuất chính xác áp suất riêng phần hơi nước $e = (RH / 100) \cdot e_s$.
  • Lý thuyết cắt lớp vô tuyến vệ tinh (Radio Occultation Inversion): Khai thác nguyên lý tia sóng GNSS truyền qua khí quyển bị uốn cong một góc uốn $\alpha(a)$ do trọng lực và gradient chiết suất, từ đó giải bài toán ngược Abel để khôi phục cấu trúc thẳng đứng của chiết suất khí quyển $n(r)$.
  • Lý thuyết tương đương Trái Đất phẳng và mô-đun khúc xạ $M$: Sử dụng phép biến đổi ánh xạ quỹ đạo cong trên Trái Đất thực sang tia sóng thẳng trên Trái Đất tương đương: $$M = N + 10^6 \frac{h}{a} \approx N + 157h \quad (\text{với } h \text{ tính bằng km})$$ cho phép nhận diện điều kiện tạo ống dẫn sóng (ducting) khi $dM/dh = G + 157 < 0\text{ M-units/km}$.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị theo quan điểm bản thể học hiện thực và nhận thức luận thực chứng (positivism/empirical stance), trong đó các quy luật vật lý của môi trường truyền sóng được kiểm chứng thông qua các phép đo đạc định lượng khách quan. Luận án áp dụng thiết kế đa tầng dữ liệu (multi-level data design) kết hợp giữa quan trắc viễn thám không gian từ quỹ đạo Trái Đất tầm thấp (LEO) và quan trắc in-situ tại mặt đất.

Phạm vi không gian được xác định chặt chẽ: tọa độ tâm trạm khí tượng Hà Nội ($21^\circ01'\text{N}, 105^\circ48'\text{E}$), mở rộng trong phạm vi bán kính $2^\circ$ kinh vĩ ($19^\circ\text{N} - 23^\circ\text{N}, 103{,}8^\circ\text{E} - 107{,}8^\circ\text{E}$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và tiền xử lý dữ liệu được thiết kế đạt tiêu chuẩn quốc tế:

  1. Dữ liệu cắt lớp vô tuyến COSMIC-1:
    • Thu thập từ Trung tâm Phân tích và Lưu trữ Dữ liệu COSMIC (CDAAC) trực thuộc UCAR/NOAA/NASA.
    • Sử dụng tệp dữ liệu mặt cắt ẩm chuyên sâu wetPrf (Level-2 data) với định dạng khoa học netCDF (Network Common Data Format) và cấu trúc ngôn ngữ mô tả dữ liệu chung CDL.
    • Các biến trường trích xuất bao gồm: Độ cao hình học ($MSL_Alt$), Áp suất ($Pres$), Nhiệt độ ($Temp$), Áp suất hơi nước ($Vap_Pres$), và Độ khúc xạ trích xuất ($Ref$). Độ phân giải thẳng đứng đạt từ $0{,}1\text{ km}$ đến $0{,}2\text{ km}$.
  2. Dữ liệu bóng thám không (Radiosonde Sounding):
    • Thu thập từ Trạm Khí tượng Thám không Hà Nội thuộc Trung tâm Khí tượng Thủy văn Quốc gia.
    • Thiết bị đo đạt chuẩn WMO, bóng mang thiết bị vô tuyến thăm dò (radiosonde Vaisala) thả định kỳ 2 lần/ngày vào các thời điểm chuẩn quốc tế 00:00 UTC (07:00 giờ Hà Nội) và 12:00 UTC (19:00 giờ Hà Nội).
    • Đo trực tiếp các thông số thẳng đứng: $P$ (hPa), $T$ ($^\circ\text{C}$), $RH$ ($%$), độ cao địa thế vị $H$ (gpm) từ mặt đất lên đến đỉnh tầng đối lưu ($> 16\text{ km}$).
  3. Phân tích sai số và lan truyền độ không đảm bảo đo: Độ chính xác của phép đo gián tiếp độ khúc xạ được mô hình hóa theo phương trình vi phân sai số toàn phần: $$\Delta N = a_i \Delta T + b_i \Delta e + c_i \Delta P$$ với các trọng số nhạy cảm theo độ cao $h$ (dựa trên khí quyển chuẩn ICAO):
    • Tại mặt đất ($h = 0\text{ km}$): $a_0 = -1{,}27$, $b_0 = 4{,}50$, $c_0 = 0{,}27$.
    • Tại $h = 10\text{ km}$: $a_{10} = -0{,}50$, $b_{10} = 7{,}52$, $c_{10} = 0{,}30$. Phân tích cho thấy sai số đo lường của cảm biến hơi nước ($b_i$) tăng mạnh theo độ cao, đòi hỏi các quy trình lọc kiểm tra chất lượng (QC) nghiêm ngặt nhằm loại bỏ các cấu hình dị thường (profile outliers).

Data và phân tích

Tác giả đã lập trình các module phần mềm chuyên dụng bằng ngôn ngữ khoa học để tự động hóa việc phân tích hàng nghìn tệp dữ liệu cấu hình ẩm:

  • Kỹ thuật trích xuất và đồng bộ: Chuyển đổi nhị phân netCDF sang cấu trúc dữ liệu phẳng đa chiều (CSV, mảng cấu trúc), nội suy đồng bộ độ cao giữa các lát cắt vệ tinh và bóng thám không.
  • Xác định chỉ số thống kê sai lệch:
    • Độ lệch tuyệt đối: $\Delta N = |N_h - N_s|\text{ (N-units)}$
    • Độ lệch tương đối: $\Delta R = 100 \cdot \frac{|N_h - N_s|}{N_s}\text{ (%)}$ trong đó $N_h$ là giá trị tính toán thực nghiệm ở độ cao $h$, và $N_s$ là giá trị mô hình chuẩn ITU-R P.453 ($N_s = 315 e^{-h/7{,}35}$).
  • Tính toán gradient thẳng đứng và kiểm tra độ bền vững (robustness checks): Tính vi phân số $G = \frac{N(h_{j+1}) - N(h_j)}{h_{j+1} - h_j}$, làm mượt bằng bộ lọc trung bình trượt và ước lượng độ tin cậy thống kê tại các khoảng tin cậy $95%$.

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

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

Dữ liệu phân tích thực nghiệm giai đoạn 2014–2018 tại khu vực Hà Nội đã mang lại 4 phát hiện khoa học quan trọng:

+-------------------------------------------------------------------------------+
|            ĐẶC ĐIỂM BIẾN THIÊN ĐỘ KHÚC XẠ N VÀ GRADIENT G TẠI HÀ NỘI         |
+-------------------------------------------------------------------------------+
| Độ cao (h)   Độ khúc xạ N (N-units)   Gradient G (N-units/km)   Trạng thái khúc xạ    |
| 0 km (Mặt đất)     375 - 410              -65 đến -95           Khúc xạ phụ / Siêu    |
| 1 km               310 - 335              -45 đến -60           Khúc xạ phụ           |
| 2 km               260 - 280              -38 đến -42           Khúc xạ chuẩn         |
| > 5 km             140 - 160              -20 đến -25           Ổn định khí quyển cao |
+-------------------------------------------------------------------------------+
  1. Giá trị độ khúc xạ bề mặt cao vượt trội so với chuẩn quốc tế: Độ khúc xạ bề mặt ($N_s$) thực tế tại Hà Nội dao động phổ biến trong khoảng từ $375\text{ N-units}$ đến $410\text{ N-units}$ vào mùa hè ẩm ướt (tháng 5–tháng 9), và từ $330\text{ N-units}$ đến $360\text{ N-units}$ vào mùa khô lạnh (tháng 12–tháng 2). Giá trị này cao hơn đáng kể so với mức tham chiếu chuẩn toàn cầu $N_0 = 315\text{ N-units}$ của ITU-R.
  2. Sai lệch cực đại tập trung tại tầng biên 0–2 km: Độ lệch tuyệt đối $\Delta N$ giữa số liệu thực nghiệm và mô hình ITU-R P.453 đạt cực đại tại độ cao sát đất từ $0\text{ km}$ đến $1{,}5\text{ km}$, với $\Delta N$ dao động từ $40\text{ N-units}$ đến trên $75\text{ N-units}$, tương ứng độ lệch tương đối $\Delta R$ lên tới $15% - 25%$. Khi độ cao đạt trên $3\text{ km}$, $\Delta R$ giảm nhanh xuống dưới $5%$ và tiến dần về $0$ ở độ cao trên $6\text{ km}$.
  3. Phân bố gradient độ dốc $G$ bất thường và xuất hiện hiện tượng siêu khúc xạ: Tại lớp khí quyển sát đất ($0 - 1\text{ km}$), giá trị trung bình của $G = dN/dh$ đạt từ $-65\text{ N-units/km}$ đến $-95\text{ N-units/km}$ (khúc xạ phụ và tiệm cận siêu khúc xạ), hoàn toàn khác biệt so với giá trị chuẩn $-40\text{ N-units/km}$. Đặc biệt, vào các thời điểm sáng sớm (00Z) trong mùa hè và mùa thu, hiện tượng nghịch nhiệt kết hợp độ ẩm ngưng tụ cục bộ đã tạo ra các lớp ống dẫn sóng khí quyển với $G < -157\text{ N-units/km}$ và $dM/dh < 0$.
  4. Tính tương thích cao giữa phương pháp COSMIC-1 và bóng thám không: Đối chiếu trực tiếp giữa cấu hình $N(h)$ tính từ vệ tinh COSMIC-1 và bóng thám không cùng thời điểm cho thấy mức độ tương đồng vượt bậc với hệ số tương quan $r > 0{,}98$, độ lệch tuyệt đối trung bình giữa hai phương pháp dưới $3{,}5\text{ N-units}$ ở độ cao từ $1\text{ km}$ đến $10\text{ km}$, khẳng định tính khả thi của việc dùng dữ liệu LEO-RO thay thế hoặc bổ trợ cho các trạm thám không truyền thống.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp bộ dữ liệu thực nghiệm chuẩn hóa đầu tiên về mặt cắt khúc xạ đối lưu khu vực Hà Nội, đặt nền móng cho việc hiệu chỉnh hệ số thực nghiệm trong mô hình ITU-R cho toàn bộ vùng nhiệt đới gió mùa Đông Nam Á.
  • Phương pháp luận: Thiết lập thành công chu trình kết hợp viễn thám vô tuyến vệ tinh (LEO-RO) với bóng thám không khí tượng, chứng minh phương pháp cắt lớp vô tuyến là công cụ tin cậy để khảo sát khí quyển diện rộng với chi phí tối ưu.
  • Ứng dụng công nghệ và viễn thông:
    • Quy hoạch tuyến viba và phát thanh truyền hình (VHF/UHF): Cần sử dụng hệ số bán kính Trái Đất hiệu dụng thực tế $k \approx 1{,}45 - 1{,}60$ cho mùa hè thay vì $k = 4/3$ mặc định, giúp tính toán chính xác khoảng cự ly thông tin tầm nhìn thẳng (LOS), giảm thiểu can nhiễu đồng kênh và tối ưu hóa độ cao an-ten.
    • Hệ thống thông tin vệ tinh và định vị GNSS: Hỗ trợ mô hình hóa chính xác độ trễ thiên đỉnh tầng đối lưu (Zenith Tropospheric Delay - ZTD), loại bỏ sai số định vị từ $0{,}5\text{ m}$ đến $2{,}0\text{ m}$ do khúc xạ khí quyển ẩm gây ra.
    • Vệ tinh đo độ cao (Satellite Altimetry): Cung cấp tham số hiệu chỉnh trễ đường truyền (Path Delay - PD) đạt độ chính xác cấp centimet ($1 - 3\text{ cm}$) cho các nhiệm vụ viễn thám hải văn và giám sát mực nước biển.
    • Hệ thống ra-đa trinh sát: Dự báo chính xác hiện tượng "vùng mù ra-đa" hoặc cự ly trinh sát vượt đường chân trời (Over-The-Horizon - OTH) khi xảy ra khúc xạ ống dẫn.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn kỹ thuật:

  1. Độ phân giải không gian và thời gian của dữ liệu COSMIC-1: Do đặc điểm quỹ đạo của chòm vệ tinh COSMIC thế hệ 1, số lượng lát cắt vô tuyến đi qua phạm vi bán kính $2^\circ$ của Hà Nội trong một ngày còn hạn chế, dẫn đến việc lấy mẫu theo chuỗi thời gian chưa đạt tính liên tục từng giờ.
  2. Hạn chế đo lường của cảm biến bóng thám không ở tầng thấp và độ trễ trôi: Bóng thám không di chuyển theo hướng gió trong quá trình bay lên, khiến quỹ đạo đo không hoàn toàn là một đường thẳng đứng cục bộ. Cảm biến độ ẩm có hệ số trễ thời gian nhất định khi đi qua các lớp nghịch nhiệt mỏng.
  3. Chưa triển khai đo đạc trực tiếp đối chứng bằng khúc xạ kế vi sóng: Do điều kiện thiết bị thực nghiệm tại chỗ, nghiên cứu chủ yếu khai thác phương pháp đo gián tiếp qua các tham số khí tượng.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) được định hình qua 4 hướng trọng tâm:

  • Mở rộng khai thác dữ liệu từ chòm vệ tinh thế hệ mới COSMIC-2 kết hợp GNSS đa hệ (GPS, Galileo, BeiDou, GLONASS) nhằm tăng mật độ lát cắt khí quyển lên gấp 4–5 lần.
  • Mở rộng phạm vi nghiên cứu từ khu vực Hà Nội ra toàn bộ 7 vùng kinh tế - xã hội của Việt Nam, xây dựng Bản đồ số khúc xạ vô tuyến quốc gia theo từng mùa.
  • Tích hợp mô hình dự báo thời tiết số phân giải cao (WRF) với dữ liệu đồng hóa cắt lớp vô tuyến để thiết lập hệ thống cảnh báo sớm thời gian thực về hiện tượng khúc xạ ống dẫn phục vụ hàng không và hải quân.
  • Nghiên cứu ảnh hưởng của sol khí (aerosol) và bụi mịn ô nhiễm đô thị đến tán xạ và suy hao tín hiệu vô tuyến băng tần cực cao (mmWave / Terahertz).

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

Nghiên cứu của tác giả Phạm Chí Công tạo ra những tác động lan tỏa sâu rộng:

  • Tác động học thuật: Đóng góp 03 công trình khoa học công bố trên các tạp chí chuyên ngành uy tín, cung cấp cơ sở dữ liệu tham chiếu quan trọng cho các nhà nghiên cứu trong lĩnh vực vật lý khí quyển, vô tuyến điện tử và viễn thông tại khu vực Đông Nam Á.
  • Chuyển đổi công nghiệp viễn thông & phát thanh truyền hình: Giúp các tập đoàn viễn thông (Viettel, VNPT, MobiFone) và Đài Truyền hình tối ưu hóa thiết kế tuyến truyền dẫn vi ba số, tính toán chính xác ngân sách đường truyền (link budget), giảm biên độ dự phòng fading không cần thiết, tiết kiệm chi phí đầu tư hạ tầng trạm thu phát.
  • An ninh quốc phòng và dẫn đường hàng không: Hỗ trợ Quân chủng Phòng không - Không quân và Cục Hàng không Dân dụng (ICAO/VATM) nắm bắt chính xác quy luật truyền sóng ra-đa, dự báo cự ly phát hiện mục tiêu bay thấp trong điều kiện thời tiết phức tạp.
  • Lợi ích kinh tế - xã hội: Nâng cao hiệu quả dự báo thời tiết cực đoan, hỗ trợ giám sát biến đổi khí hậu thông qua phân tích xu thế hơi nước tổng cộng (PWV) trích xuất từ dữ liệu trễ tầng đối lưu GNSS.

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

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Viễn thông, Vật lý Địa cầu: Tiếp cận một khung phương pháp luận hoàn chỉnh về xử lý dữ liệu viễn thám netCDF, thuật toán tính toán quang hình vô tuyến và quy trình đối soát mô hình thực nghiệm.
  • Các nhà khoa học cao cấp và chuyên gia thiết kế hệ thống: Sử dụng các công thức thực nghiệm và bảng số liệu gradient độ dốc $G$ khu vực để hiệu chỉnh các mô hình truyền sóng cục bộ thay thế cho khuyến nghị ITU-R chuẩn.
  • Kỹ sư R&D trong công nghiệp viễn thông, ra-đa và định vị vệ tinh: Ứng dụng thuật toán tính toán $k$-factor và ZTD trực tiếp vào các phần mềm quy hoạch mạng vô tuyến (như Atoll, ICS Telecom) và thiết bị định vị chính xác cao (RTK, PPP).
  • Cơ quan quản lý tần số vô tuyến (Cục Tần số Vô tuyến điện): Có căn cứ khoa học thực nghiệm vững chắc để xây dựng quy chuẩn kỹ thuật quốc gia về phân bổ phổ tần, quy định mức công suất phát và khoảng cách an toàn chống can nhiễu giữa các trạm vô tuyến cố định và lưu động.

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ì?

Đóng góp độc đáo nhất là việc thiết lập và chứng minh thực nghiệm giới hạn bất cập của mô hình tham chiếu toàn cầu ITU-R P.453-14 tại tầng biên đối lưu ($0 - 2\text{ km}$) vùng nhiệt đới gió mùa Hà Nội. Luận án đã chỉ ra rằng độ khúc xạ bề mặt thực tế cao hơn mô hình lý thuyết từ $15% - 25%$, đồng thời thiết lập mối quan hệ phi tuyến giữa gradient độ dốc thẳng đứng $G$ và hàm lượng hơi nước nhiệt đới ẩm, bổ khuyết cho lý thuyết khúc xạ sóng điện từ trong môi trường khí quyển phân tầng phức tạp.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu truyền thống tại Việt Nam vốn chỉ dựa trên số liệu quan trắc trạm khí tượng bề mặt hoặc bóng thám không đơn lẻ, luận án là công trình tiên phong kết hợp đồng thời phương pháp cắt lớp vô tuyến vệ tinh quỹ đạo thấp (COSMIC-1 LEO-RO Level-2 wetPrf) với số liệu bóng thám không chuyên sâu. Quy trình này cung cấp độ phân giải thẳng đứng siêu mịn ($0{,}1 - 0{,}2\text{ km}$), cho phép trích xuất cấu trúc không gian ba chiều của chỉ số khúc xạ với độ tin cậy tương đương phương pháp đo trực tiếp nhưng với chi phí tối ưu.

3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự xuất hiện với tần suất đáng kể của các lớp nghịch đảo gradient khúc xạ ($G < -157\text{ N-units/km}$, tương ứng $dM/dh < 0$) tại khu vực Hà Nội vào các buổi sáng sớm mùa hè và thu. Hiện tượng này tạo ra các ống dẫn sóng vô tuyến tự nhiên (tropospheric ducting) sát mặt đất, làm biến dạng hoàn toàn quỹ đạo truyền sóng LOS thông thường, có khả năng bẻ cong tia sóng đi men theo bề mặt cong của Trái Đất và gây ra can nhiễu tín hiệu nghiêm trọng ngoài tầm nhìn danh định.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án mô tả chi tiết toàn bộ quy trình công nghệ tái lập:

  • Phương thức truy xuất nguồn dữ liệu mở CDAAC netCDF và dữ liệu trạm thám không 48820.
  • Mã nguồn và thuật toán xử lý dữ liệu cấu trúc CDL/CSV.
  • Bộ công thức toán học xác định các đại lượng trung gian ($EF, e_s, e, N_h, N_s, \Delta N, \Delta R, G, k, M$).
  • Tiêu chí lọc dữ liệu nhiễu và quy chuẩn hiển thị đồ thị phân bố.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình 10 năm hướng tới việc xây dựng "Hệ thống số hóa khí tượng vô tuyến thời gian thực toàn quốc" (National Real-time Radio-meteorological Digital Twin), bao gồm: tích hợp mạng lưới vệ tinh GNSS/COSMIC-2 đa chùm; thiết lập bản đồ động 3D về hệ số $k$-factor và hiện tượng ducting; phát triển các module AI/Machine Learning dự báo suy hao truyền sóng cho mạng thông tin 6G và vệ tinh viễn thám quỹ đạo siêu thấp (VLEO).

Kết luận

  1. Xác lập cơ sở khoa học thực nghiệm toàn diện: Luận án đã làm sáng tỏ bản chất truyền sóng khúc xạ trong tầng đối lưu nhiệt đới khu vực Hà Nội, phân tích định lượng sự phụ thuộc của chỉ số khúc xạ vào nhiệt độ, áp suất và đặc biệt là áp suất hơi nước biến động theo mùa.
  2. Định lượng chính xác sai lệch của mô hình toàn cầu ITU-R P.453: Chứng minh sai số tuyệt đối $\Delta N$ tại tầng sát đất đạt tới $40 - 75\text{ N-units}$ ($\Delta R \approx 15% - 25%$), đưa ra khuyến nghị mang tính bước ngoặt về việc bắt buộc phải sử dụng số liệu thực nghiệm địa phương trong thiết kế các hệ thống vô tuyến chính xác cao.
  3. Chứng minh tính ưu việt của phương pháp cắt lớp vô tuyến COSMIC-1: Dữ liệu mặt cắt ẩm vệ tinh LEO-RO đạt độ tương thích cao ($r > 0{,}98$) với bóng thám không truyền thống, khẳng định vai trò là nguồn dữ liệu viễn thám thay thế hiệu quả cho các vùng thiếu trạm quan trắc mặt đất.
  4. Phát hiện quy luật biến thiên $k$-factor và hiện tượng khúc xạ dị thường: Xác định dải biến thiên thực tế của hệ số $k$ từ $1{,}15$ đến các trạng thái siêu khúc xạ và khúc xạ ống dẫn ($G < -157\text{ N-units/km}$), cung cấp tham số cốt lõi để tối ưu hóa cự ly truyền dẫn tầm nhìn thẳng và kiểm soát can nhiễu.
  5. Mở ra 3 hướng nghiên cứu ứng dụng đột phá: (i) Hiệu chỉnh sai số tầng đối lưu ZTD/PD cho định vị GNSS và đo cao vệ tinh Sentinel; (ii) Tối ưu hóa vùng phủ sóng và chống nghẽn cho mạng thông tin di động băng rộng; (iii) Nâng cao hiệu quả giám sát ra-đa tầm xa vượt đường chân trời cho an ninh quốc phòng.
  6. Đóng góp giá trị thực tiễn và chuyển giao công nghệ: Toàn bộ kết quả và thuật toán tính toán của luận án đã sẵn sàng tích hợp vào các mô hình dự báo chuyên dụng tại Viện Vật lý Địa cầu, tạo tiền đề vững chắc cho việc xây dựng bản đồ truyền sóng vô tuyến trên phạm vi toàn lãnh thổ Việt Nam.