Tổng quan luận án

Luận án tiến sĩ vật lý với đề tài "Nghiên cứu xây dựng hệ đo thông số quang sử dụng công nghệ xử lý tín hiệu số và ứng dụng", chuyên ngành Vật lý Vô tuyến và Điện tử (mã số: 62 44 03 01), được nghiên cứu sinh Nguyễn Tuấn Anh thực hiện tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội năm 2011, dưới sự hướng dẫn khoa học của Bạch Gia Dương và Vũ Anh Phi.

Về tính cấp thiết và khoảng trống nghiên cứu, sự bùng phát của các loài vi tảo độc, hại trong các thủy vực ven bờ biển (gây ra hiện tượng thủy triều đỏ hay sự nở hoa của nước) là nguyên nhân và hậu quả trực tiếp của tình trạng ô nhiễm môi trường thủy vực. Khi môi trường chuyển từ nghèo dưỡng sang giàu dưỡng, mật độ tảo độc có thể tăng nhanh từ vài chục, vài trăm tế bào/ml lên đến hàng triệu tế bào/ml chỉ trong vài ngày, sau đó phân hủy và giải phóng độc tố gây tổn hại cho hệ sinh thái, thủy hải sản và sức khỏe con người. Theo khuyến cáo của Tổ chức Y tế Thế giới (WHO), khi mật độ tảo lam chạm ngưỡng 1000 tế bào/ml, nguồn nước đã bị ô nhiễm nghiêm trọng. Tuy nhiên, việc phát hiện và cảnh báo sớm ở mật độ thấp dưới 1000 tế bào/ml ngoài hiện trường hiện gặp nhiều khó khăn. Các phương pháp truyền thống trong phòng thí nghiệm đòi hỏi thời gian phân tích từ vài ngày đến vài tuần; phương pháp viễn thám vệ tinh kết hợp đo đạc mặt đất có chi phí cao và dễ bị nhiễu do bức xạ khí quyển, đáy nước; trong khi các phương pháp quang phổ huỳnh quang hiện thời chỉ phân loại được khoảng 5 nhóm tảo do trùng lặp sắc tố. Khoảng trống nghiên cứu được tác giả xác định là việc thiếu hụt một hệ thống đo lưu động có khả năng xác định trực tiếp, đồng thời và nhanh chóng bộ ba thông số quang cơ bản của các phần tử lơ lửng trong môi trường nước ở mật độ thấp gồm: hệ số hấp thụ ($\mu_a$), hệ số tán xạ ($\mu_s$) và hệ số bất đối xứng ($g$).

Mục tiêu nghiên cứu của luận án được xác định gồm:

  1. Xây dựng hệ đo thông số quang của thủy vực vùng ven bờ biển Việt Nam có mật độ tảo thấp trên cơ sở nghiên cứu các mô hình phù hợp, dựa trên công nghệ tích hợp quang điện tử và kỹ thuật xử lý tín hiệu. Hệ đo hoạt động theo hai cấu hình: cấu hình đo nhanh và cấu hình đo chính xác ba thông số quang ($\mu_a$, $\mu_s$, $g$).
  2. Trên cơ sở hệ đo đã xây dựng, chế tạo thiết bị đo lưu động cho phép đo nhanh và đồng thời ba thông số quang của mẫu thủy vực có mật độ tảo thấp.

Nội dung nghiên cứu chính bao gồm:

  1. Nghiên cứu, mô tả đối tượng dưới góc độ vật lý tương tác photon - sinh học ($\mu_a$, $\mu_s$, $g$).
  2. Lựa chọn và phân tích chi tiết mô hình Kubelka-Munk và mô phỏng Monte Carlo áp dụng cho thủy vực ven bờ; kết hợp kỹ thuật quả cầu tích phân đôi.
  3. Nghiên cứu kỹ thuật xử lý tín hiệu ánh sáng: chế tạo nguồn sáng chuẩn trực có điều khiển ổn định công suất laser diode, chế tạo quang hệ và hệ hai quả cầu tích phân.
  4. Phát triển ứng dụng xử lý tín hiệu số (DSP) trên TMS320C6713 sử dụng kiến trúc song song và kỹ thuật đường ống, thiết kế chế tạo mạch ADC chuyên dụng tích hợp qua giao tiếp EMIF.
  5. Thực hiện đo đạc thực nghiệm trên các mẫu sữa tiêu chuẩn và ba mẫu tảo độc được phân lập, nuôi cấy ở các mật độ khác nhau.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu vật lý: Bộ ba thông số quang của môi trường khuếch tán ($\mu_a$, $\mu_s$, $g$) và các thành phần quang thông phản xạ ngược ($R_d$), tán xạ xuôi ($T_d$), truyền qua chuẩn trực ($T_c$).
  • Đối tượng thực nghiệm: Mẫu dung dịch sữa chuẩn có tỷ lệ chất béo 4% được pha loãng ở các nồng độ khác nhau (dùng làm môi trường chuẩn quang học mô phỏng hạt cầu lơ lửng) và 03 mẫu vi tảo độc được phân lập từ vùng biển Đồ Sơn gồm Prorocentrum rhathymum, Pseudo-nitzschia và Alexandrium minutum.
  • Phạm vi không gian và thời gian: Mẫu nước và vi tảo được khảo sát thuộc vùng thủy vực ven bờ biển Việt Nam (tập trung tại vùng biển Đồ Sơn, Hải Phòng); luận án hoàn thành năm 2011.

