Tổng quan về luận án

Ô nhiễm không khí và biến đổi khí hậu đang là thách thức an ninh môi trường cấp bách tại các đô thị lớn và khu công nghiệp ở Việt Nam. Theo số liệu công bố từ Trung tâm Quan trắc Tài nguyên và Môi trường, nồng độ bụi mịn $\text{PM}_{2.5}$ trung bình 24 giờ tại Hà Nội trong nhiều thời điểm đã vượt ngưỡng giới hạn cho phép của Quy chuẩn Kỹ thuật Quốc gia QCVN 05:2013/BTNMT. Sự biến động phức tạp của các lớp sol khí (aerosol) trong tầng đối lưu và cấu trúc lớp biên khí quyển (Atmospheric Boundary Layer - ABL / Planetary Boundary Layer - PBL) đặt ra yêu cầu phải có công cụ viễn thám quang học chủ động với độ phân giải cao theo cả không gian và thời gian.

Khoảng trống nghiên cứu (research gap) then chốt xuất phát từ thực trạng: hầu hết các hệ thống Lidar (Light Detection and Ranging) hiện đại trên thế giới như MPLNET của NASA, EARLINET của Châu Âu, hay các hệ thống thương mại của Raymetrics (Hy Lạp) và Photonics Technology Obninsk (LB Nga) đều là hệ thống kín ("hộp đen"), giữ kín bí mật công nghệ về thiết kế phần cứng, thuật toán xử lý dữ liệu và có chi phí đầu tư, vận hành rất cao. Tại Việt Nam, các đề tài tiền nhiệm (như Đề tài KC.21/06-10 do Nguyễn Thế Truyện chủ nhiệm, hay luận án của Đoàn Hải Đăng, Nguyễn Xuân Tuấn) đã bước đầu tiếp cận công nghệ Lidar nhưng hoàn toàn thiếu vắng công cụ mô phỏng định lượng thông số phần cứng dựa trên phương trình Lidar trước khi chế tạo, thiếu quy trình chuẩn trực cơ quang học vi sai để tối ưu hóa tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR), và chưa xây dựng được thuật toán tự động xác định tham số giãn nở tối ưu khi trích xuất chiều cao lớp biên khí quyển.

Luận án "Nghiên cứu thiết kế, chế tạo và tích hợp hệ Lidar quan trắc sol khí khí quyển" của nghiên cứu sinh Vũ Duy Tuấn (người hướng dẫn: TS. Nguyễn Xuân Anh, TS. Nguyễn Thế Truyện) thực hiện tại Viện Nghiên cứu Điện tử, Tin học, Tự động hóa và Viện Vật lý (Viện Hàn lâm Khoa học và Công nghệ Việt Nam) đã giải quyết toàn diện bài toán này. Nghiên cứu tập trung vào 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1 & H1: Liệu việc mô phỏng đa thông số phần cứng theo phương trình Lidar trước chế tạo có cho phép tối ưu hóa việc lựa chọn linh kiện và dự báo chính xác biên độ tán xạ ngược của sol khí theo độ cao từ 0 đến 10 km? Giả thuyết khẳng định mô hình tham số hóa sẽ giảm thiểu sai số thiết kế phần cứng xuống dưới 5%.
  • RQ2 & H2: Quy trình chuẩn trực cơ quang học đa chiều có thể nâng cao hàm chồng chập $G(R)$ tại vùng mù quang học và gia tăng tỷ số SNR cho ống nhân quang điện PMT ở bước sóng 532 nm và 1064 nm không? Giả thuyết khẳng định việc hiệu chỉnh tâm quang học hội tụ chính xác sẽ mở rộng dải động đo đạc từ cự ly sát mặt đất ($< 200\text{ m}$) lên tới trên 10 km.
  • RQ3 & H3: Thuật toán biến đổi Wavelet hiệp phương sai (Wavelet Covariance Transform - WCT) cải tiến có thể tự động xác định giá trị giãn nở tối ưu $a$ để trích xuất chiều cao lớp biên khí quyển (ABL) từ số liệu Lidar thực nghiệm với độ chính xác cao hơn phương pháp gradient truyền thống hay không?

Khung lý thuyết của nghiên cứu tích hợp lý thuyết tán xạ đàn hồi Mie - Rayleigh, nguyên lý chuyển đổi quang - điện tử của ống nhân quang điện tử PMT (kết hợp mạch nguồn cao áp Hamamatsu C4900 điều khiển điện áp dynode từ $0\text{V}$ đến $-1250\text{V}$), và lý thuyết xử lý tín hiệu Wavelet Haar đa tỷ lệ. Luận án mang lại đóng góp đột phá khi làm chủ 100% công nghệ tích hợp phần cứng và phần mềm, đạt độ phân giải không gian cực cao từ $1 - 3\text{ m}$, chu kỳ lấy mẫu thời gian thực 1 phút, cung cấp chuỗi dữ liệu thực nghiệm đo đạc khí quyển Hà Nội trong các năm 2012, 2016 và 2018 với giá thành chế tạo chỉ bằng một phần nhỏ so với thiết bị nhập ngoại.

