Khóa luận về phương pháp xác định thông số kỹ thuật của đầu dò NaI(Tl)

Luận văn tốt nghiệp kỹ thuật nghiên cứu nghiên cứu phương pháp mới xác định các thông số kỹ thuật của đầu dò nai tl, điều tra thực trạng, phân tích số liệu, đề xuất biện pháp cải

• Trường học: Trường Đại Học
• Chuyên ngành: Kỹ Thuật Hạt Nhân
• Đề tài: Nghiên Cứu Phương Pháp Xác Định Thông Số Kỹ Thuật Đầu Dò NaI(Tl)
• Loại tài liệu: Luận Văn
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỞ ĐẦU
• 1. TƯƠNG TÁC BỨC XẠ GAMMA VỚI VẬT CHẤT
• 1.1. Sự truyền bức xạ gamma qua vật chất
• 1.2. Hiệu ứng quang điện
• 1.3. Hiệu ứng Compton
• 1.4. Hiệu ứng tạo cặp
• 1.5. Các hiệu ứng xảy ra khi bức xạ truyền từ nguồn tới đầu dò
• 1.6. Phổ đo bức xạ gamma năng lượng 1408 keV
• 2. PHƯƠNG PHÁP MONTE CARLO VÀ CHƯƠNG TRÌNH MCNP5
• 2.1. Phương pháp Monte Carlo
• 2.2. Chương Trình MCNP5
• 2.3. Cấu trúc của một tập tin đầu vào (file input) trong chương trình MCNP5
• 2.3.1. Thẻ tiêu đề (Title Card)
• 2.3.2. Thẻ khai báo ô mạng (Cell Card)
• 2.3.3. Thẻ khai báo mặt (Surface Card)
• 2.3.4. Thẻ dữ liệu (Data Card)
• 2.3.4.1. Khai báo nguồn
• 3. ĐẦU DÒ NAI(TL) VÀ PHƯƠNG PHÁP XÁC ĐỊNH CÁC THÔNG SỐ TỐI ƯU CỦA ĐẦU DÒ NAI(TL)
• 3.1. Hiệu suất của đầu dò NaI(Tl)
• 3.2. Cấu hình và thông số kỹ thuật của Detector NaI(Tl)
• 3.3. Mô hình hóa hệ đo thực nghiệm trong mô phỏng MCNP5
• 3.4. Phương pháp xác định các thông số tối ưu của đầu dò NaI(Tl)
• 3.4.1. Phương pháp xác định mật độ tối ưu của lớp phản xạ Al2O3
• 3.4.2. Phương pháp xác định bán kính tối ưu của tinh thể NaI(Tl)
• 3.4.3. Phương pháp xác định chiều dài tối ưu của tinh thể NaI(Tl)
• 4. KẾT QUẢ VÀ THẢO LUẬN
• 4.1. Kết quả xác định mật độ của lớp phản xạ
• 4.2. Kết quả xác định bán kính tối ưu của tinh thể NaI(Tl)
• 4.3. Kết quả xác định chiều dài tối ưu của tinh thể NaI(Tl)

Tóm tắt

I. Tổng quan về phương pháp xác định thông số kỹ thuật đầu dò NaI Tl

Đầu dò NaI(Tl) là một trong những thiết bị quan trọng trong lĩnh vực đo lường bức xạ gamma. Việc xác định các thông số kỹ thuật của đầu dò này là cần thiết để đảm bảo hiệu suất hoạt động tối ưu. Phương pháp xác định thông số kỹ thuật đầu dò NaI(Tl) không chỉ giúp cải thiện độ chính xác trong các phép đo mà còn nâng cao hiệu suất ghi nhận bức xạ. Bài viết này sẽ trình bày chi tiết về các phương pháp xác định thông số kỹ thuật của đầu dò NaI(Tl) và ứng dụng của chúng trong thực tiễn.

1.1. Giới thiệu về đầu dò NaI Tl và ứng dụng

Đầu dò NaI(Tl) là một loại đầu dò gamma sử dụng tinh thể natri iodide doped thallium. Thiết bị này được sử dụng rộng rãi trong các ứng dụng y tế, nghiên cứu hạt nhân và bảo vệ môi trường. Đầu dò NaI(Tl) có khả năng phát hiện bức xạ gamma với độ nhạy cao, giúp người dùng theo dõi hoạt động phóng xạ trong môi trường.

1.2. Tầm quan trọng của việc xác định thông số kỹ thuật

Việc xác định thông số kỹ thuật của đầu dò NaI(Tl) là rất quan trọng để đảm bảo độ chính xác trong các phép đo bức xạ. Các thông số như độ nhạy, độ phân giải năng lượng và hiệu suất ghi nhận ảnh hưởng trực tiếp đến kết quả đo. Do đó, việc tối ưu hóa các thông số này là cần thiết để nâng cao hiệu quả sử dụng đầu dò.

II. Các thách thức trong việc xác định thông số kỹ thuật đầu dò NaI Tl

Mặc dù đầu dò NaI(Tl) có nhiều ưu điểm, nhưng việc xác định thông số kỹ thuật của nó cũng gặp phải nhiều thách thức. Các yếu tố như độ dày lớp phản xạ, bán kính tinh thể và chiều dài tinh thể đều có thể ảnh hưởng đến hiệu suất ghi nhận. Việc hiểu rõ các thách thức này sẽ giúp cải thiện quy trình xác định thông số kỹ thuật.