Tổng quan tài liệu và vị trí của luận án

Trong Chương 1, tác giả đã hệ thống hóa các hướng tiếp cận quan trắc thủy vực và xác định thông số quang học trên thế giới và tại Việt Nam:

  1. Nhóm phương pháp phân tích sinh học và hóa học truyền thống trong phòng thí nghiệm: Bao gồm phương pháp Utermöhl (sử dụng buồng lắng hình trụ kết hợp buồng phẳng để đếm dưới kính hiển vi), phương pháp buồng đếm Sedgewick Rafter (ô trung tâm chứa 100 $\text{mm}^3$, chia lưới vi thể để đếm mật độ), kỹ thuật nhuộm màu mô học (sử dụng thuốc nhuộm Carmin cho vách cellulose hoặc pectin, xanh iod cho vách thấm lignin/suberin), và phương pháp phân tích miễn dịch gắn enzyme ELISA (sử dụng kháng thể bắt giữ kháng nguyên đặc hiệu). Tác giả chỉ ra rằng các phương pháp này đòi hỏi 4 công đoạn phức tạp: thu thập mẫu tại hiện trường (dùng tàu, thuyền, lưới vớt, GPS), bảo quản mẫu trong bình chuyên dụng, vận chuyển về phòng thí nghiệm và lưu trữ mẫu. Quy trình này tiêu tốn nhiều nhân lực, thời gian kéo dài từ vài ngày đến vài tuần và ghi chép thủ công dễ phát sinh sai sót, không đáp ứng yêu cầu cảnh báo sớm hiện tượng tảo nở hoa.

  2. Nhóm phương pháp viễn thám vệ tinh kết hợp đo đạc thực địa: Sử dụng các đầu thu quang học độ phân giải phổ trung bình và thấp đặt trên vệ tinh như SeaWiFS, MODIS, MERIS để thu nhận phổ bức xạ phản xạ từ mặt biển $L_d = T_{\text{atm}}(L_w + L_s + L_b) + L_a + L_l$. Hệ số phản xạ bức xạ $R(\lambda)$ liên hệ với hệ số tán xạ ngược $b_b(\lambda)$ và hệ số hấp thụ tổng hợp $a(\lambda)$ qua công thức $R(\lambda) = f \cdot b_b(\lambda) / [a(\lambda) + b_b(\lambda)]$. Phương pháp này cho phép theo dõi trên diện tích đại dương rộng lớn nhưng chi phí vận hành rất cao, phụ thuộc vào điều kiện thời tiết, mây mù và đòi hỏi các thuật toán hiệu chỉnh khí quyển phức tạp để tách bức xạ của cột nước $L_w$ khỏi bức xạ khí quyển $L_a$ và bề mặt $L_s$.

  3. Nhóm phương pháp quang học đo huỳnh quang của tảo: Dựa trên đặc tính phát quang của Chlorophyll a (hấp thụ cực đại tại bước sóng 436 nm, phát huỳnh quang tại 685 nm), Chlorophyll b (hấp thụ tại 645 nm) và Chlorophyll c (hấp thụ tại 630 nm). Các công trình tiêu biểu gồm dự án quan trắc nguồn nước của Đại học Tổng hợp Umaine (Orono, Hoa Kỳ, khởi xướng từ tháng 3/2007) và thiết bị đo huỳnh quang hiện trường của hãng bbe Moldaenke (Cộng hòa Liên bang Đức) sử dụng 5 đèn LED kích thích tại các bước sóng đỉnh 450 nm, 525 nm, 570 nm, 590 nm và 610 nm. Tuy nhiên, do nhiều loài vi tảo có cùng nhóm sắc tố nên phương pháp huỳnh quang chỉ phân biệt được tối đa khoảng 5 nhóm tảo theo vùng phổ kích thích, không phản ánh trực tiếp các chỉ số tán xạ vật lý của tế bào.

Vị trí của luận án: Luận án lựa chọn giải quyết khoảng trống về công nghệ đo đạc trực tiếp bằng cách tích hợp quang học quả cầu tích phân đôi với bộ vi xử lý tín hiệu số tốc độ cao DSP TMS320C6713. Hệ đo khắc phục triệt để ảnh hưởng của các thành phần phản xạ khí quyển và đáy nước, đo trực tiếp tại mẫu dung dịch không qua quy trình nhuộm rửa phức tạp, và cung cấp giải pháp kép: đo nhanh bằng thuật toán giải tích Kubelka-Munk phục vụ khảo sát lưu động tại hiện trường, và đo chính xác bằng mô phỏng ngẫu nhiên Monte Carlo.

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

