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 thực hiện bởi nghiên cứu sinh Nguyễn Tuấn Anh dưới sự hướng dẫn khoa học của TS. Bạch Gia Dương và PGS. TS. Vũ Anh Phi tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội vào năm 2011.

Vấn đề ô nhiễm môi trường thủy vực ven bờ biển do sự bùng phát của các loài vi tảo độc, hại (hiện tượng nở hoa của nước hoặc thủy triều đỏ) là một nguy cơ sinh thái nghiêm trọng. Trong điều kiện tự nhiên, mật độ tế bào tảo độc, hại thường duy trì ở mức thấp, từ vài chục đến vài trăm tế bào/ml. Tuy nhiên, khi môi trường nước chuyển từ trạng thái nghèo dưỡng sang giàu dưỡng, các loài tảo này có thể bùng phát số lượng lên đến hàng triệu tế bào/ml chỉ trong vòng vài ngày. Sau giai đoạn bùng phát, tảo bước vào giai đoạn thoái hóa, phân hủy, tiêu thụ lượng lớn oxy hòa tan và giải phóng độc tố vào môi trường, đe dọa trực tiếp đến hệ sinh thái, nguồn lợi thủy sản và hoạt động du lịch. Theo khuyến cáo của Tổ chức Y tế Thế giới (WHO), khi mật độ tế bào đạt ngưỡng 1000 tế bào/ml (đối với vi khuẩn lam/tảo lam), nguồn nước được xem là bị ô nhiễm nặng nề. Do đó, việc phát hiện và quan trắc tảo ở mật độ thấp dưới ngưỡng 1000 tế bào/ml là yêu cầu cấp thiết để đưa ra cảnh báo sớm.

Khoảng trống nghiên cứu thực tế tại thời điểm thực hiện luận án là các phương pháp quan trắc hiện hành tại Việt Nam chủ yếu được tiến hành trong phòng thí nghiệm thông qua các quy trình truyền thống mất nhiều tuần lễ, không đáp ứng được yêu cầu cảnh báo nhanh ngoài hiện trường. Mặt khác, các phương pháp hiện đại trên thế giới như phân tích ảnh viễn thám vệ tinh thì chi phí đắt đỏ và phụ thuộc điều kiện thời tiết, còn phương pháp đo phổ huỳnh quang chỉ phân loại được khoảng 5 nhóm tảo theo sắc tố. Luận án lựa chọn giải quyết khoảng trống này bằng cách xây dựng một hệ đo quang học tại chỗ, cho phép xác định trực tiếp và đồng thời các thông số quang cốt lõi của thủy vực ở dải mật độ tảo thấp.

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

  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 phương pháp, mô hình phù hợp và dựa trên công nghệ tích hợp quang điện tử, kỹ thuật xử lý tín hiệu. Hệ đo làm việc theo hai cấu hình: đo nhanh và đo chính xác ba thông số quang 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$).
  2. Trên cơ sở hệ đo đã xây dựng, chế tạo một thiết bị đo lưu động, cho phép đo nhanh, đo đồng thời ba thông số quang của thủy vực có mật độ tảo thấp.

Đối tượng nghiên cứu dưới góc độ vật lý là bộ ba thông số quang của môi trường dung dịch: hệ số hấp thụ ($\mu_a$), hệ số tán xạ ($\mu_s$) và hệ số bất đối xứng ($g$). Về mặt thực nghiệm, đối tượng khảo sát gồm dung dịch sữa có tỷ lệ chất béo 4% pha loãng ở các nồng độ khác nhau (được sử dụng làm mẫu chuẩn lý sinh nhờ các hạt béo có dạng hình cầu kích thước ổn định) và 03 mẫu tảo độc đã được phân lập, nuôi cấy ở các dải mật độ khác nhau: Prorocentrum rhathymum, Pseudo-nitzschia và Alexandrium minutum.

Phạm vi nghiên cứu tập trung vào môi trường thủy vực ven bờ biển Việt Nam, cụ thể là các mẫu tảo độc hại thu thập tại vùng thủy vực Đồ Sơn (Hải Phòng), khảo sát ở dải mật độ tế bào thấp dưới 1000 tế bào/ml.

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

Phương pháp xác định đối tượng thông qua thông số quang học có lịch sử từ những năm 70 của thế kỷ 18 với việc sử dụng đĩa Secchi thả chìm dưới nước để đo độ trong thông qua ánh sáng phản xạ. Trải qua quá trình phát triển, các kỹ thuật đo thông số quang được ứng dụng rộng rãi trong công nghệ sinh học, chế biến thực phẩm và đặc biệt là chẩn đoán y tế (nghiên cứu mô, tế bào qua hệ số hấp thụ, tán xạ và bất đối xứng).

Theo tổng thuật của luận án, các nghiên cứu về tảo độc, hại trong và ngoài nước tập trung vào bốn hướng chính:

  • Điều tra, phát hiện vùng xuất hiện tảo độc, hại và các hiệu ứng của độc tố.
  • Nghiên cứu về những loài tảo gây độc và phương pháp phát hiện chúng, tác động của điều kiện môi trường lên sự nở hoa và sự tích luỹ độc tố.
  • Nghiên cứu về độc tố của tảo và các nhóm độc tố.
  • Giám sát và quản lý tảo độc, hại.

