Tổng quan nghiên cứu

Trong kỹ thuật hạt nhân và nghiên cứu an toàn bức xạ toàn cầu, việc kiểm soát hơn 440 lò phản ứng hạt nhân thương mại cùng hệ thống phân tích đồng vị đòi hỏi độ chính xác cao khi xác định độ giàu của $^{235}\text{U}$ (chiếm khoảng 0,72% trong tự nhiên) cũng như tuổi của nhiên liệu hạt nhân. Phương pháp phổ kế gamma thông thường yêu cầu xây dựng đường cong hiệu suất ghi tuyệt đối dựa trên các mẫu chuẩn vật lý có cùng hình học, kích thước và mật độ với mẫu phân tích. Tuy nhiên, trong thực tế kiểm tra không phá hủy, các mẫu đo thường có hình dạng bất kỳ, kích thước không đồng nhất hoặc đặt trong các vỏ bọc che chắn phức tạp, dẫn đến sai số thực nghiệm có thể vượt mức 15% đến 30% do hiệu ứng tự hấp thụ và biến thiên góc khối.

Nhằm giải quyết triệt để rào cản này, nghiên cứu tập trung đánh giá bằng thực nghiệm tính ưu việt của phương pháp chuẩn nội hiệu suất ghi (Internal Calibration of Detection Efficiency). Mục tiêu cụ thể là chứng minh khả năng xác định chính xác tỉ số hoạt độ phóng xạ của các đồng vị mà không cần sử dụng mẫu chuẩn và hoàn toàn không phụ thuộc vào cấu hình hình học đo; đồng thời ứng dụng phương pháp để đánh giá trạng thái cân bằng phóng xạ trong dãy $^{238}\text{U}$ và xác định tỉ số hoạt độ giữa $^{235}\text{U}$ và $^{238}\text{U}$.

Nghiên cứu được triển khai thực nghiệm tại Bộ môn Vật lý hạt nhân, Khoa Vật lý, Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội vào năm 2015. Ý nghĩa của đề tài thể hiện ở việc cắt giảm 100% sự phụ thuộc vào nguồn chuẩn ngoại lai đắt đỏ, rút ngắn thời gian chuẩn bị đo kiểm từ nhiều tuần xuống vài giờ, đồng thời duy trì hệ số tương quan làm khớp hàm $R^2 > 0,998$, mang lại giá trị ứng dụng cao cho công tác giám sát nhiên liệu hạt nhân và quan trắc phóng xạ môi trường.

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 trên nền tảng định luật phân rã phóng xạ vi phân và lý thuyết chuỗi phân rã phóng xạ liên tiếp (phương trình Bateman). Hiện tượng cân bằng phóng xạ thế kỷ (Secular Equilibrium) được áp dụng cho các dãy phóng xạ tự nhiên gồm $^{238}\text{U}$ ($T_{1/2} = 4,47 \cdot 10^9$ năm), $^{235}\text{U}$ ($T_{1/2} = 7,04 \cdot 10^8$ năm) và $^{232}\text{Th}$ ($T_{1/2} = 4,39 \cdot 10^9$ năm), nơi chu kỳ bán rã của hạt nhân mẹ đầu dãy lớn hơn rất nhiều so với các hạt nhân con cháu, dẫn đến hoạt độ phóng xạ giữa các đồng vị trong chuỗi đạt trạng thái bằng nhau sau một khoảng thời gian đủ lớn.

Mô hình phương pháp chuẩn nội hiệu suất ghi mở rộng phương trình tỷ số tốc độ đếm sang hàm hiệu suất tương đối $f(E)$. Theo mô hình này, tỷ số tốc độ đếm chia cho hệ số phân nhánh $n/Br$ của một đồng vị được so sánh trực tiếp với giá trị nội suy $f(E)$ tại cùng năng lượng được xây dựng từ các vạch gamma phát ra bởi đồng vị thứ hai trong cùng một mẫu. Các khái niệm cốt lõi bao gồm: hoạt độ phóng xạ ($H$, tính bằng Becquerel - Bq), hệ số phân nhánh bức xạ gamma ($Br$), diện tích đỉnh hấp thụ toàn phần, độ phân giải năng lượng (FWHM) và hàm chuẩn nội hiệu suất ghi.

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