Luận án xây dựng cơ sở lý thuyết dựa trên các quy luật tương tác giữa photon và môi trường sinh học lơ lửng:

  1. Quá trình hấp thụ và tán xạ ánh sáng:
  • Hấp thụ được biểu diễn theo định luật Beer-Lambert: $I = I_0 e^{-\mu_a l} = I_0 e^{-\rho \sigma_a l}$, với $\mu_a$ là hệ số hấp thụ ($\text{mm}^{-1}$), $\rho$ là mật độ hạt và $\sigma_a$ là tiết diện hấp thụ. Quãng đường tự do trung bình trước khi photon bị hấp thụ là $1/\mu_a$.
  • Tán xạ được phân loại thành tán xạ đàn hồi (năng lượng photon bảo toàn, gồm tán xạ Rayleigh cho các hạt có hệ số kích thước $x = ka \ll 1$ và tán xạ Mie cho các hạt hình cầu kích thước tương đương hoặc lớn hơn bước sóng) và tán xạ không đàn hồi (huỳnh quang với thời gian sống $10^{-9}$ đến $10^{-6}$ giây, lân quang $10^{-6}$ đến $10^2$ giây, và tán xạ Raman với các vạch phổ Stokes và đối Stokes).
  • Hệ số bất đối xứng $g = \langle \cos\theta \rangle = \int_{-1}^{1} \mu p(\mu) d\mu$, nhận giá trị từ $-1$ (tán xạ ngược hoàn toàn) đến $+1$ (tán xạ xuôi hoàn toàn), với $g = 0$ là tán xạ đẳng hướng. Hàm pha tán xạ sử dụng trong mô hình là hàm Henyey-Greenstein ($P_{\text{HG}}$): $$P_{\text{HG}}(\theta) = \frac{1 - g^2}{4\pi (1 + g^2 - 2g\cos\theta)^{3/2}}$$
  • Hệ số suy giảm tổng thể $\mu_t = \mu_a + \mu_s$ và hệ số suy giảm rút gọn $\mu_{\text{tr}} = \mu_a + \mu_s(1 - g) = \mu_a + \mu_s'$.
  1. Phương trình truyền bức xạ và mô hình Kubelka-Munk (KM):
  • Phương trình truyền bức xạ tổng quát: $$\frac{dL(\vec{r}, \hat{s})}{ds} = -\mu_a L(\vec{r}, \hat{s}) - \mu_s L(\vec{r}, \hat{s}) + \mu_s \int_{4\pi} p(\hat{s}', \hat{s}) L(\vec{r}, \hat{s}') d\omega' + S(\vec{r}, \hat{s})$$
  • Mô hình Kubelka-Munk hai dòng bức xạ: Đơn giản hóa quá trình truyền bức xạ trong môi trường nhiều lớp mỏng song song thành hai dòng bức xạ ngược chiều nhau: dòng truyền xuôi $i_t(x)$ và dòng truyền ngược $i_r(x)$. Hệ phương trình vi phân tuyến tính vi mô: $$\frac{di_r(x)}{dx} = (K + S)i_r(x) - S i_t(x)$$ $$\frac{di_t(x)}{dx} = S i_r(x) - (K + S)i_t(x)$$ Trong đó $K$ và $S$ lần lượt là hệ số hấp thụ và tán xạ theo mô hình KM. Tác giả sử dụng phép biến đổi Laplace để giải hệ phương trình vi phân với các điều kiện biên phản xạ bề mặt ($r_0, r_h$), từ đó suy ra các biểu thức tường minh liên hệ giữa độ phản xạ khuếch tán $R_d$, độ truyền qua khuếch tán $T_d$ với $K$ và $S$, cũng như mối liên hệ giữa $(K, S)$ với $(\mu_a, \mu_s)$.
  1. Lý thuyết quả cầu tích phân:
  • Bề mặt phản xạ Lambertian phân bố đều năng lượng quang thông bên trong khoang cầu rỗng.
  • Hệ số nhân của quả cầu tích phân ($M$), biểu thị sự tăng cường độ dọi do phản xạ khuếch tán lặp lại nhiều lần: $$M = \frac{\rho_0}{1 - \rho_w(1 - \sum f_i) - \sum \rho_i f_i}$$ Với $\rho_w$ là độ phản xạ của thành trong quả cầu, $f_i$ là tỷ lệ diện tích các cửa mở so với tổng diện tích mặt cầu.
  • Hằng số thời gian đáp ứng của quả cầu tích phân: $\tau = -\frac{2 D_s}{3 c \ln\bar{\rho}}$, trong đó $D_s$ là đường kính quả cầu, $c$ là tốc độ ánh sáng, $\bar{\rho}$ là hệ số phản xạ trung bình.
  • Các công thức đo đạc thực nghiệm xác định quang thông:
    • Độ truyền qua khuếch tán: $T_d = A_1 / A_2$ (với $A_1, A_2$ lần lượt là tín hiệu đo khi có mẫu và không có mẫu tại cửa vào quả cầu thứ hai).
    • Độ phản xạ khuếch tán: $R_d = (X - Y) / (Z - Y)$ (với $X$ là tín hiệu khi đặt mẫu tại cửa phản xạ, $Y$ là cường độ ánh sáng nền khi tắt laser, $Z$ là tín hiệu đo với tấm phản xạ chuẩn 100%).
    • Độ truyền qua chuẩn trực: $T_c = I / I_0$.
  1. Phương pháp mô phỏng Monte Carlo (MC):
  • Mô phỏng bước di chuyển của photon dựa trên lấy mẫu biến ngẫu nhiên $\xi \in (0, 1]$ từ phân bố xác suất quãng đường tự do: $s = -\ln(\xi) / \mu_t$.
  • Góc lệch tán xạ $\theta$ được lấy mẫu theo hàm Henyey-Greenstein: $$\cos\theta = \frac{1}{2g} \left[ 1 + g^2 - \left( \frac{1 - g^2}{1 - g + 2g\xi} \right)^2 \right] \quad (\text{khi } g \neq 0)$$ $$\cos\theta = 2\xi - 1 \quad (\text{khi } g = 0)$$
  • Góc phương vị được lấy mẫu đẳng hướng: $\psi = 2\pi\xi$.
  • Kỹ thuật mô phỏng kết hợp quả cầu tích phân đôi ghi nhận lịch sử hàng triệu photon để tính lặp ngược các giá trị $(\mu_a, \mu_s, g)$ khớp với bộ ba thực nghiệm ($R_d, T_d, T_c$).
  1. Phương pháp xử lý tín hiệu và phần cứng:
  • Bộ xử lý tín hiệu số DSP TMS320C6713 (Texas Instruments) hoạt động ở xung nhịp cao, cấu trúc VLIW hỗ trợ xử lý song song và đường ống lệnh (pipeline). Lập trình kết hợp giữa ngôn ngữ C và Assembler tối ưu hóa thời gian tính toán.
  • Bộ biến đổi tương tự - số (ADC) chuyên dụng độ phân giải cao được thiết kế riêng, ghép nối với DSP qua cổng giao tiếp bộ nhớ ngoài EMIF và mạch logic CPLD.
  • Hệ quang gồm nguồn phát laser diode chip SLD102 (bước sóng đỏ) có mạch ổn định công suất quang học tự động (APC), đầu thu quang photodiode tích hợp khuếch đại (hybrid) OPT101, và hệ quang chuẩn trực tạo chùm tia song song đi qua buồng mẫu.

