Tổng quan nghiên cứu

Trong môi trường tự nhiên và nhân tạo, sự hiện diện của hơn 60 đồng vị phóng xạ tự nhiên như Kali-40 (K-40), Radi-226 (Ra-226), Thori-232 (Th-232) cùng các nguồn phát xạ từ hoạt động y tế, công nghiệp và năng lượng nguyên tử luôn tiềm ẩn nguy cơ tác động trực tiếp đến sức khỏe con người. Tiêu chuẩn an toàn bức xạ quốc tế quy định suất liều giới hạn cho phép đối với công chúng là 1 mSv/năm và đối với nhân viên làm việc trong môi trường bức xạ là 20 mSv/năm (tính trung bình trong 5 năm). Khi cơ thể hấp thụ liều lượng từ 0.5 Sv đã xuất hiện hiện tượng suy giảm bạch cầu lympho, từ 3 Sv gây rụng tóc nghiêm trọng và khoảng 10 Sv dẫn đến nguy cơ tử vong gần như 100%.

Tại Việt Nam, các trạm quan trắc phóng xạ môi trường ban đầu tại Viện Khoa học và Kỹ thuật hạt nhân, Viện Nghiên cứu hạt nhân và Trung tâm Công nghệ xử lý môi trường đã bộc lộ những hạn chế nhất định về tính đồng bộ, độ nhạy và khả năng cảnh báo trực tuyến theo thời gian thực khi so sánh với mạng lưới 37 trạm của Nhật Bản hay 13 trạm của Hàn Quốc. Vấn đề cốt lõi đặt ra là các thiết bị quan trắc ngoài hiện trường thường chịu ảnh hưởng nặng nề bởi sự thay đổi nhiệt độ và nhiễu thống kê, làm suy giảm độ chính xác và tính ổn định.

Mục tiêu của luận văn là nghiên cứu, hoàn thiện giải pháp kỹ thuật nhằm nâng cao độ chính xác và tính ổn định cho thiết bị quan trắc, cảnh báo phóng xạ môi trường sử dụng đầu dò nhấp nháy NaI(Tl). Nghiên cứu được triển khai thực nghiệm tại Viện Khoa học và Kỹ thuật hạt nhân trong giai đoạn 2011 - 2014, tập trung làm chủ phương pháp chuyển phổ thành liều JAERI, thuật toán lọc trơn Savitzky - Golay, khớp Gauss và cơ chế ghim đỉnh tự động. Kết quả nghiên cứu giúp giảm thiểu sai số đo lường xuống dưới 5%, nâng cao độ tin cậy cho mạng lưới cảnh báo bức xạ quốc gia.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên ba trụ cột lý thuyết nền tảng trong vật lý hạt nhân và kỹ thuật đo lường bức xạ:

Thứ nhất, lý thuyết tương tác của bức xạ gamma với vật chất, bao gồm ba hiệu ứng chủ đạo: hiệu ứng quang điện chiếm ưu thế ở vùng năng lượng thấp với tiết diện tỷ lệ với $Z^5 / E^{7/2}$; hiệu ứng tán xạ Compton chiếm ưu thế ở dải năng lượng trung gian với tiết diện tỷ lệ với $Z / E$; và hiệu ứng tạo cặp electron - positron xảy ra khi năng lượng gamma vượt ngưỡng 1.022 MeV với tiết diện tỷ lệ với $Z^2$. Sự chồng chập của ba hiệu ứng này cùng bức xạ tán xạ ngược tạo nên cấu trúc phổ biên độ xung thực tế trong đầu dò.

Thứ hai, cơ chế biến đổi năng lượng bức xạ của chất nhấp nháy vô cơ NaI(Tl). Tinh thể NaI pha tạp khoảng 0.1% Thallium có mật độ khối 3.67 g/cm³, phát ra photon ánh sáng có bước sóng cực đại 4200 Å với thời gian phân rã xung khoảng 230 ns. Hiệu suất nhấp nháy đạt xấp xỉ 12% năng lượng photon tới (cao gấp 2 lần chất nhấp nháy hữu cơ Anthracene), kết hợp cùng ống nhân quang điện (PMT) có hệ số khuếch đại dòng từ $10^6$ đến $10^8$ lần để biến đổi chớp sáng thành tín hiệu điện thế tương ứng.

Thứ ba, mô hình chuyển đổi phổ thành liều của Viện Nghiên cứu Năng lượng Nguyên tử Nhật Bản (JAERI). Khác với phương pháp đếm xung thô truyền thống, phương pháp JAERI sử dụng hàm toán tử chuyển đổi $G(E)$ để tính toán trực tiếp suất liều hấp thụ $X$ thông qua tích phân phổ số đếm $N(E)$ theo biểu thức:

$$X = \int_{E_{min}}^{E_{max}} N(E) G(E) dE = \sum_{I_{min}}^{I_{max}} N(I) G(I)$$

Trong đó, hàm $G(E)$ được khai triển dưới dạng đa thức bậc $K_{max} = 20$ với tham số đặc trưng hình học tinh thể $M = 3$ và hệ số $A(K)$ chuyên biệt cho đầu dò NaI(Tl) kích thước 2.5 cm.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm thu nhận trực tiếp từ hệ phổ kế hạt nhân sử dụng đầu dò NaI(Tl) hình trụ gắn liền ống nhân quang tử, ghép nối với bộ phân tích đa kênh MCA 1024 kênh. Cỡ mẫu nghiên cứu bao gồm tập hợp hơn 50 tệp phổ năng lượng được ghi nhận liên tục trong các điều kiện đo đạc khác nhau, lưu trữ dưới định dạng chuẩn của bộ MCA tại phòng thí nghiệm Trung tâm Vật lý hạt nhân.

Phương pháp chọn mẫu là chọn mẫu thực nghiệm có chủ đích, khảo sát trên các nguồn bức xạ chuẩn như Cesium-137 (Cs-137 có năng lượng đỉnh 0.662 MeV), Coban-60 (Co-60), Kali-40 (K-40 có năng lượng 1.46 MeV) trong tự nhiên và phông phóng xạ môi trường thực tế tại Hà Nội. Đồng thời, mẫu đo được kiểm thử trong buồng biến đổi nhiệt độ từ 20°C đến 50°C.

Lý do lựa chọn các phương pháp phân tích:

  1. Thuật toán làm trơn Savitzky - Golay bậc 2 qua cửa sổ 5 điểm: Lựa chọn thuật toán này thay vì biến đổi Fourier phức tạp nhằm đảm bảo hiệu năng tính toán nhanh, triệt tiêu hiệu quả thăng giáng ngẫu nhiên Poisson mà không làm méo dạng hoặc dịch chuyển tâm đỉnh phổ, rất phù hợp khi nhúng vào vi điều khiển.
  2. Phương pháp khớp hàm Gauss (Gaussian Fitting): Giúp xác định chính xác tâm đỉnh và bề rộng nửa chiều cao (FWHM) theo công thức giải tích, loại bỏ hoàn toàn sai số chủ quan của người quan sát.
  3. Thuật toán ghim đỉnh K-40 tự nhiên kết hợp bù nhiệt vi sai: Tận dụng đồng vị K-40 luôn có sẵn trong môi trường để làm mốc chuẩn nội, giúp hệ thống tự căn chỉnh độ trôi kênh mà không cần đưa thêm nguồn phóng xạ nhân tạo vào thiết bị.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích thực nghiệm trên hệ đo và phần mềm tự động hóa mang lại ba phát hiện kỹ thuật quan trọng:

Thứ nhất, việc áp dụng bộ lọc Savitzky - Golay 5 điểm kết hợp khớp phân bố Gauss giúp nâng cao rõ rệt độ nét của phổ gamma. Thăng giáng phông ngẫu nhiên giảm trên 80%, cho phép phát hiện chính xác đỉnh quang điện của Cs-137 tại năng lượng 0.662 MeV với độ phân giải năng lượng đạt từ 7.2% đến 8.5%, khắc phục hoàn toàn hiện tượng đỉnh giả xuất hiện do nhiễu đường truyền.

Thứ hai, phương pháp chuyển phổ thành liều JAERI với hàm $G(E)$ đa thức 20 bậc đem lại độ chính xác vượt trội khi so sánh với phương pháp đếm xung truyền thống. Kết quả đo suất liều tại các vị trí có mức bức xạ chuẩn khác nhau cho thấy độ sai lệch tương đối giữa giá trị đo thực nghiệm và giá trị chuẩn luôn duy trì ở mức dưới 4.5%, trong khi phương pháp đếm mật độ xung thông thường có sai số lên đến 20% - 30% ở dải liều thấp dưới 0.1 µSv/h.

Thứ ba, nghiên cứu đã lượng hóa được quy luật trôi kênh của đầu dò theo nhiệt độ môi trường. Khi nhiệt độ tăng từ 20°C lên 26°C, vị trí đỉnh năng lượng của K-40 bị dịch chuyển từ 3 đến 5 kênh trên thang đo 1024 kênh do sự thay đổi độ nhạy của màng quang catot và hệ số nhân điện tử dynode trong ống PMT. Khi áp dụng thuật toán bù nhiệt tuyến tính và ghim đỉnh K-40, độ trôi kênh được kiểm soát chặt chẽ trong phạm vi dưới 0.5 kênh trên toàn dải nhiệt độ từ 20°C đến 50°C, nâng cao độ ổn định tổng thể của thiết bị thêm hơn 90%.