Tại Việt Nam, hoạt động quan trắc định kỳ được thực hiện bởi 4 Trung tâm quan trắc và cảnh báo môi trường (3 trung tâm nội địa tại ba miền Bắc, Trung, Nam và 1 trung tâm vùng biển) với tần suất 2 lần/năm (tháng 4-5 và tháng 9-10). Bộ thông số môi trường áp dụng theo TCVN 5943-1995 và tiêu chuẩn của Bộ Thủy sản năm 2006. Việc phân tích tảo độc, hại dựa trên các phương pháp truyền thống trong phòng thí nghiệm, bao gồm:

  • Phương pháp Utermöhl: sử dụng bình trụ kết hợp buồng chứa phẳng đáy tháo rời để lắng mẫu và quan sát dưới kính hiển vi.
  • Phương pháp buồng đếm Sedgewick Rafter: buồng chứa dung tích 100 $\text{mm}^3$ có đáy chia lưới $1 \times 1 \text{ mm}^2$ để đếm số lượng vi tảo trên mỗi microlít.
  • Kỹ thuật nhuộm màu: sử dụng phẩm nhuộm Carmin (nhận diện vách cellulose) hoặc phẩm nhuộm xanh iod, pectin (nhận diện lignin hoặc suberin) sau khi xử lý hóa chất loại bỏ tế bào chất.
  • Phương pháp miễn dịch liên kết enzyme (ELISA): sử dụng kháng thể đặc hiệu gắn đáy giếng để nhận diện kháng nguyên thông qua phản ứng cạnh tranh đổi màu.

Tác giả chỉ rõ các phương pháp truyền thống tồn tại ba nhược điểm lớn: quy trình 4 công đoạn (thu thập mẫu bằng tàu thuyền trang bị GPS, bảo quản mẫu, vận chuyển về phòng thí nghiệm, lưu trữ mẫu) tốn nhiều nhân lực và thời gian; thời gian phân tích kéo dài hàng tuần nên không thể phục vụ cảnh báo nhanh; số liệu ghi chép thủ công dễ phát sinh sai sót.

Đối với xu thế quan trắc hiện đại trên thế giới, luận án điểm lại hai hướng tiếp cận chính:

  1. Kết hợp số liệu thực địa với dữ liệu từ các cảm biến quang học đặt trên vệ tinh (như SeaWiFS, MODIS, MERIS). Phương pháp này dựa trên phân tích phổ phản xạ bức xạ $R(\lambda)$ từ bề mặt đại dương: $$R(\lambda) = f \frac{b_b(\lambda)}{a(\lambda) + b_b(\lambda)}$$ Trong đó $f$ là hằng số môi trường, $b_b(\lambda)$ là hệ số tán xạ ngược, $a(\lambda)$ là hệ số hấp thụ tổng hợp gồm 4 thành phần: nước ($a_w$), chất hữu cơ hòa tan có màu ($a_y$ - CDOM), sinh vật phù du ($C_{ph} a_{ph}$) và phần tử vô cơ ($C_s a_s$). Phương pháp này có ưu thế bao quát diện rộng nhưng chi phí vận hành cao, thuật toán xử lý ảnh phức tạp và chịu ảnh hưởng lớn bởi góc chiếu mặt trời, bức xạ khí quyển và phản xạ đáy biển.
  2. Thiết bị đo phổ huỳnh quang của tảo: dựa trên sự phát quang của các sắc tố Chlorophyll a (hấp thụ 436 nm, phát huỳnh quang 685 nm), Chlorophyll b (hấp thụ 645 nm) và Chlorophyll c (hấp thụ 630 nm). Các công trình tiêu biểu gồm dự án của Đại học Tổng hợp Umaine (Hoa Kỳ, triển khai từ 2007) và thiết bị của hãng bbe (Cộng hòa Liên bang Đức) sử dụng 5 đèn LED kích thích ở các bước sóng 450 nm, 525 nm, 570 nm, 590 nm và 610 nm kết hợp đầu thu nhân quang điện. Hạn chế cốt tử của phương pháp huỳnh quang là chỉ nhận diện được khoảng 5 nhóm tảo do nhiều loài có chung nhóm sắc tố.

Vị trí của luận án được xác lập tại điểm giao thoa giữa vật lý thực nghiệm và quan trắc môi trường: thay vì dựa vào huỳnh quang hay ảnh vệ tinh, tác giả đề xuất đo trực tiếp các thông số quang học tán xạ và hấp thụ của mẫu nước biển khi bị chiếu chùm tia laser chuẩn trực. Bằng việc kết hợp kỹ thuật quả cầu tích phân đôi với mô hình Kubelka-Munk và mô phỏng Monte Carlo, hệ đo giải quyết được bài toán đo nhanh tại hiện trường và đo chính xác các mẫu có mật độ vi tảo rất thấp.

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

Luận án dựa trên lý thuyết tương tác giữa photon và môi trường sinh học. Khi chùm sáng truyền qua môi trường chứa các hạt phân tán, cường độ bức xạ bị suy giảm bởi hai quá trình cơ bản: hấp thụ và tán xạ.

Quá trình hấp thụ ánh sáng tuân theo định luật Beer-Lambert: $$I = I_0 e^{-\mu_a l} = I_0 e^{-\rho \sigma_a l}$$ Trong đó $I_0$ là cường độ chùm sáng tới, $I$ là cường độ sau khi truyền qua lớp mẫu có độ dày $l$, $\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ụ. Đại lượng $1/\mu_a$ biểu thị quãng đường tự do trung bình của photon trước khi bị hấp thụ.

Quá trình tán xạ ánh sáng được phân thành tán xạ đàn hồi (năng lượng photon không đổi sau va chạm, gồm tán xạ Rayleigh khi hệ số kích thước $x = ka = \frac{2\pi n}{\lambda} a \ll 1$, và tán xạ Mie cho các hạt hình cầu kích thước bất kỳ) và tán xạ không đàn hồi (photon biến đổi năng lượng, gồm hiện tượng phát huỳnh quang với thời gian sống $10^{-9} - 10^{-6}$ s, phát lân quang $10^{-6} - 10^2$ s, và tán xạ Raman với vạch Stokes và đối Stokes). Luận án quy ước xem xét tán xạ trong mẫu vi tảo là tán xạ đàn hồi do thành phần phi đàn hồi chiếm tỷ trọng không đáng kể. Tán xạ được đặc trưng bởi hệ số tán xạ $\mu_s$ ($\text{mm}^{-1}$), mật độ xác suất tán xạ $p(x) = \mu_s e^{-\mu_s x}$, 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_{tr} = \mu_a + \mu_s(1-g) = \mu_a + \mu_s'$.