Literature Review và Positioning

Nghiên cứu viễn thám khí quyển toàn cầu đã chứng kiến sự phát triển mạnh mẽ của kỹ thuật Lidar qua nhiều thập kỷ. Phương pháp đo sol khí truyền thống bằng thiết bị thu mẫu tại chỗ (in-situ sampling) như máy đo nồng độ hạt, bóng thám không (stratostat), hoặc máy bay nghiên cứu khí quyển (U-2, WB57) tuy cung cấp chính xác thành phần hóa học và vi vật lý nhưng bị giới hạn nghiêm ngặt về độ bao phủ không gian và tính liên tục thời gian. Ngược lại, kỹ thuật viễn thám vệ tinh (như thiết bị MODIS trên vệ tinh Terra/Aqua, AVHRR trên NOAA-7) tuy bao phủ toàn cầu nhưng cảm biến thụ động phụ thuộc bức xạ mặt trời, độ phân giải thẳng đứng kém và chịu độ trễ quỹ đạo lớn.

Trong bối cảnh đó, Lidar nổi lên như công nghệ viễn thám chủ động ưu việt nhất:

  • Mạng lưới MPLNET (Micro-Pulse Lidar Network) của NASA (Welton et al., 2000; Campbell et al., 2002) sử dụng laser Nd:YLF phát xung công suất thấp bước sóng $532\text{ nm}$, tần số lặp lại cao $2500\text{ Hz}$, thu tín hiệu qua kính thiên văn đồng trục với độ phân giải thẳng đứng $75\text{ m}$. Hệ thống này từng vận hành thử nghiệm tại Trạm Vật lý Khí quyển Hà Nội từ 11/2010 đến 11/2012 nhưng đóng kín toàn bộ mã nguồn và thuật toán trích xuất biên dạng ABL.
  • Mạng lưới EARLINET (European Aerosol Research Lidar Network) tại Châu Âu (Bösenberg et al., 2003; Pappalardo et al., 2014) và mạng lưới trạm mặt đất AERONET (Holben et al., 1998) chuẩn hóa các phương pháp đo đạc quang học sol khí (độ sâu quang học sol khí AOD, hệ số suy hao $\alpha$, hệ số tán xạ ngược $\beta$).
  • Tại Châu Á, hệ thống Lidar Raman của Trung tâm Nghiên cứu Môi trường Academia Sinica (Đài Loan) do hãng Raymetrics (Hy Lạp) chế tạo sử dụng laser Nd:YAG năng lượng $200\text{ mJ}$ ở $532\text{ nm}$ và $70\text{ mJ}$ ở $355\text{ nm}$, kết hợp kính thiên văn Cassegrain đường kính $40\text{ cm}$ để đo tín hiệu tán xạ Raman của khí $\text{N}_2$ tại $387\text{ nm}$ và $607\text{ nm}$.
  • Tại Nga, hệ thống LSA-2C của Hãng Photonics Technology Obninsk ứng dụng laser Nd:YAG phát hai bước sóng $1064\text{ nm}$ ($100\text{ mJ}$) và $532\text{ nm}$ ($50\text{ mJ}$), kính Cassegrain $25\text{ cm}$, xử lý tín hiệu bằng phần mềm Luxel nhưng không công bố cấu trúc thuật toán tính chiều cao lớp biên.
  • Tại Đại học Manoa (Honolulu, Hawaii), hệ thống Lidar của Sharma et al. sử dụng laser Nd:YAG ba bước sóng ($1064\text{ nm}$, $532\text{ nm}$, $355\text{ nm}$), tần số $20\text{ Hz}$, phân kỳ chùm tia $1 - 1.5\text{ mrad}$, kính thu đường kính $0.4\text{ m}$, thu nhận tín hiệu qua dao động ký số Tektronix.
+-----------------------------------------------------------------------------+
|                          HỆ THỐNG QUAN TRẮC SOL KHÍ                         |
+--------------------------------------+--------------------------------------+
                                       |
       +-------------------------------+-------------------------------+
       |                                                               |
+------v------------------------------+ +------------------------------v------+
|       VIỄN THÁM THỤ ĐỘNG / TẠI CHỖ   | |       VIỄN THÁM QUANG HỌC CHỦ ĐỘNG   |
| - MODIS, AVHRR (Vệ tinh)            | |                 (LIDAR)              |
| - In-situ, Bóng thám không, Tháp đo | +------------------------------+------+
|   => Hạn chế: Kém liên tục, độ phân  |                                |
|      giải không gian thẳng đứng thấp|                                |
+-------------------------------------+                                |
                                       +-------------------------------+
                                       |
       +-------------------------------+-------------------------------+
       |                                                               |
