Kết luận chươn gI

CHƯƠNG I TỔNG QUAN CÁC PHƯƠNG PHÁP ĐO KÍCH THƯỚC HẠT MƯA

1.4. Kết luận chươn gI

Chương I của luận án đã trình bày tổng quan các phương pháp đo kích thước hạt mưa trong đó đi sâu về phân tích các phương pháp đo bằng hiệu ứng quang học. Từ những ưu nhược điểm của chúng, luận án lựa chọn phương pháp đo kích thước hạt mưa bằng hai dải sáng của D. V. Kiesewetter và V. I. Malyugin làm cơ sở để nghiên cứu nhằm nâng cao độ chính xác của kết quả đo.

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

Dựa trên những nghiên cứu tổng quan về các phương pháp đo kích thước hạt mưa ở chương I, bài toán và nội dung nghiên cứu của Luận án đã được xác định ở mục 1.3. Trong chương II, tác giả trình bày kỹ mơ hình (thiết bị) đo ban đầu của các giáo sư Nga xây dựng, chỉ rõ những hạn chế rồi sau đó đưa ra các đề xuất của tác giả để nâng cao độ chính xác của thiết bị đo. Cụ thể:

- Đề xuất biểu thức tốn học tính tốn đường kính, vận tốc hạt mưa dựa trên những phân tích về xung quang điện thu được bằng thực nghiệm.

- Đề xuất hồn thiện cơng nghệ, xây dựng hoàn thiện phần cứng thiết bị cho phép hoạt động trong môi trường thực tế.

Các kết quả của chương II đã được công bố trong các bài báo [1, 2, 3]

2.1. Mơ hình đo và những hạn chế

Năm 2016 Viện nghiên cứu Điện tử, Tin học, Tự động hóa có hợp tác với Trường Đại học Bách khoa Saint Peterburg (Nga) thực hiện đề tài hợp tác nghiên cứu để thiết kế, chế tạo thiết bị đo kích thước và vận tốc hạt mưa ứng dụng trong dự báo mưa lũ phục vụ cơng tác phịng chống biến đổi khí hậu, do đó chúng tơi đã được các Giáo sư Nga bàn giao cho mơ hình thiết bị đo và phần mềm xử lý trên máy tính để cùng nghiên cứu hồn thiện thiết kế, chế tạo được thiết bị có thể ứng dụng thực tế. Sau khi nhận được mơ hình, tác giả và nhóm nghiên cứu đã tìm hiểu, phân tích, thử nghiệm và tìm ra các hạn chế cần hồn thiện để có thể chế tạo được thiết bị ứng dụng trong thực tế.

Về khoa học, mơ hình gốc của các giáo sư Nga tính tốn đường kính hạt mưa dựa vào độ sâu điều chế M

trên mơ hình thử nghiệm. Như đã chỉ ra ở chương I, vấn đề xử lý thông qua M bị phụ thuộc vào biên độ cực đại và cực tiểu của xung quang điện thu được. Khi đó, vấn đề nhiễu ở các đỉnh xung, vấn đề điểm cực của các chồi xung không bằng nhau sẽ gây ra các sai số không mong muốn. Thêm nữa, với các hạt nhỏ hơn khoảng cách hai khe, điểm cực tiểu của khoảng trũng giữa hai chồi xung bằng khơng thì phương pháp tính này khơng dùng được nữa.

1– Nguồn sáng

2 – Khối cấp điện cho nguồn sáng 3 – Thấu kính tạo chùm song song 4 – Màn khe nhạy sáng

Fl: tiêu cự của thấu kính 3 Vd: vận tốc của hạt rơi

5 – Thấu kính hội tụ chùm song song 6 – Diode quang

7 – Khối khuếch đại và lọc 8 – Khối ADC và xử lý LS: độ rộng khoảng đo

L: khoảng cách từ hạt rơi đến màn 4

Hình 2.1. Sơ đồ kết cấu hệ đo do D. V. Kiesewetter and V. I. Malyugin đề xuất

Về công nghệ, trong mô hình đo D. V. Kiesewetter and V. I. Malyugin đề xuất (hình 2.1), nguồn sáng sử dụng

là LED nên đã giảm thiểu được nhiễu đốm trong chùm sáng do nguồn LD phát ra [62, 63, 64]. Tuy nhiên bước sóng ánh sáng trung tâm của nguồn sáng sử dụng là 850nm nằm trong dải hồng ngoại khiến cho việc hiệu chỉnh quang gặp nhiều khó khăn.

Bên cạnh đó, việc điều chỉnh nguồn sáng và điốt quang vào đúng vị trí tiêu cự của các thấu kính tương ứng là điều chỉnh cứng nên sẽ khó khăn khi cần hiệu chỉnh trục quang trong điều kiện thực tế sử dụng.