2.1. Ảnh hưởng của lớp phản xạ đến hiệu suất

Lớp phản xạ bao quanh tinh thể NaI(Tl) có thể ảnh hưởng lớn đến hiệu suất ghi nhận bức xạ. Nghiên cứu cho thấy rằng sự thay đổi độ dày của lớp phản xạ Al2O3 có thể làm thay đổi hiệu suất ghi nhận, do đó cần phải xác định mật độ tối ưu của lớp phản xạ này.

2.2. Tác động của bán kính và chiều dài tinh thể

Bán kính và chiều dài của tinh thể NaI(Tl) cũng là những yếu tố quan trọng ảnh hưởng đến hiệu suất ghi nhận. Việc xác định các thông số này một cách chính xác sẽ giúp tối ưu hóa hiệu suất của đầu dò, từ đó nâng cao độ chính xác trong các phép đo bức xạ.

III. Phương pháp xác định thông số kỹ thuật đầu dò NaI Tl

Có nhiều phương pháp khác nhau để xác định thông số kỹ thuật của đầu dò NaI(Tl). Trong đó, phương pháp Monte Carlo và mô phỏng MCNP5 là hai phương pháp phổ biến nhất. Những phương pháp này cho phép mô phỏng quá trình tương tác giữa bức xạ gamma và đầu dò, từ đó xác định các thông số kỹ thuật một cách chính xác.

3.1. Phương pháp Monte Carlo

Phương pháp Monte Carlo là một kỹ thuật mô phỏng dựa trên việc sử dụng các phép thử ngẫu nhiên để giải quyết các bài toán phức tạp. Phương pháp này đã được áp dụng thành công trong việc xác định thông số kỹ thuật của đầu dò NaI(Tl), giúp cải thiện độ chính xác trong các phép đo.

3.2. Mô phỏng MCNP5

MCNP5 là một phần mềm mô phỏng mạnh mẽ được sử dụng để mô phỏng quá trình vận chuyển hạt. Phần mềm này cho phép người dùng xác định các thông số kỹ thuật của đầu dò NaI(Tl) thông qua việc mô phỏng các tương tác giữa bức xạ gamma và đầu dò, từ đó đưa ra các kết quả chính xác.

IV. Kết quả nghiên cứu và ứng dụng thực tiễn

Kết quả từ các nghiên cứu về thông số kỹ thuật của đầu dò NaI(Tl) đã chỉ ra rằng việc tối ưu hóa các thông số này có thể nâng cao hiệu suất ghi nhận bức xạ. Các ứng dụng thực tiễn của đầu dò NaI(Tl) trong lĩnh vực y tế và bảo vệ môi trường đã chứng minh được giá trị của việc xác định thông số kỹ thuật một cách chính xác.

4.1. Kết quả xác định mật độ lớp phản xạ

Nghiên cứu đã chỉ ra rằng mật độ lớp phản xạ Al2O3 tối ưu có thể cải thiện hiệu suất ghi nhận bức xạ gamma. Kết quả cho thấy sự phù hợp giữa mô phỏng và thực nghiệm với độ lệch dưới 2% ở các mức năng lượng khác nhau.

4.2. Ứng dụng trong lĩnh vực y tế

Đầu dò NaI(Tl) được sử dụng rộng rãi trong lĩnh vực y tế để phát hiện và đo lường bức xạ gamma. Việc xác định thông số kỹ thuật chính xác giúp nâng cao độ chính xác trong các phép đo, từ đó hỗ trợ hiệu quả trong chẩn đoán và điều trị.

V. Kết luận và tương lai của nghiên cứu

Việc xác định thông số kỹ thuật của đầu dò NaI(Tl) là một lĩnh vực nghiên cứu quan trọng trong ngành đo lường bức xạ. Các phương pháp hiện tại như Monte Carlo và mô phỏng MCNP5 đã cho thấy hiệu quả trong việc tối ưu hóa các thông số này. Tương lai của nghiên cứu sẽ tiếp tục tập trung vào việc cải thiện độ chính xác và hiệu suất của đầu dò NaI(Tl).

5.1. Tương lai của đầu dò NaI Tl

Đầu dò NaI(Tl) sẽ tiếp tục được cải tiến về mặt công nghệ để nâng cao hiệu suất và độ nhạy. Các nghiên cứu mới sẽ tập trung vào việc phát triển các vật liệu mới và cải tiến quy trình sản xuất để đáp ứng nhu cầu ngày càng cao trong lĩnh vực đo lường bức xạ.

5.2. Hướng nghiên cứu tiếp theo

Các nghiên cứu tiếp theo sẽ tập trung vào việc phát triển các phương pháp mới để xác định thông số kỹ thuật của đầu dò NaI(Tl). Việc áp dụng công nghệ mới và các phương pháp mô phỏng tiên tiến sẽ giúp nâng cao độ chính xác và hiệu suất của đầu dò trong tương lai.

Nội dung chính

Tổng quan nghiên cứu