+------v------------------------------+ +------------------------------v------+
|    HỆ THỐNG QUỐC TẾ / THƯƠNG MẠI    | |   HỆ THỐNG NỘI ĐỊA TIỀN NHIỆM VÀ     |
| - MPLNET (NASA): Đóng gói NetCDF     | |       ĐỀ TÀI KC.21/06-10 / NCS      |
| - Raymetrics (Hy Lạp), LSA-2C (Nga) | - Thiếu mô phỏng tham số phần cứng   |
|   => Hạn chế: Chi phí cao, công nghệ | - Chưa chuẩn hóa quy trình chuẩn trực|
|      hộp đen, thuật toán bí mật      | - Cố định tham số giãn nở Haar WCT   |
+-------------------------------------+ +------------------------------+------+
                                                                       |
                                        +------------------------------v------+
                                        |  ĐÓNG GÓP CỦA LUẬN ÁN VŨ DUY TUẤN:  |
                                        |  1. Mô phỏng phần cứng Lidar Eq.    |
                                        |  2. Chuẩn trực vi sai quang học     |
                                        |  3. Tự động hóa tham số WCT Haar     |
                                        |  4. Độ phân giải cao 1-3m, làm chủ  |
                                        +-------------------------------------+

Về mặt lý thuyết trích xuất cấu trúc ABL, tồn tại cuộc tranh luận học thuật lớn giữa ba trường phái: (1) Phương pháp Gradient cực tiểu (Gradient Method); (2) Phương pháp Phân tích cụm (Cluster Analysis); và (3) Phương pháp biến đổi Wavelet Haar (Haar Wavelet Covariance Transform - WCT). Nghiên cứu của Vanessa Caicedo et al. (2017) và Brooks (2003) chỉ ra rằng WCT vượt trội về khả năng chống nhiễu SNR thấp. Tuy nhiên, nhược điểm chí mạng của Caicedo et al. là chọn giá trị giãn nở $a$ ngẫu nhiên hoặc cố định (ví dụ $a = 30\text{ m}$), dẫn đến sai lệch lớn khi cấu trúc sol khí phân tầng phức tạp. Luận án này định vị chính xác vào khoảng trống đó: phát triển thuật toán WCT tự thích nghi tìm giá trị giãn nở tối ưu $a_{\text{opt}}$, tích hợp trên hệ Lidar tự chế tạo hoàn chỉnh tại Việt Nam.

Đó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 và làm sâu sắc thêm các lý thuyết vật lý viễn thám và toán ứng dụng trong môi trường khí quyển nhiệt đới gió mùa:

  1. Lý thuyết tán xạ đàn hồi Mie - Rayleigh tổng quát: Luận án làm rõ ranh giới chuyển tiếp giữa tán xạ phân tử khí Rayleigh (kích thước hạt $r \ll \lambda$, tỉ lệ nghịch với $\lambda^4$) và tán xạ hạt sol khí Mie ($r \approx 10^{-8} - 10^{-2}\text{ m}$) trong điều kiện độ ẩm cao tại Hà Nội. Luận án tích hợp mô hình khí quyển chuẩn MSIS-E-90 để mô phỏng tín hiệu tán xạ Rayleigh theo độ cao, tạo đường chuẩn đối sánh chính xác nhằm tách biệt hệ số suy hao sol khí $\alpha_{\text{aer}}(R)$ và hệ số tán xạ ngược $\beta_{\text{aer}}(R)$.
  2. Lý thuyết hàm chồng chập trường nhìn quang học $G(R)$: Mở rộng mô hình hình học Lidar đơn trục (mono-static) và song trục (bi-static), chứng minh toán học sự phụ thuộc phi tuyến của biên dạng tín hiệu vào góc phân kỳ laser $\theta_L$, trường nhìn kính thiên văn $\text{FOV}_T$, cự ly tách trục $d_0$ và tiêu cự hệ thấu kính.
  3. Mở rộng lý thuyết biến đổi Wavelet hiệp phương sai (WCT): Chuyển dịch từ mô hình gán tham số tĩnh sang mô hình hàm liên tục đa độ phân giải, chứng minh rằng điểm cực đại của hệ số $W_f(a, b)$ tương ứng chính xác với vị trí đảo nhiệt và ranh giới bước nhảy sol khí tại đỉnh lớp biên khí quyển.

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

Khung phân tích của luận án tích hợp chặt chẽ 3 nền tảng: Mô hình lan truyền photon theo Phương trình viễn thám Lidar $\rightarrow$ Mô hình hàm truyền quang - điện tử phần cứng $\rightarrow$ Mô hình xử lý tín hiệu ngẫu nhiên Wavelet WCT.

$$\text{Phương trình Lidar tổng quát: } P(R) = P_0 \cdot \frac{c\tau}{2} \cdot \left[ \beta_{\text{mol}}(R) + \beta_{\text{aer}}(R) \right] \cdot \frac{A_0}{R^2} \cdot \eta \cdot G(R) \cdot \exp\left{ -2 \int_0^R \left[ \alpha_{\text{mol}}(r) + \alpha_{\text{aer}}(r) \right] dr \right}$$

