Luận văn thạc sĩ xác định hoạt độ các nguyên tố phóng xạ bằng phương pháp phổ kế gamma kết hợp labsocs và mô phỏng monte carlo geant4

Luận văn thạc sĩ: Nghiên cứu xác định hoạt độ phóng xạ bằng phổ kế gamma, kết hợp LABSOCS và mô phỏng Monte Carlo GEANT4. Phân tích định lượng chính xác.

• Tác giả: Lê Quốc Việt
• Người hướng dẫn: TS. Phan Việt Cương
• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Vật lý nguyên tử và hạt nhân
• Đề tài: Xác định hoạt độ các nguyên tố phóng xạ bằng phương pháp phổ kế gamma kết hợp LabSOCS và mô phỏng Monte-Carlo GEANT4
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2019
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• DANH MỤC CÁC KÍ HIỆU, VIẾT TẮT
• LỜI NÓI ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Phân rã phóng xạ và nguồn gốc của bức xạ gamma
• 1.1.1. Phương trình toán học biểu diễn quá trình phân rã
• 1.2. Tương tác của bức xạ gamma với vật chất
• 1.2.1. Hiệu ứng quang điện
• 1.2.2. Hiệu ứng tán xạ Compton
• 1.2.3. Hiệu ứng tạo cặp
• 1.3. Độ suy giảm cường độ của chùm gamma khi đi qua lớp vật chất
• 1.3.1. Hệ số suy giảm tuyến tính
• 1.3.2. Hệ số suy giảm khối
• 1.4. Hệ số tự hấp thụ đối với nguồn thể tích
• 1.4.1. Xác định hệ số tự hấp thụ
• 1.4.2. Xác định hệ số suy giảm tuyến tính
• 1.5. Cơ sở vật lý của phương pháp gamma xác định hoạt động phóng xạ của các đồng vị
• 1.5.1. Phân tích định tính
• 1.5.2. Phân tích định lượng
• 2. CHƯƠNG 2: PHƯƠNG PHÁP MÔ PHỎNG MONTE – CARLO GEANT4 VÀ PHẦN MỀM LABSOCS
• 2.1. Phương pháp mô phỏng Monte-Carlo GEANT4
• 2.2. Phần mềm LabSOCS
• 2.2.1. Giới thiệu về phần mềm LabSOCS
• 2.2.2. Sử dụng phần mềm LabSOCS
• 2.3. Khảo sát sai số do sự khác biệt về thành phần mẫu và mật độ khối trên LabSOCS
• 3. CHƯƠNG 3: KẾT QUẢ THỰC NGHIỆM
• 3.1. Kết quả phân tích mẫu chuẩn
• 3.1.1. Mô tả các mẫu thực nghiệm
• 3.1.2. Kết quả tính toán với mẫu đất
• 3.1.3. Kết quả tính toán với mẫu thực vật
• 3.1.4. Kết quả tính toán với mẫu nước
• 3.2. Kết quả so sánh đối chứng
• TÀI LIỆU THAM KHẢO
• Các mẫu cấu hình hộp đo có sẵn trong LabSOCS
• Thông tin các loại vật liệu có sẵn trong thư viện của LabSOCS
• Hình ảnh một số mẫu chuẩn

Tóm tắt

I. Khám phá phương pháp phổ kế gamma xác định hoạt độ phóng xạ

Phương pháp phổ kế gamma là một kỹ thuật phân tích hạt nhân hiện đại. Nó được ứng dụng rộng rãi để xác định hoạt độ các nguyên tố phóng xạ. Kỹ thuật này có ưu điểm không phá hủy mẫu, cho phép phân tích định tính và định lượng các đồng vị phát ra tia gamma trong vật liệu. Cơ sở của phương pháp dựa trên việc ghi nhận và phân tích phổ năng lượng của bức xạ gamma. Bức xạ này phát ra khi hạt nhân không bền tự phân rã để chuyển về trạng thái bền hơn. Mỗi đồng vị phóng xạ tự nhiên hay nhân tạo đều phát ra tia gamma với các mức năng lượng đặc trưng, giống như “dấu vân tay” của chúng. Hệ đo sử dụng một detector bán dẫn siêu tinh khiết, thường là phổ kế gamma HPGe (Germanium siêu tinh khiết), có độ phân giải năng lượng rất cao. Nhờ đó, hệ thống có thể tách biệt rõ ràng các vạch năng lượng gamma gần nhau, giúp nhận diện chính xác các đồng vị có trong mẫu. Việc ứng dụng phương pháp này đóng vai trò quan trọng trong nhiều lĩnh vực như vật lý hạt nhân ứng dụng, quan trắc môi trường, địa chất và an toàn bức xạ.

1.1. Cơ sở vật lý của phương pháp phân tích phổ gamma