Nghiên cứu sử dụng cỡ mẫu gồm 02 mẫu nguồn phóng xạ điển hình: nguồn chuẩn dãy Thori ký hiệu TS5 và nguồn quặng Urani tự nhiên ký hiệu US2. Phương pháp chọn mẫu có chủ đích được áp dụng nhằm chọn các mẫu phát ra tập hợp vạch phổ gamma phong phú trải rộng từ vùng năng lượng thấp đến vùng năng lượng cao, đáp ứng điều kiện làm khớp đường cong chuẩn nội.

Phương pháp phân tích quang phổ gamma sử dụng đầu dò bán dẫn Gecmani siêu tinh khiết dải rộng (Broad Energy Germanium - BEGe), model Canberra BE5030 có diện tích hiệu dụng $5000\text{ mm}^2$, dày 31 mm, hoạt động trong dải năng lượng từ 3 keV đến 3 MeV. Hệ đo được đặt trong buồng chì che chắn phông thấp dày 10 cm, khối lượng 950 kg, kết hợp khối tiền khuếch đại Canberra 2002C, khối khuếch đại phổ Canberra 2026, cao áp Canberra 3106D và bộ phân tích đa kênh MCA 16.000 kênh.

Quá trình thu nhận và xử lý số liệu diễn ra liên tục với thời gian đo thực nghiệm đạt $69.270\text{ giây}$ (đối với mẫu TS5 cấu hình không bọc chì), $83.181\text{ giây}$ (đối với mẫu TS5 bọc chì 1,5 mm) và $16.058\text{ giây}$ (đối với mẫu US2). Các phần mềm chuyên dụng gồm GammaVision 6.01, FitzPeaks 3.63, Maestro 2.2 và Origin 8.0 được ứng dụng để làm khớp đỉnh Gauss phi tuyến, hiệu chỉnh thời gian chết, hiệu ứng chồng chập xung và hiệu ứng cộng đỉnh.

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

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

Thực nghiệm đã chứng minh tính độc lập hoàn toàn của phương pháp chuẩn nội đối với hình học đo và sự suy giảm bức xạ qua lớp hấp thụ. Khi khảo sát mẫu TS5 ở hai cấu hình đo khác biệt: cấu hình 1 (nguồn không bọc chì đặt song song mặt đầu dò) và cấu hình 2 (nguồn bọc chì dày 1,5 mm), tốc độ đếm và đại lượng $n/Br$ tại từng đỉnh năng lượng sụt giảm rõ rệt. Cụ thể, tại đỉnh 338,3 keV của $^{228}\text{Ac}$, đại lượng $n/Br$ giảm 50,8% (từ 187,2 xuống 92,1); tại đỉnh 583,2 keV của $^{208}\text{Tl}$, $n/Br$ giảm 24,1% (từ 116,1 xuống 88,1). Dù dữ liệu đếm thô biến thiên mạnh, tỉ số hoạt độ $A(^{208}\text{Tl})/A(^{228}\text{Ac})$ xác định qua hàm chuẩn nội ở cả hai cấu hình đều đạt kết quả đồng nhất xấp xỉ nhau với độ tin cậy hàm khớp cực cao ($R^2 = 0,9999$ cho cấu hình song song và $R^2 = 0,99838$ cho cấu hình bọc chì).

+-----------------------------------------------------------------------------+
| SO SÁNH BIẾN THIÊN n/Br VÀ TỈ SỐ HOẠT ĐỘ GIỮA 2 CẤU HÌNH ĐO (MẪU TS5)       |
+-----------------------------------------------------------------------------+
| Đỉnh năng lượng (keV) | Cấu hình không chì (n/Br) | Cấu hình bọc chì (n/Br) |
+-----------------------+---------------------------+-------------------------+
| 338,3 keV (228Ac)     | 187,2                     | 92,1 (giảm 50,8%)       |
| 583,2 keV (208Tl)     | 116,1                     | 88,1 (giảm 24,1%)       |
+-----------------------+---------------------------+-------------------------+
| KẾT QUẢ TỈ SỐ HOẠT ĐỘ | A(208Tl) / A(228Ac)       | A(208Tl) / A(228Ac)     |
|                       | -> Đồng nhất              | -> Đồng nhất            |
+-----------------------------------------------------------------------------+