Nội dung chính theo từng chương

Chương 1: Tổng quan về phương pháp đo thông số quang

Chương này khảo sát thực trạng môi trường biển ven bờ Việt Nam với khoảng 70 loài vi tảo độc, hại đã được ghi nhận. Tác giả phân tích chi tiết quy trình quan trắc môi trường biển hiện hành theo TCVN 5943-1995 và tiêu chuẩn ngành của Bộ Thủy sản năm 2006, chỉ ra sự chậm trễ và tốn kém khi phải duy trì mạng lưới lấy mẫu thủ công tại các trạm quan trắc nội địa và hải đảo. Chương 1 điểm lại các phương pháp truyền thống (Utermöhl, buồng đếm Sedgewick Rafter, kỹ thuật nhuộm màu Carmin/pectin/iod, kỹ thuật miễn dịch ELISA) và các công nghệ hiện đại trên thế giới (viễn thám đa phổ vệ tinh MODIS/MERIS/SeaWiFS, thiết bị huỳnh quang bbe của Đức, dự án Đại học Umaine). Từ các phân tích kỹ thuật, tác giả đề xuất mô hình hệ đo mới kết hợp hai quả cầu tích phân đặt nối tiếp (quả cầu tích phân đôi) cùng ba đầu thu photodiode (PD1 đo tán xạ ngược $R_d$, PD2 đo tán xạ xuôi $T_d$, PD3 đo truyền qua chuẩn trực $T_c$), kết nối với khối xử lý tín hiệu số DSP và máy tính điều khiển PC.

Chương 2: Cơ sở lý thuyết xây dựng hệ đo thông số quang

Chương 2 trình bày có hệ thống bản chất vật lý của quá trình truyền bức xạ trong môi trường đục sinh học. Tác giả thiết lập hệ thống định nghĩa và biểu thức toán học cho các đại lượng quang học cơ bản: tiết diện hấp thụ $\sigma_a$, hệ số hấp thụ $\mu_a$, hệ số tán xạ $\mu_s$, hàm pha Henyey-Greenstein $P_{\text{HG}}(\theta)$, hệ số bất đối xứng $g$, hệ số suy giảm tổng thể $\mu_t$ và hệ số suy giảm rút gọn $\mu_{\text{tr}}$. Tác giả chứng minh phương trình truyền bức xạ Boltzmann và phân tích lý thuyết phản xạ khuếch tán Lambertian trong quả cầu tích phân. Chương này đi sâu vào việc mô tả cơ chế nhân quang thông trong quả cầu, tính toán hằng số thời gian trễ $\tau$, các công thức thực nghiệm xác định $R_d, T_d, T_c$, đồng thời thiết lập thuật toán lấy mẫu ngẫu nhiên Monte Carlo cho bước tự do $s$, góc tán xạ $\theta$ và góc phương vị $\psi$.

Chương 3: Áp dụng mô hình Kubelka-Munk và mô phỏng Monte Carlo phù hợp với kỹ thuật quả cầu tích phân đôi