Trong đó:

  • $P(R)$: Công suất quang tán xạ ngược nhận được từ cự ly $R$ (W).
  • $P_0$: Công suất đỉnh xung laser phát ($\text{J/s}$); $\tau$: Độ rộng xung laser ($5\text{ ns}$).
  • $A_0$: Diện tích hiệu dụng của kính thiên văn thu tín hiệu ($\text{m}^2$).
  • $\eta$: Hiệu suất quang học toàn phần của hệ thống.
  • $G(R)$: Hàm chồng chập hình học quang học ($0 \le G(R) \le 1$).
  • $\alpha(R), \beta(R)$: Hệ số suy giảm quang học và hệ số tán xạ ngược của phân tử khí ($\text{mol}$) và sol khí ($\text{aer}$).
+------------------------------------------------------------------------------------+
|                       KHUNG PHÂN TÍCH VẬT LÝ VÀ DỮ LIỆU TOÀN PHẦN                  |
+------------------------------------------------------------------------------------+
  [KHỐI NGUỒN PHÁT LASER]           [MÔI TRƯỜNG KHÍ QUYỂN]           [KHỐI QUANG HỌC THU]
  - Nd:YAG (532nm / 1064nm)        - Phân tử khí (Rayleigh)         - Kính Celestron C11
  - Năng lượng: 50-200 mJ/xung     - Sol khí / Khói bụi (Mie)       - Bộ tách chùm PBSH
  - Tần số xung: 20 Hz             - Hấp thụ, suy giảm quang học    - Kính lọc phổ hẹp 532nm
          |                                  |                                  |
          +--------------------------------->|<---------------------------------+
                                             |
                              [TÍN HIỆU TÁN XẠ NGƯỢC P(R)]
                                             |
                                             v
                             [BỘ ĐẦU THU QUANG ĐIỆN TỬ PMT]
                             - Hamamatsu PMT + Cao áp C4900 (0 đến -1250V)
                             - Chuyển đổi quang sang dòng electron (Hiệu suất QE)
                                             |
                                             v
                             [BỘ BIẾN ĐỔI TƯƠNG TỰ - SỐ ADC]
                             - Lấy mẫu tốc độ cao, độ phân giải 1-3m
                                             |
                                             v
                            [THUẬT TOÁN XỬ LÝ DỮ LIỆU SỐ WCT]
                             - Lọc nhiễu nền bức xạ mặt trời
                             - Chuẩn hóa tín hiệu nhân bình phương khoảng cách: P(R).R^2
                             - Biến đổi Wavelet Haar hiệp phương sai tự thích nghi
                                             |
                                             v
                             [KẾT XUẤT CÁC THÔNG SỐ VẬT LÝ]
                             - Cấu trúc & Chiều cao lớp biên khí quyển (ABL / PBL)
                             - Phân bố mật độ tầng mây, sol khí và hệ số suy giảm alpha
+------------------------------------------------------------------------------------+

Điều kiện biên (boundary conditions) của khung phân tích được xác lập nghiêm ngặt: dải độ cao khảo sát tầng đối lưu từ $0.1\text{ km}$ đến $12\text{ km}$; dải nhiệt độ môi trường từ $15^\circ\text{C}$ đến $40^\circ\text{C}$; giới hạn vận hành trong điều kiện không mưa lớn hoặc sương mù đậm đặc làm tắt nghẽn hoàn toàn chùm tia laser.

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

Thiết kế nghiên cứu

Nghiên cứu tuân theo thế giới quan thực chứng (positivism) và chủ nghĩa hiện thực phản biện (critical realism), kết hợp chặt chẽ giữa tính toán giải tích, mô phỏng số trên máy tính, kỹ thuật chế tạo cơ điện tử chính xác và thực nghiệm quan trắc thực địa. Thiết kế đa tầng bao gồm:

  1. Tầng mô phỏng: Lập trình giải phương trình Lidar trên môi trường tính toán số, nhập các tham số vật lý của laser LOTIS LS-2137U, Litron Nano S120-20, Brilliant và kính thiên văn Schmidt-Cassegrain Celestron C11 (đường kính $280\text{ mm}$, tiêu cự $2800\text{ mm}$, tỷ số khẩu độ $f/10$).
  2. Tầng cơ - quang học: Thiết kế bộ gá vi chỉnh hai thấu kính chuẩn trực (thấu kính LUD 50.9-UV-355-532 và thấu kính lưỡng lồi DCX UV), khối phân chia chùm tia phân cực PBSH-450-1300-100.
  3. Tầng điện tử & đo ghép: Thiết kế mạch nguyên lý, mạch in cung cấp điện áp phân áp dynode cho PMT, ghép nối bộ biến đổi Analog-to-Digital Converter (ADC) tốc độ cao và phát triển phần mềm giao tiếp thời gian thực.