Cơ sở vật lý của phương pháp này nằm ở sự tương tác của bức xạ gamma với vật chất. Khi tia gamma đi qua vật chất, nó tương tác với các nguyên tử thông qua ba hiệu ứng chính: hiệu ứng quang điện, tán xạ Compton và hiệu ứng tạo cặp. Trong hiệu ứng quang điện, photon gamma truyền toàn bộ năng lượng cho một electron. Trong tán xạ Compton, photon chỉ truyền một phần năng lượng. Với hiệu ứng tạo cặp, photon năng lượng cao (trên 1.022 MeV) biến thành một cặp electron-positron gần hạt nhân. Các tương tác này làm suy giảm cường độ chùm tia gamma và tạo ra các tín hiệu điện trong detector. Detector bán dẫn siêu tinh khiết HPGe hoạt động dựa trên việc thu thập các cặp electron-lỗ trống được tạo ra bởi bức xạ gamma. Số lượng cặp này tỉ lệ thuận với năng lượng của photon gamma bị hấp thụ. Tín hiệu điện sau đó được khuếch đại và xử lý để tạo thành một phổ năng lượng, trong đó các đỉnh năng lượng (photopeak) cho biết sự hiện diện của các đồng vị phóng xạ cụ thể.

1.2. Nguyên lý phân tích định tính và định lượng đồng vị

Việc phân tích phổ gamma bao gồm hai giai đoạn chính: định tính và định lượng. Phân tích định tính là xác định sự có mặt của các đồng vị phóng xạ trong mẫu. Điều này được thực hiện bằng cách so sánh vị trí các đỉnh năng lượng trên phổ đo được với thư viện năng lượng gamma chuẩn của các đồng vị đã biết. Mỗi đỉnh tương ứng với một mức năng lượng đặc trưng. Phân tích định lượng nhằm mục đích đo hoạt độ phóng xạ của từng đồng vị đã xác định. Hoạt độ (A) của một đồng vị được tính toán dựa trên diện tích đỉnh hấp thụ toàn phần (S), thời gian đo (t), cường độ phát tia gamma (Iγ) và hiệu suất ghi của detector (ε) theo công thức: A = S / (t * ε * Iγ). Yếu tố quan trọng nhất và cũng phức tạp nhất trong công thức này chính là hiệu suất ghi, đòi hỏi phải được xác định một cách chính xác.

II. Thách thức hiệu chuẩn hiệu suất ghi và sai số phép đo gamma

Để đảm bảo kết quả đo hoạt độ phóng xạ chính xác, việc hiệu chuẩn hiệu suất ghi của hệ phổ kế là bắt buộc. Hiệu suất ghi (efficiency) là tỷ lệ giữa số xung được ghi nhận tại một đỉnh năng lượng và tổng số photon gamma có cùng năng lượng đó phát ra từ mẫu. Hiệu suất này phụ thuộc vào nhiều yếu tố như năng lượng gamma, hình học giữa mẫu và detector, và đặc tính của chính mẫu đo. Phương pháp hiệu chuẩn truyền thống sử dụng các mẫu chuẩn đo lường vật lý có hoạt độ đã biết. Tuy nhiên, phương pháp này đối mặt với nhiều thách thức lớn, gây ra sai số phép đo đáng kể và làm tăng chi phí cũng như thời gian phân tích. Các khó khăn này thúc đẩy sự phát triển của các phương pháp hiệu chuẩn không dùng nguồn, dựa trên tính toán và mô phỏng, nhằm nâng cao độ chính xác và hiệu quả trong vật lý hạt nhân ứng dụng.

2.1. Hạn chế của việc sử dụng mẫu chuẩn vật lý truyền thống

Sử dụng mẫu chuẩn vật lý tồn tại nhiều hạn chế. Thứ nhất, chi phí cho các mẫu chuẩn, đặc biệt là các mẫu được chứng nhận bởi các tổ chức quốc tế như IAEA, là rất cao (khoảng 180 EUR/mẫu theo tài liệu). Thứ hai, thời gian vận chuyển và các thủ tục liên quan thường kéo dài. Thứ ba, các đồng vị có chu kỳ bán rã ngắn đòi hỏi phải thay thế mẫu chuẩn thường xuyên. Hạn chế lớn nhất là sự khác biệt về ma trận (thành phần và mật độ) giữa mẫu chuẩn và mẫu phân tích. Rất khó để có một mẫu chuẩn vật lý hoàn toàn tương đồng với mọi loại mẫu thực tế (đất, tro bay, quặng, thực vật...). Sự khác biệt này dẫn đến sai số nghiêm trọng do hiệu ứng tự hấp thụ bức xạ trong mẫu, làm cho đường cong hiệu suất xây dựng từ mẫu chuẩn không còn chính xác khi áp dụng cho mẫu đo.

2.2. Yếu tố ảnh hưởng đến độ chính xác khi đo hoạt độ phóng xạ

Nhiều yếu tố gây ra sai số phép đo trong phổ kế gamma. Yếu tố quan trọng hàng đầu là sự tự hấp thụ của tia gamma bên trong mẫu thể tích. Bức xạ gamma phát ra từ các vị trí khác nhau trong mẫu sẽ đi qua những quãng đường khác nhau trong vật liệu trước khi đến detector, dẫn đến mức độ suy giảm khác nhau. Hiện tượng này phụ thuộc vào mật độ khối và thành phần hóa học của mẫu. Ngoài ra, hình học đo (khoảng cách và vị trí tương đối giữa mẫu và detector), hình dạng và kích thước của mẫu cũng ảnh hưởng trực tiếp đến hiệu suất ghi. Nếu không hiệu chỉnh chính xác các yếu tố này, kết quả đo hoạt độ phóng xạ có thể sai lệch đáng kể, đặc biệt đối với các tia gamma năng lượng thấp.

