Tổng quan nghiên cứu

Ô nhiễm không khí và sự gia tăng của các hạt sol khí (aerosol) trong khí quyển với kích thước từ 50 nm đến vài micromet đang là thách thức nghiêm trọng đối với sức khỏe con người và biến đổi khí hậu toàn cầu. Để đánh giá chính xác động lực học khí quyển, công nghệ viễn thám chủ động Lidar (Light Detection and Ranging) đóng vai trò then chốt khi cho phép quét không gian từ xa mà không cần tiếp xúc trực tiếp với mẫu khí. Tuy nhiên, các hệ đo Lidar tán xạ đàn hồi truyền thống sử dụng phương pháp một chùm tia luôn đối mặt với bài toán toán học nan giải: một phương trình Lidar duy nhất nhưng chứa tới hai ẩn số chưa biết là hệ số tán xạ ngược và hệ số suy giảm quang học. Việc phụ thuộc vào các mô hình khí quyển chuẩn như US Standard Atmosphere 1976 để tìm nghiệm gần đúng thường gây ra sai số ước tính từ 15% đến 25% trong điều kiện thời tiết thực tế biến động.

Nhằm giải quyết triệt để rào cản trên, luận văn tập trung nghiên cứu và xây dựng cấu hình hệ đo Lidar trong phương pháp hai chùm tia (Two-beam Lidar). Mục tiêu cụ thể là phân tích các hệ thống viễn thám tiên tiến trên thế giới, thiết lập sơ đồ khối chức năng và lựa chọn các cấu kiện quang - điện tử tối ưu, phù hợp với trình độ công nghệ trong nước. Nghiên cứu được triển khai tại Hà Nội trong giai đoạn 2009-2011, kế thừa các kết quả thực nghiệm từ Trung tâm Laser và Môi trường (CODECO) cùng Bộ môn Vật lý, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội. Đề tài mang ý nghĩa khoa học sâu sắc khi triệt tiêu hoàn toàn sự phụ thuộc vào các giả định khí quyển, nâng cao độ chính xác đo đạc với độ phân giải không gian đạt 15 m và dải đo an toàn lên tới 3000 m.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng của thuyết tán xạ quang học và định luật truyền bức xạ điện từ Bouguer - Lambert - Beer ($I = I_0 e^{-\alpha L}$). Khi chùm tia laser truyền trong khí quyển, năng lượng bức xạ bị suy giảm do hai quá trình chính: tán xạ và hấp thụ bởi các phân tử khí và hạt sol khí. Quá trình tán xạ đàn hồi được phân chia thành tán xạ Rayleigh (đối với các phân tử có kích thước nhỏ hơn bước sóng, cường độ tán xạ tỷ lệ nghịch với lũy thừa bậc bốn của bước sóng $\lambda^{-4}$) và tán xạ Mie (đối với các hạt sol khí có kích thước từ 0,1 đến 10 lần bước sóng). Bên cạnh đó, các quá trình tán xạ không đàn hồi Raman và hiệu ứng quang học Doppler cũng được tổng quan để làm rõ tương tác lượng tử và động học của các khối khí.

Mô hình nghiên cứu trọng tâm là phương pháp hai chùm tia của Kunz, Hughes và Paulsons. Phương pháp này vận hành đồng thời hai hệ Lidar tán xạ đàn hồi chiếu ngược chiều nhau vào cùng một thể tích khí quyển: một hệ đo AMALi đặt trên phương tiện bay ở độ cao 3000 m chiếu xuống và một hệ đo KARL đặt tại mặt đất chiếu lên. Tín hiệu thu nhận của hai hệ được mô tả qua hai phương trình Lidar độc lập:

$$S_K(h) = P_K(h)h^2 = C_K \beta(h) T_{[0,h]}^2(h)$$

$$S_A(h) = P_A(h)(H-h)^2 = C_A \beta(h) T_{[h,H]}^2(h)$$

Bằng cách lấy đạo hàm logarit của tỷ số tín hiệu hai chùm tia theo khoảng cách $h$, hệ số suy giảm $\alpha(h)$ được xác định trực tiếp:

$$\alpha(h) = \frac{1}{4} \frac{d}{dh} \left[ \ln \left( \frac{S_A(h)}{S_K(h)} \right) \right]$$

Công thức này chứng minh hệ số suy giảm quang học và hệ số tán xạ ngược hoàn toàn được giải quyết mà không cần bất kỳ giả định gần đúng nào về khí quyển.

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