Đánh giá trạng thái cân bằng phóng xạ chuỗi $^{238}\text{U}$ trên mẫu US2 qua 7 đỉnh gamma đặc trưng của $^{214}\text{Bi}$ (từ 609,31 keV đến 1764,49 keV) cho thấy hàm chuẩn nội đa thức bậc hai đạt hệ số khớp $R = 0,99654$. Tỉ số hoạt độ tính được $A(^{214}\text{Pb})/A(^{214}\text{Bi}) \approx 1,00$ và $A(^{234m}\text{Pa})/A(^{214}\text{Bi}) \approx 1,00$ khẳng định chuỗi $^{238}\text{U}$ trong mẫu đã đạt trạng thái cân bằng thế kỷ hoàn toàn.

Phương pháp đã giải quyết thành công bài toán tách đỉnh chập phức tạp tại vùng năng lượng 186 keV (gồm đỉnh 186,21 keV của $^{226}\text{Ra}$ với $Br = 3,6%$ và đỉnh 185,75 keV của $^{235}\text{U}$ với $Br = 57,2%$) trên tổng số đếm tích lũy $205.531\text{ xung}$. Tỉ số hoạt độ thực nghiệm giữa $^{235}\text{U}$ và $^{238}\text{U}$ xác định được là $0,046$, trùng khớp chính xác 100% với giá trị tính toán lý thuyết ($0,046$).

Thảo luận kết quả

Cơ chế mang lại độ chính xác vượt trội của phương pháp chuẩn nội nằm ở việc triệt tiêu đại số thừa số góc khối $\Omega$ và hệ số tự hấp thụ của mẫu trong phương trình tỷ số hoạt độ. Khi lấy tỷ số giữa tốc độ đếm thực nghiệm và giá trị nội suy từ đường cong hiệu suất ghi $f(E)$ tại cùng một giá trị năng lượng, các yếu tố hình học và hiệu ứng suy giảm vật chất tác động đồng thời lên hai vạch bức xạ đều bị loại trừ.

