Tổng quan về luận án

Nghiên cứu về vật lý khí quyển và quan trắc môi trường không khí đang đối mặt với thách thức lớn trước tình trạng suy thoái chất lượng không khí đô thị và biến đổi khí hậu toàn cầu. Theo số liệu quan trắc môi trường không khí tại Việt Nam, nồng độ bụi mịn PM2.5 trung bình 24 giờ tại nhiều thời điểm đã vượt ngưỡng giới hạn cho phép của QCVN 05:2013/BTNMT. Khí quyển Trái Đất có tổng khối lượng ước tính khoảng $5 \times 10^{18}\text{ kg}$, trong đó ba phần tư khối lượng tập trung ở độ cao dưới $11\text{ km}$. Dù ba khí chính gồm $\text{N}_2$ ($78%$), $\text{O}_2$ ($21%$) và $\text{Ar}$ ($0{,}93%$) chiếm tới $99{,}96%$ thể tích, phần còn lại $0{,}04%$ chứa các phần tử vi lượng và sol khí (aerosol) với kích thước từ vài nanomet đến hàng chục micromet lại đóng vai trò quyết định đối với cân bằng bức xạ, vi vật lý mây và sức khỏe con người.

Mặc dù kỹ thuật viễn thám quang học chủ động Lidar (Light Detection and Ranging) đã khẳng định ưu thế vượt trội so với các phương pháp lấy mẫu trực tiếp tại trạm mặt đất, bóng thám không hay ảnh vệ tinh thụ động nhờ khả năng đo cắt dọc liên tục với độ phân giải cao, việc làm chủ quy trình tính toán thiết kế phần cứng tối ưu và giải thuật xử lý tín hiệu tại Việt Nam vẫn tồn tại khoảng trống nghiên cứu sâu sắc (research gap). Phần lớn các hệ Lidar thương mại nhập ngoại (như LSA-2C của LB Nga hay Raymetrics của Hy Lạp) đều là "hộp đen" về thông số cấu hình và thuật toán xử lý dữ liệu lõi, đồng thời chi phí bảo dưỡng rất cao. Ngược lại, các công trình trong nước trước đây (như đề tài KC.21/06-10 hay các luận án của Đặng Hải Triều, Nguyễn Xuân Tuấn) chưa xây dựng được công cụ mô phỏng tham số phần cứng tích hợp dựa trên phương trình Lidar trước khi chế tạo, thiếu quy trình chuẩn hóa căn chỉnh quang sai - quang trục để tối ưu hàm chồng chập ($G(R)$) và tỷ số tín hiệu trên nhiễu (SNR), cũng như chưa có giải thuật tự động xác định giá trị giãn nở tối ưu trong phép biến đổi Wavelet để định vị chính xác độ cao lớp biên khí quyển (Atmospheric Boundary Layer - ABL / LBKQ).

Luận án đặt ra 3 câu hỏi nghiên cứu và 3 giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa và kiểm chứng định lượng các tham số kỹ thuật phần cứng (nguồn phát laser, kính thiên văn thu, bộ tách chùm, bộ lọc quang và cảm biến PMT) dựa trên phương trình vi phân Lidar nhằm tối ưu hóa thiết kế hệ thống trước khi tích hợp vật lý?
    Giả thuyết 1 (H1): Việc thiết lập mô hình tính toán mô phỏng toàn diện tham số phần cứng theo phương trình Lidar cho phép dự báo chính xác năng lượng tán xạ ngược thu nhận, giảm thiểu sai số thiết kế và rút ngắn thời gian tích hợp thực nghiệm.
  • Câu hỏi nghiên cứu 2 (RQ2): Quy trình căn chỉnh quang học đồng trục thực nghiệm nào cho phép hội tụ chính xác chùm tia laser phát vào tâm thị trường kính thu để tối đa hóa tỷ số SNR và nâng cao độ nhạy ghi nhận tín hiệu yếu ở cự ly gần và tầng cao?
    Giả thuyết 2 (H2): Phương pháp tinh chỉnh cơ - quang đa trục có kiểm soát tín hiệu hồi tiếp thực thời sẽ tối ưu hóa hàm chồng chập quang học ($OF$), loại bỏ hiện tượng tán xạ ký sinh và gia tăng đáng kể SNR của khối thu PMT.
  • Câu hỏi nghiên cứu 3 (RQ3): Giải thuật nào cho phép trích xuất tự động và chính xác độ cao lớp biên khí quyển từ dữ liệu tán xạ Lidar trong điều kiện khí hậu nhiệt đới gió mùa có độ ẩm và mật độ mây phức tạp?
    Giả thuyết 3 (H3): Phép biến đổi Wavelet hiệp phương sai (Wavelet Covariance Transform - WCT) sử dụng hàm mẹ Haar kết hợp thuật toán tối ưu hóa hệ số giãn nở $a$ sẽ xác định chính xác ranh giới lớp biên khí quyển với độ phân giải không gian $1\text{ - }3\text{ m}$.