+-----------------------------------------------------------------------------+
|              QUY TRÌNH CHUẨN TRỰC CƠ QUANG HỌC VÀ HIỆU CHỈNH HỆ THỐNG       |
+-----------------------------------------------------------------------------+
   Bước 1: Cố định trục kính thiên văn thu tín hiệu (Celestron C11)
      |
   Bước 2: Phát laser dẫn đường công suất thấp định vị trục phát
      |
   Bước 3: Điều chỉnh cụm gá vi sai 2 thấu kính chuẩn trực (LUD & DCX UV)
      |
   Bước 4: Tinh chỉnh hệ ốc vi chỉnh 3 tọa độ (X, Y, Theta)
      |
   Bước 5: Đưa vết chùm laser vào chính tâm (trung tâm) trường nhìn FOV kính thu
      |
   Bước 6: Giám sát biên độ xung tán xạ Rayleigh tầng cao trên dao động ký
      |
   Bước 7: Cố định vị trí khi tỷ số tín hiệu trên nhiễu (SNR) đạt cực đại toàn dải
+-----------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình hiệu chuẩn và thực nghiệm được chuẩn hóa qua 5 bước nghiêm ngặt:

  • Chuẩn trực trục quang học: Để triệt tiêu góc lệch trục giữa chùm phát laser và trục thu của kính thiên văn, nghiên cứu thiết kế cụm cơ khí gá lắp với núm tinh chỉnh bước dịch micromet. Quy trình đưa chùm laser đồng trục/song song tuyệt đối vào đúng tâm kính thu, đảm bảo đạt tỷ số tín hiệu trên nhiễu $\text{SNR}$ tối đa và cực tiểu hóa vùng mù chồng chập hình học.
  • Triệt tiêu nhiễu nền quang học: Tích hợp kính lọc giao thoa quang phổ dải hẹp ($\text{FWHM} = 1.0\text{ nm}$) tại bước sóng trung tâm $532\text{ nm}$ ngay trước cửa sổ tiếp nhận photon của PMT, loại bỏ bức xạ tán xạ mặt trời ban ngày.
  • Tam giác đạc dữ liệu (Data Triangulation): Đối chiếu dữ liệu Lidar thu được đồng thời với trạm đo Lidar vi xung MPLNET của NASA đặt tại Hà Nội, trạm quan trắc mặt đất của Trung tâm Mạng lưới Khí tượng Thủy văn Quốc gia và trạm đo sol khí tự động quang kế mặt trời (Sunphotometer) thuộc mạng lưới AERONET.

Data và phân tích

Thuật toán biến đổi Wavelet Haar hiệp phương sai (WCT) được xây dựng dựa trên định nghĩa hàm mẹ Haar:

$$h\left(\frac{z - b}{a}\right) = \begin{cases} +1, & b - \frac{a}{2} \le z \le b \ -1, & b \le z \le b + \frac{a}{2} \ 0, & \text{khác} \end{cases}$$

Hàm biến đổi hiệp phương sai $W_f(a, b)$ đối với tín hiệu Lidar hiệu chỉnh khoảng cách $S(z) = P(z) \cdot z^2$ được xác định theo tích phân:

$$W_f(a, b) = \frac{1}{a} \int_{z_{\text{min}}}^{z_{\text{max}}} S(z) \cdot h\left(\frac{z - b}{a}\right) dz$$

Trong đó:

  • $z$: Độ cao quan sát (m).
  • $b$: Vị trí tịnh tiến của hàm Wavelet (tương ứng với độ cao khảo sát).
  • $a$: Hệ số giãn nở không gian (spatial dilation scale).
Dạng xung hàm Haar Wavelet:
h((z-b)/a)
   +1 ^          +---------------+
      |          |               |
      |          |               |
    0 +----------+---------------+---------------+--------> Độ cao z
      |                          |               |
      |                          |               |
   -1 v                          +---------------+
              b - a/2            b            b + a/2

Thuật toán của luận án quét tự động tham số $a$ trong dải biến thiên từ $10\text{ m}$ đến $300\text{ m}$ với bước nhảy $1\text{ m}$. Vị trí độ cao $z = b$ tại đó $W_f(a, b)$ đạt giá trị cực đại toàn cục (local maximum) chính là chiều cao đỉnh lớp biên khí quyển $z_{\text{ABL}}$:

$$z_{\text{ABL}} = \arg \max_{b} \left[ W_f(a_{\text{opt}}, b) \right]$$

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

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

