Khóa luận: Ứng dụng kỹ thuật truyền qua đo độ dày kim loại (SPKHTN)

Khóa luận sư phạm: Ứng dụng kỹ thuật truyền qua đo độ dày kim loại. Nghiên cứu sử dụng nguồn cường độ thấp, phù hợp cho giáo dục khoa học tự nhiên.

• Tác giả: Lê Mỹ Hoà
• Người hướng dẫn: Pgs. Hoàng Đức Tâm
• Trường học: Trường Đại Học Sư Phạm Thành Phố Hồ Chí Minh
• Chuyên ngành: Sư phạm Khoa Học Tự Nhiên
• Đề tài: Ứng Dụng Kĩ Thuật Truyền Qua Để Đo Độ Dày Của Một Số Vật Liệu Kim Loại Dùng Nguồn Cường Độ Thấp
• Loại tài liệu: Khóa Luận Tốt Nghiệp Đại Học
• Năm xuất bản: 2023
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• LỜI CAM ĐOAN
• DANH MỤC CÁC BẢNG BIỂU
• DANH MỤC CÁC HÌNH
• DANH MỤC NHỮNG TỪ VIẾT TẮT
• 1. CHƯƠNG 1: TƯƠNG TÁC CỦA BỨC XẠ GAMMA VỚI VẬT CHẤT
• 1.1. Tương tác của bức xạ gamma với vật chất
• 1.2. Hiệu ứng quang điện
• 1.3. Tán xạ Compton
• 1.4. Hiệu ứng tạo cặp
• 1.5. Hệ số suy giảm khối của vật liệu trong tương tác với bức xạ
• 1.6. Sự suy giảm cường độ chùm tia gamma qua bề dày vật liệu
• 2. CHƯƠNG 2: PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO VÀ CHƯƠNG TRÌNH MCNP6
• 2.1. Phương pháp Monte Carlo
• 2.2. Chương trình MCNP6
• 2.3. Cấu trúc của một tập tin đầu vào trong chương trình MCNP6
• 3. CHƯƠNG 3: BỐ TRÍ THỰC NGHIỆM VÀ MÔ PHỎNG DỮ LIỆU TRONG PHÉP ĐO BỀ DÀY KIM LOẠI
• 3.1. Vật liệu và bố trí thực nghiệm
• 3.2. Mô hình xây dựng trong chương trình MCNP6
• 4. CHƯƠNG 4: XÁC ĐỊNH BỀ DÀY MỘT SỐ KIM LOẠI BẰNG KỸ THUẬT ĐO GAMMA TRUYỀN QUA VỚI NGUỒN CƯỜNG ĐỘ THẤP
• 4.1. Đường chuẩn tuyến tính của cường độ chùm tia gamma truyền qua với bề dày của một số kim loại
• 4.2. Hệ số suy giảm khối
• 4.3. Xác định bề dày theo đường chuẩn dựa trên phép đo thực nghiệm
• KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Khám phá phương pháp đo độ dày kim loại truyền qua là gì

Phương pháp đo độ dày kim loại bằng kỹ thuật truyền qua là một trong những ứng dụng quan trọng nhất của lĩnh vực kiểm tra không phá hủy NDT (Non-Destructive Testing). Kỹ thuật này sử dụng một chùm bức xạ, thường là tia gamma hoặc tia X, chiếu xuyên qua vật liệu cần đo. Một đầu dò đặt ở phía đối diện sẽ ghi nhận cường độ chùm tia sau khi đã bị suy giảm bởi vật liệu. Dựa trên định luật suy giảm Beer-Lambert, độ dày của vật liệu có thể được xác định chính xác thông qua mức độ suy giảm cường độ này. Đây là một phương pháp cốt lõi trong ultrasonic testing UT và các kỹ thuật bức xạ, cho phép đánh giá toàn vẹn kết cấu, kiểm tra ăn mòn, và đảm bảo chất lượng sản phẩm trong các ngành công nghiệp nặng, hàng không và năng lượng. Phương pháp này đặc biệt hiệu quả vì không yêu cầu phá hủy mẫu, giúp tiết kiệm chi phí và thời gian, đồng thời cung cấp dữ liệu chính xác về tình trạng bên trong của vật liệu kim loại.

1.1. Tầm quan trọng của kiểm tra không phá hủy NDT