Nghiên cứu được triển khai dựa trên khung lý thuyết tán xạ Mie - Rayleigh, nguyên lý nhân quang điện tử dynode, mô hình truyền xạ khí quyển MSIS-E-90 và lý thuyết sóng nhỏ (Wavelet analysis). Đóng góp đột phá của luận án là chế tạo thành công hệ Lidar quang trắc sol khí tối ưu "Made in Vietnam" với độ phân giải không gian cực cao ($1\text{ - }3\text{ m}$ so với mức tiêu chuẩn $75\text{ m}$ của mạng lưới MPLNET - NASA), đồng thời giải mã toán học quy luật biến thiên của hệ số giãn nở Wavelet để xác định ranh giới ABL. Phạm vi nghiên cứu bao gồm việc phân tích chuỗi dữ liệu thực nghiệm thu thập liên tục tại khu vực Hà Nội qua các năm 2012, 2016, 2018 và đối sánh chéo với các trạm quan trắc quốc tế.

Literature Review và Positioning

Nghiên cứu quan trắc sol khí và cấu trúc lớp biên khí quyển bằng kỹ thuật quang học trải qua ba luồng học thuật chính:

  1. Luồng nghiên cứu thiết kế cấu hình phần cứng Lidar: Khởi đầu từ các cấu hình Bi-static (không đồng trục) thập niên 1980 - 1990 với nguồn sáng liên tục của Parameswaran et al. (1995) tại Viện Vật lý Ahmedabad (Ấn Độ), sau đó chuyển dịch mạnh mẽ sang cấu hình Mono-static (đồng trục hoặc hai trục song song hẹp) phát xung Nd:YAG công suất cao của Sharma et al. (1999) tại Đại học Hawaii và các sáng chế độc quyền tại Mỹ của Lee et al. (2002, 2005) hay David S. Hall (2006).
  2. Luồng nghiên cứu mạng lưới quan trắc viễn thám chuẩn hóa: Đại diện tiêu biểu là mạng lưới Lidar xung vi mô MPLNET của NASA (Welton et al., 2001, 2007) vận hành ở bước sóng $532\text{ nm}$, mạng lưới Lidar sol khí Châu Âu EARLINET (Bösenberg et al., 2003; Pappalardo et al., 2014) và mạng lưới máy đo quang phổ kế mặt trời AERONET (Holben et al., 1998).
  3. Luồng nghiên cứu giải thuật toán học trích xuất ranh giới ABL: Bao gồm phương pháp gradient đạo hàm bậc nhất (Flamant et al., 1997), phân tích cụm k-means (Cluster analysis) và phương pháp biến đổi Wavelet hiệp phương sai Haar (Brooks, 2003; Davis et al., 2000; Baars et al., 2008; Caicedo et al., 2017).

Trong y văn tồn tại những tranh luận lý thuyết và thực nghiệm gay gắt. Tranh luận thứ nhất liên quan đến việc lựa chọn cấu hình quang học: cấu hình Bi-static tách biệt bộ phát - bộ thu giúp loại bỏ hiện tượng quá tải quang học ở cự ly $0\text{ m}$ nhưng lại chịu giới hạn nghiêm ngặt về tầm đo xa và sai số góc trường nhìn ($FOV$); ngược lại, cấu hình Mono-static đồng trục của NASA MPLNET tối ưu hóa không gian lắp đặt và tầm quan trắc nhưng lại đặt ra bài toán hóc búa về hàm chồng chập $G(R)$ vùng cận trường và hiện tượng bão hòa cảm biến. Tranh luận thứ hai liên quan đến thuật toán xác định độ cao LBKQ: nhóm nghiên cứu của Caicedo et al. (2017) khi sử dụng số liệu máy đo trần mây Vaisala CL31 đã chỉ ra rằng phương pháp Haar-Wavelet vượt trội hơn phương pháp Gradient và Phân tích cụm, tuy nhiên tác giả lại chọn giá trị hệ số giãn nở $a$ một cách ngẫu nhiên cố định (ví dụ $a = 30\text{ m}$), dẫn đến việc ranh giới ABL bị nhảy bậc giả tạo khi mật độ sol khí phân tầng phức tạp.

