Tổng quan nghiên cứu

Phóng xạ tự nhiên hiện diện phổ biến trong môi trường sống với suất liều bức xạ vũ trụ trung bình khoảng 0,04 µGy/h trên bề mặt Trái Đất và liều hấp thụ hàng năm của con người đạt khoảng 0,27 mSv. Tỷ lệ này tăng gấp 2 lần sau mỗi 2.000m độ cao. Bên cạnh đó, nồng độ khí phóng xạ Radon-222 trong không gian kín có thể biến thiên mạnh từ 20 Bq/m³ đến hơn 10.000 Bq/m³, trở thành một trong những tác nhân hàng đầu gây ung thư phổi. Trong bối cảnh năng lượng hạt nhân cung cấp khoảng 16% tổng sản lượng điện toàn cầu với hơn 441 nhà máy vận hành, việc quan trắc và định lượng chính xác hoạt độ phóng xạ môi trường đóng vai trò then chốt.

Vấn đề cốt lõi trong phân tích hạt nhân hiện đại là xác định hoạt độ của các mẫu môi trường thông qua phương pháp đo tương đối dựa trên mẫu chuẩn quốc tế. Tuy nhiên, khi tăng thể tích và khối lượng mẫu nhằm cải thiện thống kê xung đếm, hiện tượng tự hấp thụ bức xạ trong nền mẫu cùng sự thay đổi góc khối hình học dẫn đến suy giảm phi tuyến của tín hiệu ghi nhận. Nếu không khảo sát tường minh quy luật suy giảm này, các tính toán hoạt độ phóng xạ sẽ gặp sai số đáng kể.

Mục tiêu cụ thể của nghiên cứu là khảo sát thực nghiệm sự suy giảm của bức xạ theo độ cao mẫu chuẩn trên hệ phổ kế gamma phông thấp. Nghiên cứu tập trung xây dựng đường cong phụ thuộc của diện tích đỉnh hấp thụ toàn phần theo các mức độ cao 11mm, 22mm, 33mm và 44mm đối với các mẫu chuẩn IAEA-RGK-1, IAEA-RGTh-1 và mẫu đa nguyên tố IAEA-375. Đề tài được thực hiện tại Phòng thí nghiệm Vật lý Hạt nhân thuộc Trường Đại học Sư phạm Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang ý nghĩa thực tiễn sâu sắc, giúp tối ưu hóa chiều cao mẫu đo trong hình học hộp 2pi, tiết kiệm 30% đến 40% lượng mẫu chuẩn đắt tiền mà vẫn duy trì hiệu suất ghi nhận trên 90% giá trị bão hòa.

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 cơ sở lý thuyết tương tác bức xạ gamma với vật chất, bao gồm ba cơ chế chủ đạo:

Thứ nhất là hiệu ứng quang điện, trong đó lượng tử gamma truyền toàn bộ năng lượng cho electron liên kết ở lớp K (chiếm khoảng 80% xác suất). Tiết diện hấp thụ quang điện tỷ lệ thuận với Z lũy thừa 5 và tỷ lệ nghịch với E lũy thừa 3,5, đóng vai trò áp đảo ở vùng năng lượng thấp dưới 100 keV.

Thứ hai là tán xạ Compton, xảy ra khi photon va chạm với electron tự do và tán xạ một phần năng lượng theo hệ thức Compton với khối lượng nghỉ của electron m0c² tương đương 0,511 MeV. Đây là cơ chế tương tác ưu thế trong dải năng lượng trung bình từ 200 keV đến 3 MeV.

Thứ ba là hiệu ứng tạo cặp electron-positron, chỉ xảy ra khi năng lượng photon gamma vượt ngưỡng 1,022 MeV trong trường Coulomb của hạt nhân. Quá trình hủy cặp sau đó phát ra 2 photon gamma 0,511 MeV theo hướng ngược nhau.