Theo các cơng trình của Illingworth and Stevens (1987) [30], Grossklaus và cộng sự (1998) [26] chỉ ra rằng gió khơng làm thay đổi đáng kể thành phần vận tốc dọc của hạt rơi xuống khi sử dụng dải sáng có hình trụ. Trong phương pháp đo nghiên cứu, dải sáng có dạng hình trụ và khe nhạy sáng bố trí song song với bề mặt của trái đất nên cũng giảm được ảnh hưởng của gió. Do đó, phổ của tín hiệu coi như khơng phụ thuộc vào tốc độ gió. Trong mơ hình đo (hình 2.1), bên khối phát quang có

dng

uồ

phần tạo dải sáng hình trụ từ nguồn sáng thơng qua thấu kính 3 nên tác giả sẽ kế thừa phần tạo chùm phát sáng.

Như vậy, trong chương II, tác giả sẽ trình bày các đề xuất nâng cao hiệu quả đo sản phẩm của cơng trình [14,16] cả về khoa học và công nghệ.

Về khoa học, thực hiện theo các nội dung sau:

-Đề xuất biểu thức tính tốn đường kính, vận tốc tương đương của hạt mưa dựa trên các xung quang điện thu được

-Đề xuất thuật tốn tính tốn đường kính và vận tốc hạt mưa.

Về cơng nghệ, thực hiện ở các mục sau:

-Thay đổi nguồn sáng để giảm thiểu nhiễu quang.

-Thay đổi cơ cấu điều chỉnh trục quang tạo điều kiện dễ dàng hơn cho điều chỉnh

-Hoàn thiện phần cứng xử lý dữ liệu hiện đang sử dụng máy tính để tiến tới thiết kế, chế tạo được thiết bị có thể hoạt động ngồi thực địa.

2.2. Đề xuất về khoa học

2.2.1. Cơ sở khoa học xây dựng biểu thức, thuật toán xử lý số liệu

2.2.1.1. Phương pháp thực nghiệm để trích xuất các xung quang điện

Chương I đã đưa ra quan điểm, kích thước hạt khi đề cập trong nghiên cứu này là kích thước tương đương của hạt mưa với hạt hình cầu có cùng đáp ứng quang điện. Bởi thế, có thể coi hạt mưa như là một vật hình cầu. Hạt mưa có tiêu cự ngắn nên có thể coi nó như một vật đen tuyệt đối. Bởi thế, có thể sử dụng viên bi sắt có kích thước biết trước để tiến hành thực nghiệm nhằm xác định các điểm đặc trưng của xung quang điện thu được rồi từ đó đưa ra các đề xuất.

Sử dụng mơ hình thực nghiệm (hình 2.2), lần lượt thả các viên bi (hạt) với kích thước khác nhau biết trước qua khoảng đo. Khi hạt rơi cắt ngang qua chùm sáng song song do nguồn sáng LED 1 và thấu kính hội tụ 3 tạo ra làm suy giảm cường độ ánh sáng nhận được trên cảm biến quang. Điều này khiến xung quang điện thu được trên photodiode 6 thay đổi theo trong suốt quá trình chuyển động của hạt qua khoảng đo có độ rộng LS. Xung quang điện sau khi được đưa qua các khối khuếch đại, lọc nhiễu 7 được đo trên oscillo để tiến hành phân

1– Nguồn sáng

2 – Khối cấp điện cho nguồn sáng 3 – Thấu kính tạo chùm song song 4 – Màn khe nhạy sáng

Fl: tiêu cự của thấu kính 3 LS: độ rộng khoảng đo

5 – Thấu kính hội tụ chùm song song 6 – Cảm biến quang

7 – Khối khuếch đại và lọc

w1, w2: lần lượt là độ rộng khe nhạy sáng 1, 2 g: khoảng cách giữa hai khe nhạy sáng

Hình 2.2. Mơ hình thực nghiệm xác định xung quang điện

Với các hạt có kích thước khác nhau nhận thấy các xung quang điện thu được cũng có những dạng khác nhau.

2.2.1.2. Phân tích các xung quang điện

 Dạng xung quang điện có hai chồi

g w1

Umax1 l1 1 Umax2 l2 l3 2 3 4 t1 t2 t3 g 2 1 4 3 w1 w2 Umin

Xung quang điện thu được trên cảm biến quang khi có hạt đi qua khoảng đo của mơ hình thực nghiệm. Giá trị nhỏ nhất của xung quang điện đạt được ở điểm hạt che giữa tâm khe nhạy sáng (điểm 2 và 4 trong hình 2.3). Giá trị lớn nhất của xung quang điện đạt được ở điểm hạt nằm ở trung tâm khoảng cách giữa hai khe nhạy sáng (điểm 3 trong hình 2.3). Biên độ của chồi xung quang điện thể hiện lượng ánh sáng bị suy hao khi hạt đi qua khoảng đo. Khoảng cách giữa các sườn xung thể hiện thời gian hạt đi qua một quãng đường nào đó nằm trong khoảng đo.