Kỹ thuật hạt nhân đóng vai trò then chốt trong y tế, năng lượng và quan trắc môi trường. Trong các hệ thống đo lường bức xạ hạt nhân, hệ phổ kế gamma sử dụng đầu dò nhấp nháy NaI(Tl) kết hợp bộ phân tích đa kênh MCA là công cụ phổ biến nhất để định tính và định lượng các đồng vị phóng xạ. Tuy nhiên, việc xác định chính xác hiệu suất đỉnh năng lượng toàn phần (FEPE) của đầu dò luôn là một thách thức kỹ thuật lớn. Hiệu suất ghi bức xạ thực tế không bao giờ đạt 100% do sự suy giảm cường độ photon khi truyền qua các lớp vỏ bọc bảo vệ, lớp đệm không khí và lớp vật liệu phản xạ.

Khi ứng dụng phần mềm mô phỏng MCNP5 dựa trên phương pháp Monte Carlo để tính toán hiệu suất ghi nhận, các nhà nghiên cứu thường gặp phải sự sai lệch đáng kể giữa kết quả mô phỏng và số liệu thực nghiệm. Việc sử dụng trực tiếp các thông số kỹ thuật do nhà sản xuất cung cấp thường dẫn đến độ lệch rất cao, điển hình như mức sai số lên tới 19,51% tại đỉnh năng lượng 31 keV của nguồn Ba-133 hoặc 13,60% tại đỉnh 1408 keV của nguồn Eu-152.

Mục tiêu trọng tâm của nghiên cứu là xây dựng quy trình độc lập nhằm hiệu chỉnh và xác định chính xác ba thông số cấu trúc của đầu dò NaI(Tl) kích thước 76,2 x 76,2 mm: mật độ lớp phản xạ nhôm oxit ($Al_2O_3$), bán kính tinh thể và chiều dài tinh thể. Nghiên cứu được triển khai trong phạm vi dải năng lượng rộng từ 31 keV đến 1408 keV với các nguồn phóng xạ chuẩn. Kết quả tối ưu hóa đã đưa độ lệch mô phỏng và thực nghiệm từ mức 19,51% giảm xuống còn 0,19% tại năng lượng 31 keV và đạt mức sai số dưới 0,03% tại năng lượng 662 keV, tạo cơ sở khoa học vững chắc cho việc hiệu chuẩn phổ kế hạt nhân không cần nguồn chuẩn phức tạp.

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 tương tác giữa bức xạ photon với vật chất và cơ chế phát quang nhấp nháy. Sự suy giảm cường độ chùm tia gamma đơn sắc khi truyền qua vật chất có bề dày $d$ tuân theo định luật hàm mũ $I = I_0 \cdot e^{-\mu d}$, trong đó $\mu$ là hệ số suy giảm tuyến tính phụ thuộc trực tiếp vào mật độ vật liệu. Ba cơ chế tương tác cơ bản của photon bao gồm:

1. Hiệu ứng quang điện: Photon truyền toàn bộ năng lượng cho electron liên kết ở lớp vỏ nguyên tử theo phương trình bảo toàn năng lượng với động năng cực đại $E_k = hf - I_{lk}$. Hiệu ứng này chiếm ưu thế tuyệt đối ở dải năng lượng thấp dưới 100 keV.
2. Tán xạ Compton: Sự va chạm giữa photon năng lượng trung bình với electron xem như tự do, làm photon bị tán xạ góc $\theta$ và mất một phần năng lượng theo biểu thức $T_e = E - E'$.
3. Hiệu ứng tạo cặp: Xảy ra khi photon có năng lượng lớn hơn ngưỡng 1022 keV tương tác với điện trường hạt nhân, biến đổi thành một cặp electron - positron, sau đó quá trình hủy hạt giải phóng hai photon năng lượng 511 keV.

Bên cạnh đó, lý thuyết vùng năng lượng chất rắn giải thích cơ chế tinh thể NaI pha tạp ion Thallium (Tl) tạo ra các mức năng lượng kích hoạt trung gian giữa vùng hóa trị và vùng dẫn. Quá trình tái tổ hợp của cặp electron - lỗ trống tại tâm kích hoạt phát ra các photon ánh sáng khả kiến tỉ lệ thuận với năng lượng hấp thụ. Hiệu suất đỉnh năng lượng toàn phần (FEPE) được xác định bằng tích của hiệu suất nội và hiệu suất hình học thông qua góc khối $\Omega$.

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

Phương pháp mô phỏng Monte Carlo vận chuyển hạt được thực hiện thông qua chương trình MCNP5 viết trên nền Fortran 90. Cấu trúc tệp đầu vào (file input) thiết lập chi tiết cấu hình hình học đa tầng của đầu dò gồm: lớp vỏ nhôm dày 0,1 cm (mật độ 2,699 g/cm³), lớp đệm không khí dày 0,1 cm, lớp phản xạ $Al_2O_3$ dày 0,16 cm và tinh thể nhấp nháy NaI(Tl). Tally F8 được chỉ định tại ô mạng chứa tinh thể để ghi nhận năng lượng photon thất thoát và tích lũy phổ xung.

Nghiên cứu triển khai hệ thống ba chuỗi thực nghiệm và mô phỏng tương ứng:

• Chuỗi 1: Nguồn Ba-133 hoạt độ 37.333 Bq được đặt trong ống chuẩn trực bằng đồng, phát chùm tia chuẩn trực 31 keV để khảo sát riêng biệt mật độ lớp phản xạ $Al_2O_3$ với các mức mô phỏng từ $10^5$ đến $10^6$ hạt photon, thời gian đo thực nghiệm đạt 25.500 giây và 89.600 giây.
• Chuỗi 2: Sử dụng các nguồn điểm phát photon năng lượng thấp (31, 32, 59, 81, 121 keV) gồm Ba-133, Cs-137, Am-241, Eu-152 đặt cách đầu dò 40 cm với số lượng hạt mô phỏng $10^7$ photon để tối ưu bán kính tinh thể.
• Chuỗi 3: Sử dụng các nguồn phát photon năng lượng cao (662, 964, 1173, 1274, 1332, 1408 keV) gồm Cs-137, Na-22, Co-60, Eu-152 đặt cách đầu dò 40 cm để tối ưu chiều dài tinh thể.

Khoảng cách từ nguồn đến đầu dò được điều khiển tự động bằng hệ cơ khí chính xác với sai số dịch chuyển chỉ 0,01 mm. Phương pháp phân tích sử dụng khai triển chuỗi MacLaurin bậc nhất cho hàm truyền qua bức xạ, kết hợp thuật toán khớp hàm tuyến tính bình phương tối thiểu có trọng số và phương pháp tính trung bình trọng số để xác định các giá trị hình học tối ưu.

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 dữ liệu mô phỏng và đối chiếu thực nghiệm đã mang lại bốn phát hiện khoa học quan trọng:

1. Mật độ lớp phản xạ $Al_2O_3$ thực tế cao hơn đáng kể so với thông số danh định: Dữ liệu khớp hàm tuyến tính giữa hiệu suất đỉnh 31 keV và mật độ lớp phản xạ cho hệ số tương quan $R^2 = 0,9978$ (nguồn đặt trước) và $R^2 = 0,9977$ (nguồn đặt bên). Giá trị mật độ lớp phản xạ nội suy từ thực nghiệm đạt trung bình 2,02 ± 0,32 g/cm³, chênh lệch gấp 3,67 lần so với thông số 0,55 g/cm³ do nhà sản xuất cung cấp.
2. Bán kính tối ưu của tinh thể NaI(Tl) nhỏ hơn kích thước thiết kế ban đầu: Khảo sát trên 5 đỉnh năng lượng thấp từ 31 keV đến 121 keV cho thấy các hàm khớp tuyến tính theo bán kính đều đạt độ tin cậy $R^2$ từ 0,9990 đến 0,9999. Bằng phương pháp trung bình trọng số, bán kính tối ưu được xác định là 3,76 ± 0,03 cm (thay vì 3,81 cm ban đầu).
3. Chiều dài hiệu dụng của tinh thể NaI(Tl) đạt 7,54 cm: Dữ liệu từ các đỉnh năng lượng cao từ 662 keV đến 1408 keV cho thấy tương quan tuyến tính rất chặt chẽ giữa chiều dài và hiệu suất ghi nhận ($R^2$ đạt từ 0,9978 đến 0,9999). Chiều dài tối ưu được xác định chính xác tại mức 7,54 ± 0,589 cm, giảm 0,08 cm so với thông số nhà sản xuất là 7,62 cm.
4. Độ lệch tổng thể giữa mô phỏng và thực nghiệm giảm triệt để trên toàn dải năng lượng: Khi áp dụng đồng bộ ba thông số tối ưu, sai số tại đỉnh 31 keV giảm từ 19,51% xuống còn 0,19%; đỉnh 32 keV giảm từ 18,87% xuống 0,25%; đỉnh 81 keV giảm từ 3,96% xuống 0,03%; đỉnh 662 keV giảm từ 0,48% xuống 0,03%; đỉnh 1173 keV giảm từ 8,95% xuống 1,91%.

Thảo luận kết quả

Cơ chế vật lý đằng sau sự phân hóa độ nhạy của các thông số hình học được giải thích thông qua đường đi tự do trung bình và độ sâu đâm xuyên của photon. Đối với photon năng lượng thấp như 31 keV, toàn bộ quá trình tương tác quang điện diễn ra ngay tại lớp bề mặt tinh thể. Do đó, thông số mật độ lớp bảo vệ $Al_2O_3$ và bán kính bề mặt đóng vai trò quyết định, trong khi chiều dài tinh thể hoàn toàn không ảnh hưởng đến số đếm đỉnh. Ngược lại, photon năng lượng cao trên 662 keV có xác suất đâm xuyên sâu xuống đáy tinh thể, khiến chiều dài tinh thể trở thành yếu tố chi phối trực tiếp đến hiệu suất thu nhận.

