Nghiên cứu sự thay đổi tỷ số Lidar theo khoảng cách

Nội dung text: Nghiên cứu sự thay đổi tỷ số Lidar theo khoảng cách

1. T¹p chÝ KHKT Má - §Þa chÊt, sè 53, 01/2015, tr.81-84 KHOA HỌC CƠ BẢN (trang 81÷84) NGHIÊN CỨU SỰ THAY ĐỔI TỶ SỐ LIDAR THEO KHOẢNG CÁCH BÙI THỊ THANH LAN, Trường Đại học Mỏ - Địa chất ĐỒNG THỊ LINH, Đại học KTCN Thái guyênN Tóm tắt: LIDAR là công nghệ mới dùng laser để khảo sát đối tượng ở phạm vi xa. Hiện nay thuật toán nghịch đảo Klett là phương pháp phổ biến để giải phương trình LIDAR nhằm xác định các đại lượng cần đo như hệ số tán xạ ngược, hệ số suy giảm. Vì phương trình LIDAR chứa hai ẩn nên cần có tỷ số LIDAR để đủ điều kiện giải phương trình. Người ta có thể sử dụng các giá trị tỷ số LIDAR cho trước nằm trong khoảng 20-90 sr. Như vậy sẽ có sai số vì trong thực tế tỷ số LIDAR không phải là hằng số mà thay đổi theo thời gian và khoảng cách. Trong phương pháp hai chùm tia, người ta tích hợp hai thiết bị LIDAR. Hai thiết bị đồng thời thực hiện hai phép đo độc lập với cùng một mẫu khí quyển ở giữa chúng và do đó cung cấp hai phương trình LIDAR độc lập với nhau. Giả thiết ánh sáng laser phát ra là đơn sắc và các quá trình tán xạ phức tạp được bỏ qua. Sử dụng phương pháp hai chùm tia xácị đ nh được hệ số suy giảm và hệ số tán xạ ngược mà không cần một giả thuyết hay một sự hiệu chuẩn nào. Hệ số suy giảm và hệ số tán xạ ngược đều phụ thuộc khoảng cách. Thực nghiệm tính sự phụ thuộc vào khoảng cách của hệ số suy giảm cho thấy có tồn tại đám sol khí ở trong vùng từ 0,3 km đến 0,5 km. Theo tính toán từ số liệu thực nghiệm hệ số tán xạ ngược giảm ở vùng từ 0,3 km đến 0,6 km. Như vậy kết quả tính của hai hệ số suy giảm và tán xạ ngược là khá phù hợp. Theo công thức tính tỉ số LIDAR không phải là hằng số. Tỉ số LIDAR phụ thuộc vào đặc trưng vật lý của hạt sol khí như loại hạt, kích thước, chiết suất. Các yếu tố trên đều thay đổi theo khoảng cách, do vậy tỉ số LIDAR cũng luôn thay đổi theo khoảng cách. Theo tính toán thực nghiệm tỷ số LIDAR phụ thuộc vào khoảng cách và có cực đại bằng 46sr tại 0,4 km, cực tiểu bằng 17,5sr tại 0,18 km phù hợp với lý thuyết sol khí vùng biển. 1. Phương pháp hai chùm tia Tuy nhiên, đây là một phương pháp khá tốn Trong phương pháp một chùm tia LIDAR, kém và mất nhiều thời gian vì phương pháp này để giải quyết bài toán một phương trình LIDAR cần thiết bị thu chính xác để đo sự dịch chuyển nhưng chứa hai ẩn số là hệ số suy giảm và hệ số bước sóng trong tán xạ Raman không đàn hồi tán xạ ngược, người ta buộc phải đưa thêm giả và tiết diện tán xạ Raman phải được xác định thiết về mối quan hệ giữa hai ẩn số này. Thông chính xác qua tính toán lý thuyết và thực thường, ta tính toán tỉ số giữa hệ số suy giảm và nghiệm đối với nhiều bước sóng laser khác hệ số tán xạ ngược (tỉ số LIDAR) ở một vị trí nhau. Vì vậy, phương pháp hai chùm tia cho cho trước sau đó giả thiết rằng tỉ số này không phép tính trực tiếp tỉ số LIDAR đã giải quyết thay đổi trong cả phạm vi đo đạc. Tuy nhiên, được những khó khăn trên. phép tính trên sẽ gặp sai số lớn khi môi trường Để thực hiện phép đo các thông số vật lý sử không đồng nhất hay phạm vi đo đạc lớn. Một dụng phương pháp hai chùm tia, người ta tích phương pháp rất tinh tế về mặt vật lý để đo tỉ số hợp hai thiết bị Airborne LiDAR ( kí hiệu là A) LIDAR là kỹ thuật đo tán xạ Raman trên phân và Ground - Based LIDAR ( kí hiệu là K). Thiết tử Nitơ do Ansmann et al (1992) áp dụng [2]. bị thứ nhất được gắn trên một máy bay bay ở độ 81
2. cao H hướng xuống phía dưới và đo sự tăng dần 1.2. Tỷ số LIDAR trong phương pháp hai của mật độ sol khí và các phân tử khí. Ngược chùm tia LIDAR lại, thiết bị thứ hai đặt trên mặt đất hướng lên Giả sử ở một độ cao Rref cho trước, ta đã trên đo sự giảm dần của mật độ sol khí theo độ biết giá trị của hệ số tán xạ ngược ref. Khi đó cao [1, 4]. Hai thiết bị có thể đặt cùng độ cao và tín hiệu thu được của hai hệ LIDAR được viết được chiếu vào môi trường sol khí từ hai phía. lại như sau [5]: Hai thiết bị trên đồng thời thực hiện hai Rref SRCRrdr()()exp(2()) . (7) phép đo độc lập với cùng một mẫu khí quyển ở KrKr efef giữa chúng và do đó cung cấp hai phương trình 0 LIDAR độc lập với nhau. Giả thiết ánh sáng Sau khi biến đổi chúng ta thu được hệ số laser phát ra là đơn sắc và các quá trình tán xạ tán xạ ngược có biểu thức như sau: phức tạp được bỏ qua [3]. Khi đó phương trình SRSRKA()() ()()RR ref , (8) LIDAR tán xạ đàn hồi mô tả tín hiệu nhận được SRSRKrAr()()efef như là một hàm của các thông số khí quyển và Giá trị ref : thông số hệ thống. SRSRKrAr()()efef . (9) 1.1. Phương trình LIDAR cho phương pháp ()Rref H hai chùm tia CCrdrexp2()  KA Từ phương trình LIDAR cơ bản với giả 0 thiết tín hiệu đo được là rõ nét nhất và hàm Như vậy, ta có tỷ lệ LIDAR: Overlap O(R) = 1, ta có phương trình LIDAR 1 d SRA () ln( ) mol R của hai hệ LIDAR như sau [5]: 4(dRSR ) aer ()R K 22 LRaer () SRPRKKK()()()() RCR TR  [0,R] , (1) aer ()R SRSRKA( )( ) 22 ()(RRrmolef ) SRPRAAA()()()() RCR TR  [H,R] , (2) SRSRKrAr()()efef R là khoảng cách giữa hệ LIDAR và mẫu (10) khí quyển cần đo, β(R) là hệ số tán xạ ngược, 2. Khảo sát tỷ số LIDAR theo khoảng cách CA và CK lần lượt là hằng số LIDAR của hai hệ 2.1. Thực nghiệm tương ứng. Hàm truyền T có dạng: Khảo sát khí quyển gần biển Point Loma R Peninsula - Mỹ [6] bằng hai máy LIDAR đặt cách . (3) TRrdr()exp(2())  nhau 0.9825 km, Hughes và Paulson đã sử dụng 0 nguồn laser Nd :YAG. Hai thiết bị LIDAR đặt ở Thay (3) vào (1) và (2), để đơn giản, ta kí độ cao khoảng 38m trên mực nước biển, phép đo hiệu rK = r: R được tiến hành trong điều kiện thời tiết tốt, không SRCRrKK()()exp(2( dr  )) . (4) khí ít ô nhiễm [5]. Kết quả thu được cho trong 2 0 bảng 1, tín hiệu thu là S(R) = ln(P(R).R ). R SRCRrdr()()exp(2())  . (5) AAAA LiDAR A H Sau một số biến đổi có hệ số suy giảm: 1 d SR() (R ) ln A . (6) (9) H 4dR SK ( R ) Hệ số suy giảm thu được bằng cách trên R Mẫu đo không cần một giả thuyết hay một sự hiệu chuẩn nào. Sai lệch của (R) phụ thuộc vào tín hiệu nhiễu của cả hai hệ LIDAR, vì vậy sẽ kết quả tính chính xác hơn nếu hai hệ LIDAR có cùng tỉ số S/N (signal to noise). O LiDAR K Hình 1. Sơ đồ bố trí hai hệ LIDAR 82
3. Bảng 1. Kết quả khảo sát khí quyển tại Point Loma Peninsula - Mỹ R(km) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 S1(r) 4,00 4,35 4,92 5,25 6,6 6,42 6,21 6,00 6,95 S2(r) 6,15 5,66 5,87 5,70 4,82 4,63 4,22 3,81 4,25 Sử dụng phương pháp hai chùm tia, chúng tôi đã xác định được hệ số suy giảm và hệ số tán xạ ngược. Kết quả chúng tôi thu được các đồ thị 2 và 3. Hình 2. Sự phụ thuộc của hệ số suy giảm theo khoảng cách Ta thấy sự phụ thuộc vào khoảng cách của hệ số suy giảm đo bằng phương pháp hai chùm tia cho giá trị cực đại ở khoảng cách 0,4 km, và giảm nhanh trong bán kính khoảng 0,1 km. Như vậy có thể kết luận là có tồn tại đám xon khí ở trong vùng từ 0,3 km đến 0,5 km. Hình 3. Sự phụ thuộc của hệ số tán xạ ngược theo khoảng cách Trên hình 3 cho thấy hệ số tán xạ ngược giảm ở vùng từ 0,3 km đến 0,6 km. Như vậy chứng tỏ có sự tồn tại của đám xon khí trong vùng này. Như vậy kết quả tính của hai hệ số suy giảm và tán xạ ngược là khá phù hợp. 2.2. Xác định tỉ số LIDAR theo khoảng cách Bảng 2. Tỉ số LIDAR L theo khoảng cách R(km) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 L (sr) 43,07 18,68 34,06 48,47 37,18 36,47 33,13 38,44 83
4. Hình 4. Sự phụ thuộc của tỉ số LIDAR vào khoảng cách Nhận xét: Theo tính toán tỷ số LIDAR phụ tính các thông số khí quyển và ứng dụng thuộc vào khoảng cách cho thấy có cực đại phương pháp hai chùm tia ở Việt Nam. bằng 46 sr tại 0,4 km và cực tiểu là 17,5 sr tại Bài báo hoàn thành nhờ sự hỗ trợ của đề tài 0,18 km. Kết quả tính toán cho thấy môi trường KHCN cấp Tỉnh mã số 10/2014/HĐ-KHCN. không khí giữa hai thiết bị LIDAR không đồng nhất. Thông thường, đối với sol khí vùng biển TÀI LIỆU THAM KHẢO và vùng dân cư, tỉ số LIDAR thường có giá trị nhỏ hơn 50 sr [1,5]. Như vậy, các số liệu tính [1]. Ackermann, 1998. The extinction to toán tỉ số LIDAR bằng phương pháp hai chùm backscatter ratio of tropospheric aerosol. A tia như trên là hoàn toàn phù hợp. numerical study, J. Atmos. Ocean. Technol., 15, Trong bảng 2 ta thấy tỉ số LIDAR không 1043–1050. phải là hằng số. Tỉ số LIDAR phụ thuộc vào [2]. Ansmann, A., Wandinger, U., Riebesell, đặc trưng vật lý của hạt sol khí như loại hạt (tan M., Weitkamp, C., and Michaelis, W, 1992. trong nước, không tan trong nước, hạt loại Independent measurements of extinction and cacbon ), phân bố kích thước, chiết suất. Các backscatter profiles in Cirrus clouds by using a yếu tố trên đều thay đổi theo khoảng cách, do combined Raman elastic-backscatter Lidar. vậy tỉ số LIDAR cũng luôn thay đổi theo Appl. Optics, 31, 7113–7131. khoảng cách. [3]. Bui Thi Thanh Lan. Introduction the 4. Kết luận LIDAR two-stream technique”, Proceeding the Hiện nay, kĩ thuật LIDAR đã và đang được conference“ Advance in mining and tunneling. ứng dụng trong các phép đo đánh giá những 8/2012. thông số của khí quyển. Các phương pháp tính [4]. I.S. Stachlewska and Ritter, 2010. On toán cũng được phát triển để cho kết quả chính retrieval of lidar extinction profiles using Two – xác hơn. Qua bài báo này, chúng tôi đã bước stream and Raman techniques. Atmos. Chem. đầu nghiên cứu phương pháp tính tỉ số LIDAR Phys, 10, PP 2813 – 2824. bằng phương pháp hai chùm tia dựa trên kết quả [5]. Hughes HG, Paulson, M.R, 1998. Double – thực nghiệm đã có. Hiện nay chúng tôi chưa có end lidar techniques for aerosol studies. Applied đủ điều kiện thực nghiệm ứng dụng phương Optics, 27, PP 2273 – 3278. pháp hai chùm tia. Đây là cơ sở lý thuyết để (xem tiếp trang 80) 84
5. ABSTRACT Study the depending of lidar ratio on distance Bui Thi Thanh Lan, Hanoi University of Mining and Geology Dong Thi Linh, Thai Nguyen University of Technology LIDAR is a new technology for remote servey. Klette inversion algorithm is commontly used to solve LIDAR equation to determind the measurement parameters, such as backscattering efficient and extinction efficient. Because of the LIDAR equation has two aerosol unknown parameters, that is why the LIDAR ratio is needed. One can choose LIDAR ratio between from 20 - 90 sr. Therefore, it makes errors because of LIDAR ratio dependent on time and location. In the two stream methode two LIDAR instruments are integrated. Two instruments independendly measure with the same aerosol sample at the same time then giving two independend LIDAR equations. Assuming that laser is monochrome and complicated scattering is overlooked. Using the two stream methode to determind the efficient extinction and back scattering cofficient without any assuming or calibration. The efficient extinction and the backscattering dependend on the distance. The experiments calculating the dependence of the efficient extinction on the distance can conclude that there is an aerosol dame in distance from 0,3 km to 0,5 km. Caculating from experimental data get the backcattering efficient reduce in the distance from 0,3 km to 0.6 km. Therefore, the results of two efficient are well fit to each other. By the fomular, the LIDAR ratio is not constant. The LIDAR ratio dependends on the characteristics of the aerosol partical such as type, dimension and refractive index. These charateristics change by distance and the LIDAR ratio also depends on the distance. Experimental results showed that the LIDAR ratio depends on the distance and has maximum value 46 sr at 0.4 km and minimum value 17.5sr at 0.18 km so that it is fiting with theory of sea aerosol. 85