Kiểm tra không phá hủy (NDT) là tập hợp các kỹ thuật phân tích được sử dụng trong khoa học và công nghiệp để đánh giá các đặc tính của vật liệu, bộ phận hoặc hệ thống mà không gây ra thiệt hại. Các phương pháp NDT như siêu âm đo độ dày hay kỹ thuật truyền qua gamma đóng vai trò sống còn trong việc đảm bảo an toàn và độ tin cậy của các công trình kỹ thuật. Chúng giúp phát hiện sớm các khuyết tật vật liệu ẩn bên trong như vết nứt, rỗ khí, hoặc sự mỏng đi do ăn mòn trước khi chúng phát triển thành các sự cố nghiêm trọng. Trong các ngành như dầu khí, hàng không vũ trụ và xây dựng, việc kiểm tra kết cấu thép và đo độ dày thành ống định kỳ bằng NDT là quy định bắt buộc để ngăn ngừa thảm họa, kéo dài tuổi thọ thiết bị và tối ưu hóa chi phí bảo trì.

1.2. Nguyên lý cốt lõi của kỹ thuật truyền qua gamma

Nguyên lý cơ bản của phương pháp truyền qua dựa trên sự suy giảm cường độ của chùm tia gamma khi đi qua vật chất. Theo định luật Beer-Lambert, cường độ chùm tia sau khi đi qua vật liệu (I) có mối quan hệ hàm mũ với cường độ ban đầu (I₀), hệ số suy giảm tuyến tính của vật liệu (µ) và độ dày vật liệu (d), biểu diễn qua công thức: I = I₀ * exp(-µd). Bằng cách đo I và I₀, và biết trước hệ số µ của vật liệu, độ dày d có thể được tính toán. Trong thực tế, hệ thống bao gồm một nguồn phát sóng siêu âm hoặc gamma, vật liệu cần đo và một đầu dò phát và thu. Nghiên cứu của Lê Mỹ Hoà (2023) đã áp dụng nguyên lý này bằng cách sử dụng nguồn gamma ¹³⁷Cs và đầu dò NaI(Tl) để xác định tỉ số R = I/I₀, từ đó xây dựng mối quan hệ tuyến tính giữa ln(R) và độ dày, cho phép xác định độ dày một cách chính xác.

1.3. So sánh kỹ thuật truyền qua và kỹ thuật tán xạ ngược

Kỹ thuật truyền qua và tán xạ ngược đều là các phương pháp NDT bức xạ nhưng hoạt động theo nguyên lý khác nhau. Kỹ thuật truyền qua yêu cầu tiếp cận cả hai phía của vật thể, với nguồn và đầu dò đặt đối diện nhau. Nó cho độ chính xác cao và xử lý số liệu đơn giản hơn. Ngược lại, kỹ thuật tán xạ ngược (backscatter) chỉ yêu cầu tiếp cận một phía, vì đầu dò và nguồn được đặt cùng một bên. Đầu dò sẽ ghi nhận các bức xạ bị tán xạ ngược lại từ vật liệu. Phương pháp này hữu ích khi không thể tiếp cận mặt sau của vật thể, ví dụ như kiểm tra đường ống ngầm hoặc thành bồn chứa lớn. Tuy nhiên, theo các nghiên cứu được trích dẫn trong khóa luận của Lê Mỹ Hoà, kỹ thuật tán xạ thường phức tạp hơn trong xử lý dữ liệu do ảnh hưởng của tán xạ đa cấp và có sai số lớn hơn, trong khi phương pháp truyền qua cho thấy độ tin cậy cao hơn, đặc biệt khi đo các vật liệu có mật độ lớn.

II. Thách thức khi đo độ dày kim loại bằng nguồn cường độ thấp

Việc sử dụng nguồn phóng xạ cường độ thấp trong kỹ thuật đo độ dày truyền qua mang lại lợi ích lớn về an toàn bức xạ và chi phí, nhưng cũng đặt ra những thách thức kỹ thuật đáng kể. Nguồn cường độ thấp tạo ra dòng photon (số đếm) ít hơn, dẫn đến tín hiệu yếu và tỉ số tín hiệu trên nhiễu (Signal-to-Noise Ratio) thấp. Điều này đòi hỏi thời gian đo dài hơn để thu thập đủ dữ liệu thống kê, làm cho quy trình trở nên kém hiệu quả đối với các ứng dụng cần tốc độ. Hơn nữa, tín hiệu yếu dễ bị ảnh hưởng bởi các yếu tố môi trường và sự bất ổn định của thiết bị siêu âm đo độ dày. Các vấn đề như sự dịch đỉnh phổ do thay đổi nhiệt độ của đầu dò, như được ghi nhận trong nghiên cứu thực nghiệm của Lê Mỹ Hoà (2023), trở nên rõ rệt hơn và có thể gây ra sai số đáng kể nếu không được xử lý đúng cách. Việc khắc phục những thách thức này đòi hỏi các phương pháp xử lý tín hiệu tiên tiến và quy trình hiệu chuẩn thiết bị đo nghiêm ngặt.

2.1. Ưu điểm và hạn chế của nguồn phóng xạ cường độ thấp