Bên cạnh đó, nghiên cứu ứng dụng lý thuyết cân bằng phóng xạ thế kỷ (secular equilibrium). Trong chuỗi phân rã Thorium-232 với thời gian bán rã xấp xỉ 1,39 x 10¹⁰ năm, hoạt độ của đồng vị mẹ và các đồng vị con đạt trạng thái cân bằng bền vững. Nhờ đó, việc định lượng Thorium-232 vốn chỉ phát tia gamma yếu 63,8 keV (xác suất 0,27%) được chuyển đổi chính xác qua phép đo đỉnh gamma 583,2 keV của đồng vị con Thallium-208. Các thông số đặc trưng của hệ đo bao gồm độ phân giải năng lượng R xác định qua bề rộng ở một nửa giá trị cực đại (FWHM), hiệu suất ghi tuyệt đối và tỷ số đỉnh trên nền Compton (P/C) lớn hơn 46:1.

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

Nguồn dữ liệu thực nghiệm được thu thập từ 3 mẫu chuẩn chuẩn hóa quốc tế do Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) chứng nhận, bao gồm: mẫu đơn nguyên Kali IAEA-RGK-1, mẫu đơn nguyên Thori IAEA-RGTh-1 và mẫu đất đa nguyên tố IAEA-375 chứa hỗn hợp Kali-40, Thallium-208 và Cesium-137.

Quy trình chọn mẫu và chuẩn bị mẫu tuân thủ nghiêm ngặt kỹ thuật phân tích hạt nhân: đất được thu thập ở độ sâu từ 25cm đến 50cm tại vùng không bị xáo trộn, sấy khô ở nhiệt độ 105°C, nghiền mịn và rây đồng nhất để triệt tiêu hiệu ứng phân tầng. Mẫu được nạp vào hộp đo tiêu chuẩn đường kính 2pi tại 4 mức độ cao thực nghiệm là 11mm, 22mm, 33mm và 44mm, hình thành 12 bộ cấu hình đo lường độc lập.

Phương pháp phân tích quang phổ gamma phông thấp sử dụng detector bán dẫn Germanium siêu tinh khiết (HPGe) kiểu đồng trục p-type, vận hành dưới điện áp cao 2400V và làm lạnh liên tục bằng nitơ lỏng ở nhiệt độ -196°C. Lý do lựa chọn detector HPGe thay vì đầu dò nhấp nháy NaI(Tl) là bởi năng lượng cần thiết để tạo một cặp electron-lỗ trống trong tinh thể Ge chỉ là 2,96 eV, mang lại độ phân giải năng lượng vượt trội với FWHM chỉ 1,8 keV đến 2,0 keV tại đỉnh 1,33 MeV của Cobalt-60.

Detector được đặt trong buồng chì che chắn dày từ 5cm đến 10cm có lót các lớp kim loại chuyển tiếp gồm đồng, cadmi và parafin nhằm triệt tiêu bức xạ phông môi trường và tia vũ trụ. Tín hiệu xung được khuếch đại tuyến tính và xử lý qua khối phân tích đa kênh (MCA) 4096 kênh điều khiển bởi phần mềm chuyên dụng Maestro-32. Thời gian đo đếm được thiết lập 3 giờ (10.800 giây) đối với mẫu chuẩn đơn nguyên và 24 giờ (86.400 giây) đối với mẫu đa nguyên IAEA-375 nhằm tích lũy đủ số đếm thống kê với sai số dưới 1%.

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

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

Phát hiện thứ nhất: Diện tích đỉnh hấp thụ toàn phần thuần (Net Peak Area) của mọi đồng vị phóng xạ đều tăng theo chiều cao mẫu nhưng theo quy luật phi tuyến suy giảm dần tốc độ tăng trưởng.

Phát hiện thứ hai: Đối với mẫu chuẩn đơn nguyên IAEA-RGTh-1 tại đỉnh năng lượng 583,2 keV của Thallium-208, diện tích đỉnh thuần đạt 9.293,25 xung ở độ cao 11mm. Khi tăng độ cao lên 22mm (tăng gấp đôi chiều dày), diện tích đỉnh đạt 14.129,25 xung, tương ứng mức tăng 52,04%. Tuy nhiên, khi tăng tiếp lên 33mm, diện tích đỉnh đạt 17.724,25 xung (tốc độ tăng giảm xuống còn 25,44%), và tại 44mm chỉ đạt 20.151,25 xung (mức tăng thu hẹp còn 13,69%).