Dữ liệu thực nghiệm thu thập qua các chiến dịch đo đạc quy mô lớn tại Hà Nội (các ngày 01/10/2012, 26/10/2016, 23/11/2016, 24/11/2016, 09/3/2018, 09/4/2018) đã mang lại 4 phát hiện khoa học mang tính bước ngoặt:

  1. Khám phá quy luật biến động nhật triều của lớp biên khí quyển Hà Nội: Dữ liệu Lidar liên tục từ trưa đến tối (chiến dịch đo 23/11/2016 và 24/11/2016 từ 13h16 đến 18h37) cho thấy chiều cao lớp biên khí quyển xáo trộn (Convective Boundary Layer) đạt đỉnh cực đại từ $1200\text{ m} - 1500\text{ m}$ vào khoảng 14h00 - 15h00 do đối lưu nhiệt mặt đất cực mạnh. Sau 16h30, khi bức xạ mặt trời suy giảm, lớp biên sụp đổ nhanh chóng xuống mức $300\text{ m} - 500\text{ m}$, hình thành lớp biên ổn định ban đêm (Stable Nocturnal Boundary Layer) bẫy toàn bộ bụi mịn $\text{PM}_{2.5}$ và sol khí sát mặt đất, lý giải hiện tượng nồng độ ô nhiễm tăng vọt về đêm tại Hà Nội.
  2. Xác lập giá trị giãn nở tối ưu $a_{\text{opt}}$ trong thuật toán WCT: Trái với các nghiên cứu quốc tế trước đây cố định $a = 30\text{ m}$, dữ liệu luận án chứng minh rằng trong điều kiện khí quyển nhiệt đới nhiều sol khí xáo trộn, giá trị giãn nở tối ưu $a$ phải nằm trong dải $120\text{ m} - 180\text{ m}$ để triệt tiêu hoàn toàn các cực đại giả do lớp sol khí cục bộ gây ra, nâng độ chính xác nhận dạng ranh giới ABL lên trên 95% ($p < 0.001$).
  3. Phát hiện cấu trúc phân tầng sol khí đa lớp (Multi-layer Aerosol Structure): Dữ liệu đo ngày 09/3/2018 và 09/4/2018 phát hiện sự tồn tại của lớp sol khí lơ lửng tồn dư (Residual Layer) tại độ cao $2.2\text{ km} - 3.5\text{ km}$, kết hợp với các tầng mây Ti (Cirrus clouds) mỏng ở độ cao $8\text{ km} - 10.5\text{ km}$ với độ phân giải không gian $1.5\text{ m}$.
  4. Kiểm chứng tính tương thích cao với chuẩn quốc tế: Tín hiệu Lidar thu được sau chuẩn trực quang học có tỷ số $\text{SNR}$ vượt trội, biên dạng tán xạ Rayleigh tại tầng cao ($> 6\text{ km}$) khớp hoàn toàn với mô hình phân tử khí chuẩn MSIS-E-90 và tương thích 92% với số liệu đo song song của trạm Lidar MPLNET (NASA) tại Hà Nội.
+-----------------------------------------------------------------------------+
|        DIỄN BIẾN ĐỘ CAO LỚP BIÊN KHÍ QUYỂN (ABL) THEO THỜI GIAN THỰC        |
+-----------------------------------------------------------------------------+
 Độ cao (m)
   ^
1500|                  * * * (Đỉnh đối lưu 14h-15h)
    |                *       *
1000|              *           *
    |            *               *
 500|  * * * * *                   * * * * * * * (Lớp biên ổn định đêm 18h-24h)
    |  (Ban sáng)
   0+------------------------------------------------------------------------->
      06:00          12:00       15:00       18:00       21:00       Thời gian
+-----------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bức tranh thực nghiệm hoàn chỉnh đầu tiên về động lực học lớp biên khí quyển và phân bố thẳng đứng của sol khí tại vùng đồng bằng Bắc Bộ, bổ sung dữ liệu thực tế cho các mô hình khí hậu khu vực Đông Nam Á.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn hóa từ mô phỏng tham số phần cứng theo phương trình Lidar, thiết kế gá lắp cơ quang học chính xác, đến thuật toán tự động xử lý tín hiệu WCT có thể chuyển giao áp dụng cho các hệ thống viễn thám quang học khác.
  • Về mặt ứng dụng thực tiễn: Hệ Lidar tự phát triển giúp các cơ quan quản lý nhà nước chủ động giám sát chất lượng không khí 24/7, phát hiện sớm các nguồn phát thải xuyên biên giới, cháy rừng hoặc đốt rơm rạ mùa vụ.
  • Về chính sách môi trường: Đóng góp cơ sở khoa học định lượng để Bộ Tài nguyên và Môi trường, Sở Tài nguyên Môi trường Hà Nội xây dựng kịch bản điều hành giao thông và sản xuất công nghiệp linh hoạt theo diễn biến nâng hạ của trần lớp biên khí quyể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 chỉ ra 4 giới hạn nghiên cứu:

  1. Giới hạn kênh đo bước sóng: Hệ thống mới tích hợp hoàn chỉnh và xử lý chuyên sâu kênh tán xạ đàn hồi Mie-Rayleigh tại bước sóng $532\text{ nm}$ và $1064\text{ nm}$, kênh tán xạ Raman đo hơi nước ($\text{H}_2\text{O}$) và khí Nitơ ($\text{N}_2$) ở $387\text{ nm}/607\text{ nm}$ chưa được vận hành đồng thời liên tục trong mọi chiến dịch.
  2. Vùng mù quang học tầm cực gần: Do cấu hình quang học hai trục song song, hàm chồng chập $G(R)$ chỉ đạt giá trị 1 tại cự ly $> 200\text{ m}$, khiến việc quan trắc nồng độ sol khí sát mặt đất ($0 - 150\text{ m}$) vẫn cần phối hợp với các trạm quan trắc in-situ.
  3. Ảnh hưởng của điều kiện thời tiết cực đoan: Tín hiệu bị suy giảm mạnh khi gặp mây dày tầng thấp (mây tích dông Cumulus) hoặc mưa rào lớn do hiện tượng đa tán xạ và suy giảm năng lượng quang học toàn phần.
  4. Chuỗi thời gian quan trắc liên tục dài hạn: Các chiến dịch đo đạc tập trung vào các đợt điển hình trong năm 2012, 2016 và 2018, chưa hình thành cơ sở dữ liệu chuỗi thời gian liên tục 10 năm để đánh giá xu thế biến đổi khí hậu dài hạn.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Tích hợp thêm kênh phân cực vi sai (Depolarization Channel) để phân biệt chính xác dạng hình học hạt sol khí (dạng cầu của sương mù/hạt ngưng kết vs dạng bất đối xứng của bụi cát sa mạc, tro bụi núi lửa).
  • Ứng dụng trí tuệ nhân tạo và mạng nơ-ron tích chập (Deep CNN) để tự động nhận dạng tầng mây, lớp biên khí quyển và phân loại sol khí từ dữ liệu ảnh phổ Lidar thời gian thực.
  • Thiết kế hệ Lidar quét 3D (Scanning Lidar) cơ động đặt trên xe chuyên dụng hoặc tàu biển để vẽ bản đồ phân bố sol khí toàn không gian đô thị.
  • Mở rộng mạng lưới trạm Lidar liên kết quốc gia kết nối với mạng lưới Lidar Châu Á (AD-Net) và MPLNET toàn cầu.

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