Ưu điểm chính của việc sử dụng nguồn cường độ thấp là an toàn. Nó giảm thiểu nguy cơ phơi nhiễm bức xạ cho người vận hành và môi trường xung quanh, đồng thời giảm yêu cầu về che chắn và các biện pháp an toàn phức tạp, giúp giảm chi phí thiết lập hệ thống. Tuy nhiên, hạn chế lớn nhất là số liệu thống kê thấp. Thời gian cần thiết để đạt được độ chính xác mong muốn có thể kéo dài hàng giờ. Điều này đặc biệt bất lợi khi đo độ dày thép hoặc các vật liệu dày, nơi sự suy giảm tín hiệu đã rất lớn. Khóa luận của Lê Mỹ Hoà đã chỉ ra rằng để có kết quả tin cậy với nguồn ¹³⁷Cs hoạt độ chỉ 0.71µCi, thời gian đo cho mỗi mẫu lên tới 14,400 giây (4 giờ).

2.2. Ảnh hưởng của nhiễu và sự bất ổn định của thiết bị đo

Với tín hiệu yếu từ nguồn cường độ thấp, sự ổn định của hệ thống đo trở thành yếu tố quyết định. Đầu dò siêu âm và các thiết bị điện tử đi kèm rất nhạy cảm với sự thay đổi của môi trường. Nghiên cứu thực nghiệm đã ghi nhận hiện tượng dịch đỉnh phổ năng lượng giữa các lần đo khác nhau. Nguyên nhân được xác định là do sự không ổn định nhiệt của đầu dò NaI(Tl). Sự thay đổi nhiệt độ dù nhỏ cũng có thể làm thay đổi đặc tính của tinh thể nhấp nháy và bộ nhân quang điện, dẫn đến sự dịch chuyển vị trí đỉnh năng lượng trên phổ đồ. Vấn đề này làm phức tạp hóa quá trình xử lý phổ, yêu cầu phải xác định lại vùng quan tâm (ROI) cho mỗi lần đo và có thể gây ra sai số trong việc tính toán diện tích đỉnh, từ đó ảnh hưởng trực tiếp đến độ chính xác của độ dày đo được.

III. Hướng dẫn đo độ dày kim loại qua mô phỏng Monte Carlo

Để khắc phục những hạn chế của thực nghiệm, đặc biệt là thời gian đo kéo dài, phương pháp mô phỏng Monte Carlo nổi lên như một công cụ cực kỳ mạnh mẽ. Bằng cách sử dụng các chương trình chuyên dụng như MCNP6, các nhà khoa học có thể tái tạo lại toàn bộ quá trình tương tác của bức xạ với vật chất trên máy tính. Mô phỏng cho phép khảo sát một dải rộng các độ dày và vật liệu khác nhau mà không cần mẫu vật lý. Quan trọng nhất, nó cho phép xây dựng các đường chuẩn lý thuyết có độ chính xác cao. Dữ liệu từ mô phỏng sau đó được sử dụng làm cơ sở để đối chiếu và hiệu chỉnh kết quả từ một số ít các phép đo thực nghiệm. Sự kết hợp này, như được trình bày trong nghiên cứu của Lê Mỹ Hoà, giúp giảm đáng kể thời gian và chi phí nghiên cứu trong khi vẫn đảm bảo độ tin cậy của phương pháp truyền qua.

3.1. Thiết lập mô hình đo trong chương trình MCNP6

Mô phỏng một hệ đo gamma truyền qua trong MCNP6 đòi hỏi việc khai báo chi tiết và chính xác về mặt hình học, vật liệu và nguồn bức xạ. Người dùng phải định nghĩa cấu trúc không gian ba chiều của toàn bộ hệ thống, bao gồm nguồn phóng xạ (ví dụ: ¹³⁷Cs), ống chuẩn trực bằng chì, tấm kim loại mẫu (nhôm, sắt) với các độ dày khác nhau, và cấu trúc chi tiết của đầu dò siêu âm (tinh thể NaI(Tl), lớp phản xạ, vỏ nhôm). Các vật liệu phải được định nghĩa bằng thành phần nguyên tố và mật độ chính xác. Nguồn bức xạ được khai báo với năng lượng phát đặc trưng (661.7 keV cho ¹³⁷Cs). Việc chạy mô phỏng với số lượng hạt lớn (hàng tỷ hạt) đảm bảo kết quả có độ chính xác thống kê cao, tạo ra các phổ năng lượng ảo có thể so sánh trực tiếp với phổ thực nghiệm.

3.2. Xây dựng đường chuẩn tuyến tính từ dữ liệu mô phỏng