Phát hiện thứ ba: Trên mẫu chuẩn đơn nguyên IAEA-RGK-1 tại đỉnh năng lượng 1461 keV của Kali-40, diện tích đỉnh thuần tại mức 11mm đạt 6.783,75 xung sau khi trừ phông buồng chì 160,25 xung, và tiếp tục tiệm cận xu thế bão hòa ở các độ cao 33mm và 44mm.

Phát hiện thứ tư: Mẫu đa nguyên tố IAEA-375 đo trong 24 giờ ghi nhận diện tích đỉnh tổng của Kali-40 tăng từ 2.769 xung tại 11mm lên 4.372 xung tại 44mm (tăng 57,89%). Trong khi đó, các đỉnh năng lượng trung bình thấp như Cesium-137 (661,66 keV) và Thallium-208 (583,2 keV) có độ suy giảm mạnh hơn so với đỉnh năng lượng cao 1461 keV, khẳng định sự phụ thuộc rõ rệt của hệ số hấp thụ vào năng lượng tia gamma.

Thảo luận kết quả

Nguyên nhân vật lý của hiện tượng suy giảm tốc độ tăng diện tích đỉnh theo độ cao xuất phát từ hai yếu tố đồng thời:

Yếu tố tự hấp thụ trong thể tích mẫu: Các lượng tử gamma phát ra từ các lớp mẫu phía trên phải truyền qua chính môi trường vật chất của lớp mẫu bên dưới trước khi tới bề mặt nhạy của detector. Theo định luật hấp thụ bức xạ, chùm photon bị suy giảm theo hàm mũ exp(-µx), trong đó µ là hệ số suy giảm tuyến tính của chất nền đất. Năng lượng gamma càng thấp thì tiết diện hấp thụ càng lớn, dẫn đến tổn hao xung đếm càng cao.

Yếu tố suy giảm góc khối hình học: Khi độ cao mẫu tăng từ 11mm lên 44mm, khoảng cách trung bình từ các phần tử phóng xạ ở lớp trên cùng đến tinh thể detector tăng lên, làm giảm góc khối nhìn detector theo quy luật nghịch đảo bình phương khoảng cách.

Khi biểu diễn dữ liệu thực nghiệm trên đồ thị tọa độ hai chiều với trục hoành là độ cao mẫu (mm) và trục tung là diện tích đỉnh xung, đường cong thực nghiệm không đi theo đường thẳng tuyến tính mà uốn cong tiệm cận một đường nằm ngang. Phương trình hồi quy dạng phi tuyến mô tả chính xác quá trình chuyển tiếp từ vùng tăng trưởng nhanh sang vùng bão hòa. So sánh với các mô hình chuẩn quốc tế, kết quả này hoàn toàn tương thích với các nghiên cứu đo đạc bức xạ môi trường hiện đại. Khi độ dày mẫu vượt quá 35mm đến 40mm, việc tăng thêm thể tích mẫu chỉ đóng góp dưới 8% số đếm thực tế nhưng lại làm tăng gấp đôi chi phí mẫu chuẩn.

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

Thứ nhất, tối ưu hóa độ cao hình học mẫu chuẩn: Các phòng thí nghiệm phân tích hạt nhân cần chuẩn hóa độ cao mẫu đất và vật liệu môi trường trong hộp đo 2pi ở mức từ 25mm đến 30mm. Giải pháp này giúp tiết kiệm từ 30% đến 40% khối lượng mẫu chuẩn phân tích so với mức 44mm truyền thống, trong khi vẫn thu nhận được khoảng 85% đến 90% hiệu suất đếm cực đại. Thời gian áp dụng khuyến nghị này có thể thực hiện ngay trong vòng 1 tháng tới.

