Tổng quan nghiên cứu

Nghiên cứu kỹ thuật đo liều bức xạ ion hóa đóng vai trò then chốt trong bảo vệ sức khỏe con người và môi trường trước tác động của phóng xạ. Năm 1990, Ủy ban An toàn Bức xạ Quốc tế đã ban hành khuyến cáo ICRP-60 hạ giới hạn liều bức xạ nghề nghiệp từ 50 mSv/năm xuống 20 mSv/năm tính trung bình trong 5 năm đối với nhân viên bức xạ, và giới hạn 1 mSv/năm đối với công chúng. Đồng thời, Ủy ban Đơn vị Bức xạ Quốc tế đưa ra khuyến cáo ICRU-39 quy định các đại lượng đo liều tương đương cá nhân mới gồm Hp(10), Hp(3) và Hs(0.07). Những thay đổi này đặt ra yêu cầu cấp thiết phải nâng cao độ chính xác, độ nhạy của các hệ thống đo liều, đặc biệt ở các dải bức xạ mức thấp.

Mục tiêu trọng tâm của nghiên cứu là xây dựng phương pháp phân tích số đường cong nhiệt phát quang, hoàn thiện quy trình chế tạo và đánh giá liều kế nhiệt phát quang VINATOM-TLD dựa trên nguồn khoáng vật fluorit (CaF2) tự nhiên của Việt Nam. Đề tài tập trung giải quyết bài toán xác định liều tương đương cá nhân Hp(d), đo liều bức xạ neutron nhanh và bức xạ beta trong các trường bức xạ hỗn hợp.

Nghiên cứu được triển khai thực nghiệm tại Trung tâm Kỹ thuật An toàn Bức xạ thuộc Viện Năng lượng Nguyên tử Việt Nam, kết hợp đối sánh quốc tế tại Viện Đo liều Bức xạ thuộc Viện Hàn lâm Khoa học Tiệp Khắc và Viện PNC Tokai Works Nhật Bản. Công trình mang ý nghĩa thực tiễn to lớn khi tạo ra sản phẩm liều kế nội địa đạt độ chính xác cao với sai số dưới 3% trong khoảng liều từ vài mGy đến 10 Gy, giúp giảm thiểu hơn 70% chi phí nhập khẩu thiết bị đo liều từ nước ngoài.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết của hiệu ứng nhiệt phát quang (TLD) dựa trên sự dịch chuyển và tích lũy các hạt tải điện (electron và lỗ trống) tại các trạng thái khuyết tật trong mạng tinh thể chất rắn khi tương tác với bức xạ ion hóa. Luận văn vận dụng ba mô hình động học nền tảng:

  • Mô hình động học bậc một của Randall và Wilkins: Giả định xác suất tái bẫy bằng 0, đường cong nhiệt phát quang có dạng bất đối xứng rõ rệt và nhiệt độ đỉnh phát quang cực đại (Tmax) không phụ thuộc vào nồng độ hạt tải điện ban đầu.
  • Mô hình động học bậc hai của Garlick và Gibson: Xét trường hợp xác suất tái bẫy tương đương xác suất tái hợp phát quang, tạo nên dạng phổ đối xứng hơn và vị trí Tmax dịch chuyển về vùng nhiệt độ thấp khi nồng độ hạt tải tăng.
  • Phương trình động học tổng quát của Chen và May: Sử dụng tham số bậc động học b biến thiên từ 1 đến 2 để mô tả chính xác mọi biến dạng phổ phát quang thực tế.

Bên cạnh đó, nghiên cứu tích hợp các khái niệm chuẩn hóa quốc tế gồm đường cong nhiệt phát quang (Glow Curve), năng lượng kích hoạt bẫy (E), hệ số tần số (S), đại lượng Kerma mô và liều tương đương cá nhân theo khuyến cáo ICRU-39 và ICRU-47.

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

Nguồn dữ liệu thực nghiệm được thu thập từ hơn 120 mẫu đo nhiệt phát quang gồm CaF2 tự nhiên Việt Nam, CaSO4:Dy, CaSO4:Mn, LiF (TLD-100, TLD-600, TLD-700), BeO và Al2O3. Phương pháp chọn mẫu là chọn mẫu thực nghiệm có đối chứng, phân tầng theo các dải năng lượng photon từ 30 keV đến 1,25 MeV của nguồn bức xạ chuẩn Cobalt-60 và Cesium-137.

Phương pháp phân tích số đường cong phát quang (Numerical Glow Curve Deconvolution) được lựa chọn thay thế phương pháp đồ thị truyền thống. Lý do lựa chọn phương pháp số là khả năng tách bạch các đỉnh nhiệt phát quang chồng lấn, loại bỏ hoàn toàn sai số chủ quan của người đọc, rút ngắn 50% thời gian xử lý dữ liệu và hỗ trợ trích xuất chính xác các thông số động học E, S, b. Quá trình nghiên cứu và chuẩn hóa được tiến hành đồng bộ trong giai đoạn từ năm 1994 đến năm 1996 với sự hỗ trợ của hệ thống máy tính tự động hóa.

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

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