Hệ số bất đối xứng $g$ (anisotropy factor) được định nghĩa qua giá trị trung bình của cosin góc tán xạ $\theta$: $$g = \int_{-1}^{1} \cos\theta , p(\cos\theta) , d(\cos\theta)$$ Luận án sử dụng hàm pha tán xạ Henyey-Greenstein ($P_{HG}$) để mô tả phân bố góc tán xạ: $$P_{HG}(\theta) = \frac{1 - g^2}{4\pi (1 + g^2 - 2g \cos\theta)^{3/2}}$$

Quá trình lan truyền ánh sáng trong môi trường sinh học được mô tả tổng quát bằng phương trình truyền bức xạ (Radiative Transfer Equation - RTE): $$\frac{dL(r, \hat{s})}{ds} = -\mu_a L(r, \hat{s}) - \mu_s L(r, \hat{s}) + \mu_s \int_{4\pi} p(\hat{s}, \hat{s}') L(r, \hat{s}') , d\omega' + S(r, \hat{s})$$ Trong đó $L(r, \hat{s})$ là độ bức xạ tại vị trí $r$ theo hướng $\hat{s}$, $S(r, \hat{s})$ là nguồn bức xạ nội tại, và $d\omega'$ là góc đặc vi phân.

Để giải phương trình truyền bức xạ trong điều kiện thực nghiệm, luận án áp dụng hai phương pháp:

  1. Mô hình Kubelka-Munk (KM) hai dòng bức xạ: giả thiết môi trường gồm các lớp mỏng phẳng song song vô hạn, môi trường đồng nhất quang học, tán xạ bán cầu ngẫu nhiên đẳng hướng. Hệ phương trình vi phân mô tả dòng bức xạ xuôi $i_t(x)$ và dòng bức xạ ngược $i_r(x)$ có dạng: $$\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)$$ Hệ phương trình này được tác giả giải chính xác bằng phương pháp biến đổi Laplace. Với điều kiện biên lớp mẫu có độ dày vô hạn, tỷ số giữa hệ số hấp thụ KM ($K$) và hệ số tán xạ KM ($S$) liên hệ với hệ số phản xạ $R_\infty$ qua hệ thức: $$\frac{K}{S} = \frac{(1 - R_\infty)^2}{2 R_\infty}$$ Từ nghiệm giải tích kết hợp điều kiện biên của mẫu mỏng thực nghiệm, các hệ số $K$ và $S$ được tính toán trực tiếp từ thành phần phản xạ khuếch tán ($R_d$) và truyền qua khuếch tán ($T_d$).

  2. Mô phỏng Monte Carlo (MC): theo dõi đường đi ngẫu nhiên của từng photon trong môi trường tán xạ và hấp thụ. Các biến số ngẫu nhiên được lấy mẫu từ biến ngẫu nhiên đồng nhất $\xi \in (0, 1]$ do máy tính tạo ra:

  • Quãng đường tự do trung bình giữa hai sự kiện tương tác: $s = -\frac{\ln(\xi)}{\mu_t}$.
  • Góc lệch tán xạ $\theta$ từ hàm pha 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ị: $\psi = 2\pi\xi$.

Kỹ thuật thực nghiệm trọng tâm là phương pháp quả cầu tích phân đôi. Quả cầu tích phân có mặt trong phủ lớp vật liệu phản xạ khuếch tán cao $\rho \approx 100%$, tạo ra nguồn sáng khuếch tán Lambertian. Tổng quang thông bên trong quả cầu sau vô hạn lần phản xạ đạt giá trị: $$\phi_s = \frac{\phi_i \rho (1 - f)}{1 - \rho(1 - f)}$$ Với $f$ là tỷ lệ diện tích khe hở trên diện tích toàn phần quả cầu, hệ số nhân của quả cầu tích phân là $M = \frac{\rho}{1 - \bar{\rho}}$, và hằng số thời gian đáp ứng là $\tau = -\frac{2}{3} \frac{D_s}{c \ln\rho}$ ($D_s$ là đường kính quả cầu, $c$ là tốc độ ánh sáng).

Hệ đo sử dụng ba đầu thu quang học để ghi nhận ba đại lượng vật lý:

  • Thành phần phản xạ ngược khuếch tán: $R_d = \frac{X - Y}{Z - Y}$ (với $X$ là cường độ khi có mẫu, $Y$ là tín hiệu nền khi laser tắt, $Z$ là cường độ khi đặt thanh chuẩn phản xạ 100%).
  • Thành phần truyền qua khuếch tán: $T_d = \frac{A_1}{A_2}$ (với $A_1, A_2$ là tín hiệu đầu thu tương ứng khi có mẫu và không có mẫu).
  • Thành phần truyền qua chuẩn trực: $T_c = \frac{I}{I_0}$.

Nguồn dữ liệu thực nghiệm bao gồm: các tín hiệu điện chuẩn (kiểm tra độ chính xác của bộ chuyển đổi ADC và DSP); dung dịch sữa chuẩn 4% chất béo pha loãng theo các nồng độ khác nhau; và 03 mẫu vi tảo biển độc hại thu thập tại vùng biển Đồ Sơn được phân lập và nuôi cấy tại phòng thí nghiệm.

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 tập trung trình bày tổng quan về đối tượng thực nghiệm và các vấn đề liên quan đến hiện tượng ô nhiễm môi trường thủy vực do vi tảo bùng phát. Tác giả khảo sát thực trạng khoảng 70 loài tảo độc, hại hiện diện tại vùng ven biển Việt Nam, đồng thời phân tích hệ thống quy chuẩn quốc gia gồm TCVN 5943-1995 và quy chuẩn ngành của Bộ Thủy sản năm 2006.