Dữ liệu nghiên cứu được biểu diễn trực quan qua các đồ thị hàm hồi quy tuyến tính của hiệu suất đỉnh theo mật độ và kích thước hình học, kết hợp bảng ma trận so sánh sai số tương đối giữa mô hình ban đầu và mô hình sau tối ưu. So với các công trình trước đây như nghiên cứu của Tam và cộng sự (chỉ tập trung thay đổi bề dày lớp phản xạ $Al_2O_3$), hướng tiếp cận của luận văn đã tách biệt thành công ảnh hưởng độc lập của từng biến số cấu trúc. Việc hiệu chỉnh đồng thời mật độ $Al_2O_3$ (2,02 g/cm³), bán kính (3,76 cm) và chiều dài (7,54 cm) giúp mô hình MCNP5 đạt độ tương thích hoàn hảo với hệ phổ kế thực nghiệm, duy trì sai số dưới 2,82% trên hầu hết các đỉnh năng lượng khảo sát.

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

Dựa trên kết quả nghiên cứu, bốn giải pháp kỹ thuật và khuyến nghị thực tiễn được đề xuất nhằm nâng cao độ chính xác cho các phép đo phổ kế hạt nhân:

1. Cập nhật tham số hình học chuẩn hóa vào các tệp đầu vào MCNP5: Các phòng thí nghiệm và cơ sở phân tích phóng xạ cần thay thế thông số danh định của đầu dò NaI(Tl) 76,2 x 76,2 mm bằng bộ thông số tối ưu gồm mật độ lớp phản xạ $Al_2O_3$ bằng 2,02 g/cm³, bán kính tinh thể 3,76 cm và chiều dài tinh thể 7,54 cm. Giải pháp này giúp đưa sai số mô phỏng về dưới ngưỡng 3% mà không cần tốn chi phí thực nghiệm lặp lại. Thời gian áp dụng có thể thực hiện ngay lập tức bởi các kỹ sư vận hành hệ phổ kế.
2. Chuẩn hóa quy trình tách biến độc lập trong hiệu chuẩn đầu dò: Các viện nghiên cứu bức xạ hạt nhân nên áp dụng quy trình chuẩn trực nguồn năng lượng thấp (dưới 35 keV) để khóa ảnh hưởng của kích thước tinh thể khi xác định lớp vỏ phản xạ, sau đó mới dùng nguồn dải rộng để xác định bán kính và chiều dài. Mục tiêu là kiểm soát độ không đảm bảo đo vị trí dưới 0,01 mm trong vòng 3 tháng tới.
3. Xây dựng cơ sở dữ liệu số hóa cho các dòng đầu dò nhấp nháy thương mại: Đề xuất các cơ quan quản lý an toàn bức xạ và viện nghiên cứu chuyên ngành chủ trì dự án đo đạc, xây dựng thư viện mô hình MCNP5 tối ưu cho các kích thước đầu dò phổ biến như 50,8 x 50,8 mm hoặc 127 x 127 mm. Dự án thực hiện theo lộ trình 6 đến 12 tháng với mục tiêu chuẩn hóa 100% các dòng detector tại Việt Nam.
4. Nâng cao số lượng lịch sử hạt trong mô phỏng vận chuyển bức xạ: Khuyến nghị người phân tích luôn thiết lập số hạt mô phỏng tối thiểu đạt từ $10^7$ photon trở lên trong chương trình MCNP5 đối với các bài toán xác định hiệu suất đỉnh, nhằm đảm bảo sai số thống kê của Tally F8 luôn nhỏ hơn 0,01%.

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

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho bốn nhóm đối tượng chính sau:

• Nghiên cứu viên và kỹ sư vật lý hạt nhân: Cung cấp phương pháp luận chặt chẽ và thuật toán tối ưu hóa thông số hình học cho các hệ detector nhấp nháy, hỗ trợ tính toán hiệu suất ghi nhận chính xác cho các phép phân tích kích hoạt nơtron hoặc quan trắc bức xạ môi trường.
• Cán bộ kỹ thuật tại các trung tâm an toàn bức xạ và y học hạt nhân: Cung cấp bộ thông số chuẩn đã được kiểm chứng thực nghiệm, giúp xây dựng các mô hình mô phỏng liều lượng bức xạ và kiểm định hoạt độ phóng xạ mẫu môi trường với độ tin cậy cao trên dải từ 31 keV đến 1408 keV.
• Giảng viên, học viên cao học và sinh viên chuyên ngành vật lý: Đóng vai trò là giáo trình thực hành mẫu về cách thức xây dựng tệp đầu vào MCNP5 (khai báo Cell, Surface, Data, Tally F8) và kỹ thuật xử lý số liệu thực nghiệm bằng các phương pháp toán học thống kê.
• Nhà sản xuất và đơn vị kiểm định thiết bị đo bức xạ: Giúp nhận diện nguyên nhân gây ra sự sai lệch giữa thông số thiết kế xuất xưởng và đặc tính vật lý thực tế của tinh thể sau quá trình gia công, đóng gói và lão hóa vật liệu.

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

1. Tại sao có sự chênh lệch lớn giữa mật độ lớp phản xạ thực tế (2,02 g/cm³) và thông số nhà sản xuất (0,55 g/cm³)?
Thông số 0,55 g/cm³ do nhà sản xuất cung cấp thường là mật độ khối biểu kiến của bột nhôm oxit khi chưa nén chặt. Trong quá trình đóng gói cơ học để bọc kín tinh thể NaI(Tl), lớp bột phản xạ bị nén lại dưới áp lực cao, làm mật độ thực tế tăng lên mức 2,02 ± 0,32 g/cm³.