So sánh với hai nghiên cứu quốc tế điển hình:

  • Hệ thống Lidar MPLNET của NASA tại Hà Nội (2011 - 2012): Sử dụng laser Nd:YLF ($532\text{ nm}$, xung $6\text{ - }8\text{ }\mu\text{J}$, tần số $2500\text{ Hz}$), xuất file NetCDF với độ phân giải thẳng đứng $75\text{ m}$ và thời gian $1\text{ phút}$. Toàn bộ tham số modul phần cứng và thuật toán lõi được bảo mật tuyệt đối.
  • Nghiên cứu của Caicedo et al. (2017) tại Mỹ: Khảo sát so sánh 3 phương pháp Gradient, Cluster và Wavelet trên thiết bị Vaisala CL31 nhưng thất bại trong việc thiết lập hàm liên tục xác định hệ số giãn nở tối ưu $a$.

Luận án này đã vượt lên các giới hạn trên bằng cách làm chủ hoàn toàn chu trình từ tính toán thiết kế mô phỏng, chế tạo module cơ - quang điện tử, thiết lập phần mềm điều khiển thu thập thời gian thực, đến phát triển giải thuật WCT tối ưu hóa tự động hệ số $a$, đạt độ phân giải không gian $1\text{ - }3\text{ m}$.

Đó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à kiểm chứng thực nghiệm sâu sắc hai lý thuyết nền tảng trong vật lý quang học khí quyển:

  1. Lý thuyết tán xạ đàn hồi Mie - Rayleigh: Luận án làm rõ ranh giới chuyển tiếp giữa tán xạ phân tử khí Rayleigh (áp dụng cho các hạt có đường kính $d \ll \lambda$, nhỏ hơn bước sóng kích thích khoảng 1000 lần) và tán xạ Mie (áp dụng cho các hạt sol khí có dải kích thước $d \approx 10^{-8}\text{ - }10^{-2}\text{ m}$). Luận án mô hình hóa thành công tương tác của chùm xung laser phát tại các bước sóng $532\text{ nm}$ và $1064\text{ nm}$ với môi trường sol khí đa thành phần (bụi khoáng sa mạc, tro đốt nông nghiệp, sol khí lục địa và sol khí biển).
  2. Lý thuyết biến đổi sóng nhỏ Haar (Haar Wavelet Mother Function): Luận án mở rộng lý thuyết WCT bằng việc chứng minh toán học rằng độ sắc nét của bước nhảy suy giảm tín hiệu Lidar tại ranh giới LBKQ là một hàm phi tuyến phụ thuộc vào tỷ số tán xạ ngược của sol khí tầng đối lưu so với tầng tự do. Luận án bác bỏ giả định chọn hệ số $a$ tĩnh ngẫu nhiên, thay vào đó đề xuất tiêu chuẩn cực đại hóa hệ số tương quan biến đổi Wavelet để xác định $a$ tối ưu cho từng profile khí quyển cụ thể.

Mô hình lý thuyết tổng quát được mô tả qua phương trình vi phân Lidar đàn hồi:

$$P(R) = P_0 \cdot \frac{c\tau}{2} \cdot A \cdot \eta \cdot \frac{G(R)}{R^2} \cdot \beta(R) \cdot \exp\left[-2 \int_0^R \alpha(r)dr\right] + P_{\text{bg}}$$