Nội dung chương phân tích chi tiết quy trình quan trắc truyền thống gồm 4 công đoạn: thu thập mẫu tại hiện trường bằng lưới vớt trên tàu thuyền định vị GPS, bảo quản mẫu trong các bình chuyên dụng, vận chuyển mẫu về phòng thí nghiệm và lưu trữ mẫu trước khi phân tích. Các phương pháp định lượng truyền thống được mổ xẻ gồm phương pháp Utermöhl, buồng đếm Sedgewick Rafter, kỹ thuật nhuộm màu (Carmin, xanh iod) và phương pháp miễn dịch ELISA. Tác giả chỉ ra rằng các kỹ thuật này đòi hỏi thời gian phân tích từ vài ngày đến vài tuần, phụ thuộc kỹ năng của kiểm nghiệm viên và dễ phát sinh sai số do ghi chép thủ công.

Tiếp đó, chương đánh giá xu thế quan trắc trên thế giới, bao gồm việc phân tích ảnh phản xạ thu từ các cảm biến vệ tinh SeaWiFS, MODIS, MERIS và các hệ đo phổ huỳnh quang tảo (như dự án quan trắc nguồn nước của Đại học Umaine năm 2007 và thiết bị 5 bước sóng của hãng bbe Đức). Trên cơ sở chỉ ra hạn chế của các phương pháp trên, tác giả đề xuất mô hình quan trắc thủy vực dự kiến sử dụng hệ hai quả cầu tích phân đặt nối tiếp nhau, một nguồn laser bán dẫn chuẩn trực và ba đầu thu photodiode (PD1 đo $R_d$, PD2 đo $T_d$, PD3 đo $T_c$). Tín hiệu quang điện sau đó được đưa vào khối xử lý số để tính toán bộ ba thông số quang $\mu_a, \mu_s, g$ theo hai chế độ: đo nhanh (dựa trên mô hình Kubelka-Munk) và đo chính xác (dựa trên mô phỏng Monte Carlo).

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

Chương 2 thiết lập toàn bộ nền tảng lý thuyết vật lý cho hệ đo. Tác giả trình bày các tương tác cơ bản giữa photon và môi trường sinh học, phân tích cơ chế dịch chuyển mức năng lượng phân tử trong quá trình hấp thụ ánh sáng, định luật Beer-Lambert, tiết diện hấp thụ $\sigma_a$ và hệ số hấp thụ $\mu_a$. Trong phần lý thuyết tán xạ, tác giả làm rõ sự khác biệt giữa tán xạ đàn hồi (Rayleigh, Mie) và tán xạ phi đàn hồi (huỳnh quang, lân quang, tán xạ Raman với mức dịch chuyển Stokes và đối Stokes).