Nội dung trọng tâm của Chương 3 là chuyển đổi các lý thuyết toán lý sang thuật toán tính toán cụ thể cho hệ quả cầu tích phân kép. Tác giả trình bày chi tiết lời giải hệ phương trình vi phân Kubelka-Munk bằng phép biến đổi Laplace, đưa ra nghiệm giải tích biểu diễn dòng bức xạ $i_t(x)$ và $i_r(x)$ theo các hàm hypebolic $\sinh(bSh)$ và $\cosh(bSh)$ với $a = (K + S)/S$ và $b = \sqrt{a^2 - 1}$. Tác giả chứng minh rằng đối với mẫu có độ dày quang học vô hạn, độ phản xạ $R_\infty$ chỉ phụ thuộc vào tỷ số $K/S$ theo công thức $(1 - R_\infty)^2 / (2 R_\infty) = K/S$. Tác giả thiết lập thuật toán chuyển đổi giữa các hệ số $(K, S)$ của mô hình KM với các hệ số quang học vật lý $(\mu_a, \mu_s)$. Đối với mô phỏng Monte Carlo, tác giả xây dựng lưu đồ thuật toán chi tiết mô phỏng đường đi của photon qua cuvet mẫu dạng phẳng đặt giữa hai quả cầu tích phân, có xét đến phản xạ và khúc xạ tại các mặt phân giới theo định luật Fresnel và góc khúc xạ Snell, cho phép tính lặp khớp dữ liệu để tìm ra bộ thông số $(\mu_a, \mu_s, g)$ với độ chính xác cao.

Chương 4: Xây dựng hệ đo, chế tạo thiết bị và đo mẫu

Chương 4 là chương thực nghiệm kỹ thuật, mô tả toàn bộ quá trình thiết kế, chế tạo phần cứng, phát triển phần mềm và thực hiện các phép đo kiểm chứng:

  1. Thiết kế và tối ưu hóa hệ thống xử lý tín hiệu số:
  • Sử dụng bo mạch DSK TMS320C6713. Tác giả tiến hành khảo sát thời gian tính toán của DSP đối với phép nhân hai ma trận $n \times n$ (Bảng 4.1), so sánh hiệu năng giữa việc có sử dụng và không sử dụng cấu trúc song song, giữa chương trình viết bằng ngôn ngữ C thuần túy và Assembler kết hợp cấu trúc song song và kỹ thuật đường ống lệnh.
  • Thiết kế mạch phần cứng ADC tốc độ cao, ghép nối đồng bộ với bus EMIF của DSP, tạo các tín hiệu điều khiển đọc thông qua CPLD.
  • Kiểm tra độ chính xác của hệ thống bằng các tín hiệu điện áp chuẩn ở đầu vào: Sai khác biên độ giữa tín hiệu chuẩn lối vào với giá trị hiển thị (Bảng 4.2) và sai khác về tần số (Bảng 4.3) đều nằm trong giới hạn sai số cho phép của thiết bị đo lường chuẩn.
  1. Chế tạo khối quang học và cơ khí:
  • Chế tạo hệ hai quả cầu tích phân bằng kim loại phủ lớp vật liệu phản xạ khuếch tán cao ở thành trong.
  • Nguồn sáng: Sử dụng chip laser bán dẫn SLD102 (Bảng 4.5) kết hợp hệ thấu kính chuẩn trực tạo tia song song và mạch điều khiển ổn định công suất laser.
  • Bộ thu tín hiệu quang: Sử dụng ba mạch tích hợp photodiode lai khuếch đại OPT101 (Bảng 4.6) có độ nhạy cao và độ tuyến tính đáp ứng tần số tốt.
  1. Chế tạo thiết bị hoàn chỉnh và giao diện phần mềm:
  • Thiết bị được đóng gói dạng hộp đo xách tay, tích hợp bộ nguồn, khối quang học quả cầu đôi, khối ADC-DSP và cổng kết nối máy tính.
  • Xây dựng phần mềm trên máy tính với giao diện cấu hình phép đo (chọn chế độ đo nhanh KM hoặc đo chính xác MC), cửa sổ hiển thị đồ thị và chức năng tự động lưu trữ dữ liệu đo.
  1. Kết quả thực nghiệm trên mẫu dung dịch sữa chuẩn và mẫu tảo độc:
  • Mẫu sữa tiêu chuẩn: Đo các thành phần $R_d, T_d, T_c$ trên các nồng độ sữa khác nhau (Bảng 4.7). Từ đó tính toán các thông số $\mu_a, \mu_s, g$ theo cả hai phương pháp KM và MC (Bảng 4.8), cho thấy sự tương thích chặt chẽ với các công bố lý sinh quốc tế đã có.
  • Mẫu vi tảo độc hại: Đo phổ huỳnh quang phụ thuộc bước sóng kích thích (Bảng 4.9) và phổ kích thích huỳnh quang (Bảng 4.10) của các loài tảo thu thập tại vùng biển Đồ Sơn. Tiến hành đo $R_d, T_d, T_c$ và tính toán tỷ lệ $\mu_s/\mu_a$ theo mật độ tế bào đối với ba loài vi tảo:
    • Tảo Prorocentrum rhathymum: Sự phụ thuộc của $R_d, T_d, T_c$ vào mật độ (Bảng 4.11) và sự thay đổi của tỷ lệ $\mu_s/\mu_a$ theo mật độ tế bào (Bảng 4.12).
    • Tảo Pseudo-nitzschia: Giá trị $R_d, T_d, T_c$ (Bảng 4.13) và tỷ lệ $\mu_s/\mu_a$ theo mật độ (Bảng 4.14).
    • Tảo Alexandrium minutum: Giá trị $R_d, T_d, T_c$ (Bảng 4.15) và tỷ lệ $\mu_s/\mu_a$ theo mật độ (Bảng 4.16).