Trong đó:

  • $P(R)$ là công suất quang tán xạ ngược thu nhận tại khoảng cách $R$.
  • $P_0$ là công suất xung laser phát xạ ban đầu; $\tau$ là độ rộng xung; $c$ là vận tốc ánh sáng.
  • $A$ là diện tích hữu dụng của gương sơ cấp kính thiên văn thu; $\eta$ là hiệu suất truyền quang của toàn hệ.
  • $G(R)$ là hàm chồng chập quang học hình học ($0 \le G(R) \le 1$).
  • $\beta(R) = \beta_{\text{mol}}(R) + \beta_{\text{aer}}(R)$ là hệ số tán xạ ngược toàn phần (phân tử khí + sol khí).
  • $\alpha(r) = \alpha_{\text{mol}}(r) + \alpha_{\text{aer}}(r)$ là hệ số suy giảm quang học (extinction coefficient).
  • $P_{\text{bg}}$ là công suất bức xạ phông nền môi trường.

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 3 trụ cột lý thuyết: (1) Quang điện tử vi mô (cơ chế khuếch đại điện tử phát xạ thứ cấp qua chuỗi dynode PMT theo mô hình phân bố Poisson); (2) Quang học sóng và hình học (lý thuyết truyền chùm Gaussian, hàm chồng chập $OF$ và góc mở trường nhìn $FOV$); (3) Xử lý tín hiệu số nâng cao (biến đổi Wavelet liên tục 1D và thuật toán lọc thích nghi phông nền).

Khung tiếp cận này xác lập rõ các điều kiện biên (boundary conditions):

  • Vùng cận trường ($R < 150\text{ m}$): Tín hiệu bị chi phối bởi hàm $G(R) < 1$ và độ bão hòa khuếch đại dynode của PMT.
  • Vùng viễn trường ($R > 10\text{ km}$): Tín hiệu tán xạ suy giảm theo quy luật $1/R^2$, giới hạn bởi ngưỡng dòng tối (dark current) và tỷ số $SNR \to 1$.

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

Thiết kế nghiên cứu

Nghiên cứu quán triệt lập trường triết học thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản biện (critical realism). Thiết kế nghiên cứu tuân thủ quy trình thực nghiệm kỹ thuật đa tầng (multi-level engineering design) gồm 3 giai đoạn hữu cơ:

[Giai đoạn 1: Mô phỏng số lý thuyết]
[Giai đoạn 2: Chế tạo cơ - quang - điện & Căn chỉnh hệ thống]
[Giai đoạn 3: Phát triển phần mềm & Đo đạc kiểm chứng thực địa]

Quy trình nghiên cứu rigorous

Quy trình thu thập và tiền xử lý dữ liệu được chuẩn hóa nghiêm ngặt:

  • Cấu hình phần cứng thu - phát: Nguồn phát laser Nd:YAG xung (tùy chọn dòng LOTIS LS-2137U, Litron Nano S120-20 hoặc Brilliant) phát tia đơn sắc bước sóng $532\text{ nm}$ và $1064\text{ nm}$. Anten thu quang học là kính thiên văn phản xạ Cassegrain model Celestron C11 đường kính gương sơ cấp $280\text{ mm}$ ($11\text{ inches}$), tiêu cự $2800\text{ mm}$ ($f/10$).
  • Khối biến đổi quang điện: Tín hiệu sau khi qua bộ tách chùm phân cực (PBSH-450-1300-100) và kính lọc giao thoa băng hẹp ($355\text{ nm}$, $532\text{ nm}$, $\Delta\lambda = 1\text{ - }3\text{ nm}$) được hội tụ vào cửa sổ quang của ống nhân quang điện PMT (đối với kênh $532\text{ nm}$) và đi-ốt quang tuyết lở APD (đối với kênh $1064\text{ nm}$). Khối PMT sử dụng bộ nguồn cao áp C4900 của Hamamatsu với các tham số chuẩn xác: điện áp nguồn cung cấp $+15\text{V}$, dải điện áp cao áp điều khiển liên tục từ $0\text{V}$ đến $-1250\text{V}$, dòng tiêu thụ định mức $0{,}6\text{ mA}$.
  • Quy trình căn chỉnh quang học độc quyền: Thiết lập hệ gá hai thấu kính LUD 50.9-UV-355-532 và DCX UV cùng cơ cấu vi chỉnh 4 vít góc, dịch chuyển trục phát laser song song và trùng khít với trục quang kính thu, loại bỏ méo dạng trường nhìn biểu kiến ($AFV$).
  • Đảm bảo độ tin cậy và giá trị đo lường: Tín hiệu Lidar thô $S(R)$ được khử dòng tối phông nền $P_{\text{bg}}$, nhân chuẩn hóa với bình phương khoảng cách $R^2$ để tạo thành tín hiệu hiệu chỉnh khoảng cách (Range Corrected Signal - RCS):

$$S_{\text{RCS}}(R) = [S(R) - P_{\text{bg}}] \cdot R^2$$

Data và phân tích

Phép biến đổi Wavelet hiệp phương sai Haar (Haar-WCT) được định nghĩa qua tích phân hàm bước nhảy mẹ $h\left(\frac{z - b}{a}\right)$:

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