Thứ hai, xây dựng ma trận hiệu chỉnh tự hấp thụ bằng thuật toán mô phỏng: Đơn vị nghiên cứu và các trung tâm hạt nhân cần triển khai mô phỏng Monte Carlo trên nền tảng MCNP hoặc GEANT4 trong vòng 6 tháng tới. Mục tiêu là thiết lập đường cong hiệu chuẩn tự hấp thụ chi tiết cho từng dải năng lượng từ 60 keV đến 2000 keV, giảm thiểu sai số hình học xuống dưới 2%.

Thứ ba, nâng cấp hệ thống che chắn phông phóng xạ: Khuyến nghị các phòng thí nghiệm bổ sung lớp lót kim loại kép gồm thiếc dày 1mm và đồng dày 1,5mm bên trong thành buồng chì 10cm. Cải tiến này giúp triệt tiêu tia X huỳnh quang chì năng lượng 75 keV đến 85 keV, cắt giảm thêm 15% đến 20% phông nền liên tục trong vòng 3 tháng.

Thứ tư, ban hành quy chuẩn kiểm soát nồng độ Radon và phóng xạ môi trường: Cơ quan quản lý an toàn bức xạ quốc gia cần sớm ban hành văn bản quy chuẩn áp dụng mức can thiệp nồng độ Radon trong nhà ở ngưỡng 200 Bq/m³ theo khuyến nghị của IAEA, đồng thời tổ chức quan trắc định kỳ 12 tháng một lần tại các khu vực mỏ khoáng sản và vùng địa chất có nền phóng xạ tự nhiên cao.

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

Thứ nhất, kỹ sư và nghiên cứu viên vật lý hạt nhân: Tiếp cận quy trình vận hành chi tiết hệ phổ kế gamma bán dẫn HPGe, kỹ thuật xử lý phổ đa kênh với phần mềm Maestro-32, phương pháp xác định diện tích đỉnh thực và hiệu chuẩn hiệu suất theo hình học nguồn khối.

Thứ hai, chuyên viên quan trắc an toàn bức xạ và môi trường: Ứng dụng phương pháp đo tương đối để định lượng chính xác hàm lượng các đồng vị phóng xạ tự nhiên (Uranium-238, Thorium-232, Kali-40) và nhân tạo (Cesium-137) trong đất, nước, thực phẩm và vật liệu xây dựng theo tiêu chuẩn kiểm soát an toàn.

Thứ ba, giảng viên và sinh viên chuyên ngành Vật lý - Kỹ thuật hạt nhân: Sử dụng làm tài liệu tham khảo giảng dạy chuyên sâu về thực nghiệm tương tác gamma với vật chất, cơ chế cân bằng phóng xạ thế kỷ và phương pháp tối ưu hóa thiết bị ghi đo bức xạ.

Thứ tư, các nhà hoạch định chính sách an toàn bức xạ và y tế công cộng: Tham khảo cơ sở dữ liệu thực nghiệm để xây dựng định mức kỹ thuật, quy chuẩn can thiệp phơi nhiễm khí Radon và tiêu chuẩn giám sát an toàn môi trường xung quanh các cơ sở ứng dụng bức xạ.

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

Vì sao diện tích đỉnh phổ gamma không tăng tuyến tính theo chiều cao mẫu chuẩn? Khi tăng chiều cao mẫu từ 11mm lên 44mm, các lớp mẫu phía trên làm tăng khoảng cách tới detector khiến góc khối thu nhận suy giảm. Đồng thời, tia gamma phát ra từ các lớp trên bị hấp thụ một phần bởi chính khối vật chất phía dưới. Thực nghiệm mẫu IAEA-RGTh-1 chứng minh tốc độ tăng diện tích đỉnh giảm mạnh từ 52,04% (ở 22mm) xuống chỉ còn 13,69% (ở 44mm).

Detector bán dẫn HPGe có ưu điểm gì vượt trội so với detector nhấp nháy NaI(Tl)? Detector HPGe có độ phân giải năng lượng cực cao với FWHM chỉ 1,8 keV tại đỉnh 1,33 MeV của Cobalt-60, vượt trội hoàn toàn so với mức phân giải 5% đến 10% của đầu dò NaI(Tl). Nhờ đó, hệ đo HPGe có khả năng tách biệt các đỉnh năng lượng nằm sát nhau, điển hình như đỉnh 185,72 keV của Uranium-235 và đỉnh 186,21 keV của Radium-226.