2. Làm thế nào để tách rời sự ảnh hưởng của mật độ lớp phản xạ khỏi các kích thước tinh thể?
Nghiên cứu sử dụng nguồn Ba-133 phát bức xạ 31 keV đặt trong ống chuẩn trực bằng đồng. Chùm photon bị giới hạn hẹp ở tâm nên không phụ thuộc vào bán kính đầu dò. Năng lượng 31 keV hấp thụ hoàn toàn ở bề mặt nên không bị ảnh hưởng bởi chiều dài tinh thể.

3. Tại sao bán kính tinh thể tối ưu (3,76 cm) lại nhỏ hơn thông số danh định (3,81 cm)?
Sự suy giảm bán kính hiệu dụng xuất phát từ hiện tượng lớp chết (dead layer) hoặc sự suy giảm khả năng thu nhận ánh sáng nhấp nháy tại các mép biên của tinh thể. Các photon tương tác gần viền ngoài không tạo ra tín hiệu quang điện trọn vẹn, dẫn đến bán kính thu nhận hiệu dụng giảm đi khoảng 0,05 cm.

4. Hiệu quả cải thiện sai số sau khi tối ưu đồng bộ cả ba thông số đạt mức nào?
Mô hình tối ưu giúp giảm triệt để sai số trên toàn phổ đo. Cụ thể, độ lệch tại đỉnh 31 keV giảm từ 19,51% xuống 0,19%, đỉnh 32 keV giảm từ 18,87% xuống 0,25%, đỉnh 662 keV giảm từ 0,48% xuống 0,03%, và các đỉnh trên 1000 keV đều duy trì sai số dưới mức chấp nhận học thuật.

5. Phương pháp này có thể áp dụng cho các loại đầu dò nhấp nháy khác như LaBr3(Ce) hay BGO không?
Hoàn toàn có thể áp dụng. Quy trình phân tách ba bước (chuẩn trực năng lượng thấp xác định lớp vỏ, nguồn dải thấp xác định bán kính, nguồn dải cao xác định chiều dài) là phương pháp tổng quát cho mọi hệ đo nhấp nháy hình trụ khi mô phỏng bằng MCNP5 hoặc Geant4.

Kết luận

• Luận văn đã giải quyết triệt để bài toán sai lệch hiệu suất ghi nhận giữa mô phỏng MCNP5 và thực nghiệm trên hệ phổ kế gamma NaI(Tl).
• Xác định chính xác bộ ba thông số thực tế của đầu dò 76,2 x 76,2 mm: mật độ $Al_2O_3$ đạt 2,02 ± 0,32 g/cm³, bán kính tinh thể đạt 3,76 ± 0,03 cm, chiều dài tinh thể đạt 7,54 ± 0,589 cm.
• Giảm độ lệch hiệu suất đỉnh năng lượng toàn phần từ mức cực đại 19,51% xuống còn 0,19% tại vùng năng lượng thấp và đạt mức 0,03% tại năng lượng 662 keV.
• Xây dựng thành công quy trình chuẩn hóa bán thực nghiệm có khả năng ứng dụng rộng rãi trong tính toán hiệu suất mà không cần nguồn chuẩn phóng xạ đắt tiền.
• Lộ trình phát triển tiếp theo bao gồm việc mở rộng thuật toán tối ưu hóa cho các dạng hình học nguồn thể tích và tích hợp vào các hệ thống quan trắc phóng xạ môi trường tự động.

Các nhà nghiên cứu và cơ sở phân tích hạt nhân quan tâm có thể áp dụng trực tiếp bộ tham số này vào các mô hình tính toán phổ kế để tối ưu hóa độ chính xác của các phép phân tích định lượng bức xạ.

Trích đoạn nội dung tài liệu

Mở đầu Ngày nay, nhiều kỹ thuật hạt nhân được ứng dụng vào đời sống đặc biệt là những kỹ thuật được ứng dụng trong các lĩnh vực y tế, năng lượng, môi trường. Việc bắt đầu sử dụng nguồn phóng xạ làm ảnh hưởng đến sức khỏe của người vận hành. Các máy đo phóng xạ môi trường trở thành những công cụ cơ bản cho phép người sử dụng kiểm tra về hoạt độ phóng xạ từ môi trường. Hệ phổ kế gamma là một trong những hệ thống phát hiện bức xạ được sử dụng rộng rãi nhất.

Trong phép đo phóng xạ cần có kiến thức chính xác về hiệu suất ghi của máy đo bởi chỉ một phần của bức xạ đi vào vật liệu tương tác bên trong nên hiệu suất ghi không đạt 100%. Một trong những đầu dò có hiệu suất cao để đo hoạt độ môi trường là đầu dò sử dụng chất nhấp nháy rắn điển hình như hệ đo phổ gamma NaI(Tl) bao gồm một đầu dò NaI(Tl) và máy phân tích đa kênh MCA, hiệu suất ghi bức xạ phụ thuộc vào nhiều tham số của đầu dò. Nhiều phần mềm đã phát triển rất sớm từ những năm 1940 cho phép người dùng tính toán phù hợp với mô hình thực nghiệm mà không cần làm việc trực tiếp với nguồn phóng xạ. Phần mềm mô phỏng MCNP5 dựa trên phương pháp Monte Carlo được xây dựng bởi các nhà khoa học tại phòng thí nghiệm quốc gia Los Alamos, MCNP5 được nhiều nhà khoa học trên thế giới sử dụng vì sự phù hợp của mô phỏng so với thực nghiệm.