Luận án khẳng định vị thế tiên phong với những tác động đa tầng rõ nét:

  • Tác động học thuật: Đặt nền móng cho kỹ nghệ thiết kế - chế tạo thiết bị viễn thám quang điện tử phức tạp tại Việt Nam; tạo ra nguồn dữ liệu Lidar độ phân giải cao phục vụ công bố trên các tạp chí quốc tế chuyên ngành khí quyển (ISI/Scopus).
  • Làm chủ công nghệ và tự chủ kinh tế: Chế tạo thành công hệ Lidar nội địa với chi phí chỉ bằng $20% - 30%$ giá thành thiết bị nhập ngoại, giúp tiết kiệm hàng chục tỷ đồng ngân sách nhà nước, đồng thời hoàn toàn làm chủ công tác bảo trì, nâng cấp phần cứng và tùy biến phần mềm.
  • An ninh quốc phòng và ứng dụng kép (Dual-use): Công nghệ Lidar phát triển trong luận án có khả năng ứng dụng trực tiếp trong quốc phòng để phát hiện khí tài bay tầm thấp không người lái (UAV), đo đạc trần mây và gió phục vụ an toàn cất hạ cánh quân sự, phát hiện tác nhân sinh - hóa học trong khí quyển.
  • Lợi ích xã hội: Nâng cao độ chính xác của các bản tin dự báo ô nhiễm không khí và cảnh báo sớm thiên tai khí tượng thủy văn cho hàng triệu người dân Thủ đô Hà Nội và các vùng lân cận.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận hệ phương pháp luận thiết kế Lidar hoàn chỉnh, phương trình viễn thám tham số hóa và mã nguồn thuật toán xử lý dữ liệu Wavelet WCT chuẩn mực.
  • Các nhà khoa học vật lý quang học & khí quyển: Sở hữu chuỗi dữ liệu thực nghiệm chuẩn xác về sol khí nhiệt đới và lớp biên khí quyển để kiểm chứng các mô hình hóa khí hậu và lan truyền chất ô nhiễm.
  • Kỹ sư R&D thiết bị quang điện tử: Tham khảo bản vẽ thiết kế cụm cơ quang vi sai, sơ đồ mạch cấp nguồn cao áp điều khiển PMT và các giải pháp lọc nhiễu quang học dải hẹp.
  • Cơ quan hoạch định chính sách & Quản lý môi trường: Ứng dụng công cụ viễn thám chủ động để giám sát tuân thủ môi trường, quy hoạch không gian đô thị và đưa ra các quyết sách giảm thiểu phát thải dựa trên bằng chứng khoa học thực chứ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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mô hình hóa tường minh sự tương tác giữa phương trình truyền xạ Lidar và hàm biến đổi Wavelet Haar hiệp phương sai (WCT) đa tỷ lệ, giải quyết triệt để bài toán xác định tham số giãn nở tối ưu $a_{\text{opt}}$. Nghiên cứu đã mở rộng lý thuyết viễn thám tán xạ Mie - Rayleigh trong điều kiện khí quyển nhiệt đới gió mùa có độ ẩm cao và mật độ hạt sol khí phức tạp, chuẩn hóa việc đối sánh tín hiệu thực nghiệm với mô hình khí quyển chuẩn MSIS-E-90.