Sau khi thu được phổ năng lượng từ mô phỏng MCNP6 cho từng độ dày kim loại, bước tiếp theo là xử lý dữ liệu để xây dựng đường chuẩn. Diện tích đỉnh năng lượng toàn phần (tương ứng với số đếm N) được tính toán. Tỉ số R được xác định bằng cách lấy N của trường hợp có mẫu chia cho N của trường hợp không có mẫu (phông). Theo lý thuyết, logarit tự nhiên của tỉ số này, ln(R), có mối quan hệ tuyến tính với độ dày d: ln(R) = - (µ/ρ) * ρ * d, trong đó (µ/ρ) là hệ số suy giảm khối. Bằng cách khớp các điểm dữ liệu (ln(R), d) từ mô phỏng bằng một hàm tuyến tính, một đường chuẩn dạng y = a + bx được tạo ra. Đường chuẩn này chính là chìa khóa để suy ra độ dày của một mẫu chưa biết khi đo tỉ số R của nó trong thực nghiệm.

3.3. Xác định hệ số suy giảm khối và so sánh với NIST

Một trong những bước quan trọng để xác thực độ tin cậy của mô hình mô phỏng là tính toán hệ số suy giảm khối (mass attenuation coefficient) và so sánh nó với các giá trị tiêu chuẩn từ cơ sở dữ liệu uy tín như NIST (Viện Tiêu chuẩn và Công nghệ Quốc gia Hoa Kỳ). Hệ số này là một đại lượng vật lý đặc trưng cho khả năng một vật liệu làm suy yếu bức xạ gamma ở một mức năng lượng nhất định. Từ đường chuẩn mô phỏng, hệ số này có thể được trích xuất từ độ dốc của đường thẳng. Trong nghiên cứu của Lê Mỹ Hoà (2023), kết quả hệ số suy giảm khối tính toán từ mô phỏng cho nhôm và sắt có độ lệch tương đối so với dữ liệu của NIST đều dưới 4%. Kết quả này khẳng định rằng mô hình mô phỏng MCNP6 đã được xây dựng một cách chính xác và đáng tin cậy.

IV. Quy trình đo độ dày kim loại truyền qua trong thực nghiệm

Quy trình thực nghiệm là bước kiểm chứng cuối cùng cho các kết quả từ mô phỏng và là nền tảng để ứng dụng phương pháp vào thực tế. Một quy trình chuẩn bao gồm các bước: bố trí hệ thống đo, hiệu chuẩn thiết bị, tiến hành đo lường, và xử lý dữ liệu thô. Sự cẩn trọng trong từng bước là yếu tố quyết định đến độ chính xác của kết quả cuối cùng. Đặc biệt với nguồn cường độ thấp, việc đảm bảo sự đồng trục giữa nguồn, mẫu và đầu dò siêu âm, duy trì điều kiện môi trường ổn định (nhiệt độ), và sử dụng chất tiếp âm siêu âm (trong trường hợp siêu âm) hoặc đảm bảo tiếp xúc tốt là cực kỳ quan trọng. Quy trình này đòi hỏi sự kiên nhẫn và tỉ mỉ, từ việc đo đạc các kích thước hình học đến việc xử lý từng phổ năng lượng thu được để loại bỏ nhiễu và các yếu tố gây sai lệch.

4.1. Bố trí hệ thống đo với đầu dò NaI Tl và nguồn ¹³⁷Cs

Hệ thống thực nghiệm được bố trí sao cho tâm của nguồn phóng xạ ¹³⁷Cs, tâm của mẫu kim loại và tâm của đầu dò NaI(Tl) nằm trên cùng một đường thẳng (đồng trục). Nguồn được đặt trong một khối chì che chắn có một ống chuẩn trực để tạo ra một chùm tia gamma hẹp. Mẫu vật liệu được đặt giữa nguồn và đầu dò. Đầu dò phát và thu tín hiệu là một tinh thể NaI(Tl) kết nối với một bộ phân tích đa kênh (MCA) để ghi nhận phổ năng lượng. Toàn bộ hệ thống được đặt trong môi trường được kiểm soát nhiệt độ (23-24°C) để hạn chế sự dịch đỉnh phổ. Khoảng cách giữa các thành phần được đo đạc và giữ cố định trong suốt quá trình thí nghiệm để đảm bảo tính nhất quán của phép đo. Mỗi phép đo, bao gồm cả đo phông (không có mẫu) và đo với mẫu, được thực hiện trong một khoảng thời gian dài (14,400 giây) để tích lũy đủ số đếm.

4.2. Kỹ thuật xử lý phổ gamma và tính toán tỉ số R

Phổ gamma thô thu được từ MCA là một biểu đồ thể hiện số đếm theo kênh năng lượng. Để tính toán chính xác số đếm trong đỉnh năng lượng toàn phần (ví dụ: 661.7 keV), cần phải xử lý phổ. Quá trình này bao gồm việc trừ phông nền Compton, là phần nền liên tục nằm dưới đỉnh phổ. Các phần mềm chuyên dụng như Colegram được sử dụng để khớp một hàm Gaussian cho đỉnh phổ và một hàm nền thích hợp. Diện tích của hàm Gaussian sau khi trừ phông chính là số đếm thực của chùm tia truyền qua (N). Sau khi có được số đếm cho trường hợp có mẫu (N_exp) và không có mẫu (N_phong), tỉ số thực nghiệm R_exp = N_exp / N_phong được tính toán. Đây là giá trị đầu vào quan trọng để sử dụng đường chuẩn mô phỏng nhằm xác định độ dày.