Bảng tổng hợp cấu trúc các khối kỹ thuật chính của thiết bị được chế tạo trong luận án:

Khối chức năng Linh kiện / Công nghệ sử dụng Vai trò kỹ thuật trong hệ đo
Nguồn phát quang Laser diode chip SLD102 Cung cấp chùm tia đơn sắc ổn định công suất, bước sóng đỏ
Hệ tạo tia Hệ thấu kính chuẩn trực Biến đổi chùm tia phát tán thành chùm sáng song song chuẩn trực
Hệ tán xạ quang học Hệ quả cầu tích phân đôi (#1 và #2) Tách và tích lũy quang thông phản xạ ngược ($R_d$), tán xạ xuôi ($T_d$)
Đầu thu quang điện Photodiode hybrid OPT101 (3 kênh PD) Biến đổi quang thông $R_d, T_d, T_c$ thành tín hiệu điện áp tương tự
Biến đổi dữ liệu Mạch tích hợp ADC chuyên dụng Lấy mẫu và biến đổi tương tự - số đa kênh đồng thời
Xử lý số lõi DSP TMS320C6713 (Texas Instruments) Chạy thuật toán song song, kỹ thuật đường ống, tính KM và MC
Ghép nối bộ nhớ Cổng giao tiếp EMIF, mạch logic CPLD Truyền dữ liệu tốc độ cao giữa ADC và bộ nhớ ngoài của DSP
Giao tiếp người dùng Máy tính PC và phần mềm chuyên dụng Cấu hình tham số, hiển thị kết quả đo và tự động sao lưu số liệu

Bảng tổng hợp kết quả khảo sát các loài vi tảo thực nghiệm trong luận án:

Tên loài vi tảo độc hại Vùng phân lập Dữ liệu quang phổ và thông số đo trong luận án Dữ liệu tỷ lệ tán xạ/hấp thụ
Prorocentrum rhathymum Vùng biển Đồ Sơn Phổ huỳnh quang, kích thích huỳnh quang, bộ ba $R_d, T_d, T_c$ (Bảng 4.11) Tỷ lệ $\mu_s/\mu_a$ biến thiên theo mật độ tế bào (Bảng 4.12)
Pseudo-nitzschia Vùng biển Đồ Sơn Phổ huỳnh quang, kích thích huỳnh quang, bộ ba $R_d, T_d, T_c$ (Bảng 4.13) Tỷ lệ $\mu_s/\mu_a$ biến thiên theo mật độ tế bào (Bảng 4.14)
Alexandrium minutum Vùng biển Đồ Sơn Phổ huỳnh quang, kích thích huỳnh quang, bộ ba $R_d, T_d, T_c$ (Bảng 4.15) Tỷ lệ $\mu_s/\mu_a$ biến thiên theo mật độ tế bào (Bảng 4.16)

Kết quả và những đóng góp mới

Luận án tự xác định các đóng góp mới về mặt khoa học lý luận và thực tiễn kỹ thuật như sau:

  1. Đóng góp về phương pháp luận và lý luận khoa học:
  • Xây dựng thành công phương pháp xác định đồng thời ba thông số quang cơ bản ($\mu_a, \mu_s, g$) của môi trường dung dịch chứa các hạt lơ lửng nồng độ thấp bằng kỹ thuật quang học quả cầu tích phân đôi kết hợp thuật toán xử lý tín hiệu số.
  • Thiết lập lời giải giải tích hoàn chỉnh cho hệ phương trình vi phân Kubelka-Munk hai dòng bức xạ bằng biến đổi Laplace, đưa ra mối liên hệ tường minh giữa các hệ số hấp thụ, tán xạ thực tế với các giá trị đo quang thông của quả cầu tích phân.
  • Xây dựng thuật toán và chương trình mô phỏng Monte Carlo tương thích cấu hình hình học của hệ quả cầu tích phân đôi, cho phép tính toán ngược trực tiếp để xác định hệ số bất đối xứng $g$ cùng lúc với $\mu_a$ và $\mu_s$ ngay trong quá trình đo mà không cần đo đạc độc lập giá trị $g$ trên một thiết bị goniometer riêng biệt.
  1. Đóng góp về mặt kỹ thuật, công nghệ và thực tiễn:
  • Thiết kế, chế tạo thành công thiết bị đo lưu động hoàn chỉnh tích hợp giữa quang học, điện tử tương tự và xử lý số. Hệ thống cho phép vận hành linh hoạt ở hai cấu hình: cấu hình đo nhanh (dựa trên mô hình Kubelka-Munk) phục vụ đánh giá nhanh ngoài hiện trường và cấu hình đo chính xác (dựa trên mô phỏng Monte Carlo).
  • Thiết kế mạch ADC chuyên dụng giao tiếp với DSP TMS320C6713 qua cổng EMIF; tối ưu hóa thuật toán tính toán ma trận và xử lý dữ liệu bằng cách lập trình phối hợp ngôn ngữ C và Assembler tận dụng triệt để kiến trúc song song và kỹ thuật đường ống, giúp rút ngắn đáng kể thời gian tính toán.
  • Đã kiểm chứng thực nghiệm hệ đo trên mẫu chuẩn dung dịch sữa 4% chất béo và thực hiện đo đạc thực tế trên 03 loài vi tảo độc (Prorocentrum rhathymum, Pseudo-nitzschia, Alexandrium minutum) ở các mật độ tế bào khác nhau, chứng minh tính khả thi của việc sử dụng bộ thông số quang và tỷ lệ $\mu_s/\mu_a$ làm chỉ số nhận dạng và quan trắc mật độ tảo ở mức thấp dưới 1000 tế bào/ml.
  1. Kiến nghị và giải pháp đề xuất:
  • Ứng dụng thiết bị đã chế tạo vào mạng lưới các trạm quan trắc tài nguyên môi trường biển ven bờ tại Việt Nam nhằm phát hiện và cảnh báo sớm hiện tượng thủy triều đỏ.
  • Mở rộng ứng dụng của hệ đo sang các lĩnh vực liên ngành: kiểm tra chất lượng nguồn nước sinh hoạt (phát hiện nhanh tạp chất, nấm, vi khuẩn lơ lửng) và lĩnh vực lý sinh y học (xét nghiệm máu phát hiện sớm tế bào lạ dựa trên biến đổi thông số tán xạ).

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Hệ thống đo mới được thử nghiệm tập trung trên đối tượng mẫu dung dịch có nồng độ thấp các hạt tán xạ (mẫu sữa chuẩn pha loãng và 3 chủng vi tảo độc được nuôi cấy trong phòng thí nghiệm); chưa thực hiện đo đạc mở rộng trên toàn bộ 70 loài vi tảo độc hại hiện diện tại vùng biển Việt Nam.
  • Thiết bị sử dụng nguồn phát laser đơn bước sóng (chip laser SLD102 bước sóng đỏ), do đó dữ liệu thông số quang thu được tại một thời điểm mới chỉ phản ánh đặc tính quang học tại một bước sóng cố định, chưa khảo sát liên tục toàn dải phổ quang học.
  • Phép đo thực tế đối với mẫu thủy vực tự nhiên có thể chịu thêm ảnh hưởng phức tạp từ các chất hữu cơ hòa tan có màu (CDOM) và các hạt vô cơ lơ lửng khác trong nước biển tự nhiên chưa được phân tách hoàn toàn.

Hướng nghiên cứu tiếp theo:

  • Nâng cấp nguồn phát quang sang hệ thống laser đa bước sóng hoặc nguồn sáng dải rộng kết hợp đơn sắc tử để thu nhận thông số quang học phụ thuộc bước sóng ($\mu_a(\lambda), \mu_s(\lambda)$).
  • Xây dựng cơ sở dữ liệu thông số quang chuẩn cho đa dạng các loài vi tảo biển tại các vùng sinh thái biển khác nhau của Việt Nam.
  • Triển khai thử nghiệm thiết bị trực tiếp trên các tàu nghiên cứu khảo sát ngoài hiện trường thực tế để hoàn thiện thiết kế cơ khí chống chịu môi trường biển.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn đối với nhiều đối tượng nghiên cứu và ứng dụng:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý Vô tuyến, Điện tử học, Kỹ thuật Quang điện tử và Kỹ thuật Viễn thông: Chương 2 và Chương 3 cung cấp tài liệu chi tiết về phương pháp giải phương trình truyền bức xạ, biến đổi Laplace cho mô hình Kubelka-Munk và kỹ thuật mô phỏng Monte Carlo quá trình di chuyển photon.
  • Kỹ sư phần cứng và lập trình nhúng: Chương 4 cung cấp sơ đồ thiết kế chi tiết, mạch nguyên lý ghép nối ADC với DSP TMS320C6713 qua giao tiếp EMIF, phương pháp lập trình tối ưu hóa song song bằng C và Assembler, cùng các kết quả khảo sát thời gian xử lý thực tế.
  • Các nhà nghiên cứu hải dương học, sinh thái học biển và cán bộ quản lý môi trường: Chương 1 và Chương 4 cung cấp dữ liệu thực nghiệm về đặc tính quang học, phổ huỳnh quang và sự biến đổi tỷ lệ tán xạ/hấp thụ theo mật độ của ba loài vi tảo độc hại vùng biển Đồ Sơn (Prorocentrum rhathymum, Pseudo-nitzschia, Alexandrium minutum), làm cơ sở xây dựng phương pháp quan trắc quang học thay thế quy trình đếm tế bào truyền thống.

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

1. Luận án giải quyết vấn đề gì trong quan trắc môi trường thủy vực ven bờ? Luận án giải quyết bài toán phát hiện sớm sự bùng phát của vi tảo độc hại ở mật độ thấp (dưới 1000 tế bào/ml). Bằng việc chế tạo thiết bị đo quang học lưu động sử dụng quả cầu tích phân đôi và DSP, hệ thống cho phép đo trực tiếp tại chỗ ba thông số quang ($\mu_a, \mu_s, g$) mà không cần trải qua quy trình lưu trữ, vận chuyển mẫu kéo dài nhiều tuần như các phương pháp truyền thống (Utermöhl, Sedgewick Rafter, ELISA), đồng thời không bị nhiễu nền từ khí quyển như phương pháp viễn thám vệ tinh.

2. Hai cấu hình làm việc của hệ đo (đo nhanh và đo chính xác) hoạt động dựa trên nguyên lý nào?

  • Cấu hình đo nhanh: Dựa trên mô hình Kubelka-Munk hai dòng bức xạ kết hợp với biến đổi Laplace. Cấu hình này sử dụng các biểu thức giải tích tường minh để tính nhanh hệ số hấp thụ $K$ và tán xạ $S$, chuyển đổi sang $\mu_a, \mu_s$, thích hợp cho việc đánh giá nhanh ngoài hiện trường.
  • Cấu hình đo chính xác: Dựa trên mô phỏng ngẫu nhiên Monte Carlo kết hợp kỹ thuật quả cầu tích phân đôi. Cấu hình này mô phỏng chi tiết từng bước di chuyển, góc lệch tán xạ của hàng triệu photon trong cuvet mẫu có xét đến phản xạ khúc xạ bề mặt, cho phép xác định chính xác đồng thời cả ba thông số $\mu_a, \mu_s$ và hệ số bất đối xứng $g$.

3. Phần cứng xử lý tín hiệu của hệ đo sử dụng những linh kiện chính nào? Hệ thống sử dụng bộ vi xử lý tín hiệu số DSP TMS320C6713 (Texas Instruments) trên bo mạch DSK TMS320C6713, tích hợp với khối biến đổi tương tự - số ADC chuyên dụng tự thiết kế thông qua bus mở rộng bộ nhớ ngoài EMIF và CPLD. Khối quang học bao gồm nguồn phát chip laser diode SLD102 có mạch điều khiển dòng ổn định công suất, hệ quang chuẩn trực, hai quả cầu tích phân chế tạo bằng kim loại phủ lớp khuếch tán trong lòng cầu, và ba đầu thu photodiode hybrid OPT101.

4. Mẫu thực nghiệm nào đã được sử dụng để kiểm chứng độ chính xác của hệ đo? Tác giả đã sử dụng hai nhóm mẫu thực nghiệm:

  • Nhóm mẫu chuẩn lý sinh: Dung dịch sữa có tỷ lệ chất béo 4% pha loãng ở các nồng độ khác nhau để tạo môi trường tán xạ hạt cầu ổn định; kết quả đo $\mu_a, \mu_s, g$ bằng KM và MC được đối chiếu khớp với các số liệu chuẩn quốc tế đã công bố.
  • Nhóm mẫu sinh học: 03 loài vi tảo độc hại được phân lập và nuôi cấy từ vùng biển Đồ Sơn gồm Prorocentrum rhathymum, Pseudo-nitzschia và Alexandrium minutum; hệ đo đã xác định chính xác các giá trị $R_d, T_d, T_c$ và quy luật biến thiên của tỷ lệ $\mu_s/\mu_a$ theo mật độ tế bào của từng loài.

5. Vì sao tác giả phải tối ưu hóa chương trình tính toán trên DSP bằng Assembler và cấu trúc song song? Thuật toán mô phỏng Monte Carlo và xử lý dữ liệu ma trận đòi hỏi khối lượng phép tính dấu phẩy động rất lớn. Bằng cách khai thác kiến trúc VLIW của DSP TMS320C6713, kết hợp kỹ thuật xử lý song song, kỹ thuật đường ống lệnh và viết mã bằng Assembler (thay vì chỉ dùng C thuần túy), thời gian tính toán của các phép nhân ma trận và vòng lặp mô phỏng đã được giảm thiểu đáng kể (được định lượng tại Bảng 4.1 và các đồ thị trong Chương 4), đảm bảo thiết bị hiển thị kết quả đo đạc gần với thời gian thực.

Kết luận

Luận án của Nguyễn Tuấn Anh đã nghiên cứu và phát triển thành công hệ đo quang học hoàn chỉnh xác định bộ ba thông số quang học ($\mu_a, \mu_s, g$) của các mẫu dung dịch có nồng độ tán xạ thấp. Công trình kết hợp chặt chẽ giữa cơ sở lý thuyết truyền bức xạ, thuật toán giải tích Kubelka-Munk, mô phỏng Monte Carlo với kỹ thuật chế tạo quang điện tử và xử lý tín hiệu số tốc độ cao trên chip DSP TMS320C6713. Các kết quả kiểm chứng thực nghiệm trên mẫu dung dịch sữa chuẩn và ba loài vi tảo độc hại vùng biển Đồ Sơn đã chứng minh tính đúng đắn của phương pháp, mở ra khả năng ứng dụng thực tế trong việc quan trắc, cảnh báo sớm hiện tượng tảo nở hoa cũng như các ứng dụng phân tích hạt trong lý sinh học và kiểm định chất lượng nước.