III. Cách mạng hóa hiệu chuẩn hiệu suất ghi bằng phần mềm LABSOCS

Để khắc phục những hạn chế của mẫu chuẩn vật lý, phần mềm LABSOCS (Laboratory Sourceless Calibration Software) ra đời như một giải pháp đột phá. Đây là một công cụ thương mại được tích hợp vào phần mềm phân tích phổ Genie 2000 của Canberra. LABSOCS cho phép xây dựng đường cong hiệu chuẩn hiệu suất ghi một cách chính xác cho hầu hết mọi dạng hình học mẫu mà không cần đến mẫu chuẩn vật lý. Nguyên lý của phần mềm dựa trên việc mô phỏng quá trình đo lường, kết hợp các thông số đặc trưng của detector đã được nhà sản xuất đo đạc chi tiết và các mô hình toán học từ chương trình MCNP (Monte Carlo N-Particle). Bằng cách nhập vào các thông tin về cấu hình đo và đặc tính của mẫu, LABSOCS có thể tính toán đường cong hiệu suất chỉ trong vài phút, giúp giảm thiểu đáng kể sai số phép đo liên quan đến sự khác biệt ma trận mẫu.

3.1. Nguyên lý và ưu điểm chính của phần mềm LABSOCS

Nguyên lý cốt lõi của phần mềm LABSOCS là sử dụng một mô hình toán học đã được xác thực để thay thế cho việc đo đạc thực nghiệm với mẫu chuẩn. Phần mềm chứa một cơ sở dữ liệu chi tiết về các đặc tính của từng detector cụ thể, được nhà sản xuất đo lường và xác nhận. Các thông số này bao gồm kích thước chính xác của tinh thể Germanium, độ dày lớp chết (dead-layer), vị trí trong vỏ bọc... Người dùng chỉ cần định nghĩa hình học đo (hình dạng hộp chứa mẫu, khoảng cách tới detector) và các thuộc tính của mẫu (thành phần vật liệu, mật độ khối). LABSOCS sẽ sử dụng các thông tin này để thực hiện một phép tính toán dựa trên phương pháp Monte Carlo, mô phỏng quá trình tương tác của bức xạ với vật chất và xác định hiệu suất ghi. Ưu điểm của nó là tốc độ, sự linh hoạt và khả năng loại bỏ sai số do sự không tương đồng giữa mẫu chuẩn và mẫu đo.

3.2. Quy trình ứng dụng LABSOCS xây dựng đường cong hiệu suất

Quy trình thực hiện hiệu chuẩn bằng LABSOCS tương đối đơn giản. Đầu tiên, người dùng khởi tạo một tệp hình học mới trong giao diện của phần mềm. Sau đó, lựa chọn đúng mã detector đang sử dụng từ danh sách. Tiếp theo, khai báo chi tiết các thông số của cấu hình đo, bao gồm loại hộp chứa (ví dụ: hình trụ, bình Marinelli), kích thước, vật liệu và độ dày của thành hộp. Quan trọng nhất là phải định nghĩa chính xác vật liệu và mật độ khối của mẫu đo. LABSOCS cung cấp một thư viện vật liệu có sẵn và cho phép người dùng tự tạo vật liệu mới. Cuối cùng, phần mềm sẽ tính toán và tạo ra đường cong hiệu suất ghi cho dải năng lượng quan tâm. Kết quả này sau đó được sử dụng trực tiếp để đo hoạt độ phóng xạ từ phổ của mẫu thực tế.

IV. Bí quyết mô phỏng detector bức xạ bằng Monte Carlo Geant4

Bên cạnh các giải pháp thương mại như LABSOCS, mô phỏng Monte Carlo Geant4 là một công cụ mã nguồn mở cực kỳ mạnh mẽ trong lĩnh vực vật lý hạt nhân ứng dụng. Geant4 (Geometry and Tracking 4) là một bộ công cụ được phát triển bởi CERN, cho phép mô phỏng chi tiết đường đi và sự tương tác của các hạt qua vật chất. Không giống LABSOCS, Geant4 đòi hỏi người dùng phải tự xây dựng toàn bộ mô hình thí nghiệm từ đầu, bao gồm cả việc mô tả hình học chi tiết của detector bán dẫn siêu tinh khiết, nguồn phóng xạ và các quá trình vật lý liên quan. Mặc dù phức tạp hơn, phương pháp này mang lại sự linh hoạt tối đa và là một công cụ nghiên cứu vô giá để hiểu sâu sắc các quá trình vật lý xảy ra trong phép đo, cũng như để kiểm chứng chéo với các phương pháp khác. Đây là một công cụ không thể thiếu cho việc mô phỏng detector bức xạ.

4.1. Cấu trúc và cơ chế làm việc của bộ công cụ Geant4