Dữ liệu thực nghiệm khi biểu diễn qua biểu đồ phân bố đường cong chuẩn nội $f(E)$ theo năng lượng $E$ thể hiện quy luật suy giảm đơn điệu và liên tục. Bảng đối chiếu số liệu giữa hai cấu hình bọc chì và không bọc chì minh chứng rõ ràng: mặc dù lớp chì 1,5 mm làm suy giảm mạnh các bức xạ năng lượng thấp, đường cong hiệu suất tương đối vẫn tái lập chính xác tỉ số hoạt độ của các hạt nhân phóng xạ. So với các phương pháp mô phỏng Monte Carlo vốn dễ mắc sai số hệ thống từ 5% đến 10% do thiếu thông số chính xác về vùng chết detector, phương pháp chuẩn nội thể hiện độ tin cậy và tính khả thi vượt trội.

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

  1. Chuẩn hóa quy trình phân tích chuẩn nội trong kiểm định an toàn hạt nhân: Bộ Khoa học và Công nghệ phối hợp với Viện Năng lượng Nguyên tử Việt Nam xây dựng tiêu chuẩn quốc gia (TCVN) về quy trình kiểm tra không phá hủy thanh nhiên liệu hạt nhân bằng phương pháp chuẩn nội, nhằm mục tiêu giảm 30% thời gian phân tích mẫu tại các cơ sở bức xạ trước năm 2028.
  2. Nâng cấp trang thiết bị và tích hợp phần mềm xử lý phổ tự động: Các viện nghiên cứu chuyên ngành và phòng thí nghiệm trọng điểm đầu tư hệ phổ kế bán dẫn BEGe thế hệ mới với độ phân giải năng lượng đạt FWHM dưới 1,8 keV tại đỉnh 1332 keV của $^{60}\text{Co}$, đồng thời lập trình tự động hóa các thuật toán tách đỉnh chập và khớp hàm chuẩn nội đạt độ tin cậy trên 99% trong giai đoạn 2026 – 2027.
  3. Mở rộng ứng dụng trong quan trắc môi trường và công nghiệp khai khoáng: Cục An toàn Bức xạ và Hạt nhân triển khai áp dụng phương pháp chuẩn nội để đánh giá nhanh tỉ số hoạt độ phóng xạ tự nhiên trong quặng đất hiếm, cát đen và vật liệu xây dựng, kiểm soát sai số đo lường dưới 5% bắt đầu từ quý 2/2027.
  4. Đào tạo và chuyển giao kỹ thuật xử lý phổ hạt nhân nâng cao: Trường Đại học Khoa học Tự nhiên phối hợp cùng các trung tâm nghiên cứu hạt nhân tổ chức các khóa đào tạo chuyên sâu về hiệu chỉnh hiệu ứng cộng đỉnh, xử lý thời gian chết và khớp hàm phi tuyến cho ít nhất 150 cán bộ kỹ thuật hạt nhân trên cả nước đến năm 2029.

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

  • Cán bộ nghiên cứu vật lý hạt nhân và kỹ thuật đo lường bức xạ: Khai thác thuật toán xây dựng hàm chuẩn nội $f(E)$, kỹ thuật làm khớp hàm phi tuyến và các giải pháp hiệu chỉnh hiệu ứng cộng đỉnh nhằm nâng cao độ chính xác phép đo phổ gamma.
  • Kỹ sư giám sát an toàn và quản lý nhiên liệu hạt nhân: Ứng dụng quy trình xác định tỉ số hoạt độ $^{235}\text{U}/^{238}\text{U}$ để kiểm tra độ giàu và tuổi nhiên liệu hạt nhân trong các thanh nhiên liệu có vỏ bọc kim loại mà không cần phá hủy mẫu.
  • Chuyên viên quan trắc phóng xạ môi trường và địa chất khoáng sản: Sử dụng phương pháp để đánh giá nhanh trạng thái cân bằng phóng xạ của các dãy $^{238}\text{U}$ và $^{232}\text{Th}$ trong mẫu đất, đá, nước ngầm có hình học phức tạp mà không tốn chi phí chế tạo mẫu chuẩn.
  • Giảng viên, học viên cao học và sinh viên ngành Vật lý: Sử dụng tài liệu như một cẩm nang thực nghiệm mẫu mực về kỹ thuật vận hành phổ kế bán dẫn BEGe Canberra, xử lý phổ đa kênh bằng GammaVision, Maestro và Origin.

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

Phương pháp chuẩn nội hiệu suất ghi có ưu điểm vượt trội gì so với phương pháp chuẩn tuyệt đối?

Phương pháp chuẩn tuyệt đối bắt buộc phải sử dụng mẫu chuẩn có hình học, mật độ và thành phần hóa học đồng nhất tuyệt đối với mẫu đo để xác định hiệu suất ghi $\varepsilon(E)$. Ngược lại, phương pháp chuẩn nội tận dụng chính các vạch bức xạ gamma của đồng vị trong mẫu để xây dựng hàm hiệu suất tương đối $f(E)$, loại bỏ hoàn toàn yêu cầu về mẫu chuẩn và cho kết quả chính xác trên mọi hình dạng mẫu đo.

Làm thế nào để chọn các vạch gamma phù hợp nhằm xây dựng đường cong chuẩn nội?

Để xây dựng đường cong $f(E)$ tin cậy, cần chọn các vạch gamma phát ra từ cùng một đồng vị (hoặc các đồng vị cân bằng phóng xạ) có dải năng lượng trải rộng, cường độ phát xạ $Br$ lớn và đã biết chính xác, đồng thời các đỉnh phổ phải độc lập, không bị chồng chập bởi bức xạ của các nguyên tố khác, đảm bảo hệ số tương quan làm khớp $R^2 \ge 0,99$.

Tại sao cần thực hiện hiệu chỉnh thời gian chết và hiệu ứng cộng đỉnh trong phép đo phổ gamma?

Thời gian chết làm mất số đếm khi bộ biến đổi tương tự số (ADC) bận xử lý xung, gây sai lệch tốc độ đếm thực tế ở các mẫu có hoạt độ cao. Hiệu ứng cộng đỉnh xảy ra khi hai photon gamma phát ra đồng thời bị đầu dò ghi nhận cùng lúc, làm suy giảm diện tích đỉnh thành phần và tạo đỉnh giả, do đó bắt buộc phải hiệu chỉnh để bảo đảm tính chính xác của tỉ số hoạt độ.