Khơng làm mất tính tổng qt, xét xung quang điện được đảo lại (hình 2.4) so với xung quang điện mơ tả ở hình 2.3.

Hình 2.4. Hình dạng xung hai chồi đảo ngược trên cảm biến quang khi hạt rơi

Xét quãng đường hạt rơi từ điểm 1  2  3  4 với điểm 1, 2 cách đều tâm khe nhạy sáng có độ rộng w1, điểm 3,4 cách đều tâm khe nhạy sáng có độ rộng w2, khoảng cách giữa hai khe là g như mơ tả trên hình 2.4.

Gọi quãng đường hạt rơi từ điểm 1 đến điểm 2 mô tả trên hai khe nhạy sáng là l12 và tương ứng với chồi xung từ điểm 1 đến điểm 2

1 1 1 1 1 1 v v    v 12

Quãng đường hạt rơi từ điểm 2 đến điểm 3 mô tả trên hai khe nhạy sáng là l23 và tương ứng với chồi xung từ điểm 2 đến điểm 3

Quãng đường hạt rơi từ điểm 3 đến điểm 4 mô tả trên hai khe nhạy sáng là l34 và tương ứng với chồi xung từ điểm 3 đến điểm 4

Quãng đường l12, l23, l34 không đổi.

- Xét với cùng một cỡ hạt đường kính D đi qua khoảng đo với các vận tốc

khác nhau.

Độ rộng khoảng đo rất nhỏ so với chiều cao mà hạt rơi nên có thể coi vận tốc hạt đi qua khoảng đo là đều. Gọi các vận tốc đó lần lượt là

khi hạt đi qua dải sáng v ,v, và v ,v, không đổi Hạt rơi với vận tốc v1 , thời gian hạt đi hết quãng đường l12, l23, l34 lần lượt là t12, t23, t34

Hạt rơi với vận tốc

v,

, thời gian hạt đi hết quãng đường l12, l23, l34 lần lượt là

t , ,t , ,t ,

12 23 34

Khi đó, thời gian hạt đi trên những quãng đường l12, l23, l34 được tính như sau: Với vận tốc v1 Với vận tốc v, t12 t23 t34  l12 v1  l23 v1  l34 v1 , l12 12 , 1 , l23 23 , 1 , l34 34 , 1 Khơng làm mất tính tổng qt, giả sử v,  1 .v Khi đó : ,  n. l12 v 1 n 1  n.t12 (2.1) t t t t

23 34 ,  n. l23 v1 ,  n. l34 v1  n.t23  n.t34 ' t , (2.2) (2.3) Ứng với từng vận tốc, xét các tỷ số: k  23 t ,  t ,  t' (2.4) k 

Thay (2.1), (2.2), (2.3) vào (2.4) thu được

12 23 34 t23 t12  t23  t34 (2.5) k' n.t 23 n.t12  n.t23  n.t34 (2.6)

Rút gọn n ở cả từ số là mẫu số thu được

 k'  t23

t12  t23  t34

(2.7)

Từ (2.5) và (2.7), có thể thấy

đổi khi vận tốc hạt thay đổi. k' 

k tức là với cùng một cỡ hạt, tỷ số k khơng

- Xét các hạt có đường kính khác nhau bay cùng vận tốc qua khoảng đo Khi đường kính D của hạt càng lớn thì lượng áng sáng mất đi càng nhiều dẫn đến giá trị lớn nhất Umax1 và Umax2 của các chồi xung càng lớn, biên độ của các điểm xung càng tăng khiến độ rộng các chồi xung tăng lên. Mà khoảng cách, độ rộng của hai khe nhạy sáng là không đổi nên khoảng trũng giữa hai chồi xung sẽ bị thu hẹp lại. Điều này có nghĩa rằng các khoảng thời gian t12, t34 tăng lên còn t23 giảm khi D tăng. Theo biểu thức (2.5) có: k  t23 t12  t23  t34 k  1 t12  t34  1 t

chia cả tử số và mẫu số cho t23 ta được

t t



Khi D tăng thì:   (t12+ t34) tăng t23 giảm   t12  t34  1 tăng    k  t23  1 giảm t12  t34  1 t23

Mà theo kết luận phần trên, ở cùng cỡ hạt, khi vận tốc thay đổi thì tỷ số k khơng thay đổi. Do đó: khi đường kính D của hạt tăng thì tỷ số k giảm.