Geant4 được xây dựng dựa trên ngôn ngữ lập trình hướng đối tượng C++. Một chương trình mô phỏng Geant4 hoàn chỉnh yêu cầu người dùng phải định nghĩa ba lớp cơ bản. Lớp DetectorConstruction dùng để mô tả toàn bộ cấu trúc hình học của hệ đo, từ tinh thể detector, vỏ bọc, đến hộp chứa mẫu và vật liệu của chúng. Lớp PhysicsList cho phép lựa chọn và khai báo các mô hình vật lý sẽ được sử dụng để mô phỏng tương tác của bức xạ với vật chất, ví dụ như các gói vật lý điện từ tiêu chuẩn. Cuối cùng, lớp PrimaryGeneratorAction định nghĩa các đặc tính ban đầu của nguồn bức xạ, như loại hạt, năng lượng, vị trí và hướng phát. Khi chương trình chạy, Geant4 sẽ theo dõi từng hạt, tính toán các tương tác của chúng và ghi lại các thông tin cần thiết, ví dụ như năng lượng lắng đọng trong thể tích hoạt động của detector.

4.2. Xây dựng mô hình hiệu chuẩn hiệu suất ghi bằng Geant4

Để hiệu chuẩn hiệu suất ghi bằng mô phỏng Geant4, bước đầu tiên là xây dựng một mô hình số hóa chính xác của hệ phổ kế gamma HPGe. Mô hình này phải tái tạo lại tất cả các thành phần quan trọng ảnh hưởng đến quá trình đo. Sau khi mô hình được xây dựng, một nguồn bức xạ ảo với số lượng lớn các phân rã đã biết được mô phỏng. Chương trình sẽ theo dõi từng photon gamma phát ra và xác định xem nó có được ghi nhận bởi detector hay không. Bằng cách đếm số photon được ghi nhận trong đỉnh hấp thụ toàn phần và chia cho tổng số photon phát ra ở năng lượng tương ứng, ta có thể tính được hiệu suất ghi tuyệt đối. Quá trình này được lặp lại cho nhiều mức năng lượng khác nhau để xây dựng một đường cong hiệu suất hoàn chỉnh. Kết quả từ Geant4 có thể được dùng để kiểm chứng và so sánh với kết quả từ LABSOCS và thực nghiệm.

V. So sánh kết quả xác định hoạt độ phóng xạ từ các phương pháp

Luận văn đã tiến hành áp dụng cả ba phương pháp để phân tích các mẫu thực tế: sử dụng mẫu chuẩn vật lý, phần mềm LABSOCS, và mô phỏng Geant4 kết hợp đo hệ số truyền qua. Đối tượng nghiên cứu là các mẫu chuẩn do IAEA cung cấp, bao gồm mẫu đất, mẫu thực vật và mẫu nước. Mục tiêu là so sánh và đánh giá độ chính xác, ưu nhược điểm của từng phương pháp trong các điều kiện đo lường khác nhau. Các kết quả thực nghiệm cho thấy sự vượt trội rõ rệt của các phương pháp không dùng nguồn, đặc biệt là trong các trường hợp mà mẫu đo và mẫu chuẩn có sự khác biệt lớn về mật độ và thành phần. Việc đối chứng kết quả giữa các phương pháp giúp khẳng định độ tin cậy của việc đo hoạt độ phóng xạ khi không có mẫu chuẩn phù hợp.

5.1. Phân tích mẫu chuẩn IAEA Đất thực vật và nước

Kết quả cho thấy đối với mẫu nước, nơi mẫu chuẩn và mẫu đo có thành phần gần như tương đồng (chủ yếu là nước), cả ba phương pháp đều cho kết quả tốt và tương đương nhau. Tuy nhiên, với mẫu đất và thực vật, sự khác biệt trở nên rõ rệt. Phương pháp dùng mẫu chuẩn vật lý cho thấy sai số phép đo lớn. Ví dụ, khi dùng mẫu chuẩn đất IAEA-375 (khối lượng ~260g) để đo mẫu đất S5 (khối lượng ~140g), độ lệch tương đối rất cao. Ngược lại, khi sử dụng phần mềm LABSOCS (giả định vật liệu là 'Drydirt' và hiệu chỉnh mật độ), sai số giảm đi đáng kể. Điều này chứng tỏ khả năng hiệu chỉnh hiệu quả của LABSOCS đối với sự khác biệt về mật độ và hiệu ứng tự hấp thụ, giúp nâng cao độ chính xác khi phân tích phổ gamma.

5.2. Đối chứng kết quả giữa LABSOCS và mô phỏng Monte Carlo

Luận văn đã thực hiện so sánh đối chứng kết quả giữa LABSOCS và phương pháp mô phỏng Geant4 kết hợp đo hệ số truyền qua. Các kết quả xác định hoạt độ của các đồng vị phóng xạ tự nhiên như K-40, Ac-228 trong mẫu đất cho thấy sự tương đồng cao giữa hai phương pháp mô phỏng. Sai khác giữa kết quả của LABSOCS và Geant4 là rất nhỏ, thường nằm trong khoảng cho phép của sai số phép đo. Sự phù hợp này khẳng định rằng cả hai công cụ mô phỏng đều có khả năng tái tạo chính xác quá trình đo lường vật lý. Việc này cung cấp một bằng chứng mạnh mẽ về độ tin cậy của các phương pháp hiệu chuẩn không dùng nguồn, mở ra khả năng thay thế cho các chuẩn đo lường vật lý đắt đỏ và kém linh hoạt.

VI. Tương lai ứng dụng LABSOCS và Geant4 trong vật lý hạt nhân