Nguồn dữ liệu của luận văn được tổng hợp từ các thông số kỹ thuật tiêu chuẩn của hệ thống Lidar di động quốc tế (như hệ AMALi của Viện Alfred Wegener) kết hợp với các phép đo kiểm chứng quang học tại các phòng thí nghiệm chuyên ngành ở Việt Nam. Cỡ mẫu khảo sát bao gồm chuỗi dữ liệu quang học thu được từ 4000 xung laser đơn lẻ, ghi nhận đồng thời ở hai bước sóng cơ bản 1064 nm và bước sóng nhân tần bậc hai 532 nm.

Phương pháp chọn mẫu dựa trên việc khảo sát các xung laser có độ rộng 5 ns với tần số lặp lại 10 Hz và 20 Hz, đảm bảo thu thập đầy đủ tín hiệu phản xạ yếu từ các lớp sol khí ở độ cao từ mặt đất đến 3 km. Phương pháp phân tích được lựa chọn là phân tích quang phổ vi phân kết hợp thuật toán xử lý tín hiệu tương tự và đếm photon. Lý do lựa chọn phương pháp phân tích này là nhằm triệt tiêu hoàn toàn hằng số hệ thống và các yếu tố hình học thông qua phép chia tỷ số quang phổ, từ đó loại bỏ các sai số hệ thống thường gặp trong phương pháp một chùm tia. Toàn bộ lộ trình khảo sát lý thuyết, thiết kế sơ đồ khối và tính toán lựa chọn thiết bị được thực hiện xuyên suốt trong 24 tháng (2009-2011).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, luận văn đã chứng minh phương pháp hai chùm tia giải quyết trọn vẹn bài toán nghịch đảo Lidar, nâng cao độ chính xác tính toán hệ số suy giảm quang học lên hơn 30% so với phương pháp Klett truyền thống của hệ một chùm tia. Phép giải tích toán học khẳng định hệ số suy giảm $\alpha(h)$ và hệ số tán xạ ngược $\beta(h)$ được tách biệt rõ ràng, không còn bị ảnh hưởng bởi tham số hiệu chuẩn ban đầu tại điểm biên.

Thứ hai, nghiên cứu đã thiết kế thành công sơ đồ khối hệ phát quang học hai bước sóng sử dụng nguồn laser Nd:YAG xung cực ngắn. Hệ phát cung cấp mức năng lượng tối đa 360 mJ tại bước sóng 1064 nm và 170 mJ tại bước sóng 532 nm, độ rộng xung đạt 5 ns. Độ phân kỳ của chùm tia được kiểm soát chặt chẽ ở mức dưới 1 mrad (đường kính vết sáng nhỏ hơn 20 mm ở khoảng cách 20 m), kết hợp kính thiên văn Cassegrain đường kính 20,96 cm, tiêu cự 200 cm với thị trường góc mở 2 mrad, giúp loại bỏ trên 85% nhiễu nền ánh sáng ban ngày và triệt tiêu rung động do gió.

Thứ ba, cấu hình hệ thu quang điện tử được xác lập tối ưu cho từng dải sóng chuyên biệt:

  • Kênh hồng ngoại 1064 nm sử dụng đầu thu quang điện thác lũ bán dẫn APD S9251 (kích thước bề mặt 1,5 mm, dải phổ 400-1200 nm, độ nhạy 10 A/W tại điện áp ngược 100-250V) hoạt động ở chế độ tương tự.
  • Kênh khả kiến 532 nm và các kênh tán xạ Raman (607 nm cho Nitơ và 660 nm cho hơi nước) sử dụng các ống nhân quang điện PMT R7400U (dải phổ 185-650 nm, cực đại tại 420 nm, điện áp 1000V) và PMT H6783-20 (dải phổ 300-900 nm, cực đại tại 630 nm). Hệ thống thu thập dữ liệu qua phần mềm chuyên dụng đạt tần số đáp ứng $\ge 250$ MHz và thời gian chết của đầu thu $\le 4$ ns.

Thảo luận kết quả

Nguyên nhân tạo nên sự vượt trội của hệ đo hai chùm tia bắt nguồn từ tính đối xứng quang học khi thực hiện đồng thời hai phép đo đối ngược từ máy bay (hướng xuống) và mặt đất (hướng lên). Khi lấy tỷ số tín hiệu giữa $S_A(h)$ và $S_K(h)$, toàn bộ các hàm truyền khí quyển trung gian và hằng số thiết bị đều được chuẩn hóa, giúp đồ thị phân bố mật độ hạt sol khí thu được phản ánh chính xác cấu trúc thực của các tầng khí quyển. So sánh với các nghiên cứu một chùm tia sử dụng mô hình khí quyển chuẩn US Standard Atmosphere 1976, phương pháp hai chùm tia loại bỏ hoàn toàn các sai số do biến thiên nhiệt độ và áp suất cục bộ gây ra.

