Tổng quan luận án
Bản luận án tiến sĩ kỹ thuật của nghiên cứu sinh Lai Thị Vân Quyên tập trung vào việc nghiên cứu, làm chủ và nâng cao độ chính xác của phương pháp đo kích thước và vận tốc hạt mưa bằng hiệu ứng quang học. Đề tài được triển khai trong khuôn khổ đề tài nghiên cứu khoa học cấp nhà nước do Bộ Công Thương tài trợ kinh phí: "Tiếp nhận, làm chủ công nghệ đo kích thước và vận tốc hạt bằng phương pháp quang học ứng dụng trong hệ thống quan trắc lượng mưa phục vụ nghiên cứu chống biến đổi khí hậu", với sự hỗ trợ về trang thiết bị và nhân lực của Viện Nghiên cứu Điện tử, Tin học, Tự động hóa phối hợp cùng các chuyên gia thuộc Trường Đại học Bách khoa Saint Petersburg (Liên bang Nga).
Tính cấp thiết và khoảng trống nghiên cứu
Đo kích thước hạt mưa và xác định phổ phân bố kích thước hạt mưa (DSD - Drop Size Distribution) giữ vai trò then chốt trong nhiều lĩnh vực: đánh giá sự suy giảm và méo tín hiệu trong lan truyền sóng điện từ, nghiên cứu vật lý khí quyển, cấu trúc tầng đối lưu, vi vật lý mây, cũng như dự báo xói mòn và trượt lở đất trong canh tác nông nghiệp. Thiết bị đo mưa tự động thông dụng hiện nay là loại chao lật chỉ đo được lượng mưa tích lũy ($R$), cường độ mưa ($I$) và thời gian mưa, không thể xác định kích thước ($D$) và vận tốc ($v$) của từng hạt mưa.
Các thiết bị đo quang học hai dải sáng hiện hành theo mô hình của D. Kiesewetter và V. I. Malyugin (2004, 2009) tính toán đường kính hạt dựa vào độ sâu điều chế $M$. Phương pháp này bộc lộ những khoảng trống kỹ thuật lớn:
- Độ sâu điều chế $M$ phụ thuộc vào giá trị cực đại và cực tiểu của tín hiệu dòng quang điện nên rất nhạy cảm với nhiễu biên độ;
- Chưa có giải pháp xử lý khi hai chồi xung quang điện có biên độ cực đại không bằng nhau;
- Khi đo các hạt có kích thước nhỏ hơn khoảng cách giữa hai khe nhạy sáng, điểm cực tiểu của đáy xung giữa hai chồi bị triệt tiêu về mức không, khiến phương pháp tính theo độ sâu điều chế $M$ không còn sử dụng được;
- Hệ thống đo nguyên mẫu sử dụng bước sóng hồng ngoại 850 nm gây khó khăn cho việc cân chỉnh quang học ngoài thực địa, cơ cấu gá cố định khó hiệu chỉnh trục và hệ thống phụ thuộc hoàn toàn vào máy tính cùng khối biến đổi ADC gắn ngoài cồng kềnh.
Mục tiêu nghiên cứu
- Về mặt khoa học: Đề xuất biểu thức toán học và thuật toán xử lý dữ liệu mới để tính toán đường kính tương đương ($D$) và vận tốc ($v$) của hạt mưa từ các xung quang điện thu được, khắc phục hạn chế của phương pháp độ sâu điều chế và nâng cao độ chính xác của phép đo.
- Về mặt công nghệ: Hoàn thiện công nghệ thông qua việc cải tiến phần cứng: thay thế nguồn sáng, thiết kế lại cơ cấu gá đỡ điều chỉnh trục quang và phát triển mạch vi xử lý nhúng tích hợp tính toán thời gian thực nhằm thay thế hoàn toàn máy tính xử lý trung tâm.
- Về mặt thực tiễn: Chế tạo thành công thiết bị đo hoàn chỉnh có khả năng hoạt động độc lập ngoài trời, tích hợp vào mạng lưới quan trắc khí tượng tự động và ứng dụng dữ liệu thu được vào bài toán đánh giá động năng mưa phục vụ tính toán xói mòn đất.
Đối tượng, phạm vi và giới hạn nghiên cứu
- Đối tượng nghiên cứu: Bài toán đo kích thước tương đương và vận tốc của hạt mưa bằng phương pháp phân tích hai xung quang điện.
- Phạm vi nghiên cứu: Tập trung vào kỹ thuật xử lý tín hiệu sau khử nhiễu từ bộ thu quang, cấu trúc khe nhạy sáng đôi và thiết kế thiết bị phần cứng nhúng kế thừa từ mô hình của D. Kiesewetter và V. I. Malyugin.
- Giới hạn bài toán: Dải đo đường kính hạt khảo sát từ 0,5 mm đến 6 mm (dải hạt mưa phổ biến trong tự nhiên); không xét đến nhiễu môi trường ngẫu nhiên như côn trùng xâm nhập khoảng đo hoặc biến đổi độ ẩm không khí cực đoan; không đánh giá định lượng sai số vận tốc do thiếu thiết bị đối chuẩn vận tốc chuyên dụng tại thời điểm thực hiện.
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa các phương pháp đo kích thước và vận tốc hạt mưa từ giai đoạn thủ công đến các thiết bị tự động hiện đại:
| Nhóm phương pháp |
Công trình tiêu biểu |
Nguyên lý hoạt động |
Ưu điểm |
Hạn chế chính |
| Thủ công |
Lowe (1892) [39]; Bentley (1904) [9]; Fuchs & Petrjanoff (1937) [21] |
Vết giấy thấm nhuộm màu, viên bột sàng mịn, phương pháp ngâm giọt trong dầu kỵ nước. |
Thiết bị đơn giản, chi phí ban đầu thấp. |
Tốn thời gian, không có tính thời gian thực, sai số lớn, không đo được vận tốc rơi. |
| Cơ học (Disdrometer tác động) |
Joss & Waldvogel (1967) [51] (JWD); Tokay và cs (2001); Krajewski và cs (2006) [34]; Thurai và cs (2013) [6, 19] |
Hạt mưa rơi tác động vào màng cảm biến làm dịch chuyển lõi từ trong cuộn cảm, biến đổi xung năng lượng điện từ. |
Cấu tạo gọn hơn hệ quang, vận hành tự động cơ bản. |
Rung liên tục khi mưa to hoặc có gió mạnh gây sai số xung nhiễu lớn; hạt vỡ khi va chạm; hoàn toàn không đo được vận tốc hạt. |
| Xử lý hình ảnh |
Kruger & Krajewski (2002) [33]; Schönhuber và cs (2007) [57] (2DVD); Newman và cs (2009) (SVI) |
Hai camera tốc độ cao ghi hình chiếu vuông góc của giọt mưa rơi qua vùng quan sát. |
Đo trực quan kích thước, vận tốc và chi tiết hình học của hạt. |
Xử lý ảnh phức tạp, thời gian phân tích dài; khó phân tách các hạt dính nhau; thiết bị ngoài trời cồng kềnh nối máy tính lớn; giá thành cao gấp 3-5 lần thiết bị quang. |
| Quang học một dải sáng (1 xung) |
Illingworth & Stevens (1987) [30]; Grossklaus và cs (1998) [26]; Löffler-Mang & Joss (2000) [38]; Bloemink & Lanzinger (2005) [8] (Parsivel, Thies LPM) |
Chiếu một chùm sáng phẳng hoặc hình khuyên/trụ; đo mức suy giảm điện áp hoặc độ rộng và biên độ cực đại của một xung quang điện. |
Thời gian đáp ứng nhanh, đo đồng thời kích thước và vận tốc, kích thước gọn nhẹ. |
Kích thước suy ra từ biên độ đỉnh đơn lẻ dễ bị nhiễu; mật độ năng lượng chùm sáng không đều khiến hạt đi qua tâm bị đo quá to, qua rìa bị đo quá nhỏ (Fransson và cs, 2011 [13]). |
| Quang học hai dải sáng (2 xung) |
Cloos (2007) [11]; Winsky (2012) [6]; Kiesewetter & Malyugin (2004, 2009) [14, 16] |
Chiếu chùm sáng qua hệ hai khe nhạy sáng hoặc tạo hai dải sáng song song, ghi nhận dạng xung hai chồi. |
Giảm thiểu sai số do vị trí rơi không đều, ít phụ thuộc vào suy hao môi trường truyền dẫn, độ tin cậy theo Luật số Lớn cao hơn. |
Mô hình của Winsky dùng laser diode (LD) bị nhiễu đốm quang học, cần máy tính xử lý thuật toán cosine phức tạp; mô hình Kiesewetter dùng độ sâu điều chế $M$ nhạy nhiễu đỉnh xung và giới hạn dải hạt nhỏ. |
Khoảng trống luận án lựa chọn giải quyết
Luận án lựa chọn kế thừa mô hình đo hai xung quang điện của D. Kiesewetter và V. I. Malyugin nhưng tập trung giải quyết triệt để các tồn tại về giải thuật và phần cứng:
- Thay thế việc xác định đường kính từ độ sâu điều chế $M = I_{max} - I_{min}$ bằng giải thuật trích xuất tham số tỷ lệ $k_{075}$ tại mức 75% của biên độ xung nhằm loại bỏ ảnh hưởng của việc đỉnh xung bị biến dạng hoặc đáy xung bị bão hòa triệt tiêu.
- Chuyển đổi mô hình thực nghiệm sử dụng máy tính và khối ADC rời sang cấu trúc bo mạch xử lý nhúng chuyên dụng, cho phép tính toán kích thước, vận tốc tại chỗ theo thời gian thực.
- Cải tiến cơ cấu cơ quang học: dùng nguồn phát LED thích hợp, thiết kế hệ gá điều chỉnh trục quang bằng cơ cấu lò xo đàn hồi vi chỉnh giúp định tâm quang học chính xác ngoài hiện trường.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết quang học và mô hình toán học tín hiệu
Khi hạt mưa bay cắt ngang qua chùm ánh sáng song song, sự suy giảm năng lượng ánh sáng tới cảm biến quang là kết quả của các hiện tượng tương tác vật lý: phản xạ, khúc xạ, tán xạ rìa và hấp thụ. Vì hạt mưa trong khí quyển có tiêu cự quang học rất ngắn, trong bài toán thu năng lượng trên cảm biến đặt tại vùng nhiễu xạ Fresnel, hạt được xem như một đĩa chắn sáng phẳng tương đương với hình cầu có cùng đáp ứng quang học.
Khoảng cách $z$ từ màn chắn khe nhạy sáng đến tâm vùng đo phải thỏa mãn điều kiện nhiễu xạ Fresnel:
$$z_{max} \ge z \ge \frac{D_{min}^2}{\lambda}$$
Trong đó:
- $\lambda$ là bước sóng của nguồn sáng phát;
- $D_{min}$ là đường kính hạt nhỏ nhất cần đo;
- $z_{max}, z_{min}$ là các giới hạn khoảng cách hình học của buồng đo.
Trong mô hình khe đôi có độ rộng hai khe bằng nhau $w_1 = w_2 = w$ và khoảng cách giữa hai khe là $g$, tín hiệu dòng quang điện $I(X_C)$ biến thiên theo tọa độ dịch chuyển $X_C$ của tâm hạt. Khi hạt có kích thước $D < w_1 + g + w_2$, tín hiệu thu được có dạng xung hai chồi riêng biệt với hai đỉnh cực đại tương ứng khi hạt đi qua lần lượt hai khe và một đáy cực tiểu khi hạt nằm ở khoảng giữa hai khe.
Phương pháp nghiên cứu thực nghiệm và trích xuất dữ liệu
- Mô hình hóa hình học giọt hạt: Sử dụng các viên bi sắt chuẩn (được kiểm tra kích thước bằng thước panme cơ khí có độ phân giải cao) để làm vật chuẩn thực nghiệm mô phỏng hạt mưa dạng vật đen tuyệt đối, loại bỏ sai số do biến dạng giọt nước trong quá trình hiệu chuẩn tham số.
- Trích xuất tín hiệu xung quang điện: Tín hiệu từ điốt quang (photodiode) sau khi qua mạch tiền xử lý (khuếch đại thuật toán và lọc thông dải để khử trôi dạt đường nền quang) được đưa vào dao động ký số (oscilloscope) để ghi nhận hình dạng xung thực tế, đo đạc các giá trị rời rạc $i_{x_y}$ (giá trị thời gian tại mức $x%$ của biên độ xung lớn nhất trên sườn xung thứ $y$).
- Phương pháp thực nghiệm đa cấp:
- Thử nghiệm kiểm tra độ phân giải với bi sắt chuẩn từ 1 mm đến 10 mm;
- Thử nghiệm với các giọt chất lỏng rơi tự do tạo ra từ các kim nhỏ giọt y tế chuẩn;
- Thử nghiệm với hệ thống vòi tạo mưa nhân tạo đối chứng với máy đo mưa chao lật chuẩn thương mại The Rain Collector II (hãng Davis, Hoa Kỳ);
- Thực nghiệm quan trắc trận mưa tự nhiên ngoài trời tại Trạm Khí tượng Hà Đông.
Nội dung chính theo từng chương
Chương I. Tổng quan các phương pháp đo kích thước hạt mưa
Chương I phân tích toàn diện lịch sử hình thành và phát triển của các công nghệ đo thông số giọt mưa. Nội dung chương chia tách và đánh giá chi tiết:
- Phân tích cơ chế suy giảm quang học: Khẳng định kích thước đo được bằng phương pháp quang học là "đường kính tương đương hình cầu" dựa trên năng lượng che khuất quang học tương đương.
- So sánh các hệ thống thương mại: Hệ thống JWD chịu ảnh hưởng nặng nề bởi tiếng ồn cơ học và dao động màng; hệ thống 2DVD và SVI đạt độ chi tiết cao về hình ảnh nhưng chi phí cao, cồng kềnh; các hệ thống một dải sáng phổ biến như OTT Parsivel (thế hệ 1 và 2) và Thies LPM gặp vấn đề cố hữu về sự bất đồng nhất của chùm tia laser trên diện tích thu nhận.
- Đánh giá sâu hai hướng tiếp cận đo bằng hai xung quang điện: Mô hình của Michael Peter Cloos (2007) và Bryson Evan Winsky (2012) tại Đại học Iowa dùng hai dải laser diode 650 nm, yêu cầu tần số lấy mẫu 66,6 kHz và thuật toán xấp xỉ hàm cosine để xử lý độ bằng đầu của xung $l_d$. Ngược lại, mô hình của D. Kiesewetter và V. I. Malyugin (2004, 2009) tính toán theo độ sâu điều chế $M$ có cấu trúc quang học tối ưu hơn nhưng còn vướng các hạn chế về mặt lý thuyết xử lý biên độ và tính khả thi ngoài hiện trường. Từ đó, chương I xác lập mục tiêu nghiên cứu và giới hạn phạm vi đề tài của luận án.
Chương II. Nghiên cứu, nâng cao độ chính xác của phương pháp đo vận tốc, kích thước hạt mưa
Chương II là phần đóng góp cốt lõi về lý thuyết toán học và giải pháp kỹ thuật công nghệ của tác giả:
1. Đề xuất biểu thức và thuật toán xử lý dữ liệu mới
Thay vì sử dụng độ sâu điều chế $M$, tác giả phân tích dạng hình học của xung hai chồi đảo ngược và đề xuất tham số hình học $k_{075}$. Tham số này được xác định từ tỷ lệ khoảng cách thời gian giữa các điểm cắt mức 75% biên độ cực đại trên các sườn xung:
$$k_{075} = f(i_{75_1}, i_{75_2}, i_{75_3}, i_{75_4})$$
Việc lấy ngưỡng tại mức 75% biên độ giúp loại bỏ hoàn toàn hiện tượng méo do nhiễu tại đỉnh xung (xảy ra ở mức 100%) và không bị phụ thuộc vào mức đáy cực tiểu (vốn bị triệt tiêu khi kích thước hạt nhỏ hơn khoảng cách khe $g$).
Biểu thức xác định vận tốc hạt $v$ dựa trên khoảng cách giữa hai khe nhạy sáng $g$ và khoảng thời gian giữa hai đỉnh chồi xung $\tau_i$:
$$v = \frac{g + w}{\tau_i}$$
Đối với các hạt kích thước lớn gây chập hai xung thành một xung duy nhất, vận tốc được tính thông qua độ rộng xung kết hợp với hệ số hiệu chỉnh hình học.
2. Cải tiến cấu trúc phần cứng và thiết kế chế tạo thiết bị
Nhằm chuyển hóa từ mô hình phòng thí nghiệm sang thiết bị công nghiệp ngoài trời, tác giả đã tiến hành các cải tiến công nghệ đồng bộ:
- Nguồn phát quang: Thay thế diode laser bằng đèn LED có bước sóng thích hợp trong dải khả kiến kết hợp thấu kính chuẩn trực, triệt tiêu hoàn toàn hiện tượng nhiễu đốm (speckle noise) và tạo điều kiện thuận lợi cho việc quan sát, cân chỉnh trục quang bằng mắt thường.
| Đặc tính quang học |
Diode Laser (LD) [63] |
Điốt phát quang (LED) [63] |
Giải pháp của luận án |
| Độ kết hợp (Coherence) |
Rất cao, dễ tạo nhiễu đốm quang |
Thấp, chùm sáng đồng đều |
Sử dụng LED để loại bỏ nhiễu đốm trên cảm biến |
| Hiệu chỉnh quang trục |
Khó quan sát nếu dùng tia cận hồng ngoại |
Thuận lợi trong dải quang phổ chọn lọc |
Cho phép cân chỉnh trực quan hiện trường |
| Độ ổn định nhiệt |
Độ nhạy nhiệt cao, trôi bước sóng |
Tương đối ổn định |
Tích hợp mạch ổn dòng cho LED |
- Cơ cấu gá đỡ cơ quang học: Bổ sung cơ cấu lò xo đàn hồi đa điểm giúp điều chỉnh vi cấp vị trí nguồn phát và cảm biến so với trục quang của thấu kính hội tụ, khắc phục tình trạng kẹt cơ khí và mất đồng trục khi nhiệt độ môi trường thay đổi.
- Phát triển phần cứng xử lý nhúng: Luận án thiết kế chi tiết 3 khối mạch in chuyên dụng (PCB):
- Mạch phát quang: Ổn định dòng cấp cho nguồn LED;
- Mạch thu quang: Sử dụng photodiode có tốc độ đáp ứng cao, mạch khuếch đại chuyển đổi dòng - áp (TIA) và tầng lọc thông dải loại bỏ can nhiễu quang nền;
- Mạch CPU xử lý trung tâm: Tích hợp bộ vi điều khiển tốc độ cao, trực tiếp số hóa tín hiệu qua ADC nội, thực thi thuật toán phân loại dạng xung, tính toán $D$ và $v$, xuất dữ liệu qua giao tiếp nối tiếp tiêu chuẩn.
- Quy trình hiệu chuẩn: Xây dựng quy trình hiệu chuẩn hai bước: hiệu chỉnh cứng (căn chỉnh cơ quang học) và hiệu chỉnh mềm (bù trừ đường cong thực nghiệm cho các hạt $D < 3,5\text{ mm}$).
Chương III. Xây dựng mô hình thực nghiệm, phân tích và đánh giá kết quả đo đạc trên mô hình đề xuất
Chương III trình bày quá trình thực nghiệm toàn diện nhằm chứng minh tính đúng đắn của giải thuật và độ tin cậy của thiết bị chế tạo:
[Thử nghiệm Bi sắt chuẩn] -> [Thử nghiệm Kim nhỏ giọt] -> [Thử nghiệm Mưa nhân tạo] -> [Thực địa Trạm Hà Đông]
(1 mm ÷ 10 mm qua Panme) (Thay đổi tốc độ rơi) (Đối chứng Davis Rain II) (Tính toán KE & USLE)
- So sánh thực nghiệm giữa phương pháp $k_{075}$ và phương pháp $M$:
Tiến hành thả các hạt mẫu chuẩn có kích thước biết trước qua buồng đo. Kết quả cho thấy phương pháp $k_{075}$ đề xuất duy trì sai số đường kính trung bình nhỏ hơn rõ rệt so với phương pháp đo theo độ sâu điều chế $M$ của nghiên cứu gốc, đặc biệt là ở dải hạt có kích thước nhỏ dưới 2 mm nơi tín hiệu $M$ bị méo dạng đáy xung.
- Kịch bản thực nghiệm 1 - Viên bi sắt chuẩn:
Sử dụng tập hợp bi sắt có đường kính xác định từ 1 mm đến 10 mm. Kết quả đo kiểm tra lặp lại nhiều lần khẳng định thiết bị đo có sai số tương đối ổn định trong ngưỡng cho phép đối với toàn bộ các cỡ hạt chuẩn.
- Kịch bản thực nghiệm 2 - Giọt chất lỏng từ ống nhỏ giọt y tế:
Kiểm tra với các đầu nhỏ giọt có kích thước kim khác nhau trong hai điều kiện: tốc độ nhỏ giọt cố định và tốc độ nhỏ giọt thay đổi; thả rơi từng giọt đơn lẻ và thả rơi đồng thời nhiều đầu giọt. Dữ liệu thực nghiệm xác nhận thiết bị nhận diện chính xác các xung quang điện liên tiếp mà không bị chồng lấn tín hiệu.
- Kịch bản thực nghiệm 3 - Mô hình mưa nhân tạo đối chứng:
Thiết bị của luận án được lắp đặt song song cùng thiết bị đo mưa chao lật thương mại The Rain Collector II (Davis Instruments - Mỹ) dưới cùng một hệ thống phun mưa nhân tạo. Dữ liệu lượng mưa tích lũy trung bình và cường độ mưa đo từ thiết bị quang học tự chế tạo có độ sai lệch tuyệt đối rất nhỏ so với thiết bị thương mại của Mỹ.
- Thử nghiệm thực tế tại Trạm Khí tượng Hà Đông và đề xuất ứng dụng:
Thiết bị được lắp đặt vận hành thực tế ngoài trời tại Trạm Khí tượng Hà Đông. Luận án trích xuất bộ dữ liệu chi tiết của một trận mưa điển hình:
- Xây dựng biểu đồ phổ phân bố kích thước hạt mưa (DSD) và mối quan hệ giữa kích thước hạt với xác suất xuất hiện;
- Thiết lập bảng quan hệ giữa cường độ mưa $I$ và tổng động năng của các hạt mưa rơi $KE_{time}(I)$ trong từng khoảng thời gian quan sát;
- Đề xuất quy trình ứng dụng trực tiếp thông số động năng $KE$ thu được từ thiết bị vào các mô hình dự báo mất đất do xói mòn như Phương trình mất đất tổng quát (USLE) và Phương trình mất đất hiệu chỉnh (RUSLE) tại Việt Nam.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học và lý luận
- Thiết lập thành công mô hình toán học giải tích mô tả dạng xung hai chồi đảo ngược và đề xuất tham số tỷ lệ xung $k_{075}$ (tại mức 75% biên độ xung lớn nhất) dùng để tính toán kích thước hạt mưa. Giải pháp này khắc phục được các hạn chế cơ bản của phương pháp đo theo độ sâu điều chế $M$ do D. Kiesewetter và V. I. Malyugin đề xuất, giúp phép đo không bị sai lệch khi hai đỉnh chồi xung mất cân bằng hoặc đáy xung bị bão hòa.
- Xây dựng thuật toán phân loại dạng xung quang điện (phân biệt xung hai chồi riêng biệt, xung hai chồi biến dạng và xung chập đơn) kết hợp công thức nội suy vận tốc và kích thước, mở rộng khả năng xử lý đối với các dòng hạt chuyển động có dải kích thước rộng.
Đóng góp mới về mặt công nghệ và kỹ thuật thực tiễn
- Làm chủ và cải tiến kết cấu phần cứng của thiết bị đo mưa hai dải sáng:
- Thay thế nguồn phát bước sóng 850 nm bằng nguồn LED thích hợp dải khả kiến kết hợp lò xo vi chỉnh quang học, giải quyết dứt điểm hiện tượng nhiễu đốm và đơn giản hóa thao tác cân chỉnh ngoài hiện trường;
- Thiết kế và chế tạo hoàn chỉnh các khối phần cứng chuyên dụng (khối phát quang, khối thu quang TIA, khối xử lý trung tâm vi điều khiển) cho phép thiết bị hoạt động độc lập theo thời gian thực mà không cần kết nối máy tính cồng kềnh;
- Xây dựng mô-đun phần mềm nhúng hoàn chỉnh tính toán trực tiếp các thông số vi vật lý của hạt mưa (đường kính $D$, vận tốc $v$) và thông số trận mưa tích lũy ($R, I, KE$).
- Chế tạo thành công thiết bị đo các thông số mưa thực tế, kiểm định hoạt động tin cậy ngoài trời tại Trạm Khí tượng Hà Đông và bước đầu đề xuất công thức tính toán động năng mưa ($KE$) ứng dụng trong đánh giá xói mòn đất đai tại Việt Nam.
Kiến nghị và đề xuất giải pháp
- Tích hợp thiết bị đo quang học chế tạo được vào mạng lưới trạm quan trắc khí tượng tự động quốc gia để cung cấp thêm thông tin phân bố hạt mưa phục vụ dự báo định lượng mưa bằng radar thời tiết.
- Sử dụng trực tiếp giá trị động năng mưa thực đo từ thiết bị để thay thế các hệ số thực nghiệm tính gián tiếp trong phương trình tính toán xói mòn đất RUSLE tại các khu vực canh tác nông nghiệp vùng núi có nguy cơ sạt lở cao.
Hạn chế và hướng nghiên cứu tiếp
Các hạn chế được ghi nhận trong luận án
- Giới hạn dải đo: Luận án tập trung khảo sát dải kích thước hạt phổ biến từ 0,5 mm đến 6 mm; các hạt mưa phùn kích thước nhỏ hơn 0,5 mm và các hạt dị thường lớn hơn 6 mm chưa được đánh giá chuyên sâu.
- Nhiễu môi trường tự nhiên: Thiết kế phần cứng và giải thuật hiện tại chưa tính đến các tình huống nhiễu môi trường ngẫu nhiên phát sinh ngoài thực địa như côn trùng bay cắt ngang vùng lấy mẫu hoặc ảnh hưởng của sương mù và hơi ẩm ngưng tụ trên bề mặt thấu kính.
- Kiểm chuẩn thông số vận tốc: Tác giả thừa nhận chưa thể đánh giá định lượng độ chính xác của phép đo vận tốc rơi do điều kiện phòng thí nghiệm chưa trang bị được thiết bị đối chứng vận tốc hạt chuyên dụng đạt chuẩn đo lường cấp cao.
Hướng nghiên cứu tiếp theo
- Hoàn thiện thuật toán xử lý tín hiệu thông minh trên vi điều khiển nhằm nhận diện và loại bỏ tự động các dạng xung dị thường do côn trùng hoặc bụi bẩn gây ra.
- Trang bị thêm cơ cấu sấy nhiệt bảo vệ bề mặt quang học tại các cửa sổ thu phát để triệt tiêu hiện tượng đọng sương khi làm việc lâu dài trong môi trường độ ẩm cao.
- Phối hợp với các phòng thí nghiệm đo lường chuyên chuẩn để kiểm định độc lập thông số vận tốc hạt; mở rộng thử nghiệm quan trắc liên tục qua nhiều mùa mưa nhằm thu thập ngân hàng dữ liệu phân bố hạt mưa phục vụ khí quyển học.
Giá trị tham khảo
Bản luận án là tài liệu kỹ thuật có giá trị tham khảo ứng dụng rõ nét cho các đối tượng chuyên môn sau:
- Nghiên cứu sinh và học viên sau đại học ngành Kỹ thuật Điện tử, Đo lường và Điều khiển, Quang điện tử: Có thể tham khảo phương pháp luận thiết kế hệ thống cảm biến quang điện tử, kỹ thuật tiền xử lý khuếch đại dòng quang điện cực nhỏ, kỹ thuật lọc nhiễu tương tự và giải thuật phân tích hình học xung số hóa.
- Kỹ sư phát triển thiết bị quan trắc môi trường và khí quyển: Cung cấp tài liệu thiết kế phần cứng hoàn chỉnh (các bản vẽ mạch nguyên lý, sơ đồ mạch in PCB của khối phát quang, thu quang, CPU nhúng trong phần Phụ lục) cùng quy trình cân chỉnh cơ quang học bằng thực nghiệm.
- Cơ quan quản lý khí tượng thủy văn và viện nghiên cứu nông nghiệp - địa chất: Cung cấp cơ sở dữ liệu thực nghiệm về phân bố hạt mưa tại khu vực Hà Nội và giải pháp công nghệ nội địa hóa thiết bị đo mưa hiện đại nhằm thay thế máy đo mưa nhập ngoại đắt tiền; cung cấp phương pháp trích xuất chỉ số động năng hạt mưa phục vụ đánh giá nguy cơ xói mòn đất theo phương trình USLE/RUSLE.
Câu hỏi thường gặp
1. Tại sao luận án thay thế phương pháp tính đường kính theo độ sâu điều chế $M$ bằng tham số $k_{075}$?
Phương pháp độ sâu điều chế gốc tính toán theo biên độ cực trị $M = I_{max} - I_{min}$. Trong thực tế, các đỉnh cực đại của xung thường xuyên bị nhiễu biên độ tần số cao gây méo và hai đỉnh chồi xung hiếm khi bằng nhau tuyệt đối do sai lệch quang cơ. Hơn nữa, khi hạt mưa có đường kính nhỏ hơn khoảng cách giữa hai khe nhạy sáng, đáy xung bị triệt tiêu (bằng không), khiến giá trị $M$ mất tính quy luật. Tham số $k_{075}$ trích xuất tại ngưỡng 75% chiều cao sườn xung giải quyết triệt để vấn đề này vì bỏ qua hoàn toàn nhiễu tại đỉnh xung 100% và không phụ thuộc vào mức độ sâu của đáy xung.
2. Tại sao luận án lại sử dụng bi sắt làm vật mẫu chuẩn trong quá trình thực nghiệm ban đầu?
Hạt mưa khi rơi trong không khí có tiêu cự quang học rất ngắn và chuyển động với vận tốc cao nên đáp ứng quang học của nó tương đương với một vật chắn sáng (vật đen tuyệt đối). Bi sắt hình cầu có kích thước cố định xác định trước (được kiểm tra chính xác bằng panme cơ khí), không bị biến dạng hình thái hoặc bay hơi như giọt nước tự do, giúp loại bỏ các biến số bất định để tập trung chuẩn hóa và kiểm chứng chính xác biểu thức toán học của xung quang điện.
3. Việc thay thế nguồn phát quang diode laser (LD) bằng đèn LED mang lại ưu thế kỹ thuật gì?
Nguồn phát laser diode (LD) có độ kết hợp quang học (coherence) quá cao, dẫn đến hiện tượng giao thoa tạo ra các đốm nhiễu quang học (speckle noise) trên bề mặt cảm biến thu, gây méo dạng sóng xung quang điện. Đèn LED có độ kết hợp thấp, tạo ra trường sáng đồng đều hơn, giúp biên dạng xung thu được mịn và ổn định. Ngoài ra, việc dùng LED có bước sóng thích hợp trong dải khả kiến giúp kỹ sư quan sát và vi chỉnh tiêu cự hệ thấu kính trực tiếp bằng mắt thường mà không cần thiết bị hiển thị hồng ngoại chuyên dụng.
4. Thiết bị của luận án tính toán động năng hạt mưa ($KE$) để phục vụ mục đích gì trong thực tiễn?
Thiết bị đo đồng thời được đường kính tương đương $D$ (từ đó suy ra thể tích và khối lượng giọt nước) và vận tốc rơi $v$ của từng hạt mưa. Từ hai thông số này, vi xử lý tính trực tiếp động năng của hạt $KE = \frac{1}{2} m v^2$ và tích lũy thành tổng động năng trận mưa theo thời gian $KE_{time}(I)$. Đây là thông số đầu vào có tính quyết định để tính toán sức công phá của hạt mưa làm phá vỡ kết cấu đất mặt, phục vụ trực tiếp việc dự báo xói mòn, rửa trôi đất trong canh tác nông nghiệp theo các mô hình USLE và RUSLE tại Việt Nam.
Kết luận
Bản luận án tiến sĩ của tác giả Lai Thị Vân Quyên đã giải quyết trọn vẹn bài toán nâng cao độ chính xác của phương pháp đo kích thước và vận tốc hạt mưa sử dụng hai xung quang điện kế thừa từ mô hình của D. Kiesewetter và V. I. Malyugin. Tác giả đã đề xuất thành công mô hình toán học dựa trên tham số tỷ lệ xung $k_{075}$, thay thế hoàn toàn phương pháp độ sâu điều chế truyền thống vốn nhạy cảm với nhiễu biên độ và hạn chế ở dải hạt nhỏ. Không chỉ dừng lại ở các đóng góp lý thuyết, công trình đã hoàn thiện công nghệ phần cứng thông qua việc cải tiến nguồn phát LED, tối ưu hóa cơ cấu định tâm trục quang bằng lò xo và chế tạo thành công thiết bị đo nhúng vi xử lý hoạt động độc lập theo thời gian thực. Các kết quả thử nghiệm đa dạng từ bi sắt chuẩn, giọt lỏng nhân tạo đến thực nghiệm hiện trường tại Trạm Khí tượng Hà Đông đã chứng minh tính khả thi, độ tin cậy và tiềm năng ứng dụng thiết thực của thiết bị trong mạng lưới quan trắc khí hậu và dự báo xói mòn đất đai.