Thực nghiệm đo liều và khảo sát đặc trưng vật liệu đã mang lại 4 phát hiện khoa học quan trọng:

  • Liều kế VINATOM-TLD chế tạo từ fluorit thiên nhiên Việt Nam có độ nhạy rất cao với bức xạ photon năng lượng thấp. Tỷ số độ nhạy tại mức năng lượng 30 keV cao gấp 15 lần so với mức năng lượng 1,2 MeV từ nguồn Cobalt-60, thể hiện sự phù hợp tối ưu cho việc đo trực tiếp liều tương đương cá nhân Hp(d).
  • Phát hiện các dị thường quang phổ tại các đỉnh phát quang nhiệt độ thấp (dưới 150°C) của mẫu CaF2 và CaSO4:Dy. Sự suy giảm theo hàm mũ của các đỉnh bẫy nông này cho phép ứng dụng phương pháp tách phổ số để xác định chính xác thời điểm xảy ra sự cố chiếu xạ với sai số thời gian dưới 5%.
  • Việc phối trộn hợp chất giàu hydro như Polyethylene (PE) vào vật liệu CaSO4:Mn, CaSO4:Dy và Al2O3 đã làm tăng độ nhạy đối với neutron nhanh lên gấp 10 đến 20 lần so với liều kế TLD nguyên bản, tạo tiền đề phân tách thành phần Kerma neutron (Kn) và Kerma gamma (Kγ) trong trường bức xạ hỗn hợp n + γ.
  • Khảo sát liều kế hình xuyến UD-100M8 xác định được mức độ suy giảm tán xạ của chùm tia beta chuẩn, duy trì độ ổn định tín hiệu với độ suy giảm dưới 10% trong 24 giờ đầu ở nhiệt độ phòng.

Thảo luận kết quả

Cơ chế vật lý của sự sai khác về độ nhạy năng lượng bắt nguồn từ hiệu ứng quang điện chiếm ưu thế ở dải năng lượng thấp đối với các vật liệu có nguyên tử số hiệu dụng lớn (Zeff của CaF2 là 16,57 so với mô sinh học là 7,4). Dữ liệu thực nghiệm được trình bày trực quan qua 36 đồ thị đường cong nhiệt phát quang biểu diễn cường độ theo nhiệt độ từ 50°C đến 400°C và 26 bảng số liệu so sánh hệ số chuẩn liều.

So sánh với các nghiên cứu quốc tế trên dòng vật liệu thương mại của hãng Harshaw, liều kế VINATOM-TLD đạt độ tuyến tính tương đương trong dải liều rộng từ 10 µGy đến 10 Gy. Kết quả thử nghiệm liên phòng tại khu vực Châu Á - Thái Bình Dương khẳng định thuật toán số xác định năng lượng hiệu dụng và liều tương đương cá nhân của tác giả đạt độ tin cậy cao, độ lệch chuẩn chỉ dao động trong mức 3% đến 6%, hoàn toàn đáp ứng các tiêu chuẩn nghiêm ngặt của IAEA.

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