Trong đó:

  • $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{vị trí khác} \end{cases}$
  • $a$ là hệ số giãn nở (dilation parameter), biểu thị độ dày của lớp chuyển tiếp sol khí.
  • $b$ là vị trí dịch chuyển tịnh tiến theo độ cao thẳng đứng.

Độ cao lớp biên khí quyển ($z_{\text{ABL}}$) chính là tọa độ $b$ tại đó hàm hiệp phương sai $W_f(a, b)$ đạt giá trị cực đại địa phương:

$$z_{\text{ABL}} = \arg\max_b \left[ W_f(a, b) \right]$$

Thuật toán của luận án quét tự động dải giá trị $a$ từ $30\text{ m}$ đến $300\text{ m}$ với bước nhảy $\Delta a = 5\text{ m}$, xác định giá trị $a_{\text{opt}}$ sao cho tỷ số giữa đỉnh cực đại $W_f$ và phương sai nhiễu nền đạt giá trị lớn nhất. Các phép kiểm tra độ vững (robustness checks) được thực hiện bằng cách so sánh profile tính toán từ hệ tự chế với dữ liệu của trạm đo viễn thám MPLNET - NASA tại Hà Nội (tháng 10/2012, 11/2016, 03/2018 và 04/2018), ghi nhận hệ số tương quan Pearson $r > 0{,}91$ ($p < 0{,}001$).

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

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

  1. Làm chủ tuyệt đối công nghệ mô phỏng phần cứng Lidar trước chế tạo: Kết quả chạy mô phỏng phương trình Lidar trên mô hình khí quyển chuẩn MSIS-E-90 cho phép xác định chính xác điện áp tối ưu cung cấp cho bộ cao áp C4900 (từ $-800\text{V}$ đến $-1050\text{V}$) tùy theo cường độ bức xạ ngày/đêm, duy trì hiệu suất lượng tử ($QE$) của PMT ở mức tối đa mà không gây bão hòa dòng anode.
  2. Nâng cao đột phá tỷ số SNR nhờ quy trình căn chỉnh đồng trục: Sau khi thực hiện quy trình căn chỉnh cơ - quang 4 bước với bộ gá hai thấu kính, biên độ tín hiệu tán xạ ngược tại tâm kính thu tăng $38%$, vùng mù quang học cận trường ($OF = 0$) rút ngắn từ $280\text{ m}$ xuống dưới $120\text{ m}$, cho phép thu nhận tín hiệu sol khí tầng thấp với SNR $> 15\text{ dB}$ ở cự ly $3\text{ km}$.
  3. Phát hiện quy luật phi tuyến của hệ số giãn nở Haar-Wavelet ($a$): Luận án chứng minh rằng giá trị giãn nở tối ưu $a$ không cố định mà biến thiên mạnh từ $60\text{ m}$ (đối với lớp biên ổn định ban đêm có độ dốc sol khí lớn) lên $180\text{ - }240\text{ m}$ (đối với lớp biên hỗn hợp ban ngày do đối lưu nhiệt mạnh). Việc cố định $a = 30\text{ m}$ như các nghiên cứu quốc tế trước đây dẫn đến sai số xác định đỉnh ABL lên tới $150\text{ - }250\text{ m}$.
  4. Định lượng chu kỳ nhật động của lớp biên khí quyển tại Hà Nội: Dữ liệu đo đạc thực tế liên tục trong các đợt quan sát (23/11/2016, 24/11/2016, 09/3/2018 và 09/4/2018) chỉ ra rằng độ cao LBKQ tại Hà Nội đạt mức thấp nhất khoảng $350\text{ - }500\text{ m}$ vào lúc $05\text{ - }07$ giờ sáng, sau đó phát triển mạnh do phát xạ mặt trời và đạt đỉnh đối lưu cực đại từ $1200\text{ - }1650\text{ m}$ vào khoảng $14\text{ - }16$ giờ chiều, trước khi suy sụp nhanh thành lớp cặn (residual layer) sau $18$ giờ.
  5. Ghi nhận hiện tượng bẫy sol khí trong các đợt nghịch nhiệt mùa đông: Dữ liệu Lidar phân giải cao $1\text{ - }3\text{ m}$ bóc tách rõ nét cấu trúc phân tầng sol khí phức tạp; trong các ngày nghịch nhiệt tháng 11/2016, sol khí PM2.5 bị giam giữ hoàn toàn dưới độ cao $600\text{ m}$, giải thích nguyên nhân gây ô nhiễm không khí nghiêm trọng tại Hà Nội.

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Cung cấp giải thuật mở toàn diện cho cộng đồng nghiên cứu quang học khí quyển, cho phép tích hợp linh hoạt mô hình Haar-WCT vào các trạm Lidar đo gió, Lidar đo ôzôn và trạm trần mây tự động.
  • Về mặt thực tiễn kinh tế - xã hội: Làm chủ công nghệ chế tạo hệ Lidar hoàn chỉnh trong nước với tổng kinh phí chỉ bằng $1/3$ đến $1/4$ giá thành thiết bị nhập khẩu từ Châu Âu hoặc Bắc Mỹ; hoàn toàn tự chủ về công tác bảo dưỡng, thay thế linh kiện và nâng cấp tính năng.
  • Về chính sách và quản lý nhà nước: Cung cấp công cụ giám sát phát thải theo phương thẳng đứng với độ chính xác cao cho Bộ Tài nguyên và Môi trường, Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, hỗ trợ đắc lực cho công tác cảnh báo sớm ô nhiễm bụi mịn đô thị và lập quy hoạch phát triển bền vững.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 hạn chế kỹ thuật:

  1. Giới hạn quan trắc trong điều kiện thời tiết cực đoan: Khi trời có sương mù dày đặc hoặc mưa rào cường độ lớn, xung laser $532\text{ nm}$ và $1064\text{ nm}$ bị suy giảm quang học cực mạnh ($Extinction \to \infty$), khiến chùm tia không thể xuyên thấu qua các lớp mây thấp dày để thu thập dữ liệu tầng trên $2\text{ km}$.
  2. Hiện tượng nhiễu bức xạ mặt trời ban ngày: Dù đã trang bị kính lọc giao thoa băng hẹp ($\Delta\lambda = 1\text{ nm}$), công suất phông nền $P_{\text{bg}}$ vào buổi trưa nắng gắt vẫn làm suy giảm tỷ số SNR khoảng $25\text{ - }30%$ so với đo đạc ban đêm, đòi hỏi phải giảm điện áp cấp cho PMT.
  3. Độ bao phủ không gian hạn chế: Hệ Lidar chế tạo hoạt động ở chế độ phát thẳng đứng cố định (Zenith pointing), chưa tích hợp đầu quét 3D (3D Scanning scanner) để quét góc phương vị và góc tà nhằm dựng bản đồ ô nhiễm không gian toàn cảnh thành phố.
  4. Chuỗi số liệu chưa bao phủ trọn vẹn 4 mùa: Các đợt đo đạc tập trung chủ yếu vào các tháng mùa thu - đông và chuyển mùa (tháng 3, 4, 10, 11), chưa có số liệu quan sát liên tục 365 ngày để đánh giá đầy đủ quy luật khí hậu học sol khí dài hạn.

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

  • Bổ sung kênh tán xạ Raman từ khí $\text{N}_2$ ($387\text{ nm}$) và hơi nước ($407\text{ nm}$) để đo trực tiếp hệ số suy giảm $\alpha(R)$ độc lập mà không cần giả định tỷ số Lidar ($S_a$).
  • Nâng cấp hệ gá quét chùm tia 3 chiều (Azimuth-Zenith Scanner) phục vụ giám sát nguồn thải công nghiệp cục bộ.
  • Tích hợp mạng nơ-ron học sâu (Deep Learning / CNN-LSTM) kết hợp Haar-WCT để tự động phân loại nguồn gốc hạt sol khí (khói sinh khối, bụi đường, bụi cát) theo thời gian thực.
  • Mở rộng mạng lưới liên trạm Lidar quan trắc tầng biên tại các đô thị trọng điểm của Việt Nam (Hà Nội, TP. Hồ Chí Minh, Đà Nẵng).

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

  • Tác động học thuật: Tạo bước nhảy vọt cho chuyên ngành Vật lý Kỹ thuật và Quang điện tử tại Việt Nam. Khung tính toán mô phỏng phần cứng và giải thuật Wavelet tối ưu của luận án là tài liệu tham khảo chuẩn mực cho các nhóm nghiên cứu viễn thám, dự báo tiềm năng trích dẫn cao trên các tạp chí quốc tế chuyên ngành (như Applied Optics, Atmospheric Environment, Remote Sensing).
  • Tác động công nghiệp và chuyển giao công nghệ: Đặt nền móng cho các viện nghiên cứu (Viện Nghiên cứu Điện tử, Tin học, Tự động hóa - VIELINA; Viện Vật lý - VAST) và doanh nghiệp công nghệ trong nước làm chủ dây chuyền sản xuất thiết bị viễn thám Lidar công nghệ cao.
  • Lợi ích xã hội và môi trường: Giúp cơ quan quản lý nhà nước định lượng chính xác độ dày tầng xáo trộn khí quyển, phục vụ dự báo chất lượng không khí (AQI) hàng ngày, góp phần bảo vệ sức khỏe cộng đồng trước hiểm họa bụi mịn PM2.5 và PM10.

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

  • Nghiên cứu sinh và Giảng viên ngành Quang tử / Vật lý Khí quyển: Tiếp cận toàn bộ phương pháp luận mô phỏng phương trình Lidar, thiết kế mạch ghép nối PMT C4900 và thuật toán WCT hoàn chỉnh không bị che giấu mã nguồn.
  • Kỹ sư R&D Thiết bị Quang điện tử: Thụ hưởng quy trình công nghệ căn chỉnh đồng trục hai thấu kính, kỹ thuật lọc nhiễu quang học và giải pháp chế tạo cơ khí chính xác cho kính thiên văn Cassegrain.
  • Chuyên gia Khí tượng Thủy văn và Môi trường: Sở hữu bộ dữ liệu cắt dọc sol khí và độ cao lớp biên khí quyển có độ phân giải siêu cao ($1\text{ - }3\text{ m}$) phục vụ kiểm chứng các mô hình dự báo thời tiết số trị (WRF, WRF-Chem).
  • Nhà hoạch định chính sách đô thị: Có cơ sở dữ liệu khoa học thực chứng để ban hành các quy chế kiểm soát khí thải giao thông và công nghiệp trong các kịch bản thời tiết nghịch nhiệt nguy hại.

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 lý thuyết độc đáo nhất là việc mở rộng Lý thuyết biến đổi Wavelet hiệp phương sai (Haar-WCT) trong việc nhận dạng cấu trúc lớp biên khí quyển phân tầng. Thay vì áp dụng hệ số giãn nở tĩnh $a$ theo kinh nghiệm như các nghiên cứu trước đây, luận án thiết lập mối quan hệ hàm số giữa hệ số giãn nở tối ưu $a_{\text{opt}}$ với gradient mật độ quang học sol khí tại vùng đệm ranh giới ($Entrainment\text{ }Zone$). Đóng góp này giải quyết triệt để sự mâu thuẫn giữa độ mịn phân giải không gian và độ ổn định kháng nhiễu trong phép biến đổi sóng nhỏ.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất hai nghiên cứu trước đây?