Việc kết hợp phương pháp phổ kế gamma với các công cụ mô phỏng như phần mềm LABSOCS và mô phỏng Geant4 đang mở ra một kỷ nguyên mới cho lĩnh vực phân tích phóng xạ. Các phương pháp này không chỉ giúp nâng cao độ chính xác mà còn tiết kiệm đáng kể thời gian và chi phí so với phương pháp truyền thống. Khả năng xử lý linh hoạt các loại mẫu đa dạng mà không cần mẫu chuẩn tương ứng là một lợi thế vượt trội. Trong tương lai, việc phát triển và ứng dụng các phần mềm mô phỏng sẽ tiếp tục là xu hướng chủ đạo. Chúng đóng vai trò then chốt trong các lĩnh vực yêu cầu độ chính xác cao như quan trắc môi trường, địa chất, khảo cổ và đặc biệt là an toàn bức xạ, góp phần bảo vệ sức khỏe con người và môi trường sống.

6.1. Tổng kết ưu điểm vượt trội của phương pháp mô phỏng

Các phương pháp mô phỏng mang lại nhiều ưu điểm. Thứ nhất, chúng loại bỏ được sai số do sự khác biệt về khối lượng, mật độ và thành phần giữa mẫu chuẩn và mẫu đo. Thứ hai, tiết kiệm chi phí mua sắm, lưu trữ và xử lý mẫu chuẩn phóng xạ. Thứ ba, rút ngắn đáng kể thời gian phân tích, vì chỉ cần đo mẫu một lần và thực hiện tính toán hiệu suất trên máy tính trong vài phút. Cuối cùng, chúng cung cấp sự linh hoạt tối đa, cho phép phân tích hầu hết mọi loại mẫu với bất kỳ cấu hình hình học nào. Những lợi ích này làm cho các phương pháp mô phỏng trở thành công cụ không thể thiếu trong các phòng thí nghiệm vật lý hạt nhân ứng dụng hiện đại.

6.2. Hướng phát triển và ứng dụng trong an toàn bức xạ

Trong tương lai, các nghiên cứu có thể tập trung vào việc cải tiến các mô hình vật lý trong Geant4 để mô tả chính xác hơn các tương tác ở năng lượng thấp. Việc kết hợp dữ liệu từ nhiều kỹ thuật đo khác nhau để xây dựng mô hình mẫu chính xác hơn cho LABSOCS cũng là một hướng đi tiềm năng. Trong lĩnh vực an toàn bức xạ, các công cụ này có thể được dùng để đánh giá nhanh hoạt độ phóng xạ tại hiện trường (sử dụng phiên bản ISOCS), lập bản đồ ô nhiễm phóng xạ, hay tính toán liều chiếu cho nhân viên bức xạ và dân chúng. Sự phát triển của các phương pháp này sẽ góp phần nâng cao năng lực giám sát và ứng phó sự cố bức xạ, đảm bảo an ninh hạt nhân quốc gia.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật phân tích hạt nhân hiện đại, hơn 85% các phép đo hoạt độ phóng xạ môi trường sử dụng hệ phổ kế gamma độ phân giải cao trang bị đầu dò bán dẫn Germanium siêu tinh khiết (HPGe). Tuy nhiên, rào cản lớn nhất của phương pháp phân tích định lượng truyền thống là sự phụ thuộc vào các mẫu chuẩn phóng xạ vật lý. Chi phí cho mỗi mẫu chuẩn nhập khẩu từ Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) lên tới khoảng 180 EUR (tương đương gần 5.000.000 VND), đi kèm thời gian vận chuyển kéo dài từ 2 đến 4 tuần và giới hạn nghiêm ngặt về chu kỳ bán rã của các đồng vị. Bên cạnh đó, sự không tương thích về hình học, mật độ khối và thành phần hóa học giữa mẫu chuẩn và mẫu đo thực tế gây ra sai số tự hấp thụ đáng kể, đặc biệt đối với các nguồn dạng thể tích.

Nhằm giải quyết triệt để vấn đề này, luận văn tập trung nghiên cứu phương pháp xác định hoạt độ các nguyên tố phóng xạ bằng phổ kế gamma kết hợp phần mềm hiệu chuẩn không dùng nguồn LabSOCS và công cụ mô phỏng Monte-Carlo GEANT4. Mục tiêu cụ thể là xây dựng đường cong hiệu suất ghi chính xác cho hệ phổ kế HPGe mà không cần sử dụng nguồn chuẩn phóng xạ, đồng thời khảo sát định lượng các yếu tố gây sai số do hình học hộp đo, mật độ khối và thành phần ma trận mẫu.