Phần mềm MCNP5 cho phép người sử dụng mô phỏng lại quá trình vận chuyển hạt từ những dữ liệu đầu vào của mô hình thực nghiệm, trong mô phỏng để tính được hiệu suất ghi của đầu dò cần xác định được số hạt để lại năng lượng trong tinh thể. Trong quá trình photon phát ra từ nguồn trên đường đi nó phải qua các vật liệu môi trường và các lớp che chắn tinh thể. Hiệu suất đỉnh năng lượng hấp thụ toàn phần được tính giữa thực nghiệm và mô phỏng bao giờ cũng có sự chệnh lệch tùy thuộc vào các thông số đầu vào. Sự phù hợp giữa tính toán hiệu suất từ mô phỏng và thực nghiệm cần có những nghiên cứu liên quan giữa các thông số đầu vào đối với kết quả tính toán.

Khi tính hiệu suất bằng phương pháp gamma truyền qua thì các yếu tố chính ảnh hưởng đến kết quả là các thông số của đầu dò NaI(Tl) được cung cấp từ nhà sản 1 xuất, việc hiệu chỉnh các yếu tố này trước khi mô phỏng sẽ cho kết quả tối ưu hơn. Sự ảnh hưởng của lớp phản xạ bao quanh tinh thể đã được nghiên cứu bởi Tam và cộng sự [5]. Kết quả mô phỏng cho thấy khi thay đổi bề dày lớp phản xạ hiệu suất cũng thay đổi và phụ thuộc tuyến tính theo bề dày lớp phản xạ Al2O3 , sự hiệu chỉnh thông số này cho thấy sự phù hợp với độ lệch dưới 2% giữa kết quả mô phỏng và thực nghiệm ở các mức năng lượng trải dài từ 88 keV- 1332 keV. Thay vì hiệu chỉnh thông số bề dày lớp phản xạ chúng tôi hiệu chỉnh các thông số khác, đồng thời đưa ra phương pháp xây dựng một quy trình để tách rời sự ảnh hưởng của từng thông số lên bài toán, sau đó đưa ra mô hình tối ưu giữa mô phỏng và thực nghiệm.

Theo những nội dung trên nên luận văn được chia thành bốn chương. Chương một trình bày cơ sở lý thuyết về tương tác giữa bức xạ gamma và vật chất, những tương các cơ bản như: quang điện, Compton và tạo cặp. Chương hai giới thiệu về phương pháp Monte Carlo và chương trình MCNP5. Chương ba trình bày về mô hình thực nghiệm và phương pháp xác định các thông số tối ưu.

Trong chương này nghiên cứu về ba phương pháp để xác định lần lượt các thông số mật độ lớp phản xạ, bán kính tinh thể NaI(Tl), và chiều dài tinh thể NaI(Tl). Chương bốn sẽ trình bày về kết quả của các thông số tối ưu của đầu dò NaI(Tl) thu được đối với từng phương pháp, từ kết quả thu được sẽ thay lại các thông số này vào mô phỏng, sau đó tiến hành so sánh hiệu suất đỉnh năng lượng hấp thụ toàn phần giữa mô phỏng và thực nghiệm và thảo luận về các thông số tối ưu với mô hình mới. TƯƠNG TÁC BỨC XẠ GAMMA VỚI VẬT CHẤT 1. Sự truyền bức xạ gamma qua vật chất Bản chất của bức xạ gamma là sóng điện từ mang năng lượng cao ứng với bước sóng nhỏ hơn 10 11 m.

Bức xạ gamma thực chất là các hạt photon, chúng có tính chất của cả sóng và hạt, khi đi vào vật liệu photon tương tác với các electron, thường xảy ra ba hiệu ứng: quang điện, Compton và tạo cặp. Do xảy ra tương tác giữa photon và electron, khi truyền qua vật liệu bia một phần cường độ của chùm tia bị suy giảm, vì vậy số đếm photon suy giảm về số lượng tùy thuộc vào độ dày vật liệu bia và năng lượng photon tới. Quy luật suy giảm của cường độ chùm tia photon đi qua vật liệu được tính theo công thức: I  I 0 e   d photon.2) trong đó: d (cm) là bề dày vật liệu.  ( cm 1 ) là hệ số suy giảm tuyến tính đối với vật liệu bia, do hệ số tuyến tính phụ thuộc vào mật độ bia nên người ta thường sử dụng hệ số suy giảm khối để mô tả sự suy giảm của cường độ photon truyền qua.cm3 ) là mật độ của bia.