Chương này trình bày chi tiết các thông số mô tả trường bức xạ khuếch tán: mật độ xác suất tán xạ $p(x)$, quãng đường tự do trung bình $1/\mu_s$, hàm tán xạ $p(e_s, e_s')$, hệ số bất đối xứng $g$, hệ số suy giảm tổng thể $\mu_t$, hệ số suy giảm rút gọn $\mu_{tr}$ và hàm pha Henyey-Greenstein. Phương trình truyền bức xạ tổng quát được dẫn giải chi tiết cho thể tích vi phân $ds , dA$.

Các mô hình gần đúng được phân tích có mô hình Pierce-Marcus, mô hình Hapke, mô hình Kubelka-Munk (KM) hai dòng bức xạ ($i_t, i_r$) và mô phỏng Monte Carlo (MC). Đồng thời, chương 2 xây dựng lý thuyết quả cầu tích phân: định luật bức xạ Lambertian, hệ số trao đổi bức xạ vi phân $dF_{d1-d2}$, phương trình bức xạ trong quả cầu tích phân, hệ số nhân $M$, hằng số thời gian $\tau$ và thiết lập sơ đồ nguyên lý đo ba thành phần quang học $R_d, T_d, T_c$. Đối với mô phỏng Monte Carlo, tác giả trình bày cơ sở toán học của việc tạo số ngẫu nhiên $\xi$, kỹ thuật lấy mẫu biến ngẫu nhiên từ hàm mật độ xác suất để xác định bước dịch chuyển $s$, góc lệch $\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

Chương 3 giải quyết bài toán cốt lõi về phương pháp luận và thuật toán tính toán. Tác giả áp dụng mô hình Kubelka-Munk cho mẫu dung dịch mỏng, sử dụng phép biến đổi Laplace để giải hệ hai phương trình vi phân tuyến tính không đồng nhất của KM: $$p I(p) - i_r(0) = (K + S)I(p) - S J(p)$$ $$p J(p) - i_t(0) = S I(p) - (K + S)J(p)$$ Bằng cách áp dụng phép biến đổi Laplace ngược và giải hệ phương trình với các điều kiện biên thực tế (xét đến hệ số phản xạ tại mặt trước $r_0$ và mặt sau $r_h$ của cuvet chứa mẫu), tác giả rút ra công thức nghiệm giải tích tường minh xác định các hệ số hấp thụ KM ($K$) và tán xạ KM ($S$). Đồng thời, chương này thiết lập các phương trình chuyển đổi giải tích từ các hệ số $K, S$ sang các đại lượng vật lý cơ bản $\mu_a, \mu_s$.

Về thiết kế phần cứng quang học, chương 3 trình bày các tính toán tối ưu hình học cho quả cầu tích phân: đường kính quả cầu, diện tích khe vào/ra để giảm thiểu tỷ lệ khe $f$, tính toán chất liệu lớp phủ phản xạ bên trong và thiết kế quang hệ tạo chùm tia laser chuẩn trực có góc phân kỳ cực tiểu.

Đối với mô phỏng Monte Carlo, tác giả xây dựng lưu đồ thuật toán chuyên biệt cho hệ quả cầu tích phân đôi: mô phỏng quá trình photon di chuyển, tán xạ và hấp thụ trong mẫu dung dịch; tính toán hiện tượng phản xạ toàn phần và khúc xạ tại ranh giới dung dịch - vách cuvet thủy tinh - không khí theo định luật Fresnel; theo dõi số lượng photon thoát ra ở nửa không gian bán cầu trước (đi vào quả cầu tích phân thứ nhất tạo nên $R_d$), photon thoát ra ở bán cầu sau (đi vào quả cầu tích phân thứ hai tạo nên $T_d$), và photon truyền thẳng không tương tác (đi tới đầu thu PD3 tạo nên $T_c$). Thuật toán này cho phép giải bài toán ngược: từ ba giá trị đo thực nghiệm $R_d, T_d, T_c$, chương trình tính lặp để tìm ra bộ ba thông số $\mu_a, \mu_s$ và $g$ của mẫu 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 trọng tâm, trình bày toàn bộ quá trình chế tạo thiết bị, thử nghiệm mạch điện tử, tối ưu hóa thuật toán và đo đạc trên mẫu thực tế.

Tác giả thiết kế sơ đồ khối phần cứng hệ đo hoàn chỉnh gồm hai cấu hình:

  • Cấu hình đo nhanh: ứng dụng mô hình Kubelka-Munk, cho phép tính toán trực tiếp kết quả thông số quang trong thời gian ngắn, thích hợp cho việc khảo sát lưu động tại hiện trường.
  • Cấu hình đo chính xác: ứng dụng mô phỏng Monte Carlo giải bài toán ngược kết hợp dữ liệu từ quả cầu tích phân đôi, cho kết quả có độ chính xác cao đối với các mẫu có nồng độ tán xạ cực thấp.

Trung tâm xử lý tín hiệu số của hệ thống là bo mạch phát triển DSK TMS320C6713 (Texas Instruments) hoạt động ở tần số xung nhịp cao, có cấu trúc song song VLIW (Very Long Instruction Word) và kỹ thuật đường ống (pipeline). Tác giả đã khảo sát khả năng xử lý số của chip thông qua phép nhân hai ma trận $n \times n$ (Bảng 4.1), so sánh thời gian thực thi giữa chương trình viết bằng C thuần túy và chương trình viết bằng Hợp ngữ (Assembler) tối ưu hóa cấu trúc song song, chứng minh việc tối ưu mã lệnh giúp giảm đáng kể thời gian tính toán.

Để thu nhận dữ liệu đa kênh đồng thời từ các cảm biến quang điện, tác giả đã thiết kế và chế tạo bo mạch biến đổi tương tự - số (ADC) chuyên dụng, tích hợp với DSK TMS320C6713 thông qua giao diện bộ nhớ ngoài (EMIF) và khối logic khả trình CPLD. Mạch ADC được kiểm tra độ chính xác với các nguồn tín hiệu điện chuẩn lối vào; kết quả đo kiểm ghi nhận sự sai khác rất nhỏ về biên độ (Bảng 4.2) và tần số (Bảng 4.3) giữa tín hiệu phát chuẩn và giá trị hiển thị trên phần mềm.

Hệ thống quang điện tử bao gồm:

  • Nguồn phát quang: sử dụng chip laser diode SLD102 (Bảng 4.5) đi kèm mạch ổn định dòng và ổn định nhiệt độ để đảm bảo công suất phát không suy giảm theo thời gian; đáp ứng tiêu chuẩn nguồn sáng cho kiểm tra chất lượng nguồn nước (Bảng 4.4).
  • Đầu thu quang điện: sử dụng 3 cảm biến photodiode hybrid tích hợp mạch khuếch đại thuật toán OPT101 (Bảng 4.6), có độ nhạy cao, độ tuyến tính tốt và dải tần đáp ứng phù hợp.
  • Buồng đo: hệ hai quả cầu tích phân chế tạo bằng hợp kim nhẹ, phủ vật liệu phản xạ khuếch tán cao ở mặt trong, tích hợp cuvet mẫu đặt giữa hai quả cầu.

Phần mềm điều khiển được xây dựng hoàn chỉnh với giao diện đồ họa cho phép thiết lập cấu hình đo, hiển thị trực quan các tín hiệu điện áp thu được, tự động tính toán $\mu_a, \mu_s, g$ theo thuật toán đã chọn và lưu trữ cơ sở dữ liệu vào máy tính.

Tác giả tiến hành đo thực nghiệm trên hai nhóm đối tượng:

  1. Mẫu dung dịch sữa chuẩn 4% chất béo: đo các thành phần $R_d, T_d, T_c$ ở nhiều nồng độ pha loãng khác nhau (Bảng 4.7). Từ đó, tính toán các thông số $\mu_a, \mu_s$ và $g$ theo cả hai phương pháp KM và MC (Bảng 4.8), đối chiếu với các công bố quốc tế về tán xạ hạt mỡ sữa để kiểm chuẩn độ tin cậy của thiết bị.
  2. Mẫu vi tảo biển độc hại: tiến hành đo phổ huỳnh quang theo 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) trên các mẫu tảo thu thập tại Đồ Sơn. Tiếp đó, hệ đo quả cầu tích phân đôi thực hiện đo các thành phần $R_d, T_d, T_c$, từ đó xác định $\mu_a, \mu_s, g$ và khảo sát sự phụ thuộc của tỷ lệ $\mu_s / \mu_a$ theo mật độ tế bào đối với 3 loài:
  • Loài Prorocentrum rhathymum: số liệu đo $R_d, T_d, T_c$ theo mật độ (Bảng 4.11) và sự biến thiên tỷ lệ $\mu_s / \mu_a$ (Bảng 4.12).
  • Loài Pseudo-nitzschia: số liệu đo $R_d, T_d, T_c$ theo mật độ (Bảng 4.13) và sự biến thiên tỷ lệ $\mu_s / \mu_a$ (Bảng 4.14).
  • Loài Alexandrium minutum: số liệu đo $R_d, T_d, T_c$ theo mật độ (Bảng 4.15) và sự biến thiên tỷ lệ $\mu_s / \mu_a$ (Bảng 4.16).