Nghiên cứu được triển khai thực nghiệm tại Trung tâm Phân tích thuộc Viện Công nghệ Xạ hiếm và Viện Vật lý, thực hiện trên 9 bộ mẫu chuẩn quốc tế của IAEA thuộc dải năng lượng từ 46,5 keV đến 1836 keV. Ý nghĩa thực tiễn của công trình thể hiện qua việc rút ngắn thời gian xây dựng đường chuẩn hiệu suất từ nhiều ngày xuống dưới 2 phút, tiết kiệm 100% chi phí mua nguồn chuẩn định kỳ và đảm bảo độ chính xác phép đo với sai lệch dưới 10% trong vùng năng lượng trên 80 keV.

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 lý thuyết phân rã phóng xạ hạt nhân và quy luật suy giảm chùm tia bức xạ điện từ theo định luật Beer-Lambert: $I = I_0 \cdot e^{-\mu x}$, trong đó cường độ bức xạ suy giảm phụ thuộc vào bề dày lớp vật chất và hệ số suy giảm tuyến tính. Khi bức xạ gamma truyền qua môi trường, ba hiệu ứng tương tác cơ bản diễn ra gồm: hiệu ứng quang điện chiếm ưu thế ở năng lượng thấp với tiết diện tỷ lệ với $Z^5 / E_\gamma^{3,5}$; hiệu ứng tán xạ Compton chiếm ưu thế ở dải năng lượng trung bình với tiết diện tỷ lệ với $Z / E_\gamma$; và hiệu ứng tạo cặp xuất hiện khi năng lượng gamma vượt ngưỡng 1022 keV với tiết diện tỷ lệ với $Z^2 \ln E_\gamma$.

Các mô hình giải tích làm khớp đường cong hiệu suất ghi của detector bán dẫn như mô hình Gray - Ahmad áp dụng cho dải năng lượng 80-1850 keV và mô hình Sanchez-Reyes cho dải 63-3050 keV được sử dụng làm cơ sở đối sánh. Bên cạnh đó, các khái niệm trọng tâm bao gồm: hiệu suất ghi tuyệt đối tại đỉnh hấp thụ toàn phần, hệ số tự hấp thụ của nguồn thể tích, hệ số suy giảm khối ($\mu/\rho$), mật độ khối ($\rho$) và công thức xác định hoạt độ phóng xạ $A = S / (\varepsilon \cdot I_\gamma \cdot t)$, trong đó $S$ là diện tích đỉnh đã trừ phông, $\varepsilon$ là hiệu suất ghi, $I_\gamma$ là cường độ phát xạ và $t$ là thời gian đo thực.

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

Nguồn dữ liệu thực nghiệm gồm 9 mẫu chuẩn được chứng nhận bởi IAEA, phân bổ trên 3 dạng ma trận môi trường phổ biến: mẫu đất (S5-Bottle 108, IAEA-447, S4-Bottle 169), mẫu thực vật (S4-Bottle 108 cỏ khô, IAEA-330 rau bina), và mẫu nước (S1, S2-Bottle 178 và S1, S2-Bottle 169). Phương pháp chọn mẫu có chủ đích nhằm đại diện cho các trạng thái vật chất khác nhau với hai cấu hình hình học hộp đo chuẩn hóa: dạng hình trụ (đường kính 85,76 mm, chiều cao 30 mm, vỏ polypropylene dày 0,82 mm) và dạng bình giếng Marinelli 500 ml bằng polystyrene.

Quy trình phân tích kết hợp hai kỹ thuật mô phỏng hiện đại: phần mềm LabSOCS tích hợp trong hệ Genie 2000 và bộ công cụ mô phỏng hạt nhân Monte-Carlo GEANT4 viết bằng C++. Lý do lựa chọn LabSOCS là khả năng nhập chính xác các thông số cấu trúc tinh thể thực của detector HPGe (bao gồm lớp chết dead-layer, góc xiên, vùng làm lạnh), cho phép loại bỏ sai số chế tạo và tự động tính toán hiệu suất ghi cùng hệ số tự hấp thụ theo mật độ khối chỉ trong 60 đến 120 giây. Đồng thời, mô phỏng GEANT4 với 3 lớp cấu trúc bắt buộc (Detector Construction, Physics List, Primary Generator Action) kết hợp đo thực nghiệm hệ số truyền qua được sử dụng làm phương pháp đối chứng độc lập nhằm kiểm chứng độ tin cậy khoa học.

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

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

Thứ nhất, việc ứng dụng LabSOCS đã loại bỏ hoàn toàn sai số do chênh lệch mật độ khối giữa mẫu chuẩn và mẫu đo. Đối với mẫu đất S5-Bottle 108 có khối lượng thực 140g, khi áp dụng phương pháp mẫu chuẩn truyền thống (sử dụng mẫu chuẩn IAEA-375 nặng 260g trong cùng hộp đo), kết quả bị sai lệch nghiêm trọng trên 35%. Tuy nhiên, khi dùng LabSOCS hiệu chỉnh mật độ khối về giá trị thực 0,7565 g/cm3, sai khác hoạt độ của Ac-228 giảm xuống còn 16%, và đồng vị Cs-137 được xác định chính xác ở mức 103,4 Bq/kg. Trên mẫu đất chuẩn IAEA-447, hoạt độ đo được của Cs-137 là 425 Bq/kg so với giá trị chuẩn 449 Bq/kg (độ sai khác chỉ 5,6%) và K-40 là 573 Bq/kg so với 550 Bq/kg (độ sai khác 4,2%).