V. Kết quả đo độ dày kim loại và ứng dụng thực tiễn

Sự kết hợp giữa mô phỏng và thực nghiệm đã chứng tỏ hiệu quả vượt trội trong việc xác định độ dày kim loại, ngay cả khi sử dụng nguồn cường độ thấp. Bằng cách thay giá trị ln(R_exp) thu được từ thực nghiệm vào phương trình đường chuẩn xây dựng từ mô phỏng, độ dày của các tấm kim loại đã được tính toán với độ chính xác cao. Kết quả nghiên cứu cho thấy, đối với các mẫu nhôm có độ dày trên 10mm, độ lệch tương đối giữa giá trị đo được và giá trị tham chiếu (đo bằng thước kẹp) là dưới 2%. Đây là một kết quả rất tốt, khẳng định tính khả thi của phương pháp. Những kết quả này mở ra tiềm năng ứng dụng rộng rãi trong công nghiệp, từ việc đo chiều dày lớp phủ bảo vệ, kiểm tra ăn mòn định kỳ trên các đường ống và bồn chứa, đến việc đảm bảo chất lượng của các kết cấu thép trong xây dựng.

5.1. Đối chiếu kết quả đo với đường chuẩn mô phỏng

Đây là bước quan trọng nhất của phương pháp kết hợp. Các giá trị tỉ số R thu được từ thực nghiệm (R_exp) cho các mẫu nhôm có độ dày khác nhau được xử lý và tính toán. Sau đó, giá trị ln(R_exp) được đưa vào phương trình đường chuẩn tuyến tính (y = a + bx) đã được thiết lập từ dữ liệu mô phỏng MCNP6, trong đó y = ln(R) và x là độ dày. Độ dày thực nghiệm (x) được giải ra từ phương trình. Kết quả từ khóa luận của Lê Mỹ Hoà (2023) cho thấy sự phù hợp tốt giữa dữ liệu thực nghiệm và đường chuẩn mô phỏng. Ví dụ, với mẫu nhôm dày 50.35mm, độ dày tính toán được là 50.28mm, với độ lệch tương đối chỉ 0.13%, cho thấy sự chính xác ấn tượng của phương pháp.

5.2. Đánh giá độ chính xác và độ lệch tương đối phép đo

Độ chính xác của phương pháp được đánh giá thông qua độ lệch tương đối (RD%) giữa độ dày đo được và độ dày tham chiếu. Kết quả cho thấy phương pháp này có độ chính xác cao đối với các mẫu tương đối dày. Tuy nhiên, đối với các mẫu mỏng (ví dụ: 4.90mm), độ lệch tương đối cao hơn (12.2%). Nguyên nhân là do với mẫu mỏng, sự suy giảm tín hiệu rất nhỏ, khiến cho tỉ số R gần bằng 1. Điều này làm cho sai số thống kê của phép đếm có ảnh hưởng lớn hơn đến kết quả cuối cùng. Điều này cho thấy có một giới hạn về độ dày tối thiểu mà phương pháp có thể đo chính xác với một nguồn cường độ thấp và thời gian đo nhất định. Để cải thiện độ chính xác cho mẫu mỏng, cần tăng thời gian đo để thu thập nhiều số liệu thống kê hơn.

VI. Kết luận Hướng phát triển kỹ thuật đo độ dày kim loại

Nghiên cứu ứng dụng kỹ thuật truyền qua gamma với nguồn cường độ thấp, kết hợp với mô phỏng Monte Carlo, đã chứng minh là một phương pháp hiệu quả và đáng tin cậy để xác định độ dày vật liệu kim loại. Phương pháp này không chỉ đảm bảo an toàn bức xạ mà còn giảm thiểu chi phí so với việc sử dụng nguồn cường độ cao. Sự phù hợp tốt giữa kết quả mô phỏng, dữ liệu NIST và các phép đo thực nghiệm đã khẳng định tính đúng đắn của phương pháp luận. Hướng phát triển trong tương lai có thể tập trung vào việc tối ưu hóa quy trình để tăng tốc độ đo và cải thiện độ chính xác cho các vật liệu mỏng. Việc sử dụng các loại đầu dò có hiệu suất cao hơn, ổn định nhiệt tốt hơn và áp dụng các thuật toán xử lý tín hiệu tiên tiến sẽ là những bước đi quan trọng tiếp theo để đưa kỹ thuật này vào ứng dụng công nghiệp rộng rãi.