Dưới đây là bảng tổng hợp các khối phần cứng chính được thiết kế, chế tạo và tích hợp trong hệ thống:

Khối chức năng Linh kiện / Thiết bị sử dụng Thông số kỹ thuật và vai trò trong hệ đo
Nguồn phát quang Laser Diode SLD102 Nguồn phát chuẩn trực, ổn định công suất, nhiệt độ phòng (Bảng 4.5)
Cảm biến quang học Photodiode Hybrid OPT101 Đầu thu tích hợp khuếch đại, độ nhạy cao, độ ồn thấp (Bảng 4.6), gồm 3 kênh
Xử lý tín hiệu số DSK TMS320C6713 (Texas Instruments) DSP dấu phẩy động 32-bit, kiến trúc VLIW song song, kỹ thuật đường ống
Chuyển đổi tương tự - số Mạch ADC chế tạo riêng + CPLD Tích hợp qua cổng bộ nhớ ngoài EMIF, đọc đồng thời tín hiệu từ các cảm biến
Hệ buồng đo quang học Quả cầu tích phân đôi Thu nhận tán xạ ngược ($R_d$), tán xạ xuôi ($T_d$) và truyền thẳng ($T_c$)
Giao diện người dùng Phần mềm trên PC (giao tiếp DSK) Thiết lập cấu hình đo, hiển thị đồ thị tín hiệu, tính toán $\mu_a, \mu_s, g$ và sao lưu

Bảng so sánh hai cấu hình làm việc của hệ đo được tác giả xây dựng:

Tiêu chí so sánh Cấu hình đo nhanh Cấu hình đo chính xác
Cơ sở mô hình toán học Mô hình Kubelka-Munk (KM) hai dòng bức xạ Mô phỏng Monte Carlo (MC) lan truyền photon
Thời gian tính toán Rất ngắn (tính toán đại số trực tiếp) Kéo dài hơn (lặp tính toán số ngẫu nhiên)
Độ chính xác đối với mẫu loãng Tương đối (do các giả thiết đơn giản hóa của KM) Rất cao (mô phỏng sát thực tế tương tác vật lý)
Nhu cầu thiết bị phụ trợ Không cần thiết bị đo $g$ trước Xác định trực tiếp đồng thời cả $\mu_a, \mu_s$ và $g$
Ứng dụng thực tế chủ yếu Khảo sát lưu động nhanh ngoài hiện trường Phân tích mẫu chính xác tại phòng thí nghiệm

Bảng tổng hợp danh mục dữ liệu thực nghiệm tiêu biểu được công bố trong Chương 4 của luận án:

Số hiệu bảng Nội dung số liệu thực nghiệm Đối tượng hoặc thông số kiểm tra
Bảng 4.1 Thời gian tính toán nhân ma trận $n \times n$ của DSP Khảo sát hiệu năng C so với Assembler trên TMS320C6713
Bảng 4.2 Sai khác biên độ giữa tín hiệu chuẩn lối vào và hiển thị Kiểm tra độ chính xác của khối mạch ADC-DSP
Bảng 4.3 Sai khác tần số giữa tín hiệu chuẩn lối vào và hiển thị Kiểm tra đáp ứng tần số của khối mạch ADC-DSP
Bảng 4.5 Thông số kỹ thuật chip laser SLD102 Kiểm tra nguồn phát laser tại nhiệt độ phòng
Bảng 4.6 Thông số kỹ thuật của đầu thu OPT101 Đáp ứng quang - điện của photodiode lai tích hợp
Bảng 4.7 Số liệu đo $R_d, T_d, T_c$ trên mẫu sữa tiêu chuẩn Mẫu sữa 4% chất béo pha ở các nồng độ khác nhau
Bảng 4.8 Số liệu xác định $\mu_a, \mu_s, g$ theo MC và KM Đối chuẩn kết quả hai phương pháp trên mẫu sữa chuẩn
Bảng 4.9 & 4.10 Phổ huỳnh quang và phổ kích thích huỳnh quang Khảo sát đặc tính quang phổ các loài vi tảo biển Đồ Sơn
Bảng 4.11 & 4.12 $R_d, T_d, T_c$ và tỷ lệ $\mu_s / \mu_a$ theo mật độ Thực nghiệm trên loài tảo Prorocentrum rhathymum
Bảng 4.13 & 4.14 $R_d, T_d, T_c$ và tỷ lệ $\mu_s / \mu_a$ theo mật độ Thực nghiệm trên loài tảo Pseudo-nitzschia
Bảng 4.15 & 4.16 $R_d, T_d, T_c$ và tỷ lệ $\mu_s / \mu_a$ theo mật độ Thực nghiệm trên loài tảo Alexandrium minutum

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