Dữ liệu thực nghiệm của hệ thống có thể được biểu diễn trực quan qua biểu đồ đường thể hiện sự suy giảm công suất thu $P(R)$ theo cự ly đo đối với các bước sóng 1064 nm, 532 nm và 355 nm. Đồng thời, bảng tổng hợp thông số kỹ thuật giữa các cấu kiện phát (laser Nd:YAG) và thu (APD S9251, PMT R7400U) cho thấy tính tương thích cao của các thiết bị thương mại được lựa chọn. Kết quả này mang ý nghĩa thực tiễn lớn, chứng minh Việt Nam hoàn toàn có thể tự chủ thiết kế, tích hợp các trạm viễn thám Lidar chất lượng cao với chi phí cạnh tranh so với các hệ thống nhập ngoại nguyên chiếc.

Đề xuất và khuyến nghị

Dựa trên kết quả thiết kế và mô phỏng hệ đo, luận văn đưa ra bốn đề xuất hành động cụ thể nhằm hoàn thiện công nghệ Lidar hai chùm tia:

  1. Gia công và tích hợp khối quang cơ tiêu chuẩn: Tiến hành lắp ráp khung giá đỡ quang học chính xác cao cho kính thiên văn Cassegrain 20,96 cm và hệ laser Nd:YAG 360 mJ, hướng tới mục tiêu hoàn thiện mô hình thử nghiệm thực địa trong thời gian 6 đến 12 tháng, do Trung tâm CODECO và Bộ môn Vật lý chủ trì thực hiện.
  2. Lập trình tự động hóa xử lý dữ liệu thời gian thực: Phát triển thuật toán nghịch đảo hai chùm tia trên môi trường LabVIEW của hãng National Instruments, đặt target metric đạt độ trễ tính toán dưới 0,5 giây cho mỗi profile quét không gian, hoàn thành trong vòng 6 tháng do nhóm kỹ sư điều khiển và xử lý tín hiệu đảm trách.
  3. Mở rộng module thu tán xạ Raman đa kênh: Bổ sung hệ kính lọc giao thoa băng hẹp dưới 0,5 nm tại bước sóng 607 nm và 660 nm để đo đồng thời mật độ hơi nước và sol khí, nâng độ nhạy phát hiện chất ô nhiễm lên 95%, triển khai trong giai đoạn 12 đến 18 tháng bởi các chuyên gia quang lượng tử.
  4. Xây dựng mạng lưới trạm quan trắc Lidar tầng đối lưu: Phối hợp lắp đặt hệ thống Lidar mặt đất kết hợp thiết bị đo trên phương tiện bay không người lái để giám sát bụi mịn PM2.5 và ô nhiễm công nghiệp tại các đô thị lớn ở độ cao dưới 3000 m, lộ trình thực hiện trong 24 tháng do Tổng cục Môi trường và các viện nghiên cứu liên quan phối hợp triển khai.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu chuyên khảo giá trị dành cho các nhóm đối tượng sau:

  1. Học viên cao học và nghiên cứu sinh ngành Quang học, Vật lý laser, Quang điện tử: Nắm vững khung toán học giải phương trình Lidar, kỹ thuật dịch chuyển tần số Raman, hiệu ứng Doppler quang học và nguyên lý phân cực bức xạ laser.
  2. Kỹ sư thiết kế thiết bị quang - cơ - điện tử và viễn thám: Khai thác tài liệu kỹ thuật chi tiết về cách lựa chọn kính thiên văn Cassegrain, gương lưỡng chiết tách chùm, bộ lọc giao thoa băng hẹp 3 nm và mạch tích hợp đầu thu APD/PMT.
  3. Chuyên gia quan trắc khí tượng, thủy văn và ô nhiễm môi trường: Ứng dụng phương pháp hai chùm tia để đo đạc chính xác nồng độ bụi sol khí, mây tầng thấp và phân bố các chất khí độc hại trong tầng đối lưu từ 0 đến 3 km.
  4. Nhà quản lý khoa học công nghệ và đơn vị nghiên cứu ứng dụng: Định hướng chiến lược đầu tư trang thiết bị đo lường quang học hiện đại, phục vụ chương trình nội địa hóa công nghệ viễn thám laser phục vụ an ninh, giao thông và bảo vệ môi trường.

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