Dựa trên kết quả nghiên cứu, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

  • Chuẩn hóa quy trình chế tạo hàng loạt liều kế VINATOM-TLD từ fluorit thiên nhiên: Đạt sản lượng mục tiêu trên 10.000 liều kế mỗi năm, kiểm soát độ đồng đều giữa các mẻ sản xuất với độ lệch nhạy dưới 5% trong vòng 12 tháng tới, do Trung tâm Kỹ thuật An toàn Bức xạ chủ trì.
  • Ứng dụng phần mềm tính toán số tự động hóa: Tích hợp thuật toán xác định đại lượng Hp(10), Hp(3) và Hs(0.07) vào hệ thống quản lý dữ liệu bức xạ quốc gia, hướng tới cắt giảm 60% thời gian lập báo cáo định kỳ cho các cơ sở y tế và công nghiệp trong vòng 6 tháng tới.
  • Xây dựng hệ thống chuẩn liều bức xạ hỗn hợp neutron - gamma và bức xạ beta: Thiết lập trường bức xạ chuẩn quốc gia với độ bất định đo dưới 3,5% cho dải liều từ 1 mGy đến 100 mGy, do Viện Năng lượng Nguyên tử Việt Nam phối hợp Tổng cục Tiêu chuẩn Đo lường Chất lượng thực hiện trong lộ trình 2 năm.
  • Ban hành quy chuẩn giám sát an toàn bức xạ cá nhân mới: Áp dụng đầy đủ khuyến cáo ICRP-60 và ICRU-39, kiểm soát nghiêm ngặt mức liều trần 20 mSv/năm cho 100% nhân viên làm việc với bức xạ ion hóa trên toàn quốc, do Bộ Khoa học và Công nghệ ban hành định kỳ hàng năm.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chuyên môn:

  • Cán bộ phụ trách an toàn bức xạ tại các bệnh viện có khoa Ung bướu, Xạ trị, Y học hạt nhân: Nắm bắt phương pháp tính liều tương đương cá nhân Hp(10) và Hp(0.07), ứng dụng kiểm soát an toàn cho hơn 500 nhân viên y tế trực tiếp vận hành máy gia tốc và nguồn phóng xạ.
  • Kỹ sư và nhà nghiên cứu công nghệ vật liệu hạt nhân: Tiếp cận công nghệ chế tạo liều kế TLD từ khoáng sản trong nước, tối ưu hóa quá trình xử lý nhiệt (annealing) và giảm thiểu 70% chi phí mua vật liệu ngoại nhập.
  • Cơ quan quản lý nhà nước về an toàn hạt nhân: Sử dụng các dữ liệu đối chuẩn quốc tế và phân tích đại lượng đo liều mới để hoàn thiện hành lang pháp lý, quy chuẩn đo lường bức xạ tại Việt Nam.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Vật lý hạt nhân, Vật lý y khoa: Khai thác 106 tài liệu tham khảo cùng hệ thống mô hình động học Randall - Wilkins, Garlick - Gibson để phục vụ công tác giảng dạy và phát triển các đề tài mô phỏng chuyên sâu.

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

  • Tại sao cần chuyển đổi cách đo liều sang đại lượng Hp(d) theo ICRU-39? Khuyến cáo ICRU-39 đưa ra đại lượng liều tương đương cá nhân Hp(d) giúp đánh giá trực tiếp liều bức xạ hấp thụ tại các độ sâu mô thực tế (0,07 mm cho da, 3 mm cho thủy tinh thể và 10 mm cho toàn thân), khắc phục nhược điểm của khái niệm liều chiếu cũ và đảm bảo an toàn dưới ngưỡng 20 mSv/năm.

  • Ưu điểm của fluorit (CaF2) tự nhiên Việt Nam trong chế tạo liều kế TLD là gì? Fluorit thiên nhiên Việt Nam có độ nhạy nhiệt phát quang cao, đặc biệt nhạy với bức xạ photon năng lượng thấp (gấp 15 lần ở mức 30 keV so với 1,2 MeV), cấu trúc bẫy ổn định, độ bền cơ học tốt và chi phí khai thác chế tạo rẻ hơn nhiều so với tinh thể nhập ngoại.

  • Làm thế nào để xác định thời điểm xảy ra sự cố bức xạ bằng liều kế TLD? Phương pháp tách phổ số phân tích tốc độ suy giảm tín hiệu của các đỉnh nhiệt độ thấp (dưới 150°C). Do các đỉnh này phân rã nhanh theo thời gian sau khi chiếu xạ, tỷ số diện tích giữa đỉnh nhiệt độ thấp và đỉnh chính ở 220°C cho phép suy ngược chính xác thời điểm xảy ra sự cố trong vòng 24 đến 48 giờ.

  • Polyethylene (PE) có vai trò gì trong việc đo liều neutron nhanh? Polyethylene chứa hàm lượng hydro rất cao. Khi neutron nhanh va chạm với hạt nhân hydro (proton) trong PE sẽ tạo ra các proton giật lùi, các proton này ion hóa vật liệu TLD xung quanh, giúp tăng độ nhạy phát quang của liều kế CaSO4 và Al2O3 lên 10 đến 20 lần đối với neutron.

  • Chế độ xử lý nhiệt (annealing) ảnh hưởng ra sao đến liều kế TLD? Quá trình nung nhiệt ở mức nhiệt độ chuẩn từ 300°C đến 400°C giúp giải phóng toàn bộ electron còn sót lại trong các bẫy sâu, tái lập trạng thái ban đầu của mạng tinh thể, đảm bảo độ nhạy không đổi và loại bỏ sai số tích lũy cho các chu kỳ đo tiếp theo.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết động học nhiệt phát quang bậc một, bậc hai và phương trình động học tổng quát, phục vụ phân tích chính xác cấu trúc bẫy năng lượng vật liệu.
  • Chế tạo thành công liều kế VINATOM-TLD từ fluorit thiên nhiên Việt Nam đạt tiêu chuẩn đo liều tương đương cá nhân với sai số đo kiểm soát dưới 3%.
  • Phát triển thành công thuật toán số hóa tự động giải đoán đường cong nhiệt phát quang và định lượng liều bức xạ photon theo các tiêu chuẩn quốc tế ICRU-39 và ICRP-60.
  • Đề xuất giải pháp pha trộn Polyethylene nâng cao độ nhạy đo neutron nhanh và ứng dụng hiệu quả liều kế UD-100M8 trong đo liều beta.
  • Mở ra hướng ứng dụng độc đáo trong việc xác định thời điểm xảy ra sự cố bức xạ thông qua phân tích các đỉnh nhiệt phát quang nhiệt độ thấp.

Trong giai đoạn 12 đến 24 tháng tới, hệ thống đo liều số hóa cần được nhân rộng tại tất cả các cơ sở ứng dụng bức xạ trên toàn quốc. Các đơn vị y tế và công nghiệp hạt nhân hãy chủ động liên hệ với Trung tâm Kỹ thuật An toàn Bức xạ để chuyển giao công nghệ liều kế VINATOM-TLD, nâng cao tiêu chuẩn an toàn bức xạ và tối ưu hóa chi phí vận hành.