Thứ hai, phương pháp đạt độ chính xác rất cao trên các ma trận thực vật và nước. Đối với mẫu rau bina IAEA-330, hoạt độ Cs-137 tính toán qua LabSOCS đạt 1248 Bq/kg so với giá trị chuẩn 1235 Bq/kg, đạt độ sai khác ấn tượng chỉ 1,1%, trong khi K-40 có độ lệch 8,3%. Đối với các mẫu nước S1 và S2 (Bottle 178) trong bình Marinelli 500 ml, hoạt độ Cs-134 đo được đạt 9,9 Bq/kg với độ sai khác so với giá trị thực nghiệm dưới 5,0%, khẳng định tính ưu việt của mô phỏng khi ma trận mẫu có tính đồng nhất cao.

Thứ ba, kết quả so sánh đối chứng giữa LabSOCS và mô phỏng GEANT4 kết hợp đo hệ số truyền qua trên mẫu đất tự nhiên (cấu hình trụ 76,00 mm) cho thấy sự tương đồng vượt trội: hoạt độ K-40 ghi nhận bởi LabSOCS là 139,9 Bq/kg, sai lệch so với mô phỏng GEANT4 chỉ dưới 4,5%.

Thảo luận kết quả

Các dữ liệu nghiên cứu được biểu diễn thông qua đường cong hiệu suất ghi (Efficiency Curve) phụ thuộc năng lượng và các bảng so sánh hoạt độ phóng xạ riêng. Đồ thị hiệu suất ghi phản ánh đúng quy luật vật lý: hiệu suất tăng nhanh từ 40 keV, đạt cực đại tại vùng 100-150 keV, sau đó giảm dần đều theo hàm logarit ở dải năng lượng cao từ 200 keV đến 1836 keV.

Nguyên nhân dẫn đến sai số ở vùng năng lượng thấp dưới 80 keV (như đỉnh 46,5 keV của Pb-210 với sai số 14-17% và 59,54 keV của Am-241) là do bề dày mẫu đo lên tới 3,0 cm gây ra hiện tượng tự hấp thụ mạnh mẽ, kết hợp với sự đóng góp phức tạp của bức xạ huỳnh quang tia X đặc trưng từ các nguyên tố cấu thành ma trận. Ngược lại, tại dải năng lượng từ 80 keV trở lên (như Cs-137 ở 661,6 keV hay K-40 ở 1460,8 keV), sai số đều nằm trong ngưỡng chấp nhận được (dưới 8,5%). Tất cả các kết quả phân tích kiểm tra tay nghề đều nằm trong khoảng giới hạn sai lệch tương đối cho phép (Maximum Acceptable Relative Bias) và được IAEA công nhận đạt chuẩn chất lượng quốc tế.

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

Thứ nhất, chuẩn hóa và áp dụng quy trình hiệu chuẩn không dùng nguồn LabSOCS vào các phép phân tích dịch vụ mẫu môi trường thường quy. Trung tâm Phân tích và các phòng thí nghiệm kiểm định bức xạ cần ban hành quy trình thao tác chuẩn trong vòng 6 tháng, đặt mục tiêu rút ngắn 90% thời gian xây dựng đường chuẩn và tiết kiệm 100% ngân sách mua mẫu chuẩn định kỳ 180 EUR mỗi mẫu.

Thứ hai, mở rộng và phát triển thư viện vật liệu chuyên biệt cho các mẫu môi trường đặc trưng tại Việt Nam. Viện Công nghệ Xạ hiếm và các cơ sở nghiên cứu cần tiến hành lập danh mục 25 loại vật liệu mới (bao gồm các loại đất đỏ bazan, quặng sa khoáng titan, tro bay nhiệt điện than và bùn lắng sông biển) với dải mật độ khối từ 0,60 đến 2,20 g/cm3 trên phần mềm Geometry Composer trong thời gian 12 tháng.

Thứ ba, tối ưu hóa cấu hình hình học hộp đo đối với các đồng vị phát gamma năng lượng thấp. Nhóm nghiên cứu thực nghiệm cần thiết kế và áp dụng các hộp đo có bề dày lớp mẫu mỏng dưới 1,0 cm (thay vì 3,0 cm như hiện nay) nhằm giảm thiểu hiệu ứng tự hấp thụ, giúp hạ mức sai số đo Pb-210 (46,5 keV) và Am-241 (59,54 keV) từ 17% xuống dưới mức 7,0% trong giai đoạn 2024-2025.

Thứ tư, xây dựng cầu nối liên thông tự động giữa phần mềm thương mại LabSOCS và mã nguồn mở GEANT4. Nhóm chuyên gia tính toán vật lý hạt nhân tại các trường đại học cần hoàn thiện module đối soát chéo dữ liệu hiệu suất trong 9 tháng, nâng cao độ tin cậy của phép đo đạt mức trên 99,5% cho các phân tích pháp y hạt nhân và an toàn bức xạ quốc gia.

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

Nhóm chuyên viên và kỹ thuật viên phân tích phóng xạ môi trường tại các trung tâm quan trắc: Tiếp cận trực tiếp quy trình thiết lập cấu hình đo trên hệ Genie 2000, giúp xử lý hàng trăm mẫu đất, nước, thực vật nhanh chóng mà không cần phụ thuộc vào nguồn chuẩn ngoại nhập.

Nhóm cán bộ nghiên cứu và quản lý an toàn bức xạ: Tham khảo hệ phương pháp đối chứng chi tiết giữa LabSOCS và GEANT4 để xây dựng phương án thẩm định thiết bị, tính toán che chắn bức xạ và đánh giá liều lượng phóng xạ cho cộng đồng.