1. Phương pháp Lidar hai chùm tia giải quyết khó khăn gì lớn nhất của Lidar truyền thống?

Phương pháp hai chùm tia khắc phục triệt để bài toán "một phương trình hai ẩn số" của Lidar tán xạ đàn hồi truyền thống. Bằng cách sử dụng hai nguồn phát chiếu ngược chiều vào cùng một thể tích khí quyển, hệ thống tạo ra hai phương trình độc lập, giúp tính trực tiếp hệ số suy giảm và tán xạ ngược mà không cần dựa vào các mô hình giả định gần đúng.

2. Vì sao hệ đo cần sử dụng kết hợp cả đầu thu APD và ống nhân quang điện PMT?

Tại bước sóng hồng ngoại 1064 nm, ống nhân quang điện thông thường không nhạy nên bắt buộc phải dùng photodiode thác lũ APD S9251 có độ nhạy 10 A/W. Ngược lại, ở bước sóng khả kiến 532 nm và tán xạ Raman (607 nm, 660 nm) có cường độ rất yếu, đầu thu PMT R7400U với chế độ đếm photon và độ nhạy cực đại tại 420 nm mang lại tỷ số tín hiệu trên nhiễu vượt trội.

3. Thông số phát của laser Nd:YAG trong luận văn có điểm gì nổi bật?

Nguồn laser Nd:YAG được lựa chọn có công suất phát lớn với năng lượng xung đạt 360 mJ ở bước sóng 1064 nm và 170 mJ ở bước sóng 532 nm. Thời gian xung cực ngắn 5 ns kết hợp tần số hoạt động 10-20 Hz cho phép đạt độ phân giải không gian cao lên tới 15 m dọc theo đường truyền tia laser.

4. Việc lựa chọn thị trường góc mở kính thiên văn 2 mrad có ý nghĩa gì?

Thị trường góc mở 2 mrad của kính thiên văn Cassegrain được tính toán để bao trọn chùm tia laser có độ phân kỳ dưới 1 mrad. Thông số này vừa đảm bảo thu trọn vẹn toàn bộ bức xạ tán xạ ngược hữu ích, vừa giảm thiểu tối đa ánh sáng nền ban ngày và hạn chế ảnh hưởng của rung động cơ học trong điều kiện gió mạnh.

5. Hệ đo Lidar hai chùm tia có khả năng hoạt động ở dải độ cao nào?

Hệ thống được thiết kế tối ưu cho tầng đối lưu ở độ cao từ mặt đất đến 3000 m (phù hợp với trần bay an toàn của thiết bị đo AMALi). Dải đo này bao phủ toàn bộ tầng biên khí quyển, nơi tập trung phần lớn các hoạt động phát thải sol khí, hạt bụi và khí nhà kính.

Kết luận

  • Đóng góp lý thuyết: Luận văn đã phân tích toàn diện cơ sở toán học của phương pháp hai chùm tia, khẳng định tính khả thi của việc loại bỏ hoàn toàn các mô hình giả định khí quyển để xác định trực tiếp hệ số suy giảm và tán xạ ngược.
  • Thiết kế khối phát chuẩn xác: Xây dựng cấu hình nguồn laser Nd:YAG phát đồng thời hai bước sóng 1064 nm (360 mJ) và 532 nm (170 mJ) với độ rộng xung 5 ns và độ phân kỳ chùm tia dưới 1 mrad.
  • Tối ưu hóa hệ thu quang điện: Lựa chọn đồng bộ kính thiên văn Cassegrain 20,96 cm, các phin lọc giao thoa băng hẹp cùng hệ đầu thu APD S9251 và PMT R7400U/H6783-20 cho phép hoạt động linh hoạt ở cả chế độ tương tự và đếm photon.
  • Tính khả thi và ứng dụng thực tiễn: Đề xuất cấu hình thiết bị thương mại có chất lượng trung bình khá so với thế giới, chi phí hợp lý và hoàn toàn tương thích với điều kiện nghiên cứu tại Việt Nam.
  • Lộ trình kế tiếp: Hoàn thiện tích hợp hệ thống phần cứng và thuật toán điều khiển trên LabVIEW trong giai đoạn 12 tháng tiếp theo để đưa hệ thống vào thử nghiệm thực tế tại các vùng khí quyển trọng điểm.