So với nghiên cứu của Caicedo et al. (2017) (chỉ sử dụng máy đo trần mây Vaisala CL31 với hệ số $a$ cố định $30\text{ m}$) và đề tài KC.21/06-10 của Nguyễn Thế Hiếu (2010) (chế tạo Lidar nhưng thiếu mô phỏng tham số trước khi gia công và thiếu thuật toán trích xuất ABL), luận án đã:

  • Xây dựng phần mềm mô phỏng tham số quang điện tử toàn diện theo phương trình vi phân Lidar trước chế tạo.
  • Thiết lập quy trình căn chỉnh đồng trục 4 bước đạt độ phân giải không gian $1\text{ - }3\text{ m}$ (vượt bậc so với độ phân giải $75\text{ m}$ của NASA MPLNET và $24\text{ m}$ của Đại học Quốc lập Trung ương Đài Loan).

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng dữ liệu kèm theo?

Phát hiện bất ngờ nhất là sự biến thiên đột ngột của hệ số giãn nở tối ưu $a_{\text{opt}}$ trong chu kỳ chuyển tiếp ngày - đêm tại Hà Nội: vào ban đêm khi lớp biên co cụm ổn định ($350\text{ - }500\text{ m}$), $a_{\text{opt}}$ đạt giá trị rất nhỏ ($60\text{ m}$); nhưng vào ban ngày khi đối lưu nhiệt phát triển mạnh đưa ABL lên $1650\text{ m}$, lớp đệm xáo trộn mở rộng khiến $a_{\text{opt}}$ bắt buộc phải tự động nới rộng lên $240\text{ m}$ thì giải thuật mới bắt đúng đỉnh cực đại WCT. Nếu giữ nguyên $a = 30\text{ m}$ như các công bố quốc tế, thuật toán sẽ nhận diện nhầm các cụm sol khí cục bộ tầng thấp thành độ cao đỉnh ABL.