Các đóng góp mới về mặt lý luận và khoa học của luận án bao gồm:

  1. Đã kết hợp thành công mô hình Kubelka-Munk giải tích (thông qua biến đổi Laplace) và mô phỏng số Monte Carlo với kỹ thuật quả cầu tích phân đôi, áp dụng cho đối tượng mẫu chất lỏng sinh học có nồng độ hạt tán xạ thấp.
  2. Xây dựng được giải pháp xác định đồng thời cả ba thông số quang học cốt lõi gồm hệ số hấp thụ ($\mu_a$), hệ số tán xạ ($\mu_s$) và đặc biệt là hệ số bất đối xứng ($g$) ngay trong một lần đo thực nghiệm, khắc phục được hạn chế cố hữu của các phương pháp trước đây vốn phải đo riêng hoặc giả định trước giá trị của $g$.

Các đóng góp mới về mặt công nghệ, kỹ thuật và thực tiễn:

  1. Thiết kế và chế tạo thành công một thiết bị đo hoàn chỉnh hoạt động lưu động, tích hợp đồng bộ các khối quang điện tử (nguồn laser diode SLD102 có mạch ổn định công suất, hệ quang học chuẩn trực chùm tia, hệ hai quả cầu tích phân ghép nối tiếp và ba đầu thu photodiode lai OPT101).
  2. Thiết kế, chế tạo thành công mạch biến đổi ADC chuyên dụng ghép nối với bo mạch DSP TMS320C6713 thông qua giao diện bộ nhớ ngoài EMIF và CPLD; tối ưu hóa phần mềm xử lý tín hiệu số bằng việc khai thác cấu trúc song song VLIW và kỹ thuật đường ống, đáp ứng yêu cầu xử lý khối lượng dữ liệu lớn từ nhiều kênh tín hiệu đồng thời.
  3. Cung cấp bộ cơ sở dữ liệu thực nghiệm về các thông số quang ($\mu_a, \mu_s, g$) và mối tương quan giữa tỷ lệ $\mu_s / \mu_a$ với mật độ tế bào của ba loài vi tảo độc hại phổ biến tại vùng biển Đồ Sơn (Prorocentrum rhathymum, Pseudo-nitzschia, Alexandrium minutum). Thực nghiệm chứng minh hệ thống có khả năng nhận biết và đo đạc chính xác các mẫu thủy vực có mật độ tảo thấp dưới ngưỡng 1000 tế bào/ml.

Về mặt kiến nghị và giải pháp đề xuất, tác giả kiến nghị đưa thiết bị vào ứng dụng tại các trạm quan trắc môi trường ven biển để phục vụ công tác cảnh báo sớm hiện tượng thủy triều đỏ. Đồng thời, đề xuất khả năng mở rộng hệ đo sang các ứng dụng kiểm tra chất lượng nguồn nước sinh hoạt (phát hiện nhanh vi nấm, tạp chất lơ lửng) và trong y sinh học (xét nghiệm máu phát hiện sớm các tế bào lạ).

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

Về mặt hạn chế, luận án thẳng thắn chỉ ra các điểm tồn tại xuất phát từ điều kiện thực nghiệm và phạm vi nghiên cứu:

  • Đối tượng thực nghiệm vi sinh vật mới chỉ dừng lại ở 03 loài tảo độc hại đơn lẻ đã được phân lập và nuôi cấy trong môi trường phòng thí nghiệm; chưa tiến hành khảo sát trên các mẫu nước biển tự nhiên phức hợp có chứa đồng thời nhiều loài vi tảo cùng các chất vẩn vô cơ lơ lửng và chất hữu cơ hòa tan (CDOM).
  • Nguồn sáng laser diode sử dụng bước sóng đơn sắc cố định (SLD102), chưa khảo sát trên dải phổ biến thiên liên tục để xây dựng dấu vết quang phổ đầy đủ cho từng loài vi tảo.
  • Mô hình Kubelka-Munk dựa trên một số giả thiết đơn giản hóa (tán xạ đẳng hướng, bỏ qua hiệu ứng phản xạ biên phức tạp) nên cấu hình đo nhanh có độ chuẩn xác hạn chế hơn cấu hình Monte Carlo khi áp dụng cho các mẫu quá loãng.

Các hướng nghiên cứu phát triển tiếp theo được đề xuất gồm:

  • Mở rộng hệ đo sử dụng nguồn sáng laser đa bước sóng hoặc nguồn sáng khả kiến có khả năng quét bước sóng để phân tích chi tiết hơn các đỉnh hấp thụ đặc trưng của từng nhóm sắc tố vi tảo.
  • Tiến hành các đợt đo đạc thử nghiệm trên các vùng thủy vực tự nhiên ven bờ vào các mùa khác nhau để hiệu chuẩn mô hình trong điều kiện môi trường nước đa thành phần.
  • Tiếp tục tối ưu hóa thuật toán Monte Carlo trên phần cứng DSP hoặc chip chuyên dụng để rút ngắn thời gian xử lý của cấu hình đo chính xác, hướng tới việc thực hiện phép đo Monte Carlo trực tiếp tại hiện trường trong thời gian thực.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị chuyên môn cao, hữu ích đối với các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên thuộc các chuyên ngành Vật lý Vô tuyến, Điện tử y sinh, Quang điện tử, Kỹ thuật đo lường và Xử lý tín hiệu số. Phần thiết kế mạch ADC tích hợp DSP qua cổng EMIF/CPLD và kỹ thuật tối ưu hóa mã lệnh VLIW (Chương 4) là tài liệu tham khảo thực tế cho các bài toán phát triển hệ thống nhúng xử lý tín hiệu đa kênh.
  • Các nhà nghiên cứu trong lĩnh vực Hải dương học, Sinh thái biển và Công nghệ môi trường. Toàn bộ nội dung Chương 2 và Chương 3 cung cấp khung lý thuyết hoàn chỉnh về tương tác photon trong môi trường dung dịch sinh học, phương pháp giải hệ phương trình Kubelka-Munk bằng biến đổi Laplace và thuật toán Monte Carlo cho hệ quả cầu tích phân đôi.
  • Cán bộ kỹ thuật và nhà quản lý tại 4 Trung tâm quan trắc và cảnh báo môi trường quốc gia. Dữ liệu thực nghiệm về các thông số $\mu_a, \mu_s, g$ của các loài vi tảo biển Đồ Sơn (Bảng 4.11 đến 4.16) là nguồn dữ liệu chuẩn hữu ích phục vụ việc đối chuẩn quang học trong quan trắc môi trường biển ven bờ.

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