6.1. Tóm tắt ưu điểm của phương pháp gamma truyền qua

Phương pháp gamma truyền qua kết hợp mô phỏng sở hữu nhiều ưu điểm nổi bật. Thứ nhất, nó là một kỹ thuật kiểm tra không phá hủy NDT, bảo toàn nguyên vẹn mẫu vật. Thứ hai, việc sử dụng nguồn cường độ thấp giúp tăng cường an toàn và giảm chi phí. Thứ ba, sự kết hợp với mô phỏng MCNP6 cho phép xây dựng đường chuẩn chính xác mà không cần nhiều mẫu vật lý, tiết kiệm thời gian và nguồn lực. Thứ tư, phương pháp cho kết quả chính xác cao, đặc biệt với các vật liệu có độ dày trung bình và lớn, với độ lệch dưới 2%. Cuối cùng, nó có khả năng ứng dụng đa dạng trong nhiều lĩnh vực công nghiệp để kiểm tra ăn mòn và đảm bảo chất lượng sản phẩm.

6.2. Kiến nghị cải tiến để tăng độ chính xác cho phép đo

Để nâng cao hơn nữa hiệu quả của phương pháp, một số cải tiến có thể được thực hiện. Đầu tiên, cần sử dụng các loại đầu dò siêu âm hoặc gamma có độ ổn định nhiệt cao hơn, ví dụ như các đầu dò được làm mát, để loại bỏ hiện tượng dịch đỉnh phổ. Thứ hai, đối với việc đo các mẫu vật liệu mỏng, cần tăng thời gian đo một cách đáng kể để cải thiện số liệu thống kê và giảm sai số tương đối. Thứ ba, có thể nghiên cứu áp dụng các thuật toán xử lý phổ tiên tiến hơn để tách đỉnh và trừ phông hiệu quả hơn, đặc biệt trong điều kiện tín hiệu yếu. Cuối cùng, việc phát triển một hệ thống đo tự động hóa có thể giúp thực hiện các phép đo dài hạn một cách hiệu quả và giảm thiểu sai sót do con người.

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

Chương 1. Tương tác bức xạ gamma với vật chất. Chương này trình bày về những tương tác của bức xạ gamma với vật chất, bao gồm hiệu ứng quang điện. tán xạ Compton, hiệu ứng tạo cặp, sự suy giảm cường độ gamma qua vật chất và hệ số suy giảm của vật liệu trong tương tác với bức xạ.

Phương pháp mô phỏng Monte Carlo và chương trình MCNP6. Chương này trình bày về đặc điểm của chương trình MCNP6 và phương pháp Monte Carlo. Bồ trí thực nghiệm và mô phỏng dit liệu trong phép đo bề dày kim loại. Trong chương nảy trình bày các thông tin cần thiết cho việc mô phỏng, cách bố trí và vật liệu thực nghiệm.

Xác định bẻ day một số kim loại bằng kỹ thuật đo gamma truyền qua với nguồn cường độ thấp. Chương này chúng tôi trình bày các kết quả đạt được từ việc mô phỏng, so sánh với cơ sở đữ liệu NIST và so với kết quả thực nghiệm. TƯƠNG TÁC CUA BỨC XA GAMMA VỚI VAT CHAT Bản chất của bức xa gamma là sóng điện từ mang năng lượng cao va có bước song rat nhỏ (cỡ 10%cm) so với kích thước nguyên tử. Bức xa gamma thực chất là các hạt photon hay còn được gọi là lượng tử gamma, có tính chất của sóng và hạt, khi đi vảo vật liệu photon tương tac với các clectron.

Co ba dang tương tác cơ bản của bức xạ gamma với electron được kẻ đến là hiệu ứng quang điện, tán xạ Compton và hiệu ứng tạo cặp. Trong chương nảy, chúng tôi sẽ trình bày các cơ chế tương tác của bức xạ gamma với vật chất. Tương tác của bức xạ gamma với vật chat 1. Hiệu ứng quang điện Hiệu ứng quang điện là quá trình tương tác của photon với electron quỹ đạo của nguyên tử, khi đó photon biển mat và truyền toàn bộ năng lượng cho electron liên kết làm cho electron bị đánh bật ra khỏi qui đạo và trở thành electron quang điện.

Theo thuyết lượng tử ánh sáng, động năng của electron quang điện [11]: K,=hv-A (1.1) Trong đó: -K : động nang của electron quang điện -hy: năng lượng photon tới - A: năng lượng liên kết của electron Hiệu ứng quang điện chủ yếu xay ra với các electron có năng lượng thấp (<1MeV). Quá trình hap thu quang điện có hiệu suất lớn nhất khi tương tác xảy ra với các electron liên kết chặt chẽ với hạt nhân nguyên tử (lớp K). Hiệu ứng quang điện không xảy ra đối với các electron tự do. Khi quá trình này xảy ra, tôn tại một lỗ trống tại vị trí xảy ra hiệu ứng.