Làm thế nào để định lượng Thorium-232 khi đồng vị này phát tia gamma cường độ quá yếu? Thorium-232 phát tia gamma 63,8 keV với hiệu suất phân rã chỉ 0,27% nên rất khó đo trực tiếp. Nhờ hiện tượng cân bằng phóng xạ thế kỷ trong chuỗi Thorium, hoạt độ của đồng vị con Thallium-208 bằng đúng hoạt độ của Thorium-232. Do đó, người phân tích có thể định lượng gián tiếp thông qua việc đo đỉnh gamma năng lượng cao 583,2 keV của Thallium-208.

Tiêu chuẩn an toàn quốc tế quy định nồng độ khí Radon trong nhà ở mức nào? Theo tiêu chuẩn an toàn bức xạ của IAEA, nồng độ khí Radon trong nhà ở không được vượt quá dải 200 Bq/m³ đến 600 Bq/m³/năm (tương đương 0,6 đến 1,7 Bq/m³/ngày). Tại Mỹ, luật môi trường quy định mức giới hạn can thiệp nghiêm ngặt hơn ở mức dưới 148 Bq/m³/năm nhằm ngăn ngừa tổn thương phế nang phổi do hạt alpha.

Lợi ích thực tiễn của việc xác định điểm bão hòa độ cao mẫu chuẩn là gì? Xác định điểm bão hòa (khoảng 30mm đến 35mm trong hộp 2pi) giúp phòng thí nghiệm chọn đúng lượng mẫu chuẩn vừa đủ để đạt trên 90% hiệu suất đếm cực đại. Điều này giúp tránh lãng phí các chất chuẩn phóng xạ đắt tiền từ IAEA và rút ngắn thời gian chuẩn bị mẫu mà vẫn đảm bảo độ chính xác phép đo.

Kết luận

  • Luận văn đã khảo sát thành công sự suy giảm bức xạ gamma theo độ cao mẫu chuẩn từ 11mm đến 44mm trên hệ phổ kế phông thấp sử dụng detector bán dẫn HPGe.
  • Chứng minh bằng số liệu thực nghiệm quy luật tăng trưởng phi tuyến của diện tích đỉnh phổ, với tốc độ tăng giảm từ 52,04% xuống 13,69% do hiệu ứng tự hấp thụ và suy giảm góc khối.
  • Chuẩn hóa quy trình phân tích định lượng chuỗi Thorium-232 thông qua trạng thái cân bằng phóng xạ thế kỷ tại đỉnh 583,2 keV của Thallium-208.
  • Đề xuất độ cao tối ưu 25mm đến 30mm cho hình học hộp đo 2pi, giúp tiết kiệm từ 30% đến 40% khối lượng mẫu chuẩn IAEA quý hiếm.
  • Cung cấp cơ sở thực nghiệm vững chắc phục vụ công tác giám sát phóng xạ môi trường, đo kiểm Radon-222 và kiểm chuẩn an toàn bức xạ hạt nhân.

Đóng góp chính của đề tài là xây dựng đường cong thực nghiệm chuẩn xác về sự phụ thuộc diện tích đỉnh theo chiều cao mẫu, giải quyết bài toán tối ưu hóa giữa hiệu suất ghi đo và chi phí mẫu chuẩn. Trong giai đoạn 12 tháng tiếp theo, nghiên cứu cần được mở rộng khảo sát trên bình đo Marinelli thể tích lớn và tích hợp thuật toán mô phỏng Monte Carlo để hoàn thiện ma trận hiệu chuẩn tự hấp thụ toàn diện. Các đơn vị nghiên cứu và phòng thí nghiệm chuyên ngành nên áp dụng ngay quy chuẩn chiều cao tối ưu này vào quy trình đo kiểm định kỳ để nâng cao độ chính xác và hiệu quả kinh tế.