Hiệu ứng quang điện Thí nghiệm nổi tiếng của Heinrich Hertz vào năm 1887 là một trong những điều kỳ lạ trong lịch sử khoa học, ông phát hiện ra sóng điện từ xác nhận lý thuyết sóng của James Maxwell, ông cũng là người khám phá ra hiệu ứng quang điện dẫn đến tính chất hạt ánh sáng [6]. Abert Einstein dựa vào lý thuyết lượng tử năng lượng của Max Plank đã giải thích thành công hiện tượng quang điện. Giả thuyết photon mang năng lượng lớn hơn năng lượng liên kết của electron với hạt nhân nguyên tử đi vào kim loại truyền hết 3 năng lượng cho electron, theo định luật bảo toàn năng lượng thì động năng cực đại của electron bứt ra khỏi bề mặt kim loại bằng hiệu năng lượng photon tới và năng lượng liên kết của electron với hạt nhân nguyên tử. Đường cong năng lượng của electron trên bề mặt kim loại, một electron ở lớp vỏ ngoài cùng hấp thụ một photon năng lượng bật ra khỏi kim loại [6].

1 2 mvmax  hf   lk (1.4) 2 trong đó: 1 2 mvmax là động năng cực đại của electron. 2 hf là năng lượng của photon tới.  lk là năng lượng liên kết giữa electron và hạt nhân. Năng lượng liên kết của electron giảm dần theo các lớp K, L, M, N… Nếu năng lượng của photon nhỏ hơn năng lượng liên kết của electron ở lớp K thì hiệu ứng quang điện chỉ xảy ra cho các electron ở lớp xa hạt nhân hơn.

Mỗi một nguyên tử có cấu trúc năng lượng ở các lớp vỏ electron khác nhau nên xác suất xảy ra hiệu ứng quang điện 4 không những phụ thuộc vào năng lượng của photon tới mà còn phụ thuộc vào số hiệu nguyên tử. Hiệu ứng Compton Hiệu ứng Compton là sự va chạm giữa photon và electron tự do, trong thực tế thì electron không tự do mà là những electron liên kết với hạt nhân trong nguyên tử môi trường. Đối với photon đi vào môi trường có năng lượng lớn hơn nhiều so với năng lượng liên kết của electron ta có thể bỏ qua năng lượng liên kết của electron và xem như bài toán va chạm giữa photon với các electron tự do. Hiệu ứng Compton[2].

Sự va chạm giữa photon với các electron ở lớp ngoài cùng của nguyên tử (xem như electron tự do) được gọi là tán xạ Compton. Sau tán xạ photon thay đổi phương chuyển động và bị mất một phần năng lượng còn electron được giải phóng ra khỏi nguyên tử. Theo định luật bảo toàn năng lượng: Te  E  E ' (1.5) Trong đó: Te là động năng cực đại của electron. E là năng lượng của photon tới.

E ' là năng lượng của photon sau tán xạ Compton. Trên cơ sở tính toán dựa trên bảo toàn năng lượng và động lượng có xét đến tương đối tính, năng lượng của photon và electron sau khi tán xạ theo góc  lần lượt được tính theo công thức (1. Hiệu ứng tạo cặp Hiệu ứng tạo cặp là sự tương tác giữa một photon có năng lượng lớn hơn 1022 keV với hạt nhân nguyên tử, kết quả là sự biến mất của photon cùng với sự xuất hiện cặp electron và positron. Positron sau đó nhanh chóng bị hủy do tương tác với các electron khác từ môi trường và sinh ra hai photon với năng lượng 511 keV.3 Hiệu ứng tạo cặp [2].

Theo định luật bảo toàn năng lượng: Te  Te  E  2(mec 2 )  E  1022 KeV (1.8) Hiêu ứng tạo cặp xảy ra bên ngoài tinh thể tạo thành photon được ghi nhận ở đỉnh năng lượng 511 keV, nếu hiệu ứng tạo cặp xảy ra bên trong tinh thể và các photon được tạo thành sau đó thoát ra thì sẽ ghi nhận các đỉnh thoát đơn hoặc thoát đôi. Các hiệu ứng xảy ra khi bức xạ truyền từ nguồn tới đầu dò [2]. Trong ghi nhận bức xạ từ phổ gamma các hiệu ứng quang điện, Compton và tạo cặp được thể hiện trong hình 1. Phổ đo bức xạ gamma năng lượng 1408 keV.

PHƯƠNG PHÁP MONTE CARLO VÀ CHƯƠNG TRÌNH MCNP5 2. Phương pháp Monte Carlo Sự ra đời của máy tính điện tử đầu những năm 1940 là một bước ngoặc quan trọng đối với sự phát triển của xã hội loài người. Kích thước của những hệ máy tính điện tử khi ấy rất cồng kềnh, điển hình như máy tính ENIAC được thiết kế bởi hai nhà khoa học Mỹ là John Mauchy và J. Bởi vì việc tính toán bằng máy tính vẫn còn phức tạp, người vận hành phải sử dụng bằng những máy điện cơ, nên các kỹ sư phải lựa chọn phương pháp giải toán phù hợp với khối lượng tính toán từ máy tính.

Một trong những phương pháp phù hợp với máy tính điện tử là giải toán bằng phương pháp số. Tuy nhiên, ý tưởng giải toán bằng phương pháp số đã xuất hiện từ rất sớm nhưng chưa được quan tâm. Năm 1777, nhà toán học người Pháp Georges- Louis Leclerc Comte de Buffon đã đưa ra một bài toán nổi tiếng về tính số  hết sức kỳ lạ.