4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ giao thức tái lập thực nghiệm:

  • Bản vẽ thiết kế cơ khí và sơ đồ mạch ghép nối bộ nguồn cao áp PMT Hamamatsu C4900 (điện áp vào $+15\text{V}$, điện áp ra $0\text{V}$ đến $-1250\text{V}$, dòng $0{,}6\text{ mA}$).
  • Thông số quang học chuẩn: Kính Celestron C11 ($D = 280\text{ mm}$, $F = 2800\text{ mm}$), lăng kính tách chùm PBSH-450-1300-100, thấu kính chuẩn trực LUD 50.9-UV-355-532.
  • Lưu đồ thuật toán và mã nguồn chương trình đo đạc, thu thập và xử lý tín hiệu Haar-WCT tối ưu hóa hệ số $a$.

5. Kế hoạch nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?

Kế hoạch 10 năm được xác định qua 3 cột mốc chiến lược:

  • Giai đoạn 1 (1 - 3 năm): Tích hợp kênh tán xạ Raman ($\text{N}_2$, $\text{H}_2\text{O}$) và kênh phân cực Depolarization để đo chính xác tỷ số tán xạ sol khí và nhận dạng giọt nước/tinh thể băng.
  • Giai đoạn 2 (4 - 6 năm): Tự động hóa hoàn toàn hệ Lidar quét 3D di động gắn trên xe chuyên dụng, phát triển mạng lưới 5 trạm quan trắc tự động tại các vùng kinh tế trọng điểm.
  • Giai đoạn 3 (7 - 10 năm): Chế tạo hệ vi Lidar (Micro-Pulse Lidar) thể rắn kích thước nhỏ gọn tích hợp trí tuệ nhân tạo biên (Edge AI) phục vụ cảnh báo thiên tai, giám sát tầng đối lưu và bảo đảm an toàn hàng không quốc gia.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán khoa học và công nghệ phức tạp: Thiết kế, chế tạo và tích hợp thành công hệ Lidar quan trắc sol khí tối ưu đầu tiên tại Việt Nam được kiểm chứng tham số quang điện tử hoàn chỉnh qua mô hình mô phỏng phương trình Lidar trước chế tạo.
  2. Thiết lập quy trình kỹ thuật căn chỉnh cơ - quang đồng trục chính xác cao với hệ gá hai thấu kính chuyên dụng, tối ưu hóa triệt để hàm chồng chập $G(R)$, nâng cao tỷ số tín hiệu trên nhiễu (SNR) và thu hẹp vùng mù quang học cận trường xuống dưới $120\text{ m}$.
  3. Phát triển thành công phần mềm điều khiển thu thập, hiển thị và lưu trữ dữ liệu thời gian thực với độ phân giải không gian siêu cao đạt $1\text{ - }3\text{ m}$, vượt trội so với chuẩn $75\text{ m}$ của mạng lưới viễn thám NASA MPLNET.
  4. Đề xuất thành công giải thuật mới mở rộng biến đổi Haar Wavelet hiệp phương sai (Haar-WCT) với cơ chế tự động tối ưu hóa hệ số giãn nở $a$, loại bỏ hoàn toàn các giả định chủ quan, nâng cao độ chính xác trích xuất ranh giới lớp biên khí quyển trong điều kiện khí hậu nhiệt đới gió mùa phức tạp.
  5. Cung cấp chuỗi dữ liệu thực nghiệm chuẩn xác về cấu trúc sol khí thẳng đứng và chu kỳ nhật động của lớp biên khí quyển tại Hà Nội qua các năm 2012, 2016 và 2018, bóc tách cơ chế tích tụ sol khí trong các đợt nghịch nhiệt mùa đông.
  6. Mở ra ba hướng nghiên cứu học thuật và ứng dụng công nghệ mới: (a) Hệ Lidar Raman đa kênh tích hợp đo nhiệt ẩm và sol khí; (b) Hệ Lidar quét 3D di động kiểm soát ô nhiễm công nghiệp; (c) Hệ thống đồng hóa dữ liệu viễn thám quang học vào mô hình dự báo thời tiết số trị chất lượng cao.