Vì sao có thể xác định hoạt độ của $^{238}\text{U}$ thông qua đỉnh gamma 1001,03 keV của $^{234m}\text{Pa}$?

Đồng vị $^{238}\text{U}$ phát hạt alpha và có chu kỳ bán rã rất dài (khoảng 4,47 tỉ năm) nên hầu như không phát tia gamma đo được trực tiếp. Tuy nhiên, các đồng vị con cháu đời đầu gồm $^{234}\text{Th}$ và $^{234m}\text{Pa}$ ($T_{1/2} = 1,175$ phút) nhanh chóng đạt trạng thái cân bằng thế kỷ với $^{238}\text{U}$, cho phép tính toán trực tiếp hoạt độ của $^{238}\text{U}$ qua vạch gamma 1001,03 keV của $^{234m}\text{Pa}$.

Phương pháp xử lý bài toán tách đỉnh chập tại vùng 186 keV giữa $^{226}\text{Ra}$ và $^{235}\text{U}$ như thế nào?

Vùng 186 keV chứa hai đỉnh năng lượng nằm sát nhau là 186,21 keV ($^{226}\text{Ra}$) và 185,75 keV ($^{235}\text{U}$). Phương pháp chuẩn nội xác định hoạt độ $^{226}\text{Ra}$ thông qua các vạch cân bằng của $^{214}\text{Bi}$, sau đó nội suy giá trị hàm $f(E)$ tại 186,21 keV để tính chính xác phần đóng góp của $^{226}\text{Ra}$, từ đó trừ ra và xác định chuẩn xác hoạt độ riêng biệt của $^{235}\text{U}$.

Kết luận

  • Đã kiểm chứng thực nghiệm xuất sắc tính ưu việt của phương pháp chuẩn nội hiệu suất ghi, khẳng định kết quả đo tỉ số hoạt độ hoàn toàn không phụ thuộc vào hình học mẫu và lớp chắn hấp thụ bên ngoài.
  • Khảo sát thành công mẫu nguồn TS5 ở hai cấu hình đo bọc chì 1,5 mm và không bọc chì, cho tỉ số hoạt độ $A(^{208}\text{Tl})/A(^{228}\text{Ac})$ đồng nhất với hệ số tương quan hàm khớp đạt mức $R^2 \ge 0,998$.
  • Đánh giá chuẩn xác trạng thái cân bằng phóng xạ thế kỷ của chuỗi phân rã $^{238}\text{U}$ trên mẫu US2 với các tỉ số hoạt độ $A(^{214}\text{Pb})/A(^{214}\text{Bi})$ và $A(^{234m}\text{Pa})/A(^{214}\text{Bi})$ đều tiệm cận 1,00.
  • Tách thành công vùng đỉnh chập 186 keV, xác định tỉ số hoạt độ $^{235}\text{U}/^{238}\text{U}$ thực nghiệm đạt $0,046$, trùng khớp tuyệt đối 100% với giá trị tính toán lý thuyết hạt nhân.
  • Khẳng định giá trị thực tiễn to lớn của phương pháp trong công tác kiểm soát độ giàu nhiên liệu hạt nhân và phân tích mẫu môi trường phức tạp tại Việt Nam.

Đóng góp quan trọng nhất của luận văn là thiết lập một quy trình thực nghiệm chuẩn tắc, tin cậy cao, giúp các nhà khoa học loại bỏ hoàn toàn sự lệ thuộc vào các mẫu chuẩn đắt đỏ. Kế hoạch nghiên cứu tiếp theo đến năm 2027 cần hướng tới việc tự động hóa thuật toán chuẩn nội trên các nền tảng số hóa để tích hợp trực tiếp vào các trạm quan trắc bức xạ tự động. Hãy liên hệ các trung tâm nghiên cứu chuyên ngành hoặc tham khảo trực tiếp toàn văn tài liệu tại Thư viện Đại học Quốc gia Hà Nội để áp dụng phương pháp khoa học tiên tiến này vào các dự án nghiên cứu thực tế.