Thảo luận kết quả

Các kết quả thực nghiệm trong luận văn được mô hình hóa và trình bày trực quan thông qua đồ thị phổ biên độ xung ở thang đo tuyến tính và thang logarit, biểu đồ đường cong phân bố Gauss sau khi khớp số liệu, cùng bảng so sánh chi tiết giữa suất liều đo được và suất liều quy chuẩn tại 5 dải giá trị từ 0.08 µSv/h đến 10 µSv/h. Đồ thị hồi quy thực nghiệm cũng thể hiện rõ nét mối tương quan giữa nhiệt độ đo và vị trí kênh đỉnh K-40.

Nguyên nhân chính dẫn đến sự cải thiện vượt bậc về độ chính xác là do hàm toán tử $G(E)$ đã tính đến trọng số năng lượng của từng photon gamma đi vào tinh thể NaI(Tl), thay vì coi mọi xung điện sinh ra đều mang mức năng lượng như nhau. Về tính ổn định, việc sử dụng chính đỉnh K-40 tự nhiên (chu kỳ bán rã 1.25 tỷ năm) làm mốc tham chiếu nội tại là một giải pháp thiết kế thông minh. Giải pháp này giúp thiết bị tự thích ứng với sự thay đổi thời tiết giữa ngày và đêm ngoài trời mà không đòi hỏi hệ thống điều hòa nhiệt độ cồng kềnh hay các nguồn phóng xạ chuẩn nhân tạo đi kèm gây phức tạp về an toàn bức xạ.

Đề xuất và khuyến nghị

Nhằm hoàn thiện và sớm đưa thiết bị quan trắc, cảnh báo phóng xạ vào ứng dụng thực tiễn rộng rãi, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Tối ưu hóa giải thuật nhúng trên vi điều khiển 32-bit: Chuyển đổi mã nguồn xử lý từ phần mềm tính toán sang ngôn ngữ C/C++ nhúng trực tiếp trên chip xử lý ARM, hoàn thiện bảng tra cứu hàm $G(E)$ rút gọn để thời gian tính toán chu kỳ phổ giảm xuống dưới 100 ms, mục tiêu hoàn thành trong 6 tháng do nhóm kỹ sư phần cứng hạt nhân phụ trách.
  2. Tích hợp cảm biến nhiệt độ kỹ thuật số độ nhạy cao: Bố trí cảm biến nhiệt tiếp xúc trực tiếp với vỏ kim loại bọc tinh thể NaI(Tl) với độ chính xác ±0.1°C, liên tục cấp số liệu nhiệt độ thực cho vi xử lý để tự động hiệu chỉnh hệ số khuếch đại điện áp cao (HV) nuôi ống PMT, triển khai thử nghiệm trong 12 tháng tại các trạm đo mẫu.
  3. Mở rộng mạng lưới trạm quan trắc tự động quốc gia: Đề xuất các cơ quan quản lý đầu tư nâng cấp từ 3 trạm hiện tại lên tối thiểu 15 - 20 trạm cảnh báo trực tuyến tại các khu vực trọng điểm biên giới, đô thị lớn và cơ sở hạt nhân trong lộ trình 3 năm tới, do Cục An toàn Bức xạ và Hạt nhân chủ trì.
  4. Xây dựng quy chuẩn hiệu chuẩn định kỳ: Thiết lập quy trình kiểm định sai số suất liều định kỳ 12 tháng/lần bằng các nguồn chuẩn Cs-137 và Co-60 đạt chuẩn quốc tế, đảm bảo độ không đảm bảo đo mở rộng luôn duy trì dưới 5% trong suốt vòng đời thiết bị.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng chính:

  1. Kỹ sư thiết kế và phát triển thiết bị đo lường hạt nhân: Nắm bắt chi tiết nguyên lý cấu tạo đầu dò NaI(Tl), mạch ghép nối ống nhân quang điện, phương pháp tính toán hàm $G(E)$ và kỹ thuật xử lý tín hiệu số trên phổ kế đa kênh MCA.
  2. Cán bộ quản lý và kiểm soát an toàn bức xạ môi trường: Tham khảo phương pháp đánh giá suất liều hấp thụ chính xác, từ đó xây dựng các kịch bản ứng phó sự cố rò rỉ phóng xạ và thiết lập ngưỡng cảnh báo sớm tại địa phương.
  3. Học viên cao học và nghiên cứu sinh ngành Vật lý hạt nhân, Vật lý nguyên tử: Sử dụng làm tài liệu nghiên cứu chuyên khảo về tương tác bức xạ gamma với vật chất, giải thuật lọc trơn Savitzky - Golay và phương pháp fit Gauss trong phân tích phổ hạt nhân.
  4. Đơn vị kỹ thuật thuộc lực lượng phòng hóa và ứng phó sự cố khẩn cấp: Khai thác phương án chế tạo các trạm quan trắc bức xạ tự động hóa, nhỏ gọn, hoạt động bền bỉ trong điều kiện dã ngoại khắc nghiệt.

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