1. Tại sao hệ đo lại tích hợp đồng thời hai cấu hình đo nhanh (Kubelka-Munk) và đo chính xác (Monte Carlo)?

Hệ đo tích hợp hai cấu hình nhằm đáp ứng các tình huống sử dụng thực tế khác nhau. Cấu hình đo nhanh dựa trên nghiệm giải tích của mô hình Kubelka-Munk hai dòng bức xạ, cho phép xử lý đại số trực tiếp và trả về kết quả gần như tức thời, rất thích hợp khi cần khảo sát diện rộng hoặc cảnh báo sơ bộ tại hiện trường. Ngược lại, cấu hình đo chính xác sử dụng thuật toán mô phỏng Monte Carlo để theo dõi từng hành trình photon và giải bài toán ngược kết hợp số liệu từ quả cầu tích phân đôi; cấu hình này khắc phục được các giả thiết đơn giản hóa của mô hình KM, mang lại độ chính xác cao đối với các mẫu có mật độ tế bào rất thấp trong phòng thí nghiệm.

2. Thiết bị giải quyết việc xác định hệ số bất đối xứng $g$ như thế nào mà không cần dùng máy đo riêng?

Trong các nghiên cứu quang học truyền thống, hệ số bất đối xứng $g$ thường phải được đo trước bằng các thiết bị đo tán xạ góc chuyên dụng hoặc giả định theo hàm pha. Trong luận án này, nhờ thiết kế quang hệ quả cầu tích phân đôi kết hợp đầu thu thứ ba (PD3) đo chùm truyền thẳng chuẩn trực, hệ đo thu nhận đồng thời ba giá trị độc lập: thành phần phản xạ ngược khuếch tán ($R_d$), thành phần truyền qua khuếch tán ($T_d$) và thành phần truyền qua chuẩn trực ($T_c$). Dựa trên ba đại lượng thực nghiệm này, thuật toán xử lý số (KM và MC) giải đồng thời hệ phương trình để trích xuất cả ba thông số $\mu_a, \mu_s$ và $g$ trong một lần đo duy nhất.

3. Những linh kiện điện tử và quang học cốt lõi nào đã được sử dụng để chế tạo thiết bị?

Phần cứng của hệ đo được xây dựng từ các khối linh kiện chuyên dụng:

  • Bộ xử lý tín hiệu số: bo mạch DSK TMS320C6713 (Texas Instruments) nền tảng DSP dấu phẩy động 32-bit.
  • Mạch ghép nối: bo mạch ADC đa kênh tự thiết kế tích hợp với DSP qua giao diện bộ nhớ ngoài EMIF và vi mạch logic CPLD.
  • Nguồn phát quang: chip laser bán dẫn SLD102 tích hợp mạch điều khiển ổn định công suất phát và nhiệt độ.
  • Đầu thu quang học: 03 cảm biến photodiode hybrid tích hợp mạch khuếch đại thuật toán OPT101 có độ nhạy cao.
  • Hệ quang: buồng đo tích hợp hai quả cầu tích phân ghép nối tiếp có mặt trong phủ lớp phản xạ khuếch tán cao.

4. Kết quả khảo sát thực nghiệm trên các mẫu vi tảo tại Đồ Sơn cho thấy điều gì?

Kết quả đo thực nghiệm trên ba loài tảo độc hại thu thập tại Đồ Sơn (Prorocentrum rhathymum, Pseudo-nitzschia, Alexandrium minutum) ở các dải mật độ khác nhau chỉ ra rằng: khi mật độ tế bào thay đổi, các thành phần quang học $R_d, T_d, T_c$ biến thiên có quy luật rõ rệt. Đáng chú ý, tỷ lệ giữa hệ số tán xạ và hệ số hấp thụ ($\mu_s / \mu_a$) biến đổi theo mật độ tế bào vi tảo, tạo thành một chỉ báo quang học tin cậy giúp định lượng mật độ vi tảo ngay cả khi mật độ ở mức thấp dưới ngưỡng 1000 tế bào/ml.

Kết luận

Luận án "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" của nghiên cứu sinh Nguyễn Tuấn Anh đã giải quyết trọn vẹn cả phương diện lý thuyết mô hình hóa và kỹ thuật chế tạo thiết bị đo quang điện tử ứng dụng trong quan trắc môi trường thủy vực. Công trình đã kết hợp thành công mô hình Kubelka-Munk và mô phỏng Monte Carlo với kỹ thuật quả cầu tích phân đôi, cho phép đo đồng thời ba thông số quang học cốt lõi ($\mu_a, \mu_s, g$) của dung dịch có nồng độ tán xạ thấp. Việc tích hợp thành công phần cứng xử lý số DSP TMS320C6713, mạch ADC-CPLD tự chế tạo, chip laser SLD102 và đầu thu OPT101 đã tạo ra một thiết bị đo lưu động hoàn chỉnh có khả năng đo nhanh và đo chính xác. Các kết quả thực nghiệm trên 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 khả thi của giải pháp trong việc phát hiện sớm tảo nở hoa ở mật độ dưới 1000 tế bào/ml, đóng góp một công cụ kỹ thuật giá trị cho hoạt động bảo vệ môi trường biển ven bờ Việt Nam.