Nhóm giảng viên, học viên cao học và sinh viên chuyên ngành Vật lý hạt nhân: Sử dụng luận văn như tài liệu học thuật chuyên sâu về tương tác bức xạ gamma với vật chất, kỹ thuật lập trình hướng đối tượng trong mô phỏng Monte-Carlo GEANT4 và các thuật toán làm khớp hiệu suất ghi.

Nhóm lãnh đạo quản lý các phòng thí nghiệm đạt chuẩn ISO/IEC 17025: Nắm bắt cơ sở khoa học để tối ưu hóa ngân sách đầu tư thiết bị, cắt giảm hàng trăm triệu đồng chi phí duy trì nguồn chuẩn và nâng cao năng lực cạnh tranh trong dịch vụ thử nghiệm hạt nhân.

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

Phần mềm LabSOCS hoạt động dựa trên cơ chế nào để loại bỏ nhu cầu sử dụng mẫu chuẩn vật lý? LabSOCS kết hợp thuật toán mô phỏng Monte-Carlo MCNP với cơ sở dữ liệu vật lý của NIST và các thông số kích thước tinh thể thực do nhà sản xuất detector cung cấp. Phần mềm tự động tính toán hiệu suất nạp đỉnh và hệ số tự hấp thụ dựa trên cấu hình hình học và mật độ khối được khai báo chỉ trong 120 giây.

Tại sao sai số của phương pháp LabSOCS lại tăng cao ở dải năng lượng dưới 80 keV? Ở mức năng lượng thấp dưới 80 keV (như Pb-210 tại 46,5 keV), tiết diện hiệu ứng quang điện tăng mạnh khiến khả năng tự hấp thụ gamma trong mẫu dày 3 cm trở nên rất lớn. Đồng thời, sự xuất hiện của các tia X đặc trưng từ thành phần nền ma trận gây nhiễu phổ, dẫn đến sai lệch khoảng 14% đến 17%.

Hiệu quả kinh tế mà phương pháp này mang lại cho phòng thí nghiệm là gì? Phương pháp giúp cắt giảm hoàn toàn chi phí mua các mẫu chuẩn từ IAEA với giá khoảng 180 EUR (gần 5.000.000 VND) cho mỗi chủng loại mẫu. Ngoài ra, cơ sở không còn phải chịu chi phí vận chuyển quốc tế phức tạp, chi phí lưu kho và xử lý chất thải phóng xạ khi nguồn chuẩn hết hạn sử dụng.

Mô phỏng Monte-Carlo GEANT4 đóng vai trò gì trong đề tài này? GEANT4 là công cụ mô phỏng độc lập viết bằng C++, cho phép tái lập chi tiết cấu trúc buồng đo và tính chất tương tác hạt nhân. Trong luận văn, GEANT4 kết hợp với phương pháp đo hệ số truyền qua được sử dụng làm chuẩn đối sánh để kiểm chứng tính đúng đắn và độ chính xác của kết quả tính từ phần mềm LabSOCS.

Phương pháp này có thể áp dụng cho mọi loại detector gamma thương mại không? LabSOCS yêu cầu tệp đặc tính hóa tinh thể chính xác từ nhà sản xuất, do đó chủ yếu tối ưu hóa cho các dòng đầu dò HPGe của Canberra/Mirion. Với các loại detector khác, người dùng cần kết hợp sử dụng mã nguồn mở GEANT4 hoặc các phần mềm như ANGLE để xây dựng mô hình mô phỏng tương đương.

Kết luận

• Luận văn đã làm chủ thành công quy trình ứng dụng phần mềm LabSOCS kết hợp mô phỏng Monte-Carlo GEANT4 trên hệ phổ kế gamma HPGe, mở ra giải pháp hiệu chuẩn không dùng nguồn toàn diện.
• Kết quả kiểm định trên 9 mẫu chuẩn IAEA cho thấy độ chính xác vượt trội với sai số dưới 5,6% trên mẫu đất, 1,1% trên mẫu thực vật và dưới 5,0% trên mẫu nước ở vùng năng lượng trên 80 keV.
• Nghiên cứu giải quyết triệt để bài toán kinh tế khi loại bỏ nhu cầu mua mẫu chuẩn 180 EUR/mẫu, đồng thời rút ngắn thời gian phân tích đường chuẩn xuống dưới 2 phút.
• Luận văn đóng góp trực tiếp vào công bố quốc tế về phương pháp xác định hoạt độ phóng xạ môi trường sử dụng một giá trị hiệu suất tuyệt đối kết hợp đường cong hiệu suất nội.
• Trong giai đoạn tới, cần tiếp tục mở rộng mô phỏng cho các cấu hình mẫu mỏng dưới 1,0 cm nhằm tối ưu hóa phép đo các đồng vị phóng xạ năng lượng thấp.

Quý độc giả và các nhà nghiên cứu quan tâm có thể khai thác toàn bộ phương pháp luận và số liệu thực nghiệm trong tài liệu để ứng dụng vào công tác phân tích phóng xạ và quản lý an toàn bức xạ chuyên nghiệp.