Tại sao đầu dò nhấp nháy NaI(Tl) lại được ưu tiên sử dụng trong quan trắc môi trường hơn đầu dò bán dẫn HPGe?

Đầu dò nhấp nháy NaI(Tl) có hiệu suất ghi nhận tia gamma cao khoảng 12%, hoạt động ổn định ở nhiệt độ phòng và chi phí chế tạo hợp lý. Trong khi đó, đầu dò bán dẫn Ge siêu tinh khiết (HPGe) tuy có độ phân giải năng lượng rất cao nhưng bắt buộc phải làm lạnh liên tục bằng nitơ lỏng ở -196°C, gây khó khăn lớn cho việc vận hành tự động ngoài hiện trường.

Phương pháp chuyển phổ thành liều JAERI ưu việt hơn phương pháp đếm xung truyền thống như thế nào?

Phương pháp đếm xung truyền thống chỉ ghi nhận tổng số lượng xung trong một đơn vị thời gian mà bỏ qua năng lượng của từng hạt, dẫn đến sai số lớn khi môi trường có bức xạ hỗn hợp. Phương pháp JAERI áp dụng toán tử $G(E)$ nhân trực tiếp với số đếm từng kênh năng lượng, giúp xác định chính xác suất liều tương đương với sai số chỉ dưới 4.5%.

Vì sao đồng vị K-40 lại được lựa chọn để ghim đỉnh và bù nhiệt độ cho thiết bị?

Đồng vị phóng xạ Kali-40 (K-40) phát ra tia gamma đơn năng lượng cao 1.46 MeV và luôn tồn tại ổn định trong môi trường đất đá với chu kỳ bán rã 1.25 tỷ năm. Nhờ vậy, thiết bị có thể sử dụng bức xạ tự nhiên này làm tín hiệu chuẩn nội để tự động căn chỉnh trôi kênh mà không cần gắn thêm nguồn phóng xạ nhân tạo nguy hiểm.

Thuật toán Savitzky - Golay làm trơn phổ gamma dựa trên nguyên lý nào?

Thuật toán làm trơn Savitzky - Golay sử dụng phương pháp bình phương tối thiểu để khớp một đa thức bậc 2 qua cửa sổ trượt gồm 5 điểm dữ liệu phổ liên tiếp. Quá trình này triệt tiêu hiệu quả các dao động thăng giáng ngẫu nhiên của nhiễu mà vẫn bảo toàn trung thực diện tích, độ cao và vị trí đỉnh năng lượng của phổ thực.

Nhiệt độ môi trường gây ảnh hưởng tiêu cực ra sao đến kết quả đo của đầu dò NaI(Tl)?

Sự biến thiên nhiệt độ làm thay đổi hiệu suất phát huỳnh quang của tinh thể NaI(Tl) và hệ số nhân điện tử bên trong các dynode của ống PMT. Thực nghiệm cho thấy khi nhiệt độ tăng từ 20°C lên 26°C, đỉnh phổ có thể bị trôi từ 3 đến 5 kênh, dẫn đến sai lệch lớn trong việc tính toán suất liều nếu không có thuật toán bù trừ nhiệt.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán nâng cao độ chính xác và tính ổn định cho thiết bị quan trắc, cảnh báo phóng xạ môi trường sử dụng đầu dò NaI(Tl).
  • Tích hợp thành công thuật toán làm trơn phổ Savitzky - Golay bậc 2 và phương pháp khớp phân bố Gauss, giúp xác định chuẩn xác các đỉnh năng lượng gamma.
  • Áp dụng chuẩn xác phương pháp JAERI với hàm chuyển đổi $G(E)$ đa thức bậc 20, hạ thấp sai số đo suất liều xuống dưới 4.5%.
  • Xây dựng cơ chế bù nhiệt vi sai và ghim đỉnh K-40 tự nhiên, kiểm soát độ trôi kênh dưới 0.5 kênh trong dải nhiệt độ 20°C - 50°C.
  • Toàn bộ thuật toán được tự động hóa hoàn toàn trên môi trường tính toán số liệu, sẵn sàng chuyển giao để nhúng vào phần cứng các trạm quan trắc thực tế.

Trong giai đoạn 12 - 24 tháng tới, hướng phát triển trọng tâm là hoàn thiện bo mạch nhúng vi điều khiển tích hợp truyền thông không dây phục vụ mạng lưới quan trắc thông minh. Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm có thể tham khảo chi tiết các thuật toán và hệ số đa thức trong toàn văn luận văn để ứng dụng vào các dự án phát triển thiết bị đo lường bức xạ chuyên dụng.