Như vậy với dạng xung hai chồi có thể dựa vào tỷ số k ở biểu thức 2.5 để tính đường kính hạt. Đây là một biến số không thay đổi khi vận tốc của một hạt rơi biến đổi và biến số này giảm khi đường kính hạt tăng.

Khi sử dụng tỷ số k để tính kích thước hạt thì chỉ có thể tính được khi có khoảng trũng giữa hai chồi xung tức có dạng xung hai chồi. Mà dạng xung hai chồi chỉ có thể có được khi t23>0. Điều này chỉ xảy ra khi D < w1+w2+g. Trong khoảng đường kính này:

Khi D  g, điểm cực tiểu bắt đầu có xu hướng tăng dần lên tức Umin tăng lên tiến dần đến Umax khi D tăng.

Khi D < g, điểm cực tiểu bằng khơng và có xu hướng đi ngang một đoạn trước khi đi lên để tạo thành chồi xung tiếp theo, độ rộng khoảng đi ngang của điểm cực tiểu tỉ lệ nghịch với đường kính D. (xem phục lục 1 hình các xung quang điện khi thả các hạt đường kính khác nhau qua khoảng đo của thiết bị thiết kế ở mục 2.3.4).

Do đó, cần lựa chọn điểm xác định tỷ số k sao cho càng gần điểm cực đại càng tốt nhưng lại không rơi vào khu vực có khả năng xảy ra nhiễu ở các đỉnh xung. Cụ thể tác giả sẽ trình bày ở mục 2.2.2.

Khi D  w1+w2+g thì t23 = 0, lúc này dạng xung trở thành xung một chồi.

 Dạng xung quang điện có một chồi xung

Xung quang điện có một chồi xung xảy ra khi hạt rơi qua khoảng đo có đường

tăng bằng và vượt điểm cực đại (hình 2.5). Khi đó, điểm cực tiểu của dạng xung hai chồi trở thành điểm cực đại của dạng xung một chồi. Thử nghiệm với tập các hạt có đường kính ở dải quan sát, biên độ cực đại của xung, diện tích bao bởi xung thay đổi đồng biến với đường kính hạt. Hình 2.5 là ví dụ xung quang điện sau tiền xử lý thu được khi đo trên oscillo với đường kính hạt trong dải kích thước này khi thả hạt mẫu qua phần cứng thiết kế ở mục 2.3.4. Với đường kính hạt D = 8mm thì biên độ xung cực đại Umax = 2,8V; với đường kính hạt D = 10mm thì biên độ xung cực đại Umax = 4,28V; với đường kính hạt D = 12mm thì biên độ xung cực đại Umax = 4,28V và bằng đầu do rơi vào trạng thái bão hịa của mạch khuếch đại. Như vậy có thể dựa vào biên độ cực đại của xung, độ bằng đầu của xung để nội suy ra đường kính của hạt khi kích thước hạt ở dải đo lớn hơn hoặc bằng tổng khoảng cách giữa các khe và chiều rộng của hai khe.

(a) (b)

(c) a. Xung quang điện của hạt có đường kính D=10mm

b. Xung quang điện của hạt có đường kính D = 8mm Xung quang điện của hạt có đường kính D=12mm

2.2.2. Đề xuất biểu thức, thuật tốn xử lý tính kích thước và tốc độ hạt2.2.2.1. Biểu thức tính kích thước hạt 2.2.2.1. Biểu thức tính kích thước hạt

+ Với các hạt tạo xung quang điện dạng hai chồi

Đường kính của các hạt và tốc độ của hạt có thể xác định chính xác nhất bằng phương pháp so sánh xung nhận được với hình dạng xung tham chiếu (hình dạng xung mẫu từ các hạt có đường kính biết trước). Phương pháp này địi hỏi phải có một tập mẫu các xung tương ứng với các kích thước hạt khác nhau. Ứng với mỗi thiết bị đo có cấu trúc tương tự nhau nhưng các tham số của các thành phần quang, điện của các thiết bị đó lại có sự sai khác nên tập mẫu này cần phải xây dựng riêng với từng thiết bị. Trong trường hợp xây dựng được tập mẫu và lưu được vào bộ nhớ thì việc so sánh xung nhận được với xung tham chiếu mẫu sẽ cho kết quả chính xác hơn tuy nhiên việc này rất tốn thời gian khi điểm so sánh nhiều, nên không thể cho kết quả theo thời gian thực. Từ những phân tích về dạng xung hai chồi ở mục 2.2.1, tác giả đề xuất chọn biến số k ở biểu thức 2.5 để tính đường kính hạt. Đây là một biến số không thay đổi khi vận tốc của một hạt rơi biến đổi và biến số này nghịch