2. Đột phá về phương pháp luận của luận án so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Vanessa Caicedo et al. (2017) trên thiết bị Vaisala CL31 (chọn $a = 30\text{ m}$ cố định) và hệ thống LSA-2C của Nga (chỉ trích xuất các hệ số quang học cơ bản), luận án đã:

  • Tự động hóa thuật toán quét tìm $a_{\text{opt}}$ trong dải $10 - 300\text{ m}$, loại bỏ hoàn toàn việc gán tham số chủ quan.
  • Tích hợp thành công phần mềm mô phỏng tham số phần cứng theo phương trình Lidar trước chế tạo, kiểm soát sai số thiết kế linh kiện quang - điện tử.
  • Đạt độ phân giải không gian vượt trội ($1 - 3\text{ m}$ so với $75\text{ m}$ của MPLNET NASA).

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

Phát hiện bất ngờ nhất là hiện tượng sụp đổ cực nhanh của lớp biên xáo trộn đối lưu tại Hà Nội vào các buổi chiều mùa đông (sau 16h30): lớp biên hạ độ cao đột ngột từ $1500\text{ m}$ xuống dưới $400\text{ m}$ chỉ trong vòng chưa đầy 90 phút. Dữ liệu Lidar chứng minh rõ nét sự tồn tại của lớp biên ổn định ban đêm có độ dày rất mỏng bẫy toàn bộ bụi mịn sát mặt đất, tạo ra đỉnh ô nhiễm $\text{PM}_{2.5}$ nghịch lý vào ban đêm dù mật độ giao thông đã hạ nhiệt.

4. Luận án có cung cấp quy trình tái lặp thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết:

  • Thông số kỹ thuật các khối phần cứng: Laser LOTIS/Litron Nano, kính thiên văn Celestron C11 ($D = 280\text{ mm}$, $f = 2800\text{ mm}$), bộ cấp nguồn Hamamatsu C4900 ($15\text{V}$, ra $0$ đến $-1250\text{V}$, dòng $0.6\text{ mA}$).
  • Bản vẽ cơ khí 2D/3D cụm giá đỡ hai thấu kính chuẩn trực LUD 50.9-UV-355-532 và DCX UV.
  • Sơ đồ mạch điện nguyên lý cấp nguồn và điều khiển hệ số khuếch đại PMT.
  • Lưu đồ thuật toán đo đạc, lưu trữ và mã giả xử lý WCT trích xuất chiều cao ABL.

5. Chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?

Chương trình 10 năm định hướng: (1) Phát triển hệ Lidar phân cực đa bước sóng (High-Spectral-Resolution Lidar - HSRL); (2) Xây dựng mạng lưới 5 trạm quan trắc Lidar tự động trên toàn lãnh thổ Việt Nam; (3) Tích hợp mô hình đồng hóa dữ liệu viễn thám Lidar vào phần mềm dự báo thời tiết khí tượng thủy văn WRF-Chem; (4) Chế tạo hệ Lidar vi xung gọn nhẹ gắn trên phương tiện bay không người lái (UAV) và vệ tinh nano.

Kết luận

Luận án "Nghiên cứu thiết kế, chế tạo và tích hợp hệ Lidar quan trắc sol khí khí quyển" của tác giả Vũ Duy Tuấn là một công trình khoa học mẫu mực, thể hiện sự kết hợp hoàn hảo giữa nghiên cứu cơ bản, kỹ thuật công nghệ cao và thực nghiệm thực tế tại hiện trường. Tóm lược 6 đóng góp cốt lõi của công trình:

  1. Làm chủ quy trình thiết kế và mô phỏng phần cứng Lidar: Ứng dụng thành công phương trình Lidar để tính toán tối ưu hóa thông số linh kiện quang - điện tử trước chế tạo, giảm thiểu rủi ro và chi phí phần cứng.
  2. Đột phá trong kỹ thuật chuẩn trực cơ quang học vi sai: Thiết kế cụm gá lắp và quy trình hiệu chỉnh trục quang học chính xác cao, tối ưu hóa hàm chồng chập $G(R)$, tối đa hóa tỷ số tín hiệu trên nhiễu SNR cho đầu thu PMT.
  3. Phát triển phần mềm điều khiển và đo đạc thời gian thực: Xây dựng hệ thống thu thập dữ liệu tự động với độ phân giải không gian cực cao ($1 - 3\text{ m}$), chu kỳ quét 1 phút, hoàn toàn tự chủ công nghệ.
  4. Cải tiến thuật toán Wavelet Haar WCT tự thích nghi: Giải quyết dứt điểm bài toán xác định tham số giãn nở tối ưu $a_{\text{opt}}$, nâng cao độ chính xác trích xuất chiều cao lớp biên khí quyển lên trên 95%.
  5. Cung cấp bức tranh thực nghiệm toàn diện về sol khí đô thị: Khám phá quy luật động lực học lớp biên và sự phân tầng sol khí tại Hà Nội qua các chuỗi số liệu thực địa tin cậy năm 2012, 2016 và 2018.
  6. Mở ra kỷ nguyên tự chủ công nghệ viễn thám quang học nội địa: Khẳng định năng lực khoa học công nghệ của Việt Nam trong việc chế tạo các thiết bị quan trắc môi trường và khí quyển phức tạp, sẵn sàng chuyển giao phục vụ an ninh quốc phòng và phát triển bền vững kinh tế - xã hội.