Các electron lớp ngoài nhanh chóng chiếm chỗ và dẫn đến quá trình tạo tia X đặc trưng hay các electron Auger. Tiết điện của hiệu ứng quang điện chủ yêu phụ thuộc vào năng lượng của chùm tia gamma tới và điện tích Z của hạt nhân nguyên tử. Hiệu ứng quang điện chiếm ưu thé ở vùng năng lượng thấp và vật liệu hấp thụ có Z lớn. Tiết điện hap thụ quang điện được tính theo công thức [12]: O45 = Z°.2) © - (quang điện tửi Hình 1.

Mô hình hiệu ứng quang điện [13] 1. Tan xạ Compton Tan xạ Compton là sự va chạm giữa photon vả clectron tự do ở lớp ngoài, trong đó photon sẽ truyền một phan năng lượng cho electron, electron hap thụ năng lượng từ photon sẽ bit ra khỏi nguyên tử, photon bị mat một phan năng lượng sau khi tương tác sẽ bị lệch phương so với quỹ đạo ban đầu. Mô hình tán xạ Compton [13] Ap dụng định luật bảo toàn năng lượng và động lượng thì động năng của electron sau tan xạ [11]: hy P 8, ) 2, (1-cos K_=E,+E, =ht—=hw`= lv~ ———— 8 -—--- 1.) ; Trong đó: * hy hay E, : Năng lượng photo ban dau *hy hay E ,: Năng lượng của photon sau khi tan xạ *Â= be mạc? =51lkeV : Năng lượng nghỉ của electron myc” «9: Góc tan xạ của chùm tia gamma, có thé thay đôi từ 0° đến 180°. Khi góc tán xạ lớn hơn 90° thì hiện tượng gamma tán xạ ngược Xảy ra.

Tiết diện tan xạ Compton chiếm ưu thé trong vùng năng lượng khoảng vai MeV. Tiết điện tán xạ Compton được xác định theo công thức [12]: ơ,~Z. Hiệu ứng tạo cặp Hiệu ứng tạo cặp là qua trình tương tac giữa một photon với hạt nhân nguyên tử. Điều kiện dé hiện tượng này xảy ra khi tia gamma tương tác phải có năng lượng tối thiểu bằng năng lượng ở trạng thái nghỉ của hai electron (tức trên 1022keV).

Kết quả là sự biến mat của photon củng với sự xuất hiện cặp eletron và positron. Positron sau đó nhanh chóng bị huỷ do tương tác với các electron khác từ môi trường va sinh ra hai photon với năng lượng Š1 IkeV. Tiết diện xảy ra hiệu ứng tạo cặp tăng khi nang lượng của chùm tia gamma tới tăng và điện tích hạt nhân Z lớn như chỉ, uranium. Peeve 1002 kev E, « — \ “A À1 ì „Hs bu ® tT e ọ es eT øu \ b we, ø a Sa t11 kev Hình 1.

Mô hình hiệu ứng tao cap [13] 1. Hệ số suy giảm khối của vật liệu trong tương tác với bức xa Dé ghi nhận được sự suy giám của bức xạ trong môi trường, ta quan tâm đến xác suất xảy ra ba quá trình (hiệu ứng quang điện, tán xạ Compton và hiệu ứng tạo cặp). Dai lượng tiết diện tương tác toàn phan o là tổng tiết điện của tất cả các quá trình tương tác giữa photon và vật chat. Tiết diện toàn phần vi mô [12]: Ø=Ơ,nto, +Ơ, (1.6) Với: ® Ø„„ là tiết điện tương tác của hiệu ứng quang điện.

là tiết diện tương tác của tán xa Compton. °c, là tiết điện tương tác của hiệu ứng tạo cặp. Khi nhân tiết diện vi mô với số nguyên tử N có trong lem? ta được hệ số suy giảm tuyến tính [12]: ua=Nơ (1.7) Với ¿ là hệ số suy giảm tuyến tính. Hệ số suy giảm tuyến tính có thể được biểu diễn qua hệ số suy giảm khối.

Hệ số suy giảm khối không phụ thuộc vào mật độ vật chất và có thé áp dụng cho bat kỳ dang vật chat nao (rắn, lỏng, khí). Nó được tinh bởi công thức [12]: tử, = L (cm*/g) (1.8) p Với p là mật độ vật chat (g/cm). Nếu xác định hệ số suy giảm tuyến tinh dựa vào số đếm. ta có công thức [12]: N=N,.9) Với: *N là số đếm dau do ghi nhận được ở trường hợp có mẫu đo.

« N, là số đếm đầu đò ghi nhận được ở trường hợp không có mẫu đo. Hệ số suy giảm khối (t„) còn được xác định theo số dém ghi nhận bởi đầu đò thông qua công thức sau: i= Lf 2 Ả px \N,} (1. Sự suy giảm cường độ chùm tia gamma qua bề dày vật liệu Theo định luật Lambert - Beer, cường độ của chùm tia gamma sẽ bị suy giảm khi đi qua một vật liệu có bé day d theo quy luật hàm mũ: f= I„exp(— ,4,) (photon.11) Trong đó, In, I lần lượt là cường độ của chim tia gamma phát trực tiếp ra từ nguồn phóng xạ và của chùm tia mà đầu đò ghi nhận khi đi qua vật liệu, /¿(em'”) là hệ số Suy giảm tuyến tính của vật liệu thứ i có be day đ (cm). Trong khoá luận, chúng tôi thực hiện phép đo cường độ chùm tia qua một tam kim loại có bề day đ,.

tr là hé số suy giảm khối của kim loại có bề day d. pix là hệ số suy giảm tuyến tính của không khí. Bên cạnh chùm tia truyền qua bị suy giảm thì dau đò cũng ghi nhận một thành phan khác là chim tia tan xa Compton với góc nhỏ khi tương tác với vật liệu. Tuy nhiên, trong trường hợp của khóa luận, cả đầu đò và nguồn đều được chuan trực với một đường kính nhỏ nên có thé bỏ qua chùm tia tán xạ này.

Như vậy, diện tích đỉnh của cham tia gamma truyền qua có thé xác định từ một phô phân bồ dang xung được ghi nhận bởi đầu đò theo công thức sau: Nx, = @6tl , = œet.13) Với ¢ là hiệu suất ghi nhận đính năng lượng toàn phan của đầu đò, œ (cm?) là tiết điện mặt cat của chùm tia, t (s) là thời gian ghi nhận của đầu dò. Trường hợp khi không có mẫu thì diện tích đính của chùm tia là: N,, = đet.14) Khi đó, tỉ số R được tính như sau: R= Này =exp(T-,đ) (1.15) kế Thay fy, = My, P. Ta có: InR =~yy, pd (1.16) cho thay giữa bề dày va ti số InR có mối quan hệ tuyến tinh với nhau. Vì vậy, có thê viết công thức (1.16) lại dưới dạng y = a + bx, trong đó các hệ số a, b phụ thuộc vao tính chất của vật liệu.

Trong khóa luận, các hệ sé này được xác định từ đường chuẩn xây dựng bằng đữ liệu mô phỏng MCNP6 và kiểm chứng lại bằng dữ liệu thực nghiệm. PHUONG PHAP MO PHONG MONTE CARLO VA CHUONG TRÌNH MCNP6 2. Phuong phap Monte Carlo Phuong pháp Monte Carlo là một lớp các thuật toán dé giai quyết nhiều bai toán trên máy tính theo kiều không dự đoán trước được kết qua từ dit liệu đầu vào. Phương pháp nay được xây dựng dựa trên nên tang các số ngẫu nhiên, luật số lớn vả định lý giới hạn trung tâm.

Tên gọi Monte Carlo đề cập đến các song bạc ở Monaco. Vào khoảng thé ki 18 phương pháp Monte Carlo đã thực hiện, tuy nhiên phương pháp này không được công bé trong bat kỳ công trình nghiên cứu nao. Phương pháp Monte Carlo chỉ được thực sự sử dụng như một công cụ nghiên cứu trong việc chế tạo bom nguyên tử vào thế chiến thứ 2. Và đến khoảng những năm 1970, các lý thuyết về phương pháp dan được hoàn thiện với độ phức tạp và độ chính xác cao hơn.

Ngày nay, theo sự phát triển của công nghệ hiện đại thì phương pháp Monte Carlo là phương pháp mô phông rat quan trọng trong nhiều lĩnh vực khoa học. Chương trình MCNP6 MCNP (Monte Carlo N-Particle) là chương trình ứng dụng phương pháp Monte Carlo để mô phỏng các quá trình vật lý hạt nhân đối với neutron, photon và electron. Chương trình cho phép người dùng sử dụng các dạng hình học phức tạp và mô phỏng dựa trên các thư viện dữ liệu tương tác hạt nhân. Chương trình điều khiển các quá trình tương tác bang cách gieo số ngau nhiên theo quy luật thong kê cho trước và mô phỏng thực hiện trên máy tính vì số lần thử cần thiết thường rat lớn.

C huong trinh duge xay dựng bởi phòng thi nghiệm quốc gia Los Alamos. Khoảng từ hai đến ba năm họ lại cho ra một phiên bản mới của chương trình và trong nghiên cứu này chúng tôi sử dụng phiên ban MCNP6.